Abundance Ratios of OH/CO and HCO+/CO as Probes of the Cosmic-Ray Ionization Rate in Diffuse Clouds

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 H2O+, 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) <= 10(-15) s(-1), which are potentially useful for modeling 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 Atacama Large Millimeter/submillimeter Array, we have derived the H-2 number densities (n(H2)) toward 4 lines of sight; assuming a kinetic temperature of T-k = 50 K, we find a range of (0.14 +/- 0.03-1.2 +/- 0.1) x 10(2) cm(-3). By comparing the observed and modeled HCO+/CO ratios, we find that zeta(2) in our diffuse gas sample is in the range of 1.0(-1.0)(+14.8) x 10(-16)-2.5(-2.4)(+1.4) x 10(-15) s(-1). This is similar to 2 times higher than the average value measured at higher extinction, supporting an attenuation of CRs as suggested by theoretical models.