Showing 1–50 of 116 results for author: Leurini, S

Massive clumps in W43-main: Structure formation in an extensively shocked molecular cloud

Authors: Yuxin Lin, Friedrich Wyrowski, Hauyu Baobab Liu, Yan Gong, Olli Sipilä, Andrés F. Izquierdo, Timea Csengeri, Adam Ginsburg, Guang-Xing Li, Silvia Spezzano, Jaime E. Pineda, Silvia Leurini, Paola Caselli, Karl M. Menten

Abstract: W43-main is a massive molecular complex located at the interaction of the Scutum arm and the Galactic bar undergoing starburst activities. We aim to investigate the gas dynamics, in particular, the prevailing shock signatures from the cloud to clump scale and assess the impact of shocks on the formation of dense gas and early-stage cores. We have carried out NOEMA and IRAM-30m observations at 3 mm… ▽ More W43-main is a massive molecular complex located at the interaction of the Scutum arm and the Galactic bar undergoing starburst activities. We aim to investigate the gas dynamics, in particular, the prevailing shock signatures from the cloud to clump scale and assess the impact of shocks on the formation of dense gas and early-stage cores. We have carried out NOEMA and IRAM-30m observations at 3 mm with an angular resolution of $\sim$0.1 pc towards five massive clumps in W43 main. We use CH$_{3}$CCH and H$_{2}$CS lines to trace the extended gas temperature and CH$_{3}$OH lines to probe the volume density of the dense gas ($\gtrsim$10$^{5}$ cm$^{-3}$). The emission of SiO (2-1) is extensive across the region ($\sim$4 pc) and is mostly contained within a low-velocity regime, hinting at a large-scale origin of the shocks. The position-velocity maps of multiple tracers show systematic spatio-kinematic offsets supporting the cloud-cloud collision/merging scenario. We identify an additional extended velocity component in CCH emission, which coincides with one of the velocity components of the larger scale $^{13}$CO (2-1) emission, likely representing an outer, less dense gas layer in the cloud merging process. We find that the V-shaped, asymmetric SiO wings are tightly correlated with localised gas density enhancements, which is direct evidence of dense gas formation and accumulation in shocks. We resolve two categories of NH$_{2}$D cores: ones exhibiting only subsonic to transonic velocity dispersion, and the others with an additional supersonic velocity dispersion. The centroid velocities of the latter cores are correlated with the shock front seen by SiO. The kinematics of the $\sim$0.1 pc NH$_{2}$D cores are heavily imprinted by shock activities, and may represent a population of early-stage cores forming around the shock interface. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Comments: 35 pages, 34 figures (including appendices); Accepted in A&A

OGHReS: Star formation in the Outer Galaxy ($\ell = 250^\circ$-$280^\circ$)

Authors: J. S. Urquhart, C. König, D. Colombo, A. Karska, F. Wyrowski, K. M. Menten, T. J. T. Moore, J. Brand, D. Elia, A. Giannetti, S. Leurini, M. Figueira, M. -Y. Lee, M. Dumke

Abstract: We have used data from the Outer Galaxy High-Resolution Survey (OGHReS) to refine the velocities, distances, and physical properties of a large sample of 3584 clumps detected in far infrared/submillimetre emission in the HiGAL survey located in the $\ell = 250^\circ-280^\circ$ region of the Galactic plane. Using $^{12}$CO and $^{13}$CO spectra, we have determined reliable velocities to 3412 clumps… ▽ More We have used data from the Outer Galaxy High-Resolution Survey (OGHReS) to refine the velocities, distances, and physical properties of a large sample of 3584 clumps detected in far infrared/submillimetre emission in the HiGAL survey located in the $\ell = 250^\circ-280^\circ$ region of the Galactic plane. Using $^{12}$CO and $^{13}$CO spectra, we have determined reliable velocities to 3412 clumps (95% of the sample). In comparison to the velocities from the HiGAL catalogue, we find good agreement for 80% of the sample (within 5 km/s). Using the higher resolution and sensitivity of OGHReS has allowed us to correct the velocity for 632 clumps and provide velocities for 687 clumps for which no velocity had been previously allocated. The velocities are used with a rotation curve to refine the distances to the clumps and to calculate the clumps' properties using a distance-dependent gas-to-dust ratio. We have determined reliable physical parameters for 3200 outer Galaxy dense clumps (~90% of the HiGAL sources in the region). We find a trend of decreasing luminosity-to-mass ratio with increasing Galactocentric distance, suggesting the star formation efficiency is lower in the outer Galaxy or that it is resulting in more lower mass stars than in the inner Galaxy. We also find a similar surface density for protostellar clumps located in the inner and outer Galaxy, revealing that the surface density requirements for star formation are the same across the Galactic disc. △ Less

Submitted 1 January, 2024; originally announced January 2024.

Comments: 15 pages, 13 figures. Accepted for publication by MNRAS

Density distributions, magnetic field structures and fragmentation in high-mass star formation

Authors: H. Beuther, C. Gieser, J. D. Soler, Q. Zhang, R. Rao, D. Semenov, Th. Henning, R. Pudritz, T. Peters, P. Klaassen, M. T. Beltran, A. Palau, T. Moeller, K. G. Johnston, H. Zinnecker, J. Urquhart, R. Kuiper, A. Ahmadi, A. Sanchez-Monge, S. Feng, S. Leurini, S. E. Ragan

Abstract: Methods: Observing the large pc-scale Stokes I mm dust continuum emission with the IRAM 30m telescope and the intermediate-scale (<0.1pc) polarized submm dust emission with the Submillimeter Array toward a sample of 20 high-mass star-forming regions allows us to quantify the dependence of the fragmentation behaviour of these regions depending on the density and magnetic field structures. Results… ▽ More Methods: Observing the large pc-scale Stokes I mm dust continuum emission with the IRAM 30m telescope and the intermediate-scale (<0.1pc) polarized submm dust emission with the Submillimeter Array toward a sample of 20 high-mass star-forming regions allows us to quantify the dependence of the fragmentation behaviour of these regions depending on the density and magnetic field structures. Results: We infer density distributions n~r^{-p} of the regions with typical power-law slopes p around ~1.5. There is no obvious correlation between the power-law slopes of the density structures on larger clump scales (~1pc) and the number of fragments on smaller core scales (<0.1pc). Comparing the large-scale single-dish density profiles to those derived earlier from interferometric observations at smaller spatial scales, we find that the smaller-scale power-law slopes are steeper, typically around ~2.0. The flattening toward larger scales is consistent with the star-forming regions being embedded in larger cloud structures that do not decrease in density away from a particular core. Regarding the magnetic field, for several regions it appears aligned with filamentary structures leading toward the densest central cores. Furthermore, we find different polarization structures with some regions exhibiting central polarization holes whereas other regions show polarized emission also toward the central peak positions. Nevertheless, the polarized intensities are inversely related to the Stokes I intensities. We estimate magnetic field strengths between ~0.2 and ~4.5mG, and we find no clear correlation between magnetic field strength and the fragmentation level of the regions. Comparison of the turbulent to magnetic energies shows that they are of roughly equal importance in this sample. The mass-to-flux ratios range between ~2 and ~7, consistent with collapsing star-forming regions. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Comments: Accepted for Astronomy & Astrophysics, 14 pages, 14 figures plus appendices, also download option at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

One, Two, Three ... An Explosive Outflow in IRAS 12326$-$6245 revealed by ALMA

Authors: Luis A. Zapata, Manuel Fernández-López, Silvia Leurini, Estrella Guzmán Ccolque, Luis F. Rodriguez, Aina Palau, Karl M. Menten, Friedrich Wyrowski

Abstract: In the last years there has been a substantial increase in the number of the reported massive and luminous star-forming regions with related explosive outflows thanks to the superb sensitivity and angular resolution provided by the new radio, infrared, and optical facilities. Here, we report one more explosive outflow related with the massive and bright star-forming region IRAS 12326$-$6245 using… ▽ More In the last years there has been a substantial increase in the number of the reported massive and luminous star-forming regions with related explosive outflows thanks to the superb sensitivity and angular resolution provided by the new radio, infrared, and optical facilities. Here, we report one more explosive outflow related with the massive and bright star-forming region IRAS 12326$-$6245 using Band 6 sensitive and high angular resolution ($\sim$0.2$"$) Atacama Large Millimeter/Submillimeter Array (ALMA) observations. We find over 10 molecular and collimated well-defined streamers, with Hubble-Lemaitre like expansion motions, and pointing right to the center of a dusty and molecular shell (reported for the first time here) localized in the northern part of the UCHII region known as G301.1A. The estimated kinematic age, and energy for the explosion are $\sim$700 yrs, and 10$^{48}$ erg, respectively. Taking into account the recently reported explosive outflows together with IRAS 12326$-$6245, we estimate an event rate of once every 90 yr in our Galaxy, similar to the formation rate of massive stars. △ Less

Submitted 20 September, 2023; originally announced September 2023.

Comments: Accepted to The Astrophysical Journal Letters

Infall and Outflow Towards High-mass Starless Clump Candidates

Authors: T. G. S. Pillai, J. S. Urquhart, S. Leurini, Q. Zhang, A. Traficante, D. Colombo, K. Wang, L. Gomez, F. Wyrowski

Abstract: The evolutionary sequence for high-mass star formation starts with massive starless clumps that go on to form protostellar, young stellar objects and then compact HII regions. While there are many examples of the three later stages, the very early stages have proved to be elusive. We follow-up a sample of 110 mid-infrared dark clumps selected from the ATLASGAL catalogue with the IRAM telescope in… ▽ More The evolutionary sequence for high-mass star formation starts with massive starless clumps that go on to form protostellar, young stellar objects and then compact HII regions. While there are many examples of the three later stages, the very early stages have proved to be elusive. We follow-up a sample of 110 mid-infrared dark clumps selected from the ATLASGAL catalogue with the IRAM telescope in an effort to identify a robust sample of massive starless clumps. We have used the HCO+ (1-0) and HNC (1-0) transitions to identify clumps associated with infall motion and the SiO (2-1) transition to identity outflow candidates. We have found blue asymmetric line profile in 65% of the sample, and have measured the infall velocities and mass infall rates (0.6-$36 \times 10^{-3}$ Msun/yr) for 33 of these clumps. We find a trend for the mass infall rate decreasing with an increase of bolometric luminosity to clump mass i.e. star formation within the clumps evolves. Using the SiO 2-1 line, we have identified good outflow candidates. Combining the infall and outflow tracers reveals that 67% of quiescent clumps are already undergoing gravitational collapse or are associated with star formation; these clumps provide us with our best opportunity to determined the initial conditions and study the earliest stages of massive star formation. Finally, we provide an overview of a systematic high-resolution ALMA study of quiescent clumps selected that allows us to develop a detailed understanding of earliest stages and their subsequent evolution. △ Less

Submitted 7 May, 2023; originally announced May 2023.

Comments: 10 pages and 7 figures

Journal ref: 2023MNRAS.522.3357P

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)

The SOUL view of IRAS20126+4104. Kinematics and variability of the H$_2$ jet from a massive protostar

Authors: F. Massi, A. Caratti o Garatti, R. Cesaroni, T. K. Sridharan, E. Ghose, E. Pinna, M. T. Beltrán, S. Leurini, L. Moscadelli, A. Sanna, G. Agapito, R. Briguglio, J. Christou, S. Esposito, T. Mazzoni, D. Miller, C. Plantet, J. Power, A. Puglisi, F. Rossi, B. Rothberg, G. Taylor, C. Veillet

