arXiv:0710.4941v1 [astro-ph] 25 Oct 2007 7MsahstsAe,Cmrde A 02139 MA, Cambridge, Ave., Massachusetts 77 aseolnt e eetd( detected lowest- the yet . to for different exoplanets relationship lead on mass-radius mass can masses the and atmosphere explore We of exist range mechanisms a loss and formation 2005). al. et (Sato envelope ascuet nttt fTcnlg,7 ascuet A Massachusetts 77 Technology, 02139 MA, of Cambridge, Institute Massachusetts t necnrcobt;adtegsrc yrdplanet hybrid gas-rock plan- the its giant and with of 149026b orbits; predominance HD eccentric exis- the in the Jupiters; ets include: hot examples 2007). of al. Notable tence et Sotin before. 2007; sumptions al. et 2007; Selsis al. 2007; et al. Fortney et 2006a,b; al. Seager plan- et (Valencia at- solid 10 hydrogen-rich to of that 5 thick assumed the a have in ra- ets with authors and Most mass rocky interpretation mosphere. a alternative the An be on planet. could based transiting a identified of dius uniquely be not mass potentially could measured planet a water identified. From low-density be a transits. radius through and is star worlds the mi- ter near for is evidence exist world strong volatiles for be water insufficient would a it because star, If gration a to 2004). water close found solid al. L´eger et of 2003; mostly (Kuchner composed worlds”, “water composition. bulk and prop- the mass planetary as other precise such determine to erties, the are begin derive planets can Transiting and can radius, we stars. because their transit important to known are 2 eaotteie htawd ag fatmospheric of range wide a that idea the adopt We 1 xpaeshv,hwvr otaitdorbscas- basic our contradicted however, have, Exoplanets can- worlds water that possibility the examine We or planets” “ocean to given been has attention Much 20 over date, to known exoplanets 250 over of Out rpittpstuigL using 2018 typeset 5, Preprint September version Draft eateto hsc,MsahstsIsiueo Technol of Institute Massachusetts Physics, of Department nsitu in eateto at,Amshrc n lntr Sciences, Planetary and Atmospheric, Earth, of Department CA LNTO HC TOPEE NTEMS-AISRELATI MASS-RADIUS THE ON ATMOSPHERE: THICK OR PLANET OCEAN aiso J46,atog l oesrqiesm /e econclud impossible We means. is other headings: H/He. alone by compositio Subject relationship some out of mass-radius ruled require variety the be wide models can on a all based example, although worlds For 436b, water pla He. GJ planet A and/or of ocean H radius value. so-called o some gas-free (a 0.1%-10% and content its to ice above core equivalent water 5-60% envelope substantial increase have gas to might a radius adding the rang that dens o causes average find the difficult same We are within the parameters—which positions. with unique planetary planets not important alone—produce is measurements of number density a average of the however, planets, h ukcmoiino neolnti omnyifre rmisave its from inferred commonly is an of composition bulk The omto.Tepooe dnicto fwa- of identification proposed The formation. M 1. ⊕ A ag aea ninfiatamount insignificant an have range INTRODUCTION T ∼ E tl mltajv 12/01/06 v. emulateapj style X 60 xrslrplanets extrasolar M ⊕ oeand core XPAESWT ASV ATMOSPHERES MASSIVE WITH EXOPLANETS ∼ .R Adams R. E. 5-20 ∼ M 30 ⊕ rf eso etme ,2018 5, September version Draft norder in ) M 1 .Seager S. , ⊕ H/He ABSTRACT ogy, ve., 1,2 .Elkins-Tanton L. , pei as to mass) spheric rsneo asv topee eepoeatmo- explore We the from from atmosphere. result ranging massive spheres that a ambiguities of presence potential identify to iinsre ydniy eslefrtemass the for solve We compo- different density. by of sorted shells sition concentric with equilibrium, pressure rniigNpuesz lntG 3b(ulre al. et (Butler the 2007b). 436b al. for GJ et compositions Gillon mass planet 2004; potential smaller Neptune-size the explore on also transiting focus We a 1995; with al. range. 1991)), et al. Podolak et 2005; Hubbard (Guillot H/He Neptune’s and nertdotad and outward, integrated hr h topeecae ob pqea iil and visible at opaque be wavelengths to infrared ceases atmosphere the where ino tt o h aeil ttebudr between boundary the At material. the layers, for state of tion o h etlyr h nerto ssopda the at stopped is integration The condition boundary layer. outer next the for n h qaino tt ftematerial. the of state of equation the and rlpesr dutdt civ h eie as the mass, cen- desired a the from for achieve Starting equations to differential specified. adjusted are pressure layers three tral to up of equilibrium, hydrostatic of equation the ue ai fetaoa planets. extrasolar of radii sured h qainfrms fashrclshell, spherical a of mass for equation the 3 u oe sue peia lnti hydrostatic in planet spherical a assumes model Our obidamdlpae,tecmoiinadmass and composition the planet, model a build To einr n aeeghdpneteet nptnilmea potential on effects wavelength-dependent any ignore We P P ( r ( and ) r ,density ), ratraieyalrerocky- large a alternatively or ) te u ieyvrigbl com- bulk varying widely but ities e ihagvnms n radius and mass given a with net dP ∼ m fcmoiin.Variations compositions. of e mosbet osri from constrain to impossible r ρ dr dm ( nesasgicn a layer gas significant a unless ( 1 1 ( r scnepanteobserved the explain can ns r h aso oi planet solid a of mass the f r dr M 3 = ) htteietfiainof identification the that e r sda h nta conditions initial the as used are ) ) . ( 2. r ∼ ⊕ ρ = ) aedniy o small For density. rage ( ρ f r teetmtdms fUranus’ of mass estimated (the P MODELS ( 4 = − 10 ,adtemperature and ), ( r P T ( scluae rmteequa- the from calculated is ) Gm r − ( ( and ) NHPFRSOLID FOR ONSHIP πr r r 3 a,approximately bar, 1 = ) M ) P , ( 2 r r ρ ⊕ 2 ) ( ( ρ r m r (10 ( ) )) r , ( ) . r , r numerically are ) × eu’atmo- Venus’ T ( r using ) m ( (2) (1) (3) r ), - 2

The temperature of a planet’s gas layer has a sig- 1987), and H2O are described in more detail in nificant effect on the radius. Instead of computing a Seager et al. (2007). cooling history to obtain T (r), we calculate T (r), sep- arately in three different regimes. Within the solid por- 3. RESULTS AND DISCUSSION tion of the planet (r ≥ rsolid), temperature has little effect on the final radius, and is assumed to be isother- Our fiducial planet consists of a 30% Fe core and a 70% mal (see Seager et al. 2007). In the deep hydrogen- MgSiO3 mantle, roughly analogous to Earth. We used layer where the pressure is greater than 1 kbar a H/He mixture with helium mass fraction Y = 0.28 (the He mass fraction of the solar nebula). We chose (rsolid ≤ r ≤ r1kbar), we assume that convection domi- nates and the temperature follows an adiabat tied to the Teq = 300 K, based on the observation that a planet around an M-dwarf at an orbital distance of 0.1 AU has entropy, S, at 1 kbar. At lower pressures (r > r1kbar), we use the radiative equilibrium gray analytical model of a similar equilibrium temperature to Earth (assuming Chevallier et al. (2007) and Hansen (2007) for irradiated similar albedos). We set Teff = 30 K, similar to Earth and Uranus4. For the atmospheric parameters, we fixed atmospheres, ◦ µ0 = cos60 and γ =0.1 to represent radiation absorbed 3 4 2 4 deep in the atmosphere5. We later investigate variations T (r)= Teff τ + + TeqF (τ,µ0,γ), (4) 4  3 on Y , Teff , Teq and γ. Figure 1 shows a plot of the mass-radius relationship for fiducial planets of masses 5, where τ is the optical depth, µ0 = cos θ0 where θ0 is the angle of incident radiation with respect to the surface 10, 15, and 20 M⊕. For each planet mass, we added normal, and γ is a parameter that accounts for the al- atmospheres ranging in mass from 0.001-1 M⊕. titude at which radiation is absorbed. We convert from A robust finding for all models is that a small amount optical depth to pressure through a constant scaling re- of gas creates a large radius increase. While this result is lation under the assumption of hydrostatic equilibrium expected, the radius increase is far more dramatic than and constant opacity per gram, following Hansen (2007) anticipated. but taking into account the planet’s surface gravity. Here For example, an H/He atmosphere of ∼ 0.001 by mass—only ten times greater than Venus’ atmospheric Teff is the effective temperature in the absence of stel- lar irradiation, representing energy from internal sources. mass fraction—is required for a noticeable radius in- crease. As seen in Figure 1, adding a hydrogen-helium Teq is the equilibrium temperature which represents heat- ing from the parent star. atmosphere with just 0.1% of the mass of a 10 M⊕ rocky We emphasize that by choosing an effective planet results in a 5% increase in the planetary radius— temperature—or adiabat—we are subsuming a cooling within a measurement precision that has been obtained calculation. This is common practice for modeling the for currently known transiting planets. interior structure of giant planets (e.g., As a second example, adding a gas layer of H/He equal Stevenson 1982; Hubbard et al. 1991; Marley et al. 1995; to 1% of the mass of our fiducial planets increases the Saumon & Guillot 2004) and was also the case even be- radius by ∼ 20% of the original planet radius, or by about 0.35 R⊕ for the planet masses we considered. fore their Teff s and gravitational moments were known or well known (Demarcus 1958; Zapolsky & Salpeter Our major finding is that that exoplanets with a sig- nificant H/He layer cannot be distinguished from water 1969). In fact, even with a known Teff cooling calculations only match observed parameters for worlds, based on Mp and Rp alone. For our fiducial solid Jupiter (e.g, Hubbard 1977); Saturn is hotter than ex- exoplanets, adding up to 5% H/He by mass (for 10 M⊕ pected (Pollack et al. 1977; Stevenson & Salpeter 1977; planets) is sufficient to push the planet’s radius through Fortney & Hubbard 2003) and Uranus and Neptune are the entire range of radii corresponding to solid planets colder than expected (Stevenson 1982; Podolak et al. with no gas, including planets with up to 100 percent 1991; Guillot 2005). An additional motivation to water composition. While we have not completed an ex- adopt a simple framework for temperature is the huge haustive study of possible compositions, we find the non- uncertainty in interior parameters for rocky planets. uniqueness of water planets to be valid for any conditions This includes the temperature-dependent equation of we investigated. state of liquids and solids at high temperature and This generic finding holds for a wide range of as- pressure, the temperature-dependent viscosity, and sumptions of assumed temperatures. Taking our fiducial the effect of tides on the cooling history for eccentric model, we vary Teq, Teff , and γ individually. For a 10 exoplanets. We further note that our choice of radiative M⊕ solid planet with an additional 0.1 M⊕ H/He atmo- equilibrium down to 1 kbar is based on the irradiated sphere, increasing Teq from 300 to 500 K increases the atmosphere/interior models by Fortney et al. (2007). radius by about 1% (Figure 2). For the same planet varying Teff from 10 K to 50 K results in an 8% increase For our fiducial model we choose Teff to be similar to Earth’s and Uranus’s. While the largest uncertainty in radius. While large, this value is comparable to ex- pected radius uncertanties for these planets (Gillon et al. in our treatment is the qualitative choice of Teff , we subsequently vary over a reasonable range of effective 2007b,a; Deming et al. 2007). Varying the altitude where temperatures (or adiabats). radiation is absorbed (specified by γ) has a much smaller We use the H/He equation of state from Saumon et al. effect on the planet radius. Varying γ from 0.1 to 10 (1995), ignoring the “plasma-phase transition” which causes the radius to decrease by 0.2%. may be a numerical artifact (D. Saumon 2007, private 4 12 Earth has 44 ×10 W (Pollack et al. 1993) and Uranus has communication). The equations of state for the solid 340 × 1012 W of energy flow (Pearl et al. 1990). materials Fe (Anderson et al. 2001), MgSiO3 perovskite 5 In comparison γ = 10 would correspond to absorption high in (Karki et al. 2000), (Mg,Fe)SiO3 (Knittle & Jeanloz the atmosphere. 3

5 MÅ core 10 MÅ core A corollary of our main result is that when a planet has 3. 3. a significant H/He atmosphere there is a wide degener- water acy in allowable internal composition. This is not just 2.5 2.5

Å Å compositional, but also relates to the trade off of tem- R water R perature and mass of H/He gas. It could be argued that rock 2. 2. specifying a planet’s composition implies a particular in- rock ternal thermal profile derived from a consistent cooling history. As addressed in §2, the many unknowns and 5. 5.5 6. 10. 10.5 11. MÅ MÅ free input parameters for rocky planet interiors—such as the possible differences between atmospheric and interior 15 MÅ core 20 MÅ core water compositions, equation of states, and the effect of tides 3. 3. water on the planet’s cooling history—prevent a self-consistent solution for the present time. 2.5 2.5 Å Å rock How could a 5–20 M⊕ exoplanet get a substantial R rock R H/He layer? Two different scenarios may produce them: 2. 2. direct capture of gas from the (possi- iron iron bly modified by the escape of some fraction of the original 15. 15.5 16. 20. 20.5 21. gas) or outgassing during . MÅ MÅ A planet may capture and retain up to 1 to 2 M⊕ of H/He if the planetary core did not grow quickly enough Fig. 1.— The increase in radius due to adding H/He to a solid to capture more before the gas in the disk evaporated planet. A H/He layer of 0.002-1 M⊕ is added to a solid planet of 5, 10, 15, or 20 M⊕, with fiducial model parameters (30% Fe and 70% (as is the paradigm for Uranus and Neptune). Alterna- MgSiO3). The black points are for atmospheres at 0.01 M⊕ and tively, for short-period exoplanets, a 1 to 2 M⊕ H/He every 0.1 M⊕ afterwards. The mass-radius relationship of solid envelope may result after substantial loss of an initially planets with no gas is plotted for comparison. The water (blue), rock (red), and iron (green) curves are taken from Seager et al. massive gas envelope from irradiated evaporation (e.g. (2007) and represent homogeneous solid planets. Intermediate Baraffe et al. 2006). Alibert et al. (2006) consider atmo- compositions for differentiated planets are, from top down: dashed- spheric evaporation during migration, and conclude that blue, 75% H2O, 22% MgSiO3, and 3% Fe; dashed-dotted-blue, 48% the 10M⊕ innermost planet in HD 69830, at 0.08 AU, H2O, 48.5% MgSiO3, and 6.5% Fe; dotted-blue, 25% H2O, 52.5% MgSiO3, and 22.5% Fe; dashed-red, 67.5% MgSiO3 and 32.5% Fe; kept ∼2 M⊕ of H/He over the 4 Gyr lifetime of the star. and dotted-red, 30% MgSiO3 and 70% Fe. In general, the addition Little attention has been given to the mass and compo- of a gas layer of up to ∼ 5% of the solid planet mass will inflate sition of exoplanet atmospheres from outgassing. Venus’ the radius of a rocky-iron planet through the range of radii cor- −4M responding to water planets with different water mass fractions. atmosphere is 10 ⊕; if Venus had a surface gravity high enough to prevent H escape, its atmosphere would −3 be over 10 M⊕. Even more massive H-rich atmospheres are possible. If a massive iron-silicate planet formed with Teq H300, 400 and 500 KL enough water, the iron may react with the water dur-

2.75 ing differentiation, liberating hydrogen gas (Ringwood 1979; Waenke & Dreibus 1994). L. Elkins-Tanton et al. 2.5 (in prep.) estimate that the maximum H component is Å R 2.25 about six percent by mass for a terrestrial-composition planet. For a 10M⊕ planet this would result in a 0.6M⊕ 2. H envelope. For short-period, low-mass planets, theoretical argu- 10. 10.1 10.2 10.3 10.4 10.5 MÅ ments of atmospheric escape may be the best way to identify a water world based on the mass and radius mea- Teff H10, 30 and 50 KL surements alone (L´eger et al. 2004; Selsis et al. 2007). 2.75 Indeed, our assumption of H/He atmospheres for exo- planets relies on the condition that atmospheric mass loss 2.5

Å has not evaporated all of the H/He. In the absence of R 2.25 hydrodynamic escape, the exospheric temperature (and not the atmospheric T ) drives the thermal Jeans es- 2. eff cape of light gases. Earth and Jupiter both have exobase

10. 10.1 10.2 10.3 10.4 10.5 temperatures of 1000 K (de Pater & Lissauer 2001), sig- Å M nificantly above their Teff of 255 K and 124.4 K respec- tively (Cox 2000). Uranus and Neptune have exobase Fig. 2.— The effects on the radius of varying the equilibrium temperature (top) and effective temperature (bottom) for a 10 M⊕ temperatures around 750 K (de Pater & Lissauer 2001). planet with otherwise fiducial parameters. Teq values of 300, 400, We note that because GJ 436 must have at least 1 M⊕ and 500 K are plotted to simulate the effect of uncertainty in the of H/He, its exospheric temperature is not too high. On orbital parameters and albedo on the expected radius. T values eff the other extreme, planets of 5M⊕ would require very low of 10, 30, and 50 K are plotted on the same scale as the Teq plot, to show uncertainties in the planet’s interior temperature. Uncer- exospheric temperatures (∼ 300 K) to retain a massive tainties in the internal energy of a planet lead to large variations atmosphere over the course of billions of years. Never- in radii for a given mass, showing how temperature is a large un- theless, a young 5M⊕ Earth-mass planet with a captured certainty in the interpretation of a planet’s internal composition. atmosphere could still have a H/He atmosphere and an old 5M⊕ planet could retain a substantial He fraction, 4

35 so capture must be at least partially invoked to explain

19 MÅ solid H30% Fe, 70% MgSiO3 L such a massive atmosphere. The second composition for 30 3.2 MÅ HHe GJ 436b we considered is for water worlds, one with 25 a 50% water mantle (by mass) and 50% silicate core, L and another with 90% water mantle and a 10% silicate cm^3

 20 g

H core. These planets also need some H/He to match the

15 known radius, 12% and 8% by mass, respectively. The

Density third model approximates planets with atmospheres cre- 10 20 MÅ solid H100% MgFeSiO3 L ated from outgassing, considering an extreme scenario 2.2 MÅ H Å where all of the available water has oxidized iron, leav- 5 20.5 M solid H90% H2O, 10% MgSiO3 L ing a 100% (Mg,Fe)SiO3 solid planet core. To match 1.7 MÅ HHe 0 the observed radius, a 22.2 M⊕ planet requires ∼ 2.2M⊕ 0 1. 2. 3. 4. of H alone, a case that assumes no initial trapping and H ÅL Radius R subsequent outgassing of He. The pure-hydrogen atmo- Fig. 3.— Density vs. radius for three different potential compo- sphere is 10% of the mass of the solid planet. This is sitions of GJ436b. From top to bottom: solid (black) curve, 19.0 M⊕ core (30% Fe, 70% MgSiO3) with 3.2 M⊕ H/He (Y=0.28); above the theoretical maximum of outgassing based on dotted (red) curve, 20.0 M⊕ core (100% (Mg,Fe)SiO3) with 2.2 observed abundances of metallic iron in chondritic mete- M⊕ H (Y=0); dashed (blue) curve, 20.5 M⊕ core (90% H2O, 10% orites from our solar system (see Elkins-Tanton et al. in MgSiO3) with 1.7 M⊕ H/He (Y=0.28). All three planets have the prep). Although not an exhaustive study, the range of in- same total radius (4.3 R⊕) and total mass (22.2 M⊕). terior compositions illustrates the variety of possibilities, making its compositional identification ambiguous. though all models require some H/He. We now turn to a qualitative study of GJ436b, to show While our study is preliminary, we make the robust that the interpretation that GJ436b is a water world point that H-rich thick atmospheres will confuse the in- akin to Uranus and Neptune (Gillon et al. 2007a) is not terpretation of planets based on a measured mass and the only possibility. We consider the GJ 436b values radius. This point is independent of the uncertainties re- Mp = 22.2M⊕ and Rp = 4.3R⊕ from Deming et al. tained by our model including Teq, Teff , the mass frac- (2007). The internal structure in Figure 3 shows how tion of H/He, and the mixing ratio of H and He. We three planets with very different internal compositions find that the identification of water worlds based on the can have the same total mass and radius. We first ex- mass-radius relationship alone is impossible unless a sig- plore a planet similar to our fiducial model: a 22.2 M⊕ nificant gas layer can be ruled out by other means. Spec- solid planet with Earth-like iron/rock mass ratio (30/70), troscopy is the likely means and may become rou- Teff = 30 K, and Teq = 600 K in rough agreement tine with transit transmission and emission spectroscopy, with the orbital parameters (assuming an albedo of 0.1). and eventually with spectroscopy by direct imaging. By adding ∼3.2M⊕ of hydrogen-helium to the 19.0 M⊕ solid planet, we are able to reproduce GJ 436b’s radius. We note that the mass of gas is 15% of the solid mass, We thank Mark Marley, Jonathan Fortney, and Wade likely too much to have originated from outgassing, and Henning for useful discussions.

REFERENCES Alibert, Y., Baraffe, I., Benz, W., Chabrier, G., Mordasini, C., Gillon, M., Pont, F., Demory, B., Mallmann, F., Mayor, M., Mazeh, Lovis, C., Mayor, M., Pepe, F., Bouchy, F., Queloz, D., & Udry, T., Queloz, D., Shporer, A., Udry, S., & Vuissoz, C. 2007b, ArXiv S. 2006, A&A, 455, L25 e-prints, 705 Anderson, O. L., Dubrovinsky, L., Saxena, S. K., & LeBihan, T. Guillot, T. 2005, Annual Review of Earth and Planetary Sciences, 2001, Geophys. Res. Lett., 28, 399 33, 493 Butler, R. P., Vogt, S. S., Marcy, G. W., Fischer, D. A., Wright, Hansen, B. M. S. 2007, submitted to ApJS J. T., Henry, G. W., Laughlin, G., & Lissauer, J. J. 2004, ApJ, Hubbard, W. B. 1977, Icarus, 30, 305 617, 580 Hubbard, W. B., Nellis, W. J., Mitchell, A. C., Holmes, N. C., Baraffe, I., Alibert, Y., Chabrier, G., & Benz, W. 2006, A&A, 450, McCandless, P. C., & Limaye, S. S. 1991, Science, 253, 648 1221 Karki, B. B., Wentzcovitch, R. M., de Gironcoli, S., & Baroni, S. Chevallier, L., Pelkowski, J., & Rutily, B. 2007, Journal of 2000, Phys. Rev. B, 62, 14750 Quantitative Spectroscopy and Radiative Transfer, 104, 357 Knittle, E. & Jeanloz, R. 1987, Science, 235, 668 Cox, A. N. 2000, Allen’s astrophysical quantities (Allen’s Kuchner, M. J. 2003, ApJ, 596, L105 astrophysical quantities, 4th ed. Publisher: New York: AIP Press; L´eger, A., Selsis, F., Sotin, C., Guillot, T., Despois, D., Mawet, D., Springer, 2000. Edited by Arthur N. Cox. ISBN: 0387987460) Ollivier, M., Lab`eque, A., Valette, C., Brachet, F., Chazelas, B., de Pater, I. & Lissauer, J. J. 2001, Planetary Sciences (Planetary & Lammer, H. 2004, Icarus, 169, 499 Sciences, by Imke de Pater and Jack J. Lissauer, pp. 544. ISBN Marley, M. S., G´omez, P., & Podolak, M. 1995, J. Geophys. Res., 0521482194. Cambridge, UK: Cambridge University Press, 100, 23349 December 2001.) Pearl, J. C., Conrath, B. J., Hanel, R. A., & Pirraglia, J. A. 1990, Demarcus, W. C. 1958, AJ, 63, 2 Icarus, 84, 12 Deming, D., Harrington, J., Laughlin, G., Seager, S., Navarro, Podolak, M., Weizman, A., & Marley, M. 1995, Planet. Space Sci., S. B., Bowman, W. C., & Horning, K. 2007, ArXiv e-prints, 43, 1517 707 Podolak, M., Hubbard, W. B., & Stevenson, D. J. Model of Uranus’ Fortney, J. J. & Hubbard, W. B. 2003, Icarus, 164, 228 interior and magnetic field (Uranus), 29–61 Fortney, J. J., Marley, M. S., & Barnes, J. W. 2007, ApJ, 659, 1661 Pollack, H. N., Hurter, S. J., & Johnson, J. R. 1993, Reviews of Gillon, M., Demory, B.-O., Barman, T., Bonfils, X., Mazeh, T., Geophysics, 31, 267 Pont, F., Udry, S., Mayor, M., & Queloz, D. 2007a, A&A, 471, Pollack, J. B., Grossman, A. S., Moore, R., & Graboske, Jr., H. C. L51 1977, Icarus, 30, 111 Ringwood, A. E. 1979, Origin of the earth and moon (New York, Springer-Verlag New York, Inc., 1979. 307 p.) 5

Sato, B., Fischer, D. A., Henry, G. W., Laughlin, G., Butler, R. P., Sotin, C., Grasset, O., & Mocquet, A. 2007, Icarus in press Marcy, G. W., Vogt, S. S., Bodenheimer, P., Ida, S., Toyota, E., Stevenson, D. J. 1982, Planet. Space Sci., 30, 755 Wolf, A., Valenti, J. A., Boyd, L. J., Johnson, J. A., Wright, Stevenson, D. J. & Salpeter, E. E. 1977, ApJS, 35, 239 J. T., Ammons, M., Robinson, S., Strader, J., McCarthy, C., Valencia, D., O’Connell, R. J., & Sasselov, D. 2006a, Icarus, 181, Tah, K. L., & Minniti, D. 2005, ApJ, 633, 465 545 Saumon, D., Chabrier, G., & van Horn, H. M. 1995, ApJS, 99, 713 Valencia, D., Sasselov, D. D., & O’Connell, R. J. 2006b, ArXiv Saumon, D. & Guillot, T. 2004, ApJ, 609, 1170 Astrophysics e-prints Seager, S., Kuchner, M., Hier-Majumder, C. A., & Militzer, B. Waenke, H. & Dreibus, G. 1994, in Deep Earth and Planetary 2007, ApJ, in press Volatiles, 46–+ Selsis, F., Chazelas, B., Borde, P., Ollivier, M., Brachet, F., Zapolsky, H. S. & Salpeter, E. E. 1969, ApJ, 158, 809 Decaudin, M., Bouchy, F., Ehrenreich, D., Griessmeier, J. ., Lammer, H., Sotin, C., Grasset, O., Moutou, C., Barge, P., Deleuil, M., Mawet, D., Despois, D., Kasting, J. F., & Leger, A. 2007, ArXiv Astrophysics e-prints