Search

The stellar initial mass function (IMF) is critical to our understanding of star formation and the effects of young stars on their environment. On large scales, it enables us to use tracers such as UV or Halpha emission to estimate the star formation rate of a system and interpret unresolved star clusters across the universe. So far, there is little firm evidence of large-scale variations of the IMF, which is thus generally considered universal. Stars form from cores and it is now possible to estimate core masses and compare the core mass function (CMF) with the IMF, which it presumably produces. The goal of the ALMA-IMF large program is to measure the core mass function at high linear resolution (2700 au) in 15 typical Milky Way protoclusters spanning a mass range of 2500 to 32700 Msun. In this work, we used two different core extraction algorithms to extract about 680 gravitationally bound cores from these 15 protoclusters. We adopt per core temperature using the temperature estimate from the PPMAP Bayesian method. A power-law fit to the CMF of the sub-sample of cores above the 1.64 Msun completeness limit, 330 cores, through the maximum likelihood estimate technique yields a slope of 1.97 +/- 0.06, significantly flatter than the 2.35 Salpeter slope. Assuming a self-similar mapping between the CMF and the IMF, this result implies that these 15 high-mass protoclusters will generate atypical IMFs. This sample is the largest to date produced and analysed self-consistently, derived at matched physical resolution, with per-core temperature estimates and cores as massive as 150 Msun. We provide the raw source extraction catalogues and the source derived size, temperature, mass, and spectral indices in the 15 protoclusters., Comment: Accepted by 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

Sites associated with high-mass star and cluster formation exhibit a so-called hot core phase, characterized by high temperatures and column densities of complex organic molecules. We built a comprehensive census of hot core candidates towards the ALMA-IMF protoclusters based on the detection of two CH3OCHO emission lines at 216.1 GHz. We used the source extraction algorithm GExt2D to identify peaks of methyl formate (CH3OCHO) emission that is a complex species commonly observed towards sites of star formation. We built up a catalog of 76 hot core candidates with masses ranging from about 0.2 to 80 Msun , of which 56 are new detections. A large majority of these objects are compact, rather circular, with deconvolved FWHM sizes of about 2300 au on average. About 30% of our sample of methyl formate sources have core masses above 8 Msun within sizes ranging from about 1000 au to 13400 au, which well correspond to archetypical hot cores. The origin of the CH3OCHO emission toward the lower-mass cores can be explained by a mixture of contribution from shocks, or may correspond to objects in a more evolved state, i.e. beyond the hot core stage. We find that the fraction of hot core candidates increases with the core mass. The large fraction of hot core candidates towards the most massive cores suggests that they rapidly enter the hot core phase and feedback effects from the forming protostar(s) impact their environment on short time-scales.

The origin of the stellar initial mass function (IMF) and its relation with the core mass function (CMF) are actively debated issues with important implications in astrophysics. Recent observations in the W43 molecular complex of top-heavy CMFs, with an excess of high-mass cores compared to the canonical mass distribution, raise questions about our understanding of the star formation processes and their evolution in space and time. We aim to compare populations of protostellar and prestellar cores in three regions imaged in the ALMA-IMF Large Program. We created an homogeneous core catalogue in W43, combining a new core extraction in W43-MM1 with the catalogue of W43-MM2&MM3 presented in a previous work. Our detailed search for protostellar outflows enabled us to identify between 23 and 30 protostellar cores out of 127 cores in W43-MM1 and between 42 and 51 protostellar cores out of 205 cores in W43-MM2&MM3. Cores with neither outflows nor hot core emission are classified as prestellar candidates. We found a similar fraction of cores which are protostellar in the two regions, about 35%. This fraction strongly varies in mass, from 15-20% at low mass, between 0.8 and 3$M_{\odot} $ up to about 80% above 16$M_{\odot}$. Protostellar cores are found to be, on average, more massive and smaller in size than prestellar cores. Our analysis also revealed that the high-mass slope of the prestellar CMF in W43, $\alpha=-1.46_{-0.19}^{+0.12}$, is consistent with the Salpeter slope, and thus the top-heavy form measured for the global CMF, $\alpha=-0.96$, is due to the protostellar core population. Our results could be explained by clump-fed models in which cores grow in mass, especially during the protostellar phase, through inflow from their environment. The difference between the slopes of the prestellar and protostellar CMFs moreover implies that high-mass cores grow more in mass than low-mass cores., Comment: 28 pages, 18 figures. Accepted for publication in A&A. Title and appendices updated

We have used the axisymmetric chemo-hydrodynamical code WALKIMYA-2D to numerically model and reproduce the physical and CO emission properties of the jet-driven outflow from the intermediate-mass protostar Cep E, which was observed at $\sim 800$au resolution in the CO $J=2\to 1$ line with the IRAM interferometer. Our simulations take into account the observational constraints available on the physical structure of the protostellar envelope to provide constraints on the dynamics of the inner protostellar environment from the study of the outflow/jet propagation away from the launch region. WALKIMYA-2D successfully reproduces the main qualitative and quantitative features of the Cep E outflow and the jet kinematics, naturally accounting for their time variability. Signatures of internal shocks are detected as knots along the jet. In the early times of the ejection process, the young emitted knots interact with the dense circumstellar envelope through high-velocity, dissociative shocks, which strongly decrease the CO gas abundance in the jet. As time proceeds, the knots propagate more smoothly through the envelope and dissociative shocks disappear after $\sim 10^3$ yr. The distribution of CO abundance along the jet shows that the latter bears memory of the early dissociative phase in the course of its propagation. Analysis of the velocity field shows that the jet material mainly consists of gas entrained from the circumstellar envelope and accelerated away from the protostar at $700$ au scale. As a result, the overall jet mass loss rate appears higher than the actual mass ejection rate by a factor $\sim 3$. Numerical modeling of the Cep E jet-driven outflow and comparison with the CO observations have allowed us to peer into the outflow formation mechanism with unprecedented detail and to retrieve the history of the mass-loss events that have shaped the outflow.

Among the most central open questions regarding the initial mass function (IMF) of stars is the impact of environment on the shape of the core mass function (CMF) and thus potentially on the IMF. The ALMA-IMF Large Program aims to investigate the variations in the core distributions with cloud characteristics, as diagnostic observables of the formation process and evolution of clouds. The present study focuses on the W43-MM2&MM3 mini-starburst, whose CMF has recently been found to be top-heavy with respect to the Salpeter slope. W43-MM2&MM3 harbors a rich cluster that contains a statistically significant number of cores, which was previously characterized in Paper III. We applied a multi-scale decomposition technique to the ALMA 1.3 mm and 3 mm continuum images to define six subregions. For each subregion we characterized the high column density probability distribution function, n-PDF, and the shape of the cloud gas using the 1.3 mm image. Using the core catalog, we investigate correlations between the CMF and cloud and core properties. We classify the subregions into different stages of evolution, from quiescent to burst to post-burst, based on the surface number density of cores, number of outflows, and UCHii presence. The high-mass end of the subregion CMFs varies from being close to the Salpeter slope (quiescent) to top-heavy (burst and post-burst). Moreover, the second tail of the n-PDF varies from steep, to flat like observed for the high mass star-forming clouds. We found that subregions with flat second n-PDF tails display top-heavy CMFs. The CMF may evolve from Salpeter to top-heavy throughout the star formation process from the quiescent to the burst phase. This scenario raises the question of if the CMF might revert again to Salpeter as the cloud approaches the end of its star formation stage, a hypothesis that remains to be tested., Comment: Accepted for publication in A&A (december, 5th 2022) after language editing; 25 pages, 14 figures

