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

This paper presents models to estimate the structure of density and magnetic field strength in spheroidal condensations, from maps of their column density and their polarization of magnetically aligned dust grains. The density model is obtained by fitting a column density map with an embedded p = 2 Plummer spheroid of any aspect ratio and inclination. The magnetic properties are based on the density model, on the Davis-Chandrasekhar-Fermi (DCF) model of Alfv\'enic fluctuations, and on the Spheroid Flux Freezing (SFF) model of mass and flux conservation in Paper I. The field strength model has the resolution of the column density map, which is finer than the resolution of the DCF estimate of field strength. The models are applied to ALMA observations of the envelope of the protostar BHR71 IRS1. Column density fits give the density model, from (2.0 +- 0.4) x 10^5 cm^-3 to (7 +- 1) x 10^7 cm^-3 . The density model predicts the field directions map, which fits the polarization map best within 1100 au, with standard deviation of angle differences 17{\deg}. In this region the DCF mean field strength is 0.7 +- 0.2 mG and the envelope mass is supercritical, with ratio of mass to magnetic critical mass 1.5 +- 0.4. The SFF field strength profile scales with the DCF field strength, from 60 x 10{\mu}G to 4+-1 mG. The spatial resolution of the SFF field strength estimate is finer than the DCF resolution by a factor ~7, and the peak SFF field strength exceeds the DCF field strength by a factor ~5.

We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particularly efficient. Methanol and sulphur monoxide are most abundant in a narrow zone that follows the eastern side of the core. This side is sheltered from the stronger external radiation field coming from the west. We show that photodissociation on the illuminated side can give rise to an asymmetric methanol distribution, but that the stark contrast observed in H-MM1 is hard to explain without assuming enhanced desorption on the shaded side. The region of the brightest emission has a wavy structure that rolls up at one end. This is the signature of Kelvin-Helmholtz instability occurring in sheared flows. We suggest that in this zone, methanol and sulphur are released as a result of grain-grain collisions induced by shear vorticity.

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.

The CS molecule is known to be absorbed onto dust in the cold and dense conditions, causing it to get significantly depleted in the central region of cores. This study is aimed to investigate the depletion of the CS molecule using the optically thin C$^{34}$S molecular line observations. We mapped five prestellar cores, L1544, L1552, L1689B, L694-2, and L1197 using two molecular lines, C$^{34}$S $(J=2-1)$ and N$_2$H$^+$ $(J=1-0)$ with the NRO 45-m telescope, doubling the number of cores where the CS depletion was probed using C$^{34}$S. In most of our targets, the distribution of C$^{34}$S emission shows features that suggest that the CS molecule is generally depleted in the center of the prestellar cores. The radial profile of the CS abundance with respect to H$_2$ directly measured from the CS emission and the Herschel dust emission indicates that the CS molecule is depleted by a factor of $\sim$3 toward the central regions of the cores with respect to their outer regions. The degree of the depletion is found to be even more enhanced by an order of magnitude when the contaminating effect introduced by the presence of CS molecules in the surrounding envelope that lie along the line-of-sight is removed. Except for L1197 which is classified as relatively the least evolved core in our targets based on its observed physical parameters, we found that the remaining four prestellar cores are suffering from significant CS depletion at their central region regardless of the relative difference in their evolutionary status.

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

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

We present initial results from the K-band focal plane array Examinations of Young STellar Object Natal Environments (KEYSTONE) survey, a large project on the 100-m Green Bank Telescope mapping ammonia emission across eleven giant molecular clouds at distances of $0.9-3.0$ kpc (Cygnus X North, Cygnus X South, M16, M17, MonR1, MonR2, NGC2264, NGC7538, Rosette, W3, and W48). This data release includes the NH$_3$ (1,1) and (2,2) maps for each cloud, which are modeled to produce maps of kinetic temperature, centroid velocity, velocity dispersion, and ammonia column density. Median cloud kinetic temperatures range from $11.4\pm2.2$ K in the coldest cloud (MonR1) to $23.0\pm6.5$ K in the warmest cloud (M17). Using dendrograms on the NH$_3$ (1,1) integrated intensity maps, we identify 856 dense gas clumps across the eleven clouds. Depending on the cloud observed, $40-100\%$ of the clumps are aligned spatially with filaments identified in H$_2$ column density maps derived from SED-fitting of dust continuum emission. A virial analysis reveals that 523 of the 835 clumps ($\sim63\%$) with mass estimates are bound by gravity alone. We find no significant difference between the virial parameter distributions for clumps aligned with the dust-continuum filaments and those unaligned with filaments. In some clouds, however, hubs or ridges of dense gas with unusually high mass and low virial parameters are located within a single filament or at the intersection of multiple filaments. These hubs and ridges tend to host water maser emission, multiple 70$\mu$m-detected protostars, and have masses and radii above an empirical threshold for forming massive stars.

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

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

"TRAO FUNS" is a project to survey Gould Belt's clouds in molecular lines. This paper presents its first results on the central region of the California molecular cloud, L1478. We performed On-The-Fly mapping observations using the Taedeok Radio Astronomy Observatory (TRAO) 14m single dish telescope equipped with a 16 multi-beam array covering $\sim$1.0 square degree area of this region using C$^{18}$O (1-0) mainly tracing low density cloud and about 460 square arcminute area using N$_{2}$H$^{+}$ (1-0) mainly tracing dense cores. CS (2-1) and SO $(3_{2}-2_{1})$ were also used simultaneously to map $\sim$440 square arcminute area of this region. We identified 10 filaments by applying the dendrogram technique to the C$^{18}$O data-cube and 8 dense N$_{2}$H$^{+}$ cores by using {\sc FellWalker}. Basic physical properties of filaments such as mass, length, width, velocity field, and velocity dispersion are derived. It is found that L1478 consists of several filaments with slightly different velocities. Especially the filaments which are supercritical are found to contain dense cores detected in N$_{2}$H$^{+}$. Comparison of non-thermal velocity dispersions derived from C$^{18}$O and N$_{2}$H$^{+}$ for the filaments and dense cores indicates that some of dense cores share similar kinematics with those of the surrounding filaments while several dense cores have different kinematics with those of their filaments. This suggests that the formation mechanism of dense cores and filaments can be different in individual filaments depending on their morphologies and environments.

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

We develop a new ''core field structure'' (CFS) model to predict the magnetic field strength and magnetic field fluctuation profile of dense cores using gas kinematics. We use spatially resolved observations of the nonthermal velocity dispersion from the Green Bank Ammonia survey along with column density maps from SCUBA-2 to estimate the magnetic field strength across seven dense cores located in the L1688 region of Ophiuchus. The CFS model predicts the profile of the relative field fluctuation, which is related to the observable dispersion in direction of the polarization vectors. Within the context of our model we find that all the cores have a transcritical mass-to-flux ratio.

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

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

A model of magnetic field structure is presented to help test the prevalence of flux freezing in star-forming clouds of various shapes, orientations, and degrees of central concentration, and to estimate their magnetic field strength. The model is based on weak-field flux freezing in centrally condensed Plummer spheres and spheroids of oblate and prolate shape. For a spheroid of given density contrast, aspect ratio, and inclination, the model estimates the local field strength and direction, and the global field pattern of hourglass shape. Comparisons with a polarization simulation indicate typical angle agreement within 1 - 10 degrees. Scalable analytic expressions are given to match observed polarization patterns, and to provide inputs to radiative transfer codes for more accurate predictions. The model may apply to polarization observations of dense cores, elongated filamentary clouds, and magnetized circumstellar disks.

This paper presents our observational attempts to precisely measure the central mass of a proto-brown dwarf candidate, L328-IRS, in order to investigate whether L328-IRS is in the substellar mass regime. Observations were made for the central region of L328-IRS with the dust continuum and CO isotopologue line emission at ALMA band 6, discovering the detailed outflow activities and a deconvolved disk structure of a size of $\sim$87 AU $\times$ 37 AU. We investigated the rotational velocities as a function of the disk radius, finding that its motions between 130 AU and 60 AU are partially fitted with a Keplerian orbit by a stellar object of $\sim0.30 ~\rm\,M_\odot$, while the motions within 60 AU do not follow any Keplerian orbit at all. This makes it difficult to lead a reliable estimation of the mass of L328-IRS. Nonetheless, our ALMA observations were useful enough to well constrain the inclination angle of the outflow cavity of L328-IRS as $\sim 66^{\circ}$, enabling us to better determine the mass accretion rate of $\sim 8.9\times 10^{-7} \rm\,M_\odot~yr^{-1}$. From assumptions that the internal luminosity of L328-IRS is mostly due to this mass accretion process in the disk, or that L328-IRS has mostly accumulated the mass through this constant accretion rate during its outflow activity, its mass was estimated to be $\sim 0.012 - 0.023\rm\,M_\odot$, suggesting L328-IRS to be a substellar object. However, we leave our identification of L328-IRS as a proto-brown dwarf to be tentative because of various uncertainties especially regarding the mass accretion rate.

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

Filamentary structures are ubiquitous in high-mass star-forming molecular clouds. Their relation with high-mass star formation is still to be understood. Here we report interferometric observations toward 8 filamentary high-mass star-forming clouds. A total of 50 dense cores are identified in these clouds, most of which present signatures of high-mass star formation. Five of them are not associated with any star formation indicators, hence are prestellar core candidates. Evolutionary phases of these cores and their linewidths, temperatures, NH$_3$ abundances, and virial parameters are found to be correlated. In a sub-sample of 4 morphologically well-defined filaments, we find that their fragmentation can not be solely explained by thermal or turbulence pressure support. We also investigate distributions of gas temperatures and non-thermal motions along the filaments, and find a spatial correlation between non-thermal linewidths and star formation activities. We find evidence of gas flows along these filaments, and derive an accretion rate along filaments of $\sim$10$^{-4}$ M$_\odot$ yr$^{-1}$. These results suggest a strong relationship between massive filaments and high-mass star formation, through i) filamentary fragmentation in very early evolutionary phases to form dense cores, ii) accretion flows along filaments that are important for the growth of dense cores and protostars, and iii) enhancement of non-thermal motion in the filaments by the feedback or accretion during star formation.

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

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

Understanding the early stages of star formation is a research field of ongoing development, both theoretically and observationally. In this context, molecular data have been continuously providing observational constraints on the gas dynamics at different excitation conditions and depths in the sources. We have investigated the Barnard 59 core, the only active site of star formation in the Pipe Nebula, to achieve a comprehensive view of the kinematic properties of the source. These information were derived by simultaneously fitting ammonia inversion transition lines (1,1) and (2,2). Our analysis unveils the imprint of protostellar feedback, such as increasing line widths, temperature and turbulent motions in our molecular data. Combined with complementary observations of dust thermal emission, we estimate that the core is gravitationally bound following a virial analysis. If the core is not contracting, another source of internal pressure, most likely the magnetic field, is supporting it against gravitational collapse and limits its star formation efficiency.

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 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 the results of our investigation of the star-forming potential in the Perseus star-forming complex. We build on previous starless core, protostellar core, and young stellar object (YSO) catalogs from Spitzer, Herschel, and SCUBA observations in the literature. We place the cores and YSOs within seven star-forming clumps based on column densities greater than 5x10^21 cm^-2. We calculate the mean density and free-fall time for 69 starless cores as 5.55x10^-19 gcm^-3 and 0.1 Myr,respectively, and we estimate the star formation rate for the near future as 150 Msun Myr^-1. According to Bonnor Ebert stability analysis, we find that majority of starless cores in Perseus are unstable. Broadly, these cores can collapse to form the next generation of stars. We found a relation between starless cores and YSOs, where the numbers of young protostars (Class 0 + Class I) are similar to the numbers of starless cores. This similarity, which shows a one-to-one relation, suggests that these starless cores may form the next generation of stars with approximately the same formation rate as the current generation, as identified by the Class 0 and Class I protostars. It follows that if such a relation between starless cores and any YSO stage exists, the SFR values of these two populations must be nearly constant. In brief, we propose that this one-to-one relation is an important factor in better understanding the star formation process within a cloud.

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.

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

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.

New models of star-forming filamentary clouds are presented, to quantify their properties and to predict their evolution. These 2D axisymmetric models describe filaments having no core, one low-mass core, and one cluster-forming core. They are based on Plummer-like cylinders and spheroids, bounded by a constant-density surface of finite extent. In contrast to 1D Plummer-like models, they have specific values of length and mass, they approximate observed column density maps, and their distributions of column density (N-pdfs) are pole-free. Each model can estimate the star-forming potential of a core-filament system, by identifying the zone of gas dense enough to form low-mass stars, and by counting the number of enclosed thermal Jeans masses. This analysis suggests that the Musca Center filament may be near the start of its star-forming life, with enough dense gas to make its first ~3 protostars, while the Coronet filament is near the midpoint of its star formation, with enough dense gas to add ~ 8 protostars to its ~20 known stars. In contrast L43 appears near the end of its star-forming life, since it lacks enough dense gas to add any new protostars to the 2 YSOs already known.

Ammonia and its deuterated isotopologues probe physical conditions in dense molecular cloud cores. With the aim of testing the current understanding of the spin-state chemistry of these molecules, we observed spectral lines of NH3, NH2D, NHD2, ND3, and N2D+ towards a dense, starless core in Ophiuchus with the APEX, GBT, and IRAM 30-m telescopes. The observations were interpreted using a gas-grain chemistry model combined with radiative transfer calculations. The chemistry model distinguishes between the different nuclear spin states of light hydrogen molecules, ammonia, and their deuterated forms. High deuterium fractionation ratios with NH2D/NH3=0.4, NHD2/NH2D=0.2, and ND3/NHD2=0.06 were found in the core. The observed ortho/para ratios of NH2D and NHD2 are close to the corresponding nuclear spin statistical weights. The chemistry model can approximately reproduce the observed abundances, but predicts uniformly too low ortho/para-NH2D, and too large ortho/para-NHD2 ratios. The longevity of N2H+ and NH3 in dense gas, which is prerequisite to their strong deuteration, can be attributed to the chemical inertia of N2 on grain surfaces. The discrepancies between the chemistry model and the observations are likely to be caused by the fact that the model assumes complete scrambling in principal gas-phase deuteration reactions of ammonia, which means that all the nuclei are mixed in reactive collisions. If, instead, these reactions occur through proton hop/hydrogen abstraction processes, statistical spin ratios are to be expected. The present results suggest that while the deuteration of ammonia changes with physical conditions and time, the nuclear spin ratios of ammonia isotopologues do not probe the evolutionary stage of a cloud.

Filamentary structures are ubiquitous in molecular clouds, and have been recently argued to play an important role in regulating the size and mass of embedded clumps through fragmentation and mass accretion. Here, we reveal the dynamical state and fragmentation of filamentary molecular gas associated with the Serpens South protocluster through analysis of wide (~4 x 4 pc) observations of NH3 (1,1) and (2,2) inversion transitions with the Green Bank Telescope. Detailed modeling of the NH3 lines reveals that the kinematics of the cluster and surrounding filaments are complex. We identify hierarchical structure using a dendrogram analysis of the NH3 emission. The distance between neighbour structures that are embedded within the same parent structure is generally greater than expected from a spherical Jeans analysis, and is in better agreement with cylindrical fragmentation models. The NH3 line width-size relation is flat, and average gas motions are sub- or trans-sonic over all physical scales observed. Subsonic regions extend far beyond the typical 0.1 pc scale previously identified in star-forming cores. As a result, we find a strong trend of decreasing virial parameter with increasing structure mass in Serpens South. Extremely low virial parameters on the largest scales probed by our data suggest that the previously observed, ordered magnetic field is insufficient to support the region against collapse, in agreement with large radial infall motions previously measured toward some of the filaments. A more complex magnetic field configuration in the dense gas, however, may be able to support the filaments.

Although surveys of infall motions in dense cores have been carried out for years, few surveys have focused on mapping infall across cores using multiple spectral line observations. To fill this gap, we present IRAM 30-m Telescope maps of N2H+(1-0), DCO+(2-1), DCO+(3-2), and HCO+(3-2) emission towards two prestellar cores (L492 and L694-2) and one protostellar core (L1521F). We find that the measured infall velocity varies with position across each core and choice of molecular line, likely as a result of radial variations in core chemistry and dynamics. Line-of-sight infall speeds estimated from DCO+(2-1) line profiles can decrease by 40-50 m/s when observing at a radial offset >= 0.04 pc from the core's dust continuum emission peak. Median infall speeds calculated from all observed positions across a core can also vary by as much as 65 m/s depending on the transition. These results show that while single-pointing, single-transition surveys of core infall velocities may be good indicators of whether a core is either contracting or expanding, the magnitude of the velocities they measure are significantly impacted by the choice of molecular line, proximity to the core center, and core evolutionary state.

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

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.

Column density (N) PDFs serve as a powerful tool to characterize the physical processes that influence the structure of molecular clouds. Star-forming clouds can best be characterized by lognormal PDFs for the lower N range and a power-law tail for higher N, commonly attributed to turbulence and self-gravity and/or pressure, respectively. We report here on PDFs obtained from observations of 12CO, 13CO, C18O, CS, and N2H+ in the Cygnus X North region and compare to a PDF derived from dust observations with the Herschel satellite. The PDF of 12CO is lognormal for Av~1-30, but is cut for higher Av due to optical depth effects. The PDFs of C18O and 13CO are mostly lognormal up for Av~1-15, followed by excess up to Av~40. Above that value, all CO PDFs drop, most likely due to depletion. The high density tracers CS and N2H+ exhibit only a power law distribution between Av~15 and 400, respectively. The PDF from dust is lognormal for Av~3-15 and has a power-law tail up to Av~500. Absolute values for the molecular line column densities are, however, rather uncertain due to abundance and excitation temperature variations. Taken the dust PDF face value, we 'calibrate' the molecular line PDF of CS to the one of the dust and determined an abundance [CS]/[H2] of 10^-9. The slopes of the power-law tails of the CS, N2H+, and dust PDFs are consistent with free-fall collapse of filaments and clumps. A quasi static configuration of filaments and clumps can possibly also account for the observed N-PDFs, as long as they have a sufficiently condensed density structure and external ram pressure by gas accretion is provided. The somehow flatter slopes of N2H+ and CS can reflect an abundance change and/or subthermal excitation at low column densities.

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.

We analyze the spatial distribution of dusty young stellar objects (YSOs) identified in the Spitzer Survey of the Orion Molecular clouds, augmenting these data with Chandra X-ray observations to correct for incompleteness in dense clustered regions. We also devise a scheme to correct for spatially varying incompleteness when X-ray data are not available. The local surface densities of the YSOs range from 1 pc$^{-2}$ to over 10,000 pc$^{-2}$, with protostars tending to be in higher density regions. This range of densities is similar to other surveyed molecular clouds with clusters, but broader than clouds without clusters. By identifying clusters and groups as continuous regions with surface densities $\ge10$ pc$^{-2}$, we find that 59% of the YSOs are in the largest cluster, the Orion Nebular Cluster (ONC), while 13% of the YSOs are found in a distributed population. A lower fraction of protostars in the distributed population is evidence that it is somewhat older than the groups and clusters. An examination of the structural properties of the clusters and groups show that the peak surface densities of the clusters increase approximately linearly with the number of members. Furthermore, all clusters with more than 70 members exhibit asymmetric and/or highly elongated structures. The ONC becomes azimuthally symmetric in the inner 0.1 pc, suggesting that the cluster is only $\sim 2$ Myr in age. We find the star formation efficiency (SFE) of the Orion B cloud is unusually low, and that the SFEs of individual groups and clusters are an order of magnitude higher than those of the clouds. Finally, we discuss the relationship between the young low mass stars in the Orion clouds and the Orion OB 1 association, and we determine upper limits to the fraction of disks that may be affected by UV radiation from OB stars or by dynamical interactions in dense, clustered regions.

