The second-order Ehrenfest method

Controlling Electronic Coherences and the Curvature Induced by the Derivative Coupling at a Conical Intersection: A Quantum Ehrenfest (QuEh) Protocol for Reaction Path Following Application to "Channel 3" Benzene Photochemistry

We report a protocol for the implementation of "reaction path following" from a transition state through a conical intersection, including both the path curvature induced by the derivative coupling and the corresponding induced electronic coherences. This protocol focuses on the "central" Gaussian wavepacket (initially unexcited) in the quantum Ehrenfest (QuEh) method. Like the reaction path following, the normal mode corresponding to the imaginary frequency at the transition state is given an initial momentum. The protocol is applied to the "channel 3" radiationless decay of benzene. We also demonstrate that one can enhance the effect of the derivative coupling and the electronic coherence with an IR pulse.

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.

Prediction of fluorescence quantum yields using the extended thawed Gaussian approximation

Spontaneous emission and internal conversion rates are calculated within harmonic approximations and compared to the results obtained within the semi-classical extended thawed Gaussian approximation (ETGA). This is the first application of the ETGA in the calculation of internal conversion and emission rates for real molecular systems, namely, formaldehyde, fluorobenzene, azulene, and a dicyano-squaraine dye. The viability of the models as black-box tools for prediction of spontaneous emission and internal conversion rates is assessed. All calculations were done using a consistent protocol in order to investigate how different methods perform without previous experimental knowledge using density functional theory (DFT) and time-dependent DFT (TD-DFT) with B3LYP, PBE0, ωB97XD, and CAM-B3LYP functionals. Contrasting the results with experimental data shows that there are further improvements required before theoretical predictions of emission and internal conversion rates can be used as reliable indicators for the photo-luminescence properties of molecules. We find that the ETGA performs rather similar to the vertical harmonical model. Including anharmonicities in the calculation of internal conversion rates has a moderate effect on the quantitative results in the studied systems. The emission rates are fairly stable with respect to computational parameters, but the internal conversion rate reveals itself to be highly dependent on the choice of the spectral line shape function, particularly the width of the Lorentzian function, associated with homogeneous broadening.

Coherent Mixing of Singlet and Triplet States in Acrolein and Ketene: A Computational Strategy for Simulating the Electron-Nuclear Dynamics of Intersystem Crossing

We present a theoretical study of intersystem crossing (ISC) in acrolein and ketene with the Ehrenfest method that can describe a superposition of singlet and triplet states. Our simulations illustrate a new mechanistic effect of ISC, namely, that a superposition of singlets and triplets yields nonadiabatic dynamics characteristic of that superposition rather than the constituent state potential energy surfaces. This effect is particularly significant in ketene, where mixing of singlet and triplet states along the approach to a singlet/singlet conical intersection occurs, with the spin-orbit coupling (SOC) remaining small throughout. In both cases, the effects require many recrossings of the singlet/triplet state crossing seam, consistent with the textbook treatment of ISC.

Nonadiabatic Molecular Dynamics Simulations Based on Time-Dependent Density Functional Tight-Binding Method

