Context. Complex organic molecules (COMs) are found to be abundant in various astrophysical environments, particularly toward star-forming regions, where they are observed both toward protostellar envelopes as well as shocked regions. The emission spectrum, especially that of heavier COMs, might consist of up to hundreds of lines, where line blending hinders the analysis. However, identifying the molecular composition of the gas that leads to the observed millimeter spectra is the first step toward a quantitative analysis. Aims. We have developed a new method based on supervised machine learning to recognize spectroscopic features of the rotational spectrum of molecules in the 3 mm atmospheric transmission band for a list of species including COMs, with the aim of obtaining a detection probability. Methods. We used local thermodynamic equilibrium (LTE) modeling to build a large set of synthetic spectra of 20 molecular species, including COMs with a range of physical conditions typical for star-forming regions. We successfully designed and trained a convolutional neural network (CNN) that provides detection probabilities of individual species in the spectra. Results. We demonstrate that the CNN model we developed has a robust performance to detect spectroscopic signatures from these species in synthetic spectra. We evaluated its ability to detect molecules according to the noise level, frequency coverage, and line-richness, as well as to test its performance for an incomplete frequency coverage with high detection probabilities for the tested parameter space, with no false predictions. Finally, we applied the CNN model to obtain predictions on observational data from the literature toward line-rich hot core-like sources, where the detection probabilities remain reasonable, with no false detections. Conclusions. We demonstrate the use of CNNs in facilitating the analysis of complex millimeter spectra both on synthetic spectra, along with the first tests performed on observational data. Further analyses on its explainability, as well as calibration using a larger observational dataset, will help improve the performance of our method for future applications.
Context. High-mass prestellar cores are extremely rare. Until recently, the search for such objects has been hampered by small sample sizes, leading to large ambiguities in their lifetimes and hence the conditions in the cores in which high-mass stars (greater than or similar to 8 M-circle dot) form. Aims. Here we leverage the large sample (similar to 580 cores) detected in the ALMA-IMF survey to identify both protostellar and prestellar cores to estimate their relative lifetimes. Methods. We used CO and SiO outflows to identify protostellar cores. We present a new automated method based on aperture line emission and background subtraction to systematically detect outflows associated with each of the 141 most massive cores. Massive cores that are not driving an outflow in either tracer are identified as prestellar. After careful scrutiny of the sample, we derived statistical lifetime estimates for the prestellar phase. Results. Our automated method allows the efficient detection of CO and SiO outflows and has a performance efficiency similar to that of more cumbersome classical techniques. We identified 30 likely prestellar cores with M greater than or similar to 8 M-circle dot, of which 12 have core masses M greater than or similar to 16 M-circle dot. The latter group contains the best candidates for high-mass star precursors. Moreover, most of these 12 high-mass prestellar cores are located inside the crowded central regions of the protoclusters, where most high-mass stars are expected to form. Using the relative ratios of prestellar to protostellar cores, and assuming a high-mass protostellar lifetime of 300 kyr, we derive a prestellar core lifetime of 120 kyr to 240 kyr for cores with masses 8 M-circle dot < M < 16 M-circle dot. For 30 M-circle dot < M < 55 M-circle dot, the lifetimes range from 50 kyr to 100 kyr. The spread in timescales reflects different assumptions for scenarios for the mass reservoir evolution. These timescales are remarkably long compared to the 4 kyr to 15 kyr free-fall time of the cores. Hence, we suggest that high-mass cores live similar to 10 to 30 free-fall times, with a tentative trend of a slight decrease with core mass. Such high ratios suggest that the collapse of massive cores is slowed down by non-thermal support of turbulent, magnetic or rotational origin at or below the observed scale.
The 'Diamond Ring' within Cygnus X, south-west of the DR21 ridge, stands out as a prominent, nearly circular structure in infrared (IR) and far-infrared (FIR) emission, spanning approximately 6 pc in diameter. It is enclosed by clumpy molecular clouds seen in CO lines and contains an H II region, visible in cm emission. It resembles a classical H II region associated with an expanding bubble seen routinely in the 158 mu m line of ionised carbon ([C II]). However, our recent observations utilising the Stratospheric Observatory for Far-Infrared Astronomy (SOFIA) under the FEEDBACK program for the spectrally resolved [C II] line have revealed a slightly tilted ring with a mass of similar to 10(3) M-circle dot, advancing at a velocity of similar to 1.3 km s(-1). The bulk emission of the gas has a line-of-sight (LOS) velocity around -2 km s(-1). The [C II] data revealed that the 'Diamond' of the Diamond Ring is an unrelated, dense gas clump at a LOS velocity of similar to 7 km s(-1). The driving source, which is also responsible for powering the associated H II region, is a B0.5e star, classified by our IR spectroscopy. This observation marks the first instance where we observe only a slowly expanding ring of [C II] emission and not an expanding 3D shell. We suggest that the H II region (along with its associated [C II] bubble), initially formed by a massive star, expanded outwards from a flat slab of molecular gas nearly in the plane of the sky. Presently, the [C II] ring is confined by the swept-up gas of the slab, while the parts of the shell moving in the directions perpendicular to the shell along the LOS have dissipated, resulting in a notable decrease in the expansion of the remaining ring. This scenario is supported by dedicated simulations that trace the evolution of the [C II] bubble. Our observations support the scenario of H II region evolution in 'flat' molecular clouds, reported earlier in the literature. In this geometry, we propose that the Diamond Ring represents the terminal phase of an expanding [C II] bubble driven by stellar winds and thermal pressure.
Recent high-angular-resolution observations indicate the need for core growth to form high-mass stars. To understand the gas dynamics at the core scale in the very early evolutionary stages before being severely affected by feedback, we have conducted Atacama Large Millimeter/submillimeter Array (ALMA) observations toward a 70 μ m dark massive clump, G337.541-00.082 as part of the Global and Local infall in the ASHES sample (GLASHES) program. Using dense gas tracers such as N 2 H + ( J = 1–0) and HNC ( J = 3–2), we find signs of infall from the position–velocity diagram and more directly from the blue asymmetry profile in addition to the clump-scale velocity gradient. We estimate infall velocities from intermediate and low-mass cores to be 0.28–1.45 km s −1 , and infall rates to be on the order of 10 −4 –10 −3 M ⊙ yr −1 , both are higher than those measured in low-mass star-forming regions by more than a factor of 5 and an order of magnitude, respectively. We find a strong correlation between the infall velocity with the nonthermal velocity dispersion, suggesting that infall may contribute significantly to the observed line width. Consistent with clump-fed scenarios, we show that the mass infall rate is larger for larger core masses and shorter distances to the clump center. Such high infall rates in cores embedded in IRDCs can be considered as strong signs of core growth, allowing high-mass star formation from intermediate-mass cores that would not initially form high-mass stars at their current mass.
Context. The ALMA-IMF Large Program imaged 15 massive protoclusters down to a resolution of similar to 2 kau scales, identifying about 10(3) 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. Aims. 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 one-sixth of the sample, we measured the bolometric luminosities, implementing deblending corrections when necessary. Methods. We used spectral energy distribution (SED) analysis obtained with the point process mapping (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. Results. For the first time, we provide data-informed estimates of dust temperatures for 882 cores identified with ALMA-IMF: 17-31 K and 28-79 K (5th and 95th percentiles, up to 127 K) for the 616 prestellar and 266 protostellar cores, respectively. We also measured protostellar luminosities spanning 20-80 000 L-circle dot. Conclusions. Dust temperatures previously estimated from SED-based analyses at a comparatively lower resolution validate our method. For hot cores, on the other hand, 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 intermediate- to 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.
Context . A crucial aspect in addressing the challenge of measuring the core mass function (CMF), that is pivotal for comprehending the origin of the initial mass function (IMF), lies in constraining the temperatures of the cores. Aims . We aim to measure the luminosity, mass, column density and dust temperature of star-forming regions imaged by the ALMA-IMF large program. These fields were chosen to encompass early evolutionary stages of massive protoclusters. High angular resolution mapping is required to capture the properties of protostellar and pre-stellar cores within these regions, and to effectively separate them from larger features, such as dusty filaments. Methods . 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. Results . 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.
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.
The ALMA-IMF Large Program provides multi-tracer observations of 15 Galactic massive protoclusters at a matched sensitivity and spatial resolution. We focus on the dense gas kinematics of the G353.41 protocluster traced by N2H+ (1-0), with a spatial resolution of similar to 0.02 pc. G353.41, at a distance of similar to 2kpc, is embedded in a larger-scale (similar to 8 pc) filament and has a mass of similar to 2.5 x 10(3) M-circle dot within 1.3 x 1.3 pc(2). We extracted the N2H+ (1-0) isolated line component and decomposed it by fitting up to three Gaussian velocity components. This allows us to identify velocity structures that are either muddled or impossible to identify in the traditional position-velocity diagram. We identify multiple velocity gradients on large (similar to 1 pc) and small scales (similar to 0.2pc). We find good agreement between the N2H+ velocities and the previously reported DCN core velocities, suggesting that cores are kinematically coupled with the dense gas in which they form. We have measured nine converging "V-shaped" velocity gradients (VGs) (similar to 20 km s(-1) pc(-1)) that are well resolved (sizes similar to 0.1 pc), mostly located in filaments, which are sometimes associated with cores near their point of convergence. We interpret these V-shapes as inflowing gas feeding the regions near cores (the immediate sites of star formation). We estimated the timescales associated with V-shapes as VG(-1), and we interpret them as inflow timescales. The average inflow timescale is similar to 67 kyr, or about twice the free-fall time of cores in the same area (similar to 33 kyr) but substantially shorter than protostar lifetime estimates (similar to 0.5 Myr). We derived mass accretion rates in the range of (0.35-8.77) x 10(-4) M-circle dot yr(-1). This feeding might lead to further filament collapse and the 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. Thus, 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's lifetime.
Context. During the process of star formation, the dense gas undergoes significant chemical evolution leading to the emergence of a rich variety of molecules associated with hot cores and hot corinos. However, the physical conditions and the chemical processes involved in this evolution are poorly constrained. In particular, the early phases, corresponding to a stage prior to the emergence of any strong ionising emission from the protostar, are still poorly studied. Aims. In this work, we provide a full inventory of the emission from complex organic molecules (COMs) to investigate the physical structure and chemical composition of six high-mass protostellar envelopes. We aim to investigate the conditions for the emergence of COMs in hot cores. Methods. We performed an unbiased spectral survey towards six infrared-quiet massive clumps between 159 GHz and 374 GHz with the APEX 12 m telescope, covering the entire atmospheric windows at 2 mm, 1.2 mm, and 0.8 mm. To identify the spectral lines, we used rotational diagrams and radiative transfer modelling assuming local thermodynamic equilibrium. Results. We detect up to 11 COMs plus three isotopologues, of which at least five COMs (CH3OH, CH3CN, CH3OCHO, CH3OCH3, and CH3CHO) are detected towards all sources. Towards all the objects, most of the COM emission is found to be cold, with respect to the typical temperatures at which COMs are found, with a temperature of 30 K and extended with a size of ~0.3 pc. Although the bulk of the gas for our sample of young massive clumps has a cold temperature, we also detect emission from COMs originating from the immediate vicinity of the protostar. This warm component of the envelope is best traced by methanol and methyl cyanide, in particular methyl cyanide traces a compact (~1″) and the hottest (T ~200 K) component of the envelope. Only three out of the six sources exhibit a robustly detected hot gas component (T > 100 K) traced by several COMs. We find a gradual emergence of the warm component in terms of size and temperature, together with an increasing molecular complexity, allowing us to establish an evolutionary sequence for our sample based on COMs. While they can already be well characterised by an emerging molecular richness, gas temperatures of COMs in the hot gas and molecular abundances suggest that COMs may become abundant in the gas phase at temperatures below the thermal desorption temperature. Conclusions. Our findings confirm that the sources of our sample of infrared-quiet massive clumps are in an early evolutionary stage during which the bulk of the gas is cold. The presence of COMs is found to be characteristic of these early evolutionary stages accompanying high-mass star and cluster formation. While the extent of the compact heated regions resembles that of hot cores, the molecular abundances, except for complex cyanides, resemble those of hot corinos and are lower than the peak COM abundances of hot cores. We suggest that the emergence of hot cores is preceded by a phase in which mostly O-bearing COMs appear first with similar abundances to hot corinos albeit with larger source sizes.
We use the H41 α recombination line to create templates of the millimeter free–free emission in the ALMA-IMF continuum maps, which allows us to separate it from dust emission. This method complements spectral-index information and extrapolation from centimeter-wavelength maps. We use the derived maps to estimate the properties of up to 34 H ii regions across the ALMA-IMF protoclusters. The hydrogen ionizing photon rate Q 0 and spectral types follow the evolutionary trend proposed by Motte et al. The youngest protoclusters lack detectable ionized gas, followed by protoclusters with increasing numbers of OB stars. The total Q 0 increases from ∼10 45 s −1 to >10 49 s −1 . We used the adjacent He41 α line to measure the relative number abundances of helium, finding values consistent with the Galactic interstellar medium, although a few outliers are discussed. A search for sites of maser amplification of the H41 α line returned negative results. We looked for possible correlations between the electron densities, emission measures, and Q 0 with H ii region size D . The latter is the best correlated, with Q 0 ∝ D 2.49 ± 0.18 . This favors interpretations in which smaller ultracompact H ii regions are not necessarily the less dynamically evolved versions of larger ones but rather are ionized by less massive stars. Moderate correlations were found between the dynamical width Δ V dyn with D and Q 0 . Δ V dyn increases from about 1 to 2 times the ionized-gas sound speed. Finally, an outlier H ii region south of W43-MM2 is discussed. We suggest that this source could harbor an embedded stellar or disk wind.
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 use the H41α recombination line to create templates of the millimeter free-free emission in the ALMA-IMF continuum maps, which allows to separate it from dust emission. This method complements spectral-index information and extrapolation from centimeter wavelength maps. We use the derived maps to estimate the properties of up to 34 HII regions across the ALMA-IMF protoclusters. The hydrogen ionizing-photon rate Q_0 and spectral types follow the evolutionary trend proposed by Motte et al. The youngest protoclusters lack detectable ionized gas, followed by protoclusters with increasing numbers of OB stars. The total Q_0 increases from ∼ 10^45 s^-1 to > 10^49 s^-1. We used the adjacent He41α line to measure the relative number abundances of helium, finding values consistent with the Galactic interstellar medium, although a few outliers are discussed. A search for sites of maser amplification of the H41α line returned negative results. We looked for possible correlations between the electron densities (n_e), emission measures (EM), and Q_0 with HII region size D. The latter are the better correlated, with Q_0 ∝ D^2.49±0.18. This favors interpretations where smaller ultracompact HII regions are not necessarily the less dynamically evolved versions of larger ones, but rather are ionized by less massive stars. Moderate correlations were found between dynamical width Δ V_dyn with D and Q_0. Δ V_dyn increases from about one to two times the ionized-gas sound speed. Finally, an outlier HII region south of W43-MM2 is discussed. We suggest that this source could harbor an embedded stellar or disk wind.
We present a catalogue of dense cores and filaments in a 3.8(degrees)x2.4(degrees) field around the TMC1 region of the Taurus molecular cloud. The catalogue was created using photometric data from the Herschel SPIRE and PACS instruments in the 70, 160, 250, 350, and 500 mu m continuum bands. Extended structure in the region was reconstructed from a Herschel column density map. Power spectra and probability density functions (PDFs) of this structure are presented. The PDF splits into lognormal and power-law forms, with the high-density power-law component associated primarily with the central part of TMC1. The total mass in the mapped region is 2000 M-circle dot, of which 34 per cent is above an extinction of A(V)similar to 3 mag - a level that appears as a break in the PDF and as the minimum column density at which dense cores are found. A total of 35 dense filaments were extracted from the column density map. These have a typical full width at half-maximum (FWHM) width of 0.07 pc, but the TMC1 filament itself has a mean FWHM of similar to 0.13 pc. The thermally supercritical filaments in the region are aligned orthogonal to the prevailing magnetic field direction. Derived properties for the supercritical TMC1 filament support the scenario of it being relatively young. A catalogue of 44 robust and candidate prestellar cores is created and is assessed to be complete down to 0.1 M-circle dot. The combined prestellar core mass function for the TMC1 and L1495 regions is well fit by a single lognormal distribution and is comparable to the standard initial mass function.
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).
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.
Context. One of the central questions in astrophysics is the origin of the initial mass function (IMF). It is intrinsically linked to the processes from which it originates, and hence its connection with the core mass function (CMF) must be elucidated. Aims. We aim 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 that are part of the ALMA-IMF program. Methods. We used observations from the ALMA-IMF large programme: ~2′ × 2′ maps of emission from the continuum and selected lines at 1.3 mm and 3 mm observed by the ALMA 12m only antennas. Our angular resolution was typically 1″, that is, ~2400 au at a distance of 2.4 kpc. The lines we analysed are CO (2–1), SiO (5–4), N2H+ (1–0), H41α as well as He41α blended with C41α. We built a census of dense cores in the region, and we measured the associated CMF based on a core-dependent temperature value. Results. We confirmed the ‘evolved’ status of W33-Main by identifiying three H II regions within the field, and to a lesser extent based on the number and extension of N2H+ filaments. We produced a filtered core catalogue 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.2 M⊙. We fitted the resulting high-mass end of the CMF with a single power law of the form N(log(M)) ∝ Mα, obtaining α = −1.44−0.22+0.16, which is slightly steeper but consistent with the Salpeter index. We categorised our cores as prestellar and protostellar, mostly based on outflow 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 H II regions and their surroundings than in the rest of the field. We also found that the cores’ masses were rather low (maximum mass of ~13 M⊙). Conclusions. We find that star formation in W33-Main could be compatible with a ‘clump-fed’ scenario of star formation in an evolved cloud characterised by stellar feedback in the form of H II regions, and under the influence of massive stars outside the field. Our results differ from those found in less evolved young star-forming regions in the ALMA-IMF program. Further investigations are needed to elucidate the evolution of late CMFs towards the IMF over statistically significant samples.
We present significant improvements to our previous work on noise reduction in Herschel observation maps by defining sparse filtering tools capable of handling, in a unified formalism, a significantly improved noise reduction as well as a deconvolution in order to reduce effects introduced by the limited instrumental response (beam). We implement greater flexibility by allowing a wider choice of parsimonious priors in the noise-reduction process. More precisely, we introduce a sparse filtering and deconvolution approach approach of type l^2-l^p, with p > 0 variable and apply it to a larger set of molecular clouds using Herschel 250 μm data in order to demonstrate their wide range of application. In the Herschel data, we are able to use this approach to highlight extremely fine filamentary structures and obtain singularity spectra that tend to show a significantly less log-normal behavior and a filamentary nature in the less dense regions. We also use high-resolution adaptive magneto-hydrodynamic simulation data to assess the quality of deconvolution in such a simulated beaming framework.
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.
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 12 m-array configurations. The data include 12 spectral windows, with eight at 1.3 mm and four at 3 mm. 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 (3–2) 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.
Context. 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. Aims. We aim to compare populations of protostellar and prestellar cores in three regions imaged in the ALMA-IMF Large Program. Methods. 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. Results. We found a similar fraction of cores which are protostellar in the two regions, about 35%. This fraction strongly varies in mass, from f(pro).similar or equal to 15-20% at low mass, between 0.8 and 3 M-circle dot up to f(pro) similar or equal to 80% above 16 M-circle dot. 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.12) (-0.19), is consistent with the Salpeter slope, and thus the top-heavy form measured for the global CMF, alpha = -0.96 +/- 0.09, is due to the protostellar core population. Conclusions. 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.