Decoherence and Its Role in Electronically Nonadiabatic Dynamics

SPADE 1.0: A Simulation Package for Non-Adiabatic Dynamics in Extended Systems

Nonadiabatic molecular dynamics (NAMD) simulations are crucial for revealing the underlying mechanisms of photochemical and photophysical processes. Typical NAMD simulation software packages rely on on-the-fly ab initio electronic structure and nonadiabatic coupling calculations, and thus become challenging when dealing with large complex systems. We here introduce a new Simulation Package for non-Adiabatic Dynamics in Extended systems (SPADE), which is designed to address the limitations of traditional surface hopping methods in dealing with these problems. By design, SPADE enables the users to define arbitrary quasi-diabatic Hamiltonians through parametrized functions and incorporates a variety of algorithms (e.g., global flux hopping probabilities, complex crossing and decoherence corrections), which can realize efficient and reliable NAMD simulations without using nonadiabatic couplings at all. All the employed methods and expressions for diabatic Hamiltonian matrix elements can be flexibly set through the input files. SPADE is mainly written in Fortran based on a modular design and has a great capacity for further implementation of new methods. SPADE can be used to simulate both model and atomistic systems as long as proper Hamiltonians are provided. As demonstrations, a series of representative models are studied to show the main features and capabilities.

Generalization of Quantum-Trajectory Surface Hopping to Multiple Quantum States

In this work, we present a generalization of the quantum trajectory surface hopping (QTSH) to multiple states and its implementation in the Libra package for nonadiabatic dynamics. In lieu of the ad hoc velocity rescaling used in many trajectory-based surface hopping approaches, QTSH utilizes quantum forces to evolve nuclear degrees of freedom continuously. It also lifts the unphysical constraint of enforcing the total energy conservation at the individual trajectory level and rather conserves the total energy at the trajectory ensemble level. Leveraging our new implementation of the multistate QTSH, we perform a comparative analysis of this method with the conventional fewest switches surface hopping approach. We combine the QTSH and decoherence corrections based on the simplified decay of mixing (SDM) and exact factorization (XF), leading to the QTSH-SDM and QTSH-XF schemes. Using the Holstein, superexchange, and phenol model Hamiltonians, we assess the relative accuracy of the resulting combined schemes in reproducing branching ratios, population, and coherence dynamics for a broad range of initial conditions. We observe that the decoherence correction in QTSH is crucial to improve energy conservation as well as the internal consistency between the population from the quantum probability and active state.

Assessing the performance of coupled-trajectory schemes on full-dimensional two-state linear vibronic coupling models

We investigate the performance of coupled-trajectory methods for nonadiabatic molecular dynamics in simulating the photodynamics of 4-(dimethylamino)benzonitrile (DMABN) and fulvene, with electronic structure provided by linear vibrational coupling models. We focus on the coupled-trajectory mixed quantum-classical (CTMQC) algorithm and on the (combined) coupled-trajectory Tully surface hopping [(C)CTTSH] in comparison to independent-trajectory approaches, such as multi-trajectory Ehrenfest and Tully surface hopping. Our analysis includes not only electronic populations but also additional electronic and nuclear properties in position and momentum space. For both DMABN and fulvene, the recently developed CCTTSH algorithm successfully resolves the internal inconsistencies of coupled-trajectory Tully surface hopping. Instead, we find that DMABN highlights a significant weakness of CTMQC, which arises when the trajectories remain for a long time in the vicinity of a region of strong nonadiabaticity.

Decoherence dynamics in molecular qubits: Exponential, Gaussian and beyond

In this work, we examine how the structure of system-bath interactions can determine commonly encountered temporal decoherence patterns, such as Gaussian and exponential decay, in molecular and other qubits coupled to a thermal bosonic bath. The analysis, based on a pure dephasing picture that admits analytical treatment, shows that decoherence, in general, is neither purely Gaussian nor exponential but rather the exponential of oscillatory functions, with periods determined by the bath's frequencies. For initially unentangled qubit-bath states, Gaussian decay is always present at early times. It becomes increasingly dominant with increasing temperature, qubit-bath interaction strength, and bath correlation time. Initial system-bath entanglement that arises due to displacement in the position of the bath states preserves the Gaussian decay. By contrast, strict exponential decay arises only in very specific models that we isolate. However, it becomes dominant for times longer than the bath correlation time or for early times when there is initial entanglement due to momentum displacement of the bath states. For molecular electronic decoherence, the long-time exponential regime plays a limited role as it emerges after most coherence is lost. Thus, the Gaussian decay provides a more suitable (albeit imperfect) model of such decoherence. Furthermore, we discuss the connection between electronic decoherence dynamics and electronic spectroscopic line shape theory, where Gaussian spectral peaks correspond to Gaussian coherence decay and Lorentzian peaks correspond to exponential coherence decay. We find that Gaussian spectral peaks, usually associated with inhomogeneous broadening, can emerge from the entangling unitary system-bath dynamics even when there is no inhomogeneity in the initial conditions.

