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

arXiv:2006.13902

Complexified phase spaces, initial value representations, and the accuracy of semiclassical propagation

Gabriel M. Lando

Published Jun 24, 2020• 0 citations
arXivHTMLPDF

Abstract

Using phase-space complexification, an Initial Value Representation (IVR) for the semiclassical propagator in position space is obtained as a composition of inverse Segal-Bargmann (S-B) transforms of the semiclassical coherent state propagator. The result is shown to be free of caustic singularities and identical to the Herman-Kluk (H-K) propagator, found ubiquitously in physical and chemical applications. We contrast the theoretical aspects of this particular IVR with the van Vleck-Gutzwiller (vV-G) propagator and one of its IVRs, often employed in order to evade the non-linear "root-search" for trajectories required by vV-G. We demonstrate that bypassing the root-search comes at the price of serious numerical instability for all IVRs except the H-K propagator. We back up our theoretical arguments with comprehensive numerical calculations performed using the homogeneous Kerr system, about which we also unveil some unexpected new phenomena, namely: (1) the observation of a clear mark of half the Ehrenfest's time in semiclassical dynamics; and (2) the accumulation of trajectories around caustics as a function of increasing time (dubbed "caustic stickiness"). We expect these phenomena to be more general than for the Kerr system alone.

DOI
10.48550/arxiv.2006.13902

Related papers

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

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.

The Journal of Chemical Physics•Dec 24, 2015•19 citations

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

Semiclassical propagation: Phase indices and the initial-value formalism

Giuseppe Campolieti, Paul W. Brumer

Initial-value phase-space integral representations for a time-dependent propagator are obtained in the coordinate and momentum representations. To do so we first derive nonuniform semiclassical propagators for the various representations, obtaining the global-time semiclassical phase indices (Maslov indices) due to caustics. Results include readily implementable general phase index formulas for any type of caustic, including cases where the Morse index theorem is inapplicable. The method of obtaining the indices is general and based simply on concatenating uniform short-time propagators which also gives rise to alternative path-integral forms. Initial-value integral representations are then derived by introducing a method of extending short-time initial-value propagator formulas to global times via a simple stationary-phase asympotoic-equivalence approach. The integrals reduce to the nonuniform semiclassical propagators within the stationary-phase approximation, are uniform about caustics, and have integrand phases which properly account for the global-time phases in terms of appropriate Maslov indices. The initial-value integrals are also consistently derived via a canonical mapping procedure on the coordinate-space path integral. Initial-value integrals for time-dependent wave-function propgation are also given. Evaluation of the initial-value integral expressions do not require trajectory root searches for propagation.

Physical Review A•Aug 1, 1994•132 citations

An Initial Value Representation with Complex Trajectories

Marcus A. M. de Aguiar, S. A. Vitiello, Adriano Grigolo

We present an Initial Value Representation for the semiclassical coherent state propagator based on complex trajectories. We map the complex phase space into a real phase space with twice as many dimensions and use a simple procedure to automatically eliminate non-contributing trajectories. The resulting semiclassical formulas do not show divergences due to caustics and provide accurate results.

arXiv:0908.1875•Aug 13, 2009•0 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

Semiclassical coherent-states propagator

Alejandro M. F. Rivas

In this work, we derived a semiclassical approximation for the matrix elements of a quantum propagator in coherent states (CS) basis that avoids complex trajectories; it only involves real ones. For that purpose, we used the symplectically invariant semiclassical Weyl propagator obtained by performing a stationary phase approximation (SPA) for the path integral in the Weyl representation. After that, for the transformation to CS representation SPA is avoided; instead a quadratic expansion of the complex exponent is used. This procedure also allows us to express the semiclassical CS propagator uniquely in terms of the classical evolution of the initial point without the need of any root search typical of van Vleck--Gutzwiller-based propagators. For the case of chaotic Hamiltonian systems, the explicit time dependence of the CS propagator has been obtained. The comparison with a realistic chaotic system that derives from a quadratic Hamiltonian, the cat map, reveals that the expression here derived is exact up to quadratic Hamiltonian systems.

arXiv:1207.0012•Physical Review A•Jul 8, 2013•2 citations
Using phase-space complexification, an Initial Value Representation (IVR) for the semiclassical propagator in position space is obtained as a composition of inverse Segal-Bargmann (S-B) transforms of the semiclassical coherent state propagator. The result is shown to be free of caustic singularities and identical to the Herman-Kluk (H-K) propagator, found ubiquitously in physical and chemical applications. We contrast the theoretical aspects of this particular IVR with the van Vleck-Gutzwiller (vV-G) propagator and one of its IVRs, often employed in order to evade the non-linear "root-search" for trajectories required by vV-G. We demonstrate that bypassing the root-search comes at the price of serious numerical instability for all IVRs except the H-K propagator. We back up our theoretical arguments with comprehensive numerical calculations performed using the homogeneous Kerr system, about which we also unveil some unexpected new phenomena, namely: (1) the observation of a clear mark of half the Ehrenfest's time in semiclassical dynamics; and (2) the accumulation of trajectories around caustics as a function of increasing time (dubbed "caustic stickiness"). We expect these phenomena to be more general than for the Kerr system alone.
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.
Initial-value phase-space integral representations for a time-dependent propagator are obtained in the coordinate and momentum representations. To do so we first derive nonuniform semiclassical propagators for the various representations, obtaining the global-time semiclassical phase indices (Maslov indices) due to caustics. Results include readily implementable general phase index formulas for any type of caustic, including cases where the Morse index theorem is inapplicable. The method of obtaining the indices is general and based simply on concatenating uniform short-time propagators which also gives rise to alternative path-integral forms. Initial-value integral representations are then derived by introducing a method of extending short-time initial-value propagator formulas to global times via a simple stationary-phase asympotoic-equivalence approach. The integrals reduce to the nonuniform semiclassical propagators within the stationary-phase approximation, are uniform about caustics, and have integrand phases which properly account for the global-time phases in terms of appropriate Maslov indices. The initial-value integrals are also consistently derived via a canonical mapping procedure on the coordinate-space path integral. Initial-value integrals for time-dependent wave-function propgation are also given. Evaluation of the initial-value integral expressions do not require trajectory root searches for propagation.
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.
In this work, we derived a semiclassical approximation for the matrix elements of a quantum propagator in coherent states (CS) basis that avoids complex trajectories; it only involves real ones. For that purpose, we used the symplectically invariant semiclassical Weyl propagator obtained by performing a stationary phase approximation (SPA) for the path integral in the Weyl representation. After that, for the transformation to CS representation SPA is avoided; instead a quadratic expansion of the complex exponent is used. This procedure also allows us to express the semiclassical CS propagator uniquely in terms of the classical evolution of the initial point without the need of any root search typical of van Vleck--Gutzwiller-based propagators. For the case of chaotic Hamiltonian systems, the explicit time dependence of the CS propagator has been obtained. The comparison with a realistic chaotic system that derives from a quadratic Hamiltonian, the cat map, reveals that the expression here derived is exact up to quadratic Hamiltonian systems.