Showing 1–50 of 59 results for author: Sipilä, O

A low cosmic-ray ionisation rate in the prestellar core Ophiuchus/H-MM1. Mapping of the molecular ions ortho-H2D+, N2H+, and DCO+

Authors: Jorma Harju, Charlotte Vastel, Olli Sipilae, Elena Redaelli, Paola Caselli, Jaime E. Pineda, Arnaud Belloche, Friedrich Wyrowski

Abstract: (abridged) We have mapped the prestellar core H-MM1 in Ophiuchus in rotational lines of ortho-H2D+ (oH2D+), N2H+, and DCO+ at the wavelength 0.8 mm with the Large APEX sub-Millimeter Array (LAsMA) multibeam receiver of the Atacama Pathfinder EXperiment (APEX) telescope. We also ran a series of chemistry models to predict the abundance distributions of the observed molecules, and to estimate the ef… ▽ More (abridged) We have mapped the prestellar core H-MM1 in Ophiuchus in rotational lines of ortho-H2D+ (oH2D+), N2H+, and DCO+ at the wavelength 0.8 mm with the Large APEX sub-Millimeter Array (LAsMA) multibeam receiver of the Atacama Pathfinder EXperiment (APEX) telescope. We also ran a series of chemistry models to predict the abundance distributions of the observed molecules, and to estimate the effect of the cosmic-ray ionisation rate on their abundances. The three line maps show different distributions. The oH2D+ map is extended and outlines the general structure of the core, while N2H+ mainly shows the density maxima, and the DCO+ emission peaks are shifted towards one edge of the core where a region of enhanced desorption has been found previously. According to the chemical simulation, the fractional oH2D+ abundance remains relatively high in the centre of the core, and its column density correlates strongly with the cosmic-ray ionisation rate. Simulated line maps constrain the cosmic-ray ionisation rate per hydrogen molecule to be low, between 5e-18/s and 1e-17/s in the H-MM1 core. This estimate agrees with the gas temperature measured in the core. Modelling line emission of oH2D+ provides a straightforward method of determining the cosmic-ray ionisation rate in dense clouds, where the primary ion, H3+, is not observable. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: accepted for publication in Astronomy & Astrophysics

Fractionation in young cores: Direct determinations of nitrogen and carbon fractionation in HCN

Authors: S. S. Jensen, S. Spezzano, P. Caselli, O. Sipilä, E. Redaelli, K. Giers, J. Ferrer Asensio

Abstract: We aim to determine the $^{14}$N/$^{15}$N and $^{12}$C/$^{13}$C ratios for HCN in six starless and prestellar cores and compare the results between the direct method using radiative transfer modeling and the indirect double isotope method assuming a fixed $^{12}$C/$^{13}$C ratio. We present IRAM 30m observations of the HCN 1-0, HCN 3-2, HC15N 1-0 and H13CN 1-0 transitions toward six embedded cores… ▽ More We aim to determine the $^{14}$N/$^{15}$N and $^{12}$C/$^{13}$C ratios for HCN in six starless and prestellar cores and compare the results between the direct method using radiative transfer modeling and the indirect double isotope method assuming a fixed $^{12}$C/$^{13}$C ratio. We present IRAM 30m observations of the HCN 1-0, HCN 3-2, HC15N 1-0 and H13CN 1-0 transitions toward six embedded cores. The ${}^{14}$N/${}^{15}$N ratio was derived using both the indirect double isotope method and directly through non-local thermodynamic equilibrium (NLTE) 1D radiative transfer modeling of the HCN emission. The latter also provides the ${}^{12}$C/${}^{13}$C ratio, which we compared to the local interstellar value. The derived ${}^{14}$N/${}^{15}$N ratios using the indirect method are generally in the range of 300-550. This result could suggest an evolutionary trend in the nitrogen fractionation of HCN between starless cores and later stages of the star formation process. However, the direct method reveals lower fractionation ratios of around $\sim$250, mainly resulting from a lower ${}^{12}$C/${}^{13}$C ratio in the range $\sim$20-40, as compared to the local interstellar medium value of 68. This study reveals a significant difference between the nitrogen fractionation ratio in HCN derived using direct and indirect methods. This can influence the interpretation of the chemical evolution and reveal the pitfalls of the indirect double isotope method for fractionation studies. However, the direct method is challenging, as it requires well-constrained source models to produce accurate results. No trend in the nitrogen fractionation of HCN between earlier and later stages of the star formation process is evident when the results of the direct method are considered. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: Accepted for publication in A&A

Probing the physics of star formation (ProPStar): I. First resolved maps of the electron fraction and cosmic-ray ionization rate in NGC 1333

Authors: Jaime E. Pineda, Olli Sipilä, Dominique M. Segura-Cox, Maria Teresa Valdivia-Mena, Roberto Neri, Michael Kuffmeier, Alexei V. Ivlev, Stella S. R. Offner, Maria Jose Maureira, Paola Caselli, Silvia Spezzano, Nichol Cunningham, Anika Schmiedeke, Mike Chen

Abstract: Electron fraction and cosmic-ray ionization rates (CRIR) in star-forming regions are important quantities in astrochemical modeling and are critical to the degree of coupling between neutrals, ions, and electrons, which regulates the dynamics of the magnetic field. However, these are difficult quantities to estimate. We aim to derive the electron fraction and CRIR maps of an active star-forming re… ▽ More Electron fraction and cosmic-ray ionization rates (CRIR) in star-forming regions are important quantities in astrochemical modeling and are critical to the degree of coupling between neutrals, ions, and electrons, which regulates the dynamics of the magnetic field. However, these are difficult quantities to estimate. We aim to derive the electron fraction and CRIR maps of an active star-forming region. We combined observations of the nearby NGC 1333 star-forming region carried out with the NOEMA interferometer and IRAM 30-m single dish to generate high spatial dynamic range maps of different molecular transitions. We used the DCO$^+$ and H$^{13}$CO$^+$ ratio (in addition to complementary data) to estimate the electron fraction and produce cosmic-ray ionization rate maps. We derived the first large-area electron fraction and CRIR resolved maps in a star-forming region, with typical values of $10^{-6.5}$ and $10^{-16.5}$ s$^{-1}$, respectively. The maps present clear evidence of enhanced values around embedded young stellar objects (YSOs). This provides strong evidence for locally accelerated cosmic rays. We also found a strong enhancement toward the northwest region in the map that might be related either to an interaction with a bubble or to locally generated cosmic rays by YSOs. We used the typical electron fraction and derived a MHD turbulence dissipation scale of 0.054 pc, which could be tested with future observations. We found a higher cosmic-ray ionization rate compared to the canonical value for $N({\rm H_2})=10^{21}-10^{23}$ cm$^{-2}$ of $10^{-17}$ s$^{-1}$ in the region, and it is likely generated by the accreting YSOs. The high value of the electron fraction suggests that new disks will form from gas in the ideal-MHD limit. This indicates that local enhancements of $ζ({\rm H_2})$, due to YSOs, should be taken into account in the analysis of clustered star formation. △ Less

Submitted 25 February, 2024; originally announced February 2024.

Comments: 13 pages, 13 Figures. Accepted to A&A. Abridged abstract. GitHub repo with code available at https://github.com/jpinedaf/NGC1333_NOEMA_Ions

Journal ref: A&A 686, A162 (2024)

Massive clumps in W43-main: Structure formation in an extensively shocked molecular cloud

Authors: Yuxin Lin, Friedrich Wyrowski, Hauyu Baobab Liu, Yan Gong, Olli Sipilä, Andrés F. Izquierdo, Timea Csengeri, Adam Ginsburg, Guang-Xing Li, Silvia Spezzano, Jaime E. Pineda, Silvia Leurini, Paola Caselli, Karl M. Menten

Abstract: W43-main is a massive molecular complex located at the interaction of the Scutum arm and the Galactic bar undergoing starburst activities. We aim to investigate the gas dynamics, in particular, the prevailing shock signatures from the cloud to clump scale and assess the impact of shocks on the formation of dense gas and early-stage cores. We have carried out NOEMA and IRAM-30m observations at 3 mm… ▽ More W43-main is a massive molecular complex located at the interaction of the Scutum arm and the Galactic bar undergoing starburst activities. We aim to investigate the gas dynamics, in particular, the prevailing shock signatures from the cloud to clump scale and assess the impact of shocks on the formation of dense gas and early-stage cores. We have carried out NOEMA and IRAM-30m observations at 3 mm with an angular resolution of $\sim$0.1 pc towards five massive clumps in W43 main. We use CH$_{3}$CCH and H$_{2}$CS lines to trace the extended gas temperature and CH$_{3}$OH lines to probe the volume density of the dense gas ($\gtrsim$10$^{5}$ cm$^{-3}$). The emission of SiO (2-1) is extensive across the region ($\sim$4 pc) and is mostly contained within a low-velocity regime, hinting at a large-scale origin of the shocks. The position-velocity maps of multiple tracers show systematic spatio-kinematic offsets supporting the cloud-cloud collision/merging scenario. We identify an additional extended velocity component in CCH emission, which coincides with one of the velocity components of the larger scale $^{13}$CO (2-1) emission, likely representing an outer, less dense gas layer in the cloud merging process. We find that the V-shaped, asymmetric SiO wings are tightly correlated with localised gas density enhancements, which is direct evidence of dense gas formation and accumulation in shocks. We resolve two categories of NH$_{2}$D cores: ones exhibiting only subsonic to transonic velocity dispersion, and the others with an additional supersonic velocity dispersion. The centroid velocities of the latter cores are correlated with the shock front seen by SiO. The kinematics of the $\sim$0.1 pc NH$_{2}$D cores are heavily imprinted by shock activities, and may represent a population of early-stage cores forming around the shock interface. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Comments: 35 pages, 34 figures (including appendices); Accepted in A&A

