Search

This work focuses on the ice features toward the binary protostellar system Ced 110 IRS 4A and 4B, and observed with JWST as part of the Early Release Science Ice Age collaboration. We aim to explore the JWST observations of the binary protostellar system Ced~110~IRS4A and IRS4B to unveil and quantify the ice inventories toward these sources. We compare the ice abundances with those found for the same molecular cloud. The analysis is performed by fitting or comparing laboratory infrared spectra of ices to the observations. Spectral fits are carried out with the ENIIGMA fitting tool that searches for the best fit. For Ced~110~IRS4B, we detected the major ice species H$_2$O, CO, CO$_2$ and NH$_3$. All species are found in a mixture except for CO and CO$_2$, which have both mixed and pure ice components. In the case of Ced~110~IRS4A, we detected the same major species as in Ced~110~IRS4B, as well as the following minor species CH$_4$, SO$_2$, CH$_3$OH, OCN$^-$, NH$_4^+$ and HCOOH. Tentative detection of N$_2$O ice (7.75~$\mu$m), forsterite dust (11.2~$\mu$m) and CH$_3^+$ gas emission (7.18~$\mu$m) in the primary source are also presented. Compared with the two lines of sight toward background stars in the Chameleon I molecular cloud, the protostar has similar ice abundances, except in the case of the ions that are higher in IRS4A. The clearest differences are the absence of the 7.2 and 7.4~$\mu$m absorption features due to HCOO$^-$ and icy complex organic molecules in IRS4A and evidence of thermal processing in both IRS4A and IRS4B as probed by the CO$_2$ ice features. We conclude that the binary protostellar system Ced~110~IRS4A and IRS4B has a large inventory of icy species. The similar ice abundances in comparison to the starless regions in the same molecular cloud suggest that the chemical conditions of the protostar were set at earlier stages in the molecular cloud., Comment: 33 pages, 19 Figures. Accepted for publication in Astronomy & Astrophysics

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

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

Context. The formation of stars has been subject to extensive studies in the past decades from molecular cloud to protoplanetary disk scales. It is still not fully understood how the surrounding material in a protostellar system, that often shows asymmetric structures with complex kinematic properties, feeds the central protostar(s) and their disk(s). Aims. We study the spatial morphology and kinematic properties of the molecular gas surrounding the IRS3A and IRS3B protostellar systems in the L1448N region located in the Perseus molecular cloud. Methods. We present 1 mm NOEMA observations of the PRODIGE large program and analyze the kinematic properties of molecular lines. Given the complexity of the spectral profiles, the lines are fitted with up to three Gaussian velocity components. The clustering algorithm DBSCAN is used to disentangle the velocity components into the underlying physical structure. Results. We discover an extended gas bridge (~3000 au) surrounding both the IRS3A and IRS3B systems in six molecular line tracers (C18O, SO, DCN, H2CO, HC3N, and CH3OH). This gas bridge is oriented along the northeast-southwest direction and shows clear velocity gradients on the order of 100 km/s/pc towards the IRS3A system. We find that the observed velocity profile is consistent with analytical streamline models of gravitational infall towards IRS3A. The high-velocity C18O (2-1) emission towards IRS3A indicates a protostellar mass of ~1.2 Msun. Conclusions. While high angular resolution continuum data often show IRS3A and IRS3B in isolation, molecular gas observations reveal that these systems are still embedded within a large-scale mass reservoir with a complex spatial morphology as well as velocity profiles. The kinematic properties of the extended gas bridge are consistent with gravitational infall toward the IRS3A protostar., Comment: 16 pages, 10 figures; accepted for publication in A&A

The Milky Way's Central Molecular Zone (CMZ) differs dramatically from our local solar neighbourhood, both in the extreme interstellar medium conditions it exhibits (e.g. high gas, stellar, and feedback density) and in the strong dynamics at play (e.g. due to shear and gas influx along the bar). Consequently, it is likely that there are large-scale physical structures within the CMZ that cannot form elsewhere in the Milky Way. In this paper, we present new results from the Atacama Large Millimeter/submillimeter Array (ALMA) large programme ACES (ALMA CMZ Exploration Survey) and conduct a multi-wavelength and kinematic analysis to determine the origin of the M0.8$-$0.2 ring, a molecular cloud with a distinct ring-like morphology. We estimate the projected inner and outer radii of the M0.8$-$0.2 ring to be 79" and 154", respectively (3.1 pc and 6.1 pc at an assumed Galactic Centre distance of 8.2 kpc) and calculate a mean gas density $> 10^{4}$ cm$^{-3}$, a mass of $\sim$ $10^6$ M$_\odot$, and an expansion speed of $\sim$ 20 km s$^{-1}$, resulting in a high estimated kinetic energy ($> 10^{51}$ erg) and momentum ($> 10^7$ M$_\odot$ km s$^{-1}$). We discuss several possible causes for the existence and expansion of the structure, including stellar feedback and large-scale dynamics. We propose that the most likely cause of the M0.8$-$0.2 ring is a single high-energy hypernova explosion. To viably explain the observed morphology and kinematics, such an explosion would need to have taken place inside a dense, very massive molecular cloud, the remnants of which we now see as the M0.8$-$0.2 ring. In this case, the structure provides an extreme example of how supernovae can affect molecular clouds., Comment: 17 pages, 13 figures, and 2 tables. Published in Astronomy & Astrophysics

