We investigate deuterium chemistry coupled with the nuclear spin-state chemistry of H$_2$ and H$_3^+$ in protoplanetary disks. Multiple paths of deuterium fractionation are found; exchange reactions with D atoms, such as HCO$^+$ + D, are effective in addition to those with HD. In a disk model with grain sizes appropriate for dark clouds, the freeze-out of molecules is severe in the outer midplane, while the disk surface is shielded from UV radiation. Gaseous molecules, including DCO$^+$, thus become abundant at the disk surface, which tends to make their column density distribution relatively flat. If the dust grains have grown to millimeter size, the freeze-out rate of neutral species is reduced, and the abundances of gaseous molecules, including DCO$^+$ and N$_2$D$^+$, are enhanced in the cold midplane. Turbulent diffusion transports D atoms and radicals at the disk surface to the midplane, and stable ice species in the midplane to the disk surface. The effects of turbulence on chemistry are thus multifold; while DCO$^+$ and N$_2$D$^+$ abundances increase or decrease depending on the regions, HCN and DCN in the gas and ice are much reduced at the innermost radii, compared with the model without turbulence. When cosmic rays penetrate the disk, the ortho-to-para ratio (OPR) of H$_2$ is found to be thermal in the disk, except in the cold ($\lesssim 10$ K) midplane. We also analyze the OPR of H$_3^+$ and H$_2$D$^+$, as well as the main reactions of H$_2$D$^+$, DCO$^+$, and N$_2$D$^+$ to analytically derive their abundances in the cold midplane.

We investigate the delivery of regular and deuterated forms of water from prestellar cores to circumstellar disks. We adopt a semi-analytical axisymmetric two-dimensional collapsing core model with post-processing gas-ice astrochemical simulations, in which a layered ice structure is considered. The physical and chemical evolutions are followed until the end of the main accretion phase. When mass averaged over the whole disk, a forming disk has a similar H2O abundance and HDO/H2O abundance ratio as their precollapse values (within a factor of 2), regardless of time in our models. Consistent with previous studies, our models suggest that interstellar water ice is delivered to forming disks without significant alteration. On the other hand, the local vertically averaged H2O ice abundance and HDO/H2O ice ratio can differ more, by up to a factor of several, depending on time and distance from a central star. Key parameters for the local variations are the fluence of stellar UV photons en route into the disk and the ice layered structure, the latter of which is mostly established in the prestellar stages. We also find that even if interstellar water ice is destroyed by stellar UV and (partly) reformed prior to disk entry, the HDO/H2O ratio in reformed water ice is similar to the original value. This finding indicates that some caution is needed in discussions on the prestellar inheritance of H2O based on comparisons between the observationally derived HDO/H2O ratio in clouds/cores and that in disks/comets. Alternatively, we propose that the ratio of D2O/HDO to HDO/H2O better probes the prestellar inheritance of H2O. It is also found that icy organics are more enriched in deuterium than water ice in forming disks. The differential deuterium fractionation in water and organics is inherited from the prestellar stages.

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 reexamined, with special attention paid to the effect of the initial value of the molecular hydrogen ortho-to-para ratio (OPR). After discussing the use of the probe for cold cores, we then 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 $\zeta$ caused by the nearby $\gamma$-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 $\zeta$, 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 to constrain the cosmic ray ionization rate in dense regions increases with source age and ionization rate as the ratio becomes far less sensitive to the initial value of the OPR.

The first hydrostatic core, also called the first Larson core, is one of the first steps in low-mass star formation, as predicted by theory. With recent and future high performance telescopes, details of these first phases become accessible, and observations may confirm theory and even bring new challenges for theoreticians. In this context, we study from a theoretical point of view the chemical and physical evolution of the collapse of prestellar cores until the formation of the first Larson core, in order to better characterize this early phase in the star formation process. We couple a state-of-the-art hydrodynamical model with full gas-grain chemistry, using different assumptions on the magnetic field strength and orientation. We extract the different components of each collapsing core (i.e., the central core, the outflow, the disk, the pseudodisk, and the envelope) to highlight their specific physical and chemical characteristics. Each component often presents a specific physical history, as well as a specific chemical evolution. From some species, the components can clearly be differentiated. The different core models can also be chemically differentiated. Our simulation suggests some chemical species as tracers of the different components of a collapsing prestellar dense core, and as tracers of the magnetic field characteristics of the core. From this result, we pinpoint promising key chemical species to be observed.