The Large Program "Astrochemical Surveys At IRAM" (ASAI) investigates the emergence of molecular complexity along the different stages of the solar-type star formation process, by carrying out unbiased line surveys of a sample of ten template sources in the range 80-272 GHz with the IRAM 30m telescope. We present here an overview of the main results of the Large Program ASAI., Comment: 6 pages, 4 figures; invited talk at the IRAM conference "Multi-line Diagnostics of the Interstellar Medium", Nice, April 2022; To be published in The European Physical Journal Web of Conferences

W43-MM1 is a young region, very rich in terms of high-mass star formation. We aim to systematically identify the massive cores which contain a hot core and compare their molecular composition. We used ALMA high-spatial resolution (2500 au) data of W43-MM1 to identify line-rich protostellar cores and make a comparative study of their temperature and molecular composition. The identification of hot cores is based on both the spatial distribution of the complex organic molecules and the contribution of molecular lines relative to the continuum intensity. We rely on the analysis of CH3CN and CH3CCH to estimate the temperatures of the selected cores. Finally, we rescale the spectra of the different hot cores based on their CH3OCHO line intensities to directly compare the detections and line intensities of the other species. W43-MM1 turns out to be a region rich in massive hot cores with at least 1 less massive and 7 massive hot cores. The excitation temperature of CH3CN is of the same order for all of them (120-160 K). There is a factor of up to 30 difference in the intensity of the complex organic molecules (COMs) lines. However the molecular emission of the hot cores appears to be the same within a factor 2-3. This points towards both a similar chemical composition and excitation of most of the COMs over these massive cores, which span about an order of magnitude in core mass. In contrast, CH3CCH emission is found to preferentially trace more the envelope, with a temperature ranging from 50 K to 90 K. Core 1, the most massive hot core of W43-MM1, shows a richer line spectrum than the other cores. In core 2, the emission of O-bearing molecules does not peak at the dust continuum core center; the blue and red shifted emission correspond to the outflow lobes, suggesting a formation via the sublimation of the ice mantles through shocks or UV irradiation on the walls of the cavity., Comment: 19 pages, 17 figures. Accepted for publication in A&A

The ALMA interferometer, with its unprecedented combination of high-sensitivity and high-angular resolution, allows for (sub-)mm wavelength mapping of protostellar systems at Solar System scales. Astrochemistry has benefited from imaging interstellar complex organic molecules in these jet-disk systems. Here we report the first detection of methanol (CH3OH) and methyl formate (HCOOCH3) emission towards the triple protostellar system VLA1623-2417 A1+A2+B, obtained in the context of the ALMA Large Program FAUST. Compact methanol emission is detected in lines from Eu = 45 K up to 61 K and 537 K towards components A1 and B, respectively. LVG analysis of the CH3OH lines towards VLA1623-2417 B indicates a size of 0.11-0.34 arcsec (14-45 au), a column density N(CH3OH) = 10^16-10^17 cm-2, kinetic temperature > 170 K, and volume density > 10^8 cm-3. An LTE approach is used for VLA1623-2417 A1, given the limited Eu range, and yields Trot < 135 K. The methanol emission around both VLA1623-2417 A1 and B shows velocity gradients along the main axis of each disk. Although the axial geometry of the two disks is similar, the observed velocity gradients are reversed. The CH3OH spectra from B shows two broad (4-5 km s-1) peaks, which are red- and blue-shifted by about 6-7 km s-1 from the systemic velocity. Assuming a chemically enriched ring within the accretion disk, close to the centrifugal barrier, its radius is calculated to be 33 au. The methanol spectra towards A1 are somewhat narrower (about 4 km s-1), implying a radius of 12-24 au., Comment: MNRAS, in press

Planetary systems such as our own are formed after a long process where matter condenses from diffuse clouds to stars, planets, asteroids, comets and residual dust, undergoing dramatic changes in physical and chemical state in less than a few million years. Several studies have shown that the chemical composition during the early formation of a Solar-type planetary system is a powerful diagnostic to track the history of the system itself. Among the approximately 270 molecules so far detected in the ISM, the so-called interstellar complex organic molecules (iCOMs) are of particular interest both because of their evolutionary diagnostic power and because they might be potential precursors of biomolecules, which are at the basis of terrestrial life. This Chapter focuses on the evolution of organic molecules during the early stages of a Solar-type planetary system, represented by the prestellar, Class 0/I and protoplanetary disk phases, and compares them with what is observed presently in Solar System comets. Our twofold goal is to review the processes at the base of organic chemistry during Solar-type star formation and, in addition, to possibly provide constraints on the early history of our own planetary system., Comment: Chapter for Protostars & Planets VII

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

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 ALMA-IMF Large Program imaged a total noncontiguous area of 53pc2, covering 15 extreme, nearby protoclusters of the Milky Way. They were selected to span relevant early protocluster evolutionary stages. Our 1.3mm and 3mm observations provide continuum images that are homogeneously sensitive to point-like cores with masses of 0.2 and 0.6Msun, respectively, with a matched spatial resolution of 2000au. We also detect lines that probe the protocluster structure, kinematics, chemistry, and feedback over scales from clouds to filaments to cores. We classify ALMA-IMF protoclusters as Young, Intermediate, or Evolved based on the amount of dense gas in the cloud that has potentially been impacted by HII regions. The ALMA-IMF catalog contains 700 cores that span a mass range of 0.15-250Msun at a typical size of 2100au. We show that this core sample has no significant distance bias and can be used to build core mass functions at similar physical scales. Significant gas motions, which we highlight here in the G353.41 region, are traced down to core scales and can be used to look for inflowing gas streamers and to quantify the impact of the possible associated core mass growth on the shape of the CMF with time. Our first analysis does not reveal any significant evolution of the matter concentration from clouds to cores or from the youngest to more evolved protoclusters, indicating that cloud dynamical evolution and stellar feedback have for the moment only had a slight effect on the structure of high-density gas in our sample. Furthermore, the first-look analysis of the line richness toward bright cores indicates that the survey encompasses several tens of hot cores, of which we highlight the most massive in the G351.77 cloud. Their homogeneous characterization can be used to constrain the emerging molecular complexity in protostars of high to intermediate masses., Comment: 27 pages, 9 figures, accepted for publication by A&A

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

