Showing 1–50 of 106 results for author: Gálvan-Madrid, R

ALMA-IMF XV: The core mass function in the high-mass star-formation regime

Authors: F. Louvet, P. Sanhueza, A. Stutz, A. Men'shchikov, F. Motte, R. Galván-Madrid, S. Bontemps, Y. Pouteau, A. Ginsburg, T. Csengeri, J. Di Francesco, P. Dell'Ova, M. González, P. Didelon, J. Braine, N. Cunningham, B. Thomasson, P. Lesaffre, P. Hennebelle, M. Bonfand, A. Gusdorf, R. H. Álverez-Gutiérrez, T. Nony, G. Busquet, F. Olguin , et al. (16 additional authors not shown)

Abstract: 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 I… ▽ More 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. △ Less

Submitted 26 July, 2024; originally announced July 2024.

Comments: Accepted by A&A

ALMA-IMF XII: Point-process mapping of 15 massive protoclusters

Authors: P. Dell'Ova, F. Motte, A. Gusdorf, Y. Pouteau, A. Men'shchikov, D. Diaz-Gonzalez, R. Galván-Madrid, P. Lesaffre, P. Didelon, A. M. Stutz, A. P. M. Towner, K. Marsh, A. Whitworth, M. Armante, M. Bonfand, T. Nony, M. Valeille-Manet, S. Bontemps, T. Csengeri, N. Cunningham, A. Ginsburg, F. Louvet, R. H. Alvarez-Gutierrez, N. Brouillet, J. Salinas , et al. (7 additional authors not shown)

Abstract: 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 captur… ▽ More 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. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 35 pages, 19 figures, 5 tables Accepted by A&A

ALMA-IMF XIV: Free-Free Templates Derived from H$41α$ and Ionized Gas Content in Fifteen Massive Protoclusters

Authors: Roberto Galván-Madrid, Daniel J. Díaz-González, Frédérique Motte, Adam Ginsburg, Nichol Cunningham, Karl M. Menten, Mélanie Armante, Mélisse Bonfand, Jonathan Braine, Timea Csengeri, Pierre Dell'Ova, Fabien Louvet, Thomas Nony, Rudy Rivera-Soto, Patricio Sanhueza, Amelia M. Stutz, Friedrich Wyrowski, Rodrigo H. Álvarez-Gutiérrez, Tapas Baug, Sylvain Bontemps, Leonardo Bronfman, Manuel Fernández-López, Antoine Gusdorf, Atanu Koley, Hong-Li Liu , et al. (3 additional authors not shown)

