Complex organic molecules with three carbon atoms are found in the earliest stages of star formation. In particular, propenal (C$_2$H$_3$CHO) is a species of interest due to its implication in the formation of more complex species and even biotic molecules. This study aims to search for the presence of C$_2$H$_3$CHO and other three-carbon species such as propylene (C$_3$H$_6$) in the hot corino region of the low-mass protostellar binary IRAS 16293--2422 to understand their formation pathways. We use ALMA observations in Band 6 and 7 from various surveys to search for the presence of C$_3$H$_6$ and C$_2$H$_3$CHO towards the protostar IRAS 16293--2422 B (IRAS 16293B). We report the detection of both C$_3$H$_6$ and C$_2$H$_3$CHO towards IRAS 16293B, however, no unblended lines were found towards the other component of the binary system, IRAS 16293A. We derive column density upper limits for C$_3$H$_8$, HCCCHO, n-C$_3$H$_7$OH, i-C$_3$H$_7$OH, C$_3$O, and cis-HC(O)CHO towards IRAS 16293B. We then use a three-phase chemical model to simulate the formation of these species in a typical prestellar environment followed by its hydrodynamical collapse until the birth of the central protostar. Different formation paths, such as successive hydrogenation and radical-radical additions on grain surfaces, are tested and compared to the observational results. The simulations reproduce the abundances within one order of magnitude from those observed towards IRAS 16293B, with the best agreement found for a rate of $10^{-12}$ cm$^3$ s$^{-1}$ for the gas-phase reaction C$_3$ + O $\rightarrow$ C$_2$ + CO. Successive hydrogenations of C$_3$, HC(O)CHO, and CH$_3$OCHO on grain surfaces are a major and crucial formation route of complex organics molecules, whereas both successive hydrogenation pathways and radical-radical addition reactions contribute to the formation of C$_2$H$_5$CHO.

Methyl carbamate CH$_3$OC(O)NH$_2$ is an isomer of glycine. Quantum chemical analyses show that methyl carbamate is more stable isomer than glycine. Because of this, there could be a higher chance for methyl carbamte to exist in the interstellar medium as compared to glycine. Despite immense searches, till now glycine has not been detected in the ISM, therefore it is worthwhile to search its isomer methyl carbamate. In this paper, we present the constraints of methyl carbamate formation under the interstellar conditions. Large complex organic molecules are favorably produced in hot-corino environments of low mass protostars. We for the first time carried out astrochemical modeling focusing on the formation of methyl carbamate in physical conditions similar to hot-corino objects. Consequently, we examined ALMA archival data for existing spectral line observations toward hot corinos NGC1333 IRAS 4A2 and IRAS 16293B. Within the common spectral range towards these sources, we found three features are possibly related to the spectral transitions of methyl carbamate and consequently estimate the upper limit of column densities. Results of chemical modeling are consistent with the observational upper limit of estimated column density/abundance toward the sources. This may hint the validation of the proposed formation mechanism. Future observations using telescope like ngVLA may confirm the presence of MC toward the hot corinos.

Sulphur is one of the most abundant elements in the Universe. Surprisingly, sulphuretted molecules are not as abundant as expected in the interstellar medium, and the identity of the main sulphur reservoir is still an open question. Our goal is to investigate the H$_{2}$S chemistry in dark clouds, as this stable molecule is a potential sulphur reservoir. Using millimeter observations of CS, SO, H$_{2}$S, and their isotopologues, we determine the physical conditions and H$_{2}$S abundances along the cores TMC 1-C, TMC 1-CP, and Barnard 1b. The gas-grain model Nautilus is then used to model the sulphur chemistry and explore the impact of photo-desorption and chemical desorption on the H$_2$S abundance. Our model shows that chemical desorption is the main source of gas-phase H$_2$S in dark cores. The measured H$_{2}$S abundance can only be fitted if we assume that the chemical desorption rate decreases by more than a factor of 10 when $n_{\rm H}>2\times10^{4}$. This change in the desorption rate is consistent with the formation of thick H$_2$O and CO ice mantles on grain surfaces. The observed SO and H$_2$S abundances are in good agreement with our predictions adopting an undepleted value of the sulphur abundance. However, the CS abundance is overestimated by a factor of $5-10$. Along the three cores, atomic S is predicted to be the main sulphur reservoir. We conclude that the gaseous H$_2$S abundance is well reproduced, assuming undepleted sulphur abundance and chemical desorption as the main source of H$_2$S. The behavior of the observed H$_{2}$S abundance suggests a changing desorption efficiency, which would probe the snowline in these cores. Our model, however, overestimates the observed gas-phase CS abundance. Given the uncertainty in the sulphur chemistry, our data are consistent with a cosmic elemental S abundance with an uncertainty of a factor of 10.

Aims. Our goal is to complete the inventory of S-bearing molecules and their abundances in the prototypical photodissociation region (PDR) the Horsehead nebula to gain insight into sulphur chemistry in UV irradiated regions. Based on the WHISPER millimeter (mm) line survey, our goal is to provide an improved and more accurate description of sulphur species and their abundances towards the core and PDR positions in the Horsehead. Methods. The Monte Carlo Markov chain (MCMC) methodology and the molecular excitation and radiative transfer code RADEX were used to explore the parameter space and determine physical conditions and beam-averaged molecular abundances. Results. A total of 13 S-bearing species (CS, SO, SO2, OCS, H2CS - both ortho and para - HDCS, C2S, HCS+, SO+, H2S, S2H, NS and NS+) have been detected in the two targeted positions. This is the first detection of SO+ in the Horsehead and the first detection of NS+ in any PDR. We find a differentiated chemical behaviour between C-S and O-S bearing species within the nebula. The C-S bearing species C2S and o-H2CS present fractional abundances a factor grater than two higher in the core than in the PDR. In contrast, the O-S bearing molecules SO, SO2, and OCS present similar abundances towards both positions. A few molecules, SO+, NS, and NS+, are more abundant towards the PDR than towards the core, and could be considered as PDR tracers. Conclusions. This is the first complete study of S-bearing species towards a PDR. Our study shows that CS, SO, and H2S are the most abundant S-bearing molecules in the PDR with abundances of a few 1E-9. We recall that SH, SH+, S, and S+ are not observable at the wavelengths covered by the WHISPER survey. At the spatial scale of our observations, the total abundance of S atoms locked in the detected species is < 1E-8, only ~0.1% of the cosmic sulphur abundance.

Molecular oxygen has been the subject of many observational searches as chemical models predicted it to be a reservoir of oxygen. Although it has been detected in two regions of the interstellar medium, its rarity is a challenge for astrochemical models. In this paper, we have combined the physical conditions computed with smoothed particle hydrodynamics (SPH) simulations with our full gas-grain chemical model Nautilus, to study the predicted O2 abundance in interstellar material forming cold cores. We thus follow the chemical evolution of gas and ices in parcels of material from the diffuse interstellar conditions to the cold dense cores. Most of our predicted O2 abundances are below 1e-8 (with respect to the total proton density) and the predicted column densities in simulated cold cores is at maximum a few 1e14 cm-2, in agreement with the non detection limits. This low O2 abundance can be explained by the fact that, in a large fraction of the interstellar material, the atomic oxygen is depleted onto the grain surface (and hydrogenated to form H2O) before O2 can be formed in the gas-phase and protected from UV photo-dissociations. We could achieve this result only because we took into account the full history of the evolution of the physical conditions from the diffuse medium to the cold cores.

The complexity of physico-chemical models of star formation is increasing, with models that take into account new processes and more realistic setups. These models allow astrochemists to compute the evolution of chemical species throughout star formation. Hence, comparing the outputs of such models to observations allows to bring new constraints on star formation. The work presented in this paper is based on the recent public release of a database of radiation hydrodynamical low-mass star formation models. We used this database as physical parameters to compute the time dependent chemical composition of collapsing cores with a 3-phase gas-grain model. The results are analyzed to find chemical tracers of the initial physical parameters of collapse such as the mass, radius, temperature, density, and free-fall time. They are also compared to observed molecular abundances of Class 0 protostars. We find numerous tracers of the initial parameters of collapse, except for the initial mass. More particularly, we find that gas phase CH3CN, NS and OCS trace the initial temperature, while H2CS trace the initial density and free-fall time of the parent cloud. The comparison of our results with a sample of 12 Class 0 low mass protostars allows us to constrain the initial parameters of collapse of low-mass prestellar cores. We find that low-mass protostars are preferentially formed within large cores with radii greater than 20000 au, masses between 2 and 4 Msol, temperatures lower or equal to 15 K, and densities between 6e4 and 2.5e5 part.cm-3, corresponding to free-fall times between 100 and 200 kyrs.

Cyanogen (NCCN) is the simplest member of the dicyanopolyynes group, and has been proposed as a major source of the CN radical observed in cometary atmospheres. Although not detected through its rotational spectrum in the cold interstellar medium, this very stable species is supposed to be very abundant. The chemistry of cyanogen in the cold interstellar medium can be investigated through its metastable isomer, CNCN (isocyanogen). Its formation may provide a clue on the widely abundant CN radical observed in cometary atmospheres. We performed an unbiased spectral survey of the L1544 proto-typical prestellar core, using the IRAM-30m and have analysed, for this paper, the nitrogen chemistry that leads to the formation of isocyanogen. We report on the first detection of CNCN, NCCNH+, C3N, CH3CN, C2H3CN, and H2CN in L1544. We built a detailed chemical network for NCCN/CNCN/HC2N2+ involving all the nitrogen bearing species detected (CN, HCN, HNC, C3N, CNCN, CH3CN, CH2CN, HCCNC, HC3N, HNC3, H2CN, C2H3CN, HCNH+, HC3NH+) and the upper limits on C4N, C2N. The main cyanogen production pathways considered in the network are the CN + HNC and N + C3N reactions. The comparison between the observations of the nitrogen bearing species and the predictions from the chemical modelling shows a very good agreement, taking into account the new chemical network. The expected cyanogen abundance is greater than the isocyanogen abundance by a factor of 100. Although cyanogen cannot be detected through its rotational spectrum, the chemical modelling predicts that it should be abundant in the gas phase and hence might be traced through the detection of isocyanogen. It is however expected to have a very low abundance on the grain surfaces compared to HCN.

Nitrogen oxides are thought to play a significant role as a nitrogen reservoir and to potentially participate in the formation of more complex species. Until now, only NO, N$_2$O and HNO have been detected in the interstellar medium. We report the first interstellar detection of nitrous acid (HONO). Twelve lines were identified towards component B of the low-mass protostellar binary IRAS~16293--2422 with the Atacama Large Millimeter/submillimeter Array, at the position where NO and N$_2$O have previously been seen. A local thermodynamic equilibrium model was used to derive the column density ($\sim$ 9 $\times$ 10$^{14}$ cm$^{-2}$ in a 0.5'' beam) and excitation temperature ($\sim$ 100 K) of this molecule. HNO, NO$_2$, NO$^+$, and HNO$_3$ were also searched for in the data, but not detected. We simulated the HONO formation using an updated version of the chemical code Nautilus and compared the results with the observations. The chemical model is able to reproduce satisfactorily the HONO, N$_2$O, and NO$_2$ abundances, but not the NO, HNO, and NH$_2$OH abundances. This could be due to some thermal desorption mechanisms being destructive and therefore limiting the amount of HNO and NH$_2$OH present in the gas phase. Other options are UV photodestruction of these species in ices or missing reactions potentially relevant at protostellar temperatures.

We have developed the first gas-grain chemical model for oxygen fractionation (also including sulphur fractionation) in dense molecular clouds, demonstrating that gas-phase chemistry generates variable oxygen fractionation levels, with a particularly strong effect for NO, SO, O2, and SO2. This large effect is due to the efficiency of the neutral 18O + NO, 18O + SO, and 18O + O2 exchange reactions. The modeling results were compared to new and existing observed isotopic ratios in a selection of cold cores. The good agreement between model and observations requires that the gas-phase abundance of neutral oxygen atoms is large in the observed regions. The S16O/S18O ratio is predicted to vary substantially over time showing that it can be used as a sensitive chemical proxy for matter evolution in dense molecular clouds.

Protoplanetary disks are challenging objects for astrochemical models due to strong density and temperature gradients and due to the UV photons 2D propagation. In this paper, we have studied the importance of several model parameters on the predicted column densities of observed species. We considered: 1) 2-phase (gas and homogeneous grains) or 3-phase (gas, surface, and bulk of grains) models, 2) several initial compositions, 3) grain growth and dust settling, and 4) several cosmic-ray ionization rates. Our main result is that dust settling is the most crucial parameter. Including this effect renders the computed column densities sensitive to all the other model parameters, except cosmic-ray ionization rate. In fact, we found almost no effect of this parameter for radii larger than 10 au (the minimum radius studied here) except for N2H+. We also compared all our models with all the column densities observed in the protoplanetary disk around DM Tau and were not able to reproduce all the observations despite the studied parameters. N2H+ seems to be the most sensitive species. Its observation in protoplanetary disks at large radius could indicate enough N2 in the gas-phase (inhibited by the 3-phase model, but boosted by the settling) and a low electron abundance (favored by low C and S elemental abundances).

Methyl cyanide (CH3CN) and propyne (CH3CCH) are two molecules commonly used as gas thermometers for interstellar gas. They are detected in several astrophysical environments and in particular towards protostars. Using data of the low-mass protostar IRAS 16293-2422 obtained with the IRAM 30m single-dish telescope, we constrained the origin of these two molecules in the envelope of the source. The line shape comparison and the results of a radiative transfer analysis both indicate that the emission of CH3CN arises from a warmer and inner region of the envelope than the CH3CCH emission. We compare the observational results with the predictions of a gas-grain chemical model. Our model predicts a peak abundance of CH3CCH in the gas-phase in the outer part of the envelope, at around 2000 au from the central star, which is relatively close to the emission size derived from the observations. The predicted CH3CN abundance only rises at the radius where the grain mantle ices evaporate, with an abundance similar to the one derived from the observations.

