Trajectory-adjusted electronic zero point energy in classical Meyer-Miller vibronic dynamics: Symmetrical quasiclassical application to photodissociation

Unsupervised Machine Learning in the Analysis of Nonadiabatic Molecular Dynamics Simulation

The all-atomic full-dimensional-level simulations of nonadiabatic molecular dynamics (NAMD) in large realistic systems has received high research interest in recent years. However, such NAMD simulations normally generate an enormous amount of time-dependent high-dimensional data, leading to a significant challenge in result analyses. Based on unsupervised machine learning (ML) methods, considerable efforts were devoted to developing novel and easy-to-use analysis tools for the identification of photoinduced reaction channels and the comprehensive understanding of complicated molecular motions in NAMD simulations. Here, we tried to survey recent advances in this field, particularly to focus on how to use unsupervised ML methods to analyze the trajectory-based NAMD simulation results. Our purpose is to offer a comprehensive discussion on several essential components of this analysis protocol, including the selection of ML methods, the construction of molecular descriptors, the establishment of analytical frameworks, their advantages and limitations, and persistent challenges.

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.

Nonadiabatic Field with Triangle Window Functions on Quantum Phase Space

Recent progress on the constraint coordinate-momentum phase space (CPS) formulation of finite-state quantum systems has revealed that the triangle window function approach is an isomorphic representation of the exact population-population correlation function of the two-state system. We use the triangle window (TW) function and the CPS mapping kernel element to formulate a novel useful representation of discrete electronic degrees of freedom (DOFs). When it is employed with nonadiabatic field (NaF) dynamics, a new variant of the NaF approach (i.e., NaF-TW) is proposed. The NaF-TW expression of the population of any adiabatic state is always positive semidefinite. Extensive benchmark tests of model systems in both the condensed phase and gas phase demonstrate that the NaF-TW approach is able to faithfully capture the dynamical interplay between electronic and nuclear DOFs in a broad region, including where the states remain coupled all the time, as well as where the bifurcation characteristic of nuclear motion is important.

Resonance Enhancement of Vibrational Polariton Chemistry Obtained from the Mixed Quantum-Classical Dynamics Simulations

We applied a variety of mixed quantum-classical (MQC) approaches to simulate the VSC-influenced reaction rate constant. All of these MQC simulations treat the key vibrational levels associated with the reaction coordinate in the quantum subsystem (as quantum states), whereas all other degrees of freedom (DOFs) are treated inside the classical subsystem. We find that, as long as we have the quantum state descriptions for the vibrational DOFs, one can correctly describe the VSC resonance condition when the cavity frequency matches the bond vibrational frequency. This correct resonance behavior can be obtained regardless of the detailed MQC methods that one uses. The results suggest that the MQC approaches can generate semiquantitative agreement with the exact results for rate constant changes when changing the cavity frequency, the light-matter coupling strength, or the cavity lifetime. The finding of this work suggests that one can use computationally economic MQC approaches to explore the collective coupling scenario when many molecules are collectively coupled to many cavity modes in the future.

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.

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.

Real-Time Dynamics and Detailed Balance in Ring Polymer Surface Hopping: The Impact of Frustrated Hops

Ring polymer surface hopping (RPSH) has been recently introduced as a well-tailored method for incorporating nuclear quantum effects, such as zero-point energy and tunneling, into nonadiabatic molecular dynamics simulations. The practical widespread usage of RPSH demands a comprehensive benchmarking of different reaction regimes and conditions with equal emphasis on demonstrating both the cons and the pros of the method. Here, we investigate the fundamental questions related to the conservation of energy and detailed balance in the context of RPSH. Using Tully's avoided crossing model as well as a 2-state quantum system coupled to a classical bath undergoing Langevin dynamics, we probe the critical problem of the proper treatment of the classically forbidden transitions stemming from the surface hopping algorithm. We show that proper treatment of these frustrated hops is key to the accurate description of real-time dynamics as well as reproducing the correct quantum Boltzmann populations.