Star formation in the central molecular zone (CMZ) is suppressed with respect to that of the Galactic disk, and this is likely related to its high turbulent environment. This turbulence impedes the potential detection of prestellar cores. We present the temperature, density, and spatial structure of the CMZ molecular cloud G+0.693-0.027, which has been proposed to host a prestellar cluster in the Sgr B2 region. We analysed multiple HC$_{3}$N rotational transitions that were observed with the IRAM 30m, APEX, Yebes 40m, and GBT radio telescopes, together with SMA+APEX spatially resolved maps. The spectral shape of HC$_{3}$N lines shows three velocity components: a broad component with a line width of 23 km s$^{-1}$ (C1), and two narrow components with line widths of 7.2 and 8.8 km s$^{-1}$ (C2 and C3). This suggests that a fraction of the molecular gas in this cloud is undergoing turbulence dissipation. From a non-local thermodynamic equilibrium analysis we have found H$_{2}$ densities of 2$\times$10$^{4}$ cm$^{-3}$, 5$\times$10$^{4}$ cm$^{-3}$, and 4$\times$10$^{5}$ cm$^{-3}$ and kinetic temperatures of 140 K, 30 K, and 80 K for C1, C2, and C3, respectively. The spatially resolved maps confirm that the colder and high-density condensations C2 and C3, which peak in the 70-85 km s$^{-1}$ velocity range, are embedded in a more diffuse and warmer gas (C1). The larger-scale structure of the Sgr B2 region shows a hole at 40-50 km s$^{-1}$ that is likely due to a small cloud that shocked the Sgr B2 region and is spatially related with a massive cloud at 60-80 km s$^{-1}$. We propose that the impacting small cloud sequentially triggered the formation of Sgr B2(M), (N), and (S) and the condensations in G+0.693-0.027 during its passage. Based on the analysis of the masses of the two condensations and on the virial parameters, C2 might expand, while C3 might further fragment or collapse., Comment: 21 pages, 16 figures, 5 tables, 4 appendix - Accepted for publication in A&A

We report the first detection of the metal-bearing molecules sodium sulfide (NaS) and magnesium sulfide (MgS) and the tentative detection of calcium monoxide (CaO) in the interstellar medium (ISM) towards the Galactic Center molecular cloud G+0.693-0.027. The derived column densities are (5.0+-1.1) x 10$^{10}$ cm$^{-2}$, (6.0+-0.6) x $^{10}$ cm$^{-2}$, and (2.0+-0.5) x $^{10}$ cm$^{-2}$, respectively. This translates into fractional abundances with respect to H$_2$ of (3.7+-1.0) x $10^{-13}$, (4.4+-0.8) x $10^{-13}$, and (1.5+-0.4) x $10^{-13}$, respectively. We have also searched for other Na-, Mg- and Ca-bearing species towards this source but none of them have been detected and thus we provide upper limits for their abundances. We discuss the possible chemical routes involved in the formation of these molecules containing metals under interstellar conditions. Finally, we compare the ratio between sulfur-bearing and oxygen-bearing molecules with and without metals, finding that metal-bearing sulfur molecules are much more abundant than metal-bearing oxygen ones, in contrast with the general trend found in the ratios between other non metal- oxygen- and sulfur-bearing molecules. This further strengthen the idea that sulfur may be little depleted in G+0.693-0.027 as a result of the low velocity shocks present in this source sputtering large amounts of material from dust grains., Comment: 18 pages, 7 tables and 7 figures, final version accepted by ApJ

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

Complex Organic Molecules (COMs) have been found toward low-mass protostars but the origins of the COM emission are still unclear. It can be associated with, for example, hot corinos, outflows, and/or accretion shock/disk atmosphere. We have conducted NOEMA observations toward SVS13A from the PROtostars & DIsks: Global Evolution (PRODIGE) program. Our previous \ce{DCN} observations reveal a possible infalling streamer, which may affect the chemistry of the central protobinary by inducing accretion outbursts and/or shocked gas. Here, we further analyze six O-bearing COMs: CH3OH, aGg'-(CH2OH)2, C2H5OH, CH2(OH)CHO, CH3CHO, and CH3OCHO. Although the COM emission is not spatially resolved, we constrain the source sizes to $\lesssim0.3-0.4$ arcsec (90$-$120 au) by conducting uv-domain Gaussian fitting. Interestingly, the high-spectral resolution data reveal complex line profiles with multiple peaks showing differences between these six O-bearing COMs. The LTE fitting unveils differences in excitation temperatures and emitting areas among these COMs. We further conduct multiple-velocity-component LTE fitting to decompose the line emission into different kinematic components. Up to 6 velocity components are found from the LTE modeling. The temperature, column density, and source size of these components from each COM are obtained. We find a variety in excitation temperatures ($100-500$ K) and source sizes (D$\sim10-70$ au) from these kinematic components from different COMs. The emission of each COM can trace several components and different COMs most likely trace different regions. Given this complex structure, we suggest that the central region is inhomogeneous and unlikely to be heated by only protostellar radiation. We conclude that accretion shocks induced by the large-scale infalling streamer likely exist and contribute to the complexity of the COM emission., Comment: 29 pages, 18 figures, accepted to A&A

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

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

