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

The sub-mm sky is a unique window for probing the architecture of the Universe and structures within it. From the discovery of dusty sub-mm galaxies, to the ringed nature of protostellar disks, our understanding of the formation, destruction, and evolution of objects in the Universe requires a comprehensive view of the sub-mm sky. The current generation single-dish sub-mm facilities have shown of the potential for discovery, while interferometers have presented a high resolution view into the finer details. However, our understanding of large-scale structure and our full use of these interferometers is now hampered by the limited sensitivity of our sub-mm view of the universe at larger scales. Thus, now is the time to start planning the next generation of sub-mm single dish facilities, to build on these revolutions in our understanding of the sub-mm sky. Here we present the case for the Atacama Large Aperture Submillimeter Telescope (AtLAST), a concept for a 50m class single dish telescope. We envision AtLAST as a facility operating as an international partnership with a suite of instruments to deliver the transformative science described in many Astro2020 science white papers. A 50m telescope with a high throughput and 1$^\circ$ FoV with a full complement of advanced instrumentation, including highly multiplexed high-resolution spectrometers, continuum cameras and Integral Field Units, AtLAST will have mapping speeds thousands of times greater than any current or planned facility. It will reach confusion limits below $L_*$ in the distant universe and resolve low-mass protostellar cores at the distance of the Galactic Center, providing synergies with upcoming facilities across the spectrum. Located on the Atacama plateau, to observe frequencies un-obtainable by other observatories, AtLAST will enable a fundamentally new understanding of the sub-mm universe at unprecedented depths.

We present a high-resolution ($\sim0.''12$, $\sim16$ au, mean sensitivity of $50~\mu$Jy~beam$^{-1}$ at 225 GHz) snapshot survey of 32 protoplanetary disks around young stars with spectral type earlier than M3 in the Taurus star-forming region using Atacama Large Millimeter Array (ALMA). This sample includes most mid-infrared excess members that were not previously imaged at high spatial resolution, excluding close binaries and highly extincted objects, thereby providing a more representative look at disk properties at 1--2 Myr. Our 1.3 mm continuum maps reveal 12 disks with prominent dust gaps and rings, 2 of which are around primary stars in wide binaries, and 20 disks with no resolved features at the observed resolution (hereafter smooth disks), 8 of which are around the primary star in wide binaries. The smooth disks were classified based on their lack of resolved substructures, but their most prominent property is that they are all compact with small effective emission radii ($R_{\rm eff,95\%} \lesssim 50$ au). In contrast, all disks with $R_{\rm eff,95\%}$ of at least 55 au in our sample show detectable substructures. Nevertheless, their inner emission cores (inside the resolved gaps) have similar peak brightness, power law profiles, and transition radii to the compact smooth disks, so the primary difference between these two categories is the lack of outer substructures in the latter. These compact disks may lose their outer disk through fast radial drift without dust trapping, or they might be born with small sizes. The compact dust disks, as well as the inner disk cores of extended ring disks, that look smooth at the current resolution will likely show small-scale or low-contrast substructures at higher resolution. The correlation between disk size and disk luminosity correlation demonstrates that some of the compact disks are optically thick at millimeter wavelengths.

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.

Young stars accrete mass from a circumstellar disk, but at the same time disk and star eject outflows and jets. These outflows have an onion-like structure where the innermost and fastest layers are surrounded by increasingly lower velocity components. The outer layers are probably photo-evaporative and magnetocentrifugally launched disk winds, but the nature of the inner winds is still uncertain. Since the fastest components carry only a small fraction of the mass, they are best observed at high-energies (X-ray and UV) as the slower, more massive components do not reach plasma temperatures sufficient for relevant X-ray or UV emission. Outflows are the most likely way in which a star or its disk can shed angular momentum and allow accretion to proceed; thus we cannot understand the accretion and the rotation rate of young stars if we cannot solve the origin of the inner jet components. Stellar jets share characteristics with their counterparts in more massive astrophysical objects, such as stellar mass black holes and AGN, with the added benefit that young stars are found at much closer distances and thus scales not accessible in other types of objects can be resolved. To understand the origin and impact of the inner jets, sub-arcsecond imaging and spectroscopy in the UV and X-rays is required, together with theory and modelling to interpret existing and future observations.

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

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

