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

Quantum Dynamics With Classical and Quantum Trajectories

Jian Liu

Published Jan 1, 2005• 0 citations
OpenAlex

Abstract

The forward-backward semiclassical dynamics (FBSD) approximation provides a rigorous and practically useful methodology for simulating the dynamics of condensed phase systems that exhibit quantum mechanical effects. It is shown that FBSD autocorrelation functions of momentum and position-dependent operators satisfy rigorously the time symmetry relation of the quantum mechanical expression. The detailed balance condition follows from this symmetry in the special case where the density is given by the Boltzmann operator. The characteristics of phase space distributions in FBSD calculations are investigated. By virtue of Liouville's theorem and energy conservation, the volume of the negative regions is rigorously conserved and the energy distribution is invariant during time evolution. Thus, while the phase space density is not invariant under FBSD, exhibiting a weak time dependence mostly in its wings, it retains its quantum mechanical characteristics and does not revert to a classical Boltzmann distribution at long times. Illustrative applications on liquid neon are presented.

Related papers

Theoretical methods for the study of quantum dynamics

Jonathan E. Chen

Several theoretical methods for the computation of quantum dynamical quantities are formulated, implemented, and applied, in the overall context of overcoming the exponential scaling property of quantum mechanics. The convergence of the approximate forward-backward semiclassical dynamics (FBSD) method has been accelerated by working with only the imaginary part of the real-time position and momentum time correlation functions. Simple manipulations allow the computationally favorable limit of a single factorization (“bead”) of the Boltzmann operator to yield greatly improved results. Numerical examples for a one-dimensional anharmonic oscillator and Lennard-Jones fluids are presented. Bohm's trajectory formulation of quantum mechanics is an exact mapping of the time-dependent Schrödinger equation onto classical fluid dynamics, plus one extra term responsible for all non-classical effects, the dreaded quantum potential. A computational procedure has been introduced to stably solve Bohm's hydrodynamic equations of motion in the forward-backward quantum dynamics (FBQD) formulation of time correlation functions, and applied to one-dimensional bound systems at zero and low temperature. A key step in the formulation of this procedure was the use of Hamilton's law of varying action, a generalization of Hamilton's principle of stationary action applicable to initial-value problems. This has resulted in the first reported synthetic computation of non-stationary quantum trajectories for a fully anharmonic oscillator. Success in one dimension has laid the groundwork for future extension to a rigorous quantum-classical scheme, which is briefly discussed. Other unpublished methods for obtaining time correlation functions and propagators in quantum dynamics are also discussed.

•May 19, 2010•0 citations

Forward–backward semiclassical dynamics in the interaction representation

Jiushu Shao, Nancy Makri

The forward–backward semiclassical dynamics methodology [J. Phys. Chem. 103, 7753, 9479 (1999)] is reformulated in the interaction representation. The new version of the method allows for a fully quantum mechanical description of a low-dimensional subsystem of interest, along with a semiclassical forward–backward treatment of the solvent coordinates and their coupling to the reference subsystem. Application to the long-time tunneling dynamics in a symmetric double-well system coupled to a harmonic bath shows that the interaction FBSD is capable of capturing quantitatively the tunneling and decoherence effects induced by weakly dissipative environments.

The Journal of Chemical Physics•Sep 1, 2000•36 citations

Nonadiabatic dynamics in open quantum-classical systems: Forward-backward trajectory solution

Chang‐Yu Hsieh, Raymond Kapral

A new approximate solution to the quantum-classical Liouville equation is derived starting from the formal solution of this equation in forward-backward form. The time evolution of a mixed quantum-classical system described by this equation is obtained in a coherent state basis using the mapping representation, which expresses N quantum degrees of freedom in a 2N-dimensional phase space. The solution yields a simple dynamics in which a set of N coherent state coordinates evolves in forward and backward trajectories, while the bath coordinates evolve under the influence of the mean potential that depends on these forward and backward trajectories. It is shown that the solution satisfies the differential form of the quantum-classical Liouville equation exactly. Relations to other mixed quantum-classical and semi-classical schemes are discussed.

arXiv:1205.1283•The Journal of Chemical Physics•Jul 18, 2012•112 citations

A quantum mechanical representation in phase space

Go. Torres-Vega, John H. Frederick

