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

Quantum wave packet dynamics with trajectories: Application to reactive scattering

Róbert E. Wyatt

Published Sep 8, 1999• 135 citations
Publisher

Abstract

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.

DOI
10.1063/1.479205

Related papers

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

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

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

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

Quantum tunneling dynamics using hydrodynamic trajectories

Eric R. Bittner

In this paper we compute quantum trajectories arising from Bohm’s causal description of quantum mechanics. Our computational methodology is based upon a finite-element moving least-squares method (MWLS) presented recently by Wyatt and co-workers [Lopreore and Wyatt, Phys. Rev. Lett. 82, 5190 (1999)]. This method treats the “particles” in the quantum Hamilton–Jacobi equation as Lagrangian fluid elements that carry the phase, S, and density, ρ, required to reconstruct the quantum wave function. Here, we compare results obtained via the MWLS procedure to exact results obtained either analytically or by numerical solution of the time-dependent Schrödinger equation. Two systems are considered: first, dynamics in a harmonic well and second, tunneling dynamics in a double well potential. In the case of tunneling in the double well potential, the quantum potential acts to lower the barrier, separating the right- and left-hand sides of the well, permitting trajectories to pass from one side to another. However, as probability density passes from one side to the other, the effective barrier begins to rise and eventually will segregate trajectories in one side from the other. We note that the MWLS trajectories exhibited long time stability in the purely harmonic cases. However, this stability was not evident in the barrier crossing dynamics. Comparisons to exact trajectories obtained via wave packet calculations indicate that the MWLS trajectories tend to underestimate the effects of constructive and destructive interference effects.

arXiv:quant-ph/0001119•The Journal of Chemical Physics•Jun 8, 2000•92 citations
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 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 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 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 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) 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.
In this paper we compute quantum trajectories arising from Bohm’s causal description of quantum mechanics. Our computational methodology is based upon a finite-element moving least-squares method (MWLS) presented recently by Wyatt and co-workers [Lopreore and Wyatt, Phys. Rev. Lett. 82, 5190 (1999)]. This method treats the “particles” in the quantum Hamilton–Jacobi equation as Lagrangian fluid elements that carry the phase, S, and density, ρ, required to reconstruct the quantum wave function. Here, we compare results obtained via the MWLS procedure to exact results obtained either analytically or by numerical solution of the time-dependent Schrödinger equation. Two systems are considered: first, dynamics in a harmonic well and second, tunneling dynamics in a double well potential. In the case of tunneling in the double well potential, the quantum potential acts to lower the barrier, separating the right- and left-hand sides of the well, permitting trajectories to pass from one side to another. However, as probability density passes from one side to the other, the effective barrier begins to rise and eventually will segregate trajectories in one side from the other. We note that the MWLS trajectories exhibited long time stability in the purely harmonic cases. However, this stability was not evident in the barrier crossing dynamics. Comparisons to exact trajectories obtained via wave packet calculations indicate that the MWLS trajectories tend to underestimate the effects of constructive and destructive interference effects.