A comprehensive study of Kepler phase curves and secondary eclipses – temperatures and albedos of confirmed Kepler giant planets Daniel Angerhausen Rensselaer Polytechnic Institute, 110 Eighth Street, Troy, NY 12180; ORAU NASA Postdoctoral Fellow, NASA Goddard Space Flight Center, Greenbelt, MD 20771 Em DeLarme Rensselaer Polytechnic Institute, 110 Eighth Street, Troy, NY 12180; UCF, Planetary Sciences Group, Department of Physics, Orlando, FL 32816 and Jon A. Morse Rensselaer Polytechnic Institute, 110 Eighth Street, Troy, NY 12180; BoldlyGo Institute, 1370 Broadway, 5th Floor Suite 572, New York, NY 10018 Received ; accepted submitted to PASP arXiv:1404.4348v3 [astro-ph.SR] 24 Aug 2015 –2– ABSTRACT We present a comprehensive study of phase curves and secondary eclipses in the Kepler data set using all data from 16 quarters that were available in 2013- 2014. Our sample consists of 20 confirmed planets with Rp > 4Re , P < 10d, Vmag < 15. Here we derive their temperatures and albedos, with an eye towards constraining models for the formation and evolution of such planets. Where there was overlap our results confirm parameters derived by previous studies, whereas we present new results for Kepler 1b-8b, 12b-15b, 17b, 40b, 41b, 43b, 44b, 76b, 77b, and 412b derived in a consistent manner. We also present lightcurve analyses for Kepler 91b and Kepler 74b, which both show extra dimmings at times other than from the expected primary and secondary eclipses. Corrected for thermal emission we find most of the massive planets from our sample to be low in albedo (< 0.1) with a few having higher albedo (> 0.1). Subject headings: planets and satellites: atmospheres planets and satellites: fundamental parameters planets and satellites: gaseous planets planets and satellites: general –3– 1. Introduction 1.1. The Kepler mission Studying extrasolar planets is one of the major frontiers of astronomy today. The field has transformed from simple identification to comprehensive categorization and characterization of exoplanets and exoplanetary systems. Analyses of data provided by NASA’s Kepler1 mission has revolutionized this field by compiling a statistically significant number of transiting planets and planetary candidates (e.g., Borucki et al. 2010b; Borucki et al. 2011; Batalha et al. 2013; Burke et al. 2014; Rowe et al. 2015; Mullally et al. 2015). For example, Kepler data allowed researchers to discover Kepler 9b (Holman et al. 2010), the first multi-planetary system outside our solar system; Kepler 10b (Batalha et al. 2011), one of the first confirmed rocky planets outside the solar system; and Kepler 16b (Doyle et al. 2011), the first circumbinary planet. More recently the Kepler team announced the discovery of potentially habitable worlds in the Kepler 62 (Borucki et al. 2013) and Kepler 69 (Barclay et al. 2013) systems, and the near-Earth-sized planets Kepler-186f (Quintana et al. 2014) and Kepler-452b (Jenkins et al. 2015) in the habitable zones of their parent stars. 1The data presented in this paper were obtained from the Mikulski Archive for Space Telescopes (MAST), which is managed by the Space Telescope Science Institute (STScI). STScI is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. Support for MAST for non-HST data is provided by the NASA Science Mission Directorate via grant NNX13AC07G and by other grants and contracts. This paper includes data collected by the Kepler mission. Funding for the Kepler mission is also provided by the NASA Science Mission Directorate. –4– Deeper analyses are possible using the exquisite Kepler data beyond merely detecting exoplanetary systems: researchers are now able to analyze large samples of planetary candidates to pin down occurrence rates such as ηearth (e.g., Howard et al. 2012; Dressing & Charbonneau 2013; Fressin et al. 2013; Petigura et al. 2013; Burke et al. 2015), find non-transiting planets via transit timing variations (Ballard et al. 2011), perform phase-curve analyses (Faigler et al. 2013), and may eventually even be able to detect exomoons (Kipping et al. 2013) and exotrojans (Hippke & Angerhausen 2015). For the close-in, and therefore hot, planets around bright, high-signal host stars in the Kepler data set, we are able to analyze secondary eclipses, i.e., the modulated flux from the star-planet system when the light (reflection and thermal emission) of the planet disappears during its passage behind the parent star. Differential measurements then help us to characterize physical parameters of the planet such as albedo and temperature. Measuring such quantities, together with complementary spectroscopic measurements, can provide contraints on models for the formation and atmospheric photochemistry of such close-in planets (e.g., Line et al. 2010; Moses et al. 2011; Visscher & Moses 2011) 1.2. Transits and eclipses Systems with transiting extrasolar planets can offer two important observational opportunities for deriving physical parameters of the planets. In primary transit the planet crosses the star. From a broadband transit-lightcurve, in this case, one can measure the 2 planetary radius Rp in units of the stellar radius R∗. The depth of the transit is ∼ (Rp/R∗) , which for a Jupiter radius planet transiting a sun-like star, is of the order of ∼ 1% (e.g., Henry et al. 2000). If the geometry (inclination, eccentricity) is right the planet also disappears behind its –5– Table 1. Kepler Quarters used in our analysis KOI LC Quarters SC Quarters 1 0,1,2,3,4,5,6,7,9,10,11,13,14,15 0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 2 0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 3 0,1,2,3,4,5,6,8,9,10,12,13,14 0,1,2,3,4,5,6,8,9,10,11,12,13,14 7 0,1,2,3,4,5,6,7,9,10,11,13,15 0,1,2,3,4,5,6,7,9,10,11,13,14,15 10 0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 1,2,3,4,5,6,7,8,9,10,11,12,13 13 0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 1,2,3,6,7,8,9,10,11,12,13,14,15 17 0,1,2,3,4,5,6,8,9,10,12,13,14 1,2,3,4,5,6,8,9,10,11,12 18 0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 1,2,3,4,5,6,7,8,9,10,11,12 20 0,1,2,3,4,5,6,7,9,10,11,13,14,15 1,2,3,4,5,6,7,8,9,10,11 97 0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 2,3,4,5,6,7,8 98 0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 2,3,4,5,6,7,8,9,10,11,12 127 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 2,3,4,5,6,7 128 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 2,3,4,5,6,7 135 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 2,3,4,5,6,7,9,10,11,12,13,14,15 196 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 3,4,5,6,7 200 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 3,4,5,6,7 202 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 3,4,5,6,7 203 1,2,3,4,5,6,8,9,10,12,13,14 3,4,5,6,7,8,9,10,11,12,13,14 204 1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 3,4,5,6,7 428 1,2,3,4,5,6,8,9,10,12,13,14 – 1658 1,2,3,4,5,7,8,9,11,12,13,15 – 2133 0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15 – Note. — For our analysis we used all data available prior to data release Q16, August 2013 –6– host star in a so-called secondary eclipse. For a 2000 K hot Jupiter-size planet, the typical flux deficit during secondary eclipse is ∼ 200 ppm at ∼ 2µm in the near-infrared and even larger at longer wavelengths, but considerably smaller at optical wavelengths (400-900 nm) at which Kepler observes. However, for host stars that are bright enough, Kepler’s outstanding sensitivity provides a direct measure of the planet’s disk-averaged day side flux for some of its targets, particularly close-in gas giants – so called ’Hot Jupiters’ and ’Hot Neptunes’. Observing secondary eclipses combined with planetary phase curves can help us to characterize the planet and its atmosphere. For example, the depth of the secondary eclipse can constrain the albedo of the planet, while the timing and width of the secondary eclipse can help determine its orbital parameters. Comparing the amplitude of the reflected light in the phase curve with the depth of the secondary eclipse can constrain the day and night side temperatures, and therefore confirm the planetary nature of a candidate that is not self-luminous, and help to understand day to night side heat exchange.