The ALMA interferometer has played a key role in revealing a new component of the Sun-like star forming process: the molecular streamers, i.e. structures up to thousands of au long funneling material non-axisymmetrically to disks. In the context of the FAUST ALMA LP, the archetypical VLA1623-2417 protostellar cluster has been imaged at 1.3 mm in the SO(5$_6$--4$_5$), SO(6$_6$--5$_5$), and SiO(5--4) line emission at the spatial resolution of 50 au. We detect extended SO emission, peaking towards the A and B protostars. Emission blue-shifted down to 6.6 km s$^{-1}$ reveals for the first time a long ($\sim$ 2000 au) accelerating streamer plausibly feeding the VLA1623 B protostar. Using SO, we derive for the first time an estimate of the excitation temperature of an accreting streamer: 33$\pm$9 K. The SO column density is $\sim$ 10$^{14}$ cm$^{-2}$, and the SO/H$_2$ abundance ratio is $\sim$ 10$^{-8}$. The total mass of the streamer is 3 $\times$ 10$^{-3}$ $Msun$, while its accretion rate is 3--5 $\times$ 10$^{-7}$ Msun yr$^{-1}$. This is close to the mass accretion rate of VLA1623 B, in the 0.6--3 $\times$ 10$^{-7}$ Msun yr$^{-1}$ range, showing the importance of the streamer in contributing to the mass of protostellar disks. The highest blue- and red-shifted SO velocities behave as the SiO(5--4) emission, the latter species detected for the first time in VLA1623-2417: the emission is compact (100-200 au), and associated only with the B protostar. The SO excitation temperature is $\sim$ 100 K, supporting the occurrence of shocks associated with the jet, traced by SiO., Comment: accepted by MNRAS

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

We report the first detection of SiC$_2$ in the interstellar medium. The molecule was identified through six rotational transitions toward G\,+0.693$-$0.027, a molecular cloud located in the Galactic center. The detection is based on a line survey carried out with the GBT, the Yebes 40m, and the IRAM 30m telescopes covering a range of frequencies from 12 to 276 GHz. We fit the observed spectra assuming local thermodynamic equilibrium and derive a column density of ($1.02\pm0.04)\times10^{13}$ cm$^{-2}$, which gives a fractional abundance of $7.5\times10^{-11}$ with respect to H$_2$, and an excitation temperature of $5.9\pm0.2$ K. We conclude that SiC$_2$ can be formed in the shocked gas by a reaction between the sputtered atomic silicon and C$_2$H$_2$, or it can be released directly from the dust grains due to disruption. We also search for other Si-bearing molecules and detect eight rotational transitions of SiS and four transitions of Si$^{18}$O. The derived fractional abundances are $3.9\times10^{-10}$ and $2.1\times10^{-11}$, respectively. All Si-bearing species toward G\,+0.693$-$0.027 show fractional abundances well below what is typically found in late-type evolved stars., Comment: 8 pages, 5 figures

We report new observations of molecular anions with the Yebes 40m and IRAM 30m telescopes toward the cold dense clouds TMC-1 CP, Lupus-1A, L1527, L483, L1495B, and L1544. We detected for the first time C3N- and C5N- in Lupus-1A and C4H- and C6H- in L483. In addition, we report new lines of C6H- toward the six targeted sources, of C4H- toward TMC-1 CP, Lupus-1A, and L1527, and of C8H- and C3N- in TMC-1 CP. Excitation calculations indicate that the lines of anions accessible to radiotelescopes run from subthermally excited to thermalized as the size of the anion increases, with the degree of departure from thermalization depending on the H2 volume density and the line frequency. We noticed that the collision rate coefficients available for the radical C6H cannot explain various observational facts, which advises for a revisitation of the collision data for this species. The observations presented here, together with observational data from the literature, are used to model the excitation of interstellar anions and to constrain their abundances. In general, the anion-to-neutral ratios derived here agree within 50 % (a factor of two at most) with literature values, when available, except for the C4H-/C4H ratio, which shows higher differences due to a revision of the dipole moment of C4H. From the set of anion-to-neutral abundance ratios derived two conclusions can be drawn. First, the C6H-/C6H ratio shows a tentative trend in which it increases with increasing H2 density, as expected from theoretical grounds. And second, it is incontestable that the higher the molecular size the higher the anion-to-neutral ratio, which supports a formation mechanism based on radiative electron attachment. Nonetheless, calculated rate coefficients for electron attachment to the medium size species C4H and C3N are probably too high and too low, respectively, by more than one order of magnitude., Comment: Accepted for publication in A&A

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

