Generalized quantum master equations in and out of equilibrium: When can one win?

Unified framework for open quantum dynamics with memory

The dynamics of quantum systems coupled to baths are typically studied using the Nakajima-Zwanzig memory kernel ( K ) or the influence functions (I), particularly when memory effects are present. Despite their significance, formal connections between the two have not been explicitly known. We establish their connections by examining the system propagator for a N-level system linearly coupled to Gaussian baths with various types of system-bath coupling. For a certain class of problems, we devised a non-perturbative, diagrammatic approach to construct K from I for (driven) systems interacting with Gaussian baths, bypassing conventional projection-free dynamics inputs. Our work provides a way to interpret approximate path integral methods in terms of approximate memory kernels. Moreover, it offers a Hamiltonian learning procedure to extract the bath spectral density from reduced system trajectories, opening new avenues in quantum sensing and engineering. The insights we provide advance our understanding of non-Markovian dynamics and will serve as a stepping stone for future theoretical and experimental developments in this area.

Semiclassical dynamics in Wigner phase space II: Nonadiabatic hybrid Wigner dynamics

We present an approximate semiclassical (SC) framework for mixed quantized dynamics in Wigner phase space in a two-part series. In the first article, we introduced the Adiabatic Hybrid Wigner Dynamics (AHWD) method that allows for a few important "system" degrees of freedom to be quantized using high-level double Herman-Kluk SC theory while describing the rest (the "bath") using classical-limit linearized SC theory. In this second article, we extend our hybrid Wigner dynamics to nonadiabatic processes. The resulting Nonadiabatic Hybrid Wigner Dynamics (NHWD) has two variants that differ in the choice of degrees of freedom to be quantized. Specifically, we introduce NHWD(E) where only the electronic state variables are quantized and the NHWD(V) where both electronic state variables and a handful of strongly coupled nuclear modes are quantized. We show that while NHWD(E) proves accurate for a wide range of scattering models and spin-boson models, systems where a few nuclear modes are strongly coupled to electronic states require NHWD(V) to accurately capture the long-time dynamics. Taken together, we show that AHWD and NHWD represent a new framework for SC simulations of high-dimensional systems with significant quantum effects.

Simulating Chemistry on Bosonic Quantum Devices

Bosonic quantum devices offer a novel approach to realize quantum computations, where the quantum two-level system (qubit) is replaced with the quantum (an)harmonic oscillator (qumode) as the fundamental building block of the quantum simulator. The simulation of chemical structure and dynamics can then be achieved by representing or mapping the system Hamiltonians in terms of bosonic operators. In this Perspective, we review recent progress and future potential of using bosonic quantum devices for addressing a wide range of challenging chemical problems, including the calculation of molecular vibronic spectra, the simulation of gas-phase and solution-phase adiabatic and nonadiabatic chemical dynamics, the efficient solution of molecular graph theory problems, and the calculations of electronic structure.

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.

Bath-induced interactions and transient dynamics in open quantum systems at strong coupling: Effective Hamiltonian approach

Understanding the dynamics of dissipative quantum systems, particularly beyond the weak coupling approximation, is central to various quantum applications. While numerically exact methods provide accurate solutions, they often lack the analytical insight provided by theoretical approaches. In this study, we employ the recently developed method dubbed the effective Hamiltonian theory to understand the dynamics of system-bath configurations without resorting to a perturbative description of the system-bath coupling energy. Through a combination of mapping steps and truncation, the effective Hamiltonian theory offers both analytical insights into signatures of strong couplings in open quantum systems and a straightforward path for numerical simulations. To validate the accuracy of the method, we apply it to two canonical models: a single spin immersed in a bosonic bath and two noninteracting spins in a common bath. In both cases, we study the transient regime and the steady state limit at nonzero temperature and spanning system-bath interactions from the weak to the strong regime. By comparing the results of the effective Hamiltonian theory with numerically exact simulations, we show that although the former overlooks non-Markovian features in the transient equilibration dynamics, it correctly captures non-perturbative bath-generated couplings between otherwise non-interacting spins, as observed in their synchronization dynamics and correlations. Altogether, the effective Hamiltonian theory offers a powerful approach for understanding strong coupling dynamics and thermodynamics, capturing the signatures of such interactions in both relaxation dynamics and in the steady state limit.

General framework for quantifying dissipation pathways in open quantum systems. I. Theoretical formulation