ALMA has observed a plethora of ring-like structures in planet forming discs at distances of 10-100 au from their host star. Although several mechanisms have been invoked to explain the origin of such rings, a common explanation is that they trace new-born planets. Under the planetary hypothesis, a natural question is how to reconcile the apparently high frequency of gap-carving planets at 10-100 au with the paucity of Jupiter mass planets observed around main sequence stars at those separations. Here, we provide an analysis of the new-born planet population emerging from observations of gaps in discs, under the assumption that the observed gaps are due to planets. We use a simple estimate of the planet mass based on the gap morphology, and apply it to a sample of gaps recently obtained by us in a survey of Taurus with ALMA. We also include additional data from recent published surveys, thus analysing the largest gap sample to date, for a total of 48 gaps. The properties of the purported planets occupy a distinctively different region of parameter space with respect to the known exo-planet population, currently not accessible through planet finding methods. Thus, no discrepancy in the mass and radius distribution of the two populations can be claimed at this stage. We show that the mass of the inferred planets conforms to the theoretically expected trend for the minimum planet mass needed to carve a dust gap. Finally, we estimate the separation and mass of the putative planets after accounting for migration and accretion, for a range of evolutionary times, finding a good match with the distribution of cold Jupiters.

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.

Earth's earliest sedimentary record contains evidence that surface temperatures were similar to, or perhaps even warmer than modern. In contrast, standard Solar models suggest the Sun was 25% less luminous at this ancient epoch, implying a cold, frozen planet-all else kept equal. This discrepancy, known as the Faint Young Sun Paradox, remains unresolved. Most proposed solutions invoke high concentrations of greenhouse gases in the early atmosphere to offset for the fainter Sun, though current geological constraints are insufficient to verify or falsify these scenarios. In this work, we examined several simple mechanisms that involve the role played by Earth's spin rate, which was significantly faster during Archean time. This faster spin rate enhances the equator-to-pole temperature gradient, facilitating a warm equator, while maintaining cold poles. Results show that such an enhanced meridional gradient augments the meridional gradient in carbonate deposition, which biases the surviving geological record away from the global mean, toward warmer waters. Moreover, using simple atmospheric models, we found that the faster-spinning Earth was less sensitive to ice-albedo feedbacks, facilitating larger meridional temperature gradients before succumbing to global glaciation. We show that within the faster-spinning regime, the greenhouse warming required to generate an ice-free Earth can differ from that required to generate an Earth with permanent ice caps by the equivalent of 1-2 orders of magnitude of pCO2. Accordingly, the resolution of the Faint Young Sun problem depends significantly on whether the early Earth was ever, or even at times, ice-free.

Outbursts due to dramatic increases in the mass accretion rate are the most extreme type of variability in young stellar objects. We searched for outbursts among 319 protostars in the Orion molecular clouds by comparing 3.6, 4.5, and 24 micron photometry from the Spitzer Space Telescope to 3.4, 4.6, and 22 micron photometry from the Wide-field Infrared Survey Explorer (WISE) obtained ~ 6.5 yr apart. Sources that brightened by more than two standard deviations above the mean variability at all three wavelengths were marked as burst candidates, and they were inspected visually to check for false positives due primarily to the reduced angular resolution of WISE compared to Spitzer. We recovered the known burst V2775 Ori (HOPS 223) as well as a previously unknown burst, HOPS 383, which we announced in an earlier paper. No other outbursts were found. With observations over 6.5 yr, we estimate an interval of about 1000 yr between bursts with a 90% confidence interval of 690 to 40,300 yr. The most likely burst interval is shorter than those found in studies of optically revealed young stellar objects, suggesting that outbursts are more frequent in protostars than in pre-main-sequence stars that lack substantial envelopes.

Far-infrared astronomy has advanced rapidly since its inception in the late 1950's, driven by a maturing technology base and an expanding community of researchers. This advancement has shown that observations at far-infrared wavelengths are important in nearly all areas of astrophysics, from the search for habitable planets and the origin of life, to the earliest stages of galaxy assembly in the first few hundred million years of cosmic history. The combination of a still developing portfolio of technologies, particularly in the field of detectors, and a widening ensemble of platforms within which these technologies can be deployed, means that far-infrared astronomy holds the potential for paradigm-shifting advances over the next decade. In this review, we examine current and future far-infrared observing platforms, including ground-based, sub-orbital, and space-based facilities, and discuss the technology development pathways that will enable and enhance these platforms to best address the challenges facing far-infrared astronomy in the 21st century.

