FU Orionis objects (FUors) are rapidly-accreting, pre-main sequence objects that are known to exhibit large outbursts at optical and near-infrared wavelengths, with post-eruption, small-scale photometric variability superimposed on longer-term trends. In contrast, little is known about the variability of FUors at longer wavelengths. To explore this further, we observed six FUor objects using the NOrthern Extended Millimeter Array (NOEMA) and for a subset of three objects we obtained coordinated observations with NOEMA and the Lowell Discovery Telescope (LDT). In combination with previously published NOEMA observations from 2014, our 2017 observations of V1735 Cyg provide the first detection of variability in an FUor object at 2.7 mm. In the absence of significant optical variability, we discount the possibility that the mm flux density changed as a result of irradiation from the central disk. In addition, a change in the dust mass due to infall is highly unlikely. A plausible explanation for the change in 2.7 mm flux density is variability in free-free emission due to changes in the object's jet/wind. Thus, it may be that free-free emission in some FUor objects is significant at $\sim$3 mm and must be considered when deriving disk masses in order to help constrain the mechanism responsible for triggering FUor outbursts.

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.

We identify 45 protostellar outflows in CO maps of the Orion A giant molecular cloud from the CARMA-NRO Orion survey. Our sample includes 11 newly detected outflows. We measure the mass and energetics of the outflows, including material at low-velocities by correcting for cloud contributions. The total momentum and kinetic energy injection rates of outflows is comparable to the turbulent dissipation rate of the cloud. We also compare the outflow position angles to the orientation of C$^{18}$O filaments. We find that the full sample of outflows is consistent with being randomly oriented with respect to the filaments. A subsample of the most reliable measurements shows a moderately perpendicular outflow-filament alignment which may reflect accretion of mass across filaments and onto the protostellar cores.

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/.

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 report new Atacama Large Millimeter/submillimeter Array Band 3 (86-100 GHz; $\sim$80 mas angular resolution) and Band 4 (146-160 GHz; $\sim$50 mas angular resolution) observations of the dust continuum emission towards the archetypal and ongoing accretion burst young stellar object FU Ori, which simultaneously covered its companion, FU Ori S. In addition, we present near-infrared (2-2.45 $\mu$m) observations of FU Ori taken with the General Relativity Analysis via VLT InTerferometrY (GRAVITY; $\sim$1 mas angular resolution) instrument on the Very Large Telescope Interferometer (VLTI). We find that the emission in both FU Ori and FU Ori S at (sub)millimeter and near infrared bands is dominated by structures inward of $\sim$10 au radii. We detected closure phases close to zero from FU Ori with VLTI/GRAVITY, which indicate the source is approximately centrally symmetric and therefore is likely viewed nearly face-on. Our simple model to fit the GRAVITY data shows that the inner 0.4 au radii of the FU Ori disk has a triangular spectral shape at 2-2.45 $\mu$m, which is consistent with the H$_{2}$O and CO absorption features in a $\dot{M}\sim$10$^{-4}$ $M_{\odot}\,yr^{-1}$, viscously heated accretion disk. At larger ($\sim$0.4-10 au) radii, our analysis shows that viscous heating may also explain the observed (sub)millimeter and centimeter spectral energy distribution when we assume a constant, $\sim$10$^{-4}$ $M_{\odot}\,yr^{-1}$ mass inflow rate in this region. This explains how the inner 0.4 au disk is replenished with mass at a modest rate, such that it neither depletes nor accumulates significant masses over its short dynamic timescale. Finally, we tentatively detect evidence of vertical dust settling in the inner 10 au of the FU Ori disk, but confirmation requires more complete spectral sampling in the centimeter bands.

During star formation, the accretion disk drives fast MHD winds which usually contain two components, a collimated jet and a radially distributed wide-angle wind. These winds entrain the surrounding ambient gas producing molecular outflows. We report recent observation of 12CO (2-1) emission of the HH 46/47 molecular outflow by the Atacama Large Millimeter/sub-millimeter Array, in which we identify multiple wide-angle outflowing shell structures in both the blue and red-shifted outflow lobes. These shells are highly coherent in position-position-velocity space, extending to >40-50 km/s in velocity and 10^4 au in space with well defined morphology and kinematics. We suggest these outflowing shells are the result of the entrainment of ambient gas by a series of outbursts from an intermittent wide-angle wind. Episodic outbursts in collimated jets are commonly observed, yet detection of a similar behavior in wide-angle winds has been elusive. Here we show clear evidence that the wide-angle component of the HH 46/47 protostellar outflows experiences similar variability seen in the collimated component.

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.

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.

The dense clusters within the Serpens Molecular Cloud are among the most active regions of nearby star formation. In this paper, we use Gaia DR2 parallaxes and proper motions to statistically measure 1167 kinematic members of Serpens, few of which were previously identified, to evaluate the star formation history of the complex. The optical members of Serpens are concentrated in three distinct groups located at 380 to 480 pc; the densest clusters are still highly obscured by optically-thick dust and have few optical members. The total population of young stars and protostars in Serpens is at least 2000 stars, including past surveys that were most sensitive to protostars and disks, and may be far higher. Distances to dark clouds measured from deficits in star counts are consistent with the distances to the optical star clusters. The Serpens Molecular Cloud is seen in the foreground of the Aquila Rift, dark clouds located at 600 to 700 pc, and behind patchy extinction, here called the Serpens Cirrus, located at ~250 pc. Based on the lack of a distributed population of older stars, the star formation rate throughout the Serpens Molecular Cloud increased by at least a factor of 20 within the past ~5 Myr. The optically bright stars in Serpens Northeast are visible because their natal molecular cloud has been eroded and not because they were flung outwards from a central factory of star formation. The separation between subclusters of 20 to 100 pc and the absence of an older population leads to speculation that an external forcing was needed to trigger the active star formation.

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.

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.

Variability of pre-main-sequence stars observed at optical wavelengths has been attributed tofluctuations in the mass accretion rate from the circumstellar disk onto the forming star.Detailed models of accretion disks suggest that young deeply embedded protostars should also exhibit variations in their accretion rates, and that these changes can be tracked indirectly by monitoring the response of the dust envelope at mid-IR to millimeter wavelengths. Interferometers such as ALMA offer the resolution and sensitivity to observe small fluctuations in brightness at the scale of the disk where episodic accretion may be driven. In this work, we present novel methods for comparing interferometric observations and apply them to CARMA and ALMA 1.3mm observations of deeply embedded protostars in Serpens taken 9 years apart. We find no brightness variation above the limits of our analysis of a factor of $\gtrsim 50\%$, due to the limited sensitivity of the CARMA observations and small number of sources common to both epochs. We further show that follow up ALMA observations with a similar sample size and sensitivity may be able to uncover variability at the level of a few percent, and discuss implications for future work.

A good constraint of when the growth of dust grains from sub-micrometer to millimeter sizes occurs, is crucial for planet formation models. This provides the first step towards the production of pebbles and planetesimals in protoplanetary disks. Currently, it is well established that Class II objects have large dust grains. However, it is not clear when in the star formation process this grain growth occurs. We use multi-wavelength millimeter observations of a Class I protostar to obtain the spectral index of the observed flux densities $\alpha_\mathrm{mm}$ of the unresolved disk and the surrounding envelope. Our goal is to compare our observational results with visibility modeling at both wavelengths simultaneously. We present data from NOEMA at 2.7 mm and SMA at 1.3 mm of the Class I protostar, Per-emb-50. We model the dust emission with a variety of parametric and radiative transfer models to deduce the grain size from the observed emission spectral index. We find a spectral index in the envelope of Per-emb-50 of $\alpha_{\rm env}$=$3.3\pm0.3$, similar to the typical ISM values. The radiative transfer modeling of the source confirms this value of $\alpha_{\rm env}$ with the presence of dust with a $a_\mathrm{max}$$\leq$100 $\mu$m. Additionally, we explore the backwarming effect, where we find that the envelope structure affects the millimeter emission of the disk. Our results reveal grains with a maximum size no larger than $100$ $\mu$m in the inner envelope of the Class I protostar Per-emb-50, providing an interesting case to test the universality of millimeter grain growth expected in these sources.

We analyze high-quality stellar catalogs for 24 young and nearby (within 1 kpc) embedded clusters and present a catalogue of 32 groups which have a high concentration of protostars. The median effective radius of these groups is 0.17 pc. The median protostellar and pre-main sequence star surface densities are 46 M_{\odot} pc^{-2} and 11 M_{\odot} pc^{-2}, respectively. We estimate the age of these groups using a model of constant birthrate and random accretion stopping and find a median value of 0.25 Myr. Some groups in Aquila and Serpens, Corona Australia and Ophichus L1688 show high protostellar surface density and high molecular gas surface density, which seem to be undergoing vigorous star formation. These groups provide an excellent opportunity to study initial conditions of clustered star formation. Comparison of protostellar and pre-main-sequence stellar surface densities reveal continuous low-mass star formation of these groups over several Myr in some clouds. For groups with typical protostellar separations of less than 0.4 pc, we find that these separations agree well with the thermal Jeans fragmentation scale. On the other hand, for groups with typical protostellar separations larger than 0.4 pc, these separations are always larger than the associated Jeans length.

Episodic accretion may be a common occurrence in the evolution of young pre-main sequence stars and has important implications for our understanding of star and planet formation. Many fundamental aspects of what drives the accretion physics, however, are still unknown. The ngVLA will be a key tool in understanding the nature of these events. The high spatial resolution, broad spectral coverage, and unprecedented sensitivity will allow for the detailed analysis of outburst systems. The proposed frequency range of the ngVLA allows for observations of the gas, dust, and non-thermal emission from the star and disk.

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.

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 present moderate (${\sim}5^{\prime\prime}$) and high angular resolution (${\sim}1^{\prime\prime}$) observations of $^{12}\rm{CO\,}(J=2-1)$ emission toward nearby, interacting galaxy NGC 3627 taken with the Submillimeter Array (SMA). These SMA mosaic maps of NGC 3627 reveal a prominent nuclear peak, inter-arm regions, and diffuse, extended emission in the spiral arms. A velocity gradient of ${\sim}400$-$450$ km s$^{-1}$ is seen across the entire galaxy with velocity dispersions ranging from $\lesssim 80$ km s$^{-1}$ toward the nuclear region to $\lesssim 15$ km s$^{-1}$ in the spiral arms. We also detect unresolved $^{13}\rm{CO\,}(J=2-1)$ line emission toward the nuclear region, southern bar end, and in a relatively isolated clump in the southern portion of the galaxy, while no $\rm{C}^{18}O(J=2-1)$ line emission is detected at a $3\sigma$ rms noise level of 42 mJy beam$^{-1}$ per 20 km s$^{-1}$ channel. Using RADEX modeling with a large velocity gradient approximation, we derive kinetic temperatures ranging from ${\sim}5$-$10$ K (in the spiral arms) to ${\sim}25$ K (at the center) and H$_2$ number densities from ${\sim}$400-1000 cm$^{-3}$ (in the spiral arms) to ${\sim}$12500 cm$^{-3}$ (at the center). From this density modeling, we find a total H$_2$ mass of $9.6\times10^9 M_{\odot}$, which is ${\sim}50\%$ higher than previous estimates made using a constant H$_2$-CO conversion factor but is largely dependent on the assumed vertical distribution of the CO gas. With the exception of the nuclear region, we also identify a tentative correlation between star formation efficiency and kinetic temperature. We derive a galactic rotation curve, finding a peak velocity of ${\sim}207$ km s$^{-1}$ and estimate a total dynamical mass of $4.94 \pm 0.70 \times 10^{10} M_{\odot}$ at a galactocentric radius of ${\sim}6.2$ kpc ($121^{\prime\prime}$).

We present near-IR imaging polarimetry of five classical FU Ori-type objects (FU Ori, V1057 Cyg, V1515 Cyg, V1735 Cyg, Z CMa) with a $\sim$0\farcs1 resolution observed using HiCIAO+AO188 at Subaru Telescope. We observed scattered light associated with circumstellar dust around four of them (i.e., all but V1515 Cyg). Their polarized intensity distribution shows a variety of morphologies with arms, tails or streams, spikes and fragmented distributions, many of which were reported in our previous paper. The morphologies of these reflection nebulae significantly differ from many other normal young stellar objects (Class I-II objects). These structures are attributed to gravitationally unstable disks, trails of clump ejections, dust blown by a wind or a jet, and a stellar companion. We can consistently explain our results with the scenario that their accretion outbursts (FUor outbursts) are triggered by gravitationally fragmenting disks, and with the hypothesis that many low-mass young stellar objects experience such outbursts.

(Abridged) Protostellar systems evolve from prestellar cores, through the deeply embedded stage and then disk-dominated stage, before they end up on the main sequence. Knowing how much time a system spends in each stage is crucial for understanding how stars and associated planetary systems form, because a key constraint is the time available to form such systems. Equally important is understanding what the spread in these time scales is. The most commonly used method for inferring protostellar ages is to assume the lifetime of one evolutionary stage, and then scale this to the relative number of protostars in the other stages, i.e., assuming steady state. This method does not account for the underlying age distribution and apparent stochasticity of star formation, nor that relative populations are not in steady state. To overcome this, we propose a new scheme where the lifetime of each protostellar stage follows a distribution based on the formalism of sequential nuclear decay. The main assumptions are: Class 0 sources follow a straight path to Class III sources, the age distribution follows a binomial distribution, and the star-formation rate is constant. The results are that the half-life of Class 0, Class I, and Flat sources are (2.4+/-0.2)%, (4.4+/-0.3)%, and (4.3+/-0.4)% of the Class II half-life, respectively, which translates to 47+/-4, 88+/-7, and 87+/-8 kyr, respectively, for a Class II half-life of 2 Myr for protostars in the Gould Belt clouds with more than 100 protostars. The mean age of these clouds is 1.2+/-0.1 Myr, and the star formation rate is (8.3+/-0.5)x10^-4 Msun/yr. The critical parameters in arriving at these numbers are the assumed half-life of the Class II stage, and the assumption that the star-formation rate and half-lives are constant. This method presents a first step in moving from steady-state to non-steady-state solutions of protostellar populations.

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.

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

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 aim at studying the causal link between the knotty jet structure in CARMA 7, a young Class 0 protostar in the Serpens South cluster, and episodic accretion in young protostellar disks. We used numerical hydrodynamics simulations to derive the protostellar accretion history in gravitationally unstable disks around solar-mass protostars. We compared the time spacing between luminosity bursts \Delta\tau_mod, caused by dense clumps spiralling on the protostar, with the differences of dynamical timescales between the knots \Delta\tau_obs in CARMA 7. We found that the time spacing between the bursts have a bi-modal distribution caused by isolated and clustered luminosity bursts. The former are characterized by long quiescent periods between the bursts with \Delta\tau_mod = a few * (10^3-10^4) yr, whereas the latter occur in small groups with time spacing between the bursts \Delta\tau_mod= a few * (10-10^2) yr. For the clustered bursts, the distribution of \Delta\tau_mod in our models can be fit reasonably well to the distribution of \Delta\tau_obs in the protostellar jet of CARMA 7, if a certain correction for the (yet unknown) inclination angle with respect to the line of sight is applied. The K-S test on the model and observational data sets suggests the best-fit values for the inclination angles of 55-80 deg., which become narrower (75-80 deg.) if only strong luminosity bursts are considered. The dynamical timescales of the knots in the jet of CARMA 7 are too short for a meaningful comparison with the long time spacings between isolated bursts in our models. The exact sequences of time spacings between the luminosity bursts in our models and knots in the jet of CARMA 7 were found difficult to match. (abridged)

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.

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.

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.