We present mid-IR (19 - 37 microns) imaging observations of S106 from SOFIA/FORCAST, complemented with IR observations from Spitzer/IRAC (3.6 - 8.0 microns), IRTF/MIRLIN (11.3 and 12.5 microns), and Herschel/PACS (70 and 160 microns). We use these observations, observations in the literature, and radiation transfer modeling to study the heating and composition of the warm (~ 100 K) dust in the region. The dust is heated radiatively by the source S106 IR, with little contributions from grain-electron collisions and Ly-alpha radiation. The dust luminosity is >~ (9.02 +/- 1.01) x 10^4 L_sun, consistent with heating by a mid- to late-type O star. We find a temperature gradient (~ 75 - 107 K) in the lobes, which is consistent with a dusty equatorial geometry around S106 IR. Furthermore, the SOFIA observations resolve several cool (~ 65 - 70 K) lanes and pockets of warmer (~ 75 - 90 K) dust in the ionization shadow, indicating that the environment is fragmented. We model the dust mass as a composition of amorphous silicates, amorphous carbon, big grains, very small grains, and PAHs. We present the relative abundances of each grain component for several locations in S106.

Observations of the isolated globule B335 with ALMA have yielded absorption features against the continuum that are redshifted from the systemic velocity in both HCN and HCO$^+$ lines. These features provide unambiguous evidence for infall toward a central luminosity source. Previously developed models of inside-out collapse can match the observed line profiles of HCN and HCO$^+$ averaged over the central 50 AU. At the new distance of 100 pc, the inferred infall radius is 0.012 pc, the mass infall rate is $3 \times 10^{-6}$ solar masses per year, the age is 50,000 years, and the accumulated mass in the central zone is 0.15 solar masses, most of which must be in the star or in parts of a disk that are opaque at 0.8 mm. The continuum detection indicates an optically thin mass (gas and dust) of only $7.5\times 10^{-4}$ solar masses in the central region, consistent with only a very small disk mass.

This paper gives a new way to diagnose the star-forming potential of a molecular cloud region from the probability density function of its column density (N-pdf). It gives expressions for the column density and mass profiles of a symmetric filament having the same N-pdf as a filamentary region. The central concentration of this characteristic filament can distinguish regions and can quantify their fertility for star formation. Profiles are calculated for N-pdfs which are pure lognormal, pure power law, or a combination. In relation to models of singular polytropic cylinders, characteristic filaments can be unbound, bound, or collapsing depending on their central concentration. Such filamentary models of the dynamical state of N-pdf gas are more relevant to star-forming regions than are models of spherical collapse. The star formation fertility of a bound or collapsing filament is quantified by its mean mass accretion rate when in radial free fall. For a given mass per length, the fertility increases with the filament mean column density and with its initial concentration. In selected regions the fertility of their characteristic filaments increases with the level of star formation.

We present Spitzer 3.6, 4.5, 5.8, 8.0, and 24 micron images of the Mon OB1 East giant molecular cloud, which contains the young star forming region NGC 2264, as well as more extended star formation. With Spitzer data and 2MASS photometry, we identify and classify young stellar objects (YSOs) with dusty circumstellar disks and/or envelopes in Mon OB1 East by their infrared-excess emission and study their distribution with respect to cloud material. We find a correlation between the local surface density of YSOs and column density of molecular gas as traced by dust extinction that is roughly described as a power law in these quantities. NGC 2264 follows a power law index of ~2.7, exhibiting a large YSO surface density for a given gas column density. Outside of NGC 2264 where the surface density of YSOs is lower, the power law is shallower and the region exhibits a larger gas column density for a YSO surface density, suggesting the star formation is more recent. In order to measure the fraction of cloud members with circumstellar disks/envelopes, we estimate the number of diskless pre-main sequence stars by statistical removal of background star detections. We find that the disk fraction of the NGC 2264 region is 45%, while the surrounding more distributed regions show a disk fraction of 19%. This may be explained by the presence an older, more dispersed population of stars. In total, the Spitzer observations provide evidence for heterogenous, non-coeval star formation throughout the Mon OB1 cloud.

We report on Submillimeter Array observations of the 870 micron continuum and CO(3-2), 13CO(2-1) and C18O(2-1) line emission of a faint object, SMM2E, near the driving source of the HH797 outflow in the IC348 cluster. The continuum emission shows an unresolved source for which we estimate a mass of gas and dust of 30 Mjup, and the CO(3-2) line reveals a compact bipolar outflow centred on SMM2E, and barely seen also in 13CO(2-1). In addition, C18O(2-1) emission reveals hints of a possible rotating envelope/disk perpendicular to the outflow, for which we infer a dynamical mass of ~16 Mjup. In order to further constrain the accreted mass of the object, we gathered data from Spitzer, Herschel, and new and archive submillimetre observations, and built the Spectral Energy Distribution (SED). The SED can be fitted with one single modified black-body from 70 micron down to 2.1 cm, using a dust temperature of ~24 K, a dust emissivity index of 0.8, and an envelope mass of ~35 Mjup. The bolometric luminosity is 0.10 Lsun, and the bolometric temperature is 35 K. Thus, SMM2E is comparable to the known Class 0 objects in the stellar domain. An estimate of the final mass indicates that SMM2E will most likely remain substellar, and the SMM2E outflow force matches the trend with luminosity known for young stellar objects. Thus, SMM2E constitutes an excellent example of a Class 0 proto-brown dwarf candidate which forms as a scaled-down version of low-mass stars. Finally, SMM2E seems to be part of a wide (~2400 AU) multiple system of Class 0 sources.

We consider the interaction between a marginally stable Bonnor-Ebert (BE) sphere and the surrounding ambient medium. In particular, we explore how the infall from an evolving ambient medium can trigger collapse of the sphere using 3-D adaptive mesh refinement simulations. We find the resulting collapse dynamics to vary considerably with ambient density. In the highest ambient density cases, infalling material drives a strong compression wave into the cloud. It is the propagation of this wave through the cloud interior that triggers the subsequent collapse. For lower ambient densities, we find the main trigger of collapse to be a quasistatic adjustment of the BE sphere to gravitational settling of the ambient gas. In all cases, we find that the classic "outside-in" collapse mode for super-critical BE spheres is recovered before a protostar (i.e. sink particle) forms. Our work supports scenarios in which BE dynamics naturally begins with either a compression-wave or infall dominated phase, and only later assumes the usual outside-in collapse behavior.

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

In an improved model of protostar mass functions (PMFs), protostars gain mass from isothermal cores in turbulent clumps. Their mass accretion rate is similar to Shu accretion at low mass, and to reduced Bondi accretion at high mass. Accretion durations follow a simple expression in which higher-mass protostars accrete for longer times. These times are set by ejections, stellar feedback, and gravitational competition, which terminate accretion and reduce its efficiency. The mass scale is the mass of a critically stable isothermal core. In steady state, the PMF approaches a power law at high mass due to competition between clump accretion and accretion stopping. The power law exponent is the ratio of the time scales of accretion and accretion stopping. The luminosity function (PLF) peaks near 1 L_Sun, due to inefficient accretion of core gas. Models fit observed PLFs in four large embedded clusters. These indicate that their underlying PMFs may be top-heavy compared to the IMF, depending on the model of protostar radius.