Context: Diffusion of species on the dust surface is a key process for determining the chemical composition of interstellar ices. On the dust surface, adsorbed species diffuse from one potential well to another and react with other adsorbed reactants, resulting in the formation of simple and complex molecules. Aims: We study the impact on the abundances of the species simulated by the chemical codes by considering the uncertainties in the diffusion energy of adsorbed species. We aim to limit the uncertainties in the abundances as calculated by chemical codes by identifying the surface species that result in a larger error because of the uncertainties in their diffusion energy. Methods: We ran various cases with 2000 to 10000 simulations in each case and varied the diffusion energies of some or all surface species randomly. We calculated Pearson correlation coefficients between the abundances and the ratio of diffusion to binding energy of adsorbed species. We identified the species that introduce maximum uncertainty in the ice and gas-phase abundances. With these species we ran three sets, with 2000 simulations in each, to quantify the uncertainties they introduce. Results: We present the abundances of various molecules in the gas phase and also on the dust surface at different time intervals during the simulation. We show which species produce a large uncertainty in the abundances. We sorted species into different groups in accordance with their importance in propagating uncertainty in the chemical network. Conclusions: We show that CO, H$_2$, O, N, and CH$_3$ are the key species for uncertainties in the abundances, while CH$_2$, HCO, S and O$_2$ come next, followed by NO, HS, and CH. We also show that by limiting the uncertainties in the ratio of diffusion to binding energy of these species, we can eliminate the uncertainties in the gas-phase abundances of almost all the species.

Studying molecular species in protoplanetary disks is very useful to characterize the properties of these objects, which are the site of planet formation. We attempt to constrain the chemistry of S-bearing molecules in the cold parts of circumstellar disk of GG Tau A. We searched for H$_2$S, CS, SO, and SO$_2$ in the dense disk around GG Tau A with the NOrthem Extended Millimeter Array (NOEMA) interferometer. We detected H$_2$S emission from the dense and cold ring orbiting around GG Tau A. This is the first detection of H$_2$S in a protoplanetary disk. We also detected HCO$^+$, H$^{13}$CO$^+$, and DCO$^+$ in the disk. Upper limits for other molecules, CCS, SO$_2$, SO, HC$_3$N, and $c$-C$_3$H$_2$ are also obtained. The observed DCO$^+$/HCO$^+$ ratio is similar to those in other disks. The observed column densities, derived using our radiative transfer code DiskFit, are then compared with those from our chemical code Nautilus. The column densities are in reasonable agreement for DCO$^{+}$, CS, CCS, and SO$_2$. For H$_2$S and SO, our predicted vertical integrated column densities are more than a factor of 10 higher than the measured values. Our results reinforce the hypothesis that only a strong sulfur depletion may explain the low observed H$_2$S column density in the disk. The H$_2$S detection in GG Tau A is most likely linked to the much larger mass of this disk compared to that in other T Tauri systems.

The study of the chemical evolution of glycine in the interstellar medium is one of challenging topics in astrochemistry. Here, we present the chemical modeling of glycine in hot cores using the state-of-the-art three-phase chemical model NAUTILUS, which is focused on the latest glycine chemistry. For the formation process of glycine on the grain surface, we obtained consistent results with previous studies that glycine would be formed via the reactions of COOH with CH$_2$NH$_2$. However, we will report three important findings regarding the chemical evolution and the detectability of interstellar glycine. First, with the experimentally obtained binding energy from the temperature programmed thermal desorption (TPD) experiment, a large proportion of glycine was destroyed through the grain surface reactions with NH or CH$_3$O radicals before it fully evaporates. As a result, the formation process in the gas phase is more important than thermal evaporation from grains. If this is the case, NH$_2$OH and CH$_3$COOH rather than CH$_3$NH$_2$ and CH$_2$NH would be the essential precursors to the gas phase glycine. Secondly, since the gas phase glycine will be quickly destroyed by positive ions or radicals, early evolutionary phase of the hot cores would be the preferable target for the future glycine surveys. Thirdly, we suggest the possibility that the suprathermal hydrogen atoms can strongly accelerate the formation of COOH radicals from CO$_2$, resulting in the dramatic increase of formation rate of glycine on grains. The efficiency of this process should be investigated in detail by theoretical and experimental studies in the future.

The different spatial distributions of N-bearing and O-bearing species, as is well known towards Orion~KL, is one of the long-lasting mysteries. We conducted a survey observation and chemical modeling study to investigate if the different distributions of O- and N-bearing species are widely recognized in general star-forming regions. First, we report our observational results of complex organic molecules (COMs) with the 45~m radio telescope at the Nobeyama Radio Observatory towards eight star-forming regions. Through our spectral survey ranging from 80 to 108~GHz, we detected CH$_3$OH, HCOOCH$_3$, CH$_3$OCH$_3$, (CH$_3$)$_2$CO, CH$_3$CHO, CH$_3$CH$_2$CN, CH$_2$CHCN, and NH$_2$CHO. Their molecular abundances were derived via the rotation diagram and the least squares methods. We found that N-bearing molecules, tend to show stronger correlations with other N-bearing molecules rather than O-bearing molecules. While G10.47+0.03 showed high fractional abundances of N-bearing species, those in NGC6334F were not so rich, being less than 0.01 compared to CH$_3$OH. Then, the molecular abundances towards these sources were evaluated by chemical modeling with NAUTILUS three-phase gas-grain chemical code. Through the simulations of time evolutions for the abundances of COMs, we suggest that observed correlations of fractional abundances between COMs can be explained by the combination of the different temperature structures inside the hot cores and the different evolutionary phase. Since our modeling could not fully explain the observed excitation temperatures, it is important to investigate the efficiency of grain surface reactions and their activation barriers, and the binding energy of COMs to further promote our understanding.

The L1544 pre-stellar core has been observed as part of the ASAI IRAM 30m Large Program as well as follow-up programs. These observations have revealed the chemical richness of the earliest phases of low-mass star-forming regions. In this paper we focus on the twenty-one sulphur bearing species (ions, isotopomers and deuteration) that have been detected in this spectral-survey through fifty one transitions: CS, CCS, C3S, SO, SO2, H2CS, OCS, HSCN, NS, HCS+, NS+ and H2S. We also report the tentative detection (4 {\sigma} level) for methyl mercaptan (CH3SH). LTE and non-LTE radiative transfer modelling have been performed and we used the nautilus chemical code updated with the most recent chemical network for sulphur to explain our observations. From the chemical modelling we expect a strong radial variation for the abundances of these species, which mostly are emitted in the external layer where non thermal desorption of other species has previously been observed. We show that the chemical study cannot be compared to what has been done for the TMC-1 dark cloud, where the abundance is supposed constant along the line of sight, and conclude that a strong sulphur depletion is necessary to fully reproduce our observations of the prototypical pre-stellar core L1544.

Ions and electrons play an important role in various stages of the star formation process. By following the magnetic field of their environment and interacting with neutral species, they slow down the gravitational collapse of the proto-star envelope. This process (known as ambipolar diffusion) depends on the ionisation degree, which can be derived from the \hco abundance. We present a study of \hco and its isotopologues (H$^{13}$CO$^+$, HC$^{18}$O$^+$, DCO$^+$, and H$^{13}$CO$^+$) in the low-mass proto-star IRAS16293$-$2422. The structure of this object is complex, and the HCO$^+$ emission arises from the contribution of a young NW-SE outflow, the proto-stellar envelope and the foreground cloud. We aim at constraining the physical parameters of these structures using all the observed transitions. For the young NW-SE outflow, we derive $T_{\rm kin}=180-220$ K and $n({\rm H_2})=(4-7)\times10^6$ cm$^{-3}$ with an HCO$^+$ abundance of $(3-5)\times10^{-9}$. Following previous studies, we demonstrate that the presence of a cold ($T_{\rm kin}$$\leqslant$30 K) and low density ($n({\rm H_2})\leqslant1\times10^4$ cm$^{-3}$) foreground cloud is also necessary to reproduce the observed line profiles. We have used the gas-grain chemical code \textsc{nautilus} to derive the HCO$^+$ abundance profile across the envelope and the external regions where X(HCO$^+$)$\gtrsim1\times10^{-9}$ dominate the envelope emission. From this, we derive an ionisation degree of $10^{-8.9}\,\lesssim\,x(e)\,\lesssim\,10^{-7.9}$. The ambipolar diffusion timescale is $\sim$5 times the free-fall timescale, indicating that the magnetic field starts to support the source against gravitational collapse and the magnetic field strength is estimated to be $6-46 \mu$G.

The protonated form of CO2, HOCO+, is assumed to be an indirect tracer of CO2 in the millimeter/submillimeter regime since CO2 lacks a permanent dipole moment. Here, we report the detection of two rotational emission lines (4 0,4-3 0,3) and (5 0,5-4 0,4) of HOCO+ in IRAS 16293-2422. For our observations, we have used EMIR heterodyne 3 mm receiver of the IRAM 30m telescope. The observed abundance of HOCO+ is compared with the simulations using the 3-phase NAUTILUS chemical model. Implications of the measured abundances of HOCO+ to study the chemistry of CO2 ices using JWST-MIRI and NIRSpec are discussed as well.

Species abundances in the interstellar medium (ISM) strongly depend on the chemistry occurring at the surfaces of the dust grains. To describe the complexity of the chemistry, various numerical models have been constructed. In most of these models, the grains are described by a single size of 0.1$\mu$m. We study the impact on the abundances of many species observed in the cold cores by considering several grain sizes in the Nautilus multi-grain model. We used grain sizes with radii in the range of $0.005\mu$m to $0.25\mu$m. We sampled this range in many bins. We used the previously published, MRN and WD grain size distributions to calculate the number density of grains in each bin. Other parameters such as the grain surface temperature or the cosmic-ray-induced desorption rates also vary with grain sizes. We present the abundances of various molecules in the gas phase and also on the dust surface at different time intervals during the simulation. We present a comparative study of results obtained using the single grain and the multi-grain models. We also compare our results with the observed abundances in TMC-1 and L134N clouds. We show that the grain size, the grain size dependent surface temperature and the peak surface temperature induced by cosmic ray collisions, play key roles in determining the ice and the gas phase abundances of various molecules. We also show that the differences between the MRN and the WD models are crucial for better fitting the observed abundances in different regions in the ISM. We show that the small grains play a very important role in the enrichment of the gas phase with the species which are mainly formed on the grain surface, as non-thermal desorption induced by collisions of cosmic ray particles is very efficient on the small grains.

Aims. We study the effect of large scale dynamics on the molecular composition of the dense interstellar medium during the transition between diffuse to dense clouds. Methods. We followed the formation of dense clouds (on sub-parsec scales) through the dynamics of the interstellar medium at galac- tic scales. We used results from smoothed particle hydrodynamics (SPH) simulations from which we extracted physical parameters that are used as inputs for our full gas-grain chemical model. In these simulations, the evolution of the interstellar matter is followed for ~50 Myr. The warm low-density interstellar medium gas flows into spiral arms where orbit crowding produces the shock formation of dense clouds, which are held together temporarily by the external pressure. Results. We show that depending on the physical history of each SPH particle, the molecular composition of the modeled dense clouds presents a high dispersion in the computed abundances even if the local physical properties are similar. We find that carbon chains are the most affected species and show that these differences are directly connected to differences in (1) the electronic fraction, (2) the C/O ratio, and (3) the local physical conditions. We argue that differences in the dynamical evolution of the gas that formed dense clouds could account for the molecular diversity observed between and within these clouds. Conclusions. This study shows the importance of past physical conditions in establishing the chemical composition of the dense medium.

We present the first detection of gas phase S2H in the Horsehead, a moderately UV-irradiated nebula. This confirms the presence of doubly sulfuretted species in the interstellar medium and opens a new challenge for sulfur chemistry. The observed S2H abundance is ~5x10$^{-11}$, only a factor 4-6 lower than that of the widespread H2S molecule. H2S and S2H are efficiently formed on the UV-irradiated icy grain mantles. We performed ice irradiation experiments to determine the H2S and S2H photodesorption yields. The obtained values are ~1.2x10$^{-3}$ and <1x10$^{-5}$ molecules per incident photon for H2S and S2H, respectively. Our upper limit to the S2H photodesorption yield suggests that photo-desorption is not a competitive mechanism to release the S2H molecules to the gas phase. Other desorption mechanisms such as chemical desorption, cosmic-ray desorption and grain shattering can increase the gaseous S2H abundance to some extent. Alternatively, S2H can be formed via gas phase reactions involving gaseous H2S and the abundant ions S+ and SH+. The detection of S2H in this nebula could be therefore the result of the coexistence of an active grain surface chemistry and gaseous photo-chemistry.

Sulphur-bearing species are often used to probe the evolution of hot cores since their abundances are particularly sensitive to physical and chemical variations. However, the chemistry of sulphur is not well understood in these regions, notably because observations of several hot cores displayed a large variety of sulphur compositions, and because the reservoir of sulphur in dense clouds, in which hot cores form, is still poorly constrained. In order to help disentangled its complexity, we present a fresh comprehensive review of sulphur chemistry in hot cores along with a study of its sensibility to temperature and pre-collapse chemical composition. In parallel, we analyse the discrepencies that result from the use of two different types of models (static and dynamic) to highlight the sensitivity to the choice of model to be used in astrochemical studies. Our results show that the pre-collapse chemical composition is a critical parameter for sulphur chemistry in hot cores and could explain the different sulphur compositions observed. We also report that differences in abundances for a given species between the static and dynamic models can reach six orders of magnitude in the hot core, which reveals the key role of the choice of model in astrochemical studies.

Molecular hydrogen is the most abundant molecule in the universe. It is the first one to form and survive photo-dissociation in tenuous environments. Its formation involves catalytic reactions on the surface of interstellar grains. The micro-physics of the formation process has been investigated intensively in the last 20 years, in parallel of new astrophysical observational and modeling progresses. In the perspectives of the probable revolution brought by the future satellite JWST, this article has been written to present what we think we know about the H$_2$ formation in a variety of interstellar environments.