Deuterium fractionation in cold dense cores in the low-mass star forming region L1688

Authors: I. V. Petrashkevich, A. F. Punanova, P. Caselli, O. Sipilä, J. E. Pineda, R. K. Friesen, M. G. Korotaeva, A. I. Vasyunin

Abstract: In this work, we study deuterium fractionation in four starless cores in the low-mass star-forming region L1688 in the Ophiuchus molecular cloud. We study how the deuterium fraction ($R_D$) changes with environment, compare deuteration of ions and neutrals, core centre and its envelope, and attempt to reproduce the observed results with a gas-grain chemical model. We chose high and low gas density… ▽ More In this work, we study deuterium fractionation in four starless cores in the low-mass star-forming region L1688 in the Ophiuchus molecular cloud. We study how the deuterium fraction ($R_D$) changes with environment, compare deuteration of ions and neutrals, core centre and its envelope, and attempt to reproduce the observed results with a gas-grain chemical model. We chose high and low gas density tracers to study both core centre and the envelope. With the IRAM 30m antenna, we mapped N$_2$H$^+$(1-0), N$_2$D$^+$(1-0), H$^{13}$CO$^+$ (1-0) and (2-1), DCO$^+$(2-1), and $p$-NH$_2$D(1$_{11}$-1$_{01}$) towards the chosen cores. The missing $p$-NH$_3$ and N$_2$H$^+$(1-0) data were taken from the literature. To measure the molecular hydrogen column density, dust and gas temperature within the cores, we used the Herschel/SPIRE dust continuum emission data, the GAS survey data (ammonia), and the COMPLETE survey data to estimate the upper limit on CO depletion. We present the deuterium fraction maps for three species towards four starless cores. Deuterium fraction of the core envelopes traced by DCO$^+$/H$^{13}$CO$^+$ is one order of magnitude lower ($\sim$0.08) than that of the core central parts traced by the nitrogen-bearing species ($\sim$0.5). Deuterium fraction increases with the gas density as indicated by high deuterium fraction of high gas density tracers and low deuterium fraction of lower gas density tracers and by the decrease of $R_D$ with core radii, consistent with the predictions of the chemical model. Our model results show a good agreement with observations for $R_D$(N$_2$D$^+$/N$_2$H$^+$) and R$_D$(DCO$^+$/HCO$^+$) and underestimate the $R_D$(NH$_2$D/NH$_3$). △ Less

Submitted 10 January, 2024; v1 submitted 9 January, 2024; originally announced January 2024.

Comments: 30 pages with 22 figures, accepted for publication in MNRAS

MSC Class: 70F15 (Primary) ACM Class: J.2

Nuclear spin ratios of deuterated ammonia in prestellar cores. LAsMA observations of H-MM1 and Oph D

Authors: Jorma Harju, Jaime E. Pineda, Olli Sipilä, Paola Caselli, Arnaud Belloche, Friedrich Wyrowski, Wiebke Riedel, Elena Redaelli, Anton I. Vasyunin

Abstract: We determine the ortho/para ratios of NH2D and NHD2 in two dense, starless cores, where their formation is supposed to be dominated by gas-phase reactions, which, in turn, is predicted to result in deviations from the statistical spin ratios. The Large APEX sub-Millimeter Array (LAsMA) multibeam receiver of the Atacama Pathfinder EXperiment (APEX) telescope was used to observe the prestellar cores… ▽ More We determine the ortho/para ratios of NH2D and NHD2 in two dense, starless cores, where their formation is supposed to be dominated by gas-phase reactions, which, in turn, is predicted to result in deviations from the statistical spin ratios. The Large APEX sub-Millimeter Array (LAsMA) multibeam receiver of the Atacama Pathfinder EXperiment (APEX) telescope was used to observe the prestellar cores H-MM1 and Oph D in Ophiuchus in the ground-state lines of ortho and para NH2D and NHD2. The fractional abundances of these molecules were derived employing 3D radiative transfer modelling, using different assumptions about the abundance profiles as functions of density. We also ran gas-grain chemistry models with different scenarios concerning proton or deuteron exchanges and chemical desorption from grains to find out if one of these models can reproduce the observed spin ratios. The observationally deduced ortho/para ratios of NH2D and NHD2 are in both cores within 10% of their statistical values 3 and 2, respectively, and taking 3-sigma limits, deviations from these of about 20% are allowed. Of the chemistry models tested here, the model that assumes proton hop (as opposed to full scrambling) in reactions contributing to ammonia formation, and a constant efficiency of chemical desorption, comes nearest to the observed abundances and spin ratios. The nuclear spin ratios derived here are in contrast with spin-state chemistry models that assume full scrambling in proton donation and hydrogen abstraction reactions leading to deuterated ammonia. The efficiency of chemical desorption influences strongly the predicted abundances of NH3, NH2D, and NHD2, but has a lesser effect on their ortho/para ratios. For these the proton exchange scenario in the gas is decisive. We suggest that this is because of rapid re-processing of ammonia and related cations by gas-phase ion-molecule reactions. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: accepted for publication in Astronomy & Astrophysics

Modelling Deuterated Isotopologues of Methanol toward the Pre-Stellar Core L1544

Authors: W. Riedel, O. Sipilä, E. Redaelli, P. Caselli, A. I. Vasyunin, F. Dulieu, N. Watanabe

Abstract: Aims. We aim to improve a previous model for the prediction of column densities and deuterium fractions of non- and singly deuterated methanol. Thereby, we try to identify crucial chemical and physical parameters, for which the study of deuteration could provide valuable additional constraints. Methods. We employed a gas-grain chemical code to devise a model that is in agreement with the observe… ▽ More Aims. We aim to improve a previous model for the prediction of column densities and deuterium fractions of non- and singly deuterated methanol. Thereby, we try to identify crucial chemical and physical parameters, for which the study of deuteration could provide valuable additional constraints. Methods. We employed a gas-grain chemical code to devise a model that is in agreement with the observed column density and deuterium fraction profiles of the innermost region of the pre-stellar core L1544. For that purpose, we developed a new treatment of reactive desorption, deriving an individual reactive desorption efficiency for every product species in a chemical reaction, that depends on the reaction enthalpy and type of underlying surface. Furthermore, we explored several options to promote the diffusion of hydrogen and deuterium atoms over the surface of interstellar dust grains, in order to increase methanol formation. Results. Our fiducial model employs diffusion by quantum tunneling of hydrogen and deuterium atoms, resulting in CH$_3$OH and CH$_2$DOH column densities that are approximately an order of magnitude lower than the observed values, which improves the results compared to the previous model by a factor 10. The $N$(CH$_2$DOH)/$N$(CH$_3$OH) ratio is reproduced within a factor of 1.2 for the centre and 1.8 for the position of the methanol peak. Given the large uncertainties that chemical models typically have, we consider our predictions to be in agreement with the observations. In general, we conclude that a diffusion process with a high diffusion rate needs to be employed to obtain methanol column densities that are in accordance with the observed values. Also, we find that the introduction of abstraction reactions into the methanol formation scheme suppresses deuteration, when used in combination with a high diffusion rate. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Journal ref: A&A 680, A87 (2023)

Combined model for $\rm ^{15}N$, $\rm ^{13}C$, and spin-state chemistry in molecular clouds

Authors: O. Sipilä, L. Colzi, E. Roueff, P. Caselli, F. Fontani, E. Wirström

