We present a new framework to characterize the occurrence rates of planet candidates identified by Kepler based on hierarchical Bayesian modeling, Approximate Bayesian Computing (ABC), and sequential importance sampling. For this study we adopt a simple 2-D grid in planet radius and orbital period as our model and apply our algorithm to estimate occurrence rates for Q1-Q16 planet candidates orbiting around solar-type stars. We arrive at significantly increased planet occurrence rates for small planet candidates ($R_p<1.25 R_{\oplus}$) at larger orbital periods ($P>80$d) compared to the rates estimated by the more common inverse detection efficiency method. Our improved methodology estimates that the occurrence rate density of small planet candidates in the habitable zone of solar-type stars is $1.6^{+1.2}_{-0.5}$ per factor of 2 in planet radius and orbital period. Additionally, we observe a local minimum in the occurrence rate for strong planet candidates marginalized over orbital period between 1.5 and 2$R_{\oplus}$ that is consistent with previous studies. For future improvements, the forward modeling approach of ABC is ideally suited to incorporating multiple populations, such as planets, astrophysical false positives and pipeline false alarms, to provide accurate planet occurrence rates and uncertainties. Furthermore, ABC provides a practical statistical framework for answering complex questions (e.g., frequency of different planetary architectures) and providing sound uncertainties, even in the face of complex selection effects, observational biases, and follow-up strategies. In summary, ABC offers a powerful tool for accurately characterizing a wide variety of astrophysical populations.

We present new multi-band transit photometry of three small (R$_{p}$ < 6 R$_{\oplus}$), short-period (P < 6 days) Kepler planet candidates acquired with the Gran Telescopio Canarias. These observations supplement the results presented in Col\'on & Ford (2011) and Col\'on et al. (2012), where we used multicolor transit photometry of five Kepler planet candidates to search for wavelength-dependent transit depths and either validate planet candidates or identify eclipsing binary false positives within our sample. In those previous studies, we provided evidence that three targets were false positives and two targets were planets. Here, we present observations that provide evidence supporting a planetary nature for KOI 439.01 and KOI 732.01, and we find that KOI 531.01, a 6 R$_{\oplus}$ planet candidate around an M dwarf, is likely a false positive. We also present a discussion of the purported "sub-Jovian desert" in the orbital period-planet radius plane, which cannot be easily explained by observational bias. Both KOI 439.01 and KOI 732.01 are likely planets located within the so-called desert and should be investigated with further follow-up observations. As only ~30 of the ~3600 currently active Kepler planet candidates are located within the sub-Jovian desert, it will be interesting to see if these candidates also survive the vetting process and fill in the gap in the period-radius plane. Confirming planets in this regime will be important for understanding planetary migration and evolution processes, and we urge additional follow-up observations of these planet candidates to confirm their nature.

We report on the orbital architectures of Kepler systems having multiple planet candidates identified in the analysis of data from the first six quarters of Kepler data and reported by Batalha et al. (2013). These data show 899 transiting planet candidates in 365 multiple-planet systems and provide a powerful means to study the statistical properties of planetary systems. Using a generic mass-radius relationship, we find that only two pairs of planets in these candidate systems (out of 761 pairs total) appear to be on Hill-unstable orbits, indicating ~96% of the candidate planetary systems are correctly interpreted as true systems. We find that planet pairs show little statistical preference to be near mean-motion resonances. We identify an asymmetry in the distribution of period ratios near first-order resonances (e.g., 2:1, 3:2), with an excess of planet pairs lying wide of resonance and relatively few lying narrow of resonance. Finally, based upon the transit duration ratios of adjacent planets in each system, we find that the interior planet tends to have a smaller transit impact parameter than the exterior planet does. This finding suggests that the mode of the mutual inclinations of planetary orbital planes is in the range 1.0-2.2 degrees, for the packed systems of small planets probed by these observations.

