A time-dependent quantum mechanical investigation of dynamical resonances in three-dimensional HeH2+ and HeHD+ systems

Time-dependent quantum mechanical wave packet dynamics

Starting from a model study of the collinear (H, H₂) exchange reaction in 1959, the time-dependent quantum mechanical wave packet (TDQMWP) method has come a long way in dealing with systems as large as Cl + CH₄. The fast Fourier transform method for evaluating the second order spatial derivative of the wave function and split-operator method or Chebyshev polynomial expansion for determining the time evolution of the wave function for the system have made the approach highly accurate from a practical point of view. The TDQMWP methodology has been able to predict state-to-state differential and integral reaction cross sections accurately, in agreement with available experimental results for three dimensional (H, H₂) collisions, and identify reactive scattering resonances too. It has become a practical computational tool in predicting the observables for many A + BC exchange reactions in three dimensions and a number of larger systems. It is equally amenable to determining the bound and quasi-bound states for a variety of molecular systems. Just as it is able to deal with dissociative processes (without involving basis set expansion), it is able to deal with multi-mode nonadiabatic dynamics in multiple electronic states with equal ease. We present an overview of the method and its strength and limitations, citing examples largely from our own research groups.

Dynamical resonances of the deuterated CH2+ complex in the electronic ground state: A quantum wavepacket study

