Philip C. Myers

From Molecular Cores to Planet‐forming Disks: An SIRTF Legacy Program

Crucial steps in the formation of stars and planets can be studied only at mid‐ to far‐infrared wavelengths, where the Space Infrared Telescope (SIRTF) provides an unprecedented improvement in sensitivity. We will use all three SIRTF instruments (Infrared Array Camera [IRAC], Multiband Imaging Photometer for SIRTF [MIPS], and Infrared Spectrograph [IRS]) to observe sources that span the evolutionary sequence from molecular cores to protoplanetary disks, encompassing a wide range of cloud masses, stellar masses, and star‐forming environments. In addition to targeting about 150 known compact cores, we will survey with IRAC and MIPS (3.6–70 μm) the entire areas of five of the nearest large molecular clouds for new candidate protostars and substellar objects as faint as 0.001 solar luminosities. We will also observe with IRAC and MIPS about 190 systems likely to be in the early stages of planetary system formation (ages up to about 10 Myr), probing the evolution of the circumstellar dust, the raw material for planetary cores. Candidate planet‐forming disks as small as 0.1 lunar masses will be detectable. Spectroscopy with IRS of new objects found in the surveys and of a select group of known objects will add vital information on the changing chemical and physical conditions in the disks and envelopes. The resulting data products will include catalogs of thousands of previously unknown sources, multiwavelength maps of about 20 deg^2 of molecular clouds, photometry of about 190 known young stars, spectra of at least 170 sources, ancillary data from ground‐based telescopes, and new tools for analysis and modeling. These products will constitute the foundations for many follow‐up studies with ground‐based telescopes, as well as with SIRTF itself and other space missions such as SIM, JWST, Herschel, and TPF/Darwin.

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

SPITZER OBSERVATIONS OF IC 348: THE DISK POPULATION AT 2-3 MILLION YEARS

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 disk-bearing stars in the cluster to be 50% ± 6%. However, only 30% ± 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% ± 8% for stars earlier than K6 to 47% ± 12% for K6-M2 stars to 28% ± 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.

FILAMENTARY STRUCTURE OF STAR-FORMING COMPLEXES

The nearest young stellar groups are associated with hubs of column density exceeding 1022 cm?2, according to recent observations. These hubs radiate multiple filaments of parsec length, having lower column density and fewer stars. Systems with many filaments tend to have parallel filaments with similar spacing. Such hub-filament structure is associated with all of the nine young stellar groups within 300 pc, forming low-mass stars. Similar properties are seen in infrared dark clouds forming more massive stars. In a new model, an initial clump in a uniform medium is compressed into a self-gravitating, modulated layer. The outer layer resembles the modulated equilibrium of Schmid-Burgk with nearly parallel filaments. The filaments converge onto the compressed clump, which collapses to form stars with high efficiency. The initial medium and condensations have densities similar to those in nearby star-forming clouds and clumps. The predicted structures resemble observed hub-filament systems in their size, shape, and column density, and in the appearance of their filaments. These results suggest that HFS associated with young stellar groups may arise from compression of clumpy gas in molecular clouds.

Early evolution of stellar groups and clusters: environmental effects on forming planetary systems

This paper studies the dynamical evolution of young groups/clusters, with N = 100-1000 members, from their embedded stage out to ages of ~10 Myr. We use N-body simulations to explore how their evolution depends on the system size N and the initial conditions. Motivated by recent observations suggesting that stellar groups begin their evolution with subvirial speeds, this study compares subvirial starting states with virial starting states. Multiple realizations of equivalent cases (100 simulations per initial condition) are used to build up a robust statistical description of these systems, e.g., the probability distribution of closest approaches, the mass profiles, and the probability distribution for the radial location of cluster members. These results provide a framework from which to assess the effects of groups/clusters on the processes of star and planet formation and to study cluster evolution. The distributions of radial positions are used in conjunction with the probability distributions of the expected far-ultraviolet (FUV) luminosities (calculated here as a function of cluster size N) to determine the radiation exposure of circumstellar disks. The distributions of closest approaches are used in conjunction with scattering cross sections (calculated here as a function of stellar mass using ~105 Monte Carlo scattering experiments) to determine the probability of disruption for newly formed solar systems. We use the nearby cluster NGC 1333 as a test case in this investigation. The main conclusion of this study is that clusters in this size range have only a modest effect on forming planetary systems. The interaction rates are low, so that the typical solar system experiences a single encounter with closest approach distance b ~ 1000 AU. The radiation exposure is also low, with median FUV flux G0 ~ 900 (1.4 ergs s-1 cm-2), so that photoevaporation of circumstellar disks is only important beyond 30 AU. Given the low interaction rates and modest radiation levels, we suggest that solar system disruption is a rare event in these clusters.

Identifying the Low-Luminosity Population of Embedded Protostars in the c2d Observations of Clouds and Cores

