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

Berichte der Bunsengesellschaft für physikalische Chemie

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

Published Apr 1, 1994• 4 citations
Publisher

Abstract

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.

DOI
10.1002/bbpc.19940980404

Related papers

Classical Hamiltonian structures in wave packet dynamics

Stephen K. Gray, John M. Verosky

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.

The Journal of Chemical Physics•Apr 1, 1994•78 citations

Numerical solutions of the quantum Hamilton-Jacobi equation and WKB like representations for one dimensional wave functions

M. Fusco Girard

By means of numerical solutions of the quantum Hamilton Jacobi equation, a general WKB-like representation for one-dimensional wave functions is obtained. This representation is unique in the classically forbidden regions, while in the allowed one, each wave function corresponds to a one parameter family of solutions of the QHJE. The method has been applied to various systems, with different energies and initial conditions. In all investigated cases, the wave functions so obtained accurately reproduce the solutions of the Schrödinger equation, analytically or numerically computed by other ways. Some results for harmonic oscillator and radial Coulomb motion are presented.

arXiv:1403.0825•Mar 4, 2014•0 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

Particle Method for the Numerical Solution of the Time-Dependent Schrödinger Equation

Jerome H. Weiner, Attila Aşkar

A method for the solution of the time-dependent Schrödinger equation in two dimensions is presented which is based on the hydrodynamic analogy to quantum mechanics. The continuum introduced in that analogy is approximated by a finite number of particles, whose trajectories are computed. The procedure is applied here to the dynamics of a Gaussian wave packet on a two-dimensional quadratic potential surface containing a saddle point. Comparison with the known analytical solution indicates that the method is capable of accurate results with substantially less computer time required than in methods previously employed. This particular check problem was chosen because of its close relationship to the wave-packet dynamics of chemical reaction rates.

The Journal of Chemical Physics•Apr 15, 1971•71 citations

Theoretical methods for the study of quantum dynamics

Jonathan E. Chen

Several theoretical methods for the computation of quantum dynamical quantities are formulated, implemented, and applied, in the overall context of overcoming the exponential scaling property of quantum mechanics. The convergence of the approximate forward-backward semiclassical dynamics (FBSD) method has been accelerated by working with only the imaginary part of the real-time position and momentum time correlation functions. Simple manipulations allow the computationally favorable limit of a single factorization (“bead”) of the Boltzmann operator to yield greatly improved results. Numerical examples for a one-dimensional anharmonic oscillator and Lennard-Jones fluids are presented. Bohm's trajectory formulation of quantum mechanics is an exact mapping of the time-dependent Schrödinger equation onto classical fluid dynamics, plus one extra term responsible for all non-classical effects, the dreaded quantum potential. A computational procedure has been introduced to stably solve Bohm's hydrodynamic equations of motion in the forward-backward quantum dynamics (FBQD) formulation of time correlation functions, and applied to one-dimensional bound systems at zero and low temperature. A key step in the formulation of this procedure was the use of Hamilton's law of varying action, a generalization of Hamilton's principle of stationary action applicable to initial-value problems. This has resulted in the first reported synthetic computation of non-stationary quantum trajectories for a fully anharmonic oscillator. Success in one dimension has laid the groundwork for future extension to a rigorous quantum-classical scheme, which is briefly discussed. Other unpublished methods for obtaining time correlation functions and propagators in quantum dynamics are also discussed.

•May 19, 2010•0 citations

Solution to the initial value problem for the quantum nonlinear Schrödinger equation

Bernard Yurke, Mary J. Potasek

The quantum nonlinear Schrödinger equation, which models the propagation of light through dispersive nonlinear waveguides, provides an integrable quantum field theory that has been solved by a number of methods. Most recently, Gutkin [ J. Func. Anal.77, 327 ( 1988)] has developed an interwining operator technique for obtaining the time evolution of the field operators. Using Gutkin’s formalism, we show how to obtain the exact time evolution of an initial Schrödinger state vector. As an example, the two-particle case is carried out explicitly. Numerical results for the time evolution of the position probability distribution are presented for the case when the two particles initially occupy the same Gaussian wave packet.

Journal of the Optical Society of America B•Jun 1, 1989•22 citations
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.
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.
By means of numerical solutions of the quantum Hamilton Jacobi equation, a general WKB-like representation for one-dimensional wave functions is obtained. This representation is unique in the classically forbidden regions, while in the allowed one, each wave function corresponds to a one parameter family of solutions of the QHJE. The method has been applied to various systems, with different energies and initial conditions. In all investigated cases, the wave functions so obtained accurately reproduce the solutions of the Schrödinger equation, analytically or numerically computed by other ways. Some results for harmonic oscillator and radial Coulomb motion are presented.
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.
A method for the solution of the time-dependent Schrödinger equation in two dimensions is presented which is based on the hydrodynamic analogy to quantum mechanics. The continuum introduced in that analogy is approximated by a finite number of particles, whose trajectories are computed. The procedure is applied here to the dynamics of a Gaussian wave packet on a two-dimensional quadratic potential surface containing a saddle point. Comparison with the known analytical solution indicates that the method is capable of accurate results with substantially less computer time required than in methods previously employed. This particular check problem was chosen because of its close relationship to the wave-packet dynamics of chemical reaction rates.
Several theoretical methods for the computation of quantum dynamical quantities are formulated, implemented, and applied, in the overall context of overcoming the exponential scaling property of quantum mechanics. The convergence of the approximate forward-backward semiclassical dynamics (FBSD) method has been accelerated by working with only the imaginary part of the real-time position and momentum time correlation functions. Simple manipulations allow the computationally favorable limit of a single factorization (“bead”) of the Boltzmann operator to yield greatly improved results. Numerical examples for a one-dimensional anharmonic oscillator and Lennard-Jones fluids are presented. Bohm's trajectory formulation of quantum mechanics is an exact mapping of the time-dependent Schrödinger equation onto classical fluid dynamics, plus one extra term responsible for all non-classical effects, the dreaded quantum potential. A computational procedure has been introduced to stably solve Bohm's hydrodynamic equations of motion in the forward-backward quantum dynamics (FBQD) formulation of time correlation functions, and applied to one-dimensional bound systems at zero and low temperature. A key step in the formulation of this procedure was the use of Hamilton's law of varying action, a generalization of Hamilton's principle of stationary action applicable to initial-value problems. This has resulted in the first reported synthetic computation of non-stationary quantum trajectories for a fully anharmonic oscillator. Success in one dimension has laid the groundwork for future extension to a rigorous quantum-classical scheme, which is briefly discussed. Other unpublished methods for obtaining time correlation functions and propagators in quantum dynamics are also discussed.
The quantum nonlinear Schrödinger equation, which models the propagation of light through dispersive nonlinear waveguides, provides an integrable quantum field theory that has been solved by a number of methods. Most recently, Gutkin [ J. Func. Anal.77, 327 ( 1988)] has developed an interwining operator technique for obtaining the time evolution of the field operators. Using Gutkin’s formalism, we show how to obtain the exact time evolution of an initial Schrödinger state vector. As an example, the two-particle case is carried out explicitly. Numerical results for the time evolution of the position probability distribution are presented for the case when the two particles initially occupy the same Gaussian wave packet.