In this study, we explore the impact of inhomogeneities in the spatial distribution of interstellar dust on spatial scales of $\le1$ au caused by ion shadowing forces on the optical properties of diffuse interstellar medium (ISM) as well as on the dust temperature. We show that recently proposed possibility that interstellar dust grains in the diffuse ISM are grouped in spherical cloudlets (clumps) may significantly affect the observed optical properties of the diffuse ISM in comparison to that calculated under the commonly accepted assumption on the uniform dust/gas mixture if the size of clumps $\gtrsim0.1$ au. We found that opacity of an arbitrary region of diffuse ISM quickly decreases with growth of dusty clumps. We also studied the dependence of opacity and dust temperature inside the dusty clumps on their size. We show that the clumps larger than 0.1 au are opaque for far ultraviolet radiation. Dust temperature exhibit a gradient inside a clump, decreasing from the edge to the center by several degrees for a clump of a size of 0.1 au and larger. We argue that dust temperatures and high opacity within clumps larger than 0.1 au may facilitate somewhat more efficient synthesis of molecules on surfaces of interstellar grains in the diffuse ISM than it was anticipated previously. On the other hand, the presence of clumps with sizes below 0.1 au makes small or negligible influence on the optical properties of the diffuse ISM in comparison to the case with uniform dust/gas mixture.

We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particularly efficient. Methanol and sulphur monoxide are most abundant in a narrow zone that follows the eastern side of the core. This side is sheltered from the stronger external radiation field coming from the west. We show that photodissociation on the illuminated side can give rise to an asymmetric methanol distribution, but that the stark contrast observed in H-MM1 is hard to explain without assuming enhanced desorption on the shaded side. The region of the brightest emission has a wavy structure that rolls up at one end. This is the signature of Kelvin-Helmholtz instability occurring in sheared flows. We suggest that in this zone, methanol and sulphur are released as a result of grain-grain collisions induced by shear vorticity.

The SOLIS (Seeds Of Life In Space) IRAM/NOEMA Large Program aims at studying a set of crucial complex organic molecules in a sample of sources, with well-known physical structure, covering the various phases of Solar-type star formation. One representative object of the transition from the prestellar core to the protostar phases has been observed toward the Very Low Luminosity Object (VeLLO) called L1521F. This type of source is important to study to make the link between prestellar cores and Class 0 sources and also to constrain the chemical evolution during the process of star formation. Two frequency windows (81.6-82.6 GHz and 96.65-97.65 GHz) were used to observe the emission from several complex organics toward the L1521F VeLLO. Only 2 transitions of methanol (A+, E2) have been detected in the narrow window centered at 96.7 GHz (with an upper limit on E1) in a very compact emission blob (~7'' corresponding to ~1000au) toward the NE of the L1521F protostar. The CS 2-1 transition is also detected within the WideX bandwidth. Consistently, with what has been found in prestellar cores, the methanol emission appears ~1000au away from the dust peak. The location of the methanol blob coincides with one of the filaments previously reported in the literature. The Tex of the gas inferred from methanol is (10$\pm$2) K, while the H2 gas density (estimated from the detected CS 2-1 emission and previous CS 5-4 ALMA obs.) is a factor >25 higher than the density in the surrounding environment (n(H2) >10$^{7}$ cm$^{-3}$). From its compactness, low excitation temperature and high gas density, we suggest that the methanol emission detected with NOEMA is either a cold and dense shock-induced blob, recently formed ($\leq$ few hundred years) by infalling gas or a cold and dense fragment that may have just been formed as a result of the intense gas dynamics found within the L1521F VeLLO system.

To understand how Phosphorus-bearing molecules are formed in star-forming regions, we have analysed ALMA observations of PN and PO towards the massive star-forming region AFGL 5142, combined with a new analysis of the data of the comet 67P/Churyumov-Gerasimenko taken with the ROSINA instrument onboard Rosetta. The ALMA maps show that the emission of PN and PO arises from several spots associated with low-velocity gas with narrow linewidths in the cavity walls of a bipolar outflow. PO is more abundant than PN in most of the spots, with the PO/PN ratio increasing as a function of the distance to the protostar. Our data favor a formation scenario in which shocks sputter phosphorus from the surface of dust grains, and gas-phase photochemistry induced by UV photons from the protostar allows efficient formation of the two species in the cavity walls. Our analysis of the ROSINA data has revealed that PO is the main carrier of P in the comet, with PO/PN>10. Since comets may have delivered a significant amount of prebiotic material to the early Earth, this finding suggests that PO could contribute significantly to the phosphorus reservoir during the dawn of our planet. There is evidence that PO was already in the cometary ices prior to the birth of the Sun, so the chemical budget of the comet might be inherited from the natal environment of the Solar System, which is thought to be a stellar cluster including also massive stars.

In this work, we reexamine sulfur chemistry occurring on and in the ice mantles of interstellar dust grains, and report the effects of two new modifications to standard astrochemical models; namely, (a) the incorporation of cosmic ray-driven radiation chemistry and (b) the assumption of fast, non-diffusive reactions for key radicals in the bulk. Results from our models of dense molecular clouds show that these changes can have a profound influence on the abundances of sulfur-bearing species in ice mantles, including a reduction in the abundance of solid-phase H$_2$S and HS, and a significant increase in the abundances of OCS, SO$_2$, as well as pure allotropes of sulfur, especially S$_8$. These pure-sulfur species - though nearly impossible to observe directly - have long been speculated to be potential sulfur reservoirs and our results represent possibly the most accurate estimates yet of their abundances in the dense ISM. Moreover, the results of these updated models are found to be in good agreement with available observational data. Finally, we examine the implications of our findings with regard to the as-yet-unknown sulfur reservoir thought to exist in dense interstellar environments.

