Search

While dust disks around optically visible, Class II protostars are found to be vertically thin, when and how dust settles to the midplane are unclear. As part of the Atacama Large Millimeter/submillimeter Array large program, Early Planet Formation in Embedded Disks, we analyze the edge-on, embedded, Class I protostar IRAS 04302+2247, also nicknamed the “Butterfly Star.” With a resolution of 0.″05 (8 au), the 1.3 mm continuum shows an asymmetry along the minor axis that is evidence of an optically thick and geometrically thick disk viewed nearly edge-on. There is no evidence of rings and gaps, which could be due to the lack of radial substructure or the highly inclined and optically thick view. With 0.″1 (16 au) resolution, we resolve the 2D snow surfaces, i.e., the boundary region between freeze-out and sublimation, for 12CO J = 2-1, 13CO J = 2-1, C18O J = 2-1, H 2CO J = 30,3-20,2, and SO J = 65-54, and constrain the CO midplane snow line to ∼130 au. We find Keplerian rotation around a protostar of 1.6 ± 0.4 M ⊙ using C18O. Through forward ray-tracing using RADMC-3D, we find that the dust scale height is ∼6 au at a radius of 100 au from the central star and is comparable to the gas pressure scale height. The results suggest that the dust of this Class I source has yet to vertically settle significantly.

Studying the physical and chemical conditions of young embedded disks is crucial to constrain the initial conditions for planet formation. Here we present Atacama Large Millimeter/submillimeter Array observations of dust continuum at ∼0.″06 (8 au) resolution and molecular line emission at ∼0.″17 (24 au) resolution toward the Class 0 protostar L1527 IRS from the Large Program eDisk (Early Planet Formation in Embedded Disks). The continuum emission is smooth without substructures but asymmetric along both the major and minor axes of the disk as previously observed. The detected lines of 12CO, 13CO, C18O, H2CO, c-C3H2, SO, SiO, and DCN trace different components of the protostellar system, with a disk wind potentially visible in 12CO. The 13CO brightness temperature and the H2CO line ratio confirm that the disk is too warm for CO freezeout, with the snowline located at ∼350 au in the envelope. Both molecules show potential evidence of a temperature increase around the disk-envelope interface. SO seems to originate predominantly in UV-irradiated regions such as the disk surface and the outflow cavity walls rather than at the disk-envelope interface as previously suggested. Finally, the continuum asymmetry along the minor axis is consistent with the inclination derived from the large-scale (100″ or 14,000 au) outflow, but opposite to that based on the molecular jet and envelope emission, suggesting a misalignment in the system. Overall, these results highlight the importance of observing multiple molecular species in multiple transitions to characterize the physical and chemical environment of young disks.

We present Atacama Large Millimeter Array band 6/7 (1.3 mm/0.87 mm) and Very Large Array Ka-band (9 mm) observations toward NGC 2071 IR, an intermediate-mass star-forming region. We characterize the continuum and associated molecular line emission toward the most luminous protostars, i.e., IRS1 and IRS3, on ∼100 au (0farcs2) scales. IRS1 is partly resolved in the millimeter and centimeter continuum, which shows a potential disk. IRS3 has a well-resolved disk appearance in the millimeter continuum and is further resolved into a close binary system separated by ∼40 au at 9 mm. Both sources exhibit clear velocity gradients across their disk major axes in multiple spectral lines including C18O, H2CO, SO, SO2, and complex organic molecules like CH3OH, 13CH3OH, and CH3OCHO. We use an analytic method to fit the Keplerian rotation of the disks and give constraints on physical parameters with a Markov Chain Monte Carlo routine. The IRS3 binary system is estimated to have a total mass of 1.4–1.5 M⊙. IRS1 has a central mass of 3–5 M⊙ based on both kinematic modeling and its spectral energy distribution, assuming that it is dominated by a single protostar. For both IRS1 and IRS3, the inferred ejection directions from different tracers, including radio jet, water maser, molecular outflow, and H2 emission, are not always consistent, and for IRS1 these can be misaligned by ∼50°. IRS3 is better explained by a single precessing jet. A similar mechanism may be present in IRS1 as well but an unresolved multiple system in IRS1 is also possible. © 2022. The Author(s). Published by the American Astronomical Society.

