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

We report ALMA and VLA continuum observations that potentially identify the four youngest protostars in the Orion Molecular Clouds taken as part of the Orion VANDAM program. These are distinguished by bright, extended, irregular emission at 0.87 mm and 8 mm and are optically thick at 0.87 mm. These structures are distinct from the disk or point-like morphologies seen toward the other Orion protostars. The 0.87 mm emission implies temperatures of 41-170 K, requiring internal heating. The bright 8 mm emission implies masses of 0.5 to 1.2 M_sun assuming standard dust opacity models. One source has a Class 0 companion, while another exhibits substructure indicating a companion-candidate. Three compact outflows are detected, two of which may be driven by companions, with dynamical times of ~300 to ~400 years. The slowest outflow may be driven by a first hydrostatic core. These protostars appear to trace an early phase when the centers of collapsing fragments become optically thick to their own radiation and compression raises the gas temperature. This phase is thought to accompany the formation of hydrostatic cores. A key question is whether these structures are evolving on free fall times of ~100 years, or whether they are evolving on Kelvin-Helmholtz times of several thousand years. The number of these sources imply a lifetime of ~6000 years, in closer agreement with the Kelvin-Helmholtz time. In this case, rotational and/or magnetic support could be slowing the collapse.

We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase.

We present ALMA (0.87~mm) and VLA (9~mm) observations toward OMC2-FIR4 and OMC2-FIR3 within the Orion integral-shaped filament that are thought to be the nearest regions of intermediate mass star formation. We characterize the continuum sources within these regions on $\sim$40~AU (0\farcs1) scales and associated molecular line emission at a factor of $\sim$30 better resolution than previous observations at similar wavelengths. We identify six compact continuum sources within OMC2-FIR4, four in OMC2-FIR3, and one additional source just outside OMC2-FIR4. This continuum emission is tracing the inner envelope and/or disk emission on less than 100~AU scales. HOPS-108 is the only protostar in OMC2-FIR4 that exhibits emission from high-excitation transitions of complex organic molecules (e.g., methanol and other lines) coincident with the continuum emission. HOPS-370 in OMC2-FIR3 with L~$\sim$~360~\lsun, also exhibits emission from high-excitation methanol and other lines. The methanol emission toward these two protostars is indicative of temperatures high enough to thermally evaporate methanol from icy dust grains; overall these protostars have characteristics similar to hot corinos. We do not identify a clear outflow from HOPS-108 in \twco, but find evidence of interaction between the outflow/jet from HOPS-370 and the OMC2-FIR4 region. The multitude of observational constraints indicate that HOPS-108 is likely a low to intermediate-mass protostar in its main mass accretion phase and it is the most luminous protostar in OMC2-FIR4. The high resolution data presented here are essential for disentangling the embedded protostars from their surrounding dusty environments and characterizing them.

Significant advances have been made over the past decade in the characterization of multiple protostar systems, enabled by the Karl G. Jansky Very Large Array (VLA), high-resolution infrared observations with the Hubble Space Telescope, and ground-based facilities. To further understand the mechanism(s) of multiple star formation, a combination of statistics, high-angular resolution radio/millimeter continuum imaging, characterization of kinematic structure, magnetic fields via polarimetry, and comparison with numerical simulations are needed. Thus, understanding the origin of stellar multiplicity in different regimes of companion separation will soon be within reach. However, to overcome challenges that studies in this field are now confronted with, a range of new capabilities are required: a new millimeter/centimeter wave facility with 10 mas resolution at {\lambda}=1 cm, space-based near to far-infrared observatories, continued development of low to high resolution spectroscopy on 3m to 10m class telescopes, and an ELT-class telescope with near to mid-infrared imaging/spectroscopic capability.

Knowledge of protostellar evolution has been revolutionized with the advent of surveys at near-infrared to submillimeter wavelengths. This has enabled the bolometric luminosities and bolometric temperatures (traditional protostellar evolution diagnostics) to be measured for large numbers of protostars. However, further progress is difficult without knowing the masses of the central protostars. Protostar masses can be most accurately determined via molecular line kinematics from millimeter interferometers (i.e., ALMA). Theoretical investigations have predicted the protostellar mass function (PMF) for various protostellar mass accretion models, and it is now imperative to observationally constrain its functional form. While ALMA has enabled protostellar mass measurements, samples approaching 100 sources are necessary to constrain the functional form of the PMF, and upgrades to ALMA and/or a new mm/cm facility will increase the feasibility of measuring such a large number of protostar masses. The masses of protostars will enable their stellar structure (radius and intrinsic luminosity), evolution, and accretion histories to be better understood. This is made more robust when effective temperatures and accretion rates can be measured via ground/space-based near to mid-infrared spectroscopy. Furthermore, access to supercomputing facilities is essential to fit the protostar masses via radiative transfer modeling and updated theoretical/numerical modeling of stellar structure may also be required.

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

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

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

We present a study of the kinematics and structure of the Cep OB3b cluster based on new spectra obtained with the Hectoschelle spectrograph on the MMT and data from Spitzer, Chandr}, and Gaia. At a distance of 819+/-16 pc, Cep OB3b is one of the closest examples of a young (~3 - 5 Myr), large (~3000 total members) cluster at the late stages of gas dispersal. The cluster is broken into two sub-clusters surrounded by a lower density halo. We fit the empirical density law of King (1962) to each sub-cluster to constrain their sizes and structure. The richer eastern sub-cluster has circular symmetry, a modest central density, and lacks molecular gas toward its core suggesting it has undergone expansion due to gas dispersal. In contrast, the western sub-cluster deviates from circular symmetry, has a smaller core size, and contains significant molecular gas near its core, suggesting that it is in an earlier phase of gas dispersal. We present posterior probability distributions for the velocity dispersions from the Hectoschelle spectra. The east will continue to expand and likely form a bound cluster with ~35% of stars remaining. The west is undergoing slower gas dispersal and will potentially form a bound cluster with ~75% of stars remaining. If the halo dissipates, this will leave two independent clusters with ~300 members; proper motions suggest that the two sub-clusters are not bound to each other.

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

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

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

We present the results of a 40 ks XMM-Newton observation centered on the variable star V1818 Ori. Using a combination of the XMM-Newton and AllWISE catalog data, we identify a group of about 31 young stellar objects around V1818 Ori. This group is coincident with the eastern edge of the dust ring surrounding Kappa Ori. Previously, we concluded that the young stellar objects on the western side of ring were formed in an episode of star formation that started 3-5 Myr ago, and are at a distance similar to that of kappa Ori (250-280 pc) and in the foreground to the Orion A cloud. Here we use the XMM-Newton observation to calculate X-ray fluxes and luminosities of the young stars around V1818 Ori. We find that their X-ray luminosity function (XLF), calculated for a distance of ~270 pc, matches the XLF of the YSOs west of Kappa Ori. We rule out that this group of young stars is associated to Mon R2 as assumed in the literature, but rather they are part of the same Kappa Ori's ring stellar population.

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

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

The progenitors of high-mass stars and clusters are still challenging to recognise. Only unbiased surveys, sensitive to compact regions of high dust column density, can unambiguously reveal such a small population of particularly massive and cold clumps. Here we study a flux limited sample of compact sources from the ATLASGAL survey to identify a sample of candidate progenitors of massive clusters in the inner Galaxy. Sensitive mid-infrared data at 21-24 $\mu$m from the WISE and MIPSGAL surveys were explored to search for embedded objects, and complementary spectroscopic data were used to investigate their stability and star formation activity. Based on such ancillary data we identify an unbiased sample of infrared-quiet massive clumps in the Galaxy that potentially represent the earliest stages of massive cluster formation. An important fraction of this sample consists of sources that have not been studied in detail before. Comparing their properties to clumps hosting more evolved embedded objects, we find that they exhibit similar physical properties in terms of mass and size, suggesting that infrared-quiet massive clumps are not only capable of forming high-mass stars, but likely also follow a single evolutionary track leading to the formation of massive clusters. The majority of the sources are not in virial-equilibrium, suggesting collapse on the clump scale. This is in line with the low number of infrared-quiet massive clumps and earlier findings that star formation, in particular for high-mass objects is a fast, dynamic process. We propose a scenario in which massive clumps start to fragment and collapse before their final mass is accumulated indicating that strong self-gravity and global collapse is needed to build up rich clusters and the most massive stars.

Infrared variability is common among young stellar objects, with surveys finding daily to weekly fluctuations of a few tenths of a magnitude. Space-based observations can produce highly sampled infrared light curves, but are often limited to total baselines of about a month due to the orientation of the spacecraft. Here we present observations of the Chameleon I cluster whose low declination makes it observable by the Spitzer space telescope over a 200 day period. We observe 30 young stellar objects with a daily cadence to better sample variability on timescales of months. We find such variability is common, occurring in ~80% of the detected cluster members. The change in [3.6]-[4.5] color over 200 days for many of the sources falls between that expected for extinction and fluctuations in disk emission. With our high cadence and long baseline we can derive power spectral density curves covering two orders of magnitude in frequency and find significant power at low frequencies, up to the boundaries of our 200 day survey. Such long timescales are difficult to explain with variations driven by the interaction between the disk and stellar magnetic field, which has a dynamical timescale of days to weeks. The most likely explanation is either structural or temperature fluctuations spread throughout the inner ~0.5 au of the disk, suggesting that the intrinsic dust structure is highly dynamic.

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

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

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

We present our investigation of 319 Class II objects in Orion A observed by $Spitzer$/IRS. We also present the follow-up observation of 120 of these Class II objects in Orion A from IRTF/SpeX. We measure continuum spectral indices, equivalent widths, and integrated fluxes that pertain to disk structure and dust composition from IRS spectra of Class II objects in Orion A. We estimate mass accretion rates using hydrogen recombination lines in the SpeX spectra of our targets. Utilizing these properties, we compare the distributions of the disk and dust properties of Orion A disks to those of Taurus disks with respect to position within Orion A (ONC and L1641) and to the sub-groups by the inferred radial structures, such as transitional disks vs. radially continuous full disks. Our main findings are as follows. (1) Inner disks evolve faster than the outer disks. (2) Mass accretion rate of transitional disks and that of radially continuous full disks are statistically significantly displaced from each other. The median mass accretion rate of radially continuous disks in ONC and L1641 is not very different from that in Taurus. (3) Less grain processing has occurred in the disks in ONC compared to those in Taurus, based on analysis of the shape index of the 10 $\mu$m silicate feature ($F_{11.3}/F_{9.8}$). (4) The 20-31 $\mu$m continuum spectral index tracks the projected distance from the most luminous Trapezium star, $\theta^{1}$ Ori C. A possible explanation is the effect of UV ablation of the outer part of the disks.

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

X-rays are a powerful probe of activity in early stages of star formation. They allow us to identify young stars even after they have lost the IR signatures of circumstellar disks and provide constraints on their distance. Here we report on XMM-Newton observations which detect 121 young stellar objects (YSOs) in two fields between L1641S and $\kappa$ Ori. These observations extend the Survey of Orion A with XMM and Spitzer (SOXS). The YSOs are contained in a ring of gas and dust apparent at millimeter wavelengths, and in far-IR and near-IR surveys. The X-ray luminosity function of the young stellar objects detected in the two fields indicates a distance of 250-280 pc, much closer than the Orion A cloud and similar to distance estimates of $\kappa$ Ori. We propose that the ring is a 5-8 pc diameter shell that has been swept up by $\kappa$ Ori. This ring contains several groups of stars detected by Spitzer and WISE including one surrounding the Herbig Ae/Be stars V1818 Ori. In this interpretation, the $\kappa$ Ori ring is one of several shells swept up by massive stars within the Orion Eridanus Superbubble, and is unrelated to the southern portion of Orion A / L1641 S.

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

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

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.

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

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

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

