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:2511.05282 · The Journal of Chemical Physics

Open quantum-classical systems: A hybrid MASH master equation

Kasra Asnaashari, Jeremy O. Richardson

quant-phphysics.chem-phSubmitted Nov 7, 2025• Announced Nov 10, 2025• 2 citations• Revised Mar 29, 2026 (v2)
v1 · Nov 7, 2025v2 · Mar 29, 2026
arXivHTMLPDF

Abstract

We propose a method which combines the quantum-classical mapping approach to surface hopping (MASH) with the dissipative quantum dynamics of the Lindblad master equation. Like conventional surface-hopping methods, our approach is based on classical trajectories coupled to the dynamics of a quantum subsystem. However, instead of evolving the subsystem wavefunction according to the time-dependent Schrödinger equation, we use stochastic quantum trajectories derived from secular Redfield theory. This enables the simulation of open quantum systems coupled simultaneously to Markovian quantum baths and anharmonic non-Markovian classical degrees of freedom. Applications to the spin--boson model and to the cavity-enhanced fluorescence of an electronically nonadiabatic molecule show excellent agreement with fully quantum-mechanical benchmarks.

Journal reference
J. Chem. Phys. 164, 024116 (2026)
DOI
10.1063/5.0311156
Comments
14 pages, 5 figures

Related papers

Quantum Trajectories for Time-Local Non-Lindblad Master Equations

Tobias Becker, Ché Netzer, André Eckardt

For the efficient simulation of open quantum systems, we often use quantum jump trajectories given by pure states that evolve stochastically to unravel the dynamics of the underlying master equation. In the Markovian regime, when the dynamics is described by a Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) master equation, this procedure is known as Monte Carlo wave function approach. However, beyond ultraweak system-bath coupling, the dynamics of the system is not described by an equation of GKSL type, but rather by the Redfield equation, which can be brought into pseudo-Lindblad form. Here, negative dissipation strengths prohibit the conventional approach. To overcome this problem, we propose a pseudo-Lindblad quantum trajectory (PLQT) unraveling. It does not require an effective extension of the state space, like other approaches, except for the addition of a single classical bit. We test the PLQT for the eternal non-Markovian master equation for a single qubit and an interacting Fermi-Hubbard chain coupled to a thermal bath and discuss its computational effort compared to solving the full master equation.

arXiv:2306.14876•Physical Review Letters•Oct 17, 2023•11 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

Stochastic Bloch-Redfield theory: Quantum jumps in a solid-state environment

Nicolas Vogt, Jan Jeske, Jared H. Cole

We discuss mapping the Bloch-Redfield master equation to Lindblad form and then unraveling the resulting evolution into a stochastic Schr\"odinger equation according to the quantum-jump method. We give two approximations under which this mapping is valid. This approach enables us to study solid-state systems of much larger sizes than is possible with the standard Bloch-Redfield master equation, while still providing a systematic method for obtaining the jump operators and corresponding rates. We also show how the stochastic unraveling of the Bloch-Redfield equations becomes the kinetic Monte Carlo algorithm in the secular approximation when the system-bath-coupling operators are given by tunneling operators between system eigenstates. The stochastic unraveling is compared to the conventional Bloch-Redfield approach with the superconducting single-electron transistor (SSET) as an example.

arXiv:1310.0559•Physical Review B•Nov 19, 2013•24 citations

Decoherence and quantum-classical master equation dynamics

Robbie Grunwald, Raymond Kapral

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.

arXiv:cond-mat/0612203•The Journal of Chemical Physics•Mar 20, 2007•23 citations

Quantum-classical Liouville dynamics in the mapping basis

Hyojoon Kim, Ali Nassimi, Raymond Kapral

The quantum-classical Liouville equation describes the dynamics of a quantum subsystem coupled to a classical environment. It has been simulated using various methods, notably, surface-hopping schemes. A representation of this equation in the mapping Hamiltonian basis for the quantum subsystem is derived. The resulting equation of motion, in conjunction with expressions for quantum expectation values in the mapping basis, provides another route to the computation of the nonadiabatic dynamics of observables that does not involve surface-hopping dynamics. The quantum-classical Liouville equation is exact for the spin-boson system. This well-known model is simulated using an approximation to the evolution equation in the mapping basis, and close agreement with exact quantum results is found.

arXiv:0903.0847•The Journal of Chemical Physics•Aug 25, 2008•156 citations

Universal Lindblad equation for open quantum systems

Frederik Nathan, Mark Spencer Rudner

