Search

Deuterium was primarily created during the Big Bang Nucleosynthesis (BBN). This fact, alongside its fractionation reactions resulting in enhanced abundances of deuterated molecules, means that these abundances can be used to better understand many processes within the interstellar medium (ISM), as well as its history. Previously, observations of deuterated molecules have been limited to the Galaxy, the Magellanic Clouds and (with respect to HD) to quasar absorption spectra. We present the first robust detection of a deuterated molecule in a starburst environment and, besides HD, the first one detected outside the Local Group. We therefore can constrain the deuterium fractionation, as observed by DCN. We observed the CMZ of the nearby starburst galaxy NGC 253 covering multiple Giant Molecular Clouds (GMC) with cloud scale observations ($\sim 30$ pc) using ALMA. Via the use of the \texttt{MADCUBA} package we were able to perform LTE analysis in order to obtain deuterium fractionation estimates. We detect DCN in the nuclear region of the starburst galaxy NGC 253 and estimate the deuterium fractionation (D/H ratio) of DCN within the GMCs of the CMZ of NGC 253. We find a range of $5 \times 10^{-4}$ to $10 \times 10^{-4}$, relatively similar to the values observed in warm Galactic star-forming regions. We also determine an upper limit of D/H of $8 \times 10^{-5}$ from DCO\plus within one region, closer to the cosmic value of D/H. Our observations of deuterated molecules within NGC 253 appear to be consistent with previous galactic studies of star forming regions. This implies that warmer gas temperatures increase the abundance of DCN relative to other deuterated species. This study also further expands the regions, particularly in the extragalactic domain, in which deuterated species have been detected., Comment: Accepted for Publication to Astronomy & Astrophysics. 15 pages, 11 figures (Including appendices)

How Supernova Remnant (SNR) shocks impact nearby molecular clouds is still poorly observationally constrained. It is unclear if SNRs can positively or negatively affect clouds star formation potential. We have studied the dense gas morphology and kinematics toward the Infrared Dark Cloud (IRDC) G034.77-00.55, shock-interacting with the SNR W44, to identify evidence of early stage star formation induced by the shock. We have used high-angular resolution N2H+(1-0) images across G034.77-00.55, obtained with ALMA. N2H+ is a well known tracer of dense and cold material, optimal to identify gas with the highest potential to harbour star formation. The N2H+ emission is distributed into two elongated structures, one toward the dense ridge at the edge of the source and one toward the inner cloud. Both elongations are spatially associated with well-defined mass-surface density features. The velocities of the gas in the two structures i.e., 38-41 km s-1 and 41-43 km s-1 are consistent with the lowest velocities of the J- and C-type parts of the SNR-driven shock, respectively. A third velocity component is present at 43-45.5 km s-1. The dense gas shows a fragmented morphology with core-like fragments of scales consistent with the Jeans lengths, masses $\sim$1-20 M$_{\odot}$, densities (n(H$_2$)$\geq$10$^5$ cm$^{-3}$) sufficient to host star formation in free-fall time scales (few 10$^4$ yr) and with virial parameters that hint toward possible collapse. The W44 driven shock may have swept up the encountered material which is now seen as a dense ridge, almost detached from the main cloud, and an elongation within the inner cloud, well constrained in both N2H+ emission and mass surface density. This shock compressed material may have then fragmented into cores that are either in a starless or pre-stellar stage. Additional observations are needed to confirm this scenario and the nature of the cores., Comment: 12 pages, 8 figures, accepted for publication in A&A

Single-dish observations suggest that the abundances of organic species in star-forming regions of the outer Galaxy, characterised by sub-Solar metallicities, are comparable to those found in the local Galaxy. To understand this counter-intuitive result, and avoid misleading interpretation due to beam dilution effects at such large distances, spatially resolved molecular emission maps are needed to link correctly measured abundances and local physical properties. We observed several organic molecules with the Atacama Large Millimeter Array towards WB89-671, the source with the largest Galactocentric distance (23.4~kpc) of the project "CHEMical complexity in star-forming regions of the OUTer Galaxy" (CHEMOUT), at a resolution of 15000~au. We compared the observed molecular abundances with chemical model predictions. We detected emission of c-C3H2, C4H, CH3OH, H2CO, HCO, H13CO+, HCS+, CS, HN13C, and SO. The emission morphology is complex, extended, and different in each tracer. The most intense emission in H13CO+, H2CO and c-C3H2 arises from two millimeter continuum, infrared-bright cores. The most intense CH3OH and SO emission arises predominantly from the part of the filament with no continuum sources. The narrow linewidths across the filament indicate quiescent gas, despite the two embedded protostars. Derived molecular column densities are comparable with those in local star-forming regions, and suggest anti-correlation between hydrocarbons, ions, HCO, and H2CO on one side, and CH3OH and SO on the other. Static chemical models that best match the observed column densities favour low energetic conditions, expected at large Galactocentric radii, but carbon elemental abundances 3 times higher than that derived extrapolating the [C/H] Galactocentric gradient at 23~kpc. This would indicate a flatter [C/H] trend at large Galactocentric radii, in line with a flat abundance of organics., Comment: accepted for publication in Astronomy & Astrophysics, 17 pages, 14 figures, 6 tables, 3 Appendices

We use data from the ALMA Evolutionary Study of High Mass Protocluster Formation in the Galaxy (ALMAGAL) survey to study 100 ALMAGAL regions at $\sim$ 1 arsecond resolution located between $\sim$ 2 and 6 kpc distance. Using ALMAGAL $\sim$ 1.3mm line and continuum data we estimate flow rates onto individual cores. We focus specifically on flow rates along filamentary structures associated with these cores. Our primary analysis is centered around position velocity cuts in H$_2$CO (3$_{0,3}$ - 2$_{0,2}$) which allow us to measure the velocity fields, surrounding these cores. Combining this work with column density estimates we derive the flow rates along the extended filamentary structures associated with cores in these regions. We select a sample of 100 ALMAGAL regions covering four evolutionary stages from quiescent to protostellar, Young Stellar Objects (YSOs), and HII regions (25 each). Using dendrogram and line analysis, we identify a final sample of 182 cores in 87 regions. In this paper, we present 728 flow rates for our sample (4 per core), analysed in the context of evolutionary stage, distance from the core, and core mass. On average, for the whole sample, we derive flow rates on the order of $\sim$10$^{-4}$ M$_{sun}$yr$^{-1}$ with estimated uncertainties of $\pm$50%. We see increasing differences in the values among evolutionary stages, most notably between the less evolved (quiescent/protostellar) and more evolved (YSO/HII region) sources. We also see an increasing trend as we move further away from the centre of these cores. We also find a clear relationship between the flow rates and core masses $\sim$M$^{2/3}$ which is in line with the result expected from the tidal-lobe accretion mechanism. Overall, we see increasing trends in the relationships between the flow rate and the three investigated parameters; evolutionary stage, distance from the core, and core mass., Comment: 11 pages, 11 figures, accepted for publication in A&A

Recent observations suggest that planets formation starts early, in protostellar disks of $\le10^5$ yrs, which are characterized by strong interactions with the environment, e.g., through accretion streamers and molecular outflows. To investigate the impact of such phenomena on disk physical and chemical properties it is key to understand what chemistry planets inherit from their natal environment. In the context of the ALMA Large Program Fifty AU STudy of the chemistry in the disk/envelope system of Solar-like protostars (FAUST), we present observations on scales from ~1500 au to ~60 au of H$_2$CO, HDCO, and D$_2$CO towards the young planet-forming disk IRS~63. H$_2$CO probes the gas in the disk as well as in a large scale streamer (~1500 au) impacting onto the South-East (SE) disk side. We detect for the first time deuterated formaldehyde, HDCO and D$_2$CO, in a planet-forming disk, and HDCO in the streamer that is feeding it. This allows us to estimate the deuterium fractionation of H$_2$CO in the disk: [HDCO]/[H$_2$CO]$\sim0.1-0.3$ and [D$_2$CO]/[H$_2$CO]$\sim0.1$. Interestingly, while HDCO follows the H$_2$CO distribution in the disk and in the streamer, the distribution of D$_2$CO is highly asymmetric, with a peak of the emission (and [D]/[H] ratio) in the SE disk side, where the streamer crashes onto the disk. In addition, D$_2$CO is detected in two spots along the blue- and red-shifted outflow. This suggests that: (i) in the disk, HDCO formation is dominated by gas-phase reactions similarly to H$_2$CO, while (ii) D$_2$CO was mainly formed on the grain mantles during the prestellar phase and/or in the disk itself, and is at present released in the gas-phase in the shocks driven by the streamer and the outflow. These findings testify on the key role of streamers in the build-up of the disk both concerning the final mass available for planet formation and its chemical composition., Comment: 12 pages, 10 figures, accepted for publication on A&A