Young stellar objects (YSOs) may undergo periods of active accretion (outbursts), during which the protostellar accretion rate is temporarily enhanced by a few orders of magnitude. Whether or not these accretion outburst YSOs possess similar dust/gas reservoirs to each other, and whether or not their dust/gas reservoirs are similar as quiescent YSOs, are issues not yet clarified. The aim of this work is to characterize the millimeter thermal dust emission properties of a statistically significant sample of long and short duration accretion outburst YSOs (i.e., FUors and EXors) and the spectroscopically identified candidates of accretion outbursting YSOs (i.e., FUor-like objects). We have carried out extensive Submillimeter Array (SMA) observations mostly at $\sim$225 GHz (1.33 mm) and $\sim$272 GHz (1.10 mm), from 2008 to 2017. We covered accretion outburst YSOs located at $<$1 kpc distances from the solar system. We analyze all the existing SMA data of such objects, both published and unpublished, in a coherent way to present a millimeter interferometric database of 29 objects. We obtained 21 detections at $>$3-$\sigma$ significance. Detected sources except for the two cases of V883 Ori and NGC 2071 MM3 were observed with $\sim$1$"$ angular resolution. Overall our observed targets show a systematically higher millimeter luminosity distribution than those of the $M_{*}>$0.3 $M_{\odot}$ Class II YSOs in the nearby ($\lesssim$400 pc) low-mass star-forming molecular clouds (e.g., Taurus, Lupus, Upp Scorpio, and Chameleon I). In addition, at 1 mm our observed confirmed binaries or triple-system sources are systematically fainter than the rest of the sources even though their 1 mm fluxes are broadly distributed. We may have detected $\sim$30-60\% millimeter flux variability from V2494 Cyg and V2495 Cyg, from the observations separated by $\sim$1 year.

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

Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion -- prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In "A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star Forming Regions", we are monitoring monthly with SCUBA-2 the sub-mm emission in eight fields within nearby (<500 pc) star forming regions to measure the accretion variability of protostars. The total survey area of ~1.6 sq.deg. includes ~105 peaks with peaks brighter than 0.5 Jy/beam (43 associated with embedded protostars or disks) and 237 peaks of 0.125-0.5 Jy/beam (50 with embedded protostars or disks). Each field has enough bright peaks for flux calibration relative to other peaks in the same field, which improves upon the nominal flux calibration uncertainties of sub-mm observations to reach a precision of ~2-3% rms, and also provides quantified confidence in any measured variability. The timescales and amplitudes of any sub-mm variation will then be converted into variations in accretion rate and subsequently used to infer the physical causes of the variability. This survey is the first dedicated survey for sub-mm variability and complements other transient surveys at optical and near-IR wavelengths, which are not sensitive to accretion variability of deeply embedded protostars.

A long-standing open issue of the paradigm of low-mass star formation is the luminosity problem: most protostars are less luminous than theoretically predicted. One possible solution is that the accretion process is episodic. FU Ori-type stars (FUors) are thought to be the visible examples for objects in the high accretion state. FUors are often surrounded by massive envelopes, which replenish the disk material and enable the disk to produce accretion outbursts. However, we have insufficient information on the envelope dynamics in FUors, about where and how mass transfer from the envelope to the disk happens. Here we present ALMA observations of the FUor-type star V346 Nor at 1.3 mm continuum and in different CO rotational lines. We mapped the density and velocity structure of its envelope and analyze the results using channel maps, position-velocity diagrams, and spectro-astrometric methods. We found that V346 Nor is surrounded by gaseous material on 10000 au scale in which a prominent outflow cavity is carved. Within the central $\sim$700 au, the circumstellar matter forms a flattened pseudo-disk where material is infalling with conserved angular momentum. Within $\sim$350 au, the velocity profile is more consistent with a disk in Keplerian rotation around a central star of 0.1 $M_{\odot}$. We determined an infall rate from the envelope onto the disk of 6$\times$10$^{-6}\,M_{\odot}$yr$^{-1}$, a factor of few higher than the quiescent accretion rate from the disk onto the star, hinting for a mismatch between the infall and accretion rates as the cause of the eruption.

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.

We present high angular resolution dust polarization and molecular line observations carried out with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the Class 0 protostar Serpens SMM1. By complementing these observations with new polarization observations from the Submillimeter Array (SMA) and archival data from the Combined Array for Research in Millimeter-wave Astronomy (CARMA) and the James Clerk Maxwell Telescopes (JCMT), we can compare the magnetic field orientations at different spatial scales. We find major changes in the magnetic field orientation between large (~0.1 pc) scales -- where the magnetic field is oriented E-W, perpendicular to the major axis of the dusty filament where SMM1 is embedded -- and the intermediate and small scales probed by CARMA (~1000 AU resolution), the SMA (~350 AU resolution), and ALMA (~140 AU resolution). The ALMA maps reveal that the redshifted lobe of the bipolar outflow is shaping the magnetic field in SMM1 on the southeast side of the source; however, on the northwestern side and elsewhere in the source, low velocity shocks may be causing the observed chaotic magnetic field pattern. High-spatial-resolution continuum and spectral-line observations also reveal a tight (~130 AU) protobinary system in SMM1-b, the eastern component of which is launching an extremely high-velocity, one-sided jet visible in both CO(2-1) and SiO(5-4); however, that jet does not appear to be shaping the magnetic field. These observations show that with the sensitivity and resolution of ALMA, we can now begin to understand the role that feedback (e.g., from protostellar outflows) plays in shaping the magnetic field in very young, star-forming sources like SMM1.

Understanding the formation of wide binary systems of very low mass stars (M $\le$ 0.1 Msun) is challenging. The most obvious route is via widely separated low-mass collapsing fragments produced through turbulent fragmentation of a molecular core. However, close binaries/multiples from disk fragmentation can also evolve to wide binaries over a few initial crossing times of the stellar cluster through tidal evolution. Finding an isolated low mass wide binary system in the earliest stage of formation, before tidal evolution could occur, would prove that turbulent fragmentation is a viable mechanism for (very) low mass wide binaries. Here we report high resolution ALMA observations of a known wide-separation protostellar binary, showing that each component has a circumstellar disk. The system is too young to have evolved from a close binary and the disk axes are misaligned, providing strong support for the turbulent fragmentation model. Masses of both stars are derived from the Keplerian rotation of the disks; both are very low mass stars.

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.

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 resolved FU Ori at 29-37 GHz using the JVLA with $\sim$0$''$.07 resolution, and performed the complementary JVLA 8-10 GHz observations, the SMA 224 GHz and 272 GHz observations, and compared with archival ALMA 346 GHz observations to obtain the SEDs. Our 8-10 GHz observations do not find evidence for the presence of thermal radio jets, and constrain the radio jet/wind flux to at least 90 times lower than the expected value from the previously reported bolometric luminosity-radio luminosity correlation. The emission at $>$29 GHz may be dominated by the two spatially unresolved sources, which are located immediately around FU Ori and its companion FU Ori S, respectively. Their deconvolved radii at 33 GHz are only a few au. The 8-346 GHz SEDs of FU Ori and FU Ori S cannot be fit by constant spectral indices (over frequency). The more sophisticated models for SEDs suggest that the $>$29 GHz emission is contributed by a combination of free-free emission from ionized gas, and thermal emission from optically thick and optically thin dust components. We hypothesize that dust in the innermost parts of the disks ($\lesssim$0.1 au) has been sublimated, and thus the disks are no more well shielded against the ionizing photons. The estimated overall gas and dust mass based on SED modeling, can be as high as a fraction of a solar mass, which is adequate for developing disk gravitational instability. Our present explanation for the observational data is that the massive inflow of gas and dust due to disk gravitational instability or interaction with a companion/intruder, was piled up at the few au scale due to the development of a deadzone with negligible ionization. The piled up material subsequently triggered the thermal and the MRI instabilities when the ionization fraction in the inner sub-au scale region exceeded a threshold value, leading to the high protostellar accretion rate.

[Abridged] We present spectroscopic observations in H$_{2}$O, CO and related species with \textit{Herschel} HIFI and PACS, as well as ground-based follow-up with the JCMT and APEX in CO, HCO$^{+}$ and isotopologues, of a sample of 49 nearby ($d<$500\,pc) candidate protostars. These data are used to study the outflow and envelope properties of these sources. We also compile their continuum SEDs in order to constrain their physical properties. Water emission is dominated by shocks associated with the outflow, rather than the cooler, slower entrained outflowing gas probed by ground-based CO observations. These shocks become less energetic as sources evolve from Class 0 to Class I. The fraction of mass in the outflow relative to the total envelope (i.e. $M_{\mathrm{out}}/M_{\mathrm{env}}$) remains broadly constant between Class 0 and I. The median value ($\sim$1$\%$) is consistent with a core to star formation efficiency on the order of 50$\%$ and an outflow duty cycle on the order of 5$\%$. Entrainment efficiency, as probed by $F_{\mathrm{CO}}/\dot{M}_{\mathrm{acc}}$, is also invariant with source properties and evolutionary stage. The median value (6.3\kms{}) suggests an entrainment efficiency of between 30 and 60$\%$ if the wind is launched at $\sim$1AU. $L$[O\,{\sc i}] is strongly correlated with $L_{\mathrm{bol}}$ but not with $M_{\mathrm{env}}$, while low-$J$ CO is more closely correlated with the latter than the former. This suggests that [O\,{\sc i}] traces the present-day accretion activity while CO traces time-averaged accretion over the dynamical timescale of the outflow. $L$[O\,{\sc i}] does not vary from Class 0 to Class I, unlike CO and H$_{2}$O. This is likely due to the ratio of atomic to molecular gas in the wind increasing as the source evolves, balancing out the decrease in mass accretion rate. Infall signatures are detected in HCO$^{+}$ and H$_{2}$O in a few sources.

We use Herschel spectrophotometry of BHR71, an embedded Class 0 protostar, to provide new constraints on its physical properties. We detect 645 (non-unique) spectral lines amongst all spatial pixels. At least 61 different spectral lines originate from the central region. A CO rotational diagram analysis shows four excitation temperature components, 43 K, 197 K, 397 K, and 1057 K. Low-J CO lines trace the outflow while the high-J CO lines are centered on the infrared source. The low-excitation emission lines of H2O trace the large-scale outflow, while the high-excitation emission lines trace a small-scale distribution around the equatorial plane. We model the envelope structure using the dust radiative transfer code, Hyperion, incorporating rotational collapse, an outer static envelope, outflow cavity, and disk. The evolution of a rotating collapsing envelope can be constrained by the far-infrared/millimeter SED along with the azimuthally-averaged radial intensity profile, and the structure of the outflow cavity plays a critical role at shorter wavelengths. Emission at 20-40 um requires a cavity with a constant-density inner region and a power-law density outer region. The best fit model has an envelope mass of 19 solar mass inside a radius of 0.315 pc and a central luminosity of 18.8 solar luminosity. The time since collapse began is 24630-44000 yr, most likely around 36000 yr. The corresponding mass infall rate in the envelope (1.2x10$^{-5}$ solar mass per year) is comparable to the stellar mass accretion rate, while the mass loss rate estimated from the CO outflow is 20% of the stellar mass accretion rate. We find no evidence for episodic accretion.

The nature of very low luminosity objects with the internal luminosity L_obj <= 0.1 Lsun is investigated by means of numerical modeling coupling the core collapse simulations with the stellar evolution calculations. The gravitational collapse of a large sample of model cores in the mass range 0.1-2.0 Msun is investigated. Numerical simulations were started at the pre-stellar phase and terminated at the end of the embedded phase when 90% of the initial core mass had been accreted onto the forming protostar plus disk system. The disk formation and evolution was studied using numerical hydrodynamics simulations, while the formation and evolution of the central star was calculated using a stellar evolution code. Three scenarios for mass accretion from the disk onto the star were considered: hybrid accretion in which a fraction of accreted energy absorbed by the protostar depends on the accretion rate, hot accretion wherein a fraction of accreted energy is constant, and cold accretion wherein all accretion energy is radiated away. Our conclusions on the nature of VeLLOs depend crucially on the character of protostellar accretion. In the hybrid accretion scenario, most VeLLOs (90.6%) are expected to be the first hydrostatic cores (FHSCs) and only a small fraction (9.4%) are true protostars. In the hot accretion scenario, all VeLLOs are FHSCs due to overly high photospheric luminosity of protostars. In the cold accretion scenario, on the contrary, the majority of VeLLOs belong to the Class I phase of stellar evolution. The reason is that the stellar photospheric luminosity, which sets the floor for the total internal luminosity of a young star, is lower in cold accretion, thus enabling more VeLLOs in the protostellar stage. VeLLOs are relatively rare objects occupying 7%-11% of the total duration of the embedded phase and their masses do not exceed 0.3 Msun. (abridged).

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

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.

Chandra observations have revealed a feature,within 16" of Cygnus X-3 which varied in phase with Cygnus X-3. This feature was shown to be a Bok globule which is along the line of sight to Cygnus X-3. We report on observations made with Submillimeter Array (SMA) to search for molecular emission from this globule, also known as Cygnus X-3's "Little Friend" (LF). We have found a counterpart in both 12CO (2-1) and 13}CO (2-1) emission. From the velocity shift of the molecular lines, we are able to find two probable distances based on the Bayesian model of Milky Way kinematics of Reid et al. 2016. For the LF velocity of -47.5 km/s, we get a distance of 6.1 +/- 0.6 kpc (62% probability) and 7.8 +/- 0.6 ~kpc (38% probability). This yields a distance to Cyg X-3 of 7.4 +/- 1.1 kpc and 10.2 +/-1.2 ~kpc, respectively. Based on the probabilities entailed, we take 7.4 +/- 1.1 kpc as the preferred distance to Cyg X-3. We also report the discovery of bipolar molecular outflow, suggesting that there is active star formation occurring within the Little Friend.