Abstract: We exploit the increased sensitivity of the recently installed AO SOUL at the LBT to obtain new high-spatial-resolution NIR images of the massive young stellar object IRAS20126+4104 and its outflow. We aim to derive the jet proper motions and kinematics, as well as to study its photometric variability by combining the novel performances of SOUL together with previous NIR images. We used both broad… ▽ More We exploit the increased sensitivity of the recently installed AO SOUL at the LBT to obtain new high-spatial-resolution NIR images of the massive young stellar object IRAS20126+4104 and its outflow. We aim to derive the jet proper motions and kinematics, as well as to study its photometric variability by combining the novel performances of SOUL together with previous NIR images. We used both broad-band ($K_{s}$, $K'$) and narrow-band (Br$γ$, H2) observations from a number of NIR cameras (UKIRT/UFTI,SUBARU/CIAO,TNG/NICS,LBT/PISCES,and LBT/LUCI1) to derive maps of the continuum and the H$_2$ emission in the 2.12 $μ$m line. Three sets of images, obtained with AO systems (CIAO,2003; FLAO,2012; SOUL,2020), allowed us to derive the proper motions of a large number of H$_2$ knots along the jet. Photometry from all images was used to study the jet variability. We derived knot proper motions in the range of 1.7-20.3 mas yr$^{-1}$ (i.e. 13-158 km s$^{-1}$ at 1.64 kpc, avg. outflow tangential velocity $\sim$ 80 km s$^{-1}$). The derived knot dynamical age spans a $\sim$ 200-4000 yr interval. A ring-like H$_2$ feature near the protostar location exhibits peculiar kinematics and may represent the outcome of a wide-angle wind impinging on the outflow cavity. Both H$_2$ geometry and velocities agree with those inferred from proper motions of the H$_2$O masers, located at a smaller distance from the protostar. Although the total H$_2$ line emission from the knots does not exhibit time variations at a $\widetilde{>}$ 0.3 mag level, we have found a clear continuum flux variation (radiation scattered by the dust in the cavity opened by the jet) which is anti-correlated between the blue-shifted and red-shifted lobes and may be periodic (with a period of $\sim$ 12-18 yr). We suggest that the continuum variability might be related to inner-disc oscillations which have also caused the jet precession. △ Less

Submitted 17 January, 2023; originally announced January 2023.

Comments: 26 pages, 22 figures, 2 mpeg files, accepted by Astronomy & Astrophysics

Journal ref: A&A 672, A113 (2023)

The Star Formation Rate of the Milky Way as seen by Herschel

Authors: D. Elia, S. Molinari, E. Schisano, J. D. Soler, M. Merello, D. Russeil, M. Veneziani, A. Zavagno, A. Noriega-Crespo, L. Olmi, M. Benedettini, P. Hennebelle, R. S. Klessen, S. Leurini, R. Paladini, S. Pezzuto, A. Traficante, D. J. Eden, P. G. Martin, M. Sormani, A. Coletta, T. Colman, R. Plume, Y. Maruccia, C. Mininni , et al. (1 additional authors not shown)

Abstract: We present a new derivation of the Milky Way's current star formation rate (SFR) based on the data of the Hi-GAL Galactic plane survey. We estimate the distribution of the SFR across the Galactic plane from the star-forming clumps identified in the Hi-GAL survey and calculate the total SFR from the sum of their contributions. The estimate of the global SFR amounts to $2.0 \pm 0.7$~M$_{\odot}$~yr… ▽ More We present a new derivation of the Milky Way's current star formation rate (SFR) based on the data of the Hi-GAL Galactic plane survey. We estimate the distribution of the SFR across the Galactic plane from the star-forming clumps identified in the Hi-GAL survey and calculate the total SFR from the sum of their contributions. The estimate of the global SFR amounts to $2.0 \pm 0.7$~M$_{\odot}$~yr$^{-1}$, of which $1.7 \pm 0.6$~M$_{\odot}$~yr$^{-1}$ coming from clumps with reliable heliocentric distance assignment. This value is in general agreement with estimates found in the literature of last decades. The profile of SFR density averaged in Galactocentric rings is found to be qualitatively similar to others previously computed, with a peak corresponding to the Central Molecular Zone and another one around Galactocentric radius $R_\mathrm{gal} \sim 5$~kpc, followed by an exponential decrease as $\log(Σ_\mathrm{SFR}/[\mathrm{M}_\odot~\mathrm{yr}^{-1}~\mathrm{kpc}^{-2}])=a\,R_\mathrm{gal}/[\mathrm{kpc}]+b $, with $a=-0.28 \pm 0.01$. In this regard, the fraction of SFR produced within and outside the Solar circle is 84\% and 16\%, respectively; the fraction corresponding to the far outer Galaxy ($R_\mathrm{gal} > 13.5$~kpc) is only 1\%. We also find that, for $R_\mathrm{gal}>3$~kpc, our data follow a power law as a function of density, similarly to the Kennicutt-Schmidt relation. Finally, we compare the distribution of the SFR density across the face-on Galactic plane and those of median parameters, such as temperature, luminosity/mass ratio and bolometric temperature, describing the evolutionary stage of Hi-GAL clumps. We found no clear correlation between the SFR and the clump evolutionary stage. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: To be published in ApJ

How magnetic field and stellar radiative feedback influences the collapse and the stellar mass spectrum of a massive star forming clump

Authors: Patrick Hennebelle, Ugo Lebreuilly, Tine Colman, Davide Elia, Gary Fuller, Silvia Leurini, Thomas Nony, Eugenio Schisano, Juan D. Soler, Alessio Traficante, Ralf S. Klessen, Sergio Molinari, Leonardo Testi

Abstract: In spite of decades of theoretical efforts, the physical origin of the stellar initial mass function (IMF) is still debated. We aim at understanding the influence of various physical processes such as radiative stellar feedback, magnetic field and non-ideal magneto-hydrodynamics on the IMF. We present a series of numerical simulations of collapsing 1000 M$_\odot$ clumps taking into account radiati… ▽ More In spite of decades of theoretical efforts, the physical origin of the stellar initial mass function (IMF) is still debated. We aim at understanding the influence of various physical processes such as radiative stellar feedback, magnetic field and non-ideal magneto-hydrodynamics on the IMF. We present a series of numerical simulations of collapsing 1000 M$_\odot$ clumps taking into account radiative feedback and magnetic field with spatial resolution down to 1 AU. Both ideal and non-ideal MHD runs are performed and various radiative feedback efficiencies are considered. We also develop analytical models that we confront to the numerical results. The sum of the luminosities produced by the stars in the calculations is computed and it compares well with the bolometric luminosities reported in observations of massive star forming clumps. The temperatures, velocities and densities are also found to be in good agreement with recent observations. The stellar mass spectrum inferred for the simulations is, generally speaking, not strictly universal and in particular varies with magnetic intensity. It is also influenced by the choice of the radiative feedback efficiency. In all simulations, a sharp drop in the stellar distribution is found at about $M_{min} \simeq$ 0.1 M$_\odot$, which is likely a consequence of the adiabatic behaviour induced by dust opacities at high densities. As a consequence, when the combination of magnetic and thermal support is not too large, the mass distribution presents a peak located at 0.3-0.5 M$_\odot$. When magnetic and thermal support are large, the mass distribution is better described by a plateau, i.e. $d N / d \log M \propto M^{-Γ}$, $Γ\simeq 0$. Abridged △ Less

Submitted 22 October, 2022; originally announced October 2022.

Comments: accepted for publication in A&A

Journal ref: A&A 668, A147 (2022)

The SEDIGISM survey: Molecular cloud morphology. II. Integrated source properties

Authors: K. R. Neralwar, D. Colombo, A. Duarte-Cabral, J. S. Urquhart, M. Mattern, F. Wyrowski, K. M. Menten, P. Barnes, A. Sanchez-Monge, A. J. Rigby, P. Mazumdar, D. Eden, T. Csengeri, C. L. Dobbs, V. S. Veena, S. Neupane, T. Henning, F. Schuller, S. Leurini, M. Wienen, A. Y. Yang, S. E. Ragan, S. Medina, Q. Nguyen-Luong

Abstract: The Structure, Excitation, and Dynamics of the Inner Galactic InterStellar Medium (SEDIGISM) survey has produced high (spatial and spectral) resolution $^{13}$CO (2-1) maps of the Milky Way. It has allowed us to investigate the molecular interstellar medium in the inner Galaxy at an unprecedented level of detail and characterise it into molecular clouds. In a previous paper, we have classified the… ▽ More The Structure, Excitation, and Dynamics of the Inner Galactic InterStellar Medium (SEDIGISM) survey has produced high (spatial and spectral) resolution $^{13}$CO (2-1) maps of the Milky Way. It has allowed us to investigate the molecular interstellar medium in the inner Galaxy at an unprecedented level of detail and characterise it into molecular clouds. In a previous paper, we have classified the SEDIGISM clouds into four morphologies. However, how the properties of the clouds vary for these four morphologies is not well understood. Here, we use the morphological classification of SEDIGISM clouds to find connections between the cloud morphologies, their integrated properties, and their location on scaling relation diagrams. We observe that ring-like clouds show the most peculiar properties, having, on average, higher masses, sizes, aspect ratios and velocity dispersions compared to other morphologies. We speculate that this is related to the physical mechanisms that regulate their formation and evolution, for example, turbulence from stellar feedback can often results in the creation of bubble-like structures. We also see a trend of morphology with virial parameter whereby ring-like, elongated, clumpy and concentrated clouds have virial parameters in a decreasing order. Our findings provide a foundation for a better understanding of the molecular cloud behaviour based on their measurable properties. △ Less

Submitted 4 May, 2022; originally announced May 2022.

Comments: Accepted for publication in A&A. 28 pages (17 of Appendices), 32 figures, 6 tables

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

The SEDIGISM survey: Molecular cloud morphology. I. Classification and star formation

Authors: K. R. Neralwar, D. Colombo, A. Duarte-Cabral, J. S. Urquhart, M. Mattern, F. Wyrowski, K. M. Menten, P. Barnes, A. Sanchez-Monge, H. Beuther, A. J. Rigby, P. Mazumdar, D. Eden, T. Csengeri, C. L. Dobbs, V. S. Veena, S. Neupane, T. Henning, F. Schuller, S. Leurini, M. Wienen, A. Y. Yang, S. E. Ragan, S. Medina, Q. Nguyen-Luong

Abstract: We present one of the very first extensive classifications of a large sample of molecular clouds based on their morphology. This is achieved using a recently published catalogue of 10663 clouds obtained from the first data release of the SEDIGISM survey. The clouds are classified into four different morphologies by visual inspection and using an automated algorithm -- J plots. The visual inspectio… ▽ More We present one of the very first extensive classifications of a large sample of molecular clouds based on their morphology. This is achieved using a recently published catalogue of 10663 clouds obtained from the first data release of the SEDIGISM survey. The clouds are classified into four different morphologies by visual inspection and using an automated algorithm -- J plots. The visual inspection also serves as a test for the J plots algorithm, as this is the first time it has been used on molecular gas. Generally, it has been found that the structure of molecular clouds is highly filamentary and our observations indeed verify that most of our molecular clouds are elongated structures. Based on our visual classification of the 10663 SEDIGISM clouds, 15% are ring-like, 57% are elongated, 15% are concentrated and 10% are clumpy clouds. The remaining clouds do not belong to any of these morphology classes and are termed unclassified. We compare the SEDIGISM molecular clouds with structures identified through other surveys, i.e. ATLASGAL elongated structures and the bubbles from Milky Way Project (MWP). We find that many of the ATLASGAL and MWP structures are velocity coherent. ATLASGAL elongated structures overlap with ~ 21% of the SEDIGISM elongated structures (elongated and clumpy clouds) and MWP bubbles overlap with ~ 25% of the SEDIGISM ring-like clouds. We also analyse the star-formation associated with different cloud morphologies using two different techniques. The first technique examines star formation efficiency (SFE) and the dense gas fraction (DGF), based on SEDIGISM clouds and ATLASGAL clumps data. The second technique uses the high-mass star formation (HMSF) threshold for molecular clouds. The results indicate that clouds with ring-like and clumpy morphologies show a higher degree of star formation. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: Accepted for publication in A&A. 40 pages (26 of Appendices), 55 figures, 13 tables. The updated SEDIGISM cloud catalogue, containing cloud morphology, will be available as part of the SEDIGISM database

Journal ref: A&A 663, A56 (2022)

NIR jets from a clustered region of massive star formation: Morphology and composition in the IRAS 18264-1152 region

Authors: A. R. Costa Silva, R. Fedriani, J. C. Tan, A. Caratti o Garatti, S. Ramsay, V. Rosero, G. Cosentino, P. Gorai, S. Leurini

Abstract: Massive stars form deeply embedded in their parental clouds, making it challenging to directly observe these stars and their immediate environments. It is known that accretion and ejection processes are intrinsically related, thus observing massive protostellar outflows can provide crucial information about the processes governing massive star formation close to the central engine. We aim to probe… ▽ More Massive stars form deeply embedded in their parental clouds, making it challenging to directly observe these stars and their immediate environments. It is known that accretion and ejection processes are intrinsically related, thus observing massive protostellar outflows can provide crucial information about the processes governing massive star formation close to the central engine. We aim to probe the IRAS 18264-1152 (G19.88-0.53) high-mass star-forming complex in the near infrared (NIR) through its molecular hydrogen (H2) jets to analyse the morphology and composition of the line emitting regions and to compare with other outflow tracers. We observed the H2 NIR jets via K-band (1.9-2.5um) observations obtained with the integral field units VLT/SINFONI and VLT/KMOS. SINFONI provides the highest NIR angular resolution achieved so far for the central region (~0.2''). We compared the geometry of the NIR outflows with that of the associated molecular outflow probed by CO (2-1) emission mapped with SMA. We identify nine point sources. Four of these display a rising continuum in the K-band and are BrG emitters, revealing that they are young, potentially jet-driving sources. The spectro-imaging analysis focusses on the H2 jets, for which we derived visual extinction, temperature, column density, area, and mass. The intensity, velocity, and excitation maps based on H2 emission strongly support the existence of a protostellar cluster, with at least two (and up to four) different large-scale outflows. The literature is in agreement with the outflow morphology found here. We derived a stellar density of ~4000 stars pc^-3. Our study reveals the presence of several outflows driven by young sources from a forming cluster of young, massive stars. The derived stellar number density together with the geometry of the outflows suggest that stars can form in a relatively ordered manner in this cluster. △ Less

Submitted 19 July, 2022; v1 submitted 8 December, 2021; originally announced December 2021.

Comments: 23 pages, 14 figures. Published in A&A. Updated to correct typos and match proofs

Journal ref: A&A 659, A23 (2022)

The evolution of temperature and density structures of OB cluster-forming molecular clumps

Authors: Yuxin Lin, Friedrich Wyrowski, Hauyu Baobab Liu, Andrés Izquierdo, Timea Csengeri, Silvia Leurini, Karl M. Menten

Abstract: OB star clusters originate from parsec-scale massive molecular clumps. We aim to understand the evolution of temperature and density structures on the intermediate-scale ($\lesssim$0.1-1 pc) extended gas of massive clumps. We performed $\sim$0.1 pc resolution observations (SMA+APEX) of multiple molecular line tracers (e.g., CH$_{3}$CCH, H$_{2}$CS, CH$_{3}$CN, CH$_{3}$OH) which cover a wide range o… ▽ More OB star clusters originate from parsec-scale massive molecular clumps. We aim to understand the evolution of temperature and density structures on the intermediate-scale ($\lesssim$0.1-1 pc) extended gas of massive clumps. We performed $\sim$0.1 pc resolution observations (SMA+APEX) of multiple molecular line tracers (e.g., CH$_{3}$CCH, H$_{2}$CS, CH$_{3}$CN, CH$_{3}$OH) which cover a wide range of excitation conditions, towards a sample of eight massive clumps. Based on various radiative transfer models, we constrain the gas temperature and density structures and establish an evolutionary picture, aided by a spatially-dependent virial analysis and abundance ratios of multiple species. We determine temperature radial profiles varying between 30-200 K over a continuous scale, from the center of the clumps out to 0.3-0.4 pc radii. The clumps' radial gas density profiles, described by radial power-laws with slopes between -0.6 and $\sim$-1.5, are steeper for more evolved sources, as suggested by results based on both dust continuum, representing the bulk of the gas ($\sim$10$^{4}$ cm$^{-3}$), and CH$_{3}$OH lines probing the dense gas ($\gtrsim$10$^{6}$-10$^{8}$ cm$^{-3}$) regime. The density contrast between the dense gas and the bulk gas increases with evolution, and may be indicative of spatially and temporally varying star formation efficiencies. The radial profiles of the virial parameter show a global variation towards a sub-virial state as the clump evolves. The line-widths decline with increasing radius around the central core region and increase in the outer envelope, with a slope shallower than the case of the supersonic turbulence ($\,\propto\,$$r^{0.5}$) and the subsonic Kolmogorov scaling ($\,\propto\,$$r^{0.33}$). In the context of clump evolution, we also find that the abundance ratios of [CCH]/[CH$_{3}$OH] and [CH$_{3}$CN]/[CH$_{3}$OH] show correlations with clump $L/M$. △ Less

Submitted 2 December, 2021; originally announced December 2021.

Comments: Accepted by A&A; 46 pages, 33 figures incl. appendices; abstract abridged to meet arxiv requirements

Journal ref: A&A 658, A128 (2022)

ATLASGAL -- Evolutionary trends in high-mass star formation

Authors: J. S. Urquhart, M. R. A. Wells, T. Pillai, S. Leurini, A. Giannetti, T. J. T. Moore, M. A. Thompson, C. Figura, D. Colombo, A. Y. Yang, C. Koenig, F. Wyrowski, K. M. Menten, A. J. Rigby, D. J. Eden, S. E. Ragan

Abstract: ATLASGAL is a 870-mircon dust survey of 420 square degrees of the inner Galactic plane and has been used to identify ~10 000 dense molecular clumps. Dedicated follow-up observations and complementary surveys are used to characterise the physical properties of these clumps, map their Galactic distribution and investigate the evolutionary sequence for high-mass star formation. The analysis of the AT… ▽ More ATLASGAL is a 870-mircon dust survey of 420 square degrees of the inner Galactic plane and has been used to identify ~10 000 dense molecular clumps. Dedicated follow-up observations and complementary surveys are used to characterise the physical properties of these clumps, map their Galactic distribution and investigate the evolutionary sequence for high-mass star formation. The analysis of the ATLASGAL data is ongoing: we present an up-to-date version of the catalogue. We have classified 5007 clumps into four evolutionary stages (quiescent, protostellar, young stellar objects and HII regions) and find similar numbers of clumps in each stage, suggesting a similar lifetime. The luminosity-to-mass (L/M) ratio curve shows a smooth distribution with no significant kinks or discontinuities when compared to the mean values for evolutionary stages indicating that the star-formation process is continuous and that the observational stages do not represent fundamentally different stages or changes in the physical mechanisms involved. We compare the evolutionary sample with other star-formation tracers (methanol and water masers, extended green objects and molecular outflows) and find that the association rates with these increases as a function of evolutionary stage, confirming that our classification is reliable. This also reveals a high association rate between quiescent sources and molecular outflows, revealing that outflows are the earliest indication that star formation has begun and that star formation is already ongoing in many of the clumps that are dark even at 70 micron. △ Less

Submitted 1 December, 2021; v1 submitted 24 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS. Consists of 20 pages, 15 figures, 8 table. The complete tables will be available from CDS and upon request

OGHReS: Large-scale filaments in the outer Galaxy

Authors: D. Colombo, C. König, J. S. Urquhart, F. Wyrowski, M. Mattern, K. M. Menten, M. -Y. Lee, J. Brand, M. Wienen, P. Mazumdar, F. Schuller, S. Leurini

Abstract: Filaments are a ubiquitous morphological feature of the molecular interstellar medium and are identified as sites of star formation. In recent years, more than 100 large-scale filaments (with a length $>10$\,pc) have been observed in the inner Milky Way. As they appear linked to Galactic dynamics, studying those structures represents an opportunity to link kiloparsec-scale phenomena to the physics… ▽ More Filaments are a ubiquitous morphological feature of the molecular interstellar medium and are identified as sites of star formation. In recent years, more than 100 large-scale filaments (with a length $>10$\,pc) have been observed in the inner Milky Way. As they appear linked to Galactic dynamics, studying those structures represents an opportunity to link kiloparsec-scale phenomena to the physics of star formation, which operates on much smaller scales. In this letter, we use newly acquired Outer Galaxy High Resolution Survey (OGHReS) $^{12}$CO(2-1) data to demonstrate that a significant number of large-scale filaments are present in the outer Galaxy as well. The 37 filaments identified appear tightly associated with inter-arm regions. In addition, their masses and linear masses are, on average, one order of magnitude lower than similar-sized molecular filaments located in the inner Galaxy, showing that Milky Way dynamics is able to create very elongated features in spite of the lower gas supply in the Galactic outskirts. △ Less

Submitted 4 November, 2021; originally announced November 2021.

Comments: 5 pages (+ appendices, 9 pages), 8 figures, 1 table, accepted for publication in A&A Letters

Journal ref: A&A 655, L2 (2021)

Astronomical source finding services for the CIRASA visual analytic platform

Authors: S. Riggi, C. Bordiu, F. Vitello, G. Tudisco, E. Sciacca, D. Magro, R. Sortino, C. Pino, M. Molinaro, M. Benedettini, S. Leurini, F. Bufano, M. Raciti, U. Becciani

Abstract: Innovative developments in data processing, archiving, analysis, and visualization are nowadays unavoidable to deal with the data deluge expected in next-generation facilities for radio astronomy, such as the Square Kilometre Array (SKA) and its precursors. In this context, the integration of source extraction and analysis algorithms into data visualization tools could significantly improve and sp… ▽ More Innovative developments in data processing, archiving, analysis, and visualization are nowadays unavoidable to deal with the data deluge expected in next-generation facilities for radio astronomy, such as the Square Kilometre Array (SKA) and its precursors. In this context, the integration of source extraction and analysis algorithms into data visualization tools could significantly improve and speed up the cataloguing process of large area surveys, boosting astronomer productivity and shortening publication time. To this aim, we are developing a visual analytic platform (CIRASA) for advanced source finding and classification, integrating state-of-the-art tools, such as the CAESAR source finder, the ViaLactea Visual Analytic (VLVA) and Knowledge Base (VLKB). In this work, we present the project objectives and the platform architecture, focusing on the implemented source finding services. △ Less

Submitted 18 October, 2021; v1 submitted 15 October, 2021; originally announced October 2021.

Comments: 16 pages, 6 figures

Disk fragmentation in high-mass star formation. High-resolution observations towards AFGL 2591-VLA 3

Authors: S. Suri, H. Beuther, C. Gieser, A. Ahmadi, Á. Sánchez-Monge, J. M. Winters, H. Linz, Th. Henning, M. T. Beltrán, F. Bosco, R. Cesaroni, T. Csengeri, S. Feng, M. G. Hoare, K. G. Johnston, P. Klaasen, R. Kuiper, S. Leurini, S. Longmore, S. Lumsden, L. Maud, L. Moscadelli, T. Möller, A. Palau, T. Peters , et al. (7 additional authors not shown)

Abstract: Increasing evidence suggests that, similar to their low-mass counterparts, high-mass stars form through a disk-mediated accretion process. At the same time, formation of high-mass stars still necessitates high accretion rates, and hence, high gas densities, which in turn can cause disks to become unstable against gravitational fragmentation. We study the kinematics and fragmentation of the disk ar… ▽ More Increasing evidence suggests that, similar to their low-mass counterparts, high-mass stars form through a disk-mediated accretion process. At the same time, formation of high-mass stars still necessitates high accretion rates, and hence, high gas densities, which in turn can cause disks to become unstable against gravitational fragmentation. We study the kinematics and fragmentation of the disk around the high-mass star forming region AFGL 2591-VLA 3 which was hypothesized to be fragmenting based on the observations that show multiple outflow directions. We use a new set of high-resolution (0.19 arcsec) IRAM/NOEMA observations at 843 micron towards VLA 3 which allow us to resolve its disk, characterize the fragmentation, and study its kinematics. In addition to the 843 micron continuum emission, our spectral setup targets warm dense gas and outflow tracers such as HCN, HC$_3$N and SO$_2$, as well as vibrationally excited HCN lines. The high resolution continuum and line emission maps reveal multiple fragments with subsolar masses within the inner 1000 AU of VLA 3. Furthermore, the velocity field of the inner disk observed at 843 micron shows a similar behavior to that of the larger scale velocity field studied in the CORE project at 1.37 mm. We present the first observational evidence for disk fragmentation towards AFGL 2591-VLA 3, a source that was thought to be a single high-mass core. While the fragments themselves are low-mass, the rotation of the disk is dominated by the protostar with a mass of 10.3$\pm 1.8~M_{\odot}$. These data also show that NOEMA Band 4 can obtain the highest currently achievable spatial resolution at (sub-)mm wavelengths in observations of strong northern sources. △ Less

Submitted 10 September, 2021; originally announced September 2021.

Comments: 16 pages, 11 figures, Accepted for publication in A&A

Journal ref: A&A 655, A84 (2021)

Establishing the evolutionary timescales of the massive star formation process through chemistry

Authors: G. Sabatini, S. Bovino, A. Giannetti, T. Grassi, J. Brand, E. Schisano, F. Wyrowski, S. Leurini, K. M. Menten

Abstract: (Abridged) Understanding the details of the formation process of massive (i.e. M<8-10M$_\odot$) stars is a long-standing problem in astrophysics. [...] We present a method to derive accurate timescales of the different evolutionary phases of the high-mass star formation process. We model a representative number of massive clumps of the ATLASGAL-TOP100 sample which cover all the evolutionary stages… ▽ More (Abridged) Understanding the details of the formation process of massive (i.e. M<8-10M$_\odot$) stars is a long-standing problem in astrophysics. [...] We present a method to derive accurate timescales of the different evolutionary phases of the high-mass star formation process. We model a representative number of massive clumps of the ATLASGAL-TOP100 sample which cover all the evolutionary stages. The models describe an isothermal collapse and the subsequent warm-up phase, for which we follow their chemical evolution. The timescale of each phase is derived by comparing the results of the models with the properties of the sources of the ATLASGAL-TOP100 sample, taking into account the mass and luminosity of the clumps, and the column densities of methyl acetylene (CH$_3$CCH), acetonitrile (CH$_3$CN), formaldehyde (H$_2$CO) and methanol (CH$_3$OH). We find that the chosen molecular tracers are affected by the thermal evolution of the clumps, showing steep ice evaporation gradients from 10$^3$ to 10$^5$ AU during the warm-up phase. We succeed in reproducing the observed column densities of CH$_3$CCH and CH$_3$CN, while H$_2$CO and CH$_3$OH show a poorer agreement with the observed values. The total (massive) star formation time is found to be $\sim5.2\times10^5$ yr, which is defined by the timescales of the individual evolutionary phases of the ATLASGAL-TOP100 sample: $\sim5\times10^4$ yr for 70-$μ$m weak, $\sim1.2\times10^5$ yr for mid-IR weak, $\sim2.4\times10^5$ yr for mid-IR bright and $\sim1.1\times10^5$ yr for HII-regions phases. Our models, with an appropriate selection of molecular tracers that can act as chemical clocks, allow to get robust estimates of the duration of the individual phases of the high-mass star formation process, with the advantage of being capable to include additional tracers aimed at increasing the accuracy of the estimated timescales. △ Less

Submitted 13 October, 2021; v1 submitted 1 June, 2021; originally announced June 2021.

Comments: Published on A&A (19 pages, 9 figures, 7 tables)

Journal ref: A&A 652, A71 (2021)

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)

