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 adopt the deep learning method CASI (Convolutional Approach to Shell Identification) and extend it to 3D (CASI-3D) to identify signatures of stellar feedback in molecular line spectra, such as 13CO. We adopt magneto-hydrodynamics simulations that study the impact of stellar winds in a turbulent molecular cloud as an input to generate synthetic observations. We apply the 3D radiation transfer code radmc-3d to model 13CO (J=1-0) line emission from the simulated clouds. We train two CASI-3d models: ME1 predicts only the position of feedback, while MF predicts the fraction of the mass coming from feedback in each voxel. We adopt 75% of the synthetic observations as the training set and assess the accuracy of the two models with the remaining data. We demonstrate that model ME1 identifies bubbles in simulated data with 95% accuracy, and model MF predicts the bubble mass within 4% of the true value. We use bubbles previously visually identified in Taurus in 13CO to validate the models and show both perform well on the highest confidence bubbles. We apply our two models on the full 98 deg2 FCRAO 13CO survey of the Taurus cloud. Models ME1 and MF predict feedback gas mass of 2894 M and 302 M, respectively. When including a correction factor for missing energy due to the limited velocity range of the 13CO data cube, model ME1 predicts feedback kinetic energies of 4.0*1e46 ergs and 1.5*1e47 ergs with/without subtracting the cloud velocity gradient. Model MF predicts feedback kinetic energy of 9.6*1e45 ergs and 2.8*1e46 ergs with/without subtracting the cloud velocity gradient. Model ME1 predicts bubble locations and properties consistent with previous visually identified bubbles. However, model MF demonstrates that feedback properties computed based on visual identifications are significantly overestimated due to line of sight confusion and contamination from background and foreground gas.

In this paper we undertake a study of the 21 square degree SMOG field, a Spitzer cryogenic mission Legacy program to map a region of the outer Milky Way towards the Perseus and Outer spiral arms with the IRAC and MIPS instruments. We identify 4648 YSOs across the field. Using the DBSCAN method we identify 68 clusters or aggregations of YSOs in the region, having 8 or more members. We identify 1197 class Is, 2632 class IIs, 819 class IIIs, of which 45 are candidate transition disk objects, utilizing the MIPS 24 photometry. The ratio of YSOs identified as members of clusters was 2872/4648, or 62%. The ratios of class I to class II YSOs in the clusters are broadly consistent with those found in the inner galactic and nearby Gould's Belt young star formation regions. The clustering properties indicate that the protostars may be more tightly bound to their natal sites than the class IIs, and the class IIIs are generally widely distributed. We further perform an analysis of the WISE data of the SMOG field to determine how the lower resolution and sensitivity of WISE affects the identification of YSOs as compared to Spitzer: we identify 931 YSOs using combined WISE and 2MASS photometry, 931/4648 or 20% of the total number identified with Spitzer. Performing the same clustering analysis finds 31 clusters which reliably trace the larger associations identified with the Spitzer data. Twelve of the clusters identified have previously measured distances from the WISE HII survey. SEDFitter modeling of these YSOs is reported, leading to an estimation of the IMF in the aggregate of these clusters which approximates that found in the inner galaxy, implying that the processes behind stellar mass distribution during star formation are not widely affected by the lower density and metallicity of the outer galaxy.

We utilize techniques from deep learning to identify signatures of stellar feedback in simulated molecular clouds. Specifically, we implement a deep neural network with an architecture similar to U-Net and apply it to the problem of identifying wind-driven shells and bubbles using data from magneto-hydrodynamic simulations of turbulent molecular clouds with embedded stellar sources. The network is applied to two tasks, dense regression and segmentation, on two varieties of data, simulated density and synthetic 12 CO observations. Our Convolutional Approach for Shell Identification (CASI) is able to obtain a true positive rate greater than 90\%, while maintaining a false positive rate of 1\%, on two segmentation tasks and also performs well on related regression tasks. The source code for CASI is available on GitLab.

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.

We present multi-epoch optical and near-infrared (NIR) photometry and spectroscopy of the spectroscopic binary T Tauri star DQ Tau. The photometric monitoring, obtained using SMARTS ANDICAM, recovers the previously-seen correlation between optical flux and the 15.8-day binary orbital period, with blue flux peaks occurring close to most observed periastron passages. For the first time, we find an even more consistent correlation between orbital period and NIR brightness and color. The onset of pulse events in the NIR on average precedes those in the optical by a few days, with the rise usually starting near apastron orbital phase. We further obtained five epochs of spectroscopy using IRTF SpeX, with a wavelength range of 0.8 to 5 microns, and derived spectra of the infrared excess emission. The shape and strength of the excess varies with time, with cooler and weaker characteristic dust emission (T ~ 1100-1300 K) over most of the binary orbit, and stronger/warmer dust emission (T ~ 1600 K, indicative of dust sublimation) just before periastron passage. We suggest our results are broadly consistent with predictions of simulations of disk structure and accretion flows around close binaries, with the varying dust emission possibly tracing the evolution of accretion streams falling inwards through a circumbinary disk cavity and feeding the accretion pulses traced by the optical photometry and NIR emission lines. However, our results also show more complicated behavior that is not fully explained by this simple picture, and will require further observations and modeling to fully interpret.

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.

Understanding the collapse of clouds and the formation of protoplanetary disks is essential to understanding the formation of stars and planets. Infall and accretion, the mass-aggregation processes that occur at envelope and disk scales, drive the dynamical evolution of protostars. While the observations of protostars at different stages constrain their evolutionary tracks, the impact of variability due to accretion bursts on dynamical and chemical evolution of the source is largely unknown. The lasting effects on protostellar envelopes and disks are tracked through multi-wavelength and time domain observational campaigns, requiring deep X-ray, infrared, and radio imaging and spectroscopy, at a sufficient level of spatial detail to distinguish contributions from the various substructures (i.e., envelope from disk from star from outflow). Protostellar models derived from these campaigns will illuminate the initial chemical state of protoplanetary disks during the epoch of giant planet formation. Insight from individual star formation in the Milky Way is also necessary to understand star formation rates in extragalactic sources. This cannot be achieved with ground-based observatories and is not covered by currently approved instrumentation. Requirements: High (v < 10 km/s for survey; v < 1 km/s for followup) spectral resolution capabilities with relatively rapid response times in the IR (3-500 um), X-ray (0.1-10 keV), and radio (cm) are critical to follow the course of accretion and outflow during an outburst. Complementary, AU-scale radio observations are needed to probe the disk accretion zone, and 10 AU-scale to probe chemical and kinematic structures of the disk-forming regions, and track changes in the dust, ice, and gas within protostellar envelopes.

Context: The study of dust emission at millimeter wavelengths is important to shed light on the dust properties and physical structure of pre-stellar cores, the initial conditions in the process of star and planet formation. Aims: Using two new continuum facilities, AzTEC at the LMT and MUSTANG-2 at the GBO, we aim to detect changes in the optical properties of dust grains as a function of radius for the well-known pre-stellar core L1544. Methods: We determine the emission profiles at 1.1 and 3.3 mm and examine whether they can be reproduced in terms of the current best physical models for L1544. We also make use of various tools to determine the radial distributions of the density, temperature, and the dust opacity in a self-consistent manner. Results: We find that our observations cannot be reproduced without invoking opacity variations. With the new data, new temperature and density profiles, as well as opacity variations across the core, have been derived. The opacity changes are consistent with the expected variations between uncoagulated bare grains, toward the outer regions of the core, and grains with thick ice mantles, toward the core center. A simple analytical grain growth model predicts the presence of grains of ~3-4 um within the central 2000 au for the new density profile.

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