HBC 722 (V2493 Cyg) is a young eruptive star in outburst since 2010. It is an FU Orionis-type object with an atypically low outburst luminosity. Because it was well characterized in the pre-outburst phase, HBC 722 is one of the few FUors where we can learn about the physical changes and processes associated with the eruption. We monitored the source in the BVRIJHKs bands from the ground, and at 3.6 and 4.5 $\mu$m from space with the Spitzer Space Telescope. We analyzed the light curves and the spectral energy distribution by fitting a series of steady accretion disk models at many epochs. We also analyzed the spectral properties of the source based on new optical and infrared spectra. We also mapped HBC 722 and its surroundings at millimeter wavelengths. From the light curve analysis we concluded that the first peak of the outburst in 2010 September was due to an abrupt increase of the accretion rate in the innermost part of the system. This was followed by a long term process, when the brightening was mainly due to a gradual increase of the accretion rate and the emitting area. Our new observations show that the source is currently in a constant plateau phase. We found that around the peak the continuum was bluer and the H$\alpha$ profile changed significantly between 2012 and 2013. The source was not detected in the millimeter continuum, but we discovered a flattened molecular gas structure with a diameter of 1700 au and mass of 0.3 M$_{\odot}$ centered on HBC 722. While the first brightness peak could be interpreted as a rapid fall of piled-up material from the inner disk onto the star, the later monotonic flux rise suggests the outward expansion of a hot component according to the theory of Bell & Lin (1994). Our study of HBC 722 demonstrated that accretion-related outbursts can occur in young stellar objects even with very low mass disks, in the late Class II phase.

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.

We present ALMA Cycle 1 observations of the HH46/47 molecular outflow using combined 12m array and ACA observations. The improved angular resolution and sensitivity of our multi-line maps reveal structures that help us study the entrainment process in much more detail and allow us to obtain more precise estimates of outflow properties than previous observations. We use 13CO(1-0) and C18O(1-0) emission to correct for the 12CO(1-0) optical depth to accurately estimate the outflow mass, momentum and kinetic energy. This correction increases the estimates of the mass, momentum and kinetic energy by factors of about 9, 5 and 2, respectively, with respect to estimates assuming optically thin emission. The new 13CO and C18O data also allow us to trace denser and slower outflow material than that traced by the 12CO maps, and they reveal an outflow cavity wall at very low velocities (as low as 0.2km/s with respect to the cores central velocity). Adding with the slower material traced only by 13CO and C18O, there is another factor of 3 increase in the mass estimate and 50% increase in the momentum estimate. The estimated outflow properties indicate that the outflow is capable of dispersing the parent core within the typical lifetime of the embedded phase of a low-mass protostar, and that it is responsible for a core-to-star efficiency of 1/4 to 1/3. We find that the outflow cavity wall is composed of multiple shells associated with a series of jet bow-shock events. Within about 3000AU of the protostar the 13CO and C18O emission trace a circumstellar envelope with both rotation and infall motions, which we compare with a simple analytic model. The CS(2-1) emission reveals tentative evidence of a slowly-moving rotating outflow, which we suggest is entrained not only poloidally but also toroidally by a disk wind that is launched from relatively large radii from the source.

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.

The protostellar outflows of wide-separation forming binaries frequently appear misaligned. We use magneto-hydrodynamic simulations to investigate the alignment of protostellar spin and molecular outflows for forming binary pairs. We show that the protostellar pairs, which form from turbulent fragmentation within a single parent core, have randomly oriented angular momenta. Although the pairs migrate to closer separations, their spins remain partially misaligned. We produce $^{12}$CO(2-1) synthetic observations of the simulations and characterize the outflow orientation in the emission maps. The CO-identified outflows exhibit a similar random distribution and are also statistically consistent with the observed distribution of molecular outflows. We conclude that observed misalignment provides a clear signature of binary formation via turbulent fragmentation. The persistence of misaligned outflows and stellar spins following dynamical evolution may provide a signature of binary origins for more evolved multiple star systems.

We report ALMA observations of a one-sided, high-velocity ($\sim$80 km s$^{-1}$) CO($J = 2 \rightarrow 1$) jet powered by the intermediate-mass protostellar source Serpens SMM1-a. The highly collimated molecular jet is flanked at the base by a wide-angle cavity; the walls of the cavity can be seen in both 4 cm free-free emission detected by the VLA and 1.3 mm thermal dust emission detected by ALMA. This is the first time that ionization of an outflow cavity has been directly detected via free-free emission in a very young, embedded Class 0 protostellar source that is still powering a molecular jet. The cavity walls are ionized either by UV photons escaping from the accreting protostellar source, or by the precessing molecular jet impacting the walls. These observations suggest that ionized outflow cavities may be common in Class 0 protostellar sources, shedding further light on the radiation, outflow, and jet environments in the youngest, most embedded forming stars.

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 a catalog of low-mass dense cores observed with the SHARC-II instrument at 350 microns. Our observations have an effective angular resolution of 10", approximately 2.5 times higher than observations at the same wavelength obtained with the Herschel Space Observatory, albeit with lower sensitivity, especially to extended emission. The catalog includes 81 maps covering a total of 164 detected sources. For each detected source, we tabulate basic source properties including position, peak intensity, flux density in fixed apertures, and radius. We examine the uncertainties in the pointing model applied to all SHARC-II data and conservatively find that the model corrections are good to within ~3", approximately 1/3 of the SHARC-II beam. We examine the differences between two array scan modes and find that the instrument calibration, beam size, and beam shape are similar between the two modes. We also show that the same flux densities are measured when sources are observed in the two different modes, indicating that there are no systematic effects introduced into our catalog by utilizing two different scan patterns during the course of taking observations. We find a detection rate of 95% for protostellar cores but only 45% for starless cores, and demonstrate the existence of a SHARC-II detection bias against all but the most massive and compact starless cores. Finally, we discuss the improvements in protostellar classification enabled by these 350 micron observations.

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.

Young stellar objects (YSOs) may not accumulate their mass steadily, as was previously thought, but in a series of violent events manifesting themselves as sharp stellar brightening. These events can be caused by fragmentation due to gravitational instabilities in massive gaseous disks surrounding young stars, followed by migration of dense gaseous clumps onto the star. We report our high angular resolution, coronagraphic near-infrared polarization imaging observations using the High Contrast Instrument for the Subaru Next Generation Adaptive Optics (HiCIAO) of the Subaru 8.2 m Telescope, towards four YSOs which are undergoing luminous accretion outbursts. The obtained infrared images have verified the presence of several hundred AUs scale arms and arcs surrounding these YSOs. In addition, our hydrodynamics simulations and radiative transfer models further demonstrate that these observed structures can indeed be explained by strong gravitational instabilities occurring at the beginning of the disk formation phase. The effect of those tempestuous episodes of disk evolution on star and planet formation remains to be understood.

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 report Submillimeter Array (SMA) 1.3 mm high angular resolution observations towards the four EXor type outbursting young stellar objects (YSOs) VY Tau, V1118 Ori, V1143 Ori, and NY Ori. The data mostly show low dust masses $M_{dust}$ in the associated circumstellar disks. Among the sources, NY Ori possesses a relatively massive disk with $M_{dust} \sim 9 \times 10^{-4}$ $M_{\odot}$. V1118 Ori has a marginal detection equivalent to $M_{dust} \sim 6 \times 10^{-5}$ $M_{\odot}$. V1143 Ori has a non-detection also equivalent to $M_{dust} < 6 \times 10^{-5}$ $M_{\odot}$. For the nearest source VY Tau, we get a surprising non-detection which provides a stringent upper limit $M_{dust} < 6 \times 10^{-6}$ $M_{\odot}$. We interpret our findings as suggesting that the gas and dust reservoirs that feed the short duration, repetitive optical outbursts seen in some EXors may be limited to the small scale, innermost region of their circumstellar disks. This hot dust may have escaped our detection limits. Follow-up, more sensitive millimeter observations are needed to improve our understanding of the triggering mechanisms of EXor type outbursts.

We present millimetre (SMA) and sub-millimetre (SCUBA-2) continuum observations of the peculiar star KIC 8462852 which displayed several deep and aperiodic dips in brightness during the Kepler mission. Our observations are approximately confusion-limited at 850 $\mu$m and are the deepest millimetre and sub-millimetre photometry of the star that has yet been carried out. No significant emission is detected towards KIC 8462852. We determine upper limits for dust between a few 10$^{-6}$ M$_{\oplus}$ and 10$^{-3}$ M$_{\oplus}$ for regions identified as the most likely to host occluding dust clumps and a total overall dust budget of $<$7.7 M$_{\oplus}$ within a radius of 200 AU. Such low limits for the inner system make the catastrophic planetary disruption hypothesis unlikely. Integrating over the Kepler lightcurve we determine that at least 10$^{-9}$ M$_{\oplus}$ of dust is required to cause the observed Q16 dip. This is consistent with the currently most favoured cometary breakup hypothesis, but nevertheless implies the complete breakup of $\sim$ 30 Comet 1/P Halley type objects. Finally, in the wide SCUBA-2 field-of-view we identify another candidate debris disc system that is potentially the largest yet discovered.

We sort $4683$ molecular clouds between $10^\circ< \ell <65^\circ$ from the Bolocam Galactic Plane Survey based on observational diagnostics of star formation activity: compact $70$ $\mu{\rm m}$ sources, mid-IR color-selected YSOs, ${\rm H_2O}$ and ${\rm CH_3OH}$ masers, and UCHII regions. We also present a combined ${\rm NH_3}$-derived gas kinetic temperature and ${\rm H_2O}$ maser catalog for $1788$ clumps from our own GBT 100m observations and from the literature. We identify a subsample of $2223$ ($47.5\%$) starless clump candidates, the largest and most robust sample identified from a blind survey to date. Distributions of flux density, flux concentration, solid angle, kinetic temperature, column density, radius, and mass show strong ($>1$ dex) progressions when sorted by star formation indicator. The median starless clump candidate is marginally sub-virial ($\alpha \sim 0.7$) with $>75\%$ of clumps with known distance being gravitationally bound ($\alpha < 2$). These samples show a statistically significant increase in the median clump mass of $\Delta M \sim 170-370$ M$_\odot$ from the starless candidates to clumps associated with protostars. This trend could be due to (i) mass growth of the clumps at $\dot{M}\sim200-440$ Msun Myr$^{-1}$ for an average free-fall $0.8$ Myr time-scale, (ii) a systematic factor of two increase in dust opacity from starless to protostellar phases, (iii) and/or a variation in the ratio of starless to protostellar clump lifetime that scales as $\sim M^{-0.4}$. By comparing to the observed number of ${\rm CH_3OH}$ maser containing clumps we estimate the phase-lifetime of massive ($M>10^3$ M$_\odot$) starless clumps to be $0.37 \pm 0.08 \ {\rm Myr} \ (M/10^3 \ {\rm M}_\odot)^{-1}$; the majority ($M<450$ M$_\odot$) have phase-lifetimes longer than their average free-fall time.

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.

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

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.

We present the full catalog of Young Stellar Objects (YSOs) identified in the 18 molecular clouds surveyed by the Spitzer Space Telescope "cores to disks" (c2d) and "Gould Belt" (GB) Legacy surveys. Using standard techniques developed by the c2d project, we identify 3239 candidate YSOs in the 18 clouds, 2966 of which survive visual inspection and form our final catalog of YSOs in the Gould Belt. We compile extinction corrected SEDs for all 2966 YSOs and calculate and tabulate the infrared spectral index, bolometric luminosity, and bolometric temperature for each object. We find that 326 (11%), 210 (7%), 1248 (42%), and 1182 (40%) are classified as Class 0+I, Flat-spectrum, Class II, and Class III, respectively, and show that the Class III sample suffers from an overall contamination rate by background AGB stars between 25% and 90%. Adopting standard assumptions, we derive durations of 0.40-0.78 Myr for Class 0+I YSOs and 0.26-0.50 Myr for Flat-spectrum YSOs, where the ranges encompass uncertainties in the adopted assumptions. Including information from (sub)millimeter wavelengths, one-third of the Class 0+I sample is classified as Class 0, leading to durations of 0.13-0.26 Myr (Class 0) and 0.27-0.52 Myr (Class I). We revisit infrared color-color diagrams used in the literature to classify YSOs and propose minor revisions to classification boundaries in these diagrams. Finally, we show that the bolometric temperature is a poor discriminator between Class II and Class III YSOs.

The Spitzer Space Telescope mapped the Perseus molecular cloud complex with IRAC and MIPS as part of the c2d Spitzer Legacy project. This paper combines the observations from both instruments giving an overview of low-mass star formation across Perseus from 3.6 to 70 micron. We provide an updated list of young stellar objects with new classifications and source fluxes from previous works, identifying 369 YSOs in Perseus with the Spitzer dataset. By synthesizing the IRAC and MIPS maps of Perseus and building on the work of previous papers in this series (Jorgensen et al. 2006, Rebull et al. 2007), we present a current census of star formation across the cloud and within smaller regions. 67% of the YSOs are associated with the young clusters NGC 1333 and IC 348. The majority of the star formation activity in Perseus occurs in the regions around the clusters, to the eastern and western ends of the cloud complex. The middle of the cloud is nearly empty of YSOs despite containing regions of high visual extinction. The western half of Perseus contains three-quarters of the total number of embedded YSOs (Class 0+I and Flat SED sources) in the cloud and nearly as many embedded YSOs as Class II and III sources. Class II and III greatly outnumber Class 0+I objects in eastern Perseus and IC 348. These results are consistent with previous age estimates for the clusters. Across the cloud, 56% of YSOs and 91% of the Class 0+I and Flat sources are in areas where Av > 5 mag, indicating a possible extinction threshold for star formation.

We present 107 maps of continuum emission at 350 microns from Galactic molecular clumps. Observed sources were mainly selected from the Bolocam Galactic Plane Survey (BGPS) catalog, with 3 additional maps covering star forming regions in the outer Galaxy. The higher resolution of the SHARC-II images (8.5" beam) compared with the 1.1 mm images from BGPS (33" beam) allowed us to identify a large population of smaller substructures within the clumps. A catalog is presented for the 1386 sources extracted from the 350 micron maps. The color temperature distribution of clumps based on the two wavelengths has a median of 13.3 K and mean of 16.3 +- 0.4 K, assuming an opacity law index of 1.7. For the structures with good determination of color temperatures, the mean ratio of gas temperature, determined from NH3 observations, to dust color temperature is 0.88 and the median ratio is 0.76. About half the clumps have more than two substructures and 22 clumps have more than 10. The fraction of the mass in dense substructures seen at 350 microns compared to the mass of their parental clump is ~0.19, and the surface densities of these substructures are, on average, 2.2 times those seen in the clumps identified at 1.1 mm. For a well-characterized sample, 88 structures (31%) exceed a surface density of 0.2 g cm^(-2), and 18 (6%) exceed 1.0 g cm^(-2), thresholds for massive star formation suggested by theorists.

