Showing 1–50 of 311 results for author: Beuther, H

A global view on star formation: The GLOSTAR Galactic plane survey X. Galactic HII region catalog using radio recombination lines

Authors: S. Khan, M. R. Rugel, A. Brunthaler, K. M. Menten, F. Wyrowski, J. S. Urquhart, Y. Gong, A. Y. Yang, H. Nguyen, R. Dokara, S. A. Dzib, S. -N. X. Medina, G. N. Ortiz-León, J. D. Pandian, H. Beuther, V. S. Veena, S. Neupane, A. Cheema, W. Reich, N. Roy

Abstract: Studies of Galactic HII regions are of crucial importance for studying star formation and the evolution of the interstellar medium. Gaining an insight into their physical characteristics contributes to a more comprehensive understanding of these phenomena. The GLOSTAR project aims to provide a GLObal view on STAR formation in the Milky Way by performing an unbiased and sensitive survey. This is ac… ▽ More Studies of Galactic HII regions are of crucial importance for studying star formation and the evolution of the interstellar medium. Gaining an insight into their physical characteristics contributes to a more comprehensive understanding of these phenomena. The GLOSTAR project aims to provide a GLObal view on STAR formation in the Milky Way by performing an unbiased and sensitive survey. This is achieved by using the extremely wideband (4{-}8 GHz) C-band receiver of the Karl G. Jansky Very Large Array and the Effelsberg 100 m telescope. Using radio recombination lines observed in the GLOSTAR survey with the VLA in D-configuration with a typical line sensitivity of 1σしぐま {\sim} 3.0 mJy beam{^-1} at {\sim} 5 km s{^-1} and an angular resolution of 25", we cataloged 244 individual Galactic HII regions and derived their physical properties. We examined the mid-infrared (MIR) morphology of these HII regions and find that a significant portion of them exhibit a bubble-like morphology in the GLIMPSE 8 μみゅーm emission. We also searched for associations with the dust continuum and sources of methanol maser emission, other tracers of young stellar objects, and find that 48\% and 14\% of our HII regions, respectively, are coextensive with those. We measured the electron temperature for a large sample of HII regions within Galactocentric distances spanning from 1.6 to 13.1 kpc and derived the Galactic electron temperature gradient as {\sim} 372 {\pm} 28 K kpc{^-1} with an intercept of 4248 {\pm} 161 K, which is consistent with previous studies. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Accepted for publication in A&A

MagMar III -- Resisting the Pressure, Is the Magnetic Field Overwhelmed in NGC6334I?

Authors: Paulo C. Cortes, Josep M. Girart, Patricio Sanhueza, Junhao Liu, Sergio Martin, Ian W. Stephens, Henrik Beuther, Patrick M. Koch, M. Fernandez-Lopez, Alvaro Sanchez-Monge, Jia-Wei Wang, Kaho Morii, Shanghuo Li, Piyali Saha, Qizhou Zhang, David Rebolledo, Luis A. Zapata, Ji-hyun Kang, Wenyu Jiao, Jongsoo Kim, Yu Cheng, Jihye Hwang, Eun Jung Chung, Spandan Choudhury, A-Ran Lyo , et al. (1 additional authors not shown)

Abstract: We report on ALMA observations of polarized dust emission at 1.2 mm from NGC6334I, a source known for its significant flux outbursts. Between five months, our data show no substantial change in total intensity and a modest 8\% variation in linear polarization, suggesting a phase of stability or the conclusion of the outburst. The magnetic field, inferred from this polarized emission, displays a pr… ▽ More We report on ALMA observations of polarized dust emission at 1.2 mm from NGC6334I, a source known for its significant flux outbursts. Between five months, our data show no substantial change in total intensity and a modest 8\% variation in linear polarization, suggesting a phase of stability or the conclusion of the outburst. The magnetic field, inferred from this polarized emission, displays a predominantly radial pattern from North-West to South-East with intricate disturbances across major cores, hinting at spiral structures. Energy analysis of CS$(J=5 \rightarrow 4)$ emission yields an outflow energy of approximately $3.5\times10^{45}$ ergs, aligning with previous interferometric studies. Utilizing the Davis-Chandrasekhar-Fermi method, we determined magnetic field strengths ranging from 1 to 11 mG, averaging at 1.9 mG. This average increases to 4 $\pm 1$ mG when incorporating Zeeman measurements. Comparative analyses using gravitational, thermal, and kinetic energy maps reveal that magnetic energy is significantly weaker, possibly explaining the observed field morphology. We also find that the energy in the outflows and the expanding cometary {\HII} region is also larger than the magnetic energy, suggesting that protostellar feedback maybe the dominant driver behind the injection of turbulence in NGC6334I at the scales sampled by our data. The gas in NGC6334I predominantly exhibits supersonic and trans-Alfvenic conditions, transitioning towards a super-Alfvenic regime, underscoring a diminished influence of the magnetic field with increasing gas density. These observations are in agreement with prior polarization studies at 220 GHz, enriching our understanding of the dynamic processes in high-mass star-forming regions. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: Accepted for Publication at the Astrophysical Journal

The Cygnus Allscale Survey of Chemistry and Dynamical Environments: CASCADE III. The large scale distribution of DCO+, DNC and DCN in the DR21 filament

Authors: I. Barlach Christensen, F. Wyrowski, V. S. Veena, H. Beuther, D. Semenov, K. M. Menten, A. M. Jacob, W. -J. Kim, N. Cunningham, C. Gieser, A. Hacar, S. Li, N. Schneider, I. Skretas, J. M. Winters

Abstract: Deuterated molecules and their molecular D/H-ratios (RD(D)) are important diagnostic tools to study the physical conditions of star-forming regions. The degree of deuteration, RD(D), can be significantly enhanced over the elemental D/H-ratio depending on physical parameters. Within the Cygnus Allscale Survey of Chemistry and Dynamical Environments (CASCADE), we aim to explore the large-scale distr… ▽ More Deuterated molecules and their molecular D/H-ratios (RD(D)) are important diagnostic tools to study the physical conditions of star-forming regions. The degree of deuteration, RD(D), can be significantly enhanced over the elemental D/H-ratio depending on physical parameters. Within the Cygnus Allscale Survey of Chemistry and Dynamical Environments (CASCADE), we aim to explore the large-scale distribution of deuterated molecules in the nearby Cygnus-X region. We focus on the analysis of large-scale structures of deuterated molecules in the filamentary region hosting the prominent Hii region DR21 and DR21(OH). Here we discuss the HCO+, HNC and HCN molecules and their deuterated isotopologues DCO+, DNC and DCN. The spatial distributions of integrated line emissions from DCO+, DNC, and DCN reveal morphological differences. DCO+ displays the most extended emission, characterized by several prominent peaks. Likewise, DNC exhibits multiple peaks, although its emission appears less extended compared to DCO+. In contrast to the extended emission of DCO+ and DNC, DCN appears the least extended, with distinct peaks. Focusing only on the regions where all three molecules are observed, the mean deuteration ratios for each species are 0.01 for both DNC and DCN, and = 0.005 for DCO+. Anti-correlations are found with deuterated molecules and dust temperature or N(H2). The strongest anti-correlation is found with RD(DCO+) and N(H2). The anti-correlation of RD(DCO+) and N(H2) is suggested to be a result of a combination of an increased photodissociation degree and shocks. A strong positive correlation between the ratio of integrated intensities of DCN and DNC with their 13C-isotopologues, are found in high column density regions. The positive relationship between the ratios implies that the D-isotopologue of the isomers could potentially serve as a tracer for the kinetic gas temperature. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 24 pages, 21 figures, accepted to A&A

Spectrum and polarization of the Galactic center radio transient ASKAP J173608.2-321635 from THOR-GC and VLITE

Authors: Kierra J. Weatherhead, Jeroen M. Stil, Michael Rugel, Wendy M. Peters, Loren Anderson, Ashley Barnes, Henrik Beuther, Tracy E. Clarke, Sergio A. Dzib, Paul Goldsmith, Karl M. Menten, Kristina E. Nyland, Mattia C. Sormani, James Urquhart

Abstract: The radio transient ASKAP J173608.2-321735, at the position (l,b)= (356.0872,-0.0390), was serendipitously observed by The HI/OH/Recombination Line Survey of the Galactic Center (THOR-GC) at three epochs in March 2020, April 2020 and February 2021. The source was detected only on 2020 April 11 with flux density 20.6 +/- 1.1 mJy at 1.23 GHz and in-band spectral index alpha = -3.1 +/- 0.2. The comme… ▽ More The radio transient ASKAP J173608.2-321735, at the position (l,b)= (356.0872,-0.0390), was serendipitously observed by The HI/OH/Recombination Line Survey of the Galactic Center (THOR-GC) at three epochs in March 2020, April 2020 and February 2021. The source was detected only on 2020 April 11 with flux density 20.6 +/- 1.1 mJy at 1.23 GHz and in-band spectral index alpha = -3.1 +/- 0.2. The commensal VLA Low-band Ionsophere and Transient Experiment (VLITE) simultaneously detected the source at 339 MHz with a flux density 122.6 +/- 20.4 mJy, indicating a spectral break below 1 GHz. The rotation measure in April 2020 was 63.9 +/- 0.3rad/m2, which almost triples the range of the variable rotation measure observed by Wang et al. (2021) to ~130 rad/m2. The polarization angle, corrected for Faraday rotation, was 97 +/- 6 degrees. The 1.23 GHz linear polarization was 76.7% +/- 3.9% with wavelength-dependent depolarization indicating Faraday depth dispersion sigma_phi = 4.8^{+0.5}_{-0.7} rad/m2. We find an upper limit to circular polarization |V|/I < 10.1%. Interpretation of the data in terms of diffractive scattering of radio waves by a plasma near the source indicates electron density and line-of-sight magnetic field strength within a factor 3 of n_e ~2 cm^{-3} and B_par ~2 x 10^5 microgauss. Combined with causality limits to the size of the source, these parameters are consistent with the low-frequency spectral break resulting from synchrotron self-absorption, not free-free absorption. A possible interpretation of the source is a highly supersonic neutron star interacting with a changing environment. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: Accepted for publication in the Astrophysical Journal

Magnetic field dragging in filamentary molecular clouds

Authors: Domitilla Tapinassi, Daniele Galli, Marco Padovani, Henrik Beuther

Abstract: Maps of polarized dust emission of molecular clouds reveal the morphology of the magnetic field associated with star-forming regions. In particular, polarization maps of hub-filament systems show the distortion of magnetic field lines induced by gas flows onto and inside filaments. We aim to understand the relation between the curvature of magnetic field lines associated with filaments in hub-fila… ▽ More Maps of polarized dust emission of molecular clouds reveal the morphology of the magnetic field associated with star-forming regions. In particular, polarization maps of hub-filament systems show the distortion of magnetic field lines induced by gas flows onto and inside filaments. We aim to understand the relation between the curvature of magnetic field lines associated with filaments in hub-filament systems and the properties of the underlying gas flows. We consider steady-state models of gas with finite electrical resistivity flowing across a transverse magnetic field. We derive the relation between the bending of the field lines and the flow parameters represented by the Alfvén Mach number and the magnetic Reynolds number. We find that, on the scale of the filaments, the relevant parameter for a gas of finite electrical resistivity is the magnetic Reynolds number, and we derive the relation between the deflection angle of the field from the initial direction (assumed perpendicular to the filament) and the value of the electrical resistivity, due to either Ohmic dissipation or ambipolar diffusion. Application of this model to specific observations of polarized dust emission in filamentary clouds shows that magnetic Reynolds numbers of a few tens are required to reproduce the data. Despite significant uncertainties in the observations (the flow speed, the geometry and orientation of the filament), and the idealization of the model, the specific cases considered show that ambipolar diffusion can provide the resistivity needed to maintain a steady state flow across magnetic fields of significant strength over realistic time scales. △ Less

Submitted 19 May, 2024; originally announced May 2024.

Comments: 9 pages, 4 figures. Accepted for publication in Astronomy & Astrophysics

3D structure of the Milky Way out to 10 kpc from the Sun. Catalogue of large molecular clouds in the Galactic Plane

Authors: Sara Rezaei Kh., Henrik Beuther, Robert A. Benjamin, Anna-Christina Eilers, Thomas Henning, Maria J. Jiménez-Donaire, Marc-Antoine Miville-Deschênes

Abstract: Understanding the 3D structure of the Milky Way is a crucial step in deriving properties of the star-forming regions, as well as the Galaxy as a whole. We present a novel 3D map of the Milky Way plane that extends to 10 kpc distance from the Sun. We leverage the wealth of information in the near-IR APOGEE dataset and combine that with our state-of-the-art 3D mapping technique using Bayesian statis… ▽ More Understanding the 3D structure of the Milky Way is a crucial step in deriving properties of the star-forming regions, as well as the Galaxy as a whole. We present a novel 3D map of the Milky Way plane that extends to 10 kpc distance from the Sun. We leverage the wealth of information in the near-IR APOGEE dataset and combine that with our state-of-the-art 3D mapping technique using Bayesian statistics and the Gaussian process to provide a large-scale 3D map of the dust in the Milky Way. Our map stretches across 10 kpc along both the X and Y axes, and 750 pc in the Z direction, perpendicular to the Galactic plane. Our results reveal multi-scale over-densities as well as large cavities in the Galactic plane and shed new light on the Galactic structure and spiral arms. We also provide a catalogue of large molecular clouds identified by our map with accurate distance and volume density estimates. Utilising volume densities derived from this map, we explore mass distribution across various Galactocentric radii. A general decline towards the outer Galaxy is observed, followed by local peaks, some aligning with established features like the Molecular Ring and segments of the spiral arms. Moreover, this work explores extragalactic observational effects on derived properties of molecular clouds by demonstrating the potential biases arising from column density measurements in inferring properties of these regions, and opens exciting avenues for further exploration and analysis, offering a deeper perspective on the complex processes that shape our galaxy and beyond. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: Revised version after first referee report; to appear in A&A. Includes SRKh's goodbye letter to Astronomy!

JWST observations of the Horsehead photon-dominated region I. First results from multi-band near- and mid-infrared imaging

Authors: A. Abergel, K. Misselt, K. D. Gordon, A. Noriega-Crespo, P. Guillard, D. Van De Putte, A. N. Witt, N. Ysard, M. Baes, H. Beuther, P. Bouchet, B. R. Brandl, M. Elyajouri, O. Kannavou, S. Kendrew, P. Klassen, B. Trahin