The origin of the chemical diversity observed around low-mass protostars probably resides in the earliest history of these systems. We aim to investigate the impact of protostellar feedback on the chemistry and grain growth in the circumstellar medium of multiple stellar systems. In the context of the ALMA Large Program FAUST, we present high-resolution (50 au) observations of CH$_3$OH, H$_2$CO, and SiO and continuum emission at 1.3 mm and 3 mm towards the Corona Australis star cluster. Methanol emission reveals an arc-like structure at $\sim$1800 au from the protostellar system IRS7B along the direction perpendicular to the major axis of the disc. The arc is located at the edge of two elongated continuum structures that define a cone emerging from IRS7B. The region inside the cone is probed by H$_2$CO, while the eastern wall of the arc shows bright emission in SiO, a typical shock tracer. Taking into account the association with a previously detected radio jet imaged with JVLA at 6 cm, the molecular arc reveals for the first time a bow shock driven by IRS7B and a two-sided dust cavity opened by the mass-loss process. For each cavity wall, we derive an average H$_2$ column density of $\sim$7$\times$10$^{21}$ cm$^{-2}$, a mass of $\sim$9$\times$10$^{-3}$ M$_\odot$, and a lower limit on the dust spectral index of $1.4$. These observations provide the first evidence of a shock and a conical dust cavity opened by the jet driven by IRS7B, with important implications for the chemical enrichment and grain growth in the envelope of Solar System analogues., Comment: 12 pages, 8 figures, 3 tables. Accepted Letter in Astronomy & Astrophysics

iCOMs are species commonly found in the interstellar medium. They are believed to be crucial seed species for the build-up of chemical complexity in star forming regions as well as our own Solar System. Thus, understanding how their abundances evolve during the star formation process and whether it enriches the emerging planetary system is of paramount importance. We use data from the ALMA Large Program FAUST to study the compact line emission towards the [BHB2007] 11 proto-binary system (sources A and B), where a complex structure of filaments connecting the two sources with a larger circumbinary disk has previously been detected. More than 45 CH3OCHO lines are clearly detected, as well as 8 CH3OCH3 transitions , 1 H2CCO transition and 4 t-HCOOH transitions. We compute the abundance ratios with respect to CH3OH for CH3OCHO, CH3OCH3, H2CCO, t-HCOOH (as well as an upper limit for CH3CHO) through a radiative transfer analysis. We also report the upper limits on the column densities of nitrogen bearing iCOMs, N(C2H5CN) and N(C2H3CN). The emission from the detected iCOMs and their precursors is compact and encompasses both protostars, which are separated by only 0.2" (~ 28 au). The integrated intensities tend to align with the Southern filament, revealed by the high spatial resolution observations of the dust emission at 1.3 mm. A PV and 2D analysis are performed on the strongest and uncontaminated CH3OCH3 transition and show three different spatial and velocity regions, two of them being close to 11B (Southern filament) and the third one near 11A. All our observations suggest that the detected methanol, as well as the other iCOMs, are generated by the shocked gas from the incoming filaments streaming towards [BHB2007] 11A and 11B, respectively, making this source one of the few where chemical enrichment of the gas caused by the streaming material is observed., Comment: Accepted in A&A

Most stars, including the Sun, are born in rich stellar clusters containing massive stars. Therefore, the study of the chemical reservoir of massive star-forming regions is crucial to understand the basic chemical ingredients available at the dawn of planetary systems. We present a detailed study of the molecular inventory of the hot molecular core G31.41+0.31 from the project GUAPOS (G31.41+0.31 Unbiased ALMA sPectral Observational Survey). We analyze 34 species for the first time plus 20 species analyzed in previous GUAPOS works, including oxygen, nitrogen, sulfur, phosphorus, and chlorine species. We compare the abundances derived in G31.41+0.31 with those observed in other chemically-rich sources that represent the initial and last stages of the formation of stars and planets: the hot corino in the Solar-like protostar IRAS 16293-2422 B, and the comets 67P/Churyumov-Gerasimenko and 46P/Wirtanen. The comparative analysis reveals that the chemical feedstock of the two star-forming regions are similar. The abundances of oxygen- and nitrogen-bearing molecules exhibit a good correlation for all pair of sources, including the two comets, suggesting a chemical heritage of these species during the process of star formation, and hence an early phase formation of the molecules. However, sulfur- and phosphorus-bearing species present worse correlations, being more abundant in comets. This suggests that while sulfur- and phosphorus-bearing species are predominantly trapped on the surface of icy grains in the hot close surroundings of protostars, they could be more easily released into gas phase in comets, allowing their cosmic abundances to be almost recovered., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)

The astrochemistry of the important biogenic element phosphorus (P) is still poorly understood, but observational evidence indicates that P-bearing molecules are likely associated with shocks. We study P-bearing molecules, as well as some shock tracers, towards one of the chemically richest hot molecular core, G31.41+0.31, in the framework of the project "G31.41+0.31 Unbiased ALMA sPectral Observational Survey" (GUAPOS), observed with the Atacama Large Millimeter Array (ALMA). We have observed the molecules PN, PO, SO, SO2, SiO, and SiS, through their rotational lines in the spectral range 84.05-115.91 GHz, covered by the GUAPOS project. PN is clearly detected while PO is tentatively detected. The PN emission arises from two regions southwest of the hot core peak, "1" and "2", and is undetected or tentatively detected towards the hot core peak. the PN and SiO lines are very similar both in spatial emission morphology and spectral shape. Region "1" is in part overlapping with the hot core and it is warmer than region "2", which is well separated from the hot core and located along the outflows identified in previous studies. The column density ratio SiO/PN remains constant in regions "1" and "2", while SO/PN, SiS/PN, and SO2/PN decrease by about an order of magnitude from region "1" to region "2", indicating that SiO and PN have a common origin even in regions with different physical conditions. Our study firmly confirms previous observational evidence that PN emission is tightly associated with SiO and it is likely a product of shock-chemistry, as the lack of a clear detection of PN towards the hot-core allows to rule out relevant formation pathways in hot gas. We propose the PN emitting region "2" as a new astrophysical laboratory for shock-chemistry studies, Comment: 15 pages, 11 figures, accepted for publication in A&A

Context. Traditionally, supersonic turbulence is considered to be one of the most likely mechanisms to slow down the gravitational collapse in dense clumps, thereby enabling the formation of massive stars. However, several recent studies have raised differing points of view based on observations carried out with sufficiently high spatial and spectral resolution. These studies call for a re-evaluation of the role turbulence plays in massive star-forming regions. Aims. Our aim is to study the gas properties, especially the turbulence, in a sample of massive star-forming regions with sufficient spatial and spectral resolution, which can both resolve the core fragmentation and the thermal line width. Methods. We observed NH3 metastable lines with the Very Large Array (VLA) to assess the intrinsic turbulence. Results. Analysis of the turbulence distribution histogram for 32 identified NH3 cores reveals the presence of three distinct components. Furthermore, our results suggest that (1) sub- and transonic turbulence is a prevalent (21 of 32) feature of massive star-forming regions and those cold regions are at early evolutionary stage. This investigation indicates that turbulence alone is insufficient to provide the necessary internal pressure required for massive star formation, necessitating further exploration of alternative candidates; and (2) studies of seven multi-core systems indicate that the cores within each system mainly share similar gas properties and masses. However, two of the systems are characterized by the presence of exceptionally cold and dense cores that are situated at the spatial center of each system. Our findings support the hub-filament model as an explanation for this observed distribution, Comment: 34 pages, 15 figures, 4 tables. Accepted for publication on A&A