Methyl isocyanate (CH3NCO) is one of the important complex organic molecules detected on the comet 67P/Churyumov-Gerasimenko by Rosetta's Philae lander. It was also detected in hot cores around high-mass protostars along with a recent detection in the solar-type protostar IRAS 16293-2422. We propose here a gas-grain chemical model to form CH3NCO after reviewing various formation pathways with quantum chemical computations. We have used NAUTILUS 3-phase gas-grain chemical model to compare observed abundances in the IRAS 16293-2422. Our chemical model clearly indicates the ice phase origin of CH3NCO.

We report the detection of linear and cyclic isomers of C3H and C3H2 towards various starless cores and review the corresponding chemical pathways involving neutral (C3Hx with x=1,2) and ionic (C3Hx+ with x = 1,2,3) isomers. We highlight the role of the branching ratio of electronic Dissociative Recombination (DR) reactions of C3H2+ and C3H3+ isomers showing that the statistical treatment of the relaxation of C3H* and C3H2* produced in these DR reactions may explain the relative c,l-C3H and c,l-C3H2 abundances. We have also introduced in the model the third isomer of C3H2 (HCCCH). The observed cyclic-to-linear C3H2 ratio vary from 110 + or - 30 for molecular clouds with a total density around 1e4 molecules.cm-3 to 30 + or - 10 for molecular clouds with a total density around 4e5 molecules.cm-3, a trend well reproduced with our updated model. The higher ratio for low molecular cloud densities is mainly determined by the importance of the H + l-C3H2 -> H + c-C3H2 and H + t-C3H2 -> H + c-C3H2 isomerization reactions.

A wide variety of molecules have recently been detected in the Horsehead nebula photodissociation region (PDR) suggesting that: (i) gas-phase and grain chemistries should both contribute to the formation of organic molecules, and (ii) far-ultraviolet (FUV) photodesorption may explain the release into the gas phase of grain surface species. In order to tackle these specific problems and more generally in order to better constrain the chemical structure of these types of environments we present a study of the Horsehead nebula gas-grain chemistry. To do so we used the 1D astrochemical gas-grain code Nautilus with an appropriate physical structure computed with the Meudon PDR Code and compared our modeled outcomes with published observations and with previously modeled results when available. The use of a large set of chemical reactions coupled with the time-dependent code Nautilus allows us to reproduce most of the observations well, including those of the first detections in a PDR of the organic molecules HCOOH, CH2CO, CH3CHO and CH3CCH, which are mostly associated with hot cores. We also provide some abundance predictions for other molecules of interest. Understanding the chemistry behind the detection of these organic molecules is crucial to better constrain the environments these molecules can probe.

Determining the gas density and temperature structures of protoplanetary disks is a fundamental task to constrain planet formation theories. This is a challenging procedure and most determinations are based on model-dependent assumptions. We attempt a direct determination of the radial and vertical temperature structure of the Flying Saucer disk, thanks to its favorable inclination of 90 degrees. We present a method based on the tomographic study of an edge-on disk. Using ALMA, we observe at 0.5$"$ resolution the Flying Saucer in CO J=2-1 and CS J=5-4. This edge-on disk appears in silhouette against the CO J=2-1 emission from background molecular clouds in $\rho$ Oph. The combination of velocity gradients due to the Keplerian rotation of the disk and intensity variations in the CO background as a function of velocity provide a direct measure of the gas temperature as a function of radius and height above the disk mid-plane. The overall thermal structure is consistent with model predictions, with a cold ($< 15-12 $~K), CO-depleted mid-plane, and a warmer disk atmosphere. However, we find evidence for CO gas along the mid-plane beyond a radius of about 200\,au, coincident with a change of grain properties. Such a behavior is expected in case of efficient rise of UV penetration re-heating the disk and thus allowing CO thermal desorption or favoring direct CO photo-desorption. CO is also detected up to 3-4 scale heights while CS is confined around 1 scale height above the mid-plane. The limits of the method due to finite spatial and spectral resolutions are also discussed. This method appears to be very promising to determine the gas structure of planet-forming disks, provided that the molecular data have an angular resolution which is high enough, of the order of $0.3 - 0.1"$ at the distance of the nearest star forming regions.

Sulphur-bearing species are often used to probe the physical structure of star forming regions of the interstellar medium, but the chemistry of sulphur in these regions is still poorly understood. In dark clouds, sulphur is supposed to be depleted under a form which is still unknown despite numerous observations and chemical modeling studies that have been performed. In order to improve the modeling of sulphur chemistry, we propose an enhancement of the sulphur chemical network using experimental and theoretical literature. We test the effect of the updated network on the outputs of a three phases gas-grain chemical model for dark cloud conditions using different elemental sulphur abundances. More particularly, we focus our study on the main sulphur reservoirs as well as on the agreement between model predictions and the abundances observed in the dark cloud TMC-1 (CP). Our results show that depending on the age of the observed cloud, the reservoir of sulphur could either be atomic sulphur in the gas-phase or HS and H2S in icy grain bulks. We also report the first chemical model able to reproduce the abundances of observed S-bearing species in TMC-1 (CP) using as elemental abundance of sulphur its cosmic value.

Cyclic-C3HD (c-C3HD) is a singly deuterated isotopologue of c-C3H2, which is one of the most abundant and widespread molecules in our Galaxy. We observed IRAS 16293-2422 in the 3 mm band with a single frequency setup using the EMIR heterodyne 3 mm receiver of the IRAM 30m telescope. We observed seven lines of c-C3HD and three lines of c-C3H2. Observed abundances are compared with astrochemical simulations using the NAUTILUS gas-grain chemical model. Our results clearly show that c-C3HD can be used as an important supplement for studying chemistry and physical conditions for cold environments. Assuming that the size of the protostellar envelope is 3000 AU and same excitation temperatures for both c-C3H2 and c-C3HD, we obtain a deuterium fraction of $14_{-3}^{+4}\%$.

Recent laboratory measurements have confirmed that chemical desorption (desorption of products due to exothermic surface reactions) can be an efficient process. The impact of including this process into gas-grain chemical models entirely depends on the formalism used and the associated parameters. Among these parameters, binding energies are probably the most uncertain ones for the moment. We propose a new model to compute binding energy of species to water ice surfaces. We have also compared the model results using either the new chemical desorption model proposed by Minissale et al. (2016) or the one of Garrod et al. (2007). The new binding energies have a strong impact on the formation of complex organic molecules. In addition, the new chemical desorption model from Minissale produces a much smaller desorption of these species and also of methanol. Combining the two effects, the abundances of CH3OH and COMs observed in cold cores cannot be reproduced by astrochemical models anymore.

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.

Recent detections of complex organic molecules in dark clouds have rekindled interest in the astrochemical modeling of these environments. Because of its relative closeness and rich molecular complexity, TMC-1 has been extensively observed to study the chemical processes taking place in dark clouds. We use local thermodynamical equilibrium radiative transfer modeling coupled with a Bayesian statistical method which takes into account outliers to analyze the data from the Nobeyama spectral survey of TMC-1 between 8 and 50 GHz. We compute the abundance relative to molecular hydrogen of 57 molecules, including 19 isotopologues in TMC-1 along with their associated uncertainty. The new results are in general agreement with previous abundance determination from Ohishi & Kaifu and the values reported in the review from Agundez & Wakelam. However, in some cases, large opacity and low signal to noise effects allow only upper or lower limits to be derived, respectively.

Protoplanetary disks are the target of many chemical studies (both observational and theoretical) as they contain the building material for planets. Their large vertical and radial gradients in density and temperature make them challenging objects for chemical models. In the outer part of these disks, the large densities and low temperatures provide a particular environment where the binding of species onto the dust grains can be very efficient and can affect the gas-phase chemical composition. We attempt to quantify to what extent the vertical abundance profiles and the integrated column densities of molecules predicted by a detailed gas-grain code are affected by the treatment of the molecular hydrogen physisorption at the surface of the grains. We performed three different models using the Nautilus gas-grain code. One model uses a H2 binding energy on the surface of water (440 K) and produces strong sticking of H2. Another model uses a small binding energy of 23 K (as if there were already a monolayer of H2), and the sticking of H$_2$ is almost negligible. Finally, the remaining model is an intermediate solution known as the encounter desorption mechanism. We show that the efficiency of molecular hydrogen binding (and thus its abundance at the surface of the grains) can have a quantitative effect on the predicted column densities in the gas phase of major species such as CO, CS, CN, and HCN.

The atomic hydrogen formation channels of the C + C2H2 reaction have been investigated using a continuous supersonic flow reactor over the 52 K 296 K temperature range. H-atoms were detected directly at 121.567 nm by vacuum ultraviolet laser induced fluorescence. Absolute H-atom yields were determined by comparison with the H-atom signal generated by the C + C2H4 reaction. The product yields agree with earlier crossed beam experiments employing universal detection methods. Incorporating these branching ratios in a gas-grain model of dense interstellar clouds increases the cC3H abundance. This reaction is a minor source of C3 containing molecules in the present simulations.

We conducted survey observations of a glycine precursor, methanimine or methylenimine (CH$_2$NH), with the NRO 45 m telescope and the SMT telescope towards 12 high-mass and two low-mass star-forming regions in order to increase number of CH$_2$NH sources and to better understand the characteristics of CH2NH sources. As a result of our survey, CH$_2$NH was detected in eight sources, including four new sources. The estimated fractional abundances were ~10$^8$ in Orion KL and G10.47+0.03, while they were ~10$^9$ towards the other sources. Our hydrogen recombination line and past studies suggest that CH$_2$NH-rich sources have less evolved HII regions. The less destruction rates by UV flux from the central star would be contributed to the high CH$_2$NH abundances towards CH$_2$NH-rich sources. Our gas-grain chemical simulations suggest that CH$_2$NH is mostly formed in the gas-phase by neutral-neutral reactions rather than the product of thermal evaporation from the dust surfaces.

We present an extended version of the 2-phase gas-grain code NAUTILUS to the 3-phase modelling of gas and grain chemistry of cold cores. In this model, both the mantle and the surface are considered as chemically active. We also take into account the competition among reaction, diffusion and evaporation. The model predictions are confronted to ice observations in the envelope of low-mass and massive young stellar objects as well as toward background stars. Modelled gas-phase abundances are compared to species observed toward TMC-1 (CP) and L134N dark clouds. We find that our model successfully reproduces the observed ice species. It is found that the reaction-diffusion competition strongly enhances reactions with barriers and more specifically reactions with H2, which is abundant on grains. This finding highlights the importance to have a good approach to determine the abundance of H2 on grains. Consequently, it is found that the major N-bearing species on grains go from NH3 to N2 and HCN when the reaction-diffusion competition is accounted. In the gas-phase and before few 10^5 yrs, we find that the 3-phase model does not have a strong impact on the observed species compared to the 2-phase model. After this time, the computed abundances dramatically decrease due to the strong accretion on dust, which is not counterbalanced by the desorption less efficient than in the 2-phase model. This strongly constrains the chemical-age of cold cores to be of the order of few 10^5 yrs.

(abridged) We used the IRAM 30-m to perform a sensitive wideband survey of 30 protoplanetary disks in the Taurus Auriga region. We simultaneously observed HCO$^+$(3-2), HCN(3-2), C$_2$H(3-2), CS(5-4), and two transitions of SO. We combine the results with a previous survey which observed $^{13}$CO (2-1), CN(2-1), two o-H$_2$CO lines and one of SO. We use available interferometric data to derive excitation temperatures of CN and C$_2$H in several sources. We determine characteristic sizes of the gas disks and column densities of all molecules using a parametric power-law disk model. Our study is mostly sensitive to molecules at 200-400 au from the stars. We compare the derived column densities to the predictions of an extensive gas-grain chemical disk model, under conditions representative of T Tauri disks. This survey provides 20 new detections of HCO$^+$ in disks, 18 in HCN, 11 in C$_2$H, 8 in CS and 4 in SO. HCO$^+$ is detected in almost all sources, and its J=3-2 line is essentially optically thick, providing good estimates of the disk radii. The other transitions are (at least partially) optically thin. Variations of the column density ratios do not correlate with any specific property of the star or disk. Disks around Herbig Ae stars appear less rich in molecules than those around T Tauri stars, although the sample remains small. SO is only found in the (presumably younger) embedded objects, perhaps reflecting an evolution of the S chemistry due to increasing depletion with time. Overall, the molecular column densities, and in particular the CN/HCN and CN/C$_2$H ratios, are well reproduced by gas-grain chemistry in cold disks. This study provides a census of simple molecules in disks of radii $> 200-300$ au. Extending that to smaller disks, or searching for less abundant or more complex molecules requires a much more sensitive facility, i.e. NOEMA and ALMA.

The product formation channels of ground state carbon atoms, C(3P), reacting with ammonia, NH3, have been investigated using two complementary experiments and electronic structure calculations. Reaction products are detected in a gas flow tube experiment (330 K, 4 Torr) using tunable VUV photoionization coupled with time of flight mass spectrometry. Temporal profiles of the species formed and photoionization spectra are used to identify primary products of the C + NH3 reaction. In addition, H-atom formation is monitored by VUV laser induced fluorescence from room temperature to 50 K in a supersonic gas flow generated by the Laval nozzle technique. Electronic structure calculations are performed to derive intermediates, transition states and complexes formed along the reaction coordinate. The combination of photoionization and laser induced fluorescence experiments supported by theoretical calculations indicate that in the temperature and pressure range investigated, the H + H2CN production channel represents 100% of the product yield for this reaction. Kinetics measurements of the title reaction down to 50 K and the effect of the new rate constants on interstellar nitrogen hydride abundances using a model of dense interstellar clouds are reported in paper II.

