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

The Journal of Chemical Physics

Initial-value semiclassical propagators for the Wigner phase space representation: Formulation based on the interpretation of the Moyal equation as a Schrödinger equation

Shin-ichi Koda

Published Dec 24, 2015• 19 citations
Publisher

Abstract

We formulate various semiclassical propagators for the Wigner phase space representation from a unified point of view. As is shown in several studies, the Moyal equation, which is an equation of motion for the Wigner distribution function, can be regarded as the Schrödinger equation of an extended Hamiltonian system where its "position" and "momentum" correspond to the middle point of two points of the original phase space and the difference between them, respectively. Then we show that various phase-space semiclassical propagators can be formulated just by applying existing semiclassical propagators to the extended system. As a result, a phase space version of the Van Vleck propagator, the initial-value Van Vleck propagator, the Herman-Kluk propagator, and the thawed Gaussian approximation are obtained. In addition, we numerically compare the initial-value phase-space Van Vleck propagator, the phase-space Herman-Kluk propagator, and the classical mechanical propagation as approximation methods for the time propagation of the Wigner distribution function in terms of both accuracy and convergence speed. As a result, we find that the convergence speed of the Van Vleck propagator is far slower than others as is the case of the Hilbert space, and the Herman-Kluk propagator keeps its accuracy for a long period compared with the classical mechanical propagation while the convergence speed of the latter is faster than the former.

DOI
10.1063/1.4938235

Related papers

Semiclassical propagation of Wigner functions

T. Dittrich, E. A. Gómez, Leonardo A. Pachón

We present a comprehensive study of semiclassical phase-space propagation in the Wigner representation, emphasizing numerical applications, in particular as an initial-value representation. Two semiclassical approximation schemes are discussed. The propagator of the Wigner function based on van Vleck's approximation replaces the Liouville propagator by a quantum spot with an oscillatory pattern reflecting the interference between pairs of classical trajectories. Employing phase-space path integration instead, caustics in the quantum spot are resolved in terms of Airy functions. We apply both to two benchmark models of nonlinear molecular potentials, the Morse oscillator and the quartic double well, to test them in standard tasks such as computing autocorrelation functions and propagating coherent states. The performance of semiclassical Wigner propagation is very good even in the presence of marked quantum effects, e.g., in coherent tunneling and in propagating Schrodinger cat states, and of classical chaos in four-dimensional phase space. We suggest options for an effective numerical implementation of our method and for integrating it in Monte-Carlo-Metropolis algorithms suitable for high-dimensional systems.

arXiv:0911.3871•The Journal of Chemical Physics•Jun 2, 2010•36 citations

General initial value form of the semiclassical propagator

Bambi Hu, Quanlin Jie, Baowen Li +1 more

We show a general initial value form of the semiclassical propagator. Similar to cellular dynamics, this formulation involves only the nearby orbits approximation: the evolution of nearby orbits is approximated by linearized dynamics. This phase space smearing formulation keeps the accuracy of the original Van Vleck--Gutzwiller propagator. As an illustration, we present a simple initial value form of the semiclassical propagator. It is nonsingular everywhere and is efficient for numeric implementation.

Physical Review A•Mar 16, 2001•7 citations

Comment on Semiclassical approximations in phase space with coherent states

Frank Großmann, Michael F. Herman

The semiclassical Herman–Kluk initial value approximation of the quantum-mechanical propagator, which is heavily used in dynamical studies of atomic and molecular systems, takes into account nonlinearity of the underlying classical dynamics by using multiple classical initial value solutions. In applications to the propagation of Gaussian wavepackets non-Gaussian distortions in the course of time and therefore a realistic description of the quantum dynamics are possible. This is in contrast to the single-trajectory approximations investigated in a recent paper by Baranger et al , one of which was erroneously termed Herman–Kluk approximation.

Journal of Physics A Mathematical and General•Oct 21, 2002•38 citations

Semiclassical dynamics with quantum trajectories: Formulation and comparison with the semiclassical initial value representation propagator

Sophya V. Garashchuk, Vitaly A. Rassolov