Learning Decoherence Time Formulas for Surface Hopping from Quantum Dynamics

Surface hopping simulations have achieved great success in many different fields, but their reliability has long been limited by the overcoherence problem. We here present a general machine learning assisted approach to identify optimal decoherence time formulas for surface hopping using exact quantum dynamics as references. In order to avoid computationally expensive force calculations, we use the nuclear kinetic energy and the adiabatic energy difference to iteratively generate the descriptor space. Through multilayer screening of the candidate descriptors and discrete optimization of the relevant parameters, we obtain new energy-based decoherence time formulas. As benchmarked in thousands of diverse multilevel systems and six standard scattering models, surface hopping with our new decoherence time formulas nearly reproduces the exact quantum dynamics while maintaining high efficiency. Thereby, our approach provides a promising avenue for systematically improving the accuracy of surface hopping simulations in complex systems from quantum dynamics data.

Theoretical Advances in Polariton Chemistry and Molecular Cavity Quantum Electrodynamics

When molecules are coupled to an optical cavity, new light-matter hybrid states, so-called polaritons, are formed due to quantum light-matter interactions. With the experimental demonstrations of modifying chemical reactivities by forming polaritons under strong light-matter interactions, theorists have been encouraged to develop new methods to simulate these systems and discover new strategies to tune and control reactions. This review summarizes some of these exciting theoretical advances in polariton chemistry, in methods ranging from the fundamental framework to computational techniques and applications spanning from photochemistry to vibrational strong coupling. Even though the theory of quantum light-matter interactions goes back to the midtwentieth century, the gaps in the knowledge of molecular quantum electrodynamics (QED) have only recently been filled. We review recent advances made in resolving gauge ambiguities, the correct form of different QED Hamiltonians under different gauges, and their connections to various quantum optics models. Then, we review recently developed ab initio QED approaches which can accurately describe polariton states in a realistic molecule-cavity hybrid system. We then discuss applications using these method advancements. We review advancements in polariton photochemistry where the cavity is made resonant to electronic transitions to control molecular nonadiabatic excited state dynamics and enable new photochemical reactivities. When the cavity resonance is tuned to the molecular vibrations instead, ground-state chemical reaction modifications have been demonstrated experimentally, though its mechanistic principle remains unclear. We present some recent theoretical progress in resolving this mystery. Finally, we review the recent advances in understanding the collective coupling regime between light and matter, where many molecules can collectively couple to a single cavity mode or many cavity modes. We also lay out the current challenges in theory to explain the observed experimental results. We hope that this review will serve as a useful document for anyone who wants to become familiar with the context of polariton chemistry and molecular cavity QED and thus significantly benefit the entire community.

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.

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.

Trajectory Propagation of Symmetrical Quasi-classical Dynamics with Meyer-Miller Mapping Hamiltonian Using Machine Learning

The long short-term memory recurrent neural network (LSTM-RNN) approach is applied to realize the trajectory-based nonadiabatic dynamics within the framework of the symmetrical quasi-classical dynamics method based on the Meyer-Miller mapping Hamiltonian (MM-SQC). After construction, the LSTM-RNN model allows us to propagate the entire trajectory evolutions of all involved degrees of freedoms (DOFs) from initial conditions. The proposed idea is proven to be reliable and accurate in the simulations of the dynamics of several site-exciton electron-phonon coupling models and three Tully's scattering models. It indicates that the LSTM-RNN model perfectly captures the dynamical information on the trajectory evolution in the MM-SQC dynamics. Our work proposes a novel machine learning approach in the simulation of trajectory-based nonadiabatic dynamic of complex systems with a large number of DOFs.

All-Atom Nonadiabatic Semiclassical Mapping Dynamics for Photoinduced Charge Transfer of Organic Photovoltaic Molecules in Explicit Solvents

