Search

The contribution of the magnetic field to the formation of high-mass stars is poorly understood. We report the high angular resolution (∼0.″3, 870 au) map of the magnetic field projected on the plane of the sky ( B _POS ) toward the high-mass star-forming region G333.46−0.16 (G333), obtained with the Atacama Large Millimeter/submillimeter Array at 1.2 mm as part of the Magnetic fields in Massive star-forming Regions survey. The B _POS morphology found in this region is consistent with a canonical “hourglass” with an embedded flattened envelope in a perpendicular direction, which suggests a dynamically important field. This region is fragmented into two protostars that appear to be gravitationally bound in a stable binary system with a separation of ∼1740 au. Interestingly, by analyzing H ^13 CO ^+ ( J = 3–2) line emission, we find no velocity gradient over the extent of the continuum, which is consistent with a strong field. We model the B _POS , obtaining a marginally supercritical mass-to-flux ratio of 1.43, suggesting an initially strongly magnetized environment. Based on the Davis–Chandrasekhar–Fermi method, the magnetic field strength toward G333 is estimated to be 5.7 mG. The absence of strong rotation and outflows toward the central region of G333 suggests strong magnetic braking, consistent with a highly magnetized environment. Our study shows that despite being a strong regulator, the magnetic energy fails to prevent the process of fragmentation, as revealed by the formation of the two protostars in the central region.

Using new continuum and molecular line data from the Atacama Large Millimeter/submillimeter Array Three-millimeter Observations of Massive Star-forming Regions (ATOMS) survey and archival Very Large Array, 4.86 GHz data, we present direct observational evidence of hierarchical triggering relating three epochs of massive star formation in a ringlike H ii region, G24.47+0.49. We find from radio flux analysis that it is excited by a massive star(s) of spectral type O8.5V–O8V from the first epoch of star formation. The swept-up ionized ring structure shows evidence of secondary collapse, and within this ring, a burst of massive star formation is observed in different evolutionary phases, which constitutes the second epoch. ATOMS spectral line (e.g., HCO ^+ (1–0)) observations reveal an outer concentric molecular gas ring expanding at a velocity of ∼9 km s ^−1 , constituting the direct and unambiguous detection of an expanding molecular ring. It harbors twelve dense molecular cores with surface mass density greater than 0.05 g cm ^−2 , a threshold typical of massive star formation. Half of them are found to be subvirial and thus in gravitational collapse making them the third epoch of potential massive star-forming sites.

Only a handful of massive starless core candidates have been discovered so far, but none of them have been fully confirmed. Within the MM1 clump in the filamentary infrared dark cloud G34.43+0.24 that was covered by the Atacama Large Millimeter/submillimeter Array (ALMA) ATOMS survey at Band 3 (∼2″, 6000 au) and the ALMA-QUARKS survey at Band 6 (∼0.″3, 900 au), two prestellar core candidates MM1-C and E1 with masses of 71 and 20 M _⊙ and radii of 2100–4400 au were discovered. The two cores show no obvious sign of star formation activities. In particular, MM1-C is a very promising massive prestellar core candidate with a total gas mass of 71 M _⊙ . Within MM1-C, we detected two extremely dense substructures, C1 and C2, as characterized by their high densities of ${n}_{{{\rm{H}}}_{2}}\sim {10}^{8\mbox{--}9}\,{\mathrm{cm}}^{-3}$ . Moreover, evidence of further fragmentation in C2 was also revealed. We have detected the primordial fragmentation in the earliest stage of massive star formation, and we speculate that MM1-C would be the birthplace of a massive multiple system. However, we cannot fully rule out the possibility that the massive prestellar core MM1-C will just form a cluster of low-mass stars if it undergoes further fragmentation.

