Massive Star-forming Regions Research Articles

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions – IV. Radio recombination lines and evolution of star formation efficiencies

ABSTRACT We report the detection of radio recombination line (RRL) H40α towards 75 sources, with data obtained from ACA (Atacama Compact 7 m Array) observations in the ATOMS (ALMA Three-millimeter Observations of Massive Star-forming regions) survey of 146 active Galactic star-forming regions. We calculated ionized gas mass and star formation rate (SFR) with H40α line emission. The mass of ionized gas is significantly smaller than molecular gas mass, indicating that ionized gas is negligible in the star-forming clumps of the ATOMS sample. The SFR estimated with RRL H40α agrees well with that calculated with the total bolometric luminosity (Lbol) when SFR ≥ 5 M⊙ Myr−1, suggesting that millimetre RRLs could well sample the upper part of the initial mass function and thus be good tracers of SFR. We also study the relationships between Lbol and the molecular line luminosities (L$^{^{\prime }}_{\rm mol}$) of CS J = 2–1 and HC3N J = 11–10 for all the 146 ATOMS sources. The $L_{\rm bol} \!-\! L^{\prime }_{\rm mol}$ correlations of both the CS J = 2–1 and HC3N J = 11–10 lines appear approximately linear and these transitions have success in predicting Lbol similar to that of more commonly used transitions. The Lbol-to-$L^{\prime }_{\rm mol}$ ratios or SFR-to-mass ratios (star formation efficiency) do not change with galactocentric distances (RGC). Sources with H40α emission (or H ii regions) show higher Lbol-to-$L^{\prime }_{\rm mol}$ ratios than those without H40α emission, which may be an evolutionary effect.

High-spatial resolution long-mid-infrared observations of massive star-forming regions: Local star formation sequence

Abstract The formation of massive stars in dense and cold molecular clouds is a fundamental problem in star formation. In this work, we studied three compact massive star-forming regions, M8E, RAFGL 6366S, and IRAS 18317−0513, to investigate the environment of massive star formation, specifically the possibility of a local feedback process. Our new mid-infrared observations using miniTAO/MAX38 at 31 and 37 μm resolved individual objects in each region and allowed us to obtain their individual luminosities and masses. Together with existence/absence of ultra-compact H ii regions, it is suggested that less-massive objects are more evolved than the more-massive objects in two out of the three regions. Because the objects that are more massive evolve faster, those that are less massive form earlier. This formation trend of local mass-sequential star was first suggested in 0.1 pc-scale massive star-forming regions. Therefore, feedback, such as outflows, jets, or radiative heating, from previously formed lower-mass young stellar objects might affect the environment of parental clouds and lead to next-generation massive young stellar objects.

Kinetic temperature of massive star-forming molecular clumps measured with formaldehyde

We mapped the kinetic temperature structure of two massive star-forming regions, N113 and N159W, in the Large Magellanic Cloud (LMC). We have used ~1.′′6 (~0.4 pc) resolution measurements of the para-H2CO JKaKc = 303–202, 322–221, and 321–220 transitions near 218.5 GHz to constrain RADEX non local thermodynamic equilibrium models of the physical conditions. The gas kinetic temperatures derived from the para-H2CO line ratios 322–221/303–202 and 321–220/303–202 range from 28 to 105 K in N113 and 29 to 68 K in N159W. Distributions of the dense gas traced by para-H2CO agree with those of the 1.3 mm dust and Spitzer 8.0 μm emission, but they do not significantly correlate with the Hα emission. The high kinetic temperatures (Tkin ≳ 50 K) of the dense gas traced by para-H2CO appear to be correlated with the embedded infrared sources inside the clouds and/or young stellar objects in the N113 and N159W regions. The lower temperatures (Tkin < 50 K) were measured at the outskirts of the H2CO-bearing distributions of both N113 and N159W. It seems that the kinetic temperatures of the dense gas traced by para-H2CO are weakly affected by the external sources of the Hα emission. The non thermal velocity dispersions of para-H2CO are well correlated with the gas kinetic temperatures in the N113 region, implying that the higher kinetic temperature traced by para-H2CO is related to turbulence on a ~0.4 pc scale. The dense gas heating appears to be dominated by internal star formation activity, radiation, and/or turbulence. It seems that the mechanism heating the dense gas of the star-forming regions in the LMC is consistent with that in Galactic massive star-forming regions located in the Galactic plane.

