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:2004.02387 · Physical Review A

Solving quantum trajectories for systems with linear Heisenberg-picture dynamics and Gaussian measurement noise

Prahlad Warszawski, Howard M. Wiseman, Andrew C. Doherty

Published Oct 12, 2020• 12 citations
arXivHTMLPDF

Abstract

We study solutions to the quantum trajectory evolution of $N$-mode open quantum systems possessing a time-independent Hamiltonian, linear Heisenberg-picture dynamics, and Gaussian measurement noise. In terms of the mode annihilation and creation operators, a system will have linear Heisenberg-picture dynamics under two conditions. First, the Hamiltonian must be quadratic. Second, the Lindblad operators describing the coupling to the environment (including those corresponding to the measurement) must be linear. In cases where we can solve the $2N$-degree polynomials that arise in our calculations, we provide an analytical solution for initial states that are arbitrary (i.e., they are not required to have Gaussian Wigner functions). The solution takes the form of an evolution operator, with the measurement-result dependence captured in $2N$ stochastic integrals over these classical random signals. The solutions also allow the positive operator-valued measure (POVM), which generates the probabilities of obtaining measurement outcomes, to be determined. To illustrate our results, we solve some single-mode example systems, with the POVMs being of practical relevance to the inference of an initial state, via quantum state tomography. Our key tool is the representation of mixed states of quantum mechanical oscillators as state vectors rather than state matrices (albeit in a larger Hilbert space). Together with methods from Lie algebra, this allows a more straightforward manipulation of the exponential operators comprising the system evolution than is possible in the original Hilbert space.

DOI
10.1103/physreva.102.042210

Related papers

Quantum-jump vs stochastic Schrödinger dynamics for Gaussian states with quadratic Hamiltonians and linear Lindbladians

Robson Christie, Jessica K. Eastman, Roman Schubert +1 more

Abstract The dynamics of Gaussian states for open quantum systems described by Lindblad equations can be solved analytically for systems with quadratic Hamiltonians and linear Lindbladians, showing the familiar phenomena of dissipation and decoherence. It is well known that the Lindblad dynamics can be expressed as an ensemble average over stochastic pure-state dynamics, which can be interpreted as individual experimental implementations, where the form of the stochastic dynamics depends on the measurement setup. Here we consider quantum-jump and stochastic Schrödinger dynamics for initially Gaussian states. While both unravellings converge to the same Lindblad dynamics when averaged, the individual dynamics can differ qualitatively. For the stochastic Schrödinger equation, Gaussian states remain Gaussian during the evolution, with stochastic differential equations governing the evolution of the phase-space centre and a deterministic evolution of the covariance matrix. In contrast to this, individual pure-state dynamics arising from the quantum-jump evolution do not remain Gaussian in general. Applying results developed in the non-Hermitian context for Hagedorn wavepackets, we formulate a method to generate quantum-jump trajectories that is described entirely in terms of the evolution of an underlying Gaussian state. To illustrate the behaviours of the different unravellings in comparison to the Lindblad dynamics, we consider two examples in detail, which can be largely treated analytically, a harmonic oscillator subject to position measurement and a damped harmonic oscillator. In both cases, we highlight the differences as well as the similarities of the stochastic Schrödinger and the quantum-jump dynamics.

arXiv:2203.11530•Journal of Physics A Mathematical and Theoretical•Oct 27, 2022•12 citations

On Non-Gaussian Linear Quantum Optics

Alexander Evgen’evich Teretenkov

Abstract This paper studies a class of single-mode bosonic models of open quantum systems in which the dynamics of the creation and annihilation operators remains linear despite the strongly non-Gaussian nature of the induced dissipation. It is shown that, in contrast to isolated bosonic systems—where linear dynamics uniquely corresponds to quadratic Hamiltonians and Gaussian states—linearity of the operator dynamics in open systems does not guarantee Gaussianity of the state evolution. We construct quantum channels and Gorini–Kossakowski–Sudarshan–Lindblad generators for which the dynamics of the creation and annihilation operators is linear and closed, while Gaussian states are driven beyond the convex hull of Gaussian states during the evolution. In particular, the emergence of negativity of the Wigner function at a finite evolution time is demonstrated.

