We present a new gas-grain chemical model to constrain the effect of grain size distribution on molecular abundances in starless and pre-stellar cores. We introduce grain-size dependence simultaneously for cosmic-ray (CR)-induced desorption efficiency and for grain equilibrium temperatures. We keep explicit track of ice abundances on a set of grain populations. We find that the size-dependent CR desorption efficiency affects ice abundances in a highly non-trivial way that depends on the molecule. Species that originate in the gas phase follow a simple pattern where the ice abundance is highest on the smallest grains (the most abundant in the distribution). Some molecules, such as HCN, are instead concentrated on large grains throughout the time evolution, while others (like $\rm N_2$) are initially concentrated on large grains, but at late times on small grains, due to grain-size-dependent competition between desorption and hydrogenation. Most of the water ice is on small grains at high medium density ($n({\rm H_2}) \gtrsim 10^6 \, \rm cm^{-3}$), where the water ice fraction, with respect to total water ice reservoir, can be as low as $\sim 10^{-3}$ on large (> 0.1 $\mu$m) grains. Allowing the grain equilibrium temperature to vary with grain size induces strong variations in relative ice abundances in low-density conditions where the interstellar radiation field and in particular its ultraviolet component are not attenuated. Our study implies consequences not only for the initial formation of ices preceding the starless core stage, but also for the relative ice abundances on the grain populations going into the protostellar stage. In particular, if the smallest grains can lose their mantles due to grain-grain collisions as the core is collapsing, the ice composition in the beginning of the protostellar stage could be very different to that in the pre-collapse phase.

C-fractionation has been studied from a theoretical point of view with different models of time-dependent chemistry, including both isotope-selective photodissociation and low-temperature isotopic exchange reactions. Recent chemical models predict that the latter may lead to a depletion of $^{13}$C in nitrile-bearing species, with $^{12}$C/$^{13}$C ratios two times higher than the elemental abundance ratio of 68 in the local ISM. Since the carbon isotopic ratio is commonly used to evaluate the $^{14}$N/$^{15}$N ratios with the double-isotope method, it is important to study C-fractionation in detail to avoid incorrect assumptions. In this work we implemented a gas-grain chemical model with new isotopic exchange reactions and investigated their introduction in the context of dense and cold molecular gas. In particular, we investigated the $^{12}$C/$^{13}$C ratios of HNC, HCN, and CN using a grid of models, with temperatures and densities ranging from 10 to 50 K and 2$\times$10$^{3}$ to 2$\times$10$^{7}$ cm$^{-3}$, respectively. We suggest a possible $^{13}$C exchange through the $^{13}$C + C$_{3}$ $\rightarrow$ $^{12}$C +$^{13}$CC$_{2}$ reaction, which does not result in dilution, but rather in $^{13}$C enhancement, for molecules formed starting from atomic carbon. This effect is efficient in a range of time between the formation of CO and its freeze-out on grains. Furthermore, we show that the $^{12}$C/$^{13}$C ratios of nitriles are predicted to be a factor 0.8-1.9 different from the local value of 68 for massive star-forming regions. This result also affects the $^{14}$N/$^{15}$N ratio: a value of 330 obtained with the double-isotope method is predicted to be 260-1150, depending on the physical conditions. Finally, we studied the $^{12}$C/$^{13}$C ratios by varying the cosmic-ray ionization rate: the ratios increase with it because of secondary photons and cosmic-ray reactions.

Deuterated interstellar molecules frequently have abundances relative to their main isotopologues much higher than the overall elemental D-to-H ratio in the cold dense interstellar medium. The H$_3^+$ and its isotopologues play a key role in the deuterium fractionation; however, the abundances of these isotopologues have not been measured empirically with respect to H$_3^+$ to date. Our aim was to constrain the relative abundances of H$_2$D$^+$ and D$_3^+$ in the cold outer envelope of the hot core CRL 2136 IRS 1. We carried out three observations targeting H$_3^+$ and its isotopologues using the spectrographs CRIRES at the VLT, iSHELL at IRTF, and EXES on board SOFIA. In addition, the CO overtone band at 2.3 $\mu$m was observed by iSHELL to characterize the gas on the line of sight. The H$_3^+$ ion was detected toward CRL 2136 IRS 1 as in previous observations. Spectroscopy of lines of H$_2$D$^+$ and D$_3^+$ resulted in non-detections. The 3$\sigma$ upper limits of $N({\rm H_2D^+})/N({\rm H_3^+})$ and $N({\rm D_3^+})/N({\rm H_3^+})$ are 0.24 and 0.13, respectively. The population diagram for CO is reproduced by two components of warm gas with the temperatures 58 K and 530 K, assuming a local thermodynamic equilibrium (LTE) distribution of the rotational levels. Cold gas ($<$20 K) makes only a minor contribution to the CO molecular column toward CRL 2136 IRS 1. The critical conditions for deuterium fractionation in a dense cloud are low temperature and CO depletion. Given the revised cloud properties, it is no surprise that H$_3^+$ isotopologues are not detected toward CRL 2136 IRS 1. The result is consistent with our current understanding of how deuterium fractionation proceeds.

