Exact factorization of the time-dependent electron-nuclear wave function

Machine Learning for Electronically Excited States of Molecules

Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.

Nuclear-electronic orbital methods: Foundations and prospects

The incorporation of nuclear quantum effects and non-Born-Oppenheimer behavior into quantum chemistry calculations and molecular dynamics simulations is a longstanding challenge. The nuclear-electronic orbital (NEO) approach treats specified nuclei, typically protons, quantum mechanically on the same level as the electrons with wave function and density functional theory methods. This approach inherently includes nuclear delocalization and zero-point energy in molecular energy calculations, geometry optimizations, reaction paths, and dynamics. It can also provide accurate descriptions of excited electronic, vibrational, and vibronic states as well as nuclear tunneling and nonadiabatic dynamics. Nonequilibrium nuclear-electronic dynamics simulations beyond the Born-Oppenheimer approximation can be used to investigate a wide range of excited state processes. This Perspective provides an overview of the foundational NEO methods and enumerates the prospects for using these methods as building blocks for future developments. The conceptual simplicity and computational efficiency of the NEO approach will enhance its accessibility and applicability to diverse chemical and biological systems.

Ab initio molecular dynamics with screened Lorentz forces. I. Calculation and atomic charge interpretation of Berry curvature

The dynamics of a molecule in a magnetic field is significantly different from its zero-field counterpart. One important difference in the presence of a field is the Lorentz force acting on the nuclei, which can be decomposed as the sum of the bare nuclear Lorentz force and a screening force due to the electrons. This screening force is calculated from the Berry curvature and can change the dynamics qualitatively. It is therefore important to include the contributions from the Berry curvature in molecular dynamics simulations in a magnetic field. In this work, we present a scheme for calculating the Berry curvature numerically using a finite-difference technique, addressing challenges related to the arbitrary global phase of the wave function. The Berry curvature is calculated as a function of bond distance for H₂ at the restricted and unrestricted Hartree-Fock levels of theory and for CH⁺ as a function of the magnetic field strength at the restricted Hartree-Fock level of theory. The calculations are carried out using basis sets of contracted Gaussian functions equipped with London phase factors (London orbitals) to ensure gauge-origin invariance. In this paper, we also interpret the Berry curvature in terms of atomic charges and discuss its convergence in basis sets with and without London phase factors. The calculation of the Berry curvature allows for its inclusion in ab initio molecular dynamics simulations in a magnetic field.

Study of the Decoherence Correction Derived from the Exact Factorization Approach for Nonadiabatic Dynamics

We present a detailed study of the decoherence correction to surface hopping that was recently derived from the exact factorization approach. Ab initio multiple spawning calculations that use the same initial conditions and the same electronic structure method are used as a reference for three molecules: ethylene, the methaniminium cation, and fulvene, for which nonadiabatic dynamics follows a photoexcitation. A comparison with the Granucci–Persico energy-based decoherence correction and the augmented fewest-switches surface-hopping scheme shows that the three decoherence-corrected methods operate on individual trajectories in a qualitatively different way, but the results averaged over trajectories are similar for these systems.

PyUNIxMD: A Python-based excited state molecular dynamics package

Theoretical/computational description of excited state molecular dynamics is nowadays a crucial tool for understanding light-matter interactions in many materials. Here we present an open-source Python-based nonadiabatic molecular dynamics program package, namely PyUNIxMD, to deal with mixed quantum-classical dynamics for correlated electron-nuclear propagation. The PyUNIxMD provides many interfaces for quantum chemical calculation methods with commercial and noncommercial ab initio and semiempirical quantum chemistry programs. In addition, the PyUNIxMD offers many nonadiabatic molecular dynamics algorithms such as fewest-switch surface hopping and its derivatives as well as decoherence-induced surface hopping based on the exact factorization (DISH-XF) and coupled-trajectory mixed quantum-classical dynamics (CTMQC) for general purposes. Detailed structures and flows of PyUNIxMD are explained for the further implementations by developers. We perform a nonadiabatic molecular dynamics simulation for a molecular motor system as a simple demonstration.

