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.

The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc- to core-scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flow along the magnetic field toward dense cores. When comparing the observed cores identified from the GBT Ammonia Survey (GAS) and Planck polarization-inferred magnetic field orientations, we find that the relative core-field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core-field orientation could be used to probe the relative significance of the magnetic field within the cloud.

We observed radio recombination lines (RRLs) toward the W51 molecular cloud complex, one of the most active star forming regions in our Galaxy. The UV radiation from young massive stars ionizes gas surrounding them to produce HII regions. Observations of the W51 IRS1 HII region were made with the Arecibo 305 m telescope. Of the full 1-10 GHz database, we have analyzed the observations between 4.5 and 5 GHz here. The steps involved in the analysis were: a) bandpass calibration using on-source/off-source observations; b) flux density calibration; c) removing spectral baselines due to errors in bandpass calibration and d) Gaussian fitting of the detected lines. We detected alpha, beta and gamma transitions of hydrogen and alpha transitions of helium. We used the observed line parameters to 1) measure the source velocity (56.6 $\pm$ 0.3 km s$^{-1}$) with respect to the Local Standard of Rest (LSR); 2) estimate the electron temperature (8500 $\pm$ 1800 K) of the HII region and 3) derive the emission measure (5.4 $\pm$ 2.7 $\times$ 10$^{6}$ pc cm$^{-6}$) of the ionized gas.

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 conducted an exploration of 12CO molecular outflows in the Orion A giant molecular cloud to investigate outflow feedback using 12CO (J = 1-0) and 13CO (J = 1-0) data obtained by the Nobeyama 45-m telescope. In the region excluding the center of OMC 1, we identified 44 12CO (including 17 newly detected) outflows based on the unbiased and systematic procedure of automatically determining the velocity range of the outflows and separating the cloud and outflow components. The optical depth of the 12CO emission in the detected outflows is estimated to be approximately 5. The total momentum and energy of the outflows, corrected for optical depth, are estimated to be 1.6 x 10 2 M km s-1 and 1.5 x 10 46 erg, respectively. The momentum and energy ejection rate of the outflows are estimated to be 36% and 235% of the momentum and energy dissipation rates of the cloud turbulence, respectively. Furthermore, the ejection rates of the outflows are comparable to those of the expanding molecular shells estimated by Feddersen et al. (2018, ApJ, 862, 121). Cloud turbulence cannot be sustained by the outflows and shells unless the energy conversion efficiency is as high as 20%.

We have investigated the formation and kinematics of sub-mm continuum cores in the Orion A molecular cloud. A comparison between sub-mm continuum and near infrared extinction shows a continuum core detection threshold of $A_V\sim$ 5-10 mag. The threshold is similar to the star formation extinction threshold of $A_V\sim$ 7 mag proposed by recent work, suggesting a universal star formation extinction threshold among clouds within 500 pc to the Sun. A comparison between the Orion A cloud and a massive infrared dark cloud G28.37+0.07 indicates that Orion A produces more dense gas within the extinction range 15 mag $\lesssim A_V \lesssim$ 60 mag. Using data from the CARMA-NRO Orion Survey, we find that dense cores in the integral-shaped filament (ISF) show sub-sonic core-to-envelope velocity dispersion that is significantly less than the local envelope line dispersion, similar to what has been found in nearby clouds. Dynamical analysis indicates that the cores are bound to the ISF. An oscillatory core-to-envelope motion is detected along the ISF. Its origin is to be further explored.

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.

Jets and outflows are ubiquitous phenomena in astrophysics, found in our Galaxy in diverse environments, from the formation of stars to late-type stellar objects. We present observations conducted with the 305m Arecibo Telescope of the pre-planetary nebula CRL 618 (Westbrook Nebula) - a well studied late-type star that has developed bipolar jets. The observations resulted in the first detection of 4765 MHz OH in a late-type stellar object. The line was narrow (FWHM ~ 0.6 km/s) and ~40 km/s blueshifted with respect to the systemic velocity, which suggests association with the expanding jets/bullets in CRL 618. We also report non-detection at Arecibo of any other OH transition between 1 and 9 GHz. The non-detections were obtained during the observations in 2008, when the 4765 MHz OH line was first discovered, and also in 2015 when the 4765 MHz OH line was not detected. Our data indicate that the 4765 MHz OH line was a variable maser. Modeling of the 4765 MHz OH detection and non-detection of the other transitions is consistent with the physical conditions expected in CRL 618. The 4765 MHz OH maser could originate from dissociation of H2O by shocks after sublimation of icy objects in this dying carbon-rich stellar system, although other alternatives such as OH in an oxygen-rich circumstellar region associated with a binary companion are also possible.

We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transition to coherence" in velocity dispersion around its periphery. The identification of these structures provides a chance to study the coherent structures in molecular clouds statistically. The identified coherent structures have a typical radius of 0.04 pc and a typical mass of 0.4 Msun, generally smaller than previously known coherent cores identified by Goodman et al. (1998), Caselli et al. (2002), and Pineda et al. (2010). We call these structures "droplets." We find that unlike previously known coherent cores, these structures are not virially bound by self-gravity and are instead predominantly confined by ambient pressure. The droplets have density profiles shallower than a critical Bonnor-Ebert sphere, and they have a velocity (VLSR) distribution consistent with the dense gas motions traced by NH3 emission. These results point to a potential formation mechanism through pressure compression and turbulent processes in the dense gas. We present a comparison with a magnetohydrodynamic simulation of a star-forming region, and we speculate on the relationship of droplets with larger, gravitationally bound coherent cores, as well as on the role that droplets and other coherent structures play in the star formation process.

We present ALMA CO(2-1) observations toward a massive infrared dark cloud G28.37+0.07. The ALMA data reveal numerous molecular (CO) outflows with a wide range of sizes throughout the cloud. Sixty-two 1.3 mm continuum cores were identified to be driving molecular outflows. We have determined the position angle in the plane-of-sky of 120 CO outflow lobes and studied their distribution. We find that the distribution of the plane-of-sky outflow position angles peaks at about 100 degree, corresponding to a concentration of outflows with an approximately east-west direction. For most outflows, we have been able to estimate the plane-of-sky angle between the outflow axis and the filament that harbors the protostar that powers the outflow. Statistical tests strongly indicate that the distribution of outflow-filament orientations is consistent with most outflow axes being mostly orthogonal to their parent filament in 3D. Such alignment may result from filament fragmentation or continuous mass transportation from filament to the embedded protostellar core. The latter is suggested by recent numerical studies with moderately strong magnetic fields.

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.

We investigate the relationship between turbulence and feedback in the Orion A molecular cloud using maps of $^{12}$CO(1-0), $^{13}$CO(1-0) and C$^{18}$O(1-0) from the CARMA-NRO Orion survey. We compare gas statistics with the impact of feedback in different parts of the cloud to test whether feedback changes the structure and kinematics of molecular gas. We use principal component analysis, the spectral correlation function, and the spatial power spectrum to characterize the cloud. We quantify the impact of feedback with momentum injection rates of protostellar outflows and wind-blown shells as well as the surface density of young stars. We find no correlation between shells or outflows and any of the gas statistics. However, we find a significant anti-correlation between young star surface density and the slope of the $^{12}$CO spectral correlation function, suggesting that feedback may influence this statistic. While calculating the principal components, we find peaks in the covariance matrix of our molecular line maps offset by 1-3 km s$^{-1}$ toward several regions of the cloud which may be produced by feedback. We compare these results to predictions from molecular cloud simulations.

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.

We present an initial overview of the filamentary structure in the Orion A molecular cloud utilizing a high angular and velocity resolution C$^{18}$O(1-0) emission map that was recently produced as part of the CARMA-NRO Orion Survey. The main goal of this study is to build a credible method to study varying widths of filaments which has previously been linked to star formation in molecular clouds. Due to the diverse star forming activities taking place throughout its $\sim$20 pc length, together with its proximity of 388 pc, the Orion A molecular cloud provides an excellent laboratory for such an experiment to be carried out with high resolution and high sensitivity. Using the widely-known structure identification algorithm, DisPerSE, on a 3-dimensional (PPV) C$^{18}$O cube, we identified 625 relatively short (the longest being 1.74 pc) filaments over the entire cloud. We study the distribution of filament widths using FilChaP, a python package that we have developed and made publicly available. We find that the filaments identified in a 2 square degree PPV cube do not overlap spatially, except for the complex OMC-4 region that shows distinct velocity components along the line of sight. The filament widths vary between 0.02 and 0.3 pc depending on the amount of substructure that a filament possesses. The more substructure a filament has, the larger is its width. We also find that despite this variation, the filament width shows no anticorrelation with the central column density which is in agreement with previous Herschel observations.