Electronically Nonadiabatic Quenching of Excited States of O2 by Collisions with O Atoms

The kinetics of electronically inelastic quenching of O2(a1Δg) and O2(b1Σg+) by collisions with O(3P) have been investigated using mixed quantum-classical trajectories governed by adiabatic potential energy surfaces and state couplings generated from a recently developed diabatic potential energy matrix (DPEM) for the 14 lowest-energy 3A' states of O3. Using the coherent switching with decay of mixing (CSDM) method, dynamics calculations were performed both with 14 coupled electronic states and with 8 coupled electronical states, and similar results were obtained. The calculated thermal quenching rate coefficients are generally small, but they increase with temperature. The positive temperature dependence is attributed to high-energy locally avoided crossings that can only be easily accessed by high collision energies. We find that, depending on the temperature, 86-97% of the b state quenches are into the a state with the remainder into the ground electronic state of O2. The calculated rate coefficients for quenching of O2(b1Σg+) and O2(a1Δg) by O(3P), coupled with the assumption that electronically nonadiabatic probabilities for collisions on the 3A' surfaces are similar to those on 3A″ surfaces, are compared with the available experimental results.

Spin-Orbit Coupling and Admixture Coefficients in SA-CASSCF and MS-CASPT2, and Triplet Excitation Yield Simulated via Trajectory Surface Hopping and Calibrated SA-CASSCF in 1,2-Dioxetane Derivatives

The energy gaps, spin-orbit coupling (SOC), and admixture coefficients over a series of the configurations are evaluated by the SA-CASSCF/6-31G, SA-CASSCF/6-31G*, SA-CASSCF/ANO-RCC-VDZP, and MS-CASPT2/ANO-RCC-VDZP to reveal the extent of the inaccuracy of the SA-CASSCF. By comparing the mean absolute errors for the energy gaps and the admixture coefficient magnitudes (ACMs) measured between the SA-CASSCF/6-31G, SA-CASSCF/6-31G*, or SA-CASSCF/ANO-RCC-VDZP and the MS-CASPT2/ANO-RCC-VDZP, the SA-CASSCF/6-31G is selected as the electronic structure method in the nonadiabatic molecular dynamics simulation. The major components of the ACMs of the SA-CASSCF/6-31G and MS-CASPT2/ANO-RCC-VDZP are identified and compared; we find that the ACMs are underestimated by the SA-CASSCF/6-31G, which is verified by the reasonable triplet quantum yield simulated by the trajectory surface hopping and the calibrated SA-CASSCF/6-31G. The magnitude of the singlet-triplet mixing positively correlates to the hopping probability between the mixed singlet and triplet states, which is confirmed by the computed S-T transition probability.

Coherence in Chemistry: Foundations and Frontiers

Coherence refers to correlations in waves. Because matter has a wave-particle nature, it is unsurprising that coherence has deep connections with the most contemporary issues in chemistry research (e.g., energy harvesting, femtosecond spectroscopy, molecular qubits and more). But what does the word "coherence" really mean in the context of molecules and other quantum systems? We provide a review of key concepts, definitions, and methodologies, surrounding coherence phenomena in chemistry, and we describe how the terms "coherence" and "quantum coherence" refer to many different phenomena in chemistry. Moreover, we show how these notions are related to the concept of an interference pattern. Coherence phenomena are indeed complex, and ambiguous definitions may spawn confusion. By describing the many definitions and contexts for coherence in the molecular sciences, we aim to enhance understanding and communication in this broad and active area of chemistry.

Simulations of photoinduced processes with the exact factorization: state of the art and perspectives