Molecular vibrational frequencies from analytic Hessian of constrained nuclear-electronic orbital density functional theory

Nuclear quantum effects are important in a variety of chemical and biological processes. The constrained nuclear-electronic orbital density functional theory (cNEO-DFT) has been developed to include nuclear quantum effects in energy surfaces. Herein, we develop the analytic Hessian for cNEO-DFT energy with respect to the change in nuclear (expectation) positions, which can be used to characterize stationary points on energy surfaces and compute molecular vibrational frequencies. This is achieved by constructing and solving the multicomponent cNEO coupled-perturbed Kohn-Sham (cNEO-CPKS) equations, which describe the response of electronic and nuclear orbitals to the displacement of nuclear (expectation) positions. With the analytic Hessian, the vibrational frequencies of a series of small molecules are calculated and compared to those from conventional DFT Hessian calculations as well as those from the vibrational second-order perturbation theory (VPT2). It is found that even with a harmonic treatment, cNEO-DFT significantly outperforms DFT and is comparable to DFT-VPT2 in the description of vibrational frequencies in regular polyatomic molecules. Furthermore, cNEO-DFT can reasonably describe the proton transfer modes in systems with a shared proton, whereas DFT-VPT2 often faces great challenges. Our results suggest the importance of nuclear quantum effects in molecular vibrations, and cNEO-DFT is an accurate and inexpensive method to describe molecular vibrations.

Light-matter quantum dynamics of complex laser-driven systems

We propose a novel general approximation to transform and simplify the description of a complex fully quantized system describing the interacting light and matter. The method has some similarities to the time-dependent Born-Oppenheimer approach: we consider a quantum description of light rather than of nuclei and follow a similar separation procedure. Our approximation allows us to obtain a decoupled system for the light-excited matter and "dressed" light connected parametrically. With these equations at hand, we study how intense light as a quantum state is affected due to the back-action of the interacting matter. We discuss and demonstrate the possibility of the light-mode entanglement and nonclassical light generation during the interaction.

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.

On-the-Fly Symmetrical Quasi-Classical Dynamics with Meyer-Miller Mapping Hamiltonian for the Treatment of Nonadiabatic Dynamics at Conical Intersections

The on-the-fly version of the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian (SQC/MM) is implemented to study the nonadiabatic dynamics at conical intersections of polyatomic systems. The current on-the-fly implementation of the SQC/MM method is based on the adiabatic representation and the dressed momentum. To include the zero-point energy (ZPE) correction of the electronic mapping variables, we employ both the γ-adjusted and γ-fixed approaches. Nonadiabatic dynamics of the methaniminium cation (CH₂NH₂⁺) and azomethane are simulated using the on-the-fly SQC/MM method. For CH₂NH₂⁺, both ZPE correction approaches give reasonable and consistent results. However, for azomethane, the γ-adjusted version of the SQC/MM dynamics behaves much better than the γ-fixed version. Further analysis indicates that it is always recommended to use the γ-adjusted SQC/MM dynamics in the on-the-fly simulation of photoinduced dynamics of polyatomic systems, particularly when the excited state is well separated from the ground state in the Franck-Condon region. This work indicates that the on-the-fly SQC/MM method is a powerful simulation protocol to deal with the nonadiabatic dynamics of realistic polyatomic systems.

Signatures of Conical Intersection Dynamics in the Time-Resolved Photoelectron Spectrum of Furan: Theoretical Modeling with an Ensemble Density Functional Theory Method