We have detected bright HC7N J = 21-20 emission toward multiple locations in the Serpens South cluster-forming region using the K-Band Focal Plane Array at the Robert C. Byrd Green Bank Telescope. HC7N is seen primarily toward cold filamentary structures that have yet to form stars, largely avoiding the dense gas associated with small protostellar groups and the main central cluster of Serpens South. Where detected, the HC7N abundances are similar to those found in other nearby star forming regions. Toward some HC7N `clumps', we find consistent variations in the line centroids relative to NH3 (1,1) emission, as well as systematic increases in the HC7N non-thermal line widths, which we argue reveal infall motions onto dense filaments within Serpens South with minimum mass accretion rates of M ~ 2-5 M_sun Myr^-1. The relative abundance of NH3 to HC7N suggests that the HC7N is tracing gas that has been at densities n ~ 10^4 cm^-3, for timescales t < 1-2 x 10^5 yr. Since HC7N emission peaks are rarely co-located with those of either NH3 or continuum, it is likely that Serpens South is not particularly remarkable in its abundance of HC7N, but instead the serendipitous mapping of HC7N simultaneously with NH3 has allowed us to detect HC7N at low abundances in regions where it otherwise may not have been looked for. This result extends the known star-forming regions containing significant HC7N emission from typically quiescent regions, like the Taurus molecular cloud, to more complex, active environments.

We present an XMM-Newton survey of the part of Orion A cloud south of the Orion Nebula. This survey includes the Lynds 1641 (L1641) dark cloud, a region of the Orion A cloud with very few massive stars and hence a relatively low ambient UV flux, and the region around the O9 III star Iota Orionis. In addition to proprietary data, we used archival XMM data of the Orion Nebula Cluster (ONC) to extend our analysis to a major fraction of the Orion A cloud. We have detected 1060 X-ray sources in L1641 and Iota Ori region. About 94% of the sources have 2MASS & Spitzer counterparts, 204 and 23 being Class II and Class I or protostars objects, respectively. In addition, we have identified 489 X-ray sources as counterparts to Class III candidates, given they are bright in X-rays and appear as normal photospheres at mid-IR wavelengths. The remaining 205 X-ray sources are likely distant AGNs or other galactic sources not related to Orion A. We find that Class III candidates appear more concentrated in two main clusters in L1641. The first cluster of Class III stars is found toward the northern part of L1641, concentrated around Iota Ori. The stars in this cluster are more evolved than those in the Orion Nebula. We estimate a distance of 300-320 pc for this cluster and thus it is closer than the Orion A cloud. Another cluster rich in Class III stars is located in L1641 South and appears to be a slightly older cluster embedded in the Orion A cloud. Furthermore, other evolved Class III stars are found north of the ONC toward NGC 1977.

One puzzle in understanding how stars form in clusters is the source of mass -- is all of the mass in place before the first stars are born, or is there an extended period when the cluster accretes material which can continuously fuel the star formation process? We use a multi-line spectral survey of the southern filament associated with the Serpens South embedded cluster-forming region in order to determine if mass is accreting from the filament onto the cluster, and whether the accretion rate is significant. Our analysis suggests that material is flowing along the filament's long axis at a rate of ~30Msol/Myr (inferred from the N2H+ velocity gradient along the filament), and radially contracting onto the filament at ~130Msol/Myr (inferred from HNC self-absorption). These accretion rates are sufficient to supply mass to the central cluster at a similar rate to the current star formation rate in the cluster. Filamentary accretion flows may therefore be very important in the ongoing evolution of this cluster.

Mass functions of starless dense cores (CMFs) may arise from contraction and dispersal of core-forming filaments. In an illustrative model, a filament contracts radially by self-gravity, increasing the mass of its cores. During this contraction, FUV photoevaporation and ablation by shocks and winds disperse filament gas and limit core growth. The stopping times of core growth are described by a waiting-time distribution. The initial filament column density profile and the resulting CMF each match recent Herschel observations in detail. Then low-mass cores have short growth ages and arise from the innermost filament gas, while massive cores have long growth ages and draw from more extended filament gas. The model fits the initial density profile and CMF best for mean core density 2 10^4 cm^-3 and filament dispersal time scale 0.5 Myr. Then the typical core mass, radius, mean column density, and contraction speed are respectively 0.8 solar masses, 0.06 pc, 6 10^21 cm^-2, and 0.07 km s^-1, also in accord with observed values.

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

We present a survey of the Orion A and B molecular clouds undertaken with the IRAC and MIPS instruments onboard Spitzer. In total, five distinct fields were mapped covering 9 sq. degrees in five mid-IR bands spanning 3-24 microns. The survey includes the Orion Nebula Cluster, the Lynds 1641, 1630 and 1622 dark clouds, and the NGC 2023, 2024, 2068 and 2071 nebulae. These data are merged with the 2MASS point source catalog to generate a catalog of eight band photometry. We identify 3479 dusty young stellar objects (YSOs) in the Orion molecular clouds by searching for point sources with mid-IR colors indicative of reprocessed light from dusty disks or infalling envelopes. The YSOs are subsequently classified on the basis of their mid-IR colors and their spatial distributions are presented. We classify 2991 of the YSOs as pre-main sequence stars with disks and 488 as likely protostars. Most of the sources were observed with IRAC in 2-3 epochs over 6 months; we search for variability between the epochs by looking for correlated variability in the 3.6 and 4.5 micron bands. We find that 50% of the dusty YSOs show variability. The variations are typically small (0.2 mag.) with the protostars showing a higher incidence of variability and larger variations. The observed correlations between the 3.6, 4.5, 5.8 and 8 micron variability suggests that we are observing variations in the heating of the inner disk due to changes in the accretion luminosity or rotating accretion hot spots.

We present Spitzer IRAC and MIPS observations of the star-forming region containing intermediate-mass young stellar object (YSO) AFGL 490. We supplement these data with near-IR 2MASS photometry and with deep SQIID observations off the central high extinction region. We have more than doubled the known membership of this region to 57 Class I and 303 Class II YSOs via the combined 1-24 um photometric catalog derived from these data. We construct and analyze the minimum spanning tree of their projected positions, isolating one locally over-dense cluster core containing 219 YSOs (60.8% of the region's members). We find this cluster core to be larger yet less dense than similarly analyzed clusters. Although the structure of this cluster core appears irregular, we demonstrate that the parsec-scale surface densities of both YSOs and gas are correlated with a power law slope of 2.8, as found for other similarly analyzed nearby molecular clouds. We also explore the mass segregation implications of AFGL 490's offset from the center of its core, finding that it has no apparent preferential central position relative to the low-mass members.

We identify protostars in Spitzer surveys of nine star-forming molecular clouds within 1 kpc: Serpens, Perseus, Ophiuchus, Chamaeleon, Lupus, Taurus, Orion, Cep OB3, and Mon R2, which combined host over 700 protostar candidates. Our diverse cloud sample allows us to compare protostar luminosity functions in these varied environments. We combine photometry from 2MASS J, H, and Ks bands and Spitzer IRAC and MIPS 24 micron bands to create 1 - 24 micron spectral energy distributions (SEDs). Using protostars from the c2d survey with well-determined bolometric luminosities (Lbol), we derive a relationship between Lbol, L_MIR (integrated from 1 - 24 microns), and SED slope. Estimations of Lbol for protostar candidates are combined to create luminosity functions for each cloud. Contamination due to edge-on disks, reddened Class II sources, and galaxies is estimated and removed from the luminosity functions. We find that luminosity functions for high mass star forming clouds peak near 1 Lsun and show a tail extending toward luminosities above 100 Lsun. The luminosity functions of the low mass star forming clouds do not exhibit a common peak, however the combined luminosity function of these regions peaks below 1 Lsun. Finally, we examine the luminosity functions as a function of the local surface density of YSOs. In the Orion molecular cloud, we find a significant difference between the luminosity functions of protostars in regions of high and low stellar density, the former of which is biased toward more luminous sources. This may be the result of primordial mass segregation, although this interpretation is not unique. We compare our luminosity functions to those predicted by models and find that our observed luminosity functions are best matched by models which invoke competitive accretion, although we do not find strong agreement of the high mass star forming clouds with any of the models.

A model of protostar mass and luminosity evolution in clusters gives new estimates of cluster age, protostar birthrate, accretion rate and mean accretion time. The model assumes constant protostar birthrate, core-clump accretion, and equally likely accretion stopping. Its parameters are set to reproduce the initial mass function, and to match protostar luminosity distributions in nearby star-forming regions. It obtains cluster ages and birthrates from the observed numbers of protostars and pre-main sequence (PMS) stars, and from the modal value of the protostar luminosity. In 31 embedded clusters and complexes the global cluster age is 1-3 Myr, matching available estimates based on optical spectroscopy and evolutionary tracks. This method of age estimation is simpler than optical spectroscopy, and is more useful for young embedded clusters where optical spectrocopy is not possible. In the youngest clusters, the protostar fraction decreases outward from the densest gas, indicating that the local star-forming age increases outward from a few 0.1 Myr in small protostar-dominated zones to a few Myr in large PMS-dominated zones.

We present high spatial (<300 AU) and spectral (0.07 km/s) resolution Submillimeter Array observations of the dense starless cluster core Oph A-N6, in the 1 mm dust continuum and the 3-2 line of N2H+ and N2D+. The dust continuum observations reveal a compact source not seen in single-dish observations, of size ~1000 AU and mass 0.005-0.01 M\odot. The combined line and single-dish observations reveal a core of size 3000 \times 1400 AU elongated in a NW-SE direction, with almost no variation in either line width or line center velocity across the map, and very small non-thermal motions. The deuterium fraction has a peak value of ~0.15 and is >0.05 over much of the core. The N2H+ column density profile across the major axis of Oph A-N6 is well represented by an isothermal cylinder, with temperature 20 K, peak density 7.1 \times 10^6 cm^{-3}, and N2H+ abundance 2.7 \times 10^{-10}. The mass of Oph A-N6 is estimated to be 0.29 M\odot, compared to a value of 0.18 M\odot from the isothermal cylinder analysis, and 0.63 M\odot for the critical mass for fragmentation of an isothermal cylinder. Compared to isolated low-mass cores, Oph A-N6 shows similar narrow line widths and small velocity variation, with a deuterium fraction similar to "evolved" dense cores. It is significantly smaller than isolated cores, with larger peak column and volume density. The available evidence suggests Oph A-N6 has formed through the fragmentation of the Oph A filament and is the precursor to a low-mass star. The dust continuum emission suggests it may already have begun to form a star.

We analyze high quality, complete stellar catalogs for four young (roughly 1 Myr) and nearby (within ~300 pc) star-forming regions: Taurus, Lupus3, ChaI, and IC348, which have been previously shown to have stellar groups whose properties are similar to those of larger clusters such as the ONC. We find that stars at higher stellar surface densities within a region or belonging to groups tend to have a relative excess of more massive stars, over a wide range of masses. We find statistically significant evidence for this result in Taurus and IC348 as well as the ONC. These differences correspond to having typically a ~10 - 20% higher mean mass in the more clustered environment. Stars in ChaI show no evidence for a trend with either surface density or grouped status, and there are too few stars in Lupus3 to make any definitive interpretation. Models of clustered star formation do not typically extend to sufficiently low masses or small group sizes in order for their predictions to be tested but our results suggest that this regime is important to consider.

A model of core-clump accretion with equally likely stopping describes star formation in the dense parts of clusters, where models of isolated collapsing cores may not apply. Each core accretes at a constant rate onto its protostar, while the surrounding clump gas accretes as a power of protostar mass. Short accretion flows resemble Shu accretion, and make low-mass stars. Long flows resemble reduced Bondi accretion and make massive stars. Accretion stops due to environmental processes of dynamical ejection, gravitational competition, and gas dispersal by stellar feedback, independent of initial core structure. The model matches the field star IMF from 0.01 to more than 10 solar masses. The core accretion rate and the mean accretion duration set the peak of the IMF, independent of the local Jeans mass. Massive protostars require the longest accretion durations, up to 0.5 Myr. The maximum protostar luminosity in a cluster indicates the mass and age of its oldest protostar. The distribution of protostar luminosities matches those in active star-forming regions if protostars have a constant birthrate but not if their births are coeval. For constant birthrate, the ratio of YSOs to protostars indicates the star-forming age of a cluster, typically ~1 Myr. The protostar accretion luminosity is typically less than its steady spherical value by a factor of ~2, consistent with models of episodic disk accretion.

We report the discovery and characterization of a power law correlation between the local surface densities of Spitzer-identified, dusty young stellar objects and the column density of gas (as traced by near-IR extinction) in eight molecular clouds within 1 kpc and with 100 or more known YSOs. This correlation, which appears in data smoothed over size scales of ~1 pc, varies in quality from cloud to cloud; those clouds with tight correlations, MonR2 and Ophiuchus, are fit with power laws of slope 2.67 and 1.87, respectively. The spread in the correlation is attributed primarily to local gas disruption by stars that formed there or to the presence of very young sub-regions at the onset of star formation. We explore the ratio of the number of Class II to Class I sources, a proxy for the star formation age of a region, as a function of gas column density; this analysis reveals a declining Class II to Class I ratio with increasing column density. We show that the observed star-gas correlation is consistent with a star formation law where the star formation rate per area varies with the gas column density squared. We also propose a simple picture of thermal fragmentation of dense gas in an isothermal, self-gravitating layer as an explanation for the power law. Finally, we briefly compare the star gas correlation and its implied star formation law with other recent proposed of star formation laws at similar and larger size scales from nearby star forming regions.

We report the detection of a compact (of order 5 arcsec; about 1800 AU projected size) CO outflow from L1148-IRS. This confirms that this Spitzer source is physically associated with the nearby (about 325 pc) L1148 dense core. Radiative transfer modeling suggests an internal luminosity of 0.08 to 0.13 L_sun. This validates L1148-IRS as a Very Low Luminosity Object (VeLLO; L < 0.1 L_sun). The L1148 dense core has unusually low densities and column densities for a star-forming core. It is difficult to understand how L1148-IRS might have formed under these conditions. Independent of the exact final mass of this VeLLO (which is likely < 0.24 M_sun), L1148-IRS and similar VeLLOs might hold some clues about the isolated formation of brown dwarfs.

New models are presented for star-forming condensations in clusters. In each model, the condensation mass increases linearly with radius on small scales, and more rapidly on large scales, as in "thermal-nonthermal" models. Spherical condensations with this structure form protostars which match the IMF if their infall is subject to equally likely stopping. However such spherical models do not match the filamentary nature of cluster gas, and they are too extended to form protostars having high mass and short spacing. Two hybrid models are presented, which are spherical on small scales and filamentary on large scales. In and around clusters, cores embedded in linear filaments match the elongation of cluster gas, and the central concentration of low-mass stars. In cluster centers, condensations require a low volume filling factor to produce massive stars with short spacing. These may have stellate shape, where cores are nodes of filamentary networks, as seen in some simulations of colliding flows and of collapsing turbulent clumps. A dense configuration of such stellate condensations may be indistinguishable from a clump forming multiple protostars via filamentary flow paths.

This paper discusses the statistics of internal motions in starless dense cores and the relation of these motions to core density and evolution. Four spectral lines from three molecular species are analyzed from single-pointing and mapped observations of several tens of starless cores. Blue asymmetric profiles are dominant, indicating that inward motions are prevalent in sufficiently dense starless cores. These blue profiles are found to be more abundant, and their asymmetry is bluer, at core positions with stronger $\rm N_2H^+$ line emission or higher column density. Thirty three starless cores are classified into four types according to the blue and red shifts of the lines in their molecular line maps. Among these cores, contracting motions dominate: 19 are classified as contracting, 3 as oscillating, 3 as expanding, and 8 as static. Contracting cores have inward motions all over the core with predominance of those motions near the region of peak density. Cores with the bluest asymmetry tend to have greater column density than other cores and all five cores with peak column density $> \rm 6\times 10^{21} cm^{-2}$ are found to be contracting. This suggests that starless cores are likely to have contracting motions if they are sufficiently condensed. Our classification of the starless cores may indicate a sequence of core evolution in the sense that column density increases from static to contracting cores: the static cores in the earliest stage, the expanding and/or the oscillating cores in the next, and the contracting cores in the latest stage.

We investigate the relation of the stellar initial mass function (IMF) and the dense core mass function (CMF), using stellar masses and positions in 14 well-studied young groups. Initial column density maps are computed by replacing each star with a model initial core having the same star formation efficiency (SFE). For each group the SFE, core model, and observational resolution are varied to produce a realistic range of initial maps. A clumpfinding algorithm parses each initial map into derived cores, derived core masses, and a derived CMF. The main result is that projected blending of initial cores causes derived cores to be too few and too massive. The number of derived cores is fewer than the number of initial cores by a mean factor 1.4 in sparse groups and 5 in crowded groups. The mass at the peak of the derived CMF exceeds the mass at the peak of the initial CMF by a mean factor 1.0 in sparse groups and 12.1 in crowded groups. These results imply that in crowded young groups and clusters, the mass distribution of observed cores may not reliably predict the mass distribution of protostars which will form in those cores.

We analyze the masses and spatial distributions of fourteen young stellar groups in Taurus, Lupus3, ChaI, and IC348. These nearby groups, which typically contain 20 to 40 members, have membership catalogs complete to ~0.02 M_sun, and are sufficiently young that their locations should be similar to where they formed. These groups show five properties seen in clusters having many more stars and much greater surface density of stars: (1) a broad range of masses, (2) a concentration of the most massive star towards the centre of the group, (3) an association of the most massive star with a high surface density of lower-mass stars, (4) a correlation of the mass of the most massive star with the total mass of the group, and (5) the distribution of a large fraction of the mass in a small fraction of the stars.

We present 3.6 to 70 {\mu}m Spitzer photometry of 154 weak-line T Tauri stars (WTTS) in the Chamaeleon, Lupus, Ophiuchus and Taurus star formation regions, all of which are within 200 pc of the Sun. For a comparative study, we also include 33 classical T Tauri stars (CTTS) which are located in the same star forming regions. Spitzer sensitivities allow us to robustly detect the photosphere in the IRAC bands (3.6 to 8 {\mu}m) and the 24 {\mu}m MIPS band. In the 70 {\mu}m MIPS band, we are able to detect dust emission brighter than roughly 40 times the photosphere. These observations represent the most sensitive WTTS survey in the mid to far infrared to date, and reveal the frequency of outer disks (r = 3-50 AU) around WTTS. The 70 {\mu}m photometry for half the c2d WTTS sample (the on-cloud objects), which were not included in the earlier papers in this series, Padgett et al. (2006) and Cieza et al. (2007), are presented here for the first time. We find a disk frequency of 19% for on-cloud WTTS, but just 5% for off- cloud WTTS, similar to the value reported in the earlier works. WTTS exhibit spectral energy distributions (SEDs) that are quite diverse, spanning the range from optically thick to optically thin disks. Most disks become more tenuous than Ldisk/L* = 2 x 10^-3 in 2 Myr, and more tenuous than Ldisk/L* = 5 x 10^-4 in 4 Myr.

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

We present N2D+ 3-2 (IRAM) and H2D+ 1_11 - 1_10 and N2H+ 4-3 (JCMT) maps of the small cluster-forming Ophiuchus B2 core in the nearby Ophiuchus molecular cloud. In conjunction with previously published N2H+ 1-0 observations, the N2D+ data reveal the deuterium fractionation in the high density gas across Oph B2. The average deuterium fractionation R_D = N(N2D+)/N(N2H+) ~ 0.03 over Oph B2, with several small scale R_D peaks and a maximum R_D = 0.1. The mean R_D is consistent with previous results in isolated starless and protostellar cores. The column density distributions of both H2D+ and N2D+ show no correlation with total H2 column density. We find, however, an anticorrelation in deuterium fractionation with proximity to the embedded protostars in Oph B2 to distances >= 0.04 pc. Destruction mechanisms for deuterated molecules require gas temperatures greater than those previously determined through NH3 observations of Oph B2 to proceed. We present temperatures calculated for the dense core gas through the equating of non-thermal line widths for molecules (i.e., N2D+ and H2D+) expected to trace the same core regions, but the observed complex line structures in B2 preclude finding a reasonable result in many locations. This method may, however, work well in isolated cores with less complicated velocity structures. Finally, we use R_D and the H2D+ column density across Oph B2 to set a lower limit on the ionization fraction across the core, finding a mean x_e, lim >= few x 10^{-8}. Our results show that care must be taken when using deuterated species as a probe of the physical conditions of dense gas in star-forming regions.

We present Chandra X-ray data of the NGC 1333 embedded cluster, combining these data with existing Chandra data, Sptizer photometry and ground based spectroscopy of both the NGC 1333 & Serpens North clusters to perform a detailed study of the X-ray properties of two of the nearest embedded clusters to the Sun. In NGC 1333, a total of 95 cluster members are detected in X-rays, of which 54 were previously identified with Spitzer. Of the Spitzer sources, we detect 23% of the Class I protostars, 53% of the Flat Spectrum sources, 52% of the Class II, and 50% of the Transition Disk YSOs. Forty-one Class III members of the cluster are identified, bringing the total identified YSO population to 178. The X-ray Luminosity Functions (XLFs) of the NGC 1333 and Serpens clusters are compared to each other and the Orion Nebula Cluster. Based on this comparison, we obtain a new distance for the Serpens cluster of 360+22/-13 pc. The X-ray luminosity was found to depend on the bolometric luminosity as in previous studies of other clusters, and that Lx depends primarily on the stellar surface area. In the NGC 1333 cluster, the Class III sources have a somewhat higher X-ray luminosity for a given surface area. We also find evidence in NGC 1333 for a jump in the X-ray luminosity between spectral types of M0 and K7, we speculate that this may result from the presence of radiative zones in the K-stars. The gas column density vs. extinction in the NGC 1333 was found to be N_H = 0.89 +/- 0.13 x 10^22 A_K, this is lower than expected of the standard ISM but similar to that found previously in the Serpens Cloud Core.

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

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

We measure the mass and size of cloud fragments in several molecular clouds continuously over a wide range of spatial scales (0.05 < r / pc < 3). Based on the recently developed "dendrogram-technique", this characterizes dense cores as well as the enveloping clouds. "Larson's 3rd Law" of constant column density, m(r) = C*r^2, is not well suited to describe the derived mass-size data. Solar neighborhood clouds not forming massive stars (< 10 M_sun; Pipe Nebula, Taurus, Perseus, and Ophiuchus) obey m(r) < 870 M_sun (r / pc)^1.33 . In contrast to this, clouds forming massive stars (Orion A, G10.15$-$0.34, G11.11$-$0.12) do exceed the aforementioned relation. Thus, this limiting mass-size relation may approximate a threshold for the formation of massive stars. Across all clouds, cluster-forming cloud fragments are found to be---at given radius---more massive than fragments devoid of clusters. The cluster-bearing fragments are found to roughly obey a mass-size law m = C*r^1.27 (where the exponent is highly uncertain in any given cloud, but is certainly smaller than 1.5).

Initial conditions for star formation in clusters are estimated for protostars whose masses follow the initial mass function (IMF) from 0.05 to 10 solar masses. Star-forming infall is assumed equally likely to stop at any moment, due to gas dispersal dominated by stellar feedback. For spherical infall, the typical initial condensation must have a steep density gradient, as in low-mass cores, surrounded by a shallower gradient, as in the clumps around cores. These properties match observed column densities in cluster-forming regions when the mean infall stopping time is 0.05 Myr and the accretion efficiency is 0.5. The infall duration increases with final protostar mass, from 0.01 to 0.3 Myr, and the mass accretion rate increases from 3 to 300 x 10^(-6) solar masses/yr. The typical spherical accretion luminosity is ~5 solar luminosities, reducing the luminosity problem to a factor ~3. The initial condensation density gradient changes from steep to shallow at radius 0.04 pc, enclosing 0.9 solar masses, with mean column density 2 x 10^(22) cm^(-2), and with effective central temperature 16 K. These initial conditions are denser and warmer than those for isolated star formation.

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 use a new contour-based map analysis technique to measure the mass and size of molecular cloud fragments continuously over a wide range of spatial scales (0.05 < r / pc < 10), i.e., from the scale of dense cores to those of entire clouds. The present paper presents the method via a detailed exploration of the Perseus Molecular Cloud. Dust extinction and emission data are combined to yield reliable scale-dependent measurements of mass. This scale-independent analysis approach is useful for several reasons. First, it provides a more comprehensive characterization of a map (i.e., not biased towards a particular spatial scale). Such a lack of bias is extremely useful for the joint analysis of many data sets taken with different spatial resolution. This includes comparisons between different cloud complexes. Second, the multi-scale mass-size data constitutes a unique resource to derive slopes of mass-size laws (via power-law fits). Such slopes provide singular constraints on large-scale density gradients in clouds.

The distribution of protostar masses is studied for core-environment systems whose duration of infall follows a waiting-time distribution. Each core-environment system has a continuous density profile with no barrier to mass flow. The core is an isothermal sphere and the environment is a filament, a layer, or a uniform medium. The infall is terminated by gas dispersal due to outflows and turbulence. The distribution of infall durations is a declining exponential, the simplest waiting-time distribution. The resulting distribution of protostar masses closely resembles the initial mass function, provided the environment density is sufficiently high, and the distribution of initial core masses is sufficiently narrow. The high-mass tail of the mass function increases strongly with environment density and weakly with environment dimension. Isolated regions of low environment density form protostars of low mass from within the parent core. In contrast, clustered regions of high environment density form protostars of low mass from core gas, and protostars of high mass from core and environment gas.

Fully sampled degree-scale maps of the 13CO 2-1 and CO 4-3 transitions toward three members of the Lupus Molecular Cloud Complex - Lupus I, III, and IV - trace the column density and temperature of the molecular gas. Comparison with IR extinction maps from the c2d project requires most of the gas to have a temperature of 8-10 K. Estimates of the cloud mass from 13CO emission are roughly consistent with most previous estimates, while the line widths are higher, around 2 km/s. CO 4-3 emission is found throughout Lupus I, indicating widespread dense gas, and toward Lupus III and IV. Enhanced line widths at the NW end and along the edge of the B228 ridge in Lupus I, and a coherent velocity gradient across the ridge, are consistent with interaction between the molecular cloud and an expanding HI shell from the Upper-Scorpius subgroup of the Sco-Cen OB Association. Lupus III is dominated by the effects of two HAe/Be stars, and shows no sign of external influence. Slightly warmer gas around the core of Lupus IV and a low line width suggest heating by the Upper-Centaurus-Lupus subgroup of Sco-Cen, without the effects of an HI shell.

The key question about early protostellar evolution is how matter is accreted from the large-scale molecular cloud, through the circumstellar disk onto the central star. A sample of 20 Class 0 and I protostars has been observed in continuum at (sub)millimeter wavelengths at high angular resolution with the Submillimeter Array. Using detailed dust radiative transfer models, we have developed a framework for disentangling the continuum emission from the envelopes and disks, and from that estimated their masses. For the Class I sources in the sample, HCO+ 3-2 line emission has furthermore been observed with the Submillimeter Array. Four of these sources show signs of Keplerian rotation, constraining the masses of the central stars. Both Class 0 and I protostars are surrounded by disks with typical masses of about 0.05 M_sun. No evidence is found for a correlation between the disk mass and evolutionary stage of the young stellar objects. This contrasts the envelope mass, which decreases sharply from 1 M_sun in the Class 0 stage to <0.1 M_sun in the Class I stage. Typically, the disks have masses that are 1-10% of the corresponding envelope masses in the Class 0 stage and 20-60% in the Class I stage. For the Class I sources for which Keplerian rotation is seen, the central stars contain 70-98% of the total mass in the star-disk-envelope system, confirming that these objects are late in their evolution through the embedded protostellar stages. Theoretical models tend to overestimate the disk masses relative to the stellar masses in the late Class I stage. The results argue in favor of a picture in which circumstellar disks are formed early during the protostellar evolution and rapidly process material accreted from the larger scale envelope onto the central star.

The nearest young stellar groups are associated with "hubs" of column density exceeding 10^22 cm^-2, according to recent observations. These hubs radiate multiple "filaments" of parsec length, having lower column density and fewer stars. Systems with many filaments tend to have parallel filaments with similar spacing. Such "hub-filament structure" is associated with all of the nine young stellar groups within 300 pc, forming low-mass stars. Similar properties are seen in infrared dark clouds forming more massive stars. In a new model, an initial clump in a uniform medium is compressed into a self-gravitating, modulated layer. The outer layer resembles the modulated equilibrium of Schmid-Burgk (1967) with nearly parallel filaments. The filaments converge onto the compressed clump, which collapses to form stars with high efficiency. The initial medium and condensations have densities similar to those in nearby star-forming clouds and clumps. The predicted structures resemble observed hub-filament systems in their size, shape, and column density, and in the appearance of their filaments. These results suggest that hub-filament structure associated with young stellar groups may arise from compression of clumpy gas in molecular clouds.

We present a uniform mid-infrared imaging and photometric survey of 36 young, nearby, star-forming clusters and groups using {\it Spitzer} IRAC and MIPS. We have confidently identified and classified 2548 young stellar objects using recently established mid-infrared color-based methods. We have devised and applied a new algorithm for the isolation of local surface density enhancements from point source distributions, enabling us to extract the overdense cores of the observed star forming regions for further analysis. We have compiled several basic structural measurements of these cluster cores from the data, such as mean surface densities of sources, cluster core radii, and aspect ratios, in order to characterize the ranges for these quantities. We find that a typical cluster core is 0.39 pc in radius, has 26 members with infrared excess in a ratio of Class II to Class I sources of 3.7, is embedded in a $A_K$=0.8 mag cloud clump, and has a surface density of 60 pc$^{-2}$. We examine the nearest neighbor distances among the YSOs in several ways, demonstrating similarity in the spacings between Class II and Class I sources but large member clusters appear more dense than smaller clusters. We demonstrate that near-uniform source spacings in cluster cores are common, suggesting that simple Jeans fragmentation of parsec-scale cloud clumps may be the dominant process governing star formation in nearby clusters and groups. Finally, we compare our results to other similar surveys in the literature and discuss potential biases in the data to guide further interpretation.

We present a comprehensive analysis of structure in the young, embedded cluster, NGC 1333 using members identified with Spitzer and 2MASS photometry based on their IR-excess emission. In total, 137 members are identified in this way, composed of 39 protostars and 98 more evolved pre-main sequence stars with disks. Of the latter class, four are transition/debris disk candidates. The fraction of exposed pre-main sequence stars with disks is 83% +/- 11%, showing that there is a measurable diskless pre-main sequence population. The sources in each of the Class I and Class II evolutionary states are shown to have very different spatial distributions relative to the distribution of the dense gas in their natal cloud. However, the distribution of nearest neighbor spacings among these two groups of sources are found to be quite similar, with a strong peak at spacings of 0.045 pc. Radial and azimuthal density profiles and surface density maps computed from the identified YSOs show that NGC 1333 is elongated and not strongly centrally concentrated, confirming previous claims in the literature. We interpret these new results as signs of a low velocity dispersion, extremely young cluster that is not in virial equilibrium.

We present spectral observations of 130 young stellar objects (YSOs) in the Serpens Cloud Core and NGC 1333 embedded clusters. The observations consist of near-IR spectra in the H and K-bands, from SpeX on the IRTF and far-red spectra (6000 - 9000 A) from Hectospec on the MMT. These YSOs were identified in previous Spitzer and Chandra observations, and the evolutionary classes of the YSOs were determined from the Spitzer mid-IR photometry. With these spectra, we search for corroborating evidence for the pre-main sequence nature of the objects, study the properties of the detected emission lines as a function of evolutionary class, and obtain spectral types for the observed YSOs. By comparing the positions of the YSOs in the HR diagrams with the pre-main sequence tracks of Baraffe (1998), we determine ages of the embedded sources and study the relative ages of the YSOs with and without optically thick circumstellar disks. The apparent isochronal ages of the YSOs in both clusters range from less than 1 Myr to 10 Myr, with most objects below 3 Myr. The observed distribution of ages for the Class II and Class III objects are statistically indistinguishable. We examine the spatial distribution and extinction of the YSOs as a function of their isochronal ages. We find the sources < 3 Myr to be concentrated in the molecular cloud gas while the older sources are spatially dispersed and are not deeply embedded. Nonetheless, the sources with isochronal ages > 3 Myr show all the characteristics of young stellar objects in their spectra, their IR spectral energy distributions, and their X-ray emission.

We present combined BIMA interferometer and IRAM 30 m Telescope data of N2H+ 1-0 line emission across the nearby dense, star forming core Ophiuchus A (Oph A) at high linear resolution (e.g., ~1000 AU). Six maxima of integrated line intensity are detected which we designate Oph A-N1 through N6. The N4 and N5 maxima are coincident with the starless continuum objects SM1 and SM2 respectively but the other maxima are not coincident with previously-identified objects. In contrast, relatively little N2H+ 1-0 emission is coincident with the starless object SM2 and the Class 0 protostar VLA 1623. The FWHM of the N2H+ 1-0 line, Delta V, varies by a factor of ~5 across Oph A. Values of Delta V < 0.3 km/s are found in 14 locations in Oph A, but only that associated with N6 is both well-defined spatially and larger than the beam size. Centroid velocities of the line, V_LSR, vary relatively little, having an rms of only \~0.17 km/s. Small-scale V_LSR gradients of <0.5 km/s over ~0.01 pc are found near SM1, SM1N, and SM2, but not N6. The low N2H+ abundances of SM2 or VLA 1623 relative to SM1, SM1N, or N6 may reflect relatively greater amounts of N2 adsorption onto dust grains in their colder and probably denser interiors. The low Delta V of N6, i.e., 0.193 km/s FWHM, is only marginally larger than the FWHM expected from thermal motions alone, suggesting turbulent motions in the Oph A core have been reduced dramatically at this location. The non-detection of N6 in previous thermal continuum maps suggests that interesting sites possibly related to star formation may be overlooked in such data.

We present combined interferometer and single dish telescope data of NH3 (J,K) = (1,1) and (2,2) emission towards the clustered star forming Ophiuchus B, C and F Cores at high spatial resolution (~1200 AU) using the Australia Telescope Compact Array, the Very Large Array, and the Green Bank Telescope. While the large scale features of the NH3 (1,1) integrated intensity appear similar to 850 micron continuum emission maps of the Cores, on 15" (1800 AU) scales we find significant discrepancies between the dense gas tracers in Oph B, but good correspondence in Oph C and F. Using the Clumpfind structure identifying algorithm, we identify 15 NH3 clumps in Oph B, and 3 each in Oph C and F. Only five of the Oph B NH3 clumps are coincident within 30" (3600 AU) of a submillimeter clump. We find v_LSR varies little across any of the Cores, and additionally varies by only ~1.5 km/s between them. The observed NH3 line widths within the Oph B and F Cores are generally large and often mildly supersonic, while Oph C is characterized by narrow line widths which decrease to nearly thermal values. We find several regions of localized narrow line emission (\Delta v < 0.4 km/s), some of which are associated with NH3 clumps. We derive the kinetic temperatures of the gas, and find they are remarkably constant across Oph B and F, with a warmer mean value (T_K = 15 K) than typically found in isolated regions and consistent with previous results in clustered regions. Oph C, however, has a mean T_K = 12 K, decreasing to a minimum T_K = 9.4 K towards the submillimeter continuum peak, similar to previous studies of isolated starless cores. There is no significant difference in temperature towards protostars embedded in the Cores. [Abridged]

We present the chemistry, temperature, and dynamical state of a sample of 193 dense cores or core candidates in the Perseus Molecular cloud and compare the properties of cores associated with young stars and clusters with those which are not. The combination of our NH3 and CCS observations with previous millimeter, sub-millimeter, and Spitzer data available for this cloud enable us both to determine core properties precisely and to accurately classify cores as starless or protostellar. The properties of cores in different cluster environments and before-and-after star formation provide important constraints on simulations of star-formation, particularly under the paradigm that the essence of star formation is set by the turbulent formation of prestellar cores. We separate the influence of stellar content from that of cluster environment and find that cores within clusters have (1) higher kinetic temperatures and (2) lower fractional abundances of CCS and NH3. Cores associated with protostars have (1) slightly higher kinetic temperatures (2) higher NH3 excitation temperatures), (3) are at higher column density, have (4) slightly more non-thermal/turbulent NH3 linewidths, have (5) higher masses and have (6) lower fractional abundance of CCS. We find that neither cluster environment nor protostellar content makes a significant difference to the dynamical state of cores as estimated by the virial parameter -- most cores in each category are gravitationally bound. Overall, cluster environment and protostellar content have a smaller influence on the properties of the cores than is typically assumed, and the variation within categories is larger than the differences between categories.

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

We present Spitzer infrared observations of the starless core L429. The IR images of this core show an absorption feature, caused by the dense core material, at wavelengths <= 70 micron. The core has a steep density profile, and reaches A_V > 35 mag near the center. We show that L429 is either collapsing or in a near-collapse state.

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

The mass of a protostar is calculated from the infall and dispersal of an isothermal sphere in a uniform background. For high contrast between peak and background densities and for short dispersal time t_d, the accretion is "self-limiting": gas beyond the core is dispersed before it accretes, and the protostar mass approaches a time-independent value of low mass. For lower density contrast and longer dispersal time, the accretion "runs away": gas accretes from beyond the core, and the protostar mass approaches massive star values. The final protostar mass is approximately the initial gas mass whose free-fall time equals t_d. This mass matches the peak of the IMF for gas temperature 10 K, peak and background densities 10^6 and 10^3 cm^-3, and for t_d comparable to the core free-fall time t_core. The accretion luminosity exceeds 1 L-Sun for 0.1 Myr, as in the "Class 0" phase. For t_d/t_core=0.4-0.8 and temperature 7-50 K, self-limiting protostar masses are 0.08-5 M-Sun. These protostar and core masses have ratio 0.4 +- 0.2, as expected if the core mass distribution and the IMF have the same shape.

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

We present an unbiased census of starless cores in Perseus, Serpens, and Ophiuchus, assembled by comparing large-scale Bolocam 1.1 mm continuum emission maps with Spitzer c2d surveys. We use the c2d catalogs to separate 108 starless from 92 protostellar cores in the 1.1 mm core samples from Enoch et al. (2006), Young et al. (2006), and Enoch et al. (2007). A comparison of these populations reveals the initial conditions of the starless cores. Starless cores in Perseus have similar masses but larger sizes and lower densities on average than protostellar cores, with sizes that suggest density profiles substantially flatter than r^-2. By contrast, starless cores in Serpens are compact and have lower masses than protostellar cores; future star formation will likely result in lower mass objects than the currently forming protostars. Comparison to dynamical masses estimated from the NH_3 survey of Perseus cores by Rosolowsky et al. (2007) suggests that most of the starless cores are likely to be gravitationally bound, and thus prestellar. The combined prestellar core mass distribution includes 108 cores and has a slope of -2.3+/-0.4 for M>0.8 Msun. This slope is consistent with recent measurements of the stellar initial mass function, providing further evidence that stellar masses are directly linked to the core formation process. We place a lower limit on the core-to-star efficiency of 25%. There are approximately equal numbers of prestellar and protostellar cores in each cloud, thus the dense prestellar core lifetime must be similar to the lifetime of embedded protostars, or 4.5x10^5 years, with a total uncertainty of a factor of two. Such a short lifetime suggests a dynamic, rather than quasi-static, core evolution scenario, at least at the relatively high mean densities (n>2x10^4 cm^-3) to which we are sensitive.

We present a census of the population of deeply embedded young stellar objects (YSOs) in the Ophiuchus molecular cloud complex based on a combination of Spitzer Space Telescope mid-infrared data from the "Cores to Disks" (c2d) legacy team and JCMT/SCUBA submillimeter maps from the COMPLETE team. We have applied a method developed for identifying embedded protostars in Perseus to these datasets and in this way construct a relatively unbiased sample of 27 candidate embedded protostars with envelopes more massive than our sensitivity limit (about 0.1 M_sun). Embedded YSOs are found in 35% of the SCUBA cores - less than in Perseus (58%). On the other hand the mid-infrared sources in Ophiuchus have less red mid-infrared colors, possibly indicating that they are less embedded. We apply a nearest neighbor surface density algorithm to define the substructure in each of the clouds and calculate characteristic numbers for each subregion - including masses, star formation efficiencies, fraction of embedded sources etc. Generally the main clusters in Ophiuchus and Perseus (L1688, NGC1333 and IC348) are found to have higher star formation efficiencies than small groups such as B1, L1455 and L1448, which on the other hand are completely dominated by deeply embedded protostars. We discuss possible explanations for the differences between the regions in Perseus and Ophiuchus, such as different evolutionary timescales for the YSOs or differences, e.g., in the accretion in the two clouds.

We study the collapse of protoclusters within a giant molecular cloud (GMC) to determine the conditions under which collapse is significantly disrupted. Motivated by observations of star forming regions which exhibit flattened cloud structures, this study considers collapsing protoclusters with disk geometries. The collapse of a 10^3 Msun protocluster initially a distance of 2-10 pc from a 10^3 - 10^6 Msun point mass is numerically calculated. Simulations with zero initial relative velocity between the two are completed as well as simulations with relative velocities consistent with those observed in GMCs. The results allow us to define the conditions under which it is safe to assume protocluster collapse proceeds as if in isolation. For instance, we find the collapse of a 10^3 Msun protocluster will be significantly disrupted if it is within 2-4 pc of a 10^4 Msun point mass. Thus, the collapse of a 10^3 Msun protocluster can be considered to proceed as if in isolation if it is more than ~ 4 pc away from a 10^4 Msun compact object. In addition, in no portion of the sampled parameter space does the gravitational interaction between the protocluster disk and the massive particle significantly disperse the disk into the background GMC. We discuss the distribution of clusters of young stellar objects within the Perseus and Mon R2 star forming regions, which are consistent with the results of our simulations and the limitations of our results in gas dominated regions such as the Orion cloud.

We present a census of circumstellar disks in the Chamaeleon I star-forming region. Using the Infrared Array Camera and the Multiband Imaging Photometer onboard the Spitzer Space Telescope, we have obtained images of Chamaeleon I at 3.6, 4.5, 5.8, 8.0, and 24 um. To search for new disk-bearing members of the cluster, we have performed spectroscopy on objects that have red colors in these data. Through this work, we have discovered four new members of Chamaeleon I with spectral types of M4, M6, M7.5, and L0. The first three objects are highly embedded (A_J~5) and reside near known protostars, indicating that they may be among the youngest low-mass sources in the cluster (<1 Myr). The L0 source is the coolest known member of Chamaeleon I. Its luminosity implies a mass of 0.004-0.01 M_sun, making it the least massive brown dwarf for which a circumstellar disk has been reliably detected. To characterize the disk population in Chamaeleon I, we have classified the infrared spectral energy distributions of the 203 known members that are encompassed by the Spitzer images. Through these classifications, we find that the disk fraction in Chamaeleon I is roughly constant at ~50% from 0.01 to 0.3 M_sun. These data are similar to the disk fraction of IC 348, which is a denser cluster at the same age as Chamaeleon I. However, the disk fraction at M>1 M_sun is significantly higher in Chamaeleon I than in IC 348 (65% vs. 20%), indicating longer disk lifetimes in Chamaeleon I for this mass range. Thus, low-density star-forming regions like Chamaeleon I may offer more time for planet formation around solar-type stars than denser clusters.

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

We present {\it Spitzer} images of the relatively sparse, low luminosity young cluster L988e, as well as complementary near-infrared (NIR) and submillimeter images of the region. The cluster is asymmetric, with the western region of the cluster embedded within the molecular cloud, and the slightly less dense eastern region to the east of, and on the edge of, the molecular cloud. With these data, as well as with extant H$\alpha$ data of stars primarily found in the eastern region of the cluster, and a molecular $^{13}$CO gas emission map of the entire region, we investigate the distribution of forming young stars with respect to the cloud material, concentrating particularly on the differences and similarities between the exposed and embedded regions of the cluster. We also compare star formation in this region to that in denser, more luminous and more massive clusters already investigated in our comprehensive multi-wavelength study of young clusters within 1 kpc of the Sun.

We present ammonia observations of 193 dense cores and core candidates in the Perseus molecular cloud made using the Robert F. Byrd Green Bank Telescope. We simultaneously observed the NH3(1,1), NH3(2,2), CCS (2_1 -> 1_0) and CC34S (2_1 -> 1_0) transitions near 23 GHz for each of the targets with a spectral resolution of dv ~ 0.024 km/s. We find ammonia emission associated with nearly all of the (sub)millimeter sources as well as at several positions with no associated continuum emission. For each detection, we have measured physical properties by fitting a simple model to every spectral line simultaneously. Where appropriate, we have refined the model by accounting for low optical depths, multiple components along the line of sight and imperfect coupling to the GBT beam. For the cores in Perseus, we find a typical kinetic temperature of T=11 K, a typical column density of N(NH3)~ 10^14.5 /cm^2 and velocity dispersions ranging from sigma_v = 0.07 km/s to 0.7 km/s. However, many cores with velocity dispersions > 0.2 km/s show evidence for multiple velocity components along the line of sight.

We present maps of 14.4 deg^2 of the Ophiuchus dark clouds observed by the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS). These high quality maps depict both numerous point sources as well as extended dust emission within the star-forming and non-star-forming portions of these clouds. Using PSF-fitting photometry, we detect 5779 sources at 24 um and 81 sources at 70 um at the 10 sigma level of significance. Three hundred twenty-three candidate young stellar objects (YSOs) were identified according to their positions on the MIPS/2MASS K versus K$-$[24] color-magnitude diagrams as compared to 24 um detections in the SWIRE extragalactic survey. We find that more than half of the YSO candidates, and almost all the ones with protostellar Class I spectral energy distributions, are confined to the known cluster and aggregates.

We present Spitzer and Chandra observations of the nearby (~260 pc) embedded stellar cluster in the Serpens Cloud Core. We observed, using Spitzer's IRAC and MIPS instruments, in six wavelength bands from 3 to 70 ${\mu}m$, to detect thermal emission from circumstellar disks and protostellar envelopes, and to classify stars using color-color diagrams and spectral energy distributions (SEDs). These data are combined with Chandra observations to examine the effects of circumstellar disks on stellar X-ray properties. Young diskless stars were also identified from their increased X-ray emission. We have identified 138 YSOs in Serpens: 22 class 0/I, 16 flat spectrum, 62 class II, 17 transition disk, and 21 class III stars; 60 of which exhibit X-ray emission. Our primary results are the following: 1.) ten protostars detected previously in the sub-millimeter are detected at lambda < 24 microns, seven at lambda < 8 microns, 2.) the protostars are more closely grouped than more evolved YSOs (median separation : ~0.024 pc, and 3.) the luminosity and temperature of the X-ray emitting plasma around these YSOs does not show any significant dependence on evolutionary class. We combine the infrared derived values of AK and X-ray values of NH for 8 class III objects and find that the column density of hydrogen gas per mag of extinctions is less than half the standard interstellar value, for AK > 1. This may be the result of grain growth through coagulation and/or the accretion of volatiles in the Serpens cloud core.

We present maps of 7.78 square degrees of the Lupus molecular cloud complex at 24, 70, and $160\:\mu$m. They were made with the Spitzer Space Telescope's Multiband Imaging Photometer for Spitzer (MIPS) instrument as part of the Spitzer Legacy Program, ``From Molecular Cores to Planet-Forming Disks'' (c2d). The maps cover three separate regions in Lupus, denoted I, III, and IV. We discuss the c2d pipeline and how our data processing differs from it. We compare source counts in the three regions with two other data sets and predicted star counts from the Wainscoat model. This comparison shows the contribution from background galaxies in Lupus I. We also create two color magnitude diagrams using the 2MASS and MIPS data. From these results, we can identify background galaxies and distinguish them from probable young stellar objects. The sources in our catalogs are classified based on their spectral energy distribution (SED) from 2MASS and Spitzer wavelengths to create a sample of young stellar object candidates. From 2MASS data, we create extinction maps for each region and note a strong corresponence between the extinction and the $160\:\mu$m emission. The masses we derived in each Lupus cloud from our extinction maps are compared to masses estimated from $^{13}$CO and C$^{18}$O and found to be similar to our extinction masses in some regions, but significantly different in others. Finally, based on our color-magnitude diagrams, we selected 12 of our reddest candidate young stellar objects for individual discussion. Five of the 12 appear to be newly-discovered YSOs.

One of the central goals of the Spitzer Legacy Project ``From Molecular Cores to Planet-forming Disks'' (c2d) is to determine the frequency of remnant circumstellar disks around weak-line T Tauri stars (wTTs) and to study the properties and evolutionary status of these disks. Here we present a census of disks for a sample of over 230 spectroscopically identified wTTs located in the c2d IRAC (3.6, 4.5, 4.8, and 8.0 um) and MIPS (24 um) maps of the Ophiuchus, Lupus, and Perseus Molecular Clouds. We find that ~20% of the wTTs in a magnitude limited subsample have noticeable IR-excesses at IRAC wavelengths indicating the presence of a circumstellar disk. The disk frequencies we find in these 3 regions are ~3-6 times larger than that recently found for a sample of 83 relatively isolated wTTs located, for the most part, outside the highest extinction regions covered by the c2d IRAC and MIPS maps. The disk fractions we find are more consistent with those obtained in recent Spitzer studies of wTTs in young clusters such as IC 348 and Tr 37. From their location in the H-R diagram, we find that, in our sample, the wTTs with excesses are among the younger part of the age distribution. Still, up to ~50% of the apparently youngest stars in the sample show no evidence of IR excess, suggesting that the circumstellar disks of a sizable fraction of pre-main-sequence stars dissipate in a timescale of ~1 Myr. We also find that none of the stars in our sample apparently older than ~10 Myrs have detectable circumstellar disks at wavelengths < 24 um. Also, we find that the wTTs disks in our sample exhibit a wide range of properties (SED morphology, inner radius, L_DISK/L*, etc) which bridge the gaps observed between the cTTs and the debris disk regimes.

We have carried out a survey toward the central regions of 85 starless cores in HCN J = 1-0 to study inward motions in the cores. Sixty-four cores were detected with HCN lines. The infall asymmetry in the HCN spectra is found to be more prevalent, and more prominent than in any other previously used infall tracers such as CS J = 2-1, DCO+ J = 2-1, and N2H+ J = 1-0. We found close relation between the intensities of the HCN and N2H+ lines. This implies that the HCN is not much depleted in the central regions of the cores. In some cores, the HCN spectra show different sign of asymmetry from other molecular lines. A few cores show various signs of asymmetry in individual HCN hyperfine lines. The distribution of the velocity shift dV of the HCN profiles with respect to the systemic velocity of the optically thin tracer is found to be more shifted toward bluer side than those of other infall tracers, indicating that the HCN traces inward motions more frequently. The dV distribution of each HCN hyperfine line for all sources is similar. Moreover the dV values obtained from different HCN hyperfine lines for each source are nearly similar. These may mean that most of starless cores are in similar kinematic states across the layers of the cores. We identify 17 infall candidates using all available indicators such as the velocity shift dV and the blue to red peak intensity ratio of double peaked profiles for HCN J = 1-0, CS J = 2-1, J = 3-2, DCO+ J = 2-1, and N2H+ J = 1-0. Four of them, L63, L492, L694-2, and L1197 are found to show higher blue to red ratio in the HCN hyperfine line along the lower opacity, suggesting that infall speed becomes higher toward the center.

We discuss the results from the combined IRAC and MIPS c2d Spitzer Legacy observations of the Serpens star-forming region. In particular we present a set of criteria for isolating bona fide young stellar objects, YSO's, from the extensive background contamination by extra-galactic objects. We then discuss the properties of the resulting high confidence set of YSO's. We find 235 such objects in the 0.85 deg^2 field that was covered with both IRAC and MIPS. An additional set of 51 lower confidence YSO's outside this area is identified from the MIPS data combined with 2MASS photometry. We describe two sets of results, color-color diagrams to compare our observed source properties with those of theoretical models for star/disk/envelope systems and our own modeling of the subset of our objects that appear to be star+disks. These objects exhibit a very wide range of disk properties, from many that can be fit with actively accreting disks to some with both passive disks and even possibly debris disks. We find that the luminosity function of YSO's in Serpens extends down to at least a few x .001 Lsun or lower for an assumed distance of 260 pc. The lower limit may be set by our inability to distinguish YSO's from extra-galactic sources more than by the lack of YSO's at very low luminosities. A spatial clustering analysis shows that the nominally less-evolved YSO's are more highly clustered than the later stages and that the background extra-galactic population can be fit by the same two-point correlation function as seen in other extra-galactic studies. We also present a table of matches between several previous infrared and X-ray studies of the Serpens YSO population and our Spitzer data set.

We present maps of 1.5 square degrees of the Serpens dark cloud at 24, 70, and 160\micron observed with the Spitzer Space Telescope MIPS Camera. More than 2400 compact sources have been extracted at 24um, nearly 100 at 70um, and 4 at 160um. We estimate completeness limits for our 24um survey from Monte Carlo tests with artificial sources inserted into the Spitzer maps. We compare source counts, colors, and magnitudes in the Serpens cloud to two reference data sets, a 0.50 deg^2 set on a low-extinction region near the dark cloud, and a 5.3 deg^2 subset of the SWIRE ELAIS N1 data that was processed through our pipeline. These results show that there is an easily identifiable population of young stellar object candidates in the Serpens Cloud that is not present in either of the reference data sets. We also show a comparison of visual extinction and cool dust emission illustrating a close correlation between the two, and find that the most embedded YSO candidates are located in the areas of highest visual extinction.

Forty molecular cloud cores in the southern hemisphere from the initial Spitzer Space Telescope Cores-to-Disks (c2d) Legacy program source list have been surveyed in 13CO(2-1), 12CO(4-3), and 12CO(7-6) with the Antarctic Submillimeter Telescope and Remote Observatory (AST/RO). The cores, ten of which contain embedded sources, are located mostly in the Vela, Ophiuchus, Lupus, Chamaeleon, Musca, and Scorpius complexes. 12CO(7-6) emission was undetected in all 40 clouds. We present data of 40 sources in 13CO(2-1) and 12CO(4-3), significant upper limits of 12CO(7-6), as well as a statistical analysis of the observed properties of the clouds. We find the typical 13CO(2-1) linewidth to be 2.0 km/s for cores with embedded stars, and 1.8 km/s for all others. The typical 12CO(4-3) linewidth is 2.6 to 3.7 km/s for cores with known embedded sources, and 1.6 to 2.3 km/s for all others. The average 13CO column density derived from the line intensities was found to be 1.9 x 10^15 cm^(-2) for cores with embedded stars, and 1.5 x 10^15 cm^(-2) for all others. The average kinetic temperature in the molecular cores, determined through a Large Velocity Gradient analysis of a set of nine cores, has an average lower limit of 16 K and an average upper limit of 26 K. The average molecular hydrogen density has an average lower limit of 10^2.9 cm^(-3) and an average upper limit of 10^3.3 cm^(-3) for all cores. For a different subset of nine cores, we have derived masses. They range from 4 to 255 M_sun. Overall, our c2d sample of southern molecular cores has a range of properties (linewidth, column density, size, mass, embedded stars) similar to those of past studies.

Using the Infrared Array Camera and the Multiband Imaging Photometer aboard the {\it Spitzer Space Telescope}, we have obtained images of the Lupus 3 star-forming cloud at 3.6, 4.5, 5.8, 8.0, and 24 \micron. We present photometry in these bands for the 41 previously known members that are within our images. In addition, we have identified 19 possible new members of the cloud based on red 3.6-8.0 \micron colors that are indicative of circumstellar disks. We have performed optical spectroscopy on 6 of these candidates, all of which are confirmed as young low-mass members of Lupus 3. The spectral types of these new members range from M4.75 to M8, corresponding to masses of 0.2-0.03 $M_\odot$ for ages of $\sim1$ Myr according to theoretical evolutionary models. We also present optical spectroscopy of a candidate disk-bearing object in the vicinity of the Lupus 1 cloud, 2M 1541-3345, which Jayawardhana & Ivanov recently classified as a young brown dwarf ($M\sim0.03$ $M_\odot$) with a spectral type of M8. In contrast to their results, we measure an earlier spectral type of M5.75$\pm$0.25 for this object, indicating that it is probably a low-mass star ($M\sim0.1$ $M_\odot$). In fact, according to its gravity-sensitive absorption lines and its luminosity, 2M 1541-3345 is older than members of the Lupus clouds ($\tau\sim1$ Myr) and instead is probably a more evolved pre-main-sequence star that is not directly related to the current generation of star formation in Lupus.

We present observations of 10.6 square degrees of the Perseus molecular cloud at 24, 70, and 160 microns with the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS). The image mosaics show prominent, complex extended emission dominated by illuminating B stars on the East side of the cloud, and by cold filaments of 160 micron emission on the West side. Of 3950 point sources identified at 24 microns, 1141 have 2MASS counterparts. A quarter of these populate regions of the Ks vs. Ks-[24] diagram that are distinct from stellar photospheres and background galaxies, and thus are likely to be cloud members with infrared excess. Nearly half (46%) of these 24 micron excess sources are distributed outside the IC 348 and NGC 1333 clusters. NGC 1333 shows the highest fraction of stars with flat or rising spectral energy distributions (28%), while Class II SEDs are most common in IC 348. These results are consistent with previous relative age determinations for the two clusters. The intercluster region contains several tightly clumped (r~0.1 pc) young stellar aggregates whose members exhibit a wide variety of infrared spectral energy distributions characteristic of different circumstellar environments. One possible explanation is a significant age spread among the aggregate members, such that some have had time to evolve more than others. Alternatively, if the aggregate members all formed at roughly the same time, then remarkably rapid circumstellar evolution would be required to account for the association of Class I and Class III sources at ages <~1 Myr. We highlight important results for several other objects as well (full abstract in the paper).

This paper presents a large spectral line and continuum survey of 8 deeply embedded, low-mass protostellar cores using the Submillimeter Array. Each source was observed in high excitation lines of some of the most common molecular species, CO, HCO+, CS, SO, H2CO, CH3OH and SiO. Line emission from 11 species originating from warm and dense gas have been imaged at high angular resolution (1-3"; typically 200-600 AU) together with continuum emission at 230 GHz (1.3 mm) and 345 GHz (0.8 mm). Compact continuum emission is observed for all sources which likely originates in marginally optically thick circumstellar disks, with typical lower limits to their masses of 0.1 M_sun (1-10% of the masses of their envelopes) and having a dust opacity law with beta approximately 1. Prominent outflows are present in CO 2-1 observations in all sources: the most diffuse outflows are found in the sources with the lowest ratios of disk-to-envelope mass, and it is suggested that these sources are in a phase where accretion of matter from the envelope has almost finished and the remainder of the envelope material is being dispersed by the outflows. Other characteristic dynamical signatures are found with inverse P Cygni profiles indicative of infalling motions seen in the 13CO 2-1 lines toward NGC1333-IRAS4A and -IRAS4B. Outflow-induced shocks are present on all scales in the protostellar environments and are most clearly traced by the emission of CH3OH in NGC1333-IRAS4A and -IRAS4B. These observations suggest that the emission of CH3OH and H2CO from these proposed "hot corinos" are related to the shocks caused by the protostellar outflows. Only one source, NGC1333-IRAS2A, has evidence for hot, compact CH3OH emission coincident with the embedded protostar.

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

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

We observed the clustered star forming complex NGC1333 with the BIMA and FCRAO telescopes in the transitions HCO+(1-0) and N2H+(1-0) over an 11'x11' area with resolution ~10" (0.015pc). The N2H+ emission follows very closely the submillimeter dust continuum emission, while HCO+ emission appears more spatially extended and also traces outflows. We have identified 93 N2H+ cores using the CLUMPFIND algorithm, and we derive N2H+ core masses between 0.05 and 2.5M_sun, with uncertainties of a factor of a few, dominated by the adopted N2H+ abundance. From a comparison with virial masses, we argue that most of these N2H+ cores are likely to be bound, even at the lowest masses, suggesting that the cores do not trace transient structures, and implies the entire mass distribution consists of objects that can potentially form stars. We find that the mass distribution of N2H+ cores resembles the field star IMF, which suggests that the IMF is locked in at the pre-stellar stage of evolution. We find that the N2H+ cores associated with stars identified from Spitzer infrared images have a flat mass distribution. This might be because lower mass cores lose a larger fraction of their mass when forming a star. Even in this clustered environment, we find no evidence for ballistic motions of the cores relative to their lower density surroundings traced by isotopic CO emission, though this conclusion must remain tentative until the surroundings are observed at the same high resolution as the N2H+.

We present a census of the population of deeply embedded young stellar objects (YSOs) in the Perseus molecular cloud complex based on a combination of Spitzer Space Telescope mid-IR data from the c2d legacy team and JCMT/SCUBA submillimeter maps from the COMPLETE team. The mid-IR sources detected at 24 micron and having [3.6]-[4.5] > 1 are located close to the center of the SCUBA cores, typically within 15" of their peaks. The narrowness of the spatial distribution of mid-IR sources around the peaks of the SCUBA cores suggests that no significant dispersal of the newly formed YSOs has occurred. This argues against the suggestion that motions of protostars regulate the time scales over which significant (Bondi-Hoyle) accretion can occur. The most deeply embedded YSOs are found in regions with high extinction, AV > 5, similar to the extinction threshold observed for the SCUBA cores. All the SCUBA cores with high concentrations have embedded YSOs, but not all cores with low concentrations are starless. An unbiased sample of 49 deeply embedded YSOs is constructed. Embedded YSOs are found in 40 of the 72 SCUBA cores with only three cores harboring multiple embedded YSOs within 15". The equal number of SCUBA cores with and without embedded YSOs suggests that the time scale for the evolution through the dense prestellar stages, where the cores are recognized in the submillimeter maps and have central densities of 5e4-1e5 cm^{-3}, is similar to the time scale for the embedded protostellar stages. The current star formation efficiency of cores is estimated to be approximately 10-15%. In contrast, the star formation efficiency averaged over the cloud life time and compared to the total cloud mass is only a few percent, reflecting also the efficiency in assembling cloud material into the dense cores forming stars.

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

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

[Abridged] We present a molecular survey of the starless cores L1498 and L1517B. These cores have been selected for their relative isolation and close-to-round shape, and they have been observed in a number of lines of 13 molecular species (4 already presented in the first part of this series): CO, CS, N2H+, NH3, CH3OH, SO, C3H2, HC3N, C2S, HCN, H2CO, HCO+, and DCO+. Using a physical model of core structure and a Monte Carlo radiative transfer code, we determine for each core a self-consistent set abundances that fits simultaneously the observed radial profile of integrated intensity and the emergent spectrum towards the core center (for abundant species, optically thin isopologues are used). From this work, we find that L1498 and L1517B have similar abundance patterns, with most species suffering a significant drop toward the core center. This occurs for CO, CS, CH3OH, SO, C3H2, HC3N, C2S, HCN, H2CO, HCO+, and DCO+, which we fit with profiles having a sharp central hole. The size of this hole varies with molecule: DCO+, HCN, and HC3N have the smallest holes while SO, C2S and CO have the largest holes. Only N2H+ and NH3 are present in the gas phase at the core centers. From the different behavior of molecules, we select SO, C2S, and CH3OH as the most sensitive tracers of molecular depletion. Comparing our abundance determinations with the predictions from current chemical models we find order of magnitude discrepancies. Finally, we show how the ``contribution function'' can be used to study the formation of line profiles from the different regions of a core.

We present the first Spitzer Infrared Spectrograph (IRS; The IRS was a collaborative venture between Cornell University and Ball Aerospace Corporation funded by NASA through the Jet Propulsion Laboratory and the Ames Research Center.) observations of the disks around classical T Tauri stars: spectra in the 5.2-30 micron range of six stars. The spectra are dominated by emission features from amorphous silicate dust, and a continuous component from 5 to 8 microns that in most cases comprises an excess above the photosphere throughout our spectral range. There is considerable variation in the silicate feature/continuum ratio, which implies variations of inclination, disk flaring, and stellar mass accretion rate. In most of our stars, structure in the silicate feature suggests the presence of a crystalline component. In one, CoKu Tau/4, no excess above the photosphere appears at wavelengths shortward of the silicate features, similar to 10 Myr old TW Hya, Hen 3-600, and HR 4796A. This indicates the optically thick inner disk is largely absent. The silicate emission features with peaks at 9.7 and 18 microns indicate small dust grains are present. The extremely low 10-20 micron color temperature of the dust excess, 135 K, indicates these grains are located more than 10 AU from the star. These features are suggestive of gravitational influence by planets or close stellar companions and grain growth in the region within 10 AU of the star, somewhat surprising for a star this young (1 Myr).

We present a comprehensive study of a low-mass star-forming region,L1251B, at wavelengths from the near-infrared to the millimeter. L1251B, where only one protostar, IRAS 22376+7455, was known previously, is confirmed to be a small group of protostars based on observations with the Spitzer Space Telescope. The most luminous source of L1251B is located 5" north of the IRAS position. A near-infrared bipolar nebula, which is not associated with the brightest object and is located at the southeast corner of L1251B, has been detected in the IRAC bands. OVRO and SMA interferometric observations indicate that the brightest source and the bipolar nebula source in the IRAC bands are deeply embedded disk sources.Submillimeter continuum observations with single-dish telescopes and the SMA interferometric observations suggest two possible prestellar objects with very high column densities. Outside of the small group, many young stellar object candidates have been detected over a larger region of 12' x 12'. Extended emission to the east of L1251B has been detected at 850 micron; this "east core" may be a site for future star formation since no point source has been detected with IRAC or MIPS. This region is therefore a possible example of low-mass cluster formation, where a small group of pre- and protostellar objects (L1251B) is currently forming, alongside a large starless core (the east core).

We present observations of 3.86 sq. deg. of the Perseus molecular cloud complex with the Spitzer Space Telescope Infrared Array Camera (IRAC). The maps show strong extended emission arising from shocked H2 in outflows in the region and from polycyclic aromatic hydrocarbon features. More than 120,000 sources are extracted toward the cloud. Based on their IRAC colors and comparison to off-cloud and extragalactic fields, we identify 400 candidate young stellar objects. About two thirds of these are associated with the young clusters IC348 and NGC1333, while the last third is distributed over the remaining cloud. We classify the young stellar objects according to the traditional scheme based on the slope of their spectral energy distributions. Significant differences are found for the numbers of embedded Class I objects relative to more evolved Class II objects in IC348, NGC1333 and the remaining cloud with the embedded Class I and "flat spectrum" YSOs constituting 14%, 36% and 47% of the total number of YSOs identified in each of these regions. The high number of Class I objects in the extended cloud (61% of the Class I objects in the entire cloud) suggests that a significant fraction of the current star formation is occuring outside the two main clusters. Finally we discuss a number of outflows and identify their driving sources, including the known deeply embedded Class 0 sources outside the two major clusters. The Class 0 objects are found to be detected by Spitzer and have very red [3.6]-[4.5] colors but do not show similarly red [5.8]-[8.0] colors. The Class 0 objects are easily identifiable in color-color diagrams plotting these two colors but are problematic to extract automatically due to the extended emission from shocked gas or scattered light in cavities related to the associated outflows.

Using the Spitzer Space Telescope, we have observed 90 weak-line and classical T Tauri stars in the vicinity of the Ophiuchus, Lupus, Chamaeleon, and Taurus star-forming regions as part of the Cores to Disks (c2d) Spitzer Legacy project. In addition to the Spitzer data, we have obtained contemporaneous optical photometry to assist in constructing spectral energy distributions. These objects were specifically chosen as solar-type young stars with low levels of H alpha emission, strong X-ray emission, and lithium absorption i.e. weak-line T Tauri stars, most of which were undetected in the mid-to-far IR by the IRAS survey. Weak-line T Tauri stars are potentially extremely important objects in determining the timescale over which disk evolution may take place. Our objective is to determine whether these young stars are diskless or have remnant disks which are below the detection threshold of previous infrared missions. We find that only 5/83 weak-line T Tauri stars have detectable excess emission between 3.6 and 70 micron which would indicate the presence of dust from the inner few tenths of an AU out to the planet-forming regions a few tens of AU from the star. Of these sources, two have small excesses at 24 microns consistent with optically thin disks; the others have optically thick disks already detected by previous IR surveys. All of the seven classical T Tauri stars show excess emission at 24 and 70 micron, although their properties vary at the shorter wavelengths. Our initial results show that disks are rare among young stars selected for their weak H alpha emission.

We examine the properties of embedded clusters within 1 kiloparsec using new data from the Spitzer Space Telescope, as well as recent results from 2MASS and other ground-based near-infrared surveys. We use surveys of entire molecular clouds to understand the range and distribution of cluster membership, size and surface density. The Spitzer data demonstrate clearly that there is a continuum of star- forming environments, from relative isolation to dense clusters. The number of members of a cluster is correlated with the cluster radius, such that the average surface density of clusters having a few to a thousand members varies by a factor of only a few. The spatial distributions of Spitzer-identified young stellar objects frequently show elongation, low density halos, and sub-clustering. The spatial distributions of protostars resemble the distribution of dense molecular gas, suggesting that their morphologies result directly from the fragmentation of the natal gas. We also examine the effects of the cluster environments on star and planet formation. Although Far-UV and Extreme-UV radiation from massive stars can truncate disks in a few million years, fewer than half of the young stars in our sample (embedded clusters within 1 kpc) are found in regions of strong FUV and EUV fields. Typical volume densities and lifetimes of the observed clusters suggest that dynamical interactions are not an important mechanism for truncating disks on solar system size scales.

Low-mass dense cores represent the state of molecular gas associated with the earliest phases of low-mass star formation. Such cores are called "protostellar" or "starless," depending on whether they do or do not contain compact sources of luminosity. In this chapter, the first half of the review of low-mass dense cores, we describe the numerous inferences made about the nature of starless cores as a result of recent observations, since these reveal the initial conditions of star formation. We focus on the identification of isolated starless cores and their internal physical and chemical properties, including morphologies, densities, temperatures, kinematics, and molecular abundances. These objects display a wide range of properties since they are each at different points on evolutionary paths from ambient molecular cloud material to cold, contracting, and centrally concentrated configurations with significant molecular depletions and, in rare cases, enhancements.

This paper studies the dynamical evolution of young stellar clusters with $N$ = 100 - 1000 members. We use N-body simulations to explore how evolution depends on system size $N$ and the initial conditions. Motivated by recent observations of extremely young systems, this study compares subvirial and virial starting states. Multiple realizations of equivalent cases (100 simulations per case) are used to build up a robust statistical description of these systems, e.g., distributions of closest approaches, mass profiles, and distributions of radial locations. These results provide a framework from which to assess the effects of these clusters on star and planet formation. The distributions of radial positions are used in conjunction with distributions of FUV luminosities (also calculated here) to determine the radiation exposure of circumstellar disks. The distributions of closest approaches are used in conjunction with scattering cross sections (calculated here from $10^5$ scattering experiments) to determine the probability of solar system disruption. We also use the nearby cluster NGC 1333 as a test case. Our main conclusion is that clusters in this size range have only a modest effect on forming planetary systems: Interaction rates are low so that the typical solar system experiences a single encounter within 1000 AU. Radiation exposure is low, with median FUV flux $G_0$ = 900, so that photoevaporation of disks is only important beyond 30 AU. Given the low interaction rates and modest radiation levels, we suggest that solar system disruption is a rare event in these clusters.

We present near and mid-infrared photometry obtained with the Spitzer Space Telescope of 300 known members of the IC348 cluster. We merge this photometry with existing ground-based optical and near-infrared photometry in order to construct optical-infrared spectral energy distributions (SEDs) for all the cluster members and present a complete atlas of these SEDs. We employ these observations to both investigate the frequency and nature of the circumstellar disk population in the cluster. The observations are sufficiently sensitive to enable the first detailed measurement of the disk frequency for very low mass stars at the peak of the stellar IMF. Using measurements of infrared excess between 3.6 and 8 microns we find the total frequency of disk-bearing stars in the cluster to be 50 +/- 6%. However, only 30 +/- 4% of the member stars are surrounded by optically thick, primordial disks, while the remaining disk-bearing stars are surrounded by what appear to be optically thin,anemic disks. The disk fraction appears to be a function of spectral type and stellar mass. The disk longevity and thus conditions for planet formation appear to be most favorable for stars which are of comparable mass to the sun. The optically thick disks around later type (> M4) stars appear to be less flared than the disks around earlier type stars. This may indicate a greater degree of dust settling and a more advanced evolutionary state for the late M disk population. Finally we find that the presence of an optically thick dust disk is correlated with gaseous accretion as measured by the strength of Halpha emission. These results suggest that it is more likely for dust disks to persist in the absence of active gaseous accretion than for active accretion to persist in the absence of dusty disks.

Using the Infrared Array Camera (IRAC) aboard the Spitzer Space Telescope, we have obtained mid-infrared photometry for 25 and 18 low-mass members of the IC 348 and Chamaeleon I star-forming clusters, respectively (>M6, M<=0.08 Msun). We find that 42+/-13% and 50+/-17% of the two samples exhibit excess emission indicative of circumstellar disks. In comparison, the disk fractions for stellar members of these clusters are 33+/-4% and 45+/-7% (M0-M6, 0.7 Msun>=M>=0.1 Msun). The similarity of the disk fractions of stars and brown dwarfs is consistent with a common formation mechanism and indicates that the raw materials for planet formation are available around brown dwarfs as often as around stars.

Recently, the Spitzer Space Telescope discovered L1014-IRS, a mid-infrared source with protostellar colors, toward the heretofore "starless" core L1014. We present deep near-infrared observations that show a scattered light nebula extending from L1014-IRS. This nebula resembles those typically associated with protostars and young stellar objects, tracing envelope cavities presumably evacuated by an outflow. The northern lobe of the nebula has an opening angle of ~100 degrees, while the southern lobe is barely detected. Its morphology suggests that the bipolar cavity and inferred protostellar disk is not inclined more than 30 degrees from an edge-on orientation. The nebula extends at least 8" from the source at Ks, strongly suggesting that L1014-IRS is embedded within L1014 at a distance of 200 pc rather than in a more distant cloud associated with the Perseus arm at 2.6 kpc. In this case, the apparently low luminosity of L1014-IRS, 0.090 Lsun, is consistent with it having a substellar mass. However, if L1014-IRS is obscured by a circumstellar disk, its luminosity and inferred mass may be greater. Using near-infrared colors of background stars, we investigate characteristics of the L1014 molecular cloud core. We determine a mass of 3.6 Msun for regions of the core with Av > 2 magnitudes. A comparison of the radial extinction profile of L1014 with other cores suggests that L1014 may be among the most centrally condensed cores known, perhaps indicative of the earliest stages of brown dwarf or star formation processes.

We present models for the chemistry in gas moving towards the ionization front of an HII region. When it is far from the ionization front, the gas is highly depleted of elements more massive than helium. However, as it approaches the ionization front, ices are destroyed and species formed on the grain surfaces are injected into the gas phase. Photodissociation removes gas phase molecular species as the gas flows towards the ionization front. We identify models for which the OH column densities are comparable to those measured in observations undertaken to study the magnetic fields in star forming regions and give results for the column densities of other species that should be abundant if the observed OH arises through a combination of the liberation of H2O from surfaces and photodissociation. They include CH3OH, H2CO, and H2S. Observations of these other species may help establish the nature of the OH spatial distribution in the clouds, which is important for the interpretation of the magnetic field results.

We report ``infall asymmetry'' in the HCO$^+$ (1--0) and (3--2) lines toward NGC 1333, extended over $\sim 0.39 {\rm pc}^2$, a larger extent than has been reported be fore, for any star-forming region. The infall asymmetry extends over a major portion of the star-forming complex, and is not limited to a single protostar, or to a single dense core, or to a single spectral line. It seems likely that the infall asymmetry represents inward motions, and that these motions are physically associated with the complex. Both blue-asymmetric and red-asymmetric lines are seen, but in both the (3--2) and (1--0) lines of HCO$^+$ the vast majority of the asymmetric lines are blue, indicating inward motions. The (3--2) line, tracing denser gas, has the spectra with the strongest asymmetry and these spectra are associated with the protostars IRAS 4A and 4B, which most likely indicates a warm central source is affecting the line profiles. The (3--2) and (1--0) lines usually have the same sense of asymmetry in common positions, but their profiles differ significantly, and the (1--0) line appears to trace motions on much larger spatial scales than does the (3--2) line. Line profile models fit the spectra well, but do not strongly constrain their parameters. The mass accretion rate of the inward motions is of order 10$^{-4}$ M$_\odot$/yr, similar to the ratio of stellar mass to cluster age.

Using the Submillimeter Array we report the discovery of a compact low mass bipolar molecular outflow from L1014-IRS and confirm its association with the L1014 dense core at 200 pc. Consequently, L1014-IRS is the lowest luminosity (L \~0.09 Lsun) and perhaps the lowest mass source known to be driving a bipolar molecular outflow, which is one of the smallest known in size (~500 AU), mass (< 10^{-4} Msun), and energetics (e.g., force < 10^{-7} Msun km/s/yr).

We present interferometric observations of N2H+(1-0) in the starless, dense core L694-2 and compare them to previously published maps of L1544. Both cores are starless, centrally condensed, and show spectral signatures of rotation and collapse. We fit radially averaged spectra using a two-layer infall model and measure the variation of opacity and infall speed in each core. Both functions increase toward the center of each core but the radial gradients are shallower, and the central values lower, in L694-2. This general behavior is predicted in models of gravitational collapse with thermal plus magnetic support and the lower values in L694-2 may be due to its lower mass or a slightly earlier evolutionary state. In either case, it appears that both cores will form stars within a few 10^4 years.

We present initial results from a survey of the Orion A and B molecular clouds made with the InfraRed Array Camera (IRAC) onboard the Spitzer Space Telescope. This survey encompasses a total of 5.6 square degrees with the sensitivity to detect objects below the hydrogen burning limit at an age of 1 Myr. These observations cover a number of known star forming regions, from the massive star forming clusters in the Orion Nebula and NGC 2024, to small groups of low mass stars in the L1641. We combine the IRAC photometry with photometry from the 2MASS point source catalog and use the resulting seven band data to identify stars with infrared excesses due to dusty disks and envelopes. Using the presence of an infrared excess as an indicator of youth, we show the distribution of young stars and protostars in the two molecular clouds. We find that roughly half of the stars are found in dense clusters surrounding the two regions of recent massive star formation in the Orion clouds, NGC 2024 and the Orion Nebula.

We present an analysis of K-band stellar distributions for the young stellar clusters GGD 12-15, IRAS 20050+2720, and NGC 7129. We find that the two deeply embedded clusters, GGD 12-15 and IRAS 20050+2720, are not azimuthally symmetric and show a high degree of structure which traces filamentary structure observed in 850 micron emission maps. In contrast, the NGC 7129 cluster is circularly symmetric, less dense, and anti-correlated to 850 micron emission, suggesting recent gas expulsion and dynamical expansion have occured. We estimate stellar volume densities from nearest neighbor distances, and discuss the impact of these densities on the evolution of circumstellar disks and protostellar envelopes in these regions.

We present spectra of four Herbig AeBe stars obtained with the Infrared Spectrograph (IRS). on the Spitzer Space Telescope. All four of the sources show strong emission from polycyclic aromatic hydrocarbons (PAHs), with the 6.2 um emission feature shifted to 6.3 um and the strongest C-C skeletal-mode feature occuring at 7.9 um instead of at 7.7 um as is often seen. Remarkably, none of the four stars have silicate emission. The strength of the 7.9 um feature varies with respect to the 11.3 um feature among the sources, indicating that we have observed PAHs with a range of ionization fractions. The ionization fraction is higher for systems with hotter and brighter central stars. Two sources, HD 34282 and HD 169142, show emission features from aliphatic hydrocarbons at 6.85 and 7.25 um. The spectrum of HD 141569 shows a previously undetected emission feature at 12.4 um which may be related to the 12.7 um PAH feature. The spectrum of HD 135344, the coolest star in our sample, shows an unusual profile in the 7-9 um region, with the peak emission to the red of 8.0 um and no 8.6 um PAH feature.

The Submillimeter Array has opened a new window to study the innermost warm and dense regions of the envelopes and disks around deeply embedded protostars. This paper presents high-angular resolution (< 2'') submillimeter observations of the class 0 young stellar object NGC1333-IRAS2A. Dust continuum emission and lines of complex organic molecules such as CH3OCH3 and CH3OCHO, high excitation CH3OH transitions, deuterated methanol CH3OD as well as lines of CO, HCN, H13CN, SO and SO2 are detected on < 200 AU scales. The observations are interpreted using detailed radiative transfer models of the physical and chemical structure, consistent with both single-dish and interferometer data. The continuum emission is explained by an extended envelope and a compact but resolved component, presumably a circumstellar disk with a diameter of 200-300 AU and a mass of a few times 0.01-0.1 M_sun. If related to the rotation of the envelope, then the size of this disk suggests a centrifugal barrier of 200-300 AU, which implies that the temperature in the envelope does not increase above 100 K. Its large size also suggests that the build-up of disks proceeds rapidly throughout the early protostellar stages. The smaller (< 100 AU) disks found around other deeply embedded protostars may be a result of tidal truncation. The high-resolution observations of SO can be explained with a simple constant abundance, ~1E-9, constrained through single-dish observations, whereas those of H13CN and the organic species require high abundances, increased by one to two orders of magnitude, or an additional compact source of emission at small scales. The compact molecular emission could originate in a hot core region of the inner envelope, but a more likely reservoir is the circumstellar disk.

Motivated by recent observations that show starless molecular cloud cores exhibit subsonic inward velocities, we revisit the collapse problem for polytropic gaseous spheres. In particular, we provide a generalized treatment of protostellar collapse in the spherical limit and find semi-analytic (self-similar) solutions, corresponding numerical solutions, and purely analytic calculations of the mass infall rates (the three approaches are in good agreement). This study focuses on collapse solutions that exhibit nonzero inward velocities at large radii, as observed in molecular cloud cores, and extends previous work in four ways: (1) The initial conditions allow nonzero initial inward velocity. (2) The starting states can exceed the density of hydrostatic equilibrium so that the collapse itself can provide the observed inward motions. (3) We consider different equations of state, especially those that are softer than isothermal. (4) We consider dynamic equations of state that are different from the effective equation of state that produces the initial density distribution. This work determines the infall rates over a wide range of parameter space, as characterized by four variables: the initial inward velocity $\vin$, the overdensity $\overdense$ of the initial state, the index $\Gamma$ of the static equation of state, and the index $\gamma$ of the dynamic equation of state. For the range of parameter space applicable to observed cores, the resulting infall rate is about a factor of two larger than found in previous theoretical studies (those with hydrostatic initial conditions and $\vin = 0$).

We present maps of over 1.5 square degrees in Chamaeleon (Cha) II at 24, 70, and 160 micron observed with the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS) and a 1.2 square degree millimeter map from SIMBA on the Swedish-ESO Submillimetre Telescope (SEST). The c2d Spitzer Legacy Team's data reduction pipeline is described in detail. Over 1500 24 micron sources and 41 70 micron sources were detected by MIPS with fluxes greater than 10-sigma. More than 40 potential YSOs are identified with a MIPS and 2MASS color-color diagram and by their spectral indices, including two previously unknown sources with 24 micron excesses. Our new SIMBA millimeter map of Cha II shows that only a small fraction of the gas is in compact structures with high column densities. The extended emission seen by MIPS is compared with previous CO observations. Some selected interesting sources, including two detected at 1 mm, associated with Cha II are discussed in detail and their SEDs presented. The classification of these sources using MIPS data is found to be consistent with previous studies.

Spitzer Space Telescope observations of a point-like source, L1014-IRS, close to the dust peak of the low-mass dense core L1014 have questioned its starless nature. The presence of an object with colors expected for an embedded protostar makes L1014-IRS the lowest luminosity isolated protostar known, and an ideal target with which to test star formation theories at the low mass end. In order to study its molecular content and to search for the presence of a molecular outflow, we mapped L1014 in at least one transition of 12CO, N2H+, HCO+, CS and of their isotopologues 13CO, C18O, C17O, N2D+ and H13CO+, using the FCRAO, the IRAM 30 meter and the CSO. The data show physical and chemical properties in L1014 typical of the less evolved starless cores: i.e. H2 central density of a few 10^5 molecules cm^-3, estimated mass of ~2M_sun, CO integrated depletion factor less than 10, N(N2H+)~6*10^12 cm^-2, N(N2D+)/N(N2H+) equal to 10% and relatively broad N2H+(1--0) lines (0.35 km/s). Infall signatures and significant velocity shifts between optically thick and optically thin tracers are not observed in the line profiles. No classical signatures of molecular outflow are found in the 12CO and 13CO observations. In particular, no high velocity wings are found, and no well-defined blue-red lobes of 12CO emission are seen in the channel maps. If sensitive, higher resolution observations confirm the absence of an outflow on a smaller scale than probed by our observations, L1014-IRS would be the only protostellar object known to be formed without driving an outflow.

Using the Infrared Array Camera (IRAC) aboard the Spitzer Space Telescope, we have obtained mid-infrared photometry of the least massive known brown dwarf in the Chamaeleon I star-forming region. For this young brown dwarf, OTS 44, we have constructed a spectral energy distribution (SED) from 0.8-8 um by combining the measurements at 3.6, 4.5, 5.8, and 8.0 um from IRAC with ground-based photometry at I, J, H, and K. The resulting SED for OTS 44 exhibits significant excess emission longward of 3 um relative to the SED expected from the photosphere of the brown dwarf. We have successfully modeled the source of this excess emission in terms of an irradiated viscous accretion disk with M'<=10^-10 M_sun/year. With a spectral type of M9.5 and a mass of ~15 M_Jup, OTS 44 is now the coolest and least massive brown dwarf observed to have a circumstellar disk. These measurements demonstrate that disks exist around brown dwarfs even down to the deuterium burning mass limit and the approximate upper mass limit of extrasolar planetary companions.

We have undertaken a survey of N2H+ and N2D+ towards 31 low-mass starless cores using the IRAM 30m telescope. Our main objective has been to determine the abundance ratio of N2D+ and N2H+ towards the nuclei of these cores and thus to obtain estimates of the degree of deuterium enrichment, a symptom of advanced chemical evolution according to current models. We find that the N(N2D+)/N(N2H+) ratio is larger in more "centrally concentrated cores" with larger peak H2 and N2H+ column density than the sample mean. The deuterium enrichment in starless cores is presently ascribed to depletion of CO in the high density (> 3*10^4 cm-3) core nucleus. To substantiate this picture, we compare our results with observations in dust emission at 1.2 mm and in two transitions of C18O. We find a good correlation between deuterium fractionation and N(C18O)/N(H2) for the nuclei of 14 starless cores. We, thus, identified a set of properties that characterize the most evolved, or "pre-stellar", starless cores. These are: higher N2H+ and N2D+ column densities, higher N(N2D+)/N(N2H+), more pronounced CO depletion, broader N2H+ lines with infall asymmetry, higher central H2 column densities and a more compact density profile than in the average core. We conclude that this combination of properties gives a reliable indication of the evolutionary state of the core. Seven cores in our sample (L1521F, OphD, L429, L694, L183, L1544 and TMC2) show the majority of these features and thus are believed to be closer to forming a protostar than are the other members of our sample. Finally, we note that the subsample of Taurus cores behaves more homogeneously than the total sample, an indication that the external environment could play an important role in the core evolution.

Models of self-gravitating gas in the early stages of pressure-free collapse are compared for initial states which are equilibrium layers, cylinders, and Bonnor-Ebert spheres. For each geometrical case the density profile has an inner region of shallow slope surrounded by an outer region of steep slope, and the profile shape during early collapse remains similar to the profile shape of the initial equilibrium. The two-slope density structure divides the spherical collapse history into a starless infall phase and a protostellar accretion phase. The similarity of density profiles implies that Bonnor-Ebert fits to observed column density maps may not distinguish spherical cores from oblate or prolate cores, and may not distinguish static cores from collapsing cores. The velocity profiles discriminate better than the density profiles between initial geometries and between collapse ages. The infall velocity generally has a subsonic maximum value, which is approximately equal to the initial velocity dispersion times the ratio of collapse age to central free-fall time. Observations of starless core line profiles constrain collapse models. Collapse from initial states which are strongly condensed and slightly prolate is consistent with infall asymmetry observed around starless cores, and is more consistent than collapse from initial states which are weakly condensed, and/or oblate. Spherical models match observed inward speeds 0.05-0.09 km/s over 0.1-0.2 pc, if the collapse has a typical age 0.3-0.5 free fall times, and if it began from a centrally condensed state which was not in stable equilibrium. In a collapsing core, optically thin line profiles should broaden and develop two-peak structure as seen in L1544.

We present a study of analytic models of starless cores whose line profiles have ``infall asymmetry,'' or blue-skewed shapes indicative of contracting motions. We compare the ability of two types of analytical radiative transfer models to reproduce the line profiles and infall speeds of centrally condensed starless cores whose infall speeds are spatially constant and range between 0 and 0.2 km s-1. The model line profiles of HCO+ (J=1-0) and HCO+ (J=3-2) are produced by a self-consistent Monte Carlo radiative transfer code. The analytic models assume that the excitation temperature in the front of the cloud is either constant (``two-layer'' model) or increases inward as a linear function of optical depth (``hill'' model). Each analytic model is matched to the line profile by rapid least-squares fitting. The blue-asymmetric line profiles with two peaks, or with a blue shifted peak and a red shifted shoulder, can be well fit by the ``HILL5'' model (a five parameter version of the hill model), with an RMS error of 0.02 km s-1. A peak signal to noise ratio of at least 30 in the molecular line observations is required for performing these analytic radiative transfer fits to the line profiles.

Recent observations suggest that the behaviour of tracer species such as N_2H+ and CS is significantly different in regions of high and low mass star formation. In the latter, N_2H+ is a good tracer of mass, while CS is not. Observations show the reverse to be true in high-mass star formation regions. We use a computational chemical model to show that the abundances of these and other species may be significantly altered by a period of accelerated collapse in high mass star forming regions. We suggest these results provide a potential explanation of the observations, and make predictions for the behaviour of other species.

We present observations of L1014, a dense core in the Cygnus region previously thought to be starless, but data from the Spitzer Space Telescope shows the presence of an embedded source. We propose a model for this source that includes a cold core, heated by the interstellar radiation field, and a low-luminosity internal source. The low luminosity of the internal source suggests a substellar object. If L1014 is representative, other "starless" cores may turn out to harbor central sources.

Sensitive 5-38 um Spitzer Space Telescope (SST) and ground based 3-5 um spectra of the embedded low mass protostars B5 IRS1 and HH46 IRS show deep ice absorption bands superposed on steeply rising mid-infrared continua. The ices likely originate in the circumstellar envelopes. The CO2 bending mode at 15 um is a particularly powerful tracer of the ice composition and processing history. Toward these protostars, this band shows little evidence for thermal processing at temperatures above 50 K. Signatures of lower temperature processing are present in the CO and OCN- bands, however. The observed CO2 profile indicates an intimate mixture with H2O, but not necessarily with CH3OH, in contrast to some high mass protostars. This is consistent with the low CH3OH abundance derived from the ground based L band spectra. The CO2/H2O column density ratios are high in both B5 IRS1 and HH46 IRS (~35%). Clearly, the SST spectra are essential to study ice evolution in low mass protostellar environments, and to eventually determine the relation between interstellar and solar system ices.

We have used the Spitzer Infrared Spectrograph to observe seven members of the TW Hya association, the nearest stellar association whose age ($\sim$ 10 Myr) is similar to the timescales thought to apply to planet formation and disk dissipation. Only two of the seven targets display infrared excess emission, indicating that substantial amounts of dust still exist closer to the stars than is characteristic of debris disks; however, in both objects we confirm an abrupt short-wavelength edge to the excess, as is seen in disks with cleared-out central regions. The mid-infrared excesses in the spectra of Hen 3-600 and TW Hya include crystalline silicate emission features, indicating that the grains have undergone significant thermal processing. We offer a detailed comparison between the spectra of TW Hya and Hen 3-600, and a model that corroborates the spectral shape and our previous understanding of the radial structure of these protoplanetary disks.

We present the analysis of seven band (1.2 to 8 micron) ground and space-based imaging of the NGC 7129 young stellar cluster from FLAMINGOS on MMT, 2MASS, and the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. An analysis of the H-[4.5] vs. J-H colors reveals 84 objects with circumstellar disks. Of these, 42 are located within the cluster core, a 0.5 pc (100'') radius region of enhanced stellar surface density. From a luminosity and extinction limited sample of the stars within the cluster core boundary we have determined that 54% +/- 14% have circumstellar disks. Finally, we report the detection of several resolved outflows in the IRAC 4.5 micron mosaic.

We present observations of the star formation region NGC 7129 taken with the Multiband Imaging Photometer for Spitzer (MIPS). A significant population of sources, likely pre-main sequence members of the young stellar cluster, is revealed outside the central photoionization region. Combining with Infrared Array Camera (IRAC) and ground-based near-infrared images, we have obtained colors and spectral energy distributions for some 60 objects. The [3.6]-[4.5] vs. [8]-[24] color-color plane shows sources clustered at several different loci, which roughly correspond to the archetypal evolutionary sequence Class 0, I, II, and III. We obtain preliminary classifications for 36 objects, and find significant numbers of both Class I and II objects. Most of the pre-main sequence candidates are associated with the densest part of the molecular cloud surrounding the photoionization region, indicating active star formation over a broad area outside the central cluster. We discuss three Class II candidates that exhibit evidence of inner disk clearing, which would be some of the youngest known examples of a transition from accretion to optically thin quiescent disks.

We report initial results from IRAC observations of four young stellar clusters. These regions are part of a larger Spitzer survey of 31 young stellar groups and clusters within 1 kpc of the Sun. In each of the four clusters, there are between 39 and 85 objects with colors inconsistent with reddened stellar photospheres. We identify these objects as young stars with significant emission from circumstellar dust. Applying an analysis developed in a companion paper (Allen et al. 2004), we classify these objects as either pre-main sequence stars with disks (class II) or protostellar objects (class I). These show that the sites of recent star formation are distributed over multi-parsec size scales. In two clusters, Cepheus C and S140, we find protostars embedded in filamentary dark clouds seen against diffuse emission in the IRAC bands.

We report spectra obtained with the Spitzer Space Telescope in the wavelength range between 14 microns and 35 microns of 19 nearby main-sequence stars with infrared excesses. The six stars with strong dust emission show no recognizable spectral features, suggesting that the bulk of the emitting particles have diameters larger than 10 microns. If the observed dust results from collisional grinding of larger solids, we infer minimum masses of the parent body population between 0.004 of the Earth's mass and 0.06 of the Earth's mass. We estimate grain production rates of 10 Gg/s around lambda Boo and HR 1570; selective accretion of this matter may help explain their peculiar surface abundances. There appear to be inner truncations in the dust clouds at 48 AU, 11 AU, 52 AU and 54 AU around HR 333, HR 506, HR 1082 and HR 3927, respectively.

We compare the infrared colors predicted by theoretical models of protostellar envelopes and protoplanetary disks with initial observations of young stellar objects made with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope (Werner et al. 2004, Fazio et al. 2004). Disk and envelope models characterized by infall and/or accretion rates found in previous studies can quantitatively account for the range of IRAC colors found in four young embedded clusters: S 140, S 171, NGC 7129, and Cep C. The IRAC color-color diagram ([3.6]-[4.5] vs. [5.8]-[8.0]) can be used to help to distinguish between young stars with only disk emission and protostars with circumstellar envelopes.

We present CS(3-2) and DCO+(2-1) observations of 94 starless cores and compare the results with previous CS(2-1) and N2H+(1-0) observations to study inward motions in starless cores. The velocity shifts of the CS(3-2) and (2-1) lines with respect to N2H+ correlate well with each other and have similar distributions. This implies that, in many cores, systematic inward motions of gaseous material may occur over a range of density of at least a factor 4. We identify 18 infall candidates based on observations of CS(3-2), CS(2-1), DCO+(2-1) and N2H+(1-0). The eight best candidates, L1355, L1498, L1521F, L1544, L158, L492, L694-2, and L1155C-1, each show at least four indications of infall asymmetry and no counter-indications. Fits of the spectra to a 2-layer radiative transfer model in ten infall candidates suggest that the median effective line-of-sight speed of the inward-moving gas is ~0.07 km/s for CS (3-2) and ~0.04 km/s for CS(2-1). Considering that the optical depth obtained from the fits is usually smaller in CS(3-2) than in (2-1) line, this may imply that CS(3-2) usually traces inner denser gas in higher inward motions than CS(2-1). However, it is also possible that this conclusion is not representative of all starless core infall candidates, due to the statistically small number analyzed here. Further line observations will be useful to test this conclusion.

We observed the pre-stellar core L1521F in dust emission at 1.2mm and in two transitions each of N2H+, N2D+, C18O, and C17O in order to increase the sample of well studied centrally concentrated and chemically evolved starless cores, likely on the verge of star formation, and to determine the initial conditions for low--mass star formation in the Taurus Molecular Cloud. We derived in this object a molecular hydrogen number density n(H2) ~ 10^6 cm-3 and a CO depletion factor, integrated along the line of sight, fD ~ 15 in the central 20'', similar to the pre-stellar core L1544. However, the N(N2D+)/N(N2H+) column density ratio is ~0.1, a factor of about 2 lower than that found in L1544. The observed relation between the deuterium fractionation and the integrated CO depletion factor across the core can be reproduced by chemical models if N2H+ is slightly (factor of ~2 in fractional abundance) depleted in the central 3000 AU. The N2H+ and N2D+ linewidths in the core center are ~0.3 km/s, significantly larger than in other more quiescent Taurus starless cores but similar to those observed in the center of L1544. The kinematical behaviour of L1521F is more complex than seen in L1544, and a model of contraction due to ambipolar diffusion is only marginally consistent with the present data. Other velocity fields, perhaps produced by unresolved substructure, are present. Both chemical and kinematical analyses suggest that L1521F is less evolved than L1544, but, in analogy with L1544, it is approaching the ``critical'' state.

We have characterized the physical structure and chemical composition of two close-to-round starless cores in Taurus-Auriga, L1498 and L1517B. Our analysis is based on high angular resolution observations in at least two transitions of NH3, N2H+, CS, C34S, C18O, and C17O, together with maps of the 1.2 mm continuum. For both cores, we derive radial profiles of constant temperature and constant turbulence, together with density distributions close to those of non-singular isothermal spheres. Using these physical conditions and a Monte Carlo radiative transfer model, we derive abundance profiles for all species and model the strong chemical differentiation of the core interiors. According to our models, the NH3 abundance increases toward the core centers by a factor of several while N2H+ has a constant abundance over most of the cores. In contrast, both C18O and CS (and isotopomers) are strongly depleted in the core interiors, most likely due to their freeze out onto grains at densities of a few 10^4 cm^{-3}. Concerning the kinematics of the dense gas, we find (in addition to constant turbulence) a pattern of internal motions at the level of 0.1 km s^{-1}. These motions seem correlated with asymmetries in the pattern of molecular depletion, and we interpret them as residuals of core contraction. Their distribution and size suggest that core formation occurs in a rather irregular manner and with a time scale of a Myr. A comparison of our derived core properties with those predicted by supersonic turbulence models of core formation shows that our Taurus cores are much more quiescent than representative predictions from these models.

We present a near-infrared extinction study of the dark globule L694-2, a starless core that shows strong evidence for inward motions in molecular line profiles. The J,H, and K band data were taken using the European Southern Observatory New Technology Telescope. The best fit simple spherical power law model has index p=2.6 +/- 0.2, over the 0.036--0.1 pc range in radius sampled in extinction. This power law slope is steeper than the value of p=2 for a singular isothermal sphere, the initial condition of the inside-out model for protostellar collapse. Including an additional extinction component along the line of sight further steepens the inferred profile. Fitting a Bonnor-Ebert sphere results in a super-critical value of the dimensionless radius xi_max=25 +/- 3. The unstable configuration of material may be related to the observed inward motions. The Bonnor-Ebert model matches the shape of the observed profile, but significantly underestimates the amount of extinction (by a factor of ~4). This discrepancy in normalization has also been found for the nearby protostellar core B335 (Harvey et al. 2001). A cylindrical density model with scale height H=0.0164+/- 0.002 pc viewed at a small inclination to the cylinder axis provides an equally good radial profile as a power law model, and reproduces the asymmetry of the core remarkably well. In addition, this model provides a basis for understanding the discrepancy in the normalization of the Bonnor-Ebert model, namely that L694-2 has prolate structure, with the full extent (mass) of the core being missed by assuming symmetry between the profiles in the plane of the sky and along the line-of-sight. If the core is sufficiently magnetized then fragmentation may be avoided, and later evolution might produce a protostar similar to B335.

We present a catalog of near-infrared surveys of young (<= a few 10^6 yr) stellar groups and clusters within 1 kpc from the Sun, based on an extensive search of the literature from the past ten years. We find 143 surveys from 69 published articles, covering 73 different regions. The number distribution of stars in a region has a median of 28 and a mean of 100. About 80% of the stars are in clusters with at least 100 members. By a rough classification of the groups and clusters based on the number of their associated stars, we show that most of the stars form in large clusters. The spatial distribution of cataloged regions in the Galactic plane shows a relative lack of observed stellar groups and clusters in the range 270 deg < l < 60 deg of Galactic longitude, reflecting our location between the Local and Sagittarius arms. This compilation is intended as a useful resource for future studies of nearby young regions of multiple star formation.

We study the density structure of candidate contracting starless core L694-2 using 1.3 mm dust continuum observations from the IRAM Plateau de Bure Interferometer and the Berkeley-Illinois-Maryland Array, probing spatial scales from 10000-500 AU. The long baseline PdBI observations detect no emission, and limit the contamination from a compact component F_c < 2.7 mJy. This limit corresponds to a very small disk mass, M_disk < 5e-4 M_sun x (60 K / T_disk), and bolsters the ``starless'' interpretation of the core. The shorter baseline BIMA data are compared to density models using a physically motivated temperature distribution with a central minimum. This analysis provides clear evidence for a turn-over from the steep outer density profile observed in dust extinction to much more shallow behavior in the inner regions (<7500 AU). The best fit Bonnor-Ebert, Plummer-like, broken power law, and end-on cylinder models produce very similar flattened profiles and cannot be distinguished. We quantify the sensitivity of the inferred structure to various uncertainties, including the temperature distribution, the central position, and the presence of a weak unresolved central component. The largest uncertainty comes from the temperature assumption; an isothermal model modifies the best fit parameters by \~2 sigma, with the inferred profiles more shallow. Dust emission and extinction profiles are reproduced by an embedded cylinder with scale height H=13.5'' inclined at a small angle to the line of sight. The flat central density profile suggests that pressure forces still support the core, despite the extended inward motions inferred from molecular line observations (Lee, Myers, & Tafalla 2001). In the context of the cylindrical density model, these inward motions may represent the contraction of a prolate core along its major axis.

We present subarcsecond resolution dust continuum observations of the protostellar collapse candidate B335 made with the IRAM Plateau de Bure Interferometer at wavelengths of 1.2 and 3.0 mm. These observations probe to < 100 AU size scales and reveal a compact source component that we identify with a circumstellar disk. We analyze these data in concert with previous lower resolution interferometer observations and find a best fit density structure for B335 that consists of a power law envelope with index p=1.55 +/- 0.04 (r < 5000 AU) together with a disk (r < 100 AU) of flux F_{1.2 mm}=21 +/-2 mJy. We estimate a systematic uncertainty in the power law index delta(p) < 0.15, where the largest error comes from the assumed form of the dust temperature falloff with radius. This determination of the inner density structure of B335 has a precision unique amongst protostellar cores, and it is consistent with the r^{-1.5} profile of gravitational free-fall, in accord with basic expectations for the formation of a star. The flux (and implied mass) of the compact component in B335 is typical of the disks around T Tauri stars.

Deep JHKs photometry was obtained towards eight dense molecular cores and J-H vs. H-Ks color-color plots are presented. Our photometry, sensitive to the detection of a 1 solar mass, 1 X 10^6 year old star through approx. 35 - 50 magnitudes of visual extinction, shows no indication of the presence of star/disk systems based on J-H vs. H-Ks colors of detected objects. The stars detected towards the cores are generally spatially anti-correlated with core centers suggesting a background origin, although we cannot preclude the possibility that some stars detected at H and Ks alone, or Ks alone, are not low mass stars or brown dwarfs (< 0.3 Solar Masses) behind substantial amounts of visual extinction (e.g. 53 magnitudes for L183B). Lower limits to optical extinctions are estimated for the detected background stars, with high extinctions being encountered, in the extreme case ranging up to at least Av = 46, and probably higher. The extinction data are used to estimate cloud masses and densities which are comparable to those determined from molecular line studies. Variations in cloud extinctions are consistent with a systematic nature to cloud density distributions and column density variations and extinctions are found to be consistent with submillimeter wave continuum studies of similar regions. The results suggest that some cores have achieved significant column density contrasts (approx. 30) on sub-core scales (approx. 0.05 pc) without having formed known stars.

We present the results of N2H+(1-0) observations of 35 dense molecular cloud cores from the northern and southern hemispheres where massive stars and star clusters are formed. Line emission has been detected in 33 sources, for 28 sources detailed maps have been obtained. The optical depth of (23-12) component toward peak intensity positions of 10 sources is ~ 0.2-1. In total, 47 clumps have been revealed in 26 sources. Integrated intensity maps with aspect ratios < 2 have been fitted with a power-law radial distribution $r^{-p}$ convolved with the telescope beam. Mean power-law index is close to unity corresponding to the $\sim r^{-2}$ density profile provided N2H+ excitation conditions do not vary inside these regions. Line widths of the cores either decrease or stay constant with distance from the center. The ratio of rotational to gravitational energy is too low for rotation to play a significant role in the dynamics of the cores. A correlation between mean line widths and sizes of clumps has been found.

As the first step of a multi-wavelength investigation into the relationship between young stellar clusters and their environment we present fully-sampled maps in the J=1--0 lines of 13CO and C18O and the J=2--1 line of C18O for a selected group of thirty young stellar groups and clusters within 1kpc of the Sun. This is the first systematic survey of these regions to date. The clusters range in size from several stars to a few hundred stars. Thirty fields ranging in size from 8'x 8' to 60'x 30' were mapped with 47'' resolution simultaneously in the J=1-0 lines with Five College Radio Observatory. Seventeen sources were mapped over fields ranging in size from 3'x 3' to 13'x 13' in the J=2--1 line with 35'' resolution with the Heinrich Hertz Telescope. We compare the cloud properties derived from each of the three tracers in order to better understand systematic uncertainties in determing masses and linewidths. Cloud masses are determined independently using the 13CO and C18O transitions, these masses range from 30 to 4000 M_sun. Finally, we present a simple morphological classification scheme which may serve as a rough indicator of cloud evolution.

We present a study of the density structure of the protostellar collapse candidate B335 using continuum observations from the IRAM Plateau de Bure Interferometer made at wavelengths of 1.2mm and 3.0mm. We analyze these data, which probe spatial scales from 5000 AU to 500 AU, directly in the visibility domain by comparison to synthetic observations constructed from models that assume different physical conditions. This approach allows for much more stringent constraints to be derived from the data than from analysis of images. A single radial power law in density provides a good description of the data, with best fit power law index p=1.65+/-0.05. Through simulations, we quantify the sensitivity of this result to various model uncertainties, including assumptions of temperature distribution, outer boundary, dust opacity spectral index, and an unresolved central component. The largest uncertainty comes from the unknown presence of a centralized point source. A point source with 1.2mm flux of F=12+/-7 mJy reduces the density index to p=1.47+/-0.07. The remaining sources of systematic uncertainty, the most important of which is the temperature distribution, likely contribute a total uncertainty of < 0.2. We therefore find strong evidence that the power law index of the density distribution within 5000 AU is significantly less than the value at larger radii, close to 2.0 from previous studies of dust emission and extinction. These results conform well to the generic paradigm of isolated, low-mass star formation which predicts a power law density index close to p=1.5 for an inner region of gravitational free fall onto the protostar.

We used the NRAO Very Large Array to search for 3.6 cm continuum emission from embedded protostars in a sample of 8 nearby ``starless'' cores that show spectroscopic evidence for infalling motions in molecular emission lines. We detect a total of 13 compact sources in the eight observed fields to 5 sigma limiting flux levels of typically 0.09 mJy. None of these sources lie within 1' of the central positions of the cores, and they are all likely background objects. Based on an extrapolation of the empirical correlation between the bolometric luminosity and 3.6 cm luminosity for the youngest protostars, these null-detections place upper limits of ~0.1 L_sun (d/140pc)^2 on the luminosities of protostellar sources embedded within these cores. These limits, together with the extended nature of the inward motions inferred from molecular line mapping (Lee et al. 2001), are inconsistent with the inside-out collapse model of singular isothermal spheres and suggest a less centrally condensed phase of core evolution during the earliest stages of star formation.

We present results of an extensive mapping survey of N2H+(1-0) in about 60 low mass cloud cores already mapped in the NH3(1,1) inversion transition line. The survey has been carried out at the FCRAO antenna with an angular resolution about 1.5 times finer than the previous ammonia observations. Cores with stars typically have map sizes about a factor of two smaller for N2H+ than for NH3, indicating the presence of denser and more centrally concentrated gas compared to starless cores. Significant correlations are found between NH3 and N2H+ column densities and excitation temperatures in starless cores, but not in cores with stars, suggesting a different chemical evolution of the two species. Velocity gradients range between 0.5 and 6 km/s/pc, similar to what has been found with NH3 data. ``Local'' velocity gradients show significant variation in both magnitude and direction, suggesting the presence of complexmotions not interpretable as simple solid body rotation. Integrated intensity profiles of starless cores present a ``central flattening'' and are consistent with a spherically symmetric density law n ~ r^{-1.2} for r < ~0.03 pc and n ~ r^{-2} at larger r. Cores with stars are better modelled with single density power laws with n ~ r^{-2}. Line widths change across the core but we did not find a general trend. The deviation in line width correlates with the mean line width, suggesting that the line of sight contains ~ 10 coherence lengths. The corresponding value of the coherence length, ~ 0.01 pc, is similar to the expected cutoff wavelength for MHD waves. This similarity may account for the increased ``coherence'' of line widths on small scales. Despite of the finer angular resolution, the majority of N2H+ and NH3 maps show a similar ``simple'' structure, with single peaks and no elongation.

(Abridged) We present evidence that low-mass starless cores, the simplest units of star formation, are systematically differentiated in their chemical composition. Molecules including CO and CS almost vanish near the core centers, where the abundance decreases by one or two orders of magnitude. At the same time, N2H+ has a constant abundance, and the fraction of NH3 increases toward the core center. Our conclusions are based on a study of 5 mostly-round starless cores (L1498, L1495, L1400K, L1517B, and L1544), which we have mappedin C18O(1-0), C17O(1-0), CS(2-1), C34S(2-1), N2H+(1-0), NH3(1,1) and (2,2), and the 1.2 mm continuum. For each core we have built a model that fits simultaneously the radial profile of all observed emission and the central spectrum for the molecular lines. The observed abundance drops of CO and CS are naturally explained by the depletion of these molecules onto dust grains at densities of 2-6 10^4 cm-3. N2H+ seems unaffected by this process up to densities of several 10^5, while the NH3 abundance may be enhanced by reactions triggered by the disappearance of CO from the gas phase. With the help of our models, we show that chemical differentiation automatically explains the discrepancy between the sizes of CS and NH3 maps, a problem which has remained unexplained for more than a decade. Our models, in addition, show that a combination of radiative transfer effects can give rise to the previously observed discrepancy in the linewidth of these two tracers. Although this discrepancy has been traditionally interpreted as resulting from a systematic increase of the turbulent linewidth with radius, our models show that it can arise in conditions of constant gas turbulence.

We present a catalogue of near-infrared photometry of young stars associated with the Ophiuchus molecular cloud, based on observations made with the Hubble Space Telescope NICMOS-3 camera at 1.1 $\mu$m and 1.6 $\mu$m. Our survey covers 0.02 square degrees centered on the dense molecular cores in Lynds 1688. We detect 165 sources at 1.6 $\mu$m and 65 sources at 1.1 $\mu$m, within our estimated completeness limits of 21.0 mag and 21.5 mag, respectively. An analysis of the cloud extinction, based on existing molecular line maps, suggests that most of the sources lying within the 40 $\rm{A_V}$ extinction contour of the cloud are probable cloud members. Approximately half (58/108) of these sources are previously unpublished. The faint embedded sources revealed by these observations are spatially concentrated in three regions of high stellar space density (N$> 10^4$ stars pc$^{-3}$). While the spatial distribution of these sources reflects that of the brighter, well--known population of young stars in Ophiuchus, it is distinctly different from the distribution of cool concentrations seen in the submillimeter. Seven new brown dwarf candidates are identified, based on their infrared colors and their projected locations on high column-density regions of the molecular cloud. Eight new candidate binary and five new candidate triple systems, having separations between 0\farcs2 to 10\arcsec (29 to 1450 AU) are reported. The spatial resolution and sensitivity of these observations reveal five apparent disk/envelope systems seen via scattered light, and four nebulous objects with complex morphologies.

The maps presented in Paper I are here used to infer the variation of the column densities of HCO+, DCO+, N2H+, and N2D+ as a function of distance from the dust peak. These results are interpreted with the aid of a crude chemical model which predicts the abundances of these species as a function of radius in a spherically symmetric model with radial density distribution inferred from the observations of dust emission at millimeter wavelengths and dust absorption in the infrared. Our main observational finding is that the N(N2D+)/N(N2H+) column density ratio is of order 0.2 towards the L1544 dust peak as compared to N(DCO+)/N(HCO+) = 0.04. We conclude that this result as well as the general finding that N2H+ and N2D+ correlate well with the dust is caused by CO being depleted to a much higher degree than molecular nitrogen in the high density core of L1544. Depletion also favors deuterium enhancement and thus N2D+, which traces the dense and highly CO-depleted core nucleus, is much more enhanced than DCO+. Our models do not uniquely define the chemistry in the high density depleted nucleus of L1544 but they do suggest that the ionization degree is a few times 10^{-9} and that the ambipolar diffusion time scale is locally similar to the free fall time. It seems likely that the lower limit which one obtains to ionization degree by summing all observable molecular ions is not a great underestimate of the true ionization degree. We predict that atomic oxygen is abundant in the dense core and, if so, H3O+ may be the main ion in the central highly depleted region of the core.

We have mapped the dense dark core L1544 in H13CO+(1-0), DCO+(2-1), DCO+(3-2), N2H+(1-0), NTH+(3-2), N2D+(2-1), N2D+(3-2), C18O(1-0), and C17O(1-0) using the IRAM 30-m telescope. We have obtained supplementary observations of HC18O+(1-0), HC17O+(1-0), and D13CO+(2-1). Many of the observed maps show a general correlation with the distribution of dust continuum emission in contrast to C18O(1-0) and C17O(1-0) which give clear evidence for depletion of CO at positions close to the continuum peak. In particular N2D+(2-1) and (3-2) and to a lesser extent N2H+(1-0) appear to be excellent tracers of the dust continuum. We find that the tracers of high density gas (in particular N2D+) show a velocity gradient along the minor axis of the L1544 core and that there is evidence for larger linewidths close to the dust emission peak. We interpret this using the model of the L1544 proposed by Ciolek & Basu (2000) and by comparing the observed velocities with those expected on the basis of their model. The results show reasonable agreement between observations and model in that the velocity gradient along the minor axis and the line broadening toward the center of L1544 are predicted by the model. This is evidence in favour of the idea that amipolar diffusion across field lines is one of the basic processes leading to gravitational collapse. However, line widths are significantly narrower than observed and are better reproduced by the Myers & Zweibel (2001) model which considers the quasistatic vertical contraction of a layer due to dissipation of its Alfvenic turbulence, indicating the importance of this process for cores in the verge of forming a star.

We present a near-infrared extinction study of the dark globule B335, a protostellar collapse candidate, using data from HST/NICMOS and the W.M. Keck Observatory. These data allow a new quantitative test of the "inside-out" collapse model previously proposed to explain molecular line profiles observed toward this region. We find that the shape of the density profile is well matched by the collapse model, but that the amount of extinction corresponds to larger column densities than predicted. An unstable Bonnor-Ebert sphere with dimensionless outer radius xi_max=12.5+/-2.6 provides an equally good description of the density profile, and is indistinguishable from the collapse model over the range in radius sampled by the data. The bipolar outflow driven by the embedded YSO has an important effect on the extinction through the core, and modeling the outflow as a hollowed-out bipolar cone of constant opening angle provides a good match to theobservations. The complete exinction map is well reproduced by a model that includes both infall and outflow, and an extra 20% dispersion that likely results from residual turbulent motions. This fitted model has an infall radius of R_inf=26"+/-3"(0.031pc for 250pc distance), and an outflow cone semi-opening angle of alpha=41+/-2 degrees. The fitted infall radius is consistent with those derived from molecular line observations and supports the inside-out collapse interpretation of the density structure. The fitted opening angle of the outflow is slightly larger than observed in CO emission, perhaps because the full extent of the outflow cone in CO becomes confused with ambient core emission at low velocities.

Millimeter wavelength observations are presented of NGC 1333 IRAS 4, a group of highly-embedded young stellar objects in Perseus, that reveal motions of infall, outflow, rotation, and turbulence in the dense gas around its two brightest continuum objects, 4A and 4B. These data have finest angular resolution of approximately 2" (0.0034 pc) and finest velocity resolution of 0.13 km/s. Infall motions are seen from inverse P-Cygni profiles observed in H2CO 3_12-2_11 toward both objects, but also in CS 3-2 and N2H+ 1-0 toward 4A, providing the least ambiguous evidence for such motions toward low-mass protostellar objects. Outflow motions are probed by bright line wings of H2CO 3_12-2_11 and CS 3-2 observed at positions offset from 4A and 4B, likely tracing dense cavity walls. Rotational motions of dense gas are traced by a systematic variation of the N2H+ line velocities, and such variations are found around 4A but not around 4B. Turbulent motions appear reduced with scale, given N2H+ line widths around both 4A and 4B that are narrower by factors of 2 or 3 than those seen from single-dish observations. Minimum observed line widths of approximately 0.2 km/s provide a new low, upper bound to the velocity dispersion of the parent core to IRAS 4, and demonstrate that turbulence within regions of clustered star formation can be reduced significantly. A third continuum object in the region, 4B', shows no detectable line emission in any of the observed molecular species.

We present the results of an extensive mapping survey of 53 `starless' cores in the optically thick line of CS 2-1 and the optically thin lines of N2H+ 1-0 and C18O 1-0. The purpose of this survey was to search for signatures of extended inward motions. This study finds 10 `strong' and 9 `probable' infall candidates, based on $\delta V_{CS}$ analysis and on the spectral shapes of CS lines. From our analysis of the blue-skewed CS spectra and the $\delta V_{CS}$ parameter, we find typical infall radii of 0.06-0.14 pc. Also, using a simple two layer radiative transfer model to fit the profiles, we derive one-dimensional infall speeds, half of whose values lie in the range of 0.05-0.09 km s$^{-1}$. These values are similar to those found in L1544 by Tafalla et al., and this result confirms that infall speeds in starless cores are generally faster than expected from ambipolar diffusion in a strongly sub-critical core. In addition, the observed infall regions are too extended to be consistent with the `inside-out' collapse model applied to a very low-mass star. In the largest cores, the spatial extent of the CS spectra with infall asymmetry is larger than the extent of the $\rm N_2H^+$ core by a factor of 2-3. All these results suggest that extended inward motions are a common feature in starless cores, and that they could represent a necessary stage in the condensation of a star-forming dense core.