Nonadiabatic excited-state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamic (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. The nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Further, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well that obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.

Incorporating Lindblad Decay Dynamics into Mixed Quantum-Classical Simulations

We derive the $\mathcal{L}$-MFE method to incorporate Lindblad jump operator dynamics into the mean-field Ehrenfest (MFE) approach. We map the density matrix evolution of Lindblad dynamics onto pure state coefficients using trajectory averages. We use simple assumptions to construct the $\mathcal{L}$-MFE method that satisfies this exact mapping. This establishes a method that uses independent trajectories which exactly reproduces Lindblad decay dynamics using a wavefunction description, with deterministic changes of the magnitudes of the quantum expansion coefficients, while only adding on a stochastic phase. We further demonstrate that when including nuclei in the Ehrenfest dynamics, the $\mathcal{L}$-MFE method gives semi-quantitatively accurate results, with the accuracy limited by the accuracy of the approximations present in the semiclassical MFE approach. This work provides a general framework to incorporate Lindblad dynamics into semiclassical or mixed quantum-classical simulations.

How electronic superpositions drive nuclear motion following the creation of a localized hole in the glycine radical cation

In this work, we have studied the nuclear and electron dynamics in the glycine cation starting from localized hole states using the quantum Ehrenfest method. The nuclear dynamics is controlled both by the initial gradient and by the instantaneous gradient that results from the oscillatory electron dynamics (charge migration). We have used the Fourier transform (FT) of the spin densities to identify the “normal modes” of the electron dynamics. We observe an isomorphic relationship between the electron dynamics normal modes and the nuclear dynamics, seen in the vibrational normal modes. The FT spectra obtained this way show bands that are characteristic of the energy differences between the adiabatic hole states. These bands contain individual peaks that are in one-to-one correspondence with atom pair (+·) ↔ (·+) resonances, which, in turn, stimulate nuclear motion involving the atom pair. With such understanding, we anticipate “designer” coherent superpositions that can drive nuclear motion in a particular direction.

Simulations of molecular photodynamics in long timescales

Nonadiabatic dynamics simulations in the long timescale (much longer than 10 ps) are the next challenge in computational photochemistry. This paper delimits the scope of what we expect from methods to run such simulations: they should work in full nuclear dimensionality, be general enough to tackle any type of molecule and not require unrealistic computational resources. We examine the main methodological challenges we should venture to advance the field, including the computational costs of the electronic structure calculations, stability of the integration methods, accuracy of the nonadiabatic dynamics algorithms and software optimization. Based on simulations designed to shed light on each of these issues, we show how machine learning may be a crucial element for long time-scale dynamics, either as a surrogate for electronic structure calculations or aiding the parameterization of model Hamiltonians. We show that conventional methods for integrating classical equations should be adequate to extended simulations up to 1 ns and that surface hopping agrees semiquantitatively with wave packet propagation in the weak-coupling regime. We also describe our optimization of the Newton-X program to reduce computational overheads in data processing and storage. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.

High-order geometric integrators for representation-free Ehrenfest dynamics

Ehrenfest dynamics is a useful approximation for ab initio mixed quantum-classical molecular dynamics that can treat electronically nonadiabatic effects. Although a severe approximation to the exact solution of the molecular time-dependent Schrödinger equation, Ehrenfest dynamics is symplectic, is time-reversible, and conserves exactly the total molecular energy as well as the norm of the electronic wavefunction. Here, we surpass apparent complications due to the coupling of classical nuclear and quantum electronic motions and present efficient geometric integrators for "representation-free" Ehrenfest dynamics, which do not rely on a diabatic or adiabatic representation of electronic states and are of arbitrary even orders of accuracy in the time step. These numerical integrators, obtained by symmetrically composing the second-order splitting method and exactly solving the kinetic and potential propagation steps, are norm-conserving, symplectic, and time-reversible regardless of the time step used. Using a nonadiabatic simulation in the region of a conical intersection as an example, we demonstrate that these integrators preserve the geometric properties exactly and, if highly accurate solutions are desired, can be even more efficient than the most popular non-geometric integrators.

Unlocking the Double Bond in Protonated Schiff Bases by Coherent Superposition of S1 and S2

The primary event occurring during the E-to-Z photoisomerization reaction of retinal protonated Schiff base (rPSB) is single-to-double bond inversion. In this work we examine the nuclear dynamics that occurs when the initial excited state is a superposition of the S₁ and S₂ electronic excited states that might be created in a laser experiment. The nuclear dynamics is dominated by double bond inversion that is parallel to the derivative coupling vector of S₁ and S₂. Thus, the molecule behaves as if it were at a conical intersection even if the states are nondegenerate.

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.

Grid-Based Ehrenfest Model To Study Electron-Nuclear Processes

The two-dimensional electron-nuclear Schrödinger equation using soft-core Coulomb potentials has been a cornerstone for modeling and predicting the behavior of one-active-electron diatomic molecules, particularly for processes where both bound and continuum states are important. The model, however, is computationally expensive to extend to more electron or nuclear coordinates. Here we propose use of the Ehrenfest approach to treat the nuclear motion, while the electronic motion is still solved by quantum propagation on a grid. In this work, we present results for a one-dimensional treatment of H₂⁺, where the quantum and semiclassical dynamics can be directly compared, showing remarkably good agreement for a variety of situations. The advantage of the Ehrenfest approach is that it can be easily extended to treat as many nuclear degrees of freedom as needed.

Nonadiabatic Quantum Molecular Dynamics in Dense Manifolds of Electronic States

Most nonadiabatic molecular dynamics methods require the determination of a basis of adiabatic or diabatic electronic states at every time step, but in dense manifolds of electronic states, such approaches become intractable. A notable exception is Ehrenfest molecular dynamics, which can be implemented without explicit determination of such a basis but suffers from unphysical behavior when propagation on a mean-field potential energy surface (PES) does not accurately reflect the true dynamics on multiple electronic states. Here we introduce the multiple cloning for dense manifolds of states (MCDMS) method, a systematically improvable approximation to the multiple cloning method. MCDMS avoids both the mean-field PES problem and the need to compute the full electronic spectrum. This is achieved by reformulating multiple cloning to use a subspace of approximate eigenstates constructed from the time-dependent Ehrenfest electronic wave function. By application to model systems, we show that this approach allows a substantial reduction in the size of the required electronic basis without significant loss in accuracy.

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.

Computational Methods for Ab Initio Molecular Dynamics

Ab initio molecular dynamics is an irreplaceable technique for the realistic simulation of complex molecular systems and processes from first principles. This paper proposes a comprehensive and self-contained review of ab initio molecular dynamics from a computational perspective and from first principles. Quantum mechanics is presented from a molecular dynamics perspective. Various approximations and formulations are proposed, including the Ehrenfest, Born–Oppenheimer, and Hartree–Fock molecular dynamics. Subsequently, the Kohn–Sham formulation of molecular dynamics is introduced as well as the afferent concept of density functional. As a result, Car–Parrinello molecular dynamics is discussed, together with its extension to isothermal and isobaric processes. Car–Parrinello molecular dynamics is then reformulated in terms of path integrals. Finally, some implementation issues are analysed, namely, the pseudopotential, the orbital functional basis, and hybrid molecular dynamics.

Assessment of Real-Time Time-Dependent Density Functional Theory (RT-TDDFT) in Radiation Chemistry: Ionized Water Dimer

Ionization in the condensed phase and molecular clusters leads to a complicated chain of processes with coupled electron-nuclear dynamics. It is difficult to describe such dynamics with conventional nonadiabatic molecular dynamics schemes since the number of states swiftly increases as the molecular system grows. It is therefore attractive to use a direct electron and nuclear propagation such as the real-time time-dependent density functional theory (RT-TDDFT). Here we report a RT-TDDFT benchmark study on simulations of singly and doubly ionized states of a water monomer and dimer as a prototype for more complex processes in a condensed phase. We employed the RT-TDDFT based Ehrenfest molecular dynamics with a generalized gradient approximate (GGA) functional and compared it with wave-function-based surface hopping (SH) simulations. We found that the initial dynamics of a singly HOMO ionized water dimer is similar for both the RT-TDDFT/GGA and the SH simulations but leads to completely different reaction channels on a longer time scale. This failure is attributed to the self-interaction error in the GGA functionals and it can be avoided by using hybrid functionals with large fraction of exact exchange (represented here by the BHandHLYP functional). The simulations of doubly ionized states are reasonably described already at the GGA level. This suggests that the RT-TDDFT/GGA method could describe processes following the autoionization processes such as Auger emission, while its applicability to more complex processes such as intermolecular Coulombic decay remains limited.

Ultrafast charge dynamics in glycine induced by attosecond pulses

The combination of attosecond pump-probe techniques with mass spectrometry methods has recently led to the first experimental demonstration of ultrafast charge dynamics in a biomolecule, the amino acid phenylalanine [Calegari et al., Science, 2014, 346, 336]. Using an extension of the static-exchange density functional theory (DFT) method, the observed dynamics was explained as resulting from the coherent superposition of ionic states produced by the broadband attosecond pulse. Here, we have used the static-exchange DFT method to investigate charge migration induced by attosecond pulses in the glycine molecule. We show that the observed dynamics follows patterns similar to those previously found in phenylalanine, namely that charge fluctuations occur all over the molecule and that they can be explained in terms of a few typical frequencies of the system. We have checked the validity of our approach by explicitly comparing with the photoelectron spectra obtained in synchrotron radiation experiments and with the charge dynamics that follows the removal of an electron from a given molecular orbital, for which fully correlated ab initio results are available in the literature. From this comparison, we conclude that our method provides an accurate description of both the coherent superposition of cationic states generated by the attosecond pulse and its subsequent time evolution. Hence, we expect that the static-exchange DFT method should perform equally well for other medium-size and large molecules, for which the use of fully correlated ab initio methods is not possible.

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.

A corpuscular picture of electrons in chemical bond

We introduce a theory of chemical bond with a corpuscular picture of electrons. It employs a minimal set of localized electron wave packets with "floating and breathing" degrees of freedom and the spin-coupling of non-orthogonal valence-bond theory. Its accuracy for describing potential energy curves of chemical bonds in ground and excited states of spin singlet and triplet is examined.

Electron and nuclear dynamics following ionisation of modified bismethylene-adamantane

We have simulated the coupled electron and nuclear dynamics using the Ehrenfest method upon valence ionisation of modified bismethylene-adamantane (BMA) molecules where there is an electron transfer between the two π bonds. We have shown that the nuclear motion significantly affects the electron dynamics after a few fs when the electronic states involved are close in energy. We have also demonstrated how the non-stationary electronic wave packet determines the nuclear motion, more precisely the asymmetric stretching of the two π bonds, illustrating “charge-directed reactivity”. Taking into account the nuclear wave packet width results in the dephasing of electron dynamics with a half-life of 8 fs; this eventually leads to the equal delocalisation of the hole density over the two methylene groups and thus symmetric bond lengths.

Geometric Rotation of the Nuclear Gradient at a Conical Intersection: Extension to Complex Rotation of Diabatic States

Nonadiabatic dynamics in the vicinity of conical intersections is of essential importance in photochemistry. It is well known that if the branching space is represented in polar coordinates, then for a geometry represented by angle θ, the corresponding adiabatic states are obtained from the diabatic states with the mixing angle θ/2. In an equivalent way, one can study the relation between the real rotation of diabatic states and the resulting nuclear gradient. In this work, we extend the concept to allow a complex rotation of diabatic states to form a nonstationary superposition of electronic states. Our main result is that this leads to an elliptical transformation of the effective potential energy surfaces; i.e., the magnitude of the initial nuclear gradient changes as well as its direction. We fully explore gradient changes that result from varying both θ and ϕ (the complex rotation angle) as a way of electronically controlling nuclear motion, through Ehrenfest dynamics simulations for benzene cation.