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:quant-ph/0208169 · Physical Review A

Perturbative approach to non-Markovian stochastic Schrödinger equations

Jay Gambetta, Howard M. Wiseman

Published Nov 18, 2002• 32 citations
arXivHTMLPDF

Abstract

In this paper we present a perturbative procedure that allows one to numerically solve diffusive non-Markovian stochastic Schr\"odinger equations, for a wide range of memory functions. To illustrate this procedure numerical results are presented for a classically driven two-level atom immersed in an environment with a simple memory function. It is observed that as the order of the perturbation is increased the numerical results for the ensemble average state ${\ensuremath{\rho}}_{\mathrm{red}}(t)$ approach the exact reduced state found via $\mathrm{Imamo}\mathrm{g}\ifmmode\bar\else\textasciimacron\fi{}{\mathrm{lu}}^{\ensuremath{'}}\mathrm{s}$ enlarged system method [Phys. Rev. A 50, 3650 (1994)].

DOI
10.1103/physreva.66.052105

Related papers

Perturbation expansions of stochastic wavefunctions for open quantum systems

Yaling Ke, Yi Zhao

Based on the stochastic unravelling of the reduced density operator in the Feynman path integral formalism for an open quantum system in touch with harmonic environments, a new non-Markovian stochastic Schrödinger equation (NMSSE) has been established that allows for the systematic perturbation expansion in the system-bath coupling to arbitrary order. This NMSSE can be transformed in a facile manner into the other two NMSSEs, i.e., non-Markovian quantum state diffusion and time-dependent wavepacket diffusion method. Benchmarked by numerically exact results, we have conducted a comparative study of the proposed method in its lowest order approximation, with perturbative quantum master equations in the symmetric spin-boson model and the realistic Fenna-Matthews-Olson complex. It is found that our method outperforms the second-order time-convolutionless quantum master equation in the whole parameter regime and even far better than the fourth-order in the slow bath and high temperature cases. Besides, the method is applicable on an equal footing for any kind of spectral density function and is expected to be a powerful tool to explore the quantum dynamics of large-scale systems, benefiting from the wavefunction framework and the time-local appearance within a single stochastic trajectory.

The Journal of Chemical Physics•Nov 8, 2017•21 citations

An alternative realization of the exact non-Markovian stochastic Schrödinger equation

Kai Song, Linze Song, Qiang Shi

Based on the path integral approach, we derive a new realization of the exact non-Markovian stochastic Schrödinger equation (SSE). The main difference from the previous non-Markovian quantum state diffusion (NMQSD) method is that the complex Gaussian stochastic process used for the forward propagation of the wave function is correlated, which may be used to reduce the amplitude of the non-Markovian memory term at high temperatures. The new SSE is then written into the recently developed hierarchy of pure states scheme, in a form that is more closely related to the hierarchical equation of motion approach. Numerical simulations are then performed to demonstrate the efficiency of the new method.

The Journal of Chemical Physics•Jun 10, 2016•26 citations

Non-Markovian stochastic Schrödinger equations in different temperature regimes: A study of the spin-boson model

Inés de Vega, Daniel Alonso, Pierre Gaspard +1 more

Stochastic Schrodinger equations are used to describe the dynamics of a quantum open system in contact with a large environment, as an alternative to the commonly used master equations. We present a study of the two main types of non-Markovian stochastic Schrodinger equations, linear and nonlinear ones. We compare them both analytically and numerically, the latter for the case of a spin-boson model. We show in this paper that two linear stochastic Schrodinger equations, derived from different perspectives by Diosi, Gisin, and Strunz [Phys. Rev. A 58, 1699 (1998)], and Gaspard and Nagaoka [J. Chem. Phys. 13, 5676 (1999)], respectively, are equivalent in the relevant order of perturbation theory. Nonlinear stochastic Schrodinger equations are in principle more efficient than linear ones, as they determine solutions with a higher weight in the ensemble average which recovers the reduced density matrix of the quantum open system. However, it will be shown in this paper that for the case of a spin-boson system and weak coupling, this improvement does only occur in the case of a bath at high temperature. For low temperatures, the sampling of realizations of the nonlinear equation is practically equivalent to the sampling of the linear ones. We study further this result by analyzing, for both temperature regimes, the driving noise of the linear equations in comparison to that of the nonlinear equations.

The Journal of Chemical Physics•Mar 22, 2005•40 citations

