Showing 1–50 of 59 results for author: Kipping, D M

LcTools: A Windows-Based Software System for Finding and Recording Signals in Lightcurves from NASA Space Missions

Authors: Allan R. Schmitt, Joel D. Hartman, David M. Kipping

Abstract: 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… ▽ More 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. △ Less

Submitted 17 October, 2019; originally announced October 2019.

Comments: 28 pages, 31 figures, arXiv-only submission

Evidence for a Large Exomoon Orbiting Kepler-1625b

Authors: Alex Teachey, David M. Kipping

Abstract: 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 t… ▽ More 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. △ Less

Submitted 4 October, 2018; originally announced October 2018.

Comments: 98 pages, 22 figures, 5 tables. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC). Select data products available at https://academiccommons.columbia.edu/doi/10.7916/D8795NHS

Journal ref: Science Advances, Vol. 4, no. 10, eaav1784 (October 3 2018)

On the detectability of transiting planets orbiting white dwarfs using LSST

Authors: Jorge Cortes, David M. Kipping

Abstract: 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 sh… ▽ More 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%. △ Less

Submitted 1 October, 2018; originally announced October 2018.

Comments: 9 pages, 6 figures, Submitted for publication in MNRAS on September 6, 2018

A second terrestrial planet orbiting the nearby M dwarf LHS 1140

Authors: Kristo Ment, Jason A. Dittmann, Nicola Astudillo-Defru, David Charbonneau, Jonathan Irwin, Xavier Bonfils, Felipe Murgas, Jose-Manuel Almenara, Thierry Forveille, Eric Agol, Sarah Ballard, Zachory K. Berta-Thompson, François Bouchy, Ryan Cloutier, Xavier Delfosse, René Doyon, Courtney D. Dressing, Gilbert A. Esquerdo, Raphaëlle D. Haywood, David M. Kipping, David W. Latham, Christophe Lovis, Elisabeth R. Newton, Francesco Pepe, Joseph E. Rodriguez , et al. (5 additional authors not shown)

Abstract: 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,… ▽ More 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. △ Less

Submitted 21 November, 2018; v1 submitted 1 August, 2018; originally announced August 2018.

Comments: 22 pages, 8 figures, accepted by AJ

Jupiter Analogues Orbit Stars with an Average Metallicity Close to that of the Sun

Authors: Lars A. Buchhave, Bertram Bitsch, Anders Johansen, David W. Latham, Martin Bizzarro, Allyson Bieryla, David M. Kipping

Abstract: 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 a… ▽ More 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. △ Less

Submitted 19 February, 2018; originally announced February 2018.

Comments: Accepted in ApJ

Finding Mountains with Molehills: The Detectability of Exotopography

Authors: Moiya A. S. McTier, David M. Kipping

Abstract: 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. T… ▽ More 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. △ Less

Submitted 17 January, 2018; originally announced January 2018.

Comments: Accepted to MNRAS

HEK VI: On the Dearth of Galilean Analogs in Kepler and the Exomoon Candidate Kepler-1625b I

Authors: Alex Teachey, David M. Kipping, Allan R. Schmitt

Abstract: 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… ▽ More 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 $ηいーた<0.38$ to 95% confidence for the 284 KOIs considered, with a 68.3% confidence interval of $ηいーた=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. △ Less

Submitted 18 October, 2017; v1 submitted 26 July, 2017; originally announced July 2017.

Comments: Accepted to The Astronomical Journal October 14, 2017. TTV posteriors are available at https://github.com/alexteachey/TTV_posteriors

Relativistic Light Sails

Authors: David M. Kipping

Abstract: 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 vel… ▽ More 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. △ Less

Submitted 22 November, 2017; v1 submitted 13 April, 2017; originally announced April 2017.

Comments: Accepted to AJ. This version corrects the comparison in Figure 4 between our prediction and that of previous works, by accounting for finite light travel time effects. Other results are unchanged

