Belavkin equation

From HandWiki

In quantum probability, the Belavkin equation, also known as Belavkin-Schrödinger equation, quantum filtering equation, stochastic master equation, is a quantum stochastic differential equation describing the dynamics of a quantum system undergoing observation in continuous time. It was derived and henceforth studied by Viacheslav Belavkin in 1988.[1][2][3]


Unlike the Schrödinger equation, which describes the deterministic evolution of the wavefunction [math]\displaystyle{ \psi(t) }[/math] of a closed system (without interaction), the Belavkin equation describes the stochastic evolution of a random wavefunction [math]\displaystyle{ \psi(t,\omega) }[/math] of an open quantum system interacting with an observer:

[math]\displaystyle{ d\psi=-\left(\frac{1}{2} L^\ast L+\frac{i}{\hbar}H\right)\psi dt+L\psi dy }[/math]

Here, [math]\displaystyle{ L }[/math] is a self-adjoint operator (or a column vector of operators) of the system coupled to the external field, [math]\displaystyle{ H }[/math] is the Hamiltonian, [math]\displaystyle{ i=\sqrt{-1} }[/math] is the imaginary unit, [math]\displaystyle{ \hbar }[/math] is the Planck constant, and [math]\displaystyle{ y(t)=\int_0^t dy(r) }[/math] is a stochastic process representing the measurement noise that is a martingale with independent increments with respect to the input probability measure [math]\displaystyle{ P(0,d\omega)=\|\psi(0,\omega)\|^2\mu(d\omega) }[/math]. Note that this noise has dependent increments with respect to the output probability measure [math]\displaystyle{ P(t,d\omega)=\|\psi(t,\omega)\|^2\mu(d\omega) }[/math] representing the output innovation process (the observation). For [math]\displaystyle{ L=0 }[/math], the equation becomes the standard Schrödinger equation.

The stochastic process [math]\displaystyle{ y(t) }[/math] can be a mixture of two basic types: the Poisson (or jump) type [math]\displaystyle{ y(t)\simeq n(t)-t }[/math], where [math]\displaystyle{ n(t) }[/math] is a Poisson process corresponding to counting observation, and the Brownian (or diffusion) type [math]\displaystyle{ y(t)\simeq w(t) }[/math], where [math]\displaystyle{ w(t) }[/math] is the standard Wiener process corresponding to continuous observation. The equations of the diffusion type can be derived as the central limit of the jump type equations with the expected rate of the jumps increasing to infinity.

The random wavefunction [math]\displaystyle{ \psi(t,\omega) }[/math] is normalized only in the mean-squared sense [math]\displaystyle{ \int\|\psi(t,\omega)\|^2\mu(d\omega)=\|\psi_0\|=1 }[/math], but generally [math]\displaystyle{ \psi(t,\omega) }[/math] fails to be normalized for each [math]\displaystyle{ \omega }[/math]. The normalization of [math]\displaystyle{ \psi(t,\omega) }[/math] for each [math]\displaystyle{ \omega }[/math] gives the random posterior state vector [math]\displaystyle{ \varphi(t,\omega)=\psi(t,\omega)/\|\psi(t,\omega)\| }[/math], the evolution of which is described by the posterior Belavkin equation, which is nonlinear, because operators [math]\displaystyle{ L }[/math] and [math]\displaystyle{ H }[/math] depend on [math]\displaystyle{ \varphi }[/math] due to normalization. The stochastic process [math]\displaystyle{ y(t) }[/math] in the posterior equation has independent increments with respect to the output probability measure [math]\displaystyle{ P(t,d\omega)=\|\psi(t,\omega)\|^2\mu(d\omega) }[/math], but not with respect to the input measure. Belavkin also derived linear equation for unnormalized density operator [math]\displaystyle{ \varrho(t,\omega)=\psi(t,\omega)\psi^\ast(t,\omega) }[/math] and the corresponding nonlinear equation for the normalized random posterior density operator [math]\displaystyle{ \rho(t,\omega)=\varphi(t,\omega)\varphi^\ast(t,\omega) }[/math]. For two types of measurement noise, this gives eight basic quantum stochastic differential equations. The general forms of the equations include all types of noise and their representations in Fock space.[4] [5]