Direct all-atom simulation of nonadiabatic dynamics in disordered condensed phases like liquid solutions and amorphous solids has been challenging. The first all-atom simulation of the photoinduced charge-transfer dynamics of a prototypical organic photovoltaic carotenoid-porphyrin-C₆₀ molecular triad in explicit tetrahydrofuran is presented. Based on the Meyer-Miller mapping Hamiltonian, various semiclassical and mixed quantum-classical dynamics are employed, including the linearized semiclassical, symmetrical quasiclassical, mean-field Ehrenfest, classical mapping model, and spin-mapping model approaches. The all-atom nonadiabatic dynamics were compared to multi-state harmonic models with a globally shared bath, and the models built using the ensemble averages on the initial electronic state could reproduce the all-atom results. The solvent effect was found to be critical for the photoinduced charge transfer, and the time-dependent solute-solvent radial distribution functions revealed that only the nonadiabatic dynamics started with the effective forces on the initial electronic state could capture the correct nuclear dynamics. The proposed strategy for modeling condensed-phase nonadiabatic dynamics with atomistic details is readily applied to complex condensed-phase systems.

Transition Path Flight Times and Nonadiabatic Electronic Transitions

Transition path flight times are studied for scattering on two electronic surfaces with a single crossing. These flight times reveal nontrivial quantum effects such as resonance lifetimes and nonclassical passage times and reveal that nonadiabatic effects often increase flight times. The flight times are computed using numerically exact time propagation and compared with results obtained from the Fewest Switches Surface Hopping (FSSH) method. Comparison of the two methods shows that the FSSH method is reliable for transition path times only when the scattering is classically allowed on the relevant adiabatic surfaces. However, where quantum effects such as tunneling and resonances dominate, the FSSH method is not adequate to accurately predict the correct times and transition probabilities. These results highlight limitations in methods which do not account for quantum interference effects, and suggest that measuring flight times is important for obtaining insights from the time-domain into quantum effects in nonadiabatic scattering.

Non-adiabatic Mapping Dynamics in the Phase Space of the ${SU}(N)$ Lie Group

We present the rigorous theoretical framework of the generalized spin mapping representation for non-adiabatic dynamics. Our work is based up a new mapping formalism recently introduced by Runeson and Richardson in [J. Chem. Phys. 152, 084110 (2020)], which uses the generators of the su(N) Lie algebra to represent N discrete electronic states, thus preserving the size of the original Hilbert space. Following this interesting idea, the Stratonovich-Weyl transform is used to map an operator in the Hilbert space to a continuous function on the SU(N) Lie group, i.e., a smooth manifold which is a phase space of continuous variables. We further use the Wigner representation to describe the nuclear degrees of freedom, and derived an exact expression of the time-correlation function as well as the exact quantum Liouvillian for the non-adiabatic system. Making the linearization approximation, this exact Liouvillian is reduced to the Liouvillian of several recently proposed methods, and the performance of this Linearized method is tested using non-adiabatic models. We envision that the theoretical work presented here provides a rigorous and unified framework to formally derive non-adiabatic quantum dynamics approaches with continuous variables and connect the previous methods in a clear and concise manner.

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.

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.

Dynamic signatures of electronically nonadiabatic coupling in sodium hydride: a rigorous test for the symmetric quasi-classical model applied to realistic, ab initio electronic states in the adiabatic representation

