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

Time-dependent approach to semiclassical dynamics

Eric J. Heller

Published Feb 15, 1975• 1,619 citations
Publisher

Abstract

In this paper we develop a new approach to semiclassical dynamics which exploits the fact that extended wavefunctions for heavy particles (or particles in harmonic potentials) may be decomposed into time−dependent wave packets, which spread minimally and which execute classical or nearly classical trajectories. A Gaussian form for the wave packets is assumed and equations of motion are derived for the parameters characterizing the Gaussians. If the potential (which may be nonseparable in many coordinates) is expanded in a Taylor series about the instantaneous center of the (many−particle) wave packet, and up to quadratic terms are kept, we find the classical parameters of the wave packet (positions, momenta) obey Hamilton’s equation of motion. Quantum parameters (wave packet spread, phase factor, correlation terms, etc.) obey similar first order quantum equations. The center of the wave packet is shown to acquire a phase equal to the action integral along the classical path. State−specific quantum information is obtained from the wave packet trajectories by use of the superposition principle and projection techniques. Successful numerical application is made to the collinear He + H2 system widely used as a test case. Classically forbidden transitions are accounted for and obtained in the same manner as the classically allowed transitions; turning points present no difficulties and flux is very nearly conserved.

DOI
10.1063/1.430620

Related papers

A strategy for time dependent quantum mechanical calculations using a Gaussian wave packet representation of the wave function

Shin‐ichi Sawada, Robert Heather, Bret Jackson +1 more

We develop methodology for performing time dependent quantum mechanical calculations by representing the wave function as a sum of Gaussian wave packets (GWP) each characterized by a set of parameters such as width, position, momentum, and phase. The problem of computing the time evolution of the wave function is thus reduced to that of finding the time evolution of the parameters in the Gaussians. This parameter motion is determined by minimizing the error made by replacing the exact wave function in the time dependent Schrödinger equation with its Gaussian representation approximant. This leads to first order differential equations for the time dependence of the parameters, and those describing the packet position and the momentum of each packet have some resemblance with the classical equations of motion. The paper studies numerically the strategy needed to achieve the best GWP representation of time dependent processes. The issues discussed are the representation of the initial wave function, the numerical stability and the solution of the differential equations giving the evolution of the parameters, and the analysis of the final wave function. Extensive comparisons are made with an approximate method which assumes that the Gaussians are independent and their width is smaller than the length scale over which the potential changes. This approximation greatly simplifies the calculations and has the advantage of a greater resemblance to classical mechanics, thus being more intuitive. We find, however, that its range of applications is limited to problems involving localized degrees of freedom that participate in the dynamic process for a very short time. Finally we give particular attention to the notion that the GWP representation of the wave function reduces the dynamics of one quantum degree of freedom to that of a set of pseudoparticles (each represented by one packet) moving according to a pseudoclassical (i.e., classical-like) mechanics whose ‘‘phase space’’ is described by a position and momentum as well as a complex phase and width.

The Journal of Chemical Physics•Sep 15, 1985•205 citations

Symplectic semiclassical wave packet dynamics

Tomoki Ohsawa, Melvin Boon-tiong Leok

The paper gives a symplectic-geometric account of semiclassical Gaussian wave packet dynamics. We employ geometric techniques to ‘strip away’ the symplectic structure behind the time-dependent Schrödinger equation and incorporate it into semiclassical wave packet dynamics. We show that the Gaussian wave packet dynamics is a Hamiltonian system with respect to the symplectic structure, apply the theory of symplectic reduction and reconstruction to the dynamics, and discuss dynamic and geometric phases in semiclassical mechanics. A simple harmonic oscillator example is worked out to illustrate the results. We show that the reduced semiclassical harmonic oscillator dynamics is completely integrable by finding the action–angle coordinates for the system, and calculate the associated dynamic and geometric phases explicitly. We also propose an asymptotic approximation of the potential term that provides a practical semiclassical correction term to the approximation by Heller. Numerical results for a simple one-dimensional example show that the semiclassical correction term realizes a semiclassical tunneling.