We constructed two new models for deuterium and spin-state chemistry for the purpose of modeling the low-temperature environment prevailing in starless and pre-stellar cores. The fundamental difference between the two models is in the treatment of ion-molecule proton-donation reactions of the form $\rm XH^+ + Y \longrightarrow X + YH^+$, which are allowed to proceed either via full scrambling or via direct proton hop, i.e., disregarding proton exchange. The choice of the reaction mechanism affects both deuterium and spin-state chemistry, and in this work our main interest is on the effect on deuterated ammonia. We applied the new models to the starless core H-MM1, where several deuterated forms of ammonia have been observed. Our investigation slightly favors the proton hop mechanism over full scrambling because the ammonia D/H ratios are better fit by the former model, although neither model can reproduce the observed $\rm NH_2D$ ortho-to-para ratio of 3 (the models predict a value of $\sim$2). Extending the proton hop scenario to hydrogen atom abstraction reactions yields a good agreement for the spin-state abundance ratios, but greatly overestimates the deuterium fractions of ammonia. However, one can find a reasonably good agreement with the observations with this model by increasing the cosmic-ray ionization rate over the commonly-adopted value of $\sim$$10^{-17}\,\rm s^{-1}$. We also find that the deuterium fractions of several other species, such as $\rm H_2CO$, $\rm H_2O$, and $\rm CH_3$, are sensitive to the adopted proton-donation reaction mechanism. Whether the full scrambling or proton hop mechanism dominates may be dependent on the reacting system, and new laboratory and theoretical studies for various reacting systems are needed to constrain chemical models.

We develop a self-consistent model for the equilibrium gas temperature and size-dependent dust temperature in cold, dense pre-stellar cores, assuming an arbitrary power-law size distribution of dust grains. Compact analytical expressions applicable to a broad range of physical parameters are derived and compared with predictions of the commonly used standard model. It is suggested that combining the theoretical results with observations should allow us to constrain the degree of dust evolution and the cosmic-ray ionization rate in dense cores, and to help in discriminating between different regimes of cosmic-ray transport in molecular clouds. In particular, assuming a canonical MRN distribution of grain sizes, our theory demonstrates that the gas temperature measurements in the pre-stellar core L1544 are consistent with an ionization rate as high as $\sim 10^{-16}$ s$^{-1}$, an order of magnitude higher than previously thought.

The chemical structure of high-mass star nurseries is important for a general understanding of star formation. Deuteration is a key chemical process in the earliest stages of star formation because its efficiency is sensitive to the environment. Using the IRAM-30 m telescope at 1.3--4.3 mm wavelengths, we have imaged two parsec-scale high-mass protostellar clumps (P1 and S) that show different evolutionary stages but are located in the same giant filamentary {infrared dark cloud} G28.34+0.06. Deep spectral images at subparsec resolution reveal the dust and gas physical structures of both clumps. We find that (1) the low-$J$ lines of $\rm N_2H^+$, HCN, HNC, and $\rm HCO^+$ isotopologues are subthermally excited; and (2) the deuteration of $\rm N_2H^+$ is more efficient than that of $\rm HCO^+$, HCN, and HNC by an order of magnitude. The deuterations of these species are enriched toward the chemically younger clump S compared with P1, indicating that this process favors the colder and denser environment ($\rm T_{kin} \sim14 K$, $\rm N(NH_3) \sim 9\times 10^{15}\,cm^{-2}$). In contrast, single deuteration of $\rm NH_3$ is insensitive to the environmental difference between P1 and S; and (3) single deuteration of $\rm CH_3OH$ ($\rm > 10\%$) is detected toward the location where CO shows a depletion of $\sim10$. This comparative chemical study between P1 and S links the chemical variations to the environmental differences and shows chemical similarities between the early phases of high- and low-mass star-forming regions.

Context. The deuterium fraction in low-mass prestellar cores is a good diagnostic indicator of the initial phases of star formations, and it is also a fundamental quantity to infer the ionisation degree in these objects. Aims. With the analysis of multiple transitions of $\rm N_2H^+$, $\rm N_2D^+$, $\rm HC^{18}O^+$ and $\rm DCO^+$ we are able to determine the molecular column density maps and the deuterium fraction in $\rm N_2H^+$ and $\rm HCO^+$ toward the prototypical prestellar core L1544. This is the preliminary step to derive the ionisation degree in the source. Methods. We use a non-local thermodynamic equilibrium (non-LTE) radiative transfer code, combined with the molecular abundances derived from a chemical model, to infer the excitation conditions of all the observed transitions, which allows us to derive reliable maps of each molecule's column density. The ratio between the column density of a deuterated species and its non-deuterated counterpart gives the searched deuteration level. Results. The non-LTE analysis confirms that, for the analysed molecules, higher-J transitions are characterised by excitation temperatures $\approx 1-2\,$K lower than the lower-J ones. The chemical model that provides the best fit to the observational data predict the depletion of $\rm N_2H^+$ and to a lesser extent of $\rm N_2D^+$ in the innermost region. The peak values for the deuterium fraction that we find are $\mathrm{D/H_{N_2H^+}} = 0.26^{+0.15}_{-0.14}$ and $\mathrm{D/H_{HCO^+}} = 0.035^{+0.015}_{-0.012}$, in good agreement with previous estimates in the source.

