We present a grid of radiation transfer models of axisymmetric young stellar objects (YSOs), covering a wide range of stellar masses (from 0.1 to 50 M ) and evolutionary stages (from the early envelope infall stage to the late disk-only stage). The grid consists of 20,000 YSO models, with spectral energy distributions (SEDs) and polarization spectra computed at 10 viewing angles for each model, resulting in a total of 200,000 SEDs. We have made a careful assessment of the theoretical and observational constraints on the physical conditions of disks and envelopes in YSOs and have attempted to fully span the corresponding regions in parameter space. These models are publicly available on a dedicated Web server. In this paper we summarize the main features of our models, as well as the range of parameters explored. Having a large grid covering reasonable regions of parameter space allows us to shed light on many trends in near-and mid-IR observations of YSOs (such as changes in the spectral indices and colors of their SEDs), linking them with physical parameters (such as disk and infalling envelope parameters). In particular, we examine the dependence of the spectral indices of the model SEDs on envelope accretion rate and disk mass. In addition, we show variations of spectral indices with stellar temperature, disk inner radius, and disk flaring power for a subset of disk-only models. We also examine how changing the wavelength range of data used to calculate spectral indices affects their values. We show sample color-color plots of the entire grid as well as simulated clusters at various distances with typical Spitzer sensitivities. We find that young embedded sources generally occupy a large region of color-color space due to inclination and stellar temperature effects. Disk sources occupy a smaller region of color-color space but overlap substantially with the region occupied by embedded sources, especially in the near-and mid-IR. We identify regions in color-color space where our models indicate that only sources at a given evolutionary stage should lie. We find that, while near-IR (such as JHK ) and mid-IR (such as IRAC) fluxes are useful in discriminating between stars and YSOs, and are useful for identifying very young sources, the addition of longer wavelength data such as MIPS 24 m is extremely valuable for determining the evolutionary stage of YSOs.

We present a method to analyze the spectral energy distributions (SEDs) of young stellar objects (YSOs). Our approach is to fit data with precomputed two-dimensional (2D) radiation transfer models spanning a large region of parameter space. This allows us to determine not only a single set of physical parameter values but the entire range of values consistent with the multiwavelength observations of a given source. In this way we hope to avoid any overinterpretation when modeling a set of data. We have constructed spectral energy distributions from optical to submillimeter wavelengths, including new Spitzer IRAC and MIPS photometry, for 30 young and spatially resolved sources in the Taurus-Auriga star-forming region. We demonstrate fitting model SEDs to these sources and find that we correctly identify the evolutionary stage and physical parameters found from previous independent studies, such as disk mass, disk accretion rate, and stellar temperature. We also explore how fluxes at various wavelengths help to constrain physical parameters and show examples of degeneracies that can occur when fitting SEDs. A Web-based version of this tool is available to the community.

We present near-and mid-infrared photometry obtained with the Spitzer Space Telescope of $300 known members of the IC 348 cluster. We merge this photometry with existing ground-based optical and near-infrared photometry in order to construct optical-infrared spectral energy distributions (SEDs) for all the cluster members and present a complete atlas of these SEDs. We employ these observations to investigate both the frequency and nature of the circumstellar disk population in the cluster. The Spitzer observations span a wavelength range between 3.6 and 24 m, corresponding to disk radii of $0.1-5 AU from the central star. The observations are sufficiently sensitive to enable the first detailed measurement of the disk frequency for very low mass stars at the peak of the stellar initial mass function. Using measurements of infrared excess between 3.6 and 8.0 m, we find the total frequency of diskbearing stars in the cluster to be 50% AE 6%. However, only 30% AE 4% of the member stars are surrounded by optically thick, primordial disks, while the remaining disk-bearing stars are surrounded by what appear to be optically thin, anemic disks. Both these values are below previous estimates for this cluster. The disk fraction appears to be a function of spectral type and stellar mass. The fraction of stars with optically thick disks ranges from 11% AE 8% for stars earlier than K6 to 47% AE 12% for K6-M2 stars to 28% AE 5% for M2-M6 stars. The disk longevity and thus conditions for planet formation appear to be most favorable for the K6-M2 stars, which are objects of comparable mass to the Sun for the age of this cluster. The optically thick disks around later type (>M4) stars appear to be less flared than the disks around earlier type stars. This may indicate a greater degree of dust settling and a more advanced evolutionary state for the late M disk population. Finally, we find that the presence of an optically thick dust disk is correlated with gaseous accretion, as measured by the strength of H emission. A large fraction of stars classified as classical T Tauri stars possess robust, optically thick disks, and very few such stars are found to be diskless. The majority (64%) of stars classified as weak-lined T Tauri stars are found to be diskless. However, a significant fraction (12%) of these stars are found to be surrounded by thick, primordial disks. These results suggest that it is more likely for dust disks to persist in the absence of active gaseous accretion than for active accretion to persist in the absence of dusty disks.

We present 2-D radiation transfer models of Class I Protostars and show the effect of including more realistic geometries on the resulting spectral energy distributions and images. We begin with a rotationally flattened infalling envelope as our comparison model, and add a flared disk and bipolar cavity. The disk affects the spectral energy distribution most strongly at edge-on inclinations, causing a broad dip at about 10 µm (independent of the silicate feature) due to high extinction and low scattering albedo in this wavelength region. The bipolar cavities allow more direct stellar+disk radiation to emerge into polar directions, and more scattering radiation to emerge into all directions. The wavelengthintegrated flux, often interpreted as luminosity, varies with viewing angle, with pole-on viewing angles seeing 2-4 times as much flux as edge-on, depending on geometry. Thus, observational estimates of luminosity should take into account the inclination of a source. The envelopes with cavities are significantly bluer in near-IR and mid-IR color-color plots than those without cavities. Using 1-D models to interpret Class I sources with bipolar cavities would lead to an underestimate of envelope mass and an overestimate of the implied evolutionary state. We compute images at near-, mid-, and far-IR wavelengths. We find that the mid-IR colors and images are sensitive to scattering albedo, and that the flared disk shadows the midplane on large size scales at all wavelengths plotted. Finally, our models produce polarization spectra which can be used to diagnose dust properties, such as albedo variations due to grain growth. Our results of polarization across the 3.1 µm ice feature agree well with observations for ice mantles covering 5% of the radius of the grains..

We present the Coordinated Synoptic Investigation of NGC 2264, a continuous 30-day multiwavelength 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α 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. Subject headings:Electronic address: amc@ipac.caltech.edu * Based on data from the Spitzer and CoRoT missions. The CoRoT space mission was developed and is operated by the French space agency CNES, with particpiation of ESA's RSSD