Klaus M. Pontoppidan

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.

The Spitzer Ice Legacy: Ice Evolution from Cores to Protostars

Ices regulate much of the chemistry during star formation and account for up to 80% of the available oxygen and carbon. In this paper, we use the Spitzer c2d Legacy ice survey, complimented with data sets on ices in cloud cores and high-mass protostars, to determine standard ice abundances and to present a coherent picture of the evolution of ices during low- and high-mass star formation. The median ice composition H_(2)O:CO:CO_2:CH_(3)OH:NH_3:CH_4:XCN is 100:29:29:3:5:5:0.3 and 100:13:13:4:5:2:0.6 toward low- and high-mass protostars, respectively, and 100:31:38:4:-:-:- in cloud cores. In the low-mass sample, the ice abundances with respect to H_(2)O of CH_4, NH_3, and the component of CO_2 mixed with H_(2)O typically vary by <25%, indicative of co-formation with H_(2)O. In contrast, some CO and CO_2 ice components, XCN, and CH3OH vary by factors 2-10 between the lower and upper quartile. The XCN band correlates with CO, consistent with its OCN– identification. The origin(s) of the different levels of ice abundance variations are constrained by comparing ice inventories toward different types of protostars and background stars, through ice mapping, analysis of cloud-to-cloud variations, and ice (anti-)correlations. Based on the analysis, the first ice formation phase is driven by hydrogenation of atoms, which results in an H_(2)O-dominated ice. At later prestellar times, CO freezes out and variations in CO freezeout levels and the subsequent CO-based chemistry can explain most of the observed ice abundance variations. The last important ice evolution stage is thermal and UV processing around protostars, resulting in CO desorption, ice segregation, and the formation of complex organic molecules. The distribution of cometary ice abundances is consistent with the idea that most cometary ices have a protostellar origin.

c2d Spitzer IRS Spectra of Disks around T Tauri Stars. I. Silicate Emission and Grain Growth

Infrared ~5-35 μm spectra for 40 solar mass T Tauri stars and 7 intermediate-mass Herbig Ae stars with circumstellar disks were obtained using the Spitzer Space Telescope as part of the c2d IRS survey. This work complements prior spectroscopic studies of silicate infrared emission from disks, which were focused on intermediate-mass stars, with observations of solar mass stars limited primarily to the 10 μm region. The observed 10 and 20 μm silicate feature strengths/shapes are consistent with source-to-source variations in grain size. A large fraction of the features are weak and flat, consistent with micron-sized grains indicating fast grain growth (from 0.1 to 1.0 μm in radius). In addition, approximately half of the T Tauri star spectra show crystalline silicate features near 28 and 33 μm, indicating significant processing when compared to interstellar grains. A few sources show large 10-to-20 μm ratios and require even larger grains emitting at 20 μm than at 10 μm. This size difference may arise from the difference in the depth into the disk probed by the two silicate emission bands in disks where dust settling has occurred. The 10 μm feature strength versus shape trend is not correlated with age or Hα equivalent width, suggesting that some amount of turbulent mixing and regeneration of small grains is occurring. The strength versus shape trend is related to spectral type, however, with M stars showing significantly flatter 10 μm features (larger grain sizes) than A/B stars. The connection between spectral type and grain size is interpreted in terms of the variation in the silicate emission radius as a function of stellar luminosity, but could also be indicative of other spectral-type-dependent factors (e.g., X-rays, UV radiation, and stellar/disk winds).

Spectroastrometric Imaging of Molecular Gas within Protoplanetary Disk Gaps

