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

Physical review. E

Semiclassical treatment of quantum propagation with nonlinear classical dynamics: A third-order thawed Gaussian approximation

Lucas Kocia, Anna Klales

Published Sep 14, 2016• 3 citations
Publisher

Abstract

The time-dependent WKB approximation for a coherent state is expanded to third order around a guiding real trajectory, allowing for the novel treatment of nonlinearity in its semiclassical dynamics, which is generally present in all physical systems far from the classical regime. The result is a closed-form solution consisting of a linear combination of Airy functions and their derivatives multiplied by an exponential. The expression's ability to capture nonlinearity is demonstrated by examining the autocorrelation of initial coherent states in anharmonic systems with few to many contributing periodic orbits. Its accuracy is compared to the quadratic expansion and found to be superior in regimes of ℏ where the curvature begins to be significant, as expected. Moreover, the expression is shown to be a real-trajectory uniformization over two coalescing saddle points that are emblematic of significant curvature. This extends real-trajectory time-dependent wave-packet semiclassical methods to highly anharmonic systems for the first time and establishes their regime of validity.

DOI
10.1103/physreve.94.032211

Related papers

A regular semiclassical approximation for the propagation of wave packets with complex trajectories

Fernando Parisio, Marcus A. M. de Aguiar

The semiclassical propagation of Gaussian wave packets by complex classical trajectories involves multiple contributing and noncontributing solutions interspersed by phase space caustics. Although the phase space caustics do not generally lie exactly on the relevant trajectories, they might strongly affect the semiclassical evolution depending on their proximity to them. In this paper, we derive a third-order regular semiclassical approximation which correctly accounts for the caustics and which is finite everywhere. We test the regular formula for the potential V ( x ) = 1/ x 2 , where the complex classical trajectories and phase space caustics can be computed analytically. We make a detailed analysis of the structure of the complex functions involved in the saddle point approximations and show how the changes in the steepest descent integration contour control both the contributing and noncontributing trajectories and the type of Airy function that appears in the regular approximation.

Journal of Physics A Mathematical and General•Oct 5, 2005•14 citations

Semiclassical propagation of coherent states using complex and real trajectories

Marcel Novaes, Marcus A. M. de Aguiar

A semiclassical approximation to the time evolution of coherent states may be derived from a saddle point approximation to the exact quantum propagator, and in general it involves complex classical dynamics. We generalize previous one-dimensional results to $d$ dimensions, and for the case $d=2$ we present several applications. We also consider other simple approximations that depend only on real classical trajectories, but are not initial value representations. These approximations are able to reproduce interference and tunneling effects and involve propagating a few classical initial conditions compatible with the quantum uncertainties.

arXiv:quant-ph/0504037•Physical Review A•Sep 12, 2005•19 citations

Generalized Gaussian wave packet dynamics: Integrable and chaotic systems

Harinder Pal, Manan Vyas, Steven Tomsovic

The ultimate semiclassical wave packet propagation technique is a complex, time-dependent Wentzel-Kramers-Brillouin method known as generalized Gaussian wave packet dynamics (GGWPD). It requires overcoming many technical difficulties in order to be carried out fully in practice. In its place roughly twenty years ago, linearized wave packet dynamics was generalized to methods that include sets of off-center, real trajectories for both classically integrable and chaotic dynamical systems that completely capture the dynamical transport. The connections between those methods and GGWPD are developed in a way that enables a far more practical implementation of GGWPD. The generally complex saddle-point trajectories at its foundation are found using a multidimensional Newton-Raphson root search method that begins with the set of off-center, real trajectories. This is possible because there is a one-to-one correspondence. The neighboring trajectories associated with each off-center, real trajectory form a path that crosses a unique saddle; there are exceptions that are straightforward to identify. The method is applied to the kicked rotor to demonstrate the accuracy improvement as a function of ℏ that comes with using the saddle-point trajectories.

arXiv:1510.08051•Physical review. E•Jan 21, 2016•22 citations

Semiclassical propagation of coherent states and wave packets: Hidden saddles