We investigate the water deuteration ratio and ortho-to-para nuclear spin ratio of H2 (OPR(H2)) during the formation and early evolution of a molecular cloud, following the scenario that accretion flows sweep and accumulate HI gas to form molecular clouds. We follow the physical evolution of post-shock materials using a one-dimensional shock model, with post-processing gas-ice chemistry simulations. This approach allows us to study the evolution of the OPR(H2) and water deuteration ratio without an arbitrary assumption concerning the initial molecular abundances, including the initial OPR(H2). When the conversion of hydrogen into H2 is almost complete, the OPR(H2) is already much smaller than the statistical value of three due to the spin conversion in the gas phase. As the gas accumulates, the OPR(H2) decreases in a non-equilibrium manner. We find that water ice can be deuterium-poor at the end of its main formation stage in the cloud, compared to water vapor observed in the vicinity of low-mass protostars where water ice is likely sublimated. If this is the case, the enrichment of deuterium in water should mostly occur at somewhat later evolutionary stages of star formation, i.e., cold prestellar/protostellar cores. The main mechanism to suppress water ice deuteration in the cloud is the cycle of photodissociation and reformation of water ice, which efficiently removes deuterium from water ice chemistry. The removal efficiency depends on the main formation pathway of water ice. The OPR(H2) plays a minor role in water ice deuteration at the main formation stage of water ice.

Water is usually the main component of ice mantles, which cover the cores of dust grains in cold portions of dense interstellar clouds. When molecular hydrogen is adsorbed onto an icy mantle through physisorption, a common assumption in gas-grain rate equation models is to use an adsorption energy for molecular hydrogen on a pure water substrate. However, at high density and low temperature, when H2 is efficiently adsorbed onto the mantle, its surface abundance can be strongly overestimated if this assumption is still used. Unfortunately, the more detailed microscopic Monte Carlo treatment cannot be used to study the abundance of H2 in ice mantles if a full gas-grain network is utilized. We present a numerical method adapted for rate-equation models that takes into account the possibility that an H2 molecule can, while diffusing on the surface, find itself bound to another hydrogen molecule, with a far weaker bond than the H2-water bond, which can lead to more efficient desorption. We label the ensuing desorption "encounter desorption". The method is implemented first in a simple system consisting only of hydrogen molecules at steady state between gas and dust using the rate-equation approach and comparing the results with the results of a microscopic Monte Carlo calculation. We then discuss the use of the rate-equation approach with encounter desorption embedded in a complete gas-grain chemical network. For both systems, the rate-equation model with encounter desorption reproduces the H2 granular coverage computed by the microscopic Monte Carlo model. The method is especially useful for dense and cold environments, and for time-dependent physical conditions, such as occur in the collapse of dense cores and the formation of protoplanetary disks. It is not significantly CPU time consuming, so can be used for example with complex 3D chemical-hydrodynamical simulations.

Understanding water deuterium fractionation is important for constraining the mechanisms of water formation in interstellar clouds. Observations of HDO and H$_2^{18}$O transitions were carried out towards the high-mass star-forming region G34.26+0.15 with the HIFI instrument onboard the Herschel Space Observatory, as well as with ground-based single-dish telescopes. Ten HDO lines and three H$_2^{18}$O lines covering a broad range of upper energy levels (22-204 K) were detected. We used a non-LTE 1D analysis to determine the HDO/H$_2$O ratio as a function of radius in the envelope. Models with different water abundance distributions were considered in order to reproduce the observed line profiles. The HDO/H$_2$O ratio is found to be lower in the hot core ($\sim$3.5 $\times$ 10$^{-4}$ - 7.5 $\times$ 10$^{-4}$) than in the colder envelope ($\sim$1.0 $\times$ 10$^{-3}$ - 2.2 $\times$ 10$^{-3}$). This is the first time that a radial variation of the HDO/H$_2$O ratio has been found to occur in a high-mass source. The chemical evolution of this source was modeled as a function of its radius and the observations are relatively well reproduced. The comparison between the chemical model and the observations leads to an age of $\sim$10$^5$ years after the infrared dark cloud stage.

