Quantum and quasi-classical calculations for the S⁺ + H₂(v,j) → SH⁺(v',j') + H reactive collisions

Diabatic Potential Energy Surfaces of the H2S+ System and the Dynamics Studies of the S+ + H2 (v0 = 2, j0 = 0) Reaction

Global diabatic potential energy surfaces (PESs) of the H2S+ system, corresponding to the 14A″ and 24A″ electronic states, were built by using the neural network method. In ab initio calculations, the aug-cc-pVTZ basis set and MRCI-F12 method were adopted. The topographic features of the new diabatic PESs were discussed and compared with the available theoretical and experimental results in detail. The spectroscopic parameters obtained from the diabatic PESs are in good agreement with previous theoretical and experimental results. Based on the newly constructed diabatic PESs, the nonadiabatic dynamics calculations of the S+ + H2(v0 = 2, j0 = 0) → SH+ + H reaction were carried out using the time-dependent wave packet method. To further understand the nonadiabatic effect, the adiabatic dynamical calculations of the title reaction were also performed based on the adiabatic PES, which was obtained by diagonalizing the diabatic PESs. The deviation between nonadiabatic results and adiabatic values is very obvious. In general, the adiabatic results underestimate the dynamics result at low collision energies and overestimate the dynamics results within high collision energy ranges. Therefore, to obtain accurate dynamics results, the nonadiabatic effect should be included in the calculation.

The Experimental Rate Constant of the S +(2 D) + H 2 Reaction

Endothermic reactions such as S +(4 S) + H 2 are not expected to play a significant role in the chemistry of the interstellar medium (ISM). However, in some specific environments, such as photon-dominated regions (PDR), UV radiation may catalyze the reaction by providing enough internal energy to reactants to overcome endothermicity. For instance, it was recently shown that the vibrational excitation of H2 greatly enhances the reactivity of C+ and S+ with H2, explaining the presence of their respective hydrides CH+ and SH+ in these regions. However, vibrational excitation of H2 is not a unique way to enhance the reactivity by UV radiation. Electronic excitation is an alternative way to effectively inject a huge amount of internal energy into the system, thus favoring reactivity. In this work, we will address how electronic excitation of the sulfur cation can strongly enhance the production of SH+. This is done by measuring experimentally the cross section of the title reaction for collision energies from 50 meV up to several eV and comparing the results with theoretical predictions in the 0.001-3 eV range. The reaction cross section is then used to derive the rate constant for a wide range of temperatures.

The temperature variation of the CH+ + H reaction rate coefficients: a puzzle finally understood?

CH+ was the first molecular ion identified in the interstellar medium and is found to be ubiquitous in interstellar clouds. However, its formation and destruction paths are not well understood, especially at low temperatures. A new theoretical approach based on the canonical variational transition state theory was used to study the H + CH+ reactive collisions. Rate coefficients for formation of C+ ions are calculated as a function of temperature. We considered the participation of a direct path and an indirect path in which the reactants should overcome an entropic barrier to form a van der Waals complex or pass through a CH2+ intermediate complex, respectively. We show that the contribution of both pathways to the formation of C+ has to be taken into account. The new reactive rate coefficients for the title reaction, complemented by reactive data for CH+/CH2+ in the H/H2/He mixture, have been used to simulate the corresponding kinetics experimentally measured using an Atomic Beam 22 Pole Trap apparatus at low temperature. A good agreement with the experimental findings was found at 50 K. At a lower temperature, the model overestimates the formation of C+. This shows that secondary reactions are not responsible for the weak C+ production in the experiments at such temperature. Then, we discuss the possible impact of non-adiabatic effects in the study of the H + CH+ reactive collisions and we found that such effects can be responsible for the decrease of the H + CH+ rate coefficients at low temperature. This study offers an explanation for the disagreement between H + CH+ theoretical and experimental rate coefficients which has been going on for 20 years and highlights the need for performing non-adiabatic studies for this simple chemical reaction.

The role of intersystem crossing in the reactive collision of S+(4S) with H2

We report a study on the reactive collision of S+(4S) with H2, HD, and D2 combining guided ion beam experiments and quantum-mechanical calculations. It is found that the reactive cross sections reflect the existence of two different mechanisms, one being spin-forbidden. Using different models, we demonstrate that the spin-forbidden pathway follows a complex mechanism involving three electronic states instead of two as previously thought. The good agreement between theory and experiment validates the methodology employed and allows us to fully understand the reaction mechanism. This study also provides new fundamental insights into the intersystem crossing process.

Neural network potential energy surface for the low temperature ring polymer molecular dynamics of the H2CO + OH reaction

A new potential energy surface (PES) and dynamical study of the reactive process of H₂CO + OH toward the formation of HCO + H₂O and HCOOH + H are presented. In this work, a source of spurious long range interactions in symmetry adapted neural network (NN) schemes is identified, which prevents their direct application for low temperature dynamical studies. For this reason, a partition of the PES into a diabatic matrix plus a NN many-body term has been used, fitted with a novel artificial neural network scheme that prevents spurious asymptotic interactions. Quasi-classical trajectory (QCT) and ring polymer molecular dynamics (RPMD) studies have been carried on this PES to evaluate the rate constant temperature dependence for the different reactive processes, showing good agreement with the available experimental data. Of special interest is the analysis of the previously identified trapping mechanism in the RPMD study, which can be attributed to spurious resonances associated with excitations of the normal modes of the ring polymer.

