Mean-field dynamics with stochastic decoherence (MF-SD): A new algorithm for nonadiabatic mixed quantum/classical molecular-dynamics simulations with nuclear-induced decoherence

Nonadiabatic Dynamics with Exact Factorization: Implementation and Assessment

In this work, we report our implementation of several independent-trajectory mixed-quantum-classical (ITMQC) nonadiabatic dynamics methods based on exact factorization (XF) in the Libra package for nonadiabatic and excited-state dynamics. Namely, the exact factorization surface hopping (SHXF), mixed quantum-classical dynamics (MQCXF), and mean-field (MFXF) are introduced. Performance of these methods is compared to that of several traditional surface hopping schemes, such as the fewest-switches surface hopping (FSSH), branching-corrected surface hopping (BCSH), and the simplified decay of mixing (SDM), as well as conventional Ehrenfest (mean-field, MF) method. Based on a comprehensive set of 1D model Hamiltonians, we find the ranking SHXF ≈ MQCXF > BCSH > SDM > FSSH ≫ MF, with the BCSH sometimes outperforming the XF methods in terms of describing coherences. Although the MFXF method can yield reasonable populations and coherences for some cases, it does not conserve the total energy and is therefore not recommended. We also find that the branching correction for auxiliary trajectories is important for the XF methods to yield accurate populations and coherences. However, the branching correction can worsen the quality of the energy conservation in the MQCXF. Finally, we find that using the time-dependent Gaussian width approximation used in the XF methods for computing decoherence correction can improve the quality of energy conservation in the MQCXF dynamics. The parameter-free scheme of Subotnik for computing the Gaussian widths is found to deliver the best performance in situations where such widths are not known a priori.

Using Machine Learning to Understand the Causes of Quantum Decoherence in Solution-Phase Bond-Breaking Reactions

Decoherence is a fundamental phenomenon that occurs when an entangled quantum state interacts with its environment, leading to collapse of the wave function. The inevitability of decoherence provides one of the most intrinsic limits of quantum computing. However, there has been little study of the precise chemical motions from the environment that cause decoherence. Here, we use quantum molecular dynamics simulations to explore the photodissociation of Na2+ in liquid Ar, in which solvent fluctuations induce decoherence and thus determine the products of chemical bond breaking. We use machine learning to characterize the solute-solvent environment as a high-dimensional feature space that allows us to predict when and onto which photofragment the bonding electron will localize. We find that reaching a requisite photofragment separation and experiencing out-of-phase solvent collisions underlie decoherence during chemical bond breaking. Our work highlights the utility of machine learning for interpreting complex solution-phase chemical processes as well as identifies the molecular underpinnings of decoherence.

Energy-Conserving and Thermally Corrected Neglect of Back-Reaction Approximation Method for Nonadiabatic Molecular Dynamics

In this work, the energy-conserving and thermally corrected neglect of the back-reaction approximation approach for nonadiabatic molecular dynamics in extended atomistic systems is developed. The new approach introduces three key corrections to the original method: (1) it enforces the total energy conservation, (2) it introduces an explicit coupling of the system to its environment, and (3) it introduces a renormalization of nonadiabatic couplings to account for a difference between the instantaneous nuclear kinetic energy and the kinetic energy of guiding trajectories. In the new approach, an auxiliary kinetic energy variable is introduced as an independent dynamical variable. The new approach produces nonzero equilibrium populations, whereas the original neglect of the back-reaction approximation method does not. It yields population relaxation time scales that are favorably comparable to the reference values, and it introduces an explicit and controllable way of dissipating energy into a bath without an assumption of the bath being at equilibrium.

A first principles derivation of energy-conserving momentum jumps in surface hopping simulations

The fewest switches surface hopping (FSSH) method proposed by Tully in 1990 [Tully, J. Chem. Phys. 93, 1061 (1990)]-along with its many later variations-forms the basis for most practical simulations of molecular dynamics with electronic transitions in realistic systems. Despite its popularity, a rigorous formal derivation of the algorithm has yet to be achieved. In this paper, we derive the energy-conserving momentum jumps employed by FSSH from the perspective of quantum trajectory surface hopping (QTSH) [Martens, J. Phys. Chem. A 123, 1110 (2019)]. In the limit of localized nonadiabatic transitions, simple mathematical and physical arguments allow the FSSH algorithm to be derived from first principles. For general processes, the quantum forces characterizing the QTSH method provide accurate results for nonadiabatic dynamics with rigorous energy conservation, at the ensemble level, within the consistency of the underlying stochastic surface hopping without resorting to the artificial momentum rescaling of FSSH.

