Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/0612203 · The Journal of Chemical Physics

Decoherence and quantum-classical master equation dynamics

Robbie Grunwald, Raymond Kapral

Published Mar 20, 2007• 23 citations
arXivHTMLPDF

Abstract

The conditions under which quantum-classical Liouville dynamics may be reduced to a master equation are investigated. Systems that can be partitioned into a quantum-classical subsystem interacting with a classical bath are considered. Starting with an exact non-Markovian equation for the diagonal elements of the density matrix, an evolution equation for the subsystem density matrix is derived. One contribution to this equation contains the bath average of a memory kernel that accounts for all coherences in the system. It is shown to be a rapidly decaying function, motivating a Markovian approximation on this term in the evolution equation. The resulting subsystem density matrix equation is still non-Markovian due to the fact that bath degrees of freedom have been projected out of the dynamics. Provided the computation of nonequilibrium average values or correlation functions is considered, the non-Markovian character of this equation can be removed by lifting the equation into the full phase space of the system. This leads to a trajectory description of the dynamics where each fictitious trajectory accounts for decoherence due to the bath degrees of freedom. The results are illustrated by computations of the rate constant of a model nonadiabatic chemical reaction.

DOI
10.1063/1.2567164

Related papers

Quantum-Classical Master Equation Dynamics: An Analysis of Decoherence and Surface-hopping Techniques

Robbie Grunwald

In this thesis quantum-classical dynamics is applied to the study of quantum condensed phase processes. This approach is based on the quantum-classical Liouville equation where the dynamics of a small subset of the degrees of freedom are treated quantum mechanically while the remaining degrees of freedom are treated by classical mechanics to a good approximation. We use this approach as it is computationally tractable, and the resulting equation of motion accurately accounts for the quantum and classical dynamics, as well as the coupling between these two components of the system. \n\nBy recasting the quantum-classical Liouville equation into the form of a generalized master equation we investigate connections to surface-hopping. The link between these approaches is decoherence arising from interaction of the subsystem with the environment. We derive an evolution equation for the subsystem which contains terms accounting for the effects of the environment. One of these terms involves a memory kernel that accounts for the coherent dynamics. If this term decays rapidly, a Markovian approximation can be made. By lifting the resulting subsystem master equation into the full phase space, we obtain a Markovian master equation that prescribes surface-hopping-like dynamics. Our analysis outlines the conditions under which such a description is valid.\n\nNext, we consider the calculation of the rate constant for a quantum mechanical barrier crossing process. Starting from the reactive-flux autocorrelation function, we derive a quantum-classical expression for the rate kernel. This expression involves quantum-classical evolution of a species operator averaged over the initial quantum equilibrium structure of the system making it possible to compute the rate constant via computer simulation. Using a simple model for a proton transfer reaction we compare the results of the rate calculation obtained by quantum-classical Liouville dynamics with that of master equation dynamics. The master equation provides a good approximation to the full quantum-classical Liouville calculation for our model and a more stable algorithm results due to the elimination of oscillating phase factors in the simulation. \n\nFinally, we make use of the theoretical framework established in this thesis to analyze some aspects of decoherence used in popular surface-hopping techniques.

•Jan 19, 2009•0 citations

Decoherence and quantum-classical dynamics in a dissipative bath

Jean Pierre Rank, Raymond Kapral

The dynamics of a mixed quantum-classical system, in which the classical subsystem interacts with a dissipative bath, is investigated. This description of the dynamics will be appropriate if the details of the bath dynamics are unimportant but its presence plays an important role in the dissipation of energy to the environment. In this dynamical description, which can be simulated employing an ensemble of stochastic surface-hopping trajectories, the strength of the dissipation is controlled by a friction coefficient. We show that if decoherence, whose effects are controlled by the bath friction, is sufficiently rapid, the equation of motion can be reduced to a master equation. Thus, decoherence and the validity of master equation models may be explored as a function of bath friction. We use this framework to study the mechanism of decoherence in a simple model nonadiabatic chemical reaction.

The Journal of Chemical Physics•Feb 17, 2010•14 citations

Surface-hopping dynamics and decoherence with quantum equilibrium structure