We study a very young star-forming region in the outer Galaxy that is the most concentrated source of outflows in the Spitzer Space Telescope GLIMPSE360 survey. This region, dubbed CMa-l224, is located in the Canis Major OB1 association. CMa-l224 is relatively faint in the mid-infrared, but it shines brightly at the far-infrared wavelengths as revealed by the Herschel Space Observatory data from the Hi-GAL survey. Using the 3.6 and 4.5 $\mu$m data from the Spitzer/GLIMPSE360 survey, combined with the JHK$_s$ 2MASS and the 70-500 $\mu$m Herschel/Hi-GAL data, we develop a young stellar object (YSO) selection criteria based on color-color cuts and fitting of the YSO candidates' spectral energy distributions with YSO 2D radiative transfer models. We identify 293 YSO candidates and estimate physical parameters for 210 sources well-fit with YSO models. We select an additional 47 sources with GLIMPSE360-only photometry as `possible YSO candidates'. The vast majority of these sources are associated with high H$_2$ column density regions and are good targets for follow-up studies. The distribution of YSO candidates at different evolutionary stages with respect to Herschel filaments supports the idea that stars are formed in the filaments and become more dispersed with time. Both the supernova-induced and spontaneous star formation scenarios are plausible in the environmental context of CMa-l224. However, our results indicate that a spontaneous gravitational collapse of filaments is a more likely scenario. The methods developed for CMa-l224 can be used for larger regions in the Galactic plane where the same set of photometry is available.

Models of the Sun's long-term evolution suggest that its luminosity was substantially reduced 2-4 billion years ago, which is inconsistent with substantial evidence for warm and wet conditions in the geological records of both ancient Earth and Mars. Typical solutions to this so-called "faint young Sun paradox" consider changes in the atmospheric composition of Earth and Mars, and while attractive, geological verification of these ideas is generally lacking-particularly for Mars. One possible underexplored solution to the faint young Sun paradox is that the Sun has simply lost a few percent of its mass during its lifetime. If correct, this would slow, or potentially even offset the increase in luminosity expected from a constant-mass model. However, this hypothesis is challenging to test. Here, we propose a novel observational proxy of the Sun's ancient mass that may be readily measured from accumulation patterns in sedimentary rocks on Earth and Mars. We show that the orbital parameters of the Solar system planets undergo quasi-cyclic oscillations at a frequency, given by secular mode g_2-g_5, that scales approximately linearly with the Sun's mass. Thus by examining the cadence of sediment accumulation in ancient basins, it is possible distinguish between the cases of a constant mass Sun and a more massive ancient Sun to a precision of greater than about 1 per cent. This approach provides an avenue toward verification, or of falsification, of the massive early Sun hypothesis.

Gap-like structures in protoplanetary disks are likely related to planet formation processes. In this paper, we present and analyze high resolution (0.17*0.11 arcsec) 1.3 mm ALMA continuum observations of the protoplanetary disk around the Herbig Ae star MWC 480. Our observations for the first time show a gap centered at ~74au with a width of ~23au, surrounded by a bright ring centered at ~98au from the central star. Detailed radiative transfer modeling of both the ALMA image and the broadband spectral energy distribution is used to constrain the surface density profile and structural parameters of the disk. If the width of the gap corresponds to 4~8 times the Hill radius of a single forming planet, then the putative planet would have a mass of 0.4~3 M_Jup. We test this prediction by performing global three-dimensional smoothed particle hydrodynamic gas/dust simulations of disks hosting a migrating and accreting planet. We find that the dust emission across the disk is consistent with the presence of an embedded planet with a mass of ~2.3 M_Jup at an orbital radius of ~78au. Given the surface density of the best-fit radiative transfer model, the amount of depleted mass in the gap is higher than the mass of the putative planet, which satisfies the basic condition for the formation of such a planet.

Rings are the most frequently revealed substructure in ALMA dust observations of protoplanetary disks, but their origin is still hotly debated. In this paper, we identify dust substructures in 12 disks and measure their properties to investigate how they form. This subsample of disks is selected from a high-resolution ($\sim0.12''$) ALMA 1.33 mm survey of 32 disks in the Taurus star-forming region, which was designed to cover a wide range of sub-mm brightness and to be unbiased to previously known substructures. While axisymmetric rings and gaps are common within our sample, spiral patterns and high contrast azimuthal asymmetries are not detected. Fits of disk models to the visibilities lead to estimates of the location and shape of gaps and rings, the flux in each disk component, and the size of the disk. The dust substructures occur across a wide range of stellar mass and disk brightness. Disks with multiple rings tend to be more massive and more extended. The correlation between gap locations and widths, the intensity contrast between rings and gaps, and the separations of rings and gaps could all be explained if most gaps are opened by low-mass planets (super-Earths and Neptunes) in the condition of low disk turbulence ($\alpha=10^{-4}$). The gap locations are not well correlated with the expected locations of CO and N$_2$ ice lines, so condensation fronts are unlikely to be a universal mechanism to create gaps and rings, though they may play a role in some cases.

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