We present velocity-resolved spectroastrometric imaging of the 4.7 µm rovibrational lines of CO gas in protoplanetary disks using the CRIRES high-resolution infrared spectrometer on the Very Large Telescope (VLT). The method as applied to three disks with known dust gaps or inner holes out to 4-45 AU (SR 21, HD 135344B, and TW Hya) achieves an unprecedented spatial resolution of 0.1-0.5 AU. While one possible gap formation mechanism is dynamical clearing by giant planets, other equally good explanations (stellar companions, grain growth, photo-evaporation) exist. One way of distinguishing between different scenarios is the presence and distribution of gas inside the dust gaps. Keplerian disk models are fit to the spectroastrometric position-velocity curves to derive geometrical parameters of the molecular gas. We determine the position angles and inclinations of the inner disks with accuracies as good as 1°-2°, as well as the radial extent of the gas emission. Molecular gas is detected well inside the dust gaps in all three disks. The gas emission extends to within a radius of 0.5 AU for HD 135344B and to 0.1 AU for TW Hya, supporting partial clearing by a < 1-10 MJup planetary body as the cause of the observed dust gaps, or removal of the dust by extensive grain coagulation and planetesimal formation. The molecular gas emission in SR 21 appears to be truncated within ~7 AU, which may be caused by complete dynamical clearing by a more massive companion. We find a smaller inclination angle of the inner disk of TW Hya than that determined for the outer disk, suggestive of a disk warp. We also detect significant azimuthal asymmetries in the SR 21 and HD 135344B inner disks.

A 3-5 mu m VLT spectroscopic survey of embedded young low mass stars I. Structure of the CO ice

Medium resolution (λ/Δλ = 5000-10000) VLT-ISAAC M-band spectra are presented of 39 young stellar objects in nearby low-mass star forming clouds showing the 4.67 μm stretching vibration mode of solid CO. By taking advantage of the unprecedentedly large sample, high S/N ratio and high spectral resolution, similarities in the ice profiles from source to source are identified. It is found that excellent fits to all the spectra can be obtained using a phenomenological decomposition of the CO stretching vibration profile at 4.67 μm into 3 components, centered on 2143.7 cm^(-1),2139.9 cm^(-1), and 2136.5 cm^(-1) with fixed widths of 3.0, 3.5 and 10.6 cm ^(-1), respectively. All observed interstellar CO profiles can thus be uniquely described by a model depending on only 3 linear fit parameters, indicating that a maximum of 3 specific molecular environments of solid CO exist under astrophysical conditions. A simple physical model of the CO ice is presented, which shows that the 2139.9 cm^(-1) component is indistinguishable from pure CO ice. It is concluded, that in the majority of the observed lines of sight, 60-90% of the CO is in a nearly pure form. In the same model the 2143.7 cm^(-1) component can possibly be explained by the longitudinal optical (LO) component of the vibrational transition in pure crystalline CO ice which appears when the background source is linearly polarised. The model therefore predicts the polarisation fraction at 4.67 μm, which can be confirmed by imaging polarimetry. The 2152 cm^(-1) feature characteristic of CO on or in an unprocessed water matrix is not detected toward any source and stringent upper limits are given. When this is taken into account, the 2136.5 cm ^(-1) component is not consistent with the available water-rich laboratory mixtures and we suggest that the carrier is not yet fully understood. A shallow absorption band centered between 2165 cm^(-1) and 2180 cn^(1) is detected towards 30 sources. For low-mass stars, this band is correlated with the CO component at 2136.5 cm^(-1), suggesting the presence of a carrier different from XCN at 2175 cm^(-1). Furthermore the absorption band from solid ^(13)CO at 2092 cm^(-1) is detected towards IRS 51 in the ρ Ophiuchi cloud complex and an isotopic ratio of ^(12)CO/^(13)CO = 68 ± 10 is derived. It is shown that all the observed solid ^(12)CO profiles, along with the solid ^(13)CO profile, are consistent with grains with an irregularly shaped CO ice mantle simulated by a Continuous Distribution of Ellipsoids (CDE), but inconsistent with the commonly used models of spherical grains in the Rayleigh limit.

An Old Disk Still Capable of Forming a Planetary System