We present combined observations of the NH3 (J,K) = (1,1) and (2,2) inversion transitions towards OMC1 in Orion A obtained by the Karl G. Jansky Very Large Array (VLA) and the 100 m Robert C. Byrd Green Bank Telescope (GBT). With an angular resolution of 6" (0.01 pc), these observations reveal with unprecedented detail the complex filamentary structure extending north of the active Orion BN/KL region in a field covering 6' x 7'. We find a 0.012 pc wide filament within OMC1, with an aspect ratio of ~37:1, that was missed in previous studies. Its orientation is directly compared to the relative orientation of the magnetic field from the James Clerk Maxwell Telescope BISTRO survey in Orion A. We find a small deviation of ~11 deg between the mean orientation of the filament and the magnetic field, suggesting that they are almost parallel to one another. The filament's column density is estimated to be 2-3 orders of magnitude larger than the filaments studied with Herschel and is possibly self-gravitating given the low values of turbulence found. We further produce maps of the gas kinematics by forward modeling the hyperfine structure of the NH3 (J,K) = (1,1) and (2,2) lines. The resulting distribution of velocity dispersions peaks at ~0.5 km/s, close to the subsonic regime of the gas. This value is about 0.2 km/s smaller than previously measured in single-dish observations of the same region, suggesting that higher angular and spectral resolution observations will identify even lower velocity dispersions that might reach the subsonic turbulence regime in dense gas filaments.

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

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 .

We present the discovery of expanding spherical shells around low to intermediate-mass young stars in the Orion A giant molecular cloud using observations of $^{12}$CO (1-0) and $^{13}$CO (1-0) from the Nobeyama Radio Observatory 45-meter telescope. The shells have radii from 0.05 to 0.85 pc and expand outward at 0.8 to 5 km/s. The total energy in the expanding shells is comparable to protostellar outflows in the region. Together, shells and outflows inject enough energy and momentum to maintain the cloud turbulence. The mass-loss rates required to power the observed shells are two to three orders of magnitude higher than predicted for line-driven stellar winds from intermediate-mass stars. This discrepancy may be resolved by invoking accretion-driven wind variability. We describe in detail several shells in this paper and present the full sample in the online journal.

We present ALMA 850 $\mu$m continuum observations of the Orion Nebula Cluster that provide the highest angular resolution ($\sim 0\rlap{.}''1 \approx 40$ AU) and deepest sensitivity ($\sim 0.1$ mJy) of the region to date. We mosaicked a field containing $\sim 225$ optical or near-IR-identified young stars, $\sim 60$ of which are also optically-identified "proplyds". We detect continuum emission at 850 $\mu$m towards $\sim 80$% of the proplyd sample, and $\sim 50$% of the larger sample of previously-identified cluster members. Detected objects have fluxes of $\sim 0.5$-80 mJy. We remove sub-mm flux due to free-free emission in some objects, leaving a sample of sources detected in dust emission. Under standard assumptions of isothermal, optically thin disks, sub-mm fluxes correspond to dust masses of $\sim 0.5$ to 80 Earth masses. We measure the distribution of disk sizes, and find that disks in this region are particularly compact. Such compact disks are likely to be significantly optically thick. The distributions of sub-mm flux and inferred disk size indicate smaller, lower-flux disks than in lower-density star-forming regions of similar age. Measured disk flux is correlated weakly with stellar mass, contrary to studies in other star forming regions that found steeper correlations. We find a correlation between disk flux and distance from the massive star $\theta^1$ Ori C, suggesting that disk properties in this region are influenced strongly by the rich cluster environment.

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.

The nearby open cluster NGC 752 presents a rare opportunity to study stellar properties at ages >1 Gyr. However, constructing a membership catalog for it is challenging; most surveys have been limited to identifying its giants and dwarf members earlier than mid-K. We supplement past membership catalogs with candidates selected with updated photometric and proper-motion criteria, generating a list of 258 members, a >50% increase over previous catalogs. Using a Bayesian framework to fit MESA Isochrones & Stellar Tracks evolutionary models to literature photometry and the Tycho-Gaia Astrometric Solution data available for 59 cluster members, we infer the age of, and distance to, NGC 752: 1.34$\pm$0.06 Gyr and 438$_{-6}^{+8}$ pc. We also report the results of our optical monitoring of the cluster using the Palomar Transient Factory. We obtain rotation periods for 12 K and M cluster members, the first periods measured for such low-mass stars with a well-constrained age >1 Gyr. We compare these new periods to data from the younger clusters Praesepe and NGC 6811, and to a theoretical model for angular-momentum loss, to examine stellar spin down for low-mass stars over their first 1.3 Gyr. While on average NGC 752 stars are rotating more slowly than their younger counterparts, the difference is not significant. Finally, we use our spectroscopic observations to measure Halpha for cluster stars, finding that members earlier than $\approx$M2 are magnetically inactive, as expected at this age. Forthcoming Gaia data should solidify and extend the membership of NGC 752 to lower masses, thereby increasing its importance for studies of low-mass stars.

We present the first results from a new, high resolution, $^{12}$CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the Orion A cloud, hereafter referred to as the CARMA-NRO Orion Survey. CARMA observations have been combined with single-dish data from the Nobeyama 45m telescope to provide extended images at about 0.01 pc resolution, with a dynamic range of approximately 1200 in spatial scale. Here we describe the practical details of the data combination in uv space, including flux scale matching, the conversion of single dish data to visibilities, and joint deconvolution of single dish and interferometric data. A $\Delta$-variance analysis indicates that no artifacts are caused by combining data from the two instruments. Initial analysis of the data cubes, including moment maps, average spectra, channel maps, position-velocity diagrams, excitation temperature, column density, and line ratio maps provides evidence of complex and interesting structures such as filaments, bipolar outflows, shells, bubbles, and photo-eroded pillars. The implications for star formation processes are profound and follow-up scientific studies by the CARMA-NRO Orion team are now underway. We plan to make all the data products described here generally accessible; some are already available at https://dataverse.harvard.edu/dataverse/CARMA-NRO-Orion

Stars are born from dense cores in molecular clouds. Observationally, it is crucial to capture the formation of cores in order to understand the necessary conditions and rate of the star formation process. The {\it Atacama Large Mm/sub-mm Array} (ALMA) is extremely powerful for identifying dense gas structures, including cores, at mm wavelengths via their dust continuum emission. Here we use ALMA to carry out a survey of dense gas and cores in the central region of the massive ($\sim10^5\:M_\odot$) Infrared Dark Cloud (IRDC) G28.37+0.07. The observation consists of a mosaic of 86 pointings of the 12m-array and produces an unprecedented view of the densest structures of this IRDC. In this first paper about this data set, we focus on a comparison between the 1.3 mm continuum emission and a mid-infrared (MIR) extinction map of the IRDC. This allows estimation of the "dense gas" detection probability function (DPF), i.e., as a function of the local mass surface density, $\Sigma$, for various choices of thresholds of mm continuum emission to define "dense gas". We then estimate the dense gas mass fraction, $f_{\rm dg}$, in the central region of the IRDC and, via extrapolation with the DPF and the known $\Sigma$ probability distribution function, to the larger-scale surrounding regions, finding values of about 5\% to 15\% for the fiducial choice of threshold. We argue that this observed dense gas is a good tracer of the protostellar core population and, in this context, estimate a star formation efficiency per free-fall time in the central IRDC region of $\epsilon_{\rm ff}\sim$10\%, with approximately a factor of two systematic uncertainties.

We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc $\lesssim R_{eff} \lesssim$ 0.1 pc) and are spatially correlated with the highest H$_2$ column density portions of the cloud. We also compare the ammonia data to a catalog of dense cores identified by higher-resolution (18.2'' FWHM) Herschel Space Observatory observations of dust continuum emission from Cepheus-L1251. Maps of kinetic gas temperature, velocity dispersion, and NH$_3$ column density, derived from detailed modeling of the NH$_3$ data, are used to investigate the stability and chemistry of the ammonia-identified and Herschel-identified structures. We show that the dust and dense gas in the structures have similar temperatures, with median $T_{dust}$ and $T_K$ measurements of 11.7 $\pm$ 1.1 K and 10.3 $\pm$ 2.0 K, respectively. Based on a virial analysis, we find that the ammonia-identified structures are gravitationally dominated, yet may be in or near a state of virial equilibrium. Meanwhile, the majority of the Herschel-identified dense cores appear to be not bound by their own gravity and instead confined by external pressure. CCS $(2_0-1_0)$ and HC$_5$N $(9-8)$ emission from the region reveal broader line widths and centroid velocity offsets when compared to the NH$_3$ (1,1) emission in some cases, likely due to these carbon-based molecules tracing the turbulent outer layers of the dense cores.

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.