A continuous supersonic flow reactor has been used to measure rate constants for the C + NH3 reaction over the temperature range 50 to 296 K. C atoms were created by the pulsed laser photolysis of CBr4. The kinetics of the title reaction were followed directly by vacuum ultra-violet laser induced fluorescence (VUV LIF) of C loss and through H formation. The experiments show unambiguously that the reaction is rapid at 296 K, becoming faster at lower temperatures, reaching a value of 1.8 10-10 cm3 molecule-1 s-1 at 50 K. As this reaction is not currently included in astrochemical networks, its influence on interstellar nitrogen hydride abundances is tested through a dense cloud model including gas-grain interactions. In particular, the effect of the ortho-to-para ratio of H2 which plays a crucial role in interstellar NH3 synthesis is examined.

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

We present an extensive review of gas phase reactions producing methylacetylene and propene showing that these relatively abundant unsaturated hydrocarbons cannot be synthesized through gas-phase reactions. We explain the formation of propene and methylacetylene through surface hydrogenation of C3 depleted onto interstellar ices, C3 being a very abundant species in the gas phase.

The nature of the main sulphur reservoir in star forming regions is a long standing mystery. The observed abundance of sulphur-bearing species in dense clouds is only about 0.1 per cent of the same quantity in diffuse clouds. Therefore, the main sulphur species in star forming regions of the interstellar medium are still unknown. IRAS 16293-2422 is one of the regions where production of S-bearing species is favourable due to its conditions which allows the evaporation of ice mantles. We carried out observations in the 3 mm band towards the solar type protostar IRAS 16293-2422 with the IRAM 30m telescope. We observed a single frequency setup with the EMIR heterodyne 3 mm receiver with an Lower Inner (LI) tuning frequency of 89.98 GHz. Several lines of the complex sulphur species CH3SH were detected. Observed abundances are compared with simulations using the NAUTILUS gas-grain chemical model. Modelling results suggest that CH3SH has the constant abundance of 4e-9 (compared to H2) for radii lower than 200 AU and is mostly formed on the surfaces. Detection of CH3SH indicates that there may be several new families of S-bearing molecules (which could form starting from CH3SH) which have not been detected or looked for yet.

We present the detection of two H2C3O isomers, propynal and cyclopropenone, toward various starless cores and molecular clouds, together with upper limits for the third isomer propadienone. We review the processes controlling the abundances of H2C3O isomers in interstellar media showing that the reactions involved are gas-phase ones. We show that the abundances of these species are controlled by kinetic rather than thermodynamic effects.

The ketenyl radical (HCCO) has recently been discovered in two cold dense clouds with a non-negligible abundance of a few 1e-11 (compared to H2) (Agundez et al. 2015). Until now, no chemical network has been able to reproduce this observation. We propose here a chemical scheme that can reproduce HCCO abundances together with HCO, H2CCO and CH3CHO in the dark clouds Lupus-1A and L486. The main formation pathway for HCCO is the OH + CCH -> HCCO + H reaction as suggested by Agundez et al. (2015) but with a much larger rate coefficient than used in current models. Since this reaction has never been studied experimentally or theoretically, this larger value is based on a comparison with other similar systems.

Hereafter we describe the activities of the $Grand \, Sud-Ouest$ Data Centre operated for INSU/CNRS by the OMP-IRAP and the Universit\'e Paul Sabatier (Toulouse), in a collaboration with the OASU-LAB (Bordeaux) and OREME-LUPM (Montpellier).

The gas mass of protoplanetary disks, and the gas-to-dust ratio, are two key elements driving the evolution of these disks and the formation of planetary system. We explore here to what extent CO (or its isotopologues) can be used as a tracer of gas mass. We use a detailed gas-grain chemical model and study the evolution of the disk composition, starting from a dense pre-stellar core composition. We explore a range of disk temperature profiles, cosmic rays ionization rates, and disk ages for a disk model representative of T Tauri stars. At the high densities that prevail in disks, we find that, due to fast reactions on grain surfaces, CO can be converted to less volatile forms (principally s-CO$_2$, and to a lesser extent s-CH$_4$) instead of being evaporated over a wide range of temperature. The canonical gas-phase abundance of 10$^{-4}$ is only reached above about 30-35 K. The dominant Carbon bearing entity depends on the temperature structure and age of the disk. The chemical evolution of CO is also sensitive to the cosmic rays ionization rate. Larger gas phase CO abundances are found in younger disks. Initial conditions, such as parent cloud age and density, have a limited impact. This study reveals that CO gas-phase abundance is heavily dependent on grain surface processes, which remain very incompletely understood so far. The strong dependence on dust temperature profile makes CO a poor tracer of the gas-phase content of disks.

We attempt to determine the molecular composition of disks around young low-mass stars in the $\rho$ Oph region and to compare our results with a similar study performed in the Taurus-Auriga region. We used the IRAM 30 m telescope to perform a sensitive search for CN N=2-1 in 29 T Tauri stars located in the $\rho$ Oph and upper Scorpius regions. $^{13}$CO J=2-1 is observed simultaneously to provide an indication of the level of confusion with the surrounding molecular cloud. The bandpass also contains two transitions of ortho-H$_2$CO, one of SO, and the C$^{17}$O J=2-1 line, which provides complementary information on the nature of the emission. Contamination by molecular cloud in $^{13}$CO and even C$^{17}$O is ubiquitous. The CN detection rate appears to be lower than for the Taurus region, with only four sources being detected (three are attributable to disks). H$_2$CO emission is found more frequently, but appears in general to be due to the surrounding cloud. The weaker emission than in Taurus may suggest that the average disk size in the $\rho$ Oph region is smaller than in the Taurus cloud. Chemical modeling shows that the somewhat higher expected disk temperatures in $\rho$ Oph play a direct role in decreasing the CN abundance. Warmer dust temperatures contribute to convert CN into less volatile forms. In such a young region, CN is no longer a simple, sensitive tracer of disks, and observations with other tracers and at high enough resolution with ALMA are required to probe the gas disk population.

Chemical models used to study the chemical composition of the gas and the ices in the interstellar medium are based on a network of chemical reactions and associated rate coefficients. These reactions and rate coefficients are partially compiled from data in the literature, when available. We present in this paper kida.uva.2014, a new updated version of the kida.uva public gas-phase network first released in 2012. In addition to a description of the many specific updates, we illustrate changes in the predicted abundances of molecules for cold dense cloud conditions as compared with the results of the previous version of our network, kida.uva.2011.

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.

We present resolved Plateau de Bure Array observations of DM Tau in lines of HCO+ (3-2), (1-0) and DCO+ (3-2). A power-law fitting approach allowed a derivation of column densities of these two molecules. A chemical inner hole of ~50 AU was found in both HCO+ and DCO+ with DCO+ emission extending to only 450 AU. An isotopic ratio of R_D = N(DCO+) / N(HCO+) was found to range from 0.1 at 50 AU and 0.2 at 450 AU. Chemical modeling allowed an exploration of the sensitivity of these molecular abundances to physical parameters out with temperature, finding that X-rays were the domination ionization source in the HCO+ molecular region and that R_D also is sensitive to the CO depletion. The ionization fraction, assuming a steady state system, was found to be x(e-) ~ 10$^{-7}$. Modeling suggests that HCO+ is the dominant charged molecule in the disk but its contribution to ionization fraction is dwarfed by atmoic ions such as C+, S+ and H+.

Recent observations have revealed the existence of Complex Organic Molecules (COMs) in cold dense cores and prestellar cores. The presence of these molecules in such cold conditions is not well understood and remains a matter of debate since the previously proposed "warm- up" scenario cannot explain these observations. In this article, we study the effect of Eley- Rideal and complex induced reaction mechanisms of gas-phase carbon atoms with the main ice components of dust grains on the formation of COMs in cold and dense regions. Based on recent experiments we use a low value for the chemical desorption efficiency (which was previously invoked to explain the observed COM abundances). We show that our introduced mechanisms are efficient enough to produce a large amount of complex organic molecules in the gas-phase at temperatures as low as 10K.

The freezeout of gas-phase species onto cold dust grains can drastically alter the chemistry and the heating-cooling balance of protostellar material. In contrast to well-known species such as carbon monoxide (CO), the freezeout of various carriers of elements with abundances $<10^{-5}$ has not yet been well studied. Our aim here is to study the depletion of chlorine in the protostellar core, OMC-2 FIR 4. We observed transitions of HCl and H2Cl+ towards OMC-2 FIR 4 using the Herschel Space Observatory and Caltech Submillimeter Observatory facilities. Our analysis makes use of state of the art chlorine gas-grain chemical models and newly calculated HCl-H$_{2}$ hyperfine collisional excitation rate coefficients. A narrow emission component in the HCl lines traces the extended envelope, and a broad one traces a more compact central region. The gas-phase HCl abundance in FIR 4 is 9e-11, a factor of only 0.001 that of volatile elemental chlorine. The H2Cl+ lines are detected in absorption and trace a tenuous foreground cloud, where we find no depletion of volatile chlorine. Gas-phase HCl is the tip of the chlorine iceberg in protostellar cores. Using a gas-grain chemical model, we show that the hydrogenation of atomic chlorine on grain surfaces in the dark cloud stage sequesters at least 90% of the volatile chlorine into HCl ice, where it remains in the protostellar stage. About 10% of chlorine is in gaseous atomic form. Gas-phase HCl is a minor, but diagnostically key reservoir, with an abundance of <1e-10 in most of the protostellar core. We find the 35Cl/37Cl ratio in OMC-2 FIR 4 to be 3.2\pm0.1, consistent with the solar system value.

IRAS 16293-2422 is a well studied low-mass protostar characterized by a strong level of deuterium fractionation. In the line of sight of the protostellar envelope, an additional absorption layer, rich in singly and doubly deuterated water has been discovered by a detailed multiline analysis of HDO. To model the chemistry in this source, the gas-grain chemical code Nautilus has been used with an extended deuterium network. For the protostellar envelope, we solve the chemical reaction network in infalling fluid parcels in a protostellar core model. For the foreground cloud, we explored several physical conditions (density, cosmic ionization rate, C/O ratio). The main results of the paper are that gas-phase abundances of H2O, HDO and D2O observed in the inner regions of IRAS16293-2422 are lower than those predicted by a 1D dynamical/chemical (hot corino) model in which the ices are fully evaporated. The abundance in the outer part of the envelope present chaotic profiles due to adsorption/evaporation competition, very different from the constant abundance assumed for the analysis of the observations. We also found that the large abundances of gas-phase H2O, HDO and D2O observed in the absorption layer are more likely explained by exothermic surface reactions rather than photodesorption processes.

We review the reactions involving HCN and HNC in dark molecular clouds to elucidate new chemical sources and sinks of these isomers. We find that the most important reactions for the HCN-HNC system are Dissociative Recombination (DR) reactions of HCNH+ (HCNH+ + e-), the ionic CN + H3+, HCN + C+, HCN and HNC reactions with H+/He+/H3+/H3O+/HCO+, the N + CH2 reaction and two new reactions: H + CCN and C + HNC. We test the effect of the new rate constants and branching ratios on the predictions of gas-grain chemical models for dark cloud conditions. The rapid C + HNC reaction keeps the HCN/HNC ratio significantly above one as long as the carbon atom abundance remains high. However, the reaction of HCN with H3+ followed by DR of HCNH+ acts to isomerize HCN into HNC when carbon atoms and CO are depleted leading to a HCN/HNC ratio close to or slightly greater than 1. This agrees well with observations in TMC-1 and L134N taking into consideration the overestimation of HNC abundances through the use of the same rotational excitation rate constants for HNC as for HCN in many radiative transfer models.

Rate constants for the C(3P) + CH3OH reaction have been measured in a continuous supersonic flow reactor over the range 50 K to 296 K. C(3P) was created by the in-situ pulsed laser photolysis of CBr4, a multiphoton process which also produced some C(1D), allowing us to investigate simultaneously the low temperature kinetics of the C(1D) + CH3OH reaction. C(1D) atoms were followed by an indirect chemiluminescent tracer method in the presence of excess CH3OH. C(3P) atoms were detected by the same chemiluminescence technique and also by direct vacuum ultra-violet laser induced fluorescence (VUV LIF). Secondary measurements of product H(2S) atom formation have been undertaken allowing absolute H atom yields to be obtained by comparison with a suitable reference reaction. In parallel, statistical calculations have been performed based on ab-initio calculations of the complexes, adducts and transition states (TSs) relevant to the title reaction. By comparison with the experimental H atom yields, the preferred reaction pathways could be determined, placing important constraints on the statistical calculations. The experimental and theoretical work are in excellent agreement, predicting a negative temperature dependence of the rate constant increasing from 2.2 x 10-11 cm3 molecule-1 s-1 at 296 K to 20.0 x 10-11 cm3 molecule-1 s-1 at 50 K. CH3 and HCO are found to be the major products under our experimental conditions. As this reaction is not considered in current astrochemical networks, its influence on interstellar methanol abundances is tested using a model of dense interstellar clouds.

While recent studies of the solar-mass protostar IRAS16293-2422 have focused on its inner arcsecond, the wealth of Herschel/HIFI data has shown that the structure of the outer envelope and of the transition region to the more diffuse ISM is not clearly constrained. We use rotational ground-state transitions of CH (methylidyne), as a tracer of the lower-density envelope. Assuming LTE, we perform a $\chi^2$ minimization of the high spectral resolution HIFI observations of the CH transitions at ~532 and ~536 GHz in order to derive column densities in the envelope and in the foreground cloud. We obtain column densities of (7.7$\pm$0.2)$\times10^{13}$ cm$^{-2}$ and (1.5$\pm$0.3)$\times10^{13}$ cm$^{-2}$, respectively. The chemical modeling predicts column densities of (0.5-2)$\times10^{13}$ cm$^{-2}$ in the envelope (depending on the cosmic-ray ionization rate), and 5$\times10^{11}$ to 2.5$\times10^{14}$ cm$^{-2}$ in the foreground cloud (depending on time). Both observed abundances are reproduced by the model at a satisfactory level. The constraints set by these observations on the physical conditions in the foreground cloud are however weak. Furthermore, the CH abundance in the envelope is strongly affected by the rate coefficient of the reaction H+CH$\rightarrow$C+H$_2$ ; further investigation of its value at low temperature would be necessary to facilitate the comparison between the model and the observations.