The SEDIGISM survey: first data release and overview of the Galactic structure

Authors: F. Schuller, J. S. Urquhart, T. Csengeri, D. Colombo, A. Duarte-Cabral, M. Mattern, A. Ginsburg, A. R. Pettitt, F. Wyrowski, L. Anderson, F. Azagra, P. Barnes, M. Beltran, H. Beuther, S. Billington, L. Bronfman, R. Cesaroni, C. Dobbs, D. Eden, M. -Y. Lee, S. -N. X. Medina, K. M. Menten, T. Moore, F. M. Montenegro-Montes, S. Ragan , et al. (35 additional authors not shown)

Abstract: The SEDIGISM (Structure, Excitation and Dynamics of the Inner Galactic Interstellar Medium) survey used the APEX telescope to map 84 deg^2 of the Galactic plane between l = -60 deg and l = +31 deg in several molecular transitions, including 13CO(2-1) and C18O(2-1), thus probing the moderately dense (~10^3 cm^-3) component of the interstellar medium. With an angular resolution of 30'' and a typical… ▽ More The SEDIGISM (Structure, Excitation and Dynamics of the Inner Galactic Interstellar Medium) survey used the APEX telescope to map 84 deg^2 of the Galactic plane between l = -60 deg and l = +31 deg in several molecular transitions, including 13CO(2-1) and C18O(2-1), thus probing the moderately dense (~10^3 cm^-3) component of the interstellar medium. With an angular resolution of 30'' and a typical 1-sigma sensitivity of 0.8-1.0 K at 0.25 km/s velocity resolution, it gives access to a wide range of structures, from individual star-forming clumps to giant molecular clouds and complexes. The coverage includes a good fraction of the first and fourth Galactic quadrants, allowing us to constrain the large scale distribution of cold molecular gas in the inner Galaxy. In this paper we provide an updated overview of the full survey and the data reduction procedures used. We also assess the quality of these data and describe the data products that are being made publicly available as part of this first data release (DR1). We present integrated maps and position-velocity maps of the molecular gas and use these to investigate the correlation between the molecular gas and the large scale structural features of the Milky Way such as the spiral arms, Galactic bar and Galactic centre. We find that approximately 60 per cent of the molecular gas is associated with the spiral arms and these appear as strong intensity peaks in the derived Galactocentric distribution. We also find strong peaks in intensity at specific longitudes that correspond to the Galactic centre and well known star forming complexes, revealing that the 13CO emission is concentrated in a small number of complexes rather than evenly distributed along spiral arms. △ Less