We use gas temperature and velocity dispersion data from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between self-gravity and the thermal and non-thermal motions present in the dense gas. Including the additional pressure binding imposed by the weight of the ambient molecular cloud material and additional smaller pressure terms, however, suggests that most of the dense cores are pressure confined.

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 an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with $A_V \gtrsim 7$ mag visible from the northern hemisphere in emission from NH$_3$ and other key molecular tracers. This first release includes the data for four regions in Gould Belt clouds: B18 in Taurus, NGC 1333 in Perseus, L1688 in Ophiuchus, and Orion A North in Orion. We compare the NH$_3$ emission to dust continuum emission from Herschel, and find that the two tracers correspond closely. NH$_3$ is present in over 60\% of lines-of-sight with $A_V \gtrsim 7$ mag in three of the four DR1 regions, in agreement with expectations from previous observations. The sole exception is B18, where NH$_3$ is detected toward ~ 40\% of lines-of-sight with $A_V \gtrsim 7$ mag. Moreover, we find that the NH$_3$ emission is generally extended beyond the typical 0.1 pc length scales of dense cores. We produce maps of the gas kinematics, temperature, and NH$_3$ column densities through forward modeling of the hyperfine structure of the NH$_3$ (1,1) and (2,2) lines. We show that the NH$_3$ velocity dispersion, ${\sigma}_v$, and gas kinetic temperature, $T_K$, vary systematically between the regions included in this release, with an increase in both the mean value and spread of ${\sigma}_v$ and $T_K$ with increasing star formation activity. The data presented in this paper are publicly available.

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.

Most massive stars form in dense clusters where gravitational interactions with other stars may be common. The two nearest forming massive stars, the BN object and Source I, located behind the Orion Nebula, were ejected with velocities of $\sim$29 and $\sim$13 km s$^{-1}$ about 500 years ago by such interactions. This event generated an explosion in the gas. New ALMA observations show in unprecedented detail, a roughly spherically symmetric distribution of over a hundred $^{12}$CO J=2$-$1 streamers with velocities extending from V$_{LSR}$ =$-$150 to +145 km s$^{-1}$. The streamer radial velocities increase (or decrease) linearly with projected distance from the explosion center, forming a `Hubble Flow' confined to within 50 arcseconds of the explosion center. They point toward the high proper-motion, shock-excited H$_2$ and [Fe ii ] `fingertips' and lower-velocity CO in the H$_2$ wakes comprising Orion's `fingers'. In some directions, the H$_2$ `fingers' extend more than a factor of two farther from the ejection center than the CO streamers. Such deviations from spherical symmetry may be caused by ejecta running into dense gas or the dynamics of the N-body interaction that ejected the stars and produced the explosion. This $\sim$10$^{48}$ erg event may have been powered by the release of gravitational potential energy associated with the formation of a compact binary or a protostellar merger. Orion may be the prototype for a new class of stellar explosion responsible for luminous infrared transients in nearby galaxies.

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.

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.

We present a 3 mm spectral line and continuum survey of L1451 in the Perseus Molecular Cloud. These observations are from the CARMA Large Area Star Formation Survey (CLASSy), which also imaged Barnard 1, NGC 1333, Serpens Main and Serpens South. L1451 is the survey region with the lowest level of star formation activity---it contains no confirmed protostars. HCO+, HCN, and N2H+ (J=1-0) are all detected throughout the region, with HCO+ the most spatially widespread, and molecular emission seen toward 90% of the area above N(H_2) column densities of 1.9x10^21 cm^-2. HCO+ has the broadest velocity dispersion, near 0.3 km/s on average, compared to ~0.15 km/s for the other molecules, thus representing a range from supersonic to subsonic gas motions. Our non-binary dendrogram analysis reveals that the dense gas traced by each molecule has similar hierarchical structure, and that gas surrounding the candidate first hydrostatic core (FHSC), L1451-mm, and other previously detected single-dish continuum clumps have similar hierarchical structure; this suggests that different sub-regions of L1451 are fragmenting on the pathway to forming young stars. We determined the three-dimensional morphology of the largest detectable dense gas structures to be relatively ellipsoidal compared to other CLASSy regions, which appeared more flattened at largest scales. A virial analysis shows the most centrally condensed dust structures are likely unstable against collapse. Additionally, we identify a new spherical, centrally condensed N2H+ feature that could be a new FHSC candidate. The overall results suggest L1451 is a young region starting to form its generation of stars within turbulent, hierarchical structures.

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.

FRIDA is a diffraction-limited imager and integral-field spectrometer that is being built for the adaptive-optics focus of the Gran Telescopio Canarias. In imaging mode FRIDA will provide scales of 0.010, 0.020 and 0.040 arcsec/pixel and in IFS mode spectral resolutions of 1500, 4000 and 30,000. FRIDA is starting systems integration and is scheduled to complete fully integrated system tests at the laboratory by the end of 2017 and to be delivered to GTC shortly thereafter. In this contribution we present a summary of its design, fabrication, current status and potential scientific applications.

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

We present SMA CO(2-1) observations toward the protostellar jet driven by SVS13A, a variable protostar in the NGC1333 star-forming region. The SMA CO(2-1) images show an extremely high-velocity jet composed of a series of molecular 'bullets'. Based on the SMA CO observations, we discover clear and large systematic velocity gradients, perpendicular to the jet axis, in the blueshifted and redshifted bullets. After discussing several alternative interpretations, such as twin-jets, jet precession, warped disk, and internal helical shock, we suggest that the systematic velocity gradients observed in the bullets result from the rotation of the SVS13A jet. From the SMA CO images, the measured rotation velocities are 11.7-13.7 km/s for the blueshifted bullet and 4.7+/-0.5 km/s for the redshifted bullet. The estimated specific angular momenta of the two bullets are comparable to those of dense cores, about 10 times larger than those of protostellar envelopes, and about 20 times larger than those of circumstellar disks. If the velocity gradients are due to the rotation of the SVS13A jet, the significant amount of specific angular momenta of the bullets indicates that the rotation of jets/outflows is a key mechanism to resolve the so-called 'angular momentum problem' in the field of star formation. The kinematics of the bullets suggests that the jet launching footprint on the disk has a radius of about 7.2-7.7 au, which appears to support the extended disk-wind model. We note that further observations are needed to comprehensively understand the kinematics of the SVS13A jet, in order to confirm the rotation nature of the bullets.

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.

Observations of nearby molecular clouds detect "shells", which are likely caused by winds from young main sequence stars. However, the progenitors of these observed features are not well characterized and the mass-loss rates inferred from the gas kinematics are several orders of magnitude greater than those predicted by atomic line-driven stellar wind models. We use magnetohydrodynamic simulations to model winds launching within turbulent molecular clouds and explore the impact of wind properties on cloud morphology and turbulence. We find that winds do not produce clear features in turbulent statistics such as the Fourier spectra of density and momentum but do impact the Fourier velocity spectrum. The density and velocity distribution functions, especially as probed by CO spectral lines, strongly indicate the presence and influence of winds. We show that stellar mass-loss rates for individual stars must be $\dot m_w \gtrsim 10^{-7}$ Msun yr$^{-1}$, similar to those estimated from observations, to reproduce shell properties. Consequently, we conclude that B and A-type main sequence stars have mass-loss rates several orders of magnitude larger that those predicted by models or that young stars are more variable than expected due to magnetic activity or accretion.

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.

Protostellar outflows provide a means to probe the accretion process of forming stars and their ability to inject energy into their surroundings. However, conclusions based on outflow observations depend upon the degree of accuracy with which their properties can be estimated. We examine the quality of Atacama Large Millimeter/submillimeter Array (ALMA) observations of protostellar outflows by producing synthetic $^{12}$CO(1-0) and $^{13}$CO(1-0) observations of numerical simulations. We use various ALMA configurations, observational parameters, and outflow inclinations to assess how accurately different assumptions and setups can recover underlying properties. We find that more compact arrays and longer observing times can improve the mass and momentum recovery by a factor of two. During the first $\sim$0.3 Myr of evolution, $^{12}$CO(1-0) is optically thick, even for velocities $|v|\ge 1$ km s$^{-1}$, and outflow mass is severely underestimated without an optical depth correction. Likewise, $^{13}$CO(1-0) is optically thick during the first $\simeq 0.1$ Myr. However, underestimation due to shorter observing time, missing flux, and optical depth are partially offset by the assumption of local thermodynamic equilibrium and higher excitation temperatures. Overall, we expect that full ALMA $^{13}$CO(1-0) observations of protostellar sources within 500 pc with observing times $\gtrsim 1$ hrs and assumed excitation temperatures of $T<20$K will reliably measure mass and line-of-sight momentum to within 20\%.

