HD-H+ collisions: statistical and quantum state-to-state studies

In the early Universe, the cooling mechanisms of the gas significantly rely on the HD abundance and excitation conditions. A proper modeling of its formation and destruction paths as well as its excitation by both radiative and collisional processes is then required to accurately describe the cooling mechanisms of the pristine gas. In such media, ion-molecule reactions are dominant. Their theoretical study is challenging and state-of-the-art quantum time-independent methods are computationally limited to collisions involving light molecules. Here, we report a state-to-state scattering study of the HD-H⁺ collisional system using two different methods: an exact quantum time-independent approach and a recently developed fast and efficient statistical method. Reactive and inelastic rate coefficients were obtained for temperatures up to 300 K. The statistical method is able to reproduce exact calculations with an accuracy reaching the astrophysical needs while drastically reducing the computational resources requirements. Such results suggest that this new statistical method should be considered to provide the astrophysical community collisional data for which quantum calculations are impossible.

Dynamics of the isotope exchange reaction of D with H3 +, H2D+, and D2H. J Chem Phys 2021;154:084307. [PMID: 33639774 DOI: 10.1063/5.0038434]

We have measured the merged-beams rate coefficient for the titular isotope exchange reactions as a function of the relative collision energy in the range of ∼3 meV-10 eV. The results appear to scale with the number of available sites for deuteration. We have performed extensive theoretical calculations to characterize the zero-point energy corrected reaction path. Vibrationally adiabatic minimum energy paths were obtained using a combination of unrestricted quadratic configuration interaction of single and double excitations and internally contracted multireference configuration interaction calculations. The resulting barrier height, ranging from 68 meV to 89 meV, together with the various asymptotes that may be reached in the collision, was used in a classical over-the-barrier model. All competing endoergic reaction channels were taken into account using a flux reduction factor. This model reproduces all three experimental sets quite satisfactorily. In order to generate thermal rate coefficients down to 10 K, the internal excitation energy distribution of each H₃ ⁺ isotopologue is evaluated level by level using available line lists and accurate spectroscopic parameters. Tunneling is accounted for by a direct inclusion of the exact quantum tunneling probability in the evaluation of the cross section. We derive a thermal rate coefficient of <1×10^-12 cm³ s^-1 for temperatures below 44 K, 86 K, and 139 K for the reaction of D with H₃ ⁺, H₂D⁺, and D₂H⁺, respectively, with tunneling effects included. The derived thermal rate coefficients exceed the ring polymer molecular dynamics prediction of Bulut et al. [J. Phys. Chem. A 123, 8766 (2019)] at all temperatures.

Kinetics Of The H + CH2 → CH + H2 Reaction At Low Temperature

A quasiclassical trajectory study of the kinetics of the title astrochemical reaction in a range of temperatures varying from 5 to 1000 K (corresponding to both the outer and the inner regions of the protostar and the circumstellar envelopes) was carried out and a clear dependence of the rate coefficient on the temperature is given, in contrast with the constant value adopted in kinetics astrochemical databases. Levering the massive nature of the performed calculations and of the detailed dynamical investigation of the reactive process, a rationalization of the temperature dependence of the released translational energy and of the rovibrational population of the CH and H₂ diatomic products is also provided. Furthermore, the effect of the initial rovibrational energy of CH₂ on the state-specific rate coefficients and cross sections is analyzed in order to single out the role played by the different regions of the potential energy surface on the dynamical outcomes and on the modeling the temperature dependence of the reactive efficiency of the investigated process. This led to a parametrization of the computed rate in terms of the following double Arrhenius expression (in cm³ s^-1), k(T) = 2.50 × 10^-10 exp(- 1.67/T) + 5.98 × 10^-11 exp(- 280.5/T), alternative to the piecewise formulation into the three subintervals of temperature in which the overall 5-1000 K interval can be divided.

The low temperature D+ + H2→ HD + H+ reaction rate coefficient: a ring polymer molecular dynamics and quasi-classical trajectory study

The reaction between D⁺ and H₂ plays an important role in astrochemistry at low temperatures and also serves as a prototype for a simple ion-molecule reaction. Its ground X[combining tilde]¹A' state has a very small thermodynamic barrier (up to 1.8 × 10^-2 eV) and the reaction proceeds through the formation of an intermediate complex lying within the potential well with a depth of at least 0.2 eV, thus representing a challenge for dynamical studies. In the present work, we analyze the title reaction within the temperature range of 20-100 K by means of ring polymer molecular dynamics (RPMD) and quasi-classical trajectory (QCT) methods over the full-dimensional global potential energy surface developed by Aguado et al. [A. Aguado, O. Roncero, C. Tablero, C. Sanz and M. Paniagua, J. Chem. Phys., 2000, 112, 1240]. The computed thermal RPMD and QCT rate coefficients are found to be almost independent of temperature and fall within the range of 1.34-2.01 × 10^-9 cm³ s^-1. They are also in very good agreement with previous time-independent quantum mechanical and statistical quantum method calculations. Furthermore, we observe that the choice of asymptotic separation distance between the reactants can markedly alter the rate coefficient in the low temperature regime (20-50 K). Therefore it is of utmost importance to correctly assign the value of this parameter for dynamical studies, particularly at very low temperatures of astrochemical importance. We finally conclude that the experimental rate measurements for the title reaction are highly desirable in future.