We analyze Herschel Space Observatory observations of 104 young stellar objects with protoplanetary disks in the ~1.5 Myr star-forming region Lynds 1641 (L1641) within the Orion A Molecular Cloud. We present spectral energy distributions from the optical to the far-infrared including new photometry from the Herschel Photodetector Array Camera and Spectrometer (PACS) at 70 microns. Our sample, taken as part of the Herschel Orion Protostar Survey, contains 24 transitional disks, eight of which we identify for the first time in this work. We analyze the full disks with irradiated accretion disk models to infer dust settling properties. Using forward modeling to reproduce the observed nKS-[70] index for the full disk sample, we find the observed disk indices are consistent with models that have depletion of dust in the upper layers of the disk relative to the mid plane, indicating significant dust settling. We perform the same analysis on full disks in Taurus with Herschel data and find that Taurus is slightly more evolved, although both samples show signs of dust settling. These results add to the growing literature that significant dust evolution can occur in disks by ~1.5 Myr.

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

We present new HST ACS observations and detailed models for a recently discovered edge-on protoplanetary disk around ESO H$\alpha$ 569 (a low-mass T Tauri star in the Cha I star forming region). Using radiative transfer models we probe the distribution of the grains and overall shape of the disk (inclination, scale height, dust mass, flaring exponent and surface/volume density exponent) by model fitting to multiwavelength (F606W and F814W) HST observations together with a literature compiled spectral energy distribution. A new tool set was developed for finding optimal fits of MCFOST radiative transfer models using the MCMC code emcee to efficiently explore the high dimensional parameter space. It is able to self-consistently and simultaneously fit a wide variety of observables in order to place constraints on the physical properties of a given disk, while also rigorously assessing the uncertainties in those derived properties. We confirm that ESO H$\alpha$ 569 is an optically thick nearly edge-on protoplanetary disk. The shape of the disk is well described by a flared disk model with an exponentially tapered outer edge, consistent with models previously advocated on theoretical grounds and supported by millimeter interferometry. The scattered light images and spectral energy distribution are best fit by an unusually high total disk mass (gas+dust assuming a ratio of 100:1) with a disk-to-star mass ratio of 0.16.

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

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

We present new SOFIA-FORCAST observations obtained in Feburary 2016 of the archetypal outbursting low mass young stellar object FU Orionis, and compare the continuum, solid state, and gas properties with mid-IR data obtained at the same wavelengths in 2004 with Spitzer-IRS. In this study, we conduct the first mid-IR spectroscopic comparison of an FUor over a long time period. Over a 12 year period, UBVR monitoring indicates that FU Orionis has continued its steady decrease in overall brightness by ~ 14%. We find that this decrease in luminosity occurs only at wavelengths < 20 microns. In particular, the continuum short ward of the silicate emission complex at 10 microns exhibits a ~ 12% (~ 3 sigma) drop in flux density, but no apparent change in slope; both the Spitzer and SOFIA spectra are consistent with a 7200 K blackbody. Additionally, the detection of water absorption is consistent with the Spitzer spectrum. The silicate emission feature at 10 microns continues to be consistent with unprocessed grains, unchanged over 12 years. We conclude that either the accretion rate in FU Orionis has decreased by ~ 12-14% over this time baseline, or that the inner disk has cooled, but the accretion disk remains in a superheated state outside of the innermost region.

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

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

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

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

With the Wide-Field Infrared Survey Explorer (WISE), we searched for young stellar objects (YSOs) in a 100 deg^2 region centered on the lightly studied Canis Major star forming region. Applying stringent magnitude cuts to exclude the majority of extragalactic contaminants, we find 144 Class I candidates and 335 Class II candidates. The sensitivity to Class II candidates is limited by their faintness at the distance to Canis Major (assumed as 1000 pc). More than half the candidates (53%) are found in 16 groups of more than four members, including four groups with more than 25 members each. The ratio of Class II to Class I objects, N_II/N_I, varies from 0.4 to 8.3 in just the largest four groups. We compare our results to those obtainable with combined Two Micron All Sky Survey (2MASS) and post-cryogenic Spitzer Space Telescope data; the latter approach recovers missing Class II sources. Via a comparison to protostars characterized with the Herschel Space Observatory, we propose new WISE color criteria for flat-spectrum and Class 0 protostars, finding 80 and seven of these, respectively. The distribution of YSOs in CMa OB1 is consistent with supernova-induced star formation, although the diverse N_II/N_I ratios are unexpected if this parameter traces age and the YSOs are due to the same supernova. Less massive clouds feature larger N_II/N_I ratios, suggesting that initial conditions play a role in determining this quantity.

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

