New transiting planet candidates are identified in 16 months (2009 May-2010 of data from the Kepler spacecraft. Nearly 5000 periodic transit-like signals are vetted against astrophysical and instrumental false positives yielding 1108 viable new planet candidates, bringing the total count up to over 2300. Improved vetting metrics are employed, contributing to higher catalog reliability. Most notable is the noise-weighted robust averaging of multiquarter photo-center offsets derived from difference image analysis that identifies likely background eclipsing binaries. Twenty-two months of photometry are used for the purpose of characterizing each of the candidates. Ephemerides (transit epoch, T 0 , and orbital period, P) are tabulated as well as the products of light curve modeling: reduced radius (R P /R ), reduced semimajor axis (d/R ), and impact parameter (b). The largest fractional increases are seen for the smallest planet candidates (201% for candidates smaller than 2 R ⊕ compared to 53% for candidates larger than 2 R ⊕ ) and those at longer orbital periods (124% for candidates outside of 50 day orbits versus 86% for candidates inside of 50 day orbits). The gains are larger than expected from increasing the observing window from 13 months (Quarters 1-5) to 16 months (Quarters 1-6) even in regions of parameter space where one would have expected the previous catalogs to be complete. Analyses of planet frequencies based on previous catalogs will be affected by such incompleteness. The fraction of all planet candidate host stars with multiple candidates has grown from 17% to 20%, and the paucity of short-period giant planets in multiple systems is still evident. The progression 1The Astrophysical Journal Supplement Series, 204:24 (21pp), 2013 February Batalha et al. toward smaller planets at longer orbital periods with each new catalog release suggests that Earth-size planets in the habitable zone are forthcoming if, indeed, such planets are abundant.

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 massradius 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 firstorder 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• , for the packed systems of small planets probed by these observations.

We report the distribution of planets as a function of planet radius, orbital period, and stellar effective temperature for orbital periods less than 50 days around Solar-type (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 R ⊕ . For each of the 156,000 target stars we assess the detectability of planets as a function of planet radius, R p , and orbital period, P , using a measure of the detection efficiency for each star. We also correct for the geometric probability of transit, R ⋆ /a. We consider first Kepler target stars within the "solar subset" having T eff = 4100-6100 K, log g = 4.0-4.9, and Kepler magnitude Kp < 15 mag, i.e. bright, main sequence GK stars. We include only those stars having photometric noise low enough to permit detection of planets down to 2 R ⊕ . 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. The resulting occurrence of planets varies by more than three orders of magnitude in the radius-orbital period plane and increases substantially down to the smallest radius (2 R ⊕ ) and out to the longest orbital period (50 days, ∼0.25 AU) in our study. For P < 50 days, the distribution of planet radii is given by a power law, df /d log R = k R R α with k R = 2.9 +0.5 −0.4 , α = −1.92 ± 0.11, and R = R p /R ⊕ . This rapid increase in planet occurrence with decreasing planet size agrees with the prediction of core-accretion formation, but disagrees with population synthesis models that predict a desert at super-Earth and Neptune sizes for close-in orbits. Planets with orbital periods shorter than 2 days are extremely rare; for R p > 2 R ⊕ we measure an occurrence of less than 0.001 planets per star. For all planets with orbital periods less than 50 days, we measure occurrence of 0.130 ± 0.008, 0.023 ± 0.003, and 0.013 ± 0.002 planets per star for planets with radii 2-4, 4-8, and 8-32 R ⊕ , in agreement with Doppler surveys. 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 a broader stellar T eff range of 3600-7100 K, spanning M0 to F2 dwarfs. Over this range, the occurrence of 2-4 R ⊕ planets in the Kepler field linearly increases with decreasing T eff , making these small planets seven times more abundant around cool stars (3600-4100 K) than the hottest stars in our sample (6600-7100 K).