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

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics

Quantum dynamics as a stochastic process

Sarah John, John W. Wilson

Published Jan 1, 1994• 4 citations
Publisher

Abstract

The quantum Liouville equation is solved in the Wigner representation using generalized Monte Carlo techniques. For small increments of time, the solution is represented as a sequential classical evolution in phase space followed by a quantum "jump" distribution in momentum space, with the latter simulated via a stochastic method. Extending the work initiated by John and Remler [Ann Phys. (N.Y.) 180, 152 (1987)] the technique is developed and validated for higher dimensions. Also, an alternative algorithm is developed and applied to study motion of a quantum system in an anharmonic quartic potential well, with significantly improved results.

DOI
10.1103/physreve.49.145

Related papers

Stochastic Solution to Quantum Dynamics

Sarah John, Walter John

The quantum Liouville equation in the Wigner representation is solved numerically by using Monte Carlo methods. For incremental time steps, the propagation is implemented as a classical evolution in phase space modified by a quantum correction. The correction, which is a momentum jump function, is simulated in the quasi-classical approximation via a stochastic process. In this paper the technique, which is developed and validated in two- and three-dimensional momentum space, extends an earlier one-dimensional work. Also, by developing a new algorithm, the application to bound state motion in an anharmonic quartic potential shows better agreement with exact solutions in two-dimensional phase space. I. Introduction Nuclear interaction theory is formulated in the language of quantum mechanics, and hence the development of general methods of solutions to quantum dynamics will provide useful tools for application to a large class of problems in nuclear many-body theory. Different approaches...

NASA STI/Recon Technical Report N•Feb 1, 2003•0 citations

Solution To Quantum Liouville Equation In Phase Space Via A Stochastic Process

Sarah John

A numerical technique for the solution of the Quantum Liouville equation is discussed which is expected to be useful for mixed states describing highly excited many body systems.;The equation is represented in phase-space, where the density operator becomes the Wigner function. The Liouville operator which determines its time evolution, separates into a classical and quantum term. When the Wigner function varies slowly in momentum, the quasi-classical approximation is valid in which one keeps only the first order quantum term in the Liouville equation.;The Wigner function is represented by a sample set of gaussian smeared points. In an iterative time development scheme, the points are propagated first along a classical phase-space trajectory, which is then followed by a quantum development. The latter is represented by a stochastic process incorporated into the numerical procedure as a generalized Markovian process.;This method is numerically tested for an anharmonic quartic potential. The contours of the Wigner function at different time intervals are compared with the purely classical and exact results. It is seen that it provides a definite improvement over the purely classical approximation. Various averages are also compared and sources of error discussed.;In addition it was discovered that intermediate mixed states show a 'quantum focusing' effect, wherein the function peaks to values greater than the original maximum. This quantum-mechanical effect did not appear in the pure case, and disappears because of Liouville's theorem for highly mixed states which are essentially classical. There are other numerical schemes, which treat the points as moving in an effective potential, wherein the points retain their identity. These again by Liouville's theorem cannot show this focusing effect.

•Jan 1, 1986•2 citations

Approximation scheme for the quantum Liouville equation using phase-space distribution functions

Marek Płoszajczak, Mark J. Rhoades-Brown

A general approximation scheme for solving the quantum Liouville equation has been developed and is presented here specifically within the Hartree approximation. This method uses the connection between the density matrix and the Wigner phase-space representations. By use of the Wigner representation the off-diagonal elements of the density matrix are expanded in powers of deviations from the mean momentum value. The time evolution of these deviations is governed by a closed self-consistent set of equations.

Physical Review Letters•Jul 8, 1985•40 citations

Quantum dynamics and Wigner representation of quantum mechanics

Vladimir S. Filinov, Yu. V. Medvedev, V.L. Kamskyi

A new method for solving the Wigner-Liouville equation allowing a treatment of quantum dynamics has been developed. The method combines two widely known approaches: molecular dynamics and Monte Carlo. Numerical results have been obtained for one-, two- and three-dimensional potential wells of realistic type. The quantum dynamics method was tested by a comparison of the numerical results with analytical estimates for such values as the average position x¯(t), average momentum p¯(t), energy Ē(t), position and momentum dispersions, and Heisenberg inequality.

