1538-3881/150/2/56

We report on the discovery and validation of Kepler-452b, a transiting planet identified by a search through the 4 years of data collected by NASA's Kepler Mission. This possibly rocky planet orbits its G2 host star every days, the longest orbital period for a small ( ) transiting exoplanet to date. The likelihood that this planet has a rocky composition lies between 49% and 62%. The star has an effective temperature of 5757 ± 85 K and a of 4.32 ± 0.09. At a mean orbital separation of AU, this small planet is well within the optimistic habitable zone of its star (recent Venus/early Mars), experiencing only 10% more flux than Earth receives from the Sun today, and slightly outside the conservative habitable zone (runaway greenhouse/maximum greenhouse). The star is slightly larger and older than the Sun, with a present radius of and an estimated age of ~6 Gyr. Thus, Kepler-452b has likely always been in the habitable zone and should remain there for another ~3 Gyr.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

The study of exoplanets began in earnest just 20 years ago with the discovery of 51 Peg b, a hot jupiter in a 4.2 day period orbit of its host star (Mayor & Queloz 1995). This discovery was followed by a trickle of hot jupiters that became a torrent of giant planets over a large range of orbital periods, chiefly by means of radial velocity observations. Transit surveys were already underway as early as 199413 (see, e.g., Doyle et al. 1996), and rapidly gained momentum once the hot jupiter population was established. These surveys did not bear fruit until 2000 when the Doppler-detected planet HD 209458b (Mazeh et al. 2000) was seen in transit (Charbonneau et al. 2000). The field has advanced at a rapid pace ever since, enabling discoveries at longer orbital periods and smaller sizes. Today ~1500 exoplanets have been found through a variety of means, including Doppler surveys, ground- and space-based transit surveys, gravitational microlensing surveys (starting in 1996 with OGLE-2003-BLG-235/MOA-2003-BLG-53, Bennett et al. 2006), and recently, direct imaging (beginning in 2011 with HR 8799, Soummer et al. 2011). As impressive as the progress has been in this nascent field of astronomy, perhaps the most exalted goal, and one that remains tantalizingly close, is the discovery of a sufficient number of (near) Earth–Sun analogs to inform estimates of the intrinsic frequency of Earth-size planets in the habitable zone of Sun-like stars. The habitable zone is an evolving concept, here taken to be that range of distances about a star permitting liquid water to pool on the surface of a small rocky planet. For our purposes, we adopt the optimistic habitable zone as per Kopparapu et al. (2013) with the insolation flux, , between that experienced by recent Venus as the inner edge, and that experienced by early Mars as the outer edge. While the exact values depend somewhat on the effective temperature of the star, the optimistic habitable zone lies within the range from ~20% to ~180% of the radiation experienced by Earth today, that is from 0.2 to 1.8 . To date, nine small planets with radii less than 2.0 have been found in the 4 years of Kepler data in or near the habitable zones of their stars. These discoveries include Kepler-62e and f (Borucki et al. 2013), Kepler-186f (Quintana et al. 2014), a few among the multitude of planets reported by Rowe et al. (2014): Kepler-283c, Kepler-296e and f; and a set of small planets announced in Torres et al. (2015): Kepler-438b, Kepler-440b, and Kepler-442b. The repurposed Kepler Mission, dubbed K2, has also added a small planet on the inside edge of its habitable zone: K2-3d (aka EPIC 201367065d; Crossfield et al. 2015). Most of these small planets are beyond the reach of Doppler surveys and cannot be confirmed independently by radial velocity detections, and therefore must be validated statistically. Follow-up observations are conducted to rule out, to a very high statistical significance, scenarios in which astrophysical signals can mimic the planetary transits under investigation. Doppler surveys have been successful for a handful of very bright (nearby), cool stars (where the habitable zone is close-in and reflex velocity amplitude is high). These possibly rocky planets ( ) include GJ 163c (Bonfils et al. 2013), GJ 667 Cc (Anglada-Escudé et al. 2013; Delfosse et al. 2013), HD 40307g (Tuomi et al. 2013), Kapteyn b (Anglada-Escudé et al. 2014), and GJ 832c (Wittenmyer et al. 2014). All of these are orbiting small ( ), cool ( K) stars. According to the recent Kepler Q1–Q16 catalog (Mullally et al. 2015) available at the NExSCI exoplanet archive,14 approximately 65 additional small planetary candidates have been identified orbiting in or near the habitable zones of stars over a much wider range of spectral types with effective temperatures from 2703 K to as high as 6640 K. All of these candidates need to be vetted thoroughly, given their importance to the prospect of reaching the primary Kepler Mission objective: to determine the frequency and distribution of terrestrial planets in the habitable zones of Sun-like stars. Here we present the detection and statistical validation of Kepler-452b, a possibly rocky 1.6 planet orbiting in the habitable zone of its G2 host star KIC 8311864 every 384.84 days, the longest orbital period for a small ( ) transiting exoplanet discovered to date. In this paper, we apply the BLENDER analysis technique (see Section 9) to validate the planetary nature of Kepler-452b, but note that similar approaches have been developed by others and are in use in transit surveys (see, e.g., Morton 2012; Barclay et al. 2013; Díaz et al. 2014). The paper is organized as follows: Section 2 describes the initial discovery of the transit signature of this super Earth in 2014 May. Section 3 presents the photometry for this transiting planet, describing the discovery of the transit signature in the data, and Section 4 describes the results of several tests performed in the Kepler science pipeline to help validate the planetary nature of the photometric signature. Section 5 describes the spectroscopic measurements obtained for this host star, and Section 6 presents the stellar properties obtained for KIC 8311864. The analytical approach used to fit the light curve and determine the planetary properties is described in Section 7. Section 8 details the adaptive optics observations made to search for close companions. Section 9 describes the results of a full BLENDER validation analysis, which indicate that the planet hypothesis is favored at an odds ratio of 424:1, providing high confidence that this is, indeed, a planetary signature. Section 10 discusses the history of Kepler-452b in the context of the evolving habitable zone of its host star. Section 11 summarizes the results. The Appendix presents a mathematical development of the bootstrap diagnostic calculation.

2. DISCOVERY

Kepler-452b was discovered in a test run of the Kepler Science Operations Center (SOC) 9.2 codebase in 2014 May when one of us (J. Twicken) inspected the planet search pipeline results to assess performance of an enhanced pipeline codebase for small, cool planets. According to the Data Validation pipeline module, the transit signature of this object featured four 10.5 hr, 200 ppm deep transits spaced 384.846 days apart, a radius of 1.1 , and an equilibrium temperature of 221 K assuming an albedo of 0.3 and full redistribution of heat. None of the Data Validation diagnostics definitively rejected the planetary nature of this signature, and while both the radius and the equilibrium temperature proved to be significantly larger once better stellar parameters were obtained, this was clearly a highly interesting threshold crossing event (TCE). This transit signature was not identified in the previous search through all 4 years of Kepler data with the SOC 9.1 codebase conducted in 2013 August (Tenenbaum et al. 2014), most likely due to overly aggressive consistency checks put into place to reduce the large number of false alarms by instrumental effects (Seader et al. 2013). After all 4 years of Kepler data had been re-processed with the deployed SOC 9.2 codebase, Kepler-452b's transit signature was re-identified in the recent run of the planet search pipeline in 2014 November (Seader et al. 2015). The re-processed light curves and calibrated pixels have been archived at the Mikulski Archive for Space Telescopes (MAST) as data release 24, and the transit search results have been archived at NASA's Exoplanet Science Institute's exoplanet archive.15 This planet has been designated as KOI-7016.01 and appears in the online KOI catalog table on NExScI with the first delivery of the Q1–Q17 KOIs. We expect that several more terrestrial planetary candidates will be identified when all the new SOC 9.2 TCEs are fully vetted. This is the first planetary signature detected for KIC 8311864, which is noteworthy, given the long orbital period and the high fraction of multiple transiting planet systems identified by Kepler (see e.g., Rowe et al. 2014). If there are planets in orbits interior to that of Kepler-452b, either they are inclined so as not to present transits or they are too small to detect against the photometric noise for this 13.4 mag star (which averaged ~40 ppm at 10.5 hr timescales across the 4 year data set). The stellar parameters for this star at the time of the test run included an effective temperature of 5578 K but a radius of only 0.79 , yielding a planet with a correspondingly small radius of 1.1 at a mean separation of 0.99 AU. The small stellar radius was due to an anomalously high of 4.58.16 Given that the star was likely to be significantly larger than 0.79 on account of its temperature, we sought reconnaissance spectra in order to better pin down the stellar parameters and eventually obtained a spectrum with the Keck HIRES spectrometer. These measurements indicated that the surface gravity of the star was overstated by the KIC, and thus imply that the planet is ~45% larger than it appeared in the Data Validation summary report produced in the 2014 May test run of the Kepler pipeline.

