Search

A complete understanding of the initial conditions of high-mass star formation and what processes determine multiplicity require the study of the magnetic field (B-field) in young, massive cores. Using ALMA 250 GHz polarization (0.3" = 1000 au) and ALMA 220 GHz high-angular resolution observations (0.05" = 160 au), we have performed a full energy analysis including the B-field at core scales and have assessed what influences the multiplicity inside a massive core previously believed to be in the prestellar phase. With 31 Msun, the G11.92 MM2 core has a young CS outflow with a dynamical time scale of a few thousand years. At high-resolution, the MM2 core fragments into a binary system with a projected separation of 505 au and a binary mass ratio of 1.14. Using the DCF method with an ADF analysis, we estimate in this core a B-field strength of 6.2 mG and a mass-to-flux ratio of 18. The MM2 core is strongly subvirialized with a virial parameter of 0.064, including the B-field. The high mass-to-flux ratio and low virial parameter indicate that this massive core is very likely undergoing runaway collapse, which is in direct contradiction with the core-accretion model. The MM2 core is embedded in a filament that has a velocity gradient consistent with infall. In line with clump-fed scenarios, the core can grow in mass at a rate of 1.9--5.6 x 10^-4 Msun/yr. In spite of the B-field having only a minor contribution to the total energy budget at core scales, it likely plays a more important role at smaller scales by setting the binary properties. Considering energy ratios and a fragmentation criterion at the core scale, the binary could have been formed by core fragmentation. The binary properties (separation and mass ratio), however, are also consistent with radiation-magnetohydrodynamic simulations with super-Alfvenic, supersonic (or sonic) turbulence that form binaries by disk fragmentation., Comment: Accepted for publications in ApJ (9 pages, 3 figures, Appendix)

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

Star formation primarily occurs in filaments where magnetic fields are expected to be dynamically important. The largest and densest filaments trace spiral structure within galaxies. Over a dozen of these dense ($\sim$10$^4$\,cm$^{-3}$) and long ($>$10\,pc) filaments have been found within the Milky Way, and they are often referred to as "bones." Until now, none of these bones have had their magnetic field resolved and mapped in their entirety. We introduce the SOFIA legacy project FIELDMAPS which has begun mapping $\sim$10 of these Milky Way bones using the HAWC+ instrument at 214\,$\mu$m and 18$\farcs$2 resolution. Here we present a first result from this survey on the $\sim$60\,pc long bone G47. Contrary to some studies of dense filaments in the Galactic plane, we find that the magnetic field is often not perpendicular to the spine (i.e., the center-line of the bone). Fields tend to be perpendicular in the densest areas of active star formation and more parallel or random in other areas. The average field is neither parallel or perpendicular to the Galactic plane nor the bone. The magnetic field strengths along the spine typically vary from $\sim$20 to $\sim$100\,$\mu$G. Magnetic fields tend to be strong enough to suppress collapse along much of the bone, but for areas that are most active in star formation, the fields are notably less able to resist gravitational collapse., Comment: Accepted to ApJL; typo in author list metadata corrected

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 present a spatio-kinematical analysis of the CO~($J$=2$\rightarrow$1) line emission, observed with the Atacama Large Millimter/submillimter Array (ALMA), of the outflow associated with the most massive core, ALMA1, in the 70 $\mu$m dark clump G010.991$-$00.082. The position-velocity (P-V) diagram of the molecular outflow exhibits a peculiar $\mathsf{S}$-shaped morphology that has not been seen in any other star forming region. We propose a spatio-kinematical model for the bipolar molecular outflow that consists of a decelerating high-velocity component surrounded by a slower component whose velocity increases with distance from the central source. The physical interpretation of the model is in terms of a jet that decelerates as it entrains material from the ambient medium, which has been predicted by calculations and numerical simulations of molecular outflows in the past. One side of the outflow is shorter and shows a stronger deceleration, suggesting that the medium through which the jet moves is significantly inhomogeneous. The age of the outflow is estimated to be $\tau$$\approx$1300 years, after correction for a mean inclination of the system of $\approx$57$^{\circ}$., Comment: 15 pages, 7 figures. Accepted for publication in ApJ

