Numerical simulation of stochastic evolution equations associated to quantum Markov semigroups

Abstract

We address the problem of approximating numerically the solutions (X t : t ∈[0, T]) of stochastic evolution equations on Hilbert spaces (Heng hooktop sign, 〈., .〉), with respect to Brownian motions, arising in the unraveling of backward quantum master equations. In particular, we study the computation of mean values of (X t , AX t ), where A is a linear operator. First, we introduce estimates on the behavior of X t . Then we characterize the error induced by the substitution of X t with the solution X t,n of a convenient stochastic ordinary differential equation. It allows us to establish the rate of convergence of E 〈X̃ t,n , AX̃ t,n 〉 to E 〈X t , AX t 〉, where X̃ t,n denotes the explicit Euler method. Finally, we consider an extrapolation method based on the Euler scheme. An application to the quantum harmonic oscillator system is included.