Abstract: The JWST has captured the sharpest IR images ever taken of the Horsehead nebula, a prototypical moderately irradiated PDR that is fully representative of most of the UV-illuminated molecular gas in the Milky Way and star-forming galaxies. We investigate the impact of FUV radiation of a molecular cloud and constrain the structure of the edge of the PDR and its illumination conditions. We used NIRCa… ▽ More The JWST has captured the sharpest IR images ever taken of the Horsehead nebula, a prototypical moderately irradiated PDR that is fully representative of most of the UV-illuminated molecular gas in the Milky Way and star-forming galaxies. We investigate the impact of FUV radiation of a molecular cloud and constrain the structure of the edge of the PDR and its illumination conditions. We used NIRCam and MIRI to obtain 17 broadband and 6 narrowband maps from 0.7 to 28 $μみゅー$m. We mapped the dust emission, scattered light, and several gas phase lines. We also used HST-WFC3 maps at 1.1 and 1. 6 $μみゅー$m, along with HST-STIS spectroscopic observations of the H$αあるふぁ$ line. We probed the structure of the edge of the Horsehead and resolved its spatial complexity. We detected a network of faint striated features extending perpendicularly to the PDR front into the H\,II region in filters sensitive to nano-grain emission and light scattered by larger grains. This may indeed figure as the first detection of the entrainment of dust particles in the evaporative flow. The map of the 1-0 S(1) line of H$_2$ presents sharp sub-structures on scales as small as 1.5 arcsec. The ionization and dissociation fronts appear at distances 1-2 arcsec behind the edge of the PDR and seem to spatially coincide, indicating a thickness of the neutral atomic layer below 100 au. All broadband maps present strong color variations which can be explained by dust attenuation. Deviations of the emissions in the H$αあるふぁ$, Pa$αあるふぁ,$ and Br$αあるふぁ$ lines also indicate dust attenuation. With a very simple model, we derive the main features of the extinction curve. A small excess of extinction at 3 $μみゅー$m may be attributed to icy H$_2$O mantles onto grains. In all lines of sight crossing the inner regions of the Horsehead, it appears that dust attenuation is non-negligible over the entire spectral range. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 25 pages, 17 figures

JWST detections of amorphous and crystalline HDO ice toward massive protostars

Authors: Katerina Slavicinska, Ewine F. van Dishoeck, Łukasz Tychoniec, Pooneh Nazari, Adam E. Rubinstein, Robert Gutermuth, Himanshu Tyagi, Yuan Chen, Nashanty G. C. Brunken, Will R. M. Rocha, P. Manoj, Mayank Narang, S. Thomas Megeath, Yao-Lun Yang, Leslie W. Looney, John J. Tobin, Henrik Beuther, Tyler L. Bourke, Harold Linnartz, Samuel Federman, Dan M. Watson, Hendrik Linz

Abstract: This work aims to utilize the increased sensitivity and resolution of the JWST to quantify the HDO/H$_{2}$O ratio in ices toward young stellar objects (YSOs) and to determine if the HDO/H$_{2}$O ratios measured in the gas phase toward massive YSOs (MYSOs) are representative of the ratios in their ice envelopes. Two protostars observed in the Investigating Protostellar Accretion (IPA) program using… ▽ More This work aims to utilize the increased sensitivity and resolution of the JWST to quantify the HDO/H$_{2}$O ratio in ices toward young stellar objects (YSOs) and to determine if the HDO/H$_{2}$O ratios measured in the gas phase toward massive YSOs (MYSOs) are representative of the ratios in their ice envelopes. Two protostars observed in the Investigating Protostellar Accretion (IPA) program using JWST NIRSpec were analyzed: HOPS 370, an intermediate-mass YSO (IMYSO), and IRAS 20126+4104, a MYSO. The HDO ice toward these sources was detected above the 3$σしぐま$ level and quantified via its 4.1 $μみゅー$m band. The contributions from the CH$_{3}$OH combination modes to the observed optical depth in this spectral region were constrained via the CH$_{3}$OH 3.53 $μみゅー$m band to ensure that the integrated optical depth of the HDO feature was not overestimated. H$_{2}$O ice was quantified via its 3 $μみゅー$m band. From these fits, ice HDO/H$_{2}$O abundance ratios of 4.6$\pm$1.8$\times$10$^{-3}$ and 2.6$\pm$1.2$\times$10$^{-3}$ are obtained for HOPS 370 and IRAS 20126+4104, respectively. The simultaneous detections of both crystalline HDO and crystalline H$_{2}$O corroborate the assignment of the observed feature at 4.1 $μみゅー$m to HDO ice. The ice HDO/H$_{2}$O ratios are similar to the highest reported gas HDO/H$_{2}$O ratios measured toward MYSOs as well as the hot inner regions of isolated low-mass protostars, suggesting that at least some of the gas HDO/H$_{2}$O ratios measured toward massive hot cores are representative of the HDO/H$_{2}$O ratios in ices. The need for an H$_{2}$O-rich CH$_{3}$OH component in the CH$_{3}$OH ice analysis supports recent experimental and observational results that indicate that some CH$_{3}$OH ice may form prior to the CO freeze-out stage in H$_{2}$O-rich ice layers. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: Accepted for publication in A&A. 23 pages, 17 figures, 10 tables

JWST/MIRI detection of suprathermal OH rotational emissions: probing the dissociation of the water by Lyman alpha photons near the protostar HOPS 370

Authors: David A. Neufeld, P. Manoj, Himanshu Tyagi, Mayank Narang, Dan M. Watson, S. Thomas Megeath, Ewine F. Van Dishoeck, Robert A. Gutermuth, Thomas Stanke, Yao-Lun Yang, Adam E. Rubinstein, Guillem Anglada, Henrik Beuther, Alessio Caratti o Garatti, Neal J. Evans II, Samuel Federman, William J. Fischer, Joel Green, Pamela Klaassen, Leslie W. Looney, Mayra Osorio, Pooneh Nazari, John J. Tobin, Lukasz Tychoniec, Scott Wolk

Abstract: Using the MIRI/MRS spectrometer on JWST, we have detected pure rotational, suprathermal OH emissions from the vicinity of the intermediate-mass protostar HOPS 370 (OMC2/FIR3). These emissions are observed from shocked knots in a jet/outflow, and originate in states of rotational quantum number as high as 46 that possess excitation energies as large as $E_U/k = 4.65 \times 10^4$ K. The relative str… ▽ More Using the MIRI/MRS spectrometer on JWST, we have detected pure rotational, suprathermal OH emissions from the vicinity of the intermediate-mass protostar HOPS 370 (OMC2/FIR3). These emissions are observed from shocked knots in a jet/outflow, and originate in states of rotational quantum number as high as 46 that possess excitation energies as large as $E_U/k = 4.65 \times 10^4$ K. The relative strengths of the observed OH lines provide a powerful diagnostic of the ultraviolet radiation field in a heavily-extinguished region ($A_V \sim 10 - 20$) where direct UV observations are impossible. To high precision, the OH line strengths are consistent with a picture in which the suprathermal OH states are populated following the photodissociation of water in its $\tilde B - X$ band by ultraviolet radiation produced by fast ($\sim 80\,\rm km\,s^{-1}$) shocks along the jet. The observed dominance of emission from symmetric ($A^\prime$) OH states over that from antisymmetric ($A^{\prime\prime}$) states provides a distinctive signature of this particular population mechanism. Moreover, the variation of intensity with rotational quantum number suggests specifically that Ly$αあるふぁ$ radiation is responsible for the photodissociation of water, an alternative model with photodissociation by a 10$^4$ K blackbody being disfavored at a high level of significance. Using measurements of the Br$αあるふぁ$ flux to estimate the Ly$αあるふぁ$ production rate, we find that $\sim 4\%$ of the Ly$αあるふぁ$ photons are absorbed by water. Combined with direct measurements of water emissions in the $νにゅー_2 = 1 -0$ band, the OH observations promise to provide key constraints on future models for the diffusion of Ly$αあるふぁ$ photons in the vicinity of a shock front. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 30 pages, 7 figures. Accepted for publication in ApJ Letters

The ALMA Survey of 70 $μみゅー$m Dark High-mass Clumps in Early Stages (ASHES). XI. Statistical Study of Early Fragmentation

Authors: Kaho Morii, Patricio Sanhueza, Qizhou Zhang, Fumitaka Nakamura, Shanghuo Li, Giovanni Sabatini, Fernando A. Olguin, Henrik Beuther, Daniel Tafoya, Natsuko Izumi, Ken'ichi Tatematsu, Takeshi Sakai

Abstract: Fragmentation during the early stages of high-mass star formation is crucial for understanding the formation of high-mass clusters. We investigated fragmentation within thirty-nine high-mass star-forming clumps as part of the Atacama Large Millimeter/submillimeter Array (ALMA) Survey of 70 $μみゅー$m Dark High-mass Clumps in Early Stages (ASHES). Considering projection effects, we have estimated core se… ▽ More Fragmentation during the early stages of high-mass star formation is crucial for understanding the formation of high-mass clusters. We investigated fragmentation within thirty-nine high-mass star-forming clumps as part of the Atacama Large Millimeter/submillimeter Array (ALMA) Survey of 70 $μみゅー$m Dark High-mass Clumps in Early Stages (ASHES). Considering projection effects, we have estimated core separations for 839 cores identified from the continuum emission and found mean values between 0.08 and 0.32 pc within each clump. We find compatibility of the observed core separations and masses with the thermal Jeans length and mass, respectively. We also present sub-clump structures revealed by the 7 m-array continuum emission. Comparison of the Jeans parameters using clump and sub-clump densities with the separation and masses of gravitationally bound cores suggests that they can be explained by clump fragmentation, implying the simultaneous formation of sub-clumps and cores within rather than a step-by-step hierarchical fragmentation. The number of cores in each clump positively correlates with the clump surface density and the number expected from the thermal Jeans fragmentation. We also find that the higher the fraction of protostellar cores, the larger the dynamic range of the core mass, implying that the cores are growing in mass as the clump evolves. The ASHES sample exhibits various fragmentation patterns: aligned, scattered, clustered, and sub-clustered. Using the Q-parameter, which can help to distinguish between centrally condensed and subclustered spatial core distributions, we finally find that in the early evolutionary stages of high-mass star formation, cores tend to follow a subclustered distribution. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: Accepted for Publication in ApJ. 19 pages, 7 figures, 4 tables

Relative alignment between gas structures and magnetic field in Orion A at different scales using different molecular gas tracers

Authors: Wenyu Jiao, Ke Wang, Fengwei Xu, Chao Wang, Henrik Beuther

Abstract: Context: Magnetic fields can play crucial roles in high-mass star formation. Nonetheless, the significance of magnetic fields at various scales and their relationship with gas structures is largely overlooked. Aims: Our goal is to examine the relationship between the magnetic field and molecular gas structures within the Orion A giant molecular cloud at different scales and density regimes. Method… ▽ More Context: Magnetic fields can play crucial roles in high-mass star formation. Nonetheless, the significance of magnetic fields at various scales and their relationship with gas structures is largely overlooked. Aims: Our goal is to examine the relationship between the magnetic field and molecular gas structures within the Orion A giant molecular cloud at different scales and density regimes. Methods: We assess the gas intensity structures and column densities in Orion A by utilizing $^{12}$CO, $^{13}$CO, and C$^{18}$O from Nobeyama observations. Through comparing Nobeyama observations with {\it{Planck}} polarization observations on large scales ($\sim0.6$ pc) and JCMT polarization observations on small scales ($\sim0.04$ pc), we investigate how the role of magnetic fields change with scale and density. Results: We find a similar trend from parallel to perpendicular alignment with increasing column densities in Orion A at both large and small spatial scales. Besides, when changing from low-density to high-density tracers, the relative orientation preference changes from random to perpendicular. The self-similar results at different scales indicate that magnetic fields are dynamically important in both cloud formation and filament formation. However, magnetic fields properties at small scales are relative complicated, and the interplay between magnetic field and star-forming activities needs to be discussed case-by-case. △ Less

Submitted 19 June, 2024; v1 submitted 7 March, 2024; originally announced March 2024.

Comments: 12 pages, 8 figures, published in A&A

JWST Observations of Young protoStars (JOYS): Linked accretion and ejection in a Class I protobinary system

Authors: Łukasz Tychoniec, Martijn L. van Gelder, Ewine F. van Dishoeck, Logan Francis, Will R. M. Rocha, Alessio Caratti o Garatti, Henrik Beuther, Caroline Gieser, Kay Justtanont, Harold Linnartz, Valentin J. M. Le Gouellec, Giulia Perotti, R. Devaraj, Benoît Tabone, Thomas P. Ray, Nashanty G. C. Brunken, Yuan Chen, Patrick J. Kavanagh, Pamela Klaassen, Katerina Slavicinska, Manuel Güdel, Goran Östlin

Abstract: Accretion and ejection sets the outcome of the star and planet formation process. The mid-infrared wavelength range offers key tracers of those processes that were difficult to detect and spatially resolve in protostars until now. We aim to characterize the interplay between accretion and ejection in the low-mass Class I protobinary system TMC1, comprising two young stellar objects: TMC1-W and TMC… ▽ More Accretion and ejection sets the outcome of the star and planet formation process. The mid-infrared wavelength range offers key tracers of those processes that were difficult to detect and spatially resolve in protostars until now. We aim to characterize the interplay between accretion and ejection in the low-mass Class I protobinary system TMC1, comprising two young stellar objects: TMC1-W and TMC1-E with 85 au separation. With the {\it James Webb} Space Telescope (JWST) - Mid-Infrared Instrument (MIRI) observations in 5 - 28 $μみゅー$m range, we measure intensities of emission lines of H$_2$, atoms and ions, e.g., [Fe II] and [Ne II], and HI recombination lines. We detect H$_2$ outflow coming from TMC1-E, with no significant H$_2$ emission from TMC1-W. The H$_2$ emission from TMC1-E outflow appears narrow and extends to wider opening angles with decreasing E$_{up}$ from S(8) to S(1) rotational transitions, indicating a disk wind origin. The outflow from TMC1-E protostar shows spatially extended emission lines of [Ne II], [Ne III], [Ar II], and [Ar III], with their line ratios consistent with UV radiation as a source of ionization. With ALMA, we detect accretion streamer infalling from $>$ 1000 au scales onto the TMC1-E component. TMC1-W protostar powers a collimated jet, detected with [Fe II] and [Ni II] consistent with energetic flow. A much weaker ionized jet is observed from TMC1-E. TMC1-W is associated with strong emission from hydrogen recombination lines, tracing the accretion onto the young star. Observations of a binary Class I protostellar system show that the two processes are clearly intertwined, with accretion from the envelope onto the disk influencing a wide-angle wind ejected on disk scales, while accretion from the protostellar disk onto the protostar is associated with the source launching a collimated high-velocity jet within the innermost regions of the disk. △ Less

