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

CRC Press eBooks

Direct Numerical Solution of the Quantum Hydrodynamic Equations of Motion

Published Apr 19, 2016• 0 citations
Publisher

Abstract

Systems .........................................................................................337 21.5 The Iterative Finite Difference Method ......................................................339 21.6 Concluding Remarks and Future Research.................................................342 Acknowledgments..................................................................................................343 Bibliography ..........................................................................................................343 The quantum hydrodynamic formulation of quantum mechanics is intuitively appealing [1-7]. The quantum potential (Q) and its associated force fq = −∇Q appear on an equal footing with the classical potential and force in the equations of motion. Thus, this approach provides a unified computational method for including both classical and quantum mechanical effects in a molecular dynamics (trajectory based) calculation. However, despite their deceptively simple form, the quantum hydrodynamic equations of motion are very challenging to solve numerically. The quantum potential is non-local (i.e., the trajectories are coupled) and it is a non-linear function of the density (i.e., it depends upon the solution itself). Thus, unlike a typical classical potential, the quantum potential is manifestly time dependent and evolves in “lock-step” fashion with the dynamics. Furthermore, the quantum potential can

DOI
10.1201/b10287-26

Related papers

Quantum Fluid Dynamics within the Framework of Density Functional Theory

Framework for QFD in 3-D Space .............................................................188 12.5 Quantum Fluid Dynamical Equations within a Density Functional Framework..................................................................................................191 12.6 Miscellaneous Aspects of the Quantum Fluid Dynamical Approach .........193 12.7 Concluding Remarks...................................................................................194 Acknowledgments..................................................................................................194 Bibliography ..........................................................................................................194 The concepts of quantum mechanics (QM) are so radically different from those of classical mechanics that soon after Schrödinger proposed the concept of quantum mechanical wavefunction, there were attempts to provide a classical interpretation of QM. In fact, the success of Madelung [1] in transforming the Schrödinger equation into a pair of hydrodynamical equations consisting of the continuity equation and an Euler type equation of hydrodynamics, giving birth to so-called quantum fluid dynamics (QFD) [2-5], in 1926, itself is testimony to the strong urge to have a classical interpretation of the strange world of QM. The concept was resurrected with full vigor more than 25 years later through the pioneering works of Bohm [6] and Takabayashi [7]. In the meantime, Wigner [8] introduced, in 1932, the concept

CRC Press eBooks•Apr 19, 2016•1 citation

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

An iterative finite difference method for solving the quantum hydrodynamic equations of motion

Brian K. Kendrick

Journal of Molecular Structure THEOCHEM•Dec 18, 2009•14 citations

Trajectory integration of the quantum hydrodynamic equations of motion

Corey Jason Trahan

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.

•Jan 1, 2003•0 citations

The Derivation and Mathematical Models of Quantum Hydrodynamic Equations

WORLD SCIENTIFIC eBooks•Jun 29, 2023•0 citations

Theoretical methods for the study of quantum dynamics

Jonathan E. Chen

Several theoretical methods for the computation of quantum dynamical quantities are formulated, implemented, and applied, in the overall context of overcoming the exponential scaling property of quantum mechanics. The convergence of the approximate forward-backward semiclassical dynamics (FBSD) method has been accelerated by working with only the imaginary part of the real-time position and momentum time correlation functions. Simple manipulations allow the computationally favorable limit of a single factorization (“bead”) of the Boltzmann operator to yield greatly improved results. Numerical examples for a one-dimensional anharmonic oscillator and Lennard-Jones fluids are presented. Bohm's trajectory formulation of quantum mechanics is an exact mapping of the time-dependent Schrödinger equation onto classical fluid dynamics, plus one extra term responsible for all non-classical effects, the dreaded quantum potential. A computational procedure has been introduced to stably solve Bohm's hydrodynamic equations of motion in the forward-backward quantum dynamics (FBQD) formulation of time correlation functions, and applied to one-dimensional bound systems at zero and low temperature. A key step in the formulation of this procedure was the use of Hamilton's law of varying action, a generalization of Hamilton's principle of stationary action applicable to initial-value problems. This has resulted in the first reported synthetic computation of non-stationary quantum trajectories for a fully anharmonic oscillator. Success in one dimension has laid the groundwork for future extension to a rigorous quantum-classical scheme, which is briefly discussed. Other unpublished methods for obtaining time correlation functions and propagators in quantum dynamics are also discussed.

