Cosmic-rays are the primary initiators of interstellar chemistry, and getting a better understanding of the varying impact they have on the chemistry of interstellar clouds throughout the Milky Way will not only expand our understanding of interstellar medium chemistry in our own galaxy, but also aid in extra-galactic studies. This work uses the ALCHEMIC astrochemical modeling code to perform numerical simulations of chemistry for a range of ionization rates. We study the impact of variations in the cosmic ray ionization rate on molecular abundances under idealized conditions, given by constant temperatures and a fixed density of 1e4 cm-3. As part of this study we examine whether observations of molecular abundances can be used to infer the cosmic ray ionization rate in such a simplified case. We find that intense cosmic-ray ionisation results in molecules, in particular the large and complex ones, being largely dissociated, and the medium becoming increasingly atomic. Individual species have limitations in their use as probes of the cosmic ray ionization rate. At early time (<1 Myrs) ions such as N2H+ and HOC+ make the best probes, while at later times, neutral species such as HNCO and SO stand out, in particular due to their large abundance variations. It is however by combining species into pairs that we find the best probes. Molecular ions such as N2H+ combined with different neutral species can provide probe candidates that outmatch individual species, in particular N2H+/C4H, N2H+/C2H, HOC+/O and HOC+/HNCO. These still have limitations to their functional range, but are more functional as probes than individual species previously used.

Aims: Our aim is to determine the critical parameters in water chemistry and the contribution of water to the oxygen budget by observing and modelling water gas and ice for a sample of eleven low-mass protostars, for which both forms of water have been observed. Methods: A simplified chemistry network, which is benchmarked against more sophisticated chemical networks, is developed that includes the necessary ingredients to determine the water vapour and ice abundance profiles in the cold, outer envelope in which the temperature increases towards the protostar. Comparing the results from this chemical network to observations of water emission lines and previously published water ice column densities, allows us to probe the influence of various agents (e.g., FUV field, initial abundances, timescales, and kinematics). Results: The observed water ice abundances with respect to hydrogen nuclei in our sample are 30-80ppm, and therefore contain only 10-30% of the volatile oxygen budget of 320ppm. The keys to reproduce this result are a low initial water ice abundance after the pre-collapse phase together with the fact that atomic oxygen cannot freeze-out and form water ice in regions with T(dust)>15 K. This requires short prestellar core lifetimes of less than about 0.1Myr. The water vapour profile is shaped through the interplay of FUV photodesorption, photodissociation, and freeze-out. The water vapour line profiles are an invaluable tracer for the FUV photon flux and envelope kinematics. Conclusions: The finding that only a fraction of the oxygen budget is locked in water ice can be explained either by a short pre-collapse time of less than 0.1 Myr at densities of n(H)~1e4 cm-3, or by some other process that resets the initial water ice abundance for the post-collapse phase. A key for the understanding of the water ice abundance is the binding energy of atomic oxygen on ice.

The chemistry in the diffuse interstellar medium initiates the gradual increase of molecular complexity during the life cycle of matter. A key molecule that enables build-up of new molecular bonds and new molecules via proton-donation is H3+. Its evolution is tightly related to molecular hydrogen and thought to be well understood. However, recent observations of ortho and para lines of H2 and H3+ in the diffuse ISM showed a puzzling discrepancy in nuclear spin excitation temperatures and populations between these two key species. H3+, unlike H2, seems to be out of thermal equilibrium, contrary to the predictions of modern astrochemical models. We conduct the first time-dependent modeling of the para-fractions of H2 and H3+ in the diffuse ISM and compare our results to a set of line-of-sight observations, including new measurements presented in this study. We isolate a set of key reactions for H3+ and find that the destruction of the lowest rotational states of H3+ by dissociative recombination largely control its ortho/para ratio. A plausible agreement with observations cannot be achieved unless a ratio larger than 1:5 for the destruction of (1,1)- and (1,0)-states of H3+ is assumed. Additionally, an increased CR ionization rate to 10(-15) 1/s further improves the fit whereas variations of other individual physical parameters, such as density and chemical age, have only a minor effect on the predicted ortho/para ratios. Thus our study calls for new laboratory measurements of the dissociative recombination rate and branching ratio of the key ion H3+ under interstellar conditions.

Formation and evolution of water in the Solar System and the origin of water on Earth constitute one of the most interesting questions in astronomy. The prevailing hypothesis for the origin of water on Earth is by delivery through water-rich small Solar system bodies. In this paper, the isotopic and chemical evolution of water during the early history of the solar nebula, before the onset of planetesimal formation, is studied. A gas-grain chemical model that includes multiply-deuterated species and nuclear spin-states is combined with a steady-state solar nebula model. To calculate initial abundances, we simulated 1 Myr of evolution of a cold and dark TMC1-like prestellar core. Two time-dependent chemical models of the solar nebula are calculated over 1 Myr: (1) a laminar model and (2) a model with 2D turbulent mixing. We find that the radial outward increase of the H2O D/H ratio is shallower in the chemo-dynamical nebular model compared to the laminar model. This is related to more efficient de-fractionation of HDO via rapid gas-phase processes, as the 2D mixing model allows the water ice to be transported either inward and thermally evaporated or upward and photodesorbed. The laminar model shows the Earth water D/H ratio at r ~<2.5 AU, while for the 2D chemo-dynamical model this zone is larger, r ~<9 AU. Similarly, the water D/H ratios representative of the Oort-family comets, ~2.5-10 x 10-4, are achieved within ~2-6 AU and ~2-20 AU in the laminar and the 2D model, respectively. We find that with regards to the water isotopic composition and the origin of the comets, the mixing model seems to be favored over the laminar model.

Observations of deuterated species are useful in probing the temperature, ionization level, evolutionary stage, chemistry, and thermal history of astrophysical environments. The analysis of data from ALMA and other new telescopes requires an elaborate model of deuterium fractionation. This paper presents a publicly available chemical network with multi-deuterated species and an extended, up-to-date set of gas-phase and surface reactions. To test this network, we simulate deuterium fractionation in diverse interstellar sources. Two cases of initial abundances are considered: i) atomic except for H2 and HD, and ii) molecular from a prestellar core. We reproduce the observed D/H ratios of many deuterated molecules, and sort the species according to their sensitivity to temperature gradients and initial abundances. We find that many multiply-deuterated species produced at 10 K retain enhanced D/H ratios at temperatures $\la 100$ K. We study how recent updates to reaction rates affect calculated D/H ratios, and perform a detailed sensitivity analysis of the uncertainties of the gas-phase reaction rates in the network. We find that uncertainties are generally lower in dark cloud environments than in warm IRDCs and that uncertainties increase with the size of the molecule and number of D-atoms. A set of the most problematic reactions is presented. We list potentially observable deuterated species predicted to be abundant in low- and high-mass star-formation regions.