Angular momentum-scattering angle quantum correlation: a generalized deflection function

A natural generalization of the classical deflection function, the functional dependence of the deflection angle on the angular momentum (or the impact parameter), is the joint probability density function of these two quantities, revealing the correlation between them. It provides, at a glance, detailed information about the reaction mechanisms and how changes in the impact parameter affect the product angular distribution. It is also useful to predict the presence of quantum phenomena such as interference. However, the classical angular momentum-scattering angle correlation function has a limited use whenever quantum effects become important. Rigorously speaking, there is not a quantum equivalent of the classical joint distribution, as the differential cross section depends on the coherences between the different values of J caused by the cross terms in the expansion of partial waves. In this article, we present a simple method to calculate a quantum analog of this correlation, a generalized deflection function that can shed light onto the reaction mechanism using just quantum mechanical results. Our results show that there is a very good agreement between the quantum and classical correlation functions as long as quantum effects are not all relevant. When this is not the case, it will also be shown that the quantum correlation function is most useful to observe the extent of quantum effects such as interference among different reaction mechanisms.

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+.

A comparative account of quantum dynamics of the H⁺ + H₂ reaction at low temperature on two different potential energy surfaces

Rotationally resolved reaction probabilities, integral cross sections, and rate constant for the H(+) + H2 (v = 0, j = 0 or 1) → H2 (v' = 0, j') + H(+) reaction are calculated using a time-independent quantum mechanical method and the potential energy surface of Kamisaka et al. [J. Chem. Phys. 116, 654 (2002)] (say KBNN PES). All partial wave contributions of the total angular momentum, J, are included to obtain converged cross sections at low collision energies and rate constants at low temperatures. In order to test the accuracy of the KBNN PES, the results obtained here are compared with those obtained in our earlier work [P. Honvault et al., Phys. Rev. Lett. 107, 023201 (2011)] using the accurate potential energy surface of Velilla et al. [J. Chem. Phys. 129, 084307 (2008)]. Integral cross sections and rate constants obtained on the two potential energy surfaces considered here show remarkable differences in terms of magnitude and dependence on collision energy (or temperature) which can be attributed to the differences observed in the topography of the surfaces near to the entrance channel. This clearly shows the inadequacy of the KBNN PES for calculations at low collision energies.

The D(+) + H2 reaction: differential and integral cross sections at low energy and rate constants at low temperature

The D(+) + H2 reaction is investigated by means of a time independent quantum mechanical (TIQM) and statistical quantum mechanical (SQM) methods. Differential cross sections and product rotational distributions obtained with these two theoretical approaches for collision energies between 1 meV and 0.1 eV are compared to analyze the dynamics of the process. The agreement observed between the TIQM differential cross sections and the SQM predictions as the energy increases revealed the role played by the complex-forming mechanism. The importance of a good description of the asymptotic regions is also investigated by calculating rate constants for the title reaction at low temperature.

Dynamics of the D(+) + H2 → HD + H(+) reaction at the low energy regime by means of a statistical quantum method

The D(+) +H2(v = 0, j = 0, 1) → HD+H(+) reaction has been investigated at the low energy regime by means of a statistical quantum mechanical (SQM) method. Reaction probabilities and integral cross sections (ICSs) between a collisional energy of 10(-4) eV and 0.1 eV have been calculated and compared with previously reported results of a time independent quantum mechanical (TIQM) approach. The TIQM results exhibit a dense profile with numerous narrow resonances down to Ec ~ 10(-2) eV and for the case of H2(v = 0, j = 0) a prominent peak is found at ~2.5 × 10(-4) eV. The analysis at the state-to-state level reveals that this feature is originated in those processes which yield the formation of rotationally excited HD(v' = 0, j' > 0). The statistical predictions reproduce reasonably well the overall behaviour of the TIQM ICSs at the larger energy range (Ec ≥ 10(-3) eV). Thermal rate constants are in qualitative agreement for the whole range of temperatures investigated in this work, 10-100 K, although the SQM values remain above the TIQM results for both initial H2 rotational states, j = 0 and 1. The enlargement of the asymptotic region for the statistical approach is crucial for a proper description at low energies. In particular, we find that the SQM method leads to rate coefficients in terms of the energy in perfect agreement with previously reported measurements if the maximum distance at which the calculation is performed increases noticeably with respect to the value employed to reproduce the TIQM results.

