Reconciling semiclassical and Bohmian mechanics. II. Scattering states for discontinuous potentials

In a previous paper [B. Poirier, J. Chem. Phys. 121, 4501 (2004)] a unique bipolar decomposition, psi = psi1 + psi2, was presented for stationary bound states Psi of the one-dimensional Schrodinger equation, such that the components psi1 and psi2 approach their semiclassical WKB analogs in the large action limit. Moreover, by applying the Madelung-Bohm ansatz to the components rather than to Psi itself, the resultant bipolar Bohmian mechanical formulation satisfies the correspondence principle. As a result, the bipolar quantum trajectories are classical-like and well behaved, even when psi has many nodes or is wildly oscillatory. In this paper, the previous decomposition scheme is modified in order to achieve the same desirable properties for stationary scattering states. Discontinuous potential systems are considered (hard wall, step potential, and square barrier/well), for which the bipolar quantum potential is found to be zero everywhere, except at the discontinuities. This approach leads to an exact numerical method for computing stationary scattering states of any desired boundary conditions, and reflection and transmission probabilities. The continuous potential case will be considered in a companion paper [C. Trahan and B. Poirier, J. Chem. Phys. 124, 034116 (2006), following paper].     

Reconciling semiclassical and Bohmian mechanics. V. Wavepacket dynamics

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.     

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.     

Rovibrational spectroscopy calculations of neon dimer using a phase space truncated Weyl-Heisenberg wavelet basis