Abstract: We use the H$41α$ 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 i… ▽ More We use the H$41α$ 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 $\sim 10^{45}$ s$^{-1}$ to $> 10^{49}$ s$^{-1}$. We used the adjacent He$41α$ 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 H$41α$ 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 \propto D^{2.49\pm0.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_\mathrm{dyn}$ with $D$ and $Q_0$. $ΔV_\mathrm{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. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: Accepted to AAS Journals (ApJS). Data available at https://dataverse.harvard.edu/dataverse/h41a-freefree

Next Generation Very Large Array Memo #122: Characterization of the synthesized beam with and without MID antennas in Mexico

Authors: Alfonso Trejo-Cruz, Roberto Galván-Madrid, Carlos Carrasco-González, Eric F. Jiménez-Andrade, Stan Kurtz, Jesús M. Jáquez-Domínguez, Alice Pasetto, Luis A. Zapata

Abstract: Synthesized beam (PSF) synthetic observations with and without the antennas in Mexico are analyzed. For a simple continuum observing setup, we generated visibility files and their associated PSF images for a grid of parameters (robust weighting, tapering, and declination). The tests were done for both the MID and MID+Spiral+Core configurations and their cropped versions without antennas in Mexico.… ▽ More Synthesized beam (PSF) synthetic observations with and without the antennas in Mexico are analyzed. For a simple continuum observing setup, we generated visibility files and their associated PSF images for a grid of parameters (robust weighting, tapering, and declination). The tests were done for both the MID and MID+Spiral+Core configurations and their cropped versions without antennas in Mexico. We show that the performance of the Array, in terms of the beam properties, is in general significantly better when both the MID array antennas are present in Northern Mexico and observations target southern sources. At a declination of -40 deg, there are increments in the ellipticity of at least ~1.3X and 1.2X for a tapering of 3.0 and 4.0 mas, if the antennas in Mexico are not included. For the parameter space tested, the changes in ellipticity of the MID and MID+Spiral+Core configurations differ by ~10%. Larger tapering values help to reduce the ellipticity for cropped configurations at all declinations, but it will impose more constraints in terms of angular resolution. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 10 pages, 8 figures

ALMA-IMF XIII: N$_2$H$^+$ kinematic analysis on the intermediate protocluster G353.41

Authors: R. H. Álvarez-Gutiérrez, A. M. Stutz, N. Sandoval-Garrido, F. Louvet, F. Motte, R. Galván-Madrid, N. Cunningham, P. Sanhueza, M. Bonfand, S. Bontemps, A. Gusdorf, A. Ginsburg, T. Csengeri, S. D. Reyes, J. Salinas, T. Baug, L. Bronfman, G. Busquet, D. J. Díaz-González, M. Fernandez-Lopez, A. Guzmán, A. Koley, H. -L. Liu, F. A. Olguin, M. Valeille-Manet , et al. (1 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 13 June, 2024; v1 submitted 10 April, 2024; originally announced April 2024.

Comments: Accepted in A&A, 22 pages, 21 figures, 1 interactive figure, 4 tables

INvestigations of massive Filaments ANd sTar formation (INFANT). I. Core Identification and Core Mass Function

Authors: Yu Cheng, Xing Lu, Patricio Sanhueza, Hauyu Baobab Liu, Qizhou Zhang, Roberto Galván-Madrid, Ke Wang, Fumitaka Nakamura, Tie Liu, Siyi Feng, Shanghuo Li, Sihan Jiao, Kei E. I. Tanaka, Xunchuan Liu, Pak Shing Li, Qiuyi Luo, Qilao Gu, Yuxin Lin, András E. Guzmán

Abstract: Filamentary structures are ubiquitously found in high-mass star-forming clouds. To investigate the relationship between filaments and star formation, we carry out the INFANT (INvestigations of massive Filaments ANd sTar formation) survey, a multi-scale, multi-wavelength survey of massive filamentary clouds with ALMA band 3/band 6 and VLA K band. In this first paper, we present the ALMA band 6 cont… ▽ More Filamentary structures are ubiquitously found in high-mass star-forming clouds. To investigate the relationship between filaments and star formation, we carry out the INFANT (INvestigations of massive Filaments ANd sTar formation) survey, a multi-scale, multi-wavelength survey of massive filamentary clouds with ALMA band 3/band 6 and VLA K band. In this first paper, we present the ALMA band 6 continuum observations toward a sample of 8 high-mass star forming filaments. We covered each target with approximately rectangular mosaic field of view with two 12-m array configurations, achieving an angular resolution of $\sim$0.6" (2700 AU at 4.5 kpc) and a continuum rms of $\sim$0.1 mJy/beam ($\sim$0.06 Msun in gas mass assuming 15 K). We identify cores using the getsf and astrodendro and find the former is more robust in terms of both identification and measuring flux densities. We identify in total 183 dense cores (15--36 cores in each cloud) and classify their star formation states via outflow and warm gas tracers. The protostellar cores are statistically more massive than the prestellar cores, possibly indicating further accretion onto cores after formation of protostars. For the high-mass end ($M_\text{core}$ $>$ 1.5 Msun) of the core mass function (CMF) we derive a power-law index of $-$1.15 $\pm$ 0.12 for the whole sample, and $-$1.70 $\pm$ 0.25 for the prestellar population. We also find a steepening trend in CMF with cloud evolution ($-$0.89 $\pm$ 0.15 for the young group v.s. $-$1.44 $\pm$ 0.25 for the evolved group) and discuss its implication for cluster formation. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: 25 pages, 8 figures, accepted for ApJ

Core to ultracompact HII region evolution in the W49A massive protocluster

Authors: T. Nony, R. Galván-Madrid, N. Brouillet, G. Suárez, F. Louvet, C. G. De Pree, M. Juárez-Gama, A. Ginsburg, K. Immer, Y. Lin, H. B. Liu, C. G. Román-Zúñiga, Q. Zhang

Abstract: 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 contin… ▽ More 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$α$ 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$α$. 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. △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: 15 pages, 7 figures. Accepted for publication in A&A

Journal ref: A&A 687, A84 (2024)

ALMA-IMF XI: The sample of hot core candidates A rich population of young high-mass proto-stars unveiled by the emission of methyl formate

Authors: M. Bonfand, T. Csengeri, S. Bontemps, N. Brouillet, F. Motte, F. Louvet, A. Ginsburg, N. Cunningham, R. Galván-Madrid, F. Herpin, F. Wyrowski, M. Valeille-Manet, A. M. Stutz, J. Di Francesco, A. Gusdorf, M. Fernández-López, B. Lefloch, H-L. Liu, P. Sanhueza, R. H. Álvarez-Gutiérrez, F. Olguin, T. Nony, A. Lopez-Sepulcre, P. Dell'Ova, Y. Pouteau , et al. (6 additional authors not shown)

Abstract: 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 pea… ▽ More 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. △ Less

Submitted 22 February, 2024; originally announced February 2024.

Journal ref: A&A 687, A163 (2024)

ALMA-IMF X -- The core population in the evolved W33-Main (G012.80) protocluster

Authors: M. Armante, A. Gusdorf, F. Louvet, F. Motte, Y. Pouteau, P. Lesaffre, R. Galván-Madrid, P. Dell'Ova, M. Bonfand, T. Nony, N. Brouillet, N. Cunningham, A. Ginsburg, A. Men'shchikov, S. Bontemps, D. Díaz González, T. Csengeri, M. Fernández-López, M. González, F. Herpin, H. -L. Liu, P. Sanhueza, A. M. Stutz, M. Valeille-Manet

Abstract: 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 reso… ▽ More 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). △ Less

Submitted 17 January, 2024; originally announced January 2024.

Journal ref: A&A 686, A122 (2024)

Shedding light on the ejection history of molecular outflows: Multiple velocity modes and precession

Authors: Veronica Lora, Thomas Nony, Alejandro Esquivel, Roberto Galván-Madrid

Abstract: Variable accretion has been well studied in evolved stages of low-mass stars formation. However, the accretion history in the initial phases of star formation is still a seldom studied topic. The outflows and jets emerging from protostellar objects could shed some light onto their accretion history. We consider the recently studied case of W43-MM1, a protocluster containing 46 outflows driven by 2… ▽ More Variable accretion has been well studied in evolved stages of low-mass stars formation. However, the accretion history in the initial phases of star formation is still a seldom studied topic. The outflows and jets emerging from protostellar objects could shed some light onto their accretion history. We consider the recently studied case of W43-MM1, a protocluster containing 46 outflows driven by 27 protostellar cores. The outflow kinematics of the individual cores and associated knots in W43-MM1 indicate episodic protostellar ejection. We took the observed parameters of an individual core system (core #8) and performed 3D hydrodynamic simulations of such system including episodic changes in the velocity of the outflow. The simulations consist of a collimated jet emerging from a core taking into account one- and two-velocity modes in the variation of the ejection velocity of the jet. In addition, we investigated the effect of including the precession of the jet to the one- and two-velocity mode models. From the simulations, we constructed position-velocity diagrams and compared them with the observations. We find that including a second mode in the ejection velocity, as well as the precession, are required to explain the positions of the outflow knots and other position-velocity features observed in the outflow of core #8 in W43-MM1. △ Less

Submitted 20 December, 2023; originally announced December 2023.

Comments: 13 pages, 8 figures. Accepted for publication in ApJ

ALMA-IMF VIII -- Combination of Interferometric Continuum Images with Single-Dish Surveys and Structural Analysis of Six Protoclusters

Authors: Daniel J. Díaz-González, Roberto Galván-Madrid, Adam Ginsburg, Frédérique Motte, Pierre Dell'Ova, Stan Kurtz, Nichol Cunningham, Amelia M. Stutz, Fabien Louvet, Timea Csengeri, Manuel Fernández-López, Patricio Sanhueza, Thomas Nony, Rudy Rivera-Soto, Rodrigo H. Álvarez-Gutiérrez, Melanie Armante, Melisse Bonfand, Sylvain Bontemps, Antoine Gusdorf, Hong-Li Liu

Abstract: We present the combination of ALMA-IMF and single-dish continuum images from the Mustang-2 Galactic Plane Survey (MGPS90) at 3 millimeters and the Bolocam Galactic Plane Survey (BGPS) at 1 millimeter. Six and ten out of the fiffteen ALMA-IMF fields are combined with MGPS90 and BGPS, respectively. The combination is made via the feathering technique. We used the dendrogram algorithm throughout the… ▽ More We present the combination of ALMA-IMF and single-dish continuum images from the Mustang-2 Galactic Plane Survey (MGPS90) at 3 millimeters and the Bolocam Galactic Plane Survey (BGPS) at 1 millimeter. Six and ten out of the fiffteen ALMA-IMF fields are combined with MGPS90 and BGPS, respectively. The combination is made via the feathering technique. We used the dendrogram algorithm throughout the combined images, and performed further analysis in the six fields with combination in both bands (G012.80, W43-MM1, W43-MM2, W43-MM3, W51-E, W51-IRS2). In these fields, we calculated spectral index maps and used them to separate regions dominated by dust or free-free emission, and then performed further structural analysis. We report the basic physical parameters of the dust-dominated (column densities, masses) and ionized (emission measures, hydrogen ionization photon rates) structures. We also searched for multi-scale relations in the dust-dominated structures across the analyzed fields, finding that the fraction of mass in dendrogram leaves (which we label as "Leaf Mass Eficiency", LME) as a function of molecular gas column density follows a similar trend: a rapid, exponential-like growth, with maximum values approaching 100% in most cases. The observed behaviour of the LME with gas column is tentatively interpreted as an indicator of large star formation activity within the ALMA-IMF protoclusters. W51-E and G012.80 stand out as cases with comparatively large and reduced potential for further star formation, respectively. △ Less

Submitted 20 October, 2023; originally announced October 2023.

Comments: Accepted to The Astrophysical Journal Supplement

ALMA-IMF IX: Catalog and Physical Properties of 315 SiO Outflow Candidates in 15 Massive Protoclusters

Authors: A. P. M. Towner, A. Ginsburg, P. Dell'Ova, A. Gusdorf, S. Bontemps, T. Csengeri, R. Galván-Madrid, F. K. Louvet, F. Motte, P. Sanhueza, A. M. Stutz, J. Bally, T. Baug, H. R. V. Chen, N. Cunningham, M. Fernández-López, H. -L. Liu, X. Lu, T. Nony, M. Valeille-Manet, B. Wu, R. H. Álvarez-Gutiérrez, M. Bonfand, J. Di Francesco, Q. Nguyen-Luong , et al. (2 additional authors not shown)

Abstract: 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,00… ▽ More 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$σ$), and no correlations above 3$σ$ 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. △ Less

Submitted 27 October, 2023; v1 submitted 19 October, 2023; originally announced October 2023.

Comments: 46 pages, 14 figures, 10 tables. This publication has an associated Zenodo entry, which can be found here: https://zenodo.org/records/8350595

A high-resolution radio study of the L1551 IRS 5 and L1551 NE jets

Authors: A. Feeney-Johansson, S. J. D. Purser, T. P. Ray, C. Carrasco-González, A. Rodríguez-Kamenetzky, J. Eislöffel, J. Lim, R. Galván-Madrid, S. Lizano, L. F. Rodríguez, H. Shang, P. Ho, M. Hoare

Abstract: 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… ▽ More 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. △ Less

Submitted 24 July, 2023; v1 submitted 24 July, 2023; originally announced July 2023.

Comments: 13 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 677, A97 (2023)

ALMA-IMF. VII. First release of the full spectral line cubes: Core kinematics traced by DCN J=(3-2)

Authors: N. Cunningham, A. Ginsburg, R. Galván-Madrid, F. Motte, T. Csengeri, A. M. Stutz, M. Fernández-López, R. H. Álvarez-Gutiérrez, M. Armante, T. Baug, M. Bonfand, S. Bontemps, J. Braine, N. Brouillet, G. Busquet, D. J. Díaz-González, J. Di Francesco, A. Gusdorf, F. Herpin, H. Liu, A. López-Sepulcre, F. Louvet, X. Lu, L. Maud, T. Nony , et al. (8 additional authors not shown)

Abstract: 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), prod… ▽ More 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. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: 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

Journal ref: A&A 678, A194 (2023)

Kinematics and stability of high-mass protostellar disk candidates at sub-arcsecond resolution -- Insights from the IRAM NOEMA large program CORE

Authors: Aida Ahmadi, H. Beuther, F. Bosco, C. Gieser, S. Suri, J. C. Mottram, R. Kuiper, Th. Henning, Á. Sánchez-Monge, H. Linz, R. E. Pudritz, D. Semenov, J. M. Winters, T. Möller, M. T. Beltrán, T. Csengeri, R. Galván-Madrid, K. G. Johnston, E. Keto, P. D. Klaassen, S. Leurini, S. N. Longmore, S. L. Lumsden, L. T. Maud, L. Moscadelli , et al. (6 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 3 May, 2023; v1 submitted 28 April, 2023; originally announced May 2023.

Comments: 27 pages, 12 figures, 6 appendices - accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 677, A171 (2023)

Simulated observations of star formation regions: infrared evolution of globally collapsing clouds

Authors: Jesús M. Jáquez-Domínguez, Roberto Galván-Madrid, Jacopo Fritz, Manuel Zamora-Avilés, Peter Camps, Gustavo Bruzual, Maarten Baes, Yuxin Lin, Enrique Vázquez-Semadeni

Abstract: The direct comparison between hydrodynamical simulations and observations is needed to improve the physics included in the former and test biases in the latter. Post-processing radiative transfer and synthetic observations are now the standard way to do this. We report on the first application of the \texttt{SKIRT} radiative transfer code to simulations of a star-forming cloud. The synthetic obser… ▽ More The direct comparison between hydrodynamical simulations and observations is needed to improve the physics included in the former and test biases in the latter. Post-processing radiative transfer and synthetic observations are now the standard way to do this. We report on the first application of the \texttt{SKIRT} radiative transfer code to simulations of a star-forming cloud. The synthetic observations are then analyzed following traditional observational workflows. We find that in the early stages of the simulation, stellar radiation is inefficient in heating dust to the temperatures observed in Galactic clouds, thus the addition of an interstellar radiation field is necessary. The spectral energy distribution of the cloud settles rather quickly after $\sim3$ Myr of evolution from the onset of star formation, but its morphology continues to evolve for $\sim8$ Myr due to the expansion of \textsc{Hii} regions and the respective creation of cavities, filaments, and ridges. Modeling synthetic \textit{Herschel} fluxes with 1- or 2-component modified black bodies underestimates total dust masses by a factor of $\sim2$. Spatially-resolved fitting recovers up to about $70\%$ of the intrinsic value. This ``missing mass'' is located in a very cold dust component with temperatures below $10$ K, which does not contribute appreciably to the far-infrared flux. This effect could bias real observations if such dust exists in large amounts. Finally, we tested observational calibrations of the SFR based on infrared fluxes and concluded that they are in agreement when compared to the intrinsic SFR of the simulation averaged over $\sim100$ Myr. △ Less

Submitted 10 April, 2023; originally announced April 2023.

ALMA-IMF. V. Prestellar and protostellar core populations in the W43 cloud complex

Authors: T. Nony, R. Galvan-Madrid, F. Motte, Y. Pouteau, N. Cunningham, F. Louvet, A. M. Stutz, B. Lefloch, S. Bontemps, N. Brouillet, A. Ginsburg, I. Joncour, F. Herpin, P. Sanhueza, T. Csengeri, A. P. M. Towner, M. Bonfand, M. Fernández-López, T. Baug, L. Bronfman, G. Busquet, J. Di Francesco, A. Gusdorf, X. Lu, F. Olguin , et al. (2 additional authors not shown)

Abstract: 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… ▽ More 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, $α=-1.46_{-0.19}^{+0.12}$, is consistent with the Salpeter slope, and thus the top-heavy form measured for the global CMF, $α=-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. △ Less

Submitted 15 March, 2023; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: 28 pages, 18 figures. Accepted for publication in A&A. Title and appendices updated

Journal ref: A&A 674, A75 (2023)

ALMA-IMF VI -- Investigating the origin of stellar masses: Core mass function evolution in the W43-MM2&MM3 mini-starburst

Authors: Y. Pouteau, F. Motte, T. Nony, M. Gonzalez, I. Joncour, J. -F. Robitaille, G. Busquet, R. Galvan-Madrid, A. Gusdorf, P. Hennebelle, A. Ginsburg, T. Csengeri, P. Sanhueza, P. Dell'Ova, A. M. Stutz, A. P. M. Towner, N. Cunningham, F. Louvet, A. Men'shchikov, M. Fernandez-Lopez, N. Schneider, M. Armante, J. Bally, T. Baug, M. Bonfand , et al. (13 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 20 February, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: Accepted for publication in A&A (december, 5th 2022) after language editing; 25 pages, 14 figures

Journal ref: A&A 674, A76 (2023)

Clustered Formation of Massive Stars within an Ionized Rotating Disk

Authors: Roberto Galván-Madrid, Qizhou Zhang, Andrés Izquierdo, Charles J. Law, Thomas Peters, Eric Keto, Hauyu Baobab Liu, Paul T. P. Ho, Adam Ginsburg, Carlos Carrasco-González

Abstract: We present ALMA observations with a 800 au resolution and radiative-transfer modelling of the inner part ($r\approx6000$ au) of the ionized accretion flow around a compact star cluster in formation at the center of the luminous ultra-compact (UC) HII region G10.6-0.4. We modeled the flow with an ionized Keplerian disk with and without radial motions in its outer part, or with an external Ulrich en… ▽ More We present ALMA observations with a 800 au resolution and radiative-transfer modelling of the inner part ($r\approx6000$ au) of the ionized accretion flow around a compact star cluster in formation at the center of the luminous ultra-compact (UC) HII region G10.6-0.4. We modeled the flow with an ionized Keplerian disk with and without radial motions in its outer part, or with an external Ulrich envelope. The MCMC fits to the data give total stellar masses $M_\star$ from 120 to $200~M_\odot$, with much smaller ionized-gas masses $M_\mathrm{ion-gas} = 0.2$ to $0.25~M_\odot$. The stellar mass is distributed within the gravitational radius $R_g\approx 1000$ to 1500 au, where the ionized gas is bound. The viewing inclination angle from the face-on orientation is $i = 49$ to $56~°$. Radial motions at radii $r > R_g$ converge to $v_{r,0} \approx 8.7$ km/s, or about the speed of sound of ionized gas, indicating that this gas is marginally unbound at most. From additional constraints on the ionizing-photon rate and far-IR luminosity of the region, we conclude that the stellar cluster consists of a few massive stars with $M_\mathrm{star} = 32$ to $60~M_\odot$, or one star in this range of masses accompanied by a population of lower-mass stars. Any active accretion of ionized gas onto the massive (proto)stars is residual. The inferred cluster density is very large, comparable to that reported at similar scales in the Galactic Center. Stellar interactions are likely to occur within the next Myr. △ Less

Submitted 11 December, 2022; originally announced December 2022.

Comments: Accepted to ApJ Letters, December 7 2022

Time-Variable Jet Ejections from RW Aur A, RY Tau and DG Tau

Authors: Michihiro Takami, Hans Moritz Guenther, P. Christian Schneider, Tracy L. Beck, Jennifer L. Karr, Youichi Ohyama, Roberto Galvan-Madrid, Taichi Uyama, Marc White, Konstantin Grankin, Deirdre Coffey, Chun-Fan Liu, Misato Fukagawa, Nadine Manset, Wen-Ping Chen, Tae-Soo Pyo, Hsien Shang, Thomas P. Ray, Masaaki Otsuka, Mei-Yin Chou

Abstract: We present Gemini-NIFS, VLT-SINFONI and Keck-OSIRIS observations of near-infrared [Fe II] emission associated with the well-studied jets from three active T Tauri stars; RW Aur A, RY Tau and DG Tau taken from 2012-2021. We primarily covered the redshifted jet from RW Aur A, and the blueshifted jets from RY Tau and DG Tau, to investigate long-term time variabilities potentially related to the activ… ▽ More We present Gemini-NIFS, VLT-SINFONI and Keck-OSIRIS observations of near-infrared [Fe II] emission associated with the well-studied jets from three active T Tauri stars; RW Aur A, RY Tau and DG Tau taken from 2012-2021. We primarily covered the redshifted jet from RW Aur A, and the blueshifted jets from RY Tau and DG Tau, to investigate long-term time variabilities potentially related to the activities of mass accretion and/or the stellar magnetic fields. All of these jets consist of several moving knots with tangential velocities of 70-240 km s-1, ejected from the star with different velocities and at irregular time intervals. Via comparison with literature, we identify significant differences in tangential velocities for the DG Tau jet between 1985-2008 and 2008-2021. The sizes of the individual knots appear to increase with time, and in turn, their peak brightnesses in the 1.644-micron emission decreased up to a factor of ~30 during the epochs of our observations. A variety of the decay timescales measured in the [Fe II] 1.644 micron emission can be attributed to different pre-shock conditions if the moving knots are unresolved shocks. However, our data do not exclude the possibility that these knots are due to non-uniform density/temperature distributions with another heating mechanism, or in some cases due to stationary shocks without proper motions. Spatially resolved observations of these knots with significantly higher angular resolutions are necessary to better understand their physical nature. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: 24 pages, 11 figures, 8 tables, accepted for ApJS

ALMA-IMF IV -- A comparative study of the main hot cores in W43-MM1: detection, temperature and molecular composition

Authors: N. Brouillet, D. Despois, J. Molet, T. Nony, F. Motte, A. Gusdorf, F. Louvet, S. Bontemps, F. Herpin, M. Bonfand, T. Csengeri, A. Ginsburg, N. Cunningham, R. Galvan-Madrid, L. Maud, G. Busquet, L. Bronfman, M. Fernandez-Lopez, D. L. Jeff, B. Lefloch, Y. Pouteau, P. Sanhueza, A. M. Stutz, M. Valeille-Manet

Abstract: 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… ▽ More 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. △ Less

Submitted 7 July, 2022; originally announced July 2022.

Comments: 19 pages, 17 figures. Accepted for publication in A&A

Journal ref: A&A 665, A140 (2022)

Resolving the collimation zone of an intermediate-mass protostar

Authors: Adriana R. Rodríguez-Kamenetzky, Carlos Carrasco-González, Luis Felipe Rodríguez Jorge, Tom P. Ray, Alberto Sanna, Luca Moscadelli, Melvin Hoare, Roberto Galván-Madrid, Hsien Shang, Susana Lizano, Jochen Eislöffel, Jeremy Lim, José M. Torrelles, Paul Ho, Anton Feeney-Johansson

Abstract: We report new VLA and e-MERLIN high resolution and sensitivity images of the Triple Radio continuum Source in the Serpens star forming region. These observations allowed us to perform a deep multi-frequency, multi-epoch study by exploring the innermost regions (<~100 au) of an intermediate-mass YSO for the first time, with a physical resolution of ~15 au. The kinematic analysis of knots recently e… ▽ More We report new VLA and e-MERLIN high resolution and sensitivity images of the Triple Radio continuum Source in the Serpens star forming region. These observations allowed us to perform a deep multi-frequency, multi-epoch study by exploring the innermost regions (<~100 au) of an intermediate-mass YSO for the first time, with a physical resolution of ~15 au. The kinematic analysis of knots recently ejected by the protostar indicates that the jet is undergoing episodic variations in velocity. In addition, our multi-frequency images reveal striking characteristics, e.g., a highly collimated ionized stream that would be launched at a radial distance of ~0.4 au from the protostar, and a narrow (~28 au wide) ionized cavity that would be excited by the interaction of a wide-angle component with the surrounding toroid of infalling material. In light of these results, we propose the scenario in which both a highly-collimated jet and a wide-angle wind coexist to be the most plausible to explain our observations, either launched by the X-wind or X- plus Disk-wind mechanism. △ Less

Submitted 20 May, 2022; originally announced May 2022.

Comments: 13 pages, 3 figures. Accepted for publication in APJL

MSC Class: 85-02 ACM Class: J.2

ALMA-IMF III -- Investigating the origin of stellar masses: Top-heavy core mass function in the W43-MM2&MM3 mini-starburst

Authors: Y. Pouteau, F. Motte, T. Nony, R. Galván-Madrid, A. Men'shchikov, S. Bontemps, J. -F. Robitaille, F. Louvet, A. Ginsburg, F. Herpin, A. López-Sepulcre, P. Dell'Ova, A. Gusdorf, P. Sanhueza, A. M. Stutz, N. Brouillet, B. Thomasson, M. Armante, T. Baug, G. Busquet, T. Csengeri, N. Cunningham, M. Fernández-López, H. -L. Liu, F. Olguin , et al. (13 additional authors not shown)

Abstract: 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 p… ▽ More 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. △ Less

Submitted 28 April, 2022; v1 submitted 7 March, 2022; originally announced March 2022.

Comments: 28 pages, 16 figures, 6 tables, A&A accepted on 02/08/22, revised following the language referee report

Journal ref: A&A 664, A26 (2022)

ALMA-IMF II -- investigating the origin of stellar masses: Continuum Images and Data Processing

Authors: A. Ginsburg, T. Csengeri, R. Galván-Madrid, N. Cunningham, R. H. Álvarez-Gutiérrez, T. Baug, M. Bonfand, S. Bontemps, G. Busquet, D. J. Díaz-González, M. Fernández-López, A. Guzmán, F. Herpin, H. Liu, A. López-Sepulcre, F. Louvet, L. Maud, F. Motte, F. Nakamura, T. Nony, F. A. Olguin, Y. Pouteau, P. Sanhueza, A. M. Stutz, A. P. M. Towner , et al. (27 additional authors not shown)

Abstract: 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 Ar… ▽ More 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. △ Less

Submitted 13 May, 2023; v1 submitted 15 December, 2021; originally announced December 2021.

Comments: Data released on Zenodo at https://zenodo.org/record/5702966. Accepted for publication in A&A

Journal ref: A&A 662, A9 (2022)

ALMA-IMF I -- Investigating the origin of stellar masses: Introduction to the Large Program and first results

Authors: F. Motte, S. Bontemps, T. Csengeri, Y. Pouteau, F. Louvet, A. M. Stutz, N. Cunningham, A. López-Sepulcre, N. Brouillet, R. Galván-Madrid, A. Ginsburg, L. Maud, A. Men'shchikov, F. Nakamura, T. Nony, P. Sanhueza, R. H. Álvarez-Gutiérrez, M. Armante, T. Baug, M. Bonfand, G. Busquet, E. Chapillon, D. Díaz-González, M. Fernández-López, A. E. Guzmán , et al. (39 additional authors not shown)

Abstract: 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 2… ▽ More 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. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: 27 pages, 9 figures, accepted for publication by A&A

Journal ref: A&A 662, A8 (2022)

Clustered star formation at early evolutionary stages. Physical and chemical analysis of the young star-forming regions ISOSS J22478+6357 and ISOSS J23053+5953

Authors: C. Gieser, H. Beuther, D. Semenov, S. Suri, J. D. Soler, H. Linz, J. Syed, Th. Henning, S. Feng, T. Möller, A. Palau, J. M. Winters, M. T. Beltrán, R. Kuiper, L. Moscadelli, P. Klaassen, J. S. Urquhart, T. Peters, S. N. Longmore, Á. Sánchez-Monge, R. Galván-Madrid, R. E. Pudritz, K. G. Johnston

Abstract: 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 t… ▽ More 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. △ Less

Submitted 14 October, 2021; v1 submitted 5 October, 2021; originally announced October 2021.

Comments: 34 pages, 17 figures, accepted for publication in A&A

Journal ref: A&A 657, A3 (2022)

A Core Mass Function Indistinguishable from the Salpeter Stellar Initial Mass Function Using 1000 au Resolution ALMA Observations

Authors: Genaro Suárez, Roberto Galván-Madrid, Luis Aguilar, Adam Ginsburg, Sundar Srinivasan, Hauyu Baobab Liu, Carlos G. Román-Zúñiga

Abstract: We present the core mass function (CMF) of the massive star-forming clump G33.92+0.11 using 1.3 mm observations obtained with the Atacama Large Millimeter/submillimeter Array (ALMA). With a resolution of 1000 au, this is one of the highest resolution CMF measurements to date. The CMF is corrected by flux and number incompleteness to obtain a sample that is complete for gas masses… ▽ More We present the core mass function (CMF) of the massive star-forming clump G33.92+0.11 using 1.3 mm observations obtained with the Atacama Large Millimeter/submillimeter Array (ALMA). With a resolution of 1000 au, this is one of the highest resolution CMF measurements to date. The CMF is corrected by flux and number incompleteness to obtain a sample that is complete for gas masses $M\gtrsim2.0\ M_\odot$. The resulting CMF is well represented by a power-law function ($dN/d\log M\propto M^Γ$), whose slope is determined using two different approaches: $i)$ by least-squares fitting of power-law functions to the flux- and number-corrected CMF, and $ii)$ by comparing the observed CMF to simulated samples with similar incompleteness. We provide a prescription to quantify and correct a flattening bias affecting the slope fits in the first approach, which is caused by small-sample or edge effects when the data is represented by either classical histograms or a kernel density estimate, respectively. The resulting slopes from both approaches are in good agreement each other, with $Γ=-1.11_{-0.11}^{+0.12}$ being our adopted value. Although this slope appears to be slightly flatter than the Salpeter slope $Γ=-1.35$ for the stellar initial mass function (IMF), we find from Monte Carlo simulations that the CMF in G33.92+0.11 is statistically indistinguishable from the Salpeter representation of the stellar IMF. Our results are consistent with the idea that the form of the IMF is inherited from the CMF, at least at high masses and when the latter is observed at high-enough resolution. △ Less

Submitted 29 July, 2021; originally announced July 2021.

Comments: Accepted for publication in ApJ

Zooming into the Collimation Zone in a Massive Protostellar Jet

Authors: Carlos Carrasco-González, Alberto Sanna, Adriana Rodríguez-Kamenetzky, Luca Moscadelli, Melvin Hoare, José M. Torrelles, Roberto Galván-Madrid, Andrés F. Izquierdo

Abstract: Protostellar jets have a fundamental role at the earliest evolution of protostars of all masses. In the case of low-mass (<8 Msun) protostars, strong observational evidence exists that the launching and collimation is due to the X- and/or disk-wind mechanisms. In these models, it is the protostar/disk system that creates all the necessary conditions to launch and collimate the jets near the protos… ▽ More Protostellar jets have a fundamental role at the earliest evolution of protostars of all masses. In the case of low-mass (<8 Msun) protostars, strong observational evidence exists that the launching and collimation is due to the X- and/or disk-wind mechanisms. In these models, it is the protostar/disk system that creates all the necessary conditions to launch and collimate the jets near the protostar via strong magnetic fields. The origin of jets from more massive protostars has been investigated much less, in part because of the difficulty of resolving the collimation zone in these more distant objects. Here we present the highest angular resolution observations of a jet powered by a massive protostar, the Cep A HW2 radio jet. We imaged the radio emission at projected distances of only ~20 au from the protostar, resolving the innermost 100 au of a massive protostellar jet for the first time. The morphology of the radio jet emission in this massive object is very different than what is usually observed in jets from low-mass protostars. We found that the outflowing material in HW2 has two components: a wide-angle wind launched from the protostar/disk system, and a highly collimated jet starting at 20-30 au from the protostar. We discuss two possible scenarios: an extension of the classical disk-wind to a massive protostar, or external collimation of a wide-angle wind. These results have important consequences for our understanding of how stars of different masses are formed. △ Less

Submitted 2 June, 2021; originally announced June 2021.

Comments: Accepted for publication in ApJ Letters

The physical and chemical structure of high-mass star-forming regions. Unraveling chemical complexity with the NOEMA large program "CORE"

Authors: C. Gieser, H. Beuther, D. Semenov, A. Ahmadi, S. Suri, T. Möller, M. T. Beltran, P. Klaassen, Q. Zhang, J. S. Urquhart, Th. Henning, S. Feng, R. Galván-Madrid, V. de Souza Magalhães, L. Moscadelli, S. Longmore, S. Leurini, R. Kuiper, T. Peters, K. M. Menten, T. Csengeri, G. Fuller, F. Wyrowski, S. Lumsden, Á. Sánchez-Monge , et al. (8 additional authors not shown)

Abstract: 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 profil… ▽ More 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. △ Less

Submitted 23 February, 2021; originally announced February 2021.

Comments: 29 pages, 14 figures, accepted for publication in A&A

Journal ref: A&A 648, A66 (2021)

Multi-scale view of star formation in IRAS 21078+5211: From clump fragmentation to disk wind

Authors: L. Moscadelli, H. Beuther, A. Ahmadi, C. Gieser, F. Massi, R. Cesaroni, Á. Sánchez-Monge, F. Bacciotti, M. T. Beltrán, T. Csengeri, R. Galván-Madrid, Th. Henning, P. D. Klaassen, R. Kuiper, S. Leurini, S. N. Longmore, L. T. Maud, T. Möller, A. Palau, T. Peters, R. E. Pudritz, A. Sanna, D. Semenov, J. S. Urquhart, J. M. Winters , et al. (1 additional authors not shown)

Abstract: 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 deri… ▽ More 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. △ Less

Submitted 9 February, 2021; originally announced February 2021.

Comments: 25 pages, 16 figures

Journal ref: A&A 647, A114 (2021)

Sub-arcsecond Imaging of the Complex Organic Chemistry in Massive Star-forming Region G10.6-0.4

Authors: Charles J. Law, Qizhou Zhang, Karin I. Öberg, Roberto Galván-Madrid, Eric Keto, Hauyu Baobab Liu, Paul T. P. Ho

Abstract: Massive star-forming regions exhibit an extremely rich and diverse chemistry, which in principle provides a wealth of molecular probes, as well as laboratories for interstellar prebiotic chemistry. Since the chemical structure of these sources displays substantial spatial variation among species on small scales (${\lesssim}10^4$ au), high angular resolution observations are needed to connect chemi… ▽ More Massive star-forming regions exhibit an extremely rich and diverse chemistry, which in principle provides a wealth of molecular probes, as well as laboratories for interstellar prebiotic chemistry. Since the chemical structure of these sources displays substantial spatial variation among species on small scales (${\lesssim}10^4$ au), high angular resolution observations are needed to connect chemical structures to local environments and inform astrochemical models of massive star formation. To address this, we present ALMA 1.3 mm observations toward OB cluster-forming region G10.6-0.4 (hereafter "G10.6") at a resolution of 0.14$^{\prime\prime}$ (700 au). We find highly-structured emission from complex organic molecules (COMs) throughout the central 20,000 au, including two hot molecular cores and several shells or filaments. We present spatially-resolved rotational temperature and column density maps for a large sample of COMs and warm gas tracers. These maps reveal a range of gas substructure in both O- and N-bearing species. We identify several spatial correlations that can be explained by existing models of COM formation, including NH$_2$CHO/HNCO and CH$_3$OCHO/CH$_3$OCH$_3$, but also observe unexpected distributions and correlations which suggest that our current understanding of COM formation is far from complete. Importantly, complex chemistry is observed throughout G10.6, rather than being confined to hot cores. The COM composition appears to be different in the cores compared to the more extended structures, which illustrates the importance of high spatial resolution observations of molecular gas in elucidating the physical and chemical processes associated with massive star formation. △ Less

Submitted 19 January, 2021; originally announced January 2021.

Comments: 41 pages, 25 figures, accepted for publication in ApJ

Formation of the Hub-Filament System G33.92+0.11: Local Interplay between Gravity, Velocity, and Magnetic Field

Authors: Jia-Wei Wang, Patrick M. Koch, Roberto Galván-Madrid, Shih-Ping Lai, Hauyu Baobab Liu, Sheng-Jun Lin, Kate Pattle

Abstract: The formation of filaments in molecular clouds is an important process in star formation. Hub-filament systems (HFSs) are a transition stage connecting parsec-scale filaments and proto-clusters. Understanding the origin of HFSs is crucial to reveal how star formation proceeds from clouds to cores. Here, we report JCMT POL-2 850 $μ$m polarization and IRAM 30-m C$^{18}$O (2-1) line observations towa… ▽ More The formation of filaments in molecular clouds is an important process in star formation. Hub-filament systems (HFSs) are a transition stage connecting parsec-scale filaments and proto-clusters. Understanding the origin of HFSs is crucial to reveal how star formation proceeds from clouds to cores. Here, we report JCMT POL-2 850 $μ$m polarization and IRAM 30-m C$^{18}$O (2-1) line observations toward the massive HFS G33.92+0.11. The 850 $μ$m continuum map reveals four major filaments converging to the center of G33.92+0.11 with numerous short filaments connecting to the major filaments at local intensity peaks. We estimate the local orientations of filaments, magnetic field, gravity, and velocity gradients from observations, and we examine their correlations based on their local properties. In the high-density areas, our analysis shows that the filaments tend to align with the magnetic field and local gravity. In the low-density areas, we find that the local velocity gradients tend to be perpendicular to both the magnetic field and local gravity, although the filaments still tend to align with local gravity. A global virial analysis suggests that the gravitational energy overall dominates the magnetic and kinematic energy. Combining local and global aspects, we conclude that the formation of G33.92+0.11 is predominantly driven by gravity, dragging and aligning the major filaments and magnetic field on the way to the inner dense center. Traced by local velocity gradients in the outer diffuse areas, ambient gas might be accreted onto the major filaments directly or via the short filaments. △ Less

Submitted 3 November, 2020; originally announced November 2020.

Comments: 32 pages, 23 figures, accepted for publication in ApJ

Time-Variable Radio Recombination Line Emission in W49A

Authors: C. G. De Pree, D. J. Wilner, L. E. Kristensen, R. Galván-Madrid, W. M. Goss, R. S. Klessen, M. -M. Mac Low, T. Peters, A. Robinson, S. Sloman, M. Rao

Abstract: 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 elonga… ▽ More 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 $Δ$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. △ Less

Submitted 9 October, 2020; originally announced October 2020.

Comments: 18 pages, 7 figures, Accepted for publication in The Astronomical Journal

Possible Time Correlation Between Jet Ejection and Mass Accretion for RW Aur A

Authors: Michihiro Takami, Tracy L. Beck, P. Christian Schneider, Hans Moritz Guenther, Marc White, Konstantin Grankin, Jennifer L. Karr, Youichi Ohyama, Deirdre Coffey, Hauyu Baobab Liu, Roberto Galvan-Madrid, Chun-Fan Liu, Misato Fukagawa, Nadine Manset, Wen-Ping Chen, Tae-Soo Pyo, Hsien Shang, Thomas P. Ray, Masaaki Otsuka, Mei-Yin Chou

Abstract: For the active T-Taur star RW Aur A we have performed long-term (~10 yr) monitoring observations of (1) jet imaging in the [Fe II] 1.644-micron emission line using Gemini-NIFS and VLT-SINFONI; (2) optical high-resolution spectroscopy using CFHT-ESPaDOnS; and (3) V-band photometry using the CrAO 1.25-m telescope and AAVSO. The latter two observations confirm the correlation of time variabilities be… ▽ More For the active T-Taur star RW Aur A we have performed long-term (~10 yr) monitoring observations of (1) jet imaging in the [Fe II] 1.644-micron emission line using Gemini-NIFS and VLT-SINFONI; (2) optical high-resolution spectroscopy using CFHT-ESPaDOnS; and (3) V-band photometry using the CrAO 1.25-m telescope and AAVSO. The latter two observations confirm the correlation of time variabilities between (A) the Ca II 8542 A and O I 7772 A line profiles associated with magnetospheric accretion, and (B) optical continuum fluxes. The jet images and their proper motions show that four knot ejections occurred at the star over the past ~15 years with an irregular interval of 2-6 years. The time scale and irregularity of these intervals are similar to those of the dimming events seen in the optical photometry data. Our observations show a possible link between remarkable (Delta_V < -1 mag.) photometric rises and jet knot ejections. Observations over another few years may confirm or reject this trend. If confirmed, this would imply that the location of the jet launching region is very close to the star (r <<0.1 au) as predicted by some jet launching models. Such a conclusion would be crucial for understanding disk evolution within a few au of the star, and therefore possible ongoing planet formation at these radii. △ Less

Submitted 31 July, 2020; originally announced July 2020.

Comments: 14 pages, 5 figures, Accepted by Astrophysical Journal

Recombination Lines and Molecular Gas from Hypercompact HII regions in W51 A

Authors: Rudy Rivera-Soto, Roberto Galván-Madrid, Adam Ginsburg, Stan Kurtz

Abstract: We present a detailed characterization of the population of compact radio-continuum sources in W51 A using subarcsecond VLA and ALMA observations. We analyzed their 2-cm continuum, the recombination lines (RL's) H77$α$ and H30$α$, and the lines of $\rm H_{2}CO(3_{0,3}-2_{0,2})$, $\rm H_{2}CO(3_{2,1}-2_{2,0})$, and $\rm SO(6_{5}-5_{4})$. We derive diameters for 10/20 sources in the range… ▽ More We present a detailed characterization of the population of compact radio-continuum sources in W51 A using subarcsecond VLA and ALMA observations. We analyzed their 2-cm continuum, the recombination lines (RL's) H77$α$ and H30$α$, and the lines of $\rm H_{2}CO(3_{0,3}-2_{0,2})$, $\rm H_{2}CO(3_{2,1}-2_{2,0})$, and $\rm SO(6_{5}-5_{4})$. We derive diameters for 10/20 sources in the range $D \sim 10^{-3}$ to $\sim 10^{-2}$ pc, thus placing them in the regime of hypercompact HII regions (HC HII's). Their continuum-derived electron densities are in the range $n_{\rm e} \sim 10^4$ to $10^5$ cm$^{-3}$, lower than typically considered for HC HII's. We combined the RL measurements and independently derived $n_{\rm e}$, finding the same range of values but significant offsets for individual measurements between the two methods. We found that most of the sources in our sample are ionized by early B-type stars, and a comparison of $n_{\rm e}$ vs $D$ shows that they follow the inverse relation previously derived for ultracompact (UC) and compact HII's. When determined, the ionized-gas kinematics is always (7/7) indicative of outflow. Similarly, 5 and 3 out of the 8 HC HII's still embedded in a compact core show evidence for expansion and infall motions in the molecular gas, respectively. We hypothesize that there could be two different types of $hypercompact$ ($D< 0.05$ pc) HII regions: those that essentially are smaller, expanding UC HII's; and those that are also $hyperdense$ ($n_{\rm e} > 10^6$ cm$^{-3}$), probably associated with O-type stars in a specific stage of their formation or early life. △ Less

Submitted 23 July, 2020; originally announced July 2020.

Comments: The Astrophysical Journal, in press. 23 pages, 9 figures

The MUSTANG-2 Galactic Plane Survey (MGPS90) pilot

Authors: Adam Ginsburg, L. D. Anderson, Simon Dicker, Charles Romero, Brian Svoboda, Mark Devlin, Roberto Galván-Madrid, Remy Indebetouw, Hauyu Baobab Liu, Brian Mason, Tony Mroczkowski, W. P. Armentrout, John Bally, Crystal Brogan, Natalie Butterfield, Todd R. Hunter, Erik D. Reese, Erik Rosolowsky, Craig Sarazin, Yancy Shirley, Jonathan Sievers, Sara Stanchfield

Abstract: We report the results of a pilot program for a Green Bank Telescope (GBT) MUSTANG Galactic Plane survey at 3 mm (90 GHz), MGPS90. The survey achieves a typical $1σ$ depth of $1-2$ mJy beam$^{-1}$ with a 9" beam. We describe the survey parameters, quality assessment process, cataloging, and comparison with other data sets. We have identified 709 sources over seven observed fields selecting some of… ▽ More We report the results of a pilot program for a Green Bank Telescope (GBT) MUSTANG Galactic Plane survey at 3 mm (90 GHz), MGPS90. The survey achieves a typical $1σ$ depth of $1-2$ mJy beam$^{-1}$ with a 9" beam. We describe the survey parameters, quality assessment process, cataloging, and comparison with other data sets. We have identified 709 sources over seven observed fields selecting some of the most prominent millimeter-bright regions between $0°< \ell < 50°$ (total area $\approx 7.5 °^2$). The majority of these sources have counterparts at other wavelengths. By applying flux selection criteria to these sources, we successfully recovered several known hypercompact HII (HCHII) regions, but did not confirm any new ones. We identify 126 sources that have mm-wavelength counterparts but do not have cm-wavelength counterparts and are therefore candidate HCHII regions; of these, 10 are morphologically compact and are strong candidates for new HCHII regions. Given the limited number of candidates in the extended area in this survey compared to the relatively large numbers seen in protoclusters W51 and W49, it appears that most HCHII regions exist within dense protoclusters. Comparing the counts of HCHII to ultracompact HII (UCHII) regions, we infer the HCHII region lifetime is 16-46% that of the UCHII region lifetime. We additionally separated the 3 mm emission into dust and free-free emission by comparing with archival 870 $μ$m and 20 cm data. In the selected pilot fields, most ($\gtrsim80$%) of the 3 mm emission comes from plasma, either through free-free or synchrotron emission. △ Less

Submitted 22 April, 2020; v1 submitted 20 April, 2020; originally announced April 2020.

Comments: Accepted to ApJS. Data publicly released, with links in paper. Reposted to fix missing figure labels

A Massive Young Runaway Star in W49 North

Authors: Luis F. Rodriguez, Roberto Galvan-Madrid, Joel Sanchez-Bermudez, Christopher G. De Pree

Abstract: We analyzed high angular resolution 45.5 GHz images of the W49 North massive star forming region obtained in 1998 and 2016 with the Very Large Array. Most of the ultracompact HII regions show no detectable changes over the time interval of the observations. However, subcomponents B1, B2, G2a and G2c have increased its peak flux densities by values in the range of 3.8 to 21.4 \%. Most interestingly… ▽ More We analyzed high angular resolution 45.5 GHz images of the W49 North massive star forming region obtained in 1998 and 2016 with the Very Large Array. Most of the ultracompact HII regions show no detectable changes over the time interval of the observations. However, subcomponents B1, B2, G2a and G2c have increased its peak flux densities by values in the range of 3.8 to 21.4 \%. Most interestingly, the cometary region C clearly shows proper motions that at the distance of the region are equivalent to a velocity of 76$\pm$6 km s$^{-1}$ in the plane of the sky. We interpret this region as the ionized bowshock produced by a runaway O6 ZAMS star that was ejected from the eastern edge of Welch's ring about 6,400 years ago. △ Less

Submitted 7 February, 2020; originally announced February 2020.

Comments: 11 pages, 6 figures,to appear in ApJ

Modeling the accretion disk around the high-mass protostar GGD 27-MM1

Authors: N. Añez-Lopez, M. Osorio, G. Busquet, J. M. Girart, E. Macıas, C. Carrasco-Gonzalez, S. Curiel, R. Estalella, M. Fernandez-Lopez, R. Galvan-Madrid, J. Kwon, J. M. Torrelles

Abstract: 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 r… ▽ More 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). △ Less

Submitted 27 November, 2019; originally announced November 2019.

Comments: 20 pages, 15 figures

Chemical complexity in high-mass star formation: An observational and modeling case study of the AFGL 2591 VLA 3 hot core

Authors: C. Gieser, D. Semenov, H. Beuther, A. Ahmadi, J. C. Mottram, Th. Henning, M. Beltran, L. T. Maud, F. Bosco, S. Leurini, T. Peters, P. Klaassen, R. Kuiper, S. Feng, J. S. Urquhart, L. Moscadelli, T. Csengeri, S. Lumsden, J. M. Winters, S. Suri, Q. Zhang, R. Pudritz, A. Palau, K. M. Menten, R. Galvan-Madrid , et al. (8 additional authors not shown)

Abstract: 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 observ… ▽ More 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. △ Less

Submitted 11 October, 2019; originally announced October 2019.

Comments: 22 pages, 13 figures, accepted for publication in A&A

Journal ref: A&A 631, A142 (2019)

The Radial Distribution of Dust Particles in the HL Tau Disk from ALMA and VLA Observations

Authors: Carlos Carrasco-González, Anibal Sierra, Mario Flock, Zhaohuan Zhu, Thomas Henning, Claire Chandler, Roberto Galván-Madrid, Enrique Macías, Guillem Anglada, Hendrik Linz, Mayra Osorio, Luis F. Rodríguez, Leonardo Testi, José M. Torrelles, Laura Pérez, Yao Liu

Abstract: Understanding planet formation requires to discern how dust grows in protoplanetary disks. An important parameter to measure in disks is the maximum dust grain size present. This is usually estimated through measurements of the dust opacity at different millimeter wavelengths assuming optically thin emission and dust opacity dominated by absorption. However, ALMA observations have shown that these… ▽ More Understanding planet formation requires to discern how dust grows in protoplanetary disks. An important parameter to measure in disks is the maximum dust grain size present. This is usually estimated through measurements of the dust opacity at different millimeter wavelengths assuming optically thin emission and dust opacity dominated by absorption. However, ALMA observations have shown that these assumptions might not be correct in the case of protoplanetary disks, leading to overestimation of particle sizes and to underestimation of the disk's mass. Here, we present an analysis of high quality ALMA and VLA images of the HL Tau protoplanetary disk, covering a wide range of wavelengths, from 0.8 mm to 1 cm, and with a physical resolution of $\sim$7.35 au. We describe a procedure to analyze a set of millimeter images without any assumption about the optical depth of the emission, and including the effects of absorption and scattering in the dust opacity. This procedure allows us to obtain the dust temperature, the dust surface density and the maximum particle size at each radius. In the HL Tau disk, we found that particles have already grown up to a few millimeters in size. We detect differences in the dust properties between dark and bright rings, with dark rings containing low dust density and small dust particles. Different features in the HL Tau disk seem to have different origins. Planet-disk interactions can explain substructure at the external half of the disk, but the internal rings seem to be associated to the presence of snow lines of several molecules. △ Less

Submitted 24 September, 2019; v1 submitted 19 August, 2019; originally announced August 2019.

Diagnosing 0.1-10 au scale morphology of the FU\,Ori disk using ALMA and VLTI/GRAVITY

Authors: Hauyu Baobab Liu, Antoine Mérand, Joel D. Green, Sebastián Pérez, Antonio S. Hales, Yao-Lun Yang, Michael M. Dunham, Yasuhiro Hasegawa, Thomas Henning, Roberto Galván-Madrid, Ágnes Kóspál, Michihiro Takami, Eduard I. Vorobyov, Zhaohuan Zhu

Abstract: We report new Atacama Large Millimeter/submillimeter Array Band 3 (86-100 GHz; $\sim$80 mas angular resolution) and Band 4 (146-160 GHz; $\sim$50 mas angular resolution) observations of the dust continuum emission towards the archetypal and ongoing accretion burst young stellar object FU Ori, which simultaneously covered its companion, FU Ori S. In addition, we present near-infrared (2-2.45 $μ$m)… ▽ More We report new Atacama Large Millimeter/submillimeter Array Band 3 (86-100 GHz; $\sim$80 mas angular resolution) and Band 4 (146-160 GHz; $\sim$50 mas angular resolution) observations of the dust continuum emission towards the archetypal and ongoing accretion burst young stellar object FU Ori, which simultaneously covered its companion, FU Ori S. In addition, we present near-infrared (2-2.45 $μ$m) observations of FU Ori taken with the General Relativity Analysis via VLT InTerferometrY (GRAVITY; $\sim$1 mas angular resolution) instrument on the Very Large Telescope Interferometer (VLTI). We find that the emission in both FU Ori and FU Ori S at (sub)millimeter and near infrared bands is dominated by structures inward of $\sim$10 au radii. We detected closure phases close to zero from FU Ori with VLTI/GRAVITY, which indicate the source is approximately centrally symmetric and therefore is likely viewed nearly face-on. Our simple model to fit the GRAVITY data shows that the inner 0.4 au radii of the FU Ori disk has a triangular spectral shape at 2-2.45 $μ$m, which is consistent with the H$_{2}$O and CO absorption features in a $\dot{M}\sim$10$^{-4}$ $M_{\odot}\,yr^{-1}$, viscously heated accretion disk. At larger ($\sim$0.4-10 au) radii, our analysis shows that viscous heating may also explain the observed (sub)millimeter and centimeter spectral energy distribution when we assume a constant, $\sim$10$^{-4}$ $M_{\odot}\,yr^{-1}$ mass inflow rate in this region. This explains how the inner 0.4 au disk is replenished with mass at a modest rate, such that it neither depletes nor accumulates significant masses over its short dynamic timescale. Finally, we tentatively detect evidence of vertical dust settling in the inner 10 au of the FU Ori disk, but confirmation requires more complete spectral sampling in the centimeter bands. △ Less

Submitted 8 August, 2019; originally announced August 2019.

Comments: 16 pages, 9 figures; Accepted to ApJ

Substructures in the Keplerian disc around the O-type (proto)star G17.64+0.16

Authors: L. T. Maud, R. Cesaroni, M. S. N. Kumar, V. M. Rivilla, A. Ginsburg, P. D. Klaassen, D. Harsono, A. Sanchez-Monge, A. Ahmadi, V. Allen, M. T. Beltran, H. Beuther, R. Galvan-Madrid, C. Goddi, M. G. Hoare, M. R. Hogerheijde, K. G. Johnston, R. Kuiper, L. Moscadelli, T. Peters, L. Testi, F. F. S. van der Tak, W. J. de Wit

Abstract: 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… ▽ More 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. △ Less

Submitted 15 June, 2019; originally announced June 2019.

Comments: 9 Total, inc references, appendix, long author list, 4 figures (inc appendix)

Journal ref: A&A 627, L6 (2019)

IRAS23385+6053: An embedded massive cluster in the making

Authors: R. Cesaroni, H. Beuther, A. Ahmadi, M. T. Beltran, T. Csengeri, R. Galvan-Madrid, C. Gieser, T. Henning, K. G. Johnston, P. D. Klaassen, R. Kuiper, S. Leurini, H. Linz, S. Longmore, S. L. Lumsden, L. T. Maud, L. Moscadelli, J. C. Mottram, A. Palau, T. Peters, R. E. Pudritz, A. Sanchez-Monge, P. Schilke, D. Semenov, S. Suri , et al. (4 additional authors not shown)

Abstract: 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 additi… ▽ More 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. △ Less

Submitted 27 May, 2019; originally announced May 2019.

Journal ref: A&A 627, A68 (2019)

ALMA survey of Class II protoplanetary disks in Corona Australis: a young region with low disk masses

Authors: P. Cazzoletti, C. F. Manara, H. B. Liu, E. F. van Dishoeck, S. Facchini, J. M. Alcalà, M. Ansdell, L. Testi, J. P. Williams, C. Carrasco-González, R. Dong, J. Forbrich, M. Fukagawa, R. Galván-Madrid, N. Hirano, M. Hogerheijde, Y. Hasegawa, T. Muto, P. Pinilla, M. Takami, M. Tamura, M. Tazzari, J. P. Wisniewski

Abstract: 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 obs… ▽ More 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. △ Less

Submitted 4 April, 2019; originally announced April 2019.

Comments: 20 pages, 8 Figures, accepted for publication on A&A. The abstract has been shortened in order to fit arXiv requirements

Journal ref: A&A 626, A11 (2019)

Unveiling a cluster of protostellar disks around the massive protostar GGD27 MM1

Authors: G. Busquet, J. M. Girart, R. Estalella, M. Fernández-López, R. Galván-Madrid, G. Anglada, C. Carrasco-González, N. Añez-López, S. Curiel, M. Osorio, L. F. Rodríguez, J. M. Torrelles

Abstract: 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 indicat… ▽ More 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. △ Less

Submitted 20 February, 2019; originally announced February 2019.

Comments: Accepted for publication in Astronomy & Astrophysics Letters

Journal ref: A&A 623, L8 (2019)

A VLBI study of the wind-wind collision region in the massive multiple HD 167971

Authors: J. Sanchez-Bermudez, A. Alberdi, R. Schödel, W. Brandner, R. Galván-Madrid, J. C. Guirado, R. Herrero-Illana, C. A. Hummel, J. M. Marcaide, M. A. Pérez-Torres

Abstract: 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 b… ▽ More 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. △ Less

Submitted 29 January, 2019; originally announced January 2019.

Comments: Accepted for publication in Astronomy and Astrophysics, 7 pages, 6 figures

Journal ref: A&A 624, A55 (2019)

Extreme fragmentation and complex kinematics at the center of the L1287 cloud

Authors: Carmen Juárez, Hauyu-Baobab Liu, Josep M. Girart, Aina Palau, Gemma Busquet, Roberto Galván-Madrid, Naomi Hirano, Yuxin Lin

Abstract: The filamentary infrared dark cloud L1287 is actively forming a dense cluster of low-mass YSOs at its inner $\sim$0.1 pc region. To help understand the origin of this low-mass YSO cluster, the present work aims at resolving the gas structures and kinematics. We have performed $\sim$1$"$ angular resolution ($\sim$930 AU) SMA observations at $\sim$1.3 mm. From a $\sim$2$"$ resolution 1.3 mm continuu… ▽ More The filamentary infrared dark cloud L1287 is actively forming a dense cluster of low-mass YSOs at its inner $\sim$0.1 pc region. To help understand the origin of this low-mass YSO cluster, the present work aims at resolving the gas structures and kinematics. We have performed $\sim$1$"$ angular resolution ($\sim$930 AU) SMA observations at $\sim$1.3 mm. From a $\sim$2$"$ resolution 1.3 mm continuum image we identified six dense cores, namely SMA1-6 with masses in the range of $\sim0.4-4$ M$_\odot$. From a $\sim$1$"$ resolution 1.3 mm continuum image, we find a high fragmentation level, with 14 compact millimeter sources within 0.1 pc (two of them associated with the known accretion outburst YSOs RNO 1C and RNO 1B). The dense gas tracer DCN (3--2) traces well the dust continuum emission and shows a clear velocity gradient along the NW-SE direction centered at SMA3. There is another velocity gradient with opposite direction around the most luminous YSO IRAS 00338+6312. The fragmentation within 0.1 pc in L1287 is very high compared to other regions at the same spatial scales. The incoherent motions of dense gas flows are sometimes interpreted by being influenced by (proto)stellar feedback (e.g., outflows), which is not yet ruled out in this particular target source. The directions of the velocity gradients traced by DCN are approximately perpendicular to those of the dominant CO outflow(s). Therefore, we alternatively hypothesize that the velocity gradients revealed by DCN trace the convergence from the $\gtrsim$0.1 pc scales infalling motion towards the rotational motions around the more compact ($\sim0.02$ pc) sources. This global molecular gas converging flow may feed the formation of the dense low-mass YSO cluster. IRAS 00338+6312 is the most likely powering source of the dominant CO outflow. A compact blue-shifted outflow from RNO 1C is also identified. △ Less

Submitted 20 November, 2018; originally announced November 2018.

Comments: 14 pages, 11 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 621, A140 (2019)

Chasing discs around O-type (proto)stars - ALMA evidence for an SiO disc and disc wind from G17.64+0.16

Authors: L. T. Maud, R. Cesaroni, M. S. N. Kumar, F. F. S. van der Tak, V. Allen, M. G. Hoare, P. D. Klaassen, D. Harsono, M. R. Hogerheijde, Á. Sánchez-Monge, P. Schilke, A. Ahmadi, M. T. Beltrán, H. Beuther, T. Csengeri, S. Etoka, G. Fuller, R. Galván-Madrid, C. Goddi, Th. Henning, K. G. Johnston, R. Kuiper, S. Lumsden, L. Moscadelli, J. C. Mottram , et al. (6 additional authors not shown)

Abstract: 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 peak… ▽ More 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. △ Less

Submitted 9 October, 2018; originally announced October 2018.

Comments: 23 pages, 13 fig (inc 3 Appendix figures)

Journal ref: A&A 620, A31 (2018)

On the Effects of Self-Obscuration in the (Sub-)Millimeter Spectral Indices and Appearance of Protostellar Disks

Authors: Roberto Galván-Madrid, Hauyu Baobab Liu, Andrés F. Izquierdo, Anna Miotello, Bo Zhao, Carlos Carrasco-González, Susana Lizano, Luis F. Rodríguez

Abstract: In this paper we explore the effects of self-obscuration in protostellar disks with a radially decreasing temperature gradient and a colder midplane. We are motivated by recent reports of resolved dark lanes (`hamburgers') and (sub)mm spectral indices systematically below the ISM value for optically thin dust $α_{\rm ISM} =3.7$. We explore several model grids, scaling disk mass and varying inclina… ▽ More In this paper we explore the effects of self-obscuration in protostellar disks with a radially decreasing temperature gradient and a colder midplane. We are motivated by recent reports of resolved dark lanes (`hamburgers') and (sub)mm spectral indices systematically below the ISM value for optically thin dust $α_{\rm ISM} =3.7$. We explore several model grids, scaling disk mass and varying inclination angle $i$ and observing frequency $ν$ from the VLA Ka band ($\sim 37$ GHz) to ALMA Band 8 ($\sim 405$ GHz). We also consider the effects of decreasing the index of the (sub-)mm dust opacity power law $β$ from 1.7 to 1. We find that a distribution of disk masses in the range $M_{\rm disk} = 0.01-2~M_\odot$ is needed to reproduce the observed distribution of spectral indices, and that assuming a fixed $β=1.7$ gives better results than $β=1$. A wide distribution of disk masses is also needed to produce some cases with $α<2$, as reported for some sources in the literature. Such extremely low spectral indices arise naturally when the selected observing frequencies sample the appropriate change in the temperature structure of the optically thick model disk. Our results show that protostellar disk masses could often be underestimated by $> \times10$, and are consistent with recent hydrodynamical simulations. Although we do not rule out the possibility of some grain growth occurring within the short protostellar timescales, we conclude that self-obscuration needs to be taken into account. △ Less

Submitted 1 October, 2018; originally announced October 2018.

Comments: Accepted to The Astrophysical Journal Supplement

Investigating fragmentation of gas structures in OB cluster-forming molecular clump G33.92+0.11 with 1000 AU resolution observations of ALMA

Authors: Hauyu Baobab Liu, Huei-Ru Vivien Chen, Carlos G. Román-Zúñiga, Roberto Galván-Madrid, Adam Ginsburg, Paul T. P. Ho, Young Chol Minh, Izaskun Jiménez-Serra, Leonardo Testi, Qizhou Zhang

Abstract: We report new, $\sim$1000 AU spatial resolution observations of 225 GHz dust continuum emission towards the OB cluster-forming molecular clump G33.92+0.11. On parsec scales, this molecular clump presents a morphology with several arm-like dense gas structures surrounding the two central massive ($\gtrsim$100 $M_{\odot}$) cores. From the new, higher resolution observations, we identified 28 localiz… ▽ More We report new, $\sim$1000 AU spatial resolution observations of 225 GHz dust continuum emission towards the OB cluster-forming molecular clump G33.92+0.11. On parsec scales, this molecular clump presents a morphology with several arm-like dense gas structures surrounding the two central massive ($\gtrsim$100 $M_{\odot}$) cores. From the new, higher resolution observations, we identified 28 localized, spatially compact dust continuum emission sources, which may be candidates of young stellar objects. Only one of them is not embedded within known arm-like (or elongated) dense gas structures. The spatial separations of these compact sources can be very well explained by Jeans lengths. We found that G33.92+0.11 may be consistently described by a marginally centrifugally supported, Toomre unstable accretion flow which is approximately in a face-on projection. The arm-like overdensities are natural consequence of the Toomre instability, which can fragment to form young stellar objects in shorter time scales than the timescale of the global clump contraction. On our resolved spatial scales, there is not yet evidence that the fragmentation is halted by turbulence, magnetic field, or stellar feedback. △ Less

Submitted 4 December, 2018; v1 submitted 23 August, 2018; originally announced August 2018.

Comments: 24 pages, 18 figures. Accepted to publish on December 04, 2018; updated to arXiv on December 05, 2018