Assessing Nonadiabatic Dynamics Methods in Long Timescales.

In this work, the recently introduced quantum trajectory mean-field (QTMF) approach is implemented and employed to explore photodissociation dynamics of diazirinone (N2CO), which are based on the high-level ab initio calculation. For comparison, the photodissociation process has been simulated as well with the fewest-switches surface hopping (FSSH) and the ab initio multiple spawning (AIMS) methods. Overall, the dynamical behavior predicted by the three methods is consistent. The N2CO photodissociation at λ > 335 nm is an ultrafast process and the two C-N bonds are broken in a stepwise way, giving birth to CO and N2 as the final products in the ground state. Meanwhile, some noticeable differences were found in the QTMF, FSSH, and AIMS simulated time constants for fission of the C-N bonds, excited-state lifetime, and nonadiabatic transition ratios in different intersection regions. These have been discussed in detail. The present study provides a clear evidence that direct ab initio QTMF approach is one of the reliable tools for simulating nonadiabatic dynamics processes.

Accurate quantum dynamics simulations of nonadiabatic processes are important for studies of electron transfer, energy transfer, and photochemical reactions in complex systems. In this comparative study, we benchmark various approximate nonadiabatic dynamics methods with mapping variables against numerically exact calculations based on the tensor-train (TT) representation of high-dimensional arrays, including TT-KSL for zero-temperature dynamics and TT-thermofield dynamics for finite-temperature dynamics. The approximate nonadiabatic dynamics methods investigated include mixed quantum-classical Ehrenfest mean-field and fewest-switches surface hopping, linearized semiclassical mapping dynamics, symmetrized quasiclassical dynamics, the spin-mapping method, and extended classical mapping models. Different model systems were evaluated, including the spin-boson model for nonadiabatic dynamics in the condensed phase, the linear vibronic coupling model for electronic transition through conical intersections, the photoisomerization model of retinal, and Tully's one-dimensional scattering models. Our calculations show that the optimal choice of approximate dynamical method is system-specific, and the accuracy is sensitively dependent on the zero-point-energy parameter and the initial sampling strategy for the mapping variables.

We perform on-the-fly non-adiabatic molecular dynamics simulations using the recently developed spin-mapping formalism. Two quantum dynamics approaches based on this mapping formalism, (i) the fully linearized Spin-LSC and (ii) the partially linearized Spin-PLDM, are explored using the quasi-diabatic propagation scheme. We have performed dynamics simulations in four abinitio molecular models for which benchmark abinitio multiple spawning (AIMS) data have been published. We find that the spin-LSC and the previously reported symmetric quasi-classical (SQC) approaches provide nearly equivalent population dynamics. While we expected the more involved spin-PLDM method to provide superior accuracy compared to the other mapping-based approaches, SQC and spin-LSC, we found that it performed with equivalent accuracy compared to the AIMS benchmark results. We further explore the underpinnings of the spin-PLDM correlation function by decomposing its N2 density matrix-focused initial conditions, where N is the number of states in the quantum subsystem. Finally, we found an approximate form of the spin-PLDM correlation function, which simplifies the simulation and reduces the computational costs from N2 to N.

The recent improvement in experimental capabilities for interrogating and controlling molecular systems with ultrafast coherent light sources calls for the development of theoretical approaches that can accurately and efficiently treat electronic coherence. However, the most popular and practical nonadiabatic molecular dynamics techniques, Tully's fewest-switches surface hopping and Ehrenfest mean-field dynamics, are unable to describe the dynamics proceeding from an initial electronic coherence. While such issues are not encountered with the analogous coupled-trajectory algorithms or numerically exact quantum dynamics methods, applying such techniques necessarily comes with a higher computational cost. Here we show that a correct description of initial electronic coherence can indeed be achieved using independent-trajectory methods derived from the semiclassical mapping formalism. The key is the introduction of an initial sampling over the electronic phase space and a means of incorporating phase interference between trajectories, both of which are naturally achieved when working within the semiclassical mapping framework.