Grain-surface reactions play an essential role in interstellar chemistry, since dust grain catalyses reactions at its surface allowing for the formation of molecules. We used a chemical model in which both gas-phase and grain-surface reactions occur and studied particularly the diffusion mechanisms on the surface of the grains. Surface reactions can occur via thermal hopping when species cross over a potential barrier or via quantum tunneling when species cross through this barrier. We show that the thermal diffusion (hopping) can be much more efficient after a cosmic ray particle collides with a dust grain, heating it to a peak temperature of 70K. We present here the results of numerical simulations after including the quantum tunneling mechanism for species H, H2 and O and considering the effect of cosmic ray particle collision on the surface reactions. As a consequence, the gas-phase and grain-surface abundances are affected and we show that more complex molecules can be formed in molecular clouds.

Protoplanetary disks composed of dust and gas are ubiquitous around young stars and are commonly recognized as nurseries of planetary systems. Their lifetime, appearance, and structure are determined by an interplay between stellar radiation, gravity, thermal pressure, magnetic field, gas viscosity, turbulence, and rotation. Molecules and dust serve as major heating and cooling agents in disks. Dust grains dominate the disk opacities, reprocess most of the stellar radiation, and shield molecules from ionizing UV/X-ray photons. Disks also dynamically evolve by building up planetary systems which drastically change their gas and dust density structures. Over the past decade significant progress has been achieved in our understanding of disk chemical composition thanks to the upgrade or advent of new millimeter/Infrared facilities (SMA, PdBI, CARMA, Herschel, e-VLA, ALMA). Some major breakthroughs in our comprehension of the disk physics and chemistry have been done since PPV. This review will present and discuss the impact of such improvements on our understanding of the disk physical structure and chemical composition.

We review the reactions between carbon chain molecules and radicals, namely Cn, CnH, CnH2, C2n+1O, CnN, HC2n+1N, with C, N and O atoms. Rate constants and branching ratios for these processes have been re-evaluated using experimental and theoretical literature data. In total 8 new species have been introduced, 41 new reactions have been proposed and 122 rate coefficients from kida.uva.2011 (Wakelam et al. 2012) have been modified. We test the effect of the new rate constants and branching ratios on the predictions of gas-grain chemical models for dark cloud conditions using two different C/O elemental ratios. We show that the new rate constants produce large differences in the predicted abundances of carbon chains since the formation of long chains is less effective. The general agreement between the model predictions and observed abundances in the dark cloud TMC-1 (CP) is improved by the new network and we find that C/O ratios of 0.7 and 0.95 both produce a similar agreement for different times. The general agreement for L134N (N) is not significantly changed. The current work specifically highlights the importance of O + CnH and N + CnH reactions. As there are very few experimental or theoretical data for the rate constants of these reactions we highlight the need for experimental studies of the O + CnH and N + CnH reactions, particularly at low temperature.