The formation of the massive stars, and in particular, the role that the magnetic fields play in their early evolutionary phase is still far from being completely understood. Here, we present the Atacama Large Millimeter/submillimeter Array 1.2 mm full polarized continuum and H ^13 CO ^+ (3−2), CS(5−4), and HN ^13 C(3−2) line observations with a high angular resolution (∼0.″4 or 1100 au). In the 1.2 mm continuum emission, we reveal a dusty envelope surrounding the massive protostars, IRAS16547-E and IRAS16547-W, with dimensions of ∼10,000 au. This envelope has a biconical structure likely carved by the powerful thermal radio jet present in region. The magnetic field vectors follow very well the biconical envelope. The polarization fraction is ∼2.0% in this region. Some of these vectors seem to converge to IRAS 16547-E and IRAS 16547-W, the most massive protostars. Moreover, the velocity fields revealed from the spectral lines H ^13 CO ^+ (3−2) and HN ^13 C(3−2) show velocity gradients with a good correspondence with the magnetic fields, which maybe are tracing the cavities of molecular outflows or maybe infalling in some parts. We derived a magnetic field strength in some filamentary regions that goes from 2 to 6.1 mG. We also find that the CS(5−4) molecular line emission reveals multiple outflow cavities or bow shocks with different orientations, some of which seem to follow the NW–SE radio thermal jet.

Fragmentation contributes to the formation and evolution of stars. Observationally, high-mass stars are known to form multiple-star systems, preferentially in cluster environments. Theoretically, Jeans instability has been suggested to determine characteristic fragmentation scales, and thermal or turbulent motion in the parental gas clump mainly contributes to the instability. To search for such a characteristic fragmentation scale, we have analyzed Atacama Large Millimeter/submillimeter Array (ALMA) 1.33 mm continuum observations toward 30 high-mass star-forming clumps taken by the Digging into the Interior of Hot Cores with ALMA survey. We have identified 573 cores using the dendrogram algorithm and measured the separation of cores by using the Minimum Spanning Tree technique. The core separation corrected by projection effects has a distribution peaked around 5800 au. In order to remove biases produced by different distances and sensitivities, we further smooth the images to a common physical scale and perform completeness tests. Our careful analysis finds a characteristic fragmentation scale of ∼7000 au, comparable to the thermal Jeans length of the clumps. We conclude that thermal Jeans fragmentation plays a dominant role in determining the clump fragmentation in high-mass star-forming regions, without the need to invoke turbulent Jeans fragmentation.

We report on Atacama Large Millimeter/submillimeter Array 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 northwest to southeast with intricate disturbances across major cores, hinting at spiral structures. Energy analysis of CS ( J = 5 → 4) emission yields an outflow energy of approximately 3.5 × 10 ^45 erg, 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 ± 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 may be 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.

Filamentary structures are ubiquitously found in high-mass star-forming clouds. To investigate the relationship between filaments and star formation, we carry out the INFANT (Investigations of Massive Filaments and Sar Formation) survey, a multiscale, multiwavelength survey of massive filamentary clouds with Atacama Large Millimeter/submillimeter Array (ALMA) band 3/band 6 and Very Large Array K band. In this first paper, we present the ALMA band 6 continuum observations toward a sample of eight high-mass star-forming filaments. We covered each target with an approximately rectangular mosaic field of view with two 12 m array configurations, achieving an angular resolution of ∼0.″6 (2700 au at 4.5 kpc) and a continuum rms of ∼0.1 mJy beam ^−1 (∼0.06 M _⊙ in gas mass assuming 15 K). We identify cores using the getsf and astrodendro and find the former is more robust in terms of both identification and measuring flux densities. We identify in total 183 dense cores (15–36 cores in each cloud) and classify their star formation states via outflow and warm gas tracers. The protostellar cores are statistically more massive than the prestellar cores, possibly indicating further accretion onto cores after the formation of protostars. For the high-mass end ( ${M}_{\mathrm{core}}$ > 1.5 M _⊙ ) of the core mass function (CMF) we derive a power-law index of −1.15 ± 0.12 for the whole sample and −1.70 ± 0.25 for the prestellar population. We also find a steepening trend in CMF with cloud evolution (−0.89 ± 0.15 for the young group versus. −1.44 ± 0.25 for the evolved group) and discuss its implications for cluster formation.

Fragmentation during the early stages of high-mass star formation is crucial for understanding the formation of high-mass clusters. We investigated fragmentation within 39 high-mass star-forming clumps as part of the Atacama Large Millimeter/submillimeter Array Survey of 70 μ m Dark High-mass Clumps in Early Stages (ASHES) survey. 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 subclump structures revealed by the 7 m array continuum emission. Comparison of the Jeans parameters using clump and subclump densities with the separation and masses of gravitationally bound cores suggests that they can be explained by clump fragmentation, implying the simultaneous formation of subclumps 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 subclustered. Using the ${ \mathcal Q }$ -parameter, which can help 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.

