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

Journal of Physics A Mathematical and Theoretical

Effective classical Liouville-like evolution equation for the quantum phase-space dynamics

Omar Morandi

Published Jul 23, 2010• 22 citations
Publisher

Abstract

A model for the evolution of a quantum particle in the phase plane is derived. The quasi-distribution function is decomposed into a suitable over-complete basis of coherent states. In this approach, a coarse grain length scale is introduced by using the spread of the coherent states. The obtained evolution system is completely equivalent to the Wigner formulation of the quantum mechanics. It shows some analogy with the von Neumann approach and the particle in the cell method in classical plasma physics. Differing from the usual formulation of the quantum phase space, given in terms of infinite-order partial differential equations, in this model, the evolution equations of the second differential order are expressed by a hierarchy of coupled functions (the first term being the Husimi function). The resulting formulation reveals itself to be particularly close to the classical description of the particles motion. This formal analogy is useful to gain new physical insights and to profit from numerical methods developed for classical systems.

DOI
10.1088/1751-8113/43/36/365302

Related papers

Phase space formulation of the quantum many-body problem

Paul H. Levine

By means of a quantum mechanical phase space distribution function introduced by von Roos, the Schroedinger equation for a non-relativistic system of N identical particles with scalar interactions is transformed into a quantum mechanical generalization of the Liouville equation, thereby formulating the problem in terms of a generalized density in phase space, a quantity of primary interest in most treatments of the corresponding classical system (or plasma). This transformation permits a parallel development of the theories of classical and quantum plasmas and thus allows the quantum many-body problem to be discussed virtually completely in classical terms. In particular, a kinetic theory of quantum plasmas is obtained by deriving the quantum analogue of the BBGKY hierarchy, and applying thereto approximation techniques similar to those of Rostocker and Rosenbluth, and Bogoliubov. The point of departure from similar previous studies based on the Wigner distribution function is that the proper exchange symmetry can be tractably introduced into the formalism. Attention is first focused on the Hartree and Hartree-Fock approximations, in which case the quantum BBGKY system reduces to a simple quantum generalization of the Vlasov equation. This equation is used to study the response of spatially homogeneous systems to weak external forces, and the associated problems of plasmon and spin-wave excitations. It is also used to derive the quantum and exchange corrected equations of inviscid hydrodynamical transport which are then applied to the problem of sound propagation in the degenerate electron gas. The second part of the study is concerned with the theory of the many-electron atom in the Hartree and Hartree-Fock approximations. The relevant quantum Vlasov equations lead naturally to a statistical theory of the atom which reduces to the Thomas-Fermi-Amaldi and Thomas Fermi models (respectively) as ħ → 0. For ħ ≠ 0, the quantum and exchange corrections to these models are simultaneously generated. The quantum hydrodynamical theory developed earlier is used to determine the influence of these corrections on the boundary conditions of the model, and a theory of the compressed atom is consequently obtained. Considered in somewhat less detail are the effects of non-zero temperature, net orbital angular momentum, relativity and correlations, as well as time dependent processes. The final part deals with the problem of the degenerate electron gas with a uniform neutralizing background. Going beyond the Hartree-Fock approximation, the pair functions for particles with parallel and anti-parallel spin are obtained by neglecting three particle correlations. From these functions, a quantum-mechanical collision integral is derived which differs from that obtained by Silin and Guernsey and conjectured by Wyld and Pines in that dynamical exchange effects are included. Also obtained from the pair function is an expression for the correlation which reduces in the high density limit to the result of Gell-Mann and Brueckner. At intermediate densities an additional term appears in the energy due to the screening of the exchange interaction by the dielectric properties of the medium. It is evaluated in the high density limit and found to be -0.151 r_s ln r_s Ryd/electron in marked disagreement with the corresponding value obtained by DuBois.

•Jan 1, 1963•0 citations

The Quantum Mechanical Density Operator And Its Time Evolution: Quantum Dynamics Using The Quantum Liouville Equation

Abraham Nitzan

The starting point of the classical description of motion is the Newton equations that yield a phase space trajectory (rN (t), pN (t)) for a given initial condition (rN (0), pN (0)). Alternatively one may describe classical motion in the framework of the Liouville equation (Section (1.2.2)) that describes the time evolution of the phase space probability density f (rN , pN ; t). For a closed system fully described in terms of a well specified initial condition, the two descriptions are completely equivalent. Probabilistic treatment becomes essential in reduced descriptions that focus on parts of an overall system, as was demonstrated in Sections 5.1–5.3 for equilibrium systems, and in Chapters 7 and 8 that focus on the time evolution of classical systems that interact with their thermal environments. This chapter deals with the analogous quantum mechanical problem. Within the limitations imposed by its nature as expressed, for example, by Heisenbergtype uncertainty principles, the Schrödinger equation is deterministic. Obviously it describes a deterministic evolution of the quantum mechanical wavefunction. The analog of the phase space probability density f (rN , pN ; t) is now the quantum mechanical density operator (often referred to as the “density matrix”), whose time evolution is determined by the quantum Liouville equation. Again, when the system is fully described in terms of a well specified initial wavefunction, the two descriptions are equivalent. The density operator formalism can, however, be carried over to situations where the initial state of the system is not well characterized and/or a reduced description of part of the overall system is desired. Such situations are considered later in this chapter.

