Eugene Chiang

Spectral Energy Distributions of T Tauri Stars with Passive Circumstellar Disks

We derive hydrostatic, radiative equilibrium models for passive disks surrounding T Tauri stars. Each disk is encased by an optically thin layer of superheated dust grains. This layer reemits directly to space about half the stellar energy it absorbs. The other half is emitted inward and regulates the interior temperature of the disk. The heated disk flares. As a consequence, it absorbs more stellar radiation, especially at large radii, than a flat disk would. The portion of the spectral energy distribution contributed by the disk is fairly flat throughout the thermal infrared. At fixed frequency, the contribution from the surface layer exceeds that from the interior by about a factor 3 and is emitted at more than an order of magnitude greater radius. Spectral features from dust grains in the superheated layer appear in emission if the disk is viewed nearly face-on.

Optical Images of an Exosolar Planet 25 Light Years from Earth

Fomalhaut, a bright star 7.7 parsecs (25 light-years) from Earth, harbors a belt of cold dust with a structure consistent with gravitational sculpting by an orbiting planet. Here, we present optical observations of an exoplanet candidate, Fomalhaut b. Fomalhaut b lies about 119 astronomical units (AU) from the star and 18 AU of the dust belt, matching predictions of its location. Hubble Space Telescope observations separated by 1.73 years reveal counterclockwise orbital motion. Dynamical models of the interaction between the planet and the belt indicate that the planets mass is at most three times that of Jupiter; a higher mass would lead to gravitational disruption of the belt, matching predictions of its location. The flux detected at 0.8 μm is also consistent with that of a planet with mass no greater than a few times that of Jupiter. The brightness at 0.6 μm and the lack of detection at longer wavelengths suggest that the detected flux may include starlight reflected off a circumplanetary disk, with dimension comparable to the orbits of the Galilean satellites. We also observe variability of unknown origin at 0.6 μm.

Atmospheric Escape from Hot Jupiters

Photoionization heating from ultraviolet (UV) radiation incidents on the atmospheres of hot Jupiters may drive planetary mass loss. Observations of stellar Lyman-α (Lyα) absorption have suggested that the hot Jupiter HD 209458b is losing atomic hydrogen. We construct a model of escape that includes realistic heating and cooling, ionization balance, tidal gravity, and pressure confinement by the host star wind. We show that mass loss takes the form of a hydrodynamic (Parker) wind, emitted from the planets dayside during lulls in the stellar wind. When dayside winds are suppressed by the confining action of the stellar wind, nightside winds might pick up if there is sufficient horizontal transport of heat. A hot Jupiter loses mass at maximum rates of ~2 × 1012 g s–1 during its host stars pre-main-sequence phase and ~2 × 1010 g s–1 during the stars main-sequence lifetime, for total maximum losses of ~0.06% and ~0.6% of the planets mass, respectively. For UV fluxes F UV 104 erg cm–2 s–1, the mass-loss rate is approximately energy limited and scales as . For larger UV fluxes, such as those typical of T Tauri stars, radiative losses and plasma recombination force to increase more slowly as F 0.6 UV. Dayside winds are quenched during the T Tauri phase because of confinement by overwhelming stellar wind pressure. During this early stage, nightside winds can still blow if the planet resides outside the stellar Alfven radius; otherwise, even nightside winds are stifled by stellar magnetic pressure, and mass loss is restricted to polar regions. We conclude that while UV radiation can indeed drive winds from hot Jupiters, such winds cannot significantly alter planetary masses during any evolutionary stage. They can, however, produce observable signatures. Candidates for explaining why the Lyman-α photons of HD 209458 are absorbed at Doppler-shifted velocities of ±100 km s–1 include charge-exchange in the shock between the planetary and stellar winds.

The minimum-mass extrasolar nebula: in situ formation of close-in super-Earths

Close-in super-Earths, with radii R = 2-5 R_Earth and orbital periods P < 100 days, orbit more than half, and perhaps nearly all Sun-like stars in the universe. We use this omnipresent population to construct the minimum-mass extrasolar nebula (MMEN), the circumstellar disk of solar-composition solids and gas from which such planets formed, if they formed near their current locations and did not migrate. In a series of back-of-the-envelope calculations, we demonstrate how in-situ formation in the MMEN is fast, efficient, and can reproduce many of the observed properties of close-in super-Earths, including their gas-to-rock fractions. Testable predictions are discussed.

Forming Planetesimals in Solar and Extrasolar Nebulae

