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.
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.
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.
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.
The quantum Hamilton-Jacobi equation for the action function is approximately solved by propagating individual Lagrangian quantum trajectories in complex-valued phase space. Equations of motion for these trajectories are derived through use of the derivative propagation method (DPM), which leads to a hierarchy of coupled differential equations for the action function and its spatial derivatives along each trajectory. In this study, complex-valued classical trajectories (second order DPM), along which is transported quantum phase information, are used to study low energy barrier transmission for a model two-dimensional system involving either an Eckart or Gaussian barrier along the reaction coordinate coupled to a harmonic oscillator. The arrival time for trajectories to reach the transmitted (product) region is studied. Trajectories launched from an "equal arrival time surface," defined as an isochrone, all reach the real-valued subspace in the transmitted region at the same time. The Rutherford-type diffraction of trajectories around poles in the complex extended Eckart potential energy surface is described. For thin barriers, these poles are close to the real axis and present problems for computing the transmitted density. In contrast, for the Gaussian barrier or the thick Eckart barrier where the poles are further from the real axis, smooth transmitted densities are obtained. Results obtained using higher-order quantum trajectories (third order DPM) are described for both thick and thin barriers, and some issues that arise for thin barriers are examined.
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) 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.
In this study, an analysis of the one-dimensional Eckart and Gaussian barrier scattering problems is undertaken using approximate quantum trajectories. Individual quantum trajectories are computed using the derivative propagation method (DPM). Both real-valued and complex-valued DPM quantum trajectories are employed. Of interest are the deep tunneling and the higher energy barrier scattering problems in cases in which the scattering barrier is "thick" by comparison to the width of the initial wave packet. For higher energy scattering problems, it is found that real-valued DPM trajectories very accurately reproduce the transmitted probability densities at low orders when compared to large fixed-grid calculations. However, higher orders must be introduced to obtain good probabilities for deep tunneling problems. Complex-valued DPM is found to accurately reproduce transmitted probability densities at low order for both the deep tunneling and the higher energy scattering problems. Of particular note, complex-classical trajectories are found to very nearly give the exact result for the deep barrier tunneling scattering problem, and the complex DPM converges well at high orders for these thick barrier scattering problems. A variety of analyses are performed to elucidate the dynamics of complex-valued DPM trajectories. The complex-extended barrier potentials are examined in detail, including an analysis of the complex force. Of particular interest are initial conditions for complex-valued DPM trajectories known as isochrones. All trajectories launched from an isochrone arrive on the real axis on the transmitted side of the barrier at the same time. The computation and properties of isochrones as well as the behavior of the initial wave packet in the complex plane are also examined.