With the Stratospheric Observatory for Infrared Astronomy (SOFIA) routinely operating science flights, we demonstrate that observations with the Faint Object infraRed CAmera for the SOFIA Telescope (FORCAST) can provide reliable estimates of the internal luminosities, $L_{\rm int}$, of protostars. We have developed a technique to estimate $L_{\rm int}$ using a pair of FORCAST filters: one "short-wavelength" filter centered within 19.7-25.3 $\mu$m, and one "long-wavelength" filter within 31.5-37.1 $\mu$m. These $L_{\rm int}$ estimates are reliable to within 30-40% for 67% of protostars and to within a factor of 2.3-2.6 for 99% of protostars. The filter pair comprised of F25.3$\mu$m and F37.1$\mu$m achieves the best sensitivity and most constrained results. We evaluate several assumptions that could lead to systematic uncertainties. The OH5 dust opacity matches observational constraints for protostellar environments best, though not perfectly; we find that any improved dust model will have a small impact of 5-10% on the $L_{\rm int}$ estimates. For protostellar envelopes, the TSC84 model yields masses that are twice those of the Ulrich model, but we conclude this mass difference does not significantly impact results at the mid-infrared wavelengths probed by FORCAST. Thus, FORCAST is a powerful instrument for luminosity studies targeting newly discovered protostars or suspected protostars lacking detections longward of 24 $\mu$m. Furthermore, with its dynamic range and greater angular resolution, FORCAST may be used to characterize protostars that were either saturated or merged with other sources in previous surveys using the Spitzer Space Telescope or Herschel Space Observatory.

Kilonovae represent an important electromagnetic counterpart for compact binary mergers, which could become the most commonly detected gravitational wave (GW) source. Follow-up observations, triggered by GW events, of kilonovae are nevertheless difficult due to poor localization by GW detectors and due to their faint near-infrared peak emission that has limited observational capability. We show that the Near-Infrared Camera (NIRCam) on the James Webb Space Telescope (JWST) will be able to detect kilonovae within the relevant GW-detection range of $\sim$ 200 Mpc in short ($\lesssim$ 12-second) exposure times for a week following the merger. Despite this sensitivity, a kilonova search fully covering a fiducial localized area of $10$ $\mbox{deg}^2$ will not be viable with NIRCam due to its limited field of view. However, targeted surveys may be developed to optimize the likelihood of discovering kilonovae efficiently within limited observing time. We estimate that a survey of $10$ $\mbox{deg}^2$ focused on galaxies within 200 Mpc would require about 13 hours, dominated by overhead times; a survey further focused on galaxies exhibiting high star-formation rates would require $\sim$ 5 hours. The characteristic time may be reduced to as little as $\sim$4 hours, without compromising the likelihood of detecting kilonovae, by surveying sky areas associated with 50%, rather than 90%, confidence regions of 3 GW events, rather than a single event. On detection and identification of a kilonova, a limited number of NIRCam follow-up observations could constrain the properties of matter ejected by the binary and the equation of state of dense nuclear matter.

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

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

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

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

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

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

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

The relation between ices in the envelopes and disks surrounding YSOs and those in the quiescent interstellar medium is investigated. For a sample of 31 stars behind isolated dense cores, ground-based and Spitzer spectra and photometry in the 1-25 um wavelength range are combined. The baseline for the broad and overlapping ice features is modeled, using calculated spectra of giants, H2O ice and silicates. The adopted extinction curve is derived empirically. Its high resolution allows for the separation of continuum and feature extinction. The extinction between 13-25 um is ~50% relative to that at 2.2 um. The strengths of the 6.0 and 6.85 um absorption bands are in line with those of YSOs. Thus, their carriers, which, besides H2O and CH3OH, may include NH4+, HCOOH, H2CO and NH3, are readily formed in the dense core phase, before stars form. The 3.53 um C-H stretching mode of solid CH3OH was discovered. The CH3OH/H2O abundance ratios of 5-12% are larger than upper limits in the Taurus molecular cloud. The initial ice composition, before star formation occurs, therefore depends on the environment. Signs of thermal and energetic processing that were found toward some YSOs are absent in the ices toward background stars. Finally, the peak optical depth of the 9.7 um band of silicates relative to the continuum extinction at 2.2 um is significantly shallower than in the diffuse interstellar medium. This extends the results of Chiar et al. (2007) to a larger sample and higher extinctions.

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