Journal ref: AJ, 153, 277 (2017)

Masses of Kepler-46b, c from Transit Timing Variations

Authors: Ximena Saad-Olivera, David Nesvorný, David M. Kipping, Fernando Roig

Abstract: 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… ▽ More 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. △ Less

Submitted 5 April, 2017; originally announced April 2017.

Comments: 12 pages, 6 figures

Journal ref: The Astronomical Journal, volume 153, number 4, page 198, 2017

Transit Clairvoyance: Enhancing TESS follow-up using artificial neural networks

Authors: David M. Kipping, Christopher Lam

Abstract: 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)… ▽ More 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. △ Less

Submitted 15 November, 2016; originally announced November 2016.

Comments: Accepted to MNRAS. Grid of results and interpolating code available at https://github.com/cl3425/Clairvoyance

No Conclusive Evidence for Transits of Proxima b in MOST photometry

Authors: David M. Kipping, Chris Cameron, Joel D. Hartman, James R. A. Davenport, Jaymie M. Matthews, Dimitar Sasselov, Jason Rowe, Robert J. Siverd, Jingjing Chen, Emily Sandford, Gáspár Á. Bakos, Andres Jordan, Daniel Bayliss, Thomas Henning, Luigi Mancini, Kaloyan Penev, Zoltan Csubry, Waqas Bhatti, Joao Da Silva Bento, David B. Guenther, Rainer Kuschnig, Anthony F. J. Moffat, Slavek M. Rucinski, Werner W. Weiss

Abstract: The analysis of Proxima Centauri's radial velocities recently led Anglada-Escudé 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 observ… ▽ More The analysis of Proxima Centauri's radial velocities recently led Anglada-Escudé 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. △ Less

Submitted 1 December, 2016; v1 submitted 27 September, 2016; originally announced September 2016.

Comments: Accepted to ApJ. Posterior samples, MOST photometry and HATSouth photometry are all available at https://github.com/CoolWorlds/Proxima

MOST Observations of our Nearest Neighbor: Flares on Proxima Centauri

Authors: James R. A. Davenport, David M. Kipping, Dimitar Sasselov, Jaymie M. Matthews, Chris Cameron

Abstract: 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… ▽ More 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. △ Less

Submitted 28 September, 2016; v1 submitted 23 August, 2016; originally announced August 2016.

Comments: 6 pages, 4 figures, ApJ Letters accepted

A Cloaking Device for Transiting Planets

Authors: David M. Kipping, Alex Teachey

Abstract: 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 deliberate… ▽ More 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. △ Less

Submitted 29 March, 2016; originally announced March 2016.

Comments: 10 pages, 6 figures, accepted in MNRAS

Probabilistic Forecasting of the Masses and Radii of Other Worlds

Authors: Jingjing Chen, David M. Kipping

Abstract: 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… ▽ More 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. △ Less

Submitted 3 November, 2016; v1 submitted 28 March, 2016; originally announced March 2016.

Comments: Accepted in ApJ. New robustness test added, results unchanged. Forecaster code is available at https://github.com/chenjj2/forecaster

Observational biases for transiting planets

Authors: David M. Kipping, Emily Sandford

Abstract: 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 ef… ▽ More 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. △ Less

Submitted 17 August, 2016; v1 submitted 17 March, 2016; originally announced March 2016.

Comments: 10 pages, 9 figures, accepted for publication in MNRAS

A Transiting Jupiter Analog

Authors: David M. Kipping, Guillermo Torres, Chris Henze, Alex Teachey, Howard T. Isaacson, Erik A. Petigura, Geoffrey W. Marcy, Lars A. Buchhave, Jingjing Chen, Steve T. Bryson, Emily Sandford

Abstract: 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 met… ▽ More 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. △ Less

Submitted 29 February, 2016; originally announced March 2016.