Magnetic Fields in the Massive Star-forming Region GL 437

Abstract In this note, I report a R-band polarimetric study on a massive star-forming region GL 437 (a.k.a. AFGL 437) to infer the magnetic field (B-field) geometry of this region. The mean B-field in the diffused region, is found to be parallel to the extended cloud structure and the Galactic plane. I compared present optical polarization results with archival submm results adapted from Curran & Chrysostomou. A flipped field geometry is seen in denser parts of GL 437 which is an interesting finding. An infrared source WK 34 embedded in the cloud is associated with a poorly collimated outflow. These outflows are parallel to the inner field orientation but misaligned with outer B-fields. An investigation of gas kinematics of this region will be carried out to understand if the gas motion have some role in changing the field geometry.

ATOMS: ALMA Three-millimeter Observations of Massive Star-forming regions–VI. On the formation of the ‘L’ type filament in G286.21+0.17

ABSTRACT Filaments play an important role in star formation, but the formation process of filaments themselves is still unclear. The high-mass star-forming clump G286.21+0.17 (G286 for short) that contains an ‘L’ type filament was thought to undergo global collapse. Our high-resolution ALMA band 3 observations resolve the gas kinematics of G286 and reveal two sub-clumps with very different velocities inside it. We find that the ‘blue profile’ (an indicator of gas infall) of HCO+ lines in single dish observations of G286 is actually caused by gas emission from the two sub-clumps rather than gas infall. We advise great caution in interpreting gas kinematics (e.g. infall) from line profiles towards distant massive clumps in single dish observations. Energetic outflows are identified in G286 but the outflows are not strong enough to drive expansion of the two sub-clumps. The two parts of the ‘L’ type filament (‘NW–SE’ and ‘NE–SW’ filaments) show prominent velocity gradients perpendicular to their major axes, indicating that they are likely formed due to large-scale compression flows. We argue that the large-scale compression flows could be induced by the expansion of nearby giant H ii regions. The ‘NW–SE’ and ‘NE–SW’ filaments seem to be in collision, and a large amount of gas has been accumulated in the junction region where the most massive core G286c1 forms.

Evolutions of CH3CN abundance in molecular clumps

To investigate the effects of massive star evolution on surrounding molecules, we select nine massive clumps previously observed with the Atacama Pathfinder Experiment (APEX) telescope and the Submillimeter Array (SMA) telescope. Based on the observations of APEX, we obtain luminosity to mass ratios L clump/M clump that range from 10 to 154 L ⊙/M ⊙, where some of them embedded ultra compact (UC) HII region. Using the SMA, CH3CN (12K–11K) transitions were observed toward nine massive star-forming regions. We derive the CH3CN rotational temperature and column density using the XCLASS program, and calculate its fractional abundance. We find that CH3CN temperature seems to increase with the increase of L clump/M clump when the ratio is between 10 to 40 L ⊙/M ⊙, then decrease when L clump/M clump ≥ 40 L ⊙/M ⊙. Assuming that the CH3CN gas is heated by radiation from the central star, the effective distance of CH3CN relative to the central star is estimated. The distance ranges from ∼ 0.003 to ∼ 0.083 pc, which accounts for ∼ 1/100 to ∼ 1/1000 of clump size. The effective distance increases slightly as L clump/M clump increases (R eff ∼ (L clump/M clump)0.5±0.2). Overall, the CH3CN abundance is found to decrease as the clumps evolve, e.g., X CH3CN ∼ (L clump/M clump)−1.0 ± 0.7. The steady decline of CH3CN abundance as the clumps evolution can be interpreted as a result of photodissociation.

Calibrating the Davis–Chandrasekhar–Fermi Method with Numerical Simulations: Uncertainties in Estimating the Magnetic Field Strength from Statistics of Field Orientations

