AU Microscopii (AU Mic) is the second closest pre main sequence star, at a distance of 9.79 parsecs and with an age of 22 million years. AU Mic possesses a relatively rare and spatially resolved3 edge-on debris disk extending from about 35 to 210 astronomical units from the star, and with clumps exhibiting non-Keplerian motion. Detection of newly formed planets around such a star is challenged by the presence of spots, plage, flares and other manifestations of magnetic activity on the star. Here we report observations of a planet transiting AU Mic. The transiting planet, AU Mic b, has an orbital period of 8.46 days, an orbital distance of 0.07 astronomical units, a radius of 0.4 Jupiter radii, and a mass of less than 0.18 Jupiter masses at 3 sigma confidence. Our observations of a planet co-existing with a debris disk offer the opportunity to test the predictions of current models of planet formation and evolution.

Life emerged on the Earth within the first quintile of its habitable window, but a technological civilization did not blossom until its last. Efforts to infer the rate of abiogenesis, based on its early emergence, are frustrated by the selection effect that if the evolution of intelligence is a slow process, then life's early start may simply be a prerequisite to our existence, rather than useful evidence for optimism. In this work, we interpret the chronology of these two events in a Bayesian framework, extending upon previous work by considering that the evolutionary timescale is itself an unknown that needs to be jointly inferred, rather than fiducially set. We further adopt an objective Bayesian approach, such that our results would be agreed upon even by those using wildly different priors for the rates of abiogenesis and evolution - common points of contention for this problem. It is then shown that the earliest microfossil evidence for life indicates that the rate of abiogenesis is at least 2.8 times more likely to be a typically rapid process, rather than a slow one. This modest limiting Bayes factor rises to 8.7 if we accept the more disputed evidence of C13 depleted zircon deposits (Bell et al. 2015). For intelligence evolution, it is found that a rare-intelligence scenario is slightly favored at 3:2 betting odds. Thus, if we re-ran Earth's clock, one should statistically favor life to frequently re-emerge, but intelligence may not be as inevitable.

The Galactic bulge is a tumultuous dense region of space, packed with stars separated by far smaller distances than those in the Solar neighborhood. A quantification of the frequency and proximity of close stellar encounters in this environment dictates the exchange of material, disruption of planetary orbits, and threat of sterilizing energetic events. We present estimated encounter rates for stars in the Milky Way bulge found using a combination of numerical and analytical methods. By integrating the orbits of bulge stars with varying orbital energy and angular momentum to find their positions over time, we were able to estimate how many close stellar encounters the stars should experience as a function of orbit shape. We determined that ~80% of bulge stars have encounters within 1000 AU and that half of bulge stars will have >35 such encounters, both over a gigayear. Our work has interesting implications for the long-term survivability of planets in the Galactic bulge.

Exomoons are predicted to produce transit timing variations (TTVs) upon their host planet. Unfortunately, so are many other astrophysical phenomena - most notably other planets in the system. In this work, an argument of reductio ad absurdum is invoked, by deriving the transit timing effects that are impossible for a single exomoon to produce. Our work derives three key analytic tests. First, one may exploit the fact that a TTV signal from an exomoon should be accompanied by transit duration variations (TDVs), and that one can derive a TDV floor as a minimum expected level of variability. Cases for which the TDV upper limit is below this floor can thus be killed as exomoon candidates. Second, formulae are provided for estimating whether moons are expected to be 'killable' when no TDVs presently exist, thus enabling the community to estimate whether it's even worth deriving TDVs in the first place. Third, a TTV ceiling is derived, above which exomoons should never be able to produce TTV amplitudes. These tools are applied to a catalog of TTVs and TDVs for two and half thousand Kepler Objects Interest, revealing over two hundred cases that cannot be due to a moon. These tests are also applied to the exomoon candidate Kepler-1625b i, which comfortably passes the criteria. These simple analytic results should provide a means of rapidly rejecting putative exomoons and streamlining the search for satellites.

Photometric phase curves provide an important window onto exoplanetary atmospheres and potentially even their surfaces. With similar amplitudes to occultations but far longer baselines, they have a higher sensitivity to planetary photons at the expense of a more challenging data reduction in terms of long-term stability. In this work, we introduce a novel non-parametric algorithm dubbed "phasma" to produce clean, robust exoplanet phase curves and apply it to 115 Neptunian and 50 Terran exoplanets observed by Kepler. We stack the signals to further improve signal-to-noise, and measure an average Neptunian albedo of Ag < 0.23 to 95% confidence, indicating a lack of bright clouds consistent with theoretical models. Our Terran sample provides the first constraint on the ensemble albedo of exoplanets which are most likely solid, constraining Ag < 0.42 to 95% confidence. In agreement with our constraint on the greenhouse effect, our work implies that Kepler's solid planets are unlikely to resemble cloudy Venusian analogs, but rather dark Mercurian rocks.

We report the detection of the full orbital phase curve and occultation of the hot-Jupiter WASP-100b using TESS photometry. The phase curve is isolated by suppressing low frequency stellar and instrumental modes using both a non-parametric harmonic notch filter (phasma) and semi-sector long polynomials. This yields a phase curve signal of (73 +/- 9) ppm amplitude, preferred over a null-model by deltaBIC = 25, indicating very strong evidence for an observed effect. We recover the occultation event with a suite of five temporally localized tools, including Gaussian processes and cosine filtering. This allows us to infer an occultation depth of (100 +/- 14) ppm, with an additional +/- 16 ppm systematic error from the differences between methods. We regress a model including atmospheric reflection, emission, ellipsoidal variations and Doppler beaming to the combined phase curve and occultation data. This allows us to infer that WASP-100b has a geometric albedo of A_g = 0.16 (+0.04, -0.03) in the TESS bandpass, with a maximum dayside brightness temperature of (2710 +/- 100) K and a warm nightside temperature of (2380 [+170, -200]) K. Additionally, we find evidence that WASP-100b has a high thermal redistribution efficiency, manifesting as a substantial eastward hotspot offset of (71 [+2, -4]) degrees. These results present the first measurement of a thermal phase shift among the phase curves observed by TESS so far, and challenge the predicted efficiency of heat transport in the atmospheres of ultra-hot Jupiters.

The claim of an exomoon candidate in the Kepler-1625b system has generated substantial discussion regarding possible alternative explanations for the purported signal. In this work we examine in detail these possibilities. First, the effect of more flexible trend models is explored and we show that sufficiently flexible models are capable of attenuating the signal, although this is an expected byproduct of invoking such models. We also explore trend models using X and Y centroid positions and show that there is no data-driven impetus to adopt such models over temporal ones. We quantify the probability that the 500 ppm moon-like dip could be caused by a Neptune-sized transiting planet to be < 0.75%. We show that neither autocorrelation, Gaussian processes nor a Lomb-Scargle periodogram are able to recover a stellar rotation period, demonstrating that K1625 is a quiet star with periodic behavior < 200 ppm. Through injection and recovery tests, we find that the star does not exhibit a tendency to introduce false-positive dip-like features above that of pure Gaussian noise. Finally, we address a recent re-analysis by Kreidberg et al (2019) and show that the difference in conclusions is not from differing systematics models but rather the reduction itself. We show that their reduction exhibits i) slightly higher intra-orbit and post-fit residual scatter, ii) $\simeq$ 900 ppm larger flux offset at the visit change, iii) $\simeq$ 2 times larger Y-centroid variations, and iv) $\simeq$ 3.5 times stronger flux-centroid correlation coefficient than the original analysis. These points could be explained by larger systematics in their reduction, potentially impacting their conclusions.

Since 2009, the Kepler, K2, and TESS missions have produced a vast number of lightcurves for public use. To assist citizen scientists in processing those lightcurves, the LcTools software system was developed. The system provides a set of tools to efficiently search for signals of interest in large sets of lightcurves using automated and manual (visual) techniques. At the heart of the system is a multipurpose lightcurve viewer and signal processor with advanced navigation and display capabilities to facilitate the search for signals. Other applications in the system are available for building lightcurve files in bulk, finding periodic signals automatically, and generating signal reports. This paper describes each application in the system and the methods by which the software can be used to detect and record signals. The software is free and can be obtained from the lead author by request.

Distant starlight passing through the Earth's atmosphere is refracted by an angle of just over one degree near the surface. This focuses light onto a focal line starting at an inner (and chromatic) boundary out to infinity - offering an opportunity for pronounced lensing. It is shown here that the focal line commences at ~85% of the Earth-Moon separation, and thus placing an orbiting detector between here and one Hill radius could exploit this refractive lens. Analytic estimates are derived for a source directly behind the Earth (i.e. on-axis) showing that starlight is lensed into a thin circular ring of thickness $W H_{\Delta}/R$, yielding an amplification of $8 H_{\Delta}/W$, where $H_{\Delta}$ is the Earth's refractive scale height, $R$ is its geopotential radius and $W$ is the detector diameter. These estimates are verified through numerical ray-tracing experiments from optical to 30 micron light with standard atmospheric models. The numerical experiments are extended to include extinction from both a clear atmosphere and one with clouds. It is found that a detector at one Hill radius is least affected by extinction since lensed rays travel no deeper than 13.7 km, within the stratosphere and above most clouds. Including extinction, a 1 metre Hill radius 'terrascope' is calculated to produce an amplification of ~45,000 for a lensing timescale of ~20 hours. In practice, the amplification is likely halved in order to avoid daylight scattering i.e. 22,500 ($\Delta$mag=10.9) for $W=$1 metre, or equivalent to a 150 metre optical/infrared telescope.

We demonstrate that planet occurrence does not depend on stellar galactocentric velocity in the Solar neighborhood. Using Gaia DR2 astrometry and radial velocity data, we calculate 3D galactocentric velocities for 197,090 Kepler field stars, 1647 of which are confirmed planet hosts. When we compare the galactocentric velocities of planet hosts to those of the entire field star sample, we observe a statistically significant (KS p-value = $10^{-70}$) distinction, with planet hosts being apparently slower than field stars by $\sim$40 km/s. We explore some potential explanations for this difference and conclude that it is not a consequence of the planet-metallicity relation or distinctions in the samples' thin/thick disk membership, but rather an artefact of Kepler's selection function. Non Kepler-host stars that have nearly identical distances, temperatures, surface gravities, and Kepler magnitudes to the confirmed planet hosts also have nearly identical velocity distributions. Using one of these identical non-host samples, we consider that the probability of a star with velocity $v_{\mathrm{tot}}$ hosting a planet can be described by an exponential function proportional to $e^{(-v_{\mathrm{tot}}/v_0)}$. Using a Markov Chain Monte Carlo sampler, we determine that $v_0 > $976 km/s to 99% confidence, which implies that planets in the Solar neighborhood are just as likely to form around fast stars as they are around slow stars. Our work highlights the subtle ways in which selection biases can create strong correlations without physical underpinnings.

The true multiplicity distribution of transiting planet systems is obscured by strong observational biases, leading low-multiplicity systems to be overrepresented in the observed sample. Using the Kepler FGK planet hosts, we employ approximate Bayesian computation to infer the multiplicity distribution by comparing simulated catalogs to the observed one. After comparing a total of ten different multiplicity distributions, half of which were two-population models, to the observed data, we find that a single-population model following a Zipfian distribution is able to explain the Kepler data as well as any of the dichotomous models we test. Our work provides another example of a way to explain the observed Kepler multiplicities without invoking a dichotomous planet population. Using our preferred Zipfian model, we estimate that an additional $2393_{-717}^{+904}$ planets likely reside in the 1537 FGK Kepler systems studied in this work, which would increase the planet count by a factor of $2.22_{-0.36}^{+0.46}$. Of these hidden worlds, $663_{-151}^{+158}$ are expected to reside in ostensibly single-transiting-planet systems, meaning that an additional planet(s) is expected for approximately 1-in-2 such Kepler systems.

Kepler, K2, TESS, and similar time-domain photometric projects, while designed with exoplanet detection in mind, are also well-suited projects for searches for large artificial structures orbiting other stars in the Galaxy. An effort to examine these data sets with an eye towards non-spherical or otherwise anomalous transit events, and a robust follow-up program to understand the stars and occulters that generate them, would enable the first robust upper limits on such megastructures in terms of their sizes, occurrence rates, and orbital properties. Such work also has the ancillary benefit of improving our understanding of stellar photometric variability and orbital and physical parameter estimation of exoplanets from photometric time series, and may lead to the identification of new, unexpected classes of stellar variables and exoplanets. Ultimately, searching for the most unusual and anomalous signatures benefits not only the search for technologies, but also the entire astronomical community by uncovering new mysteries to advance our understanding of the Universe.

We report Spitzer Space Telescope observations during predicted transits of the exoplanet Proxima Centauri b. As the nearest terrestrial habitable-zone planet we will ever discover, any potential transit of Proxima b would place strong constraints on its radius, bulk density, and atmosphere. Subsequent transmission spectroscopy and secondary-eclipse measurements could then probe the atmospheric chemistry, physical processes, and orbit, including a search for biosignatures. However, our photometric results rule out planetary transits at the 200~ppm level at 4.5$~{\mu}m$, yielding a 3$\sigma$ upper radius limit of 0.4~$R_\rm{\oplus}$ (Earth radii). Previous claims of possible transits from optical ground- and space-based photometry were likely correlated noise in the data from Proxima Centauri's frequent flaring. Follow-up observations should focus on planetary radio emission, phase curves, and direct imaging. Our study indicates dramatically reduced stellar activity at near-to-mid infrared wavelengths, compared to the optical. Proxima b is an ideal target for space-based infrared telescopes, if their instruments can be configured to handle Proxima's brightness.

