The impact of hydrogen atom tunneling on aromatic chemistry in TMC-1

TL;DR

This study investigates how hydrogen atom tunneling influences aromatic molecule chemistry in interstellar environments, identifying key reactions and assessing their impact on aromatic abundances through detailed kinetic analyses.

Contribution

The paper provides new potential energy surfaces and rate coefficients for hydrogen abstraction reactions relevant to astrochemistry, highlighting the significance of tunneling effects in the ISM.

Findings

Reactions with C₂H, OH, and CN are competitive in the ISM despite low rate coefficients.

The NH₂ reaction is inefficient at typical interstellar conditions.

Reaction involving c-C₆H₅⁺ significantly affects aromatic abundance predictions.

Abstract

Hydrogen atom tunneling likely plays a substantial role in the gas-phase chemistry of astrochemical environments. To determine the potential effect that it has on the chemical modeling of aromatic molecules, we screened the kida.uva.2024 network, and our own expanded network to find reactions which could be significantly accelerated by hydrogen atom tunneling in the ISM. In total, 64 reactions were identified. The hydrogen abstraction reactions from H$_{2}$ to four key interstellar radicals (C$_{2}$H, OH, CN, and NH$_{2}$) were studied further using newly calculated potential energy surfaces and RRKM analyses to determine rate coefficients for a temperature of 10 K and a density of 2 $\times$ 10$^{4}$ cm$^{-3}$. Despite having low rate coefficients of 1.66 $\times$ 10$^{-15}$, 8.17 $\times$ 10$^{-16}$ and 3.15 $\times$ 10$^{-16}$ $\mathrm{cm^{3}\,s^{-1}}$ the C$_{2}$H, OH, and CN reactions are competitive in the ISM, due to large overall rates caused by the high abundance of molecular hydrogen. The calculated value for the NH$_{2}$ reaction, however, was much smaller and found to be inefficient at ISM conditions. The possible effects of all other considered reactions were studied with simulations using calculated collision limit rate coefficients. Upper and lower bounds were then placed on modeled aromatic abundances using the most significant reactions. Due to the dependence of calculated aromatic abundances on reactions involving c-C$_{6}$H$_{5}^{+}$ and the recent questions surrounding its reactivity, we also explored the abundance variations caused by reactions leading to or involving c-C$_{6}$H$_{5}^{+}$.