To understand the chemistry of sulphur (S) in the interstellar medium, models need to be tested by observations of S-bearing molecules in different physical conditions. We analyse observations obtained with the IRAM 30m telescope towards 15 well-known cores classified in the three main evolutionary stages of the high-mass star-formation process: high-mass starless cores,high-mass protostellar objects, and ultracompact HII regions. We detected rotational lines of SO, SO+, NS, C34S, 13CS, SO2, CCS, H2S, HCS+, OCS, H2CS, and CCCS. We also analyse for the first time lines of the NO molecule to complement the analysis. From a local thermodynamic equilibrium approach, we have derived column densities of each species and excitation temperatures. Based on a statistical analysis on the line widths and the excitation temperatures, we find that: NS, C34S, 13CS, CCS, and HCS+ trace cold, quiescent, and likely extended material; OCS, and SO2 trace warmer, more turbulent, and likely denser and more compact material; SO and perhaps SO+ trace both quiescent and turbulent material depending on the target. The abundances of SO, SO2, and H2S show the strongest positive correlations with the kinetic temperature, believed to be an evolutionary indicator. Moreover, the sum of all molecular abundances show an enhancement of gaseous S from the less evolved to the more evolved stages. These trends could be due to the increasing amount of S sputtered from dust grains owing to the increasing protostellar activity with evolution. The average abundances in each evolutionary group increase especially in the oxygen-bearing molecules, perhaps due to the increasing abundance of atomic oxygen with evolution owing to photodissociation of water in gas phase. Our observational work represents a test-bed for theoretical studies aimed at modelling the chemistry of sulphur during the evolution of high-mass star-forming cores, Comment: 16 pages. 3 appendices. Accepted for publication in A&A

