About this Journal Submit a Manuscript Table of Contents
The Scientific World Journal
Volume 2013 (2013), Article ID 424137, 6 pages
Research Article

Simulation of Quantum Dynamics Based on the Quantum Stochastic Differential Equation

School of Automation, Guangdong University of Technology, No. 100 Waihuan Xi Road, Guangzhou Higher Education Mega Center, Pan Yu District, Guangzhou, Guangdong 510006, China

Received 11 April 2013; Accepted 28 April 2013

Academic Editors: C. Chen, D. Dong, M. Jiang, and L.-C. Wang

Copyright © 2013 Ming Li. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


The quantum stochastic differential equation derived from the Lindblad form quantum master equation is investigated. The general formulation in terms of environment operators representing the quantum state diffusion is given. The numerical simulation algorithm of stochastic process of direct photodetection of a driven two-level system for the predictions of the dynamical behavior is proposed. The effectiveness and superiority of the algorithm are verified by the performance analysis of the accuracy and the computational cost in comparison with the classical Runge-Kutta algorithm.

1. Introduction

Since Nelson successfully described the kinematics law of the quantum fluctuations by the Itô equation [1] and the Schrödinger equation was derived from Newtonian mechanics; the stochastic interpretation of quantum mechanics was established, in which a diffusion process was used to analyze the quantum fluctuation instead of the wave function. Then, the stochastic mechanics has gradually drawn much attention with research fields ranging from atomic and optical physics to condensed matter physics and quantum information science [2]. It becomes clear that a deep understanding of the effects of environments on a quantum system such as the mechanisms of decoherence and the dynamics of entanglement in the framework of quantum open systems is both of fundamental interest in quantum foundation issues and of practical importance in quantum information sciences. Many scholars have made thorough research on the quantum diffusion movement based on the basic theory and achieved fruitful results. For example, the normative structure of the dynamics equation of the particle fluctuation and its stability analysis methods were determined, followed by the unified interpretation of the Brown motion and the basic equation of the quantum mechanics [35]. The quantum stochastic dynamics elaborated the organic link between the microscopic behavior and macroevolution of the system, by which the details of the system evolution from any initial state to the final state can be analyzed. It has been successfully applied to the proliferation of microscopic particles, the molecular motors, the quantum chaos, and so forth [6].

Though the deterministic differential equations of quantum stochastic mechanics are relatively complex, it can be seen that with the development of the computer technology [7], the evolution of a microsystem can be analyzed using the numerical simulation method [810]. The quantum trajectory method is a typical one. It can be used for a wide range of open quantum systems to solve the master equation by unraveling the density operator evolution into individual stochastic trajectories in Hilbert space [11]. Over the last twenty years the theory of quantum trajectories has been developed by many researchers for a variety of purposes, including modeling continuously monitored open systems, improving numerical calculation and investigating the problem of quantum measurement [12, 13].

In this paper, based on the quantum stochastic dynamics, the master equation describing the time evolution law of the quantum state and its reduced density operator are investigated, and the effect of nonunitary operators on the evolution of the system is analyzed. Then the quantum stochastic differential equation is established to describe the microkinetic characteristics of the system, and a numerical iterative algorithm for the simulation of the system evolution is proposed. The practicality and advantages of the algorithm are verified by comparison with the classical Runge-Kutta numerical iterative algorithm, which is followed by further discussions on the convergence of the algorithm.

2. Methods

The quantum state diffusion theory replaces the deterministic evolution of the density operator representing an ensemble of open systems [14] by a unique stochastic diffusion of a quantum state, representing an individual system of the ensemble in interaction with its environment. is the Hamiltonian, and are a set of environment operators which represent the collective effects of interaction with the environment.

However, for some complicated systems, it can be very difficult to get either the analytical or the numerical solution of (1). In that case, it is often advantageous to take alternative ways considering an unraveling of the master equation into individual quantum trajectories. Quantum state diffusion (QSD) is one of these unraveling techniques. The corresponding quantum state diffusion equation is a stochastic differential equation for the normalized state vector representing the pure state of the system that evolves according to the QSD equation. The differerential Itô form is [15] where are defined by

In the QSD equation, the standard normalized terms represent independent complex Wiener processes and satisfy the relations where denotes an ensemble average of the noise. QSD reproduces the master equation in the mean That is to say, the reduced state of the system, , is obtained as an ensemble average. And this is what is meant by an unraveling of the master equation. Expectation values for operators obey a similar relationship: The use of QSD as a practical algorithm to solve master equations has been widely investigated [16, 17]. This includes calculations of output spectra in quantum optics [18]. As a practical method of computation, QSD gains over the direct solution of the master equation, because of a basis of states, QSD needs a computer store with elements, and the time of computation is also proportional to . For the direct solution these are proportional to .