Submitted 4 June, 2024; v1 submitted 6 February, 2024; originally announced February 2024.

Comments: 18 pages, 13 figures, accepted for publication in Astronomy and Astrophysics

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

JWST observations of $^{13}$CO$_{2}$ ice: Tracing the chemical environment and thermal history of ices in protostellar envelopes

Authors: Nashanty G. C. Brunken, Will R. M. Rocha, Ewine F. van Dishoeck, Robert Gutermuth, Himanshu Tyagi, Katerina Slavicinska, Pooneh Nazari, S. Thomas Megeath, Neal J. Evans II, Mayank Narang, P. Manoj, Adam E. Rubinstein, Dan M. Watson, Leslie W. Looney, Harold Linnartz, Alessio Caratti o Garatti, Henrik Beuther, Hendrik Linz, Pamela Klaassen, Charles A. Poteet, Samuel Federman, Guillem Anglada, Prabhani Atnagulov, Tyler L. Bourke, William J. Fischer , et al. (16 additional authors not shown)

Abstract: The structure and composition of simple ices can be modified during stellar evolution by protostellar heating. Key to understanding the involved processes are thermal and chemical tracers that can diagnose the history and environment of the ice. The 15.2 $μみゅー$m bending mode of $^{12}$CO$_2$ has proven to be a valuable tracer of ice heating events but suffers from grain shape and size effects. A viab… ▽ More The structure and composition of simple ices can be modified during stellar evolution by protostellar heating. Key to understanding the involved processes are thermal and chemical tracers that can diagnose the history and environment of the ice. The 15.2 $μみゅー$m bending mode of $^{12}$CO$_2$ has proven to be a valuable tracer of ice heating events but suffers from grain shape and size effects. A viable alternative tracer is the weaker $^{13}$CO$_2$ isotopologue band at 4.39 $μみゅー$m which has now become accessible at high S/N with the $\textit{James Webb}$ Space Telescope (JWST). We present JWST NIRSpec observations of $^{13}$CO$_2$ ice in five deeply embedded Class 0 sources spanning a wide range in luminosities (0.2 - 10$^4$ L$_{\odot}$ ) taken as part of the Investigating Protostellar Accretion Across the Mass Spectrum (IPA) program. The band profiles vary significantly, with the most luminous sources showing a distinct narrow peak at 4.38 $μみゅー$m. We first apply a phenomenological approach and show that a minimum of 3-4 Gaussian profiles are needed to fit the $^{13}$CO$_2$ absorption feature. We then combine these findings with laboratory data and show that a 15.2 $μみゅー$m $^{12}$CO$_2$ band inspired five-component decomposition can be applied for the isotopologue band where each component is representative of CO$_2$ ice in a specific molecular environment. The final solution consists of cold mixtures of CO$_2$ with CH$_3$OH, H$_2$O and CO as well as segregated heated pure CO$_2$ ice. Our results are in agreement with previous studies of the $^{12}$CO$_2$ ice band, further confirming that $^{13}$CO$_{2}$ is a useful alternative tracer of protostellar heating events. We also propose an alternative solution consisting only of heated CO$_2$:CH$_3$OH and CO$_2$:H$_2$O ices and warm pure CO$_2$ ice for decomposing the ice profiles of the two most luminous sources in our sample. △ Less

Submitted 7 March, 2024; v1 submitted 6 February, 2024; originally announced February 2024.

Journal ref: A&A 685, A27 (2024)

Hunt for complex cyanides in protostellar ices with JWST: Tentative detection of CH$_3$CN and C$_2$H$_5$CN

Authors: P. Nazari, W. R. M. Rocha, A. E. Rubinstein, K. Slavicinska, M. G. Rachid, E. F. van Dishoeck, S. T. Megeath, R. Gutermuth, H. Tyagi, N. Brunken, M. Narang, P. Manoj, D. M. Watson, N. J. Evans II, S. Federman, J. Muzerolle Page, G. Anglada, H. Beuther, P. Klaassen, L. W. Looney, M. Osorio, T. Stanke, Y. -L. Yang

Abstract: Nitrogen-bearing complex organic molecules have been commonly detected in the gas phase but not yet in interstellar ices. This has led to the long-standing question of whether these molecules form in the gas phase or in ices. $\textit{James Webb}$ Space Telescope ($\textit{JWST}$) offers the sensitivity, spectral resolution, and wavelength coverage needed to detect them in ices and investigate whe… ▽ More Nitrogen-bearing complex organic molecules have been commonly detected in the gas phase but not yet in interstellar ices. This has led to the long-standing question of whether these molecules form in the gas phase or in ices. $\textit{James Webb}$ Space Telescope ($\textit{JWST}$) offers the sensitivity, spectral resolution, and wavelength coverage needed to detect them in ices and investigate whether their abundance ratios are similar in gas and ice. We report the first tentative detection of CH$_3$CN, C$_2$H$_5$CN, and the simple molecule, N$_2$O, based on the CN-stretch band in interstellar ices toward three (HOPS 153, HOPS 370, and IRAS 20126+4104) out of the five protostellar systems observed as part of the Investigating Protostellar Accretion (IPA) GO program with $\textit{JWST}$-NIRSpec. We also provide upper limits for the two other sources with smaller luminosities in the sample. We detect OCN$^-$ in the ices of all sources with typical CH$_3$CN/OCN$^-$ ratios of around 1. Ice and gas column density ratios of the nitrogen-bearing species with respect to each other are better matched than those with respect to methanol, which are a factor of ${\sim}5$ larger in the ices than the gas. We attribute the elevated ice column densities with respect to methanol to the difference in snowline locations of nitrogen-bearing molecules and of methanol, biasing the gas-phase observations toward fewer nitrogen-bearing molecules. Moreover, we find tentative evidence for enhancement of OCN$^-$, CH$_3$CN, and C$_2$H$_5$CN in warmer ices, although formation of these molecules likely starts along with methanol in the cold prestellar phase. Future surveys combining NIRSpec and MIRI, and additional laboratory spectroscopic measurements of C$_2$H$_5$CN ice, are necessary for robust detection and conclusions on the formation history of complex cyanides. △ Less

Submitted 15 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A

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

JOYS: MIRI/MRS spectroscopy of gas-phase molecules from the high-mass star-forming region IRAS 23385+6053

Authors: L. Francis, M. L. van Gelder, E. F. van Dishoeck, C. Gieser, H. Beuther, L. Tychoniec, G. Perotti, A. Caratti o Garatti, P. J. Kavanagh, T. Ray, P. Klaassen, K. Justtanont, H. Linnartz, W. R. M. Rocha, K. Slavicinska, M. Güdel, T. Henning, P. O. Lagage, G. Östlin

Abstract: Space-based mid-IR spectroscopy provides tracers of warm gas in star-forming regions that are inaccessible from the ground. Past mid-IR spectra of bright high-mass protostars in the hot-core phase typically showed strong absorption features from molecules such as CO$_2$, C$_2$H$_2$, and HCN. However, little is known about their fainter counterparts at earlier stages. We thus aim to characterize th… ▽ More Space-based mid-IR spectroscopy provides tracers of warm gas in star-forming regions that are inaccessible from the ground. Past mid-IR spectra of bright high-mass protostars in the hot-core phase typically showed strong absorption features from molecules such as CO$_2$, C$_2$H$_2$, and HCN. However, little is known about their fainter counterparts at earlier stages. We thus aim to characterize the gas-phase molecular features in JWST MIRI/MRS observations of the young high-mass star-forming region IRAS 23385+6053. Spectra were extracted from two mid-IR sources and three H$_2$ bright outflow knots in the MIRI/MRS field of view. Rich molecular spectra with emission from CO, H$_2$, HD, H$_2$O, C$_2$H$_2$, HCN, CO$_2$, and OH are detected towards the two mid-IR sources. However, only CO and OH are seen towards the brightest H$_2$ knots, suggesting that the majority of the observed species are associated with disks or hot core regions rather than outflows. Simple Local thermodynamic equilibrium (LTE) slab models were used to fit the observed molecular features. The LTE model fits to $^{12}$CO$_{2}$, C$_{2}$H$_{2}$, and HCN emission suggest warm $120-200$ K emission arising from a disk surface around one or both protostars. Weak $\sim500$ K H$_2$O emission at $\sim$ 6-7 $μみゅー$m is detected towards one mid-IR source, whereas $250-1050$ K H$_2$O absorption is found in the other. The H$_2$O absorption may occur in the disk atmosphere due to strong accretion-heating of the midplane, or in a disk wind viewed at an ideal angle for absorption. CO emission may originate in the hot inner disk or outflow shocks. OH emission is likely excited in a non-LTE manner through water photodissociation or chemical formation. The observations are consistent with disks having already formed in the young IRAS 23385+6053 system, but further observations are needed to disentangle the effects of geometry and evolution. △ Less

Submitted 12 January, 2024; originally announced January 2024.

Comments: 21 Pages, 16 Figures. Accepted for publication in Astronomy and Astrophysics

Observations of high-order multiplicity in a high-mass stellar protocluster

Authors: Shanghuo Li, Patricio Sanhueza, Henrik Beuther, Huei-Ru Vivien Chen, Rolf Kuiper, Fernando A. Olguin, Ralph E. Pudritz, Ian W. Stephens, Qizhou Zhang, Fumitaka Nakamura, Xing Lu, Rajika L. Kuruwita, Takeshi Sakai, Thomas Henning, Kotomi Taniguchi, Fei Li

Abstract: The dominant mechanism forming multiple stellar systems in the high-mass regime (M$_\ast \gtrsim $ 8 $M_{\odot}$) remained unknown because direct imaging of multiple protostellar systems at early phases of high-mass star formation is very challenging. High-mass stars are expected to form in clustered environments containing binaries and higher-order multiplicity systems. So far only a few high-mas… ▽ More The dominant mechanism forming multiple stellar systems in the high-mass regime (M$_\ast \gtrsim $ 8 $M_{\odot}$) remained unknown because direct imaging of multiple protostellar systems at early phases of high-mass star formation is very challenging. High-mass stars are expected to form in clustered environments containing binaries and higher-order multiplicity systems. So far only a few high-mass protobinary systems, and no definitive higher-order multiples, have been detected. Here we report the discovery of one quintuple, one quadruple, one triple and four binary protostellar systems simultaneously forming in a single high-mass protocluster, G333.23--0.06, using Atacama Large Millimeter/submillimeter Array high-resolution observations. We present a new example of a group of gravitationally bound binary and higher-order multiples during their early formation phases in a protocluster. This provides the clearest direct measurement of the initial configuration of primordial high-order multiple systems, with implications for the in situ multiplicity and its origin. We find that the binary and higher-order multiple systems, and their parent cores, show no obvious sign of disk-like kinematic structure. We conclude that the observed fragmentation into binary and higher-order multiple systems can be explained by core fragmentation, indicating its crucial role in establishing the multiplicity during high-mass star cluster formation. △ Less

Submitted 12 January, 2024; originally announced January 2024.

Comments: Published in Nature Astronomy at https://www.nature.com/articles/s41550-023-02181-9. 27 pages, 12 figures, 1 table

IPA. Class 0 Protostars Viewed in CO Emission Using JWST/NIRSpec

Authors: Adam E. Rubinstein, Himanshu Tyagi, Pooneh Nazari, Robert Gutermuth, Samuel Federman, Mayank Narang, Will R. M. Rocha, Nashanty Brunken, Katie Slavicinska, Neal J. Evans II, Joel D. Green, Dan M. Watson, Henrik Beuther, Tyler Bourke, Alessio Caratti o Garatti, Lee Hartmann, Pamela Klaassen, Hendrik Linz, Leslie W. Looney, Puravankara Manoj, S. Thomas Megeath, James Muzerolle Page, Thomas Stanke, John J. Tobin, Ewine F. van Dishoeck , et al. (2 additional authors not shown)

Abstract: We investigate the bright CO fundamental emission in the central regions of five Class 0 protostars using the JWST's Near-Infrared Spectrograph (NIRSpec) and provide clues to what processes excite the gas. CO line emission images are extracted for a forest of $\sim$150 ro-vibrational transitions from two vibrational bands, $v=1-0$ and $v=2-1$. However, ${}^{13}$CO is not detected, and thus we can… ▽ More We investigate the bright CO fundamental emission in the central regions of five Class 0 protostars using the JWST's Near-Infrared Spectrograph (NIRSpec) and provide clues to what processes excite the gas. CO line emission images are extracted for a forest of $\sim$150 ro-vibrational transitions from two vibrational bands, $v=1-0$ and $v=2-1$. However, ${}^{13}$CO is not detected, and thus we can only statistically constrain the ${}^{12}$CO optical depth. Using noise measurements to determine upper limits to the ${}^{13}$CO emission, the flux ratio of ${}^{12}$CO/${}^{13}$CO indicates that the ${}^{12}$CO emission itself is not optically thick for ro-vibrational transitions with upper state rotational quantum number $J_u \geq 15$. We construct population diagrams to estimate the rotational temperature and number of molecules from extinction-corrected CO line fluxes assuming CO emission is optically thin. Two different temperature components are required for $v=1$ ($\sim600-1000$ K and $\sim1500-3500$ K), while one hotter component is required for $v=2$ ($\sim2000-6000$ K). The vibrational temperature is $\sim 900$ K among our sources and shows no trend with luminosity. Using vibrational temperatures and the inferred total amount of CO molecules for our sources, the total warm gas mass correlates strongly with luminosity ranging from $\sim$0.1 $\rm M_{Earth}$ for the low-mass protostars to $\sim$1 M$_{\rm sun}$ for the high-mass protostars. Interpreting the distribution of gas column densities and temperatures depends on radiative and chemical processes affecting CO. The presence of a $v=2$ population may indicate CO gas radiatively excited. Selective UV photodissociation of CO isotopologues around our high-mass sources may explain their depletion of ${}^{13}$CO. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: 30 pages, 6 figures, 5 tables, submitted to ApJ

JWST Observations of Young protoStars (JOYS+): Detection of icy complex organic molecules and ions. I. CH$_4$, SO$_2$, HCOO$^-$, OCN$^-$, H$_2$CO, HCOOH, CH$_3$CH$_2$OH, CH$_3$CHO, CH$_3$OCHO, CH$_3$COOH

Authors: W. R. M. Rocha, E. F. van Dishoeck, M. E. Ressler, M. L. van Gelder, K. Slavicinska, N. G. C. Brunken, H. Linnartz, T. P. Ray, H. Beuther, A. Caratti o Garatti, V. Geers, P. J. Kavanagh, P. D. Klaassen, K. Justannont, Y. Chen, L. Francis, C. Gieser, G. Perotti, Ł. Tychoniec, M. Barsony, L. Majumdar, V. J. M. le Gouellec, L. E. U. Chu, B. W. P. Lew, Th. Henning , et al. (1 additional authors not shown)

