Robert Allen Gutermuth

Spitzer Observations of NGC 1333: A Study of Structure and Evolution in a Nearby Embedded Cluster

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. A total of 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 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.

Infrared Extinction toward Nearby Star-forming Regions

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 24 μm MIPS band. The source data for these measurements are observations of five nearby star-forming regions: the Orion A cloud, NGC 2068/2071, NGC 2024/2023, Serpens, and Ophiuchus. Color excess ratios E/E 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 24 μm for Serpens and NGC 2068/2071 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 to 8 μm than the previously assumed standard, as noted by all of these recent studies. The extinction law at 24 μm is consistent with previous measurements and models, although there are relatively large uncertainties.

Infrared Array Camera (IRAC) Colors of Young Stellar Objects

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. 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: S140, S171, NGC 7129, and Cep C. The IRAC color-color diagram ([3.6]� [4.5] vs. [5.8]� [8.0]) can be used to help distinguish between young stars with only disk emission and protostars with circumstellar envelopes. Subject heading gs: infrared: stars — stars: formation — stars: pre–main-sequence

The spatial distribution of star formation in the solar neighbourhood: do all stars form in dense clusters?

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

The Disk Population of the Chamaeleon I Star-forming Region*

We present a census of circumstellar disks in the Chamaeleon I star-forming region. Using the Infrared Array Camera and the Multiband Imaging Photometer on board the Spitzer Space Telescope, we have obtained images of Chamaeleon I at 3.6, 4.5, 5.8, 8.0, and 24 ?m. 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 ( -->AJ ~ 5) and reside near known protostars, indicating that they may be among the youngest low-mass sources in the cluster ( -->? M 1 M? 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.

A Spitzer View of Protoplanetary Disks in the γ Velorum Cluster

We present new Spitzer Space Telescope observations of stars in the young (~5 Myr)γ 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 five debris disks around A-type stars and five to six 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 infrared 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 γ Velorum stellar cluster is slightly older than reported ages and the low frequency of primordial disks reflects the fast disk dissipation observed at ~5 Myr.

THE SPITZER SPACE TELESCOPE SURVEY OF THE ORION A AND B MOLECULAR CLOUDS. I. A CENSUS OF DUSTY YOUNG STELLAR OBJECTS AND A STUDY OF THEIR MID-INFRARED VARIABILITY

We present a survey of the Orion A and B molecular clouds undertaken with the IRAC and MIPS instruments on board Spitzer. In total, five distinct fields were mapped, covering 9 deg^2 in five mid-IR bands spanning 3-24 μm. 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 Two Micron All Sky Survey 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 two to three epochs over six months; we search for variability between the epochs by looking for correlated variability in the 3.6 and 4.5 μm 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 μm 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.

The Initial Configuration of Young Stellar Clusters: A K-Band Number Counts Analysis of the Surface Density of Stars

We present an analysis of stellar distributions for the young stellar clusters GGD 12-15, IRAS 20050+2720, and NGC 7129, which range in far-IR luminosity from 227 to 5:68 ; 10 3 Land are all still associated with their natal molecular clouds. The data used for this analysis include near-IR data obtained with FLAMINGOS on the MMTand newlyobtainedwide-field 850 � memissionmaps from SCUBA ontheJCMT.Clustersizeandazimuthal asymmetry are measured via azimuthal and radial averaging methods, respectively. To quantify the deviation of the distribution of stars from circular symmetry, we define an azimuthal asymmetry parameter, and we investigate the statistical properties of this parameter through Monte Carlo simulations. The distribution of young stars is compared to the morphology of the molecular gas using stellar surface density maps and the 850 � m maps. We find that two of the clusters are not azimuthally symmetric and show a high degree of structure. The GGD 12-15 cluster is elongated and is aligned with newly detected filamentary structure at 850 � m. IRAS 20050+2720 is composed of a chain of three subclusters, in agreement with Chen and coworkers, although our results show that two of the subclusters appear to overlap. Significant 850 � m emission is detected toward two of the subclusters but is not detected toward the central subcluster, suggesting that the dense gas may already be cleared there. In contrast to these two highly embedded subclusters, wefind an anticorrelation of the stars and dust in NGC 7129, indicating that much of the parental gas and dust has been dispersed. The NGC 7129 cluster exhibits a higher degree of azimuthal symmetry, a lower stellar sur- face density, and a larger size than the other two clusters, suggesting that the cluster may be dynamically expanding following the recent dispersal of natal molecular gas. These analyses are further evidence that embedded, forming clusters are often not spherically symmetric structures but can be elongated and clumpy and that these morphologies may reflect the initial structure of the dense molecular gas. Furthermore, this work suggests that gas expulsion by stellar feedback results in significant dynamical evolution within the first 3 Myr of cluster evolution. We estimate peak stellar volume densities and discuss the impact of these densities on the evolution of circumstellar disks and protostellar envelopes.

YSOVAR: THE FIRST SENSITIVE, WIDE-AREA, MID-INFRARED PHOTOMETRIC MONITORING OF THE ORION NEBULA CLUSTER

We present initial results from time-series imaging at infrared wavelengths of 0.9 deg^2 in the Orion Nebula Cluster (ONC). During Fall 2009 we obtained 81 epochs of Spitzer 3.6 and 4.5 μm 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 (I_c) and/or near-infrared (JK_s ) wavelengths. Our data set can be mined to determine stellar rotation periods, identify new pre-main-sequence 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 and 4.5 μm 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 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.

Spitzer Observations of the Orion OB1 Association: Disk Census in the Low-Mass Stars

We present new Spitzer observations of two fields in the Orion OB1 association. We report 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: those with optically thick disks (class II systems), with an inner transitional disks (transitional disk candidates), and with evolved disks; the last exhibit smaller IRAC/MIPS excesses than class II systems. In all, we identify one transitional disk candidate in the 25 Orionis aggregate and three 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; larger degrees of dust settling are necessary to explain these objects.