Nonadiabatic simulations of photoisomerization and dissociation in ethylene using ab initio classical trajectories

We simulate the nonadiabatic dynamics of photo-induced isomerization and dissociation in ethylene using ab initio classical trajectories in an extended phase space of nuclear and electronic variables. This is achieved by employing the linearized semiclassical initial value representation method for nonadiabatic dynamics, where discrete electronic states are mapped to continuous classical variables using either the Meyer-Miller-Stock-Thoss representation or a more recently introduced spin mapping approach. Trajectory initial conditions are sampled by constraining electronic state variables to a single initial excited state and by drawing nuclear phase space configurations from a Wigner distribution at a finite temperature. An ensemble of classical ab initio trajectories is then generated to compute thermal population correlation functions and analyze the mechanisms of isomerization and dissociation. Our results serve as a demonstration that this parameter-free semiclassical approach is computationally efficient and accurate, identifying mechanistic pathways in agreement with previous theoretical studies and also uncovering dissociation pathways observed experimentally.

Zeno and Anti-Zeno Effects in Nonadiabatic Molecular Dynamics

Decoherence plays an important role in nonadiabatic (NA) molecular dynamics (MD) simulations because it provides a physical mechanism for trajectory hopping and can alter transition rates by orders of magnitude. Generally, decoherence effects slow quantum transitions, as exemplified by the quantum Zeno effect: in the limit of infinitely fast decoherence, the transitions stop. If the measurements are not sufficiently frequent, an opposite quantum anti-Zeno effect occurs, in which the transitions are accelerated with faster decoherence. Using two common NA-MD approaches, fewest switches surface hopping and decoherence-induced surface hopping, combined with analytic examination, we demonstrate that including decoherence into NA-MD slows down NA transitions; however, many realistic systems operate in the anti-Zeno regime. Therefore, it is important that NA-MD methods describe both Zeno and anti-Zeno effects. Numerical simulations of charge trapping and relaxation in graphitic carbon nitride suggest that time-dependent NA Hamiltonians encountered in realistic systems produce robust results with respect to errors in the decoherence time, a favorable feature for NA-MD simulations.

Nonadiabatic Dynamics in a Continuous Circularly Polarized Laser Field with Floquet Phase-Space Surface Hopping

Nonadiabatic chemical reactions involving continuous circularly polarized light (cw CPL) have not attracted as much attention as dynamics in unpolarized/linearly polarized light. However, including circularly (in contrast to linearly) polarized light allows one to effectively introduce a complex-valued time-dependent Hamiltonian, which offers a new path for control or exploration through the introduction of Berry forces. Here, we investigate several inexpensive semiclassical approaches for modeling such nonadiabatic dynamics in the presence of a time-dependent complex-valued Hamiltonian, beginning with a straightforward instantaneous adiabatic fewest-switches surface hopping (IA-FSSH) approach (where the electronic states depend on position and time), continuing to a standard Floquet fewest switches surface hopping (F-FSSH) approach (where the electronic states depend on position and frequency), and ending with an exotic Floquet phase-space surface hopping (F-PSSH) approach (where the electronic states depend on position, frequency, and momentum). Using a set of model systems with time-dependent complex-valued Hamiltonians, we show that the Floquet phase-space adiabats are the optimal choice of basis as far as accounting for Berry phase effects and delivering accuracy. Thus, the F-PSSH algorithm sets the stage for future modeling of nonadiabatic dynamics under strong externally pumped circular polarization.

Pedagogical Overview of the Fewest Switches Surface Hopping Method

