Search

ALMA-IMF imaged 15 massive protoclusters down to a resolution of of 2 kau scales, identifying about 1000 star-forming cores. The mass and luminosity of these cores, which are fundamental physical characteristics, are difficult to determine, a problem greatly exacerbated at the distances >2 kpc of ALMA-IMF protoclusters. We combined new datasets and radiative transfer modeling to characterize these cores. We estimated their mass-averaged temperature and the masses these estimates imply. For 1/6 of the sample, we measured the bolometric luminosities, implementing deblending corrections when necessary. We used spectral energy distribution (SED) analysis obtained with the PPMAP Bayesian procedure, which aims to preserve the best angular resolution of the input data. We extrapolated the luminosity and dust temperature images provided by PPMAP at 2.5" resolution to estimate those of individual cores, which were identified at higher angular resolution. To do this, we applied approximate radiative transfer relationships between the luminosity of a protostar and the temperature of its surrounding envelope and between the external heating of prestellar cores and their temperatures. For the first time, we provide data-informed estimates of dust temperatures for 883 cores identified with ALMA-IMF: 17-31 K and 28-79 K (5th and 95th percentiles, up to 127 K) for the 617 prestellar and 266 protostellar cores, respectively. We also measured protostellar luminosities spanning 20-80 000 Lsun. For hot cores, we estimated systematically lower temperatures than studies based on complex organic molecules. We established a mass-luminosity evolutionary diagram, for the first time at the core spatial resolution and for a large sample of high-mass protostellar cores. The ALMA-IMF data favor a scenario in which protostars accrete their mass from a larger mass reservoir than their host cores., Comment: 42 pages, 13 figures, accepted by A&A

ALMA-IMF observed 15 massive protoclusters capturing multiple spectral lines and the continuum emission. We focus on the G351.77 protocluster ($\sim$ 2500 M$_{\odot}$, estimated from single-dish continuum observations) located at 2 kpc. We trace the dense gas emission and kinematics with N$_2$H$^+$ (1-0) at $\sim$ 4 kau resolution. We estimate an N$_2$H$^+$ relative abundance $\sim (1.7 \pm 0.5) \times 10^{-10}$. We decompose the N$_2$H$^+$ emission into up to two velocity components, highlighting the kinematic complexity in the dense gas. By examining the position-velocity (PV) diagrams on small scales, we observe clear inflow signatures (V-shapes) associated with dense cores. The most prominent V-shape has a mass inflow rate of $\sim 9.8 \times 10^{-4}$ M$_{\odot}$ yr$^{-1}$ and a short timescale of $\sim$ 15.6 kyr. We also observe V-shapes without associated cores. This suggests both that cores or centers of accretion exist below the 1.3 mm detection limit, and that the V-shapes may be viable tracers of very early accretion and star formation on $\sim$ 4 kau scales. The large-scale PV diagram shows that the protocluster is separated into 2 principal velocity structures. Combined with smaller scale DCN and H$_2$CO emission, we propose a scenario of larger scale slow contraction with rotation in the center based on simple toy models. This scenario leads the suggestion of outside-in evolution of the protocluster as it collapses. The gas depletion times implied by the V-shapes are short ($\sim$ 0.3 Myr), requiring either very fast cluster formation, and/or continuous mass feeding of the protocluster. The latter is possible via the Mother Filament G351.77 is forming out of. The similarities in the properties of G351.77 and the recently published work in G353.41 indicate that many of the physical conditions inferred via the ALMA-IMF N$_2$H$^+$ observations may be generic to protoclusters., Comment: Submitted in A&A, 28 pages, 31 figures, 4 interactive figures, 7 tables

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

A crucial aspect in addressing the challenge of measuring the core mass function, that is pivotal for comprehending the origin of the initial mass function, lies in constraining the temperatures of the cores. We aim to measure the luminosity, mass, column density and dust temperature of star-forming regions imaged by the ALMA-IMF large program. High angular resolution mapping is required to capture the properties of protostellar and pre-stellar cores and to effectively separate them from larger features, such as dusty filaments. We employed the point process mapping (PPMAP) technique, enabling us to perform spectral energy distribution fitting of far-infrared and submillimeter observations across the 15 ALMA-IMF fields, at an unmatched 2.5" angular resolution. By combining the modified blackbody model with near-infrared data, we derived bolometric luminosity maps. We estimated the errors impacting values of each pixel in the temperature, column density, and luminosity maps. Subsequently, we employed the extraction algorithm getsf on the luminosity maps in order to detect luminosity peaks and measure their associated masses. We obtained high-resolution constraints on the luminosity, dust temperature, and mass of protoclusters, that are in agreement with previously reported measurements made at a coarser angular resolution. We find that the luminosity-to-mass ratio correlates with the evolutionary stage of the studied regions, albeit with intra-region variability. We compiled a PPMAP source catalog of 313 luminosity peaks using getsf on the derived bolometric luminosity maps. The PPMAP source catalog provides constraints on the mass and luminosity of protostars and cores, although one source may encompass several objects. Finally, we compare the estimated luminosity-to-mass ratio of PPMAP sources with evolutionary tracks and discuss the limitations imposed by the 2.5" beam., Comment: 35 pages, 19 figures, 5 tables Accepted by A&A

The ALMA-IMF Large Program provides multi-tracer observations of 15 Galactic massive protoclusters at matched sensitivity and spatial resolution. We focus on the dense gas kinematics of the G353.41 protocluster traced by N$_2$H$^+$ (1$-$0), with an spatial resolution $\sim$0.02 pc. G353.41, at a distance of $\sim$2 kpc, has a mass of $\sim$2500 M$_{\odot}$ within $1.3\times1.3$ pc$^2$. We extract the N$_2$H$^+$ isolated line component and we decompose it by fitting up to 3 Gaussian velocity components. This allows us to identify velocity structures that are impossible to identify in the traditional position-velocity diagram. We identify multiple velocity gradients on large ($\sim$1 pc) and small scales ($\sim$0.2 pc). We find good agreement between the N$_2$H$^+$ velocities and the previously reported DCN core velocities, suggesting that cores are kinematically coupled to the dense gas in which they form. We measure 9 converging ``V-shaped'' velocity gradients ($\sim20$ km/s/pc), located in filaments, which are sometimes associated with cores near their point of convergence. The average inflow timescale is $\sim67$ kyr, or about twice the free-fall time of cores in the same area ($\sim33$ kyr) but substantially shorter than protostar lifetime estimates ($\sim$0.5 Myr). We derive mass accretion rates in the range of $(0.35-8.77)\,\times\,10^{-4}$ M$_{\odot}$/yr. This feeding might lead to further filament collapse and formation of new cores. We suggest that the protocluster is collapsing on large scales, but the velocity signature of collapse is slow compared to pure free-fall. These data are consistent with a comparatively slow global protocluster contraction under gravity, and faster core formation within, suggesting the formation of multiple generations of stars over the protocluster lifetime., Comment: Accepted in A&A, 22 pages, 21 figures, 1 interactive figure, 4 tables