(Abridged) ASHES has been designed to systematically characterize the earliest stages and to constrain theories of high-mass star formation. A total of 12 massive (>500 $M_{\odot}$), cold (<15 K), 3.6-70 $\mu$m dark prestellar clump candidates, embedded in IRDCs, were carefully selected in the pilot survey to be observed with ALMA. We mosaiced each clump (~1 arcmin^2) in dust and line emission with the 12m/7m/TP arrays at 224 GHz, resulting in ~1.2" resolution (~4800 AU). As the first paper of the series, we concentrate on the dust emission to reveal the clump fragmentation. We detect 294 cores, from which 84 (29%) are categorized as protostellar based on outflow activity or 'warm core' line emission. The remaining 210 (71%) are considered prestellar core candidates. The number of detected cores is independent of the mass sensitivity range of the observations. On average, more massive clumps tend to form more cores. We find a large population of low-mass (<1 M) cores and no high-mass (>30 $M_{\odot}$) prestellar cores. The most massive prestellar core has a mass of 11 $M_{\odot}$. From the prestellar CMF, we derive a power law index of 1.17+-0.1, slightly shallower than Salpeter (1.35). We use the MST technique to characterize the separation between cores and their spatial distribution, and derive mass segregation ratios. While there is a range of core masses and separations detected in the sample, the mean separation and mass of cores are well explained by thermal fragmentation and are inconsistent with turbulent Jeans fragmentation. The core spatial distribution is well described by hierarchical subclustering rather than centrally peaked clustering. There is no conclusive evidence of mass segregation. We test several theoretical conditions, and conclude that overall, competitive accretion and global hierarchical collapse scenarios are favored over the turbulent core accretion scenario., Comment: Accepted for publication in ApJ (September 17, 2019). 37 pages, 19 figures, and 7 tables

We have observed complex organic molecules (COMs) and deuterated species toward a hot core/corino (HC) associated with the infrared dark cloud (IRDC) clump G34.43+00.24 MM3 with the Atacama Large Millimeter/submillimeter Array (ALMA). We have detected six normal-COMs (CH3OH, CH3CHO, CH3CH2CN, CH3OCH3, HCOOCH3, and NH2CHO), one deuterated-COM (CH2DCN), and two deuterated fundamental molecules (D2CO and DNC) toward G34.43+00.24 MM3 HC. None of these lines, except for CH3OH, are detected toward the shocked regions in our data, which suggests that COMs do not originate in shocks. The abundance of the COMs relative to CH3OH in G34.43+00.24 MM3 HC is found to be similar to those in high-mass hot cores, rather than those in hot corinos in low-mass star-forming regions. This result suggests that the physical conditions of the warm-up phase of G34.43+00.24 MM3 HC are similar to those of high-mass sources. On the other hand, the D2CO abundance relative to CH$_3$OH in G34.43+00.24 MM3 HC is higher than that of other hot cores, and seems to be comparable to that of hot corinos. The relatively high D2CO/CH3OH ratio of G34.43+00.24 MM3 HC implies a long cold starless phase of G34.43+00.24 MM3 HC., Comment: accepted to ApJ