At the low temperatures ($\sim$10 K) and high densities ($\sim$100,000 H$_2$ molecules per cc) of molecular cloud cores and protostellar envelopes, a large amount of molecular species (in particular those containing C and O) freeze-out onto dust grain surfaces. It is in these regions that the deuteration of H$_3^+$ becomes very efficient, with a sharp abundance increase of H$_2$D$^+$ and D$_2$H$^+$. The multi-deuterated forms of H$_3^+$ participate in an active chemistry: (i) their collision with neutral species produces deuterated molecules such as the commonly observed N$_2$D$^+$, DCO$^+$ and multi-deuterated NH$_3$; (ii) their dissociative electronic recombination increases the D/H atomic ratio by several orders of magnitude above the D cosmic abundance, thus allowing deuteration of molecules (e.g. CH$_3$OH and H$_2$O) on the surface of dust grains. Deuterated molecules are the main diagnostic tools of dense and cold interstellar clouds, where the first steps toward star and protoplanetary disk formation take place. Recent observations of deuterated molecules are reviewed and discussed in view of astrochemical models inclusive of spin-state chemistry. We present a new comparison between models based on complete scrambling (to calculate branching ratio tables for reactions between chemical species that include protons and/or deuterons) and models based on non-scrambling (proton hop) methods, showing that the latter best agree with observations of NH$_3$ deuterated isotopologues and their different nuclear spin symmetry states.

We carried out a parameter-space exploration of the ammonia abundance in the pre-stellar core L1544, where it has been observed to increase toward the center of the core with no signs of freeze-out onto grain surfaces. We considered static and dynamical physical models coupled with elaborate chemical and radiative transfer calculations, and explored the effects of varying model parameters on the (ortho+para) ammonia abundance profile. None of our models are able to reproduce the inward-increasing tendency in the observed profile; ammonia depletion always occurs in the center of the core. In particular, our study shows that including the chemical desorption process, where exothermic association reactions on the grain surface can result in the immediate desorption of the product molecule, leads to ammonia abundances that are over an order of magnitude above the observed level in the innermost 15000 au of the core - at least when one employs a constant efficiency for the chemical desorption process irrespective of the ice composition. Our results seemingly constrain the chemical desorption efficiency of ammonia on water ice to below 1%. It is increasingly evident that time-dependent effects must be considered so that the results of chemical models can be reconciled with observations.

We use the deuteration of $c$-$\mathrm{C_{3}H_{2}}$ to probe the physical parameters of starless and protostellar cores, related to their evolutionary states, and compare it to the $\mathrm{N_{2}H^{+}}$-deuteration in order to study possible differences between the deuteration of C- and N-bearing species. We observed the main species $c$-$\mathrm{C_{3}H_{2}}$, the singly and doubly deuterated species $c$-$\mathrm{C_{3}HD}$ and $c$-$\mathrm{C_{3}D_{2}}$, as well as the isotopologue $c$-$\mathrm{{H^{13}CC_{2}H}}$ toward 10 starless cores and 5 protostars in the Taurus and Perseus Complexes. We examined the correlation between the $N$($c$-$\mathrm{C_{3}HD}$)/$N$($c$-$\mathrm{C_{3}H_{2}}$) ratio and the dust temperature along with the $\mathrm{H_2}$ column density and the CO depletion factor. The resulting $N$($c$-$\mathrm{C_{3}HD}$)/$N$($c$-$\mathrm{C_{3}H_{2}}$) ratio is within the error bars consistent with $10\%$ in all starless cores with detected $c$-$\mathrm{C_{3}HD}$. This also accounts for the protostars except for the source HH211, where we measure a high deuteration level of $23\%$. The deuteration of $\mathrm{N_{2}H^{+}}$ follows the same trend but is considerably higher in the dynamically evolved core L1544. Toward the protostellar cores the coolest objects show the largest deuterium fraction in $c$-$\mathrm{C_{3}H_{2}}$. We show that the deuteration of $c$-$\mathrm{C_{3}H_{2}}$ can trace the early phases of star formation and is comparable to that of $\mathrm{N_{2}H^{+}}$. However, the largest $c$-$\mathrm{C_{3}H_{2}}$ deuteration level is found toward protostellar cores, suggesting that while $c$-$\mathrm{C_{3}H_{2}}$ is mainly frozen onto dust grains in the central regions of starless cores, active deuteration is taking place on ice.

We quantify if the chemical abundance gradients given by a dynamical model of core collapse including time-dependent changes in density and temperature differ greatly from abundances derived from static models, where the density and temperature structures of the core are kept fixed as the chemistry evolves. For this study we developed a new one-dimensional spherically symmetric hydrodynamics code that couples the hydrodynamics equations with a comprehensive time-dependent gas-grain chemical model, including deuterium and spin-state chemistry, and radiative transfer calculations to derive self-consistent time-dependent chemical abundance gradients. We applied the code to model the collapse of a starless core up to the point when the infall flow becomes supersonic. The abundances predicted by the dynamical and static models are almost identical during the quiescent phase of core evolution, but the results start to diverge after the onset of core collapse, where the static model underestimates abundances at high medium density (inner core) and underestimates them at low density (outer core), and this is clearly reflected in simulated lines. The static model generally overestimates deuteration, which is increasingly evident the more D atoms are substituted in the molecule. We also find that using a limited chemical network, or a limited set of cooling molecules, may lead to an overestimate of the collapse timescale, and in some cases may prevent the collapse altogether. In our model, most of the line cooling near the center of the core is due to HCN, CO, and NO. In conclusion, the use of a static physical model is not a reliable method of simulating chemical abundances in starless cores after the onset of gravitational collapse. The adoption of complex chemistry and a comprehensive set of cooling molecules is necessary to model the collapse adequately.