We aim to identify and characterize cores in the high-mass proto-cluster W49, determine their evolutionary stages and measure the associated lifetimes. We built a catalog of 129 cores extracted from an ALMA 1.3 mm continuum image at 0.26" (2900 au) angular resolution. The association between cores and Hyper/Ultra Compact HII (H/UC HII) regions was established from the analysis of VLA 3.3 cm continuum and H30$\alpha$ line observations. We also looked for emission of hot molecular cores (HMCs) using the methyl formate doublet at 218.29 GHz. We identified 40 cores associated with an H/UC HII region and 19 HMCs over the ALMA mosaic. The 52 cores with an H/UC HII region and/or a HMC are assumed to be high-mass protostellar cores, while the rest of the core population likely consists in prestellar cores and low-mass protostellar cores. We found a good agreement between the two tracers of ionized gas, with 23 common detections and only four cores detected at 3.3 cm and not in H30$\alpha$. The spectral indexes from 3.3 cm to 1.3 mm range from 1, for the youngest cores with partially optically thick free-free emission, to about -0.1, that is optically thin free-free emission obtained for cores likely more evolved. Using as a reference the H/UC HII regions, we found the statistical lifetimes of the HMC and massive protostellar phases in W49N to be about $6\times10^4$ yr and $1.4\times10^5$ yr, respectively. We also showed that HMC can co-exist with H/UC HII regions during a short fraction of the core lifetime, about $2\times10^4$ yr. This indicates a rapid dispersal of the inner molecule envelope once the HC HII is formed., Comment: 15 pages, 7 figures. Accepted for publication in A&A

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

We aimed to measure the CMF in the evolved W33-Main star-forming protocluster to compare it with CMF recently obtained in other Galactic star-forming regions, including the ones included in the ALMA-IMF program. We used observations from the ALMA-IMF large program: 2'x2' maps of emission from the continuum and selected lines at 1.3mm and 3mm observed by the ALMA 12m only antennas. Our angular resolution was typically 1'', that is 2400au at a distance of 2.4kpc. The lines we analysed are CO(2-1), SiO(5-4), N2H+(1-0), H41alpha as well as He41alpha blended with C41alpha. We built a census of dense cores in the region, and we measured the associated CMF based on a core-dependent temperature value. We confirmed the 'evolved' status of W33-Main by identifiying three HII regions within the field, and to a lesser extent based on the number and extension of N2H+ filaments. We produced a filtered core catalog of 94 candidates, that we refined to take into account the contamination of the continuum by free-free and line emission, obtaining 80 cores with masses that range from 0.03 to 13.2Msun. We fitted the resulting high-mass end of the CMF with a single power law of the form N(log(M)) ~ M^alpha, obtaining alpha = -1.44(+0.16)(-0.22), slightly steeper but consistent with the Salpeter index. We categorized our cores in pre- and protostellar, mostly based on outlow activity and hot core nature. We found the prestellar CMF to be steeper than a Salpeter-like distribution, and the protostellar CMF to be slightly top heavy. We found a higher proportion of cores within the HII regions and their surroundings than in the rest of the field. We also found that the cores' masses were rather low (maximum mass of 13Msun).

We present a catalog of 315 protostellar outflow candidates detected in SiO J=5-4 in the ALMA-IMF Large Program, observed with ~2000 au spatial resolution, 0.339 km/s velocity resolution, and 2-12 mJy/beam (0.18-0.8 K) sensitivity. We find median outflow masses, momenta, and kinetic energies of ~0.3 M$_{\odot}$, 4 M$_{\odot}$ km/s, and 10$^{45}$ erg, respectively. Median outflow lifetimes are 6,000 years, yielding median mass, momentum, and energy rates of $\dot{M}$ = 10$^{-4.4}$ M$_{\odot}$ yr$^{-1}$, $\dot{P}$ = 10$^{-3.2}$ M$_{\odot}$ km/s yr$^{-1}$, and $\dot{E}$ = 1 L$_{\odot}$. We analyze these outflow properties in the aggregate in each field. We find correlations between field-aggregated SiO outflow properties and total mass in cores (~3$-$5$\sigma$), and no correlations above 3$\sigma$ with clump mass, clump luminosity, or clump luminosity-to-mass ratio. We perform a linear regression analysis and find that the correlation between field-aggregated outflow mass and total clump mass - which has been previously described in the literature - may actually be mediated by the relationship between outflow mass and total mass in cores. We also find that the most massive SiO outflow in each field is typically responsible for only 15-30% of the total outflow mass (60% upper limit). Our data agree well with the established mechanical force-bolometric luminosity relationship in the literature, and our data extend this relationship up to L $\geq$ 10$^6$ L$_{\odot}$ and $\dot{P}$ $\geq$ 1 M$_{\odot}$ km/s yr$^{-1}$. Our lack of correlation with clump L/M is inconsistent with models of protocluster formation in which all protostars start forming at the same time., Comment: 46 pages, 14 figures, 10 tables. This publication has an associated Zenodo entry, which can be found here: https://zenodo.org/records/8350595

Using observations with e-MERLIN and the VLA, together with archival data from ALMA, we obtain high-resolution radio images of two binary YSOs: L1551 IRS 5 and L1551 NE, covering a wide range of frequencies from 5 - 336 GHz, and resolving emission from the radio jet on scales of only ~15 au. By comparing these observations to those from a previous epoch, it is shown that there is a high degree of variability in the free-free emission from the jets of these sources. In particular, the northern component of L1551 IRS 5 shows a remarkable decline in flux density of a factor of ~5, suggesting that the free-free emission of this source has almost disappeared. By fitting the spectra of the sources, the ionised mass-loss rates of the jets are derived and it is shown that there is significant variability of up to a factor of ~6 on timescales of ~20 years. Using radiative transfer modelling, we also obtained a model image for the jet of the southern component of L1551 IRS 5 to help study the inner region of the ionised high-density jet. The findings favour the X-wind model launched from a very small innermost region., Comment: 13 pages, 7 figures, accepted for publication in A&A

