Researcher profile

David M. Kipping

David M. Kipping contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
39works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

39 published item(s)

preprint2016arXiv

A Cloaking Device for Transiting Planets

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.

preprint2016arXiv

A Transiting Jupiter Analog

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.

preprint2016arXiv

Observational biases for transiting planets

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.

preprint2016arXiv

Probabilistic Forecasting of the Masses and Radii of Other Worlds

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

preprint2015arXiv

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

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.

preprint2015arXiv

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

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.

preprint2015arXiv

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

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.

preprint2015arXiv

Validation of Twelve Small Kepler Transiting Planets in the Habitable Zone

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.

preprint2014arXiv

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

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$%.

preprint2014arXiv

Bayesian priors for the eccentricity of transiting planets

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.

preprint2014arXiv

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

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.

preprint2014arXiv

Characterizing Distant Worlds with Asterodensity Profiling

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.

preprint2014arXiv

Discovery of a Transiting Planet Near the Snow-Line

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.

preprint2014arXiv

Flicker as a tool for characterizing planets through Asterodensity Profiling

Variability in the time series brightness of a star on a timescale of 8 hours, known as &#39;flicker&#39;, 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&#39;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.

preprint2014arXiv

In Search of Exomoons

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 &#34;exomoons&#34;. 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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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.

preprint2014arXiv

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

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 &#39;Hunt for Exomoons with Kepler&#39; (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&#39;s pixels, which we confirm using a novel way of examining pixel-level data which we dub the &#39;transit centroid&#39;. 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.

preprint2013arXiv

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

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 &#39;minimum atmospheric height&#39; (R_MAH) using only a planet&#39;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.

preprint2013arXiv

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

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.

preprint2013arXiv

Parametrizing the exoplanet eccentricity distribution with the Beta distribution

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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 &#39;Hunt for Exomoons with Kepler&#39; (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.

preprint2012arXiv

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

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.

preprint2012arXiv

LUNA: An algorithm for generating dynamic planet-moon transits

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.

preprint2012arXiv

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

The Kepler Mission is monitoring the brightness of ~150,000 stars searching for evidence of planetary transits. As part of the &#34;Hunt for Exomoons with Kepler&#34; (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.

preprint2012arXiv

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

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 &#34;The Hunt for Exomoons with Kepler&#34; (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.

preprint2012arXiv

The Transits of Extrasolar Planets with Moons

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 &#34;exomoons&#34; 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.

preprint2011arXiv

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

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 &#34;Multibody Asterodensity Profiling&#34; (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.

preprint2011arXiv

An Independent Analysis of Kepler-4b through Kepler-8b

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

preprint2011arXiv

Analysis of Kepler&#39;s short-cadence photometry for TrES-2b

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&#39;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.

preprint2011arXiv

Detection of visible light from the darkest world

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&#39; ~ 0.3-0.4) provide a compatible explanation to the data.

preprint2011arXiv

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

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.

preprint2010arXiv

How to Weigh a Star Using a Moon

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&#39;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&#39;s composition.

preprint2010arXiv

Investigations of approximate expressions for the transit duration

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.

preprint2010arXiv

Nightside Pollution of Exoplanet Transit Depths

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

preprint2009arXiv

On the detectability of habitable exomoons with Kepler-class photometry

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

preprint2009arXiv

Transit timing effects due to an exomoon II

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