We present the results of a search for all embedded protostars with internal luminosities ≤1.0 L☉ in the full sample of nearby, low-mass star-forming regions surveyed by the Spitzer Space Telescope Legacy Project From Molecular Cores to Planet Forming Disks (c2d). The internal luminosity of a source, Lint, is the luminosity of the central source and excludes luminosity arising from external heating. On average, the Spitzer c2d data are sensitive to embedded protostars with -->Lint ≥ 4 × 10−3(d/140 pc)2 L☉, a factor of 25 better than the sensitivity of the Infrared Astronomical Satellite (IRAS) to such objects. We present a set of selection criteria used to identify candidates from the Spitzer data and examine complementary data to decide whether each candidate is truly an embedded protostar. We find a tight correlation between the 70 μm flux and internal luminosity of a protostar, an empirical result based on both observations and detailed two-dimensional radiative transfer models of protostars. We identify 50 embedded protostars with -->Lint ≤ 1.0 L☉; 15 have -->Lint ≤ 0.1 L☉. The intrinsic distribution of source luminosities increases to lower luminosities. While we find sources down to the above sensitivity limit, indicating that the distribution may extend to luminosities lower than probed by these observations, we are able to rule out a continued rise in the distribution below -->Lint = 0.1 L☉. Between 75% and 85% of cores classified as starless prior to being observed by Spitzer remain starless to our luminosity sensitivity; the remaining 15%-25% harbor low-luminosity, embedded protostars. We compile complete spectral energy distributions for all 50 objects and calculate standard evolutionary signatures (Lbol, Tbol, and Lbol/Lsmm) and argue that these objects are inconsistent with the simplest picture of star formation, wherein mass accretes from the core onto the protostar at a constant rate.

X-ray sources in regions of star formation. III. Naked T Tauri stars associated with the Taurus-Auriga complex

Ground-based and IRAS optical and IR spectroscopic and photometric observations are reported for 90 stars in or near 59 Einstein Observatory X-ray error circles in the Tau-Aur region. The data are presented in extensive tables and sample spectra and characterized in detail, with particular attention to 28 newly discovered naked T Tau stars, which are shown to be normal stars with no significant IR or UV excess and ages of 1-40 Myr. These stars are found to outnumber normal T Tau stars by a factor of 10 in an area near the Tau-Aur dark clouds, and it is argued that their evolution toward the ZAMS is typical for low-mass stars. The implications of this finding for the time scales of circumstellar-disk dissipation and planet formation are discussed. 100 references.

THE SPITZER c2d SURVEY OF LARGE, NEARBY, INTERSTELLAR CLOUDS. IX. THE SERPENS YSO POPULATION AS OBSERVED WITH IRAC AND MIPS

We discuss the combined IRAC/MIPS c2d Spitzer Legacy observations of the Serpens star-forming region. We describe criteria for isolating bona fide YSOs from the extensive background of extragalactic objects. We then discuss the properties of the resulting high-confidence set of 235 YSOs. An additional 51 lower confidence YSOs outside this area are identified from the MIPS data and 2MASS photometry. We present color-color diagrams to compare our observed source properties with those of theoretical models for star/disk/envelope systems and our own modeling of the objects that are well represented by a stellar photosphere plus circumstellar disk. These objects exhibit a wide range of disk properties, from many with actively accreting disks to some with both passive disks and even possibly debris disks. The YSO luminosity function extends down to at least a few times 10^(-3) L_☉ or lower. The lower limit may be set more by our inability to distinguish YSOs from extragalactic sources than by the lack of YSOs at very low luminosities. We find no evidence for variability in the shorter IRAC bands between the two epochs of our data set, Δt ~ 6 hr. A spatial clustering analysis shows that the nominally less evolved YSOs are more highly clustered than the later stages. The background extragalactic population can be fitted by the same two-point correlation function as seen in other extragalactic studies. We present a table of matches between several previous infrared and X-ray studies of the Serpens YSO population and our Spitzer data set. The clusters in Serpens have a very high surface density of YSOs, primarily with SEDs suggesting extreme youth. The total number of YSOs, mostly Class II, is greater outside the clusters.

A Catalog of Optically Selected Cores

We present a new catalog of 406 dense cores optically selected by using the STScI Digitized Sky Survey (DSS). In this catalog 306 cores have neither an embedded young stellar object (EYSO) nor a pre-main-sequence (PMS) star, 94 cores have EYSOs (one core has both an EYSO and a PMS star), and six cores have PMS stars only. Our sample of dense cores in the catalog is fairly complete within a category of northern Lynds class 5 and 6 clouds and southern Hartley et al. class A clouds, providing a database useful for the systematic study of dense cores. Most of the cores listed in the catalog have diameters between 0.05 and 0.36 pc with a mean of ~0.24 pc. The sizes (~0.33 pc in the mean) of cores with EYSOs are found to be usually larger than the sizes (~0.22 pc in the mean) of starless cores. The typical mean gas density of the cores is ~7 ? 103 cm-3. Most of the cores are more likely elongated than spherical (mean aspect ratio: ~2.4). The ratio of the number of cores with EYSOs to the number of starless cores for our sample is about 0.3, suggesting that the typical lifetime of starless cores is 0.3-1.6 Myr, about 3 times longer than the duration of the class 0 and class I phases. This lifetime is shorter than expected from models of ambipolar diffusion, by factors of 2-44.