Theories of massive star formation predict that massive protostars accrete gas through circumstellar disks. Although several cases have been found already thanks to high-angular resolution interferometry, it remains unknown the internal physical structure of these disks and, in particular, whether they present warps or internal holes as observed in low-mass proto-planetary disks. Here, we report very high angular resolution observations of the H21alpha radio recombination line carried out in Band 9 with the Atacama Large Millimeter/submillimeter Array (beam of 80 mas x 60 mas, or 70 au x 50 au) toward the IRS2 massive young stellar object in the Monoceros R2 star-forming cluster. The H21alpha line shows maser amplification, which allows us to study the kinematics and physical structure of the ionised gas around the massive protostar down to spatial scales of ~1-2 au. Our ALMA images and 3D radiative transfer modelling reveal that the ionized gas around IRS2 is distributed in a Keplerian circumstellar disk and an expanding wind. The H21alpha emission centroids at velocities between -10 and 20 km s-1 deviate from the disk plane, suggesting a warping for the disk. This could be explained by the presence of a secondary object (a stellar companion or a massive planet) within the system. The ionized wind seems to be launched from the disk surface at distances ~11 au from the central star, consistent with magnetically-regulated disk wind models. This suggests a similar wind launching mechanism to that recently found for evolved massive stars such as MWC349A and MWC922.

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.

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.

Small imines containing up to three carbon atoms are present in the interstellar medium. As alkynyl compounds are abundant in this medium, propargylimine thus represents a promising candidate for a new interstellar detection. The goal of the present work is to perform a comprehensive laboratory investigation of the rotational spectrum of propargylimine in its ground vibrational state in order to obtain a highly precise set of rest-frequencies and to search it in space. The rotational spectra of $E$ and $Z$ geometrical isomers of propargylimine have been recorded in laboratory in the 83-500 GHz frequency interval. The measurements have been performed using a source-modulation millimetre-wave spectrometer equipped with a pyrolysis system for the production of unstable species. High-level ab initio calculations were performed to assist the analysis and to obtain reliable estimates for an extended set of spectroscopic quantities. The improved spectral data allow us to perform a successful search for this new imine in the quiescent G+0.693-0.027 molecular cloud. Eighteen lines of $Z$-propargylimine have been detected at level $>2.5 \sigma$, resulting in a column density estimate of $N = (0.24\pm 0.02)\times 10^{14}$ cm$^{-2}$. An upper limit was retrieved for the higher-energy $E$ isomer, which was not detected in the data. The fractional abundance (w.r.t. H$_2$) derived for $Z$-propargylimine is $1.8\times 10^{-10}$. We discuss the possible formation routes by comparing the derived abundance with those measured in the source for possible chemical precursors.

In this work, we present the results of our investigation into the chemistry of Z- and E-cyanomethanimine (HNCHCN), both of which are possible precursors to the nucleobase adenine. Ab initio quantum chemical calculations for a number of reactions with atomic hydrogen were carried out. We find that the reaction H + Z/E-HNCHCN leading both to H-addition as well as H$_2$-abstraction proceed via similar short-range barriers with bimolecular rate coefficients on the order of $\sim10^{-17}$ cm$^{3}$ s$^{-1}$. These results were then incorporated into astrochemical models and used in simulations of the giant molecular cloud G+0.693. The calculated abundances obtained from these models were compared with previous observational data and found to be in good agreement, with a predicted [Z/E] ratio of $\sim3$ - somewhat smaller than the previously derived value of $6.1\pm2.4$. We find that the [Z/E] ratio in our simulations is due mostly to ion-molecule destruction rates driven by the different permanent dipoles of the two conformers. Based on these results, we propose a general rule-of-thumb for estimating the abundances of isomers in interstellar environments.

The lines of HOC$^+$, HCO and CO$^+$ are considered good tracers of photon-dominated regions (PDRs) and X-ray dominated regions. We study these tracers towards regions of the Sgr B2 cloud selected to be affected by different heating mechanisms. We find the lowest values of the column density ratios of HCO$^+$ versus HOC$^+$, HCO and CO$^+$ in dense HII gas, where UV photons dominate the heating and chemistry of gas. HOC$^+$, HCO and CO$^+$ abundances and the above ratios are compared with those of chemical modeling, finding that high temperature chemistry, a cosmic-ray ionization rate of 10$^{-16}$ s$^{-1}$ and timescales $>$10$^{5.0}$ years explain well the HOC$^+$ abundances in quiescent Sgr B2 regions, while shocks are also needed to explain the highest HCO abundances derived for these regions. CO$^+$ is mainly formed in PDRs since the highest CO$^+$ abundances of $\sim$(6-10)$\times$10$^{-10}$ are found in HII regions with electron densities $>$540 cm$^{-3}$ and that CO$^+$ emission is undetected in quiescent gas. Between the ratios, the HCO$^+$/HCO ratio is sensitive to the electron density as it shows different values in dense and diffuse HII regions. We compare SiO J=2-1 emission maps of Sgr B2 with X-ray maps from 2004 and 2012. One known spot shown on the 2012 X-ray map is likely associated with molecular gas at velocities of 15-25 km s$^{-1}$. We also derive the X-ray ionization rate of $\sim$10$^{-19}$ s$^{-1}$ for Sgr B2 regions pervaded by X-rays in 2004, which is quite low to affect the chemistry of the molecular gas.

