A combined theoretical and experimental investigation of the kinetics and dynamics of the O(1D) + D2 reaction at low temperature

An Experimental and Theoretical Investigation of the Gas-Phase C(3P) + N2O Reaction. Low Temperature Rate Constants and Astrochemical Implications

The reaction between atomic carbon in its ground electronic state, C(³P), and nitrous oxide, N₂O, has been studied below room temperature due to its potential importance for astrochemistry, with both species considered to be present at high abundance levels in a range of interstellar environments. On the experimental side, we measured rate constants for this reaction over the 50-296 K range using a continuous supersonic flow reactor. C(³P) atoms were generated by the pulsed photolysis of carbon tetrabromide at 266 nm and were detected by pulsed laser-induced fluorescence at 115.8 nm. Additional measurements allowing the major product channels to be elucidated were also performed. On the theoretical side, statistical rate theory was used to calculate low temperature rate constants. These calculations employed the results of new electronic structure calculations of the ³A″ potential energy surface of CNNO and provided a basis to extrapolate the measured rate constants to lower temperatures and pressures. The rate constant was found to increase monotonically as the temperature falls (k_C(³P)+N2O (296 K) = (3.4 ± 0.3) × 10^-11 cm³ s^-1), reaching a value of k_C(³P)+N2O (50 K) = (7.9 ± 0.8) × 10^-11 cm³ s^-1 at 50 K. As current astrochemical models do not include the C + N₂O reaction, we tested the influence of this process on interstellar N₂O and other related species using a gas-grain model of dense interstellar clouds. These simulations predict that N₂O abundances decrease significantly at intermediate times (10³ - 10⁵ years) when gas-phase C(³P) abundances are high.

Experimental and theoretical studies of the gas-phase reactions of O(1D) with H2O and D2O at low temperature

Here we report the results of an experimental and theoretical study of the gas-phase reactions between O(¹D) and H₂O and O(¹D) and D₂O at room temperature and below. On the experimental side, the kinetics of these reactions have been investigated over the 50-127 K range using a continuous flow Laval nozzle apparatus, coupled with pulsed laser photolysis and pulsed laser induced fluorescence for the production and detection of O(¹D) atoms respectively. Experiments were also performed at 296 K in the absence of a Laval nozzle. On the theoretical side, the existing full-dimensional ground X ¹A potential energy surface for the H₂O₂ system involved in this process has been reinvestigated and enhanced to provide a better description of the barrierless H-atom abstraction pathway. Based on this enhanced potential energy surface, quasiclassical trajectory calculations and ring polymer molecular dynamics simulations have been performed to obtain low temperature rate constants. The measured and calculated rate constants display similar behaviour above 100 K, showing little or no variation as a function of temperature. Below 100 K, the experimental rate constants increase dramatically, in contrast to the essentially temperature independent theoretical values. The possible origins of the divergence between experiment and theory at low temperatures are discussed.

Benchmarking an improved statistical adiabatic channel model for competing inelastic and reactive processes

Inelastic collisions and elementary chemical reactions proceeding through the formation and subsequent decay of an intermediate collision complex, with an associated deep well on the potential energy surface, pose a challenge for accurate fully quantum mechanical approaches, such as the close-coupling method. In this study, we report on the theoretical prediction of temperature-dependent state-to-state rate coefficients for these complex-mode processes, using a statistical quantum method. This statistical adiabatic channel model is benchmarked by a direct comparison using accurate rate coefficients from the literature for a number of systems (H₂ + H⁺, HD + H⁺, SH⁺ + H, and CH⁺ + H) of interest in astrochemistry and astrophysics. For all of the systems considered, an error of less than factor 2 was found, at least for the dominant transitions and at low temperatures, which is sufficiently accurate for applications in the above mentioned disciplines.

Atom-Diatom Reactive Scattering Collisions in Protonated Rare Gas Systems

The study of the dynamics of atom-diatom reactions involving two rare gas (Rg) atoms and protons is of crucial importance given the astrophysical relevance of these processes. In a series of previous studies, we have been investigating a number of such Rg(1)+ Rg(2)H+→ Rg(2)+ Rg(1)H+ reactions by means of different numerical approaches. These investigations comprised the construction of accurate potential energy surfaces by means of ab initio calculations. In this work, we review the state-of-art of the study of these protonated Rg systems making special emphasis on the most relevant features regarding the dynamical mechanisms which govern these reactive collisions. The aim of this work therefore is to provide an as complete as possible description of the existing information regarding these processes.

Theoretical H + O3 rate coefficients from ring polymer molecular dynamics on an accurate global potential energy surface: assessing experimental uncertainties

Thermal rate coefficients and kinetic isotope effects have been calculated for an important atmospheric reaction H/D + O3 → OH/OD + O2 based on an accurate permutation invariant polynomial-neural network potential energy surface, using ring polymer molecular dynamics (RPMD), quasi-classical trajectory (QCT) and variational transition-state theory (VTST) with multidimensional tunneling. The RPMD approach yielded results that are generally in better agreement with experimental rate coefficients than the VTST and QCT ones, especially at low temperatures, attributable to its capacity to capture quantum effects such as tunneling and zero-point energy. The theoretical results support one group of existing experiments over the other. In addition, rate coefficients for the D + O3 → OD + O2 reaction are also reported using the same methods, which will allow a stringent assessment of future experimental measurements, thus helping to reduce the uncertainty in the recommended rate coefficients of this reaction.

A kinetic study of the N(2D) + C2H4 reaction at low temperature

