In recent years evidence has been building that planet formation starts early, in the first $\sim$ 0.5 Myr. Studying the dust masses available in young disks enables understanding the origin of planetary systems since mature disks are lacking the solid material necessary to reproduce the observed exoplanetary systems, especially the massive ones. We aim to determine if disks in the embedded stage of star formation contain enough dust to explain the solid content of the most massive exoplanets. We use Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 observations of embedded disks in the Perseus star-forming region together with Very Large Array (VLA) Ka-band (9 mm) data to provide a robust estimate of dust disk masses from the flux densities. Using the DIANA opacity model including large grains, with a dust opacity value of $\kappa_{\rm 9\ mm}$ = 0.28 cm$^{2}$ g$^{-1}$, the median dust masses of the embedded disks in Perseus are 158 M$_\oplus$ for Class 0 and 52 M$_\oplus$ for Class I from the VLA fluxes. The lower limits on the median masses from ALMA fluxes are 47 M$_\oplus$ and 12 M$_\oplus$ for Class 0 and Class I, respectively, obtained using the maximum dust opacity value $\kappa_{\rm 1.3mm}$ = 2.3 cm$^{2}$ g$^{-1}$. The dust masses of young Class 0 and I disks are larger by at least a factor of 10 and 3, respectively, compared with dust masses inferred for Class II disks in Lupus and other regions. The dust masses of Class 0 and I disks in Perseus derived from the VLA data are high enough to produce the observed exoplanet systems with efficiencies acceptable by planet formation models: the solid content in observed giant exoplanets can be explained if planet formation starts in Class 0 phase with an efficiency of $\sim$ 15%. Higher efficiency of $\sim$ 30% is necessary if the planet formation is set to start in Class I disks.

We explore the relation between the stellar mass surface density and the mass surface density of molecular hydrogen gas in twelve nearby molecular clouds that are located at $<$1.5 kpc distance. The sample clouds span an order of magnitude range in mass, size, and star formation rates. We use thermal dust emission from $Herschel$ maps to probe the gas surface density and the young stellar objects from the most recent $Spitzer$ Extended Solar Neighborhood Archive (SESNA) catalog to probe the stellar surface density. Using a star-sampled nearest neighbor technique to probe the star-gas surface density correlations at the scale of a few parsecs, we find that the stellar mass surface density varies as a power-law of the gas mass surface density, with a power-law index of $\sim$2 in all the clouds. The consistent power-law index implies that star formation efficiency is directly correlated with gas column density, and no gas column density threshold for star formation is observed. We compare the observed correlations with the predictions from an analytical model of thermal fragmentation, and with the synthetic observations of a recent hydrodynamic simulation of a turbulent star-forming molecular cloud. We find that the observed correlations are consistent for some clouds with the thermal fragmentation model and can be reproduced using the hydrodynamic simulations.

Recent observations have suggested that circumstellar disks may commonly form around young stellar objects. Although the formation of circumstellar disks can be a natural result of the conservation of angular momentum in the parent cloud, theoretical studies instead show disk formation to be difficult from dense molecular cores magnetized to a realistic level, owing to efficient magnetic braking that transports a large fraction of the angular momentum away from the circumstellar region. We review recent progress in the formation and early evolution of disks around young stellar objects of both low-mass and high-mass, with an emphasis on mechanisms that may bridge the gap between observation and theory, including non-ideal MHD effects and asymmetric perturbations in the collapsing core (e.g., magnetic field misalignment and turbulence). We also address the associated processes of outflow launching and the formation of multiple systems, and discuss possible implications in properties of protoplanetary disks.

The collapse of the protostellar envelope results in the growth of the protostar and the development of a protoplanetary disk, playing a critical role during the early stages of star formation. Characterizing the gas infall in the envelope constrains the dynamical models of star formation. We present unambiguous signatures of infall, probed by optically thick molecular lines, toward an isolated embedded protostar, BHR 71 IRS1. The three dimensional radiative transfer calculations indicate that a slowly rotating infalling envelope model following the "inside-out" collapse reproduces the observations of both HCO$^{+}$ $J=4\rightarrow3$ and CS $J=7\rightarrow6$ lines, and the low velocity emission of the HCN $J=4\rightarrow3$ line. The envelope has a model-derived age of 12000$\pm$3000 years after the initial collapse. The envelope model underestimates the high velocity emission at the HCN $J=4\rightarrow3$ and H$^{13}$CN $J=4\rightarrow3$ lines, where outflows or a Keplerian disk may contribute. The ALMA observations serendipitously discover the emission of complex organic molecules (COMs) concentrated within a radius of 100 au, indicating that BHR 71 IRS1 harbors a hot corino. Eight species of COMs are identified, including CH$_{3}$OH and CH$_{3}$OCHO, along with H$_{2}$CS, SO$_{2}$ and HCN $v_{2}=1$. The emission of methyl formate and $^{13}$C-methanol shows a clear velocity gradient within a radius of 50 au, hinting at an unresolved Keplerian rotating disk.

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.

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.

We present and release the full dataset for the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. This survey used the Submillimeter Array (SMA) to image the 74 known protostars within the Perseus molecular cloud. The SMA was used in two array configurations to capture outflows for scales $>$30$^{\prime\prime}$ ($>$9000 au) and to probe scales down to $\sim$1$^{\prime\prime}$ ($\sim$300 au). The protostars were observed with the 1.3 mm and 850 $\mu$m receivers simultaneously to detect continuum at both wavelengths and molecular line emission from CO(2-1), $^{13}$CO(2-1), C$^{18}$O(2-1), N$_2$D$^+$(3-2), CO(3-2), HCO$^+$(4-3), and H$^{13}$CO$^+$(4-3). Some of the observations also used the SMA's recently upgraded correlator, SWARM, whose broader bandwidth allowed for several more spectral lines to be observed (e.g., SO, H$_2$CO, DCO$^+$, DCN, CS, CN). Of the main continuum and spectral tracers observed, 84% of the images and cubes had emission detected. The median C$^{18}$O(2-1) linewidth is $\sim$1.0 km s$^{-1}$, which is slightly higher than those measured with single-dish telescopes at scales of 3000-20000 au. Of the 74 targets, six are suggested to be first hydrostatic core candidates, and we suggest that L1451-mm is the best candidate. We question a previous continuum detection toward L1448 IRS2E. In the SVS13 system, SVS13A certainly appears to be the most evolved source, while SVS13C appears to be hotter and more evolved than SVS13B. The MASSES survey is the largest publicly available interferometric continuum and spectral line protostellar survey to date, and is largely unbiased as it only targets protostars in Perseus. All visibility ($uv$) data and imaged data are publicly available at https://dataverse.harvard.edu/dataverse/full_MASSES/.

The fastest molecular component to the protostellar outflows -- extremely high-velocity (EHV) molecular jets -- are still puzzling since they are seen only rarely. The first aim is to analyze the interaction between the EHV jet and the slow outflow by comparing their outflow force content. The second aim is to analyze the chemical composition of the different outflow velocity components and to reveal the spatial location of molecules. ALMA 3 mm and 1.3 mm observations of five outflow sources at 130 -- 260 au resolution in the Serpens Main cloud are presented. Observations of CO, SiO, H$_2$CO and HCN reveal the kinematic and chemical structure of those flows. Three velocity components are distinguished: the slow and the fast wing, and the EHV jet. Out of five sources, three have the EHV component. Comparison of outflow forces reveals that only the EHV jet in the youngest source Ser-emb 8 (N) has enough momentum to power the slow outflow. The SiO abundance is generally enhanced with velocity, while HCN is present in the slow and the fast wing, but disappears in the EHV jet. For Ser-emb 8 (N), HCN and SiO show a bow-shock shaped structure surrounding one of the EHV peaks suggesting sideways ejection creating secondary shocks upon interaction with the surroundings. Also, the SiO abundance in the EHV gas decreases with distance from this protostar, whereas that in the fast wing increases. H$_2$CO is mostly associated with low-velocity gas but also appears surprisingly in one of the bullets in the Ser-emb~8~(N) EHV jet. The high detection rate suggests that the presence of the EHV jet may be more common than previously expected. The origin and temporal evolution of the abundances of SiO, HCN and H$_2$CO through high-temperature chemistry are discussed. The data are consistent with a low C/O ratio in the EHV gas versus high C/O ratio in the fast and slow wings.

We present 1.3 mm ALMA observations of polarized dust emission toward the wide-binary protostellar system BHR 71 IRS1 and IRS2. IRS1 features what appears to be a natal, hourglass-shaped magnetic field. In contrast, IRS2 exhibits a magnetic field that has been affected by its bipolar outflow. Toward IRS2, the polarization is confined mainly to the outflow cavity walls. Along the northern edge of the redshifted outflow cavity of IRS2, the polarized emission is sandwiched between the outflow and a filament of cold, dense gas traced by N$_2$D$^+$, toward which no dust polarization is detected. This suggests that the origin of the enhanced polarization in IRS2 is the irradiation of the outflow cavity walls, which enables the alignment of dust grains with respect to the magnetic field -- but only to a depth of ~300 au, beyond which the dust is cold and unpolarized. However, in order to align grains deep enough in the cavity walls, and to produce the high polarization fraction seen in IRS2, the aligning photons are likely to be in the mid- to far-infrared range, which suggests a degree of grain growth beyond what is typically expected in very young, Class 0 sources. Finally, toward IRS1 we see a narrow, linear feature with a high (10-20%) polarization fraction and a well ordered magnetic field that is not associated with the bipolar outflow cavity. We speculate that this feature may be a magnetized accretion streamer; however, this has yet to be confirmed by kinematic observations of dense-gas tracers.

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.

We present the first results of high-spectral resolution (0.023 km/s) N$_2$H$^+$ observations of dense gas dynamics at core scales (~0.01 pc) using the recently commissioned Argus instrument on the Green Bank Telescope (GBT). While the fitted linear velocity gradients across the cores measured in our targets nicely agree with the well-known power-law correlation between the specific angular momentum and core size, it is unclear if the observed gradients represent core-scale rotation. In addition, our Argus data reveal detailed and intriguing gas structures in position-velocity (PV) space for all 5 targets studied in this project, which could suggest that the velocity gradients previously observed in many dense cores actually originate from large-scale turbulence or convergent flow compression instead of rigid-body rotation. We also note that there are targets in this study with their star-forming disks nearly perpendicular to the local velocity gradients, which, assuming the velocity gradient represents the direction of rotation, is opposite to what is described by the classical theory of star formation. This provides important insight on the transport of angular momentum within star-forming cores, which is a critical topic on studying protostellar disk formation.

Observations of young stellar objects (YSOs) in centimeter bands can probe the continuum emission from growing dust grains, ionized winds, and magnetospheric activity, which are intimately connected to the evolution of protoplanetary disks and the formation of planets. We have carried out sensitive continuum observations toward the Ophiuchus A star-forming region using the Karl G. Jansky Very Large Array (VLA) at 10 GHz over a field-of-view of 6$'$ with a spatial resolution of $\theta_{maj}$ $\times$ $\theta_{min}$ $\sim$ 0.4$''$ $\times$ 0.2$''$. We achieved a 5 $\mu$Jy beam$^{-1}$ root-mean-square noise level at the center of our mosaic field of view. Among the eighteen sources we detected, sixteen are YSOs (three Class 0, five Class I, six Class II, and two Class III) and two are extragalactic candidates. We find that thermal dust emission generally contributes less that 30% of the emission at 10 GHz. The radio emission is dominated by other types of emission such as gyro-synchrotron radiation from active magnetospheres, free-free emission from thermal jets, free-free emission from the outflowing photo-evaporated disk material, and/or synchrotron emission from accelerated cosmic-rays in jet or protostellar surface shocks. These different types of emission could not be clearly disentangled. Our non-detections towards Class II/III disks suggest that extreme UV-driven photoevaporation is insufficient to explain the disk dispersal, assuming that the contribution of UV photoevaporating stellar winds to radio flux does not evolve with time. The sensitivity of our data cannot exclude photoevaporation due to X-ray photons as an efficient mechanism for disk dispersal. Deeper surveys with the Square Kilometre Array will be able to provide strong constraints on disk photoevaporation.

We present 0.25 arcsec resolution (35 au) ALMA 1.3 mm dust polarization observations for 37 young stellar objects (YSOs) in the Ophiuchus cloud. These data encompass all the embedded protostars in the cloud and several Flat and Class II objects to produce the largest, homogeneous study of dust polarization on disk scales to date. The goal of this study is to study dust polarization down to disk scales. We find that 14/37 (38%) of the YSOs are detected in polarization. Nine of these sources have uniform polarization angles and four sources have azimuthal polarization. The sources with uniform polarization tend to have steeper inclinations (> 60 degree) than those with azimuthal polarization (< 60 degree). The majority (9/14) of the detected sources have polarization morphologies and disk properties consistent with dust self-scattering in optically thick disks. The remaining sources may be instead tracing magnetic fields. Their inferred field directions from rotating the polarization vectors by 90 degree are mainly poloidal or hourglass shaped. We find no evidence of a strong toroidal field component toward any of our disks. For the 23 YSOs that are undetected in polarization, roughly half of them have 3-sigma upper limits of < 2%. These sources also tend to have inclinations < 60 degree and they are generally compact. Since lower inclination sources tend to have azimuthal polarization, these YSOs may be undetected in polarization due to unresolved polarization structure within our beam. We propose that disks with inclinations > 60 degree are the best candidates for future polarization studies of dust self-scattering as these systems will generally show uniform polarization vectors that do not require very high resolution to resolve. We release the continuum and polarization images for all the sources with this publication. Data from the entire survey can be obtained from Dataverse.

The determination of the specific angular momentum radial profile, $j(r)$, in the early stages of star formation is crucial to constrain star and circumstellar disk formation theories. The specific angular momentum is directly related to the largest Keplerian disk possible, and it could constrain the angular momentum removal mechanism. We determine $j(r)$ towards two Class 0 objects and a first hydrostatic core candidate in the Perseus cloud, which is consistent across all three sources and well fit with a single power-law relation between 800 and 10,000\,au: $j_{fit}(r)=10^{-3.60\pm0.15}\left(r/\textrm{1,000au}\right)^{1.80\pm0.04}$ km s$^{-1}$ pc. This power-law relation is in between solid body rotation ($\propto r^2$) and pure turbulence ($\propto r^{1.5}$). This strongly suggests that even at 1,000\,au, the influence of the dense core's initial level of turbulence or the connection between core and the molecular cloud is still present. The specific angular momentum at 10,000\,au is $\approx3\times$ higher than previously estimated, while at 1,000\,au it is lower by $2\times$. We do not find a region of conserved specific angular momentum, although it could still be present at a smaller radius. We estimate an upper limit to the largest Keplerian disk radius of 60\,au, which is small but consistent with published upper limits. Finally, these results suggest that more realistic initial conditions for numerical simulations of disk formation are needed. Some possible solutions include: a) use a larger simulation box to include some level of driven turbulence or connection to the parental cloud, or b) incorporate the observed $j(r)$ to setup the dense core kinematics initial conditions.

We have carried out polarimetric observations with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the Class 0 protostellar system L1448 IRS 2, which is a proto-binary embedded within a flattened, rotating structure, and for which a hint of a central disk has been suggested, but whose magnetic fields are aligned with the bipolar outflow on the cloud core scale. Our high sensitivity and high resolution ($\sim 100$ au) observations show a clear hourglass magnetic field morphology centered on the protostellar system, but the central pattern is consistent with a toroidal field indicative of a circumstellar disk, although other interpretations are also possible, including field lines dragged by an equatorial accretion flow into a configuration parallel to the midplane. If a relatively large disk does exist, it would suggest that the magnetic braking catastrophe is averted in this system, not through a large misalignment between the magnetic and rotation axes, but rather through some other mechanisms, such as non-ideal magneto-hydrodynamic effects and/or turbulence. We have also found a relationship of decreasing polarization fractions with intensities and the various slopes of this relationship can be understood as multiple polarization mechanisms and/or depolarization from a changing field morphology. In addition, we found a prominent clumpy depolarization strip crossing the center perpendicular to the bipolar outflow. Moreover, a rough estimate of the magnetic field strength indicates that the field is strong enough to hinder formation of a rotationally supported disk, which is inconsistent with the feature of a central toroidal field.

Context. The bipolar outflows associated with low-mass protostars create cavities in the infalling envelope. These cavities are illuminated by the central protostar and inner disk, creating a bipolar scattered light nebula at near-infrared and mid-infrared wavelengths. The variability of the scattered light nebula in both total intensity and intensity as a function of position in the scattered light nebula can provide important insights into the structure of the inner disk that cannot be spatially resolved. Aims. We aim to determine the likelihood that a warped inner disk is the origin of the surface brightness variability in the bipolar scattered light nebula associated with L1527 IRS. Methods. We present results from near-IR imaging conducted over the course of seven years, with periods of monthly cadence monitoring. We used Monte Carlo radiative transfer models to interpret the observations. Results. We find a time varying, asymmetrical brightness in the scattered light nebulae within the outflow cavities of the protostar. Starting in 2007, the surface brightnesses of the eastern and western outflow cavities were roughly symmetric. Then, in 2009, the surface brightnesses of the cavities were found to be asymmetric, with a substantial increase in surface brightness and a larger increase in the eastern outflow cavity. More regular monitoring was conducted from 2011 to 2014, revealing a rotating pattern of surface brightness variability in addition to a slow change of the eastern and western outflow cavities toward symmetry, but still not as symmetric as observed in 2007. We find that an inner disk warp is a feasible mechanism to produce the rotating pattern of surface brightness variability.

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.

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.

Background: low-mass stars are the dominant product of the star formation process, and they trace star formation over the full range of environments, from isolated globules to clusters in the central molecular zone. In the past two decades, our understanding of the spatial distribution and properties of young low-mass stars and protostars has been revolutionized by sensitive space-based observations at X-ray and IR wavelengths. By surveying spatial scales from clusters to molecular clouds, these data provide robust measurements of key star formation properties. Goal: with their large numbers and their presence in diverse environments, censuses of low mass stars and protostars can be used to measure the dependence of star formation on environmental properties, such as the density and temperature of the natal gas, strengths of the magnetic and radiation fields, and the density of stars. Here we summarize how such censuses can answer three basic questions: i.) how is the star formation rate influenced by environment, ii.) does the IMF vary with environment, and iii.) how does the environment shape the formation of bound clusters? Answering these questions is an important step toward understanding star and cluster formation across the extreme range of environments found in the Universe. Requirements: sensitivity and angular resolution improvements will allow us to study the full range of environments found in the Milky Way. High spatial dynamic range (< 1arcsec to > 1degree scales) imaging with space-based telescopes at X-ray, mid-IR, and far-IR and ground-based facilities at near-IR and sub-mm wavelengths are needed to identify and characterize young stars.

The majority of the ultimate main-sequence mass of a star is assembled in the protostellar phase, where a forming star is embedded in an infalling envelope and encircled by a protoplanetary disk. Studying mass accretion in protostars is thus a key to understanding how stars gain their mass and ultimately how their disks and planets form and evolve. At this early stage, the dense envelope reprocesses most of the luminosity generated by accretion to far-infrared and submillimeter wavelengths. Time-domain photometry at these wavelengths is needed to probe the physics of accretion onto protostars, but variability studies have so far been limited, in large part because of the difficulty in accessing these wavelengths from the ground. We discuss the scientific progress that would be enabled with far-infrared and submillimeter programs to probe protostellar variability in the nearest kiloparsec.

