Elementary processes associated with ionization of liquid water provide a framework for understanding radiation-matter interactions in chemistry and biology. Although numerous studies have been conducted on the dynamics of the hydrated electron, its partner arising from ionization of liquid water, H2O+, remains elusive. We used tunable femtosecond soft x-ray pulses from an x-ray free electron laser to reveal the dynamics of the valence hole created by strong-field ionization and to track the primary proton transfer reaction giving rise to the formation of OH. The isolated resonance associated with the valence hole (H2O+/OH) enabled straightforward detection. Molecular dynamics simulations revealed that the x-ray spectra are sensitive to structural dynamics at the ionization site. We found signatures of hydrated-electron dynamics in the x-ray spectrum.

X-ray free-electron lasers enable the investigation of the structure and dynamics of diverse systems, including atoms, molecules, nanocrystals and single bioparticles, under extreme conditions. Many imaging applications that target biological systems and complex materials use hard X-ray pulses with extremely high peak intensities (exceeding 10 watts per square centimetre). However, fundamental investigations have focused mainly on the individual response of atoms and small molecules using soft X-rays with much lower intensities. Studies with intense X-ray pulses have shown that irradiated atoms reach a very high degree of ionization, owing to multiphoton absorption, which in a heteronuclear molecular system occurs predominantly locally on a heavy atom (provided that the absorption cross-section of the heavy atom is considerably larger than those of its neighbours) and is followed by efficient redistribution of the induced charge. In serial femtosecond crystallography of biological objects-an application of X-ray free-electron lasers that greatly enhances our ability to determine protein structure-the ionization of heavy atoms increases the local radiation damage that is seen in the diffraction patterns of these objects and has been suggested as a way of phasing the diffraction data. On the basis of experiments using either soft or less-intense hard X-rays, it is thought that the induced charge and associated radiation damage of atoms in polyatomic molecules can be inferred from the charge that is induced in an isolated atom under otherwise comparable irradiation conditions. Here we show that the femtosecond response of small polyatomic molecules that contain one heavy atom to ultra-intense (with intensities approaching 10 watts per square centimetre), hard (with photon energies of 8.3 kiloelectronvolts) X-ray pulses is qualitatively different: our experimental and modelling results establish that, under these conditions, the ionization of a molecule is considerably enhanced compared to that of an individual heavy atom with the same absorption cross-section. This enhancement is driven by ultrafast charge transfer within the molecule, which refills the core holes that are created in the heavy atom, providing further targets for inner-shell ionization and resulting in the emission of more than 50 electrons during the X-ray pulse. Our results demonstrate that efficient modelling of X-ray-driven processes in complex systems at ultrahigh intensities is feasible.

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.

We present a theoretical investigation of x-ray multiphoton ionization dynamics of polyatomic molecules, based on the rate equation model and molecular electronic structure calculations. An efficient numerical procedure is developed to calculate photoionization cross sections, Auger rates, and fluorescence rates for all possible electronic multiple-hole configurations of molecules. We investigate the charge-state distribution of a water molecule after interaction with an intense x-ray pulse and discuss its dependence on the fluence and the pulse duration of the x-ray beam. Our results demonstrate that a water molecule exposed to an intense x-ray pulse is more ionized than what would be expected within the independent-atom picture.

Ultrafast proton migration and isomerization are key processes for acetylene and its ions. However, the mechanism for ultrafast isomerization of acetylene [HCCH]2+ to vinylidene [H2CC]2+ dication remains nebulous. Theoretical studies show a large potential barrier ( > 2 eV) for isomerization on low-lying dicationic states, implying picosecond or longer isomerization timescales. However, a recent experiment at a femtosecond X-ray free-electron laser suggests sub-100 fs isomerization. Here we address this contradiction with a complete theoretical study of the dynamics of acetylene dication produced by Auger decay after X-ray photoionization of the carbon atom K shell. We find no sub-100 fs isomerization, while reproducing the salient features of the time-resolved Coulomb imaging experiment. This work resolves the seeming contradiction between experiment and theory and also calls for careful interpretation of structural information from the widely applied Coulomb momentum imaging method.

In many cases fragmentation of molecules upon inner-shell ionization is very unspecific with respect to the initially localized ionization site. Often this finding is interpreted in terms of an equilibration of internal energy into vibrational degrees of freedom after Auger decay. We investigate the X-ray photofragmentation of ethyl trifluoroacetate upon core electron ionization at environmentally distinct carbon sites using photoelectron-photoion-photoion coincidence measurements and ab initio electronic structure calculations. For all four carbon ionization sites, the Auger decay weakens the same bonds and transfers the two charges to opposite ends of the molecule, which leads to a rapid dissociation into three fragments, followed by further fragmentation steps. The lack of site specificity is attributed to the character of the dicationic electronic states after Auger decay instead of a fast equilibration of internal energy.

Determination of the electronic structure of mass-selected transient molecular ions which can be considered as building blocks of biomolecules.

There has been considerable debate on the existence of a low-barrier hydrogen bond (LBHB) in the photoactive yellow protein (PYP). The debate was initially triggered by the neutron diffraction study of Yamaguchi et al. ( Proc. Natl. Acad. Sci., U. S. A. , 2009 , 106 , 440 - 444 ) who suggested a model in which a neutral Arg52 residue triggers the formation of the LBHB in PYP. Here, we present an alternative model that is consistent within the error margins of the Yamaguchi structure factors. The model explains an increased hydrogen bond length without nuclear quantum effects and for a protonated Arg52. We tested both models by calculations under crystal, solution, and vacuum conditions. Contrary to the common assumption in the field, we found that a single PYP in vacuum does not provide an accurate description of the crystal conditions but instead introduces strong artifacts, which favor a LBHB and a large H NMR chemical shift. Our model of the crystal environment was found to stabilize the two Arg52 hydrogen bonds and crystal water positions for the protonated Arg52 residue in free MD simulations and predicted an Arg52 pK upshift with respect to PYP in solution. The crystal and solution environments resulted in almost identical H chemical shifts that agree with NMR solution data. We also calculated the effect of the Arg52 protonation state on the LBHB in 3D nuclear equilibrium density calculations. Only the charged crystal structure in vacuum supports a LBHB if Arg52 is neutral in PYP at the previously reported level of theory ( J. Am. Chem. Soc. , 2014 , 136 , 3542 - 3552 ). We attribute the anomalies in the interpretation of the neutron data to a shift of the potential minimum, which does not involve nuclear quantum effects and is transferable beyond the Yamaguchi structure.