ABSTRACT
The Kepler mission has discovered more than 2500 exoplanet candidates in the first two years of spacecraft data, with approximately 40% of those in candidate multi-planet systems. The high rate of multiplicity combined with the low rate of identified false positives indicates that the multiplanet systems contain very few false positive signals due to other systems not gravitationally bound to the target star. False positives in the multi-planet systems are identified and removed, leaving behind a residual population of candidate multi-planet transiting systems expected to have a false positive rate less than 1%. We present a sample of 340 planetary systems that contain 851 planets that are validated to substantially better than the 99% confidence level; the vast majority of these have not been previously verified as planets. We expect ∼two unidentified false positives making our sample of planet very reliable. We present fundamental planetary properties of our sample based on a comprehensive analysis of Kepler light curves, ground-based spectroscopy, and high-resolution imaging. Since we do not require spectroscopy or high-resolution imaging for validation, some of our derived parameters for a planetary system may be systematically incorrect due to dilution from light due to additional stars in the photometric aperture. Nonetheless, our result nearly doubles the number verified exoplanets.
Export citation and abstract BibTeX RIS
1. INTRODUCTION
Data from the first two years of Kepler spacecraft operations have identified 3670 target stars with periodic or transit-like signatures indicative of transiting planets or eclipsing binary stars. Approximately 50% of these targets have signatures that can be attributed to false positives (FPs), primarily eclipsing binaries (EBs) centered on the target star, a chance alignment of a distant EB within the photometric aperture, or flux bleeds into the photometric aperture. The remaining 2530 systems are composed of primarily exoplanetary systems with an expected FP rate of approximately 10% due to photometric blends (Morton & Johnson 2011; Fressin et al. 2013; Santerne et al. 2013). However, a subset of 457 systems show more than one candidate transiting planet candidate (PC); we refer to these candidates as "multis." FPs should be nearly randomly distributed among Kepler targets, whereas if flat multi-planet systems are common, then many targets should have multiple transiting planets (Lissauer et al. 2011). The large number of multis observed thus implies a high reliability rate, as quantified by (Lissauer et al. 2012, 2013, henceforth Paper I and Paper II, respectively). The small number of FPs found in multis observed by Kepler (Latham et al. 2011, and Section 5.9 of this paper) reinforces our confidence in the high reliability of the PCs remaining in multis (see Appendix C of Paper II for details). After making selections to minimize the odds of blend scenarios, we find 340 systems containing a total of 851 planets that can be validated to better than the 99th percentile, with 768 planets across 306 systems being newly validated. Some of these systems have also been confirmed from radial velocity (RV) detection (Marcy et al. 2013; Gautier et al. 2012), transit timing variations (Ford et al. 2011; Steffen et al. 2012), and planet validation techniques such as BLENDER (Torres et al. 2011; Fressin et al. 2012) and now through multiplicity boost (Paper I). We increase the known number of exoplanets from 94216 to 1710.
With excellent precision and high-duty cycle, Kepler observations of transiting exoplanet systems provide photometric data that can be used to measure fundamental transit properties such as transit duration and observed transit depth. These properties are determined by the relative sizes of the planet, host star and additional flux sources (other stars) within the photometric aperture. Since observations provide nearly full coverage of the planetary orbit, the resulting photometric phase curves enable a useful diagnostic for the identification of astrophysical FPs that can mimic an exoplanet transit signature. In Section 4, we examine the light curves of the Kepler sample and describe the nature of the planetary systems and how they are identified.
As an imaging instrument, Kepler also provides time series measurements of the centroid of the photometric signal. When multiple sources are present in the photometric aperture, the photometric centroid can move in response to flux changes from any of the sources. This property allows Kepler pixel-level data to be used to search for scenarios where a planetary transit-like event is produced by a diluted background eclipsing binary star (Batalha et al. 2010).
Paper II presents a theoretical exploration of the expected and predicted FP rate for transiting multi-planet systems. From the ∼190,000 targets observed by Kepler, there are roughly 2500 transit-like patterns of events in the Kepler Object of Interest (KOI) catalog, split between PCs and EBs. Thus, for a Kepler target chosen at random, it is unlikely that a transit event will be present. One of the most common classes of FP event is caused by a background eclipsing binary (BEB) in the aperture of the target star. The source of a BEB can either be an additional star that is found within or close to the photometric aperture or a bright star within the Kepler field of view that introduces flux in the photometric aperture due to optical ghosting, such as mirror images or electronic interference, such as CCD crosstalk (J. Coughlin et al. 2013, in preparation). The occurrence rate of FPs is largely independent of the target stars and thus it is far more likely that an FP source will produce a single transit-like event as opposed to a photometric light curve containing transit-like signals from multi-planet sources. It is important to distinguish between FPs produced by background stars from those caused by instrumental effects. We identify the latter as period/phase (P/T0) collisions to indicate that the period and epoch solution for the candidate event are not due to a unique event. Rather, the transit signal from one target is seen on another target as well. P/T0 collisions account for most of the identified instrumental FPs. A full description and catalog of P/T0 collisions can be found in Coughlin et al. (2014). Figure 1 of Paper II shows the galactic distribution of Kepler targets, PCs and FPs. The results indicate that only a small fraction of the remaining Kepler PCs are likely to be BEB FPs.
The process of identifying and cataloging FPs from the KOI list continues to evolve. Thus, the FP list from various iterations of the KOI catalog (Borucki et al. 2011; Batalha et al. 2011, 2013) have different reliability rates. This paper describes in detail the steps taken to develop a reliable and uniform classification scheme and its application to the multi-planet sample to classify KOIs as false alarms (FAs), FPs, or PCs. FAs are transit event candidates with a signal-to-noise ratio (S/N) below a S/N threshold of 7.1 (see Section 5.2) or a transit candidate mimicked by stellar variability or an instrumental artifact. We find that the FP rate for multi-planet systems is low and consistent with the predictions from Appendix C of Paper II based on a statistical analysis of the FP rate found in the single planet population. The predicted FP rate allows us to conclude that the PCs that pass our FP and morphological tests would misclassify only ∼2 PCs, allowing us to claim that 768 PCs are bona fide planets with a confidence level greater than 99%.
As an example, consider only the non-transiting, transiting single planet and transiting double planet systems. Here we ignore the case of FP+FP, where two FPs are associated with the same target and estimate the number of FP for the double planet systems in the spirit of Paper II. There are approximately 140,000 stars that contribute to the transiting planet population (Section 4 of Paper II). After removing FAs and T/P0 collisions there are 2182 systems that show one transiting body and 284 systems that show two transiting bodies. From the single-planet systems, 662 systems were identified as FPs, which provides an estimate of the FP rate of 0.44%. Thus, from the 1500 good single-planet candidates, 7 PC+FP systems are predicted. From the sample of 295 that have two transiting candidates, there are 6 that were identified as a PC+FP combination. Good agreement is found even in our simplified case. Of course, one needs to properly account for EBs and the entire range of multi-planet candidates and the multitude of PC and FP mixes that can be produced. These considerations are the basis of Papers I and II, which provide predictions of FPs rates that are verified in this paper.
The combined analysis reported in Paper II and this manuscript validates more than 300 new Kepler multi-planet systems. Paper II introduces the binary star planet hosts Kepler-132, where one star hosts two transiting planets and its companion hosts one transiting planet and Kepler-296, a pair of small stars with a total of five transiting planets, the multi-resonant four-planet Kepler-223 system and two additional planets in the Kepler-80 = KOI-500 system that includes two three-body resonances, as well as several high-multiplicity systems, including the new five-planet systems Kepler-102 = KOI-82, Kepler-169 = KOI-505, Kepler-238 = KOI-834, and Kepler-292 = KOI-1364, three new planets orbiting Kepler-84 = KOI-1589 (bringing the total count to 5) and partial validations of the five-candidate systems Kepler-122 = KOI-232 (four planets validated) and Kepler-154 = KOI-435 (two planets validated). Hundreds of new planetary systems are announced herein, with special attention given to four new planets with radii roughly twice that of Earth located in or near the habitable zones of their host stars.
This paper is organized as follows: in Section 2, a description of the planetary sample is presented. The adopted stellar parameters are discussed in Section 3. Detailed descriptions of the transit models and lightcurve analysis are described in Section 4. The process of identifying FPs can be found in Section 5. Section 5.6 covers centroid measurements that provide a clean sample of highly probable transiting multi-planet systems which we demonstrate in Section 5.9 are genuine extrasolar planets. We conclude with a discussion of the multi-planet population in Section 6.
2. PLANET CANDIDATE SAMPLE
Photometric surveys for extrasolar planets are contaminated by FPs that are caused by eclipsing stellar binaries and transits and eclipses of stars that are spatially offset from the target stars. The KOI catalog is an inhomogeneous working list used to track transit candidates of interest identified from Kepler photometric light curves. The FP to PC ratio of the raw KOI catalog is approximately 0.39 (Burke et al. 2014).17 A quick survey of KOI dispositions (available on the Kepler exoplanet archive; Akeson et al. 2013) shows most of the FP occurrences are linked to events that show a single pattern of periodic transits. They rarely occur when there is evidence of multiple transiting objects, a fact first noted by Latham et al. (2011). While the vast majority of the catalog is dedicated to exoplanet candidates and eclipsing binary stars, the list contains some astrophysically interesting light curves that do not show transits, such as Heartbeat Stars (Thompson et al. 2012). We have excluded events that have been classified as a non-transit event. Similarly, FAs were also excluded based on a transit S/N of 7.1. The Kepler transit detection pipeline was designed to identify events that have at least three transits. However, through human examination of the photometric time series, some deep, single occurrence events are noted and cataloged. The orbital periods for some of these events have now been measured based on the detection of a third transit. However, there are 18 PCs for which only 1 or 2 transits have been observed. It is difficult to assess the FP nature of these candidates as one cannot reliably estimate the orbital period. These KOIs are also excluded from our sample.
Given that the Kepler pipeline is continually undergoing substantial improvements, we restrict ourselves to exoplanet candidates found with Q1–Q8 light curves. This includes the Q1–Q5 (Borucki et al. 2011) and Q1–Q6 (Batalha et al. 2013) lists and a subset of the Q1–Q8 catalog (Burke et al. 2014). The Q1–Q8 candidate list includes data products based on Q1–Q10 data. During this process additional multi-planet candidates were discovered and cataloged, including PCs found with independent methods such as QATS (Carter & Agol 2013). An example is KOI-351 (Kepler-90), which had six PCs found in the Q1–Q8 sample but is now known to contain at least seven candidates (Paper II and E. Agol et al. 2013, in preparation). We have excluded these Q1–Q10 discoveries to avoid a bias in our statistics that would overestimate the quantity of multi-planet systems for validation. Any multi-planet candidate found in Burke et al. (2014) that is not listed in our sample is a Q1–Q10 discovery. We used 3737 KOIs associated with 3008 stars in our analysis; 1210 PCs in multi-planet systems and 2527 PCs in single planet systems.
Table 1 lists the number of systems and planets that have been considered in our sample. Each row lists the number of systems that pass various tests, such as our FA test (Section 5.2), P/T0 collisions, and transit characterization from models (Section 4.1). The columns indicate the number of PCs found in each system with S1 indicating that one PC was found. M1 means that a light curve was identified to have multiple transiting candidates, but after cuts the number of remaining candidates has fallen to one. The columns labeled 2 through 6 indicate the number of systems with the corresponding number of planets that pass criteria indicated in the last column. Each test was chosen to eliminate and mitigate FPs (see Section 3 of Paper II), but only candidates labeled as FP are definitive FPs. For example, we do not validate candidates that show "V"-shaped light curves based on transit model fits. This distinctive transit shape can be produced by either a transit duration that is similar to the photometric cadence, a grazing planet with an impact parameter near 1, or stellar binary with a large value of Rp/R⋆. Transit-modeling of a "V"-shaped transit leads to a large uncertainty on the measurement of quantities such as Rp/R⋆, making it difficult to assess the properties of the PC and identify blended BEBs. These candidates are not labeled as FPs, but they are also not validated as planets.
Table 1. Multi-planet Counts
S1 | M1 | 2 | 3 | 4 | 5 | 6 | # Multi-pl | # New Multi-pl | Cuts |
---|---|---|---|---|---|---|---|---|---|
2527 | 0 | 313 | 107 | 45 | 13 | 3 | 1210 | 1092 | No Cuts |
2482 | 17 | 300 | 104 | 43 | 13 | 3 | 1167 | 1049 | FA |
2161 | 21 | 295 | 102 | 43 | 13 | 3 | 1151 | 1033 | FA, Col |
1861 | 27 | 293 | 100 | 42 | 10 | 3 | 1122 | 1011 | FA, Col, P |
1496 | 24 | 284 | 102 | 43 | 13 | 3 | 1129 | 1011 | FA, Col, FPs |
1377 | 29 | 285 | 100 | 42 | 10 | 3 | 1106 | 995 | FA, Col, FPs, P |
1492 | 32 | 279 | 100 | 43 | 13 | 2 | 1107 | 989 | FA, Col, FPs, SN |
1366 | 41 | 273 | 98 | 40 | 13 | 2 | 1077 | 964 | FA, Col, FPs, SN, b |
1374 | 37 | 281 | 96 | 43 | 10 | 2 | 1084 | 973 | FA, Col, FPs, SN, P |
1257 | 46 | 272 | 96 | 40 | 10 | 2 | 1054 | 947 | FA, Col, FPs, SN, P, b |
⋅⋅⋅ | 66 | 236 | 81 | 32 | 12 | 1 | 909 | 821 | FA, centr |
⋅⋅⋅ | 66 | 236 | 81 | 32 | 12 | 1 | 909 | 821 | FA, Col, centr |
⋅⋅⋅ | 67 | 235 | 81 | 32 | 12 | 1 | 907 | 819 | FA, Col, FP, centr |
⋅⋅⋅ | 67 | 235 | 81 | 33 | 11 | 1 | 906 | 818 | FA, Col, FP, centr, SN |
⋅⋅⋅ | 72 | 233 | 79 | 34 | 9 | 1 | 890 | 803 | FA, Col, FP, centr, SN, P |
⋅⋅⋅ | 79 | 227 | 78 | 32 | 9 | 1 | 867 | 784 | FA, Col, FP, centr, SN, P, b |
⋅⋅⋅ | 83 | 221 | 78 | 31 | 9 | 1 | 851 | 768 | Dynamical and SP Cuts |
Notes. Number of multi-planet systems after various cuts. See Section 5.9 for detailed definition of the cuts.
Download table as: ASCIITypeset image
Approximately one-third of the entries in the KOI catalog have been identified as EBs either through light curve or centroid analysis. Criteria for KOI promotion have not been uniformly applied, and there have been some systematic biases in categorization based upon multiplicity of sets of transit signatures. For example, a low amplitude sequence of dips with alternating depths is generally categorized as an eclipsing binary and not given a KOI number, but if such a signature is seen for a target that has already has been classified as a KOI based upon another set of dips, it is given a KOI number and immediately classified as an FP, thereby providing a biased increase in the number of FPs in multis. Searching for additional candidates is terminated around targets that have been identified as FPs, although there is often a significant amount of time between initial identification as a KOI and labeling as a FP. If a target is first identified as an EB, then the pipeline transit search is terminated, and since most such identifications are done rapidly, the distribution is more strongly biased against combinations of EBs and PCs. To improve our estimate of the true FP rate, we supplement our FP list with the eclipsing binary catalog.18 Since we are interested in EBs that roughly match the signal from a transiting extrasolar planet we only considered EBs that are detached based on having a morphological classification criteria less than 0.5 (Matijevic et al. 2012) and have a primary eclipse depth less than 2%. These criteria remove contact binaries, which have a characteristic w shape. This step increases the total number of FPs by 138.
3. CHARACTERISTICS OF PLANET-HOSTING STARS
We used a diverse set of measurements to estimate the properties of each stellar host of the Kepler multis that we validate as planets. Our goal is to obtain the best classification of each planet-hosting star given all of the information available to us rather than to produce a homogeneous data set. We prioritized classification in the following order, choosing for each target the first available option.
- 1.Combined asteroseismology + spectroscopy analysis (Huber et al. 2013).
- 2.Spectrometry Made Easy (SME; Valenti et al. 1996) analysis using spectra taken at the Keck I telescope.
- 3.SpecMatch fitting (see below and Petigura et al. 2013) using spectra taken at the Keck I telescope.
- 4.Stellar Parameter Classification (SPC) analysis of spectra taken at various telescopes (Buchhave et al. 2012).
- 5.Modified KIC photometric classification from adjustments to original KIC values of Teff, log g and [Fe/H] to match Yale-Yonsei (Demarque et al. 2004) stellar evolution models.
Note that by using this heterogeneous set of stellar characterization techniques, we sacrifice uniformity for accuracy, and care should be taken in performing statistical studies based on fit parameters.
To support the SME and SpecMatch analyses, high resolution spectra were taken of multi-planet candidate host stars with HIRES spectrometer on the Keck I telescope using the observing setup of the CPS group (Marcy et al. 2008). We acquired spectra with a resolution of R = 55,000 and a wavelength coverage of 360–800 nm, which have a S/N per pixel of 40 (or better in some cases that were used for the SME analysis) at 550 nm, corresponding to a S/N = 85 per resolution element. The spectra were observed without the iodine cell in the light path. Using the C2 decker, which projects to 087 × 14
0 on the sky, we removed the signal from moonlight that otherwise could contaminate the stellar spectra at the level of a few percent.
When determining atmospheric parameters of the planet host stars using SpecMatch, we compared each spectrum to a library of 800 spectra having Teff = 3500–7500 K and log g = 2.0–5.0, which spans the FGK and early M type main sequence and subgiant stars. All library stars have accurate parallax measurements, allowing for good estimates of stellar mass and radius for each. We then compared the observed spectrum with that of each library star. The spectrum is placed on a common wavelength scale and normalized in intensity. The χ2 value is calculated as the sum of the squares of the differences between the observed spectrum and each library spectrum. The final stellar properties, listed in Table 2, are determined by the weighted mean of the ten library spectra with the lowest χ2 values. We adopted errors in each parameter by comparing results to a range of standard stars.
Table 2. Stellar Parameters
KOI | Kepler-ID | KID | Teff | Teffσ | log g | log gσ | [Fe/H] | [Fe/H]σ | R⋆ | R⋆σ | ρ⋆ | ρ⋆σ | Flag | Blend |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
(K) | (K) | (cgs) | (cgs) | (R☉) | (R☉) | (g cm−3) | (g cm−3) | |||||||
41 | Kepler-100 | 6521045 | 5825 | 75 | 4.125 | 0.045 | 0.02 | 0.10 | 1.490 | 0.035 | 0.457 | 0.013 | 5 | 3 |
46 | Kepler-101 | 10905239 | 5570 | 134 | 4.065 | 0.240 | 0.30 | 0.10 | 1.666 | 0.415 | 0.351 | 0.300 | 3 | 4 |
70 | Kepler-20 | 6850504 | 5443 | 74 | 4.398 | 0.100 | 0.00 | 0.07 | 0.986 | 0.095 | 1.304 | 0.400 | 4 | 3 |
72 | Kepler-10 | 11904151 | 5627 | 44 | 4.342 | 0.046 | −0.15 | 0.04 | 1.056 | 0.021 | 1.068 | 0.008 | 5 | 3 |
82 | Kepler-102 | 10187017 | 4908 | 74 | 4.640 | 0.100 | 0.08 | 0.07 | 0.716 | 0.032 | 3.132 | 0.304 | 4 | 3 |
85 | Kepler-65 | 5866724 | 6169 | 50 | 4.236 | 0.035 | 0.09 | 0.08 | 1.424 | 0.024 | 0.621 | 0.011 | 5 | 3 |
89 | 8056665 | 6688 | 342 | 4.059 | 0.150 | −0.21 | 0.10 | 1.773 | 0.357 | 0.329 | 0.186 | 3 | 3 | |
94 | Kepler-89 | 6462863 | 6184 | 83 | 4.196 | 0.068 | 0.11 | 0.07 | 1.486 | 0.139 | 0.543 | 0.127 | 4 | 3 |
102 | 8456679 | 5705 | 100 | 4.311 | 0.150 | 0.18 | 0.10 | 1.199 | 0.219 | 0.867 | 0.419 | 3 | 3 | |
108 | Kepler-103 | 4914423 | 5845 | 88 | 4.162 | 0.051 | 0.07 | 0.11 | 1.436 | 0.039 | 0.513 | 0.020 | 5 | 3 |
Notes. Flag: −1—solar parameters assumed; 0—original KIC; 1—revised KIC; 2—SPC; 3—SpecMatch; 4—SME; 5—asteroseismology blend flag: 0—nearby star detected that may produce blend, 1—no measurement, 2—has speckle, 3—has spectral matching, 4—has both.
Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.
Download table as: DataTypeset image
Stellar parameters are derived by matching atmospheric parameters (Teff, log g, and [Fe/H]) to stellar evolution models (M⋆, Age, and Z). Atmospheric parameters are based on SME (Valenti et al. 1996), SpecMatch, SPC (Buchhave et al. 2012), asteroseismology (Huber et al. 2013), or the KIC (Brown et al. 2011) including the revision of Teff by Pinsonneault et al. (2012). For SME parameters, we added 59 K to Teff and 0.062 dex to [Fe/H] in quadrature following the Torres et al. (2012) recommendation. For SpecMatch and SPC results, we adopt uncertainties as reported, with preference given to SME and then SpecMatch. For stars without asteroseismic or spectroscopic constraints, we adopt Teff from Table 7 of Pinsonneault et al. (2012) and log g and [Fe/H] as given in the KIC. For uncertainties, we adopt values of 200 K in Teff, 0.3 dex in log g and 0.5 dex in [Fe/H], in agreement with typical residuals of KIC values to stellar properties determined from asteroseismology and spectroscopy (e.g., Bruntt et al. 2012). We adopted the Yonsei-Yale stellar evolution models (Demarque et al. 2004) to determine stellar parameters. The model matching was done by varying the stellar mass, age, and Z and comparing the model-derived values of Teff, log g, and [Fe/H] with the spectroscopic values with a chi-square statistic. An initial match was found by scanning in mass increments of 0.1 M☉ restricting ages from 0 to 14 Gyr, and identifying a best matching model. A Markov Chain Monte Carlo (MCMC) routine was then seeded with this trial value of stellar mass, age, and Z to determine posterior distributions. All stellar models with ages greater than 14 Gyr were excluded. In total 100,000 chain elements were generated for each star. The models were also used to determine posterior distributions for the stellar radius, luminosity, and mean stellar density. The resulting stellar parameters are listed in Table 2.
4. LIGHT CURVE ANALYSIS
Kepler photometry was used to both identify FPs and to characterize the transiting planets. We used Kepler Q1 to Q10, long cadence, simple aperture, photometric observations19 gathered every 29.4 minutes over a time span of 868 days. These measurements do not account for the effects of dilution from the addition of stars near or in the photometric aperture, thus, there is a bias in our measured planetary parameters toward underestimating planetary radii. Measuring dilution and determining corrections are difficult tasks and outside the scope of our goal to validate hundreds of extrasolar planets. However, there are rough estimates of the dilution based on the KIC, retrieved from MAST, from which this bias can be estimated and used to place conservative upper limits on systematics introduced from contamination. The mean value of light contamination for validated Kepler planets is 5%. As the transit depth is proportional to (Rp/R⋆)2, a 5% dilution translates into a 2.6% systematic bias in the planetary radius, which is small compared to the uncertainty in the stellar radius. From our validated planet sample, the largest light contamination was found to be 20% of the total light for KOI-907 (Kepler-251), which translates into an error on the planetary radii of 11.8%.
Although Kepler had a high duty cycle, some transits were missed due to a variety of logistical details such as sky location, data downlink, spacecraft safe modes, and a dead module. An extreme example is KOI-94 (Kepler-89) (Weiss et al. 2013), which has an effective duty cycle of less than 50% due to its location on the edge of one of the CCD detectors. The spacecraft rotates each quarter and this target lies in the gap between detectors during two of every four consecutive quarters. For a majority of our targets, the effective duty fraction after removal of flagged20 data was approximately 86%.
We filtered the data to remove instrumental and astrophysical signatures that are independent of the planetary transit as follows: each observation was corrected by fitting a cubic polynomial to a segment of the time series photometric measurements centered on the time of measurement. A segment is defined by selecting observations that were taken within one day of the measurement. We also require that the time series not contain any gaps longer than five cadences (∼2.5 hr). If such a gap is encountered, the data collected near that gap are not considered. Such gaps were commonly produced by the monthly data downlinks. The removed data dropped the duty cycle by ∼1%. After repointing the spacecraft, there was usually a photometric offset produced due to thermal changes in the telescope. Thus, astrophysical signals with timescales of approximately two days are strongly filtered by this process. The filter is destructive to the shape of a planetary transit. Thus, we exclude any measurement taken within one transit duration of the measured center of the transit time and use an extrapolation of the polynomial fit to estimate corrections during transits. The transit duration is defined as the time from first to last contact, Tdur. The segment is fit with a cubic polynomial and used to measure the photometric offset, which is then removed. We repeat the process for each observation to produce a detrended time series. When significant transit timing variations are detected, we rebuild the detrended time series using the updated center of transit times.
An initial multi-planet photometric model was fit to each detrended light curve. The photometric model assumes non-interacting circular orbits and used the quadratic limb darkening transit model of Mandel & Agol (2002). We used limb-darkening parameters from Claret & Bloemen (2011), which were fixed for each target based on our stellar classification (Teff, log g, and [Fe/H]). The model was parameterized by the mean-stellar density (ρ⋆), photometric zero point, and for each planet (n) an epoch (T0n), period (Pn), scaled planetary radius (Rp/R⋆n), and impact parameter (bn). The semi-major axis for each PC is estimated by
![Equation (1)](https://content.cld.iop.org/journals/0004-637X/784/1/45/revision1/apj490915fd1.gif)
where the assumption was made that the sum of the planetary masses is much less than the mass of the host star. For a Jupiter-mass companion of a Sun-like star, a systematic error of 0.1% is incurred on the determination of ρ⋆. To account for the ∼30 minute integration of Kepler observations, the transit model was sampled 11 times temporally with equal spacings within the integration window. The 11 separate models were then averaged. A best fit model was calculated by a Levenberg–Marquardt chi-square minimization routine (More et al. 1980). This model was used primarily to seed our MCMC routines to measure fundamental physical properties of each planet.
4.1. Measuring Planet Parameters
Our main objective is to identify FPs and to select candidates found in multi-transiting systems that have a very high probability of being bona fide extrasolar planets. Our strategy was to examine each photometric light curve for signatures of stellar binarity: secondary eclipses, phase curve variations and a comparison of the transit model determination of ρ⋆ to our classification and modeling of the host star. We also examine the populations of stellar and planetary systems to establish regions of parameter space, namely orbital period and impact parameter, that are most susceptible to contamination from FPs.
Our measured planetary parameters are listed in Table 3 and are based on a transit model fit similar to the description given at the beginning of Section 4, except that we have modeled each PC in a system independently. We start by using the best fit model from the multi-planet model to remove the photometric signature of all transiting candidates except the one we wish to measure. We assumed a circular orbit and fit for ρ⋆, T0, P, b, Rp/R⋆, and ρc, where ρc is the value of ρ⋆ when a circular orbit is assumed. Thus, each PC provides an independent measurement of ρc. If the value of ρc is statistically the same for each PC, then the planetary system is consistent with each planet being in a circular orbit around the same host star. We examine the distribution of transit-determined values of ρc in Appendix A.
Table 3. Transit Model Parameters
KOI | Kepler ID | P | Rp | S | b | Rp/R⋆ | ρc | Tdep | Tdur | T0 | S/N | fp | DynTest |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
(days) | (R⊕) | (S⊕) | (g cm−3) | (ppm) | (hr) | BJD-2454900 | |||||||
41.01 | Kepler-100 c | 12.815842 | 2.25 | 189.1 | 0.38 | 0.01384 | 0.3520 | 221.8 | 6.344 | 55.94755 | 123.4 | 5 | 1 |
0.000029 | 0.11 | 18.6 | +0.01–0.38 | 0.00059 | 0.0678 | 2.5 | 0.042 | 0.00135 | |||||
41.02 | Kepler-100 b | 6.887067 | 1.47 | 429.3 | 0.87 | 0.00901 | 0.1008 | 72.6 | 4.352 | 66.17732 | 48.2 | 5 | 1 |
0.000036 | 0.12 | 42.3 | +0.00–0.60 | 0.00073 | 0.2565 | 2.3 | 0.080 | 0.00305 | |||||
41.03 | Kepler-100 d | 35.333072 | 1.55 | 49.0 | 0.59 | 0.00952 | 0.8392 | 100.8 | 5.826 | 86.98063 | 33.8 | 5 | 1 |
0.000319 | 0.06 | 4.8 | +0.34–0.36 | 0.00027 | 0.6508 | 4.1 | 0.151 | 0.00562 | |||||
46.01 | Kepler-101 b | 3.487691 | 5.87 | 1037.8 | 0.02 | 0.03226 | 0.5797 | 1263.5 | 3.834 | 103.93126 | 360.2 | 7 | 1 |
0.000002 | 1.46 | 625.0 | +0.23–0.02 | 0.00018 | 0.0492 | 5.1 | 0.013 | 0.00031 | |||||
46.02 | Kepler-101 c | 6.029809 | 1.33 | 517.1 | 0.66 | 0.00734 | 0.3429 | 57.8 | 4.030 | 65.48205 | 11.3 | 6 | 1 |
0.000115 | 0.34 | 311.4 | +0.05–0.49 | 0.00047 | 0.3619 | 6.7 | 0.339 | 0.01146 |
Notes. fp: −1—false-alarm, 0—false positive or false-alarm, 1—period/epoch collision, 2—not clean centroid, 3—nearby stars makes unclean, 4—unsaturated manual centroid pass, 5—saturated pass, 6—Q1-Q12 autopass, 7—Q1-Q15 autopass. DynTest: 0—failed dynamic test, 1—dynamic test not performed, 2—passes dynamic test. Periods marked with "tt" indicated that transit timing variations were accounted for in the transit models. Radii marked with * indicate systems known to be planet hosting stellar binaries. As dilution is not accounted for, the actual planet radii are larger.
Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.
Download table as: DataTypeset image
To estimate the posterior distribution on each fitted parameter, we use a MCMC approach similar to the procedure outlined in Ford (2005). To account for the strong correlation between ρc, b, and Rp/R⋆, we use a Gibbs sampler to shuffle the value of parameters for each step of the MCMC procedure and use a control set of parameters to approximate the scale and orientation for the jumping distribution of correlated parameters as outlined in Gregory (2011). An initial control set consists of 2000 chains was generated by a MCMC run where the width of the Gaussian proposal distributions was adjusted to achieve a success rate of ∼25%. Once the success rate for a jump was between 20% and 30%, the width of the Gaussian was fixed for the duration of the calculations. The control set is updated during the MCMC run by adding every second accepted jump proposal parameter set and removing the oldest element. The control set is fixed once an acceptance rate between 20% and 25% is achieved. Any chain that was generated before the proposal sample was fixed was discarded. We found that this method allows the MCMC approach to efficiently sample parameter space even with highly correlated model parameters. We generated four 1,000,000 Markov chains for each PC. The first 20% of each chain was discarded and the remaining sets were combined and used to calculate the median; standard deviation; and 1σ, 2σ, 3σ bounds of the distribution centered on the median of each modeled parameter. Our model fits and uncertainties are reported in Table 3. We use the Markov chains to derive model dependent measurements of the transit depth (Tdep) and transit duration (Tdur). We also convolve the transit model parameters with the stellar parameters (see Section 3) to compute the planetary radius, Rp, and the flux received by the planet relative to the Earth (S). To compute the transit duration, we used Equation (3) from Seager & Mallen-Ornelas (2003) for a circular orbit,
![Equation (2)](https://content.cld.iop.org/journals/0004-637X/784/1/45/revision1/apj490915fd2.gif)
which defines the transit duration as the time from first to last contact. We estimate the ratio of incident flux received by the planet relative to the Earth's incident flux,
![Equation (3)](https://content.cld.iop.org/journals/0004-637X/784/1/45/revision1/apj490915fd3.gif)
where Teff is the effective temperature of host star, Teff☉ is the temperature of the Sun and a is the semi-major axis calculated with Kepler's Second Law using the measured orbital period and estimated stellar mass.
4.2. Transit Timing Variations
We estimate transit timing variations (TTVs) for each light curve using the best fit models from Sections 4 and 4.1 as a template. Center-of-transit times are measured by selecting data obtained within one transit duration of the predicted center of transit time (thus the time series has a length that is twice the transit duration). If the transit duration is less than two hours, then we select data within two hours of the center-of-transit time. We then refit the transit model but we allow only T0 to vary. The measured center-of-transit time is then compared to the predicted time to produce the observed minus calculated transit time (TTVn) for each transit, n. If significant variations are detected, we improve the transit model template by compressing and expanding the time interval between measurements by linearly interpolating between timing offsets observed for each transit. The improved template is then used to redetermine the transit times. We report our measured transit-timing variations in Table 4. If fewer than four observations were selected for fitting, we do not report TTVn. Note that here, tn is the measured transit time and not the prediction of a linear transit ephemeris (unlike Ford et al. 2011).
Table 4. TTV Measurements
n | tn | TTVn | TTVnσ |
---|---|---|---|
(days) | (days) | ||
KOI-89.01 | |||
1 | 83.56887677 | −0.00469193 | 0.00463798 |
2 | 168.25615788 | 0.00778280 | 0.00516824 |
3 | 252.94343900 | −0.00271064 | 0.00555464 |
4 | 337.63072011 | −0.00016954 | 0.00476691 |
5 | 422.31800122 | 0.00327206 | 0.00417621 |
6 | 507.00528233 | −0.00463520 | 0.00479403 |
7 | 591.69256345 | 0.01299045 | 0.00805245 |
8 | 676.37984456 | 0.00975161 | 0.00425755 |
9 | 761.06712567 | 0.02165216 | 0.00395487 |
10 | 845.75440678 | −0.04881705 | 0.01201981 |
11 | 930.44168789 | 0.00158643 | 0.00540031 |
12 | 1015.12896901 | 0.00220596 | 0.00437501 |
13 | 1099.81625012 | 0.00631963 | 0.00419635 |
KOI-89.02 | |||
1 | 222.88301163 | −0.02540113 | 0.00275003 |
2 | 430.46918841 | 0.02681994 | 0.00248308 |
3 | 638.05536519 | 0.00738233 | 0.00321930 |
4 | 845.64154197 | 0.01459505 | 0.00787256 |
5 | 1053.22771875 | −0.02928310 | 0.00312264 |
KOI-94.01 | |||
1 | 65.74076388 | 0.00225177 | 0.00036658 |
2 | 88.08374268 | −0.00049457 | 0.00028418 |
11 | 289.17055180 | −0.00025281 | 0.00049360 |
12 | 311.51353060 | 0.00114436 | 0.00270664 |
13 | 333.85650939 | 0.00046567 | 0.00432967 |
14 | 356.19948818 | 0.00106531 | 0.00057328 |
15 | 378.54246697 | 0.00003749 | 0.00055978 |
16 | 400.88544576 | −0.00020364 | 0.00056129 |
17 | 423.22842456 | 0.00156312 | 0.00027064 |
18 | 445.57140335 | −0.00175704 | 0.00029028 |
19 | 467.91438214 | 0.00051162 | 0.00059437 |
28 | 669.00119127 | 0.00237404 | 0.00036207 |
29 | 691.34417006 | −0.00064501 | 0.00285220 |
30 | 713.68714885 | 0.00142408 | 0.00038260 |
32 | 758.37310644 | −0.00142232 | 0.00126874 |
33 | 780.71608523 | 0.00043150 | 0.00025590 |
34 | 803.05906402 | −0.00101907 | 0.00030884 |
35 | 825.40204281 | 0.00169851 | 0.00036495 |
45 | 1048.83183073 | −0.00113621 | 0.00041459 |
... |
Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.
Download table as: DataTypeset image
Although the vast majority of the TTVs were processed in bulk, some KOIs with large TTVs received individual attention. When the center-of-transit time was shifted substantially away from the predicted transit time, the fitting process failed. An example is KOI-142 (Kepler-88), where the transit times shift by ∼20 hr. For such cases, the previous two transit timing measurements were used to linear extrapolate an estimate of the next transit time to initialize the fitter.
5. PLANET DISPOSITIONS AND FALSE POSITIVE IDENTIFICATION
The adopted Q1–Q8 dispositions were produced by a combination of work developed for the general KOI catalog (Burke et al. 2014) and the multi-planet population listed in this paper. The end result is a set of dispositions shared between the two papers. Each PC was subjected to tests described below: lightcurve inspection by eye, a S/N threshold, searches for secondary events, phase-linked variations, odd–even numbered transit comparison and centroid motion during transit. The disposition of KOIs presented here and the underlying statistics presented in Paper II sets the stage for the validation of a large number of multi-planet candidates at greater than the 99th percentile; specifically, we expect ∼2 FPs from the 851 planets validated in Section 5.9.
One of the major results from Paper II (see Section 4 therein) is that the FP rate in multi-planet systems must be low. The predictions are that ∼27 FPs should have been detected in the multi-planet candidate sample, and that ∼2 FPs have been missed. Demonstration of the accuracy of the first of these predictions builds a strong case that currently viable PCs in multi-planet systems are bona fide planets. The types of FPs that can be searched for include: a planet transiting a star not physically bound to the Kepler target star; or an eclipsing binary star system or other astrophysical phenomenon. If a bound stellar companion is found, it is sometimes possible to determine which star is the source of the transits, but Kepler data are not sensitive to isolating the transit host in most bound systems.
This section is dedicated to describing the tests that were carried out to identify FPs in our multi-planet sample. These tests include searches for secondary events and classifying the event as a planetary occultation or stellar eclipse (Section 5.3). Tidal interactions and motion of the host star around the center-of-mass produce variations related to the orbital period. When present, the amplitude of these variations reveals the masses of the binary components (Section 5.4). In the cases where the primary and secondary eclipses of a stellar binary with nearly equal mass stars are reported, the orbital period may be incorrect by a factor of two. To test for this scenario we compare the depths of the odd- and even-numbered transit events (Section 5.5). One of the most powerful tests is the use of pixel-level Kepler data to focus on finding clean targets for validation by localizing the source of the transit on the detector (Section 5.6). A common source of FPs are centroid offsets due to motion in the difference of in- and out-of-transit combined images across a transit event. The most frequent source of centroid offsets are BEBs that track the spatial density of background stars (see Figure 1 of Paper II).
5.1. Quality of Model Fit
We calculated the reduced chi-square for each transit model. If the value was greater than 2 or less than 0.5, then the fit and photometry were visually inspected. In most cases, it was found that a transit overlapped an instrumental effect, the most common effect being photometric deviations observed after the instrument returned to nominal operations that involved a reorientation of the spacecraft. In these cases, the offending segment of data was excluded and the model fits were repeated. Other cases include models that produced a poor transit fit from convergence to a local minimum, excess scatter from stellar variability, and evidence of a stellar binary in the light curve shape from the presence of a strong occultation or ellipsoidal variations. This level of vetting was performed for both the single and multi-planet population.
From our inspection we discovered that KOI-1134.01 and 1134.02 were both tracing the same EB that had a period of 100 days but had transit depths that were heavily modulated due to third-light contamination that appeared as independent transit candidates in early analysis using only a few months of observations. We have labeled KOI-1134.01 as an EB FP and 1134.02 was labeled as an FA. KOI-1792.02 was found to be an FP with stellar variability mimicking a transit signal. KOI-2048.02 was identified as residuals of the transit fit to KOI-2048.01, thus 2048.02 was labeled as an FA and removed.
5.2. Signal-to-Noise Ratio
The S/N was calculated from depth of the transit using the transit model and noise was estimated by the standard deviation of observations obtained outside of transit and then scaled with a geometric sum to match the transit duration, yielding
![Equation (4)](https://content.cld.iop.org/journals/0004-637X/784/1/45/revision1/apj490915fd4.gif)
where Tdep is the transit depth, NT is the number observations obtained during transit and σOT is the standard deviation of out-of-transit observations. The estimation of the S/N assumes that the depth of the transit is uniform, which is a good approximation for small Earth-sized planets with central transits (b = 0). For relatively large planet-to-star radius ratio and/or large impact parameters, our technique will overestimate the S/N, but this has minimal impact on our assessment of PCs. The impact parameter is not well defined for low S/N events, but the transit depth and duration are measurable quantities.
We used S/N estimates for two purposes: to identify FAs and to determine a threshold for planet validation. The KOI catalog has an adopted S/N limit of 7.1 to classify a target as a KOI. FAs are present in the KOI catalog as initially the transit signal was estimated to have a S/N greater than 7.1 (Jenkins et al. 2002), and then as additional observations were gathered the S/N dropped below 7.1. These types of events inform us that validation of transiting planets with a S/N near the KOI threshold has a risk of introducing FAs, which we will now assess.
Using a S/N cut of 7.1 based on the transit depth and visual inspection, 26 KOIs in multis were classified as FAs: KOI-111.04, 439.02, 489.02, 966.02, 989.01, 1070.03, 1134.02, 1198.04, 1312.01, 1316.02, 1408.02, 1576.03, 1639.01, 1639.02, 1792.02, 1940.02, 1961.02, 2048.02, 2160.02, 2188.02, 2224.02, 2261.02, 2339.02, 2473.02, 2533.02, and 2586.02. For 15 systems, only a single PC remained, as indicated by the M1 column in Table 1 when the S/N cut is applied. The objects were considered single-planet systems for statistical counts. The rate of FAs from the single-planet and multi-planet population were both found to be ∼2%. As FAs do not represent real detections (quite the opposite), there is no reason to expect the rates to be predictable or reliable. The KOI creation process has been very inhomogeneous. This has caused the introduction of biases that favor finding and identifying additional PCs once the first candidate in a system has been found. This is especially true because of the notion that the FP rate for multi-planet systems is low. Quantifying this human bias is difficult and is part of our motivation to choose a larger transit S/N cut of 10 for planet validation.
The distribution of the transit S/N is shown in Figure 1. There is a rise in the observed number of PCs from a S/N of 50 to ∼15. The increase is driven by the increase in the number of PCs toward smaller radii and the increase in Kepler targets toward fainter magnitudes. A sharp drop is observed at S/N below 15, which marks the transition where the KOI catalog becomes significantly incomplete. We also inspected all transit candidates with a S/N less than 15 and found convincing transit signals for all candidates with a S/N greater than 10. Based on our observation of the S/N distribution shown in Figure 1 and inspection of the observed transit, we only validate PCs with a S/N greater than 10. We expect a large number of lower S/N candidates in the range of 7.1 to 10 to still be good PCs.
Figure 1. Histogram of the S/N for the Q1–Q8 KOI (first row of Table 1 sample). Approximately one-third of the sample has a S/N greater than 50.
Download figure:
Standard image High-resolution image5.3. Occultation/Secondary Eclipse Search
The primary signature of an EB relative to a transiting planet is the presence of eclipses of different depths due to the difference in surface brightness of the two orbiting stars. The change in depth can be quite dramatic depending on the nature of the two stars. As Kepler photometry has high precision and a spectral bandpass extending to 850 nm, occultations or eclipses can reliably be found for companions with radii similar to Jupiter with temperatures greater than approximately 2000 K. For bright host stars (Kepmag ∼ 10), this limit can be pushed to even cooler temperatures (e.g., TReS-2b; Barclay et al. 2012). As such, a secondary event can be due to a secondary eclipse from a stellar binary, or an occultation when a planet is blocked by the host star. To distinguish between secondary eclipses and occultations, we estimated the expected equilibrium temperature (Teq) for an orbiting body heated by incident stellar flux and compared it to an estimate of the temperature (Teffp) based on the depth of the occultation. The expectation is that a star, which is self luminous from nuclear fusion, will have a temperature Teffp that is much larger than Teq. We also test whether the depth of the occultation is consistent with reflected light from a planet by computing the geometric albedo, Ag in the Kepler bandpass. The secondary event is inconsistent with the planet hypothesis if Ag is significantly greater than unity.
Although visual inspection reveals some obvious occultations present in the data, we performed a more thorough search to identify occultations and eclipses. To search for secondary events, the light curve was phased to the orbital period and for each phase point the mean was calculated. Observations that occurred within one transit duration were compared to mean values computed at phases within ±1 transit duration. The difference divided by the standard deviation of observations at all phases was computed and used to identify occultations at any phase outside of transit.
To distinguish between planetary occultations and stellar eclipses, we compared the event depth with the occultation depth expected by a highly radiated exoplanet. We estimate the equilibrium temperature by
![Equation (5)](https://content.cld.iop.org/journals/0004-637X/784/1/45/revision1/apj490915fd5.gif)
where R⋆ and Teff are the stellar radius and temperature, a is the semi-major axis, AB is the planet's Bond albedo, and f is a proxy for atmospheric thermal circulation. To calculate the mean, we assume AB = 0.1 for highly irradiated planets (Rowe et al. 2006) and f = 1 for efficient heat distribution to the night side. The occultation depth was used to estimate the temperature of the companion (Teffp) using our best estimate of the stellar parameters,
![Equation (6)](https://content.cld.iop.org/journals/0004-637X/784/1/45/revision1/apj490915fd6.gif)
where Fp/F⋆ is the ratio of the companion and stellar flux and is equal to the depth of the occultation. We are assuming that the occultation depth observed over the Kepler bandpass is a proxy for the true bolometric flux ratio. We estimate uncertainties in Teq and Teffp by propagating our determined errors in the stellar parameters from Table 2. We also estimate whether the occultation could be due to reflection rather than thermal emission by estimating the geometric albedo,
![Equation (7)](https://content.cld.iop.org/journals/0004-637X/784/1/45/revision1/apj490915fd7.gif)
In the case that Teffp is greater than Teq at the 99.7 percentile (3σ) and Ag is greater than 1, we identify the event as a stellar eclipse and the candidate as an EB FP. While unexpected, such a test may classify self-luminous planets (e.g., youth or external forces) as FPs.
A number of FPs were detected through the identification of secondary eclipses (see Section 5.8). The only planet in a multi-planet system with a detected occultation was Kepler-10b with Ag < 1 and Teffp ∼ Teq. The lack of detected occultations in the multi-planet population is a consequence of the dearth of large, highly irradiated planets in these systems (see Figure 3). For the single-planet population, it is likely a handful of EBs that show occultations are classified as close-in Jupiter-sized PCs heated to ∼2000 K because such planets have an occultation and transit depth similar to an eclipsing low-mass star. The philosophy for the KOI catalog has been to keep a candidate classified as a PC until strong evidence is presented that shows the FP nature of the candidate. A consequence is that a handful of FPs will be misclassified as large (Jupiter-sized) candidates, which has no impact on our validation of multi-planet systems.
5.4. Phase Linked Variations
The orbital motion of a companion is imprinted in the photometric light curve due to day-night effects (thermal emission and reflectivity), ellipsoidal variations from gravity darkening due to tidal forces, and Doppler boosting from orbital motion. The latter two effects are dependent on the mass of the companion and, when present, can be used to estimate the mass of the companion (e.g., Mazeh et al. 2013).
There are cases when a companion interacts with the stellar surface through magnetic fields and produces star spots that could be misinterpreted as phase linked variations due to a massive companion (e.g., Tau Bootis; Walker et al. 2008). Thus our analysis would label a planet such as Tau Bootis as an FP.
To search for phase-curve variations (only relevant for short period systems), we filtered the data using the same procedure described in Section 4 except we changed the timescale of the polynomial fitter to five days instead of two days as the filter is destructive to astrophysical signals with a similar or longer timescale. This means our initial search is not sensitive to phase-linked variations on these longer timescales. The search was performed by calculating the occultation depth statistic introduced in Section 5.3, which is equivalent to average filtered data with a width equal to the transit duration. The standard deviation of the set of occultation measurements was calculated. This value was compared to the standard deviation of the data. This test determines whether the scatter is Gaussian on transit-duration timescales. When the ratio was found to be greater than two, variability in the phased light curve is detected. In these cases we inspected the light curves and found evidence of coherent EB effects as well as variable stars with fast pulsation or rotation timescales.
As was the case with the occultation search, it is expected that phase-linked variations will be rare because highly irradiated Jupiter-sized objects are rare in the Kepler multi-planet sample. KOI-1731.02 and 1447.02 were found to show phase curves and labeled as FPs. From the single-planet population, phase curves were discovered in KOI-23, 130, 143, 631, 636, 681, and 699. If ellipsoidal variations are source of the signal, then the mass of the companion was estimated, which ranged from ∼0.1 to 0.23 M☉, indicating that these sources are EB FPs.
5.5. Odd–Even Metric
The occultation/eclipse search is effective when the orbital period is correctly estimated. In cases when EB eclipses are similar in depth, it is common to have the period off by a factor of two. For most EBs, the depths of the alternating (odd- and even-numbered) transits differ. We search for this odd–even effect by separately modeling the odd- and even-numbered transits where we only allowed Rp/R⋆ to vary and the other parameters were fixed to their global solution. We used the change in Rp/R⋆ as a proxy for a change in transit depth. When the change was greater than 3σ, we inspected the transit light curves to insure that the effect was real. For the cases where we noticed spot-crossing-induced variations in the transit depths, the systems were retained as candidates.
From the multi-planet sample, only KOI-966.01 was found to exhibit an odd–even transit effect. Thus, the true orbital period is double the reported KOI value. From the single-planet population, 102 candidates had detected odd–even effects, although this count is incomplete as an FP is not always searched for additional effects or FP signatures in the light curve.
5.6. Centroid Analysis
A dominant source of FP planetary transit detections is EBs, or giant planet transits, on background stars that are captured in the aperture of the target star (BEB). These background signals are diluted by the target star and can have the appearance of small-planet transits. In this subsection, we describe the method we use to find KOIs with "clean" centroids, where the measurement is of high quality and there is no indication that the transit is not on the target star. This "clean" centroid standard, described in detail below, is a more stringent centroid standard than that used for PC status (see, for example, Batalha et al. 2013), and gives us confidence that the centroid signal is coincident on the sky with the target star.
We use centroid analysis to identify KOIs that are not clearly on the target star. The centroid method we use is the fit of a point response function (PRF) to the pixel difference image constructed by subtracting an average in-transit image from an average out-of-transit image (Bryson et al. 2013). This centroid method provides an offset from the target star position for each quarter, and the final offset for the KOI is a robust average of the quarterly offsets. We also use data quality metrics that indicate whether the data support the centroid offset measurement. We do not validate a KOI if the centroid offset suggests that the transit is not very close to the target star or if the data quality does not provide confidence that the transit signal is on the target star. This centroiding method does not work for highly saturated targets. Our treatment of saturated targets is described in Section 5.6.1.
We provide here a brief overview of the PRF-difference image centroid method. For details see Bryson et al. (2013). PRF-based centroids are measured on both out-of-transit and difference images quarter by quarter. When the target star is isolated, the centroid of the out-of-transit image gives the position of the target star. Assuming that the transit source is the only source of variability in the aperture, the centroid of the difference image gives the location of the transit signal source. These quarterly centroid measurements are robustly averaged to estimate the target star and centroid signal locations on the sky. These average locations are differenced to provide the average offset of the transit signal from the target. The robust average also provides a 1σ uncertainty per quarter, which is propagated through the robust average and offset calculation to provide an offset uncertainty. An alternative method for estimating the centroid uncertainty is via a bootstrap, using a resampling with replacement of the quarterly centroid measurements. The bootstrap-estimated uncertainty is also propagated through the robust average and offset calculation. We choose the larger of the two offset uncertainties when performing the cuts described below.
Centroid measurements are subject to several systematic errors, caused primarily by errors in the measurement of the Kepler PRF and crowding by background stars. The systematic error due to PRF error is mitigated by computing the offset as the difference between the PRF centroid of the out-of-transit image and the PRF centroid of the transit signal in the difference image. Because the transit source and the target star are near each other on the Kepler focal plane, their PRF errors are very similar so centroid systematic errors due to PRF error approximately cancel. The residual PRF error systematic varies from quarter to quarter and is statistically zero mean, so averaging over quarters further reduces PRF-error-driven centroid systematics. The residual systematics have a statistical standard deviation of less than 01. To account for this systematic error, a constant 0
067 is added in quadrature to the final offset uncertainty. This added constant does not, however, eliminate all apparently significant offsets due to systematic error, so we pass any KOI with offsets less than 0
3, even if that offset is formally statistically significant.
In a few cases, there is a field star in the target's aperture that is brighter than the target. In this case, the centroid of the out-of-transit image is strongly biased by the bright star, and the centroid offsets are invalid. We detect such cases by computing the offset of the out-of-transit image centroid from the catalog position of the target star, and declare the centroid measurement to be invalid if the offset of the out-of-transit image centroid from the target star catalog position is ⩾15.
We classify a KOI as having "clean centroids" if it passes three criteria, described in more detail below: (1) it has a good centroid measurement, (2) that centroid measurement indicates small offsets from the target star, and (3) there is at least a 99% probability that the transit signal is on the target star rather than another known star.
Good centroid measurement. The quality of a centroid measurement is determined by several factors, most notably the transit S/N and systematic error. We do not validate KOIs as planets for which difference images are not available. There are three ways in which a KOI can fail to have a good measurement.
- 1.When the S/N is very low, the measured offset uncertainty can be too large to sufficiently localize the transit signal. When the offset uncertainty is ⩾1
5, we say that the KOI does not have a good measurement.
- 2.The measured offset of the out-of-transit centroid from the target star's catalog position is ⩾1
5, indicating that it is likely that the out-of-transit centroid measurement is strongly biased by crowding.
- 3.The quality of the difference image in a quarter is determined by measuring the correlation of the difference image pixels with the fit PRF. If the correlation is less than 0.7, we consider the signal in the difference image too weak to trust the centroid value; otherwise we say that quarter has a good PRF quality. We demand that there be at least three quarters with good PRF quality, or that with two-thirds of the observed quarters have good PRF quality, otherwise we say the KOI does not have a good measurement.
Small offsets. We demand that the measured offsets be close to the target, satisfying both of the following criteria:
- 1.The offset is statistically close, that is the offset is <3σ, or the offset is <0
3 to allow for small systematic error.
- 2.The offset is smaller than 4''.
Probability ⩾99%. The systematic due to crowding is addressed via forward modeling of the observed pixels based on catalogs and the Kepler PRF (S. T. Bryson & T. D. Morton 2013, in preparation). A synthetic pixel scene is created for each quarter by placing a flux-scaled PRF at the pixel location of every known star close enough to contribute flux to the observed aperture. In this way, a synthetic pixel image modeling the average out-of-transit image is created for each quarter. A synthetic in-transit pixel image is created for each star in the aperture by reducing the flux of that star by the transit depth that best reproduces the overall observed transit depth, accounting for dilution. These images are analyzed for each star via difference-image PRF centroiding just like the observed pixels. The resulting offsets provide a prediction for the transit signal offset from the target star under the hypothesis that the transit occurs on each star in the aperture. The predicted offsets are compared with the observed offsets by inferring the underlying probability distributions. For each star in the aperture, the normalized integral of the product of the observed and modeled distributions provides a relative probability that the transit signal is at the same location as that star, when the modeled depth on that star is less than 100%. An unknown background source is also included as an alternative hypothesis. For details, see S. T. Bryson & T. D. Morton (2013, in preparation). For this paper we assume an underlying Gaussian distribution, which is characterized by the mean and uncertainty of the offset averages. We say a KOI is not clean if its relative probability is less than 0.99.
The KOI is considered clean if the measurement quality, small centroid offset, and probability criteria are all satisfied.
5.6.1. Manual KOI Inspection
Some KOIs considered in this paper do not have well-computed centroids, either because the KOIs are on saturated target stars or because the centroid did not satisfy the "good measurement" criterion. Some of these KOIs were subject to manual inspection based on the criteria described in this section. If they pass inspection, they are considered "good." We do not consider a "good" classification as strong as a "clean" classification, but as described in Paper II the multiple planet probability boost allows us to validate "good" KOIs.
Saturated target stars. When the target star is saturated or near saturation (Kp < 12), centroiding methods based on the PRF are no longer valid. In these cases the transit signal has a distinctive pattern in the difference image (Bryson et al. 2013). Visual examination of the difference image in each quarter provides a qualitative indication that the transit source is in the same pixel as the target star. Specifically, when the transit is on the target star, the transit signal appears at the end of the saturated columns, as well as in the non-saturated wings of the PRF. We pass a saturated KOI as "good" when the transit signal visibly appears as expected at the end of the saturated column and the transit signal wings match the non-saturated wings of the target star. All of the saturated multis considered in this paper for which there are difference images and which are not already confirmed planets passed this test.
Bad centroid measurement. When there are two "clean" KOIs in a system and additional KOIs that fail the "good measurement" criteria, manual inspection of the pixel data was performed to see if there is any indication that these additional KOIs are not at the target star location. The typical situation is that the difference images were too noisy to support a high-quality centroid measurement. In this case, when manual inspection indicates that the transit signal is on the same pixel as the target star, and that there is no significant signal in the difference image away from the target star, we consider the KOI to be "good."
5.7. KOIs with Validation Issues from Imaging
Table 5 lists candidates that have newly detected companions inside the photometric aperture by Adams et al. (2012, 2013) and one of us (S.H.). Because these companions were not in the catalog used to compute the probability criterion described in Section 5.6, we give these special attention. Table 5 gives the observed offset of the newly found companion from the target star and the offset of the companion from the measured transit source in units of the centroid uncertainty. When the companion star is more than 3σ from the measured transit source, we consider that companion as ruled out as a source of the transit signal. Because our validation criteria includes the requirement that the centroid source be no more than 3σ from the target star, companions outside of 4σ will not reduce the KOI's probability of being on the target star to less than 99%. We do not consider whether or not the companion is gravitationally bound to the primary star.
Table 5. Close Neighbors Discovered with High-resolution Imaging
KOI | Companion Offset | Transit Position σ | Companion Offset from Transit Position in σ | Disposition |
---|---|---|---|---|
(arcsec) | (arcsec) | |||
112.01 | 0.1 | 0.1840 | 1.0 | fail |
112.02 | 0.1 | 0.2900 | 0.4 | fail |
270.01 | 0.05 | saturated | <3 | fail |
270.02 | 0.05 | saturated | <3 | fail |
285.01 | 1.44 | saturated | <3 | fail |
285.02 | 1.44 | saturated | <3 | fail |
279.01 | 0.922 | saturated | <3 | fail |
279.02 | 0.922 | saturated | <3 | fail |
307.01 | 0.080 | 0.2745 | 1.1, 1.1 | fail |
307.02 | 0.080 | 0.2333 | 0.6, 1.3 | fail |
102.01 | 2.76, 5.45 | 0.7886 | 6.0 | pass |
102.02 | 2.76, 5.45 | 0.7886 | 6.0 | pass |
123.01 | 2.03, 5.27 | 0.1950 | 11.3 | pass |
123.02 | 2.03, 5.27 | 0.2290 | 9.2 | pass |
124.01 | 2.4 | 0.1945 | 12.7 | pass |
124.02 | 2.4 | 0.2187 | 11.8 | pass |
153.01 | 5.14 | 0.1104 | 53.7 | pass |
153.02 | 5.14 | 0.1186 | 61.8 | pass |
251.01 | 3.45, 4.76 | 0.1026 | 44.5 | pass |
251.02 | 3.45, 4.76 | 0.6723 | 9.6 | pass |
283.01 | 5.96 | saturated | >3 | pass |
283.02 | 5.96 | saturated | >3 | pass |
555.01 | 4.01 | 0.3456 | 25 | pass |
555.02 | 4.01 | 0.1852 | 12.5 | pass |
298.02 | 0.825 | 0.9002 | 1.4 | fail |
Notes. KOIs with otherwise clean centroids that have close companions revealed by high-resolution imaging. For non-saturated targets, the companion offset in σ is companion's offset divided by the uncertainty in the centroid measurement for that KOI. When the companion is within 3σ of the measured transit location we do not validate the target. Two distances in σ are given for KOI-307 because there is a 180 deg ambiguity in the position angle of the companion.
Download table as: ASCIITypeset image
When the companion stars are within 3σ of the transit position, the transit signal is not necessarily an FP. However, this indicates that we did not determine which star was the source of the transits. We do not validate such candidates unless we have strong evidence that the nearby star is a bound companion to the Kepler target as described in Section 9 of Paper II.
5.8. Identified FPs Among the Multis
All of the FP tests described above have been applied to the original sample of 1212 PCs identified as potential multi-planet systems. For each FP, a brief description of the types of FPs detected in multi-planet transiting systems is presented below. The FP disposition was used for a comparison of the single planet, multi-planet, EB, and FP samples. There are three classes of FPs used to describe the nature of the transiting object: (1) Period and epoch (P/T0) collisions where multiple sources show the same orbital period and transit times. Such events can be produced by direct PRF contamination, optical reflections or electronic interference, such as crosstalk (Coughlin et al. 2014); (2) flux FPs, where the photometric light curve shows evidence of an EB; and (3) active-Pixel-Offsets (APOs) FPs, where centroid measurements of the photometric aperture indicate that the source of the transits is due to a source offset from the Kepler target. Categories (2) and (3) are not mutually exclusive. Table 1 gives a breakdown of the PCs for various cuts and number of systems with six candidates, five candidates, four candidates, and so forth.
There were 12 P/T0 collisions detected in the following KOIs: 376.01, 489.01, 989.02, 1119.01, 1196.01, 1231.01, 1231.02, 1747.02, 1803.02, 1806.01, 1944.02, and 2188.01. A secondary eclipse was detected for most of these systems, indicating that the primary sources of P/T0 collisions are EBs. The distribution of P/T0 collisions will not favor transit candidate targets, thus it is expected that the rate of P/T0 collisions is lower for multi-planet sample relative to the single planet sample. KOI-489.02 was flagged a FA with a transit S/N of 6.6, thus the KOI-489 system does not count as a multi-planet system in any of our statistical counts.
KOI-199.02 shows a secondary eclipse with a depth of 50 ppm and an orbital period of 8.8 days. The depth of the secondary eclipse is inconsistent with the planet hypothesis. Derived values from the occultation are Ag ∼ 22, Teff ∼ 4500 K and Teq = 1200 K.
KOI-376.01 is a P/T0 collision and shows strong quarterly depth variations due to quarterly variation of contamination. The second candidate, KOI-376.02, has a period of 1.4 days and shows a strong secondary eclipse with an observed depth of 260 ppm. Since the signals observed are likely heavily diluted, the true eclipse depths are likely much deeper.
KOI-379.01 was found to have an additional star within the photometric aperture with a separation of 1''. Centroid analysis points toward the fainter star as the source of the transits, thus this candidate has been flagged as an APO. Centroids analysis of KOI-379.02 is inconclusive to determine which star is the source of the transits and is kept as a PC.
KOI-414.01 shows a secondary eclipse with a depth of 400 ppm. The location of the secondary eclipse indicates that the orbit is non-circular. KOI-414.02 shows a clear centroid offset and was flagged as an APO.
KOI-2671.01 was marked as an FP as a secondary eclipse was detected. The occultation shows the orbit to be eccentric (35 hr offset). KOI-2671.02 has a centroid offset and was labeled as an APO.
KOI-989.01 and KOI-989.02 are the same event with the periods being integer multiples of one another. These two candidates were also flagged as a P/T0 collision. The confusion of these two candidates arose because of strong variations in the quarter transit depths from quarterly dependence of dilution. KOI-989.03 remains as a PC.
KOI-549.01, 549.02, 1196.01, 1231.01, 1378.02, and 2007.01 are flagged as APOs from centroid analysis. KOI-1119.01 is a P/T0 collision and 1119.02 shows strong centroid offsets. KOI-1342.01 shows a small offset of 09 with a significance of 4.1σ. It is therefore considered to be an APO. KOI-2159.02 shows centroid offsets and a secondary eclipse. KOI-1731.02 shows an occultation and phase-linked variations and was labeled as an FP because its transit depth appears to be heavily diluted.
KOI-966.01 shows an odd–even transit effect, but KOI-966.02 was labeled as an FA due the low S/N of the transit event. Thus, the KOI-966 system does not count as a multi for our statistics.
KOI-1447.01 showed a "V" shaped transit event with a depth greater that 15%. While transit-depth is not an indication of the FP nature of the candidate, KOI-1447.02 shows large amplitude phase-linked variations. Thus, KOI-1447.02 is a clear FP, which removes the KOI-1447 system as a multi-planet system. Due to the large transit depth, we classified 1447.01 as an FP.
From the single-planet population of 2482 PCs, 976 were classified as FPs resulting in an FP rate of ∼40%, which, as expected, is in stark contrast to the multi-planet FP rate. In total we found 26 FPs (including P/T0 collisions) in the multi-planet sample of 1167 PCs remaining after the removal of FAs and single planets. A few of the classified FPs had an associated FA, such as KOI-1231.02, and are included in the FP totals for the single planet population. The 26 FPs include cases of two FPs associated with the same target. Candidates that were flagged as having not clean centroids in Table 3 are not FPs and remain as unvalidated PCs. There is no strong evidence to suggest that any one target with not clean centroids is a blend; however, the probability of blends existing within the population is large enough that we cannot validate this sample at the 99th percentile. The 26 FPs are found around 20 systems, with 6 double FP systems, 12 cases of FP associated with a single PC, and 2 cases where 2 PCs and 1 FP where associated with the same target.
The results of the multi-planet disposition are summarized in Table 4 of Paper II, which gives a comparison and breakdown of the expected FP rate. The agreement between the observations and predictions is very good which leads to our conclusion that a vast majority of transiting candidates found in multi-planet systems are genuine planets. After removal of FPs, there are 1129 remaining multi-planet candidates. The next step is to explore this large population of candidates and set additional criteria to reduce the chances that undetected blends still exist. There will be a population of FPs from blends that exist in the Kepler transit sample, but cannot be detected via our methods, for example, a blend from a BEB where the separation from the target source is too small to be detected by centroid motion.
5.9. Validation of Multi-planet Candidates
To reduce the number of potential FPs in our validated list of planets in multi-planet systems, we eliminate regions of phase space where we have reduced confidence. For example, we have reduced confidence in the validity of a PC if centroids cannot localize the position of transit to eliminate the chance of a background blend at the 99th percentile. The first requirement is that a candidate has a S/N > 10 as established in Section 5.2. This insures that the multi-planet sample is free of FAs and removes an additional 21 candidates after the removal of FPs and P/T0 collisions.
Using the analysis for Section 3, Section 4, and this section, we are able to use the criteria set out in Paper II to select a population of multi-planet transiting systems that have an FP rate substantially less than 1% (additional details below). Table 1 lists the number of PCs after various cuts and tests are applied. The various cuts are: FA, where either a transit candidate has insufficient S/N (<7.1) or was labeled as a non-transit event such as stellar variability, or was observed to have less than three transits. Col indicates a P/T0 collision. These sources are non-unique by nature so they are classified separately. FP indicates when an FP is identified that is a not a P/T0 collision. An FP can either be an EB masquerading as a PC or a diluted signal where the source of the transits has been localized off the Kepler target. SN: the transit models were used to determine the S/N of the phase folded transit for each candidate. We adopted a threshold of S/N > 10 to consider a transit event. P marks period cuts. We require the orbital period to be greater than 1.6 days due to the increased rate of FP found with shorter orbital periods (see Figure 2 in Section 5.8 and Appendix A of Paper II). b marks when cuts are made based on the measured impact parameter. When a transit is "V"-shaped, there is a larger chance that a FP has been identified. This does not mean that "V"-shaped events are FPs, only that we have less confidence in declaring such objects as planets. The fraction of "V"-shaped signatures that are produced by EBs as opposed to transiting planets is far larger than that for "U" shaped profiles. Our criteria for "V"-shaped transits is that b + bσ + Rp/R⋆ > 1.00. centr indicates the centroid test as outlined in Section 5.6. A target that does not pass our centroid test is not a statement that the object is an FP/APO, but rather that we had insufficient information to localize the source of the transits on the Kepler target. The column number of multi-pl indicates the total number of multi-planet systems that pass the indicated test. The column number of new multi-pl indicated the total number of multi-planet systems that pass the indicated test and have been previously verified (already assigned a Kepler ID).
Figure 2. Cumulative distribution of orbital periods for multi-planets (black), single KOIs (red), FP KOIs (green), P/T0 collisions (cyan), and EB systems (blue). All periods used in normalization. Other timescales are shown in Figure 2 of Paper II.
Download figure:
Standard image High-resolution imageFigure 2 shows the cumulative distribution of orbital periods for candidate multi-planet systems (black), candidate single planet systems (red), FPs KOIs (green), P/T0 collisions (cyan), and EBs from the Kepler Eclipsing Binary catalog (blue).21 FA are not considered. The distribution of EBs shows a large population of short period events due to the inclusion of contact binaries. The sample of FPs from the KOI catalog also shows a relatively strong population of FPs at short orbital periods, which is different from the EB population, but there is still a much larger fraction of FPs at shorter orbital periods compared to the PC populations. There are two reasons for this difference. (1) KOIs are selected based on a visual inspection of the photometric transit event, with a requirement that the event has an appearance of a planetary transit. This process heavily reduces the number of contact binaries in the KOI catalog, as a distinct transit that shows a clear ingress and egress are not present. (2) The Kepler transit search algorithm does not conduct searches for events with periods less than 0.5 d, so unless a strong harmonic of the orbital period is detected at a longer period, many short period events will also be missed. Similarly, the distribution of P/T0 collisions, which is dominated by EBs, shows a larger fraction of events with short orbital periods relative to the PC distributions. As articulated in Appendix A of Paper II, the expected abundance of unidentified FPs in multis to planets in multis is far larger at small orbital periods than at large ones. Thus we do not validate PCs with orbital periods less than 1.6 days. We also excluded systems with an orbital period less than 4 days and a S/N < 15, due to concerns of an increased FP rate and the lack of well constrained transit model parameters due to the low S/N of the transit. Eliminating candidates with a P < 1.6 days reduced the sample of 1107 PCs in multis that passed our FP and S/N tests to 1084 PCs as shown by row 6 of Table 1.
Dispositioning of PCs relies on transit models to characterize the orbiting companion. A transit can be "V"-shaped in appearance when the transit duration is comparable to the photometric cadence or we have a grazing transit (b + Rp/R⋆ ∼ 1). The transit model incorporates the cadence time to convolve the synthetic lightcurve to match observations and allows a quantitative assessment as to whether a grazing transit is observed. A grazing transit also results in increased uncertainty in transit model parameters, which makes assessment of the transit event difficult. From the single-transit population, after the removal of FAs, 16% have grazing or close to grazing transits (b + bσ + Rp/R⋆ > 1.00), which drops to 8% after the removal of FPs and P/T0 collisions, which is larger than one would expect based on an isotropic distribution of orbital planes. Almost half of the grazing-transit single PCs have a radius greater than 15 R⊕, which indicates that a large fraction of this population is likely FPs. However, planetary radius is not a criterion to label a candidate as FP because an upper limit on the radius of a planet is not well established, and measurement uncertainty on "V"-shaped transits precludes making a definitive statement regarding the absolute radius of the transiting object. From the multi-planet population, after the removal of FAs, 3.8% are measured to be grazing or near grazing, which drops to 3.1% after the removal of FPs and T/P0 collisions. Only 1 (5.9%) of the grazing multi-planet candidates has a large inferred planetary radius, KOI-1477.01, which is associated with an EB (KOI-1477.02). It is very likely that KOI-1477.01 is also an EB FP. The lower rate of grazing transits in the multi-planet population leads us to conclude that a large fraction of the grazing transits in the single-planet population are FPs. This means that there is a higher probability of a FP being found when a grazing transit is present, thus we do not validate multi-planet candidates that have b + bσ + Rp/R⋆ > 1.00. These KOIs are still good PCs. Application of FP, P/T0, S/N, and cuts based on impact parameter reduces the set of multi-planet candidates for validation to 1054.
Cuts based on S/N, period, and transit shape use the transit models and comparison with the single-planet population to identify regions of model parameter space with reduced confidence in the validity of a PC when the rate of FPs is observed to be larger relative to the multi-planet population. Cuts based on methodology are presented in Section 5.6, where PRF models were used to identify which multi-planet candidates have at least a 99% probability that the transit signal is on the target star rather than another known or unknown star. This statement means that there is a low probability that a blended background transit event is present. After application of our centroid criteria, the number of multi-planet candidates still considered for validation is 851.
As previously mentioned, planet radius is not a criterion for classification of a transiting candidate as an FP. A problem thus arises for large Jupiter-sized PCs, as there is a degeneracy in radius for planets, brown dwarfs, and low-mass stars. Additional dynamical tests were applied to the multi-planet candidates: using Hill's criterion to test for stability of neighboring pairs of PCs, and when , dynamical fits to transit times observed in Q1–Q14 Kepler long cadence data were conducted to determine whether the giant candidates can have masses exceeding 13 MJ, assuming all of the candidates orbit the Kepler target. Details of both tests can be found in Appendix C of Paper II. All candidates that passed the above tests, apart from the special case of KOI-284 (Kepler-132, addressed in Section 9.1 of Paper II), passed the stability tests. Table 3 categorizes candidates as having passed, failed, or being too small to have been tested for mass limits large enough to be stars. The last row of Table 1 gives the number of planetary candidates after all tests have been applied and gives us a total of 851 planets that we validate. From this sample, 60 have been previously validated via other methods, thus we are able to introduce 768 newly validated planets, which roughly doubles the current number of confirmed and validated planets. Planets discovered and confirmed by the Kepler mission currently account for more than half of the known and validated extrasolar planets.
6. POPULATION OF VALIDATED PLANETS AND DISCUSSION
After the application of all the tests listed above, we validate 851 extrasolar planets associated with 340 planetary systems. These systems are expected to have an FP rate that is significantly less than 1% due to the reasons listed in Section 5.9 and the theoretical framework laid out in Papers I and II. From this sample, 768 candidates in 306 systems have not been previously validated, but are now extrasolar planets validated above a confidence level of 99%. Thus, we introduce Kepler-100 through 405. From this population there are 106 new planets that have a radius less than 1.25 R⊕, compared to 16 that have been previously validated. There are 6 planets with incident solar flux S less than 1.5 times that of the Earth, including four planets: KOI-518.03, 1422.04, 1430.03, and 1596.02 (Kepler-174 d, Kepler-296 f, Kepler-298 d, and Kepler-309 c, respectively) that are new validations that we discuss below. Figure 3 plots incident flux versus radius and displays our new validations as filled circles. Multi-planet candidates that pass all tests except our centroid criterion are plotted as open circles and single PCs after the exclusion of FPs are plotted as small dots. The falloff in the number of planets below 1 R⊕ is driven by incompleteness due to insufficient S/N. The falloff in the number of planets with S < 1.5 is due to decreasing transiting probability and incompleteness to longer period events (>150 d) as our sample is based on Q1–Q8 photometry (∼two years). As noted by Latham et al. (2011), there is a lack of hot-Jupiters in multi-planet systems verified in Figure 3 as a deficit of planets with S > 200 and Rp > 7 R⊕ relative to the single-planet population.
Figure 3. Multi-planet population showing incident flux S (relative to the flux received by the Earth) vs. planetary radius (R⊕). Planets validated in this paper are marked with filled circles, unvalidated planets in multi systems are marked with open circles, and planet candidates found in single systems are plotted with dots. All planets and planet candidates shown have P > 1.6 days, b + bσRp/R⋆ < 1.00, and have not been identified as a FA, FP, or P/T0 collision. The two colored bands display estimates of the inner and outer habitable zone based on the Recent Venus and Early Mars models from Kopparapu et al. (2013). The colors display the range each boundary as a function of Teff from 3000 K (red) to 7000 K (light blue).
Download figure:
Standard image High-resolution imageThe FP tests presented in Section 5 are not sensitive to most hierarchical blends. A hierarchical blend is a bound stellar binary with a transiting planet orbiting one of the stellar components. It is not known if widely separated binaries host planetary systems with orbital planes aligned with the stellar orbital plane. If the alignment distribution is isotropic, then hierarchal, transiting triples in the multis may be rare, however, if alignment is common, say because of star-planet formation processes that favor aligned systems, then the rate of hierarchal, transiting triples could be much larger. For an isotropic distribution it was shown in Section 5 of Paper II that four or five systems are likely to have a PC around each stellar component and it is very unlikely that there is more than one hierarchical triple multi (three bound stars, one hosting a transiting planet and the other two in a eclipsing configuration). If the orbital planes of planets around both components of a stellar binary are aligned, then we might expect to find a greater number of blends. In Appendix A, we develop a synthetic population model to test whether hierarchical blends could contribute a large fraction of the observed multi-planet population by comparing the measured value of ρ⋆ from our transit models. We find that the multi-planet population is not dominated by hierarchical blends, but the strongest constraints come from KECK HIRES observations to search for spectroscopic blends (Section 8.1 of Paper II). Together, it appears the rate of hierarchical blends is low. It is important to note that even if any of our validated planets are found to be orbiting a fainter and bound star, they are still planets; however, the stellar parameters listed in Table 2 will need to be revised.
The single- and multi-planet populations also appear to have different fractions of planets at longer orbital periods. The relative cumulative distribution of the multiple planet systems overtakes the single planet populations at a period of ∼25 days. There are 1027 and 897 single- and multi-planet candidates with periods between 5 and 150 days. If we separate these samples into short periods, 5–25 days, and long periods, 25–150 days, we find 689 and 627 short period planets and 338 and 270 long period planets for the single- and multi-planet populations, respectively. Thus, 55% of the multi-planet sample are found in the short period bin compared to 45% for the single-planet sample. Explaining this difference seems counterintuitive, as alignment of orbital planes in multi-planet systems would make it more likely to find longer period planets relative to the single-planet population under the assumption that long period planets are equally common in multi- and single-planet systems. However, there are strong biases in the detection process that generates the raw KOI list. In particular, the candidates are found via different numerical and inspection methods. For example, there are a number of long period single-transit candidates that were found through identification of a single transit event and then the candidate was continuously monitored for additional transits.
6.1. New Planetary Systems with Planets in or near the Habitable Zone
We discuss herein, validated multiple-planet systems that contain a planet that is in the nominal habitable zone of their star. The location of the habitable zone depends on stellar luminosity (and the orbital period range also depends on stellar mass), so we introduce only those planets whose host stars have been characterized spectroscopically in this section.
As we do not have information regarding either the albedo or atmospheric characteristics of these planets (nor of any moons that they might have), we can only make reasonable estimates of the flux of stellar radiation that they intercept, i.e., the amount of insolation that they receive. We therefore quote results in terms of the average solar flux intercepted by Earth, S, which is generally referred to as the solar constant. For the purposes of our tabulation, we list objects that intercept flux less than 1.5 S (comparable to the flux Venus received 1 billion years ago, see Kopparapu et al. 2013 and references therein). For comparison, Venus receives 1.91 S. We don't specify an outer boundary to the HZ because few of the Kepler planets that we have validated have significantly smaller insolation than does Earth, but note that Mars receives an average of 0.43 S. Orbital eccentricity, e, which generally is unknown, affects the flux of stellar radiation that a planet receives, but the change that it induces in the annual average insolation is roughly quadratic in e, and as few of the planets in Kepler's multis have large eccentricities, the magnitude of the change in mean insolation resulting from planetary eccentricity is likely to be small. The spectrum of stellar radiation received by a planet also affects atmospheric and surface temperatures (Kasting et al. 1993), but these variations are small compared to uncertainties in estimated insolation and in-atmospheric properties. Nonetheless, we note the effective temperatures of the stellar hosts to aid in investigations by other researchers. Only six of the planets that we validate intercept less than 150% of the radiation flux encountered by Earth. Two of these orbit KOI-701 (Kepler-62) and have been analyzed in detail by Borucki et al. (2013). The transits for the four new planets that receive less than 1.5 S, KOI-518.03 (Kepler-174 d), 1422.04 (Kepler-296 f) (see Paper II for more detailed discussion), 1430.03 (Kepler-298 d), and 1596.02 (Kepler-309 c), are shown in Figure 4. The planets have nominal radii of 2.19, 1.79, 2.50, and 2.51 R⊕.
Figure 4. Transits of all of the validated planets in systems with a newly validated planet with S < 1.5. The beauty shots in the upper panels display the sizes of the stars and planet candidates to a uniform scale. The color of the stars and the impact parameters of the planetary transits reflect estimates of stellar and transit characteristics given in Tables 2 and 3. Verified planets are shown in black while other candidates are green. Planets and candidates are displayed with distance below the middle of the star corresponding to the transit impact parameter. The lower panels show the detrended Kepler flux from the host star phased at the period of each transit signal and zoomed to a region around mid-transit, shown in order of increasing orbital period. Black dots represent individual Kepler long cadence observations. The blue bars are the data binned 30 minutes in phase with 1σ uncertainties. The colored curves show the model transit fits, with colors corresponding to the last two digits of KOI designators as follows: red = .01, green = .02, blue = .03, cyan = .04. In each panel, the best-fit model for the other planet candidates was removed before plotting. All panels for a given system have an identical vertical scale to show the relative depths, and identical horizontal scale to show the relative durations, but scales differ between systems. The successive panels show KOIs 518, 1422, 1430, and 1596.
Download figure:
Standard image High-resolution image6.2. Conclusions
Our work provides a substantial increase in the number of verified exoplanetary systems and demonstrates the ability of the Kepler mission to probe the statistics of exoplanetary systems with a sample that is relativity clean of FPs. Both transit models and centroid models are used to characterize the photometric data and various tests were used to identify FPs. The rate of FPs was found to be low relative to the single transiting planet population, in quantitative agreement with theory (Papers I and II). This result demonstrates that 851 PCs in multi-planet systems are valid planets at greater than the 99% level. In Appendix A, the multi-planet population was used to investigate the rate of hierarchical blends in multi-planet systems; while no limits on the occurrence rate can be currently set, it was found that the eccentricity distribution of transiting multi-planet systems found by Kepler is significantly different from the planetary distribution found by RV surveys. The list of validated planets presented is reliable, but the sample suffers from both incompleteness and strong biases. Many of the candidates that were not validated in this study are still excellent planetary candidates.
Funding for this Discovery mission is provided by NASA's Science Mission Directorate. We are indebted to the entire Kepler team for all the hard work and dedication that has made such discoveries possible. J.F.R. is partially supported by a Kepler Participating Scientist grant (NNX12AD21G).
Facility: Kepler - The Kepler Mission
APPENDIX A: RATE OF FALSE MULTIS FROM TRANSIT DURATIONS
The FP tests presented in Section 5 are not sensitive to all types of hierarchical blends. We consider a hierarchical blend to be a gravitationally bound stellar binary with a transiting planet orbiting one of the stellar components. In this section we lay out the framework for estimating the number of blends in the Kepler multi-planet sample. Our aim was to set limits on the rate of hierarchical blends, which could be a large source of error in the Kepler exoplanet database. We present two types of analysis: (1) comparison of measured values of ρc for each planet pair within a multi-planet candidate system, and (2) comparison of ρc to ρ⋆ derived from stellar theory.
Our transit models provide a measurement of ρc, which is the measurement of the mean stellar density, ρ⋆, if the transiting planet travels on a circular orbit. The value of ρc is strongly correlated to the transit duration, which in turn depends on the planet's semi-major axis and the impact parameter for a circular orbit. Each transiting planet in a hypothetical planetary system with planets in perfect circular orbits will produce the same measurement of ρc. Variations in the measured value of ρc can be produced in three ways. (1) Eccentric orbits change the transit duration depending on the star–planet separation at the time of the transit in accordance with Kepler's Second Law. Comparison of a population of transits provides some insight on the eccentricity distribution because ρc and ρ⋆ will differ. (2) Unresolved, gravitationally bound stars that host transiting planets that have diluted transits as both stars are observed within the same photometric aperture. The unseen companion is the fainter and, generally, smaller relative to the target star. Thus, hierarchical triples produce systematically larger values of ρc. When the planet system is bound to the fainter star, the measurement of ρc can disagree with the estimate of ρ⋆ from stellar classification and modeling. If transiting planets are found around both components of the stellar binary, then ρc will not agree planet to planet. (3) Measurement error introduces scatter that produces ρc values that are equally too small or large. The population of Kepler multi-planet systems was used to place limits on the rate of hierarchal blend occurrences as described below.
The difference between values of ρc for each KOI was computed for each PC in the system. For example, if a system has three planets (P1, P2, and P3), then we would compare the difference in ρc based on P1-P2, P1-P3, and P2-P3 relative to the sum of ρc for each pair,
![Equation (A1)](https://content.cld.iop.org/journals/0004-637X/784/1/45/revision1/apj490915fd8.gif)
Since the distribution is symmetric, we used the absolute value. The binned distribution of dρc, i − j measurements is shown in black in Figure 5. The sample included 1158 PCs after the removal of FAs, FPs, and P/T0 collisions. The observed distribution is broad but peaked toward zero. The next step was to construct a synthetic population to reproduce the observed distribution that accounted for eccentricity, binarity, and measurement error.
Figure 5. Distribution of differences in ρ⋆ as measured by two planets orbiting the same star. The observed Kepler multi-planet population is shown in black with 1σ Poisson uncertainties. Other lines show synthetic population models. All models incorporate the measurement uncertainty of ρ⋆ from Table 3. The cyan line shows a synthetic model population with circular orbits and no hierarchical blends. The red line shows a synthetic population with an eccentricity distribution based on RV detected planets. The green line shows a synthetic population with eccentricity and a hierarchical blend rate of 0.5. The stellar companion masses are drawn from a 1/q distribution. The blue line shows a synthetic population similar to the green line except the stellar companion masses are drawn from a uniform distribution.
Download figure:
Standard image High-resolution imageA.1. Constructing Synthetic Populations
To construct a synthetic population, ρ⋆ for the primary star was adopted from Table 2 and matching orbital periods for the orbiting planets from Table 3. Thus, a synthetic transiting planet is generated for each planet in our transiting multi-planet sample. A co-eval binary companion is generated for each star and later we decide whether the primary or secondary component is hosting the planet. The companion is not used to estimate the binary fraction, but to determine the change in ρc and estimate the number of hierarchical blends. A hierarchical blend occurs when transiting planets are observed around both components of a stellar binary. Masses for a bound companion (M2) were chosen to be greater than 0.1 M☉ and less than the primary (M1) and have a mass ratio (q) distribution,
![Equation (A2)](https://content.cld.iop.org/journals/0004-637X/784/1/45/revision1/apj490915fd9.gif)
where n = −1 would produce a 1/q distribution matching the distribution observed in RV surveys (Trimble 1990). The transit depth from Table 3 together with an estimate of the luminosity of the primary (L1) and secondary (L2) stellar components was used to check that the undiluted transit depth around the fainter secondary star would not exceed 50%. This sets a lower limit on the luminosity of the secondary,
![Equation (A3)](https://content.cld.iop.org/journals/0004-637X/784/1/45/revision1/apj490915fd10.gif)
where Tdep is the transit depth from Table 3. If L2 < L2, min, then we choose another mass for the stellar companion and repeated until a suitable choice is found for this check. To determine which stellar component the planet would be orbiting in our model, a fitted parameter, binfrac, was used to represent the fraction of planets that are orbiting the fainter, bound companion. For each transiting planet in the system we drew a uniform random number. If that number was less than binfrac, then we adopted ρ⋆ of the bound companion.
A system has now been constructed that consists of two bound stars, each having a probability of having a transiting planet. To account for eccentricity in the synthetic population, a two-parameter model of the eccentricity distribution based on the beta function, as described in Kipping (2013), was adopted,
![Equation (A4)](https://content.cld.iop.org/journals/0004-637X/784/1/45/revision1/apj490915fd11.gif)
where e was the eccentricity of the planet and x and y were fitted parameters. Other distributions, such as a Rayleigh distribution, could also be used to produce similar results. The distribution Pβ was used to draw a value of e for each transiting planet. With the orbital period and ρ⋆ selected, a/R⋆ was calculated using Equation (1). The argument of periapsis (ω) was then randomly selected from a uniform distribution and used to determine the star–planet separation, d/R⋆, during transit. The transit probability was then calculated,
![Equation (A5)](https://content.cld.iop.org/journals/0004-637X/784/1/45/revision1/apj490915fd12.gif)
and a uniform random number from 0 to 1 was drawn. If the random number was greater than Tprob then the choice of ω was rejected and a new value was drawn and the exercise repeated. This process insures that transits occurring near periastron are preferred. The estimate of d/R⋆ was substituted in Equation (1) to estimate ρc,
![Equation (A6)](https://content.cld.iop.org/journals/0004-637X/784/1/45/revision1/apj490915fd13.gif)
Measurement error was incorporated by choosing a model solution from the MCMC analysis for the corresponding KOI and comparing ρc to the median value of ρc from all the chains. The difference was added to the synthetic value of ρc. To investigate the dependance of the synthetic model on reliability of estimating uncertainties in transit parameters, we used a nuisance parameter, errfrac, to scale the errors on ρ⋆ as measured by the transit when assuming a circular orbit.
A.2. Results
In Figure 5, we plot the binned distribution of dρc, i − j based on various synthetic populations for comparison to the observations. The cyan line shows a population of planets with only circular orbits and no hierarchical blends; only incorporating measurement error. This model does not match the observations shown in black. The red line was produced by incorporating measurement error and an eccentric distribution of planets from RV planets (Wright et al. 2011) with best-fit parameters from Kipping (2013). This model produced a better fit, but not an ideal fit to the observations. The green line shows a population produced using measurement error, circular orbits and a hierarchical blend rate of 0.5 and N(q)∝q−1. In this scenario, half of the planets are transiting primarily low-mass stellar companions. The blue line shows a population with circular orbits and a hierarchical blend rate of 0.5 and a uniform distribution of q. For both cases with hierarchical blends, we see an overabundance of mismatches in ρc between planet pairs, with the strength of the mismatches modulated by the distribution of q.
To measure posterior distributions of the parameters to describe the planet population we use a MCMC routine that uses methods similar to the description found in Section 4.1. In Figure 6 we show distributions for binfrac, n, a, b, and errfrac based on 48,000 chains. Parameters where restricted to binfrac = {0, 0.5}, n = { − 2: 2}, a, b = {0, 10} and errfrac > 0. It is immediately clear that posterior distributions for each fitted parameter are quite broad, but we can draw a few conclusions.
Figure 6. Plots of MCMC distributions for parameters describing a synthetic population of planet systems to measure the rate of hierarchical blends. The diagonal panels show histograms for each parameter, other panels show scatter plots of two parameters against one another.
Download figure:
Standard image High-resolution imageAn extensive search for blended companions based on KECK high-resolution spectra was described in Section 8.1 of Paper II. The sample included 270 multi-planet candidate systems and would be sensitive to blends due to companions that are 2%–3% as bright as the companion and show a RV difference of ∼10 km s−1. From this sample, only one clear blend was found, and in that case (KOI-2311.02) the S/N of the transit was found to be very low, so we do not even consider that system to be a multi-planet candidate. However, based on the one potential blend detection, an estimate of the blend rate is 0.004 ± 0.004 for companion stars that are 2%–3% as bright as the primary. Beyond ∼5 AU, the RV component of the stellar binary will be too small to allow reliable detection of stellar blends, which would account for less than half of companions found in solar-type stars (Raghavan et al. 2010). We double the potential blend rate and take the 3σ upper limit to get a rough estimate of the number of hierarchical blends that could exist. The synthetic population has 1158 planets in 460 systems, so observations suggest that no more than 14 blended systems could exist and either be missed by the spectroscopic survey or have separations large enough that a companion would not be detected. The inclusion of high-resolution imaging observations would allow additional constraints on the number of companions detected at large separations. In Figure 6, all simulations that have less than 14 blended systems (e.g., at least one planet around each stellar component) have been marked in red.
The rate of hierarchical blends, binfrac, was found to be dependent on how well uncertainties are determined for ρc. The transit model used a fixed set of limb-darkening parameters, thus we expect our uncertainties to be both underestimated and potentially systematically biased. We do not expect our uncertainties to be underestimated. There is weak (2σ) evidence that there are zero hierarchical blends in our sample based on the measurement uncertainties in Table 3. Limits on the blend rate from spectroscopic analysis also suggests that the uncertainties on ρc have been underestimated.
The eccentricity distribution as parameterized by Equation (A4) is relatively independent of binfrac, errfrac, and n. Figure 7 shows the range of allowed eccentricity distributions in black based on the synthetic populations from the MCMC analysis with 1σ uncertainties. The blue line shows the eccentricity distribution based on RV surveys (Wright et al. 2011) based on an analysis by Kipping (2013). A chi-square test of the two distributions gives χ2 = 46.2 for 10 samples, which indicates that the two distributions are different with high confidence. The multi-planet eccentricity distribution is more sharply peaked toward zero, thus highly eccentric planets in multi-planet systems are relativity rare compared to the RV sample. It would be interesting to further break down the RV sample to compare the eccentricity distributions of single and multi-planet samples, but it outside the scope of the initial analysis presented here.
Figure 7. Comparison of the eccentricity distribution of multi-planet systems (black), based on matching synthetic population models to the transiting multi-planet sample, to the eccentricity distribution of RV planets (blue; Wright et al. 2011; Kipping 2013). The distributions are different at high significance. The relative fraction of low eccentricity planets was found to be large for the transiting multi-planet population.
Download figure:
Standard image High-resolution imageThe difference in ρc when hierarchical blends are present will be strongly dependent on probability distribution of the underlying mass function. As shown in Figure 5, the distribution of dρc, i − j will contain an increased number of large mismatches as n decreases. In the case of equal or nearly equal-mass binaries the difference in ρc will be indistinguishable from measurement error. The synthetic population model shows that as errfrac is reduced, n pushes toward large values that produce a large number of equal mass binaries. For errfrac = 1, there is no strong measurement of the mass fraction distribution.
The value of ρc can also be directly compared to ρ⋆ (i.e., Tingley et al. 2011) from Table 2. Eccentric planets will be seen when the transit occurs close to pericenter, which decreases the transit duration relative to a circular orbit. This pushes ρc toward larger values (a denser star). As stated above, hierarchical blends will also produce a bias toward larger values of ρc as Kepler planet host stars are typically close to the main sequence. Measurement error should not introduce any bias. A comparison of ρ⋆ and ρc in a similar manor to the comparison of ρc for planet pairs could be carried out, but a strict requirement is that ρ⋆ is a good estimate of the true mean stellar density, which is not true. In particular, the use of Yale-Yonsei evolution models are known to produce radii too large and hence, densities that are systematically too small for low-mass stars (Plavchan et al. 2014). However, when ρ⋆ is based on asteroseismology (Huber et al. 2013), such biases are likely better controlled. Figure 8 shows the difference in ρ⋆ and ρc scaled by the uncertainty versus ρ⋆. The asteroseismology sample is limited to solar-like and evolved stars as the amplitudes of p-mode oscillations scale proportionally to stellar luminosity. The bias toward smaller values of ρ⋆ for low mass stars with ρ⋆ based on Yonsei-Yale models can be seen for stars with ρ⋆ >3 g cm−3. From this small sample, there is evidence of a bias of ρc being larger that ρ⋆; however, the sample is too small to draw any inferences on the underlying eccentricity and hierarchical blend population.
Figure 8. Difference between ρ⋆ and ρc is plotted vs. ρ⋆ to compare the mean stellar density derived from stellar modeling (ρ⋆) to the geometrical estimate of the mean stellar density when a circular orbit is assumed (ρc). Values of ρ⋆ derived from asteroseismology are shown with filled circles, values of ρ⋆ derived from fitting stellar evolution models to spectroscopic classification (see Section 3) are shown with dots. The population of dots shifts toward positive values of ρ⋆–ρc as ρ⋆ increases due to a poor match of our adopted stellar evolution models at low stellar mass.
Download figure:
Standard image High-resolution imageCurrently, the sample of multi-planet systems and well characterized stars is too small to draw strong conclusions about the number of hierarchical blends. The strongest constraints currently come from observations that attempt to directly detect a nearby companion that may be gravitationally bound. When additional multi-planet candidates and observations become available, model fits can be repeated to determine the rate of hierarchical blends and, in turn, establish whether orbital planes in binary stars systems are aligned.
APPENDIX B: LIST OF SYMBOLS AND ABBREVIATIONS
- 1.MJ, RJ–mass and radius in Jupiter units.
- 2.S/N—signal-to-noise ratio.
- 3.KOI—Kepler Object of Interest.
- 4.FA—false alarms are candidates with a S/N < 7.1.
- 5.FP—astrophysical false positive.
- 6.P/T0—period/T0 collision false positive produced by the instrument.
- 7.PC—planetary candidate.
- 8.EB—eclipsing binary.
- 9.BEB—background eclipsing binary.
- 10.S—ratio of incident flux relative to the Earth.
- 11.ρ⋆—mean stellar density.
- 12.ρc—transiting derived mean stellar density for circular orbits.
Footnotes
- 16
Based on the NASA Exoplanet Archive 2013/11/12.
- 17
Based on KOIs 1–3149 having an FP or PC status on 2013/08/01.
- 18
V3 retrieved 2013/04/24 http://keplerebs.villanova.edu/.
- 19
Observations labeled as SAP_FLUX from FITS files retrieved from The Barbara A. Mikulski Archive for Space Telescopes (MAST).
- 20
Observations with SAP_QUALITY=0 from FITS files retrieved from MAST.
- 21