We present dust column densities and dust temperatures for $\sim3000$ young high-mass molecular clumps from the Millimetre Astronomy Legacy Team 90 GHz (MALT90) survey, derived from adjusting single temperature dust emission models to the far-infrared intensity maps measured between 160 and 870 \micron\ from the Herschel/Hi-Gal and APEX/ATLASGAL surveys. We discuss the methodology employed in analyzing the data, calculating physical parameters, and estimating their uncertainties. The population average dust temperature of the clumps are: $16.8\pm0.2$ K for the clumps that do not exhibit mid-infrared signatures of star formation (Quiescent clumps), $18.6\pm0.2$ K for the clumps that display mid-infrared signatures of ongoing star formation but have not yet developed an HII region (Protostellar clumps), and $23.7\pm0.2$ and $28.1\pm0.3$ K for clumps associated with HII and photo-dissociation regions, respectively. These four groups exhibit large overlaps in their temperature distributions, with dispersions ranging between 4 and 6 K. The median of the peak column densities of the Protostellar clump population is $0.20\pm0.02$ gr cm$^{-2}$, which is about 50% higher compared to the median of the peak column densities associated with clumps in the other evolutionary stages. We compare the dust temperatures and column densities measured toward the center of the clumps with the mean values of each clump. We find that in the Quiescent clumps the dust temperature increases toward the outer regions and that they are associated with the shallowest column density profiles. In contrast, molecular clumps in the Protostellar or HII region phase have dust temperature gradients more consistent with internal heating and are associated with steeper column density profiles compared with the Quiescent clumps., Comment: 41 pages, 9 figures, accepted in ApJ

We have observed the clump G34.43+00.24 MM3 associated with an infrared dark cloud in DNC $J$=3--2, HN$^{13}$C $J$=3--2, and N$_2$H$^+$ $J$=3--2 with the Atacama Large Millimeter/submillimeter Array (ALMA). The N$_2$H$^+$ emission is found to be relatively weak near the hot core and the outflows, and its distribution is clearly anti-correlated with the CS emission. This result indicates that a young outflow is interacting with cold ambient gas. The HN$^{13}$C emission is compact and mostly emanates from the hot core, whereas the DNC emission is extended around the hot core. Thus, the DNC and HN$^{13}$C emission traces warm regions near the protostar differently. The DNC emission is stronger than the HN$^{13}$C emission toward most parts of this clump. The DNC/HNC abundance ratio averaged within a $15^{\prime\prime} \times 15^{\prime\prime}$ area around the phase center is higher than 0.06. This ratio is much higher than the value obtained by the previous single-dish observations of DNC and HN$^{13}$C $J$=1--0 ($\sim$0.003). It seems likely that the DNC and HNC emission observed with the single-dish telescope traces lower density envelopes, while that observed with ALMA traces higher density and highly deuterated regions. We have compared the observational results with chemical-model results in order to investigate the behavior of DNC and HNC in the dense cores. Taking these results into account, we suggest that the low DNC/HNC ratio in the high-mass sources obtained by the single-dish observations are at least partly due to the low filling factor of the high density regions., Comment: accepted to ApJ