A quantum mechanical representation suitable for studying the time evolution of quantum densities in phase space is proposed and examined in detail. This representation on ℒ2 (2) phase space is based on definitions of the operators P̂ and Q̂ in phase space that satisfy various correspondences for the Liouville equation in classical and quantum phase space, as well as quantum position and momentum ℒ2 (1) spaces. The definitions presented here, P̂=p/2−iℏ∂/∂q and Q̂=q/2+iℏ∂/∂p, are related to definitions that have been recently proposed [J. Chem. Phys. 93, 8862 (1990)]. The resulting quantum phase space representation shares many of the mathematical properties of usual representations in coordinate and momentum spaces. Within this representation, time evolution equations for complex-valued functions (wave functions) and their square magnitudes (distribution functions) are derived, and it is shown that the coordinate and momentum space time evolution equations can be recovered by a simple Fourier projection. The phase space quantum probability conservation equation obtained is a good illustration of the quantization rule that requires one to replace the classical Poisson bracket between the Hamiltonian and the probability density with the quantum commutator between the corresponding operators. The possible classical analogs to quantum probabilities densities are also considered and some of the present results are illustrated for the dynamics of the coherent state.

The Journal of Chemical Physics•Feb 15, 1993•95 citations

An approach for generating trajectory-based dynamics which conserves the canonical distribution in the phase space formulation of quantum mechanics. II. Thermal correlation functions

Jian Liu, William Hughes Miller

We show the exact expression of the quantum mechanical time correlation function in the phase space formulation of quantum mechanics. The trajectory-based dynamics that conserves the quantum canonical distribution-equilibrium Liouville dynamics (ELD) proposed in Paper I is then used to approximately evaluate the exact expression. It gives exact thermal correlation functions (of even nonlinear operators, i.e., nonlinear functions of position or momentum operators) in the classical, high temperature, and harmonic limits. Various methods have been presented for the implementation of ELD. Numerical tests of the ELD approach in the Wigner or Husimi phase space have been made for a harmonic oscillator and two strongly anharmonic model problems, for each potential autocorrelation functions of both linear and nonlinear operators have been calculated. It suggests ELD can be a potentially useful approach for describing quantum effects for complex systems in condense phase.

The Journal of Chemical Physics•Mar 8, 2011•38 citations

Nonadiabatic dynamics in open quantum-classical systems: Forward-backward trajectory solution

Chang‐Yu Hsieh, Raymond Kapral

A new approximate solution to the quantum-classical Liouville equation is derived starting from the formal solution of this equation in forward-backward form. The time evolution of a mixed quantum-classical system described by this equation is obtained in a coherent state basis using the mapping representation, which expresses N quantum degrees of freedom in a 2N-dimensional phase space. The solution yields a simple non-Hamiltonian dynamics in which a set of N coherent state coordinates evolve in forward and backward trajectories while the bath coordinates evolve under the influence of the mean potential that depends on these forward and backward tra-jectories. It is shown that the solution satisfies the differential form of the quantum-classical Liouville equation exactly. Relations to other mixed quantum-classical and semi-classical schemes are discussed. 1 ar