The fewest switches surface hopping method continues to grow in popularity to capture electronic nonadiabaticity and quantum nuclear effects due to its simplicity and accuracy. Knowing the basics of the method is essential for the correct implementation and interpretation of results. This review covers the fundamentals of the fewest switches surface hopping method with a detailed discussion of the nuances such as decoherence schemes and frustrated hops and the correct approach to calculating populations. The consequences of incorrect implementation are further discussed toward calculating kinetic and thermodynamic properties. Some tips for practitioners and a step-by-step algorithm for developers are provided. Finally, some of the finer technicalities of the fewest switches surface hopping method that are buried deep in the literature are pointed out to help graduate students better appreciate this method.

Modeling Spin-Dependent Nonadiabatic Dynamics with Electronic Degeneracy: A Phase-Space Surface-Hopping Method

Nuclear Berry curvature effects emerge from electronic spin degeneracy and can lead to nontrivial spin-dependent (nonadiabatic) nuclear dynamics. However, such effects are not captured fully by any current mixed quantum-classical method such as fewest-switches surface hopping. In this work, we present a phase-space surface-hopping (PSSH) approach to simulate singlet-triplet intersystem crossing dynamics. We show that with a simple pseudodiabatic ansatz, a PSSH algorithm can capture the relevant Berry curvature effects and make predictions in agreement with exact quantum dynamics for a simple singlet-triplet model Hamiltonian. Thus, this approach represents an important step toward simulating photochemical and spin processes concomitantly, as relevant to intersystem crossing and spin-lattice relaxation dynamics.

Electronic Absorption Spectra from Off-Diagonal Quantum Master Equations

Quantum master equations (QMEs) provide a general framework for describing electronic dynamics within a complex molecular system. Off-diagonal QMEs (OD-QMEs) correspond to a family of QMEs that describe the electronic dynamics in the interaction picture based on treating the off-diagonal coupling terms between electronic states as a small perturbation within the framework of second-order perturbation theory. The fact that OD-QMEs are given in terms of the interaction picture makes it non-trivial to obtain Schrodinger picture electronic coherences from them. A key experimental quantity that relies on the ability to obtain accurate Schrodinger picture electronic coherences is the absorption spectrum. In this paper, we propose using a recently introduced procedure for extracting Schrodinger picture electronic coherences from interaction picture inputs to calculate electronic absorption spectra from electronic dynamics generated by OD-QMEs. The accuracy of the absorption spectra obtained in this way is studied in the context of a biexciton benchmark model, by comparing spectra calculated based on time-local and time-nonlocal OD-QMEs to spectra calculated based on a Redfield-type QME and the non-perturbative and quantum-mechanically exact hierarchical equations of motion (HEOM) method.

Multilayer Subsystem Surface Hopping Method for Large-Scale Nonadiabatic Dynamics Simulation with Hundreds of Thousands of States

We present a multilayer subsystem surface hopping (MSSH) method to deal with nonadiabatic dynamics in large-scale systems. A small subsystem instead of the full system is adopted for surface hopping and is updated on-the-fly to achieve a reliable description of important adiabatic states and the wave function evolution. Additional subsystems for molecular dynamics and statistical description are introduced to further improve the simulation reliability. The global flux hopping probabilities with optimal state assignments are utilized to treat the complex surface crossings. As demonstrated in a series of one- and two-dimensional Holstein models with up to hundreds of thousands of states, MSSH shows weak parameter dependence in all investigated systems. Especially, the computational costs are reduced by 2-6 orders of magnitude compared to traditional surface hopping simulations in full systems, and size-independent results are achieved with a large time-step size of 2-5 fs. The new method is compatible with different decoherence correction strategies and achieves a much better balance between efficiency and reliability, thus promising for applications in general charge and exciton dynamics simulations.

A Mixed Deterministic-Stochastic Algorithm of the Branching Corrected Mean Field Method for Nonadiabatic Dynamics

We present a new algorithm of the branching corrected mean field (BCMF) method for nonadiabatic dynamics [J. Xu and L. Wang, J. Phys. Chem. Lett. 11, 8283 (2020)], which combines the key advantages of the two existed algorithms, i.e., the deterministic BCMF algorithm based on weights of trajectory branches (BCMF-w) and the stochastic BCMF algorithm with random collapse of the electronic wavefunction (BCMF-s). The resulting mixed deterministic-stochastic BCMF algorithm (BCMF-ws) is benchmarked in a series of standard scattering problems with potential wells on the excited-state surfaces, which are common in realistic systems. In all investigated cases, BCMF-ws holds the same high accuracy while the computational time is reduced about two orders of magnitude compared to the original BCMF-w and BCMF-s algorithms, thus promising for nonadiabatic dynamics simulations of general systems.