Understanding the collapse of clouds and the formation of protoplanetary disks is essential to understanding the formation of stars and planets. Infall and accretion, the mass-aggregation processes that occur at envelope and disk scales, drive the dynamical evolution of protostars. While the observations of protostars at different stages constrain their evolutionary tracks, the impact of variability due to accretion bursts on dynamical and chemical evolution of the source is largely unknown. The lasting effects on protostellar envelopes and disks are tracked through multi-wavelength and time domain observational campaigns, requiring deep X-ray, infrared, and radio imaging and spectroscopy, at a sufficient level of spatial detail to distinguish contributions from the various substructures (i.e., envelope from disk from star from outflow). Protostellar models derived from these campaigns will illuminate the initial chemical state of protoplanetary disks during the epoch of giant planet formation. Insight from individual star formation in the Milky Way is also necessary to understand star formation rates in extragalactic sources. This cannot be achieved with ground-based observatories and is not covered by currently approved instrumentation. Requirements: High (v < 10 km/s for survey; v < 1 km/s for followup) spectral resolution capabilities with relatively rapid response times in the IR (3-500 um), X-ray (0.1-10 keV), and radio (cm) are critical to follow the course of accretion and outflow during an outburst. Complementary, AU-scale radio observations are needed to probe the disk accretion zone, and 10 AU-scale to probe chemical and kinematic structures of the disk-forming regions, and track changes in the dust, ice, and gas within protostellar envelopes.

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.

In the standard picture for low-mass star formation, a dense molecular cloud undergoes gravitational collapse to form a protostellar system consisting of a new central star, a circumstellar disk, and a surrounding envelope of remaining material. The mass distribution of the system evolves as matter accretes from the large-scale envelope through the disk and onto the protostar. While this general picture is supported by simulations and indirect observational measurements, the specific timescales related to disk growth and envelope dissipation remain poorly constrained. In this paper we conduct a rigorous test of a method introduced by J{\o}rgensen et al. (2009) to obtain mass measurements of disks and envelopes around embedded protostars with observations that do not resolve the disk (resolution of $\sim$1000\,AU). Using unresolved data from the recent Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, we derive disk and envelope mass estimates for $59$ protostellar systems in the Perseus molecular cloud. We compare our results to independent disk mass measurements from the VLA Nascent Disk and Multiplicity (VANDAM) survey and find a strong linear correlation, suggesting that accurate disk masses can be measured from unresolved observations. Then, leveraging the size of the MASSES sample, we find no significant trend in protostellar mass distribution as a function of age, as approximated from bolometric temperatures. These results may indicate that the disk mass of a protostar is set near the onset of the Class 0 protostellar stage and remains roughly constant throughout the Class I protostellar stage.

Since the discovery of the multi-ring structure of the HL Tau disk, ALMA data suggest that the dust continuum emission of many, if not all, protoplanetary disks consists of rings and gaps, no matter their spectral type or age. The origin of these gaps so far remains unclear. We present a sample study of 16 disks with multiple ring-like structures in the continuum, using published ALMA archival data, to compare their morphologies and gap locations in a systematic way. The 16 targets range from early to late type stars, from <0.5 Myr to >10 Myr, from ~0.2 to 40 L_Sun and include both full and transitional disks with cleared inner dust cavities. Stellar ages are revised using new Gaia distances. Gap locations are derived using a simple radial fit to the intensity profiles. Using a radiative transfer model, the temperature profiles are computed. The gap radii generally do not correspond to the orbital radii of snow lines of the most common molecules. A snow line model can likely be discarded as a common origin of multi-ring systems. In addition, there are no systematic trends in the gap locations that could be related to resonances of planets. Finally, the outer radius of the disks decreases for the oldest disks in the sample, indicating that if multi-ring disks evolve in a similar way, outer dust rings either dissipate with the gas or grow into planetesimal belts.

We present high-resolution (~ 35 au) ALMA Band 6 1.3 mm dust polarization observations of IRAS 16293. These observations spatially resolve the dust polarization across the two protostellar sources and toward the filamentary structures between them. The dust polarization and inferred magnetic field have complicated structures throughout the region. In particular, we find that the magnetic field is aligned parallel to three filamentary structures. We characterize the physical properties of the filamentary structure that bridges IRAS 16293A and IRAS 16293B and estimate a magnetic field strength of 23-78 mG using the Davis-Chandrasekhar-Fermi method. We construct a toy model for the bridge material assuming that the young stars dominate the mass and gravitational potential of the system. We find that the expected gas flow to each star is of comparable order to the Alfven speed, which suggests that the field may be regulating the gas flow. We also find that the bridging material should be depleted in ~ 1000 yr. If the bridge is part of the natal filament that formed the stars, then it must have accreted new material. Alternatively, the bridge could be a transient structure. Finally, we show that the 1.3 mm polarization morphology of the optically thick IRAS 16293B system is qualitatively similar to dust self-scattering. Based on similar polarization measurements at 6.9 mm, we propose that IRAS 16293B has produced a substantial population of large dust grains with sizes between 200 and 2000 um.

We present Gaia-DR2 astrometry of a sample of YSO candidates in Ophiuchus, Serpens Main and Serpens South/W40 in the Aquila Rift, which had been mainly identified by their infrared excess with Spitzer. We compare the Gaia-DR2 parallaxes against published and new parallaxes obtained from our Very Long Baseline Array (VLBA) program GOBELINS. We obtain consistent results between Gaia and the VLBA for the mean parallaxes in each of the regions analyzed here. We see small offsets, when comparing mean values, of a few tens of micro-arcseconds in the parallaxes, which are either introduced by the Gaia zero-point error or due to a selection effect by Gaia toward the brightest, less obscured stars. Gaia-DR2 data alone conclusively places Serpens Main and Serpens South at the same distance, as we first inferred from VLBA data alone in a previous publication. Thus, Serpens Main, Serpens South and W40 are all part of the same complex of molecular clouds, located at a mean distance of 436+/-9 pc. In Ophiuchus, both Gaia and VLBA suggest a small parallax gradient across the cloud, and the distance changes from 144.2+/-1.3 pc to 138.4+/-2.6 pc when going from L1689 to L1688.

We present a characterization of the binary protostar system that is forming within a dense core in the isolated dark cloud BHR71. The pair of protostars, IRS1 and IRS2, are both in the Class 0 phase, determined from observations that resolve the sources from 1 um out to 250 um and from 1.3 mm to 1.3cm. The resolved observations enable the luminosities of IRS1 and IRS2 to be independently measured (14.7 and 1.7L_sun, respectively), in addition to the bolometric temperatures 68~K, and 38~K, respectively. The surrounding core was mapped in NH3 (1,1) with the Parkes radio telescope, and followed with higher-resolution observations from ATCA in NH3 (1,1) and 1.3cm continuum. The protostars were then further characterized with ALMA observations in the 1.3~mm continuum along with N2D+ (J=3-2), 12CO, 13CO, and C18O (J=2-1) molecular lines. The Parkes observations find evidence for a velocity gradient across the core surrounding the two protostars, while ATCA reveals more complex velocity structure toward the protostars within the large-scale gradient. The ALMA observations then reveal that the two protostars are at the same velocity in C18O, and N2H+ exhibits a similar velocity structure as NH3. However, the C18O kinematics reveal that the rotation on scales $<$1000~AU around IRS1 and IRS2 are in opposite directions. Taken with the lack of a systematic velocity difference between the pair, it is unlikely that their formation resulted from rotational fragmentation. We instead conclude that the binary system most likely formed via turbulent fragmentation of the core.

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.

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.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of multiple protostar systems in the Perseus molecular cloud previously detected by the Karl G. Jansky Very Large Array (VLA). We observed 17 close ($<$600~AU separation) multiple systems at 1.3~mm in continuum and five molecular lines (i.e., \twco, \cateo, \thco, H$_2$CO, SO) to characterize the circum-multiple environments in which these systems are forming. We detect at least one component in the continuum for the 17 multiple systems. In three systems, one companion is not detected, and for two systems the companions are unresolved at our observed resolution. We also detect circum-multiple dust emission toward 8 out of 9 Class 0 multiples. Circum-multiple dust emission is not detected toward any of the 8 Class I multiples. Twelve systems are detected in the dense gas tracers toward their disks/inner envelopes. For these 12 systems, we use the dense gas observations to characterize their formation mechanism. The velocity gradients in the circum-multiple gas are clearly orthogonal to the outflow directions in 8 out of the 12 systems, consistent with disk fragmentation. Moreover, only two systems with separations $<$200~AU are \textit{inconsistent} with disk fragmentation, in addition to the two widest systems ($>$500~AU). Our results suggest that disk fragmentation via gravitational instability is an important formation mechanism for close multiple systems, but further statistics are needed to better determine the relative fraction formed via this method.

Context: Simulations suggest that gas heating due to radiative feedback is a key factor in whether or not multiple protostellar systems will form. Chemistry is a good tracer of the physical structure of a protostellar system, since it depends on the temperature structure. Aims: We aim to study the relationship between envelope gas temperature and protostellar multiplicity. Methods: Single dish observations of various molecules that trace the cold, warm and UV-irradiated gas are used to probe the temperature structure of multiple and single protostellar systems on 7000 AU scales. Results: Single, close binary, and wide multiples present similar current envelope gas temperatures, as estimated from H$_2$CO and DCO$^+$ line ratios. The temperature of the outflow cavity, traced by c-C$_3$H$_2$, on the other hand, shows a relation with bolometric luminosity and an anti-correlation with envelope mass. Although the envelope gas temperatures are similar for all objects surveyed, wide multiples tend to exhibit a more massive reservoir of cold gas compared to close binary and single protostars. Conclusions: Although the sample of protostellar systems is small, the results suggest that gas temperature may not have a strong impact on fragmentation. We propose that mass, and density, may instead be key factors in fragmentation.

We present the full disk-fit results VANDAM survey of all Class 0 and I protostars in the Perseus molecular cloud. We have 18 new protostellar disk candidates around Class 0 and I sources, which are well described by a simple, parametrized disk model fit to the 8 mm VLA dust-continuum observations. 33% of Class 0 protostars and just 11% of Class I protostars have candidate disks, while 78% of Class 0 and I protostars do not have signs of disks within our 12 AU disk diameter resolution limit, indicating that at 8 mm most disks in the Class 0 and I phases are <10 AU in radius. These small radii may be a result of surface brightness sensitivity limits. Modeled 8 mm radii are similar to the radii of known Class 0 disks with detected Keplerian rotation. Since our 8 mm data trace a population of larger dust grains which radially drift towards the protostar and are lower limits on true disk sizes, large disks at early times do not seem to be particularly rare. We find statistical evidence that Class 0 and I disks are likely drawn from the same distribution, meaning disk properties may be defined early in the Class 0 phase and do not undergo large changes through the Class I phase. By combining our candidate disk properties with previous polarization observations, we find a qualitative indication that misalignment between inferred envelope-scale magnetic fields and outflows may indicate disks on smaller scales in Class 0 sources.

We report the detection of methanol in the disk around the young outbursting star V883 Ori with the Atacama Large Millimeter/submillimeter Array (ALMA). Four transitions are observed with upper level energies ranging between 115 and 459 K. The emission is spatially resolved with the 0.14" beam and follows the Keplerian rotation previously observed for C$^{18}$O. Using a rotational diagram analysis, we find a disk-averaged column density of $\sim10^{17}$ cm$^{-2}$ and a rotational temperature of $\sim90-100$ K, suggesting that the methanol has thermally desorbed from the dust grains. We derive outer radii between 120 and 140 AU for the different transitions, compared to the 360 AU outer radius for C$^{18}$O. Depending on the exact physical structure of the disk, the methanol emission could originate in the surface layers beyond the water snowline. Alternatively, the bulk of the methanol emission originates inside the water snowline, which can then be as far out as ~100 AU, instead of 42 AU as was previously inferred from the continuum opacity. In addition, these results show that outbursting young stars like V883 Ori are good sources to study the ice composition of planet forming material through thermally desorbed complex molecules, which have proven to be hard to observe in more evolved protoplanetary disks.

We derive the distance and structure of the Perseus molecular cloud by combining trigonometric parallaxes from Very Long Baseline Array (VLBA) observations, taken as part of the GOBELINS survey, and Gaia Data Release 2. Based on our VLBA astrometry, we obtain a distance of 321+/-10 pc for IC 348. This is fully consistent with the mean distance of 320+/-26 measured by Gaia. The VLBA observations toward NGC 1333 are insufficient to claim a successful distance measurement to this cluster. Gaia parallaxes, on the other hand, yield a mean distance of 293+/-22 pc. Hence, the distance along the line of sight between the eastern and western edges of the cloud is ~30 pc, which is significantly smaller than previously inferred. We use Gaia proper motions and published radial velocities to derive the spatial velocities of a selected sample of stars. The average velocity vectors with respect to the LSR are (u,v,w) = (-6.1+/-1.6, 6.8+/-1.1, -0.9+/-1.2) and (-6.4+/-1.0, 2.1+/-1.4, -2.4+/-1.0) km/s for IC 348 and NGC 1333, respectively. Finally, our analysis of the kinematics of the stars has shown that there is no clear evidence of expansion, contraction, or rotational motions within the clusters.

We present the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey, which uses the Submillimeter Array (SMA) interferometer to map the continuum and molecular lines for all 74 known Class 0/I protostellar systems in the Perseus molecular cloud. The primary goal of the survey is to observe an unbiased sample of young protostars in a single molecular cloud so that we can characterize the evolution of protostars. This paper releases the MASSES 1.3 mm data from the subcompact configuration ($\sim$4$^{\prime\prime}$ or $\sim$1000 au resolution), which is the SMA's most compact array configuration. We release both $uv$ visibility data and imaged data for the spectral lines CO(2-1),$^{13}$CO(2-1), C$^{18}$O(2-1), and N$_2$D$^+$(3-2), as well as for the 1.3 mm continuum. We identify the tracers that are detected toward each source. We also show example images of continuum and CO(2-1) outflows, analyze C$^{18}$O(2-1) spectra, and present data from the SVS 13 star-forming region. The calculated envelope masses from the continuum show a decreasing trend with bolometric temperature (a proxy for age). Typical C$^{18}$O(2-1) linewidths are 1.45 km s$^{-1}$, which is higher than the C$^{18}$O linewidths detected toward Perseus filaments and cores. We find that N$_2$D$^+$(3-2) is significantly more likely to be detected toward younger protostars. We show that the protostars in SVS 13 are contained within filamentary structures as traced by C$^{18}$O(2-1) and N$_2$D$^+$(3-2). We also present the locations of SVS 13A's high velocity (absolute line-of-sight velocities $>$150 km s$^{-1}$) red and blue outflow components. Data can be downloaded from https://dataverse.harvard.edu/dataverse/MASSES .