In the past decade, Astrochemistry has witnessed an impressive increase in the number of detections of complex organic molecules. Some of these species are of prebiotic interest such as glycolaldehyde, the simplest sugar, or amino acetonitrile, a possible precursor of glycine. Recently, we have reported the detection of two new nitrogen-bearing complex organics, glycolonitrile and Z-cyanomethanimine, known to be intermediate species in the formation process of ribonucleotides within theories of a primordial ribonucleic acid (RNA)-world for the origin of life. In this paper, we present deep and high-sensitivity observations toward two of the most chemically rich sources in the Galaxy: a Giant Molecular Cloud in the center of the Milky Way (G+0.693-0.027) and a proto-Sun (IRAS16293-2422 B). Our aim is to explore whether the key precursors considered to drive the primordial RNA-world chemistry, are also found in space. Our high-sensitivity observations reveal that urea is present in G+0.693-0.027 with an abundance of about 5x10-11. This is the first detection of this prebiotic species outside a star-forming region. Urea remains undetected toward the proto-Sun IRAS16293-2422 B (upper limit to its abundance of less than 2x10-11). Other precursors of the RNA-world chemical scheme such as glycolaldehyde or cyanamide are abundant in space, but key prebiotic species such as 2- amino-oxazole, glyceraldehyde or dihydroxyacetone are not detected in either source. Future more sensitive observations targeting the brightest transitions of these species will be needed to disentangle whether these large prebiotic organics are certainly present in space.

Using $0.2^{\prime \prime}$ ($\sim3$ pc) ALMA images of vibrationally excited HC$_3$N emission (HC$_3$N$^*$) we reveal the presence of $8$ unresolved Super Hot Cores (SHCs) in the inner $160$ pc of NGC\,253. Our LTE and non-LTE modelling of the HC$_3$N$^*$ emission indicate that SHCs have dust temperatures of $200-375$ K, relatively high H$_2$ densities of $1-6\times 10^{6}$ cm$^{-3}$ and high IR luminosities of $0.1-1\times 10^8$ L$_\odot$. As expected from their short lived phase ($\sim 10^4$ yr), all SHCs are associated with young Super Star Clusters (SSCs). We use the ratio of luminosities form the SHCs (protostar phase) and from the free-free emission (ZAMS star phase), to establish the evolutionary stage of the SSCs. The youngest SSCs, with the larges ratios, have ages of a few $10^4$ yr (proto-SSCs) and the more evolved SSCs are likely between $10^5$ and $10^6$ yr (ZAMS-SSCs). The different evolutionary stages of the SSCs are also supported by the radiative feedback from the UV radiation as traced by the HNCO/CS ratio, with this ratio being systematically higher in the young proto-SSCs than in the older ZAMS-SSCs. We also estimate the SFR and the SFE of the SSCs. The trend found in the estimated SFE ($\sim40\%$ for proto-SSCs and $>85\%$ for ZAMS-SSCs) and in the gas mass reservoir available for star formation, one order of magnitude higher for proto-SSCs, suggests that star formation is still going on in proto-SSCs. We also find that the most evolved SSCs are located, in projection, closer to the center of the galaxy than the younger proto-SSCs, indicating an inside-out SSC formation scenario.