The nonlinear equation describing observation of position of a free particle, which is a special case of the posterior Belavkin equation of the diffusion type, was also obtained by Diosi[6] and appeared in the works of Gisin,[7] Ghirardi, Pearle and Rimini,[8] although with a rather different motivation or interpretation. Similar nonlinear equations for posterior density operators were postulated (although without derivation) in quantum optics and the quantum trajectories theory,[9] where they are called stochastic master equations. The averaging of the equations for the random density operators [math]\displaystyle{ \rho(t,\omega) }[/math] over all random trajectories [math]\displaystyle{ \omega }[/math] leads to the Lindblad equation,[10] which is deterministic.

The nonlinear Belavkin equations for posterior states play the same role as the Stratonovich–Kushner equation in classical probability, while the linear equations correspond to the Zakai equation.[11] The Belavkin equations describe continuous-time decoherence of initially pure state [math]\displaystyle{ \rho(0)=\psi_0\psi_0^\ast }[/math] into a mixed posterior state [math]\displaystyle{ \rho(t)=\int\varphi(t,\omega)\varphi^\ast(t,\omega)P(t,d\omega) }[/math] giving a rigorous description of the dynamics of the wavefunction collapse due to an observation or measurement.[12] [13] [14]

Non-demolition measurement and quantum filtering

Noncommutativity presents a major challenge for probabilistic interpretation of quantum stochastic differential equations due to non-existence of conditional expectations for general pairs of quantum observables. Belavkin resolved this issue by discovering the error-perturbation uncertainty relation and formulating the non-demolition principle of quantum measurement.[13][15] In particular, if the stochastic process [math]\displaystyle{ dy(t) }[/math] corresponds to the error [math]\displaystyle{ e(t)dt=dw }[/math] (white noise in the diffusive case) of a noisy observation [math]\displaystyle{ Y(t)=X(t)+e(t) }[/math] of operator [math]\displaystyle{ X(t)=\lambda[L(t)+L^\ast(t)] }[/math] with the accuracy coefficient [math]\displaystyle{ \lambda\gt 0 }[/math], then the indirect observation perturbs the dynamics of the system by a stochastic force [math]\displaystyle{ f(t) }[/math], called the Langevin force, which is another white noise of intensity [math]\displaystyle{ (\hbar\lambda)^2 }[/math] that does not commute with the error [math]\displaystyle{ e(t) }[/math]. The result of such a perturbation is that the output process [math]\displaystyle{ Y(t) }[/math] is commutative [math]\displaystyle{ [Y(s),Y(t)]=0 }[/math], and hence [math]\displaystyle{ \int_0^tY(r)dr }[/math] corresponds to a classical observation, while the system operators [math]\displaystyle{ X }[/math] satisfy the non-demolition condition: all future observables must commute with the past observations (but not with the future observations): [math]\displaystyle{ [X(s),Y(t)]=0 }[/math] for all [math]\displaystyle{ s\geq t }[/math] (but not [math]\displaystyle{ s\lt t }[/math]). Note that commutation of [math]\displaystyle{ X(s) }[/math] with [math]\displaystyle{ Y(t) }[/math] and another operator [math]\displaystyle{ Z(s) }[/math] with [math]\displaystyle{ Y(t) }[/math] does not imply commutation of [math]\displaystyle{ X(s) }[/math] with [math]\displaystyle{ Z(s) }[/math], so that the algebra of future observables is still non-commutative. The non-demolition condition is necessary and sufficient for the existence of conditional expectations [math]\displaystyle{ \mathbb{E}\{X(s)\mid Y(t)\} }[/math], which makes the quantum filtering possible.[16]

Posterior state equations

Counting observation

Let [math]\displaystyle{ n(t) }[/math] be a Poisson process with forward increments [math]\displaystyle{ dn(t)=0 }[/math] almost everywhere and [math]\displaystyle{ dn(t)=1 }[/math] otherwise and having the property [math]\displaystyle{ (dn)^2=dn }[/math]. The expected number of events is [math]\displaystyle{ \mathbb{E}\{n(t)\}=\nu t }[/math], where [math]\displaystyle{ \nu }[/math] is the expected rate of jumps. Then substituting [math]\displaystyle{ m(t)=\frac{n(t)-\nu t}{\sqrt{\nu}} }[/math] for the stochastic process [math]\displaystyle{ y(t) }[/math] gives the linear Belavkin equation for the unnormalized random wavefunction [math]\displaystyle{ \psi(t,\omega) }[/math] undergoing counting observation. Substituting [math]\displaystyle{ L=\sqrt\nu(C-I) }[/math], where [math]\displaystyle{ C }[/math] is the collapse operator, and [math]\displaystyle{ H=E+i\hbar\frac{\nu}{2}(C-C^\ast) }[/math], where [math]\displaystyle{ E }[/math] is the energy operator, this equation can be written in the following form