Sodium hydride (NaH) in the gas phase presents a seemingly simple electronic structure making it a potentially tractable system for the detailed investigation of nonadiabatic molecular dynamics from both computational and experimental standpoints. The single vibrational degree of freedom, as well as the strong nonadiabatic coupling that arises from the excited electronic states taking on considerable ionic character, provides a realistic chemical system to test the accuracy of quasi-classical methods to model population dynamics where the results are directly comparable against quantum mechanical benchmarks. Using a simulated pump-probe type experiment, this work presents computational predictions of population transfer through the avoided crossings of NaH via symmetric quasi-classical Meyer-Miller (SQC/MM), Ehrenfest, and exact quantum dynamics on realistic, ab initio potential energy surfaces. The main driving force for population transfer arises from the ground vibrational level of the D¹Σ⁺ adiabatic state that is embedded in the manifold of near-dissociation C¹Σ⁺ vibrational states. When coupled through a sharply localized first-order derivative coupling most of the population transfers between t = 15 and t = 30 fs depending on the initially excited vibronic wavepacket. While quantum mechanical effects are expected due to the reduced mass of NaH, predictions of the population dynamics from both the SQC/MM and Ehrenfest models perform remarkably well against the quantum dynamics benchmark. Additionally, an analysis of the vibronic structure in the nonadiabatically coupled regime is presented using a variational eigensolver methodology.

Path Integrals for Nonadiabatic Dynamics: Multistate Ring Polymer Molecular Dynamics

This review focuses on a recent class of path-integral-based methods that simulate nonadiabatic dynamics in the condensed phase using only classical molecular dynamics trajectories in an extended phase space. Specifically, a semiclassical mapping protocol is used to derive an exact, continuous, Cartesian variable path-integral representation for the canonical partition function of a system in which multiple electronic states are coupled to nuclear degrees of freedom. Building on this exact statistical foundation, multistate ring polymer molecular dynamics methods are developed for the approximate calculation of real-time thermal correlation functions. The remarkable promise of these multistate ring polymer methods, their successful applications, and their limitations are discussed in detail.Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Calculating Multidimensional Optical Spectra from Classical Trajectories

Multidimensional optical spectra are measured from the response of a material system to a sequence of laser pulses and have the capacity to elucidate specific molecular interactions and dynamics whose influences are absent or obscured in a conventional linear absorption spectrum. Interpretation of complex spectra is supported by theoretical modeling of the spectroscopic observable, requiring implementation of quantum dynamics for coupled electrons and nuclei. Performing numerically correct quantum dynamics in this context may pose computational challenges, particularly in the condensed phase. Semiclassical methods based on calculating classical trajectories offer a practical alternative. Here I review the recent application of some semiclassical, trajectory-based methods to nonlinear molecular vibrational and electronic spectra.

Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Unified Formulation of Phase Space Mapping Approaches for Nonadiabatic Quantum Dynamics

Nonadiabatic dynamical processes are one of the most important quantum mechanical phenomena in chemical, materials, biological, and environmental molecular systems, where the coupling between different electronic states is either inherent in the molecular structure or induced by the (intense) external field. The curse of dimensionality indicates the intractable exponential scaling of calculation effort with system size and restricts the implementation of "numerically exact" approaches for realistic large systems. The phase space formulation of quantum mechanics offers an important theoretical framework for constructing practical approximate trajectory-based methods for quantum dynamics. This Account reviews our recent progress in phase space mapping theory: a unified framework for constructing the mapping Hamiltonian on phase space for coupled F-state systems where the renowned Meyer-Miller Hamiltonian model is a special case, a general phase space formulation of quantum mechanics for nonadiabatic systems where the electronic degrees of freedom are mapped onto constraint space and the nuclear degrees of freedom are mapped onto infinite space, and an isomorphism between the mapping phase space approach for nonadiabatic systems and that for nonequilibrium electron transport processes. While the zero-point-energy parameter is conventionally assumed to be positive, we show that the constraint implied in the conventional Meyer-Miller mapping Hamiltonian requires that such a parameter can be negative as well and lies in (-1/F, +∞) for each electronic degree of freedom. More importantly, the zero-point-energy parameter should be interpreted as a special case of a commutator matrix in the comprehensive phase space mapping Hamiltonian for nonadiabatic systems. From the rigorous formulation of mapping phase space, we propose approximate but practical trajectory-based nonadiabatic dynamics methods. The applications to both gas phase and condensed phase problems include the spin-boson model for condensed phase dissipative two-state systems, the three-state photodissociation models, the seven-site model of the Fenna-Matthews-Olson monomer in photosynthesis of green sulfur bacteria, the strongly coupled molecular/atomic matter-optical cavity systems designed for controlling and manipulating chemical dynamical processes, and the Landauer model for a quantum dot state coupled with two electrodes. In these applications the overall performance of our phase space mapping dynamics approach is superior to two prevailing trajectory-based methods, Ehrenfest dynamics and fewest switches surface hopping.