We present the results of a new survey of 23 molecular clouds for the Zeeman effect in OH undertaken with the ATNF Parkes 64-m radio telescope and the NRAO Green Bank 43-m radio telescope. The Zeeman effect was clearly detected in the cloud associated with the HII region RCW 38, with a field strength of 38+/-3 micro-Gauss, and possibly detected in a cloud associated with the HII region RCW 57, with a field strength of -203+/-24 micro-Gauss. The remaining 21 measurements give formal upper limits to the magnetic field strength, with typical 1-sigma sensitivities <20 micro-Gauss. For 22 of the molecular clouds we are also able to determine thecolumn density of the gas in which we have made a sensitive search for the Zeeman effect. We combine these results with previous Zeeman studies of 29 molecular clouds, most of which were compiled by Crutcher (1999), for a comparsion of theoretical models with the data. This comparison implies that if the clouds can be modeled as initially spherical with uniform magnetic fields and densities that evolve to their final equilibrium state assuming flux-freezing then the typical cloud is magnetically supercritical, as was found by Crutcher (1999). If the clouds can be modeled as highly flattened sheets threaded by uniform perpendicular fields, then the typical cloud is approximately magnetically critical, in agreement with Shu et al. (1999), but only if the true values of the field for the non-detections are close to the 3-sigma upper limits. If instead these values are significantly lower (for example, similar to the 1-sigma limits), then the typical cloud is generally magnetically supercritical.