We provide updates to the Kepler planet candidate sample based upon nearly two years of high-precision photometry (i.e., Q1-Q8). From an initial list of nearly 13,400 Threshold Crossing Events (TCEs), 480 new host stars are identified from their flux time series as consistent with hosting transiting planets. Potential transit signals are subjected to further analysis using the pixel-level data, which allows background eclipsing binaries to be identified through small image position shifts during transit. We also re-evaluate Kepler Objects of Interest (KOI) 1-1609, which were identified early in the mission, using substantially more data to test for background false positives and to find additional multiple systems. Combining the new and previous KOI samples, we provide updated parameters for 2,738 Kepler planet candidates distributed across 2,017 host stars. From the combined Kepler planet candidates, 472 are new from the Q1-Q8 data examined in this study. The new Kepler planet candidates represent ~40% of the sample with Rp~1 Rearth and represent ~40% of the low equilibrium temperature (Teq<300 K) sample. We review the known biases in the current sample of Kepler planet candidates relevant to evaluating planet population statistics with the current Kepler planet candidate sample.

We confirm 27 planets in 13 planetary systems by showing the existence of statistically significant anti-correlated transit timing variations (TTVs), which demonstrates that the planet candidates are in the same system, and long-term dynamical stability, which places limits on the masses of the candidates---showing that they are planetary. %This overall method of planet confirmation was first applied to \kepler systems 23 through 32. All of these newly confirmed planetary systems have orbital periods that place them near first-order mean motion resonances (MMRs), including 6 systems near the 2:1 MMR, 5 near 3:2, and one each near 4:3, 5:4, and 6:5. In addition, several unconfirmed planet candidates exist in some systems (that cannot be confirmed with this method at this time). A few of these candidates would also be near first order MMRs with either the confirmed planets or with other candidates. One system of particular interest, Kepler-56 (KOI-1241), is a pair of planets orbiting a 12th magnitude, giant star with radius over three times that of the Sun and effective temperature of 4900 K---among the largest stars known to host a transiting exoplanetary system.

Using the OSIRIS instrument installed on the 10.4-m Gran Telescopio Canarias (GTC) we acquired multi-color transit photometry of four small (Rp < 5 R_Earth) short-period (P < 6 days) planet candidates recently identified by the Kepler space mission. These observations are part of a program to constrain the false positive rate for small, short-period Kepler planet candidates. Since planetary transits should be largely achromatic when observed at different wavelengths (excluding the small color changes due to stellar limb darkening), we use the observed transit color to identify candidates as either false positives (e.g., a blend with a stellar eclipsing binary either in the background/foreground or bound to the target star) or validated planets. Our results include the identification of KOI 225.01 and KOI 1187.01 as false positives and the tentative validation of KOI 420.01 and KOI 526.01 as planets. The probability of identifying two false positives out of a sample of four targets is less than 1%, assuming an overall false positive rate for Kepler planet candidates of 10% (as estimated by Morton & Johnson 2011). Therefore, these results suggest a higher false positive rate for the small, short-period Kepler planet candidates than has been theoretically predicted by other studies which consider the Kepler planet candidate sample as a whole. Furthermore, our results are consistent with a recent Doppler study of short-period giant Kepler planet candidates (Santerne et al. 2012). We also investigate how the false positive rate for our sample varies with different planetary and stellar properties. Our results suggest that the false positive rate varies significantly with orbital period and is largest at the shortest orbital periods (P < 3 days), where there is a corresponding rise in the number of detached eclipsing binary stars... (truncated)

We present a statistical analysis that demonstrates that the overwhelming majority of Kepler candidate multiple transiting systems (multis) indeed represent true, physically-associated transiting planets. Binary stars provide the primary source of false positives among Kepler planet candidates, implying that false positives should be nearly randomly-distributed among Kepler targets. In contrast, true transiting planets would appear clustered around a smaller number of Kepler targets if detectable planets tend to come in systems and/or if the orbital planes of planets encircling the same star are correlated. There are more than one hundred times as many Kepler planet candidates in multi-candidate systems as would be predicted from a random distribution of candidates, implying that the vast majority are true planets. Most of these multis are multiple planet systems orbiting the Kepler target star, but there are likely cases where (a) the planetary system orbits a fainter star, and the planets are thus significantly larger than has been estimated, or (b) the planets orbit different stars within a binary/multiple star system. We use the low overall false positive rate among Kepler multis, together with analysis of Kepler spacecraft and ground-based data, to validate the closely-packed Kepler-33 planetary system, which orbits a star that has evolved somewhat off of the main sequence. Kepler-33 hosts five transiting planets with periods ranging from 5.67 to 41 days.

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.