Multi-state harmonic models with globally shared bath for nonadiabatic dynamics in the condensed phase

Model Hamiltonians constructed from quantum chemistry calculations and molecular dynamics simulations are widely used for simulating nonadiabatic dynamics in the condensed phase. The most popular two-state spin-boson model could be built by mapping the all-atom anharmonic Hamiltonian onto a two-level system bilinearly coupled to a harmonic bath using the energy gap time correlation function. However, for more than two states, there lacks a general strategy to construct multi-state harmonic (MSH) models since the energy gaps between different pairs of electronic states are not entirely independent and need to be considered consistently. In this paper, we extend the previously proposed approach for building three-state harmonic models for photoinduced charge transfer to the arbitrary number of electronic states with a globally shared bath and the system-bath couplings are scaled differently according to the reorganization energies between each pair of states. We demonstrate the MSH model construction for an organic photovoltaic carotenoid-porphyrin-C₆₀ molecular triad dissolved in explicit tetrahydrofuran solvent. Nonadiabatic dynamics was simulated using mixed quantum-classical techniques, including the linearized semiclassical and symmetrical quasiclassical dynamics with the mapping Hamiltonians, mean-field Ehrenfest, and mixed quantum-classical Liouville dynamics in two-state, three-state, and four-state harmonic models of the triad system. The MSH models are shown to provide a general and flexible framework for simulating nonadiabatic dynamics in complex systems.

Ab initio symmetric quasi-classical approach to investigate molecular Tully models

We perform on-the-fly non-adiabatic molecular dynamics simulations using the symmetrical quasi-classical (SQC) approach with the recently suggested molecular Tully models: ethylene and fulvene. We attempt to provide benchmarks of the SQC methods using both the square and triangle windowing schemes as well as the recently proposed electronic zero-point-energy correction scheme (the so-called γ correction). We use the quasi-diabatic propagation scheme to directly interface the diabatic SQC methods with adiabatic electronic structure calculations. Our results showcase the drastic improvement of the accuracy by using the trajectory-adjusted γ-corrections, which outperform the widely used trajectory surface hopping method with decoherence corrections. These calculations provide useful and non-trivial tests to systematically investigate the numerical performance of various diabatic quantum dynamics approaches, going beyond simple diabatic model systems that have been used as the major workhorse in the quantum dynamics field. At the same time, these available benchmark studies will also likely foster the development of new quantum dynamics approaches based on these techniques.

Commutator Matrix in Phase Space Mapping Models for Nonadiabatic Quantum Dynamics

We show that a novel, general phase space mapping Hamiltonian for nonadiabatic systems, which is reminiscent of the renowned Meyer-Miller mapping Hamiltonian, involves a commutator variable matrix rather than the conventional zero-point-energy parameter. In the exact mapping formulation on constraint space for phase space approaches for nonadiabatic dynamics, the general mapping Hamiltonian with commutator variables can be employed to generate approximate trajectory-based dynamics. Various benchmark model tests, which range from gas phase to condensed phase systems, suggest that the overall performance of the general mapping Hamiltonian is better than that of the conventional Meyer-Miller Hamiltonian.

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

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

Two-dimensional vibrational-electronic spectra with semiclassical mechanics