Submitted 2 December, 2020; originally announced December 2020.

Comments: MNRAS, in press

The SEDIGISM survey: Molecular clouds in the inner Galaxy

Authors: A. Duarte-Cabral, D. Colombo, J. S. Urquhart, A. Ginsburg, D. Russeil, F. Schuller, L. D. Anderson, P. J. Barnes, M. T. Beltran, H. Beuther, S. Bontemps, L. Bronfman, T. Csengeri, C. L. Dobbs, D. Eden, A. Giannetti, J. Kauffmann, M. Mattern, S. -N. X. Medina, K. M. Menten, M. -Y. Lee, A. R. Pettitt, M. Riener, A. J. Rigby, A. Trafficante , et al. (35 additional authors not shown)

Abstract: We use the 13CO(2-1) emission from the SEDIGISM high-resolution spectral-line survey of the inner Galaxy, to extract the molecular cloud population with a large dynamic range in spatial scales, using the SCIMES algorithm. This work compiles a cloud catalogue with a total of 10663 molecular clouds, 10300 of which we were able to assign distances and compute physical properties. We study some of the… ▽ More We use the 13CO(2-1) emission from the SEDIGISM high-resolution spectral-line survey of the inner Galaxy, to extract the molecular cloud population with a large dynamic range in spatial scales, using the SCIMES algorithm. This work compiles a cloud catalogue with a total of 10663 molecular clouds, 10300 of which we were able to assign distances and compute physical properties. We study some of the global properties of clouds using a science sample, consisting of 6664 well resolved sources and for which the distance estimates are reliable. In particular, we compare the scaling relations retrieved from SEDIGISM to those of other surveys, and we explore the properties of clouds with and without high-mass star formation. Our results suggest that there is no single global property of a cloud that determines its ability to form massive stars, although we find combined trends of increasing mass, size, surface density and velocity dispersion for the sub-sample of clouds with ongoing high-mass star formation. We then isolate the most extreme clouds in the SEDIGISM sample (i.e. clouds in the tails of the distributions) to look at their overall Galactic distribution, in search for hints of environmental effects. We find that, for most properties, the Galactic distribution of the most extreme clouds is only marginally different to that of the global cloud population. The Galactic distribution of the largest clouds, the turbulent clouds and the high-mass star-forming clouds are those that deviate most significantly from the global cloud population. We also find that the least dynamically active clouds (with low velocity dispersion or low virial parameter) are situated further afield, mostly in the least populated areas. However, we suspect that part of these trends may be affected by some observational biases, and thus require further follow up work in order to be confirmed. △ Less

Submitted 2 December, 2020; originally announced December 2020.

Comments: 25 pages (+ appendices, 15 pages), 26 figures, MNRAS

SEDIGISM-ATLASGAL: Dense Gas Fraction and Star Formation Efficiency Across the Galactic Disk

Authors: J. S. Urquhart, C. Figura, J. R. Cross, M. R. A. Wells, T. J. T. Moore, D. J. Eden, S. E. Ragan, A. R. Pettitt, A. Duarte-Cabral, D. Colombo, F. Schuller, T. Csengeri, M. Mattern, H. Beuther, K. M. Menten, F. Wyrowski, L. D. Anderson, P. J. Barnes, M. T. Beltrán, S. J. Billington, L. Bronfman, A. Giannetti, J. Kainulainen, J. Kauffmann, M. -Y. Lee , et al. (10 additional authors not shown)

Abstract: By combining two surveys covering a large fraction of the molecular material in the Galactic disk we investigate the role the spiral arms play in the star formation process. We have matched clumps identified by ATLASGAL with their parental GMCs as identified by SEDIGISM, and use these giant molecular cloud (GMC) masses, the bolometric luminosities, and integrated clump masses obtained in a concurr… ▽ More By combining two surveys covering a large fraction of the molecular material in the Galactic disk we investigate the role the spiral arms play in the star formation process. We have matched clumps identified by ATLASGAL with their parental GMCs as identified by SEDIGISM, and use these giant molecular cloud (GMC) masses, the bolometric luminosities, and integrated clump masses obtained in a concurrent paper to estimate the dense gas fractions (DGF$_{\rm gmc}=\sum M_{\rm clump}/M_{\rm gmc}$) and the instantaneous star forming efficiencies (i.e., SFE$_{\rm gmc} = \sum L_{\rm clump}/M_{\rm gmc}$). We find that the molecular material associated with ATLASGAL clumps is concentrated in the spiral arms ($\sim$60% found within $\pm$10 km s$^{-1}$ of an arm). We have searched for variations in the values of these physical parameters with respect to their proximity to the spiral arms, but find no evidence for any enhancement that might be attributable to the spiral arms. The combined results from a number of similar studies based on different surveys indicate that, while spiral-arm location plays a role in cloud formation and HI to H$_2$ conversion, the subsequent star formation processes appear to depend more on local environment effects. This leads us to conclude that the enhanced star formation activity seen towards the spiral arms is the result of source crowding rather than the consequence of a any physical process. △ Less

Submitted 2 December, 2020; originally announced December 2020.

Comments: 15 pages, 12 figures. Full list of affiliations can be found at the end of the paper

ATLASGAL -- Relationship between dense star forming clumps and interstellar masers

Authors: S. J. Billington, J. S. Urquhart, C. König, H. Beuther, S. L. Breen, K. M. Menten, J. Campbell-White, S. P. Ellingsen, M. A. Thompson, T. J. T. Moore, D. J. Eden, W. -J. Kim, S. Leurini

Abstract: We have used catalogues from several Galactic plane surveys and dedicated observations to investigate the relationship between various maser species and Galactic star forming clumps, as identified by the ATLASGAL survey. The maser transitions of interest are the 6.7 & 12.2 GHz methanol masers, 22.2 GHz water masers, and the masers emitting in the four ground-state hyperfine structure transitions o… ▽ More We have used catalogues from several Galactic plane surveys and dedicated observations to investigate the relationship between various maser species and Galactic star forming clumps, as identified by the ATLASGAL survey. The maser transitions of interest are the 6.7 & 12.2 GHz methanol masers, 22.2 GHz water masers, and the masers emitting in the four ground-state hyperfine structure transitions of hydroxyl. We find clump association rates for the water, hydroxyl and methanol masers to be 56, 39 and 82 per cent respectively, within the Galactic longitude range of 60° > $l$ > -60°. We investigate the differences in physical parameters between maser associated clumps and the full ATLASGAL sample, and find that clumps coincident with maser emission are more compact with increased densities and luminosities. However, we find the physical conditions within the clumps are similar for the different maser species. A volume density threshold of $n$(H$_{2}$) > 10$^{4.1}$ cm$^{-3}$ for the 6.7 GHz methanol maser found in our previous study is shown to be consistent across for all maser species investigated. We find limits that are required for the production of maser emission to be 500 L$_{\odot}$ and 6 M$_{\odot}$ respectively. The evolutionary phase of maser associated clumps is investigated using the L/M ratio of clumps coincident with maser emission, and these have similar L/M ranges (~10$^{0.2}$ - 10$^{2.7}$ L$_{\odot}$/M$_{\odot}$) regardless of the associated transitions. This implies that the conditions required for the production of maser emission only occur during a relatively narrow period during a star's evolution. Lower limits of the statistical lifetimes for each maser species are derived, ranging from ~0.4 - 2 x 10$^{4}$ yrs and are in good agreement with the "straw man" evolutionary model previously presented. △ Less

Submitted 30 September, 2020; originally announced September 2020.

Comments: Accepted for publication in the Monthly Notices of the Royal Astronomical Society, 17 pages, 14 figures, 5 tables

36 GHz methanol lines from nearby galaxies: maser or quasi-thermal emission?

Authors: P. Humire, C. Henkel, Y. Gong, S. Leurini, R. Mauersberger, S. A. Levshakov, B. Winkel, A. Tarchi, P. Castangia, A. Malawi, H. Asiri, S. P. Ellingsen, T. P. McCarthy, X. Chen, X. Tang

Abstract: Methanol (CH3OH) is one of the most abundant interstellar molecules, offering a vast number of transitions to be studied, including many maser lines. While the strongest Galactic CH3OH lines, the so-called class II masers, show no indications for the presence of superluminous counterparts in external galaxies, the less luminous Galactic class I sources appear to be different. Here we report class… ▽ More Methanol (CH3OH) is one of the most abundant interstellar molecules, offering a vast number of transitions to be studied, including many maser lines. While the strongest Galactic CH3OH lines, the so-called class II masers, show no indications for the presence of superluminous counterparts in external galaxies, the less luminous Galactic class I sources appear to be different. Here we report class I 36GHz CH3OH 4(-1) - 3(0) E line emission from the nearby galaxies Maffei2 and IC342, measured with the 100-m telescope at Effelsberg at three different epochs within a time span of about five weeks. The 36GHz methanol line of Maffei2 is the second most luminous among the sources detected with certainty outside the Local Group of galaxies. This is not matched by the moderate infrared luminosity of Maffei2. Higher resolution data are required to check whether this is related to its prominent bar and associated shocks. Upper limits for M82, NGC4388, NGC5728 and Arp220 are also presented. The previously reported detection of 36GHz maser emission in Arp220 is not confirmed. Non-detections are reported from the related class I 44GHz methanol transition towards Maffei2 and IC342, indicating that this line is not stronger than its 36GHz counterpart. In contrast to the previously detected 36GHz CH3OH emission in NGC253 and NGC4945, our 36GHz profiles towards Maffei2 and IC342 are similar to those of previously detected non-masing lines from other molecular species. However, by analogy to our Galactic center region, it may well be possible that the 36GHz methanol lines in Maffei~2 and IC~342 are composed of a large number of faint and narrow maser features that remain spatially unresolved. In view of this, a search for a weak broad 36GHz line component would also be desirable in NGC253 and NGC4945. △ Less

