sign in sign up
<<
Home , Hot Jupiter

arXiv:astro-ph/0407620v1 29 Jul 2004 rpittpstuigL using typeset Preprint l o ial ersra lnt ofr Amtg 2003) (Armitage form impossi to it terrestrial making sizable re-form, for to ble later, time strongl enough occurs not without migration is destroyed million if populati planets) However, first terrestrial the the of pleted. within lifetime, blocks occurs disk building the (the migration of this (Myr) If years Richard- & Bodenheimer 1996). (Lin, gravitation disk son via gas massive inward a migrated with and torques out farther formed have oa ytmhv riswti . srnmcluis(AU) units astronomical central 0.5 their within of orbits have System Solar fapae riigti eyyugsa (Forrest young presence very the this is orbiting explanation One a AU. an revealed 10 of have inside Telescope dust Space of Spitzer absence the with 4 Tau Coku nie01A.I htnmnltr,oemgtcl lnt a planets call might one nomenclature, that In AU. 0.1 inside Bieoe l20) bevtoso the of Observations with d 2001). agreement circumstellar in al of less, lifetime et or year (Briceño Myr million 1-10 1 observed, in can the planets mechanism giant this that via shown have form 2004) Benz (Alib & formation Mordasini during migration ert, and 2003) Armitage & (Rice ftegatpae (Lufkin planet 2003). giant forma Henning the the of during, & even Kley, after, immediately (D’Angelo, begins planets Migration mass - to ed ntems ftepae n h asadvsoiyof viscosity and 10 mass than less the typically and is planet and disk, the gaseous of the mass the on pends gi fast for evidence observational be formation. planet would this correct, If crto cnro(Pollack scenario Rice 2002; hc irto tp sntwl nesod n a involv may and understood, well not is stops migration which ieclsvagaiainlclas Bs 97 Mayer 1997; (Boss collapse gravitational via timescales D ATVERSION RAFT ogl n hr ftegatpaesdsoee usd th outside discovered planets giant the of third one Roughly h iecl o h nadmgaino in lntde- planet giant a of migration inward the for timescale The eetrslsso htgatpaescnfr nvr short very on form can planets giant that show results Recent 4 3 2 1 euetetr o uie oapyt a in lnt ihs with planets giant gas to apply to jupiter hot term the use We http://www.exoplanets.org from Data ua n lntr aoaoy h nvriyo Arizona, of University The Laboratory, Planetary and Lunar eateto srnm,Uiest fWsigo,Bx351 Box Washington, of University Astronomy, of Department ujc headings: e Subject may belts . th that hot show speculate terrestria We We the of timescales. content suppressed. water delivery is the formation enhance potential not their AU) does greater 0.25 Jupiter conte 0.15, or water AU (e.g. have 0.5 jupiters can of hot and of jupiters, radii hot orbital with stars pla around terrestrial form of radii orbital simulations having of as results defined lifetime a present the If We in re-generate disk. 2003). to gas time massive has a population with etesimal interactions via inward migrate tal. et H OMTO N AIAIIYO ERSRA LNT NTH IN PLANETS TERRESTRIAL OF HABITABILITY AND FORMATION THE Htjptr, in lnt ihobt eycoet thei to close very orbits with planets giant jupiters,’ ‘Hot A UGUST A 03.Nwsmltoso h tnad core- standard, the of simulations New 2003). , T E tl mltajv 04/03/99 v. emulateapj style X 0 2018 30, 3 hs htjupiters” “hot These . 1. INTRODUCTION tal. et tal. et lntr omto xrslrpaes–cosmochemistry – planets extrasolar – formation planetary 04.Temcaimby mechanism The 2004). , 96 nldn turbulence including 1996) , S EAN ≤ .R N. . U eso httretilpaessmlrt hs nt in those to similar planets terrestrial that show We AU. 0.5 ∼ 5 AYMOND 4 er o - for years y l star old Myr 1 r huh to thought are . U“amjupiters.” “warm AU 0.5 t tal. et usn Z827 ([email protected]) 85287. AZ Tucson, rf eso uut3,2018 30, August version Draft mmjraxes emimajor 1 8,Sate A915([email protected] tr ([email protected]; 98195 WA Seattle, 580, T , 2004). , tal. et HOMAS . de- y isks tion ABSTRACT ant are JUPITERS on al e e - - , Q 1 a UINN ≤ lnt,btrte erae hi omto n water and formation their decreases rather but planets, l yhbtbepaescnfr,btfrsm-ao axes semi-major for but form, can planets habitable ly fteds n ersra lnt a om(Armitage form may planets terrestrial and disk the of e omto ntepeec fhtjptr,broadly jupiters, hot of presence the in formation net so ie”weewtrcudeita c ntelwpressur low the the past in ice belt, (Morbidelli as disk asteroid exist protoplanetary outer could formation water the where during in line,” water “snow formed its which of bodies from acquired Earth The 2004). t ihdfeetcmoiin seRaymond (see compositions different pr with proto- of habita stirring ets a the radial of significant in formation requires The out therefore line. planet T farther snow the much past content. lies disk, water planetary however, substantial water, a of needs source its also in it reside A to zone, needs AU 1993). habitable only 1.37 Reynolds not & - planet Whitmire th 0.95 habitable (Kasting, potentially roughly on System is exist Solar and to our planet, water in Earth-like liquid an for of right surface is temperature the which 04 n h rsneo n rmr in lnt (Cham- Raymond planets 2003; giant Agnor more & Levison or 2002; one Cassen of & bers presence Lunine the & and Quinn Raymond, 2004) 2002; (Wet Cassen material & solid Chambers of 1996; amount ill the by affected strongly is ets iin n ueia ehd.Scin3peet u resul our presents 3 cases Section initial some our methods. in outlines numerical 2 and and Section AU, AU. jupiter ditions 0.5 5.2 hot and at planets 0.25 include giant We 0.15, outer of planet. radii orbital giant both jupiter, with exterior hot an a of without presence and the in formation planet trial discussion). a lnt a omi tnad otmu aho ntepres the in fashion jupiter. hot bottom-up a standard, of a ence in form can 2000). Papaloizou planets & (Nelson star Mat- the 2002; fact into in way Lecar may the & planets all Many migrate (Kuchner 2003). disk Murray evac- & inner Johnstone, an suyama, the or in 2003) region (Terquem uated fields magnetic with interactions . U hsdprsfo eti sso h em hc rese which term, the of uses certain from departs This AU. 0.5 ititro otetretilpaesi ytm with systems in planets terrestrial the to interior xist aetsas r huh ofr ate wyand away farther form to thought are stars, parent r 1 h aial oeaon tri enda h nuu in annulus the as defined is star a around zone habitable The eew rsn eut fdnmclsmltoso terres- of simulations dynamical of results present we Here h hrce n opsto fasse ftretilpla terrestrial of system a of composition and terrestrial character The that expect we arguments, above the on Based in lntfrsadmgae ucl,teplan- the quickly, migrates and forms planet giant t qa oo ihrta h at’.Frsmall For Earth’s. the than higher or to equal nts J & , ttepeec fa ue in lntsc as such planet giant outer an of presence the at ONATHAN .L I. exobiology – UNINE RSNEO HOT OF PRESENCE E 2 tal. et eSlrSse can System Solar he [email protected]) 2000). , tal. et v tfrplanets for it rve 04for 2004 , otoplan- tal. et star’s with con- her- ble he ts, n- e e s - , 2 which are discussed in section 4. Figure 2 shows the final state of twelve simulations (out of twenty), with the included for scale. Grey circles 2. METHOD represent the positions of giant planets in each simulation and are not on the same scale as the terrestrial bodies. The eccen- A simulation begins with a disk of which re- tricity of each body is shown beneath it by its radial excursion flects the minimum mass solar nebula model (Hayashi 1981). over the course of one orbit. Terrestrial planets can form in Planetary embryos have densities of 3 gcm−3 and are placed the habitable zone in the presence of a , often with from the hot jupiter out to 5.2 AU. These are randomly spaced substantial water contents. The possibility of a potentially hab- by 3-6 mutual Hill radii assuming the surface density of solids itable planet forming depends on the location of the hot jupiter. scales with heliocentric distance r as r−1.5. The surface density −2 In most cases, no planets more massive than 0.2 Earth masses is normalized to 10 gcm at 1 AU, with each disk of embryos form within a factor of 3 in period to the hot jupiter, roughly a containing 6-7 earth masses of material inside 5 AU. The dis- factor of two in semi-major axis. If a planet forms in the habit- covered giant planets are found to preferentially orbit stars with able zone with a hot jupiter at 0.5 AU, it is the innermostterres- metallicities higher than the Sun’s (Laws et al. , 2003), indi- trial planet and tends to be relatively small and dry. Water-rich cating that they likely contain a large amount of solid material planets form readily in the habitable zone with a hot jupiter at with which to build terrestrial planets. Our chosen value for the 0.15 or 0.25 AU. surface density is therefore quite low, and accounts for some In our Solar System planet formation was suppressed in the depletion during hot jupiter migration. All hot jupiters have by the gravitational effects of Jupiter. This is seen masses of 0.5 Jupiter masses and all outer giant planets are 1 in Fig. 2, as no terrestrial planets form within a factor of 3-4 . in period to a hot jupiter or an outer . We speculate We assign protoplanets an initial distribution of water con- that as this gap is filled with the remnants of terrestrial bodies tent which reflects the distribution in chondritic meteorites (see in our solar system, systems with hot jupiters may contain as- Fig. 2 from Raymond et al. 2004), such that the inner bodies teroid belts interior to the terrestrial planets. The resolution of are dry, past 2 AU planetary embryos contain 0.1% water, and current simulations is too low to test this hypothesis. past 2.5 AU embryos contain 5% water. Their iron distribution Simulations without a exterior to the terrestrial is interpolated between the content of the planets and chondritic region form planets of substantial mass in the asteroid belt and asteroid classes, ignoring the planet Mercury. These range from beyond on time scales of hundreds of Myr. Indeed, the systems 0.40 (40% iron by mass) at 0.2 AU to 0.15 at 5 AU. Each em- in Fig. 2 with no outer giant planet have not yet finished accret- bryo is given a small initial inclination ( 1◦) and eccentricity < ing. Simulations23 and 24 were runfor 200 Myr, sims 9 and 10 ( 0.02). < for 500 Myr, and sims 13 and 14 for 800+ Myr. A comparison Each simulation is evolved for at least 200 million years us- between the outer regions of these demonstrates the long for- ing a hybrid integrator called Mercury (Chambers 1999), which mation timescales. An outer gas giant clears the asteroid belt evolves the orbits of all bodies and keeps track of collisions. of protoplanets quickly, although the water content of terres- The hybrid scheme in Mercury uses a symplectic algorithm to trial planets is roughly the same in the absence of an outer giant evolve orbits of bodies unless they are involved in a close en- planet. In all cases, terrestrial planets in the habitable zone form counter, in which case it switches to a Bulirsch-Stoer method. morequickly in the presenceof an outergiant planet and are de- Collisions are treated as inelastic mergers which conserve mass livered water at earlier times than with no outer gas giant. This and water content. The time step in each simulation is chosen suggests that an outer giant planet’s net effect is to clear ma- to be less than 1/20 of the of the innermost body terial from the asteroid belt and to accelerate in the simulation, and ranges from 1 day for a hot jupiter at 0.15 formation. Its role in delivering water to the terrestrial planets AU to 6 days for a hot jupiter at 0.5 AU. Each simulation con- 5 is not a vital one in terms of quantity. served energyto at least one part in 10 , and took between three The amount of material ejected from the system is also a weeks and two months to complete on a desktop PC. function of the number and configuration of giant planets. An outer giant planet ejects approximately one half of the total ter- 3. RESULTS restrial mass in the system, while a hot jupiter can remove up Figure 1 shows the time evolution of one simulation which to one third. In the case of a hot jupiter at 0.5 AU and an outer formed a planet in the habitable zone, with a hot jupiter at 0.25 giant planet, the terrestrial planets comprise only one quarter AU and an outer giant planet at 5.2 AU (not shown). Planetary of the initial mass. These planets are systematically depleted embryos are dynamically excited by the giant planets and their in iron, because the inner, iron-rich material has been largely mutual gravitation, increasing their eccentricities and causing removed by the hot jupiter. their orbits to cross. This results in both accretional impacts We have run three simulations for one billion years or more and close encounters with giant planets, which eject roughly to test the long term stability of terrestrial planets in the pres- half of the terrestrial bodies. By the end of a simulation onlya ence of hot jupiters. The short dynamical timescales in the in- few terrestrial planets remain. In this case four terrestrial bodies ner disk result in a fast clearing of unstable objects, so a longer have formed including two planets inside 2 AU, one of which integration produces no change. The asteroid belt is slowly lies in the habitable zone at 1.06 AU with 1.68 times the mass cleared by an outer giant planet, but all planets which are well- of Earth with water content higher than the Earth’s. As our separated from a giant planet (by a factor of 3-4 or more in simulations do not account for water loss during impacts, wa- orbital period) are stable for long timescales. ter content values are upper limits. However, we do not sim- We have run two simulations under the assumption that the ulate the secondary delivery of volatiles from farther out in hot jupiter’s migration took place later in the lifetime of the the disk (“late veneer”) which would increase the water con- . The surface density of solid material was tent.although likely not by more than 10% if it proceeds as in reduced by a factor of five, and we included a hot jupiter at our Solar system (Morbidelli et al. , 2000). 0.25 AU and an outer giant planet. After 200 Myr of evolution, 3 these systems formed no planets more massive than 0.16 earth jupiter. Based on this and previous work it is possible to pre- masses and left a large number of small bodies in the terrestrial dict the character of the terrestrial planets around a star, from region, reminiscent of a large asteroid belt. observables such as the orbit and mass of a giant planet and the metallicity of the star. Our predictions will be testable in the 4. DISCUSSION near future with upcoming space missions such as Kepler5 and All the simulations presented here contain giant planets on COROT6, that will detect giant and (hopefully) terrestrial plan- circular orbits with fixed masses. The observed hot jupiters ets around other stars. Longer-term missions like Terrestrial inside 0.1 AU (51 Peg-type hot jupiters) tend to have circu- Planet Finder7 and Darwin8 hope to obtain spectra of terres- lar orbits due to tidal interactions with the central star. More trial planets and search for signs of water and life. We suggest distant giant planets can have a large range in eccentricity and that stars with hot jupiters may be a good place to look for extra mass. The effects of these parameters can be extrapolated using solar terrestrial planets. previous results. An eccentric giant planet preferentially ejects This result can also be applied to constrain the location of the water-rich material from the rather than scat- (Gonzalez, Brownlee & Ward 2002; tering it inward, which results in dry terrestrial planets (Cham- Lineweaver, Fenner & Gibson 2004). This is defined as the re- bers & Cassen 2002; Raymond et al. , 2004) with large eccen- gion in the galaxy in which various factors conspire to make the tricities, located far from the giant planet. A more massive giant area suitable for life (e.g. the average metallicity of stars, the planet or a higher surface density of solid material results in a rate of supernovae,time needed for life to evolve). In particular, smaller number of more massive terrestrial planets (Wetherill Lineweaver et al. , (2004) assume (from Lineweaver 2001) that 1996; Chambers & Cassen 2002; Raymond et al. , 2004). We the probability of a star to host a potentially habitable planet apply this to a known planetary system, 55 Cancri (Marcy et drops precipitously if its metallicity is higher than 0.2-0.3 dex al. , 2002), which contains two hot jupiters at 0.115 and 0.241 (solar metallicity is defined to be 0.0). This is based on the AU and an exterior giant planet at 5.9 AU. The hot jupiters are fact that higher metallicity stars are more likely to have hot close to being in 3:1 resonance and the less massive one has an jupiters (Laws et al. , 2003), and the assumption that any mi- eccentricity of 0.33. The outer giant planet’s eccentricityis 0.16 gration event would preclude the formation of terrestrial plan- and it is four times as massive as Jupiter. By our previous argu- ets in the system. Our result, that potentially habitable planets ments, we expect a small number of terrestrial planets to form can exist around stars with hot jupiters, effectively widens the in 55 Cancri far away from the hot jupiters as well as from the Galactic Habitable Zone to include regions at small galactocen- outer giant. The high eccentricities should strongly deplete the tric distances and recent times (“too metal rich” regions in Figs. solid material, resulting in low-mass planets. Simulations have 3 and 4 of Lineweaver et al. , 2004). shown this to be the case, with at most two terrestrial planets forming in 55 Cancri, with masses no greater than 0.6 earth 5. ACKNOWLEDGMENTS masses (Raymond & Barnes 2004). We have argued that terrestrial planets can form in the pres- We thank Lucio Mayer, Chris Laws, and Graeme Lufkin for ence of hot jupiters. We have shown that potentially habitable helpful discussions. This work was funded by NASA’s Astrobi- planets with orbits in the habitable zone and substantial water ology Institute and NASA Planetary . These sim- contents can form in such conditions. We hypothesize that as- ulations were run under Condor, which is publicly available at teroid belts may exist between the terrestrial planets and a hot http://www.cs.wisc.edu/condor.

REFERENCES

Alibert, Y., Mordasini, C., & Benz, W., 2004. Migration and giant planet Kuchner, M. J. & Lecar, M., 2002. Halting Planet Migration in the Evacuated formation. A&A, 417, L25. Centers of Protoplanetary Disks. ApJ, 574, L87. Armitage, P. J., 2003. A Reduced Efficiency of Terrestrial Planet Formation Laws, C., Gonzalez, G., Walker, K. M., Tyagi, S., Dodsworth, J., Snider, K., & following Giant Planet Migration. ApJ582, L47. Suntzeff, N. B. 2003, AJ, 125, 2664 Boss, A. P., 1997. Giant planet formation by gravitational instability. Science Levison, H. F. & Agnor, C., 2003. The Role of Giant Planets in Terrestrial 276, 1836. Planet Formation. AJ125, 2692. Briceño, C., Vivas, A. K., Calvet, N., Hartmann, L., Pacheco, R., Herrera, D., Lin, D. N. C., Bodenheimer, P., & Richardson, D. C., 1996. Orbital migration Romero, L., Berlind, P., Sanchez, G., Snyder, J. A., Andrews, P., 2001. The of the planetary companion of to its present location. , 380, CIDA-QUEST Large-Scale Survey of Orion OB1: Evidence for Rapid Disk 606. Dissipation in a Dispersed Stellar Population. Science 291, 93-97. Lineweaver, C. H., 2001. An Estimate of the Age Distribution of Terrestrial Chambers, J. E., 1999. A Hybrid Symplectic Integrator that Permits Close Planets in the Universe: Quantifying Metallicity as a Selection Effect. Icarus, Encounters between Massive Bodies. MNRAS, 304, 793-799. 151, 307. Chambers, J. E. & Cassen, P., 2002. The effects of nebula surface density Lineweaver, C. H., Fenner, Y., & Gibson, B. K., 2004. The Galactic Habitable profile and giant-planet eccentricities on planetary accretion in the inner solar Zone and the Age Distribution of Complex Life in the Milky Way. Science, system. Meteoritics and , 37, 1523. 303, 59. D’Angelo, G., Kley, W., & Henning, T., 2003. Orbital Migration and Mass Lufkin, G., Quinn, T., Wadsley, J., Stadel, J., & Governato, F., 2004. Accretion of Protoplanets in Three-dimensional Global Computations with Simulations of gaseous disc-embedded planet interaction. MNRAS, 347, Nested Grids. ApJ, 586, 540, and references therein 421. Forrest, W. J., et al. , 2004. Mid-Infrared Spectroscopy of Disks around Marcy, G. W., Butler, R. P., Fischer, D. A., Laughlin, G., Vogt, S. S., Henry, Classical T Tauri Stars. ApJS, in press. G. W., & Pourbaix, D., 2002. ApJ, 581, 1375. Gonzalez, G., Brownlee, D., & Ward, P., 2001. The Galactic Habitable Zone: Matsuyama, I., Johnstone, D., & Murray, N., 2003. Halting Planet Migration Galactic Chemical Evolution. Icarus, 152, 185. by Photoevaporation from the Central Source. ApJ, 585, L143. Hayashi, C. 1981. Prog. Theor. Phys. Suppl., 70, 35 Mayer, L., Quinn, T., Wadsley, J., & Stadel, J., 2002. Formation of Giant Planets Kasting, J. F., Whitmire, D. P., and Reynolds, R. T., 1993. Habitable zones by Fragmentation of Protoplanetary Disks. Science, 298, 1756. around main sequence stars. Icarus 101, 108-128.

5http://www.kepler.arc.nasa.gov 6http://www.astrsp-mrs.fr/projets/corot 7http://planetquest.jpl.nasa.gov/TPF 8http://ast.star.rl.ac.uk/darwin 4

Morbidelli, A., Chambers, J., Lunine, J. I., Petit, J. M., Robert, F., Valsecchi, Raymond, S. N., Quinn, T., & Lunine, J. I., 2004. Making other earths: G. B., and Cyr, K. E., 2000. Source regions and timescales for the delivery dynamical simulations of terrestrial planet formation and water delivery. of water on Earth. Meteoritics and Planetary Science 35, 1309-1320. Icarus, 168, 1. Nelson, R. P., Papaloizou, J. C. B., Masset, F., & Kley, W., 2000. The migration Rice, W. K. M., Armitage, P. J., 2003. On the Formation Timescale and Core and growth of protoplanets in protostellar discs. MNRAS, 318, 18. Masses of Gas Giant Planets. ApJ, 598, L54. Pollack, J. B., Hubickyj, O., Bodenheimer, P., Lissauer, J. J., Podolak, M., Rice, W. K. M., Armitage, P. J., Bonnell, I. A., Bate, M. R., Jeffers, S. V., & & Greenzweig, Y., 1996. Formation of the Giant Planets by Concurrent Vine, S. G., 2003. Substellar companions and isolated planetary-mass objects Accretion of Solids and Gas. Icarus, 124, 62. from protostellar disc fragmentation. MNRAS, 346, L36. Raymond, S. N. & Barnes, R., 2004. Predicting Planets in Known Extra-Solar Terquem, C. E. J. M. L. J., 2003. Stopping inward by a Planetary Systems III: Forming Terrestrial Planets. ApJ, submitted, astro- toroidal magnetic field. MNRAS, 341, 1157. ph/0404212 Wetherill, G. W., 1996. The Formation and Habitability of Extra-Solar Planets. Icarus, 119, 219

FIG. 1.— Six snapshots in time from a simulation with two giant planets (not shown): a 0.5 Jupiter-mass hot jupiter at 0.25 AU and a Jupiter-mass planet at 5.2 AU, both on circular, coplanar orbits. Each panel plots the eccentricity and semi-major axis of each surviving body in the simulation. The size of a body is proportional to its mass(1/3), and the dark region in the center represents the size of its iron core, on the same scale. The color corresponds to the water mass fraction, which ranges initially from 10−5 to 0.05. 5

FIG. 2.— Final configurations of twelve simulations, with the Solar System shown for scale. The gray circles represent the giant planets in each simulation and are not to the same scale as the terrestrial bodies. The eccentricity of each body is represented by its excursion in heliocentric distance over an orbit. The x axis is on a logarithmic scale such that a given separation corresponds to a fixed ratio of orbital periods, shown in the scale bar on the top left. The dashed vertical lines represent the boundaries of the habitable zone (Kasting et al. , 1993). Simulations 23 and 24 were run for 200 Myr, simulations 9 and 10 for 500 Myr, and simulations 13 and 14 for 800+ Myr. A comparison shows the long accretion timescales in the outer terrestrial region. Note the presence of protoplanets in 1:1 resonance with a giant planet in some cases.