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 new method for solving the quantum hydrodynamic equations of motion

Brian K. Kendrick

Published Sep 8, 2003• 100 citations
Publisher

Abstract

The quantum hydrodynamic equations associated with the de Broglie–Bohm formulation of quantum mechanics are solved using a meshless method based on a moving least squares approach. An arbitrary Lagrangian–Eulerian frame of reference is used which significantly improves the accuracy and stability of the method when compared to an approach based on a purely Lagrangian frame of reference. A regridding algorithm is implemented which adds and deletes points when necessary in order to maintain accurate and stable calculations. It is shown that unitarity in the time evolution of the quantum wave packet is significantly improved by propagating using averaged fields. As nodes in the reflected wave packet start to form, the quantum potential and force become very large and numerical instabilities occur. By introducing artificial viscosity into the equations of motion, these instabilities can be avoided and the stable propagation of the wave packet for very long times becomes possible. Results are presented for the scattering of a wave packet from a repulsive Eckart barrier.

DOI
10.1063/1.1601217

Related papers

A new method for solving the quantum hydrodynamic equations of motion: Application to two-dimensional reactive scattering

Denise K. Pauler, Brian K. Kendrick

The de Broglie-Bohm hydrodynamic equations of motion are solved using a meshless method based on a moving least squares approach and an arbitrary Lagrangian-Eulerian frame of reference. A regridding algorithm adds and deletes computational points as needed in order to maintain a uniform interparticle spacing, and unitary time evolution is obtained by propagating the wave packet using averaged fields. The numerical instabilities associated with the formation of nodes in the reflected portion of the wave packet are avoided by adding artificial viscosity to the equations of motion. The methodology is applied to a two-dimensional model collinear reaction with an activation barrier. Reaction probabilities are computed as a function of both time and energy, and are in excellent agreement with those based on the quantum trajectory method.

The Journal of Chemical Physics•Dec 24, 2003•52 citations

Quantum hydrodynamics: Capturing a reactive scattering resonance

Sean Derrickson, Eric R. Bittner, Brian K. Kendrick

The hydrodynamic equations of motion associated with the de Broglie-Bohm formulation of quantum mechanics are solved using a meshless method based upon a moving least-squares approach. An arbitrary Lagrangian-Eulerian frame of reference and a regridding algorithm which adds and deletes computational points are used to maintain a uniform and nearly constant interparticle spacing. The methodology also uses averaged fields to maintain unitary time evolution. The numerical instabilities associated with the formation of nodes in the reflected portion of the wave packet are avoided by adding artificial viscosity to the equations of motion. A new and more robust artificial viscosity algorithm is presented which gives accurate scattering results and is capable of capturing quantum resonances. The methodology is applied to a one-dimensional model chemical reaction that is known to exhibit a quantum resonance. The correlation function approach is used to compute the reactive scattering matrix, reaction probability, and time delay as a function of energy. Excellent agreement is obtained between the scattering results based upon the quantum hydrodynamic approach and those based upon standard quantum mechanics. This is the first clear demonstration of the ability of moving grid approaches to accurately and robustly reproduce resonance structures in a scattering system.

The Journal of Chemical Physics•Aug 1, 2005•20 citations

Quantum Wave Packet Dynamics with Trajectories

Courtney L. Lopreore, Róbert E. Wyatt

The de Broglie--Bohm causual (hydrodynamic) formulation of quantum mechanics is computationally implemented in the Lagrangian (moving with the fluid) viewpoint. The quantum potential and force are accurately evaluated with a moving weighted least squares algorithm. The quantum trajectory method is then applied to barrier tunneling on smooth potential surfaces. Analysis of the tunneling mechanism leads to a novel and accurate approximation: shortly after the wave packet is launched, completely neglect all quantum terms in the dynamical equations for motion along the tunneling coordinate.

Physical Review Letters•Jun 28, 1999•441 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

Quantum wave packet dynamics with trajectories: Implementation with adaptive Lagrangian grids

Róbert E. Wyatt, Eric R. Bittner

