Ionization-induced proton and energy transfer in liquid water

We report computational simulation results addressing the ionization response of liquid water upon valence ionization. The simulations cover ionizations in the whole valence-orbital range of liquid water, i.e., vacancies in 1b1, 3a1, 1b2, and 2a1 orbitals. It is found that ionization in any of these valence orbitals leads to rapid proton-transfer dynamics. The timescale on which the proton transfer occurs depends on which type of orbital is ionized. For ionization in the 2a1 orbitals, the proton transfer takes place in about 22 fs, competing with the intermolecular Coulombic decay mechanism that takes place on a similar timescale. This result is discussed in the context of earlier experimental results (Richter et al., Nat. Commun. 9, 4988) regarding the intermolecular Coulombic decay in water. For ionization in the outer-valence orbitals (1b1, 3a1, 1b2), we see rapid internal conversion via non-adiabatic transitions to the electronic ground state. The proton transfer occurs 46, 70, and 91 fs after the initial ionization from a 1b1, 3a1, and 1b2 orbital, respectively. The initial valence ionization induces strong vibrational excitations in the surrounding water molecules, leading to a considerable increase in the local effective temperature. The created heat diffuses into the liquid environment on a timescale of several hundred femtoseconds. We compare the results using two different embedding schemes, subtractive and electrostatic embedding, and find overall very similar dynamics.

Ultrafast time-resolved x-ray absorption spectroscopy of ionized urea and its dimer through ab initio nonadiabatic dynamics

Investigating the early dynamics of chemical systems following ionization is essential for our understanding of radiation damage. However, experimental as well as theoretical investigations are very challenging due to the complex nature of these processes. Time-resolved x-ray absorption spectroscopy on a femtosecond timescale, in combination with appropriate simulations, is able to provide crucial insights into the ultrafast processes that occur upon ionization due to its element-specific probing nature. In this theoretical study, we investigate the ultrafast dynamics of valence-ionized states of urea and its dimer employing Tully's fewest switches surface hopping approach using Koopmans' theorem to describe the ionized system. We demonstrate that following valence ionization through a pump pulse, the time-resolved x-ray absorption spectra at the carbon, nitrogen, and oxygen K-edges reveal rich insights into the dynamics. Excited states of the ionized system give rise to time-delayed blueshifts in the x-ray absorption spectra as a result of electronic relaxation dynamics through nonadiabatic transitions. Moreover, our statistical analysis reveals specific structural dynamics in the molecule that induce time-dependent changes in the spectra. For the urea monomer, we elucidate the possibility to trace effects of specific molecular vibrations in the time-resolved x-ray absorption spectra. For the urea dimer, where ionization triggers a proton transfer reaction, we show how the x-ray absorption spectra can reveal specific details on the progress of proton transfer.

Non-adiabatic quantum dynamics without potential energy surfaces based on second-quantized electrons: Application within the framework of the MCTDH method

A first principles quantum formalism to describe the non-adiabatic dynamics of electrons and nuclei based on a second quantization representation (SQR) of the electronic motion combined with the usual representation of the nuclear coordinates is introduced. This procedure circumvents the introduction of potential energy surfaces and non-adiabatic couplings, providing an alternative to the Born-Oppenheimer approximation. An important feature of the molecular Hamiltonian in the mixed first quantized representation for the nuclei and the SQR representation for the electrons is that all degrees of freedom, nuclear positions and electronic occupations, are distinguishable. This makes the approach compatible with various tensor decomposition Ansätze for the propagation of the nuclear-electronic wavefunction. Here, we describe the application of this formalism within the multi-configuration time-dependent Hartree framework and its multilayer generalization, corresponding to Tucker and hierarchical Tucker tensor decompositions of the wavefunction, respectively. The approach is applied to the calculation of the photodissociation cross section of the HeH⁺ molecule under extreme ultraviolet irradiation, which features non-adiabatic effects and quantum interferences between the two possible fragmentation channels, He + H⁺ and He⁺ + H. These calculations are compared with the usual description based on ab initio potential energy surfaces and non-adiabatic coupling matrix elements, which fully agree. The proof-of-principle calculations serve to illustrate the advantages and drawbacks of this formalism, which are discussed in detail, as well as possible ways to overcome them. We close with an outlook of possible application domains where the formalism might outperform the usual approach, for example, in situations that combine a strong static correlation of the electrons with non-adiabatic electronic-nuclear effects.

Opportunities for Two-Color Experiments in the Soft X-ray Regime at the European XFEL

X-ray pump/X-ray probe applications are made possible at X-ray Free Electron Laser (XFEL) facilities by generating two X-ray pulses with different wavelengths and controllable temporal delay. In order to enable this capability at the European XFEL, an upgrade project to equip the soft X-ray SASE3 beamline with a magnetic chicane is underway. In the present paper we describe the status of the project, its scientific focus and expected performance, including start-to-end simulations of the photon beam transport up to the sample, as well as recent experimental results demonstrating two-color lasing at photon energies of 805 eV + 835 eV and 910 eV + 950 eV. Additionally, we discuss methods to analyze the spectral properties and the intensity of the generated radiation to provide on-line diagnostics for future user experiments.

Dyson orbitals within the fc-CVS-EOM-CCSD framework: theory and application to X-ray photoelectron spectroscopy of ground and excited states