We present a general and practical theoretical framework to investigate how energy is dissipated in open quantum system dynamics. This is performed by quantifying the contributions of individual bath components to the overall dissipation of the system. The framework is based on the Nakajima-Zwanzig projection operator technique, which allows us to express the rate of energy dissipation into a specific bath degree of freedom by using traces of operator products. The approach captures system-bath interactions to all orders, but is based on second-order perturbation theory on the off-diagonal subsystem's couplings and a Markovian description of the bath. The usefulness of our theory is demonstrated by applying it to various models of open quantum systems involving harmonic oscillators or spin baths and connecting the outcomes to existing results such as our previously reported formula derived for locally coupled harmonic baths [Kim and Franco, J. Chem. Phys. 154, 084109 (2021)]. We also prove that the dissipation calculated by our theory rigorously satisfies thermodynamic principles such as energy conservation and detailed balance. Overall, the strategy can be used to develop the theory and simulation of dissipation pathways to interpret and engineer the dynamics of open quantum systems.

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.

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.

Predicting rate kernels via dynamic mode decomposition

Simulating dynamics of open quantum systems is sometimes a significant challenge, despite the availability of various exact or approximate methods. Particularly when dealing with complex systems, the huge computational cost will largely limit the applicability of these methods. In this work, we investigate the usage of dynamic mode decomposition (DMD) to evaluate the rate kernels in quantum rate processes. DMD is a data-driven model reduction technique that characterizes the rate kernels using snapshots collected from a small time window, allowing us to predict the long-term behaviors with only a limited number of samples. Our investigations show that whether the external field is involved or not, the DMD can give accurate prediction of the result compared with the traditional propagations, and simultaneously reduce the required computational cost.

Integrative generalized master equation: A method to study long-timescale biomolecular dynamics via the integrals of memory kernels

The generalized master equation (GME) provides a powerful approach to study biomolecular dynamics via non-Markovian dynamic models built from molecular dynamics (MD) simulations. Previously, we have implemented the GME, namely the quasi Markov State Model (qMSM), where we explicitly calculate the memory kernel and propagate dynamics using a discretized GME. qMSM can be constructed with much shorter MD trajectories than the MSM. However, since qMSM needs to explicitly compute the time-dependent memory kernels, it is heavily affected by the numerical fluctuations of simulation data when applied to study biomolecular conformational changes. This can lead to numerical instability of predicted long-time dynamics, greatly limiting the applicability of qMSM in complicated biomolecules. We present a new method, the Integrative GME (IGME), in which we analytically solve the GME under the condition when the memory kernels have decayed to zero. Our IGME overcomes the challenges of the qMSM by using the time integrations of memory kernels, thereby avoiding the numerical instability caused by explicit computation of time-dependent memory kernels. Using our solutions of the GME, we have developed a new approach to compute long-time dynamics based on MD simulations in a numerically stable, accurate and efficient way. To demonstrate its effectiveness, we have applied the IGME in three biomolecules: the alanine dipeptide, FIP35 WW-domain, and Taq RNA polymerase. In each system, the IGME achieves significantly smaller fluctuations for both memory kernels and long-time dynamics compared to the qMSM. We anticipate that the IGME can be widely applied to investigate biomolecular conformational changes.

A deep learning approach to the measurement of long-lived memory kernels from generalized Langevin dynamics

Memory effects are ubiquitous in a wide variety of complex physical phenomena, ranging from glassy dynamics and metamaterials to climate models. The Generalized Langevin Equation (GLE) provides a rigorous way to describe memory effects via the so-called memory kernel in an integro-differential equation. However, the memory kernel is often unknown, and accurately predicting or measuring it via, e.g., a numerical inverse Laplace transform remains a herculean task. Here, we describe a novel method using deep neural networks (DNNs) to measure memory kernels from dynamical data. As a proof-of-principle, we focus on the notoriously long-lived memory effects of glass-forming systems, which have proved a major challenge to existing methods. In particular, we learn the operator mapping dynamics to memory kernels from a training set generated with the Mode-Coupling Theory (MCT) of hard spheres. Our DNNs are remarkably robust against noise, in contrast to conventional techniques. Furthermore, we demonstrate that a network trained on data generated from analytic theory (hard-sphere MCT) generalizes well to data from simulations of a different system (Brownian Weeks-Chandler-Andersen particles). Finally, we train a network on a set of phenomenological kernels and demonstrate its effectiveness in generalizing to both unseen phenomenological examples and supercooled hard-sphere MCT data. We provide a general pipeline, KernelLearner, for training networks to extract memory kernels from any non-Markovian system described by a GLE. The success of our DNN method applied to noisy glassy systems suggests that deep learning can play an important role in the study of dynamical systems with memory.

