A partially linearized spin-mapping approach for nonadiabatic dynamics. II. Analysis and comparison with related approaches

Benchmarking various nonadiabatic semiclassical mapping dynamics methods with tensor-train thermo-field dynamics

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.

Consistent Construction of the Density Matrix from Surface Hopping Trajectories

Proper construction of the density matrix based on surface hopping trajectories remains a difficult problem. Due to the well-known overcoherence in traditional surface hopping simulations, the electronic wave function cannot be used directly. In this work, we propose a consistent density matrix construction method, which takes the advantage of occupation of active states to rescale the coherence calculated by wave functions and ensures the intrinsic consistency of the density matrix. This new trajectory analysis method can be used for both Tully's fewest switches surface hopping (FSSH) and our recently proposed branching corrected surface hopping (BCSH). As benchmarked in both one- and two-dimensional standard scattering models, the new approach combined with BCSH trajectories achieves highly accurate time-dependent spatial distributions of adiabatic populations and coherence compared to exact quantum results.

Detailed balance in mixed quantum-classical mapping approaches

The violation of detailed balance poses a serious problem for the majority of current quasiclassical methods for simulating nonadiabatic dynamics. In order to analyze the severity of the problem, we predict the long-time limits of the electronic populations according to various quasiclassical mapping approaches by applying arguments from classical ergodic theory. Our analysis confirms that regions of the mapping space that correspond to negative populations, which most mapping approaches introduce in order to go beyond the Ehrenfest approximation, pose the most serious issue for reproducing the correct thermalization behavior. This is because inverted potentials, which arise from negative electronic populations entering the nuclear force, can result in trajectories unphysically accelerating off to infinity. The recently developed mapping approach to surface hopping (MASH) provides a simple way of avoiding inverted potentials while retaining an accurate description of the dynamics. We prove that MASH, unlike any other quasiclassical approach, is guaranteed to describe the exact thermalization behavior of all quantum-classical systems, confirming it as one of the most promising methods for simulating nonadiabatic dynamics in real condensed-phase systems.

Which Algorithm Best Propagates the Meyer-Miller-Stock-Thoss Mapping Hamiltonian for Non-Adiabatic Dynamics?

A common strategy to simulate mixed quantum-classical dynamics is by propagating classical trajectories with mapping variables, often using the Meyer-Miller-Stock-Thoss (MMST) Hamiltonian or the related spin-mapping approach. When mapping the quantum subsystem, the coupled dynamics reduce to a set of equations of motion to integrate. Several numerical algorithms have been proposed, but a thorough performance comparison appears to be lacking. Here, we compare three time-propagation algorithms for the MMST Hamiltonian: the Momentum Integral (MInt) (J. Chem. Phys., 2018, 148, 102326), the Split-Liouvillian (SL) (Chem. Phys., 2017, 482, 124-134), and the algorithm in J. Chem. Phys., 2012, 136, 084101 that we refer to as the Degenerate Eigenvalue (DE) algorithm due to the approximation required during derivation. We analyze the accuracy of individual trajectories, correlation functions, energy conservation, symplecticity, Liouville's theorem, and the computational cost. We find that the MInt algorithm is the only rigorously symplectic algorithm. However, comparable accuracy at a lower computational cost can be obtained with the SL algorithm. The approximation implicitly made within the DE algorithm conserves energy poorly, even for small timesteps, and thus leads to slightly different results. These results should guide future mapping-variable simulations.

A mapping approach to surface hopping

We present a nonadiabatic classical-trajectory approach that offers the best of both worlds between fewest-switches surface hopping (FSSH) and quasiclassical mapping dynamics. This mapping approach to surface hopping (MASH) propagates the nuclei on the active adiabatic potential-energy surface, such as in FSSH. However, unlike in FSSH, transitions between active surfaces are deterministic and occur when the electronic mapping variables evolve between specified regions of the electronic phase space. This guarantees internal consistency between the active surface and the electronic degrees of freedom throughout the dynamics. MASH is rigorously derivable from exact quantum mechanics as a limit of the quantum-classical Liouville equation (QCLE), leading to a unique prescription for momentum rescaling and frustrated hops. Hence, a quantum-jump procedure can, in principle, be used to systematically converge the accuracy of the results to that of the QCLE. This jump procedure also provides a rigorous framework for deriving approximate decoherence corrections similar to those proposed for FSSH. We apply MASH to simulate the nonadiabatic dynamics in various model systems and show that it consistently produces more accurate results than FSSH at a comparable computational cost.

An accurate and efficient Ehrenfest dynamics approach for calculating linear and nonlinear electronic spectra