Phosphorus-bearing species are an essential key to form life on Earth, however they have barely been detected in the interstellar medium. Since only PN and PO have been identified so far towards star-forming regions, the chemical formation pathways of P-bearing molecules are not easy to constrain and are thus highly debatable. An important factor still missing in the chemical models is the initial elemental abundance of phosphorus, i.e. the depletion level of P. In order to overcome this problem, we study P-bearing species in diffuse/translucent clouds. In these objects phosphorus is mainly in the gas phase and therefore the elemental initial abundance needed in our chemical simulations corresponds to the cosmic one and is thus well constrained. An advanced chemical model with an updated P-chemistry network has been used. Single-pointing observations were performed with the IRAM 30m telescope towards the line of sight to the blazar B0355+508 aiming for the (2-1) transitions of PN, PO, HCP and CP. This line of sight incorporates five diffuse/translucent clouds. The (2-1) transitions of the PN, PO, HCP and CP were not detected. We report detections of the (1-0) lines of $\mathrm{^{13}CO}$, HNC and CN along with a first detection of $\mathrm{C^{34}S}$. We have reproduced the observations of HNC, CN, CS and CO in every cloud by applying typical conditions for diffuse/translucent clouds. According to our best-fit model, the most abundant P-bearing species are HCP and CP ($\sim10^{-10}$), followed by PN, PO and $\mathrm{PH_3}$ ($\sim10^{-11}$). We show that the production of P-bearing species is favoured towards translucent rather than diffuse clouds, where the environment provides a stronger shielding from the interstellar radiation. Based on our improved model, the (1-0) transitions of HCP, CP, PN and PO are expected to be detectable with estimated intensities up to $\sim200$ mK.

L1521E is a dense starless core in Taurus that was found to have relatively low molecular depletion by earlier studies, thus suggesting a recent formation. We aim to characterize the chemical structure of L1521E and compare it to the more evolved L1544 pre-stellar core. We have obtained $\sim$2.5$\times$2.5 arcminute maps toward L1521E using the IRAM-30m telescope in transitions of various species. We derived abundances for the species and compared them to those obtained toward L1544. We estimated CO depletion factors. Similarly to L1544, $c$-C$_3$H$_2$ and CH$_3$OH peak at different positions. Most species peak toward the $c$-C$_3$H$_2$ peak. The CO depletion factor derived toward the $Herschel$ dust peak is 4.3$\pm$1.6, which is about a factor of three lower than that toward L1544. The abundances of sulfur-bearing molecules are higher toward L1521E than toward L1544 by factors of $\sim$2-20. The abundance of methanol is similar toward the two cores. The higher abundances of sulfur-bearing species toward L1521E than toward L1544 suggest that significant sulfur depletion takes place during the dynamical evolution of dense cores, from the starless to pre-stellar stage. The CO depletion factor measured toward L1521E suggests that CO is more depleted than previously found. Similar CH$_3$OH abundances between L1521E and L1544 hint that methanol is forming at specific physical conditions in Taurus, characterized by densities of a few $\times$10$^4$ cm$^{-3}$ and $N$(H$_2$)$\gtrsim$10$^{22}$ cm$^{-2}$, when CO starts to catastrophically freeze-out, while water can still be significantly photodissociated, so that the surfaces of dust grains become rich in solid CO and CH$_3$OH, as already found toward L1544. Methanol can thus provide selective crucial information about the transition region between dense cores and the surrounding parent cloud.

Many astrochemical models today explicitly consider the species that comprise the bulk of interstellar dust grain ice-mantles separately from those in the top few monolayers. Bombardment of these ices by ionizing radiation - whether in the form of cosmic rays, stellar winds, or radionuclide emission - represents an astrochemically viable means of driving a rich chemistry even in the bulk of the ice-mantle, now supported by a large body of work in laboratory astrophysics. In this study, using an existing rate equation-based astrochemical code modified to include a method of considering radiation chemistry recently developed by us, we attempted to simulate two such studies in which (a) pure O$_2$ ice at 5 K and, (b) pure H$_2$O ice at 16 K and 77 K, were bombarded by keV H$^+$ ions. Our aims are twofold: (1) to test the capability of our newly developed method to replicate the results of ice-irradiation experiments, and (2) to determine in such a well-constrained system how bulk chemistry is best handled using the same gas-grain codes that are used to model the interstellar medium (ISM). We find that our modified astrochemical model is able to reproduce both the abundance of O$_3$ in the 5 K pure O$_2$ ice, as well as both the abundance of H$_2$O$_2$ in the 16 K water ice and the previously noted decrease of hydrogen peroxide at higher temperatures. However, these results require the assumption that radicals and other reactive species produced via radiolysis react quickly and non-diffusively with neighbors in the ice.