In this paper, we take advantage of the system simulation method to simulate the evolutionary behavior of the open quantum systems and thus calculate and analyze various physical properties of the ensemble of open quantum systems. However, noting that we cannot usually get the analytical solution of (2), an alternative way is to find the numerical solution of the system evolution and investigate various control algorithms and control strategies based on the simulation method, which is a powerful tool built on the systems science, system identification, control theory, and computer technology for the analysis and synthesis of complex systems, especially large-scale systems [19].

In the system simulation, we should pay attention to the problem that the physical description of the stochastic process is relied on the master equation in the Lindblad operator form. If we want to get the numerical solution, the Lindblad operator must be explicitly quantified. Fortunately, existing results have summarized various forms of decoherence Lindblad operator in open quantum systems for reference; a general form of (2) can be written as [20] in which is a Wiener incremental process and the random items characterizing the system decoherence are

According to the above system model, in the given interval , a sample of realization can be generated by the following algorithm [20].(1) At the initial time , the initial state of the process is determined by the initial distribution.(2)It is assumed that at time , the normalized state is reached through a quantum jump; then we set . (3)Determine a random waiting time . This can be done, for example, by drawing a random number which is uniformly distributed over the interval and by determining from the equation: First we define the defect of the waiting time distribution by the identity

For , a unique solution can be obtained. If , we set in which case there will be no further jumps. Within the time interval the realization follows the deterministic time evolution: (4) At time (if is finite and ), one of the possible jumps labeled by the index occurs according to Then we select a specific jump of type with probability And then we update the state of system as (5) Repeat steps (1) to (4) until the desired final time is reached, which yields the realization over the whole time interval . Once a sample of realizations , has been generated according to this algorithm, any statistical quantity can be estimated through an appropriate ensemble average.

According to the above framework, an iterative algorithm is described as follows: where

3. Results and Discussion

We consider the process describing the direct photodetection of a driven two-level system, and the piecewise deterministic process is given by the following equation [21, 22]:

The Poisson increments satisfy

The corresponding stochastic Schrodinger equation takes the form where

In the numerical simulations, it is assumed that the atom is in its ground state and the probability of finding the atom in the excited state can be calculated as follows:

From a sample of realizations this probability is estimated by determining the average

An appropriate estimator for the corresponding statistical errors in the determination of from a finite sample of size is given by

According [21, 22], the analytical solution of the process is where

Hence, it is possible to compare the numerical results with the analytical results.

3.1. Simulation Results

A sample of realizations for calculated from (20) is as the dashed line shows in Figure 1 with the following parameters: step size , , and , while the smooth line gives the analytical solution.

Figure 1: The comparison between the numerical solution of and the analytical solution direct photodetection of a driven two-level system with the parameters: , , and .

In order to discuss the performance of the algorithm, we introduced the classic Runge-Kutta iterative algorithm for generating the sample of realizations as follows [23]: where

In this paper, the computer simulation platform is the Intel(R) Core (TM)2 Duo CPU E7500 @ 2.93 GHz under the Windows XP operating system with the numerical calculation software Matlab7. By selecting different simulation step sizes , we can get set of results according to the corresponding with different algorithms:(1)the estimated means of minus the exact values obtained by the analytical solution of for different step sizes and methods, , as Table 1 shows,(2)the normalized CPU time versus the absolute error for the data points, as Table 2 shows.

Table 1: The estimated means of minus the exact values obtained by the analytical solution of , that is, , for different step sizes and methods.
Table 2: The normalized CPU time for different step sizes and methods.

As can be seen from Table 1, when the simulation step size is increased, there is little change in the errors in the proposed algorithm, while the errors using the Runge-Kutta method increase linearly. That is, in a certain step length range, the proposed iterative algorithm can generate more accurate approximation numerical solution than the Runge-Kutta algorithm in comparison to the analytical solution.

At the same time, we can draw a conclusion from Table 2: when obtaining a more accurate numerical solution using the proposed algorithm, the computational time is on the same order of magnitude as the Runge-Kutta method consumes. And it is not difficult to find that, when the step length is gradually reduced, accompanied by improving the accuracy, the computational time of the Runge-Kutta algorithm grows faster.

Taking the above two advantages comparing to the classical Runge-Kutta algorithm, the proposed algorithm can generate a more accurate sample of realizations while paying lower computational costs, which reflects the superiority and practicality of the proposed algorithm.