From the masses of the planets orbiting the Sun, and the abundance of elements relative to hydrogen, it is estimated that when the Solar System formed, the circumstellar disk must have had a minimum mass of around 0.01 solar masses within about 100 astronomical units of the star. (One astronomical unit is the Earth–Sun distance.) The main constituent of the disk, gaseous molecular hydrogen, does not efficiently emit radiation from the disk mass reservoir, and so the most common measure of the disk mass is dust thermal emission and lines of gaseous carbon monoxide. Carbon monoxide emission generally indicates properties of the disk surface, and the conversion from dust emission to gas mass requires knowledge of the grain properties and the gas-to-dust mass ratio, which probably differ from their interstellar values. As a result, mass estimates vary by orders of magnitude, as exemplified by the relatively old (3–10 million years) star TW Hydrae, for which the range is 0.0005–0.06 solar masses. Here we report the detection of the fundamental rotational transition of hydrogen deuteride from the direction of TW Hydrae. Hydrogen deuteride is a good tracer of disk gas because it follows the distribution of molecular hydrogen and its emission is sensitive to the total mass. The detection of hydrogen deuteride, combined with existing observations and detailed models, implies a disk mass of more than 0.05 solar masses, which is enough to form a planetary system like our own.

H2O AND OH GAS IN THE TERRESTRIAL PLANET-FORMING ZONES OF PROTOPLANETARY DISKS

We present detections of numerous 10-20 μm H_2O emission lines from two protoplanetary disks around the T Tauri stars AS 205A and DR Tau, obtained using the InfraRed Spectrograph on the Spitzer Space Telescope. Follow-up 3-5 μm Keck NIRSPEC data confirm the presence of abundant water and spectrally resolve the lines. We also detect the P4.5 (2.934 μm) and P9.5 (3.179 μm) doublets of OH and ^(12)CO/^(13)CO v = 1 → 0 emission in both sources. Line shapes and LTE models suggest that the emission from all three molecules originates between ~0.5 and 5 AU, and so will provide a new window for understanding the chemical environment during terrestrial planet formation. LTE models also imply significant columns of H_2O and OH in the inner disk atmospheres, suggesting physical transport of volatile ices either vertically or radially, while the significant radial extent of the emission stresses the importance of a more complete understanding of nonthermal excitation processes.

A Spitzer Survey of Mid-infrared Molecular Emission from Protoplanetary Disks. I. Detection Rates

We present a Spitzer InfraRed Spectrometer search for 10-36 μm molecular emission from a large sample of protoplanetary disks, including lines from H_(2)O, OH, C_(2)H_2, HCN, and CO2. This paper describes the sample and data processing and derives the detection rate of mid-infrared molecular emission as a function of stellar mass. The sample covers a range of spectral type from early M to A, and is supplemented by archival spectra of disks around A and B stars. It is drawn from a variety of nearby star-forming regions, including Ophiuchus, Lupus, and Chamaeleon. Spectra showing strong emission lines are used to identify which lines are the best tracers of various physical and chemical conditions within the disks. In total, we identify 22 T Tauri stars with strong mid-infrared H2O emission. Integrated water line luminosities, where water vapor is detected, range from 5 × 10^(–4) to 9 × 10^(–3) L_☉ sun, likely making water the dominant line coolant of inner disk surfaces in classical T Tauri stars. None of the five transitional disks in the sample show detectable gaseous molecular emission with Spitzer upper limits at the 1% level in terms of line-to-continuum ratios (apart from H_2), but the sample is too small to conclude whether this is a general property of transitional disks. We find a strong dependence on detection rate with spectral type; no disks around our sample of 25 A and B stars were found to exhibit water emission, down to 1%-2% line-to-continuum ratios, in the mid-infrared, while more than half of disks around late-type stars (M-G) show sufficiently intense water emission to be detected by Spitzer, with a detection rate approaching 2/3 for disks around K stars. Some Herbig Ae/Be stars show tentative H_(2)O/OH emission features beyond 20 μm at the 1%-2% level, however, and one of them shows CO_2 in emission. We argue that the observed differences between T Tauri disks and Herbig Ae/Be disks are due to a difference in excitation and/or chemistry depending on spectral type and suggest that photochemistry may be playing an important role in the observable characteristics of mid-infrared molecular line emission from protoplanetary disks.

