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

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, Horia I. Metiu

Published Sep 15, 1985• 205 citations
Publisher

Abstract

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.

DOI
10.1063/1.449204

Related papers

Evolution of the wave function's shape in a time-dependent harmonic potential

Etera R. Livine

Abstract An effective operational approach to quantum mechanics is to focus on the evolution of wave packets, for which the wave function can be seen in the semi-classical regime as representing a classical motion dressed with extra degrees of freedom describing the shape of the wave packet and its fluctuations. These quantum dressing are independent degrees of freedom, mathematically encoded in the higher moments of the wave function. We review how to extract the effective dynamics for Gaussian wave packets evolving according to the Schrödinger equation with time-dependent potential in a 1 + 1-dimensional spacetime, and derive the equations of motion for the quadratic uncertainty. We then show how to integrate the evolution of all the higher moments for a general wave function in a time-dependent harmonic potential.

arXiv:2305.03847•Europhysics Letters (EPL)•Aug 1, 2023•2 citations

Generalized Gaussian wave packet dynamics, Schrödinger equation, and stationary phase approximation

Daniel Huber, Eric J. Heller, Robert G. Littlejohn

We have recently published a new semiclassical method, generalized Gaussian wave packet dynamics, which extends Gaussian wave packet dynamics into complex phase space. Although we were able to give an accurate formulation of the method, we had at the time of writing that paper only an intuitive, heuristic understanding of the deeper causes which make the method work. A more mathematical understanding was needed. To close this gap we show in this paper the equivalence of the new method with a first order expansion of ℏ of the Schrödinger equation. We further prove that the new method is equivalent to the stationary phase approximation, using the usual WKB formula for the propagator. The latter equivalence enables us to show that all the symmetry properties of time-dependent quantum mechanics also hold in the new semiclassical theory. Finally, we provide some elaboration of the method, and clarify several issues that were not discussed before. With this new insight we are able to formulate a simple rule for the calculation of semiclassical wave functions that contain contributions from more than one branch. This corrects for the divergence of semiclassical wave functions near caustics, a problem that we encountered in the preceding paper.

The Journal of Chemical Physics•Aug 15, 1988•118 citations

The Time-Dependent Schrödinger Equation

Reinhold A. Bertlmann, Nicolai Friis

Abstract The second chapter introduces the notion of wave functions and their equation of motion proposed by Erwin Schrödinger. The interpretation and normalization of the wave function are discussed, leading to the continuity equation. The chapter also introduces basic mathematical operations that can be carried out for wave functions and operators, describing e.g., how to calculate scalar products and commutators, as well as expectation values and variances of observables. With these tools at hand, Heisenberg’s uncertainty principle is re-examined, focusing on the example of a Gaussian wave packet, which is studied both in position space and momentum space. We also discuss the time evolution of general states and of the exemplary wave packet. Finally, we will see how classical physics is recovered as a limiting case of the quantum-mechanical description

•Sep 7, 2023•0 citations

Time-dependent approach to semiclassical dynamics

Eric J. Heller

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.

The Journal of Chemical Physics•Feb 15, 1975•1,619 citations

Lewis-Riesenfeld approach to the solutions of the Schrödinger equation in the presence of a time-dependent linear potential

Pi‐Gang Luan, Chi‐Shung Tang

We reexamine the general solution of a Schr\"odinger equation in the presence of a time-dependent linear potential in configuration space based on the Lewis-Riesenfeld framework. For comparison, we also solve the problem in momentum space and then Fourier transform the solution to get the general wave function. Appropriately choosing the weight function in the latter method, we can obtain the same wave function as the former method. It is found that a non-Hermitian time-dependent linear invariant can be used to obtain Gaussian-type wave-packet solutions of the time-dependent system. This operator is a specific linear combination of the initial momentum and initial position operators. This fact indicates that the constants of integration such as the initial position and initial momentum that determine the classical motion play important roles in the time-dependent quantum system. The eigenfunction of the linear invariant is interpreted as a wave packet with a ``center of mass'' moving along the classical trajectory, while the ratio between the coefficients of the initial position and initial momentum determines the width of the wave packet.

arXiv:quant-ph/0309174•Physical Review A•Jan 5, 2005•34 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
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 An effective operational approach to quantum mechanics is to focus on the evolution of wave packets, for which the wave function can be seen in the semi-classical regime as representing a classical motion dressed with extra degrees of freedom describing the shape of the wave packet and its fluctuations. These quantum dressing are independent degrees of freedom, mathematically encoded in the higher moments of the wave function. We review how to extract the effective dynamics for Gaussian wave packets evolving according to the Schrödinger equation with time-dependent potential in a 1 + 1-dimensional spacetime, and derive the equations of motion for the quadratic uncertainty. We then show how to integrate the evolution of all the higher moments for a general wave function in a time-dependent harmonic potential.
We have recently published a new semiclassical method, generalized Gaussian wave packet dynamics, which extends Gaussian wave packet dynamics into complex phase space. Although we were able to give an accurate formulation of the method, we had at the time of writing that paper only an intuitive, heuristic understanding of the deeper causes which make the method work. A more mathematical understanding was needed. To close this gap we show in this paper the equivalence of the new method with a first order expansion of ℏ of the Schrödinger equation. We further prove that the new method is equivalent to the stationary phase approximation, using the usual WKB formula for the propagator. The latter equivalence enables us to show that all the symmetry properties of time-dependent quantum mechanics also hold in the new semiclassical theory. Finally, we provide some elaboration of the method, and clarify several issues that were not discussed before. With this new insight we are able to formulate a simple rule for the calculation of semiclassical wave functions that contain contributions from more than one branch. This corrects for the divergence of semiclassical wave functions near caustics, a problem that we encountered in the preceding paper.
Abstract The second chapter introduces the notion of wave functions and their equation of motion proposed by Erwin Schrödinger. The interpretation and normalization of the wave function are discussed, leading to the continuity equation. The chapter also introduces basic mathematical operations that can be carried out for wave functions and operators, describing e.g., how to calculate scalar products and commutators, as well as expectation values and variances of observables. With these tools at hand, Heisenberg’s uncertainty principle is re-examined, focusing on the example of a Gaussian wave packet, which is studied both in position space and momentum space. We also discuss the time evolution of general states and of the exemplary wave packet. Finally, we will see how classical physics is recovered as a limiting case of the quantum-mechanical description
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 reexamine the general solution of a Schr\"odinger equation in the presence of a time-dependent linear potential in configuration space based on the Lewis-Riesenfeld framework. For comparison, we also solve the problem in momentum space and then Fourier transform the solution to get the general wave function. Appropriately choosing the weight function in the latter method, we can obtain the same wave function as the former method. It is found that a non-Hermitian time-dependent linear invariant can be used to obtain Gaussian-type wave-packet solutions of the time-dependent system. This operator is a specific linear combination of the initial momentum and initial position operators. This fact indicates that the constants of integration such as the initial position and initial momentum that determine the classical motion play important roles in the time-dependent quantum system. The eigenfunction of the linear invariant is interpreted as a wave packet with a ``center of mass'' moving along the classical trajectory, while the ratio between the coefficients of the initial position and initial momentum determines the width of the wave packet.
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.