Rate coefficients from quantum and quasi-classical cumulative reaction probabilities for the S(1D) + H2 reaction

Cumulative reaction probabilities (CRPs) at various total angular momenta have been calculated for the barrierless reaction S((1)D) + H(2) → SH + H at total energies up to 1.2 eV using three different theoretical approaches: time-independent quantum mechanics (QM), quasiclassical trajectories (QCT), and statistical quasiclassical trajectories (SQCT). The calculations have been carried out on the widely used potential energy surface (PES) by Ho et al. [J. Chem. Phys. 116, 4124 (2002)] as well as on the recent PES developed by Song et al. [J. Phys. Chem. A 113, 9213 (2009)]. The results show that the differences between these two PES are relatively minor and mostly related to the different topologies of the well. In addition, the agreement between the three theoretical methodologies is good, even for the highest total angular momenta and energies. In particular, the good accordance between the CRPs obtained with dynamical methods (QM and QCT) and the statistical model (SQCT) indicates that the reaction can be considered statistical in the whole range of energies in contrast with the findings for other prototypical barrierless reactions. In addition, total CRPs and rate coefficients in the range of 20-1000 K have been calculated using the QCT and SQCT methods and have been found somewhat smaller than the experimental total removal rates of S((1)D).

State specific stabilization of H+ + H2(j) collision complexes

Stabilization of H3(+) collision complexes has been studied at nominal temperatures between 11 and 33 K using a 22-pole radio frequency (rf) ion trap. Apparent binary rate coefficients, k(*) = kr + k3[H2], have been measured for para- and normal-hydrogen at number densities between some 10(11) and 10(14) cm(-3). The state specific rate coefficients extracted for radiative stabilization, kr(T;j), are all below 2 × 10(-16) cm(3) s(-1). There is a slight tendency to decrease with increasing temperature. In contrast to simple expectations, kr(11 K;j) is for j = 0 a factor of 2 smaller than for j = 1. The ternary rate coefficients for p-H2 show a rather steep T-dependence; however, they are increasing with temperature. The state specific ternary rate coefficients, k3(T;j), measured for j = 0 and derived for j = 1 from measurements with n-H2, differ by an order of magnitude. Most of these surprising observations are in disagreement with predictions from standard association models, which are based on statistical assumptions and the separation of complex formation and competition between stabilization and decay. Most probably, the unexpected collision dynamics are due to the fact that, at the low translational energies of the present experiment, only a small number of partial waves participate. This should make exact quantum mechanical calculations of kr feasible. More complex is three-body stabilization, because it occurs on the H5(+) potential energy surface.

State-to-state quantum mechanical calculations of rate coefficients for the D+ + H2 → HD + H+ reaction at low temperature

The dynamics of the D(+) + H2 → HD + H(+) reaction on a recent ab initio potential energy surface (Velilla, L.; Lepetit, B.; Aguado, A.; Beswick, J. A.; Paniagua, M. J. Chem. Phys. 2008, 129, 084307) has been investigated by means of a time-independent quantum mechanical approach. Cross-sections and rate coefficients are calculated, respectively, for collision energies below 0.1 eV and temperatures up to 100 K for astrophysical application. An excellent accord is found for collision energy above 5 meV, while a disagreement between theory and experiment is observed below this energy. We show that the rate coefficients reveal a slightly temperature-dependent behavior in the upper part of the temperature range considered here. This is in agreement with the experimental data above 80 K, which give a temperature independent value. However, a significant decrease is found at temperatures below 20 K. This decrease can be related to quantum effects and the decay back to the reactant channel, which are not considered by simple statistical approaches, such as the Langevin model. Our results have been fitted to appropriate analytical expressions in order to be used in astrochemical and cosmological models.

Stabilization of H+-H2 collision complexes between 11 and 28 K

Formation of H(3)(+) via association of H(+) with H(2) has been studied at low temperatures using a 22-pole radiofrequency trap. Operating at hydrogen number densities from 10(11) to 10(14) cm(-3), the contributions of radiative, k(r), and ternary, k(3), association have been extracted from the measured apparent binary rate coefficients, k*=k(r)+k(3)[H(2)]. Surprisingly, k(3) is constant between 11 and 22 K, (2.6±0.8)×10(-29) cm(6) s(-1), while radiative association decreases from k(r)(11 K)=(1.6±0.3)×10(-16) cm(3) s(-1) to k(r)(28 K)=(5±2)×10(-17) cm(3) s(-1). These results are in conflict with simple association models in which formation and stabilization of the complex are treated separately. Tentative explanations are based on the fact that, at low temperatures, only few partial waves contribute to the formation of the collision complex and that ternary association with H(2) may be quite inefficient because of the 'shared proton' structure of H(5)(+).