The non-adiabatic dynamics of furan excited in the ππ* state (S2 in the Franck-Condon geometry) was studied using non-adiabatic molecular dynamics simulations in connection with an ensemble density functional method. The time-resolved photoelectron spectra were theoretically simulated in a wide range of electron binding energies that covered the valence as well as the core electrons. The dynamics of the decay (rise) of the photoelectron signal were compared with the excited-state population dynamics. It was observed that the photoelectron signal decay parameters at certain electron binding energies displayed a good correlation with the events occurring during the excited-state dynamics. Thus, the time profile of the photoelectron intensity of the K-shell electrons of oxygen (decay constant of 34 ± 3 fs) showed a reasonable correlation with the time of passage through conical intersections with the ground state (47 ± 2 fs). The ground-state recovery constant of the photoelectron signal (121 ± 30 fs) was in good agreement with the theoretically obtained excited-state lifetime (93 ± 9 fs), as well as with the experimentally estimated recovery time constant (ca. 110 fs). Hence, it is proposed to complement the traditional TRPES observations with the trXPS (or trNEXAFS) measurements to obtain more reliable estimates of the most mechanistically important events during the excited-state dynamics.

Validating fewest-switches surface hopping in the presence of laser fields

The capability of fewest-switches surface hopping (FSSH) to describe non-adiabatic dynamics under explicit excitation with external fields is evaluated. Different FSSH parameters are benchmarked against multi-configurational time dependent Hartree (MCTDH) reference calculations using SO₂ and 2-thiocytosine as model, yet realistic, molecular systems. Qualitatively, FSSH is able to reproduce the trends in the MCTDH dynamics with (also without) an explicit external field; however, no set of FSSH parameters is ideal. The adequate treatment of the overcoherence in FSSH is revealed as the driving factor to improve the description of the excitation process with respect to the MCTDH reference. Here, two corrections were tested: the augmented-FSSH (AFSSH) correction and the energy-based decoherence correction. A dependence on the employed basis is detected in AFSSH, performing better when spin-orbit and external laser field couplings are treated as off-diagonal elements instead of projecting them onto the diagonal of the Hamilton operator. In the presence of an electric field, the excited state dynamics was found to depend strongly on the vector used to rescale the kinetic energy along after a transition between surfaces. For SO₂, recurrence of the excited wave packet throughout the duration of the applied laser pulse is observed for laser pulses (>100 fs), resulting in additional interferences missed by FSSH and only visible in variational multi-configurational Gaussian when utilizing a large number of Gaussian basis functions. This feature vanishes when going toward larger molecules, such as 2-thiocytosine, where this effect is barely visible in a laser pulse 200 fs long.

Coupled nuclear and electron dynamics in the vicinity of a conical intersection

Ultrafast optical techniques allow us to study ultrafast molecular dynamics involving both nuclear and electronic motion. To support interpretation, theoretical approaches are needed that can describe both the nuclear and electron dynamics. Hence, we revisit and expand our ansatz for the coupled description of the nuclear and electron dynamics in molecular systems (NEMol). In this purely quantum mechanical ansatz, the quantum-dynamical description of the nuclear motion is combined with the calculation of the electron dynamics in the eigenfunction basis. The NEMol ansatz is applied to simulate the coupled dynamics of the molecule NO₂ in the vicinity of a conical intersection (CoIn) with a special focus on the coherent electron dynamics induced by the non-adiabatic coupling. Furthermore, we aim to control the dynamics of the system when passing the CoIn. The control scheme relies on the carrier envelope phase of a few-cycle IR pulse. The laser pulse influences both the movement of the nuclei and the electrons during the population transfer through the CoIn.

Charge-Transfer Steps in Density Functional Theory from the Perspective of the Exact Electron Factorization