Ionization energies and Dyson orbitals within frozen-core core–valence separated equation-of-motion coupled cluster singles and doubles (fc-CVS-EOM-CCSD) enable efficient and reliable calculations of standard XPS and of UV-pump/XPS probe spectra.

Spin adapted implementation of EOM-CCSD for triplet excited states: Probing intersystem crossings of acetylacetone at the carbon and oxygen K-edges

We present an equation of motion coupled cluster singles and doubles approach for computing transient absorption spectra from a triplet excited state. The implementation determines the left and right excitation vectors by explicitly spin-adapting the triplet excitation space. As an illustrative application, we compute transient state X-ray absorption spectra at the carbon and oxygen K-edges for the acetylacetone molecule.

Hole dynamics in a photovoltaic donor-acceptor couple revealed by simulated time-resolved X-ray absorption spectroscopy

Theoretical and experimental methodologies that can characterize electronic and nuclear dynamics, and the coupling between the two, are needed to understand photoinduced charge transfer in molecular building blocks used in organic photovoltaics. Ongoing developments in ultrafast pump-probe techniques such as time-resolved X-ray absorption spectroscopy, using an X-ray free electron laser in combination with an ultraviolet femtosecond laser, present desirable probes of coupled electronic and nuclear dynamics. In this work, we investigate the charge transfer dynamics of a donor-acceptor pair, which is widely used as a building block in low bandgap block copolymers for organic photovoltaics. We simulate the dynamics of the benzothiadiazole-thiophene molecule upon photoionization with a vacuum ultraviolet (VUV) pulse and study the potential of probing the subsequent charge dynamics using time-resolved X-ray absorption spectroscopy. The photoinduced dynamics are calculated using on-the-fly nonadiabatic molecular dynamics simulations based on Tully's Fewest Switches Surface Hopping approach. We calculate the X-ray absorption spectrum as a function of time after ionization at the Hartree-Fock level. The changes in the time-resolved X-ray absorption spectrum at the sulfur K-edge reveal the ultrafast charge carrier dynamics in the molecule occurring on a femtosecond time scale. These theoretical findings anticipate that ultrafast time-resolved X-ray absorption spectroscopy using an X-ray probe in combination with a VUV pump offers a new approach to investigate the detailed dynamics of organic photovoltaic materials.

New and Efficient Equation-of-Motion Coupled-Cluster Framework for Core-Excited and Core-Ionized States

We present a fully analytical implementation of the core-valence separation (CVS) scheme for the equation-of-motion (EOM) coupled-cluster singles and doubles (CCSD) method for calculations of core-level states. Inspired by the CVS idea as originally formulated by Cederbaum, Domcke, and Schirmer, pure valence excitations are excluded from the EOM target space and the frozen-core approximation is imposed on the reference-state amplitudes and multipliers. This yields an efficient, robust, practical, and numerically balanced EOM-CCSD framework for calculations of excitation and ionization energies as well as state and transition properties (e.g., spectral intensities, natural transition, and Dyson orbitals) from both the ground and excited states. The errors in absolute excitation/ionization energies relative to the experimental reference data are on the order of 0.2-3.0 eV, depending on the K-edge considered and on the basis set used, and the shifts are systematic for each edge. Compared to a previously proposed CVS scheme where CVS was applied as a posteriori projection only during the solution of the EOM eigenvalue equations, the new scheme is computationally cheaper. It also achieves better cancellation of errors, yielding similar spectral profiles but with absolute core excitation and ionization energies that are systematically closer to the corresponding experimental data. Among the presented results are calculations of transient-state X-ray absorption spectra, relevant for interpretation of UV-pump/X-ray probe experiments.

Challenges in XUV Photochemistry Simulations: A Case Study on Ultrafast Fragmentation Dynamics of the Benzene Radical Cation

The challenges of simulating extreme ultraviolet (XUV)-induced dissociation dynamics of organic molecules on a multitude of coupled potential energy surfaces are discussed for the prototypical photoionization of benzene. The prospects of Koopmans' theorem-based electronic structure calculations in combination with classical trajectories and Tully's fewest switches surface hopping are explored. It is found that a Koopmans' theorem-based approach overestimates the CH dissociation barrier and thus underestimates the fragmentation yield. However, the nonadiabatic population dynamics are in good agreement with previous approaches, indicating that the Koopmans' theorem based potentials are well described around the Franck-Condon point. This is explicitly tested for the ground state potential of the benzene cation employing CASPT2 calculations, for which very good agreement is found. This work highlights the need for efficient electronic structure approaches that can treat medium-sized organic molecules with a multitude of coupled excited states and several dissociation channels.

Correlated proton-electron hole dynamics in protonated water clusters upon extreme ultraviolet photoionization

The ultrafast nuclear and electronic dynamics of protonated water clusters H(+)(H2O) n after extreme ultraviolet photoionization is investigated. In particular, we focus on cluster cations with n = 3, 6, and 21. Upon ionization, two positive charges are present in the cluster related to the excess proton and the missing electron, respectively. A correlation is found between the cluster's geometrical conformation and initial electronic energy with the size of the final fragments produced. For situations in which the electron hole and proton are initially spatially close, the two entities become correlated and separate in a time-scale of 20 to 40 fs driven by strong non-adiabatic effects.