Aims. We study whether rotational excitation makes a difference to the abundances of the $\rm H_3^+$ isotopologs, including spin states, in physical conditions corresponding to starless cores and protostellar envelopes. Methods. We developed a new rate coefficient set for the $\rm H_3^+$ isotopologs, allowing for rotational excitation, using the state-to-state rate coefficients from Hugo et al. These new so-called species-to-species rate coefficients are compared with previously-used ground state-to-species rate coefficients. Results. The species-to-species and ground state-to-species model results differ at high density and toward increasing temperatures ($T > 10$ K). The species-to-species model predicts a lower $\rm H_3^+$ deuteration degree at high density owing to an increase of the rate coefficients of endothermic reactions that decrease deuteration. At 20 K the ground state-to-species model overestimates the abundance of $\rm H_2D^+$ by a factor of about two while the abundance of $\rm D_3^+$ can differ by an order of magnitude between the models. Spin-state abundance ratios are also affected, and the new model better reproduces recent observations of ortho and para $\rm H_2D^+$ and $\rm D_2H^+$. The applicability regime of the new rate coefficients depends on the critical densities of the various rotational transitions. Conclusions. The difference in the abundances of the $\rm H_3^+$ isotopologs predicted by the two models is negligible at 10 K but excited states are very important in studies of deuteration at higher temperatures, for example in protostellar envelopes. The species-to-species rate coefficients provide a more realistic approach to the chemistry of the $\rm H_3^+$ isotopologs than the ground state-to-species rate coefficients do, and so the former should be adopted in chemical models describing the chemistry of the $\rm H_3^+ + H_2$ reacting system.

We report on the detection of the ground-state rotational line of ortho-D2H+ at 1.477 THz (203 micron) using the German REceiver for Astronomy at Terahertz frequencies (GREAT) onboard the Stratospheric Observatory For Infrared Astronomy (SOFIA). The line is seen in absorption against far-infrared continuum from the protostellar binary IRAS 16293-2422 in Ophiuchus. The para-D2H+ line at 691.7 GHz was not detected with the APEX telescope toward this position. These D2H+ observations complement our previous detections of para-H2D+ and ortho-H2D+ using SOFIA and APEX. By modeling chemistry and radiative transfer in the dense core surrounding the protostars, we find that the ortho-D2H+ and para-H2D+ absorption features mainly originate in the cool (T<18 K) outer envelope of the core. In contrast, the ortho-H2D+ emission from the core is significantly absorbed by the ambient molecular cloud. Analyses of the combined D2H+ and H2D+ data result in an age estimate of ~500 000 yr for the core, with an uncertainty of ~200 000 yr. The core material has probably been pre-processed for another 500 000 years in conditions corresponding to those in the ambient molecular cloud. The inferred time scale is more than ten times the age of the embedded protobinary. The D2H+ and H2D+ ions have large and nearly equal total (ortho+para) fractional abundances of ~$10^{-9}$ in the outer envelope. This confirms the central role of H3+ in the deuterium chemistry in cool, dense gas, and adds support to the prediction of chemistry models that also D3+ should be abundant in these conditions.

HOCO$^+$ is a polar molecule that represents a useful proxy for its parent molecule CO$_2$, which is not directly observable in the cold interstellar medium. This cation has been detected towards several lines of sight, including massive star forming regions, protostars, and cold cores. Despite the obvious astrochemical relevance, protonated CO$_2$ and its deuterated variant, DOCO$^+$, still lack an accurate spectroscopic characterisation. The aim of this work is to extend the study of the ground-state pure rotational spectra of HOCO$^+$ and DOCO$^+$ well into the sub-millimetre region. Ground-state transitions have been recorded in the laboratory using a frequency-modulation absorption spectrometer equipped with a free-space glow-discharge cell. The ions were produced in a low-density, magnetically-confined plasma generated in a suitable gas mixture. The ground-state spectra of HOCO$^+$ and DOCO$^+$ have been investigated in the 213-967 GHz frequency range, with the detection of 94 new rotational transitions. Additionally, 46 line positions taken from the literature have been accurately remeasured. The newly-measured lines have significantly enlarged the available data sets for HOCO$^+$ and DOCO$^+$, thus enabling the determination of highly accurate rotational and centrifugal distortion parameters. Our analysis showed that all HOCO$^+$ lines with Ka $\geq$ 3 are perturbed by a ro-vibrational interaction that couples the ground state with the v$_ 5$ = 1 vibrationally-excited state. This resonance has been explicitly treated in the analysis in order to obtain molecular constants with clear physical meaning. The improved sets of spectroscopic parameters provide enhanced lists of very accurate, sub-millimetre rest-frequencies of HOCO$^+$ and DOCO$^+$ for astrophysical applications. These new data challenges a recent tentative identification of DOCO$^+$ toward a pre-stellar core.

