Impact of cosmic-ray propagation on the chemistry and ionisation fraction of dark clouds

TL;DR

This paper introduces an efficient numerical framework to model cosmic-ray propagation in dark clouds, revealing its significant impact on chemical species and ionisation, and providing practical formulas for estimating ionisation rates from observations.

Contribution

The authors develop a novel, computationally efficient method to simulate cosmic-ray propagation and its effects on cloud chemistry, improving upon previous models.

Findings

Cosmic-ray propagation significantly influences deuterated and non-deuterated species.

An analytical formula links ionisation fraction to observable tracers.

A linear fit estimates cosmic-ray ionisation rate from local H2 density.

Abstract

A proper modelling of the cosmic-ray ionisation rate within gas clouds is crucial to describe their chemical evolution accurately. However, this modelling is computationally demanding because it requires the propagation of cosmic rays throughout the cloud over time. We present a more efficient approach that simultaneously guarantees a reliable estimate of the cosmic-ray impact on the chemistry of prestellar cores. We introduce a numerical framework that mimics the cosmic-ray propagation within gas clouds and applies it to magnetohydrodynamic simulations performed with the code GIZMO. It simulates the cosmic-ray attenuation by computing the effective column density of H$_2$ that is traversed, which is estimated using the same kernel weighting approach as employed in the simulation. The obtained cosmic-ray ionisation rate is then used in post-processing to study the chemical evolution of the clouds. We found that cosmic-ray propagation affects deuterated and non-deuterated species significantly and that it depends on the assumed cosmic-ray spectrum. We explored correlations between the electron abundance, the cosmic-ray ionisation rate, and the abundance of the most relevant ions (HCO$^+$, N$_2$H$^+$, DCO$^+$, N$_2$D$^+$, and o-H$_2$D$^+$), with the purpose of finding simple expressions that link them. We provide an analytical formula to estimate the ionisation fraction, X(e$^-$), from observable tracers and applied it to existing observations of high-mass clumps. We obtained values of about 10$^{-8}$, which is in line with previous works and with expectations for dense clouds. We also provide a linear fit to calculate the cosmic-ray ionisation rate from the local H$_2$ density, which is to be employed in three-dimensional simulations that do not include cosmic-ray propagation.