Robbie Grunwald, Hyojoon Kim, Raymond Kapral

In open quantum systems, decoherence occurs through interaction of a quantum subsystem with its environment. The computation of expectation values requires a knowledge of the quantum dynamics of operators and sampling from initial states of the density matrix describing the subsystem and bath. We consider situations where the quantum evolution can be approximated by quantum-classical Liouville dynamics and examine the circumstances under which the evolution can be reduced to surface-hopping dynamics, where the evolution consists of trajectory segments exclusively evolving on single adiabatic surfaces, with probabilistic hops between these surfaces. The justification for the reduction depends on the validity of a Markovian approximation on a bath averaged memory kernel that accounts for quantum coherence in the system. We show that such a reduction is often possible when initial sampling is from either the quantum or classical bath initial distributions. If the average is taken only over the quantum dispersion that broadens the classical distribution, then such a reduction is not always possible.

arXiv:0801.3406•The Journal of Chemical Physics•Apr 24, 2008•18 citations

A Study of Quantum-classical Dynamics in the Mapping Basis

Ali M. Nassimi

Solving quantum dynamics is an exponentially difficult problem. Thus, an exact numerical solution is inaccessible for any condensed matter system.\nA promising approach is to divide the system into a quantum subsystem containing degrees of freedom which are of greater interest or those which have more profound quantum character (e.g., have smaller mass) and a classical bath containing the rest of the system.\nImposing such a partition and treating the bath classically results in quantum-classical dynamics. The quantum-classical Liouville equation is a general equation in the Hilbert space of quantum degrees of freedom while it resides in the phase space of the classical degrees of freedom.\nAny numerical solution to this equation requires representation of the quantum subsystem in some basis. Solutions to this equation have been already proposed in the subsystem, adiabatic and force bases, each with its own cons and pros.\n\nIn this work, the quantum-classical equations of motion are cast in the subsystem basis and subsequently mapped to a number of fictitious harmonic oscillators.\nThe result is quantum-classical dynamics in the mapping basis which treats both quantum and classical degrees of freedom on the same footing, i.e., in phase space. Neglecting a portion of the back reaction of the quantum-subsystem to classical bath results in an expression for the time evolution of an operator (density matrix) equal to its Poisson bracket with the Hamiltonian.\nThis equation can be solved in terms of characteristics to provide a computationally tractable method for calculating quantum-classical dynamical properties. The expressions for expectation values and correlation functions in this formalism are derived.\nCalculations on spin-boson system, barrier crossing models---the so called Tully models---and the Fenna-Mathews-Olson pigments show very good agreement between the results of this method and numerical solutions to the Schrödinger equation.

•Aug 31, 2011•0 citations

Decoherence in open quantum systems: influence of the intrinsic bath dynamics

V. V. Ignatyuk, V. G. Morozov

The non-Markovian master equation for open quantum systems is obtained by generalization of the standard Zwanzig-Nakajima (ZN) projection technique. To this end, a coupled chain of equations is written for the reduced density matrices of the bath ϱB(t) and of the system ϱS(t). A formal solution of the equation for ϱB(t) in the 2nd approximation in interaction yields a specific extra term related to the intrinsic bath dynamics. This term is nonlinear in the reduced density matrix ϱS(t), and vanishes in the Markovian limit. To verify the consistence and robustness of our approach, we apply the generalized ZN projection scheme to a simple dephasing model. We study the obtained kinetic equation both in the Markovian approximation and beyond it (for the term related to the intrinsic bath dynamics) and compare the results with the exact ones.

arXiv:2112.05595•Condensed Matter Physics•Jan 1, 2022•4 citations

Master equation for open quantum systems: Zwanzig-Nakajima projection technique and the intrinsic bath dynamics

V. V. Ignatyuk, V. G. Morozov