Bulletin of the Russian Academy of Sciences Physics•Jul 1, 2026•1 citation

Bohmian Trajectories Within Hilbert Space Based Quantum Mechanics. Solution of the Measurement Problem

Tulsi Dass

de Broglie-Bohm theory (dBBT), treating quantum particles as point objects moving along well defined (Bohmian) trajectories, offers an appealing solution of the measurement problem in quantum mechanics; it has, however, problems relating to spin, relativity and lack of proper integration with the Hilbert space based framework. In this work, we present a consistent formalism which has the traditional state-observable framework integrated with the desirable features of dBBT. We adopt ensemble interpretation for the Schrodinger wave function $ψ$. Given a Schrodinger wave function $ψ$, we use its value $ψ_0$ at some fixed time (say, $t = 0$) to define the probability measure $|ψ_0|^2 {\rm d}x$ on the system configuration space $M$ ($=\mathbb{R}^n$). On the resulting probability space $\mathcal{M}_0$, we introduce a stochastic process $ξ(t)$ corresponding to the Heisenberg position operator $X_H(t)$ such that, in the Heisenberg state $|ψ_h\rangle$ corresponding to $ψ_0$, the expectation value of $X_H(t)$ equals that of $ξ(t)$ in $\mathcal{M}_0$. This condition leads to the de Broglie-Bohm guidance equation for the sample paths of the process $ξ(t)$ which are, therefore, Bohmian trajectories supposedly representing time-evolutions of individual members of the $ψ_0$-ensemble. Stochastic processes and Bohmian trajectories corresponding to observables with discrete eigenvalues (in particular spin) are treated by extending the configuration space to the spectral space of the commutative algebra obtained by adding appropriate discrete observables to the position observables. Pauli's equation is treated as an example. A straightforward derivation of von Neumann's projection rule employing the Schrodinger-Bohm evolution of individual systems along their Bohmian trajectories is given. Some comments on the potential application of the formalism developed here to quantum mechanics of the universe are included.

quant-pharXiv:2512.07007•Dec 9, 2025•0 citations

Diffusive Dynamics and Continuous Monitoring

Alessio Serafini