The YSOVAR (Young Stellar Object VARiability) Spitzer Space Telescope observing program obtained the first extensive mid-infrared (3.6 & 4.5 um) time-series photometry of the Orion Nebula Cluster plus smaller footprints in eleven other star-forming cores (AFGL490, NGC1333, MonR2, GGD 12-15, NGC2264, L1688, Serpens Main, Serpens South, IRAS 20050+2720, IC1396A, and Ceph C). There are ~29,000 unique objects with light curves in either or both IRAC channels in the YSOVAR data set. We present the data collection and reduction for the Spitzer and ancillary data, and define the "standard sample" on which we calculate statistics, consisting of fast cadence data, with epochs about twice per day for ~40d. We also define a "standard sample of members", consisting of all the IR-selected members and X-ray selected members. We characterize the standard sample in terms of other properties, such as spectral energy distribution shape. We use three mechanisms to identify variables in the fast cadence data--the Stetson index, a chi^2 fit to a flat light curve, and significant periodicity. We also identified variables on the longest timescales possible of ~6 years, by comparing measurements taken early in the Spitzer mission with the mean from our YSOVAR campaign. The fraction of members in each cluster that are variable on these longest timescales is a function of the ratio of Class I/total members in each cluster, such that clusters with a higher fraction of Class I objects also have a higher fraction of long-term variables. For objects with a YSOVAR-determined period and a [3.6]-[8] color, we find that a star with a longer period is more likely than those with shorter periods to have an IR excess. We do not find any evidence for variability that causes [3.6]-[4.5] excesses to appear or vanish within our data; out of members and field objects combined, at most 0.02% may have transient IR excesses.

Variability in the infrared emission from disks around pre-main sequence stars over the course of days to weeks appears to be common, but the physical cause of the changes in disk structure are not constrained. Here we present coordinated monitoring of one young cluster with the Spitzer and Chandra space telescopes aimed at studying the physical source of the variability. In fall 2011 we obtained ten epochs of Chandra ACIS photometry over a period of 30 days with a roughly 3 day cadence contemporaneous with 20 epochs of Spitzer [3.6],[4.5] photometry over 40 days with a roughly 2 day cadence of the IC 348 cluster. This cadence allows us to search for week to month long responses of the infrared emission to changes in the high-energy flux. We find no strong evidence for a direct link between the X-ray and infrared variability on these timescales among 39 cluster members with circumstellar disks. There is no significant correlation between the shape of the infrared and X-ray light curves, or between the size of the X-ray and infrared variability. Among the stars with an X-ray flare none showed evidence of a correlated change in the infrared photometry on timescales of days to weeks following the flare. This lack of connection implies that X-ray heating of the planet forming region of the disk is not significant, although we cannot rule out rapid or instantaneous changes in infrared emission.

The Cygnus-X star-forming complex is one of the most active regions of low and high mass star formation within 2 kpc of the Sun. Using mid-infrared photometry from the IRAC and MIPS Spitzer Cygnus-X Legacy Survey, we have identified over 1800 protostar candidates. We compare the protostellar luminosity functions of two regions within Cygnus-X: CygX-South and CygX-North. These two clouds show distinctly different morphologies suggestive of dissimilar star-forming environments. We find the luminosity functions of these two regions are statistically different. Furthermore, we compare the luminosity functions of protostars found in regions of high and low stellar density within Cygnus-X and find that the luminosity function in regions of high stellar density is biased to higher luminosities. In total, these observations provide further evidence that the luminosities of protostars depend on their natal environment. We discuss the implications this dependence has for the star formation process.

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

We determine the extinction law through Cep OB3b, a young cluster of 3000 stars undergoing gas dispersal. The extinction is measured toward 76 background K giants identified with MMT/Hectospec spectra. Color excess ratios were determined toward each of the giants using $V$ and $R$ photometry from the literature, $g$,$r$,$i$ and $z$ photometry from SDSS and $J$, $H$, and $K_{s}$ photometry from 2MASS. These color excess ratios were the used to construct the extinction law through the dusty material associated with Cep OB3b. The extinction law through Cep OB3b is intermediate between the $R_{V} = 3.1$ and $R_{V} = 5$ laws commonly used for the diffuse atomic ISM and dense molecular clouds, respectively. The dependence of the extinction law on line-of-sight $A_{V}$ is investigated and we find the extinction law becomes shallower for regions with $A_{V} > 2.5$ magnitudes. We speculate that the intermediate dust law results from dust processing during the dispersal of the molecular cloud by the cluster.

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

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

[Context] Two competing models describe the formation of massive stars in objects like the Orion Trapezium. In the turbulent core accretion model, the resulting stellar masses are directly related to the mass distribution of the cloud condensations. In the competitive accretion model, the gravitational potential of the protocluster captures gas from the surrounding cloud for which the individual cluster members compete. [Aims] With high resolution submillimeter observations of the structure, kinematics, and chemistry of the proto-Trapezium cluster W3 IRS5, we aim to determine which mode of star formation dominates. [Methods] We present 354 GHz Submillimeter Array observations at resolutions of 1"-3" (1800-5400 AU) of W3 IRS5. ...... [Results] The observations show five emission peaks (SMM1-5). SMM1 and SMM2 contain massive embedded stars (~20 Msun); SMM3-5 are starless or contain low-mass stars (<8 Msun). The inferred densities are high, >= 10^7 cm^-3, but the core masses are small, 0.2-0.6 Msun. The detected molecular emission reveals four different chemical zones. ...... [Conclusions] The proto-Trapezium cluster W3 IRS5 is an ideal test case to discriminate between models of massive star formation. Either the massive stars accrete locally from their local cores; in this case the small core masses imply that W3 IRS5 is at the very end stages (1000 yr) of infall and accretion, or the stars are accreting from the global collapse of a massive, cluster forming core. We find that the observed masses, densities and line widths observed toward W3 IRS 5 and the surrounding cluster forming core are consistent with the competitive accretion of gas at rates of Macc~10^-4 Msun yr^-1 by the massive young forming stars. ......

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.

We report the detection of a unique CO2 ice band toward the deeply embedded, low-mass protostar HOPS-68. Our spectrum, obtained with the Infrared Spectrograph onboard the Spitzer Space Telescope, reveals a 15.2 micron CO2 ice bending mode profile that cannot modeled with the same ice structure typically found toward other protostars. We develop a modified CO2 ice profile decomposition, including the addition of new high-quality laboratory spectra of pure, crystalline CO2 ice. Using this model, we find that 87-92% of the CO2 is sequestered as spherical, CO2-rich mantles, while typical interstellar ices show evidence of irregularly-shaped, hydrogen-rich mantles. We propose that (1) the nearly complete absence of unprocessed ices along the line-of-sight is due to the flattened envelope structure of HOPS-68, which lacks cold absorbing material in its outer envelope, and possesses an extreme concentration of material within its inner (10 AU) envelope region and (2) an energetic event led to the evaporation of inner envelope ices, followed by cooling and re-condensation, explaining the sequestration of spherical, CO2 ice mantles in a hydrogen-poor mixture. The mechanism responsible for the sublimation could be either a transient accretion event or shocks in the interaction region between the protostellar outflow and envelope. The proposed scenario is consistent with the rarity of the observed CO2 ice profile, the formation of nearly pure CO2 ice, and the production of spherical ice mantles. HOPS-68 may therefore provide a unique window into the protostellar feedback process, as outflows and heating shape the physical and chemical structure of protostellar envelopes and molecular clouds.

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

We report on synoptic observations at 3.6 and 4.5micron of young stellar objects in IC 348 with 38 epochs covering 40 days. We find that among the detected cluster members, 338 at [3.6] and 269 at both [3.6] and [4.5], many are variable on daily to weekly timescales with typical fluctuations of ~0.1 mag. The fraction of variables ranges from 20% for the diskless pre-main sequence stars to 60% for the stars still surrounded by infalling envelopes. We also find that stars in the exposed cluster core are less variable than the stars in the dense, slightly younger, south-western ridge. This trend persists even after accounting for the underlying correlation with infrared SED type, suggesting that the change in variable fraction is not simply a reflection of the change in relative fraction of class I vs. class II sources across the cloud, but instead reflects a change in variability with age. We also see a strong correlation between infrared variability and X-ray luminosity among the class II sources. The observed variability most likely reflects large changes in the structure of the inner wall located at the dust sublimation radius. We explore the possibility that these structural perturbations could be caused by a hot spot on the star heating dust above the sublimation temperature, causing it to evaporate rapidly, and increasing the inner radius for a portion of the disk. Under a number of simplifying assumptions we show that this model can reproduce the size and timescale of the 3.6 and 4.5micron fluctuations. Regardless of its source, the infrared variability indicates that the inner disk is not a slowly evolving entity, but instead is a bubbling, warped, dented mass of gas and dust whose global size and shape fluctuate in a matter of days.

We extend our previous study of the stellar population of L1641, the lower-density star-forming region of the Orion A cloud south of the dense Orion Nebula Cluster (ONC), with the goal of testing whether there is a statistically significant deficiency of high-mass stars in low-density regions. Previously, we compared the observed ratio of low-mass stars to high-mass stars with theoretical models of the stellar initial mass function (IMF) to infer a deficiency of the highest-mass stars in L1641. We expand our population study to identify the intermediate mass (late B to G) L1641 members in an attempt to make a more direct comparison with the mass function of the nearby ONC. The spectral type distribution and the K-band luminosity function of L1641 are similar to those of the ONC (Hillenbrand 1997; Muench et al. 2002), but problems of incompleteness and contamination prevent us from making a detailed test for differences. We limit our analysis to statistical tests of the ratio of high-mass to low-mass stars, which indicate a probability of only 3% that the ONC and the southern region of L1641 were drawn from the same population, supporting the hypothesis that the upper mass end of the IMF is dependent on environmental density.

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

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

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 infrared observations of the ultra-compact H II region W3(OH) made by the FORCAST instrument aboard SOFIA and by Spitzer/IRAC. We contribute new wavelength data to the spectral energy distribution, which constrains the optical depth, grain size distribution, and temperature gradient of the dusty shell surrounding the H II region. We model the dust component as a spherical shell containing an inner cavity with radius ~ 600 AU, irradiated by a central star of type O9 and temperature ~ 31,000 K. The total luminosity of this system is 71,000 L_solar. An observed excess of 2.2 - 4.5 microns emission in the SED can be explained by our viewing a cavity opening or clumpiness in the shell structure whereby radiation from the warm interior of the shell can escape. We claim to detect the nearby water maser source W3 (H2O) at 31.4 and 37.1 microns using beam deconvolution of the FORCAST images. We constrain the flux densities of this object at 19.7 - 37.1 microns. Additionally, we present in situ observations of four young stellar and protostellar objects in the SOFIA field, presumably associated with the W3 molecular cloud. Results from the model SED fitting tool of Robitaille et al. (2006, 2007} suggest that two objects (2MASS J02270352+6152357 and 2MASS J02270824+6152281) are intermediate-luminosity (~ 236 - 432 L_solar) protostars; one object (2MASS J02270887+6152344) is either a high-mass protostar with luminosity 3000 L_solar or a less massive young star with a substantial circumstellar disk but depleted envelope; and one object (2MASS J02270743+6152281) is an intermediate-luminosity (~ 768 L_solar) protostar nearing the end of its envelope accretion phase or a young star surrounded by a circumstellar disk with no appreciable circumstellar envelope.