We present a near-IR survey for the visual multiples in the Orion molecular clouds region at separations between 100 and 1000 AU. These data were acquired at 1.6~$\mu$m with the NICMOS and WFC3 cameras on the Hubble Space Telescope. Additional photometry was obtained for some of the sources at 2.05~$\mu$m with NICMOS and in the $L'$-band with NSFCAM2 on the IRTF. Towards 129 protostars and 197 pre-main sequence stars with disks observed with WFC3, we detect 21 and 28 candidate companions between the projected separations of 100---1000 AU, of which less than 5 and 8, respectively, are chance line of sight coincidences. The resulting companion fraction ($CF$) after the correction for the line of sight contamination is 14.4$^{+1.1}_{-1.3}$% for protostars and 12.5$^{+1.2}_{-0.8}$% for the pre-main sequence stars. These values are similar to those found for main sequence stars, suggesting that there is little variation in the $CF$ with evolution, although several observational biases may mask a decrease in the $CF$ from protostars to the main sequence stars. After segregating the sample into two populations based on the surrounding surface density of YSOs, we find that the $CF$ in the high stellar density regions ($\Sigma_{YSO} > 45$~pc$^{-2}$) is approximately 50% higher than that found in the low stellar density regions ($\Sigma_{YSO} < 45$~pc$^{-2}$). We interpret this as evidence for the elevated formation of companions at 100 to 1000 AU in the denser environments of Orion. We discuss possible reasons for this elevated formation.

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

While searches for young stellar objects (YSOs) with the Spitzer Space Telescope focused on known molecular clouds, photometry from the Wide-field Infrared Survey Explorer (WISE) can be used to extend the search to the entire sky. As a precursor to more expansive searches, we present results for a 100 square degree region centered on the Canis Major clouds.

During the "WISE at 5: Legacy and Prospects" conference in Pasadena, CA -- which ran from February 10 - 12, 2015 -- attendees were invited to engage in an interactive session exploring the future uses of the Wide-field Infrared Survey Explorer (WISE) data. The 65 participants -- many of whom are extensive users of the data -- brainstormed the top questions still to be answered by the mission, as well as the complementary current and future datasets and additional processing of WISE/NEOWISE data that would aid in addressing these most important scientific questions. The results were mainly bifurcated between topics related to extragalactic studies (e.g. AGN, QSOs) and substellar mass objects. In summary, participants found that complementing WISE/NEOWISE data with cross-correlated multiwavelength surveys (e.g. SDSS, Pan-STARRS, LSST, Gaia, Euclid, etc.) would be highly beneficial for all future mission goals. Moreover, developing or implementing machine-learning tools to comb through and understand cross-correlated data was often mentioned for future uses. Finally, attendees agreed that additional processing of the data such as co-adding WISE and NEOWISE and extracting a multi-epoch photometric database and parallax and proper motion catalog would greatly improve the scientific results of the most important projects identified. In that respect, a project such as MaxWISE which would execute the most important additional processing and extraction as well as make the data and catalogs easily accessible via a public portal was deemed extremely important.

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

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

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

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

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

In accreting young stars one of the prominent spectral features in the near infrared is the Paschen and Brackett series in emission. We examine hydrogen line ratios for 16 classical T Tauri stars from SpeX spectra and and assess the trends with veiling and accretion. The observed line ratios are compared to two theoretical models for line formation: (1) Baker and Menzel's (1938) Case B for radiative ionization and recombination and (2) a set of local line excitation calculations designed to replicate the conditions in T Tauri winds and magnetic accretion columns (Kwan & Fischer 2011). While the comparison between Case B and observed line ratios implies a wide range in electron density and temperature among the hydrogen line formation regions in T Tauri stars, the predictions of the local line excitation models give consistent results across multiple diagnostics. Under the assumptions of the local line excitation calculations, we find that n_H in the hydrogen line formation region is constrained to 2x10^10 -- 2x10^11 cm^-3, where stars with higher accretion rates have densities at the higher end of this range. Because of uncertainties in extinction, temperature is not well delineated but falls within the range expected for collisional excitation to produce the line photons. We introduce new diagnostics for assessing extinction based on near infrared hydrogen line ratios from the local line excitation calculations.

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