We here investigate the effects of isotopic substituents on the vibrational energy levels of the CH₂⁺ complex, supported by the electronic ground (1 ²A') potential energy surface (PES) of the H + CH⁺ reaction. We calculate the transition state spectrum by Fourier transforming the time-autocorrelation function of the initial wavepacket (WP) chosen in the interaction region of the PES. Using the time-dependent WP approach, the dynamical resonances are identified as bound and quasibound in nature, and they are characterized in terms of the eigenfunctions and lifetimes. The present work on the isotopic variants [CHD⁺(CDH⁺) and CD₂⁺] is compared with our earlier work [P. Sundaram et al., Phys. Chem. Chem. Phys. 19, 20172 (2017)] on the parent CH₂⁺ species. The isotopic variants reveal a large number of peaks in the spectra and the eigenfunctions exhibit the systematic nodal progressions and periodic orbits, the same as in CH₂⁺. While the CD₂⁺ complex exactly mimics the resonance behaviors (local and hyperspherical modes) of the bound and quasibound CH₂⁺ complex, the CHD⁺(CDH⁺) complex reveals only the local mode behaviors at low energies and significantly less number of resonance structures at high energies. Lifetime analysis of the isotopic variants implies that the CD₂⁺ complex survives much longer than the CHD⁺(CDH⁺) complex and concludes the work by noting the following order in the decay profile of the deuterated CH₂⁺ resonances as CH₂⁺>CHD⁺(CDH⁺) >CD₂⁺.

Time-dependent wave packet quantum and quasi-classical trajectory study of He + H₂⁺, D₂⁺ → HeH⁺ + H, HeD⁺ + D reaction on an accurate FCI potential energy surface

The quantum scattering dynamics and quasi-classical trajectory (QCT) calculations have been carried out for the title reaction on an accurate potential energy surface (PES) computed using the full configuration interaction (FCI). On the basis of the PES, the integral cross-sections of He + H₂⁺ (v = 0-3, j = 1) → HeH⁺ + H reaction have been calculated, and the results are generally agreed with the experimental cross-sections obtained by Tang et al. [J. Chem. Phys. 2005, 122, 164301] after taking into account the experimental uncertainties, which proves the reliability of implementing dynamics calculations on the FCI PES. The reaction probability of He + D₂⁺ (v = 0-2, j = 0) → HeD⁺ + D reactions for total angular momentum J = 0 and the integral cross-section (ICS) have been calculated. The significant quantum effect has been explored by the comparison between the QCT reaction probabilities (or ICS) and the quantum mechanical (QM) reaction probabilities (or ICS), which may be attributed to the deep well in the PES of this light atoms system. Furthermore, the role of Coriolis coupling (CC) effects has also been found not important by the comparison between the CC calculation and the centrifugal sudden (CS) approximation calculation, except that the CC total cross-sections for the v = 1 and 2 states show the collision energy-dependent behaviors in the low-energy area, which are different from those based on the CS calculation.

Influence of Reagent Rotation on (H-, D2) and (D-, H2) Collisions: A Quantum Mechanical Study

Time-independent quantum mechanical (TIQM) approach (helicity basis truncated at k = 2) has been used for computing differential and integral cross sections for the exchange reaction H- + D2 (v = 0, j = 0-4) --> HD + D- and D- + H2 (v = 0, j = 0-3) --> HD + H- in three dimensions on an accurate ab initio potential energy surface. It is shown that the j-weighted differential reaction cross section values are in good agreement with the experimental results reported by Zimmer and Linder at four different relative translational energies (Etrans = 0.55, 0.93, 1.16 and 1.48 eV) for (H-, D2) and at one relative translational energy (Etrans = 0.6 eV) by Haufler et al. for both (H-, D2) and (D-, H2) collisions. The j-weighted integral reaction cross section values are in good agreement with the crossed beam measurements by Zimmer and Linder in the Etrans range 0.5-1.5 eV and close to the guided ion beam results by Haufler et al. for (H-, D2) in the range 0.8-1.2 eV. Time-dependent quantum mechanical (TDQM) results obtained using centrifugal sudden approximation are reported in the form of integral reaction cross section values as a function of Etrans in the range 0.3-3.0 eV for both reactions in three dimensions on the same potential energy surface. The TDQM reaction cross section values decline more sharply than the TIQM results with increase in the initial rotational quantum number (j) for the D2 molecules in their ground vibrational state (v = 0) for (H-, D2) collisions. The computed j-weighted reaction cross section values are in good agreement with the experimental results reported by Zimmer and Linder for (H-, D2) collisions and guided ion beam results by Haufler et al. for both (H-, D2) and (D-, H2) collisions for energies below the threshold for electron detachment channel.

The study of state-selected ion-molecule reactions using the vacuum ultraviolet pulsed field ionization-photoion technique

This paper presents the methodology to generate beams of ions in single quantum states for bimolecular ion-molecule reaction dynamics studies using pulsed field ionization (PFI) of atoms or molecules in high-n Rydberg states produced by vacuum ultraviolet (VUV) synchrotron or laser photoexcitation. Employing the pseudocontinuum high-resolution VUV synchrotron radiation at the Advanced Light Source as the photoionization source, PFI photoions (PFI-PIs) in selected rovibrational states have been generated for ion-molecule reaction studies using a fast-ion gate to pass the PFI-PIs at a fixed delay with respect to the detection of the PFI photoelectrons (PFI-PEs). The fast ion gate provided by a novel interleaved comb wire gate lens is the key for achieving the optimal signal-to-noise ratio in state-selected ion-molecule collision studies using the VUV synchrotron based PFI-PE secondary ion coincidence (PFI-PESICO) method. The most recent development of the VUV laser PFI-PI scheme for state-selected ion-molecule collision studies is also described. Absolute integral cross sections for state-selected H2+ ions ranging from v+ = 0 to 17 in collisions with Ar, Ne, and He at controlled translational energies have been obtained by employing the VUV synchrotron based PFI-PESICO scheme. The comparison between PFI-PESICO cross sections for the H2+(HD+)+Ne and H2+(HD+)+He proton-transfer reactions and theoretical cross sections based on quasiclassical trajectory (QCT) calculations and three-dimensional quantum scattering calculations performed on the most recently available ab initio potential energy surfaces is highlighted. In both reaction systems, quantum scattering resonances enhance the integral cross sections significantly above QCT predictions at low translational and vibrational energies. At higher energies, the agreement between experiment and quasiclassical theory is very good. The profile and magnitude of the kinetic energy dependence of the absolute integral cross sections for the H2+(v+ = 0-2,N+ = 1)+He proton-transfer reaction unambiguously show that the inclusion of Coriolis coupling is important in quantum dynamics scattering calculations of ion-molecule collisions.

A pulsed-field ionization photoelectron secondary ion coincidence study of the H2+ (X,upsilon+=0-15,N+=1)+He proton transfer reaction

The endothermic proton transfer reaction, H2+(upsilon+)+He-->HeH+ + H(DeltaE=0.806 eV), is investigated over a broad range of reactant vibrational levels using high-resolution vacuum ultraviolet to prepare reactant ions either through excitation of autoionization resonances, or using the pulsed-field ionization-photoelectron-secondary ion coincidence (PFI-PESICO) approach. In the former case, the translational energy dependence of the integral reaction cross sections are measured for upsilon+=0-3 with high signal-to-noise using the guided-ion beam technique. PFI-PESICO cross sections are reported for upsilon+=1-15 and upsilon+=0-12 at center-of-mass collision energies of 0.6 and 3.1 eV, respectively. All ion reactant states selected by the PFI-PESICO scheme are in the N+=1 rotational level. The experimental cross sections are complemented with quasiclassical trajectory (QCT) calculations performed on the ab initio potential energy surface provided by Palmieri et al. [Mol. Phys. 98, 1839 (2000)]. The QCT cross sections are significantly lower than the experimental results near threshold, consistent with important contributions due to resonances observed in quantum scattering studies. At total energies above 2 eV, the QCT calculations are in excellent agreement with the present results. PFI-PESICO time-of-flight (TOF) measurements are also reported for upsilon+=3 and 4 at a collision energy of 0.6 eV. The velocity inverted TOF spectra are consistent with the prevalence of a spectator-stripping mechanism.

Time-dependent quantum mechanical wave packet study of the He+H(2) (+)(v,j)-->HeH(+)+H reaction

A detailed three-dimensional time-dependent quantum dynamical study of the He+H(2) (+)(v=0-3,j=0)-->HeH(+)+H reaction is reported for different vibrational v states of H(2) (+) in its ground rotational (j=0) state over a range of translational E(trans) energies on an accurate ab initio potential energy surface published by Palmieri et al. Plots of reaction probability as a function of total energy E reveal a large number of oscillations indicating the presence of a number of reactive scattering resonances. When averaged over total angular momentum J, some of the oscillations survive, indicating that they may be amenable to experimental observation. A comparison of our present results with our earlier results on the McLaughlin-Thompson-Joseph-Sathyamurthy surface and the experimental results from different research groups reveal a good deal of agreement as well as some discrepancies between theory and experiment at the level of state-selected gas phase dynamics.

Elastic and charge transfer processes in H+ + CO collisions

Proton and hydrogen atom time-of-flight spectra in collision energy range of E(trans) = 9.5-30 eV show that the endoergic charge transfer process in the H+ + CO system is almost an order of magnitude less probable than the elastic scattering [G. Niedner-Schatteburg and J. P. Toennies, Adv. Chem. Phys. LXXXII, 553 (1992)]. Ab initio computations at the multireference configuration interaction level have been performed to obtain the ground- and several low-lying excited electronic state potential energy curves in three different molecular orientations namely, H+ approaching the O-end and the C-end (collinear), and H+ approaching the CO molecule in perpendicular configuration with fixed CO internuclear distance. Nonadiabatic coupling terms between the ground electronic state (H+ + CO) and the three low-lying excited electronic states (H + CO+) have been computed and the corresponding diabatic potentials have been obtained. A time-dependent wavepacket dynamics study is modeled first involving only the ground and the first excited states and then involving the ground and the three lowest excited states at the collision energy of 9.5 eV. The overall charge transfer probability have been found to be approximately 20%-30% which is in qualitative agreement with the experimental findings.

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.

Three Dimensional Quantum Dynamics of (H-, H2) and Its Isotopic Variants

We present the results of a time-dependent quantum mechanical investigation using centrifugal sudden approximation in the form of reaction probability as a function of collision energy (E(trans)) in the range 0.3-3.0 eV for a range of total angular momentum (J) values and the excitation function sigma(E(trans)) for the exchange reaction H(-) + H(2) (v = 0, j = 0) --> H(2) + H(-) and its isotopic variants in three dimensions on an accurate ab initio potential energy surface published recently (J. Chem. Phys. 2004, 121, 9343). The excitation function results are shown to be in excellent agreement with those obtained from crossed beam measurements by Zimmer and Linder for H(-) + D(2) collisions for energies below the threshold for electron detachment channel and somewhat larger than the most recent results of Haufler et al. for (H(-), D(2)) and (D(-), H(2)) collisions.

Resonances in three-dimensional H+HLi scattering: A time-dependent wave packet dynamical study

This paper examines the resonances in H + HLi scattering. The signature of these resonances is obtained from the oscillations in its reaction probability versus energy curves. They are identified here from a set of pseudospectra calculated for different initial locations of a stationary Gaussian wave packet on the ab initio potential energy surface (PES) reported by Dunne, Murrel, and Jemmer. The nuclear motion on this PES is monitored with the aid of a time-dependent wave packet method and the pseudospectrum are calculated by Fourier transforming the time autocorrelation function of the initial wave packet. The resonances are further examined and assigned by computing their eigenfunctions through spectral quantization algorithm. Both the linewidth as well as decay lifetimes of the resonances are reported.

Quantum scattering calculations on chemical reactions

This review discusses recent quantum scattering calculations on bimolecular chemical reactions in the gas phase. This theory provides detailed and accurate predictions on the dynamics and kinetics of reactions containing three atoms. In addition, the method can now be applied to reactions involving polyatomic molecules. Results obtained with both time-independent and time-dependent quantum dynamical methods are described. The review emphasises the recent development in time-dependent wave packet theories and the applications of reduced dimensionality approaches for treating polyatomic reactions. Calculations on over 40 different reactions are described.