Infrared extinction maps and submillimeter dust continuum maps are powerful probes of the density structure in the envelope of star-forming cores. We make a direct comparison between infrared and submillimeter dust continuum observations of the low-mass Class 0 core, B335, to constrain the ratio of submillimeter to infrared opacity (\kaprat) and the submillimeter opacity power-law index ($\kappa \propto \lambda^{-\beta}$). Using the average value of theoretical dust opacity models at 2.2 \micron, we constrain the dust opacity at 850 and 450 \micron . Using new dust continuum models based upon the broken power-law density structure derived from interferometric observations of B335 and the infall model derived from molecular line observations of B335, we find that the opacity ratios are $\frac{\kappa_{850}}{\kappa_{2.2}} = (3.21 - 4.80)^{+0.44}_{-0.30} \times 10^{-4}$ and $\frac{\kappa_{450}}{\kappa_{2.2}} = (12.8 - 24.8)^{+2.4}_{-1.3} \times 10^{-4}$ with a submillimeter opacity power-law index of $\beta_{smm} = (2.18 - 2.58)^{+0.30}_{-0.30}$. The range of quoted values are determined from the uncertainty in the physical model for B335. For an average 2.2 \micron\ opacity of $3800 \pm 700$ cm$^2$g$^{-1}$, we find a dust opacity at 850 and 450 \micron\ of $\kappa_{850} = (1.18 - 1.77)^{+0.36}_{-0.24}$ and $\kappa_{450} = (4.72 - 9.13)^{+1.9}_{-0.98}$ cm$^2$g$^{-1}$ of dust. These opacities are from $(0.65 - 0.97) \kappa^{\rm{OH}5}_{850}$ of the widely used theoretical opacities of Ossenkopf and Henning for coagulated ice grains with thin mantles at 850\micron.

We investigate the relation between star formation rate (SFR) and gas surface densities in Galactic star forming regions using a sample of YSOs and massive clumps. Our YSO sample consists of objects located in 20 molecular clouds from the Spitzer cores to disks and Gould's Belt surveys. We estimate the gas surface density (Sigma_gas) from Av maps and YSO SFR surface densities (Sigma_SFR) from the number of YSOs, assuming a mean mass and lifetime. We also divide the clouds into contour levels of Av, counting only the youngest Class I and Flat SED YSOs. For a sample of massive star forming clumps, we derive SFRs from the infrared luminosity and use HCN gas maps to estimate Sigma_gas. We find that Galactic clouds lie above the extragalactic relations (e.g., Kennicutt-Schmidt Law) by factors up to 17. Cloud regions with high Sigma_gas lie above extragalactic relations up to a factor of 54 and overlap with massive clumps. We use 12CO and 13CO gas maps of the Perseus and Ophiuchus clouds to estimate Sigma_gas and compare to Sigma_gas from Av maps. We find that 13CO, underestimates the Av-based mass by factors of 4-5. 12CO may underestimate the total gas mass at Sigma_gas > 200 Msun pc^-2 by > 30%;however, this does not explain the large discrepancy between Galactic and extragalactic relations. We find evidence for a threshold of star formation (Sigma_th) at 129+-14 Msun pc^-2. At Sigma_gas > Sigma_th, the Galactic relation is linear. A possible reason for the difference between Galactic and extragalactic relations is that all the CO-emitting gas, including Sigma_gas below Sigma_th, is measured in extragalactic studies. If the Kennicutt-Schmidt relation (Sigma_SFR Sigma_gas^1.4) and a linear relation between dense gas and star formation is assumed, the fraction of dense star forming gas (f_dense) increases as Sigma_gas^0.4. When Sigma_gas reaches ~300Sigma_th, f_dense is 1. (Abridged)

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.