We aim to reveal the physical properties and chemical composition of the cores in the California molecular cloud (CMC), so as to better understand the initial conditions of star formation. We made a high-resolution column density map (18.2") with Herschel data, and extracted a complete sample of the cores in the CMC with the \textsl{fellwalker} algorithm. We performed new single-pointing observations of molecular lines near 90 GHz with the IRAM 30m telescope along the main filament of the CMC. In addition, we also performed a numerical modeling of chemical evolution for the cores under the physical conditions. We extracted 300 cores, of which 33 are protostellar and 267 are starless cores. About 51\% (137 of 267) of the starless cores are prestellar cores. Three cores have the potential to evolve into high-mass stars. The prestellar core mass function (CMF) can be well fit by a log-normal form. The high-mass end of the prestellar CMF shows a power-law form with an index $\alpha=-0.9\pm 0.1$ that is shallower than that of the Galactic field stellar mass function. Combining the mass transformation efficiency ($\varepsilon$) from the prestellar core to the star of $15\pm 1\%$ and the core formation efficiency (CFE) of 5.5\%, we suggest an overall star formation efficiency of about 1\% in the CMC. In the single-pointing observations with the IRAM 30m telescope, we find that 6 cores show blue-skewed profile, while 4 cores show red-skewed profile. [$\rm {HCO}^{+}$]/[HNC] and [$\rm {HCO}^{+}$]/$\rm [N_{2}H^{+}]$ in protostellar cores are higher than those in prestellar cores; this can be used as chemical clocks. The best-fit chemical age of the cores with line observations is $\sim 5\times 10^4$~years.

The ketenyl radical, HCCO, has recently been detected in the ISM for the first time. Further astronomical detections of HCCO will help us understand its gas-grain chemistry, and subsequently revise the oxygen-bearing chemistry towards dark clouds. Moreover, its deuterated counterpart, DCCO, has never been observed in the ISM. HCCO and DCCO still lack a broad spectroscopic investigation, although they exhibit a significant astrophysical relevance. In this work we aim to measure the pure rotational spectra of the ground state of HCCO and DCCO in the millimeter and submillimeter region, considerably extending the frequency range covered by previous studies. The spectral acquisition was performed using a frequency-modulation absorption spectrometer between 170 and 650 GHz. The radicals were produced in a low-density plasma generated from a select mixture of gaseous precursors. For each isotopologue we were able to detect and assign more than 100 rotational lines. The new lines have significantly enhanced the previous data set allowing the determination of highly precise rotational and centrifugal distortion parameters. In our analysis we have taken into account the interaction between the ground electronic state and a low-lying excited state (Renner-Teller pair) which enables the prediction and assignment of rotational transitions with $K_a$ up to 4. The present set of spectroscopic parameters provides highly accurate, millimeter and submillimeter rest-frequencies of HCCO and DCCO for future astronomical observations. We also show that towards the pre-stellar core L1544, ketenyl peaks in the region where $c$-$\mathrm{C_3H_2}$ peaks, suggesting that HCCO follows a predominant hydrocarbon chemistry, as already proposed by recent gas-grain chemical models.

The formyl radical HCO has been proposed as the basic precursor of many complex organic molecules such as methanol (CH$_3$OH) or glycolaldehyde (CH$_2$OHCHO). Using ALMA, we have mapped, for the first time at high angular resolution ($\sim$1$^{\prime\prime}$, $\sim$140 au), HCO towards the Solar-type protostellar binary IRAS 16293$-$2422, where numerous complex organic molecules have been previously detected. We also detected several lines of the chemically related species H$_2$CO, CH$_3$OH and CH$_2$OHCHO. The observations revealed compact HCO emission arising from the two protostars. The line profiles also show redshifted absorption produced by foreground material of the circumbinary envelope that is infalling towards the protostars. Additionally, IRAM 30m single-dish data revealed a more extended HCO component arising from the common circumbinary envelope. The comparison between the observed molecular abundances and our chemical model suggests that whereas the extended HCO from the envelope can be formed via gas-phase reactions during the cold collapse of the natal core, the HCO in the hot corinos surrounding the protostars is predominantly formed by the hydrogenation of CO on the surface of dust grains and subsequent thermal desorption during the protostellar phase. The derived abundance of HCO in the dust grains is high enough to produce efficiently more complex species such as H$_2$CO, CH$_3$OH, and CH$_2$OHCHO by surface chemistry. We found that the main formation route of CH$_2$OHCHO is the reaction between HCO and CH$_2$OH.

Aims. In this paper we investigate the detectability of the molecular oxygen in icy dust grain mantles towards astronomical objects. Methods. We present a systematic set of experiments with O2-H2O ice mixtures designed to disentangle how the molecular ratio affects the O2 signature in the mid- and near-infrared spectral regions. All the experiments were conducted in a closed-cycle helium cryostat coupled to a Fourier transform infrared spectrometer. The ice mixtures comprise varying thicknesses from 8 $\times$ 10$^{-3}$ to 3 $\mu$m. The absorption spectra of the O2-H2O mixtures are also compared to the one of pure water. In addition, the possibility to detect the O2 in icy bodies and in the interstellar medium is discussed. Results. We are able to see the O2 feature at 1551 cm$^{-1}$ even for the most diluted mixture of H2O : O2 = 9 : 1, comparable to a ratio of O2/H2O = 10 % which has already been detected in situ in the coma of the comet 67P/Churyumov-Gerasimenko. We provide an estimate for the detection of O2 with the future mission of the James Webb Space Telescope (JWST).

The study of deuteration in pre-stellar cores is important to understand the physical and chemical initial conditions in the process of star formation. In particular, observations toward pre-stellar cores of methanol and deuterated methanol, solely formed on the surface of dust grains, may provide useful insights on surface processes at low temperatures. Here we analyze maps of CO, methanol, formaldehyde and their deuterated isotopologues toward a well-known pre-stellar core. This study allows us to test current gas-dust chemical models. Single-dish observations of CH$_3$OH, CH$_2$DOH, H$_2$CO, H$_2\,^{13}$CO, HDCO, D$_2$CO and C$^{17}$O toward the prototypical pre-stellar core L1544 were performed at the IRAM 30 m telescope. We analyze their column densities, distributions, and compare these observations with gas-grain chemical models. The maximum deuterium fraction derived for methanol is [CH$_2$DOH]/[CH$_3$OH] $\sim$ 0.08$\pm$0.02, while the measured deuterium fractions of formaldehyde at the dust peak are [HDCO]/[H$_2$CO] $\sim$ 0.03$\pm$0.02, [D$_2$CO]/[H$_2$CO] $\sim$ 0.04$\pm$0.03 and [D$_2$CO]/[HDCO] $\sim$ 1.2$\pm$0.3. Observations differ significantly from the predictions of models, finding discrepancies between a factor of 10 and a factor of 100 in most cases. It is clear though that to efficiently produce methanol on the surface of dust grains, quantum tunneling diffusion of H atoms must be switched on. It also appears that the currently adopted reactive desorption efficiency of methanol is overestimated and/or that abstraction reactions play an important role. More laboratory work is needed to shed light on the chemistry of methanol, an important precursor of complex organic molecules in space.

Phosphorus is a crucial element in prebiotic chemistry, especially the P$-$O bond, which is key for the formation of the backbone of the deoxyribonucleic acid. So far, PO had only been detected towards the envelope of evolved stars, and never towards star-forming regions. We report the first detection of PO towards two massive star-forming regions, W51 e1/e2 and W3(OH), using data from the IRAM 30m telescope. PN has also been detected towards the two regions. The abundance ratio PO/PN is 1.8 and 3 for W51 and W3(OH), respectively. Our chemical model indicates that the two molecules are chemically related and are formed via gas-phase ion-molecule and neutral-neutral reactions during the cold collapse. The molecules freeze out onto grains at the end of the collapse and desorb during the warm-up phase once the temperature reaches $\sim$35 K. The observed molecular abundances of 10$^{-10}$ are predicted by the model if a relatively high initial abundance of phosphorus, 5$\times$10$^{-9}$, is assumed.

The formation of compact dusty clouds in diffuse interstellar medium (ISM) has been recently proposed and studied by Tsytovich et al. (2014). In the present paper, an effect of the clouds on the rate of H$\to$H$_2$ transition in the ISM is examined. We discuss the mechanisms leading to the formation of the clouds and the creation of gaseous clumps around them, analyze the main processes determining the efficiency of the recombination of atomic hydrogen on dust grains, and estimate the expected enhancement of the global H$_2$ formation due to the presence of the clouds. In conclusion, we argue that the compact dusty clouds provide optimum conditions for the atomic recombination in diffuse ISM, and point out some astrophysical implications of the resulting H$_2$ formation enhancement.

We present multi-transition observations of PN towards a sample of nine massive dense cores in different evolutionary stages. Using transitions with different excitation conditions, we have found for the first time that the excitation temperatures of PN are in the range 5-30 K. To investigate the main chemical route for the PN formation (surface-chemistry vs. gas-phase chemistry), and the dominant desorption mechanism (thermal vs. shock), we have compared our results with those obtained from molecules tracing different chemical and physical conditions (SiO, SO, CH3OH, and N2H+). We have found that the PN line profiles are very well correlated with those of SiO and SO in six out of the nine targets, which indicates that PN may be released by sputtering of dust grains due to shocks. This finding is corroborated by a faint but statistically significant positive trend between the PN abundance and those of SiO and SO. However, in three objects the PN lines have no hints of high velocity wings, which indicates an alternative origin of PN. Overall, our results indicate that the origin of PN is not unique, as it can be formed in protostellar shocks, but also in colder and more quiescent gas through alternative pathways.

Towards the pre-stellar core L1544, the methanol (CH$_3$OH) emission forms an asymmetric ring around the core centre, where CH$_3$OH is mostly in solid form, with a clear peak 4000~au to the north-east of the dust continuum peak. As part of the NOEMA Large Project SOLIS (Seeds of Life in Space), the CH$_3$OH peak has been spatially resolved to study its kinematics and physical structure and to investigate the cause behind the local enhancement. We find that methanol emission is distributed in a ridge parallel to the main axis of the dense core. The centroid velocity increases by about 0.2~km~s$^{-1}$ and the velocity dispersion increases from subsonic to transonic towards the central zone of the core, where the velocity field also shows complex structure. This could be indication of gentle accretion of material onto the core or interaction of two filaments, producing a slow shock. We measure the rotational temperature and show that methanol is in local thermodynamic equilibrium (LTE) only close to the dust peak, where it is significantly depleted. The CH$_3$OH column density, $N_{tot}({\rm CH_3OH})$, profile has been derived with non-LTE radiative transfer modelling and compared with chemical models of a static core. The measured $N_{tot}({\rm CH_3OH})$ profile is consistent with model predictions, but the total column densities are one order of magnitude lower than those predicted by models, suggesting that the efficiency of reactive desorption or atomic hydrogen tunnelling adopted in the model may be overestimated; or that an evolutionary model is needed to better reproduce methanol abundance.

Complex organic molecules have been observed for decades in the interstellar medium. Some of them might be considered as small bricks of the macromolecules at the base of terrestrial life. It is hence particularly important to understand organic chemistry in Solar-like star forming regions. In this article, we present a new observational project: SOLIS (Seeds Of Life In Space). This is a Large Project at the IRAM-NOEMA interferometer, and its scope is to image the emission of several crucial organic molecules in a sample of Solar-like star forming regions in different evolutionary stage and environments. Here, we report the first SOLIS results, obtained from analysing the spectra of different regions of the Class 0 source NGC1333-IRAS4A, the protocluster OMC-2 FIR4, and the shock site L1157-B1. The different regions were identified based on the images of formamide (NH2CHO) and cyanodiacetylene (HC5N) lines. We discuss the observed large diversity in the molecular and organic content, both on large (3000-10000 au) and relatively small (300-1000 au) scales. Finally, we derive upper limits to the methoxy fractional abundance in the three observed regions of the same order of magnitude of that measured in few cold prestellar objects, namely ~10^-12-10^-11 with respect to H2 molecules.

Protoplanetary disk mass is a key parameter controlling the process of planetary system formation. CO molecular emission is often used as a tracer of gas mass in the disk. In this study we consider the ability of CO to trace the gas mass over a wide range of disk structural parameters and search for chemical species that could possibly be used as alternative mass tracers to CO. Specifically, we apply detailed astrochemical modeling to a large set of models of protoplanetary disks around low-mass stars, to select molecules with abundances correlated with the disk mass and being relatively insensitive to other disk properties. We do not consider sophisticated dust evolution models, restricting ourselves with the standard astrochemical assumption of $0.1~\mu $m dust. We find that CO is indeed the best molecular tracer for total gas mass, despite the fact that it is not the main carbon carrier, provided reasonable assumptions about CO abundance in the disk are used. Typically, chemical reprocessing lowers the abundance of CO by a factor of 3, compared to the case of photo-dissociation and freeze-out as the only ways of CO depletion. On average only 13% C-atoms reside in gas-phase CO, albeit with variations from 2 to 30%. CO$_2$, H$_2$O and H$_2$CO can potentially serve as alternative mass tracers, the latter two being only applicable if disk structural parameters are known.

Context: Modern versions of the Miller-Urey experiment claim that formamide (NH$_2$CHO) could be the starting point for the formation of metabolic and genetic macromolecules. Intriguingly, formamide is indeed observed in regions forming Solar-type stars as well as in external galaxies. Aims: How NH$_2$CHO is formed has been a puzzle for decades: our goal is to contribute to the hotly debated question of whether formamide is mostly formed via gas-phase or grain surface chemistry. Methods: We used the NOEMA interferometer to image NH$_2$CHO towards the L1157-B1 blue-shifted shock, a well known interstellar laboratory, to study how the components of dust mantles and cores released into the gas phase triggers the formation of formamide. Results: We report the first spatially resolved image (size $\sim$ 9", $\sim$ 2300 AU) of formamide emission in a shocked region around a Sun-like protostar: the line profiles are blueshifted and have a FWHM $\simeq$ 5 km s$^{-1}$. A column density of $N_{\rm NH_2CHO}$ = 8 $\times$ 10$^{12}$ cm$^{-1}$, and an abundance (with respect to H-nuclei) of 4 $\times$ 10$^{-9}$ are derived. We show a spatial segregation of formamide with respect to other organic species. Our observations, coupled with a chemical modelling analysis, indicate that the formamide observed in L1157-B1 is formed by gas-phase chemical process, and not on grain surfaces as previously suggested. Conclusions: The SOLIS interferometric observations of formamide provide direct evidence that this potentially crucial brick of life is efficiently formed in the gas-phase around Sun-like protostars.

The interstellar delivery of carbon atoms locked into molecules might be one of the key ingredients for the emergence of life. Cyanopolyynes are carbon chains delimited at their two extremities by an atom of hydrogen and a cyano group, so that they might be excellent reservoirs of carbon. The simplest member, HC3N, is ubiquitous in the galactic interstellar medium and found also in external galaxies. Thus, understanding the growth of cyanopolyynes in regions forming stars similar to our Sun, and what affects it, is particularly relevant. In the framework of the IRAM/NOEMA Large Program SOLIS (Seeds Of Life In Space), we have obtained a map of two cyanopolyynes, HC3N and HC5N, in the protocluster OMC2-FIR4. Because our Sun is thought to be born in a rich cluster, OMC2-FIR4 is one of the closest and best known representatives of the environment in which the Sun may have been born. We find a HC3N/HC5N abundance ratio across the source in the range ~ 1 - 30, with the smallest values (< 10) in FIR5 and in the Eastern region of FIR4. The ratios < 10 can be reproduced by chemical models only if: (1) the cosmic-ray ionisation rate $\zeta$ is ~ $4 \times 10^{-14}$ s$^{-1}$; (2) the gaseous elemental ratio C/O is close to unity; (3) O and C are largely depleted. The large $\zeta$ is comparable to that measured in FIR4 by previous works and was interpreted as due to a flux of energetic (> 10 MeV) particles from embedded sources. We suggest that these sources could lie East of FIR4 and FIR5. A temperature gradient across FIR4, with T decreasing by about 10 K, could also explain the observed change in the HC3N/HC5N line ratio, without the need of a cosmic ray ionisation rate gradient. However, even in this case, a high constant cosmic-ray ionisation rate (of the order of $10^{-14}$ s$^{-1}$) is necessary to reproduce the observations.

We present the results of chemical modeling of complex organic molecules (COMs) under conditions typical for prestellar cores. We utilize an advanced gas-grain astrochemical model with updated gas-phase chemistry, with a multilayer approach to ice-surface chemistry and an up-to-date treatment of reactive desorption based on recent experiments of Minissale et al. (2016). With the chemical model, radial profiles of molecules including COMs are calculated for the case of the prototypical prestellar core L1544 at the timescales when the modeled depletion factor of CO becomes equal to that observed. We find that COMs can be formed efficiently in L1544 up to the fractional abundances of 10(-10) wrt. total hydrogen nuclei. Abundances of many COMs such as CH3OCH3, HCOOCH3, and others peak at similar radial distances of ~2000-4000 AU. Gas-phase abundances of COMs depend on the efficiency of reactive desorption, which in turn depends on the composition of the outer monolayers of icy mantles. In prestellar cores, the outer monolayers of mantles likely include large fractions of CO and its hydrogenation products, which may increase the efficiency of reactive desorption according to Minissale et al. (2016), and makes the formation of COMs efficient under conditions typical for prestellar cores, although this assumption is yet to be confirmed experimentally. The hydroxyl radical (OH) appears to play an important role in gas-phase chemistry of COMs, which makes it deserving further detailed studies.

The detection of complex organic molecules (COMs) toward cold sources such as pre-stellar cores (with T<10 K), has challenged our understanding of the formation processes of COMs in the interstellar medium. Recent modelling on COM chemistry at low temperatures has provided new insight into these processes predicting that COM formation depends strongly on parameters such as visual extinction and the level of CO freeze out. We report deep observations of COMs toward two positions in the L1544 pre-stellar core: the dense, highly-extinguished continuum peak with Av>=30 mag within the inner 2700 au; and a low-density shell with average Av~7.5-8 mag located at 4000 au from the core's center and bright in CH3OH. Our observations show that CH3O, CH3OCH3 and CH3CHO are more abundant (by factors ~2-10) toward the low-density shell than toward the continuum peak. Other COMs such as CH3OCHO, c-C3H2O, HCCCHO, CH2CHCN and HCCNC show slight enhancements (by factors <=3) but the associated uncertainties are large. This suggests that COMs are actively formed and already present in the low-density shells of pre-stellar cores. The modelling of the chemistry of O-bearing COMs in L1544 indicates that these species are enhanced in this shell because i) CO starts freezing out onto dust grains driving an active surface chemistry; ii) the visual extinction is sufficiently high to prevent the UV photo-dissociation of COMs by the external interstellar radiation field; and iii) the density is still moderate to prevent severe depletion of COMs onto grains.

Phosphorus is a crucial element in biochemistry, especially the P-O bond, which is key for the formation of the backbone of the deoxyribonucleic acid. So far, PO has only been detected towards the envelope of evolved stars, and never towards star-forming regions. We report the first detections of PO towards two massive star-forming regions, W51 e1/e2 and W3(OH), using data from the IRAM 30m telescope. PN has also been detected towards the two regions. The abundance ratio PO/PN is 1.8 and 3 for W51 and W3(OH), respectively. Our chemical model indicates that the two molecules are chemically related and are formed via gas-phase ion-molecule and neutral-neutral reactions during the cold collapse. The molecules freeze out onto grains at the end of the collapse and desorb during the warm-up phase once the temperature reaches 35 K. Similar abundances of the two species are expected during a period of 5x10^{4} yr at the early stages of the warm-up phase, when the temperature is in the range 35-90 K. The observed molecular abundances of 10^{-10} are predicted by the model if a relatively high initial abundance of 5x10^{-9} of depleted phosphorus is assumed.

Phosphorus is a crucial element for the development of life, but so far P-bearing molecules have been detected only in a few astrophysical objects, hence its interstellar chemistry is almost totally unknown. Here we show new detections of phosphorus nitride in a sample of dense cores in different evolutionary stages of the intermediate- and high-mass star formation process: starless, with protostellar objects, and with ultracompact HII regions. All detected PN line widths are smaller than ~5 km/s , and they arise from regions associated with kinetic temperatures smaller than 100 K. Because the few previous detections reported in the literature are associated with warmer and more turbulent sources, the results of this work show that PN can arise from relatively quiescent and cold gas. This information is challenging for theoretical models that invoke either high desorption temperatures or grain sputtering from shocks to release phosphorus into the gas phase. Derived column densities are of the order of 10^{11-12} cm^{-2}, marginally lower than the values derived in the few high-mass star forming regions detected so far. New constraints on the column density upper limit of phosphorus monoxide, the fundamental unit of biologically relevant molecules, are also given.

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.

The problem of impulsive heating of dust grains in cold, dense interstellar clouds is revisited theoretically, with the aim to better understand leading mechanisms of the explosive desorption of icy mantles. It is rigorously shown that if the heating of a reactive medium occurs within a sufficiently localized spot (e.g., heating of mantles by cosmic rays), then the subsequent thermal evolution is characterized by a single dimensionless number $\lambda$. This number identifies a bifurcation between two distinct regimes: When $\lambda$ exceeds a critical value (threshold), the heat equation exhibits the explosive solution, i.e., the thermal (chemical) explosion is triggered. Otherwise, thermal diffusion causes the deposited heat to spread over the entire grain -- this regime is commonly known as the whole-grain heating. The theory allows us to find a critical combination of the physical parameters that govern the explosion of icy mantles due to impulsive spot heating. In particular, the calculations suggest that heavy cosmic ray species (e.g., iron ions) colliding with dust are able to trigger the explosion. Based on the recently calculated local cosmic-ray spectra, the expected rate of the explosive desorption is estimated. The efficiency of the desorption, which affects all solid species independent of their binding energy, is shown to be comparable with other cosmic-ray desorption mechanisms typically considered in the literature. Also, the theory allows us to estimate maximum abundances of reactive species that may be stored in the mantles, which provides important constraints on available astrochemical models.

Core accretion models of massive star formation require the existence of stable massive starless cores, but robust observational examples of such objects have proven elusive. We report subarcsecond-resolution SMA 1.3 mm, 1.1 mm, and 0.88 mm and VLA 1.3 cm observations of an excellent massive starless core candidate, G11.92-0.61-MM2, initially identified in the course of studies of GLIMPSE Extended Green Objects (EGOs). Separated by ~7.2" from the nearby MM1 protostellar hot core, MM2 is a strong, compact dust continuum source (submillimeter spectral index alpha=2.6+/-0.1), but is devoid of star formation indicators. In contrast to MM1, MM2 has no masers, no centimeter continuum, and no (sub)millimeter wavelength line emission in ~24 GHz of bandwidth observed with the SMA, including N2H+(3-2), HCO+(3-2), and HCN(3-2). Additionally, there is no evidence for an outflow driven by MM2. The (sub)millimeter spectral energy distribution (SED) of MM2 is best fit with a dust temperature of ~17-19 K and luminosity of ~5-7 L_sun. The combined physical properties of MM2, as inferred from its dust continuum emission, are extreme: M>30 M_sun within a radius<1000 AU, N(H2)>10^25 cm^-2 and n(H2)>10^9 cm^-3. Comparison of the molecular abundance limits derived from our SMA observations with gas-grain chemical models indicates that extremely dense (n(H)>>10^8 cm^-3), cold (<20 K) conditions are required to explain the lack of observed (sub)millimeter line emission, consistent with the dust continuum results. Our data suggest that G11.92-0.61-MM2 is the best candidate for a bonafide massive prestellar core found to date, and a promising target for future, higher-sensitivity observations.

HNC and HCN, typically used as dense gas tracers in molecular clouds, are a pair of isomers that have great potential as a temperature probe because of temperature dependent, isomer-specific formation and destruction pathways. Previous observations of the HNC/HCN abundance ratio show that the ratio decreases with increasing temperature, something that standard astrochemical models cannot reproduce. We have undertaken a detailed parameter study on which environmental characteristics and chemical reactions affect the HNC/HCN ratio and can thus contribute to the observed dependence. Using existing gas and gas-grain models updated with new reactions and reaction barriers, we find that in static models the H + HNC gas-phase reaction regulates the HNC/HCN ratio under all conditions, except for very early times. We quantitively constrain the combinations of H abundance and H + HNC reaction barrier that can explain the observed HNC/HCN temperature dependence and discuss the implications in light of new quantum chemical calculations. In warm-up models, gas-grain chemistry contributes significantly to the predicted HNC/HCN ratio and understanding the dynamics of star formation is therefore key to model the HNC/HCN system.

[Abridged] Ethylene oxide and its isomer acetaldehyde are important complex organic molecules because of their potential role in the formation of amino acids. Despite the fact that acetaldehyde is ubiquitous in the interstellar medium, ethylene oxide has not yet been detected in cold sources. We aim to understand the chemistry of the formation and loss of ethylene oxide in hot and cold interstellar objects (i) by including in a revised gas-grain network some recent experimental results on grain surfaces and (ii) by comparison with the chemical behaviour of its isomer, acetaldehyde. We test the code for the case of a hot core. The model allows us to predict the gaseous and solid ethylene oxide abundances during a cooling-down phase prior to star formation and during the subsequent warm-up phase. We can therefore predict at what temperatures ethylene oxide forms on grain surfaces and at what temperature it starts to desorb into the gas phase. The model reproduces the observed gaseous abundances of ethylene oxide and acetaldehyde towards high-mass star-forming regions. In addition, our results show that ethylene oxide may be present in outer and cooler regions of hot cores where its isomer has already been detected. Despite their different chemical structures, the chemistry of ethylene oxide is coupled to that of acetaldehyde, suggesting that acetaldehyde may be used as a tracer for ethylene oxide towards cold cores.

(Abridged) Protoplanetary disks are vital objects in star and planet formation, possessing all the material which may form a planetary system orbiting the new star. We investigate the synthesis of complex organic molecules (COMs) in disks to constrain the achievable chemical complexity and predict species and transitions which may be observable with ALMA. We have coupled a 2D model of a protoplanetary disk around a T Tauri star with a gas-grain chemical network including COMs. We compare compare synthesised line intensities and calculated column densities with observations and determine those COMs which may be observable in future. COMs are efficiently formed in the disk midplane via grain-surface chemical reactions, reaching peak grain-surface fractional abundances 1e-6 - 1e-4 that of the H nuclei number density. COMs formed on grain surfaces are returned to the gas phase via non-thermal desorption; however, gas-phase species reach lower fractional abundances than their grain-surface equivalents, 1e-12 - 1e-7. Including the irradiation of grain mantle material helps build further complexity in the ice through the replenishment of grain-surface radicals which take part in further grain-surface reactions. There is reasonable agreement with several line transitions of H2CO observed towards several T Tauri star-disk systems. The synthesised line intensities for CH3OH are consistent with upper limits determined towards all sources. Our models suggest CH3OH should be readily observable in nearby protoplanetary disks with ALMA; however, detection of more complex species may prove challenging. Our grain-surface abundances are consistent with those derived from cometary comae observations providing additional evidence for the hypothesis that comets (and other planetesimals) formed via the coagulation of icy grains in the Sun's natal disk.

It is currently assumed that infrared dark clouds (IRDCs) represent the earliest evolutionary stages of high-mass stars ($>$ 8 M$_{\odot}$). Submillimeter and millimeter-wave studies performed over the past 15 years show that IRDCs possess a broad variety of properties, and hence a wide range of problems and questions that can be tackled. In this paper, we report an investigation of the molecular composition and chemical processes in two groups of IRDCs. Using the Mopra, APEX, and IRAM radio telescopes over the last four years, we have collected molecular line data for CO, H$_2$CO, HNCO, CH$_3$CCH, CH$_3$OH, CH$_3$CHO, CH$_3$OCHO, and CH$_3$OCH$_3$. For all of these species we estimated molecular abundances. We then undertook chemical modeling studies, concentrating on the source IRDC028.34+0.06, and compared observed and modeled abundances. This comparison showed that to reproduce observed abundances of complex organic molecules (COMs), a 0-D gas-grain model with constant physical conditions is not sufficient. We achieved greater success with the use of a warm-up model, in which warm-up from 10 K to 30 K occurs following a cold phase.

The chemical changes of high-mass star-forming regions provide a potential method for classifying their evolutionary stages and, ultimately, ages. In this study, we search for correlations between molecular abundances and the evolutionary stages of dense molecular clumps associated with high-mass star formation. We use the molecular line maps from Year 1 of the Millimetre Astronomy Legacy Team 90 GHz (MALT90) Survey. The survey mapped several hundred individual star-forming clumps chosen from the ATLASGAL survey to span the complete range of evolution, from prestellar to protostellar to H II regions. The evolutionary stage of each clump is classified using the Spitzer GLIMPSE/MIPSGAL mid-IR surveys. Where possible, we determine the dust temperatures and H2 column densities for each clump from Herschel Hi-GAL continuum data. From MALT90 data, we measure the integrated intensities of the N2H+, HCO+, HCN and HNC (1-0) lines, and derive the column densities and abundances of N2H+ and HCO+. The Herschel dust temperatures increase as a function of the IR-based Spitzer evolutionary classification scheme, with the youngest clumps being the coldest, which gives confidence that this classification method provides a reliable way to assign evolutionary stages to clumps. Both N2H+ and HCO+ abundances increase as a function of evolutionary stage, whereas the N2H+ (1-0) to HCO+ (1-0) integrated intensity ratios show no discernable trend. The HCN (1-0) to HNC(1-0) integrated intensity ratios show marginal evidence of an increase as the clumps evolve.

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.

The recent discovery of terrestrial-type organic species such as methyl formate and dimethyl ether in the cold interstellar gas has proved that the formation of organic matter in the Galaxy begins at a much earlier stage of star formation than was thought before. This discovery represents a challenge for astrochemical modelers. The abundances of these molecules cannot be explained by the previously developed "warm-up" scenario, in which organic molecules are formed via diffusive chemistry on surfaces of interstellar grains starting at 30 K, and then released to the gas at higher temperatures during later stages of star formation. In this article, we investigate an alternative scenario in which complex organic species are formed via a sequence of gas-phase reactions between precursor species formed on grain surfaces and then ejected into the gas via efficient reactive desorption, a process in which non-thermal desorption occurs as a result of conversion of the exothermicity of chemical reactions into the ejection of products from the surface. The proposed scenario leads to reasonable if somewhat mixed results at temperatures as low as 10 K and may be considered as a step towards the explanation of abundances of terrestrial-like organic species observed during the earliest stages of star formation.

We present a self-consistent model of a protoplanetary disk: 'ANDES' ('AccretioN disk with Dust Evolution and Sedimentation'). ANDES is based on a flexible and extendable modular structure that includes 1) a 1+1D frequency-dependent continuum radiative transfer module, 2) a module to calculate the chemical evolution using an extended gas-grain network with UV/X-ray-driven processes surface reactions, 3) a module to calculate the gas thermal energy balance, and 4) a 1+1D module that simulates dust grain evolution. For the first time, grain evolution and time-dependent molecular chemistry are included in a protoplanetary disk model. We find that grain growth and sedimentation of large grains to the disk midplane lead to a dust-depleted atmosphere. Consequently, dust and gas temperatures become higher in the inner disk (R < 50 AU) and lower in the outer disk (R > 50 AU), in comparison with the disk model with pristine dust. The response of disk chemical structure to the dust growth and sedimentation is twofold. First, due to higher transparency a partly UV-shielded molecular layer is shifted closer to the dense midplane. Second, the presence of big grains in the disk midplane delays the freeze-out of volatile gas-phase species such as CO there, while in adjacent upper layers the depletion is still effective. Molecular concentrations and thus column densities of many species are enhanced in the disk model with dust evolution, e.g., CO2, NH2CN, HNO, H2O, HCOOH, HCN, CO. We also show that time-dependent chemistry is important for a proper description of gas thermal balance.

