Search

The composition of protoplanetary disks, and hence the initial conditions of planet formation, may be strongly influenced by the infall and thermal processing of material during the protostellar phase. Composition of dust and ice in protostellar envelopes, shaped by energetic processes driven by the protostar, serves as the fundamental building material for planets and complex organic molecules. As part of the JWST GO program, "Investigating Protostellar Accretion" (IPA), we observed an intermediate-mass protostar HOPS 370 (OMC2-FIR3) using NIRSpec/IFU and MIRI/MRS. This study presents the gas and ice phase chemical inventory revealed with the JWST in the spectral range of $\sim$2.9 to 28 $\mu$m and explores the spatial variation of volatile ice species in the protostellar envelope. We find evidence for thermal processing of ice species throughout the inner envelope. We present the first high-spatial resolution ($\sim 80$ au) maps of key volatile ice species H$_{2}$O, CO$_{2}$, $^{13}$CO$_2$, CO, and OCN$^-$, which reveal a highly structured and inhomogeneous density distribution of the protostellar envelope, with a deficiency of ice column density that coincides with the jet/outflow shocked knots. Further, we observe high relative crystallinity of H$_{2}$O ice around the shocked knot seen in the H$_2$ and OH wind/outflow, which can be explained by a lack of outer colder material in the envelope along the line of sight due to the irregular structure of the envelope. These observations show clear evidence of thermal processing of the ices in the inner envelope, close to the outflow cavity walls, heated by the luminous protostar., Comment: Submitted to ApJ. Main text:16 pages with 11 figures

We present ALMA Band 7 molecular line observations of the protostars within the VLA 1623 system. We map C$^{17}$O (3 - 2) in the circumbinary disk around VLA 1623A and the outflow cavity walls of the collimated outflow. We further detect red-shifted and blue-shifted velocity gradients in the circumstellar disks around VLA 1623B and VLA 1623W that are consistent with Keplerian rotation. We use the radiative transfer modeling code, pdspy, and simple flared disk models to measure stellar masses of $0.27 \pm 0.03$ M$_\odot$, $1.9^{+0.3}_{-0.2}$ M$_\odot$, and $0.64 \pm 0.06$ M$_\odot$ for the VLA 1623A binary, VLA 1623B, and VLA 1623W, respectively. These results represent the strongest constraints on stellar mass for both VLA 1623B and VLA 1623W, and the first measurement of mass for all stellar components using the same tracer and methodology. We use these masses to discuss the relationship between the young stellar objects (YSOs) in the VLA 1623 system. We find that VLA 1623W is unlikely to be an ejected YSO, as has been previously proposed. While we cannot rule out that VLA 1623W is a unrelated YSO, we propose that it is a true companion star to the VLA 1623A/B system and that the these stars formed in situ through turbulent fragmentation and have had only some dynamical interactions since their inception., Comment: Accepted to A&A; 16 pages, 12 figures

The envelopes and disks that surround protostars reflect the initial conditions of star and planet formation and govern the assembly of stellar masses. Characterizing these structures requires observations that span the near-infrared to centimeter wavelengths. Consequently, the past two decades have seen progress driven by numerous advances in observational facilities across this spectrum, including the \textit{Spitzer Space Telescope}, \textit{Herschel Space Observatory}, the Atacama Large Millimeter/submillimeter Array, and a host of other ground-based interferometers and single-dish radio telescopes. Nearly all protostars appear to have well-formed circumstellar disks that are likely to be rotationally-supported; the ability to detect a disk around a protostar is more a question of spatial resolution than whether or not a disk is present. The disks around protostars have inherently higher millimeter/submillimeter luminosities as compared to disks around more-evolved pre-main sequence stars, though there may be systematic variations between star forming regions. The envelopes around protostars are inherently asymmetric and streamers emphasize that mass flow through the envelopes to the disks may not be homogeneous. The current mass distribution of protostars may be impacted by selection bias given that it is skewed toward solar-mass protostars, inconsistent with the stellar initial mass function., Comment: To appear in Annual Reviews in Astronomy and Astrophysics; 42 pages, 12 Figures, 1 Table

The structure and composition of simple ices can be modified during stellar evolution by protostellar heating. Key to understanding the involved processes are thermal and chemical tracers that can diagnose the history and environment of the ice. The 15.2 $\mu$m bending mode of $^{12}$CO$_2$ has proven to be a valuable tracer of ice heating events but suffers from grain shape and size effects. A viable alternative tracer is the weaker $^{13}$CO$_2$ isotopologue band at 4.39 $\mu$m which has now become accessible at high S/N with the $\textit{James Webb}$ Space Telescope (JWST). We present JWST NIRSpec observations of $^{13}$CO$_2$ ice in five deeply embedded Class 0 sources spanning a wide range in luminosities (0.2 - 10$^4$ L$_{\odot}$ ) taken as part of the Investigating Protostellar Accretion Across the Mass Spectrum (IPA) program. The band profiles vary significantly, with the most luminous sources showing a distinct narrow peak at 4.38 $\mu$m. We first apply a phenomenological approach and show that a minimum of 3-4 Gaussian profiles are needed to fit the $^{13}$CO$_2$ absorption feature. We then combine these findings with laboratory data and show that a 15.2 $\mu$m $^{12}$CO$_2$ band inspired five-component decomposition can be applied for the isotopologue band where each component is representative of CO$_2$ ice in a specific molecular environment. The final solution consists of cold mixtures of CO$_2$ with CH$_3$OH, H$_2$O and CO as well as segregated heated pure CO$_2$ ice. Our results are in agreement with previous studies of the $^{12}$CO$_2$ ice band, further confirming that $^{13}$CO$_{2}$ is a useful alternative tracer of protostellar heating events. We also propose an alternative solution consisting only of heated CO$_2$:CH$_3$OH and CO$_2$:H$_2$O ices and warm pure CO$_2$ ice for decomposing the ice profiles of the two most luminous sources in our sample.