About one-third of the ~1200 transiting planet candidates detected in the first four months of \ik data are members of multiple candidate systems. There are 115 target stars with two candidate transiting planets, 45 with three, 8 with four, and one each with five and six. We characterize the dynamical properties of these candidate multi-planet systems. The distribution of observed period ratios shows that the vast majority of candidate pairs are neither in nor near low-order mean motion resonances. Nonetheless, there are small but statistically significant excesses of candidate pairs both in resonance and spaced slightly too far apart to be in resonance, particularly near the 2:1 resonance. We find that virtually all candidate systems are stable, as tested by numerical integrations that assume a nominal mass-radius relationship. Several considerations strongly suggest that the vast majority of these multi-candidate systems are true planetary systems. Using the observed multiplicity frequencies, we find that a single population of planetary systems that matches the higher multiplicities underpredicts the number of singly-transiting systems. We provide constraints on the true multiplicity and mutual inclination distribution of the multi-candidate systems, revealing a population of systems with multiple super-Earth-size and Neptune-size planets with low to moderate mutual inclinations.

We report the distribution of planets as a function of planet radius (R_p), orbital period (P), and stellar effective temperature (Teff) for P < 50 day orbits around GK stars. These results are based on the 1,235 planets (formally "planet candidates") from the Kepler mission that include a nearly complete set of detected planets as small as 2 Earth radii (Re). For each of the 156,000 target stars we assess the detectability of planets as a function of R_p and P. We also correct for the geometric probability of transit, R*/a. We consider first stars within the "solar subset" having Teff = 4100-6100 K, logg = 4.0-4.9, and Kepler magnitude Kp < 15 mag. We include only those stars having noise low enough to permit detection of planets down to 2 Re. We count planets in small domains of R_p and P and divide by the included target stars to calculate planet occurrence in each domain. Occurrence of planets varies by more than three orders of magnitude and increases substantially down to the smallest radius (2 Re) and out to the longest orbital period (50 days, ~0.25 AU) in our study. For P < 50 days, the radius distribution is given by a power law, df/dlogR= k R^\alpha. This rapid increase in planet occurrence with decreasing planet size agrees with core-accretion, but disagrees with population synthesis models. We fit occurrence as a function of P to a power law model with an exponential cutoff below a critical period P_0. For smaller planets, P_0 has larger values, suggesting that the "parking distance" for migrating planets moves outward with decreasing planet size. We also measured planet occurrence over Teff = 3600-7100 K, spanning M0 to F2 dwarfs. The occurrence of 2-4 Re planets in the Kepler field increases with decreasing Teff, making these small planets seven times more abundant around cool stars than the hottest stars in our sample. [abridged]

NASA's Kepler Mission uses transit photometry to determine the frequency of earth-size planets in or near the habitable zone of Sun-like stars. The mission reached a milestone toward meeting that goal: the discovery of its first rocky planet, Kepler-10b. Two distinct sets of transit events were detected: 1) a 152 +/- 4 ppm dimming lasting 1.811 +/- 0.024 hours with ephemeris T[BJD]=2454964.57375+N*0.837495 days and 2) a 376 +/- 9 ppm dimming lasting 6.86 +/- 0.07 hours with ephemeris T[BJD]=2454971.6761+N*45.29485 days. Statistical tests on the photometric and pixel flux time series established the viability of the planet candidates triggering ground-based follow-up observations. Forty precision Doppler measurements were used to confirm that the short-period transit event is due to a planetary companion. The parent star is bright enough for asteroseismic analysis. Photometry was collected at 1-minute cadence for >4 months from which we detected 19 distinct pulsation frequencies. Modeling the frequencies resulted in precise knowledge of the fundamental stellar properties. Kepler-10 is a relatively old (11.9 +/- 4.5 Gyr) but otherwise Sun-like Main Sequence star with Teff=5627 +/- 44 K, Mstar=0.895 +/- 0.060 Msun, and Rstar=1.056 +/- 0.021 Rsun. Physical models simultaneously fit to the transit light curves and the precision Doppler measurements yielded tight constraints on the properties of Kepler-10b that speak to its rocky composition: Mpl=4.56 +/- 1.29 Mearth, Rpl=1.416 +/- 0.036 Rearth, and density=8.8 +/- 2.9 gcc. Kepler-10b is the smallest transiting exoplanet discovered to date.

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 NASA's 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.