We have comprehensively studied the multiscale physical properties of the massive infrared dark cloud G28.34 (the Dragon cloud) with dust polarization and molecular line data from Planck, FCRAO-14 m, James Clerk Maxwell Telescope, and Atacama Large Millimeter/submillimeter Array. We find that the averaged magnetic fields of clumps tend to be either parallel with or perpendicular to the cloud-scale magnetic fields, while the cores in clump MM4 tend to have magnetic fields aligned with the clump fields. Implementing the relative orientation analysis (for magnetic fields, column density gradients, and local gravity), velocity gradient technique, and modified Davis–Chandrasekhar–Fermi analysis, we find that G28.34 is located in a trans-to-sub-Alfvénic environment; the magnetic field is effectively resisting gravitational collapse in large-scale diffuse gas, but is distorted by gravity within the cloud and affected by star formation activities in high-density regions, and the normalized mass-to-flux ratio tends to increase with increasing density and decreasing radius. Considering the thermal, magnetic, and turbulent supports, we find that the environmental gas of G28.34 is in a supervirial (supported) state, the infrared dark clumps may be in a near-equilibrium state, and core MM4-core4 is in a subvirial (gravity-dominant) state. In summary, we suggest that magnetic fields dominate gravity and turbulence in the cloud environment at large scales, resulting in relatively slow cloud formation and evolution processes. Within the cloud, gravity could overwhelm both magnetic fields and turbulence, allowing local dynamical star formation to happen.

Massive infrared dark clouds (IRDCs) are considered to host the earliest stages of high-mass star formation. In particular, 70 μ m dark IRDCs are the colder and more quiescent clouds. At a scale of about 5000 au using formaldehyde (H _2 CO) emission, we investigate the kinetic temperature of dense cores in 12 IRDCs obtained from the pilot Atacama Large Millimeter/submillimeter Array Survey of 70 μ m dark High-mass clumps in Early Stages (ASHES). Compared to the 1.3 mm dust continuum and other molecular lines, such as C ^18 O and deuterated species, we find that H _2 CO is mainly sensitive to low-velocity outflow components rather than to quiescent gas expected in the early phases of star formation. The kinetic temperatures of these components range from 26 to 300 K. The Mach number reaches about 15 with an average value of about 4, suggesting that the velocity distribution of gas traced by H _2 CO is significantly influenced by a supersonic nonthermal component. In addition, we detect warm line emission from HC _3 N and OCS in 14 protostellar cores, which requires high excitation temperatures ( E _u / k ∼ 100 K). These results show that some of the embedded cores in the ASHES fields are in an advanced evolutionary stage, previously unexpected for 70 μ m dark IRDCs.

We conducted a Ka -band (26.1–35 GHz) line survey toward Orion KL using the TianMa 65 m Radio Telescope (TMRT). It is the first blind line survey in the Ka band and achieves a sensitivity at the mK level (1–3 mK at a spectral resolution of ∼1 km s ^−1 ). In total, 592 Gaussian features are extracted. Among them, 257 radio recombination lines (RRLs) are identified. The maximum Δ n of RRLs of H, He, and C are 20, 15, and 5, respectively. Through stacking, we have detected the β lines of ion RRLs (RRLs of C ^+ with the possible contribution of other ions like O ^+ ) for the first time, and a tentative signal of the γ lines of ion RRLs can also be seen on the stacked spectrum. Besides this, 318 other line features were assigned to 37 molecular species, and 10 of these species were not detected in the Q -band survey of TMRT. The vibrationally excited states of nine species were also detected. The emission of most species can be modeled under LTE. A number of transitions of E-CH3OH ( J _2 − J _1 ) display maser effects, which are confirmed by our modeling, and besides the bumping peak at J ∼ 6, there is another peak at J ∼ 13. Methylcyanoacetylene (CH _3 C _3 N) is detected in Orion KL for the first time. This work emphasizes that the Ka band, which was long ignored for spectral line surveys, is very useful for surveying RRLs and molecular lines simultaneously.