Abstract: Complex organic molecules (COMs) detected in the gas phase are thought to be mostly formed on icy grains, but no unambiguous detection of icy COMs larger than CH3OH has been reported so far. Exploring this matter in more detail has become possible with the JWST the critical 5-10 $μみゅー$m range. In the JOYS+ program, more than 30 protostars are being observed with the MIRI/MRS. This study explores the… ▽ More Complex organic molecules (COMs) detected in the gas phase are thought to be mostly formed on icy grains, but no unambiguous detection of icy COMs larger than CH3OH has been reported so far. Exploring this matter in more detail has become possible with the JWST the critical 5-10 $μみゅー$m range. In the JOYS+ program, more than 30 protostars are being observed with the MIRI/MRS. This study explores the COMs ice signatures in the low and high-mass protostar, IRAS 2A and IRAS 23385, respectively. We fit continuum and silicate subtracted observational data with IR laboratory ice spectra. We use the ENIIGMA fitting tool to find the best fit between the lab data and the observations and to performs statistical analysis of the solutions. We report the best fits for the spectral ranges between 6.8 and 8.6 $μみゅー$m in IRAS 2A and IRAS 23385, originating from simple molecules, COMs, and negative ions. The strongest feature in this range (7.7 $μみゅー$m) is dominated by CH4 and has contributions of SO2 and OCN-. Our results indicate that the 7.2 and 7.4 $μみゅー$m bands are mostly dominated by HCOO-. We find statistically robust detections of COMs based on multiple bands, most notably CH3CHO, CH3CH2OH, and CH3OCHO. The likely detection of CH3COOH is also reported. The ice column density ratios between CH3CH2OH and CH3CHO of IRAS 2A and IRAS 23385, suggests that these COMs are formed on icy grains. Finally, the derived ice abundances for IRAS 2A correlate well with those in comet 67P/GC within a factor of 5. Based on the MIRI/MRS data, we conclude that COMs are present in interstellar ices, thus providing additional proof for a solid-state origin of these species in star-forming regions. The good correlation between the ice abundances in comet 67P and IRAS 2A is in line with the idea that cometary COMs can be inherited from the early protostellar phases. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: Accepted for publication in A&A

JOYS+: mid-infrared detection of gas-phase SO$_2$ emission in a low-mass protostar. The case of NGC 1333 IRAS2A: hot core or accretion shock?

Authors: M. L. van Gelder, M. E. Ressler, E. F. van Dishoeck, P. Nazari, B. Tabone, J. H. Black, Ł. Tychoniec, L. Francis, M. Barsony, H. Beuther, A. Caratti o Garatti, Y. Chen, C. Gieser, V. J. M. le Gouellec, P. J. Kavanagh, P. D. Klaassen, B. W. P. Lew, H. Linnartz, L. Majumdar, G. Perotti, W. R. M. Rocha

Abstract: JWST/MIRI has sharpened our infrared eyes toward the star formation process. This paper presents the first mid-infrared detection of gaseous SO$_2$ emission in an embedded low-mass protostellar system. MIRI-MRS observations of the low-mass protostellar binary NGC 1333 IRAS2A are presented from the JWST Observations of Young protoStars (JOYS+) program, revealing emission from the SO$_2~νにゅー_3$ asymmet… ▽ More JWST/MIRI has sharpened our infrared eyes toward the star formation process. This paper presents the first mid-infrared detection of gaseous SO$_2$ emission in an embedded low-mass protostellar system. MIRI-MRS observations of the low-mass protostellar binary NGC 1333 IRAS2A are presented from the JWST Observations of Young protoStars (JOYS+) program, revealing emission from the SO$_2~νにゅー_3$ asymmetric stretching mode at 7.35 micron. The results are compared to those derived from high-angular resolution SO$_2$ data obtained with ALMA. The SO$_2$ emission from the $νにゅー_3$ band is predominantly located on $\sim50-100$ au scales around the main component of the binary, IRAS2A1. A rotational temperature of $92\pm8$ K is derived from the $νにゅー_3$ lines. This is in good agreement with the rotational temperature derived from pure rotational lines in the vibrational ground state (i.e., $νにゅー=0$) with ALMA ($104\pm5$ K). However, the emission of the $νにゅー_3$ lines is not in LTE given that the total number of molecules predicted by a LTE model is found to be a factor $2\times10^4$ higher than what is derived for the $νにゅー=0$ state. This difference can be explained by a vibrational temperature that is $\sim100$ K higher than the derived rotational temperature of the $νにゅー=0$ state. The brightness temperature derived from the continuum around the $νにゅー_3$ band of SO$_2$ is $\sim180$ K, which confirms that the $νにゅー_3=1$ level is not collisionally populated but rather infrared pumped by scattered radiation. This is also consistent with the non-detection of the $νにゅー_2$ bending mode at 18-20 micron. Given the rotational temperature, the extent of the emission ($\sim100$ au in radius), and the narrow line widths in the ALMA data (3.5 km/s), the SO$_2$ in IRAS2A likely originates from ice sublimation in the central hot core around the protostar rather than from an accretion shock at the disk-envelope boundary. △ Less

Submitted 28 November, 2023; originally announced November 2023.

Comments: 19 pages, 17 figures, accepted for publication in A&A, abstract abbreviated

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

A Collection of German Science Interests in the Next Generation Very Large Array

Authors: M. Kadler, D. A. Riechers, J. Agarwal, A. -K. Baczko, H. Beuther, F. Bigiel, T. Birnstiel, B. Boccardi, D. J. Bomans, L. Boogaard, T. T. Braun, S. Britzen, M. Brüggen, A. Brunthaler, P. Caselli, D. Elsässer, S. von Fellenberg, M. Flock, C. M. Fromm, L. Fuhrmann, P. Hartogh, M. Hoeft, R. P. Keenan, Y. Kovalev, K. Kreckel , et al. (66 additional authors not shown)

Abstract: The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy com… ▽ More The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy community under the lead of the National Radio Astronomy Observatory (NRAO), and is anticipated to be built between 2027 and 2037. Two workshops have been held in 2022 and 2023 with the goal to discuss and consolidate the scientific interests in the ngVLA within the German astronomical community. This community paper constitutes a collection of 48 science ideas which the German community aims to pursue with the ngVLA in the 2030s. This is not a complete list and the ideas are not developed at the level of a "Science Book", such that the present document is mainly meant provide a basis for further discussion within the community. As such, additional contributions are welcome, and will be considered for inclusion in future revisions. △ Less

Submitted 18 June, 2024; v1 submitted 16 November, 2023; originally announced November 2023.

Comments: Version 2.0 (status June 18, 2024): 169 pages, comments and future contributions welcome [v2.0: 7 new science cases added, some minor revisions to other chapters]

Discovery of a collimated jet from the low luminosity protostar IRAS 16253$-$2429 in a quiescent accretion phase with the JWST

Authors: Mayank Narang, Manoj P., Himanshu Tyagi, Dan M. Watson, S. Thomas Megeath, Samuel Federman, Adam E. Rubinstein, Robert Gutermuth, Alessio Caratti o Garatti, Henrik Beuther, Tyler L. Bourke, Ewine F. Van Dishoeck, Neal J. Evans II, Guillem Anglada, Mayra Osorio, Thomas Stanke, James Muzerolle, Leslie W. Looney, Yao-Lun Yang, John J. Tobin, Pamela Klaassen, Nicole Karnath, Prabhani Atnagulov, Nashanty Brunken, William J. Fischer , et al. (14 additional authors not shown)

Abstract: Investigating Protostellar Accretion (IPA) is a JWST Cycle~1 GO program that uses NIRSpec IFU and MIRI MRS to obtain 2.9--28~$μみゅー$m spectral cubes of young, deeply embedded protostars with luminosities of 0.2 to 10,000~L$_{\odot}$ and central masses of 0.15 to 12~M$_{\odot}$. In this Letter, we report the discovery of a highly collimated atomic jet from the Class~0 protostar IRAS~16253$-$2429, the l… ▽ More Investigating Protostellar Accretion (IPA) is a JWST Cycle~1 GO program that uses NIRSpec IFU and MIRI MRS to obtain 2.9--28~$μみゅー$m spectral cubes of young, deeply embedded protostars with luminosities of 0.2 to 10,000~L$_{\odot}$ and central masses of 0.15 to 12~M$_{\odot}$. In this Letter, we report the discovery of a highly collimated atomic jet from the Class~0 protostar IRAS~16253$-$2429, the lowest luminosity source ($L_\mathrm{bol}$ = 0.2 $L_\odot$) in the IPA program. The collimated jet is detected in multiple [Fe~II] lines, [Ne~II], [Ni~II], and H~I lines, but not in molecular emission. The atomic jet has a velocity of about 169~$\pm$~15~km\,s$^{-1}$, after correcting for inclination. The width of the jet increases with distance from the central protostar from 23 to~60 au, corresponding to an opening angle of 2.6~$\pm$~0.5\arcdeg. By comparing the measured flux ratios of various fine structure lines to those predicted by simple shock models, we derive a shock {speed} of 54~km\,s$^{-1}$ and a preshock density of 2.0$\times10^{3}$~cm$^{-3}$ at the base of the jet. {From these quantities and using a suite of jet models and extinction laws we compute a mass loss rate between $0.4 -1.1\times10^{-10}~M_{\odot}$~yr~$^{-1}$.} The low mass loss rate is consistent with simultaneous measurements of low mass accretion rate ($2.4~\pm~0.8~\times~10^{-9}~M_{\odot}$~yr$^{-1}$) for IRAS~16253$-$2429 from JWST observations (Watson et al. in prep), indicating that the protostar is in a quiescent accretion phase. Our results demonstrate that very low-mass protostars can drive highly collimated, atomic jets, even during the quiescent phase. △ Less

Submitted 11 January, 2024; v1 submitted 21 October, 2023; originally announced October 2023.

Comments: Accepted to ApJL. Comments and feedback welcome

Far-Infrared Luminosity Bursts Trace Mass Accretion onto Protostars

Authors: William J. Fischer, Cara Battersby, Doug Johnstone, Rachel Lee, Marta Sewilo, Henrik Beuther, Yasuhiro Hasegawa, Adam Ginsburg, Klaus Pontoppidan

Abstract: Evidence abounds that young stellar objects undergo luminous bursts of intense accretion that are short compared to the time it takes to form a star. It remains unclear how much these events contribute to the main-sequence masses of the stars. We demonstrate the power of time-series far-infrared (far-IR) photometry to answer this question compared to similar observations at shorter and longer wave… ▽ More Evidence abounds that young stellar objects undergo luminous bursts of intense accretion that are short compared to the time it takes to form a star. It remains unclear how much these events contribute to the main-sequence masses of the stars. We demonstrate the power of time-series far-infrared (far-IR) photometry to answer this question compared to similar observations at shorter and longer wavelengths. We start with model spectral energy distributions that have been fit to 86 Class 0 protostars in the Orion molecular clouds. The protostars sample a broad range of envelope densities, cavity geometries, and viewing angles. We then increase the luminosity of each model by factors of 10, 50, and 100 and assess how these luminosity increases manifest in the form of flux increases over wavelength ranges of interest. We find that the fractional change in the far-IR luminosity during a burst more closely traces the change in the accretion rate than photometric diagnostics at mid-infrared and submillimeter wavelengths. We also show that observations at far-IR and longer wavelengths reliably track accretion changes without confusion from large, variable circumstellar and interstellar extinction that plague studies at shorter wavelengths. We close by discussing the ability of a proposed far-IR surveyor for the 2030s to enable improvements in our understanding of the role of accretion bursts in mass assembly. △ Less

Submitted 22 December, 2023; v1 submitted 19 October, 2023; originally announced October 2023.

Comments: Replaced with version accepted by AAS Journals

A global view on star formation: The GLOSTAR Galactic plane survey. IX. Radio Source Catalog III: 2<l<28, 36<l<40, 56<l<60 and |b|<1, VLA B-configuration

Authors: A. Y. Yang, S. A. Dzib, J. S. Urquhart, A. Brunthaler, S. -N. X. Medina, K. M. Menten, F. Wyrowski, G. N. Ortiz-León, W. D. Cotton, Y. Gong, R. Dokara, M. R. Rugel, H. Beuther, J. D. Pandian, T. Csengeri, V. S. Veena, N. Roy, H. Nguyen, B. Winkel, J. Ott, C. Carrasco-Gonzalez, S. Khan, A. Cheema

Abstract: As part of the GLOSTAR (GLObal view of STAR formation in the Milky Way) survey, we present the high-resolution continuum source catalog for the regions (l = 2-28, 36-40, 56-60, &|b|<1.0), observed with the Karl G. Jansky Very Large Array (VLA) in its B-configuration. The continuum images are optimized to detect compact sources on angular scales up to 4", and have a typical noise level of 1sigma ~… ▽ More As part of the GLOSTAR (GLObal view of STAR formation in the Milky Way) survey, we present the high-resolution continuum source catalog for the regions (l = 2-28, 36-40, 56-60, &|b|<1.0), observed with the Karl G. Jansky Very Large Array (VLA) in its B-configuration. The continuum images are optimized to detect compact sources on angular scales up to 4", and have a typical noise level of 1sigma ~ 0.08mJy/beam for an angular resolution of 1", which makes GLOSTAR currently the highest resolution as well as the most sensitive radio survey of the northern Galactic plane at 4-8GHz. We extracted 13354 sources above a threshold of 5sigma and 5437 sources above 7sigma that represent the high-reliability catalog. We determined the in-band spectral index (alpha) for the sources in the 7sigma-threshold catalog. The mean value is alpha=-0.6, which indicates that the catalog is dominated by sources emitting non-thermal radio emission. We identified the most common source types detected in radio surveys: 251 HII region candidates (113 new), 282 planetary nebulae (PNe) candidates (127 new), 784 radio star candidates (581 new), and 4080 extragalactic radio source candidates (2175 new). A significant fraction of HII regions and PNe candidates have alpha<-0.1 indicating that these candidates could contain radio jets, winds or outflows from high-mass and low-mass stellar objects. We identified 245 variable radio sources by comparing the flux densities of compact sources from the GLOSTAR survey and the Co-Ordinated Radio `N' Infrared Survey for High-mass star formation (CORNISH), and find that most of them are infrared quiet. The catalog is typically 95% complete for point sources at a flux density of 0.6 mJy (i.e. typical 7sigma level) and the systematic positional uncertainty is <= 0.1". The GLOSTAR data and catalogs are available online at https://glostar.mpifr-bonn.mpg.de. △ Less