Modeling Non-adiabatic Dynamics in Nanoscale and Condensed Matter Systems

Rapid, far-from-equilibrium processes involving excitation of electronic, vibrational, spin, photon, topological, and other degrees of freedom form the basis of modern technologies, including electronics and optoelectronics, solar energy harvesting and conversion to electrical and chemical energy, quantum information processing, spin- and valleytronics, chemical detection, and medical therapies. Such processes are studied experimentally with various time-resolved spectroscopies that allow scientists to track system's evolution on ultrafast time scales and at close to atomistic level of detail. The availability of various forms of lasing has made such measurements easily accessible to many experimental groups worldwide, to study atoms and small molecules, nanoscale and condensed matter systems, proteins, cells, and mesoscopic materials. The experimental work necessitates parallel theoretical efforts needed to interpret the experiments and to provide insights that cannot be gained through measurements due to experimental limitations.Non-adiabatic (NA) molecular dynamics (MD) allows one to study processes at the atomistic level and in the time domain most directly mimicking the time-resolved experiments. Atomistic modeling takes full advantage of chemical intuition and principles that guide design and fabrication of molecules and materials. It provides atomistic origins of quasi-particles, such as holes, excitons, trions, plasmons, phonons, polarons, polaritons, spin-waves, momentum-resolved and topological states, electrically and magnetically polarized structures, and other abstract concepts. An atomistic description enables one to study realistic aspects of materials, which necessarily contain defects, dopants, surfaces, interfaces, passivating ligands, and solvent layers. Often, such realistic features govern material properties and are hard to account for phenomenologically. NA-MD requires few approximations and assumptions. It does not need to assume that atomic motions are harmonic, that electrons are Drude oscillators, that coupling between different degrees of freedom is weak, that dynamics is Markovian or has short memory, or that evolution occurs by exponential kinetics of transitions between few states. The classical or semiclassical treatment of atomic motions constitutes the main approximation of NA-MD and is used because atoms are 3-5 orders of magnitude heavier than electrons. NA-MD is limited by system size, typically hundreds or thousands of atoms, and time scale, picoseconds to nanoseconds. The quality of NA-MD simulations depends on the electronic structure method used to obtain excited state energies and NA couplings.NA-MD has been largely popularized and advanced in the chemistry community that focuses on molecules. Modeling far-from-equilibrium dynamics in nanoscale and condensed matter systems often has to account for other types of physics. At the same time, condensed phase NA-MD allows for approximations that may not work in molecules. Focusing on the recent NA-MD developments aimed at studying excited state processes in nanoscale and condensed phases, this Account considers how the phenomena important on the nanoscale can be incorporated into NA-MD and what approximations can be made to increase its efficiency with complex systems and processes.

An accurate, non-empirical method for incorporating decoherence into Ehrenfest dynamics

In mixed quantum-classical nonadiabatic molecular dynamics methods, the anchoring of the electronic wave function to a single nuclear geometry results in both quantitative and qualitative errors in the dynamics. In the context of both Ehrenfest and trajectory surface hopping methods, methods for incorporating decoherence are widely used to eliminate these errors. However, the accuracy of these methods often depends strongly on the parameterization of the decoherence time and/or other related quantities. Here, we present a refinement of the recently introduced collapse to a block (TAB) scheme for incorporating decoherence into Ehrenfest dynamics. The proposed approach incorporates an approximation to the history of the population dynamics and treats the coherence decay as Gaussian, rather than exponential. This method uses parameters that can be obtained from first principles, rather than empirical fitting. Application to one-dimensional models indicates excellent agreement with numerically exact simulations. We also introduce a second refinement to the TAB method: a robust linear least-squares algorithm for determining collapse probabilities.

On simulating the dynamics of electronic populations and coherences via quantum master equations based on treating off-diagonal electronic coupling terms as a small perturbation