This paper argues that star forming environments should be classified into finer divisions than the traditional isolated and clustered modes. Using the observed set of galactic open clusters and theoretical considerations regarding cluster formation, we estimate the fraction of star formation that takes place within clusters. We find that less than 10% of the stellar population originates from star forming regions destined to become open clusters, confirming earlier estimates. The smallest clusters included in the observational surveys (having at least N=100 members) roughly coincide with the smallest stellar systems that are expected to evolve as clusters in a dynamical sense. We show that stellar systems with too few members N < N_\star have dynamical relaxation times that are shorter than their formation times (1-2 Myr), where the critical number of stars N_\star \approx 100. Our results suggest that star formation can be characterized by (at least) three principal modes: I. isolated singles and binaries, II. groups (N<N_\star), and III. clusters (N>N_\star). Many -- if not most -- stars form through the intermediate mode in stellar groups with 10<N<100. Such groups evolve and disperse much more rapidly than open clusters; groups also have a low probability of containing massive stars and are unaffected by supernovae and intense ultraviolet radiation fields. Because of their short lifetimes and small stellar membership, groups have relatively little effect on the star formation process (on average) compared to larger open clusters.

We present molecular-line observations toward the dense core in L483 and its bipolar outflow powered by the Class 0 object IRAS 18148-0440. CO (carbon monoxide) maps show that the outflow is well collimated and asymmetric, and that its gas is warmer than the surrounding cloud by at least a factor of 2. H2CO (formaldehyde) lines toward the outflow show prominent high-velocity wings and evidence for an H2CO abundance enhancement of a factor of 20. At ambient velocities, these lines show strong self-absorption and a brighter blue peak, a characteristic signature of inward motions. Finally, and in contrast with the outflows from other Class 0 objects, the CH3OH (methanol) lines in L483 do not show high velocity wings, and no evidence for abundance enhancement is found in this molecule or in SiO (silicon monoxide). Comparing the physical and chemical properties of the outflow in L483 with those of other outflows from Class 0 and Class I sources, we find that the L483 outflow is somewhat intermediate between these types. This suggests that the L483 central source has already started its transition between Class 0 and Class I, and that its mixed properties illustrate how this transition occurs.