Context. As a building block for amino acids, formamide (NH$_2$CHO) is an important molecule in astrobiology and astrochemistry, but its formation path in the interstellar medium is not understood well. Aims. We aim to find empirical evidence to support the chemical relationships of formamide to HNCO and H$_2$CO. Methods. We examine high angular resolution (~0.2") Atacama Large Millimeter/submillimeter Array (ALMA) maps of six sources in three high-mass star-forming regions and compare the spatial extent, integrated emission peak position, and velocity structure of HNCO and H$_2$CO line emission with that of NH$_2$CHO by using moment maps. Through spectral modeling, we compare the abundances of these three species. Results. In these sources, the emission peak separation and velocity dispersion of formamide emission is most often similar to HNCO emission, while the velocity structure is generally just as similar to H$_2$CO and HNCO (within errors). From the spectral modeling, we see that the abundances between all three of our focus species are correlated, and the relationship between NH$_2$CHO and HNCO reproduces the previously demonstrated abundance relationship. Conclusions. In this first interferometric study, which compares two potential parent species to NH$_2$CHO, we find that all moment maps for HNCO are more similar to NH$_2$CHO than H$_2$CO in one of our six sources (G24 A1). For the other five sources, the relationship between NH$_2$CHO, HNCO, and H$_2$CO is unclear as the different moment maps for each source are not consistently more similar to one species as opposed to the other.

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

Phosphorus-bearing species are an essential key to form life on Earth, however they have barely been detected in the interstellar medium. Since only PN and PO have been identified so far towards star-forming regions, the chemical formation pathways of P-bearing molecules 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, i.e. the depletion level of P. In order to overcome this problem, we study P-bearing species in diffuse/translucent clouds. In these objects phosphorus is mainly in the gas phase and therefore the elemental initial abundance needed in our chemical simulations corresponds to the cosmic one and is thus well constrained. An advanced chemical model with an updated P-chemistry network has been used. Single-pointing observations were performed with the IRAM 30m telescope towards the line of sight to the blazar B0355+508 aiming for the (2-1) transitions of PN, PO, HCP and CP. This line of sight incorporates five diffuse/translucent clouds. The (2-1) transitions of the PN, PO, HCP and CP were not detected. We report detections of the (1-0) lines of $\mathrm{^{13}CO}$, HNC and CN along with a first detection of $\mathrm{C^{34}S}$. We have reproduced the observations of HNC, CN, CS and CO in every cloud by applying typical conditions for diffuse/translucent clouds. According to our best-fit model, the most abundant P-bearing species are HCP and CP ($\sim10^{-10}$), followed by PN, PO and $\mathrm{PH_3}$ ($\sim10^{-11}$). We show that the production of P-bearing species is favoured towards translucent rather than diffuse clouds, where the environment provides a stronger shielding from the interstellar radiation. Based on our improved model, the (1-0) transitions of HCP, CP, PN and PO are expected to be detectable with estimated intensities up to $\sim200$ mK.

In this paper we present the detailed formalism at the core of the Spectral Line Identification and Modelling (SLIM) within the MAdrid Data CUBe Analysis (MADCUBA) package and their main data handling functionalities. These tools have been developed to visualize, analyze and model large spectroscopic data cubes. We present the highly interactive on-the-fly visualization and modelling tools of MADCUBA and SLIM, which includes an stand-alone spectroscopic database. The parameters stored therein are used to solve the full radiative transfer equation under Local Thermodynamic Equilibrium (LTE). SLIM provides tools to generate synthetic LTE model spectra based on input physical parameters of column density, excitation temperature, velocity, line width and source size. SLIM also provides an automatic fitting algorithm to obtain the physical parameters (with their associated errors) better fitting the observations. Synthetic spectra can be overlayed in the data cubes/spectra to easy the task of multi-molecular line identification and modelling.We present the Java-based MADCUBA and its internal module SLIM packages which provide all the necessary tools for manipulation and analysis of spectroscopic data cubes. We describe in detail the spectroscopic fitting equations and make use of this tool to explore the breaking conditions and implicit errors of commonly used approximations in the literature. Easy-to-use tools like MADCUBA allow the users to derive the physical information from spectroscopic data without the need of resourcing to simple approximations. SLIM allows to use the full radiative transfer equation, and to interactively explore the space of physical parameters and associated uncertainties from observational data.

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.

We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation.

Theories of a pre-RNA world suggest that glycolonitrile (HOCH$_2$CN) is a key species in the process of ribonucleotide assembly, which is considered as a molecular precursor of nucleic acids. In this Letter, we report the first detection of this pre-biotic molecule in the interstellar medium (ISM) by using ALMA data obtained at frequencies between 86.5$\,$GHz and 266.5$\,$GHz toward the Solar-type protostar IRAS16293-2422 B. A total of 15 unblended transitions of HOCH$_2$CN were identified. Our analysis indicates the presence of a cold (T$\rm _{ex}$=24$\pm$8$\,$K) and a warm (T$\rm _{ex}$=158$\pm$38$\,$K) component meaning that this molecule is present in both the inner hot corino and the outer cold envelope of IRAS16293 B. The relative abundance with respect to H$_2$ is (6.5$\pm$0.6)$\times$10$^{-11}$ and $\geq$(6$\pm$2)$\times$10$^{-10}$ for the warm and cold components respectively. Our chemical modelling seems to underproduce the observed abundance for both the warm and cold component under various values of the cosmic-ray ionisation rate ($\zeta$). Key gas phase routes for the formation of this molecule might be missing in our chemical network.

We previously observed the star-forming region G16.59$-$0.05 through interferometric observations of both thermal and maser lines, and identified a high-mass young stellar object (YSO) which is surrounded by an accretion disk and drives a nonthermal radio jet. We performed high-angular-resolution (beam FWHM ~0.15") 1.2-mm continuum and line observations towards G16.59$-$0.05 with the Atacama Large Millimeter Array (ALMA). The main dust clump, with size ~10$^4$ au, is resolved into four relatively compact (diameter ~2000 au) millimeter (mm) sources. The source harboring the high-mass YSO is the most prominent in molecular emission. By fitting the emission profiles of several unblended and optically thin transitions of CH$_3$OCH$_3$ and CH$_3$OH, we derived gas temperatures inside the mm-sources in the range 42--131 K, and calculated masses of 1--5 $M_{\odot}$. A well-defined Local Standard of Rest velocity (Vlsr) gradient is detected in most of the high-density molecular tracers at the position of the high-mass YSO, pinpointed by compact 22-GHz free-free emission. This gradient is oriented along a direction forming a large (~70 degree) angle with the radio jet, traced by elongated 13-GHz continuum emission. The butterfly-like shapes of the P-V plots and the linear pattern of the emission peaks of the molecular lines at high velocity confirm that this Vlsr gradient is due to rotation of the gas in the disk surrounding the high-mass YSO. The disk radius is ~500 au, and the Vlsr distribution along the major axis of the disk is well reproduced by a Keplerian profile around a central mass of 10$\pm$2 $M_{\odot}$. The position of the YSO is offset by >~ 0.1" from the axis of the radio jet and the dust emission peak. To explain this displacement we argue that the high-mass YSO could have moved from the center of the parental mm source owing to dynamical interaction with one or more companions.

Orion KL is one of the prime templates of astrochemical and prebiotic chemical studies. We wish to explore more organic molecules with increasing complexity in this region. In particular, we have searched for one of the most complex organic molecules detected in space so far, ethyl formate (C$_{2}$H$_{5}$OCHO). This species is the next step in chemical complexity after the simplest member of esters (methyl formate, CH$_{3}$OCHO). The mechanisms leading to its formation are still poorly known. We have used high angular resolution ($\sim$ 1.$^{\prime\prime}$5) ALMA observations covering a large bandwidth from 214 to 247 GHz. We have detected 82 unblended lines of C$_{2}$H$_{5}$OCHO (49 and 33 of the trans and gauche conformers, respectively). The line images showed that C$_{2}$H$_{5}$OCHO arises mainly from the compact ridge and the hot core-southwest regions. The derived rotational temperatures and column densities are 122 $\pm$ 34 K, (0.9 $\pm$ 0.3) $\times$ 10$^{16}$ cm$^{-2}$ for the hot core-SW, and 103 $\pm$ 13 K, (0.6 $\pm$ 0.3) $\times$ 10$^{16}$ cm$^{-2}$ for the compact ridge. The comparison of spatial distribution and abundance ratios with chemically related molecules (methyl formate, ethanol and formic acid) indicates that C$_{2}$H$_{5}$OCHO is likely formed on the surface of dust grains by addition of CH$_{3}$ to functional-group radicals (CH$_{2}$OCHO) derived from methyl formate (CH$_{3}$OCHO).

We report the first detection in the interstellar medium of the Z-isomer of cyanomethanimine (HNCHCN), an HCN dimer proposed as precursor of adenine. We identified six transitions of Z-cyanomethanimine, along with five transitions of E-cyanomethanimine, using IRAM 30m observations towards the Galactic Center quiescent molecular cloud G+0.693. The Z-isomer has a column density of (2.0$\pm$0.6)$\times$10$^{14}$ cm$^{-2}$ and an abundance of 1.5$\times$10$^{-9}$. The relative abundance ratio between the isomers is [Z/E]$\sim$6. This value cannot be explained by the two chemical formation routes previously proposed (gas-phase and grain surface), which predicts abundances ratios between 0.9 and 1.5. The observed [Z/E] ratio is in good agreement with thermodynamic equilibrium at the gas kinetic temperature (130$-$210 K). Since isomerization is not possible in the ISM, the two species may be formed at high temperature. New chemical models, including surface chemistry on dust grains and gas-phase reactions, should be explored to explain our findings. Whatever the formation mechanism, the high abundance of Z-HNCHCN shows that precursors of adenine are efficiently formed in the ISM.

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.

We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources.

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).

Phosphorus is a crucial element in prebiotic chemistry, especially the P$-$O bond, which is key for the formation of the backbone of the deoxyribonucleic acid. So far, PO had only been detected towards the envelope of evolved stars, and never towards star-forming regions. We report the first detection of PO towards two massive star-forming regions, W51 e1/e2 and W3(OH), using data from the IRAM 30m telescope. PN has also been detected towards the two regions. The abundance ratio PO/PN is 1.8 and 3 for W51 and W3(OH), respectively. Our chemical model indicates that the two molecules are chemically related and are formed via gas-phase ion-molecule and neutral-neutral reactions during the cold collapse. The molecules freeze out onto grains at the end of the collapse and desorb during the warm-up phase once the temperature reaches $\sim$35 K. The observed molecular abundances of 10$^{-10}$ are predicted by the model if a relatively high initial abundance of phosphorus, 5$\times$10$^{-9}$, is assumed.

Hot molecular cores (HMCs), the cradles of massive stars, are the most chemically rich sources in the Galaxy. The typical masses of these cores (few hundreds of solar masses) make them the most important reservoirs of complex organic molecules (COMs), including key species for prebiotic processes. This rich chemistry is thought to be the result of the evaporation of dust grain mantles by the strong radiation of the deeply embedded early-type star(s). Our own Sun may have been born in a high-mass star-forming region, so our Earth may have inherited the primordial chemical composition of its parental hot core region, as suggested by recent studies of oxygen and sulfur chemistry in comets. In this chapter, we discuss how the next generation Very Large Array (ngVLA) can help us to study the emission of heavy COMs in both low- and high-mass star-forming regions. The emission of COMs is important not only because it allows us to understand how chemistry may have developed to eventually form life in our Earth, but also because COMs are a powerful tool for studying the physical properties and kinematics of the dense regions very close to the central protostars.

As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22" (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE-SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up towards the center. Spectral lines towards the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating towards the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core.

The role of accretion disks in the formation of low-mass stars has been well assessed by means of high angular resolution observations at various wavelengths. These findings confirm the prediction that conservation of angular momentum during the collapse leading to the formation of a star is bound to produce flattening and rotation of the collapsing core. What about high-mass stars? At present, several authors have reported on detections of disks around high-mass YSOs. Notwithstanding these important results, the presence of disks rotating about high-mass stars is not sufficient by itself to prove unambiguously the accretion model: what is needed is iron-clad evidence of infall. Such evidence is very difficult to find, as the free-fall velocity becomes significant only very close to the accreting star, i.e., over a region of a few 0.01 pc ($\sim$2000 au), which is very difficult to access and disentangle from the surrounding quiescent or rotating material. In this chapter we discuss how to characterize the infall of material in a sample of 36 high-mass accretion disk candidates covering a broad range of luminosities, from 10$^3$ $L_\odot$ to 10$^6$ $L_\odot$, compiled by Beltr\'an & de Wit (2016) with the next generation Very Large Array (ngVLA).

We present an unbiased spectral line survey toward the Galactic Centre (GC) quiescent giant molecular cloud (QGMC), G+0.693 using the GBT and IRAM 30$\,$ telescopes. Our study highlights an extremely rich organic inventory of abundant amounts of nitrogen (N)-bearing species in a source without signatures of star formation. We report the detection of 17 N-bearing species in this source, of which 8 are complex organic molecules (COMs). A comparison of the derived abundances relative to H$_2$ is made across various galactic and extragalactic environments. We conclude that the unique chemistry in this source is likely to be dominated by low-velocity shocks with X-rays/cosmic rays also playing an important role in the chemistry. Like previous findings obtained for O-bearing molecules, our results for N-bearing species suggest a more efficient hydrogenation of these species on dust grains in G+0.693 than in hot cores in the Galactic disk, as a consequence of the low dust temperatures coupled with energetic processing by X-ray/cosmic ray radiation in the GC.

The fractionation of nitrogen (N) in star-forming regions is a poorly understood process. To put more stringent observational constraints on the N-fractionation, we have observed with the IRAM-30m telescope a large sample of 66 cores in massive star-forming regions. We targeted the (1-0) rotational transition of HN$^{13}$C, HC$^{15}$N, H$^{13}$CN and HC$^{15}$N, and derived the $^{14}$N/$^{15}$N ratio for both HCN and HNC. We have completed this sample with that already observed by Colzi et al. (2018), and thus analysed a total sample of 87 sources. The $^{14}$N/$^{15}$N ratios are distributed around the Proto-Solar Nebula value with a lower limit near the terrestrial atmosphere value ($\sim$272). We have also derived the $^{14}$N/$^{15}$N ratio as a function of the Galactocentric distance and deduced a linear trend based on unprecedented statistics. The Galactocentric dependences that we have found are consistent, in the slope, with past works but we have found a new local $^{14}$N/$^{15}$N value of $\sim$400, i.e. closer to the Prosolar Nebula value. A second analysis was done, and a parabolic Galactocentric trend was found. Comparison with Galactic chemical evolution models shows that the slope until 8 kpc is consistent with the linear analysis, while the flattening trend above 8 kpc is well reproduced by the parabolic analysis.

We present multi-transition observations of PN towards a sample of nine massive dense cores in different evolutionary stages. Using transitions with different excitation conditions, we have found for the first time that the excitation temperatures of PN are in the range 5-30 K. To investigate the main chemical route for the PN formation (surface-chemistry vs. gas-phase chemistry), and the dominant desorption mechanism (thermal vs. shock), we have compared our results with those obtained from molecules tracing different chemical and physical conditions (SiO, SO, CH3OH, and N2H+). We have found that the PN line profiles are very well correlated with those of SiO and SO in six out of the nine targets, which indicates that PN may be released by sputtering of dust grains due to shocks. This finding is corroborated by a faint but statistically significant positive trend between the PN abundance and those of SiO and SO. However, in three objects the PN lines have no hints of high velocity wings, which indicates an alternative origin of PN. Overall, our results indicate that the origin of PN is not unique, as it can be formed in protostellar shocks, but also in colder and more quiescent gas through alternative pathways.

Glycolaldehyde (HOCH$_2$CHO) and ethylene glycol ((CH$_2$OH)$_2$) are two complex organic molecules detected in the hot cores and hot corinos of several star-forming regions. The ethylene glycol/glycolaldehyde abundance ratio seems to show an increase with the source luminosity. In the literature, several surface-chemistry formation mechanisms have been proposed for these two species. With the UCLCHEM chemical code, we explored the different scenarios and compared the predictions for a range of sources of different luminosities with the observations. None of the scenarios reproduce perfectly the trend. A better agreement is, however, found for a formation through recombination of two HCO radicals followed by successive hydrogenations. The reaction between HCO and CH$_2$OH could also contribute to the formation of glycolaldehyde in addition to the hydrogenation pathway. The predictions are improved when a trend of decreasing H$_2$ density within the core region with T $\geq$ 100 K as a function of luminosity, is included in the model. Destruction reactions of complex organic molecules in the gas phase would also need to be investigated, since they can affect the abundance ratios once the species have desorbed in the warm inner regions of the star-forming regions.

Phosphorus (P) is one of the essential elements for life due to its central role in biochemical processes. Recent searches have shown that P-bearing molecules (in particular PN and PO) are present in star-forming regions, although their formation routes remain poorly understood. In this Letter, we report observations of PN and PO towards seven molecular clouds located in the Galactic Center, which are characterized by different types of chemistry. PN is detected in five out of seven sources, whose chemistry is thought to be shock-dominated. The two sources with PN non-detections correspond to clouds exposed to intense UV/X-rays/cosmic-ray radiation. PO is detected only towards the cloud G+0.693$-$0.03, with a PO/PN abundance ratio of $\sim$1.5. We conclude that P-bearing molecules likely form in shocked gas as a result of dust grain sputtering, while are destroyed by intense UV/X-ray/cosmic ray radiation.

Young stellar objects are known to exhibit strong radio variability on timescales of weeks to months, and a few reports have documented extreme radio flares, with at least an order of magnitude change in flux density on timescales of hours to days. However, there have been few constraints on the occurrence rate of such radio flares or on the correlation with pre-main sequence X-ray flares, although such correlations are known for the Sun and nearby active stars. Here we report simultaneous deep VLA radio and Chandra X-ray observations of the Orion Nebula Cluster, targeting hundreds of sources to look for the occurrence rate of extreme radio variability and potential correlation with the most extreme X-ray variability. We identify 13 radio sources with extreme radio variability, with some showing an order of magnitude change in flux density in less than 30 minutes. All of these sources show X-ray emission and variability, but only on timescales <1h do we find clear correlations with extreme radio flaring. Strong X-ray variability does not predict the extreme radio sources and vice versa. Radio flares thus provide us with a new perspective on high-energy processes in YSOs and the irradiation of their protoplanetary disks. Finally, our results highlight implications for interferometric imaging of sources violating the constant-sky assumption.

We report the detection of the prebiotic molecule CH3NCO in a solar-type protostar, IRAS16293-2422 B. A significant abundance of this species on the surface of the comet 67P/Churyumov-Gerasimenko has been proposed, and it has recently been detected in hot cores around high-mass protostars. We observed IRAS16293-2422 B with ALMA in the 90 GHz to 265 GHz range, and detected 8 unblended transitions of CH3NCO. From our Local Thermodynamic Equilibrium analysis we derived an excitation temperature of 110+-19 K and a column density of (4.0+-0.3)x10^15 cm^-2 , which results in an abundance of <=(1.4+-0.1)x10^-10 with respect to molecular hydrogen. This implies a CH3NCO/HNCO and CH3NCO/NH2CHO column density ratios of ~0.08. Our modelling of the chemistry of CH3NCO suggests that both ice surface and gas phase formation reactions of this molecule are needed to explain the observations.

The detection of organic molecules with increasing complexity and potential biological relevance is opening the possibility to understand the formation of the building blocks of life in the interstellar medium. One of the families of molecules with astrobiological interest are the esters, whose simplest member, methyl formate, is rather abundant in star-forming regions. The next step in the chemical complexity of esters is ethyl formate, C$_2$H$_5$OCHO. Only two detections of this species have been reported so far, which strongly limits our understanding of how complex molecules are formed in the interstellar medium. We have searched for ethyl formate towards the W51 e2 hot molecular core, one of the most chemically rich sources in the Galaxy and one of the most promising regions to study prebiotic chemistry, especially after the recent discovery of the P$-$O bond, key in the formation of DNA. We have analyzed a spectral line survey towards the W51 e2 hot molecular core, which covers 44 GHz in the 1, 2 and 3 mm bands, carried out with the IRAM 30m telescope. We report the detection of the trans and gauche conformers of ethyl formate. A Local Thermodynamic Equilibrium analysis indicates that the excitation temperature is 78$\pm$10 K and that the two conformers have similar source-averaged column densities of (2.0$\pm$0.3)$\times$10$^{16}$ cm$^{-2}$ and an abundance of $\sim$10$^{-8}$. We compare the observed molecular abundances of ethyl formate with different competing chemical models based on grain surface and gas-phase chemistry. We propose that grain-surface chemistry may have a dominant role in the formation of ethyl formate (and other complex organic molecules) in hot molecular cores, rather than reactions in the gas phase.

We study the molecular abundance and spatial distribution of the simplest sugar alcohol, ethylene glycol (EG), the simplest sugar glycoladehyde (GA), and other chemically related complex organic species towards the massive star-forming region G31.41+0.31. We have analyzed multiple single dish and interferometric data, and obtained excitation temperatures and column densities using an LTE analysis. We have reported for the first time the presence of EG towards G31.41+0.31, and we have also detected multiple transitions of other complex organic molecules such as GA, methyl formate (MF), dimethyl ether (DME) and ethanol (ET). The high angular resolution images show that the EG emission is very compact, peaking towards the maximum of the continuum. These observations suggest that low abundance complex organic molecules, like EG or GA, are good probes of the gas located closer to the forming stars. Our analysis confirms that EG is more abundant than GA in G31.41+0.31, as previously observed in other interstellar regions. Comparing different star-forming regions we find evidence of an increase of the EG/GA abundance ratio with the luminosity of the source. The DME/MF and EG/ET ratios are nearly constant with luminosity. We have also found that the abundance ratios of pairs of isomers GA/MF and ET/DME decrease with the luminosity of the sources. The most likely explanation for the behavior of the EG/GA ratio is that these molecules are formed by different chemical formation routes not directly linked; although warm-up timescales effects and different formation and destruction efficiencies in the gas phase cannot be ruled out. The most likely formation route of EG is by combination of two CH$_{2}$OH radicals on dust grains. We also favor that GA is formed via the solid-phase dimerization of the formyl radical HCO, and a chemical link between MF and DME.

Phosphorus is a crucial element in biochemistry, especially the P-O bond, which is key for the formation of the backbone of the deoxyribonucleic acid. So far, PO has only been detected towards the envelope of evolved stars, and never towards star-forming regions. We report the first detections of PO towards two massive star-forming regions, W51 e1/e2 and W3(OH), using data from the IRAM 30m telescope. PN has also been detected towards the two regions. The abundance ratio PO/PN is 1.8 and 3 for W51 and W3(OH), respectively. Our chemical model indicates that the two molecules are chemically related and are formed via gas-phase ion-molecule and neutral-neutral reactions during the cold collapse. The molecules freeze out onto grains at the end of the collapse and desorb during the warm-up phase once the temperature reaches 35 K. Similar abundances of the two species are expected during a period of 5x10^{4} yr at the early stages of the warm-up phase, when the temperature is in the range 35-90 K. The observed molecular abundances of 10^{-10} are predicted by the model if a relatively high initial abundance of 5x10^{-9} of depleted phosphorus is assumed.

Phosphorus is a crucial element for the development of life, but so far P-bearing molecules have been detected only in a few astrophysical objects, hence its interstellar chemistry is almost totally unknown. Here we show new detections of phosphorus nitride in a sample of dense cores in different evolutionary stages of the intermediate- and high-mass star formation process: starless, with protostellar objects, and with ultracompact HII regions. All detected PN line widths are smaller than ~5 km/s , and they arise from regions associated with kinetic temperatures smaller than 100 K. Because the few previous detections reported in the literature are associated with warmer and more turbulent sources, the results of this work show that PN can arise from relatively quiescent and cold gas. This information is challenging for theoretical models that invoke either high desorption temperatures or grain sputtering from shocks to release phosphorus into the gas phase. Derived column densities are of the order of 10^{11-12} cm^{-2}, marginally lower than the values derived in the few high-mass star forming regions detected so far. New constraints on the column density upper limit of phosphorus monoxide, the fundamental unit of biologically relevant molecules, are also given.

We present a deep centimeter-wavelength catalog of the Orion Nebula Cluster (ONC), based on a 30h single-pointing observation with the Karl G. Jansky Very Large Array in its high-resolution A-configuration using two 1 GHz bands centered at 4.7 GHz and 7.3 GHz. A total of 556 compact sources were detected in a map with a nominal rms noise of 3 muJy/bm, limited by complex source structure and the primary beam response. Compared to previous catalogs, our detections increase the sample of known compact radio sources in the ONC by more than a factor of seven. The new data show complex emission on a wide range of spatial scales. Following a preliminary correction for the wideband primary-beam response, we determine radio spectral indices for 170 sources whose index uncertainties are less than +/-0.5. We compare the radio to the X-ray and near-infrared point-source populations, noting similarities and differences.

The Premiale Project "Science and Technology in Italy for the upgraded ALMA Observatory - iALMA" has the goal of strengthening the scientific, technological and industrial Italian contribution to the Atacama Large Millimeter/submillimeter Array (ALMA), the largest ground based international infrastructure for the study of the Universe in the microwave. One of the main objectives of the Science Working Group (SWG) inside iALMA, the Work Package 1, is to develop the Italian contribution to the Science Case for the ALMA Band 2 or Band 2+3 receiver. ALMA Band 2 receiver spans from ~67 GHz (bounded by an opaque line complex of ozone lines) up to 90 GHz which overlaps with the lower frequency end of ALMA Band 3. Receiver technology has advanced since the original definition of the ALMA frequency bands. It is now feasible to produce a single receiver which could cover the whole frequency range from 67 GHz to 116 GHz, encompassing Band 2 and Band 3 in a single receiver cartridge, a so called Band 2+3 system. In addition, upgrades of the ALMA system are now foreseen that should double the bandwidth to 16 GHz. The science drivers discussed below therefore also discuss the advantages of these two enhancements over the originally foreseen Band 2 system.

We used the Very Large Array to carry out a multi-epoch radio continuum monitoring of the Orion Nebula Cluster and Orion Molecular Cloud. Our observations reveal the presence of 19 sources. With the exception of the sources BN and C the sources show variability between the different epochs. We have found tentative evidence of variability in the massive object related with source I. Our observations also confirm radio flux density variations of a factor >2 on timescales of hours to days in 5 sources. One of these flaring sources, OHC-E, has been detected for the first time. We conclude that the radio emission arises from: i) highly-variable non-thermal gyrosynchrotron emission produced by electrons accelerated in the magnetospheres of pre-main sequence stars; ii) thermal emission from ionized gas and/or heated dust around massive objects and proplyds. Combining our sample with other radio monitoring and a X-ray catalog, we have studied the properties of 51 radio/X-ray stars. We have found severals hints of a direct relation between the X-ray activity and the mechanisms responsible for (at least some fraction of) the radio emission. We have estimated a radio flaring rate of 0.14 flares day-1 in the densest stellar cluster of the region, suggesting that radio flares from young stars are more common events than previously thought.