Aims. The brown dwarf (BD) formation process has not yet been completely understood. To shed more light on the differences and similarities between star and BD formation processes, we study and compare the disk fraction among both kinds of objects over a large angular region in the Taurus cloud. In addition, we examine the spatial distribution of stars and BD relative to the underlying molecular gas Methods. In this paper, we present new and updated photometry data from the Infrared Array Camera (IRAC) aboard the Spitzer Space Telescope on 43 BDs in the Taurus cloud, and recalculate of the BD disk fraction in this region. We also useed recently available CO mm data to study the spatial distribution of stars and BDs relative to the cloud's molecular gas. Results. We find that the disk fraction among BDs in the Taurus cloud is 41 \pm 12%, a value statistically consistent with the one among TTS (58 \pm 9%). We find that BDs in transition from a state where they have a disk to a diskless state are rare, and we study one isolated example of a transitional disk with an inner radius of \approx 0.1 AU (CFHT BD Tau 12, found via its relatively small mid-IR excess compared to most members of Taurus that have disks. We find that BDs are statistically found in regions of similar molecular gas surface density to those associated with stars. Furthermore, we find that the gas column density distribution is almost identical for stellar and substellar objects with and without disks.

Using data from the SCUBA Legacy Catalogue (850 um) and Spitzer Space Telescope (3.6 - 70 um), we explore dense cores in the Ophiuchus, Taurus, Perseus, Serpens, and Orion molecular clouds. We develop a new method to discriminate submillimeter cores found by SCUBA as starless or protostellar, using point source photometry from Spitzer wide field surveys. First, we identify infrared sources with red colors associated with embedded young stellar objects (YSOs). Second, we compare the positions of these YSO-candidates to our submillimeter cores. With these identifications, we construct new, self-consistent starless and protostellar core mass functions (CMFs) for the five clouds. We find best fit slopes to the high-mass end of the CMFs of -1.26 +/- 0.20, -1.22 +/- 0.06, -0.95 +/- 0.20, and -1.67 +/- 0.72 for Ophiuchus, Taurus, Perseus, and Orion, respectively. Broadly, these slopes are each consistent with the -1.35 power-law slope of the Salpeter IMF at higher masses, but suggest some differences. We examine a variety of trends between these CMF shapes and their parent cloud properties, potentially finding a correlation between the high-mass slope and core temperature. We also find a trend between core mass and effective size, but we are very limited by sensitivity. We make similar comparisons between core mass and size with visual extinction (for A_V >= 3) and find no obvious trends. We also predict the numbers and mass distributions of cores that future surveys with SCUBA-2 may detect in each of these clouds.

We report on the properties of pre-main-sequence objects in the Taurus molecular clouds as observed in 7 mid- and far-infrared bands with the Spitzer Space Telescope. There are 215 previously-identified members of the Taurus star-forming region in our ~44 square degree map; these members exhibit a range of Spitzer colors that we take to define young stars still surrounded by circumstellar dust (noting that ~20% of the bonafide Taurus members exhibit no detectable dust excesses). We looked for new objects in the survey field with similar Spitzer properties, aided by extensive optical, X-ray, and ultraviolet imaging, and found 148 candidate new members of Taurus. We have obtained follow-up spectroscopy for about half the candidate sample, thus far confirming 34 new members, 3 probable new members, and 10 possible new members, an increase of 15-20% in Taurus members. Of the objects for which we have spectroscopy, 7 are now confirmed extragalactic objects, and one is a background Be star. The remaining 93 candidate objects await additional analysis and/or data to be confirmed or rejected as Taurus members. Most of the new members are Class II M stars and are located along the same cloud filaments as the previously-identified Taurus members. Among non-members with Spitzer colors similar to young, dusty stars are evolved Be stars, planetary nebulae, carbon stars, galaxies, and AGN.