ALMA-IMF is an Atacama Large Millimeter/submillimeter Array (ALMA) Large Program designed to measure the core mass function (CMF) of 15 protoclusters chosen to span their early evolutionary stages. It further aims to understand their kinematics, chemistry, and the impact of gas inflow, accretion, and dynamics on the CMF. We present here the first release of the ALMA-IMF line data cubes (DR1), produced from the combination of two ALMA 12m-array configurations. The data include 12 spectral windows, with eight at 1.3mm and four at 3mm. The broad spectral coverage of ALMA-IMF (~6.7 GHz bandwidth coverage per field) hosts a wealth of simple atomic, molecular, ionised, and complex organic molecular lines. We describe the line cube calibration done by ALMA and the subsequent calibration and imaging we performed. We discuss our choice of calibration parameters and optimisation of the cleaning parameters, and we demonstrate the utility and necessity of additional processing compared to the ALMA archive pipeline. As a demonstration of the scientific potential of these data, we present a first analysis of the DCN (3-2) line. We find that DCN traces a diversity of morphologies and complex velocity structures, which tend to be more filamentary and widespread in evolved regions and are more compact in the young and intermediate-stage protoclusters. Furthermore, we used the DCN (3-2) emission as a tracer of the gas associated with 595 continuum cores across the 15 protoclusters, providing the first estimates of the core systemic velocities and linewidths within the sample. We find that DCN (3-2) is detected towards a higher percentage of cores in evolved regions than the young and intermediate-stage protoclusters and is likely a more complete tracer of the core population in more evolved protoclusters. The full ALMA 12m-array cubes for the ALMA-IMF Large Program are provided with this DR1 release., Comment: 75 pages (21 main body; 54 appendix), 37 figures. The ALMA-IMF DR1 line release is hosted at https://dataverse.harvard.edu/dataverse/alma-imf