Abstract The Davis–Chandrasekhar–Fermi (DCF) method is widely used to indirectly estimate the magnetic field strength from the plane-of-sky field orientation. In this work, we present a set of 3D MHD simulations and synthetic polarization images using radiative transfer of clustered massive star-forming regions. We apply the DCF method to the synthetic polarization maps to investigate its reliability in high-density molecular clumps and dense cores where self-gravity is significant. We investigate the validity of the assumptions of the DCF method step by step and compare the model and estimated field strength to derive the correction factors for the estimated uniform and total (rms) magnetic field strength at clump and core scales. The correction factors in different situations are cataloged. We find the DCF method works well in strong-field cases. However, the magnetic field strength in weak-field cases could be significantly overestimated by the DCF method when the turbulent magnetic energy is smaller than the turbulent kinetic energy. We investigate the accuracy of the angular dispersion function (ADF; a modified DCF method) method on the effects that may affect the measured angular dispersion and find that the ADF method correctly accounts for the ordered field structure, the beam smoothing, and the interferometric filtering, but may not be applicable to account for the signal integration along the line of sight in most cases. Our results suggest that the DCF methods should be avoided to be applied below ∼0.1 pc scales if the effect of line-of-sight signal integration is not properly addressed.

Singly and doubly deuterated formaldehyde in massive star-forming regions

Context. Deuterated molecules are good tracers of the evolutionary stage of star-forming cores. During the star formation process, deuterated molecules are expected to be enhanced in cold, dense pre-stellar cores and to deplete after protostellar birth. Aims. In this paper, we study the deuteration fraction of formaldehyde in high-mass star-forming cores at different evolutionary stages to investigate whether the deuteration fraction of formaldehyde can be used as an evolutionary tracer. Methods. Using the APEX SEPIA Band 5 receiver, we extended our pilot study of the J = 3 →2 rotational lines of HDCO and D2CO to eleven high-mass star-forming regions that host objects at different evolutionary stages. High-resolution follow-up observations of eight objects in ALMA Band 6 were performed to reveal the size of the H2CO emission and to give an estimate of the deuteration fractions HDCO/H2CO and D2CO/HDCO at scales of ~6″ (0.04–0.15 pc at the distance of our targets). Results. Our observations show that singly and doubly deuterated H2CO are detected towards high-mass protostellar objects (HMPOs) and ultracompact H II regions (UC H II regions), and the deuteration fraction of H2CO is also found to decrease by an order of magnitude from the earlier HMPO phases to the latest evolutionary stage (UC H II), from ~0.13 to ~0.01. We have not detected HDCO and D2CO emission from the youngest sources (i.e. high-mass starless cores or HMSCs). Conclusions. Our extended study supports the results of the previous pilot study: the deuteration fraction of formaldehyde decreases with the evolutionary stage, but higher sensitivity observations are needed to provide more stringent constraints on the D/H ratio during the HMSC phase. The calculated upper limits for the HMSC sources are high, so the trend between HMSC and HMPO phases cannot be constrained.

The JCMT BISTRO Survey: An 850/450 μm Polarization Study of NGC 2071IR in Orion B

Abstract We present the results of simultaneous 450 μm and 850 μm polarization observations toward the massive star-forming region NGC 2071IR, a target of the BISTRO (B-fields in STar-forming Region Observations) Survey, using the POL-2 polarimeter and SCUBA-2 camera mounted on the James Clerk Maxwell Telescope. We find a pinched magnetic field morphology in the central dense core region, which could be due to a rotating toroidal disklike structure and a bipolar outflow originating from the central young stellar object IRS 3. Using the modified Davis–Chandrasekhar–Fermi method, we obtain a plane-of-sky magnetic field strength of 563 ± 421 μG in the central ∼0.12 pc region from 850 μm polarization data. The corresponding magnetic energy density of 2.04 × 10−8 erg cm−3 is comparable to the turbulent and gravitational energy densities in the region. We find that the magnetic field direction is very well aligned with the whole of the IRS 3 bipolar outflow structure. We find that the median value of polarization fractions is 3.0% at 450 μm in the central 3′ region, which is larger than the median value of 1.2% at 850 μm. The trend could be due to the better alignment of warmer dust in the strong radiation environment. We also find that polarization fractions decrease with intensity at both wavelengths, with slopes, determined by fitting a Rician noise model of 0.59 ± 0.03 at 450 μm and 0.36 ± 0.04 at 850 μm, respectively. We think that the shallow slope at 850 μm is due to grain alignment at the center being assisted by strong radiation from the central young stellar objects.

