Arxiv:1404.4368V2 [Astro-Ph.EP] 5 Sep 2014 Univ
Total Page:16
File Type:pdf, Size:1020Kb
Draft version September 8, 2014 Preprint typeset using LATEX style emulateapj v. 5/2/11 FORMATION, TIDAL EVOLUTION AND HABITABILITY OF THE KEPLER-186 SYSTEM Emeline Bolmont1,2,*, Sean N. Raymond1,2, Philip von Paris3,+, Franck Selsis1,2, Franck Hersant1,2, Elisa V. Quintana4,5 and Thomas Barclay5,6 Draft version September 8, 2014 ABSTRACT The Kepler-186 system consists of five planets orbiting an early M dwarf. The planets have physical radii of 1.0{1.50 R⊕ and orbital periods of 4{130 days. The 1:1R⊕ Kepler-186f with a period of 130 days is of particular interest. Its insolation of roughly 0:32 S⊕ places it within the surface liquid water habitable zone (HZ). We present a multifaceted study of the Kepler-186 system, using two sets of parameters which are consistent with the data and also self-consistent. First, we show that the distribution of planet masses can be roughly reproduced if the planets were accreted from a high surface density disk presumably sculpted by an earlier phase of migration. However, our simulations predict the existence of one to two undetected planets between planets e and f. Next, we present a dynamical analysis of the system including the effect of tides. The timescale for tidal evolution is short enough that the four inner planets must have small obliquities and near-synchronous rotation rates. The tidal evolution of Kepler-186f is slow enough that its current spin state depends on a combination of its initial spin state, its dissipation rate, and the stellar age. Finally, we study the habitability of Kepler-186f with a one-dimensional climate model. The planet's surface temperature can be raised above 273 K with 0.5{5 bars of CO2, depending on the amount of N2 present. Kepler-186f represents a case study of an Earth-sized planet in the cooler regions of the HZ of a cool star. Subject headings: methods: numerical { planets and satellites: atmospheres { planets and satellites: dynamical evolution and stability { planets and satellites: formation { stars: indi- vidual (Kepler-186, KIC 8120608) 1. INTRODUCTION sustain surface liquid water (e.g., Wordsworth et al. The Kepler mission (Borucki et al. 2010) has made 2010). Kepler-186f receives a comparable or perhaps key discoveries on the road to finding Earth-like planets slightly higher stellar flux than GJ 581, placing it more (e.g., Batalha et al. 2011; Borucki et al. 2012; Fressin comfortably within the HZ. et al. 2012; Borucki et al. 2013). The recent detection of Here, we use the definition of the classical HZ (HZ, an Earth-sized planet in the HZ of an M star (i.e., the e.g., Dole 1964; Hart 1978; Kasting et al. 1993; Selsis Kepler-186 system, Quintana et al. 2014) brings us a step et al. 2007a; Kopparapu et al. 2013). Acknowledging the closer to finding a true Earth twin. fact that all terrestrial life needs liquid water, the HZ is The Kepler-186 planetary system hosts five known defined as the region around a star where a terrestrial planets including Kepler-186f, an Earth-sized planet in planet can host liquid water on the surface. The extent the HZ (Selsis et al. 2007a; Kopparapu et al. 2013). Fig- of this HZ naturally depends on the atmospheric con- ure 1 shows a comparison between the Kepler-186 sys- ditions (composition, pressure) as well as on the prop- tem, the solar system, and two other systems with poten- erties of the central star. Many more factors influence tially habitable planets: Kepler-62 (Borucki et al. 2013) the width of the HZ, such as the geological activity (e.g., and GJ 581 (Udry et al. 2007; Mayor et al. 2009). Cli- Lammer et al. 2010), the biosphere itself (e.g., Grenfell mate models have shown that GJ 581d, a super-Earth et al. 2010), or the dynamical environment of the plan- near the outer edge of the HZ of its host M star, could etary system (e.g., Menou & Tabachnik 2003; Barnes & Raymond 2004; Jones et al. 2006; S´andoret al. 2007; 1 arXiv:1404.4368v2 [astro-ph.EP] 5 Sep 2014 Univ. Bordeaux, Laboratoire d'Astrophysique de Bordeaux, Kopparapu & Barnes 2010). UMR 5804, F-33270 Floirac, France 2 CNRS, Laboratoire d'Astrophysique de Bordeaux, UMR We present a three-pronged study of the Kepler-186 5804, F-33270 Floirac, France system. We first try to reproduce the orbital architecture 3 Institut f¨ur Planetenforschung, Deutsches Zentrum f¨ur of the system using simple accretion simulations (Section Luft- und Raumfahrt (DLR), Rutherfordstr. 2, D-12489 Berlin, Germany 3). We show that certain features of the system{such as 4 SETI Institute, 189 Bernardo Ave, Suite 100, Mountain the large gap between planets e and f{are hard to ex- View, CA 94043, USA plain. We next briefly discuss the long-term dynamical 5 NASA Ames Research Center, Moffett Field, CA 94035, stability of the system (Section 4). In Section 5, we study USA 6 Bay Area Environmental Research Institute/NASA Ames the long-term dynamical, tidal, and spin evolution of the Research Center, Moffett Field, CA 94035, USA system. We use both simple tidal models and N-body + Current address: Univ. Bordeaux, Laboratoire simulations which include both tides and general relativ- d'Astrophysique de Bordeaux, UMR 5804, F-33270 Floirac, France and CNRS, Laboratoire d'Astrophysique de Bordeaux, ity. Next, we study the atmospheric conditions needed UMR 5804, F-33270 Floirac, France to bring Kepler-186f's surface temperature into the liq- * Author to whom all correspondence should be addressed. uid water range (Section 6). We discuss our findings and email: [email protected] conclude in Section 7. 2 Bolmont et al. Astronomical Units (AU) forming MCMC realizations. The stellar properties -0.4 -0.2 0.0 0.2 0.4 of this set correspond to those of table S1 in Quin- Orbits of tana et al. (2014). Kepler-186’s planets Set A and B are both valid sets of parameters, meaning that they are consistent with the data and are also self- consistent. c Line of sight The stellar mass is 0:5359 M in set A and 0:478 M to Earth in set B (Table 1 and Quintana et al. 2014). This pro- b f duces different values for the planets' semi-major axes d and insolations (Table 2). The same is true for their e planetary radii, as these are derived from the detected transit depth and the stellar radius. The fact that there is such a difference between the parameters values of each Habitable Zone set illustrate the limited knowledge we have of the sys- Potentially habitable tem. We expect the real parameters of the system to planets be within the range of sets A and B, so we performed Planetary system Too hot for life dynamical simulations for both these sets. However, av- comparison eraging between set A and B is not an option in so far Solar System as it would not correspond to anything realistic, the self- consistency would be lost (as in Quintana et al. 2014). Kepler-62 That is why we chose a different approach than in Quin- tana et al. (2014) where semi-major axes are not chosen Kepler-186 to match the derived stellar properties, but are rather determined by the transit model only. All five planets in the Kepler-186 system have radii GJ 581 between 1.0 and 1.5 R⊕ for both parameter sets. Given that low-density, gas-dominated planets tend to be larger than 1:5{2R⊕ (Weiss et al. 2013a; Weiss & Marcy 2014; 100 10 5 2 1 0.5 0.25 Lopez & Fortney 2013; Jontof-Hutter et al. 2014; Marcy Incident flux (normalized to Earth) et al. 2014b,a), all five Kepler-186 planets are probably Fig. 1.| Orbital configuration of the Kepler-186 system. The top rocky or at least solid. The planets' masses have not part shows a top-down view of the system, assuming orbits from been constrained with radial velocity or transit timing set A. The habitable zone boundaries are from Kopparapu et al. measurements (Quintana et al. 2014). Table 3 shows the (2013): the inner boundaries are the moist/runaway greenhouse limits and the outer boundaries are the maximum greenhouse and range of plausible planetary masses assuming a range early Mars limits. The sizes of the symbols are not to scale with the of compositions: 100% ice, 50% ice/ 50% rock, Earth- planetary orbits. The bottom part of the plot shows a comparison like composition, and 100% iron (following Fortney et al. between four different planetary systems which contain planets in 2007). the HZ: the solar system, Kepler-62 (Borucki et al. 2013), Kepler- 186 (Quintana et al. 2014), and GJ 581 (Udry et al. 2007; Mayor 3. FORMATION et al. 2009). Note that the inner moist and runaway greenhouse limits of the habitable zone are the same for Kepler-62, Kepler- 3.1. Formation Models 186, and GJ 581. Given the consistent insolation scaling, the x axis is linear in orbital distance but the scale is different for each At least six candidate mechanisms have been proposed system. The planets' relative sizes are correct, although for GJ 581 to explain the origin of close-in low-mass planets (Ray- the planetary radii were calculated as R = [Msin(i)]2:06, following mond et al. 2008a, 2013). In theory, these mechanisms Lissauer et al. (2011). can be distinguished using two observable quantities: 2. MODEL INPUT PARAMETER the inner planetary system architecture and the mean The stellar properties and planetary parameter given planet density.