Protostellar outflows and jets are almost ubiquitous characteristics during the mass accretion phase and encode the history of stellar accretion, complex organic molecule (COM) formation, and planet formation. Episodic jets are likely connected to episodic accretion through the disk. Despite the importance, studies on episodic accretion and ejection links have not been done yet in a systematic fashion using high-sensitivity and high-resolution observations. To explore episodic accretion mechanisms and the chronologies of episodic events, we investigated 39 fields containing protostars with Atacama Large Millimeter/submillimeter Array observations of CO, SiO, and 1.3 mm continuum emission. We detected SiO emission in 19 fields, where 17 sources are driving molecular jets. Jet velocities, mass-loss rates, mass accretion rates, and periods of accretion events appear to have some dependence on the driving forces of the jet (e.g., bolometric luminosity, envelope mass). Next, velocities and mass-loss rates appear to be somewhat correlated with the surrounding envelope mass, suggesting that the presence of high mass around protostars increases the ejection–accretion activity. We determine mean periods of ejection events of 20–175 yr for our sample, which could be associated with perturbation zones of ∼2−25 au extent around the protostars. In addition, mean ejection periods show an apparent anticorrelation with the envelope mass, where high accretion rates may trigger more frequent ejection events. The observed periods of outburst/ejection are much shorter than the freezeout timescale of the simplest COMs like CH _3 OH, suggesting that episodic events could affect the ice–gas balance inside and around the snowline.

The internal structure of the prestellar core G208.68-19.02-N2 (G208-N2) in the Orion Molecular Cloud 3 (OMC-3) region has been studied with the Atacama Large Millimeter/submillimeter Array. The dust continuum emission revealed a filamentary structure with a length of ∼5000 au and an average H _2 volume density of ∼6 × 10 ^7 cm ^−3 . At the tip of this filamentary structure, there is a compact object, which we call a nucleus , with a radius of ∼150–200 au and a mass of ∼0.1 M _⊙ . The nucleus has a central density of ∼2 × 10 ^9 cm ^−3 with a radial density profile of r ^−1.87±0.11 . The density scaling of the nucleus is ∼3.7 times higher than that of the singular isothermal sphere (SIS). This as well as the very low virial parameter of 0.39 suggests that the gravity is dominant over the pressure everywhere in the nucleus. However, there is no sign of CO outflow localized to this nucleus. The filamentary structure is traced by the N _2 D ^+ 3–2 emission, but not by the C ^18 O 2–1 emission, implying the significant CO depletion due to high density and cold temperature. Toward the nucleus, the N _2 D ^+ also shows the signature of depletion. This could imply either the depletion of the parent molecule, N _2 , or the presence of the embedded very-low luminosity central source that could sublimate the CO in the very small area. The nucleus in G208-N2 is considered to be a prestellar core on the verge of first hydrostatic core (FHSC) formation or a candidate for the FHSC.

We observed the high-mass star-forming core G336.01–0.82 at 1.3 mm and 0.″05 (∼150 au) angular resolution with the Atacama Large Millimeter/submillimeter Array (ALMA) as part of the Digging into the Interior of Hot Cores with ALMA survey. These high-resolution observations reveal two spiral streamers feeding a circumstellar disk at opposite sides in great detail. Molecular line emission from CH _3 OH shows velocity gradients along the streamers consistent with infall. Similarly, a flattened envelope model with rotation and infall implies a mass larger than 10 M _☉ for the central source and a centrifugal barrier of 300 au. The location of the centrifugal barrier is consistent with local peaks in the continuum emission. We argue that gas brought by the spiral streamers is accumulating at the centrifugal barrier, which can result in future accretion burst events. A total high infall rate of ∼4 × 10 ^−4 M _☉ yr ^−1 is derived by matching models to the observed velocity gradient along the streamers. Their contribution accounts for 20%–50% the global infall rate of the core, indicating streamers play an important role in the formation of high-mass stars.