We present new 2 cm and 6 cm maps of H2CO, radio recombination lines, and the radio continuum in the W51 star forming complex acquired with Arecibo and the Green Bank Telescope at ~50" resolution. We use H2CO absorption to determine the relative line-of-sight positions of molecular and ionized gas. We measure gas densities using the H2CO densitometer, including continuous measurements of the dense gas mass fraction (DGMF) over the range $10^4$ cm$^{-3}$ < n(H$_2$) < $10^6$ cm$^{-3}$- this is the first time a dense gas mass fraction has been measured over a range of densities with a single data set. The DGMF in W51A is high,f >~70% above $n>10^4$ cm$^{-3}$, while it is low, f<20%, in W51 B. We did not detect any H2CO emission throughout the W51 GMC; all gas dense enough to emit under normal conditions is in front of bright continuum sources and therefore is seen in absorption instead. The data set has been made public at http://dx.doi.org/10.7910/DVN/26818. Conclusions. (1) The dense gas fraction in the W51 A and B clouds shows that W51 A will continue to form stars vigorously, while star formation has mostly ended in W51 B. The lack of dense, star-forming gas around W51 C indicates that collect-and-collapse is not acting or is inefficient in W51. (2) Ongoing high-mass star formation is correlated with n ~ 1x10$^5$ cm$^{-3}$ gas. Gas with n > 10$^4$ cm$^{-3}$ is weakly correlated with low and moderate mass star formation, but does not strongly correlate with high-mass star formation. (3) The nondetection of H$_2$CO emission implies that the emission detected in other galaxies, e.g. Arp 220, comes from high-density gas that is not directly affiliated with already-formed massive stars. Either the non-star-forming ISM of these galaxies is very dense, implying the star formation density threshold is higher, or H II regions have their emission suppressed.

The initial velocity dispersion of newborn stars is a major unconstrained aspect of star formation theory. Using near-infrared spectra obtained with the APOGEE spectrograph, we show that the velocity dispersion of young (1-2 Myr) stars in NGC 1333 is 0.92+/-0.12 km/s after correcting for measurement uncertainties and the effect of binaries. This velocity dispersion is consistent with the virial velocity of the region and the diffuse gas velocity dispersion, but significantly larger than the velocity dispersion of the dense, star-forming cores, which have a sub-virial velocity dispersion of 0.5 km/s. Since the NGC 1333 cluster is dynamically young and deeply embedded, this measurement provides a strong constraint on the initial velocity dispersion of newly-formed stars. We propose that the difference in velocity dispersion between stars and dense cores may be due to the influence of a 70 micro-Gauss magnetic field acting on the dense cores, or be the signature of a cluster with initial sub-structure undergoing global collapse.

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

Open clusters are collections of stars with a single, well-determined age, and can be used to investigate the connections between angular-momentum evolution and magnetic activity over a star's lifetime. We present the results of a comparative study of the relationship between stellar rotation and activity in two benchmark open clusters: Praesepe and the Hyades. As they have the same age and roughly solar metallicity, these clusters serve as an ideal laboratory for testing the agreement between theoretical and empirical rotation-activity relations at $\approx$600 Myr. We have compiled a sample of 720 spectra --- more than half of which are new observations --- for 516 high-confidence members of Praesepe; we have also obtained 139 new spectra for 130 high-confidence Hyads. We have collected rotation periods ($P_{rot}$) for 135 Praesepe members and 87 Hyads. To compare $H\alpha$ emission, an indicator of chromospheric activity, as a function of color, mass, and Rossby number $R_o$, we first calculate an expanded set of $\chi$ values, with which we can obtain the $H\alpha$ to bolometric luminosity ratio, $L_{H\alpha}/L_{bol}$, even when spectra are not flux-calibrated and/or stars lack reliable distances. Our $\chi$ values cover a broader range of stellar masses and colors (roughly equivalent to spectral types from K0 to M9), and exhibit better agreement between independent calculations, than existing values. We find no difference between the two clusters in their $H\alpha$ equivalent width or $L_{H\alpha}/L_{bol}$ distributions, and therefore take the merged $H\alpha$ and $P_{rot}$ data to be representative of 600-Myr-old stars. Our analysis shows that $H\alpha$ activity in these stars is saturated for $R_o\leq0.11^{+0.02}_{-0.03}$. Above that value activity declines as a power-law with slope $\beta=-0.73^{+0.16}_{-0.12}$, before dropping off rapidly at $R_o\approx0.4$...

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

We have used deep near-infrared observations with adaptive optics to discover a distributed population of low-mass protostars within the filamentary Infrared Dark Cloud G34.43+00.24. We use maps of dust emission at multiple wavelengths to determine the column density structure of the cloud. In combination with an empirically-verified model of the magnitude distribution of background stars, this column density map allows us to reliably determine overdensities of red sources that are due to embedded protostars in the cloud. We also identify protostars through their extended emission in K-band which comes from excited H2 in protostellar outflows or reflection nebulosity. We find a population of distributed low-mass protostars, suggesting that low-mass protostars may form earlier than, or contemporaneously with, high-mass protostars in such a filament. The low-mass protostellar population may also produce the narrow linewidth SiO emission observed in some clouds without high-mass protostars. Finally, we use a molecular line map of the cloud to determine the virial parameter per unit length along the filament and find that the highest mass protostars form in the most bound portion of the filament, as suggested by theoretical models.

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

We present line and continuum observations made with the Submillimeter Array (SMA) of the young stellar object Barnard 5 - IRS1 located in the Perseus molecular cloud. Our $^{12}$CO(2-1) line observations resolve the high-velocity bipolar northeast-southwest outflow associated with this source. We find that the outflowing gas shows different structures at three different velocity regimes, in both lobes, resulting in a spider-like morphology. In addition to the low-velocity, cone-like (wide-angle) lobes that have previously been observed, we report the presence of intermediate-velocity parabolic shells emerging very close to the Class I protostar, as well as high velocity molecular bullets that appear to be associated to the optical/IR jet emanating from this source. These compact high-velocity features reach radial velocities of about 50 km s$^{-1}$ away from the cloud velocity. We interpret the peculiar spider-like morphology is a result of the molecular material being entrained by a wind with both a collimated jet-like component and a wide-angle component. We suggest the outflow is in a transitional evolutionary phase between a mostly jet-driven flow and an outflow in which the entrainment is dominated by the wide-angle wind component. We also detect 1300 $\mu$m continuum emission at the position of the protostar, which likely arises from the dusty envelope and disk surrounding the protostar. Finally, we report the detection of $^{13}$CO(2-1) and SO(6$_5$-5$_4$) emission arising from the outflow and the location of the young stellar object.

In this review we focus on the role jets and outflows play in the star and planet formation process. Our essential question can be posed as follows: are jets/outflows merely an epiphenomenon associated with star formation or do they play an important role in mediating the physics of assembling stars both individually and globally? We address this question by reviewing the current state of observations and their key points of contact with theory. Our review of jet/outflow phenomena is organized into three length-scale domains: Source and Disk Scales ($0.1-10^2$ au) where the connection with protostellar and disk evolution theories is paramount; Envelope Scales ($10^2-10^5$ au) where the chemistry and propagation shed further light on the jet launching process, its variability and its impact on the infalling envelope; Parent Cloud Scales ($10^5-10^6$ au) where global momentum injection into cluster/cloud environments become relevant. Issues of feedback are of particular importance on the smallest scales where planet formation regions in a disk may be impacted by the presence of disk winds, irradiation by jet shocks or shielding by the winds. Feedback on envelope scales may determine the final stellar mass (core-to-star efficiency) and envelope dissipation. Feedback also plays an important role on the larger scales with outflows contributing to turbulent support within clusters including alteration of cluster star formation efficiencies (feedback on larger scales currently appears unlikely). A particularly novel dimension of our review is that we consider results on jet dynamics from the emerging field of High Energy Density Laboratory Astrophysics (HEDLA). HEDLA is now providing direct insights into the 3-D dynamics of fully magnetized, hypersonic, radiative outflows.