Radiatively Cooling Superwinds in Ultracompact Hii Regions

AbstractUltracompact Hii regions (UC-HII) are the young, very dense cores of massive star-forming regions in dwarf galaxies, where newly formed massive OB stars are surrounded by natal molecular clouds. Thermal energy deposited by mechanical feedback from a cluster of massive OB stars can form a superwind, which may lead to a wind-blown bubble as well as radiative cooling. We investigate the formation of radiatively cooling superwinds in UC-HII using a radiative cooling module in the hydrodynamics program. We built a grid of hydrodynamic simulations to determine the dependence of radiative cooling on the cluster radius, mass-deposition rate, wind velocity, and ambient medium in UC-HII. Our findings could help to better understand star formation in massive star-forming regions, where cool superwinds could trigger the formation of molecular clumpy regions.

Classification of Filament Formation Mechanisms in Magnetized Molecular Clouds

Abstract Recent observations of molecular clouds show that dense filaments are the sites of present-day star formation. Thus, it is necessary to understand the filament formation process because these filaments provide the initial condition for star formation. Theoretical research suggests that shock waves in molecular clouds trigger filament formation. Since several different mechanisms have been proposed for filament formation, the formation mechanism of the observed star-forming filaments requires clarification. In the present study, we perform a series of isothermal magnetohydrodynamics simulations of filament formation. We focus on the influences of shock velocity and turbulence on the formation mechanism and identified three different mechanisms for the filament formation. The results indicate that when the shock is fast, at shock velocity v sh ≃ 7 km s−1, the gas flows driven by the curved shock wave create filaments irrespective of the presence of turbulence and self-gravity. However, at a slow shock velocity v sh ≃ 2.5 km s−1, the compressive flow component involved in the initial turbulence induces filament formation. When both the shock velocities and turbulence are low, the self-gravity in the shock-compressed sheet becomes important for filament formation. Moreover, we analyzed the line-mass distribution of the filaments and showed that strong shock waves can naturally create high-line-mass filaments such as those observed in the massive star-forming regions in a short time. We conclude that the dominant filament formation mode changes with the velocity of the shock wave triggering the filament formation.

C II] 158 μm line emission from Orion A I. A template for extragalactic studies?

Context. The [C II] 158 μm fine-structure line is one of the dominant coolants of the neutral interstellar medium. It is hence one of the brightest far-infrared (FIR) emission lines and can be observed not only in star-forming regions throughout the Galaxy, but also in the diffuse interstellar medium and in distant galaxies. [C II] line emission has been suggested to be a powerful tracer of star formation. Aims. We aim to understand the origin of [C II] emission and its relation to other tracers of interstellar gas and dust. This includes a study of the heating efficiency of interstellar gas as traced by the [C II] line to test models of gas heating. Methods. We made use of a one-square-degree map of velocity-resolved [C II] line emission toward the Orion Nebula complex, including M 43 and NGC 1977. We employed Herschel FIR photometric images to determine dust properties. Moreover, we compared with Hα emission from the ionized gas, Spitzer mid-infrared photometry to trace hot dust and large polycyclic aromatic hydrocarbons (PAHs), and velocity-resolved IRAM 30m CO(2–1) observations of the molecular gas. Results. The [C II] intensity is tightly correlated with PAH emission in the IRAC 8 μm band and FIR emission from warm dust. However, the [C II] intensity depends less than linearly on the 8 μm and FIR intensity, while 8 μm and FIR intensities are approximately linearly correlated. The correlation between [C II] and CO(2–1) does not show a clear trend and is affected by the detailed geometry of the region. We find particularly low [C II]-over-FIR intensity ratios toward large columns of (warm and cold) dust, which suggest the interpretation of the “[C II] deficit” in terms of a “FIR excess”. Conclusions. In terms of the [C II] deficit, we find clear evidence in our data for the importance of [O I] 63 μm emission in the photodissociation regions (PDRs) associated with the Huygens region. A smaller contribution is made by a decreased heating efficiency in regions of high UV irradiation. FIR emission from deeply embedded protostars leads to palpably deficient [C II]/FIR intensity ratios. The [C II] directly associated with the M 42, M 43, and NGC 1977 regions underestimates the star formation rate derived from extragalactic scaling relations. We ascribe this to the importance of [C II] emission from low surface brightness PDR surfaces of molecular clouds which are not included in our survey. Future studies of more active regions of massive star formation will be instrumental in validating the general applicability of these conclusions.