The observed gas-phase molecular inventory of hot cores is believed to be significantly impacted by the products of chemistry in interstellar ices. In this study, we report the construction of a full macroscopic Monte Carlo model of both the gas-phase chemistry and the chemistry occurring in the icy mantles of interstellar grains. Our model treats icy grain mantles in a layer-by-layer manner, which incorporates laboratory data on ice desorption correctly. The ice treatment includes a distinction between a reactive ice surface and an inert bulk. The treatment also distinguishes between zeroth and first order desorption, and includes the entrapment of volatile species in more refractory ice mantles. We apply the model to the investigation of the chemistry in hot cores, in which a thick ice mantle built up during the previous cold phase of protostellar evolution undergoes surface reactions and is eventually evaporated. For the first time, the impact of a detailed multilayer approach to grain mantle formation on the warm-up chemistry is explored. The use of a multilayer ice structure has a mixed impact on the abundances of organic species formed during the warm-up phase. For example, the abundance of gaseous HCOOCH3 is lower in the multilayer model than in previous grain models that do not distinguish between layers (so-called "two phase" models). Other gaseous organic species formed in the warm-up phase are affected slightly. Finally, we find that the entrapment of volatile species in water ice can explain the two-jump behavior of H2CO previously found in observations of protostars.

We simulate the chemistry of infrared dark clouds (IRDCs) with a model in which the physical conditions are homogeneous and time-independent. The chemistry is solved as a function of time with three networks: one purely gas-phase, one that includes accretion and desorption, and one, the complete gas-grain network, that includes surface chemistry in addition. We compare our results with observed molecular abundances for two representative IRDCs -- IRDC013.90-1 and IRDC321.73-1 -- using the molecular species N$_2$H$^+$, HC$_3$N, HNC, HCO$^+$, HCN, C$_2$H, NH$_3$ and CS. IRDC013.90-1 is a cold IRDC, with a temperature below 20 K, while IRDC321.73-1 is somewhat warmer, in the range 20 - 30 K. We find that the complete gas-grain model fits the data very well, but that the goodness-of-fit is not sharply peaked at a particular temperature. Surface processes are important for the explanation of the high gas-phase abundance of N$_2$H$^+$ in IRDC321.73-1. The general success of the 0-D model in reproducing single-dish observations of our limited sample of 8 species shows that it is probably sufficient for an explanation of this type of data. To build and justify more complicated models, including spatial temperature and density structure, contraction, and heating, we require high-resolution interferometric data.