We report extensive new photometry and spectroscopy of the highly variable young stellar object PTF 10nvg including optical and near-infrared time series data as well as mid-infrared and millimeter data. Following the previously reported 2010 rise, during 2011 and 2012 the source underwent additional episodes of brightening and dimming events including prolonged faint states. The observed high-amplitude variations are largely consistent with extinction changes having a 220 day quasi-periodic signal. Spectral evolution includes not only changes in the spectral slope but correlated variation in the prominence of TiO/VO/CO bands and atomic line emission, as well as anticorrelated variation in forbidden line emission which, along with H_2, dominates optical and infrared spectra at faint epochs. Neutral and singly-ionized atomic species are likely formed in an accretion flow and/or impact while the origin of zero-velocity atomic LiI 6707 in emission is unknown. Forbidden lines, including several rare species, exhibit blueshifted emission profiles and likely arise from an outflow/jet. Several of these lines are also seen spatially offset from the continuum source position, presumably in a shocked region of an extended jet. CARMA maps resolve on larger scales a spatially extended outflow in mm-wavelength CO. We attribute the observed photometric and spectroscopic behavior in terms of occultation of the central star as well as the bright inner disk and the accretion/outflow zones that renders shocked gas in the inner part of the jet amenable to observation at the faint epochs. We discuss PTF 10nvg as a source exhibiting both accretion-driven (perhaps analogous to V1647 Ori) and extinction-driven (perhaps analogous to UX Ori or GM Cep) high-amplitude variability phenomena.

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

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

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

We examine eight young stellar objects in the OMC-2 star forming region based on observations from the SOFIA/FORCAST early science phase, the Spitzer Space Telescope, the Herschel Space Observatory, 2MASS, APEX, and other results in the literature. We show the spectral energy distributions of these objects from near-infrared to millimeter wavelengths, and compare the SEDs with those of sheet collapse models of protostars and circumstellar disks. Four of the objects can be modelled as protostars with infalling envelopes, two as young stars surrounded by disks, and the remaining two objects have double-peaked SEDs. We model the double-peaked sources as binaries containing a young star with a disk and a protostar. The six most luminous sources are found in a dense group within a 0.15 x 0.25 pc region; these sources have luminosities ranging from 300 L_sun to 20 L_sun. The most embedded source (OMC-2 FIR 4) can be fit by a class 0 protostar model having a luminosity of ~50 L_sun and mass infall rate of ~10^-4 solar masses per year.

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

We present the first characterization of the excess continuum emission of accreting T Tauri stars between optical and near-infrared wavelengths. With nearly simultaneous spectra from 0.48 to 2.4 microns acquired with HIRES and NIRSPEC on Keck and SpeX on the IRTF, we find significant excess continuum emission throughout this region, including the I, Y, and J bands, which are usually thought to diagnose primarily photospheric emission. The IYJ excess correlates with the excess in the V band, attributed to accretion shocks in the photosphere, and the excess in the K band, attributed to dust in the inner disk near the dust sublimation radius, but it is too large to be an extension of the excess from these sources. The spectrum of the excess emission is broad and featureless, suggestive of blackbody radiation with a temperature between 2200 and 5000 K. The luminosity of the IYJ excess is comparable to the accretion luminosity inferred from modeling the blue and ultraviolet excess emission and may require reassessment of disk accretion rates. The source of the IYJ excess is unclear. In stars of low accretion rate, the size of the emitting region is consistent with cooler material surrounding small hot accretion spots in the photosphere. However, for stars with high accretion rates, the projected area is comparable to or exceeds that of the stellar surface. We suggest that at least some of the IYJ excess emission arises in the dust-free gas inside the dust sublimation radius in the disk.