IRAS 20050+2720 is young star forming region at a distance of 700 pc without apparent high mass stars. We present results of our multiwavelength study of IRAS 20050+2720 which includes observations by Chandra and Spitzer, and 2MASS and UBVRI photometry. In total, about 300 YSOs in different evolutionary stages are found. We characterize the distribution of young stellar objects (YSOs) in this region using a minimum spanning tree (MST) analysis. We newly identify a second cluster core, which consists mostly of class II objects, about 10 arcmin from the center of the cloud. YSOs of earlier evolutionary stages are more clustered than more evolved objects. The X-ray luminosity function (XLF) of IRAS 20050+2720 is roughly lognormal, but steeper than the XLF of the more massive Orion nebula complex. IRAS 20050+2720 shows a lower N_H/A_K ratio compared with the diffuse ISM.

We have obtained time series observations of the Orion Nebula Cluster at 70 microns and 160 microns from the Herschel/PACS Photometer. This represents the first wide-field far-infrared photometric monitoring of a young star forming region. The acquired 35'x35' maps show complex extended structures, with unprecedented details, that trace the interaction between the molecular gas and the young hot stars. We detect 43 protostars, most of which are situated along the integral-shaped filament extending from the Orion nebula, through OMC2 and to OMC3. We present high-reliability light curves for some of these objects using the first six epochs of our observing program spread over 6 weeks. We find amplitude variations in excess of 20% for a fraction of the detected protostars over periods as short as a few weeks. This is inconsistent with the dynamical time-scales of cool far-IR emitting material that orbits at hundreds of AU from the protostar, and it suggests that the mechanism(s) responsible for the observed variability originates from the inner region of the protostars, likely driven by variable mass accretion.

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

We present results from an optical photometric and spectroscopic survey of the young stellar population in L1641, the low-density star-forming region of the Orion A cloud south of the Orion Nebula Cluster (ONC). Our goal is to determine whether L1641 has a large enough low-mass population to make the known lack of high-mass stars a statistically-significant demonstration of environmental dependence of the upper mass stellar initial mass function (IMF). Our spectroscopic sample consists of IR-excess objects selected from the Spitzer/IRAC survey and non-excess objects selected from optical photometry. We have spectral confirmation of 864 members, with another 98 probable members; of the confirmed members, 406 have infrared excesses and 458 do not. Assuming the same ratio of stars with and without IR excesses in the highly-extincted regions, L1641 may contain as many as ~1600 stars down to ~0.1 solar mass, comparable within a factor of two to the the ONC. Compared to the standard models of the IMF, L1641 is deficient in O and early B stars to a 3-4 sigma significance level, assuming that we know of all the massive stars in L1641. With a forthcoming survey of the intermediate-mass stars, we will be in a better position to make a direct comparison with the neighboring, dense ONC, which should yield a stronger test of the dependence of the high-mass end of the stellar initial mass function upon environment.

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.

We map the full extent of a rich massive young cluster in the Cep OB3b association with the IRAC and MIPS instruments aboard the {\it Spitzer} Space Telescope and the ACIS instrument aboard the $\it{Chandra}$ X-Ray Observatory. At 700 pc, it is revealed to be the second nearest large ($>1000$ member), young ($< 5$ Myr) cluster known. In contrast to the nearest large cluster, the Orion Nebula Cluster, Cep OB3b is only lightly obscured and is mostly located in a large cavity carved out of the surrounding molecular cloud. Our infrared and X-ray datasets, as well as visible photometry from the literature, are used to take a census of the young stars in Cep OB3b. We find that the young stars within the cluster are concentrated in two sub-clusters; an eastern sub-cluster, near the Cep B molecular clump, and a western sub-cluster, near the Cep F molecular clump. Using our census of young stars, we examine the fraction of young stars with infrared excesses indicative of circumstellar disks. We create a map of the disk fraction throughout the cluster and find that it is spatially variable. Due to these spatial variations, the two sub-clusters exhibit substantially different average disk fractions from each other: $32% \pm 4%$ and $50% \pm 6%$. We discuss whether the discrepant disk fractions are due to the photodestruction of disks by the high mass members of the cluster or whether they result from differences in the ages of the sub-clusters. We conclude that the discrepant disk fractions are most likely due to differences in the ages.

