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

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

Jerome H. Weiner, Attila Aşkar

Published Apr 15, 1971• 71 citations
Publisher

Abstract

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.

DOI
10.1063/1.1675377

Related papers

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.

International Journal of Quantum Chemistry•Jun 18, 2009•0 citations

Wave Packet Dynamics on Two-Dimensional Quadratic Potential Surfaces

Attila Aşkar, Jerome H. Weiner

Computer-generated figures within the framework of the hydrodynamic analogy to the time-dependent Schrödinger equation are presented for the motion of Gaussian wave packets on two-dimensional potential surfaces. The three cases treated illustrate (a) motion in a potential well, (b) tunneling through a saddle point, and (c) scattering by a potential barrier.

American Journal of Physics•Oct 1, 1971•26 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

The Quantum Hydrodynamic Description of Tunneling

Brian K. Kendrick

The quantum hydrodynamic approach is based on the de Broglie-Bohm formulation of quantum mechanics. The resulting fluid-like equations of motion describe the flow of probability and an accurate solution to these equations is equivalent to solving the time-dependent Schroedinger equation. Furthermore, the hydrodynamic approach provides new insight into the mechanisms as well as an alternative computational approach for treating tunneling phenomena. New concepts include well-defined 'quantum trajectories', 'quantum potential', and 'quantum force' all of which have classical analogues. The quantum potential and its associated force give rise to all quantum mechanical effects such as zero point energy, tunneling, and interference. A new numerical approach called the Iterative Finite Difference Method (IFDM) will be discussed. The IFDM is used to solve the set of non-linear coupled hydrodynamic equations. It is 2nd-order accurate in both space and time and exhibits exponential convergence with respect to the iteration count. The stability and computational efficiency of the IFDM is significantly improved by using a 'smart' Eulerian grid which has the same computational advantages as a Lagrangian or Arbitrary Lagrangian Eulerian (ALE) grid. The IFDM is also capable of treating anharmonic potentials. Example calculations using the IFDM will be presented which include: a one-dimensional Gaussian wave packet tunneling through an Eckart barrier, a one-dimensional bound-state Morse oscillator, and a two-dimensional (2D) model collinear reaction using an anharmonic potential energy surface. Approximate treatments of the quantum hydrodynamic equations will also be discussed which could allow scaling of the calculations to hundreds of degrees of freedom which is important for treating tunneling phenomena in condensed phase systems.

•Jun 15, 2012•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

Quantum wave packet dynamics with trajectories: Application to reactive scattering

Róbert E. Wyatt

The quantum trajectory method (QTM) for time-dependent wave packet dynamics involves integration of the de Broglie–Bohm hydrodynamic equations for the evolving probability fluid [C. Lopreore and R. E. Wyatt, Phys. Rev. Lett. 82, 5190 (1999)]. The equations of motion for discretized elements of the probability fluid (particles) are integrated in the Lagrangian, moving with the fluid, picture. These fluid elements move under the influence of both the usual potential energy function and the quantum potential, which involves the curvature of the quantum amplitude. The quantum potential and the quantum force are evaluated using a moving weighted least squares algorithm. As a demonstration of applicability, the QTM is applied to a model collinear reaction with an activation barrier. The reaction probabilities are in good agreement with exact quantum results, even for a relatively small number of particles in the ensemble. The QTM accurately describes tunneling using only real valued trajectories. In addition to the reaction probability, plots are presented to show the probability density and the flux distributions at several time steps during the reaction.

The Journal of Chemical Physics•Sep 8, 1999•135 citations
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.
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.
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 quantum hydrodynamic approach is based on the de Broglie-Bohm formulation of quantum mechanics. The resulting fluid-like equations of motion describe the flow of probability and an accurate solution to these equations is equivalent to solving the time-dependent Schroedinger equation. Furthermore, the hydrodynamic approach provides new insight into the mechanisms as well as an alternative computational approach for treating tunneling phenomena. New concepts include well-defined 'quantum trajectories', 'quantum potential', and 'quantum force' all of which have classical analogues. The quantum potential and its associated force give rise to all quantum mechanical effects such as zero point energy, tunneling, and interference. A new numerical approach called the Iterative Finite Difference Method (IFDM) will be discussed. The IFDM is used to solve the set of non-linear coupled hydrodynamic equations. It is 2nd-order accurate in both space and time and exhibits exponential convergence with respect to the iteration count. The stability and computational efficiency of the IFDM is significantly improved by using a 'smart' Eulerian grid which has the same computational advantages as a Lagrangian or Arbitrary Lagrangian Eulerian (ALE) grid. The IFDM is also capable of treating anharmonic potentials. Example calculations using the IFDM will be presented which include: a one-dimensional Gaussian wave packet tunneling through an Eckart barrier, a one-dimensional bound-state Morse oscillator, and a two-dimensional (2D) model collinear reaction using an anharmonic potential energy surface. Approximate treatments of the quantum hydrodynamic equations will also be discussed which could allow scaling of the calculations to hundreds of degrees of freedom which is important for treating tunneling phenomena in condensed phase systems.
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 quantum trajectory method (QTM) for time-dependent wave packet dynamics involves integration of the de Broglie–Bohm hydrodynamic equations for the evolving probability fluid [C. Lopreore and R. E. Wyatt, Phys. Rev. Lett. 82, 5190 (1999)]. The equations of motion for discretized elements of the probability fluid (particles) are integrated in the Lagrangian, moving with the fluid, picture. These fluid elements move under the influence of both the usual potential energy function and the quantum potential, which involves the curvature of the quantum amplitude. The quantum potential and the quantum force are evaluated using a moving weighted least squares algorithm. As a demonstration of applicability, the QTM is applied to a model collinear reaction with an activation barrier. The reaction probabilities are in good agreement with exact quantum results, even for a relatively small number of particles in the ensemble. The QTM accurately describes tunneling using only real valued trajectories. In addition to the reaction probability, plots are presented to show the probability density and the flux distributions at several time steps during the reaction.