The quantum trajectory method was recently developed to solve the hydrodynamic equations of motion in the Lagrangian, moving-with-the-fluid, picture. In this approach, trajectories are integrated for fluid elements (“particles”) moving under the influence of the combined force from the potential surface and the quantum potential. To accurately compute the quantum potential and the quantum force, it is necessary to obtain the derivatives of a function given only the values on the unstructured mesh defined by the particle locations. However, in some regions of space–time, the particle mesh shows compression and inflation associated with regions of large and small density, respectively. Inflation is especially severe near nodes in the wave function. In order to circumvent problems associated with highly nonuniform grids defined by the particle locations, adaptation of moving grids is introduced in this study. By changing the representation of the wave function in these local regions (which can be identified by diagnostic tools), propagation is possible to much longer times. These grid adaptation techniques are applied to the reflected portion of a wave packet scattering from an Eckart potential.

The Journal of Chemical Physics•Nov 22, 2000•122 citations

An arbitrary Lagrangian–Eulerian approach to solving the quantum hydrodynamic equations of motion: Equidistribution with “smart” springs

Corey Jason Trahan, Róbert E. Wyatt

Recently, the quantum trajectory method (QTM) has been utilized in solving several quantum mechanical wave packet scattering problems including barrier transmission and electronic nonadiabatic dynamics. By propagating the real-valued action and amplitude functions in the Lagrangian frame, only a fraction of the grid points needed for Eulerian fixed-grid methods are used while still obtaining accurate solutions. Difficulties arise, however, near wave function nodes and in regions of sharp oscillatory features, and because of this many quantum mechanical problems have not yet been amenable to solution with the QTM. This study proposes a hybrid of both the Lagrangian and Eulerian techniques in what is termed the arbitrary Lagrangian–Eulerian method (ALE). In the ALE method, an additional equation of motion governing the momentum of the grid points is coupled into the quantum hydrodynamic equations. These new “quasi-” Bohmian trajectories can be dynamically adapted to the emergent features of the time evolving hydrodynamic fields and are non-Lagrangian. In this study it is shown that the ALE method applied to an uphill ramp potential that was previously unsolvable by the current Lagrangian QTM not only yields stable transmission probabilities with accuracies comparable to that of a high resolution Eulerian method, but does so with a small number of grid points and for extremely long propagation times. To determine the grid point positions at each new time, an equidistribution method is used that is constructed similar to the stiffness matrix of a classical spring system in equilibrium. Each “smart” spring is dependent on a local function M(x) called the monitor function which can sense gradients or curvatures of the fields surrounding its position. To constrain grid points from having zero separation and possible overlap, a new system of equations is derived that includes a minimum separation parameter which prevents this from occurring.