HCN and its isomer HNC play an important role in molecular cloud chemistry and the formation of more complex molecules. We investigate here the impact of protostellar shocks on the HCN and HNC abundances from high-sensitivity IRAM 30m observations of the prototypical shock region L1157-B1 and the envelope of the associated Class 0 protostar, as a proxy for the pre-shock gas. The isotopologues H$^{12}$CN, HN$^{12}$C, H$^{13}$CN, HN$^{13}$C, HC$^{15}$N, H$^{15}$NC, DCN and DNC were all detected towards both regions. Abundances and excitation conditions were obtained from radiative transfer analysis of molecular line emission under the assumption of Local Thermodynamical Equilibrium. In the pre-shock gas, the abundances of the HCN and HNC isotopologues are similar to those encountered in dark clouds, with a HCN/HNC abundance ratio $\approx 1$ for all isotopologues. A strong D-enrichment (D/H$\approx 0.06$) is measured in the pre-shock gas. There is no evidence of $^{15}$N fractionation neither in the quiescent nor in the shocked gas. At the passage of the shock, the HCN and HNC abundances increase in the gas phase in different manners so that the HCN/HNC relative abundance ratio increases by a factor 20. The gas-grain chemical and shock model UCLCHEM allows us to reproduce the observed trends for a C-type shock with pre-shock density $n$(H)= $10^5$cm$^{-3}$ and shock velocity $V_s= 40$km/s. We conclude that the HCN/HNC variations across the shock are mainly caused by the sputtering of the grain mantle material in relation with the history of the grain ices., Comment: 14 pages, 7 figures. Accepted for publication in MNRAS

Using the Yebes 40m and IRAM 30m radiotelescopes, we detected two series of harmonically related lines in space that can be fitted to a symmetric rotor. The lines have been seen towards the cold dense cores TMC-1, L483, L1527, and L1544. High level of theory ab initio calculations indicate that the best possible candidate is the acetyl cation, CH3CO+, which is the most stable product resulting from the protonation of ketene. We have produced this species in the laboratory and observed its rotational transitions Ju = 10 up to Ju = 27. Hence, we report the discovery of CH3CO+ in space based on our observations, theoretical calculations, and laboratory experiments. The derived rotational and distortion constants allow us to predict the spectrum of CH3CO+ with high accuracy up to 500 GHz. We derive an abundance ratio N(H2CCO)/N(CH3CO+) = 44. The high abundance of the protonated form of H2CCO is due to the high proton affinity of the neutral species. The other isomer, H2CCOH+, is found to be 178.9 kJ/mol above CH3CO+. The observed intensity ratio between the K=0 and K=1 lines, 2.2, strongly suggests that the A and E symmetry states have suffered interconversion processes due to collisions with H and/or H2, or during their formation through the reaction of H3+ with H2CCO., Comment: Accepted for publication in A&A Letters

As a part of the CALYPSO large programme, we constrain the properties of protostellar jets and outflows in a sample of 21 Class 0 protostars with internal luminosities, Lint, from 0.035 to 47 Lsun. We analyse high angular resolution (~0.5"-1") IRAM PdBI observations in CO (2-1), SO ($5_6-4_5$), and SiO (5-4). CO (2-1), which probes outflowing gas, is detected in all the sources (for the first time in SerpS-MM22 and SerpS-MM18b). Collimated high-velocity jets in SiO (5-4) are detected in 67% of the sources (for the first time in IRAS4B2, IRAS4B1, L1448-NB, SerpS-MM18a), and 77% of these also show jet/outflow emission in SO ($5_6-4_5$). In 5 sources (24% of the sample) SO ($5_6-4_5$) probes the inner envelope and/or the disk. The CALYPSO survey shows that the outflow phenomenon is ubiquitous and that the detection rate of high-velocity jets increases with protostellar accretion, with at least 80% of the sources with Lint>1 Lsun driving a jet. The protostellar flows exhibit an onion-like structure, where the SiO jet (opening angle ~10$^o$) is nested into a wider angle SO (~15$^o$) and CO (~25$^o$) outflow. On scales >300 au the SiO jets are less collimated than atomic jets from Class II sources (~3$^o$). Velocity asymmetry between the two jet lobes are detected in one third of the sources, similarly to Class II atomic jets, suggesting that the same launching mechanism is at work. Most of the jets are SiO rich (SiO/H2 from >2.4e-7 to >5e-6), which indicates efficient release of >1%-10% of silicon in gas phase likely in dust-free winds, launched from inside the dust sublimation radius. The mass-loss rates (from ~7e-8 to ~3e-6 Msun/yr) are larger than what was measured for Class II jets. Similarly to Class II sources, the mass-loss rates are ~1%-50% of the mass accretion rates suggesting that the correlation between ejection and accretion in young stars holds from 1e4 yr up to a few Myr., Comment: Accepted for publication by A&A, 40 pages, 7 figures + appendixes

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

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

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

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

The accretion history of protostars remains widely mysterious even though it represents one of the best ways to understand the protostellar collapse that leads to the formation of stars. Molecular outflows are here used to characterize the protostellar accretion phase in W43-MM1. The W43-MM1 protocluster host a sufficient number of protostars to statistically investigate molecular outflows in a single, homogeneous region. We used the CO(2-1) and SiO(5-4) line datacubes, taken as part of an ALMA mosaic with a 2000 AU resolution, to search for protostellar outflows, evaluate the influence that the environment has on these outflows' characteristics and put constraints on outflow variability in W43-MM1. We discovered a rich cluster of 46 outflow lobes, driven by 27 protostars with masses of 1-100 Msun. The complex environment inside which these outflow lobes develop has a definite influence on their length, limiting the validity of using outflows' dynamical timescales as a proxy of the ejection timescale in clouds with high dynamics and varying conditions. We performed a detailed study of Position-Velocity (PV) diagrams of outflows that revealed clear events of episodic ejection. The time variability of W43-MM1 outflows is a general trend and is more generally observed than in nearby, low- to intermediate-mass star-forming regions. The typical timescale found between two ejecta, about 500 yr, is consistent with that found in nearby protostars. If ejection episodicity reflects variability in the accretion process, either protostellar accretion is more variable or episodicity is easier to detect in high-mass star-forming regions than in nearby clouds. The timescale found between accretion events could be resulting from disk instabilities, associated with bursts of inflowing gas arising from the dynamical environment of high-mass star-forming cores., Comment: Accepted for publication in Astronomy & Astrophysics (23 pages, 15 figures)