This perspective offers an overview of the applications of the exact factorization of the electron-nuclear wavefunction to the domain of theoretical photochemistry, where the aim is to gain insights into the ultrafast dynamics of molecular systems via simulations of their excited-state dynamics beyond the Born-Oppenheimer approximation. The exact factorization offers an alternative viewpoint to the Born-Huang representation for the interpretation of dynamical processes involving the electronic ground and excited states as well as their coupling through the nuclear motion. Therefore, the formalism has been used to derive algorithms for quantum molecular-dynamics simulations where the nuclear motion is treated using trajectories and the electrons are treated quantum mechanically. These algorithms have the characteristic features of being based on coupled and on auxiliary trajectories, and have shown excellent performance in describing a variety of excited-state processes, as this perspective illustrates. We conclude with a discussion on the authors' point of view on the future of the exact factorization.

Extending Non-Perturbative Simulation Techniques for Open-Quantum Systems to Excited-State Proton Transfer and Ultrafast Non-Adiabatic Dynamics

Excited state proton transfer is an ubiquitous phenomenon in biology and chemistry, spanning from the ultrafast reactions of photobases and acids to light-driven, enzymatic catalysis and photosynthesis. However, the simulation of such dynamics involves multiple challenges, since high-dimensional, out-of-equilibrium vibronic states play a crucial role, while a fully quantum description of the proton's dissipative, real-space dynamics is also required. In this work, we extend the powerful matrix product state approach to open quantum systems (TEDOPA) to study these demanding dynamics, and also more general nonadiabatic processes that can appear in complex photochemistry subject to strong laser driving. As an illustration, we initially consider an open model of a four-level electronic system interacting with hundreds of intramolecular vibrations that drive ultrafast excited state proton transfer, as well as an explicit photonic environment that allows us to directly monitor the resulting dual fluorescence in this system. We then demonstrate how to include a continuous "reaction coordinate" of the proton transfer that allows numerically exact simulations that can be understood, visualized and interpreted in the familiar language of diabatic and adiabatic dynamics on potential surfaces, while also retaining an exact quantum treatment of dissipation and driving effects that could be used to study diverse problems in ultrafast photochemistry.

Machine Learning-Assisted Mixed Quantum-Classical Dynamics without Explicit Nonadiabatic Coupling: Application to the Photodissociation of Peroxynitric Acid

We have devised a hybrid quantum-classical scheme utilizing machine-learned potential energy surfaces (PES), which circumvents the need for explicit computation of nonadiabatic coupling elements. The quantities necessary to account for the nonadiabatic effects are directly obtained from the PESs. The simulation of dynamics is based on the fewest-switches surface-hopping method. We applied this scheme to model the photodissociation of both N-O and O-O bonds in a conformer of peroxynitric acid (HO2NO2). Adiabatic PES data for the six lowest states of this molecule were computed at the CASSCF level for various nuclear configurations. These served as the training data for the machine-learning models for the PESs. The dynamics simulation was initiated on the lowest optically bright singlet excited state (S4) and propagated along the two Jacobi coordinates J → 1 and J → 2 while accounting for the nonadiabatic effects through transitions between PESs. Our analysis revealed that there is a very high chance of dissociation of the N-O bond leading to the HO2 and NO2 fragments.

Time-Dependent Extension of Grimme's Continuous Chirality Measure for Electronic Chirality Flips in Femto- and Attosecond Time Domains

Grimme's Continuous Chirality Measure (C C M ${CCM}$ ) was developed for comparisons of the chirality of the electronic wave functions of molecules, typically in their ground states. For example,C C M = 14 . 5 ${CCM=14.5}$ ,1 . 2 ${1.2}$ and0 . 0 ${0.0}$ for alanine, hydrogen-peroxide, and for achiral molecules, respectively. Well-designed laser pulses can excite achiral molecules from the electronic ground state to time-dependent chiral superposition states, with chirality flips in the femto- or even attosecond (fs or as) time domains. Here we provide a time-dependent extensionC C M t ${CCM\left(t\right)}$ of Grimme'sC C M ${CCM}$ for trailing the electronic chirality flips. As examples, we consider two laser driven electronic wavefunctions which represent flips between opposite electronic enantiomers of oriented NaK within4 . 76 f s ${4.76\ {\rm f}{\rm s}}$ and433 a s ${433\ {\rm a}{\rm s}}$ . The correspondingC C M t ${CCM\left(t\right)}$ vary respectively from14 . 5 ${14.5}$ or from13 . 3 ${13.3}$ to0 . 0 ${0.0}$ , and back.

Quantum Dynamics from Classical Trajectories