We present a statistical characterization of circumstellar disk orientations toward 12 protostellar multiple systems in the Perseus molecular cloud using the Atacama Large Millimeter/submillimeter Array at Band 6 (1.3 mm) with a resolution of 25 mas (8 au). This exquisite resolution enabled us to resolve the compact inner disk structures surrounding the components of each multiple system and to determine the projected 3-D orientation of the disks (position angle and inclination) to high precision. We performed a statistical analysis on the relative alignment of disk pairs to determine whether the disks are preferentially aligned or randomly distributed. We considered three subsamples of the observations selected by the companion separations, a <100 au, a >500 au, and a < 10,000 au. We found for the compact (< 100 au) subsample, the distribution of orientation angles is best described by an underlying distribution of preferentially aligned sources (within 30deg) but does not rule out distributions with 40% misaligned sources. The wide companion (>500 au) subsample appears to be consistent with a distribution of 40%-80% preferentially aligned sources. Similarly, the full sample of systems with companions (a< 10, 000 au) is most consistent with a fractional ratio of at most 80% preferentially aligned source and rules out purely randomly aligned distributions. Thus our results imply the compact sources (<100 au) and the wide companions (>500 au) are statistically different.

Investigating Protostellar Accretion (IPA) is a JWST Cycle~1 GO program that uses NIRSpec IFU and MIRI MRS to obtain 2.9--28~$\mu$m spectral cubes of young, deeply embedded protostars with luminosities of 0.2 to 10,000~L$_{\odot}$ and central masses of 0.15 to 12~M$_{\odot}$. In this Letter, we report the discovery of a highly collimated atomic jet from the Class~0 protostar IRAS~16253$-$2429, the lowest luminosity source ($L_\mathrm{bol}$ = 0.2 $L_\odot$) in the IPA program. The collimated jet is detected in multiple [Fe~II] lines, [Ne~II], [Ni~II], and H~I lines, but not in molecular emission. The atomic jet has a velocity of about 169~$\pm$~15~km\,s$^{-1}$, after correcting for inclination. The width of the jet increases with distance from the central protostar from 23 to~60 au, corresponding to an opening angle of 2.6~$\pm$~0.5\arcdeg. By comparing the measured flux ratios of various fine structure lines to those predicted by simple shock models, we derive a shock {speed} of 54~km\,s$^{-1}$ and a preshock density of 2.0$\times10^{3}$~cm$^{-3}$ at the base of the jet. {From these quantities and using a suite of jet models and extinction laws we compute a mass loss rate between $0.4 -1.1\times10^{-10}~M_{\odot}$~yr~$^{-1}$.} The low mass loss rate is consistent with simultaneous measurements of low mass accretion rate ($2.4~\pm~0.8~\times~10^{-9}~M_{\odot}$~yr$^{-1}$) for IRAS~16253$-$2429 from JWST observations (Watson et al. in prep), indicating that the protostar is in a quiescent accretion phase. Our results demonstrate that very low-mass protostars can drive highly collimated, atomic jets, even during the quiescent phase., Comment: Accepted to ApJL. Comments and feedback welcome

The Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy (SCALES) instrument is a lenslet-based integral field spectrograph that will operate at 2 to 5 microns, imaging and characterizing colder (and thus older) planets than current high-contrast instruments. Its spatial resolution for distant science targets and/or close-in disks and companions could be improved via interferometric techniques such as sparse aperture masking. We introduce a nascent Python package, NRM-artist, that we use to design several SCALES masks to be non-redundant and to have uniform coverage in Fourier space. We generate high-fidelity mock SCALES data using the scalessim package for SCALES' low spectral resolution modes across its 2 to 5 micron bandpass. We include realistic noise from astrophysical and instrument sources, including Keck adaptive optics and Poisson noise. We inject planet and disk signals into the mock datasets and subsequently recover them to test the performance of SCALES sparse aperture masking and to determine the sensitivity of various mask designs to different science signals.

SCALES (Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy) is a 2 - 5 micron high-contrast lenslet-based integral field spectrograph (IFS) designed to characterize exoplanets and their atmospheres. The SCALES medium-spectral-resolution mode uses a lenslet subarray with a 0.34 x 0.36 arcsecond field of view which allows for exoplanet characterization at increased spectral resolution. We explore the sensitivity limitations of this mode by simulating planet detections in the presence of realistic noise sources. We use the SCALES simulator scalessim to generate high-fidelity mock observations of planets that include speckle noise from their host stars, as well as other atmospheric and instrumental noise effects. We employ both angular and reference differential imaging as methods of disentangling speckle noise from the injected planet signals. These simulations allow us to assess the feasibility of speckle deconvolution for SCALES medium resolution data, and to test whether one approach outperforms another based on planet angular separations and contrasts.

The Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy (SCALES) is a $2-5~\mu$m, high-contrast integral field spectrograph (IFS) currently being built for Keck Observatory. With both low ($R\lesssim250$) and medium ($R\sim3500-7000$) spectral resolution IFS modes, SCALES will detect and characterize significantly colder exoplanets than those accessible with near-infrared ($\sim1-2~\mu$m) high-contrast spectrographs. This will lead to new progress in exoplanet atmospheric studies, including detailed characterization of benchmark systems that will advance the state of the art of atmospheric modeling. SCALES' unique modes, while designed specifically for direct exoplanet characterization, will enable a broader range of novel (exo)planetary observations as well as galactic and extragalactic studies. Here we present the science cases that drive the design of SCALES. We describe an end-to-end instrument simulator that we use to track requirements, and show simulations of expected science yields for each driving science case. We conclude with a discussion of preparations for early science when the instrument sees first light in $\sim2025$., Comment: 10 pages, 16 figures, submitted to Proceedings of the SPIE

