A Nonperturbative Methodology for Simulating Multidimensional Spectra of Multiexcitonic Molecular Systems via Quasiclassical Mapping Hamiltonian Methods

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.

Coarse-Grained Approach to Simulate Signatures of Excitation Energy Transfer in Two-Dimensional Electronic Spectroscopy of Large Molecular Systems

Two-dimensional electronic spectroscopy (2DES) has proven to be a highly effective technique in studying the properties of excited states and the process of excitation energy transfer in complex molecular assemblies, particularly in biological light-harvesting systems. However, the accurate simulation of 2DES for large systems still poses a challenge because of the heavy computational demands it entails. In an effort to overcome this limitation, we devised a coarse-grained 2DES method. This method encompasses the treatment of the entire system by dividing it into distinct weakly coupled segments, which are assumed to communicate predominantly through incoherent exciton transfer. We first demonstrate the efficiency of this method through simulation on a model dimer system, which demonstrates a marked improvement in calculation efficiency, with results that exhibit good concordance with reference spectra calculated with less approximate methods. Additionally, the application of this method to the light-harvesting antenna 2 (LH2) complex of purple bacteria showcases its advantages, accuracy, and limitations. Furthermore, simulating the anisotropy decay in LH2 induced by energy transfer and its comparison with experiments confirm that the method is capable of accurately describing dynamical processes in a biologically relevant system. This method presented lends itself to an extension that accounts for the effect of intrasegment relaxation processes on the 2DES spectra, which for computational efficiency are ignored in the implementation reported here. It is envisioned that the method will be employed in the future to accurately and efficiently calculate 2D spectra of more extensive systems, such as photosynthetic supercomplexes.

Generalized quantum master equations can improve the accuracy of semiclassical predictions of multitime correlation functions

Multitime quantum correlation functions are central objects in physical science, offering a direct link between the experimental observables and the dynamics of an underlying model. While experiments such as 2D spectroscopy and quantum control can now measure such quantities, the accurate simulation of such responses remains computationally expensive and sometimes impossible, depending on the system's complexity. A natural tool to employ is the generalized quantum master equation (GQME), which can offer computational savings by extending reference dynamics at a comparatively trivial cost. However, dynamical methods that can tackle chemical systems with atomistic resolution, such as those in the semiclassical hierarchy, often suffer from poor accuracy, limiting the credence one might lend to their results. By combining work on the accuracy-boosting formulation of semiclassical memory kernels with recent work on the multitime GQME, here we show for the first time that one can exploit a multitime semiclassical GQME to dramatically improve both the accuracy of coarse mean-field Ehrenfest dynamics and obtain orders of magnitude efficiency gains.

Efficient formulation of multitime generalized quantum master equations: Taming the cost of simulating 2D spectra

Modern 4-wave mixing spectroscopies are expensive to obtain experimentally and computationally. In certain cases, the unfavorable scaling of quantum dynamics problems can be improved using a generalized quantum master equation (GQME) approach. However, the inclusion of multiple (light-matter) interactions complicates the equation of motion and leads to seemingly unavoidable cubic scaling in time. In this paper, we present a formulation that greatly simplifies and reduces the computational cost of previous work that extended the GQME framework to treat arbitrary numbers of quantum measurements. Specifically, we remove the time derivatives of quantum correlation functions from the modified Mori-Nakajima-Zwanzig framework by switching to a discrete-convolution implementation inspired by the transfer tensor approach. We then demonstrate the method's capabilities by simulating 2D electronic spectra for the excitation-energy-transfer dimer model. In our method, the resolution of data can be arbitrarily coarsened, especially along the t2 axis, which mirrors how the data are obtained experimentally. Even in a modest case, this demands O(103) fewer data points. We are further able to decompose the spectra into one-, two-, and three-time correlations, showing how and when the system enters a Markovian regime where further measurements are unnecessary to predict future spectra and the scaling becomes quadratic. This offers the ability to generate long-time spectra using only short-time data, enabling access to timescales previously beyond the reach of standard methodologies.