We investigate protostellar outflow evolution, gas entrainment, and star formation efficiency using radiation-hydrodynamic simulations of isolated, turbulent low-mass cores. We adopt an X-wind launching model, in which the outflow rate is coupled to the instantaneous protostellar accretion rate and evolution. We vary the outflow collimation angle from $\theta$=0.01-0.1 and find that even well collimated outflows effectively sweep up and entrain significant core mass. The Stage 0 lifetime ranges from 0.14-0.19 Myr, which is similar to the observed Class 0 lifetime. The star formation efficiency of the cores spans 0.41-0.51. In all cases, the outflows drive strong turbulence in the surrounding material. Although the initial core turbulence is purely solenoidal by construction, the simulations converge to approximate equipartition between solenoidal and compressive motions due to a combination of outflow driving and collapse. When compared to a simulation of a cluster of protostars, which is not gravitationally centrally condensed, we find that the outflows drive motions that are mainly solenoidal. The final turbulent velocity dispersion is about twice the initial value of the cores, indicating that an individual outflow is easily able to replenish turbulent motions on sub-parsec scales. We post-process the simulations to produce synthetic molecular line emission maps of $^{12}$CO, $^{13}$CO, and C$^{18}$O and evaluate how well these tracers reproduce the underlying mass and velocity structure.

Stars do not generally form in isolation. Instead, they form in clusters, and in these clustered environments newborn stars can have profound effects on one another and on their parent gas clouds. Feedback from clustered stars is almost certainly responsible for a number of otherwise puzzling facts about star formation: that it is an inefficient process that proceeds slowly when averaged over galactic scales; that most stars disperse from their birth sites and dissolve into the galactic field over timescales $\ll 1$ Gyr; and that newborn stars follow an initial mass function (IMF) with a distinct peak in the range $0.1 - 1$ $M_\odot$, rather than an IMF dominated by brown dwarfs. In this review we summarize current observational constraints and theoretical models for the complex interplay between clustered star formation and feedback.

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 morphology, kinematics and entrainment mechanism of the HH 46/47 molecular outflow were studied using new ALMA Cycle 0 observations. Results show that the blue and red lobes are strikingly different. We argue that these differences are partly due to contrasting ambient densities that result in different wind components having a distinct effect on the entrained gas in each lobe. A 29-point mosaic, covering the two lobes at an angular resolution of about 3", detected outflow emission at much higher velocities than previous observations, resulting in significantly higher estimates of the outflow momentum and kinetic energy than previous studies of this source, using the CO(1-0) line. The morphology and the kinematics of the gas in the blue lobe are consistent with models of outflow entrainment by a wide-angle wind, and a simple model describes the observed structures in the position-velocity diagram and the velocity-integrated intensity maps. The red lobe exhibits a more complex structure, and there is evidence that this lobe is entrained by a wide-angle wind and a collimated episodic wind. Three major clumps along the outflow axis show velocity distribution consistent with prompt entrainment by different bow shocks formed by periodic mass ejection episodes which take place every few hundred years. Position-velocity cuts perpendicular to the outflow cavity show gradients where the velocity increases towards the outflow axis, inconsistent with outflow rotation. Additionally, we find evidence for the existence of a small outflow driven by a binary companion.

We present observations of outflows in the star-forming region NGC 1333 using the Combined Array for Research in Millimeter-Wave Astronomy (CARMA). We combined the 12CO and 13CO (1-0) CARMA mosaics with data from the 14-m Five College Radio Astronomy Observatory (FCRAO) to probe the central, most dense and active region of this protostellar cluster at scales from 5'' to 7' (or 1000 AU to 0.5 pc at a distance of 235 pc). We map and identify 12CO outflows, and along with 13CO data we estimate their mass, momentum and energy. Within the 7'x7' map, the 5'' resolution allows for a detailed study of morphology and kinematics of outflows and outflow candidates, some of which were previously confused with other outflow emission in the region. In total, we identify 22 outflow lobes, as well as 9 dense circumstellar envelopes marked by continuum emission, of which 6 drive outflows. We calculate a total outflow mass, momentum and energy within the mapped region of 6 Msun, 19 Msun km/s, and 7x10^44 erg, respectively. Within this same region, we compare outflow kinematics with turbulence and gravitational energy, and we suggest that outflows are likely important agents for the maintenance of turbulence in this region. In the earliest stages of star formation, outflows do not yet contribute enough energy to totally disrupt the clustered region where most star formation is happening, but have the potential to do so as the protostellar sources evolve. Our results can be used to constrain outflow properties, such as outflow strength, in numerical simulations of outflow-driven turbulence in clusters.

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 use two 4.5micron Spitzer (IRAC) maps of the NGC 1333 region taken over approx. 7 yr interval to determine proper motions of its associated outflows. This is a first, successful attempt at obtaining proper motions of stellars outflow from Spitzer observations. For the outflow formed by the Herbig-Haro objects HH7, 8 and 10, we find proper motions of approx. 9-13 km/s, which are consistent with previously determined optical proper motions of these objects. We determine proper motions for a total of 8 outflows, ranging from approx. 10 to 100 km/s. The derived proper motions show that out of these 8 outflows, 3 have tangential velocities less or equal to 20 km/s. This result shows that a large fraction of the observed outflows have low intrinsic velocities, and that the low proper motions are not merely a projection effect.

We investigate the shape of the extinction law in two 1-degree square fields of the Perseus Molecular Cloud complex. We combine deep red-optical (r, i, and z-band) observations obtained using Megacam on the MMT with UKIDSS near-infrared (J, H, and K-band) data to measure the colours of background stars. We develop a new hierarchical Bayesian statistical model, including measurement error, intrinsic colour variation, spectral type, and dust reddening, to simultaneously infer parameters for individual stars and characteristics of the population. We implement an efficient MCMC algorithm utilising generalised Gibbs sampling to compute coherent probabilistic inferences. We find a strong correlation between the extinction (Av) and the slope of the extinction law (parameterized by Rv). Because the majority of the extinction toward our stars comes from the Perseus molecular cloud, we interpret this correlation as evidence of grain growth at moderate optical depths. The extinction law changes from the diffuse value of Rv = 3 to the dense cloud value of Rv = 5 as the column density rises from Av = 2 mags to Av = 10 mags. This relationship is similar for the two regions in our study, despite their different physical conditions, suggesting that dust grain growth is a fairly universal process.

A survey of the Milky Way disk and the Magellanic System at the wavelengths of the 21-cm atomic hydrogen (HI) line and three 18-cm lines of the OH molecule will be carried out with the Australian Square Kilometre Array Pathfinder telescope. The survey will study the distribution of HI emission and absorption with unprecedented angular and velocity resolution, as well as molecular line thermal emission, absorption, and maser lines. The area to be covered includes the Galactic plane (|b|< 10deg) at all declinations south of delta = +40deg, spanning longitudes 167deg through 360deg to 79deg at b=0deg, plus the entire area of the Magellanic Stream and Clouds, a total of 13,020 square degrees. The brightness temperature sensitivity will be very good, typically sigma_T ~ 1 K at resolution 30arcsec and 1 km/s. The survey has a wide spectrum of scientific goals, from studies of galaxy evolution to star formation, with particular contributions to understanding stellar wind kinematics, the thermal phases of the interstellar medium, the interaction between gas in the disk and halo, and the dynamical and thermal states of gas at various positions along the Magellanic Stream.

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 present results from three-dimensional, self-gravitating, radiation-hydrodynamic simulations of low-mass protostellar outflows. We construct synthetic observations in 12CO in order to compare with observed outflows and evaluate the effects of beam resolution and outflow orientation on inferred outflow properties. To facilitate the comparison, we develop a quantitative prescription for measuring outflow opening angles. Using this prescription, we demonstrate that, in both simulations and synthetic observations, outflow opening angles broaden with time similarly to observed outflows. However, the interaction between the outflowing gas and the turbulent core envelope produces significant asymmetry between the red and blue shifted outflow lobes. We find that applying a velocity cutoff may result in outflow masses that are underestimated by a factor 5 or more, and masses derived from optically thick CO emission further underpredict the mass of the high-velocity gas by a factor of 5-10. Derived excitation temperatures indicate that outflowing gas is hotter than the ambient gas with temperature rising over time, which is in agreement with the simulation gas temperatures. However, excitation temperatures are otherwise not well correlated with the actual gas temperature.

