# Abundance ratios of OH/CO and HCO+/CO as probes of the cosmic ray ionization rate in diffuse clouds

## TL;DR

This paper introduces a new method to estimate the cosmic-ray ionization rate in diffuse clouds using the abundance ratios of OH/CO and HCO+/CO, providing an alternative to traditional tracers that are difficult to detect.

## Contribution

It proposes an analytical approach linking molecular abundance ratios to the cosmic-ray ionization rate, validated with ALMA observations, enhancing modeling of interstellar chemistry.

## Key findings

- Abundance ratios correlate with cosmic-ray ionization rate $\
- Derived CRIR values are higher than those at higher extinction, indicating CR attenuation.
- Analytical expressions accurately predict OH abundance for $\

## Abstract

The cosmic-ray ionization rate (CRIR, $\zeta_2$) is one of the key parameters controlling the formation and destruction of various molecules in molecular clouds. However, the current most commonly used CRIR tracers, such as H$_3^+$, OH$^+$, and H$_2$O$^+$, are hard to detect and require the presence of background massive stars for absorption measurements. In this work, we propose an alternative method to infer the CRIR in diffuse clouds using the abundance ratios of OH/CO and HCO$^+$/CO. We have analyzed the response of chemical abundances of CO, OH, and HCO$^+$ on various environmental parameters of the interstellar medium in diffuse clouds and found that their abundances are proportional to $\zeta_2$. Our analytic expressions give an excellent calculation of the abundance of OH for $\zeta_2$ $\leq$10$^{-15}$ s$^{-1}$, which are potentially useful for modelling chemistry in hydrodynamical simulations. The abundances of OH and HCO$^+$ were found to monotonically decrease with increasing density, while the CO abundance shows the opposite trend. With high-sensitivity absorption transitions of both CO (1--0) and (2--1) lines from ALMA, we have derived the H$_2$ number densities ($n_{\rm H_2}$) toward 4 line-of-sights (LOSs); assuming a kinetic temperature of $T_{\rm k}=50\,{\rm K}$, we find a range of (0.14$\pm$0.03--1.2$\pm$0.1)$\times$10$^2$ cm$^{-3}$}. By comparing the observed and modelled HCO$^+$/CO ratios, we find that $\zeta_2$ in our diffuse gas sample is in the { range of $1.0_{-1.0}^{+14.8}$ $\times$10$^{-16}- 2.5_{-2.4}^{+1.4}$ $\times$10$^{-15}$ s$^{-1}$. This is $\sim$2 times higher than the average value measured at higher extinction, supporting an attenuation of CRs as suggested by theoretical models.