Simulating Attochemistry: Which Dynamics Method to Use?

Attochemistry aims to exploit the properties of coherent electronic wavepackets excited via attosecond pulses to control the formation of photoproducts. Such molecular processes can, in principle, be simulated with various nonadiabatic dynamics methods, yet the impact of the approximations underlying the methods is rarely assessed. The performances of widely used mixed quantum-classical approaches, Tully surface hopping, and classical Ehrenfest methods are evaluated against the high-accuracy DD-vMCG quantum dynamics. This comparison is conducted for the valence ionization of fluorobenzene. Analyzing the nuclear motion induced in the branching space of the nearby conical intersection, the results show that the mixed quantum-classical methods reproduce quantitatively the average motion of a quantum wavepacket when initiated on a single electronic state. However, they fail to properly capture the nuclear motion induced by an electronic wavepacket along the derivative coupling, the latter originating from the quantum electronic coherence property, key to attochemistry.

Nonadiabatic dynamics with classical trajectories: The problem of an initial coherent superposition of electronic states

Advances in coherent light sources and development of pump-probe techniques in recent decades have opened the way to study electronic motion in its natural time scale. When an ultrashort laser pulse interacts with a molecular target, a coherent superposition of electronic states is created and the triggered electron dynamics is coupled to the nuclear motion. A natural and computationally efficient choice to simulate this correlated dynamics is a trajectory-based method where the quantum-mechanical electronic evolution is coupled to a classical-like nuclear dynamics. These methods must approximate the initial correlated electron-nuclear state by associating an initial electronic wavefunction to each classical trajectory in the ensemble. Different possibilities exist that reproduce the initial populations of the exact molecular wavefunction when represented in a basis. We show that different choices yield different dynamics and explore the effect of this choice in Ehrenfest, surface hopping, and exact-factorization-based coupled-trajectory schemes in a one-dimensional two-electronic-state model system that can be solved numerically exactly. This work aims to clarify the problems that standard trajectory-based techniques might have when a coherent superposition of electronic states is created to initialize the dynamics, to discuss what properties and observables are affected by different choices of electronic initial conditions and to point out the importance of quantum-momentum-induced electronic transitions in coupled-trajectory schemes.

Competing quantum effects in heavy-atom tunnelling through conical intersections

Thermally activated chemical reactions are typically understood in terms of overcoming potential-energy barriers. However, standard rate theories break down in the presence of a conical intersection (CI) because these processes are inherently nonadiabatic, invalidating the Born-Oppenheimer approximation. Moreover, CIs give rise to intricate nuclear quantum effects such as tunnelling and the geometric phase, which are neglected by standard trajectory-based simulations and remain largely unexplored in complex molecular systems. We present new semiclassical transition-state theories based on an extension of golden-rule instanton theory to describe nonadiabatic tunnelling through CIs and thus provide an intuitive picture for the reaction mechanism. We apply the method in conjunction with first-principles electronic-structure calculations to the electron transfer in the bis(methylene)-adamantyl cation. Our study reveals a strong competition between heavy-atom tunnelling and geometric-phase effects.

Photo-ionization initiated differential ultrafast charge migration: impacts of molecular symmetries and tautomeric forms

Photo-ionization induced ultrafast electron dynamics is considered as a precursor for the slower nuclear dynamics associated with molecular dissociation. Here, using the ab initio multielectron wave-packet propagation method, we study the overall many-electron dynamics, triggered by ionizing the outer-valence orbitals of different tautomers for a prototype molecule with more than one symmetry element. From the time evolution of the initially created averaged hole density of each system, we identify distinctly different charge dynamics responses in the tautomers. We observe that the keto form shows a charge migration direction away from the nitrogen bonded with hydrogen, while in enol-U - away from oxygen bonded to hydrogen. Additionally, the dynamics following the ionization of molecular orbitals with different symmetries reveals that a' orbitals show a fast and highly delocalized charge density in comparison to a'' symmetry. These observations indicate why different tautomers respond differently to an XUV ionization, and might explain the subsequent different fragmentation pathways. An experimental schematics allowing the detection and reconstruction of such charge dynamics is also proposed. Although the present study uses a simple, prototypical bio-relevant molecule, it reveals the explicit role of molecular symmetry and tautomerism in the ionization-triggered charge migration that controls many ultrafast physical, chemical, and biological processes, making tautomeric forms a promising tool of molecular design for desired charge migration.