We present Spitzer images of the Taurus Complex (TC) and take advantage of the sensitivity and spatial resolution of the observations to characterize the diffuse IR emission across the cloud. This work highlights evidence of dust evolution within the translucent sections of the archetype reference for studies of quiescent molecular clouds. We combine Spitzer 160 um and IRAS 100 um observations to produce a dust temperature map and a far-IR dust opacity map at 5' resolution. The average dust temperature is about 14.5K with a dispersion of +/-1K across the cloud. The far-IR dust opacity is a factor 2 larger than the average value for the diffuse ISM. This opacity increase and the attenuation of the radiation field (RF) both contribute to account for the lower emission temperature of the large grains. The structure of the TC significantly changes in the mid-IR images that trace emission from PAHs and VSGs. We focus our analysis of the mid-IR emission to a range of ecliptic latitudes where the zodiacal light residuals are small. Within this cloud area, there are no 8 and 24 um counterparts to the brightest 160 um emission features. Conversely, the 8 and 24 um images reveal filamentary structure that is strikingly inconspicuous in the 160 um and extinction maps. The IR colors vary over sub-parsec distances across this filamentary structure. We compare the observed colors with model calculations quantifying the impact of the RF intensity and the abundance of stochastically heated particles on the dust SED. To match the range of observed colors, we have to invoke variations by a factor of a few of both the interstellar RF and the abundance of PAHs and VSGs. We conclude that within this filamentary structure a significant fraction of the dust mass cycles in and out the small size end of the dust size distribution.

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

We combine the results from several multiplicity surveys of pre-main-sequence stars located in four nearby star-forming regions with Spitzer data from three different Legacy Projects. This allows us to construct a sample of 349 targets, including 125 binaries, which we use to to investigate the effect of companions on the evolution of circumstellar disks. We find that the distribution of projected separations of systems with Spitzer excesses is significantly different (P ~2.4e-5, according to the KS test for binaries with separations < 400 AU) from that of systems lacking evidence for a disk. As expected, systems with projected separations < 40 AU are half as likely to retain at least one disk than are systems with projected separations in the 40-400 AU range. These results represent the first statistically significant evidence for a correlation between binary separation and the presence of an inner disk (r ~ 1 AU). Several factors (e.g., the incompleteness of the census of close binaries, the use of unresolved disk indicators, and projection effects) have previously masked this correlation in smaller samples. We discuss the implications of our findings for circumstellar disk lifetimes and the formation of planets in multiple systems.

We investigate the environment of the very low-luminosity object L1521F IRS using data from the Taurus Spitzer Legacy Survey. The MIPS 160 micron image shows both extended emission from the Taurus cloud as well as emission from multiple cold cores over a 1 X 2 deg region. Analysis shows that the cloud dust temperature is 14.2 +- 0.4 K and the extinction ratio is A_160/A_K = 0.010 +- 0.001 up to A_V ~ 4 mag. We find kappa_160 = 0.23 +- 0.046 cm^2/g for the specific opacity of the gas-dust mixture. Therefore, for dust in the Taurus cloud we find the 160 um opacity is significantly higher than that measured for the diffuse ISM, but not too different from dense cores, even at modest extinction values. Furthermore, the 160 um image shows features that do not appear in the IRAS 100 um image. We identify six regions as cold cores, i.e. colder than 14.2 K, all of which have counterparts in extinction maps or C18O maps. We compare the effects of L1521F IRS on its natal core and find there is no evidence for dust heating at 160 or 100 um by the embedded source. From the infrared luminosity L_TIR = 0.024 Lo we find L_bol = 0.034 - 0.046L_o, thus confirming the source's low-luminosity. Comparison of L1521F IRS with theoretical simulations for the very early phases of star formation appears to rule out the first core collapse phase. The evolutionary state appears similar to or younger than the class 0 phase, and the estimated mass is likely to be substellar.

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

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