arXiv:1302.1139•Journal of Physics A Mathematical and Theoretical•Sep 19, 2013•31 citations

Complex trajectory method in semiclassical propagation of wave packets

Marc Boiron, M. Lombardi

We propose a semiclassical wave packet propagation method relying on classical trajectories in a complex phase space. It is based on the Schrödinger wave equation and the usual expansion with respect to ℏ, except that the amplitude of the wave packet is taken into account at the very zeroth order, unlike in the usual WKB method where it is treated as a corrective or first order term. Formally, it amounts to making both the wavelength and the width of the wave packet tend to zero with ℏ. The action and consequently the classical trajectories derived are complex. This method is tested successfully in many cases, analytically or numerically, including the bounce and even the splitting of the wave packet. Our method appears to be much more accurate than the WKB method while less computationally demanding than the Van-Vleck formula. Moreover, it has a particularly interesting property: the singularities (caustics) of the usual semiclassical theories do not appear in this formalism in all cases tested.

The Journal of Chemical Physics•Mar 1, 1998•32 citations

Hamiltonian formulation of quantum mechanics with semiclassical implications. II. Variational treatment

Kenneth G. Kay

It is shown that application of the time-dependent variational principle to trial wave functions of a particular form results in a description of the quantum system as a classical system consisting of two mutually interacting parts: the ordinary classical analog of the quantum system and a (generally finite) number of additional degrees of freedom. This holds out the prospect of simulating quantum dynamics by coupling the corresponding classical system to additional degrees of freedom, and it is pointed out that such calculations have, in effect, already been performed. Since the trial wave functions that qualify for this treatment may be arbitrarily accurate, there is no limitation to the accuracy of the resulting classical description. In a particular case, the canonical coordinates and momenta for all degrees of freedom of the overall classical system are identified and investigated. In general, certain parameters appearing in the trial wave functions treated in this work may become redundant at certain times, leading to singularities in the Hamiltonian equations of motion. It is shown, however, that in certain such cases, a classical formalism for the quantum dynamics remains valid if the Hamiltonian function is suitably modified. Calculations applying the present Hamiltonian formalism are carried out that demonstrate that it has better convergence properties than an alternative, nonvariational, Hamiltonian technique suggested by previous work. These numerical studies also suggest that the observed spreading of wave packets and development of complicated structure in wave functions may be interpreted classically in terms of energy transfer among the various degrees of freedom in the overall classical system.

Physical Review A•Aug 1, 1992•24 citations

Generalized Gaussian wave packet dynamics

Daniel Huber, Eric J. Heller

We propose an extension of the semiclassical Gaussian wave packet dynamics to eliminate the three main restrictions of this method. The first restriction is that the wave packet is forced to remain Gaussian. This is correct only for quadratic, linear, or constant potentials. The second restriction is that the method is, in general, not able to treat most classically forbidden processes. The third restriction is that the norm is conserved only for Gaussian wave packets. For a superposition of Gaussians this is no longer true. We can eliminate these restrictions by an extension of the method into complex phase space, keeping time real.

The Journal of Chemical Physics•Nov 1, 1987•179 citations

Classical mechanics and the spreading of localized wave packets in condensed phase molecular systems

Michael F. Herman, D. F. Coker

The relationship between the diverging of classical trajectories in chaotic many-body systems, the spreading of quantum wave packets, and the validity and use of classical molecular dynamics is explored. This analysis, which is based on the semiclassical description of wave function propagation in terms of a weighted integration over a traveling fixed width coherent state basis, suggests that the exponential divergence of nearby classical trajectories in chaotic many-body systems should result in the rapid delocalization of an initially localized quantum wave packet describing the state of the system. Thus the justification for the use of classical molecular dynamics procedures for these supposedly classical systems cannot be based on the picture of the system wave function remaining localized as its center follows a nearly classical trajectory. The quantum evolution of the system density, on the other hand, requires two propagators, and each of these propagators is represented as an integration over trajectories in the semiclassical picture. The interference between the contributions from these two integrations over classical trajectories focuses the analysis on the most important points in this trajectory pair space, which are shown to occur when both trajectories in the pair are the same. Given reasonable assumptions for the initial density for a system that is expected to be well described by classical molecular dynamics, and given an appropriate choice for the width of the coherent state basis which is employed in the semiclassical description, it is shown that the semiclassical expressions for time dependent observables and correlation functions reduce the purely classical expressions, despite the fact that an initially localized wave packet would rapidly delocalize for the same system.