Investigating Protostellar Accretion is a Cycle 1 JWST program using the NIRSpec+MIRI integral field units to obtain 2.9--28 $\mu$m spectral cubes of five young protostars with luminosities of 0.2-10,000 L$_{\odot}$ in their primary accretion phase. This paper introduces the NIRSpec 2.9--5.3 $\mu$m data of the inner 840-9000 au with spatial resolutions from 28-300 au. The spectra show rising continuum emission; deep ice absorption; emission from H$_{2}$, H~I, and [Fe~II]; and the CO fundamental series in emission and absorption. Maps of the continuum emission show scattered light cavities for all five protostars. In the cavities, collimated jets are detected in [Fe~II] for the four $< 320$~L$_{\odot}$ protostars, two of which are additionally traced in Br-$\alpha$. Knots of [Fe~II] emission are detected toward the most luminous protostar, and knots of [FeII] emission with dynamical times of $< 30$~yrs are found in the jets of the others. While only one jet is traced in H$_2$, knots of H$_2$ and CO are detected in the jets of four protostars. H$_2$ is seen extending through the cavities, showing that they are filled by warm molecular gas. Bright H$_2$ emission is seen along the walls of a single cavity, while in three cavities narrow shells of H$_2$ emission are found, one of which has an [Fe~II] knot at its apex. These data show cavities containing collimated jets traced in atomic/ionic gas surrounded by warm molecular gas in a wide-angle wind and/or gas accelerated by bow shocks in the jets., Comment: 26 pages, 12 figures

High-resolution, millimeter observations of disks at the protoplanetary stage reveal substructures such as gaps, rings, arcs, spirals, and cavities. While many protoplanetary disks host such substructures, only a few at the younger protostellar stage have shown similar features. We present a detailed search for early disk substructures in ALMA 1.3 and 0.87~mm observations of ten protostellar disks in the Ophiuchus star-forming region. Of this sample, four disks have identified substructure, two appear to be smooth disks, and four are considered ambiguous. The structured disks have wide Gaussian-like rings ($\sigma_R/R_{\mathrm{disk}}\sim0.26$) with low contrasts ($C<0.2$) above a smooth disk profile, in comparison to protoplanetary disks where rings tend to be narrow and have a wide variety of contrasts ($\sigma_R/R_{\mathrm{disk}}\sim0.08$ and $C$ ranges from $0-1$). The four protostellar disks with the identified substructures are among the brightest sources in the Ophiuchus sample, in agreement with trends observed for protoplanetary disks. These observations indicate that substructures in protostellar disks may be common in brighter disks. The presence of substructures at the earliest stages suggests an early start for dust grain growth and, subsequently, planet formation. The evolution of these protostellar substructures is hypothesized in two potential pathways: (1) the rings are the sites of early planet formation, and the later observed protoplanetary disk ring-gap pairs are secondary features, or (2) the rings evolve over the disk lifetime to become those observed at the protoplanetary disk stage., Comment: Accepted by ApJ, 22 pages, 10 figures

The first step toward planet formation is grain growth from (sub-)micrometer to millimeter/centimeter sizes. Grain growth has been reported not only in Class II protoplanetary disks but also in Class 0/I protostellar envelopes. However, early-stage grain growth occurring in Class 0/I stages has rarely been observed on the protostellar disk scale. Here we present the results from the ALMA Band 3 ($\lambda$ = 3.1 mm) and 7 ($\lambda$ = 0.87 mm) archival data of the Class I protostellar disk WL 17 in the $\rho$ Ophiuchus molecular cloud. Disk substructures are found in both bands, but they are different: while a central hole and a symmetric ring appear in Band 3, an off-center hole and an asymmetric ring are shown in Band 7. Furthermore, we obtain an asymmetric spectral index map with a low mean value of $\alpha$ = 2.28 $\pm$ 0.02, suggestive of grain growth and dust segregation on the protostellar disk scale. Our radiative transfer modeling verifies these two features by demonstrating that 10 cm-sized large grains are symmetrically distributed, whereas 10 $\mu$m-sized small grains are asymmetrically distributed. Also, the analysis shows that the disk is expected to be massive and gravitationally unstable. We thus suggest a single Jupiter-mass protoplanet formed by gravitational instability as the origin of the ring-like structure, grain growth, and dust segregation identified in WL 17., Comment: 23 pages, 7 figures; to be published in The Astrophysical Journal

We present observations of the Class 0 protostar IRAS 16544-1604 in CB 68 from the ''Early Planet Formation in Embedded Disks (eDisk)'' ALMA Large program. The ALMA observations target continuum and lines at 1.3-mm with an angular resolution of $\sim$5 au. The continuum image reveals a dusty protostellar disk with a radius of $\sim$30 au seen close to edge-on, and asymmetric structures both along the major and minor axes. While the asymmetry along the minor axis can be interpreted as the effect of the dust flaring, the asymmetry along the major axis comes from a real non-axisymmetric structure. The C$^{18}$O image cubes clearly show the gas in the disk that follows a Keplerian rotation pattern around a $\sim$0.14 $M_{\odot}$ central protostar. Furthermore, there are $\sim$1500 au-scale streamer-like features of gas connecting from North-East, North-North-West, and North-West to the disk, as well as the bending outflow as seen in the $^{12}$CO (2-1) emission. At the apparent landing point of NE streamer, there are SO (6$_5$-5$_4$) and SiO (5-4) emission detected. The spatial and velocity structure of NE streamer can be interpreted as a free-falling gas with a conserved specific angular momentum, and the detection of the SO and SiO emission at the tip of the streamer implies presence of accretion shocks. Our eDisk observations have unveiled that the Class 0 protostar in CB 68 has a Keplerian rotating disk with flaring and non-axisymmetric structure associated with accretion streamers and outflows., Comment: 30 pages, 24 figures, accepted for publication in The Astrophysical Journal as one of the first-look papers of the eDisk ALMA Large Program

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 (ALMA) 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 which 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 $^{12}$CO $J$=2--1, $^{13}$CO $J$=2--1, C$^{18}$O $J$=2--1, $H_{2}$CO $J$=$3_{0,3}$--$2_{0,2}$, and SO $J$=$6_{5}$--$5_{4}$, and constrain the CO midplane snow line to $\sim 130$ au. We find Keplerian rotation around a protostar of $1.6 \pm 0.4 M_{\odot}$ using C$^{18}$O. Through forward ray-tracing using RADMC-3D, we find that the dust scale height is $\sim 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., Comment: 33 pages, 21 figures. Accepted for publication in ApJ as one of the first-look papers of the eDisk ALMA Large Program

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 (ALMA) observations of dust continuum at $\sim$0.06" (8 au) resolution and molecular line emission at $\sim$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 $^{12}$CO, $^{13}$CO, C$^{18}$O, H$_2$CO, c-C$_3$H$_2$, SO, SiO, and DCN trace different components of the protostellar system, with a disk wind potentially visible in $^{12}$CO. The $^{13}$CO brightness temperature and the H$_2$CO line ratio confirm that the disk is too warm for CO freeze out, with the snowline located at $\sim$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., Comment: 27 pages, 16 figures, 2 tables, 10 pages appendix with 12 figures. Accepted for publication in ApJ as one of the first-look papers of the eDisk ALMA Large Program