Attochemistry: Is Controlling Electrons the Future of Photochemistry?

Controlling matter with light has always been a great challenge, leading to the ever-expanding field of photochemistry. In addition, since the first generation of light pulses of attosecond (1 as = 10^-18 s) duration, a great deal of effort has been devoted to observing and controlling electrons on their intrinsic time scale. Because of their short duration, attosecond pulses have a large spectral bandwidth populating several electronically excited states in a coherent manner, i.e., an electronic wavepacket. Because of interference, such a wavepacket has a new electronic distribution implying a potentially different and totally new reactivity as compared to traditional photochemistry, leading to the novel concept of "attochemistry". This nascent field requires the support of theory right from the start. In this Perspective, we discuss the opportunities offered by attochemistry, the related challenges, and the current and future state-of-the-art developments in theoretical chemistry needed to model it accurately.

Quantum electronic coherences by attosecond transient absorption spectroscopy: ab initio B-spline RCS-ADC study

Here I present a fully ab initio time-resolved study of X-ray attosecond transient absorption spectroscopy (ATAS) in a prototypical polyatomic molecule, pyrazine, and demonstrate the possibility of retrieving the many-electron quantum ionic coherences arising in attosecond molecular photoionisation and pre-determining the subsequent charge-directed photochemical reactivity. Advanced first-principles many-electron simulations are performed, within a hybrid XUV pump/X-ray probe setup, to describe the interaction of pyrazine with both XUV pump and X-ray probe pulses, and study the triggered correlated many-electron dynamics. The calculations are carried out by means of the recently-developed ab initio method for many-electron dynamics in polyatomic molecules, the time-dependent (TD) B-spline Restricted Correlation Space-Algebraic Diagrammatic Construction (RCS-ADC). RCS-ADC simulates molecular ionisation from first principles, combining the accurate description of electron correlation of quantum chemistry with the full account of the continuum dynamics of the photoelectron. Complete theoretical characterisation of the atto-ionised many-electron state and photo-induced attosecond charge dynamics is achieved by calculating the reduced ionic density matrix (R-IDM) of the bipartite ion-photoelectron system, with full inclusion of the correlated shakeup states. Deviations from the sudden approximation picture of photoionisation, in the low-photon-energy limit, are presented. The effect of the multi-channel interaction between the parent-ion and the emitted photoelectron on the onset of the quantum electronic coherences is analysed. Moreover, I show how the Schmidt decomposition of the R-IDM unravels the many-electron dynamics triggered by the pump, allowing for the identification of the key channels involved. Finally, I calculate the X-ray attosecond transient absorption spectra of XUV-ionised pyrazine. The results unveil the mapping of the ATAS measurement onto the quantum electronic coherences, and related non-zero R-IDM matrix elements, produced upon ionisation by the XUV pump laser pulse.

Control of nuclear dynamics in the benzene cation by electronic wavepacket composition

The study of coupled electron-nuclear dynamics driven by coherent superpositions of electronic states is now possible in attosecond science experiments. The objective is to understand the electronic control of chemical reactivity. In this work we report coherent 8-state non-adiabatic electron-nuclear dynamics simulations of the benzene radical cation. The computations were inspired by the extreme ultraviolet (XUV) experimental results in which all 8 electronic states were prepared with significant population. Our objective was to study the nuclear dynamics using various bespoke coherent electronic state superpositions as initial conditions in the Quantum-Ehrenfest method. The original XUV measurements were supported by Multi-configuration time-dependent Hartree (MCTDH) simulations, which suggested a model of successive passage through conical intersections. The present computations support a complementary model where non-adiabatic events are seen far from a conical intersection and are controlled by electron dynamics involving non-adjacent adiabatic states. It proves to be possible to identify two superpositions that can be linked with two possible fragmentation paths.