Submitted 23 October, 2023; v1 submitted 15 October, 2023; originally announced October 2023.

Comments: 25pages, 21 figures, has been accepted for publication in Astronomy & Astrophysics (A&A)

Journal ref: A&A 680, A92 (2023)

Investigating Protostellar Accretion-Driven Outflows Across the Mass Spectrum: JWST NIRSpec IFU 3-5~$μみゅー$m Spectral Mapping of Five Young Protostars

Authors: Samuel Federman, S. Thomas Megeath, Adam E. Rubinstein, Robert Gutermuth, Mayank Narang, Himanshu Tyagi, P. Manoj, Guillem Anglada, Prabhani Atnagulov, Henrik Beuther, Tyler L. Bourke, Nashanty Brunken, Alessio Caratti o Garatti, Neal J. Evans II, William J. Fischer, Elise Furlan, Joel Green, Nolan Habel, Lee Hartmann, Nicole Karnath, Pamela Klaassen, Hendrik Linz, Leslie W. Looney, Mayra Osorio, James Muzerolle Page , et al. (13 additional authors not shown)

Abstract: Investigating Protostellar Accretion is a Cycle 1 JWST program using the NIRSpec+MIRI integral field units to obtain 2.9--28 $μみゅー$m spectral cubes of five young protostars with luminosities of 0.2-10,000 L$_{\odot}$ in their primary accretion phase. This paper introduces the NIRSpec 2.9--5.3 $μみゅー$m data of the inner 840-9000 au with spatial resolutions from 28-300 au. The spectra show rising continuum… ▽ More Investigating Protostellar Accretion is a Cycle 1 JWST program using the NIRSpec+MIRI integral field units to obtain 2.9--28 $μみゅー$m spectral cubes of five young protostars with luminosities of 0.2-10,000 L$_{\odot}$ in their primary accretion phase. This paper introduces the NIRSpec 2.9--5.3 $μみゅー$m data of the inner 840-9000 au with spatial resolutions from 28-300 au. The spectra show rising continuum emission; deep ice absorption; emission from H$_{2}$, H~I, and [Fe~II]; and the CO fundamental series in emission and absorption. Maps of the continuum emission show scattered light cavities for all five protostars. In the cavities, collimated jets are detected in [Fe~II] for the four $< 320$~L$_{\odot}$ protostars, two of which are additionally traced in Br-$αあるふぁ$. Knots of [Fe~II] emission are detected toward the most luminous protostar, and knots of [FeII] emission with dynamical times of $< 30$~yrs are found in the jets of the others. While only one jet is traced in H$_2$, knots of H$_2$ and CO are detected in the jets of four protostars. H$_2$ is seen extending through the cavities, showing that they are filled by warm molecular gas. Bright H$_2$ emission is seen along the walls of a single cavity, while in three cavities narrow shells of H$_2$ emission are found, one of which has an [Fe~II] knot at its apex. These data show cavities containing collimated jets traced in atomic/ionic gas surrounded by warm molecular gas in a wide-angle wind and/or gas accelerated by bow shocks in the jets. △ Less

Submitted 24 April, 2024; v1 submitted 5 October, 2023; originally announced October 2023.

Comments: 26 pages, 12 figures

Journal ref: ApJ 966 41 (2024)

Turbulent Structure In Supernova Remnants G46.8-0.3 And G39.2-0.3 From THOR Polarimetry

Authors: Russell Shanahan, Jeroen Stil, Loren Anderson, Henrik Beuther, Paul Goldsmith, Ralf S. Klessen, Michael Rugel, Juan D. Soler

Abstract: We present the continued analysis of polarization and Faraday rotation for the supernova remnants (SNRs) G46.8-0.3 and G39.2-0.3 in L-band (1-2 GHz) radio continuum in The HI/OH/Recombination line (THOR) survey. In this work, we present our investigation of Faraday depth fluctuations from angular scales comparable to the size of the SNRs down to scales less than our 16" beam (<~0.7 pc) from Farada… ▽ More We present the continued analysis of polarization and Faraday rotation for the supernova remnants (SNRs) G46.8-0.3 and G39.2-0.3 in L-band (1-2 GHz) radio continuum in The HI/OH/Recombination line (THOR) survey. In this work, we present our investigation of Faraday depth fluctuations from angular scales comparable to the size of the SNRs down to scales less than our 16" beam (<~0.7 pc) from Faraday dispersion (sigma_phi). From THOR, we find median sigma_phi of 15.9 +/- 3.2 rad m^-2 for G46.8-0.3 and 17.6 +/- 1.6 rad m^-2 for G39.2-0.3. When comparing to polarization at 6cm, we find evidence for sigma_phi >~ 30 rad m^-2 in localized regions where we detect no polarization in THOR. We combine Faraday depth dispersion with the rotation measure (RM) structure function (SF) and find evidence for a break in the SF on scales less than the THOR beam. We estimate the RM SF of the foreground interstellar medium (ISM) using the SF of extra-Galactic radio sources (EGRS) and pulsars to find that the RM fluctuations we measure originate within the SNRs for all but the largest angular scales. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: Has been accepted by ApJ for publication. Figures 3 and 4 did do not render well when using an internet browser to view them, but when the pdf file is downloaded these figures look as they should

Cold atomic gas identified by HI self-absorption. Cold atomic clouds toward giant molecular filaments

Authors: J. Syed, H. Beuther, P. F. Goldsmith, Th. Henning, M. Heyer, R. S. Klessen, J. M. Stil, J. D. Soler, L. D. Anderson, J. S. Urquhart, M. R. Rugel, K. G. Johnston, A. Brunthaler

Abstract: Stars form in the dense interiors of molecular clouds. The dynamics and physical properties of the atomic interstellar medium (ISM) set the conditions under which molecular clouds and eventually stars will form. It is, therefore, critical to investigate the relationship between the atomic and molecular gas phase to understand the global star formation process. Using the high angular resolution dat… ▽ More Stars form in the dense interiors of molecular clouds. The dynamics and physical properties of the atomic interstellar medium (ISM) set the conditions under which molecular clouds and eventually stars will form. It is, therefore, critical to investigate the relationship between the atomic and molecular gas phase to understand the global star formation process. Using the high angular resolution data from The HI/OH/Recombination line survey of the Milky Way (THOR), we aim to constrain the kinematic and physical properties of the cold atomic hydrogen gas phase toward the inner Galactic plane. HI self-absorption (HISA) has proven to be a viable method to detect cold atomic hydrogen clouds in the Galactic plane. With the help of a newly developed self-absorption extraction routine (astroSABER), we build upon previous case studies to identify HI self-absorption toward a sample of Giant Molecular Filaments (GMFs). We find the cold atomic gas to be spatially correlated with the molecular gas on a global scale. The column densities of the cold atomic gas traced by HISA are usually of the order of $10^{20}\rm\,cm^{-2}$ whereas those of molecular hydrogen traced by $\rm^{13}CO$ are at least an order of magnitude higher. The HISA column densities are attributed to a cold gas component that accounts for a fraction of $\sim$5% of the total atomic gas budget within the clouds. The HISA column density distributions show pronounced log-normal shapes that are broader than those traced by HI emission. The cold atomic gas is found to be moderately supersonic with Mach numbers of a $\sim$few. In contrast, highly supersonic dynamics drive the molecular gas within most filaments. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: 41 pages, 28 figures, accepted for publication in A&A

Journal ref: A&A 679, A130 (2023)

JOYS: Disentangling the warm and cold material in the high-mass IRAS 23385+6053 cluster

Authors: C. Gieser, H. Beuther, E. F. van Dishoeck, L. Francis, M. L. van Gelder, L. Tychoniec, P. J. Kavanagh, G. Perotti, A. Caratti o Garatti, T. P. Ray, P. Klaassen, K. Justtanont, H. Linnartz, W. R. M. Rocha, K. Slavicinska, L. Colina, M. Güdel, Th. Henning, P. -O. Lagage, G. Östlin, B. Vandenbussche, C. Waelkens, G. Wright

Abstract: (abridged) We study and compare the warm (>100 K) and cold (<100 K) material toward the high-mass star-forming region IRAS 23385+6053 (IRAS 23385 hereafter) combining high angular resolution observations in the mid-infrared (MIR) with the JWST Observations of Young protoStars (JOYS) project and with the NOEMA at mm wavelengths at angular resolutions of 0.2"-1". The spatial morphology of atomic a… ▽ More (abridged) We study and compare the warm (>100 K) and cold (<100 K) material toward the high-mass star-forming region IRAS 23385+6053 (IRAS 23385 hereafter) combining high angular resolution observations in the mid-infrared (MIR) with the JWST Observations of Young protoStars (JOYS) project and with the NOEMA at mm wavelengths at angular resolutions of 0.2"-1". The spatial morphology of atomic and molecular species is investigated by line integrated intensity maps. The temperature and column density of different gas components is estimated using H2 transitions (warm and hot component) and a series of CH3CN transitions as well as 3 mm continuum emission (cold component). Toward the central dense core in IRAS 23385 the material consists of relatively cold gas and dust (~50 K), while multiple outflows create heated and/or shocked H2 and show enhanced temperatures (~400 K) along the outflow structures. An energetic outflow with enhanced emission knots of [Fe II] and [Ni II] hints at J-type shocks, while two other outflows have enhanced emission of only H2 and [S I] caused by C-type shocks. The latter two outflows are also more prominent in molecular line emission at mm wavelengths (e.g., SiO, SO, H2CO, and CH3OH). Even higher angular resolution data are needed to unambiguously identify the outflow driving sources given the clustered nature of IRAS 23385. While most of the forbidden fine structure transitions are blueshifted, [Ne II] and [Ne III] peak at the source velocity toward the MIR source A/mmA2 suggesting that the emission is originating from closer to the protostar. △ Less

Submitted 19 September, 2023; originally announced September 2023.

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

Journal ref: A&A 679, A108 (2023)

The Cygnus Allscale Survey of Chemistry and Dynamical Environments: CASCADE. II. A detailed kinematic analysis of the DR21 Main outflow

Authors: I. M. Skretas, A. Karska, F. Wyrowski, K. M. Menten, H. Beuther, A. Ginsburg, A. Hernández-Gómez, C. Gieser, S. Li, W. -J. Kim, D. A. Semenov, L. Bouscasse, I. B. Christensen, J. M. Winters, A. Hacar

Abstract: Molecular outflows are believed to be a key ingredient in the process of star formation. The molecular outflow associated with DR21 Main in Cygnus-X is one of the most extreme, in mass and size, molecular outflows in the Milky Way. The outflow is suggested to belong to a rare class of explosive outflows which are formed by the disintegration of protostellar systems.We aim to explore the morphology… ▽ More Molecular outflows are believed to be a key ingredient in the process of star formation. The molecular outflow associated with DR21 Main in Cygnus-X is one of the most extreme, in mass and size, molecular outflows in the Milky Way. The outflow is suggested to belong to a rare class of explosive outflows which are formed by the disintegration of protostellar systems.We aim to explore the morphology, kinematics,and energetics of the DR21 Main outflow, and compare those properties to confirmed explosive outflows to unravel the underlying driving mechanism behind DR21. Line and continuum emission are studied at a wavelength of 3.6\,mm with IRAM 30 m and NOEMA telescopes as part of the Cygnus Allscale Survey of Chemistry and Dynamical Environments (CASCADE) program. The spectra include ($J= 1-0$) transitions of HCO$^+$, HCN, HNC, N$_2$H$^+$, H$_2$CO, CCH tracing different temperature and density regimes of the outflowing gas at high-velocity resolution ($\sim$ 0.8 km s$^{-1}$). The map encompasses the entire DR21 Main outflow and covers all spatial scales down to a resolution of ~3" ($\sim$ 0.02 pc). Integrated intensity maps of the HCO$^+$ emission reveal a strongly collimated bipolar outflow with significant overlap of the blue- and red-shifted emission. The opening angles of both outflow lobes decrease with velocity, from $\sim80$ to 20$^{\circ}$ for the velocity range from 5 to 45 km s$^{-1}$ relative to the source velocity. No evidence is found for the presence of elongated, "filament-like" structures expected in explosive outflows. N$_2$H$^+$ emission near the western outflow lobe reveals the presence of a dense molecular structure which appears to be interacting with the DR21 Main outflow. The overall morphology as well as the detailed kinematics of the DR21 Main outflow is more consistent with that of a typical bipolar outflow instead of an explosive counterpart. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: 26 pages, 37 figures. Accepted for publication in A&A

Journal ref: A&A 679, A66 (2023)

Metrewave Galactic Plane with the uGMRT (MeGaPluG) Survey: Lessons from the Pilot Study

Authors: Rohit Dokara, Nirupam Roy, Karl Menten, Sarita Vig, Prasun Dutta, Henrik Beuther, Jagadheep D. Pandian, Michael Rugel, Md Rashid, Andreas Brunthaler

Abstract: Context. The advent of wide-band receiver systems on interferometer arrays enables one to undertake high-sensitivity and high-resolution radio continuum surveys of the Galactic plane in a reasonable amount of telescope time. However, to date, there are only a few such studies of the first quadrant of the Milky Way that have been carried out at frequencies below 1 GHz. The Giant Metrewave Radio Tel… ▽ More Context. The advent of wide-band receiver systems on interferometer arrays enables one to undertake high-sensitivity and high-resolution radio continuum surveys of the Galactic plane in a reasonable amount of telescope time. However, to date, there are only a few such studies of the first quadrant of the Milky Way that have been carried out at frequencies below 1 GHz. The Giant Metrewave Radio Telescope (GMRT) has recently upgraded its receivers with wide-band capabilities (now called the uGMRT) and provides a good opportunity to conduct high resolution surveys, while also being sensitive to the extended structures. Aims. We wish to assess the feasibility of conducting a large-scale snapshot survey, the Metrewave Galactic Plane with the uGMRT Survey (MeGaPluG), to simultaneously map extended sources and compact objects at an angular resolution lower than $10''$ and a point source sensitivity of 0.15 mJy/beam. Methods. We performed an unbiased survey of a small portion of the Galactic plane, covering the W43/W44 regions ($l=29^\circ-35^\circ$ and $|b|<1^\circ$) in two frequency bands: 300$-$500 MHz and 550$-$750 MHz. The 200 MHz wide-band receivers on the uGMRT are employed to observe the target field in several pointings, spending nearly 14 minutes on each pointing in two separate scans. We developed an automated pipeline for the calibration, and a semi-automated self-calibration procedure is used to image each pointing using multi-scale CLEAN and outlier fields. Results. We produced continuum mosaics of the surveyed region at a final common resolution of $25''$ in the two bands that have central frequencies of 400 MHz and 650 MHz, with a point source sensitivity better than 5 mJy/beam. We plan to cover a larger footprint of the Galactic plane in the near future based on the lessons learnt from this study. (Abridged) △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: To be published in A&A. 13 pages, 10 figures

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