We present the discovery of a pair of transiting giant planets, TOI-216b and c, using four sectors of TESS photometry. TOI-216 is a $0.87 M_{\odot}$ dwarf orbited by two transiters with radii of $8.2 R_{\oplus}$ and $11.3 R_{\oplus}$, and periods of $17.01$d and $34.57$d, respectively. Anti-correlated TTVs are clearly evident indicating that the transiters orbit the same star and interact via a near 2:1 mean motion resonance. By fitting the TTVs with a dynamical model, we infer masses of $26_{-11}^{+24} M_{\oplus}$ and $190_{-80}^{+220} M_{\oplus}$, establishing that the objects are planetary in nature and have likely sub-Kronian and Kronian densities. TOI-216 lies close to the southern ecliptic pole and thus will be observed by TESS throughout the first year, providing an opportunity for continuous dynamical monitoring and considerable refinement of the dynamical masses presented here. TOI-216 closely resembles Kepler-9 in architecture, and we hypothesize that in such systems these Saturn-analogs failed to fully open a gap and thus migrated far deeper into the system before becoming trapped into resonance, which would imply that future detections of new analogs may also have sub-Jupiter masses.

By design, model-based approaches for flagging transiting exoplanets in light curves, such as boxed least squares, excel at detecting planets with low S/N at the expense of finding signals that are not well described by the assumed model, such as self-lensing binaries, disintegrating or evaporating planets, or planets with large rings. So far, such signals have typically been found through visual searches by professional or citizen scientists, or by inspection of the photometric power-spectra. We present a nonparametric detection algorithm, for short duty-cycle periodic signals in photometric time series based on phase dispersion minimization. We apply our code to 161,786 Kepler sources and detect 18 new periodic signals consistent with heartbeat binaries/planets, 4 new singly-transiting systems, and 2 new doubly-transiting systems. We show that our code is able to recover the majority of known Kepler objects of interest (KOIs) to high confidence, as well as more unusual events such as Boyajian's star and a comet passing through the Kepler field. Nonparametric signal-flagging techniques, such as the one presented here, will become increasingly valuable with the coming data from TESS and future transit surveys as the volume of data available to us exceeds that which can be feasibly examined manually.

Gravitational slingshots around a neutron star in a compact binary have been proposed as a means of accelerating large masses to potentially relativistic speeds. Such a slingshot is attractive since fuel is not expended for the acceleration, however it does entail a spacecraft diving into close proximity of the binary, which could be hazardous. It is proposed here that such a slingshot can be performed remotely using a beam of light which follows a boomerang null geodesic. Using a moving black hole as a gravitational mirror, kinetic energy from the black hole is transferred to the beam of light as a blueshift and upon return the recycled photons not only accelerate, but also add energy to, the spacecraft. It is shown here that this gained energy can be later expended to reach a terminal velocity of approximately 133% the velocity of the black hole. A civilization could exploit black holes as galactic way points but would be difficult to detect remotely, except for an elevated binary merger rate and excess binary eccentricity.

We consider the problem of inferring the shape of a transiting object's silhouette from its light curve alone, without assuming a physical model for the object. We model the object as a grid of pixels which transits the star; each pixel has an opacity, ranging from transparent to opaque, which we infer from the light curve. We explore three interesting degeneracies inherent to this problem, in which markedly different transiting shapes can produce identical light curves: (i) the "flip" degeneracy, by which two pixels transiting at the same impact parameter on opposite sides of the star's horizontal midplane generate the same light curve; (ii) the "arc" degeneracy, by which opacity can be redistributed along the semicircular arc of pixels which undergoes ingress or egress at the same time without consequence to the light curve, and (iii) the "stretch" degeneracy, by which a wide shape moving fast can produce the same light curve as a narrow shape moving more slowly. By understanding these degeneracies and adopting some additional assumptions, we are able to numerically recover informative shadow images of transiting objects, and we explore a number of different algorithmic approaches to this problem. We apply our methods to real data, including the TRAPPIST-1c,e,f triple transit and two dips of Boyajian's Star. We provide Python code to calculate the transit light curve of any grid and, conversely, infer the image grid which generates any light curve in a software package accompanying this paper, EightBitTransit.

LHS 1140 is a nearby mid-M dwarf known to host a temperate rocky super-Earth (LHS 1140 b) on a 24.737-day orbit. Based on photometric observations by MEarth and Spitzer as well as Doppler spectroscopy from HARPS, we report the discovery of an additional transiting rocky companion (LHS 1140 c) with a mass of $1.81\pm0.39~{\rm M_{Earth}}$ and a radius of $1.282\pm0.024~{\rm R_{Earth}}$ on a tighter, 3.77795-day orbit. We also obtain more precise estimates of the mass and radius of LHS 1140 b to be $6.98\pm0.89~{\rm M_{Earth}}$ and $1.727\pm0.032~{\rm R_{Earth}}$. The mean densities of planets b and c are $7.5\pm1.0~\rm{g/cm^3}$ and $4.7\pm1.1~\rm{g/cm^3}$, respectively, both consistent with the Earth's ratio of iron to magnesium silicate. The orbital eccentricities of LHS 1140 b and c are consistent with circular orbits and constrained to be below 0.06 and 0.31, respectively, with 90% confidence. Because the orbits of the two planets are co-planar and because we know from previous analyses of Kepler data that compact systems of small planets orbiting M dwarfs are commonplace, a search for more transiting planets in the LHS 1140 system could be fruitful. LHS 1140 c is one of the few known nearby terrestrial planets whose atmosphere could be studied with the upcoming James Webb Space Telescope.

Exomoons are the natural satellites of planets orbiting stars outside our solar system, of which there are currently no confirmed examples. We present new observations of a candidate exomoon associated with Kepler-1625b using the Hubble Space Telescope to validate or refute the moon's presence. We find evidence in favor of the moon hypothesis, based on timing deviations and a flux decrement from the star consistent with a large transiting exomoon. Self-consistent photodynamical modeling suggests that the planet is likely several Jupiter masses, while the exomoon has a mass and radius similar to Neptune. Since our inference is dominated by a single but highly precise Hubble epoch, we advocate for future monitoring of the system to check model predictions and confirm repetition of the moon-like signal.

White dwarfs are one of the few types of stellar objects for which we know almost nothing about the possible existence of companion planets. Recent evidence for metal contaminated atmospheres, circumstellar debris disks and transiting planetary debris all indicate that planets may be likely. However, white dwarf transit surveys are challenging due to the intrinsic faintness of such objects, the short timescale of the transits and the low transit probabilities due to their compact radii. The Large Synoptic Survey Telescope (LSST) offers a remedy to these problems as a deep, half-sky survey with fast exposures encompassing approximately 10 million white dwarfs with $r<24.5$ apparent magnitude. We simulate LSST photometric observations of 3.5 million white dwarfs over a ten-year period and calculate the detectability of companion planets with $P<10$ d via transits. We find typical detection rates in the range of $5 \times 10^{-6}$ to $4 \times 10^{-4}$ for Ceres-sized bodies to Earth-sized worlds, yielding $\sim 50$ to $4\,000$ detections for a 100% occurrence rate of each. For terrestrial planets in the continuously habitable zone, we find detection rates of $\sim 10^{-3}$ indicating that LSST would reveal hundreds of such worlds for occurrence rates in the range of 1% to 10%.

The TRAPPIST-1 system provides an exquisite laboratory for advancing our understanding exoplanetary atmospheres, compositions, dynamics and architectures. A remarkable aspect of TRAPPIST-1 is that it represents the longest known resonance chain, where all seven planets share near mean motion resonances with their neighbors. Prior to the measurement of 1h's period, Luger et al. (2017) showed that six possible and highly precise periods for 1h were expected, assuming it also participated in the resonant chain. We show here that combining this argument with a Titius-Bode law fit of the inner six worlds narrows the choices down to a single precise postdiction for 1h's period, which is ultimately the correct period. But a successful postdiction is never as convincing as a successful prediction, and so we take the next step and apply this argument to a hypothetical TRAPPIST-1i. We find two possible periods predicted by this argument, either 25.345 or 28.699 days. If successful, this may provide the basis for planet prediction in compact resonant chain systems. If falsified, this would indicate that the argument lacks true predictive power and may not be worthwhile pursuing further in our efforts to build predictive models for planetary systems.

Life appears to have emerged relatively quickly on the Earth, a fact sometimes used to justify a high rate of spontaneous abiogenesis ($\lambda$) among Earth-like worlds. Conditioned upon a single datum - the time of earliest evidence for life ($t_{\mathrm{obs}}$) - previous Bayesian formalisms for the posterior distribution of $\lambda$ have demonstrated how inferences are highly sensitive to the priors. Rather than attempt to infer the true $\lambda$ posterior, we here compute the relative change to $\lambda$ when new experimental/observational evidence is introduced. By simulating posterior distributions and resulting entropic information gains, we compare three experimental pressures on $\lambda$: 1) evidence for an earlier start to life; $t_{\mathrm{obs}}$; 2) constraints on spontaneous abiogenesis from the lab; and 3) an exoplanet survey for biosignatures. First, we find that experiments 1 and 2 can only yield lower limits on $\lambda$, unlike 3. Second, evidence for an earlier start to life can yield negligible information on $\lambda$ if $t_{\mathrm{obs}} \ll \lambda_{\mathrm{max}}^{-1}$. Vice versa, experiment 2 is uninformative when $\lambda_{\mathrm{max}} \gg t_{\mathrm{obs}}^{-1}$. Whilst experiment 3 appears the most direct means of measuring $\lambda$, we highlight that early starts inform us of the conditions of abiogenesis, and that lab experiments could succeed in building new life. Together then, the three experiments are complimentary and we encourage activity in all to solve this grand challenge.

An exoplanet's optical phase curve constrains the thermal emission and albedo of the planet's surface and/or atmosphere, as well as potentially constraining the mass via gravitational influences on the host star. Recently, Jansen & Kipping (2018) demonstrated that exoplanets with precisely constrained orbital periods, as is typical of transiting planets, enable one to exploit a Fourier-argument to non-parametrically separate optical phase curves from long-term photometric time series in the presence of realistic noise structures. This algorithm, dubbed phasma, was applied to 477 confirmed exoplanets in the original paper. In this research note, we extend the sample to all non-false positive Kepler Objects of Interest with periods less than ten days. Our sample of 1998 detrended phase curves are made publicly available to help the community search for various effects. A Python port of the phasma code is also made available.

Studying photometric time series in the frequency domain can serve as a means of detecting rotational modulations, measuring asteroseismic modes and even detecting short-period transiting planets. To our knowledge, there is no prior archive of the NASA Kepler Mission's power spectra and so we present one here to aid the community in searching for such effects. Using DR25 PDC long-cadence Kepler photometry, 2,594,616 individual periodograms are computed using Welch's method with a Nuttall window, where we provide a unique periododogram for each quarter (up to 16) of each star (196,791 in total). Additionally, we normalize the periodograms in the high-frequency end and combine them into channel- and quarter-averaged power spectra to track common instrumental modes occurring onboard the telescope, with a particularly notable feature at ~80 minutes (~200 $\mu$Hz) observed.

The TRAPPIST-1 system provides an exquisite laboratory for understanding exoplanetary atmospheres and interiors. Their mutual gravitational interactions leads to transit timing variations, from which Grimm et al. (2018) recently measured the planetary masses with precisions ranging from 5% to 12%. Using these masses and the <5% radius measurements on each planet, we apply the method described in Suissa et al. (2018) to infer the minimum and maximum CRF (core radius fraction) of each planet. Further, we modify the maximum limit to account for the fact that a light volatile envelope is excluded for planets b through f. Only planet e is found to have a significant probability of having a non-zero minimum CRF, with a 0.7% false-alarm probability it has no core. Our method further allows us to measure the CRF of planet e to be greater than (49 +/- 7)% but less than (72 +/- 2)%, which is compatible with that of the Earth. TRAPPIST-1e therefore possess a large iron core similar to the Earth, in addition to being Earth-sized and located in the temperature zone.

The interior structure of an exoplanet is hidden from direct view yet likely plays a crucial role in influencing the habitability of the Earth analogues. Inferences of the interior structure are impeded by a fundamental degeneracy that exists between any model comprising more than two layers and observations constraining just two bulk parameters: mass and radius. In this work, we show that although the inverse problem is indeed degenerate, there exists two boundary conditions that enables one to infer the minimum and maximum core radius fraction, CRFmin and CRFmax. These hold true even for planets with light volatile envelopes, but require the planet to be fully differentiated and that layers denser than iron are forbidden. With both bounds in hand, a marginal CRF can also be inferred by sampling in-between. After validating on the Earth, we apply our method to Kepler-36b and measure CRFmin = (0.50 +/- 0.07), CRFmax = (0.78 +/- 0.02), and CRFmarg = (0.64 +/- 0.11), broadly consistent with the Earth's true CRF value of 0.55. We apply our method to a suite of hypothetical measurements of synthetic planets to serve as a sensitivity analysis. We find that CRFmin and CRFmax have recovered uncertainties proportional to the relative error on the planetary density, but CRFmarg saturates to between 0.03 and 0.16 once delta rho/rho drops below 1-2 per cent. This implies that mass and radius alone cannot provide any better constraints on internal composition once bulk density constraints hit around a per cent, providing a clear target for observers.

The search for extrasolar planets in the past decades has shown that planets abound in the Solar neighborhood. While we are still missing an Earth twin, the forthcoming space missions and ground-based instrumentation are already driven to achieve this goal. But, in order to fully understand the conditions for life appearing in the Solar System, we still miss some pieces of the planetary system jigsaw puzzle, namely a deeper understanding of the minor bodies. Trojans, moons, and comets are tracers of the formation and evolution processes of planetary systems. These missing pieces are also critical to understand the emergence and evolution of life over millions of years. With the large crop of planetary systems discovered so far and yet to be detected with the forthcoming missions, the hunt for minor bodies in extrasolar systems is a natural continuation of our search for real Solar System- and, in particular, Earth- analogs. This white paper is focused on detection of these minor components and their relevance in the emergence, evolution and survival of life.