•May 19, 2010•0 citations
Systems .........................................................................................337 21.5 The Iterative Finite Difference Method ......................................................339 21.6 Concluding Remarks and Future Research.................................................342 Acknowledgments..................................................................................................343 Bibliography ..........................................................................................................343 The quantum hydrodynamic formulation of quantum mechanics is intuitively appealing [1-7]. The quantum potential (Q) and its associated force fq = −∇Q appear on an equal footing with the classical potential and force in the equations of motion. Thus, this approach provides a unified computational method for including both classical and quantum mechanical effects in a molecular dynamics (trajectory based) calculation. However, despite their deceptively simple form, the quantum hydrodynamic equations of motion are very challenging to solve numerically. The quantum potential is non-local (i.e., the trajectories are coupled) and it is a non-linear function of the density (i.e., it depends upon the solution itself). Thus, unlike a typical classical potential, the quantum potential is manifestly time dependent and evolves in “lock-step” fashion with the dynamics. Furthermore, the quantum potential can
Framework for QFD in 3-D Space .............................................................188 12.5 Quantum Fluid Dynamical Equations within a Density Functional Framework..................................................................................................191 12.6 Miscellaneous Aspects of the Quantum Fluid Dynamical Approach .........193 12.7 Concluding Remarks...................................................................................194 Acknowledgments..................................................................................................194 Bibliography ..........................................................................................................194 The concepts of quantum mechanics (QM) are so radically different from those of classical mechanics that soon after Schrödinger proposed the concept of quantum mechanical wavefunction, there were attempts to provide a classical interpretation of QM. In fact, the success of Madelung [1] in transforming the Schrödinger equation into a pair of hydrodynamical equations consisting of the continuity equation and an Euler type equation of hydrodynamics, giving birth to so-called quantum fluid dynamics (QFD) [2-5], in 1926, itself is testimony to the strong urge to have a classical interpretation of the strange world of QM. The concept was resurrected with full vigor more than 25 years later through the pioneering works of Bohm [6] and Takabayashi [7]. In the meantime, Wigner [8] introduced, in 1932, the concept
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.
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.
Several theoretical methods for the computation of quantum dynamical quantities are formulated, implemented, and applied, in the overall context of overcoming the exponential scaling property of quantum mechanics. The convergence of the approximate forward-backward semiclassical dynamics (FBSD) method has been accelerated by working with only the imaginary part of the real-time position and momentum time correlation functions. Simple manipulations allow the computationally favorable limit of a single factorization (“bead”) of the Boltzmann operator to yield greatly improved results. Numerical examples for a one-dimensional anharmonic oscillator and Lennard-Jones fluids are presented. Bohm's trajectory formulation of quantum mechanics is an exact mapping of the time-dependent Schrödinger equation onto classical fluid dynamics, plus one extra term responsible for all non-classical effects, the dreaded quantum potential. A computational procedure has been introduced to stably solve Bohm's hydrodynamic equations of motion in the forward-backward quantum dynamics (FBQD) formulation of time correlation functions, and applied to one-dimensional bound systems at zero and low temperature. A key step in the formulation of this procedure was the use of Hamilton's law of varying action, a generalization of Hamilton's principle of stationary action applicable to initial-value problems. This has resulted in the first reported synthetic computation of non-stationary quantum trajectories for a fully anharmonic oscillator. Success in one dimension has laid the groundwork for future extension to a rigorous quantum-classical scheme, which is briefly discussed. Other unpublished methods for obtaining time correlation functions and propagators in quantum dynamics are also discussed.