Abstract: We present a new gas-grain chemical model for the combined isotopic fractionation of carbon and nitrogen in molecular clouds, in which the isotope chemistry of carbon and nitrogen is coupled with a time-dependent description of spin-state chemistry. We updated the rate coefficients of some isotopic exchange reactions considered in the literature, and present here a set of new exchange reactions in… ▽ More We present a new gas-grain chemical model for the combined isotopic fractionation of carbon and nitrogen in molecular clouds, in which the isotope chemistry of carbon and nitrogen is coupled with a time-dependent description of spin-state chemistry. We updated the rate coefficients of some isotopic exchange reactions considered in the literature, and present here a set of new exchange reactions involving molecules substituted in $\rm ^{13}C$ and $\rm ^{15}N$ simultaneously. We apply the model to a series of zero-dimensional simulations representing a set of physical conditions across a prototypical prestellar core, exploring the deviations of the isotopic abundance ratios in the various molecules from the elemental isotopic ratios as a function of physical conditions and time. We find that the $\rm ^{12}C/^{13}C$ ratio can deviate from the elemental ratio by up to a factor of several depending on the molecule, and that there are highly time-dependent variations in the ratios. The $\rm HCN/H^{13}CN$ ratio, for example, can obtain values of less than 10 depending on the simulation time. The $\rm ^{14}N/^{15}N$ ratios tend to remain close to the assumed elemental ratio within $\sim$ ten per cent, with no clear trends as a function of the physical conditions. Abundance ratios between $\rm ^{13}C$-containing molecules and $\rm ^{13}C$+$\rm ^{15}N$-containing molecules show somewhat increased levels of fractionation due to the newly included exchange reactions, though still remaining within a few tens of per cent of the elemental $\rm ^{14}N/^{15}N$ ratio. Our results imply the existence of gradients in isotopic abundance ratios across prestellar cores, suggesting that detailed simulations are required to interpret observations of isotopically substituted molecules correctly, especially given that the various isotopic forms of a given molecule do not necessarily trace the same gas layers. △ Less

Submitted 17 October, 2023; v1 submitted 5 September, 2023; originally announced September 2023.

Comments: Accepted to A&A; abstract abridged to meet arXiv requirements

Journal ref: A&A 678, A120 (2023)

Similar levels of deuteration in the pre-stellar core L1544 and the protostellar core HH211

Authors: K. Giers, S. Spezzano, P. Caselli, E. Wirström, O. Sipilä, J. E. Pineda, E. Redaelli, C. T. Bop, F. Lique

Abstract: In the centre of pre-stellar cores, deuterium fractionation is enhanced due to the low temperatures and high densities. Therefore, the chemistry of deuterated molecules can be used to study the earliest stages of star formation. We analyse the deuterium fractionation of simple molecules, comparing the level of deuteration in the envelopes of the pre-stellar core L1544 in Taurus and the protostella… ▽ More In the centre of pre-stellar cores, deuterium fractionation is enhanced due to the low temperatures and high densities. Therefore, the chemistry of deuterated molecules can be used to study the earliest stages of star formation. We analyse the deuterium fractionation of simple molecules, comparing the level of deuteration in the envelopes of the pre-stellar core L1544 in Taurus and the protostellar core HH211 in Perseus. We used single-dish observations of CCH, HCN, HNC, HCO$^+$, and their $^{13}$C-, $^{18}$O- and D-bearing isotopologues, detected with the Onsala 20m telescope. We derived the column densities and the deuterium fractions of the molecules. Additionally, we used radiative transfer simulations and results from chemical modelling to reproduce the observed molecular lines. We used new collisional rate coefficients for HNC, HN$^{13}$C, DNC, and DCN that consider the hyperfine structure of these molecules. We find high levels of deuteration for CCH (10%) in both sources, consistent with other carbon chains, and moderate levels for HCN (5-7%) and HNC (8%). The deuterium fraction of HCO$^+$ is enhanced towards HH211, most likely caused by isotope-selective photodissociation of C$^{18}$O. Similar levels of deuteration show that the process is likely equally efficient towards both cores, suggesting that the protostellar envelope still retains the chemical composition of the original pre-stellar core. The fact that the two cores are embedded in different molecular clouds also suggests that environmental conditions do not have a significant effect on the deuteration within dense cores. Radiative transfer modelling shows that it is necessary to include the outer layers of the cores to consider the effects of extended structures. Besides HCO$^+$ observations, HCN observations towards L1544 also require the presence of an outer diffuse layer where the molecules are relatively abundant. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Comments: 27 pages, 17 figures, accepted for publication in A&A

Journal ref: A&A 676, A78 (2023)

First detection of deuterated methylidyne (CD) in the interstellar medium

Authors: Arshia M. Jacob, Karl M. Menten, Friedrich Wyrowski, Olli Sipilä