Memory Unlocks the Future of Biomolecular Dynamics: Transformative Tools to Uncover Physical Insights Accurately and Efficiently

Conformational changes underpin function and encode complex biomolecular mechanisms. Gaining atomic-level detail of how such changes occur has the potential to reveal these mechanisms and is of critical importance in identifying drug targets, facilitating rational drug design, and enabling bioengineering applications. While the past two decades have brought Markov state model techniques to the point where practitioners can regularly use them to glimpse the long-time dynamics of slow conformations in complex systems, many systems are still beyond their reach. In this Perspective, we discuss how including memory (i.e., non-Markovian effects) can reduce the computational cost to predict the long-time dynamics in these complex systems by orders of magnitude and with greater accuracy and resolution than state-of-the-art Markov state models. We illustrate how memory lies at the heart of successful and promising techniques, ranging from the Fokker-Planck and generalized Langevin equations to deep-learning recurrent neural networks and generalized master equations. We delineate how these techniques work, identify insights that they can offer in biomolecular systems, and discuss their advantages and disadvantages in practical settings. We show how generalized master equations can enable the investigation of, for example, the gate-opening process in RNA polymerase II and demonstrate how our recent advances tame the deleterious influence of statistical underconvergence of the molecular dynamics simulations used to parameterize these techniques. This represents a significant leap forward that will enable our memory-based techniques to interrogate systems that are currently beyond the reach of even the best Markov state models. We conclude by discussing some current challenges and future prospects for how exploiting memory will open the door to many exciting opportunities.

A derivation of the conditions under which bosonic operators exactly capture fermionic structure and dynamics

The dynamics of many-body fermionic systems are important in problems ranging from catalytic reactions at electrochemical surfaces to transport through nanojunctions and offer a prime target for quantum computing applications. Here, we derive the set of conditions under which fermionic operators can be exactly replaced by bosonic operators that render the problem amenable to a large toolbox of dynamical methods while still capturing the correct dynamics of n-body operators. Importantly, our analysis offers a simple guide on how one can exploit these simple maps to calculate nonequilibrium and equilibrium single- and multi-time correlation functions essential in describing transport and spectroscopy. We use this to rigorously analyze and delineate the applicability of simple yet effective Cartesian maps that have been shown to correctly capture the correct fermionic dynamics in select models of nanoscopic transport. We illustrate our analytical results with exact simulations of the resonant level model. Our work provides new insights as to when one can leverage the simplicity of bosonic maps to simulate the dynamics of many-electron systems, especially those where an atomistic representation of nuclear interactions becomes essential.

Tensor-Train Thermo-Field Memory Kernels for Generalized Quantum Master Equations

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.

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.

Compact and complete description of non-Markovian dynamics

Generalized master equations provide a theoretically rigorous framework to capture the dynamics of processes ranging from energy harvesting in plants and photovoltaic devices to qubit decoherence in quantum technologies and even protein folding. At their center is the concept of memory. The explicit time-nonlocal description of memory is both protracted and elaborate. When physical intuition is at a premium, one would desire a more compact, yet complete, description. Here, we demonstrate how and when the time-convolutionless formalism constitutes such a description. In particular, by focusing on the dissipative dynamics of the spin-boson and Frenkel exciton models, we show how to: easily construct the time-local generator from reference reduced dynamics, elucidate the dependence of its existence on the system parameters and the choice of reduced observables, identify the physical origin of its apparent divergences, and offer analysis tools to diagnose their severity and circumvent their deleterious effects. We demonstrate that, when applicable, the time-local approach requires as little information as the more commonly used time-nonlocal scheme, with the important advantages of providing a more compact description, greater algorithmic simplicity, and physical interpretability. We conclude by introducing the discrete-time analog and a straightforward protocol to employ it in cases where the reference dynamics have limited resolution. The insights we present here offer the potential for extending the reach of dynamical methods, reducing both their cost and conceptual complexity.

Quasiclassical approaches to the generalized quantum master equation

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

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

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

Generalized master equation for charge transport in a molecular junction: Exact memory kernels and their high order expansion