3. KEPLER PHOTOMETRY

The Kepler telescope obtained nearly continuous photometry in a 115 deg2 field of view (FOV) near Cygnus and Lyra for nearly 4 years (Borucki et al. 2010). During this period (2009 May 13 through 2013 May 11), Kepler observed 111,800 stars nearly continuously and altogether observed a total of more than 200,000 stars over its operational mission. The Kepler observations were broken nominally into 3 month (93 day) segments, dubbed quarters, due to the need to rotate the spacecraft by 90° to keep its solar arrays pointed toward the Sun and its sunshade properly oriented (Koch et al. 2010). The star Kepler-452 (Kepler Input Catalog (KIC) = 8311864, α = , δ = , J2000, KIC mag) was observed during each of the 17 quarters (Q1, Q2, ..., Q17) returned by Kepler.17 The first and last quarters were very much abbreviated compared to the others: Q1 lasted 34 days and took place immediately after commissioning was complete (Haas et al. 2010), while Q17 lasted only 32.4 days, being truncated by the failure of a second reaction wheel, thus ending the data collection phase of the Kepler Mission. Figure 1 shows all 17 quarters of the systematic error-corrected flux time series data for this star, with the locations of the transits identified by triangle markers.18 The first transit takes place during Q3 at a barycentric-corrected Kepler Julian date (BKJD)19 of , and three additional transits occur every days, one each in Q7, Q11 and Q15. The transits each last approximately 10.5 hr and are 200 ppm deep. The scatter in the photometric data in Figure 1 is quite uniform in all quarters and there is no obvious evidence of the electronic image artifacts that plague a subset of the Kepler target stars as they rotate through the handful of affected CCD readout channels, resulting in a large number of false alarm TCEs (van Cleve & Caldwell 2009; Tenenbaum et al. 2013, 2014). Indeed, this star does not fall on any of the CCD channels associated with these image artifacts in any observing season. Figure 1. Detrended flux time series for Kepler-452 vs. barycentric-corrected Kepler Julian day (BKJD—2,454,833.0). The triangle symbols indicate the locations of each of the four transits which occur in quarters 3, 7, 11 and 15. All 4 years of Kepler data are plotted here with vertical dashed red lines indicating the boundaries between consecutive quarters. The detrending was accomplished by subtracting the result of passing the systematic error-corrected light curve through a moving median filter with a window width of five times the ~10.5 hr pulse duration. The light curve residuals are remarkably uniform across all 17 quarterly data segments, showing no evidence of electronic image artifacts which are responsible for a large number of long period false positive transit-like features in the Kepler light curves. Download figure: Standard image High-resolution image Export PowerPoint slide The planet was detected at the level, well above the threshold and passed the other consistency checks exercised by the Transiting Planet Search (TPS) component of the SOC pipeline (Jenkins 2002; Jenkins et al. 2010; Seader et al. 2013; Tenenbaum et al. 2013, 2014). The light curve with its putative transit signature was fitted with a Mandel & Agol (2002) limb-darkened transit waveform and subjected to a suite of diagnostic tests by the Data Validation module (Tenenbaum et al. 2010; Wu et al. 2010). Figure 2 shows the phase-folded, detrended light curve, zoomed on the transit epoch. The individual transits are significant at the ~4.75 level, so that each can be credibly seen by eye in the long cadence (LC) data at a 29.4 minute resolution. The DV diagnostics appeared to be consistent with the planetary hypothesis, so we coordinated follow-up observations to help refine the stellar (and hence planetary) properties and to further bolster the case for the planetary nature of this object. Figure 2. Phase-folded, detrended flux time series for Kepler-452 vs. phase in hours. The flux time series plotted here is the same as that of Figure 1, except that it has been phase folded at the period of the planet, 384.84 days. The triangle symbol indicates the location of the four transits at zero phase. The full resolution data are represented by black dots. The cyan filled circles represent the results of binning the data in phase with a bin width of 1/5 of the DV-fitted duration of the transit, 10.5 hr. Download figure: Standard image High-resolution image Export PowerPoint slide The following section describes the results of the validation tests performed in DV.

4. DATA VALIDATION TESTS