Huichao Wang, Steven Tomsovic

Semiclassical methods are extremely important in the subjects of wave-packet and coherent-state dynamics. Unfortunately, these essentially saddle-point approximations are considered nearly impossible to carry out in detail for systems with multiple degrees of freedom due to the difficulties of solving the resulting two-point boundary-value problems. However, recent developments have extended the applicability to a broader range of systems and circumstances. The most important advances are first to generate a set of real reference trajectories using appropriately reduced dimensional spaces of initial conditions, and second to feed that set into a Newton-Raphson search scheme to locate the exposed complex saddle trajectories. The arguments for this approach were based mostly on intuition and numerical verification. In this paper, the methods are put on a firmer theoretical foundation and then extended to incorporate saddles hidden from Newton-Raphson searches initiated with real trajectories. This hidden class of saddles is relevant to tunneling-type processes, but a hidden saddle can sometimes contribute just as much as or more than an exposed one. The distinctions between hidden and exposed saddles clarifies the interpretation of what constitutes tunneling for wave packets and coherent states in the time domain.

arXiv:2107.08799•Physical review. E•May 17, 2022•1 citation

Semiclassical Analysis of Wave Dynamics via Wigner Transforms: Applications to the Quantum Harmonic, Anharmonic, and Dirac Oscillator

August Burg

In this Master's thesis, we give a pedagogical introduction to the semiclassical analysis of wave dynamics in single-band and multi-band systems. We show how the semiclassical evolutions of Wigner functions can approximate the dynamics of wave packet centers and wave packet deformations. Three systems of increasing complexity are considered: the quantum harmonic oscillator and the quantum anharmonic oscillator as single-band systems, and the Dirac oscillator as a two-band system. We explain the origins of the first-order semiclassical corrections (Moyal, Littlejohn-Flynn, and Berry) in the equations of motion and examine their effects on predictive accuracy. Another focus is set on the conditions for validity of the semiclassical analysis. For both the quantum anharmonic and the Dirac oscillators, a dimensionless parameter controlling the accuracy of the approximation is identified. In multi-band systems, the formalism applies only to states within a single adiabatic subspace, which, for the Dirac oscillator, corresponds to an energy branch. This motivates the construction of single-branch almost coherent states: approximately Gaussian wave packets formed by superpositions of states from a single energy branch. Their centers undergo oscillations with an aperiodic envelope. These states enable comparison between semiclassical Wigner function evolution and the commonly used Sundaram-Niu formalism. We show that semiclassical Wigner function evolution captures deformation-related dynamics, including the initial decay of oscillations, which the Sundaram-Niu formalism cannot describe by construction. Additionally, exact quantum evolution reveals partial relocalization of these single-branch almost coherent states, a phenomenon that cannot be captured by the semiclassical evolution.

cond-mat.mes-hallCross-listedarXiv:2610.22414•Oct 16, 2026•0 citations

Semiclassical approximations based on complex trajectories

A. D. Ribeiro, Marcus A. M. de Aguiar, Michel Baranger