3.2. Convergence Analysis

Generally speaking, the smaller the step size is, the closer the numerical solution matches the true solution, and consequently the convergence seems to take place. We denote as the true analytical solution and as the numerical approximation. Noting that and are random variables, we can measure the difference using , where represents the expected value. Then a method is said to have strong order of convergence equal to if there exists a constant such that for any fixed and sufficiently small.

The strong order of convergence (29) measures the rate at which the “mean of the error” decays as . A less demanding alternative is to measure the rate of decay of the “error of the means.” This leads to the concept of weak convergence. A method is said to have weak order of convergence equal to if there exists a constant such that for all functions in some class for any fixed and sufficiently small [24].

According the above theory, we consider a numerical algorithm to integrate (7) in the time period . Such an algorithm will generate a discrete time approximate realizations for the exact process at the given time where . Before our discussion, several conventions should be made as follows.(1)In this section, always represents a numerical approximation, while stands for the exact process described by (7). (2)Similarly, we define the discrete time approximation of the density matrix as , whereas is the exact density matrix satisfying the Lindblad equation  (9).(3)We always set a deterministic initial state .

In order to illustrate the numerical convergence of the algorithm, we will compare the Taylor expansion of the exact density matrix which is given by (31) with the generated approximation :

We compare the numerical simulation approximate obtained by the proposed algorithm to the true evolution . The single-step error of a certain numerical scheme may then be expressed through the difference which means that the strategy reproduces the Taylor expansion of including terms of order in .

Thus, it is direct to prove that the integration over a finite time period decreases the convergence order by one since time steps are needed to compute the density at time , that is, with . If (33) can be satisfied, the numerical strategy is defined to be a strategy of order [15]. It should be noted that (33) is a special case of weak convergence of order describing the degree of proximity of the probability distributions of and which is a much weaker criterion comparing to strong convergence defined in (29). In actual applications, one tends to care about this weaker form of convergence especially when considering the approximation of functionals of the stochastic variable.

When investigating the stochastic differential equations and numerical simulation solution process, if the issues are related to the numerical simulation of the stochastic process, the evaluation criteria of the solutions convergence are usually defined as the numerical approximation curve. It must be sufficiently close to the real trace of the evolution. That is to say, the higher the convergence order is, the closer the distribution of the numerical solution is with the analytical solution of the distribution. So it requires that the probability distribution obtained by the simulation iterative algorithm is close enough to the probability distribution determined by the quantum master equation [24].

Performing the error analysis according to (16) one can find it converging with order in contrast to the classical Runge-Kutta method with . Thus, it is really a higher-order strategy in the weaker convergence sense.

4. Conclusions

The decoherence of open quantum systems usually makes the system evolve from the initial pure state to mixed states (in some cases, may also be mixed state into a pure state). Being a powerful tool for investigating the open quantum systems, the quantum master equation can give a quantitative description of the transition, dissipation, and decoherence caused by the interaction between the closed system and the environment. Taking this as the starting point of our research, in order to obtain the evolution of the open quantum systems according to its dynamic characteristics, we used the system simulation method to get the numerical solution to the reduced density operator of a typical open quantum system. And its effectiveness and superiority were verified in comparison with the classical algorithm. Further research includes the control scheme [25, 26] for quantum manipulation based on the characteristics of quantum dynamics.