Jupiter played an important role in determining the structure and configuration of the Solar System. Whereas hot-Jupiter type exoplanets preferentially form around metal-rich stars, the conditions required for the formation of planets with masses, orbits and eccentricities comparable to Jupiter (Jupiter analogues) are unknown. Using spectroscopic metallicities, we show that stars hosting Jupiter analogues have an average metallicity close to solar, in contrast to their hot-Jupiter and eccentric cool Jupiter counterparts, which orbit stars with super-solar metallicities. Furthermore, the eccentricities of Jupiter analogues increase with host star metallicity, suggesting that planet-planet scatterings producing highly eccentric cool Jupiters could be more common in metal-rich environments. To investigate a possible explanation for these metallicity trends, we compare the observations to numerical simulations, which indicate that metal-rich stars typically form multiple Jupiters, leading to planet-planet interactions and, hence, a prevalence of either eccentric cool Jupiters or hot-Jupiters with circularized orbits. Although the samples are small and exhibit variations in their metallicities, suggesting that numerous processes other than metallicity affect the formation of planetary systems, the data in hand suggests that Jupiter analogues and terrestrial-sized planets form around stars with average metallicities close to solar, whereas high metallicity systems preferentially host eccentric cool Jupiter or hot-Jupiters, indicating higher metallicity systems may not be favorable for the formation of planetary systems akin to the Solar System.

Long-period circumbinary planets appear to be as common as those orbiting single stars and have been found to frequently have orbital radii just beyond the critical distance for dynamical stability. Assessing the stability is typically done either through N-body simulations or using the stability criterion first considered by Dvorak and later developed by Holman and Wiegert: a second-order polynomial calibrated to broadly match numerical simulations. However, the polynomial is unable to capture islands of instability introduced by mean motion resonances, causing the accuracy of the criterion to approach that of a random coin-toss when close to the boundary. We show how a deep neural network (DNN) trained on N-body simulations generated with REBOUND is able to significantly improve stability predictions for circumbinary planets on initially coplanar, circular orbits. Specifically, we find that the accuracy of our DNN never drops below 86%, even when tightly surrounding the boundary of instability, and is fast enough to be practical for on-the-fly calls during likelihood evaluations typical of modern Bayesian inference. Our binary classifier DNN is made publicly available at https://github.com/CoolWorlds/orbital-stability

Mountain ranges, volcanoes, trenches, and craters are common on rocky bodies throughout the Solar System, and we might we expect the same for rocky exoplanets. With ever larger telescopes under design and a growing need to not just detect planets but also to characterize them, it is timely to consider whether there is any prospect of remotely detecting exoplanet topography in the coming decades. To test this, we devised a novel yet simple approach to detect and quantify topographical features on the surfaces of exoplanets using transit light curves. If a planet rotates as it transits its parent star, its changing silhouette yields a time-varying transit depth, which can be observed as an apparent and anomalous increase in the photometric scatter. Using elevation data for several rocky bodies in our solar system, we quantify each world's surface integrated relief with a "bumpiness" factor, and calculate the corresponding photometric scatter expected during a transit. Here we describe the kinds of observations that would be necessary to detect topography in the ideal case of Mars transiting a nearby white dwarf star. If such systems have a conservative occurrence rate of 10%, we estimate that the upcoming Colossus or OWL telescopes would be able to detect topography with <20 hours of observing time, which corresponds to ~400 transits with a duration of 2 minutes and orbital period of ~10 hours.

The Characterizing Exoplanets Satellite (CHEOPS) mission is planned for launch next year with a major objective being to search for transits of known RV planets, particularly those orbiting bright stars. Since the radial velocity method is only sensitive to planetary mass, the radii, transit depths and transit signal-to-noise values of each RV planet are, a-priori, unknown. Using an empirically calibrated probabilistic mass-radius relation, forecaster (Chen & Kipping 2017a), we address this by predicting a catalog of homogeneous credible intervals for these three keys terms for 468 planets discovered via radial velocities. Of these, we find that the vast majority should be detectable with CHEOPS, including terrestrial bodies, if they have the correct geometric alignment. In particular, we predict that 22 mini-Neptunes and 82 Neptune-sized planets would be suitable for detection and that more than 80% of these will have apparent magnitude of V < 10, making them highly suitable for follow-up characterization work. Our work aims to assist the CHEOPS team in scheduling efforts and highlights the great value of quantifiable, statistically robust estimates for upcoming exoplanetary missions.

One proposed method for spacecraft to reach nearby stars is by accelerating sails using either solar radiation pressure or directed energy. This idea constitutes the thesis behind the Breakthrough Starshot project, which aims to accelerate a gram-mass spacecraft up to one-fifth the speed of light towards Proxima Centauri. For such a case, the combination of the sail's low mass and relativistic velocity render previous treatments formally incorrect, including that of Einstein himself in his seminal 1905 paper introducing special relativity. To address this, we present formulae for a sail's acceleration, first in response to a single photon and then extended to an ensemble. We show how the sail's motion in response to an ensemble of incident photons is equivalent to that of a single photon of energy equal to that of the ensemble. We use this 'principle of ensemble equivalence' for both perfect and imperfect mirrors, enabling a simple analytic prediction of the sail's velocity curve. Using our results and adopting putative parameters for Starshot, we estimate that previous relativistic treatments underestimate the spacecraft's terminal velocity by ~50m/s for the same incident energy, sufficient to miss a target by several Earth radii. Additionally, we use a simple model to predict the sail's temperature and diffraction beam losses during the laser firing period, allowing us to estimate that for firing times of a few minutes and operating temperatures below 300C (573K), Starshot will require a sail of which absorbs less than 1 in 260,000 photons.

The properties of a transiting planet's host star are written in its transit light curve. The light curve can reveal the stellar density and the limb darkening profile in addition to the characteristics of the planet and its orbit. For planets with strong prior constraints on orbital eccentricity, we may measure these stellar properties directly from the light curve; this method promises to aid greatly in the characterization of transiting planet host stars targeted by the upcoming NASA TESS mission and any long-period, singly-transiting planets discovered in the same systems. Using Bayesian inference, we fit a transit model, including a nonlinear limb darkening law, to 66 Kepler transiting planet hosts to measure their stellar properties. We present posterior distributions of stellar density, limb-darkening coefficients, and other system parameters for these stars. We measure densities to within 5% for the majority of our target stars, with the dominant precision-limiting factor being the signal-to-noise ratio of the transits. 95 % of our measured stellar densities are in 3-sigma or better agreement with previously published literature values. We make posterior distributions for all of our target KOIs available for online download.

Exomoons represent an outstanding challenge in modern astronomy, with the potential to provide rich insights into planet formation theory and habitability. In this work, we stack the phase-folded transits of 284 viable moon hosting Kepler planetary candidates, in order to search for satellites. These planets range from Earth-to-Jupiter sized and from ${\sim}$0.1-to-1.0 AU in separation - so-called "warm" planets. Our data processing includes two-pass harmonic detrending, transit timing variations, model selection and careful data quality vetting to produce a grand light curve with an r.m.s. of 5.1 ppm. We find that the occurrence rate of Galilean-analog moon systems for planets orbiting between ${\sim}$0.1 and 1.0 AU can be constrained to be $\eta<0.38$ to 95% confidence for the 284 KOIs considered, with a 68.3% confidence interval of $\eta=0.16_{-0.10}^{+0.13}$. A single-moon model of variable size and separation locates a slight preference for a population of short-period moons with radii ${\sim}0.5$ $R_{\oplus}$ orbiting at 5-10 planetary radii. However, we stress that the low Bayes factor of just 2 in this region means it should be treated as no more than a hint at this time. Splitting our data into various physically-motivated subsets reveals no strong signal. The dearth of Galilean-analogs around warm planets places the first strong constraint on exomoon formation models to date. Finally, we report evidence for an exomoon candidate Kepler-1625b I, which we briefly describe ahead of scheduled observations of the target with the Hubble Space Telescope.

Recent transit surveys have discovered thousands of planetary candidates with directly measured radii, but only a small fraction have measured masses. Planetary mass is crucial in assessing the feasibility of numerous observational signatures, such as radial velocities (RVs), atmospheres, moons and rings. In the absence of a direct measurement, a data-driven, probabilistic forecast enables observational planning and so here we compute posterior distributions for the forecasted mass of approximately seven thousand Kepler Objects of Interest (KOIs). Our forecasts reveal that the predicted RV amplitudes of Neptunian planets are relatively consistent, as a result of transit survey detection bias, hovering around the few m/s level. We find that mass forecasts are unlikely to improve through more precise planetary radii, with the error budget presently dominated by the intrinsic model uncertainty. Our forecasts identify a couple of dozen KOIs near the Terran-Neptunian divide with particularly large RV semi-amplitudes which could be promising targets to follow-up, particularly in the near-IR. With several more transit surveys planned in the near-future, the need to quickly forecast observational signatures is likely to grow and the work here provides a template example of such calculations.

One of the most directly observable features of a transiting multi-planet system is their size-ordering when ranked in orbital separation. Kepler has revealed a rich diversity of outcomes, from perfectly ordered systems, like Kepler-80, to ostensibly disordered systems, like Kepler-20. Under the hypothesis that systems are born via preferred formation pathways, one might reasonably expect non-random size-orderings reflecting these processes. However, subsequent dynamical evolution, often chaotic and turbulent in nature, may erode this information and so here we ask - do systems remember how they formed? To address this, we devise a model to define the entropy of a planetary system's size-ordering, by first comparing differences between neighboring planets and then extending to accommodate differences across the chain. We derive closed-form solutions for many of the micro state occupancies and provide public code with look-up tables to compute entropy for up to ten-planet systems. All three proposed entropy definitions exhibit the expected property that their credible interval increases with respect to a proxy for time. We find that the observed Kepler multis display a highly significant deficit in entropy compared to a randomly generated population. Incorporating a filter for systems deemed likely to be dynamically packed, we show that this result is robust against the possibility of missing planets too. Put together, our work establishes that Kepler systems do indeed remember something of their younger years and highlights the value of information theory for exoplanetary science.

We use 16 quarters of the \textit{Kepler} mission data to analyze the transit timing variations (TTVs) of the extrasolar planet Kepler-46b (KOI-872). Our dynamical fits confirm that the TTVs of this planet (period $P=33.648^{+0.004}_{-0.005}$ days) are produced by a non-transiting planet Kepler-46c ($P=57.325^{+0.116}_{-0.098}$ days). The Bayesian inference tool \texttt{MultiNest} is used to infer the dynamical parameters of Kepler-46b and Kepler-46c. We find that the two planets have nearly coplanar and circular orbits, with eccentricities $\simeq 0.03$ somewhat higher than previously estimated. The masses of the two planets are found to be $M_{b}=0.885^{+0.374}_{-0.343}$ and $M_{c}=0.362^{+0.016}_{-0.016}$ Jupiter masses, with $M_{b}$ being determined here from TTVs for the first time. Due to the precession of its orbital plane, Kepler-46c should start transiting its host star in a few decades from now.

The analysis of Proxima Centauri's radial velocities recently led Anglada-Escud\'e et al. (2016) to claim the presence of a low mass planet orbiting the Sun's nearest star once every 11.2 days. Although the a-priori probability that Proxima b transits its parent star is just 1.5%, the potential impact of such a discovery would be considerable. Independent of recent radial velocity efforts, we observed Proxima Centauri for 12.5 days in 2014 and 31 days in 2015 with the MOST space telescope. We report here that we cannot make a compelling case that Proxima b transits in our precise photometric time series. Imposing an informative prior on the period and phase, we do detect a candidate signal with the expected depth. However, perturbing the phase prior across 100 evenly spaced intervals reveals one strong false-positive and one weaker instance. We estimate a false-positive rate of at least a few percent and a much higher false-negative rate of 20-40%, likely caused by the very high flare rate of Proxima Centauri. Comparing our candidate signal to HATSouth ground-based photometry reveals that the signal is somewhat, but not conclusively, disfavored (1-2 sigmas) leading us to argue that the signal is most likely spurious. We expect that infrared photometric follow-up could more conclusively test the existence of this candidate signal, owing to the suppression of flare activity and the impressive infrared brightness of the parent star.

The upcoming TESS mission is expected to find thousands of transiting planets around bright stars, yet for three-quarters of the fields observed the temporal coverage will limit discoveries to planets with orbital periods below 13.7 days. From the Kepler catalog, the mean probability of these short-period transiting planets having additional longer period transiters (which would be missed by TESS) is 18%, a value ten times higher than the average star. In this work, we show how this probability is not uniform but functionally dependent upon the properties of the observed short-period transiters, ranging from less than 1% up to over 50%. Using artificial neural networks (ANNs) trained on the Kepler catalog and making careful feature selection to account for the differing sensitivity of TESS, we are able to predict the most likely short-period transiters to be accompanied by additional transiters. Through cross-validation, we predict that a targeted, optimized TESS transit and/or radial velocity follow-up program using our trained ANN would have a discovery yield improved by a factor of two. Our work enables a near-optimal follow-up strategy for surveys following TESS targets for additional planets, improving the science yield derived from TESS and particularly beneficial in the search for habitable-zone transiting worlds.

