A Habitable Zone Census via Transit Timing and the Imperative for Continuing to Observe the Kepler Field

We propose a scientific program to complete a census of planets, characterizing their masses, orbital properties, and dynamical histories using continued observations of the Kepler field of view with the Kepler spacecraft in a two reaction wheel mission. Even with a significantly reduced photometric precision, extending time-domain observations of this field is uniquely capable of pursuing several critical science goals: 1) measuring the architectures of planetary systems by identifying non-transiting planets interleaved among known transiting planets, 2) establishing the mass-radius relationship for planets in the important transition region between small, gas-rich sub-Neptune planets and large, rocky super-Earths, and 3) uncovering dynamical evidence of the formation and evolution of the inner regions of planetary systems. To meet these objectives, the unique multi-object observing capabilities of Kepler will be used in a set of concurrent campaigns with specific motivations. These campaigns focus largely on the ability to interpret Transit Timing Variations (TTVs) that result from dynamical interactions among planets in a system and include: 1) observations of systems that exhibit large TTVs and are particularly rich in dynamical information, 2) observations of systems where additional transit times will yield mass measurements of the constituent planets, 3) observations of systems where the TTV signal evolves over very long timescales, and 4) observations of systems with long-period planet candidates where additional transits will remove orbital period ambiguities caused by gaps in the original Kepler data.

On The Relative Sizes of Planets Within Kepler Multiple Candidate Systems

We present a study of the relative sizes of planets within the multiple candidate systems discovered with the $Kepler$ mission. We have compared the size of each planet to the size of every other planet within a given planetary system after correcting the sample for detection and geometric biases. We find that for planet-pairs for which one or both objects is approximately Neptune-sized or larger, the larger planet is most often the planet with the longer period. No such size--location correlation is seen for pairs of planets when both planets are smaller than Neptune. Specifically, if at least one planet in a planet-pair has a radius of $gtrsim 3R_oplus$, $68pm 6%$ of the planet pairs have the inner planet smaller than the outer planet, while no preferred sequential ordering of the planets is observed if both planets in a pair are smaller than $lesssim3 R_oplus$.

Transit Timing Observations from Kepler: IV. Confirmation of 4 Multiple Planet Systems by Simple Physical Models

Eighty planetary systems of two or more planets are known to orbit stars other than the Sun. For most, the data can be sufficiently explained by non-interacting Keplerian orbits, so the dynamical interactions of these systems have not been observed. Here we present 4 sets of lightcurves from the Kepler spacecraft, which each show multiple planets transiting the same star. Departure of the timing of these transits from strict periodicity indicates the planets are perturbing each other: the observed timing variations match the forcing frequency of the other planet. This confirms that these objects are in the same system. Next we limit their masses to the planetary regime by requiring the system remain stable for astronomical timescales. Finally, we report dynamical fits to the transit times, yielding possible values for the planets masses and eccentricities. As the timespan of timing data increases, dynamical fits may allow detailed constraints on the systems architectures, even in cases for which high-precision Doppler follow-up is impractical.

Transit Timing Observations from Kepler: II. Confirmation of Two Multiplanet Systems via a Non-parametric Correlation Analysis

We present a new method for confirming transiting planets based on the combination of transit timingn variations (TTVs) and dynamical stability. Correlated TTVs provide evidence that the pair of bodies are in the same physical system. Orbital stability provides upper limits for the masses of the transiting companions that are in the planetary regime. This paper describes a non-parametric technique for quantifying the statistical significance of TTVs based on the correlation of two TTV data sets. We apply this method to an analysis of the transit timing variations of two stars with multiple transiting planet candidates identified by Kepler. We confirm four transiting planets in two multiple planet systems based on their TTVs and the constraints imposed by dynamical stability. An additional three candidates in these same systems are not confirmed as planets, but are likely to be validated as real planets once further observations and analyses are possible. If all were confirmed, these systems would be near 4:6:9 and 2:4:6:9 period commensurabilities. Our results demonstrate that TTVs provide a powerful tool for confirming transiting planets, including low-mass planets and planets around faint stars for which Doppler follow-up is not practical with existing facilities. Continued Kepler observations will dramatically improve the constraints on the planet masses and orbits and provide sensitivity for detecting additional non-transiting planets. If Kepler observations were extended to eight years, then a similar analysis could likely confirm systems with multiple closely spaced, small transiting planets in or near the habitable zone of solar-type stars.

Transit Timing Observations from Kepler: VI. Transit Timing Variation Candidates in the First Seventeen Months from Polynomial Models

Transit timing variations provide a powerful tool for confirming and characterizing transiting planets, as well as detecting non-transiting planets. We report the results an updated TTV analysis for 1481 planet candidates (Borucki et al. 2011; Batalha et al. 2012) based on transit times measured during the first sixteen months of Kepler observations. We present 39 strong TTV candidates based on long-term trends (2.8% of suitable data sets). We present another 136 weaker TTV candidates (9.8% of suitable data sets) based on excess scatter of TTV measurements about a linear ephemeris. We anticipate that several of these planet candidates could be confirmed and perhaps characterized with more detailed TTV analyses using publicly available Kepler observations. For many others, Kepler has observed a long-term TTV trend, but an extended Kepler mission will be required to characterize the system via TTVs. We find that the occurrence rate of planet candidates that show TTVs is significantly increased (~68%) for planet candidates transiting stars with multiple transiting planet candidate when compared to planet candidates transiting stars with a single transiting planet candidate.

Transit Timing Observations from Kepler: VI. Potentially interesting candidate systems from Fourier-based statistical tests