A General First- and Second-Order Numerical Solver for Non-Markovian Quantum State Diffusion

Zhenning Cai, Quanhui Zhu

The numerical simulation of non-Markovian open quantum systems based on the non-Markovian quantum state diffusion (NMQSD) equation is complicated by functional derivatives with respect to the stochastic process. A general numerical framework that directly treats these functional derivatives without relying on prescribed decompositions of the bath correlation function is still lacking. In this work, we derive an analytical solution of the linear NMQSD equation that reveals three elementary structures of the non-Markovian stochastic dynamics: stochastic propagation, functional-derivative insertion, and memory pairing. Based on this structure, we construct a general auxiliary-state framework for arbitrary bath correlation functions. The framework separates the numerical construction into time discretization, memory quadrature, and hierarchy truncation. We then construct first- and second-order schemes and provide diagrammatic transition rules for their explicit implementation. Numerical results verify the expected temporal accuracy and demonstrate the applicability of the proposed methods to different bath correlation functions and multi-level quantum systems.

quant-pharXiv:2607.25427•Jul 29, 2026•0 citations

Markovian unraveling of non-Markovian stochastic Schr\"odinger equations

Antoine Tilloy

Non-Markovian stochastic Schrodinger equations (NMSSE) are important tools in quantum mechanics, from the theory of open systems to foundations. Yet, in general, they are but formal objects: their solution can be computed numerically only in some specific cases or perturbatively. This article is focused on the NMSSE themselves rather than on the open-system evolution they unravel and aims at making them less abstract. Namely, we propose to write the stochastic realizations of linear NMSSE as averages over the solutions of an auxiliary equation with an additional random field. Our method yields a non-perturbative numerical simulation algorithm for generic linear NMSSE that can be made arbitrarily accurate for reasonably short times. For isotropic complex noises, the method extends from linear to non-linear NMSSE and allows to sample the solutions of norm-preserving NMSSE directly.

•Jan 8, 2017•0 citations

Time-local unraveling of non-Markovian stochastic Schrödinger equations

Antoine Tilloy

Non-Markovian stochastic Schrödinger equations (NMSSE) are important tools in quantum mechanics, from the theory of open systems to foundations. Yet, in general, they are but formal objects: their solution can be computed numerically only in some specific cases or perturbatively. This article is focused on the NMSSE themselves rather than on the open-system evolution they unravel and aims at making them less abstract. Namely, we propose to write the stochastic realizations of linear NMSSE as averages over the solutions of an auxiliary equation with an additional random field. Our method yields a non-perturbative numerical simulation algorithm for generic linear NMSSE that can be made arbitrarily accurate for reasonably short times. For isotropic complex noises, the method extends from linear to non-linear NMSSE and allows to sample the solutions of norm-preserving NMSSE directly.