Complex organic molecules (COMs) have been detected in a few Class 0 protostars but their origin is not well understood. Going beyond studies of individual objects, we want to investigate the origin of COMs in young protostars on a statistical basis. We use the CALYPSO survey performed with the IRAM PdBI to search for COMs at high angular resolution in a sample of 26 solar-type protostars, including 22 Class 0 and four Class I objects. Methanol is detected in 12 sources and tentatively in one source, which represents half of the sample. Eight sources (30%) have detections of at least three COMs. We find a strong chemical differentiation in multiple systems with five systems having one component with at least three COMs detected but the other component devoid of COM emission. The internal luminosity is found to be the source parameter impacting the most the COM chemical composition of the sources, while there is no obvious correlation between the detection of COM emission and that of a disk-like structure. A canonical hot-corino origin may explain the COM emission in four sources, an accretion-shock origin in two or possibly three sources, and an outflow origin in three sources. The CALYPSO sources with COM detections can be classified into three groups on the basis of the abundances of oxygen-bearing molecules, cyanides, and CHO-bearing molecules. These chemical groups correlate neither with the COM origin scenarii, nor with the evolutionary status of the sources if we take the ratio of envelope mass to internal luminosity as an evolutionary tracer. We find strong correlations between molecules that are a priori not related chemically (for instance methanol and methyl cyanide), implying that the existence of a correlation does not imply a chemical link. [abridged], Comment: Accepted for publication in Astronomy and Astrophysics

Aim: In the past, observations of protostellar shocks have been able to set constraints on the formation route of formamide (NH2CHO), exploiting its observed spatial distribution and comparison with astrochemical model predictions. In this work, we follow the same strategy to study the case of acetaldehyde (CH3CHO). Method: To this end, we used the data obtained with the IRAM-NOEMA interferometer in the framework of the Large Program SOLIS to image the B0 and B1 shocks along the L1157 blueshifted outflow in methanol (CH3OH) and acetaldehyde line emission. Results: We imaged six CH3OH and eight CH3CHO lines which cover upper level energies up to 30 K. Both species trace the B0 molecular cavity as well as the northern B1 portion, i.e. the regions where the youngest shocks (1000 yr) occurred. The CH$_3$OH and CH$_3$CHO emission peaks towards the B1b clump, where we measured the following column densities and relative abundances: 1.3 x 10^16 cm-2 and 6.5 x 10-6 (methanol), and 7 x 10^13 cm-2 and 3.5 x 10-8 (acetaldehyde). We carried out a non-LTE LVG analysis of the observed CH3OH line: the average kinetic temperature and density of the emitting gas are Tkin = 90 K and nH2 = 4 x 10^5 cm-3, respectively. The CH3OH and CH3CHO abundance ratio towards B1b is 190, varying by less than a factor 3 throughout the whole B0-B1 structure. Conclusions: The comparison of astrochemical model predictions with the observed methanol and acetaldehyde spatial distribution does not allow to distinguish whether acetaldehyde is formed on the grain mantles or rather on the gas-phase, as its gas-phase formation, dominated by the reaction of ethyl radical (CH3CH2) with atomic oxygen, is very fast. Observations of acetaldehyde in younger shocks, e.g. 10^2 yr old, or/and of the ethyl radical, whose frequencies are not presently available, are necessary to settle the issue., Comment: A&A, in press

The chemical composition of protostellar jets and its origin are still badly understood. More observational constraints are needed to make progress. With that objective, we have carried out a systematic search for molecular species in the jet of Cep E-mm, a template for intermediate-mass Class 0 protostars, associated with a luminous, high-velocity outflow. We made use of an unbiased spectral line survey in the range 72-350 GHz obtained with the IRAM 30m telescope, complementary observations of the CO $J$=3-2 transition with the JCMT, and observations at 1" angular resolution of the CO $J$=2-1 transition with the IRAM Plateau de Bure interferometer. In addition to CO, we have detected rotational transitions from SiO, SO, H$_2$CO, CS, HCO$^{+}$ and HCN. A strong chemical differentiation is observed in the southern and northern lobes of the jet. Radiative transfer analysis in the Large Velocity Gradient approximation yields typical molecular abundances of the order of $10^{-8}$ for all molecular species other than CO. Overall, the jets exhibit an unusual chemical composition, as CS, SO and H$_2$CO are found to be the most abundant species, with a typical abundance of (3-4)$\times 10^{-8}$. The transverse size of the CO jet emission estimated from interferometric observations is about 1000 au, suggesting that we are detecting emission from a turbulent layer of gas entrained by the jet in its propagation and not the jet itself. We propose that some molecular species could be the signatures of the specific photochemistry driven by the UV radiation field generated in the turbulent envelope., Comment: 14 pages, 11 figures, accepted to publication on MNRAS

We present the first census of the interstellar Complex Organic Molecules (iCOMs) in the low-mass Class I protostar SVS13-A, obtained by analysing data from the IRAM-30m Large Project ASAI (Astrochemical Surveys At IRAM). They consist of an high-sensitivity unbiased spectral survey at the 1mm, 2mm and 3mm IRAM bands. We detected five iCOMs: acetaldehyde (CH$_3$CHO), methyl formate (HCOOCH$_3$), dimethyl ether (CH$_3$OCH$_3$), ethanol (CH$_3$CH$_2$OH) and formamide (NH$_2$CHO). In addition we searched for other iCOMs and ketene (H$_2$CCO), formic acid (HCOOH) and methoxy (CH$_3$O), whose only ketene was detected. The numerous detected lines, from 5 to 37 depending on the species, cover a large upper level energy range, between 15 and 254 K. This allowed us to carry out a rotational diagram analysis and derive rotational temperatures between 35 and 110 K, and column densities between $3\times 10^{15}$ and $1\times 10^{17}$ cm$^{-2}$ on the 0."3 size previously determined by interferometric observations of glycolaldehyde. These new observations clearly demonstrate the presence of a rich chemistry in the hot corino towards SVS13-A. The measured iCOMs abundances were compared to other Class 0 and I hot corinos, as well as comets, previously published in the literature. We find evidence that (i) SVS13-A is as chemically rich as younger Class 0 protostars, and (ii) the iCOMs relative abundances do not substantially evolve during the protostellar phase., Comment: 24 pages, MNRAS in press