•Aug 17, 2016•70 citations
The forward-backward semiclassical dynamics (FBSD) approximation provides a rigorous and practically useful methodology for simulating the dynamics of condensed phase systems that exhibit quantum mechanical effects. It is shown that FBSD autocorrelation functions of momentum and position-dependent operators satisfy rigorously the time symmetry relation of the quantum mechanical expression. The detailed balance condition follows from this symmetry in the special case where the density is given by the Boltzmann operator. The characteristics of phase space distributions in FBSD calculations are investigated. By virtue of Liouville's theorem and energy conservation, the volume of the negative regions is rigorously conserved and the energy distribution is invariant during time evolution. Thus, while the phase space density is not invariant under FBSD, exhibiting a weak time dependence mostly in its wings, it retains its quantum mechanical characteristics and does not revert to a classical Boltzmann distribution at long times. Illustrative applications on liquid neon are presented.
Several theoretical methods for the computation of quantum dynamical quantities are formulated, implemented, and applied, in the overall context of overcoming the exponential scaling property of quantum mechanics. The convergence of the approximate forward-backward semiclassical dynamics (FBSD) method has been accelerated by working with only the imaginary part of the real-time position and momentum time correlation functions. Simple manipulations allow the computationally favorable limit of a single factorization (“bead”) of the Boltzmann operator to yield greatly improved results. Numerical examples for a one-dimensional anharmonic oscillator and Lennard-Jones fluids are presented. Bohm's trajectory formulation of quantum mechanics is an exact mapping of the time-dependent Schrödinger equation onto classical fluid dynamics, plus one extra term responsible for all non-classical effects, the dreaded quantum potential. A computational procedure has been introduced to stably solve Bohm's hydrodynamic equations of motion in the forward-backward quantum dynamics (FBQD) formulation of time correlation functions, and applied to one-dimensional bound systems at zero and low temperature. A key step in the formulation of this procedure was the use of Hamilton's law of varying action, a generalization of Hamilton's principle of stationary action applicable to initial-value problems. This has resulted in the first reported synthetic computation of non-stationary quantum trajectories for a fully anharmonic oscillator. Success in one dimension has laid the groundwork for future extension to a rigorous quantum-classical scheme, which is briefly discussed. Other unpublished methods for obtaining time correlation functions and propagators in quantum dynamics are also discussed.
The forward–backward semiclassical dynamics methodology [J. Phys. Chem. 103, 7753, 9479 (1999)] is reformulated in the interaction representation. The new version of the method allows for a fully quantum mechanical description of a low-dimensional subsystem of interest, along with a semiclassical forward–backward treatment of the solvent coordinates and their coupling to the reference subsystem. Application to the long-time tunneling dynamics in a symmetric double-well system coupled to a harmonic bath shows that the interaction FBSD is capable of capturing quantitatively the tunneling and decoherence effects induced by weakly dissipative environments.
A new approximate solution to the quantum-classical Liouville equation is derived starting from the formal solution of this equation in forward-backward form. The time evolution of a mixed quantum-classical system described by this equation is obtained in a coherent state basis using the mapping representation, which expresses N quantum degrees of freedom in a 2N-dimensional phase space. The solution yields a simple dynamics in which a set of N coherent state coordinates evolves in forward and backward trajectories, while the bath coordinates evolve under the influence of the mean potential that depends on these forward and backward trajectories. It is shown that the solution satisfies the differential form of the quantum-classical Liouville equation exactly. Relations to other mixed quantum-classical and semi-classical schemes are discussed.
A quantum mechanical representation suitable for studying the time evolution of quantum densities in phase space is proposed and examined in detail. This representation on ℒ2 (2) phase space is based on definitions of the operators P̂ and Q̂ in phase space that satisfy various correspondences for the Liouville equation in classical and quantum phase space, as well as quantum position and momentum ℒ2 (1) spaces. The definitions presented here, P̂=p/2−iℏ∂/∂q and Q̂=q/2+iℏ∂/∂p, are related to definitions that have been recently proposed [J. Chem. Phys. 93, 8862 (1990)]. The resulting quantum phase space representation shares many of the mathematical properties of usual representations in coordinate and momentum spaces. Within this representation, time evolution equations for complex-valued functions (wave functions) and their square magnitudes (distribution functions) are derived, and it is shown that the coordinate and momentum space time evolution equations can be recovered by a simple Fourier projection. The phase space quantum probability conservation equation obtained is a good illustration of the quantization rule that requires one to replace the classical Poisson bracket between the Hamiltonian and the probability density with the quantum commutator between the corresponding operators. The possible classical analogs to quantum probabilities densities are also considered and some of the present results are illustrated for the dynamics of the coherent state.
We show the exact expression of the quantum mechanical time correlation function in the phase space formulation of quantum mechanics. The trajectory-based dynamics that conserves the quantum canonical distribution-equilibrium Liouville dynamics (ELD) proposed in Paper I is then used to approximately evaluate the exact expression. It gives exact thermal correlation functions (of even nonlinear operators, i.e., nonlinear functions of position or momentum operators) in the classical, high temperature, and harmonic limits. Various methods have been presented for the implementation of ELD. Numerical tests of the ELD approach in the Wigner or Husimi phase space have been made for a harmonic oscillator and two strongly anharmonic model problems, for each potential autocorrelation functions of both linear and nonlinear operators have been calculated. It suggests ELD can be a potentially useful approach for describing quantum effects for complex systems in condense phase.
A new approximate solution to the quantum-classical Liouville equation is derived starting from the formal solution of this equation in forward-backward form. The time evolution of a mixed quantum-classical system described by this equation is obtained in a coherent state basis using the mapping representation, which expresses N quantum degrees of freedom in a 2N-dimensional phase space. The solution yields a simple non-Hamiltonian dynamics in which a set of N coherent state coordinates evolve in forward and backward trajectories while the bath coordinates evolve under the influence of the mean potential that depends on these forward and backward tra-jectories. It is shown that the solution satisfies the differential form of the quantum-classical Liouville equation exactly. Relations to other mixed quantum-classical and semi-classical schemes are discussed. 1 ar