Search

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 Atacama Large Millimeter/Submillimeter Array (ALMA) 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 ($\sim$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 $\sim$10,000 au. This envelope has a bi-conical structure likely carved by the powerful thermal radio jet present in region. The magnetic fields vectors follow very-well the bi-conical envelope. The polarization fraction is $\sim$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, that maybe are tracing the cavities of molecular outflows or maybe in some parts infall. 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., Comment: Accepted by the Astrophysical Journal, 13 pages

The contribution of the magnetic field to the formation of high-mass stars is poorly understood. We report the high-angular resolution ($\sim0.3^{\prime\prime}$, 870 au) map of the magnetic field projected on the plane of the sky (B$_\mathrm{POS}$) towards the high-mass star forming region G333.46$-$0.16 (G333), obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.2 mm as part of the Magnetic Fields in Massive Star-forming Regions (MagMaR) survey. The B$_\mathrm{POS}$ morphology found in this region is consistent with a canonical ``hourglass'' which suggest a dynamically important field. This region is fragmented into two protostars separated by $\sim1740$ au. Interestingly, by analysing H$^{13}$CO$^{+}$ ($J=3-2$) line emission, we find no velocity gradient over the extend of the continuum which is consistent with a strong field. We model the B$_\mathrm{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 towards G333 is estimated to be 5.7 mG. The absence of strong rotation and outflows towards 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.

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 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 (DIHCA) survey. We have identified 573 cores using the dendrogram algorithm and measured the separation of cores by using the Minimum Spanning Tree (MST) 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 $\sim$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 of invoking turbulent Jeans fragmentation., Comment: 30 pages, 18 figures, Accepted in ApJ

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., Comment: Accepted for Publication at the Astrophysical Journal

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 $\mu$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., Comment: Accepted for Publication in ApJ. 19 pages, 7 figures, 4 tables

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., Comment: Published in Nature Astronomy at https://www.nature.com/articles/s41550-023-02181-9. 27 pages, 12 figures, 1 table

Massive infrared dark clouds (IRDCs) are considered to host the earliest stages of high-mass star formation. In particular, 70 $\mu$m dark IRDCs are the colder and more quiescent clouds. At a scale of about 5000 au using formaldehyde (H2CO) emission, we investigate the kinetic temperature of dense cores in 12 IRDCs obtained from the pilot ALMA Survey of 70 $\mu$m dark High-mass clumps in Early Stages (ASHES). Compared to 1.3 mm dust continuum and other molecular lines, such as C18O and deuterated species, we find that H2CO 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 H2CO is significantly influenced by a supersonic non-thermal component. In addition, we detect warm line emission from HC3N and OCS in 14 protostellar cores, which requires high excitation temperatures (Eu/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 $\mu$m dark IRDCs., Comment: Accepted for Publication in ApJ. 39 pages, 22 figures, 6 tables

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 $\mu$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., Comment: Accepted for Publication in ApJ. 54 pages, 40 figures, 5 tables

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 $\sim 0.3''$ obtained by the Atacama Large Millimeter/submillimeter Array (ALMA) 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 ($\rho=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., Comment: Accepted for The Astrophysical Journal. 27 pages, 10 tables, and 13 figures

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., Comment: accepted to ApJ

We observed the high-mass protostellar core G335.579-0.272 ALMA1 at ${\sim}200$ au (0.05") resolution with the Atacama Large Millimeter/submillimeter Array (ALMA) at 226 GHz (with a mass sensitivity of $5\sigma=0.2$ M$_\odot$ at 10 K). We discovered that at least a binary system is forming inside this region, with an additional nearby bow-like structure (${\lesssim}1000$ au) that could add an additional member to the stellar system. These three sources are located at the center of the gravitational potential well of the ALMA1 region and the larger MM1 cluster. The emission from CH$_3$OH (and many other tracers) is extended ($>1000$ au), revealing a common envelope toward the binary system. We use CH$_2$CHCN line emission to estimate an inclination angle of the rotation axis of $26^\circ$ with respect to the line of sight based on geometric assumptions and derive a kinematic mass of the primary source (protostar+disk) of 3.0 M$_\odot$ within a radius of 230 au. Using SiO emission, we find that the primary source drives the large scale outflow revealed by previous observations. Precession of the binary system likely produces a change in orientation between the outflow at small scales observed here and large scales observed in previous works. The bow structure may have originated by entrainment of matter into the envelope due to widening or precession of the outflow, or, alternatively, an accretion streamer dominated by the gravity of the central sources. An additional third source, forming due to instabilities in the streamer, cannot be ruled out as a temperature gradient is needed to produce the observed absorption spectra., Comment: 23 pages, 16 figures, 4 tables, 2 appendices. Accepted for publication in ApJ

Here, we report ALMA detections of polarized emission from dust, CS($J=5 \rightarrow 4$), and C$^{33}$S($J=5 \rightarrow 4$) toward the high-mass star-forming region NGC6334I(N). A clear ``hourglass'' magnetic field morphology was inferred from the polarized dust emission which is also directly seen from the polarized CS emission across velocity, where the polarization appears to be parallel to the field. By considering previous findings, the field retains a pinched shape which can be traced to clump length-scales from the envelope scales traced by ALMA, suggesting that the field is dynamically important across multiple length-scales in this region. The CS total intensity emission is found to be optically thick ($\tau_{\mathrm{CS}} = 32 \pm 12$) while the C$^{33}$S emission appears to be optically thin ($\tau_{\mathrm{C^{33}S}} = 0.1 \pm 0.01$). This suggests that sources of anisotropy other than large velocity gradients, i.e. anisotropies in the radiation field are required to explain the polarized emission from CS seen by ALMA. By using four variants of the Davis-Chandrasekhar-Fermi technique and the angle dispersion function methods (ADF), we obtain an average of estimates for the magnetic field strength onto the plane of the sky of $\left< \mathrm{B}_{\mathrm{pos}} \right> = 16$ mG from the dust and $\left< \mathrm{B}_{\mathrm{pos}} \right> \sim 2$ mG from the CS emission, where each emission traces different molecular hydrogen number densities. This effectively enables a tomographic view of the magnetic field within a single ALMA observation., Comment: Accepted for publication in ApJ

With a mass of $\sim$1000 $M_\odot$ and a surface density of $\sim$0.5 g cm$^{-2}$, G023.477+0.114 also known as IRDC 18310-4 is an infrared dark cloud (IRDC) that has the potential to form high-mass stars and has been recognized as a promising prestellar clump candidate. To characterize the early stages of high-mass star formation, we have observed G023.477+0.114 as part of the ALMA Survey of 70 $\mu$m Dark High-mass Clumps in Early Stages (ASHES). We have conducted $\sim$1."2 resolution observations with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.3 mm in dust continuum and molecular line emission. We identified 11 cores, whose masses range from 1.1 $M_\odot$ to 19.0 $M_\odot$. Ignoring magnetic fields, the virial parameters of the cores are below unity, implying that the cores are gravitationally bound. However, when magnetic fields are included, the prestellar cores are close to virial equilibrium, while the protostellar cores remain sub-virialized. Star formation activity has already started in this clump. Four collimated outflows are detected in CO and SiO. H$_2$CO and CH$_3$OH emission coincide with the high-velocity components seen in the CO and SiO emission. The outflows are randomly oriented for the natal filament and the magnetic field. The position-velocity diagrams suggest that episodic mass ejection has already begun even in this very early phase of protostellar formation. The masses of the identified cores are comparable to the expected maximum stellar mass that this IRDC could form (8-19 $M_\odot$). We explore two possibilities on how IRDC G023.477+0.114 could eventually form high-mass stars in the context of theoretical scenarios., Comment: Accepted for publication in ApJ (September 1, 2021). 33 pages, 20 figures, and 5 tables

A full understanding of high-mass star formation requires the study of one of the most elusive components of the energy balance in the interstellar medium: magnetic fields. We report ALMA 1.2 mm, high-resolution (700 au) dust polarization and molecular line observations of the rotating hot molecular core embedded in the high-mass star-forming region IRAS 18089-1732. The dust continuum emission and magnetic field morphology present spiral-like features resembling a whirlpool. The velocity field traced by the H13CO+ (J=3-2) transition line reveals a complex structure with spiral filaments that are likely infalling and rotating, dragging the field with them. We have modeled the magnetic field and find that the best model corresponds to a weakly magnetized core with a mass-to-magnetic-flux ratio (lambda) of 8.38. The modeled magnetic field is dominated by a poloidal component, but with an important contribution from the toroidal component that has a magnitude of 30% of the poloidal component. Using the Davis-Chandrasekhar-Fermi method, we estimate a magnetic field strength of 3.5 mG. At the spatial scales accessible to ALMA, an analysis of the energy balance of the system indicates that gravity overwhelms turbulence, rotation, and the magnetic field. We show that high-mass star formation can occur in weakly magnetized environments, with gravity taking the dominant role., Comment: Accepted for publications in ApJL (12 pages, 5 figures, Appendix)

We observed the high-mass star-forming region G335.579-0.292 with the Atacama Large Millimeter/submillimeter Array (ALMA) at 226 GHz with an angular resolution of 0.3'' ($\sim 1000$ au resolution at the source distance). G335.579-0.292 hosts one of the most massive cores in the Galaxy (G335-MM1). The continuum emission shows that G335-MM1 fragments into at least five sources, while molecular line emission is detected in two of the continuum sources (ALMA1 and ALMA3). We found evidence of large and small scale infall in ALMA1 revealed by an inverse P-Cygni profile and the presence of a blue-shifted spot at the center of the first moment map of the CH$_3$CN emission. In addition, hot gas expansion in the innermost region is unveiled by a red-shifted spot in the first moment map of HDCO and (CH$_3$)$_2$CO (both with $E_u > 1100$ K). Our modeling reveals that this expansion motion originates close to the central source, likely due to reversal of the accretion flow induced by the expansion of the HII region, while infall and rotation motions originate in the outer regions. ALMA3 shows clear signs of rotation, with a rotation axis inclination with respect to the line of sight close to $90^\circ$, and a system mass (disk + star) in the range of 10-30 M$_\odot$., Comment: 15 pages, 12 figures, 3 tables. Accepted for publication in ApJ