The non-Markovian master equation for open quantum systems is obtained by generalization of the ordinary Zwanzig-Nakajima (ZN) projection technique. To this end, a coupled chain of equations for the reduced density matrices of the bath $\varrho_{B}(t)$ and the system $\varrho_{S}(t)$ are written down. A formal solution of the equation for $\varrho_{B}(t)$, having been inserted in the equation for the reduced density matrix of the system, in the 2-nd approximation in interaction yields a very specific extra term in the generalized master equation. This term, being nonlinear in $\varrho_{S}(t)$, is related to the intrinsic bath dynamics and vanishes in the Markovian limit. To verify the consistence and robustness of our approach, we applied the generalized ZN projection scheme to a simple dephasing model. It is shown that consideration of the lowest order in interaction is insufficient to describe time evolution of the system coherence adequately. We explain this fact by analyzing the exact and approximate forms of $\varrho_{B}(t)$ and give some hints how to take the dynamic correlations (which originate from the spin-bath coupling) into account.

arXiv:1911.12312•Nov 27, 2019•0 citations
The conditions under which quantum-classical Liouville dynamics may be reduced to a master equation are investigated. Systems that can be partitioned into a quantum-classical subsystem interacting with a classical bath are considered. Starting with an exact non-Markovian equation for the diagonal elements of the density matrix, an evolution equation for the subsystem density matrix is derived. One contribution to this equation contains the bath average of a memory kernel that accounts for all coherences in the system. It is shown to be a rapidly decaying function, motivating a Markovian approximation on this term in the evolution equation. The resulting subsystem density matrix equation is still non-Markovian due to the fact that bath degrees of freedom have been projected out of the dynamics. Provided the computation of nonequilibrium average values or correlation functions is considered, the non-Markovian character of this equation can be removed by lifting the equation into the full phase space of the system. This leads to a trajectory description of the dynamics where each fictitious trajectory accounts for decoherence due to the bath degrees of freedom. The results are illustrated by computations of the rate constant of a model nonadiabatic chemical reaction.
In this thesis quantum-classical dynamics is applied to the study of quantum condensed phase processes. This approach is based on the quantum-classical Liouville equation where the dynamics of a small subset of the degrees of freedom are treated quantum mechanically while the remaining degrees of freedom are treated by classical mechanics to a good approximation. We use this approach as it is computationally tractable, and the resulting equation of motion accurately accounts for the quantum and classical dynamics, as well as the coupling between these two components of the system. \n\nBy recasting the quantum-classical Liouville equation into the form of a generalized master equation we investigate connections to surface-hopping. The link between these approaches is decoherence arising from interaction of the subsystem with the environment. We derive an evolution equation for the subsystem which contains terms accounting for the effects of the environment. One of these terms involves a memory kernel that accounts for the coherent dynamics. If this term decays rapidly, a Markovian approximation can be made. By lifting the resulting subsystem master equation into the full phase space, we obtain a Markovian master equation that prescribes surface-hopping-like dynamics. Our analysis outlines the conditions under which such a description is valid.\n\nNext, we consider the calculation of the rate constant for a quantum mechanical barrier crossing process. Starting from the reactive-flux autocorrelation function, we derive a quantum-classical expression for the rate kernel. This expression involves quantum-classical evolution of a species operator averaged over the initial quantum equilibrium structure of the system making it possible to compute the rate constant via computer simulation. Using a simple model for a proton transfer reaction we compare the results of the rate calculation obtained by quantum-classical Liouville dynamics with that of master equation dynamics. The master equation provides a good approximation to the full quantum-classical Liouville calculation for our model and a more stable algorithm results due to the elimination of oscillating phase factors in the simulation. \n\nFinally, we make use of the theoretical framework established in this thesis to analyze some aspects of decoherence used in popular surface-hopping techniques.
The dynamics of a mixed quantum-classical system, in which the classical subsystem interacts with a dissipative bath, is investigated. This description of the dynamics will be appropriate if the details of the bath dynamics are unimportant but its presence plays an important role in the dissipation of energy to the environment. In this dynamical description, which can be simulated employing an ensemble of stochastic surface-hopping trajectories, the strength of the dissipation is controlled by a friction coefficient. We show that if decoherence, whose effects are controlled by the bath friction, is sufficiently rapid, the equation of motion can be reduced to a master equation. Thus, decoherence and the validity of master equation models may be explored as a function of bath friction. We use this framework to study the mechanism of decoherence in a simple model nonadiabatic chemical reaction.
In open quantum systems, decoherence occurs through interaction of a quantum subsystem with its environment. The computation of expectation values requires a knowledge of the quantum dynamics of operators and sampling from initial states of the density matrix describing the subsystem and bath. We consider situations where the quantum evolution can be approximated by quantum-classical Liouville dynamics and examine the circumstances under which the evolution can be reduced to surface-hopping dynamics, where the evolution consists of trajectory segments exclusively evolving on single adiabatic surfaces, with probabilistic hops between these surfaces. The justification for the reduction depends on the validity of a Markovian approximation on a bath averaged memory kernel that accounts for quantum coherence in the system. We show that such a reduction is often possible when initial sampling is from either the quantum or classical bath initial distributions. If the average is taken only over the quantum dispersion that broadens the classical distribution, then such a reduction is not always possible.
Solving quantum dynamics is an exponentially difficult problem. Thus, an exact numerical solution is inaccessible for any condensed matter system.\nA promising approach is to divide the system into a quantum subsystem containing degrees of freedom which are of greater interest or those which have more profound quantum character (e.g., have smaller mass) and a classical bath containing the rest of the system.\nImposing such a partition and treating the bath classically results in quantum-classical dynamics. The quantum-classical Liouville equation is a general equation in the Hilbert space of quantum degrees of freedom while it resides in the phase space of the classical degrees of freedom.\nAny numerical solution to this equation requires representation of the quantum subsystem in some basis. Solutions to this equation have been already proposed in the subsystem, adiabatic and force bases, each with its own cons and pros.\n\nIn this work, the quantum-classical equations of motion are cast in the subsystem basis and subsequently mapped to a number of fictitious harmonic oscillators.\nThe result is quantum-classical dynamics in the mapping basis which treats both quantum and classical degrees of freedom on the same footing, i.e., in phase space. Neglecting a portion of the back reaction of the quantum-subsystem to classical bath results in an expression for the time evolution of an operator (density matrix) equal to its Poisson bracket with the Hamiltonian.\nThis equation can be solved in terms of characteristics to provide a computationally tractable method for calculating quantum-classical dynamical properties. The expressions for expectation values and correlation functions in this formalism are derived.\nCalculations on spin-boson system, barrier crossing models---the so called Tully models---and the Fenna-Mathews-Olson pigments show very good agreement between the results of this method and numerical solutions to the Schrödinger equation.
The non-Markovian master equation for open quantum systems is obtained by generalization of the standard Zwanzig-Nakajima (ZN) projection technique. To this end, a coupled chain of equations is written for the reduced density matrices of the bath ϱB(t) and of the system ϱS(t). A formal solution of the equation for ϱB(t) in the 2nd approximation in interaction yields a specific extra term related to the intrinsic bath dynamics. This term is nonlinear in the reduced density matrix ϱS(t), and vanishes in the Markovian limit. To verify the consistence and robustness of our approach, we apply the generalized ZN projection scheme to a simple dephasing model. We study the obtained kinetic equation both in the Markovian approximation and beyond it (for the term related to the intrinsic bath dynamics) and compare the results with the exact ones.
The non-Markovian master equation for open quantum systems is obtained by generalization of the ordinary Zwanzig-Nakajima (ZN) projection technique. To this end, a coupled chain of equations for the reduced density matrices of the bath ϱB(t)\varrho_{B}(t)ϱB​(t) and the system ϱS(t)\varrho_{S}(t)ϱS​(t) are written down. A formal solution of the equation for ϱB(t)\varrho_{B}(t)ϱB​(t), having been inserted in the equation for the reduced density matrix of the system, in the 2-nd approximation in interaction yields a very specific extra term in the generalized master equation. This term, being nonlinear in ϱS(t)\varrho_{S}(t)ϱS​(t), is related to the intrinsic bath dynamics and vanishes in the Markovian limit. To verify the consistence and robustness of our approach, we applied the generalized ZN projection scheme to a simple dephasing model. It is shown that consideration of the lowest order in interaction is insufficient to describe time evolution of the system coherence adequately. We explain this fact by analyzing the exact and approximate forms of ϱB(t)\varrho_{B}(t)ϱB​(t) and give some hints how to take the dynamic correlations (which originate from the spin-bath coupling) into account.