C2D Spitzer-IRS spectra of disks around T Tauri stars - IV. Crystalline silicates

Aims. Dust grains in the planet-forming regions around young stars are expected to be heavily processed due to coagulation, fragmentation, and crystallization. This paper focuses on the crystalline silicate dust grains in protoplanetary disks for a statistically significant number of TTauri stars (96). Methods. As part of the cores to disks (c2d) legacy program, we obtained more than a hundred Spitzer/IRS spectra of TTauri stars, over a spectral range of 5-35 μm where many silicate amorphous and crystalline solid-state features are present. At these wavelengths, observations probe the upper layers of accretion disks up to distances of a dozen AU from the central object. Results. More than 3/4 of our objects show at least one crystalline silicate emission feature that can be essentially attributed to Mg-rich silicates. The Fe-rich crystalline silicates are largely absent in the c2d IRS spectra. The strength and detection frequency of the crystalline features seen at λ > 20 μm correlate with each other, while they are largely uncorrelated with the observational properties of the amorphous silicate 10 μm feature. This supports the idea that the IRS spectra essentially probe two independent disk regions: a warm zone (≤1 AU) emitting at ~ 10 μm and a much colder region emitting at λ > 20 μm (≤10 AU). We identify a crystallinity paradox, as the long-wavelength (λ > 20 m) crystalline silicate features are detected 3.5 times more frequently (~55% vs. ~15%) than the crystalline features arising from much warmer disk regions (λ ~ 10 μm). This suggests that the disk has an inhomogeneous dust composition within ~10 AU. The analysis of the shape and strength of both the amorphous 10 μm feature and the crystalline feature around 23 μm provides evidence for the prevalence of μm-sized (amorphous and crystalline) grains in upper layers of disks. Conclusions. The abundant crystalline silicates found far from their presumed formation regions suggest efficient outward radial transport mechanisms in the disks around TTauri stars. The presence of μm-sized grains in disk atmospheres, despite the short timescales for settling to the midplane, suggests efficient (turbulent) vertical diffusion, probably accompanied by grain-grain fragmentation to balance the expected efficient growth. In this scenario, the depletion of submicron-sized grains in the upper layers of the disks points toward removal mechanisms such as stellar winds or radiation pressure.

Probing dust grain evolution in IM Lupi's circumstellar disc Multi-wavelength observations and modelling of the dust disc

Aims. We present a panchromatic study, involving a multiple technique approach, of the circumstellar disc surrounding the T Tauri star IM Lupi (Sz 82). Methods. We have undertaken a comprehensive observational study of IM Lupi using photometry, spectroscopy, millimetre interferometry and multi-wavelength imaging. For the first time, the disc is resolved from optical and near-infrared wavelengths in scattered light, to the millimetre regime in thermal emission. Our data-set, in conjunction with existing photometric data, provides an extensive coverage of the spectral energy distribution, including a detailed spectrum of the silicate emission bands. We have performed a simultaneous modelling of the various observations, using the radiative transfer code MCFOST, and analysed a grid of models over a large fraction of the parameter space via Bayesian inference. Results. We have constructed a model that can reproduce all of the observations of the disc. Our analysis illustrates the importance of combining a wide range of observations in order to fully constrain the disc model, with each observation providing a strong constraint only on some aspects of the disc structure and dust content. Quantitative evidence of dust evolution in the disc is obtained: grain growth up to millimetre-sized particles, vertical stratification of dust grains with micrometric grains close to the disc surface and larger grains which have settled towards the disc midplane, and possibly the formation of fluffy aggregates and/or ice mantles around grains.