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

Trajectory integration of the quantum hydrodynamic equations of motion

Corey Jason Trahan

Published Jan 1, 2003• 0 citations
OpenAlex

Abstract

Recently, in an effort to solve more realistic problems in quantum dynamics, much attention has been directed into numerically integrating the quantum hydrodynamic equations of motions (QHEM), as opposed to directly solving the time-dependent Schrödinger equation (TDSE). Such efforts have been provoked by the many numerical drawbacks encountered when solving the TDSE on a fixed-grid. In this dissertation, one trajectory method for integrating the QHEM is reviewed, and two novel trajectories methods are described. The first of these, the quantum trajectory method (QTM), was introduced in 1999 and has been used to solve many problems in quantum dynamics since then. However, severe numerical problems are encountered when this method is applied to problems that form wave function nodes. To get around this problem, new methods for numerically integrating the QHEM are needed. In the first novel method described, the arbitrary Lagrangian-Eulerian (ALE) method, particle trajectories are governed by a predetermined equation of motion that is user-supplied. The ALE method remedies inflation and compression problems encountered in the pure Lagrangian QTM. In the second new method discussed, the derivative propagating method (DPM), single quantum trajectories can be calculated one at a time, as opposed to the ensemble propagation of the QTM and ALE method. Using these two methods, new solutions to the QHEM are obtained where the QTM fails. In addition to solving the QHEM, the DPM is also used to solve the classical Klein-Kramers equation in this dissertation. This equation governs the Markovian phase space evolution of a system coupled to an environment such as a heat bath. This marks the first time single trajectories have been used to solve both the QHEM and the Klein-Kramers equations.

Related papers

Multidimensional quantum trajectories: Applications of the derivative propagation method

Corey Jason Trahan, Róbert E. Wyatt, Bill Poirier

In a previous publication [J. Chem. Phys. 118, 9911 (2003)], the derivative propagation method (DPM) was introduced as a novel numerical scheme for solving the quantum hydrodynamic equations of motion (QHEM) and computing the time evolution of quantum mechanical wave packets. These equations are a set of coupled, nonlinear partial differential equations governing the time evolution of the real-valued functions C and S in the complex action, S=C(r,t) + iS(r,t)/Planck's over 2pi, where Psi(r,t)=exp(S). Past numerical solutions to the QHEM were obtained via ensemble trajectory propagation, where the required first- and second-order spatial derivatives were evaluated using fitting techniques such as moving least squares. In the DPM, however, equations of motion are developed for the derivatives themselves, and a truncated set of these are integrated along quantum trajectories concurrently with the original QHEM equations for C and S. Using the DPM quantum effects can be included at various orders of approximation; no spatial fitting is involved; there is no basis set expansion; and single, uncoupled quantum trajectories can be propagated (in parallel) rather than in correlated ensembles. In this study, the DPM is extended from previous one-dimensional (1D) results to calculate transmission probabilities for 2D and 3D wave packet evolution on coupled Eckart barrier/harmonic oscillator surfaces. In the 2D problem, the DPM results are compared to standard numerical integration of the time-dependent Schrodinger equation. Also in this study, the practicality of implementing the DPM for systems with many more degrees of freedom is discussed.

The Journal of Chemical Physics•Apr 22, 2005•40 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

Evolution of classical and quantum phase-space distributions: A new trajectory approach for phase space hydrodynamics

Corey Jason Trahan, Róbert E. Wyatt

Recently, Donoso and Martens described a method for evolving both classical and quantum phase-space distribution functions, W(q,p,t), that involves the propagation of an ensemble of correlated trajectories. The trajectories are linked into a unified whole by spatial and momentum derivatives of density dependent terms in the equations of motion. On each time step, these nonlocal terms were evaluated by fitting the density around each trajectory to an assumed functional form. In the present study, we develop a different trajectory method for propagating phase-space distribution functions. A hierarchy of coupled analytic equations of motion are derived for the q and p derivatives of the density and a truncated set of these are integrated along each trajectory concurrently with the equation of motion for the density. The advantage of this approach is that individual trajectories can be propagated, one at a time, and function fitting is not required to evaluate the nonlocal terms. Regional nonlocality can be incorporated at various levels of approximation to “dress” what would otherwise be “thin” locally propagating trajectories. This derivative propagation method is used to obtain trajectory solutions for the Klein–Kramers equation, the Husimi equation, and for a smoothed version of the Caldeira–Leggett equation derived by the Diosi. Trajectory solutions are obtained for the relaxation of an oscillator in contact with a thermal bath and for the decay of a metastable state.

The Journal of Chemical Physics•Oct 3, 2003•43 citations

Electronic transitions with quantum trajectories

Róbert E. Wyatt, Courtney L. Lopreore, Gérard Parlant