We present a study on the shells (and bubbles) in the Perseus molecular cloud using the COMPLETE survey large-scale 12CO(1-0) and 13CO(1-0) maps. The twelve shells reported here are spread throughout most of the Perseus cloud and have circular or arc-like morphologies with a range in radius of about 0.1 to 3 pc. Most of them have not been detected before most likely as maps of the region lacked the coverage and resolution needed to distinguish them. The majority of the shells are coincident with infrared nebulosity of similar shape and have a candidate powering source near the center. We suggest they are formed by the interaction of spherical or very wide-angle winds powered by young stars inside or near the Perseus molecular cloud -a cloud that is commonly considered to be mostly forming low-mass stars. Two of the twelve shells are powered by high-mass stars close to the cloud, while the others appear to be powered by low or intermediate mass stars in the cloud. We argue that winds with a mass loss rate of about 10^-8 to 10^-6 M_sun/yr are required to produce the observed shells. Our estimates indicate that the energy input rate from these stellar winds is similar to the turbulence dissipation rate. We conclude that in Perseus the total energy input from both collimated protostellar outflows and powerful spherical winds from young stars is sufficient to maintain the turbulence in the molecular cloud. Large scale molecular line and IR continuum maps of a sample of clouds will help determine the frequency of this phenomenon in other star forming regions.

We present the detection of a dust continuum source at 3-mm (CARMA) and 1.3-mm (SMA), and 12CO(2-1) emission (SMA) towards the L1451-mm dense core. These detections suggest a compact object and an outflow where no point source at mid-infrared wavelengths is detected using Spitzer. An upper limit for the dense core bolometric luminosity of 0.05 Lsun is obtained. By modeling the broadband SED and the continuum interferometric visibilities simultaneously, we confirm that a central source of heating is needed to explain the observations. This modeling also shows that the data can be well fitted by a dense core with a YSO and disk, or by a dense core with a central First Hydrostatic Core (FHSC). Unfortunately, we are not able to decide between these two models, which produce similar fits. We also detect 12CO(2-1) emission with red- and blue-shifted emission suggesting the presence of a slow and poorly collimated outflow, in opposition to what is usually found towards young stellar objects but in agreement with prediction from simulations of a FHSC. This presents the best candidate, so far, for a FHSC, an object that has been identified in simulations of collapsing dense cores. Whatever the true nature of the central object in L1451-mm, this core presents an excellent laboratory to study the earliest phases of low-mass star formation.

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.

We present a ~6.5'x8' Expanded Very Large Array (EVLA) mosaic observations of the NH3 (1,1) emission in the Barnard 5 region in Perseus, with an angular resolution of 6". This map covers the coherent region, where the dense gas presents subsonic non-thermal motions (as seen from single dish observations with the Green Bank Telescope, GBT). The combined EVLA and GBT observations reveal, for the first time, a striking filamentary structure (20" wide or 5,000 AU at the distance of Perseus) in this low-mass star forming region. The integrated intensity profile of this structure is consistent with models of an isothermal filament in hydrostatic equilibrium. The observed separation between the B5-IRS1 young stellar object (YSO), in the central region of the core, and the northern starless condensation matches the Jeans length of the dense gas. This suggests that the dense gas in the coherent region is fragmenting. The region observed displays a narrow velocity dispersion, where most of the gas shows evidence for subsonic turbulence, and where little spatial variations are present. It is only close to the YSO where an increase in the velocity dispersion is found, but still displaying subsonic non-thermal motions

We present a study of the structure and dynamics of the dense gas surrounding the HH 211-mm source, using VLA observations of the ammonia (1,1) and (2,2) inversion transitions. We find the envelope around this Class 0 source has an elongated geometry, extending about 104 AU in the direction perpendicular to the well-known HH 211 jet, and exhibits a velocity distribution consistent with rotation along the major axis. Our VLA observations indicate that the envelope is mostly in virial equilibrium. However, comparing our data with results from previous studies it appears that the gas within approximately 0.005 pc of the central protostar is undergoing dynamical collapse. The size of this collapsing radius may constrain the amount of mass that can eventually infall into the forming star. We also find that the envelope is mostly internally heated, most probably by radiation from the central protostar. In addition, we detect evidence of outflow-envelope interaction in the ammonia data. These include a velocity gradient in the dense gas along the outflow axis and significant line broadening that is spatially correlated with the jet and it might be the result of outflow-induced turbulence in the envelope.

We report the discovery of multiple condensations in the prestellar core candidate SMM1A in the R~CrA cloud, which may represent the earliest phase of core fragmentation observed thus far. The separation between the condensations is between 1000 and 2100 AU, and their masses range from about 0.1 to 0.2 M_sun. We find that the three condensations have extremely low bolometric luminosities (< 0.1 L_sun) and temperatures (< 20 K), indicating that these are young sources that have yet to form protostars. We suggest that these sources were formed through the fragmentation of an elongated prestellar core. Our results, in concert with other observed protostellar binary systems with separations in the scale of 1000 AU, support the scenario that prompt fragmentation in the isothermal collapse phase is an efficient mechanism for wide binary star formation, while the fragmentation in the subsequent adiabatic phase may be an additional mechanism for close (< 100 AU) binary star formation.

We present two Very Large Array observations of the pre-main-sequence star PV Cephei, taken with a separation of 10.5 years. These data show that the proper motions of this star are $\mu_\alpha \cos \delta = +10.9 \pm 3.0$ mas yr$^{-1}$; $\mu_\delta = +0.2 \pm 1.8$ mas yr$^{-1}$, very similar to those -previously known- of HD 200775, the B2Ve star that dominates the illumination of the nearby reflection nebula NGC 7023. This result suggests that PV Cephei is not a rapidly moving run-away star as suggested by previous studies. The large velocity of PV Cephei had been inferred from the systematic eastward displacement of the bisectors of successive pairs of Herbig Haro knots along its flow. These systematic shifts might instead result from an intrinsic dissymmetry in the ejection mechanisms, or from an asymmetric distribution of the circumstellar material.

We present a study on the impact of molecular outflows in the Perseus molecular cloud complex using the COMPLETE survey large-scale 12CO(1-0) and 13CO(1-0) maps. We used three-dimensional isosurface models generated in RA-DEC-Velocity space to visualize the maps. This rendering of the molecular line data allowed for a rapid and efficient way to search for molecular outflows over a large (~ 16 sq. deg.) area. Our outflow-searching technique detected previously known molecular outflows as well as new candidate outflows. Most of these new outflow-related high-velocity features lie in regions that have been poorly studied before. These new outflow candidates more than double the amount of outflow mass, momentum, and kinetic energy in the Perseus cloud complex. Our results indicate that outflows have significant impact on the environment immediately surrounding localized regions of active star formation, but lack the energy needed to feed the observed turbulence in the entire Perseus complex. This implies that other energy sources, in addition to protostellar outflows, are responsible for turbulence on a global cloud scale in Perseus. We studied the impact of outflows in six regions with active star formation within Perseus of sizes in the range of 1 to 4 pc. We find that outflows have enough power to maintain the turbulence in these regions and enough momentum to disperse and unbind some mass from them. We found no correlation between outflow strength and star formation efficiency for the six different regions we studied, contrary to results of recent numerical simulations. The low fraction of gas that potentially could be ejected due to outflows suggests that additional mechanisms other than cloud dispersal by outflows are needed to explain low star formation efficiencies in clusters.

Intermediate between the prestellar and Class 0 protostellar phases, the first core is a quasi-equilibrium hydrostatic object with a short lifetime and an extremely low luminosity. Recent MHD simulations suggest that the first core can even drive a molecular outflow before the formation of the second core (i.e., protostar). Using the Submillimeter Array and the Spitzer Space Telescope, we present high angular resolution observations towards the embedded dense core IRS2E in L1448. We find that source L1448 IRS2E is not visible in the sensitive Spitzer infrared images (at wavelengths from 3.6 to 70 um), and has weak (sub-)millimeter dust continuum emission. Consequently, this source has an extremely low bolometric luminosity (< 0.1 L_sun). Infrared and (sub-)millimeter observations clearly show an outflow emanating from this source; L1448 IRS2E represents thus far the lowest luminosity source known to be driving a molecular outflow. Comparisons with prestellar cores and Class 0 protostars suggest that L1448 IRS2E is more evolved than prestellar cores but less evolved than Class 0 protostars, i.e., at a stage intermediate between prestellar cores and Class 0 protostars. All these results are consistent with the theoretical predictions of the radiative/magneto hydrodynamical simulations, making L1448 IRS2E the most promising candidate of the first hydrostatic core revealed so far.