Molecular Physics•Jul 1, 1995•59 citations

Quantum dynamics in Wigner and tomography representations

Vladimir S. Filinov, G. Schubert, Pavel R. Levashov +3 more

In the standard formulation of quantum mechanics states of a system are described by wave functions or density operators . However quantum description of the system can be given in many other ways , for example, in Wigner-Moyal or Feynman formulations . All these representation are equivalent in the the sense of physical results, but the form of presentation of quantum mechanics is different as wave functions, Wigner-Liouville functions or other representations of density operators differ in essential way from each other. In this work we are dealing with Wigner and probability representation of quantum mechanics. The probability representations or tomography representation was recently proposed in terms of marginal distribution functions (MDF). MDF describing the quantum states are positive distribution functions connected with wave functions or density matrices by known integral transformation. Sign conservation of the of the MDF can be valuable in computer simulations to overcome the "sign problem". Important advantage of the probability representation is that the quantum transitions between quantum states are descried by non-negative probabilities. This work is devoted to the development of a new stochastic approaches to numerical solution of the evolution equations for Wigner-Liouviile (WL) and marginal distribution functions. To solve the WL equation we combine both molecular dynamics and Monte Carlo methods and compute traces of the dynamical operators. The results obtained for a system of electrons and random scatterers clearly demonstrate that the many-particle interaction between the electrons leads to an enhancement of the conductivity compared to the case of noninteracting electrons. To obtain evolution of MDF we developed new stochastic approach to solution of the generalized Langevin equation (GLE). GLE is derived from Kolmogorov equations for Green function of evolution equation for MDF. We discuss the basic relations and main ideas of this approach, compare obtained numerical results with results of independent finite-difference calculations for quantum oscillator and quantum particles crossing the finite well and tunneling through the Gaussian barrier.

IEEE International Conference on Plasma Science (ICOPS)•Jun 1, 2008•0 citations

Stochastic path-integral simulation of quantum scattering

Thomas C. Schmidt, K. Möhring

We present a path-integral solution for the exact propagation of the Wigner distribution in phase space, which is an improved version of results obtained [Maslov, Bertrand, Combe, and co-workers, J. Sov. Math. 13, 315 (1980); 19, 55 (1982); Lett. Math. Phys. 7, 327 (1983); Physica 124A, 561 (1984)] and is suitable for Monte Carlo simulations. The trajectories involved are purely stochastic processes exhibiting discontinuous jumps in momentum. The applicability of the method is demonstrated for scattering at a one-dimensional barrier.