This work was supported by the National Natural Science Foundation of China (61203062) and the Doctor Start Foundation of Guangdong University of Technology (405105028).


  1. E. Nelson, “Derivation of the Schrodinger equation from Newtonian dynamics,” Physical Review, vol. 150, no. 4, pp. 1079–1085, 1966. View at Publisher · View at Google Scholar
  2. E. Nelson, “Review of stochastic mechanics,” Journal of Physics: Conference Series, vol. 361, no. 1, Article ID 012011, 2012. View at Google Scholar
  3. F. Guerra and L. M. Morato, “Quantization of dynamical systems and stochastic control theory,” Physical Review D, vol. 27, no. 8, pp. 1774–1786, 1983. View at Publisher · View at Google Scholar · View at Scopus
  4. E. A. Carlen, “Conservative diffusion,” Communications in Mathematical Physics, vol. 94, pp. 293–315, 1984. View at Google Scholar
  5. M. I. Loffredo and L. M. Morato, “Lagrangian variational principle in stochastic mechanics: gauge structure and stability,” Journal of Mathematical Physics, vol. 30, no. 2, pp. 354–360, 1989. View at Google Scholar · View at Scopus
  6. M. A. Taylor, J. Janousek, V. Daria, et al., “Biological measurement beyond the quantum limit,” Nature Photonics, vol. 7, pp. 229–233, 2013. View at Publisher · View at Google Scholar
  7. J. I. Cirac and P. Zoller, “Goals and opportunities in quantum simulation,” Nature Physics, vol. 8, pp. 264–266, 2012. View at Google Scholar
  8. J. Cai, A. Retzker, F. Jelezko, et al., “A large-scale quantum simulator on a diamond surface at room temperature,” Nature Physics, vol. 9, pp. 168–173, 2013. View at Publisher · View at Google Scholar
  9. K. C. Young and M. D. Grace, “Simulation of stochastic quantum systems using polynomial chaos expansions,” Physical Review Letters, vol. 110, Article ID 110402, 2013. View at Publisher · View at Google Scholar
  10. K. Ichikawa, M. Fukuda, and A. Tachibana, “Study of simulation method of time evolution in rigged QED,” International Journal of Quantum Chemistry, vol. 113, no. 3, pp. 190–202, 2013. View at Publisher · View at Google Scholar
  11. H. M. Wiseman and H. B. Sun, “Continuous quantum measurement of two coupled quantum dots using a point contact: a quantum trajectory approach,” Physical Review B, vol. 63, no. 12, Article ID 125326, 2001. View at Google Scholar · View at Scopus
  12. T. A. Brun, “A simple model of quantum trajectories,” American Journal of Physics, vol. 70, no. 7, pp. 719–737, 2002. View at Publisher · View at Google Scholar · View at Scopus
  13. O. V. Zhirov and D. L. Shepelyansky, “Dissipative decoherence in the Grover algorithm,” European Physical Journal D, vol. 38, no. 2, pp. 405–408, 2006. View at Publisher · View at Google Scholar · View at Scopus
  14. T. A. Brun and N. Gisin, “Quantum state diffusion and time correlation functions,” 2008, http://arxiv.org/abs/quant-ph/9512013v3.
  15. I. Percival, Quantum State Diffusion, Cambridge University Press, Cambridge, UK, 1998.
  16. A. Puga and B. N. Miller, “Solution of the quantum initial value problem with transparent boundary conditions,” American Journal of Physics, vol. 81, no. 1, pp. 50–55, 2013. View at Google Scholar
  17. V. E. Tarasov, “Fractional diffusion equations for open quantum system,” Nonlinear Dynamics, vol. 71, no. 4, pp. 663–670, 2013. View at Google Scholar
  18. T. Peyronel, O. Firstenberg, Q. Liang, et al., “Quantum nonlinear optics with single photons enabled by strongly interacting atoms,” Nature, vol. 488, pp. 57–60, 2012. View at Google Scholar
  19. A. A. Houck, H. E. Türeci, and J. Koch, “On-chip quantum simulation with superconducting circuits,” Nature Physics, no. 8, pp. 292–299, 2012. View at Google Scholar
  20. H. P. Breuer and F. Petruccione, The Theory of Open Quantum Systems, Oxford University Press, Oxford, UK, 2003.
  21. H. M. Wiseman and G. J. Milburn, “Quantum theory of field-quadrature measurements,” Physical Review A, vol. 47, no. 1, pp. 642–662, 1993. View at Publisher · View at Google Scholar · View at Scopus
  22. H. M. Wiseman and G. J. Milburn, “Interpretation of quantum jump and diffusion processes illustrated on the Bloch sphere,” Physical Review A, vol. 47, no. 3, pp. 1652–1666, 1993. View at Publisher · View at Google Scholar · View at Scopus
  23. J. Roden, W. T. Strunz, and A. Eisfeld, “Non-Markovian quantum state diffusion for absorption spectra of molecular aggregates,” Journal of Chemical Physics, vol. 134, Article ID 034902, 2011. View at Google Scholar
  24. D. J. Higham, “An algorithmic introduction to numerical simulation of stochastic differential equations,” SIAM Review, vol. 43, no. 3, pp. 525–546, 2001. View at Google Scholar · View at Scopus
  25. D. Dong and I. R. Petersen, “Quantum control theory and applications: a survey,” IET Control Theory and Applications, vol. 4, no. 12, pp. 2651–2671, 2010. View at Publisher · View at Google Scholar · View at Scopus
  26. D. Dong and I. R. Peterson, “Sliding mode control of quantum systems,” New Journal of Physics, vol. 11, Article ID 105033, 2009. View at Google Scholar