One of the most poorly understood aspects of low-mass star formation is how multiple-star systems are formed. Here we present the results of Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 observations toward a forming quadruple protostellar system, G206.93-16.61E2, in the Orion B molecular cloud. ALMA 1.3 mm continuum emission reveals four compact objects, of which two are Class I young stellar objects and the other two are likely in prestellar phase. The 1.3 mm continuum emission also shows three asymmetric ribbon-like structures that are connected to the four objects, with lengths ranging from ∼500 to ∼2200 au. By comparing our data with magnetohydrodynamic simulations, we suggest that these ribbons trace accretion flows and also function as gas bridges connecting the member protostars. Additionally, ALMA CO J = 2−1 line emission reveals a complicated molecular outflow associated with G206.93-16.61E2, with arc-like structures suggestive of an outflow cavity viewed pole-on.

The Atacama Large Millimeter/submillimeter Array (ALMA) Survey of Star Formation and Evolution in Massive Protoclusters with Blue Profiles (ASSEMBLE) aims to investigate the process of mass assembly and its connection to high-mass star formation theories in protoclusters in a dynamic view. We observed 11 massive ( M _clump ≳ 10 ^3 M _⊙ ), luminous ( L _bol ≳ 10 ^4 L _⊙ ), and blue-profile (infall signature) clumps by ALMA with resolution of ∼2200–5500 au (median value of 3500 au) at 350 GHz (870 μ m). We identified 248 dense cores, including 106 cores showing protostellar signatures and 142 prestellar core candidates. Compared to early stage infrared dark clouds (IRDCs) by ASHES, the core mass and surface density within the ASSEMBLE clumps exhibited a significant increment, suggesting concurrent core accretion during the evolution of the clumps. The maximum mass of prestellar cores was found to be 2 times larger than that in IRDCs, indicating that evolved protoclusters have the potential to harbor massive prestellar cores. The mass relation between clumps and their most massive core (MMCs) is observed in ASSEMBLE but not in IRDCs, which is suggested to be regulated by multiscale mass accretion. The mass correlation between the core clusters and their MMCs has a steeper slope compared to that observed in stellar clusters, which can be due to fragmentation of the MMC and stellar multiplicity. We observe a decrease in core separation and an increase in central concentration as protoclusters evolve. We confirm primordial mass segregation in the ASSEMBLE protoclusters, possibly resulting from gravitational concentration and/or gas accretion.

Hot corinos are of great interest due to their richness in interstellar complex organic molecules (COMs) and the consequent potential prebiotic connection to solar-like planetary systems. Recent surveys have reported an increasing number of detected hot corinos in Class 0/I protostars; however, the relationships between their physical properties and the hot-corino signatures remain elusive. In this study, our objective is to establish a general picture of the detectability of hot corinos by identifying the origins of the hot-corino signatures in the sample of young stellar objects (YSOs) obtained from the Atacama Large Millimeter/submillimeter Array Survey of Orion Planck Galactic Cold Clumps project. We apply spectral energy distribution modeling to our sample and identify the physical parameters of the modeled YSOs directly, linking the detection of hot-corino signatures to the envelope properties of the YSOs. Imaging simulations of the methanol emission further support this scenario. We therefore posit that the observed COM emission originates from the warm inner envelopes of the sample YSOs, based on both the warm region size and the envelope density profile. The former is governed by the source luminosity and is additionally affected by the disk and cavity properties, while the latter is related to the evolutionary stages. This scenario provides a framework for detecting hot-corino signatures toward luminous Class 0 YSOs, with fewer detections being observed toward similarly luminous Class I sources.

We studied the unique kinematic properties in massive filament G352.63-1.07 at 10 ^3 au spatial scale with the dense molecular tracers observed with the Atacama Large Millimeter/submillimeter Array. We find the central massive core M1 (12 M _⊙ ) being separated from the surrounding filament with a velocity difference of $v-{\overline{v}}_{\mathrm{sys}}=-2\,\mathrm{km}\,{{\rm{s}}}^{-1}$ and a transverse separation within 3″. Meanwhile, as shown in multiple dense-gas tracers, M1 has a spatial extension closely aligned with the main filament and is connected to the filament toward both its ends. M1 thus represents a very beginning state for a massive, young star-forming core escaping from the parental filament, within a timescale of ∼4000 yr. Based on its kinetic energy (3.5 × 10 ^44 erg), the core escape is unlikely solely due to the original filament motion or magnetic field but requires more energetic events such as a rapid intense anisotropic collapse. The released energy also seems to noticeably increase the environmental turbulence. This may help the filament to become stabilized again.