Abstract: While the abundance of elemental deuterium is relatively low (D/H ~ a few 1E-5), orders of magnitude higher D/H abundance ratios have been found for many interstellar molecules, enhanced by deuterium fractionation. In cold molecular clouds (T < 20K) deuterium fractionation is driven by the H2D+ ion, whereas at higher temperatures (T > 20-30K) gas-phase deuteration is controlled by reactions with C… ▽ More While the abundance of elemental deuterium is relatively low (D/H ~ a few 1E-5), orders of magnitude higher D/H abundance ratios have been found for many interstellar molecules, enhanced by deuterium fractionation. In cold molecular clouds (T < 20K) deuterium fractionation is driven by the H2D+ ion, whereas at higher temperatures (T > 20-30K) gas-phase deuteration is controlled by reactions with CH2D+ and C2HD+. While the role of H2D+ in driving cold interstellar deuterium chemistry is well understood, thanks to observational constraints from direct measurements of H2D+, deuteration stemming from CH2D+ is far less understood, caused by the absence of direct observational constraints of its key ions. Therefore, making use of chemical surrogates is imperative for exploring deuterium chemistry at intermediate temperatures. Formed at an early stage of ion-molecule chemistry, directly from the dissociative recombination of CH3+ (CH2D+), CH (CD) is an ideal tracer for investigating deuterium substitution initiated by reactions with CH2D+. This paper reports the first detection of CD in the interstellar medium, carried out using the APEX 12m telescope toward the widely studied low-mass protostellar system IRAS 16293-2422. Gas-phase chemical models reproducing the observed CD/CH abundance ratio of 0.016 suggests that it reflects `warm deuterium chemistry' (which ensues in moderately warm conditions of the interstellar medium) and illustrates the potential use of the CD/CH ratio in constraining the gas temperatures of the envelope gas clouds it probes. △ Less

Submitted 11 May, 2023; originally announced May 2023.

Comments: Accepted for publication in Astronomy & Astrophysics (14 pages, 11 figures and 7 tables including Appendix)

Journal ref: A&A 675, A69 (2023)

A large (~1 pc) contracting envelope around the prestellar core L1544

Authors: E. Redaelli, A. Chacón-Tanarro, P. Caselli, M. Tafalla, J. E. Pineda, S. Spezzano, O. Sipilä

Abstract: Prestellar cores, the birthplace of Sun-like stars, form from the fragmentation of the filamentary structure that composes molecular clouds, from which they must inherit at least partially the kinematics. Furthermore, when they are on the verge of gravitational collapse, they show signs of subsonic infall motions. How extended these motions are, which depends on how the collapse occurs, remains la… ▽ More Prestellar cores, the birthplace of Sun-like stars, form from the fragmentation of the filamentary structure that composes molecular clouds, from which they must inherit at least partially the kinematics. Furthermore, when they are on the verge of gravitational collapse, they show signs of subsonic infall motions. How extended these motions are, which depends on how the collapse occurs, remains largely unknown. We want to investigate the kinematics of the envelope that surrounds the prototypical prestellar core L1544, studying the cloud-core connection. To our aims, we observed the $\rm HCO^+$(1-0) transition in a large map. \hcop is expected to be abundant in the envelope, making it an ideal probe of the large-scale kinematics in the source. We modelled the spectrum at the dust peak by means of a non local-thermodynamical-equilibrium radiative transfer. In order to reproduce the spectrum at the dust peak, a large ($\sim 1\, \rm pc$) envelope is needed, with low density (tens of $\rm cm^{-3}$ at most) and contraction motions, with an inward velocity of $\approx 0.05\,\rm km \, s^{-1}$. We fitted the data cube using the Hill5 model, which implements a simple model {for the optical depth and excitation temperature profiles along the line-of-sight,} in order to obtain a map of the infall velocity. This shows that the infall motions are extended, with typical values in the range $0.1-0.2\,\rm km \, s^{-1}$. Our results suggest that the contraction motions extend in the diffuse envelope surrounding the core, which is consistent with recent magnetic field measurements in the source, which showed that the envelope is magnetically supercritical. △ Less

Submitted 24 October, 2022; originally announced October 2022.

Comments: Accepted for publication on ApJ, 24 Oct. 2022

Chemistry and dynamics of the prestellar core L1544

Authors: O. Sipilä, P. Caselli, E. Redaelli, S. Spezzano

Abstract: We aim to quantify the effect of chemistry on the infall velocity in the prestellar core L1544. Previous observational studies have found evidence for double-peaked line profiles for the rotational transitions of several molecules, which cannot be accounted for with the models presently available for the physical structure of the source, without ad hoc up-scaling of the infall velocity. We ran one… ▽ More We aim to quantify the effect of chemistry on the infall velocity in the prestellar core L1544. Previous observational studies have found evidence for double-peaked line profiles for the rotational transitions of several molecules, which cannot be accounted for with the models presently available for the physical structure of the source, without ad hoc up-scaling of the infall velocity. We ran one-dimensional hydrodynamical simulations of the collapse of a core with L1544-like properties (in terms of mass and outer radius), using a state-of-the-art chemical model with a very large chemical network combined with an extensive description of molecular line cooling, determined via radiative transfer simulations, with the aim of determining whether these expansions of the simulation setup (as compared to previous models) can lead to a higher infall velocity. After running a series of simulations where the simulation was sequentially simplified, we found that the infall velocity is almost independent of the size of the chemical network or the approach to line cooling. We conclude that chemical evolution does not have a large impact on the infall velocity, and that the higher infall velocities that are implied by observations may be the result of the core being more dynamically evolved than what is now thought, or alternatively the average density in the simulated core is too low. However, chemistry does have a large influence on the lifetime of the core, which varies by about a factor of two across the simulations and grows longer when the chemical network is simplified. Therefore, although the model is subject to several sources of uncertainties, the present results clearly indicate that the use of a small chemical network leads to an incorrect estimate of the core lifetime, which is naturally a critical parameter for the development of chemical complexity in the precollapse phase. △ Less

Submitted 21 October, 2022; v1 submitted 28 September, 2022; originally announced September 2022.

Comments: 12 pages incl. Appendix; 8 figures; accepted to A&A; replacement to reflect the change of the title of the paper post-acceptance

Journal ref: A&A 668, A131 (2022)

Tracing the contraction of the pre-stellar core L1544 with HC$^{17}$O$^+$ $J$ = 1-0 emission

Authors: J. Ferrer Asensio, S. Spezzano, P. Caselli, F. O. Alves, O. Sipilä, E. Redaelli, L. Bizzocchi, F. Lique, A. Mullins

Abstract: Spectral line profiles of several molecules observed towards the pre-stellar core L1544 appear double-peaked. For abundant molecular species this line morphology has been linked to self-absorption. However, the physical process behind the double-peaked morphology for less abundant species is still under debate. In order to understand the cause behind the double-peaked spectra of optically thin tra… ▽ More Spectral line profiles of several molecules observed towards the pre-stellar core L1544 appear double-peaked. For abundant molecular species this line morphology has been linked to self-absorption. However, the physical process behind the double-peaked morphology for less abundant species is still under debate. In order to understand the cause behind the double-peaked spectra of optically thin transitions and their link to the physical structure of pre-stellar cores, we present high-sensitivity and high-spectral resolution HC$^{17}$O$^+$ $J =$1-0 observations towards the dust peak in L1544. We observed the HC$^{17}$O$^+$ (1-0) spectrum with the Institut de Radioastronomie Millimétrique (IRAM) 30m telescope. By using new state-of-the-art collisional rate coefficients, a physical model for the core and the fractional abundance profile of HC$^{17}$O$^+$, the hyperfine structure of this molecular ion is modelled for the first time with the radiative transfer code LOC applied to the predicted chemical structure of a contracting pre-stellar core. We applied the same analysis to the chemically related C$^{17}$O molecule. The observed HC$^{17}$O$^+$(1-0) and C$^{17}$O(1-0) lines have been successfully reproduced with a non-local thermal equilibrium (LTE) radiative transfer model applied to chemical model predictions for a contracting pre-stellar core. An upscaled velocity profile (by 30%) is needed to reproduce the HC$^{17}$O$^+$(1-0) observations. The double peaks observed in the HC$^{17}$O$^+$(1-0) hyperfine components are due to the contraction motions at densities close to the critical density of the transition ($\sim$10$^{5}$ cm$^{-3}$) and to the fact that the HCO$^{+}$ fractional abundance decreases toward the centre. △ Less

Submitted 6 September, 2022; originally announced September 2022.

A Survey of Deuterated Ammonia in the Cepheus Star-Forming Region L1251

Authors: Maria Galloway-Sprietsma, Yancy L. Shirley, James Di Francesco, Jared Keown, Samantha Scibelli, Olli Sipilä, Rachel Smullen

Abstract: Understanding the chemical processes during starless core and prestellar core evolution is an important step in understanding the initial stages of star and disk formation. This project is a study of deuterated ammonia, o-NH$_2$D, in the L1251 star-forming region toward Cepheus. Twenty-two dense cores (twenty of which are starless or prestellar, and two of which have a protostar), previously ident… ▽ More Understanding the chemical processes during starless core and prestellar core evolution is an important step in understanding the initial stages of star and disk formation. This project is a study of deuterated ammonia, o-NH$_2$D, in the L1251 star-forming region toward Cepheus. Twenty-two dense cores (twenty of which are starless or prestellar, and two of which have a protostar), previously identified by p-NH$_3$ (1,1) observations, were targeted with the 12m Arizona Radio Observatory telescope on Kitt Peak. o-NH$_2$D J$_{\rm{K_a} \rm{K_c}}^{\pm} =$ $1_{11}^{+} \rightarrow 1_{01}^{-}$ was detected in 13 (59\%) of the NH$_3$-detected cores with a median sensitivity of $σ_{T_{mb}} = 17$ mK. All cores detected in o-NH$_2$D at this sensitivity have p-NH$_3$ column densities $> 10^{14}$ cm$^{-2}$. The o-NH$_2$D column densities were calculated using the constant excitation temperature (CTEX) approximation while correcting for the filling fraction of the NH$_3$ source size. The median deuterium fraction was found to be 0.11 (including 3$σ$ upper limits). However, there are no strong, discernible trends in plots of deuterium fraction with any physical or evolutionary variables. If the cores in L1251 have similar initial chemical conditions, then this result is evidence of the cores physically evolving at different rates. △ Less

Submitted 21 July, 2022; originally announced July 2022.

Comments: 18 pages, 14 figures, MNRAS, in press

Nitrogen fractionation towards a pre-stellar core traces isotope-selective photodissociation

Authors: Silvia Spezzano, Paola Caselli, Olli Sipilä, Luca Bizzocchi

Abstract: Isotopologue abundance ratios are important to understand the evolution of astrophysical objects and ultimately the origins of a planetary system like our own. Being nitrogen a fundamental ingredient of pre-biotic material, understanding its chemistry and inheritance is of fundamental importance to understand the formation of the building blocks of life. We present here single-dish observations of… ▽ More Isotopologue abundance ratios are important to understand the evolution of astrophysical objects and ultimately the origins of a planetary system like our own. Being nitrogen a fundamental ingredient of pre-biotic material, understanding its chemistry and inheritance is of fundamental importance to understand the formation of the building blocks of life. We present here single-dish observations of the ground state rotational transitions of the $^{13}$C and $^{15}$N isotopologues of HCN, HNC and CN with the IRAM 30m telescope. We analyse their column densities and compute the $^{14}$N/$^{15}$N ratio map for HCN. The $^{15}$N-fractionation of CN and HNC is computed towards different offsets across L1544. The $^{15}$N-fractionation map of HCN shows a clear decrease of the $^{14}$N/$^{15}$N ratio towards the southern edge of L1544, where carbon chain molecules present a peak, strongly suggesting that isotope-selective photodissociation has a strong effect on the fractionation of nitrogen across pre-stellar cores. The $^{14}$N/$^{15}$N ratio in CN measured towards four positions across the core also shows a decrease towards the South-East of the core, while HNC shows opposite behaviour. The uneven illumination of the pre-stellar core L1544 provides clear evidence that $^{15}$N-fractionation of HCN and CN is enhanced toward the region more exposed to the interstellar radiation field. Isotope-selective photodissociation of N$_2$ is then a crucial process to understand $^{15}$N fractionation, as already found in protoplanetary disks. Therefore, the $^{15}$N-fractionation in pre-stellar material is expected to change depending on the environment within which pre-stellar cores are embedded. The $^{12}$CN/$^{13}$CN ratio also varies across the core, but its variation does not affect our conclusions on the effect of the environment on the fractionation of nitrogen. △ Less

Submitted 13 July, 2022; originally announced July 2022.

Comments: Accepted for publication in A&A Letters

Journal ref: A&A 664, L2 (2022)

Deuteration of c-C$_3$H$_2$ towards the pre-stellar core L1544

Authors: K. Giers, S. Spezzano, F. Alves, P. Caselli, E. Redaelli, O. Sipilä, M. Ben Khalifa, L. Wiesenfeld, S. Brünken, L. Bizzocchi

Abstract: Context: In the centre of pre-stellar cores, the deuterium fractionation is enhanced due to the cold temperatures and high densities. Therefore, the chemistry of deuterated molecules can be used to probe the evolution and the kinematics in the earliest stages of star formation. Aims: We analyse emission maps of cyclopropenylidene, c-C$_3$H$_2$, to study the distribution of the deuteration througho… ▽ More Context: In the centre of pre-stellar cores, the deuterium fractionation is enhanced due to the cold temperatures and high densities. Therefore, the chemistry of deuterated molecules can be used to probe the evolution and the kinematics in the earliest stages of star formation. Aims: We analyse emission maps of cyclopropenylidene, c-C$_3$H$_2$, to study the distribution of the deuteration throughout the prototypical pre-stellar core L1544. Methods: We use single-dish observations of c-C$_3$H$_2$, c-H$^{13}$CC$_2$H, c-C$_3$HD, and c-C$_3$D$_2$ towards the pre-stellar core L1544, performed at the IRAM 30m telescope. We derive the column density and deuterium fraction maps, and compare these observations with non-LTE radiative transfer simulations. Results: The highest deuterium fractions are found close to the dust peak at the centre of L1544, where the increased abundance of H$_2$D$^+$ ions drives the deuteration process. The peak values are N(c-C$_3$HD)/N(c-C$_3$H$_2)=0.17\pm0.01$, N(c-C$_3$D$_2$)/N(c-C$_3$H$_2)=0.025\pm0.003$ and N(c-C$_3$D$_2$)/N(c-C$_3$HD$)=0.16\pm0.03$, which is consistent with previous single point observations. The distributions of c-C$_3$HD and c-C$_3$D$_2$ indicate that the deuterated forms of c-C$_3$H$_2$ in fact trace the dust peak and not the c-C$_3$H$_2$ peak. Conclusions: The N(c-C$_3$D$_2$)/N(c-C$_3$HD) map confirms that the process of deuteration is more efficient towards the centre of the core and demonstrates that carbon-chain molecules are still present at high densities. This is likely caused by an increased abundance of He$^+$ ions destroying CO, which increases the amount of carbon atoms in the gas phase. △ Less

Submitted 30 May, 2022; originally announced May 2022.

Comments: 13 pages, 9 figures

Journal ref: A&A 664, A119 (2022)

Multi-line observations of CH$_{3}$OH, c-C$_{3}$H$_{2}$ and HNCO towards L1544: Dissecting the core structure with chemical differentiation

Authors: Yuxin Lin, Silvia Spezzano, Olli Sipilä, Anton Vasyunin, Paola Caselli

Abstract: Pre-stellar cores are the basic unit for the formation of stars and stellar systems. The anatomy of the physical and chemical structures of pre-stellar cores is critical for understanding the star formation process. L1544 is a prototypical pre-stellar core, which shows significant chemical differentiation surrounding the dust peak. We aim to constrain the physical conditions at the different molec… ▽ More Pre-stellar cores are the basic unit for the formation of stars and stellar systems. The anatomy of the physical and chemical structures of pre-stellar cores is critical for understanding the star formation process. L1544 is a prototypical pre-stellar core, which shows significant chemical differentiation surrounding the dust peak. We aim to constrain the physical conditions at the different molecular emission peaks. This study allows us to compare the abundance profiles predicted from chemical models together with the classical density structure of Bonnor-Ebert (BE) sphere. We conducted multi-transition pointed observations of CH$_{3}$OH, c-C$_{3}$H$_{2}$ and HNCO with the IRAM 30m telescope, towards the dust peak and the respective molecular peaks of L1544. With non-LTE radiative transfer calculations and a 1-dimensional model, we revisit the physical structure of L1544, and benchmark with the abundance profiles from current chemical models. We find that the HNCO, c-C$_{3}$H$_{2}$ and CH$_{3}$OH lines in L1544 are tracing progressively higher density gas, from $\sim$10$^{4}$ to several times 10$^{5}$ cm$^{-3}$. Particularly, we find that to produce the observed intensities and ratios of the CH$_{3}$OH lines, a local gas density enhancement upon the BE sphere is required. This suggests that the physical structure of an early-stage core may not necessarily follow a smooth decrease of gas density profile locally, but can be intercepted by clumpy substructures surrounding the gravitational center. Multiple transitions of molecular lines from different molecular species can provide a tomographic view of the density structure of pre-stellar cores. The local gas density enhancement deviating from the BE sphere may reflect the impact of accretion flows that appear asymmetric and are enhanced at the meeting point of large-scale cloud structures. △ Less

Submitted 19 May, 2022; originally announced May 2022.

Comments: accepted by A&A; 22 pages, 22 figures incl. appendices

Journal ref: A&A 665, A131 (2022)

An Interferometric View of H-MM1. I. Direct Observation of NH3 Depletion

Authors: Jaime E. Pineda, Jorma Harju, Paola Caselli, Olli Sipilä, Mika Juvela, Charlotte Vastel, Erik Rosolowsky, Andreas Burkert, Rachel K. Friesen, Yancy Shirley, María José Maureira, Spandan Choudhury, Dominique M. Segura-Cox, Rolf Güsten, Anna Punanova, Luca Bizzocchi, Alyssa A. Goodman

Abstract: Spectral lines of ammonia, NH$_3$, are useful probes of the physical conditions in dense molecular cloud cores. In addition to advantages in spectroscopy, ammonia has also been suggested to be resistant to freezing onto grain surfaces, which should make it a superior tool for studying the interior parts of cold, dense cores. Here we present high-resolution NH$_3$ observations with the Very Large A… ▽ More Spectral lines of ammonia, NH$_3$, are useful probes of the physical conditions in dense molecular cloud cores. In addition to advantages in spectroscopy, ammonia has also been suggested to be resistant to freezing onto grain surfaces, which should make it a superior tool for studying the interior parts of cold, dense cores. Here we present high-resolution NH$_3$ observations with the Very Large Array (VLA) and Green Bank Telescope (GBT) towards a prestellar core. These observations show an outer region with a fractional NH$_3$ abundance of X(NH$_3$) = (1.975$\pm$0.005)$\times 10^{-8}$ ($\pm 10\%$ systematic), but it also reveals that after all, the X(NH$_3$) starts to decrease above a H$_2$ column density of $\approx 2.6 \times 10^{22}$ cm$^{-2}$. We derive a density model for the core and find that the break-point in the fractional abundance occurs at the density n(H$_2$) $\sim 2\times10^5$ cm$^{-3}$, and beyond this point the fractional abundance decreases with increasing density, following the power law $n^{-1.1}$. This power-law behavior is well reproduced by chemical models where adsorption onto grains dominates the removal of ammonia and related species from the gas at high densities. We suggest that the break-point density changes from core to core depending on the temperature and the grain properties, but that the depletion power law is anyway likely to be close to $n^{-1}$ owing to the dominance of accretion in the central parts of starless cores. △ Less

Submitted 2 May, 2022; originally announced May 2022.

Comments: Accepted for publication in AJ. 18 Pages, 16 Figures, 3 Tables

H2CS deuteration maps towards the pre-stellar core L1544

Authors: S. Spezzano, O. Sipilä, P. Caselli, S. S. Jensen, S. Czakli, L. Bizzocchi, J. Chantzos, G. Esplugues, A. Fuente, F. Eisenhauer

Abstract: Deuteration is a crucial tool to understand the complexity of interstellar chemical processes, especially when they involve the interplay of gas-phase and grain-surface chemistry. In the case of multiple deuteration, comparing observation with the results of chemical modelling is particularly effective to study how molecules are inherited in the different stages within the process of star and plan… ▽ More Deuteration is a crucial tool to understand the complexity of interstellar chemical processes, especially when they involve the interplay of gas-phase and grain-surface chemistry. In the case of multiple deuteration, comparing observation with the results of chemical modelling is particularly effective to study how molecules are inherited in the different stages within the process of star and planet formation. We aim to study the the D/ H ratio in H2CS across the prototypical pre-stellar core L1544. This study allows us to test current gas-dust chemical models involving sulfur in dense cores. We present here single-dish observations of H2CS, HDCS and D2CS with the IRAM 30m telescope. We analyse their column densities and distributions, and compare these observations with gas-grain chemical models. The deuteration maps of H2CS in L1544 are compared with the deuteration maps of methanol, H2CO, N2H+ and HCO+ towards the same source. Furthermore, the single and double deuteration of H2CS towards the dust peak of L1544 is compared with H2CO and c-C3H2. The difference between the deuteration of these molecules in L1544 is discussed and compared with the prediction of chemical models. The maximum deuterium fractionation for the first deuteration of H2CS is N(HDCS)/N(H2CS)$\sim$30$\%$ and is located towards the north-east at a distance of about 10000 AU from the dust peak. While for c-C3H2 the first and second deuteration have a similar efficiency, for H2CS and H2CO the second deuteration is more efficient, leading to D2CX/HDCX$\sim$100$\%$ (with X= O or S). Our results imply that the large deuteration of H2CO and H2CS observed in protostellar cores as well as in comets is likely inherited from the pre-stellar phase. However, the comparison with state-of-the-art chemical models suggests that the reaction network for the formation of the doubly deuterated H2CS and H2CO it is not complete yet. △ Less

Submitted 15 March, 2022; originally announced March 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 661, A111 (2022)

The Central 1000 au of a Pre-stellar Core Revealed with ALMA. II. Almost Complete Freeze-out

Authors: Paola Caselli, Jaime E. Pineda, Olli Sipilä, Bo Zhao, Elena Redaelli, Silvia Spezzano, Maria José Maureira, Felipe Alves, Luca Bizzocchi, Tyler L. Bourke, Ana Chacón-Tanarro, Rachel Friesen, Daniele Galli, Jorma Harju, Izaskun Jiménez-Serra, Eric Keto, Zhi-Yun Li, Marco Padovani, Anika Schmiedeke, Mario Tafalla, Charlotte Vastel

Abstract: Pre-stellar cores represent the initial conditions in the process of star and planet formation. Their low temperatures ($<$10 K) allow the formation of thick icy dust mantles, which will be partially preserved in the future protoplanetary disks, ultimately affecting the chemical composition of planetary systems. Previous observations have shown that carbon- and oxygen-bearing species, in particula… ▽ More Pre-stellar cores represent the initial conditions in the process of star and planet formation. Their low temperatures ($<$10 K) allow the formation of thick icy dust mantles, which will be partially preserved in the future protoplanetary disks, ultimately affecting the chemical composition of planetary systems. Previous observations have shown that carbon- and oxygen-bearing species, in particular CO, are heavily depleted in pre-stellar cores due to the efficient molecular freeze-out onto the surface of cold dust grains. However, N-bearing species such as NH$_3$ and, in particular, its deuterated isotopologues, appear to maintain high abundances where CO molecules are mainly in solid phase. Thanks to ALMA, we present here the first clear observational evidence of NH$_2$D freeze-out toward the L1544 pre-stellar core, suggestive of the presence of a"complete-depletion zone" within a $\simeq$1800 au radius, in agreement with astrochemical pre-stellar core model predictions. Our state-of-the-art chemical model coupled with a non-LTE radiative transfer code demonstrates that NH$_2$D becomes mainly incorporated in icy mantles in the central 2000 au and starts freezing-out already at $\simeq$7000 au. Radiative transfer effects within the pre-stellar core cause the NH$_2$D(1$_{11}$-1$_{01}$) emission to appear centrally concentrated, with a flattened distribution within the central $\simeq$3000 au, unlike the 1.3 mm dust continuum emission which shows a clear peak within the central $\simeq$1800 au. This prevented NH$_2$D freeze-out to be detected in previous observations, where the central 1000 au cannot be spatially resolved. △ Less

Submitted 27 February, 2022; originally announced February 2022.

Comments: accepted for publication in The Astrophysical Journal

The cosmic-ray ionisation rate in the pre-stellar core L1544

Authors: Elena Redaelli, Olli Sipilä, Marco Padovani, Paola Caselli, Daniele Galli, Alexei V. Ivlev

Abstract: Context. Cosmic rays (CRs) play an important role in the chemistry and dynamics of the interstellar medium. In dense environments, they represent the main ionising agent, driving the rich chemistry of molecular ions and determining the ionisation fraction, which regulates the degree of coupling between the gas and magnetic fields. Estimates of the CR ionisation rate ($ζ_2$) span several orders of… ▽ More Context. Cosmic rays (CRs) play an important role in the chemistry and dynamics of the interstellar medium. In dense environments, they represent the main ionising agent, driving the rich chemistry of molecular ions and determining the ionisation fraction, which regulates the degree of coupling between the gas and magnetic fields. Estimates of the CR ionisation rate ($ζ_2$) span several orders of magnitude, depending on the targeted sources and on the used method. Aims. Recent theoretical models have characterised the CR attenuation with increasing density. We aim to test these models for the attenuation of CRs in the low-mass pre-stellar core L1544. Methods. We use a state-of-the-art gas-grain chemical model, which accepts the CR ionisation rate profile as input, to predict the abundance profiles of four ions: $\rm N_2H^+$, $\rm N_2D^+$, $\rm HC^{18}O^+$, and $\rm DCO^+$. Non-LTE radiative transfer is performed to produce synthetic spectra based on the derived abundances. These are compared with observations obtained with the Institut de Radioastronomie Millimétrique (IRAM) 30m telescope. Results. Our results indicate that a model with $ζ_2 > 10^{-16} \rm \, s^{-1}$ is excluded by the observations. Also the model with the standard $ζ_2 = 1.3 \times 10^{-17} \rm \, s^{-1}$ produces a worse agreement with respect to the attenuation model based on Voyager observations, which has an average $ζ_2 = 3 \times 10^{-17} \rm \, s^{-1}$ at the column densities typical of L1544. The single-dish data, however, are not sensitive to the attenuation of the CR profile, which changes only by a factor of two in the range of column densities spanned by the core model. Interferometric observations at higher spatial resolution, combined with observations of transitions with lower critical density are needed to observe a decrease of $ζ_2$ with density. △ Less

Submitted 16 September, 2021; originally announced September 2021.

Comments: Accepted in A&A on 15 Sept. 2021

Journal ref: A&A 656, A109 (2021)

A Revised Description of the Cosmic Ray-Induced Desorption of Interstellar Ices

Authors: O. Sipilä, K. Silsbee, P. Caselli

Abstract: Non-thermal desorption of ices on interstellar grains is required to explain observations of molecules that are not synthesized efficiently in the gas phase in cold dense clouds. Perhaps the most important non-thermal desorption mechanism is one induced by cosmic rays (CRs), which, when passing through a grain, heat it transiently to a high temperature - the grain cools back to its original equili… ▽ More Non-thermal desorption of ices on interstellar grains is required to explain observations of molecules that are not synthesized efficiently in the gas phase in cold dense clouds. Perhaps the most important non-thermal desorption mechanism is one induced by cosmic rays (CRs), which, when passing through a grain, heat it transiently to a high temperature - the grain cools back to its original equilibrium temperature via the (partial) sublimation of the ice. Current cosmic-ray-induced desorption (CRD) models assume a fixed grain cooling time. In this work we present a revised description of CRD in which the desorption efficiency depends dynamically on the ice content. We apply the revised desorption scheme to two-phase and three-phase chemical models in physical conditions corresponding to starless and prestellar cores, and to molecular cloud envelopes. We find that inside starless and prestellar cores, introducing dynamic CRD can decrease gas-phase abundances by up to an order of magnitude in two-phase chemical models. In three-phase chemical models our model produces very similar results to the static cooling scheme - when only one monolayer of ice is considered active. Ice abundances are generally insensitive to variations in the grain cooling time. Further improved CRD models need to take into account additional effects in the transient heating of the grains, introduced for example by the adoption of a spectrum of CR energies. △ Less

Submitted 3 September, 2021; v1 submitted 8 June, 2021; originally announced June 2021.

Comments: Accepted to ApJ; note significant changes to the model results compared to the version of the paper previously posted on arXiv, especially for low volume densities and for the three-phase ice models, induced by the addition of an upper limit to the CRD rate

First survey of HCNH$^+$ in high-mass star-forming cloud cores

Authors: F. Fontani, L. Colzi, E. Redaelli, O. Sipilä, P. Caselli

Abstract: Most stars in the Galaxy, including the Sun, were born in high-mass star-forming regions. It is hence important to study the chemical processes in these regions to better understand the chemical heritage of both the Solar System and most stellar systems in the Galaxy. The molecular ion HCNH+ is thought to be a crucial species in ion-neutral astrochemical reactions, but so far it has been detected… ▽ More Most stars in the Galaxy, including the Sun, were born in high-mass star-forming regions. It is hence important to study the chemical processes in these regions to better understand the chemical heritage of both the Solar System and most stellar systems in the Galaxy. The molecular ion HCNH+ is thought to be a crucial species in ion-neutral astrochemical reactions, but so far it has been detected only in a handful of star-forming regions, and hence its chemistry is poorly known. We have observed with the IRAM-30m Telescope 26 high-mass star-forming cores in different evolutionary stages in the J=3-2 rotational transition of HCNH+. We report the detection of HCNH+ in 16 out of 26 targets. This represents the largest sample of sources detected in this molecular ion so far. The fractional abundances of HCNH+, [HCNH+], w.r.t. H2, are in the range 0.9 - 14 X $10^{-11}$, and the highest values are found towards cold starless cores. The abundance ratios [HCNH+]/[HCN] and [HCNH+]/[HCO+] are both < 0.01 for all objects except for four starless cores, for which they are well above this threshold. These sources have the lowest gas temperature in the sample. We run two chemical models, a "cold" one and a "warm" one, which attempt to match as much as possible the average physical properties of the cold(er) starless cores and of the warm(er) targets. The reactions occurring in the latter case are investigated in this work for the first time. Our predictions indicate that in the warm model HCNH+ is mainly produced by reactions with HCN and HCO+, while in the cold one the main progenitor species of HCNH+ are HCN+ and HNC+. The results indicate that the chemistry of HCNH+ is different in cold/early and warm/evolved cores, and the abundance ratios [HCNH+]/[HCN] and [HCNH+]/[HCO+] is a useful astrochemical tool to discriminate between different evolutionary phases in the process of star formation. △ Less

Submitted 18 May, 2021; originally announced May 2021.

Comments: 13 pages, 3 tables, 12 figures (+3 figures in appendix). Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 651, A94 (2021)

The chemical structure of young high-mass star-forming clumps: (II) parsec-scale CO depletion and deuterium fraction of $\rm HCO^+$

Authors: S. Feng, D. Li, P. Caselli, F. Du, Y. Lin, O. Sipilä, H. Beuther, Patricio Sanhueza, K. Tatematsu, S. Y. Liu, Q. Zhang, Y. Wang, T. Hogge, I. Jimenez-Serra, X. Lu, T. Liu, K. Wang, Z. Y. Zhang, S. Zahorecz, G. Li, H. B. Liu, J. Yuan

Abstract: The physical and chemical properties of cold and dense molecular clouds are key to understanding how stars form. Using the IRAM 30 m and NRO 45 m telescopes, we carried out a Multiwavelength line-Imaging survey of the 70 $μ$m dark and bright clOuds (MIAO). At a linear resolution of 0.1--0.5 pc, this work presents a detailed study of parsec-scale CO depletion and $\rm HCO^+$ deuterium (D-) fraction… ▽ More The physical and chemical properties of cold and dense molecular clouds are key to understanding how stars form. Using the IRAM 30 m and NRO 45 m telescopes, we carried out a Multiwavelength line-Imaging survey of the 70 $μ$m dark and bright clOuds (MIAO). At a linear resolution of 0.1--0.5 pc, this work presents a detailed study of parsec-scale CO depletion and $\rm HCO^+$ deuterium (D-) fractionation toward four sources (G11.38+0.81, G15.22-0.43, G14.49-0.13, and G34.74-0.12) included in our full sample. In each source with $\rm T<20$ K and $n_{\rm H}\rm\sim10^4$--$\rm 10^5 cm^{-3}$, we compared pairs of neighboring 70 $μ$m bright and dark clumps and found that (1) the $\rm H_2$ column density and dust temperature of each source show strong spatial anticorrelation; (2) the spatial distribution of CO isotopologue lines and dense gas tracers, such as 1--0 lines of $\rm H^{13}CO^+$ and $\rm DCO^+$, are anticorrelated; (3) the abundance ratio between $\rm C^{18}O$ and $\rm DCO^+$ shows a strong correlation with the source temperature; (4) both the $\rm C^{18}O$ depletion factor and D-fraction of $\rm HCO^+$ show a robust decrease from younger clumps to more evolved clumps by a factor of more than 3; and (5) preliminary chemical modeling indicates chemical ages of our sources are ${\sim}8\times10^4$ yr, which is comparable to their free-fall timescales and smaller than their contraction timescales, indicating that our sources are likely dynamically and chemically young. △ Less

Submitted 27 September, 2020; v1 submitted 8 August, 2020; originally announced August 2020.

Comments: accepted by ApJ

First detection of NHD and ND$_2$ in the interstellar medium

Authors: M. Melosso, L. Bizzocchi, O. Sipilä, B. M. Giuliano, L. Dore, F. Tamassia, M. -A. Martin-Drumel, O. Pirali, E. Redaelli, P. Caselli

Abstract: Deuterium fractionation processes in the interstellar medium (ISM) have been shown to be highly efficient in the family of nitrogen hydrides. To date, observations were limited to ammonia (NH$_2$D, NHD$_2$, ND$_3$) and imidogen radical (ND) isotopologues. We want to explore the high frequency windows offered by the \emph{Herschel Space Observatory} to search for deuterated forms of amidogen radica… ▽ More Deuterium fractionation processes in the interstellar medium (ISM) have been shown to be highly efficient in the family of nitrogen hydrides. To date, observations were limited to ammonia (NH$_2$D, NHD$_2$, ND$_3$) and imidogen radical (ND) isotopologues. We want to explore the high frequency windows offered by the \emph{Herschel Space Observatory} to search for deuterated forms of amidogen radical NH$_2$ and to compare the observations against the predictions of our comprehensive gas-grain chemical model. Making use of the new molecular spectroscopy data recently obtained at high frequencies for NHD and ND$_2$, both isotopologues have been searched for in the spectral survey towards the class 0 IRAS 16293-2422, a source in which NH$_3$, NH and their deuterated variants have been previously detected. We used the observations carried out with HIFI (Heterodyne Instrument for the Far Infrared) in the framework of the key program "Chemical Herschel surveys of star forming regions" (CHESS). We report the first detection of interstellar NHD and ND$_2$. Both species are observed in absorption against the continuum of the protostar. From the analysis of their hyperfine structure, accurate excitation temperature and column density values have been determined. The latter were combined with the column density of the parent species NH$_2$ to derive the deuterium fractionation in amidogen. The amidogen D/H ratio measured in the low-mass protostar IRAS 16293-2422 is comparable to the one derived for the related species imidogen and much higher than that observed for ammonia. Additional observations of these species will give more insights into the mechanism of ammonia formation and deuteration in the ISM. We finally indicate the current possibilities to further explore these species at submillimeter wavelengths. △ Less

Submitted 15 July, 2020; originally announced July 2020.

Comments: 11 pages, 5 figures, 7 tables. Accepted for publication in A&A

Journal ref: A&A 641, A153 (2020)

Rapid Elimination of Small Dust Grains in Molecular Clouds

Authors: Kedron Silsbee, Alexei Ivlev, Paola Caselli, Olli Sipila, Bo Zhao

Abstract: We argue that impact velocities between dust grains with sizes less than $\sim 0.1$ $μm$ in molecular cloud cores are dominated by drift arising from ambipolar diffusion. This effect is due to the size dependence of the dust coupling to the magnetic field and the neutral gas. Assuming perfect sticking in collisions up to $\approx 50$ m/s, we show that this effect causes rapid depletion of small gr… ▽ More We argue that impact velocities between dust grains with sizes less than $\sim 0.1$ $μm$ in molecular cloud cores are dominated by drift arising from ambipolar diffusion. This effect is due to the size dependence of the dust coupling to the magnetic field and the neutral gas. Assuming perfect sticking in collisions up to $\approx 50$ m/s, we show that this effect causes rapid depletion of small grains - consistent with starlight extinction and IR/microwave emission measurements, both in the core center ($n \sim 10^{6}$ cm$^{-3}$) and envelope ($n \sim 10^{4}$ cm$^{-3}$). The upper end of the size distribution does not change significantly if only velocities arising from this effect are considered. We consider the impact of an evolved dust size distribution on the gas temperature, and argue that if the depletion of small dust grains occurs as would be expected from our model, then the cosmic ray ionization rate must be well below $10^{-16}$ s$^{-1}$ at a number density of $10^{5}$ cm$^{-3}$. △ Less

Submitted 29 September, 2020; v1 submitted 23 June, 2020; originally announced June 2020.

Journal ref: A&A 641, A39 (2020)

Effect of grain size distribution and size-dependent grain heating on molecular abundances in starless and pre-stellar cores

Authors: O. Sipilä, B. Zhao, P. Caselli

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

Submitted 17 June, 2020; originally announced June 2020.

Comments: Accepted by A&A; 16 pages incl. appendices; abstract abridged to meet ArXiV requirements

Journal ref: A&A 640, A94 (2020)

Carbon isotopic fractionation in molecular clouds

Authors: L. Colzi, O. Sipilä, E. Roueff, P. Caselli, F. Fontani

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

Submitted 2 October, 2020; v1 submitted 5 June, 2020; originally announced June 2020.

Comments: The published exponent in the A&A journal of the second rate coefficient of Table 1 is a misprint

Journal ref: A&A 640, A51 (2020)

Search for H$_3^+$ isotopologues toward CRL 2136 IRS 1

Authors: Miwa Goto, T. R. Geballe, Jorma Harju, Paola Caselli, Olli Sipilä, Karl M. Menten, Tomonori Usuda

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

Submitted 9 October, 2019; originally announced October 2019.

Comments: 9 pages, 3 tables, 8 figures, accepted for publication at A&A

Journal ref: A&A 632, A29 (2019)

Modeling deuterium chemistry in starless cores: full scrambling versus proton hop

Authors: O. Sipilä, P. Caselli, J. Harju

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

Submitted 10 September, 2019; originally announced September 2019.

Comments: Accepted to A&A; 14 pages, 13 figures

Journal ref: A&A 631, A63 (2019)

Gas and dust temperature in pre-stellar cores revisited: New limits on cosmic-ray ionization rate

Authors: Alexei V. Ivlev, Kedron Silsbee, Olli Sipilä, Paola Caselli

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

Submitted 9 September, 2019; originally announced September 2019.

Comments: Accepted to ApJ

The chemical structure of young high-mass star-forming clumps: (I) Deuteration

Authors: S. Feng, P. Caselli, K. Wang, Y. Lin, H. Beuther, O. Sipilä

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

Submitted 31 August, 2019; originally announced September 2019.

High sensitivity maps of molecular ions in L1544: I. Deuteration of N2H+ and HCO+ and first evidence of N2D+ depletion

Authors: Elena Redaelli, Luca Bizzocchi, Paola Caselli, Olli Sipilä, Valerio Lattanzi, Barbara Michela Giuliano, Silvia Spezzano

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

Submitted 18 July, 2019; originally announced July 2019.

Comments: Accepted for publication on A&A on 18/7/19

Journal ref: A&A 629, A15 (2019)

Deuterated forms of H${_3^+}$ and their importance in astrochemistry

Authors: Paola Caselli, Olli Sipilä, Jorma Harju

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

Submitted 21 May, 2019; originally announced May 2019.

Comments: accepted for publication in the Philosophical Transactions A of the Royal Society

Why does ammonia not freeze out in the center of pre-stellar cores?

Authors: O. Sipilä, P. Caselli, E. Redaelli, M. Juvela, L. Bizzocchi

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

Submitted 7 May, 2019; originally announced May 2019.

Comments: Accepted to MNRAS

Efficient methanol production on the dark side of a prestellar core

Authors: Jorma Harju, Jaime E. Pineda, Anton I. Vasyunin, Paola Caselli, Stella S. R. Offner, Alyssa A. Goodman, Mika Juvela, Olli Sipilae, Alexandre Faure, Romane Le Gal, Pierre Hily-Blant, Joao Alves, Luca Bizzocchi, Andreas Burkert, Hope Chen, Rachel K. Friesen, Rolf Guesten, Philip C. Myers, Anna Punanova, Claire Rist, Erik Rosolowsky, Stephan Schlemmer, Yancy Shirley, Silvia Spezzano, Charlotte Vastel , et al. (1 additional authors not shown)

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

Submitted 4 May, 2020; v1 submitted 27 March, 2019; originally announced March 2019.

Comments: Accepted for publication in the Astrophysical Journal

Mapping deuterated methanol toward L1544: I. Deuterium fraction and comparison with modeling

Authors: A. Chacon-Tanarro, P. Caselli, L. Bizzocchi, J. E. Pineda, O. Sipila, A. Vasyunin, S. Spezzano, A. Punanova, B. M. Giuliano, V. Lattanzi

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

Submitted 29 August, 2018; originally announced August 2018.

Comments: Accepted for publication in A&A

Journal ref: A&A 622, A141 (2019)

A study of the $c$-$\mathrm{C_{3}HD}$/$c$-$\mathrm{C_{3}H_{2}}$ ratio in low-mass star forming regions

Authors: J. Chantzos, S. Spezzano, P. Caselli, A. Chacón-Tanarro, L. Bizzocchi, O. Sipilä, B. M. Giuliano

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

Submitted 12 July, 2018; originally announced July 2018.

Hydrodynamics with gas-grain chemistry and radiative transfer: comparing dynamical and static models

Authors: O. Sipilä, P. Caselli

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

Submitted 24 January, 2018; originally announced January 2018.

Comments: Accepted by A&A; 15 pages, 12 figures, 5 tables; arXiv abstract heavily modified and redacted from the original

Journal ref: A&A 615, A15 (2018)

Species-to-species rate coefficients for the $\rm H_3^+ + H_2$ reacting system

Authors: O. Sipilä, J. Harju, P. Caselli

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

Submitted 11 July, 2017; originally announced July 2017.

Comments: accepted to A&A; 22 pages, 8 figures, 2 appendices; arXiv abstract heavily redacted and somewhat modified

Journal ref: A&A 607, A26 (2017)

NH_3(1_0-0_0) in the pre-stellar core L1544

Authors: P. Caselli, L. Bizzocchi, E. Keto, O. Sipila, M. Tafalla, L. Pagani, L. E. Kristensen, F. F. S. van der Tak, C. M. Walmsley, C. Codella, B. Nisini, Y. Aikawa, A. Faure, E. F. van Dishoeck

Abstract: Pre-stellar cores represent the initial conditions in the process of star and planet formation, therefore it is important to study their physical and chemical structure. Because of their volatility, nitrogen-bearing molecules are key to study the dense and cold gas present in pre-stellar cores. The NH_3 rotational transition detected with Herschel-HIFI provides a unique combination of sensitivity… ▽ More Pre-stellar cores represent the initial conditions in the process of star and planet formation, therefore it is important to study their physical and chemical structure. Because of their volatility, nitrogen-bearing molecules are key to study the dense and cold gas present in pre-stellar cores. The NH_3 rotational transition detected with Herschel-HIFI provides a unique combination of sensitivity and spectral resolution to further investigate physical and chemical processes in pre-stellar cores. Here we present the velocity-resolved Herschel-HIFI observations of the ortho-NH_3(1_0-0_0) line at 572 GHz and study the abundance profile of ammonia across the pre-stellar core L1544 to test current theories of its physical and chemical structure. Recently calculated collisional coefficients have been included in our non-LTE radiative transfer code to reproduce Herschel observations. A gas-grain chemical model, including spin-state chemistry and applied to the (static) physical structure of L1544 is also used to infer the abundance profile of ortho-NH_3 . The hyperfine structure of ortho-NH_3(1_0-0_0) is resolved for the first time in space. All the hyperfine components are strongly self-absorbed. The profile can be reproduced if the core is contracting in quasi-equilibrium, consistent with previous work, and if the NH_3 abundance is slightly rising toward the core centre, as deduced from previous interferometric observations of para-NH_3(1,1). The chemical model overestimates the NH_3 abundance at radii between ~ 4000 and 15000 AU by about two orders of magnitude and underestimates the abundance toward the core centre by more than one order of magnitude. Our observations show that chemical models applied to static clouds have problems in reproducing NH_3 observations. △ Less

Submitted 9 June, 2017; originally announced June 2017.

Comments: accepted for publication in A&A Letters

Detection of interstellar ortho-D2H+ with SOFIA

Authors: Jorma Harju, Olli Sipilä, Sandra Brünken, Stephan Schlemmer, Paola Caselli, Mika Juvela, Karl M. Menten, Jürgen Stutzki, Oskar Asvany, Tomasz Kaminski, Yoko Okada, Ronan Higgins

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

Submitted 8 April, 2017; originally announced April 2017.

Comments: 19 pages, 12 figures, accepted for publication in The Astrophysical Journal

Accurate sub-millimetre rest-frequencies for HOCO$^+$ and DOCO$^+$ ions

Authors: Luca Bizzocchi, Valerio Lattanzi, Jacob Laas, Silvia Spezzano, Barbara Michela Giuliano, Domenico Prudenzano, Christian Endres, Olli Sipilä, Paola Caselli

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

Submitted 24 March, 2017; originally announced March 2017.

Comments: 7 pages, 1 figure, accepted for publication in A&A

Journal ref: A&A 602, A34 (2017)

On the stability of nonisothermal Bonnor-Ebert spheres. III. The role of chemistry in core stabilization

Authors: O. Sipilä, P. Caselli, M. Juvela

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

Submitted 21 March, 2017; originally announced March 2017.

Comments: 8 pages, 10 figures; accepted for publication in A&A

Journal ref: A&A 601, A113 (2017)

Chemistry of TMC-1 with multiply deuterated species and spin chemistry of H2, H2+, H3+ and their isotopologues

Authors: L. Majumdar, P. Gratier, M. Ruaud, V. Wakelam, C. Vastel, O. Sipilä, F. Hersant, A. Dutrey, S. Guilloteau

Abstract: 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 ch… ▽ More 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. △ Less

Submitted 22 December, 2016; originally announced December 2016.

Comments: Accepted for publication in MNRAS on Dec. 22nd

Understanding the $\rm C_3H_2$ cyclic-to-linear ratio in L1544

Authors: O. Sipilä, S. Spezzano, P. Caselli

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

Submitted 24 May, 2016; originally announced May 2016.

Comments: 5 pages, 2 figures; accepted as an A&A Letter

Journal ref: A&A 591, L1 (2016)

Deuteration of ammonia in the starless core Ophiuchus/H-MM1

Authors: Jorma Harju, Fabien Daniel, Olli Sipilä, Paola Caselli, Jaime E. Pineda, Rachel K. Friesen, Anna Punanova, Rolf Güsten, Laurent Wiesenfeld, Philip C. Myers, Alexandre Faure, Pierre Hily-Blant, Claire Rist, Erik Rosolowsky, Stephan Schlemmer, Yancy L. Shirley

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

Submitted 22 December, 2016; v1 submitted 19 April, 2016; originally announced April 2016.

Comments: to appear in Astronomy & Astrophysics

Journal ref: A&A 600, A61 (2017)

The effect of multilayer ice chemistry on gas-phase deuteration in starless cores

Authors: O. Sipilä, P. Caselli, V. Taquet

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

Submitted 18 April, 2016; originally announced April 2016.

Comments: 13 pages, 11 figures; accepted by A&A; abstract redacted and modified to fit arXiv format

Journal ref: A&A 591, A9 (2016)

The NH$_2$D hyperfine structure revealed by astrophysical observations

Authors: F. Daniel, L. H. Coudert, A. Punanova, J. Harju, A. Faure, E. Roueff, O. Sipilä, P. Caselli, R. Güsten, A. Pon, J. E. Pineda

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

Submitted 2 January, 2016; originally announced January 2016.

Journal ref: A&A 586, L4 (2016)

On the stability of non-isothermal Bonnor-Ebert spheres. II. The effect of gas temperature on the stability

Authors: O. Sipilä, J. Harju, M. Juvela

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

Submitted 25 August, 2015; originally announced August 2015.

Comments: 7 pages, 5 figures; accepted for publication in A&A; abstract (heavily) abridged for arXiv