Belavkin equation

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

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

The stochastic process ${\displaystyle y(t)}$ can be a mixture of two basic types: the Poisson (or jump) type ${\displaystyle y(t)\simeq n(t)-t}$, where ${\displaystyle n(t)}$ is a Poisson process corresponding to counting observation, and the Brownian (or diffusion) type ${\displaystyle y(t)\simeq w(t)}$, where ${\displaystyle w(t)}$ 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 ${\displaystyle \psi (t,\omega )}$ is normalized only in the mean-squared sense ${\displaystyle \int \|\psi (t,\omega )\|^{2}\mu (d\omega )=\|\psi _{0}\|=1}$, but generally ${\displaystyle \psi (t,\omega )}$ fails to be normalized for each ${\displaystyle \omega }$. The normalization of ${\displaystyle \psi (t,\omega )}$ for each ${\displaystyle \omega }$ gives the random posterior state vector ${\displaystyle \varphi (t,\omega )=\psi (t,\omega )/\|\psi (t,\omega )\|}$, the evolution of which is described by the posterior Belavkin equation, which is nonlinear, because operators ${\displaystyle L}$ and ${\displaystyle H}$ depend on ${\displaystyle \varphi }$ due to normalization. The stochastic process ${\displaystyle y(t)}$ in the posterior equation has independent increments with respect to the output probability measure ${\displaystyle P(t,d\omega )=\|\psi (t,\omega )\|^{2}\mu (d\omega )}$, but not with respect to the input measure. Belavkin also derived linear equation for unnormalized density operator ${\displaystyle \varrho (t,\omega )=\psi (t,\omega )\psi ^{\ast }(t,\omega )}$ and the corresponding nonlinear equation for the normalized random posterior density operator ${\displaystyle \rho (t,\omega )=\varphi (t,\omega )\varphi ^{\ast }(t,\omega )}$. 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 ${\displaystyle \rho (t,\omega )}$ over all random trajectories ${\displaystyle \omega }$ 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 ${\displaystyle \rho (0)=\psi _{0}\psi _{0}^{\ast }}$ into a mixed posterior state ${\displaystyle \rho (t)=\int \varphi (t,\omega )\varphi ^{\ast }(t,\omega )P(t,d\omega )}$ giving a rigorous description of the dynamics of the wavefunction collapse due to an observation or measurement.[12] [13] [14]

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 ${\displaystyle dy(t)}$ corresponds to the error ${\displaystyle e(t)dt=dw}$ (white noise in the diffusive case) of a noisy observation ${\displaystyle Y(t)=X(t)+e(t)}$ of operator ${\displaystyle X(t)=\lambda [L(t)+L^{\ast }(t)]}$ with the accuracy coefficient ${\displaystyle \lambda >0}$, then the indirect observation perturbs the dynamics of the system by a stochastic force ${\displaystyle f(t)}$, called the Langevin force, which is another white noise of intensity ${\displaystyle (\hbar \lambda )^{2}}$ that does not commute with the error ${\displaystyle e(t)}$. The result of such a perturbation is that the output process ${\displaystyle Y(t)}$ is commutative ${\displaystyle [Y(s),Y(t)]=0}$, and hence ${\displaystyle \int _{0}^{t}Y(r)dr}$ corresponds to a classical observation, while the system operators ${\displaystyle X}$ satisfy the non-demolition condition: all future observables must commute with the past observations (but not with the future observations): ${\displaystyle [X(s),Y(t)]=0}$ for all ${\displaystyle s\geq t}$ (but not ${\displaystyle s<t}$). Note that commutation of ${\displaystyle X(s)}$ with ${\displaystyle Y(t)}$ and another operator ${\displaystyle Z(s)}$ with ${\displaystyle Y(t)}$ does not imply commutation of ${\displaystyle X(s)}$ with ${\displaystyle Z(s)}$, 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 ${\displaystyle \mathbb {E} \{X(s)\mid Y(t)\}}$, which makes the quantum filtering possible.[16]

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

${\displaystyle L={\hat {x}},\qquad H={\frac {{\hat {p}}^{2}}{2m}}}$

the posterior stochastic equation becomes

${\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]}$

where ${\displaystyle \langle {\hat {x}}\rangle (t)=\mathrm {Tr} \left\{{\hat {x}}\rho (t)\right\}}$ is the posterior expectation of ${\displaystyle {\hat {x}}}$. Motivated by the spontaneous collapse theory rather than the filtering theory, this equation was also obtained by Diosi,[17] showing that the measurement noise ${\displaystyle dy-2\langle {\hat {x}}\rangle dt}$ is the increment ${\displaystyle d\xi }$ 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 ${\displaystyle \psi _{0}}$ these solutions correspond to optimal quantum linear filter.[15] Solutions to the Belavkin equation show that in the limit ${\displaystyle t\to \infty }$ the wavefunction has finite dispersion,[20] therefore resolving the quantum Zeno effect.[11]