We present the first high-angular resolution study of the MonR2 star-forming complex carried out with the Submillimeter Array at (sub-)millimeter wavelengths. We image the continuum and molecular line emission toward the young stellar objects in MonR2 at 0.85mm and 1.3mm, with resolutions ranging from 0.5" to ~3". While free-free emission dominates the IRS1 and IRS2 continuum, dust thermal emission prevails for IRS3 and IRS5, giving envelope masses of ~0.1-0.3 M_Sun. IRS5 splits into at least two sub-arcsecond scale sources, IRS5B and the more massive IRS5A. Our 12CO(2-1) images reveal 11 previously unknown molecular outflows in the MonR2 clump. Comparing these outflows with known IR sources in the IRS5 and IRS3 subclusters allows for tentative identification of driving stars. Line images of molecular species such as CH3CN or CH3OH show that, besides IRS3 (a well-known hot molecular core), IRS5 is also a chemically active source in the region. The gas excitation temperature derived from CH3CN lines toward IRS5 is 144 \pm 15 K, indicating a deeply embedded protostar at the hot-core evolutionary stage. SED fitting of IRS5 gives a mass of ~7 M_Sun and a luminosity of 300 L_Sun for the central source. The derived physical properties of the CO outflows suggest that they contribute to the turbulent support of the MonR2 complex and to the gas velocity dispersion in the clump's center. The detection of a large number of CO outflows widespread across the region supports the competitive accretion scenario as origin of the MonR2 star cluster.