Mixed Quantum-Classical Liouville Equation Treatment of Electronic Spectroscopy of Condensed Systems: Harmonic and Anharmonic Electron-Phonon Coupling

This Review integrates the use of electronic optical response function theory and the mixed quantum-classical (MQC) Liouville equation (MQCLE), thereby leading to electronic spectroscopy in MQC media. It further sheds light on the applicability, utility, and efficiency of the mixed quantum-classical dynamics (MQCD) formalism, which starts off with the MQCLE, in probing spectroscopy and dynamics of condensed systems, whereby quantum mechanics and classical mechanics are combined systematically. The author has been exploring and implementing MQCD to investigate electron-phonon coupling effects on electronic dephasing in harmonic and anharmonic systems by calculating linear and nonlinear optical transition analytically and numerically dipole moment time correlation functions in an MQC environment, thereby presenting an in depth spectral profile analysis and their shape and symmetry. The distinctive capability of the MQC time correlation functions is that ergodicity and stationarity properties are inherently satisfied as part of the mixed quantum-classical dynamics (MQCD) framework, unlike classical correlation functions. While some research groups have applied MQCLE to calculate vibrational spectra to study hydrogen-bonded complexes in a MQC environment and other groups calculated Optical response function to probe electron transfer dynamics using the basis mapping technique, the approach, purpose, rigor, applications, and path to the end results reported herein are different. Finally, the same framework is employed to study dissipative systems in the MQC limit, whereby the zero-phonon line adopts the correct width and eliminates its asymmetry. While the full quantum mechanical model, like the multimode Brownian oscillator (MBO) model, yields the correct width and inaccurate shape in the low-temperature limit, the MQCD formalism seems to produce an accurate zero-phonon profile. Nonlinear optical signals are also reviewed in MQC media to show the applicability and utility of this approach. The vibronic optical response functions developed here will account for geometry change, frequency change, and anharmonicity upon electronic excitation to accurately probe electronic dephasing, electron-phonon coupling, shape, and symmetry of profiles and present differences and similarities to the MBO model on pure electronic dephasing. Frequency change and anharmonicity are vitally crucial for accurately assessing electron-phonon coupling upon electronic excitation. This is an additional unique result obtained by the author to further demonstrate the applicability and utility of this approach over other approximation schemes in probing electronic dephasing including that of the MBO model.

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.

Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals

Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.

Effects of Heterogeneous Protein Environment on Excitation Energy Transfer Dynamics in the Fenna-Matthews-Olson Complex

The Fenna-Matthews-Olson (FMO) complex of green sulfur bacteria has been serving as a prototypical light-harvesting protein for studying excitation energy transfer (EET) dynamics in photosynthesis. The most widely used Frenkel exciton model for FMO complex assumes that each excited bacteriochlorophyll site couples to an identical and isolated harmonic bath, which does not account for the heterogeneous local protein environment. To better describe the realistic environment, we propose to use the recently developed multistate harmonic (MSH) model, which contains a globally shared bath that couples to the different pigment sites according to the atomistic quantum mechanics/molecular mechanics simulations with explicit protein scaffold and solvent. In this work, the effects of heterogeneous protein environment on EET in FMO complexes from Prosthecochloris aestuarii and Chlorobium tepidum, specifically including realistic spectral density, site-dependent reorganization energies, and system-bath couplings are investigated. Semiclassical and mixed quantum-classical mapping dynamics were applied to obtain the nonadiabatic EET dynamics in several models ranging from the Frenkel exciton model to the MSH model and their variants. The MSH model with realistic spectral density and site-dependent system-bath couplings displays slower EET dynamics than the Frenkel exciton model. Our comparative study shows that larger average reorganization energy, heterogeneity in spectral densities, and low-frequency modes could facilitate energy dissipation, which is insensitive to the static disorder in reorganization energies. The effects of the spectral densities and system-bath couplings along with the MSH model can be used to optimize EET dynamics for artificial light-harvesting systems.

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.

Imaginary-Time Open-Chain Path-Integral Approach for Two-State Time Correlation Functions and Applications in Charge Transfer

