2016ApJ...830..151S

The chemistry of dense interstellar regions was analyzed using a time-dependent gas-grain astrochemical simulation and a new chemical network that incorporates deuterated chemistry, taking into account nuclear spin states for the hydrogen chemistry and its deuterated isotopologues. With this new network, the utility of the [HCO⁺]/[DCO⁺] abundance ratio as a probe of the cosmic-ray ionization rate has been re-examined, with special attention paid to the effect of the initial value of the ortho-to-para ratio (OPR) of molecular hydrogen. After discussing the use of the probe for cold cores, we compare our results with previous theoretical and observational results for a molecular cloud close to the supernova remnant W51C, which is thought to have an enhanced cosmic-ray ionization rate ζ caused by the nearby γ-ray source. In addition, we attempt to use our approach to estimate the cosmic-ray ionization rate for L1174, a dense core with an embedded star. Beyond the previously known sensitivity of [HCO⁺]/[DCO⁺] to ζ, we demonstrate its additional dependence on the initial OPR and, secondarily, on the age of the source, its temperature, and its density. We conclude that the usefulness of the [HCO⁺]/[DCO⁺] abundance ratio in constraining the cosmic-ray ionization rate in dense regions increases with the age of the source and the ionization rate as the ratio becomes far less sensitive to the initial value of the OPR.