Two-dimensional vibrational-electronic (2DVE) spectra probe the effects on vibronic spectra of initial vibrational excitation in an electronic ground state. The optimized mean trajectory (OMT) approximation is a semiclassical method for computing nonlinear spectra from response functions. Ensembles of classical trajectories are subject to semiclassical quantization conditions, with the radiation-matter interaction inducing discontinuous transitions. This approach has been previously applied to two-dimensional infrared and electronic spectra and is extended here to 2DVE spectra. For a system including excitonic coupling, vibronic coupling, and interaction of a chromophore vibration with a resonant environment, the OMT method is shown to well approximate exact quantum dynamics.

Negative Zero-Point-Energy Parameter in the Meyer-Miller Mapping Model for Nonadiabatic Dynamics

The celebrated Meyer-Miller mapping model has been a useful approach for generating practical trajectory-based nonadiabatic dynamics methods. It is generally assumed that the zero-point-energy (ZPE) parameter is positive. The constraint implied in the conventional Meyer-Miller mapping Hamiltonian for an F-electronic-state system actually requires γ∈(-1/F, ∞) for the ZPE parameter for each electronic degree of freedom. Both negative and positive values are possible for such a parameter. We first establish a rigorous formulation to construct exact mapping models in the Cartesian phase space when the constraint is applied. When nuclear dynamics is approximated by the linearized semiclassical initial value representation, a negative ZPE parameter could lead to reasonably good performance in describing dynamic behaviors in typical spin-boson models for condensed-phase two-state systems, even at challenging zero temperature.

Spectroscopic response theory with classical mapping Hamiltonians

Exact quantum dynamics with a time-independent Hamiltonian in a discrete state space can be computed using classical mechanics through the classical Meyer-Miller-Stock-Thoss mapping Hamiltonian. In order to compute quantum response functions from classical dynamics, we extend this mapping to a quantum Hamiltonian with time-dependence arising from a classical field. This generalization requires attention to time-ordering in quantum and classical propagators. Quantum response theory with the original quantum Hamiltonian is equivalent to classical response theory with the classical mapping Hamiltonian. We elucidate the structure of classical response theory with the mapping Hamiltonian, thereby generating classical versions of the two-sided quantum density operator diagrams conventionally used to describe spectroscopic processes. This formal development can provide a foundation for new semiclassical approximations to spectroscopic observables for models in which classical nuclear degrees of freedom are introduced into a mapping Hamiltonian describing electronic states. Calculations of the temperature-dependence of two-dimensional electronic spectra for an exciton dimer using two semiclassical approaches are compared with benchmark calculations using the hierarchical equations of motion method.

Study of the exciton dynamics in perylene bisimide (PBI) aggregates with symmetrical quasiclassical dynamics based on the Meyer–Miller mapping Hamiltonian

The exciton dynamics in one-dimensional stacked PBI (Perylene Bisimide) aggregates was studied with SQC-MM dynamics (Symmetrical Quasiclassical Dynamics based on the Meyer–Miller mapping Hamiltonian).

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.

Quasi-Diabatic Scheme for Nonadiabatic On-the-Fly Simulations

We use the quasi-diabatic (QD) propagation scheme to perform on-the-fly nonadiabatic simulations of the photodynamics of ethylene. The QD scheme enables a seamless interface between accurate diabatic-based quantum dynamics approaches and adiabatic electronic structure calculations, explicitly avoiding any efforts to construct global diabatic states or reformulate the diabatic dynamics approach to the adiabatic representation. Using the partial linearized path-integral approach and the symmetrical quasi-classical approach as the diabatic dynamics methods, the QD propagation scheme enables direct nonadiabatic simulation with complete active space self-consistent field on-the-fly electronic structure calculations. The population dynamics obtained from both approaches are in a close agreement with the quantum wavepacket-based method and outperform the widely used trajectory surface-hopping approach. Further analysis of the ethylene photodeactivation pathways demonstrates the correct predictions of competing processes of nonradiative relaxation mechanism through various conical intersections. This work provides the foundation of using accurate diabatic dynamics approaches and on-the-fly adiabatic electronic structure information to perform ab initio nonadiabatic simulation.