Cyanides, ranging from three carbon atoms to PAHs, and alkenyl compounds are abundant in the interstellar medium (ISM). Aminoacrylonitrile (3-Amino-2-propenenitrile, H$_{2}$N-CH=CH-CN), an alkenyl cyanide, thus represents a promising candidate for new interstellar detection. A comprehensive spectroscopic laboratory investigation of aminoacrylonitrile in its rotational ground vibrational state has been herein performed. The measurements carried out up to the THz regime made it possible to generate a precise set of reliable rest frequencies for its search in space up to sub-millimetre wavelengths. The $Z$-aminoacrylonitrile ($Z$-apn) isomer spectrum has been recorded employing a source-modulated sub-millimetre spectrometer, from 80 GHz to 1 THz. A combination of Doppler and sub-Doppler measurement regimes allowed to record 600 new lines. The collected data have enabled the characterisation of a set of spectroscopic parameters up to decic centrifugal distortion constants. The catalogue generated from the improved spectral data has been used for the search of $Z$-apn in the spectral survey of the G+0.693-0.027 molecular cloud located in the central molecular zone, in the proximity of the Galactic centre.

Interstellar amides have attracted significant attentions as they are potential precursors for a wide variety of organics essential to life. However, our current understanding of their formation in space is heavily based on observations in star-forming regions and hence the chemical networks lack the constraints on their early origin. In this work, unbiased sensitive spectral surveys with IRAM 30m and Yebes 40m telescopes are used to systematically study a number of amides towards a quiescent Galactic Centre molecular cloud, G+0.693-0.027. We report the first detection of acetamide (CH3C(O)NH2) and trans-N-methylformamide (CH3NHCHO) towards this cloud. In addition, with the wider frequency coverage of the survey, we revisited the detection of formamide (NH2CHO) and urea (carbamide; NH2C(O)NH2), which had been reported previously towards G+0.693-0.027. Our results are compared with those present in the literature including recent laboratory experiments and chemical models. We find constant abundance ratios independently of the evolutionary stages, suggesting that amides related chemistry is triggered in early evolutionary stages of molecular cloud and remain unaffected by the warm-up phase during the star formation process. Although a correlation between more complex amides and NH2CHO have been suggested, alternative formation routes involving other precursors such as acetaldehyde (CH3CHO), methyl isocyanate (CH3NCO) and methylamine (CH3NH2) may also contribute to the production of amides. Observations of amides together with these species towards a larger sample of sources can help to constrain the amide chemistry in the interstellar medium., Comment: 16 pages, 8 figures, 3 tables, accepted for publication in Monthly Notices of the Royal Astronomical Society

Methylamine has been the only simple alkylamine detected in the interstellar medium for a long time. With the recent secure and tentative detections of vinylamine and ethylamine, respectively, dimethylamine has become a promising target for searches in space. Its rotational spectrum, however, has been known only up to 45 GHz until now. Here we investigate the rotation-tunneling spectrum of dimethylamine in selected regions between 76 and 1091 GHz using three different spectrometers in order to facilitate its detection in space. The quantum number range is extended to $J = 61$ and $K_a = 21$, yielding an extensive set of accurate spectroscopic parameters. To search for dimethylamine, we refer to the spectral line survey ReMoCA carried out with the Atacama Large Millimeter/submillimeter Array toward the high-mass star-forming region Sagittarius B2(N) and a spectral line survey of the molecular cloud G+0.693$-$0.027 employing the IRAM 30 m and Yebes 40 m radio telescopes. We report nondetections of dimethylamine toward the hot molecular cores Sgr B2(N1S) and Sgr B2(N2b) as well as G+0.693$-$0.027 which imply that dimethylamine is at least 14, 4.5 and 39 times less abundant than methylamine toward these sources, respectively. The observational results are compared to computational results from a gas-grain astrochemical model. The modeled methylamine to dimethylamine ratios are compatible with the observational lower limits. However, the model produces too much ethylamine compared with methylamine which could mean that the already fairly low levels of dimethylamine in the models may also be too high., Comment: Mon. Not. R. Astron. Soc., accepted. 33 pages including tables, figures, and appendix