The Journal of Chemical Physics•Feb 28, 2003•58 citations
The quantum hydrodynamic equations associated with the de Broglie–Bohm formulation of quantum mechanics are solved using a meshless method based on a moving least squares approach. An arbitrary Lagrangian–Eulerian frame of reference is used which significantly improves the accuracy and stability of the method when compared to an approach based on a purely Lagrangian frame of reference. A regridding algorithm is implemented which adds and deletes points when necessary in order to maintain accurate and stable calculations. It is shown that unitarity in the time evolution of the quantum wave packet is significantly improved by propagating using averaged fields. As nodes in the reflected wave packet start to form, the quantum potential and force become very large and numerical instabilities occur. By introducing artificial viscosity into the equations of motion, these instabilities can be avoided and the stable propagation of the wave packet for very long times becomes possible. Results are presented for the scattering of a wave packet from a repulsive Eckart barrier.
The de Broglie-Bohm hydrodynamic equations of motion are solved using a meshless method based on a moving least squares approach and an arbitrary Lagrangian-Eulerian frame of reference. A regridding algorithm adds and deletes computational points as needed in order to maintain a uniform interparticle spacing, and unitary time evolution is obtained by propagating the wave packet using averaged fields. The numerical instabilities associated with the formation of nodes in the reflected portion of the wave packet are avoided by adding artificial viscosity to the equations of motion. The methodology is applied to a two-dimensional model collinear reaction with an activation barrier. Reaction probabilities are computed as a function of both time and energy, and are in excellent agreement with those based on the quantum trajectory method.
The hydrodynamic equations of motion associated with the de Broglie-Bohm formulation of quantum mechanics are solved using a meshless method based upon a moving least-squares approach. An arbitrary Lagrangian-Eulerian frame of reference and a regridding algorithm which adds and deletes computational points are used to maintain a uniform and nearly constant interparticle spacing. The methodology also uses averaged fields to maintain unitary time evolution. The numerical instabilities associated with the formation of nodes in the reflected portion of the wave packet are avoided by adding artificial viscosity to the equations of motion. A new and more robust artificial viscosity algorithm is presented which gives accurate scattering results and is capable of capturing quantum resonances. The methodology is applied to a one-dimensional model chemical reaction that is known to exhibit a quantum resonance. The correlation function approach is used to compute the reactive scattering matrix, reaction probability, and time delay as a function of energy. Excellent agreement is obtained between the scattering results based upon the quantum hydrodynamic approach and those based upon standard quantum mechanics. This is the first clear demonstration of the ability of moving grid approaches to accurately and robustly reproduce resonance structures in a scattering system.
The de Broglie--Bohm causual (hydrodynamic) formulation of quantum mechanics is computationally implemented in the Lagrangian (moving with the fluid) viewpoint. The quantum potential and force are accurately evaluated with a moving weighted least squares algorithm. The quantum trajectory method is then applied to barrier tunneling on smooth potential surfaces. Analysis of the tunneling mechanism leads to a novel and accurate approximation: shortly after the wave packet is launched, completely neglect all quantum terms in the dynamical equations for motion along the tunneling coordinate.
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.
The quantum trajectory method was recently developed to solve the hydrodynamic equations of motion in the Lagrangian, moving-with-the-fluid, picture. In this approach, trajectories are integrated for fluid elements (“particles”) moving under the influence of the combined force from the potential surface and the quantum potential. To accurately compute the quantum potential and the quantum force, it is necessary to obtain the derivatives of a function given only the values on the unstructured mesh defined by the particle locations. However, in some regions of space–time, the particle mesh shows compression and inflation associated with regions of large and small density, respectively. Inflation is especially severe near nodes in the wave function. In order to circumvent problems associated with highly nonuniform grids defined by the particle locations, adaptation of moving grids is introduced in this study. By changing the representation of the wave function in these local regions (which can be identified by diagnostic tools), propagation is possible to much longer times. These grid adaptation techniques are applied to the reflected portion of a wave packet scattering from an Eckart potential.
Recently, the quantum trajectory method (QTM) has been utilized in solving several quantum mechanical wave packet scattering problems including barrier transmission and electronic nonadiabatic dynamics. By propagating the real-valued action and amplitude functions in the Lagrangian frame, only a fraction of the grid points needed for Eulerian fixed-grid methods are used while still obtaining accurate solutions. Difficulties arise, however, near wave function nodes and in regions of sharp oscillatory features, and because of this many quantum mechanical problems have not yet been amenable to solution with the QTM. This study proposes a hybrid of both the Lagrangian and Eulerian techniques in what is termed the arbitrary Lagrangian–Eulerian method (ALE). In the ALE method, an additional equation of motion governing the momentum of the grid points is coupled into the quantum hydrodynamic equations. These new “quasi-” Bohmian trajectories can be dynamically adapted to the emergent features of the time evolving hydrodynamic fields and are non-Lagrangian. In this study it is shown that the ALE method applied to an uphill ramp potential that was previously unsolvable by the current Lagrangian QTM not only yields stable transmission probabilities with accuracies comparable to that of a high resolution Eulerian method, but does so with a small number of grid points and for extremely long propagation times. To determine the grid point positions at each new time, an equidistribution method is used that is constructed similar to the stiffness matrix of a classical spring system in equilibrium. Each “smart” spring is dependent on a local function M(x) called the monitor function which can sense gradients or curvatures of the fields surrounding its position. To constrain grid points from having zero separation and possible overlap, a new system of equations is derived that includes a minimum separation parameter which prevents this from occurring.