In the last chapter, we dealt with the wide class of physical operations that preserve the Gaussian character of the initial state, treating them as discrete transformations. Here, we describe what could loosely be considered the continuous-time version of such transformations. This chapter is intended to provide the reader with the tools to describe a wide range of notable quantum dynamics, with a very broad applicative scope. After a note on the action of linear and quadratic Hamiltonians, whose description follows directly from the theory of Gaussian states presented in Chapter 3 but was not explicitly clarified so far, we shall address systems linearly coupled to a reservoir of modes that are delta-correlated in time (a situation commonly referred to as ‘white noise’), and thus give rise to Markovian dynamics which we shall refer to as ‘diffusive’ for reasons that will be made clear in due time. Our story will commence with a derivation of the Gaussian dynamics resulting from such an environmental coupling, in the form of dynamical equations for the first and second statistical moments. Besides their interest per se , the latter will be instrumental in obtaining an unconventional derivation of the master equations that govern general, not necessarily Gaussian, quantum states under the same conditions. The corresponding quantum Langevin equations for the canonical operators in the Heisenberg picture will also be derived. Next, we will turn our attention to the study of the system's dynamics that result from the measurement of the environment. This is a topic with a longstanding tradition in quantum mechanics, grounded in the formalism of quantum stochastic calculus that emerged within the mathematical physics community as early as the seventies, and then re-introduced as the quantum trajectories approach in quantum optics almost twenty years later. The ‘weak’ POVMs resulting from observing a weakly coupled environment (as opposed to ‘strong’, projective measurements), carried out continuously in time (a prescription customarily subsumed by the term ‘monitoring’), leads to a stochastic dynamics for the system, where the randomness of the measurement outcomes is compactly and most elegantly represented by infinitesimal stochastic increments, giving rise to continuous but not differentiable quantum dynamics. Existing approaches to this subject tend to be quite impervious to the uninitiated. In the quantum optics literature, this is largely due to the technicalities involved in incorporating the discontinuous update of the quantum state prescribed by the von Neumann postulate. More mathematically oriented treatments have instead typically preferred an algebraic approach where the usual difficulty of describing measurements in the Heisenberg picture is compounded by the adoption of a rather abstract hierarchy of filters, in part inspired by classical control theory, to describe different possible updates of the quantum state. This state of affairs is further complicated by the added difficulties that come with the necessary introduction of some form of quantum stochastic calculus. Our approach to monitoring will therefore take a different route. Like in our derivation of diffusion and master equations, we shall keep restricting to linear couplings with a white noise environment, forfeiting full generality This situation is well modelled by the so called input-output formalism, whose physical origin will be discussed in detail in Chapter 9 , which approximates reliably, among other instances, the interaction between internal resonating modes and electromagnetic radiation leaking in and out of an optical cavity. First, we shall introduce stochastic Wiener increments and the Itô rule, which will arise naturally from the identification of the continuum field's modes subject to monitoring. We shall deliberately keep our arsenal of stochastic methods to these bare essentials, sacrificing to some extent both depth and generality, but still providing the reader with a powerful and seminal demonstration of the formalism. We shall then restrict to Gaussian general-dyne detections and apply the exact formalism developed in the previous chapter to the continuous filtering of the system's state conditioned by such measurements' outcomes. Next, as in our derivation of the unconditional master equations, we shall use the Gaussian dynamics as a key to unlock the full stochastic master equations governing general states. Hence, we will build on the quantum theory of filtering in the Gaussian regime to derive dynamical equations for linear feedback control and also apply our theoretical tools to the enhanced production of squeezing. Finally, we will turn to the closely related paradigm of coherent feedback, where environmental modes do not undergo measurements but are fed back to interact coherently with the system, see in detail how it modifies Gaussian diffusive dynamics and assess its potential for control by reconsidering the production of squeezing assisted by interferometric coherent feedback loops.

CRC Press eBooks•Jun 1, 2023•0 citations

On the Unraveling of Open Quantum Dynamics

Brecht Donvil, Paolo Muratore-Ginanneschi

It is well known that the state operator of an open quantum system can be generically represented as the solution of a time-local equation — a quantum master equation. Unraveling in quantum trajectories offers a picture of open system dynamics dual to solving master equations. In the unraveling picture, physical indicators are computed as Monte Carlo averages over a stochastic process valued in the Hilbert space of the system. This approach is particularly adapted to simulate systems in large Hilbert spaces. We show that the dynamics of an open quantum system generically admits an unraveling in the Hilbert space of the system described by a Markov process generated by ordinary stochastic differential equations for which rigorous concentration estimates are available. The unraveling can be equivalently formulated in terms of norm-preserving state vectors or in terms of linear “ostensible” processes trace preserving only on average. We illustrate the results in the case of a two level system in a simple boson environment. Next, we derive the state-of-the-art form of the Diósi-Gisin-Strunz Gaussian random ostensible state equation in the context of a model problem. This equation provides an exact unraveling of open systems in Gaussian environments. We compare and contrast the two unravelings and their potential for applications to quantum error mitigation.

arXiv:2309.13408•Open Systems & Information Dynamics•Sep 1, 2023•5 citations

Measurement-based Initial Point Smoothing and Control Approach to Quantum Memory Systems

Igor G. Vladimirov, Ian R. Petersen, Guodong Shi