We present a time-dependent semiclassical method based on quantum trajectories. Quantum-mechanical effects are described via the quantum potential computed from the wave function density approximated as a linear combination of Gaussian fitting functions. The number of the fitting functions determines the accuracy of the approximate quantum potential (AQP). One Gaussian fit reproduces time-evolution of a Gaussian wave packet in a parabolic potential. The limit of the large number of fitting Gaussians and trajectories gives the full quantum-mechanical result. The method is systematically improvable from classical to fully quantum. The fitting procedure is implemented as a gradient minimization. We also compare AQP method to the widely used semiclassical propagator of Herman and Kluk by computing energy-resolved transmission probabilities for the Eckart barrier from the wave packet time-correlation functions. We find the results obtained with the Herman–Kluk propagator to be essentially equivalent to those of AQP method with a one-Gaussian density fit for several barrier widths.

The Journal of Chemical Physics•Jan 30, 2003•54 citations

Generalization of the coherent-state path integrals and systematic derivation of semiclassical propagators

Shin-ichi Koda, Kazuo Takatsuka

The coherent path integral is generalized such that the identity operator represented in a complete (actually overcomplete) set of the coherent states with the ``time-variable'' exponents are inserted between two consecutive short-time propagators. Since such a complete set of any given exponent can constitute the identity operator, the exponent may be varied from time to time without loss of generality as long as it is set common to all the Gaussians. However, a finite truncation of the coherent state expansion should result in different values of the propagator depending on the choice of the exponents. Furthermore, approximation methodology to treat with the exact propagator can also depend on this choice, and thereby many different semiclassical propagators may emerge from these combinations. Indeed, we show that the well-known semiclassical propagators such as those of Van Vleck, Herman-Kluk, Heller's thawed Gaussian, and many others can be derived in a systematic manner, which enables one to comprehend these semiclassical propagators from a unified point of view. We are particularly interested in our generalized form of the Herman-Kluk propagator, since the relative accuracy of this propagator has been well established by Kay, and since, nevertheless, its derivation was not necessarily clear. Thus our generalized Herman-Kluk propagator replaces the classical Hamiltonian with a Gaussian averaged quantum Hamiltonian, generating non-Newtonian trajectories. We perform a numerical test to assess the quality of such a family of generalized Herman-Kluk propagators and find that the original Herman-Kluk gives an accurate result. The reason why this has come about is also discussed.

Physical Review A•Mar 21, 2011•6 citations

Semiclassical approximations in phase space with coherent states

Michel Baranger, Marcus A. M. de Aguiar, F Keck +2 more

We present a complete derivation of the semiclassical limit of the coherent-state propagator in one dimension, starting from path integrals in phase space. We show that the arbitrariness in the path integral representation, which follows from the overcompleteness of the coherent states, results in many different semiclassical limits. We explicitly derive two possible semiclassical formulae for the propagator, we suggest a third one, and we discuss their relationships. We also derive an initial-value representation for the semiclassical propagator, based on an initial Gaussian wavepacket. It turns out to be related to, but different from, Heller's thawed Gaussian approximation. It is very different from the Herman-Kluk formula, which is not a correct semiclassical limit. We point out errors in two derivations of the latter. Finally we show how the semiclassical coherent-state propagators lead to WKB-type quantization rules and to approximations for the Husimi distributions of stationary states.