Mass and radius are two of the most fundamental properties of an astronomical object. Increasingly, new planet discoveries are being announced with a measurement of one of these terms, but not both. This has led to a growing need to forecast the missing quantity using the other, especially when predicting the detectability of certain follow-up observations. We present am unbiased forecasting model built upon a probabilistic mass-radius relation conditioned on a sample of 316 well-constrained objects. Our publicly available code, Forecaster, accounts for observational errors, hyper-parameter uncertainties and the intrinsic dispersions observed in the calibration sample. By conditioning our model upon a sample spanning dwarf planets to late-type stars, Forecaster can predict the mass (or radius) from the radius (or mass) for objects covering nine orders-of-magnitude in mass. Classification is naturally performed by our model, which uses four classes we label as Terran worlds, Neptunian worlds, Jovian worlds and stars. Our classification identifies dwarf planets as merely low-mass Terrans (like the Earth), and brown dwarfs as merely high-mass Jovians (like Jupiter). We detect a transition in the mass-radius relation at $2.0_{-0.6}^{+0.7} M_\oplus$, which we associate with the divide between solid, Terran worlds and Neptunian worlds. This independent analysis adds further weight to the emerging consensus that rocky Super-Earths represent a narrower region of parameter space than originally thought. Effectively, then, the Earth is the Super-Earth we have been looking for.

We present a study of white light flares from the active M5.5 dwarf Proxima Centauri using the Canadian microsatellite MOST. Using 37.6 days of monitoring data from 2014 and 2015, we have detected 66 individual flare events, the largest number of white light flares observed to date on Proxima Cen. Flare energies in our sample range from $10^{29}$-$10^{31.5}$ erg. The flare rate is lower than that of other classic flare stars of similar spectral type, such as UV Ceti, which may indicate Proxima Cen had a higher flare rate in its youth. Proxima Cen does have an unusually high flare rate given its slow rotation period, however. Extending the observed power-law occurrence distribution down to $10^{28}$ erg, we show that flares with flux amplitudes of 0.5% occur 63 times per day, while superflares with energies of $10^{33}$ erg occur ~8 times per year. Small flares may therefore pose a great difficulty in searches for transits from the recently announced 1.27 M_earth Proxima b, while frequent large flares could have significant impact on the planetary atmosphere.

Observational biases distort our view of nature, such that the patterns we see within a surveyed population of interest are often unrepresentative of the truth we seek. Transiting planets currently represent the most informative data set on the ensemble properties of exoplanets within 1 AU of their star. However, the transit method is inherently biased due to both geometric and detection-driven effects. In this work, we derive the overall observational biases affecting the most basic transit parameters from first principles. By assuming a trapezoidal transit and using conditional probability, we infer the expected distribution of these terms both as a joint distribution and in a marginalized form. These general analytic results provide a baseline against which to compare trends predicted by mission-tailored injection/recovery simulations and offer a simple way to correct for observational bias. Our results explain why the observed population of transiting planets displays a non-uniform impact parameter distribution, with a bias towards near-equatorial geometries. We also find that the geometric bias towards observed planets transiting near periastron is attenuated by the longer durations which occur near apoastron. Finally, we predict that the observational bias with respect to ratio-of-radii is super-quadratic, scaling as $(R_P/R_{\star})^{5/2}$, driven by an enhanced geometric transit probability and modestly longer durations.