The water snowline in circumstellar disks is a crucial component in planet formation, but direct observational constraints on its location remain sparse owing to the difficulty of observing water in both young embedded and mature protoplanetary disks. Chemical imaging provides an alternative route to locate the snowline, and HCO+ isotopologues have been shown to be good tracers in protostellar envelopes and Herbig disks. Here we present similar to 0.'' 5 resolution (similar to 35 au radius) Atacama Large Millimeter/submillimeter Array (ALMA) observations of HCO+ J = 4 - 3 and (HCO+)-C-13 J = 3 - 2 toward the young (Class 0/I) disk L1527 IRS. Using a source-specific physical model with the midplane snowline at 3.4 au and a small chemical network, we are able to reproduce the HCO+ and (HCO+)-C-13 emission, but for HCO+ only when the cosmic-ray ionization rate is lowered to 10(-18) s(-1). Even though the observations are not sensitive to the expected HCO+ abundance drop across the snowline, the reduction in HCO+ above the snow surface and the global temperature structure allow us to constrain a snowline location between 1.8 and 4.1 au. Deep observations are required to eliminate the envelope contribution to the emission and to derive more stringent constraints on the snowline location. Locating the snowline in young disks directly with observations of H2O isotopologues may therefore still be an alternative option. With a direct snowline measurement, HCO+ will be able to provide constraints on the ionization rate.

We characterize protostellar multiplicity in the Orion molecular clouds using Atacama Large Millimeter/submillimeter Array 0.87 mm and Very Large Array 9 mm continuum surveys toward 328 protostars. These observations are sensitive to projected spatial separations as small as ∼20 au, and we consider source separations up to 104 au as potential companions. The overall multiplicity fraction (MF) and companion fraction (CF) for the Orion protostars are 0.30 ± 0.03 and 0.44 ± 0.03, respectively, considering separations from 20 to 104 au. The MFs and CFs are corrected for potential contamination by unassociated young stars using a probabilistic scheme based on the surface density of young stars around each protostar. The companion separation distribution as a whole is double peaked and inconsistent with the separation distribution of solar-type field stars, while the separation distribution of Flat Spectrum protostars is consistent solar-type field stars. The multiplicity statistics and companion separation distributions of the Perseus star-forming region are consistent with those of Orion. Based on the observed peaks in the Class 0 separations at ∼100 au and ∼103 au, we argue that multiples with separations <500 au are likely produced by both disk fragmentation and turbulent fragmentation with migration, and those at ≳103 au result primarily from turbulent fragmentation. We also find that MFs/CFs may rise from Class 0 to Flat Spectrum protostars between 100 and 103 au in regions of high young stellar object density. This finding may be evidence for the migration of companions from >103 au to <103 au, and that some companions between 103 and 104 au must be (or become) unbound. © 2022. The Author(s). Published by the American Astronomical Society.

We present Markov Chain Monte Carlo radiative transfer modeling of a joint ALMA 345 GHz and spectral energy distribution dataset for a sample of 97 protostellar disks from the VLA and ALMA Nascent Disk and Multiplicity Survey of Orion Protostars. From this modeling, we derive disk and envelope properties for each protostar, allowing us to examine the bulk properties of a population of young protostars. We find that disks are small, with a median dust radius of $ 29.4^{+ 4.1}_{- 2.7}$ au and a median dust mass of $ 5.8^{+ 4.6}_{- 2.7}$ M$_{\oplus}$. We find no statistically significant difference between most properties of Class 0, I, and Flat Spectrum sources with the exception of envelope dust mass and inclination. The distinction between inclination is an indication that the Class 0/I/Flat Spectrum system may be difficult to tie uniquely to the evolutionary state of protostars. When comparing with Class II disk dust masses in Taurus from similar radiative transfer modeling, we further find that the trend of disk dust mass decreasing from Class 0 to Class II disks is no longer present, though it remains unclear whether such a comparison is fair due to differences in star forming region and modeling techniques. Moreover, the disks we model are broadly gravitationally stable. Finally, we compare disk masses and radii with simulations of disk formation and find that magnetohydrodynamical effects may be important for reproducing the observed properties of disks., Comment: 31 pages, 15 figures, 2 tables, accepted for publication in ApJ