An outstanding question of astrobiology is the link between the chemical composition of planets, comets, and other solar system bodies and the molecules formed in the interstellar medium. Understanding the chemical and physical evolution of the matter leading to the formation of protoplanetary disks is an important step for this. We provide some new clues to this long-standing problem using three-dimensional chemical simulations of the early phases of disk formation: we interfaced the full gas-grain chemical model Nautilus with the radiation-magnetohydrodynamic model RAMSES, for different configurations and intensities of magnetic field. Our results show that the chemical content (gas and ices) is globally conserved during the collapsing process, from the parent molecular cloud to the young disk surrounding the first Larson core. A qualitative comparison with cometary composition suggests that comets are constituted of different phases, some molecules being direct tracers of interstellar chemistry, while others, including complex molecules, seem to have been formed in disks, where higher densities and temperatures allow for an active grain surface chemistry. The latter phase, and its connection with the formation of the first Larson core, remains to be modeled.

Many chemical models of dense interstellar clouds predict that the majority of gas-phase elemental nitrogen should be present as N2, with an abundance approximately five orders of magnitude less than that of hydrogen. As a homonuclear diatomic molecule, N2 is difficult to detect spectroscopically through infrared or millimetre-wavelength transitions so its abundance is often inferred indirectly through its reaction product N2H+. Two main formation mechanisms each involving two radical-radical reactions are the source of N2 in such environments. Here we report measurements of the low temperature rate constants for one of these processes, the N + CN reaction down to 56 K. The effect of the measured rate constants for this reaction and those recently determined for two other reactions implicated in N2 formation are tested using a gas-grain model employing a critically evaluated chemical network. We show that the amount of interstellar nitrogen present as N2 depends on the competition between its gas-phase formation and the depletion of atomic nitrogen onto grains. As the reactions controlling N2 formation are inefficient, we argue that N2 does not represent the main reservoir species for interstellar nitrogen. Instead, elevated abundances of more labile forms of nitrogen such as NH3 should be present on interstellar ices, promoting the eventual formation of nitrogen-bearing organic molecules.

Context: Dark cloud chemical models usually predict large amounts of O2, often above observational limits. Aims: We investigate the reason for this discrepancy from a theoretical point of view, inspired by the studies of Jenkins and Whittet on oxygen depletion. Methods: We use the gas-grain code Nautilus with an up-to-date gas-phase network to study the sensitivity of the molecular oxygen abundance to the oxygen elemental abundance. We use the rate coefficient for the reaction O + OH at 10 K recommended by the KIDA (KInetic Database for Astrochemistry) experts. Results: The updates of rate coefficients and branching ratios of the reactions of our gas-phase chemical network, especially N + CN and H3+ + O, have changed the model sensitivity to the oxygen elemental abundance. In addition, the gas-phase abundances calculated with our gas-grain model are less sensitive to the elemental C/O ratio than those computed with a pure gas-phase model. The grain surface chemistry plays the role of a buffer absorbing most of the extra carbon. Finally, to reproduce the low abundance of molecular oxygen observed in dark clouds at all times, we need an oxygen elemental abundance smaller than 1.6E-4. Conclusion: The chemistry of molecular oxygen in dense clouds is quite sensitive to model parameters that are not necessarily well constrained. That O2 abundance may be sensitive to nitrogen chemistry is an indication of the complexity of interstellar chemistry.