Understanding the formation mechanisms of protoplanetary disks and multiple systems, and their pristine properties, is a key question for modern astrophysics. The properties of the youngest disks, embedded in rotating infalling protostellar envelopes, have largely remained unconstrained up to now. In the framework of the IRAM-PdBI CALYPSO survey, we have obtained sub-arcsecond observations of the dust continuum emission at 231 GHz and 94 GHz, for a sample of 16 solar-type Class 0 protostars. In an attempt to identify disk-like structures embedded at small scales in the protostellar envelopes, we model the dust continuum emission visibility profiles using both Plummer-like envelope models and envelope models including additional Gaussian disk-like components. Our analysis shows that in the CALYPSO sample, 11 of the 16 Class 0 protostars are better reproduced by models including a disk-like dust continuum component contributing to the flux at small scales, but less than 25% of these candidate protostellar disks are resolved at radii > 60 au. Including all available literature constraints on Class 0 disks at subarcsecond scales, we show that our results are representative: most (> 72% in a sample of 26 protostars) Class 0 protostellar disks are small and emerge only at radii < 60 au. Our multiplicity fraction at scales 100-5000 au is in global agreement with the multiplicity properties of Class I protostars at similar scales. We confront our observational constraints on the disk size distribution in Class 0 protostars to the typical disk properties from protostellar formation models. Because they reduce the centrifugal radius, and produce a disk size distribution peaking at radii <100 au during the main accretion phase, magnetized models of rotating protostellar collapse are favored by our observations., Comment: Accepted for publication in Astronomy \& Astrophysics. 24 pages of main text, 22 pages of additional content (under the form of three appendices) - revised version including language and typos corrections

Intermediate-mass (IM) protostars provide a bridge between the low- and high-mass protostars. Despite their relevance, little is known about their chemical diversity. We want to investigate the molecular richness towards the envelope of IM protostars and to compare their properties with other sources. We have selected the isolated IM Class 0 protostar CepE-mm to carry out an unbiased molecular survey with the IRAM 30m telescope between 72 and 350 GHz with an angular resolution lying in the range 7-34". These data were complemented with NOEMA observations of the spectral bands 85.9 - 89.6 GHz and 216.8 - 220.4 GHz at an angular resolution of 2.3 and 1.4", respectively. 30m spectra show bright emission of O- and N-bearing COMs. We identify three components in the spectral signature: an extremely broad line component associated with the outflowing gas, a narrow line component associated with the cold envelope, and a broad line component tracing the hot corino. NOEMA observations reveal CepE-mm as a binary protostellar system CepE-A and B, separated by 1.7". CepE-A dominates the core continuum emission and powers the high-velocity jet associated with HH377. Our interferometric maps show that COMs arises from a region of 0.7" size around CepE-A. Rotational temperatures were found to lie in the range 20-40 K with column densities ranging from a few 10^15 cm-2 for O-bearing species, down to a few 10^14 cm-2 for N-bearing species. Molecular abundances are similar to those measured towards other low- and intermediate-mass protostars. High-mass hot cores are significantly less abundant in methanol and N-bearing species are more abundant by 3 orders of magnitude.CepE-mm reveals itself as a binary protostellar system with a strong chemical differentiation. The brightest component of the binary is associated with a hot corino. Its properties are similar to those of low-mass hot corinos.

The HH 212 protostellar system, in Orion B, has been mapped thanks to ALMA-Band 7 Cycle 1 and Cycle 4 observations of dueterated water (HDO) and acetaldehyde (CH$_3$CHO) emission with an angular resolution down to $\sim$0.15 arcsec (60 au). Many emission lines due to 14 CH$_3$CHO and 1 HDO transitions at high excitation ($E_{\rm u}$ between 163 K and 335 K) have been imaged in the inner $\sim$ 70 au region. The local thermal equilibrium analysis of the CH$_3$CHO emission leads to a temperature of 78$\pm$14 K and a column density of 7.6$\pm$3.2 $\times$ 10$^{15}$ cm$^{-2}$, which, when $N_{\rm H_2}$ of 10$^{24}$ cm$^{-2}$ is assumed, leads to an abundance of $X_{\rm CH_3CHO}$ $\simeq$ 8 $\times$ 10$^{-9}$. The large velocity gradient analysis of the HDO emission also places severe constraints on the volume density, n$_{\rm H_2}$ $\geq$ 10$^8$ cm$^{-3}$. The line profiles are 5--7 km s$^{-1}$ wide, and CH$_3$CHO and HDO both show a $\pm$ 2 km s$^{-1}$ velocity gradient over a size of $\sim$ 70 au (blue-shifted emission towards the north-west and red-shifted emission towards the south-east) along the disc equatorial plane, in agreement with what was found so far using other molecular tracers. The kinematics of CH$_3$CHO and HDO are consistent with the occurrence of a centrifugal barrier, that is, the infalling envelope-rotating disc ring, which is chemically enriched through low-velocity accretion shocks. The emission radius is $\sim$ 60 au, in good agreement with what was found before for another interstellar complex organic molecule such as NH$_2$CHO. We support a vertical structure for the centrifugal barrier, suggesting the occurrence of two outflowing, expanding, and rotating rings above and below (of about 40-45 au) the optically thick equatorial disc plane. It is tempting to speculate that these rings could probe the basis of a wind launched from this region., Comment: Astronomy & Astrophysics, in press

The origin of the heavily fractionated reservoir of nitrogen in comets remains an issue in the theory of their formation and hence of the solar system. Whether the fractionated reservoir traced by comets is inherited from the interstellar cloud or is the product of processes taking place in the protostar, or in the protoplanetary disk, remains unclear. So far, observations of nitrogen isotopic ratios in protostars or prestellar cores have not securely identified such a fractionated reservoir owing to the intrinsic difficulty of direct isotopic ratios measurements. In this article, we report the detection of 5 rotational lines of \ce{HC3N}, {including the weaker components of the hyperfine multiplets}, and two rotational lines of its $^{15}$N isotopologue, towards the L1544 prestellar core. Based on a MCMC/non-LTE multi-line analysis at the hyperfine level, we derive the column densities of \ce{HC3N} ($8.0\pm0.4\tdix{13}$\cc) and \ce{HC3^{15}N} ($2.0\pm0.4\tdix{11}$\cc) and derive an isotopic ratio of 400$\pm$20(1$\sigma$). This value suggests that \ce{HC3N} is slightly depleted in $^{15}$N in L1544 with respect to the elemental $^{14}$N/$^{15}$N ratio {of $\approx$330} in the present-day local interstellar medium. Our study also stresses the need for radiative calculations at the hyperfine level. Finally, the comparison of the derived ratio with those obtained in CN and HCN in the same core seems to favor CN+C$_2$H$_2$ as the dominant formation route to HC$_3$N. However, uncertainties in the isotopic ratios preclude definitive conclusions., Comment: Accepted for publication in MNRAS