arXiv:1701.01948•Quantum•Sep 19, 2017•10 citations
In this paper we present a perturbative procedure that allows one to numerically solve diffusive non-Markovian stochastic Schr"odinger equations, for a wide range of memory functions. To illustrate this procedure numerical results are presented for a classically driven two-level atom immersed in an environment with a simple memory function. It is observed that as the order of the perturbation is increased the numerical results for the ensemble average state \ensuremathρred(t){\ensuremath{\rho}}_{\mathrm{red}}(t)\ensuremathρred​(t) approach the exact reduced state found via Imamog\ifmmode\elseˉ\textasciimacron\filu\ensuremath′s\mathrm{Imamo}\mathrm{g}\ifmmode\bar\else\textasciimacron\fi{}{\mathrm{lu}}^{\ensuremath{'}}\mathrm{s}Imamog\ifmmode\elseˉ​\textasciimacron\filu\ensuremath′s enlarged system method [Phys. Rev. A 50, 3650 (1994)].
Based on the stochastic unravelling of the reduced density operator in the Feynman path integral formalism for an open quantum system in touch with harmonic environments, a new non-Markovian stochastic Schrödinger equation (NMSSE) has been established that allows for the systematic perturbation expansion in the system-bath coupling to arbitrary order. This NMSSE can be transformed in a facile manner into the other two NMSSEs, i.e., non-Markovian quantum state diffusion and time-dependent wavepacket diffusion method. Benchmarked by numerically exact results, we have conducted a comparative study of the proposed method in its lowest order approximation, with perturbative quantum master equations in the symmetric spin-boson model and the realistic Fenna-Matthews-Olson complex. It is found that our method outperforms the second-order time-convolutionless quantum master equation in the whole parameter regime and even far better than the fourth-order in the slow bath and high temperature cases. Besides, the method is applicable on an equal footing for any kind of spectral density function and is expected to be a powerful tool to explore the quantum dynamics of large-scale systems, benefiting from the wavefunction framework and the time-local appearance within a single stochastic trajectory.
Based on the path integral approach, we derive a new realization of the exact non-Markovian stochastic Schrödinger equation (SSE). The main difference from the previous non-Markovian quantum state diffusion (NMQSD) method is that the complex Gaussian stochastic process used for the forward propagation of the wave function is correlated, which may be used to reduce the amplitude of the non-Markovian memory term at high temperatures. The new SSE is then written into the recently developed hierarchy of pure states scheme, in a form that is more closely related to the hierarchical equation of motion approach. Numerical simulations are then performed to demonstrate the efficiency of the new method.
Stochastic Schrodinger equations are used to describe the dynamics of a quantum open system in contact with a large environment, as an alternative to the commonly used master equations. We present a study of the two main types of non-Markovian stochastic Schrodinger equations, linear and nonlinear ones. We compare them both analytically and numerically, the latter for the case of a spin-boson model. We show in this paper that two linear stochastic Schrodinger equations, derived from different perspectives by Diosi, Gisin, and Strunz [Phys. Rev. A 58, 1699 (1998)], and Gaspard and Nagaoka [J. Chem. Phys. 13, 5676 (1999)], respectively, are equivalent in the relevant order of perturbation theory. Nonlinear stochastic Schrodinger equations are in principle more efficient than linear ones, as they determine solutions with a higher weight in the ensemble average which recovers the reduced density matrix of the quantum open system. However, it will be shown in this paper that for the case of a spin-boson system and weak coupling, this improvement does only occur in the case of a bath at high temperature. For low temperatures, the sampling of realizations of the nonlinear equation is practically equivalent to the sampling of the linear ones. We study further this result by analyzing, for both temperature regimes, the driving noise of the linear equations in comparison to that of the nonlinear equations.
The numerical simulation of non-Markovian open quantum systems based on the non-Markovian quantum state diffusion (NMQSD) equation is complicated by functional derivatives with respect to the stochastic process. A general numerical framework that directly treats these functional derivatives without relying on prescribed decompositions of the bath correlation function is still lacking. In this work, we derive an analytical solution of the linear NMQSD equation that reveals three elementary structures of the non-Markovian stochastic dynamics: stochastic propagation, functional-derivative insertion, and memory pairing. Based on this structure, we construct a general auxiliary-state framework for arbitrary bath correlation functions. The framework separates the numerical construction into time discretization, memory quadrature, and hierarchy truncation. We then construct first- and second-order schemes and provide diagrammatic transition rules for their explicit implementation. Numerical results verify the expected temporal accuracy and demonstrate the applicability of the proposed methods to different bath correlation functions and multi-level quantum systems.
Non-Markovian stochastic Schrodinger equations (NMSSE) are important tools in quantum mechanics, from the theory of open systems to foundations. Yet, in general, they are but formal objects: their solution can be computed numerically only in some specific cases or perturbatively. This article is focused on the NMSSE themselves rather than on the open-system evolution they unravel and aims at making them less abstract. Namely, we propose to write the stochastic realizations of linear NMSSE as averages over the solutions of an auxiliary equation with an additional random field. Our method yields a non-perturbative numerical simulation algorithm for generic linear NMSSE that can be made arbitrarily accurate for reasonably short times. For isotropic complex noises, the method extends from linear to non-linear NMSSE and allows to sample the solutions of norm-preserving NMSSE directly.
Non-Markovian stochastic Schrödinger equations (NMSSE) are important tools in quantum mechanics, from the theory of open systems to foundations. Yet, in general, they are but formal objects: their solution can be computed numerically only in some specific cases or perturbatively. This article is focused on the NMSSE themselves rather than on the open-system evolution they unravel and aims at making them less abstract. Namely, we propose to write the stochastic realizations of linear NMSSE as averages over the solutions of an auxiliary equation with an additional random field. Our method yields a non-perturbative numerical simulation algorithm for generic linear NMSSE that can be made arbitrarily accurate for reasonably short times. For isotropic complex noises, the method extends from linear to non-linear NMSSE and allows to sample the solutions of norm-preserving NMSSE directly.