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

International Journal of Quantum Chemistry

Particle method for the numerical solution of the time-dependent schrödinger equation

Jerome H. Weiner, Attila Aşkar

Published Jun 18, 2009• 0 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.

DOI
10.1002/qua.560050830

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. 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

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

Computing the wavefunction from trajectories: particle and wave pictures in quantum mechanics and their relation

Peter Holland

We discuss the particle method in quantum mechanics which provides an exact scheme to calculate the time-dependent wavefunction from a single-valued continuum of trajectories where two spacetime points are linked by at most a single orbit. A natural language for the theory is offered by the hydrodynamic analogy, in which wave mechanics corresponds to the Eulerian picture and the particle theory to the Lagrangian picture. The Lagrangian model for the quantum fluid may be developed from a variational principle. The Euler-Lagrange equations imply a fourth-order nonlinear partial differential equation to calculate the trajectories of the fluid particles as functions of their initial coordinates using as input the initial wavefunction. The admissible solutions are those consistent with quasi-potential flow. The effect of the superposition principle is represented via a nonclassical force on each particle. The wavefunction is computed via the standard map between the Lagrangian coordinates and the Eulerian fields, which provides the analogue in this model of Huygens principle in wave mechanics. The method is illustrated by calculating the time-dependence of a free Gaussian wavefunction. The Eulerian and Lagrangian pictures are complementary descriptions of a quantum process in that they have associated Hamiltonian formulations that are connected by a canonical transformation. The de Broglie-Bohm interpretation, which employs the same set of trajectories, should not be conflated with the Lagrangian version of the hydrodynamic interpretation. The theory implies that the mathematical results of the de Broglie-Bohm model may be regarded as statements about quantum mechanics itself rather than about its interpretation.

•Jan 1, 2004•51 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

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

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
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.
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.
We discuss the particle method in quantum mechanics which provides an exact scheme to calculate the time-dependent wavefunction from a single-valued continuum of trajectories where two spacetime points are linked by at most a single orbit. A natural language for the theory is offered by the hydrodynamic analogy, in which wave mechanics corresponds to the Eulerian picture and the particle theory to the Lagrangian picture. The Lagrangian model for the quantum fluid may be developed from a variational principle. The Euler-Lagrange equations imply a fourth-order nonlinear partial differential equation to calculate the trajectories of the fluid particles as functions of their initial coordinates using as input the initial wavefunction. The admissible solutions are those consistent with quasi-potential flow. The effect of the superposition principle is represented via a nonclassical force on each particle. The wavefunction is computed via the standard map between the Lagrangian coordinates and the Eulerian fields, which provides the analogue in this model of Huygens principle in wave mechanics. The method is illustrated by calculating the time-dependence of a free Gaussian wavefunction. The Eulerian and Lagrangian pictures are complementary descriptions of a quantum process in that they have associated Hamiltonian formulations that are connected by a canonical transformation. The de Broglie-Bohm interpretation, which employs the same set of trajectories, should not be conflated with the Lagrangian version of the hydrodynamic interpretation. The theory implies that the mathematical results of the de Broglie-Bohm model may be regarded as statements about quantum mechanics itself rather than about its interpretation.
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.
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.
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.