Evidence is mounting that the small bodies of our Solar System, such as comets and asteroids, have at least partially inherited their chemical composition from the first phases of the Solar System formation. It then appears that the molecular complexity of these small bodies is most likely related to the earliest stages of star formation. It is therefore important to characterize and to understand how the chemical evolution changes with solar-type protostellar evolution. We present here the Large Program "Astrochemical Surveys At IRAM" (ASAI). Its goal is to carry out unbiased millimeter line surveys between 80 and 272 GHz of a sample of ten template sources, which fully cover the first stages of the formation process of solar-type stars, from prestellar cores to the late protostellar phase. In this article, we present an overview of the surveys and results obtained from the analysis of the 3 mm band observations. The number of detected main isotopic species barely varies with the evolutionary stage and is found to be very similar to that of massive star-forming regions. The molecular content in O- and C- bearing species allows us to define two chemical classes of envelopes, whose composition is dominated by either a) a rich content in O-rich complex organic molecules, associated with hot corino sources, or b) a rich content in hydrocarbons, typical of Warm Carbon Chain Chemistry sources. Overall, a high chemical richness is found to be present already in the initial phases of solar-type star formation., Comment: 19 pages, 7 figures, 5 Tables

Although chemical models predict that the deuterium fractionation in N$_2$H$^+$ is a good evolutionary tracer in the star formation process, the fractionation of nitrogen is still a poorly understood process. Recent models have questioned the similar evolutionary trend expected for the two fractionation mechanisms in N$_2$H$^+$, based on a classical scenario in which ion-neutral reactions occurring in cold gas should have caused an enhancement of the abundance of N$_2$D$^+$, $^{15}$NNH$^+$, and N$^{15}$NH$^+$. In the framework of the ASAI IRAM-30m large program, we have investigated the fractionation of deuterium and $^{15}$N in N$_2$H$^+$ in the best known representatives of the different evolutionary stages of the Sun-like star formation process. The goal is to ultimately confirm (or deny) the classical "ion-neutral reactions" scenario that predicts a similar trend for D and $^{15}$N fractionation. We do not find any evolutionary trend of the $^{14}$N/$^{15}$N ratio from both the $^{15}$NNH$^+$ and N$^{15}$NH$^+$ isotopologues. Therefore, our findings confirm that, during the formation of a Sun-like star, the core evolution is irrelevant in the fractionation of $^{15}$N. The independence of the $^{14}$N/$^{15}$N ratio with time, found also in high-mass star-forming cores, indicates that the enrichment in $^{15}$N revealed in comets and protoplanetary disks is unlikely to happen at core scales. Nevertheless, we have firmly confirmed the evolutionary trend expected for the H/D ratio, with the N$_2$H$^+$/N$_2$D$^+$ ratio decreasing before the pre-stellar core phase, and increasing monotonically during the protostellar phase. We have also confirmed clearly that the two fractionation mechanisms are not related., Comment: 9 pages, 2 figures, accepted for publication in MNRAS

We report the detection in space of a new molecular species which has been characterized spectroscopically and fully identified from astrophysical data. The observations were carried out with the 30m IRAM telescope. The molecule is ubiquitous as its $J$=2$\rightarrow$1 transition has been found in cold molecular clouds, prestellar cores, and shocks. However, it is not found in the hot cores of Orion-KL and in the carbon-rich evolved star IRC+10216. Three rotational transitions in perfect harmonic relation J'=2/3/5 have been identified in the prestellar core B1b. The molecule has a 1Sigma electronic ground state and its J=2-1 transition presents the hyperfine structure characteristic of a molecule containing a nucleus with spin 1. A careful analysis of possible carriers shows that the best candidate is NS+. The derived rotational constant agrees within 0.3-0.7 % with ab initio calculations. NS+ was also produced in the laboratory to unambiguously validate the astrophysical assignment. The observed rotational frequencies and determined molecular constants confirm the discovery of the nitrogen sulfide cation in space. The chemistry of NS+ and related nitrogen-bearing species has been analyzed by means of a time-dependent gas phase model. The model reproduces well the observed NS/NS+ abundance ratio, in the range 30-50, and indicates that NS+ is formed by reactions of the neutral atoms N and S with the cations SH+ and NH+, respectively., Comment: Accepted for publication in ApJ Letters

We present here a systematic search for cyanopolyynes in the shock region L1157-B1 and its associated protostar L1157-mm in the framework of the Large Program "Astrochemical Surveys At IRAM" (ASAI), dedicated to chemical surveys of solar-type star forming regions with the IRAM 30m telescope. Observations of the millimeter windows between 72 and 272 GHz permitted the detection of HC$_3$N and its $^{13}$C isotopologues, and HC$_5$N (for the first time in a protostellar shock region). In the shock, analysis of the line profiles shows that the emission arises from the outflow cavities associated with L1157-B1 and L1157-B2. Molecular abundances and excitation conditions were obtained from analysis of the Spectral Line Energy Distributions under the assumption of Local Thermodynamical Equilibrium or using a radiative transfer code in the Large Velocity Gradient approximation. Towards L1157mm, the HC$_3$N emission arises from the cold envelope ($T_{rot}=10$ K) and a higher-excitation region ($T_{rot}$= $31$ K) of smaller extent around the protostar. We did not find any evidence of $^{13}$C or D fractionation enrichment towards L1157-B1. We obtain a relative abundance ratio HC$_3$N/HC$_5$N of 3.3 in the shocked gas. We find an increase by a factor of 30 of the HC$_3$N abundance between the envelope of L1157-mm and the shock region itself. Altogether, these results are consistent with a scenario in which the bulk of HC$_3$N was produced by means of gas phase reactions in the passage of the shock. This scenario is supported by the predictions of a parametric shock code coupled with the chemical model UCL_CHEM., Comment: 13 pages, 12 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