Quantum master equations provide a general framework for describing the dynamics of electronic observables within a complex molecular system. One particular family of such equations is based on treating the off-diagonal coupling terms between electronic states as a small perturbation within the framework of second-order perturbation theory. In this paper, we show how different choices of projection operators, as well as whether one starts out with the time-convolution or the time-convolutionless forms of the generalized quantum master equation, give rise to four different types of such off-diagonal quantum master equations (OD-QMEs), namely, time-convolution and time-convolutionless versions of a Pauli-type OD-QME for only the electronic populations and an OD-QME for the full electronic density matrix (including both electronic populations and coherences). The fact that those OD-QMEs are given in terms of the interaction picture makes it non-trivial to obtain Schrödinger picture electronic coherences from them. To address this, we also extend a procedure for extracting Schrödinger picture electronic coherences from interaction picture populations recently introduced by Trushechkin in the context of time-convolutionless Pauli-type OD-QME to the other three types of OD-QMEs. The performance of the aforementioned four types of OD-QMEs is explored in the context of the Garg-Onuchic-Ambegaokar benchmark model for charge transfer in the condensed phase across a relatively wide parameter range. The results show that time-convolution OD-QMEs can be significantly more accurate than their time-convolutionless counterparts, particularly in the case of Pauli-type OD-QMEs, and that rather accurate Schrödinger picture coherences can be obtained from interaction picture electronic inputs.

Nonadiabatic Dynamics with Coupled Trajectories

In this paper, we discuss coupled-trajectory schemes for molecular-dynamics simulations of excited-state processes. New coupled-trajectory strategies to capture decoherence effects, revival of coherence and nonadiabatic interferences in long-time dynamics are proposed, and compared to independent-trajectory schemes. The working framework is provided by the exact factorization of the electron-nuclear wave function, and it exploits ideas emanating from various surface-hopping schemes. The new coupled-trajectory algorithms are tested on a one-dimensional two-state system using different model parameters which allow one to induce different dynamics. The benchmark is provided by the numerically exact solution of the time-dependent Schrödinger equation.

Semiclassical description of nuclear dynamics moving through complex-valued single avoided crossings of two electronic states

The standard fewest-switches surface hopping (FSSH) approach fails to model nonadiabatic dynamics when the electronic Hamiltonian is complex-valued and there are multiple nuclear dimensions; FSSH does not include geometric magnetic effects and does not have access to a gauge independent direction for momentum rescaling. In this paper, for the case of a Hamiltonian with two electronic states, we propose an extension of Tully's FSSH algorithm, which includes geometric magnetic forces and, through diabatization, establishes a well-defined rescaling direction. When combined with a decoherence correction, our new algorithm shows satisfying results for a model set of two-dimensional single avoided crossings.

Branching and phase corrected surface hopping: A benchmark of nonadiabatic dynamics in multilevel systems

Since the seminal work of Tully [J. Chem. Phys. 93, 1061 (1990)], two-level scattering models have been extensively adopted as the standard benchmark systems to assess the performance of different trajectory surface hopping methods for nonadiabatic dynamics simulations. Here, we extend the branching and phase corrections to multilevel systems and combine them with both the traditional fewest switches surface hopping (FSSH) and its variant global flux surface hopping (GFSH) algorithms. To get a comprehensive evaluation of the proposed methods, we construct a series of more challenging and diverse three-level and four-level scattering models and use exact quantum solutions as references. Encouragingly, both FSSH and GFSH with the branching and phase corrections produce excellent and nearly identical results in all investigated systems, indicating that the new surface hopping methods are robust to describe multilevel problems and the reliability is insensitive to the definition of self-consistent hopping probabilities in the adiabatic representation. Furthermore, the branching correction is found to be especially important when dealing with strongly repulsive potential energy surfaces, which are common in realistic systems, thus promising for general applications.

Simulating energy transfer dynamics in the Fenna-Matthews-Olson complex via the modified generalized quantum master equation

