Low temperature rate coefficients of the H + CH(+) → C(+) + H2 reaction: New potential energy surface and time-independent quantum scattering

The observed abundances of the methylidyne cation, CH(+), in diffuse molecular clouds can be two orders of magnitude higher than the prediction of the standard gas-phase models which, in turn, predict rather well the abundances of neutral CH. It is therefore necessary to investigate all the possible formation and destruction processes of CH(+) in the interstellar medium with the most abundant species H, H2, and e(-). In this work, we address the destruction process of CH(+) by hydrogen abstraction. We report a new calculation of the low temperature rate coefficients for the abstraction reaction, using accurate time-independent quantum scattering and a new high-level ab initio global potential energy surface including a realistic model of the long-range interaction between the reactants H and CH(+). The calculated thermal rate coefficient is in good agreement with the experimental data in the range 50 K-800 K. However, at lower temperatures, the experimental rate coefficient takes exceedingly small values which are not reproduced by the calculated rate coefficient. Instead, the latter rate coefficient is close to the one given by the Langevin capture model, as expected for a reaction involving an ion and a neutral species. Several recent theoretical works have reported a seemingly good agreement with the experiment below 50 K, but an analysis of these works show that they are based on potential energy surfaces with incorrect long-range behavior. The experimental results were explained by a loss of reactivity of the lowest rotational states of the reactant; however, the quantum scattering calculations show the opposite, namely, a reactivity enhancement with rotational excitation.

Spin-orbit quenching of the C+(2P) ion by collisions with para- and ortho-H2

Spin-orbit (de-)excitation of C(+)((2)P) by collisions with H2, a key process for astrochemistry, is investigated. Quantum-mechanical calculations of collisions between C(+) ions and para- and ortho-H2 have been performed in order to determine the cross section for the C(+) (2)P3∕2 → (2)P1∕2 fine-structure transition at low and intermediate energies. The calculation are based on new ab initio potential energy surfaces obtained using the multireference configuration interaction method. Corresponding rate coefficients were obtained for temperatures ranging from 5 to 500 K. These rate coefficients are compared to previous estimations, and their impact is assessed through radiative transfer computation. They are found to increase the flux of the (2)P3∕2 → (2)P1∕2 line at 158 μm by up to 30% for typical diffuse interstellar cloud conditions.