Emission from protostars at centimeter radio wavelengths has been shown to trace the free-free emission arising from ionizing shocks as a result of jets and outflows driven by protostars. Therefore, measuring properties of protostars at radio frequencies can provide valuable insights into the nature of their outflows and jets. We present a C-band (4.1 cm and 6.4 cm) survey of all known protostars (Class 0 and Class I) in Perseus as part of the VLA Nascent Disk and Multiplicity (VANDAM) Survey. We examine the known correlations between radio flux density and protostellar parameters such as bolometric luminosity and outflow force, for our sample. We also investigate the relationship between radio flux density and far-infrared line luminosities from Herschel. We show that free-free emission originates most likely from J-type shocks; however, the large scatter indicates that those two types of emission probe different time and spatial scales. Using C-band fluxes, we removed an estimation of free-free contamination from the corresponding Ka-band (9 mm) flux densities that primarily probe dust emission from embedded disks. We find that the compact ($<$~1") dust emission is lower for Class I sources (median dust mass 96 M$_{\oplus}$) relative to Class 0 (248 M$_{\oplus}$), but several times higher than in Class II (5-15 M$_{\oplus}$). If this compact dust emission is tracing primarily the embedded disk, as is likely for many sources, this result provides evidence for decreasing disk masses with protostellar evolution, with sufficient mass for forming giant planet cores primarily at early times.

We present high resolution (~ 30 au) ALMA Band 6 dust polarization observations of VLA 1623. The VLA 1623 data resolve compact ~ 40 au inner disks around the two protobinary sources, VLA 1623-A and VLA 1623-B, and also an extended ~ 180 au ring of dust around VLA 1623-A. This dust ring was previously identified as a large disk in lower-resolution observations. We detect highly structured dust polarization toward the inner disks and the extended ring with typical polarization fractions ~ 1.7% and ~ 2.4%, respectively. The two components also show distinct polarization morphologies. The inner disks have uniform polarization angles aligned with their minor axes. This morphology is consistent with expectations from dust scattering. By contrast, the extended dust ring has an azimuthal polarization morphology not previously seen in lower-resolution observations. We find that our observations are well fit by a static, oblate spheroid model with a flux-frozen, poloidal magnetic field. We propose that the polarization traces magnetic grain alignment likely from flux freezing on large scales and magnetic diffusion on small scales. Alternatively, the azimuthal polarization may be attributed to grain alignment by the anisotropic radiation field. If the grains are radiatively aligned, then our observations indicate that large (~ 100 um) dust grains grow quickly at large angular extents. Finally, we identify significant proper motion of VLA 1623 using our observations and those in the literature. This result indicates that the proper motion of nearby systems must be corrected for when combining ALMA data from different epochs.

We present new trigonometric parallaxes and proper motions of young stellar objects in the Taurus molecular cloud complex from observations collected with the Very Long Baseline Array as part of the Gould's Belt Distances Survey (GOBELINS). We detected 26 young stellar objects and derived trigonometric parallaxes for 18 stars with an accuracy of 0.3$\%$ to a few percent. We modeled the orbits of six binaries and determined the dynamical masses of the individual components in four of these systems (V1023 Tau, T Tau S, V807 Tau and V1000 Tau). Our results are consistent with the first trigonometric parallaxes delivered by the Gaia satellite and reveal the existence of significant depth effects. We find that the central portion of the dark cloud Lynds 1495 is located at $d=129.5\pm 0.3$ pc while the B 216 clump in the filamentary structure connected to it is at $d=158.1\pm 1.2$ pc. The closest and remotest stars in our sample are located at $d=126.6\pm 1.7$ pc and $d=162.7\pm 0.8$ pc yielding a distance difference of about 36 pc. We also provide a new distance estimate for HL Tau that was recently imaged. Finally, we compute the spatial velocity of the stars with published radial velocity and investigate the kinematic properties of the various clouds and gas structures in this region.

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.

[Abridged] Protoplanetary disks have been studied extensively, both physically and chemically, to understand the environment in which planets form. However, the first steps of planet formation are likely to occur already when the protostar and disk are still embedded in their natal envelope. The initial conditions for planet formation may thus be provided by these young embedded disks, of which the physical and chemical structure is poorly characterized. We aim to constrain the midplane temperature structure, one of the critical unknowns, of the embedded disk around L1527. In particular, we set out to determine whether there is an extended cold outer region where CO is frozen out, as is the case for Class II disks. We use archival ALMA data to directly observe the midplane of the near edge-on L1527 disk. Optically thick $^{13}$CO ($J=2-1$) and C$^{18}$O ($J=2-1$) emission is observed throughout the disk and inner envelope, while N$_2$D$^+ (J=3-2$), which can only be abundant when CO is frozen out, is not detected. Both CO isotopologues have brightness temperatures $\gtrsim$ 25 K along the midplane. Disk and envelope emission can be disentangled kinematically, because the largest velocities are reached in the disk. A power law radial temperature profile constructed using the highest midplane temperature at these velocities suggest that the temperature is above 20 K out to at least 75 AU, and possibly throughout the entire 125 AU disk. Radiative transfer models show that a model without CO freeze-out in the disk matches the C$^{18}$O observations better than a model with the CO snowline at $\sim$70 AU. In addition, there is no evidence for a large (order of magnitude) depletion of CO. The disk around L1527 is likely to be warm enough to have CO present in the gas phase throughout the disk, suggesting that young embedded disks can indeed be warmer than the more evolved Class II disks.

We present a characterization of the protostar embedded within the BHR7 dark cloud, based on both photometric measurements from the near-infrared to millimeter and interferometric continuum and molecular line observations at millimeter wavelengths. We find that this protostar is a Class 0 system, the youngest class of protostars, measuring its bolometric temperature to be 50.5~K, with a bolometric luminosity of 9.3~L$_{\odot}$. The near-infrared and \textit{Spitzer} imaging show a prominent dark lane from dust extinction separating clear bipolar outflow cavities. Observations of $^{13}$CO ($J=2\rightarrow1$), C$^{18}$O ($J=2\rightarrow1$), and other molecular lines with the Submillimeter Array (SMA) exhibit a clear rotation signature on scales $<$1300~AU. The rotation can be traced to an inner radius of $\sim$170~AU and the rotation curve is consistent with an R$^{-1}$ profile, implying that angular momentum is being conserved. Observations of the 1.3~mm dust continuum with the SMA reveal a resolved continuum source, extended in the direction of the dark lane, orthogonal to the outflow. The deconvolved size of the continuum indicates a radius of $\sim$100~AU for the continuum source at the assumed distance of 400~pc. The visibility amplitude profile of the continuum emission cannot be reproduced by an envelope alone and needs a compact component. Thus, we posit that the resolved continuum source could be tracing a Keplerian disk in this very young system. If we assume that the continuum radius traces a Keplerian disk (R$\sim$120~AU) the observed rotation profile is consistent with a protostar mass of 1.0~$M_{\odot}$.

We present 870 $\mu$m ALMA dust polarization observations of 10 young Class 0/I protostars in the Perseus Molecular Cloud. At $\sim$ 0.35$"$ (80 au) resolution, all of our sources show some degree of polarization, with most (9/10) showing significantly extended emission in the polarized continuum. Each source has incredibly intricate polarization signatures. In particular, all three disk-candidates have polarization vectors roughly along the minor axis, which is indicative of polarization produced by dust scattering. On $\sim$ 100 au scales, the polarization is at a relatively low level ($\lesssim 1\%$) and is quite ordered. In sources with significant envelope emission, the envelope is typically polarized at a much higher ($\gtrsim 5\%$) level and has a far more disordered morphology. We compute the cumulative probability distributions for both the small (disk-scale) and large (envelope-scale) polarization percentage. We find that the two are intrinsically different, even after accounting for the different detection thresholds in the high/low surface brightness regions. We perform Kolmogorov-Smirnov and Anderson-Darling tests on the distributions of angle offsets of the polarization from the outflow axis. We find disk-candidate sources are different from the non-disk-candidate sources. We conclude that the polarization on the 100 au scale is consistent with the signature of dust scattering for disk-candidates and that the polarization on the envelope-scale in all sources may come from another mechanism, most likely magnetically aligned grains.

We present a study of hierarchical structure in the Perseus molecular cloud, from the scale of the entire cloud ($\gtrsim$10 pc) to smaller clumps ($\sim$1 pc), cores ($\sim$0.05-0.1 pc), envelopes ($\sim$300-3000 AU) and protostellar objects ($\sim$15 AU). We use new observations from the Submillimeter Array (SMA) large project "Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES)" to probe the envelopes, and recent single-dish and interferometric observations from the literature for the remaining scales. This is the first study to analyze hierarchical structure over five scales in the same cloud complex. We compare the number of fragments with the number of Jeans masses in each scale to calculate the Jeans efficiency, or the ratio of observed to expected number of fragments. The velocity dispersion is assumed to arise either from purely thermal motions, or from combined thermal and non-thermal motions inferred from observed spectral line widths. For each scale, thermal Jeans fragmentation predicts more fragments than observed, corresponding to inefficient thermal Jeans fragmentation. For the smallest scale, thermal plus non-thermal Jeans fragmentation also predicts too many protostellar objects. However at each of the larger scales thermal plus non-thermal Jeans fragmentation predicts fewer than one fragment, corresponding to no fragmentation into envelopes, cores, and clumps. Over all scales, the results are inconsistent with complete Jeans fragmentation based on either thermal or thermal plus non-thermal motions. They are more nearly consistent with inefficient thermal Jeans fragmentation, where the thermal Jeans efficiency increases from the largest to the smallest scale.

The nature of the far-IR source HOPS 108 has been a matter of debate in the last years. Previous radio observations detected a 3.6 cm source (VLA 12), coincident with HOPS 108, that was interpreted as a radio jet from this protostar. We present new multi-wavelength (0.7-5 cm), multi-configuration VLA observations as well as archive data (3.6 cm) that reveal VLA 12 as three knots of non-thermal emission, with HOPS 108 close to the central knot. We show that these knots have not been ejected by HOPS 108. We propose that the VLA 12 knots are actually part of a radio jet driven by VLA 11 (HOPS 370), a strong nearby source clearly elongated in the direction of the knots. The position of HOPS 108 in the path of the VLA 11-VLA 12 jet suggests an appealing new scenario: the triggered formation of HOPS 108 by the interaction of the jet with the surrounding medium.

We present intensity-corrected Herschel maps at 100 um, 160 um, 250 um, 350 um, and 500 um for 56 isolated low-mass clouds. We determine the zero-point corrections for Herschel PACS and SPIRE maps from the Herschel Science Archive (HSA) using Planck data. Since these HSA maps are small, we cannot correct them using typical methods. Here, we introduce a technique to measure the zero-point corrections for small Herschel maps. We use radial profiles to identify offsets between the observed HSA intensities and the expected intensities from Planck. Most clouds have reliable offset measurements with this technique. In addition, we find that roughly half of the clouds have underestimated HSA-SPIRE intensities in their outer envelopes relative to Planck, even though the HSA-SPIRE maps were previously zero-point corrected. Using our technique, we produce corrected Herschel intensity maps for all 56 clouds and determine their line-of-sight average dust temperatures and optical depths from modified black body fits. The clouds have typical temperatures of ~ 14-20 K and optical depths of ~ 1e-5 - 1e-3. Across the whole sample, we find an anti-correlation between temperature and optical depth. We also find lower temperatures than what was measured in previous Herschel studies, which subtracted out a background level from their intensity maps to circumvent the zero-point correction. Accurate Herschel observations of clouds are key to obtain accurate density and temperature profiles. To make such future analyses possible, intensity-corrected maps for all 56 clouds are publicly available in the electronic version.

Low-mass stars build a significant fraction of their total mass during short outbursts of enhanced accretion known as FUor and EXor outbursts. FUor objects are characterized by a sudden brightening of ~5 magnitudes at visible wavelengths within one year and remain bright for decades. EXor objects have lower amplitude outbursts on shorter timescales. Here we discuss a 1.3 mm ALMA mini-survey of eight outbursting sources (three FUor, four EXor, and the borderline object V1647 Ori) in the Orion Molecular Cloud. While previous papers in this series discuss the remarkable molecular outflows observed in the three FUor objects and V1647 Ori, here we focus on the continuum data and the differences and similarities between the FUor and EXor populations. We find that FUor discs are significantly more massive (~80-600 M_JUP) than the EXor objects (~0.5-40 M_JUP ). We also report that the EXor sources lack the prominent outflows seen in the FUor population. Even though our sample is small, the large differences in disc masses and outflow activity suggest that the two types of objects represent different evolutionary stages. The FUor sources seem to be rather compact (Rc < 20-40 au) and to have a smaller characteristic radius for a given disc mass when compared to T Tauri stars. V1118 Ori, the only known close binary systemin our sample, is shown to host a disc around each one of the stellar components. The disc around HBC 494 is asymmetric, hinting at a structure in the outer disc or the presence of a second disc.

Centimeter continuum emission from protostars offers insight into the innermost part of the outflows, as shock-ionized gas produces free-free emission. We observed a complete population of Class 0 and I protostars in the Perseus molecular cloud at 4.1 cm and 6.4 cm with resolution and sensitivity superior to previous surveys. From a total of 71 detections, 8 sources exhibit resolved emission at 4.1 cm and/or 6.4 cm. In this paper we focus on this sub-sample, analyzing their spectral indices along the jet, and their alignment with respect to the large-scale molecular outflow. Spectral indices for fluxes integrated toward the position of the protostar are consistent with free-free thermal emission. The value of the spectral index along a radio jet decreases with distance from the protostar. For six sources, emission is well aligned with the outflow central axis, showing that we observe the ionized base of the jet. This is not the case for two sources, where we note misalignment of the emission with respect to the large-scale outflow. This might indicate that the emission does not originate in the radio jet, but rather in an ionized outflow cavity wall or disk surface. For five of the sources, the spectral indices along the jet decrease well below the thermal free-free limit of -0.1 with $>2\sigma$ significance. This is indicative of synchrotron emission, meaning that high energy electrons are being produced in the outflows close to the disk. This result can have far-reaching implications for the chemical composition of the embedded disks.

We present a high angular resolution ($\sim 0.2^{\prime\prime}$), high sensitivity ($\sigma \sim 0.2$ mJy) survey of the 870 $\mu$m continuum emission from the circumstellar material around 49 pre-main sequence stars in the $\rho$ Ophiuchus molecular cloud. Because most millimeter instruments have resided in the northern hemisphere, this represents the largest high-resolution, millimeter-wave survey of the circumstellar disk content of this cloud. Our survey of 49 systems comprises 63 stars; we detect disks associated with 29 single sources, 11 binaries, 3 triple systems and 4 transition disks. We present flux and radius distributions for these systems; in particular, this is the first presentation of a reasonably complete probability distribution of disk radii at millimeter-wavelengths. We also compare the flux distribution of these protoplanetary disks with that of the disk population of the Taurus-Auriga molecular cloud. We find that disks in binaries are both significantly smaller and have much less flux than their counterparts around isolated stars. We compute truncation calculations on our binary sources and find that these disks are too small to have been affected by tidal truncation and posit some explanations for this. Lastly, our survey found 3 candidate gapped disks, one of which is a newly identified transition disk with no signature of a dip in infrared excess in extant observations.

The youngest low-mass protostars are known to be chemically rich, accreting matter most vigorously, and producing the most powerful outflows. Molecules are unique tracers of these phenomena. We use ALMA to study several outflow sources in the Serpens Main region. The most luminous source, Ser-SMM1, shows the richest chemical composition, but some complex molecules are also present in S68N. No emission from complex organics is detected toward Ser-emb 8N, which is the least luminous in the sample. We discuss whether these differences reflect an evolutionary effect or whether they are due to different physical structures. We also analyze the outflow structure from these young protostars by comparing emission of CO and SiO. EHV molecular jets originating from SMM1-a,b and Ser-emb 8N contrast with no such activity from S68N, which on the other hand presents a complex outflow structure.

Planets form in disks around young stars. The planet formation process may start when the protostar and disk are still deeply embedded within their infalling envelope. However, unlike more evolved protoplanetary disks, the physical and chemical structure of these young embedded disks are still poorly constrained. We have analyzed ALMA data for $^{13}$CO, C$^{18}$O and N$_2$D$^+$ to constrain the temperature structure, one of the critical unknowns, in the disk around L1527. The spatial distribution of $^{13}$CO and C$^{18}$O, together with the kinetic temperature derived from the optically thick $^{13}$CO emission and the non-detection of N$_2$D$^+$, suggest that this disk is warm enough ($\gtrsim$ 20 K) to prevent CO freeze-out.

We present Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of the star-forming environment surrounding V1647 Ori, an outbursting FUor/EXor pre-MS star. Dust continuum and the (J = 2 - 1) $^{12}$CO, $^{13}$CO, C$^{18}$O molecular emission lines were observed to characterize the V1647 Ori circumstellar disc and any large scale molecular features present. We detect continuum emission from the circumstellar disc and determine a radius r = 40 au, inclination i = 17$^{\circ}$$^{+6}_{-9}$ and total disc mass of M$_{\mathrm{disk}}$ of ~0.1 M$_{\odot}$. We do not identify any disc structures associated with nearby companions, massive planets or fragmentation. The molecular cloud environment surrounding V1647 Ori is both structured and complex. We confirm the presence of an excavated cavity north of V1647 Ori and have identified dense material at the base of the optical reflection nebula (McNeil's Nebula) that is actively shaping its surrounding environment. Two distinct outflows have been detected with dynamical ages of ~11,700 and 17,200 years. These outflows are misaligned suggesting disc precession over ~5500 years as a result of anisotropic accretion events is responsible. The collimated outflows exhibit velocities of ~2 km s$^{-1}$, similar in velocity to that of other FUor objects presented in this series but significantly slower than previous observations and model predictions. The V1647 Ori system is seemingly connected by an "arm" of material to a large unresolved structure located ~20$"$ to the west. The complex environment surrounding V1647 Ori suggests it is in the early stages of star formation which may relate to its classification as both an FUor and EXor type object.

The kinematics and dynamics of young stellar populations enable us to test theories of star formation. With this aim, we continue our analysis of the SDSS-III/APOGEE IN-SYNC survey, a high resolution near infrared spectroscopic survey of young clusters. We focus on the Orion A star-forming region, for which IN-SYNC obtained spectra of $\sim2700$ stars. In Paper IV we used these data to study the young stellar population. Here we study the kinematic properties through radial velocities ($v_r$). The young stellar population remains kinematically associated with the molecular gas, following a $\sim10\:{\rm{km\:s}}^{-1}$ gradient along filament. However, near the center of the region, the $v_r$ distribution is slightly blueshifted and asymmetric; we suggest that this population, which is older, is slightly in foreground. We find evidence for kinematic subclustering, detecting statistically significant groupings of co-located stars with coherent motions. These are mostly in the lower-density regions of the cloud, while the ONC radial velocities are smoothly distributed, consistent with it being an older, more dynamically evolved cluster. The velocity dispersion $\sigma_v$ varies along the filament. The ONC appears virialized, or just slightly supervirial, consistent with an old dynamical age. Here there is also some evidence for on-going expansion, from a $v_r$--extinction correlation. In the southern filament, $\sigma_v$ is $\sim2$--$3$ times larger than virial in the L1641N region, where we infer a superposition along the line of sight of stellar sub-populations, detached from the gas. On the contrary, $\sigma_v$ decreases towards L1641S, where the population is again in agreement with a virial state.

We present new Submillimeter Array (SMA) observations of CO(2-1) outflows toward young, embedded protostars in the Perseus molecular cloud as part of the Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES) survey. For 57 Perseus protostars, we characterize the orientation of the outflow angles and compare them with the orientation of the local filaments as derived from $Herschel$ observations. We find that the relative angles between outflows and filaments are inconsistent with purely parallel or purely perpendicular distributions. Instead, the observed distribution of outflow-filament angles are more consistent with either randomly aligned angles or a mix of projected parallel and perpendicular angles. A mix of parallel and perpendicular angles requires perpendicular alignment to be more common by a factor of $\sim$3. Our results show that the observed distributions probably hold regardless of the protostar's multiplicity, age, or the host core's opacity. These observations indicate that the angular momentum axis of a protostar may be independent of the large-scale structure. We discuss the significance of independent protostellar rotation axes in the general picture of filament-based star formation.

The Herschel Orion Protostar Survey obtained well-sampled 1.2 - 870 micron spectral energy distributions (SEDs) of over 300 protostars in the Orion molecular clouds, home to most of the young stellar objects (YSOs) in the nearest 500 pc. We plot the bolometric luminosities and temperatures for 330 Orion YSOs, 315 of which have bolometric temperatures characteristic of protostars. The histogram of bolometric temperature is roughly flat; 29% of the protostars are in Class 0. The median luminosity decreases by a factor of four with increasing bolometric temperature; consequently, the Class 0 protostars are systematically brighter than the Class I protostars, with a median luminosity of 2.3 L_sun as opposed to 0.87 L_sun. At a given bolometric temperature, the scatter in luminosities is three orders of magnitude. Using fits to the SEDs, we analyze how the luminosities corrected for inclination and foreground reddening relate to the mass in the inner 2500 AU of the best-fit model envelopes. The histogram of envelope mass is roughly flat, while the median corrected luminosity peaks at 15 L_sun for young envelopes and falls to 1.7 L_sun for late-stage protostars with remnant envelopes. The spread in luminosity at each envelope mass is three orders of magnitude. Envelope masses that decline exponentially with time explain the flat mass histogram and the decrease in luminosity, while the formation of a range of stellar masses explains the dispersion in luminosity.

Context: Understanding how accretion proceeds is a key question of star formation, with important implications for both the physical and chemical evolution of young stellar objects. In particular, very little is known about the accretion variability in the earliest stages of star formation. Aims: To characterise protostellar accretion histories towards individual sources by utilising sublimation and freeze-out chemistry of CO. Methods: A sample of 24 embedded protostars are observed with the Submillimeter Array (SMA) in context of the large program "Mass Assembly of Stellar Systems and their Evolution with the SMA" (MASSES). The size of the C$^{18}$O emitting region, where CO has sublimated into the gas-phase, is measured towards each source and compared to the expected size of the region given the current luminosity. The SMA observations also include 1.3 mm continuum data, which are used to investigate whether a link can be established between accretion bursts and massive circumstellar disks. Results: Depending on the adopted sublimation temperature of the CO ice, between 20% and 50% of the sources in the sample show extended C$^{18}$O emission indicating that the gas was warm enough in the past that CO sublimated and is currently in the process of refreezing; something which we attribute to a recent accretion burst. Given the fraction of sources with extended C$^{18}$O emission, we estimate an average interval between bursts of 20000-50000 yr, which is consistent with previous estimates. No clear link can be established between the presence of circumstellar disks and accretion bursts, however the three closest known binaries in the sample (projected separations <20 AU) all show evidence of a past accretion burst, indicating that close binary interactions may also play a role in inducing accretion variability.

We carried out multiwavelength (0.7-5 cm), multiepoch (1994-2015) Very Large Array (VLA) observations toward the region enclosing the bright far-IR sources FIR 3 (HOPS 370) and FIR 4 (HOPS 108) in OMC-2. We report the detection of 10 radio sources, seven of them identified as young stellar objects. We image a well-collimated radio jet with a thermal free-free core (VLA 11) associated with the Class I intermediate-mass protostar HOPS 370. The jet presents several knots (VLA 12N, 12C, 12S) of non-thermal radio emission (likely synchrotron from shock-accelerated relativistic electrons) at distances of ~7,500-12,500 au from the protostar, in a region where other shock tracers have been previously identified. These knots are moving away from the HOPS 370 protostar at ~ 100 km/s. The Class 0 protostar HOPS 108, which itself is detected as an independent, kinematically decoupled radio source, falls in the path of these non-thermal radio knots. These results favor the previously proposed scenario where the formation of HOPS 108 has been triggered by the impact of the HOPS 370 outflow with a dense clump. However, HOPS 108 presents a large proper motion velocity of ~ 30 km/s, similar to that of other runaway stars in Orion, whose origin would be puzzling within this scenario. Alternatively, an apparent proper motion could result because of changes in the position of the centroid of the source due to blending with nearby extended emission, variations in the source shape, and /or opacity effects.

We present Atacama Large Millimeter/ sub-millimeter Array (ALMA) observations of V883 Ori, an FU Ori object. We describe the molecular outflow and envelope of the system based on the $^{12}$CO and $^{13}$CO emissions, which together trace a bipolar molecular outflow. The C$^{18}$O emission traces the rotational motion of the circumstellar disk. From the $^{12}$CO blue-shifted emission, we estimate a wide opening angle of $\sim$ 150$^{^{\circ}}$ for the outflow cavities. Also, we find that the outflow is very slow (characteristic velocity of only 0.65 km~s$^{-1}$), which is unique for an FU Ori object. We calculate the kinematic properties of the outflow in the standard manner using the $^{12}$CO and $^{13}$CO emissions. In addition, we present a P Cygni profile observed in the high-resolution optical spectrum, evidence of a wind driven by the accretion and being the cause for the particular morphology of the outflows. We discuss the implications of our findings and the rise of these slow outflows during and/or after the formation of a rotationally supported disk.

Compact substructure is expected to arise in a starless core as mass becomes concentrated in the central region likely to form a protostar. Additionally, multiple peaks may form if fragmentation occurs. We present ALMA Cycle 2 observations of 60 starless and protostellar cores in the Ophiuchus molecular cloud. We detect eight compact substructures which are >15 arcsec from the nearest Spitzer YSO. Only one of these has strong evidence for being truly starless after considering ancillary data, e.g., from Herschel and X-ray telescopes. An additional extended emission structure has tentative evidence for starlessness. The number of our detections is consistent with estimates from a combination of synthetic observations of numerical simulations and analytical arguments. This result suggests that a similar ALMA study in the Chamaeleon I cloud, which detected no compact substructure in starless cores, may be due to the peculiar evolutionary state of cores in that cloud.

We present Atacama Large Millimeter/sub-millimeter Array (ALMA) Cycle-2 observations of the HBC 494 molecular outflow and envelope. HBC 494 is an FU Ori-like object embedded in the Orion A cloud and is associated with the reflection nebulae Re50 and Re50N. We use $^{12}$CO, $^{13}$CO and C$^{18}$O spectral line data to independently describe the outflow and envelope structures associated with HBC 494. The moment-1 map of the $^{12}$CO emission shows the widest outflow cavities in a Class I object known to date (opening angle $\sim$ 150$^{^{\circ}}$). The morphology of the wide outflow is likely to be due to the interaction between winds originating in the inner disc and the surrounding envelope. The low-velocity blue- and red-shifted $^{13}$CO and C$^{18}$O emission trace the rotation and infall motion of the circumstellar envelope. Using molecular line data and adopting standard methods for correcting optical depth effects, we estimate its kinematic properties, including an outflow mass on the order of 10$^{-1}$ M$_{\odot}$. Considering the large estimated outflow mass for HBC 494, our results support recent theoretical work suggesting that wind-driven processes might dominate the evolution of protoplanetary discs via energetic outflows.

We present the first results of the Gould's Belt Distances Survey (GOBELINS), a project aimed at measuring the proper motion and trigonometric parallax of a large sample of young stars in nearby regions using multi-epoch Very Long Baseline Array (VLBA) radio observations. Enough VLBA detections have now been obtained for 16 stellar systems in Ophiuchus to derive their parallax and proper motion. This leads to distance determinations for individual stars with an accuracy of 0.3 to a few percent. In addition, the orbits of 6 multiple systems were modelled by combining absolute positions with VLBA (and in some cases, near infrared) angular separations. Twelve stellar systems are located in the dark cloud Lynds 1688, the individual distances for this sample are highly consistent with one another, and yield a mean parallax for Lynds 1688 of $\varpi=7.28\pm0.06$ mas, corresponding to a distance $d=137.3\pm1.2$ pc. This represents an accuracy better than 1%. Three systems for which astrometric elements could be measured are located in the eastern streamer (Lynds 1689) and yield an estimate of $\varpi=6.79\pm0.16$ mas, corresponding to a distance $d=147.3\pm3.4$ pc. This suggests that the eastern streamer is located about 10 pc farther than the core, but this conclusion needs to be confirmed by observations (currently being collected) of additional sources in the eastern streamer. From the measured proper motions, we estimate the one-dimensional velocity dispersion in Lynds 1688 to be 2.8$\pm$1.8 and 3.0$\pm$2.0 ${\rm km~s}^{-1}$, in R.A. and DEC., respectively, these are larger than, but still consistent within $1\sigma$, with those found in other studies.

Using multi-epoch VLA observations, covering a time baseline of 29.1 years, we have measured the proper motions of 88 young stars with compact radio emission in the core of the Orion Nebula Cluster (ONC) and the neighboring BN/KL region. Our work increases the number of young stars with measured proper motion at radio frequencies by a factor of 2.5 and enables us to perform a better statistical analysis of the kinematics of the region than was previously possible. Most stars (79 out of 88) have proper motions consistent with a Gaussian distribution centered on $\overline{\mu_{\alpha}\cos{\delta}}=1.07\pm0.09\quad{\rm mas\,yr^{-1}}$, and $\overline{\mu_{\delta}}=-0.84\pm0.16\quad{\rm mas\,yr^{-1}}$, with velocity dispersions of $\sigma_{\alpha}=1.08\pm0.07\quad{\rm mas\,\,yr^{-1}},$ $\sigma_{\delta}=1.27\pm0.15\quad{\rm mas\,\,yr^{-1}}$. We looked for organized movements of these stars but found no clear indication of radial expansion/contraction or rotation. The remaining nine stars in our sample show peculiar proper motions that differ from the mean proper motions of the ONC by more than 3-$\sigma$. One of these stars, V 1326 Ori, could have been expelled from the Orion Trapezium 7,000 years ago. Two could be related to the multi-stellar disintegration in the BN/KL region, in addition to the previously known sources BN, I and n. The others either have high uncertainties (so their anomalous proper motions are not firmly established) or could be foreground objects.

As part of an ALMA survey to study the origin of episodic accretion in young eruptive variables, we have observed the circumstellar environment of the star V2775 Ori. This object is a very young, pre-main sequence object which displays a large amplitude outburst characteristic of the FUor class. We present Cycle-2 band 6 observations of V2775 Ori with a continuum and CO (2-1) isotopologue resolution of 0.25\as (103 au). We report the detection of a marginally resolved circumstellar disc in the ALMA continuum with an integrated flux of $106 \pm 2$ mJy, characteristic radius of $\sim$ 30 au, inclination of $14.0^{+7.8}_{-14.5}$ deg, and is oriented nearly face-on with respect to the plane of the sky. The \co~emission is separated into distinct blue and red-shifted regions that appear to be rings or shells of expanding material from quasi-episodic outbursts. The system is oriented in such a way that the disc is seen through the outflow remnant of V2775 Ori, which has an axis along our line-of-sight. The $^{13}$CO emission displays similar structure to that of the \co, while the C$^{18}$O line emission is very weak. We calculated the expansion velocities of the low- and medium-density material with respect to the disc to be of -2.85 km s$^{-1}$ (blue), 4.4 km s$^{-1}$ (red) and -1.35 and 1.15 km s$^{-1}$ (for blue and red) and we derived the mass, momentum and kinetic energy of the expanding gas. The outflow has an hourglass shape where the cavities are not seen. We interpret the shapes that the gas traces as cavities excavated by an ancient outflow. We report a detection of line emission from the circumstellar disc and derive a lower limit of the gas mass of 3 \MJup.

Binary and multiple star systems are a frequent outcome of the star formation process, and as a result, almost half of all sun-like stars have at least one companion star. Theoretical studies indicate that there are two main pathways that can operate concurrently to form binary/multiple star systems: large scale fragmentation of turbulent gas cores and filaments or smaller scale fragmentation of a massive protostellar disk due to gravitational instability. Observational evidence for turbulent fragmentation on scales of $>$1000~AU has recently emerged. Previous evidence for disk fragmentation was limited to inferences based on the separations of more-evolved pre-main sequence and protostellar multiple systems. The triple protostar system L1448 IRS3B is an ideal candidate to search for evidence of disk fragmentation. L1448 IRS3B is in an early phase of the star formation process, likely less than 150,000 years in age, and all protostars in the system are separated by $<$200~AU. Here we report observations of dust and molecular gas emission that reveal a disk with spiral structure surrounding the three protostars. Two protostars near the center of the disk are separated by 61 AU, and a tertiary protostar is coincident with a spiral arm in the outer disk at a 183 AU separation. The inferred mass of the central pair of protostellar objects is $\sim$1 M$_{sun}$, while the disk surrounding the three protostars has a total mass of $\sim$0.30 M$_{\sun}$. The tertiary protostar itself has a minimum mass of $\sim$0.085 M$_{sun}$. We demonstrate that the disk around L1448 IRS3B appears susceptible to disk fragmentation at radii between 150~AU and 320~AU, overlapping with the location of the tertiary protostar. This is consistent with models for a protostellar disk that has recently undergone gravitational instability, spawning one or two companion stars.

We report on new distances and proper motions to seven stars across the Serpens/Aquila complex. The observations were obtained as part of the Gould's Belt Distances Survey (GOBELINS) project between September 2013 and April 2016 with the Very Long Baseline Array (VLBA). One of our targets is the proto-Herbig AeBe object EC 95, which is a binary system embedded in the Serpens Core. For this system, we combined the GOBELINS observations with previous VLBA data to cover a total period of ~8 years, and derive the orbital elements and an updated source distance. The individual distances to sources in the complex are fully consistent with each other, and the mean value corresponds to a distance of $436.0\pm9.2$~pc for the Serpens/W40 complex. Given this new evidence, we argue that Serpens Main, W40 and Serpens South are physically associated and form a single cloud structure.

We present the results of the Gould's Belt Distances Survey (GOBELINS) of young star forming regions towards the Orion Molecular Cloud Complex. We detected 36 YSOs with the Very Large Baseline Array (VLBA), 27 of which have been observed in at least 3 epochs over the course of 2 years. At least half of these YSOs belong to multiple systems. We obtained parallax and proper motions towards these stars to study the structure and kinematics of the Complex. We measured a distance of 388$\pm$5 pc towards the Orion Nebula Cluster, 428$\pm$10 pc towards the southern portion L1641, 388$\pm$10 pc towards NGC 2068, and roughly $\sim$420 pc towards NGC 2024. Finally, we observed a strong degree of plasma radio scattering towards $\lambda$ Ori.

Protostellar evolution, following the formation of the protostar is becoming reasonably well characterized, but the evolution from a prestellar core to a protostar is not well known, although the first hydrostatic core (FHSC) must be a pivotal step. NGC 1333 IRAS 4C is a potentially very young object, that we directly compare with the nearby Class 0 IRAS 4A and IRAS 4B. Observational constraints are provided by spectral imaging from the JCMT Spectral Legacy Survey (330-373 GHz) and continuum and line observations from CARMA. We present integrated intensity and velocity maps of several species, including CO, H2CO and CH3OH. The velocity of an observed outflow, the degree of CO depletion, the deuteration of DCO+/HCO+ and gas kinetic temperatures are observational signatures that we present. We report differences between the three sources in four aspects: a) the kinetic temperature is much lower towards IRAS 4C, b) the line profiles of the detected species show strong outflow activity towards IRAS 4A and IRAS 4B but not towards IRAS 4C, c) the HCN/HNC is < 1 towards IRAS 4C, which confirms the cold nature of the source, d) the degree of CO depletion and the deuteration are the lowest towards the warmest of the sources, IRAS 4B. IRAS 4C seems to be in a different evolutionary state than the IRAS 4A and IRAS 4B sources. We can probably exclude the FHSC stage due to the relatively low Lsmm/Lbol (~6\%) and we investigate the earliest accretion phase of Class 0 stage and the transition between Class 0 to Class I. Our results do not show a consistent scenario for either case, with the major issue being the absence of outflow activity and the cold nature of IRAS 4C. The number of FHSC candidates in Perseus is ~10 times higher than current models predict, which suggests that the lifespan of these objects is greater than 1000 years, possibly due to an accretion rate lower than 4x10^-5 msun/yr.