[math]\displaystyle{ d\psi=-\left(\frac{\nu}{2}(C^\ast C-I)+\frac{i}{\hbar} E\right)\psi dt+(C-I)\psi dn }[/math]

Normalized wavefunction [math]\displaystyle{ \varphi(t,\omega)=\psi(t,\omega)/\|\psi(t,\omega)\| }[/math] is called the posterior state vector, the evolution of which is described by the following nonlinear equation

[math]\displaystyle{ d\varphi=-\left(\frac{\nu}{2}(C^\ast C-\|C\varphi\|^2)+\frac{i}{\hbar} E\right)\varphi dt+\left(\frac{C}{\|C\varphi\|}-I\right)\varphi d\tilde n }[/math]

where [math]\displaystyle{ \tilde n(t) }[/math] has expectation [math]\displaystyle{ \mathbb{E}\{\tilde n(t)\}=\nu\|C\varphi\|^2t }[/math]. The posterior equation can be written in the standard form

[math]\displaystyle{ d\varphi=-\left(\frac{1}{2} \tilde L^\ast\tilde L+\frac{i}{\hbar}\tilde H\right)\varphi dt+\tilde L\varphi d\tilde m }[/math]

with [math]\displaystyle{ \tilde L=\sqrt\nu(C-\|C\varphi\|) }[/math], [math]\displaystyle{ \tilde H=E+i\hbar\frac{\nu}{2}(C-C^\ast)\|C\varphi\| }[/math], and [math]\displaystyle{ \tilde m(t)=\frac{\tilde n(t)-\nu\|C\varphi\|^2 t}{\sqrt\nu\|C\varphi\|} }[/math]. The corresponding equations for the unnormalized random density operator [math]\displaystyle{ \varrho(t,\omega)=\psi(t,\omega)\psi^\ast(t,\omega) }[/math] and for the normalized random posterior density operator [math]\displaystyle{ \rho(t,\omega)=\varphi(t,\omega)\varphi^\ast(t,\omega) }[/math] are as follows

[math]\displaystyle{ \begin{align} d\varrho&=-[G\varrho+\varrho G^\ast-\nu\varrho]dt+[C\varrho C^\ast-\varrho]dn\\ d\rho&=-[G\rho+\rho G^\ast-\nu\rho\mathrm{Tr}\left\{C\rho C^\ast\right\}]dt+\left(\frac{C\rho C^\ast}{\mathrm{Tr}\left\{C\rho C^\ast\right\}}-\rho\right)d\tilde n \end{align} }[/math]

where [math]\displaystyle{ G=\frac{\nu}{2} C^\ast C+\frac{i}{\hbar} E }[/math]. Note that the latter equation is nonlinear.

Continuous observation

Stochastic process [math]\displaystyle{ m(t) }[/math], defined in the previous section, has forward increments [math]\displaystyle{ (dm)^2=\nu^{-1/2}dm+dt }[/math], which tend to [math]\displaystyle{ dt }[/math] as [math]\displaystyle{ \nu\to\infty }[/math]. Therefore, [math]\displaystyle{ m(t) }[/math] becomes standard Wiener process with respect to the input probability measure. Substituting [math]\displaystyle{ w(t) }[/math] for [math]\displaystyle{ y(t) }[/math] gives the linear Belavkin equation for the unnormalized random wavefunction [math]\displaystyle{ \psi(t,\omega) }[/math] undergoing continuous observation. The output process [math]\displaystyle{ \tilde m(t) }[/math] becomes the diffusion innovation process [math]\displaystyle{ \tilde w(t) }[/math] with increments [math]\displaystyle{ d\tilde w(t,\omega)=dw(t,\omega)-2\mathrm{Re}\langle\varphi(t,\omega)\mid L\varphi(t,\omega)\rangle dt }[/math]. The nonlinear Belavkin equation of the diffusion type for the posterior state vector [math]\displaystyle{ \varphi(t,\omega)=\psi(t,\omega)/\|\psi(t,\omega)\| }[/math] is

