We infer the number of planets-per-star as a function of orbital period and planet size using $Kepler$ archival data products with updated stellar properties from the $Gaia$ Data Release 2. Using hierarchical Bayesian modeling and Hamiltonian Monte Carlo, we incorporate planet radius uncertainties into an inhomogeneous Poisson point process model. We demonstrate that this model captures the general features of the outcome of the planet formation process around GK stars, and provides an infrastructure to use the $Kepler$ results to constrain analytic planet distribution models. We report an increased mean and variance in the marginal posterior distributions for the number of planets per $GK$ star when including planet radius measurement uncertainties. We estimate the number of planets-per-$GK$ star between 0.75 and 2.5 $R_{\oplus}$ and 50 to 300 day orbital periods to have a $68\%$ credible interval of $0.49$ to $0.77$ and a posterior mean of $0.63$. This posterior has a smaller mean and a larger variance than the occurrence rate calculated in this work and in Burke et al. (2015) for the same parameter space using the $Q1-Q16$ (previous $Kepler$ planet candidate and stellar catalog), and a larger mean and variance than when using the $DR25$ (latest $Kepler$ planet candidate and stellar catalog). We find that the accuracy and precision of our hierarchical Bayesian model posterior distributions are less sensitive to the total number of planets in the sample, and more so on the characteristics of the catalog completeness and reliability and the span of the planet parameter space.

We present the Kepler Object of Interest (KOI) catalog of transiting exoplanets based on searching four years of Kepler time series photometry (Data Release 25, Q1-Q17). The catalog contains 8054 KOIs of which 4034 are planet candidates with periods between 0.25 and 632 days. Of these candidates, 219 are new and include two in multi-planet systems (KOI-82.06 and KOI-2926.05), and ten high-reliability, terrestrial-size, habitable zone candidates. This catalog was created using a tool called the Robovetter which automatically vets the DR25 Threshold Crossing Events (TCEs, Twicken et al. 2016). The Robovetter also vetted simulated data sets and measured how well it was able to separate TCEs caused by noise from those caused by low signal-to-noise transits. We discusses the Robovetter and the metrics it uses to sort TCEs. For orbital periods less than 100 days the Robovetter completeness (the fraction of simulated transits that are determined to be planet candidates) across all observed stars is greater than 85%. For the same period range, the catalog reliability (the fraction of candidates that are not due to instrumental or stellar noise) is greater than 98%. However, for low signal-to-noise candidates between 200 and 500 days around FGK dwarf stars, the Robovetter is 76.7% complete and the catalog is 50.5% reliable. The KOI catalog, the transit fits and all of the simulated data used to characterize this catalog are available at the NASA Exoplanet Archive.

We characterize the eccentricity distribution of a sample of ~50 short-period planet candidates using transit and occultation measurements from NASA's Kepler Mission. First, we evaluate the sensitivity of our hierarchical Bayesian modeling and test its robustness to model misspecification using simulated data. When analyzing actual data assuming a Rayleigh distribution for eccentricity, we find that the posterior mode for the dispersion parameter is $\sigma=0.081 \pm^{0.014}_{0.003}$. We find that a two-component Gaussian mixture model for $e \cos \omega$ and $e \sin \omega$ provides a better model than either a Rayleigh or Beta distribution. Based on our favored model, we find that $\sim90\%$ of planet candidates in our sample come from a population with an eccentricity distribution characterized by a small dispersion ($\sim0.01$), and $\sim10\%$ come from a population with a larger dispersion ($\sim0.22$). Finally, we investigate how the eccentricity distribution correlates with selected planet and host star parameters. We find evidence that suggests systems around higher metallicity stars and planet candidates with smaller radii come from a more complex eccentricity distribution.

The meteorology of hot Jupiters has been characterized primarily with thermal measurements, but recent observations suggest the possibility of directly detecting the winds by observing the Doppler shift of spectral lines seen during transit. Motivated by these observations, we show how Doppler measurements can place powerful constraints on the meteorology. We show that the atmospheric circulation--and Doppler signature--of hot Jupiters splits into two regimes. Under weak stellar insolation, the day-night thermal forcing generates fast zonal jet streams from the interaction of atmospheric waves with the mean flow. In this regime, air along the terminator (as seen during transit) flows toward Earth in some regions and away from Earth in others, leading to a Doppler signature exhibiting superposed blueshifted and redshifted components. Under intense stellar insolation, however, the strong thermal forcing damps these planetary-scale waves, inhibiting their ability to generate jets. Strong frictional drag likewise damps these waves and inhibits jet formation. As a result, this second regime exhibits a circulation dominated by high-altitude, day-to-night airflow, leading to a predominantly blueshifted Doppler signature during transit. We present state-of-the-art circulation models including nongray radiative transfer to quantify this regime shift and the resulting Doppler signatures; these models suggest that cool planets like GJ 436b lie in the first regime, HD 189733b is transitional, while planets hotter than HD 209458b lie in the second regime. Moreover, we show how the amplitude of the Doppler shifts constrains the strength of frictional drag in the upper atmospheres of hot Jupiters. If due to winds, the ~2-km/sec blueshift inferred on HD 209458b may require drag time constants as short as 10^4-10^6 seconds, possibly the result of Lorentz-force braking on this planet's hot dayside.