The "Square Kilometre Array" (SKA) is a large international radio telescope project characterised, as suggested by its name, by a total collecting area of approximately one square kilometre, and consisting of several interferometric arrays to observe at metric and centimetric wavelengths. The deployment of the SKA will take place in two sites, in South Africa and Australia, and in two successive phases. From its Phase 1, the SKA will be one of the most formidable scientific machines ever deployed by mankind, and by far the most impressive in terms of data throughput and required computing power. With the participation of almost 200 authors from forty research institutes and six private companies, the publication of this French SKA white paper illustrates the strong involvement in the SKA project of the French astronomical community and of a rapidly growing number of major scientific and technological players in the fields of Big Data, high performance computing, energy production and storage, as well as system integration., Comment: Editor in chief: C. Ferrari; Editors: M. Alves, S. Bosse, S. Corbel, A. Ferrari, K. Ferri\`ere, S. Gauffre, E. Josselin, G. Lagache, S. Lambert, G. Marquette, J.-M. Martin, M.-A. Miville-Desch\^enes, L. Montier, B. Semelin, G. Theureau, S. Vergani, N. Vilmer, P. Zarka; Original file with high resolution figures at SKA-France link: https://ska-france.oca.eu/images/SKA-France-Media/FWB_051017.pdf

The high-sensitivity of the IRAM 30-m ASAI unbiased spectral survey in the mm-window allows us to detect NO emission towards both the Class I object SVS13-A and the protostellar outflow shock L1157-B1. We detect the hyperfine components of the $^2\Pi_{\rm 1/2}$ $J$ = 3/2 $\to$ 1/2 (at 151 GHz) and the $^2\Pi_{\rm 1/2}$ $J$ = 5/2 $\to$ 3/2 (250 GHz) spectral pattern. The two objects show different NO profiles: (i) SVS13-A emits through narrow (1.5 km s$^{-1}$) lines at the systemic velocity, while (ii) L1157-B1 shows broad ($\sim$ 5 km s$^{-1}$) blue-shifted emission. For SVS13-A the analysis leads to $T_{\rm ex}$ $\geq$ 4 K, $N(\rm NO)$ $\leq$ 3 $\times$ 10$^{15}$ cm$^{-2}$, and indicates the association of NO with the protostellar envelope. In L1157-B1, NO is tracing the extended outflow cavity: $T_{\rm ex}$ $\simeq$ 4--5 K, and $N(\rm NO)$ = 5.5$\pm$1.5 $\times$ 10$^{15}$ cm$^{-2}$. Using C$^{18}$O, $^{13}$C$^{18}$O, C$^{17}$O, and $^{13}$C$^{17}$O ASAI observations we derive an NO fractional abundance less than $\sim$ 10$^{-7}$ for the SVS13-A envelope, in agreement with previous measurements towards extended PDRs and prestellar objects. Conversely, a definite $X(NO)$ enhancement is measured towards L1157-B1, $\sim$ 6 $\times$ 10$^{-6}$, showing that the NO production increases in shocks. The public code UCLCHEM was used to interpret the NO observations, confirming that the abundance observed in SVS13-A can be attained in an envelope with a gas density of 10$^5$ cm$^{-3}$ and a kinetic temperature of 40 K. The NO abundance in L1157-B1 is reproduced with pre-shock densities of 10$^5$ cm$^{-3}$ subjected to a $\sim$ 45 km s$^{-1}$ shock., Comment: MNRAS, in press

HCCNC and HNC3 are less commonly found isomers of cyanoacetylene, HC3N, a molecule that is widely found in diverse astronomical sources. We want to know if HNC3 is present in sources other than the dark cloud TMC-1 and how its abundance is relative to that of related molecules. We used the ASAI unbiased spectral survey at IRAM 30m towards the prototypical prestellar core L1544 to search for HNC3 and HCCNC which are by-product of the HC3NH+ recombination, previously detected in this source. We performed a combined analysis of published HNC3 microwave rest frequencies with thus far unpublished millimeter data because of issues with available rest frequency predictions. We determined new spectroscopic parameters for HNC3, produced new predictions and detected it towards L1544. We used a gas-grain chemical modelling to predict the abundances of N-species and compare with the observations. The modelled abundances are consistent with the observations, considering a late stage of the evolution of the prestellar core. However the calculated abundance of HNC3 was found 5-10 times higher than the observed one. The HC3N, HNC3 and HCCNC versus HC3NH+ ratios are compared in the TMC-1 dark cloud and the L1544 prestellar core., Comment: Accepted in MNRAS letters. 5 pages plus 2 additional pages for the on-line material

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

Context: Methanol is thought to be mainly formed during the prestellar phase and its deuterated form keeps memory of the conditions at that epoch. Thanks to the unique combination of high angular resolution and sensitivity provided by ALMA, we wish to measure methanol deuteration in the planet formation region around a Class 0 protostar and to understand its origin. Aims: We mapped both the $^{13}$CH$_3$OH and CH$_2$DOH distribution in the inner regions ($\sim$100 au) of the HH212 system in Orion B. To this end, we used ALMA Cycle 1 and Cycle 4 observations in Band 7 with angular resolution down to $\sim$0.15$"$. Results: We detected 6 lines of $^{13}$CH$_3$OH and 13 lines of CH$_2$DOH with upper level energies up to 438 K in temperature units. We derived a rotational temperature of (171 $\pm$ 52) K and column densities of 7$\times$10$^{16}$ cm$^{-2}$ ($^{13}$CH$_3$OH) and 1$\times$10$^{17}$ cm$^{-2}$ (CH$_2$DOH), respectively. Consequently, the D/H ratio is (2.4 $\pm$ 0.4)$\times$10$^{-2}$, a value lower by an order of magnitude with respect to what was previously measured using single dish telescopes toward protostars located in Perseus. Our findings are consistent with the higher dust temperatures in Orion B with respect to that derived for the Perseus cloud. The emission is tracing a rotating structure extending up to 45 au from the jet axis and elongated by 90 au along the jet axis. So far, the origin of the observed emission appears to be related with the accretion disk. Only higher spatial resolution measurements however, will be able to disentangle between different possible scenarios: disk wind, disk atmosphere, or accretion shocks., Comment: Astronomy & Astrophysics Letter, in press

The early stages of low-mass star formation are likely to be subject to intense ionization by protostellar energetic MeV particles. As a result, the surrounding gas is enriched in molecular ions, such as HCO$^{+}$ and N$_{2}$H$^{+}$. Nonetheless, this phenomenon remains poorly understood for Class 0 objects. Recently, based on Herschel observations taken as part of the key program Chemical HErschel Surveys of Star forming regions (CHESS), a very low HCO$^{+}$/N$_{2}$H$^{+}$ abundance ratio of about 3-4, has been reported toward the protocluster OMC-2 FIR4. This finding suggests a cosmic-ray ionization rate in excess of 10$^{-14}$ s$^{-1}$, much higher than the canonical value of $\zeta$ = 3$\times$10$^{-17}$ s$^{-1}$ (value expected in quiescent dense clouds). To assess the specificity of OMC-2 FIR4, we have extended this study to a sample of sources in low- and intermediate mass. More specifically, we seek to measure the HCO$^{+}$/N$_2$H$^{+}$ abundance ratio from high energy lines (J $\ge$ 6) toward this source sample in order to infer the flux of energetic particles in the warm and dense gas surrounding the protostars. We use observations performed with the Heterodyne Instrument for the FarInfrared spectrometer on board the Herschel Space Observatory toward a sample of 9 protostars. We report HCO$^{+}$/N$_2$H$^{+}$ abundance ratios in the range of 5 up to 73 toward our source sample. The large error bars do not allow us to conclude whether OMC-2~FIR4 is a peculiar source. Nonetheless, an important result is that the measured HCO$^{+}$/N$_2$H$^{+}$ ratio does not vary with the source luminosity. At the present time, OMC-2 FIR4 remains the only source where a high flux of energetic particles is clearly evident. More sensitive and higher angular resolution observations are required to further investigate this process., Comment: Accepted for publication in A&A 18 pages, 12 figures, 15 tables

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

We used the Atacama Large (sub-)Millimeter Array (ALMA) and the IRAM Plateau de Bure Interferometer (PdBI) to image, with an angular resolution of 0.5$''$ (120 au) and 1$''$ (235 au), respectively, the emission from 11 different organic molecules in the protostellar binary NGC1333 IRAS 4A. We clearly disentangled A1 and A2, the two protostellar cores present. For the first time, we were able to derive the column densities and fractional abundances simultaneously for the two objects, allowing us to analyse the chemical differences between them. Molecular emission from organic molecules is concentrated exclusively in A2 even though A1 is the strongest continuum emitter. The protostellar core A2 displays typical hot corino abundances and its deconvolved size is 70 au. In contrast, the upper limits we placed on molecular abundances for A1 are extremely low, lying about one order of magnitude below prestellar values. The difference in the amount of organic molecules present in A1 and A2 ranges between one and two orders of magnitude. Our results suggest that the optical depth of dust emission at these wavelengths is unlikely to be sufficiently high to completely hide a hot corino in A1 similar in size to that in A2. Thus, the significant contrast in molecular richness found between the two sources is most probably real. We estimate that the size of a hypothetical hot corino in A1 should be less than 12 au. Our results favour a scenario in which the protostar in A2 is either more massive and/or subject to a higher accretion rate than A1, as a result of inhomogeneous fragmentation of the parental molecular clump. This naturally explains the smaller current envelope mass in A2 with respect to A1 along with its molecular richness., Comment: Accepted in Astronomy and Astrophysics

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

L1157-B1 is the brightest shocked region of the large-scale molecular outflow, considered the prototype of chemically rich outflows, being the ideal laboratory to study how shocks affect the molecular gas. Several deuterated molecules have been previously detected with the IRAM 30m, most of them formed on grain mantles and then released into the gas phase due to the shock. We aim to observationally investigate the role of the different chemical processes at work that lead to formation the of DCN and test the predictions of the chemical models for its formation. We performed high-angular resolution observations with NOEMA of the DCN(2-1) and H13CN(2-1) lines to compute the deuterated fraction, Dfrac(HCN). We detected emission of DCN(2-1) and H13CN(2-1) arising from L1157-B1 shock. Dfrac(HCN) is ~4x10$^{-3}$ and given the uncertainties, we did not find significant variations across the bow-shock. Contrary to HDCO, whose emission delineates the region of impact between the jet and the ambient material, DCN is more widespread and not limited to the impact region. This is consistent with the idea that gas-phase chemistry is playing a major role in the deuteration of HCN in the head of the bow-shock, where HDCO is undetected as it is a product of grain-surface chemistry. The spectra of DCN and H13CN match the spectral signature of the outflow cavity walls, suggesting that their emission result from shocked gas. The analysis of the time dependent gas-grain chemical model UCL-CHEM coupled with a C-type shock model shows that the observed Dfrac(HCN) is reached during the post-shock phase, matching the dynamical timescale of the shock. Our results indicate that the presence of DCN in L1157-B1 is a combination of gas-phase chemistry that produces the widespread DCN emission, dominating in the head of the bow-shock, and sputtering from grain mantles toward the jet impact region., Comment: Accepted for publication in A&A. 7 pages, 5 Figures, 1 Table

The shock L1157-B1 driven by the low-mass protostar L1157-mm is an unique environment to investigate the chemical enrichment due to molecules released from dust grains. IRAM-30m and Plateau de Bure Interferometer observations allow a census of Si-bearing molecules in L1157-B1. We detect SiO and its isotopologues and, for the first time in a shock, SiS. The strong gradient of the [SiO/SiS] abundance ratio across the shock (from >=180 to ~25) points to a different chemical origin of the two species. SiO peaks where the jet impacts the cavity walls ([SiO/H2] ~ 1e-6), indicating that SiO is directly released from grains or rapidly formed from released Si in the strong shock occurring at this location. In contrast, SiS is only detected at the head of the cavity opened by previous ejection events ([SiS/H2] ~ 2e-8). This suggests that SiS is not directly released from the grain cores but instead should be formed through slow gas-phase processes using part of the released silicon. This finding shows that Si-bearing molecules can be useful to distinguish regions where grains or gas-phase chemistry dominates., Comment: 5 pages, 4 figure, 1 page of appendices

We report on a systematic search for oxygen-bearing Complex Organic Molecules (COMs) in the Solar-like protostellar shock region L1157-B1, as part of the IRAM Large Program "Astrochemical Surveys At IRAM" (ASAI). Several COMs are unambiguously detected, some for the first time, such as ketene H$_2$CCO, dimethyl ether (CH$_3$OCH$_3$) and glycolaldehyde (HCOCH$_2$OH), and others firmly confirmed, such as formic acid (HCOOH) and ethanol (C$_2$H$_5$OH). Thanks to the high sensitivity of the observations and full coverage of the 1, 2 and 3mm wavelength bands, we detected numerous (10--125) lines from each of the detected species. Based on a simple rotational diagram analysis, we derive the excitation conditions and the column densities of the detected COMs. Combining our new results with those previously obtained towards other protostellar objects, we found a good correlation between ethanol, methanol and glycolaldehyde. We discuss the implications of these results on the possible formation routes of ethanol and glycolaldehyde., Comment: 5 pages; 1 table; 3 figures; to appear in MNRAS

Protostellar jets and outflows are key features of the star-formation process, and primary processes of the feedback of young stars on the interstellar medium. Understanding the underlying shocks is necessary to explain how jets and outflows are launched, and to quantify their chemical and energetic impacts on the surrounding medium. We performed a high-spectral resolution study of the [OI]$_{\rm 63 \mu m}$ emission in the outflow of the intermediate-mass Class 0 protostar Cep E-mm. We present observations of the OI $^3$P$_1 \rightarrow$ $^3$P$_2$, OH between $^2\Pi_{1/2}$ $J = 3/2$ and $J = 1/2$ at 1837.8 GHz, and CO (16-15) lines with SOFIA-GREAT at three positions in the Cep E outflow: mm (the driving protostar), BI (in the southern lobe), and BII (the terminal position in the southern lobe). The CO line is detected at all three positions. The OI line is detected in BI and BII, whereas the OH line is not detected. In BII, we identify three kinematical components in OI and CO, already detected in CO: the jet, the HH377 terminal bow-shock, and the outflow cavity. The OI column density is higher in the outflow cavity than in the jet, which itself is higher than in the terminal shock. The terminal shock is where the abundance ratio of OI to CO is the lowest (about 0.2), whereas the jet component is atomic (ratio $\sim$2.7). In the jet, we compare the OI observations with shock models that successfully fit the integrated intensity of 10 CO lines: these models do not fit the OI data. The high intensity of OI emission points towards the propagation of additional dissociative or alternative FUV-irradiated shocks, where the illumination comes from the shock itself. From the sample of low-to-high mass protostellar outflows where similar observations have been performed, the effects of illumination seem to increase with the mass of the protostar., Comment: 9 pages, 4 figures, 6 tables, Astronomy and astrophysics, in press