When a molecule dissociates, the exact Kohn-Sham (KS) and Pauli potentials may form step structures. Reproducing these steps correctly is central for the description of dissociation and charge-transfer processes in density functional theory (DFT): The steps align the KS eigenvalues of the dissociating subsystems relative to each other and determine where electrons localize. While the step height can be calculated from the asymptotic behavior of the KS orbitals, this provides limited insight into what causes the steps. We give an explanation of the steps with an exact mapping of the many-electron problem to a one-electron problem, the exact electron factorization (EEF). The potentials appearing in the EEF have a clear physical meaning that translates to the DFT potentials by replacing the interacting many-electron system with the KS system. With a simple model of a diatomic, we illustrate that the steps are a consequence of spatial electron entanglement and are the result of a charge transfer. From this mechanism, the step height can immediately be deduced. Moreover, two methods to approximately reproduce the potentials during dissociation are proposed. One is based on the states of the dissociated system, while the other one is based on an analogy to the Born-Oppenheimer treatment of a molecule. The latter method also shows that the steps connect adiabatic potential energy surfaces. The view of DFT from the EEF thus provides a better understanding of how many-electron effects are encoded in a one-electron theory and how they can be modeled.

How important are the residual nonadiabatic couplings for an accurate simulation of nonadiabatic quantum dynamics in a quasidiabatic representation?

Diabatization of the molecular Hamiltonian is a standard approach to remove the singularities of nonadiabatic couplings at conical intersections of adiabatic potential energy surfaces. In general, it is impossible to eliminate the nonadiabatic couplings entirely-the resulting "quasidiabatic" states are still coupled by smaller but nonvanishing residual nonadiabatic couplings, which are typically neglected. Here, we propose a general method for assessing the validity of this potentially drastic approximation by comparing quantum dynamics simulated either with or without the residual couplings. To make the numerical errors negligible to the errors due to neglecting the residual couplings, we use the highly accurate and general eighth-order composition of the implicit midpoint method. The usefulness of the proposed method is demonstrated on nonadiabatic simulations in the cubic Jahn-Teller model of nitrogen trioxide and in the induced Renner-Teller model of hydrogen cyanide. We find that, depending on the system, initial state, and employed quasidiabatization scheme, neglecting the residual couplings can result in wrong dynamics. In contrast, simulations with the exact quasidiabatic Hamiltonian, which contains the residual couplings, always yield accurate results.

Quantum-classical nonadiabatic dynamics of Floquet driven systems

We develop a trajectory-based approach for excited-state molecular dynamics simulations of systems subject to an external periodic drive. We combine the exact-factorization formalism, allowing us to treat electron-nuclear systems in nonadiabatic regimes, with the Floquet formalism for time-periodic processes. The theory is developed starting with the molecular time-dependent Schrödinger equation with the inclusion of an external periodic drive that couples to the system dipole moment. With the support of the Floquet formalism, quantum dynamics is approximated by combining classical-like, trajectory-based, nuclear evolution with electronic dynamics represented in the Floquet basis. The resulting algorithm, which is an extension of the coupled-trajectory mixed quantum-classical scheme for periodically driven systems, is applied to a model study, exactly solvable, with different field intensities.

Recent progress in approximate quantum dynamics methods for the study of proton-coupled electron transfer reactions

Proton-coupled electron transfer (PCET) reactions are ubiquitous natural processes at the heart of energy conversion reactions in photosynthesis and respiration, DNA repair, and diverse enzymatic reactions. Theoretical formulation and computational method developments have eyed modeling of thermal and photoinduced PCET for the last three decades. The accumulation of these studies, collected in dozens of reviews, accounts, and perspectives, has firmly established the influence of quantum effects, including non-adiabatic electronic transitions, vibrational relaxation, zero-point energy, and proton tunneling, on the rate and mechanism of PCET reactions. Here, we focus on some recently-developed methods, spanning the last eight years, that can quantitatively capture these effects in the PCET context and provide efficient means for their qualitative description in complex systems. The theoretical background of each method and their accuracy with respect to exact results are discussed and the results of relevant PCET simulations based on each method are presented.

Relaxation dynamics through a conical intersection: Quantum and quantum-classical studies