Formation of stars is now believed to be tightly linked to the dynamical evolution of interstellar filaments in which they form. In this paper we analyze the density structure and kinematics of a small network of infrared dark filaments, SDC13, observed in both dust continuum and molecular line emission with the IRAM 30m telescope. These observations reveal the presence of 18 compact sources amongst which the two most massive, MM1 and MM2, are located at the intersection point of the parsec-long filaments. The dense gas velocity and velocity dispersion observed along these filaments show smooth, strongly correlated, gradients. We discuss the origin of the SDC13 velocity field in the context of filament longitudinal collapse. We show that the collapse timescale of the SDC13 filaments (from 1 Myr to 4 Myr depending on the model parameters) is consistent with the presence of Class I sources in them, and argue that, on top of bringing more material to the centre of the system, collapse could generate additional kinematic support against local fragmentation, helping the formation of starless super-Jeans cores.

We have studied the young low-mass pre-main sequence (PMS) stellar population associated with the massive star-forming region DR 21 by using archival X-ray Chandra observations and by complementing them with existing optical and IR surveys. The Chandra observations have revealed for the first time a new highly extincted population of PMS low-mass stars previously missed in observations at other wavelengths. The X-ray population exhibits three main stellar density peaks, coincident with the massive star-forming regions, being the DR 21 core the main peak. The cross-correlated X-ray/IR sample exhibits a radial "Spokes-like" stellar filamentary structure that extends from the DR 21 core towards the northeast. The near IR data reveal a centrally peaked structure for the extinction, which exhibits its maximum in the DR 21 core and gradually decreases with the distance to the N-S cloud axis and to the cluster center. We find evidence of a global mass segregation in the full low-mass stellar cluster, and of an stellar age segregation, with the youngest stars still embedded in the N-S cloud, and more evolved stars more spatially distributed. The results are consistent with the scenario where an elongated overall potential well created by the full low-mass stellar cluster funnels gas through filaments feeding stellar formation. Besides the full gravitational well, smaller-scale local potential wells created by dense stellar sub-clusters of low-mass stars are privileged in the competition for the gas of the common reservoir, allowing the formation of massive stars. We also discuss the possibility that a stellar collision in the very dense stellar cluster revealed by Chandra in the DR 21 core is the origin of the large-scale and highly-energetic outflow arising from this region.