The generalized quantum master equation (GQME) provides a general and formally exact framework for simulating the reduced dynamics of open quantum systems. The recently introduced modified approach to the GQME (M-GQME) corresponds to a specific implementation of the GQME that is geared toward simulating the dynamics of the electronic reduced density matrix in systems governed by an excitonic Hamiltonian. Such a Hamiltonian, which is often used for describing energy and charge transfer dynamics in complex molecular systems, is given in terms of diabatic electronic states that are coupled to each other and correspond to different nuclear Hamiltonians. Within the M-GQME approach, the effect of the nuclear degrees of freedom on the time evolution of the electronic density matrix is fully captured by a memory kernel superoperator, which can be obtained from short-lived (compared to the time scale of energy/charge transfer) projection-free inputs. In this paper, we test the ability of the M-GQME to predict the energy transfer dynamics within a seven-state benchmark model of the Fenna-Matthews-Olson (FMO) complex, with the short-lived projection-free inputs obtained via the Ehrenfest method. The M-GQME with Ehrenfest-based inputs is shown to yield accurate results across a wide parameter range. It is also found to dramatically outperform the direct application of the Ehrenfest method and to provide better-behaved convergence with respect to memory time in comparison to an alternative implementation of the GQME approach previously applied to the same FMO model.

Benchmarking the Surface Hopping Method to Include Nuclear Quantum Effects

We have benchmarked the surface hopping method to capture nuclear quantum effects in the spin-Boson model in the deep tunneling regime. The thermal populations and the rate constants calculated using the surface hopping method are compared with those calculated using Boltzmann theory and Fermi's golden rule, respectively. Additionally, we have proposed a simple kinetic model that partially includes nuclear quantum effects within Marcus theory, and the results of the surface hopping method are analyzed under the framework of this simple kinetic model. A broad range of parameters are investigated to identify the regimes for the successes and failures of the surface hopping method. This work shows that with the accurate treatment of decoherence and velocity reversal, surface hopping can generally capture the nuclear quantum effects in the deep tunneling and weak diabatic coupling regime.

Effect of varying the TD-lc-DFTB range-separation parameter on charge and energy transfer in a model pentacene/buckminsterfullerene heterojunction

Atomistic modeling of energy and charge transfer at the heterojunction of organic solar cells is an active field with many remaining outstanding questions owing, in part, to the difficulties in performing reliable photodynamics calculations on very large systems. One approach to being able to overcome these difficulties is to design and apply an appropriate simplified method. Density-functional tight binding (DFTB) has become a popular form of approximate density-functional theory based on a minimal valence basis set and neglect of all but two center integrals. We report the results of our tests of a recent long-range correction (lc) [A. Humeniuk and R. Mitrić, J. Chem. Phys. 143, 134120 (2015)] for time-dependent (TD) lc-DFTB by carrying out TD-lc-DFTB fewest switches surface hopping calculations of energy and charge transfer times using the relatively new DFTBABY [A. Humeniuk and R. Mitrić, Comput. Phys. Commun. 221, 174 (2017)] program. An advantage of this method is the ability to run enough trajectories to get meaningful ensemble averages. Our interest in the present work is less in determining exact energy and charge transfer rates than in understanding how the results of these calculations vary with the value of the range-separation parameter (R_lc = 1/μ) for a model organic solar cell heterojunction consisting of a gas-phase van der Waals complex P/F made up of a single pentacene (P) molecule together with a single buckminsterfullerene (F) molecule. The default value of R_lc = 3.03 a₀ is found to be much too small as neither energy nor charge transfer is observed until R_lc ≈ 10 a₀. Tests at a single geometry show that the best agreement with high-quality ab initio spectra is obtained in the limit of no lc (i.e., very large R_lc). A plot of energy and charge transfer rates as a function of R_lc is provided, which suggests that a value of R_lc ≈ 15 a₀ yields the typical literature (condensed-phase) charge transfer time of about 100 fs. However, energy and charge transfer times become as high as ∼300 fs for R_lc ≈ 25 a₀. A closer examination of the charge transfer process P^*/F → P⁺/F^- shows that the initial electron transfer is accompanied by a partial delocalization of the P hole onto F, which then relocalizes back onto P, consistent with a polaron-like picture in which the nuclei relax to stabilize the resultant redistribution of charges.

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.

Nonadiabatic dynamics simulation of photoinduced ring-opening reaction of 2(5H)-thiophenone with internal conversion and intersystem crossing