We present a multi-wavelength observational study of a low-mass star-forming region, L1251-C, with observational results at wavelengths from the near-infrared to the millimeter. Spitzer Space Telescope observations confirmed that IRAS 22343+7501 is a small group of protostellar objects. The extended emission to east-west direction with its intensity peak at the center of L1251A has been detected at 350 and 850 micron with the CSO and JCMT telescopes, tracing dense envelope materials around L1251A. The single-dish data from the KVN and TRAO telescopes show inconsistencies between the intensity peaks of several molecular line emission and that of the continuum emission, suggesting complex distributions of molecular abundances around L1251A. The SMA interferometer data, however, show intensity peaks of CO 2-1 and 13CO 2-1 located at the position of IRS 1, which is both the brightest source in IRAC image and the weakest source in the 1.3 mm dust continuum map. IRS 1 is the strongest candidate for the driving source of the newly detected compact CO 2-1 outflow. Over the whole region (14' by 14') of L125l-C, 3 Class I and 16 Class II sources have been detected, including three YSOs in L1251A. A comparison with the average projected distance among 19 YSOs in L1251-C and that among 3 YSOs in L1251A suggests L1251-C is an example of low-mass cluster formation, where protostellar objects are forming in a small group.

Molecular outflows driven by protostellar cluster members likely impact their surroundings and contribute to turbulence, affecting subsequent star formation. The very young Serpens South cluster consists of a particularly high density and fraction of protostars, yielding a relevant case study for protostellar outflows and their impact on the cluster environment. We combined CO $J=1-0$ observations of this region using the Combined Array for Research in Millimeter-wave Astronomy (CARMA) and the Institut de Radioastronomie Millim\'{e}trique (IRAM) 30 m single dish telescope. The combined map allows us to probe CO outflows within the central, most active region at size scales of 0.01 pc to 0.8 pc. We account for effects of line opacity and excitation temperature variations by incorporating $^{12}$CO and $^{13}$CO data for the $J=1-0$ and $J=3-2$ transitions (using Atacama Pathfinder Experiment and Caltech Submillimeter Observatory observations for the higher CO transitions), and we calculate mass, momentum, and energy of the molecular outflows in this region. The outflow mass loss rate, force, and luminosity, compared with diagnostics of turbulence and gravity, suggest that outflows drive a sufficient amount of energy to sustain turbulence, but not enough energy to substantially counter the gravitational potential energy and disrupt the clump. Further, we compare Serpens South with the slightly more evolved cluster NGC 1333, and we propose an empirical scenario for outflow-cluster interaction at different evolutionary stages.

We present SCUBA-2 450micron and 850micron observations of the Serpens MWC 297 region, part of the JCMT Gould Belt Survey of nearby star-forming regions. Simulations suggest that radiative feedback influences the star-formation process and we investigate observational evidence for this by constructing temperature maps. Maps are derived from the ratio of SCUBA-2 fluxes and a two component model of the JCMT beam for a fixed dust opacity spectral index of beta = 1.8. Within 40 of the B1.5Ve Herbig star MWC 297, the submillimetre fluxes are contaminated by free-free emission with a spectral index of 1.03+-0.02, consistent with an ultra-compact HII region and polar winds/jets. Contamination accounts for 73+-5 per cent and 82+-4 per cent of peak flux at 450micron and 850micron respectively. The residual thermal disk of the star is almost undetectable at these wavelengths. Young Stellar Objects are confirmed where SCUBA-2 850micron clumps identified by the fellwalker algorithm coincide with Spitzer Gould Belt Survey detections. We identify 23 objects and use Tbol to classify nine YSOs with masses 0.09 to 5.1 Msun. We find two Class 0, one Class 0/I, three Class I and three Class II sources. The mean temperature is 15+-2K for the nine YSOs and 32+-4K for the 14 starless clumps. We observe a starless clump with an abnormally high mean temperature of 46+-2K and conclude that it is radiatively heated by the star MWC 297. Jeans stability provides evidence that radiative heating by the star MWC 297 may be suppressing clump collapse.

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 couple non-magnetic, hydrodynamical simulations of collapsing protostellar cores with radiative transfer evolutionary models to generate synthetic observations. We then use these synthetic observations to investigate the extent to which a simple method for measuring protostellar disc masses used in the literature recovers the intrinsic masses of the discs formed in the simulations. We evaluate the effects of contamination from the surrounding core, partially resolving out the disc, optical depth, fixed assumed dust temperatures, inclination, and the dust opacity law. We show that the combination of these effects can lead to disc mass underestimates by up to factors of 2-3 at millimeter wavelengths and up to an order of magnitude or larger at submillimeter wavelengths. The optically thin portions of protostellar discs are generally cooler in the Class I stage than the Class 0 stage since Class I discs are typically larger and more optically thick, and thus more shielded. The observed disc mass distribution closely resembles the intrinsic distribution if this effect is taken into account, especially at millimeter wavelengths where optical depth effects are minimized. Approximately 50%-70% of protostellar discs observed to date with this method are consistent with the masses of the gravitationally unstable discs formed in the simulations, suggesting that at least some protostellar discs are likely sufficiently massive to fragment. We emphasize key future work needed to confirm these results, including assembling larger, less biased samples, and using molecular line observations to distinguish between rotationally supported, Keplerian discs and magnetically supported pseudodiscs.

We explore the relationship between gas and dust in massive star-forming regions by comparing physical properties derived from each. We compare the temperatures and column densities in a massive star-forming Infrared Dark Cloud (IRDC, G32.02+0.05), which shows a range of evolutionary states, from quiescent to active. The gas properties were derived using radiative transfer modeling of the (1,1), (2,2), and (4,4) transitions of NH3 on the Karl G. Jansky Very Large Array (VLA), while the dust temperatures and column densities were calculated using cirrus-subtracted, modified blackbody fits to Herschel data. We compare the derived column densities to calculate an NH3 abundance, 4.6 x 10^-8. In the coldest star-forming region, we find that the measured dust temperatures are lower than the measured gas temperatures (mean and standard deviations T_dust ~ 11.6 +/- 0.2 K vs. T_gas ~ 15.2 +/- 1.5 K), which may indicate that the gas and dust are not well-coupled in the youngest regions (~0.5 Myr) or that these observations probe a regime where the dust and/or gas temperature measurements are unreliable. Finally, we calculate millimeter fluxes based on the temperatures and column densities derived from NH3 which suggest that millimeter dust continuum observations of massive star-forming regions, such as the Bolocam Galactic Plane Survey or ATLASGAL, can probe hot cores, cold cores, and the dense gas lanes from which they form, and are generally not dominated by the hottest core.

We present new NH3 (1,1), (2,2), and (4,4) observations from the Karl G. Jansky Very Large Array (VLA) compiled with work in the literature to explore the range of conditions observed in young, massive star-forming regions. To sample the effects of evolution independent from those of distance/resolution, abundance, and large-scale environment, we compare clumps in different evolutionary stages within a single Infrared Dark Cloud (IRDC), G32.02+0.06. We find that the early stages of clustered star formation are characterized by dense, parsec-scale filamentary structures interspersed with complexes of dense cores (<0.1 pc cores clustered in complexes separated by ~1 pc) with masses from about 10 to 100 Msun. The most quiescent core is the most extended while the star-forming cores are denser and more compact, showing very similar column density structure before and shortly after the onset of massive star formation, with peak surface densities >~ 1 g cm^-2. Quiescent cores and filaments show smoothly varying temperatures from 10-20 K, rising to over 40 K in star-forming cores. We calculate the virial parameters for 16 cores and find that the level of support provided by turbulence is generally insufficient to support them against gravitational collapse (alpha_vir ~ 0.6). The star-forming filaments show smooth velocity fields, punctuated by discontinuities at the sites of active star formation. We discuss the Massive Molecular Filament (MMF; M > 10^5 Msun, l > 60 pc) hosting the IRDC, hypothesizing that it may have been shaped by previous generations of massive stars.

We present observations of the Auriga-California Molecular Cloud (AMC) at 3.6, 4.5, 5.8, 8.0, 24, 70 and 160 micron observed with the IRAC and MIPS detectors as part of the Spitzer Gould Belt Legacy Survey. The total mapped areas are 2.5 sq-deg with IRAC and 10.47 sq-deg with MIPS. This giant molecular cloud is one of two in the nearby Gould Belt of star-forming regions, the other being the Orion A Molecular Cloud (OMC). We compare source counts, colors and magnitudes in our observed region to a subset of the SWIRE data that was processed through our pipeline. Using color-magnitude and color-color diagrams, we find evidence for a substantial population of 166 young stellar objects (YSOs) in the cloud, many of which were previously unknown. Most of this population is concentrated around the LkHalpha 101 cluster and the filament extending from it. We present a quantitative description of the degree of clustering and discuss the fraction of YSOs in the region with disks relative to an estimate of the diskless YSO population. Although the AMC is similar in mass, size and distance to the OMC, it is forming about 15 - 20 times fewer stars.

In the last twenty years, the topic of episodic accretion has gained significant interest in the star formation community. It is now viewed as a common, though still poorly understood, phenomenon in low-mass star formation. The FU Orionis objects (FUors) are long-studied examples of this phenomenon. FUors are believed to undergo accretion outbursts during which the accretion rate rapidly increases from typically $10^{-7}$ to a few $10^{-4}$ $M_\odot$ yr$^{-1}$, and remains elevated over several decades or more. EXors, a loosely defined class of pre-main sequence stars, exhibit shorter and repetitive outbursts, associated with lower accretion rates. The relationship between the two classes, and their connection to the standard pre-main sequence evolutionary sequence, is an open question: do they represent two distinct classes, are they triggered by the same physical mechanism, and do they occur in the same evolutionary phases? Over the past couple of decades, many theoretical and numerical models have been developed to explain the origin of FUor and EXor outbursts. In parallel, such accretion bursts have been detected at an increasing rate, and as observing techniques improve each individual outburst is studied in increasing detail. We summarize key observations of pre-main sequence star outbursts, and review the latest thinking on outburst triggering mechanisms, the propagation of outbursts from star/disk to disk/jet systems, the relation between classical EXors and FUors, and newly discovered outbursting sources -- all of which shed new light on episodic accretion. We finally highlight some of the most promising directions for this field in the near- and long-term.

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.

We present a survey of 28 molecular outflows driven by low-mass protostars, all of which are sufficiently isolated spatially and/or kinematically to fully separate into individual outflows. Using a combination of new and archival data from several single-dish telescopes, 17 outflows are mapped in CO (2-1) and 17 are mapped in CO (3-2), with 6 mapped in both transitions. For each outflow, we calculate and tabulate the mass, momentum, kinetic energy, mechanical luminosity, and force assuming optically thin emission in LTE at an excitation temperature of 50 K. We show that all of the calculated properties are underestimated when calculated under these assumptions. Taken together, the effects of opacity, outflow emission at low velocities confused with ambient cloud emission, and emission below the sensitivities of the observations increase outflow masses and dynamical properties by an order of magnitude, on average, and factors of 50-90 in the most extreme cases. Different (and non-uniform) excitation temperatures, inclination effects, and dissociation of molecular gas will all work to further increase outflow properties. Molecular outflows are thus almost certainly more massive and energetic than commonly reported. Additionally, outflow properties are lower, on average, by almost an order of magnitude when calculated from the CO (3-2) maps compared to the CO (2-1) maps, even after accounting for different opacities, map sensitivities, and possible excitation temperature variations. It has recently been argued in the literature that the CO (3-2) line is subthermally excited in outflows, and our results support this finding.

The Bolocam Galactic Plane Survey (BGPS) is a 1.1 mm continuum survey of dense clumps of dust throughout the Galaxy covering 170 square degrees. We present spectroscopic observations using the Heinrich Hertz Submillimeter Telescope of the dense gas tracers, HCO+ and N2H+ 3-2, for all 6194 sources in the Bolocam Galactic Plane Survey v1.0.1 catalog between 7.5 <= l <= 194 degrees. This is the largest targeted spectroscopic survey of dense molecular gas in the Milky Way to date. We find unique velocities for 3126 (50.5%) of the BGPS v1.0.1 sources observed. Strong N2H+ 3-2 emission (T_{mb} > 0.5 K) without HCO+ 3-2 emission does not occur in this catalog. We characterize the properties of the dense molecular gas emission toward the entire sample. HCO+ is very sub-thermally populated and the 3-2 transitions are optically thick toward most BGPS clumps. The median observed line width is 3.3 km/s consistent with supersonic turbulence within BGPS clumps. We find strong correlations between dense molecular gas integrated intensities and 1.1 mm peak flux and the gas kinetic temperature derived from previously published NH3 observations. These intensity correlations are driven by the sensitivity of the 3-2 transitions to excitation conditions rather than by variations in molecular column density or abundance. We identify a subset of 113 sources with stronger N2H+ than HCO+ integrated intensity, but we find no correlations between the N2H+ / HCO+ ratio and 1.1 mm continuum flux density, gas kinetic temperature, or line width. Self-absorbed profiles are rare (1.3%).

We present a re-reduction and expansion of the Bolocam Galactic Plane Survey, first presented by Aguirre et al. (2011) and Rosolowsky et al. (2010). The BGPS is a 1.1 mm survey of dust emission in the Northern galactic plane, covering longitudes -10 < \ell < 90 and latitudes |b| < 0.5 with a typical 1-\sigma RMS sensitivity of 30-100 mJy in a 33" beam. Version 2 of the survey includes an additional 20 square degrees of coverage in the 3rd and 4th quadrants and 2 square degrees in the 1st quadrant. The new data release has improved angular recovery, with complete recovery out to 80" and partial recovery to 300", and reduced negative bowls around bright sources resulting from the atmospheric subtraction process. We resolve the factor of 1.5 flux calibration offset between the v1.0 data release and other data sets and determine that there is no offset between v2.0 and other data sets. The v2.0 pointing accuracy is tested against other surveys and demonstrated to be accurate and an improvement over v1.0. We present simulations and tests of the pipeline and its properties, including measurements of the pipeline's angular transfer function. The Bolocat cataloging tool was used to extract a new catalog, which includes 8594 sources, with 591 in the expanded regions. We have demonstrated that the Bolocat 40" and 80" apertures are accurate even in the presence of strong extended background emission. The number of sources is lower than in v1.0, but the amount of flux and area included in identified sources is larger.