A global view on star formation: The GLOSTAR Galactic plane survey VIII. Formaldehyde absorption in Cygnus~X

Authors: Y. Gong, G. N. Ortiz-León, M. R. Rugel, K. M. Menten, A. Brunthaler, F. Wyrowski, C. Henkel, H. Beuther, S. A. Dzib, J. S. Urquhart, A. Y. Yang, J. D. Pandian, R. Dokara, V. S. Veena, H. Nguyen, S. -N. X. Medina, W. D. Cotton, W. Reich, B. Winkel, P. Müller, I. Skretas, T. Csengeri, S. Khan, A. Cheema

Abstract: Cygnus X is one of the closest and most active high-mass star-forming regions in our Galaxy, making it one of the best laboratories for studying massive star formation. As part of the GLOSTAR Galactic plane survey, we performed large scale simultaneous H$_{2}$CO (1$_{1,0}$-1$_{1,1}$) spectral line and radio continuum imaging observations toward Cygnus X at $λらむだ\sim$6 cm with the Karl G. Jansky Very… ▽ More Cygnus X is one of the closest and most active high-mass star-forming regions in our Galaxy, making it one of the best laboratories for studying massive star formation. As part of the GLOSTAR Galactic plane survey, we performed large scale simultaneous H$_{2}$CO (1$_{1,0}$-1$_{1,1}$) spectral line and radio continuum imaging observations toward Cygnus X at $λらむだ\sim$6 cm with the Karl G. Jansky Very Large Array and the Effelsberg-100 m radio telescope. Our Effelsberg observations reveal widespread H$_{2}$CO (1$_{1,0}$-1$_{1,1}$) absorption with a spatial extent of $\gtrsim$50 pc in Cygnus~X for the first time. On large scales of 4.4 pc, the relative orientation between local velocity gradient and magnetic field tends to be more parallel at H$_{2}$ column densities of $\gtrsim$1.8$\times 10^{22}$~cm$^{-2}$. On the smaller scale of 0.17 pc, our VLA+Effelsberg combined data reveal H$_{2}$CO absorption only toward three bright H{\scriptsize II} regions. Our observations demonstrate that H$_{2}$CO (1$_{1,0}$-1$_{1,1}$) is commonly optically thin. Kinematic analysis supports the assertion that molecular clouds generally exhibit supersonic motions on scales of 0.17-4.4 pc. We show a non-negligible contribution of the cosmic microwave background radiation in producing extended absorption features in Cygnus X. Our observations suggest that H$_{2}$CO ($1_{1,0}-1_{1,1}$) can trace molecular gas with H$_{2}$ column densities of $\gtrsim 5 \times 10^{21}$ cm$^{-2}$. The ortho-H$_{2}$CO fractional abundance with respect to H$_{2}$ has a mean value of 7.0$\times 10^{-10}$. A comparison of velocity dispersions on different linear scales suggests that the dominant $-3$ km s$^{-1}$ velocity component in the prominent DR21 region has nearly identical velocity dispersions on scales of 0.17-4.4 pc, which deviates from the expected behavior of classic turbulence. △ Less

Submitted 2 August, 2023; originally announced August 2023.

Comments: 27 pages, 23 figures, accepted for publication in A&A

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

The diverse chemistry of protoplanetary disks as revealed by JWST

Authors: Ewine F. van Dishoeck, S. Grant, B. Tabone, M. van Gelder, L. Francis, L. Tychoniec, G. Bettoni, A. M. Arabhavi, D. Gasman, P. Nazari, M. Vlasblom, P. Kavanagh, V. Christiaens, P. Klaassen, H. Beuther, Th. Henning, I. Kamp

Abstract: Early results from the JWST-MIRI guaranteed time programs on protostars (JOYS) and disks (MINDS) are presented. Thanks to the increased sensitivity, spectral and spatial resolution of the MIRI spectrometer, the chemical inventory of the planet-forming zones in disks can be investigated with unprecedented detail across stellar mass range and age. Here data are presented for five disks, four around… ▽ More Early results from the JWST-MIRI guaranteed time programs on protostars (JOYS) and disks (MINDS) are presented. Thanks to the increased sensitivity, spectral and spatial resolution of the MIRI spectrometer, the chemical inventory of the planet-forming zones in disks can be investigated with unprecedented detail across stellar mass range and age. Here data are presented for five disks, four around low-mass stars and one around a very young high-mass star. The mid-infrared spectra show some similarities but also significant diversity: some sources are rich in CO2, others in H2O or C2H2. In one disk around a very low-mass star, booming C2H2 emission provides evidence for a ``soot'' line at which carbon grains are eroded and sublimated, leading to a rich hydrocarbon chemistry in which even di-acetylene (C4H2) and benzene (C6H6) are detected (Tabone et al. 2023). Together, the data point to an active inner disk gas-phase chemistry that is closely linked to the physical structure (temperature, snowlines, presence of cavities and dust traps) of the entire disk and which may result in varying CO2/H2O abundances and high C/O ratios >1 in some cases. Ultimately, this diversity in disk chemistry will also be reflected in the diversity of the chemical composition of exoplanets. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: 17 pages, 8 figures. Author's version of paper submitted to Faraday Discussions January 18 2023, Accepted March 16 2023

Journal ref: Faraday Discussion "Astrochemistry at high resolution", 2023

Large-scale Velocity-coherent Filaments in the SEDIGISM Survey: Association with Spiral Arms and Fraction of Dense Gas

Authors: Y. Ge, K. Wang, A. Duarte-Cabral, A. R. Pettitt, C. L. Dobbs, Á. Sánchez-Monge, K. R. Neralwar, J. S. Urquhart, D. Colombo, E. Durán-Camacho, H. Beuther, L. Bronfman, A. J. Rigby, D. Eden, S. Neupane, P. Barnes, T. Henning, A. Y. Yang

Abstract: Context. Filamentary structures in the interstellar medium are closely related to star formation. Dense gas mass fraction (DGMF) or clump formation efficiency in large-scale filaments possibly determine their hosting star formation activities. Aims. We aim to automatically identify large-scale filaments, characterize them, investigate their association with Galactic structures, and study their DGM… ▽ More Context. Filamentary structures in the interstellar medium are closely related to star formation. Dense gas mass fraction (DGMF) or clump formation efficiency in large-scale filaments possibly determine their hosting star formation activities. Aims. We aim to automatically identify large-scale filaments, characterize them, investigate their association with Galactic structures, and study their DGMFs. Methods. We use a modified minimum spanning tree (MST) algorithm to chain parsec-scale 13CO clumps previously extracted from the SEDIGISM (Structure, Excitation, and Dynamics of the Inner Galactic InterStellar Medium) survey. The MST connects nodes in a graph such that the sum of edge lengths is minimum. Modified MST also ensures velocity coherence between nodes, so the identified filaments are coherent in position-position-velocity (PPV) space. Results. We generate a catalog of 88 large-scale ($>10pc$) filaments in the inner Galactic plane (with $-60^\circ < l < 18^\circ and $|b| < 0.5^\circ$). These SEDIGISM filaments are larger and less dense than MST filaments previously identified from the BGPS and ATLASGAL surveys. We find that eight of the filaments run along spiral arms and can be regarded as "bones" of the Milky Way. We also find three bones associated with the Local Spur in PPV space. By compiling 168 large-scale filaments with available DGMF across the Galaxy, an order of magnitude more than previously investigated, we find that DGMFs do not correlate with Galactic location, but bones have higher DGMFs than other filaments. △ Less

Submitted 12 May, 2023; originally announced May 2023.

Comments: Accepted for publication in A&A. 30 pages, 8 figures, 3 tables

Journal ref: A&A 675, A119 (2023)

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

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

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

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

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

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

Physical and chemical complexity in high-mass star-forming regions with ALMA. I. Overview and evolutionary trends of physical properties

Authors: Caroline Gieser, Henrik Beuther, Dmitry Semenov, Aida Ahmadi, Thomas Henning, Molly R. A. Wells

Abstract: In this study, we investigate how physical properties, such as the density and temperature profiles, evolve on core scales through the evolutionary sequence during high-mass star formation ranging from protostars in cold infrared dark clouds to evolved UCHII regions. We observed 11 high-mass star-forming regions with ALMA at 3 mm wavelengths. Based on the 3 mm continuum morphology and recombinatio… ▽ More In this study, we investigate how physical properties, such as the density and temperature profiles, evolve on core scales through the evolutionary sequence during high-mass star formation ranging from protostars in cold infrared dark clouds to evolved UCHII regions. We observed 11 high-mass star-forming regions with ALMA at 3 mm wavelengths. Based on the 3 mm continuum morphology and recombination line emission, tracing locations with free-free (ff) emission, the fragmented cores analyzed in this study are classified into either dust or dust+ff cores. In addition, we resolve three cometary UCHII regions with extended 3 mm emission that is dominated by free-free emission. The temperature structure and radial profiles (T~r^-q ) are determined by modeling molecular emission of CH3CN and CH313CN with XCLASS and by using the HCN-to- HNC intensity ratio as probes for the gas kinetic temperature. The density profiles (n~r^-p ) are estimated from the 3 mm continuum visibility profiles. The masses M and H2 column densities N(H2) are then calculated from the 3 mm dust continuum emission. Results. We find a large spread in mass and peak H2 column density in the detected sources ranging from 0.1-150 Msun and 10^23 - 10^26 cm-2 , respectively. Including the results of the CORE and CORE-extension studies (Gieser et al. 2021, 2022) to increase the sample size, we find evolutionary trends on core scales for the temperature power-~law index q increasing from 0.1 to 0.7 from infrared dark clouds to UCHII regions, while for the the density power-law index p on core scales, we do not find strong evidence for an evolutionary trend. However, we find that on the larger clump scales throughout these evolutionary phases the density profile flattens from p = 2.2 to p = 1.2. (abridged) △ Less

Submitted 24 April, 2023; v1 submitted 14 April, 2023; originally announced April 2023.

Comments: 21 pages, 10 figures, accepted for publication in A&A

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

The ALMA Survey of 70 $μみゅー$m Dark High-mass Clumps in Early Stages (ASHES). IX. Physical Properties and Spatial Distribution of Cores in IRDCs

Authors: Kaho Morii, Patricio Sanhueza, Fumitaka Nakamura, Qizhou Zhang, Giovanni Sabatini, Henrik Beuther, Xing Lu, Shanghuo Li, Guido Garay, James M. Jackson, Fernando A. Olguin, Daniel Tafoya, Ken'ichi Tatematsu, Natsuko Izumi, Takeshi Sakai, Andrea Silva