We develop a Markovian master equation in the Lindblad form that enables the efficient study of a wide range of open quantum many-body systems that would be inaccessible with existing methods. The validity of the master equation is based entirely on properties of the bath and the system-bath coupling, without any requirements on the level structure within the system itself. The master equation is derived using a Markov approximation that is distinct from that used in earlier approaches. We provide a rigorous bound for the error induced by this Markov approximation; the error is controlled by a dimensionless combination of intrinsic correlation and relaxation timescales of the bath. Our master equation is accurate on the same level of approximation as the Bloch-Redfield equation. In contrast to the Bloch-Redfield approach, our approach ensures preservation of the positivity of the density matrix. As a result, our method is robust, and can be solved efficiently using stochastic evolution of pure states (rather than density matrices). We discuss how our method can be applied to static or driven quantum many-body systems, and illustrate its power through numerical simulation of a spin chain that would be challenging to treat by existing methods.

arXiv:2004.01469•Physical review. B./Physical review. B•Sep 3, 2020•212 citations
We propose a method which combines the quantum-classical mapping approach to surface hopping (MASH) with the dissipative quantum dynamics of the Lindblad master equation. Like conventional surface-hopping methods, our approach is based on classical trajectories coupled to the dynamics of a quantum subsystem. However, instead of evolving the subsystem wavefunction according to the time-dependent Schrödinger equation, we use stochastic quantum trajectories derived from secular Redfield theory. This enables the simulation of open quantum systems coupled simultaneously to Markovian quantum baths and anharmonic non-Markovian classical degrees of freedom. Applications to the spin--boson model and to the cavity-enhanced fluorescence of an electronically nonadiabatic molecule show excellent agreement with fully quantum-mechanical benchmarks.
For the efficient simulation of open quantum systems, we often use quantum jump trajectories given by pure states that evolve stochastically to unravel the dynamics of the underlying master equation. In the Markovian regime, when the dynamics is described by a Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) master equation, this procedure is known as Monte Carlo wave function approach. However, beyond ultraweak system-bath coupling, the dynamics of the system is not described by an equation of GKSL type, but rather by the Redfield equation, which can be brought into pseudo-Lindblad form. Here, negative dissipation strengths prohibit the conventional approach. To overcome this problem, we propose a pseudo-Lindblad quantum trajectory (PLQT) unraveling. It does not require an effective extension of the state space, like other approaches, except for the addition of a single classical bit. We test the PLQT for the eternal non-Markovian master equation for a single qubit and an interacting Fermi-Hubbard chain coupled to a thermal bath and discuss its computational effort compared to solving the full master equation.
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.
We discuss mapping the Bloch-Redfield master equation to Lindblad form and then unraveling the resulting evolution into a stochastic Schr\"odinger equation according to the quantum-jump method. We give two approximations under which this mapping is valid. This approach enables us to study solid-state systems of much larger sizes than is possible with the standard Bloch-Redfield master equation, while still providing a systematic method for obtaining the jump operators and corresponding rates. We also show how the stochastic unraveling of the Bloch-Redfield equations becomes the kinetic Monte Carlo algorithm in the secular approximation when the system-bath-coupling operators are given by tunneling operators between system eigenstates. The stochastic unraveling is compared to the conventional Bloch-Redfield approach with the superconducting single-electron transistor (SSET) as an example.
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.
The quantum-classical Liouville equation describes the dynamics of a quantum subsystem coupled to a classical environment. It has been simulated using various methods, notably, surface-hopping schemes. A representation of this equation in the mapping Hamiltonian basis for the quantum subsystem is derived. The resulting equation of motion, in conjunction with expressions for quantum expectation values in the mapping basis, provides another route to the computation of the nonadiabatic dynamics of observables that does not involve surface-hopping dynamics. The quantum-classical Liouville equation is exact for the spin-boson system. This well-known model is simulated using an approximation to the evolution equation in the mapping basis, and close agreement with exact quantum results is found.
We develop a Markovian master equation in the Lindblad form that enables the efficient study of a wide range of open quantum many-body systems that would be inaccessible with existing methods. The validity of the master equation is based entirely on properties of the bath and the system-bath coupling, without any requirements on the level structure within the system itself. The master equation is derived using a Markov approximation that is distinct from that used in earlier approaches. We provide a rigorous bound for the error induced by this Markov approximation; the error is controlled by a dimensionless combination of intrinsic correlation and relaxation timescales of the bath. Our master equation is accurate on the same level of approximation as the Bloch-Redfield equation. In contrast to the Bloch-Redfield approach, our approach ensures preservation of the positivity of the density matrix. As a result, our method is robust, and can be solved efficiently using stochastic evolution of pure states (rather than density matrices). We discuss how our method can be applied to static or driven quantum many-body systems, and illustrate its power through numerical simulation of a spin chain that would be challenging to treat by existing methods.