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

Recent (sub)millimeter polarization observations of protoplanetary disks reveal toroidally aligned, effectively prolate dust grains large enough (at least ~100 $\mu$m) to efficiently scatter millimeter light. The alignment mechanism for these grains remains unclear. We explore the possibility that gas drag aligns grains through gas-dust relative motion when the grain's center of mass is offset from its geometric center, analogous to a badminton birdie's alignment in flight. A simple grain model of two non-identical spheres illustrates how a grain undergoes damped oscillations from flow-induced restoring torques which align its geometric center in the flow direction relative to its center of mass. Assuming specular reflection and subsonic flow, we derive an analytical equation of motion for spheroids where the center of mass can be shifted away from the spheroid's geometric center. We show that a prolate or an oblate grain can be aligned with the long axis parallel to the gas flow when the center of mass is shifted along that axis. Both scenarios can explain the required effectively prolate grains inferred from observations. Application to a simple disk model shows that the alignment timescales are shorter than or comparable to the orbital time. The grain alignment direction in a disk depends on the disk (sub-)structure and grain Stokes number (St) with azimuthal alignment for large St grains in sub-Keplerian smooth gas disks and for small St grains near the gas pressure extrema, such as rings and gaps., Comment: 21 pages, 12 figures, submitted to MNRAS

We use the H$41\alpha$ recombination line to create templates of the millimeter free-free emission in the ALMA-IMF continuum maps, which allows to separate it from dust emission. This method complements spectral-index information and extrapolation from centimeter wavelength maps. We use the derived maps to estimate the properties of up to 34 HII regions across the ALMA-IMF protoclusters. The hydrogen ionizing-photon rate $Q_0$ and spectral types follow the evolutionary trend proposed by Motte et al. The youngest protoclusters lack detectable ionized gas, followed by protoclusters with increasing numbers of OB stars. The total $Q_0$ increases from $\sim 10^{45}$ s$^{-1}$ to $> 10^{49}$ s$^{-1}$. We used the adjacent He$41\alpha$ line to measure the relative number abundances of helium, finding values consistent with the Galactic interstellar medium, although a few outliers are discussed. A search for sites of maser amplification of the H$41\alpha$ line returned negative results. We looked for possible correlations between the electron densities ($n_e$), emission measures (EM), and $Q_0$ with HII region size $D$. The latter are the better correlated, with $Q_0 \propto D^{2.49\pm0.18}$. This favors interpretations where smaller ultracompact HII regions are not necessarily the less dynamically evolved versions of larger ones, but rather are ionized by less massive stars. Moderate correlations were found between dynamical width $\Delta V_\mathrm{dyn}$ with $D$ and $Q_0$. $\Delta V_\mathrm{dyn}$ increases from about one to two times the ionized-gas sound speed. Finally, an outlier HII region south of W43-MM2 is discussed. We suggest that this source could harbor an embedded stellar or disk wind., Comment: Accepted to AAS Journals (ApJS). Data available at https://dataverse.harvard.edu/dataverse/h41a-freefree

Millimeter and sub-millimeter observations of continuum linear dust polarization provide insight into dust grain growth in protoplanetary disks, which are the progenitors of planetary systems. We present the results of the first survey of dust polarization in protoplanetary disks at 870 $\mu$m and 3 mm. We find that protoplanetary disks in the same molecular cloud at similar evolutionary stages can exhibit different correlations between observing wavelength and polarization morphology and fraction. We explore possible origins for these differences in polarization, including differences in dust populations and protostar properties. For RY Tau and MWC 480, which are consistent with scattering at both wavelengths, we present models of the scattering polarization from several dust grain size distributions. These models aim to reproduce two features of the observational results for these disks: (1) both disks have an observable degree of polarization at both wavelengths and (2) the polarization fraction is higher at 3 mm than at 870 $\mu$m in the centers of the disks. For both disks, these features can be reproduced by a power-law distribution of spherical dust grains with a maximum radius of 200 $\mu$m and high optical depth. In MWC 480, we can also reproduce features (1) and (2) with a model containing large grains ($a_{max}$ = 490 $\mu$m ) near the disk midplane and small grains ($a_{max}$ = 140 $\mu$m) above and below the midplane., Comment: 22 pages, 12 figures

We investigate the crescent-shaped dust trap in the transition disk, Oph IRS 48, using well-resolved (sub)millimeter polarimetric observations at ALMA Band 7 (870 $\mu$m). The dust polarization map reveals patterns consistent with dust scattering-induced polarization. There is a relative displacement between the polarized flux and the total flux, which holds the key to understanding the dust scale heights in this system. We model the polarization observations, focusing on the effects of dust scale heights. We find that the interplay between the inclination-induced polarization and the polarization arising from radiation anisotropy in the crescent determines the observed polarization; the anisotropy is controlled by the dust optical depth along the midplane, which is, in turn, determined by the dust scale height in the vertical direction. We find that the dust grains can neither be completely settled nor well mixed with the gas. The completely settled case produces little radial displacement between the total and polarized flux, while the well-mixed case produces an azimuthal pattern in the outer (radial) edge of the crescent that is not observed. Our best model has a gas-to-dust scale height ratio of 2, and can reproduce both the radial displacement and the azimuthal displacement between the total and polarized flux. We infer an effective turbulence $\alpha$ parameter of approximately $0.0001-0.005$. The scattering-induced polarization provides insight into a turbulent vortex with a moderate level of dust settling in the IRS 48 system, which is hard to achieve otherwise., Comment: 9 pages in the main text, 5 figures, accepted by ApJ