CII] emission properties of the massive star-forming region RCW 36 in a filamentary molecular cloud

Aims. We investigate the properties of [C II] 158 μm emission of RCW 36 in a dense filamentary cloud. Methods. [C II] observations of RCW 36, covering an area of ~30′ × 30′, were carried out with a Fabry-Pérot spectrometer on board a 100-cm balloon-borne far-infrared (IR) telescope with an angular resolution of 90′′. Using AKARI and Herschel images, we compared the spatial distribution of the [C II] intensity with the emission from the large grains and polycyclic aromatic hydrocarbon (PAH). Results. The [C II] emission is in good spatial agreement with shell-like structures of a bipolar lobe observed in IR images, which extend along the direction perpendicular to the direction of cold dense filament. We found that the [C II]–160 μm relation for RCW 36 shows a higher brightness ratio of [C II]/160 μm than that for RCW 38, while the [C II]–9 μm relation for RCW 36 is in good agreement with that for RCW 38. Conclusions. Via a spectral decomposition analysis on a pixel-by-pixel basis using IR images, the [C II] emission is spatially well correlated with PAH and cold dust emissions. This means that the observed [C II] emission predominantly comes from photo-dissociation regions. Moreover, the L[C II]∕LFIR ratio shows large variation (10−2–10−3), as compared with the L[C II]/LPAH ratio. In view of the observed tight correlation between L[C II]∕LFIR and the optical depth at λ = 160 μm, the large variation in L[C II]∕LFIR can be simply explained by the geometrical effect, that is, LFIR has contributions from the entire dust-cloud column along the line of sight, while L[C II] has contributions from far-UV illuminated cloud surfaces. Based on the picture of the geometry effect, the enhanced brightness ratio of [C II]/160 μm is attributed to the difference in gas structures where massive stars are formed: filamentary (RCW 36) and clumpy (RCW 38) molecular clouds; thus suggesting that RCW 36 is dominated by far-UV illuminated cloud surfaces, as compared with RCW 38.

A high pitch angle structure in the Sagittarius Arm

Context. In spiral galaxies, star formation tends to trace features of the spiral pattern, including arms, spurs, feathers, and branches. However, in our own Milky Way, it has been challenging to connect individual star-forming regions to their larger Galactic environment owing to our perspective from within the disk. One feature in nearly all modern models of the Milky Way is the Sagittarius Arm, located inward of the Sun with a pitch angle of ∼12°. Aims. We map the 3D locations and velocities of star-forming regions in a segment of the Sagittarius Arm using young stellar objects (YSOs) from the Spitzer/IRAC Candidate YSO (SPICY) catalog to compare their distribution to models of the arm. Methods. Distances and velocities for these objects are derived from Gaia EDR3 astrometry and molecular line surveys. We infer parallaxes and proper motions for spatially clustered groups of YSOs and estimate their radial velocities from the velocities of spatially associated molecular clouds. Results. We identify 25 star-forming regions in the Galactic longitude range ℓ ∼ 4.​ ° 0–18.​ ° 5 arranged in a narrow, ∼1 kpc long linear structure with a high pitch angle of ψ = 56° and a high aspect ratio of ∼7:1. This structure includes massive star-forming regions such as M8, M16, M17, and M20. The motions in the structure are remarkably coherent, with velocities in the direction of Galactic rotation of |Vφ|≈240 ± 3 km s−1 (slightly higher than average) and slight drifts inward (VR ≈ −4.3 km s−1) and in the negative Z direction (VZ ≈ −2.9 km s−1). The rotational shear experienced by the structure is ΔΩ = 4.6 km s−1 kpc−1. Conclusions. The observed 56° pitch angle is remarkably high for a segment of the Sagittarius Arm. We discuss possible interpretations of this feature as a substructure within the lower pitch angle Sagittarius Arm, as a spur, or as an isolated structure.