Gravitational wave (GW) bursts (short duration signals) are expected to be associated with highly energetic astrophysical processes. With such high energies present, it is likely these astrophysical events will have signatures in the EM spectrum as well as in gravitational radiation. We have initiated a program, "Locating and Observing Optical Counterparts to Unmodeled Pulses in Gravitational Waves" (LOOC UP) to promptly search for counterparts to GW burst candidates. The proposed method analyzes near real-time data from the LIGO-Virgo network, and then uses a telescope network to seek optical-transient counterparts to candidate GW signals. We carried out a pilot study using S5/VSR1 data from the LIGO-Virgo network to develop methods and software tools for such a search. We will present the method, with an emphasis on the potential for such a search to be carried out during the next science run of LIGO and Virgo, expected to begin in 2009.

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 c2d Spitzer/IRAC observations of the Lupus I, III and IV dark clouds and discuss them in combination with optical and near-infrared and c2d MIPS data. With the Spitzer data, the new sample contains 159 stars, 4 times larger than the previous one. It is dominated by low- and very-low mass stars and it is complete down to M $\approx$ 0.1M$_\odot$. We find 30-40 % binaries with separations between 100 to 2000 AU with no apparent effect in the disk properties of the members. A large majority of the objects are Class II or Class III objects, with only 20 (12%) of Class I or Flat spectrum sources. The disk sample is complete down to ``debris''-like systems in stars as small as M $\approx$ 0.2 M$_\odot$ and includes sub-stellar objects with larger IR excesses. The disk fraction in Lupus is 70 -- 80%, consistent with an age of 1 -- 2 Myr. However, the young population contains 20% optically thick accretion disks and 40% relatively less flared disks. A growing variety of inner disk structures is found for larger inner disk clearings for equal disk masses. Lupus III is the most centrally populated and rich, followed by Lupus I with a filamentary structure and by Lupus IV, where a very high density core with little star-formation activity has been found. We estimate star formation rates in Lupus of 2 -- 10 M$_\odot$ Myr$^{-1}$ and star formation efficiencies of a few percent, apparently correlated with the associated cloud masses.

We discuss the results of the optical spectroscopic follow-up of pre-main sequence (PMS) objects and candidates selected in the Chamaeleon II dark cloud based on data from the Spitzer Legacy survey "From Molecular Cores to Planet Forming Disks" (c2d) and from previous surveys. Our sample includes both objects with infrared excess selected according to c2d criteria and referred to as Young Stellar Objects and other cloud members and candidates selected from complementary optical and near-infrared data. We characterize the sample of objects by deriving their physical parameters. The vast majority of objects have masses < 1 solar mass and ages < 6 Myr. Several of the PMS objects and candidates lie very close to or below the Hydrogen-burning limit. A first estimate of the slope of the Initial Mass Function in Cha II is consistent with that of other T associations. The star formation efficiency in the cloud (1-4%) is consistent with our own estimates for Taurus and Lupus, but significantly lower than for Cha I. This might mean that different star-formation activities in the Chamaeleon clouds may reflect a different history of star formation. We also find that the Cha II cloud is turning some 8 solar masses into stars every Myr, which is less than the star formation rate in the other c2d clouds. However, the star formation rate is not steady and evidence is found that the star formation in Cha II might have occurred very rapidly. The H_alpha emission of the Cha II PMS objects, as well as possible correlations between their stellar and disk properties, are also investigated.

We discuss the results from the combined IRAC and MIPS c2d Spitzer Legacy survey observations and complementary optical and near infrared data of the Chamaeleon II (Cha II) dark cloud. We perform a census of the young population of Cha II, in a mapped area of ~1.75 square degrees, and study the spatial distribution and properties of the cloud members and candidate pre-main sequence (PMS) objects and their circumstellar matter. From the analysis of the volume density of the PMS objects and candidates we find two tight groups of objects with volume densities higher than 25 solar masses per cubic parsec and 5-10 members each. These groups correlate well in space with the regions of high extinction. A multiplicity fraction of about 13% is observed for objects with separations between 0.8" and 6.0". Using the results of masses and ages from a companion paper, we estimate the star formation efficiency to be 1-4% significantly lower than for Cha I. This might mean that different star-formation activities in the Chamaeleon clouds reflect a different history of star formation. We also find that the Cha II cloud is turning some 6-7 solar masses into stars every Myr, which is low in comparison with the star formation rate in other c2d clouds. On the other hand, the disk fraction of 70-80% that we estimate in Cha II is much higher than in other star forming regions and indicates that the population in this cloud is dominated by objects with active accretion. Finally, the Cha II outflows are discussed, with particular regard to the discovery of a new Herbig-Haro outflow, HH 939, driven by the classical T Tauri star Sz 50.