The semiclassical limit of the coherent state propagator $⟨{\mathbf{z}}^{\ensuremath{''}}\ensuremath{\mid}{e}^{\ensuremath{-}i\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}T∕\ensuremath{\hbar}}\ensuremath{\mid}{\mathbf{z}}^{\ensuremath{'}}⟩$ involves complex classical trajectories of the Hamiltonian $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{H}(\mathbf{u},\mathbf{v})=⟨\mathbf{v}\ensuremath{\mid}\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}\ensuremath{\mid}\mathbf{u}⟩$ satisfying $\mathbf{u}(0)={\mathbf{z}}^{\ensuremath{'}}$ and $\mathbf{v}(T)={\mathbf{z}}^{\ensuremath{''}}*$. In this work we study mostly the case ${\mathbf{z}}^{\ensuremath{'}}={\mathbf{z}}^{\ensuremath{''}}$. The propagator is then the return probability amplitude of a wave packet. We show that a plot of the exact return probability brings out the quantal images of the classical periodic orbits. Then we compare the exact return probability with its semiclassical approximation for a soft chaotic system with two degrees of freedom. We find two situations where classical trajectories satisfying the correct boundary conditions must be excluded from the semiclassical formula. The first occurs when the contribution of the trajectory to the propagator becomes exponentially large as $\ensuremath{\hbar}$ goes to zero. The second occurs when the contributing trajectories undergo bifurcations. Close to the bifurcation the semiclassical formula diverges. More interestingly, in the example studied, after the bifurcation, where more than one trajectory satisfying the boundary conditions exist, only one of them in fact contributes to the semiclassical formula, a phenomenon closely related to Stokes lines. When the contributions of these trajectories are filtered out, the semiclassical results show excellent agreement with the exact calculations.

Physical Review E•Jun 4, 2004•48 citations
The time-dependent WKB approximation for a coherent state is expanded to third order around a guiding real trajectory, allowing for the novel treatment of nonlinearity in its semiclassical dynamics, which is generally present in all physical systems far from the classical regime. The result is a closed-form solution consisting of a linear combination of Airy functions and their derivatives multiplied by an exponential. The expression's ability to capture nonlinearity is demonstrated by examining the autocorrelation of initial coherent states in anharmonic systems with few to many contributing periodic orbits. Its accuracy is compared to the quadratic expansion and found to be superior in regimes of ℏ where the curvature begins to be significant, as expected. Moreover, the expression is shown to be a real-trajectory uniformization over two coalescing saddle points that are emblematic of significant curvature. This extends real-trajectory time-dependent wave-packet semiclassical methods to highly anharmonic systems for the first time and establishes their regime of validity.
The semiclassical propagation of Gaussian wave packets by complex classical trajectories involves multiple contributing and noncontributing solutions interspersed by phase space caustics. Although the phase space caustics do not generally lie exactly on the relevant trajectories, they might strongly affect the semiclassical evolution depending on their proximity to them. In this paper, we derive a third-order regular semiclassical approximation which correctly accounts for the caustics and which is finite everywhere. We test the regular formula for the potential V ( x ) = 1/ x 2 , where the complex classical trajectories and phase space caustics can be computed analytically. We make a detailed analysis of the structure of the complex functions involved in the saddle point approximations and show how the changes in the steepest descent integration contour control both the contributing and noncontributing trajectories and the type of Airy function that appears in the regular approximation.
A semiclassical approximation to the time evolution of coherent states may be derived from a saddle point approximation to the exact quantum propagator, and in general it involves complex classical dynamics. We generalize previous one-dimensional results to ddd dimensions, and for the case d=2d=2d=2 we present several applications. We also consider other simple approximations that depend only on real classical trajectories, but are not initial value representations. These approximations are able to reproduce interference and tunneling effects and involve propagating a few classical initial conditions compatible with the quantum uncertainties.
The ultimate semiclassical wave packet propagation technique is a complex, time-dependent Wentzel-Kramers-Brillouin method known as generalized Gaussian wave packet dynamics (GGWPD). It requires overcoming many technical difficulties in order to be carried out fully in practice. In its place roughly twenty years ago, linearized wave packet dynamics was generalized to methods that include sets of off-center, real trajectories for both classically integrable and chaotic dynamical systems that completely capture the dynamical transport. The connections between those methods and GGWPD are developed in a way that enables a far more practical implementation of GGWPD. The generally complex saddle-point trajectories at its foundation are found using a multidimensional Newton-Raphson root search method that begins with the set of off-center, real trajectories. This is possible because there is a one-to-one correspondence. The neighboring trajectories associated with each off-center, real trajectory form a path that crosses a unique saddle; there are exceptions that are straightforward to identify. The method is applied to the kicked rotor to demonstrate the accuracy improvement as a function of ℏ that comes with using the saddle-point trajectories.
Semiclassical methods are extremely important in the subjects of wave-packet and coherent-state dynamics. Unfortunately, these essentially saddle-point approximations are considered nearly impossible to carry out in detail for systems with multiple degrees of freedom due to the difficulties of solving the resulting two-point boundary-value problems. However, recent developments have extended the applicability to a broader range of systems and circumstances. The most important advances are first to generate a set of real reference trajectories using appropriately reduced dimensional spaces of initial conditions, and second to feed that set into a Newton-Raphson search scheme to locate the exposed complex saddle trajectories. The arguments for this approach were based mostly on intuition and numerical verification. In this paper, the methods are put on a firmer theoretical foundation and then extended to incorporate saddles hidden from Newton-Raphson searches initiated with real trajectories. This hidden class of saddles is relevant to tunneling-type processes, but a hidden saddle can sometimes contribute just as much as or more than an exposed one. The distinctions between hidden and exposed saddles clarifies the interpretation of what constitutes tunneling for wave packets and coherent states in the time domain.
In this Master's thesis, we give a pedagogical introduction to the semiclassical analysis of wave dynamics in single-band and multi-band systems. We show how the semiclassical evolutions of Wigner functions can approximate the dynamics of wave packet centers and wave packet deformations. Three systems of increasing complexity are considered: the quantum harmonic oscillator and the quantum anharmonic oscillator as single-band systems, and the Dirac oscillator as a two-band system. We explain the origins of the first-order semiclassical corrections (Moyal, Littlejohn-Flynn, and Berry) in the equations of motion and examine their effects on predictive accuracy. Another focus is set on the conditions for validity of the semiclassical analysis. For both the quantum anharmonic and the Dirac oscillators, a dimensionless parameter controlling the accuracy of the approximation is identified. In multi-band systems, the formalism applies only to states within a single adiabatic subspace, which, for the Dirac oscillator, corresponds to an energy branch. This motivates the construction of single-branch almost coherent states: approximately Gaussian wave packets formed by superpositions of states from a single energy branch. Their centers undergo oscillations with an aperiodic envelope. These states enable comparison between semiclassical Wigner function evolution and the commonly used Sundaram-Niu formalism. We show that semiclassical Wigner function evolution captures deformation-related dynamics, including the initial decay of oscillations, which the Sundaram-Niu formalism cannot describe by construction. Additionally, exact quantum evolution reveals partial relocalization of these single-branch almost coherent states, a phenomenon that cannot be captured by the semiclassical evolution.
The semiclassical limit of the coherent state propagator $⟨{\mathbf{z}}^{\ensuremath{''}}\ensuremath{\mid}{e}^{\ensuremath{-}i\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}T∕\ensuremath{\hbar}}\ensuremath{\mid}{\mathbf{z}}^{\ensuremath{'}}⟩$ involves complex classical trajectories of the Hamiltonian $\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{H}(\mathbf{u},\mathbf{v})=⟨\mathbf{v}\ensuremath{\mid}\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{H}\ensuremath{\mid}\mathbf{u}⟩$ satisfying $\mathbf{u}(0)={\mathbf{z}}^{\ensuremath{'}}$ and $\mathbf{v}(T)={\mathbf{z}}^{\ensuremath{''}}*$. In this work we study mostly the case ${\mathbf{z}}^{\ensuremath{'}}={\mathbf{z}}^{\ensuremath{''}}$. The propagator is then the return probability amplitude of a wave packet. We show that a plot of the exact return probability brings out the quantal images of the classical periodic orbits. Then we compare the exact return probability with its semiclassical approximation for a soft chaotic system with two degrees of freedom. We find two situations where classical trajectories satisfying the correct boundary conditions must be excluded from the semiclassical formula. The first occurs when the contribution of the trajectory to the propagator becomes exponentially large as $\ensuremath{\hbar}$ goes to zero. The second occurs when the contributing trajectories undergo bifurcations. Close to the bifurcation the semiclassical formula diverges. More interestingly, in the example studied, after the bifurcation, where more than one trajectory satisfying the boundary conditions exist, only one of them in fact contributes to the semiclassical formula, a phenomenon closely related to Stokes lines. When the contributions of these trajectories are filtered out, the semiclassical results show excellent agreement with the exact calculations.