Comments: 14 pages, 10 figures. Accepted to ApJ. Posteriors available at https://github.com/CoolWorlds/Kepler-167-Posteriors

Efficient, uninformative sampling of limb-darkening coefficients for a three-parameter law

Authors: David M. Kipping

Abstract: 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-darkenin… ▽ More 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. △ Less

Submitted 30 October, 2015; v1 submitted 11 September, 2015; originally announced September 2015.

Comments: 12 pages, 8 figures. Accepted to MNRAS. Code and LaTeX source at https://github.com/davidkipping/LDC3

The Hunt for Exomoons with Kepler (HEK): V. A Survey of 41 Planetary Candidates for Exomoons

Authors: David M. Kipping, Allan R. Schmitt, Chelsea X. Huang, Guillermo Torres, David Nesvorny, Lars A. Buchhave, Joel Hartman, Gáspár Á. Bakos

Abstract: 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 exh… ▽ More 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. △ Less

Submitted 24 September, 2015; v1 submitted 18 March, 2015; originally announced March 2015.

Comments: 18 pages, 9 figures, 4 tables. Accepted in ApJ

Validation of Twelve Small Kepler Transiting Planets in the Habitable Zone

Authors: Guillermo Torres, David M. Kipping, Francois Fressin, Douglas A. Caldwell, Joseph D. Twicken, Sarah Ballard, Natalie M. Batalha, Stephen T. Bryson, David R. Ciardi, Christopher E. Henze, Steve B. Howell, Howard T. Isaacson, Jon M. Jenkins, Philip S. Muirhead, Elisabeth R. Newton, Erik A. Petigura, Thomas Barclay, William J. Borucki, Justin R. Crepp, Mark E. Everett, Elliott P. Horch, Andrew W. Howard, Rea Kolbl, Geoffrey W. Marcy, Sean McCauliff , et al. (1 additional authors not shown)

Abstract: 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 regi… ▽ More 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. △ Less

Submitted 7 January, 2015; v1 submitted 6 January, 2015; originally announced January 2015.

Comments: 27 pages in emulateapj format, including tables and figures. To appear in The Astrophysical Journal

The Possible Moon of Kepler-90g is a False Positive

Authors: David M. Kipping, Chelsea X. Huang, David Nesvorny, Guillermo Torres, Lars A. Buchhave, Gáspár Á. Bakos, Allan R. Schmitt

Abstract: 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 Kepl… ▽ More 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. △ Less

Submitted 25 November, 2014; originally announced November 2014.

Comments: 5 pages, 4 figures. Accepted for publication in The Astrophysical Journal Letters

Kepler-432: a red giant interacting with one of its two long period giant planets

Authors: Samuel N. Quinn, Timothy R. White, David W. Latham, William J. Chaplin, Rasmus Handberg, Daniel Huber, David M. Kipping, Matthew J. Payne, Chen Jiang, Victor Silva Aguirre, Dennis Stello, David H. Sliski, David R. Ciardi, Lars A. Buchhave, Timothy R. Bedding, Guy R. Davies, Saskia Hekker, Hans Kjeldsen, Mark E. Everett, Steve B. Howell, Sarbani Basu, Tiago L. Campante, Jørgen Christensen-Dalsgaard, Yvonne P. Elsworth, Christoffer Karoff , et al. (6 additional authors not shown)

Abstract: 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… ▽ More 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えーゆー. △ Less

Submitted 20 April, 2015; v1 submitted 17 November, 2014; originally announced November 2014.

Comments: 22 pages, 19 figures, 5 tables. Accepted to ApJ on Jan 24, 2015 (submitted Nov 11, 2014). Updated with minor changes to match published version

Journal ref: 2015, ApJ, 803, 49

Bayesian priors for the eccentricity of transiting planets

Authors: David M. Kipping

Abstract: 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 pro… ▽ More 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. △ Less

Submitted 6 August, 2014; originally announced August 2014.

