State-to-state quasi-classical trajectory study of the D + H2 collision for high temperature astrophysical applications.

We report state-to-state quasi-classical trajectory rate constants for the D + H2 reactive collision, using the accurate H3 global potential energy surface of Mielke et al. [J. Chem. Phys. 116, 4142 (2002)]. High relative collision energies (up to ≈56 000 K) and high rovibrational levels of H2 (up to ≈50 000 K), relevant to various astrophysical media, are considered. The HD product molecule is formed in highly excited rovibrational states, over a wide collision energy range. The collision-induced dissociation channel (often overlooked in fully quantum reaction dynamics calculations) is found to be significantly populated, even at collision energies as low as 1500 K.

Communication: Adiabatic quantum trajectory capture for cold and ultra-cold chemical reactions.

The Langevin capture model is often used to describe barrierless reactive collisions. At very low temperatures, quantum effects may alter this simple capture image and dramatically affect the reaction probability. In this paper, we use the trajectory-ensemble reformulation of quantum mechanics, as recently proposed by one of the authors (Poirier) to compute adiabatic-channel capture probabilities and cross-sections for the highly exothermic reaction Li + CaH(v = 0, j = 0) → LiH + Ca, at low and ultra-low temperatures. Each captured quantum trajectory takes full account of tunneling and quantum reflection along the radial collision coordinate. Our approach is found to be very fast and accurate, down to extremely low temperatures. Moreover, it provides an intuitive and practical procedure for determining the capture distance (i.e., where the capture probability is evaluated), which would otherwise be arbitrary.

Quantum trajectory calculations for bipolar wavepacket dynamics in one dimension.

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.

Reconciling semiclassical and Bohmian mechanics: IV. Multisurface dynamics.

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.

A hybrid hydrodynamic-Liouvillian approach to mixed quantum-classical dynamics: application to tunneling in a double well.

The hybrid quantum-classical approach of Burghardt and Parlant [Burghardt, I.; Parlant, G. J. Chem. Phys. 2004, 120, 3055], referred to here as the quantum-classical moment (QCM) approach, is demonstrated for the dynamics of a quantum double well coupled to a classical harmonic coordinate. The approach combines the quantum hydrodynamic and classical Liouvillian representations by the construction of a particular type of moments (that is, partial hydrodynamic moments) whose evolution is determined by a hierarchy of coupled equations. For pure states, which are at the center of the present study, this hierarchy terminates at the first order. In the Lagrangian picture, the deterministic trajectories result in dynamics which is Hamiltonian in the classical subspace, while the projection onto the quantum subspace evolves under a generalized hydrodynamic force. Importantly, this force also depends upon the classical (Q, P) variables. The present application demonstrates the tunneling dynamics in both the Eulerian and Lagrangian representations. The method is exact if the classical subspace is harmonic, as is the case for the systems studied here.

Combining fixed- and moving-grid methods to study direct dissociation processes involving nonadiabatic transitions.

We present a novel quantum-dynamics approach suitable for computing direct dissociation processes, including electronic transitions. This approach combines quantum trajectories in the Lagrangian reference frame with standard fixed-grid wave packets in order to overcome the limitations and difficulties of both techniques. As a model application, we consider the ultrafast photodissociation of H2 excited by a femtosecond extreme UV laser pulse.

On the dynamics of coupled Bohmian and phase-space variables: a new hybrid quantum-classical approach.

A new approach to the coupling of quantum and classical dynamics is developed, by combining a hydrodynamic, Bohmian description for the quantum subsystem with a Liouville-space description for the classical subsystem. To this end, partial hydrodynamic moments are introduced, the dynamics of which is determined by a hierarchy of equations derived from the quantum Liouville equation. We focus on pure states (wave functions) and introduce a trajectory representation in a hybrid hydrodynamic-Liouvillian phase space. The interleaved trajectory dynamics is guided by a new type of quantum force. For illustration, we consider a pair of bilinearly coupled harmonic oscillators, for which the method is exact.