We have mapped the dense gas distribution and dynamics in the NW region of the Serpens molecular cloud in the CS(2-1) and N2H+(1-0) lines and 3 mm continuum using the FCRAO telescope and BIMA interferometer. 7 continuum sources are found. The N2H+ spectra are optically thin and fits to the 7 hyperfine components are used to determine the distribution of velocity dispersion. 8 cores, 2 with continuum sources, 6 without, lie at a local linewidth minimum and optical depth maximum. The CS spectra are optically thick and generally self-absorbed over the full 0.2 pc extent of the map. We use the line wings to trace outflows around at least 3, and possibly 4, of the continuum sources, and the asymmetry in the self-absorption as a diagnostic of relative motions between core centers and envelopes. The quiescent regions with low N2H+ linewidth tend to have more asymmetric CS spectra than the spectra around the continuum sources indicating higher infall speeds. These regions have typical sizes ~5000 AU, linewidths ~0.5 km/s, and infall speeds ~0.05 km/s. The correlation of CS asymmetry with N2H+ velocity dispersion suggests that the inward flows of material that build up pre-protostellar cores are driven at least partly by a pressure gradient rather than by gravity alone. We discuss a scenario for core formation and eventual star forming collapse through the dissipation of turbulence.

We have observed HCO+ J=3-2 toward 16 Class I sources and 18 Class 0 sources, many of which were selected from Mardones et al. (1997). Eight sources have profiles significantly skewed to the blue relative to optically thin lines. We suggest six sources as new infall candidates. We find an equal "blue excess" among Class 0 and Class I sources after combining this sample with that of Gregersen et al. (1997). We used a Monte Carlo code to simulate the temporal evolution of line profiles of optically thick lines of HCO+, CS and H2CO in a collapsing cloud and found that HCO+ had the strongest asymmetry at late times. If a blue-peaked line profile implies infall, then the dividing line between the two classes does not trace the end of the infall stage.