We present a 1.1~mm census of dense cores in the Mon~R2 Giant Molecular Cloud with the AzTEC instrument on the Large Millimeter Telescope (LMT). We detect 295 cores (209 starless, and 86 with protostars) in a two square degree shallow survey. We also carry out a deep follow-up survey of 9 regions with low to intermediate ($3<A_V<7$) gas column densities and detect 60 new cores in the deeper survey which allows us to derive a completeness limit. After performing corrections for low signal-to-noise cores, we find a median core mass of $\sim 2.1 \text{M}_{\odot}$ and a median size of $ 0.08$~pc. $46\%$ of the cores (141) have masses exceeding the local Bonor-Ebert mass for cores with T=12K, suggesting that in the absence of supporting non-thermal pressure, these regions are unstable to gravitational collapse. We present the core mass function (CMF) for various subdivisions of the core sample. We find that cores with masses $>$10~$M_{\odot}$ are exclusively found in regions with high core number densities and that the CMF of the starless cores has an excess of low-mass cores ($<$5~$M_{\odot}$) compared to the CMF of protostellar cores. We report a power law correlation of index $1.99 \pm 0.03$ between local core mass density and gas column density (as traced by Herschel) over a wide range of size scales (0.3-5~pc). This power law is consistent with that predicted for thermal fragmentation of a self-gravitating sheet.

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 present a Chandra study of the deeply embedded Serpens South star-forming region, examining cluster structure and disk properties at the earliest stages. In total, 152 X-ray sources are detected. Combined with Spitzer and 2MASS photometry, 66 X-ray sources are reliably matched to an IR counterpart. We identify 21 class I, 6 flat spectrum, 16 class II, and 18 class III young stars; 5 were unclassified. Eighteen sources were variable in X-rays, 8 exhibiting flare-like emission, and one periodic source. The cluster X-ray luminosity distance was estimated, the best match was to the nearer distance of 260pc for the front of the Aquila Rift complex. The $N_{H}$ vs. $A_{K}$ ratio is found to be $\sim$0.68x10$^{22}$, similar to that measured in other young low mass regions, but lower than that measured in the ISM and high mass clusters ($\sim$1.6-2x10$^{22}$). We find the spatial distribution closely follows that of the dense filament from which the stars have formed, with the class II population still strongly associated with the filament. There are four sub-clusters in the field, with three forming knots in the filament, and a fourth to the west, which may not be associated but may be contributing to the distributed class III population. A high percentage of diskless class IIIs (upper limit 30% of classified X-ray sources) in such a young cluster could indicate that processing of disks is influenced by the cluster environment and is not solely time-scale dependent.

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

Sensitive, imaging observations of the 1.1 mm dust continuum emission from a 1 deg^2 area collected with the AzTEC bolometer camera on the Large Millimeter Telescope are presented. A catalog of 1545 compact sources is constructed based on a Wiener-optimization filter. These sources are linked to larger clump structures identified in the Bolocam Galactic Plane Survey. Hydrogen column densities are calculated for all sources and mass and mean volume densities are derived for the subset of sources for which kinematic distances can be assigned. The AzTEC sources are localized, high density peaks within the massive clumps of molecular clouds and comprise 5-15% of the clump mass. We examine the role of the gravitational instability in generating these fragments by comparing the mass of embedded AzTEC sources to the Jeans' mass of the parent BGPS object. For sources with distances less than 6 kpc the fragment masses are comparable to the clump Jeans' mass, despite having isothermal Mach numbers between 1.6 and 7.2. AzTEC sources linked to ultra-compact HII regions have mass surface densities greater than the critical value implied by the mass-size relationship of infrared dark clouds with high mass star formation while AzTEC sources associated with Class II methanol masers have mass surface densities greater than 0.7 g cm^{-2} that approaches the proposed threshold required to form massive stars.

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

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 new deep UBVRI images and high-resolution multi-object optical spectroscopy of the young (~ 6 - 10 Myr old), relatively nearby (800 pc) open cluster IC 2395. We identify nearly 300 cluster members and use the photometry to estimate their spectral types, which extend from early B to middle M. We also present an infrared imaging survey of the central region using the IRAC and MIPS instruments on board the Spitzer Space Telescope, covering the wavelength range from 3.6 to 24 microns. Our infrared observations allow us to detect dust in circumstellar disks originating over a typical range of radii ~ 0.1 to ~ 10AU from the central star. We identify 18 Class II, 8 transitional disk, and 23 debris disk candidates, respectively 6.5%, 2.9%, and 8.3% of the cluster members with appropriate data. We apply the same criteria for transitional disk identification to 19 other stellar clusters and associations spanning ages from ~ 1 to ~ 18 Myr. We find that the number of disks in the transitional phase as a fraction of the total with strong 24 micron excesses ([8] - [24] > 1.5) increases from 8.4 +\- 1.3% at ~ 3 Myr to 46 +\- 5% at ~ 10 Myr. Alternative definitions of transitional disks will yield different percentages but should show the same trend.

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

Theoretical models and spectroscopic observations of newborn stars suggest that protoplantary disks have an inner "wall" at a distance set by the disk interaction with the star. Around T Tauri stars, the size of this disk hole is expected to be on a 0.1-AU scale that is unresolved by current adaptive optics imaging, though some model-dependent constraints have been obtained by near-infrared interferometry. Here we report the first measurement of the inner disk wall around a solar-mass young stellar object, YLW 16B in the {\rho} Ophiuchi star-forming region, by detecting the light travel time of the variable radiation from the stellar surface to the disk. Consistent time lags were detected on two nights, when the time series in H (1.6 {\mu}m) and K (2.2 {\mu}m) bands were synchronized while the 4.5 {\mu}m emission lagged by 74.5 +/- 3.2 seconds. Considering the nearly edge-on geometry of the disk, the inner rim should be 0.084 AU from the protostar on average, with an error of order 0.01 AU. This size is likely larger than the range of magnetospheric truncations, and consistent with an optically and geometrically thick disk front at the dust sublimation radius at ~1500 K. The widths of the cross-correlation functions between the data in different wavebands place possible new constraints on the geometry of the disk.

We examine the light and color evolution of the T Tauri binary KH 15D through photometry obtained at wavelengths between 0.55 and 8.0 $\mu$m. The data were collected with ANDICAM on the 1.3 m SMARTS telescope at Cerro-Tololo Inter-American Observatory and with IRAC on the Spitzer Space Telescope. We show that the system's circumbinary ring, which acts as a screen that covers and uncovers different portions of the binary orbit as the ring precesses, has reached an orientation where the brighter component (star B) fully or nearly fully emerges during each orbital cycle. The fainter component (star A) remains fully occulted by the screen at all phases. The leading and trailing edges of the screen move across the sky at the same rate of $\sim$15 meters per second, consistent with expectation for a ring with a radius and width of $\sim$4 AU and a precession period of $\sim$6500 years. Light and color variations continue to indicate that the screen is sharp edged and opaque at \emph{VRIJH} wavelengths. However, we find an increasing transparency of the ring edge at 2.2, 3.6, and 4.5 $\mu$m. Reddening seen at the beginning of the eclipse that occurred during the CSI 2264 campaign particularly suggests selective extinction by a population of large dust grains. Meanwhile, the gradual bluing observed while star B is setting is indicative of forward scattering effects at the edge of the ring. The SED of the system at its bright phase shows no evidence of infrared excess emission that can be attributed to radiation from the ring or other dust component out to 8 microns.

Newborn stars form within the localized, high density regions of molecular clouds. The sequence and rate at which stars form in dense clumps and the dependence on local and global environments are key factors in developing descriptions of stellar production in galaxies. We seek to observationally constrain the rate and latency of star formation in dense massive clumps that are distributed throughout the Galaxy and to compare these results to proposed prescriptions for stellar production. A sample of 24 micron-based Class~I protostars are linked to dust clumps that are embedded within molecular clouds selected from the APEX Telescope Large Area Survey of the Galaxy. We determine the fraction of star-forming clumps, f*, that imposes a constraint on the latency of star formation in units of a clump's lifetime. Protostellar masses are estimated from models of circumstellar environments of young stellar objects from which star formation rates are derived. Physical properties of the clumps are calculated from 870 micron dust continuum emission and NH_3 line emission. Linear correlations are identified between the star formation rate surface density, Sigma_{SFR}, and the quantities Sigma_{H2}/tau_{ff} and Sigma_{H2}/tau_{cross}, suggesting that star formation is regulated at the local scales of molecular clouds. The measured fraction of star forming clumps is 23%. Accounting for star formation within clumps that are excluded from our sample due to 24 micron saturation, this fraction can be as high as 31%. Dense, massive clumps form primarily low mass (< 1-2 msun) stars with emergent 24 micron fluxes below our sensitivity limit or are incapable of forming any stars for the initial 70% of their lifetimes. The low fraction of star forming clumps in the Galactic center relative to those located in the disk of the Milky Way is verified.

We provide CoRoT and Spitzer light curves, as well as broad-band multi-wavelength photometry and high resolution, multi- and single-epoch spectroscopy for 17 classical T Tauris in NGC 2264 whose CoRoT light curves (LCs) exemplify the "stochastic" LC class as defined in Cody et al. (2014). The most probable physical mechanism to explain the optical variability in this LC class is time-dependent mass accretion onto the stellar photosphere, producing transient hot spots. As evidence in favor of this hypothesis, multi-epoch high resolution spectra for a subset of these stars shows that their veiling levels also vary in time and that this veiling variability is consistent in both amplitude and timescale with the optical LC morphology. Furthermore, the veiling variability is well-correlated with the strength of the HeI 6678A emission line, a feature predicted by models to arise in accretion shocks on or near the stellar photosphere. Stars with accretion burst LC morphology (Stauffer et al. 2014) are also attributed to variable mass accretion. Both the stochastic and accretion burst LCs can be explained by a simple model of randomly occurring flux bursts, with the stochastic LC class having a higher frequency of lower amplitude events. Based on their UV excesses, veiling, and mean Ha equivalent widths, members of the stochastic LC class have only moderate time-averaged mass accretion rates. The most common feature of their Ha profiles is for them to exhibit blue-shifted absorption features, most likely originating in a disk wind. The lack of periodic signatures in the LCs suggests that little of the variability is due to long-lived hot spots rotating into or out of our line of sight; instead, the primary driver of the observed photometric variability is likely to be instabilities in the inner disk that lead to variable mass accretion.