We study the relaxation process through a conical intersection of a photo-excited retinal chromophore model. The analysis is based on a two-electronic-state two-dimensional Hamiltonian developed by Hahn and Stock [J. Phys. Chem. B 104 1146 (2000)] to reproduce, with a minimal model, the main features of the 11-cis to all-trans isomerization of the retinal of rhodopsin. In particular, we focus on the performance of various trajectory-based schemes to nonadiabatic dynamics, and we compare quantum-classical results to the numerically exact quantum vibronic wavepacket dynamics. The purpose of this work is to investigate, by analyzing electronic and nuclear observables, how the sampling of initial conditions for the trajectories affects the subsequent dynamics.

Case studies of the time-dependent potential energy surface for dynamics in cavities

The exact time-dependent potential energy surface driving the nuclear dynamics was recently shown to be a useful tool to understand and interpret the coupling of nuclei, electrons, and photons in cavity settings. Here, we provide a detailed analysis of its structure for exactly solvable systems that model two phenomena: cavity-induced suppression of proton-coupled electron-transfer and its dependence on the initial state, and cavity-induced electronic excitation. We demonstrate the inadequacy of simply using a weighted average of polaritonic surfaces to determine the dynamics. Such a weighted average misses a crucial term that redistributes energy between the nuclear and the polaritonic systems, and this term can in fact become a predominant term in determining the nuclear dynamics when several polaritonic surfaces are involved. Evolving an ensemble of classical trajectories on the exact potential energy surface reproduces the nuclear wavepacket quite accurately, while evolving on the weighted polaritonic surface fails after a short period of time. The implications and prospects for application of mixed quantum-classical methods based on this surface are discussed.

The role of electron-nuclear coupling on multi-state photoelectron spectra, scattering processes and phase transitions

We present first principle based beyond Born-Oppenheimer (BBO) theory and its applications on various models as well as realistic spectroscopic and scattering processes, where the Jahn-Teller (JT) theory is brought in conjunction with the BBO approach on the phase transition of lanthanide complexes. Over one and half decades, our development of BBO theory is demonstrated with ab initio calculations on representative molecules of spectroscopic interest (NO2 radical, Na3 and K3 clusters, NO3 radical, C6H6+ and 1,3,5-C6H3F3+ radical cations) as well as triatomic reactive scattering processes (H+ + H2 and F + H2). Such an approach exhibits the effect of JT, Renner-Teller (RT) and pseudo Jahn-Teller (PJT) type of interactions. While implementing the BBO theory, we generate highly accurate diabatic potential energy surfaces (PESs) to carry out quantum dynamics calculation and find excellent agreement with experimental photoelectron spectra of spectroscopic systems and cross-sections/rate constants of scattering processes. On the other hand, such electron-nuclear couplings incorporated through JT theory play a crucial role in dictating higher energy satellite transitions in the dielectric function spectra of the LaMnO3 complex. Overall, this article thoroughly sketches the current perspective of the BBO approach and its connection with JT theory with various applications on physical and chemical processes.

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.

Quantum and Quantum-Classical Studies of the Photoisomerization of a Retinal Chromophore Model

We report an in-depth analysis of the photo-induced isomerization of the 2-cis-penta-2,4-dieniminium cation: a minimal model of the 11-cis retinal protonated Schiff base chromophore of the dim-light photoreceptor rhodopsin. Based on recently developed three-dimensional potentials parametrized on ab initio multi-state multi-configurational second-order perturbation theory data, we perform quantum-dynamical studies. In addition, simulations based on various quantum-classical methods, among which Tully surface hopping and the coupled-trajectory approach derived from the exact factorization, allow us to validate their performance against vibronic wavepacket propagation and, therefore, a purely quantum treatment. Quantum-dynamics results uncover qualitative differences with respect to the two-dimensional Hahn-Stock potentials, widely used as model potentials for the isomerization of the same chromophore, due to the increased dimensionality and three-mode correlation. Quantum-classical simulations show, instead, that three-dimensional model potentials are capable of capturing a number of features revealed by atomistic simulations and experimental observations. In particular, a recently reported vibrational phase relationship between double-bond torsion and hydrogen-out-of-plane modes critical for rhodopsin isomerization efficiency is correctly reproduced.