Linear and nonlinear electronic spectra provide an important tool to probe the absorption and transfer of electronic energy. Here, we introduce a pure state Ehrenfest approach to obtain accurate linear and nonlinear spectra that is applicable to systems with large numbers of excited states and complex chemical environments. We achieve this by representing the initial conditions as sums of pure states and unfolding multi-time correlation functions into the Schrödinger picture. By doing this, we show that one can obtain significant improvements in accuracy over the previously used projected Ehrenfest approach and that these benefits are particularly pronounced in cases where the initial condition is a coherence between excited states. While such initial conditions do not arise when calculating linear electronic spectra, they play a vital role in capturing multidimensional spectroscopies. We demonstrate the performance of our method by showing that it is able to quantitatively capture the exact linear, 2D electronic spectroscopy, and pump-probe spectra for a Frenkel exciton model in slow bath regimes and is even able to reproduce the main spectral features in fast bath regimes.

On detailed balance in nonadiabatic dynamics: From spin spheres to equilibrium ellipsoids

Trajectory-based methods that propagate classical nuclei on multiple quantum electronic states are often used to simulate nonadiabatic processes in the condensed phase. A long-standing problem of these methods is their lack of detailed balance, meaning that they do not conserve the equilibrium distribution. In this article, we investigate ideas for restoring detailed balance in mixed quantum-classical systems by tailoring the previously proposed spin-mapping approach to thermal equilibrium. We find that adapting the spin magnitude can recover the correct long-time populations but is insufficient to conserve the full equilibrium distribution. The latter can however be achieved by a more flexible mapping of the spin onto an ellipsoid, which is constructed to fulfill detailed balance for arbitrary potentials. This ellipsoid approach solves the problem of negative populations that has plagued previous mapping approaches and can therefore be applied also to strongly asymmetric and anharmonic systems. Because it conserves the thermal distribution, the method can also exploit efficient sampling schemes used in standard molecular dynamics, which drastically reduces the number of trajectories needed for convergence. The dynamics does however still have mean-field character, as is observed most clearly by evaluating reaction rates in the golden-rule limit. This implies that although the ellipsoid mapping provides a rigorous framework, further work is required to find an accurate classical-trajectory approximation that captures more properties of the true quantum dynamics.

Quasiclassical approaches to the generalized quantum master equation

The formalism of the generalized quantum master equation (GQME) is an effective tool to simultaneously increase the accuracy and the efficiency of quasiclassical trajectory methods in the simulation of nonadiabatic quantum dynamics. The GQME expresses correlation functions in terms of a non-Markovian equation of motion, involving memory kernels that are typically fast-decaying and can therefore be computed by short-time quasiclassical trajectories. In this paper, we study the approximate solution of the GQME, obtained by calculating the kernels with two methods: Ehrenfest mean-field theory and spin-mapping. We test the approaches on a range of spin-boson models with increasing energy bias between the two electronic levels and place a particular focus on the long-time limits of the populations. We find that the accuracy of the predictions of the GQME depends strongly on the specific technique used to calculate the kernels. In particular, spin-mapping outperforms Ehrenfest for all the systems studied. The problem of unphysical negative electronic populations affecting spin-mapping is resolved by coupling the method with the master equation. Conversely, Ehrenfest in conjunction with the GQME can predict negative populations, despite the fact that the populations calculated from direct dynamics are positive definite.

Incorporating Lindblad Decay Dynamics into Mixed Quantum-Classical Simulations

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

2D electronic-vibrational spectroscopy with classical trajectories

Two-dimensional electronic-vibrational (2DEV) spectra have the capacity to probe electron–nuclear interactions in molecules by measuring correlations between initial electronic excitations and vibrational transitions at a later time. The trajectory-based semiclassical optimized mean trajectory approach is applied to compute 2DEV spectra for a system with excitonically coupled electronic excited states vibronically coupled to a chromophore vibration. The chromophore mode is in turn coupled to a bath, inducing redistribution of vibrational populations. The lineshapes and delay-time dynamics of the resulting spectra compare well with benchmark calculations, both at the level of the observable and with respect to contributions from distinct spectroscopic processes.

Explaining the Efficiency of Photosynthesis: Quantum Uncertainty or Classical Vibrations?

Photosynthetic organisms are known to use a mechanism of vibrationally assisted exciton energy transfer to efficiently harvest energy from light. The importance of quantum effects in this mechanism is a long-standing topic of debate, which has traditionally focused on the role of excitonic coherences. Here, we address another recent claim: that the efficient energy transfer in the Fenna-Matthews-Olson complex relies on nuclear quantum uncertainty and would not function if the vibrations were classical. We present a counter-example to this claim, showing by trajectory-based simulations that a description in terms of quantum electrons and classical nuclei is indeed sufficient to describe the funneling of energy to the reaction center. We analyze and compare these findings to previous classical-nuclear approximations that predicted the absence of an energy funnel and conclude that the key difference and the reason for the discrepancy is the ability of the trajectories to properly account for Newton's third law.

Two-dimensional vibronic spectroscopy with semiclassical thermofield dynamics

