Quantum trajectory methods (QTMs) hold great promise as a potential means of obtaining dynamical insight and computational scaling similar to classical trajectory simulations but in an exact quantum dynamical context. To date, the development of QTMs has been stymied by the "node problem"--highly nonclassical and numerically unstable trajectories that arise when the wavepacket density |psi|2 exhibits substantial interference oscillations. In a recent paper, however [B. Poirier, J. Chem. Phys. 128, 164115 (2008)], a "bipolar decomposition," psi=psi+(+)(psi)psi(-), was introduced for one-dimensional (1D) wavepacket dynamics calculations such that the component densities |psi(+)(-)|2 are slowly varying and otherwise interference-free, even when |psi|2 itself is highly oscillatory. The bipolar approach is thus ideally suited to a QTM implementation, as is demonstrated explicitly in this paper. Two model 1D benchmark systems exhibiting substantial interference are considered--one with more "quantum" system parameters and the other more classical-like. For the latter, more challenging application, synthetic QTM results are obtained and found to be extremely accurate, as compared to a corresponding fixed-grid calculation. Ramifications of the bipolar QTM approach for the classical limit and also for multidimensional applications, are discussed.
Quantum trajectory methods (QTMs) hold great promise as a potential means of obtaining dynamical insight and computational scaling similar to classical trajectory simulations but in an exact quantum dynamical context. To date, the development of QTMs has been stymied by the "node problem"--highly nonclassical and numerically unstable trajectories that arise when the wavepacket density |psi|2 exhibits substantial interference oscillations. In a recent paper, however [B. Poirier, J. Chem. Phys. 128, 164115 (2008)], a "bipolar decomposition," psi=psi+(+)(psi)psi(-), was introduced for one-dimensional (1D) wavepacket dynamics calculations such that the component densities |psi(+)(-)|2 are slowly varying and otherwise interference-free, even when |psi|2 itself is highly oscillatory. The bipolar approach is thus ideally suited to a QTM implementation, as is demonstrated explicitly in this paper. Two model 1D benchmark systems exhibiting substantial interference are considered--one with more "quantum" system parameters and the other more classical-like. For the latter, more challenging application, synthetic QTM results are obtained and found to be extremely accurate, as compared to a corresponding fixed-grid calculation. Ramifications of the bipolar QTM approach for the classical limit and also for multidimensional applications, are discussed.
In a previous paper [Park K, Poirier B, Parlant G, J Chem Phys129:194112, 2008], a synthetic quantum trajectory method (QTM) was successfully implemented for wave-packet dynamics in a one-dimensional (1D) symmetric Eckart barrier system, utilizing a "double-wavepacket" version of the bipolar decomposition, ψ = ψ++ ψ-= (ψ1++ ψ2+) + (ψ1-+ ψ2-), to avoid a technical difficulty involving negligible initial ψ-density. In this paper, we develop a new synthetic algorithm which overcomes this difficulty directly, utilizing the original "single-wavepacket" version of the bipolar decomposition, ψ =ψ++ ψ-, and also show that the initial propagation of ψ-is mainly governed by probability transfer from ψ+, rather than by the given initial conditions for ψ-. The new algorithm makes it possible to apply the synthetic bipolar QTM to asymptotically asymmetric as well as symmetric potential systems. Successful application results for both symmetric and asymmetric Eckart barrier systems in 1D are presented.
In previous articles [B. Poirier J. Chem. Phys. 121, 4501 (2004); C. Trahan and B. Poirier, ibid. 124, 034115 (2006); 124, 034116 (2006); B. Poirier and G. Parlant, J. Phys. Chem. A 111, 10400 (2007)] a bipolar counterpropagating wave decomposition, psi = psi(+) + psi(-), was presented for stationary states psi of the one-dimensional Schrodinger equation, such that the components psi(+/-) approach their semiclassical Wentzel-Kramers-Brillouin analogs in the large action limit. The corresponding bipolar quantum trajectories are classical-like and well behaved, even when psi has many nodes, or is wildly oscillatory. In this paper, the method is generalized for time-dependent wavepacket dynamics applications and applied to several benchmark problems, including multisurface systems with nonadiabatic coupling.
In previous articles [J. Chem. Phys. 121, 4501 (2004); J. Chem. Phys. 124, 034115 (2006); J. Chem. Phys. 124, 034116 (2006); J. Phys. Chem. A 111, 10400 (2007); J. Chem. Phys. 128, 164115 (2008)] an exact quantum, bipolar wave decomposition, psi=psi(+)+psi(-), was presented for one-dimensional stationary state and time-dependent wavepacket dynamics calculations, such that the components psi(+/-) approach their semiclassical WKB analogs in the large action limit. The corresponding bipolar quantum trajectories are classical-like and well behaved, even when psi has many nodes or is wildly oscillatory. In this paper, both the stationary state and wavepacket dynamics theories are generalized for multidimensional systems and applied to several benchmark problems, including collinear H+H(2).
In this work we present a numerical analysis of the so-called ‘interacting quantum trajectory method’ (IQTM), which is based on a formal theory developed by Poirier [Bohmian mechanics without pilot waves, Chem. Phys. 370 (2010), pp. 4–14.], and later extended by Schiff and Poirier [Communication: Quantum mechanics without wavefunctions, J. Chem. Phys. 136 (2012), p. 031102.], for representing the quantum mechanical state of a particle in terms of an ensemble of Bohmian quantum trajectories. This has led to a number of numerical algorithms over the years, mostly for time-independent applications. A numerically robust, exact time-dependent wavepacket propagation scheme would have to await the refinements of Dupuy et al. [L. Dupuy, F. Talotta, F. Agostini, D. Lauvergnat, B. Poirier, and Y. Scribano, Adiabatic and nonadiabatic dynamics with interacting quantum trajectories, J. Chem. Theory Comput. 18 (2022), pp. 6447–6462.], which serves as the starting point for our continued development in the present study. Here we will discuss the numerical protocols comprising IQTM, and present detailed error analysis for several 1D model problems involving simple free particles and oscillators. These studies reveal some key advantages of IQTM; however, they also highlight problems for which the algorithms need further work.
In a previous paper [B. Poirier, J. Chem. Phys. 121, 4501 (2004)] a unique bipolar decomposition Ψ=Ψ1+Ψ2 was presented for stationary bound states Ψ of the one-dimensional Schrödinger equation, such that the components Ψ1 and Ψ2 approach their semiclassical WKB analogs in the large-action limit. The corresponding bipolar quantum trajectories, as defined in the usual Bohmian mechanical formulation, are classical-like and well behaved, even when Ψ has many nodes or is wildly oscillatory. A modification for discontinuous potential stationary scattering states was presented in a second, companion paper [C. Trahan and B. Poirier, J. Chem. Phys.124, 034115 (2006), previous paper], whose generalization for continuous potentials is given here. The result is an exact quantum scattering methodology using classical trajectories. For additional convenience in handling the tunneling case, a constant-velocity-trajectory version is also developed.
In previous articles (J. Chem. Phys. 2004, 121, 4501; 2006, 124, 034115; 2006, 124, 034116) a bipolar counter-propagating wave decomposition, Psi = Psi+ + Psi-, was presented for stationary states Psi of the one-dimensional Schrödinger equation, such that the components Psi+/- approach their semiclassical WKB analogs in the large action limit. The corresponding bipolar quantum trajectories are classical-like and well-behaved, even when Psi has many nodes or is wildly oscillatory. In this paper, the method is generalized for multisurface scattering applications and applied to several benchmark problems. A natural connection is established between intersurface transitions and (+ <--> -) transitions.