Quantum time correlation functions (TCFs) involving two states are important for describing nonadiabatic dynamical processes such as charge transfer. Based on a previous single-state method, we propose an imaginary-time open-chain path-integral (OCPI) approach for evaluating the two-state symmetrized TCFs. Expressing the forward and backward propagation on different electronic potential energy surfaces as a complex-time path integral, we then transform the path variables to average and difference variables such that the integration over the difference variables up to the second order can be performed analytically. The resulting expression for the symmetrized TCF is equivalent to sampling the open-chain configurations in an effective potential that corresponds to the average surface. Using importance sampling over the extended OCPI space via open path integral molecular dynamics, we tested the resulting path-integral approximation by calculating the Fermi's golden rule charge transfer rate constant within a widely-used spin-boson model. Comparing with the real-time linearized semiclassical method and analytical result, we show that the imaginary-time OCPI provides an accurate two-state symmetrized TCF and rate constant in the typical turnover region. It is shown that the first bead of the open chain corresponds to physical zero-time, and the endpoint bead corresponds to final time t; oscillations of the end-to-end distance perfectly match the nuclear mode frequency. The two-state OCPI scheme is seen to capture the tested model's electronic quantum coherence and nuclear quantum effects accurately.

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.

Equilibrium-nonequilibrium ring-polymer molecular dynamics for nonlinear spectroscopy

Two-dimensional Raman and hybrid terahertz-Raman spectroscopic techniques provide invaluable insight into molecular structures and dynamics of condensed-phase systems. However, corroborating experimental results with theory is difficult due to the high computational cost of incorporating quantum-mechanical effects in the simulations. Here, we present the equilibrium-nonequilibrium ring-polymer molecular dynamics (RPMD), a practical computational method that can account for nuclear quantum effects on the two-time response function of nonlinear optical spectroscopy. Unlike a recently developed approach based on the double Kubo transformed (DKT) correlation function, our method is exact in the classical limit, where it reduces to the established equilibrium-nonequilibrium classical molecular dynamics method. Using benchmark model calculations, we demonstrate the advantages of the equilibrium-nonequilibrium RPMD over classical and DKT-based approaches. Importantly, its derivation, which is based on the nonequilibrium RPMD, obviates the need for identifying an appropriate Kubo transformed correlation function and paves the way for applying real-time path-integral techniques to multidimensional spectroscopy.

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.

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.

Ultrafast Internal Conversion Dynamics through the on-the-Fly Simulation of Transient Absorption Pump-Probe Spectra with Different Electronic Structure Methods

An on-the-fly surface-hopping simulation protocol is developed for the evaluation of transient absorption (TA) pump-probe (PP) signals of molecular systems exhibiting internal conversion to the electronic ground state. We study the nonadiabatic dynamics of azomethane and the associating TA PP spectra at three levels of the electronic-structure theory, OM2/MRCI, SA-CASSCF, and XMS-CASPT2. The impact of these methods on the population dynamics and time-resolved TA PP signals is substantially different. This difference is attributed to the strong non-Condon effects that must be taken into account for the proper understanding and interpretation of time-resolved TA PP signals of nonadiabatic polyatomic systems. This shows that the combination of the dynamical and spectral simulations definitely provides more accurate and detailed information on the microscopic mechanisms of photophysical and photochemical processes. Hence the simulation of time-resolved spectroscopic signals provides another important dimension to examine the accuracy of quantum chemistry methods.

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.

Forecasting nonadiabatic dynamics using hybrid convolutional neural network/long short-term memory network

Modeling nonadiabatic dynamics in complex molecular or condensed-phase systems has been challenging, especially for the long-time dynamics. In this work, we propose a time series machine learning scheme based on the hybrid convolutional neural network/long short-term memory (CNN-LSTM) framework for predicting the long-time quantum behavior, given only the short-time dynamics. This scheme takes advantage of both the powerful local feature extraction ability of CNN and the long-term global sequential pattern recognition ability of LSTM. With feature fusion of individually trained CNN-LSTM models for the quantum population and coherence dynamics, the proposed scheme is shown to have high accuracy and robustness in predicting the linearized semiclassical and symmetrical quasiclassical mapping dynamics as well as the mixed quantum-classical Liouville dynamics of various spin-boson models with learning time up to 0.3 ps. Furthermore, if the hybrid network has learned the dynamics of a system, this knowledge is transferable that could significantly enhance the accuracy in predicting the dynamics of a similar system. The hybrid CNN-LSTM network is thus believed to have high predictive power in forecasting the nonadiabatic dynamics in realistic charge and energy transfer processes in photoinduced energy conversion.

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.

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

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

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.

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.