We present the initial results of our investigation of the star-forming complex W49, one of the youngest and most luminous massive star forming regions in our Galaxy. We used Spitzer/Infrared Array Camera (IRAC) data to investigate massive star formation with the primary objective to locate a representative set of protostars and the clusters of young stars that are forming around them. We present our source catalog with the mosaics from the IRAC data. In this study we used a combination of IRAC, MIPS, Two Micron All Sky Survey (2MASS) and UKIRT Deep Infrared Sky Survey (UKIDSS) data to identify and classify the Young Stellar Objects (YSOs). We identified 232 Class 0/I YSOs, 907 Class II YSOs, and 74 transition disk candidate objects using color-color and color-magnitude diagrams. In addition, to understand the evolution of star formation in W49 we analysed the distribution of YSOs in the region to identify clusters using a minimal spanning tree method. The fraction of YSOs that belong to clusters with >7 members is found to be 52% for a cut-off distance of 96" and the ratio of Class II/I objects is 2.1. We compared the W49 region to the G305 and G333 star forming regions and concluded that the W49 has the richest population with 7 subclusters of YSOs.

We present deep radio continuum observations of the star-forming core of the Serpens South Infrared Dark Cloud with the Karl G. Jansky Very Large Array (VLA). Observations were conducted in two bands centered at 7.25 GHz (4.14 cm) and 4.75 GHz (6.31 cm) with an rms of 8.5 and 11.1 microJy/beam, respectively. We also use 2MASS, Spitzer and Herschel data to put our radio observations in the context of young stellar populations characterized by near and far infrared observations. Within a 5 arcmin x 5 arcmin region of interest around the central cluster, we detect roughly eighteen radio sources, seven of which we determine are protostellar in nature due to their radio spectral indices and their association with infrared sources. We find evidence for a previously undetected embedded Class 0 protostar and reaffirm Class 0 protostellar classifications determined by previous millimeter wavelength continuum studies. We use our infrared data to derive mid-infrared luminosities for three of our protostellar sources and find relative agreement between the known YSO radio luminosity vs bolometric luminosity correlation. Lastly, we marginally detect an additional six radio sources at the 2-3 sigma level that lie within two arcseconds of infrared YSO candidates, providing motivation for higher sensitivity studies to clarify the nature of these sources and further probe embedded and/or low luminosity YSOs in Serpens South.

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.

LRLL 54361 is an infrared source located in the star forming region IC 348 SW. Remarkably, its infrared luminosity increases by a factor of 10 during roughly one week every 25.34 days. To understand the origin of these remarkable periodic variations, we obtained sensitive 3.3 cm JVLA radio continuum observations of LRLL 54361 and its surroundings in six different epochs: three of them during the IR-on state and three during the IR-off state. The radio source associated with LRLL 54361 remained steady and did not show a correlation with the IR variations. We suggest that the IR is tracing the results of fast (with a timescale of days) pulsed accretion from an unseen binary companion, while the radio traces an ionized outflow with an extent of $\sim$100 AU that smooths out the variability over a period of order a year. The average flux density measured in these 2014 observations, 27$\pm$5 $\mu$Jy, is about a factor of two less than that measured about 1.5 years before, $53\pm$11 $\mu$Jy, suggesting that variability in the radio is present, but over larger timescales than in the IR. We discuss other sources in the field, in particular two infrared/X-ray stars that show rapidly varying gyrosynchrotron emission.

As part of the Young Stellar Object VARiability (YSOVAR) program, we monitored NGC 1333 for ~35 days at 3.6 and 4.5 um using the Spitzer Space Telescope. We report here on the mid-infrared variability of the point sources in the ~10x~20arcmin area centered on 03:29:06, +31:19:30 (J2000). Out of 701 light curves in either channel, we find 78 variables over the YSOVAR campaign. About half of the members are variable. The variable fraction for the most embedded SEDs (Class I, flat) is higher than that for less embedded SEDs (Class II), which is in turn higher than the star-like SEDs (Class III). A few objects have amplitudes (10-90th percentile brightness) in [3.6] or [4.5]>0.2 mag; a more typical amplitude is 0.1-0.15 mag. The largest color change is >0.2 mag. There are 24 periodic objects, with 40% of them being flat SED class. This may mean that the periodic signal is primarily from the disk, not the photosphere, in those cases. We find 9 variables likely to be 'dippers', where texture in the disk occults the central star, and 11 likely to be 'bursters', where accretion instabilities create brightness bursts. There are 39 objects that have significant trends in [3.6]-[4.5] color over the campaign, about evenly divided between redder-when-fainter (consistent with extinction variations) and bluer-when-fainter. About a third of the 17 Class 0 and/or jet-driving sources from the literature are variable over the YSOVAR campaign, and a larger fraction (~half) are variable between the YSOVAR campaign and the cryogenic-era Spitzer observations (6-7 years), perhaps because it takes time for the envelope to respond to changes in the central source. The NGC 1333 brown dwarfs do not stand out from the stellar light curves in any way except there is a much larger fraction of periodic objects (~60% of variable brown dwarfs are periodic, compared to ~30% of the variables overall).

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

We present an IR-monitoring survey with the $Spitzer$ Space Telescope of the star forming region GGD 12-15. Over 1000 objects were monitored including about 350 objects within the central 5 arcminutes which is found to be especially dense in cluster members. The monitoring took place over 38 days and is part of the Young Stellar Object VARiability (YSOVAR) project. The region was also the subject of a contemporaneous 67ks $Chandra$ observation. The field includes 119 previously identified pre-main sequence star candidates. X-rays are detected from 164 objects, 90 of which are identified with cluster members. Overall, we find that about half the objects in the central 5 arcminutes are young stellar objects based on a combination of their spectral energy distribution, IR variability and X-ray emission. Most of the stars with IR excess relative to a photosphere show large amplitude (>0.1 mag) mid-IR variability. There are 39 periodic sources, all but one of these is found to be a cluster member. Almost half of the periodic sources do not show IR excesses. Overall, more than 85% of the Class I, flat spectrum, and Class II sources are found to vary. The amplitude of the variability is larger in more embedded young stellar objects. Most of the Class~I/II objects exhibit redder colors in a fainter state, compatible with time-variable extinction. A few become bluer when fainter, which can be explained with significant changes in the structure of the inner disk. A search for changes in the IR due to X-ray events is carried out, but the low number of flares prevented an analysis of the direct impact of X-ray flares on the IR lightcurves. However, we find that X-ray detected Class II sources have longer timescales for change in the mid-IR than a similar set of non-X-ray detected Class IIs.

We present a time-variability study of young stellar objects in the cluster IRAS 20050+2720, performed at 3.6 and 4.5 micron with the Spitzer Space Telescope; this study is part of the Young Stellar Object VARiability project (YSOVAR). We have collected light curves for 181 cluster members over 40 days. We find a high variability fraction among embedded cluster members of ca. 70%, whereas young stars without a detectable disk display variability less often (in ca. 50% of the cases) and with lower amplitudes. We detect periodic variability for 33 sources with periods primarily in the range of 2-6 days. Practically all embedded periodic sources display additional variability on top of their periodicity. Furthermore, we analyze the slopes of the tracks that our sources span in the color-magnitude diagram (CMD). We find that sources with long variability time scales tend to display CMD slopes that are at least partially influenced by accretion processes, while sources with short variability time scales tend to display extinction-dominated slopes. We find a tentative trend of X-ray detected cluster members to vary on longer time scales than the X-ray undetected members.