There is growing evidence that high-mass star formation and hub-filament systems (HFS) are intricately linked. The gas kinematics along the filaments and the forming high-mass star(s) in the central hub are in excellent agreement with the new generation of global hierarchical high-mass star formation models. In this paper, we present an observational investigation of a typical HFS cloud, G310.142+0.758 (G310 hereafter), which reveals unambiguous evidence of mass inflow from the cloud scale via the filaments onto the forming protostar(s) at the hub conforming with the model predictions. Continuum and molecular line data from the ATOMS and MALT90 surveys that cover different spatial scales are used. Three filaments (with a total mass of 5.7 ± 1.1 × 10 ^3 M _⊙ ) are identified converging toward the central hub region where several signposts of high-mass star formation have been observed. The hub region contains a massive clump (1280 ± 260 M _⊙ ) harboring a central massive core. Additionally, five outflow lobes are associated with the central massive core implying a forming cluster. The observed large-scale, smooth, and coherent velocity gradients from the cloud down to the core scale, and the signatures of infall motion seen in the central massive clump and core, clearly unveil a nearly continuous, multi-scale mass accretion/transfer process at a similar mass infall rate of ∼10 ^−3 M _⊙ yr ^−1 over all scales, feeding the central forming high-mass protostar(s) in the G310 HFS cloud.

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 12 and 7 m arrays from the Atacama Large Millimeter/submillimeter Array. The targets are 70 μ m dark massive (220–4900 M _⊙ ), dense (>10 ^4 cm ^−3 ), and cold (∼10–20 K) 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 _⊙ using 1.3 mm dust continuum emission. About 55% of the cores are low-mass (

We have analyzed the NH _2 CHO, HNCO, H _2 CO, and CH _3 CN ( ^13 CH _3 CN) molecular lines at an angular resolution of ∼0.″3 obtained by the Atacama Large Millimeter/submillimeter Array Band 6 toward 30 high-mass star-forming regions. The NH _2 CHO emission has been detected in 23 regions, while the other species have been detected toward 29 regions. A total of 44 hot molecular cores (HMCs) have been identified using the moment 0 maps of the CH _3 CN line. The fractional abundances of the four species have been derived at each HMC. In order to investigate pure chemical relationships, we have conducted a partial correlation test to exclude the effect of temperature. Strong positive correlations between NH _2 CHO and HNCO ( ρ = 0.89) and between NH _2 CHO and H _2 CO (0.84) have been found. These strong correlations indicate their direct chemical links; dual-cyclic hydrogen addition and abstraction reactions between HNCO and NH _2 CHO and gas-phase formation of NH _2 CHO from H _2 CO. Chemical models including these reactions can reproduce the observed abundances in our target sources.

We present dynamical properties of 294 cores embedded in twelve IRDCs observed as part of the ASHES Survey. Protostellar cores have higher gas masses, surface densities, column densities, and volume densities than prestellar cores, indicating core mass growth from the prestellar to the protostellar phase. We find that ∼80% of cores with virial parameter ( α ) measurements are gravitationally bound ( α < 2). We also find an anticorrelation between the mass and the virial parameter of cores, with massive cores having on average lower virial parameters. Protostellar cores are more gravitationally bound than prestellar cores, with an average virial parameter of 1.2 and 1.5, respectively. The observed nonthermal velocity dispersion (from N _2 D ^+ or DCO ^+ ) is consistent with simulations in which turbulence is continuously injected, whereas the core-to-core velocity dispersion is neither in agreement with driven nor decaying turbulence simulations. We find a not significant increment in the line velocity dispersion from prestellar to protostellar cores, suggesting that the dense gas within the core traced by these deuterated molecules is not yet severely affected by turbulence injected from outflow activity at the early evolutionary stages traced in ASHES. The most massive cores are strongly self-gravitating and have greater surface density, Mach number, and velocity dispersion than cores with lower masses. Dense cores do not have significant velocity shifts relative to their low-density envelopes, suggesting that dense cores are comoving with their envelopes. We conclude that the observed core properties are more in line with the predictions of clump-fed scenarios rather than with those of core-fed scenarios.

We use molecular line data from the Atacama Large Millimeter/submillimeter Array, Submillimeter Array, James Clerk Maxwell Telescope, and NANTEN2 to study the multiscale (∼15–0.005 pc) velocity statistics in the massive star formation region NGC 6334. We find that the nonthermal motions revealed by the velocity dispersion function (VDF) stay supersonic over scales of several orders of magnitude. The multiscale nonthermal motions revealed by different instruments do not follow the same continuous power law, which is because the massive star formation activities near central young stellar objects have increased the nonthermal motions in small-scale and high-density regions. The magnitudes of VDFs vary in different gas materials at the same scale, where the infrared dark clump N6334S in an early evolutionary stage shows a lower level of nonthermal motions than other more evolved clumps due to its more quiescent star formation activity. We find possible signs of small-scale-driven (e.g., by gravitational accretion or outflows) supersonic turbulence in clump N6334IV with a three-point VDF analysis. Our results clearly show that the scaling relation of velocity fields in NGC 6334 deviates from a continuous and universal turbulence cascade due to massive star formation activities.

We present ALMA dust polarization and molecular line observations toward four clumps (I(N), I, IV, and V) in the massive star-forming region NGC 6334. In conjunction with large-scale dust polarization and molecular line data from JCMT, Planck, and NANTEN2, we make a synergistic analysis of relative orientations between magnetic fields ( θ _B ), column density gradients ( θ _NG ), local gravity ( θ _LG ), and velocity gradients ( θ _VG ) to investigate the multi-scale (from ∼30 to 0.003 pc) physical properties in NGC 6334. We find that the relative orientation between θ _B and θ _NG changes from statistically more perpendicular to parallel as column density ( ${N}_{{{\rm{H}}}_{2}}$ ) increases, which is a signature of trans-to-sub-Alfvénic turbulence at complex/cloud scales as revealed by previous numerical studies. Because θ _NG and θ _LG are preferentially aligned within the NGC 6334 cloud, we suggest that the more parallel alignment between θ _B and θ _NG at higher ${N}_{{{\rm{H}}}_{2}}$ is because the magnetic field line is dragged by gravity. At even higher ${N}_{{{\rm{H}}}_{2}}$ , the angle between θ _B and θ _NG or θ _LG transits back to having no preferred orientation, or statistically slightly more perpendicular, suggesting that the magnetic field structure is impacted by star formation activities. A statistically more perpendicular alignment is found between θ _B and θ _VG throughout our studied ${N}_{{{\rm{H}}}_{2}}$ range, which indicates a trans-to-sub-Alfvénic state at small scales as well, and this signifies that magnetic field has an important role in the star formation process in NGC 6334. The normalized mass-to-flux ratio derived from the polarization-intensity gradient (KTH) method increases with ${N}_{{{\rm{H}}}_{2}}$ , but the KTH method may fail at high ${N}_{{{\rm{H}}}_{2}}$ due to the impact of star formation feedback.

Starless cores represent the initial stage of evolution toward (proto)star formation, and a subset of them, known as prestellar cores, with high density (∼ 10 ^6 cm ^−3 or higher) and being centrally concentrated are expected to be embryos of (proto)stars. Determining the density profile of prestellar cores therefore provides an important opportunity to gauge the initial conditions of star formation. In this work, we perform rigorous modeling to estimate the density profiles of three nearly spherical prestellar cores among a sample of five highly dense cores detected by our recent observations. We employed multiscale observational data of the (sub)millimeter dust continuum emission, including those obtained by SCUBA-2 on the James Clerk Maxwell Telescope with a resolution of ∼ 5600 au and by multiple Atacama Large Millimeter/submillimeter Array observations with a resolution as high as ∼ 480 au. We are able to consistently reproduce the observed multiscale dust continuum images of the cores with a simple prescribed density profile, which bears an inner region of flat density and an r ^−2 profile toward the outer region. By utilizing the peak density and the size of the inner flat region as a proxy for the dynamical stage of the cores, we find that the three modeled cores are most likely unstable and prone to collapse. The sizes of the inner flat regions, as compact as ∼ 500 au, signify them as being the highly evolved prestellar cores rarely found to date.