The Kepler DR25 planet candidate catalog was produced using an automated method of planet candidate identification based on various tests. These tests were tuned to obtain a reasonable but arbitrary balance between catalog completeness and reliability. We produce new catalogs with differing balances of completeness and reliability by varying these tests, and study the impact of these alternative catalogs on occurrence rates. We find that if there is no correction for reliability, different catalogs give statistically inconsistent occurrence rates, while if we correct for both completeness and reliability, we get statistically consistent occurrence rates. This is a strong indication that correction for completeness and reliability is critical for the accurate computation of occurrence rates. Additionally, we find that this result is the same whether using Bayesian Poisson likelihood MCMC or Approximate Bayesian Computation methods. We also examine the use of a Robovetter disposition score cut as an alternative to reliability correction, and find that while a score cut does increase the reliability of the catalog, it is not as accurate as performing a full reliability correction. We get the same result when performing a reliability correction with and without a score cut. Therefore removing low-score planets removes data without providing any advantage, and should be avoided when possible. We make our alternative catalogs publicly available, and propose that these should be used as a test of occurrence rate methods, with the requirement that a method should provide statistically consistent occurrence rates for all these catalogs.

We present Kepler exoplanet occurrence rates inferred with approximate Bayesian computation (ABC). By using the same planet catalogue, stellar sample, and characterization of completeness and reliability as Bryson et al. (2020), we are able to provide the first direct comparison of results from ABC to those derived with the popular Poisson-likelihood method. For planets with orbital periods between 50 and 400 days and radii between 0.75 and 2.5 $R_{\oplus}$, we find an integrated occurrence rate $F_{0} = 0.596_{-0.099}^{+0.092}$ planets per GK dwarf star. After correcting for reliability against astrophysical false positives and false alarms, we find $F_{0} = 0.421_{-0.072}^{+0.086}$. Our findings agree within 1$\sigma$ of Bryson et al. (2020), indicating that the results are robust and not method-dependent.

We present exoplanet occurrence rates estimated with approximate Bayesian computation for planets with radii between 0.5 and 16 $R_{\bigoplus}$ and orbital periods between 0.78 and 400 days, orbiting FGK dwarf stars. We base our results on an independent planet catalogue compiled from our search of all ~200,000 stars observed over the Kepler mission, with precise planetary radii supplemented by Gaia DR2-incorporated stellar radii. We take into account detection and vetting efficiency, planet radius uncertainty, and reliability against transit-like noise signals in the data. By analyzing our FGK occurrence rates as well as those computed after separating F-, G-, and K-type stars, we explore dependencies on stellar effective temperature, planet radius, and orbital period. We reveal new characteristics of the photoevaporation-driven "radius gap" between ~1.5 and 2 $R_{\bigoplus}$, indicating that the bimodal distribution previously revealed for $P$ < 100 days exists only over a much narrower range of orbital periods, above which sub-Neptunes dominate and below which super-Earths dominate. Finally, we provide several estimates of the "eta-Earth" value -- the frequency of potentially habitable, rocky planets orbiting Sun-like stars. For planets with sizes 0.75 - 1.5 $R_{\bigoplus}$ orbiting in a conservatively defined habitable zone (0.99 - 1.70 AU) around G-type stars, we place an upper limit (84.1th percentile) of <0.18 planets per star.

We present the results of an independent search of all ~200,000 stars observed over the four year Kepler mission (Q1-Q17) for multiplanet systems, using a three-transit minimum detection criteria to search orbital periods up to hundreds of days. We incorporate both automated and manual triage, and provide estimates of the completeness and reliability of our vetting pipeline. Our search returned 17 planet candidates (PCs) in addition to thousands of known Kepler Objects of Interest (KOIs), with a 98.8% recovery rate of already confirmed planets. We highlight the discovery of one candidate, KIC-7340288 b, that is both rocky (radius < 1.6 $R_{\bigoplus}$) and in the Habitable Zone (insolation between 0.25 and 2.2 times the Earth's insolation). Another candidate is an addition to the already known KOI-4509 system. We also present adaptive optics imaging follow-up for six of our new PCs, two of which reveal a line-of-sight stellar companion within 4''.

We study the geometry of the AGN obscurer in IRAS 09104+4109, an IR-luminous, radio-intermediate FR-I source at $z=0.442$, using infrared data from Spitzer and Herschel, X-ray data from Nustar, Swift, Suzaku, and Chandra, and an optical spectrum from Palomar. The infrared data imply a total rest-frame 1-1000$\mu$m luminosity of $5.5\times10^{46}$erg s$^{-1}$ and require both an AGN torus and starburst model. The AGN torus has an anisotropy-corrected IR luminosity of $4.9\times10^{46}$erg s$^{-1}$, and a viewing angle and half opening angle both of approximately $36$ degrees from pole-on. The starburst has a star formation rate of $(110\pm34)$M$_{\odot}$ yr$^{-1}$ and an age of $<50$Myr. These results are consistent with two epochs of luminous activity in 09104: one approximately $150$Myr ago, and one ongoing. The X-ray data suggest a photon index of $\Gamma \simeq 1.8$ and a line-of-sight column of $N_{\rm H} \simeq 5\times10^{23}$cm$^{-2}$. This argues against a reflection-dominated hard X-ray spectrum, which would have implied a much higher $N_{\rm H}$ and luminosity. The X-ray and infrared data are consistent with a bolometric AGN luminosity of $L_{\rm bol}\sim(0.5-2.5)\times10^{47}$erg s$^{-1}$. The X-ray and infrared data are further consistent with coaligned AGN obscurers in which the line of sight "skims" the torus. This is also consistent with the optical spectra, which show both coronal iron lines and broad lines in polarized but not direct light. Combining constraints from the X-ray, optical, and infrared data suggests that the AGN obscurer is within a vertical height of $20$pc, and a radius of $125$pc, of the nucleus.