More than 50% of solar-mass stars form in multiple systems. It is therefore crucial to investigate how multiplicity affects the star and planet formation processes at the protostellar stage. We report continuum and C$^{18}$O (2-1) observations of the VLA 1623-2417 protostellar system at 50 au angular resolution as part of the ALMA Large Program FAUST. The 1.3 mm continuum probes the disks of VLA 1623A, B, and W, and the circumbinary disk of the A1+A2 binary. The C$^{18}$O emission reveals, for the first time, the gas in the disk-envelope of VLA 1623W. We estimate the dynamical mass of VLA 1623W, $M_{\rm dyn}=0.45\pm0.08$ M$_{\odot}$, and the mass of its disk, $M_{\rm disk}\sim6\times10^{-3}$ M$_{\odot}$. C$^{18}$O also reveals streamers that extend up to 1000 au, spatially and kinematically connecting the envelope and outflow cavities of the A1+A2+B system with the disk of VLA 1623W. The presence of the streamers, as well as the spatial ($\sim$1300 au) and velocity ($\sim$2.2 km/s) offset of VLA 1623W suggest that either sources W and A+B formed in different cores, interacting between them, or that source W has been ejected from the VLA 1623 multiple system during its formation. In the latter case, the streamers may funnel material from the envelope and cavities of VLA 1623AB onto VLA 1623W, thus concurring to set its final mass and chemical content., Comment: 10 pages, 5 figures, accepted by MNRAS

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

Icy grain mantles are the main reservoir of the volatile elements that link chemical processes in dark, interstellar clouds with the formation of planets and composition of their atmospheres. The initial ice composition is set in the cold, dense parts of molecular clouds, prior to the onset of star formation. With the exquisite sensitivity of JWST, this critical stage of ice evolution is now accessible for detailed study. Here we show the first results of the Early Release Science program "Ice Age" that reveal the rich composition of these dense cloud ices. Weak ices, including, $^{13}$CO$_2$, OCN$^-$, $^{13}$CO, OCS, and COMs functional groups are now detected along two pre-stellar lines of sight. The $^{12}$CO$_2$ ice profile indicates modest growth of the icy grains. Column densities of the major and minor ice species indicate that ices contribute between 2 and 19% of the bulk budgets of the key C, O, N, and S elements. Our results suggest that the formation of simple and complex molecules could begin early in a water-ice rich environment., Comment: To appear in Nature Astronomy on January 23rd, 2023. 33 pages, 16 figures, 3 tables; includes extended and supplemental data sections. Part of the JWST Ice Age Early Release Science program's science enabling products. Enhanced spectra downloadable on Zenodo at the following DOI: 10.5281/zenodo.7501239

Gas phase Elemental abundances in molecular CloudS (GEMS) is an IRAM 30m large program aimed at determining the elemental abundances of carbon (C), oxygen (O), nitrogen (N), and sulfur (S) in a selected set of prototypical star-forming filaments. In particular, the elemental abundance of S remains uncertain by several orders of magnitude and its determination is one of the most challenging goals of this program. We have carried out an extensive chemical modeling of the fractional abundances of CO, HCO$^+$, HCN, HNC, CS, SO, H$_2$S, OCS, and HCS$^+$ to determine the sulfur depletion toward the 244 positions in the GEMS database. These positions sample visual extinctions from A$_V$ $\sim$ 3 mag to $>$50 mag, molecular hydrogen densities ranging from a few 10$^3$~cm$^{-3}$ to 3$\times$10$^6$~cm$^{-3}$, and T$_k$ $\sim$ 10$-$35 K. Most of the positions in Taurus and Perseus are best fitted assuming early-time chemistry, t=0.1 Myr, $\zeta_{H_2}$$\sim$ (0.5$-$1)$\times$10$^{-16}$ s$^{-1}$, and [S/H]$\sim$1.5$\times$10$^{-6}$. On the contrary, most of the positions in Orion are fitted with t=1~Myr and $\zeta_{H_2}$$\sim$ 10$^{-17}$ s$^{-1}$. Moreover, $\sim$40% of the positions in Orion are best fitted assuming the undepleted sulfur abundance, [S/H]$\sim$1.5$\times$10$^{-5}$. Our results suggest that sulfur depletion depends on the environment. While the abundances of sulfur-bearing species are consistent with undepleted sulfur in Orion, a depletion factor of $\sim$20 is required to explain those observed in Taurus and Perseus. We propose that differences in the grain charge distribution in the envelopes of the studied clouds might explain these variations. The shocks associated with past and ongoing star formation could also contribute to enhance [S/H] in Orion., Comment: 22 pages, 15 figures, Astronomy and Astrophysics, in press

Context. In the past few years, there has been a rise in the detection of streamers, asymmetric flows of material directed toward the protostellar disk with material from outside the star's natal core. It is unclear how they affect the process of mass accretion, in particular beyond the Class 0 phase. Aims. We investigate the gas kinematics around Per-emb-50, a Class I source in the crowded star-forming region NGC 1333. Our goal is to study how the mass infall proceeds from envelope to disk scales in this source. Results. We discover a streamer delivering material toward Per-emb-50 in H$_2$CO and C$^{18}$O emission. The streamer's emission can be well described by the analytic solutions for an infalling parcel of gas along a streamline with conserved angular momentum, both in the image plane and along the line of sight velocities. The streamer has a mean infall rate of $1.3 \times 10^{ -6}$ M$_{ \odot}$ yr$^{ -1}$, $5 -10$ times higher than the current accretion rate of the protostar. SO and SO$_2$ emission reveal asymmetric infall motions in the inner envelope, additional to the streamer around Per-emb-50. Furthermore, the presence of SO$_2$ could mark the impact zone of the infalling material. Conclusions. The streamer delivers sufficient mass to sustain the protostellar accretion rate and might produce an accretion burst, which would explain the protostar's high luminosity with respect to other Class I sources. Our results highlight the importance of late infall for protostellar evolution: streamers might provide a significant amount of mass for stellar accretion after the Class 0 phase., Comment: 20 pages, 14 figures, accepted for publication in A&A

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