To distinguish between the different theories proposed to explain massive star formation, it is crucial to establish the distribution, the extinction, and the density of low-mass stars in massive star-forming regions. We analyzed deep X-ray observations of the Orion massive star-forming region using the Chandra Orion Ultradeep Project (COUP) catalog. We found that pre-main sequence (PMS) low-mass stars cluster toward the three massive star-forming regions: the Trapezium Cluster (TC), the Orion Hot Core (OHC), and OMC1-S. We derived low-mass stellar densities of 10^{5} stars pc^{-3} in the TC and OMC1-S, and of 10^{6} stars pc^{-3} in the OHC. The close association between the low-mass star clusters with massive star cradles supports the role of these clusters in the formation of massive stars. The X-ray observations show for the first time in the TC that low-mass stars with intermediate extinction are clustered toward the position of the most massive star, which is surrounded by a ring of non-extincted low-mass stars. Our analysis suggests that at least two basic ingredients are needed in massive star formation: the presence of dense gas and a cluster of low-mass stars. The scenario that better explains our findings assumes high fragmentation in the parental core, accretion at subcore scales that forms a low-mass stellar cluster, and subsequent competitive accretion.

Outflows arising from very young stars affect their surroundings and influence the star formation in the parental core. Multiple molecular outflows and Herbig-Haro (HH) objects have been observed in Orion, many of them originating from the embedded massive star-forming region known as OMC1-S. The detection of the outflow driving sources is commonly difficult, because they are still hidden behind large extinction, preventing their direct observation at optical and even near and mid-IR wavelengths. With the aim of improving the identification of the driving sources of the multiple outflows detected in OMC1-S, we used the catalog provided by deep X-ray observations, which have unveiled the very embedded population of pre-main sequence stars. We compared the position of stars observed by the Chandra Orion Ultra Deep project (COUP) in OMC1-S with the morphology of the molecular outflows and the directions of measured proper motions of HH optical objects. We find that 6 out of 7 molecular outflows reported in OMC1-S (detection rate of 86 %) have an extincted X-ray COUP star located at the expected position of the driving source. In several cases, X-rays detected the possible driving sources for the first time. This clustered embedded population revealed by Chandra is very young, with an estimated average age of few 10^{5} yr. It is also likely responsible for the multiple HH objects, which are the optical correspondence of flows arising from the cloud. We show that the molecular outflows driven by the members of the OMC1-S cluster can account for the observed turbulence at core-scales and regulate the star formation efficiency. We discuss the effects of outflow feedback in the formation of massive stars, concluding that the injected turbulence in OMC1-S is compatible with a competitive accretion scenario.