In this paper, we present the first results of a survey of 220 starless cores selected primarily by their optical obscuration (Lee & Myers 1999) and observed in CS (2-1), N$_2$H$^+$ (1-0), and C$^{18}$O (1-0) using the NEROC Haystack 37-m telescope, to quantify the incidence of infall motions toward starless cores and find infall candidates for further study in the future.

We present combined single-dish and interferometric CS(2--1) and N2H+(1--0) observations of a compact core in the NW region of the Serpens molecular cloud. The core is starless according to observations from optical to millimeter wavelengths and its lines have turbulent widths and ``infall asymmetry''. Line profile modeling indicates supersonic inward motions v_in>0.34 km/s over an extended region L>12000AU. The high infall speed and large extent exceeds the predictions of most thermal ambipolar diffusion models and points to a more dynamical process for core formation. A short (dynamic) timescale, ~1e5 yr=L/v_in, is also suggested by the low N2H+ abundance ~1e-10.

We present a new catalogue of 406 dense cores optically selected by using the STScI Digitized Sky Survey (DSS). In this catalogue 306 cores have neither an Embedded YSO (EYSO) nor a Pre-Main-Sequence (PMS) star, 94 cores have EYSOs (1 core has both an EYSO and a PMS star), and 6 cores have PMS star only. Our sample of dense cores in the catalogue is fairly complete within a category of northern Lynds class 5, 6 clouds, and southern Hartley et al. (1986)'s class A clouds, providing a database useful for the systematic study of dense cores. Most of the cores listed in the catalogue have diameters between $0.05 - 0.36$ pc with a mean of $\sim 0.24$ pc. The sizes ($\sim 0.33$ pc in the mean) of cores with EYSOs are found to be usually larger than the sizes ($\sim 0.22 $ pc in the mean) of starless cores. The typical mean gas density of the cores is $\sim7\times 10^3 cm^{-3}$. Most of the cores are more likely elongated than spherical (mean aspect ratio: $\sim 2.4$). The ratio of the number of cores with EYSOs to the number of starless cores for our sample is about 0.3, suggesting that the typical lifetime of starless cores is $0.3-1.6$ Myr, about 3 times longer than the duration of the Class 0 and Class I phases. This lifetime is shorter than expected from models of ambipolar diffusion, by factors of 2-44.