We present new mid-infrared images of the central region of the Orion Nebula using the newly commissioned SOFIA airborne telescope and its 5 -- 40 micron camera FORCAST. The 37.1 micron images represent the highest resolution observations (<4") ever obtained of this region at these wavelengths. After BN/KL (which is described in a separate letter in this issue), the dominant source at all wavelengths except 37.1 micron is the Ney-Allen Nebula, a crescent-shaped extended source associated with theta 1D. The morphology of the Ney-Allen nebula in our images is consistent with the interpretation that it is ambient dust swept up by the stellar wind from theta 1D, as suggested by Smith et al. (2005). Our observations also reveal emission from two "proplyds" (proto-planetary disks), and a few embedded young stellar objects (YSOs; IRc9, and OMC1S IRS1, 2, and 10). The spectral energy distribution for IRc9 is presented and fitted with standard YSO models from Robitaille et al. (2007) to constrain the total luminosity, disk size, and envelope size. The diffuse, nebular emission we observe at all FORCAST wavelengths is most likely from the background photodissociation region (PDR) and shows structure that coincides roughly with H_alpha and [N II] emission. We conclude that the spatial variations in the diffuse emission are likely due to undulations in the surface of the background PDR.

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

We present the results of a multi-wavelength multi-epoch survey of five evolved protoplanetary disks in the IC 348 cluster that show significant infrared variability. Using 3-8micron and 24micron photometry along with 5-40micron spectroscopy from the Spitzer Space Telescope, as well as ground-based 0.8-5micron spectroscopy, optical spectroscopy and near-infrared photometry, covering timescales of days to years, we examine the variability in the disk, stellar and accretion flux. We find substantial variations (10-60%) at all infrared wavelengths on timescales of weeks to months for all of these young stellar objects. This behavior is not unique when compared to other cluster members and is consistent with changes in the structure of the inner disk, most likely scale height fluctuations on a dynamical timescale. Previous observations, along with our near-infrared photometry, indicate that the stellar fluxes are relatively constant; stellar variability does not appear to drive the large changes in the infrared fluxes. Based on our near-infrared spectroscopy of the Pa-beta and Br-gamma lines we find that the accretion rates are variable in most of the evolved disks but the overall rates are probably too small to cause the infrared variability. We discuss other possible physical causes for the variability, including the influence of a companion, magnetic fields threading the disk, and X-ray flares.

We report Spitzer observations of five newly identified bow shocks in the massive star-forming region RCW 38. Four are visible at IRAC wavelengths, the fifth is visible only at 24 microns. Chandra X-ray emission indicates that winds from the central O5.5 binary, IRS~2, have caused an outflow to the NE and SW of the central subcluster. The southern lobe of hot ionised gas is detected in X-rays; shocked gas and heated dust from the shock-front are detected with Spitzer at 4.5 and 24 microns. The northern outflow may have initiated the present generation of star formation, based on the filamentary distribution of the protostars in the central subcluster. Further, the bow-shock driving star, YSO 129, is photo-evaporating a pillar of gas and dust. No point sources are identified within this pillar at near- to mid-IR wavelengths. We also report on IRAC 3.6 & 5.8 micron observations of the cluster DBS2003-124, NE of RCW 38, where 33 candidate YSOs are identified. One star associated with the cluster drives a parsec-scale jet. Two candidate HH objects associated with the jet are visible at IRAC and MIPS wavelengths. The jet extends over a distance of ~3 pc. Assuming a velocity of 100 km/s for the jet material gives an age of about 30,000 years, indicating that the star (and cluster) are likely to be very young, with a similar or possibly younger age than RCW 38, and that star formation is ongoing in the extended RCW 38 region.

We present a study of the structure of the high mass star-forming region RCW~38 and the spatial distribution of its young stellar population. Spitzer IRAC photometry 3-8um are combined with 2MASS near-IR data to identify young stellar objects by IR-excess emission from their circumstellar material. Chandra X-ray data are used to identify class III pre-main sequence stars lacking circumstellar material. We identify 624 YSOs: 23 class 0/I and 90 flat spectrum protostars, 437 Class II stars, and 74 Class III stars. We also identify 29 (27 new) O star candidates over the IRAC field. Seventy-two stars exhibit IR-variability, including seven class 0/I and 12 flat spectrum YSOs. A further 177 tentative candidates are identified by their location in the IRAC [3.6] vs. [3.6]-[5.8] cmd. We find strong evidence of subclustering in the region. Three subclusters were identified surrounding the central cluster, with massive and variable stars in each subcluster. The central region shows evidence of distinct spatial distributions of the protostars and pre-main sequence stars. A previously detected IR cluster, DB2001_Obj36, has been established as a subcluster of RCW 38. This suggests that star formation in RCW 38 occurs over a more extended area than previously thought. The gas to dust ratio is examined using the X-ray derived hydrogen column density, N_H and the K-band extinction, and found to be consistent with the diffuse ISM, in contrast with Serpens & NGC1333. We posit that the high photoionising flux of massive stars in RCW 38 affects the agglomeration of the dust grains.

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

We describe extensive synoptic multi-wavelength observations of the transition disk LRLL 31 in the young cluster IC 348. We combined four epochs of IRS spectra, nine epochs of MIPS photometry, seven epochs of cold-mission IRAC photometry and 36 epochs of warm mission IRAC photometry along with multi-epoch near-infrared spectra, optical spectra and polarimetry to explore the nature of the rapid variability of this object. We find that the inner disk, as traced by the 2-5micron excess stays at the dust sublimation radius while the strength of the excess changes by a factor of 8 on weekly timescales, and the 3.6 and 4.5micron photometry shows a drop of 0.35 magnitudes in one week followed by a slow 0.5 magnitude increase over the next three weeks. The accretion rate, as measured by PaBeta and BrGamma emission lines, varies by a factor of five with evidence for a correlation between the accretion rate and the infrared excess. While the gas and dust in the inner disk are fluctuating the central star stays relatively static. Our observations allow us to put constraints on the physical mechanism responsible for the variability. The variabile accretion, and wind, are unlikely to be causes of the variability, but both are effects of the same physical process that disturbs the disk. The lack of periodicity in our infrared monitoring indicates that it is unlikely that there is a companion within ~0.4 AU that is perturbing the disk. The most likely explanation is either a companion beyond ~0.4 AU or a dynamic interface between the stellar magnetic field and the disk leading to a variable scale height and/or warping of the inner disk.

We present a global study of low mass, young stellar object (YSO) surface densities in nearby (< 500 pc) star forming regions based on a comprehensive collection of Spitzer Space Telescope surveys. We show that the distribution of YSO surface densities in the solar neighbourhood is a smooth distribution, being adequately described by a lognormal function from a few to 10^3 YSOs per pc^2, with a peak at 22 stars/pc^2 and a dispersion of 0.85. We do not find evidence for multiple discrete modes of star-formation (e.g. clustered and distributed). Comparing the observed surface density distribution to previously reported surface density threshold definitions of clusters, we find that the fraction of stars in clusters is crucially dependent on the adopted definitions, ranging from 40 to 90%. However, we find that only a low fraction (< 26%) of stars are formed in dense environments where their formation/evolution (along with their circumstellar disks and/or planets) may be affected by the close proximity of their low-mass neighbours.

We present new images and photometry of the massive star forming complex Cygnus X obtained with IRAC and MIPS on the Spitzer Space Telescope. A combination of IRAC, MIPS, UKIDSS, and 2MASS data are used to identify and classify young stellar objects. Of the 8,231 sources detected exhibiting infrared excess in Cygnus X North, 670 are classified as Class I and 7,249 are classified as Class II. Using spectra from the FAST spectrograph at the Fred L. Whipple Observatory and Hectospec on the MMT, we spectrally typed 536 sources in the Cygnus X complex to identify the massive stars. We find that YSOs tend to be grouped in the neighborhoods of massive B stars (spectral types B0 to B9). We present a minimal spanning tree analysis of clusters in two regions in Cygnus X North. The fraction of infrared excess sources that belong to clusters with >10 members is found to be 50--70%. Most Class II objects lie in dense clusters within blown out HII regions, while Class I sources tend to reside in more filamentary structures along the bright-rimmed clouds, indicating possible triggered star formation.

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.

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

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

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

We present observations of newly discovered 24 micron circumstellar structures detected with the Multiband Imaging Photometer for Spitzer (MIPS) around three evolved stars in the Cygnus-X star forming region. One of the objects, BD+43 3710, has a bipolar nebula, possibly due to an outflow or a torus of material. A second, HBHA 4202-22, a Wolf-Rayet candidate, shows a circular shell of 24 micron emission suggestive of either a limb-brightened shell or disk seen face-on. No diffuse emission was detected around either of these two objects in the Spitzer 3.6-8 micron Infrared Array Camera (IRAC) bands. The third object is the luminous blue variable candidate G79.29+0.46. We resolved the previously known inner ring in all four IRAC bands. The 24 micron emission from the inner ring extends ~1.2 arcmin beyond the shorter wavelength emission, well beyond what can be attributed to the difference in resolutions between MIPS and IRAC. Additionally, we have discovered an outer ring of 24 micron emission, possibly due to an earlier episode of mass loss. For the two shell stars, we present the results of radiative transfer models, constraining the stellar and dust shell parameters. The shells are composed of amorphous carbon grains, plus polycyclic aromatic hydrocarbons in the case of G79.29+0.46. Both G79.29+0.46 and HBHA 4202-22 lie behind the main Cygnus-X cloud. Although G79.29+0.46 may simply be on the far side of the cloud, HBHA 4202-22 is unrelated to the Cygnus-X star formation region.

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.

[abridged] Protoplanetary disks with AU-scale inner clearings, often referred to as transitional disks, provide a unique sample for understanding disk dissipation mechanisms and possible connections to planet formation. Observations of young stellar clusters with the Spitzer Space Telescope have amassed mid-infrared spectral energy distributions for thousands of star-disk systems from which transition disks can be identified. From a sample of 8 relatively nearby young regions (d <= 400 pc), we have identified about 20 such objects, which we term "classical" transition disks, spanning a wide range of stellar age and mass. We also identified two additional categories representing more ambiguous cases: "warm excess" objects with transition-like spectral energy distributions but moderate excess at 5.8 microns, and "weak excess" objects with smaller 24 micron excess that may be optically thin or exhibit advanced dust grain growth and settling. From existing Halpha emission measurements, we find evidence for different accretion activity among the three categories, with a majority of the classical and warm excess transition objects still accreting gas through their inner holes and onto the central stars, while a smaller fraction of the weak transition objects are accreting at detectable rates. We find a possible age dependence to the frequency of classical transition objects, with fractions relative to the total population of disks in a given region of a few percent at 1-2 Myr rising to 10-20% at 3-10 Myr. The trend is even stronger if the weak and warm excess objects are included. Classical transition disks appear to be less common, and weak transition disks more common, around lower-mass stars (M <= 0.3 Msun).

We present multi-epoch Spitzer Space Telescope observations of the transitional disk LRLL 31 in the 2-3 Myr-old star forming region IC 348. Our measurements show remarkable mid-infrared variability on timescales as short as one week. The infrared continuum emission exhibits systematic wavelength-dependent changes that suggest corresponding dynamical changes in the inner disk structure and variable shadowing of outer disk material. We propose several possible sources for the structural changes, including a variable accretion rate or a stellar or planetary companion embedded in the disk. Our results indicate that variability studies in the infrared can provide important new constraints on protoplanetary disk behavior.

We present results of a sensitive Chandra X-ray observation and Spitzer mid-IR observations of the infrared cluster lying north of the NGC 2071 reflection nebula in the Orion B molecular cloud. We focus on the dense cluster core known as NGC 2071-IR which contains at least nine IR sources within a 40 x 40 arcsecond region. This region shows clear signs of active star formation including powerful molecular outflows, Herbig-Haro objects, and both OH and H2O masers. We use Spitzer IRAC images to aid in X-ray source identification and to determine YSO classes using mid-IR colors. Spitzer IRAC colors show that the luminous source IRS 1 is a class I protostar. IRS 1 is believed to be driving a powerful bipolar molecular outflow and may be an embedded B-type star or its progenitor. Its X-ray spectrum reveals a fluorescent Fe emission line at 6.4 keV, arising in cold material near the protostar. The line is present even in the absence of large flares, raising questions about the nature of the ionizing mechanism responsible for producing the 6.4 keV fluorescent line. Chandra also detects X-ray sources at or near the positions of IRS 2, IRS 3, IRS 4, and IRS 6 and a variable X-ray source coincident with the radio source VLA 1, located just 2 arcsec north of IRS 1. No IR data are yet available to determine a YSO classification for VLA 1, but its high X-ray absorption shows that it is even more deeply-embedded than IRS 1, suggesting that it could be an even younger, less-evolved protostar.

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.

Fomalhaut and epsilon Eridani are two young, nearby stars that possess extended debris disks whose structures suggest the presence of perturbing planetary objects. With its high sensitivity and stable point spread function, Spitzer/IRAC is uniquely capable of detecting cool, Jupiter-like planetary companions whose peak emission is predicted to occur near 4.5 um. We report on deep IRAC imaging of these two stars, taken at 3.6 and 4.5 um using subarray mode and in all four channels in wider-field full array mode. Observations acquired at two different telescope roll angles allowed faint surrounding objects to be separated from the stellar diffraction pattern. No companion candidates were detected at the reported position of Fomalhaut b with 3 sigma model-dependent mass upper limits of 3 MJ (for an age of 200 Myr). Around epsilon Eridani we instead set a limit of 4 and <1 MJ (1 Gyr model age) at the inner and outer edge of the sub-millimeter debris ring, respectively. These results are consistent with non-detections in recent near-infrared imaging searches, and set the strongest limits to date on the presence of planets outside epsilon Eridani sub-millimeter ring.

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 adaptive optics (AO) near-infrared (JHKs) observations of the deeply embedded massive cluster RCW 38 using NACO on the VLT. Narrowband AO observations centered at wavelengths of 1.28, 2.12, and 2.17 micron were also obtained. The area covered by these observations is about 0.5 pc square, centered on the O star RCW 38 IRS2. We use the JHKs colors to identify young stars with infrared excess. Through a detailed comparison to a nearby control field, we find that most of the 337 stars detected in all three infrared bands are cluster members (~317), with essentially no contamination due to background or foreground sources. Five sources have colors suggestive of deeply embedded protostars, while 53 sources are detected at Ks only; their spatial distribution with respect to the extinction suggests they are highly reddened cluster members. Detectable Ks-band excess is found toward 29 +/- 3 % of the stars. For comparison to a similar area of Orion observed in the near-infrared, mass and extinction cuts are applied, and the excess fractions redetermined. The resulting excesses are then 25 +/- 5 % for RCW 38, and 42 +/- 8 % for Orion. RCW 38 IRS2 is shown to be a massive star binary with a projected separation of ~500 AU. Two regions of molecular hydrogen emission are revealed through the 2.12 micron imaging. One shows a morphology suggestive of a protostellar jet, and is clearly associated with a star only detected at H and Ks, previously identified as a highly obscured X-ray source. Three spatially extended cometary-like objects, suggestive of photoevaporating disks, are identified, but only one is clearly directly influenced by RCW 38 IRS2. A King profile provides a reasonable fit to the cluster radial density profile and a nearest neighbor distance analysis shows essentially no sub-clustering.

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 a kinematic study of the Orion Nebula Cluster based upon radial velocities measured by multi-fiber echelle spectroscopy at the 6.5 meter MMT and Magellan telescopes. Velocities are reported for 1613 stars, with multi-epoch data for 727 objects as part of our continuing effort to detect and analyze spectroscopic binaries. We confirm and extend the results of Furesz et al. showing that the ONC is not relaxed, consistent with its youth, and that the stars generally follow the position-velocity structure of the moderate density gas in the region, traced by $^{13}$CO. The additional radial velocities we have measured enable us to probe some discrepancies between stellar and gaseous structure which can be attributed to binary motion and the inclusion of non-members in our kinematic sample. Our multi-epoch data allow us to identify 89 spectroscopic binaries; more will be found as we continue monitoring. Our results reinforce the idea that the ONC is a cluster in formation, and thus provides a valuable testing ground for theory. In particular, our observations are not consistent with the quasi-equilibrium or slow contraction models of cluster formation, but are consistent with cold collapse models.

Most stars are born in rich young stellar clusters (YSCs) embedded in giant molecular clouds. The most massive stars live out their short lives there, profoundly influencing their natal environments by ionizing HII regions, inflating wind-blown bubbles, and soon exploding as supernovae. Thousands of lower-mass pre-main sequence stars accompany the massive stars, and the expanding HII regions paradoxically trigger new star formation as they destroy their natal clouds. While this schematic picture is established, our understanding of the complex astrophysical processes involved in clustered star formation have only just begun to be elucidated. The technologies are challenging, requiring both high spatial resolution and wide fields at wavelengths that penetrate obscuring molecular material and remove contaminating Galactic field stars. We outline several important projects for the coming decade: the IMFs and structures of YSCs; triggered star formation around YSC; the fate of OB winds; the stellar populations of Infrared Dark Clouds; the most massive star clusters in the Galaxy; tracing star formation throughout the Galactic Disk; the Galactic Center region and YSCs in the Magellanic Clouds. Programmatic recommendations include: developing a 30m-class adaptive optics infrared telescope; support for high-resolution and wide field X-ray telescopes; large-aperture sub-millimeter and far-infrared telescopes; multi-object infrared spectrographs; and both numerical and analytical theory.

The giant molecular cloud G216-2.5, also known as Maddalena's cloud or the Maddalena-Thaddeus cloud, is distinguished by an unusual combination of high gas mass (1-6 x 10^5) solar masses, low kinetic temperatures (10 K), and the lack of bright far infrared emission. Although star formation has been detected in neighboring satellite clouds, little evidence for star formation has been found in the main body of this cloud. Using a combination of mid-infrared observations with the IRAC and MIPS instruments onboard the Spitzer space telescope, and near-IR images taken with the Flamingos camera on the KPNO 2.1-meter, we identify a population of 41 young stars with disks and 33 protostars in the center of the cloud. Most of the young stellar objects are coincident with a filamentary structure of dense gas detected in CS (2-1). These observations show that the main body of G216 is actively forming stars, although at a low stellar density comparable to that found in the Taurus cloud.

At ~ 400 pc, the Horsehead Nebula (B33) is the closest radiatively-sculpted pillar to the Sun, but the state and extent of star formation in this structure is not well understood. We present deep near-infrared (IRSF/SIRIUS JHKs) and mid-infrared (Spitzer/IRAC) observations of the Horsehead Nebula in order to characterize the star forming properties of this region and to assess the likelihood of triggered star formation. Infrared color-color and color-magnitude diagrams are used to identify young stars based on infrared excess emission and positions to the right of the Zero-Age Main Sequence, respectively. Of the 45 sources detected at both near- and mid-infrared wavelengths, three bona fide and five candidate young stars are identified in this 7' by 7' region. Two bona fide young stars have flat infrared SEDs and are located at the western irradiated tip of the pillar. The spatial coincidence of the protostars at the leading edge of this elephant trunk is consistent with the Radiation-Driven Implosion (RDI) model of triggered star formation. There is no evidence, however, for sequential star formation within the immediate ~ 1.5' (0.17 pc) region from the cloud/H II region interface.

We present wide-field near-IR images of Orion A. K and H2 1-0S(1) images of a contiguous 8 sqr degree region are compared to photometry from Spitzer and dust-continuum maps obtained with MAMBO and SCUBA. We also measure proper motions for H2 features in 33 outflows. We increase the number of known H2 outflows in Orion A to 116. A total of 111 H2 flows were observed with Spitzer; outflow sources are identified for at least 72 of them. The MAMBO 1200 micron maps cover 97 H2 flows; 57 of them are associated with Spitzer sources and dust cores or extended emission. The H2 jets are widely distributed and randomly orientated; the jets do not appear to be orthogonal to large-scale filaments or even to the small-scale cores. Moreover, H2 jet lengths and opening angles are not obviously correlated with indicators of outflow source age - source spectral index or (sub)millimetre core flux. We demonstrate that H2 jet sources are predominantly protostellar with flat or positive spectral indices, rather than disk-excess (or T Tauri) stars. Most protostars in molecular cores drive H2 outflows. However, not all molecular cores are associated with protostars or H2 jets. On statistical grounds, the H2 jet phase may be marginally shorter than the protostellar phase, though must be considerably shorter than the prestellar phase. In terms of their spectral index, H2 jet sources are indistinguishable from protostars. The few true protostars without H2 jets are almost certainly more evolved than their H2-jet-driving counterparts. We also find that protostars that power molecular outflows are no more (nor no less) clustered than protostars that do not. The H2 emission regions in outflows from young stars clearly weaken and fade very quickly, before the source evolves from protostar to pre-main-sequence star.

We present several corrections for point source photometry to be applied to data from the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. These corrections are necessary because of characteristics of the IRAC arrays and optics and the way the instrument is calibrated in-flight. When these corrections are applied, it is possible to achieve a ~2% relative photometric accuracy for sources of adequate signal to noise in an IRAC image.

We present Spitzer Space Telescope time series photometry of the exoplanet system HD 189733 spanning two times of secondary eclipse, when the planet passes out of view behind the parent star. We estimate the relative eclipse depth in 5 distinct bands and find the planet-to-star flux ratio to be 0.256 +/- 0.014% (3.6 microns), 0.214 +/- 0.020% (4.5 microns), 0.310 +/- 0.034% (5.8 microns), 0.391 +/- 0.022% (8.0 microns), and 0.598 +/- 0.038% (24 microns). For consistency, we re-analyze a previously published time series to deduce a contrast ratio in an additional band, 0.519 +/- 0.020% (16 microns). Our data are strongly inconsistent with a Planck spectrum, and we clearly detect emission near 4 microns as predicted by published theoretical models in which this feature arises from a corresponding opacity window. Unlike recent results for the exoplanet HD 209458b, we find that the emergent spectrum from HD 189733b is best matched by models that do not include an atmospheric temperature inversion. Taken together, these two studies provide initial observational support for the idea that hot Jupiter atmospheres diverge into two classes, in which a thermal inversion layer is present for the more strongly irradiated objects.

Forty new low mass members with spectral types ranging from M4-M9 have been confirmed in the Orion Molecular Cloud 2/3 region. Through deep, I, z', J, H, K photometry of a 20' x 20' field in OMC 2/3, we selected brown dwarf candidates for follow-up spectroscopy. Low resolution far-red and near-infrared spectra were obtained for the candidates, and 19 young brown dwarfs in the OMC 2/3 region are confirmed. They exhibit spectral types of M6.5-M9, corresponding to approximate masses of 0.075-0.015 M_solar using the evolutionary models of Baraffe et al. (1998). At least one of these bona fide young brown dwarfs has strong Halpha emission, indicating that it is actively accreting. In addition, we confirm 21 new low mass members with spectral types of M4-M6, corresponding to approximate masses of 0.35-0.10 M_solar in OMC 2/3. By comparing pre-main sequence tracks to the positions of the members in the H-R diagram, we find that most of the brown dwarfs are less than 1 Myr, but find a number of low mass stars with inferred ages greater than 3 Myr. The discrepancy in the stellar and substellar ages is due to our selection of only low luminosity sources; however, the presence of such objects implies the presence of an age spread in the OMC 2/3 region. We discuss possible reasons for this apparent age spread.

We present Spitzer Space Telescope IRAC and MIPS observations toward a sample of nine high-mass star forming regions at a distance of around 2 kpc. Based on IRAC and MIPS 24 $\mu$m photometric results and 2MASS JHKs data, we carry out a census of young stellar objects (YSOs) in a 5' by 5' field toward each region. Toward seven out of the nine regions, we detect parsec sized clusters with around 20 YSOs surrounded by a more extended and sparse distribution of young stars and protostars. For the other two regions, IRAS 20126+4104 and IRAS 22172+5549, the former has the lowest number of YSOs in the sample and shows no obvious cluster, and the latter appears to be part of a larger, potentially more evolved cluster. The deep IRAC imaging reveals at least twelve outflows in eight out of the nine regions, with nine outflows prominent in the 4.5 $\mu$m band most probably attributed to shocked H$_2$ emission, two outflows dominated by scattered light in the 3.6 and 4.5 $\mu$m bands, and one outflow standing out from its hydrocarbon emission in the 8.0 $\mu$m band. In comparison with previous ground-based observations, our IRAC observations reveal new outflow structures in five regions. The dramatically different morphologies of detected outflows can be tentatively interpreted in terms of possible evolution of massive outflows. The driving sources of these outflows are deeply embedded in dense dusty cores revealed by previous millimeter interferometric observations. We detect infrared counterparts of these dusty cores in the IRAC or MIPS 24 $\mu$m bands. Reflection nebulae dominated by the emission from UV heated hydrocarbons in the 8 $\mu$m band can be found in most regions and they may imply the presence of young B stars.

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.

Aims: Investigation of the dense gas, the outflows and the continuum emission from the massive twin cores NGC6334I and I(N) at high spatial resolution. Methods: We imaged the region with the Australia Telescope Compact Array (ATCA) at 3.4mm wavelength in continuum as well as CH3CN(5_K-4_K) and HCN(1-0) spectral line emission. Results: While the continuum emission in NGC6334I mainly traces the UCHII region, toward NGC6334I(N) we detect line emission from four of the previously identified dust continuum condensations that are of protostellar or pre-stellar nature. The CH3CN(5_K-4_K) lines are detected in all K-components up to energies of 128K above ground toward two protostellar condensations in both regions. We find line-width increasing with increasing K for all sources, which indicates a higher degree of internal motions closer to the central protostars. Toward the main mm and CH3CN source in NGC6334I we identify a velocity gradient approximately perpendicular to the large-scale molecular outflow. This may be interpreted as a signature of an accretion disk, although other scenarios, e.g., an unresolved double source, could produce a similar signature as well. No comparable signature is found toward any of the other sources. HCN does not trace the dense gas well but it is dominated by the molecular outflows. While the outflow in NGC6334I exhibits a normal Hubble-law like velocity structure, the data indicate a precessing outflow close to the plane of the sky for NGC6334I(N). Furthermore, we observe a wide (~15.4km/s) HCN absorption line, much broader than the previously observed CH3OH and NH3 absorption lines. Several explanations for the difference are discussed.

During the last two decades, the focus of star formation research has shifted from understanding the collapse of a single dense core into a star to studying the formation hundreds to thousands of stars in molecular clouds. In this chapter, we overview recent observational and theoretical progress toward understanding star formation on the scale of molecular clouds and complexes, i.e the macrophysics of star formation. We begin with an overview of recent surveys of young stellar objects (YSOs) in molecular clouds and embedded clusters, and we outline an emerging picture of cluster formation. We then discuss the role of turbulence to both support clouds and create dense, gravitationally unstable structures, with an emphasis on the role of magnetic fields (in the case of distributed stars) and feedback (in the case of clusters) to slow turbulent decay and mediate the rate and density of star formation. The discussion is followed by an overview of how gravity and turbulence may produce observed scaling laws for the properties of molecular clouds, stars and star clusters, and how the observed, low star formation rate may result from self regulated star formation. We end with some concluding remarks, including a number of questions to be addressed by future observations and simulations.

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 results from 1351 high resolution spectra of 1215 stars in the Orion Nebula Cluster (ONC) and the surrounding Orion 1c association, obtained with the Hectochelle multiobject echelle spectrograph on the 6.5m MMT. We confirmed 1111 stars as members, based on their radial velocity and/or H-alpha emission. The radial velocity distribution of members shows a dispersion of 3.1 km/s. We found a substantial north-south velocity gradient and spatially coherent structure in the radial velocity distribution, similar to that seen in the molecular gas in the region. We also identified several binary and high velocity stars, a region exhibiting signs of triggered star formation, and a possible foreground population of stars somewhat older than the ONC. Stars without infrared excesses (as detected with the IRAC instrument on the Spitzer Space Telescope) exhibit a wider spread in radial velocity than the infrared excess stars; this spread is mostly due to a blue-shifted population of stars that may constitute a foreground population. We also identify some accreting stars, based on H-alpha, that do not have detectable infrared excesses with IRAC, and thus are potential transitional disk systems (objects with inner disk holes). We propose that the substructure seen both in stellar and gaseous component is the result of non-uniform gravitational collapse to a filamentary distribution of gas. The spatial and kinematic correlation between the stellar and gaseous components suggests the region is very young, probably only ~1 crossing time old or less to avoid shock dissipation and gravitational interactions which would tend to destroy the correlation between stars and gas.

We estimate the strength of the bandpass-integrated thermal emission from the extrasolar planet HD 209458b at 3.6, 4.5, 5.8, and 8.0 microns using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We observe a single secondary eclipse simultaneously in all four bandpasses and find relative eclipse depths of 0.00094 +/- 0.00009, 0.00213 +/- 0.00015, 0.00301 +/- 0.00043, and 0.00240 +/- 0.00026, respectively. These eclipse depths reveal that the shape of the inferred emission spectrum for the planet differs significantly from the predictions of standard atmosphere models; instead the most plausible explanation would require the presence of an inversion layer high in the atmosphere leading to significant water emission in the 4.5 and 5.8 micron bandpasses. This is the first clear indication of such a temperature inversion in the atmosphere of a hot Jupiter, as previous observations of other planets appeared to be in reasonably good agreement with the predictions of models without such an inversion layer.

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.

"Hot Jupiter" extrasolar planets are expected to be tidally locked because they are close (<0.05 astronomical units, where 1 AU is the average Sun-Earth distance) to their parent stars, resulting in permanent daysides and nightsides. By observing systems where the planet and star periodically eclipse each other, several groups have been able to estimate the temperatures of the daysides of these planets. A key question is whether the atmosphere is able to transport the energy incident upon the dayside to the nightside, which will determine the temperature at different points on the planet's surface. Here we report observations of HD 189733, the closest of these eclipsing planetary systems, over half an orbital period, from which we can construct a 'map' of the distribution of temperatures. We detected the increase in brightness as the dayside of the planet rotated into view. We estimate a minimum brightness temperature of 973 +/- 33 K and a maximum brightness temperature of 1212 +/- 11 K at a wavelength of 8 microns, indicating that energy from the irradiated dayside is efficiently redistributed throughout the atmosphere, in contrast to a recent claim for another hot Jupiter. Our data indicate that the peak hemisphere-integrated brightness occurs 16$\pm$6 degrees before opposition, corresponding to a hot spot shifted east of the substellar point. The secondary eclipse (when the planet moves behind the star) occurs 120 +/- 24 s later than predicted, which may indicate a slightly eccentric orbit.

During the last three decades, evidence has mounted that star and planet formation is not an isolated process, but is influenced by current and previous generations of stars. Although cool stars form in a range of environments, from isolated globules to rich embedded clusters, the influences of other stars on cool star and planet formation may be most significant in embedded clusters, where hundreds to thousands of cool stars form in close proximity to OB stars. At the cool stars 14 meeting, a splinter session was convened to discuss the role of environment in the formation of cool stars and planetary systems; with an emphasis on the ``hot'' environment found in rich clusters. We review here the basic results, ideas and questions presented at the session. We have organized this contribution into five basic questions: what is the typical environment of cool star formation, what role do hot star play in cool star formation, what role does environment play in planet formation, what is the role of hot star winds and supernovae, and what was the formation environment of the Sun? The intention is to review progress made in addressing each question, and to underscore areas of agreement and contention.

We present an independent estimate of the interstellar extinction law for the Spitzer IRAC bands as well as a first attempt at extending the law to the 24micron MIPS band. The source data for these measurements are observations of five nearby star-forming regions: the Orion A cloud, NGC 2068/71, NGC 2024/23, Serpens and Ophiuchus. Color excess ratios E(H-Ks)/E(Ks-[lambda]) were measured for stars without infrared excess dust emission from circumstellar disks/envelopes. For four of these five regions, the extinction laws are similar at all wavelengths and differ systematically from a previous determination of the extinction law, which was dominated by the diffuse ISM, derived for the IRAC bands. This difference could be due to the difference in the dust properties of the dense molecular clouds observed here and those of the diffuse ISM. The extinction law at longer wavelengths toward the Ophiuchus region lies between that to the other four regions studied here and that for the ISM. In addition, we extended our extinction law determination to 24micron for Serpens and NGC 2068/71 using Spitzer MIPS data. We compare these results against several ISO extinction law determinations, although in each case there are assumptions which make absolute comparison uncertain. However, our work confirms a relatively flatter extinction curve from 4 - 8micron than the previously assumed standard, as noted by all of these recent studies. The extinction law at 24micron is consistent with previous measurements and models, although there are relatively large uncertainties.

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.

Aims: The massive twin cores NGC6334I and I(N) are in different evolutionary stages and hence ideal targets to study evolutionary variations within the same larger-scale environment. Here, we study the warm, compact gas components. Methods: We imaged the two regions with the Australia Telescope Compact Array (ATCA) at high angular resolution in the NH3(3,3) to (6,6) inversion lines. Results: Compact emission is detected toward both regions in all observed inversion lines with energy levels up to 407K above ground. This is particularly surprising for NGC6334I(N) since it lacks bright infrared emission and is considered a massive cold core at an early evolutionary stage. High optical depth and multiply-peaked line profiles complicate rotation temperature estimates, and we can only conclude that gas components with temperatures >100K are present in both regions. Toward NGC6334I, we confirm previous reports of NH3(3,3) maser emission toward the outflow bow-shocks. Furthermore, we report the first detection of an NH3(6,6) maser toward the central region of NGC6334I. This maser is centered on the second millimeter (mm) peak and elongated along the outflow axis, indicating that this mm continuum core harbors the driving source of the molecular outflow. Toward the main mm peak in NGC6334I(N), we detect a double-horn line profile in the NH3(6,6) transition. The current data do not allow us to differentiate whether this double-horn profile is produced by multiple gas components along the line of sight, or whether it may trace a potential underlying massive accretion disk. The data to Figures 3 to 7 are also available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/.

We report the discovery of T dwarf companions to the nearby stars HN Peg (G0V, 18.4 pc, ~0.3 Gyr) and HD 3651 (K0V, 11.1 pc, ~7 Gyr). During an ongoing survey of 5'x5' fields surrounding stars in the solar neighborhood with IRAC aboard the Spitzer Space Telescope, we identified these companions as candidate T dwarfs based on their mid-IR colors. Using near-IR spectra obtained with SpeX at the NASA IRTF, we confirm the presence of methane absorption that characterizes T dwarfs and measure spectral types of T2.5+/-0.5 and T7.5+/-0.5 for HN Peg B and HD 3651 B, respectively. By comparing our Spitzer data to images from 2MASS obtained several years earlier, we find that the proper motions of HN Peg B and HD 3651 B are consistent with those of the primaries, confirming their companionship. HN Peg B and HD 3651 B have angular separations of 43.2" and 42.9" from their primaries, which correspond to projected physical separations of 795 and 476 AU, respectively. A comparison of their luminosities to the values predicted by theoretical evolutionary models implies masses of 0.021+/-0.009 and 0.051+/-0.014 Msun for HN Peg B and HD 3651 B. In addition, the models imply an effective temperature for HN Peg B that is significantly lower than the values derived for other T dwarfs at similar spectral types, which is the same behavior reported by Metchev & Hillenbrand for the young late-L dwarf HD 203030 B. Thus, the temperature of the L/T transition appears to depend on surface gravity. Meanwhile, HD 3651 B is the first substellar companion directly imaged around a star that is known to harbor a close-in planet from RV surveys. The discovery of this companion supports the notion that the high eccentricities of close-in planets like the one near HD 3651 may be the result of perturbations by low-mass companions at wide separations.

We have conducted a photometric monitoring program of 3 field late-L brown dwarfs looking for evidence of non-axisymmetric structure or temporal variability in their photospheres. The observations were performed using Spitzer/IRAC 4.5 and 8 micron bandpasses and were designed to cover at least one rotational period of each object. One-sigma RMS (root mean squared) uncertainties of less than 3 mmag at 4.5 micron and around 9 mmag at 8 micron were achieved. Two out of the three objects studied exhibit some modulation in their light curves at 4.5 micron - but not 8 micron - with periods of 7.4 hr and 4.6 hr and peak-to-peak amplitudes of 10 mmag and 8 mmag. Although the lack of detectable 8 micron variation suggests an instrumental origin for the detected variations, the data may nevertheless still be consistent with intrinsic variability since the shorter wavelength IRAC bandpasses probe more deeply into late L dwarf atmospheres than the longer wavelengths. A cloud feature occupying a small percentage (1-2 %) of the visible hemisphere could account for the observed amplitude of variation. If, instead, the variability is indeed instrumental in origin, then our non-variable L dwarfs could be either completely covered with clouds or objects whose clouds are smaller and uniformly distributed. Such scenarios would lead to very small photometric variations. Followup IRAC photometry at 3.6 and 5.8 micron bandpasses should distinguish between the two cases. In any event, the present observations provide the most sensitive search to date for structure in the photospheres of late-L dwarfs at mid-IR wavelengths, and our photometry provides stringent upper limits to the extent to which the photospheres of these transition L dwarfs are structured.

We present the results of a program to acquire photometry for eighty-six late-M, L, and T dwarfs using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We examine the behavior of these cool dwarfs in various color-color and color-magnitude diagrams composed of near-IR and IRAC data. The T dwarfs exhibit the most distinctive positions in these diagrams. In M_5.8 versus [5.8]-[8.0], the IRAC data for T dwarfs are not monotonic in either magnitude or color, giving the clearest indication yet that the T dwarfs are not a one parameter family in Teff. Because metallicity does not vary enough in the solar neighborhood to act as the second parameter, the most likely candidate then is gravity, which in turn translates to mass. Among objects with similar spectral type, the range of mass suggested by our sample is about a factor of five (~70 M_Jup to ~15 M_Jup), with the less massive objects making up the younger members of the sample. We also find the IRAC 4.5 micron fluxes to be lower than expected, from which we infer a stronger CO fundamental band at ~4.67 microns. This suggests that equilibrium CH_4/CO chemistry underestimates the abundance of CO in T dwarf atmospheres, confirming earlier results based on M-band observations from the ground. In combining IRAC photometry with near-IR JHK photometry and parallax data, we find the combination of Ks, IRAC 3.6 micron, and 4.5 micron bands to provide the best color-color discrimination for a wide range of M, L, and T dwarfs. Also noteworthy is the M_Ks versus Ks-[4.5] relation, which shows a smooth progression over spectral type and splits the M, L, and T types cleanly.

Using the Submillimeter Array (SMA), we have imaged the 1.3 millimeter continuum emission at the center of the massive star-forming regions NGC 6334 I and I(N). In both regions, the SMA observations resolve the emission into multiple millimeter sources, with most of the sources clustered into areas only 10,000 AU in diameter. Toward NGC 6334 I, we find four compact sources: the two brightest (I-SMA1 and I-SMA2) are associated with previously-known ammonia cores; I-SMA3 is coincident with the peak of the compact HII region (NGC 6334 F); and I-SMA4 is a newly-discovered object. While I-SMA3 exhibits a mixture of free-free and dust emission, the rest of the objects are dust cores. Toward NGC 6334 I(N), seven compact dust cores are found, one of which is associated with a faint centimeter source. With the exception of I-SMA3, none of the millimeter sources have infrared counterparts in Spitzer Space Telescope 3-8 micron images. Using a simple physical model for the dust continuum emission, the estimated mass of interstellar material toward each of these compact objects is in the range of 3 to 66 solar masses. The total mass in the compact objects appears to be similar in I and I(N). The small size of these groups of sources suggest that these objects are proto-Trapezia forming in the centers of clusters of low to intermediate-mass stars.

We have used the InfraRed Array Camera (IRAC) onboard the Spitzer Space telescope to search for low mass companions of the nearby debris disk star epsilon Eridani. The star was observed in two epochs 39 days apart, with different focal plane rotation to allow the subtraction of the instrumental Point Spread Function, achieving a maximum sensitivity of 0.01 MJy/sr at 3.6 and 4.5 um, and 0.05 MJy/sr at 5.8 and 8.0 um. This sensitivity is not sufficient to directly detect scattered or thermal radiation from the epsilon Eridani debris disk. It is however sufficient to allow the detection of Jovian planets with mass as low as 1 MJ in the IRAC 4.5 um band. In this band, we detected over 460 sources within the 5.70 arcmin field of view of our images. To test if any of these sources could be a low mass companion to epsilon Eridani, we have compared their colors and magnitudes with models and photometry of low mass objects. Of the sources detected in at least two IRAC bands, none fall into the range of mid-IR color and luminosity expected for cool, 1 Gyr substellar and planetary mass companions of epsilon Eridani, as determined by both models and observations of field M, L and T dwarf. We identify three new sources which have detections at 4.5 um only, the lower limit placed on their [3.6]-[4.5] color consistent with models of planetary mass objects. Their nature cannot be established with the currently available data and a new observation at a later epoch will be needed to measure their proper motion, in order to determine if they are physically associated to epsilon Eridani.

Early observations of T Tauri stars suggested that stars with evidence of circumstellar accretion disks rotated slower than stars without such evidence, but more recent results are not as clear. Near-IR circumstellar disk indicators, though the most widely available, are subject to uncertainties that can result from inner disk holes and/or the system inclination. Mid-infrared observations are less sensitive to such effects, but until now, these observations have been difficult to obtain. The Spitzer Space Telescope now easily enables mid-infrared measurements of large samples of PMS stars covering a broad mass range in nearby star-forming regions. Megeath and collaborators surveyed the Orion Molecular Clouds (~1 Myr) with the IRAC instrument (3.6, 4.5, 5.8, 8 microns) as part of a joint IRAC and MIPS GTO program. We examine the relationship between rotation and Spitzer mid-IR fluxes for ~900 stars in Orion for stars between 3 and 0.1 Msun. We find in these Spitzer data the clearest indication to date that stars with longer periods are more likely than those with short periods to have IR excesses suggestive of disks.

We present observations of six late-type members of the young cluster IC 348 detected at 24 microns with the Multiband Imaging Photometer for Spitzer(MIPS). At least four of the objects are probably substellar. Combining these data with ground-based optical and near-infrared photometry and complementary observations with the Infrared Array Camera (IRAC), we have modeled the spectral energy distributions using detailed models of irradiated accretion disks. We are able to fit the observations with models using a range of maximum grain sizes from ISM-type dust to grains as large as 1 millimeter. Two objects show a lack of excess emission at wavelengths shortward of 5.8-8 microns but significant excess at longer wavelengths, indicative of large optically thin or evacuated inner holes. Our models indicate a inner hole of radius ~ 0.5-0.9 AU for the brown dwarf L316; this is the first brown dwarf with evidence for an AU-scale inner disk hole. We examine several possible mechanisms for the inner disk clearing in this case, including photoevaporation and planet formation.

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.

Using the Hubble Space Telescope, the 4 m Blanco telescope at the Cerro Tololo Inter-American Observatory, and the Spitzer Space Telescope, we have performed deep imaging from 0.8 to 8 um of the southern subcluster in the Chamaeleon I star-forming region. In these data, we have discovered an object, Cha 110913-773444, whose colors and magnitudes are indicative of a very low-mass brown dwarf with a circumstellar disk. In a near-infrared spectrum of this source obtained with the Gemini Near-Infrared Spectrograph, the presence of strong steam absorption confirms its late-type nature (>=M9.5) while the shapes of the H- and K-band continua and the strengths of the Na I and K I lines demonstrate that it is a young, pre-main-sequence object rather than a field dwarf. A comparison of the bolometric luminosity of Cha 110913-773444 to the luminosities predicted by the evolutionary models of Chabrier and Baraffe and Burrows and coworkers indicates a mass of 8+7/-3 M_Jup, placing it fully within the mass range observed for extrasolar planetary companions (M<=15 M_Jup). The spectral energy distribution of this object exhibits mid-infrared excess emission at >5 um, which we have successfully modeled in terms of an irradiated viscous accretion disk with M'<=10e-12 M_sun/year. Cha 110913-773444 is now the least massive brown dwarf observed to have a circumstellar disk, and indeed is one of the least massive free-floating objects found to date. These results demonstrate that the raw materials for planet formation exist around free-floating planetary-mass bodies.

We present new Spitzer Space Telescope observations of the young cluster NGC2264. Observations at 24 micron with the Multiband Imaging Photometer has enabled us to identify the most highly embedded and youngest objects in NGC2264. This letter reports on one particular region of NGC2264 where bright 24 micron sources are spatially configured in curious linear structures with quasi-uniform separations. The majority of these sources (~60% are found to be protostellar in nature with Class I spectral energy distributions. Comparison of their spatial distribution with sub-millimeter data from Wolf-Chase (2003) and millimeter data from Peretto et al. (2005) shows a close correlation between the dust filaments and the linear spatial configurations of the protostars, indicating that star formation is occurring primarily within dense dusty filaments. Finally, the quasi-uniform separations of the protostars are found to be comparable in magnitude to the expected Jeans length suggesting thermal fragmentation of the dense filamentary material.

We present IRAC 3.6, 4.5, 5.8 and 8 micron photometry for the 17 A, K and M type members of the Eta Chameleontis association. These data show infrared excesses toward six of the 15 K and M stars, indicating the presence of circumstellar disks around 40% of the stars with masses of 0.1-1 solar mass. The two A-stars show no infrared excesses. The excess emission around one of the stars is comparable to the median excess for classical T Tauri stars in the Taurus association; the remaining five show comparatively weak excess emission. Taking into account published Halpha spectroscopy that shows that five of the six stars are accreting, we argue that the disks with weak mid-infrared excesses are disks in which the inner disks have been largely depleted of small grains by grain growth, or, in one case, the small grains have settled to the midplane. This suggests that Eta Cha has a much higher fraction of disks caught in the act of transitioning into optically thin disks than that measured in younger clusters and associations.

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.

We are conducting a survey of several regions of high-mass star formation to assess their content and structure. The observations include Spitzer observations, ground-based optical and near-IR imaging surveys, and optical and IR spectra of objects and locations in the molecular clouds. The goal of the survey is to gain a better understanding of the processes involved in high mass star formation by determining the characteristics of the stars detected in these regions and investigating the properties of the interstellar medium (ISM) environment in which these stars form. In this contribution, we present results on the identification and spatial analysis of young stars in three clusters, W5/AFGL 4029, S255, and S235. First we show how the IRAC data are used to roughly segregate young stars according to their mid-infrared colors, into two groups corresponding the SED Class I and Class II young stellar objects. Then using the IRAC data in combination with 2MASS, we show how more young stars can be identified. Finally, we examine the spatial distributions of young stars in these clusters and find a range of morphologies and of peak surface densities.

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.

The Infrared Array Camera (IRAC) on the Spitzer Space Telescope is absolutely calibrated by comparing photometry on a set of A stars near the north ecliptic pole to predictions based on ground-based observations and a stellar atmosphere model. The brightness of point sources is calibrated to an accuracy of 3%, relative to models for A star stellar atmospheres, for observations performed and analyzed in the same manner as the calibration stars. This includes corrections for location of the star in the array and the location if the centroid within the peak pixel. Long-term stability of the IRAC photometry was measured by monitoring the brightness of A dwarfs and K giants (near the north ecliptic pole) observed several times per month; the photometry is stable to 1.5% (rms) over a year. Intermediate-time-scale stability of the IRAC photometry was measured by monitoring at least one secondary calibrator (near the ecliptic plane) every 12 hr while IRAC is in nominal operations; the intermediate-term photometry is stable with a 1% dispersion (rms). One of the secondary calibrators was found to have significantly time-variable (5%) mid-infrared emission, with period (7.4 days) matching the optical light curve; it is possibly a Cepheid variable.

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 infrared photometry obtained with the IRAC camera on the Spitzer Space Telescope of a sample of 82 pre-main sequence stars and brown dwarfs in the Taurus star-forming region. We find a clear separation in some IRAC color-color diagrams between objects with and without disks. A few ``transition'' objects are noted, which correspond to systems in which the inner disk has been evacuated of small dust. Separating pure disk systems from objects with remnant protostellar envelopes is more difficult at IRAC wavelengths, especially for objects with infall at low rates and large angular momenta. Our results generally confirm the IRAC color classification scheme used in previous papers by Allen et al. and Megeath et al. to distinguish between protostars, T Tauri stars with disks, and young stars without (inner) disks. The observed IRAC colors are in good agreement with recent improved disk models, and in general accord with models for protostellar envelopes derived from analyzing a larger wavelength region. We also comment on a few Taurus objects of special interest. Our results should be useful for interpreting IRAC results in other, less well-studied star-forming regions.

We present high spatial resolution mid- and far-infrared images of the Vega debris disk obtained with the Multiband Imaging Photometer for Spitzer (MIPS). The disk is well resolved and its angular size is much larger than found previously. The radius of the disk is at least 43" (330 AU), 70"(543 AU), and 105" (815 AU) in extent at 24, 70 and 160 um, respectively. The disk images are circular, smooth and without clumpiness at all three wavelengths. The radial surface brightness profiles imply an inner boundary at a radius of 11"+/-2" (86 AU). Assuming an amalgam of amorphous silicate and carbonaceous grains, the disk can be modeled as an axially symmetric and geometrically thin disk, viewed face-on, with the surface particle number density following an r^-1 power law. The disk radiometric properties are consistent with a range of models using grains of sizes ~1 to ~50 um. We find that a ring, containing grains larger than 180 um and at radii of 86-200 AU from the star, can reproduce the observed 850 um flux, while its emission does not violate the observed MIPS profiles. This ring could be associated with a population of larger asteroidal bodies analogous to our own Kuiper Belt. Cascades of collisions starting with encounters amongthese large bodies in the ring produce the small debris that is blown outward by radiation pressure to much larger distances where we detect its thermal emission. The dust production rate is >~10^15 g/s based on the MIPS results. This rate would require a very massive asteroidal reservoir for the dust to be produced in a steady state throughout Vega's life. Instead, we suggest that the disk we imaged is ephemeral and that we are witnessing the aftermath of a large and relatively recent collisional event, and subsequent collisional cascade.

Molecular line observations of NH3 (J,K)=(1,1), (2,2) and CH3OH at 24.93GHz taken with the Australian Telescope Compact Array (ATCA) toward the massive twin cores NGC6334 I & I(N) reveal significant variations in the line emission between the two massive cores. The UCHII region/hot core NGC6334 I exhibits strong thermal NH3 and CH3OH emission adjacent to the UCHII region and coincident with two mm continuum peaks observed by Hunter et al. (in prep.). In contrast, we find neither compact NH3 nor thermal CH3OH line emission toward NGC6334 I(N). There, the NH3 emission is distributed over a broad region (>1') without a clear peak, and we find Class I CH3OH maser emission with peak brightness temperatures up to 7000K. The maser emission peaks appear to be spatially associated with the interfaces between the molecular outflows and the ambient dense gas. Peak NH3(1,1) line brightness temperatures >= 70K in both regions indicate gas temperatures of the same order. NH3 emission is also detected toward the outflow in NGC6334 I resulting in an estimated rotational temperature of Trot~19K. Furthermore, we observe CH3OH and NH3 absorption toward the UCHII region, the velocity structure is consistent with expanding molecular gas around the UCHII region. Thermal and kinematic effects possibly imposed from the UCHII region on the molecular core are also discussed.

We present Spitzer Space Telescope infrared photometric time series of the transiting extrasolar planet system TrES-1. The data span a predicted time of secondary eclipse, corresponding to the passage of the planet behind the star. In both bands of our observations, we detect a flux decrement with a timing, amplitude, and duration as predicted by published parameters of the system. This signal represents the first direct detection of (i.e. the observation of photons emitted by) a planet orbiting another star. The observed eclipse depths (in units of relative flux) are 0.00066 +/- 0.00013 at 4.5um and 0.00225 +/- 0.00036 at 8.0um. These estimates provide the first observational constraints on models of the thermal emission of hot Jupiters. Assuming that the planet emits as a blackbody, we estimate an effective temperature of T_p=1060 +/- 50 K. Under the additional assumptions that the planet is in thermal equilibrium with the radiation from the star and emits isotropically, we find a Bond albedo of A = 0.31 +/- 0.14. This would imply that the planet absorbs the majority of stellar radiation incident upon it, a conclusion of significant impact to atmospheric models of these objects. We compare our data to a previously-published model of the planetary thermal emission, which predicts prominent spectral features in our observational bands due to water and carbon monoxide. Based on the time of secondary eclipse, we present an upper limit on the orbital eccentricity that is sufficiently small that we conclude that tidal dissipation is unlikely to provide a significant source of energy interior to the planet.(abridged)

We present Hubble Space Telescope NICMOS imaging of W3 IRS 5, a binary high-mass protostar. In addition to the two protostars, NICMOS images taken in the F222M and F160W filters show three new 2.22 micron sources with very red colors; these sources fall within a region 5600 AU in diameter, and are coincident with a 100 solar mass dense molecular clump. Two additional point sources are found within 0.4'' (800 AU) of one of the high-mass protostars; these may be stellar companions or unresolved emission knots from an outflow. We propose that these sources constitute a nascent Trapezium system in the center of the W3 IRS 5 cluster containing as many as five proto OB stars. This would be the first identification of a Trapezium still deeply embedded in its natal gas.

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 present Spitzer Space Telescope observations of V1647 Ori, the outbursting source lighting McNeil's nebula, taken near the optical peak of the outburst in early March 2004. The source is easily detected in all Spitzer imaging bands from 3.6 - 70 microns. The fluxes at all wavelengths are roughly a factor of 15 brighter than pre-outburst levels; we measure a bolometric luminosity of 44 Lsun. We posit that this event is due to an increase in the accretion luminosity of the source. Simple models of an accretion disk plus tenuous envelope can qualitatively explain the observed pre- and post-outburst spectral energy distributions. The accretion activity implied by our results indicates that the outburst may be intermediate between FUor and EXor-type events. We also report the discovery of a previously unknown mid-infrared counterpart to the nearby Herbig-Haro object HH 22.

We present spectroscopy of two candidate substellar members of the Chamaeleon I star-forming region. The candidates, which were identified photometrically by Oasa, Tamura, & Sugitani, have been observed at 1-2.5 micron during commissioning of the Gemini Near-Infrared Spectrograph. The late-type nature of one of the candidates, OTS 44, is confirmed through the detection of strong steam absorption bands. The other object, OTS 7, exhibits no late-type features and is likely a background star or galaxy. The gravity-sensitive shape of the H- and K-band continua demonstrate that OTS 44 is a young, pre-main-sequence object rather than a field dwarf. We measure a spectral type of M9.5 for OTS 44 based on a comparison of its spectrum to data for optically-classified young late-type objects. Because OTS 44 is the coolest and faintest object with confirmed membership in Chamaeleon I, it is very likely the least massive known member of the cluster. By comparing the position of OTS 44 on the H-R diagram to the evolutionary models of Chabrier & Baraffe, we infer a mass of ~0.015 M_sun. Although this estimate is uncertain by at least a factor of two, OTS 44 is nevertheless one of the least massive free-floating brown dwarfs confirmed spectroscopically to date.

We analyze multi-band photometry of a subsample of low mass stars in the associations Ori OB1a and 1b discovered during the CIDA Orion Variability Survey, which have ages of 7 - 10 Myr and 3 - 5 Myr, respectively. We obtained UBVRcIc photometry at Mt. Hopkins for 6 Classical T Tauri stars (CTTS) and 26 Weak T Tauri stars (WTTS) in Ori OB1a, and for 21 CTTS and 2 WTTS in Ori OB1b. We also obtained L band photometry for 14 CTTS at Mt. Hopkins, and 10um and 18um photometry with OSCIR at Gemini for 6 CTTS; of these, all 6 were detected at 10um while only one was detected at 18um. We estimate mass accretion rates from the excess luminosity at U, and find that they are consistent with determinations for a number of other associations, with or without high mass star formation. The observed decrease of mass accretion rate with age is qualitatively consistent with predictions of viscous evolution of accretion disks. We find an overall decrease of disk emission from Taurus to Ori OB1b to Ori OB1a. This decrease implies that significant grain growth and settling towards the midplane has taken place in the inner disks of Ori OB1. We compare the SED of the star detected at both 10um and 18um with disk models for similar stellar and accretion parameters. We find that the low <= 18 um fluxes of this Ori OB1b star cannot be due to the smaller disk radius expected from viscous evolution in the presence of the FUV radiation fields from the OB stars in the association. Instead, we find that the disk of this star is essentially a flat disk, with little if any flaring, indicating a a significant degree of dust settling towards the midplane, as expected from dust evolution in 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 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.

High-resolution images of the hypercompact HII regions (HCHII) in W3 IRS 5 taken with the Very Large Array (VLA) at 1.3 and 0.7 cm are presented. Four HCHII regions were detected with sufficient signal-to-noise ratios to allow the determination of relevant parameters such as source position, size and flux density. The sources are slightly extended in our ~0.2 arcsecond beams; the deconvolved radii are less than 240 AU. A comparison of our data with VLA images taken at epoch 1989.1 shows proper motions for sources IRS 5a and IRS 5f. Between 1989.1 and 2002.5, we find a proper motion of 210 mas at a position angle of 12 deg for IRS 5f and a proper motion of 190 mas at a position angle of 50 deg for IRS 5a. At the assumed distance to W3 IRS 5, 1.83 +/- 0.14 kpc, these offsets translate to proper motions of ~135 km/s and ~122 km/s$ respectively. These sources are either shock ionized gas in an outflow or ionized gas ejected from high mass stars. We find no change in the positions of IRS 5d1/d2 and IRS 5b; and we show through a comparison with archival NICMOS 2.2 micron images that these two radio sources coincide with the infrared double constituting W3 IRS 5. These sources contain B or perhaps O stars. The flux densities of the four sources have changed compared to the epoch 1989.1 results. In our epoch 2002.5 data, none of the spectral indicies obtained from flux densities at 1.3 and 0.7 cm are consistent with optically thin free-free emission; IRS 5d1/d2 shows the largest increase in flux density from 1.3 cm to 0.7 cm. This may be an indication of free-free optical depth within an ionized wind, a photoevaporating disk, or an accretion flow. It is less likely that this increase is caused by dust emission at 0.7 cm.

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.

Element-by-element we have combined the optical components in the three 2MASS cameras, and incorporated detector quantum efficiency curves and site-specific atmospheric transmissions, to create three relative spectral response curves (RSRs). We provide absolute 2MASS attributes associated with "zero magnitude" in the JHKs bands so that these RSRs may be used for synthetic photometry. The RSRs tie 2MASS to the Cohen-Walker-Witteborn framework of absolute photometry and spectra for the purpose of using 2MASS data to support the development of absolute calibrators for IRAC and pairwise cross-calibrators between all three SIRTF instruments. We examine the robustness of these RSRs to changes in water vapor within a night. We compare the observed 2MASS magnitudes of 33 stars (converted from the precision optical calibrators of Landolt and Carter-Meadows into absolute infrared calibrators from 1.2-35 um with our predictions, thereby deriving 2MASS "zero point offsets" from the ensemble. These offsets are the final ingredients essential to merge 2MASS $JHK_{s}$ data with other absolutely calibrated bands and stellar spectra, and to support the creation of faint calibration stars for SIRTF.

We describe the technique used to develop on-orbit calibrators for IRAC and validate the method for either K0-M0III or A0-A5V stars. For SIRTF, the approach uses all available optical, NIR, and MIR photometry, and yields absolute spectra from UV to MIR. One set of stars is from Landolt's network of optical standards, the other from Carter-Meadows IR standards. Traceability to the Cohen-Walker- Witteborn framework of absolute photometry and spectra is assured. The method is based on using either "supertemplates" to represent the intrinsic shapes of the spectra of K0-M0IIIs from 1150A to 35 um, or Kurucz synthetic spectra for A0-5V stars. Each supertemplate/model is reddened according to a star's extinction and normalized by characterized optical photometry. This paper tests our ability to predict NIR (JHK) magnitudes from supertemplates or models constrained by Hippa- rcos/Tycho or precision ground-based optical data. We offer absolute calibrated spectra of 33 optical standards to demonstrate the viability of this technique for a set of IR calibrators 100-1000 times fainter than we have previously publ- ished. We calculate the absolute uncertainties associated with predicting IRAC mags for the faintest cool giant and hot dwarf in this new set of calibrators.

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.

In order to obtain a census of the pre-stellar and star-forming molecular cores, we have begun an unbiased survey in CS and N2H+ of the L1630 and L1641 molecular clouds. The use of these two molecular species enables us to quantify and disentangle the effects of depletion often seen in CS observations of dense cores. The spectral line data will provide essential kinematical information not present in similar studies of the sub-millimeter dust-continuum, enabling us to examine the overall core to core velocity dispersion and study the effects of infall and outflows around known sub-mm and infra-red sources. Here we present our initial observations of part of L1630, taken during the commissioning phase of the FCRAO On-the-Fly Mapping system in January 2002.

The connection between high-mass stars and young stellar clusters has been well established by near-IR observations showing young massive stars in HII regions surrounded by clusters of lower mass stars. Mm-wave observations show that these clusters form in approx. 1000 Msol dense cores. Interestingly, the mm-wave observations of very active cluster-forming regions also reveal the presence of seemingly quiescent dense massive cores, which contain few signs of ongoing star formation. In the W3 Main region we find an approx. 1000 Msol mass dense core which contains no mid-IR or cm-continuum sources. Near-IR imaging has failed to detect any embedded sources. Yet, we suggest that this region is at an early stage of cluster formation. This is given support by comparative observations of the NGC 6334 GMC. Here, we find another quiescent massive 1000 Msol core without mid-IR sources, HII regions, or a detectable near-IR cluster. Motivated by the presence of CH3OH and H2O masers in this IR-quiet core, we have searched for and detected a bipolar outflow driven by a young massive star deeply embedded in the core. The presence of an embedded high-mass star strengthens our arguments that IR-quiet cores are at an early stage of cluster formation. We propose that masers and outflows from very deeply embedded (and consequently IR-quiet) high-mass stars are among the first incipient signposts of cluster formation in massive cores. Accordingly, we have undertaken a SEST survey of 46 CH3OH maser clusters which are not coincident with IRAS sources or HII regions. From this sample, we detected outflows in transitions of SiO toward 15 of the maser clusters. These sources have now been observed at 1200 micron with SIMBA, the 37-multichannel bolometer at SEST. Each one is associated with a bright 1 mm point source.

Massive OB stars play an important role in the evolution of molecular clouds and star forming regions. The OB stars both photo--ionize molecular gas as well as sweep up and compress interstellar gas through winds, ionization fronts, and supernovae. In this contribution, we examine the effect massive stars have on the formation of young stellar clusters. We first discuss the processes by which OB stars destroy cluster--forming molecular cores, and hence terminate star formation. We overview observational evidence that OB stars forming in young stellar clusters destroy their parental cores on a timescale of 0.1 Myr, and we discuss some of the implications of this result. We then summarize extensive observations of the NGC 281 molecular cloud complex, and present evidence that two types of triggered star formation are occurring in this complex. Our goal is to underscore the impact massive stars have on cluster formation over distances ranging from 0.1 pc to 300 pc.

We find strong evidence for outflows originating in the dense molecular core NGC 6334 I(North): a 1000 Msol molecular core distinguished by its lack of HII regions and mid-IR emission. New observations were obtained of the SiO 2-1 and 5-4 lines with the SEST 15-m telescope and the H2 (1-0) S(1) line with the ESO 2.2-m telescope. The line profiles of the SiO transitions show broad wings extending from -50 to 40 km/s, and spatial maps of the line wing emission exhibit a bipolar morphology with the peaks of the red and blue wing separated by 30". The estimated mass loss rate of the outflow is comparable to those for young intermediate to high-mass stars. The near-IR images show eight knots of H2 emission. Five of the knots form a linear chain which is displaced from the axis of the SiO outflow; these knots may trace shock excited gas along the path of a second, highly collimated outflow. We propose that I(N) is a rare example of a molecular core in an early stage of cluster formation.