This section describes the most salient validation tests performed inside of the Data Validation component of the SOC pipeline. The DV report resulting from the official SOC 9.2 processing is available at http://exoplanetarchive.ipac.caltech.edu/data/KeplerData/008/008311/008311864/dv/kplr008311864-20141002224145_dvr.pdf.20 4.1. Discriminating Against Background Eclipsing Binaries Data validation applies several tests for the presence of background eclipsing binaries. The first is the odd/even test: DV fits a limb-darkened transit model to the odd transits and an independent model to the even transits, and then compares the depth fitted to the odd transits to that fitted to the even transits. This results in a test statistic of with a corresponding significance of 89.0%.21 In addition, DV examines the ratio of the temporal offsets between the odd and even transits compared to the fitted orbital period, the so-called "epoch" test, giving a result of and a corresponding significance of 92%. Figure 3 illustrates the odd/even test by plotting the light curve data phase folded separately for the odd and even transits. No other TCEs were generated for this star and no evidence was seen for weak secondary eclipses. Based on these tests, we see no evidence for systematic differences in the odd/even transit depths and durations which often result in the case of background eclipsing binaries, and the transit times are consistent with a strict periodic spacing. Figure 3. Phase-folded, detrended flux time series for Kepler-452 vs. phase in hours as in Figure 2, with the data for the even transits separated from that of the odd transits. The full resolution data are represented by black dots. The cyan fllled circles represent the results of binning the data in phase with a bin width of of the DV-fitted duration of the transit, 10.5 hr. The red horizontal lines are the fitted depths for the odd transits (left) and the even transits (right) and the dashed horizontal lines represent the error bars on the fitted depths. Download figure: Standard image High-resolution image Export PowerPoint slide 4.2. Statistical Bootstrap Test In order to better determine the statistical significance of the transit sequence, we examined the bootstrap distribution for the population of null statistics for the out-of-transit data (Jenkins et al. 2002; J. Jenkins & S. Seader 2015, in preparation). This test relaxes the assumption that the observation noise is broadband colored Gaussian noise to establish a robust significance level for the detection statistic of the planet candidate. The Appendix contains a mathematical derivation of the statistical approach and a brief description of its implementation. The false alarm rate is estimated using a bootstrap statistical analysis of the out-of-transit data for this star. The single event statistics time series, representing the likelihood of a transit of the fitted duration as a function of time, are drawn at random (with replacement) to construct artificial multiple event statistics corresponding to the same number of transits that occurred for the TCE in question. The distribution of all such bootstrap samples provides a means to assess the statistical significance of the maximum multiple event statistic, the maxMES, that resulted in the identification of the TCE in the first place. This methodology relaxes the assumption of perfection on the part of the adaptive, power spectral density estimator in TPS that informs the pre-whitener, allowing for the detection signal-to-noise ratio (S/N) to be evaluated against the distribution of null statistics drawn directly from the observation noise of the light curve. The bootstrap false alarm rate for the transit signature of Kepler-452b is securely established as statistically significant by this test for the detection of this planet, at a level of . This figure is well below that required by the Kepler experiment, , to limit the number of false positives to no more than one over the full mission. Figure 4 illustrates the bootstrap test for this planet. Figure 4. Bootstrap false alarm rate (FAR) vs. maximum multiple event statistic for the light curve for Kepler-452 (blue curve), and that expected for ideal Gaussian noise fitted to the empirical FAR between 1 × 10−4 and 1 × 10−13 (dashed green curve). The Gaussian fit allows us to extrapolate the FAR out to the maximum MES provided by the planetary signature as the highest MES probed by the bootstrap distribution falls short of the 9.7σ MES of this transiting planet in the Kepler data. The red star indicates the estimated bootstrap FAR obtained from the Gaussian fit. Download figure: Standard image High-resolution image Export PowerPoint slide 4.3. Centroid Analysis Data validation performs two types of pixel-level centroid analysis (Bryson et al. 2013) to determine the location of the transit source relative to KIC 8311864. This analysis constrains the likelihood that transits observed in a light curve may be due to background eclipsing binaries by providing a "confusion radius," inside which a background transit source cannot be distinguished from the target star. Sources outside this confusion radius are ruled out as possible causes of the transit signal. This radius is passed to the BLENDER analysis of Section 9. The first analysis technique uses the change in the photocenter (centroid shift) of the flux on pixels when the transit occurs to infer the location of the source of the transit signal. The row and column of the photocenter of pixels collected for this target are measured for every LC in which the target was observed. These row and column time series are then fitted to the transit model (normalized to unit depth) derived from the light curve. The fitted amplitude of the row and column fits provides the apparent centroid shift on the sky. The use of every observed LC in this fit allows a very high precision measurement of this centroid shift. Dividing this centroid shift by the transit depth gives an estimate for the transit source location relative to the target star. For KIC 8311864, the source offset was estimated as 1 832 at , consistent with the transit source being co-located with KIC 8311864. This estimate, however, tends to be biased so as to overestimate the source offset relative to the mean photocenter of the pixels. It also tends to be more vulnerable to noise, instrumental systematics, and other flux variations in the photometric mask than estimates obtained via the second technique. The second method consists of performing difference image analysis on in-transit pixel data compared with neighboring out-of-transit pixel data in each quarter. The average of the in-transit images is subtracted from the averaged out-of-transit images. When flux variation in the photometric mask other than the transit is negligible, the transit source appears as a star in the difference image. Fitting the quarterly difference and out-of-transit images with the Kepler point-spread function (PSF) then provides a direct measurement of the offset of the transit source from KIC 8311864 per quarter. The quarterly offsets are combined via robust averaging to provide an estimate of transit source offset. The difference image analysis generally tends to be less vulnerable to noise and other flux variation than the brightness weighted centroid results. For shallow transits which occur only once per quarter, such as Kepler-452b, however, instrumental systematics can result in very poor (un-star-like) difference images whose PSF fit will result in significantly biased position estimates. When the systematics are sufficiently strong, the fit can fail to converge. For KIC 8311864, all four quarters containing transits had successful fits, shown as the green crosses in Figure 5. Two of the quarters, however, had very poor, non-star-like difference images and resulted in the two locations to the northwest of KIC 8311864. The other two quarters had good, star-like images, resulting in offsets that straddle KIC 8311864, slightly offset to the west. Because the two non-star-like difference images are not validly fit by the Kepler PSF (the difference images in these quarters do not contain anything that looks like a star), we exclude them from our analysis. The remaining two quarters give a transit source offset from KIC 8311864 of 0.782 ± 0 33, or . This measurement formally has a confusion radius of 1 0. Figure 5. Difference image centroiding results for Kepler-452b with a J-band infrared image from the UK Infrared Telescope Survey (UKIRT; Lawrence et al. 2007) overlaid. There is a faint, 18th mag star offset by ~4'' to the west of KIC 8311864, but neither the difference image centroids in this figure, nor the centroid shift analysis provide any evidence suggesting that this star is the source of the transit-like features. The four green crosses indicate the difference image centroids for each of the four transits of this super Earth. The blue circle delimits the 3σ confusion radius, which contains the target star, KIC 8311864, indicating that the target star is the most likely source of the transit features in the Kepler light curve, to the limits of the resolution of the UKIRT images. Download figure: Standard image High-resolution image Export PowerPoint slide To address the concern that a measurement using only two data points can be misleading, we estimate that a Gaussian uncertainty of 0 50 would result in the observed uncertainty of 0 33 half the time. We therefore conservatively set the 3σ confusion radius to 1 5. Potential transit sources outside this 1 5 radius are ruled out. Possible unknown sources inside 1 5 are addressed in the BLENDER analysis of Section 9. 4.4. Contamination by Other Sources In addition to the tests conducted by the pipeline module Data Validation on the Kepler data themselves, we also performed manual checks against contamination by other astrophysical sources in the Kepler FOV. The two principal sources of contamination are (1) optical ghosting from eclipsing binaries well offset from the target star due to reflections off of optical elements in the Kepler telescope, and (2) video crosstalk between CCD channels that are read out simultaneously. Optical ghosting, or direct PSF contamination, was the subject of a study by Coughlin et al. (2014) identifying 685 TCEs in the Q1–Q12 vetting process as false positives due to these and other effects. To determine whether optical ghosting was an issue, we looked for period-epoch collisions between Kepler-452b and the TCEs identified in both the 2012 August search through Q1–Q16 (Tenenbaum et al. 2013) and those from the recent 2014 September search through Q1–Q17 (Seader et al. 2015). None were found. Nor were any TCEs found either at one half or one third of the orbital period of Kepler-452b with transit epochs matching the epochs of Kepler-452b. In addition, we queried the general catalogue of variable stars (GCVS) (Samus et al. 2009) for variable stars with periods within 10 hr of Kepler-452b. Only two stars in the GCVS had periods close to 384.84 days that were within ~6° of KIC 8311864. Both of these are pulsating variables, however, not eclipsing binaries, and therefore cannot explain the transits of Kepler-452b. One of them is FX CYG aka KIC 4671784, which was observed in LC mode throughout Kepler's mission. The other variable star is V1299 Cyg which is more than from the Kepler FOV and thus cannot be the source of any contamination. No variable stars were found in the GCVS in or near Kepler's FOV with periods near one half or one third that of Kepler-452b. The likelihood that optical ghosting could explain the transits of KIC 8311864 is further diminished by the fact that most optical ghosts tend to create multiple "children" signals in the Kepler light curves as the ghost images are defocussed or severely scattered. Video crosstalk occurs when analog voltage signals in one video circuit impress themselves onto other video circuits in the readout electronics boards (van Cleve & Caldwell 2009, p. 71), contaminating the other CCD readout values. This represents an additive noise source and the contamination is an attenuated version of the original signal, but can either be positive or negative, with a typical value of ~4400 for the attenuation factor (or a gain factor of ). While video crosstalk does occur between channels read out at the same time on different CCD modules, the attenuation is so high that the only appreciable video crosstalk occurs within each pair of CCDs mounted together on the same module. To rule out contamination by video crosstalk, we examined the quarter 3 month 2 full frame image (FFI) kplr2009322233047_ffi-cal.fits obtained on UT 2009 November 18 for evidence of bright stars read out on the same pixels as KIC 8311864 on the adjacent CCD channels.22 For a target star of brightness B experiencing video crosstalk from a perfectly situated eclipsing binary (or transiting giant planet) of brightness b and an eclipse (or transit) depth of on an adjacent CCD readout channel, the apparent transit depth exhibited by the target star is given by where the second line arises for the case of interest here as . All four transits of KIC 8311864 occur on CCD module/output 7.3, as the orbital period is close to one year and the first transit occurs sufficiently early in Q3 so that the subsequent three transits take place when the Kepler photometer is in the same seasonal orientation. Thus, only stars on the adjacent channels 7.1, 7.2, and 7.4 could produce video crosstalk contamination for Kepler-452b. Of these, only channel 7.1 has a readily visible star ( , ) near the same row/column coordinates as for KIC 8311864 on channel 7.3. Its peak pixel actually is just 1 row below the optimal aperture for Q3, so only flux from the wing of the PSF could be getting into the aperture of KIC 8311864. No Kepler light curves exist for this star since it was not selected as a target during the 4 year Kepler Mission. This star on channel 7.1 is about 3 mag fainter than KIC 8311864 and the video crosstalk coefficient for channel 7.1 on channel 7.3 is . Hence, a full eclipse (50% drop in brightness) on this star on channel 7.1 would produce an apparent transit depth of ~0.75 ppm. A 1% giant planet transit on this star would impute an even smaller apparent transit depth of ~0.015 ppm. We can thus eliminate this star as a potential source of confusion as it cannot credibly cause a transit-like feature deep enough to explain the 199 ppm transits exhibited by Kepler-452b. Channels 7.2 and 7.4 are also incapable of producing a false positive here through video crosstalk. Channel 7.2 has a negative crosstalk coefficient with 7.3 (−1.9759 ×10−4) and 7.4 has a crosstalk coefficient of , but no apparent star in the pixels corresponding to the aperture of KIC 8311864. We conclude that contamination through either optical ghosting or video crosstalk is highly unlikely for Kepler-452b.

5. SPECTROSCOPY