Molecular Photochemistry: Recent Developments in Theory

Photochemistry is a fascinating branch of chemistry that is concerned with molecules and light. However, the importance of simulating light-induced processes is reflected also in fields as diverse as biology, material science, and medicine. This Minireview highlights recent progress achieved in theoretical chemistry to calculate electronically excited states of molecules and simulate their photoinduced dynamics, with the aim of reaching experimental accuracy. We focus on emergent methods and give selected examples that illustrate the progress in recent years towards predicting complex electronic structures with strong correlation, calculations on large molecules, describing multichromophoric systems, and simulating non-adiabatic molecular dynamics over long time scales, for molecules in the gas phase or in complex biological environments.

Effect of many modes on self-polarization and photochemical suppression in cavities

The standard description of cavity-modified molecular reactions typically involves a single (resonant) mode, while in reality, the quantum cavity supports a range of photon modes. Here, we demonstrate that as more photon modes are accounted for, physicochemical phenomena can dramatically change, as illustrated by the cavity-induced suppression of the important and ubiquitous process of proton-coupled electron-transfer. Using a multi-trajectory Ehrenfest treatment for the photon-modes, we find that self-polarization effects become essential, and we introduce the concept of self-polarization-modified Born-Oppenheimer surfaces as a new construct to analyze dynamics. As the number of cavity photon modes increases, the increasing deviation of these surfaces from the cavity-free Born-Oppenheimer surfaces, together with the interplay between photon emission and absorption inside the widening bands of these surfaces, leads to enhanced suppression. The present findings are general and will have implications for the description and control of cavity-driven physical processes of molecules, nanostructures, and solids embedded in cavities.

Ultrafast Spectroscopy of Photoactive Molecular Systems from First Principles: Where We Stand Today and Where We Are Going

Computational spectroscopy is becoming a mandatory tool for the interpretation of the complex, and often congested, spectral maps delivered by modern non-linear multi-pulse techniques. The fields of Electronic Structure Methods, Non-Adiabatic Molecular Dynamics, and Theoretical Spectroscopy represent the three pillars of the virtual ultrafast optical spectrometer, able to deliver transient spectra in silico from first principles. A successful simulation strategy requires a synergistic approach that balances between the three fields, each one having its very own challenges and bottlenecks. The aim of this Perspective is to demonstrate that, despite these challenges, an impressive agreement between theory and experiment is achievable now regarding the modeling of ultrafast photoinduced processes in complex molecular architectures. Beyond that, some key recent developments in the three fields are presented that we believe will have major impacts on spectroscopic simulations in the very near future. Potential directions of development, pending challenges, and rising opportunities are illustrated.

Quantum Trajectories for the Dynamics in the Exact Factorization Framework: A Proof-of-Principle Test

In the framework of the exact factorization of the time-dependent electron-nuclear wave function, we investigate the possibility of solving the nuclear time-dependent Schrödinger equation based on trajectories. The nuclear equation is separated in a Hamilton-Jacobi equation for the phase of the wave function, and a continuity equation for its (squared) modulus. For illustrative adiabatic and nonadiabatic one-dimensional models, we implement a procedure to follow the evolution of the nuclear density along the characteristics of the Hamilton-Jacobi equation. Those characteristics are referred to as quantum trajectories, since they are generated via ordinary differential equations similar to Hamilton's equations, but including the so-called quantum potential, and they can be used to reconstruct exactly the quantum-mechanical nuclear wave function, provided infinite initial conditions are propagated in time.

Internal Conversion and Intersystem Crossing with the Exact Factorization