Oxford University Press eBooks•Apr 6, 2006•1 citation

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

Perturbation Theory in Terms of a Generalized Phase-Space Quantization Procedure

Omar Morandi, Luigi Barletti, Giovanni Frosali

A new approach to perturbation theory in the quantum phase-space formalism is proposed, in order to devise some approximated description of the quantum phase-space dynamics, which is not directly related to the usual semi-classical approximation. A general class of equivalent quasi-distribution functions based on the Wigner-Moyal formalism is considered and a first-order invariant formulation of the dynamics is obtained. The relationship between the various phase-space representation is expressed in term of a pseudo-differential operator defined by the Moyal product. In particular, our theory is applied to the sub-class of representations obtained by a first order perturbation of the Wigner representation. Finally the connection of our approach with some well established gauge-invariant formulation of the Wigner dynamics in the presence of an external magnetic field is investigated.

Bollettino Della Unione Matematica Italiana•Jan 1, 2011•57 citations

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

Is the Moyal equation for the Wigner function a quantum analogue of the Liouville equation?

E. E. Perepelkin, Б. И. Садовников, Наталья Германовна Иноземцева +2 more

Abstract The Moyal equation describes the evolution of the Wigner function of a quantum system in the phase space. The right-hand side of the equation contains an infinite series with coefficients proportional to powers of the Planck constant. There is an interpretation of the Moyal equation as a quantum analogue of the classical Liouville equation. Indeed, if one uses the notion of the classical passage to the limit as the Planck constant tends to zero, then formally the right-hand side of the Moyal equation tends to zero. As a result, the Moyal equation becomes the classical Liouville equation for the distribution function. In this paper, we show that the right side of the Moyal equation does not explicitly depend on the Planck constant, and all terms of the series can make a significant contribution. The transition between the classical and quantum descriptions is related not to the Planck constant, but to the spatial scale. For a model quantum system with a potential in the form of a «quadratic funnel», an exact 3D solution of the Schrödinger equation is found and the corresponding Wigner function is constructed in the paper. Using trajectory analysis in the phase space, based on the representation of the right-hand side of the Moyal equation, it is shown that on the spatial microscale there is an infinite number of «trajectories» of the particle motion (thereby the concept of a trajectory is indefinite), and when passing to the macroscale, all «trajectories» concentrate around the classical trajectory.