We present the first detection of a jet in the far-IR [O I] lines from an intermediate mass protostar. We have carried out a Herschel/PACS spectral mapping study in the [O I] lines of OMC-2 FIR 3 and FIR 4, two of the most luminous protostars in Orion outside of the Orion Nebula. The spatial morphology of the fine structure line emission reveals the presence of an extended photodissociation region (PDR) and a narrow, but intense jet connecting the two protostars. The jet seen in [O I] emission is spatially aligned with the Spitzer/IRAC 4.5 micron jet and the CO (6-5) molecular outflow centered on FIR 3. The mass loss rate derived from the total [O I] 63 micron line luminosity of the jet is 7.7 x 10^-6 M_sun/yr, more than an order of magnitude higher than that measured for typical low mass class 0 protostars. The implied accretion luminosity is significantly higher than the observed bolometric luminosity of FIR 4, indicating that the [O I] jet is unlikely to be associated with FIR 4. We argue that the peak line emission seen toward FIR 4 originates in the terminal shock produced by the jet driven by FIR 3. The higher mass-loss rate that we find for FIR 3 is consistent with the idea that intermediate mass protostars drive more powerful jets than their low-mass counterparts. Our results also call into question the nature of FIR 4.

We investigate the evolution of far-IR CO emission from protostars observed with Herschel/PACS for 50 sources from the combined sample of HOPS and DIGIT Herschel key programs. From the uniformly sampled spectral energy distributions, we computed $L_{\rm{bol}}$, $T_{\rm{bol}}$ and $L_{\rm {bol}}/L_{\rm {smm}}$ for these sources to search for correlations between far-IR CO emission and protostellar properties. We find a strong and tight correlation between far-IR CO luminosity ($L^{\rm fir}_{\rm CO}$) and the bolometric luminosity ($L_{\rm{bol}}$) of the protostars with $L^{\rm fir}_{\rm CO}$ $\propto$ $L_{\rm{bol}}^{0.7}$. We, however, do not find a strong correlation between $L^{\rm fir}_{\rm CO}$ and protostellar evolutionary indicators, $T_{\rm{bol}}$ and $L_{\rm {bol}}/L_{\rm {smm}}$. FIR CO emission from protostars traces the currently shocked gas by jets/outflows, and $L^{\rm fir}_{\rm CO}$ is proportional to the instantaneous mass loss rate, $\dot{M}_{\rm{out}}$. The correlation between $L^{\rm fir}_{\rm CO}$ and $L_{\rm{bol}}$ is indicative of instantaneous $\dot{M}_{\rm{out}}$ tracking instantaneous $\dot{M}_{\rm{acc}}$. The lack of correlation between $L^{\rm fir}_{\rm CO}$ and evolutionary indicators $T_{\rm{bol}}$ and $L_{\rm {bol}}/L_{\rm {smm}}$ suggests that $\dot{M}_{\rm{out}}$ and, therefore, $\dot{M}_{\rm{acc}}$ do not show any clear evolutionary trend. These results are consistent with mass accretion/ejection in protostars being episodic. Taken together with the previous finding that the time-averaged mass ejection/accretion rate declines during the protostellar phase (e.g., Bontemps et al. 1996), our results suggest that the instantaneous accretion/ejection rate of protostars is highly time variable and episodic, but the amplitude and/or frequency of this variability decreases with time such that the time averaged accretion/ejection rate declines with system age.

A snow-line is the region of a protoplanetary disk at which a major volatile, such as water or carbon monoxide, reaches its condensation temperature. Snow-lines play a crucial role in disk evolution by promoting the rapid growth of ice-covered grains. Signatures of the carbon monoxide snow-line (at temperatures of around 20 kelvin) have recently been imaged in the disks surrounding the pre-main-sequence stars TW Hydra and HD163296, at distances of about 30 astronomical units (au) from the star. But the water snow-line of a protoplanetary disk (at temperatures of more than 100 kelvin) has not hitherto been seen, as it generally lies very close to the star (less than 5 au away for solar-type stars). Water-ice is important because it regulates the efficiency of dust and planetesimal coagulation, and the formation of comets, ice giants and the cores of gas giants. Here we report ALMA images at 0.03-arcsec resolution (12 au) of the protoplanetary disk around V883 Ori, a protostar of 1.3 solar masses that is undergoing an outburst in luminosity arising from a temporary increase in the accretion rate. We find an intensity break corresponding to an abrupt change in the optical depth at about 42 au, where the elevated disk temperature approaches the condensation point of water, from which we conclude that the outburst has moved the water snow-line. The spectral behaviour across the snow-line confirms recent model predictions: dust fragmentation and the inhibition of grain growth at higher temperatures results in soaring grain number densities and optical depths. As most planetary systems are expected to experience outbursts caused by accretion during their formation our results imply that highly dynamical water snow-lines must be considered when developing models of disk evolution and planet formation.