The role of massive ($\geq$ 8M$_{\odot}$) stars in defining the energy budget and chemical enrichment of the interstellar medium in their host galaxy is significant. In this first paper from the Tracing Evolution in Massive Protostellar Objects (TEMPO) project we introduce a colour-luminosity selected (L$_*$ $\sim$ 3$\times10^3$ to 1$\times10^5$ L$_{\odot}$) sample of 38 massive star forming regions observed with ALMA at 1.3mm and explore the fragmentation, clustering and flux density properties of the sample. The TEMPO sample fields are each found to contain multiple fragments (between 2-15 per field). The flux density budget is split evenly (53%-47%) between fields where emission is dominated by a single high flux density fragment and those in which the combined flux density of fainter objects dominates. The fragmentation scales observed in most fields are not comparable with the thermal Jeans length, $\lambda_J$, being larger in the majority of cases, suggestive of some non-thermal mechanism. A tentative evolutionary trend is seen between luminosity of the clump and the `spectral line richness' of the TEMPO fields; with 6.7GHz maser associated fields found to be lower luminosity and more line rich. This work also describes a method of line-free continuum channel selection within ALMA data and a generalised approach used to distinguishing sources which are potentially star-forming from those which are not, utilising interferometric visibility properties., Comment: 42 pages, 17 figures, 4 Tables including 3 Appendices (with additional tables and figures). Accepted for publication in MNRAS

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., Comment: Accepted to A&A; abstract abridged to meet arXiv requirements

The G31.41+0.31 Unbiased ALMA sPectral Observational Survey (GUAPOS) project targets the hot molecular core (HMC) G31.41+0.31 (G31), to unveil the complex chemistry of one of the most chemically rich high-mass star-forming regions outside the Galactic Center (GC). In the third paper of the project, we present a study of nine O-bearing (CH$_3$OH, $^{13}$CH$_3$OH, CH$_3^{18}$OH, CH$_3$CHO, CH$_3$OCH$_3$, CH$_3$COCH$_3$ , C$_2$H$_5$OH, aGg'-(CH$_2$OH)$_2$, and gGg'-(CH$_2$OH)$_2$) and six N-bearing (CH$_3$CN, $^{13}$CH$_3$CN, CH$_3^{13}$CN, C$_2$H$_3$CN, C$_2$H$_5$CN, and C$_2$H$_5^{13}$CN) complex organic molecules toward G31. The aim of this work is to characterize the abundances in one of the most chemically-rich hot molecular cores outside the GC and to search for a possible chemical segregation between O-bearing and N-bearing species in G31, which hosts four compact sources as seen with higher angular resolution data. Observations were carried out with the interferometer ALMA and covered the entire Band 3 from 84 to 116 GHz ($\sim 32$ GHz bandwidth) with an angular resolution of $1.2''$ ($\sim4400\,\mathrm{au}$). The spectrum has been analyzed with the tool SLIM of MADCUBA to determine the physical parameters of the emitting gas. Moreover, we have analyzed the morphology of the emission of the molecular species. We have compared the abundances w.r.t methanol of COMs in G31 with other twenty-seven sources, including other hot molecular cores inside and outside the Galactic Center, hot corinos, shocked regions, envelopes around young stellar objects, and quiescent molecular clouds, and with chemical models. Different species peak at slightly different positions, and this, together with the different central velocities of the lines obtained from the spectral fitting, point to chemical differentiation of selected O-bearing species., Comment: accepted for publication in A&A, 39 pages

Supernova remnants (SNRs) may regulate star formation in galaxies. For example, SNR-driven shocks may form new molecular gas or compress pre-existing clouds and trigger the formation of new stars. To test this scenario, we measure the deuteration of $N_2H^+$, $D_{frac}^{N_2H^+}$, a well-studied tracer of pre-stellar cores, across the Infrared Dark Cloud (IRDC) G034.77-00.55, known to be experiencing a shock interaction with the SNR W44. We use N$_2$H$^+$ and N$_2$D$^+$ J=1-0 single pointing observations obtained with the 30m antenna at the Instituto de Radioastronomia Millimetrica to infer $D_{frac}^{N_2H^+}$ toward five positions across the cloud, namely a massive core, different regions across the shock front, a dense clump and ambient gas. We find $D_{frac}^{N_2H^+}$ in the range 0.03-0.1, several orders of magnitude larger than the cosmic D/H ratio ($\sim$10$^{-5}$). Across the shock front, $D_{frac}^{N_2H^+}$ is enhanced by more than a factor of 2 ($D_{frac}^{N_2H^+}\sim$0.05-0.07) with respect to the ambient gas ($\leq$0.03) and similar to that measured generally in pre-stellar cores. Indeed, in the massive core and dense clump regions of this IRDC we measure $D_{frac}^{N_2H^+}$}$\sim$0.1. We find enhanced deuteration of $N_2H^+$ across the region of the shock, at a level that is enhanced with respect to regions of unperturbed gas. It is possible that this has been induced by shock compression, which would then be indirect evidence that the shock is triggering conditions for future star formation. However, since unperturbed dense regions also show elevated levels of deuteration, further, higher-resolution studies are needed to better understand the structure and kinematics of the deuterated material in the shock region, e.g., if it still in relatively diffuse form or already organised in a population of low-mass pre-stellar cores., Comment: Accepted for publication on A&A; 8 pages, 5 figures

Context: Core accretion models of massive star formation require the existence of massive, starless cores within molecular clouds. Yet, only a small number of candidates for such truly massive, monolithic cores are currently known. Aims: Here we analyse a massive core in the well-studied infrared-dark cloud (IRDC) called the 'dragon cloud' (also known as G028.37+00.07 or 'Cloud C'). This core (C2c1) sits at the end of a chain of a roughly equally spaced actively star-forming cores near the centre of the IRDC. Methods: We present new high-angular resolution 1 mm ALMA dust continuum and molecular line observations of the massive core. Results: The high-angular resolution observations show that this region fragments into two cores C2c1a and C2c1b, which retain significant background-subtracted masses of 23 Msun and 2 Msun (31 Msun and 6 Msun without background subtraction), respectively. The cores do not appear to fragment further on the scales of our highest angular resolution images (0.200 arcsec, 0.005 pc ~ 1000 AU). We find that these cores are very dense (nH2 > 10^6 cm-3) and have only trans-sonic non-thermal motions (Ms ~ 1). Together the mass, density and internal motions imply a virial parameter of < 1, which suggests the cores are gravitationally unstable, unless supported by strong magnetic fields with strengths of ~ 1 - 10 mG. From CO line observations, we find that there is tentative evidence for a weak molecular outflow towards the lower-mass core, and yet the more massive core remains devoid of any star formation indicators. Conclusions: We present evidence for the existence of a massive, pre-stellar core, which has implications for theories of massive star formation. This source warrants follow-up higher-angular-resolution observations to further assess its monolithic and pre-stellar nature., Comment: 8+4 pages, 4+2 Figures, 2 Tables, accepted for publication in Astronomy & Astrophysics

The study of the early phases of star and planet formation is important to understand the physical and chemical history of stellar systems such as our own. In particular, protostars born in rich clusters are prototypes of the young Solar System. In the framework of the Seeds Of Life In Space (SOLIS) large observational project, the aim of the present work is to investigate the origin of the previously inferred high flux of energetic particles in the protocluster FIR4 of the Orion Molecular Cloud 2 (OMC-2), which appears asymmetric within the protocluster itself. Interferometric observations carried out with the IRAM NOEMA interferometer were used to map the silicon monoxide (SiO) emission around the FIR4 protocluster. Complementary archival data from the ALMA interferometer were also employed to help constrain excitation conditions. A physical-chemical model was implemented to characterise the particle acceleration along the protostellar jet candidate, along with a non-LTE analysis of the SiO emission along the jet. The emission morphology of the SiO rotational transitions hints for the first time at the presence of a collimated jet originating very close to the brightest protostar in the cluster, HOPS-108. The NOEMA observations unveiled a possible jet in the OMC-2 FIR4 protocluster propagating towards a previously measured enhanced cosmic-ray ionisation rate. This suggests that energetic particle acceleration by the jet shock close to the protostar might be at the origin of the enhanced cosmic-ray ionisation rate, as confirmed by modelling the protostellar jet., Comment: 8 pages, 7 figures, accepted for publications in A&A

Nitrogen isotopic ratios are a key tool for tracing Galactic stellar nucleosynthesis. We present the first study of the $^{14}$N/$^{15}$N abundance ratio in the outer regions of the Milky Way (namely, for galactocentric distances, $R_{\rm GC}$, from 12 kpc up to 19 kpc), with the aim to study the stellar nucleosynthesis effects in the global Galactic trend. We analysed IRAM 30m observations towards a sample of 35 sources in the context of the CHEMical complexity in star-forming regions of the OUTer Galaxy (CHEMOUT) project. We derived the $^{14}$N/$^{15}$N ratios from HCN and HNC for 14 and 3 sources, respectively, using the $J$ = 1-0 rotational transition of HN$^{13}$C, H$^{15}$NC, H$^{13}$CN, and HC$^{15}$N. The results found in the outer Galaxy have been combined with previous measurements obtained in the inner Galaxy. We find an overall linear decreasing H$^{13}$CN/HC$^{15}$N ratio with increasing $R_{\rm GC}$. This translates to a parabolic $^{14}$N/$^{15}$N ratio with a peak at 11 kpc. Updated Galactic chemical evolution models have been taken into account and compared with the observations. The parabolic trend of the $^{14}$N/$^{15}$N ratio with $R_{\rm GC}$ can be naturally explained (i) by a model that assumes novae as the main $^{15}$N producers on long timescales ($\ge$1 Gyr) and (ii) by updated stellar yields for low- and intermediate-mass stars., Comment: 19 pages, 8 figures, 6 tables, 4 appendix - Accepted for publication in A&A

Star-forming molecular clouds are characterised by the ubiquity of intertwined filaments. The filaments have been observed in both high- and low-mass star-forming regions, and are thought to split into collections of sonic fibres. The locations where filaments converge are termed hubs, and these are associated with the young stellar clusters. However, the observations of filamentary structures within hubs at distances require a high angular resolution that limits the number of such studies conducted so far. The integral shaped filament of the Orion A molecular cloud is noted for harbouring several hubs within which no filamentary structures have been observed so far. The goal of our study is to investigate the nature of the filamentary structures within one of these hubs, which is the chemically rich hub OMC-2 FIR 4, and to analyse their emission with high density and shock tracers. We observed the OMC-2 FIR 4 proto-cluster using Band 6 of the ALMA in Cycle 4 with an angular resolution of ~0.26"(100 au). We analysed the spatial distribution of dust, the shock tracer SiO, and dense gas tracers (i.e., CH$_{3}$OH, CS, and H$^{13}$CN). We also studied gas kinematics using SiO and CH3OH maps. Our observations for the first time reveal interwoven filamentary structures within OMC-2 FIR 4 that are probed by several tracers. Each filamentary structure is characterised by a distinct velocity as seen from the emission peak of CH$_{3}$OH lines. They also show transonic and supersonic motions. SiO is associated with filaments and also with multiple bow-shock features. In addition, for the first time, we reveal a highly collimated SiO jet (~1$^{\circ}$) with a projected length of ~5200 au from the embedded protostar VLA15. Our study shows that multi-scale observations of these regions are crucial for understanding the accretion processes and flow of material that shapes star formation.

Methanol is a ubiquitous species commonly found in the molecular interstellar medium. It is also a crucial seed species for the building-up of the chemical complexity in star forming regions. Thus, understanding how its abundance evolves during the star formation process and whether it enriches the emerging planetary system is of paramount importance. We used new data from the ALMA Large Program FAUST (Fifty AU STudy of the chemistry in the disk/envelope system of Solar-like protostars) to study the methanol line emission towards the [BHB2007] 11 protobinary system (sources A and B), where a complex structure of filaments connecting the two sources with a larger circumbinary disk has been previously detected. Twelve methanol lines have been detected with upper energies in the range [45-537] K along with one 13CH3OH transition. The methanol emission is compact and encompasses both protostars, separated by only 28 au and presents three velocity components, not spatially resolved by our observations, associated with three different spatial regions, with two of them close to 11B and the third one associated with 11A. A non-LTE radiative transfer analysis of the methanol lines concludes that the gas is hot and dense and highly enriched in methanol with an abundance as high as 1e-5. Using previous continuum data, we show that dust opacity can potentially completely absorb the methanol line emission from the two binary objects. Although we cannot firmly exclude other possibilities, we suggest that the detected hot methanol is resulting from the shocked gas from the incoming filaments streaming towards [BHB2007] 11 A and B, respectively. Higher spatial resolution observations are necessary to confirm this hypothesis., Comment: 17 pages, accepted in A&A

The outer Galaxy is an environment with metallicity lower than the Solar one and, because of this, the formation and survival of molecules in star-forming regions located in the inner and outer Galaxy is expected to be different. To gain understanding on how chemistry changes throughout the Milky Way, it is crucial to observe outer Galaxy star-forming regions to constrain models adapted for lower metallicity environments. The project "chemical complexity in star-forming regions of the outer Galaxy" (CHEMOUT) aims to address this problem observing a sample of 35 high-mass star-forming cores at Galactocentric distances up to ~23 kpc with the IRAM 30m telescope in various 3mm and 2mm bands. In this work we analyse observations of methanol (CH3OH), one of the simplest complex organic molecules crucial for organic chemistry in star-forming regions, and of two chemically related species, HCO and formaldehyde (H2CO), towards 15 out of the 35 targets of the CHEMOUT sample. In fact, only targets previously detected in both HCO and H2CO, both precursors of methanol, were considered. We detected CH3OH in all 15 targets. Using a Local Thermodynamic Equilibrium approach, we derive CH3OH excitation temperatures in the range 7 - 16 K and line widths smaller than 4 km/s, consistent with emission from a cold and quiescent envelope. The CH3OH fractional abundances w.r.t. H2 range between ~0.6 x 10^{-9} and ~7.4 x 10^{-9}. These values are comparable to those found in star-forming regions in the inner and local Galaxy. Our results have important implications in the organic, and possibly pre-biotic, chemistry occurring in the outermost star-forming regions of the Galaxy, and can help setting the frontiers of the Galactic habitable zone., Comment: 13 pages, 5 figures, 2 appendices, accepted for publication in A&A

Protostellar jets are an important agent of star formation feedback, tightly connected with the mass-accretion process. The history of jet formation and mass-ejection provides constraints on the mass accretion history and the nature of the driving source. We want to characterize the time-variability of the mass-ejection phenomena at work in the Class 0 protostellar phase, in order to better understand the dynamics of the outflowing gas and bring more constraints on the origin of the jet chemical composition and the mass-accretion history. We have observed the emission of the CO 2-1 and SO N_J=5_4-4_3 rotational transitions with NOEMA, towards the intermediate-mass Class 0 protostellar system Cep E. The CO high-velocity jet emission reveals a central component associated with high-velocity molecular knots, also detected in SO, surrounded by a collimated layer of entrained gas. The gas layer appears to accelerate along the main axis over a length scale delta_0 ~700 au, while its diameter gradually increases up to several 1000au at 2000au from the protostar. The jet is fragmented into 18 knots of mass ~10^-3 Msun, unevenly distributed between the northern and southern lobes, with velocity variations up to 15 km/s close to the protostar, well below the jet terminal velocities. The knot interval distribution is approximately bimodal with a scale of ~50-80yr close to the protostar and ~150-200yr at larger distances >12". The mass-loss rates derived from knot masses are overall steady, with values of 2.7x10^-5 Msun/yr (8.9x10^-6 Msun/yr) in the northern (southern) lobe. The interaction of the ambient protostellar material with high-velocity knots drives the formation of a molecular layer around the jet, which accounts for the higher mass-loss rate in the north. The jet dynamics are well accounted for by a simple precession model with a period of 2000yr and a mass-ejection period of 55yr., Comment: 13 pages, 9 figures, 3 table. Accepted in A&A

The outer Galaxy is an environment with metallicity lower than the Solar one. Because of this, the formation and survival of molecules in star-forming regions located in the inner and outer Galaxy is expected to be different. To gain understanding on how chemistry changes throughout the Milky Way, it is crucial to observe outer Galaxy star-forming regions to constrain models adapted for lower metallicity environments. In this paper we present a new observational project: chemical complexity in star-forming regions of the outer Galaxy (CHEMOUT). The goal is to unveil the chemical composition in 35 dense molecular clouds associated with star-forming regions of the outer Galaxy through observations obtained with the IRAM 30m telescope. In this first paper, we present the sample, and report the detection at 3~mm of simple organic species HCO+, H13CO+, HCN, c-C3H2, HCO, C4H, and HCS+, of the complex hydrocarbon CH3CCH, and of SiO, CCS and SO. From c-C3H2, we estimate new kinematic heliocentric and Galactocentric distances based on an updated rotation curve of the Galaxy. The detection of the molecular tracers does not seem to have a clear dependence on the Galactocentric distance. We also analyse the HCO+ line profiles. We find high velocity wings in ~71% of the targets, and their occurrence does not depend on the Galactocentric distance. Our results, confirmed by a statistical analysis, show that the presence of organic molecules and tracers of protostellar activity is ubiquitous in the low-metallicity environment of the outer Galaxy. Based on this, and on the additional evidence that small, terrestrial planets are omnipresent in the Galaxy, we support previous claims that the definition of Galactic Habitable Zone should be rediscussed in view of the ubiquitous capacity of the interstellar medium to form organic molecules., Comment: 11 pages, 5 figures, 3 appendices, accepted for publication in A&A

Supernova remnants (SNRs) contribute to regulate the star formation efficiency and evolution of galaxies. As they expand into the interstellar medium (ISM), they transfer vast amounts of energy and momentum that displace, compress and heat the surrounding material. Despite the extensive work in galaxy evolution models, it remains to be observationally validated to what extent the molecular ISM is affected by the interaction with SNRs. We use the first results of the ESO-ARO Public Spectroscopic Survey SHREC, to investigate the shock interaction between the SNR IC443 and the nearby molecular clump G. We use high sensitivity SiO(2-1) and H$^{13}$CO$^+$(1-0) maps obtained by SHREC together with SiO(1-0) observations obtained with the 40m telescope at the Yebes Observatory. We find that the bulk of the SiO emission is arising from the ongoing shock interaction between IC443 and clump G. The shocked gas shows a well ordered kinematic structure, with velocities blue-shifted with respect to the central velocity of the SNR, similar to what observed toward other SNR-cloud interaction sites. The shock compression enhances the molecular gas density, n(H$_2$), up to $>$10$^5$ cm$^{-3}$, a factor of >10 higher than the ambient gas density and similar to values required to ignite star formation. Finally, we estimate that up to 50\% of the momentum injected by IC443 is transferred to the interacting molecular material. Therefore the molecular ISM may represent an important momentum carrier in sites of SNR-cloud interactions., Comment: Accepted for Publication in MNRAS; 11 pages, 7 figures

New-generation spectroscopic surveys of the Milky Way plane have been revealing the structure of the interstellar medium, allowing the simultaneous study of dense structures from single star-forming objects or systems to entire spiral arms. We present the catalogue of molecular clouds extracted from the $^{13}$CO(1-0) data cubes of the Forgotten Quadrant Survey, which mapped the Galactic plane in the range 220\deg

Aims: The goal is to obtain a census of S-bearing species using interferometric images, towards SVS13-A, a Class I object associated with a hot corino rich in interstellar complex organic molecules. Methods: We used data at 3mm and 1.4mm obtained with IRAM-NOEMA in the framework of the Large Program SOLIS. Results: We imaged the spatial distribution of the line emission of 32SO, 34SO, C32}S, C34S, C33S, OCS, H2C32S, H2C34S, and NS. The low excitation (9 K) 32SO line is tracing the fast collimated jet driven by the nearby SVS13-B. Conversely, the rest of the lines are confined in the inner SVS13-A region, where complex organics have been previously imaged. The non-LTE LVG analysis of SO, SO2, and H2CS indicates a hot corino origin (60-120 au). Temperatures between 50 K and 300 K, and volume densities larger than 10^5 cm-3 have been derived. The abundances are in the following ranges: 0.3-6 10^-6 (CS), 7 10^-9} - 1 10^-7 (SO), 1-10 10^-7 (SO2), a few 10^-10 (H2CS and OCS), and 10^{-10} - 10^{-9}(NS). The N(NS)/N(NS^+) ratio is larger than 10, supporting that the NS^+ ion is mainly formed in the extended envelope. Conclusions: The [H2CS]/[H2CO] ratio increases with time (from Class 0 to Class II objects) by more than one order of magnitude. This suggests that [S]/[O] changes along the Sun-like star forming process. The estimate of the [S]/[H] budget in SVS13-A is 2%-17% of the Solar System value (1.8 10^-5), being consistent with what was previously measured towards Class 0 objects (1%-8%). This supports that the enrichment of the sulphuretted species with respect to dark clouds keeps constant from the Class 0 to the Class I stages of low-mass star formation. The present findings stress the importance of investigating the chemistry of star forming regions using large observational surveys as well as sampling regions on a Solar System scale., Comment: A&A, in press

Peptide-like bond molecules, which can take part to the formation of proteins in a primitive Earth environment, have been detected up to now only towards a few sources. We present a study of HNCO, HC(O)NH$_{2}$, CH$_{3}$NCO, CH$_{3}$C(O)NH$_{2}$, CH$_{3}$NHCHO, CH$_{3}$CH$_{2}$NCO, NH$_{2}$C(O)NH$_{2}$, NH$_{2}$C(O)CN, and HOCH$_{2}$C(O)NH$_{2}$ towards the hot core G31.41+0.31. We have used the spectrum obtained from the ALMA 3mm spectral survey GUAPOS, with an angular resolution of 1.2"$\times$1.2" ($\sim$4500 au), to derive column densities of all the molecular species, together with other 0.2"$\times$0.2" ($\sim$750 au) ALMA observations to study the morphology of HNCO, HC(O)NH$_{2}$ and CH$_{3}$C(O)NH$_{2}$. We have detected HNCO, HC(O)NH$_{2}$, CH$_{3}$NCO, CH$_{3}$C(O)NH$_{2}$, and CH$_{3}$NHCHO, for the first time all together outside the Galactic center. We have obtained molecular fractional abundances with respect to H$_{2}$ from 10$^{-7}$ down to a few 10$^{-9}$ and with respect to CH$_{3}$OH from 10$^{-3}$ to $\sim$4$\times$10$^{-2}$. From the comparison with other sources, we find that regions in an earlier stage of evolution, such as pre-stellar cores, show abundances at least two orders of magnitude lower than those in hot cores, hot corinos or shocked regions. Moreover, molecular abundance ratios towards different sources are found to be consistent between them within $\sim$1 order of magnitude, regardless of the physical properties (e.g. different masses and luminosities), or the source position throughout the Galaxy. New correlations between pairs of molecular abundances have also been found. These results suggest that all these species are formed on grain surfaces in early evolutionary stages of molecular clouds, and that they are subsequently released back to the gas-phase through thermal desorption or shock-triggered desorption., Comment: 51 pages, 24 figures, 5 tables, 5 appendix - Accepted for publication in A&A

The TOPG\"ot project studies a sample of 86 high-mass star-forming regions in different evolutionary stages from starless cores to ultra compact HII regions. The aim of the survey is to analyze different molecular species in a statistically significant sample to study the chemical evolution in high-mass star-forming regions, and identify chemical tracers of the different phases. The sources have been observed with the IRAM 30m telescope in different spectral windows at 1, 2, and 3 mm. In this first paper, we present the sample and analyze the spectral energy distributions (SEDs) of the TOPG\"ot sources to derive physical parameters. We use the MADCUBA software to analyze the emission of methyl cyanide (CH$_3$CN), a well-known tracer of high-mass star formation. The emission of the $\rm{CH_3CN(5_{K}-4_{K})}$ K-transitions has been detected towards 73 sources (85% of the sample), with 12 non-detections and one source not observed in the frequency range of $\rm{CH_3CN(5_{K}-4_{K})}$. The emission of CH$_3$CN has been detected towards all evolutionary stages, with the mean abundances showing a clear increase of an order of magnitude from high-mass starless-cores to later evolutionary stages. We found a conservative abundance upper limit for high-mass starless cores of $X_{\rm CH_3CN}<4.0\times10^{-11}$, and a range in abundance of $4.0\times10^{-11}

Deuterated molecules are good tracers of the evolutionary stage of star-forming cores. During the star formation process, deuterated molecules are expected to be enhanced in cold, dense pre-stellar cores and to deplete after protostellar birth. In this paper we study the deuteration fraction of formaldehyde in high-mass star-forming cores at different evolutionary stages to investigate whether the deuteration fraction of formaldehyde can be used as an evolutionary tracer. Using the APEX SEPIA Band 5 receiver, we extended our pilot study of the $J$=3$\rightarrow$2 rotational lines of HDCO and D$_2$CO to eleven high-mass star-forming regions that host objects at different evolutionary stages. High-resolution follow-up observations of eight objects in ALMA Band 6 were performed to reveal the size of the H$_2$CO emission and to give an estimate of the deuteration fractions HDCO/H$_2$CO and D$_2$CO/HDCO at scales of $\sim$6" (0.04-0.15 pc at the distance of our targets). Our observations show that singly- and doubly deuterated H$_2$CO are detected toward high-mass protostellar objects (HMPOs) and ultracompact HII regions (UCHII regions), the deuteration fraction of H$_2$CO is also found to decrease by an order of magnitude from the earlier HMPO phases to the latest evolutionary stage (UCHII), from $\sim$0.13 to $\sim$0.01. We have not detected HDCO and D$_2$CO emission from the youngest sources (high-mass starless cores, HMSCs). Our extended study supports the results of the previous pilot study: the deuteration fraction of formaldehyde decreases with evolutionary stage, but higher sensitivity observations are needed to provide more stringent constraints on the D/H ratio during the HMSC phase. The calculated upper limits for the HMSC sources are high, so the trend between HMSC and HMPO phases cannot be constrained., Comment: 15 pages, 4 figures, 4 tables, accepted for publication in A&A

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., Comment: 13 pages, 3 tables, 12 figures (+3 figures in appendix). Accepted for publication in Astronomy & Astrophysics

Infrared dark clouds (IRDCs) are potential hosts of the elusive early phases of high-mass star formation (HMSF). Here we conduct an in-depth analysis of the fragmentation properties of a sample of 10 IRDCs, which have been highlighted as some of the best candidates to study HMSF within the Milky Way. To do so, we have obtained a set of large mosaics covering these IRDCs with ALMA at band 3 (or 3mm). These observations have a high angular resolution (~3arcsec or ~0.05pc), and high continuum and spectral line sensitivity (~0.15mJy/beam and ~0.2K per 0.1km/s channel at the N2H+(1-0) transition). From the dust continuum emission, we identify 96 cores ranging from low- to high-mass (M = 3.4 to 50.9Msun) that are gravitationally bound (alpha_vir = 0.3 to 1.3) and which would require magnetic field strengths of B = 0.3 to 1.0mG to be in virial equilibrium. We combine these results with a homogenised catalogue of literature cores to recover the hierarchical structure within these clouds over four orders of magnitude in spatial scale (0.01pc to 10pc). Using supplementary observations at an even higher angular resolution, we find that the smallest fragments (<0.02pc) within this hierarchy do not currently have the mass and/or the density required to form high-mass stars. Nonetheless, the new ALMA observations presented in this paper have facilitated the identification of 19 (6 quiescent and 13 star-forming) cores that retain >16Msun without further fragmentation. These high-mass cores contain trans-sonic non-thermal motions, are kinematically sub-virial, and require moderate magnetic field strengths for support against collapse. The identification of these potential sites of high-mass star formation represents a key step in allowing us to test the predictions from high-mass star and cluster formation theories., Comment: 21(+6 in appendix) pages, 12(+3) figures, 3(+3) tables. Machine-readable versions of Table A1, A2 and A3 are available online. Accepted for publication in MNRAS

The 14N/15N ratio in molecules exhibits a large variation in star-forming regions, especially when measured from N2H+ isotopologues. However, there are only a few studies performed at high-angular resolution. We present the first interferometric survey of the 14N/15N ratio in N2H+ obtained with the Atacama Large Millimeter Array towards four infrared-dark clouds harbouring 3~mm continuum cores associated with different physical properties. We detect N15NH+ (1-0) in about 20-40% of the cores, depending on the host cloud. The 14N/15N values measured towards the millimeter continuum cores range from a minimum of 80 up to a maximum of 400. The spread of values is narrower than that found in any previous single-dish survey of high-mass star-forming regions, and than that obtained using the total power data only. This suggests that the 14N/15N ratio is on average higher in the diffuse gaseous envelope of the cores, and stresses the need for high-angular resolution maps to measure correctly the 14N/15N ratio in dense cores embedded in IRDCs. The average 14N/15N ratio of 210 is also lower than the interstellar value at the Galactocentric distance of the clouds (300-330), although the sensitivity of our observations does not allow us to unveil 14N/15N ratios higher than 400. No clear trend is found between the 14N/15N ratio and the core physical properties. We find only a tentative positive trend between 14N/15N and the H2 column density. However, firmer conclusions can be drawn only with higher sensitivity measurements., Comment: 19 pages, 7 figures, 6 tables, 1 appendix, accepted for publication in MNRAS

The isotopic ratio of nitrogen presents a wide range of values in the Solar System and in star forming system whose origin is still unclear. Chemical reactions in the gas phase are one of the possible processes that could modify the $^{14}$N/$^{15}$N ratio. We aim at investigating if and how the passage of a shock wave in the interstellar medium, can affect the relative fraction of nitrogen isotopes. The ideal place for such a study is the L1157 outflow, where several shocked clumps are present. We present the first measurement of the $^{14}$N/$^{15}$N ratio in the two shocked clumps, B1 and B0, of the protostellar outflow L1157, derived from the interferomteric maps of the H$^{13}$CN(1-0) and the HC$^{15}$N(1-0) lines. In B1, we find that the H$^{13}$CN(1-0) and HC$^{15}$N(1-0) emission traces the front of the clump, with averaged column density of $N$(H$^{13}$CN) $\sim$ 7$\times$10$^{12}$ cm$^{-2}$ and $N$(HC$^{15}$N) $\sim$ 2$\times$10$^{12}$ cm$^{-2}$. In this region the ratio H$^{13}$CN(1-0)/HC$^{15}$N(1-0) is quite uniform with an average value of $\sim$ 5$\pm$1. The same average value is also measured in the smaller clump B0e. Assuming the standard $^{12}$C/$^{13}$C = 68, we obtain $^{14}$N/$^{15}$N = 340$\pm$70, similar to those usually found in prestellar cores and protostars. We analysed the prediction of a chemical shock model for several shock conditions and we found that the nitrogen and carbon fractionations do not vary much for the first period after the shock. The observed H$^{13}$CN/HC$^{15}$N can be reproduced by a non-dissociative, C-type shock with parameters in agreement with previous modelling of L1157-B1. Both observations and chemical models indicate that the shock propagation does not affect the nitrogen isotopic ratio that remains similar to that measured in lower temperature gas in prestellar cores and in protostellar envelopes., Comment: 7 pages, 6 figures

Infrared Dark Clouds (IRDCs) are very dense and highly extincted regions that host the initial conditions of star and stellar cluster formation. It is crucial to study the kinematics and molecular content of IRDCs to test their formation mechanism and ultimately characterise these initial conditions. We have obtained high-sensitivity Silicon Monoxide, SiO(2-1), emission maps toward the six IRDCs, G018.82$-$00.28, G019.27+00.07, G028.53$-$00.25, G028.67+00.13, G038.95$-$00.47 and G053.11+00.05 (cloud A, B, D, E, I and J, respectively), using the 30-m antenna at the Instituto de Radioastronom\'{i}a Millim\'{e}trica (IRAM30m). We have investigated the SiO spatial distribution and kinematic structure across the six clouds to look for signatures of cloud-cloud collision events that may have formed the IRDCs and triggered star formation within them. Toward clouds A, B, D, I and J we detect spatially compact SiO emission with broad line profiles which are spatially coincident with massive cores. Toward the IRDCs A and I, we report an additional SiO component that shows narrow line profiles and that is widespread across quiescent regions. Finally, we do not detect any significant SiO emission toward cloud E. We suggest that the broad and compact SiO emission detected toward the clouds is likely associated with ongoing star formation activity within the IRDCs. However, the additional narrow and widespread SiO emission detected toward cloud A and I may have originated from the collision between the IRDCs and flows of molecular gas pushed toward the clouds by nearby HII regions., Comment: 17 pages, 9 figures

Understanding the degree of chemical complexity that can be reached in star-forming regions, together with the identification of precursors of the building blocks of life in the interstellar medium, is one of the goals of astrochemistry. Unbiased spectral surveys with large bandwidth and high spectral resolution are thus needed, to resolve line blending in chemically rich sources and identify complex organic molecules. This kind of observations has been successfully carried out, mainly towards the Galactic Center, a region that shows peculiar environmental conditions. We present an unbiased spectral survey at 3mm of one of the most chemically rich hot molecular cores located outside the Galactic Center, in the high-mass star-forming region G31.41+0.31. In this first paper, we present the survey and discuss the detection of the 3 isomers of C$_{2}$H$_{4}$O$_{2}$: methyl formate, glycolaldehyde and acetic acid. Observations were carried out with ALMA and cover the entire Band 3 from 84 to 116 GHz with an angular resolution of $1.2^{''}$x$1.2^{''}$ and a spectral resolution of $\sim0.488$ MHz. The transitions of the 3 molecules have been analyzed with the software XCLASS. All three isomers were detected and methyl formate and acetic acid abundances in G31 are the highest detected up to now, if compared to sources in literature. The size of the emission varies among the three isomers with acetic acid showing the most compact emission while methyl formate the most extended. The comparison with chemical models suggests the necessity of grain-surface routes for the formation of methyl formate in G31, while for glycolaldehyde both scenarios could be feasible. Proposed grain-surface reaction for acetic acid is not able to reproduce the observed abundance in this work, while gas-phase scenario should be further tested due to large uncertainties., Comment: 49 pages, 29 figures, accepted for publication in A&A

The study of hot corinos in Solar-like protostars has been so far mostly limited to the Class 0 phase, hampering our understanding of their origin and evolution. In addition, recent evidence suggests that planet formation starts already during Class I phase, which, therefore, represents a crucial step in the future planetary system chemical composition. Hence, the study of hot corinos in Class I protostars has become of paramount importance. Here we report the discovery of a hot corino towards the prototypical Class I protostar L1551 IRS5, obtained within the ALMA Large Program FAUST. We detected several lines from methanol and its isopotologues ($^{13}$CH$_{\rm 3}$OH and CH$_{\rm 2}$DOH), methyl formate and ethanol. Lines are bright toward the north component of the IRS5 binary system, and a possible second hot corino may be associated with the south component. The methanol lines non-LTE analysis constrains the gas temperature ($\sim$100 K), density ($\geq$1.5$\times$10$^{8}$ cm$^{-3}$), and emitting size ($\sim$10 au in radius). All CH$_{\rm 3}$OH and $^{13}$CH$_{\rm 3}$OH lines are optically thick, preventing a reliable measure of the deuteration. The methyl formate and ethanol relative abundances are compatible with those measured in Class 0 hot corinos. Thus, based on the present work, little chemical evolution from Class 0 to I hot corinos occurs., Comment: 6 pages, 2 figures

We have studied four complex organic molecules (COMs), methyl formate ($CH_3OCHO$), dimethyl ether ($CH_3OCH_3$), formamide ($NH_2CHO$), and ethyl cyanide ($C_2H_5CN$), towards a large sample of 39 high-mass star-forming regions representing different evolutionary stages, from early to evolved phases. We aim to identify potential correlations between the molecules and to trace their evolutionary sequence through the star formation process. We analysed spectra obtained at 3, 2, and 0.9 mm with the IRAM-30m telescope. We derived the main physical parameters for each species by fitting the molecular lines. We compared them and evaluated their evolution, also taking several other interstellar environments into account. We report detections in 20 sources, revealing a clear dust absorption effect on column densities. Derived abundances are ~$10^{-10}-10^{-7}$ for $CH_3OCHO$ and $CH_3OCH_3$, ~$10^{-12}-10^{-10}$ for $NH_2CHO$, and ~$10^{-11}-10^{-9}$ for $C_2H_5CN$. The abundances of $CH_3OCHO$, $CH_3OCH_3$, and $C_2H_5CN$ are very strongly correlated (r>0.92) across ~4 orders of magnitude. $CH_3OCHO$ and $CH_3OCH_3$ show the strongest correlations in most parameters, and a nearly constant ratio (~1) over a remarkable ~9 orders of magnitude in luminosity for a wide variety of sources: pre-stellar to evolved cores, low- to high-mass objects, shocks, Galactic clouds, and comets. This indicates that COMs chemistry is likely early developed and then preserved through evolved phases. Moreover, the molecular abundances clearly increase with evolution. We consider $CH_3OCHO$ and $CH_3OCH_3$ to be most likely chemically linked: they could e.g. share a common precursor, or be formed one from the other. We propose a general scenario for all COMs, involving a formation in the cold, earliest phases of star formation and a following increasing desorption with the progressive heating of the evolving core., Comment: 27 pages, 14 figures, 20 tables - Accepted for publication in A&A on 25 June 2020

Aims: A unique environment to study how interstellar Complex Organic Molecules (iCOMs) can be formed is the shocked gas along low-mass protostellar outflows, as the dust mantles composition is sputtered into the gas phase. The chemical richness in these environments has been so far studied only in the L1157 blue shifted outflow. Methods: To understand if the L1157-B1 case is unique, we imaged the NGC 1333 IRAS 4A outflows using the NOEMA (NOrthern Extended Millimeter Array) interferometer as part of the IRAM SOLIS (Seeds Of Life in Space) Large Program and compared the observations with the GRAINOBLE+ gas phase astrochemical model. Results: Several iCOMs were detected in the IRAS 4A outflows: methanol (CH$_3$OH), acetaldehyde (CH$_3$CHO), formamide (NH$_2$CHO) and dimethyl ether (CH$_3$OCH$_3$), all sampling upper excitation energy up to $\sim$30 K. We found a significant chemical differentiation between the IRAS 4A1 outflow, showing a richer molecular content, and the IRAS 4A2 one. The CH$_3$OH/CH$_3$CHO abundance ratio is lower by a factor $\sim$4 in the former; furthermore the ratio in both outflows is lower by a factor $\sim$10 with respect to hot corinos values. Conclusions: After L1157-B1, IRAS 4A outflow is now the second outflow to show an evident chemical complexity. Given that CH$_3$OH is a grain surface species, GRAINOBLE+ reproduced our observations assuming acetaldehyde formation in gas phase by the reaction of ethyl radical (CH$_3$CH$_2$) with atomic oxygen. Moreover, the chemical differentiation between the two outflows suggests that the IRAS 4A1 outflow is likely younger than the IRAS 4A2 one. Further investigation is needed to constrain the age of the outflow and observations of even younger shocks are necessary and future spectroscopic studies on CH$_3$CH$_2$ are needed to be able to observe this species and provide strong constraints on the CH$_3$CHO formation., Comment: 12 pages, 7 figures, 5 tables

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., Comment: The published exponent in the A&A journal of the second rate coefficient of Table 1 is a misprint

We present the study of deuteration of cyanoacetylene (HC$_3$N) towards a sample of 28 high-mass star-forming cores divided into different evolutionary stages, from starless to evolved protostellar cores. We report for the first time the detection of DC$_3$N towards 15 high-mass cores. The abundance ratios of DC$_3$N with respect HC$_3$N range in the interval 0.003$-$0.022, lower than those found in low-mas protostars and dark clouds. No significant trend with the evolutionary stage, or with the kinetic temperature of the region, has been found. We compare the level of deuteration of HC$_3$N with those of other molecules towards the same sample, finding weak correlation with species formed only or predominantly in gas phase (N$_2$H$^+$ and HNC, respectively), and no correlation with species formed only or predominantly on dust grains (CH$_3$OH and NH$_3$, respectively). We also present a single-dish map of DC$_3$N towards the protocluster IRAS 05358+3543, which shows that DC$_3$N traces an extended envelope ($\sim$0.37 pc) and peaks towards two cold condensations separated from the positions of the protostars and the dust continuum. The observations presented in this work suggest that deuteration of HC$_3$N is produced in the gas of the cold outer parts of massive star-forming clumps, giving us an estimate of the deuteration factor prior to the formation of denser gas., Comment: Accepted in Monthly Notices of the Royal Astronomical Society -- 11 pages, 7 Figures, 2 Tables. Version with some typos corrected

We present observations of SO and $\rm SO_2$ lines toward the shocked regions along the L1157 chemically rich outflow, taken in the context of the Seeds Of Life In Space IRAM-NOrthern Extended Millimeter Array Large Program, and supported by data from Submillimeter Array and IRAM-30 m telescope at 1.1--3.6 mm wavelengths. We simultaneously analyze, for the first time, all of the brightest shocks in the blueshifted lobe, namely, B0, B1, and B2. We found the following. (1) SO and $\rm SO_2$ may trace different gas, given that the large(-scale) velocity gradient analysis indicates for $\rm SO_2$ a volume density ($\rm 10^5\text{--}10^6\,cm^{-3}$) denser than that of the gas emitting in SO by a factor up to an order of magnitude. (2) Investigating the 0.1 pc scale field of view, we note a tentative gradient along the path of the precessing jet. More specifically, $\rm \chi({SO/SO_2})$ decreases from the B0-B1 shocks to the older B2. (3) At a linear resolution of 500--1400 au, a tentative spatial displacement between the two emitting molecules is detected, with the SO peak closer (with respect to $\rm SO_2$) to the position where the recent jet is impinging on the B1 cavity wall. Our astrochemical modeling shows that the SO and $\rm SO_2$ abundances evolve on timescales less than about 1000 years. Furthermore, the modeling requires high abundances ($2\times10^{-6}$) of both $\rm H_2S/H$ and S/H injected in the gas phase due to the shock occurrence, so pre-frozen OCS only is not enough to reproduce our new observations., Comment: 13 pages, 9 figures, accepted by ApJ on 2020 Apr. 7

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

We show the first interferometric maps of the 14N/15N ratio obtained with the Atacama Large Millimeter Array (ALMA) towards the Solar-like forming protocluster OMC-2 FIR4. We observed N2H+, 15NNH+, N15NH+ (1-0), and N2D+ (2-1), from which we derive the isotopic ratios 14N/15N and D/H. The target is one of the closest analogues of the environment in which our Sun may have formed. The ALMA images, having synthesised beam corresponding to ~600 au, show that the emission of the less abundant isotopologues is distributed in several cores of ~10" (i.e. ~0.02~pc or 4000 au) embedded in a more extended N2H+ emission. Overall, our results indicate that: (1) 14N/15N does not change across the region at core scales, and (2) 14N/15N does not depend on temperature variations. Our findings also suggest that the 14N/15N variations found in pristine Solar System objects are likely not inherited from the protocluster stage, and hence their reason has to be found elsewhere., Comment: 10 pages, 5 figures, accepted for publication by MNRAS

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

We present the Forgotten Quadrant Survey (FQS), an ESO large project that used the 12m antenna of the Arizona Radio Observatory to map the Galactic Plane in the range 220\deg$

Observations of young stellar objects (YSOs) in centimeter bands can probe the continuum emission from growing dust grains, ionized winds, and magnetospheric activity, which are intimately connected to the evolution of protoplanetary disks and the formation of planets. We have carried out sensitive continuum observations toward the Ophiuchus A star-forming region using the Karl G. Jansky Very Large Array (VLA) at 10 GHz over a field-of-view of 6$'$ with a spatial resolution of $\theta_{maj}$ $\times$ $\theta_{min}$ $\sim$ 0.4$''$ $\times$ 0.2$''$. We achieved a 5 $\mu$Jy beam$^{-1}$ root-mean-square noise level at the center of our mosaic field of view. Among the eighteen sources we detected, sixteen are YSOs (three Class 0, five Class I, six Class II, and two Class III) and two are extragalactic candidates. We find that thermal dust emission generally contributes less that 30% of the emission at 10 GHz. The radio emission is dominated by other types of emission such as gyro-synchrotron radiation from active magnetospheres, free-free emission from thermal jets, free-free emission from the outflowing photo-evaporated disk material, and/or synchrotron emission from accelerated cosmic-rays in jet or protostellar surface shocks. These different types of emission could not be clearly disentangled. Our non-detections towards Class II/III disks suggest that extreme UV-driven photoevaporation is insufficient to explain the disk dispersal, assuming that the contribution of UV photoevaporating stellar winds to radio flux does not evolve with time. The sensitivity of our data cannot exclude photoevaporation due to X-ray photons as an efficient mechanism for disk dispersal. Deeper surveys with the Square Kilometre Array will be able to provide strong constraints on disk photoevaporation., Comment: Accepted in A&A

Phosphorus nitride (PN) is the P-bearing species with the highest number of detections in star-forming regions. Multi-line studies of the molecule have shown that the excitation temperature of PN is usually lower than the gas kinetic temperature, suggesting that PN is likely in conditions of sub-thermal excitation. We present an analysis of PN which takes the possible sub-thermal excitation conditions into account in a sample of 24 massive star-forming regions. We observed PN (2-1), (3-2), (4-3), and (6-5) with the IRAM-30m and APEX telescopes and detected PN lines in 15 of them. Together with 9 similar sources detected in PN in previous works, we have analysed the largest sample of star-forming regions to date, made of 33 sources with 24 detections in total (among which 13 are new detections). Hence, we have increased the number of star-forming regions detected in PN by more than a factor 2. Our analysis indicates that the PN lines are indeed sub-thermally excited, but well described by a single excitation temperature. We have compared line profiles and fractional abundances of PN and SiO, a typical shock tracer, and found that almost all objects detected in PN have high-velocity SiO wings. Moreover, the SiO and PN abundances with respect to H2 are correlated over several orders of magnitude, and uncorrelated with gas temperature. This clearly shows that the production of PN is strongly linked to the presence of shocked gas, and rules out alternative scenarios based on thermal evaporation from iced grain mantles., Comment: 14 pages, 10 figures, accepted for publication in MNRAS

Nitrogen (N) fractionation is used as a tool to search for a link between the chemical history of the Solar System and star-forming regions. A large variation of $^{14}$N/$^{15}$N is observed towards different astrophysical sources, and current chemical models cannot reproduce it. With the advent of high angular resolution radiotelescopes it is now possible to search for N-fractionation at core scales. We present IRAM NOEMA observations of the J=1-0 transition of N$_{2}$H$^{+}$, $^{15}$NNH$^{+}$ and N$^{15}$NNH$^{+}$ towards the high-mass protocluster IRAS 05358+3543. We find $^{14}$N/$^{15}$N ratios that span from $\sim$100 up to $\sim$220 and these values are lower or equal than those observed with single-dish observations towards the same source. Since N-fractionation changes across the studied region, this means that it is regulated by local environmental effects. We find also the possibility, for one of the four cores defined in the protocluster, to have a more abundant $^{15}$NNH$^{+}$ with respect to N$^{15}$NNH$^{+}$. This is another indication that current chemical models may be missing chemical reactions or may not take into account other mechanisms, like photodissociation or grain surface chemistry, that could be important., Comment: 19 pages, 8 figures, 6 tables, 3 appendices Accepted in Monthly Notices of the Royal Astronomical Society Letters

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

Interstellar CO2 is an important reservoir of carbon and oxygen, and one of the major constituents of the icy mantles of dust grains, but it is not observable directly in the cold gas because has no permanent dipole moment. Its protonated form, HOCO+, is believed to be a good proxy for gaseous CO2. However, it has been detected in only a few star-forming regions so far, so that its interstellar chemistry is not well understood. We present new detections of HOCO+ lines in 11 high-mass star-forming clumps. Our observations increase by more than three times the number of detections in star-forming regions so far. We have derived beam-averaged abundances relative to H2 in between 0.3 and 3.8 x 10^{-11}. We have compared these values with the abundances of H13CO+, a possible gas-phase precursor of HOCO+, and CH3OH, a product of surface chemistry. We have found a positive correlation with H13CO+, while with CH3OH there is no correlation. We suggest that the gas-phase formation route starting from HCO+ plays an important role in the formation of HOCO+, perhaps more relevant than protonation of CO2 (upon evaporation of this latter from icy dust mantles)., Comment: 5 pages, 4 figures, 1 table, accepted for publication in MNRAS