We present the first results from the MALT-45 (Millimetre Astronomer's Legacy Team - 45 GHz) Galactic Plane survey. We have observed 5 square-degrees ($l = 330 - 335$, $b = \pm0.5$) for spectral lines in the 7 mm band (42-44 and 48-49 GHz), including $\text{CS}$ $(1-0)$, class I $\text{CH}_3\text{OH}$ masers in the $7(0,7)-6(1,6)$ $\text{A}^{+}$ transition and $\text{SiO}$ $(1-0)$ $v=0,1,2,3$. MALT-45 is the first unbiased, large-scale, sensitive spectral line survey in this frequency range. In this paper, we present data from the survey as well as a few intriguing results; rigorous analyses of these science cases are reserved for future publications. Across the survey region, we detected 77 class I $\text{CH}_3\text{OH}$ masers, of which 58 are new detections, along with many sites of thermal and maser $\text{SiO}$ emission and thermal $\text{CS}$. We found that 35 class I $\text{CH}_3\text{OH}$ masers were associated with the published locations of class II $\text{CH}_3\text{OH}$, $\text{H}_2\text{O}$ and $\text{OH}$ masers but 42 have no known masers within 60 arcsec. We compared the MALT-45 $\text{CS}$ with $\text{NH}_3$ (1,1) to reveal regions of $\text{CS}$ depletion and high opacity, as well as evolved star-forming regions with a high ratio of $\text{CS}$ to $\text{NH}_3$. All $\text{SiO}$ masers are new detections, and appear to be associated with evolved stars from the $\it{Spitzer}$ Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE). Generally, within $\text{SiO}$ regions of multiple vibrational modes, the intensity decreases as $v=1,2,3$, but there are a few exceptions where $v=2$ is stronger than $v=1$., Comment: Accepted by Monthly Notices of the Royal Astronomical Society on 2015 January 23. Appendices will be published with MNRAS

The Millimetre Astronomy Legacy Team 90 GHz (MALT90) survey has detected high-mass star-forming clumps with anomalous N$_2$H$^+$/HCO$^+$(1-0) integrated intensity ratios that are either unusually high ("N$_2$H$^+$ rich") or unusually low ("N$_2$H$^+$ poor"). With 3 mm observations from the Australia Telescope Compact Array (ATCA), we imaged two N$_2$H$^+$ rich clumps, G333.234-00.061 and G345.144-00.216, and two N$_2$H$^+$ poor clumps, G351.409+00.567 and G353.229+00.672. In these clumps, the N$_2$H$^+$ rich anomalies arise from extreme self-absorption of the HCO$^+$ line. G333.234-00.061 contains two of the most massive protostellar cores known with diameters of less than 0.1 pc, separated by a projected distance of only 0.12 pc. Unexpectedly, the higher mass core appears to be at an earlier evolutionary stage than the lower mass core, which may suggest that two different epochs of high-mass star formation can occur in close proximity. Through careful analysis of the ATCA observations and MALT90 clumps (including the G333, NGC 6334, and NGC 6357 star formation regions), we find that N$_2$H$^+$ poor anomalies arise at clump-scales and are caused by lower relative abundances of N$_2$H$^+$ due to the distinct chemistry of H II regions or photodissociation regions., Comment: 18 pages, 18 figures, accepted to ApJ

We have observed a molecular clump (MM3) associated with the infrared dark cloud G34.43+00.24 in the CH3OH J_K=9_{-1}--8_0 E, 5_0-4_0 E, and 5_{-1}-4_{-1} E lines at sub-arcsecond resolution by using the Atacama Large Millimeter/submillimeter Array. By comparing the CH3OH J_K=9_{-1}-8_0 E emission with the CH3OH 5_0-4_0 E and 5_{-1}-4_{-1} E emission, we have found that the CH3OH J_K=9_{-1}-8_0 E emission is masing. We have clearly shown that the CH3OH J_K=9_{-1}-8_0 masers arise from the post shocked gas in the interacting regions between the outflows and ambient dense gas. Toward the strongest peak of the CH3OH maser emission, SiO J=6-5 emission is very weak. This indicates that the CH3OH maser emission traces relatively old shocks or weak shocks., Comment: 15 pages, 4 figures, accepted to ApJL

We have used deep near-infrared observations with adaptive optics to discover a distributed population of low-mass protostars within the filamentary Infrared Dark Cloud G34.43+00.24. We use maps of dust emission at multiple wavelengths to determine the column density structure of the cloud. In combination with an empirically-verified model of the magnitude distribution of background stars, this column density map allows us to reliably determine overdensities of red sources that are due to embedded protostars in the cloud. We also identify protostars through their extended emission in K-band which comes from excited H2 in protostellar outflows or reflection nebulosity. We find a population of distributed low-mass protostars, suggesting that low-mass protostars may form earlier than, or contemporaneously with, high-mass protostars in such a filament. The low-mass protostellar population may also produce the narrow linewidth SiO emission observed in some clouds without high-mass protostars. Finally, we use a molecular line map of the cloud to determine the virial parameter per unit length along the filament and find that the highest mass protostars form in the most bound portion of the filament, as suggested by theoretical models., Comment: 15 pages, 16 figures. This version add internal hypertext links to references, figures and tables. Full quality version at http://www.astro.yale.edu/foster/papers/Foster_2014_DistributedSF.pdf

We have observed a cluster forming clump (MM3) associated with the infrared dark cloud G34.43+00.24 in the 1.3 mm continuum and the CH3OH, CS, 13CS, SiO, CH3CH2CN, and HCOOCH3 lines with the Atacama Large Millimeter/submillimeter Array and in K-band with the Keck telescope. We have found a young outflow toward the center of this clump in the SiO, CS, and CH3OH lines. This outflow is likely driven by a protostar embedded in a hot core, which is traced by the CH3CH2CN, HCOOCH3, 13CS, and high excitation CH3OH lines. The size of the hot core is about 800 x 300 AU in spite of its low mass (<1.1 M_sun), suggesting a high accretion rate or the presence of multiple star system harboring a few hot corinos. The outflow is highly collimated, and the dynamical timescale is estimated to be less than 740 yr. In addition, we have also detected extended emission of SiO, CS, and CH3OH, which is not associated with the hot core and the outflow. This emission may be related to past star formation activity in the clump. Although G34.43+00.24 MM3 is surrounded by a dark feature in infrared, it has already experienced active formation of low-mass stars in an early stage of clump evolution., Comment: 15 pages, 4 figures, accepted to ApJ

We have observed the IRDC G028.23-00.19 at 3.3 mm using CARMA. In its center, the IRDC hosts one of the most massive (~1520 Msun) quiescent, cold (12 K) clumps known (MM1). The low temperature, high NH_2D abundance, narrow molecular line widths, and absence of embedded IR sources (from 3.6 to 70 mu) indicate that the clump is likely prestellar. Strong SiO emission with broad line widths (6-9 km s^-1) and high abundances (0.8-4 x 10^-9) is detected in the northern and southern regions of the IRDC, unassociated with MM1. We suggest that SiO is released to the gas phase from the dust grains through shocks produced by outflows from undetected intermediate-mass stars or clusters of low-mass stars deeply embedded in the IRDC. A weaker SiO component with narrow line widths (~2 km s^-1) and low abundances (4.3 x 10^-11) is detected in the center-west region, consistent with either a "subcloud-subcloud" collision or an unresolved population of a few low-mass stars. We report widespread CH_3OH emission throughout the whole IRDC and the first detection of extended narrow methanol emission (~2 km s^-1) in a cold, massive prestellar clump (MM1). We suggest that the most likely mechanism releasing methanol into the gas phase in such a cold region is the exothermicity of grain-surface reactions. HN^13C reveals that the IRDC is actually composed of two distinct substructures ("subclouds") separated in velocity space by ~1.4 km s^-1. The narrow SiO component arises where the subclouds overlap. The spatial distribution of C_2H resembles that of NH_2D, which suggests that C_2H also traces cold gas in this IRDC., Comment: Accepted to ApJ. 16 pages

We present high spatial resolution mid-IR images of the ring of UCHII regions in W49A obtained at Gemini North, allowing us to identify the driving source of its powerful H2O maser outflow. These data also confirm our previous report that several radio sources in the ring are undetected in the mid-IR because they are embedded deep inside the cloud core. We locate the source of the water maser outflow at the position of the compact mid-IR peak of source G (source G:IRS1). This IR source is not coincident with any identified compact radio continuum source, but is coincident with a hot molecular core, so we propose that G:IRS1 is a hot core driving an outflow analogous to the wide-angle bipolar outflow in OMC-1. G:IRS1 is at the origin of a larger bipolar cavity and CO outflow. The water maser outflow is orthogonal to the bipolar CO cavity, so the masers probably reside near its waist in the cavity walls. Models of the IR emission require a massive protostar of 45Msun, 3e5Lsun, and an effective envelope accretion rate of 1e-3Msun/yr. Feedback from the central star could potentially drive the H2O maser outflow, but it has insufficient radiative momentum to have driven the large-scale CO outflow, requiring that this massive star had an active accretion disk over the past 10^4 yr. Combined with the spatialy resolved morphology in IR images, G:IRS1 in W49 provides compelling evidence for a massive protostar that formed by accreting from a disk, accompanied by a bipolar outflow., Comment: 14 pages, MNRAS accepted