Nonadiabatic molecular dynamics plays an essential role in exploring the time evolution of molecular systems. Various methods have been developed for this study, with varying accuracy and computational cost. One very successful among them is trajectory surface hopping, which propagates nuclei as classical trajectories using forces from a quantum description of the electrons and incorporates nonadiabatic effects through stochastic state changes during each trajectory propagation. A statistical analysis of an ensemble of the independent trajectories recovers the simulated system's behavior. This approach can give good results, but it is known to overlook nuclear quantum effects, leading to inaccurate predictions. Here, we present quantum dynamics from classical trajectories (QDCT), a new protocol to recover the quantum wavepacket from the classical trajectories generated by surface hopping. In this first QDCT implementation, we apply it to recover results at the multiple spawning level from postprocessing surface hopping precomputed trajectories. With a series of examples, we demonstrate QDCT's potential to improve the accuracy of the dynamics, correct decoherence effects, and diagnose problems or increase confidence in surface hopping results. All that comes at virtually no computational cost since no new electronic calculation is required.

Gravitationally-induced wave function collapse time for molecules

The Diósi-Penrose model states that the wave function collapse ending a quantum superposition occurs due to the instability of coexisting gravitational potentials created by distinct geometric conformations of the system in different states. The Heisenberg time-energy principle can be invoked to estimate the collapse time for the energy associated with this instability, the gravitational self-energy. This paper develops atomistic models to calculate the Diósi-Penrose collapse time. It applies them to a range of systems, from small molecules to large biological structures and macroscopic systems. An experiment is suggested to test the Diósi-Penrose hypothesis, and we critically examine the model, highlighting challenges from an atomistic perspective, such as gravitational self-energy saturation and limited extensivity.

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.

Generalized Semiclassical Ehrenfest Method: A Route to Wave Function-Free Photochemistry and Nonadiabatic Dynamics with Only Potential Energies and Gradients

We reconsider recent methods by which direct dynamics calculations of electronically nonadiabatic processes can be carried out while requiring only adiabatic potential energies and their gradients. We show that these methods can be understood in terms of a new generalization of the well-known semiclassical Ehrenfest method. This is convenient because it eliminates the need to evaluate electronic wave functions and their matrix elements along the mixed quantum-classical trajectories. The new approximations and procedures enabling this advance are the curvature-driven approximation to the time-derivative coupling, the generalized semiclassical Ehrenfest method, and a new gradient correction scheme called the time-derivative matrix (TDM) scheme. When spin-orbit coupling is present, one can carry out dynamics calculations in the fully adiabatic basis using potential energies and gradients calculated without spin-orbit coupling plus the spin-orbit coupling matrix elements. Even when spin-orbit coupling is neglected, the method is useful because it allows calculations by electronic structure methods for which nonadiabatic coupling vectors are unavailable. In order to place the new considerations in context, the article starts out with a review of background material on trajectory surface hopping, the semiclassical Ehrenfest scheme, and methods for incorporating decoherence. We consider both internal conversion and intersystem crossing. We also review several examples from our group of successful applications of the curvature-driven approximation.

Prediction Challenge: Simulating Rydberg photoexcited cyclobutanone with surface hopping dynamics based on different electronic structure methods

This research examines the nonadiabatic dynamics of cyclobutanone after excitation into the n → 3s Rydberg S2 state. It stems from our contribution to the Special Topic of the Journal of Chemical Physics to test the predictive capability of computational chemistry against unseen experimental data. Decoherence-corrected fewest-switches surface hopping was used to simulate nonadiabatic dynamics with full and approximated nonadiabatic couplings. Several simulation sets were computed with different electronic structure methods, including a multiconfigurational wavefunction [multiconfigurational self-consistent field (MCSCF)] specially built to describe dissociative channels, multireference semiempirical approach, time-dependent density functional theory, algebraic diagrammatic construction, and coupled cluster. MCSCF dynamics predicts a slow deactivation of the S2 state (10 ps), followed by an ultrafast population transfer from S1 to S0 (<100 fs). CO elimination (C3 channel) dominates over C2H4 formation (C2 channel). These findings radically differ from the other methods, which predicted S2 lifetimes 10-250 times shorter and C2 channel predominance. These results suggest that routine electronic structure methods may hold low predictive power for the outcome of nonadiabatic dynamics.

Dissociation Time, Quantum Yield, and Dynamic Reaction Pathways in the Thermolysis of trans-3,4-Dimethyl-1,2-dioxetane