With the goal to study the physical and chemical evolution of ices in solar-mass systems, a spectral survey is conducted of a sample of 41 low luminosity YSOs using 3-38 um Spitzer and ground-based spectra. The long-known 6.0 and 6.85 um bands are detected toward all sources, with the Class 0-type YSOs showing the deepest bands ever observed. In almost all sources the 6.0 um band is deeper than expected from the bending mode of pure solid H2O. The depth and shape variations of the remaining 5-7 um absorption indicate that it consists of 5 independent components, which, by comparison to laboratory studies, must be from at least 8 different carriers. Simple species are responsible for much of the absorption in the 5-7 um region, at abundances of 1-30% for CH3OH, 3-8% for NH3, 1-5% for HCOOH, ~6% for H2CO, and ~0.3% for HCOO- with respect to solid H2O. The 6.85 um band likely consists of one or two carriers, of which one is less volatile than H2O because its abundance relative to H2O is enhanced at lower H2O/tau_9.7 ratios. It does not survive in the diffuse interstellar medium (ISM), however. The similarity of the 6.85 um bands for YSOs and background stars indicates that its carrier(s) must be formed early in the molecular cloud evolution. If an NH4+ salt is the carrier its abundance with respect to solid H2O is typically 7%, and low temperature acid-base chemistry or cosmic ray induced reactions must have been involved in its formation. Possible origins are discussed for the carrier of an enigmatic, very broad absorption between 5 and 8 um. Finally, all the phenomena observed for ices toward massive YSOs are also observed toward low mass YSOs, indicating that processing of the ices by internal ultraviolet radiation fields is a minor factor in the early chemical evolution of the ices. [abridged]

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 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 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 undertaken an unprecedentedly large 1.1 millimeter continuum survey of three nearby star forming clouds using Bolocam at the Caltech Submillimeter Observatory. We mapped the largest areas in each cloud at millimeter or submillimeter wavelengths to date: 7.5 sq. deg in Perseus (Paper I), 10.8 sq. deg in Ophiuchus (Paper II), and 1.5 sq. deg in Serpens with a resolution of 31", detecting 122, 44, and 35 cores, respectively. Here we report on results of the Serpens survey and compare the three clouds. Average measured angular core sizes and their dependence on resolution suggest that many of the observed sources are consistent with power-law density profiles. Tests of the effects of cloud distance reveal that linear resolution strongly affects measured source sizes and densities, but not the shape of the mass distribution. Core mass distribution slopes in Perseus and Ophiuchus (alpha=2.1+/-0.1 and alpha=2.1+/-0.3) are consistent with recent measurements of the stellar IMF, whereas the Serpens distribution is flatter (alpha=1.6+/-0.2). We also compare the relative mass distribution shapes to predictions from turbulent fragmentation simulations. Dense cores constitute less than 10% of the total cloud mass in all three clouds, consistent with other measurements of low star-formation efficiencies. Furthermore, most cores are found at high column densities; more than 75% of 1.1 mm cores are associated with Av>8 mag in Perseus, 15 mag in Serpens, and 20-23 mag in Ophiuchus.

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.

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

Aims : We want to investigate whether brown dwarfs (BDs) form like stars or are ejected embryos. We study the presence of disks around BDs in the Taurus cloud, and discuss implications for substellar formation models. Methods : We use photometric measurements from the visible to the far infrared to determine the spectral energy distributions (SEDs) of Taurus BDs. Results: We use Spitzer color indices, Halpha as an accretion indicator, and models fit to the SEDs in order to estimate physical parameters of the disks around these BDs. We study the spatial distribution of BDs with and without disks across the Taurus aggregates, and we find that BDs with and without disks are not distributed regularly across the Taurus cloud. Conclusions: We find that 48%+/- 14% of Taurus BDs have a circumstellar disk signature, a ratio similar to recent results from previous authors in other regions. We fit the SEDs and find that none of the disks around BDs in Taurus can be fitted convincingly with a flaring index beta = 0, indicating that heating by the central object is efficient and that the disks we observe retain a significant amount of gas. We find that BDs with disks are proportionally more numerous in the northern Taurus filament, possibly the youngest filament. We do not find such a clear segregation for classical T Tauri stars (CTTS) and weak-lined T Tauri stars (WTTS), suggesting that, in addition to the effects of evolution, any segregation effects could be related to the mass of the object. A by-product of our study is to propose a recalibration of the Barrado y Navascues & Martin (2003) accretion limit in the substellar domain. The global shape of the limit fits our data points if it is raised by a factor 1.25-1.30.

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

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

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

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