We present NH3 observations of the B5 region in Perseus obtained with the Green Bank Telescope (GBT). The map covers a region large enough (~11'x14') that it contains the entire dense core observed in previous dust continuum surveys. The dense gas traced by NH3(1,1) covers a much larger area than the dust continuum features found in bolometer observations. The velocity dispersion in the central region of the core is small, presenting subsonic non-thermal motions which are independent of scale. However, it is thanks to the coverage and high sensitivity of the observations that we present the detection, for the first time, of the transition between the coherent core and the dense but more turbulent gas surrounding it. This transition is sharp, increasing the velocity dispersion by a factor of 2 in less than 0.04 pc (the 31" beam size at the distance of Perseus, ~250 pc). The change in velocity dispersion at the transition is ~3 km/s/pc. The existence of the transition provides a natural definition of dense core: the region with nearly-constant subsonic non-thermal velocity dispersion. From the analysis presented here we can not confirm nor rule out a corresponding sharp density transition.

We report the detection of complex organic molecules in the young protostellar outflow L1157. We identify lines from HCOOCH3, CH3CN, HCOOH and C2H5OH at the position of the B1 shock in the blueshifted lobe, making it the first time that complex species have been detected towards a molecular outflow powered by a young low-mass protostar. The time scales associated with the warm outflow gas (< 2,000 yr) are too short for the complex molecules to have formed in the gas phase after the shock-induced sputtering of the grain mantles. It is more likely that the complex species formed in the surface of grains and were then ejected from the grain mantles by the shock. The formation of complex molecules in the grains of low-mass star forming regions must be relatively efficient, and our results show the importance of considering the impact of outflows when studying complex molecules around protostars. The relative abundance with respect to methanol of most of the detected complex molecules is similar to that of hot cores and molecular clouds in the galactic center region, which suggests that the mantle composition of the dust in the L1157 dark cloud is similar to dust in those regions.

We have mapped the starless core TMC-1C in a variety of molecular lines with the IRAM 30m telescope. High density tracers show clear signs of self-absorption and sub-sonic infall asymmetries are present in N2H+ (1-0) and DCO+ (2-1) lines. The inward velocity profile in N2H+ (1-0) is extended over a region of about 7,000 AU in radius around the dust continuum peak, which is the most extended ``infalling'' region observed in a starless core with this tracer. The kinetic temperature (~12 K) measured from C17O and C18O suggests that their emission comes from a shell outside the colder interior traced by the mm continuum dust. The C18O (2-1) excitation temperature drops from 12 K to ~10 K away from the center. This is consistent with a volume density drop of the gas traced by the C18O lines, from ~4x10^4 cm^-3 towards the dust peak to ~6x10^3 cm^-3 at a projected distance from the dust peak of 80" (or 11,000 AU). The column density implied by the gas and dust show similar N2H+ and CO depletion factors (f_D < 6). This can be explained with a simple scenario in which: (i) the TMC-1C core is embedded in a relatively dense environment (H2 ~10^4 cm^-3), where CO is mostly in the gas phase and the N2H+ abundance had time to reach equilibrium values; (ii) the surrounding material (rich in CO and N2H+) is accreting onto the dense core nucleus; (iii) TMC-1C is older than 3x10^5 yr, to account for the observed abundance of N2H+ across the core (~10^-10 w.r.t. H2); and (iv) the core nucleus is either much younger (~10^4 yr) or ``undepleted'' material from the surrounding envelope has fallen towards it in the past 10,000 yr.

We report on initial results of our Spitzer Cycle 2 program to observe the young massive star R Mon and its associated HH 39 Herbig-Haro object in the mid-infrared. Our program used all instruments on-board Spitzer to obtain deep images with IRAC of the HH 39 complex and of R Mon and its surroundings, a deep image of HH 39 at 24 and 70 $\mu$m with MIPS, and mid-infrared spectra with the SH, LH, and LL modules of IRS. The aim of this program is to study the physical links in a young massive star between accretion disk, outflows and jets, and sh ocks in the associated HH object. Our preliminary analysis reveals that several knots of HH 39 are clearly detected in most IRAC bands. In IRAC4 (8 $\mu$m), diffuse emission, probably from PAHs, appears as foreground emission covering the HH 39 emission. The HH 39 knots are detected at 24 microns, despite the fact that dust continuum emission covers the knots and shows the same structure as observed with IRAC4. The IRS spectra of HH 39 show weak evidence of [Ne II] 12.8 $\mu$m and 0--0 S(1) H$_2$ 17.0 $\mu$m lines. A more detailed analysis is, however, required due to the faintness of the Herbig-Haro knots. Finally, we obtained the SH and MIPS SED spectra of R Mon. A PAH emission feature at 11.3 $\mu$m is detected on top of the strong continuum; although no strong emission or absorption lines are observed, we will seek to detect faint lines. The combined IRAC, IRS, and MIPS data of the R Mon/HH 39 system will help us to understand circumstellar disk processing, and the connection between jets, outflows, and HH objects.

We present multi-line and continuum observations of the circumstellar environment within 10^4 AU of a sample of protostars to investigate how the effects of outflows on their immediate environment changes over time. 12CO(1-0) emission probes the high-velocity molecular outflows near the protostars and demonstrate that the outflow opening angle widens as the nascent star evolves. Maps of the 13CO(1-0) and HCO+(1-0) outflow emission show that protostellar winds erode the circumstellar envelope through the entrainment of the outer envelope gas. The spatial and velocity distribution of the dense circumstellar envelope, as well as its mass, is traced by the C18O(1-0) emission and also displays evolutionary changes. We show that outflows are largely responsible for these changes, and propose an empirical model for the evolution of outflow-envelope interactions. In addition, some of the outflows in our sample appear to affect the chemical composition of the surrounding environment, enhancing the HCO+ abundance. Overall, our results confirm that outflows play a major role in the star formation process through their strong physical and chemical impacts on the environments of the young protostars.

We review the known properties of molecular outflows from low- and high-mass young stars. General trends among outflows are identified, and the most recent studies on the morphology, kinematics, energetics, and evolution of molecular outflows are discussed, focusing on results from high-resolution millimeter observations. We review the existing four broad classes of outflow models and compare numerical simulations with the observational data. A single class of models cannot explain the range of morphological and kinematic properties that are observed, and we propose a possible solution. The impact of outflows on their cloud is examined, and we review how outflows can disrupt their surrounding environment, through the clearing of gas and the injection of momentum and energy onto the gas at distances from their powering sources from about 0.01 to a few pc. We also discuss the effects of shock-induced chemical processes on the ambient medium, and how these processes may act as a chemical clock to date outflows. Lastly, future outflow research with existing and planned millimeter and submillimeter instruments is presented.

We present an overview of data available for the Ophiuchus and Perseus molecular clouds from ``Phase I'' of the COMPLETE Survey of Star-Forming Regions. This survey provides a range of data complementary to the Spitzer Legacy Program ``From Molecular Cores to Planet Forming Disks.'' Phase I includes: Extinction maps derived from 2MASS near-infrared data using the NICER algorithm; extinction and temperature maps derived from IRAS 60 and 100um emission; HI maps of atomic gas; 12CO and 13CO maps of molecular gas; and submillimetre continuum images of emission from dust in dense cores. Not unexpectedly, the morphology of the regions appears quite different depending on the column-density tracer which is used, with IRAS tracing mainly warmer dust and CO being biased by chemical, excitation and optical depth effects. Histograms of column-density distribution are presented, showing that extinction as derived from 2MASS/NICER gives the closest match to a log-normal distribution as is predicted by numerical simulations. All the data presented in this paper, and links to more detailed publications on their implications are publically available at the COMPLETE website.

Three independent lines of evidence imply that the young star PV Ceph is moving at roughly 20 km/s through the interstellar medium. The first, and strongest, suggestion of motion comes from the geometry of the HH knots in the "giant" Herbig-Haro (HH) flow associated with PV Ceph. Bisectors of lines drawn between pairs of knots at nearly equal distances from PV Ceph imply an E-W motion of the source, and a plasmon model fit to the knot positions gives a good fit of 22 km/s motion for the star. The second bit of damning evidence comes from a redshifted "trail" of molecular gas, pointing in the same E-W direction implied by the HH knot geometry. The third exhibit we offer in accusing PV Ceph of speeding involves the tilt apparent in the high-velocity molecular jet now emanating from the star. This tilt is best explained if the true, current, jet direction is N-S, as it is in HST images, and the star is moving at roughly 20 km/s. Tracing the motion of PV Ceph backward in time, to the nearest cluster from which it might have been ejected, we find that it is very likely to have been thrown out of the massive star-forming cluster NGC 7023 (more than 10 pc away). We propose that PV Ceph was ejected, at a speed large enough to escape NGC 7023, at least 100,000 years ago, but that it did not enter the molecular cloud in which it now finds itself until more like 35,000 years ago. Our calculations show that the currently-observable molecular outflow associated with PV Ceph is about 10,000 years old, so that the flow has had plenty of time to form while in its current molecular cloud. But, the question of what PV Ceph was doing, and what gas/disk it took along with it in the time it was traveling through the low-density region between NGC 7023 and its current home is an open question.

We present a new set of high-resolution molecular line maps of the gas immediately surrounding various Herbig-Haro (HH) knots of the giant HH flow HH 315, from the young star PV Cephei. The observations, aimed at studying the entrainment mechanism of the 2.6 pc-long HH 315 flow, include IRAM 30m maps of the 12CO(2-1), 12CO(1-0), and 13CO(1-0) lines, with beam sizes of 11'', 21'', and 22'', respectively. We compare the morphology and the kinematics of the outflow gas, as well as the temperature and momentum distribution of the molecular outflow with those predicted by different entrainment models. With our detailed study we are able to conclude that jet bow shock entrainment by an episodic stellar wind, with a time-varying axis, produces most of the high-velocity molecular outflow observed far from the source. In addition, near PV Cephei we find evidence for a poorly collimated, wide-angle, molecular outflow and a collimated wiggling jet-like molecular outflow. We propose that the poorly collimated component is entrained by a wide-angle wind, and the collimated component is entrained by a variable jet with internal working surfaces. If this picture is true, then a stellar wind model which allows for the coexistence of a wide-angle component and a collimated (jet-like) stellar wind component is needed to explain the observed properties of the PV Ceph outflow. The wiggling axis of the redshifted molecular outflow lobe indicates that the outflow ejection axis is changing over time. We find that the time-scale of the axis variation shown by the molecular outflow lobe is about a factor of 10 less than that shown by the large-scale optical HH knots.

We present a set of detailed molecular line maps of the region associated with the giant Herbig-Haro flow HH 315, from the young star PV Cephei, aimed at studying the outflow-cloud interaction. Our study clearly shows that the HH 315 flow is effecting the kinematics of its surrounding medium, and has been able to redistribute considerable amounts of the surrounding medium-density gas in its star-forming core as well at parsec-scale distances from the source. The giant molecular outflow HH 315 is a highly asymmetric bipolar flow with a projected linear extent of about 2 pc. Our results indicate that the two outflow lobes are each interacting with the ambient medium in different ways. The southern (redshifted) lobe interacts with a dense ambient medium, very close to the young stellar outflow source, and its kinetic energy is comparable to both the turbulent and gravitational binding energy of its host cloud. In addition, we find evidence that the southern lobe is responsible for the creation of a cavity in the 13CO emission. In contrast, the northern (mainly blueshifted) outflow lobe extends farther from PV Ceph and interacts with ambient gas much less dense than the southern lobe. There is very little 13CO emission north of the outflow source, and the only prominent 13CO emission is a shell-like structure coincident with the outer edge of the northern lobe, about 1.2 pc northwest of PV Ceph. It appears that the northern lobe of the HH 315 outflow has been able to ``push'' aside a substantial fraction of the gas in the area, piling it in a dense shell-like structure at its edges. In addition, we find that the northern outflow lobe is responsible for a velocity gradient in the ambient gas.

While observational evidence for the episodic nature of young stellar outflows continues to mount, existing numerical and theoretical models of molecular outflows assume they are formed by the interaction of a non-episodic wind from a young stellar object with an ambient cloud. In this Letter we estimate and discuss the effects of episodicity on the mass-velocity and position-velocity relations observed in molecular outflows. We explain how many recent observational results disagree with the predictions of non-episodic outflow models, and we offer simple explanations for the discrepancies. In particular, we discuss how an episodic stellar outflow can steepen the power-law slope of the mass-velocity relation in a molecular outflow. And, we illustrate how an episodic outflow can produce multiple "Hubble-wedges'' in the position-velocity distribution of a molecular outflow. With a little more information than we have now, it may be possible to use the "fossil record" embedded in a molecular outflow's mass-velocity and position-velocity relations to reconstruct the history of a young stellar object's mass ejection episodes.

We present the first set of detailed molecular line maps of the region associated with the giant Herbig-Haro flow HH 300, from the young star IRAS 04239+2436. Our results indicate that the red lobe of the HH 300 flow is depositing a fair amount of momentum and kinetic energy over a notable volume of its host dark cloud. This makes HH 300 a key player in the evolution and fate of its parent cloud. The redshifted molecular outflow lobe of HH 300 is 1.1 pc long and 0.3 pc wide, and has a very clumpy structure. The density, velocity, and momentum distributions in the outflow indicate that the observed clumps arise from the prompt entrainment of ambient gas. Bow shock-shaped structures are associated with the clumps, and we hypothesize that the shocks are produced by different mass ejection episodes. Lines drawn from IRAS 04239+2436 to each of these clumps have different orientations on the plane of the sky, and we conclude that HH 300 is a precessing and episodic outflow. The observations include a map of the red lobe in the 12CO(2-1) line, with a beam size of 27'', and more extended maps of the outflow region in the 12CO(1-0) and 13CO(1-0) lines, with 45'' and 47'' beam sizes, respectively. We obtain a steep power-law mass spectrum for HH 300, which we believe is best explained by the evolution of the outflow mass kinematics. In addition, our 13CO observations show that the HH 300 flow has been able to redistribute (in space and velocity) considerable amounts of its surrounding medium-density gas.

We present a study of the detailed distribution of extinction in a region of the Taurus dark cloud complex. Our study uses new BVR images of the region, spectral classification data for 95 stars, and IRAS Sky Survey Atlas (ISSA) 60 and 100 micron images. We study the extinction of the region in four different ways, and we present the first inter-comparison of all these methods, which are: 1) using the color excess of background stars for which spectral types are known; 2) using the ISSA 60 and 100 micron images; 3) using star counts; and 4) using an optical (V and R) version of the average color excess method used by Lada et al. (1994). We find that all four methods give generally similar results, with important exceptions. To study the structure in the dust distribution, we compare the ISSA extinction and the extinction measured for individual stars. From the comparison, we conclude that in the relatively low extinction regions studied, with 0.9 < A_V < 3.0 mag (away from filamentary dark clouds and IRAS cores), there are no fluctuations in the dust column density greater than 45% (at the 99.7% confidence level), on scales smaller than 0.2 pc. We also report the discovery of a previously unknown stellar cluster behind the Taurus dark cloud near R.A 4h19m00s, Dec. 27:30:00 (B1950)