arXiv:2307.16316•Journal of Statistical Mechanics Theory and Experiment•Sep 1, 2023•8 citations
A model for the evolution of a quantum particle in the phase plane is derived. The quasi-distribution function is decomposed into a suitable over-complete basis of coherent states. In this approach, a coarse grain length scale is introduced by using the spread of the coherent states. The obtained evolution system is completely equivalent to the Wigner formulation of the quantum mechanics. It shows some analogy with the von Neumann approach and the particle in the cell method in classical plasma physics. Differing from the usual formulation of the quantum phase space, given in terms of infinite-order partial differential equations, in this model, the evolution equations of the second differential order are expressed by a hierarchy of coupled functions (the first term being the Husimi function). The resulting formulation reveals itself to be particularly close to the classical description of the particles motion. This formal analogy is useful to gain new physical insights and to profit from numerical methods developed for classical systems.
By means of a quantum mechanical phase space distribution function introduced by von Roos, the Schroedinger equation for a non-relativistic system of N identical particles with scalar interactions is transformed into a quantum mechanical generalization of the Liouville equation, thereby formulating the problem in terms of a generalized density in phase space, a quantity of primary interest in most treatments of the corresponding classical system (or plasma). This transformation permits a parallel development of the theories of classical and quantum plasmas and thus allows the quantum many-body problem to be discussed virtually completely in classical terms. In particular, a kinetic theory of quantum plasmas is obtained by deriving the quantum analogue of the BBGKY hierarchy, and applying thereto approximation techniques similar to those of Rostocker and Rosenbluth, and Bogoliubov. The point of departure from similar previous studies based on the Wigner distribution function is that the proper exchange symmetry can be tractably introduced into the formalism. Attention is first focused on the Hartree and Hartree-Fock approximations, in which case the quantum BBGKY system reduces to a simple quantum generalization of the Vlasov equation. This equation is used to study the response of spatially homogeneous systems to weak external forces, and the associated problems of plasmon and spin-wave excitations. It is also used to derive the quantum and exchange corrected equations of inviscid hydrodynamical transport which are then applied to the problem of sound propagation in the degenerate electron gas. The second part of the study is concerned with the theory of the many-electron atom in the Hartree and Hartree-Fock approximations. The relevant quantum Vlasov equations lead naturally to a statistical theory of the atom which reduces to the Thomas-Fermi-Amaldi and Thomas Fermi models (respectively) as ħ → 0. For ħ ≠ 0, the quantum and exchange corrections to these models are simultaneously generated. The quantum hydrodynamical theory developed earlier is used to determine the influence of these corrections on the boundary conditions of the model, and a theory of the compressed atom is consequently obtained. Considered in somewhat less detail are the effects of non-zero temperature, net orbital angular momentum, relativity and correlations, as well as time dependent processes. The final part deals with the problem of the degenerate electron gas with a uniform neutralizing background. Going beyond the Hartree-Fock approximation, the pair functions for particles with parallel and anti-parallel spin are obtained by neglecting three particle correlations. From these functions, a quantum-mechanical collision integral is derived which differs from that obtained by Silin and Guernsey and conjectured by Wyld and Pines in that dynamical exchange effects are included. Also obtained from the pair function is an expression for the correlation which reduces in the high density limit to the result of Gell-Mann and Brueckner. At intermediate densities an additional term appears in the energy due to the screening of the exchange interaction by the dielectric properties of the medium. It is evaluated in the high density limit and found to be -0.151 r_s ln r_s Ryd/electron in marked disagreement with the corresponding value obtained by DuBois.
The starting point of the classical description of motion is the Newton equations that yield a phase space trajectory (rN (t), pN (t)) for a given initial condition (rN (0), pN (0)). Alternatively one may describe classical motion in the framework of the Liouville equation (Section (1.2.2)) that describes the time evolution of the phase space probability density f (rN , pN ; t). For a closed system fully described in terms of a well specified initial condition, the two descriptions are completely equivalent. Probabilistic treatment becomes essential in reduced descriptions that focus on parts of an overall system, as was demonstrated in Sections 5.1–5.3 for equilibrium systems, and in Chapters 7 and 8 that focus on the time evolution of classical systems that interact with their thermal environments. This chapter deals with the analogous quantum mechanical problem. Within the limitations imposed by its nature as expressed, for example, by Heisenbergtype uncertainty principles, the Schrödinger equation is deterministic. Obviously it describes a deterministic evolution of the quantum mechanical wavefunction. The analog of the phase space probability density f (rN , pN ; t) is now the quantum mechanical density operator (often referred to as the “density matrix”), whose time evolution is determined by the quantum Liouville equation. Again, when the system is fully described in terms of a well specified initial wavefunction, the two descriptions are equivalent. The density operator formalism can, however, be carried over to situations where the initial state of the system is not well characterized and/or a reduced description of part of the overall system is desired. Such situations are considered later in this chapter.
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 approach to perturbation theory in the quantum phase-space formalism is proposed, in order to devise some approximated description of the quantum phase-space dynamics, which is not directly related to the usual semi-classical approximation. A general class of equivalent quasi-distribution functions based on the Wigner-Moyal formalism is considered and a first-order invariant formulation of the dynamics is obtained. The relationship between the various phase-space representation is expressed in term of a pseudo-differential operator defined by the Moyal product. In particular, our theory is applied to the sub-class of representations obtained by a first order perturbation of the Wigner representation. Finally the connection of our approach with some well established gauge-invariant formulation of the Wigner dynamics in the presence of an external magnetic field is investigated.
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...
Abstract The Moyal equation describes the evolution of the Wigner function of a quantum system in the phase space. The right-hand side of the equation contains an infinite series with coefficients proportional to powers of the Planck constant. There is an interpretation of the Moyal equation as a quantum analogue of the classical Liouville equation. Indeed, if one uses the notion of the classical passage to the limit as the Planck constant tends to zero, then formally the right-hand side of the Moyal equation tends to zero. As a result, the Moyal equation becomes the classical Liouville equation for the distribution function. In this paper, we show that the right side of the Moyal equation does not explicitly depend on the Planck constant, and all terms of the series can make a significant contribution. The transition between the classical and quantum descriptions is related not to the Planck constant, but to the spatial scale. For a model quantum system with a potential in the form of a «quadratic funnel», an exact 3D solution of the Schrödinger equation is found and the corresponding Wigner function is constructed in the paper. Using trajectory analysis in the phase space, based on the representation of the right-hand side of the Moyal equation, it is shown that on the spatial microscale there is an infinite number of «trajectories» of the particle motion (thereby the concept of a trajectory is indefinite), and when passing to the macroscale, all «trajectories» concentrate around the classical trajectory.