Using radiation hydrodynamics simulations, we explore the evolution of circumplanetary disks around wide-orbit proto-gas giants. At large distances from the star (~100 AU), gravitational instability followed by disk fragmentation can form low-mass substellar companions (massive gas giants and/or brown dwarfs) that are likely to host large disks. We examine the initial evolution of these subdisks and their role in regulating the growth of their substellar companions, as well as explore consequences of their interactions with circumstellar material. We find that subdisks that form in the context of GIs evolve quickly from a very massive state. Long-term accretion rates from the subdisk onto the proto-gas giant reach ~0.3 Jupiter masses per kyr. We also find consistency with previous simulations, demonstrating that subdisks are truncated at ~1/3 of the companion's Hill radius and are thick, with (h/r) of great than or equal to 0.2. The thickness of subdisks draws to question the use of thin-disk approximations for understanding the behavior of subdisks, and the morphology of subdisks has implications for the formation and extent of satellite systems. These subdisks create heating events in otherwise cold regions of the circumstellar disk, and serve as planet formation beacons that can be detected by instruments such as ALMA.

The dynamical evolution of planetary systems leaves observable signatures in debris disks. Optical images trace micron-sized grains, which are strongly affected by stellar radiation and need not coincide with their parent body population. Observations of mm-size grains accurately trace parent bodies, but previous images lack the resolution and sensitivity needed to characterize the ring's morphology. Here we present ALMA 350 GHz observations of the Fomalhaut debris ring. These observations demonstrate that the parent body population is 13-19 AU wide with a sharp inner and outer boundary. We discuss three possible origins for the ring, and suggest that debris confined by shepherd planets is the most consistent with the ring's morphology.

We report observations of HD 80606 using the 10.4-m Gran Telescopio Canarias (GTC) and the OSIRIS tunable filter imager. We acquired very-high-precision, narrow-band photometry in four bandpasses around the K I absorption feature during the January 2010 transit of HD 80606b and during out-of-transit observations conducted in January and April of 2010. We obtained differential photometric precisions of \sim 2.08e-4 for the in-transit flux ratio measured at 769.91-nm, which probes the K I line core. We find no significant difference in the in-transit flux ratio between observations at 768.76 and 769.91 nm. Yet, we find a difference of \sim 8.09 \pm 2.88e-4 between these observations and observations at a longer wavelength that probes the K I wing (777.36 nm). While the presence of red noise in the transit data has a non-negligible effect on the uncertainties in the flux ratio, the 777.36-769.91 nm colour during transit shows no effects from red noise and also indicates a significant colour change, with a mean value of \sim 8.99\pm0.62e-4. This large change in the colour is equivalent to a \sim 4.2% change in the apparent planetary radius with wavelength, which is much larger than the atmospheric scale height. This implies the observations probed the atmosphere at very low pressures as well as a dramatic change in the pressure at which the slant optical depth reaches unity between \sim770 and 777 nm. We hypothesize that the excess absorption may be due to K I in a high-speed wind being driven from the exoplanet's exosphere. We discuss the viability of this and alternative interpretations, including stellar limb darkening, starspots, and effects from Earth's atmosphere. We strongly encourage follow-up observations of HD 80606b to confirm the signal measured here. Finally, we discuss the future prospects for exoplanet characterization using tunable filter spectrophotometry.

We have developed 1-D photochemical and thermochemical kinetics and diffusion models for the transiting exoplanets HD 189733b and HD 209458b to study the effects of disequilibrium chemistry on the atmospheric composition of "hot Jupiters." Here we investigate the coupled chemistry of neutral carbon, hydrogen, oxygen, and nitrogen species, and we compare the model results with existing transit and eclipse observations. We find that the vertical profiles of molecular constituents are significantly affected by transport-induced quenching and photochemistry, particularly on cooler HD 189733b; however, the warmer stratospheric temperatures on HD 209458b can help maintain thermochemical equilibrium and reduce the effects of disequilibrium chemistry. For both planets, the methane and ammonia mole fractions are found to be enhanced over their equilibrium values at pressures of a few bar to less than a mbar due to transport-induced quenching, but CH4 and NH3 are photochemically removed at higher altitudes. Atomic species, unsaturated hydrocarbons (particularly C2H2), some nitriles (particularly HCN), and radicals like OH, CH3, and NH2 are enhanced overequilibrium predictions because of quenching and photochemistry. In contrast, CO, H2O, N2, and CO2 more closely follow their equilibrium profiles, except at pressures < 1 microbar, where CO, H2O, and N2 are photochemically destroyed and CO2 is produced before its eventual high-altitude destruction. The enhanced abundances of HCN, CH4, and NH3 in particular are expected to affect the spectral signatures and thermal profiles HD 189733b and other, relatively cool, close-in transiting exoplanets. We examine the sensitivity of our results to the assumed temperature structure and eddy diffusion coefficientss and discuss further observational consequences of these models.