Due to the curiously small stellar radius of KIC 8311864 for the effective temperature available at the time of the 2014 May planet search, it was clear that reconnaissance spectroscopy was called for in order to obtain a reasonable interpretation of the nature of the source of the transit-like features in the Kepler light curve. We obtained one spectrum from the McDonald Observatory in 2014 May and two additional spectra from the Fred L. Whipple Observatory in 2014 June. To improve the stellar parameters derived from spectroscopy, we obtained a spectrum from the Keck I observatory with the HIRES instrument in 2014 July. The HIRES data were analyzed with three independent analysis packages in order to better understand the effects of systematics in the software codes and to determine the most reasonable stellar parameters to use for interpreting the photometry in this system. 5.1. Reconnaisance Spectroscopy with the Tull spectrograph at the McDonald Observatory 2.7 m Telescope We obtained a reconnaissance spectrum for KIC 8311864 using the Tull Coudé spectrograph (Tull et al. 1995) at the Harlan J. Smith 2.7 m telescope at McDonald Observatory. We used the standard instrumental setup that covers the entire optical spectrum and uses a 1 2 slit that yields a resolving power of R = 60,000 with 2 pixel sampling. We observed the star in the night of UT 2014, May 18. The exposure time was 1720 s and the resulting spectrum has a S/N of around 31:1 (per resolution element) at 5650 Å. We retrieved stellar parameters using our new spectral fitting tool Kea. Kea was specifically developed to determine reliable stellar parameters from the lower S/N reconnaissance spectra of our Kepler follow-up program at McDonald Observatory. We used the MOOG (Sneden 1973) stellar spectral code in its spectral synthesis mode with the Kurucz (1993) stellar atmospheric models to construct a grid of synthetic stellar spectra. We used atomic line parameters obtained from the Vienna Atomic Line Database (VALD). We included molecular opacities for MgH (Hinkle et al. 2013), TiO (Plez 1998) and CN (Sneden et al. 2014). The spectral grid covers a range of effective temperature from 3500 to 10,000 K in 100 K steps, a stellar metallicity of [Fe/H] from −1.0 to +0.5 dex in 0.25 dex steps, and a surface gravity log (g) of 1.0 to 5.0 dex in steps of 0.25 dex. We used the MOOG "weedout" feature to eliminate atomic and molecular lines in the line list with a ratio of line to continuum opacity of less than 0.001. Kea compares 12 orders of the continuum-normalized Tull spectrum with the synthetic grid using statistics. For this comparison we convolve the synthetic spectra with a Gaussian-shaped PSF to assure the same spectral resolution of model and data. We also rotationally broaden the model spectra to find a best match for the (projected) of the star. The final stellar parameters are obtained from the mean and standard deviation of the best-fit models of the spectral orders that Kea is using. A more detailed description of this method will be presented in a forthcoming publication (M. Endl & W. Cochran 2015, in preparation). From the KIC 8311864 reconnaissance spectrum, we determined an effective temperature of K, a metallicity of dex, a surface gravity of , and a rotational broadening of km s−1 for KIC 8311864. 5.2. Reconnaisance Spectra Obtained from TRES We obtained two spectra using the Tillinghast Reflector Echelle Spectrograph (TRES; Fűresz 2008), on the 1.5 m telescope at the Fred Lawrence Whipple Observatory (FLWO) on Mt. Hopkins, Arizona, which has a resolving power of R = 44,000 and a wavelength coverage of ~3900–9100 Å. The spectra were taken on UT 2014, June 09 and UT 2014, June 14 and were extracted as described by Buchhave et al. (2010). We used the Spectral Parameter Classification (SPC) (Buchhave et al. 2012) technique to determine stellar parameters. SPC cross correlates an observed spectrum against a grid of synthetic spectra using the correlation peak heights to determine the best fit parameters. The synthetic spectra are based on Kurucz model atmospheres (Kurucz 1993). We allowed , , [Fe/H], and to float as free parameters. The S/N per resolution element on the TRES spectra were slightly lower than usual when running SPC so these results should be viewed with caution. A weighted average of the results from the two spectra gave = 5751 ± 55 K, = 4.43 ± 0.10, [Fe/H] = 0.40 ± 0.08, and = 4.2 ± 0.5 km s−1. The maximum height of the normalized cross correlation function (CCF) was 0.916 and the combined S/N of the spectra was 23.5. 5.3. HIRES Spectroscopy On UT 2014, July 3, we acquired a high resolution spectrum of KIC 8311864 with the Keck I Telescope and HIRES spectrometer. We used the standard setup of the California Planet Search (Howard et al. 2010) without the iodine cell, resulting in a spectral resolution of 60,000, and S/N, at 5500 Å, of 40 per pixel. SpecMatch compares a target stellar spectrum to a library of 800 spectra from stars that span the HR diagram ( = 3500–7500 K; = 2.0–5.0). Parameters for the library stars are determined from LTE spectral modeling. Once the target spectrum and library spectrum are placed on the same wavelength scale, SpecMatch computes the sum of the squares of the pixel-by-pixel differences in normalized intensity. The weighted mean of the ten spectra with the lowest values is taken as the final value for the effective temperature, stellar surface gravity, and metallicity. SpecMatch-derived stellar radii are uncertain to 10% rms, based on tests of stars having known radii from high resolution spectroscopy and asteroseismology (Huber et al. 2013). The stellar parameters determined using SpecMatch were = 5757 ± 60 K, = 4.32 ± 0.07 dex and [Fe/H] = 0.21 ± 0.04. Two other analyses were carried out on the HIRES spectrum to better understand the role of instrumental systematic errors between the three sets of observations, and to better determine how to set the error bars on the resulting stellar parameters. An SPC analysis of the HIRES spectrum gave the following stellar parameters: = 5740 ± 50 K, = 4.42 ± 0.10, [Fe/H] = 0.22 ± 0.08, = 0.8 ± 0.5 km s−1, maximum height of the CCF = 0.988, and S/N = 58.8 per resolution element. We also performed a standard EW analysis of Fe i and Fe ii lines using MOOG and Kurucz's odfnew models. We determined two sets of , , [Fe/H], and (microturbulence) values, one using only the Keck spectrum and a second using a solar spectrum from MIKE/Magellan for line-by-line differential analysis. The solution based solely on the Keck spectrum gave = 5740 ± 35 K, = 4.230 ± 0.072, A(Fe) = 7.649 ± 0.033 (which would correspond to [Fe/H] = +0.2 for an absolute solar iron abundance of A(Fe) = 7.45, as inferred by numerous 1D-LTE analyses), and = 1.03 ± 0.08 km s−1. The solution with the Keck spectrum + MIKE/Magellan solar spectrum template gave = 5818 ± 21 K, = 4.33 ± 0.05, [Fe/H] = 0.24 ± 0.02, and micro turbulent velocities, = 1.07 ± 0.05 km s−1. As an additional search for unseen companions, we searched the HIRES spectrum for secondary spectral lines caused by a second star that is angularly close enough to the primary to fall into the HIRES spectrometer slit (0 87 wide). The observed spectrum is first compared to a library of representative main sequence stars. The best match library spectrum is subtracted from the observed spectrum, and the residuals are then inspected for evidence of a secondary spectrum (Kolbl et al. 2015). We found no secondary spectrum brighter than 1% of the primary, separated in radial velocity by more than 10 km s−1. Stars separated by less than 10 km s−1 cannot be detected because they overlap with the primary spectrum. See Section 9 for more details. Table 1 contains the results for the spectroscopy of KIC 8311864 and attendant analyses. We considered the scatter and precision in the spectroscopic results in order to determine the most reasonable stellar parameters for the recovery of planetary parameters, as described in the next section. Table 1. Spectroscopy and Analysis Results for KIC 8311864 Parameter Value Notes Tull/McDonald Observatory Effective temperature (K) 5650 ± 108 a Surface gravity 4.45 ± 0.16 a Metallicity [Fe/H](dex) 0.21 ± 0.07 a Projected rotation Velocity −1) 4.3 ± 0.5 a TRES/Whipple Observatory Effective temperature (K) 5751 ± 55 b Surface gravity 4.43 ± 0.10 b Metallicity [Fe/H](dex) 0.40 ± 0.08 b Projected rotation Velocity −1) 4.2 ± 0.5 b HIRES/Keck Observatory Effective temperature (K) 5757 ± 60 c Surface gravity 4.32 ± 0.07 c Metallicity [Fe/H](dex) 0.21 ± 0.04 c Projected rotation Velocity −1) <1 c Effective temperature (K) 5740 ± 50 b Surface gravity 4.42 ± 0.10 b Metallicity [Fe/H](dex) 0.22 ± 0.08 b Projected rotation Velocity −1) 0.8 ± 0.5 b Effective temperature (K) 5818 ± 21 d Surface gravity 4.33 ± 0.05 d Metallicity [Fe/H](dex) 0.24 ± 0.02 d Notes. a: Analyzed with Kea. b: Analyzed with SPC. c: Analyzed with SpecMatch. d: Analyzed with MOOG. Download table as: ASCIITypeset image

6. HOST STAR PROPERTIES

