Quantum Dynamics of the Reaction H(2S) + HeH+(X1Σ+) → H2+(X2Σg+) + He(1S) from Cold to Hyperthermal Energies: Time-Dependent Wavepacket Study and Comparison with Time-Independent Calculations

Experimental study of the proton-transfer reaction C + H2+ → CH+ + H and its isotopic variant (D2+)

We report absolute integral cross section (ICS) measurements using a dual-source merged-fast-beams apparatus to study the titular reactions over the relative translational energy range of Er ∼ 0.01-10 eV. We used photodetachment of C- to produce a pure beam of atomic C in the ground electronic 3P term, with statistically populated fine-structure levels. The H2+ and D2+ were formed in an electron impact ionization source, with well known vibrational and rotational distributions. The experimental work is complemented by a theoretical study of the CH2+ electronic system in the reactant and product channels, which helps to clarify the possible reaction mechanisms underlying the ICS measurements. Our measurements provide evidence that the reactions are barrierless and exoergic. They also indicate the apparent absence of an intermolecular isotope effect, to within the total experimental uncertainties. Capture models, taking into account either the charge-induced dipole interaction potential or the combined charge-quadrupole and charge-induced dipole interaction potentials, produce reaction cross sections that lie a factor of ∼4 above the experimental results. Based on our theoretical study, we hypothesize that the reaction is most likely to proceed adiabatically through the 14A' and 14A'' states of CH2+via the reaction C(3P) + H2+(2Σ+g) → CH+(3Π) + H(2S). We also hypothesize that at low collision energies only H2+(v ≤ 2) and D2+(v ≤ 3) contribute to the titular reactions, due to the onset of dissociative charge transfer for higher vibrational v levels. Incorporating these assumptions into the capture models brings them into better agreement with the experimental results. Still, for energies ⪅0.1 eV where capture models are most relevant, the modified charge-induced dipole model yields reaction cross sections with an incorrect energy dependence and lying ∼10% below the experimental results. The capture cross section obtained from the combined charge-quadrupole and charge-induced dipole model better matches the measured energy dependence but lies ∼30-50% above the experimental results. These findings provide important guidance for future quasiclassical trajectory and quantum mechanical treatments of this reaction.

Dynamics of H + HeH+(v = 0, j = 0) → H2+ + He: Insight on the Possible Complex-Forming Behavior of the Reaction

The H + HeH⁺→ He + H₂⁺ reaction has been studied by means of a combination of theoretical approaches: a statistical quantum method (SQM), ring polymer molecular dynamics (RPMD), and the quasiclassical trajectory (QCT) method. Cross sections and rate constants have been calculated in an attempt to investigate the dynamics of the process. The comparison with previous calculations and experimental results reveals that despite the fact that statistical predictions seem to reproduce some of the overall observed features, the analysis at a more detailed state-to-state level shows noticeable deviations from a complex-forming dynamics. We find some differences in cross sections and rate constants obtained in the QCT calculation with a Gaussian binning procedure with respect to previous works in which the standard histogram binning was employed.

Scattering study of the Ne + NeH(+)(v0 = 0, j0 = 0) → NeH(+) + Ne reaction on an ab initio based analytical potential energy surface

Initial state selected dynamics of the Ne + NeH(+)(v0 = 0, j0 = 0) → NeH(+) + Ne reaction is investigated by quantum and statistical quantum mechanical (SQM) methods on the ground electronic state. The three-body ab initio energies on a set of suitably chosen grid points have been computed at CCSD(T)/aug-cc-PVQZ level and analytically fitted. The fitting of the diatomic potentials, computed at the same level of theory, is performed by spline interpolation. A collinear [NeHNe](+) structure lying 0.72 eV below the Ne + NeH(+) asymptote is found to be the most stable geometry for this system. Energies of low lying vibrational states have been computed for this stable complex. Reaction probabilities obtained from quantum calculations exhibit dense oscillatory structures, particularly in the low energy region and these get partially washed out in the integral cross section results. SQM predictions are devoid of oscillatory structures and remain close to 0.5 after the rise at the threshold thus giving a crude average description of the quantum probabilities. Statistical cross sections and rate constants are nevertheless in sufficiently good agreement with the quantum results to suggest an important role of a complex-forming dynamics for the title reaction.

Complementarity between Quantum and Classical Mechanics in Chemical Modeling. The H + HeH+ → H2 + + He Reaction: A Rigourous Test for Reaction Dynamics Methods

In this work we present a dynamical study of the H + HeH+ → H2+ + He reaction in a collision energy range from 0.1 meV to 10 eV, suitable to be used in applicative models. The paper extends and complements a recent work [ Phys. Chem. Chem. Phys. 2014, 16, 11662] devoted to the characterization of the reactivity from the ultracold regime up to the three-body dissociation breakup. In particular, the accuracy of the quasi-classical trajectory method below the three-body dissociation threshold has been assessed by a detailed comparison with previous calculations performed with different reaction dynamics methods, whereas the reliability of the results in the high energy range has been checked by a direct comparison with the available experimental data. Integral cross sections for several HeH+ roto-vibrational states have been analyzed and used to understand the extent of quantum effects in the reaction dynamics. By using the quasi-classical trajectory method and quantum mechanical close coupling data, respectively, in the high and low collision energy ranges, we obtain highly accurate thermal rate costants until 15 000 K including all (178) the roto-vibrational bound and quasi-bound states of HeH+. The role of the collision-induced dissociation is also discussed and explicitly calculated for the ground roto-vibrational state of HeH+.

Quantum, Statistical, and Quasiclassical Trajectory Studies For the Ne + HeH(+) → NeH(+) + He Reaction on the Ground Electronic State

Real wave packet, statistical quantum, and quasiclassical trajectory methods were employed to study the dynamics of Ne + HeH(+)(v0,j0) → He + NeH(+) reaction on an ab initio potential energy surface [J. Phys. Chem. A 2013, 117, 13070-13078]. Quantum and statistical quantum calculations were performed within the centrifugal sudden (CS) approximation as well as including the Coriolis coupling (CC). Dense oscillatory structures of the quantum reaction probabilities and fair agreement between quantum and statistical cross sections suggest a complex forming mechanism for the reaction. No significant differences between cross sections obtained within the CS and CC approaches are observed. Quasiclassical trajectory results give an excellent average description of the quantum CC results. At low collision energies, there is a substantial decrease in reactivity for the reaction upon rovibrational excitation. Initial state selected rate constants for the title reaction are calculated between 20 and 1000 K, and the calculated value at 300 K agrees quite well with the available experimental result. Reaction cross sections and rate constants are also compared with those calculated via the Langevin capture model for exothermic reactions.

Born–Oppenheimer and Renner–Teller coupled-channel quantum reaction dynamics of O(3P) + H2+(X2Σg+) collisions

We present Born–Oppenheimer (BO) and Renner–Teller (RT) time dependent quantum dynamics studies of the reactions O(³P) + H₂⁺(X²Σ_g⁺) → OH⁺(X³Σ⁻) + H(²S) and OH(X²Π) + H⁺.