The thermolysis of trans-3,4-dimethyl-1,2-dioxetane is studied by trajectory surface hopping. The significant difference between long and short dissociation times is rationalized by frustrated dissociations and the time spent in triplet states. If the C-C bond breaks through an excited state channel, then the trajectory passes over a ridge of the potential energy surface of that state. The calculated triplet quantum yields match the experimental results. The dissociation half-times and quantum yields follow the same ascending order as per the product states, justifying the conjecture that the longer dissociation time leads to a higher quantum yield, proposed in the context of the methylation effect. The populations of the molecular Coulomb Hamiltonian and diagonal states reach equilibrium, but the triplet populations with different Sz components fluctuate indefinitely. Certain initial velocities, leading the trajectories to given product states, can be identified as the most characteristic features for sorting trajectories according to their product states.

Total angular momentum conservation in Ehrenfest dynamics with a truncated basis of adiabatic states

We show that standard Ehrenfest dynamics does not conserve linear and angular momentum when using a basis of truncated adiabatic states. However, we also show that previously proposed effective Ehrenfest equations of motion [M. Amano and K. Takatsuka, "Quantum fluctuation of electronic wave-packet dynamics coupled with classical nuclear motions," J. Chem. Phys. 122, 084113 (2005) and V. Krishna, "Path integral formulation for quantum nonadiabatic dynamics and the mixed quantum classical limit," J. Chem. Phys. 126, 134107 (2007)] involving the non-Abelian Berry force do maintain momentum conservation. As a numerical example, we investigate the Kramers doublet of the methoxy radical using generalized Hartree-Fock with spin-orbit coupling and confirm that angular momentum is conserved with the proper equations of motion. Our work makes clear some of the limitations of the Born-Oppenheimer approximation when using ab initio electronic structure theory to treat systems with unpaired electronic spin degrees of freedom, and we demonstrate that Ehrenfest dynamics can offer much improved, qualitatively correct results.

Investigating the Photodynamics of trans-Azobenzene with Coupled Trajectories

In this work, we present the first implementation of coupled-trajectory Tully surface hopping (CT-TSH) suitable for applications to molecular systems. We combine CT-TSH with the semiempirical floating occupation molecular orbital-configuration interaction electronic structure method to investigate the photoisomerization dynamics of trans-azobenzene. Our study shows that CT-TSH can capture correctly decoherence effects in this system, yielding consistent electronic and nuclear dynamics in agreement with (standard) decoherence-corrected TSH. Specifically, CT-TSH is derived from the exact factorization and the electronic coefficients' evolution is directly influenced by the coupling of trajectories, resulting in the improvement of internal consistency if compared to standard TSH.

From chiral laser pulses to femto- and attosecond electronic chirality flips in achiral molecules

Chirality is an important topic in biology, chemistry and physics. Here we show that ultrashort circularly polarized laser pulses, which are chiral, can be fired on achiral oriented molecules to induce chirality in their electronic densities, with chirality flips within femtoseconds or even attoseconds. Our results, obtained by quantum dynamics simulations, use the fact that laser pulses can break electronic symmetry while conserving nuclear symmetry. Here two laser pulses generate a superposition of three electronic eigenstates. This breaks all symmetry elements of the electronic density, making it chiral except at the periodic rare events of the chirality flips. As possible applications, we propose the combination of the electronic chirality flips with Chiral Induced Spin Selectivity.

Mapping electronic decoherence pathways in molecules

Establishing the fundamental chemical principles that govern molecular electronic quantum decoherence has remained an outstanding challenge. Fundamental questions such as how solvent and intramolecular vibrations or chemical functionalization contribute to the decoherence remain unanswered and are beyond the reach of state-of-the-art theoretical and experimental approaches. Here we address this challenge by developing a strategy to isolate electronic decoherence pathways for molecular chromophores immersed in condensed phase environments that enables elucidating how electronic quantum coherence is lost. For this, we first identify resonance Raman spectroscopy as a general experimental method to reconstruct molecular spectral densities with full chemical complexity at room temperature, in solvent, and for fluorescent and non-fluorescent molecules. We then show how to quantitatively capture the decoherence dynamics from the spectral density and identify decoherence pathways by decomposing the overall coherence loss into contributions due to individual molecular vibrations and solvent modes. We illustrate the utility of the strategy by analyzing the electronic decoherence pathways of the DNA base thymine in water. Its electronic coherences decay in [Formula: see text]30 fs. The early-time decoherence is determined by intramolecular vibrations while the overall decay by solvent. Chemical substitution of thymine modulates the decoherence with hydrogen-bond interactions of the thymine ring with water leading to the fastest decoherence. Increasing temperature leads to faster decoherence as it enhances the importance of solvent contributions but leaves the early-time decoherence dynamics intact. The developed strategy opens key opportunities to establish the connection between molecular structure and quantum decoherence as needed to develop chemical strategies to rationally modulate it.