We report the first detection of interstellar mercapto radicals, obtained along the sight-line to the submillimeter continuum source W49N. We have used the GREAT instrument on SOFIA to observe the 1383 GHz Doublet Pi 3/2 J = 5/2 - 3/2 lambda doublet in the upper sideband of the L1 receiver. The resultant spectrum reveals SH absorption in material local to W49N, as well as in foreground gas, unassociated with W49N, that is located along the sight-line. For the foreground material at velocities in the range 37 - 44 km/s with respect to the local standard of rest, we infer a total SH column density ~ 2.6 E+12 cm-2, corresponding to an abundance of ~ 7 E-9 relative to H2, and yielding an SH/H2S abundance ratio ~ 0.13. The observed SH/H2S abundance ratio is much smaller than that predicted by standard models for the production of SH and H2S in turbulent dissipation regions and shocks, and suggests that the endothermic neutral-neutral reaction SH + H2 -> H2S + H must be enhanced along with the ion-neutral reactions believed to produce CH+ and SH+ in diffuse molecular clouds.

We study details of the UV radiation transfer in a protoplanetary disk, paying attention to the influence of dust growth and sedimentation on the disk density and temperature. Also, we show how the dust evolution affects photoreaction rates of key molecules, like CN and CS.

We study the impact of dust evolution in a protoplanetary disk around a T Tauri star on the disk chemical composition. For the first time we utilize a comprehensive model of dust evolution which includes growth, fragmentation and sedimentation. Specific attention is paid to the influence of grain evolution on the penetration of the UV field in the disk. A chemical model that includes a comprehensive set of gas phase and grain surface chemical reactions is used to simulate the chemical structure of the disk. The main effect of the grain evolution on the disk chemical composition comes from sedimentation, and, to a lesser degree, from the reduction of the total grain surface area. The net effect of grain growth is suppressed by the fragmentation process which maintains a population of small grains, dominating the total grain surface area. We consider three models of dust properties. In model GS both growth and sedimentation are taken into account. In models A5 and A4 all grains are assumed to have the same size (10(-5) cm and 10(-4) cm, respectively) with constant gas-to-dust mass ratio of 100. Like in previous studies, the "three-layer" pattern (midplane, molecular layer, hot atmosphere) in the disk chemical structure is preserved in all models, but shifted closer to the midplane in models with increased grain size (GS and A4). Unlike other similar studies, we find that in models GS and A4 column densities of most gas-phase species are enhanced by 1-3 orders of magnitude relative to those in a model with pristine dust (A5), while column densities of their surface counterparts are decreased. We show that column densities of certain species, like C2H, HC(2n+1)N (n=0-3), H2O and some other molecules, as well as the C2H2/HCN abundance ratio which are accessible with Herschel and ALMA can be used as observational tracers of early stages of the grain evolution process in protoplanetary disks.