Comments: 9 pages, 5 figures. Accepted to MNRAS. Code available at http://www.cfa.harvard.edu/~dkipping/ECCSAMPLES.html

Discovery of a Transiting Planet Near the Snow-Line

Authors: David M. Kipping, Guillermo Torres, Lars A. Buchhave, Scott J. Kenyon, Christopher E. Henze, Howard Isaacson, Rea Kolbl, Geoff W. Marcy, Stephen T. Bryson, Keivan G. Stassun, Fabienne A. Bastien

Abstract: 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 dw… ▽ More 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. △ Less

Submitted 17 July, 2014; originally announced July 2014.

Comments: 14 pages, 10 figures, 3 tables. Accepted in ApJ

spotrod: a semi-analytic model for transits of spotted stars

Authors: Bence Béky, David M. Kipping, Matthew J. Holman

Abstract: 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… ▽ More 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. △ Less

Submitted 16 July, 2014; originally announced July 2014.

Comments: 16 pages, 15 figures. Accepted for publication in MNRAS volume 442 issue 4

In Search of Exomoons

Authors: David M. Kipping

Abstract: 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 det… ▽ More 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. △ Less

Submitted 6 May, 2014; originally announced May 2014.

Comments: 18 pages, 5 figures; to appear in the proceedings for the Frank N. Bash Symposium 2013: New Horizons in Astronomy, held October 6-8, 2013 in Austin, TX

Stellar rotation--planetary orbit period commensurability in the HAT-P-11 system

Authors: Bence Béky, Matthew J. Holman, David M. Kipping, Robert W. Noyes

Abstract: 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 trans… ▽ More 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. △ Less

Submitted 28 March, 2014; originally announced March 2014.

Comments: Accepted for publication in ApJ. 10 pages, 7 figures

Flicker as a tool for characterizing planets through Asterodensity Profiling

Authors: David M. Kipping, Fabienne A. Bastien, Keivan G. Stassun, William J. Chaplin, Daniel Huber, Lars A. Buchhave

