The high intensity of free electron lasers opens up the possibility to perform single-shot molecule scattering experiments. However, even for small molecules, radiation damage induced by absorption of high intense x-ray radiation is not yet fully understood. One of the striking effects which occurs under intense x-ray illumination is the creation of double core ionized molecules in considerable quantity. To provide insight into this process, we have studied the dynamics of water molecules in single and double core ionized states by means of electronic transition rate calculations and ab initio molecular dynamics (MD) simulations. From the MD trajectories, photoionization and Auger transition rates were computed based on electronic continuum wavefunctions obtained by explicit integration of the coupled radial Schrödinger equations. These rates served to solve the master equations for the populations of the relevant electronic states. To account for the nuclear dynamics during the core hole lifetime, the calculated electron emission spectra for different molecular geometries were incoherently accumulated according to the obtained time-dependent populations, thus neglecting possible interference effects between different decay pathways. We find that, in contrast to the single core ionized water molecule, the nuclear dynamics for the double core ionized water molecule during the core hole lifetime leaves a clear fingerprint in the resulting electron emission spectra. The lifetime of the double core ionized water was found to be significantly shorter than half of the single core hole lifetime.

Organisms have evolved a wide variety of mechanisms to utilize and respond to light. In many cases, the biological response is mediated by structural changes that follow photon absorption in a protein complex. The initial step in such cases is normally the photoisomerization of a highly conjugated prosthetic group. To understand better the factors controlling the isomerization, we perform atomistic molecular dynamics simulations. In this perspective article we briefly review the key theoretical concepts of photochemical reactions and present a practical simulation scheme for simulating photochemical reactions in biomolecular systems. In our scheme, a multi-configurational quantum mechanical description is used to model the electronic rearrangement for those parts of the system that are involved in the photon absorption. For the remainder, typically consisting of the apo-protein and the solvent, a simple force field model is used. The interactions in the systems are thus computed within a hybrid quantum/classical framework. Forces are calculated on-the-fly, and a diabatic surface hopping procedure is used to model the excited-state decay. To demonstrate how this method is used we review our studies on photoactivation of the photoactive yellow protein, a bacterial photoreceptor. We will show what information can be obtained from the simulations, and, by comparing to recent experimental findings, what the limitations of our simulations are.

We performed ab initio excited-state molecular dynamics simulations of a gas-phase photoexcited protonated Schiff base (C1-N2═C3-C4═C5-C6) to search for control mechanisms of its photoisomerization. The excited molecule twists by ∼90° around either the N2C3 bond or the C4C5 bond and relaxes to the ground electronic state through a conical intersection with either a trans or cis outcome. We show that a large initial distortion of several dihedral angles and a specific normal vibrational mode combining pyramidalization and double-bond twisting can lead to a preferential rotation of atoms around the C4C5 bond. We also show that selective pretwisting of several dihedral angles in the initial ground state thermal ensemble (by analogy to a protein pocket) can significantly increase the fraction of photoreactive (cis → trans) trajectories. We demonstrate that new ensembles with higher degrees of control over the photoisomerization reaction can be obtained by a computational directed evolution approach on the ensembles of molecules with the pretwisted geometries.

Erratum: "Development of a force field for zeolitic imidazolate framework-8 with structural flexibility" [J. Chem.

No abstract