We introduce an efficient configuration interaction (CI) method for the calculation of mixed quantum and classical nonadiabatic molecular dynamics for multiple electrons. For any given realization of the classical degrees of freedom (e.g., a solvent), the method uses a novel real-space quadrature to efficiently compute the Coulomb and exchange interactions between electrons. We also introduce an approximation whereby the classical molecular dynamics is propagated for several time steps on electronic potential energy surfaces generated using only a particularly important subset of the CI basis states. By only updating the important-states subset periodically, we achieve significant reductions in the computational cost of solving the multielectron quantum problem. We test the real-space quadrature for the cases of two electrons confined in a cubic box with infinitely repulsive walls and two electrons dissolved in liquid water that occupy a single cavity, so-called hydrated dielectrons. We then demonstrate how to perform mixed quantum and classical nonadiabatic dynamics by combining these computational techniques with the mean-field with surface hopping algorithm of Prezhdo and Rossky [J. Chem. Phys. 107, 825 (1997)]. Finally, we illustrate the practicality of the approach to multielectron nonadiabatic dynamics by examining the nonadiabatic relaxation dynamics of both spin singlet and spin triplet hydrated dielectrons following excitation from the ground to the first excited state.

A mechanism explaining the photobasicity of 5-methoxyquinoline (5-MeOQ) is proposed on the basis of nonadiabatic molecular dynamics simulations using time-dependent hybrid density functional theory (TDDFT) and fewest switches surface hopping (SH) and analysis of existing ultrafast spectroscopy experiments. According to the TDDFT-SH simulations, the rate-determining step is hole transfer from photoexcited 5-MeOQ to adjacent water molecules within ∼5 ps followed by rapid electron-coupled proton transfer and deactivation to the ground state. This fast redox-catalyzed proton transfer mechanism is consistent with simple thermodynamic arguments, correlated wavefunction calculations, and recent isotope substitution time-resolved fluorescence experiments. Although energetically feasible, direct protonation of 5-MeOQ in the S1 state appears to be too slow to contribute significantly to 5-MeOQ photobasicity in aqueous solution.

A new scheme is proposed for modeling molecular nonadiabatic dynamics near metal surfaces. The charge-transfer character of such dynamics is exploited to construct an efficient reduced representation for the electronic structure. In this representation, the fewest switches surface hopping (FSSH) approach can be naturally modified to include electronic relaxation (ER). The resulting FSSH-ER method is valid across a wide range of coupling strengths as supported by tests applied to the Anderson-Holstein model for electron transfer. Future work will combine this scheme with ab initio electronic structure calculations.

Photochemical and photophysical properties of cis-stilbene molecule by electronic structure calculations and nonadiabatic surface-hopping dynamics 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.

The photochemistry of cis-stilbene proceeds through two pathways: cis-trans isomerization and ring closure to 4a,4b-dihydrophenanthrene (DHP). Despite serving for many decades as a model system for photoisomerization, the photodynamics of cis-stilbene is still not fully understood. We use ab initio multiple spawning on a SA-2-CASSCF(2,2) potential energy surface to simulate the nonadiabatic dynamics of isolated cis-stilbene. We find the cyclization (to DHP and cis-stilbene) and isomerization (to trans- and cis-stilbene) reaction coordinates to be orthogonal; branching between the two pathways is determined on the S1 excited state within 150 fs of photoexcitation. Trajectory basis functions (TBFs) undergoing cyclization decay rapidly to the ground state in 250 fs, while TBFs moving along the isomerization coordinate remain on the excited state longer, with the majority decaying between 300 and 500 fs. We observe three avoided crossing regions in the dynamics: two along the isomerization coordinate (displaying pyramidalization and migration of an ethylenic hydrogen or phenyl group), and one DHP-like conical intersection along the cyclization coordinate. The isomeric form of the vibrationally hot photoproducts (as determined by measurement 2 ps after photoexcitation) is determined within less than 50 fs of decay to the ground state mediated by passage through a conical intersection. Excess vibrational energy of ground state cis- and trans-stilbene is channelled into phenyl torsions (with mostly opposing directionality). Our simulations are validated by direct comparison to experiment for the absorption spectrum, branching ratio of the three photoproducts (44:52:4 cis-stilbene:trans-stilbene:DHP), and excited state lifetime (520 ± 40 fs).

