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:2501.09544 · Physical Review Research

Unifying quantum stochastic methods using Wick's theorem on the Keldysh contour

Vasco Cavina, Antonio D’Abbruzzo, Vittorio Giovannetti

Published Oct 22, 2025• 2 citations
arXivHTMLPDF

Abstract

We present a method, based on the Keldysh formalism, for deriving stochastic master equations that describe the non-Markovian dynamics of a quantum system coupled to a Gaussian environment. This approach yields a compact expression for the system's propagator, which we show to be equivalent to existing formulations, such as the stochastic von Neumann equation (SVNE). A key advantage of our method is its generality: It can be extended to describe any open-system evolution defined on a suitable ordering contour. As a result, we adapt it to derive generalized versions of the SVNE that account for initial system-environment correlations, as well as stochastic equations that incorporate information about the statistics of energy flows in the environment. The insights offered by our technique further allow us to examine the nature of the noise processes appearing in the SVNE. We prove that its solution can be expressed in terms of a single physical noise, without any loss of information. Finally, we propose a semiclassical scenario in which this noise can be interpreted as arising from an initial measurement process on the environment.

DOI
10.1103/cnvm-w8cy

Related papers

Exact Non-Markovian Master Equations: A Generalized Derivation for Gaussian Systems

Antonio D’Abbruzzo, Vittorio Giovannetti, Vasco Cavina

We derive an exact master equation that captures the dynamics of a quadratic quantum system linearly coupled to a Gaussian environment of the same statistics: the Gaussian master equation (GME). Unlike previous approaches, our formulation applies universally to both bosonic and fermionic setups, and remains valid even in the presence of initial system-environment correlations, allowing for the exact computation of the system's reduced density matrix across all parameter regimes. Remarkably, the GME shares the same operatorial structure as the Redfield equation and depends on a single kernel: a dressed environment correlation function accounting for all virtual interactions between the system and the environment. This simple structure grants a clear physical interpretation and makes the GME easy to simulate numerically, as we show by applying it to an open system based on two fermions coupled via superconductive pairing.

arXiv:2502.14364•Physical Review Letters•Nov 12, 2025•3 citations

A unified stochastic formulation of dissipative quantum dynamics. I. Generalized hierarchical equations

Chang‐Yu Hsieh, Jianshu Cao

We extend a standard stochastic theory to study open quantum systems coupled to a generic quantum environment. We exemplify the general framework by studying a two-level quantum system coupled bilinearly to the three fundamental classes of non-interacting particles: bosons, fermions, and spins. In this unified stochastic approach, the generalized stochastic Liouville equation (SLE) formally captures the exact quantum dissipations when noise variables with appropriate statistics for different bath models are applied. Anharmonic effects of a non-Gaussian bath are precisely encoded in the bath multi-time correlation functions that noise variables have to satisfy. Starting from the SLE, we devise a family of generalized hierarchical equations by averaging out the noise variables and expand bath multi-time correlation functions in a complete basis of orthonormal functions. The general hierarchical equations constitute systems of linear equations that provide numerically exact simulations of quantum dynamics. For bosonic bath models, our general hierarchical equation of motion reduces exactly to an extended version of hierarchical equation of motion which allows efficient simulation for arbitrary spectral densities and temperature regimes. Similar efficiency and flexibility can be achieved for the fermionic bath models within our formalism. The spin bath models can be simulated with two complementary approaches in the present formalism. (I) They can be viewed as an example of non-Gaussian bath models and be directly handled with the general hierarchical equation approach given their multi-time correlation functions. (II) Alternatively, each bath spin can be first mapped onto a pair of fermions and be treated as fermionic environments within the present formalism.

arXiv:1701.05709•The Journal of Chemical Physics•Jan 2, 2018•85 citations

Stochastic wave-function unraveling of the generalized Lindblad master equation

Mervlyn Moodley, Francesco Petruccione

