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.
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.
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.
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.
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 Caldeira-Leggett master equation, expressed in Lindblad form, has been used in the numerical study of the effect of a thermal environment on the dynamics of the scattering of a wave packet from a repulsive Eckart barrier. The dynamics are studied in terms of phase space trajectories associated with the distribution function, W(q,p,t). The equations of motion for the trajectories include quantum terms that introduce nonlocality into the motion, which imply that an ensemble of correlated trajectories needs to be propagated. However, use of the derivative propagation method (DPM) allows each trajectory to be propagated individually. This is achieved by deriving equations of motion for the partial derivatives of W(q,p,t) that appear in the master equation. The effects of dissipation on the trajectories are studied and results are shown for the transmission probability. On short time scales, decoherence is demonstrated by a swelling of trajectories into momentum space. For a nondissipative system, a comparison is made of the DPM with the "exact" transmission probability calculated from a fixed grid calculation.
In recent years trajectory based methodologies have become increasingly popular for evaluating the time evolution of quantum systems. A revival of the de Broglie Bohm interpretation of quantum mechanics has spawned several such techniques for examining quantum dynamics from a hydrodynamic perspective. Using techniques similar to those found in computational fluid dynamics one can construct the wavefunction of a quantum system at any time from the trajectories of a discrete ensemble of hydrodynamic fluid elements (Bohm particles) which evolve according to non-classical equations of motion. Until very recently these schemes have been limited to closed systems. In this paper we present our methodology for including the effects of a thermal environment into the hydrodynamic formulation of quantum dynamics. We derive hydrodynamic equations of motion from the Caldeira-Leggett master equation for the reduced density matrix and give a brief overview of our computational scheme which incorporates an adaptive Lagrangian mesh. Our applications focus upon the dissipative dynamics of open unbounded quantum systems. Using both the Wigner phase space representation and the linear entropy we probe the breakdown of the Markov approximation of the bath dynamics at low temperatures. We suggest a criteria for rationalizing the validity of the Markov approximation in open unbound systems and discuss decoherence, energy relaxation, and quantum/classical correspondence in the context of the Bohmian paths.
In recent years trajectory based methodologies have become increasingly popular for evaluating the time evolution of quantum systems. A revival of the de Broglie--Bohm interpretation of quantum mechanics has spawned several such techniques for examining quantum dynamics from a hydrodynamic perspective. Using techniques similar to those found in computational fluid dynamics one can construct the wave function of a quantum system at any time from the trajectories of a discrete ensemble of hydrodynamic fluid elements (Bohm particles) which evolve according to nonclassical equations of motion. Until very recently these schemes have been limited to conservative systems. In this paper, we present our methodology for including the effects of a thermal environment into the hydrodynamic formulation of quantum dynamics. We derive hydrodynamic equations of motion from the Caldeira-Leggett master equation for the reduced density matrix and give a brief overview of our computational scheme that incorporates an adaptive Lagrangian mesh. Our applications focus upon the dissipative dynamics of open unbounded quantum systems. Using both the Wigner phase space representation and the linear entropy, we probe the breakdown of the Markov approximation of the bath dynamics at low temperatures. We suggest a criteria for rationalizing the validity of the Markov approximation in open unbound systems and discuss decoherence, energy relaxation, and quantum/classical correspondence in the context of the Bohmian paths.