To distinguish between the different theories proposed to explain massive star formation, it is crucial to establish the distribution, the extinction, and the density of low-mass stars in massive star-forming regions. We analyze deep X-ray observations of the Orion massive star-forming region using the Chandra Orion Ultradeep Project (COUP) catalog. We studied the stellar distribution as a function of extinction, with cells of 0.03 pc x 0.03 pc, the typical size of protostellar cores. We derived stellar density maps and calculated cluster stellar densities. We found that low-mass stars cluster toward the three massive star-forming regions: the Trapezium Cluster (TC), the Orion Hot Core (OHC), and OMC1-S. We derived low-mass stellar densities of 10^{5} stars pc^{-3} in the TC and OMC1-S, and of 10^{6} stars pc^{-3} in the OHC. The close association between the low-mass star clusters with massive star cradles supports the role of these clusters in the formation of massive stars. The X-ray observations show for the first time in the TC that low-mass stars with intermediate extinction are clustered toward the position of the most massive star, which is surrounded by a ring of non-extincted low-mass stars. This 'envelope-core' structure is also supported by infrared and optical observations. Our analysis suggests that at least two basic ingredients are needed in massive star formation: the presence of dense gas and a cluster of low-mass stars. The scenario that better explains our findings assumes high fragmentation in the parental core, accretion at subcore scales that forms a low-mass stellar cluster, and subsequent competitive accretion. Finally, although coalescence does not seem a common mechanism for building up massive stars, we show that a single stellar merger may have occurred in the evolution of the OHC cluster, favored by the presence of disks, binaries, and gas accretion.