ALMA View of the ρ Ophiuchi A PDR with a 360 au Beam: The [C i] Emission Originates from the Plane-parallel PDR and Extended Gas

Abstract We present the results of data analysis of the [C i] (3 P 1–3 P 0) emission from the ρ Ophiuchi A photon-dominated region (PDR) obtained in the ALMA ACA standalone mode with a spatial resolution of 2.″6 (360 au). The [C i] emission shows filamentary structures with a width of ∼1000 au, which are adjacent to the shell structure seen in the 4.5 μm map. We found that the 4.5 μm emission, C0, and CO are distributed in this order from the excitation star (S1) in a complementary pattern. These results indicate that [C i] is emitted from a thin layer in the PDR generated by the excitation star, as predicted in the plane-parallel PDR model. In addition, extended [C i] emission was also detected, which shows nearly uniform integrated intensity over the entire field of view (1.′6 × 1.′6). The line profile of the extended component is different from that of the above shell component. The column density ratio of C0 to CO in the extended component was ∼2, which is significantly higher than those of Galactic massive star-forming regions (0.1–0.2). These results suggest that [C i] is emitted also from the extended gas with a density of cm−3, which is not greatly affected by the excitation star.

SOFIA FEEDBACK Survey: Exploring the Dynamics of the Stellar Wind–Driven Shell of RCW 49

Abstract We unveil the stellar wind–driven shell of the luminous massive star-forming region of RCW 49 using SOFIA FEEDBACK observations of the [C ii] 158 μm line. The complementary data set of the 12CO and 13CO J = 3 → 2 transitions is observed by the APEX telescope and probes the dense gas toward RCW 49. Using the spatial and spectral resolution provided by the SOFIA and APEX telescopes, we disentangle the shell from a complex set of individual components of gas centered around RCW 49. We find that the shell of radius ∼6 pc is expanding at a velocity of 13 km s−1 toward the observer. Comparing our observed data with the ancillary data at X-ray, infrared, submillimeter, and radio wavelengths, we investigate the morphology of the region. The shell has a well-defined eastern arc, while the western side is blown open and venting plasma further into the west. Though the stellar cluster, which is ∼2 Myr old, gave rise to the shell, it only gained momentum relatively recently, as we calculate the shell’s expansion lifetime of ∼0.27 Myr, making the Wolf–Rayet star WR 20a a likely candidate responsible for the shell’s reacceleration.

Probing the hadronic nature of the gamma-ray emission associated with Westerlund 2

ABSTRACT Star-forming regions have been proposed as potential Galactic cosmic ray accelerators for decades. Cosmic ray acceleration can be probed through observations of gamma-rays produced in inelastic proton–proton collisions at GeV and TeV energies. In this paper, we analyse more than 11 yr of Fermi–LAT data from the direction of Westerlund 2, one of the most massive and best-studied star-forming regions in our Galaxy. In particular, we investigate the characteristics of the bright pulsar PSR J1023–5746 that dominates the gamma-ray emission below a few GeV at the position of Westerlund 2 and the underlying extended source FGES J1023.3–5747. The analysis results in a clear identification of FGES J1023.3–5747 as the GeV counterpart of the TeV source HESS J1023-575, through its morphological and spectral properties. This identification provides new clues about the origin of the HESS J1023-575 gamma-ray emission, favouring a hadronic origin of the emission, powered by Westerlund 2, rather than a leptonic origin related to either the pulsar wind nebula associated with PSR J1023–5746 or the cluster itself. This result indirectly supports the hypothesis that star-forming regions can contribute to the cosmic ray sea observed in our Galaxy.