We test the recently published all-sky reddening map of Schlegel, Finkbeiner & Davis (1998 [SFD]) using the extinction study of a region in the Taurus dark cloud complex by Arce & Goodman (1999 [AG]). In their study, AG use four different techniques to measure the amount and structure of the extinction toward Taurus, and all four techniques agree very well. Thus we believe that the AG results are a truthful representation of the extinction in the region and can be used to test the reliability of the SFD reddening map. The results of our test show that the SFD all-sky reddening map, which is based on data from COBE/DIRBE and IRAS/ISSA, overestimates the reddening by a factor of 1.3 to 1.5 in regions of smooth extinction with A_V > 0.5 mag. In some regions of steep extinction gradients the SFD map underestimates the reddening value, probably due to its low spatial resolution. We expect that the astronomical community will be using the SFD reddening map extensively. We offer this Letter as a cautionary note about using the SFD map in regions of high extinction (A_V > 0.5 mag), as it might not be giving accurate reddening values there.

We present a study of the polarizing power of the dust in cold dense regions (dark clouds) compared to that of dust in the general interstellar medium (ISM). Our study uses new polarimetric, optical, and spectral classification data for 36 stars to carefully study the relation between polarization percentage (p) and extinction (A_V) in the Taurus dark cloud complex. We find two trends in our p-A_V study: (1) stars background to the warm ISM show an increase in p with A_V; and (2) the percentage of polarization of stars background to cold dark clouds does not increase with extinction. We detect a break in the p-A_V relation at an extinction 1.3 +/- 0.2 mag, which we expect corresponds to a set of conditions where the polarizing power of the dust associated with the Taurus dark clouds drops precipitously. This breakpoint places important restrictions on the use of polarimetry in studying interstellar magnetic fields.

We report the serendipitous discovery of a large HI cloud with an associated HI mass of $6(\pm1.5)\times 10^9 h^{-2}$ M$_\odot$ and a heliocentric velocity 8800 \kms, located near the periphery of the cluster of galaxies Abell 2634. Its velocity field appears to be very quiescent, as no gradients in the peak velocity are seen over its extent of 143$ h^{-1}$ by 103$ h^{-1}$ kpc. The distribution of gas is poorly resolved spatially, and it is thus difficult at this time to ascertain the nature of the cloud. At least two relatively small, actively star--forming galaxies appear to be embedded in the HI gas, which may (a) be an extended gaseous envelope surrounding one or both galaxies, (b) have been spread over a large region by a severe episode of tidal disruption or (c) have been affected by the ram pressure resulting from its motion through the intracluster gas of A2634.