The quantum trajectory method (QTM) is extended to the dynamics of electronic nonadiabiatic collisions. Equations of motion are first derived for the probability density, velocity, and action function for wave packets moving on each of the coupled electronic potential surfaces. These discretized equations are solved in the Lagrangian (moving with the fluid) picture to give the trajectory dynamics of fluid elements evolving on each potential surface. This trajectory method is fully quantum mechanical and does not involve “trajectory surface hopping.” The method is applied to nonadiabiatic collision models involving two coupled electronic states. The quantum trajectory results are in excellent agreement with solutions computed (using space-fixed grid methods) directly from the time-dependent Schrödinger equation.

The Journal of Chemical Physics•Mar 22, 2001•112 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
Recently, in an effort to solve more realistic problems in quantum dynamics, much attention has been directed into numerically integrating the quantum hydrodynamic equations of motions (QHEM), as opposed to directly solving the time-dependent Schrödinger equation (TDSE). Such efforts have been provoked by the many numerical drawbacks encountered when solving the TDSE on a fixed-grid. In this dissertation, one trajectory method for integrating the QHEM is reviewed, and two novel trajectories methods are described. The first of these, the quantum trajectory method (QTM), was introduced in 1999 and has been used to solve many problems in quantum dynamics since then. However, severe numerical problems are encountered when this method is applied to problems that form wave function nodes. To get around this problem, new methods for numerically integrating the QHEM are needed. In the first novel method described, the arbitrary Lagrangian-Eulerian (ALE) method, particle trajectories are governed by a predetermined equation of motion that is user-supplied. The ALE method remedies inflation and compression problems encountered in the pure Lagrangian QTM. In the second new method discussed, the derivative propagating method (DPM), single quantum trajectories can be calculated one at a time, as opposed to the ensemble propagation of the QTM and ALE method. Using these two methods, new solutions to the QHEM are obtained where the QTM fails. In addition to solving the QHEM, the DPM is also used to solve the classical Klein-Kramers equation in this dissertation. This equation governs the Markovian phase space evolution of a system coupled to an environment such as a heat bath. This marks the first time single trajectories have been used to solve both the QHEM and the Klein-Kramers equations.
In a previous publication [J. Chem. Phys. 118, 9911 (2003)], the derivative propagation method (DPM) was introduced as a novel numerical scheme for solving the quantum hydrodynamic equations of motion (QHEM) and computing the time evolution of quantum mechanical wave packets. These equations are a set of coupled, nonlinear partial differential equations governing the time evolution of the real-valued functions C and S in the complex action, S=C(r,t) + iS(r,t)/Planck's over 2pi, where Psi(r,t)=exp(S). Past numerical solutions to the QHEM were obtained via ensemble trajectory propagation, where the required first- and second-order spatial derivatives were evaluated using fitting techniques such as moving least squares. In the DPM, however, equations of motion are developed for the derivatives themselves, and a truncated set of these are integrated along quantum trajectories concurrently with the original QHEM equations for C and S. Using the DPM quantum effects can be included at various orders of approximation; no spatial fitting is involved; there is no basis set expansion; and single, uncoupled quantum trajectories can be propagated (in parallel) rather than in correlated ensembles. In this study, the DPM is extended from previous one-dimensional (1D) results to calculate transmission probabilities for 2D and 3D wave packet evolution on coupled Eckart barrier/harmonic oscillator surfaces. In the 2D problem, the DPM results are compared to standard numerical integration of the time-dependent Schrodinger equation. Also in this study, the practicality of implementing the DPM for systems with many more degrees of freedom is discussed.
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.
Recently, Donoso and Martens described a method for evolving both classical and quantum phase-space distribution functions, W(q,p,t), that involves the propagation of an ensemble of correlated trajectories. The trajectories are linked into a unified whole by spatial and momentum derivatives of density dependent terms in the equations of motion. On each time step, these nonlocal terms were evaluated by fitting the density around each trajectory to an assumed functional form. In the present study, we develop a different trajectory method for propagating phase-space distribution functions. A hierarchy of coupled analytic equations of motion are derived for the q and p derivatives of the density and a truncated set of these are integrated along each trajectory concurrently with the equation of motion for the density. The advantage of this approach is that individual trajectories can be propagated, one at a time, and function fitting is not required to evaluate the nonlocal terms. Regional nonlocality can be incorporated at various levels of approximation to “dress” what would otherwise be “thin” locally propagating trajectories. This derivative propagation method is used to obtain trajectory solutions for the Klein–Kramers equation, the Husimi equation, and for a smoothed version of the Caldeira–Leggett equation derived by the Diosi. Trajectory solutions are obtained for the relaxation of an oscillator in contact with a thermal bath and for the decay of a metastable state.
The quantum trajectory method (QTM) is extended to the dynamics of electronic nonadiabiatic collisions. Equations of motion are first derived for the probability density, velocity, and action function for wave packets moving on each of the coupled electronic potential surfaces. These discretized equations are solved in the Lagrangian (moving with the fluid) picture to give the trajectory dynamics of fluid elements evolving on each potential surface. This trajectory method is fully quantum mechanical and does not involve “trajectory surface hopping.” The method is applied to nonadiabiatic collision models involving two coupled electronic states. The quantum trajectory results are in excellent agreement with solutions computed (using space-fixed grid methods) directly from the time-dependent Schrödinger equation.
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.