Thermofield dynamics is an exactly correct formulation of quantum mechanics at finite temperature in which a wavefunction is governed by an effective temperature-dependent quantum Hamiltonian. The optimized mean trajectory (OMT) approximation allows the calculation of spectroscopic response functions from trajectories produced by the classical limit of a mapping Hamiltonian that includes physical nuclear degrees of freedom and other effective degrees of freedom representing discrete vibronic states. Here, we develop a thermofield OMT (TF-OMT) approach in which the OMT procedure is applied to a temperature-dependent classical Hamiltonian determined from the thermofield-transformed quantum mapping Hamiltonian. Initial conditions for bath nuclear degrees of freedom are sampled from a zero-temperature distribution. Calculations of two-dimensional electronic spectra and two-dimensional vibrational–electronic spectra are performed for models that include excitonically coupled electronic states. The TF-OMT calculations agree very closely with the corresponding OMT results, which, in turn, represent well benchmark calculations with the hierarchical equations of motion method.

A partially linearized spin-mapping approach for simulating nonlinear optical spectra

We present a partially linearized method based on spin-mapping for computing both linear and nonlinear optical spectra. As observables are obtained from ensembles of classical trajectories, the approach can be applied to the large condensed-phase systems that undergo photosynthetic light-harvesting processes. In particular, the recently derived spin partially linearized density matrix method has been shown to exhibit superior accuracy in computing population dynamics compared to other related classical-trajectory methods. Such a method should also be ideally suited to describing the quantum coherences generated by interaction with light. We demonstrate that this is, indeed, the case by calculating the nonlinear optical response functions relevant for the pump-probe and 2D photon-echo spectra for a Frenkel biexciton model and the Fenna-Matthews-Olsen light-harvesting complex. One especially desirable feature of our approach is that the full spectrum can be decomposed into its constituent components associated with the various Liouville-space pathways, offering a greater insight beyond what can be directly obtained from experiments.

Quantum Entanglement from Classical Trajectories

A long-standing challenge in mixed quantum-classical trajectory simulations is the treatment of entanglement between the classical and quantal degrees of freedom. We present a novel approach that describes the emergence of entangled states entirely in terms of independent and deterministic Ehrenfest-like classical trajectories. For a two-level quantum system in a classical environment, this is derived by mapping the quantum system onto a path-integral representation of a spin 1/2. We demonstrate that the method correctly accounts for coherence and decoherence and thus reproduces the splitting of a wave packet in a nonadiabatic scattering problem. This discovery opens up a new class of simulations as an alternative to stochastic surface-hopping, coupled-trajectory, or semiclassical approaches.

Generalized discrete truncated Wigner approximation for nonadiabatic quantum-classical dynamics

Nonadiabatic molecular dynamics occur in a wide range of chemical reactions and femtochemistry experiments involving electronically excited states. These dynamics are hard to treat numerically as the system's complexity increases, and it is thus desirable to have accurate yet affordable methods for their simulation. Here, we introduce a linearized semiclassical method, the generalized discrete truncated Wigner approximation (GDTWA), which is well-established in the context of quantum spin lattice systems, into the arena of chemical nonadiabatic systems. In contrast to traditional continuous mapping approaches, e.g., the Meyer-Miller-Stock-Thoss and the spin mappings, GDTWA samples the electron degrees of freedom in a discrete phase space and thus forbids an unphysical unbounded growth of electronic state populations. The discrete sampling also accounts for an effective reduced but non-vanishing zero-point energy without an explicit parameter, which makes it possible to treat the identity operator and other operators on an equal footing. As numerical benchmarks on two linear vibronic coupling models and Tully's models show, GDTWA has a satisfactory accuracy in a wide parameter regime, independent of whether the dynamics is dominated by relaxation or by coherent interactions. Our results suggest that the method can be very adequate to treat challenging nonadiabatic dynamics problems in chemistry and related fields.

Non-adiabatic Matsubara dynamics and non-adiabatic ring-polymer molecular dynamics

We present the non-adiabatic Matsubara dynamics, a general framework for computing the time-correlation function (TCF) of electronically non-adiabatic systems. This new formalism is derived based on the generalized Kubo-transformed TCF using the Wigner representation for both the nuclear degrees of freedom and the electronic mapping variables. By dropping the non-Matsubara nuclear normal modes in the quantum Liouvillian and explicitly integrating these modes out from the expression of the TCF, we derived the non-adiabatic Matsubara dynamics approach. Further making the approximation to drop the imaginary part of the Matsubara Liouvillian and enforce the nuclear momentum integral to be real, we arrived at the non-adiabatic ring-polymer molecular dynamics (NRPMD) approach. We have further justified the capability of NRPMD for simulating the non-equilibrium TCF. This work provides the rigorous theoretical foundation for several recently proposed state-dependent RPMD approaches and offers a general framework for developing new non-adiabatic quantum dynamics methods in the future.