When an extrasolar planet passes in front of its star (transits), its radius can be measured from the decrease in starlight and its orbital period from the time between transits. Multiple planets transiting the same star reveal more: period ratios determine stability and dynamics, mutual gravitational interactions reflect planet masses and orbital shapes, and the fraction of transiting planets observed as multiples has implications for the planarity of planetary systems. But few stars have more than one known transiting planet, and none has more than three. Here we report Kepler spacecraft observations of a single Sun-like star that reveal six transiting planets, five with orbital periods between 10 and 47 days plus a sixth one with a longer period. The five inner planets are among the smallest whose masses and sizes have both been measured, and these measurements imply substantial envelopes of light gases. The degree of coplanarity and proximity of the planetary orbits imply energy dissipation near the end of planet formation.

We present and discuss five candidate exoplanetary systems identified with the Kepler spacecraft. These five systems show transits from multiple exoplanet candidates. Should these objects prove to be planetary in nature, then these five systems open new opportunities for the field of exoplanets and provide new insights into the formation and dynamical evolution of planetary systems. We discuss the methods used to identify multiple transiting objects from the Kepler photometry as well as the false-positive rejection methods that have been applied to these data. One system shows transits from three distinct objects while the remaining four systems show transits from two objects. Three systems have planet candidates that are near mean motion commensurabilities---two near 2:1 and one just outside 5:2. We discuss the implications that multitransiting systems have on the distribution of orbital inclinations in planetary systems, and hence their dynamical histories; as well as their likely masses and chemical compositions. A Monte Carlo study indicates that, with additional data, most of these systems should exhibit detectable transit timing variations (TTV) due to gravitational interactions---though none are apparent in these data. We also discuss new challenges that arise in TTV analyses due to the presence of more than two planets in a system.

We present ground-based optical observations of the September 2009 and January 2010 transits of HD 80606b. Based on 3 partial light curves of the September 2009 event, we derive a midtransit time of T_c [HJD] = 2455099.196 +- 0.026, which is about 1 sigma away from the previously predicted time. We observed the January 2010 event from 9 different locations, with most phases of the transit being observed by at least 3 different teams. We determine a midtransit time of T_c [HJD] = 2455210.6502 +- 0.0064, which is within 1.3 sigma of the time derived from a Spitzer observation of the same event.

We report the discovery of a low-mass companion orbiting the metal-rich, main sequence F star TYC 2949-00557-1 during the MARVELS (Multi-object APO Radial Velocity Exoplanet Large-area Survey) Pilot Project. The host star has an effective temperature T_eff = 6135 +/- 40 K, log(g) = 4.4 +/- 0.1 and [Fe/H] = 0.32 +/- 0.01, indicating a mass of M = 1.25 +/- 0.09 M_\odot and R = 1.15 +/- 0.15 R_\odot. The companion has an orbital period of 5.69449 +/- 0.00023 days and straddles the hydrogen burning limit with a minimum mass of 64 M_J, and may thus be an example of the rare class of brown dwarfs orbiting at distances comparable to those of "Hot Jupiters." We present relative photometry that demonstrates the host star is photometrically stable at the few millimagnitude level on time scales of hours to years, and rules out transits for a companion of radius greater than 0.8 R_J at the 95% confidence level. Tidal analysis of the system suggests that the star and companion are likely in a double synchronous state where both rotational and orbital synchronization have been achieved. This is the first low-mass companion detected with a multi-object, dispersed, fixed-delay interferometer.