[math]\displaystyle{ d\varphi=-\left(\frac{1}{2} \tilde L^\ast\tilde L+\frac{i}{\hbar}\tilde H\right)\varphi dt+\tilde L\varphi d\tilde w }[/math]

with [math]\displaystyle{ \tilde L=L-\mathrm{Re}\langle\varphi\mid L\varphi\rangle }[/math] and [math]\displaystyle{ \tilde H=H+i\hbar\frac{1}{2}(L-L^\ast)\mathrm{Re}\langle\varphi\mid L\varphi\rangle }[/math]. The corresponding equations for the unnormalized random density operator [math]\displaystyle{ \varrho(t,\omega)=\psi(t,\omega)\psi^\ast(t,\omega) }[/math] and for the normalized random posterior density operator [math]\displaystyle{ \rho(t,\omega)=\varphi(t,\omega)\varphi^\ast(t,\omega) }[/math] are as follows

[math]\displaystyle{ \begin{align} d\varrho&=-[K\varrho+\varrho K^\ast-L\varrho L^\ast]dt+[L\varrho+\varrho L^\ast]dw\\ d\rho&=-[K\rho+\rho K^\ast-L\rho L^\ast]dt+[L\rho+\rho L^\ast-\rho\mathrm{Tr}\left\{(L+L^\ast)\rho\right\}]d\tilde w \end{align} }[/math]

where [math]\displaystyle{ K=\frac{1}{2} L^\ast L+\frac{i}{\hbar} H }[/math]. The second equation is nonlinear due to normalization. Because [math]\displaystyle{ \mathbb{E}\{dw(t)\}=0 }[/math], taking the average of these stochastic equations over all [math]\displaystyle{ \omega }[/math] leads to the Lindblad equation

[math]\displaystyle{ \frac{d\rho}{dt}=-[K\rho+\rho K^\ast]+L\rho L^\ast }[/math]

Example: continuous observation of position of a free particle

Consider a free particle of mass [math]\displaystyle{ m }[/math]. The position and momentum observables correspond respectively to operators [math]\displaystyle{ \hat x }[/math] of multiplication by [math]\displaystyle{ x }[/math] and [math]\displaystyle{ \hat p=(\hbar/i)d/dx }[/math]. Making the following substitutions in the Belavkin equation

[math]\displaystyle{ L=\hat x,\qquad H=\frac{\hat p^2}{2m} }[/math]

the posterior stochastic equation becomes

[math]\displaystyle{ d\varphi=-\left(\frac{1}{2}(\hat x-\langle\hat x\rangle)^2+\frac{i}{\hbar}\frac{\hat p^2}{2m}\right)\varphi dt+\Bigl(\hat x-\langle\hat x\rangle\Bigr)\varphi[dy-2\langle\hat x\rangle dt] }[/math]