We present a large-scale millimeter continuum map of the Ophiuchus molecular cloud. Nearly 11 square degrees, including all of the area in the cloud with visual extinction more than 3 magnitudes, was mapped at 1.1 mm with Bolocam on the Caltech Submillimeter Observatory (CSO). By design, the map also covers the region mapped in the infrared with the Spitzer Space Telescope. We detect 44 definite sources, and a few likely sources are also seen along a filament in the eastern streamer. The map indicates that dense cores in Ophiuchus are very clustered and often found in filaments within the cloud. Most sources are round, as measured at the half power point, but elongated when measured at lower contour levels, suggesting spherical sources lying within filaments. The masses, for an assumed dust temperature of 10 K, range from 0.24 to 3.9 solar masses, with a mean value of 0.96 solar masses. The total mass in distinct cores is 42 solar masses, 0.5 to 2% of the total cloud mass, and the total mass above 4 sigma is about 80 solar masses. The mean densities in the cores are quite high, with an average of 1.6 x 10^6 per cc, suggesting short free-fall times. The core mass distribution can be fitted with a power law with slope of 2.1 plus or minus 0.3 for M>0.5 solar masses, similar to that found in other regions, but slightly shallower than that of some determinations of the local IMF. In agreement with previous studies, our survey shows that dense cores account for a very small fraction of the cloud volume and total mass. They are nearly all confined to regions with visual extinction at least 9 mag, a lower threshold than found previously.

We have completed a 1.1 mm continuum survey of 7.5 sq deg of the Perseus Molecular Cloud using Bolocam at the Caltech Submillimeter Observatory. This represents the largest millimeter or submillimeter continuum map of Perseus to date. Our map covers more than 30,000 31" (FWHM) resolution elements to a 1 sigma RMS of 15 mJy/beam. We detect a total of 122 cores above a 5 sigma point source mass detection limit of 0.18 M_sun, assuming a dust temperature of 10 K, 60 of which are new millimeter or submillimeter detections. The 1.1 mm mass function is consistent with a broken power law of slope -1.3 (0.5 M_sun<M<2.5 M_sun) and -2.6 (M>2.5 M_sun), similar to the local initial mass function slope. No more than 5% of the total cloud mass is contained in discrete 1.1 mm cores, which account for a total mass of 285 M_sun. We suggest an extinction threshold for millimeter cores of Av~5 mag, based on our calculation of the probability of finding a 1.1 mm core as a function of Av. Much of the cloud is devoid of compact millimeter emission; despite the significantly greater area covered compared to previous surveys, only 5-10 of the newly identified sources lie outside previously observed areas. The two-point correlation function confirms that dense cores in the cloud are highly structured, with significant clustering on scales as large as 2e5 AU. These 1.1 mm results, especially when combined with recently acquired c2d Spitzer Legacy data, will provide a census of dense cores and protostars in Perseus and improve our understanding of the earliest stages of star formation in molecular clouds.

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.

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

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.

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.

We present a new infrared extinction study of Globule 2, the most opaque molecular cloud core in the Coalsack complex. Using deep near-infrared imaging observations obtained with the ESO NTT we are able to examine the structure of the globule in significantly greater detail than previously possible. We find the most prominent structural feature of this globule to be a strong central ring of dust column density which was not evident in lower resolution studies of this cloud. This ring represents a region of high density and pressure that is likely a transient structure. For a spherical cloud geometry the ring would correspond to a dense inner shell of high pressure that could not be in dynamical equilibrium with its surroundings since there appear to be no sources of pressure in the central regions of the cloud that could support the shell against gravity and prevent its inward implosion. The timescale for the inward collapse of the ring would be less than 2 x 10^5 years, suggesting that this globule is in an extremely early stage of evolution, and perhaps caught in the process of forming a centrally condensed dense core or Bok globule. Outside its central regions the globule displays a well-behaved density profile whose shape is very similar to that of a stable Bonnor-Ebert sphere. Using SEST we also obtained a C18O spectrum toward the center of the cloud. The CO observation indicates that the globule is a gravitationally bound object. Analysis of the CO line profile reveals significant non-thermal gas motions likely due to turbulence. As a whole the globule may be evolving to a global state of quasi-static dynamical equilibrium in which thermal and turbulent pressure balance gravity.