Physical Review A•Nov 1, 1993•16 citations
The quantum Liouville equation is solved in the Wigner representation using generalized Monte Carlo techniques. For small increments of time, the solution is represented as a sequential classical evolution in phase space followed by a quantum "jump" distribution in momentum space, with the latter simulated via a stochastic method. Extending the work initiated by John and Remler [Ann Phys. (N.Y.) 180, 152 (1987)] the technique is developed and validated for higher dimensions. Also, an alternative algorithm is developed and applied to study motion of a quantum system in an anharmonic quartic potential well, with significantly improved results.
The quantum Liouville equation in the Wigner representation is solved numerically by using Monte Carlo methods. For incremental time steps, the propagation is implemented as a classical evolution in phase space modified by a quantum correction. The correction, which is a momentum jump function, is simulated in the quasi-classical approximation via a stochastic process. In this paper the technique, which is developed and validated in two- and three-dimensional momentum space, extends an earlier one-dimensional work. Also, by developing a new algorithm, the application to bound state motion in an anharmonic quartic potential shows better agreement with exact solutions in two-dimensional phase space. I. Introduction Nuclear interaction theory is formulated in the language of quantum mechanics, and hence the development of general methods of solutions to quantum dynamics will provide useful tools for application to a large class of problems in nuclear many-body theory. Different approaches...
A numerical technique for the solution of the Quantum Liouville equation is discussed which is expected to be useful for mixed states describing highly excited many body systems.;The equation is represented in phase-space, where the density operator becomes the Wigner function. The Liouville operator which determines its time evolution, separates into a classical and quantum term. When the Wigner function varies slowly in momentum, the quasi-classical approximation is valid in which one keeps only the first order quantum term in the Liouville equation.;The Wigner function is represented by a sample set of gaussian smeared points. In an iterative time development scheme, the points are propagated first along a classical phase-space trajectory, which is then followed by a quantum development. The latter is represented by a stochastic process incorporated into the numerical procedure as a generalized Markovian process.;This method is numerically tested for an anharmonic quartic potential. The contours of the Wigner function at different time intervals are compared with the purely classical and exact results. It is seen that it provides a definite improvement over the purely classical approximation. Various averages are also compared and sources of error discussed.;In addition it was discovered that intermediate mixed states show a 'quantum focusing' effect, wherein the function peaks to values greater than the original maximum. This quantum-mechanical effect did not appear in the pure case, and disappears because of Liouville's theorem for highly mixed states which are essentially classical. There are other numerical schemes, which treat the points as moving in an effective potential, wherein the points retain their identity. These again by Liouville's theorem cannot show this focusing effect.
A general approximation scheme for solving the quantum Liouville equation has been developed and is presented here specifically within the Hartree approximation. This method uses the connection between the density matrix and the Wigner phase-space representations. By use of the Wigner representation the off-diagonal elements of the density matrix are expanded in powers of deviations from the mean momentum value. The time evolution of these deviations is governed by a closed self-consistent set of equations.
A new method for solving the Wigner-Liouville equation allowing a treatment of quantum dynamics has been developed. The method combines two widely known approaches: molecular dynamics and Monte Carlo. Numerical results have been obtained for one-, two- and three-dimensional potential wells of realistic type. The quantum dynamics method was tested by a comparison of the numerical results with analytical estimates for such values as the average position x¯(t), average momentum p¯(t), energy Ē(t), position and momentum dispersions, and Heisenberg inequality.
In the standard formulation of quantum mechanics states of a system are described by wave functions or density operators . However quantum description of the system can be given in many other ways , for example, in Wigner-Moyal or Feynman formulations . All these representation are equivalent in the the sense of physical results, but the form of presentation of quantum mechanics is different as wave functions, Wigner-Liouville functions or other representations of density operators differ in essential way from each other. In this work we are dealing with Wigner and probability representation of quantum mechanics. The probability representations or tomography representation was recently proposed in terms of marginal distribution functions (MDF). MDF describing the quantum states are positive distribution functions connected with wave functions or density matrices by known integral transformation. Sign conservation of the of the MDF can be valuable in computer simulations to overcome the "sign problem". Important advantage of the probability representation is that the quantum transitions between quantum states are descried by non-negative probabilities. This work is devoted to the development of a new stochastic approaches to numerical solution of the evolution equations for Wigner-Liouviile (WL) and marginal distribution functions. To solve the WL equation we combine both molecular dynamics and Monte Carlo methods and compute traces of the dynamical operators. The results obtained for a system of electrons and random scatterers clearly demonstrate that the many-particle interaction between the electrons leads to an enhancement of the conductivity compared to the case of noninteracting electrons. To obtain evolution of MDF we developed new stochastic approach to solution of the generalized Langevin equation (GLE). GLE is derived from Kolmogorov equations for Green function of evolution equation for MDF. We discuss the basic relations and main ideas of this approach, compare obtained numerical results with results of independent finite-difference calculations for quantum oscillator and quantum particles crossing the finite well and tunneling through the Gaussian barrier.
We present a path-integral solution for the exact propagation of the Wigner distribution in phase space, which is an improved version of results obtained [Maslov, Bertrand, Combe, and co-workers, J. Sov. Math. 13, 315 (1980); 19, 55 (1982); Lett. Math. Phys. 7, 327 (1983); Physica 124A, 561 (1984)] and is suitable for Monte Carlo simulations. The trajectories involved are purely stochastic processes exhibiting discontinuous jumps in momentum. The applicability of the method is demonstrated for scattering at a one-dimensional barrier.