During a synoptic survey of the North American Nebula region, the Palomar Transient Factory (PTF) detected an optical outburst (dubbed PTF10nvg) associated with the previously unstudied flat or rising spectrum infrared source IRAS 20496+4354. The PTF R-band light curve reveals that PTF10nvg brightened by more than 5 mag during the current outburst, rising to a peak magnitude of R~13.5 in 2010 Sep. Follow-up observations indicate PTF10nvg has undergone a similar ~5 mag brightening in the K band, and possesses a rich emission-line spectrum, including numerous lines commonly assumed to trace mass accretion and outflows. Many of these lines are blueshifted by ~175 km/s from the North American Nebula's rest velocity, suggesting that PTF10nvg is driving an outflow. Optical spectra of PTF10nvg show several TiO/VO bandheads fully in emission, indicating the presence of an unusual amount of dense (> 10^10 cm^-3), warm (1500-4000 K) circumstellar material. Near-infrared spectra of PTF10nvg appear quite similar to a spectrum of McNeil's Nebula/V1647 Ori, a young star which has undergone several brightenings in recent decades, and 06297+1021W, a Class I protostar with a similarly rich near--infrared emission line spectrum. While further monitoring is required to fully understand this event, we conclude that the brightening of PTF10nvg is indicative of enhanced accretion and outflow in this Class-I-type protostellar object, similar to the behavior of V1647 Ori in 2004-2005.

We perform local excitation calculations to obtain line opacities and emissivity ratios and compare them with observed properties of H, He I, O I, Ca II, and Na I lines to determine the density, temperature, and photon ionization rate. We find that UV photoionization is the most probable excitation mechanism for generating the He I 10830 opacities that produce all the associated absorption features. We also calculate the specific line flux at an observed velocity of v_obs = +/- 150 km/s for both radial wind and infall models. All the model results, together with observed correlations between absorption and emission features and between narrow and broad emission components, are used to deduce the origins of the strong H, He I, and Ca II broad line emission. We conclude that the first two arise primarily in a radial outflow that is highly clumpy. The bulk of the wind volume is filled by gas at a density ~10^9 cm^-3 and optically thick to He I 10830 and H alpha, but optically thin to He I 5876, Pa gamma, and the Ca II infrared triplet. The optically thick He I 5876 emission occurs mostly in regions of density greater than or equal to 10^11 cm^-3 and temperature greater than or equal to 1.5x10^4 K, while the optically thick H alpha and Pa gamma emission occur mostly in regions of density around 10^11 cm^-3 and temperature between 8750 and 1.25x10^4 K. In producing the observed line fluxes at a given v_obs, the covering factor of these emission clumps is sufficiently small not to incur significant absorption of the stellar and veiling continua in either He I or H lines. The strong Ca II broad line emission likely arises in both the magnetospheric accretion flow and the disk boundary layer, where the gases dissipate part of their rotational energies before infalling along magnetic field lines. The needed density and temperature are ~10^12 cm^-3 and less than or equal to 7500 K.

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

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

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

We take advantage of good atmospheric transparency and the availability of high quality instrumentation in the 1 um near-infrared atmospheric window to present a grid of F, G, K, and M spectral standards observed at high spectral resolution (R ~ 25,000). In addition to a spectral atlas, we present a catalog of atomic line absorption features in the 0.95-1.11 um range. The catalog includes a wide range of line excitation potentials, from 0-13 eV, arising from neutral and singly ionized species, most frequently those of Fe I and Ti I at low excitation, Cr I, Fe I, and Si I at moderate excitation, and C I, S I, and Si I having relatively high excitation. The spectra also include several prominent molecular bands from CN and FeH. For the atomic species, we analyze trends in the excitation potential, line depth, and equivalent width across the grid of spectroscopic standards to identify temperature and surface gravity diagnostics near 1 um. We identify the line ratios that appear especially useful for spectral typing as those involving Ti I and C I or S I, which are temperature sensitive in opposite directions, and Sr II, which is gravity sensitive at all spectral types. Ascii versions of all spectra are available to download with the electronic version of the journal.

We probe the geometry of magnetospheric accretion in classical T Tauri stars by modeling red absorption at He I 10830 via scattering of the stellar and veiling continua. Under the assumptions that the accretion flow is an azimuthally symmetric dipole and helium is sufficiently optically thick that all incident 1-micron radiation is scattered, we illustrate the sensitivity of He I 10830 red absorption to both the size of the magnetosphere and the filling factor of the hot accretion shock. We compare model profiles to those observed in 21 CTTS with subcontinuum redshifted absorption at He I 10830 and find that about half of the stars have red absorptions and 1-micron veilings that are consistent with dipole flows of moderate width with accretion shock filling factors matching the size of the magnetospheric footpoints. However, the remaining 50% of the profiles, with a combination of broad, deep absorption and low 1-micron veiling, require very wide flows where magnetic footpoints are distributed over 10-20% of the stellar surface but accretion shock filling factors are < 1%. We model these profiles by invoking large magnetospheres dilutely filled with accreting gas, leaving the disk over a range of radii in many narrow "streamlets" that fill only a small fraction of the entire infall region. In some cases accreting streamlets need to originate in the disk between several stellar radii and at least the corotation radius. A few stars have such deep absorption at velocities greater than half the stellar escape velocity that flows near the star with less curvature than a dipole trajectory seem to be required.

