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

Classical Hamiltonian structures in wave packet dynamics

Stephen K. Gray, John M. Verosky

Published Apr 1, 1994• 78 citations
Publisher

Abstract

The general, N state matrix representation of the time-dependent Schrödinger equation is equivalent to an N degree of freedom classical Hamiltonian system. We describe how classical mechanical methods and ideas can be applied towards understanding and modeling exact quantum dynamics. Two applications are presented. First, we illustrate how qualitative insights may be gained by treating the two state problem with a time-dependent coupling. In the case of periodic coupling, Poincaré surfaces of section are used to view the quantum dynamics, and features such as the Floquet modes take on interesting interpretations. The second application illustrates computational implications by showing how Liouville’s theorem, or more generally the symplectic nature of classical Hamiltonian dynamics, provides a new perspective for carrying out numerical wave packet propagation. We show how certain simple and explicit symplectic integrators can be used to numerically propagate wave packets. The approach is illustrated with an application to the problem of a diatomic molecule interacting with a laser, although it and related approaches may be useful for describing a variety of problems.

DOI
10.1063/1.467219

Related papers

Classical numerical integrators for wave-packet dynamics

Jesús María Sanz-Serna, Ana M. Portillo

Gray and Verosky have recently studied the reformulation of the N-state matrix representation of the time-dependent Schrödinger equation as an N-degrees of freedom classical Hamiltonian system. This opens the possibility of using in quantum dynamics numerical integrators originally devised for classical mechanics. When the Hamiltonian matrix is time-dependent, Gray and Verosky suggest the use of a Magnus approximation before reducing the quantum system to its classical format. We show that Magnus approximations are not necessary and suggest an alternative technique. With the new technique it is possible to obtain simple integrators of arbitrarily high orders of accuracy that can be applied to all matrix Schrödinger problems with a, possibly time-dependent, real Hamiltonian matrix. The connection between the new approach and high-order split-operator methods is studied.

The Journal of Chemical Physics•Feb 8, 1996•61 citations

Hamilton‐jacobi dynamics for the solution of time dependent quantum problems I. Formalism and wave packet propagation in one dimension

Ersin Yurtsever, Jürgen Brickmann

Abstract Two methods for the numerical integration of the time‐dependent Schrödinger equation with given initial conditions (initial wave packet) are presented. The first method (method A) is based on the Schrödinger representation of the quantum‐dynamical system while the second one (method B) is based upon the intermediate representation. In both cases the quantum dynamical equation is transformed into a system of Hamilton‐Jacobi type equations of motion as occurring in multi particle classical dynamics, i.e. standard molecular dynamics techniques can be applied for the integration. The dynamics of a minimum uncertainty Gaussian wave packet in a strongly anharmonic oscillator is taken as an example.

Berichte der Bunsengesellschaft für physikalische Chemie•Apr 1, 1994•4 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

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

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

Trajectories and the perception of classical motion in the free propagation of wave packets

J S Briggs

Abstract The free propagation in time of a normalisable wave packet is the oldest problem of continuum quantum mechanics. Its motion from microscopic to macroscopic distance is the way in which most quantum systems are detected experimentally. Although much studied and analysed since 1927 and presented in many textbooks, here the problem is re‐appraised from the standpoint of semi‐classical mechanics. Particular aspects are the emergence of deterministic trajectories of particles emanating from a region of atomic dimension and the interpretation of the wave function as describing a single particle or an ensemble of identical particles. Of possible wave packets, that of Gaussian form is most studied due to the simple exact form of the time‐dependent solution in real and in momentum space. Furthermore, this form is important in laser optics. Here the equivalence of the time‐dependent Schrödinger equation to the paraxial equation for the propagation of light is demonstrated explicitly. This parallel helps to understand the relevance of trajectory concepts and the conditions necessary for the perception of quantum motion as classical. Key Points The paper addresses a problem of fundamental quantum mechanics, important in both physics and chemistry. Within physics, the subject is relevant to the disciplines of both quantum mechanics and classical optics. The matter presented is topical also for researchers in the history and philosophy of physics.

Natural Sciences•Apr 1, 2022•13 citations
The general, N state matrix representation of the time-dependent Schrödinger equation is equivalent to an N degree of freedom classical Hamiltonian system. We describe how classical mechanical methods and ideas can be applied towards understanding and modeling exact quantum dynamics. Two applications are presented. First, we illustrate how qualitative insights may be gained by treating the two state problem with a time-dependent coupling. In the case of periodic coupling, Poincaré surfaces of section are used to view the quantum dynamics, and features such as the Floquet modes take on interesting interpretations. The second application illustrates computational implications by showing how Liouville’s theorem, or more generally the symplectic nature of classical Hamiltonian dynamics, provides a new perspective for carrying out numerical wave packet propagation. We show how certain simple and explicit symplectic integrators can be used to numerically propagate wave packets. The approach is illustrated with an application to the problem of a diatomic molecule interacting with a laser, although it and related approaches may be useful for describing a variety of problems.
Gray and Verosky have recently studied the reformulation of the N-state matrix representation of the time-dependent Schrödinger equation as an N-degrees of freedom classical Hamiltonian system. This opens the possibility of using in quantum dynamics numerical integrators originally devised for classical mechanics. When the Hamiltonian matrix is time-dependent, Gray and Verosky suggest the use of a Magnus approximation before reducing the quantum system to its classical format. We show that Magnus approximations are not necessary and suggest an alternative technique. With the new technique it is possible to obtain simple integrators of arbitrarily high orders of accuracy that can be applied to all matrix Schrödinger problems with a, possibly time-dependent, real Hamiltonian matrix. The connection between the new approach and high-order split-operator methods is studied.
Abstract Two methods for the numerical integration of the time‐dependent Schrödinger equation with given initial conditions (initial wave packet) are presented. The first method (method A) is based on the Schrödinger representation of the quantum‐dynamical system while the second one (method B) is based upon the intermediate representation. In both cases the quantum dynamical equation is transformed into a system of Hamilton‐Jacobi type equations of motion as occurring in multi particle classical dynamics, i.e. standard molecular dynamics techniques can be applied for the integration. The dynamics of a minimum uncertainty Gaussian wave packet in a strongly anharmonic oscillator is taken as an example.
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.
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 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.
Abstract The free propagation in time of a normalisable wave packet is the oldest problem of continuum quantum mechanics. Its motion from microscopic to macroscopic distance is the way in which most quantum systems are detected experimentally. Although much studied and analysed since 1927 and presented in many textbooks, here the problem is re‐appraised from the standpoint of semi‐classical mechanics. Particular aspects are the emergence of deterministic trajectories of particles emanating from a region of atomic dimension and the interpretation of the wave function as describing a single particle or an ensemble of identical particles. Of possible wave packets, that of Gaussian form is most studied due to the simple exact form of the time‐dependent solution in real and in momentum space. Furthermore, this form is important in laser optics. Here the equivalence of the time‐dependent Schrödinger equation to the paraxial equation for the propagation of light is demonstrated explicitly. This parallel helps to understand the relevance of trajectory concepts and the conditions necessary for the perception of quantum motion as classical. Key Points The paper addresses a problem of fundamental quantum mechanics, important in both physics and chemistry. Within physics, the subject is relevant to the disciplines of both quantum mechanics and classical optics. The matter presented is topical also for researchers in the history and philosophy of physics.