In this study we demonstrate for the first time that the unified Monte Carlo approach can be applied to model gas-grain chemistry in large reaction networks. Specifically, we build a time-dependent gas-grain chemical model of the interstellar medium, involving about 6000 gas-phase and 200 grain surface reactions. This model is used to test the validity of the standard and modified rate equation methods in models of dense and translucent molecular clouds and to specify under which conditions the use of the stochastic approach is desirable. We found that at temperatures 25--30 K gas-phase abundances of H$_2$O, NH$_3$, CO and many other gas-phase and surface species in the stochastic model differ from those in the deterministic models by more than an order of magnitude, at least, when tunneling is accounted for and/or diffusion energies are 3x lower than the binding energies. In this case, surface reactions, involving light species, proceed faster than accretion of the same species. In contrast, in the model without tunneling and with high binding energies, when the typical timescale of a surface recombination is greater than the timescale of accretion onto the grain, we obtain almost perfect agreement between results of Monte Carlo and deterministic calculations in the same temperature range. At lower temperatures ($\sim10$ K) gaseous and, in particular, surface abundances of most important molecules are not much affected by stochastic processes.

We study how uncertainties in the rate coefficients of chemical reactions in the RATE06 database affect abundances and column densities of key molecules in protoplanetary disks. We randomly varied the gas-phase reaction rates within their uncertainty limits and calculated the time-dependent abundances and column densities using a gas-grain chemical model and a flaring steady-state disk model. We find that key species can be separated into two distinct groups according to the sensitivity of their column densities to the rate uncertainties. The first group includes CO, C$^+$, H$_3^+$, H$_2$O, NH$_3$, N$_2$H$^+$, and HCNH$^+$. For these species, the column densities are not very sensitive to the rate uncertainties but the abundances in specific regions are. The second group includes CS, CO$_2$, HCO$^+$, H$_2$CO, C$_2$H, CN, HCN, HNC and other, more complex species, for which high abundances and abundance uncertainties co-exist in the same disk region, leading to larger scatters in the column densities. However, even for complex and heavy molecules, the dispersion in their column densities is not more than a factor of ~4. We perform a sensitivity analysis of the computed abundances to rate uncertainties and identify those reactions with the most problematic rate coefficients. We conclude that the rate coefficients of about a hundred of chemical reactions need to be determined more accurately in order to greatly improve the reliability of modern astrochemical models. This improvement should be an ultimate goal for future laboratory studies and theoretical investigations.

With the chemical reaction rate database UMIST95 (Millar et al. 1997) we analyze how uncertainties in rate constants of gas-phase chemical reactions influence the modelling of molecular abundances in the interstellar medium. Random variations are introduced into the rate constants to estimate the scatter in theoretical abundances. Calculations are performed for dark and translucent molecular clouds where gas phase chemistry is adequate. Similar approach was used by Pineau des Forets & Roueff (2000) for the study of chemical bistability. All the species are divided into 6 sensitivity groups according to the value of the scatter in their model abundances computed with varied rate constants. It is shown that the distribution of species within these groups depends on the number of atoms in a molecule and on the adopted physical conditions. The simple method is suggested which allows to single out reactions that are most important for the evolution of a given species.