The relations between observable stellar parameters are usually assumed to be deterministic. That is, given an infinitely precise measurement of independent variable, `$x$', and some model, the value of dependent variable, `$y$' can be known exactly. In practise this assumption is rarely valid and intrinsic stochasticity means that two stars with exactly the same `$x$', will have slightly different `$y$'s. The relation between short-timescale brightness fluctuations (flicker) of stars and both surface gravity and stellar density are two such stochastic relations that have, until now, been treated as deterministic ones. We recalibrate these relations in a probabilistic framework, using Hierarchical Bayesian Modelling (HBM) to constrain the intrinsic scatter in the relations. We find evidence for additional scatter in the relationships, that cannot be accounted for by the observational uncertainties alone. The scatter in surface gravity and stellar density does not depend on flicker, suggesting that using flicker as a proxy for $\log g$ and $\rho_\star$ is equally valid for dwarf and giant stars, despite the fact that the observational uncertainties tend to be larger for dwarfs.

The transit method is presently the most successful planet discovery and characterization tool at our disposal. Other advanced civilizations would surely be aware of this technique and appreciate that their home planet's existence and habitability is essentially broadcast to all stars lying along their ecliptic plane. We suggest that advanced civilizations could cloak their presence, or deliberately broadcast it, through controlled laser emission. Such emission could distort the apparent shape of their transit light curves with relatively little energy, due to the collimated beam and relatively infrequent nature of transits. We estimate that humanity could cloak the Earth from Kepler-like broadband surveys using an optical monochromatic laser array emitting a peak power of about 30 MW for roughly 10 hours per year. A chromatic cloak, effective at all wavelengths, is more challenging requiring a large array of tunable lasers with a total power of approximately 250 MW. Alternatively, a civilization could cloak only the atmospheric signatures associated with biological activity on their world, such as oxygen, which is achievable with a peak laser power of just around 160 kW per transit. Finally, we suggest that the time of transit for optical SETI is analogous to the water-hole in radio SETI, providing a clear window in which observers may expect to communicate. Accordingly, we propose that a civilization may deliberately broadcast their technological capabilities by distorting their transit to an artificial shape, which serves as both a SETI beacon and a medium for data transmission. Such signatures could be readily searched in the archival data of transit surveys.

Decadal-long radial velocity surveys have recently started to discover analogs to the most influential planet of our solar system, Jupiter. Detecting and characterizing these worlds is expected to shape our understanding of our uniqueness in the cosmos. Despite the great successes of recent transit surveys, Jupiter analogs represent a terra incognita, owing to the strong intrinsic bias of this method against long orbital periods. We here report on the first validated transiting Jupiter analog, Kepler-167e (KOI-490.02), discovered using Kepler archival photometry orbiting the K4-dwarf KIC-3239945. With a radius of $(0.91\pm0.02)$ $R_{\mathrm{Jup}}$, a low orbital eccentricity ($0.06_{-0.04}^{+0.10}$) and an equilibrium temperature of $(131\pm3)$ K, Kepler-167e bears many of the basic hallmarks of Jupiter. Kepler-167e is accompanied by three Super-Earths on compact orbits, which we also validate, leaving a large cavity of transiting worlds around the habitable-zone. With two transits and continuous photometric coverage, we are able to uniquely and precisely measure the orbital period of this post snow-line planet ($1071.2323\pm0.0006$ d), paving the way for follow-up of this $K=11.8$ mag target.

Determining which small exoplanets have stony-iron compositions is necessary for quantifying the occurrence of such planets and for understanding the physics of planet formation. Kepler-10 hosts the stony-iron world Kepler-10b (K10b), and also contains what has been reported to be the largest solid silicate-ice planet, Kepler-10c (K10c). Using 220 radial velocities (RVs), including 72 precise RVs from Keck-HIRES of which 20 are new from 2014-2015, and 17 quarters of Kepler photometry, we obtain the most complete picture of the Kepler-10 system to date. We find that K10b (Rp=1.47 Re) has mass 3.72$\pm$0.42 Me and density 6.46$\pm$0.73 g/cc. Modeling the interior of K10b as an iron core overlaid with a silicate mantle, we find that the iron core constitutes 0.17$\pm$0.11 of the planet mass. For K10c (Rp=2.35 Re) we measure Mp=13.98$\pm$1.79 Me and $\rho$=5.94$\pm$0.76 g/cc, significantly lower than the mass computed in Dumusque et al. (2014, 17.2$\pm$1.9 Me). Internal compositional modeling reveals that at least $10\%$ of the radius of Kepler-10c is a volatile envelope composed of hydrogen-helium ($0.2\%$ of the mass, $16\%$ of the radius) or super-ionic water ($28\%$ of the mass, $29\%$ of the radius). Analysis of only HIRES data yields a higher mass for K10b and a lower mass for K10c than does analysis of the HARPS-N data alone, with the mass estimates for K10c formally inconsistent by 3$\sigma$. Splitting the RVs from each instrument leads to inconsistent measurements for the mass of planet c in each data set. This suggests that time-correlated noise is present and that the uncertainties in the planet masses (especially K10c) exceed our formal estimates. Transit timing variations (TTVs) of K10c indicate the likely presence of a third planet in the system, KOI-72.X. The TTVs and RVs are consistent with KOI-72.X having an orbital period of 24, 71, or 101 days, and a mass from 1-7 Me.

Stellar limb-darkening impacts a wide range of astronomical measurements. The accuracy to which it is modelled limits the accuracy in any covariant parameters of interest, such as the radius of a transiting planet. With the ever growing availability of precise observations and the importance of robust estimates of astrophysical parameters, an emerging trend has been to freely fit the limb-darkening coefficients (LDCs) describing a limb-darkening law of choice, in order to propagate our ignorance of the true intensity profile. In practice, this approach has been limited to two-parameter limb-darkening laws, such as the quadratic law, due to the relative ease of sampling the physically allowed range of LDCs. Here, we provide a highly efficient method for sampling LDCs describing a more accurate three-parameter non-linear law. We first derive analytic criteria which can quickly test if a set of LDCs are physical, although naive sampling with these criteria leads to an acceptance rate less than 1%. We then show that the loci of allowed LDCs can be transformed into a cone-like volume, from which we are able to draw uniform samples. We show that samples drawn uniformly from the conal region are physically valid in 97.3% of realizations and encompass 94.4% of the volume of allowed parameter space. We provide Python and Fortran code (LDC3) to sample from this region (and perform the reverse calculation) at https://github.com/davidkipping/LDC3, which also includes a subroutine to efficiently test whether a sample is physically valid or not.

White dwarfs are the end state of most stars, including the Sun, after they exhaust their nuclear fuel. Between 1/4 and 1/2 of white dwarfs have elements heavier than helium in their atmospheres, even though these elements should rapidly settle into the stellar interiors unless they are occasionally replenished. The abundance ratios of heavy elements in white dwarf atmospheres are similar to rocky bodies in the Solar system. This and the existence of warm dusty debris disks around about 4% of white dwarfs suggest that rocky debris from white dwarf progenitors' planetary systems occasionally pollute the stars' atmospheres. The total accreted mass can be comparable to that of large asteroids in the solar system. However, the process of disrupting planetary material has not yet been observed. Here, we report observations of a white dwarf being transited by at least one and likely multiple disintegrating planetesimals with periods ranging from 4.5 hours to 4.9 hours. The strongest transit signals occur every 4.5 hours and exhibit varying depths up to 40% and asymmetric profiles, indicative of a small object with a cometary tail of dusty effluent material. The star hosts a dusty debris disk and the star's spectrum shows prominent lines from heavy elements like magnesium, aluminium, silicon, calcium, iron, and nickel. This system provides evidence that heavy element pollution of white dwarfs can originate from disrupted rocky bodies such as asteroids and minor planets.

We present a survey of 41 Kepler Objects of Interest (KOIs) for exomoons using Bayesian photodynamics, more than tripling the number of KOIs surveyed with this technique. We find no compelling evidence for exomoons although thirteen KOIs yield spurious detections driven by instrumental artifacts, stellar activity and/or perturbations from unseen bodies. Regarding the latter, we find seven KOIs exhibiting >5 sigma evidence of transit timing variations, including the 'mega-Earth' Kepler-10c, likely indicating an additional planet in that system. We exploit the moderately large sample of 57 unique KOIs surveyed to date to infer several useful statistics. For example, although there is a diverse range in sensitivities, we find that we are sensitive to Pluto-Charon mass-ratio systems for ~40% of KOIs studied and Earth-Moon mass-ratios for 1 in 8 cases. In terms of absolute mass, our limits probe down to 1.7 Ganymede masses, with a sensitivity to Earth-mass moons for 1 in 3 cases studied and to the smallest moons capable of sustaining an Earth-like atmosphere (0.3 Earth masses) for 1 in 4. Despite the lack of positive detections to date, we caution against drawing conclusions yet, since our most interesting objects remain under analysis. Finally, we point out that had we searched for the photometric transit signals of exomoons alone, rather than using photodynamics, we estimate that 1 in 4 KOIs would have erroneously been concluded to harbor exomoons due to residual time correlated noise in the Kepler data, posing a serious problem for alternative methods.

We report the discovery of Kepler-432b, a giant planet ($M_b = 5.41^{+0.32}_{-0.18} M_{\rm Jup}, R_b = 1.145^{+0.036}_{-0.039} R_{\rm Jup}$) transiting an evolved star $(M_\star = 1.32^{+0.10}_{-0.07} M_\odot, R_\star = 4.06^{+0.12}_{-0.08} R_\odot)$ with an orbital period of $P_b = 52.501129^{+0.000067}_{-0.000053}$ days. Radial velocities (RVs) reveal that Kepler-432b orbits its parent star with an eccentricity of $e = 0.5134^{+0.0098}_{-0.0089}$, which we also measure independently with asterodensity profiling (AP; $e=0.507^{+0.039}_{-0.114}$), thereby confirming the validity of AP on this particular evolved star. The well-determined planetary properties and unusually large mass also make this planet an important benchmark for theoretical models of super-Jupiter formation. Long-term RV monitoring detected the presence of a non-transiting outer planet (Kepler-432c; $M_c \sin{i_c} = 2.43^{+0.22}_{-0.24} M_{\rm Jup}, P_c = 406.2^{+3.9}_{-2.5}$ days), and adaptive optics imaging revealed a nearby (0\farcs87), faint companion (Kepler-432B) that is a physically bound M dwarf. The host star exhibits high signal-to-noise asteroseismic oscillations, which enable precise measurements of the stellar mass, radius and age. Analysis of the rotational splitting of the oscillation modes additionally reveals the stellar spin axis to be nearly edge-on, which suggests that the stellar spin is likely well-aligned with the orbit of the transiting planet. Despite its long period, the obliquity of the 52.5-day orbit may have been shaped by star-planet interaction in a manner similar to hot Jupiter systems, and we present observational and theoretical evidence to support this scenario. Finally, as a short-period outlier among giant planets orbiting giant stars, study of Kepler-432b may help explain the distribution of massive planets orbiting giant stars interior to 1 AU.

The discovery of rings around extrasolar planets ("exorings") is one of the next breakthroughs in exoplanetary research. Previous studies have explored the feasibility of detecting exorings with present and future photometric sensitivities by seeking anomalous deviations in the residuals of a standard transit light curve fit, at the level of ~100 ppm for Kronian rings. In this work, we explore two much larger observational consequences of exorings: (1) the significant increase in transit depth that may lead to the misclassification of ringed planetary candidates as false-positives and/or the underestimation of planetary density; and (2) the so-called "photo-ring" effect, a new asterodensity profiling effect, revealed by a comparison of the light curve derived stellar density to that measured with independent methods (e.g., asteroseismology). While these methods do not provide an unambiguous detection of exorings, we show that the large amplitude of these effects combined with their relatively simple analytic description, makes them highly suited to large-scale surveys to identify candidate ringed planets worthy of more detailed investigation. Moreover, these methods lend themselves to ensemble analyses seeking to uncover evidence of a population of ringed planets. We describe the method in detail, develop the basic underlying formalism and test it in the parameter space of rings and transit config- uration. We discuss the prospects of using this method for the first systematic search of exoplanetary rings in the Kepler database and provide a basic computational code for implementing it.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

We present an investigation of twelve candidate transiting planets from Kepler with orbital periods ranging from 34 to 207 days, selected from initial indications that they are small and potentially in the habitable zone (HZ) of their parent stars. Few of these objects are known. The expected Doppler signals are too small to confirm them by demonstrating that their masses are in the planetary regime. Here we verify their planetary nature by validating them statistically using the BLENDER technique, which simulates large numbers of false positives and compares the resulting light curves with the Kepler photometry. This analysis was supplemented with new follow-up observations (high-resolution optical and near-infrared spectroscopy, adaptive optics imaging, and speckle interferometry), as well as an analysis of the flux centroids. For eleven of them (KOI-0571.05, 1422.04, 1422.05, 2529.02, 3255.01, 3284.01, 4005.01, 4087.01, 4622.01, 4742.01, and 4745.01) we show that the likelihood they are true planets is far greater than that of a false positive, to a confidence level of 99.73% (3 sigma) or higher. For KOI-4427.01 the confidence level is about 99.2% (2.6 sigma). With our accurate characterization of the GKM host stars, the derived planetary radii range from 1.1 to 2.7 R_Earth. All twelve objects are confirmed to be in the HZ, and nine are small enough to be rocky. Excluding three of them that have been previously validated by others, our study doubles the number of known rocky planets in the HZ. KOI-3284.01 (Kepler-438b) and KOI-4742.01 (Kepler-442b) are the planets most similar to the Earth discovered to date when considering their size and incident flux jointly.

The discovery of an exomoon would provide deep insights into planet formation and the habitability of planetary systems, with transiting examples being particularly sought after. Of the hundreds of Kepler planets now discovered, the seven-planet system Kepler-90 is unusual for exhibiting an unidentified transit-like signal in close proximity to one of the transits of the long-period gas-giant Kepler-90g, as noted by Cabrera et. al. (2014). As part of the 'Hunt for Exomoons with Kepler' (HEK) project, we investigate this possible exomoon signal and find it passes all conventional photometric, dynamical and centroid diagnostic tests. However, pixel-level light curves indicate that the moon-like signal occurs on nearly all of the target's pixels, which we confirm using a novel way of examining pixel-level data which we dub the 'transit centroid'. This test reveals that the possible exomoon to Kepler-90g is likely a false positive, perhaps due to a cosmic ray induced Sudden Pixel Sensitivity Dropout (SPSD). This work highlights the extreme care required for seeking non-periodic low-amplitude transit signals, such as exomoons.

The diversity and quantity of moons in the Solar System suggest a manifold population of natural satellites exist around extrasolar planets. Of peculiar interest from an astrobiological perspective, the number of sizable moons in the stellar habitable zones may outnumber planets in these circumstellar regions. With technological and theoretical methods now allowing for the detection of sub-Earth-sized extrasolar planets, the first detection of an extrasolar moon appears feasible. In this review, we summarize formation channels of massive exomoons that are potentially detectable with current or near-future instruments. We discuss the orbital effects that govern exomoon evolution, we present a framework to characterize an exomoon's stellar plus planetary illumination as well as its tidal heating, and we address the techniques that have been proposed to search for exomoons. Most notably, we show that natural satellites in the range of 0.1 - 0.5 Earth mass (i) are potentially habitable, (ii) can form within the circumplanetary debris and gas disk or via capture from a binary, and (iii) are detectable with current technology.

Planets on eccentric orbits have a higher geometric probability of transiting their host star. By application of Bayes' theorem, we reverse this logic to show that the eccentricity distribution of transiting planets is positively biased. Adopting the flexible Beta distribution as the underlying prior for eccentricity, we derive the marginalized transit probability as well as the a-priori joint probability distribution of eccentricity and argument of periastron, given that a planet is known to transit. These results allow to demonstrate that most planet occurrence rate calculations using Kepler data have overestimated the prevalence of planets by ~10%. Indeed, the true occurrence of planets from transit surveys is fundamentally intractable without a prior assumption for the eccentricity distribution. Further more, we show that previously extracted eccentricity distributions using Kepler data are positively biased. In cases where one wishes to impose an informative eccentricity prior, we provide a recursive algorithm to apply inverse transform sampling of our joint prior probability distribution. Computer code of this algorithm, ECCSAMPLES, is provided to enable the community to sample directly from the prior.

In most theories of planet formation, the snow-line represents a boundary between the emergence of the interior rocky planets and the exterior ice giants. The wide separation of the snow-line makes the discovery of transiting worlds challenging, yet transits would allow for detailed subsequent characterization. We present the discovery of Kepler-421b, a Uranus-sized exoplanet transiting a G9/K0 dwarf once every 704.2 days in a near-circular orbit. Using public Kepler photometry, we demonstrate that the two observed transits can be uniquely attributed to the 704.2 day period. Detailed light curve analysis with BLENDER validates the planetary nature of Kepler-421b to >4 sigmas confidence. Kepler-421b receives the same insolation as a body at ~2AU in the Solar System and for a Uranian albedo would have an effective temperature of ~180K. Using a time-dependent model for the protoplanetary disk, we estimate that Kepler-421b's present semi-major axis was beyond the snow-line after ~3Myr, indicating that Kepler-421b may have formed at its observed location.

The Hubble Space Telescope (HST) and the Kepler space mission observed a large number of planetary transits showing anomalies due to starspot eclipses, with more such observations expected in the near future by the K2 mission and the Transiting Exoplanet Survey Satellite (TESS). To facilitate analysis of this phenomenon, we present spotrod, a model for planetary transits of stars with an arbitrary limb darkening law and a number of homogeneous, circular spots on their surface. A free, open source implementation written in C, ready to use in Python, is available for download. We analyze Kepler observations of the planetary host star HAT-P-11, and study the size and contrast of more than two hundred starspots. We find that the flux ratio of spots ranges at least from 0.6 to 0.9, corresponding to an effective temperature approximately 100 to 450 K lower than the stellar surface, although it is possible that some spots are darker than 0.5. The largest detected spots have a radius less than approximately 0.2 stellar radii.

Eclipsing systems, such as transiting exoplanets, allow one to measure the mean stellar density of the host star under various idealized assumptions. Asterodensity Profiling (AP) compares this density to an independently determined value in order to check the validity of the assumptions and ultimately derive useful parameters. Several physical effects can cause said assumptions to become invalid, with the most well-known example being the so-called photo-eccentric effect. In this work, we provide analytic expressions for five other effects which induce AP deviations: the photo-blend, -spot, -timing, -duration and -mass effects. We find that these effects can easily reproduce large AP deviations and so we caution that extracting the eccentricity distribution is only viable with careful consideration of the prior distributions for these other effects. We also re-investigate the photo-eccentric effect and derive a single-domain minimum eccentricity expression and the parameter range for which analytic formulae are valid. This latter result shows that the assumptions underlying the analytic model for the photo-eccentric effect break down for close-in, highly-eccentric planets, meaning that extreme care must be taken in this regime. Finally, we demonstrate that contaminated light fraction can be solved for, indicating that AP could be a potent tool for planet validation.

KOI-227, KOI-319 and KOI-884 are identified here as (at least) two planet systems. For KOI-319 and KOI-884, the observed Transit Timing Variations (TTVs) of the inner transiting planet are used to detect an outer non-transiting planet. The outer planet in KOI-884 is 2.6 Jupiter masses and has the orbital period just narrow of the 3:1 resonance with the inner planet (orbital period ratio 2.93). The distribution of parameters inferred from KOI-319.01's TTVs is bimodal with either a 1.6 Neptune-mass (M_N) planet wide of the 5:3 resonance (period 80.1 d) or a 1 Saturn-mass planet wide of the 7:3 resonance (period 109.2 d). The radial velocity measurements can be used in this case to determine which of these parameter modes is correct. KOI-227.01's TTVs with large 10 hour amplitude can be obtained for planetary-mass companions in various major resonances. Based on the Bayesian evidence, the current TTV data favor the outer 2:1 resonance with a companion mass 1.5 M_N, but this solution implies a very large density of KOI-227.01. The inner and outer 3:2 resonance solutions with sub-Neptune-mass companions are physically more plausible, but will need to be verified.

Asterodensity Profiling (AP) is a relatively new technique for studying transit light curves. By comparing the mean stellar density derived from the transit light curve to that found through some independent method, AP provides information on several useful properties such as orbital eccentricity and blended light. We present an AP survey of 41 Kepler Objects of Interest (KOIs), with a single transiting candidate, for which the target star's mean stellar density has been measured using asteroseismology. The ensemble distribution of the AP measurements for the 31 dwarf stars in our sample shows excellent agreement with the spread expected if the KOIs were genuine and have realistic eccentricities. In contrast, the same test for the 10 giants in our sample reveals significant incompatibility at 4$\sigma$ confidence. Whilst extreme eccentricities could be invoked, this hypothesis requires four of the KOIs to contact their host star at periastron passage, including the recently claimed confirmation of Kepler-91b. After carefully examining several hypotheses, we conclude that the most plausible explanation is that the transiting objects orbit a different star to that measured with asteroseismology - cases we define as false-positives. Based on the AP distribution, we estimate a false positive rate for Kepler's giant stars with a single transiting object of FPR$\simeq70\pm30$%.

Two decades ago, astronomers began detecting planets orbiting stars other than our Sun, so-called exoplanets. Since that time, the rate of detections and the sensitivity to ever-smaller planets has improved dramatically with several Earth-sized planets now known. As our sensitivity dives into the terrestrial regime, increasingly the community has wondered if the moons of exoplanets may also be detectable, so-called "exomoons". Their detection represents an outstanding challenge in modern astronomy and would provide deep insights into the uniqueness of our Solar System and perhaps even expand the definition of habitability. Here, I will briefly review theoretical studies exploring the formation and evolution of exomoons, which serve to guide observational searches and provide testable hypotheses. Next, I will outline the different methods which have been proposed to accomplish this challenging feat and their respective merits. Finally, initial results from observational efforts will be summarized with a view to future prospects as well.

We explore how finite integration times or equivalently temporal binning induces morphological distortions to the transit light-curve. These distortions, if uncorrected for, lead to the retrieval of erroneous system parameters and may even lead to some planetary candidates being rejected as ostensibly unphysical. We provide analytic expressions for estimating the disturbance to the various light-curve parameters as a function of the integration time. These effects are particularly crucial in light of the long-cadence photometry often used for discovering new exoplanets by, for example, Convection Rotation and Planetary Transits (COROT) and the Kepler Mission (8.5 and 30 min). One of the dominant effects of long integration times is a systematic underestimation of the light-curve-derived stellar density, which has significant ramifications for transit surveys. We present a discussion of numerical integration techniques to compensate for the effects and produce expressions to quickly estimate the errors of such techniques, as a function of integration time and numerical resolution. This allows for an economic choice of resolution before attempting fits of long-cadence light-curves. We provide a comparison of the short- and long-cadence light-curves of TrES-2b and show that the retrieved transit parameters are consistent using the techniques discussed here.

A number of planet-host stars have been observed to rotate with a period equal to an integer multiple of the orbital period of their close planet. We expand this list by analyzing Kepler data of HAT-P-11 and finding a period ratio of 6:1. In particular, we present evidence for a long-lived spot on the stellar surface that is eclipsed by the planet in the same position four times, every sixth transit. We also identify minima in the out-of-transit lightcurve and confirm that their phase with respect to the stellar rotation is mostly stationary for the 48-month timeframe of the observations, confirming the proposed rotation period. For comparison, we apply our methods to Kepler-17 and confirm the findings of Bonomo & Lanza (2012) that the period ratio is not exactly 8:1 in that system. Finally, we provide a hypothesis on how interactions between a star and its planet could possibly result in an observed commensurability for systems where the stellar differential rotation profile happens to include a period at some latitude which is commensurable to the planetary orbit.

Variability in the time series brightness of a star on a timescale of 8 hours, known as 'flicker', has been previously demonstrated to serve as a proxy for the surface gravity of a star by Bastien et al. (2013). Although surface gravity is crucial for stellar classification, it is the mean stellar density which is most useful when studying transiting exoplanets, due to its direct impact on the transit light curve shape. Indeed, an accurate and independent measure of the stellar density can be leveraged to infer subtle properties of a transiting system, such as the companion's orbital eccentricity via asterodensity profiling. We here calibrate flicker to the mean stellar density of 439 Kepler targets with asteroseismology, allowing us to derive a new empirical relation given by $\log_{10}(\rho_{\star}\,[\mathrm{kg}\,\mathrm{m}^{-3}]) = 5.413 - 1.850 \log_{10}(F_8\,[\mathrm{ppm}])$. The calibration is valid for stars with $4500$K$<T_{\mathrm{eff}}<6500$K, $K_P<14$ and flicker estimates corresponding to stars with $3.25<\log g_{\star}<4.43$. Our relation has a model error in the stellar density of 31.7% and so has $\sim8$ times lower precision than that from asteroseismology but is applicable to a sample $\sim40$ times greater. Flicker therefore provides an empirical method to enable asterodensity profiling on hundreds of planetary candidates from present and future missions.

With their smaller radii and high cosmic abundance, transiting planets around cool stars hold a unique appeal. As part of our on-going project to measure the occurrence rate of extrasolar moons, we here present results from a survey focussing on eight Kepler planetary candidates associated with M-dwarfs. Using photodynamical modeling and Bayesian multimodal nested sampling, we find no compelling evidence for an exomoon in these eight systems. Upper limits on the presence of such bodies probe down to $\sim0.4M_{\oplus}$ in the best case. For KOI-314, we are able to confirm the planetary nature of two out of the three known transiting candidates using transit timing variations. Of particular interest is KOI-314c, which is found to have a mass of $1.0_{-0.3}^{+0.4}M_{\oplus}$, making it the lowest mass transiting planet discovered to date. With a radius of $1.61_{-0.15}^{+0.16}R_{\oplus}$, this Earth-mass world is likely enveloped by a significant gaseous envelope comprising $\geq17_{-13}^{+12}$% of the planet by radius. We find evidence to support the planetary nature of KOI-784 too via transit timing, but we advocate further observations to verify the signals. In both systems, we infer that the inner planet has a higher density than the outer world, which may be indicative of photo-evaporation. These results highlight both the ability of Kepler to search for sub-Earth mass moons and the exciting ancillary science which often results from such efforts.

Kepler-22b is the first transiting planet to have been detected in the habitable-zone of its host star. At 2.4 Earth radii, Kepler-22b is too large to be considered an Earth-analog, but should the planet host a moon large enough to maintain an atmosphere, then the Kepler-22 system may yet possess a telluric world. Aside from being within the habitable-zone, the target is attractive due to the availability of previously measured precise radial velocities and low intrinsic photometric noise, which has also enabled asteroseismology studies of the star. For these reasons, Kepler-22b was selected as a target-of-opportunity by the 'Hunt for Exomoons with Kepler' (HEK) project. In this work, we conduct a photodynamical search for an exomoon around Kepler-22b leveraging the transits, radial velocities and asteroseismology plus several new tools developed by the HEK project to improve exomoon searches. We find no evidence for an exomoon around the planet and exclude moons of mass >0.5 Earth masses to 95% confidence. By signal injection and blind retrieval, we demonstrate that an Earth-like moon is easily detected for this planet even when the time-correlated noise of the data set is taken into account. We provide updated parameters for the planet Kepler-22b including a revised mass of <53 Earth masses to 95% confidence and an eccentricity of 0.13(-0.13)(+0.36) by exploiting Single-body Asterodensity Profiling (SAP). Finally, we show that Kepler-22b has a >95% probability of being within the empirical habitable-zone but a <5% probability of being within the conservative habitable-zone.

Stellar limb darkening affects a wide range of astronomical measurements and is frequently modelled with a parametric model using polynomials in the cosine of the angle between the line of sight and the emergent intensity. Two-parameter laws are particularly popular for cases where one wishes to fit freely for the limb darkening coefficients (i.e. an uninformative prior) due to the compact prior volume and the fact that more complex models rarely obtain unique solutions with present data. In such cases, we show that the two limb darkening coefficients are constrained by three physical boundary conditions, describing a triangular region in the two-dimensional parameter space. We show that uniformly distributed samples may be drawn from this region with optimal efficiency by a technique developed by computer graphical programming: triangular sampling. Alternatively, one can make draws using a uniform, bivariate Dirichlet distribution. We provide simple expressions for these parametrizations for both techniques applied to the case of quadratic, square-root and logarithmic limb darkening laws. For example, in the case of the popular quadratic law, we advocate fitting for q_1 = (u_1+u_2)^2 and q_2 = 0.5u_1(u_1+u_2)^{-1} with uniform priors in the interval [0,1] to implement triangular sampling easily. Employing these parametrizations allows one to derive model parameters which fully account for our ignorance about the intensity profile, yet never explore unphysical solutions, yielding robust and realistic uncertainty estimates. Furthermore, in the case of triangular sampling with the quadratic law, our parametrization leads to significantly reduced mutual correlations and provides an alternative geometric explanation as to why naively fitting the quadratic limb darkening coefficients precipitates strong correlations in the first place.

It is suggested that the distribution of orbital eccentricities for extrasolar planets is well-described by the Beta distribution. Several properties of the Beta distribution make it a powerful tool for this purpose. For example, the Beta distribution can reproduce a diverse range of probability density functions (PDFs) using just two shape parameters (a and b). We argue that this makes it ideal for serving as a parametric model in Bayesian comparative population analysis. The Beta distribution is also uniquely defined over the interval zero to unity, meaning that it can serve as a proper prior for eccentricity when analysing the observations of bound extrasolar planets. Using nested sampling, we find that the distribution of eccentricities for 396 exoplanets detected through radial velocity with high signal-to-noise is well-described by a Beta distribution with parameters a = 0.867+/-0.044 and b = 3.03+/-0.17. The Beta distribution is shown to be 3.7 times more likely to represent the underlying distribution of exoplanet eccentricities than the next best model: a Rayleigh + exponential distribution. The same data are also used in an example population comparison utilizing the Beta distribution, where we find that the short- and long-period planets are described by distinct Beta distributions at a confidence of 11.6 sigma and display a signature consistent with the effects of tidal circularization.

Amongst the many hundreds of transiting planet candidates discovered by the Kepler Mission, one finds a large number of candidates with sizes between that of the Earth and Neptune. The composition of these worlds is not immediately obvious with no Solar System analog to draw upon and there exists some ambiguity as to whether a given candidate is a rocky Super-Earth or a gas-enveloped Mini-Neptune. The potential scientific value and observability of the atmospheres of these two classes of worlds varies significantly and given the sheer number of candidates in this size-range, there is evidently a need for a quick, simple metric to rank whether the planets have an extended atmosphere or not. In this work, we propose a way to calculate the 'minimum atmospheric height' (R_MAH) using only a planet's radius and mass as inputs. We assume and exploit the boundary condition that the bulk composition of a solid/liquid Super-Earth cannot be composed of a material lighter than that of water. Mass-radius loci above a pure-water composition planet correspond to R_MAH>0. The statistical confidence of a planet maintaining an extended atmosphere can be therefore easily calculated to provide a simple ranking of target planets for follow-up observations. We also discuss how this metric can be useful in the interpretation of the spectra of observed planetary atmospheres.

With the detection of extrasolar moons (exomoons) on the horizon, it is important to consider their potential for habitability. If we consider the circumstellar Habitable Zone (HZ, often described in terms of planet semi-major axis and orbital eccentricity), we can ask, "How does the HZ for an Earth-like exomoon differ from the HZ for an Earth-like exoplanet?" For the first time, we use 1D latitudinal energy balance modelling to address this question. The model places an Earth-like exomoon in orbit around a Jupiter mass planet, which in turn orbits a Sun-like star. The exomoon's surface temperature is evolved under the action of stellar insolation, atmospheric cooling, heat diffusion, eclipses and tidal heating. We use this model to carry out two separate investigations. In the first, four test cases are run to investigate in detail the dependence of the exomoon climate on orbital direction, orbital inclination, and on the frequency of stellar eclipse by the host planet. We find that lunar orbits which are retrograde to the planetary orbit exhibit greater climate variations than prograde orbits, with global mean temperatures around 0.1 K warmer due to the geometry of eclipses. If eclipses become frequent relative to the atmospheric thermal inertia timescale, climate oscillations become extremely small. In the second investigation, we carry out an extensive parameter study, running the model many times to study the habitability of the exomoon in the 4-dimensional space composed of the planet semi-major axis and eccentricity, and the moon semi-major axis and eccentricity. We find that for zero moon eccentricity, frequent eclipses allow the moon to remain habitable in regions of high planet eccentricity, but tidal heating severely constrains habitability in the limit of high moon eccentricity, making the habitable zone a sensitive function of moon semi-major axis.

The Transit Timing Variations (TTVs) can be used as a diagnostic of gravitational interactions between planets in a multi-planet system. Many Kepler Objects of Interest (KOIs) exhibit significant TTVs, but KOI-142.01 stands out among them with an unrivaled, 12-hour TTV amplitude. Here we report a thorough analysis of KOI-142.01's transits. We discover periodic Transit Duration Variations (TDVs) of KOI-142.01 that are nearly in phase with the observed TTVs. We show that KOI-142.01's TTVs and TDVs uniquely detect a non-transiting companion with a mass 0.7 that of Jupiter (KOI-142c). KOI-142.01's mass inferred from the transit variations is consistent with the measured transit depth, suggesting a Neptune class planet (KOI-142b). The orbital period ratio P_c/P_b=2.03 indicates that the two planets are just wide of the 2:1 resonance. The present dynamics of this system, characterized here in detail, can be used to test various formation theories that have been proposed to explain the near-resonant pairs of exoplanets.

From the list of 2321 transiting planet candidates announced by the Kepler Mission, we select seven targets with favorable properties for the capacity to dynamically maintain an exomoon and present a detectable signal. These seven candidates were identified through our automatic target selection (TSA) algorithm and target selection prioritization (TSP) filtering, whereby we excluded systems exhibiting significant time-correlated noise and focussed on those with a single transiting planet candidate of radius less than 6 Earth radii. We find no compelling evidence for an exomoon around any of the seven KOIs but constrain the satellite-to-planet mass ratios for each. For four of the seven KOIs, we estimate a 95% upper quantile of M_S/M_P<0.04, which given the radii of the candidates, likely probes down to sub-Earth masses. We also derive precise transit times and durations for each candidate and find no evidence for dynamical variations in any of the KOIs. With just a few systems analyzed thus far in the in-going HEK project, projections on eta-moon would be premature, but a high frequency of large moons around Super-Earths/Mini-Neptunes would appear to be incommensurable with our results so far.

Photometric rotational modulations due to starspots remain the most common and accessible way to study stellar activity. In the Kepler-era, there now exists precise, continuous photometry of ~150,000 stars presenting an unprecedented opportunity for statistical analyses of these modulations. Modelling rotational modulations allows one to invert the observations into several basic parameters, such as the rotation period, spot coverage, stellar inclination and differential rotation rate. The most widely used analytic model for this inversion comes from Budding (1977) and Dorren (1987), who considered circular, grey starspots for a linearly limb darkened star. In this work, we extend the model to be more suitable in the analysis of high precision photometry, such as that by Kepler. Our new freely available Fortran code, macula, provides several improvements, such as non-linear limb darkening of the star and spot, a single-domain analytic function, partial derivatives for all input parameters, temporal partial derivatives, diluted light compensation, instrumental offset normalisations, differential rotation, starspot evolution and predictions of transit depth variations due to unocculted spots. Through numerical testing, we find that the inclusion of non-linear limb darkening means macula has a maximum photometric error an order-of-magnitude less than that of Dorren (1987), for Sun-like stars observed in the Kepler-bandpass. The code executes three orders-of-magnitude faster than comparable numerical codes making it well-suited for inference problems.

The Kepler Mission is monitoring the brightness of ~150,000 stars searching for evidence of planetary transits. As part of the "Hunt for Exomoons with Kepler" (HEK) project, we report a planetary system with two confirmed planets and one candidate planet discovered using the publicly available data for KOI-872. Planet b transits the host star with a period P_b=33.6d and exhibits large transit timing variations indicative of a perturber. Dynamical modeling uniquely detects an outer nontransiting planet c near the 5:3 resonance (P_c=57.0d) of mass 0.37 times that of Jupiter. Transits of a third planetary candidate are also found: a 1.7-Earth radius super-Earth with a 6.8d period. Our analysis indicates a system with nearly coplanar and circular orbits, reminiscent of the orderly arrangement within the solar system.

We present multi-band optical photometry of 94 spectroscopically-confirmed Type Ia supernovae (SN Ia) in the redshift range 0.0055 to 0.073, obtained between 2006 and 2011. There are a total of 5522 light curve points. We show that our natural system SN photometry has a precision of roughly 0.03 mag or better in BVr'i', 0.06 mag in u', and 0.07 mag in U for points brighter than 17.5 mag and estimate that it has a systematic uncertainty of 0.014, 0.010, 0.012, 0.014, 0.046, and 0.073 mag in BVr'i'u'U, respectively. Comparisons of our standard system photometry with published SN Ia light curves and comparison stars reveal mean agreement across samples in the range of ~0.00-0.03 mag. We discuss the recent measurements of our telescope-plus-detector throughput by direct monochromatic illumination by Cramer et al (in prep.). This technique measures the whole optical path through the telescope, auxiliary optics, filters, and detector under the same conditions used to make SN measurements. Extremely well-characterized natural-system passbands (both in wavelength and over time) are crucial for the next generation of SN Ia photometry to reach the 0.01 mag accuracy level. The current sample of low-z SN Ia is now sufficiently large to remove most of the statistical sampling error from the dark energy error budget. But pursuing the dark-energy systematic errors by determining highly-accurate detector passbands, combining optical and near-infrared (NIR) photometry and spectra, using the nearby sample to illuminate the population properties of SN Ia, and measuring the local departures from the Hubble flow will benefit from larger, carefully measured nearby samples.

The search for extrasolar planets is strongly motivated by the goal of characterizing how frequent habitable worlds and life may be within the Galaxy. Whilst much effort has been spent on searching for Earth-like planets, large moons may also be common, temperate abodes for life as well. The methods to detect extrasolar moons, or "exomoons" are more subtle than their planetary counterparts and in this thesis I aim to provide a method to find such bodies in transiting systems, which offer the greatest potential for detection. Before one can search for the tiny perturbations to the planetary signal, an understanding of the planetary transit must be established. Therefore, in Chapters 3 to 5 I discuss the transit model and provide several new insights. Chapter 4 presents new analytic expressions for the times of transit minima and the transit duration, which will be critical in the later search for exomoons. Chapter 5 discusses two sources of distortion to the transit signal, namely blending (with a focus on the previously unconsidered self-blending scenario) and light curve smearing due to long integration times. I provide methods to compensate for both of these effects, thus permitting for the accurate modelling of the planetary transit light curve. In Chapter 6, I discuss methods to detect exomoons through their gravitational influence on the host planet, giving rise to transit timing and duration variations (TTV and TDV). The previously known TTV effect is updated with a new model and the associated critical problems are outlined. I then predict a new effect, TDV, which solves these problems, making exomoon detection viable. Chapter 7 presents a feasibility study for detecting habitable-zone exomoons with Kepler, where it is found that moons down to 0.2 Earth masses are detectable. Finally, conclusions and future work are discussed in Chapter 8.

It has been previously shown that moons of extrasolar planets may be detectable with the Kepler Mission, for moon masses above ~0.2 Earth masses Kipping et al. 2009c. Transit timing effects have been formerly identified as a potent tool to this end, exploiting the dynamics of the system. In this work, we explore the simulation of transit light curves of a planet plus a single moon including not only the transit timing effects but also the light curve signal of the moon itself. We introduce our new algorithm, LUNA, which produces transit light curves for both bodies, analytically accounting for shadow overlaps, stellar limb darkening and planet-moon dynamical motion. By building the dynamics into the core of LUNA, the routine automatically accounts for transit timing/duration variations and ingress/egress asymmetries for not only the planet, but also the moon. We then generate some artificial data for two feasibly detectable hypothetical systems of interest: a i) prograde and ii) retrograde Earth-like moon around a habitable-zone Neptune for a M-dwarf system. We fit the hypothetical systems using LUNA and demonstrate the feasibility of detecting these cases with Kepler photometry.

Two decades ago, empirical evidence concerning the existence and frequency of planets around stars, other than our own, was absent. Since this time, the detection of extrasolar planets from Jupiter-sized to most recently Earth-sized worlds has blossomed and we are finally able to shed light on the plurality of Earth-like, habitable planets in the cosmos. Extrasolar moons may also be frequent habitable worlds but their detection or even systematic pursuit remains lacking in the current literature. Here, we present a description of the first systematic search for extrasolar moons as part of a new observational project called "The Hunt for Exomoons with Kepler" (HEK). The HEK project distills the entire list of known transiting planet candidates found by Kepler (2326 at the time of writing) down to the most promising candidates for hosting a moon. Selected targets are fitted using a multimodal nested sampling algorithm coupled with a planet-with-moon light curve modelling routine. By comparing the Bayesian evidence of a planet-only model to that of a planet-with-moon, the detection process is handled in a Bayesian framework. In the case of null detections, upper limits derived from posteriors marginalised over the entire prior volume will be provided to inform the frequency of large moons around viable planetary hosts, eta-moon. After discussing our methodologies for target selection, modelling, fitting and vetting, we provide two example analyses.

We present a novel method to determine eccentricity constraints of extrasolar planets in systems with multiple transiting planets through photometry alone. Our method is highly model independent, making no assumptions about the stellar parameters and requiring no radial velocity, transit timing or occultation events. Our technique exploits the fact the light curve derived stellar density must be the same for all planets transiting a common star. Assuming a circular orbit, the derived stellar density departs from the true value by a predictable factor, Psi, which contains information on the eccentricity of the system. By comparing multiple stellar densities, any differences must be due to eccentricity and thus meaningful constraints can be placed in the absence of any other information. The technique, dubbed "Multibody Asterodensity Profiling" (MAP), is a new observable which can be used alone or in combination with other observables, such as radial velocities and transit timing variations. An eccentricity prior can also be included as desired. MAP is most sensitive to the minimum pair-combined eccentricity e.g. (e1 + e2)[min]. Individual eccentricity constraints are less stringent but an empirical eccentricity posterior is always derivable and is freely available from transit photometry alone. We present a description of our technique using both analytic and numerical implementations, followed by two example analyses on synthetic photometry as a proof of principle. We point out that MAP has the potential to constrain the eccentricity, and thus habitability, of Earth-like planets in the absence of radial velocity data, which is likely for terrestrial-mass objects.

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

We present the detection of visible light from the planet TrES-2b, the darkest exoplanet currently known. By analysis of the orbital photometry from publicly available Kepler data (0.4-0.9 microns), we determine a day-night contrast amplitude of (6.5 +/- 1.9 ppm), constituting the lowest amplitude orbital phase variation discovered. The signal is detected to 3.7 sigma confidence and persists in six different methods of modelling the data and thus appears robust. In contrast, we are unable to detect ellipsoidal variations or beaming effects, but we do provide confidence intervals for these terms. If the day-night contrast is interpreted as being due to scattering, it corresponds to a geometric albedo of Ag = 0.0253 +/- 0.0072. However, our models indicate that there is a significant emission component to day-side brightness, and the true albedo is even lower (<1%). By combining our measurement with Spitzer and ground-based data, we show that a model with moderate redistribution (Pn ~ 0.3) and moderate extra optical opacity (kappa' ~ 0.3-0.4) provide a compatible explanation to the data.

We report the discovery of HAT-P-31b, a transiting exoplanet orbiting the V=11.660 dwarf star GSC 2099-00908. HAT-P-31b is the first HAT planet discovered without any follow-up photometry, demonstrating the feasibility of a new mode of operation for the HATNet project. The 2.17 Mj, 1.1Rj planet has a period P = 5.0054 days and maintains an unusually high eccentricity of e = 0.2450+/-0.0045, determined through Keck, FIES and Subaru high precision radial velocities. Detailed modeling of the radial velocities indicates an additional quadratic residual trend in the data detected to very high confidence. We interpret this trend as a long-period outer companion, HAT-P-31c, of minimum mass 3.4Mj and period >2.8 years. Since current RVs span less than half an orbital period, we are unable to determine the properties of HAT-P-31c to high confidence. However, dynamical simulations of two possible configurations show that orbital stability is to be expected. Further, if HAT-P-31c has non-zero eccentricity, our simulations show that the eccentricity of HAT-P-31b is actively driven by the presence of c, making HAT-P-31 a potentially intriguing dynamical laboratory.

We present an analysis of 18 short-cadence (SC) transit lightcurves of TrES-2b using quarter 0 (Q0) and quarter 1 (Q1) from the Kepler Mission. The photometry is of unprecedented precision, 237ppm per minute, allowing for the most accurate determination of the transit parameters yet obtained for this system. Global fits of the transit photometry, radial velocities and known transit times are used to obtain a self-consistent set of refined parameters for this system, including updated stellar and planetary parameters. Special attention is paid to fitting for limb darkening and eccentricity. We place an upper limit on the occultation depth to be <72.9ppm to 3-sigma confidence, indicating TrES-2b has the lowest determined geometric albedo for an exoplanet, of Ag<0.146. We also produce a transit timing analysis using Kepler's short-cadence data and demonstrate exceptional timing precision at the level of a few seconds for each transit event. With 18 fully-sampled transits at such high precision, we are able to produce stringent constraints on the presence of perturbing planets, Trojans and extrasolar moons. We introduce the novel use of control data to identify phasing effects. We also exclude the previously proposed hypotheses of short-period TTV and additional transits but find the hypothesis of long-term inclination change is neither supported nor refuted by our analysis.

We present an analysis of seven primary transit observations of the hot Neptune GJ436b at 3.6, 4.5 and $8~\mu$m obtained with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. After correcting for systematic effects, we fitted the light curves using the Markov Chain Monte Carlo technique. Combining these new data with the EPOXI, HST and ground-based $V, I, H$ and $K_s$ published observations, the range $0.5-10~\mu$m can be covered. Due to the low level of activity of GJ436, the effect of starspots on the combination of transits at different epochs is negligible at the accuracy of the dataset. Representative climate models were calculated by using a three-dimensional, pseudo-spectral general circulation model with idealised thermal forcing. Simulated transit spectra of GJ436b were generated using line-by-line radiative transfer models including the opacities of the molecular species expected to be present in such a planetary atmosphere. A new, ab-initio calculated, linelist for hot ammonia has been used for the first time. The photometric data observed at multiple wavelengths can be interpreted with methane being the dominant absorption after molecular hydrogen, possibly with minor contributions from ammonia, water and other molecules. No clear evidence of carbon monoxide and dioxide is found from transit photometry. We discuss this result in the light of a recent paper where photochemical disequilibrium is hypothesised to interpret secondary transit photometric data. We show that the emission photometric data are not incompatible with the presence of abundant methane, but further spectroscopic data are desirable to confirm this scenario.

We present two independent, homogeneous, global analyses of the transit light curves, radial velocities and spectroscopy of Kepler-4b through Kepler-8b, with numerous differences over the previous methods. These include: 1) improved decorrelated parameter fitting set, 2) new limb darkening coefficients, 3) time-stamps modified to BJD for consistency with RV, 4) two different methods for compensating for the integration-time of Kepler LC-data, 5) best-fit secondary-eclipse depths and excluded upper limits, 6) fitted mid-transit times, durations, depths and baseline fluxes for individual transits. We make several new determinations: 1) We detect a secondary eclipse for Kepler-7b of depth (47+/-14)ppm and significance 3.5-sigma. We conclude reflected light is a more plausible origin than thermal emission and determine a geometric albedo of Ag=(0.38+/-0.12). 2) An eccentric-orbit model for the Neptune-mass planet Kepler-4b is marginally detected to 2-sigma confidence with e=(0.25+/-0.12). If confirmed, this would place Kepler-4b in a similar category as GJ436b and HAT-P-11b as an eccentric, Neptune-mass planet. 3) We find marginal (2-sigma) evidence for Kepler-5b's secondary eclipse of depth (26+/-17)ppm. The most plausible explanation is reflected light caused by a planet with Ag=(0.15+/-0.10). 4) A 2.6-sigma peak in Kepler-6b TTV periodogram is detected and is not easily explained as an aliased frequency. We find that a resonant/non-resonant perturber, Trojan or exomoon all provide inadequate explanations and the most likely source is stellar rotation. 5) We find different impact parameters relative to the discovery papers in most cases, but generally self-consistent. 6) We constrain the presence of resonant planets, exomoons and Trojans using transit timing.

We report the discovery of HAT-P-24b, a transiting extrasolar planet orbiting the moderately bright V=11.818 F8 dwarf star GSC 0774-01441, with a period P = 3.3552464 +/- 0.0000071 d, transit epoch Tc = 2455216.97669 +/- 0.00024 (BJD_UTC), and transit duration 3.653 +/- 0.025 hours. The host star has a mass of 1.191 +/- 0.042 Msun, radius of 1.317 +/- 0.068 Rsun, effective temperature 6373 +/- 80 K, and a low metallicity of [Fe/H] = -0.16 +/- 0.08. The planetary companion has a mass of 0.681 +/- 0.031 MJ, and radius of 1.243 +/- 0.072 RJ yielding a mean density of 0.439 +/- 0.069 g cm-3 . By repeating our global fits with different parameter sets, we have performed a critical investigation of the fitting techniques used for previous HAT planetary discoveries. We find that the system properties are robust against the choice of priors. The effects of fixed versus fitted limb darkening are also examined. HAT-P-24b probably maintains a small eccentricity of e = 0.052 +0.022 -0.017, which is accepted over the circular orbit model with false alarm probability 5.8%. In the absence of eccentricity pumping, this result suggests HAT-P-24b experiences less tidal dissipation than Jupiter. Due to relatively rapid stellar rotation, we estimate that HAT-P-24b should exhibit one of the largest known Rossiter-McLaughlin effect amplitudes for an exoplanet (deltaVRM ~ 95 m/s) and thus a precise measurement of the sky-projected spin-orbit alignment should be possible.

We show that for a transiting exoplanet accompanied by a moon which also transits, the absolute masses and radii of the star, planet and moon are determinable. For a planet-star system, it is well known that the density of the star is calculable from the lightcurve by manipulation of Kepler's Third Law. In an analogous way, the planetary density is calculable for a planet-moon system which transits a star, and thus the ratio-of-densities is known. By combining this ratio with the observed ratio-of-radii and the radial velocity measurements of the system, we show that the absolute dimensions of the star and planet are determinable. This means such systems could be used as calibrators of stellar evolution. The detection of dynamical effects, such as transit timing variations, allows the absolute mass of the moon to be determined as well, which may be combined with the radius to infer the satellite's composition.

We report the discovery of HAT-P-26b, a transiting extrasolar planet orbiting the moderately bright V=11.744 K1 dwarf star GSC 0320-01027, with a period P = 4.234516 +- 0.000015 d, transit epoch Tc = 2455304.65122 +- 0.00035 (BJD), and transit duration 0.1023 +- 0.0010 d. The host star has a mass of 0.82 +- 0.03 Msun, radius of 0.79 + 0.10 - 0.04 Rsun, effective temperature 5079 +- 88 K, and metallicity [Fe/H] = -0.04 +- 0.08. The planetary companion has a mass of 0.059 +- 0.007 MJ, and radius of 0.565 + 0.072 - 0.032 RJ yielding a mean density of 0.40 +- 0.10 g cm-3. HAT-P-26b is the fourth Neptune-mass transiting planet discovered to date. It has a mass that is comparable to those of Neptune and Uranus, and slightly smaller than those of the other transiting Super-Neptunes, but a radius that is ~65% larger than those of Neptune and Uranus, and also larger than those of the other transiting Super-Neptunes. HAT-P-26b is consistent with theoretical models of an irradiated Neptune-mass planet with a 10 Mearth heavy element core that comprises >~ 50% of its mass with the remainder contained in a significant hydrogen-helium envelope, though the exact composition is uncertain as there are significant differences between various theoretical models at the Neptune-mass regime. The equatorial declination of the star makes it easily accessible to both Northern and Southern ground-based facilities for follow-up observations.

We report the discovery of four relatively massive (2-7MJ) transiting extrasolar planets. HAT-P-20b orbits a V=11.339 K3 dwarf star with a period P=2.875317+/-0.000004d. The host star has a mass of 0.760+/-0.03 Msun, radius of 0.690+/-0.02 Rsun, Teff=4595+/-80 K, and metallicity [Fe/H]=+0.35+/-0.08. HAT-P-20b has a mass of 7.246+/-0.187 MJ, and radius of 0.867+/-0.033 RJ yielding a mean density of 13.78+/-1.50 gcm^-3 , which is the second highest value among all known exoplanets. HAT-P-21b orbits a V=11.685 G3 dwarf on an eccentric (e=0.2280+/-0.016) orbit, with a period of P=4.1244810+/-000007d. The host star has a mass of 0.95+/-0.04Msun, radius of 1.10+/-0.08Rsun, Teff=5588+/-80K, and [Fe/H]=+0.01+/-0.08. HAT-P-21b has a mass of 4.063+/-0.161MJ, and radius of 1.024+/-0.092RJ. HAT-P-22b orbits the V=9.732 G5 dwarf HD233731, with P=3.2122200+/-0.000009d. The host star has a mass of 0.92+/-0.03Msun, radius of 1.04+/-0.04Rsun, Teff=5302+/-80K, and metallicity of +0.24+/-0.08. The planet has a mass of 2.147+/-0.061 MJ, and compact radius of 1.080+/-0.058RJ. The host star also harbors an M-dwarf companion at a wide separation. Finally, HAT-P-23b orbits a V=12.432 G0 dwarf star, with a period P=1.212884+/-0.000002d. The host star has a mass of 1.13+/-0.04sun, radius of 1.20+/-0.07Rsun, Teff=5905+/-80K, and [Fe/H]=+0.15+/-0.04. The planetary companion has a mass of 2.090+/-0.111MJ, and radius of 1.368+/-0.090RJ (abridged).

Out of the known transiting extrasolar planets, the majority are gas giants orbiting their host star at close proximity. Both theoretical and observational studies support the hypothesis that such bodies emit significant amounts of flux relative to the host star, increasing towards infrared wavelengths. For the dayside of the exoplanet, this phenomenon typically permits detectable secondary eclipses at such wavelengths, which may be used to infer atmospheric composition. In this paper, we explore the effects of emission from the nightside of the exoplanet on the primary transit lightcurve, which is essentially a self-blend. Allowing for nightside emission, an exoplanet's transit depth is no longer exclusively a function of the ratio-of-radii. The nightside of an exoplanet is emitting flux and the contrast to the star's emission is of the order of ~10^(-3) for hot-Jupiters. Consequently, we show that the transit depth in the mid-infrared will be attenuated due to flux contribution from the nightside emission by ~10^(-4). We show how this effect can be compensated for in the case where exoplanet phase curves have been measured, in particular for HD 189733b. For other systems, it may be possible to make a first-order correction by using temperature estimates of the planet. Unless the effect is accounted for, transmission spectra will also be polluted by nightside emission and we estimate that a Spitzer broadband spectrum on a bright target is altered at the 1-sigma level. Using archived Spitzer measurements, we show that the effect respectively increases the 8.0um and 24.0um transit depths by 1-sigma and 0.5-sigma per transit for HD 189733b. Consequently, we estimate that this would be 5-10 sigma effect for near-future JWST observations.

In this work, we investigate the accuracy of various approximate expressions for the transit duration of a detached binary against the exact solution, found through solving a quartic equation. Additionally, a new concise approximation is derived, which offers more accurate results than those currently in the literature. Numerical simulations are performed to test the accuracy of the various expressions. We find that our proposed expression show yields a >200% improvement in accuracy relative to the most previously employed expression. We derive a new set of equations for retrieving the lightcurve parameters and consider the effect of falsely using circular expressions for eccentric orbits, with particularly important consequences for transit surveys. The new expression also allows us to propose a new lightcurve fitting parameter set, which minimizes the mutual correlations and thus improves computational efficiency. The equation is also readily differentiated to provide analytic expressions for the transit duration variation (TDV) due to secular variations in the system parameters, for example due to apsidal precession induced by perturbing planets.

We report the discovery of HAT-P-14b, a fairly massive transiting extrasolar planet orbiting the moderately bright star GSC 3086-00152 (V = 9.98), with a period of P = 4.627669 +/- 0.000005 days. The transit is close to grazing (impact parameter 0.891 +0.007/-0.008) and has a duration of 0.0912 +/- 0.0017 days, with a reference epoch of mid transit of Tc = 2454875.28938 +/- 0.00047 (BJD). The orbit is slightly eccentric (e = 0.107 +/- 0.013), and the orientation is such that occultations are unlikely to occur. The host star is a slightly evolved mid-F dwarf with a mass of 1.386 +/- 0.045 M(Sun), a radius of 1.468 +/- 0.054 R(Sun) effective temperature 6600 +/- 90 K, and a slightly metal-rich composition corresponding to [Fe/H] = +0.11 +/- 0.08. The planet has a mass of 2.232 +/- 0.059 M(Jup) and a radius of 1.150 +/- 0.052 R(Jup), implying a mean density of 1.82 +/- 0.24 g/cm3. Its radius is well reproduced by theoretical models for the 1.3 Gyr age of the system if the planet has a heavy-element fraction of about 50 M(Earth) (7% of its total mass). The brightness, near-grazing orientation, and other properties of HAT-P-14 make it a favorable transiting system to look for changes in the orbital elements or transit timing variations induced by a possible second planet, and also to place meaningful constraints on the presence of sub-Earth mass or Earth mass exomoons, by monitoring it for transit duration variations.

The hot Jupiter HD 209458b was observed during primary transit at 3.6, 4.5, 5.8 and 8.0 microns using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We detail here the procedures we adopted to correct for the systematic trends present in the IRAC data. The light curves were fitted including limb darkening effects and fitted using Markov Chain Monte Carlo and prayer-bead Monte Carlo techniques, finding almost identical results. The final depth measurements obtained by a combined Markov Chain Monte Carlo fit are at 3.6 microns, 1.469 +- 0.013 % and 1.448 +- 0.013 %; at 4.5 microns, 1.478 +- 0.017 % ; at 5.8 microns, 1.549 +- 0.015 % and at 8.0 microns 1.535 +- 0.011 %. Our results clearly indicate the presence of water in the planetary atmosphere. Our broad band photometric measurements with IRAC prevent us from determining the additional presence of other other molecules such as CO, CO2 and methane for which spectroscopy is needed. While water vapour with a mixing ratio of 10^-4-10^-3 combined with thermal profiles retrieved from the day-side may provide a very good fit to our observations, this data set alone is unable to resolve completely the degeneracy between water abundance and atmospheric thermal profile.

The search for life outside of the Solar System should not be restricted to exclusively planetary bodies; large moons of extrasolar planets may also be common habitable environments throughout the Galaxy. Extrasolar moons, or exomoons, may be detected through transit timing effects induced onto the host planet as a result of mutual gravitational interaction. In particular, transit timing variations (TTV) and transit duration variations (TDV) are predicted to produce a unique exomoon signature, which is not only easily distinguished from other gravitational perturbations, but also provides both the period and mass of an exomoon. Using these timing effects, photometry greater or equal to that of the Kepler Mission is readily able to detect habitable-zone exomoons down to 0.2 Earth masses and could survey up to 25,000 stars for 1 Earth-mass satellites. We discuss future possibilities for spectral retrieval of such bodies and show that transmission spectroscopy with JWST should be able to detect molecular species with ~30 transit events, in the best cases.

In this paper we investigate the detectability of a habitable-zone exomoon around various configurations of exoplanetary systems with the Kepler Mission or photometry of approximately equal quality. We calculate both the predicted transit timing signal amplitudes and the estimated uncertainty on such measurements in order to calculate the confidence in detecting such bodies across a broad spectrum of orbital arrangements. The effects of stellar variability, instrument noise and photon noise are all accounted for in the analysis. We validate our methodology by simulating synthetic lightcurves and performing a Monte Carlo analysis for several cases of interest. We find that habitable-zone exomoons down to 0.2 Earth masses may be detected and ~25,000 stars could be surveyed for habitable-zone exomoons within Kepler's field-of-view. A Galactic Plane survey with Kepler-class photometry could potentially survey over one million stars for habitable-zone exomoons. In conclusion, we propose that habitable exomoons will be detectable should they exist in the local part of the galaxy.

In our previous paper, we evaluated the transit duration variation (TDV) effect for a co-aligned planet-moon system at an orbital inclination of i=90 degrees. Here, we will consider the effect for the more general case of i <= 90 degrees and an exomoon inclined from the planet-star plane by Euler rotation angles $\alpha$, $\beta$ and $\gamma$. We find that the TDV signal has two major components, one due to the velocity variation effect described in our first paper and one new component due to transit impact parameter variation. By evaluating the dominant terms, we find the two effects are additive for prograde exomoon orbits, and deductive for retrograde orbits. This asymmetry could allow for future determination of the orbital sense of motion. We re-evaluate the ratio of TDV and TTV effects, $\eta$, in the more general case of an inclined planetary orbit with a circular orbiting moon and find that it is still possible to directly determine the moon's orbital separation from just the ratio of the two amplitudes, as first proposed in our previous paper.

We report the detection of a transit egress by the ~ 3.9-Jupiter-mass planet HD 80606b, an object in a highly-eccentric orbit (e ~ 0.93) about its parent star of approximately solar type. The astrophysical reality of the signal of variability in HD 80606 is confirmed by observation with two independent telescope systems, and checks against several reference stars in the field. Differential photometry with respect to the nearby comparison star HD 80607 provides a precise light curve. Modelling of the light curve with a full eccentric-orbit model indicates a planet/star-radius ratio of 0.1057 +/- 0.0018, corresponding to a planet radius of 1.029 R_J for a solar-radius parent star; and a precise orbital inclination of 89.285 +/- 0.023 degrees, giving a total transit duration of 12.1 +/- 0.4 hours. The planet hence joins HD 17156b in a class of highly eccentric transiting planets, in which HD 80606b has both the longest period and most eccentric orbit. The recently reported discovery of a secondary eclipse of HD 80606b by the Spitzer Space Observatory permits a combined analysis with the mid-time of primary transit in which the orbital parameters of the system can be tightly constrained. We derive a transit ephemeris of T_tr = HJD (2454876.344 +/- 0.011) + (111.4277 +/- 0.0032) E.

As the number of known exoplanets continues to grow, the question as to whether such bodies harbour satellite systems has become one of increasing interest. In this paper, we explore the transit timing effects that should be detectable due to an exomoon and predict a new observable. We first consider transit time variation (TTV), where we update the model to include the effects of orbital eccentricity. We draw two key conclusions: 1) In order to maintain Hill stability, the orbital frequency of the exomoon will always be higher than the sampling frequency. Therefore, the period of the exomoon cannot be reliably determined from TTV, only a set of harmonic frequencies. 2) The TTV amplitude is proportional to M_S a_S where M_S is the exomoon mass and a_S is the semi-major axis of the moon's orbit. Therefore, M_S and a_S cannot be separately determined. We go on to predict a new observable due to exomoons - transit duration variation (TDV). We derive the TDV amplitude and conclude that its amplitude is not only detectable, but the TDV signal will provide two robust advantages: 1) The TDV amplitude is proportional to M_S a_S^{-1/2} and therefore the ratio of TTV and TDV allows for a separate determination of M_S and a_S. 2) TDV has a 90 degrees phase difference to the TTV signal, making it an excellent complementary technique.

Transiting planet lightcurves have historically been used predominantly for measuring the depth and hence ratio of the planet-star radii, p. Equations have been previously presented by Seager & Mallen-Ornelas (2003) for the analysis of the total and trough transit lightcurve times to derive the ratio of semi-major axis to stellar radius, a/R*, in the case of circular orbits. Here, a new analytic model is proposed which operates for the more general case of an eccentric orbit. We aim to investigate three major effects our model predicts: i) the degeneracy in transit lightcurve solutions for eccentricity, e>0 ii) the asymmetry of the lightcurve and the resulting shift in the mid-transit time, Tmid iii) the effect of eccentricity on the ingress and egress slopes. It is shown that a system with changing eccentricity and inclination may produce a long period transit time variation (LTTV). Furthermore, we use our model in a reanalysis of HD 209458 b archived data by Richardson et al. (2006), where we include the confirmed non-zero eccentricity and derive a 24 micron planetary radius of R_P = 1.275 +- 0.082 R_J (where R_J = 1 Jovian radius), which is 1% larger than is we assume a circular orbit.