We analyze the deviations of transit times from a linear ephemeris for the Kepler Objects of Interest (KOI) through Quarter six (Q6) of science data. We conduct two statistical tests for all KOIs and a related statistical test for all pairs of KOIs in multi-transiting systems. These tests identify several systems which show potentially interesting transit timing variations (TTVs). Strong TTV systems have been valuable for the confirmation of planets and their mass measurements. Many of the systems identified in this study should prove fruitful for detailed TTV studies.

The Distribution of Transit Durations for Kepler Planet Candidates and Implications for their Orbital Eccentricities

Doppler planet searches have discovered that giant planets follow orbits with a wide range of orbital eccentricities, revolutionizing theories of planet formation. The discovery of hundreds of exoplanet candidates by NASAs Kepler mission enables astronomers to characterize the eccentricity distribution of small exoplanets. Measuring the eccentricity of individual planets is only practical in favorable cases that are amenable to complementary techniques (e.g., radial velocities, transit timing variations, occultation photometry). Yet even in the absence of individual eccentricities, it is possible to study the distribution of eccentricities based on the distribution of transit durations (relative to the maximum transit duration for a circular orbit). We analyze the transit duration distribution of Kepler planet candidates. We find that for host stars with T_eff > 5100 K we cannot invert this to infer the eccentricity distribution at this time due to uncertainties and possible systematics in the host star densities. With this limitation in mind, we compare the observed transit duration distribution with models to rule out extreme distributions. If we assume a Rayleigh eccentricity distribution for Kepler planet candidates, then we find best-fits with a mean eccentricity of 0.1-0.25 for host stars with T_eff < 5100 K. We compare the transit duration distribution for different subsets of Kepler planet candidates and discuss tentative trends with planetary radius and multiplicity. High-precision spectroscopic follow-up observations for a large sample of host stars will be required to confirm which trends are real and which are the results of systematic errors in stellar radii. Finally, we identify planet candidates that must be eccentric or have a significantly underestimated stellar radius.

Transit Timing Observations from Kepler: I. Statistical Analysis of the First Four Months

The architectures of multiple planet systems can provide valuable constraints on models of planet formation, including orbital migration, and excitation of orbital eccentricities and inclinations. NASAs Kepler mission has identified 1235 transiting planet candidates (Borcuki et al 2011). The method of transit timing variations (TTVs) has already confirmed 7 planets in two planetary systems (Holman et al. 2010; Lissauer et al. 2011a). We perform a transit timing analysis of the Kepler planet candidates. We find that at least ~12% of planet candidates currently suitable for TTV analysis show evidence suggestive of TTVs, representing at least ~65 TTV candidates. In all cases, the time span of observations must increase for TTVs to provide strong constraints on planet masses and/or orbits, as expected based on n-body integrations of multiple transiting planet candidate systems (assuming circular and coplanar orbits). We find that the fraction of planet candidates showing TTVs in this data set does not vary significantly with the number of transiting planet candidates per star, suggesting significant mutual inclinations and that many stars with a single transiting planet should host additional non-transiting planets. We anticipate that Kepler could confirm (or reject) at least ~12 systems with multiple transiting planet candidates via TTVs. Thus, TTVs will provide a powerful tool for confirming transiting planets and characterizing the orbital dynamics of low-mass planets. If Kepler observations were extended to at least six years, then TTVs would provide much more precise constraints on the dynamics of systems with multiple transiting planets and would become sensitive to planets with orbital periods extending into the habitable zone of solar-type stars.

The California Planet Survey III. A Possible 2:1 Resonance in the Exoplanetary Triple System HD 37124

We present new radial velocities from Keck Observatory and both Newtonian and Keplerian solutions for the triple-planet system orbiting HD 37124. The orbital solution for this system has improved dramatically since the third planet was first reported in Vogt et al. 2005 with an ambiguous orbital period. We have resolved this ambiguity, and the outer two planets have an apparent period commensurability of 2:1. A dynamical analysis finds both resonant and non-resonant configurations consistent with the radial velocity data, and constrains the mutual inclinations of the planets to be less than about 30 degrees. We discuss HD 37124 in the context of the other 19 exoplanetary systems with apparent period commenserabilities, which we summarize in a table. We show that roughly one in three well-characterized multiplanet systems has a apparent low-order period commensuribility, which is more than would naively be expected if the periods of exoplanets in known multiplanet systems were drawn randomly from the observed distribution of planetary orbital periods.

The Formation Mechanism of Gas Giants on Wide Orbits

The recent discoveries of massive planets on ultra-wide orbits of HR 8799 (Marois et al. 2008) and Fomalhaut (Kalas et al. 2008) present a new challenge for planet formation theorists. Our goal is to figure out which of three giant planet formation mechanisms--core accretion (with or without migration), scattering from the inner disk, or gravitational instability--could be responsible for Fomalhaut b, HR 8799 b, c and d, and similar planets discovered in the future. This paper presents the results of numerical experiments comparing the long-period planet formation efficiency of each possible mechanism in model A star, G star and M star disks. First, a simple core accretion simulation shows that planet cores forming beyond 35 AU cannot reach critical mass, even under the most favorable conditions one can construct. Second, a set of N-body simulations demonstrates that planet-planet scattering does not create stable, wide-orbit systems such as HR 8799. Finally, a linear stability analysis verifies previous work showing that global spiral instabilities naturally arise in high-mass disks. We conclude that massive gas giants on stable orbits with semimajor axes greater than 35 AU form by gravitational instability in the disk. We recommend that observers examine the planet detection rate as a function of stellar age, controlling for planet dimming with time. If planet detection rate is found to be independent of stellar age, it would confirm our prediction that gravitational instability is the dominant mode of producing detectable planets on wide orbits. We also predict that the occurrence ratio of long-period to short-period gas giants should be highest for M dwarfs due to the inefficiency of core accretion and the expected small fragment mass in their disks.