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

Theoretical investigation of the H + HD → D + H2 chemical reaction for astrophysical applications: A state-to-state quasi-classical study

We report a large set of state-to-state rate constants for the H + HD reactive collision, using Quasi-Classical Trajectory (QCT) simulations on the accurate H₃ 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 HD (up to ≈50 000 K), relevant to various non thermal equilibrium astrophysical media, are considered. We have validated the accuracy of our QCT calculations with a new efficient adaptation of the Multi Configuration Time Dependent Hartree (MCTDH) method to compute the reaction probability of a specific reactive channel. Our study has revealed that the high temperature regime favors the production of H₂ in its highly rovibrationnally excited states, which can de-excite radiatively (cooling the gas) or collisionally (heating the gas). Those new state-to-state QCT reaction rate constants represent a significant improvement in our understanding of the possible mechanisms leading to the destruction of HD by its collision with a H atom.

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.