where [math]\displaystyle{ \langle\hat x\rangle(t)=\mathrm{Tr}\left\{\hat x\rho(t)\right\} }[/math] is the posterior expectation of [math]\displaystyle{ \hat x }[/math]. Motivated by the spontaneous collapse theory rather than the filtering theory, this equation was also obtained by Diosi,[17] showing that the measurement noise [math]\displaystyle{ dy-2\langle\hat x\rangle dt }[/math] is the increment [math]\displaystyle{ d\xi }[/math] of a standard Wiener process. There are closed-form solutions to this equation,[18] as well as equations for a particle in a linear or quadratic potentials.[1][3][19] For a Gaussian initial state [math]\displaystyle{ \psi_0 }[/math] these solutions correspond to optimal quantum linear filter.[15] Solutions to the Belavkin equation show that in the limit [math]\displaystyle{ t\to\infty }[/math] the wavefunction has finite dispersion,[20] therefore resolving the quantum Zeno effect.[11]


  1. 1.0 1.1 Belavkin, V.P. (1988). "Nondemolition measurements, nonlinear filtering and dynamic programming of quantum stochastic processes". in A. Blaquiere. Proc of Bellmann Continuum Workshop `Modelling and Control of Systems'. Lecture notes in Control and Inform Sciences. 121. Sophia-Antipolis: Springer-Verlag. pp. 245–265. 
  2. Belavkin, V.P. (1989). "A continuous counting observation and posterior quantum dynamics". J Phys A 22 (23): L1109–L1114. doi:10.1088/0305-4470/22/23/006. Bibcode1989JPhA...22L1109B. 
  3. 3.0 3.1 Belavkin, V.P. (1989). "A new wave equation for a continuous nondemolition measurement". Physics Letters A 140 (7–8): 355–358. doi:10.1016/0375-9601(89)90066-2. Bibcode1989PhLA..140..355B. 
  4. Belavkin, V.P. (1995). "On stochastic generators of completely positive cocycles". Russ Journ of Math Phys 3 (4): 523–528. 
  5. Belavkin, V.P. (1997). "Quantum stochastic positive evolutions: characterization, construction, dilation". Commun. Math. Phys. 184 (3): 533–566. doi:10.1007/s002200050072. Bibcode1997CMaPh.184..533B. 
  6. Di'osi, L. (1989). "Models for universal reduction of macroscopic quantum fluctuations". Physical Review A 40 (3): 1165–1174. doi:10.1103/PhysRevA.40.1165. PMID 9902248. Bibcode1989PhRvA..40.1165D. 
  7. Gisin, N. (1989). "Stochastic quantum dynamics and relativity". Helvetica Physica Acta 62: 363–371. 
  8. Ghirardi, G.C.; Pearle, P.; Rimini, A. (1990). "Markov processes in Hilbert space and continuous spontaneous localization of systems of identical particles". Phys. Rev. A 42 (1): 78–89. doi:10.1103/PhysRevA.42.78. PMID 9903779. Bibcode1990PhRvA..42...78G. 
  9. Carmichael, H.J. (1993). An Open Systems Approach to Quantum Optics. Springer-Verlag. 
  10. Smolyanov, O.; Truman, A. (1999). "Schrödinger-Belavkin equations and associated Kolmogorov and Lindblad equations". Theoretical and Mathematical Physics 120 (2): 973–984. doi:10.1007/BF02557405. Bibcode1999TMP...120..973S. 
  11. 11.0 11.1 Holevo, A.S. (1991). "Quantum probability and quantum statistics". in Prokhorov, Y.V. (in ru). Probability Theory~8. Itogi Nauki i Tekhniki. 83. VINITI. pp. 5–132. 
  12. Belavkin, V.P. (1990). "A quantum posterior stochastics and spontaneous collapse". in Truman, A.. Stochastics and Quantum Mechanics. World Scientific. pp. 40–68. 
  13. 13.0 13.1 Belavkin, V.P. (1992). "Quantum continual measurements and a posteriori collapse on CCR.". Comm. Math. Phys. 146 (3): 611–635. doi:10.1007/BF02097018. Bibcode1992CMaPh.146..611B. 
  14. Belavkin, V.P.; Melsheimer, O. (1995). "A hamiltonian solution to quantum collapse, state diffusion and spontaneous localization". Quantum Communications and Measurement. Plenum Publisher. pp. 201–222. doi:10.1007/978-1-4899-1391-3_20. ISBN 978-1-4899-1393-7. 
  15. 15.0 15.1 Belavkin, V.P. (1980). "Optimal filtering of Markov signals with quantum white noise". Radio Eng Electron Physics 25: 1445–1453. doi:10.1007/978-1-4899-1391-3_37. 
  16. Bouten, L.; van Handel, R.; James, M.R. (2009). "A discrete invitation to quantum filtering and feedback control". SIAM Review 51 (2): 239–316. doi:10.1137/060671504. Bibcode2009SIAMR..51..239B. 
  17. Diosi, L. (1988). "Continuous quantum measurement and Itô formalism". Phys Lett A 129 (8–9): 419–423. doi:10.1016/0375-9601(88)90309-X. Bibcode1988PhLA..129..419D. 
  18. Diosi, L. (1988). "Localized solution of simple nonlinear quantum Langevin-equation". Phys Lett A 132 (5): 233–236. doi:10.1016/0375-9601(88)90555-5. Bibcode1988PhLA..132..233D. 
  19. Belavkin, V.P.; Staszewski, P. (1992). "Nondemolition observation of a free quantum particle". Physical Review A 45 (3): 1347–1357. doi:10.1103/PhysRevA.45.1347. PMID 9907114. Bibcode1992PhRvA..45.1347B. 
  20. Kolokol'tsov1, V.N. (1995). "Scattering theory for the Belavkin equation describing a quantum particle with continuously observed coordinate". Journal of Mathematical Physics 36 (6): 2741–2760. doi:10.1063/1.531063. Bibcode1995JMP....36.2741K.