Observed changes in protostellar brightness can be complicated to interpret. In our JCMT~Transient monitoring survey, we discovered that a young binary protostar, HOPS 373, is undergoing a modest $30\%$ brightness increase at 850 $\mu$m, caused by a factor of 1.8$-$3.3 enhancement in the accretion rate. The initial burst occurred over a few months, with a sharp rise and then shallower decay. A second rise occurred soon after the decay, and the source is still bright one year later. The mid-IR emission, the small-scale CO outflow mapped with ALMA, and the location of variable maser emission indicate that the variability is associated with the SW component. The near-infrared and NEOWISE $W1$ and $W2$ emission is located along the blueshifted CO outflow, spatially offset by $\sim3$ to $4^{\prime\prime}$ from the SW component. The $K$-band emission imaged by UKIRT shows a compact H$_2$ emission source at the edge of the outflow, with a tail tracing the outflow back to the source. The $W1$ emission, likely dominated by scattered light, brightens by 0.7 mag, consistent with expectations based on the sub-mm lightcurve. The signal of continuum variability in $K$-band and $W2$ is masked by stable H$_2$ emission, as seen in our Gemini/GNIRS spectrum, and perhaps by CO emission. These differences in emission sources complicate infrared searches for variability of the youngest protostars., Comment: 29 pages, 14 figures, accepted for publication in ApJ

We present VLA observations toward the Class 0 protostar L1157 MMS at 6.8 mm and 9 mm with a resolution of ~0.04" (14 au). We detect two sources within L1157 MMS and interpret these sources as a binary protostar with a separation of ~16 au. The material directly surrounding the binary system within the inner 50 au radius of the system has an estimated mass of 0.11 M_sun, calculated from the observed dust emission. We interpret the observed binary system in the context of previous observations of its flattened envelope structure, low rates of envelope rotation from 5000 to 200 au scales, and an ordered, poloidal magnetic field aligned with the outflow. Thus, L1157 MMS is a prototype system for magnetically-regulated collapse and the presence of a compact binary within L1157 MMS demonstrates that multiple star formation can still occur within envelopes that likely have dynamically important magnetic fields., Comment: 21 pages, 6 figures, 4 tables, accepted to ApJ

At least half of a protostar's mass is accreted in the Class 0 phase, when the central protostar is deeply embedded in a dense, infalling envelope. We present the first systematic search for outbursts from Class 0 protostars in the Orion clouds. Using photometry from Spitzer/IRAC spanning 2004 to 2017, we detect three outbursts from Class 0 protostars with $\ge 2$ mag changes at 3.6 or 4.5 $\mu$m. This is comparable to the magnitude change of a known protostellar FU Ori outburst. Two are newly detected bursts from the protostars HOPS 12 and 124. The number of detections implies that Class 0 protostars burst every 438 yr, with a 95% confidence interval of 161 to 1884 yr. Combining Spitzer and WISE/NEOWISE data spanning 2004-2019, we show that the bursts persist for more than nine years with significant variability during each burst. Finally, we use $19-100$ $\mu$m photometry from SOFIA, Spitzer and Herschel to measure the amplitudes of the bursts. Based on the burst interval, a duration of 15 yr, and the range of observed amplitudes, 3-100% of the mass accretion during the Class 0 phase occurs during bursts. In total, we show that bursts from Class 0 protostars are as frequent, or even more frequent, than those from more evolved protostars. This is consistent with bursts being driven by instabilities in disks triggered by rapid mass infall. Furthermore, we find that bursts may be a significant, if not dominant, mode of mass accretion during the Class 0 phase., Comment: Accepted to ApJL

We present high-resolution Karl G. Jansky Very Large Array (VLA) observations of the protostar L1527 IRS at 7 mm, 1.3 cm, and 2 cm wavelengths. We detect the edge-on dust disk at all three wavelengths and find that it is asymmetric, with the southern side of the disk brighter than the northern side. We confirm this asymmetry through analytic modeling and also find that the disk is flared at 7 mm. We test the data against models including gap features in the intensity profile, and though we cannot rule such models out, they do not provide a statistically significant improvement in the quality of fit to the data. From these fits, we can, however, place constraints on allowed properties of any gaps that could be present in the true, underlying intensity profile. The physical nature of the asymmetry is difficult to associate with physical features owing to the edge-on nature of the disk, but it could be related to spiral arms or asymmetries seen in other imaging of more face-on disks.