In recent times, large organic molecules of exceptional complexity have been found in diverse regions of the interstellar medium. In this context, we aim to provide accurate frequencies of the ground vibrational state of two key aliphatic aldehydes, n-butanal and its branched-chain isomer, i-butanal, to enable their eventual detection in the interstellar medium. We employ a frequency modulation millimeter-wave absorption spectrometer to measure the rotational features of n- and i-butanal. We use the spectral line survey ReMoCA performed with the Atacama Large Millimeter/submillimeter Array to search for n- and i-butanal toward the star-forming region Sgr B2(N). We also search for both aldehydes toward the molecular cloud G+0.693-0.027 with IRAM 30 m and Yebes 40 m observations. Several thousand rotational transitions belonging to the lowest-energy conformers have been assigned in the laboratory spectra up to 325 GHz. A precise set of the relevant rotational spectroscopic constants has been determined for each structure. We report non-detections of n- and i-butanal toward both sources, Sgr B2(N1S) and G+0.693-0.027. We find that n- and i-butanal are at least 2-6 and 6-18 times less abundant than acetaldehyde toward Sgr B2(N1S), respectively, and that n-butanal is at least 63 times less abundant than acetaldehyde toward G+0.693-0.027. Comparison with astrochemical models indicates good agreement between observed and simulated abundances (where available). Grain-surface chemistry appears sufficient to reproduce aldehyde ratios in G+0.693-0.027; gas-phase production may play a more active role in Sgr B2(N1S). Our astronomical results indicate that the family of interstellar aldehydes in the Galactic center region is characterized by a drop of one order of magnitude in abundance at each incrementation in the level of molecular complexity.

The Central Molecular Zone (CMZ) contains most of the mass of our Galaxy but its star formation rate is one order of magnitude lower than in the Galactic disc. This is likely related to the fact that the bulk of the gas in the CMZ is in a warm ($>$100 K) and turbulent phase with little material in the pre-stellar phase. We present in this Letter observations of deuterium fractionation (D/H ratios) of HCN, HNC, HCO$^{+}$, and N$_{2}$H$^{+}$ towards the CMZ molecular cloud G+0.693-0.027. These observations clearly show, for the first time, the presence of a colder, denser, and less turbulent narrow component, with a line width of $\sim$9 km s$^{-1}$, in addition to the warm, less dense and turbulent broad component with a line width of $\sim$20 km s$^{-1}$. The very low D/H ratio $\le$6$\times$10$^{-5}$ for HCO$^{+}$ and N$_{2}$H$^{+}$, close to the cosmic value ($\sim$2.5$\times$10$^{-5}$), and the high D/H ratios $>$4$\times$10$^{-4}$ for HCN and HNC derived for the broad component, confirm the presence of high-temperatures deuteration routes for nitriles. For the narrow component we have derived D/H ratios $>$10$^{-4}$ and excitation temperatures of $7$ K for all molecules, suggesting kinetic temperatures $\le$30 K and H$_2$ densities $\ge$5$\times$10$^{4}$ cm$^{-3}$, at least one order of magnitude larger than for the broad component. The method presented in this Letter allows to identify clouds on the verge of star formation, i.e. under pre-stellar conditions, towards the CMZ. This method can also be used for the identification of such clouds in external galaxies., Comment: 11 pages, 4 figures, 2 tables, 2 appendix - Accepted for publication in ApJ Letters

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

We studied theoretically the cis-trans isomerization reactions of two astrophysically relevant acids, formic acid (HCOOH) and thioformic acid (HC(O)SH), where the latter has recently been discovered in space. We also searched for these molecules towards the hot core G31.41+0.31 to compare their abundances with the expected theoretical isomerization results. Our results demonstrate that these isomerizations are viable in the conditions of the ISM due to ground-state tunneling effects, which allow the system to reach the thermodynamic equilibrium at moderately low temperatures. At very low temperatures (Tkin 10 K), the reaction rate constants for the cis to trans isomerizations are very small, which implies that the cis isomers should not be formed under cold ISM conditions. This is in disagreement with observations of the cis/trans isomers of HCOOH in cold cores where the cis isomer is found to be 5-6% the trans isomer. At high temperatures (150-300 K), our theoretical data not only match the observed behaviour of the trans/cis abundance ratios for HCOOH (the cis form is undetected), but they support our tentative detection of the trans and, for the first time in the insterstellar medium, of the cis isomer of HC(O)SH towards the hot molecular core G31.41+0.31 (with a measured trans/cis abundance ratio of 3.7). While the trans/cis ratio for HC(O)SH in the ISM depends on the relative stability of the isomers, the trans/cis ratio for HCOOH cannot be explained by isomerization, and it is determined by other competitive chemical processes, Comment: 14 pages, 7 figures, 5 tables