Recently a generalized master equation was derived that extends the Lindblad theory to highly non-Markovian quantum processes [H.-P. Breuer, Phys. Rev. A 75, 022103 (2007)]. We perform a stochastic unraveling of this master equation by considering $n$ random state vectors that satisfy the corresponding stochastic differential equation for a piecewise deterministic process. As an application we consider a two-state system randomly coupled to an environment consisting of two energy bands with finite number of levels. Our numerical results are compared to results obtained from the time-convolutionless projection operator method using correlated projection superoperators and the exact solution of the Schr\"odinger equation for this system.

arXiv:0903.5382•Physical Review A•Apr 22, 2009•33 citations

Non-Markovian dynamics of open quantum systems

Christen H. Fleming

An open quantum system is a quantum system that interacts with some environment whose degrees of freedom have been coarse grained away. This model describes non-equilibrium processes more general than scattering-matrix formulations. Furthermore, the microscopically-derived environment provides a model of noise, dissipation and decoherence far more general than Markovian (white noise) models. The latter are fully characterized by Lindblad equations and can be motivated phenomenologically. Non-Markovian processes consistently account for backreaction with the environment and can incorporate effects such as finite temperature and spatial correlations. We consider linear systems with bilinear coupling to the environment, or quantum Brownian motion, and nonlinear systems with weak coupling to the environment. For linear systems we provide exact solutions with analytical results for a variety of spectral densities. Furthermore, we point out an important mathematical subtlety which led to incorrect master-equation coefficients in earlier derivations, given nonlocal dissipation. For nonlinear systems we provide perturbative solutions by translating the formalism of canonical perturbation theory into the context of master equations. It is shown that unavoidable degeneracy causes an unfortunate reduction in accuracy between perturbative master equations and their solutions. We also extend the famous theorem of Lindblad, Gorini, Kossakowski and Sudarshan on completely positivity to non-Markovian master equations. Our application is primarily to model atoms interacting via a common electromagnetic field. The electromagnetic field contains correlations in both space and time, which are related to its relativistic (photon-mediated) nature. As such, atoms residing in the same field experience different environmental effects depending upon their relative position and orientation. Our more accurate solutions were necessary to assess sudden death of entanglement at zero temperature. In contrast to previous claims, we found that all initial states of two-level atoms undergo finite-time disentanglement. We were also able to access regimes which cannot be described by Lindblad equations and other simpler methods, such as near resonance. Finally we revisit the infamous Abraham-Lorentz force, wherein a single particle in motion experiences backreaction from the electromagnetic field. This leads to a number of well-known problems including pre-acceleration and runaway solutions. We found a more a more-suitable open-system treatment of the nonrelativistic particle to be perfectly causal and dissipative without any extraneous requirements for finite size of the particle, weak coupling to the field, etc..

•Jan 1, 2011•6 citations

Non-Markovian reduced propagator, multiple-time correlation functions, and master equations with general initial conditions in the weak-coupling limit

Inés de Vega, Daniel Alonso

In this paper we derive the evolution equation for the reduced propagator, an object that evolves vectors of the Hilbert space of a system $\mathcal{S}$ interacting with an environment $\mathcal{B}$ in a non-Markovian way. This evolution is conditioned to certain initial and final states of the environment. Once an average over these environmental states is made, reduced propagators permit the evaluation of multiple-time correlation functions of system observables. When this average is done stochastically the reduced propagator evolves according to a stochastic Schr\"odinger equation. In addition, it is possible to obtain the evolution equations of the multiple-time correlation functions which generalize the well-known quantum regression theorem to the non-Markovian case. Here, both methods, stochastic and evolution equations, are described by assuming a weak coupling between system and environment. Finally, we show that reduced propagators can be used to obtain a master equation with general initial conditions, and not necessarily an initial vacuum state for the environment. We illustrate the theory with several examples.

Physical Review A•Feb 2, 2006•44 citations

Parametrization and Optimization of Gaussian Non-Markovian Unravelings for Open Quantum Dynamics

Nina Megier, Walter T. Strunz, Carlos Viviescas +1 more