A new global potential energy surface of the SH2+(X4A'') system and quantum calculations for the S+ + H2(v = 0-3, j = 0) reaction

A new global potential energy surface (PES) for the ground state of the SH2+(X4A'') system is constructed using a permutation invariant polynomial neural network method. In ab initio calculations, the MRCI-F12 method with the AVTZ basis set is used. Furthermore, the dynamics calculations of the S+ + H2(v = 0-3, j = 0) → SH+ + H reaction are carried out based on the new PES. The reaction probabilities and integral cross sections are compared with available theoretical calculations. Present values are in general good agreement with the previous theoretical studies. However, some discrepancies can still be found due to different PESs used in the calculation. Furthermore, the vibrational energy of the reactant molecule can significantly enhance the reactivity compared to the translational energy. The differential cross sections indicated that the reaction mechanism is changed from the "head-on" rebound mechanism to the tripping mechanism with the increasing number of initial vibrational excitation state.

Signature of a conical intersection in the dissociative photoionization of formaldehyde

Electron/ion coincidence experiments and ab initio calculations of the dissociative photoionization of formaldehyde reveal the presence of a conical intersection controlling the dynamics and favoring dissociation into the molecular channel, CO⁺ + H₂.

Formation and Destruction of SiS in Space

The presence of SiS in space seems to be restricted to a few selected types of astronomical environments. It is long known to be present in circumstellar envelopes around evolved stars and it has also been detected in a handful of star-forming regions with evidence of outflows, like Sgr B2, Orion KL and more recently L1157-B1. The kinetics of reactions involving SiS is very poorly known and here we revisit the chemistry of SiS in space by studying some potentially important reactions of formation and destruction of this molecule. We calculated ab initio potential energy surfaces of the SiOS system and computed rate coefficients in the temperature range 50-2500 K for the reaction of destruction of SiS, in collisions with atomic O, and of its formation, through the reaction between Si and SO. We find that both reactions are rapid, with rate coefficients of a few times 10-10 cm3 s-1, almost independent of temperature. In the reaction between Si and SO, SiO production is 5-7 times more efficient than SiS formation. The reaction of SiS with O atoms can play an important role in destroying SiS in envelopes around evolved stars. We built a simple chemical model of a postshock gas to study the chemistry of SiS in protostellar outflows and we found that SiS forms with a lower abundance and later than SiO, that SiS is efficiently destroyed through reaction with O, and that the main SiS-forming reactions are Si + SO and Si + SO2.

Full dimensional potential energy surface and low temperature dynamics of the H2CO + OH → HCO + H2O reaction

A new method is proposed to analytically represent the potential energy surface of reactions involving polyatomic molecules capable of accurately describing long-range interactions and saddle points, needed to describe low-temperature collisions. It is based on two terms, a reactive force field term and a many-body term. The reactive force field term accurately describes the fragments, long-range interactions among them and the saddle points for reactions. The many-body term increases the desired accuracy everywhere else. This method has been applied to the OH + H2CO → H2O + HCO reaction, giving a barrier of 27.4 meV. The simulated classical rate constants with this potential are in good agreement with recent experimental results [Ocaña et al., Astrophys. J., 2017, submitted], showing an important increase at temperatures below 100 K. The reaction mechanism is analyzed in detail here, and explains the observed behavior at low energy by the formation of long-lived collision complexes, with roaming trajectories, with a capture observed for very long impact parameters, >100 a.u., determined by the long-range dipole-dipole interaction.

State-to-state chemistry and rotational excitation of CH+ in photon-dominated regions

We present a detailed theoretical study of the rotational excitation of CH+ due to reactive and nonreactive collisions involving C+(2P), H2, CH+, H and free electrons. Specifically, the formation of CH+ proceeds through the reaction between C+(2P) and H2(νH2 = 1, 2), while the collisional (de)excitation and destruction of CH+ is due to collisions with hydrogen atoms and free electrons. State-to-state and initial-state-specific rate coefficients are computed in the kinetic temperature range 10-3000 K for the inelastic, exchange, abstraction and dissociative recombination processes using accurate potential energy surfaces and the best scattering methods. Good agreement, within a factor of 2, is found between the experimental and theoretical thermal rate coefficients, except for the reaction of CH+ with H atoms at kinetic temperatures below 50 K. The full set of collisional and chemical data are then implemented in a radiative transfer model. Our Non-LTE calculations confirm that the formation pumping due to vibrationally excited H2 has a substantial effect on the excitation of CH+ in photon-dominated regions. In addition, we are able to reproduce, within error bars, the far-infrared observations of CH+ toward the Orion Bar and the planetary nebula NGC 7027. Our results further suggest that the population of νH2 = 2 might be significant in the photon-dominated region of NGC 7027.