The Journal of Chemical Physics•Aug 1, 1999•43 citations
In this paper we develop a new approach to semiclassical dynamics which exploits the fact that extended wavefunctions for heavy particles (or particles in harmonic potentials) may be decomposed into time−dependent wave packets, which spread minimally and which execute classical or nearly classical trajectories. A Gaussian form for the wave packets is assumed and equations of motion are derived for the parameters characterizing the Gaussians. If the potential (which may be nonseparable in many coordinates) is expanded in a Taylor series about the instantaneous center of the (many−particle) wave packet, and up to quadratic terms are kept, we find the classical parameters of the wave packet (positions, momenta) obey Hamilton’s equation of motion. Quantum parameters (wave packet spread, phase factor, correlation terms, etc.) obey similar first order quantum equations. The center of the wave packet is shown to acquire a phase equal to the action integral along the classical path. State−specific quantum information is obtained from the wave packet trajectories by use of the superposition principle and projection techniques. Successful numerical application is made to the collinear He + H2 system widely used as a test case. Classically forbidden transitions are accounted for and obtained in the same manner as the classically allowed transitions; turning points present no difficulties and flux is very nearly conserved.
We develop methodology for performing time dependent quantum mechanical calculations by representing the wave function as a sum of Gaussian wave packets (GWP) each characterized by a set of parameters such as width, position, momentum, and phase. The problem of computing the time evolution of the wave function is thus reduced to that of finding the time evolution of the parameters in the Gaussians. This parameter motion is determined by minimizing the error made by replacing the exact wave function in the time dependent Schrödinger equation with its Gaussian representation approximant. This leads to first order differential equations for the time dependence of the parameters, and those describing the packet position and the momentum of each packet have some resemblance with the classical equations of motion. The paper studies numerically the strategy needed to achieve the best GWP representation of time dependent processes. The issues discussed are the representation of the initial wave function, the numerical stability and the solution of the differential equations giving the evolution of the parameters, and the analysis of the final wave function. Extensive comparisons are made with an approximate method which assumes that the Gaussians are independent and their width is smaller than the length scale over which the potential changes. This approximation greatly simplifies the calculations and has the advantage of a greater resemblance to classical mechanics, thus being more intuitive. We find, however, that its range of applications is limited to problems involving localized degrees of freedom that participate in the dynamic process for a very short time. Finally we give particular attention to the notion that the GWP representation of the wave function reduces the dynamics of one quantum degree of freedom to that of a set of pseudoparticles (each represented by one packet) moving according to a pseudoclassical (i.e., classical-like) mechanics whose ‘‘phase space’’ is described by a position and momentum as well as a complex phase and width.
The paper gives a symplectic-geometric account of semiclassical Gaussian wave packet dynamics. We employ geometric techniques to ‘strip away’ the symplectic structure behind the time-dependent Schrödinger equation and incorporate it into semiclassical wave packet dynamics. We show that the Gaussian wave packet dynamics is a Hamiltonian system with respect to the symplectic structure, apply the theory of symplectic reduction and reconstruction to the dynamics, and discuss dynamic and geometric phases in semiclassical mechanics. A simple harmonic oscillator example is worked out to illustrate the results. We show that the reduced semiclassical harmonic oscillator dynamics is completely integrable by finding the action–angle coordinates for the system, and calculate the associated dynamic and geometric phases explicitly. We also propose an asymptotic approximation of the potential term that provides a practical semiclassical correction term to the approximation by Heller. Numerical results for a simple one-dimensional example show that the semiclassical correction term realizes a semiclassical tunneling.
We propose a semiclassical wave packet propagation method relying on classical trajectories in a complex phase space. It is based on the Schrödinger wave equation and the usual expansion with respect to ℏ, except that the amplitude of the wave packet is taken into account at the very zeroth order, unlike in the usual WKB method where it is treated as a corrective or first order term. Formally, it amounts to making both the wavelength and the width of the wave packet tend to zero with ℏ. The action and consequently the classical trajectories derived are complex. This method is tested successfully in many cases, analytically or numerically, including the bounce and even the splitting of the wave packet. Our method appears to be much more accurate than the WKB method while less computationally demanding than the Van-Vleck formula. Moreover, it has a particularly interesting property: the singularities (caustics) of the usual semiclassical theories do not appear in this formalism in all cases tested.
It is shown that application of the time-dependent variational principle to trial wave functions of a particular form results in a description of the quantum system as a classical system consisting of two mutually interacting parts: the ordinary classical analog of the quantum system and a (generally finite) number of additional degrees of freedom. This holds out the prospect of simulating quantum dynamics by coupling the corresponding classical system to additional degrees of freedom, and it is pointed out that such calculations have, in effect, already been performed. Since the trial wave functions that qualify for this treatment may be arbitrarily accurate, there is no limitation to the accuracy of the resulting classical description. In a particular case, the canonical coordinates and momenta for all degrees of freedom of the overall classical system are identified and investigated. In general, certain parameters appearing in the trial wave functions treated in this work may become redundant at certain times, leading to singularities in the Hamiltonian equations of motion. It is shown, however, that in certain such cases, a classical formalism for the quantum dynamics remains valid if the Hamiltonian function is suitably modified. Calculations applying the present Hamiltonian formalism are carried out that demonstrate that it has better convergence properties than an alternative, nonvariational, Hamiltonian technique suggested by previous work. These numerical studies also suggest that the observed spreading of wave packets and development of complicated structure in wave functions may be interpreted classically in terms of energy transfer among the various degrees of freedom in the overall classical system.
We propose an extension of the semiclassical Gaussian wave packet dynamics to eliminate the three main restrictions of this method. The first restriction is that the wave packet is forced to remain Gaussian. This is correct only for quadratic, linear, or constant potentials. The second restriction is that the method is, in general, not able to treat most classically forbidden processes. The third restriction is that the norm is conserved only for Gaussian wave packets. For a superposition of Gaussians this is no longer true. We can eliminate these restrictions by an extension of the method into complex phase space, keeping time real.
The relationship between the diverging of classical trajectories in chaotic many-body systems, the spreading of quantum wave packets, and the validity and use of classical molecular dynamics is explored. This analysis, which is based on the semiclassical description of wave function propagation in terms of a weighted integration over a traveling fixed width coherent state basis, suggests that the exponential divergence of nearby classical trajectories in chaotic many-body systems should result in the rapid delocalization of an initially localized quantum wave packet describing the state of the system. Thus the justification for the use of classical molecular dynamics procedures for these supposedly classical systems cannot be based on the picture of the system wave function remaining localized as its center follows a nearly classical trajectory. The quantum evolution of the system density, on the other hand, requires two propagators, and each of these propagators is represented as an integration over trajectories in the semiclassical picture. The interference between the contributions from these two integrations over classical trajectories focuses the analysis on the most important points in this trajectory pair space, which are shown to occur when both trajectories in the pair are the same. Given reasonable assumptions for the initial density for a system that is expected to be well described by classical molecular dynamics, and given an appropriate choice for the width of the coherent state basis which is employed in the semiclassical description, it is shown that the semiclassical expressions for time dependent observables and correlation functions reduce the purely classical expressions, despite the fact that an initially localized wave packet would rapidly delocalize for the same system.