The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars ($M\gtrsim 8M_\odot$) are still not well understood. To this end, we have undertaken a large observational program (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous ($L>10^4L_\odot$) protostellar objects in the 1.37 mm wavelength regime in both continuum and line emission, reaching $\sim$0.4" resolution (800 au at 2 kpc). Using the dense gas tracer CH$_3$CN, we find velocity gradients across 13 cores perpendicular to the directions of bipolar molecular outflows, making them excellent disk candidates. Specific angular momentum ($j$) radial profiles are on average $\sim10^{-3}$ km /s pc and follow $j \propto r^{1.7}$, consistent with a poorly resolved rotating and infalling envelope/disk model. Fitting the velocity profiles with a Keplerian model, we find protostellar masses in the range of $\sim 10-25$ $M_\odot$. Modelling the level population of CH$_3$CN lines, we present temperature maps and find median gas temperatures in the range $70-210$ K. We create Toomre $Q$ maps to study the stability of the disks and find almost all (11 of 13) disk candidates to be prone to fragmentation due to gravitational instabilities at the scales probed by our observations. In particular, disks with masses greater than $\sim10-20\%$ of the mass of their host (proto)stars are Toomre unstable, and more luminous protostellar objects tend to have disks that are more massive and hence more prone to fragmentation. Our finings show that most disks around high-mass protostars are prone to disk fragmentation early in their formation due to their high disk to stellar mass ratio. This impacts the accretion evolution of high-mass protostars which will have significant implications for the formation of the most massive stars., Comment: 27 pages, 12 figures, 6 appendices - accepted for publication in Astronomy and Astrophysics

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

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

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

The ALMA-IMF Large Program observed the W43-MM2-MM3 ridge, whose 1.3mm and 3mm ALMA 12m array continuum images reach a 2500au spatial resolution. We used both the best-sensitivity and the line-free ALMA-IMF images, reduced the noise with the multi-resolution segmentation technique MnGSeg, and derived the most complete and most robust core catalog possible. Using two different extraction software packages, getsf and GExt2D, we identified 200 compact sources, whose 100 common sources have on average fluxes consistent to within 30%. We filtered sources with non-negligible free-free contamination and corrected fluxes from line contamination, resulting in a W43-MM2-MM3 catalog of 205 getsf cores. With a median deconvolved FWHM size of 3400au, core masses range from 0.1Msun to 70Msun and the getsf catalog is 90% complete down to 0.8Msun. The high-mass end of the core mass function (CMF) of W43-MM2-MM3 is top-heavy compared to the canonical IMF. Fitting the cumulative CMF with a single power law of the form N(>logM)\propto M^a, we measured a=-0.95\pm0.04, compared to the canonical a=-1.35 Salpeter IMF slope. The slope of the CMF is robust with respect to map processing, extraction software package, and reasonable variations in the assumptions taken to estimate core masses. We explore several assumptions on how cores transfer their mass to stars and sub-fragment to predict the IMF resulting from the W43-MM2-MM3 CMF. In stark contrast to the commonly accepted paradigm, our result argues against the universality of the CMF shape. More robust functions of the star-formation efficiency and core sub-fragmentation are required to better predict the resulting IMF, here suggested to remain top-heavy at the end of the star-formation phase. If confirmed, the IMFs emerging from starburst events could inherit their top-heavy shape from their parental CMFs, challenging the IMF universality., Comment: 28 pages, 16 figures, 6 tables, A&A accepted on 02/08/22, revised following the language referee report

We present the first data release of the ALMA-IMF Large Program, which covers the 12m-array continuum calibration and imaging. The ALMA-IMF Large Program is a survey of fifteen dense molecular cloud regions spanning a range of evolutionary stages that aims to measure the core mass function (CMF). We describe the data acquisition and calibration done by the Atacama Large Millimeter/submillimeter Array (ALMA) observatory and the subsequent calibration and imaging we performed. The image products are combinations of multiple 12m array configurations created from a selection of the observed bandwidth using multi-term, multi-frequency synthesis imaging and deconvolution. The data products are self-calibrated and exhibit substantial noise improvements over the images produced from the delivered data. We compare different choices of continuum selection, calibration parameters, and image weighting parameters, demonstrating the utility and necessity of our additional processing work. Two variants of continuum selection are used and will be distributed: the "best-sensitivity" data, which include the full bandwidth, including bright emission lines that contaminate the continuum, and "cleanest", which select portions of the spectrum that are unaffected by line emission. We present a preliminary analysis of the spectral indices of the continuum data, showing that the ALMA products are able to clearly distinguish free-free emission from dust emission, and that in some cases we are able to identify optically thick emission sources. The data products are made public with this release., Comment: Data released on Zenodo at https://zenodo.org/record/5702966. Accepted for publication in A&A

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

We aim to characterize the physical and chemical properties of fragmented cores during the earliest evolutionary stages in the very young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953. NOEMA 1.3 mm data are used in combination with archival mid- and far-infrared observations to construct and fit the SEDs of individual fragmented cores. The radial density profiles are inferred from the 1.3 mm continuum visibility profiles and the radial temperature profiles are estimated from H2CO rotation temperature maps. Molecular column densities are derived with the line fitting tool XCLASS. The physical and chemical properties are combined by applying the physical-chemical model MUSCLE in order to constrain the chemical timescales of a few line-rich cores. The morphology and spatial correlations of the molecular emission are analyzed using the HOG method. The mid-infrared data show that both regions contain a cluster of young stellar objects. Bipolar molecular outflows are observed in the CO 2-1 transition toward the strong mm cores indicating protostellar activity. We find strong molecular emission of SO, SiO, H2CO, and CH3OH in locations which are not associated with the mm cores. These shocked knots can be either associated with the bipolar outflows or, in the case of ISOSS J23053+5953, with a colliding flow that creates a large shocked region between the mm cores. The mean chemical timescale of the cores is lower (20 000 yr) compared to that of the sources of the more evolved CORE sample (60 000 yr). With the HOG method, we find that the spatial emission of species tracing the extended emission and of shock-tracing molecules are well correlated within transitions of these groups., Comment: 34 pages, 17 figures, accepted for publication in A&A

We use sub-arcsecond resolution ($\sim$0.4$''$) observations with NOEMA at 1.37 mm to study the dust emission and molecular gas of 18 high-mass star-forming regions. We combine the derived physical and chemical properties of individual cores in these regions to estimate their ages. The temperature structure of these regions are determined by fitting H2CO and CH3CN line emission. The density profiles are inferred from the 1.37 mm continuum visibilities. The column densities of 11 different species are determined by fitting the emission lines with XCLASS. Within the 18 observed regions, we identify 22 individual cores with associated 1.37 mm continuum emission and with a radially decreasing temperature profile. We find an average temperature power-law index of q = 0.4$\pm$0.1 and an average density power-law index of p = 2.0$\pm$0.2 on scales on the order of several 1 000 au. Comparing these results with values of p derived in the literature suggest that the density profiles remain unchanged from clump to core scales. The column densities relative to N(C18O) between pairs of dense gas tracers show tight correlations. We apply the physical-chemical model MUSCLE to the derived column densities of each core and find a mean chemical age of $\sim$60 000 yrs and an age spread of 20 000-100 000 yrs. With this paper we release all data products of the CORE project available at https://www.mpia.de/core. The CORE sample reveals well constrained density and temperature power-law distributions. Furthermore, we characterize a large variety in molecular richness that can be explained by an age spread confirmed by our physical-chemical modeling. The hot molecular cores show the most emission lines, but we also find evolved cores at an evolutionary stage, in which most molecules are destroyed and thus the spectra appear line-poor again., Comment: 29 pages, 14 figures, accepted for publication in A&A

In the massive star-forming region IRAS 21078+5211, a highly fragmented cluster (0.1~pc in size) of molecular cores is observed, located at the density peak of an elongated (1~pc in size) molecular cloud. A small (1~km/s per 0.1~pc) LSR velocity (Vlsr) gradient is detected across the axis of the molecular cloud. Assuming we are observing a mass flow from the harboring cloud to the cluster, we derive a mass infall rate of about 10^{-4}~M_{sun}~yr^{-1}. The most massive cores (labeled 1, 2, and 3) are found at the center of the cluster, and these are the only ones that present a signature of protostellar activity in terms of emission from high-excitation molecular lines or a molecular outflow. We reveal an extended (size about 0.1~pc), bipolar collimated molecular outflow emerging from core 1. We believe this is powered by a (previously discovered) compact (size <= 1000~au) radio jet, ejected by a YSO embedded in core 1 (named YSO-1), since the molecular outflow and the radio jet are almost parallel and have a comparable momentum rate. By means of high-excitation lines, we find a large (14~km/s over 500~au) Vlsr gradient at the position of YSO-1, oriented approximately perpendicular to the radio jet. Assuming this is an edge-on, rotating disk and fitting a Keplerian rotation pattern, we determine the YSO-1 mass to be 5.6+/-2.0~M_{sun}. The water masers (previously observed with VLBI) emerge within 100-300~au from YSO-1 and are unique tracers of the jet kinematics. Their three-dimensional (3D) velocity pattern reveals that the gas flows along, and rotates about, the jet axis. We show that the 3D maser velocities are fully consistent with the magneto-centrifugal disk-wind models predicting a cylindrical rotating jet. Under this hypothesis, we determine the jet radius to be about 16~au and the corresponding launching radius and terminal velocity to be about 2.2~au and 200~km/s, respectively., Comment: 25 pages, 16 figures

We present new Jansky Very Large Array (VLA) images of the central region of the W49A star-forming region at 3.6~cm and at 7~mm at resolutions of 0\farcs15 (1650 au) and 0\farcs04 (440 au), respectively. The 3.6~cm data reveal new morphological detail in the ultracompact \ion{H}{2} region population, as well as several previously unknown and unresolved sources. In particular, source A shows elongated, edge-brightened, bipolar lobes, indicative of a collimated outflow, and source E is resolved into three spherical components. We also present VLA observations of radio recombination lines at 3.6~cm and 7~mm, and IRAM Northern Extended Millimeter Array (NOEMA) observations at 1.2~mm. Three of the smallest ultracompact \ion{H}{2} regions (sources A, B2 and G2) all show broad kinematic linewidths, with $\Delta$V$_{FWHM}\gtrsim$40~km~s$^{-1}$. A multi-line analysis indicates that broad linewidths remain after correcting for pressure broadening effects, suggesting the presence of supersonic flows. Substantial changes in linewidth over the 21 year time baseline at both 3.6 cm and 7 mm are found for source G2. At 3.6 cm, the linewidth of G2 changed from 31.7$\pm$1.8 km s$^{-1}$ to 55.6$\pm$2.7 km s$^{-1}$, an increase of $+$23.9$\pm$3.4 km s$^{-1}$. The G2 source was previously reported to have shown a 3.6~cm continuum flux density decrease of 40\% between 1994 and 2015. This source sits near the center of a very young bipolar outflow whose variability may have produced these changes., Comment: 18 pages, 7 figures, Accepted for publication in The Astronomical Journal

Recent high-angular resolution (40 mas) ALMA observations at 1.14 mm resolve a compact (R~200 au) flattened dust structure perpendicular to the HH 80-81 jet emanating from the GGD 27-MM1 high-mass protostar, making it a robust candidate for a true accretion disk. The jet/disk system (HH 80-81 / GGD 27-MM1) resemble those found in association with low- and intermediate-mass protostars. We present radiative transfer models that fit the 1.14 mm ALMA dust image of this disk which allow us to obtain its physical parameters and predict its density and temperature structure. Our results indicate that this accretion disk is compact (Rdisk~170 au) and massive (5Msun), about 20% of the stellar mass of 20 Msun. We estimate the total dynamical mass of the star-disk system from the molecular line emission finding a range between 21 and 30 Msun, which is consistent with our model. We fit the density and temperature structures found by our model with power law functions. These results suggest that accretion disks around massive stars are more massive and hotter than their low-mass siblings, but they still are quite stable. We also compare the temperature distribution in the GGD 27-MM1 disk with that found in low- and intermediate-mass stars and discuss possible implications on the water snow line. We have also carried about a study of the distance based on Gaia DR2 data and the population of young stellar objects (YSOs) in this region, and from the extinction maps. We conclude that the source distance is in within 1.2 and 1.4 kpc, closer than what was derived in previous studies (1.7 kpc)., Comment: 20 pages, 15 figures

We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observed molecular abundances with the chemical evolution code MUSCLE (MUlti Stage ChemicaL codE). Results. With the kinetic temperature tracers CH3CN and H2CO we observe a temperature distribution with a power-law index of q = 0.41+-0.08. Using the visibilities of the continuum emission we derive a density structure with a power-law index of p = 1.7+-0.1. The hot core spectra reveal high molecular abundances and a rich diversity in complex molecules. The majority of the molecules have an asymmetric spatial distribution around the forming protostar(s), which indicates a complex physical structure on scales < 1 400 au. Using MUSCLE, we are able to explain the observed molecular abundance of 10 out of 14 modeled species at an estimated hot core chemical age of ~21 100 years. In contrast to the observational analysis, our chemical modeling predicts a lower density power-law index of p < 1.4. Reasons for this discrepancy are discussed. Conclusions. Combining high spatial resolution observations with detailed chemical modeling allows us to derive a concise picture of the physical and chemical structure of the famous AFGL 2591 hot core. The next steps are to conduct a similar analysis for the whole CORE sample, and then use this analysis to constrain the chemical diversity in high-mass star formation to a much greater depth., Comment: 22 pages, 13 figures, accepted for publication in A&A

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., Comment: 9 Total, inc references, appendix, long author list, 4 figures (inc appendix)

This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additional single-dish maps, merged with the interferometric data to recover the extended emission. Our correlator setup covered a number of lines from well-known hot core tracers and a few outflow tracers. The angular (~0.45"$-$0.9") and spectral (0.5 km/s) resolutions were sufficient to resolve the clump in IRAS23385+6053 and investigate the existence of large-scale motions due to rotation, infall, or expansion. We find that the clump splits into six distinct cores when observed at sub-arcsecond resolution. These are identified through their 1.4 mm continuum and molecular line emission. We produce maps of the velocity, line width, and rotational temperature from the methanol and methyl cyanide lines, which allow us to investigate the cores and reveal a velocity and temperature gradient in the most massive core. We also find evidence of a bipolar outflow, possibly powered by a low-mass star. We present the tentative detection of a circumstellar self-gravitating disk lying in the most massive core and powering a large-scale outflow previously known in the literature. In our scenario, the star powering the flow is responsible for most of the luminosity of IRAS23385+6053 (~$3000~L_\odot$). The other cores, albeit with masses below the corresponding virial masses, appear to be accreting material from their molecular surroundings and are possibly collapsing or on the verge of collapse. We conclude that we are observing a sample of star-forming cores that is bound to turn into a cluster of massive stars.

In recent years, the disk populations in a number of young star-forming regions have been surveyed with ALMA. Understanding the disk properties and their correlation with those of the central star is critical to understand planet formation. In particular, a decrease of the average measured disk dust mass with the age of the region has been observed. We conducted high-sensitivity continuum ALMA observations of 43 Class II young stellar objects in CrA at 1.3 mm (230 GHz). The typical spatial resolution is 0.3". The continuum fluxes are used to estimate the dust masses of the disks, and a survival analysis is performed to estimate the average dust mass. We also obtained new VLT/X-Shooter spectra for 12 of the objects in our sample. 24 disks are detected, and stringent limits have been put on the average dust mass of the non-detections. Accounting for the upper limits, the average disk mass in CrA is $6\pm3\,\rm M_\oplus$, significantly lower than that of disks in other young (1-3 Myr) star forming regions (e.g. Lupus) and appears consistent with the 5-10 Myr old Upper Sco. The position of the stars in our sample on the HR diagram, however, seems to confirm that that CrA has age similar to Lupus. Neither external photoevaporation nor a lower than usual stellar mass distribution can explain the low disk masses. On the other hand, a low-mass disk population could be explained if the disks are small, which could happen if the parent cloud has a low temperature or intrinsic angular momentum, or if the the angular momentum of the cloud is removed by some physical mechanism such as magnetic braking. In order to fully explain and understand the dust mass distribution of protoplanetary disks and their evolution, it may also be necessary to take into consideration the initial conditions of star and disk formation process, which may vary from region to region, and affect planet formation., Comment: 20 pages, 8 Figures, accepted for publication on A&A. The abstract has been shortened in order to fit arXiv requirements

We used ALMA to observe the star-forming region GGD27 at 1.14 mm with an unprecedented angular resolution, 40 mas (56 au) and sensitivity (0.002 Msun). We detected a cluster of 25 continuum sources, most of which are likely tracing disks around Class 0/I protostars. Excluding the two most massive objects, disks masses are in the range 0.003-0.05 Msun. The analysis of the cluster properties indicates that GGD27 displays moderate subclustering. This result combined with the dynamical timescale of the radio jet (10000 years) suggests the youthfulness of the cluster. The lack of disk mass segregation signatures may support this too. We found a clear paucity of disks with Rdisk >100 au. The median value of the radius is 34 au, smaller than the median of 92 au for Taurus but comparable to the value found in Ophiuchus and in the Orion Nebula Cluster. In GGD27 there is no evidence of a distance-dependent disk mass distribution (i. e., disk mass depletion due to external photoevaporation), most likely due to the cluster youth. There is a clear deficit of disks for distances <0.02 pc. Only for distances >0.04 pc stars can form larger and more massive disks, suggesting that dynamical interactions far from the cluster center are weaker, although the small disks found could be the result of disk truncation. This work demonstrates the potential to characterize disks from low-mass YSOs in distant and massive (still deeply embedded) clustered environments., Comment: Accepted for publication in Astronomy & Astrophysics Letters

Context. Colliding winds in massive binaries are able to accelerate particles up to relativistic speeds as the result of the interaction between the winds of the different stellar components. HD 167971 exhibits this phenomenology which makes it a strong radio source. Aims. We aim at characterizing the morphology of the radio emission and its dependence on the orbital motion, traced independently by NIR-interferometry, of the spectroscopic binary and the tertiary component that conforms HD 167971. Methods. We analyze 2006 and 2016 very long baseline interferometric data at C and X bands. We complement our analysis with a geometrical model of the wind-wind collision region, and with an astrometric description of the system. Results. We confirm that the detected non-thermal radio emission is associated with the wind-wind collision region of the spectroscopic binary and the tertiary component in HD 167971. The wind-wind collision region changes orientation in agreement with the orbital motion of the tertiary around the spectroscopic binary. The total intensity also changes between the two observing epochs in a way inversely proportional to the separation between the SB and T, with a negative-steep spectral index typical of an optically thin synchrotron emission possibly steepened by an inverse Compton cooling effect. The wind-wind collision bow-shock shape and its position with respect to the stars indicates that the wind momentum from the spectroscopic binary is stronger than that of the tertiary. Finally, the astrometric solution derived for the stellar system and the wind-wind collision region is consistent with independent Gaia data., Comment: Accepted for publication in Astronomy and Astrophysics, 7 pages, 6 figures

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., Comment: 23 pages, 13 fig (inc 3 Appendix figures)

The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.37 mm wavelength regime in both line and dust continuum emission. We present the spectral line setup of the CORE survey and a case study for W3(H2O). At ~0.35" (700 AU at 2 kpc) resolution, the W3(H2O) clump fragments into two cores (West and East), separated by ~2300 AU. Velocity shifts of a few km/s are observed in the dense-gas tracer, CH3CN, across both cores, consistent with rotation and perpendicular to the directions of two bipolar outflows, one emanating from each core. The kinematics of the rotating structure about W3(H2O) W shows signs of differential rotation of material, possibly in a disk-like object. The observed rotational signature around W3(H2O) E may be due to a disk-like object, an unresolved binary (or multiple) system, or a combination of both. We fit the emission of CH3CN (12-11) K = 4-6 and derive a gas temperature map with a median temperature of ~165 K across W3(H2O). We create a Toomre Q map to study the stability of the rotating structures against gravitational instability. The rotating structures appear to be Toomre unstable close to their outer boundaries, with a possibility of further fragmentation in the differentially-rotating core W3(H2O) W. Rapid cooling in the Toomre-unstable regions supports the fragmentation scenario. Combining millimeter dust continuum and spectral line data toward the famous high-mass star-forming region W3(H2O), we identify core fragmentation on large scales, and indications for possible disk fragmentation on smaller spatial scales., Comment: 23 pages, 26 figures, accepted for publication in Astronomy and Astrophysics

A number of ultracompact H II regions in Galactic star forming environments have been observed to vary significantly in radio flux density on timescales of 10-20 years. Theory predicted that such variations should occur when the accretion flow that feeds a young massive star becomes unstable and clumpy. We have targeted the massive star-forming region W49A with the Karl G. Jansky Very Large Array (VLA) for observations at 3.6 cm with the B-configuration at 0.8'' resolution, to compare to nearly identical observations taken almost 21 years earlier (February 2015 and August 1994). Most of the sources in the crowded field of ultracompact and hypercompact H II regions exhibit no significant changes over this time period. However, one source, W49A/G2, decreased by 20% in peak intensity (from 71+/-4 mJy/beam to 57+/-3 mJy/beam), and 40% in integrated flux (from 0.109+/-0.011 Jy to 0.067+/-0.007 Jy), where we cite 5 sigma errors in peak intensity, and 10% errors in integrated flux. We present the radio images of the W49A region at the two epochs, the difference image that indicates the location of the flux density decrease, and discuss explanations for the flux density decrease near the position of W49A/G2., Comment: 9 pages, 2 figures, accepted for publication, The Astrophysical Journal Letters, July 26, 2018

Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation., Comment: 29 pages, 17 figures, Astronomy & Astrophysics in press, for a higher-resolution version of the paper see http://www.mpia.de/homes/beuther/papers.html

Here we present deep (16 mumJy), very high (40 mas) angular resolution 1.14 mm, polarimetric, Atacama Large Millimeter/submillimeter Array (ALMA) observations towards the massive protostar driving the HH 80-81 radio jet. The observations clearly resolve the disk oriented perpendicular to the radio jet, with a radius of ~0.171 arcsec (~291 au at 1.7 kpc distance). The continuum brightness temperature, the intensity profile, and the polarization properties clearly indicate that the disk is optically thick for a radius of R<170 au. The linear polarization of the dust emission is detected almost all along the disk and its properties suggest that dust polarization is produced mainly by self-scattering. However, the polarization pattern presents a clear differentiation between the inner (optically thick) part of the disk and the outer (optically thin) region of the disk, with a sharp transition that occurs at a radius of 0.1 arcsec (~170 au). The polarization characteristics of the inner disk suggest that dust settling has not occurred yet with a maximum dust grain size between 50 and 500 mum. The outer part of the disk has a clear azimuthal pattern but with a significantly higher polarization fraction compared to the inner disk. This pattern is broadly consistent with self-scattering of a radiation field that is beamed radially outward, as expected in the optically thin outer region, although contribution from non-spherical grains aligned with respect to the radiative flux cannot be excluded., Comment: Accepted for publication in the Astrophysical Journal Letters

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., Comment: 28 pages, 16 figures, and 7 tables. Accepted for publication in A&A

We present subarcsecond angular resolution observations carried out with the Submillimeter Array (SMA) at 880 $\mu$m centered at the B0-type protostar GGD27~MM1, the driving source of the parsec scale HH 80-81 jet. We constrain its polarized continuum emission to be $\lesssim0.8\%$ at this wavelength. Its submm spectrum is dominated by sulfur-bearing species tracing a rotating disk--like structure (SO and SO$_2$ isotopologues mainly), but also shows HCN-bearing and CH$_3$OH lines, which trace the disk and the outflow cavity walls excavated by the HH 80-81 jet. The presence of many sulfurated lines could indicate the presence of shocked gas at the disk's centrifugal barrier or that MM1 is a hot core at an evolved stage. The resolved SO$_2$ emission traces very well the disk kinematics and we fit the SMA observations using a thin-disk Keplerian model, which gives the inclination (47$^{\circ}$), the inner ($\lesssim170$ AU) and outer ($\sim950-1300$~AU) radii and the disk's rotation velocity (3.4 km s$^{-1}$ at a putative radius of 1700 AU). We roughly estimate a protostellar dynamical mass of 4-18\msun. MM2 and WMC cores show, comparatively, an almost empty spectra suggesting that they are associated with extended emission detected in previous low-angular resolution observations, and therefore indicating youth (MM2) or the presence of a less massive object (WMC)., Comment: 11 pages, 8 figures. Accepted for publication in the Astrophysical Journal

We targeted the massive star forming region W33A using the Atacama Large Sub/Millimeter Array (ALMA) in band 6 (230 GHz) and 7 (345 GHz) to search for a sub-1000au disc around the central O-type massive young stellar object (MYSO) W33A MM1-Main. Our data achieves a resolution of ~0.2" (~500au) and resolves the central core, MM1, into multiple components and reveals complex and filamentary structures. There is strong molecular line emission covering the entire MM1 region. The kinematic signatures are inconsistent with only Keplerian rotation although we propose that the shift in the emission line centroids within ~1000au of MM1-Main could hint at an underlying compact disc with Keplerian rotation. We cannot however rule out the possibility of an unresolved binary or multiple system. A putative smaller disc could be fed by the large scale spiral `feeding filament' we detect in both gas and dust emission. We also discuss the nature of the now-resolved continuum sources.

Yellow hypergiants are rare and represent a fast evolutionary stage of massive evolved stars. That evolutionary phase is characterised by a very intense mass loss, the understanding of which is still very limited. Here we report ALMA Compact Array observations of a 50$"$-mosaic toward the Fried Egg nebula, around one of the few Galactic yellow hypergiants IRAS 17163-3907. The emission from the $^{12}$CO J=2-1 line, H30$\alpha$ recombination line, and continuum is imaged at a resolution of $\sim$8$"$, revealing the morphology of the molecular environment around the star. The continuum emission is unresolved and peaks at the position of the star. The radio recombination line H30$\alpha$ shows unresolved emission at the star, with an approximately gaussian spectrum centered on a velocity of 21$\pm$3~km/s with a width of 57$\pm$6~km/s. In contrast, the CO 2-1 emission is complex and decomposes into several components beyond the contamination from interstellar gas in the line of sight. The CO spectrum toward the star is a broad plateau, centered at the systemic velocity of +18 km/s and with an expansion velocity of 100$\pm$10 km/s. Assuming isotropic and constant mass-loss, we estimate a mass-loss rate of 8$\pm$1.5 $\times10^{-5}$~M$_\odot$ yr$^{-1}$. At a radius of 25$"$ from the star, we detect CO emission associated with the dust ring previously imaged by {\it Herschel}. The kinematics of this ring, however, is not consistent with an expanding shell, but show a velocity gradient of $v_{sys} \pm$20 km/s. In addition, we find a puzzling bright feature radially connecting the star to the CO ring, at a velocity of +40 km/s relative to the star. This spur feature may trace a unidirectional ejection event from the star. Our ACA observations reveal the complex morphology around IRAS 17163 and illustrate the breakthroughs that ALMA will bring to the field of massive stellar evolution., Comment: Accepted for publication in A&A. 8 pages, 7 figures, plus appendix

Young massive clusters (YMCs) have central stellar mass surface densities exceeding $10^{4} M_{\odot} pc^{-2}$. It is currently unknown whether the stars formed at such high (proto)stellar densities. We compile a sample of gas clouds in the Galaxy which have sufficient gas mass within a radius of a few parsecs to form a YMC, and compare their radial gas mass distributions to the stellar mass distribution of Galactic YMCs. We find that the gas in the progenitor clouds is distributed differently than the stars in YMCs. The mass surface density profiles of the gas clouds are generally shallower than the stellar mass surface density profiles of the YMCs, which are characterised by prominent dense core regions with radii ~ 0.1 pc, followed by a power-law tail. On the scale of YMC core radii, we find that there are no known clouds with significantly more mass in their central regions when compared to Galactic YMCs. Additionally, we find that models in which stars form from very dense initial conditions require surface densities that are generally higher than those seen in the known candidate YMC progenitor clouds. Our results show that the quiescent, less evolved clouds contain less mass in their central regions than in the highly star-forming clouds. This suggests an evolutionary trend in which clouds continue to accumulate mass towards their centres after the onset of star formation. We conclude that a conveyor-belt scenario for YMC formation is consistent with the current sample of Galactic YMCs and their progenitor clouds., Comment: 11 pages, 4 figures, 4 tables. Accepted for publication in MNRAS

We have recently published observations of significant flux density variations at 1.3 cm in HII regions in the star forming regions Sgr B2 Main and North (De Pree et al. 2014). To further study these variations, we have made new 7 mm continuum and recombination line observations of Sgr B2 at the highest possible angular resolution of the Karl G. Jansky Very Large Array (VLA). We have observed Sgr B2 Main and North at 42.9 GHz and at 45.4 GHz in the BnA configuration (Main) and the A configuration (North). We compare these new data to archival VLA 7 mm continuum data of Sgr B2 Main observed in 2003 and Sgr B2 North observed in 2001. We find that one of the 41 known ultracompact and hypercompact HII regions in Sgr B2 (K2-North) has decreased $\sim$27% in flux density from 142$\pm$14 mJy to 103$\pm$10 mJy (2.3$\sigma$) between 2001 and 2012. A second source, F3c-Main has increased $\sim$30% in flux density from 82$\pm$8 mJy to 107 $\pm$11 mJy (1.8$\sigma$) between 2003 and 2012. F3c-Main was previously observed to increase in flux density at 1.3 cm over a longer time period between 1989 and 2012 (De Pree et al. 2014). An observation of decreasing flux density, such as that observed in K2-North, is particularly significant since such a change is not predicted by the classical hypothesis of steady expansion of HII regions during massive star accretion. Our new observations at 7 mm, along with others in the literature, suggest that the formation of massive stars occurs through time-variable and violent accretion., Comment: 20 pages, 5 figures, Accepted for publication in The Astrophysical Journal

The formation of massive stars and their arrival on the zero-age main-sequence occurs hidden behind dense clouds of gas and dust. In the giant Hii region NGC 3603, the radiation of a young cluster of OB stars has dispersed dust and gas in its vicinity. At a projected distance of 2:5 pc from the cluster, a bright mid-infrared (mid-IR) source (IRS 9A) had been identified as a massive young stellar object (MYSO), located on the side of a molecular clump (MM2) of gas facing the cluster. We investigated the physical conditions in MM2, based on APEX sub-mm observations using the SABOCA and SHFI instruments, and archival ATCA 3 mm continuum and CS spectral line data. We resolved MM2 into several compact cores, one of them closely associated with IRS 9A. These are likely infrared dark clouds as they do not show the typical hot-core emission lines and are mostly opaque against the mid-IR background. The compact cores have masses of up to several hundred times the solar mass and gas temperatures of about 50 K, without evidence of internal ionizing sources. We speculate that IRS 9A is younger than the cluster stars, but is in an evolutionary state after that of the compact cores.

The object W33 is a giant molecular cloud that contains star forming regions at various evolutionary stages from quiescent clumps to developed H II regions. Since its star forming regions are located at the same distance and the primary material of the birth clouds is probably similar, we conducted a comparative chemical study to trace the chemical footprint of the different phases of evolution. We observed six clumps in W33 with the Atacama Pathfinder Experiment (APEX) telescope at 280 GHz and the Submillimeter Array (SMA) at 230 GHz. We detected 27 transitions of 10 different molecules in the APEX data and 52 transitions of 16 different molecules in the SMA data. The chemistry on scales larger than $\sim$0.2 pc, which are traced by the APEX data, becomes more complex and diverse the more evolved the star forming region is. On smaller scales traced by the SMA data, the chemical complexity and diversity increase up to the hot core stage. In the H II region phase, the SMA spectra resemble the spectra of the protostellar phase. Either these more complex molecules are destroyed or their emission is not compact enough to be detected with the SMA. Synthetic spectra modelling of the H$_{2}$CO transitions, as detected with the APEX telescope, shows that both a warm and a cold component are needed to obtain a good fit to the emission for all sources except for W33 Main1. The temperatures and column densities of the two components increase during the evolution of the star forming regions. The integrated intensity ratios N$_{2}$H$^{+}$(3$-$2)/CS(6$-$5) and N$_{2}$H$^{+}$(3$-$2)/H$_{2}$CO(4$_{2,2}$$-$3$_{2,1}$) show clear trends as a function of evolutionary stage, luminosity, luminosity-to-mass ratio, and H$_{2}$ peak column density of the clumps and might be usable as chemical clocks., Comment: 66 pages, 28 figures, 8 tables, accepted for publication at A&A

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants, Comment: Replaced figures 2 and 3: corrected position SRT. Corrected minor typo in 5.1

Accretion flows onto massive stars must transfer mass so quickly that they are themselves gravitationally unstable, forming dense clumps and filaments. These density perturbations interact with young massive stars, emitting ionizing radiation, alternately exposing and confining their HII regions. As a result, the HII regions are predicted to flicker in flux density over periods of decades to centuries rather than increasing monotonically in size as predicted by simple Spitzer solutions. We have recently observed the Sgr B2 region at 1.3 cm with the VLA in its three hybrid configurations (DnC, CnB and BnA) at a resolution of 0.25''. These observations were made to compare in detail with matched continuum observations from 1989. At 0.25'' resolution, Sgr B2 contains 41 UC HII regions, 6 of which are hypercompact. The new observations of Sgr B2 allow comparison of relative peak flux densites for the HII regions in Sgr B2 over a 23 year time baseline (1989-2012) in one of the most source-rich massive star forming regions in the Milky Way. The new 1.3 cm continuum images indicate that four of the 41 UC HII regions exhibit significant changes in their peak flux density, with one source (K3) dropping in peak flux density, and the other 3 sources (F10.303, F1 and F3) increasing in peak flux density. The results are consistent with statistical predictions from simulations of high mass star formation, suggesting that they offer a solution to the lifetime problem for ultracompact HII regions., Comment: 12 pages, 3 figures, Accepted for publication in the Astrophysical Journal Letters

The Multi-Scale Continuum and Line Exploration of W49 (MUSCLE W49) is a comprehensive gas and dust survey of the giant molecular cloud (GMC) of W49A, the most luminous star-formation region in the Milky Way. It covers the entire GMC at different scales and angular resolutions. In this paper we present: 1) an all-configuration SMA mosaic in the 230-GHz band covering the central 3 arcmin (10 pc, known as W49N), with most of the embedded massive stars; and 2) PMO 14m telescope observations in the 90-GHz band, covering the entire GMC with maps up to 35 arcmin in size, or 113 pc. We also make use of archival data from the VLA, JCMT-SCUBA, IRAM 30m, and the CSO BOLOCAM GPS. Our main findings are: 1) The W49 GMC is one of the most massive in the Galaxy, with a total mass ~1.1x10^6 Msun within a radius of 60 pc. Within a radius of 6 pc, the total gas mass is ~2x10^5 Msun. At these scales only 1% of the material is photoionized. The mass reservoir is sufficient to form several young massive clusters (YMCs) as massive as a globular cluster. 2) The mass of the GMC is distributed in a hierarchical network of filaments. At scales <10 pc, a triple, centrally condensed structure peaks toward the ring of HC HII regions in W49N. This structure extends to scales from ~10 to 100 pc. The W49A starburst most likely formed from global gravitational contraction with localized collapse in a "hub-filament" geometry. 3) Currently, feedback from the central YMCs (with a present mass Mcl > 5x10^4 Msun) is still not enough to entirely disrupt the GMC, but further stellar mass growth could be enough to allow radiation pressure to clear the cloud and halt star formation. 4) The resulting stellar content will probably remain as a gravitationally bound massive star cluster, or a small system of bound clusters. (ABRIDGED), Comment: Author version of paper published in ApJ. For high-quality figures please see the Journal version