Little is known about the chemistry of isocyanates (compounds with the functional group R-N=C=O) in the interstellar medium, as only four of them have been detected so far: isocyanate radical (NCO), isocyanic acid (HNCO), N-protonated isocyanic acid (H$_2$NCO$^+$) and methyl isocyanate (CH$_3$NCO). The molecular cloud G+0.693-0.027, located in the Galactic Centre, represents an excellent candidate to search for new isocyanates since it exhibits high abundances of the simplest ones, HNCO and CH$_3$NCO. After CH$_3$NCO, the next complex isocyanates are ethyl isocyanate (C$_2$H$_5$NCO) and vinyl isocyanate (C$_2$H$_3$NCO). Their detection in the ISM would enhance our understanding of the formation of these compounds in space. We have detected C$_2$H$_5$NCO and H$_2$NCO$^+$ towards G+0.693-0.027 (the former for the first time in the interstellar medium) with molecular abundances of (4.7$-$7.3)$\times$10$^{-11}$ and (1.0$-$1.5)$\times$10$^{-11}$, respectively. A ratio CH$_3$NCO / C$_2$H$_5$NCO = 8$\pm$1 is obtained; therefore the relative abundance determined for HNCO:CH$_3$NCO:C$_2$H$_5$NCO is 1:1/55:1/447, which implies a decrease by more than one order of magnitude going progressively from HNCO to CH$_3$NCO and to C$_2$H$_5$NCO. This is similar to what has been found for e.g. alcohols and thiols and suggests that C$_2$H$_5$NCO is likely formed on the surface of dust grains. In addition, we have obtained column density ratios of HNCO / NCO > 269, HNCO / H$_2$NCO$^+$ $\sim$ 2100 and C$_2$H$_3$NCO / C$_2$H$_5$NCO~<~4. A comparison of the Methyl~/~Ethyl ratios for isocyanates (-NCO), alcohols (-OH), formiates (HCOO-), nitriles (-CN) and thiols (-SH) is performed and shows that ethyl-derivatives may be formed more efficiently for the N-bearing molecules than for the O- and S-bearing molecules., Comment: 5 pages, 5 figures; accepted in A&A

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

The chemical and physical evolution of starless and pre-stellar cores are of paramount importance to understanding the process of star formation. The Taurus Molecular Cloud cores TMC 1-C and TMC 1-CP share similar initial conditions and provide an excellent opportunity to understand the evolution of the pre-stellar core phase. We investigated the evolutionary stage of starless cores based on observations towards the prototypical dark cores TMC 1-C and TMC 1-CP, mapping them in the CS $3\rightarrow 2$, C$^{34}$S $3\rightarrow 2$, $^{13}$CS $2\rightarrow 1$, DCN $1\rightarrow 0$, DCN $2\rightarrow 1$, DNC $1\rightarrow 0$, DNC $2\rightarrow 1$, DN$^{13}$C $1\rightarrow 0$, DN$^{13}$C $2\rightarrow 1$, N$_2$H$^+$ $1\rightarrow 0$, and N$_2$D$^+$ $1\rightarrow 0$ transitions. We performed a multi-transitional study of CS and its isotopologs, DCN, and DNC lines to characterize the physical and chemical properties of these cores. We studied their chemistry using the state-of-the-art gas-grain chemical code Nautilus and pseudo time-dependent models to determine their evolutionary stage. Observational diagnostics seem to indicate that TMC 1-C is in a later evolutionary stage than TMC 1-CP, with a chemical age $\sim$1 Myr. TMC 1-C shows signs of being an evolved core at the onset of star formation, while TMC 1-CP appears to be in an earlier evolutionary stage due to a more recent formation or, alternatively, a collapse slowed down by a magnetic support., Comment: 26 pages, 10 figures

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

We report here the first detection in the interstellar medium of the cyanomidyl radical (HNCN). Using the Yebes 40m and the IRAM 30m telescopes, we have targeted the doublets of the $N$=2$-$1, 4$-$3, 5$-$4, 6$-$5, and 7$-$6 transitions of HNCN toward the molecular cloud G+0.693-0.027. We have detected three unblended lines of HNCN, these are the $N$=6$-$5 doublet and one line of the $N$=4$-$3 transition. Additionally we present one line of the $N$=5$-$4 transition partially blended with emission from other species. The Local Thermodynamic Equilibrium best fit to the data gives a molecular abundance of (0.91$\pm$0.05)$\times$10$^{-10}$ with respect to H$_2$. The relatively low abundance of this species in G+0.693-0.027, and its high reactivity, suggest that HNCN is possibly produced by gas-phase chemistry. Our work shows that this highly reactive molecule is present in interstellar space, and thus it represents a plausible precursor of larger prebiotic molecules with the NCN backbone such as cyanamide (NH$_2$CN), carbodiimide (HNCNH) and formamidine (NH$_2$CHNH)., Comment: Accepted in Monthly Notices of the Royal Astronomical Society Letters

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

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