We present 870 um polarimetric observations toward 61 protostars in the Orion molecular clouds, with ~400 au (1") resolution using the Atacama Large Millimeter/submillimeter Array. We successfully detect dust polarization and outflow emission in 56 protostars, in 16 of them the polarization is likely produced by self-scattering. Self-scattering signatures are seen in several Class 0 sources, suggesting that grain growth appears to be significant in disks at earlier protostellar phases. For the rest of the protostars, the dust polarization traces the magnetic field, whose morphology can be approximately classified into three categories: standard-hourglass, rotated-hourglass (with its axis perpendicular to outflow), and spiral-like morphology. 40.0% (+-3.0%) of the protostars exhibit a mean magnetic field direction approximately perpendicular to the outflow on several 100--1000 au scales. However, in the remaining sample, this relative orientation appears to be random, probably due to the complex set of morphologies observed. Furthermore, we classify the protostars into three types based on the C17O (3--2) velocity envelope's gradient: perpendicular to outflow, non-perpendicular to outflow, and unresolved gradient (<1.0~km/s/arcsec). In protostars with a velocity gradient perpendicular to outflow, the magnetic field lines are preferentially perpendicular to outflow, most of them exhibit a rotated hourglass morphology, suggesting that the magnetic field has been overwhelmed by gravity and angular momentum. Spiral-like magnetic fields are associated with envelopes having large velocity gradients, indicating that the rotation motions are strong enough to twist the field lines. All of the protostars with a standard-hourglass field morphology show no significant velocity gradient due to the strong magnetic braking., Comment: 12 pages, 5 figures (16 pages and 7 figures for the supplementary materials) Accepted for publication on ApJ Letters

Polarized (sub)millimeter emission from dust grains in circumstellar disks was initially thought to be due to grains aligned with the magnetic field. However, higher resolution multi-wavelength observations along with improved models found that this polarization is dominated by self-scattering at shorter wavelengths (e.g., 870 $\mu$m) and by grains aligned with something other than magnetic fields at longer wavelengths (e.g., 3 mm). Nevertheless, the polarization signal is expected to depend on the underlying substructure, and observations hitherto have been unable to resolve polarization in multiple rings and gaps. HL Tau, a protoplanetary disk located 147.3 $\pm$ 0.5 pc away, is the brightest Class I or Class II disk at millimeter/submillimeter wavelengths. Here we show deep, high-resolution 870 $\mu$m polarization observations of HL Tau, resolving polarization in both the rings and gaps. We find that the gaps have polarization angles with a significant azimuthal component and a higher polarization fraction than the rings. Our models show that the disk polarization is due to both scattering and emission from aligned effectively prolate grains. The intrinsic polarization of aligned dust grains is likely over 10%, which is much higher than what was expected in low resolution observations (~1%). Asymmetries and dust features are revealed in the polarization observations that are not seen in non-polarimetric observations., Comment: Published in Nature

We present the combination of ALMA-IMF and single-dish continuum images from the Mustang-2 Galactic Plane Survey (MGPS90) at 3 millimeters and the Bolocam Galactic Plane Survey (BGPS) at 1 millimeter. Six and ten out of the fiffteen ALMA-IMF fields are combined with MGPS90 and BGPS, respectively. The combination is made via the feathering technique. We used the dendrogram algorithm throughout the combined images, and performed further analysis in the six fields with combination in both bands (G012.80, W43-MM1, W43-MM2, W43-MM3, W51-E, W51-IRS2). In these fields, we calculated spectral index maps and used them to separate regions dominated by dust or free-free emission, and then performed further structural analysis. We report the basic physical parameters of the dust-dominated (column densities, masses) and ionized (emission measures, hydrogen ionization photon rates) structures. We also searched for multi-scale relations in the dust-dominated structures across the analyzed fields, finding that the fraction of mass in dendrogram leaves (which we label as "Leaf Mass Eficiency", LME) as a function of molecular gas column density follows a similar trend: a rapid, exponential-like growth, with maximum values approaching 100% in most cases. The observed behaviour of the LME with gas column is tentatively interpreted as an indicator of large star formation activity within the ALMA-IMF protoclusters. W51-E and G012.80 stand out as cases with comparatively large and reduced potential for further star formation, respectively., Comment: Accepted to The Astrophysical Journal Supplement

In the last years there has been a substantial increase in the number of the reported massive and luminous star-forming regions with related explosive outflows thanks to the superb sensitivity and angular resolution provided by the new radio, infrared, and optical facilities. Here, we report one more explosive outflow related with the massive and bright star-forming region IRAS 12326$-$6245 using Band 6 sensitive and high angular resolution ($\sim$0.2$"$) Atacama Large Millimeter/Submillimeter Array (ALMA) observations. We find over 10 molecular and collimated well-defined streamers, with Hubble-Lemaitre like expansion motions, and pointing right to the center of a dusty and molecular shell (reported for the first time here) localized in the northern part of the UCHII region known as G301.1A. The estimated kinematic age, and energy for the explosion are $\sim$700 yrs, and 10$^{48}$ erg, respectively. Taking into account the recently reported explosive outflows together with IRAS 12326$-$6245, we estimate an event rate of once every 90 yr in our Galaxy, similar to the formation rate of massive stars., Comment: Accepted to The Astrophysical Journal Letters

Polarization is a unique tool to study the properties of dust grains of protoplanetary disks and detail the initial conditions of planet formation. Polarization around HL Tau was previously imaged using the Atacama Large Millimeter/submillimeter Array (ALMA) at Bands 3 (3.1 mm), 6 (1.3 mm), and 7 (0.87 mm), showing that the polarization orientation changes across wavelength $\lambda$. The polarization morphology at Band 7 is predominantly parallel to the disk minor axis but appears azimuthally oriented at Band 3, with the morphology at Band 6 in between the two. We present new ~0.2" (29 au) polarization observations at Q-Band (7.0 mm) using the Karl G. Jansky Very Large Array (VLA) and at Bands 4 (2.1 mm), 5 (1.5 mm), and 7 using ALMA, consolidating HL Tau's position as the protoplanetary disk with the most complete wavelength coverage in dust polarization. The polarization patterns at Bands 4 and 5 continue to follow the morphological transition with wavelength previously identified in Bands 3, 6, and 7. Based on the azimuthal variation, we decompose the polarization into contributions from scattering ($s$) and thermal emission ($t$). We find that $s$ decreases slowly with increasing $\lambda$, and $t$ increases more rapidly with $\lambda$ which are expected from optical depth effects of toroidally aligned, scattering prolate grains. The relatively weak $\lambda$ dependence of $s$ is consistent with large, porous grains. The sparse polarization detections from the Q-band image are also consistent with toroidally aligned prolate grains., Comment: 22 pages, 13 figures, submitted to MNRAS

Crescent-shaped structures in transition disks hold the key to studying the putative companions to the central stars. The dust dynamics, especially that of different grain sizes, is important to understanding the role of pressure bumps in planet formation. In this work, we present deep dust continuum observation with high resolution towards the Oph IRS 48 system. For the first time, we are able to significantly trace and detect emission along $95\%$ of the ring crossing the crescent-shaped structure. The ring is highly eccentric with an eccentricity of $0.27$. The flux density contrast between the peak of the flux and its counter part along the ring is $\sim 270$. In addition, we detect a compact emission toward the central star. If the emission is an inner circumstellar disk inside the cavity, it has a radius of at most a couple of astronomical units with a dust mass of $1.5\times 10^{-8}\rm\, M_\odot$, or $0.005\rm\, M_\oplus$. We also discuss the implications of the potential eccentric orbit on the proper motion of the crescent, the putative secondary companion, and the asymmetry in velocity maps., Comment: 9 pages, 5 figures, accepted for publication on AAS Journals

We present 0.9 mm continuum and CO (3-2) line emission observations retrieved from the Atacama Large Millimeter/submillimeter Array (ALMA) archive toward the high-mass star formation region IRAS 16076-5134. We identify fourteen dense cores with masses between 0.3 to 22 M$_{\odot}$. We find an ensemble of filament-like CO (3-2) ejections from -62 to +83 km s$^{-1}$ that appear to arise radially from a common central position, close to the dense core MM8. The ensemble of filaments, has a quasi-isotropic distribution in the plane of the sky. The radial velocity of several filaments follow a linear velocity gradient, incresing from a common origin. Considering the whole ensemble of filaments, we estimate its total mass to be 138 and 216 M$_{\odot}$ from its CO emission, for 70 K and 140 K respectively. Also, assuming a constant velocity expansion of the filaments (of 83 km s$^{-1}$) we estimate the dynamical age of the outflowing material (3500 years), its momentum (~10$^{4}$ M$_{\odot}$ km s$^{-1}$) and its kinetic energy (~10$^{48-49}$ erg). The morphology and kinematics presented by the filaments suggest the presence of a dispersal outflow with explosive characteristics in IRAS 16076-5134. In addition, we make a raw estimate of the lower limit of the frequency rate of the explosive dispersal outflows in the Galaxy (one every 110 years) considering constant star formation rate and efficiency with respect to the galactocentric radius of the Galaxy. This may imply a comparable rate of dispersal outflows and supernovae (approximately one every 50 years), which may be important for the energy budget of the Interstellar Medium and the link between dispersal outflows and high-mass star formation., Comment: 13 pages, 8 figures. ApJ accepted

The size of dust grains, $a$, is key to the physical and chemical processes in circumstellar disks, but observational constraints of grain size remain challenging. (Sub)millimeter continuum observations often show a percent-level polarization parallel to the disk minor axis, which is generally attributed to scattering by $\sim 100\mu$m-sized spherical grains (with a size parameter $x \equiv 2\pi a / \lambda < 1$, where $\lambda$ is the wavelength). Larger spherical grains (with $x$ greater than unity) would produce opposite polarization direction. However, the inferred size is in tension with the opacity index $\beta$ that points to larger mm/cm-sized grains. We investigate the scattering-produced polarization by large irregular grains with a range of $x$ greater than unity with optical properties obtained from laboratory experiments. Using the radiation transfer code, RADMC-3D, we find that large irregular grains still produce polarization parallel to the disk minor axis. If the original forsterite refractive index in the optical is adopted, then all samples can produce the typically observed level of polarization. Accounting for the more commonly adopted refractive index using the DSHARP dust model, only grains with $x$ of several (corresponding to $\sim$mm-sized grains) can reach the same polarization level. Our results suggest that grains in disks can have sizes in the millimeter regime, which may alleviate the tension between the grain sizes inferred from scattering and other means. Additionally, if large irregular grains are not settled to the midplane, their strong forward scattering can produce asymmetries between the near and far side of an inclined disk, which can be used to infer their presence., Comment: 14 pages, 9 figures, published by MNRAS

The nature and origin of the molecular gas component located in the circumstellar vicinity of $\eta$ Carinae are still far from being completely understood. Here, we present Atacama Large Millimeter/Submillimeter Array (ALMA) CO(3$-$2) observations with a high angular resolution ($\sim$0.15$''$), and a great sensitivity that are employed to reveal the origin of this component in $\eta$ Carinae. These observations reveal much higher velocity ($-$300 to $+$270 km s$^{-1}$) blue and redshifted molecular thermal emission than previously reported, which we associate with the lobes of the Homunculus Nebula, and that delineates very well the innermost contours of the red- and blue-shifted lobes likely due by limb brightening. The inner contour of the redshifted emission was proposed to be a {\it disrupted torus}, but here we revealed that it is at least part of the molecular emission originated from the lobes and/or the expanding equatorial skirt. On the other hand, closer to systemic velocities ($\pm$100 km s$^{-1}$), the CO molecular gas traces an inner butterfly-shaped structure that is also revealed at NIR and MIR wavelengths as the region in which the shielded dust resides. The location and kinematics of the molecular component indicate that this material has formed after the different eruptions of $\eta$ Carinae., Comment: Accepted to Astrophysical Journal. 3D animation -> https://drive.google.com/file/d/1MRWaTkij-Lf_uzTfamJfHeGvgzjDe6Oq/view?usp=sharing

Telescopes are now able to resolve dust polarization across circumstellar disks at multiple wavelengths, allowing the study of the polarization spectrum. Most disks show clear evidence of dust scattering through their unidirectional polarization pattern typically at the shorter wavelength of $\sim 870 \mu$m. However, certain disks show an elliptical pattern at $\sim 3$mm, which is likely due to aligned grains. With HL Tau, its polarization pattern at $\sim 1.3$mm shows a transition between the two patterns making it the first example to reveal such transition. We use the T-matrix method to model elongated dust grains and properly treat scattering of aligned non-spherical grains with a plane-parallel slab model. We demonstrate that a change in optical depth can naturally explain the polarization transition of HL Tau. At low optical depths, the thermal polarization dominates, while at high optical depths, dichroic extinction effectively takes out the thermal polarization and scattering polarization dominates. Motivated by results from the plane-parallel slab, we develop a simple technique to disentangle thermal polarization of the aligned grains $T_{0}$ and polarization due to scattering $S$ using the azimuthal variation of the polarization fraction. We find that, with increasing wavelength, the fractional polarization spectrum of the scattering component $S$ decreases, while the thermal component $T_{0}$ increases, which is expected since the optical depth decreases. We find several other sources similar to HL Tau that can be explained by azimuthally aligned scattering prolate grains when including optical depth effects. In addition, we explore how spirally aligned grains with scattering can appear in polarization images., Comment: 27 pages, 21 figures, accepted by MNRAS

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

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)

While magnetic fields likely play an important role in driving the evolution of protoplanetary disks through angular momentum transport, observational evidence of magnetic fields has only been found in a small number of disks. Although dust continuum linear polarization has been detected in an increasing number of disks, its pattern is more consistent with that from dust scattering than from magnetically aligned grains in the vast majority of cases. Continuum linear polarization from dust grains aligned to a magnetic field can reveal information about the magnetic field's direction, but not its strength. On the other hand, observations of circular polarization in molecular lines produced by Zeeman splitting offer a direct measure of the line-of-sight magnetic field strength in disks. We present upper limits on the net toroidal and vertical magnetic field strengths in the protoplanetary disk AS 209 derived from Zeeman splitting observations of the CN 2-1 line. The 3$\sigma$ upper limit on the net line-of-sight magnetic field strength in AS 209 is 5.0 mG on the redshifted side of the disk and 4.2 mG on the blueshifted side of the disk. Given the disk's inclination angle, we set a 3$\sigma$ upper limit on the net toroidal magnetic field strength of 8.7 and 7.3 mG for the red and blue sides of the disk, respectively, and 6.2 and 5.2 mG on the net vertical magnetic field on the red and blue sides of the disk. If magnetic disk winds are a significant mechanism of angular momentum transport in the disk, magnetic fields of a strength close to the upper limits would be sufficient to drive accretion at the rate previously inferred for regions near the protostar., Comment: June 22, 2021: Corrected disk size scale bar in Figures 1, 2 and 5. Corrected percent polarization scale bar label in Figure 5

The explosive molecular outflow detected decades ago in the Orion BN/KL region of massive star formation was considered to be a bizarre event. This belief was strengthened by the non detection of similar cases over the years with the only exception of the marginal case of DR21. Here, we confim a similar explosive outflow associated with the UCH$_{\rm II}$ region G5.89$-$0.39 that indicates that this phenomenon is not unique to Orion or DR21. Sensitive and high angular resolution ($\sim$ 0.1$''$) ALMA CO(2$-$1) and SiO(5$-$4) observations show that the molecular outflow in the massive star forming region G5.89$-$0.39 is indeed an explosive outflow with an age of about 1000 yrs and a liberated kinetic energy of 10$^{46-49}$ erg. Our new CO(2$-$1) ALMA observations revealed over 30 molecular filaments, with Hubble-like expansion motions, pointing to the center of UCH$_{\rm II}$ region. In addition, the SiO(5$-$4) observations reveal warmer and strong shocks very close to the origin of the explosion, confirming the true nature of the flow. A simple estimation for the occurrence of these explosive events during the formation of the massive stars indicates an event rate of once every $\sim$100 yrs, which is close to the supernovae rate., Comment: Accepted by ApJ Letters. The animation is located here: https://drive.google.com/file/d/1a1c9JMwm0Mf8CxlfuEzyh6XmCOsgOyMA/view?usp=sharing

Annular structures, or rings and gaps, in disks around pre-main sequence stars have been detected in abundance towards Class II objects ~1,000,000 years in age. These structures are often interpreted as evidence of planet formation, with planet-mass bodies carving rings and gaps in the disk. This implies that planet formation may already be underway in even younger disks in the Class I phase, when the protostar is still embedded in a larger scale dense envelope of gas and dust. While younger disks likely play an important role in the onset of planet formation, only within the past decade have detailed properties of disks in the youngest star-forming phases begun to be observed. Here we present 1.3 mm dust emission observations with 5 au resolution that show four annular substructures in the disk of the young (<500,000 years) protostar IRS 63. IRS 63, a single Class I source located in the nearby Ophiuchus molecular cloud (at a distance of 144 pc), is one of the brightest Class I protostars at (sub)millimeter wavelengths that also has a relatively large disk (>50 au). Multiple annular substructures observed towards disks at young times can act as an early foothold for dust grain growth, a prerequisite of planet formation. Whether planets already exist or not in the disk of IRS 63, it is clear that the planet formation process begins in the young protostellar phases, earlier than predicted by current planet formation theories., Comment: Accepted for publication in Nature, appearing in the 8 Oct 2020 issue, DOI: 10.1038/s41586-020-2779-6

Magnetic fields are expected to play an important role in accretion processes for circumstellar disks. Measuring the magnetic field morphology is difficult, especially since polarimetric (sub)millimeter continuum observations may not trace fields in most disks. The Goldreich-Kylafis (GK) effect suggests that line polarization is perpendicular or parallel to the magnetic field direction. We attempt to observe CO(2-1), $^{13}$CO(2-1), and C$^{18}$O(2-1) line polarization toward HD 142527 and IM Lup, which are large, bright protoplanetary disks. We use spatial averaging and spectral integration to search for signals in both disks, and detect a potential CO(2-1) Stokes $Q$ signal toward both disks. The total CO(2-1) polarization fractions are 1.57 $\pm$ 0.18% and 1.01 $\pm$ 0.10% for HD 142527 and IM Lup, respectively. Our Monte Carlo simulations indicate that these signals are marginal. We also stack Stokes parameters based on the Keplerian rotation, but no signal was found. Across the disk traced by dust of HD 142527, the 3$\sigma$ upper limits for $P_{\text{frac}}$ at 0.5$^{\prime\prime}$ ($\sim$80 au) resolution are typically less than 3% for CO(2-1) and $^{13}$CO(2-1) and 4% for C$^{18}$O(2-1). For IM Lup, the 3$\sigma$ upper limits for these three lines are typically less than 3%, 4%, and 12%, respectively. Upper limits based on our stacking technique are up to a factor of $\sim$10 lower, though stacking areas can potentially average out small-scale polarization structure. We also compare our continuum polarization at 1.3 mm to observations at 870 $\mu$m from previous studies. The polarization in the northern dust trap of HD 142527 shows a significant change in morphology and an increase in $P_{\text{frac}}$ as compared to 870 $\mu$m. For IM Lup, the 1.3 mm polarization may be more azimuthal and has a higher $P_{\text{frac}}$ than at 870 $\mu$m., Comment: Accepted to ApJ

We present sensitive and high angular resolution ($\sim$0.2-0.3$''$) (sub)millimeter (230 and 345 GHz) continuum and CO(2$-$1)/CO(3$-$2) line archive observations of the disk star system in UX Tauri carried out with ALMA (The Atacama Large Millimeter/Submillimeter Array). These observations reveal the gas and dusty disk surrounding the young star UX Tauri A with a large signal-to-noise ratio ($>$400 in the continuum and $>$50 in the line), and for the first time is detected the molecular gas emission associated with the disk of UX Tauri C (with a size for the disk of $<$56 au). No (sub)millimeter continuum emission is detected at 5$\sigma$-level (0.2 mJy at 0.85 mm) associated with UX Tauri C. For the component UX Tauri C, we estimate a dust disk mass of $\leq$ 0.05 M$_\oplus$. Additionally, we report a strong tidal disk interaction between both disks UX Tauri A/C, separated 360 au in projected distance. The CO line observations reveal marked spiral arms in the disk of UX Tauri A and an extended redshifted stream of gas associated with the UX Tauri C disk. No spiral arms are observed in the dust continuum emission of UX Tauri A. Assuming a Keplerian rotation we estimate the enclosed masses (disk$+$star) from their radial velocities in 1.4 $\pm$ 0.6 M$_\odot$ for UX Tauri A, and 70 $\pm$ 30 / $\sin i$ Jupiter masses for UX Tauri C (the latter coincides with the mass upper limit value for a brown dwarf). The observational evidence presented here lead us to propose that UX Tauri C is having a close approach of a possible wide, evolving and eccentric orbit around the disk of UX Tauri A causing the formation of spiral arms and the stream of molecular gas falling towards UX Tauri C., Comment: Accepted by The Astrophysical Journal

The explosive outflows are a newly-discovered family of molecular outflows associated with high-mass star forming regions. Such energetic events are possibly powered by the release of gravitational energy related with the formation of a (proto)stellar merger or a close stellar binary. Here, we present sensitive and high angular resolution observations (0.85$''$) archival CO(J=3-2) observations carried out with the Submillimeter Array (SMA) of the high-mass star forming region G5.89$-$0.39 that reveal the possible presence of an explosive outflow. We find six well-defined and narrow straight filament-like ejections pointing back approximately to the center of an expanding molecular and ionized shell located at the center of this region. These high velocity ($-$120 to $+$100 km s$^{-1}$) filaments follow a Hubble-like velocity law with the radial velocities increasing with the projected distance. The estimated kinematical age of the filaments is about of 1000 yrs, a value similar to the dynamical age found for the expanding ionized shell. G5.89 is the thus the third explosive outflow reported in the galaxy (together with Orion BN-KL and DR21) and argues in favor of the idea that this is a frequent phenomenon. In particular, explosive outflows, in conjunction with runaway stars, demonstrate that dynamical interactions in such groups are a very important ingredient in star formation., Comment: Accepted by MNRAS Letters

Polarized dust emission outside of disks reveal the magnetic field morphology of molecular clouds. Within disks, however, polarized dust emission can arise from very different mechanisms (e.g., self-scattering), and each of them are useful for constraining physical properties in the disk. For example, these mechanisms allow us to constrain the disk grain size distributions and grain/disk geometries, independent from current methods of measuring these parameters. To accurately model these features and disentangle the various polarization mechanisms, multiwavelength observations at very high resolution and sensitivity are required. With significant upgrades to current interferometric facilities, we can understand how grains evolve in disks during the planet formation process., Comment: Science white paper submitted to the Astro2020 Decadal Survey

During the last decades, a great interest has emerged to know if even the most massive stars in our galaxy (namely the spectral O-type stars) are formed in a similar manner as the low- and intermediate-mass stars, that is, through the presence of accreting disks and powerful outflows. Here, using sensitive observations of the Atacama Large Millimeter/Submillimeter Array (ALMA), we report a resolved Keplerian disk (with fifteen synthesized beams across its major axis) surrounding the deeply embedded O-type protostar IRAS16547$-$4247. The disk shows some asymmetries that could arise because of the disk is unstable and fragmenting or maybe because of different excitation conditions within the disk. The enclosed mass estimated from the disk Keplerian radial velocities is 25$\pm$3 M$_\odot$. The molecular disk is at the base of an ionized thermal radio jet and is approximately perpendicular to the jet axis orientation. We additionally find the existence of a binary system of compact dusty objects at the center of the accreting disk, which indicates the possible formation of an O-type star and a companion of lower mass. This is not surprising due to the high binary fraction reported in massive stars. Subtracting the contribution of the dusty disk plus the envelope and the companion, we estimated a mass of 20 M$_\odot$ for the central star., Comment: Accepted by The Astrophysical Journal

We present high-sensitivity ($\sigma_I\sim0.2-0.5$ mJy, $\sigma_{QU}\sim0.05$ mJy), high-resolution (0.12-0.2 arcsec) observations of polarized 872 $\mu$m dust emission from the young multiple system VLA 1623 in $\rho$ Ophiuchus and the protostar L1527 in Taurus. We detect the circumstellar material of VLA 1623A, the extended Keplerian disk surrounding VLA 1623A which we call VLA 1623CBdisk, VLA 1623B, VLA 1623W, and L1527 strongly in the polarized emission, at the $\sim1-3\%$ level. We spatially resolve VLA 1623A into two sources, VLA 1623Aa and VLA 1623Ab, separated by $\sim30$ au and located within a cavity of radius $\sim50$ au within the circumbinary Keplerian disk, as well as the edge-on disk of VLA 1623W. The polarization angle of the emission is uniform across each protostellar source and nearly coincides with each disk's minor axis. The offsets between the minor axis and the polarization angle are not uniformly distributed at the $P\lesssim2\times10^{-4}$ level. The circumbinary disk surrounding VLA 1623Aab is azimuthally symmetrically polarized. Each compact source's emission is partially optically thick ($\tau\gtrsim1$) at 872 $\mu$m, complicating interpretations of polarization involving aligned grains. We find evidence against alignment by radiative flux in each source, particularly in the edge-on VLA 1623W and L1527. We detect astrometric offsets between the polarized emission and the total intensity in VLA 1623Aa, VLA 1623Ab, and VLA 1623B, as predicted if self-scattering in the optically thick limit operates. We conclude that self-scattering is likely responsible for disk-scale polarization at 872 $\mu$m in these systems., Comment: 22 pages, 7 figures, accepted to Astrophysical Journal

Aims: We investigated the masses and spatial distributions of pre-stellar and protostellar candidates in the young, low-mass star forming region Serpens South, where active star formation is known to occur along a predominant filamentary structure. Previous observations used to study these distributions have been limited by two important observational factors: (1) sensitivity limits that leave the lowest-mass sources undetected, or (2) resolution limits that cannot distinguish binaries and/or cluster members in close proximity. Methods: Recent millimeter-wavelength interferometry observations can now uncover faint and/or compact sources in order to study a more complete population of protostars, especially in nearby ($D<500$ pc) clusters. Here we present ALMA observations of 1 mm (Band 6) continuum in a $3 \times 2$ arcminutes region at the center of Serpens South. Our angular resolution of $\sim1$ arcsec is equivalent to $\sim400$ au, corresponding to scales of envelopes and/or disks of protostellar sources. Results: We detect 52 sources with 1 mm continuum, and we measure masses of $0.002 - 0.9$ solar masses corresponding to gas and dust in the disk and/or envelope of the protostellar system. For the deeply embedded (youngest) sources with no IR counterparts, we find evidence of mass segregation and clustering according to: the Minimum Spanning Tree method, distribution of projected separations between unique sources, and concentration of higher-mass sources near to the dense gas at the cluster center. Conclusions: The mass segregation of the mm sources is likely primordial rather than dynamical given the young age of this cluster, compared with segregation time. This is the first case to show this for mm sources in a low-mass protostellar cluster environment., Comment: Accepted to A&A, FITS data to be available at CDS

For a binary protostellar outflow system in which its members are located so close to each other (the separation being smaller than the addition of the widths of the flows) and with large opening angles, the collision seems unavoidable regardless of the orientation of the outflows. This is in contrast to the current observational evidence of just a few regions with indications of colliding outflows. Here, using sensitive observations of the Atacama Large Millimeter/Submillimeter Array (ALMA), we report resolved images of carbon monoxide (CO) towards the binary flows associated with the BHR71 protostellar system. These images reveal for the first time solid evidence that their flows are partially colliding, increasing the brightness of the CO, the dispersion of the velocities in the interaction zone, and changing part of the orientation in one of the flows. Additionally, this impact opened the possibility of knowing the 3D geometry of the system, revealing that one of its components (IRS2) should be closer to us., Comment: Submitted to AJ, comments are welcome

The mechanism for producing polarized emission from protostellar disks at (sub)millimeter wavelengths is currently uncertain. Classically, polarization is expected from non-spherical grains aligned with the magnetic field. Recently, two alternatives have been suggested. One polarization mechanism is caused by self-scattering from dust grains of sizes comparable to the wavelength while the other mechanism is due to grains aligned with their short axes along the direction of radiation anisotropy. The latter has recently been shown as a likely mechanism for causing the dust polarization detected in HL Tau at 3.1 mm. In this paper, we present ALMA polarization observations of HL Tau for two more wavelengths: 870 $\mu$m and 1.3 mm. The morphology at 870 $\mu$m matches the expectation for self-scattering, while that at 1.3 mm shows a mix between self-scattering and grains aligned with the radiation anisotropy. The observations cast doubt on the ability of (sub)millimeter continuum polarization to probe disk magnetic fields for at least HL Tau. By showing two distinct polarization morphologies at 870 $\mu$m and 3.1 mm and a transition between the two at 1.3 mm, this paper provides definitive evidence that the dominant (sub)millimeter polarization mechanism transitions with wavelength. In addition, if the polarization at 870 $\mu$m is due to scattering, the lack of polarization asymmetry along the minor axis of the inclined disk implies that the large grains responsible for the scattering have already settled into a geometrically thin layer, and the presence of asymmetry along the major axis indicates that the HL Tau disk is not completely axisymmetric., Comment: Accepted to ApJ, typos corrected from proof

We use CARMA 3mm continuum and molecular lines (NH2D, N2H+, HCO+, HCN and CS) at ~1000 au resolution to characterize the structure and kinematics of the envelope surrounding the deeply embedded first core candidate Per-Bolo 58. The line profile of the observed species shows two distinct peaks separated by 0.4-0.6 km/s, most likely arising from two different optically thin velocity components rather than the product of self-absorption in an optically thick line. The two velocity components, each with a mass of ~0.5-0.6 Msun, overlap spatially at the position of the continuum emission, and produce a general gradient along the outflow direction. We investigate whether these observations are consistent with infall in a turbulent and magnetized envelope. We compare the morphology and spectra of the N2H+(1-0) with synthetic observations of an MHD simulation that considers the collapse of an isolated core that is initially perturbed with a turbulent field. The proposed model matches the data in the production of two velocity components, traced by the isolated hyperfine line of the N2H+(1-0) spectra and shows a general agreement in morphology and velocity field. We also use large maps of the region to compare the kinematics of the core with that of the surrounding large-scale filamentary structure and find that accretion from the large-scale filament could also explain the complex kinematics exhibited by this young dense core.

We use 3mm multi-line and continuum CARMA observations towards the first hydrostatic core (FHSC) candidate L1451-mm to characterize the envelope kinematics at 1000 AU scales and investigate its evolutionary state. We detect evidence of infall and rotation in the N2H+(1-0), NH2D(1(1,1)-1(0,1)) and HCN(1-0) molecular lines. We compare the position velocity diagram of the NH2D line with a simple kinematic model and find that it is consistent with an envelope that is both infalling and rotating while conserving angular momentum around a central mass of about 0.06 Msun. The N2H+(1-0) LTE mass of the envelope along with the inferred infall velocity leads to a mass infall rate of approximately 6e-6 Msun/yr, implying a young age of 10,000 years for this FHSC candidate. Assuming that the accretion onto the central object is the same as the infall rate we obtain that the minimum source size is 1.5-5 AU consistent with the size expected for a first core. We do not see any evidence of outflow motions or signs of outflow-envelope interaction at scales > 2000 AU. This is consistent with previous observations that revealed a very compact outflow (<500 AU). We conclude that L1451-mm is indeed at a very early stage of evolution, either a first core or an extremely young Class 0 protostar. Our results provide strong evidence that L1451-mm is the best candidate for being a bonafide first core., Comment: 16 pages, 11 figures

We present a 3 mm spectral line and continuum survey of L1451 in the Perseus Molecular Cloud. These observations are from the CARMA Large Area Star Formation Survey (CLASSy), which also imaged Barnard 1, NGC 1333, Serpens Main and Serpens South. L1451 is the survey region with the lowest level of star formation activity---it contains no confirmed protostars. HCO+, HCN, and N2H+ (J=1-0) are all detected throughout the region, with HCO+ the most spatially widespread, and molecular emission seen toward 90% of the area above N(H_2) column densities of 1.9x10^21 cm^-2. HCO+ has the broadest velocity dispersion, near 0.3 km/s on average, compared to ~0.15 km/s for the other molecules, thus representing a range from supersonic to subsonic gas motions. Our non-binary dendrogram analysis reveals that the dense gas traced by each molecule has similar hierarchical structure, and that gas surrounding the candidate first hydrostatic core (FHSC), L1451-mm, and other previously detected single-dish continuum clumps have similar hierarchical structure; this suggests that different sub-regions of L1451 are fragmenting on the pathway to forming young stars. We determined the three-dimensional morphology of the largest detectable dense gas structures to be relatively ellipsoidal compared to other CLASSy regions, which appeared more flattened at largest scales. A virial analysis shows the most centrally condensed dust structures are likely unstable against collapse. Additionally, we identify a new spherical, centrally condensed N2H+ feature that could be a new FHSC candidate. The overall results suggest L1451 is a young region starting to form its generation of stars within turbulent, hierarchical structures., Comment: Accepted to The Astrophysical Journal (ApJ), 45 pages, 24 figures (some with reduced resolution in this preprint); Project website is at http://carma.astro.umd.edu/classy

The loss of mass from protostars, in the form of a jet or outflow, is a necessary counterpart to protostellar mass accretion. Outflow ejection events probably vary in their velocity and/or in the rate of mass loss. Such `episodic' ejection events have been observed during the Class 0 protostellar phase (the early accretion stage), and continue during the subsequent class I phase that marks the first one million years of star formation. Previously observed episodic-ejection sources were relatively isolated; however, the most common sites of star formation are clusters. Outflows link protostars with their environment and provide a viable source of turbulence that is necessary for regulating star formation in clusters, but it is not known how an accretion-driven jet or outflow in a clustered environment manifests itself in its earliest stage. This early stage is important in establishing the initial conditions for momentum and energy transfer to the environment as the protostar and cluster evolve. Here we report that an outflow from a very young class 0 protostar, at the hub of the very active and filamentary Serpens South protostellar cluster, shows unambiguous episodic events. The $^{12}$CO (J=2-1) emission from the protostar reveals 22 distinct features of outflow ejecta, the most recent having the highest velocity. The outflow forms bipolar lobes --- one of the first detectable signs of star formation --- which originate from the peak of 1-mm continuum emission. Emission from the surrounding C$^{18}$O envelope shows kinematics consistent with rotation and an infall of material onto the protostar. The data suggest that episodic accretion-driven outflow begins in the earliest phase of protostellar evolution, and that the outflow remains intact in a very clustered environment, probably providing efficient momentum transfer for driving turbulence., Comment: Accepted for publication as a "Letter" in Nature. 6 figures

A massive envelope and a strong bipolar outflow are the two main structures characterizing the youngest protostellar systems. In order to understand the physical properties of a bipolar outflow and the relationship with those of the envelope, we obtained a mosaic map covering the whole bipolar outflow of the youngest protostellar system L1157 with about $5"$ angular resolution in CO J=2-1 using the Combined Array for Research in Millimeter-wave Astronomy. By utilizing these observations of the whole bipolar outflow, we estimate its physical properties and show that they are consistent with multiple jets. We also constrain a preferred precession direction. In addition, we observed the central envelope structure with $2"$ resolution in the $\lambda=1.3$ and 3 mm continua and various molecular lines: C$^{17}$O, C$^{18}$O, $^{13}$CO, CS, CN, N$_2$H$^+$, CH$_3$OH, H$_2$O, SO, and SO$_2$. All the CO isotopes and CS, CN, and N$_2$H$^+$ have been detected and imaged. We marginally detected the features that can be interpreted as a rotating inner envelope in C$^{17}$O and C$^{18}$O and as an infalling outer envelope in N$_2$H$^+$. We also estimated the envelope and central protostellar masses and found that the dust opacity spectral index changes with radius., Comment: 43 pages, 13 figures, 3 tables, to be published in ApJ

We present subarcsecond (~0.35") resolved observations of the 1.3 mm dust polarization from the edge-on circumstellar disk around the Class 0 protostar L1527. The inferred magnetic field is consistent with a dominantly toroidal morphology; there is no significantly detected vertical poloidal component to which observations of an edge-on disk are most sensitive. This suggests that angular momentum transport in Class 0 protostars (when large amounts of material are fed down to the disk from the envelope and accreted onto the protostar) is driven mainly by magnetorotational instability rather than magnetocentrifugal winds at 50 AU scales. In addition, with the data to date there is an early, tentative trend that R>30 AU disks have so far been found in Class 0 systems with average magnetic fields on the 1000 AU scale strongly misaligned with the rotation axis. The absence of such a disk in the aligned case could be due to efficient magnetic braking that disrupts disk formation. If this is the case, this implies that candidate Class 0 disk systems could be identified by the average magnetic field direction at ~1000 AU spatial scales., Comment: 6 pages, 3 figures, 2 tables, accepted by ApJL