We present Gran Telescopio Canarias (GTC) optical transit narrow-band photometry of the hot-Jupiter exoplanet XO-2b using the OSIRIS instrument. This unique instrument has the capabilities to deliver high cadence narrow-band photometric lightcurves, allowing us to probe the atmospheric composition of hot Jupiters from the ground. The observations were taken during three transit events which cover four wavelengths at spectral resolutions near 500, necessary for observing atmospheric features, and have near-photon limited sub-mmag precisions. Precision narrow-band photometry on a large aperture telescope allows for atmospheric transmission spectral features to be observed for exoplanets around much fainter stars than those of the well studied targets HD209458b and HD189733b, providing access to the majority of known transiting planets. For XO-2b, we measure planet-to-star radius contrasts of R_pl/R_star=0.10508+/-0.00052 at 6792 Ang, 0.10640+/-0.00058 at 7582 Ang, and 0.10686+/-0.00060 at 7664.9 Ang, and 0.10362+/-0.00051 at 8839 Ang. These measurements reveal significant spectral features at two wavelengths, with an absorption level of 0.067+/-0.016% at 7664.9 Ang due to atmospheric potassium in the line core (a 4.1-sigma significance level), and an absorption level of 0.058+/-0.016% at 7582 Ang, (a 3.6-sigma significance level). When comparing our measurements to hot-Jupiter atmospheric models, we find good agreement with models which are dominated in the optical by alkali metals. This is the first evidence for potassium in an extrasolar planet, an element that has long been theorized along with sodium to be a dominant source of opacity at optical wavelengths for hot Jupiters.

We explore the transmission spectrum of the Neptune-class exoplanet GJ 436b, including the possibility that its atmospheric opacity is dominated by a variety of non- equilibrium chemical products. We also validate our transmission code by demonstrating close agreement with analytic models that use only Rayleigh scattering or water vapor opacity. We find broad disagreement with radius variations predicted by another published model. For GJ 436b, the relative coolness of the planet's atmosphere, along with its implied high metallicity, may make it dissimilar in character compared to "hot Jupiters." Some recent observational and modeling efforts suggest low relative abundances of H2O and CH4 present in GJ 436b's atmosphere, compared to calculations from equilibrium chemistry. We include these characteristics in our models and examine the effects of absorption from methane-derived higher order hydrocarbons. Significant absorption from HCN and C2H2 are found throughout the infrared, while C2H4 and C2H6 are less easily seen. We perform detailed simulations of JWST observations, including all likely noise sources, and find that we will be able to constrain chemical abundance regimes from this planet's transmission spectrum. For instance, the width of the features at 1.5, 3.3, and 7 microns indicates the amount of HCN versus C2H2 present. The NIRSpec prism mode will be useful due to its large spectral range and the relatively large number of photo-electrons recorded per spectral resolution element. However, extremely bright host stars like GJ 436 may be better observed with a higher spectroscopic resolution mode in order to avoid detector saturation. We find that observations with the MIRI low resolution spectrograph should also have high signal-to-noise in the 5 - 10 micron range due to the brightness of the star and the relatively low spectral resolution (R ~ 100) of this mode.

We compute models of the transmission spectra of planets HD 209458b, HD 189733b, and generic hot Jupiters. We examine the effects of temperature, surface gravity, and metallicity for the generic planets as a guide to understanding transmission spectra in general. We find that carbon dioxide absorption at 4.4 and 15 microns is prominent at high metallicity, and is a clear metallicity indicator. For HD 209458b and HD 189733b, we compute spectra for both one-dimensional and three-dimensional model atmospheres and examine the differences between them. The differences are usually small, but can be large if atmospheric temperatures are near important chemical abundance boundaries. The calculations for the 3D atmospheres, and their comparison with data, serve as constraints on these dynamical models that complement the secondary eclipse and light curve data sets. For HD 209458b, even if TiO and VO gases are abundant on the day side, their abundances can be considerably reduced on the cooler planetary limb. However, given the predicted limb temperatures and TiO abundances, the model's optical opacity is too high. For HD 189733b we find a good match with some infrared data sets and constrain the altitude of a postulated haze layer. For this planet, substantial differences can exist between the transmission spectra of the leading and trailing hemispheres, which is an excellent probe of carbon chemistry. In thermochemical equilibrium, the cooler leading hemisphere is methane-dominated, and the hotter trailing hemisphere is CO-dominated, but these differences may be eliminated by non-equilibrium chemistry due to vertical mixing. It may be possible to constrain the carbon chemistry of this planet, and its spatial variation, with JWST.