Planets are built from planetesimals: solids larger than a kilometer that grow by colliding in pairs. Planetesimals themselves are unlikely to form by twobody collisions alone; subkilometer objects have gravitational fields individually too weak, and electrostatic attraction is too feeble for growth beyond a few centimeters. We review the possibility that planetesimals form when self-gravity brings together vast ensembles of small particles. Even when self-gravity is weak, aerodynamic processes can accumulate solids relative to gas, paving the way for gravitational collapse. Particles pile up as they drift radially inward. Gas turbulence stirs particles but can also seed collapse by clumping them. Whereas the feedback of solids on gas triggers vertical-shear instabilities that obstruct self-gravity, this same feedback triggers streaming instabilities that strongly concentrate particles. Numerical simulations find that solids ∼10–100 cm in size gravitationally collapse in turbulent disks. We outline areas for progress, including the possibility that still smaller objects self-gravitate.

Discovery and spectroscopy of the young jovian planet 51 Eri b with the Gemini Planet Imager

An exoplanet extracted from the bright Direct imaging of Jupiter-like exoplanets around young stars provides a glimpse into how our solar system formed. The brightness of young stars requires the use of next-generation devices such as the Gemini Planet Imager (GPI). Using the GPI, Macintosh et al. discovered a Jupiter-like planet orbiting a young star, 51 Eridani (see the Perspective by Mawet). The planet, 51 Eri b, has a methane signature and is probably the smallest exoplanet that has been directly imaged. These findings open the door to understanding solar system origins and herald the dawn of a new era in next-generation planetary imaging. Science, this issue p. 64; see also p. 39 The Gemini Planet Imager detects a Jupiter-like exoplanet orbiting the young star 51 Eridani. [Also see Perspective by Mawet] Directly detecting thermal emission from young extrasolar planets allows measurement of their atmospheric compositions and luminosities, which are influenced by their formation mechanisms. Using the Gemini Planet Imager, we discovered a planet orbiting the ~20-million-year-old star 51 Eridani at a projected separation of 13 astronomical units. Near-infrared observations show a spectrum with strong methane and water-vapor absorption. Modeling of the spectra and photometry yields a luminosity (normalized by the luminosity of the Sun) of 1.6 to 4.0 × 10−6 and an effective temperature of 600 to 750 kelvin. For this age and luminosity, “hot-start” formation models indicate a mass twice that of Jupiter. This planet also has a sufficiently low luminosity to be consistent with the “cold-start” core-accretion process that may have formed Jupiter.

FOMALHAUT'S DEBRIS DISK AND PLANET: CONSTRAINING THE MASS OF FOMALHAUT B FROM DISK MORPHOLOGY

Following the optical imaging of exoplanet candidate Fomalhaut b (Fom b), we present a numerical model of how Fomalhauts debris disk is gravitationally shaped by a single interior planet. The model is simple, adaptable to other debris disks, and can be extended to accommodate multiple planets. If Fom b is the dominant perturber of the belt, then to produce the observed disk morphology it must have a mass M pl 101.5 AU, and an orbital eccentricity e pl = 0.11-0.13. These conclusions are independent of Fom bs photometry. To not disrupt the disk, a greater mass for Fom b demands a smaller orbit farther removed from the disk; thus, future astrometric measurement of Fom bs orbit, combined with our model of planet-disk interaction, can be used to determine the mass more precisely. The inner edge of the debris disk at a ≈ 133 AU lies at the periphery of Fom bs chaotic zone, and the mean disk eccentricity of e ≈ 0.11 is secularly forced by the planet, supporting predictions made prior to the discovery of Fom b. However, previous mass constraints based on disk morphology rely on several oversimplifications. We explain why our constraint is more reliable. It is based on a global model of the disk that is not restricted to the planets chaotic zone boundary. Moreover, we screen disk parent bodies for dynamical stability over the system age of ~ 100 Myr, and model them separately from their dust grain progeny; the latters orbits are strongly affected by radiation pressure and their lifetimes are limited to ~ 0.1 Myr by destructive grain-grain collisions. The single planet model predicts that planet and disk orbits be apsidally aligned. Fomalhaut bs nominal space velocity does not bear this out, but the astrometric uncertainties may be large. If the apsidal misalignment proves real, our calculated upper mass limit of 3M J still holds. If the orbits are aligned, our model predicts M pl = 0.5M J, a pl = 115 AU, and e pl = 0.12. Parent bodies are evacuated from mean-motion resonances with Fom b; these empty resonances are akin to the Kirkwood gaps opened by Jupiter. The belt contains at least 3M ⊕ of solids that are grinding down to dust, their velocity dispersions stirred so strongly by Fom b that collisions are destructive. Such a large mass in solids is consistent with Fom b having formed in situ.

The Deep Ecliptic Survey: A Search for Kuiper Belt Objects and Centaurs. II. Dynamical Classification, the Kuiper Belt Plane, and the Core Population