We present a detailed derivation of the generalized coupled-trajectory mixed quantum-classical (G-CT-MQC) algorithm based on the exact-factorization equations. The ultimate goal is to propose an algorithm that can be employed for molecular dynamics simulations of nonradiative phenomena, as the spin-allowed internal conversions and the spin-forbidden intersystem crossings. Internal conversions are nonadiabatic processes driven by the kinetic coupling between electronic states, whereas intersystem crossings are mediated by the spin-orbit coupling. In this paper, we discuss computational issues related to the suitable representation for electronic dynamics and the different natures of kinetic and spin-orbit coupling. Numerical studies on model systems allow us to test the performance of the G-CT-MQC algorithm in different situations.

On the proper derivation of the Floquet-based quantum classical Liouville equation and surface hopping describing a molecule or material subject to an external field

We investigate different approaches to derive the proper Floquet-based quantum-classical Liouville equation (F-QCLE) for laser-driven electron-nuclear dynamics. The first approach projects the operator form of the standard QCLE onto the diabatic Floquet basis and then transforms to the adiabatic representation. The second approach directly projects the QCLE onto the Floquet adiabatic basis. Both approaches yield a form that is similar to the usual QCLE with two modifications: (1) The electronic degrees of freedom are expanded to infinite dimension and (2) the nuclear motion follows Floquet quasi-energy surfaces. However, the second approach includes an additional cross derivative force due to the dual dependence on time and nuclear motion of the Floquet adiabatic states. Our analysis and numerical tests indicate that this cross derivative force is a fictitious artifact, suggesting that one cannot safely exchange the order of Floquet state projection with adiabatic transformation. Our results are in accord with similar findings by Izmaylov et al., [J. Chem. Phys. 140, 084104 (2014)] who found that transforming to the adiabatic representation must always be the last operation applied, although now we have extended this result to a time-dependent Hamiltonian. This paper and the proper derivation of the F-QCLE should lay the basis for further improvements of Floquet surface hopping.

The quantum-Ehrenfest method with the inclusion of an IR pulse: Application to electron dynamics of the allene radical cation

We describe the implementation of a laser control pulse in the quantum-Ehrenfest method, a molecular quantum dynamics method that solves the time-dependent Schrödinger equation for both electrons and nuclei. The oscillating electric field-dipole interaction is incorporated directly in the one-electron Hamiltonian of the electronic structure part of the algorithm. We then use the coupled electron-nuclear dynamics of the π-system in the allene radical cation (•CH₂=C=CH₂)⁺ as a simple model of a pump-control experiment. We start (pump) with a two-state superposition of two cationic states. The resulting electron dynamics corresponds to the rapid oscillation of the unpaired electron between the two terminal methylenes. This electron dynamics is, in turn, coupled to the torsional motion of the terminal methylenes. There is a conical intersection at 90° twist, where the electron dynamics collapses because the adiabatic states become degenerate. After passing the conical intersection, the electron dynamics revives. The IR pulse (control) in our simulations is timed to have its maximum at the conical intersection. Our simulations show that the effect of the (control) pulse is to change the electron dynamics at the conical intersection and, as a consequence, the concomitant nuclear dynamics, which is dominated by the change in the torsional angle.

Embedding via the Exact Factorization Approach

We present a quantum electronic embedding method derived from the exact factorization approach to calculate static properties of a many-electron system. The method is exact in principle but the practical power lies in utilizing input from a low-level calculation on the entire system in a high-level method computed on a small fragment, as in other embedding methods. Here, the exact factorization approach defines an embedding Hamiltonian on the fragment. Various Hubbard models demonstrate that remarkably accurate ground-state energies are obtained over the full range of weak to strongly correlated systems.

Modeling nonadiabatic dynamics in condensed matter materials: some recent advances and applications

This review focuses on recent developments in the field of nonadiabatic molecular dynamics (NA-MD), with particular attention given to condensed-matter systems. NA-MD simulations for small molecular systems can be performed using high-level electronic structure (ES) calculations, methods accounting for the quantization of nuclear motion, and using fewer approximations in the dynamical methodology itself. Modeling condensed-matter systems imposes many limitations on various aspects of NA-MD computations, requiring approximations at various levels of theory-from the ES, to the ways in which the coupling of electrons and nuclei are accounted for. Nonetheless, the approximate treatment of NA-MD in condensed-phase materials has gained a spin lately in many applied studies. A number of advancements of the methodology and computational tools have been undertaken, including general-purpose methods, as well as those tailored to nanoscale and condensed matter systems. This review summarizes such methodological and software developments, puts them into the broader context of existing approaches, and highlights some of the challenges that remain to be solved.

