arXiv:0906.3744v1 [astro-ph.EP] 19 Jun 2009 ohse,Rcetr Y14627 NY Rochester, Rochester, ioe D228 [email protected] 21218; MD timore, ihcl pcrleeg itiuin(E)clrtem- color (SED) ( distribution disks energy peratures debris spectral some cold that with suggests imaging surveys. photometric resolution band High broad from inferred currently 2007). al. et o e nw.A h urn ie nyteexo-asteroid the is only time, rare around current or belt the At common known. our are yet not belts in lower whether asteroid Bombardment however, steady-state, 2007); Heavy mass, Late al. et of (Wyatt events system period solar stochastic the massive, disks as or self-stirred such 2004) that in 5% Bromley disks formation planet & dusty remaining of (Kenyon warm, result The the massive, be com- possess may stars. are disks (Meyer sequence belts debris Belt of exo-Kuiper main Kuiper that to around our of suggesting similar in mon 2007), 95% dust temperatures that al. cold color et and for have 2000) expected disks that al. debris et these (Lagrange in- excess surveys frared far-infared and that mid- dicate been have Space-based infrared rings, thermal via dust emission. to stars sequence belt anal- believed main Kuiper around bodies, disks, detected and parent asteroid Dusty between to collisions ogous belt. quantified. by rare generated Kuiper been or be yet common a is not and architecture has planets, planetary Jovian this four Whether belt, asteroid eg R49A oahu) ohgoere a pos- may geometries Both edges inner Fomalhaut). sharp 4796A, with HR (e.g. belts narrow law (e.g. (2) power and decaying Mic) a AU and Pic, rising a a by of approximated well combination are that sur- distributions with density disks continuous face (1) architectures: thermal dust and distinct of light images scattered emission resolved Spatially ponents. 3 2 am seodlds a emr rvln than prevalent more be may dust asteroidal Warm, u oa ytmpsessfu ersra lnt,an planets, terrestrial four possesses system solar Our 1 rpittpstuigL using 2018 typeset 21, Preprint November version Draft OO 5 ot hryAeu,Tco,A 85726 AZ Tucson, Avenue, Cherry North 950 NOAO, eateto hsc n srnm,Uiest of University Astronomy, and Physics of Department pc eecp cec nttt,30 a atnD. Bal Dr., Martin San 3700 Institute, Science Telescope Space hitn .Chen H. Christine ohteosre ctee ih n hra msin ial,ca Finally, cold, the emission. if thermal the and grains, infer light recen we porous scattered the light, or observed to scattered icy the exterior in with Belt both imaged consistent Kuiper been grains, a has and scattering which to 139664, interior HD belt For asteroid an of h eta trde o cu rmadsydisk. dusty a from occur headings: not indicatin Subject clearing, does central star a central possesses the system this that suggests msinfaue ept eaieyhtinrds opnn ( F component dust components. inner body hot relatively black a temperature despite single features two emission using modeled well λ ehv rereduced have We o,H 364 n R89.W n htta h hra emissio thermal the that that find We 8799. HR and 139664, HD Boo, T ∼ ζ gr 03%o ansqec tr oss in- possess stars sequence main of 10-30% e a ensailyrsle (Moerchen resolved spatially been has Lep < OA YTMAAOSAROUND ANALOGS SYSTEM SOLAR 0 )myas oss amcom- warm possess also may K) 100 1. ∼ A 5rsle ersdssrva two reveal disks debris resolved 15 INTRODUCTION T E tl mltajv 08/22/09 v. emulateapj style X 1 tr:cruselrmte—paeaysses formation systems: planetary matter— circumstellar stars: arc Sheehan Patrick , Spitzer R pcr n enlzdteSDso he erydbi disks debris nearby three of SED’s the reanalyzed and spectra IRS rf eso oebr2,2018 21, November version Draft 2 a .Watson M. Dan , ABSTRACT β - IRAS Lue l 04.W aebe odciga conducting been have star We the of study to 2004). spectroscopic close IRS al. component et dust may (Liu observations unresolved (this interferometry an dust systems: nulling reveal warm (2) debris of and component cold study) additional in an of dust presence warm of light (1) angular presence shed may their the techniques or Additional on faint small. because too too resolved are are spatially sizes brightnesses be currently surface not their can 1998; far thus al. et 2004). (Koerner al. emission et scat- thermal in Stapelfeldt resolved and are light that rings distinct tered dust spatially narrow cold, are unre- the that of from detected components imaging have dust also Similar warm may solved dust. Fomalhaut Wahhaj warm and 2005; of 4796A al. indicative HR et 2003) (Telesco al. star et central distances the from from originating away excess infrared of nent mid- of Ground-based imaging infrared temperatures. terrestrial at dust sess fmgiuebihe hnepce rmterstellar their 24 Recent from MIPS 1993). Paresce expected & (Backman than photospheres brighter magnitude of oitdwith sociated 24 unresolved of presence H 364,adH 79wr rtietfidbased identified 60 first of were presence 8799 the HR on and 139664), (HD dust of 139664. properties al- the HD SED on around infrared constraints grains the new and place to 2006) scattered us resolved al. lows of et system synthesis (Kalas The this images of own. light our architecture of planetary that allows the 2009) with compare al. et to (Marois us 8799 HR around planets giant o,H 85(D196) R89)wt amdust discovered warm originally was with that 8799) dust using cold HR our the to 139664), of addition ( (HD in reanalysis 5825 objects and HR four Rereduction Boo, finds observations 2006). previously-published al. et (Chen notntl,temjrt ersdssdiscovered disks debris majority the Unfortunately, h ersdssaon h stars the around disks debris The 2 .MnjPuravankara Manoj P. , Spitzer DSOEE ERSDISKS DEBRIS -DISCOVERED IRAS htslcieaceino oisonto solids of accretion selective that g µ n 70 and m T h eetdrc eeto ftreorbiting three of detection direct recent The . .-40 5.5- gr λ 5 )–w ne h presence the infer we – K) 150 = eu nlsso the of analysis reful ue ikcmoetgenerates component disk outer l icvrdobtn planets. orbiting discovered tly o S ta.20)adH 139664 HD and 2006) al. et (Su Boo rH 79–wt osilicate no with – 8799 HR or 2. µ rsneo togyforward strongly of presence β E bevtosmyrva the reveal may observations SED m µ i eetda neovdcompo- unresolved an detected Pic OBSERVATIONS htmtyhscnre the confirmed has photometry m µ rmteeojcsis objects these from n m µ IRAS n/r70 and/or m IRAS 2 onR Najita R. Joan , dsoee ersdisks debris -discovered xess orders 2 - 1 excesses, λ o SED Boo λ µ o,H 5825 HR Boo, xessas- excesses m : 3 < Spitzer Spitzer AU 5 λ 2

(Beichman et al. 2006). The 7′′ and 21′′ angular resolu- mid- to far-infared thermal emission is gravitationally tions of Spitzer at 24 and 70 µm, respectively, are a sig- bound to the system. In this case, we can estimate the nificant improvement from the 0.75′×4.6′ and 1.5′×4.7′ minimum grain size by balancing the force due to radi- angular resolutions of IRAS at 25 and 60 µm, respec- ation pressure with the force due to gravity. For small tively, and better exclude confusion from the background grains with radius, a, the force due to radiation pressure and/or other nearby sources. Chen et al. (2006) mea- overcomes gravity for solid particles smaller than sured the shape of the hot dust continuum emitted from 3L∗ these objects from low resolution Spitzer IRS spectra and amin = (1) inferred the presence of warm dust with Tgr = 100 K - 8πGM∗cρ 150 K, hotter than estimated from photometric observa- where L∗ is the stellar luminosity. For HR 8799 with tions alone. L∗ = 6.1 L⊙ and M∗ = 1.6 M⊙ (Chen et al. 2006), we We present updated IRS spectra, extracted from S15.3 estimate amin = 1.7 and 4.7 µm, assuming that the grains pipeline products using the IRS team’s SMART pro- are either an amorphous carbon/silicate mixture or water gram (Higdon et al. 2004), overlaid on Kurucz stellar −3 2 ice (densities ρ = 2.5 and 0.91 g cm ), respectively. photosphere models that are minimum χ fit to TD1, The 150 K temperature of the warm dust indicates that Johnson UBVRI, and 2MASS stellar photometry (where this population is unlikely to be composed of water ice. available) in Figure 1. Our new reduction improves the The sublimation lifetime of 4.7 µm water ice grains is signal:noise ratio of the spectra modestly and extends re- expected to be liable continuum measurements to 38 µm, beyond which 1/2 T /T second-order light decreases the SNR. In Table 1, we list aρTgr e subl gr the fluxes of our objects in two photometric bands that tsubl = (2) have been used to search for excess from silicates (8.5-13 σ˙o µm) and cold grains (30-34 µm). The calibration un- (Jura et al. 1998) whereσ ˙o is the mass rate per surface ∼ 8 −2 −1 1/2 certainty in the fluxes is 5% and the measured sta- area (= 3.8 × 10 g cm s K , Tsubl = 5530 K; tistical uncertainties are listed in Table 1. We fit the Ford & Neufeld 2001). For HR 8799, we estimate that combined IRS, MIPS, and IRAS SEDs using two tem- 4.7 µm solid-ice grains sublimate in 40 minutes. Analy- perature black body models with the parameters listed sis of the 10 µm and 20 µm silicate features emitted by in Table 2. We estimate the temperature of the warm T-Tauri disks indicates that the warm dust in these sys- 2 component by minimum χ fitting the IRS excess at tems is predominately composed of silicates (Sargent et wavelengths shorter than 30 µm and the temperature of al. 2006). However, silicate grains with Tgr = 150 K, are the cool component by estimating the color temperature expected to produce a strong 20 µm silicate emission fea- from the remaining 30-35 µm and IRAS-60 or MIPS-70 ture that should be detected using IRS because the dust µm excesses. For HR 8799 with published submillimeter is optically thin. The lack of a 20 µm silicate emission photometry (Williams & Andrews 2006), cold dust with feature may indicate that the warm grains may be sig- β an emissivity, κν ∝ ν where β = 1, is needed to fit the nificantly larger than the minimum grain size and can be far-infrared and submillimeter SED. modeled assuming that they are large, with emissivities Direct imaging of HD 139664 and HR 8799 has pro- that are constant as a function of wavelength. Another vided additional information about their circumstellar possible reason for the lack of silicate emission features environments. High contrast Keck and Gemini AO imag- is that the grains are composed of amorphous carbon. ing has detected three planetary mass objects orbiting Amorphous carbon has a flat emissivity and therefore HR 8799 with projected separations of 24, 38, and 68 lacks spectral features. Disentangling the contribution AU and inferred masses 7, 10, and 10 MJup, inferred of carbon to this grain composition is challenging. How- from the measured K-band companion luminosities and ever, since the warm dust component of T-Tauri disks the 160 Myr central star isochronal age (Marois et al. is principally composed of silicates and icy bodies such 2008). High contrast HST ACS scattered light imaging as comets possess silicate emission features (Min et al. has resolved dust in an edge-on ring around HD 139664 2005), we believe that the composition in the HR 8799 with a possible inner edge at 60 AU, a dust peak at 83 disk is most likely silicate-rich rather than carbon-rich. AU, and a sharp outer boundary at 109 AU away from The IRAS 60 µm and SCUBA 850 µm fluxes can only the central star (Kalas et al. 2006). Kalas et al. spec- be modeled using grains with an emissivity that is in- ulate on the presence of (1) planetary embryos forming versely proportional to wavelength, suggesting that the at 83 AU that dynamically stir smaller parent bodies, cold grains are small (a << λ/2π = 9.5 µm), possibly producing a collisional cascade at this radius, (2) a plan- below the sub-blow out size of 1 µm. etary body either inside or outside of the disk that traps Black bodies (with 2πa >> λ) in radiative equilibrium migrating dust grains into mean motion resonances, also around a star with a luminosity, L∗, and grain temper- producing a dusty ring, or (3) planetary bodies inside 60 ature, Tgr, are expected to possess a distance, D, from AU and/or beyond 109 AU that confine dust into the the central star narrow ring. 1/4 L∗ 3. HR 8799 T = (3) gr 16πσD2  HR 8799 was known to possess a debris disk well be- fore the three companions orbiting the central star were For HR 8799 with Tgr(warm) = 150 K, we estimate min- discovered. We examine the grain properties to elucidate imum grain distances Din = 8.7 AU, placing the warm the planetary system architecture. dust well within the orbits of the three companions, anal- We assume that the dust detected in this system via ogous to asteroidal dust in our solar system. Small grains 3

(with 2πa << λ) in radiative equilibrium around a star = 3.6×10−5 from our fit, higher than observed toward with a luminosity, L∗, and grain temperature, Tgr, are the zodiacal dust and the dust in the Fomalhaut disk expected to possess a distance, D, from the central star (ω=0.02, Chiang et al. 2009) but similar to that of β Pic 1/5 moving group member HD 181327 (Chen et al. 2008). L∗T∗ We further estimate the ratio of the 0.6 µm scattering T = (4) gr 16πσD2  and 70 µm absorption coefficients. In our simple model, a single population of N grains with radius, a, emits a For HR 8799 with Tgr(cold) = 40 K, we estimate mini- thermal flux, mum grain distances Dout = 2000 AU, placing the cold dust well outside the orbits of the three companions, πa2 F = NQ B (T (5) analogous to Kuiper belt dust in our solar system. The ν d2 abs ν gr) warm dust temperature is higher than that of water ice sublimation, suggesting that the inner belt may possess where d (=17.5 pc) is the distance from the observer to rocky grains, while the cold dust temperature is not, sug- the central star and Qabs is the absorption coefficient. gesting that the outer belt may be icy. This same grain population is also expected to produce At the current time, there are few detailed constraints a scattered light surface brightness, on the orbits of each of the planetary-mass objects. To Lν,∗ 2 first order, the system appears face-on with circular or- SB = f(θ)Qscaπa n (6) 4πD2 bits because the stellar v sin i (= 40 km s−1) is statis- tically low for a late A-/early F-type star, the observed where n is the grain column density, f(θ) is the scatter- orbital motions for the b and c components are orthog- ing phase function for θ degrees of light deflection from onal to their position vectors (relative to the star) and forward scattering, and Lν,∗ and Qsca are the specific have magnitudes that are commensurate with that ex- stellar luminosity and scattering coefficient at the scat- pected for face-on circular orbits (Marois et al. 2008). tered light wavelength observed (0.5 µm). We estimate Using the formalism of Chiang et al. (2009), we esti- the number of dust grains in the two ansa from the dust mate that the chaotic zone associated with component column density along our line-of-sight, n, assuming a pro- b extends from 68 AU outward to ∼100 AU and that of jected area on the sky, 2∆r∆h, where ∆r and ∆h are component d extends inward from 24 AU to 12 AU. The the measured radial width and the height of the ring, lack of detected dust at distances closer to the orbit of respectively. We derive the following constraint on the the innermost putative planet may suggest that either scattering and absorption coefficients, there are additional planets in this system that have not 2 yet been detected or that the orbits of the detected plan- Qsca D Bν (Tgr)SB(gr) 2∆r∆h = 4 (7) ets are slightly eccentric, producing large regions in the Qabs d Fν,∗f(π/2) Fν (gr) disk where circumstellar grains might experience secular perturbations (Moro-Martin et al. 2007). For HD 139664 with ∆r ∼ ∆h ∼ 83 AU, SB(gr) = 11.0 −2 Reidemeister et al. (2009) have recently studied com- µJy arcsec (Kalas et al. 2006) and Fν,∗ = 52.6 Jy at −4 prehensively the properties of the dust and planets V-band, Qsca(0.5µm)/Qabs(70µm) = 7.3×10 /f(π/2). around HR 8799, reaching many of the same conclusions; Typically, Qsca/Qabs is of order unity for large grains, however, they assert that the HR 8799 IRS spectrum may suggesting that the scattering phase function at 90◦ show weak silicate emission and model the cold compo- scattering, f(π/2), must be very small. Assuming the nent using large grains. Their interpretation including a Henyey-Greenstein phase function, we can rewrite the 10 µm silicate emission feature may arise from the pres- constraint on the scattering and absorption coefficients ence of apparent absorption features in the low resolution as a function of g, the asymmetric scattering parame- spectra at 8.5 and 15 µm that are probably artifacts from ter, f(π/2) = (1 − g)2/4π(1 + g)3/2 (Henyey & Green- the edges of the SL2 and LL2 spectra. stein 1941) and search for materials that satisfy this con- straint. We calculate scattering and absorption coeffi- 4. HD 139664 cients and asymmetric scattering parameters using lab- We can infer the dust architecture in the HD 139664 oratory measurements of the optical constants for amor- system assuming that the grains are large and well- phous olivine (MgFeSiO4; Dorschner et al. 1995) solid described by black bodies. For HD 139664 with L∗ = and fluffy grains (assuming a vacuum volume fraction of 3.6 L⊙ and M∗ = 2.7 M⊙ (Chen et al. 2006), we es- 0.9) using Bruggeman effective medium and Mie theory timate amin = 1.2 µm (assuming a grain density ρs = (Bohren & Huffman 1983). We find that large (a > 40 2.5 cm−3). From our SED fit, we measure grain tem- µm), porous olivine grains can satisfy our grain proper- peratures, Tgr(warm) = 80 K and Tgr(cold) = 50 K, ties constraint while solid grains can not. Large, porous corresponding to minimum distances, Din = 21 AU and grains have been used to explain the strong forward scat- Dout = 61 AU. The similarity between the black body tering and high polarization fraction observed in the AU distance for the cold dust grains and the observed scat- Mic disk using HST ACS (Graham et al. 2007). tered light distance suggests that the thermal emission 5. and scattered light are co-spatial, originating from the λ BOO same grain population. The star λ Boo has become the archetype of a class We can infer additional grain properties by compar- of Population I late B- to early F-type stars with mod- ing the scattered light and thermal emission from the erate to extreme (up to a factor of 100) surface under- cold dust population. We estimate a Bond albedo, ω abundances of Fe-peak elements and solar abundances of = 0.3, assuming ωτ = 10−5 (Kalas et al. 2006) and τ lighter elements (C, N, O, and S) (Paunzen 2004). One 4 possible hypothesis for the abundance anomalies seen to- the asteroid and Kuiper belts in our Solar System (See ward λ Boo stars is that they selectively accrete mate- Table 2). These dusty bands apparently occur around rial from their circumstellar environments. In this sce- stars with a wide range of spectral types and ages. Some nario, grains with icy mantles and refractory cores slowly are identified using SED modeling while others are dis- spiral in under Poynting-Robertson drag, sublimating or covered via high resolution imaging. High resolution photo-desorbing water ice as they migrate inward. Once imaging has discovered unresolved excesses (consistent their volatile mantles have been depleted, their refrac- with asteroidal dust) around objects whose SEDs are tory cores are expelled via radiation pressure and the well fit using a single cold dust population, indicating gas is accreted onto the star. This hypothesis predicts that multiple dust belts may be common in debris disks. that λ Boo stars should possess debris disks with con- The Spitzer Formation and Evolution of Planetary Sys- stant surface densities, governed by Poynting-Robertson tems (FEPS) team surveyed 328 solar-like stars with −1 Drag, and SED’s with Fν ∝ ν (Jura et al. 2998). Pre- ages 0.003 - 3 Gyr quantifying the evolution of dust liminary SED analysis of 6 debris disks around A-type and gas around main sequence stars using IRAC, IRS, stars with IRS excess indicated the possibility that two λ and MIPS, finding that ∼10% solar-like stars posses cold Bootis stars (λ Boo and HR 1570) possessed disks with debris disks. Of the dusty disks, fully 1/3 are appar- −1 with Fν ∝ ν and therefore uniform surface densities ently better fit using multiple temperature components in contrast to debris disks around four normal A-type than single temperature black bodies (Hillenbrand et al. stars that possessed disks with cleared inner regions, bet- 2008), also suggesting that extended disks and multiple ter fit by single temperature black bodies (Jura et al. dust components may be common in debris disks. 2004). At the time, this observational result was con- The present day asteroid and Kuiper Belts possess an sistent with the idea that the abundance anomalies of λ estimated 0.0005 and 0.01 M⊕ in large bodies, respec- Boo stars might be generated by pollution of their stel- tively, significantly smaller than originally inferred to ex- lar atmospheres via selective accretion of volatiles from ist in these regions (∼2-3 and ∼35 M⊕, respectively). inward migrating circumstellar dust grains. Since this The migration of the outer planets and the subsequent result was published, more data have been obtained and crossing of the 2:1 resonance by Jupiter and Saturn is data reduction of spectra obtained early in the Spitzer believed to have ejected or scattered the majority of the mission has improved. initial mass in the asteroid and Kuiper Belts (Levison The connection between λ Boo stars and circumstellar et al. 2007). We estimate the mass for each debris and/or interstellar dust is highly uncertain today. Su belt observed toward λ Boo, HD 139664, and HR 8799 et al. (2008) carried out a MIPS 24 and 70 µm survey (see Table 2) assuming that the current mass in parent searching for infrared excess associated with ∼160 A-type bodies is at least as large as the amount of dust that stars, including 15 λ Boo stars (R. Gray’s λ Boo website: has been removed from the system thus far, MP B ≥ 4 2 http://www1.appstate.edu/dept/physics/spectrum/lamboo.txtLIR) tage/c (Chen & Jura 2001), assuming stellar ages, with an average age of ∼300 Myr. Their observations tage = 310, 300 (Kalas et al. 2006), and 160 Myr (Marois suggest that the 24 and 70 µm excess rate for λ Boo et al. 2008), respectively. These estimates are a lower stars is higher than average, even though the average bounds because radiation pressure (and not Poynting- age for the aggregate sample is similar ∼270 Myr. We Robertson Drag) probably dominates the grain removal estimate that 8/15 (53%) and 7/9 (78%) of the λ Boo process. We find that (1) the minimum parent body stars possess 24 and 70 µm excesses, respectively, com- masses in the warm and cold belts are 0.0002-0.003 M⊕ pared with 53/160 (33%) and 39/69 (57%) of the total and 0.002-0.004 M⊕, respectively, similar to the present sample, suggesting that the λ Boo phenomenon may be day asteroid and Kuiper belts and (2) the cold dust belts associated with circumstellar dust. Low resolution IRS are 1-10× more massive that the warm ones, suggesting spectra of 6 λ Boo stars have been published thusfar: that they have more similar masses that the present day HD 30422, HR 1570, HD 74873, HD 110411, λ Boo, asteroid and Kuiper belts. and HD 183324. New models for the IRS spectra of HD 7. 30422 and HD 110411 indicate that these objects are CONCLUSIONS better fits using a power law, consistent with a uniform We have rereduced the IRS spectra of three debris disks surface density disk with an inner radius coincident and reanalyzed their SEDs. We conclude the following: with the stellar radius (Morales et al. 2009); however, 1. The infrared SEDs of λ Boo, HD 139664, and HR revised models for λ Boo and HR 1570 (Chen et al. 8799 are well-fit using two single temperature black bod- 2006) spectra and newly published models for HD 74873 ies indicating the presence of two dust populations. For and HD 183324 (Morales et al. 2009) suggest that these HR 8799, the warmer dust is interior to the newly dis- objects are better fit using a black body, suggesting that covered orbiting planets in an exo-zodiacal dust belt and the disks around these stars possess central clearings. the cooler dust is exterior to the planets in an exo-Kuiper Therefore, the mechanism by which grains material is belt. selectively accreted may be diverse and may need to be 2. The cool grains in the HD 139664 must be strongly, revised. In any case, none of the objects presented here forward scattering, consistent with large (a > 40 µm), can be modeled using the selective accretion hypothesis porous grains if it produces both the thermal emission −1 because their SEDs are not well-fit assuming Fν ∝ ν . detected using IRAS and the scattered light disk resolved HST 6. using ACS at visual wavelengths. DISCUSSION 3. The debris disks around the majority of observed At the present time, a dozen systems are known to λ Boo stars apparently possess inner clearings similar to possess multiple dust belts, suggesting the presence of those around normal A-type stars, suggesting that sim- multiple parent body populations, perhaps analogous to ple dusty disks delivering grains from the outer regions 5

TABLE 1 Measured Fluxes

Name Spectral Distance Fν (8.5-13 µm) Fν (30-34 µm) Type (pc) (mJy) (mJy)

λ Boo A0p 30 934±5 258±7 HD 139664 F5IV-V 18 1256±7 211±5 HR 8799 A5V 40 289±3 65±6 of disks onto central stars are not responsible for the pe- Forrest, J. Graham, M. Jura, B. Macintosh, and J. Pa- culiar stellar abundances observed. tience for their helpful comments and suggestions.

We would like to thank J. Debes, M. Fitzgerald, W.

REFERENCES Backman, D. E., & Paresce, F. 1993, in Protostars and Planets Marois, C., Macintosh, B., Barman, T., Zuckerman, B., Song, I., III, eds. E. Levy and J. I. Lunine (Tuscon: University of Patience, J., Lafreniere, D., & Doyon, R. 2008, Science, 322, Arizona Press), 1253 1348 Backman, D. E., Marengo, M., Stapelfeldt, K., Su, K., Wilner, Min, M., Hovenier, J. W., de Koter, A., Waters, L. B. F. M., & D., Dowell, C. D., Watson, D., Stansberry, J. et al. 2009, ApJ, Dominik, C. 2005, Icarus 179, 158 690, 1522 Moerchen, M. M., Telesco, C. M., Packham, C., & Kehoe, T. J. J. Beichman, C. A., Bryden, G., Stapelfeldt, K. R., Gautier, T. N., 2007, ApJ, 655, L109 Grogan, K., Shao, M., Velusamy, T., Lawler, S. M. et al. 2006, Morales, F. Y., Werner, M. W., Bryden, G., Plavchan, P., ApJ, 652, 1674 Stapelfeldt, K. R., Beichman, C., Chen, Grogan, K. et al. 2009, Bohren C. F., Huffman D. R., 1983, Absorption and Scattering of ApJ, submitted Light by Small Particles (New York, NY: John Wily & Sons, Meyer, M. R., Backman, D. E., Weinberger, A. J., & Wyatt, Inc.) M. C. 2007, in Protostars and Planets V, eds. B. Reipurth, D. Chen, C. H., & Jura, M. 2001, ApJ, 166, 35, 560, L171 Jewitt, and K. Keil (Tuscon: University of Arizona Press), 573 Chen, C. H., Sargent, B. A., Bohac, C., Kim, K. H., Moro-Martin, A., Malhotra, R., Carpenter, J. M., Hillenbrand, Leibensperger, E., Jura, M., Najita, J., Forrest, W. J. et al. L. A., Wolf, S., Meyer, M. R., Hollenbach, D., Najita, J. et al. 2006, ApJS, 166, 351 2007, ApJ, 668, 1165 Chen, C. H., Fitzgerald, M. P., & Smith, P. S. 2008, ApJ, 689, 539 Reidemeister, M., Krivov, A. V., Schmidt, T. O. B., Fiedler, S., Chiang, E., Kite, E., Kalas, P., Graham, J. R., & Clampin, M. M¨uller, S., L¨ohne, T. & Neuh¨auser, R. 2009, Icarus, in press 2009, ApJ, 693, 734 Sargent, B., Forrest, W. J., D’Alessio, P., Najita, J., Li, A., Dorschner J., Begemann B., Henning Th., Jaeger C. Mutschke H., Calvet, N., Furlan, E., Green, J. D. et al. 2006, ApJ, 645, 395 1995, A&A 300, 500 Schneider, G., Weinberger, A. J., Becklin, E. E., Debes, J. H., & Ford, K. E., & Neufeld, D. A. 2001, ApJ, 557, L113 Smith, B. A. 2009, AJ, 137, 53 Graham, J. R., Kalas, P., & Matthews, B. 2007, ApJ, 654, 595 Stapelfeldt, K. R., Holmes, E. K., Chen, C., Rieke, G. H., Su, K. Jura, M., Chen, C. H., Furlan, E., Green, J., Sargent, B., Forrest, Y. L., Hines, D. C., Werner, M. W., Beichman, C. A. et al. W. J., Watson, D. M. et al. 2004, ApJS, 154, 453 2004, ApJS, 154, 458 Jura, M., Malkan, M., White, R., Telesco, C., Pina, R., & Fisher, Su, K. Y. L., Rieke, G. H., Stansberry, J. A., Bryden, G., R. S., 1998, ApJ, 505, 897 Stapelfeldt, K. R., Trilling, D. E., Muzerolle, J., Beichman, C. Kalas, P., Graham, J. R., & Clampin, M. C. 2005, Nature, 435, A. et al. 2006, ApJ, 653, 675 1067 Telesco, C. T., Fisher, R. S., Wyatt, M. C., Dermott, S. F., Kalas, P., Graham, J. R., Clampin, M. C., & Fitzgerald, M. P. Kehoe, T. J. J., Novotny, S., Marinas, N., Radomski, J. T. et 2006, ApJ, 637, L57 al. 2005, Nature, 433, 133 Kenyon, S. J., & Bromley, B. C. 2004, AJ, 127, 513 Wahhaj, Z., Koerner, D. W., Ressler, M. E., Werner, M. W., Koerner, D. W., Ressler, M. E., Werner, M. W., & Backman, Backman, D. E., & Sargent, A. I.2003, ApJ, 584, L27 D. E. 1998, ApJ, 503, L83 Williams, J. P., & Andrews, S. M. 2006, ApJ, 653, 1480 Levison, H. F., Morbidelli, A., Gomes, R., & Backman, D. 2007, Wyatt, M. C., Greaves, J. S., Dent, W. R. F., & Coulson, I. M. in Protostars and Planets V, eds. B. Reipurth, D. Jewitt, and 2005, ApJ, 620, 492 K. Keil (Tucson: University of Arizona Press), 669 Wyatt, M. C., Smith, R., Greaves, J. S., Beichman, C. A., Lisse, C. M., Chen, C. H., Wyatt, M. C., & Morlok, A. 2008, Bryden, G., & Lisse, C. M. 2007, ApJ, 658, 569 ApJ, 673, 1106 Liu, W. M., Hinz, P. M., Hoffman, W. F., Brusa, G., Wildi, F., Miller, D., Lloyd-Hart, M., Kenworthy, M. A. et al. 2004, ApJ, 610, L125 6

TABLE 2 Dust Properties

Warm Dust Cold Dust Name β Tgr D LIR/L∗ MP B β Tgr D LIR/L∗ MP B (K) (AU) M⊕ (K) (AU) M⊕

λ Boo 0 110 27 5.8×10−6 0.003 0 60 115 7.6×10−6 0.004 HD 139664 0 80 21 4.2×10−6 0.0004 0 50 61 2.4×10−5 0.003 HR 8799 0 150 9 1.9×10−6 0.0002 1 40 2000 2.5×10−5 0.002

TABLE 3 Multiple Dust Belt Systems

HR HD Name Spectral Age Dhot Dcold References Type (Myr) (AU) (AU)

1084 22049 ǫ Eri K2V 850 ∼3 35-90 1 2020 39060 β Pic A6V 12 <5 100-1500 10, 11 3927 86087 A0V 50 ∼7 ∼50 2 4775 109085 η Crv F2V 1000 ∼2 130-170 2, 12 4796 109573 A0V 8 <20 60-87 5, 8 113766 F3/F5V 16 ∼4 ∼30-80 2 5351 125162 λ Boo A0p 310 ∼24 ∼115 this work 139664 F5IV-V 200 ∼15 60-110 4, this work 7329 181296 η Tel A0Vn 12 ∼5 ∼30 2 8728 216956 Fomalhaut A4V 200 <10 130-160 3, 9 8799 218396 A5V 160 ∼8 ∼2000 7, this work

References. — (1) Backman et al. 2009; (2) Chen et al. 2006; (3) Kalas et al. 2005; (4) Kalas et al. 2006; (5) Koerner et al. 1998; (6) Lisse et al. 2008; (7) Marois et al. 2008; (8) Schneider et al. 2009; (9) Stapelfeldt et al. 2004; (10) Telesco et al. 2005; (11) Wahhaj et al. 2003; (12) Wyatt et al. 2005 7

Fig. 1.— SEDs for λ Boo, HD 139664, and HR 8799: TD1 fluxes are plotted with magenta error bars; General Catalogue of Photometric Data mean UBV or Johnson et al. fluxes are plotted with yellow error bars; and 2MASS JHK fluxes are plotted with cyan error bars. Color-corrected IRS, MIPS, and IRAS data, where available, are shown with red, green, and blue error bars. 8

2 Fig. 2.— Photosphere-subtracted IRS spectra with Fν plotted as a function of wavelength. The minimum χ fits for the two single temperature black body component are overlaid with the warm dust component shown as a dashed line, the cold dust component shown as a dashed-dotted line, and the sum of the two component as a solid line.