The high opacity of He I 10830 makes it an exceptionally sensitive probe of the inner wind geometry of accreting T Tauri stars. In this line blueshifted absorption below the continuum results from simple scattering of stellar photons, a situation which is readily modeled without definite knowledge of the physical conditions and recourse to multi-level radiative transfer. We present theoretical line profiles for scattering in two possible wind geometries, a disk wind and a wind emerging radially from the star, and compare them to observed He I 10830 profiles from a survey of classical T Tauri stars. The comparison indicates that subcontinuum blueshifted absorption is characteristic of disk winds in ~30% of the stars and of stellar winds in ~40%. We further conclude that for many stars the emission profile of helium likely arises in stellar winds, increasing the fraction of accreting stars inferred to have accretion-powered stellar winds to ~60%. Stars with the highest disk accretion rates are more likely to have stellar wind than disk wind signatures and less likely to have redshifted absorption from magnetospheric funnel flows. This suggests the possibility that when accretion rates are high, disks can extend closer to the star, magnetospheric accretion zones can be reduced in size and conditions arise that favor radially outflowing stellar winds.

In a high dispersion 1 micron survey of 39 classical T Tauri stars veiling is detected in 80%, and He I 10830 and Pgamma line emission in 97% of the stars. On average, the 1 micron veiling exceeds the level expected from previously identified sources of excess emission, suggesting the presence of an additional contributor to accretion luminosity in the star-disk interface region. Strengths of both lines correlate with veiling, and at Pgamma there is a systematic progression in profile morphology with veiling. He I 10830 has an unprecedented sensitivity to inner winds, showing blueshifted absorption below the continuum in 71% of the CTTS compared to 0% at pgamma. This line is also sensitive to magnetospheric accretion flows, with redshifted absorption below the continuum found in 47% of the CTTS compared to 24% at Pgamma. The blueshifted absorption at 10830 shows considerable diversity in its breadth and penetration depth into the continuum, indicating that a range of inner wind conditions exist in accreting stars. We interpret the broadest and deepest blue absorptions as formed from scattering of the 1 micron continuum by outflowing gas whose full acceleration region envelopes the star, suggesting radial outflow from the star. In contrast, narrow blue absorption with a range of radial velocities more likely arises via scattering of the 1 micron continuum by a wind emerging from the inner disk. Both stellar and disk winds are accretion powered since neither is seen in non-accreting WTTS and among the CTTS helium strength correlates with veiling.

He I 10830 profiles acquired with Keck's NIRSPEC for 6 young low mass stars with high disk accretion rates (AS 353A, DG Tau, DL Tau, DR Tau, HL Tau and SVS 13) provide new insight into accretion-driven winds. In 4 stars the profiles have the signature of resonance scattering, and possess a deep and broad blueshifted absorption that penetrates more than 50% into the 1 micron continuum over a continuous range of velocities from near the stellar rest velocity to the terminal velocity of the wind, unlike inner wind signatures seen in other spectral features. This deep and broad absorption provides the first observational tracer of the acceleration region of the inner wind and suggests that this acceleration region is situated such that it occults a significant portion of the stellar disk. The remaining 2 stars also have blue absorption extending below the continuum although here the profiles are dominated by emission, requiring an additional source of helium excitation beyond resonant scattering. This is likely the same process that produces the emission profiles seen at He I 5876.

We present deep infrared ($2.2 \mu$m) imaging of the Galactic microquasars 1E1740-2942 and GRS 1758-258 using the Keck-I 10-meter telescope in June 1998. The observations were taken under excellent seeing conditions ($\sim 0.45 \arcsec$ full-width half-maximum), making them exceptionally deep for these crowded fields. We used the USNO-A2.0 catalog to astrometrically calibrate the infrared images (along with an optical CCD image in the case of GRS 1758-258), providing independent frame ties to the known radio positions of the objects. For 1E1740-2942, we confirm potential candidates for the microquasar previously identified by Marti et al., and show that none of the objects near the microquasar have varied significantly from 1998 to 1999. For GRS 1758-258, our astrometry indicates a position shifted from previous reports of candidates for the microquasar. We find no candidates inside our 90% confidence radius to a $2 \sigma$ limiting magnitude of $K_s = 20.3$ mag. We discuss the implications of these results for the nature of the microquasar binary systems.