Abstract: 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 tra… ▽ More 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}(ρろー_{\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. △ Less

Submitted 20 March, 2014; originally announced March 2014.

Comments: 6 pages, 3 figures, 1 table. Accepted to ApJ Letters. Code available at https://www.cfa.harvard.edu/~dkipping/flicker.html

Journal ref: Astrophys.J.Let.785:L32-L36,2014

The Hunt for Exomoons with Kepler (HEK): IV. A Search for Moons around Eight M-Dwarfs

Authors: David M. Kipping, David Nesvorny, Lars A. Buchhave, Joel Hartman, Gáspár Á. Bakos, Allan R. Schmitt

Abstract: 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 e… ▽ More 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. △ Less

Submitted 26 January, 2014; v1 submitted 6 January, 2014; originally announced January 2014.

Comments: 15 pages, 13 figures, 6 tables. Accepted in ApJ

Journal ref: Astrophys.J.784:28-41,2014

A High False Positive Rate for Kepler Planetary Candidates of Giant Stars using Asterodensity Profiling

Authors: David H. Sliski, David M. Kipping

Abstract: 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 tran… ▽ More 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$σしぐま$ 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$%. △ Less

Submitted 27 May, 2014; v1 submitted 6 January, 2014; originally announced January 2014.

Comments: 11 pages, 5 figures, 2 tables. Accepted in ApJ

Characterizing Distant Worlds with Asterodensity Profiling

Authors: David M. Kipping

Abstract: 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,… ▽ More 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. △ Less

Submitted 9 June, 2014; v1 submitted 5 November, 2013; originally announced November 2013.

Comments: 27 pages, 11 figures, accepted to MNRAS (fixed broken figures). Hi-res cheat-sheet available at http://www.cfa.harvard.edu/~dkipping/AP_Cheat_Sheet.pdf

Journal ref: Mon. Not. R. Astron. Soc., 2014, Volume 440, Issue 3, pp. 2164-2184

Efficient, uninformative sampling of limb darkening coefficients for two-parameter laws

Authors: David M. Kipping

Abstract: 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 v… ▽ More 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. △ Less

Submitted 14 August, 2013; v1 submitted 31 July, 2013; originally announced August 2013.

Comments: 11 pages, 5 figures, 2 tables, accepted in MNRAS

Journal ref: Mon. Not. R. Astron. Soc., 2013, Volume 435, Issue 3, pp. 2152-2160

Parametrizing the exoplanet eccentricity distribution with the Beta distribution

Authors: David M. Kipping

Abstract: 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 f… ▽ More 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. △ Less

Submitted 26 June, 2013; v1 submitted 20 June, 2013; originally announced June 2013.

Comments: 5 pages, 4 figures, 2 tables, accepted in MNRAS Letters. Code available at https://www.cfa.harvard.edu/~dkipping/betaprior.html

Journal ref: Mon. Not. R. Astron. Soc., 2013, Volume 434, Issue 1, pp. L51-L55

A simple, quantitative method to infer the minimum atmospheric height of small exoplanets

Authors: David M. Kipping, David S. Spiegel, Dimitar D. Sasselov

Abstract: 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.… ▽ More 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. △ Less

Submitted 25 June, 2013; v1 submitted 13 June, 2013; originally announced June 2013.

Comments: 7 pages, 3 figures, 3 tables, accepted in MNRAS. Code available at http://www.cfa.harvard.edu/~dkipping/mah.html

Journal ref: Mon. Not. R. Astron. Soc., 2013, Volume 434, Issue 3, pp. 1883-1888

The Hunt for Exomoons with Kepler (HEK): III. The First Search for an Exomoon around a Habitable-Zone Planet

Authors: David M. Kipping, Duncan Forgan, Joel Hartman, David Nesvorny, Gáspár Á. Bakos, Allan R. Schmitt, Lars A. Buchhave

Abstract: 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 avai… ▽ More 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. △ Less

Submitted 9 September, 2013; v1 submitted 6 June, 2013; originally announced June 2013.

Comments: 19 pages, 11 figures, 7 tables. Accepted in ApJ. Planet-moon transit animations available at https://www.cfa.harvard.edu/~dkipping/kepler22.html

Journal ref: Astrophys.J.777:134-150,2013

The Hunt for Exomoons with Kepler (HEK): II. Analysis of Seven Viable Satellite-Hosting Planet Candidates

Authors: David M. Kipping, Joel Hartman, Lars A. Buchhave, Allan R. Schmitt, Gáspár Á. Bakos, David Nesvorny

Abstract: 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 exhibi… ▽ More 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. △ Less

Submitted 6 March, 2013; v1 submitted 9 January, 2013; originally announced January 2013.

Comments: 32 pages, 11 figures, 23 tables, Accepted to ApJ

Journal ref: Astrophys.J.770:101-132,2013

An Analytic Model for Rotational Modulations in the Photometry of Spotted Stars

Authors: David M. Kipping

Abstract: 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… ▽ More 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. △ Less

Submitted 1 November, 2012; v1 submitted 13 September, 2012; originally announced September 2012.

Comments: 30 pages, 6 figures, 2 tables, accepted in MNRAS. Error corrected in transit depth variations section. Code available at http://www.cfa.harvard.edu/~dkipping/macula.html

Journal ref: Mon. Not. R. Astron. Soc., 2012, Volume 427, Issue 3, pp. 2487-2511

The Detection and Characterization of a Nontransiting Planet by Transit Timing Variations

Authors: David Nesvorny, David M. Kipping, Lars A. Buchhave, Gáspár Á. Bakos, Joel Hartman, Allan Schmitt

Abstract: 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 timi… ▽ More 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. △ Less

Submitted 4 August, 2012; originally announced August 2012.

Comments: Accepted for publication in Science. Published online on May 10, 2012. Main Text and supplemental information included in a single merged file, 58 pages

Journal ref: Science, Volume 336, Issue 6085, pp. 1133-1136 (2012)

The Hunt for Exomoons with Kepler (HEK): I. Description of a New Observational Project

Authors: David M. Kipping, Gáspár Á. Bakos, Lars A. Buchhave, David Nesvorny, Allan Schmitt

Abstract: 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 habit… ▽ More 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. △ Less

Submitted 15 March, 2012; v1 submitted 3 January, 2012; originally announced January 2012.

Comments: 21 pages, 8 figures, 4 tables, accepted in ApJ

Journal ref: Astrophys.J.750:115-134,2012

A Novel Method to Photometrically Constrain Orbital Eccentricities: Multibody Asterodensity Profiling (MAP)

Authors: David M. Kipping, William R. Dunn, Jamie M. Jasinski, Varun P. Manthri

Abstract: 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 t… ▽ More 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. △ Less

Submitted 12 December, 2011; originally announced December 2011.

Comments: 25 pages, 14 figures. Accepted in MNRAS

Journal ref: Mon. Not. R. Astron. Soc., 2012, Volume 421, Issue 2, pp. 1166-1188

Detection of visible light from the darkest world

Authors: David M. Kipping, David S. Spiegel

Abstract: 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 si… ▽ More 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. △ Less

Submitted 16 August, 2011; v1 submitted 10 August, 2011; originally announced August 2011.

Comments: 5 pages, 3 figures. Accepted in MNRAS Letters. Minor typos corrected from last version

Journal ref: Mon. Not. R. Astron. Soc., 2011, Volume 417, Issue 1, pp. L88-L92

HAT-P-31b,c: A Transiting, Eccentric, Hot Jupiter and a Long-Period, Massive Third-Body

Authors: David M. Kipping, Joel Hartman, Gáspár Á. Bakos, Guillermo Torres, David W. Latham, Daniel Bayliss, László L. Kiss, Bun'ei Sato, Bence Béky, Géza Kovács, Sam N. Quinn, Lars A. Buchhave, Jens Andersen, Geoff W. Marcy, Andrew W. Howard, Debra A. Fischer, John A. Johnson, Robert W. Noyes, Dimitar D. Sasselov, Robert P. Stefanik, József Lázár, István Papp, Pál Sári, Gabor Furesz

Abstract: 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, determ… ▽ More 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. △ Less

Submitted 7 July, 2011; v1 submitted 6 June, 2011; originally announced June 2011.

Comments: 12 pages, 7 figures, 7 tables. Accepted in AJ. Minor changes over previous version

Journal ref: Astronomical Journal 142 95, 2011

LUNA: An algorithm for generating dynamic planet-moon transits

Authors: David M. Kipping

Abstract: 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 on… ▽ More 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. △ Less

Submitted 23 April, 2012; v1 submitted 17 May, 2011; originally announced May 2011.

Comments: Accepted in MNRAS, 2011 May 16. Minor typos corrected (thanks to S. Awiphan)

Journal ref: Mon. Not. R. Astron. Soc., 2011, Volume 416, Issue 1, pp. 689-709

The Transits of Extrasolar Planets with Moons

Authors: David M. Kipping

Abstract: 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… ▽ More 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. △ Less

Submitted 23 April, 2012; v1 submitted 16 May, 2011; originally announced May 2011.

Comments: PhD thesis, University College London (March 14, 2011), 211 pages, 42 figures, 10 tables. Further minor typos corrected in revised version

How to Weigh a Star Using a Moon

Authors: David M. Kipping

Abstract: 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 transit… ▽ More 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. △ Less

Submitted 12 October, 2010; originally announced October 2010.

Comments: 5 pages; Accepted in MNRAS

Journal ref: Mon. Not. R. Astron. Soc., 2010, Volume 409, Issue 1, pp. L119-L123

HAT-P-26b: A Low-Density Neptune-Mass Planet Transiting a K Star

Authors: J. D. Hartman, G. Á. Bakos, D. M. Kipping, G. Torres, G. Kovács, R. W. Noyes, D. W. Latham, A. W. Howard, D. A. Fischer, J. A. Johnson, G. W. Marcy, H. Isaacson, S. N. Quinn, L. A. Buchhave, B. Béky, D. D. Sasselov, R. P. Stefanik, G. A. Esquerdo, M. Everett, G. Perumpilly, J. Lázár, I. Papp, P. Sári

Abstract: 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 metal… ▽ More 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. △ Less

Submitted 5 October, 2010; originally announced October 2010.

Comments: 16 pages, 9 figures, 5 tables, submitted to ApJ

HAT-P-24b: An inflated hot-Jupiter on a 3.36d period transiting a hot, metal-poor star

Authors: D. M. Kipping, G. A. Bakos, J. Hartman, G. Torres, A. Shporer, D. W. Latham, Geza Kovacs, R. W. Noyes, A. W. Howard, D. A. Fischer, J. A. Johnson, G. W. Marcy, B. Beky, G. Perumpilly, G. A. Esquerdo, D. D. Sasselov, R. P. Stefanik, J. Lazar, I. Papp, P. Sari

Abstract: 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… ▽ More 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. △ Less

Submitted 13 October, 2010; v1 submitted 19 August, 2010; originally announced August 2010.

Comments: 13 pages with 4 figures and 8 tables in emulateapj format. Minor changes. Accepted in The Astrophysical Journal

Journal ref: Astrophys.J.725:2017-2028,2011

Methane in the atmosphere of the transiting hot Neptune GJ436b?

Authors: J. -P. Beaulieu, G. Tinetti, D. M. Kipping, I. Ribas, R. J. Barber, J. Y-K. Cho, I. Polichtchouk, J. Tennyson, S. N. Yurchenko, C. A. Griffith, V. Batista, I. Waldmann, S. Miller, S. Carey, O. Mousis, S. J. Fossey, A. Aylward

Abstract: We present an analysis of seven primary transit observations of the hot Neptune GJ436b at 3.6, 4.5 and $8~μみゅー$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 obser… ▽ More We present an analysis of seven primary transit observations of the hot Neptune GJ436b at 3.6, 4.5 and $8~μみゅー$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~μみゅー$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. △ Less

Submitted 14 February, 2011; v1 submitted 2 July, 2010; originally announced July 2010.

Comments: 19 pages, 10 figures, 1 table, Astrophysical Journal in press

Analysis of Kepler's short-cadence photometry for TrES-2b

Authors: David M. Kipping, Gáspár Á. Bakos

Abstract: 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… ▽ More 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. △ Less

Submitted 10 March, 2011; v1 submitted 29 June, 2010; originally announced June 2010.

Comments: Accepted in ApJ. 19 pages, 12 figures, 5 tables, emulateapj

Journal ref: Astrophys.J.733:36-52,2011

Investigations of approximate expressions for the transit duration

Authors: David M. Kipping

Abstract: 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 expre… ▽ More 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. △ Less

Submitted 21 April, 2010; originally announced April 2010.

Comments: Accepted in MNRAS

Journal ref: Mon. Not. R. Astron. Soc., 2010, Volume 407, Issue 1, pp. 301-313

Binning is sinning: morphological light-curve distortions due to finite integration time

Authors: David M. Kipping

Abstract: 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-c… ▽ More 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. △ Less

Submitted 14 April, 2014; v1 submitted 21 April, 2010; originally announced April 2010.

Comments: Long delayed upload of the MNRAS accepted version, 10 pages, 3 figures

Journal ref: Mon. Not. R. Astron. Soc., 2010, Volume 408, Issue 3, pp. 1758-1769