ATOMS: ALMA three-millimeter observations of massive star-forming regions – III. Catalogues of candidate hot molecular cores and hyper/ultra compact H ii regions

ABSTRACT We have identified 453 compact dense cores in 3 mm continuum emission maps in the ALMA Three-millimetre Observations of Massive Star-forming regions survey, and compiled three catalogues of high-mass star-forming cores. One catalogue, referred to as hyper/ultra compact (H/UC)-H ii catalogue, includes 89 cores that enshroud H/UC H ii regions as characterized by associated compact H40α emission. A second catalogue, referred to as pure s-cHMC, includes 32 candidate hot molecular cores (HMCs) showing rich spectra (N ≥ 20 lines) of complex organic molecules (COMs) and not associated with H/UC-H ii regions. The third catalogue, referred to as pure w-cHMC, includes 58 candidate HMCs with relatively low levels of COM richness and not associated with H/UC-H ii regions. These three catalogues of dense cores provide an important foundation for future studies of the early stages of high-mass star formation across the Milky Way. We also find that nearly half of H/UC-H ii cores are candidate HMCs. From the number counts of COM-containing and H/UC-H ii cores, we suggest that the duration of high-mass protostellar cores showing chemically rich features is at least comparable to the lifetime of H/UC-H ii regions. For cores in the H/UC-H ii catalogue, the width of the H40α line increases as the core size decreases, suggesting that the non-thermal dynamical and/or pressure line-broadening mechanisms dominate on the smaller scales of the H/UC-H ii cores.

Magnetic Fields in Massive Star-forming Regions (MagMaR). I. Linear Polarized Imaging of the Ultracompact H ii Region G5.89–0.39

Abstract We report 1.2 mm polarized continuum emission observations carried out with the Atacama Large Millimeter/submillimeter Array toward the high-mass star formation region G5.89–0.39. The observations show a prominent 0.2 pc north–south filamentary structure. The ultracompact H ii region in G5.89–0.39 breaks the filament into two pieces. Its millimeter emission shows a dusty belt with a mass of 55–115 M ⊙ and 4500 au in radius, surrounding an inner part comprising mostly ionized gas, with dust emission only accounting for about 30% of the total millimeter emission. We also found a lattice of convex arches that may be produced by dragged dust and gas from the explosive dispersal event involving the O5 Feldt’s star. The north–south filament has a mass between 300 and 600 M ⊙ and harbors a cluster of about 20 mm envelopes with a median size and mass of 1700 au and 1.5 M ⊙, respectively, some of which are already forming protostars. We interpret the polarized emission in the filament as mainly coming from magnetically aligned dust grains. The polarization fraction is ∼4.4% in the filaments and 2.1% at the shell. The magnetic fields are along the North Filament and perpendicular to the South Filament. In the Central Shell, the magnetic fields are roughly radial in a ring surrounding the dusty belt between 4500 and 7500 au, similar to the pattern recently found in the surroundings of Orion BN/KL. This may be an independent observational signpost of explosive dispersal outflows and should be further investigated in other regions.

A Low-mass Cold and Quiescent Core Population in a Massive Star Protocluster

Abstract Pre-stellar cores represent the initial conditions of star formation. Although these initial conditions in nearby low-mass star-forming regions have been investigated in detail, such initial conditions remain vastly unexplored for massive star-forming regions. We report the detection of a cluster of low-mass starless and pre-stellar core candidates in a massive star protocluster-forming cloud, NGC 6334S. With the Atacama Large Millimeter/submillimeter Array (ALMA) observations at a ∼0.02 pc spatial resolution, we identified 17 low-mass starless core candidates that do not show any evidence of protostellar activity. These candidates present small velocity dispersions, high fractional abundances of NH2D, high NH3 deuterium fractionations, and are completely dark in the infrared wavelengths from 3.6 up to 70 μm. Turbulence is significantly dissipated and the gas kinematics are dominated by thermal motions toward these candidates. Nine out of the 17 cores are gravitationally bound, and therefore are identified as pre-stellar core candidates. The embedded cores of NGC 6334S show a wide diversity in masses and evolutionary stages.