Nonadiabatic Dynamics in a Laser Field: Using Floquet Fewest Switches Surface Hopping To Calculate Electronic Populations for Slow Nuclear Velocities

We investigate two well-known approaches for extending the fewest switches surface hopping (FSSH) algorithm to periodic time-dependent couplings. The first formalism acts as if the instantaneous adiabatic electronic states were standard adiabatic states, which just happen to evolve in time. The second formalism replaces the role of the usual adiabatic states by the time-independent adiabatic Floquet states. For a set of modified Tully model problems, the Floquet FSSH (F-FSSH) formalism gives a better estimate for both transmission and reflection probabilities than the instantaneous adiabatic FSSH (IA-FSSH) formalism, especially for slow nuclear velocities. More importantly, only F-FSSH predicts the correct final scattering momentum. Finally, in order to use Floquet theory accurately, we find that it is crucial to account for the interference between wavepackets on different Floquet states. Our results should be of interest to all those interested in laser-induced molecular dynamics.

Spin-Orbit Interactions in Ultrafast Molecular Processes

We investigate spin-orbit interactions in ultrafast molecular processes employing the exact factorization of the electron-nuclear wave function. We revisit the original derivation by including spin-orbit coupling, and show how the dynamics driven by the time-dependent potential energy surface alleviates inconsistencies arising from different electronic representations. We propose a novel trajectory-based scheme to simulate spin-forbidden non-radiative processes, and we show its performance in the treatment of excited-state dynamics where spin-orbit effects couple different spin multiplets.

How to calculate charge mobility in molecular materials from surface hopping non-adiabatic molecular dynamics - beyond the hopping/band paradigm

Charge transport in high mobility organic semiconductors is in an intermediate regime between small polaron hopping and band transport limits. We have recently shown that surface hopping non-adiabatic molecular dynamics is a powerful method for prediction of charge transport mechanisms in organic materials and for near-quantitative prediction of charge mobilities at room temperature where the effects of nuclear zero-point motion and tunneling are still relatively small [S. Giannini et al., Nat. Commun., 2019, 10, 3843]. Here we assess and critically discuss the extensions to Tully's original method that have led to this success: (i) correction for missing electronic decoherence, (ii) detection of trivial crossings and (iii) removal of decoherence correction-induced spurious charge transfer. If any one of these corrections is not included, the charge mobility diverges with system size, each for different physical reasons. Yet if they are included, convergence with system size, detailed balance and good internal consistency are achieved.

Hovering States of Ammonia in a High-Intensity, High-Frequency Oscillating Field: Trapped into Planarity by Laser-Induced Hybridization

A high-intensity, high-frequency laser can create an oscillating induced dipole moment in a molecule. At high laser frequencies with a long pulse width, a stable non-ionizing state with a laser-induced hybridization of the electrons is formed. For ammonia, aligned with the linear polarization direction of the laser, such stable states can be realized. Electronic hybridization in the presence of the high-frequency field is such that the lone pair propensity is dynamically equalized on either side of ammonia. This leads to a destabilization of pyramidal ammonia and hovering states with the electron density flipping to either side of the geometry. Electronic structure calculations in an oscillating frame of reference anticipate this effect with a predicted classical quiver distance of 0.1 Å. Electronic dynamics at a laser intensity of 1.14 × 10¹³ W/cm² and a frequency of 8.16 eV predicts negligible ionization for the planar geometry. Approximate nuclear wave packet dynamics in the oscillating potential energy generated by the electrons predicts a trapping of ammonia in its planar transition state geometry.