What Controls the Quality of Photodynamical Simulations? Electronic Structure Versus Nonadiabatic Algorithm

The field of nonadiabatic dynamics has matured over the last decade with a range of algorithms and electronic structure methods available at the moment. While the community currently focuses more on developing and benchmarking new nonadiabatic dynamics algorithms, the underlying electronic structure controls the outcome of nonadiabatic simulations. Yet, the electronic-structure sensitivity analysis is typically neglected. In this work, we present a sensitivity analysis of the nonadiabatic dynamics of cyclopropanone to electronic structure methods and nonadiabatic dynamics algorithms. In particular, we compare wave function-based CASSCF, FOMO-CASCI, MS- and XMS-CASPT2, density-functional REKS, and semiempirical MRCI-OM3 electronic structure methods with the Landau-Zener surface hopping, fewest switches surface hopping, and ab initio multiple spawning with informed stochastic selection algorithms. The results clearly demonstrate that the electronic structure choice significantly influences the accuracy of nonadiabatic dynamics for cyclopropanone even when the potential energy surfaces exhibit qualitative and quantitative similarities. Thus, selecting the electronic structure solely on the basis of the mapping of potential energy surfaces can be misleading. Conversely, we observe no discernible differences in the performance of the nonadiabatic dynamics algorithms across the various methods. Based on the above results, we discuss the present-day practice in computational photodynamics.

Nonadiabatic Dynamics Simulations for Photoinduced Processes in Molecules and Semiconductors: Methodologies and Applications

Nonadiabatic dynamics (NAMD) simulations have become powerful tools for elucidating complicated photoinduced processes in various systems from molecules to semiconductor materials. In this review, we present an overview of our recent research on photophysics of molecular systems and periodic semiconductor materials with the aid of ab initio NAMD simulation methods implemented in the generalized trajectory surface-hopping (GTSH) package. Both theoretical backgrounds and applications of the developed NAMD methods are presented in detail. For molecular systems, the linear-response time-dependent density functional theory (LR-TDDFT) method is primarily used to model electronic structures in NAMD simulations owing to its balanced efficiency and accuracy. Moreover, the efficient algorithms for calculating nonadiabatic coupling terms (NACTs) and spin-orbit couplings (SOCs) have been coded into the package to increase the simulation efficiency. In combination with various analysis techniques, we can explore the mechanistic details of the photoinduced dynamics of a range of molecular systems, including charge separation and energy transfer processes in organic donor-acceptor structures, ultrafast intersystem crossing (ISC) processes in transition metal complexes (TMCs), and exciton dynamics in molecular aggregates. For semiconductor materials, we developed the NAMD methods for simulating the photoinduced carrier dynamics within the framework of the Kohn-Sham density functional theory (KS-DFT), in which SOC effects are explicitly accounted for using the two-component, noncollinear DFT method. Using this method, we have investigated the photoinduced carrier dynamics at the interface of a variety of van der Waals (vdW) heterojunctions, such as two-dimensional transition metal dichalcogenides (TMDs), carbon nanotubes (CNTs), and perovskites-related systems. Recently, we extended the LR-TDDFT-based NAMD method for semiconductor materials, allowing us to study the excitonic effects in the photoinduced energy transfer process. These results demonstrate that the NAMD simulations are powerful tools for exploring the photodynamics of molecular systems and semiconductor materials. In future studies, the NAMD simulation methods can be employed to elucidate experimental phenomena and reveal microscopic details as well as rationally design novel photofunctional materials with desired properties.

Nonadiabatic Coupling in Trajectory Surface Hopping: Accurate Time Derivative Couplings by the Curvature-Driven Approximation