In a series of earlier articles [B. Poirier J. Theor. Comput. Chem. 2, 65 (2003); B. Poirier and A. Salam J. Chem. Phys. 121, 1690 (2004); B. Poirier and A. Salam J. Chem. Phys. 121, 1740 (2004)], a new method was introduced for performing exact quantum dynamics calculations in a manner that formally defeats exponential scaling with system dimensionality. The method combines an optimally localized, orthogonal Weyl-Heisenberg wavelet basis set with a simple phase space truncation scheme, and has already been applied to model systems up to 17 degrees of freedom (DOF's). In this paper, the approach is applied for the first time to a real molecular system (neon dimer), necessitating the development of an efficient numerical scheme for representing arbitrary potential energy functions in the wavelet representation. All bound rovibrational energy levels of neon dimer are computed, using both one DOF radial coordinate calculations and a three DOF Cartesian coordinate calculation. Even at such low dimensionalities, the approach is found to be competitive with another state-of-the-art method applied to the same system [J. Montgomery and B. Poirier J. Chem. Phys. 119, 6609 (2003)].     

Accurate quantum calculation of the bound and resonant rovibrational states of Li-(H2)

In a recent paper [B. Poirier, Chem. Phys. 308, 305 (2005)] a full-dimensional quantum method for computing the rovibrational dynamics of triatomic systems was presented, incorporating three key features: (1) exact analytical treatment of Coriolis coupling, (2) three-body "effective potential," and (3) a single bend angle basis for all rotational states. In this paper, these ideas are applied to the Li-(H2) electrostatic complex, to compute all of the rovibrational bound state energies, and a number of resonance energies and widths, to very high accuracy (thousandths of a wave number). This application is very challenging, owing to the long-range nature of the interaction and to narrow level spacings near dissociation. Nevertheless, by combining the present method with a G4 symmetry-adapted phase-space-optimized representation, only modest basis sizes are required for which the matrices are amenable to direct diagonalization. Several new bound levels are reported, as compared with a previous calculation [D. T. Chang, G. Surratt, G. Ristroff, and G. I. Gellene, J. Chem. Phys. 116, 9188 (2002)]. The resonances exhibit a clear-cut separation into shape and Feshbach varieties, with the latter characterized by extremely long lifetimes (microseconds or longer).     

Simulation of vibrational dephasing of I(2) in solid Kr using the semiclassical Liouville method

In this paper, we present simulations of the decay of quantum coherence between vibrational states of I(2) in its ground (X) electronic state embedded in a cryogenic Kr matrix. We employ a numerical method based on the semiclassical limit of the quantum Liouville equation, which allows the simulation of the evolution and decay of quantum vibrational coherence using classical trajectories and ensemble averaging. The vibrational level-dependent interaction of the I(2)(X) oscillator with the rare-gas environment is modeled using a recently developed method for constructing state-dependent many-body potentials for quantum vibrations in a many-body classical environment [J. M. Riga, E. Fredj, and C. C. Martens, J. Chem. Phys. 122, 174107 (2005)]. The vibrational dephasing rates gamma(0n) for coherences prepared between the ground vibrational state mid R:0 and excited vibrational state mid R:n are calculated as a function of n and lattice temperature T. Excellent agreement with recent experiments performed by Karavitis et al. [Phys. Chem. Chem. Phys. 7, 791 (2005)] is obtained.     

Unimolecular rovibrational bound and resonance states for large angular momentum: J=20 calculations for HO2

We explore the calculation of unimolecular bound states and resonances for deep-well species at large angular momentum using a Chebychev filter diagonalization scheme incorporating doubling of the autocorrelation function as presented recently by Neumaier and Mandelshtam [Phys. Rev. Lett. 86, 5031 (2001)]. The method has been employed to compute the challenging J=20 bound and resonance states for the HO2 system. The methodology has firstly been tested for J=2 in comparison with previous calculations, and then extended to J=20 using a parallel computing strategy. The quantum J-specific unimolecular dissociation rates for HO2-->H+O2 in the energy range from 2.114 to 2.596 eV have been reported for the first time, and comparisons with the results of Troe and co-workers [J. Chem. Phys. 113, 11019 (2000) Phys. Chem. Chem. Phys. 2, 631 (2000)] from statistical adiabatic channel method/classical trajectory calculations have been made. For most of the energies, the reported statistical adiabatic channel method/classical trajectory rate constants agree well with the average of the fluctuating quantum-mechanical rates. Near the dissociation threshold, quantum rates fluctuate more severely, but their average is still in agreement with the statistical adiabatic channel method/classical trajectory results.     

Nonadiabatic quantum dynamics in O(3P)+H2→OH+H: a revisited study

To investigate the extent of nonadiabatic effects in the title reaction, quasi-classical trajectory and nonadiabatic quantum scattering as well as the nonadiabatic quantum-classical trajectory calculations were performed on the accurate ab initio benchmark potential energy surfaces of the lowest (3)A' and (3)A" electronic states [Rogers et al., J Phys Chem A 2000, 104, 2308], together with the spin-orbit coupling matrix [Maiti and Schatz, J Chem Phys 2003, 119, 12360] and the lowest singlet (1) A' potential energy surface [Dobby and Knowles, Faraday Discuss 1998, 110, 247]. Comparison of the calculated total cross sections from both adiabatic and nonadiabatic calculations has demonstrated that for adiabatic channels including (3)A'→(3)A' and (3)A"→(3)A", difference does exist between the two kinds of adiabatic and nonadiabatic calculations, showing nonadiabatic effects to some extent. Such nonadiabatic effects tend to become more conspicuous at high collision energies and are found to be more pronounced with trajectories/quantum wave packet initiated on (3)A' than on (3)A". Furthermore, the present study also showed that nonadiabatic effects can bring the component of forward-scattering in the product angular distributions.     

Quantum-classical dynamics of scattering processes in adiabatic and diabatic representations

We demonstrate the workability of a TDDVR based [J. Chem. Phys. 118, 5302 (2003)], novel quantum-classical approach, for simulating scattering processes on a quasi-Jahn-Teller model [J. Chem. Phys. 105, 9141 (1996)] surface. The formulation introduces a set of DVR grid points defined by the Hermite part of the basis set in each dimension and allows the movement of grid points around the central trajectory. With enough trajectories (grid points), the method converges to the exact quantum formulation whereas with only one grid point, we recover the conventional molecular dynamics approach. The time-dependent Schrodinger equation and classical equations of motion are solved self-consistently and electronic transitions are allowed anywhere in the configuration space among any number of coupled states. Quantum-classical calculations are performed on diabatic surfaces (two and three) to reveal the effects of symmetry on inelastic and reactive state-to-state transition probabilities, along with calculations on an adiabatic surface with ordinary Born-Oppenheimer approximation. Excellent agreement between TDDVR and DVR results is obtained in both the representations.     

Interpretation of Adiabatic and Diabatic Populations from Trajectories of Branching Corrected Surface Hopping

The branching corrected surface hopping (BCSH) has been demonstrated as a robust approach to improve the performance of the traditional fewest switches surface hopping (FSSH) for nonadiabatic dynamics simulations of standard scattering problems [J. Chem. Phys. 150 , 164101 (2019)]. Here, we study how reliable populations of both adiabatic and diabatic states can be interpreted from BCSH trajectories. Using exact quantum solutions and FSSH results as references, we investigate a series of one-dimensional two-level scattering models and illustrate that excellent time-dependent populations can be obtained by BCSH. Especially, we show that different trajectory analysis strategies produce noticeable differences in different representations. Namely, the method based on active states performs better to get populations of adiabatic states, while the method based on wavefunctions produces more reliable results for populations of diabatic states.     

A quantum wave-packet study of intersystem crossing effects in the O(3P2,1,0,1D2)+H2 reaction

We present for the first time an exact quantum study of spin-orbit-induced intersystem crossing effects in the title reaction. The time-dependent wave-packet method, combined with an extended split operator scheme, is used to calculate the fine-structure resolved cross section. The calculation involves four electronic potential-energy surfaces of the 1A' state [J. Dobbyn and P. J. Knowles, Faraday Discuss. 110, 247 (1998)], the 3A' and the two degenerate 3A" states [S. Rogers, D. Wang, A. Kuppermann, and S. Wald, J. Phys. Chem. A 104, 2308 (2000)], and the spin-orbit couplings between them [B. Maiti, and G. C. Schatz, J. Chem. Phys. 119, 12360 (2003)]. Our quantum dynamics calculations clearly demonstrate that the spin-orbit coupling between the triplet states of different symmetries has the greatest contribution to the intersystem crossing, whereas the singlet-triplet coupling is not an important effect. A branch ratio of the spin state Pi32 to Pi12 of the product OH was calculated to be approximately 2.75, with collision energy higher than 0.6 eV, when the wave packet was initially on the triplet surfaces. The quantum calculation agrees quantitatively with the previous quasiclassical trajectory surface hopping study.     

Multidimensional quantum trajectories: applications of the derivative propagation method

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.     

Pair-correlated speed distributions for the OH+CH4/CD4 reactions: Further remarks on their classical trajectory calculations in a quantum spirit

Ten years ago, Liu and co-workers measured pair-correlated speed distributions for OH+CH₄/CD₄ reactions by means of velocity map imaging (VMI) techniques at a collision energy of ∼10 kcal/mol [B. Zhang, W. Shiu, J. J. Lin and K. Liu, J. Chem. Phys 122, 131102 (2005); B. Zhang, W. Shiu and K. Liu, J. Phys. Chem. A 109, 8989 (2005)]. Recently, two of us could semi-quantitatively reproduce these measurements by performing full-dimensional quasi-classical trajectory calculations in a quantum spirit on an ab-initio potential energy surface of their own [J. Espinosa-Garcia and J. C. Corchado, Theor. Chem. Acc. (2015) 134: 6; J. Phys. Chem. B 120, 1446 (2016)]. The goal of the present work is to show that these results can be significantly improved by adding a few more constraints in order to better comply with the restrictions imposed by VMI. Overall, the level of agreement between theory and experiment is remarkable owing to the large dimensionality of the reactions under scrutiny. This is an encouraging result considering the computational challenges of quantum scattering calculations for such large processes.     

Quantum dynamics calculations using symmetrized, orthogonal Weyl-Heisenberg wavelets with a phase space truncation scheme. II. Construction and optimization

In this paper, we extend and elaborate upon a wavelet method first presented in a previous publication [B. Poirier, J. Theo. Comput. Chem. 2, 65 (2003)]. In particular, we focus on construction and optimization of the wavelet functions, from theoretical and numerical viewpoints, and also examine their localization properties. The wavelets used are modified Wilson-Daubechies wavelets, which in conjunction with a simple phase space truncation scheme, enable one to solve the multidimensional Schrodinger equation. This approach is ideally suited to rovibrational spectroscopy applications, but can be used in any context where differential equations are involved.     

Quasi-classical trajectories study of Ne79Br2(B) vibrational predissociation

A full-dimensional quasi-classical trajectories study on the vibrational predissociation (VP) of the Ne79Br2(B) complex is presented. Following the most recent experiments, the Br2(B) vibrational levels v'=16-29 were explored. The total angular momentum, J, was taken to be zero, and a semiclassical Franck-Condon model to compute initial conditions from quantum distributions was employed. Predissociation lifetimes were extracted from Ne79Br2 population decay by using two different exponential laws. Predicted lifetimes are in excellent agreement with the last experimental results [J. A. Cabrera, C. R. Bieler, B. C. Olbricht, W. E. van der Veer and K. C. Janda, J. Chem. Phys., 2005, 123, 054311]. The Br2 fragment ro-vibrational distributions resulting from the VP of the molecule were obtained from the statistics of classical magnitudes using the standard binning procedure. Computed rotational distributions (for the Deltav'=-1, -2 channels) are also in very good agreement with the experimental results [M. Nejad-Sattari and T. A. Stephenson, J. Chem. Phys., 1997, 106 5454]. The influence of two quantum effects-the closing of the Deltav'=-1 dissociation channel and the intramolecular vibrational relaxation (IVR) mechanism-on the agreement with experimental rotational distributions, is discussed. Due to the classical character of our calculations and the binning procedure we used, the agreement of computed vibrational distributions with experimental and quantum theoretical is qualitative. For instance, for v'=28-for which the Deltav'=-1 channel is experimentally found to be closed-the Deltav'=-2 channel becomes statistically more significant. A discussion on the viability of similar quasi-classical methods to model the VP dynamics of analogous clusters is presented.     