We present wide-field near-infrared (JHK) images of the young, 2 Myr IC 348 cluster taken with FLAMINGOS. We use these new data to construct an infrared census of sources, which is sensitive enough to detect a 10 Mjup brown dwarf seen through an extinction of Av=7mag. We examine the cluster's structure and relationship to the molecular cloud and construct the cluster's K band luminosity function. Using our model luminosity function algorithm, we derive the cluster's initial mass function throughout the stellar and substellar regimes and find that the IC 348 IMF is very similar to that found for the Trapezium Cluster with both cluster IMFs having a mode between 0.2 - 0.08 Msun. In particular we find that, similar to our results for the Trapezium, brown dwarfs constitute only 1 in 4 of the sources in the IC 348 cluster. We show that a modest secondary peak forms in the substellar IC 348 KLF, corresponding to the same mass range responsible for a similar KLF peak found in the Trapezium. We interpret this KLF peak as either evidence for a corresponding secondary IMF peak at the deuterium burning limit, or as arising from a feature in the substellar mass-luminosity relation that is not predicted by current theoretical models. Lastly, we find that IC 348 displays radial variations of its sub-solar (0.5 - 0.08 Msun) IMF on a parsec scale. Whatever mechanism that is breaking the universality of the IMF on small spatial scales in IC 348 does not appear to be acting upon the brown dwarf population, whose relative size does not vary with distance from the cluster center.

We report sensitive, high resolution molecular-line observations of the dark cloud Barnard 68 obtained with the IRAM 30-m telescope. We analyze spectral-line observations of C18O, CS(2--1), C34S(2--1), and N2H+(1--0) in order to investigate the kinematics and dynamical state of the cloud. We find extremely narrow linewidths in the central regions of the cloud. These narrow lines are consistent with thermally broadened profiles for the measured gas temperature of 10.5 K. We determine the thermal pressure to be a factor 4 -- 5 times greater than the non-thermal (turbulent) pressure in the central regions of the cloud, indicating that thermal pressure is the primary source of support against gravity in this cloud. This confirms the inference of a thermally supported cloud drawn previously from deep infrared extinction measurements. The rotational kinetic energy is found to be only a few percent of the gravitational potential energy, indicating that the contribution of rotation to the overall stability of the cloud is insignificant. Finally, our observations show that CS line is optically thick and self-reversed across nearly the entire projected surface of the cloud. The shapes of the self-reversed profiles are asymmetric and are found to vary across the cloud in such a manner that the presence of both inward and outward motions are observed within the cloud. Moreover, these motions appear to be globally organized in a clear and systematic alternating spatial pattern which is suggestive of a small amplitude, non-radial oscillation or pulsation of the outer layers of the cloud about an equilibrium configuration.

We present sensitive, high angular resolution molecular-line observations of C18O and N2H+ toward the dark globule B68. We directly compare these data with the near-infrared extinction measurements of Alves, Lada, & Lada (2001) to derive the first evidence for the depletion of N2H+, and by inference N2, in a pre-stellar dark cloud. We also find widespread C18O depletion throughout the centrally condensed core of the B68 cloud. Specifically, we find the N2H+ emission to peak in a shell partially surrounding the peak of dust extinction. Moreover, the N2H+ peaks inside the much larger C18O depletion hole and has a smaller depletion zone, confirming theoretical predictions. These data are analyzed through a direct coupling of time dependent chemical models to a radiation transfer code. This analysis highlights the importance of photodissociation at cloud edges and suggests that the CO abundance declines by two orders of magnitude from edge to center. In contrast N2H+ declines in abundance, at minimum, by at least a factor of two. Indeed it is entirely possible that both N2H+ and N2 are completely absent from the central regions of the B68 core. The depletion of N2H+, and its parent molecule N2, opens the possibility that the centers of dense cores, prior to the formation of a star, may evade detection by conventional methods of probing cores using molecular emission. Under these conditions H2D+ may be the sole viable molecular probe of the innermost regions of star forming cores.