The "Hole" Story in Ionized Water from the Perspective of Ehrenfest Dynamics

The radiolysis of liquid water and the radiation-matter interactions that happen in aqueous environments are important to the fields of chemistry, materials, and environmental sciences, as well as the biological and physiological response to extreme conditions and medical treatments. The initial stage of radiolysis is the ultrafast response, or hole dynamics, that triggers chemical processes within complex energetic landscapes that may include reactivity. A fundamental understanding necessitates the use of theoretical methods that are capable of simulating both ultrafast coherence and non-adiabatic energy transfer pathways. In this work, we carry out an ab initio Ehrenfest dynamics study to provide a more complete description of the ultrafast dynamics and reactive events initiated by photoionization of water. After sudden ionization, a range of processes, including hole trapping and transfer, large OH oscillations, proton transfer and subsequent relay, formation of the metastable Zundel complex, and long-lived coherence, are identified and new insight into their driving forces is elucidated.

Onset of ionic coherence and ultrafast charge dynamics in attosecond molecular ionisation

Here is presented a fully ab initio theoretical framework for simulating the correlated many-electron dynamics occurring during and emerging from molecular ionisation by attosecond laser pulses. This is based on the time-dependent (TD) version of the B-spline restricted correlation space (RCS)-algebraic diagrammatic construction (ADC) method, with the full description of the photoelectron and inclusion of electron correlation effects, such as shakeup processes and inter-channel couplings. The nature of the ultrafast charge dynamics in the molecular ion is elucidated by quantitatively predicting the degree of electronic coherence and eigenstate content of the prepared molecular cationic state, beyond the commonly used sudden approximation. The results presented here for the acetylene and ethylene molecules show that even in the high photon energy regime the simulated hole dynamics is quantitatively different from the prediction of the sudden approximation. Moreover, for high-bandwidth ionising pulse, the residual interaction between the cation, in highly-excited shake-up states, and the emitted slow photoelectron gives rise to a loss of coherence in the ionic system which can persist for the first few femtoseconds after ionisation.

The Ehrenfest method with fully quantum nuclear motion (Qu-Eh): Application to charge migration in radical cations

An algorithm is described for quantum dynamics where an Ehrenfest potential is combined with fully quantum nuclear motion (Quantum-Ehrenfest, Qu-Eh). The method is related to the single-set variational multi-configuration Gaussian approach (vMCG) but has the advantage that only a single quantum chemistry computation is required at each time step since there is only a single time-dependent potential surface. Also shown is the close relationship to the "exact factorization method." The quantum Ehrenfest method is compared with vMCG for study of electron dynamics in a modified bismethylene-adamantane cation system. Illustrative examples of electron-nuclear dynamics are presented for a distorted allene system and for HCCI⁺ where one has a degenerate Π system.

Electron Dynamics upon Ionization of Polyatomic Molecules: Coupling to Quantum Nuclear Motion and Decoherence

Knowledge about the electronic motion in molecules is essential for our understanding of chemical reactions and biological processes. The advent of attosecond techniques opens up the possibility to induce electronic motion, observe it in real time, and potentially steer it. A fundamental question remains the factors influencing electronic decoherence and the role played by nuclear motion in this process. Here, we simulate the dynamics upon ionization of the polyatomic molecules paraxylene and modified bismethylene-adamantane, with a quantum mechanical treatment of both electron and nuclear dynamics using the direct dynamics variational multiconfigurational Gaussian method. Our simulations give new important physical insights about the expected decoherence process. We have shown that the decoherence of electron dynamics happens on the time scale of a few femtoseconds, with the interplay of different mechanisms: the dephasing is responsible for the fast decoherence while the nuclear overlap decay may actually help maintain it and is responsible for small revivals.