A coherent state approach to semiclassical nonadiabatic dynamics

A semiclassical (SC) approximation to the quantum mechanical propagator for nonadiabatic systems is derived. Our derivation starts with an exact path integral expression that uses canonical coherent states for the nuclear degrees of freedom and spin coherent states for the electronic degrees of freedom. A stationary path approximation (SPA) is then applied to the path integral to obtain the SC approximation. The SPA results in complex classical trajectories of both nuclear and electronic degrees of freedom and a double ended boundary condition. The root search problem is solved using the previously proposed "real trajectory local search" algorithm. The SC approximation is tested on three simple one dimensional two-state systems proposed by Tully [J. Chem. Phys. 93, 1061 (1990)], and the SC results are compared to Ehrenfest and surface hopping predictions. Excellent agreement with quantum results is reached when the SC trajectory is far away from caustics. We discuss the origin of caustics in this SC formalism and the strengths and weaknesses of this approach.     

A time-dependent wave-packet quantum scattering study of the reaction H2+(v = 0-2,4,6;j = 1) + He--> HeH+ + H

The quantum scattering dynamics calculation was carried out for the titled reaction in the collision energy range of 0.0-2.4 eV with reactant H(2) (+) in the rotational state j = 1 and vibrational states v = 0-2, 4, and 6. The present time-dependent wave-packet calculation takes into account the Coriolis coupling (CC) and uses the accurate ab initio potential-energy surface of Palmieri et al. [Mol. Phys. 98, 1835 (2000)]. The importance of including the CC quantum scattering calculation has been revealed by the comparison between the CC calculation and the previous coupled state (CS) calculation. The CC total cross sections for the v = 2, 4, and 6 states show collision energy-dependent behaviors different from those based on the CS calculation. Furthermore, the collision energy dependence of the total cross sections obtained in the present CC calculation only exhibits minor oscillations, indicating that the chance is slim for reactive resonances in total cross sections to survive through the partial-wave averaging. The magnitude and profile of the CC total cross sections for v = 0-2 in the collision energy range of 0.0-2.5 eV are found to be consistent with experimental cross sections obtained recently by Tang et al. [J. Chem. Phys. 122, 164301 (2005)] after taking into account the experimental uncertainties.     

The integral encounter theory of multistage reactions containing association-dissociation reaction stages. III. Taking account of quantum states of reactants

The formalism developed in Part I [K. L. Ivanov, N. N. Lukzen, A. A. Kipriyanov, and A. B. Doktorov, Phys. Chem. Chem. Phys. 6, 1706 (2004)] of the present contribution is extended to treat the reacting particles with internal quantum states. Initial spatial correlations of reactants are considered in the framework of this formalism as well.