We report the detection of a new radio recombination line (RRL) maser object toward the IRS2 source in the MonR2 ultracompact HII region. The continuum emission at 1.3mm and 0.85mm and the H30a and H26a lines were observed with the Submillimeter Array (SMA) at angular resolutions of about 0.5"-3". The SMA observations show that the MonR2-IRS2 source is very compact and remains unresolved at spatial scales <=400AU. Its continuum power spectrum at millimeter wavelengths is almost flat (alpha=-0.16, with S_nu proportional to nu^alpha), indicating that this source is dominated by optically thin free-free emission. The H30a and H26a RRL emission is also compact and peaks toward the position of the MonR2-IRS2 source. The measured RRL profiles are double-peaked with the H26a line showing a clear asymmetry in its spectrum. Since the derived line-to-continuum flux ratios (80 and 180kms-1 for H30a and H26a, respectively) exceed the LTE predictions, the RRLs toward MonR2-IRS2 are affected by maser amplification. The amplification factors are however smaller than those found toward the emission line star MWC349A, indicating that MonR2-IRS2 is a weakly amplified maser. Radiative transfer modelling of the RRL emission toward this source shows that the RRL masers arise from a dense and collimated jet embedded in a cylindrical ionized wind, oriented nearly along the direction of the line-of-sight. High-angular resolution observations at sub-millimeter wavelengths are needed to unveil weakly amplified RRL masers in very young massive stars.