We present interferometric observations of N2H+(1--0) in the starless, dense core L1544 in Taurus. Red-shifted self-absorption, indicative of inward motions, is found toward the center of an elongated core. The data are fit by a non-spherical model consisting of two isothermal, rotating, centrally condensed layers. Through a hybrid global-individual fit to the spectra, we map the variation of infall speed at scales ~1400AU and find values ~0.08 km/s around the core center. The inward motions are small in comparison to thermal, rotational, and gravitational speeds but are large enough to suggest that L1544 is very close to forming a star.

We propose that inward, subsonic flows arise from the local dissipation of turbulent motions in molecular clouds. Such "turbulent cooling flows" may account for recent observations of spatially extended inward motions towards dense cores. These pressure-driven flows may arise from various types of turbulence and dissipation mechanisms. For the example of MHD waves and turbulence damped by ion-neutral friction, sustained cooling flow requires that the outer gas be sufficiently turbulent, that the inner gas have marginal field-neutral coupling, and that this coupling decrease sufficiently rapidly with increasing density. These conditions are most likely met at the transition between outer regions ionized primarily by UV photons and inner regions ionized primarily by cosmic rays. If so, turbulent cooling flows can help form dense cores, with speeds faster than expected for ambipolar diffusion. Such motions could reduce the time needed for dense core formation and could precede and enhance the motions of star-forming gravitational infall.

We present an observational and theoretical study of the ionization fraction in several massive cores located in regions that are currently forming stellar clusters. Maps of the emission from the J = 1-> O transitions of C18O, DCO+, N2H+, and H13CO+, as well as the J = 2 -> 1 and J = 3 -> 2 transitions of CS, were obtained for each core. Core densities are determined via a large velocity gradient analysis with values typically 10^5 cm^-3. With the use of observations to constrain variables in the chemical calculations we derive electron fractions for our overall sample of 5 cores directly associated with star formation and 2 apparently starless cores. The electron abundances are found to lie within a small range, -6.9 < log10(x_e) < -7.3, and are consistent with previous work. We find no difference in the amount of ionization fraction between cores with and without associated star formation activity, nor is any difference found in electron abundances between the edge and center of the emission region. Thus our models are in agreement with the standard picture of cosmic rays as the primary source of ionization for molecular ions. With the addition of previously determined electron abundances for low mass cores, and even more massive cores associated with O and B clusters, we systematically examine the ionization fraction as a function of star formation activity. This analysis demonstrates that the most massive sources stand out as having the lowest electron abundances (x_e < 10^-8).

We report observations of 47 candidate protostars in two optically thick lines (H2CO 212-111 and CS 2-1) and one optically thin line (N2H+ 1-0) using the IRAM 30-m, SEST 15-m, and Haystack 37-m radio telescopes. The sources were selected from the redness of their spectra (Tbol < 200 K) and their near distance (d < 400 pc). Most of the sources have asymmetric optically thick lines. The observed distribution of velocity differences, deltaV = (Vthick-Vthin)/FWHMthin, is skewed toward negative (blue-shifted) velocities for both the H2CO and CS samples. This excess is much more significant for Class 0 than for Class I sources, suggesting that we detect infall motions toward Class 0 and not toward Class I sources. This indicates a difference in the physical conditions in the circumstellar envelopes around Class I and Class 0 sources, but does not rule out the presence of infall onto Class I sources by e.g. lower opacity gas. Bipolar outflows alone, or rotation alone, cannot reproduce these statistics if the sample of sources has randomly oriented symmetry axes. We identify 15 spectroscopic infall candidates, of which 6 are new. Most of these infall candidates have primarily turbulent rather than thermal motions, and are associated with clusters rather than being isolated.

A quantitative analysis of grain alignment in the filamentary dark cloud L1755 in Ophiuchus is presented. We show, that the observed decrease of the polarization to extinction ratio for the inner parts of this quiescent dark cloud can be explained as a result of the decrease of the efficiency of grain alignment.We make quantitative estimates of grain alignment efficiency for six mechanisms, involving grains with either thermal or suprathermal rotation, interacting with either magnetic field or gaseous flow. We also make semiquantitative estimates of grain alignment by radiative torques.We show that in conditions typical of dark cloud interiors, all known major mechanisms of grain alignment fail. All the studied mechanisms predict polarization at least an order of magnitude below the currently detectable levels of order of 1%. On the contrary, in the dark cloud environments where A_v<1 the grain alignment can be much more efficient. There the alignment of suprathermally rotating grains with superparamagentic inclusions, and possibly also radiative torques account for observed polarization. These results apply to L1755, which we model in detail and probably also to B216 and other similar dark clouds. Our study suggests an explanation for the difference in results obtained through polarimetry of background starlight and polarized thermal emission from the dust itself. We conjecture that the emission polarimetry selectively reveals aligned grains in the environment far from thermodynamic equilibrium as opposed to starlight polarization studie that probe the alignment of grains all the way along the line of sight, including the interiors of dark quiescent clouds, where no alignment is possible.

Trying to obtain a more complete picture of star forming regions in the Taurus Molecular Cloud, four newly discovered dense cores (L1521D, L1521F, L1524, L1507A) are identified and mapped through the ammonia (J,K)=(1,1) and (2,2) rotational inversion lines. These cores have sizes of 0.06-0.09 pc, hydrogen densities of order 10000 - 100000 cm-3 and kinetic temperatures between 8 and 10 K. The masses range from 0.2 - 1.0 solar masses, placing the cores with the lowest M values at the lower edge of the mass distribution for ammonia cores in the TMC. Turbulent, thermal and gravitational energies have been estimated. A comparison between these energy terms and considerations related to thermal and turbulent line broadening indicate that these cores are close to equilibrium. Moreover we report the detection of the J=9-8 transition of HC5N towards the four ammonia peak positions. The HC5N column densities are in agreement with values derived for other molecular cores in the TMC.