Gas phase Elemental abundances in Molecular CloudS (GEMS) is an IRAM 30m Large Program designed to estimate the S, C, N, and O depletions and gas ionization degree, X(e-), in a set of star-forming filaments of Taurus, Perseus and Orion. Our immediate goal is to build up a complete database of molecular abundances that can serve as an observational basis for estimating X(e-) and the C, O, N, and S depletions through chemical modeling. We observed and derived the abundances of 14 species (13CO, C18O, HCO+, H13CO+, HC18O+, HCN, H13CN, HNC, HCS+, CS, SO, 34SO, H2S, and OCS) in 244 positions, covering the AV 3 to 100 mag, n(H2) a few 10$^{3}$ to 10$^6$ cm$^{-3}$, and Tk 10 to 30 K ranges in these clouds, avoiding protostars, HII regions, and outflows. A statistical analysis is carried out to identify general trends between different species and with physical parameters. Relations between molecules reveal strong linear correlations which define three different families: (1) 13CO and C18O; (2) H13CO+, HC18O+, H13CN, and HNC; and (3) the S-bearing molecules. The abundances of the CO isotopologs increase with the gas kinetic temperature until TK 15 K. For higher temperatures, the abundance remains constant with a scatter of a factor of 3. The abundances of H13CO+, HC18O+, H13CN, and HNC are well correlated with each other, and all of them decrease with molecular hydrogen density, following the law n(H2)$^{-0.8\pm0.2}$. The abundances of S-bearing species also decrease with n(H2) at a rate of (S-bearing/H)gas n(H2)$^{-0.6\pm0.1}$. The abundances of molecules belonging to groups 2 and 3 do not present any clear trend with gas temperature. At scales of molecular clouds, the C18O abundance is the quantity that better correlates with the cloud mass. We discuss the utility of the 13CO/C18O, HCO+/H13CO+, and H13CO+/H13CN abundance ratios as chemical diagnostics of star formation in external galaxies., Comment: 28 pages, 22 figures. Accepted for publication in A&A

Methanol and complex organic molecules have been found in cold starless cores, where a standard warm-up scenario would not work because of the absence of heat sources. A recent chemical model attributed the presence of methanol and large organics to the efficient chemical desorption and a class of neutral-neutral reactions that proceed fast at low temperatures in the gas phase. The model calls for a high abundance of methanol ice at the edge of the CO freeze-out zone in cold cloud cores. We performed medium resolution spectroscopy toward 3 field stars behind the starless core L1544 at 3 $\mu$m to constrain the methanol ice abundance and compare it with the model predictions. One of the field stars shows a methanol-ice abundance of 11% with respect to water ice. This is higher than the typical methanol abundance previously found in cold cloud cores (4%), but is 4.5 times smaller than predicted. The reason for the disagreement between the observations and the model calculations is not yet understood., Comment: 9 pages, 6 figures, accepted to A&A

CS is among the most abundant gas-phase S-bearing molecules in cold dark molecular clouds. It is easily observable with several transitions in the millimeter wavelength range, and has been widely used as a tracer of the gas density in the interstellar medium in our Galaxy and external galaxies. Chemical models fail to account for the observed CS abundances when assuming the cosmic value for the elemental abundance of sulfur. The CS+O -> CO + S reaction has been proposed as a relevant CS destruction mechanism at low temperatures, and could explain the discrepancy between models and observations. Its reaction rate has been experimentally measured at temperatures of 150-400 K, but the extrapolation to lower temperatures is doubtful. Here we calculate the CS+O reaction rate at temperatures <150 K which are prevailing in the interstellar medium. We performed ab initio calculations to obtain the three lowest PES of the CS+O system. These PESs are used to study the reaction dynamics, using several methods to eventually calculate the CS+O thermal reaction rates. We compare the results of our theoretical calculations for 150-400 K with those obtained in the laboratory. Our detailed theoretical study on the CS+O reaction, which is in agreement with the experimental data obtained at 150-400 K, demonstrates the reliability of our approach. After a careful analysis at lower temperatures, we find that the rate constant at 10 K is negligible, which is consistent with the extrapolation of experimental data using the Arrhenius expression. We use the updated chemical network to model the sulfur chemistry in TMC1 based on molecular abundances determined from GEMS project observations. In our model, we take into account the expected decrease of the cosmic ray ionization rate along the cloud. The abundance of CS is still overestimated when assuming the cosmic value for the sulfur abundance., Comment: 11 pages, 10 figures, accepted in Astronomy & Astrophysics

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

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

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