There is substantial evidence for disk formation taking place during the early stages of star formation and for most stars being born in multiple systems; however, protostellar multiplicity and disk searches have been hampered by low resolution, sample bias, and variable sensitivity. We have conducted an unbiased, high-sensitivity Karl G. Jansky Very Large Array (VLA) survey toward all known protostars (n = 94) in the Perseus molecular cloud (d~230 pc), with a resolution of ~15 AU (0.06") at 8 mm. We have detected candidate protostellar disks toward 17 sources (with 12 of those in the Class 0 stage) and we have found substructure on < 50AU scales for three Class 0 disk candidates, possibly evidence for disk fragmentation. We have discovered 16 new multiple systems (or new components) in this survey; the new systems have separations < 500 AU and 3 by < 30 AU. We also found a bi-modal distribution of separations, with peaks at ~75 AU and ~3000 AU, suggestive of formation through two distinct mechanisms: disk and turbulent fragmentation. The results from this survey demonstrate the necessity and utility of uniform, unbiased surveys of protostellar systems at millimeter and centimeter wavelengths.

We present CARMA CO (J=1-0) observations and Herschel PACS spectroscopy, characterizing the outflow properties toward extremely young and deeply embedded protostars in the Orion molecular clouds. The sample comprises a subset of the Orion protostars known as the PACS Bright Red Sources (PBRS) (Stutz et al. 2013). We observed 14 PBRS with CARMA and 8 of these 14 with Herschel, acquiring full spectral scans from 55 micron to 200 micron. Outflows are detected in CO (J=1-0) from 8 of 14 PBRS, with two additional tentative detections; outflows are also detected from the outbursting protostar HOPS 223 (V2775 Ori) and the Class I protostar HOPS 68. The outflows have a range of morphologies, some are spatially compact, <10000 AU in extent, while others extend beyond the primary beam. The outflow velocities and morphologies are consistent with being dominated by intermediate inclination angles (80 deg > i > 20 deg). This confirms the interpretation of the very red 24 micron to 70 micron colors of the PBRS as a signpost of high envelope densities, with only one (possibly two) cases of the red colors resulting from edge-on inclinations. We detect high-J (J_up > 13) CO lines and/or H_2O lines from 5 of 8 PBRS and only for those with detected CO outflows. The far-infrared CO rotation temperatures of the detected PBRS are marginally colder (~230 K) than those observed for most protostars (~300 K), and only one of these 5 PBRS has detected [OI] 63 micron emission. The high envelope densities could be obscuring some [OI] emission and cause a ~20 K reduction to the CO rotation temperatures.

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.

We present a far-IR survey of the entire Mon R2 GMC with $Herschel-SPIRE$ cross-calibrated with $Planck-HFI$ data. We fit the SEDs of each pixel with a greybody function and an optimal beta value of 1.8. We find that mid-range column densities obtained from far-IR dust emission and near-IR extinction are consistent. For the entire GMC, we find that the column density histogram, or N-PDF, is lognormal below $\sim$10$^{21}$ cm$^{-2}$. Above this value, the distribution takes a power law form with an index of -2.16. We analyze the gas geometry, N-PDF shape, and YSO content of a selection of subregions in the cloud. We find no regions with pure lognormal N-PDFs. The regions with a combination of lognormal and one power law N-PDF have a YSO cluster and a corresponding centrally concentrated gas clump. The regions with a combination of lognormal and two power law N-PDF have significant numbers of typically younger YSOs but no prominent YSO cluster. These regions are composed of an aggregate of closely spaced gas filaments with no concentrated dense gas clump. We find that for our fixed scale regions, the YSO count roughly correlates with the N-PDF power law index. The correlation appears steeper for single power law regions relative to two power law regions with a high column density cut-off, as a greater dense gas mass fraction is achieved in the former. A stronger correlation is found between embedded YSO count and the dense gas mass among our regions.

Dust polarization in millimeter (and centimeter) has been mapped in disks around an increasing number of young stellar objects. It is usually thought to come from emission by magnetically aligned (non-spherical) grains, but can also be produced by dust scattering. We present a semi-analytic theory of disk polarization that includes both the direction emission and scattering, with an emphasis on their relative importance and how they are affected by the disk inclination. For face-on disks, both emission and scattering tend to produce polarization in the radial direction, making them difficult to distinguish, although the scattering-induced polarization can switch to the azimuthal direction if the incident radiation is beamed strongly enough in the radial direction in the disk plane. Disk inclination affects the polarizations from emission and scattering differently, especially on the major axis where, in the edge-on limit, the former vanishes while the latter reaches a polarization fraction as large as 1/3. The polarizations from the two competing mechanisms tend to cancel each other on the major axis, producing two low polarization "holes" (one on each side of the center) under certain conditions. We find tantalizing evidence for at least one such "hole" in NGC1333 IRAS4A1, whose polarization observed at 8 mm on the 100 AU scale is indicative of a pattern dominated by scattering close to the center and by direction emission in the outer region. If true, it would imply not only that a magnetic field exists on the disk scale, but that it is strong enough to align large, possibly mm-sized, grains.

Multiplicity is common in field stars and among protostellar systems. Models suggest two paths of formation: turbulent fragmentation and protostellar disk fragmentation. We attempt to find whether or not the coevality frequency of multiple protostellar systems can help to better understand their formation mechanism. The coevality frequency is determined by constraining the relative evolutionary stages of the components in a multiple system. SEDs for known multiple protostars in Perseus were constructed from literature data. Herschel PACS photometric maps were used to sample the peak of the SED for systems with separations >7", a crucial aspect in determining the evolutionary stage of a protostellar system. Inclination effects and the surrounding envelope and outflows were considered to decouple source geometry from evolution. This together with the shape and derived properties from the SED was used to determine each system's coevality as accurately as possible. SED models were used to examine the frequency of non-coevality that is due to geometry. We find a non-coevality frequency of 33+/-10% from the comparison of SED shapes of resolved multiple systems. Other source parameters suggest a somewhat lower frequency of non-coevality. The frequency of apparent non-coevality that is due to random inclination angle pairings of model SEDs is 17+/-0.5%. Observations of the outflow of resolved multiple systems do not suggest significant misalignments within multiple systems. Effects of unresolved multiples on the SED shape are also investigated. We find that 1/3 of the multiple protostellar systems sampled here are non-coeval, which is more than expected from random geometric orientations. The other 2/3 are found to be coeval. Higher order multiples show a tendency to be non-coeval. The frequency of non-coevality found here is most likely due to formation and enhanced by dynamical evolution.

We present an Atacama Large Millimeter/submillimeter Array (ALMA) 106 GHz (Band 3) continuum survey of the complete population of dense cores in the Chamaeleon I molecular cloud. We detect a total of 24 continuum sources in 19 different target fields. All previously known Class 0 and Class I protostars in Chamaeleon I are detected, whereas all of the 56 starless cores in our sample are undetected. We show that the Spitzer+Herschel census of protostars in Chamaeleon I is complete, with the rate at which protostellar cores have been misclassified as starless cores calculated as <1/56, or < 2%. We use synthetic observations to show that starless cores collapsing following the turbulent fragmentation scenario are detectable by our ALMA observations when their central densities exceed ~10^8 cm^-3, with the exact density dependent on the viewing geometry. Bonnor-Ebert spheres, on the other hand, remain undetected to central densities at least as high as 10^10 cm^-3. Our starless core non-detections are used to infer that either the star formation rate is declining in Chamaeleon I and most of the starless cores are not collapsing, matching the findings of previous studies, or that the evolution of starless cores are more accurately described by models that develop less substructure than predicted by the turbulent fragmentation scenario, such as Bonnor-Ebert spheres. We outline future work necessary to distinguish between these two possibilities.

We present key results from the Herschel Orion Protostar Survey (HOPS): spectral energy distributions (SEDs) and model fits of 330 young stellar objects, predominantly protostars, in the Orion molecular clouds. This is the largest sample of protostars studied in a single, nearby star-formation complex. With near-infrared photometry from 2MASS, mid- and far-infrared data from Spitzer and Herschel, and sub-millimeter photometry from APEX, our SEDs cover 1.2-870 $\mu$m and sample the peak of the protostellar envelope emission at ~100 $\mu$m. Using mid-IR spectral indices and bolometric temperatures, we classify our sample into 92 Class 0 protostars, 125 Class I protostars, 102 flat-spectrum sources, and 11 Class II pre-main-sequence stars. We implement a simple protostellar model (including a disk in an infalling envelope with outflow cavities) to generate a grid of 30400 model SEDs and use it to determine the best-fit model parameters for each protostar. We argue that far-IR data are essential for accurate constraints on protostellar envelope properties. We find that most protostars, and in particular the flat-spectrum sources, are well-fit. The median envelope density and median inclination angle decrease from Class 0 to Class I to flat-spectrum protostars, despite the broad range in best-fit parameters in each of the three categories. We also discuss degeneracies in our model parameters. Our results confirm that the different protostellar classes generally correspond to an evolutionary sequence with a decreasing envelope infall rate, but the inclination angle also plays a role in the appearance, and thus interpretation, of the SEDs.

(Abridged) The physical structure of deeply-embedded low-mass protostars (Class 0) on scales of less than 300 AU is still poorly constrained. Determining this is crucial for understanding the physical and chemical evolution from cores to disks. In this study two models of the emission, a Gaussian disk intensity distribution and a parametrized power-law disk model, are fitted to sub-arcsecond resolution interferometric continuum observations of five Class 0 sources, including one source with a confirmed Keplerian disk. For reference, a spherically symmetric single power-law envelope is fitted to the larger scale ($\sim$1000 AU) emission and investigated further for one of the sources on smaller scales. A thin disk model can approximate the emission and physical structure in the inner few 100 AU scales of the studied deeply-embedded low-mass protostars and paves the way for analysis of a larger sample with ALMA. While the disk radii agree with previous estimates the masses are different for some of the sources studied. Assuming a typical temperature distribution, the fractional amount of mass in the disk above 100 K varies in between 7% to 30%. Kinematic data are needed to determine the presence of any Keplerian disk. Using previous observations of p-H$_2^{18}$O, we estimate the relative gas phase water abundances roughly an order of magnitude higher than previously inferred when both warm and cold H$_2$ was used as reference. A spherically symmetric single power-law envelope model fails to simultaneously reproduce both the small and large scale emission.

We investigate the alignment between outflow axes in nine of the youngest binary/multiple systems in the Perseus Molecular Cloud. These systems have typical member spacing larger than 1000 AU. For outflow identification, we use 12CO(2-1) and 12CO(3-2) data from a large survey with the Sub- millimeter Array: Mass Assembly of Stellar Systems and their Evolution with the SMA (MASSES). The distribution of outflow orientations in the binary pairs is consistent with random or preferentially anti-aligned distributions, demonstrating that these outflows are misaligned. This result suggests that these systems are possibly formed in environments where the distribution of angular momentum is complex and disordered, and these systems do not come from the same co-rotating structures or from an initial cloud with aligned vectors of angular momentum.

We examine the spectroscopic binary population for two massive nearby regions of clustered star formation, the Orion Nebula Cluster and NGC 2264, supplementing the data presented by Tobin et al. (2009, 2015) with more recent observations and more extensive analysis. The inferred multiplicity fraction up to 10 AU based on these observations is $5.3\pm 1.2$% for NGC 2264 and $5.8\pm 1.1$% for the ONC; they are consistent with the distribution of binaries in the field in the relevant parameter range. Eight of the multiple systems in the sample have enough epochs to make an initial fit for the orbital parameters. Two of these sources are double-lined spectroscopic binaries; for them we determine the mass ratio. Our reanalysis of the distribution of stellar radial velocities towards these clusters presents a significantly better agreement between stellar and gas kinematics than was previously thought.

We present a multiplicity study of all known protostars (94) in the Perseus molecular cloud from a Karl G. Jansky Very Large Array (VLA) survey at Ka-band (8 mm and 1 cm) and C-band (4 cm and 6.6 cm). The observed sample has a bolometric luminosity range between 0.1 L$_{\odot}$ and $\sim$33 L$_{\odot}$, with a median of 0.7 L$_{\odot}$. This multiplicity study is based on the Ka-band data, having a best resolution of $\sim$0.065" (15 AU) and separations out to $\sim$43" (10000 AU) can be probed. The overall multiplicity fraction (MF) is found to be of 0.40$\pm$0.06 and the companion star fraction (CSF) is 0.71$\pm$0.06. The MF and CSF of the Class 0 protostars are 0.57$\pm$0.09 and 1.2$\pm$0.2, and the MF and CSF of Class I protostars are both 0.23$\pm$0.08. The distribution of companion separations appears bi-modal, with a peak at $\sim$75 AU and another peak at $\sim$3000 AU. Turbulent fragmentation is likely the dominant mechanism on $>$1000 AU scales and disk fragmentation is likely to be the dominant mechanism on $<$200 AU scales. Toward three Class 0 sources we find companions separated by $<$30 AU. These systems have the smallest separations of currently known Class 0 protostellar binary systems. Moreover, these close systems are embedded within larger (50 AU to 400 AU) structures and may be candidates for ongoing disk fragmentation.

We present the first dust emission results toward a sample of seven protostellar disk candidates around Class 0 and I sources in the Perseus molecular cloud from the VLA Nascent Disk and Multiplicity (VANDAM) survey with ~0.05'' or 12 AU resolution. To examine the surface brightness profiles of these sources, we fit the Ka-band 8 mm dust-continuum data in the u,v-plane to a simple, parametrized model based on the Shakura-Sunyaev disk model. The candidate disks are well-fit by a model with a disk-shaped profile and have masses consistent with known Class 0 and I disks. The inner-disk surface densities of the VANDAM candidate disks have shallower density profiles compared to disks around more evolved Class II systems. The best-fit model radii of the seven early-result candidate disks are R_c > 10 AU; at 8 mm, the radii reflect lower limits on the disk size since dust continuum emission is tied to grain size and large grains radially drift inwards. These relatively large disks, if confirmed kinematically, are inconsistent with theoretical models where the disk size is limited by strong magnetic braking to < 10 AU at early times.

We present the results of the SDSS APOGEE INfrared Spectroscopy of Young Nebulous Clusters program (IN-SYNC) survey of the Orion A molecular cloud. This survey obtained high resolution near infrared (NIR) spectroscopy of about 2700 young pre-main sequence stars throughout the region, acquired across five distinct fields spanning 6deg field of view (FOV). With these spectra, we have measured accurate stellar parameters (T_eff, log g, v sin i) and extinctions, and placed the sources in the Hertzsprung-Russel Diagram (HRD). We have also extracted radial velocities for the kinematic characterization of the population. We compare our measurements with literature results for a sub-sample of targets in order to assess the performances and accuracy of the survey. Source extinction shows evidence for dust grains that are larger than those in the diffuse interstellar medium (ISM): we estimate an average R_V=5.5 in the region. Importantly, we find a clear correlation between HRD inferred ages and spectroscopic surface-gravity inferred ages. This clearly indicates a real spread of stellar radii at fixed temperature, and together with additional correlations with extinction and with disk presence, strongly suggests a real spread of ages large than a few Myr. Focussing on the young population around NGC1980 iota Ori, which has previously been suggested to be a separate, foreground, older cluster, we confirm its older (5Myr) age and low A_V, but considering that its radial velocity distribution is indistinguishable from the Orion A's population, we suggest that NGC1980 is part of Orion A's star formation activity. Based on their stellar parameters and kinematic properties, we identify 383 new candidate members of Orion A, most of which are diskless sources in areas of the region poorly studied by previous works.

We present deep radio continuum observations of the star-forming core of the Serpens South Infrared Dark Cloud with the Karl G. Jansky Very Large Array (VLA). Observations were conducted in two bands centered at 7.25 GHz (4.14 cm) and 4.75 GHz (6.31 cm) with an rms of 8.5 and 11.1 microJy/beam, respectively. We also use 2MASS, Spitzer and Herschel data to put our radio observations in the context of young stellar populations characterized by near and far infrared observations. Within a 5 arcmin x 5 arcmin region of interest around the central cluster, we detect roughly eighteen radio sources, seven of which we determine are protostellar in nature due to their radio spectral indices and their association with infrared sources. We find evidence for a previously undetected embedded Class 0 protostar and reaffirm Class 0 protostellar classifications determined by previous millimeter wavelength continuum studies. We use our infrared data to derive mid-infrared luminosities for three of our protostellar sources and find relative agreement between the known YSO radio luminosity vs bolometric luminosity correlation. Lastly, we marginally detect an additional six radio sources at the 2-3 sigma level that lie within two arcseconds of infrared YSO candidates, providing motivation for higher sensitivity studies to clarify the nature of these sources and further probe embedded and/or low luminosity YSOs in Serpens South.

We present multi-epoch, large-scale ($\sim$ 2000 arcmin${}^2$), fairly deep ($\sim$ 16 $\mu$Jy), high-resolution ($\sim$ 1") radio observations of the Perseus star-forming complex obtained with the Karl G. Jansky Very Large Array at frequencies of 4.5 GHz and 7.5 GHz. These observations were mainly focused on the clouds NGC 1333 and IC 348, although we also observed several fields in other parts of the Perseus complex. We detect a total of 206 sources, 42 of which are associated with young stellar objects (YSOs). The radio properties of about 60% of the YSOs are compatible with a non-thermal radio emission origin. Based on our sample, we find a fairly clear relation between the prevalence of non-thermal radio emission and evolutionary status of the YSOs. By comparing our results with previously reported X-ray observations, we show that YSOs in Perseus follow a G\"udel-Benz relation with $\kappa$ = 0.03 consistent with other regions of star formation. We argue that most of the sources detected in our observations but not associated with known YSOs are extragalactic, but provide a list of 20 unidentified radio sources whose radio properties are consistent with being YSO candidates. Finally we also detect 5 sources with extended emission features which can clearly be associated with radio galaxies.

We have surveyed 84 Class 0, Class I, and flat-spectrum protostars in mid-infrared [Si II], [Fe II] and [S I] line emission, and 11 of these in far-infrared [O I] emission. We use the results to derive their mass outflow rates. Thereby we observe a strong correlation of mass outflow rates with bolometric luminosity, and with the inferred mass accretion rates of the central objects, which continues through the Class 0 range the trend observed in Class II young stellar objects. Along this trend from large to small mass-flow rates, the different classes of young stellar objects lie in the sequence Class 0 -- Class I/flat-spectrum -- Class II, indicating that the trend is an evolutionary sequence in which mass outflow and accretion rates decrease together with increasing age, while maintaining rough proportionality. The survey results include two which are key tests of magnetocentrifugal outflow-acceleration mechanisms: the distribution of the outflow/accretion branching ratio b, and limits on the distribution of outflow speeds. Neither rule out any of the three leading outflow-acceleration, angular-momentum-ejection mechanisms, but they provide some evidence that disk winds and accretion-powered stellar winds (APSWs) operate in many protostars. An upper edge observed in the branching-ratio distribution is consistent with the upper bound of b = 0.6 found in models of APSWs, and a large fraction (0.31) of the sample have branching ratio sufficiently small that only disk winds, launched on scales as large as several AU, have been demonstrated to account for them.

We present continuum and molecular line observations at 230 GHz and 345 GHz from the Sub-millimeter Array (SMA) toward three protostars in the Perseus L1448N region. The data are from the large project "Mass Assembly of Stellar Systems and their Evolution with the SMA" (MASSES). Three dust continuum sources, Source B, Source NW, and Source A, are detected at both frequencies. These sources have corresponding emission peaks in C18O (J=2-1), 13CO (J=2-1), and HCO+ (J=4-3), and have offsets with N2D+ (J=3-2) peaks. High angular resolution data from a complimentary continuum survey with the Karl G. Jansky Very Large Array show that Source B is associated with three 8 mm continuum objects, Source NW with two, and Source A remains single. These results suggest that multiplicity in L1448N exists at different spatial scales from a few thousand AU to < 100 AU. Velocity gradients in each source obtained from two-dimensional fits to the SMA C18O emission are found to be perpendicular to within 20 degrees of the outflow directions as revealed by 12CO (J=2-1). We have observed that Sources B and NW with multiplicity have higher densities than Source A without multiplicity. This suggests that thermal Jeans fragmentation can be relevant in the fragmentation process. However, we have not observed a difference in the ratio between rotational and gravitational energy between sources with and without multiplicity. We also have not observed a trend between non-thermal velocity dispersions and the level of fragmentation. Our study has provided the first direct and comprehensive comparison between multiplicity and core properties in low-mass protostars, although based on small number statistics.

Magnetic fields can regulate disk formation, accretion and jet launching. Until recently, it has been difficult to obtain high resolution observations of the magnetic fields of the youngest protostars in the critical region near the protostar. The VANDAM survey is observing all known protostars in the Perseus Molecular Cloud. Here we present the polarization data of IRAS 4A. We find that with ~ 0.2'' (50 AU) resolution at {\lambda} = 8.1 and 10.3 mm, the inferred magnetic field is consistent with a circular morphology, in marked contrast with the hourglass morphology seen on larger scales. This morphology is consistent with frozen-in field lines that were dragged in by rotating material entering the infall region. The field morphology is reminiscent of rotating circumstellar material near the protostar. This is the first polarization detection of a protostar at these wavelengths. We conclude from our observations that the dust emission is optically thin with {\beta} ~ 1.3, suggesting that mm/cm-sized grains have grown and survived in the short lifetime of the protostar.

This paper discusses compelling science cases for a future long-baseline interferometer operating at millimeter and centimeter wavelengths, like the proposed Next Generation Vary Large Array (ngVLA). We report on the activities of the Cradle of Life science working group, which focused on the formation of low- and high-mass stars, the formation of planets and evolution of protoplanetary disks, the physical and compositional study of Solar System bodies, and the possible detection of radio signals from extraterrestrial civilizations. We propose 19 scientific projects based on the current specification of the ngVLA. Five of them are highlighted as possible Key Science Projects: (1) Resolving the density structure and dynamics of the youngest HII regions and high-mass protostellar jets, (2) Unveiling binary/multiple protostars at higher resolution, (3) Mapping planet formation regions in nearby disks on scales down to 1 AU, (4) Studying the formation of complex molecules, and (5) Deep atmospheric mapping of giant planets in the Solar System. For each of these projects, we discuss the scientific importance and feasibility. The results presented here should be considered as the beginning of a more in-depth analysis of the science enabled by such a facility, and are by no means complete or exhaustive.

Azimuthally asymmetric dust distributions observed with ALMA in transition disks have been interpreted as dust traps. We present VLA Ka band (34 GHz or 0.9 cm) and ALMA Cycle 2 Band 9 (680 GHz or 0.45 mm) observations at 0.2" resolution of the Oph IRS 48 disk, which suggest that larger particles could be more azimuthally concentrated than smaller dust grains, assuming an axisymmetric temperature field or optically thin 680 GHz emission. Fitting an intensity model to both data demonstrates that the azimuthal extent of the millimeter emission is 2.3 $\pm0.9$ times as wide as the centimeter emission, marginally consistent with the particle trapping mechanism under the above assumptions. The 34 GHz continuum image also reveals evidence for ionized gas emission from the star. Both the morphology and the spectral index variations are consistent with an increase of large particles in the center of the trap, but uncertainties remain due to the continuum optical depth at 680 GHz. Particle trapping has been proposed in planet formation models to allow dust particles to grow beyond millimeter sizes in the outer regions of protoplanetary disks. The new observations in the Oph IRS 48 disk provide support for the dust trapping mechanism for centimeter-sized grains, although additional data is required for definitive confirmation.

We present deep ($\sim 17~\mu$Jy) radio continuum observations of the Serpens molecular cloud, the Serpens south cluster, and the W40 region obtained using the Very Large Array in its A configuration. We detect a total of 146 sources, 29 of which are young stellar objects (YSOs), 2 are BV stars and 5 more are associated with phenomena related to YSOs. Based on their radio variability and spectral index, we propose that about 16 of the remaining 110 unclassified sources are also YSOs. For approximately 65% of the known YSOs detected here as radio sources, the emission is most likely non-thermal, and related to stellar coronal activity. As also recently observed in Ophiuchus, our sample of YSOs with X-ray counterparts lies below the fiducial G\"udel & Benz relation. Finally, we analyze the proper motions of 9 sources in the W40 region. This allows us to better constrain the membership of the radio sources in the region.

We present a CARMA 1.3 mm continuum survey toward 9 Class 0 protostars in the Perseus molecular cloud at $\sim$0.3$^{\prime\prime}$ (70 AU) resolution. This study approximately doubles the number of Class 0 protostars observed with spatial resolutions $<$ 100 AU at millimeter wavelengths, enabling the presence of protostellar disks and proto-binary systems to be probed. We detect flattened structures with radii $>$ 100 AU around 2 sources (L1448 IRS2 and Per-emb-14) and these sources may be strong disk candidates. Marginally-resolved structures with position angles within 30$^{\circ}$ of perpendicular to the outflow are found toward 3 protostars (L1448 IRS3C, IRAS 03282+3035, and L1448C) and are considered disk candidates. Two others (L1448 IRS3B and IRAS 03292+3039) have resolved structure, possibly indicative of massive inner envelopes or disks; L1448 IRS3B also has a companion separated by 0.9$^{\prime\prime}$ ($\sim$210 AU). IC348-MMS does not have well-resolved structure and the candidate first hydrostatic core L1451-MMS is marginally resolved on 1$^{\prime\prime}$ scales. The strong disk candidate sources were followed-up with C$^{18}$O ($J=2\rightarrow1$) observations, detecting velocity gradients consistent with rotation, but it is unclear if the rotation is Keplerian. We compare the observed visibility amplitudes to radiative transfer models, finding that visibility amplitude ratios suggest a compact component (possibly a disk) is necessary for 5 of 9 Class 0 sources; envelopes alone may explain the other 4 systems. We conclude that there is evidence for the formation of large disks in the Class 0 phase with a range of radii and masses dependent upon their initial formation conditions.

We present an expanded kinematic study of the young cluster NGC 2264 based upon optical radial velocities measured using multi-fiber echelle spectroscopy at the 6.5 meter MMT and Magellan telescopes. We report radial velocities for 695 stars, of which approximately 407 stars are confirmed or very likely members. Our results more than double the number of members with radial velocities from F{\H u}r{\'e}sz et al., resulting in a much better defined kinematic relationship between the stellar population and the associated molecular gas. In particular, we find that there is a significant subset of stars that are systematically blueshifted with respect to the molecular ($^{13}$CO) gas. The detection of Lithium absorption and/or infrared excesses in this blue-shifted population suggests that at least some of these stars are cluster members; we suggest some speculative scenarios to explain their kinematics. Our results also more clearly define the redshifted population of stars in the northern end of the cluster; we suggest that the stellar and gas kinematics of this region are the result of a bubble driven by the wind from O7 star S Mon. Our results emphasize the complexity of the spatial and kinematic structure of NGC 2264, important for eventually building up a comprehensive picture of cluster formation.

We report the dramatic mid-infrared brightening between 2004 and 2006 of HOPS 383, a deeply embedded protostar adjacent to NGC 1977 in Orion. By 2008, the source became a factor of 35 brighter at 24 microns with a brightness increase also apparent at 4.5 microns. The outburst is also detected in the submillimeter by comparing APEX/SABOCA to SCUBA data, and a scattered-light nebula appeared in NEWFIRM K_s imaging. The post-outburst spectral energy distribution indicates a Class 0 source with a dense envelope and a luminosity between 6 and 14 L_sun. Post-outburst time-series mid- and far-infrared photometry shows no long-term fading and variability at the 18% level between 2009 and 2012. HOPS 383 is the first outbursting Class 0 object discovered, pointing to the importance of episodic accretion at early stages in the star formation process. Its dramatic rise and lack of fading over a six-year period hint that it may be similar to FU Ori outbursts, although the luminosity appears to be significantly smaller than the canonical luminosities of such objects.

We present a multi-epoch radio study of the Taurus-Auriga star-forming complex made with the Karl G. Jansky Very Large Array at frequencies of 4.5 GHz and 7.5 GHz. We detect a total of 610 sources, 59 of which are related to young stellar objects and 18 to field stars. The properties of 56\% of the young stars are compatible with non-thermal radio emission. We also show that the radio emission of more evolved young stellar objects tends to be more non-thermal in origin and, in general, that their radio properties are compatible with those found in other star forming regions. By comparing our results with previously reported X-ray observations, we notice that young stellar objects in Taurus-Auriga follow a G\"{u}del-Benz relation with $\kappa$=0.03, as we previously suggested for other regions of star formation. In general, young stellar objects in Taurus-Auriga and in all the previous studied regions seem to follow this relation with a dispersion of $\sim1$ dex. Finally, we propose that most of the remaining sources are related with extragalactic objects but provide a list of 46 unidentified radio sources whose radio properties are compatible with a YSO nature.

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.

For the Class 0 protostar, L1527, we compare 131 polarization vectors from SCUPOL/JCMT, SHARP/CSO and TADPOL/CARMA observations with the corresponding model polarization vectors of four ideal-MHD, non-turbulent, cloud core collapse models. These four models differ by their initial magnetic fields before collapse; two initially have aligned fields (strong and weak) and two initially have orthogonal fields (strong and weak) with respect to the rotation axis of the L1527 core. Only the initial weak orthogonal field model produces the observed circumstellar disk within L1527. This is a characteristic of nearly all ideal-MHD, non-turbulent, core collapse models. In this paper we test whether this weak orthogonal model also has the best agreement between its magnetic field structure and that inferred from the polarimetry observations of L1527. We found that this is not the case; based on the polarimetry observations the most favored model of the four is the weak aligned model. However, this model does not produce a circumstellar disk, so our result implies that a non-turbulent, ideal-MHD global collapse model probably does not represent the core collapse that has occurred in L1527. Our study also illustrates the importance of using polarization vectors covering a large area of a cloud core to determine the initial magnetic field orientation before collapse; the inner core magnetic field structure can be highly altered by a collapse and so measurements from this region alone can give unreliable estimates of the initial field configuration before collapse.

We present CARMA 2.9 mm dust continuum emission observations of a sample of 14 Herschel-detected Class 0 protostars in the Orion A and B molecular clouds, drawn from the PACS Bright Red Sources (PBRS) sample (Stutz et al.). These objects are characterized by very red 24 \micron\ to 70 \micron\ colors and prominent submillimeter emission, suggesting that they are very young Class 0 protostars embedded in dense envelopes. We detect all of the PBRS in 2.9 mm continuum emission and emission from 4 protostars and 1 starless core in the fields toward the PBRS; we also report 1 new PBRS source. The ratio of 2.9 mm luminosity to bolometric luminosity is higher by a factor of $\sim$5 on average, compared to other well-studied protostars in the Perseus and Ophiuchus clouds. The 2.9 mm visibility amplitudes for 6 of the 14 PBRS are very flat as a function of uv-distance, with more than 50\% of the source emission arising from radii $<$ 1500 AU. These flat visibility amplitudes are most consistent with spherically symmetric envelope density profiles with $\rho$~$\propto$~R$^{-2.5}$. Alternatively, there could be a massive unresolved structure like a disk or a high-density inner envelope departing from a smooth power-law. The large amount of mass on scales $<$ 1500 AU (implying high average central densities) leads us to suggest that that the PBRS with flat visibility amplitude profiles are the youngest PBRS and may be undergoing a brief phase of high mass infall/accretion and are possibly among the youngest Class 0 protostars. The PBRS with more rapidly declining visibility amplitudes still have large envelope masses, but could be slightly more evolved.

We are conducting a Jansky VLA Ka-band (8 mm and 1 cm) and C-band (4 cm and 6.4 cm) survey of all known protostars in the Perseus Molecular Cloud, providing resolution down to $\sim$0.06'' and $\sim$0.35" in Ka-band and C-band, respectively. Here we present first results from this survey that enable us to examine the source NGC 1333 IRAS2A in unprecedented detail and resolve it into a proto-binary system separated by 0.621"$\pm$0.006" ($\sim$143 AU) at 8 mm, 1 cm, and 4 cm. These 2 sources (IRAS2A VLA1 and VLA2) are likely driving the two orthogonal outflows known to originate from IRAS2A. The brighter source IRAS2A VLA1 is extended perpendicular to its outflow in the VLA data, with a deconvolved size of 0.055" ($\sim$13 AU), possibly tracing a protostellar disk. The recently reported candidate companions (IRAS2A MM2 and MM3) are not detected in either our VLA data, CARMA 1.3 mm data, or SMA 850 $\mu$m data. SMA CO ($J=3\rightarrow2$), CARMA CO ($J=2\rightarrow1$), and lower resolution CARMA CO ($J=1\rightarrow0$) observations are used to examine the outflow origins and the nature of the candidate companions to IRAS2A VLA1. The CO ($J=3\rightarrow2$) and ($J=2\rightarrow1$) data show that IRAS2A MM2 is coincident with a bright CO emission spot in the east-west outflow, and IRAS2A MM3 is within the north-south outflow. In contrast, IRAS2A VLA2 lies at the east-west outflow symmetry point. We propose that IRAS2A VLA2 is the driving source of the East-West outflow and a true companion to IRAS2A VLA1, whereas IRAS2A MM2 and MM3 may not be protostellar.

We present observations of N2H+(1-0), HCO+(1-0), and HCN(1-0) toward the Serpens Main molecular cloud from the CARMA Large Area Star Formation Survey (CLASSy). We mapped 150 square arcminutes of Serpens Main with an angular resolution of 7 arcsecs. The gas emission is concentrated in two subclusters (the NW and SE subclusters). The SE subcluster has more prominent filamentary structures and more complicated kinematics compared to the NW subcluster. The majority of gas in the two subclusters has subsonic to sonic velocity dispersions. We applied a dendrogram technique with N2H+(1-0) to study the gas structures; the SE subcluster has a higher degree of hierarchy than the NW subcluster. Combining the dendrogram and line fitting analyses reveals two distinct relations: a flat relation between nonthermal velocity dispersion and size, and a positive correlation between variation in velocity centroids and size. The two relations imply a characteristic depth of 0.15 pc for the cloud. Furthermore, we have identified six filaments in the SE subcluster. These filaments have lengths of 0.2 pc and widths of 0.03 pc, which is smaller than a characteristic width of 0.1 pc suggested by Herschel observations. The filaments can be classified into two types based on their properties. The first type, located in the northeast of the SE subcluster, has larger velocity gradients, smaller masses, and nearly critical mass-per-unit-length ratios. The other type, located in the southwest of the SE subcluster, has the opposite properties. Several YSOs are formed along two filaments which have supercritical mass per unit length ratios, while filaments with nearly critical mass-per-unit-length ratios are not associated with YSOs, suggesting that stars are formed on gravitationally unstable filaments.

We present details of the CARMA Large Area Star Formation Survey (CLASSy), while focusing on observations of Barnard 1. CLASSy is a CARMA Key Project that spectrally imaged N2H+, HCO+, and HCN (J=1-0 transitions) across over 800 square arcminutes of the Perseus and Serpens Molecular Clouds. The observations have angular resolution near 7" and spectral resolution near 0.16 km/s. We imaged ~150 square arcminutes of Barnard 1, focusing on the main core, and the B1 Ridge and clumps to its southwest. N2H+ shows the strongest emission, with morphology similar to cool dust in the region, while HCO+ and HCN trace several molecular outflows from a collection of protostars in the main core. We identify a range of kinematic complexity, with N2H+ velocity dispersions ranging from ~0.05-0.50 km/s across the field. Simultaneous continuum mapping at 3 mm reveals six compact object detections, three of which are new detections. A new non-binary dendrogram algorithm is used to analyze dense gas structures in the N2H+ position-position-velocity (PPV) cube. The projected sizes of dendrogram-identified structures range from about 0.01-0.34 pc. Size-linewidth relations using those structures show that non-thermal line-of-sight velocity dispersion varies weakly with projected size, while rms variation in the centroid velocity rises steeply with projected size. Comparing these relations, we propose that all dense gas structures in Barnard 1 have comparable depths into the sky, around 0.1-0.2 pc; this suggests that over-dense, parsec-scale regions within molecular clouds are better described as flattened structures rather than spherical collections of gas. Science-ready PPV cubes for Barnard 1 molecular emission are available for download.

We present the N2H+(J=1-0) map of the Serpens South molecular cloud obtained as part of the CARMA Large Area Star Formation Survey (CLASSy). The observations cover 250 square arcminutes and fully sample structures from 3000 AU to 3 pc with a velocity resolution of 0.16 km/s, and they can be used to constrain the origin and evolution of molecular cloud filaments. The spatial distribution of the N2H+ emission is characterized by long filaments that resemble those observed in the dust continuum emission by Herschel. However, the gas filaments are typically narrower such that, in some cases, two or three quasi-parallel N2H+ filaments comprise a single observed dust continuum filament. The difference between the dust and gas filament widths casts doubt on Herschel ability to resolve the Serpens South filaments. Some molecular filaments show velocity gradients along their major axis, and two are characterized by a steep velocity gradient in the direction perpendicular to the filament axis. The observed velocity gradient along one of these filaments was previously postulated as evidence for mass infall toward the central cluster, but these kind of gradients can be interpreted as projection of large-scale turbulence.

We present results from a high-sensitivity (60 $\mu$Jy), large-scale (2.26 square degree) survey obtained with the Karl G. Jansky Very Large Array as part of the Gould's Belt Survey program. We detected 374 and 354 sources at 4.5 and 7.5 GHz, respectively. Of these, 148 are associated with previously known Young Stellar Objects (YSOs). Another 86 sources previously unclassified at either optical or infrared wavelengths exhibit radio properties that are consistent with those of young stars. The overall properties of our sources at radio wavelengths such as their variability and radio to X-ray luminosity relation are consistent with previous results from the Gould's Belt Survey. Our detections provide target lists for followup VLBA radio observations to determine their distances as YSOs are located in regions of high nebulosity and extinction, making it difficult to measure optical parallaxes.

We present {\lambda}1.3 mm CARMA observations of dust polarization toward 30 star-forming cores and 8 star-forming regions from the TADPOL survey. We show maps of all sources, and compare the ~2.5" resolution TADPOL maps with ~20" resolution polarization maps from single-dish submillimeter telescopes. Here we do not attempt to interpret the detailed B-field morphology of each object. Rather, we use average B-field orientations to derive conclusions in a statistical sense from the ensemble of sources, bearing in mind that these average orientations can be quite uncertain. We discuss three main findings: (1) A subset of the sources have consistent magnetic field (B-field) orientations between large (~20") and small (~2.5") scales. Those same sources also tend to have higher fractional polarizations than the sources with inconsistent large-to-small-scale fields. We interpret this to mean that in at least some cases B-fields play a role in regulating the infall of material all the way down to the ~1000 AU scales of protostellar envelopes. (2) Outflows appear to be randomly aligned with B-fields; although, in sources with low polarization fractions there is a hint that outflows are preferentially perpendicular to small-scale B-fields, which suggests that in these sources the fields have been wrapped up by envelope rotation. (3) Finally, even at ~2.5" resolution we see the so-called "polarization hole" effect, where the fractional polarization drops significantly near the total intensity peak. All data are publicly available in the electronic edition of this article.

We use mid-infrared to submillimeter data from the Spitzer, Herschel, and APEX telescopes to study the bright sub-mm source OMC-2 FIR 4. We find a point source at 8, 24, and 70 $\mu$m, and a compact, but extended source at 160, 350, and 870 $\mu$m. The peak of the emission from 8 to 70 $\mu$m, attributed to the protostar associated with FIR 4, is displaced relative to the peak of the extended emission; the latter represents the large molecular core the protostar is embedded within. We determine that the protostar has a bolometric luminosity of 37 Lsun, although including more extended emission surrounding the point source raises this value to 86 Lsun. Radiative transfer models of the protostellar system fit the observed SED well and yield a total luminosity of most likely less than 100 Lsun. Our models suggest that the bolometric luminosity of the protostar could be just 12-14 Lsun, while the luminosity of the colder (~ 20 K) extended core could be around 100 Lsun, with a mass of about 27 Msun. Our derived luminosities for the protostar OMC-2 FIR 4 are in direct contradiction with previous claims of a total luminosity of 1000 Lsun (Crimier et al 2009). Furthermore, we find evidence from far-infrared molecular spectra (Kama et al. 2013, Manoj et al. 2013) and 3.6 cm emission (Reipurth et al 1999) that FIR 4 drives an outflow. The final stellar mass the protostar will ultimately achieve is uncertain due to its association with the large reservoir of mass found in the cold core.

Stars form from the gravitational collapse of dense molecular cloud cores. In the protostellar phase, mass accretes from the core onto a protostar, likely through an accretion disk, and it is during this phase that the initial masses of stars and the initial conditions for planet formation are set. Over the past decade, new observational capabilities provided by the Spitzer Space Telescope and Herschel Space Observatory have enabled wide-field surveys of entire star-forming clouds with unprecedented sensitivity, resolution, and infrared wavelength coverage. We review resulting advances in the field, focusing both on the observations themselves and the constraints they place on theoretical models of star formation and protostellar evolution. We also emphasize open questions and outline new directions needed to further advance the field.

Edge-on protostars are valuable for understanding the disk and envelope properties of embedded young stellar objects, since the disk, envelope, and envelope cavities are all distinctly visible in resolved images and well constrained in modeling. Comparing 2MASS, WISE, Spitzer, Herschel, APEX, and IRAM photometry and limits from 1.25 to 1200 microns, Spitzer spectroscopy from 5 to 40 microns, and high-resolution Hubble imaging at 1.60 and 2.05 microns to radiative transfer modeling, we determine envelope and disk properties for the Class I protostar HOPS 136, an edge-on source in Orion's Lynds 1641 region. The source has a bolometric luminosity of 0.8 L_sun, a bolometric temperature of 170 K, and a ratio of submillimeter to bolometric luminosity of 0.8%. Via modeling, we find a total luminosity of 4.7 L_sun (larger than the observed luminosity due to extinction by the disk), an envelope mass of 0.06 M_sun, and a disk radius and mass of 450 AU and 0.002 M_sun. The stellar mass is highly uncertain but is estimated to fall between 0.4 and 0.5 M_sun. To reproduce the flux and wavelength of the near-infrared scattered-light peak in the spectral energy distribution, we require 5.4 x 10^-5 M_sun of gas and dust in each cavity. The disk has a large radius and a mass typical of more evolved T Tauri stars in spite of the significant remaining envelope. HOPS 136 appears to be a key link between the protostellar and optically revealed stages of star formation.

We present observations of three Class 0/I protostars (L1157-mm, CB230 IRS1, and L1165-SMM1) using the Karl G. Jansky Very Large Array (VLA) and observations of two (L1165-SMM1 and CB230 IRS1) with the Combined Array for Research in Millimeter-wave Astronomy (CARMA). The VLA observations were taken at wavelengths of $\lambda = 7.3$ mm, 1.4 cm, 3.3 cm, 4.0 cm, and 6.5 cm with a best resolution of $\sim$0\farcs06 (18 AU) at 7.3 mm. The L1165-SMM1 CARMA observations were taken at $\lambda = 1.3$ mm with a best resolution of $\sim0\farcs3$ (100 AU), and the CB230 IRS1 observations were taken at $\lambda = 3.4$ mm with a best resolution of $\sim$3\arcsec\ (900 AU). We find that L1165-SMM1 and CB230 IRS1 have probable binary companions at separations of $\sim$0\farcs3 (100 AU) from detections of secondary peaks at multiple wavelengths. The position angles of these companions are nearly orthogonal to the direction of the observed bipolar outflows, consistent with the expected protostellar disk orientations. We suggest that these companions may have formed from disk fragmentation; turbulent fragmentation would not preferentially arrange the binary companions to be orthogonal to the outflow direction. For L1165-SMM1, both the 7.3 mm and 1.3 mm emission show evidence of a large (R $>$ 100 AU) disk. For the L1165-SMM1 primary protostar and the CB230 IRS1 secondary protostar, the 7.3 mm emission is resolved into structures consistent with $\sim20$ AU radius disks. For the other protostars, including L1157-mm, the emission is unresolved, suggesting disks with radii $< 20$ AU.

We present large-scale ($\sim$ 2000 square arcminutes), deep ($\sim$ 20 $\mu$Jy), high-resolution ($\sim$ 1$''$) radio observations of the Ophiuchus star-forming complex obtained with the Karl G.\ Jansky Very Large Array at $\lambda$ = 4 and 6 cm. In total, 189 sources were detected, 56 of them associated with known young stellar sources, and 4 with known extragalactic objects; the other 129 remain unclassified, but most of them are most probably background quasars. The vast majority of the young stars detected at radio wavelengths have spectral types K or M, although we also detect 4 objects of A/F/B types and 2 brown dwarf candidates. At least half of these young stars are non-thermal (gyrosynchrotron) sources, with active coronas characterized by high levels of variability, negative spectral indices, and (in some cases) significant circular polarization. As expected, there is a clear tendency for the fraction of non-thermal sources to increase from the younger (Class 0/I or flat spectrum) to the more evolved (Class III or weak line T Tauri) stars. The young stars detected both in X-rays and at radio wavelengths broadly follow a G\"udel-Benz relation, but with a different normalization than the most radio-active types of stars. Finally, we detect a $\sim$ 70 mJy compact extragalactic source near the center of the Ophiuchus core, which should be used as gain calibrator for any future radio observations of this region.

We present high-resolution sub/millimeter interferometric imaging of the Class 0 protostar L1527 IRS (IRAS 04368+2557) at 870 micron and 3.4 mm from the Submillimeter Array (SMA) and Combined Array for Research in Millimeter Astronomy (CARMA). We detect the signature of an edge-on disk surrounding the protostar with an observed diameter of 180 AU in the sub/millimeter images. The mass of the disk is estimated to be 0.007 M_sun, assuming optically thin, isothermal dust emission. The millimeter spectral index is observed to be quite shallow at all the spatial scales probed; alpha ~ 2, implying a dust opacity spectral index beta ~ 0. We model the emission from the disk and surrounding envelope using Monte Carlo radiative transfer codes, simultaneously fitting the sub/millimeter visibility amplitudes, sub/millimeter images, resolved L\arcmin\ image, spectral energy distribution, and mid-infrared spectrum. The best fitting model has a disk radius of R = 125 AU, is highly flared (H ~ R^1.3), has a radial density profile rho ~ R^-2.5, and has a mass of 0.0075 M_sun. The scale height at 100 AU is 48 AU, about a factor of two greater than vertical hydrostatic equilibrium. The resolved millimeter observations indicate that disks may grow rapidly throughout the Class 0 phase. The mass and radius of the young disk around L1527 is comparable to disks around pre-main sequence stars; however, the disk is considerably more vertically extended, possibly due to a combination of lower protostellar mass, infall onto the disk upper layers, and little settling of ~1 micron-sized dust grains.

We present results of 1.3 mm dust polarization observations toward 16 nearby, low-mass protostars, mapped with ~2.5" resolution at CARMA. The results show that magnetic fields in protostellar cores on scales of ~1000 AU are not tightly aligned with outflows from the protostars. Rather, the data are consistent with scenarios where outflows and magnetic fields are preferentially misaligned (perpendicular), or where they are randomly aligned. If one assumes that outflows emerge along the rotation axes of circumstellar disks, and that the outflows have not disrupted the fields in the surrounding material, then our results imply that the disks are not aligned with the fields in the cores from which they formed.

We perform a census of the reddest, and potentially youngest, protostars in the Orion molecular clouds using data obtained with the PACS instrument onboard the Herschel Space Observatory and the LABOCA and SABOCA instruments on APEX as part of the Herschel Orion Protostar Survey (HOPS). A total of 55 new protostar candidates are detected at 70 um and 160 um that are either too faint (m24 > 7 mag) to be reliably classified as protostars or undetected in the Spitzer/MIPS 24 um band. We find that the 11 reddest protostar candidates with log (lambda F_lambda 70) / (lambda F_lambda 24) > 1.65 are free of contamination and can thus be reliably explained as protostars. The remaining 44 sources have less extreme 70/24 colors, fainter 70 um fluxes, and higher levels of contamination. Taking the previously known sample of Spitzer protostars and the new sample together, we find 18 sources that have log (lambda F_lambda 70) / (lambda F_lambda 24) > 1.65; we name these sources "PACS Bright Red sources", or PBRs. Our analysis reveals that the PBRs sample is composed of Class 0 like sources characterized by very red SEDs (T_bol < 45 K) and large values of sub-millimeter fluxes (L_smm/L_bol > 0.6%). Modified black-body fits to the SEDs provide lower limits to the envelope masses of 0.2 M_sun to 2 M_sun and luminosities of 0.7 L_sun to 10 L_sun. Based on these properties, and a comparison of the SEDs with radiative transfer models of protostars, we conclude that the PBRs are most likely extreme Class 0 objects distinguished by higher than typical envelope densities and hence, high mass infall rates.

The Wide-field Infrared Survey Explorer has revealed a T8.5 brown dwarf (WISE J111838.70+312537.9) that exhibits common proper motion with a solar-neighborhood (8 pc) quadruple star system - Xi Ursae Majoris. The angular separation is 8.5 arc-min, and the projected physical separation is about 4000 AU. The sub-solar metallicity and low chromospheric activity of Xi UMa A argue that the system has an age of at least 2 Gyr. The infrared luminosity and color of the brown dwarf suggests the mass of this companion ranges between 14 and 38 Jupiter masses for system ages of 2 and 8 Gyr respectively.

We present a study on the spatial distribution of N2D+ and N2H+ in thirteen protostellar systems. Eight of thirteen objects observed with the IRAM 30m telescope show relative offsets between the peak N2D+ (J=2-1) and N2H+ (J=1-0) emission. We highlight the case of L1157 using interferometric observations from the Submillimeter Array and Plateau de Bure Interferometer of the N2D+ (J=3-2) and N2H+ (J=1-0) transitions respectively. Depletion of N2D+ in L1157 is clearly observed inside a radius of ~2000 AU (7") and the N2H+ emission is resolved into two peaks at radii of ~1000 AU (3.5"), inside the depletion region of N2D+. Chemical models predict a depletion zone in N2H+ and N2D+ due to destruction of H2D+ at T ~ 20 K and the evaporation of CO off dust grains at the same temperature. However, the abundance offsets of 1000 AU between the two species are not reproduced by chemical models, including a model that follows the infall of the protostellar envelope. The average abundance ratios of N2D+ to N2H+ have been shown to decrease as protostars evolve by Emprechtinger et al., but this is the first time depletion zones of N2D+ have been spatially resolved. We suggest that the difference in depletion zone radii for N2H+ and N2D+ is caused by either the CO evaporation temperature being above 20 K or an H2 ortho-to-para ratio gradient in the inner envelope.

Surveys with the Spitzer and Herschel space observatories are now enabling the discovery and characterization of large samples of protostars in nearby molecular clouds, providing the observational basis for a detailed understanding of star formation in diverse environments. We are pursuing this goal with the Herschel Orion Protostar Survey (HOPS), which targets 328 Spitzer-identified protostars in the Orion molecular clouds, the largest star-forming region in the nearest 500 pc. The sample encompasses all phases of protostellar evolution and a wide range of formation environments, from dense clusters to relative isolation. With a grid of radiative transfer models, we fit the 1-870 micron spectral energy distributions (SEDs) of the protostars to estimate their envelope densities, cavity opening angles, inclinations, and total luminosities. After correcting the bolometric luminosities and temperatures of the sources for foreground extinction and inclination, we find a spread of several orders of magnitude in luminosity at all evolutionary states, a constant median luminosity over the more evolved stages, and a possible deficit of high-inclination, rapidly infalling envelopes among the Spitzer-identified sample. We have detected over 100 new sources in the Herschel images; some of them may fill this deficit. We also report results from modeling the pre- and post-outburst 1-870 micron SEDs of V2775 Ori (HOPS 223), a known FU Orionis outburster in the sample. It is the least luminous FU Ori star with a protostellar envelope.

We extend our previous study of the stellar population of L1641, the lower-density star-forming region of the Orion A cloud south of the dense Orion Nebula Cluster (ONC), with the goal of testing whether there is a statistically significant deficiency of high-mass stars in low-density regions. Previously, we compared the observed ratio of low-mass stars to high-mass stars with theoretical models of the stellar initial mass function (IMF) to infer a deficiency of the highest-mass stars in L1641. We expand our population study to identify the intermediate mass (late B to G) L1641 members in an attempt to make a more direct comparison with the mass function of the nearby ONC. The spectral type distribution and the K-band luminosity function of L1641 are similar to those of the ONC (Hillenbrand 1997; Muench et al. 2002), but problems of incompleteness and contamination prevent us from making a detailed test for differences. We limit our analysis to statistical tests of the ratio of high-mass to low-mass stars, which indicate a probability of only 3% that the ONC and the southern region of L1641 were drawn from the same population, supporting the hypothesis that the upper mass end of the IMF is dependent on environmental density.

In their earliest stages, protostars accrete mass from their surrounding envelopes through circumstellar disks. Until now, the smallest observed protostar/envelope mass ratio was ~2.1. The protostar L1527 IRS is thought to be in the earliest stages of star formation. Its envelope contains ~1 solar mass of material within a ~0.05 pc radius, and earlier observations suggested the presence of an edge-on disk. Here we report observations of dust continuum emission and 13CO (J=2-1) line emission from the disk around L1527, from which we determine a protostellar mass of M = 0.19 +/- 0.04 solar masses and a protostar/envelope mass ratio of ~0.2. We conclude that most of the luminosity is generated through the accretion process, with an accretion rate of ~6.6 x 10^-7 solar masses per year. If it has been accreting at that rate through much of its life, its age is ~300,000 yr, though theory suggests larger accretion rates earlier, so it may be younger. The presence of a rotationally--supported disk is confirmed and significantly more mass may be added to its planet-forming region as well as the protostar itself.

We present far-IR (57-196 mu) spectra of 21 protostars in the Orion molecular clouds, obtained with the Photodetector Array Camera and Spectrometer (PACS) onboard the Herschel Space observatory, as part of the Herschel Orion Protostar Survey (HOPS) program. We analyzed the CO emission lines (J_up = 14-46) in the PACS spectra, extracted within a projected distance of <= 2000 AU centered on the protostar. The total luminosity of the CO lines observed with PACS (L(CO)) is found to increase with increasing L_bol. The CO rotational temperature implied by the line ratios increases with J, and at least 3-4 rotational temperature components are required to fit the observed rotational diagram. The rotational temperature components are remarkably invariant between protostars and show no dependence on L_bol, T_bol or envelope density, implying that if the emitting gas is in LTE, the CO emission must arise in multiple temperature components that remain independent of L_bol over two orders of magnitudes. The observed CO emission can also be modeled as arising from a single temperature gas component or from a medium with a power-law temperature distribution; both of these require sub-thermally excited molecular gas at low densities (n(H_2) <= 10^6 cm^-3) and high temperatures (T >= 2000 K). Our results suggest that the contribution from PDRs along the envelope cavity walls is unlikely to be the dominant component of the CO emission observed with PACS. Instead, the "universality" of the rotational temperatures and the observed correlation between L(CO) and L_bol can most easily be explained if the observed CO emission originates in shock-heated, hot (T >= 2000 K), sub-thermally excited (n(H_2) <= 10^6 cm^-3) molecular gas. Post-shock gas at these densities is more likely to be found within the outflow cavities along the molecular outflow or along the cavity walls at radii >= several 100-1000 AU.

Filamentary structures are ubiquitous from large-scale molecular clouds (few parsecs) to small-scale circumstellar envelopes around Class 0 sources (~1000 AU to ~0.1 pc). In particular, recent observations with the Herschel Space Observatory emphasize the importance of large-scale filaments (few parsecs) and star formation. The small-scale flattened envelopes around Class 0 sources are reminiscent of the large-scale filaments. We propose an observationally derived scenario for filamentary star formation that describes the evolution of filaments as part of the process for formation of cores and circumstellar envelopes. If such a scenario is correct, small-scale filamentary structures (0.1 pc in length) with higher densities embedded in starless cores should exist, although to date almost all the interferometers have failed to observe such structures. We perform synthetic observations of filaments at the prestellar stage by modeling the known Class 0 flattened envelope in L1157 using both the Combined Array for Research in Millimeter-wave Astronomy (CARMA) and the Atacama Large Millimeter/Submillimeter Array (ALMA). We show that with reasonable estimates for the column density through the flattened envelope, the CARMA D-array at 3mm wavelengths is not able to detect such filamentary structure, so previous studies would not have detected them. However, the substructures may be detected with CARMA D+E array at 3 mm and CARMA E array at 1 mm as a result of more appropriate resolution and sensitivity. ALMA is also capable of detecting the substructures and showing the structures in detail compared to the CARMA results with its unprecedented sensitivity. Such detection will confirm the new proposed paradigm of non-spherical star formation.

We present dual-wavelength observations and modeling of the nearly edge-on Class 0 young stellar object L1157-mm. Using the Combined Array for Research in Millimeter-wave Astronomy, a nearly spherical structure is seen from the circumstellar envelope at the size scale of 10^2 to 10^3 AU in both 1 mm and 3 mm dust emission. Radiative transfer modeling is performed to compare data with theoretical envelope models, including a power-law envelope model and the Terebey-Shu-Cassen model. Bayesian inference is applied for parameter estimation and information criteria is used for model selection. The results prefer the power-law envelope model against the Terebey-Shu-Cassen model. In particular, for the power-law envelope model, a steep density profile with an index of ~2 is inferred. Moreover, the dust opacity spectral index (beta) is estimated to be ~0.9, implying that grain growth has started at L1157-mm. Also, the unresolved disk component is constrained to be < 40 AU in radius and < 4-25 M_Jup in mass. However, the estimate of the embedded disk component relies on the assumed envelope model.

Individual outbursting young stars are important laboratories for studying the physics of episodic accretion and the extent to which this phenomenon can explain the luminosity distribution of protostars. We present new and archival data for V2775 Ori (HOPS 223), a protostar in the L 1641 region of the Orion molecular clouds that was discovered by Caratti o Garatti et al. (2011) to have recently undergone an order-of-magnitude increase in luminosity. Our near-infrared spectra of the source have strong blueshifted He I 10830 absorption, strong H2O and CO absorption, and no H I emission, all typical of FU Orionis sources. With data from IRTF, 2MASS, HST, Spitzer, WISE, Herschel, and APEX that span from 1 to 70 microns pre-outburst and from 1 to 870 microns post-outburst, we estimate that the outburst began between 2005 April and 2007 March. We also model the pre- and post-outburst spectral energy distributions of the source, finding it to be in the late stages of accreting its envelope with a disk-to-star accretion rate that increased from about 2x10^-6 M_sun/yr to about 10^-5 M_sun/yr during the outburst. The post-outburst luminosity at the epoch of the FU Orionis-like near-IR spectra is 28 L_sun, making V2775 Ori the least luminous documented FU Orionis outburster with a protostellar envelope. The existence of low-luminosity outbursts supports the notion that a range of episiodic accretion phenomena can partially explain the observed spread in protostellar luminosities.

We present results from an optical photometric and spectroscopic survey of the young stellar population in L1641, the low-density star-forming region of the Orion A cloud south of the Orion Nebula Cluster (ONC). Our goal is to determine whether L1641 has a large enough low-mass population to make the known lack of high-mass stars a statistically-significant demonstration of environmental dependence of the upper mass stellar initial mass function (IMF). Our spectroscopic sample consists of IR-excess objects selected from the Spitzer/IRAC survey and non-excess objects selected from optical photometry. We have spectral confirmation of 864 members, with another 98 probable members; of the confirmed members, 406 have infrared excesses and 458 do not. Assuming the same ratio of stars with and without IR excesses in the highly-extincted regions, L1641 may contain as many as ~1600 stars down to ~0.1 solar mass, comparable within a factor of two to the the ONC. Compared to the standard models of the IMF, L1641 is deficient in O and early B stars to a 3-4 sigma significance level, assuming that we know of all the massive stars in L1641. With a forthcoming survey of the intermediate-mass stars, we will be in a better position to make a direct comparison with the neighboring, dense ONC, which should yield a stronger test of the dependence of the high-mass end of the stellar initial mass function upon environment.

We present an interferometric kinematic study of morphologically complex protostellar envelopes based on observations of the dense gas tracers N2H+ and NH3. The strong asymmetric nature of most envelopes in our sample leads us to question the common interpretation of velocity gradients as rotation, given the possibility of projection effects in the observed velocities. Several "idealized" sources with well-ordered velocity fields and envelope structures are now analyzed in more detail. We compare the interferometric data to position-velocity diagrams of kinematic models for spherical rotating collapse and filamentary rotating collapse. For this purpose, we developed a filamentary parametrization of the rotating collapse model to explore the effects of geometric projection on the observed velocity structures. We find that most envelopes in our sample have PV structures that can be reproduced by an infalling filamentary envelope projected at different angles within the plane of the sky. The infalling filament produces velocity shifts across the envelope that can mimic rotation, especially when viewed at single-dish resolutions and the axisymmetric rotating collapse model does not uniquely describe any dataset. Furthermore, if the velocities are assumed to reflect rotation, then the inferred centrifugal radii are quite large in most cases, indicating significant fragmentation potential or more likely another component to the line-center velocity. We conclude that ordered velocity gradients cannot be interpreted as rotation alone when envelopes are non-axisymmetric and that projected infall velocities likely dominate the velocity field on scales larger than 1000 AU.

We present a study of dense molecular gas kinematics in seventeen nearby protostellar systems using single-dish and interferometric molecular line observations. The non-axisymmetric envelopes around a sample of Class 0/I protostars were mapped in the N2H+ (J=1-0) tracer with the IRAM 30m, CARMA and PdBI as well as NH3 (1,1) with the VLA. The molecular line emission is used to construct line-center velocity and linewidth maps for all sources to examine the kinematic structure in the envelopes on spatial scales from 0.1 pc to ~1000 AU. The direction of the large-scale velocity gradients from single-dish mapping is within 45 degrees of normal to the outflow axis in more than half the sample. Furthermore, the velocity gradients are often quite substantial, the average being ~2.3 km\s\pc. The interferometric data often reveal small-scale velocity structure, departing from the more gradual large-scale velocity gradients. In some cases, this likely indicates accelerating infall and/or rotational spin-up in the inner envelope; the median velocity gradient from the interferometric data is ~10.7 km/s/pc. In two systems, we detect high-velocity HCO+ (J=1-0) emission inside the highest-velocity \nthp\ emission. This enables us to study the infall and rotation close to the disk and estimate the central object masses. The velocity fields observed on large and small-scales are more complex than would be expected from rotation alone, suggesting that complex envelope structure enables other dynamical processes (i.e. infall) to affect the velocity field.

We present the Spitzer Space Telescope Infrared Spectrograph spectrum of the Orion A protostar HOPS-68. The mid-infrared spectrum reveals crystalline substructure at 11.1, 16.1, 18.8, 23.6, 27.9, and 33.6 microns superimposed on the broad 9.7 and 18 micron amorphous silicate features; the substructure is well matched by the presence of the olivine end-member forsterite. Crystalline silicates are often observed as infrared emission features around the circumstellar disks of Herbig Ae/Be stars and T Tauri stars. However, this is the first unambiguous detection of crystalline silicate absorption in a cold, infalling, protostellar envelope. We estimate the crystalline mass fraction along the line-of-sight by first assuming that the crystalline silicates are located in a cold absorbing screen and secondly by utilizing radiative transfer models. The resulting crystalline mass fractions of 0.14 and 0.17, respectively, are significantly greater than the upper limit found in the interstellar medium (< 0.02-0.05). We propose that the amorphous silicates were annealed within the hot inner disk and/or envelope regions and subsequently transported outward into the envelope by entrainment in a protostellar outflow

We present high-resolution L'-band imaging of the inner scattered light structure of Class 0 protostar L1527 IRS (IRAS 04368+2557) taken with the Gemini North telescope. The central point-source like feature seen in Spitzer Space Telescope IRAC images is resolved in the Gemini image into a compact bipolar structure with a narrow dark lane in the center. Two scattered light lobes are extended ~1.8" (200 AU) perpendicular to the direction of the outflow and ~2.5" (350 AU) along the outflow axis; the narrow dark lane between the scattered light lobes is ~0.45" (60 AU) thick. The observations are consistent with our initial modeling of a bright inner cavity separated by a dark lane due to extinction along the line of sight of the central protostar by the disk (Tobin et al. 2008). The bright, compact scattered light might be due to complex inner structure generated by the outflow, as suggested in our first paper, or it may more likely be the upper layers of the disk forming from infalling matter.

Extinction maps at 8 micron from the Spitzer Space Telescope show that many Class 0 protostars exhibit complex, irregular, and on-axisymmetric structure within the densest regions of their dusty envelopes. Many of the systems have highly irregular and on-axisymmetric morphologies on scales $\sim$1000 AU, with a quarter of the sample exhibiting filamentary or flattened dense structures. Complex envelope structure is observed in regions spatially distinct from outflow cavities, and the densest structures often show no systematic alignment perpendicular to the cavities. We suggest that the observed envelope complexity is the result of collapse from protostellar cores with initially non-equilibrium structures. The striking non-axisymmetry in many envelopes could provide favorable conditions for the formation of binary systems. We then show that the kinematics around L1165 as probed with N2H+ are indicative of asymmetric infall; the velocity gradient is not perpendicular to the outflow.

We present the final results from a high sampling rate, multi-month, spectrophotometric reverberation mapping campaign undertaken to obtain either new or improved Hbeta reverberation lag measurements for several relatively low-luminosity AGNs. We have reliably measured thetime delay between variations in the continuum and Hbeta emission line in six local Seyfert 1 galaxies. These measurements are used to calculate the mass of the supermassive black hole at the center of each of these AGNs. We place our results in context to the most current calibration of the broad-line region (BLR) R-L relationship, where our results remove outliers and reduce the scatter at the low-luminosity end of this relationship. We also present velocity-resolved Hbeta time delay measurements for our complete sample, though the clearest velocity-resolved kinematic signatures have already been published.

Upcoming surveys with the Herschel Space Observatory will yield far-IR photometry of large samples of young stellar objects, which will require careful interpretation. We investigate the color and luminosity diagnostics based on Herschel broad-band filters to identify and discern the properties of low-mass protostars. We compute a grid of 2,016 protostars in various physical congurations, present the expected flux densities and flux density ratios for this grid of protostars, and compare Herschel observations of three protostars to the model results. These provide useful constraints on the range of colors and fluxes of protostar in the Herschel filters. We find that Herschel data alone is likely a useful diagnostic of the envelope properties of young stars

We present 70 and 160 micron Herschel science demonstration images of a field in the Orion A molecular cloud that contains the prototypical Herbig-Haro objects HH 1 and 2, obtained with the Photodetector Array Camera and Spectrometer (PACS). These observations demonstrate Herschel's unprecedented ability to study the rich population of protostars in the Orion molecular clouds at the wavelengths where they emit most of their luminosity. The four protostars previously identified by Spitzer 3.6-40 micron imaging and spectroscopy are detected in the 70 micron band, and three are clearly detected at 160 microns. We measure photometry of the protostars in the PACS bands and assemble their spectral energy distributions (SEDs) from 1 to 870 microns with these data, Spitzer spectra and photometry, 2MASS data, and APEX sub-mm data. The SEDs are fit to models generated with radiative transfer codes. From these fits we can constrain the fundamental properties of the protostars. We find luminosities in the range 12-84 L_sun and envelope densities spanning over two orders of magnitude. This implies that the four protostars have a wide range of envelope infall rates and evolutionary states: two have dense, infalling envelopes, while the other two have only residual envelopes. We also show the highly irregular and filamentary structure of the cold dust and gas surrounding the protostars as traced at 160 microns.

The NGC 1999 reflection nebula features a dark patch with a size of ~10,000 AU, which has been interpreted as a small, dense foreground globule and possible site of imminent star formation. We present Herschel PACS far-infrared 70 and 160mum maps, which reveal a flux deficit at the location of the globule. We estimate the globule mass needed to produce such an absorption feature to be a few tenths to a few Msun. Inspired by this Herschel observation, we obtained APEX LABOCA and SABOCA submillimeter continuum maps, and Magellan PANIC near-infrared images of the region. We do not detect a submillimer source at the location of the Herschel flux decrement; furthermore our observations place an upper limit on the mass of the globule of ~2.4x10^-2 Msun. Indeed, the submillimeter maps appear to show a flux depression as well. Furthermore, the near-infrared images detect faint background stars that are less affected by extinction inside the dark patch than in its surroundings. We suggest that the dark patch is in fact a hole or cavity in the material producing the NGC 1999 reflection nebula, excavated by protostellar jets from the V 380 Ori multiple system.

We use archived IRAC images from the Spitzer Space Telescope to show that many Class 0 protostars exhibit complex, irregular, and non-axisymmetric structure within their dusty envelopes. Our 8 $\mu$m extinction maps probe some of the densest regions in these protostellar envelopes. Many of the systems are observed to have highly irregular and non-axisymmetric morphologies on scales >= 1000 AU, with a quarter of the sample exhibiting filamentary or flattened dense structures. Complex envelope structure is observed in regions spatially distinct from outflow cavities, and the densest structures often show no systematic alignment perpendicular to the cavities. These results indicate that mass ejection is not responsible for much of the irregular morphologies we detect; rather, we suggest that the observed envelope complexity is mostly the result of collapse from protostellar cores with initially non-equilibrium structures. The striking non-axisymmetry in many envelopes could provide favorable conditions for the formation of binary systems. We also note that protostars in the sample appear to be formed preferentially near the edges of clouds or bends in filaments, suggesting formation by gravitational focusing.

We present observations of the Class 0 protostar L1157-mm using the Combined Array for Research in Millimeter-wave Astronomy (CARMA) in 3 mm dust continuum and N2H+ line emission. In the N2H+ line, we detect a large-scale envelope extended over a linear size of ~20,000AU flattened in the direction perpendicular to the outflow. This N2H+ feature coincides with the outer envelope seen in the 8 micron extinction by Looney et al. Meanwhile, the dust continuum traces the compact, nearly spherical structure of the inner envelope, where N2H+ becomes depleted. This highly flattened N2H+ envelope also shows dynamical signatures consistent with gravitational infall in the inner region, but a slow, solid-body rotation at large scales. This flattened structure is not a rotationally supported circumstellar disk; instead, it resembles a prestellar core both morphologically and kinematically, representing the early phase of a Class 0 system. In this paper, we construct a simple model to interpret both the dust continuum and N2H+ emission and suggest a possible dynamical scenario for the overall properties of the envelope.

A detailed analysis of the data from a high sampling rate, multi-month reverberation mapping campaign, undertaken primarily at MDM Observatory with supporting observations from telescopes around the world, reveals that the Hbeta emission region within the broad line regions (BLRs) of several nearby AGNs exhibit a variety of kinematic behaviors. While the primary goal of this campaign was to obtain either new or improved Hbeta reverberation lag measurements for several relatively low luminosity AGNs (presented in a separate work), we were also able to unambiguously reconstruct velocity-resolved reverberation signals from a subset of our targets. Through high cadence spectroscopic monitoring of the optical continuum and broad Hbeta emission line variations observed in the nuclear regions of NGC 3227, NGC 3516, and NGC 5548, we clearly see evidence for outflowing, infalling, and virialized BLR gas motions, respectively.

We present the first results from a high sampling rate, multi-month reverberation mapping campaign undertaken primarily at MDM Observatory with supporting observations from telescopes around the world. The primary goal of this campaign was to obtain either new or improved Hbeta reverberation lag measurements for several relatively low luminosity AGNs. We feature results for NGC 4051 here because, until now, this object has been a significant outlier from AGN scaling relationships, e.g., it was previously a ~2-3sigma outlier on the relationship between the broad-line region (BLR) radius and the optical continuum luminosity - the R_BLR-L relationship. Our new measurements of the lag time between variations in the continuum and Hbeta emission line made from spectroscopic monitoring of NGC 4051 lead to a measured BLR radius of R_BLR = 1.87 (+0.54 -0.50) light days and black hole mass of M_BH = 1.73 (+0.55 -0.52) x 10^6 M_sun. This radius is consistent with that expected from the R_BLR-L relationship, based on the present luminosity of NGC 4051 and the most current calibration of the relation by Bentz et al. (2009a). We also present a preliminary look at velocity-resolved Hbeta light curves and time delay measurements, although we are unable to reconstruct an unambiguous velocity-resolved reverberation signal.

Motivated by kinematic observations of young embedded clusters, this paper explores possible kinematic signatures produced by asphericity and departures from initial virial equilibrium in these systems. Specifically, the kinematic quantity that is measured and calculated in this study is the distribution of the line-of-sight velocities as a function of position along the cluster. Although clusters are found within a wide range of sizes, we focus on the regime with stellar membership $N \sim 10^3$. The gravitational potential of these young clusters is dominated by the gas, and the geometry of the gas distribution is generalized to include axisymmetric (and triaxial) forms. With this loss of symmetry, the kinematic results thus depend on viewing angle. This work also considers the stars to begin their trajectories with subvirial speeds, as indicated by observations of core motions in such clusters. Our results determine the conditions necessary for the kinematic signature to display interesting structure, i.e., a non-spherical potential, a viewing angle that is not along one of the principal axes, and subvirial starting conditions. We characterize the effects on this signature due to projection angle, initial stellar velocities, cluster elongation, and star formation efficiency. Finally, we compare our theoretical results to recent kinematic observations of the Orion Nebula Cluster; we find that the observations can be explained provided that the cluster is non-spherical, starts with subvirial initial velocities, and is not viewed along a principal axis.

We present a kinematic study of the Orion Nebula Cluster based upon radial velocities measured by multi-fiber echelle spectroscopy at the 6.5 meter MMT and Magellan telescopes. Velocities are reported for 1613 stars, with multi-epoch data for 727 objects as part of our continuing effort to detect and analyze spectroscopic binaries. We confirm and extend the results of Furesz et al. showing that the ONC is not relaxed, consistent with its youth, and that the stars generally follow the position-velocity structure of the moderate density gas in the region, traced by $^{13}$CO. The additional radial velocities we have measured enable us to probe some discrepancies between stellar and gaseous structure which can be attributed to binary motion and the inclusion of non-members in our kinematic sample. Our multi-epoch data allow us to identify 89 spectroscopic binaries; more will be found as we continue monitoring. Our results reinforce the idea that the ONC is a cluster in formation, and thus provides a valuable testing ground for theory. In particular, our observations are not consistent with the quasi-equilibrium or slow contraction models of cluster formation, but are consistent with cold collapse models.

We model Spitzer Space Telescope observations of the Taurus Class 0 protostar L1527 IRS (IRAS 04368+2557) to provide constraints on its protostellar envelope structure. The nearly edge-on inclination of L1527 IRS, coupled with the highly spatially-resolved near to mid-infrared images of this object and the detailed IRS spectrum, enable us to constrain the outflow cavity geometry quite well, reducing uncertainties in the other derived parameters. The mid-infrared scattered light image shows a bright central source within a dark lane; the aspect ratio of this dark lane is such that it appears highly unlikely to be a disk shadow. In modeling this dark lane, we conclude that L1527 IRS is probably not described by a standard TSC envelope with simple bipolar cavities. We find it necessary to model the dark lane and central source as a modified inner envelope structure. This structure may be due either to a complex wind-envelope interaction or induced by the central binary. To fit the overall SED, we require the central source to have a large near to mid-infrared excess, suggesting substantial disk accretion. Our model reproduces the overall morphology and surface brightness distribution of L1527 IRS fairly well, given the limitations of using axisymmetric models to fit the non-axisymmetric real object, and the derived envelope infall rates are in reasonable agreement with some other investigations. IRAC observations of L1527 IRS taken 12 months apart show variability in total flux and variability in the opposing bipolar cavities, suggesting asymmetric variations in accretion. We also provide model images at high resolution for comparison to future observations with current ground-based instrumentation and future space-based telescopes.

Deep Spitzer IRAC images of L1157 reveal many of the details of the outflow and the circumstellar environment of this Class 0 protostar. In IRAC band 4, 8 microns, there is a flattened structure seen in absorption against the background emission. The structure is perpendicular to the outflow and is extended to a diameter of 2 arcminutes. This structure is the first clear detection of a flattened circumstellar envelope or pseudo-disk around a Class 0 protostar. Such a flattened morphology is an expected outcome for many collapse theories that include magnetic fields or rotation. We construct an extinction model for a power-law density profile, but we do not constrain the density power-law index.

We present deep IRAC images that highlight the scattered light emission around many of the youngest protostars, the so-called Class 0 sources, in L1448. By comparison of the data with a Monte Carlo radiative transfer code \citep{whitney2003a}, we demonstrate for the first time that the observed infrared light from these objects is consistent with scattered light from the central protostar. The scattered light escapes out the cavity, carved by molecular outflows, in the circumstellar envelope. In particular, we observe prominent scattered light nebulae associated with the Class 0 sources: L1448-mm, L1448 IRS 2, and 3B, as well as a Class I source: IRS 3A. We use a grid of models with probable protostellar properties to generate model spectral energy distributions (SEDs) and images for bands sensitive to this scattered light: J, H, Ks, and \textit{Spitzer} IRAC bands. By simultaneously fitting SEDs and images of the outflow cavities, we are able to model geometric parameters, i.e. inclination angle and opening angle, and loosely constrain physical parameters. The opening angle may be an important indicator of the evolutionary state of a source. We compare our results for Class 0 sources to similar studies of Class I sources. There may be a transition phase from Class 0 to Class I when a source has an opening angle between 20$^o$ to 30$^o$. It is important to note that while the best fit model image and SED do not fully describe the sources, the fits generally describe the circumstellar structure of Class 0 sources in L1448.

We construct a simple model for radioisotopic enrichment of the protosolar nebula by injection from a nearby supernova, based on the inverse square law for ejecta dispersion. We find that the presolar radioisotopes abundances (i.e., in solar masses) demand a nearby supernova: its distance can be no larger than 66 times the size of the protosolar nebula, at a 90% confidence level, assuming 1 solar mass of protosolar material. The relevant size of the nebula depends on its state of evolution at the time of radioactivity injection. In one scenario, a collection of low-mass stars, including our sun, formed in a group or cluster with an intermediate- to high-mass star that ended its life as a supernova while our sun was still a protostar, a starless core, or perhaps a diffuse cloud. Using recent observations of protostars to estimate the size of the protosolar nebula constrains the distance of the supernova at 0.02 to 1.6 pc. The supernova distance limit is consistent with the scales of low-mass stars formation around one or more massive stars, but it is closer than expected were the sun formed in an isolated, solitary state. Consequently, if any presolar radioactivities originated via supernova injection, we must conclude that our sun was a member of such a group or cluster that has since dispersed, and thus that solar system formation should be understood in this context. In addition, we show that the timescale from explosion to the creation of small bodies was on the order of 1.8 Myr (formal 90% confidence range of 0 to 2.2 Myr), and thus the temporal choreography from supernova ejecta to meteorites is important. Finally, we can not distinguish between progenitor masses from 15 to 25 solar masses in the nucleosynthesis models; however, the 20 solar mass model is somewhat preferred.