We identify nine young stellar objects (YSOs) in the NGC 2264 star-forming region with optical {\em CoRoT} light curves exhibiting short-duration, shallow, periodic flux dips. All of these stars have infrared (IR) excesses that are consistent with their having inner disk walls near the Keplerian co-rotation radius. The repeating photometric dips have FWHM generally less than one day, depths almost always less than 15%, and periods (3<P<11 days) consistent with dust near the Keplerian co-rotation period. The flux dips vary considerably in their depth from epoch to epoch, but usually persist for several weeks and, in two cases, were present in data collected on successive years. For several of these stars, we also measure the photospheric rotation period and find that the rotation and dip periods are the same, as predicted by standard "disk-locking" models. We attribute these flux dips to clumps of material in or near the inner disk wall, passing through our line of sight to the stellar photosphere. In some cases, these dips are also present in simultaneous {\em Spitzer} IRAC light curves at 3.6 and 4.5 microns. We characterize the properties of these dips, and compare the stars with light curves exhibiting this behavior to other classes of YSO in NGC 2264. A number of physical mechanisms could locally increase the dust scale height near the inner disk wall, and we discuss several of those mechanisms; the most plausible mechanisms are either a disk warp due to interaction with the stellar magnetic field or dust entrained in funnel-flow accretion columns arising near the inner disk wall.

We present SCUBA-2 450micron and 850micron observations of the Serpens MWC 297 region, part of the JCMT Gould Belt Survey of nearby star-forming regions. Simulations suggest that radiative feedback influences the star-formation process and we investigate observational evidence for this by constructing temperature maps. Maps are derived from the ratio of SCUBA-2 fluxes and a two component model of the JCMT beam for a fixed dust opacity spectral index of beta = 1.8. Within 40 of the B1.5Ve Herbig star MWC 297, the submillimetre fluxes are contaminated by free-free emission with a spectral index of 1.03+-0.02, consistent with an ultra-compact HII region and polar winds/jets. Contamination accounts for 73+-5 per cent and 82+-4 per cent of peak flux at 450micron and 850micron respectively. The residual thermal disk of the star is almost undetectable at these wavelengths. Young Stellar Objects are confirmed where SCUBA-2 850micron clumps identified by the fellwalker algorithm coincide with Spitzer Gould Belt Survey detections. We identify 23 objects and use Tbol to classify nine YSOs with masses 0.09 to 5.1 Msun. We find two Class 0, one Class 0/I, three Class I and three Class II sources. The mean temperature is 15+-2K for the nine YSOs and 32+-4K for the 14 starless clumps. We observe a starless clump with an abnormally high mean temperature of 46+-2K and conclude that it is radiatively heated by the star MWC 297. Jeans stability provides evidence that radiative heating by the star MWC 297 may be suppressing clump collapse.

In this contribution, we describe the applied methods to construct a 24 micron based point source catalog derived from the image data of the MIPSGAL 24 micron Galactic Plane Survey and the corresponding data products. The high quality catalog product contains 933,818 sources, with a total of 1,353,228 in the full archive catalog. The source tables include positional and photometric informa- tion derived from the 24 micron images, source quality and confusion flags and coun- terpart photometry from matched 2MASS, GLIMPSE, and WISE point sources. Completeness decay data cubes are constructed at 1 arcminute angular resolution that describe the varying background levels over the MIPSGAL field and the ability to extract sources of a given magnitude from this background. The completeness decay cubes are included in the set of data products. We present the results of our efforts to verify the astrometric and photometric calibration of the catalog, and present several analyses of minor anomalies in these measurements to justify adopted mitigation strategies.

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.

The emission from young stellar objects (YSOs) in the mid-IR is dominated by the inner rim of their circumstellar disks. We present an IR-monitoring survey of about 800 objects in the direction of the Lynds 1688 (L1688) star forming region over four visibility windows spanning 1.6 years using the \emph{Spitzer} space telescope in its warm mission phase. Among all lightcurves, 57 sources are cluster members identified based on their spectral-energy distribution and X-ray emission. Almost all cluster members show significant variability. The amplitude of the variability is larger in more embedded YSOs. Ten out of 57 cluster members have periodic variations in the lightcurves with periods typically between three and seven days, but even for those sources, significant variability in addition to the periodic signal can be seen. No period is stable over 1.6 years. Non-periodic lightcurves often still show a preferred timescale of variability which is longer for more embedded sources. About half of all sources exhibit redder colors in a fainter state. This is compatible with time-variable absorption towards the YSO. The other half becomes bluer when fainter. These colors can only be explained with significant changes in the structure of the inner disk. No relation between mid-IR variability and stellar effective temperature or X-ray spectrum is found.

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.

We report Herschel/PACS photometric observations at 70 {\mu}m and 160 {\mu}m of LRLL54361 - a suspected binary protostar that exhibits periodic (P=25.34 days) flux variations at shorter wavelengths (3.6 {\mu}m and 4.5 {\mu}m) thought to be due to pulsed accretion caused by binary motion. The PACS observations show unprecedented flux variation at these far-infrared wavelengths that are well cor- related with the variations at shorter wavelengths. At 70 {\mu}m the object increases its flux by a factor of six while at 160{\mu}m the change is about a factor of two, consistent with the wavelength dependence seen in the far-infrared spectra. The source is marginally resolved at 70 {\mu}m with varying FWHM. Deconvolved images of the sources show elongations exactly matching the outflow cavities traced by the scattered light observations. The spatial variations are anti-correlated with the flux variation indicating that a light echo is responsible for the changes in FWHM. The observed far-infrared flux variability indicates that the disk and en- velope of this source is periodically heated by the accretion pulses of the central source, and suggests that such long wavelength variability in general may provide a reasonable proxy for accretion variations in protostars.

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 present the Coordinated Synoptic Investigation of NGC 2264, a continuous 30-day multi-wavelength photometric monitoring campaign on more than 1000 young cluster members using 16 telescopes. The unprecedented combination of multi-wavelength, high-precision, high-cadence, and long-duration data opens a new window into the time domain behavior of young stellar objects. Here we provide an overview of the observations, focusing on results from Spitzer and CoRoT. The highlight of this work is detailed analysis of 162 classical T Tauri stars for which we can probe optical and mid-infrared flux variations to 1% amplitudes and sub-hour timescales. We present a morphological variability census and then use metrics of periodicity, stochasticity, and symmetry to statistically separate the light curves into seven distinct classes, which we suggest represent different physical processes and geometric effects. We provide distributions of the characteristic timescales and amplitudes, and assess the fractional representation within each class. The largest category (>20%) are optical "dippers" having discrete fading events lasting ~1-5 days. The degree of correlation between the optical and infrared light curves is positive but weak; notably, the independently assigned optical and infrared morphology classes tend to be different for the same object. Assessment of flux variation behavior with respect to (circum)stellar properties reveals correlations of variability parameters with H$\alpha$ emission and with effective temperature. Overall, our results point to multiple origins of young star variability, including circumstellar obscuration events, hot spots on the star and/or disk, accretion bursts, and rapid structural changes in the inner disk.

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

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.

We examine the effects of pre-processing across the Coma Supercluster, including 3505 galaxies over 500 sq deg, by quantifying the degree to which star-forming (SF) activity is quenched as a function of environment. We characterise environment using the complementary techniques of Voronoi Tessellation, to measure the density field, and the Minimal Spanning Tree, to define continuous structures, and so we measure SF activity as a function of local density and the type of environment (cluster, group, filament, and void), and quantify the degree to which environment contributes to quenching of SF activity. Our sample covers over two orders of magnitude in stellar mass (10^8.5 to 10^11 Msun), and consequently we trace the effects of environment on SF activity for dwarf and massive galaxies, distinguishing so-called "mass quenching" from "environment quenching". Environmentally-driven quenching of SF activity, measured relative to the void galaxies, occurs to progressively greater degrees in filaments, groups, and clusters, and this trend holds for dwarf and massive galaxies alike. A similar trend is found using g-r colours, but with a more significant disparity between galaxy mass bins driven by increased internal dust extinction in massive galaxies. The SFR distributions of massive SF galaxies have no significant environmental dependence, but the distributions for dwarf SF galaxies are found to be statistically distinct in most environments. Pre-processing plays a significant role at low redshift, as environmentally-driven galaxy evolution affects nearly half of the galaxies in the group environment, and a significant fraction of the galaxies in the more diffuse filaments. Our study underscores the need for sensitivity to dwarf galaxies to separate mass-driven from environmentally-driven effects, and the use of unbiased tracers of SF activity.

Based on more than four weeks of continuous high cadence photometric monitoring of several hundred members of the young cluster NGC 2264 with two space telescopes, NASA's Spitzer and the CNES CoRoT (Convection, Rotation, and planetary Transits), we provide high quality, multi-wavelength light curves for young stellar objects (YSOs) whose optical variability is dominated by short duration flux bursts, which we infer are due to enhanced mass accretion rates. These light curves show many brief -- several hour to one day -- brightenings at optical and near-infrared (IR) wavelengths with amplitudes generally in the range 5-50% of the quiescent value. Typically, a dozen or more of these bursts occur in a thirty day period. We demonstrate that stars exhibiting this type of variability have large ultraviolet (UV) excesses and dominate the portion of the u-g vs. g-r color-color diagram with the largest UV excesses. These stars also have large Halpha equivalent widths, and either centrally peaked, lumpy Halpha emission profiles or profiles with blue-shifted absorption dips associated with disk or stellar winds. Light curves of this type have been predicted for stars whose accretion is dominated by Rayleigh-Taylor instabilities at the boundary between their magnetosphere and inner circumstellar disk, or where magneto-rotational instabilities modulate the accretion rate from the inner disk. Amongst the stars with the largest UV excesses or largest Halpha equivalent widths, light curves with this type of variability greatly outnumber light curves with relatively smooth sinusoidal variations associated with long-lived hot spots. We provide quantitative statistics for the average duration and strength of the accretion bursts and for the fraction of the accretion luminosity associated with these bursts.

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

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