This paper is concerned with a quantum memory system for storing quantum information in the form of its initial dynamic variables in the presence of environmental noise. In order to compensate for the deviation from the initial conditions, the classical parameters of the system Hamiltonian are affected by the actuator output of a measurement-based classical controller. The latter uses an observation process produced by a measuring apparatus from the quantum output field of the memory system. The underlying system is modelled as an open quantum harmonic oscillator whose Heisenberg evolution is governed by linear Hudson-Parthasarathy quantum stochastic differential equations. The controller is organised as a classical linear time-varying system, so that the resulting closed-loop system has quantum and classical dynamic variables. We apply linear-quadratic-Gaussian control and fixed-point smoothing at the level of the first two moments and consider controllers with a separation structure which involve a continuously updated estimate for the initial quantum variables. The initial-point smoother is used for actuator signal formation so as to minimise the sum of a mean-square deviation of the quantum memory system variables at a given time horizon from their initial values and an integral quadratic penalty on the control signal.

quant-pharXiv:2512.05586•Dec 8, 2025•0 citations
We study solutions to the quantum trajectory evolution of NNN-mode open quantum systems possessing a time-independent Hamiltonian, linear Heisenberg-picture dynamics, and Gaussian measurement noise. In terms of the mode annihilation and creation operators, a system will have linear Heisenberg-picture dynamics under two conditions. First, the Hamiltonian must be quadratic. Second, the Lindblad operators describing the coupling to the environment (including those corresponding to the measurement) must be linear. In cases where we can solve the 2N2N2N-degree polynomials that arise in our calculations, we provide an analytical solution for initial states that are arbitrary (i.e., they are not required to have Gaussian Wigner functions). The solution takes the form of an evolution operator, with the measurement-result dependence captured in 2N2N2N stochastic integrals over these classical random signals. The solutions also allow the positive operator-valued measure (POVM), which generates the probabilities of obtaining measurement outcomes, to be determined. To illustrate our results, we solve some single-mode example systems, with the POVMs being of practical relevance to the inference of an initial state, via quantum state tomography. Our key tool is the representation of mixed states of quantum mechanical oscillators as state vectors rather than state matrices (albeit in a larger Hilbert space). Together with methods from Lie algebra, this allows a more straightforward manipulation of the exponential operators comprising the system evolution than is possible in the original Hilbert space.
Abstract The dynamics of Gaussian states for open quantum systems described by Lindblad equations can be solved analytically for systems with quadratic Hamiltonians and linear Lindbladians, showing the familiar phenomena of dissipation and decoherence. It is well known that the Lindblad dynamics can be expressed as an ensemble average over stochastic pure-state dynamics, which can be interpreted as individual experimental implementations, where the form of the stochastic dynamics depends on the measurement setup. Here we consider quantum-jump and stochastic Schrödinger dynamics for initially Gaussian states. While both unravellings converge to the same Lindblad dynamics when averaged, the individual dynamics can differ qualitatively. For the stochastic Schrödinger equation, Gaussian states remain Gaussian during the evolution, with stochastic differential equations governing the evolution of the phase-space centre and a deterministic evolution of the covariance matrix. In contrast to this, individual pure-state dynamics arising from the quantum-jump evolution do not remain Gaussian in general. Applying results developed in the non-Hermitian context for Hagedorn wavepackets, we formulate a method to generate quantum-jump trajectories that is described entirely in terms of the evolution of an underlying Gaussian state. To illustrate the behaviours of the different unravellings in comparison to the Lindblad dynamics, we consider two examples in detail, which can be largely treated analytically, a harmonic oscillator subject to position measurement and a damped harmonic oscillator. In both cases, we highlight the differences as well as the similarities of the stochastic Schrödinger and the quantum-jump dynamics.
Abstract This paper studies a class of single-mode bosonic models of open quantum systems in which the dynamics of the creation and annihilation operators remains linear despite the strongly non-Gaussian nature of the induced dissipation. It is shown that, in contrast to isolated bosonic systems—where linear dynamics uniquely corresponds to quadratic Hamiltonians and Gaussian states—linearity of the operator dynamics in open systems does not guarantee Gaussianity of the state evolution. We construct quantum channels and Gorini–Kossakowski–Sudarshan–Lindblad generators for which the dynamics of the creation and annihilation operators is linear and closed, while Gaussian states are driven beyond the convex hull of Gaussian states during the evolution. In particular, the emergence of negativity of the Wigner function at a finite evolution time is demonstrated.
de Broglie-Bohm theory (dBBT), treating quantum particles as point objects moving along well defined (Bohmian) trajectories, offers an appealing solution of the measurement problem in quantum mechanics; it has, however, problems relating to spin, relativity and lack of proper integration with the Hilbert space based framework. In this work, we present a consistent formalism which has the traditional state-observable framework integrated with the desirable features of dBBT. We adopt ensemble interpretation for the Schrodinger wave function ψψψ. Given a Schrodinger wave function ψψψ, we use its value ψ0ψ_0ψ0​ at some fixed time (say, t=0t = 0t=0) to define the probability measure ∣ψ0∣2dx|ψ_0|^2 {\rm d}x∣ψ0​∣2dx on the system configuration space MMM (=Rn=\mathbb{R}^n=Rn). On the resulting probability space M0\mathcal{M}_0M0​, we introduce a stochastic process ξ(t)ξ(t)ξ(t) corresponding to the Heisenberg position operator XH(t)X_H(t)XH​(t) such that, in the Heisenberg state ∣ψh⟩|ψ_h\rangle∣ψh​⟩ corresponding to ψ0ψ_0ψ0​, the expectation value of XH(t)X_H(t)XH​(t) equals that of ξ(t)ξ(t)ξ(t) in M0\mathcal{M}_0M0​. This condition leads to the de Broglie-Bohm guidance equation for the sample paths of the process ξ(t)ξ(t)ξ(t) which are, therefore, Bohmian trajectories supposedly representing time-evolutions of individual members of the ψ0ψ_0ψ0​-ensemble. Stochastic processes and Bohmian trajectories corresponding to observables with discrete eigenvalues (in particular spin) are treated by extending the configuration space to the spectral space of the commutative algebra obtained by adding appropriate discrete observables to the position observables. Pauli's equation is treated as an example. A straightforward derivation of von Neumann's projection rule employing the Schrodinger-Bohm evolution of individual systems along their Bohmian trajectories is given. Some comments on the potential application of the formalism developed here to quantum mechanics of the universe are included.
In the last chapter, we dealt with the wide class of physical operations that preserve the Gaussian character of the initial state, treating them as discrete transformations. Here, we describe what could loosely be considered the continuous-time version of such transformations. This chapter is intended to provide the reader with the tools to describe a wide range of notable quantum dynamics, with a very broad applicative scope. After a note on the action of linear and quadratic Hamiltonians, whose description follows directly from the theory of Gaussian states presented in Chapter 3 but was not explicitly clarified so far, we shall address systems linearly coupled to a reservoir of modes that are delta-correlated in time (a situation commonly referred to as ‘white noise’), and thus give rise to Markovian dynamics which we shall refer to as ‘diffusive’ for reasons that will be made clear in due time. Our story will commence with a derivation of the Gaussian dynamics resulting from such an environmental coupling, in the form of dynamical equations for the first and second statistical moments. Besides their interest per se , the latter will be instrumental in obtaining an unconventional derivation of the master equations that govern general, not necessarily Gaussian, quantum states under the same conditions. The corresponding quantum Langevin equations for the canonical operators in the Heisenberg picture will also be derived. Next, we will turn our attention to the study of the system's dynamics that result from the measurement of the environment. This is a topic with a longstanding tradition in quantum mechanics, grounded in the formalism of quantum stochastic calculus that emerged within the mathematical physics community as early as the seventies, and then re-introduced as the quantum trajectories approach in quantum optics almost twenty years later. The ‘weak’ POVMs resulting from observing a weakly coupled environment (as opposed to ‘strong’, projective measurements), carried out continuously in time (a prescription customarily subsumed by the term ‘monitoring’), leads to a stochastic dynamics for the system, where the randomness of the measurement outcomes is compactly and most elegantly represented by infinitesimal stochastic increments, giving rise to continuous but not differentiable quantum dynamics. Existing approaches to this subject tend to be quite impervious to the uninitiated. In the quantum optics literature, this is largely due to the technicalities involved in incorporating the discontinuous update of the quantum state prescribed by the von Neumann postulate. More mathematically oriented treatments have instead typically preferred an algebraic approach where the usual difficulty of describing measurements in the Heisenberg picture is compounded by the adoption of a rather abstract hierarchy of filters, in part inspired by classical control theory, to describe different possible updates of the quantum state. This state of affairs is further complicated by the added difficulties that come with the necessary introduction of some form of quantum stochastic calculus. Our approach to monitoring will therefore take a different route. Like in our derivation of diffusion and master equations, we shall keep restricting to linear couplings with a white noise environment, forfeiting full generality This situation is well modelled by the so called input-output formalism, whose physical origin will be discussed in detail in Chapter 9 , which approximates reliably, among other instances, the interaction between internal resonating modes and electromagnetic radiation leaking in and out of an optical cavity. First, we shall introduce stochastic Wiener increments and the Itô rule, which will arise naturally from the identification of the continuum field's modes subject to monitoring. We shall deliberately keep our arsenal of stochastic methods to these bare essentials, sacrificing to some extent both depth and generality, but still providing the reader with a powerful and seminal demonstration of the formalism. We shall then restrict to Gaussian general-dyne detections and apply the exact formalism developed in the previous chapter to the continuous filtering of the system's state conditioned by such measurements' outcomes. Next, as in our derivation of the unconditional master equations, we shall use the Gaussian dynamics as a key to unlock the full stochastic master equations governing general states. Hence, we will build on the quantum theory of filtering in the Gaussian regime to derive dynamical equations for linear feedback control and also apply our theoretical tools to the enhanced production of squeezing. Finally, we will turn to the closely related paradigm of coherent feedback, where environmental modes do not undergo measurements but are fed back to interact coherently with the system, see in detail how it modifies Gaussian diffusive dynamics and assess its potential for control by reconsidering the production of squeezing assisted by interferometric coherent feedback loops.
It is well known that the state operator of an open quantum system can be generically represented as the solution of a time-local equation — a quantum master equation. Unraveling in quantum trajectories offers a picture of open system dynamics dual to solving master equations. In the unraveling picture, physical indicators are computed as Monte Carlo averages over a stochastic process valued in the Hilbert space of the system. This approach is particularly adapted to simulate systems in large Hilbert spaces. We show that the dynamics of an open quantum system generically admits an unraveling in the Hilbert space of the system described by a Markov process generated by ordinary stochastic differential equations for which rigorous concentration estimates are available. The unraveling can be equivalently formulated in terms of norm-preserving state vectors or in terms of linear “ostensible” processes trace preserving only on average. We illustrate the results in the case of a two level system in a simple boson environment. Next, we derive the state-of-the-art form of the Diósi-Gisin-Strunz Gaussian random ostensible state equation in the context of a model problem. This equation provides an exact unraveling of open systems in Gaussian environments. We compare and contrast the two unravelings and their potential for applications to quantum error mitigation.
This paper is concerned with a quantum memory system for storing quantum information in the form of its initial dynamic variables in the presence of environmental noise. In order to compensate for the deviation from the initial conditions, the classical parameters of the system Hamiltonian are affected by the actuator output of a measurement-based classical controller. The latter uses an observation process produced by a measuring apparatus from the quantum output field of the memory system. The underlying system is modelled as an open quantum harmonic oscillator whose Heisenberg evolution is governed by linear Hudson-Parthasarathy quantum stochastic differential equations. The controller is organised as a classical linear time-varying system, so that the resulting closed-loop system has quantum and classical dynamic variables. We apply linear-quadratic-Gaussian control and fixed-point smoothing at the level of the first two moments and consider controllers with a separation structure which involve a continuously updated estimate for the initial quantum variables. The initial-point smoother is used for actuator signal formation so as to minimise the sum of a mean-square deviation of the quantum memory system variables at a given time horizon from their initial values and an integral quadratic penalty on the control signal.