Full multiple spawning offers an in principle exact framework for excited-state dynamics, where nuclear wavefunctions in different electronic states are represented by a set of coupled trajectory basis functions that follow classical trajectories. The couplings between trajectory basis functions can be approximated to treat molecular systems, leading to the ab initio multiple spawning method which has been successfully employed to study the photochemistry and photophysics of several molecules. However, a detailed investigation of its approximations and their consequences is currently missing in the literature. In this work, we simulate the explicit photoexcitation and subsequent excited-state dynamics of a simple system, LiH, and we analyze (i) the effect of the ab initio multiple spawning approximations on different observables and (ii) the convergence of the ab initio multiple spawning results towards numerically exact quantum dynamics upon a progressive relaxation of these approximations. We show that, despite the crude character of the approximations underlying ab initio multiple spawning for this low-dimensional system, the qualitative excited-state dynamics is adequately captured, and affordable corrections can further be applied to ameliorate the coupling between trajectory basis functions.

Simulations of photochemical reaction dynamics have been a challenge to the theoretical chemistry community for some time. In an effort to determine the predictive character of current approaches, we predict the results of an upcoming ultrafast diffraction experiment on the photodynamics of cyclobutanone after excitation to the lowest lying Rydberg state (S2). A picosecond of nonadiabatic dynamics is described with abinitio multiple spawning. We use both time dependent density functional theory (TDDFT) and equation-of-motion coupled cluster singles and doubles (EOM-CCSD) theory for the underlying electronic structure theory. We find that the lifetime of the S2 state is more than a picosecond (with both TDDFT and EOM-CCSD). The predicted ultrafast electron diffraction spectrum exhibits numerous structural features, but weak time dependence over the course of the simulations.

Upon photoexcitation by a short light pulse, molecules can reach regions of the configuration space characterized by strong nonadiabaticity, where the motion of the nuclei is strongly coupled to the motion of the electrons. The subtle interplay between the nuclear and electronic degrees of freedom in such situations is rather challenging to capture by state-of-the-art nonadiabatic dynamics approaches, limiting therefore their predictive power. The Exact Factorization of the molecular wavefunction, though, offers new perspectives in the solution of this longstanding issue. Here, we investigate the performance of a mixed quantum/classical (MQC) limit of this theory, named Coupled Trajectory-MQC, which was shown to reproduce the excited-state dynamics of small systems accurately. The method is applied to the study of the photoinduced ring opening of oxirane and the results are compared with two other nonadiabatic approaches based on different Ansatze for the molecular wavefunction, namely Ehrenfest dynamics and Ab Initio Multiple Spawning (AIMS). All simulations were performed using linear-response time-dependent density functional theory. We show that the CT-MQC method can capture the (de)coherence effects resulting from the dynamics through conical intersections, in good agreement with the results obtained with AIMS and in contrast with ensemble Ehrenfest dynamics.

Fewest-switches surface hopping (FSSH) is a popular nonadiabatic dynamics method which treats nuclei with classical mechanics and electrons with quantum mechanics. In order to simulate the motion of a wave packet as accurately as possible, standard FSSH requires a stochastic sampling of the trajectories over a distribution of initial conditions corresponding, e.g., to the Wigner distribution of the initial quantum wave packet. Although it is well-known that FSSH does not properly account for decoherence effects, there is some confusion in the literature about whether or not this averaging over a distribution of initial conditions can approximate some of the effects of decoherence. In this paper, we not only show that averaging over initial conditions does not generally account for decoherence, but also why it fails to do so. We also show how an apparent improvement in accuracy can be obtained for a fortuitous choice of model problems, even though this improvement is not possible, in general. For a basic set of one-dimensional and two-dimensional examples, we find significantly improved results using our recently introduced augmented FSSH algorithm.

The 3D photodissociation dynamics of LiNe system is investigated by quantum calculations using the multi-configuration time-dependent Hartree (MCTDH) method and by classical simulations with the trajectory surface hopping (TSH) approach. Six electronic states of A’ symmetry and two states of A” symmetry are involved in the process. Couplings in the excitation region and two conical intersections in the vicinity of the Franck–Condon zone control the non-adiabatic nuclear dynamics. A diabatic representation including all the states and the couplings is determined. Diabatic and adiabatic populations calculated for initial excitation to pure diabatic and adiabatic states lead to a clear understanding of the mechanisms governing the non-adiabatic photodissociation process. The classical and quantum photodissociation cross-sections for absorption in two adiabatic states of the A’ symmetry are calculated. A remarkable agreement between quantum and classical results is obtained regarding the populations and the absorption cross-sections.