We derive a set of generalized master equations (GMEs) to study charge transport dynamics in molecular junctions using the Nakajima-Zwanzig-Mori projection operator approach. In the new GME, time derivatives of population on each quantum state of the molecule, as well as the tunneling current, are calculated as the convolution of time non-local memory kernels with populations on all system states. The non-Markovian memory kernels are obtained by combining the hierarchical equations of motion (HEOM) method and a previous derived Dyson relation for the exact kernel. A perturbative expansion of these memory kernels is then calculated using the extended HEOM developed in our previous work [M. Xu et al., J. Chem. Phys. 146, 064102 (2017)]. By using the resonant level model and the Anderson impurity model, we study properties of the exact memory kernels and analyze convergence properties of their perturbative expansions with respect to the system-bath coupling strength and the electron-electron repulsive energy. It is found that exact memory kernels calculated from HEOM exhibit short memory times and decay faster than the population and current dynamics. The high order perturbation expansion of the memory kernels can give converged results in certain parameter regimes. The Padé and Landau-Zener resummation schemes are also found to give improved results over low order perturbation theory.

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.

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.

A Road Map to Various Pathways for Calculating the Memory Kernel of the Generalized Quantum Master Equation

The generalized quantum master equation (GQME) provides a powerful framework for simulating electronic energy, charge, and coherence transfer dynamics in molecular systems. 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. However, the actual memory kernel depends on the choice of projection operator and is therefore not unique. Furthermore, calculating the memory kernel can be done in multiple ways that use different forms of projection-free inputs. Although the electronic dynamics is invariant to those choices when quantum-mechanically exact projection-free inputs are used, this is not the case when they are obtained via more feasible semiclassical or mixed quantum-classical approximate methods. Furthermore, the accuracy and numerical stability of the resulting electronic dynamics has been observed to be sensitive to the above-mentioned choices when approximate methods are used to calculate the projection-free inputs. In this article, we provide a systematic road map to 30 possible pathways for calculating the memory kernel and highlight how they are related as well as the ways in which they differ. We also compare the performance of different pathways in the context of the spin-boson benchmark model, with the projection-free inputs obtained via a mapping Hamiltonian linearized semiclassical method. In this case, we find that expressing the memory kernel with an exponential operator where the projection operator precedes the Liouvillian yields the most accurate and most numerically stable results.

Constructing tensor network influence functionals for general quantum dynamics

We describe an iterative formalism to compute influence functionals that describe the general quantum dynamics of a subsystem beyond the assumption of linear coupling to a quadratic bath. We use a space-time tensor network representation of the influence functional and investigate its approximability in terms of its bond dimension and time-like entanglement in the tensor network description. We study two numerical models, the spin-boson model and a model of interacting hard-core bosons in a 1D harmonic trap. We find that the influence functional and the intermediates involved in its construction can be efficiently approximated by low bond dimension tensor networks in certain dynamical regimes, which allows the quantum dynamics to be accurately computed for longer times than with direct time evolution methods. However, as one iteratively integrates out the bath, the correlations in the influence functional can first increase before decreasing, indicating that the final compressibility of the influence functional is achieved via non-trivial cancellation.

Simulating energy transfer dynamics in the Fenna-Matthews-Olson complex via the modified generalized quantum master equation

The generalized quantum master equation (GQME) provides a general and formally exact framework for simulating the reduced dynamics of open quantum systems. The recently introduced modified approach to the GQME (M-GQME) corresponds to a specific implementation of the GQME that is geared toward simulating the dynamics of the electronic reduced density matrix in systems governed by an excitonic Hamiltonian. Such a Hamiltonian, which is often used for describing energy and charge transfer dynamics in complex molecular systems, is given in terms of diabatic electronic states that are coupled to each other and correspond to different nuclear Hamiltonians. Within the M-GQME approach, the effect of the nuclear degrees of freedom on the time evolution of the electronic density matrix is fully captured by a memory kernel superoperator, which can be obtained from short-lived (compared to the time scale of energy/charge transfer) projection-free inputs. In this paper, we test the ability of the M-GQME to predict the energy transfer dynamics within a seven-state benchmark model of the Fenna-Matthews-Olson (FMO) complex, with the short-lived projection-free inputs obtained via the Ehrenfest method. The M-GQME with Ehrenfest-based inputs is shown to yield accurate results across a wide parameter range. It is also found to dramatically outperform the direct application of the Ehrenfest method and to provide better-behaved convergence with respect to memory time in comparison to an alternative implementation of the GQME approach previously applied to the same FMO model.

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

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

Quantum Dynamics of Exciton Transport and Dissociation in Multichromophoric Systems