Simulation of Nonlinear Femtosecond Signals at Finite Temperature via a Thermo Field Dynamics-Tensor Train Method: General Theory and Application to Time- and Frequency-Resolved Fluorescence of the Fenna-Matthews-Olson Complex

Addressing needs of contemporary nonlinear femtosecond optical spectroscopy, we have developed a fully quantum, numerically accurate wave function-based approach for the calculation of third-order spectroscopic signals of polyatomic molecules and molecular aggregates at finite temperature. The systems are described by multimode nonadiabatic vibronic-coupling Hamiltonians, in which diagonal terms are treated in harmonic approximation, while off-diagonal interstate couplings are assumed to be coordinate independent. The approach is based on the Thermo Field Dynamics (TFD) representation of quantum mechanics and tensor-train (TT) machinery for efficient numerical simulation of quantum evolution of systems with many degrees of freedom. The developed TFD-TT approach is applied to the calculation of time- and frequency-resolved fluorescence spectra of the Fenna-Matthews-Olson (FMO) antenna complex at room temperature taking into account finite time-frequency resolution in fluorescence detection, orientational averaging, and static disorder.

Simulation of Time- and Frequency-Resolved Four-Wave-Mixing Signals at Finite Temperatures: A Thermo-Field Dynamics Approach

We propose a new approach to simulate four-wave-mixing signals of molecular systems at finite temperatures by combining the multiconfigurational Ehrenfest method with the thermo-field dynamics theory. In our approach, the four-time correlation functions at finite temperatures are mapped onto those at zero temperature in an enlarged Hilbert space with twice the vibrational degrees of freedom. As an illustration, we have simulated three multidimensional spectroscopic signals, time- and frequency-resolved fluorescence spectra, transient-absorption pump-probe spectra, and electronic two-dimensional (2D) spectra at finite temperatures, for a conical intersection-mediated singlet fission model of a rubrene crystal. It is shown that a detailed dynamical picture of the singlet fission process can be extracted from the three spectroscopic signals. An increasing temperature leads to lower intensities of the signals and broadened vibrational peaks, which can be attributed to faster singlet-triplet population transfer and stronger bath-induced electronic dephasing at higher temperatures.

Evaluating the anharmonicity contributions to the molecular excited state internal conversion rates with finite temperature TD-DMRG

In this work, we propose a new method to calculate molecular nonradiative electronic relaxation rates based on the numerically exact time-dependent density matrix renormalization group theory. This method could go beyond the existing frameworks under the harmonic approximation (HA) of the potential energy surface (PES) so that the anharmonic effect could be considered, which is of vital importance when the electronic energy gap is much larger than the vibrational frequency. We calculate the internal conversion (IC) rates in a two-mode model with Morse potential to investigate the validity of HA. We find that HA is unsatisfactory unless only the lowest several vibrational states of the lower electronic state are involved in the transition process when the adiabatic excitation energy is relatively low. As the excitation energy increases, HA first underestimates and then overestimates the IC rates when the excited state PES shifts toward the dissociative side of the ground state PES. On the contrary, HA slightly overestimates the IC rates when the excited state PES shifts toward the repulsive side. In both cases, a higher temperature enlarges the error of HA. As a real example to demonstrate the effectiveness and scalability of the method, we calculate the IC rates of azulene from S₁ to S₀ on the ab initio anharmonic PES approximated by the one-mode representation. The calculated IC rates of azulene under HA are consistent with the analytically exact results. The rates on the anharmonic PES are 30%-40% higher than the rates under HA.

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.

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.