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

Coping with the node problem in quantum hydrodynamics: The covering function method

Dmytro Babyuk, Róbert E. Wyatt

Published Nov 4, 2004• 51 citations
Publisher

Abstract

A conceptually simple approach, the covering function method (CFM), is developed to cope with the node problem in the hydrodynamic formulation of quantum mechanics. As nodes begin to form in a scattering wave packet (detected by a monitor function), a nodeless covering wave function is added to it yielding a total function that is also nodeless. Both local and global choices for the covering function are described. The total and covering functions are then propagated separately in the hydrodynamic picture. At a later time, the actual wave function is recovered from the two propagated functions. The results obtained for Eckart barrier scattering in one dimension are in excellent agreement with exact results, even for very long propagation times t=1.2 ps. The capability of the CFM is also demonstrated for multidimensional propagation of a vibrationally excited wave packet.

DOI
10.1063/1.1801211

Related papers

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 hydrodynamics: Application to N-dimensional reactive scattering

Brian K. Kendrick

The quantum hydrodynamic equations associated with the de Broglie-Bohm formulation of quantum mechanics are solved using a new methodology which gives an accurate, unitary, and stable propagation of a time dependent quantum wave packet [B. K. Kendrick, J. Chem. Phys. 119, 5805 (2003)]. The methodology is applied to an N-dimensional model chemical reaction with an activation barrier. A parallel version of the methodology is presented which is designed to run on massively parallel supercomputers. The computational scaling properties of the parallel code are investigated both as a function of the number of processors and the dimension N. A decoupling scheme is introduced which decouples the multidimensional quantum hydrodynamic equations into a set of uncoupled one-dimensional problems. The decoupling scheme dramatically reduces the computation time and is highly parallelizable. Furthermore, the computation time is shown to scale linearly with respect to the dimension N=2,...,100.

The Journal of Chemical Physics•Jul 23, 2004•39 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

Scattered-wave-packet formalism with applications to barrier scattering and quantum transistors

Chia‐Chun Chou, Róbert E. Wyatt

The scattered wave formalism developed for a quantum subsystem interacting with reservoirs through open boundaries is applied to one- or two-dimensional barrier scattering and quantum transistors. The total wave function is divided into incident and scattered components. Markovian outgoing wave boundary conditions are imposed on the scattered or total wave function by either the ratio or polynomial methods. For barrier scattering problems, accurate time-dependent transmission probabilities are obtained through the integration of the modified time-dependent Schrödinger equations for the scattered wave function. For quantum transistors, the time-dependent transport is studied for a quantum wave packet propagating through the conduction channel of a field effect transistor. This study shows that the scattered wave formalism significantly reduces computational effort relative to other open boundary methods and demonstrates wide applications to quantum dynamical processes.

Physical Review E•Nov 23, 2011•3 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

A new method for wave packet dynamics: Derivative propagation along quantum trajectories

Corey Jason Trahan, Keith H. Hughes, Róbert E. Wyatt

A new method is proposed for computing the time evolution of quantum mechanical wave packets. Equations of motion for the real-valued functions C and S in the complex action S=C(r,t)+iS(r,t)/ℏ, with ψ(r,t)=exp(S), involve gradients and curvatures of C and S. In previous implementations of the hydrodynamic formulation, various time-consuming fitting techniques of limited accuracy were used to evaluate these derivatives around each fluid element in an evolving ensemble. In this study, equations of motion are developed for the spatial derivatives themselves and a small set of these are integrated along quantum trajectories concurrently with the equations for C and S. Significantly, quantum effects can be included at various orders of approximation, no spatial fitting is involved, there are no basis set expansions, and single quantum trajectories (rather than correlated ensembles) may be propagated, one at a time. Excellent results are obtained when the derivative propagation method is applied to anharmonic potentials involving barrier transmission.

The Journal of Chemical Physics•May 22, 2003•85 citations
A conceptually simple approach, the covering function method (CFM), is developed to cope with the node problem in the hydrodynamic formulation of quantum mechanics. As nodes begin to form in a scattering wave packet (detected by a monitor function), a nodeless covering wave function is added to it yielding a total function that is also nodeless. Both local and global choices for the covering function are described. The total and covering functions are then propagated separately in the hydrodynamic picture. At a later time, the actual wave function is recovered from the two propagated functions. The results obtained for Eckart barrier scattering in one dimension are in excellent agreement with exact results, even for very long propagation times t=1.2 ps. The capability of the CFM is also demonstrated for multidimensional propagation of a vibrationally excited wave packet.
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 quantum hydrodynamic equations associated with the de Broglie-Bohm formulation of quantum mechanics are solved using a new methodology which gives an accurate, unitary, and stable propagation of a time dependent quantum wave packet [B. K. Kendrick, J. Chem. Phys. 119, 5805 (2003)]. The methodology is applied to an N-dimensional model chemical reaction with an activation barrier. A parallel version of the methodology is presented which is designed to run on massively parallel supercomputers. The computational scaling properties of the parallel code are investigated both as a function of the number of processors and the dimension N. A decoupling scheme is introduced which decouples the multidimensional quantum hydrodynamic equations into a set of uncoupled one-dimensional problems. The decoupling scheme dramatically reduces the computation time and is highly parallelizable. Furthermore, the computation time is shown to scale linearly with respect to the dimension N=2,...,100.
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 scattered wave formalism developed for a quantum subsystem interacting with reservoirs through open boundaries is applied to one- or two-dimensional barrier scattering and quantum transistors. The total wave function is divided into incident and scattered components. Markovian outgoing wave boundary conditions are imposed on the scattered or total wave function by either the ratio or polynomial methods. For barrier scattering problems, accurate time-dependent transmission probabilities are obtained through the integration of the modified time-dependent Schrödinger equations for the scattered wave function. For quantum transistors, the time-dependent transport is studied for a quantum wave packet propagating through the conduction channel of a field effect transistor. This study shows that the scattered wave formalism significantly reduces computational effort relative to other open boundary methods and demonstrates wide applications to quantum dynamical processes.
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.
A new method is proposed for computing the time evolution of quantum mechanical wave packets. Equations of motion for the real-valued functions C and S in the complex action S=C(r,t)+iS(r,t)/ℏ, with ψ(r,t)=exp(S), involve gradients and curvatures of C and S. In previous implementations of the hydrodynamic formulation, various time-consuming fitting techniques of limited accuracy were used to evaluate these derivatives around each fluid element in an evolving ensemble. In this study, equations of motion are developed for the spatial derivatives themselves and a small set of these are integrated along quantum trajectories concurrently with the equations for C and S. Significantly, quantum effects can be included at various orders of approximation, no spatial fitting is involved, there are no basis set expansions, and single quantum trajectories (rather than correlated ensembles) may be propagated, one at a time. Excellent results are obtained when the derivative propagation method is applied to anharmonic potentials involving barrier transmission.