The Deep Ecliptic Survey (DES)?a search optimized for the discovery of Kuiper belt objects (KBOs) with the Blanco and Mayall 4 m telescopes at the Cerro Tololo Inter-American Observatory and Kitt Peak National Observatory?has covered 550 deg2 from its inception in 1998 through the end of 2003. This survey has a mean 50% sensitivity at VR magnitude 22.5. We report here the discoveries of 320 designated KBOs and Centaurs for the period 2000 March through 2003 December and describe improvements to our discovery and recovery procedures. Our data and the data products needed to reproduce our analyses in this paper are available through the NOAO survey database. Here we present a dynamical classification scheme, based on the behavior of orbital integrations over 10 Myr. The dynamical classes, in order of testing, are Resonant, Centaur, Scattered-Near, Scattered-Extended, and Classical. (These terms are capitalized when referring to our rigorous definitions.) Of the 382 total designated KBOs discovered by the DES, a subset of 196 objects have sufficiently accurate orbits for dynamical classification. Summary information is given for an additional 240 undesignated objects also discovered by the DES from its inception through the end of 2003. The number of classified DES objects (uncorrected for observational bias) are Classical, 96; Resonant, 54; Scattered-Near, 24; Scattered-Extended, 9; and Centaur, 13. We use subsets of the DES objects (which can have observational biases removed) and larger samples to perform dynamical analyses on the Kuiper belt. The first of these is a determination of the Kuiper belt plane (KBP), for which the Classical objects with inclinations less than 5? from the mean orbit pole yield a pole at R.A. = 27392 ? 062 and decl. = 6670 ? 020 (J2000), consistent with the invariable plane of the solar system. A general method for removing observational biases from the DES data set is presented and used to find a provisional magnitude distribution and the distribution of orbital inclinations relative to the KBP. A power-law model fit to the cumulative magnitude distribution of all KBOs discovered by the DES in the VR filter yields an index of 0.86 ? 0.10 (with the efficiency parameters for the DES fitted simultaneously with the population power law). With the DES sensitivity parameters fixed, we derive power-law indices of 0.74 ? 0.05, 0.52 ? 0.08, and 0.74 ? 0.15, respectively, for the Classical, Resonant, and Scattered classes. Plans for calibration of the DES detection efficiency function and DES magnitudes are discussed. The inclination distribution confirms the presence of hot and cold populations; when the geometric sin i factor is removed from the inclination distribution function, the cold population shows a concentrated core with a full width at half-maximum of approximately 46, while the hot population appears as a halo, extending beyond 30?. The inclination distribution is used to infer the KBO distribution in the sky, as a function of latitude relative to the KBP. This inferred latitude distribution is reasonably consistent with the latitude distribution derived from direct observation, but the agreement is not perfect. We find no clear boundary between the Classical and Scattered classes either in their orbital inclinations with respect to the KBP or in their power-law indices in their respective magnitude distributions. This leaves open the possibility that common processes have shaped the distribution of orbital parameters for the two classes.

On the Location of the Snow Line in a Protoplanetary Disk

In a protoplanetary disk, the inner edge of the region where the temperature falls below the condensation temperature of water is referred to as the snow line. Outside the snow line, water ice increases the surface density of solids by a factor of 4. The mass of the fastest growing planetesimal (the isolation mass) scales as the surface density to the 3/2 power. It is thought that ice-enhanced surface densities are required to make the cores of the gas giants (Jupiter and Saturn) before the disk gas dissipates. Observations of our solar systems asteroid belt suggest that the snow line occurred near 2.7 AU. In this paper we revisit the theoretical determination of the snow line. In a minimum-mass disk characterized by conventional opacities and a mass accretion rate of 10-8 M☉ yr-1, the snow line lies at 1.6-1.8 AU, just past the orbit of Mars. The minimum-mass disk, with a mass of 0.02 M☉, has a lifetime of 2 million years with the assumed accretion rate. Moving the snow line past 2.7 AU requires that we increase the disk opacity, accretion rate, and/or disk mass by factors ranging up to an order of magnitude above our assumed baseline values.

Particle Pileups and Planetesimal Formation

Solid particles in protoplanetary disks that are sufficiently supersolar in metallicity overcome turbulence generated by vertical shear to gravitationally condense into planetesimals. Supersolar metallicities result if solid particles pile up as they migrate starward as a result of aerodynamic drag. Previous analyses of aerodynamic drift rates that account for mean flow differences between gas and particles yield particle pileups. We improve on these studies not only by accounting for the collective inertia of solids relative to that of gas, but also by including the transport of angular momentum by turbulent stresses within the particle layer. These turbulent stresses are derived in a physically self-consistent manner from the structure of marginally Kelvin-Helmholtz turbulent flows. They are not calculated using the usual plate drag formulae, whose use we explain is inappropriate. Accounting for the relative inertia of solids to gas retards, but does not prevent, particle pileups and generates more spatially extended regions of metal enrichment. Turbulent transport hastens pileups. We conclude that particle pileup is a robust outcome in sufficiently passive protoplanetary disks. Connections to observations of circumstellar disks, including the Kuiper Belt, and the architectures of planetary systems are made.