Predictions of astrochemical models depend strongly on the reaction rate coefficients used in the simulations. We reviewed a number of key reactions for the chemistry of nitrogen-bearing species in the dense interstellar medium and proposed new reaction rate coefficients for those reactions. The details of the reviews are given in the form of a datasheet associated with each reaction. The new recommended rate coefficients are given with an uncertainty and a temperature range of validity and will be included in KIDA (http://kida.obs.u-bordeaux1.fr).

Chemical models have been developed over the years by astrophysicists to study the pro- cesses at play in the various environments of the interstellar medium (ISM) that define the chemical composition of the gas and the dust. These qualitative aspects of the model predictions have been improved from a chemical point of view thanks to many recent developments of the experimental technics and theoretical methods that aim at studying the individual reactions in conditions as close to the ISM conditions as possible and characterize the rate constants of their efficiency. These models have also been more and more associated with dynamical evolution of the ISM physical conditions (for star forming regions for instance) since the chemical composition is far from steady-state in such regions. In this paper, we try to assess the state of the art concerning the chemical modeling of dark clouds, the initial step for the formation of stars and disks.

We investigate water and deuterated water chemistry in turbulent protoplanetary disks. Chemical rate equations are solved with the diffusion term, mimicking turbulent mixing in vertical direction. Water near the midplane is transported to the disk atmosphere by turbulence and destroyed by photoreactions to produce atomic oxygen, while the atomic oxygen is transported to the midplane and reforms water and/or other molecules. We find that this cycle significantly decreases column densities of water ice at r < 30 AU, where dust temperatures are too high to reform water ice effectively. The radial extent of such region depends on the desorption energy of atomic hydrogen. Our model indicates that water ice could be deficient even outside the sublimation radius. Outside this radius, the cycle decreases the D/H ratio of water ice from 2x10^-2, which is set by the collapsing core model, to 10^-4-10^-2 in 10^6 yr, without significantly decreasing the water ice column density. The resultant D/H ratios depend on the strength of mixing and the radial distance from the central star. Our finding suggests that the D/H ratio of cometary water (10^-3-10^-4) could be established (i.e. cometary water could be formed) in the solar nebula, even if the D/H ratio of water ice delivered to the disk was very high (10^-2).

We discuss models that astrochemists have developed to study the chemical composition of the interstellar medium. These models aim at computing the evolution of the chemical composition of a mixture of gas and dust under as- trophysical conditions. These conditions, as well as the geometry and the physical dynamics, have to be adapted to the objects being studied because different classes of objects have very different characteristics (temperatures, densities, UV radia- tion fields, geometry, history etc); e.g., proto-planetary disks do not have the same characteristics as protostellar envelopes. Chemical models are being improved continually thanks to comparisons with observations but also thanks to laboratory and theoretical work in which the individual processes are studied.

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

Rate constants for the potentially important interstellar N(4S) + CH(X2{\Pi}r) reaction have been measured in a continuous supersonic flow reactor over the range 56 K < T < 296 K using the relative rate technique employing both the N(4S) + OH(X2{\Pi}i) and N(4S) + CN(X2{\Sigma}+) reactions as references. Excess concentrations of atomic nitrogen were produced by the microwave discharge method upstream of the Laval nozzle and CH and OH radicals were created by the in-situ pulsed laser photolysis of suitable precursor molecules. In parallel, quantum dynamics calculations of the title reaction have been performed based on accurate global potential energy surfaces for the 13A' and 13A" states of HCN and HNC, brought about through a hierarchical construction scheme. Both adiabatic potential energy surfaces are barrierless, each one having two deep potential wells suggesting that this reaction is dominated by a complex-forming mechanism. The experimental and theoretical work are inexcellent agreement, predicting a positive temperature dependence of the rate constant, in contrast to earlier experimental work at low temperature. The effects of the new low temperature rate constants on interstellar N2 formation are tested using a dense cloud model, yielding N2 abundances 10-20 % lower than previously predicted.

The aim of this paper is to provide a new set of branching ratios for interstellar and planetary chemical networks based on a semi empirical model. We applied, instead of zero order theory (i.e. only the most exoergic decaying channel is considered), a statistical microcanonical model based on the construction of breakdown curves and using experimental high velocity collision branching ratios for their parametriza- tion. We applied the model to ion-molecule, neutral-neutral, and ion-pair reactions implemented in the few popular databases for astrochemistry such as KIDA, OSU and UMIST. We studied the reactions of carbon and hydrocarbon species with electrons, He+, H+, CH+, CH, C, and C+ leading to intermediate complexes of the type Cn=2,10, Cn=2,4 H, C3 H2, C+n=2,10, Cn=2,4 H+, or C3 H+2 . Comparison of predictions with measurements supports the validity of the model. Huge deviations with respect to database values are often obtained. Effects of the new branching ratios in time dependant chemistry for dark clouds and for photodissociation region chemistry with conditions similar to those found in the Horsehead Nebula are discussed.

Context: Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the environments of low-mass star-forming regions and, in particular, the Class 0 protostar IRAS 16293-2422. Aims: The key program Chemical HErschel Surveys of Star forming regions (CHESS) aims at studying the molecular complexity of the interstellar medium. The high sensitivity and spectral resolution of the Herschel/HIFI instrument provide a unique opportunity to observe the fundamental 1_{1,1}-0_{0,0} transition of ortho-D2O at 607 GHz and the higher energy 2_{1,2}-1_{0,1} transition of para-D2O at 898 GHz, both of which are inaccessible from the ground. Methods: The ortho-D2O transition at 607 GHz was previously detected. We present in this paper the first tentative detection for the para-D2O transition at 898 GHz. The spherical Monte Carlo radiative transfer code RATRAN was used to reproduce the observed line profiles of D2O with the same method that was used to reproduce the HDO and H2-18O line profiles in IRAS 16293-2422. Results: As for HDO, the absorption component seen on the D2O lines can only be reproduced by adding an external absorbing layer, possibly created by the photodesorption of the ices at the edges of the molecular cloud. The D2O column density is found to be about 2.5e12 cm^{-2} in this added layer, leading to a D2O/H2O ratio of about 0.5%. At a 3 sigma uncertainty, upper limits of 0.03% and 0.2% are obtained for this ratio in the hot corino and the colder envelope of IRAS 16293-2422, respectively. Conclusions: The deuterium fractionation derived in our study suggests that the ices present in IRAS 16293-2422 formed on warm dust grains (~15-20 K) in dense (~1e4-5e4 cm^{-3}) translucent clouds. These results allow us to address the earliest phases of star formation and the conditions in which ices form.

Turbulence is thought to be a key driver of the evolution of protoplanetary disks, regulating the mass accretion process, the transport of angular momentum, and the growth of dust particles. We intend to determine the magnitude of the turbulent motions in the outer parts (> 100 AU) of the disk surrounding DM Tau. Turbulent motions can be constrained by measuring the nonthermal broadening of line emission from heavy molecules. We used the IRAM Plateau de Bure interferometer to study emission from the CS molecule in the disk of DM Tau. High spatial (1.4 x 1 ") and spectral resolution (0.126 km/s) CS J=3-2 images provide constraints on the molecule distribution and velocity structure of the disk. A low sensitivity CS J=5-4 image was used in conjunction to evaluate the excitation conditions. We analyzed the data in terms of two parametric disk models, and compared the results with detailed time-dependent chemical simulations. The CS data confirm the relatively low temperature suggested by observations of other simple molecules. The intrinsic linewidth derived from the CS J=3-2 data is much larger than expected from pure thermal broadening. The magnitude of the derived nonthermal component depends only weakly on assumptions about the location of the CS molecules with respect to the disk plane. Our results indicate turbulence with a Mach number around 0.4 - 0.5 in the molecular layer. Geometrical constraints suggest that this layer is located near one scale height, in reasonable agreement with chemical model predictions.

We investigate the molecular evolution and D/H abundance ratios that develop as star formation proceeds from a dense-cloud core to a protostellar core, by solving a gas-grain reaction network applied to a 1-D radiative hydrodynamic model with infalling fluid parcels. Spatial distributions of gas and ice-mantle species are calculated at the first-core stage, and at times after the birth of a protostar. Gas-phase methanol and methane are more abundant than CO at radii $r\lesssim 100$ AU in the first-core stage, but gradually decrease with time, while abundances of larger organic species increase. The warm-up phase, when complex organic molecules are efficiently formed, is longer-lived for those fluid parcels in-falling at later stages. The formation of unsaturated carbon chains (warm carbon-chain chemistry) is also more effective in later stages; C$^+$, which reacts with CH$_4$ to form carbon chains, increases in abundance as the envelope density decreases. The large organic molecules and carbon chains are strongly deuterated, mainly due to high D/H ratios in the parent molecules, determined in the cold phase. We also extend our model to simulate simply the chemistry in circumstellar disks, by suspending the 1-D infall of a fluid parcel at constant disk radii. The species CH$_3$OCH$_3$ and HCOOCH$_3$ increase in abundance in $10^4-10^5$ yr at the fixed warm temperature; both also have high D/H ratios.

We investigate molecular evolution from a molecular cloud core to a first hydrostatic core in three spatial dimensions. We perform a radiation hydrodynamic simulation in order to trace fluid parcels, in which molecular evolution is investigated, using a gas-phase and grain-surface chemical reaction network. We derive spatial distributions of molecular abundances and column densities in the core harboring the first core. We find that the total of gas and ice abundances of many species in a cold era (10 K) remain unaltered until the temperature reaches ~500 K. The gas abundances in the warm envelope and the outer layer of the first core (T < 500 K) are mainly determined via the sublimation of ice-mantle species. Above 500 K, the abundant molecules, such as H2CO, start to be destroyed, and simple molecules, such as CO, H2O and N2 are reformed. On the other hand, some molecules are effectively formed at high temperature; carbon-chains, such as C2H2 and cyanopolyynes, are formed at the temperature of >700 K. We also find that large organic molecules, such as CH3OH and HCOOCH3, are associated with the first core (r < 10 AU). Although the abundances of these molecules in the first core stage are comparable or less than in the protostellar stage (hot corino), reflecting the lower luminosity of the central object, their column densities in our model are comparable to the observed values toward the prototypical hot corino, IRAS 16293-2422. We propose that these large organic molecules can be good tracers of the first cores.

The form of depleted sulphur in dense clouds is still unknown. Until now, only two molecules, OCS and SO2, have been detected in interstellar ices but cannot account for the elemental abundance of sulphur observed in diffuse medium. Chemical models suggest that solid H2S is the main form of sulphur in denser sources but observational constraints exist that infirm this hypothesis. We have used the Nautilus gas-grain code in which new chemical reactions have been added, based on recent experiments of H2S ice irradiation with UV photons and high energy protons. In particular, we included the new species Sn, H2Sn and C2S. We found that at the low temperature observed in dense clouds, i.e. 10 K, these new molecules are not efficiently produced and our modifications of the network do not change the previous pre- dictions. At slightly higher temperature, 20 K in less dense clouds or in the proximity of protostars, H2S abundance on the surfaces is strongly decreased in favor of the polysulfanes H2S3. Such a result can also be obtained if the diffusion barriers on the grains are less im- portant. In the context of the life cycle of interstellar clouds and the mixing between diffuse and denser parts of the clouds, the depletion of sulphur in the form of polysulfanes or other sulphur polymers, may have occurred in regions where the temperature is slightly higher than the cold inner parts of the clouds.

Molecular line emission from protoplanetary disks is a powerful tool to constrain their physical and chemical structure. Nevertheless, only a few molecules have been detected in disks so far. We take advantage of the enhanced capabilities of the IRAM 30m telescope by using the new broad band correlator (FTS) to search for so far undetected molecules in the protoplanetary disks surrounding the TTauri stars DM Tau, GO Tau, LkCa 15 and the Herbig Ae star MWC 480. We report the first detection of HC3N at 5 sigma in the GO Tau and MWC 480 disks with the IRAM 30-m, and in the LkCa 15 disk (5 sigma), using the IRAM array, with derived column densities of the order of 10^{12}cm^{-2}. We also obtain stringent upper limits on CCS (N < 1.5 x 10^{12} cm^{-3}). We discuss the observational results by comparing them to column densities derived from existing chemical disk models (computed using the chemical code Nautilus) and based on previous nitrogen and sulfur-bearing molecule observations. The observed column densities of HC3N are typically two orders of magnitude lower than the existing predictions and appear to be lower in the presence of strong UV flux, suggesting that the molecular chemistry is sensitive to the UV penetration through the disk. The CCS upper limits reinforce our model with low elemental abundance of sulfur derived from other sulfur-bearing molecules (CS, H2S and SO).

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

We present a novel chemical database for gas-phase astrochemistry. Named the KInetic Database for Astrochemistry (KIDA), this database consists of gas-phase reactions with rate coefficients and uncertainties that will be vetted to the greatest extent possible. Submissions of measured and calculated rate coefficients are welcome, and will be studied by experts before inclusion into the database. Besides providing kinetic information for the interstellar medium, KIDA is planned to contain such data for planetary atmospheres and for circumstellar envelopes. Each year, a subset of the reactions in the database (kida.uva) will be provided as a network for the simulation of the chemistry of dense interstellar clouds with temperatures between 10 K and 300 K. We also provide a code, named Nahoon, to study the time-dependent gas-phase chemistry of 0D and 1D interstellar sources.

The association reaction S + CO {\to} OCS + hnu has been identified as being particularly important for the prediction of gas-phase OCS abundances by chemical models of dark clouds. We performed detailed ab-initio calculations for this process in addition to undertaking an extensive review of the neutral-neutral reactions involving this species which might be important in such environments. The rate constant for this association reaction was estimated to be several orders of magnitude smaller than the one present in current astrochemical databases. The new rate for this reaction and the introduction of other processes, notably OH + CS {\to} OCS + H and C + OCS {\to} CO + CS, dramatically changes the OCS gas-phase abundance predicted by chemical models for dark clouds. The disagreement with observations in TMC-1 (CP) and L134N (N), suggests that OCS may be formed on grain surfaces as is the case for methanol. The observation of solid OCS on interstellar ices supports this hypothesis.

We study the content in S-bearing molecules of protoplanetary disks around low-mass stars. We used the new IRAM 30-m receiver EMIR to perform simultaneous observations of the $1_{10}-1_{01}$ line of H$_2$S at 168.8 GHz and $2_{23}-1_{12}$ line of SO at 99.3 GHz. We compared the observational results with predictions coming from the astrochemical code NAUTILUS, which has been adapted to protoplanetary disks. The data were analyzed together with existing CS J=3-2 observations. We fail to detect the SO and H$_2$S lines, although CS is detected in LkCa15, DM\,Tau, and GO\,Tau but not in MWC\,480. However, our new upper limits are significantly better than previous ones and allow us to put some interesting constraints on the sulfur chemistry. Our best modeling of disks is obtained for a C/O ratio of 1.2, starting from initial cloud conditions of H density of $2\times 10^5$ cm$^{-3}$ and age of $10^6$ yr. The results agree with the CS data and are compatible with the SO upper limits, but fail to reproduce the H$_2$S upper limits. The predicted H$_2$S column densities are too high by at least one order of magnitude. H$_2$S may remain locked onto grain surfaces and react with other species, thereby preventing the desorption of H$_2$S.

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

While unbiased surveys observable from ground-based telescopes have previously been obtained towards several high mass protostars, very little exists on low mass protostars. To fill up this gap, we carried out a complete spectral survey of the bands at 3, 2, 1 and 0.8 mm towards the solar type protostar IRAS16293-2422. The observations covered about 200\,GHz and were obtained with the IRAM-30m and JCMT-15m telescopes. Particular attention was devoted to the inter-calibration of the obtained spectra with previous observations. All the lines detected with more than 3 sigma and free from obvious blending effects were fitted with Gaussians to estimate their basic kinematic properties. More than 4000 lines were detected (with sigma \geq 3) and identified, yielding a line density of approximatively 20 lines per GHz, comparable to previous surveys in massive hot cores. The vast majority (~2/3) of the lines are weak and due to complex organic molecules. The analysis of the profiles of more than 1000 lines belonging 70 species firmly establishes the presence of two distinct velocity components, associated with the two objects, A and B, forming the IRAS16293-2422 binary system. In the source A, the line widths of several species increase with the upper level energy of the transition, a behavior compatible with gas infalling towards a ~1 Mo object. The source B, which does not show this effect, might have a much lower central mass of ~0.1 Mo. The difference in the rest velocities of both objects is consistent with the hypothesis that the source B rotates around the source A. This spectral survey, although obtained with single-dish telescope with a low spatial resolution, allows to separate the emission from 2 different components, thanks to the large number of lines detected. The data of the survey are public and can be retrieved on the web site http://www-laog.obs.ujf-grenoble.fr/heberges/timasss.

The main sulfur-bearing molecules OCS, H2S, SO, SO2, and CS have been observed in four high mass dense cores (W43-MM1, IRAS 18264, IRAS 05358, and IRAS 18162). Our goal is to put some constraints on the relative evolutionary stage of these sources by comparing these observations with time-dependent chemical modeling. We used the chemical model Nahoon, which computes the gas-phase chemistry and gas-grain interactions of depletion and evaporation. Mixing of the different chemical compositions shells in a 1D structure through protostellar envelope has been included since observed lines suggest nonthermal supersonic broadening. Observed radial profiles of the temperature and density are used to compute the chemistry as a function of time. With our model, we underproduce CS by several orders of magnitude compared to the other S-bearing molecules, which seems to contradict observations, although some uncertainties in the CS abundance observed at high temperature remain. The OCS/SO2, SO/SO2, and H2S/SO2 abundance ratios could in theory be used to trace the age of these massive protostars since they show a strong dependence with time, but the sources are too close in age compared to the accuracy of chemical models and observations. Our comparison between observations and modeling may, however, indicate that W43-MM1 could be chemically younger than the three other sources. Turbulent diffusivity through the protostellar envelopes has to be less efficient than 2e14 cm2s-1. Otherwise, it would have smoothed out the abundance profiles, and this would have been observed. The sulfur chemistry depends strongly on the 1D physical conditions. In our case, no conclusion can be given on the relative age of IRAS 18264, IRAS 18162 and IRAS 05358 except that they are very close. W43-MM1 seems younger than the other sources.

We survey the current situation regarding chemical modelling of the synthesis of molecules in the interstellar medium. The present state of knowledge concerning the rate coefficients and their uncertainties for the major gas-phase processes -- ion-neutral reactions, neutral-neutral reactions, radiative association, and dissociative recombination -- is reviewed. Emphasis is placed on those reactions that have been identified, by sensitivity analyses, as 'crucial' in determining the predicted abundances of the species observed in the interstellar medium. These sensitivity analyses have been carried out for gas-phase models of three representative, molecule-rich, astronomical sources: the cold dense molecular clouds TMC-1 and L134N, and the expanding circumstellar envelope IRC +10216. Our review has led to the proposal of new values and uncertainties for the rate coefficients of many of the key reactions. The impact of these new data on the predicted abundances in TMC-1 and L134N is reported. Interstellar dust particles also influence the observed abundances of molecules in the interstellar medium. Their role is included in gas-grain, as distinct from gas-phase only, models. We review the methods for incorporating both accretion onto, and reactions on, the surfaces of grains in such models, as well as describing some recent experimental efforts to simulate and examine relevant processes in the laboratory. These efforts include experiments on the surface-catalysed recombination of hydrogen atoms, on chemical processing on and in the ices that are known to exist on the surface of interstellar grains, and on desorption processes, which may enable species formed on grains to return to the gas-phase.

Fragmentation branching ratios of electronically excited molecular species are of first importance for the modeling of gas phase interstellar chemistry. Despite experimental and theoretical efforts that have been done during the last two decades there is still a strong lack of detailed information on those quantities for many molecules such as Cn, CnH or C3H2. Our aim is to provide astrochemical databases with more realistic branching ratios for Cn (n=2 to 10), CnH (n=2 to 4), and C3H2 molecules that are electronically excited either by dissociative recombination, photodissociation, or cosmic ray processes, when no detailed calculations or measurements exist in literature. High velocity collision in an inverse kinematics scheme was used to measure the complete fragmentation pattern of electronically excited Cn (n=2 to 10), CnH (n=2 to 4), and C3H2 molecules. Branching ratios of dissociation where deduced from those experiments. The full set of branching ratios was used as a new input in chemical models and branching ratio modification effects observed in astrochemical networks that describe the dense cold Taurus Molecular Cloud-1 and the photon dominated Horse Head region. The comparison between the branching ratios obtained in this work and other types of experiments showed a good agreement. It was interpreted as the signature of a statistical behavior of the fragmentation. The branching ratios we obtained lead to an increase of the C3 production together with a larger dispersion of the daughter fragments. The introduction of these new values in the photon dominated region model of the Horse Head nebula increases the abundance of C3 and C3H, but reduces the abundances of the larger Cn and hydrocarbons at a visual extinction Av smaller than 4.

In contrast with numerous studies on the physical and chemical structure of low- and high-mass protostars, much less is known about their intermediate-mass counterparts, a class of objects that could help to elucidate the mechanisms of star formation on both ends of the mass range. We present the first results from a rich HIFI spectral dataset on an intermediate-mass protostar, OMC2-FIR4, obtained in the CHESS (Chemical HErschel SurveyS of star forming regions) key programme. The more than 100 methanol lines detected between 554 and 961 GHz cover a range in upper level energy of 40 to 540 K. Our physical interpretation focusses on the hot core, but likely the cold envelope and shocked regions also play a role in reality, because an analysis of the line profiles suggests the presence of multiple emission components. An upper limit of 10^-6 is placed on the methanol abundance in the hot core, using a population diagram, large-scale source model and other considerations. This value is consistent with abundances previously seen in low-mass hot cores. Furthermore, the highest energy lines at the highest frequencies display asymmetric profiles, which may arise from infall around the hot core.

In contrast to extensively studied dense star-forming cores, little is known about diffuse gas surrounding star-forming regions. We study molecular gas in the high-mass star-forming region NGC6334I, which contains diffuse, quiescent components that are inconspicuous in widely used molecular tracers such as CO. We present Herschel/HIFI observations of CH toward NGC6334I observed as part of the CHESS key program. HIFI resolves the hyperfine components of its J=3/2-1/2 transition, observed in both emission and absorption. The CH emission appears close to the systemic velocity of NGC6334I, while its measured linewidth of 3 km/s is smaller than previously observed in dense gas tracers such as NH3 and SiO. The CH abundance in the hot core is 7 10^-11, two to three orders of magnitude lower than in diffuse clouds. While other studies find distinct outflows in, e.g., CO and H2O toward NGC6334I, we do not detect outflow signatures in CH. To explain the absorption signatures, at least two absorbing components are needed at -3.0 and +6.5 km/s with N(CH)=7 10^13 and 3 10^13 cm^-2. Two additional absorbing clouds are found at +8.0 and 0.0 km/s, both with N(CH)=2 10^13 cm^-2. Turbulent linewidths for the four absorption components vary between 1.5 and 5.0 km/s in FWHM. We constrain physical properties of our CH clouds by matching our CH absorbers with other absorption signatures. In the hot core, molecules such as H2O and CO trace gas that is heated and dynamically influenced by outflow activity, whereas CH traces more quiescent material. The four CH absorbers have column densities and turbulent properties consistent with diffuse clouds: two are located near NGC6334, and two are unrelated foreground clouds. Local density and dynamical effects influence the chemical composition of physical components of NGC6334, causing some components to be seen in CH but not in other tracers, and vice versa.

In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate species for the ammonia synthesis in the gas phase. Herschel/HIFI represents a unique opportunity to study the deuteration and formation mechanisms of such species, which are not observable from the ground. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration mechanism of this species. We observed the solar-mass Class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamental transition in absorption against the continuum background emitted from the nascent protostar. The 3 groups of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an efficient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N+ with HD, or deuteron/proton exchange with NH.

Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the environments of low-mass star forming regions, and in particular the Class 0 protostar IRAS 16293-2422. The CHESS (Chemical HErschel Surveys of Star forming regions) Key Program aims at studying the molecular complexity of the interstellar medium. The high sensitivity and spectral resolution of the HIFI instrument provide a unique opportunity to observe the fundamental 1,1,1 - 0,0,0 transition of the ortho-D2O molecule, inaccessible from the ground, and to determine the ortho-to-para D2O ratio. We have detected the fundamental transition of the ortho-D2O molecule at 607.35 GHz towards IRAS 16293-2422. The line is seen in absorption with a line opacity of 0.62 +/- 0.11 (1 sigma). From the previous ground-based observations of the fundamental 1,1,0 - 1,0,1 transition of para-D2O seen in absorption at 316.80 GHz we estimate a line opacity of 0.26 +/- 0.05 (1 sigma). We show that the observed absorption is caused by the cold gas in the envelope of the protostar. Using these new observations, we estimate for the first time the ortho to para D2O ratio to be lower than 2.6 at a 3 sigma level of uncertainty, to be compared with the thermal equilibrium value of 2:1.

We present observations of twelve rotational transitions of H2O-16, H2O-18, and H2O-17 toward the massive star-forming region NGC 6334 I, carried out with Herschel/HIFI as part of the guaranteed time key program Chemical HErschel Surveys of Star forming regions (CHESS). We analyze these observations to obtain insights into physical processes in this region. We identify three main gas components (hot core, cold foreground, and outflow) in NGC 6334 I and derive the physical conditions in these components. The hot core, identified by the emission in highly excited lines, shows a high excitation temperature of 200 K, whereas water in the foreground component is predominantly in the ortho- and para- ground states. The abundance of water varies between 4 10^-5 (outflow) and 10^-8 (cold foreground gas). This variation is most likely due to the freeze-out of water molecules onto dust grains. The H2O-18/H2O-17 abundance ratio is 3.2, which is consistent with the O-18/O-17 ratio determined from CO isotopologues. The ortho/para ratio in water appears to be relatively low 1.6(1) in the cold, quiescent gas, but close to the equilibrium value of three in the warmer outflow material (2.5(0.8)).

Abridged: We detail and benchmark two sophisticated chemical models developed by the Heidelberg and Bordeaux astrochemistry groups. The main goal of this study is to elaborate on a few well-described tests for state-of-the-art astrochemical codes covering a range of physical conditions and chemical processes, in particular those aimed at constraining current and future interferometric observations of protoplanetary disks. We consider three physical models: a cold molecular cloud core, a hot core, and an outer region of a T Tauri disk. Our chemical network (for both models) is based on the original gas-phase osu_03_2008 ratefile and includes gas-grain interactions and a set of surface reactions for the H-, O-, C-, S-, and N-bearing molecules. The benchmarking is performed with the increasing complexity of the considered processes: (1) the pure gas-phase chemistry, (2) the gas-phase chemistry with accretion and desorption, and (3) the full gas-grain model with surface reactions. Using atomic initial abundances with heavily depleted metals and hydrogen in its molecular form, the chemical evolution is modeled within 10^9 years. The time-dependent abundances calculated with the two chemical models are essentially the same for all considered physical cases and for all species, including the most complex polyatomic ions and organic molecules. This result however required a lot of efforts to make all necessary details consistent through the model runs, e.g. definition of the gas particle density, density of grain surface sites, the strength and shape of the UV radiation field, etc. The reference models and the benchmark setup, along with the two chemical codes and resulting time-dependent abundances are made publicly available in the Internet: http://www.mpia.de/homes/semenov/Chemistry_benchmark/home.html

We report the first detection of chloronium, H$_2$Cl$^+$, in the interstellar medium, using the HIFI instrument aboard the \emph{Herschel} Space Observatory. The $2_{12}-1_{01}$ lines of ortho-H$_2^{35}$Cl$^+$ and ortho-H$_2^{37}$Cl$^+$ are detected in absorption towards NGC~6334I, and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ is detected in absorption towards NGC~6334I and Sgr~B2(S). The H$_2$Cl$^+$ column densities are compared to those of the chemically-related species HCl. The derived HCl/H$_2$Cl$^+$ column density ratios, $\sim$1--10, are within the range predicted by models of diffuse and dense Photon Dominated Regions (PDRs). However, the observed H$_2$Cl$^+$ column densities, in excess of $10^{13}$~cm$^{-2}$, are significantly higher than the model predictions. Our observations demonstrate the outstanding spectroscopic capabilities of HIFI for detecting new interstellar molecules and providing key constraints for astrochemical models.

Gas and dust dissipation processes of proto-planetary disks are hardly known. Transition disks between Class II (proto-planetary disks) and Class III (debris disks) remain difficult to detect. We investigate the carbon chemistry of the peculiar CQ Tau gas disk. It is likely a transition disk because it exhibits weak CO emission with a relatively strong millimeter continuum, indicating that the disk might be currently dissipating its gas content. We used APEX to observe the two CI lines at 492GHz and 809 GHz in the disk orbiting CQ Tau. We compare the observations to several chemical model predictions. We focus our study on the influence of the stellar UV radiation shape and gas-to-dust ratio. We did not detect the CI lines. However, our upper limits are deep enough to exclude high-CI models. The only available models compatible with our limits imply very low gas-to-dust ratio, of the order of a few, only. These observations strengthen the hypothesis that CQ Tau is likely a transition disk and suggest that gas disappears before dust.

The outflow driven by the low-mass class 0 protostar L1157 is the prototype of the so-called chemically active outflows. The bright bowshock B1 in the southern outflow lobe is a privileged testbed of magneto-hydrodynamical (MHD) shock models, for which dynamical and chemical processes are strongly interdependent. We present the first results of the unbiased spectral survey of the L1157-B1 bowshock, obtained in the framework of the key program "Chemical Herschel Surveys of Star Forming Regions" (CHESS). The main aim is to trace the warm and chemically enriched gas and to infer the excitation conditions in the shock region. The CO 5-4 and H2O lines have been detected at high-spectral resolution in the unbiased spectral survey of the HIFI-Band 1b spectral window (555-636 GHz), presented by Codella et al. in this volume. Complementary ground-based observations in the submm window help establish the origin of the emission detected in the main-beam of HIFI, and the physical conditions in the shock.}{Both lines exhibit broad wings, which extend to velocities much higher than reported up to now. We find that the molecular emission arises from two regions with distinct physical conditions: an extended, warm (100K), dense (3e5 cm-3) component at low-velocity, which dominates the water line flux in Band~1; a secondary component in a small region of B1 (a few arcsec) associated with high-velocity, hot (> 400 K) gas of moderate density ((1.0-3.0)e4 cm-3), which appears to dominate the flux of the water line at 179mu observed with PACS. The water abundance is enhanced by two orders of magnitude between the low- and the high-velocity component, from 8e-7 up to 8e-5. The properties of the high-velocity component agree well with the predictions of steady-state C-shock models.

We identify a prominent absorption feature at 1115 GHz, detected in first HIFI spectra towards high-mass star-forming regions, and interpret its astrophysical origin. The characteristic hyperfine pattern of the H2O+ ground-state rotational transition, and the lack of other known low-energy transitions in this frequency range, identifies the feature as H2O+ absorption against the dust continuum background and allows us to derive the velocity profile of the absorbing gas. By comparing this velocity profile with velocity profiles of other tracers in the DR21 star-forming region, we constrain the frequency of the transition and the conditions for its formation. In DR21, the velocity distribution of H2O+ matches that of the [CII] line at 158\mu\m and of OH cm-wave absorption, both stemming from the hot and dense clump surfaces facing the HII-region and dynamically affected by the blister outflow. Diffuse foreground gas dominates the absorption towards Sgr B2. The integrated intensity of the absorption line allows us to derive lower limits to the H2O+ column density of 7.2e12 cm^-2 in NGC 6334, 2.3e13 cm^-2 in DR21, and 1.1e15 cm^-2 in Sgr B2.

Because of new telescopes that will dramatically improve our knowledge of the interstellar medium, chemical models will have to be used to simulate the chemistry of many regions with diverse properties. To make these models more robust, it is important to understand their sensitivity to a variety of parameters. In this article, we report a study of the sensitivity of a chemical model of a cold dense core, with homogeneous and time-independent physical conditions, to variations in the following parameters: initial chemical inventory, gas temperature and density, cosmic-ray ionization rate, chemical reaction rate coefficients, and elemental abundances. From the results of the parameter variations, we can quantify the sensitivity of the model to each parameter as a function of time. Our results can be used in principle with observations to constrain some parameters for different cold clouds. We also attempted to use the Monte Carlo approach with all parameters varied collectively. Within the parameter ranges studied, the most critical parameters turn out to be the reaction rate coefficients at times up to 4e5 yr and elemental abundances at later times. At typical times of best agreement with observation, models are sensitive to both of these parameters. The models are less sensitive to other parameters such as the gas density and temperature. The improvement of models will require that the uncertainties in rate coefficients of important reactions be reduced. As the chemistry becomes better understood and more robust, it should be possible to use model sensitivities concerning other parameters, such as the elemental abundances and the cosmic ray ionization rate, to yield detailed information on cloud properties and history. Nevertheless, at the current stage, we cannot determine the best values of all the parameters simultaneously based on purely observational constraints.

We studied several representative circumstellar disks surrounding the Herbig Ae star MWC 480 and the T Tauri stars LkCa 15 and DM Tau at (sub-)millimeter wavelengths in lines of CCH. Our aim is to characterize photochemistry in the heavily UV-irradiated MWC 480 disk and compare the results to the disks around cooler T Tauri stars. We detected and mapped CCH in these disks with the IRAM Plateau de Bure Interferome- ter in the C- and D-configurations in the (1-0) and (2-1) transitions. Using an iterative minimization technique, the CCH column densities and excitation conditions are con- strained. Very low excitation temperatures are derived for the T Tauri stars. These values are compared with the results of advanced chemical modeling, which is based on a steady-state flared disk structure with a vertical temperature gradient, and a gas- grain chemical network with surface reactions. Both model and observations suggest that CCH is a sensitive tracer of the X-ray and UV irradiation. The predicted radial dependency and source to source variations of CCH column densities qualitatively agree with the observed values, but the predicted column densities are too low by a factor of several. The chemical model fails to reproduce high concentrations of CCH in very cold disk midplane as derived from the observed low excitation condition for both the (1-0) and (2-1) transitions.

Chemical composition of the massive cores forming high-mass stars can put some constrains on the time scale of the massive star formation: sulphur chemistry is of specific interest due to its rapid evolution in warm gas and because the abundance of sulphur bearing species increases significantly with the temperature. Two mid-infrared quiet and two brighter massive cores are observed in various transitions (E_up up to 289K) of CS, OCS, H2S, SO, SO2 and of their isotopologues at mm wavelengths with the IRAM 30m and CSO telescopes. 1D modeling of the dust continuum is used to derive the density and temperature laws, which are then applied in the RATRAN code to model the observed line emission, and to derive the relative abundances of the molecules. All lines, except the highest energy SO2 transition, are detected. Infall (up to 2.9km/s) may be detected towards the core W43MM1. The inferred mass rate is 5.8-9.4 10^{-2} M_{\odot}/yr. We propose an evolutionary sequence of our sources (W43MM1-IRAS18264-1152-IRAS05358+3543-IRAS18162-2048), based on the SED analysis. The analysis of the variations in abundance ratios from source to source reveals that the SO and SO2 relative abundances increase with time, while CS and OCS decrease. Molecular ratios, such as [OCS/H2S], [CS/H2S], [SO/OCS], [SO2/OCS], [CS/SO] and [SO2/SO] may be good indicators of evolution depending on layers probed by the observed molecular transitions. Observations of molecular emission from warmer layers, hence involving higher upper energy levels are mandatory to include.

Chemical networks used for models of interstellar clouds contain many reactions, some of them with poorly determined rate coefficients and/or products. In this work, we report a method for improving the predictions of molecular abundances using sensitivity methods and ab initio calculations. Based on the chemical network osu.2003, we used two different sensitivity methods to determine the most important reactions as a function of time for models of dense cold clouds. Of these reactions, we concentrated on those between C and C3 and between C and C5, both for their effect on specific important species such as CO and for their general effect on large numbers of species. We then used ab initio and kinetic methods to determine an improved rate coefficient for the former reaction and a new set of products, plus a slightly changed rate coefficient for the latter. Putting our new results in a pseudo-time-dependent model of cold dense clouds, we found that the abundances of many species are altered at early times, based on large changes in the abundances of CO and atomic C. We compared the effect of these new rate coefficients/products on the comparison with observed abundances and found that they shift the best agreement from 3e4 yr to (1-3)e5 yr.

Aims. We attempt to understand the presence of gas phase CO below its freezing temperature in circumstellar disks. We study two promising mechanisms to explain this phenomenon: turbulent mixing and photodesorption. Methods. We compute the chemical evolution of circumstellar disks including grain surface reactions with and without turbulent mixing and CO photodesorption. Results. We show that photodesorption significantly enhances the gas phase CO abundance, by extracting CO from the grains when the visual extinction remains below about 5 magnitudes. However the resulting dependence of column density on radial distance is not consistent with observations so far. We propose that this inconsistency could be the result of grain growth. On the other hand, the influence of turbulent mixing is not found to be straightforward. The efficiency of turbulent mixing depends upon a variety of parameters, including the disk structure. For the set of parameters we chose, turbulent mixing is not found to have any significant influence on the CO column density.

Dielectronic recombination (DR) of singly charged ions is a reaction pathway that is commonly neglected in chemical models of molecular clouds. In this study we include state-of-the-art DR data for He$^+$, C$^+$, N$^+$, O$^+$, Na$^+$, and Mg$^+$ in chemical models used to simulate dense molecular clouds, protostars, and diffuse molecular clouds. We also update the radiative recombination (RR) rate coefficients for H$^+$, He$^+$, C$^+$, N$^+$, O$^+$, Na$^+$, and Mg$^+$ to the current state-of-the-art values. The new RR data has little effect on the models. However, the inclusion of DR results in significant differences in gas-grain models of dense, cold molecular clouds for the evolution of a number of surface and gas-phase species. We find differences of a factor of 2 in the abundance for 74 of the 655 species at times of $10^4$--$10^6$ years in this model when we include DR. Of these 74 species, 16 have at least a factor of 10 difference in abundance. We find the largest differences for species formed on the surface of dust grains. These differences are due primarily to the addition of C$^+$ DR, which increases the neutral C abundance, thereby enhancing the accretion of C onto dust. These results may be important for the warm-up phase of molecular clouds when surface species are desorbed into the gas phase. We also note that no reliable state-of-the-art RR or DR data exist for Si$^+$, P$^+$, S$^+$, Cl$^+$, and Fe$^+$. Modern calculations for these ions are needed to better constrain molecular cloud models.

We study the formation conditions of icy planetesimals in protoplanetary disks in order to determine the composition of ices in small and cold extrasolar planets. Assuming that ices are formed from hydrates, clathrates, and pure condensates, we calculate their mass fractions with respect to the total quantity of ices included in planetesimals, for a grid of disk models. We find that the composition of ices weakly depends on the adopted disk thermodynamic conditions, and is rather influenced by the initial composition of the gas phase. The use of a plausible range of molecular abundance ratios and the variation of the relative elemental carbon over oxygen ratio in the gas phase of protoplanetary disks, allow us to apply our model to a wide range of planetary systems. Our results can thus be used to constrain the icy/volatile phase composition of cold planets evidenced by microlensing surveys, hypothetical ocean-planets and carbon planets, which could be detected by Corot or Kepler.

Virtually all detailed gas-phase models of the chemistry of dense interstellar clouds exclude polycyclic aromatic hydrocarbons (PAH's). This omission is unfortunate because from the few studies that have been done on the subject, it is known that the inclusion of PAH's can affect the gas-phase chemistry strongly. We have added PAH's to our network to determine the role they play in the chemistry of cold dense cores. In the models presented here, we include radiative attachment to form PAH-, mutual neutralization between PAH anions and small positively-charged ions, and photodetachment. We also test the sensitivity of our results to changes in the size and abundance of the PAH's. Our results confirm that the inclusion of PAH's changes many of the calculated abundances of smaller species considerably. In TMC-1, the general agreement with observations is significantly improved contrary to L134N. This may indicate a difference in PAH properties between the two regions. With the inclusion of PAH's in dense cloud chemistry, high-metal elemental abundances give a satisfactory agreement with observations. As a result, we do not need to decrease the observed elemental abundances of all metals and we do not need to vary the elemental C/O ratio in order to produce large abundances of carbon species in TMC-1 (CP).

We investigate molecular evolution in a star-forming core that is initially a hydrostatic starless core and collapses to form a low-mass protostar. The results of a one-dimensional radiation-hydrodynamics calculation are adopted as a physical model of the core. We first derive radii at which CO and large organic species sublimate. CO sublimation in the central region starts shortly before the formation of the first hydrostatic core. When the protostar is born, the CO sublimation radius extends to 100 AU, and the region inside $\lesssim 10$ AU is hotter than 100 K, at which some large organic species evaporate. We calculate the temporal variation of physical parameters in infalling shells, in which the molecular evolution is solved using an updated gas-grain chemical model to derive the spatial distribution of molecules in a protostellar core. The shells pass through the warm region of $10 -100$ K in several $\times$ $10^4$ yr, and fall into the central star $\sim 100$ yr after they enter the region where $T \gtrsim 100$ K. We find that large organic species are formed mainly via grain-surface reactions at temperatures of $20 -40$ K and then desorbed into the gas-phase at their sublimation temperatures. Carbon-chain species can be formed by a combination of gas-phase reactions and grain-surface reactions following the sublimation of CH$_4$. Our model also predicts that CO$_2$ is more abundant in isolated cores, while gas-phase large organic species are more abundant in cores embedded in ambient clouds.

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.

Aims: The gas-phase abundance of methanol in dark quiescent cores in the interstellar medium cannot be explained by gas-phase chemistry. In fact, the only possible synthesis of this species appears to be production on the surfaces of dust grains followed by desorption into the gas. Yet, evaporation is inefficient for heavy molecules such as methanol at the typical temperature of 10 K. It is necessary then to consider non-thermal mechanisms for desorption. But, if such mechanisms are considered for the production of methanol, they must be considered for all surface species. Methods: Our gas-grain network of reactions has been altered by the inclusion of a non-thermal desorption mechanism in which the exothermicity of surface addition reactions is utilized to break the bond between the product species and the surface. Our estimated rate for this process derives from a simple version of classical unimolecular rate theory with a variable parameter only loosely contrained by theoretical work. Results: Our results show that the chemistry of dark clouds is altered slightly at times up to 10^6 yr, mainly by the enhancement in the gas-phase abundances of hydrogen-rich species such as methanol that are formed on grain surfaces. At later times, however, there is a rather strong change. Instead of the continuing accretion of most gas-phase species onto dust particles, a steady-state is reached for both gas-phase and grain-surface species, with significant abundances for the former. Nevertheless, most of the carbon is contained in an undetermined assortment of heavy surface hydrocarbons. Conclusions: The desorption mechanism discussed here will be better constrained by observational data on pre-stellar cores, where a significant accretion of species such as CO has already occurred.

Aim: To determine whether or not gas-phase chemical models with homogeneous and time-independent physical conditions explain the many observed molecular abundances in astrophysical sources, it is crucial to estimate the uncertainties in the calculated abundances and compare them with the observed abundances and their uncertainties. Non linear amplification of the error and bifurcation may limit the applicability of chemical models. Here we study such effects on dense cloud chemistry. Method: Using a previously studied approach to uncertainties based on the representation of rate coefficient errors as log normal distributions, we attempted to apply our approach using as input a variety of different elemental abundances from those studied previously. In this approach, all rate coefficients are varied randomly within their log normal (Gaussian) distribution, and the time-dependent chemistry calculated anew many times so as to obtain good statistics for the uncertainties in the calculated abundances. Results: Starting with so-called ``high-metal'' elemental abundances, we found bimodal rather than Gaussian like distributions for the abundances of many species and traced these strange distributions to an extreme sensitivity of the system to changes in the ratio of the cosmic ray ionization rate zeta\_He for He and that for molecular hydrogen zeta\_H2. The sensitivity can be so extreme as to cause a region of bistability, which was subsequently found to be more extensive for another choice of elemental abundances. To the best of our knowledge, the bistable solutions found in this way are the same as found previously by other authors, but it is best to think of the ratio zeta\_He/zeta\_H2 as a control parameter perpendicular to the ''standard'' control parameter zeta/n\_H.

The gas-phase chemistry of dark clouds has been studied with a treatment of uncertainties caused both by errors in individual rate coefficients and uncertainties in physical conditions. Moreover, a sensitivity analysis has been employed to attempt to determine which reactions are most important in the chemistry of individual species. The degree of overlap between calculated errors in abundances and estimated observational errors has been used as an initial criterion for the goodness of the model and the determination of a best 'chemical' age of the source. For the well-studied sources L134N and TMC-1CP, best agreement is achieved at so-called "early times" ~10$^{5}$ yr, in agreement with previous calculations but here put on a firmer statistical foundation. A more detailed criterion for agreement, which takes into account the degree of disagreement, is also proposed. Poorly understood but critical classes of reactions are delineated, especially reactions between ions and polar neutrals. Such reactions will have to be understood better before the chemistry can be made more secure. Nevertheless, the level of agreement is low enough to indicate that a static picture of physical conditions without consideration of interactions with grain surfaces is inappropriate for a complete understanding of the chemistry.

We report a detailed attempt to model the ortho-to-para abundance ratio of c-C3H2 so as to reproduce observed values in the cores of the well-known source TMC-1. According to observations, the ortho-to-para ratios vary, within large uncertainties, from a low of near unity to a high of approximately three depending on the core. We used the osu.2003 network of gas-phase chemical reactions augmented by reactions that specifically consider the formation, depletion, and interconversion of the ortho and para forms of the c-C3H2 and its precursor ion c-C3H3+. We investigated the sensitivity of the calculated ortho-to-para ratio for c-C3H2 to a large number of factors. For the less evolved cores C, CP, and D, we had no difficulty reproducing the observed ortho-to-para ratios of 1-2.

It is not common to consider the role of uncertainties in the rate coefficients used in interstellar gas-phase chemical models. In this paper, we report a new method to determine both the uncertainties in calculated molecular abundances and their sensitivities to underlying uncertainties in the kinetic data utilized. The method is used in hot core models to determine if previous analyses of the age and the applicable cosmic-ray ionization rate are valid. We conclude that for young hot cores ($\le 10^4$ yr), the modeling uncertainties related to rate coefficients are reasonable so that comparisons with observations make sense. On the contrary, the modeling of older hot cores is characterized by strong uncertainties for some of the important species. In both cases, it is crucial to take into account these uncertainties to draw conclusions from the comparison of observations with chemical models.

We present the results of a multiline survey at mm-wavelengths of the Cepheus A star forming region. Four main flows have been identified: three pointing in the SW, NE, and SE directions and accelerating high density CS clumps. The fourth outflow, revealed by high-sensitivity HDO observations, is pointing towards South and is associated with conditions particularly favourable to a chemical enrichment. The analysis of the line profiles shows that the SiO molecule dominates at the highest velocities and at the highest excitation conditions, confirming its close association with shocks. H2S, SO2, and SO preferentially trace more quiescent regions than SiO, and in particular a lack of bright H2S emission at the highest velocities is found. OCS and H2CS emit at quite high velocities, where the abundances of three shock tracers like SiO, CH3OH, and HDO are higher. These results may indicate that H2S is not the only major sulphur carrier in the grain mantles, and that OCS and H2CS may probably play an important role on the grains; or that alternatively they rapidly form once the mantle is evaporated after the passage of a shock. Finally, the outflow peak emission has been compared with recent time-dependent sulphur chemistry models: the results indicate that, if associated with accurate measurements of the physical conditions, the CH3OH/H2CS column density ratio can be used as an effective chemical clock to date the age of shocked gas.

We present SO and SO2 observations in the region of the low mass protostar IRAS2/NGC1333. The East-West outflow originating from this source has been mapped in four transitions of both SO and SO2. In addition, CS observations published in the literature have been used. We compute the SO, SO2 and CS column densities and the physical conditions at several positions of the outflow using LTE and a non-LTE LVG approximations. The SO2/SO and CS/SO abundance ratios are compared with the theoretical predictions of a chemical model adapted to the physical conditions in the IRAS2 outflow. The SO2/SO abundance ratios are constant in the two lobes of the outflow whereas CS/SO is up to 6 times lower in the shocked gas of the East lobe than in the West one. The comparison with the chemical model allows us to constrain the age of the outflow produced by IRAS2 to >5e3yr. We find low densities and temperatures for the outflow of IRAS2 (< 1e6 cm-3 and >70 K) from SO and SO2 emission probably because the two molecules trace the cooled entrained material. The East lobe of the outflow shows denser gas compared to the West lobe. We also discuss some constraints on the depleted form of sulphur.

We present IRAM 30m and JCMT observations of HDO lines towards the solar-type protostar IRAS 16293-2422. Five HDO transitions have been detected on-source, and two were unfruitfully searched for towards a bright spot of the outflow of IRAS 16293-2422. We interpret the data by means of the Ceccarelli, Hollenbach and Tielens (1996) model, and derive the HDO abundance in the warm inner and cold outer parts of the envelope. The emission is well explained by a jump model, with an inner abundance of 1e-7 and an outer abundance lower than 1e-9 (3 sigma). This result is in favor of HDO enhancement due to ice evaporation from the grains in theinner envelope. The deuteration ratio HDO/H2O is found to be f_in=3% and f_out < 0.2% (3 sigma) in the inner and outer envelope respectively and therefore, the fractionation also undergoes a jump in the inner part of the envelope. These results are consistent with the formation of water in the gas phase during the cold prestellar core phase and storage of the molecules on the grains, but do not explain why observations of H2O ices consistently derive a H2O ice abundance of several 1e-5 to 1e-4, some two orders of magnitude larger than the gas phase abundance of water in the hot core around IRAS 16293-2422.

We report a theoretical study of sulphur chemistry, as applied to hot cores, where S-bearing molecular ratios have been previously proposed and used as chemical clocks. As in previous models, we follow the S-bearing molecular composition after the injection of grain mantle components into the gas phase. For this study, we developed a time-dependent chemical model with up-to-date reaction rate coefficients. We ran several cases, using different realistic chemical compositions for the grain mantles and for the gas prior to mantle evaporation. The modeling shows that S-bearing molecular ratios depend very critically on the gas temperature and density, the abundance of atomic oxygen, and, most importantly, on the form of sulphur injected in the gas phase, which is very poorly known. Consequently, ratios of S-bearing molecules cannot be easily used as chemical clocks. However, detailed observations and careful modeling of both physical and chemical structure can give hints on the source age and constrain the mantle composition (i.e. the form of sulphur in cold molecular clouds) and, thus, help to solve the mystery of the sulphur depletion. We analyse in detail the cases of Orion and IRAS16293-2422. The comparison of the available observations with our model suggests that the majority of sulphur released from the mantles is mainly in, or soon converted into, atomic form.

We report observations of the expected main S-bearing species (SO, SO2 and H2S) in the low-mass star forming region L1689N. We obtained large scale (~300''x200'') maps of several transitions from these molecules with the goal to study the sulphur chemistry, i.e. how the relative abundances change in the different physical conditions found in L1689N. We identified eight interesting regions, where we carried out a quantitative comparative study: the molecular cloud (as reference position), five shocked regions caused by the interaction of the molecular outflows with the cloud, and the two protostars IRAS16293-2422 and 16293E. In the cloud we carefully computed the gas temperature and density by means of a non-LTE LVG code, while in other regions we used previous results. We hence derived the column density of SO, SO2 and H2S, together with SiO and H2CO - which were observed previously - and their relevant abundance ratios. We find that SiO is the molecule that shows the largest abundance variations in the shocked regions, whereas S-bearing molecules show more moderate variations. Remarkably, the region of the brightest SiO emission in L1689N is undetected in SO2, H2S and H2CO and only marginally detected in SO. In the other weaker SiO shocks, SO2 is enhanced with respect to SO. We propose a schema in which the different molecular ratios correspond to different ages of the shocks. Finally, we find that SO, SO2 and H2S have significant abundance jumps in the inner hot core of IRAS16293-2422 and discuss the implications of the measured abundances.