Due to the subtle interplay of site-to-site electronic couplings, exciton delocalization, nonadiabatic effects, and vibronic couplings, quantum dynamical studies are needed to elucidate the details of ultrafast photoinduced energy and charge transfer events in organic multichromophoric systems. In this vein, we review an approach that combines first-principles parameterized lattice Hamiltonians with accurate quantum dynamical simulations using advanced multiconfigurational methods. Focusing on the elementary transfer steps in organic functional materials, we address coherent exciton migration and creation of charge transfer excitons in homopolymers, notably representative of the poly(3-hexylthiophene) material, as well as exciton dissociation at polymer:fullerene heterojunctions. We emphasize the role of coherent transfer, trapping effects due to high-frequency phonon modes, and thermal activation due to low-frequency soft modes that drive a diffusive dynamics.

On the advantages of exploiting memory in Markov state models for biomolecular dynamics

Biomolecular dynamics play an important role in numerous biological processes. Markov State Models (MSMs) provide a powerful approach to study these dynamic processes by predicting long time scale dynamics based on many short molecular dynamics (MD) simulations. In an MSM, protein dynamics are modeled as a kinetic process consisting of a series of Markovian transitions between different conformational states at discrete time intervals (called "lag time"). To achieve this, a master equation must be constructed with a sufficiently long lag time to allow interstate transitions to become truly Markovian. This imposes a major challenge for MSM studies of proteins since the lag time is bound by the length of relatively short MD simulations available to estimate the frequency of transitions. Here, we show how one can employ the generalized master equation formalism to obtain an exact description of protein conformational dynamics both at short and long time scales without the time resolution restrictions imposed by the MSM lag time. Using a simple kinetic model, alanine dipeptide, and WW domain, we demonstrate that it is possible to construct these quasi-Markov State Models (qMSMs) using MD simulations that are 5-10 times shorter than those required by MSMs. These qMSMs only contain a handful of metastable states and, thus, can greatly facilitate the interpretation of mechanisms associated with protein dynamics. A qMSM opens the door to the study of conformational changes of complex biomolecules where a Markovian model with a few states is often difficult to construct due to the limited length of available MD simulations.

Improved population operators for multi-state nonadiabatic dynamics with the mixed quantum-classical mapping approach

Application to the 7-state Frenkel-exciton Hamiltonian for the Fenna–Matthews–Olson complex shows that using a different representation of the electronic population operators can drastically improve the accuracy of the quasiclassical mapping approach without increasing the computational effort.

Exciton dissociation and charge separation at donor–acceptor interfaces from quantum-classical dynamics simulations

Nonadiabatic dynamics simulations based on the quantum-classical Liouville equation are employed to study the real-time dynamics of exciton dissociation and charge separation at a model donor–acceptor interface.

Benchmarking semiclassical and perturbative methods for real-time simulations of cavity-bound emission and interference

We benchmark a selection of semiclassical and perturbative dynamics techniques by investigating the correlated evolution of a cavity-bound atomic system to assess their applicability to study problems involving strong light-matter interactions in quantum cavities. The model system of interest features spontaneous emission, interference, and strong coupling behavior and necessitates the consideration of vacuum fluctuations and correlated light-matter dynamics. We compare a selection of approximate dynamics approaches including fewest switches surface hopping (FSSH), multitrajectory Ehrenfest dynamics, linearized semiclassical dynamics, and partially linearized semiclassical dynamics. Furthermore, investigating self-consistent perturbative methods, we apply the Bogoliubov-Born-Green-Kirkwood-Yvon hierarchy in the second Born approximation. With the exception of fewest switches surface hopping, all methods provide a reasonable level of accuracy for the correlated light-matter dynamics, with most methods lacking the capacity to fully capture interference effects.

A hybrid memory kernel approach for condensed phase non-adiabatic dynamics

The spin-boson model is a simplified Hamiltonian often used to study non-adiabatic dynamics in large condensed phase systems, even though it has not been solved in a fully analytic fashion. Herein, we present an exact analytic expression for the dynamics of the spin-boson model in the infinitely slow-bath limit and generalize it to approximate dynamics for faster baths. We achieve the latter by developing a hybrid approach that combines the exact slow-bath result with the popular non-interacting blip approximation (NIBA) method to generate a memory kernel that is formally exact to second-order in the diabatic coupling but also contains higher-order contributions approximated from the second-order term alone. This kernel has the same computational complexity as the NIBA, but is found to yield dramatically superior dynamics in regimes where the NIBA breaks down-such as systems with large diabatic coupling or energy bias. This indicates that this hybrid approach could be used to cheaply incorporate higher-order effects into second-order methods and could potentially be generalized to develop alternate kernel resummation schemes.