Whether high mass stars continue to accrete material beyond the formation of an HII region is still an open question. Ionized infall and outflow have been seen in some sources, but their ties to the surrounding molecular gas are not well constrained. We aim to quantify the ionized and molecular gas dynamics in a high mass star forming region (K3-50A) and their interaction. We present CARMA observations of the 3mm continuum, H41\alpha, and HCO+ emission, and VLA continuum observations at 23 GHz and 14.7 GHz to quantify the gas and its dynamics in K3-50A. We find large scale dynamics consistent with previous observations. On small scales, we find evidence for interaction between the ionized and molecular gas which suggests the ionized outflow is entraining the molecular one. This is the first time such an outflow entrained by photo ionized gas has been observed. Accretion may be ongoing in K3-50A because an ionized bipolar outflow is still being powered, which is in turn entraining part of the surrounding molecular gas. This outflow scenario is similar to that predicted by ionization feedback models., Comment: 7 pages, 4 Figures. Accepted by Astronomy and Astrophysics

Context. This study is part of a large project to study the physics of accretion and molecular outflows towards a selected sample of high-mass star-forming regions that show evidence of infall and rotation from previous studies. Aims. We wish to make a thorough study at high-angular resolution of the structure and kinematics of the HMCs and corresponding molecular outflows in the high-mass star-forming region G24.78+0.08. Methods. We carried out SMA and IRAM PdBI observations at 1.3 and 1.4 mm, respectively, of dust and of typical high-density and molecular outflow tracers with resolutions of <1". Complementary IRAM 30-m 12CO and 13CO observations were carried out to recover the short spacing information of the molecular outflows. Results. The millimeter continuum emission towards cores G24 A1 and A2 has been resolved into 3 and 2 cores, respectively, and named A1, A1b, A1c, A2, and A2b. All these cores are aligned in a southeast-northwest direction coincident with that of the molecular outflows detected in the region, which suggests a preferential direction for star formation in this region. The masses of the cores range from 7 to 22 Msun, and the rotational temperatures from 128 to 180 K. The high-density tracers have revealed the existence of 2 velocity components towards A1, one of them peaks close to the position of the millimeter continuum peak and of the HC HII region, and is associated with the velocity gradient seen in CH3CN towards this core, while the other one peaks southwest of core A1 and is not associated with any millimeter continuum emission peak. The position-velocity plots along outflow A and the 13CO averaged blueshifted and redshifted emission indicate that this outflow is driven by core A2., Comment: 18 pages, 3 tables, 15 figures. Accepted for publication in A&A