In the present work, the quantum trajectory mean-field approach, which is able to overcome the overcoherence problem, was generalized to simulate internal conversion and intersystem crossing processes simultaneously. The photoinduced ring-opening and subsequent rearrangement reactions of isolated 2(5H)-thiophenone were studied based on geometry optimizations on critical structures and nonadiabatic dynamics simulations using this method. Upon 267 nm irradiation, the molecule is initially populated in the 1ππ* state. After a sudden rupture of one C-S bond within 100 fs in this state, the lowest two singlet excited states and the lowest two triplet excited states become quasi-degenerated, and then the intersystem crossing processes between singlet and triplet states accompanied by rearrangement reactions can be observed several times. Compared with our previous nonadiabatic simulations in the absence of intersystem crossing (ChemPhotoChem, 2019, 3, 897-906), some new nonadiabatic relaxation pathways involving triplet states and different ring-opening products were identified. The present work provides new mechanistic insights into the photoinduced ring-opening of thio-substituted heterocyclic molecules and reveals the importance of nonadiabatic dynamics simulation that is able to deal with multiple electronic states with different spin multiplicities.

Evaluation of mixed quantum-classical molecular dynamics on cis-azobenzene photoisomerization

The quantitative prediction of nonadiabatic transitions between different electronic states is important to understand ultrafast processes in photochemistry. A variety of mixed quantum-classical molecular dynamics methods such as surface hopping and Ehrenfest mean-field have been developed. However, how to choose an appropriate one from a wide diversity of dynamics algorithms to study a realistic photochemical process is still unclear. In this work, we implemented 30 combinations of different mixed quantum-classical dynamics methods, including 24 surface hopping models with 8 decoherence corrections and 3 momentum rescaling strategies as well as 6 mean-field models. Then we performed numerical investigations by simulating the photoisomerization of cis-azobenzene combined with on-the-fly electronic structure calculations. Predictions of the S1 lifetime and the quantum yield of the photoproduct using different models are distinct. Surface hopping is more robust than mean-field in our test system. Moreover, the choice of momentum rescaling methods in surface hopping brings more significant changes than decoherence corrections, while a large discrepancy between simulation results with different mean-field algorithms has been observed.

Experimental and theoretical study of energy transfer in a chromophore triad: What makes modeling dynamics successful?

Simulation of electronic dynamics in realistically large molecular systems is a demanding task that has not yet achieved the same level of quantitative prediction already realized for its static counterpart. This is particularly true for processes occurring beyond the Born-Oppenheimer regime. Non-adiabatic molecular dynamics (NAMD) simulations suffer from two convoluted sources of error: numerical algorithms for dynamics and electronic structure calculations. While the former has gained increasing attention, particularly addressing the validity of ad hoc methodologies, the effect of the latter remains relatively unexplored. Indeed, the required accuracy for electronic structure calculations to reach quantitative agreement with experiment in dynamics may be even more strict than that required for static simulations. Here, we address this issue by modeling the electronic energy transfer in a donor-acceptor-donor (D-A-D) molecular light harvesting system using fewest switches surface hopping NAMD simulations. In the studied system, time-resolved experimental measurements deliver complete information on spectra and energy transfer rates. Subsequent modeling shows that the calculated electronic transition energies are "sufficiently good" to reproduce experimental spectra but produce over an order of magnitude error in simulated dynamical rates. We further perform simulations using artificially shifted energy gaps to investigate the complex relationship between transition energies and modeled dynamics to understand factors affecting non-radiative relaxation and energy transfer rates.

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.

PySpawn: Software for Nonadiabatic Quantum Molecular Dynamics

The ab initio multiple spawning (AIMS) method enables nonadiabatic quantum molecular dynamics simulations in an arbitrary number of dimensions, with potential energy surfaces provided by electronic structure calculations performed on-the-fly. However, the intricacy of the AIMS algorithm complicates software development, deployment on modern shared computer resources, and postsimulation data analysis. PySpawn is a nonadiabatic molecular dynamics software package that addresses these issues. The program is designed to be easily interfaced with electronic structure software, and an interface to the TeraChem software package is described here. PySpawn introduces a task-based reorganization of the AIMS algorithm, allowing fine-grained restart capability and setting the stage for efficient parallelization in a future release. PySpawn includes a user-friendly and interactive Python analysis module that will enable novice users to painlessly adopt AIMS. As a demonstration of PySpawn's simulation capability and analysis module, we report complete active space self-consistent field-based AIMS simulations of the 1,2-dithienyl-1,2-dicyanoethene molecule, a promising molecular photoswitch.

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.