The formation of stars in massive clusters is one of the main modes of the star formation process. However, the study of massive star forming regions is hampered by their typically large distances to the Sun. One exception to this is the massive star forming region Cygnus OB2 in the Cygnus X region, at the distance of about 1400 pc. Cygnus OB2 hosts very rich populations of massive and low-mass stars, being the best target in our Galaxy to study the formation of stars, circumstellar disks, and planets in presence of massive stars. In this paper we combine a wide and deep set of photometric data, from the r band to 24 micron, in order to select the disk bearing population of stars in Cygnus OB2 and identify the class I, class II, and stars with transition and pre-transition disks. We selected 1843 sources with infrared excesses in an area of 1 degree x 1 degree centered on Cyg OB2 in several evolutionary stages: 8.4% class I, 13.1% flat-spectrum sources, 72.9% class II, 2.3% pre-transition disks, and 3.3% transition disks. The spatial distribution of these sources shows a central cluster surrounded by a annular overdensity and some clumps of recent star formation in the outer region. Several candidate subclusters are identified, both along the overdensity and in the rest of the association.

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.

Infrared stellar photometry from 2MASS and spectral line imaging observations of 12CO and 13CO J = 1-0 line emission from the FCRAO 14m telescope are analysed to assess the variation of the CO abundance with physical conditions throughout the Orion A and Orion B molecular clouds. Three distinct Av regimes are identified in which the ratio between the 13CO column density and visual extinction changes corresponding to the photon dominated envelope, the strongly self-shielded interior, and the cold, dense volumes of the clouds. Within the strongly self-shielded interior of the Orion A cloud, the 13CO abundance varies by 100% with a peak value located near regions of enhanced star formation activity. The effect of CO depletion onto the ice mantles of dust grains is limited to regions with AV > 10 mag and gas temperatures less than 20 K as predicted by chemical models that consider thermal-evaporation to desorb molecules from grain surfaces. Values of the molecular mass of each cloud are independently derived from the distributions of Av and 13CO column densities with a constant 13CO-to-H2 abundance over various extinction ranges. Within the strongly self-shielded interior of the cloud (Av > 3 mag), 13CO provides a reliable tracer of H2 mass with the exception of the cold, dense volumes where depletion is important. However, owing to its reduced abundance, 13CO does not trace the H2 mass that resides in the extended cloud envelope, which comprises 40-50% of the molecular mass of each cloud. The implied CO luminosity to mass ratios, M/L_{CO}, are 3.2 and 2.9 for Orion A and Orion B respectively, which are comparable to the value (2.9), derived from gamma-ray observations of the Orion region. Our results emphasize the need to consider local conditions when applying CO observations to derive H2 column densities.

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

Periodic increases in luminosity arising from variable accretion rates have been predicted for some close pre-main sequence binary stars as they grow from circumbinary disks. The phenomenon is known as "pulsed accretion" and can affect the orbital evolution and mass distribution of young binaries, as well as the potential for planet formation in the circumbinary environment. Accretion variability is a common feature of young stars, with a large range of amplitudes and timescales as measured from multi-epoch observations at optical and infrared wavelengths. Periodic variations consistent with pulsed accretion have been seen in only a few young binaries via optical accretion tracers, albeit intermittently with accretion luminosity variations ranging from 0-50 percent from orbit to orbit. Here we report on a young protostar (age ~10^5 yr) that exhibits periodic variability in which the infrared luminosity increases by a factor of 10 in roughly one week every 25.34 days. We attribute this to pulsed accretion associated with an unseen binary companion. The strength and regularity of this accretion signal is surprising; it may be related to the very young age of the system, which is a factor of ten younger than the other pulsed accretors previously studied.

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

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

Recent studies have identified several young stellar objects (YSOs) which exhibit significant mid-infrared (mid-IR) variability. A wide range of physical mechanisms may be responsible for these variations, including changes in a YSO's accretion rate or in the extinction or emission from the inner disk. We have obtained and analyzed multi-epoch near-infrared (NIR) spectra for five actively accreting YSOs in the $\rho$ Oph star-forming region along with contemporaneous mid-IR light curves obtained as part of the YSOVAR Spitzer/IRAC survey. Four of the five YSOs exhibit mid-IR light curves with modest ($\sim 0.2$--0.4 mag) but statistically significant variations over our 40-day observation window. Measuring the strengths of prominent photospheric absorption lines and accretion sensitive \ion{H}{1} and \ion{He}{1} lines in each NIR spectrum, we derive estimates of each YSO's spectral type, effective temperature ({\Teff}), and $H$ band extinction ($A_H$), and analyze the time evolution of their NIR veiling ($r_H$ and $r_K$) and mass accretion rates ({\Macc}). Defining a YSO's evolutionary stage such that heavily veiled, high accretion rate objects are less evolved than those with lower levels of veiling and ongoing accretion, we infer that GY 314 is the most evolved YSO in our sample, with GY 308 and GY 292 at progressively earlier evolutionary stages. Leveraging our multi-epoch, multi-wavelength dataset, we detect significant variations in mass accretion rates over timescales of days to weeks, but find that extinction levels in these YSOs remain relatively constant. We find no correlation between these YSO mid-IR light curves and time-resolved veiling or mass accretion rates, such that we are unable to link their mid-IR variability with physical processes localized near the inner edge of the circumstellar disk or within regions which are directly responsive to mass accretion...

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.

Eclipsing binaries (EBs) provide critical laboratories for empirically testing predictions of theoretical models of stellar structure and evolution. Pre-main-sequence (PMS) EBs are particularly valuable, both due to their rarity and the highly dynamic nature of PMS evolution, such that a dense grid of PMS EBs is required to properly calibrate theoretical PMS models. Analyzing multi-epoch, multi-color light curves for 2400 candidateOrion Nebula Cluster (ONC) members from our Warm Spitzer Exploration Science Program YSOVAR, we have identified 12 stars whose light curves show eclipse features. Four of these 12 EBs are previously known. Supplementing our light curves with follow-up optical and near-infrared spectroscopy, we establish two of the candidates as likely field EBs lying behind the ONC. We confirm the remaining six candidate systems, however, as newly identified ONC PMS EBs. These systems increase the number of known PMS EBs by over 50% and include the highest mass (Theta1 Ori E, for which we provide a complete set of well-determined parameters including component masses of 2.807 and 2.797 solar masses) and longest period (ISOY J053505.71-052354.1, P \sim 20 days) PMS EBs currently known. In two cases (Theta1 Ori E and ISOY J053526.88-044730.7), enough photometric and spectroscopic data exist to attempt an orbit solution and derive the system parameters. For the remaining systems, we combine our data with literature information to provide a preliminary characterization sufficient to guide follow-up investigations of these rare, benchmark systems.

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

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