Since KIC 8311864 is too faint for direct observations such as interferometry or asteroseismology, the derivation of stellar mass, radius and density relies on matching atmospheric properties to interior models. Figure 6 compares the spectroscopic solutions described in the previous section. The values are generally within 1 σ and no strong correlation of with or [Fe/H] is apparent, indicating that the differences are dominated by differences in the methods rather than degeneracies in atmospheric properties (e.g., Torres et al. 2012). On average, KIC 8311864 is slightly cooler, slightly larger, and about 60% more metal-rich than the Sun. Figure 6. Spectroscopically derived surface gravity as a function of effective temperature (top panel) and metallicity (bottom panel) for KIC 8311864, derived using different methods and spectra (see text for details). The position of the Sun is also shown. Download figure: Standard image High-resolution image Export PowerPoint slide To ensure self-consistency, we adopted a single spectroscopic solution for the remainder of the paper. We chose the SpecMatch solution, which yields an intermediate temperature, conservative surface gravity, and metallicity consistent with most other methods. To account for the systematic differences between methods, we added the standard deviation of all methods for a given parameter in quadrature to the formal errors provided by SpecMatch. The final adopted atmospheric properties are listed in Table 2. Table 2. Stellar Properties of KIC 8311864 Parameter Value Derived from high-resolution spectroscopy using SpecMatch Effective Temperature (K) 5757 ± 85 log [Surface gravity] (dex) 4.32 ± 0.09 Metallicity [Fe/H] (dex) 0.21 ± 0.09 Derived from fitting Radius ( Mass ( Mean Density (g cm−3) Age (Gyr) Download table as: ASCIITypeset image To derive interior properties, we fitted , , and [Fe/H] to a grid of Dartmouth isochrones (Dotter et al. 2008). The grid was calculated using solar-scaled alpha-element abundances, and interpolated to a stepsize of 0.01 in mass and 0.02 dex metallicity. We fitted the input values to the model grid using a Markov chain Monte Carlo (MCMC) algorithm which, at each step, draws samples in mass, age, and metallicity. For ease of computation, the samples in age and metallicity were drawn in discrete steps corresponding to the sampling of the model grid (0.5 Gyr in age and 0.02 dex in [Fe/H]). For each sample of fixed age and metallicity, we linearly interpolated the grid in mass to derive and , which were then used to evaluate the likelihood functions. We adopted uniform priors in mass, age, and metallicity. We calculated 106 iterations and discarded the first 10% as burn-in. Figure 7 shows the marginalized and full posterior distributions, as well as the position of KIC 8311864 on a radius-temperature diagram. As expected, KIC 8311864 is slightly evolved off the zero-age main-sequence. The 68% confidence interval for the stellar radius ranges from 1.02–1.26 , with a most probable value of 1.11 . Figure 7. Marginalized and full posterior distributions of Markov chain Monte Carlo chains based on using spectroscopic , and [Fe/H] to constrain a grid of Dartmouth isochrones. A random set of 750 samples (~0.1% of all iterations) is plotted for each parameter for clarity. Vertical dashed and dotted lines mark the posterior mode and 68% confidence interval around the mode for each parameter. Note that age and metallicity are discretely sampled with the resolution of the model grid. The top right panel shows the position of KIC 8311864 in a radius temperature diagram, with solid black lines showing evolutionary tracks with 0.9–1.2 in steps of 0.05 for [Fe/H] = 0.21. The gray dashed line marks the zero-age main-sequence. Download figure: Standard image High-resolution image Export PowerPoint slide As a consistency check, we also computed the age, mass, and radius of this star based on the MOOG analysis of the HIRES spectrum using Yonsei-Yale isochrones (Demarque et al. 2004; Yi et al. 2004) and the approach detailed in Section 4.5 of Ramírez et al. (2014), obtaining an age of 5.2 (3.4–7.5) Gyr, a mass of 1.07 (1.02–1.12) , and a radius of 1.12 (1.03–1.29) . The values in parentheses represent approximate ranges of values. Note that the agreement in these derived quantities with those derived for the SpecMatch results is very high, giving us confidence that we have a reasonable degree of control over the systematics in the different analysis techniques.

7. PROPERTIES OF THE PLANET KEPLER-452B

To estimate the posterior distribution on each fitted parameter in our transit model, we used an MCMC approach similar to the procedure outlined in Ford (2005) and implemented in Rowe et al. (2014). Our algorithm used both a Gibbs sampler and a control set of parameters to enable a vector jumping distribution to efficiently handle the sampling of correlated parameters as outlined in Gregory (2011). We directly fitted the following transit and orbital parameters: orbital period (P), epoch, scaled planet radius ( ), impact parameter (b), and eccentricity parameters ( and ). We generated ten sets of 1,000,000-element Markov chains, then discarded the first 20% of each chain to establish burn-in. The remaining sets were combined and used to calculate the mode, standard deviation, and bounds of the distribution centered on the mode of each model parameter. Our model fits and uncertainties are reported in Table 3. The eccentricity of Kepler-452b was not well constrained by the light curve as was expected in the absence of a radial velocity detection of the planet. Table 3. Planet Parameters for Kepler-452b Parameter Value Notes Transit and orbital parameters Orbital period P (day) a, b Epoch (BJD—2454833) a, b Scaled planet radius a, b Impact parameter a, b Orbital inclination i (deg) a Transit depth a Transit duration a Eccentricity a, b Eccentricity a, b Planetary parameters Radius a Orbital semimajor axis a (AU) a Equilibrium temperature c Insolation relative to Earth d Notes. a: Based on the photometry. b: Directly fitted parameter. c: Assumes Bond albedo = 0.3 and complete redistribution. d: Based on Dartmouth isochrones. Download table as: ASCIITypeset image We used the Markov chains to derive model dependent measurements of the transit depth, , and transit duration, . The transit depth is defined as the transit-model value when the projected distance between the star and planet is minimized. The transit duration is the time from first to last contact. We also convolved the transit model parameters with the stellar parameters (see Section 6) to compute the planetary radius, = 1.6 , and the flux received by the planet relative to the Earth ( ). We also performed a fit forcing the orbit to be circular in our transit model and obtained thereby an estimate for the density of the star of = 1.1 ± 0.3 g cm−3. While this value is perfectly consistent with the spectroscopic work, it suggests that both the star and the planet may be somewhat smaller than indicated by the SpecMatch analysis.

8. ADAPTIVE OPTICS OBSERVATIONS

KIC 8311864 has one companion visible in the UKIRT image (see Figure 5) at a separation of ~4''. This dim = 18.435 star is well outside the 3 σ confusion radius and thus cannot be the source of the transit-like features in the Kepler photometry, according to the centroid arguments presented in Section 4.3. However, it is important to probe the scene around KIC 8311864 for evidence of nearby companions that would not be identifiable at the spatial resolutions in either the Kepler or the UKIRT image data. We obtained near-infrared adaptive optics (AO) images in the J (1.248 μm) and Br-γ (2.157 μm) filters using the NIRC2 imager behind the natural guide star AO system (Wizinowich et al. 2004; Johansson et al. 2008) at the Keck II telescope on UT 2014, June 13. Data were acquired in a three-point dither pattern to avoid the lower left corner of the NIRC2 array, which is noisier than the other three quadrants. The dither pattern was performed three times with step sizes of 3 0, 2 5, and 2 0 for a total of nine frames; individual frames had an integration time of 60 s (30 s times 2 coadds) in both bands. The individual frames were flat-fielded, sky-subtracted (sky frames were constructed from a median average of the 9 individual source frames), and co-added. The average AO-corrected seeing (FWHM) was 0 04635 and 0 05229 at J- and K-bands, respectively. Figure 8 shows the Keck AO images obtained for this star. KIC 8311864 appears to be a well isolated star in both bands with no obvious background stars, or resolved, associated companions within ±2'' that would present blends at the 4'' pixel scale of Kepler's camera. Figure 9 shows the 5 σ sensitivity limits for rejecting the presence of background confusion sources as a function of radial distance from KIC 8311864. The sensitivity achieved in both bands establishes rejection of background stellar sources from ~0 5 to 2 0 within mag of the central target star, in the bands of observation. Using typical stellar Kepler-infrared colors (Howell et al. 2012), the limitations of the observations translate into the Kepler bandpass sensitivity limits near mag. Given that the transits of Kepler-452b are 199 ppm deep, a background eclipsing binary exhibiting 50% eclipses would need to be within 8.66 mag of Kepler-452. The AO images and the difference image analysis carried out with the Kepler data allow us to validate the planetary nature of the transit source with high confidence, as detailed in the next section. Figure 8. Images obtained with the NIRC instrument at the Keck II telescope on 2014 June 13 in J band (left) and K band (right). Neither band indicates the presence of any background stars that might be the source of the transits observed in the Kepler photometry of Kepler-452. Download figure: Standard image High-resolution image Export PowerPoint slide Figure 9. 5 σ sensitivity limits on the rejection of background confusion sources in terms of as a function of radial offset from KIC 8311864. Left: the dashed–dotted blue curve is the measured in-band sensitivities at , while the solid green curve indicates the corresponding estimated sensitivity limits converted to the Kepler bandpass. Right: same sensitivity curves for . Download figure: Standard image High-resolution image Export PowerPoint slide

9. BLENDER ANALYSIS

In this section, we examine the possibility that the Kepler-452b signal may be due to an astrophysical false positive. Typical examples of such phenomena include a background or foreground eclipsing binary along the line of sight to the target (i.e., within the photometric aperture of Kepler), a chance alignment with a background or foreground single star transited by a (larger) planet, or a physically associated star transited by a planetary or stellar companion (hierarchical triple). We refer to the first and second of these scenarios as "BEB" and "BP," respectively, and to the third as "HTP" or "HTS," depending on whether the object producing the eclipses is a planet or a star.23 Other possible false positive scenarios include the transit of a white dwarf, or the occultation of a brown dwarf in an eccentric orbit with an orientation such that the transit does not occur (see, e.g., Santerne et al. 2013). The white dwarf scenario can be ruled out because it would result in a very large microlensing signal of ~2400 ppm, as per Sahu & Gilliland (2003), that is not observed. A non-transiting, occulting white dwarf is unlikely to produce a 200 ppm deep transit-like feature as the white dwarf would need to be about 7000 K to produce an occultation of the correct depth, and white dwarfs cool from 10,000 to 6000 K in Gyr (see, e.g., Goldsbury et al. 2012). Thus, there would be only a very brief opportunity for a white dwarf to provide 200 ppm deep occultations of KIC 8311864 as the uncertainty in the transit depth is ~10%. Further, we expect white dwarf/main sequence pairs at a period of 385 days to be relatively rare: Farmer & Agol (2003) predicted that only five microlensing white dwarf/main sequence binaries would be discovered in Kepler's FOV out to periods of ~100 days, and indeed, only one such white dwarf has been found to date (Kruse & Agol 2014) in an 88 day period orbit of a solar-like star exhibiting ~1000 ppm microlensing events and (coincidentally) ~1000 ppm deep secondary occupations. The brown dwarf scenario is very unlikely to begin with (much more so than the BEB, BP, or HTP/HTS scenarios) because of the known paucity of brown dwarfs as spectroscopic companions to main-sequence stars within hundreds of AU, and especially G-type stars (see, e.g., Marcy & Butler 2000), whereas binaries with main-sequence companions are far more common. Only eight transiting brown dwarf/main sequence binaries have been discovered to date (Moutou et al. 2013), four of which are in the Kepler FOV (Bouchy et al. 2011; Johnson et al. 2011; Díaz et al. 2013; Moutou et al. 2013). All but one of these have orbital periods under 20 days, and KOI-415b has a period of 167 days (Moutou et al. 2013). Further, even the most massive brown dwarf would cool below a luminosity ratio of 100 ppm with the Sun in 1.5 Gyr, and below 25 ppm by 6 Gyr (Burrows et al. 2001). Even if such a case were to occur, the shape of the transit-like signals (particularly the ingress/egress) would not match what we observe because of the much larger size of brown dwarfs compared to terrestrial planets like Kepler-452b. Thus we restrict our analysis to the above BEB, BP, and HTP/HTS configurations. The general procedure we followed consists of generating realistic transit-like light curves for a large number of these false positive scenarios ("blends"), and comparing each of them with the Kepler photometry for Kepler-452 to see which provide a fit indistinguishable from that of a model of a planet transiting the target star. This comparison allowed us to place constraints on the various physical parameters of the blends, including the sizes and masses of the bodies involved, their overall color and brightness, the linear distance between the background/foreground eclipsing pair and the target, etc. We then used a Monte Carlo approach to estimate the expected rate of occurrence of each type of blend (BEB, BP, HTP/HTS) from our simulations, and we compared the sum with the expected rate of occurrence of transiting planets similar to the one implied by the Kepler-452 signal. If the a priori planet occurrence is much larger than the anticipated total blend frequency, then we would consider the candidate to be statistically validated. The algorithm used to generate blend light curves and compare them to the observed photometry in a sense is known as BLENDER (Torres et al. 2004, 2011; Fressin et al. 2012), which has been applied in the past to the validation of many of the most interesting Kepler discoveries (see, e.g., Ballard et al. 2013; Barclay et al. 2013; Borucki et al. 2013; Meibom et al. 2013; Kipping et al. 2014), including, most recently, a sample of twelve small planets in the habitable zone of their parent stars (Torres et al. 2015). For full details of the technique, we refer the reader to the latter publication. That work also explains the calculation of the blend frequencies through a Monte Carlo procedure, making use of the constraints provided by BLENDER along with other constraints available from follow-up observations. These include the sensitivity to companions at close angular separations from our AO imaging (Section 8), limits on the brightness of unseen companions at even closer separations derived from our high-resolution spectroscopy (Section 5.3), constraints on blended objects from our difference image analysis (specifically, the size of the angular 3 σ exclusion region; Section 4.3), and also the brightness measurements for KIC 8311864 as listed in the KIC (which allowed us to rule out blends whose combined light is either too red or too blue compared to the observed color of the target). The time series used for the BLENDER analysis is the detrended Q1–Q17 LC systematic error-corrected light curve for KIC 8311864. The stellar parameters adopted for the target were those given in Section 6. The parameters of the simulated stars in the various false positive configurations were drawn from model isochrones from the Dartmouth series (Dotter et al. 2008), and rely also on the known distributions of binary properties such as the orbital periods, mass ratios, and eccentricities (Raghavan et al. 2010; Tokovinin 2014). Other ingredients used for the blend frequency calculation included the number density of stars in the vicinity of the target from the Besançon Galactic structure model (Robin et al. 2003, 2012),24 and the rates of occurrence of eclipsing binaries (Slawson et al. 2011) and of transiting planets from the Kepler Mission itself (NASA Exoplanet Archive, consulted on 2014 March 16).25 The BLENDER constraints for the BEB, BP, and HTP scenarios are illustrated in Figures 10–12. Our simulations indicate that eclipsing binaries orbiting the target (HTS scenario) invariably produce light curves that have the wrong shape for a transit, or feature noticeable secondary eclipses that are not seen in the KIC 8311864 photometry. Similarly, blends involving background giant stars yield light curves that also have the wrong shape. We therefore restricted our analysis to main-sequence stars. Figure 10 displays a sample cross-section of the blend parameter space for the BEB scenario, with the horizontal axis showing the mass of the primary star in the background eclipsing binary, and the vertical axis representing the relative distance between the background binary and the target. Blends with parameters interior to the white contour give a goodness of fit to the Kepler photometry that we considered indistinguishable from that of a transiting planet fit (to within 3 σ). We consider these as viable false positives. Also indicated in the figure are regions in which blends have an color index that does not match that of the target (see figure caption), or where blends are not allowed because the intruding EB would have been detected spectroscopically. Similar diagrams for the BP and HTP configurations are shown in Figures 11 and 12. Further explanations regarding the BLENDER constraints may be found in the recent work by Torres et al. (2015). Figure 10. Map of the surface (goodness of fit) for KIC 8311864 corresponding to blends involving background eclipsing binaries. On the vertical axis we represent the linear distance between the BEB and the target , cast for convenience in terms of the difference in distance modulus, . Only blends within the solid white contour (darker colors) provide fits to the Kepler light curve that are within acceptable limits ( , where σ is the significance level of the difference compared to a transiting planet model fit; see Fressin et al. 2012). Other concentric colored areas (lighter colors) represent fits that are increasingly worse ( , , etc.), which we considered to be ruled out. The blue cross-hatched areas correspond to regions of parameter space where the blends are either too red (left) or too blue (right) compared to the measured color of the target, by more than three times the measurement uncertainty. The dashed green line labeled = 8.4 is tangent to the white contour from above, and corresponds to the faintest viable blends. The green line labeled = 5.0 represents the spectroscopic limit on faint background stars. All simulated blends below this line (green hatched region) are brighter and are generally excluded if the BEB is angularly close enough to the target to fall within the slit of the spectrograph. Download figure: Standard image High-resolution image Export PowerPoint slide Figure 11. Similar to Figure 10 (and with the same color scheme) for blends involving background or foreground stars transited by a planet (BP scenario). The faintest blends giving acceptable fits have = 8.0 relative to the target (dashed green line). Blends below the solid green ( = 5.0, hatched green area) are generally excluded by the spectroscopic constraint. Download figure: Standard image High-resolution image Export PowerPoint slide Figure 12. Similar to Figure 10 for the case of physical companions to KIC 8311864 that are transited by a planet (HTP). Only companion stars with masses larger than about 0.96 or smaller than about 0.56 yield blend colors that are consistent with the measured index of KIC 8311864. However, some of these blends are eliminated by the spectroscopic constraint, as they are brighter than = 5.0 (green hatched area). The only false positives that would remain undetected (viable) have companion stars smaller than 0.50 and are transited by planets with sizes that place them within the white contours. Download figure: Standard image High-resolution image Export PowerPoint slide Based on the above analysis, and including also the information derived from the follow-up observations mentioned earlier, we arrived at estimates of the blend frequencies for the BEB, BP, and HTP scenarios of , , and , respectively. These estimates make use of the frequencies of eclipsing binaries generated by the Kepler team (Slawson et al. 2011) as well as the frequencies of transiting planets of the specific size required to produce the above blends, based on the available sample of KOIs from the NASA Exoplanet Archive. The dominant contribution to the total blend frequency, as is often the case for Kepler candidates, was from false positives involving an unseen companion to the target star that is in a wide orbit and is transited by a larger planet (HTP). On the other hand, the expected rate of occurrence of transiting planets of the size of Kepler-452b, based again on the available sample of KOIs from the NASA Exoplanet Archive, is . The "odds ratio" is therefore . Consistent with previous applications of BLENDER , we adopted a threshold for validation of , corresponding to a 3 σ significance. Any odds ratio above this has an even greater confidence level. This means that Kepler-452b is firmly validated at the 99.76% level.

10. HABITABILITY AND COMPOSITION OF KEPLER-452B

Here we examine the probability that Kepler-452b orbits in the habitable zone of its star, consider the effects of the unknown eccentricity on the habitability, and estimate the likelihood that it is rocky. Of chief interest is whether the insolation flux experienced by Kepler-452b would permit water to exist in liquid form on its surface. Secondarily, we wished to determine the likelihood that it possesses a rocky surface upon which liquid water could pool. 10.1. Habitability Following the statistical approach taken by Torres et al. (2015), we took the insolation flux levels defined in Kopparapu et al. (2013) for recent Venus and for early Mars to be the inner and outer boundaries of the wide or optimistic habitable zone, and the insolation flux levels for the runaway greenhouse and a maximum greenhouse as the inner and outer boundaries for the narrow or conservative habitable zone. This was somewhat more conservative, as Torres et al. (2015) took the inner boundary of the habitable zone to correspond to that for dry desert worlds from Zsom et al. (2013), which permits habitable worlds as close as 0.38 AU from a solar-like host star if the albedo is high and the relative humidity is low (~1%). The 900,000 MCMC chains for the stellar parameters discussed in Section 6 were convolved with realizations of the measured orbital period distribution to construct a population of insolation flux levels relative to Earth's insolation, calculated as where = 5778 K is the solar temperature, is the effective temperature of the star, a is the mean separation between the star and the planet in AU, and R * is the radius of the star. For each realized effective temperature value, the insolation flux values corresponding to the inner boundaries of the optimistic and conservative habitable zones were calculated as per Kopparapu et al. (2013) and the sample was marked as habitable if the insolation was less than that for recent Venus/runaway greenhouse. Figure 13 shows a false color image of the joint distribution between insolation and planetary radius for all 900,000 realizations, where the planetary radius samples were generated by combining the stellar radius samples with a realization for the transit depth, or ratio of the planetary radius to that of the parent star, . Marginalizing the joint distribution over , we find the likelihood that Kepler-452b falls within the wide or optimistic habitable zone is 96.8%, while the likelihood is 28.0% that it falls within the narrow or conservative habitable zone. Figure 13. The joint density of the planetary radius, , of Kepler-452b and its insolation flux, , in Earth insolation units, displayed as a false color image. The vertical lines show the insolation flux values for the modes of the MCMC realizations for the SpecMatch stellar parameters presented in Section 6 as developed by Kopparapu et al. (2013): the optimistic or wide habitable zone is bounded on the inside by recent Venus (solid line) and on the outside by early Mars (dash dotted line), and the conservative or narrow habitable zone is bounded on the inside by the runaway greenhouse effect (dashed line) and on the outside by the maximum greenhouse effect (dotted line). Download figure: Standard image High-resolution image Export PowerPoint slide Given that the SpecMatch results are somewhat conservative with respect to the surface gravity, we repeated the calculations above for the SPC stellar parameters obtained for the HIRES spectrum as a means by which to explore the reasonable non-conservative limits. The SPC parameters yielded realizations that are within the wide and the narrow habitable zone 99.0% and 59.3% of the time, respectively. 10.2. Effect of Eccentricity on Habitability The eccentricity of the orbit of Kepler-452b is poorly constrained by the photometric measurements, at best. With an orbital period of 384.84 days, and a mass of ~5 (see Section 10.3), Kepler-452b should exhibit a reflex velocity of ~45 cm s−1. As the host star is relatively faint for such work ( = 13.43), there is little chance of measuring its orbit through radial velocity observations in the near future. To investigate the effect of the uncertainties of the eccentricity of the orbit on the habitability of Kepler-452b, we examined the variation of the mean insolation flux as a function of eccentricity. Figure 14 shows the mean, maximum and minimum insolation flux for realizations of an orbit at the semimajor axis and orbital period of Kepler-452b as a function of the eccentricity from 0 to 1. While the maximum insolation flux exceeds 4 for , the mean flux does not exceed the maximum permitted by the optimistic habitable zone (recent Venus) until e exceeds 0.8. The probability density for eccentricity developed by Rowe et al. (2014) based on the population of Kepler multiple transiting planets is also shown on this figure, illustrating the fact that transiting planets exhibit very low eccentricities in general. Rowe et al. (2014) modeled the eccentricity distribution as a beta distribution with x = 0.4497 and y = 1.7938. The eccentricity distribution-weighted mean insolation flux for Kepler-452b is 1.17 , only 9% higher than that for e = 0. Thus we conclude that the uncertainties in the eccentricity of Kepler-452b do not significantly affect the habitability of its orbit. Figure 14. The insolation flux, , in Earth insolation units received by the planet Kepler-452b as a function of orbital eccentricity, along with a prior distribution for eccentricity for transiting planets as per Rowe et al. (2014). The mean, max, and min insolation flux values are indicated by the solid blue, the dashed green and the dashed–dotted red curves, respectively. The optimistic and conservative habitable zones, according to Kopparapu et al. (2013), are indicated by the light and dark green shaded regions, respectively. The eccentricity prior is indicated by the solid magenta curve with the square symbols. Download figure: Standard image High-resolution image Export PowerPoint slide 10.3. Composition Because Kepler-452b does not have a mass measurement, estimating the probability that it is rocky necessarily requires the application of statistical mass–radius (M–R) relationships. One such relation is provided by Weiss & Marcy (2014), where the mass of the planet is a deterministic function of its radius. Alternatively, the analysis of Wolfgang et al. (2015) yields a probabilistic M–R relation that allows a range of masses to exist at a given radius, as we would expect if there is a distribution of compositions among the population. Either can be used to produce pairs that can then be compared to the theoretical M–R curves of Fortney et al. (2007): If the radius is smaller at the assigned mass than that for a 100% silicate planet, as determined by the 100% rock mass fraction in the rock-iron composition curve (Equation (8) of that paper), then that pair is consistent with being rocky. We applied both the deterministic M–R relationship of Weiss & Marcy (2014) and the probabilistic formulation of Wolfgang et al. (2015) to investigate the likelihood that Kepler-452b is rocky. Kepler-452b's radius is actually a distribution which incorporates measurement uncertainty reflected in the posterior distribution of the MCMC realizations. To calculate the probability that it is rocky, we drew a set of radii from this distribution, applied the M–R relations to each to obtain mass estimates, and measured the fraction of draws which were consistent with a rocky planet as defined above. Here the radius distribution was determined by both the posterior produced by the lightcurve fitting of Section 7, and the stellar radius posterior of Section 6 and Figure 7; drawing randomly from each of these sets of posterior samples and multiplying them then gave the set of planet radius posterior samples used for this analysis. Because the Weiss & Marcy (2014) M–R relation is a deterministic one, applying it to this problem is straightforward. However, it does not explicitly capture the intrinsic variability in planet masses that we both expect and observe, so the rocky probabilities that it produces may be inaccurate and biased. To improve on this, we computed the posterior predictive distribution for Kepler-452b's mass given the probabilistic relationship of Wolfgang et al. (2015). This involved drawing posterior samples from the M–R relation hyperparameters (i.e., the power law constant, index, and the mass dispersion around this power law) that were produced by evaluating the hierarchical Bayesian model in this sample (see Wolfgang et al. 2015, for more details). This draw defined the M–R relation which then yielded the mass corresponding to the drawn radius. The resulting posterior predictive distribution therefore accounts for both the intrinsic, expected dispersion in planet compositions in this small-radius range and for the uncertainties in the parameter values themselves. As a result, the Wolfgang et al. (2015) approach most accurately incorporates the degree of our uncertainty in this radius-to-mass "conversion." Specifically, we took the deterministic mass–density relation from Equations (1)–(3) in Weiss & Marcy (2014): where and are the planetary mass and radius, respectively, in Earth units. The probabilistic M–R relation of Wolfgang et al. (2015) has the form: where denotes a normal (Gaussian) distribution with mean, μ, and standard deviation, σ, and where and are in Earth units. The radius of a pure rock composition planet was calculated as per Equation (8) in Fortney et al. (2007) for the estimated planetary mass: where is the rock mass fraction (set to 1 for this analysis). We computed the probabilities that Kepler-452b is rocky for both the conservative SpecMatch stellar parameters and for the SPC results. Applying the M–R relation of Weiss & Marcy (2014) gave rocky probabilities of 40% and 64%, respectively. The posterior predictive distribution displayed in Figure 15 obtained using the probabilistic M–R relation of Wolfgang et al. (2015) gave 49% and 62%, respectively. The good agreement between the results for both M–R relationships lends confidence in the results and we estimate the likelihood that Kepler-452b is rocky as between that for the SpecMatch and the SPC parameters, namely, between 49% and 62%. Figure 15. The joint density of the planetary radius, , of Kepler-452b and its mass, , in Earth units, displayed as a scatter plot with marginal histograms for each quantity for the SpecMatch derived parameters. The planetary mass realizations were generated using a hierarchical Bayesian model as per Wolfgang et al. (2015) which explicitly models the planetary mass as a distribution for a given planetary radius. The data and histograms are colored red for those realizations that are denser than predictions for 100% silicate planets as per Fortney et al. (2007), and black for those realizations that are less dense and likely to be non-rocky. The fraction of realizations allows us to estimate the probability to be 49% that Kepler-452b is rocky for this set of stellar and planetary parameters. Download figure: Standard image High-resolution image Export PowerPoint slide We note that it is unlikely that Kepler-452b has an Earth-like composition: we computed the density of Earth-like composition planets for the SpecMatch realizations according to Fortney et al. (2007, i.e., a rock-iron mixing ratio of 2/3 in Equation (8) therein) and found that only 16% and 22% of the Weiss & Marcy (2014) and of the Wolfgang et al. (2015) realizations, respectively, were at least as dense as that predicted for Earth-composition planets. 10.4. History of the Habitability of Kepler-452b Given the similarities between the Sun and KIC 8311864 with respect to mass and effective temperature, and the fact that the orbital distance of Kepler-452b is essentially 1 AU, it is interesting to consider the habitability of this planet over the history of this planetary system. It is important to note that statements comparing the age of KIC 8311864 to that of the Sun are tentative, given that the estimated age for KIC 8311864 is heavily model dependent and the 2 Gyr uncertainty estimated in Section 6 does not take into account variations of poorly constrained model physics such as convection (mixing length parameter) and microscopic diffusion. In this section, we use our best estimate for the stellar age to consider a case study of the evolution of the habitable zone over the main-sequence lifetime for stars such as KIC 8311864 and our own Sun. We traced the evolution of the effective temperature and the luminosity of KIC 8311864 using the Dartmouth stellar evolution models (Dotter et al. 2008)26 appropriate to its mass (1.037 ) and metallicity ([Fe/H] = 0.21). The evolving effective temperature allowed us to determine the inner and outer radii of the conservative and optimistic habitable zone per Kopparapu et al. (2013) as a function of time. The luminosity history obtained from the Dartmouth models allowed us to estimate the radiation received by this super Earth over the main-sequence lifetime of its star. At ~6 Gyr, KIC 8311864 and its planet are about 1.5 Gyr older than the solar system, and the host star's radius is ~10% larger than that of the Sun. This planet bathes in a stream of radiation in its orbit at a level only ~10% higher than that experienced by Earth today. However, since the mass of the host star is only ~4% greater than that of the Sun, Kepler-452b spent the first ~5 Gyr in the conservative habitable zone according to Kopparapu et al. (2013), as illustrated in Figure 16. This planet will remain in the optimistic habitable zone for another ~3.5 Gyr before the star leaves the main-sequence and expands rapidly in the ascending red giant branch phase. Figure 16. A history of the habitable zone of KIC 8311864 as delineated by the amount of flux received by a planet at the orbital distance of Kepler-452b (1.05 AU—essentially the same as that of the Earth) over the main-sequence lifetime of its star. This figure plots the radiation received by Kepler-452b from 0.25 to 11 Gyr as a thick blue curve. The optimistic and conservative habitable zones, according to Kopparapu et al. (2013), are indicated by the light and dark green shaded regions, respectively. The planet is currently experiencing radiation at a level 10% above that received by Earth. However, over the first 5.2 Gyr of its history, Kepler-452b was in the conservative habitable zone, receiving less than 1.0369 of the flux received by contemporary Earth. The red and purple stippled curves represent the flux received by planets in Mars' and Venus' orbits (1.52 and 0.72 AU), respectively. Download figure: Standard image High-resolution image Export PowerPoint slide If this exoplanet system is a future version of an analog to our own solar system, with additional planets in the orbits corresponding to those of Mars and Venus, their insolation history would be very similar to their namesakes in the solar system. An exo-Venus would have spent only the first ~3 Gyr in the optimistic habitable zone before exiting it for good. In contrast, an exo-Mars would be located in the conservative habitable zone during the entire ~11 Gyr spent by the host star on the main-sequence. At present, Kepler-452b is the only known planet orbiting this star. Assuming that the dispersion in orbital inclination angles is the same in the KIC 8311864 system as it is for the solar system (~1 9), Kepler would have had only a ~10% chance of observing a transiting exo-Venus in this system as well as Kepler-452b. The probability of observing a few transits of an exo-Mars in addition to those of Kepler-452b is ~6%, but it would have to have been a rather large planet in order for its transits to be readily identifiable in the light curve by eye since it would have provided fewer than three transits over the course of the Kepler Mission. No such transits have been seen. Kepler had only a ~2% chance of seeing transits of both an exo-Venus and an exo-Mars in addition to those of the super Earth at 1.05 AU. Until such time as other planets are discovered by other exoplanet surveys, or a SETI observation of this planetary system reveals signatures of extraterrestrial technologies, we are left to speculate on the fate of an ancient civilization that may have developed first on Kepler-452b (or on a moon orbiting it). For example, it may have subsequently migrated to an as-yet-undetected outer planet to escape the inevitable loss of most of the intrinsic water inventory after the moist runaway greenhouse effect took hold at Kepler-452b's orbital distance approximately 800 Myr ago.

11. CONCLUSIONS

The identification of a 1.6 planet in the habitable zone of a G2 star at a mean separation of 1.05 AU represents the closest analog to the Earth–Sun system discovered in the Kepler data set to date with respect to the orbit and spectral type of the host star. Moreover, the small radius of this planet provides a reasonable chance between 49% and 62% that it is rocky, although in this case it is unlikely to have an iron core of significant mass. Given our current best estimates the host star is 10% larger and 1.5 Gyr older than the Sun; hence, Kepler-452b may represent the future environment that the Earth will experience as our own Sun evolves toward the red giant phase. Though it is currently just beyond the moist runaway greenhouse level for the incident radiation, this planet is still within the optimistic habitable zone as per Kopparapu et al. (2013), and will remain so for another ~3 Gyr. However, Kepler-452b likely spent the first ~5 Gyr in the conservative habitable zone, below the radiation level experienced by Earth today. The paucity of small planets in or near the habitable zone of solar-like stars, particularly G stars, makes the determination of the intrinsic frequency of Earth-size planets in the habitable zone of Sun-like stars difficult to determine robustly without resorting to extrapolation. Kepler-452b represents an important addition to the list of small planets in or near the habitable zone for its proximity to a G2 star. We are hopeful that this is only the first of many such planetary systems yet to be discovered in the fullness of time.

Kepler was competitively selected as the 10th Discovery mission. Funding for this mission is provided by NASA's Science Mission Directorate. The authors acknowledge the efforts of the Kepler Mission team in obtaining the light curve data and data validation products used in this publication. These data were generated by the Kepler Mission science pipeline through the efforts of the Kepler Science Operations Center and Science Office. The Kepler Mission is led by the project office at NASA Ames Research Center. Ball Aerospace built the Kepler photometer and spacecraft which is operated by the mission operations center at LASP. The Kepler light curves are archived at the Mikulski Archive for Space Telescopes and the Data Validation products are archived at the NASA Exoplanet Science Institute. Some of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership of the California Institute of Technology, the University of California, and NASA. The Keck Observatory was made possible by the generous financial support of the W. M. Keck Foundation. The authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Mauna Kea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain. M. E. acknowledges support by NASA under grant NNX14AB86G issued through the Kepler Participating Scientist Program. D. W. L. acknowledges partial support from NASA's Kepler Mission under Cooperative Agreements NNX13AB58A with the Smithsonian Astrophysical Observatory. G. T. acknowledges support for this work from NASA under grant NNX14AB83G (Kepler Participating Scientist Program). D. H. acknowledges support by the Australian Research Council's Discovery Projects funding scheme (project number DE140101364) and support by the National Aeronautics and Space Administration under Grant NNX14AB92G issued through the Kepler Participating Scientist Program. A. W.ʼs funding was provided by the National Science Foundation Graduate Research Fellowship under Grant No. 0809125, and by the UC Santa Cruz Graduate Division's Eugene Cota-Robles Fellowship. J. M. J. acknowledges 20 years of unwavering support from Renée M. Schell and excellent editorial suggestions for this paper. We dedicate this discovery to the memory of Kepler's late Deputy PI, Dr. David G. Koch, who represented the epitome of grace, intellect and civility, and who would have been deliriously delighted by the discovery of this planet. Facilities: Kepler.

APPENDIX: A COMPUTATIONALLY EFFICIENT STATISTICAL BOOTSTRAP TEST FOR TRANSITING PLANETS