Trajectory surface hopping (TSH) is a widely used mixed quantum-classical dynamics method that is used to simulate molecular dynamics with multiple electronic states. In TSH, time-derivative coupling is employed to propagate the electronic coefficients and in that way to determine when the electronic state on which the nuclear trajectory is propagated switches. In this work, we discuss nonadiabatic TSH dynamics algorithms employing the curvature-driven approximation and overlap-based time derivative couplings, and we report test calculations on six photochemical reactions where we compare the results to one another and to calculations employing analytic nonadiabatic coupling vectors. We correct previous published results thanks to a bug found in the software. We also provide additional, more detailed studies of the time-derivative couplings. Our results show good agreement between curvature-driven algorithms and overlap-based algorithms.

Charge Migration in HCCI Cations Probed by Strong Field Ionization: Time-Dependent Configuration Interaction and Vibrational Wavepacket Simulations

Strong field ionization of neutral iodoacetylene (HCCI) can produce a coherent superposition of the X and A cations and results in charge migration between the CC π orbital and the iodine π-type lone pair. This charge migration causes oscillations in the rate of strong field ionization of the cation to the dication that can be monitored using intense few-cycle probe pulses. The dynamics and strong field ionization of the coherent superposition the X and A states of HCCI⁺ have been modeled by time-dependent configuration interaction (TDCI) simulations. When the nuclei are allowed to move, the electronic wavefunctions need to be multiplied by vibrational wavefunctions. Nuclear motion has been modeled by vibrational packets moving on quadratic approximations to the potential energy surfaces for the X and A states of the cation. The overlap of the vibrational wavepackets decays in about 10-15 fs. Consequently, the oscillations in the strong field ionization decay on the same time scale. A revival of the vibrational overlap and in the oscillations of the strong field ionization is seen at 60-110 fs. TDCI simulations show that the decay and revival of the charge migration can be monitored by strong field ionization with intense 2- and 4-cycle linearly polarized 800 nm pulses. The revival is also seen with 7-cycle pulses.

New Gradient Correction Scheme for Electronically Nonadiabatic Dynamics Involving Multiple Spin States

It has been recommended that the best representation to use for trajectory surface hopping (TSH) calculations is the fully adiabatic basis in which the Hamiltonian is diagonal. Simulations of intersystem crossing processes with conventional TSH methods require an explicit computation of nonadiabatic coupling vectors (NACs) in the molecular-Coulomb-Hamiltonian (MCH) basis, also called the spin-orbit-free basis, in order to compute the gradient in the fully adiabatic basis (also called the diagonal representation). This explicit requirement destroys some of the advantages of the overlap-based algorithms and curvature-driven algorithms that can be used for the most efficient TSH calculations. Therefore, although these algorithms allow one to perform NAC-free simulations for internal conversion processes, one still requires NACs for intersystem crossing. Here, we show that how the NAC requirement is circumvented by a new computation scheme called the time-derivative-matrix scheme.

Direct Nonadiabatic Dynamics of Ammonia with Curvature-Driven Coherent Switching with Decay of Mixing and with Fewest Switches with Time Uncertainty: An Illustration of Population Leaking in Trajectory Surface Hopping Due to Frustrated Hops

Mixed quantum-classical nonadiabatic dynamics is a widely used approach to simulate molecular dynamics involving multiple electronic states. There are two main categories of mixed quantum-classical nonadiabatic dynamics algorithms, namely, trajectory surface hopping (TSH) in which the trajectory propagates on a single potential energy surface, interrupted by hops, and self-consistent-potential (SCP) methods, such as semiclassical Ehrenfest, in which propagation occurs on a mean-field surface without hops. In this work, we will illustrate an example of severe population leaking in TSH. We emphasize that such leaking is a combined effect of frustrated hops and long-time simulations that drive the final excited-state population toward zero as a function of time. We further show that such leaking can be alleviated-but not eliminated-by the fewest switches with time uncertainty TSH algorithm (here implemented in the SHARC program); the time uncertainty algorithm slows down the leaking process by a factor of 4.1. The population leaking is not present in coherent switching with decay of mixing (CSDM), which is an SCP method with non-Markovian decoherence included. Another result in this paper is that we find very similar results with the original CSDM algorithm, with time-derivative CSDM (tCSDM), and with curvature-driven CSDM (κCSDM). Not only do we find good agreement for electronically nonadiabatic transition probabilities but also we find good agreement of the norms of the effective nonadiabatic couplings (NACs) that are derived from the curvature-driven time-derivative couplings as implemented in κCSDM with the time-dependent norms of the nonadiabatic coupling vectors computed by state-averaged complete-active-space self-consistent field theory.