State-pairwise decoherence times for nonadiabatic dynamics on more than two electronic states

Independent trajectory (IT) nonadiabatic molecular dynamics simulation methods are powerful tools for modeling processes involving transitions between electronic states. Incorporation and refinement of decoherence corrections into popular IT methods, e.g., Ehrenfest dynamics and trajectory surface hopping, is an important means of improving their accuracies. In this work, we identify a new challenge in the development of such decoherence corrections; when a system exists in a coherent superposition of three or more electronic states, coherences may decay unphysically when the decoherence correction is based on decoherence times assigned on a state-wise basis. As a solution, we introduce decoherence corrected Ehrenfest schemes based on decoherence times assigned on a state-pairwise basis. By application of these methods to a set of very simple one-dimensional model problems, we show that one of these state-pairwise methods ("collapse to a block") correctly describes the loss of coherence between all pairs of states in our multistate model problems, whereas a method based on a state-wise description of coherence loss does not. The new one-dimensional models introduced here can serve as useful tests for other decoherence correction schemes.

Identification of a receiver triplet state in the ultrafast intersystem crossing of carbonylpyrenes

The near-unity triplet quantum yield of photoexcited carbonyl functionalized pyrenes is theoretically investigated. The estimated energetics of singlet-triplet manifolds and relevant spin-orbit coupling parameters strongly suggest triplet state formation via the S₁→ T₄/T₅ pathway. Quantum wavepacket dynamics of triplet manifolds within the linear vibronic coupling approach reveal that the receiver triplet state would undergo rapid internal conversion decay to the lower triplet state(s), facilitating efficient triplet generation by minimizing the reverse intersystem crossing possibilities. On the basis of these results, a unified mechanism is proposed to describe the ultrafast intersystem crossing process in these molecules.

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.

Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials

Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.

A molecular perspective on Tully models for nonadiabatic dynamics

We present a series of standardized molecular tests for nonadiabatic dynamics, reminiscent of the one-dimensional Tully models proposed in 1990.

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.

Initial sampling in symmetrical quasiclassical dynamics based on Li-Miller mapping Hamiltonian

A symmetrical quasiclassical (SQC) dynamics approach based on the Li-Miller (LM) mapping Hamiltonian (SQC-LM) was employed to describe nonadiabatic dynamics. In principle, the different initial sampling procedures may be applied in the SQC-LM dynamics, and the results may be dependent on different initial sampling. We provided various initial sampling approaches and checked their influence. We selected two groups of models including site-exciton models for exciton dynamics and linear vibronic coupling models for conical intersections to test the performance of SQC-LM dynamics with the different initial sampling methods. The results were examined with respect to those of the accurate multiconfigurational time-dependent Hartree (MCTDH) quantum dynamics. For both the models, the SQC-LM method more-or-less gives a reasonable description of the population dynamics, while the influence of the initial sampling approaches on the final results is noticeable. It seems that the suitable initial sampling methods should be determined by the system under study. This indicates that the combination of the SQC-LM method with a suitable sampling approach may be a potential method in the description of nonadiabatic dynamics.

Quantum Trajectory Mean-Field Method for Nonadiabatic Dynamics in Photochemistry

The mixed quantum-classical dynamical approaches have been widely used to study nonadiabatic phenomena in photochemistry and photobiology, in which the time evolutions of the electronic and nuclear subsystems are treated based on quantum and classical mechanics, respectively. The key issue is how to deal with coherence and decoherence during the propagation of the two subsystems, which has been the subject of numerous investigations for a few decades. A brief description on Ehrenfest mean-field and surface-hopping (SH) methods is first provided, and then different algorithms for treatment of quantum decoherence are reviewed in the present paper. More attentions were paid to quantum trajectory mean-field (QTMF) method under the picture of quantum measurements, which is able to overcome the overcoherence problem. Furthermore, the combined QTMF and SH algorithm is proposed in the present work, which takes advantages of the QTMF and SH methods. The potential to extend the applicability of the QTMF method was briefly discussed, such as the generalization to other type of nonadiabatic transitions, the combination with multiscale computational models, and possible improvements on its accuracy and efficiency by using machine-learning techniques.