Aims. We investigate the effect of chemistry on the stability of starless cores against gravitational collapse. Methods. We combine chemical and radiative transfer simulations in the context of a modified Bonnor-Ebert sphere to model the effect of chemistry on the gas temperature, and study the effect of temperature changes on core stability. Results. We find that chemistry has in general very little effect on the nondimensional radius $\xi_{\rm out}$ which parametrizes the core stability. Cores that are initially stable or unstable tend to stay near their initial states, in terms of stability (i.e., $\xi_{\rm out} \sim$ constant), as the chemistry develops. This result is independent of the initial conditions. We can however find solutions where $\xi_{\rm out}$ decreases at late times ($t \gtrsim 10^6 \, \rm yr$) which correspond to increased stabilization caused by the chemistry. Even though the core stability is unchanged by the chemistry in most of the models considered here, we cannot rule out the possibility that a core can evolve from an unstable to a stable state owing to chemical evolution. The reverse case, where an initially stable core becomes ultimately unstable, seems highly unlikely. Conclusions. Our results indicate that chemistry should be properly accounted for in studies of star-forming regions, and that further investigations of core stability especially with hydrodynamical models are warranted.

Deuterated species are unique and powerful tools in astronomy since they can probe the physical conditions, chemistry, and ionization level of various astrophysical media. Recent observations of several deuterated species along with some of their spin isomeric forms have rekindled the interest for more accurate studies on deuterium fractionation. This paper presents the first publicly available chemical network of multiply deuterated species along with spin chemistry implemented on the latest state-of-the-art gas-grain chemical code `NAUTILUS'. D/H ratios for all deuterated species observed at different positions of TMC-1 are compared with the results of our model, which considers multiply deuterated species along with the spin chemistry of light hydrogen bearing species H2, H2+, H3+ and their isotopologues. We also show the differences in the modeled abundances of non-deuterated species after the inclusion of deuteration and spin chemistry in the model. Finally, we present a list of potentially observable deuterated species in TMC-1 awaiting detection.

Ammonia and its deuterated isotopologues probe physical conditions in dense molecular cloud cores. With the aim of testing the current understanding of the spin-state chemistry of these molecules, we observed spectral lines of NH3, NH2D, NHD2, ND3, and N2D+ towards a dense, starless core in Ophiuchus with the APEX, GBT, and IRAM 30-m telescopes. The observations were interpreted using a gas-grain chemistry model combined with radiative transfer calculations. The chemistry model distinguishes between the different nuclear spin states of light hydrogen molecules, ammonia, and their deuterated forms. High deuterium fractionation ratios with NH2D/NH3=0.4, NHD2/NH2D=0.2, and ND3/NHD2=0.06 were found in the core. The observed ortho/para ratios of NH2D and NHD2 are close to the corresponding nuclear spin statistical weights. The chemistry model can approximately reproduce the observed abundances, but predicts uniformly too low ortho/para-NH2D, and too large ortho/para-NHD2 ratios. The longevity of N2H+ and NH3 in dense gas, which is prerequisite to their strong deuteration, can be attributed to the chemical inertia of N2 on grain surfaces. The discrepancies between the chemistry model and the observations are likely to be caused by the fact that the model assumes complete scrambling in principal gas-phase deuteration reactions of ammonia, which means that all the nuclei are mixed in reactive collisions. If, instead, these reactions occur through proton hop/hydrogen abstraction processes, statistical spin ratios are to be expected. The present results suggest that while the deuteration of ammonia changes with physical conditions and time, the nuclear spin ratios of ammonia isotopologues do not probe the evolutionary stage of a cloud.

Aims. We aim to understand the high cyclic-to-linear $\rm C_3H_2$ ratio ($32 \pm 4$) observed toward L1544 by Spezzano et al. (2016). Methods. We combine a gas-grain chemical model with a physical model for L1544 to simulate the column densities of cyclic and linear $\rm C_3H_2$ observed toward L1544. The most important reactions for the formation and destruction of both forms of $\rm C_3H_2$ are identified, and their relative rate coefficients are varied to find the best match to the observations. Results. We find that the ratio of the rate coefficients of $\rm C_3H_3^+ + e^- \longrightarrow C_3H_2 + H$ for cyclic and linear $\rm C_3H_2$ must be $\sim 20$ in order to reproduce the observations, depending on the branching ratios assumed for the $\rm C_3H_3^+ + e^- \longrightarrow C_3H + H_2$ reaction. In current astrochemical networks it is assumed that cyclic and linear $\rm C_3H_2$ are formed in a 1:1 ratio in the aforementioned reactions. Laboratory studies and/or theoretical calculations are needed to confirm the results of our chemical modeling based on observational constraints.

Aims. We aim to investigate whether a multilayer ice model can be as successful as a bulk ice model in reproducing the observed abundances of various deuterated gas-phase species toward starless cores. Methods. We calculate abundances for various deuterated species as functions of time adopting fixed physical conditions. We also estimate abundance gradients by adopting a modified Bonnor-Ebert sphere as a core model. In the multilayer ice scenario, we consider desorption from one or several monolayers on the surface. Results. We find that the multilayer model predicts abundances of $\rm DCO^+$ and $\rm N_2D^+$ that are about an order of magnitude lower than observed, caused by the trapping of CO and $\rm N_2$ into the grain mantle. As a result of the mantle trapping, deuteration efficiency in the gas phase increases and we find stronger deuterium fractionation in ammonia than what has been observed. Another distinguishing feature of the multilayer model is that $\rm D_3^+$ becomes the main deuterated ion at high density. The bulk ice model is generally easily reconciled with observations. Conclusions. Our results underline that more theoretical and experimental work is needed to understand the composition and morphology of interstellar ices, and the desorption processes that can act on them. With the current constraints, the bulk ice model appears to be better in reproducing observations than the multilayer ice model. According to our results, the $\rm H_2D^+$ to $\rm N_2D^+$ abundance ratio is higher than 100 in the multilayer model, while only a few $\times$ 10 in the bulk model, and so observations of this ratio could provide information on the ice morphology in starless cores. Observations of the abundance of $\rm D_3^+$ compared to $\rm H_2D^+$ and $\rm D_2H^+$ would provide additional constraints for the models.

The 1$_{11}$-1$_{01}$ lines of ortho and para--NH$_2$D (o/p-NH$_2$D), respectively at 86 and 110 GHz, are commonly observed to provide constraints on the deuterium fractionation in the interstellar medium. In cold regions, the hyperfine structure due to the nitrogen ($^{14}$N) nucleus is resolved. To date, this splitting is the only one which is taken into account in the NH$_2$D column density estimates. We investigate how the inclusion of the hyperfine splitting caused by the deuterium (D) nucleus affects the analysis of the rotational lines of NH$_2$D. We present 30m IRAM observations of the above mentioned lines, as well as APEX o/p-NH$_2$D observations of the 1$_{01}$-0$_{00}$ lines at 333 GHz. The hyperfine spectra are first analyzed with a line list that only includes the hyperfine splitting due to the $^{14}$N nucleus. We find inconsistencies between the line widths of the 1$_{01}$-0$_{00}$ and 1$_{11}$-1$_{01}$ lines, the latter being larger by a factor of $\sim$1.6$\pm0.3$. Such a large difference is unexpected given the two sets of lines are likely to originate from the same region. We next employ a newly computed line list for the o/p-NH$_2$D transitions, where the hyperfine structure induced by both nitrogen and deuterium nuclei is included. With this new line list, the analysis of the previous spectra leads to linewidths which are compatible. Neglecting the hyperfine structure owing to D leads to overestimate the linewidths of the o/p-NH$_2$D lines at 3 mm. The error for a cold molecular core is about 50\%. This error propagates directly to the column density estimate. It is therefore recommended to take into account the hyperfine splittings caused by both the $^{14}$N and D nuclei in any analysis relying on these lines.

Aims. We investigate the stability of non-isothermal Bonnor-Ebert spheres with a model that includes a self-consistent calculation of the gas temperature. This way we can discard the assumption of equality between the dust and gas temperatures, and study the stability as the gas temperature changes with chemical evolution of the gas. Methods. We use a gas-grain chemical model including a time-dependent treatment of depletion onto grain surfaces, which strongly influences the gas temperature as the main coolant, CO, depletes from the gas. Dust and gas temperatures are solved with radiative transfer. For comparison with previous work, we assume that the cores are deeply embedded in a larger external structure, corresponding to visual extinction $A_{\rm V}^{\rm ext}=10$ mag. Results. We find that the critical non-dimensional radius $\xi_1$ derived here is similar to our previous work where we assumed $T_{\rm dust}=T_{\rm gas}$; the $\xi_1$ values lie below the isothermal critical value $\xi_0\sim6.45$, but the difference is less than 10%. Chemical evolution does not affect notably the stability condition of low-mass cores (<0.75 $M_\odot$). For higher masses the decrease of cooling owing to CO depletion causes substantial temporal changes in the temperature and density profiles of the cores. In the mass range 1-2 $M_\odot$ , $\xi_1$ decreases with chemical evolution, whereas above 3 $M_\odot$ , $\xi_1$ instead increases. We also find that decreasing $A_{\rm V}^{\rm ext}$ increases the gas temperature especially when the gas is chemically old, causing $\xi_1$ to increase with respect to models with higher $A_{\rm V}^{\rm ext}$. The derived $\xi_1$ values are close to $\xi_0$. The density contrast between the core center and edge varies between 8 to 16 depending on core mass and the chemical age of the gas, compared to the constant value $\sim$ 14.1 for the isothermal BES.

Aims. We aim to develop a chemical model that contains a consistent description of spin-state chemistry in reactions involving chemical species with multiple deuterons. We apply the model to the specific case of deuterated ammonia, to derive values for the various spin-state ratios. Methods. We apply symmetry rules in the complete scrambling assumption to calculate branching ratio tables for reactions between chemical species that include multiple protons and/or deuterons. Reaction sets for both gas-phase and grain-surface chemistry are generated using an automated routine that forms all possible spin-state variants of any given reaction with up to six H/D atoms. Single-point and modified Bonnor-Ebert models are used to study the density and temperature dependence of ammonia and its isotopologs, and the associated spin-state ratios. Results. We find that the spin-state ratios of the ammonia isotopologs are, at late times, very different from their statistical values. The ratios are rather insensitive to variations in the density, but present strong temperature dependence. We derive high peak values ($\sim$ 0.1) for the deuterium fraction in ammonia, in agreement with previous (gas-phase) models. The deuterium fractionation is strongest at high density, corresponding to a high degree of depletion, and also presents temperature dependence. We find that in the temperature range 5 to 20 K, the deuterium fractionation peaks at $\sim$ 15 K while most of the ortho/para (and meta/para for $\rm ND_3$) ratios present a minimum at 10 K (ortho/para $\rm NH_2D$ has instead a maximum at this temperature). Conclusions. Owing to the density and temperature dependence found in the abundances and spin-state ratios of ammonia and its isotopologs, it is evident that observations of ammonia and its deuterated forms can provide important constraints on the physical structure of molecular clouds.

The deuterium fraction [N$_2$D$^+$]/[N$_2$H$^+$], may provide information about the ages of dense, cold gas structures, important to compare with dynamical models of cloud core formation and evolution. Here we introduce a complete chemical network with species containing up to three atoms, with the exception of the Oxygen chemistry, where reactions involving H$_3$O$^+$ and its deuterated forms have been added, significantly improving the consistency with comprehensive chemical networks. Deuterium chemistry and spin states of H$_2$ and H$_3^+$ isotopologues are included in this primarily gas-phase chemical model. We investigate dependence of deuterium chemistry on model parameters: density ($n_{\rm H}$), temperature, cosmic ray ionization rate, and gas-phase depletion factor of heavy elements ($f_{\rm D}$). We also explore the effects of time-dependent freeze-out of gas-phase species and dynamical evolution of density at various rates relative to free-fall collapse. For a broad range of model parameters, the timescales to reach large values of $D_{\rm frac}^{\rm N_2H^+} \gtrsim 0.1$, observed in some low- and high-mass starless cores, are relatively long compared to the local free-fall timescale. These conclusions are unaffected by introducing time-dependent freeze-out and considering models with evolving density, unless the initial $f_{\rm D} \gtrsim$ 10. For fiducial model parameters, achieving $D_{\rm frac}^{\rm N_2H^+} \gtrsim 0.1$ requires collapse to be proceeding at rates at least several times slower than that of free-fall collapse, perhaps indicating a dynamically important role for magnetic fields in the support of starless cores and thus the regulation of star formation.

Aims. We aim to present simulated chemical abundance profiles for a variety of important species, with special attention given to spin-state chemistry, in order to provide reference results against which present and future models can be compared. Methods. We employ gas-phase and gas-grain models to investigate chemical abundances in physical conditions corresponding to starless cores. To this end, we have developed new chemical reaction sets for both gas-phase and grain-surface chemistry, including the deuterated forms of species with up to six atoms and the spin-state chemistry of light ions and of the species involved in the ammonia and water formation networks. The physical model is kept simple in order to facilitate straightforward benchmarking of other models against the results of this paper. Results. We find that the ortho/para ratios of ammonia and water are similar in both gas-phase and gas-grain models, at late times in particular, implying that the ratios are determined by gas-phase processes. We derive late-time ortho/para ratios of ~0.5 and ~1.6 for ammonia and water, respectively. We find that including or excluding deuterium in the calculations has little effect on the abundances of non-deuterated species and on the ortho/para ratios of ammonia and water, especially in gas-phase models where deuteration is naturally hindered owing to the presence of abundant heavy elements. Although we study a rather narrow temperature range (10-20 K), we find strong temperature dependence in, e.g., deuteration and nitrogen chemistry. For example, the depletion timescale of ammonia is significantly reduced when the temperature is increased from 10 to 20 K; this is because the increase in temperature translates into increased accretion rates, while the very high binding energy of ammonia prevents it from being desorbed at 20 K.

Extremely large deuteration of several molecules has been observed towards prestellar cores and low-mass protostars for a decade. New observations performed towards low-mass protostars suggest that water presents a lower deuteration in the warm inner gas than in the cold external envelope. We coupled a gas-grain astrochemical model with a one-dimension model of collapsing core to properly follow the formation and the deuteration of interstellar ices as well as their subsequent evaporation in the low-mass protostellar envelopes with the aim of interpreting the spatial and temporal evolutions of their deuteration. The astrochemical model follows the formation and the evaporation of ices with a multilayer approach and also includes a state-of-the-art deuterated chemical network by taking the spin states of H$_2$ and light ions into account. Because of their slow formation, interstellar ices are chemically heterogeneous and show an increase of their deuterium fractionation towards the surface. The differentiation of the deuteration in ices induces an evolution of the deuteration within protostellar envelopes. The warm inner region is poorly deuterated because it includes the whole molecular content of ices while the deuteration predicted in the cold external envelope scales with the highly deuterated surface of ices. We are able to reproduce the observed evolution of water deuteration within protostellar envelopes but we are still unable to predict the super-high deuteration observed for formaldehyde and methanol. Finally, the extension of this study to the deuteration of complex organics (COMs), important for the prebiotic chemistry, shows a good agreement with the observations, suggesting that we can use the deuteration to retrace their mechanisms and their moments of formation.

Aims: We aim to investigate the abundances of light deuterium-bearing species such as HD, H2D+ and D2H+ in a gas-grain chemical model including an extensive description of deuterium and spin state chemistry, in physical conditions appropriate to the very centers of starless cores. Methods: We combine a gas-grain chemical model with radiative transfer calculations to simulate density and temperature structure in starless cores. The chemical model includes deuterated forms of species with up to 4 atoms and the spin states of the light species H2, H2+ and H3+ and their deuterated forms. Results: We find that HD eventually depletes from the gas phase because deuterium is efficiently incorporated to grain-surface HDO, resulting in inefficient HD production on grains. HD depletion has consequences not only on the abundances of e.g. H2D+ and D2H+, whose production depends on the abundance of HD, but also on the spin state abundance ratios of the various light species, when compared with the complete depletion model where heavy elements do not influence the chemistry. Conclusions: While the eventual HD depletion leads to the disappearance of light deuterium-bearing species from the gas phase in a relatively short timescale at high density, we find that at late stages of core evolution the abundances of H2D+ and D2H+ increase toward the core edge and the disributions become extended. The HD depletion timescale increases if less oxygen is initially present in the gas phase, owing to chemical interaction between the gas and the dust predecing the starless core phase. Our results are greatly affected if H2 is allowed to tunnel on grain surfaces, and therefore more experimental data not only on tunneling but also on the O + H2 surface reaction in particular is needed.

Aims: We aim to simulate radial profiles of molecular abundances and the gas temperature in cold and heavily shielded starless cores by combining chemical and radiative transfer models. Methods: A determination of the dust temperature in a modified Bonnor-Ebert sphere is used to calculate initial radial molecular abundance profiles. The abundances of selected cooling molecules corresponding to two different core ages are then extracted to determine the gas temperature at two time steps. The calculation is repeated in an iterative process yielding molecular abundances consistent with the gas temperature. Line emission profiles for selected substances are calculated using simulated abundance profiles. Results: The gas temperature is a function of time; the gas heats up as the core gets older because the cooling molecules are depleted onto grain surfaces. The contributions of the various cooling molecules to the total cooling power change with time. Radial chemical abundance profiles are non-trivial: different species present varying degrees of depletion and in some cases inward-increasing abundances profiles, even at t > 10^5 years. Line emission simulations indicate that cores of different ages can present significantly different line emission profiles, depending on the tracer species considered. Conclusions: Chemical abundances and the associated line cooling power change as a function of time. Most chemical species are depleted onto grain surfaces at densities exceeding ~10^5 cm^-3. Notable exceptions are NH_3 and N2H^+; the latter is largely undepleted even at n_H~10^6 cm-3. On the other hand, chemical abundances are not significantly developed in regions of low gas density even at t~10^5 years, revealed by inward-increasing abundance gradients. The gas temperature can be significantly different from the dust temperature; this may have implications on core stability.

Aims: We aim to derive a stability condition for non-isothermal Bonnor-Ebert spheres and compare the physical properties of critical non-isothermal and isothermal gas spheres. These configurations can serve as models for prestellar cores before gravitational collapse. Methods: A stability condition for non-isothermal spheres is derived by constructing an expression for the derivative of boundary pressure with respect to core volume. The temperature distribution is determined by means of radiative transfer calculations. Based on the stability analysis, we derive the physical parameters of critical cores for the mass range 0.1 - 5.0 M_sun. In addition, the properties of roughly Jupiter-mass cores are briefly examined. Results: At the low-mass end the critical non-isothermal sphere has lower central density and a slightly larger physical radius than the corresponding isothermal sphere (i.e. one with the same mass and average temperature). The temperature decrease towards the core centre becomes steeper towards smaller masses as the central density becomes higher. The slope depends on the adopted dust model. We find that the critical dimensionless radius increases above the isothermal value xi_0 = 6.45 for very low-mass cores (< 0.2 M_sun). However, in the mass-range studied here the changes are within 5% from the isothermal value. Conclusions: The density structures of non-isothermal and isothermal Bonnor-Ebert spheres for a given mass are fairly similar. However, the present models predict clear differences in the average temperatures for the same physical radius. Especially for low-mass cores, the temperature gradient probably has implications on the chemistry and the observed line emission. We also find that hydrostatic Jupiter-mass cores with radii less than 100 AU would have very high boundary pressures compared with typical pressures in the interstellar space.

Context: The depletion of heavy elements in cold cores of interstellar molecular clouds can lead to a situation where deuterated forms of H_3^+ are the most useful spectroscopic probes of the physical conditions. Aims: The aim is to predict the observability of the rotational lines of H_2D^+ and D_2H^+ from prestellar cores. Methods: Recently derived rate coefficients for the H_3^+ + H_2 isotopic system were applied to the "complete depletion" reaction scheme to calculate abundance profiles in hydrostatic core models. The ground-state lines of H_2D^+(o) (372 GHz) and D_2H^+(p) (692 GHz) arising from these cores were simulated. The excitation of the rotational levels of these molecules was approximated by using the state-to-state coefficients for collisions with H_2. We also predicted line profiles from cores with a power-law density distribution advocated in some previous studies. Results: The new rate coefficients introduce some changes to the complete depletion model, but do not alter the general tendencies. One of the modifications with respect to the previous results is the increase of the D_3^+ abundance at the cost of other isotopologues. Furthermore, the present model predicts a lower H_2D^+ (o/p) ratio, and a slightly higher D_2H^+ (p/o) ratio in very cold, dense cores, as compared with previous modelling results. These nuclear spin ratios affect the detectability of the submm lines of H_2D^+(o) and D_2H^+(p). The previously detected H_2D^+ and D_2H^+ lines towards the core I16293E, and the H_2D^+ line observed towards Oph D can be reproduced using the present excitation model and the physical models suggested in the original papers.