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

Wave packet dynamics on adaptive moving grids

Róbert E. Wyatt

Published Nov 14, 2002• 46 citations
Publisher

Abstract

A new hybrid moving grid approach to wave packet dynamics is described. Exterior points within a nonrigid frame follow Lagrangian trajectories obtained by solving the hydrodynamic equations of motion. Internal grid points within one or more open windows follow non-Lagrangian adaptive paths. Within these windows, problems encountered with quantum trajectories near wave function nodes are circumvented by directly solving the moving path transform of the Schrödinger equation. Excellent results are obtained for evolution of the density in a double well potential even though multiple ripples develop in the density.

DOI
10.1063/1.1517045

Related papers

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

Quantum wave packet dynamics with trajectories: Implementation with distributed approximating functionals

Róbert E. Wyatt, Donald Jack Kouri, David K. Hoffman

The quantum trajectory method (QTM) 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 N fluid elements (particles) moving under the influence of both the force from the potential surface and from the quantum potential. In this study, distributed approximating functionals (DAFs) are used on a uniform grid to compute the necessary derivatives in the equations of motion. Transformations between the physical grid where the particle coordinates are defined and the uniform grid are handled through a Jacobian, which is also computed using DAFs. A difficult problem associated with computing derivatives on finite grids is the edge problem. This is handled effectively by using DAFs within a least squares approach to extrapolate from the known function region into the neighboring regions. The QTM–DAF is then applied to wave packet transmission through a one-dimensional Eckart potential. Emphasis is placed upon computation of the transmitted density and wave function. A problem that develops when part of the wave packet reflects back into the reactant region is avoided in this study by introducing a potential ramp to sweep the reflected particles away from the barrier region.

The Journal of Chemical Physics•Jun 22, 2000•54 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

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

Electronic transitions with quantum trajectories. II

Courtney L. Lopreore, Róbert E. Wyatt

The quantum trajectory method (QTM) is applied to nonadiabatic electronic transitions. Equations of motion in a Lagrangian framework are derived for the probability density, velocity, position, and action functions for a discretized wave packet moving along coupled potential energy surfaces. In solving these equations of motion, we obtain agreement with exact quantum results computed by solving the time-dependent Schrödinger equation on a space-fixed grid. On each of the coupled potential energy surfaces, the dynamics of the trajectories is fully quantum mechanical, i.e., there are no “surface–hopping transitions.” We present a detailed analysis of the QTM results including density changes, complex oscillations of the wave functions, phase space analysis, and a detailed discussion of the forces that contribute to the evolution the trajectories.

The Journal of Chemical Physics•Jan 22, 2002•78 citations

A new method for solving the quantum hydrodynamic equations of motion

Brian K. Kendrick

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 Journal of Chemical Physics•Sep 8, 2003•100 citations
A new hybrid moving grid approach to wave packet dynamics is described. Exterior points within a nonrigid frame follow Lagrangian trajectories obtained by solving the hydrodynamic equations of motion. Internal grid points within one or more open windows follow non-Lagrangian adaptive paths. Within these windows, problems encountered with quantum trajectories near wave function nodes are circumvented by directly solving the moving path transform of the Schrödinger equation. Excellent results are obtained for evolution of the density in a double well potential even though multiple ripples develop in the density.
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 quantum trajectory method (QTM) 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 N fluid elements (particles) moving under the influence of both the force from the potential surface and from the quantum potential. In this study, distributed approximating functionals (DAFs) are used on a uniform grid to compute the necessary derivatives in the equations of motion. Transformations between the physical grid where the particle coordinates are defined and the uniform grid are handled through a Jacobian, which is also computed using DAFs. A difficult problem associated with computing derivatives on finite grids is the edge problem. This is handled effectively by using DAFs within a least squares approach to extrapolate from the known function region into the neighboring regions. The QTM–DAF is then applied to wave packet transmission through a one-dimensional Eckart potential. Emphasis is placed upon computation of the transmitted density and wave function. A problem that develops when part of the wave packet reflects back into the reactant region is avoided in this study by introducing a potential ramp to sweep the reflected particles away from the barrier region.
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.
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 quantum trajectory method (QTM) is applied to nonadiabatic electronic transitions. Equations of motion in a Lagrangian framework are derived for the probability density, velocity, position, and action functions for a discretized wave packet moving along coupled potential energy surfaces. In solving these equations of motion, we obtain agreement with exact quantum results computed by solving the time-dependent Schrödinger equation on a space-fixed grid. On each of the coupled potential energy surfaces, the dynamics of the trajectories is fully quantum mechanical, i.e., there are no “surface–hopping transitions.” We present a detailed analysis of the QTM results including density changes, complex oscillations of the wave functions, phase space analysis, and a detailed discussion of the forces that contribute to the evolution the trajectories.
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.