Abstract: The initial conditions found in infrared dark clouds (IRDCs) provide insights on how high-mass stars and stellar clusters form. We have conducted high-angular resolution and high-sensitivity observations toward thirty-nine massive IRDC clumps, which have been mosaicked using the 12m and 7m arrays from the Atacama Large Millimeter/submillimeter Array (ALMA). The targets are 70 $μみゅー$m dark massive (22… ▽ More The initial conditions found in infrared dark clouds (IRDCs) provide insights on how high-mass stars and stellar clusters form. We have conducted high-angular resolution and high-sensitivity observations toward thirty-nine massive IRDC clumps, which have been mosaicked using the 12m and 7m arrays from the Atacama Large Millimeter/submillimeter Array (ALMA). The targets are 70 $μみゅー$m dark massive (220-4900 $M_\odot$), dense ($>$10$^4$ cm$^{-3}$), and cold ($\sim$10-20K) clumps located at distances between 2 and 6 kpc. We identify an unprecedented number of 839 cores, with masses between 0.05 and 81 $M_\odot$ using 1.3 mm dust continuum emission. About 55% of the cores are low-mass ($<$1 $M_\odot$), whereas $\lesssim$1% (7/839) are high-mass ($\gtrsim$27 $M_\odot$). We detect no high-mass prestellar cores. The most massive cores (MMC) identified within individual clumps lack sufficient mass to form high-mass stars without additional mass feeding. We find that the mass of the MMCs is correlated with the clump surface density, implying denser clumps produce more massive cores and a larger number of cores. There is no significant mass segregation except for a few tentative detections. In contrast, most clumps show segregation once the clump density is considered instead of mass. Although the dust continuum emission resolves clumps in a network of filaments, some of which consist of hub-filament systems, the majority of the MMCs are not found in the hubs. Our analysis shows that high-mass cores and MMCs have no preferred location with respect to low-mass cores at the earliest stages of high-mass star formation. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Comments: Accepted for Publication in ApJ. 54 pages, 40 figures, 5 tables

JOYS: JWST Observations of Young protoStars: Outflows and accretion in the high-mass star-forming region IRAS23385+605

Authors: H. Beuther, E. F. van Dishoeck, L. Tychoniec, C. Gieser, P. J. Kavanagh, G. Perotti, M. L. van Gelder, P. Klaassen, A. Caratti o Garatti, L. Francis, W. R. M. Rocha, K. Slavicinska, T. Ray, K. Justtanont, H. Linnartz, C. Weakens, L. Colina, T. Greve, M. Guedel, T. Henning, P. O. Lagage, B. Vandenbussche, G. Oestlin, G. Wright

Abstract: Aims: The JWST program JOYS (JWST Observations of Young protoStars) aims at characterizing the physical and chemical properties of young high- and low-mass star-forming regions, in particular the unique mid-infrared diagnostics of the warmer gas and solid-state components. We present early results from the high-mass star formation region IRAS23385+6053. Methods: The JOYS program uses the MIRI MRS… ▽ More Aims: The JWST program JOYS (JWST Observations of Young protoStars) aims at characterizing the physical and chemical properties of young high- and low-mass star-forming regions, in particular the unique mid-infrared diagnostics of the warmer gas and solid-state components. We present early results from the high-mass star formation region IRAS23385+6053. Methods: The JOYS program uses the MIRI MRS with its IFU to investigate a sample of high- and low-mass star-forming protostellar systems. Results: The 5 to 28mum MIRI spectrum of IRAS23385+6053 shows a plethora of features. While the general spectrum is typical for an embedded protostar, we see many atomic and molecular gas lines boosted by the higher spectral resolution and sensitivity compared to previous space missions. Furthermore, ice and dust absorption features are also present. Here, we focus on the continuum emission, outflow tracers like the H2, [FeII] and [NeII] lines as well as the potential accretion tracer Humphreys alpha HI(7--6). The short-wavelength MIRI data resolve two continuum sources A and B, where mid-infrared source A is associated with the main mm continuum peak. The combination of mid-infrared and mm data reveals a young cluster in its making. Combining the mid-infrared outflow tracer H2, [FeII] and [NeII] with mm SiO data shows a complex interplay of at least three molecular outflows driven by protostars in the forming cluster. Furthermore, the Humphreys alpha line is detected at a 3-4sigma level towards the mid-infrared sources A and B. Following Rigliaco et al. (2015), one can roughly estimate accretion luminosities and corresponding accretion rates between ~2.6x10^-6 and ~0.9x10^-4 M_sun/yr. This is discussed in the context of the observed outflow rates. Conclusions: The analysis of the MIRI MRS observations for this young high-mass star-forming region reveals connected outflow and accretion signatures. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Comments: 12 pages, 9 figures, accepted for Astronomy & Astrophysics, the paper is also available at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

Journal ref: A&A 673, A121 (2023)

A global view on star formation: The GLOSTAR Galactic plane survey. VII. Supernova remnants in the Galactic longitude range $28^\circ<l<36^\circ$

Authors: R. Dokara, Y. Gong, W. Reich, M. Rugel, A. Brunthaler, K. Menten, W. Cotton, S. Dzib, S. Khan, S. Medina, H. Nguyen, G. Ortiz-León, J. Urquhart, F. Wyrowski, A. Yang, L. D. Anderson, H. Beuther, T. Csengeri, P. Müller, J. Ott, J. D. Pandian, N. Roy

Abstract: Context. While over 1000 supernova remnants (SNRs) are estimated to exist in the Milky Way, only less than 400 have been found to date. In the context of this apparent deficiency, more than 150 SNR candidates were recently identified in the D-configuration Very Large Array (VLA-D) continuum images of the 4--8 GHz global view on star formation (GLOSTAR) survey, in the Galactic longitude range… ▽ More Context. While over 1000 supernova remnants (SNRs) are estimated to exist in the Milky Way, only less than 400 have been found to date. In the context of this apparent deficiency, more than 150 SNR candidates were recently identified in the D-configuration Very Large Array (VLA-D) continuum images of the 4--8 GHz global view on star formation (GLOSTAR) survey, in the Galactic longitude range $-2^\circ<l<60^\circ$. Aims. We attempt to find evidence of nonthermal synchrotron emission from 35 SNR candidates in the region of Galactic longitude range $28^\circ<l<36^\circ$, and also to study the radio continuum emission from the previously confirmed SNRs in this region. Methods. Using the short-spacing corrected GLOSTAR VLA-D+Effelsberg images, we measure ${\sim}6$ GHz total and linearly polarized flux densities of the SNR candidates and the SNRs that were previously confirmed. We also attempt to determine the spectral indices by measuring flux densities from complementary Galactic plane surveys and from the temperature-temperature plots of the GLOSTAR-Effelsberg images. Results. We provide evidence of nonthermal emission from four candidates that have spectral indices and polarization consistent with a SNR origin, and, considering their morphology, we are confident that three of these (G28.36+0.21, G28.78-0.44, and G29.38+0.10) are indeed SNRs. However, about $25\%$ of the candidates have spectral index measurements that indicate thermal emission, and the rest of them are too faint to have a good constraint on the spectral index yet. Conclusions. Additional observations at longer wavelengths and higher sensitivities will shed more light on the nature of these candidates. A simple Monte-Carlo simulation reiterates the view that future studies must persist with the current strategy of searching for SNRs with small angular size to solve the problem of the Milky Way's missing SNRs. △ Less

Submitted 24 November, 2022; originally announced November 2022.

Comments: To be published in A&A. 21 pages, 15 figures

Journal ref: A&A 671, A145 (2023)

Salt-bearing disk candidates around high-mass young stellar objects

Authors: Adam Ginsburg, Brett A. McGuire, Patricio Sanhueza, Fernando Olguin, Luke T Maud, Kei E. I. Tanaka, Yichen Zhang, Henrik Beuther, Nick Indriolo

Abstract: Molecular lines tracing the orbital motion of gas in a well-defined disk are valuable tools for inferring both the properties of the disk and the star it surrounds. Lines that arise only from a disk, and not also from the surrounding molecular cloud core that birthed the star or from the outflow it drives, are rare. Several such emission lines have recently been discovered in one example case, tho… ▽ More Molecular lines tracing the orbital motion of gas in a well-defined disk are valuable tools for inferring both the properties of the disk and the star it surrounds. Lines that arise only from a disk, and not also from the surrounding molecular cloud core that birthed the star or from the outflow it drives, are rare. Several such emission lines have recently been discovered in one example case, those from NaCl and KCl salt molecules. We studied a sample of 23 candidate high-mass young stellar objects (HMYSOs) in 17 high-mass star-forming regions to determine how frequently emission from these species is detected. We present five new detections of water, NaCl, KCl, PN, and SiS from the innermost regions around the objects, bringing the total number of known briny disk candidates to nine. Their kinematic structure is generally disk-like, though we are unable to determine whether they arise from a disk or outflow in the sources with new detections. We demonstrate that these species are spatially coincident in a few resolved cases and show that they are generally detected together, suggesting a common origin or excitation mechanism. We also show that several disks around HMYSOs clearly do not exhibit emission in these species. Salty disks are therefore neither particularly rare in high-mass disks, nor are they ubiquitous. △ Less

Submitted 4 November, 2022; originally announced November 2022.

Comments: accepted to ApJ

Properties of atomic hydrogen gas in the Galactic plane from THOR 21-cm absorption spectra: a comparison with the high latitude gas

Authors: Arghyadeep Basu, Nirupam Roy, Henrik Beuther, Jonas Syed, Jürgen Ott, Juan D. Soler, Jeroen Stil, Michael R. Rugel

Abstract: The neutral hydrogen 21 cm line is an excellent tracer of the atomic interstellar medium in the cold and the warm phases. Combined 21 cm emission and absorption observations are very useful to study the properties of the gas over a wide range of density and temperature. In this work, we have used 21 cm absorption spectra from recent interferometric surveys, along with the corresponding emission sp… ▽ More The neutral hydrogen 21 cm line is an excellent tracer of the atomic interstellar medium in the cold and the warm phases. Combined 21 cm emission and absorption observations are very useful to study the properties of the gas over a wide range of density and temperature. In this work, we have used 21 cm absorption spectra from recent interferometric surveys, along with the corresponding emission spectra from earlier single dish surveys to study the properties of the atomic gas in the Milky Way. In particular, we focus on a comparison of properties between lines of sight through the gas disk in the Galactic plane and high Galactic latitude lines of sight through more diffuse gas. As expected, the analysis shows a lower average temperature for the gas in the Galactic plane compared to that along the high latitude lines of sight. The gas in the plane also has a higher molecular fraction, showing a sharp transition and flattening in the dust - gas correlation. On the other hand, the observed correlation between 21 cm brightness temperature and optical depth indicates some intrinsic difference in spin temperature distribution and a fraction of gas in the Galactic plane having intermediate optical depth (for 0.02 < $τたう$ < 0.2) but higher spin temperature, compared to that of the diffuse gas at high latitude with the same optical depth. This may be due to a small fraction of cold gas with slightly higher temperature and lower density present on the Galactic plane. △ Less

Submitted 20 October, 2022; originally announced October 2022.

Comments: 6 pages, 5 figures, submitted and accepted for publication in MNRAS

A global view on star formation: The GLOSTAR Galactic plane survey. VI. Radio Source Catalog II: $28^\circ < \ell < 36^\circ$ and $|b| < 1^\circ$, VLA B-configuration

Authors: S. A. Dzib, A. Y. Yang, J. S. Urquhart, S. -N. X. Medina, A. Brunthaler, K. M. Menten, F. Wyrowski, W. D. Cotton, R. Dokara, G. N. Ortiz-León, M. R. Rugel, H. Nguyen, Y. Gong, A. Chakraborty, H. Beuther, S. J. Billington, C. Carrasco-Gonzalez, T. Csengeri, P. Hofner, J. Ott, J. D. Pandian, N. Roy, V. Yanza

Abstract: As part of the GLOSTAR survey we have used the VLA in its B-configuration to observe the part of the Galactic plane between longitudes of 28d and 36d and latitudes from -1d to +1d at the C-band (4--8 GHz). To reduce the contamination of extended sources that are not well recovered by our coverage of the (u, v)-plane we discarded short baselines that are sensitive to emission on angular scales… ▽ More As part of the GLOSTAR survey we have used the VLA in its B-configuration to observe the part of the Galactic plane between longitudes of 28d and 36d and latitudes from -1d to +1d at the C-band (4--8 GHz). To reduce the contamination of extended sources that are not well recovered by our coverage of the (u, v)-plane we discarded short baselines that are sensitive to emission on angular scales $<4"$. The resulting radio continuum images have an angular resolution of 1.0", and sensitivity of $\sim60 μみゅー$Jy~beam$^{-1}$; making it the most sensitive radio survey covering a large area of the Galactic plane with this angular resolution. An automatic source extraction algorithm was used in combination with visual inspection to identify a total of 3325 radio sources. A total of 1457 radio sources are $\geq7σしぐま$ and comprise our highly reliable catalog; 72 of these are grouped as 22 fragmented sources, e.g., multiple components of an extended and resolved source. To explore the nature of the catalogued radio sources we searched for counterparts at millimeter and infrared wavelengths. Our classification attempts resulted in 93 HII region candidates, 104 radio stars, 64 planetary nebulae, while most of the remaining radio sources are suggested to be extragalactic sources. We investigated the spectral indices ($αあるふぁ$, $S_νにゅー\proptoνにゅー^αあるふぁ$) of radio sources classified as HII region candidates and found that many have negative values. This may imply that these radio sources represent young stellar objects that are members of the star clusters around the high mass stars that excite the HII regions, but not these HII regions themselves. By comparing the peak flux densities from the GLOSTAR and CORNISH surveys we have identified 49 variable radio sources, most of them with an unknown nature. Additionally, we provide the list of 1866 radio sources detected within 5 to 7$σしぐま$ levels. △ Less

Submitted 6 October, 2022; v1 submitted 2 October, 2022; originally announced October 2022.

Comments: 25 pages, 15 figures, and 6 tables, Accepted for publication in the Astronomy \& Astrophysics

Journal ref: A&A 670, A9 (2023)

Polarized Emission From Four Supernova Remnants In The THOR Survey

Authors: Russell Shanahan, Jeroen M. Stil, Loren Anderson, Henrik Beuther, Paul Goldsmith, Juergen Ott, Michael Rugel, Juan Soler, Jonas Syed

Abstract: We present polarization and Faraday rotation for the supernova remnants (SNRs) G46.8-0.3, G43.3-0.2, G41.1-0.3, and G39.2-0.3 in L-band (1-2 GHz) radio continuum in The HI/OH/Recombination line (THOR) survey. We detect polarization from G46.8-0.3, G43.3-0.2 and G39.2-0.3 but find upper limits at the 1% level of Stokes I for G41.1-0.3. For G46.8-0.3 and G39.2-0.3 the fractional polarization varies… ▽ More We present polarization and Faraday rotation for the supernova remnants (SNRs) G46.8-0.3, G43.3-0.2, G41.1-0.3, and G39.2-0.3 in L-band (1-2 GHz) radio continuum in The HI/OH/Recombination line (THOR) survey. We detect polarization from G46.8-0.3, G43.3-0.2 and G39.2-0.3 but find upper limits at the 1% level of Stokes I for G41.1-0.3. For G46.8-0.3 and G39.2-0.3 the fractional polarization varies on small scales from 1% to ~6%. G43.3-0.2 is less polarized with fractional polarization <~3%. We find upper limits at the 1% level for the brighter regions in each SNR with no evidence for associated enhanced Faraday depolarization. We observe significant variation in Faraday depth and fractional polarization on angular scales down to the resolution limit of 16". Approximately 6% of our polarization detections from G46.8-0.3 and G39.2-0.3 exhibit two-component Faraday rotation and 14% of polarization detections in G43.3-0.2 are multi-component. For G39.2-0.3 we find a bimodal Faraday depth distribution with a narrow peak and a broad peak for all polarization detections as well as for the subset with two-component Faraday rotation. We identify the narrow peak with the front side of the SNR and the broad peak with the back side. Similarly, we interpret the observed Faraday depth distribution of G46.8-0.3 as a superposition of the distributions from the front side and the back side. We interpret our results as evidence for a partially filled shell with small-scale magnetic field structure and internal Faraday rotation. △ Less

Submitted 27 September, 2022; originally announced September 2022.

Comments: Accepted for publication in the Astrophysical Journal. Figures 6, 19 and 20 may not be displayed in the browser pdf viewer, but the downloaded pdf is complete

Snapshot of a magnetohydrodynamic disk wind traced by water maser observations

Authors: Luca Moscadelli, Alberto Sanna, Henrik Beuther, André Oliva, Rolf Kuiper

Abstract: The formation of astrophysical objects of different nature and size, from black holes to gaseous giant planets, involves a disk-jet system, where the disk drives the mass accretion onto a central compact object and the jet is a fast collimated ejection along the disk rotation axis. Magnetohydrodynamic disk winds can provide the link between mass accretion and ejection, which is essential to ensure… ▽ More The formation of astrophysical objects of different nature and size, from black holes to gaseous giant planets, involves a disk-jet system, where the disk drives the mass accretion onto a central compact object and the jet is a fast collimated ejection along the disk rotation axis. Magnetohydrodynamic disk winds can provide the link between mass accretion and ejection, which is essential to ensure that the excess angular momentum is removed from the system and accretion onto the central object can proceed. However, up to now, we have been lacking direct observational proof of disk winds. This work presents a direct view of the velocity field of a disk wind around a forming massive star. Achieving a very high spatial resolution of ~0.05 au, our water maser observations trace the velocities of individual streamlines emerging from the disk orbiting the forming star. We find that, at low elevation above the disk midplane, the flow co-rotates with its launch point in the disk, in agreement with magneto-centrifugal acceleration where the gas is flung away along the magnetic field line anchored to the disk. Beyond the co-rotation point, the flow rises spiraling around the disk rotation axis along a helical magnetic field. We have performed (resistive-radiative-gravito-) magnetohydrodynamic simulations of the formation of a massive star and record the development of a magneto-centrifugally launched jet presenting many properties in agreement with our observations. △ Less

Submitted 1 September, 2022; originally announced September 2022.

Comments: Main paper: 15 pages, 6 figures, 2 tables; Appendix: 21 pages, 7 figures, 1 table

Journal ref: published on-line in Nature Astronomy on August 11th 2022; view-only at: https://rdcu.be/cTukQ ; You Tube video: https://www.youtube.com/watch?v=EcmmvuVWa3Q

The Cygnus Allscale Survey of Chemistry and Dynamical Environments: CASCADE: Overview and first results toward DR20 from the Max Planck IRAM Observatory program (MIOP)

Authors: H. Beuther, F. Wyrowski, K. M. Menten, J. M. Winters, S. Suri, W. -J. Kim, L. Bouscasse, C. Gieser, M. Sawczuck, I. B. Christensen, I. M. Skretas

Abstract: Context: While star formation on large molecular cloud scales and on small core and disk scales has been investigated intensely over the past decades, the connection of the large-scale interstellar material with the densest small-scale cores has been a largely neglected field. Methods: Using NOEMA and the IRAM 30\,m telescope, we mapped large areas (640\,arcmin$^2$) of the archetypical star form… ▽ More Context: While star formation on large molecular cloud scales and on small core and disk scales has been investigated intensely over the past decades, the connection of the large-scale interstellar material with the densest small-scale cores has been a largely neglected field. Methods: Using NOEMA and the IRAM 30\,m telescope, we mapped large areas (640\,arcmin$^2$) of the archetypical star formation complex Cygnus X at 3.6\,mm wavelengths in line and continuum emission. Results: The scope and outline of The Cygnus Allscale Survey of Chemistry and Dynamical Environments (CASCADE) is presented. We then focus on the first observed subregion in Cygnus X, namely the DR20 star formation site, which comprises sources in a range of evolutionary stages from cold pristine gas clumps to more evolved ultracompact H{\sc ii} regions. The data covering cloud to cores scales at a linear spatial resolution of $<5000$\,au reveal several kinematic cloud components that are likely part of several large-scale flows onto the central cores. The temperature structure of the region is investigated by means of the HCN/HNC intensity ratio and compared to dust-derived temperatures. We find that the deuterated DCO$^+$ emission is almost exclusively located toward regions at low temperatures below 20\,K. Investigating the slopes of spatial power spectra of dense gas tracer intensity distributions (HCO$^+$, H$^{13}$CO$^+$, and N$_2$H$^+$), we find comparatively flat slopes between $-2.9$ and $-2.6$, consistent with high Mach numbers and/or active star formation in DR20. △ Less

Submitted 22 July, 2022; originally announced July 2022.

Comments: 18 pages, 16 figures, accepted to A&A, also available at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

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

A Global View on Star Formation: The GLOSTAR Galactic Plane Survey V. 6.7 GHz Methanol Maser Catalogue

Authors: H. Nguyen, M. R. Rugel, C. Murugeshan, K. M. Menten, A. Brunthaler, J. S. Urquhart, R. Dokara, S. A. Dzib, Y. Gong, S. Khan, S-N. X. Medina, G. N. Ortiz-Leon, W. Reich, F. Wyrowski, A. Y. Yang, H. Beuther, W. D. Cotton, J. D. Pandian

Abstract: Class II methanol (CH$_{3}$OH) masers are amongst the clearest signposts of recent high-mass star formation (HMSF). A complete catalogue outlines the distribution of star formation in the Galaxy, the number of young star-forming cores, and the physical conditions of their environment. The Global View on Star Formation (GLOSTAR) survey, which is a blind survey in the radio regime of 4$-$8 GHz, maps… ▽ More Class II methanol (CH$_{3}$OH) masers are amongst the clearest signposts of recent high-mass star formation (HMSF). A complete catalogue outlines the distribution of star formation in the Galaxy, the number of young star-forming cores, and the physical conditions of their environment. The Global View on Star Formation (GLOSTAR) survey, which is a blind survey in the radio regime of 4$-$8 GHz, maps the Galactic mid-plane in the radio continuum, 6.7 GHz methanol line, the 4.8 GHz formaldehyde line, and several radio recombination lines. We present the analysis of the observations of the 6.7 GHz CH$_{3}$OH maser transition using data from the D-configuration of the Very Large Array (VLA). We analyse the data covering Galactic longitudes from $-2^{\circ}< l <60^{\circ}$ and Galactic latitudes of $|\textit{b}|<1^{\circ}$. We detect a total of 554 methanol masers, out of which 84 are new, and catalogue their positions, velocity components, and integrated fluxes. With a typical noise level of $\sim$18 mJy beam$^{-1}$, this is the most sensitive unbiased methanol survey for methanol masers to date. We search for dust continuum and radio continuum associations, and find that 97% of the sources are associated with dust, and 12% are associated with radio continuum emission. △ Less

Submitted 21 July, 2022; originally announced July 2022.

Comments: Accepted in A&A July 18, 2022

Journal ref: A&A 666, A59 (2022)

The Galactic dynamics revealed by the filamentary structure in atomic hydrogen emission

Authors: Juan D. Soler, Marc-Antoine Miville-Deschênes, Sergio Molinari, Ralf S. Klessen, Patrick Hennebelle, Leonardo Testi, Naomi M. McClure-Griffiths, Henrik Beuther, Davide Elia, Eugenio Schisano, Alessio Traficante, Philipp Girichidis, Simon C. O. Glover, Rowan J. Smith, Mattia Sormani, Robin Treß

Abstract: We present a study of the filamentary structure in the atomic hydrogen (HI) emission at the 21 cm wavelength toward the Galactic plane using the observations in the HI4PI survey. Using the Hessian matrix method across radial velocity channels, we identified the filamentary structures and quantified their orientations using circular statistics. We found that the regions of the Milky Way's disk beyo… ▽ More We present a study of the filamentary structure in the atomic hydrogen (HI) emission at the 21 cm wavelength toward the Galactic plane using the observations in the HI4PI survey. Using the Hessian matrix method across radial velocity channels, we identified the filamentary structures and quantified their orientations using circular statistics. We found that the regions of the Milky Way's disk beyond 10 kpc and up to roughly 18 kpc from the Galactic center display HI filamentary structures predominantly parallel to the Galactic plane. For regions at lower Galactocentric radii, we found that the HI filaments are mostly perpendicular or do not have a preferred orientation with respect to the Galactic plane. We interpret these results as the imprint of supernova feedback in the inner Galaxy and Galactic rotation in the outer Milky Way. We found that the HI filamentary structures follow the Galactic warp and that they highlight some of the variations interpreted as the effect of the gravitational interaction with satellite galaxies. In addition, the mean scale height of the filamentary structures is lower than that sampled by the bulk of the HI emission, thus indicating that the cold and warm atomic hydrogen phases have different scale heights in the outer galaxy. Finally, we found that the fraction of the column density in HI filaments is almost constant up to approximately 18 kpc from the Galactic center. This is possibly a result of the roughly constant ratio between the cold and warm atomic hydrogen phases inferred from the HI absorption studies. Our results indicate that the HI filamentary structures provide insight into the dynamical processes shaping the Galactic disk. Their orientations record how and where the stellar energy input, the Galactic fountain process, the cosmic ray diffusion, and the gas accretion have molded the diffuse interstellar medium in the Galactic plane. △ Less

Submitted 20 May, 2022; originally announced May 2022.

Comments: 33 pages. 39 figures. Accepted for publication in Astronomy & Astrophysics (06MAY2022)

Report number: AA/2022/43334

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

The sharp ALMA view of infall and outflow in the massive protocluster G31.41+0.31

Authors: M. T. Beltrán, V. M. Rivilla, R. Cesaroni, D. Galli, L. Moscadelli, A. Ahmadi, H. Beuther, S. Etoka, C. Goddi, P. D. Klaassen, R. Kuiper, M. S. N. Kumar, A. Lorenzani, T. Peters, Á. Sánchez-Monge, P. Schilke, F. van der Tak, S. Vig

Abstract: Context. To better understand the formation of high-mass stars, it is fundamental to investigate how matter accretes onto young massive stars, how it is ejected, and how all this differs from the low-mass case. The massive protocluster G31.41+0.31 is the ideal target to study all these processes because observations at millimeter and centimeter wavelengths have resolved the emission of the Main co… ▽ More Context. To better understand the formation of high-mass stars, it is fundamental to investigate how matter accretes onto young massive stars, how it is ejected, and how all this differs from the low-mass case. The massive protocluster G31.41+0.31 is the ideal target to study all these processes because observations at millimeter and centimeter wavelengths have resolved the emission of the Main core into at least four massive dust continuum sources, named A, B, C, and D, within 1" or 0.018 pc, and have identified signatures of infall and several outflows associated with the core. Aims. We study the interplay between infall and outflow in G31.41+0.31 by investigating at a spatial resolution of a few 100 au their properties and their possible impact on the core. Methods. We carried out molecular line observations of typical high-density tracers, such as CH3CN or H2CO, and shock and outflow tracers, such as SiO, with ALMA at 1.4 mm that achieved an angular resolution of 0.09" (340 au). Results. The observations have revealed inverse P-Cygni profiles in CH3CN and H2CO toward the four sources embedded in the Main core, suggesting that all of them are undergoing collapse. The infall rates, estimated from the red-shifted absorption are on the order of 1E-2 Msun/yr. The individual infall rates imply that the accretion timescale of the Main core is an order of magnitude smaller than its rotation timescale. This confirms that rotating toroids such as the G31 Main core are non-equilibrium, transient collapsing structures that need to be constantly replenished with fresh material from a large-scale reservoir. For sources B, C, and D, the infall could be accelerating inside the sources, while for source A, the presence of a second emission component complicates the interpretation. The SiO observations have revealed the presence of at least six outflows in the G31.41+0.31 star-forming region, ... △ Less

Submitted 25 January, 2022; originally announced January 2022.

Comments: 17 pages, 9 figures, 2 tables. Accepted by A&A

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

FEEDBACK from the NGC7538 HII region

Authors: H. Beuther, N. Schneider, R. Simon, S. Suri, V. Ossenkopf-Okada, S. Kabanovic, M. Roellig, C. Guevara, A. G. G. M. Tielens, G. Sandell, C. Buchbender, O. Ricken, R. Guesten

Abstract: Context: How do expanding HII regions interact with their environmental cloud? This is one of the central questions driving the SOFIA legacy program FEEDBACK. Here, we present a case study toward the prototypical H{\sc ii} region NGC7538. Methods: With SOFIA we mapped an area of ~210'^2 around NGC7538 in the [CII] line at 1.9THz. Complementary observed atomic carbon [CI] and high-J CO(8-7) data as… ▽ More Context: How do expanding HII regions interact with their environmental cloud? This is one of the central questions driving the SOFIA legacy program FEEDBACK. Here, we present a case study toward the prototypical H{\sc ii} region NGC7538. Methods: With SOFIA we mapped an area of ~210'^2 around NGC7538 in the [CII] line at 1.9THz. Complementary observed atomic carbon [CI] and high-J CO(8-7) data as well as archival NIR/FIR, cm continuum, CO(3-2) and HI data are folded into the analysis. Results: While the overall [CII] morphology follows the general ionized gas, the channel maps show multiple bubble-like structures with sizes on the order of ~80-100" (~1.0-1.28pc). While at least one of them may be an individual feedback bubble driven by the main exciting sources of the region, the other bubble-morphologies may also be due to the intrinsically porous structure of the HII region. An analysis of the expansion velocities around 10km s^{-1} indicates that thermal expansion is not sufficient but that wind-driving from the central O-stars is required. The most blue-shifted [CII] component has barely any molecular or atomic counterparts. At the interface to the molecular cloud, we find a typical photon-dominated region (PDR) with a bar-shape. Ionized, atomic and molecular carbon show a layered structure in this PDR. The carbon in the PDR is dominated by its ionized form with atomic and molecular masses of ~0.45+-0.1M_{\odot} and ~1.2+-0.1M_{\odot}, respectively, compared to the ionized carbon in the range of 3.6-9.7M_{\odot}. Conclusions: The NGC7538 HII region exhibits a diverse set of sub-structures that interact with each other as well as with the adjacent cloud. Compared to other recent [CII] observations of HII regions (e.g., Orion Veil, RCW120, RCW49), bubble-shape morphologies revealed in [CII] emission, indicative of expanding shells, are recurring structures of PDRs. △ Less

Submitted 14 January, 2022; originally announced January 2022.

Comments: 21 pages, 20 figures, accepted for Astronomy & Astrophysics, a high-resolution version can be found at https://www2.mpia-hd.mpg.de/homes/beuther/papers.html

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

ALMA observations of NGC 6334S. II. Subsonic and Transonic Narrow Filaments in a High-mass Star Formation Cloud

Authors: Shanghuo Li, Patricio Sanhueza, Chang Won Lee, Qizhou Zhang, Henrik Beuther, Aina Palau, Hong-Li Liu, Howard Smith, Hauyu Baobab Liu, Izaskun, Jiménez-Serra, Kee-Tae Kim, Siyi Feng, Josep Miquel. Girart, Tie Liu, Junzhi Wang, Di Li, Keping Qiu, Xing Lu, Ke Wang, Fei Li, Juan Li, Yue Cao, Shinyoung Kim, Shaye Strom

Abstract: We present a study of narrow filaments toward a massive infrared dark cloud, NGC 6334S, using the Atacama Large Millimeter/submillimeter Array (ALMA). Thirteen gas filaments are identified using the H$^{13}$CO$^{+}$ line, while a single continuum filament is revealed by the continuum emission. The filaments present a compact radial distribution with a median filament width of $\sim$0.04 pc narrowe… ▽ More We present a study of narrow filaments toward a massive infrared dark cloud, NGC 6334S, using the Atacama Large Millimeter/submillimeter Array (ALMA). Thirteen gas filaments are identified using the H$^{13}$CO$^{+}$ line, while a single continuum filament is revealed by the continuum emission. The filaments present a compact radial distribution with a median filament width of $\sim$0.04 pc narrower than the previously proposed `quasi-universal' 0.1~pc filament width. The higher spatial resolution observations and higher-density gas tracer tend to identify even narrower and lower mass filaments. The filament widths are roughly twice the size of embedded cores. The gas filaments are largely supported by thermal motions. The nonthermal motions are predominantly subsonic and transonic in both identified gas filaments and embedded cores, which may imply that stars are likely born in environments of low turbulence. A fraction of embedded objects show a narrower velocity dispersion compared with their corresponding natal filaments, which may indicate that the turbulent dissipation is taking place in these embedded cores. The physical properties (mass, mass per unit length, gas kinematics, and width) of gas filaments are analogous to those of narrow filaments found in low- to high-mass star-forming regions. The more evolved sources are found to be farther away from the filaments, a situation that may have resulted from the relative motions between the YSOs and their natal filaments. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Comments: 28 pages, 12 figures, 1 table. Accepted for publication in ApJ