arXiv:quant-ph/0105153•Journal of Physics A Mathematical and General•Sep 3, 2001•176 citations
We formulate various semiclassical propagators for the Wigner phase space representation from a unified point of view. As is shown in several studies, the Moyal equation, which is an equation of motion for the Wigner distribution function, can be regarded as the Schrödinger equation of an extended Hamiltonian system where its "position" and "momentum" correspond to the middle point of two points of the original phase space and the difference between them, respectively. Then we show that various phase-space semiclassical propagators can be formulated just by applying existing semiclassical propagators to the extended system. As a result, a phase space version of the Van Vleck propagator, the initial-value Van Vleck propagator, the Herman-Kluk propagator, and the thawed Gaussian approximation are obtained. In addition, we numerically compare the initial-value phase-space Van Vleck propagator, the phase-space Herman-Kluk propagator, and the classical mechanical propagation as approximation methods for the time propagation of the Wigner distribution function in terms of both accuracy and convergence speed. As a result, we find that the convergence speed of the Van Vleck propagator is far slower than others as is the case of the Hilbert space, and the Herman-Kluk propagator keeps its accuracy for a long period compared with the classical mechanical propagation while the convergence speed of the latter is faster than the former.
We present a comprehensive study of semiclassical phase-space propagation in the Wigner representation, emphasizing numerical applications, in particular as an initial-value representation. Two semiclassical approximation schemes are discussed. The propagator of the Wigner function based on van Vleck's approximation replaces the Liouville propagator by a quantum spot with an oscillatory pattern reflecting the interference between pairs of classical trajectories. Employing phase-space path integration instead, caustics in the quantum spot are resolved in terms of Airy functions. We apply both to two benchmark models of nonlinear molecular potentials, the Morse oscillator and the quartic double well, to test them in standard tasks such as computing autocorrelation functions and propagating coherent states. The performance of semiclassical Wigner propagation is very good even in the presence of marked quantum effects, e.g., in coherent tunneling and in propagating Schrodinger cat states, and of classical chaos in four-dimensional phase space. We suggest options for an effective numerical implementation of our method and for integrating it in Monte-Carlo-Metropolis algorithms suitable for high-dimensional systems.
We show a general initial value form of the semiclassical propagator. Similar to cellular dynamics, this formulation involves only the nearby orbits approximation: the evolution of nearby orbits is approximated by linearized dynamics. This phase space smearing formulation keeps the accuracy of the original Van Vleck--Gutzwiller propagator. As an illustration, we present a simple initial value form of the semiclassical propagator. It is nonsingular everywhere and is efficient for numeric implementation.
The semiclassical Herman–Kluk initial value approximation of the quantum-mechanical propagator, which is heavily used in dynamical studies of atomic and molecular systems, takes into account nonlinearity of the underlying classical dynamics by using multiple classical initial value solutions. In applications to the propagation of Gaussian wavepackets non-Gaussian distortions in the course of time and therefore a realistic description of the quantum dynamics are possible. This is in contrast to the single-trajectory approximations investigated in a recent paper by Baranger et al , one of which was erroneously termed Herman–Kluk approximation.
We present a time-dependent semiclassical method based on quantum trajectories. Quantum-mechanical effects are described via the quantum potential computed from the wave function density approximated as a linear combination of Gaussian fitting functions. The number of the fitting functions determines the accuracy of the approximate quantum potential (AQP). One Gaussian fit reproduces time-evolution of a Gaussian wave packet in a parabolic potential. The limit of the large number of fitting Gaussians and trajectories gives the full quantum-mechanical result. The method is systematically improvable from classical to fully quantum. The fitting procedure is implemented as a gradient minimization. We also compare AQP method to the widely used semiclassical propagator of Herman and Kluk by computing energy-resolved transmission probabilities for the Eckart barrier from the wave packet time-correlation functions. We find the results obtained with the Herman–Kluk propagator to be essentially equivalent to those of AQP method with a one-Gaussian density fit for several barrier widths.
The coherent path integral is generalized such that the identity operator represented in a complete (actually overcomplete) set of the coherent states with the ``time-variable'' exponents are inserted between two consecutive short-time propagators. Since such a complete set of any given exponent can constitute the identity operator, the exponent may be varied from time to time without loss of generality as long as it is set common to all the Gaussians. However, a finite truncation of the coherent state expansion should result in different values of the propagator depending on the choice of the exponents. Furthermore, approximation methodology to treat with the exact propagator can also depend on this choice, and thereby many different semiclassical propagators may emerge from these combinations. Indeed, we show that the well-known semiclassical propagators such as those of Van Vleck, Herman-Kluk, Heller's thawed Gaussian, and many others can be derived in a systematic manner, which enables one to comprehend these semiclassical propagators from a unified point of view. We are particularly interested in our generalized form of the Herman-Kluk propagator, since the relative accuracy of this propagator has been well established by Kay, and since, nevertheless, its derivation was not necessarily clear. Thus our generalized Herman-Kluk propagator replaces the classical Hamiltonian with a Gaussian averaged quantum Hamiltonian, generating non-Newtonian trajectories. We perform a numerical test to assess the quality of such a family of generalized Herman-Kluk propagators and find that the original Herman-Kluk gives an accurate result. The reason why this has come about is also discussed.
We present a complete derivation of the semiclassical limit of the coherent-state propagator in one dimension, starting from path integrals in phase space. We show that the arbitrariness in the path integral representation, which follows from the overcompleteness of the coherent states, results in many different semiclassical limits. We explicitly derive two possible semiclassical formulae for the propagator, we suggest a third one, and we discuss their relationships. We also derive an initial-value representation for the semiclassical propagator, based on an initial Gaussian wavepacket. It turns out to be related to, but different from, Heller's thawed Gaussian approximation. It is very different from the Herman-Kluk formula, which is not a correct semiclassical limit. We point out errors in two derivations of the latter. Finally we show how the semiclassical coherent-state propagators lead to WKB-type quantization rules and to approximations for the Husimi distributions of stationary states.