Molecular globules and pillars are spectacular features, found only in the interface region between a molecular cloud and an HII-region. Impacting Far-ultraviolet (FUV) radiation creates photon dominated regions (PDRs) on their surfaces that can be traced by typical cooling lines. With the GREAT receiver onboard SOFIA we mapped and spectrally resolved the [CII] 158 micron atomic fine-structure line and the highly excited 12CO J=11-10 molecular line from three objects in Cygnus X (a pillar, a globule, and a strong IRAS source). We focus here on the globule and compare our data with existing Spitzer data and recent Herschel Open-Time PACS data. Extended [CII] emission and more compact CO-emission was found in the globule. We ascribe this emission mainly to an internal PDR, created by a possibly embedded star-cluster with at least one early B-star. However, external PDR emission caused by the excitation by the Cyg OB2 association cannot be fully excluded. The velocity-resolved [CII] emission traces the emission of PDR surfaces, possible rotation of the globule, and high-velocity outflowing gas. The globule shows a velocity shift of ~2 km/s with respect to the expanding HII-region, which can be understood as the residual turbulence of the molecular cloud from which the globule arose. This scenario is compatible with recent numerical simulations that emphazise the effect of turbulence. It is remarkable that an isolated globule shows these strong dynamical features traced by the [CII]-line, but it demands more observational studies to verify if there is indeed an embedded cluster of B-stars.

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 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 the discovery of ten proplyd-like objects in the vicinity of the massive OB association Cygnus OB2. They were discovered in IPHAS H-Alpha images and are clearly resolved in broad-band HST/ACS, near-IR and Spitzer mid-IR images. All exhibit the familiar tadpole shape seen in photoevaporating objects such as the Orion proplyds, with a bright ionization front at the head facing the central cluster of massive stars, and a tail stretching in the opposite direction. Many also show secondary ionization fronts, complex tail morphologies or multiple heads. We consider the evidence that these are either proplyds or `evaporating gaseous globules' (EGGs) left over from a fragmenting molecular cloud, but find that neither scenario fully explains the observations. Typical sizes are 50,000--100,000 AU, larger than the Orion proplyds, but in agreement with the theoretical scaling of proplyd size with distance from the ionizing source. These objects are located at projected separations of 6-14pc from the OB association, compared to 0.1pc for the Orion proplyds, but are clearly being photoionized by the 65 O-type stars in Cyg OB2. Central star candidates are identified in near- and mid-IR images, supporting the proplyd scenario, though their large sizes and notable asymmetries is more consistent with the EGG scenario. A third possibility is therefore considered, that these are a unique class of photoevaporating partially-embedded young stellar objects that have survived the destruction of their natal molecular cloud. This has implications for the properties of stars that form in the vicinity of massive stars.

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 present 1.1 mm observations for a sample of 16 powerful radio galaxies at 0.5<z<5.2 and a radio quiet quasar at z=6.3, obtained using the AzTEC bolometer array mounted on the ASTE or the JCMT. This paper more than doubles the number of high-z radio galaxies imaged at millimetre/sub-millimetre wavelengths. We detect probable millimetre-wave counterparts for 11 of the active galaxies. The 6 active galaxies which do not have a probable millimetre counterpart in our images nevertheless have one or more likely associated millimetric source. Thus, we conclude that powerful (radio-loud) active galaxies at high-z are beacons for finding luminous millimetre/sub-millimetre galaxies at high-z. The flux densities of our AzTEC counterparts imply star formation rates ranging from <200 to ~1300 M./yr. In addition, we find that for the radio galaxoes the 1.1 mm flux density is anticorrelated with the largest angular size of the radio source. We also present new Spitzer imaging observations of several active galaxies in our sample. Combining these with archival data, we examine the mid-infrared colours of our sample. We find that radio galaxies for which we have detected a probable 1.1 mm counterpart have mid-infrared colours consistent with dusty starbursts, and are usually bluer than high-z Spitzer-selected active galaxies. In addition, we find arcs of 24 micron sources extending across ~200-500 kpc, apparently associated with three of the radio galaxies.

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 results of CO ($J=3-2$) and HCO$^+$ ($J=4-3$) mapping observations toward a nearby embedded cluster, Serpens South, using the ASTE 10 m telescope. Our CO ($J=3-2$) map reveals that many outflows are crowded in the dense cluster-forming clump that can be recognized as a HCO$^+$ clump with a size of $\sim$ 0.2 pc and mass of $\sim$ 80 M$_\odot$. The clump contains several subfragments with sizes of $\sim$ 0.05 pc. By comparing the CO ($J=3-2$) map with the 1.1 mm dust continuum image taken by AzTEC on ASTE, we find that the spatial extents of the outflow lobes are sometimes anti-correlated with the distribution of the dense gas and some of the outflow lobes apparently collide with the dense gas. The total outflow mass, momentum, and energy are estimated at 0.6 $M_\odot$, 8 $M_\odot$ km s$^{-1}$, and 64 $M_\odot$ km$^2$ s$^{-2}$, respectively. The energy injection rate due to the outflows is comparable to the turbulence dissipation rate in the clump, implying that the protostellar outflows can maintain the supersonic turbulence in this region. The total outflow energy seems only about 10 percent the clump gravitational energy. We conclude that the current outflow activity is not enough to destroy the whole cluster-forming clump, and therefore star formation is likely to continue for several or many local dynamical times.

The Serpens South embedded cluster, which is located at the constricted part in a long filamentary infrared dark cloud, is believed to be in very early stage of cluster formation. We present results of near-infrared (JHKs) polarization observations toward the filamentary cloud. Our polarization measurements of near-infrared point sources indicate a well-ordered global magnetic field that is perpendicular to the main filament, implying that the magnetic field is likely to have controlled the formation of the main filament. On the other hand, the sub-filaments, which converge on the central part of the cluster, tend to run along the magnetic field. The global magnetic field appears to be curved in the southern part of the main filament. Such morphology is consistent with the idea that the global magnetic field is distorted by gravitational contraction along the main filament toward the northern part that contains larger mass. Applying the Chandrasekhar-Fermi method, the magnetic field strength is roughly estimated to be a few x 100 microgauss, suggesting that the filamentary cloud is close to magnetically critical as a whole.

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

We present initial results from time series imaging at infrared wavelengths of 0.9 sq. degrees in the Orion Nebula Cluster (ONC). During Fall 2009 we obtained 81 epochs of Spitzer 3.6 and 4.5 micron data over 40 consecutive days. We extracted light curves with ~3% photometric accuracy for ~2000 ONC members ranging from several solar masses down to well below the hydrogen burning mass limit. For many of the stars, we also have time-series photometry obtained at optical (Ic) and/or near-infrared (JKs) wavelengths. Our data set can be mined to determine stellar rotation periods, identify new pre-main-sequence (PMS) eclipsing binaries, search for new substellar Orion members, and help better determine the frequency of circumstellar disks as a function of stellar mass in the ONC. Our primary focus is the unique ability of 3.6 & 4.5 micron variability information to improve our understanding of inner disk processes and structure in the Class I and II young stellar objects (YSOs). In this paper, we provide a brief overview of the YSOVAR Orion data obtained in Fall 2009, and we highlight our light curves for AA-Tau analogs - YSOs with narrow dips in flux, most probably due to disk density structures passing through our line of sight. Detailed follow-up observations are needed in order to better quantify the nature of the obscuring bodies and what this implies for the structure of the inner disks of YSOs.

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 (Sigma) 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 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. The observed lognormal Sigma is consistent with predictions of hierarchically structured star-formation at scales below 10 pc, arising from the molecular cloud structures. We do not find evidence for multiple discrete modes of star-formation (e.g. clustered and distributed). Comparing the observed Sigma distribution to previous Sigma threshold definitions of clusters show that they are arbitrary. 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 early results of our multiwavelength study of the star forming region IRAS 20050+2720. While we use X-rays and IR to classify young stellar objects, the optical data can be used to exclude foreground objects. The dataset set contains 57 class I sources, 183 class II sources and 183 X-ray sources, of which 140 are class III candidates. Within IRAS 20050+2720 four subclusters are found. Subcluster A and B are the densest regions, which contain stars of all evolutionary stages. Subcluster C is much younger than the other two. It has not formed any class III objects yet. We newly identify a fourth subcluster, which consists mostly of class II objects and is located about 10 arcmin from the center of the cloud.

We present a study of the variation of spatial structure of stellar populations within dwarf galaxies as a function of the population age. We use deep Hubble Space Telescope/Advanced Camera for Surveys imaging of nearby dwarf galaxies in order to resolve individual stars and create composite colour-magnitude diagrams (CMDs) for each galaxy. Using the obtained CMDs, we select Blue Helium Burning stars (BHeBs), which can be unambiguously age-dated by comparing the absolute magnitude of individual stars with stellar isochrones. Additionally, we select a very young (<10 Myr) population of OB stars for a subset of the galaxies based on the tip of the young main-sequence. By selecting stars in different age ranges we can then study how the spatial distribution of these stars evolves with time. We find, in agreement with previous studies, that stars are born within galaxies with a high degree of substructure which is made up of a continuous distribution of clusters, groups and associations from parsec to hundreds of parsec scales. These structures disperse on timescales of tens to hundreds of Myr, which we quantify using the two-point correlation function and the Q-parameter developed by Cartwright & Whitworth (2004). On galactic scales, we can place lower limits on the time it takes to remove the original structure (i.e., structure survives for at least this long), tevo, which varies between ~100~Myr (NGC~2366) and ~350 Myr (DDO~165). This is similar to what we have found previously for the SMC (~80~Myr) and the LMC (~175 Myr). We do not find any strong correlations between tevo and the luminosity of the host galaxy.

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 IRAC/MIPS Spitzer Space Telescope observations of the solar type and the low mass stellar population of the young (~5 Myr) Lambda Orionis cluster. Combining optical and 2MASS photometry, we identify 436 stars as probable members of the cluster. Given the distance (450 pc) and the age of the cluster, our sample ranges in mass from 2 solar mass to objects below the substellar limit. With the addition of the Spitzer mid-infrared data, we have identified 49 stars bearing disks in the stellar cluster. Using spectral energy distribution (SED) slopes, we place objects in several classes: non-excess stars (diskless), stars with optically thick disks, stars with ``evolved disks''( with smaller excesses than optically thick disk systems), and ``transitional disks'' candidates (in which the inner disk is partially or fully cleared). The disk fraction depends on the stellar mass, ranging from ~6% for K type stars (Rc-J<2) to ~27% for stars with spectral type M5 or later (Rc-J>4). We confirm the dependence of disk fraction on stellar mass in this age range found in other studies. Regarding clustering levels, the overall fraction of disks in the Lambda Orionis cluster is similar to those reported in other stellar groups with ages normally quoted as ~5 Myr.

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.

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 IRAC/MIPS Spitzer observations of intermediate-mass stars in the 5 Myr old Lambda Orionis cluster. In a representative sample of stars earlier than F5 (29 stars), we find a population of 9 stars with a varying degree of moderate 24um excess comparable to those produced by debris disks in older stellar groups. As expected in debris disks systems, those stars do not exhibit emission lines in their optical spectra. We also include in our study the star HD 245185, a known Herbig Ae object which displays excesses in all Spitzer bands and shows emission lines in its spectrum. We compare the disk population in the Lambda Orionis cluster with the disk census in other stellar groups studied using similar methods to detect and characterize their disks and spanning a range of ages from 3 Myr to 10 Myr. We find that for stellar groups of 5 Myr or older the observed disk frequency in intermediate mass stars (with spectral types from late B to early F) is higher than in low mass stars (with spectral types K and M). This is in contradiction with the observed trend for primordial disks evolution, in which stars with higher stellar masses dissipate their primordial disks faster. At 3 Myr the observed disk frequency in intermediate mass stars is still lower than for low mass stars indicating that second generation dusty disks start to dominate the disk population at 5 Myr for intermediate mass stars. This result agrees with recent models of evolution of solids in the region of the disk where icy objects form (>30 AU), which suggest that at 5-10 Myr collisions start to produce large amount of dust during the transition from runaway to oligarchic growth (reaching sizes of ~500 km) and then dust production peaks at 10-30 Myr, when objects reach their maximum sizes (>1000 km)

We present results of optical spectroscopic and BVR_CI_C photometric observations of 77 pre-main sequence (PMS) stars in the Cepheus flare region. A total of 64 of these are newly confirmed PMS stars, originally selected from various published candidate lists. We estimate effective temperatures and luminosities for the PMS stars, and comparing the results with pre-main sequence evolutionary models we estimate stellar masses of 0.2-2.4M_sun and stellar ages of 0.1-15 Myr. Among the PMS stars, we identify 15 visual binaries with separations of 2-10 arcsec. From archival IRAS, 2MASS, and Spitzer data, we construct their spectral energy distributions and classify 5% of the stars as Class I, 10% as Flat SED, 60% as Class II, and 3% as Class III young stellar objects (YSOs). We identify 12 CTTS and 2 WTTS as members of NGC 7023, with mean age of 1.6 Myr. The 13 PMS stars associated with L1228 belong to three small aggregates: RNO 129, L1228A, and L1228S. The age distribution of the 17 PMS stars associated with L1251 suggests that star formation has propagated with the expansion of the Cepheus flare shell. We detect sparse aggregates of 6-7 Myr old PMS stars around the dark clouds L1177 and L1219, at a distance of 400 pc. Three T Tauri stars appear to be associated with the Herbig Ae star SV Cep at a distance of 600 pc. Our results confirm that the molecular complex in the Cepheus flare region contains clouds of various distances and star forming histories.

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

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

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

We present 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 have conducted a survey of a sample of infrared-dark clouds (IRDCs) with the Spitzer Space Telescope in order to explore their mass distribution. We present a method for tracing mass using dust absorption against the bright Galactic background at 8 microns. The IRDCs in this sample are comprised of tens of clumps, ranging in sizes from 0.02 to 0.3 pc in diameter and masses from 0.5 to a few 10 Msun, the broadest dynamic range in any clump mass spectrum study to date. Structure with this range in scales confirms that IRDCs are the the precursors to stellar clusters in an early phase of fragmentation. Young stars are distributed in the vicinity of the IRDCs, but the clumps are typically not associated with stars and appear pre-stellar in nature. We find an IRDC clump mass spectrum with a slope of 1.76 +/- 0.05 for masses from 30 to 3000 Msun. This slope is consistent with numerous studies, culled from a variety of observational techniques, of massive star formation regions and is close to the mass function of Galactic stellar clusters and star clusters in other galaxies. We assert that the shape of the mass function is an intrinsic and universal feature of massive star formation regions, that are the birth sites of stellar clusters. As these clouds evolve and their constituent clumps fragment, the mass spectrum will steepen and eventually assume the form of the core mass function that is observed locally.

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

We present an analysis of the spatial distribution of various stellar populations within the Large Magellanic Cloud. We combine mid-infrared selected young stellar objects, optically selected samples with mean ages between ~9 and ~1000 Myr, and existing stellar cluster catalogues to investigate how stellar structures form and evolve within the LMC. For the analysis we use Fractured Minimum Spanning Trees, the statistical Q parameter, and the two-point correlation function. Restricting our analysis to young massive (OB) stars we confirm our results obtained for M33, namely that the luminosity function of the groups is well described by a power-law with index -2, and that there is no characteristic length-scale of star-forming regions. We find that stars in the LMC are born with a large amount of substructure, consistent with a 2D fractal distribution with dimension ~1.8 and evolve towards a uniform distribution on a timescale of ~175 Myr. This is comparable to the crossing time of the galaxy and we suggest that stellar structure, regardless of spatial scale, will be eliminated in a crossing time. This may explain the smooth distribution of stars in massive/dense young clusters in the Galaxy, while other, less massive, clusters still display large amounts of structure at similar ages. By comparing the stellar and star cluster distributions and evolving timescales, we show that infant mortality of clusters (or 'popping clusters') have a negligible influence on galactic structure. Finally, we quantify the influence of the elongation, differential extinction, and contamination of a population on the measured Q value.

We present an analysis of the spatial distribution of various stellar populations within the Large and Small Magellanic Clouds. We use optically selected stellar samples with mean ages between ~9 and ~1000 Myr, and existing stellar cluster catalogues to investigate how stellar structures form and evolve within the LMC/SMC. We use two statistical techniques to study the evolution of structure within these galaxies, the $Q$-parameter and the two-point correlation function (TPCF). In both galaxies we find the stars are born with a high degree of substructure (i.e. are highly fractal) and that the stellar distribution approaches that of the 'background' population on timescales similar to the crossing times of the galaxy (~80/150 Myr for the SMC/LMC respectively). By comparing our observations to simple models of structural evolution we find that 'popping star clusters' do not significantly influence structural evolution in these galaxies. Instead we argue that general galactic dynamics are the main drivers, and that substructure will be erased in approximately the crossing time, regardless of spatial scale, from small clusters to whole galaxies. This can explain why many young Galactic clusters have high degrees of substructure, while others are smooth and centrally concentrated. We conclude with a general discussion on cluster 'infant mortality', in an attempt to clarify the time/spatial scales involved.

We present images and initial results from our extensive Spitzer Space Telescope imaging survey of the W5 H II region with the Infrared Array Camera (IRAC) and Multiband Imaging Photometer for Spitzer (MIPS). We detect dense clusters of stars, centered on the O stars: HD 18326, BD +60 586, HD 17505 and HD 17520. At 24 microns substantial extended emission is visible, presumably from heated dust grains that survive in the strongly ionizing environment of the H {\sc ii} region. With photometry of more than 18000 point sources, we analyze the clustering properties of objects classified as young stars by their IR spectral energy distributions (a total of 2064 sources) across the region using a minimal-spanning-tree algorithm. We find ~40--70% of infrared excess sources belong to clusters with >10 members. We find that within the evacuated cavities of the H II regions that make up W5, the ratio of Class II to Class I sources is ~7 times higher than for objects coincident with molecular gas as traced by 12CO emission and near-IR extinction maps. We attribute this contrast to an age difference between the two locations, and postulate that at least two distinct generations of star formation are visible across W5. Our preliminary analysis shows that triggering is a plausible mechanism to explain the multiple generations of star formation in W5, and merits further investigation.

We present new Spitzer Space Telescope observations of stars in the young ~5 Myr gamma Velorum stellar cluster. Combining optical and 2MASS photometry, we have selected 579 stars as candidate members of the cluster. With the addition of the Spitzer mid-infrared data, we have identified 5 debris disks around A-type stars, and 5-6 debris disks around solar-type stars, indicating that the strong radiation field in the cluster does not completely suppress the production of planetesimals in the disks of cluster members. However, we find some evidence that the frequency of circumstellar primordial disks is lower, and the IR flux excesses are smaller than for disks around stellar populations with similar ages. This could be evidence for a relatively fast dissipation of circumstellar dust by the strong radiation field from the highest mass star(s) in the cluster. Another possibility is that gamma Velorum stellar cluster is slightly older than reported ages and the the low frequency of primordial disks reflects the fast disk dissipation observed at ~5 Myr.

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.

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 report on the results of an approximately 90 ks Chandra observation of a complex region that hosts multiple sites of recent and active star formation in ARA OB1a. The field is centered on the embedded cluster RCW 108-IR and includes and a large portion of the open cluster NGC 6193. We detect over 420 X-ray sources in the field and combined these data with deep near-IR, Spitzer/IRAC and MSX mid-IR data. We find about 360 of the X-ray sources have near--IR counterparts. We divide the region into 5 parts based on the X-ray point source characteristics and extended 8 micron emission. The most clearly defined regions are the central region - identified by embedded sources with high luminosities in the both the near-IR and X-ray as well as high X-ray temperatures (about 3 keV) and the eastern region - identified by low extinction and 1 keV X-ray temperatures. Other regions, identified by their directional relationship to RCW 108-IR are less uniform - representing combinations of the first two regions, independent star formation epochs, or both. Over 18% percent of the cluster members with over 100 counts exhibit flares. Overall about 50% of the stars appear to have optically thick disks when IRAC data are employed. The largest fraction of X-ray sources are best described as possessing some disk material via a more detailed extinction fitting. We estimate that the total number of pre--main sequence stars in the field is about 1600. Approximately 800 are confined to (1.1 pc) central region.

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 new Spitzer Space Telescope observations of two fields in the Orion OB1 association. We report here IRAC/MIPS observations for 115 confirmed members and 41 photometric candidates of the ~10 Myr 25 Orionis aggregate in the OB1a subassociation, and 106 confirmed members and 65 photometric candidates of the 5 Myr region located in the OB1b subassociation. The 25 Orionis aggregate shows a disk frequency of 6% while the field in the OB1b subassociation shows a disk frequency of 13%. Combining IRAC, MIPS and 2MASS photometry we place stars bearing disks in several classes: stars with optically thick disks (class II systems), stars with an inner transitional disks (transitional disk candidates) and stars with "evolved disks"; the last exhibit smaller IRAC/MIPS excesses than class II systems. In all, we identify 1 transitional disk candidate in the 25 Orionis aggregate and 3 in the OB1b field; this represents ~10% of the disk bearing stars, indicating that the transitional disk phase can be relatively fast. We find that the frequency of disks is a function of the stellar mass, suggesting a maximum around stars with spectral type M0. Comparing the infrared excess in the IRAC bands among several stellar groups we find that inner disk emission decays with stellar age, showing a correlation with the respective disk frequencies. The disk emission at the IRAC and MIPS bands in several stellar groups indicates that disk dissipation takes place faster in the inner region of the disks. Comparison with models of irradiated accretion disks, computed with several degrees of settling, suggests that the decrease in the overall accretion rate observed in young stellar groups is not sufficient to explain the weak disk emission observed in the IRAC bands for disk bearing stars with ages 5 Myr or older.

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.

Star-formation within galaxies appears on multiple scales, from spiral structure, to OB associations, to individual star clusters, and often sub-structure within these clusters. This multitude of scales calls for objective methods to find and classify star-forming regions, regardless of spatial size. To this end, we present an analysis of star-forming groups in the local group spiral galaxy M33, based on a new implementation of the Minimum Spanning Tree (MST) method. Unlike previous studies which limited themselves to a single spatial scale, we study star-forming structures from the effective resolution limit (~20pc) to kpc scales. We find evidence for a continuum of star-forming group sizes, from pc to kpc scales. We do not find a characteristic scale for OB associations, unlike that found in previous studies, and we suggest that the appearance of such a scale was caused by spatial resolution and selection effects. The luminosity function of the groups is found to be well represented by a power-law with an index, -2, similar to that found for clusters and GMCs. Additionally, the groups follow a similar mass-radius relation as GMCs. The size distribution of the groups is best described by a log-normal distribution and we show that within a hierarchical distribution, if a scale is selected to find structure, the resulting size distribution will have a log-normal distribution. We find an abrupt drop of the number of groups outside a galactic radius of ~4kpc, suggesting a change in the structure of the star-forming ISM, possibly reflected in the lack of GMCs beyond this radius. (abridged)

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.

We report new Spitzer Space Telescope observations from the IRAC and MIPS instruments of the young (~ 3 Myr) sigma Orionis cluster. We identify 336 stars as members of the cluster using optical and near-infrared color magnitude diagrams. Using the spectral energy distribution (SED) slopes in the IRAC spectral range, we place objects in several classes: non-excess stars, stars with optically thick disks(like classical T Tauri stars), class I (protostellar) candidates, and stars with ``evolved disks''; the last exhibit smaller IRAC excesses than optically thick disk systems. In general, this classification agrees with the location expected in IRAC-MIPS color-color diagrams for these objects. We find that the evolved disk systems are mostly a combination of objects with optically thick but non-flared disks, suggesting grain growth and/or settling, and transition disks, systems in which the inner disk is partially or fully cleared of small dust. In all, we identify 7 transition disk candidates and 3 possible debris disk systems. As in other young stellar populations, the fraction of disks depends on the stellar mass, ranging from ~10% for stars in the Herbig Ae/Be mass range (>2 msun) to ~35% in the T Tauri mass range (1-0.1 msun). We find that the disk fraction does not decrease significantly toward the brown dwarf candidates (<0.1 msun). The IRAC infrared excesses found in stellar clusters and associations with and without central high mass stars are similar, suggesting that external photoevaporation is not very important in many clusters. Finally, we find no correlation between the X-ray luminosity and the disk infrared excess, suggesting that the X-rays are not strongly affected by disk accretion.

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

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

Outflows from young stellar objects have been identified as a possible source of turbulence in molecular clouds. To investigate the relationship between outflows, cloud dynamics and turbulence, we compare the kinematics of the molecular gas associated with NGC 1333, traced in 13CO(1-0), with the distribution of young stellar objects (YSOs) within. We find a velocity dispersion of ~ 1-1.6 km/s in 13CO that does not significantly vary across the cloud, and is uncorrelated with the number of nearby young stellar outflows identified from optical and submillimeter observations. However, from velocity channel maps we identify about 20 cavities or depressions in the 13CO intensity of scales > 0.1-0.2 pc and velocity widths 1-3 km/s. The cavities exhibit limb brightened rims in both individual velocity channel maps and position velocity diagrams, suggesting that they are slowly expanding. We interpret these cavities to be remnants of past YSO outflow activity: If these cavities are presently empty, they would fill in on time scales of a million years. This can exceed the lifetime of a YSO outflow phase, or the transit time of the central star through the cavity, explaining the the absence of any clear correlation between the cavities and YSO outflows. We find that the momentum and energy deposition associated with the expansion of the cavities is sufficient to power the turbulence in the cloud. In this way we conclude that the cavities are an important intermediary step between the conversion of YSO outflow energy and momentum into cloud turbulent motions.

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

In the inner 3kpc of M51 we find that logarithmic spirals provide good fits to the peak intensities in molecular gas observed by BIMA in the CO (J=1-0) emission line along the spiral arms. However, we measure significant asymmetries between the location and density of the arms observed on one side of the galaxy compared to those on the opposite side. Between a radius of 1 and 2.2 kpc the gas distribution traced in CO is lopsided with densities twice as large in one arm compared to the opposite one. In the same region, the spiral arms are offset by 20 degrees from the position of the arm on the opposite side of the galaxy after a rotation of 180 degrees. We use the ratio of CO emission to that in the Pa alpha hydrogen recombination line to estimate the efficiency of star formation along each arm. Except for a narrow region at about 2kpc where star formation is surprisingly inefficient, the gas depletion time is approximately 200 million years despite large variations in the molecular gas density in the arms. We account for the deviations from two armed bi-symmetric structure with a simple model that includes an additional three-armed spiral density wave. This model accounts for the angular offset between the arm on one side compared to its opposite, the lopsided distribution in molecular gas, and inter-arm star formation to the north east of the galaxy nucleus. Since the star formation efficiency is unaffected by the variations in gas density, and the variations in gas density can be accounted for by the presence of an additional three-armed density wave, we suggest that the star formation rate is variable, and is highest where and when the maxima of the two spiral density waves coincide or constructively add.