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

NASA STI/Recon Technical Report N

Stochastic Solution to Quantum Dynamics

Sarah John, Walter John

Published Feb 1, 2003• 0 citations
OpenAlex

Abstract

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...

Related papers

Quantum dynamics as a stochastic process

Sarah John, John W. Wilson

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.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Jan 1, 1994•4 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

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

Simulation of quantum processes using entangled trajectory molecular dynamics

Arnaldo Donoso, Yujun Zheng, Craig C. Martens

In this paper, we describe a new method for simulating quantum processes using classical-like molecular dynamics. The approach is based on solving the quantum Liouville equation in the Wigner representation using ensembles of classical trajectories in phase space. The nonlocality of quantum mechanics is incorporated in the trajectory representation as nonclassical interactions between the members of the ensemble, leading to an entanglement of their evolution. The statistical independence of the individual trajectories making up an ensemble in the classical limit is lost when quantum effects are included, and the entire state of the system must be propagated as a unified whole. We develop the formalism and its numerical implementation, and illustrate its application on two model problems of quantum mechanical tunneling: escape from a metastable well and wave packet penetration of the Eckhart barrier.

The Journal of Chemical Physics•Aug 21, 2003•85 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

Wigner approach to quantum statistical mechanics and quantum generalization molecular dynamics method. Part 2

Vladimir S. Filinov

A new numerical approach for the calculation of quantum time correlation functions in the canonical and micro-canonical ensembles has been developed within the Wigner formulation of quantum statistical mechanics. For electrons in disordered systems of scatterers the numerical results for the tensor of static electrical conductivity have been obtained according to the quantum Kubo formula in the one-electron approximation in the canonical ensemble. The second quantity considered is the flux position time correlation function characterizing the momentum flow for tunnelling particles at the barrier crossing from one potential well to another. The difficulties with studies of exponentially small effects in the temperature averaging for the canonical ensemble has been overcome by a transition to the micro-canonical ensemble with a fixed total energy. The related integral Wigner-Liouville type equations have been obtained and a numerical approach combining both molecular dynamics and Monte Carlo methods for solving this equation has also been developed. Numerical results have been obtained for the momentum flow of a particle tunnelling through the barrier dividing two potential wells in two- and three-dimensional cases.

Molecular Physics•Aug 1, 1996•38 citations
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...
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.
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 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 this paper, we describe a new method for simulating quantum processes using classical-like molecular dynamics. The approach is based on solving the quantum Liouville equation in the Wigner representation using ensembles of classical trajectories in phase space. The nonlocality of quantum mechanics is incorporated in the trajectory representation as nonclassical interactions between the members of the ensemble, leading to an entanglement of their evolution. The statistical independence of the individual trajectories making up an ensemble in the classical limit is lost when quantum effects are included, and the entire state of the system must be propagated as a unified whole. We develop the formalism and its numerical implementation, and illustrate its application on two model problems of quantum mechanical tunneling: escape from a metastable well and wave packet penetration of the Eckhart barrier.
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.
A new numerical approach for the calculation of quantum time correlation functions in the canonical and micro-canonical ensembles has been developed within the Wigner formulation of quantum statistical mechanics. For electrons in disordered systems of scatterers the numerical results for the tensor of static electrical conductivity have been obtained according to the quantum Kubo formula in the one-electron approximation in the canonical ensemble. The second quantity considered is the flux position time correlation function characterizing the momentum flow for tunnelling particles at the barrier crossing from one potential well to another. The difficulties with studies of exponentially small effects in the temperature averaging for the canonical ensemble has been overcome by a transition to the micro-canonical ensemble with a fixed total energy. The related integral Wigner-Liouville type equations have been obtained and a numerical approach combining both molecular dynamics and Monte Carlo methods for solving this equation has also been developed. Numerical results have been obtained for the momentum flow of a particle tunnelling through the barrier dividing two potential wells in two- and three-dimensional cases.