We present an overview of the Large Program, ``Early Planet Formation in Embedded Disks (eDisk)'', conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The ubiquitous detections of substructures, particularly rings and gaps, in protoplanetary disks around T Tauri stars raise the possibility that at least some planet formation may have already started during the embedded stages of star formation. In order to address exactly how and when planet formation is initiated, the program focuses on searching for substructures in disks around 12 Class 0 and 7 Class I protostars in nearby ($< $200 pc) star-forming regions through 1.3 mm continuum observations at a resolution of $\sim7$ au (0.04"). The initial results show that the continuum emission, mostly arising from dust disks around the sample protostars, has relatively few distinctive substructures, such as rings and spirals, in marked contrast to Class II disks. The dramatic difference may suggest that substructures quickly develop in disks when the systems evolve from protostars to Class II sources or alternatively that high optical depth of the continuum emission could obscure internal structures. Kinematic information obtained through CO isotopologue lines and other lines reveals the presence of Keplerian disks around protostars, providing us with crucial physical parameters, in particular, the dynamical mass of the central protostars. We describe the background of the eDisk program, the sample selection and their ALMA observations, the data reduction, and also highlight representative first-look results., Comment: This is a publication of a series of eDisk ALMA large program first-look papers

We present Atacama Large Millimeter Array (ALMA) observations of the young, gas-rich debris disk around HD110058 at 0.3-0.6\arcsec resolution. The disk is detected in the 0.85 and 1.3~mm continuum, as well as in the J=2-1 and J=3-2 transitions of $^{12}$CO and $^{13}$CO. The observations resolve the dust and gas distributions and reveal that this is the smallest debris disk around stars of similar luminosity observed by ALMA. The new ALMA data confirm the disk is very close to edge-on, as shown previously in scattered light images. We use radiative transfer modeling to constrain the physical properties of dust and gas disks. The dust density peaks at around 31~au and has a smooth outer edge that extends out to $\sim70$~au. Interestingly, the dust emission is marginally resolved along the minor axis, which indicates that it is vertically thick if truly close to edge-on with an aspect ratio between 0.13 and 0.28. We also find that the CO gas distribution is more compact than the dust \ah{(similarly to the disk around 49 Ceti)}, which could be due to a low viscosity and a higher gas release rate at small radii. Using simulations of the gas evolution taking into account the CO photodissociation, shielding, and viscous evolution, we find that HD~110058's CO gas mass and distribution are consistent with a secondary origin scenario. Finally, we find that the gas densities may be high enough to cause the outward drift of small dust grains in the disk.

We present an 870 $\mu$m continuum survey of 300 protostars from the Herschel Orion Protostar Survey using the Atacama Compact Array (ACA). These data measure protostellar flux densities on envelope scales $\leq$8000 au (20") and resolve the structure of envelopes with 1600 au (4") resolution, a factor of 3-5 improvement in angular resolution over existing single-dish 870 $\mu$m observations. We compare the ACA observations to Atacama Large Millimeter/submillimeter Array 12 m array observations at 870 $\mu$m with $\sim$0.1 (40 au) resolution. Using the 12 m data to measure the fluxes from disks and the ACA data within 2500 au to measure the combined disk plus envelope fluxes, we calculate the 12 m/ACA 870 $\mu$m flux ratios. Our sample shows a clear evolution in this ratio. Class 0 protostars are mostly envelope-dominated with ratios $<$0.5. In contrast, Flat Spectrum protostars are primarily disk-dominated with ratios near 1, although with a number of face-on protostars dominated by their envelopes. Class I protostars span the range from envelope to disk-dominated. The increase in ratio is accompanied by a decrease in the envelope fluxes and estimated mass infall rates. We estimate that 80$\%$ of the mass is accreted during the envelope-dominated phase. We find that the 12 m/ACA flux ratio is an evolutionary indicator that largely avoids the inclination and foreground extinction dependence of spectral energy distribution-based indicators.

VLA 1623 West is an ambiguous source that has been described as a shocked cloudlet as well as a protostellar disk. We use deep ALMA 1.3 and 0.87 millimeter observations to constrain its shape and structure to determine its origins better. We use a series of geometric models to fit the uv visibilities at both wavelengths with GALARIO. Although the Real visibilities show structures similar to what has been identified as gaps and rings in protoplanetary disks, we find that a modified Flat-Topped Gaussian at high inclination provides the best fit to the observations. This fit agrees well with expectations for an optically thick, highly inclined disk. Nevertheless, we find that the geometric models consistently yield positive residuals at the four corners of the disk at both wavelengths. We interpret these residuals as evidence that the disk is flared in the millimeter dust. We use a simple toy model for an edge-on flared disk and find that the residuals best match a disk with flaring that is mainly restricted to the outer disk at $R \gtrsim 30$ au. Thus, VLA 1623W may represent a young protostellar disk where the large dust grains have not yet had enough time to settle into the mid-plane. This result may have implications for how disk evolution and vertical dust settling impact the initial conditions leading to planet formation., Comment: Accepted by ApJ, 12 pages, 13 figures

We present the design of SCALES (Slicer Combined with Array of Lenslets for Exoplanet Spectroscopy) a new 2-5 micron coronagraphic integral field spectrograph under construction for Keck Observatory. SCALES enables low-resolution (R~50) spectroscopy, as well as medium-resolution (R~4,000) spectroscopy with the goal of discovering and characterizing cold exoplanets that are brightest in the thermal infrared. Additionally, SCALES has a 12x12" field-of-view imager that will be used for general adaptive optics science at Keck. We present SCALES's specifications, its science case, its overall design, and simulations of its expected performance. Additionally, we present progress on procuring, fabricating and testing long lead-time components., Comment: 19 pages, 12 Figures, SPIE Astronomical Telescopes + Instrumentation, (#12184-18), 2022, Montreal, Canada

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 due to the edge-on nature of the disk, but could be related to spiral arms or asymmetries seen in other imaging of more face-on disks., Comment: 20 pages, 5 figures, 3 tables; accepted for publication in the Astrophysical Journal

We present ALMA band 6/7 (1.3 mm/0.87 mm) and VLA Ka band (9 mm) observations toward NGC 2071 IR, an intermediate-mass star forming region. We characterize the continuum and associated molecular line emission towards the most luminous protostars, i.e., IRS1 and IRS3, on ~100 au (0. 2") scales. IRS1 is partly resolved in millimeter and centimeter continuum, which shows a potential disk. IRS3 has a well resolved disk appearance in 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 MCMC routine. The IRS3 binary system is estimated to have a total mass of 1.4-1.5$M_\odot$. IRS1 has a central mass of 3-5$M_\odot$ 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$^{\circ}$. 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., Comment: 36 pages, 21 figures, accepted by ApJ

We report an Atacama Large Millimeter/submillimeter Array 0.88 mm (Band 7) continuum detection of the accretion disk around SR 12 c, an $\sim$11 $M_{\rm Jup}$ planetary-mass companion (PMC) orbiting its host binary at 980 au. This is the first submillimeter detection of a circumplanetary disk around a wide PMC. The disk has a flux density of $127 \pm14~\mu$Jy and is not resolved by the $\sim$0.1" beam, so the dust disk radius is likely less than 5 au and can be much smaller if the dust continuum is optically thick. If, however, the dust emission is optically thin, then the SR 12 c disk has a comparable dust mass to the circumplanetary disk around PDS 70 c but is about five times lower than that of the $\sim$12 $M_{\rm Jup}$ free-floating OTS 44. This suggests that disks around bound and unbound planetary-mass objects can span a wide range of masses. The gas mass estimated with an accretion rate of $10^{-11}~M_\odot$ yr$^{-1}$ implies a gas-to-dust ratio higher than 100. If cloud absorption is not significant, a nondetection of ${}^{12}$CO(3-2) implies a compact gas disk around SR 12 c. Future sensitive observations may detect more PMC disks at 0.88 mm flux densities of $\lesssim$100 $\mu$Jy., Comment: Published in ApJL

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

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 characterize protostellar multiplicity in the Orion molecular clouds using ALMA 0.87~mm and VLA 9~mm continuum surveys toward 328 protostars. These observations are sensitive to projected spatial separations as small as $\sim$20~au, and we consider source separations up to 10$^4$~au as potential companions. The overall multiplicity fraction (MF) and companion fraction (CF) for the Orion protostars are 0.30$\pm$0.03 and 0.44$\pm$0.03, respectively, considering separations from 20 to 10$^4$~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 $\sim$100~au and $\sim$10$^3$~au, we argue that multiples with separations $<$500~au are likely produced by both disk fragmentation and turbulent fragmentation with migration, and those at $\ga$10$^3$~au result primarily from turbulent fragmentation. We also find that MFs/CFs may rise from Class 0 to Flat Spectrum protostars between 100 and 10$^3$~au in regions of high YSO density. This finding may be evidence for migration of companions from $>$10$^3$~au to $<$10$^3$~au, and that some companions between 10$^3$ and 10$^4$~au must be (or become) unbound., Comment: 76 pages, 20 Figures, 10 Tables, accepted by the Astrophysical Journal

WL 17 is a young transition disk in the Ophiuchus L1688 molecular cloud complex. Even though WL 17 is among the brightest disks in L1688 and massive enough to expect dust self-scattering, it was undetected in polarization down to ALMA's instrument sensitivity limit. Such low polarization fractions could indicate unresolved polarization within the beam or optically thin dust emission. We test the latter case by combining the high sensitivity 233 GHz Stokes I data from the polarization observations with previous ALMA data at 345 GHz and 100 GHz. We use simple geometric shapes to fit the observed disk visibilities in each band. Using our simple models and assumed dust temperature profiles, we estimate the optical depth in all three bands. The optical depth at 233 GHz peaks at $\tau_{233} \sim 0.3$, which suggests the dust emission may not be optically thick enough for dust self-scattering to be efficient. We also find the higher sensitivity 233 GHz data show substructure in the disk for the first time. The substructure appears as brighter lobes along the major axis, on either side of the star. We attempt to fit the lobes with a simple geometric model, but they are unresolved in the 233 GHz data. We propose that the disk may be flared at 1 mm such that there is a higher column of dust along the major axis than the minor axis when viewed at an inclination. These observations highlight the strength of high sensitivity continuum data from dust polarization observations to study disk structures., Comment: Accepted to ApJ, 13 pages, 6 figures

Stellar masses are fundamental but often difficult to measure. Thanks to the Atacama Large Millimeter Array (ALMA) and $Gaia$, dynamical masses of pre-main sequence stars can be precisely measured using the Keplerian rotation of protoplanetary disks. We used ALMA CO(2-1) and CO(3-2) observations of CT Cha and DS Tau to determine their masses by modeling the geometry, kinematics, and physical properties of their disks with a Bayesian-based radiative transfer modeling code $pdspy$. We found that the posterior distributions of the masses from the two transitions are inconsistent at the 2-4 $\sigma$ level. These systematic errors may originate from assumptions in the disk model, or perhaps the modest spatial or spectral resolutions used in this study. Regardless, this indicates that dynamical mass measurements using disk kinematics should be treated with caution when using only a single transition line because of these systematic errors., Comment: Published on RNAAS

We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations towards a compact (230~au separation) triple protostar system, L1448 IRS3B, at 879~\micron with \contbeam~resolution. Spiral arm structure within the circum-multiple disk is well resolved in dust continuum toward IRS3B, and we detect the known wide (2300~au) companion, IRS3A, also resolving possible spiral substructure. Using dense gas tracers, C17O, H13CO$+$, and H13CN, we resolve the Keplerian rotation for both the circum-triple disk in IRS3B and the disk around IRS3A. Furthermore, we use the molecular line kinematic data and radiative transfer modeling of the molecular line emission to confirm that the disks are in Keplerian rotation with fitted masses of $1.19^{+0.13}_{-0.07}$ for IRS3B-ab, $1.51^{+0.06}_{-0.07}$~Msun for IRS3A, and place an upper limit on the central protostar mass for the tertiary IRS3B-c of 0.2~Msun. We measure the mass of the fragmenting disk of IRS3B to be 0.29~Msun from the dust continuum emission of the circum-multiple disk and estimate the mass of the clump surrounding IRS3B-c to be 0.07~Msun. We also find that the disk around IRS3A has a mass of 0.04~Msun. By analyzing the Toomre~Q parameter, we find the IRS3A circumstellar disk is gravitationally stable (Q$>$5), while the IRS3B disk is consistent with a gravitationally unstable disk (Q$<$1) between the radii 200-500~au. This coincides with the location of the spiral arms and the tertiary companion IRS3B-c, supporting the hypothesis that IRS3B-c was formed in situ via fragmentation of a gravitationally unstable disk.

We present ALMA (0.87~mm and 1.3~mm) and VLA (9~mm) observations toward the candidate intermediate-mass protostar OMC2-FIR3 (HOPS-370; L$_{bol}$~314~L$_{\odot}$) at $\sim$0.1" (40~au) resolution for the continuum emission and ~0.25" (100 au) resolution of nine molecular lines. The dust continuum observed with ALMA at 0.87~mm and 1.3~mm resolve a near edge-on disk toward HOPS-370 with an apparent radius of ~100 au. The VLA observations detect both the disk in dust continuum and free-free emission extended along the jet direction. The ALMA observations of molecular lines (H$_2$CO, SO, CH$_3$OH, $^{13}$CO, C$^{18}$O, NS, and H$^{13}$CN) reveal rotation of the apparent disk surrounding HOPS-370 orthogonal to the jet/outflow direction. We fit radiative transfer models to both the dust continuum structure of the disk and molecular line kinematics of the inner envelope and disk for the H$_2$CO, CH$_3$OH, NS, and SO lines. The central protostar mass is determined to be $\sim$2.5 M_sun with a disk radius of $\sim$94~au, when fit using combinations of the H$_2$CO, CH$_3$OH, NS, and SO lines, consistent with an intermediate-mass protostar. Modeling of the dust continuum and spectral energy distribution (SED) yields a disk mass of 0.035~M$_{\odot}$ (inferred dust+gas) and a dust disk radius of 62~au, thus the dust disk may have a smaller radius than the gas disk, similar to Class II disks. In order to explain the observed luminosity with the measured protostar mass, HOPS-370 must be accreting at a rate between 1.7 and 3.2$\times$10$^{-5}$~M$_{\odot}$~yr$^{-1}$., Comment: Accepted to ApJ; 51 pages, 12 Figures, 7 Tables

We present observations toward HOPS 383, the first known outbursting Class 0 protostar located within the Orion molecular cloud using ALMA, VLA, and SMA. The SMA observations reveal envelope scale continuum and molecular line emission surrounding HOPS 383 at 0.85 mm, 1.1 mm, and 1.3 mm. The images show that HCO$^+$ and H$^{13}$CO$^+$ peaks on or near the continuum, while N$_2$H$^+$ is reduced at the same position. This reflects the underlying chemistry where CO evaporating close to the protostar destroys N$_2$H$^+$ while forming HCO$^+$. We also observe the molecular outflow traced by $^{12}$CO ($J = 2 \rightarrow 1$) and ($J = 3 \rightarrow 2$). A disk is resolved in the ALMA 0.87 mm dust continuum, orthogonal to the outflow direction, with an apparent radius of $\sim$62 AU. Radiative transfer modeling of the continuum gives disk masses of 0.02 M$_{\odot}$ when fit to the ALMA visibilities. The models including VLA 8 mm data indicate that the disk mass could be up to a factor of 10 larger due to lower dust opacity at longer wavelengths. The disk temperature and surface density profiles from the modeling, and an assumed protostar mass of 0.5 M$_{\odot}$ suggest that the Toomre $Q$ parameter $< 1$ before the outburst, making gravitational instability a viable mechanism to explain outbursts at an early age if the disk is sufficiently massive., Comment: Accepted by ApJ

The prevalence of substructures in $\sim1-10$ Myr old protoplanetary disks, which are often linked to planet formation, has raised the question of how early such features form, and as a corollary, how early planet formation begins. Here we present observations of seven protostellar disks (aged $\sim0.1-1$ Myr) from the VLA/ALMA Nascent Disk and Multiplicity Survey of Orion Protostars (VANDAM: Orion) that show clear substructures, thereby demonstrating that these features can form early in the lifetimes of disks. We use simple analytic models as well as detailed radiative transfer modeling to characterize their structure. In particular we show that at least four of the sources have relatively massive envelopes, indicating that they are particularly young, likely the youngest disks with substructures known to-date. Several of these disks also have emission from an inner disk that is offset from the center of the ring structure. Given the size of the cleared out regions of the disk, it is unclear, however, whether these features are related to planet formation, or rather if they are signposts of close-separation binary formation at early times., Comment: 21 pages, 11 figures, 4 tables, accepted for publication in ApJ

We present Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm observations of four young, eruptive star-disk systems at 0.4" resolution: two FUors (V582 Aur and V900 Mon), one EXor (UZ Tau E) and one source with an ambiguous FU/EXor classification (GM Cha). The disks around GM Cha, V900 Mon and UZ Tau E are resolved. These observations increase the sample of FU/EXors observed at sub-arcsecond resolution by 15%. The disk sizes and masses of FU/EXors objects observed by ALMA so far suggest that FUor disks are more massive than Class 0/I disks in Orion and Class II disks in Lupus of similar size. EXor disks in contrast do not seem to be distinguishable from these two populations. We reach similar conclusions when comparing the FU/EXor sample to the Class I and Class II disks in Ophiuchus. FUor disks around binaries are host to more compact disks than those in single-star systems, similar to non-eruptive young disks. We detect a wide-angle outflow around GM Cha in $^{12}$CO emission, wider than typical Class I objects and more similar to those found around some FUor objects. We use radiative transfer models to fit the continuum and line data of the well-studied disk around UZ Tau E. The line data is well described by a keplerian disk, with no evidence of outflow activity (similar to other EXors). The detection of wide-angle outflows in FUors and not in EXors support to the current picture in which FUors are more likely to represent an accretion burst in the protostellar phase (Class I), while EXors are smaller accretion events in the protoplanetary (Class II) phase., Comment: 22 pages, 14 figures. Accepted for publication in ApJ

Characterizing the physical properties and compositions of circumplanetary disks can provide important insights into the formation of giant planets and satellites. We report ALMA 0.88 mm (Band 7) continuum observations of six planetary-mass (10-20 $M_{\rm Jup}$) companions: CT Cha b, 1RXS 1609 b, ROXs 12 b, ROXs 42B b, DH Tau b, and FU Tau b. No continuum sources are detected at the locations of the companions down to 3$\sigma$ limits of 120-210 $\mu$Jy. Given these non-detections, it is not clear whether disks around planetary-mass companions indeed follow the disk flux-host mass trend in the stellar regime. The faint radio brightness of these companion disks may result from a combination of fast radial drift and a lack of dust traps. Alternatively, as disks in binary systems are known to have significantly lower millimeter fluxes due to tidal interactions, these companion disks may instead follow the relationship of moderate-separation binary stars. This scenario can be tested with sensitive continuum imaging at rms levels of $\lesssim$10 $\mu$Jy., Comment: Accepted for publication in AJ

We report ALMA and VLA continuum observations that potentially identify the four youngest protostars in the Orion Molecular Clouds taken as part of the Orion VANDAM program. These are distinguished by bright, extended, irregular emission at 0.87 mm and 8 mm and are optically thick at 0.87 mm. These structures are distinct from the disk or point-like morphologies seen toward the other Orion protostars. The 0.87 mm emission implies temperatures of 41-170 K, requiring internal heating. The bright 8 mm emission implies masses of 0.5 to 1.2 M_sun assuming standard dust opacity models. One source has a Class 0 companion, while another exhibits substructure indicating a companion-candidate. Three compact outflows are detected, two of which may be driven by companions, with dynamical times of ~300 to ~400 years. The slowest outflow may be driven by a first hydrostatic core. These protostars appear to trace an early phase when the centers of collapsing fragments become optically thick to their own radiation and compression raises the gas temperature. This phase is thought to accompany the formation of hydrostatic cores. A key question is whether these structures are evolving on free fall times of ~100 years, or whether they are evolving on Kelvin-Helmholtz times of several thousand years. The number of these sources imply a lifetime of ~6000 years, in closer agreement with the Kelvin-Helmholtz time. In this case, rotational and/or magnetic support could be slowing the collapse., Comment: accepted by ApJ, 22 pages, 11 figures, 5 tables

We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase., Comment: 107 pages, 21 Figures, 11 Tables, accepted to ApJ. Version with all source figures: https://www.cv.nrao.edu/~jtobin/Orion-disks-fullfigs.pdf Reduced data available from https://dataverse.harvard.edu/dataverse/VANDAMOrion

We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC2-FIR4, four in OMC2-FIR3, and one additional source just outside OMC2-FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100~AU scales. HOPS-108 is the only protostar in OMC2-FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC2-FIR3 with L~$\sim$~360~\lsun, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate methanol from icy dust grains; overall these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in \twco, but find evidence of interaction between the outflow/jet from HOPS-370 and the OMC2-FIR4 region. The multitude of observational constraints indicate that HOPS-108 is likely a low to intermediate-mass protostar in its main mass accretion phase and it is the most luminous protostar in OMC2-FIR4. The high resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them., Comment: 32 Pages, 10 Figures, 5 Tables, accepted to ApJ

We present new 0.6 - 4 $\mu$m imaging of the SR 21 transition disk from Keck/NIRC2 and Magellan/MagAO. The protoplanetary disk around SR 21 has a large (~ 30 - 40 AU) clearing first inferred from its spectral energy distribution and later detected in sub-millimeter imaging. Both the gas and small dust grains are known to have a different morphology, with an inner truncation in CO at ~7 AU, and micron-sized dust detected within the millimeter clearing. Previous near-infrared imaging could not distinguish between an inner dust disk with a truncation at ~7 AU or one that extended to the sublimation radius. The imaging data presented here require an inner dust disk radius of a few AU, and complex structure such as a warp or spiral. We present a parametric warped disk model that can reproduce the observations. Reconciling the images with the spectral energy distribution gathered from the literature suggests grain growth to ~2 - 5 $\mu$m within the sub-millimeter clearing. The complex disk structure and possible grain growth can be connected to dynamical shaping by a giant-planet mass companion, a scenario supported by previous observational and theoretical studies., Comment: 22 pages, 14 figures, accepted for publication in ApJ

Search engines are some of the most popular and profitable intelligent technologies in existence. Recent research, however, has suggested that search engines may be surprisingly dependent on user-created content like Wikipedia articles to address user information needs. In this paper, we perform a rigorous audit of the extent to which Google leverages Wikipedia and other user-generated content to respond to queries. Analyzing results for six types of important queries (e.g. most popular, trending, expensive advertising), we observe that Wikipedia appears in over 80% of results pages for some query types and is by far the most prevalent individual content source across all query types. More generally, our results provide empirical information to inform a nascent but rapidly-growing debate surrounding a highly-consequential question: Do users provide enough value to intelligent technologies that they should receive more of the economic benefits from intelligent technologies?, Comment: This version includes a bibliography entry that was missing from the first version of the text due to a processing error. This is a preprint of a paper accepted at ICWSM 2019. Please cite that version instead

Knowledge of protostellar evolution has been revolutionized with the advent of surveys at near-infrared to submillimeter wavelengths. This has enabled the bolometric luminosities and bolometric temperatures (traditional protostellar evolution diagnostics) to be measured for large numbers of protostars. However, further progress is difficult without knowing the masses of the central protostars. Protostar masses can be most accurately determined via molecular line kinematics from millimeter interferometers (i.e., ALMA). Theoretical investigations have predicted the protostellar mass function (PMF) for various protostellar mass accretion models, and it is now imperative to observationally constrain its functional form. While ALMA has enabled protostellar mass measurements, samples approaching 100 sources are necessary to constrain the functional form of the PMF, and upgrades to ALMA and/or a new mm/cm facility will increase the feasibility of measuring such a large number of protostar masses. The masses of protostars will enable their stellar structure (radius and intrinsic luminosity), evolution, and accretion histories to be better understood. This is made more robust when effective temperatures and accretion rates can be measured via ground/space-based near to mid-infrared spectroscopy. Furthermore, access to supercomputing facilities is essential to fit the protostar masses via radiative transfer modeling and updated theoretical/numerical modeling of stellar structure may also be required., Comment: Science white paper submitted to the Astro2020 US Decadal Survey. 7 pages and 3 figures

Significant advances have been made over the past decade in the characterization of multiple protostar systems, enabled by the Karl G. Jansky Very Large Array (VLA), high-resolution infrared observations with the Hubble Space Telescope, and ground-based facilities. To further understand the mechanism(s) of multiple star formation, a combination of statistics, high-angular resolution radio/millimeter continuum imaging, characterization of kinematic structure, magnetic fields via polarimetry, and comparison with numerical simulations are needed. Thus, understanding the origin of stellar multiplicity in different regimes of companion separation will soon be within reach. However, to overcome challenges that studies in this field are now confronted with, a range of new capabilities are required: a new millimeter/centimeter wave facility with 10 mas resolution at {\lambda}=1 cm, space-based near to far-infrared observatories, continued development of low to high resolution spectroscopy on 3m to 10m class telescopes, and an ELT-class telescope with near to mid-infrared imaging/spectroscopic capability., Comment: Science white paper submitted to the Astro2020 US Decadal Survey. 7 pages and 2 figures

The Keplerian rotation in protoplanetary disks can be used to robustly measure stellar masses at very high precision if the source distance is known. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of spatially and spectrally resolved $^{12}$CO (2--1) emission towards the disks around 2MASS J16262774--2527247 (the tertiary companion to ROXs 12 at 5100 au), CT Cha, and DH Tau. We employ detailed modeling of the Keplerian rotation profile, coupled with accurate distances from $Gaia$, to directly measure the stellar masses with $\sim$2\% precision. We also compare these direct mass measurements with the masses inferred from evolutionary models, determined in a statistically rigorous way. We find that 2MASS J16262774--2527247 has a mass of $0.535^{+0.006}_{-0.007}~M_{\odot}$ and CT Cha has a mass of $0.796^{+0.015}_{-0.014}~M_{\odot}$, broadly consistent with evolutionary models, although potentially significant differences remain. DH Tau has a mass of $0.101^{+0.004}_{-0.003}~M_{\odot}$, but it suffers from strong foreground absorption that may affect our mass estimate. The combination of ALMA, $Gaia$, and codes like \texttt{pdspy}, presented here, can be used to infer the dynamical masses for large samples of young stars and substellar objects, and place constraints on evolutionary models., Comment: 22 pages, 16 figures, 3 tables, accepted for publication in ApJ

The formation and evolution of disks early in the protostellar phase is an area of study in which the ngVLA is poised to make significant breakthroughs. The high-sensitivity and resolution at wavelengths of 3~mm and longer will enable forming disks to be examined with unprecedented detail. The need to observe dust emission at wavelengths of 3~mm and longer is motivated by the fact that dust emission at these wavelengths is more likely to be optically thin, which is essential to understanding the structure of these disks. We explore the feasibility of detecting and resolving protostellar disks with a variety of radii, masses, and distances, out to distances as large as 1.5~kpc using radiative transfer models and simulations with the proposed ngVLA configuration. We also examine the potential for the ngVLA to enable studies of grain growth and radial migration of dust particles early in the protostellar phase with the broad multi-wavelength coverage. Studies of grain growth will require wavelength coverage extending at least to $\sim$4~cm to characterize and quantify the location and intensity of free-free emission, which is expected to be generated at $<$10 AU scales., Comment: To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA); 11 pages, 4 figures

The ngVLA will enable significant advances in our understanding of the formation and evolution of multiple star systems in the protostellar phase, building upon the breakthroughs enabled by the VLA. The high-sensitivity and resolution at 3~mm wavelengths and longer will enable closer multiple systems to be discovered in the nearby star forming regions. The ngVLA is incredibly important for multiplicity studies because dust opacity at short wavelengths ($<$3~mm) can hide multiplicity and the long wavelengths are needed to reveal forming multiples in the youngest systems. The samples sizes can be expanded to encompass star forming regions at distances of at least 1.5~kpc, enabling statistical studies that are on par with studies of field star multiplicity. We verify the capability of the ngVLA to detect and resolve multiple star systems at distances out to 1.5~kpc using empirical examples of systems detected by the VLA and scaling them to greater distances. We also use radiative transfer models and simulations to verify that the ngVLA can resolve close binary systems from their dust emission at these distances. The ngVLA will also have excellent imaging capability and the circum-multiple environments can also be examined in great detail., Comment: To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA); 11 pages, 5 figures