We derive a family of Gaussian non-Markovian stochastic Schrödinger equations for the dynamics of open quantum systems. The different unravelings correspond to different choices of squeezed coherent states, reflecting different measurement schemes on the environment. Consequently, we are able to give a single shot measurement interpretation for the stochastic states and microscopic expressions for the noise correlations of the Gaussian process. By construction, the reduced dynamics of the open system does not depend on the squeezing parameters. They determine the non-Hermitian Gaussian correlation, a wide range of which are compatible with the Markov limit. We demonstrate the versatility of our results for quantum information tasks in the non-Markovian regime. In particular, by optimizing the squeezing parameters, we can tailor unravelings for improving entanglement bounds or for environment-assisted entanglement protection.

arXiv:1710.08359•Physical Review Letters•Apr 13, 2018•21 citations
We present a method, based on the Keldysh formalism, for deriving stochastic master equations that describe the non-Markovian dynamics of a quantum system coupled to a Gaussian environment. This approach yields a compact expression for the system's propagator, which we show to be equivalent to existing formulations, such as the stochastic von Neumann equation (SVNE). A key advantage of our method is its generality: It can be extended to describe any open-system evolution defined on a suitable ordering contour. As a result, we adapt it to derive generalized versions of the SVNE that account for initial system-environment correlations, as well as stochastic equations that incorporate information about the statistics of energy flows in the environment. The insights offered by our technique further allow us to examine the nature of the noise processes appearing in the SVNE. We prove that its solution can be expressed in terms of a single physical noise, without any loss of information. Finally, we propose a semiclassical scenario in which this noise can be interpreted as arising from an initial measurement process on the environment.
We derive an exact master equation that captures the dynamics of a quadratic quantum system linearly coupled to a Gaussian environment of the same statistics: the Gaussian master equation (GME). Unlike previous approaches, our formulation applies universally to both bosonic and fermionic setups, and remains valid even in the presence of initial system-environment correlations, allowing for the exact computation of the system's reduced density matrix across all parameter regimes. Remarkably, the GME shares the same operatorial structure as the Redfield equation and depends on a single kernel: a dressed environment correlation function accounting for all virtual interactions between the system and the environment. This simple structure grants a clear physical interpretation and makes the GME easy to simulate numerically, as we show by applying it to an open system based on two fermions coupled via superconductive pairing.
We extend a standard stochastic theory to study open quantum systems coupled to a generic quantum environment. We exemplify the general framework by studying a two-level quantum system coupled bilinearly to the three fundamental classes of non-interacting particles: bosons, fermions, and spins. In this unified stochastic approach, the generalized stochastic Liouville equation (SLE) formally captures the exact quantum dissipations when noise variables with appropriate statistics for different bath models are applied. Anharmonic effects of a non-Gaussian bath are precisely encoded in the bath multi-time correlation functions that noise variables have to satisfy. Starting from the SLE, we devise a family of generalized hierarchical equations by averaging out the noise variables and expand bath multi-time correlation functions in a complete basis of orthonormal functions. The general hierarchical equations constitute systems of linear equations that provide numerically exact simulations of quantum dynamics. For bosonic bath models, our general hierarchical equation of motion reduces exactly to an extended version of hierarchical equation of motion which allows efficient simulation for arbitrary spectral densities and temperature regimes. Similar efficiency and flexibility can be achieved for the fermionic bath models within our formalism. The spin bath models can be simulated with two complementary approaches in the present formalism. (I) They can be viewed as an example of non-Gaussian bath models and be directly handled with the general hierarchical equation approach given their multi-time correlation functions. (II) Alternatively, each bath spin can be first mapped onto a pair of fermions and be treated as fermionic environments within the present formalism.
Recently a generalized master equation was derived that extends the Lindblad theory to highly non-Markovian quantum processes [H.-P. Breuer, Phys. Rev. A 75, 022103 (2007)]. We perform a stochastic unraveling of this master equation by considering nnn random state vectors that satisfy the corresponding stochastic differential equation for a piecewise deterministic process. As an application we consider a two-state system randomly coupled to an environment consisting of two energy bands with finite number of levels. Our numerical results are compared to results obtained from the time-convolutionless projection operator method using correlated projection superoperators and the exact solution of the Schr\"odinger equation for this system.
An open quantum system is a quantum system that interacts with some environment whose degrees of freedom have been coarse grained away. This model describes non-equilibrium processes more general than scattering-matrix formulations. Furthermore, the microscopically-derived environment provides a model of noise, dissipation and decoherence far more general than Markovian (white noise) models. The latter are fully characterized by Lindblad equations and can be motivated phenomenologically. Non-Markovian processes consistently account for backreaction with the environment and can incorporate effects such as finite temperature and spatial correlations. We consider linear systems with bilinear coupling to the environment, or quantum Brownian motion, and nonlinear systems with weak coupling to the environment. For linear systems we provide exact solutions with analytical results for a variety of spectral densities. Furthermore, we point out an important mathematical subtlety which led to incorrect master-equation coefficients in earlier derivations, given nonlocal dissipation. For nonlinear systems we provide perturbative solutions by translating the formalism of canonical perturbation theory into the context of master equations. It is shown that unavoidable degeneracy causes an unfortunate reduction in accuracy between perturbative master equations and their solutions. We also extend the famous theorem of Lindblad, Gorini, Kossakowski and Sudarshan on completely positivity to non-Markovian master equations. Our application is primarily to model atoms interacting via a common electromagnetic field. The electromagnetic field contains correlations in both space and time, which are related to its relativistic (photon-mediated) nature. As such, atoms residing in the same field experience different environmental effects depending upon their relative position and orientation. Our more accurate solutions were necessary to assess sudden death of entanglement at zero temperature. In contrast to previous claims, we found that all initial states of two-level atoms undergo finite-time disentanglement. We were also able to access regimes which cannot be described by Lindblad equations and other simpler methods, such as near resonance. Finally we revisit the infamous Abraham-Lorentz force, wherein a single particle in motion experiences backreaction from the electromagnetic field. This leads to a number of well-known problems including pre-acceleration and runaway solutions. We found a more a more-suitable open-system treatment of the nonrelativistic particle to be perfectly causal and dissipative without any extraneous requirements for finite size of the particle, weak coupling to the field, etc..
In this paper we derive the evolution equation for the reduced propagator, an object that evolves vectors of the Hilbert space of a system S\mathcal{S}S interacting with an environment B\mathcal{B}B in a non-Markovian way. This evolution is conditioned to certain initial and final states of the environment. Once an average over these environmental states is made, reduced propagators permit the evaluation of multiple-time correlation functions of system observables. When this average is done stochastically the reduced propagator evolves according to a stochastic Schr\"odinger equation. In addition, it is possible to obtain the evolution equations of the multiple-time correlation functions which generalize the well-known quantum regression theorem to the non-Markovian case. Here, both methods, stochastic and evolution equations, are described by assuming a weak coupling between system and environment. Finally, we show that reduced propagators can be used to obtain a master equation with general initial conditions, and not necessarily an initial vacuum state for the environment. We illustrate the theory with several examples.
We derive a family of Gaussian non-Markovian stochastic Schrödinger equations for the dynamics of open quantum systems. The different unravelings correspond to different choices of squeezed coherent states, reflecting different measurement schemes on the environment. Consequently, we are able to give a single shot measurement interpretation for the stochastic states and microscopic expressions for the noise correlations of the Gaussian process. By construction, the reduced dynamics of the open system does not depend on the squeezing parameters. They determine the non-Hermitian Gaussian correlation, a wide range of which are compatible with the Markov limit. We demonstrate the versatility of our results for quantum information tasks in the non-Markovian regime. In particular, by optimizing the squeezing parameters, we can tailor unravelings for improving entanglement bounds or for environment-assisted entanglement protection.