Electronically excited nitrogen atoms N(²D) are important species in the photochemistry of N₂ based planetary atmospheres such as Titan. Despite this, few N(²D) reactions have been studied over the appropriate low temperature range. During the present work, rate constants were measured for the N(²D) + ethene (C₂H₄) reaction using a supersonic flow reactor at temperatures between 50 K and 296 K. Here, a chemical reaction was used to generate N(²D) atoms, which were detected directly by laser induced fluorescence in the vacuum ultraviolet wavelength region. The measured rate constants displayed very little variation as a function of temperature, with substantially larger values than those obtained in previous work. Indeed, considering an average temperature of 170 K for the atmosphere of Titan leads to a rate constant that is almost seven times larger than the currently recommended value. In parallel, electronic structure calculations were performed to provide insight into the reactive process. While earlier theoretical work at a lower level predicted the presence of a barrier for the N(²D) + C₂H₄ reaction, the present calculations demonstrate that two of the five doublet potential energy surfaces correlating with reagents are likely to be attractive, presenting no barriers for the perpendicular approach of the N atom to the C[double bond, length as m-dash]C bond of ethene. The measured rate constants and new product channels taken from recent dynamical investigations of this process are included in a 1D coupled ion-neutral model of Titan's atmosphere. These simulations indicate that the modeled abundances of numerous nitrogen bearing compounds are noticeably affected by these changes.

Experimental and theoretical studies of the N(2D) + H2 and D2 reactions

This study reports the first kinetic measurements of the N(²D) + H₂, D₂ reactions below 200 K.

Experimental and Theoretical Study of the O(1D) + HD Reaction

This work addresses the kinetics and dynamics of the gas-phase reaction between O(¹D) and HD molecules down to low temperature. Here, measurements were performed by using a supersonic flow (Laval nozzle) reactor coupled with pulsed laser photolysis for O(¹D) production and pulsed-laser-induced fluorescence for O(¹D) detection to obtain rate constants over the 50-300 K range. Additionally, temperature-dependent branching ratios (OD + H/OH + D) were obtained experimentally by comparison of the H/D atom atom yields with those of a reference reaction. In parallel, theoretical rate constants and branching ratios were calculated by using three different techniques; mean potential phase space theory (MPPST), the statistical quantum mechanical method (SQM), and ring polymer molecular dynamics (RPMD). Although the agreement between experimental and theoretical rate constants is reasonably good, with differences not exceeding 30% over the entire temperature range, the theoretical branching ratios derived by the MPPST and SQM methods are as much as 50% larger than the experimental ones. These results are presented in the context of earlier work, while the possible origins of the discrepancies between experiment and theory are discussed.

Theoretical Study of Barrierless Chemical Reactions Involving Nearly Elastic Rebound: The Case of S(1D) + X2, X = H, D

For some values of the total angular momentum consistent with reaction, the title processes involve nonreactive trajectories proceeding through a single rebound mechanism during which the internal motion of the reagent diatom is nearly unperturbed. When such paths are in a significant amount, the classical reaction probability is found to be markedly lower than the quantum mechanical one. This finding was recently attributed to an unusual quantum effect called diffraction-mediated trapping, and a semiclassical correction was proposed in order to take into account this effect in the classical trajectory method. In the present work, we apply the resulting approach to the calculation of opacity functions as well as total and state-resolved integral cross sections (ICSs) and compare the values obtained with exact quantum ones, most of which are new. As the title reactions proceed through a deep insertion well, mean potential statistical calculations are also presented. Seven values of the collision energy, ranging from 30 to 1127 K, are considered. Two remarkable facts stand out: (i) The corrected classical treatment strongly improves the accuracy of the opacity function as compared to the usual classical treatment. When the entrance transition state is tight, however, those trajectories crossing it with a bending vibrational energy below the zero point energy must be discarded. (ii) The quantum opacity function, particularly its cutoff, is finely reproduced by the statistical approach. Consequently, the total ICS is also very well described by the two previous approximate methods. These, however, do not predict state-resolved ICSs with the same accuracy, proving thereby that (i) one or several genuine quantum effects involved in the dynamics are missed by the corrected classical treatment and (ii) the dynamics are not fully statistical.

A low temperature investigation of the N(2D) + CH4, C2H6 and C3H8 reactions

Low temperature rate constants for the N(²D) + C₂H₆, C₃H₈ reactions are shown to be much smaller than previously thought.

Rate constants for the N(2D) + C2H2 reaction over the 50–296 K temperature range

The reactions of metastable atomic nitrogen N(²D) are important processes in the gas-phase chemistry of several planetary atmospheres.

A Ring Polymer Molecular Dynamics Approach to Study the Transition between Statistical and Direct Mechanisms in the H2 + H3+ → H3+ + H2 Reaction

Because of its fundamental importance in astrochemistry, the H₂ + H₃⁺ → H₃⁺ + H₂ reaction has been studied experimentally in a wide temperature range. Theoretical studies of the title reaction significantly lag primarily because of the challenges associated with the proper treatment of the zero-point energy (ZPE). As a result, all previous theoretical estimates for the ratio between a direct proton-hop and indirect exchange (via the H₅⁺ complex) channels deviate from the experiment, in particular, at lower temperatures where the quantum effects dominate. In this work, the ring polymer molecular dynamics (RPMD) method is applied to study this reaction, providing very good agreement with the experiment. RPMD is immune to the shortcomings associated with the ZPE leakage and is able to describe the transition from direct to indirect mechanisms below room temperature. We argue that RPMD represents a useful tool for further studies of numerous ZPE-sensitive chemical reactions that are of high interest in astrochemistry.