Search

There is evidence that the star formation process is linked to the intricate net of filaments in molecular clouds, which may be also due to gas compression from external triggers. We studied the southern region of the Perseus NGC 1333 molecular cloud, known to be heavily shaped by similar external triggers, to shed light on the process that perturbed the filament where the Class 0 IRAS4 protostars lie. We use new IRAM-NOEMA observations of SiO and CH3OH, both known to trace violent events as shocks, toward IRAS 4A as part of the Large Program Seeds Of Life in Space (SOLIS). We detected three parallel elongated ($>$6000 au) structures, called fingers, with narrow line profiles (~1.5 $km s^{-1}$) peaked at the cloud systemic velocity, tracing gas with high density (5-20 $10^5 cm^{-3}$) and high temperature (80-160 K). They are chemically different, with the northern finger traced by both SiO and CH3OH ([CH3OH]/[SiO]~160-300), while the other two only by SiO ([CH3OH]/[SiO]$$5000 yr, would be consistent with the observations if a substantial fraction of silicon, frozen onto the grain mantles, is released by the shocks.We suggest that the shock train is due to an expanding gas bubble, coming behind NGC 1333 from the southwest and clashing against the filament, where IRAS 4A lies. Finally, we propose a solution to the two-decades long debate on the nature and origin of the widespread narrow SiO emission observed in the south part of NGC 1333, namely that it is due to unresolved trains of shocks., Comment: accepted to MNRAS, stac083

Innovative technological approaches are crucial to enhance naturalistic museum collections and develop information repositories of relevant interest to science, such as threatened animal taxa. In this context, museomics is an emerging discipline that provides a novel approach to the enhancement and exploitation of these collections. In the present study, the discovery of a neglected lion skeleton in the Museum of Zoology “Pietro Doderlein” of the University of Palermo (Italy) offered the opportunity to undertake a multidisciplinary project. The aims of the study consisted of the following: (i) adding useful information for museographic strategies, (ii) obtaining a new genetic data repository from a vulnerable species, (iii) strengthening public awareness of wildlife conservation, and (iv) sharing new learning material. The remains of the lion were examined with a preliminary osteological survey, then they were restored by means of 3D printing of missing skeletal fragments. Phylogenetic analyses based on cytochrome b sequence clearly indicate that the specimen belongs to the Central Africa mitochondrial clade. At the end of the study, the complete and restored skeleton was exhibited, along with all of the information and data available from this project. This study shows a useful approach for the restoration and enhancement of a museum specimen, with important opportunities for preserving biodiversity and driving specific conservation policies, but also for providing Life Science learning material.

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 (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., 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

Imai et al., The chemical diversity of low-mass protostellar sources has so far been recognized, and environmental effects are invoked as its origin. In this context, observations of isolated protostellar sources without the influence of nearby objects are of particular importance. Here, we report the chemical and physical structures of the low-mass Class 0 protostellar source IRAS 16544−1604 in the Bok globule CB 68, based on 1.3 mm Atacama Large Millimeter/submillimeter Array observations at a spatial resolution of ∼70 au that were conducted as part of the large program FAUST. Three interstellar saturated complex organic molecules (iCOMs), CH3OH, HCOOCH3, and CH3OCH3, are detected toward the protostar. The rotation temperature and the emitting region size for CH3OH are derived to be 131 ± 11 K and ∼10 au, respectively. The detection of iCOMs in close proximity to the protostar indicates that CB 68 harbors a hot corino. The kinematic structure of the C18O, CH3OH, and OCS lines is explained by an infalling–rotating envelope model, and the protostellar mass and the radius of the centrifugal barrier are estimated to be 0.08–0.30 M⊙ and, This paper makes use of the following ALMA data: ADS/JAO.ALMA#2016.1.01376.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), NSC and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The authors acknowledge financial support from JSPS and MAEE under the Japan–France Integrated Action Program (SAKURA: 25765VC). D.J. is supported by NRC Canada and by an NSERC Discovery Grant. E.B. and G.B. acknowledge funding from the European Research Council under the European Union's Horizon 2020 research and innovation program, for the project "The Dawn of Organic Chemistry," grant agreement No. 741002. G.B. acknowledges support from the PID2020-117710GB-I00 grant funded by MCIN/AEI/10.13039/501100011033 and from the Unit of Excellence María de Maeztu 2020–2023 award to the Institute of Cosmos Sciences (CEX2019-00918-M). I.J.-S. has received partial support from the Spanish State Research Agency (AEI) through project Nos. PID2019-105552RB-C41 and MDM-2017-0737 Unidad de Excelencia "Maria de Maeztu" - Centro de Astrobiologia (CSIC-INTA). This study is supported by grants-in-aid from the Ministry of Education, Culture, Sports, Science, and Technology of Japan (18H05222, 18J11010, 19H05069, 19K14753, and 21K13954).

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

We have observed the very low-mass Class 0 protostar IRAS 15398−3359 at scales ranging from 50 to 1800 au, as part of the Atacama Large Millimeter/Submillimeter Array Large Program FAUST. We uncover a linear feature, visible in H2CO, SO, and C18O line emission, which extends from the source in a direction almost perpendicular to the known active outflow. Molecular line emission from H2CO, SO, SiO, and CH3OH further reveals an arc-like structure connected to the outer end of the linear feature and separated from the protostar, IRAS 15398−3359, by 1200 au. The arc-like structure is blueshifted with respect to the systemic velocity. A velocity gradient of 1.2 km s−1 over 1200 au along the linear feature seen in the H2CO emission connects the protostar and the arc-like structure kinematically. SO, SiO, and CH3OH are known to trace shocks, and we interpret the arc-like structure as a relic shock region produced by an outflow previously launched by IRAS 15398−3359. The velocity gradient along the linear structure can be explained as relic outflow motion. The origins of the newly observed arc-like structure and extended linear feature are discussed in relation to turbulent motions within the protostellar core and episodic accretion events during the earliest stage of protostellar evolution. The authors thank the anonymous reviewer for invaluable comments. They are also grateful to Hauyu Baobab Liu for useful discussions. This paper makes use of the following ALMA data set: ADS/JAO.ALMA#2018.1.01205.L (PI: Satoshi Yamamoto). ALMA is a partnership of the ESO (representing its member states), the NSF (USA) and NINS (Japan), together with the NRC (Canada) and the NSC and ASIAA (Taiwan), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by the ESO, the AUI/NRAO, and the NAOJ. The authors thank to the ALMA staff for their excellent support. This work is supported by the projects PRIN-INAF-MAIN-STREAM 2017 Protoplanetary disks seen through the eyes of new-generation instruments, ARIEL and the astrochemical link between circumstellar disks and planets, and the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programmes: The Dawn of Organic Chemistry (DOC), grant agreement No. 741002, and Astro-Chemistry Origins (ACO), grant No. 811312. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. I.J.-S. has received partial support from the Spanish FEDER (project number ESP2017-86582-C4-1-R) and the State Research Agency (AEI; project number PID2019-105552RB-C41). D.J. is supported by NRC Canada and by an NSERC Discovery Grant. This project is also supported by a Grant-in-Aid from Japan Society for the Promotion of Science (KAKENHI: Nos. 18H05222, 19H05069, 19K14753). Yuki Okoda thanks the Advanced Leading Graduate Course for Photon Science (ALPS) and Japan Society for the Promotion of Science (JSPS) for financial support. Peerreview

We study star formation in the Center Ridge 1 (CR1) clump in the Vela C giant molecular cloud, selected as a high column density region that shows the lowest level of dust continuum polarization angle dispersion, likely indicating that the magnetic field is relatively strong. We observe the source with the ALMA 7m-array at 1.05~mm and 1.3~mm wavelengths, which enable measurements of dust temperature, core mass and astrochemical deuteration. A relatively modest number of eleven dense cores are identified via their dust continuum emission, with masses spanning from 0.17 to 6.7 Msun. Overall CR1 has a relatively low compact dense gas fraction compared with other typical clouds with similar column densities, which may be a result of the strong magnetic field and/or the very early evolutionary stage of this region. The deuteration ratios, Dfrac, of the cores, measured with N2H+(3-2) and N2D+(3-2) lines, span from 0.011 to 0.85, with the latter being one of the highest values yet detected. The level of deuteration appears to decrease with evolution from prestellar to protostellar phase. A linear filament, running approximately parallel with the large scale magnetic field orientation, is seen connecting the two most massive cores, each having CO bipolar outflows aligned orthogonally to the filament. The filament contains the most deuterated core, likely to be prestellar and located midway between the protostars. The observations permit measurement of the full deuteration structure of the filament along its length, which we present. We also discuss the kinematics and dynamics of this structure, as well as of the dense core population., 25 pages, 17 figures, accepted by ApJ

Roccapelago (MO) is a small village located in the Northern Central Apennines, with a population of 31 inhabitants (2014). In 2010, more than 400 individuals dated between the end of the 16th and the 18th century, many of which partially mummified, were discovered in the crypt of the church. This small village, because of its geographical location and surrounding environment, seems to possess the characteristics of a genetic isolate, useful for population genetics and genealogical analyses. Thus, a diachronic study of DNA aimed at investigating the structure and dynamics of the population of Roccapelago over the about 4 centuries, was conducted by analyzing ancient and modern inhabitants of the village. The 14 modern samples were selected by considering both the founder surnames of the village, identified thanks to the study of parish registers, and the grandparent’s criterion. From 25 ancient mummies, morphologically assigned to male individuals, the petrous bone, that harbors high DNA amounts, was selected for the DNA extraction. The quantification and qualitative assessment of total human male DNA were evaluated by a real-time PCR assay using the Quantifiler Trio DNA Quantification Kit and multiplex PCR of 27 Y-chromosome short tandem repeat (Y-STR) markers included in the Yfiler Plus PCR Amplification Kit, with seven rapidly mutating Y-STR loci for improving discrimination of male lineages, was performed to genotype the samples. Y-STRs were analyzed according to the criteria of ancient DNA (aDNA) analysis to ensure that authentic DNA typing results were obtained from these ancient samples. The molecular analysis showed the usefulness of the Y chromosome to identify historically relevant remains and discover patterns of relatedness in communities moving from anthropology to genetic genealogy and forensics.

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., 7 pages, 6 figures

Background: Phoenician and Punic expansions have been protagonists of intense trade networks and settlements in the Mediterranean Sea. Aims: The maternal genetic variability of ancient Punic samples from the Sardinian necropolis of Tharros was analysed, with the aim to explore genetic interactions and signatures of past population events. Subjects and methods: The mtDNA HVS-I and coding region SNPs were analysed in 14 Punic samples and 74 modern individuals from Cabras and Belvì (for which the HVS-II region was also analysed). The results were compared with 5,590 modern Euro-Mediterranean sequences and 127 ancient samples. Results: While contemporary groups fall within the genetic variability of other modern Sardinians, our Punic samples reveal proximity to present-day North-African and Iberian populations. Furthermore, Cabras and Belvì cluster mainly with pre-Phoenician groups, while samples from Tharros project with other Punic Sardinian individuals. Conclusion: This study provides the first preliminary insights into the population dynamics of the Punic site of Tharros. While the number of currently available samples does not allow definitive investigation of the connection with indigenous Sardinian groups, our results seem to confirm internal migratory phenomena in the central-western Mediterranean and female participation in the Punic mobility.

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, Comment: 14 pages, 11 figures, accepted by Astronomy and Astrophysics

The Italian peninsula was host to a strong history of migration processes that shaped its genomic variability since prehistoric times. During the Metal Age, Sicily and Southern Italy were the protagonists of intense trade networks and settlements along the Mediterranean. Nonetheless, ancient DNA studies in Southern Italy are, at present, still limited to prehistoric and Roman Apulia. Here, we present the first mitogenomes from a Middle Bronze Age cave burial in Calabria to address this knowledge gap. We adopted a hybridization capture approach, which enabled the recovery of one complete and one partial mitochondrial genome. Phylogenetic analysis assigned these two individuals to the H1e and H5 subhaplogroups, respectively. This preliminary phylogenetic analysis supports affinities with coeval Sicilian populations, along with Linearbandkeramik and Bell Beaker cultures maternal lineages from Central Europe and Iberia. Our work represents a starting point which contributes to the comprehension of migrations and population dynamics in Southern Italy, and highlights this knowledge gap yet to be filled by genomic studies.

Context. Recent results in astrochemistry have revealed that some molecules, such as interstellar complex organic species and deuterated species, can serve as valuable tools in the investigation of star-forming regions. Sulphuretted species can also be used to follow the chemical evolution of the early stages of a Sun-like star formation process. 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 that is rich in interstellar complex organic molecules. Methods. To this end, we used the NGC 1333 SVS13-A data at 3 mm and 1.4 mm obtained with the IRAM-NOEMA interferometer in the framework of the SOLIS (Seeds of Life in Space) Large Program. The line emission of S-bearing species was imaged and analyzed using local thermodynamic equilibrium (LTE) and large velocity gradient (LVG) approaches. Results. We imaged the spatial distribution on ≤300 au scale of the line emission of 32SO, 34SO, C32S, C34S, C33S, OCS, H2C32S, H2C34S, and NS. The low excitation (9 K) 32SO line traces: (i) the low-velocity SVS13-A outflow and (ii) the fast (up to 100 km s−1 away from the systemic velocity) collimated jet driven by the nearby SVS13-B Class 0 object. Conversely, the rest of the lines are confined in the inner SVS13-A region, where complex organics were previously imaged. More specifically, the non-LTE LVG analysis of SO, SO2, and H2CS indicates a hot corino origin (size in the 60–120 au range). Temperatures between 50 K and 300 K, as well as volume densities larger than 105 cm−3 have been derived. The abundances of the sulphuretted 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 the assessment that the NS+ ion is mainly formed in the extended envelope. Conclusions. The [H2CS]/[H2CO] ratio, once measured at high-spatial resolutions, increases with time (from Class 0 to Class II objects) by more than one order of magnitude (from ≤10−2 to a few 10−1). This suggests that [S]/[O] changes along the process of Sun-like star formation. Finally, the estimate of the [S]/[H] budget in SVS13-A is 2–17% of the Solar System value (1.8 × 10−5), which is consistent with what was previously measured towards Class 0 objects (1–8%). This finding supports the notion that the enrichment of the sulphuretted species with respect to dark clouds remains 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 the scale of the Solar System.

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

Bianchi, E. et al., 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 (Atacama Large Millimeter/submillimeter Array) Large Program FAUST (Fifty AU STudy of the chemistry in the disc/envelope system of solar-like protostars).We detected several lines from methanol and its isotopologues (13CH3OH and CH2DOH), methyl formate, and ethanol. Lines are bright towards 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 (100 K), density (1.5 × 108 cm3), and emitting size (10 au in radius). All CH3OH and 13CH3OH 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 this work, little chemical evolution from Class 0 to I hot corinos occurs., With funding from the Spanish government through the "María de Maeztu Unit of Excellence" accreditation (MDM-2017-0737)

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., 27 pages, 14 figures, 20 tables - Accepted for publication in A&A on 25 June 2020

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

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

In the interstellar medium carbon exists in the form of two stable isotopes $^{12}$C and $^{13}$C and their ratio is a good indicator of nucleosynthesis in galaxies. However, chemical fractionation can potentially significantly alter this ratio and in fact observations of carbon fractionation within the same galaxy has been found to vary from species to species. In this paper we theoretically investigate the carbon fractionation for selected abundant carbon-bearing species in order to determine the conditions that may lead to a spread of the $^{12}$C/$^{13}$C ratio in external galaxies. We find that carbon fractionation is sensitive to almost all the physical conditions we investigated, it strongly varies with time for all species but CO, and shows pronounced differences across species. Finally we discuss our theoretical results in the context of the few observations of the $^{12}$C/$^{13}$C ratio in both local and higher redshift galaxies., MNRAS, in press

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 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., 10 pages, 5 figures, accepted for publication by MNRAS

Aims. We present the Forgotten Quadrant Survey (FQS), an ESO large project that used the 12 m antenna of the Arizona Radio Observatory to map the Galactic plane in the range 220° < l < 240° and −2.°5 < b < 0°, both in 12CO (1–0) and 13CO (1–0), at a spectral resolution of 0.65 km s−1 and 0.26 km s−1. Methods. We used the (1–0) transition of carbon monoxide to trace the molecular component of the interstellar medium. Our data set allows us to easily identify how the molecular dense gas is organised at different spatial scales: from the giant clouds with their denser filamentary networks, down to the clumps and cores that host the new-born stars and to obtain reliable estimates of their key physical parameters such as size and mass. Results. We present the first release of the data of the FQS survey and discuss their quality. Spectra with 0.65 km s−1 velocity channels have noise ranging from 0.8 K to 1.3 K for 12CO (1–0) and from 0.3 K to 0.6 K for 13CO (1–0). In this first paper, we used the 12CO (1–0) spectral cubes to produce a catalogue of 263 molecular clouds. The clouds are grouped in three main structures corresponding to the Local, Perseus, and Outer arms up to a distance of ∼8.6 kpc from the Sun. This is the first self-consistent statistical catalogue of molecular clouds of the outer Galaxy obtained with a subarcminute spatial resolution and therefore able to detect not only the classical giant molecular clouds, but also the small clouds and to resolve the cloud structure at the sub-parsec scale up to a distance of a few kiloparsec. We found two classes of objects: structures with sizes above a few parsecs that are typical molecular clouds and may be self-gravitating, and subparsec structures that cannot be in gravitational equilibrium and are likely transient or confined by external pressure. We used the ratio between the Herschel H2 column density and the integrated intensity of the CO lines to calculate the CO conversion factor and we found mean values of (3.3 ± 1.4) × 1020 cm−2 (K km s−1)−1 and (1.2 ± 0.4) × 1021 cm−2 (K km s−1)−1, for 12CO (1–0) and 13CO (1–0), respectively. Conclusions. The FQS contributes to the general effort in producing a new generation of high-quality spectroscopic data for the Galactic plane in the less-studied third Galactic quadrant toward the outer Galaxy. The FQS has produced a data set of great legacy value, largely improving the data quality both in terms of sensitivity and spatial resolution over previous data sets.

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

Context. Phosphorus-bearing species are essential to the formation of life on Earth, however they have barely been detected in the interstellar medium. In particular, towards star-forming regions only PN and PO have been identified so far. Since only a small number of detections of P-bearing molecules are available, their chemical formation pathways are not easy to constrain and are thus highly debatable. An important factor still missing in the chemical models is the initial elemental abundance of phosphorus, that is, the depletion level of P at the start of chemical models of dense clouds. Aims. In order to overcome this problem, we study P-bearing species in diffuse and translucent clouds. In these objects phosphorus is expected to be mainly in the gas phase and therefore the elemental initial abundance needed in our chemical simulations corresponds to the cosmic one and is well constrained. Methods. For the study of P-bearing chemistry we used an advanced chemical model. We updated and significantly extended the P-chemistry network based on chemical databases and previous literature. We performed single-pointing observations with the IRAM 30 m telescope in the 3 mm range towards the line of sight to the strong continuum source B0355+508 aiming for the (2-1) transitions of PN, PO, HCP, and CP. This line of sight incorporates five diffuse and/or translucent clouds. Results. The (2-1) transitions of the PN, PO, HCP, and CP were not detected. We report high signal-to-noise-ratio detections of the (1-0) lines of 13CO, HNC, and CN along with a first detection of C34S towards this line of sight. We have attempted to reproduce the observations of HNC, CN, CS, and CO in every cloud with our model by applying typical physical conditions for diffuse or translucent clouds. We find that towards the densest clouds with vLSR = -10, - 17 km s-1 the best-fit model is given by the parameters (n(H), AV, Tgas) = (300 cm-3, 3 mag, 40 K). Conclusions. According to our best-fit model, the most abundant P-bearing species are HCP and CP (~10-10). The molecules PN, PO, and PH3 also show relatively high predicted abundances of ~10-11. We show that the abundances of these species are sensitive to visual extinction, cosmic-ray ionization rate, and the diffusion-to-desorption energy ratio on dust grains. The production of P-bearing species is favored towards translucent rather than diffuse clouds, where the environment provides a stronger shielding from the interstellar radiation. Based on our improved model, we show that the (1-0) transitions of HCP, CP, PN, and PO are expected to be detectable with estimated intensities of up to ~200 mK. © J. Chantzos et al. 2020. Horizon 2020 Framework Programme, H2020: 66493 lzp-2018/1-0170 We thank the anonymous referee for his/her comments that significantly improved the present manuscript. The authors also wish to thank the IRAM Granada staff for their help during the observations. V.M.R. has received funding from the European Unions Horizon 2020 research and innovation programme under the Marie Skodowska-Curie grant agreement No 66493. Work by A.V. is supported by the Latvian Science Council via the project lzp-2018/1-0170. J.C. acknowledges Dr. J. C. Laas for his support with the Python programming.

Context. One of the goals of astrochemistry is to understand the degree of chemical complexity that can be reached in star-forming regions, along with the identification of precursors of the building blocks of life in the interstellar medium. To answer such questions, unbiased spectral surveys with large bandwidth and high spectral resolution are needed, in particular, to resolve line blending in chemically rich sources and identify each molecule (especially for complex organic molecules). These kinds of observations have already been successfully carried out, primarily towards the Galactic Center, a region that shows peculiar environmental conditions. Aims. We present an unbiased spectral survey 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. The aim of this 3mm spectral survey is to identify and characterize the physical parameters of the gas emission in different molecular species, focusing on complex organic molecules. In this first paper, we present the survey and discuss the detection and relative abundances of the three isomers of C2H4O2: methyl formate, glycolaldehyde, and acetic acid. Methods. Observations were carried out with the ALMA interferometer, covering all of band 3 from 84 to 116 GHz (∼32 GHz bandwidth) with an angular resolution of 1.2″ × 1.2″ (∼ 4400 au × 4400 au) and a spectral resolution of ∼0.488 MHz (∼1.3-1.7 km s-1). The transitions of the three molecules have been analyzed with the software XCLASS to determine the physical parameters of the emitted gas. Results. All three isomers were detected with abundances of (2 ± 0.6) × 10-7, (4.3-8) × 10-8, and (5.0 ± 1.4) × 10-9 for methyl formate, acetic acid, and glycolaldehyde, respectively. Methyl formate and acetic acid abundances are the highest detected up to now, if compared to sources in the literature. The size of the emission varies among the three isomers with acetic acid showing the most compact emission while methyl formate exhibits the most extended emission. Different chemical pathways, involving both grain-surface chemistry and cold or hot gas-phase reactions, have been proposed for the formation of these molecules, but the small number of detections, especially of acetic acid and glycolaldehyde, have made it very difficult to confirm or discard the predictions of the models. The comparison with chemical models in literature suggests the necessity of grain-surface routes for the formation of methyl formate in G31, while for glycolaldehyde both scenarios could be feasible. The proposed grain-surface reaction for acetic acid is not capable of reproducing the observed abundance in this work, while the gas-phase scenario should be further tested, given the large uncertainties involved., With funding from the Spanish government through the "María de Maeztu Unit of Excellence" accreditation (MDM-2017-0737)

One of the most important tools to investigate the chemical history of our Galaxy and our own Solar System is to measure the isotopic fractionation of chemical elements. This is the process that distributes the less abundant stable isotopes of an element in different molecules. The isotopic ratios are governed by two main processes: 1. chemical evolution of the whole Galaxy due to stellar nucleosynthesis; 2. local fractionation effects. In this Proceeding we report some results highlighting both processes towards massive star-forming regions.

Context. Isotopic fractionation is an important tool for investigating the chemical history of our Solar System. In particular, the isotopic fraction of nitrogen (14N/15N) is lower in comets and other pristine Solar System bodies with respect to the value measured for the protosolar nebula, suggesting a local chemical enrichment of 15N during the formation of the Solar System. Therefore, interferometric studies of nitrogen fractionation in Solar System precursors are needed for us to obtain clues about our astrochemical origins. Aims. In this work we have investigated the variation in the 14N/15N ratio in one of the closest analogues of the environment in which the Solar System was born: the protocluster OMC-2 FIR4. We present the first comparison at high angular resolution between HCN and N2H+ using interferometric data. Methods. We analysed observations of the HCN isotopologues H13CN and HC15N in the OMC-2 FIR4 protocluster. Specifically, we observed the transitions H13CN (1−0) and HC15N (1−0) with the NOrthern Extended Millimeter Array (NOEMA) within the context of the IRAM Seeds Of Life In Space (SOLIS) Large Program. We combined our results with analysis of archival data obtained with the Atacama Large Millimeter Array of N2H+ and its 15N isotopologues. Results. Our results show a small regional variation in the 14N/15N ratio for HCN, from ~250 to 500. The ratios in the central regions of FIR4, where the candidate protostars are located, are largely consistent with one another and within that range (~300). They also show little variation from the part of the protocluster known to harbour a high cosmic-ray ionisation rate to the portion with a lower rate. We also found a small variation in the 14N/15N ratio of N2H+ across different regions, from ~200 to ~400. Conclusions. These results suggest that local changes in the physical parameters occurring on the small linear scales probed by our observations in the protocluster do not seem to affect the 14N/15N ratio in either HCN or N2H+ and hence that this is independent of the molecule used. Moreover, the high level of irradiation due to cosmic rays does not affect the N fractionation either.

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)., 5 pages, 4 figures, 1 table, accepted for publication in MNRAS

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

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

Many pieces of evidence indicate that the Solar System youth was marked by violent processes: among others, high fluxes of energetic particles (greater than or equal to 10 MeV) are unambiguously recorded in meteoritic material, where an overabundance of the short-lived 10 Be products is measured. Several hypotheses have been proposed to explain from where these energetic particles originate, but there is no consensus yet, mostly because of the scarcity of complementary observational constraints. In general, the reconstruction of the past history of the Solar System is best obtained by simultaneously considering what we know of it and of similar systems nowadays in formation. However, when it comes to studying the presence of energetic particles in young forming stars, we encounter the classical problem of the impossibility of directly detecting them toward the emitting source (analogously to what happens to galactic cosmic rays). Yet, exploiting the fact that energetic particles, such as cosmic rays, create H 3 + and that an enhanced abundance of H 3 + causes dramatic changes on the overall gas chemical composition, we can indirectly estimate the flux of energetic particles. This contribution provides an overview of the search for solar-like protostars permeated by energetic particles and the discovery of a protocluster, OMC-2 FIR4, where the phenomenon is presently occurring. This article is part of a discussion meeting issue ‘Advances in hydrogen molecular ions: H 3 + , H 5 + and beyond’.

In star forming regions in our own Galaxy, the 14N/15N ratio is found to vary from $\sim$ 100 in meteorites, comets and protoplanetary disks up to $\sim$ 1000 in pre-stellar and star forming cores, while in external galaxies the very few single-dish large scale measurements of this ratio lead to values of 100-450. The extent of the contribution of isotopic fractionation to these variations is, to date, unknown. In this paper we present a theoretical chemical study of nitrogen fractionation in external galaxies in order to determine the physical conditions that may lead to a spread of the 14N/15N ratio from the solar value of $\sim$440 and hence evaluate the contribution of chemical reactions in the ISM to nitrogen fractionation. We find that the main cause of ISM enrichment of nitrogen fractionation is high gas densities, aided by high fluxes of cosmic rays., Accepted by MNRAS

We present ALMA CO(2-1) observations toward a massive infrared dark cloud G28.37+0.07. The ALMA data reveal numerous molecular (CO) outflows with a wide range of sizes throughout the cloud. Sixty-two 1.3 mm continuum cores were identified to be driving molecular outflows. We have determined the position angle in the plane-of-sky of 120 CO outflow lobes and studied their distribution. We find that the distribution of the plane-of-sky outflow position angles peaks at about 100 degree, corresponding to a concentration of outflows with an approximately east-west direction. For most outflows, we have been able to estimate the plane-of-sky angle between the outflow axis and the filament that harbors the protostar that powers the outflow. Statistical tests strongly indicate that the distribution of outflow-filament orientations is consistent with most outflow axes being mostly orthogonal to their parent filament in 3D. Such alignment may result from filament fragmentation or continuous mass transportation from filament to the embedded protostellar core. The latter is suggested by recent numerical studies with moderately strong magnetic fields., Comment: 4 figures, 1 table, accepted by ApJ

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

Context. C2O and C3O belong to the carbon chain oxides family. Both molecules have been detected in the gas phase towards several star-forming regions, and to explain the observed abundances, ion-molecule gas-phase reactions have been invoked. On the other hand, laboratory experiments have shown that carbon chain oxides are formed after energetic processing of CO-rich solid mixtures. Therefore, it has been proposed that they are formed in the solid phase in dense molecular clouds after cosmic ion irradiation of CO-rich icy grain mantles and released in the gas phase after their desorption. Aims. In this work, we contribute to the understanding of the role of both gas-phase reactions and energetic processing in the formation of simple carbon chain oxides that have been searched for in various low-mass star-forming regions. Methods. We present observations obtained with the Noto-32m and IRAM-30 m telescopes towards star-forming regions. We compare these with the results of a gas-phase model that simulates C2O and C3O formation and destruction, and laboratory experiments in which both molecules are produced after energetic processing (with 200 keV protons) of icy grain mantle analogues. Results. New detections of both molecules towards L1544, L1498, and Elias 18 are reported. The adopted gas phase model is not able to reproduce the observed C2O/C3O ratios, while laboratory experiments show that the ion bombardment of CO-rich mixtures produces C2O/C3O ratios that agree with the observed values. Conclusions. Based on the results obtained here, we conclude that the synthesis of both species is due to the energetic processing of CO-rich icy grain mantles. Their subsequent desorption because of non-thermal processes allows the detection in the gas-phase of young star-forming regions. In more evolved objects, the non-detection of both C2O and C3O is due to their fast destruction in the warm gas.

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

Recently, there has been an increased interest in the study of the generation of low-energy cosmic rays (CRs; < 1 TeV) in shocks situated on the surface of a protostar or along protostellar jets. These locally accelerated CRs offer an attractive explanation for the high levels of non-thermal emission and ionisation rate, $\zeta$, observed close to these sources. The high $\zeta$ observed in some protostellar sources is generally attributed to shock-generated UV photons. The aim of this article is to show that when synchrotron emission and a high $\zeta$ are measured in the same spatial region, a locally shock-accelerated CR flux is sufficient to explain both phenomena. We assume that relativistic particles are accelerated according to the first-order Fermi acceleration mechanism and compute $\zeta$ and the non-thermal emission at cm wavelengths. We then apply our model to the star-forming region OMC-2 FIR 3/FIR 4. Using a Bayesian analysis, we constrain the parameters of the model and estimate the spectral indices of the non-thermal radio emission. We demonstrate that the local CR acceleration model makes it possible to simultaneously explain the synchrotron emission along the HOPS 370 jet within the FIR 3 region and $\zeta$ observed near the FIR 4 protocluster. Our model constrains the magnetic field strength (~250-450$~\mu$G), its turbulent component (~20-40$~\mu$G), and the jet velocity in the shock reference frame for the three non-thermal sources of the HOPS 370 jet (~350-1000 km s$^{-1}$). Beyond the modelling of the OMC-2 FIR 3/FIR 4 system, we show how the combination of continuum observations at cm wavelengths and molecular transitions is a powerful new tool for the analysis of star-forming regions: these two types of observations can be simultaneously interpreted by invoking only the presence of locally accelerated CRs, without having to resort to shock-generated UV photons., Comment: 10 pages, 6 figures, accepted by A&A

Dog domestication is still largely unresolved due to time-gaps in the sampling of regions. Ancient Italian canids are particularly understudied, currently represented by only a few specimens. In the present study, we sampled 27 canid remains from Northern Italy dated between the Late Pleistocene and Bronze Age to assess their genetic variability, and thus add context to dog domestication dynamics. They were targeted at four DNA fragments of the hypervariable region 1 of mitochondrial DNA. A total of 11 samples had good DNA preservation and were used for phylogenetic analyses. The dog samples were assigned to dog haplogroups A, C and D, and a Late Pleistocene wolf was set into wolf haplogroup 2. We present our data in the landscape of ancient and modern dog genetic variability, with a particular focus on the ancient Italian samples published thus far. Our results suggest there is high genetic variability within ancient Italian canids, where close relationships were evident between both a similar to 24,700 years old Italian canid, and Iberian and Bulgarian ancient dogs. These findings emphasize that disentangling dog domestication dynamics benefits from the analysis of specimens from Southern European regions.

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., 9 pages, 2 figures, accepted for publication in MNRAS

Infrared dark clouds (IRDCs) are thought to be potential hosts of the elusive early phases of high-mass star formation. Here, we conduct an in-depth kinematic analysis of one such IRDC, G034.43+00.24 (Cloud F), using high sensitivity and high spectral resolution IRAM-30m N2H+ (1-0) and C18O (1-0) observations. To disentangle the complex velocity structure within this cloud, we use Gaussian decomposition and hierarchical clustering algorithms. We find that four distinct coherent velocity components are present within Cloud F. The properties of these components are compared to those found in a similar IRDC, G035.39-00.33 (Cloud H). We find that the components in both clouds have high densities (inferred by their identification in N2H+), trans-to-supersonic non-thermal velocity dispersions with Mach numbers of ∼1.5-4, a separation in velocity of ∼3 km s-1, and a mean red-shift of ∼0.3 km s-1 between the N2H+ (dense gas) and C18O emission (envelope gas). The latter of these could suggest that these clouds share a common formation scenario. We investigate the kinematics of the larger-scale Cloud F structures, using lower-density-tracing 13CO(1-0) observations. A good correspondence is found between the components identified in the IRAM-30m observations and the most prominent component in the 13CO data. We find that the IRDC Cloud F is only a small part of a much larger structure, which appears to be an inter-arm filament of the Milky Way.

We report new interferometric images of cyclopropenylidene, c-C$_3$H$_2$, towards the young protocluster OMC-2 FIR\,4. The observations were performed at 82 and 85 GHz with the NOrthern Extended Millimeter Array (NOEMA) as part of the project Seeds Of Life In Space (SOLIS). In addition, IRAM-30m data observations were used to investigate the physical structure of OMC-2 FIR\,4. We find that the c-C$_3$H$_2$ gas emits from the same region where previous SOLIS observations showed bright HC$_5$N emission. From a non-LTE analysis of the IRAM-30m data, the c-C$_3$H$_2$ gas has an average temperature of $\sim$40K, a H$_2$ density of $\sim$3$\times$10$^{5}$~cm$^{-3}$, and a c-C$_3$H$_2$ abundance relative to H$_2$ of ($7\pm1$)$\times$10$^{-12}$. In addition, the NOEMA observations provide no sign of significant c-C$_3$H$_2$ excitation temperature gradients across the region (about 3-4 beams), with T$_{ex}$ in the range 8$\pm$3 up to 16$\pm$7K. We thus infer that our observations are inconsistent with a physical interaction of the OMC-2 FIR\,4 envelope with the outflow arising from OMC-2 FIR\,3, as claimed by previous studies. The comparison of the measured c-C$_3$H$_2$ abundance with the predictions from an astrochemical PDR model indicates that OMC-2 FIR\,4 is irradiated by a FUV field $\sim$1000 times larger than the interstellar one, and by a flux of ionising particles $\sim$4000 times larger than the canonical value of $1\times10^{-17}$~s$^{-1}$ from the Galaxy cosmic rays, which is consistent with our previous HC$_5$N observations. This provides an important and independent confirmation of other studies that one or more sources inside the OMC-2 FIR\,4 region emit energetic ($\geq10$~MeV) particles., Accepted for publication in ApJ - 14 pages, 10 figures, 4 tables