Submitted 2 December, 2019; v1 submitted 15 November, 2019; originally announced November 2019.

Comments: Accepted by Astronomy & Astrophysics, 8 pages, 1 table, 3 figures, slight changes w.r.t. version 1 due to language editing

Journal ref: A&A 633, A106 (2020)

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)

On the size of the CO-depletion radius in the IRDC G351.77-0.51

Authors: G. Sabatini, A. Giannetti, S. Bovino, J. Brand, S. Leurini, E. Schisano, T. Pillai, K. M. Menten

Abstract: An estimate of the degree of CO-depletion ($f_D$) provides information on the physical conditions occurring in the innermost and densest regions of molecular clouds. A key parameter in these studies is the size of the depletion radius, i.e. the radius within which the C-bearing species, and in particular CO, are largely frozen onto dust grains. A strong depletion state (i.e. $f_D>10$, as assumed i… ▽ More An estimate of the degree of CO-depletion ($f_D$) provides information on the physical conditions occurring in the innermost and densest regions of molecular clouds. A key parameter in these studies is the size of the depletion radius, i.e. the radius within which the C-bearing species, and in particular CO, are largely frozen onto dust grains. A strong depletion state (i.e. $f_D>10$, as assumed in our models) is highly favoured in the innermost regions of dark clouds, where the temperature is $<20$ K and the number density of molecular hydrogen exceeds a few $\times$10$^{4}$ cm$^{-3}$. In this work, we estimate the size of the depleted region by studying the Infrared Dark Cloud (IRDC) G351.77-0.51. Continuum observations performed with the $Herschel$ $Space$ $Observatory$ and the $LArge$ $APEX$ $BOlometer$ $CAmera$, together with APEX C$^{18}$O and C$^{17}$O J=2$\rightarrow$1 line observations, allowed us to recover the large-scale beam- and line-of-sight-averaged depletion map of the cloud. We built a simple model to investigate the depletion in the inner regions of the clumps in the filament and the filament itself. The model suggests that the depletion radius ranges from 0.02 to 0.15 pc, comparable with the typical filament width (i.e.$\sim$0.1 pc). At these radii, the number density of H$_2$ reaches values between 0.2 and 5.5$\times$10$^{5}$ cm$^{-3}$. These results provide information on the approximate spatial scales on which different chemical processes operate in high-mass star-forming regions and also suggest caution when using CO for kinematical studies in IRDCs. △ Less

Submitted 7 October, 2019; originally announced October 2019.

Comments: accepted to MNRAS (14 pages, 10 figures, 1 table)

Fragmentation, rotation and outflows in the high-mass star-forming region IRAS 23033+5951. A case study of the IRAM NOEMA large program CORE

Authors: F. Bosco, H. Beuther, A. Ahmadi, J. C. Mottram, R. Kuiper, H. Linz, L. Maud, J. M. Winters, T. Henning, S. Feng, T. Peters, D. Semenov, P. D. Klaassen, P. Schilke, J. S. Urquhart, M. T. Beltrán, S. L. Lumsden, S. Leurini, L. Moscadelli, R. Cesaroni, Á. Sánchez-Monge, A. Palau, R. Pudritz, F. Wyrowski, S. Longmore , et al. (1 additional authors not shown)

Abstract: The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. U… ▽ More The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. Using the Northern Extended Millimeter Array (NOEMA) in three configurations and the IRAM 30-m single-dish telescope at 220GHz, we probe the gas and dust emission at an angular resolution of ~0.45arcsec, corresponding to 1900au. In the mm continuum emission, we identify a protostellar cluster with at least four mm-sources, where three of them show a significantly higher peak intensity well above a signal-to-noise ratio of 100. Hierarchical fragmentation from large to small spatial scales is discussed. Two fragments are embedded in rotating structures and drive molecular outflows, traced by $^{13}$CO (2-1) emission. The velocity profiles across two of the cores are similar to Keplerian but are missing the highest velocity components close to the center of rotation, which is a common phenomena from observations like these, and other rotation scenarios are not excluded entirely. Position-velocity diagrams suggest protostellar masses of ~6 and 19M$_\sun$. Rotational temperatures from fitting CH$_3$CN ($12_K-11_K$) spectra are used for estimating the gas temperature and by that the disk stability against gravitational fragmentation, utilizing Toomre's $Q$ parameter. [We] identify only one candidate disk to be unstable against gravitational instability caused by axisymmetric perturbations. The dominant sources cover different evolutionary stages within the same maternal gas clump. The appearance of rotation and outflows of the cores are similar to those found in low-mass star-forming regions. △ Less

Submitted 9 July, 2019; originally announced July 2019.

Comments: 23 pages, 16 figures, 5 tables, accepted for publication in A&A

Journal ref: A&A 629, A10 (2019)

ATLASGAL -- physical parameters of dust clumps associated with 6.7 GHz methanol masers

Authors: S. J. Billington, J. S. Urquhart, C. Konig, T. J. T. Moore, D. J. Eden, S. L. Breen, W. -J. Kim, M. A. Thompson, S. P. Ellingsen, K. M. Menten, F. Wyrowski, S. Leurini

Abstract: We have constructed the largest sample of dust-associated class II 6.7 GHz methanol masers yet obtained. New measurements from the the Methanol MultiBeam (MMB) Survey were combined with the 870 $μ$m APEX Telescope Large Area Survey of the Galaxy (ATLASGAL) and the 850 $μ$m JCMT Plane Survey (JPS). Together with two previous studies we have now identified the host clumps for 958 methanol masers acr… ▽ More We have constructed the largest sample of dust-associated class II 6.7 GHz methanol masers yet obtained. New measurements from the the Methanol MultiBeam (MMB) Survey were combined with the 870 $μ$m APEX Telescope Large Area Survey of the Galaxy (ATLASGAL) and the 850 $μ$m JCMT Plane Survey (JPS). Together with two previous studies we have now identified the host clumps for 958 methanol masers across the Galactic Plane, covering approximately 99% of the MMB catalogue and increasing the known sample of dust-associated masers by over 30%. We investigate correlations between the physical properties of the clumps and masers using distances and luminosities drawn from the literature. Clumps hosting methanol masers are significantly more compact and have higher volume densities than the general population of clumps. We determine a minimum volume density threshold of $n$(H$_2$) $\geq 10^4$ cm$^{-3}$ for the efficient formation of intermediate- and high-mass stars. We find 6.7 GHz methanol masers are associated with a distinct part of the evolutionary process ($L_{\rm bol}$/$M_{\rm fwhm}$ ratios of between 10$^{0.6}$ and 10$^{2.2}$) and have well defined turning on and termination points. We estimate the lower limit for the mass of embedded objects to be $\geq$ 6 M$_{\odot}$ and the statistical lifetime of the methanol maser stage to be $\sim$ 3.3$\times$10$^{4}$ yrs. This suggests that methanol masers are indeed reliable tracers of high mass star formation, and indicates that the evolutionary period traced by this marker is relatively rapid. △ Less

Submitted 23 September, 2019; v1 submitted 1 July, 2019; originally announced July 2019.

Comments: Accepted for publication in the Monthly Notices of the Royal Astronomical Society, 20 pages, 23 figures, 6 tables

Warm gas in protostellar outflows. II. EHV jet and outflows from OMC-2/3

Authors: A. I. Gómez-Ruiz, A. Gusdorf, S. Leurini, K. M. Menten, S. Takahashi, F. Wyrowski, R. Güsten

Abstract: OMC-2/3 is one of the nearest embedded cluster-forming region that includes intermediate-mass protostars at early stages of evolution. A previous CO (3--2) mapping survey towards this region revealed outflow activity related to sources at different evolutionary phases. The present work aims to study the warm gas in the high-velocity emission from several outflows found in CO (3--2) emission by pre… ▽ More OMC-2/3 is one of the nearest embedded cluster-forming region that includes intermediate-mass protostars at early stages of evolution. A previous CO (3--2) mapping survey towards this region revealed outflow activity related to sources at different evolutionary phases. The present work aims to study the warm gas in the high-velocity emission from several outflows found in CO (3--2) emission by previous observations, determine their physical conditions, and make comparison with previous results in low-mass star-forming regions. We use the CHAMP+ heterodyne array on the APEX telescope to map the CO (6--5) and CO (7--6) emission in the OMC-2 FIR 6 and OMC-3 MMS 1-6 regions, and to observe $^{13}$CO (6--5) at selected positions. We analyze these data together with previous CO (3--2) observations. In addition, we mapped the SiO (5--4) emission in OMC-2 FIR 6. The CO (6--5) emission was detected in most of the outflow lobes in the mapped regions, while the CO (7--6) was found mostly from the OMC-3 outflows. In the OMC-3 MMS 5 outflow, a previously undetected extremely high velocity gas was found in CO (6--5). This extremely high velocity emission arises from the regions close to the central object MMS 5. Radiative transfer models revealed that the high-velocity gas from MMS 5 outflow consists of gas with $n_{\rm H_2}=$10$^4$--10$^{5}$ cm$^{-3}$ and $T>$ 200 K, similar to what is observed in young Class 0 low-mass protostars. For the other outflows, values of $n_{\rm H_2}>$10$^{4}$ cm$^{-3}$ were found. The physical conditions and kinematic properties of the young intermediate-mass outflows presented here are similar to those found in outflows from Class 0 low-mass objects. Due to their excitation requirements, mid$-J$ CO lines are good tracers of extremely high velocity gas in young outflows likely related to jets. △ Less

Submitted 8 June, 2019; originally announced June 2019.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 629, A77 (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)

First interferometric study of enhanced N-fractionation in N$_{2}$H$^{+}$: the high-mass star-forming region IRAS 05358+3543

Authors: L. Colzi, F. Fontani, P. Caselli, S. Leurini, L. Bizzocchi, G. Quaia

Abstract: Nitrogen (N) fractionation is used as a tool to search for a link between the chemical history of the Solar System and star-forming regions. A large variation of $^{14}$N/$^{15}$N is observed towards different astrophysical sources, and current chemical models cannot reproduce it. With the advent of high angular resolution radiotelescopes it is now possible to search for N-fractionation at core sc… ▽ More Nitrogen (N) fractionation is used as a tool to search for a link between the chemical history of the Solar System and star-forming regions. A large variation of $^{14}$N/$^{15}$N is observed towards different astrophysical sources, and current chemical models cannot reproduce it. With the advent of high angular resolution radiotelescopes it is now possible to search for N-fractionation at core scales. We present IRAM NOEMA observations of the J=1-0 transition of N$_{2}$H$^{+}$, $^{15}$NNH$^{+}$ and N$^{15}$NNH$^{+}$ towards the high-mass protocluster IRAS 05358+3543. We find $^{14}$N/$^{15}$N ratios that span from $\sim$100 up to $\sim$220 and these values are lower or equal than those observed with single-dish observations towards the same source. Since N-fractionation changes across the studied region, this means that it is regulated by local environmental effects. We find also the possibility, for one of the four cores defined in the protocluster, to have a more abundant $^{15}$NNH$^{+}$ with respect to N$^{15}$NNH$^{+}$. This is another indication that current chemical models may be missing chemical reactions or may not take into account other mechanisms, like photodissociation or grain surface chemistry, that could be important. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: 19 pages, 8 figures, 6 tables, 3 appendices Accepted in Monthly Notices of the Royal Astronomical Society Letters

Massive and low-mass protostars in massive "starless" cores

Authors: Thushara Pillai, Jens Kauffmann, Qizhou Zhang, Patricio Sanhueza, Silvia Leurini, Ke Wang, T. K. Sridharan, Carsten König

Abstract: The infrared dark clouds (IRDCs) G11.11$-$0.12 and G28.34$+$0.06 are two of the best-studied IRDCs in our Galaxy. These two clouds host clumps at different stages of evolution, including a massive dense clump in both clouds that is dark even at 70 and 100$μ$m. Such seemingly quiescent massive dense clumps have been speculated to harbor cores that are precursors of high-mass stars and clusters. We… ▽ More The infrared dark clouds (IRDCs) G11.11$-$0.12 and G28.34$+$0.06 are two of the best-studied IRDCs in our Galaxy. These two clouds host clumps at different stages of evolution, including a massive dense clump in both clouds that is dark even at 70 and 100$μ$m. Such seemingly quiescent massive dense clumps have been speculated to harbor cores that are precursors of high-mass stars and clusters. We observed these two "prestellar" regions at 1mm with the Submillimeter Array (SMA) with the aim of characterizing the nature of such cores. We show that the clumps fragment into several low- to high-mass cores within the filamentary structure of the enveloping cloud. However, while the overall physical properties of the clump may indicate a starless phase, we find that both regions host multiple outflows. The most massive core though 70 $μ$m dark in both clumps is clearly associated with compact outflows. Such low-luminosity, massive cores are potentially the earliest stage in the evolution of a massive protostar. We also identify several outflow features distributed in the large environment around the most massive core. We infer that these outflows are being powered by young, low-mass protostars whose core mass is below our detection limit. These findings suggest that low-mass protostars have already formed or are coevally formed at the earliest phase of high-mass star formation. △ Less

Submitted 22 January, 2019; originally announced January 2019.

Comments: in print at A&A

ATLASGAL --- Molecular fingerprints of a sample of massive star forming clumps

Authors: J. S. Urquhart, C. Figura, F. Wyrowski, A. Giannetti, W. -J. Kim, M. Wienen, S. Leurini, T. Pillai, T. Csengeri, S. J. Gibson, K. Menten, T. J. T. Moore, M. A. Thompson

Abstract: We have conducted a 3-mm molecular-line survey towards 570 high-mass star-forming clumps, using the Mopra telescope. The sample is selected from the 10,000 clumps identified by the ATLASGAL survey and includes all of the most important embedded evolutionary stages associated with massive star formation, classified into five distinct categories (quiescent, protostellar, young stellar objects, \hii\… ▽ More We have conducted a 3-mm molecular-line survey towards 570 high-mass star-forming clumps, using the Mopra telescope. The sample is selected from the 10,000 clumps identified by the ATLASGAL survey and includes all of the most important embedded evolutionary stages associated with massive star formation, classified into five distinct categories (quiescent, protostellar, young stellar objects, \hii\ regions and photo-dominated regions). The observations were performed in broadband mode with frequency coverage of 85.2 to 93.4\,GHz and a velocity resolution of $\sim$0.9\,\kms, detecting emission from 26 different transitions. We find significant evolutionary trends in the detection rates, integrated line intensities, and abundances of many of the transitions and also identify a couple of molecules that appear to be invariant to changes in the dust temperature and evolutionary stage (N$_2$H$^+$\,(1-0) and HN$^{13}$C\,(1-0)). We use the K-ladders for CH$_3$C$_2$H\,(5-4) and CH$_3$CH\,(5-4) to calculate the rotation temperatures and find $\sim$1/3 of the quiescent clumps have rotation temperatures that suggest the presence of an internal heating source. These sources may constitute a population of very young protostellar objects that are still dark at 70\,\mum\ and suggest that the fraction of truly quiescent clumps may only be a few per cent. We also identify a number of line ratios that show a strong correlation with the evolutionary stage of the embedded objects and discuss their utility as diagnostic probes of evolution. △ Less

Submitted 11 January, 2019; originally announced January 2019.

Comments: Accepted for publication in MNRAS; consists of 28 pages, 22 figures and 10 tables. For full versions of Tables 1, 5 and 10 will only be available via CDS

A timeline for massive star-forming regions via combined observation of o-H$_2$D$^+$ and N$_2$D$^+$

Authors: A. Giannetti, S. Bovino, P. Caselli, S. Leurini, D. R. G. Schleicher, B. Koertgen, K. M. Menten, T. Pillai, F. Wyrowski

Abstract: Context: In cold and dense gas prior to the formation of young stellar objects, heavy molecular species (including CO) are accreted onto dust grains. Under these conditions H$_3^+$ and its deuterated isotopologues become more abundant, enhancing the deuterium fraction of molecules such as N$_2$H$^+$ that are formed via ion-neutral reactions. Because this process is extremely temperature sensitive,… ▽ More Context: In cold and dense gas prior to the formation of young stellar objects, heavy molecular species (including CO) are accreted onto dust grains. Under these conditions H$_3^+$ and its deuterated isotopologues become more abundant, enhancing the deuterium fraction of molecules such as N$_2$H$^+$ that are formed via ion-neutral reactions. Because this process is extremely temperature sensitive, the abundance of these species is likely linked to the evolutionary stage of the source. Aims: We investigate how the abundances of o-H$_2$D$^+$ and N$_2$D$^+$ vary with evolution in high-mass clumps. Methods: We observed with APEX the ground-state transitions of o-H$_2$D$^+$ near 372 GHz, and N$_2$D$^+$(3-2) near 231 GHz for three massive clumps in different evolutionary stages. The sources were selected within the G351.77-0.51 complex to minimise the variation of initial chemical conditions, and to remove distance effects. We modelled their dust continuum emission to estimate their physical properties, and also modelled their spectra under the assumption of local thermodynamic equilibrium to calculate beam-averaged abundances. Results: We find an anticorrelation between the abundance of o-H$_2$D$^+$ and that of N$_2$D$^+$, with the former decreasing and the latter increasing with evolution. With the new observations we are also able to provide a qualitative upper limit to the age of the youngest clump of about 10$^5$ yr, comparable to its current free-fall time. Conclusions: We can explain the evolution of the two tracers with simple considerations on the chemical formation paths, depletion of heavy elements, and evaporation from the grains. We therefore propose that the joint observation and the relative abundance of o-H$_2$D$^+$ and N$_2$D$^+$ can act as an efficient tracer of the evolutionary stages of the star-formation process. △ Less

Submitted 7 January, 2019; originally announced January 2019.

ATLASGAL-selected high-mass clumps in the inner Galaxy. VII. Characterisation of mid-J CO emission

Authors: Felipe Navarete, Silvia Leurini, Andrea Giannetti, Friedrich Wyrowski, James S. Urquhart, Carsten Koenig, Timea Csengeri, Rolf Guesten, Augusto Damineli, Karl M. Menten

Abstract: High-mass stars are formed within massive molecular clumps, where a large number of stars form close together. The evolution of the clumps with different masses and luminosities is mainly regulated by its high-mass stellar content and the formation of such objects is still not well understood. In this work, we characterise the mid-J CO emission in a statistical sample of 99 clumps (Top100) selecte… ▽ More High-mass stars are formed within massive molecular clumps, where a large number of stars form close together. The evolution of the clumps with different masses and luminosities is mainly regulated by its high-mass stellar content and the formation of such objects is still not well understood. In this work, we characterise the mid-J CO emission in a statistical sample of 99 clumps (Top100) selected from the ATLASGAL survey that are representative of the Galactic proto-cluster population. High-spatial resolution APEX-CHAMP+ maps of the CO(6-5) and CO(7-6) transitions were obtained and combined with additional single-pointing APEX-FLASH+ spectra of the CO(4-3) line. We study the correlations of the CO line luminosities and profiles for the three CO transitions with the clump properties and investigate if and how they change as a function of the evolution. All sources were detected above 3-$σ$ in all three CO transitions and most of the sources exhibit broad CO emission likely associated with molecular outflows. We found that the extension of the mid-J CO emission is correlated with the size of the dust emission traced by the Herschel-PACS 70 $μ$m maps. The CO line luminosity is correlated with the luminosity and mass of the clumps. However, it does not correlate with the L/M ratio. The dependency of the CO luminosity with the properties of the clumps is steeper for higher-J transitions. Our data seem to exclude that this trend is biased by self-absorption features in the CO emission, but rather suggest that different J transitions arise from different regions of the inner envelope. Moreover, high-mass clumps show similar trends in CO luminosity as lower mass clumps, but are systematically offset towards larger values, suggesting that higher column density and/or temperature (of unresolved) CO emitters are found inside high-mass clumps. △ Less

Submitted 6 January, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: Accepted by A&A: 62 pages, 230 figures. ArXiV: 37 pages, 40 figures

Journal ref: A&A 622, A135 (2019)

Characterising the high-mass star forming filament G351.776--0.527 with Herschel and APEX dust continuum and gas observations

Authors: S. Leurini, E. Schisano, T. Pillai, A. Giannetti, J. Urquhart, T. Csengeri, S. Casu, M. Cunningham, D. Elia, P. A. Jones, C. Koenig, S. Molinari, T. Stanke, L. Testi, F. Wyrowski, K. M. Menten

Abstract: G351.776-0.527 is among the most massive, closest, and youngest filaments in the inner Galactic plane and therefore it is an ideal laboratory to study the kinematics of dense gas and mass replenishment on a large scale. In this paper, we present far-infrared (FIR) and submillimetre wavelength continuum observations combined with spectroscopic C$^{18}$O (2-1) data of the entire region to study its… ▽ More G351.776-0.527 is among the most massive, closest, and youngest filaments in the inner Galactic plane and therefore it is an ideal laboratory to study the kinematics of dense gas and mass replenishment on a large scale. In this paper, we present far-infrared (FIR) and submillimetre wavelength continuum observations combined with spectroscopic C$^{18}$O (2-1) data of the entire region to study its temperature, mass distribution, and kinematics. The structure is composed of a main elongated region with an aspect ratio of $\sim 23$, which is associated with a network of filamentary structures. The main filament has a remarkably constant width of 0.2 pc. The total mass of the network (including the main filament) is $\geq 2600$ M$_\odot$, while we estimate a mass of $\sim 2000$ M$_\odot$ for the main structure. Therefore, the network harbours a large reservoir of gas and dust that could still be accreted onto the main structure. From the analysis of the gas kinematics, we detect two velocity components in the northern part of the main filament. The data also reveal velocity oscillations in C$^{18}$O along the spine in the main filament and in at least one of the branches. Considering the region as a single structure, we find that it is globally close to virial equilibrium indicating that the entire structure is approximately in a stable state. △ Less

Submitted 3 December, 2018; originally announced December 2018.

Comments: 18 pages, 19 figures, accepted for publication in A&A

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

The Formation Conditions of the Wide Binary Class 0 Protostars within BHR 71

Authors: John Tobin, Tyler Bourke, Stacy Mader, Lars Kristensen, Hector Arce, Frederic Gueth, Antoine Gusdorf, Claudio Codella, Silvia Leurini, Xuepeng Chen

Abstract: We present a characterization of the binary protostar system that is forming within a dense core in the isolated dark cloud BHR71. The pair of protostars, IRS1 and IRS2, are both in the Class 0 phase, determined from observations that resolve the sources from 1 um out to 250 um and from 1.3 mm to 1.3cm. The resolved observations enable the luminosities of IRS1 and IRS2 to be independently measured… ▽ More We present a characterization of the binary protostar system that is forming within a dense core in the isolated dark cloud BHR71. The pair of protostars, IRS1 and IRS2, are both in the Class 0 phase, determined from observations that resolve the sources from 1 um out to 250 um and from 1.3 mm to 1.3cm. The resolved observations enable the luminosities of IRS1 and IRS2 to be independently measured (14.7 and 1.7L_sun, respectively), in addition to the bolometric temperatures 68~K, and 38~K, respectively. The surrounding core was mapped in NH3 (1,1) with the Parkes radio telescope, and followed with higher-resolution observations from ATCA in NH3 (1,1) and 1.3cm continuum. The protostars were then further characterized with ALMA observations in the 1.3~mm continuum along with N2D+ (J=3-2), 12CO, 13CO, and C18O (J=2-1) molecular lines. The Parkes observations find evidence for a velocity gradient across the core surrounding the two protostars, while ATCA reveals more complex velocity structure toward the protostars within the large-scale gradient. The ALMA observations then reveal that the two protostars are at the same velocity in C18O, and N2H+ exhibits a similar velocity structure as NH3. However, the C18O kinematics reveal that the rotation on scales $<$1000~AU around IRS1 and IRS2 are in opposite directions. Taken with the lack of a systematic velocity difference between the pair, it is unlikely that their formation resulted from rotational fragmentation. We instead conclude that the binary system most likely formed via turbulent fragmentation of the core. △ Less

Submitted 7 November, 2018; originally announced November 2018.

Comments: 51 pages, 25 Figures, 2 Tables, Accepted to ApJ

SEDIGISM: The kinematics of ATLASGAL filaments

Authors: M. Mattern, J. Kauffmann, T. Csengeri, J. S. Urquhart, S. Leurini, F. Wyrowski, A. Giannetti, P. J. Barnes, H. Beuther, L. Bronfman, A. Duarte-Cabral, T. Henning, J. Kainulainen, K. M. Menten, E. Schisano, F. Schuller

Abstract: Analysing the kinematics of filamentary molecular clouds is a crucial step towards understanding their role in the star formation process. Therefore, we study the kinematics of 283 filament candidates in the inner Galaxy, that were previously identified in the ATLASGAL dust continuum data. The $^{13}$CO(2 - 1) and C$^{18}$O(2 - 1) data of the SEDIGISM survey (Structure, Excitation, and Dynamics of… ▽ More Analysing the kinematics of filamentary molecular clouds is a crucial step towards understanding their role in the star formation process. Therefore, we study the kinematics of 283 filament candidates in the inner Galaxy, that were previously identified in the ATLASGAL dust continuum data. The $^{13}$CO(2 - 1) and C$^{18}$O(2 - 1) data of the SEDIGISM survey (Structure, Excitation, and Dynamics of the Inner Galactic Inter Stellar Medium) allows us to analyse the kinematics of these targets and to determine their physical properties at a resolution of 30 arcsec and 0.25 km/s. To do so, we developed an automated algorithm to identify all velocity components along the line-of-sight correlated with the ATLASGAL dust emission, and derive size, mass, and kinematic properties for all velocity components. We find two-third of the filament candidates are coherent structures in position-position-velocity space. The remaining candidates appear to be the result of a superposition of two or three filamentary structures along the line-of-sight. At the resolution of the data, on average the filaments are in agreement with Plummer-like radial density profiles with a power-law exponent of p = 1.5 +- 0.5, indicating that they are typically embedded in a molecular cloud and do not have a well-defined outer radius. Also, we find a correlation between the observed mass per unit length and the velocity dispersion of the filament of $m \sim σ_v^2$. We show that this relation can be explained by a virial balance between self-gravity and pressure. Another possible explanation could be radial collapse of the filament, where we can exclude infall motions close to the free-fall velocity. △ Less

Submitted 22 August, 2018; originally announced August 2018.

Comments: 45 pages, 20 pages for Table of derived parameters

Journal ref: A&A 619, A166 (2018)

Core fragmentation and Toomre stability analysis of W3(H2O): A case study of the IRAM NOEMA large program CORE

Authors: A. Ahmadi, H. Beuther, J. C. Mottram, F. Bosco, H. Linz, Th. Henning, J. M. Winters, R. Kuiper, R. Pudritz, Á. Sánchez-Monge, E. Keto, M. Beltran, S. Bontemps, R. Cesaroni, T. Csengeri, S. Feng, R. Galvan-Madrid, K. G. Johnston, P. Klaassen, S. Leurini, S. N. Longmore, S. Lumsden, L. T. Maud, K. M. Menten, L. Moscadelli , et al. (8 additional authors not shown)

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

Submitted 1 August, 2018; originally announced August 2018.

Comments: 23 pages, 26 figures, accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 618, A46 (2018)

Water and interstellar complex organics associated with the HH 212 protostellar disc - On disc atmospheres, disc winds, and accretion shocks

Authors: C. Codella, E. Bianchi, B. Tabone, C. -F. Lee, S. Cabrit, C. Ceccarelli, L. Podio, F. Bacciotti, R. Bachiller, E. Chapillon, F. Gueth, A. Gusdorf, B. Lefloch, S. Leurini, G. Pineau des Forets, K. L. J. Rygl, M. Tafalla

Abstract: The HH 212 protostellar system, in Orion B, has been mapped thanks to ALMA-Band 7 Cycle 1 and Cycle 4 observations of dueterated water (HDO) and acetaldehyde (CH$_3$CHO) emission with an angular resolution down to $\sim$0.15 arcsec (60 au). Many emission lines due to 14 CH$_3$CHO and 1 HDO transitions at high excitation ($E_{\rm u}$ between 163 K and 335 K) have been imaged in the inner $\sim$ 70… ▽ More The HH 212 protostellar system, in Orion B, has been mapped thanks to ALMA-Band 7 Cycle 1 and Cycle 4 observations of dueterated water (HDO) and acetaldehyde (CH$_3$CHO) emission with an angular resolution down to $\sim$0.15 arcsec (60 au). Many emission lines due to 14 CH$_3$CHO and 1 HDO transitions at high excitation ($E_{\rm u}$ between 163 K and 335 K) have been imaged in the inner $\sim$ 70 au region. The local thermal equilibrium analysis of the CH$_3$CHO emission leads to a temperature of 78$\pm$14 K and a column density of 7.6$\pm$3.2 $\times$ 10$^{15}$ cm$^{-2}$, which, when $N_{\rm H_2}$ of 10$^{24}$ cm$^{-2}$ is assumed, leads to an abundance of $X_{\rm CH_3CHO}$ $\simeq$ 8 $\times$ 10$^{-9}$. The large velocity gradient analysis of the HDO emission also places severe constraints on the volume density, n$_{\rm H_2}$ $\geq$ 10$^8$ cm$^{-3}$. The line profiles are 5--7 km s$^{-1}$ wide, and CH$_3$CHO and HDO both show a $\pm$ 2 km s$^{-1}$ velocity gradient over a size of $\sim$ 70 au (blue-shifted emission towards the north-west and red-shifted emission towards the south-east) along the disc equatorial plane, in agreement with what was found so far using other molecular tracers. The kinematics of CH$_3$CHO and HDO are consistent with the occurrence of a centrifugal barrier, that is, the infalling envelope-rotating disc ring, which is chemically enriched through low-velocity accretion shocks. The emission radius is $\sim$ 60 au, in good agreement with what was found before for another interstellar complex organic molecule such as NH$_2$CHO. We support a vertical structure for the centrifugal barrier, suggesting the occurrence of two outflowing, expanding, and rotating rings above and below (of about 40-45 au) the optically thick equatorial disc plane. It is tempting to speculate that these rings could probe the basis of a wind launched from this region. △ Less

Submitted 20 June, 2018; originally announced June 2018.

Comments: Astronomy & Astrophysics, in press

Discovery of a sub-Keplerian disk with jet around a 20Msun young star. ALMA observations of G023.01-00.41

Authors: A. Sanna, A. Koelligan, L. Moscadelli, R. Kuiper, R. Cesaroni, T. Pillai, K. M. Menten, Q. Zhang, A. Caratti o Garatti, C. Goddi, S. Leurini, C. Carrasco-Gonzalez

Abstract: It is well established that Solar-mass stars gain mass via disk accretion, until the mass reservoir of the disk is exhausted and dispersed, or condenses into planetesimals. Accretion disks are intimately coupled with mass ejection via polar cavities, in the form of jets and less collimated winds, which allow mass accretion through the disk by removing a substantial fraction of its angular momentum… ▽ More It is well established that Solar-mass stars gain mass via disk accretion, until the mass reservoir of the disk is exhausted and dispersed, or condenses into planetesimals. Accretion disks are intimately coupled with mass ejection via polar cavities, in the form of jets and less collimated winds, which allow mass accretion through the disk by removing a substantial fraction of its angular momentum. Whether disk accretion is the mechanism leading to the formation of stars with much higher masses is still unclear. Here, we are able to build a comprehensive picture for the formation of an O-type star, by directly imaging a molecular disk which rotates and undergoes infall around the central star, and drives a molecular jet which arises from the inner disk regions. The accretion disk is truncated between 2000-3000au, it has a mass of about a tenth of the central star mass, and is infalling towards the central star at a high rate (6x10^-4 Msun/yr), as to build up a very massive object. These findings, obtained with the Atacama Large Millimeter/submillimeter Array at 700au resolution, provide observational proof that young massive stars can form via disk accretion much like Solar-mass stars. △ Less

Submitted 15 January, 2019; v1 submitted 24 May, 2018; originally announced May 2018.

Comments: 12 pages, 9 figures, 2 tables, minor revision included, accepted by Astronomy & Astrophysics

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

Fragmentation and disk formation during high-mass star formation: The IRAM NOEMA (Northern Extended Millimeter Array) large program CORE

Authors: H. Beuther, J. C. Mottram, A. Ahmadi, F. Bosco, H. Linz, Th. Henning, P. Klaassen, J. M. Winters, L. T. Maud, R. Kuiper, D. Semenov, C. Gieser, T. Peters, J. S. Urquhart, R. Pudritz, S. E. Ragan, S. Feng, E. Keto, S. Leurini, R. Cesaroni, M. Beltran, A. Palau, A. Sanchez-Monge, R. Galvan-Madrid, Q. Zhang , et al. (8 additional authors not shown)

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

Submitted 3 May, 2018; originally announced May 2018.

Comments: 29 pages, 17 figures, Astronomy & Astrophysics in press, for a higher-resolution version of the paper see http://www.mpia.de/homes/beuther/papers.html

Journal ref: A&A 617, A100 (2018)

The search for high-mass protostars with ALMA revealed up to kilo-parsec scales (SPARKS): I. Indication for a centrifugal barrier in the environment of a single high-mass envelope

Authors: T. Csengeri, S. Bontemps, F. Wyrowski, A. Belloche, K. M. Menten, S. Leurini, H. Beuther, L. Bronfman, B. Commerccon, E. Chapillon, S. Longmore, A. Palau, J. C. Tan, J. S. Urquhart

Abstract: The formation of the most massive O-type stars is poorly understood. We present a case study of a young massive clump from the ATLASGAL survey, G328.2551-0.5321. It exhibits a bolometric luminosity of 1.3$\times$10$^4$ L$_{\odot}$ corresponding to a current protostellar mass of $\sim$11 and 16 M$_{\odot}$. We analyze high angular-resolution observations with ALMA at $\sim$0.17" corresponding a phy… ▽ More The formation of the most massive O-type stars is poorly understood. We present a case study of a young massive clump from the ATLASGAL survey, G328.2551-0.5321. It exhibits a bolometric luminosity of 1.3$\times$10$^4$ L$_{\odot}$ corresponding to a current protostellar mass of $\sim$11 and 16 M$_{\odot}$. We analyze high angular-resolution observations with ALMA at $\sim$0.17" corresponding a physical scale of $\sim$400 au in dust continuum and molecular lines. The dust continuum emission reveals a single high-mass protostellar envelope and shows evidence for a marginally resolved continuum source. We detect a rotational line of CH$_3$OH within its $v_{\rm t}$=1 torsionally excited state revealing two bright peaks of emission spatially offset from the dust continuum peak, and exhibiting a distinct velocity component $\pm$4.5 km s$^{-1}$ offset compared to the source $v_{\rm lsr}$. Local thermodynamic equilibrium analysis suggests N(CH$_3$OH)=1.2$-$2$\times$10$^{19}$ cm$^{-2}$, and kinetic temperatures of 160$-$200 K at the position of these peaks. Their velocity shifts correspond well to the expected Keplerian velocity around a central object with 15M$_{\odot}$ consistent with the mass estimate based on the source's bolometric luminosity. We propose a picture where the CH$_3$OH emission peaks trace the accretion shocks around the centrifugal barrier, pinpointing the interaction region between the collapsing envelope and an accretion disk. Because the HC$_3$N $v_{\rm 7}$=1e ($J$=38-37) line shows compact emission, and a velocity pattern consistent with models of Keplerian rotation, we suggest that this could be a new tracer for compact accretion disks around high-mass protostars. The estimated physical properties of the accretion disk suggest a specific angular momentum several times larger than typically observed towards low-mass protostars. △ Less

Submitted 17 April, 2018; originally announced April 2018.

Comments: Manuscript submitted to A&A, and revised after first referee report (abstract abbreviated)

Journal ref: A&A 617, A89 (2018)

ATLASGAL-selected massive clumps in the inner Galaxy: VI. Kinetic temperature and spatial density measured with formaldehyde

Authors: X. D. Tang, C. Henkel, F. Wyrowski, A. Giannetti, K. M. Menten, T. Csengeri, S. Leurini, J. S. Urquhart, C. Koenig, R. Guesten, Y. X. Lin, X. W. Zheng, J. Esimbek, J. J. Zhou

Abstract: We aim to directly determine the kinetic temperature and spatial density with formaldehyde for the $\sim$100 brightest ATLASGAL-selected clumps at 870 $μ$m representing various evolutionary stages of high-mass star formation. Ten transitions ($J$ = 3-2 and 4-3) of ortho- and para-H$_2$CO near 211, 218, 225, and 291 GHz were observed with the APEX 12 m telescope. Using non-LTE models with RADEX, we… ▽ More We aim to directly determine the kinetic temperature and spatial density with formaldehyde for the $\sim$100 brightest ATLASGAL-selected clumps at 870 $μ$m representing various evolutionary stages of high-mass star formation. Ten transitions ($J$ = 3-2 and 4-3) of ortho- and para-H$_2$CO near 211, 218, 225, and 291 GHz were observed with the APEX 12 m telescope. Using non-LTE models with RADEX, we derive the gas kinetic temperature and spatial density using the measured p-H$_2$CO 3$_{21}$-2$_{20}$/3$_{03}$-2$_{02}$, 4$_{22}$-3$_{21}$/4$_{04}$-3$_{03}$, and 4$_{04}$-3$_{03}$/3$_{03}$-2$_{02}$ ratios. The gas kinetic temperatures derived from the p-H$_2$CO 3$_{21}$-2$_{20}$/3$_{03}$-2$_{02}$ and 4$_{22}$-3$_{21}$/4$_{04}$-3$_{03}$ line ratios are high, ranging from 43 to $>$300 K with an unweighted average of 91 $\pm$ 4 K. Deduced $T_{\rm kin}$ values from the $J$ = 3-2 and 4-3 transitions are similar. Spatial densities of the gas derived from the p-H$_2$CO 4$_{04}$-3$_{03}$/3$_{03}$-2$_{02}$ line ratios yield 0.6-8.3 $\times$ 10$^6$ cm$^{-3}$ with an unweighted average of 1.5 ($\pm$0.1) $\times$ 10$^6$ cm$^{-3}$. A comparison of kinetic temperatures derived from p-H$_2$CO, NH$_3$, and the dust emission indicates that p-H$_2$CO traces a distinctly higher temperature than the NH$_3$ (2,2)/(1,1) transitions and the dust, tracing heated gas more directly associated with the star formation process. The H$_2$CO linewidths are found to be correlated with bolometric luminosity and increase with the evolutionary stage of the clumps, which suggests that higher luminosities tend to be associated with a more turbulent molecular medium. It seems that the spatial densities measured with H$_2$CO do not vary significantly with the evolutionary stage of the clumps. However, averaged gas kinetic temperatures derived from H$_2$CO increase with time through the evolution of the clumps. △ Less

Submitted 27 November, 2017; originally announced November 2017.

Comments: Accepted for publication in A&A

Journal ref: A&A 611, A6 (2018)

Galactocentric variation of the gas-to-dust ratio and its relation with metallicity

Authors: A. Giannetti, S. Leurini, C. Koenig, J. S. Urquhart, T. Pillai, J. Brand, J. Kauffmann, F. Wyrowski, K. M. Menten

Abstract: (Abridged) Context: The assumption of a gas-to-dust mass ratio (γ) is a common approach to estimate the basic properties of molecular clouds, such as total mass and column density of molecular hydrogen, from (sub)mm continuum observations of the dust. In the Milky Way a single value is used at all galactocentric radii, independently of the observed metallicity gradients. Both models and extragalac… ▽ More (Abridged) Context: The assumption of a gas-to-dust mass ratio (γ) is a common approach to estimate the basic properties of molecular clouds, such as total mass and column density of molecular hydrogen, from (sub)mm continuum observations of the dust. In the Milky Way a single value is used at all galactocentric radii, independently of the observed metallicity gradients. Both models and extragalactic observations suggest that this quantity increases for decreasing metallicity Z, typical of the outer regions in disks, where fewer heavy elements are available to form dust grains. Aims: We aim to investigate the variation of the gas-to-dust ratio as a function of galactocentric radius and metallicity, to allow a more accurate characterisation of the quantity of molecular gas across the galactic disk, as derived from observations of the dust. Methods: Observations of the optically thin C\$^{18}\$O (2-1) transition were obtained with the APEX telescope for a sample of 23 massive and dense star-forming regions in the far outer Galaxy (galactocentric distance greater than 14 kpc). From the modelling of this line and of the spectral energy distribution of the selected clumps we computed the gas-to-dust ratio and compared it to that of well-studied sources from the ATLASGAL TOP100 sample in the inner galactic disk. Results: The gradient in gas-to-dust ratio is found to be 0.087 dex/kpc (or equivalently γ proportional to Z\$^{-1.4}\$). The dust-to-metal ratio, decreases with galactocentric radius, which is the most common situation also for external late-type galaxies. This suggests that grain growth dominates over destruction. The predicted gas-to-dust ratio is in excellent agreement with the estimates in Magellanic clouds, for the appropriate value of Z. △ Less

Submitted 16 October, 2017; originally announced October 2017.

Comments: Accepted for publication on A&A

Journal ref: A&A 606, L12 (2017)

Deuterated methanol on Solar System scale around the HH212 protostar

Authors: E. Bianchi, C. Codella, C. Ceccarelli, V. Taquet, S. Cabrit, F. Bacciotti, R. Bachiller, E. Chapillon, F. Gueth, A. Gusdorf, B. Lefloch, S. Leurini, L. Podio, K. L. J. Rygl, B. Tabone, M. Tafalla

Abstract: Context: Methanol is thought to be mainly formed during the prestellar phase and its deuterated form keeps memory of the conditions at that epoch. Thanks to the unique combination of high angular resolution and sensitivity provided by ALMA, we wish to measure methanol deuteration in the planet formation region around a Class 0 protostar and to understand its origin. Aims: We mapped both the… ▽ More Context: Methanol is thought to be mainly formed during the prestellar phase and its deuterated form keeps memory of the conditions at that epoch. Thanks to the unique combination of high angular resolution and sensitivity provided by ALMA, we wish to measure methanol deuteration in the planet formation region around a Class 0 protostar and to understand its origin. Aims: We mapped both the $^{13}$CH$_3$OH and CH$_2$DOH distribution in the inner regions ($\sim$100 au) of the HH212 system in Orion B. To this end, we used ALMA Cycle 1 and Cycle 4 observations in Band 7 with angular resolution down to $\sim$0.15$"$. Results: We detected 6 lines of $^{13}$CH$_3$OH and 13 lines of CH$_2$DOH with upper level energies up to 438 K in temperature units. We derived a rotational temperature of (171 $\pm$ 52) K and column densities of 7$\times$10$^{16}$ cm$^{-2}$ ($^{13}$CH$_3$OH) and 1$\times$10$^{17}$ cm$^{-2}$ (CH$_2$DOH), respectively. Consequently, the D/H ratio is (2.4 $\pm$ 0.4)$\times$10$^{-2}$, a value lower by an order of magnitude with respect to what was previously measured using single dish telescopes toward protostars located in Perseus. Our findings are consistent with the higher dust temperatures in Orion B with respect to that derived for the Perseus cloud. The emission is tracing a rotating structure extending up to 45 au from the jet axis and elongated by 90 au along the jet axis. So far, the origin of the observed emission appears to be related with the accretion disk. Only higher spatial resolution measurements however, will be able to disentangle between different possible scenarios: disk wind, disk atmosphere, or accretion shocks. △ Less

Submitted 14 September, 2017; originally announced September 2017.

Comments: Astronomy & Astrophysics Letter, in press

Journal ref: A&A 606, L7 (2017)

Fast deuterium fractionation in magnetized and turbulent filaments

Authors: Bastian Körtgen, Stefano Bovino, Dominik R. G. Schleicher, Amelia Stutz, Robi Banerjee, Andrea Giannetti, Silvia Leurini

Abstract: Deuterium fractionation is considered as an important process to infer the chemical ages of prestellar cores in filaments. We present here the first magneto-hydrodynamical simulations including a chemical network to study deuterium fractionation in magnetized and turbulent filaments and their substructures. The filaments typically show widespread deuterium fractionation with average values… ▽ More Deuterium fractionation is considered as an important process to infer the chemical ages of prestellar cores in filaments. We present here the first magneto-hydrodynamical simulations including a chemical network to study deuterium fractionation in magnetized and turbulent filaments and their substructures. The filaments typically show widespread deuterium fractionation with average values $\gtrsim0.01$. For individual cores of similar age, we observe the deuteration fraction to increase with time, but also to be independent of their average properties such as density, virial or mass-to-magnetic flux ratio. We further find a correlation of the deuteration fraction with core mass, average H$_2$ density and virial parameter only at late evolutionary stages of the filament and attribute this to the lifetime of the individual cores. Specifically, chemically old cores reveal higher deuteration fractions. Within the radial profiles of selected cores, we notice differences in the structure of the deuteration fraction or surface density, which we can attribute to their different turbulent properties. High deuteration fractions of the order $0.01-0.1$ may be reached within approximately $200$~kyrs, corresponding to two free-fall times, as defined for cylindrical systems, of the filaments △ Less

Submitted 12 September, 2017; originally announced September 2017.

Comments: submitted to MNRAS. Comments welcome

Testing the variability of the proton-to-electron mass ratio from observations of methanol in the dark cloud core L1498

Authors: M. Daprà, C. Henkel, S. A. Levshakov, K. M. Menten, S. Muller, H. L. Bethlem, S. Leurini, A. V. Lapinov, W. Ubachs

Abstract: The dependence of the proton-to-electron mass ratio, mu, on the local matter density was investigated using methanol emission in the dense dark cloud core L1498. Towards two different positions in L1498, five methanol transitions were detected and an extra line was tentatively detected at a lower confidence level in one of the positions. The observed centroid frequencies were then compared with th… ▽ More The dependence of the proton-to-electron mass ratio, mu, on the local matter density was investigated using methanol emission in the dense dark cloud core L1498. Towards two different positions in L1498, five methanol transitions were detected and an extra line was tentatively detected at a lower confidence level in one of the positions. The observed centroid frequencies were then compared with their rest frame frequencies derived from least-squares fitting to a large data set. Systematic effects, as the underlying methanol hyperfine structure and the Doppler tracking of the telescope, were investigated and their effects were included in the total error budget. The comparison between the observations and the rest frame frequencies constrains potential mu variation at the level of Dmu/mu < 6 x 10^(-8), at a 3 sigma confidence level. For the dark cloud we determine a total CH3OH (A+E) beam averaged column density of 3-4 x 10^(12) cm(-2) (within roughly a factor of two), an E- to A-type methanol column density ratio of N(A-CH3OH)/N(E-CH3OH) = 1.00 +/- 0.15, a density of n(H2) = 3 x 10^5 cm^(-3) (again within a factor of two), and a kinetic temperature of Tkin = 6 +/- 1 K. In a kinetic model including the line intensities observed for the methanol lines, the n(H2) density is higher and the temperature is lower than that derived in previous studies based on different molecular species; the intensity of the 1_0 --> 1_-1 E line strength is not well reproduced. △ Less

Submitted 10 September, 2017; originally announced September 2017.

Comments: 10 pages, 6 figures, accepted for publication in MNRAS

ATLASGAL --- properties of a complete sample of Galactic clumps

Authors: J. S. Urquhart, C. Koenig, A. Giannetti, S. Leurini, T. J. T. Moore, D. J. Eden, T. Pillai, M. A. Thompson, C. Braiding, M. G. Burton, T. Csengeri, J. T. Dempsey, C. Figura, D. Froebrich, K. M. Menten, F. Schuller, M. D. Smith, F. Wyrowski

Abstract: Abridged: ATLASGAL is an unbiased 870 micron submillimetre survey of the inner Galactic plane. It provides a large and systematic inventory of all massive, dense clumps in the Galaxy (>1000 Msun) and includes representative samples of all embedded stages of high-mass star formation. Here we present the first detailed census of the properties (velocities, distances, luminosities and masses) and spa… ▽ More Abridged: ATLASGAL is an unbiased 870 micron submillimetre survey of the inner Galactic plane. It provides a large and systematic inventory of all massive, dense clumps in the Galaxy (>1000 Msun) and includes representative samples of all embedded stages of high-mass star formation. Here we present the first detailed census of the properties (velocities, distances, luminosities and masses) and spatial distribution of a complete sample of ~8000 dense clumps located in the Galactic disk. We derive highly reliable velocities and distances to ~97% of the sample and use mid- and far-infrared survey data to develop an evolutionary classification scheme that we apply to the whole sample. Comparing the evolutionary subsamples reveals trends for increasing dust temperatures, luminosities and line-widths as a function of evolution indicating that the feedback from the embedded proto-clusters is having a significant impact on the structure and dynamics of their natal clumps. We find 88\,per\,cent are already associated with star formation at some level. We also find the clump mass to be independent of evolution suggesting that the clumps form with the majority of their mass in-situ. We estimate the statistical lifetime of the quiescent stage to be ~5 x 10^4 yr for clump masses ~1000 Msun decreasing to ~1 x 10^4 yr for clump masses >10000 Msun. We find a strong correlation between the fraction of clumps associated with massive stars and peak column density. The fraction is initially small at low column densities but reaching 100\,per\,cent for column densities above 10^{23} cm^{-2}; there are no clumps with column density clumps above this value that are not already associated with massive star formation. All of the evidence is consistent with a dynamic view of star formation wherein the clumps form rapidly and are initially very unstable so that star formation quickly ensues. △ Less

Submitted 5 September, 2017; v1 submitted 1 September, 2017; originally announced September 2017.

Comments: Accepted for publication in MNRAS. Consists of 45 pages, including 10 pages of appendices. Full versions of Tables 2, 5, 8 and A2 will be made available via CDS once the paper has been published as will full versions of Figures 6, 8 and A1. The quality of the images has been reduced but a high-resolution version is available on request. V2: some minor corrections made