sign in sign up
<<
Home , Planetary migration, Roche limit

arXiv:1704.01688v2 [astro-ph.EP] 7 Apr 2017 h eetydsoee RPIT1sse provides system TRAPPIST-1 discovered recently The ( Kepler-11 are measured systems these of 2013 stability things, orbital other ( amongst the allow, of masses exploration planetary the of 2005 plan- the system of compact precision the the of and ets period For- short velocities. the radial tunately precision from not measured could masses be that easily meant where also stars parameter-space host faint quently of region the Kepler in exist systems 2013 aitos(Ts ( Timing Transit (TTVs) through Variations masses the of measurement the aefo h al eut fthe of results early systems the planetary from compact of came windfall or- sys- initial eccentric planetary The compact in tems. and planets, These circumbinary massive system. bits, planetary , own hot our which in include for analog architectures no have orbital we of shortage no provided iso h rni ehdtwrssotobtlperiods orbital short towards ( method transit the of bias amn ta.2009 al. et Raymond 2008 Braun von & Kane [email protected] oal xmlso opc lntr ytm with systems planetary compact of examples Notable rpittpstuigL using 2017 typeset 10, Preprint April version Draft vrteps e eae,tefil feolnt has of field the decades, few past the Over OLSWTOTMOS XMO OSRIT O OPC PLAN COMPACT FOR CONSTRAINTS : WITHOUT WORLDS n elr8 ( Kepler-80 and ) ; ; ad ta.2014 al. et Hands eateto hsc srnm,SnFacsoSaeUniv State Francisco San , & Physics of Department omn&Mra 2005 Murray & Holman a h togs estvt.Hwvr h fre- the However, sensitivity. strongest the had n ftepiaysrrsso xpae eetoshsbe the been has within detections reside planets of numerous surprises whereby primary systems, the of One oto h lnt r eylkl ob olswtotmoons. without worlds be to likely TRAPPIST-1 very Keywords: the are syst to planets compact f constraints the in pool these of planets tidal most apply for We of th moons rate of considerations. of su the limit Zone presence a increase Habitable the of and the presence on variations the constraints within tilt is lie axial to stabilize the of can planets habitability h the the the of to which contributors more for systems or those one the are that enabling interest indicating particular years, Of recent in tecture. discovered been have systems 1. srbooy–paeaysses–sas niiul(TRAPPIST-1 individual stars: – systems planetary – INTRODUCTION A iad-su´ 2002 Miralda-Escud´e i 2013 Xie T ; E ; tl ATX .1.0 v. AASTeX6 style X uke l 2010 al. et Funk ekr&Aas2017 Adams & Becker en htteecompact these that meant ) Kepler ; .Temeasurements The ). aDnl ta.2016 al. et MacDonald isure l 2011 al. et Lissauer Kepler htmtyallowed photometry ; iso.The mission. ; hn tal. et Zhang gle al. et Agol ). tpe .Kane R. Stephen ABSTRACT ). , ∼ riy 60Hloa vne a rnic,C 43,USA 94132, CA Francisco, San Avenue, Holloway 1600 ersity, . Uo h otsa.Mn fteeknsof kinds these of Many star. host the of AU 0.5 rmtcciaecags td by study A in resulting changes. obliquity climate in dramatic variations extreme have would aiaiiytruholqiysaiiain al work Early Earth’s by in stabilization. role obliquity key through a habitability played substantial with eido h ue lntta lcdteobtna the near ( the line placed snow that TRAPPIST-1 outer the of period yicue h ffcso h ia xaso fthe of expansion by tidal simulations the Further of effects moon. the included by and ( 2015 2013 2013 Barnes & Heller Kipping & Forgan ( habitability 1987 potential their and moons the the by provided of photometry extension of baseline time longer three of dwarf uncovered ultracool (HZ) observations faint ( Zone the Initial star transiting Habitable planets the star. terrestrial occupy host plan- to the the of appear some where ets architecture orbital compact a h lnt ob eemnd( determined be of masses to the planets allowed transit- that the terrestrial TTVs including four planets, additional ing an revealed Spitzer at osntncsaiypeld habitability. preclude moonless necessarily a not thus does and Earth the amplitude at not determined though previously obliquity, Earth’s the stabilize indeed 1996 uhhsbe rte bu h eeto fexo- of detection the about written been has Much eew xlr h oeo ria rhtcueand architecture orbital of role the explore we Here akre al. et Lasker ie al. et Li a of presence the that held commonly also is It ). ; xaddtecatcolqiycalculations obliquity chaotic the expanded ) ilne l 2016 al. et Gillon ilase l 1997 al. et Williams s tri o-as hspotentially thus low-mass, is star ost r arycmo ria archi- orbital common fairly a are y m ae nHl ai n Roche and radii Hill on based ems sata onwoetdleffects tidal whose moon bstantial ( 2014 icvr fcmatplanetary compact of discovery rain eew xlr the explore we Here ormation. ( 1993 ytmaddmntaethat demonstrate and system Kepler eosrtdta h ondoes Moon the that demonstrated ) ; niae htamols Earth moonless a that indicated ) elre l 2014 al. et Heller .Cniudosrain using observations Continued ). otsa.Oeo the of One star. host e ) iso ofimdteorbital the confirmed mission ue ta.2017 al. et Luger ; ; ede on 2010 Horne & Weidner ikl&Kn 2013 Kane & Hinkel TR SYSTEMS ETARY ilne l 2017 al. et Gillon isure al. et Lissauer oael tal. et Tomasella ; enlse al. et Reynolds ipn tal. et Kipping ). ( 2012 .A ). K2 ) ; ; 2 Stephen R. Kane planetary masses and radii on the ability of the planets to harbor . We show that gravitational in- fluence of planets in compact planetary systems places severe constraints on the presence of exomoons and we use the TRAPPIST-1 system as an example. In Sec- tion 2, we quantify the limits imposed on exomoon or- bits by the Hill radius and Roche limits of a planet. In Section 3, we construct analytical limits for exomoon stability regions in the context of compact planetary sys- tem architectures. Section 4 applies these constraints to the TRAPPIST-1 system and demonstrates why many of the TRAPPIST-1 planets cannot harbor moons. We provide discussion and concluding remarks in Section 5.

2. HILL RADIUS AND ROCHE LIMIT Figure 1. The locations of equality between the Hill ra- The two main considerations for the region in which dius and Roche limit (see Equation 5) for five different mean a moon can be maintained in orbit around a planet are densities of a potential moon. Stable moon cannot exist above the curve for a given density. For the pur- the Hill radius, which defines the outer limit of the re- poses of these calculations, we adopt a reduction factor of gion, and the Roche limit, which defines the inner limit χ = 1/3. The horizontal dashed lines indicate the stellar (Kipping 2009; Heller 2012). The Hill radius, RH , is masses for two example compact planetary systems, Kepler- 80 and TRAPPIST-1, with the locations of the planets shown given by 1 on each line. 3 Mp RH = ap (1) 3M size of the allowed region for long-term moon orbital sta-  ⋆  bility is the evolution of the orbit through tidal interac- where ap is the planetary semi-major axis, Mp is the tions with the planet (Sasaki & Barnes 2014), effectively , and M⋆ is the stellar mass. The Roche resulting in an exchange of . The limit, RR, is given by inward or outward direction of the moon migration de- 1 3 ρp pends upon the ratio of the orbital period to the rotation RR ≃ 2.44Rp (2) ρm period of the planet (Barnes & O’Brien 2002) and also   if the orbit is prograde or retrograde (Domingos et al. where R is the planetary radius, ρ is the planetary p p 2006). Note that the semi-major axis of the Moon’s or- density, and ρ is the moon density. For example, the m bit is ∼26% of the Earth’s Hill radius, which is close Hill radius of the Earth is 0.01 AU and the Roche limit to the χ ≈ 1/3 criteria mentioned above. For compari- is 1.26 × 10−4 AU, a factor of ∼78.5 difference. For an son, the outermost of the , , has Earth-analog in a compact system with a semi-major a semi-major axis that is ∼3% of ’s Hill radius. axis of 0.05 AU, the Hill radius shrinks to 6 × 10−4 AU Clearly the effective size of the Hill for planets in which results in a factor of ∼4.5 difference with the compact systems can be greatly reduced by the proper- Roche limit. ties of the local environment. Given that the radius of the is a theoret- ical approximation, the practical outer limit of a stable 3. EXOMOON CONSTRAINTS FOR COMPACT orbit for a moon can be significantly smaller due other PLANETARY SYSTEMS perturbation effects, such as those originating from the Here we consider the Hill and Roche limit radii de- host star. As such, the outer limit may be described as scribed in Section 2 in the extreme environments of com- χRH where χ is a reduction factor due to the above de- pact planetary systems. Within the orbital regime of scribed effects (Holman & Wiegert 1999). Further per- these systems, the gravitational perturbations are dom- turbation sources are other bodies in the system, partic- inated by the host star rather than the planets, en- ularly planetary bodies (Gong et al. 2013; Payne et al. abling stability through the shrinking of the planetary 2013). Furthermore, formation processes play a role Hill (Namouni 2010). An upper limit of the in the destabilization of moons, such as their poten- presence of moons may be determined by evaluated con- tial for resonant removal during planetary migration ditions where χR = R . Using Equations 1 and 2, we (Spalding et al. 2016). For example, Barnes & O’Brien H R find the following relationship: (2002) and Kipping (2009) estimate χ ≈ 1/3, effec- 1 tively reducing the outer limit of a stable moon orbit 1 3 M 3 1 M × p p to 1/3 RH . An additional process that will reduce the χap =2.44Rp 4 3 (3) 3M⋆ ρm πR   3 p ! Worlds Without Moons 3 which is simplified as exist within regions of parameter space that lie above the curve. For example, no planet in the Kepler-80 system M⋆ ρm = 10.4 3 (4) within 0.035 AU of the host star can ever host a moon (χap) with a density ≤ 3 g/cm3. The cubic nature of the re- A more convenient form of the relationship is as follows: lationship described in Equation 5 does ensure that the −3 concise exclusion of moons predominantly affects those ρm −6 −3 M⋆ ap 3 =6.18 × 10 χ (5) planets in the very small semi-major axis regime. An g/cm ! M⊙ 1 AU     interesting possibility is that the moon constraints for Equation 5 is plotted in Figure 1 for densities ranging compact systems may still allow the presence of ring from 1 g/cm3 to 5 g/cm3, where we have assumed a systems (Barnes & Fortney 2004; Zuluaga et al. 2015), Hill radius reduction factor of χ = 1/3. Horizontal particularly considering the relatively long persistence dashed lines indicate the stellar masses for the Kepler- of debris disks around low-mass stars (Plavchan et al. 80 (0.73 M⊙) and TRAPPIST-1 (0.0802 M⊙) host stars, 2009). Such rings around terrestrial planets are likely the planets of which largely fall within the range of semi- to have low optical depth and not readily detectable major axes shown (MacDonald et al. 2016; Gillon et al. from transit photometry when they are coplanar with 2017). Since the solid curves represent the boundary the inclination of the planetary orbits. where χRH = RR, moons of the given density cannot

Table 1. TRAPPIST-1 planetary Hill radii and Roche limits

−3 −3 Planet Mp (M⊕) Rp (R⊕) a (AU) RH (10 AU) RR (10 AU) RH /RR TRAPPIST-1 b 0.85 1.086 0.011 0.244 0.120 2.04 TRAPPIST-1 c 1.38 1.056 0.015 0.393 0.141 2.79 TRAPPIST-1 d 0.41 0.772 0.021 0.370 0.094 3.94 TRAPPIST-1 e 0.62 0.918 0.028 0.557 0.108 5.17 TRAPPIST-1 f 0.68 1.045 0.037 0.756 0.111 6.80 TRAPPIST-1 g 1.34 1.127 0.045 1.154 0.139 8.28 TRAPPIST-1 h 0.31 0.715 0.060 0.936 0.086 10.86

Note—Planetary masses, radii, and semi-major axes are from Gillon et al. (2017).

4. NO MOONS IN THE TRAPPIST-1 SYSTEM? of χ = 1, the difference between the Hill radii and Roche The TRAPPIST-1 system is an exceptional case of limits for the inner planets is remarkably small, similar compact planetary systems due to the extremely low to the case for the Earth-analog at 0.05 AU described in mass of the host star. The combination of transits and Section 2. TTVs provide radii and masses for the planets, allow- Shown in Figure 2 are moon exclusion boundaries ing us to calculate the Hill radii and Roche limits for (solid lines) for the TRAPPIST-1 system as a function each planet (see Section 2). Currently the outer planet of mean density and semi-major axis, with vertical dot- (planet h) does not have a reliable mass estimate, and ted lines indicating the location of the six innermost so we used the mean density of the other six planets planets. The shaded regions show the extent of the HZ for TRAPPIST-1, where we used the TRAPPIST- of ∼0.84 ρ⊕ to determine a mass estimate of 0.31 M⊕. The resulting calculations for all seven planets are shown 1 luminosity (0.000524 L⊙) and effective temperature in Table 1. Also shown are the ratios of the Hill radii (2559 K) provided by Gillon et al. (2017), along with and Roche limits, where a Hill radius reduction factor the HZ boundary equations of Kopparapu et al. (2013, of unity has been used and a moon density of 3 g/cm3 2014). The shaded regions correspond to the “conser- adopted for the Roche limit. For comparison, the mean vative” (light-gray, 0.024–0.049 AU) and “optimistic” density of the Galilean moons is 2.6 g/cm3 and the mean (dark-gray, 0.019–0.051 AU) HZ boundaries, the defi- density of the Earth’s moon is 3.3 g/cm3. The ratios in nitions of which depends upon assumptions regarding the final column of Table 1 show that, even for the case the longevity of liquid water on the surfaces of and . Although the full extent of the perturbations 4 Stephen R. Kane

(Kane et al. 2013). An aspect of particular interest is the presence of compact systems around low-mass stars where the chances of temperate terrestrial planets are greatly enhanced. The TRAPPIST-1 system is an ideal example of such a system, with three planets within the conservative HZ and four within the optimistic HZ. This has resulted in the postulated habitability of the TRAPPIST-1 planets (Gillon et al. 2017), including de- tailed models of the planetary atmospheres (Wolf 2017). The presence of moons can have a positive impact on habitability, such as promoting tidal pools and stabi- lizing rotational obliquity. Insofar as habitability is re- duced by a lack of moons, compact planetary systems such as TRAPPIST-1 may suffer due to the constraints Figure 2. The allowed moon density as a function of semi- on the presence of moons outlined in this work. As major axis for the TRAPPIST-1 system and for Hill radii re- shown, a dominant deciding factor is the mean density of duction factors of χ = 1/3 and χ = 1/4. Stable moon orbits the moon which depends on whether they formed in situ cannot exist below each of the curves. The vertical dotted lines indicate the location of the six innermost TRAPPIST-1 around the planet and the location of planet formation. planets. The light-gray shaded region shows the extent of the Planetary migration of planets near or beyond the snow conservative HZ, and the dark-gray shaded region represents line may result in a lower than expected moon density the optimistic extension to the HZ. which will in turn ensure their removal according to the calculations described herein. to potential moons is not quantified, the existence of The presence of moons, or lack thereof, in the sys- resonant chains of planetary orbits (Gillon et al. 2017) tem is a testable hypothesis given the techniques avail- is expected to increase the amplitude of such pertur- able from transit variations. A moon will produce bations. Therefore, we have used Hill radii reduction TTV and/or Transit Duration Variation (TDV) effects factors of χ = 1/3 and χ = 1/4 in the figure. In this in the timing and shape of the transit photometry, case, the exclusion boundaries imply that the parame- depending on the mass and separation of the moon ter space below the curve show where stable orbits of (Sartoretti & Schneider 1999; Kipping 2009). Addition- moons are not possible. For example, the χ =1/4 curve ally, the size of the moon will play a major role in its indicates that planets b–e cannot host moons of reason- detectability in the photometric data. For example, the able (. 3 g/cm3) mean densities. As noted at the end of radius ratio of the Earth’s moon and the Earth is 0.272; Section 2, the Hill radii reduction factors do not include thus the Earth’s moon has a transit depth that is 7.4% the effects of moon orbital migration due to a tidal ex- as large as the Earth’s. Such a deviation would have change of angular momentum (Sasaki et al. 2012). Thus been detected in the precision Spitzer photometry pre- the region for which each planet can sustain a moon sented by Gillon et al. (2017), and so moons of compara- in a stable orbit over long time periods is further re- ble size ratios are likely not present, consistent with the duced such that it is highly unlikely that any of the findings of this work. Moons as small as the minimum planets represented in Figure 2 host moons. Properties size to be round (radius of 200–300 kms) are unlikely of the host star, such as activity and rotation rate, sug- to be detected without exceptional circunstances, such gest that the may be relatively young as phase folding planetary transits for a low-mass star (Gillon et al. 2016). The outer planets may therefore with relatively low stellar activity (Heller 2014). Long- currently hold moons in temporarily stable locations al- term monitoring with current and future facilities, such though, as noted earlier, the process of migration into as thirty-meter class ground-based telescopes and the regions of resonance trapping may have stripped them James Webb Space Telescope, will place further con- already (Spalding et al. 2016). straints on the presence of moons in these system.

5. CONCLUSIONS ACKNOWLEDGEMENTS Although the discovery of compact planetary systems The author would like to thank Duncan Forgan, Ren´e was initially a surprising outcome of exoplanet obser- Heller, and Natalie Hinkel for their insightful feedback vations, it now seems as though they may be a quite on this work. This research has made use of the NASA common occurrence. In many cases, their architec- Exoplanet Archive, which is operated by the Califor- tures have been compared to that of the Jovian system nia Institute of Technology, under contract with the Na- leading to speculation of similar formation mechanisms tional Aeronautics and Space Administration under the Worlds Without Moons 5

Exoplanet Exploration Program. The results reported Science (NExSS) research coordination network spon- herein benefited from collaborations and/or information sored by NASA’s Science Mission Directorate. exchange within NASA’s Nexus for Exoplanet System

REFERENCES

Agol, E., Steffen, J., Sari, R., Clarkson, W. 2005, MNRAS, 359, Lasker, J., Joutel, F., Robutel, P. 1993, Nature, 361, 615 567 Li, G., Batygin, K. 2014, ApJ, 790, 69 Barnes, J.W., O’Brien, D.P. 2002, ApJ, 575, 1087 Lissauer, J.J., Fabrycky, D.C., Ford, E.B., et al. 2011, Nature, Barnes, J.W., Fortney, J.J. 2004, ApJ, 616, 1193 470, 53 Becker, J.C., Adams, F.C. 2017, MNRAS, in press Lissauer, J.J., Barnes, J.W., Chambers, J.E. 2012, Icarus, 217, 77 (arXiv:1702.07714) Lissauer, J.J., Jontof-Hutter, D., Rowe, J.F., et al. 2013, ApJ, Domingos, R.C., Winter, O.C., Yokoyama, T. 2006, MNRAS, 770, 131 373, 1227 Luger, R., Sestovic, M., Kruse, E., et al. 2017, NatAs, submitted Forgan, D., Kipping, D. 2013, MNRAS, 432, 2994 (arXiv:1703.04166) Funk, B., Wuchterl, G., Schwarz, R., Pilat-Lohinger, E., Eggl, S. MacDonald, M.G., Ragozzine, D., Fabrycky, D.C., et al. 2016, 2010, A&A, 516, 82 AJ, 152, 105 Gillon, M., Jehin, E., Lederer, S.M., et al. 2016, Nature, 533, 221 Miralda-Escud´e, J. 2002, ApJ, 564, 1019 Gillon, M., Triaud, A.H.M.J., Demory, B.-O., et al. 2017, Namouni, F. 2010, ApJ, 719, L145 Nature, 542, 456 Payne, M.J., Deck, K.M., Holman, M.J., Perets, H.B. 2013, ApJ, Gong, Y.-X., Zhou, J.-L., Xie, J.-W., Wu, X.-M. 2013, ApJ, 769, 775, L44 L14 Plavchan, P., Werner, M.W., Chen., C.H., et al. 2009, ApJ, 698, Hands, T.O., Alexander, R.D., Dehnen, W. 2014, MNRAS, 445, 1068 749 Raymond, S.N., Barnes, R., Veras, D., et al. 2009, ApJ, 696, L98 Heller, R. 2012, A&A, 545, L8 Reynolds, R.T., McKay, C.P., Kasting, J.F. 1987, AdSpR, 7, 125 Heller, R., Barnes, R. 2013, AsBio, 13, 18 Sasaki, T., Barnes, J.W., O’Brien, D.P. 2012, ApJ, 754, 51 Heller, R. 2014, ApJ, 787, 14 Sartoretti, P., Schneider, J. 1999, A&AS, 134, 553 Heller, R., Williams, D., Kipping, D., et al. 2014, AsBio, 14, 798 Sasaki, T., Barnes, J.W. 2014, IJAsB, 13, 324 Hinkel, N.R., Kane., S.R. 2013, ApJ, 774, 27 Spalding, C., Batygin, K., Adams, F.C. 2016, ApJ, 817, 18 Holman, M.J., Wiegert, P.A. 1999, AJ, 117, 621 Tomasella, L, Marzari, F., Vanzani, V. 1996, P&SS, 44, 427 Holman, M.J., Murray, N.W. 2005, Science, 307, 1288 Weidner, C., Horne, K. 2010, A&A, 521, A76 Kane, S.R., von Braun, K. 2008, ApJ, 689, 492 Williams, D.M., Kasting, J.F., Wade, R.A. 1997, Nature, 385, Kane, S.R., Hinkel, N.R., Raymond, S.N. 2013, AJ, 146, 122 234 Kipping, D.M. 2009, MNRAS, 392, 181 Wolf, E.T. 2017, ApJL, submitted (arXiv:1703.05815) Kipping, D.M., Schmitt, A.R., Huang, X., et al. 2015, ApJ, 813, Xie, J.-W. 2013, ApJS, 208, 22 14 Zhang, K., Hamilton, D.P., Matsumura, S. 2013, ApJ, 778, 6 Kopparapu, R.K., Ramirez, R., Kasting, J.F., et al. 2013, ApJ, Zuluaga, J.I., Kipping, D.M., Sucerquia, M., Alvarado, J.A. 765, 131 2015, ApJ, 803, L14 Kopparapu, R.K., Ramirez, R.M., SchottelKotte, J., et al. 2014, ApJ, 787, L29