Using numerical hydrodynamic simulations, we study the gravitational fragmentation of an unstable protostellar disc formed during the collapse of a pre-stellar core with a mass of 1.2 M_sun. The forming fragments span a mass range from about a Jupiter mass to very-low-mass protostars and are located at distances from a few tens to a thousand AU, with a dearth of objects at < 100 AU. We explore the possibility of observational detection of the fragments in discs viewed through the outflow cavity at a distance of 250 pc. We demonstrate that one hour of integration time with the Atacama Large Millimeter/sub-millimeter Array (ALMA) is sufficient to detect the fragments with masses as low as 1.5 M_Jup at orbital distances up to 800 AU from the protostar. The ALMA resolution sets the limit on the minimum orbital distance of detectable fragments. For the adopted resolution of our simulated ALMA images of 0.1", the fragments can be detected at distances down to 50 AU. At smaller distances, the fragments usually merge with the central density peak. The likelihood for detecting the fragments reduces significantly for a lower resolution of 0.5". Some of the most massive fragments, regardless of their orbital distance, can produce characteristic peaks at approximately 5 micron and hence their presence can be indirectly inferred from the observed spectral energy distributions of protostars.

We present Herschel-HIFI, SPIRE, and PACS 50-670 {\mu}m imaging and spectroscopy of six FU Orionis-type objects and candidates (FU Orionis, V1735 Cyg, V1515 Cyg, V1057 Cyg, V1331 Cyg, and HBC 722), ranging in outburst date from 1936-2010, from the "FOOSH" (FU Orionis Objects Surveyed with Herschel) program, as well as ancillary results from Spitzer-IRS and the Caltech Submillimeter Observatory. In their system properties (Lbol, Tbol, line emission), we find that FUors are in a variety of evolutionary states. Additionally, some FUors have features of both Class I and II sources: warm continuum consistent with Class II sources, but rotational line emission typical of Class I, far higher than Class II sources of similar mass/luminosity. Combining several classification techniques, we find an evolutionary sequence consistent with previous mid-IR indicators. We detect [O I] in every source at luminosities consistent with Class 0/I protostars, much greater than in Class II disks. We detect transitions of 13CO (J_up of 5 to 8) around two sources (V1735 Cyg and HBC 722) but attribute them to nearby protostars. Of the remaining sources, three (FU Ori, V1515 Cyg, and V1331 Cyg) exhibit only low-lying CO, but one (V1057 Cyg) shows CO up to J = 23 - 22 and evidence for H2O and OH emission, at strengths typical of protostars rather than T Tauri stars. Rotational temperatures for "cool" CO components range from 20-81 K, for ~ 10^50 total CO molecules. We detect [C I] and [N II] primarily as diffuse emission.

We present 50-210 um spectral scans of 30 Class 0/I protostellar sources, obtained with Herschel-PACS, and 0.5-1000 um SEDs, as part of the Dust, Ice, and Gas in Time (DIGIT) Key Program. Some sources exhibit up to 75 H2O lines ranging in excitation energy from 100-2000 K, 12 transitions of OH, and CO rotational lines ranging from J=14-13 up to J=40-39. [O I] is detected in all but one source in the entire sample; among the sources with detectable [O I] are two Very Low Luminosity Objects (VeLLOs). The mean 63/145 um [O I] flux ratio is 17.2 +/- 9.2. The [O I] 63 um line correlates with Lbol, but not with the time-averaged outflow rate derived from low-J CO maps. [C II] emission is in general not local to the source. The sample Lbol increased by 1.25 (1.06) and Tbol decreased to 0.96 (0.96) of mean (median) values with the inclusion of the Herschel data. Most CO rotational diagrams are characterized by two optically thin components (<N> = (0.70 +/- 1.12) x 10^49 total particles). N_CO correlates strongly with Lbol, but neither Trot nor N_CO(warm)/N_CO(hot) correlates with Lbol, suggesting that the total excited gas is related to the current source luminosity, but that the excitation is primarily determined by the physics of the interaction (e.g., UV- heating/shocks). Rotational temperatures for H2O (<Trot> = 194 +/- 85 K) and OH (<Trot> = 183 +/- 117 K) are generally lower than for CO, and much of the scatter in the observations about the best fit is attributed to differences in excitation conditions and optical depths amongst the detected lines.

We present high angular resolution 1.3 mm and 850 um dust continuum data obtained with the Submillimeter Array toward 33 Class 0 protostars in nearby clouds (distance < 500 pc), which represents so far the largest survey toward protostellar binary/multiple systems. The median angular resolution in the survey is 2.5 arcsec, while the median linear resolution is approximately 600 AU. Compact dust continuum emission is observed from all sources in the sample. Twenty-one sources in the sample show signatures of binarity/multiplicity, with separations ranging from 50 to 5000 AU. The numbers of singles, binaries, triples, and quadruples in the sample are 12, 14, 5, and 2, respectively. The derived multiplicity frequency (MF) and companion star fraction (CSF) for Class 0 protostars are 0.64+/-0.08 and 0.91+/-0.05, respectively, with no correction for completeness. The derived MF and CSF in this survey are approximately two times higher than the values found in the binary surveys toward Class I YSOs, and approximately three (for MF) and four (for CSF) times larger than the values found among MS stars, with a similar range of separations. Furthermore, the observed fraction of high order multiple systems to binary systems in Class 0 protostars (0.50+/-0.09) is also larger than the fractions found in Class I YSOs (0.31+/-0.07) and MS stars (< 0.2). These results suggest that binary properties evolve as protostars evolve, as predicted by numerical simulations. The distribution of separations for Class 0 protostellar binary systems shows a general trend in which companion star fraction increases with decreasing companion separation. We find that 67%+/-8% of the protobinary systems have circumstellar mass ratios below 0.5, implying that unequal-mass systems are preferred in the process of binary star formation. We suggest an empirical sequential fragmentation picture for binary star formation.

Motivated by the long-standing "luminosity problem" in low-mass star formation whereby protostars are underluminous compared to theoretical expectations, we identify 230 protostars in 18 molecular clouds observed by two Spitzer Space Telescope Legacy surveys of nearby star-forming regions. We compile complete spectral energy distributions, calculate Lbol for each source, and study the protostellar luminosity distribution. This distribution extends over three orders of magnitude, from 0.01 Lsun - 69 Lsun, and has a mean and median of 4.3 Lsun and 1.3 Lsun, respectively. The distributions are very similar for Class 0 and Class I sources except for an excess of low luminosity (Lbol < 0.5 Lsun) Class I sources compared to Class 0. 100 out of the 230 protostars (43%) lack any available data in the far-infrared and submillimeter (70 um < wavelength < 850 um) and have Lbol underestimated by factors of 2.5 on average, and up to factors of 8-10 in extreme cases. Correcting these underestimates for each source individually once additional data becomes available will likely increase both the mean and median of the sample by 35% - 40%. We discuss and compare our results to several recent theoretical studies of protostellar luminosities and show that our new results do not invalidate the conclusions of any of these studies. As these studies demonstrate that there is more than one plausible accretion scenario that can match observations, future attention is clearly needed. The better statistics provided by our increased dataset should aid such future work.

We present a systematic single-dish search for molecular outflows toward a sample of 9 candidate low-luminosity protostars and 30 candidate Very Low Luminosity Objects (VeLLOs; L_int < 0.1 L_sun). The sources are identified using data from the Spitzer Space Telescope catalogued by Dunham et al. toward nearby (D < 400 pc) star forming regions. Each object was observed in 12CO and 13CO J = 2-1 simultaneously using the sideband separating ALMA Band-6 prototype receiver on the Heinrich Hertz Telescope at 30 arcsecond resolution. Using 5-point grid maps we identify five new potential outflow candidates and make on-the-fly maps of the regions surrounding sources in the dense cores B59, L1148, L1228, and L1165. Of these new outflow candidates, only the map of B59 shows a candidate blue outflow lobe associated with a source in our survey. We also present larger and more sensitive maps of the previously detected L673-7 and the L1251-A IRS4 outflows and analyze their properties in comparison to other outflows from VeLLOs. The accretion luminosities derived from the outflow properties of the VeLLOs with detected CO outflows are higher than the observed internal luminosity of the protostars, indicating that these sources likely had higher accretion rates in the past. The known L1251-A IRS3 outflow is detected but not remapped. We do not detect clear, unconfused signatures of red and blue molecular wings toward the other 31 sources in the survey indicating that large-scale, distinct outflows are rare toward this sample of candidate protostars. Several potential outflows are confused with kinematic structure in the surrounding core and cloud. Interferometric imaging is needed to disentangle large-scale molecular cloud kinematics from these potentially weak protostellar outflows.

We present Spitzer IRS spectroscopy of CO2 ice bending mode spectra at 15.2 micrometer toward 19 young stellar objects with luminosity lower than 1 Lsun (3 with luminosity lower than 0.1 Lsun). Ice on dust grain surfaces can encode the history of heating because pure CO2 ice forms only at elevated temperature, T > 20 K, and thus around protostars of higher luminosity. Current internal luminosities of YSOs with L < 1 Lsun do not provide the conditions needed to produce pure CO2 ice at radii where typical envelopes begin. The presence of detectable amounts of pure CO2 ice would signify a higher past luminosity. Many of the spectra require a contribution from a pure, crystalline CO2 component, traced by the presence of a characteristic band splitting in the 15.2 micrometer bending mode. About half of the sources (9 out of 19) in the low luminosity sample have evidence for pure CO2 ice, and six of these have significant double-peaked features, which are very strong evidence of pure CO2 ice. The presence of the pure CO2 ice component indicates that the dust temperature, and hence luminosity of the central star/accretion disk system, must have been higher in the past. An episodic accretion scenario, in which mixed CO-CO2 ice is converted to pure CO2 ice during each high luminosity phase, explains the presence of pure CO2 ice, the total amount of CO2 ice, and the observed residual C18O gas.

