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.
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.
The generalized quantum master equation (GQME) provides a powerful framework for simulating electronically nonadiabatic molecular dynamics. Within this framework, the effect of the nuclear degrees of freedom on the time evolution of the electronic reduced density matrix is fully captured by a memory kernel superoperator. In this paper, we consider two different procedures for calculating the memory kernel of the GQME from projection-free inputs obtained via the combination of the mapping Hamiltonian (MH) approach and the linearized semiclassical (LSC) approximation. The accuracy and feasibility of the two procedures are demonstrated on the spin-boson model. We find that although simulating the electronic dynamics by direct application of the two LSC-based procedures leads to qualitatively different results that become increasingly less accurate with increasing time, restricting their use to calculating the memory kernel leads to an accurate description of the electronic dynamics. Comparison with a previously proposed procedure for calculating the memory kernel via the Ehrenfest method reveals that MH/LSC methods produce memory kernels that are better behaved at long times and lead to more accurate electronic dynamics.
This thesis describes the development of practical and efficient computational approaches to the quantum dynamics of complex systems. Most of the work presented here relies on the generalized quantum master equation (GQME) formalism, which provides a simple equation of motion of reduced dimensionality for a set of dynamical quantities, e.g., nonequilibrium averages and equilibrium time correlation functions. The reduced dimensionality of the GQME comes at a cost: the introduction of the memory kernel, which accounts for the influence of all ``excluded'' degrees of freedom. Focusing first on the second-order perturbative treatment of the memory kernel known as Redfield theory, I present a collaborative effort to extend its applicability into highly non-Markovian regions via a mode freezing approach. In this method, a portion of bath modes characterized by low frequencies are treated as sources of static disorder and used to calculate modified Redfield dynamics. Application of the method to the spin-boson and FMO complex models indicates that the Redfield+frozen modes scheme consistently produces dynamics that are as good or better than bare Redfield dynamics. Next, we explore GQME approach coupled to the self-consistent solution of the memory kernel, which requires the calculation of auxiliary kernels. Previous implementations of the method had shown impressive boosts in efficiency and, when approximate methods were used to calculate the auxiliary kernels, accuracy over direct calculation of nonequilibrium averages. We show that this method, when formulated from the Mori perspective, is equally applicable to nonequilibrium averages and equilibrium correlation functions. In addition, we examine the dependence of the improvements afforded by the GQME framework on the choice projection operator and kernel closure. In particular, we demonstrate that improvements in efficiency, which rely on short memory lifetimes, are sensitively dependent on the choice of projection operator, and that the choice of kernel closure directly dictates the improvements in accuracy. In addition, we present evidence that indicates that the success of the GQME formalism when the auxiliary kernels are calculated via semi- and quasi-classical methods is largely due to the exact sampling of bath operators at t=0 required by the calculation of specific kernel closures. Next, we provide analytical arguments that delineate when the GQME framework coupled to the self-consistent solution of the memory kernel is likely to provide improvements in efficiency and accuracy. Finally, we present a path integral framework that can efficiently render the partially Wigner-transformed canonical density operator for systems coupled linearly to harmonic baths. This approach permits the direct calculation of any thermodynamic quantity and can be integrated into dynamical schemes like the Ehrenfest, surface hopping, or linearized semi-classical initial value representation methods to calculate equilibrium correlation functions.
We present a modified approach for simulating electronically nonadiabatic dynamics based on the Nakajima-Zwanzig generalized quantum master equation (GQME). The modified approach utilizes the fact that the Nakajima-Zwanzig formalism does not require casting the overall Hamiltonian in system-bath form, which is arguably neither natural nor convenient in the case of the Hamiltonian that governs nonadiabatic dynamics. Within the modified approach, the effect of the nuclear degrees of freedom on the time evolution of the electronic reduced density operator is fully captured by a memory kernel super-operator. A methodology for calculating the memory kernel from projection-free inputs is developed. Simulating the electronic dynamics via the modified approach, with a memory kernel obtained using exact or approximate methods, can be more cost effective and/or lead to more accurate results than direct application of those methods. The modified approach is compared to previously proposed GQME-based approaches, and its robustness and accuracy are demonstrated on a benchmark spin-boson model with a memory kernel which is calculated within the Ehrenfest method.
The generalized quantum master equation (GQME) approach provides a powerful general-purpose framework for simulating the inherently quantum mechanical dynamics of a subset of electronic reduced density matrix elements of interest in complex molecular systems. Previous studies have found that combining the GQME approach with quasiclassical mapping Hamiltonian (QC/MH) methods can dramatically improve the accuracy of electronic populations obtained via those methods. In this paper, we perform a complimentary study of the advantages offered by the GQME approach for simulating the dynamics of electronic coherences, which play a central role in optical spectroscopy, quantum information science, and quantum technology. To this end, we focus on cases where the electronic coherences predicted for the spin-boson benchmark model by direct application of various QC/MH methods are inaccurate. We find that similar to the case of electronic populations, combining the QC/MH methods with the GQME approach can dramatically improve the accuracy of the electronic coherences obtained via those methods. We also provide a comprehensive analysis of how the performance of GQMEs depends on the choice of projection operator and electronic basis and show that the accuracy and feasibility of the GQME approach can benefit from casting the GQME in terms of the eigen-basis of the observable of interest.
The generalized quantum master equation (GQME) approach provides a rigorous framework for deriving the exact equation of motion for any subset of electronic reduced density matrix elements (e.g., the diagonal elements). In the context of electronic dynamics, the memory kernel and inhomogeneous term of the GQME introduce the implicit coupling to nuclear motion and dynamics of electronic density matrix elements that are projected out (e.g., the off-diagonal elements), allowing for efficient quantum dynamics simulations. Here, we focus on benchmark quantum simulations of electronic dynamics in a spin-boson model system described by various types of GQMEs. Exact memory kernels and inhomogeneous terms are obtained from short-time quantum-mechanically exact tensor-train thermo-field dynamics (TT-TFD) simulations and are compared with those obtained from an approximate linearized semiclassical method, allowing for assessment of the accuracy of these approximate memory kernels and inhomogeneous terms. Moreover, we have analyzed the computational cost of the full and reduced-dimensionality GQMEs. The scaling of the computational cost is dependent on several factors, sometimes with opposite scaling trends. The TT-TFD memory kernels can provide insights on the main sources of inaccuracies of GQME approaches when combined with approximate input methods and pave the road for the development of quantum circuits that implement GQMEs on digital quantum computers.
The Nakajima-Zwanzig generalized quantum master equation (GQME) provides a general, and formally exact, prescription for simulating the reduced dynamics of a quantum system coupled to a, possibly anharmonic, quantum bath. In this equation, a memory kernel superoperator accounts for the influence of the bath on the dynamics of the system. In a previous paper [Q. Shi and E. Geva, J. Chem. Phys. 119, 12045 (2003)] we proposed a new approach to calculating the memory kernel, in the case of arbitrary system-bath coupling. Within this approach, the memory kernel is obtained by solving a set of two integral equations, which requires a new type of two-time system-dependent bath correlation functions as input. In the present paper, we consider the application of the linearized semiclassical (LSC) approximation for calculating those correlation functions, and subsequently the memory kernel. The new approach is tested on a benchmark spin-boson model. Application of the LSC approximation for calculating the relatively short-lived memory kernel, followed by a numerically exact solution of the GQME, is found to provide an accurate description of the relaxation dynamics. The success of the proposed LSC-GQME methodology is contrasted with the failure of both the direct application of the LSC approximation and the weak coupling treatment to provide an accurate description of the dynamics, for the same model, except at very short times. The feasibility of the new methodology to anharmonic systems is also demonstrated in the case of a two level system coupled to a chain of Lennard-Jones atoms.