We present Spitzer IRAC (2.1 sq. deg.) and MIPS (6.5 sq. deg.) observations of star formation in the Ophiuchus North molecular clouds. This fragmentary cloud complex lies on the edge of the Sco-Cen OB association, several degrees to the north of the well-known rho Oph star-forming region, at an approximate distance of 130 pc. The Ophiuchus North clouds were mapped as part of the Spitzer Gould Belt project under the working name `Scorpius'. In the regions mapped, selected to encompass all the cloud with visual extinction AV>3, eleven Young Stellar Object (YSO) candidates are identified, eight from IRAC/MIPS colour-based selection and three from 2MASS K/MIPS colours. Adding to one source previously identified in L43 (Chen et al. 2009), this increases the number of YSOcs identified in Oph N to twelve. During the selection process, four colour-based YSO candidates were rejected as probable AGB stars and one as a known galaxy. The sources span the full range of YSOc classifications from Class 0/1 to Class III, and starless cores are also present. Twelve high-extinction (AV>10) cores are identified with a total mass of approx. 100 solar masses. These results confirm that there is little ongoing star formation in this region (instantaneous star formation efficiency <0.34%) and that the bottleneck lies in the formation of dense cores. The influence of the nearby Upper Sco OB association, including the 09V star zeta Oph, is seen in dynamical interactions and raised dust temperatures but has not enhanced levels of star formation in Ophiuchus North.

We present 230 GHz Submillimeter Array continuum and molecular line observations of the newly discovered FUor candidate HBC722. We report the detection of seven 1.3 mm continuum sources in the vicinity of HBC722, none of which correspond to HBC722 itself. We compile infrared and submillimeter continuum photometry of each source from previous studies and conclude that three are Class 0 embedded protostars, one is a Class I embedded protostar, one is a Class I/II transition object, and two are either starless cores or very young, very low luminosity protostars or first hydrostatic cores. We detect a northwest-southeast outflow, consistent with the previous detection of such an outflow in low-resolution, single-dish observations, and note that its axis may be precessing. We show that this outflow is centered on and driven by one of the nearby Class 0 sources rather than HBC722, and find no conclusive evidence that HBC722 itself is driving an outflow. The non-detection of HBC722 in the 1.3 mm continuum observations suggests an upper limit of 0.02 solar masses for the mass of the circumstellar disk. This limit is consistent with typical T Tauri disks and with a disk that provides sufficient mass to power the burst.

We present high angular resolution SMA and Spitzer observations toward the Bok globule CB17. SMA 1.3mm dust continuum images reveal within CB17 two sources with an angular separation of about 21" (about 5250 AU at a distance of 250 pc). The northwestern continuum source, referred to as CB17 IRS, dominates the infrared emission in the Spitzer images, drives a bipolar outflow extending in the northwest-southeast direction, and is classified as a low luminosity Class0/I transition object (L_bol ~ 0.5 L_sun). The southeastern continuum source, referred to as CB17 MMS, has faint dust continuum emission in the SMA 1.3mm observations (about 6 sigma detection; ~3.8 mJy), but is not detected in the deep Spitzer infrared images at wavelengths from 3.6 to 70 micron. Its bolometric luminosity and temperature, estimated from its spectral energy distribution, are less than 0.04 L_sun and 16 K, respectively. The SMA CO(2-1) observations suggest that CB17 MMS may drive a low-velocity molecular outflow (about 2.5 km/s), extending in the east-west direction. Comparisons with prestellar cores and Class0 protostars suggest that CB17 MMS is more evolved than prestellar cores but less evolved than Class0 protostars. The observed characteristics of CB17 MMS are consistent with the theoretical predictions from radiative/magneto hydrodynamical simulations of a first hydrostatic core, but there is also the possibility that CB17 MMS is an extremely low luminosity protostar deeply embedded in an edge-on circumstellar disk. Further observations are needed to study the properties of CB17 MMS and to address more precisely its evolutionary stage.

We present high angular resolution observations of the Class 0 protostar IRAM04191+1522, using the Submillimeter Array (SMA). The SMA 1.3 mm continuum images reveal within IRAM04191+1522 two distinct sources with an angular separation of 7.8\,$\pm$\,0.2$"$. The two continuum sources are located in the southeast-northwest direction, with total gas masses of about 0.011 M_sun and about 0.005 M_sun, respectively. The southeastern source, associated with an infrared source seen in the Spitzer images, is the well-known Class 0 protostar with a bolometric luminosity of about 0.08 L_sun. The newly-discovered northwestern continuum source is not visible in the Spitzer images at wavelengths from 3.6 to 70 micron, and has an extremely low bolometric luminosity (< 0.03 L_sun). Complementary IRAM N2H+(1-0) data that probe the dense gas in the common envelope suggest that the two sources were formed through the rotational fragmentation of an elongated dense core. Furthermore, comparisons between IRAM04191+1522 and other protostars suggest that most cores with binary systems formed therein have ratios of rotational energy to gravitational energy $\beta_{\rm rot}$ > 1%. This is consistent with theoretical simulations and indicates that the level of rotational energy in a dense core plays an important role in the fragmentation process.

We determine the observational signatures of protostellar cores by coupling two-dimensional radiative transfer calculations with numerical hydrodynamical simulations that predict accretion rates that both decline with time and feature short-term variability and episodic bursts caused by disk gravitational instability and fragmentation. We calculate the radiative transfer of the collapsing cores throughout the full duration of the collapse, using as inputs the core, disk, and protostellar masses, radii, and mass accretion rates predicted by the hydrodynamical simulations. From the resulting spectral energy distributions, we calculate standard observational signatures (bolometric luminosity, bolometric temperature, ratio of bolometric to submillimeter luminosity) to directly compare to observations. We show that the accretion process predicted by these models reproduces the full spread of observed protostars in both Lbol - Tbol and Lbol - core mass space, including very low luminosity objects, provides a reasonable match to the observed protostellar luminosity distribution, and resolves the long-standing luminosity problem. These models predict an embedded phase duration shorter than recent observationally determined estimates (0.12 Myr vs. 0.44 Myr), and a fraction of total time spent in Stage 0 of 23%, consistent with the range of values determined by observations. On average, the models spend 1.3% of their total time in accretion bursts, during which 5.3% of the final stellar mass accretes, with maximum values being 11.8% and 35.5% for the total time and accreted stellar mass, respectively. Time-averaged models that filter out the accretion variability and bursts do not provide as good of a match to the observed luminosity problem, suggesting that the bursts are required.

We present the results of a GBT survey of NH3(1,1), (2,2), (3,3) lines towards 631 Bolocam Galactic Plane Survey (BGPS) sources at a range of Galactic longitudes in the inner Galaxy. We have detected the NH3(1,1) line towards 72% of our targets (456), demonstrating that the high column density features identified in the BGPS and other continuum surveys accurately predict the presence of dense gas. We have determined kinematic distances and resolved the distance ambiguity for all BGPS sources detected in NH3. The BGPS sources trace the locations of the Scutum and Sagittarius spiral arms, with the number of sources peaking between Galactocentric radii of 4-5 kpc. We measure the physical properties of each source and find that depending on the distance, BGPS sources are primarily clumps, with some cores and clouds. We have examined the physical properties as a function of Galactocentric distance, and find a mean gas kinetic temperature of 15.6 K, and that the NH3 column density and abundance decrease by nearly an order of magnitude between Galactocentric radii of 3 - 11 kpc. Comparing sources at similar distances demonstrates that the physical properties are indistinguishable, which suggests a similarity in clump structure across the Galactic disk. We have also compared the BGPS sources to criteria for efficient star formation presented independently by Heiderman et al. and Lada et al., and massive star formation presented by Kauffmann et al. 48% of our sample should be forming stars (including massive stars) with high efficiency, and 87% contain subregions that should be efficiently forming stars. Indeed, we find that 67% of the sample exhibit signs of star formation activity based on an association with a mid-IR source.

We present new 230 GHz Submillimeter Array observations of the candidate first hydrostatic core Per-Bolo 58. We report the detection of a 1.3 mm continuum source and a bipolar molecular outflow, both centered on the position of the candidate first hydrostatic core. The continuum detection has a total flux density of 26.6 +/- 4.0 mJy, from which we calculate a total (gas and dust) mass of 0.11 +/- 0.05 Msun and a mean number density of 2.0 +/- 1.6 X 10^7 cm-3. There is some evidence for the existence of an unresolved component in the continuum detection, but longer-baseline observations are required in order to confirm the presence of this component and determine whether its origin lies in a circumstellar disk or in the dense inner envelope. The bipolar molecular outflow is observed along a nearly due east-west axis. The outflow is slow (characteristic velocity of 2.9 km/s), shows a jet-like morphology (opening semi-angles ~8 degrees for both lobes), and extends to the edges of the primary beam. We calculate the kinematic and dynamic properties of the outflow in the standard manner and compare them to several other protostars and candidate first hydrostatic cores with similarly low luminosities. We discuss the evidence both in support of and against the possibility that Per-Bolo 58 is a first hydrostatic core, and we outline future work needed to further evaluate the evolutionary status of this object.

The Red MSX Source (RMS) survey has identified a sample of ~1200 massive young stellar objects (MYSOs), compact and ultra compact HII regions from a sample of ~2000 MSX and 2MASS colour selected sources. We have used the 100 m Green Bank telescope to search for 22-24 GHz water maser and ammonia (1,1), (2,2) and (3,3) emission towards ~600 RMS sources located within the northern Galactic plane. We have identified 308 H2O masers which corresponds to an overall detection rate of ~50%. Abridged: We detect ammonia emission towards 479 of these massive young stars, which corresponds to ~80%. Ammonia is an excellent probe of high density gas allowing us to measure key parameters such as gas temperatures, opacities, and column densities, as well as providing an insight into the gas kinematics. The average kinetic temperature, FWHM line width and total NH3 column density for the sample are approximately 22 K, 2 km/s and 2x10^{15} cm^{-2}, respectively. We find that the NH3 (1,1) line width and kinetic temperature are correlated with luminosity and finding no underlying dependence of these parameters on the evolutionary phase of the embedded sources, we conclude that the observed trends in the derived parameters are more likely to be due to the energy output of the central source and/or the line width-clump mass relationship. The velocities of the peak H2O masers and the NH3 emission are in excellent agreement with each other, which would strongly suggest an association between the dense gas and the maser emission. Moreover, we find the bolometric luminosity of the embedded source and the isotropic luminosity of the H2O maser are also correlated. We conclude from the correlations of the cloud and water maser velocities and the bolometric and maser luminosity that there is a strong dynamical relationship between the embedded young massive star and the H2O maser.

We present the results of observations of 1882 sources in the Bolocam Galactic Plane Survey (BGPS) at 1.1 mm with the 10m Heinrich Hertz Telescope simultaneously in HCO+ J=3-2 and N2H+ J=3-2. We detect 77% of these sources in HCO^+ and 51% in N2H+ at greater than 3$\sigma$. We find a strong correlation between the integrated intensity of both dense gas tracers and the 1.1 mm dust emission of BGPS sources. We determine kinematic distances for 529 sources (440 in the first quadrant breaking the distance ambiguity and 89 in the second quadrant) We derive the size, mass, and average density for this subset of clumps. The median size of BGPS clumps is 0.75 pc with a median mass of 330 M$_{\odot}$ (assuming T_{Dust}=20 K). The median HCO+ linewidth is 2.9 km s$^{-1}$ indicating that BGPS clumps are dominated by supersonic turbulence or unresolved kinematic motions. We find no evidence for a size-linewidth relationship for BGPS clumps. We analyze the effects of the assumed dust temperature on the derived clump properties with a Monte Carlo simulation and we find that changing the temperature distribution will change the median source properties (mass, volume-averaged number density, surface density) by factors of a few. The observed differential mass distribution has a power-law slope that is intermediate between that observed for diffuse CO clouds and the stellar IMF. BGPS clumps represent a wide range of objects (from dense cores to more diffuse clumps) and are typically characterized by larger sizes and lower densities than previously published surveys of high-mass star-forming regions. This collection of objects is a less-biased sample of star-forming regions in the Milky Way that likely span a wide range of evolutionary states.

We present Spitzer Space Telescope IRAC and MIPS observations of a 0.85 deg^2 field including the Corona Australis (CrA) star-forming region. At a distance of 130 pc, CrA is one of the closest regions known to be actively forming stars, particularly within its embedded association, the Coronet. Using the Spitzer data, we identify 51 young stellar objects (YSOs) in CrA which include sources in the well-studied Coronet cluster as well as distributed throughout the molecular cloud. Twelve of the YSOs discussed are new candidates, one of which is located in the Coronet. Known YSOs retrieved from the literature are also added to the list, and a total of 116 candidate YSOs in CrA are compiled. Based on these YSO candidates, the star formation rate is computed to be 12 M_o Myr^-1, similar to that of the Lupus clouds. A clustering analysis was also performed, finding that the main cluster core, consisting of 68 members, is elongated (having an aspect ratio of 2.36), with a circular radius of 0.59 pc and mean surface density of 150 pc^-2. In addition, we analyze outflows and jets in CrA by means of new CO and H_2 data. We present 1.3 mm interferometric continuum observations made with the Submillimeter Array (SMA) covering R CrA, IRS 5, IRS 7, and IRAS 18595-3712 (IRAS 32). We also present multi-epoch H_2 maps and detect jets and outflows, study their proper motions, and identify exciting sources. The Spitzer and ISAAC/VLT observations of IRAS 32 show a bipolar precessing jet, which drives a CO (2-1) outflow detected in the SMA observations. There is also clear evidence for a parsec-scale precessing outflow, E-W oriented, and originating in the SMA 2 region, likely driven by SMA 2 or IRS 7A.

We analyze the submillimeter emission surrounding the new FU Orionis-type object, HBC 722. We present the first epoch of observations of the active environs of HBC 722, with imaging and spectroscopy from PACS, SPIRE, and HIFI aboard the Herschel Space Observatory, as well as CO J= 2-1 and 350 um imaging (SHARC-II) with the Caltech Submillimeter Observatory. The primary source of submillimeter continuum emission in the region -- 2MASS 20581767+4353310 -- is located 16$\arcsec$ south-southeast of the optical flaring source while the optical and near-IR emission is dominated by HBC 722. A bipolar outflow extends over HBC 722; the most likely driver is the submillimeter source. We detect warm (100 K) and hot (246 K) CO emission in the surrounding region, evidence of outflow-driven heating in the vicinity. The region around HBC 722 itself shows little evidence of heating driven by the new outbursting source itself.

We present the results of a search for mid-infrared signs of star formation activity in the 1.1 mm sources in the Bolocam Galactic Plane Survey (BGPS). We have correlated the BGPS catalog with available mid-IR Galactic plane catalogs based on the Spitzer Space Telescope GLIMPSE legacy survey and the Midcourse Space Experiment (MSX) Galactic plane survey. We find that 44% (3,712 of 8,358) of the BGPS sources contain at least one mid-IR source, including 2,457 of 5,067 (49%) within the area where all surveys overlap (10 deg < l < 65 deg). Accounting for chance alignments between the BGPS and mid-IR sources, we conservatively estimate that 20% of the BPGS sources within the area where all surveys overlap show signs of active star formation. We separate the BGPS sources into four groups based on their probability of star formation activity. Extended Green Objects (EGOs) and Red MSX Sources (RMS) make up the highest probability group, while the lowest probability group is comprised of "starless" BGPS sources which were not matched to any mid-IR sources. The mean 1.1 mm flux of each group increases with increasing probability of active star formation. We also find that the "starless" BGPS sources are the most compact, while the sources with the highest probability of star formation activity are on average more extended with large skirts of emission. A subsample of 280 BGPS sources with known distances demonstrates that mass and mean H_2 column density also increase with probability of star formation activity.

We present Gould's Belt (GB) Spitzer IRAC and MIPS observations of the Lupus V and VI clouds and discuss them in combination with near-infrared (2MASS) data. Our observations complement those obtained for other Lupus clouds within the frame of the Spitzer "Core to Disk" (c2d) Legacy Survey. We found 43 Young Stellar Object (YSO) candidates in Lupus V and 45 in Lupus VI, including 2 transition disks, using the standard c2d/GB selection method. None of these sources was classified as a pre-main sequence star from previous optical, near-IR and X-ray surveys. A large majority of these YSO candidates appear to be surrounded by thin disks (Class III; ~79% in Lupus V and ~87% in Lupus VI). These Class III abundances differ significantly from those observed for the other Lupus clouds and c2d/GB surveyed star-forming regions, where objects with optically thick disks (Class II) dominate the young population. We investigate various scenarios that can explain this discrepancy. In particular, we show that disk photo-evaporation due to nearby OB stars is not responsible for the high fraction of Class III objects. The gas surface densities measured for Lupus V and VI lies below the star-formation threshold (AV {\eqsim}8.6 mag), while this is not the case for other Lupus clouds. Thus, few Myrs older age for the YSOs in Lupus V and VI with respect to other Lupus clouds is the most likely explanation of the high fraction of Class III objects in these clouds, while a higher characteristic stellar mass might be a contributing factor. Better constraints on the age and binary fraction of

We present new observations of the CB130 region, composed of three separate cores. Using the \textit{Spitzer Space Telescope} we detected a Class 0 and a Class II object in one of these, CB130-1. The observed photometric data from \textit{Spitzer} and ground-based telescopes are used to establish the physical parameters of the Class 0 object. SED fitting with a radiative transfer model shows that the luminosity of the Class 0 object is 0.14 $-$ 0.16 L$_{\odot}$, which is a low luminosity for a protostellar object. In order to constrain the chemical characteristics of the core having the low luminosity object, we compare our molecular line observations to models of lines including abundance variations. We tested both ad hoc step function abundance models and a series of self-consistent chemical evolution models. In the chemical evolution models, we consider a continuous accretion model and an episodic accretion model to explore how variable luminosity affects the chemistry. The step function abundance models can match observed lines reasonably well. The best fitting chemical evolution model requires episodic accretion and the formation of CO$_2$ ice from CO ice during the low luminosity periods. This process removes C from the gas phase, providing a much improved fit to the observed gas-phase molecular lines and the CO$_2$ ice absorption feature. Based on the chemical model result, the low luminosity of CB130-1 is explained better as a quiescent stage between episodic accretion bursts rather than being at the first hydrostatic core stage.

The Bolocam Galactic Plane Survey (BGPS) data for a six square degree region of the Galactic plane containing the Galactic center is analyzed and compared to infrared and radio continuum data. The BGPS 1.1 mm emission consists of clumps interconnected by a network of fainter filaments surrounding cavities, a few of which are filled with diffuse near-IR emission indicating the presence of warm dust or with radio continuum characteristic of HII regions or supernova remnants. New 350 {\mu}m images of the environments of the two brightest regions, Sgr A and B, are presented. Sgr B2 is the brightest mm-emitting clump in the Central Molecular Zone and may be forming the closest analog to a super star cluster in the Galaxy. The Central Molecular Zone (CMZ) contains the highest concentration of mm and sub-mm emitting dense clumps in the Galaxy. Most 1.1 mm features at positive longitudes are seen in silhouette against the 3.6 to 24 {\mu}m background observed by the Spitzer Space Telescope. However, only a few clumps at negative longitudes are seen in absorption, confirming the hypothesis that positive longitude clumps in the CMZ tend to be on the near-side of the Galactic center, consistent with the suspected orientation of the central bar in our Galaxy. Some 1.1 mm cloud surfaces are seen in emission at 8 {\mu}m, presumably due to polycyclic aromatic hydrocarbons (PAHs). A ~0.2\degree (~30 pc) diameter cavity and infrared bubble between l \approx 0.0\degree and 0.2\degree surrounds the Arches and Quintuplet clusters and Sgr A. The bubble contains several clumpy dust filaments that point toward Sgr A\ast; its potential role in their formation is explored. [abstract truncated]

We present the Bolocam Galactic Plane Survey (BGPS), a 1.1 mm continuum survey at 33" effective resolution of 170 square degrees of the Galactic Plane visible from the northern hemisphere. The survey is contiguous over the range -10.5 < l < 90.5, |b| < 0.5 and encompasses 133 square degrees, including some extended regions |b| < 1.5. In addition to the contiguous region, four targeted regions in the outer Galaxy were observed: IC1396, a region towards the Perseus Arm, W3/4/5, and Gem OB1. The BGPS has detected approximately 8400 clumps over the entire area to a limiting non-uniform 1-sigma noise level in the range 11 to 53 mJy/beam in the inner Galaxy. The BGPS source catalog is presented in a companion paper (Rosolowsky et al. 2010). This paper details the survey observations and data reduction methods for the images. We discuss in detail the determination of astrometric and flux density calibration uncertainties and compare our results to the literature. Data processing algorithms that separate astronomical signals from time-variable atmospheric fluctuations in the data time-stream are presented. These algorithms reproduce the structure of the astronomical sky over a limited range of angular scales and produce artifacts in the vicinity of bright sources. Based on simulations, we find that extended emission on scales larger than about 5.9' is nearly completely attenuated (> 90%) and the linear scale at which the attenuation reaches 50% is 3.8'. Comparison with other millimeter-wave data sets implies a possible systematic offset in flux calibration, for which no cause has been discovered. This presentation serves as a companion and guide to the public data release through NASA's Infrared Processing and Analysis Center (IPAC) Infrared Science Archive (IRSA). New data releases will be provided through IPAC IRSA with any future improvements in the reduction.

The first hydrostatic core (FHSC) represents a very early phase in the low-mass star formation process, after collapse of the parent core has begun but before a true protostar has formed. This large (few AU), cool (100 K), pressure supported core of molecular hydrogen is expected from theory, but has yet to be observationally verified. Here we present observations of an excellent candidate for the FHSC phase: Per-Bolo 58, a dense core in Perseus that was previously believed to be starless. The 70 micron flux of 65 mJy, from new deep Spitzer MIPS observations, is consistent with that expected for the FHSC. A low signal-to-noise detection at 24 micron leaves open the possibility that Per-Bolo 58 could be a very low luminosity protostar, however. We utilize radiative transfer models to determine the best-fitting FHSC and protostar models to the spectral energy distribution and 2.9 mm visibilities of Per-Bolo 58. The source is consistent with a FHSC with some source of lower opacity through the envelope allowing 24 micron emission to escape; a small outflow cavity or a cavity in the envelope are both possible. While we are unable to rule out the presence of a protostar, if present it would be one of the lowest luminosity protostellar objects yet observed, with an internal luminosity of approximately 0.01 Lsun.

We present new infrared, submillimeter, and millimeter observations of the dense core L673-7 and report the discovery of a low-luminosity, embedded Class 0 protostar driving a molecular outflow. L673-7 is seen in absorption against the mid-infrared background in 5.8, 8, and 24 micron Spitzer images, allowing for a derivation of the column density profile and total enclosed mass of L673-7, independent of dust temperature assumptions. Estimates of the core mass from these absorption profiles range from 0.2-4.5 solar masses. Millimeter continuum emission indicates a mass of about 2 solar masses, both from a direct calculation assuming isothermal dust and from dust radiative transfer models constrained by the millimeter observations. We use dust radiative transfer models to constrain the internal luminosity of L673-7, defined to be the luminosity of the central source and excluding the luminosity from external heating, to be 0.01-0.045 solar luminosities, with 0.04 solar luminosities the most likely value. L673-7 is thus classified as a very low luminosity object (VeLLO), and is among the lowest luminosity VeLLOs yet studied. We calculate the kinematic and dynamic properties of the molecular outflow in the standard manner, and we show that the expected accretion luminosity based on these outflow properties is greater than or equal to 0.36 solar luminosities. The discrepancy between this expected accretion luminosity and the internal luminosity derived from dust radiative transfer models indicates that the current accretion rate is much lower than the average rate over the lifetime of the outflow. Although the protostar embedded within L673-7 is consistent with currently being substellar, it is unlikely to remain as such given the substantial mass reservoir remaining in the core.

We present the 1.1 millimeter Bolocam Galactic Plane Survey (BGPS) observations of the Gemini OB1 molecular cloud complex, and targeted ammonia observations of the BGPS sources. When paired with molecular spectroscopy of a dense gas tracer, millimeter observations yield physical properties such as masses, radii, mean densities, kinetic temperatures and line widths. We detect 34 distinct BGPS sources above 5-sigma=0.37 Jy/beam with corresponding 5-sigma detections in the ammonia (1,1) transition. Eight of the objects show water maser emission (20%). We find a mean millimeter source FWHM of 1.12 pc, and a mean kinetic temperature of 20 K for the sample of 34 BGPS sources. The observed ammonia line widths are dominated by non-thermal motions, typically found to be a few times the thermal sound speed expected for the derived kinetic temperature. We calculate the mass for each source from the millimeter flux assuming the sources are isothermal and find a mean isothermal mass within a 120" aperture of 230 +/- 180 solar masses. We find a total mass of 8,400 solar masses for all BGPS sources in the Gemini OB1 molecular cloud, representing 6.5% of the cloud mass. By comparing the millimeter isothermal mass to the virial mass within a radius equal to the mm source size calculated from the ammonia line widths, we find a mean virial parameter (M_vir/M_iso) of 1.0 +/- 0.9 for the sample. We find mean values for the distributions of column densities of 10^22 cm^-2 for H_2, and 3.0x10^14 cm^-2 for ammonia, giving a mean ammonia abundance of 3.0x10^-8 relative to H_2. We find volume-averaged densities on the order of 10^3-10^4 cm^-3. The sizes and densities suggest that in the Gem OB1 region the BGPS is detecting the clumps from which stellar clusters form, rather than smaller, higher density cores where single stars or small multiple systems form.

We present far-infrared spectroscopic observations, taken with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory, of the protoplanetary disk around the pre-main-sequence star HD 100546. These observations are the first within the DIGIT Herschel key program, which aims to follow the evolution of dust, ice, and gas from young stellar objects still embedded in their parental molecular cloud core, through the final pre-main-sequence phases when the circumstellar disks are dissipated. Our aim is to improve the constraints on temperature and chemical composition of the crystalline olivines in the disk of HD 100546 and to give an inventory of the gas lines present in its far-infrared spectrum. The 69 \mu\m feature is analyzed in terms of position and shape to derive the dust temperature and composition. Furthermore, we detected 32 emission lines from five gaseous species and measured their line fluxes. The 69 \mu\m emission comes either from dust grains with ~70 K at radii larger than 50 AU, as suggested by blackbody fitting, or it arises from ~200 K dust at ~13 AU, close to the midplane, as supported by radiative transfer models. We also conclude that the forsterite crystals have few defects and contain at most a few percent iron by mass. Forbidden line emission from [CII] at 157 \mu\m and [OI] at 63 and 145 \mu\m, most likely due to photodissociation by stellar photons, is detected. Furthermore, five H2O and several OH lines are detected. We also found high-J rotational transition lines of CO, with rotational temperatures of ~300 K for the transitions up to J=22-21 and T~800 K for higher transitions.

DK Cha is an intermediate-mass star in transition from an embedded configuration to a star plus disk stage. We aim to study the composition and energetics of the circumstellar material during this pivotal stage. Using the Range Scan mode of PACS on the Herschel Space Observatory, we obtained a spectrum of DK Cha from 55 to 210 micron as part of the DIGIT Key Program. Almost 50 molecular and atomic lines were detected, many more than the 7 lines detected in ISO-LWS. Nearly the entire ladder of CO from J=14-13 to 38-37 (E_u/k = 4080 K), water from levels as excited as E_u/k = 843 K, and OH lines up to E_u/k = 290 K were detected. The continuum emission in our PACS SED scan matches the flux expected from a model consisting of a star, a surrounding disk of 0.03 Solar mass, and an envelope of a similar mass, supporting the suggestion that the object is emerging from its main accretion stage. Molecular, atomic, and ionic emission lines in the far-infrared reveal the outflow's influence on the envelope. The inferred hot gas can be photon-heated, but some emission could be due to C-shocks in the walls of the outflow cavity.

A long-standing problem in low-mass star formation is the "luminosity problem," whereby protostars are underluminous compared to the accretion luminosity expected both from theoretical collapse calculations and arguments based on the minimum accretion rate necessary to form a star within the embedded phase duration. Motivated by this luminosity problem, we present a set of evolutionary models describing the collapse of low-mass, dense cores into protostars, using the Young & Evans (2005) model as our starting point. We calculate the radiative transfer of the collapsing cores throughout the full duration of the collapse in two dimensions. From the resulting spectral energy distributions, we calculate standard observational signatures to directly compare to observations. We incorporate several modifications and additions to the original Young & Evans model in an effort to better match observations with model predictions. We find that scattering, 2-D geometry, mass-loss, and outflow cavities all affect the model predictions, as expected, but none resolve the luminosity problem. A cycle of episodic mass accretion, however, can resolve this problem and bring the model predictions into better agreement with observations. Standard assumptions about the interplay between mass accretion and mass loss in our model give star formation efficiencies consistent with recent observations that compare the core mass function (CMF) and stellar initial mass function (IMF). The combination of outflow cavities and episodic mass accretion reduce the connection between observational Class and physical Stage to the point where neither of the two common observational signatures (bolometric temperature and ratio of bolometric to submillimeter luminosity) can be considered reliable indicators of physical Stage.

A long infrared jet has been discovered by the Spitzer c2d legacy program in core A of L1251. It is associated with a very embedded Class 0 object with an accretion luminosity of about 0.9 Lsun derived by radiative transfer model fitting to the observed SED. Comparing the observed IRAC colors along the infrared jet with those calculated from a model of an admixture of gas with a power-law temperature distribution indicates that the jet is possibly created by a paraboloidal bow shock propagating into the ambient medium of n(H_2)=10^5 cm^{-3}. In addition, the variation of the power-law index along the jet suggests that the portion of hot gas decreases with distance from the jet engine. The molecular outflow in this region has been mapped for the first time using CO data. From the calculated outflow momentum flux, a very strong lower limit to the average accretion luminosity is 3.6 sin i/cos^3 i Lsun, indicative of a decrease in the accretion rate with time.

We report observations of 15 high redshift (z = 1-5) galaxies at 350 microns using the Caltech Submillimeter Observatory and SHARC-II array detector. Emission was detected from eight galaxies, for which far-infrared luminosities, star formation rates, total dust masses, and minimum source size estimates are derived. These galaxies have star formation rates and star formation efficiencies comparable to other high redshift molecular emission line galaxies. The results are used to test the idea that star formation in these galaxies occurs in a large number of basic units, the units being similar to star-forming clumps in the Milky Way. The luminosity of these extreme galaxies can be reproduced in a simple model with (0.9-30) *10^6 dense clumps, each with a luminosity of 5 *10^5 Lsun, the mean value for such clumps in the Milky Way. Radiative transfer models of such clumps can provide reasonable matches to the overall SEDs of the galaxies. They indicate that the individual clumps are quite opaque in the far-infrared. Luminosity to mass ratios vary over two orders of magnitude, correlating strongly with the dust temperature derived from simple fits to the SED. The gas masses derived from the dust modeling are in remarkable agreement with those from CO luminosities, suggesting that the assumptions going into both calculations are reasonable.

We present the first results of a program to characterize the disk and envelope structure of typical Class 0 protostars in nearby low-mass star forming regions. We use Spitzer IRS mid-infrared spectra, high resolution CARMA 230 GHz continuum imaging, and 2-D radiative transfer models to constrain the envelope structure, as well as the size and mass of the circum-protostellar disk in Serpens FIRS 1. The primary envelope parameters (centrifugal radius, outer radius, outflow opening angle, and inclination) are well constrained by the spectral energy distribution (SED), including Spitzer IRAC and MIPS photometry, IRS spectra, and 1.1 mm Bolocam photometry. These together with the excellent uv-coverage (4.5-500 klam) of multiple antenna configurations with CARMA allow for a robust separation of the envelope and a resolved disk. The SED of Serpens FIRS 1 is best fit by an envelope with the density profile of a rotating, collapsing spheroid with an inner (centrifugal) radius of approximately 600 AU, and the millimeter data by a large resolved disk with Mdisk~1.0 Msun and Rdisk~300 AU. These results suggest that large, massive disks can be present early in the main accretion phase. Results for the larger, unbiased sample of Class~0 sources in the Perseus, Serpens, and Ophiuchus molecular clouds are needed to determine if relatively massive disks are typical in the Class 0 stage.

We have conducted a search for ionized gas at 3.6 cm, using the Very Large Array, towards 31 Galactic intermediate- and high-mass clumps detected in previous millimeter continuum observations. In the 10 observed fields, 35 HII regions are identified, of which 20 are newly discovered. Many of the HII regions are multiply peaked indicating the presence of a cluster of massive stars. We find that the ionized gas tends to be associated towards the millimeter clumps; of the 31 millimeter clumps observed, 9 of these appear to be physically related to ionized gas, and a further 6 have ionized gas emission within 1'. For clumps with associated ionized gas, the combined mass of the ionizing massive stars is compared to the clump masses to provide an estimate of the instantaneous star formation efficiency. These values range from a few percent to 25%, and have an average of 7 +/- 8%. We also find a correlation between the clump mass and the mass of the ionizing massive stars within it, which is consistent with a power law. This result is comparable to the prediction of star formation by competitive accretion that a power law relationship exists between the mass of the most massive star in a cluster and the total mass of the remaining stars.

We present a catalog of 8358 sources extracted from images produced by the Bolocam Galactic Plane Survey (BGPS). The BGPS is a survey of the millimeter dust continuum emission from the northern Galactic plane. The catalog sources are extracted using a custom algorithm, Bolocat, which was designed specifically to identify and characterize objects in the large-area maps generated from the Bolocam instrument. The catalog products are designed to facilitate follow-up observations of these relatively unstudied objects. The catalog is 98% complete from 0.4 Jy to 60 Jy over all object sizes for which the survey is sensitive (<3.5'). We find that the sources extracted can best be described as molecular clumps -- large dense regions in molecular clouds linked to cluster formation. We find the flux density distribution of sources follows a power law with dN/dS ~S^(-2.4 +/- 0.1) and that the mean Galactic latitude for sources is significantly below the midplane: <b>=(-0.095 +/- 0.001) deg.

Observations of Lynds Dark Nebula 1221 from the Spitzer Space Telescope are presented. These data show three candidate protostars towards L1221, only two of which were previously known. The infrared observations also show signatures of outflowing material, an interpretation which is also supported by radio observations with the Very Large Array. In addition, molecular line maps from the Five College Radio Astronomy Observatory are shown. One-dimensional dust continuum modelling of two of these protostars, IRS1 and IRS3, is described. These models show two distinctly different protostars forming in very similar environments. IRS1 shows a higher luminosity and larger inner radius of the envelope than IRS3. The disparity could be caused by a difference in age or mass, orientation of outflow cavities, or the impact of a binary in the IRS1 core.

Circumstellar discs are expected to be the nursery of planets. Grain growth within such discs is the first step in the planet formation process. The Bok globule CB 26 harbours such a young disc. We present a detailed model of the edge-on circumstellar disc and its envelope in the Bok globule CB 26. The model is based on HST near-infrared maps in the I, J, H, and K bands, OVRO and SMA radio maps at 1.1mm, 1.3mm and 2.7mm, and the spectral energy distribution (SED) from 0.9 microns to 3mm. New photometric and spectroscopic data from the Spitzer Space Telescope and the Caltech Submilimeter Observatory have been obtained and are part of our analysis. Using the self-consistent radiative transfer code MC3D, the model we construct is able to discriminate parameter sets and dust properties of both its parts, namely envelope and disc. We find that the disc has an inner hole with a radius of 45 +/- 5 AU. Based on a dust model including silicate and graphite the maximum grain size needed to reproduce the spectral millimetre index is 2.5 microns. Features seen in the near-infrared images, dominated by scattered light, can be described as a result of a rotating envelope. Successful employment of ISM dust in both the disc and envelope hint that grain growth may not yet play a significant role for the appearance of this system. A larger inner hole gives rise to the assumption that CB 26 is a circumbinary disc.

Molecular outflows provide an alternative method of identifying protostellar cores, complementary to recent mid-infrared studies. Continuing our studies of Perseus, we investigate whether all Spitzer-identified protostars, and particularly those with low luminosities, drive outflows, and if any new protostellar cores (perhaps harbouring low-mass sources) can be identified via their outflows alone. We have used the heterodyne array receiver HARP on JCMT to make deep 12CO 3-2 maps of submm cores in Perseus, extending and deepening our earlier study with RxB and bringing the total number of SCUBA cores studied up to 83. Our survey includes 23/25 of the Dunham et al. (2008) Spitzer low-luminosity objects believed to be embedded protostars, including three VeLLOs. All but one of the cores identified as harbouring embedded YSOs have outflows, confirming outflow detections as a good method for identifying protostars. We detect outflows from 20 Spitzer low-luminosity objects. We do not conclusively detect any outflows from IR-quiet cores, though confusion in clustered regions such as NGC1333 makes it impossible to identify all the individual driving sources. This similarity in detection rates despite the difference in search methods and detection limits suggests either that the sample of protostars in Perseus is now complete, or that the existence of an outflow contributes to the Spitzer detectability, perhaps through the contribution of shocked H2 emission in the IRAC bands. For five of the low-luminosity sources, there is no protostellar envelope detected at 350 microns and the Spitzer emission is entirely due to shocks. Additionally, we detect the outflow from IRAS 03282+3035 at 850 microns with SCUBA due to CO line contamination in the continuum passband.

We present Spitzer IRAC (~2 deg^2) and MIPS (~8 deg^2) observations of the Cepheus Flare which is associated with the Gould Belt, at an approximate distance of ~300 pc. Around 6500 sources are detected in all four IRAC bands, of which ~900 have MIPS 24 micron detections. We identify 133 YSO candidates using color-magnitude diagram techniques, a large number of the YSO candidates are associated with the NGC 7023 reflection nebula. Cross identifications were made with the Guide Star Catalog II and the IRAS Faint Source Catalog, and spectral energy distributions (SED) were constructed. SED modeling was conducted to estimate the degree of infrared excess. It was found that a large majority of disks were optically thick accreting disks, suggesting that there has been little disk evolution in these sources. Nearest-neighbor clustering analysis identified four small protostellar groups (L1228, L1228N, L1251A, and L1251B) with 5-15 members each and the larger NGC 7023 association with 32 YSO members. The star formation efficiency for cores with clusters of protostars and for those without clusters was found to be ~8% and ~1% respectively. The cores L1155, L1241, and L1247 are confirmed to be starless down to our luminosity limit of L_bol=0.06 L_sol.

This White Paper, submitted to the National Academy of Sciences' Astro2010 Decadal Review Committee, focuses on 2 central themes in the study of young brown dwarfs -- their formation mechanism and disk characteristics -- which are of direct relevance to fundamental questions of stellar and planetary origins and properties.

We report on diffraction-limited observations in the far-infrared and sub- millimeter of the Cluster B region of Serpens (G3-G6 Cluster) and of the Herbig Be star to the south, VV Ser. The observations were made with the Spitzer MIPS instrument in fine-scale mode at 70um, in normal mapping mode at 160um (VV Ser only), and the CSO SHARC-II camera at 350um (Cluster B only). We use these data to define the spectral energy distributions of the tightly grouped members of Cluster B, many of whose SEDs peak in the far-infrared. We compare our results to those of the c2d survey of Serpens and to published models for the far-infrared emission from VV Ser. We find that values of Lbol and Tbol calculated with our new photometry show only modest changes from previous values, and that most source SED classifications remain unchanged.

This paper reports the discovery of a Very Low Luminosity Object (VeLLO) in the "starless" dense core L328, using the Spitzer Space Telescope and ground based observations from near-infrared to millimeter wavelengths. The Spitzer 8 micron image indicates that L328 consists of three subcores of which the smallest one may harbor a source, L328-IRS while two other subcores remain starless. L328-IRS is a Class 0 protostar according to its bolometric temperature (44 K) and the high fraction ~72 % of its luminosity emitted at sub-millimeter wavelengths. Its inferred "internal luminosity" (0.04 - 0.06 Lsun) using a radiative transfer model under the most plausible assumption of its distance as 200 pc is much fainter than for a typical protostar, and even fainter than other VeLLOs studied previously. Note, however, that its inferred luminosity may be uncertain by a factor of 2-3 if we consider two extreme values of the distance of L328-IRS (125 or 310 pc). Low angular resolution observations of CO do not show any clear evidence of a molecular outflow activity. But broad line widths toward L328, and Spitzer and near-infrared images showing nebulosity possibly tracing an outflow cavity, strongly suggest the existence of outflow activity. Provided that an envelope of at most ~0.1 Msunis the only mass accretion reservoir for L328-IRS, and the star formation efficiency is close to the canonical value ~30%, L328-IRS has not yet accreted more than 0.05 Msun. At the assumed distance of 200 pc, L328-IRS is destined to be a brown dwarf.

(Abridged) The c2d Spitzer Legacy project obtained images and photometry with both IRAC and MIPS instruments for five large, nearby molecular clouds. This paper combines information drawn from studies of individual clouds into a combined and updated statistical analysis of star formation rates and efficiencies, numbers and lifetimes for SED classes, and clustering properties. Current star formation efficiencies range from 3% to 6%. Taken together, the five clouds are producing about 260 solar masses of stars per Myr. The star formation surface density is more than an order of magnitude larger than would be predicted from the Kennicutt relation used in extragalactic studies. Measured against the dense gas probed by the maps of dust continuum emission, the efficiencies are much higher, and the current stock of dense cores would be exhausted in 1.8 Myr on average. The derived lifetime for the Class I phase is 0.44 to 0.54 Myr, considerably longer than some estimates. Similarly, the lifetime for the Class 0 SED class, 0.10 to 0.16 Myr, is longer than early estimates. The great majority (90%) of young stars lie within loose clusters with at least 35 members and a stellar density of 1 solar mass per cubic pc. Accretion at the sound speed from an isothermal sphere over the lifetime derived for the Class I phase could build a star of about 0.25 solar masses, given an efficiency of 0.3. Our data confirm and aggravate the "luminosity problem" for protostars. Our results strongly suggest that accretion is time variable, with prolonged periods of very low accretion. Based on a very simple model and this sample of sources, half the mass of a star would be accreted during only 7% of the Class I lifetime, as represented by the eight most luminous objects.

We present a search for all embedded protostars with internal luminosities < 1 solar luminosity in the sample of nearby, low-mass star-forming regions surveyed by the Spitzer Space Telescope c2d Legacy Project. The internal luminosity (Lint) of a source is the luminosity of the central source and excludes luminosity arising from external heating. On average, the c2d data are sensitive to embedded protostars with Lint > 4E-3 (d/140 pc)^2 solar luminosities, a factor of 25 better than the sensitivity of IRAS to such objects. We present selection criteria used to identify candidates from the Spitzer data and examine complementary data to decide whether each candidate is truly an embedded protostar. We find a tight correlation between the 70 micron flux and internal luminosity of a protostar, an empirical result based on observations and two-dimensional radiative transfer models of protostars. We identify 50 embedded protostars with Lint < 1 solar luminosities; 15 have Lint < 0.1 solar luminosities. The intrinsic distribution of source luminosities increases to lower luminosities. While we find sources down to the above sensitivity limit, indicating that the distribution may extend to luminosities lower than probed by these observations, we are able to rule out a continued rise in the distribution below 0.1 solar luminosities. Between 75-85% of cores classified as starless prior to being observed by Spitzer remain starless to our luminosity sensitivity; the remaining 15-25% harbor low-luminosity, embedded protostars. We compile complete Spectral Energy Distributions for all 50 objects and calculate standard evolutionary signatures, and argue that these objects are inconsistent with the simplest picture of star formation wherein mass accretes from the core onto the protostar at a constant rate.

We report the discovery of a nearby, embedded cluster of young stellar objects, associated filamentary infrared dark cloud, and 4.5 micron shock emission knots from outflows detected in Spitzer/IRAC mid-infrared imaging of the Serpens-Aquila Rift obtained as part of the Spitzer Gould Belt Legacy Survey. We also present radial velocity measurements of the region from molecular line observations obtained with the Submillimeter Array (SMA) that suggest the cluster is co-moving with the Serpens Main embedded cluster 3 degrees to the north. We therefore assign it the same distance, 260 pc. The core of the new cluster, which we call Serpens South, is composed of an unusually large fraction of protostars (77%) at high mean surface density (>430 pc^-2) and short median nearest neighbor spacing (3700 AU). We perform basic cluster structure characterization using nearest neighbor surface density mapping of the YSOs and compare our findings to other known clusters with equivalent analyses available in the literature.

We present the results of a submillimeter survey of 53 low-mass dense cores with the Submillimeter High Angular Resolution Camera II (SHARC-II). The survey is a follow-up project to the Spitzer Legacy Program ``From Molecular Cores to Planet-Forming Disks'', with the purpose being to create a complete data set of nearby low-mass dense cores from the infrared to the millimeter. We present maps of 52 cores at 350 microns and three cores at 450 microns, two of which were observed at both wavelengths. Of these 52 cores, 41 were detected by SHARC-II: 32 contained one submillimeter source while 9 contained multiple sources. For each submillimeter source detected, we report various source properties including source position, fluxes in various apertures, size, aspect ratio, and position angle. For the 12 cores that were not detected we present upper limits. The sources detected by SHARC-II have, on average, smaller sizes at the 2sigma contours than those derived from longer-wavelength bolometer observations. We conclude that this is not caused by a failure to integrate long enough to detect the full extent of the core; instead it arises primarily from the fact that the observations presented in this survey are insensitive to smoothly varying extended emission. We find that SHARC-II observations of low-mass cores are much better suited to distinguishing between starless and protostellar cores than observations at longer wavelengths. Very Low Luminosity Objects, a new class of objects being discovered by the Spitzer Space Telescope in cores previously classified as starless, look very similar at 350 microns to other cores with more luminous protostars.

We present IRAC (3.6, 4.5, 5.8, and 8.0 micron) observations of the Chamaeleon II molecular cloud. The observed area covers about 1 square degree defined by $A_V >2$. Analysis of the data in the 2005 c2d catalogs reveals a small number of sources (40) with properties similar to those of young stellaror substellar objects (YSOs). The surface density of these YSO candidates is low, and contamination by background galaxies appears to be substantial, especially for sources classified as Class I or flat SED. We discuss this problem in some detail and conclude that very few of the candidate YSOs in early evolutionary stages are actually in the Cha II cloud. Using a refined set of criteria, we define a smaller, but more reliable, set of 24 YSO candidates.

Infrared images of the dark cloud core B59 were obtained with the Spitzer Space Telescope as part of the "Cores to Disks" Legacy Science project. Photometry from 3.6-70 microns indicates at least 20 candidate low-mass young stars near the core, more than doubling the previously known population. Out of this group, 13 are located within about 0.1 pc in projection of the molecular gas peak, where a new embedded source is detected. Spectral energy distributions span the range from small excesses above photospheric levels to rising in the mid-infrared. One other embedded object, probably associated with the millimeter source B59-MMS1, with a bolometric luminosity L(bol) roughly 2 L(sun), has extended structure at 3.6 and 4.5 microns, possibly tracing the edges of an outflow cavity. The measured extinction through the central part of the core is A(V) greater than of order 45 mag. The B59 core is producing young stars with a high efficiency.

We present Spitzer Space Telescope observations of the "evolved starless core" L1521F which reveal the presence of a very low luminosity object (L < 0.07 Lsun). The object, L1521F-IRS, is directly detected at mid-infrared wavelengths (>5 micron) but only in scattered light at shorter infrared wavelengths, showing a bipolar nebula oriented east-west which is probably tracing an outflow cavity. The nebula strongly suggests that L1521F-IRS is embedded in the L1521F core. Thus L1521F-IRS is similar to the recently discovered L1014-IRS and the previously known IRAM 04191 in its substellar luminosity and dense core environment. However these objects differ significantly in their core density, core chemistry, and outflow properties, and some may be destined to be brown dwarfs rather than stars.

We report the first detections of the Class 0 protostellar source IRAM 04191+1522 at wavelengths shortward of 60 microns with the Spitzer Space Telescope. We see extended emission in the Spitzer images that suggests the presence of an outflow cavity in the circumstellar envelope. We combine the Spitzer observations with existing data to form a complete dataset ranging from 3.6 to 1300 microns and use these data to construct radiative transfer models of the source. We conclude that the internal luminosity of IRAM 04191+1522, defined to be the sum of the luminosity from the internal sources (a star and a disk), is L_int = 0.08 +/- 0.04 L_sun, placing it among the lowest luminosity protostars known. Though it was discovered before the launch of the Spitzer Space Telescope, IRAM 04191+1522 falls within a new class of Very Low Luminosity Objects being discovered by Spitzer. Unlike the two other well-studied objects in this class, which are associated either with weak, compact outflows or no outflows at all, IRAM 04191+1522 has a well-defined molecular outflow with properties consistent with those expected based on relations derived from higher luminosity (L_int > 1 L_sun) protostars. We discuss the difficulties in understanding IRAM 04191+1522 in the context of the standard model of star formation, and suggest a possible explanation for the very low luminosity of this source.