Source author record

David Kipping

David Kipping appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

26 published item(s)

preprint2026arXiv

A JWST Transit of a Jupiter Analog I: Constraints on the Oblateness of Kepler-167e

In October 2024 JWST observed a transit of Kepler-167e, a Jupiter-analog planet on a 1000+ day orbit. These observations, recorded over a long baseline of nearly 60 hours, were designed to search for signatures of planetary oblateness and/or exomoons comparable to Ganymede. In this first in a series of studies analyzing these data we report on constraints on Kepler-167e's oblateness. We explored a large grid of data reduction pipelines and modeling choices, including a new entirely independent reduction pipeline ("katahdin") and two new treatments for limb darkening. We find that under a Bayesian model comparison framework the data are fit equally well by both spherical and oblate planet models, and that our ability to constrain the oblateness is negatively impacted by the influence of exposure-long trends. Using the most conservative of our posteriors, we place a 95% upper bound on the projected oblateness of $f<0.097$, which corresponds to a rotation period of $P\geq7.11$ hours if the planet's spin axis is aligned with the sky plane. We note, however, that the final bound depends on the choice of reduction pipeline and systematics model, and that our suite of end-to-end analyses produced bounds as low as $f<0.065$ at 95%. We conclude that leveraging JWST to make tighter constraints on planetary oblateness will require further investigation into mitigating exposure-long trends and correlated noise.

preprint2024arXiv

Giant Outer Transiting Exoplanet Mass (GOT 'EM) Survey. IV. Long-term Doppler Spectroscopy for 11 Stars Thought to Host Cool Giant Exoplanets

Discovering and characterizing exoplanets at the outer edge of the transit method's sensitivity has proven challenging owing to geometric biases and the practical difficulties associated with acquiring long observational baselines. Nonetheless, a sample of giant exoplanets on orbits longer than 100 days has been identified by transit hunting missions. We present long-term Doppler spectroscopy for 11 such systems with observation baselines spanning a few years to a decade. We model these radial velocity observations jointly with transit photometry to provide initial characterizations of these objects and the systems in which they exist. Specifically, we make new precise mass measurements for four long-period giant exoplanets (Kepler-111 c, Kepler-553 c, Kepler-849 b, and PH-2 b), we place new upper limits on mass for four others (Kepler-421 b, KOI-1431.01, Kepler-1513 b, and Kepler-952 b), and we show that several "confirmed" planets are in fact not planetary at all. We present these findings to complement similar efforts focused on closer-in short-period giant planets, and with the hope of inspiring future dedicated studies of cool giant exoplanets.

preprint2022arXiv

An Exomoon Survey of 70 Cool Giant Exoplanets and the New Candidate Kepler-1708 b-i

Exomoons represent a crucial missing puzzle piece in our efforts to understand extrasolar planetary systems. To address this deficiency, we here describe an exomoon survey of 70 cool, giant transiting exoplanet candidates found by Kepler. We identify only one which exhibits a moon-like signal that passes a battery of vetting tests: Kepler-1708 b. We show that Kepler-1708 b is a statistically validated Jupiter-sized planet orbiting a Sun-like quiescent star at ~1.6AU. The signal of the exomoon candidate, Kepler-1708 b-i, is a 4.8-sigma effect and is persistent across different instrumental detrending methods, with a 1% false-positive probability via injection-recovery. Kepler-1708 b-i is ~2.6 Earth radii and is located in an approximately coplanar orbit at ~12 planetary radii from its ~1.6AU Jupiter-sized host. Future observations will be necessary to validate or reject the candidate.

preprint2022arXiv

Could the "Wow" signal have originated from a stochastic repeating beacon?

The famous "Wow" signal detected in 1977 remains arguably the most compelling SETI signal ever found. The original Big Ear data requires that the signal turned on/off over the span of ~3 minutes (time difference between the dual antennae), yet persisted for 72 seconds (duration of a single beam sweep). Combined with the substantial and negative follow-up efforts, these observations limit the allowed range of signal repeat schedules, to the extent that one might question the credibility of the signal itself. Previous work has largely excluded the hypothesis of a strictly periodic repeating source, for periods shorter than 40 hours. However, a non-periodic, stochastic repeater remains largely unexplored. Here, we employ a likelihood emulator using the Big Ear observing logs to infer the probable signal properties under this hypothesis. We find that the maximum a-posteriori solution has a likelihood of 32.3%, highly compatible with the Big Ear data, with a broad 2 $σ$ credible interval of signal duration 72 secs < T < 77 mins and mean repeat rate 0.043 1/days < $λ$ < 59.8 1/days. We extend our analysis to include 192 hours of subsequent observations from META, Hobart and ATA, which drops the peak likelihood to 1.78%, and thus in tension with the available data at the 2.4 $σ$ level. Accordingly, the Wow signal cannot be excluded as a stochastic repeater with available data, and we estimate that 62 days of accumulated additional observations would be necessary to surpass 3 $σ$ confidence.

preprint2022arXiv

Kepler-1708 b-i is likely undetectable with HST

The exomoon candidate Kepler-1708 b-i was recently reported using two transits of Kepler data. Supported by a 1% false-positive probability, the candidate is promising but requires follow-up observations to confirm/reject its validity. In this paper, we consider the detectability of the exomoon candidate's transit, most specifically in the next window (March 2023) using the WFC3 instrument aboard the Hubble Space Telescope (HST). Using realistic noise estimates, accounting for the visit-long trends, and propagating the model posteriors derived using the Kepler data, we perform 75 injection-recovery trials with Bayesian model selection. Defining a successful detection as one which meets thresholds of the Bayes factor, AIC, and error of the retrieved parameters, only 7 of our 75 injections were recovered when considering HST data alone. This implies a true-positive probability (TPP) of $10\pm3$%. Despite HST's superior aperture to Kepler, both instrumental systematics and the compactness of the candidate exomoon's orbit typically obfuscate a strong detection. Although the noise properties of the James Webb Space Telescope (JWST) have not yet been characterized in flight, we estimate the signal would be easily recovered using NIRSpec operating in its Bright Object Time Series mode.

preprint2022arXiv

Mathematical encoding within multi-resonant planetary systems as SETI beacons

How might an advanced alien civilization manipulate the orbits within a planetary system to create a durable signpost that communicates its existence? While it is still debated whether such a purposeful advertisement would be prudent and wise, we propose that mean-motion resonances between neighboring planets -- with orbital periods that form integer ratios -- could in principle be used to encode simple sequences that one would not expect to form in nature. In this Letter we build four multi-resonant planetary systems and test their long-term orbital stability. The four systems each contain 6 or 7 planets and consist of: (i) consecutive integers from 1 to 6; (ii) prime numbers from 2 to 11; (iii) the Fibonacci sequence from 1 to 13; and (iv) the Lazy Caterer sequence from 1 to 16. We built each system using N-body simulations with artificial migration forces. We evaluated the stability of each system over the full 10 Gyr integration of the Sun's main sequence phase. We then tested the stability of these systems for an additional 10 Gyr, during and after post-main sequence evolution of the central stars (assumed to be Sun-like) to their final, white dwarf phase. The only system that was destabilized was the consecutive integer sequence (system i). The other three sequences therefore represent potential SETI beacons.

preprint2022arXiv

The Number of Transits Per Epoch for Transiting Misaligned Circumbinary Planets

The growing catalog of circumbinary planets strengthens the notion that planets form in a diverse range of conditions across the cosmos. Transiting circumbinary planets yield especially important insights and many examples are now known, in broadly coplanar obits with respect to their binary. Studies of circumbinary disks suggest misaligned transiting examples could also plausibly exist, but their existence would exacerbate the already challenging feat of automatic detection. In this work, we synthesize populations of such planets and consider the number of transits per epoch they produce, forming integer sequences. For isotropic distributions, such sequences will appear foreign to conventional expectation, rarely (~1%) producing the signature double-transits we have come to expect for circumbinaries, instead producing sparse sequences dominated by zero-transit epochs (~80%). Despite their strangeness, we demonstrate that these sequences will be non-random and that the two preceding epochs predict the next to high accuracy. Additionally, we show that even when clustering the transits into grouped epochs, they often appear unphysical if erroneously assuming a single star, due to the missing epochs. Crucially, missing epochs mean highly isotropic populations can trick the observer into assigning the wrong period in up to a quarter of cases, adding further confusion. Finally, we show that the transit sequences encode the inclination distribution and demonstrate a simple inference method that successfully matches the injected truth. Our work highlights how the simple act of flagging transits can be used to provide an initial, vetting-level analysis of misaligned transiting circumbinary planets.

preprint2021arXiv

Contact Inequality -- First Contact Will Likely Be With An Older Civilization

First contact with another civilization, or simply another intelligence of some kind, will likely be quite different depending on whether that intelligence is more or less advanced than ourselves. If we assume that the lifetime distribution of intelligences follows an approximately exponential distribution, one might naively assume that the pile-up of short-lived entities dominates any detection or contact scenario. However, it is argued here that the probability of contact is proportional to the age of said intelligence (or possibly stronger), which introduces a selection effect. We demonstrate that detected intelligences will have a mean age twice that of the underlying (detected + undetected) population, using the exponential model. We find that our first contact will most likely be with an older intelligence, provided that the maximum allowed mean lifetime of the intelligence population, $τ_{\rm max}$, is $>e$ times larger than our own. Older intelligences may be rare but they disproportionality contribute to first contacts, introducing what we call a 'contact inequality', analogous to wealth inequality. This reasoning formalizes intuitional arguments and highlights that first contact would likely be one-sided, with ramifications for how we approach SETI.

preprint2021arXiv

Validation of HD 183579b using archival radial velocities: a warm-neptune orbiting a bright solar analog

As exoplanetary science matures into its third decade, we are increasingly offered the possibility of pre existing, archival observations for newly detected candidates. This is particularly poignant for the TESS mission, whose survey spans bright, nearby dwarf stars in both hemispheres, which are precisely the types of sources targeted by previous radial velocity (RV) surveys. On this basis, we investigated whether any of the TESS Objects of Interest (TOIs) coincided with such observations, from which we find 18 single planet candidate systems. Of these, one exhibits an RV signature that has the correct period and phase matching the transiting planetary candidate with a false alarm probability of less than 1 percent. After further checks, we exploit this fact to validate HD 183579b (TOI-1055b). This planet is less than 4 Earth Radii and has better than 33 percent planetary mass measurements, thus advancing the TESS primary objective of finding 50 such worlds. We find that this planet is amongst the most accessible small transiting planets for atmospheric characterization. Our work highlights that the efforts to confirm and even precisely measure the masses of new transiting planet candidates need not always depend on acquiring new observations - that in some instances these tasks can be completed with existing data.

preprint2020arXiv

8 in 10 Stars in the Milky Way Bulge Experience Stellar Encounters Within 1000 AU in a Gigayear

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

preprint2020arXiv

A Bayesian Approach to the Simulation Argument

The Simulation Argument posed by Bostrom (2003) suggests that we may be living inside a sophisticated computer simulation. If post-human civilizations eventually have both the capability and desire to generate such Bostrom-like simulations, then the number of simulated realities would greatly exceed the one base reality, ostensibly indicating a high probability that we do not live in said base reality. In this work, it is argued that since the hypothesis that such simulations are technically possible remains unproven, then statistical calculations need to consider not just the number of state spaces, but the intrinsic model uncertainty. This is achievable through a Bayesian treatment of the problem, which is presented here. Using Bayesian model averaging, it is shown that the probability that we are sims is in fact less than 50%, tending towards that value in the limit of an infinite number of simulations. This result is broadly indifferent as to whether one conditions upon the fact that humanity has not yet birthed such simulations, or ignore it. As argued elsewhere, it is found that if humanity does start producing such simulations, then this would radically shift the odds and make it very probable that we are in fact sims.

preprint2020arXiv

A planet within the debris disk around the pre-main-sequence star AU Microscopii

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

preprint2020arXiv

An Independent Analysis of the Six Recently Claimed Exomoon Candidates

It has been recently claimed that KOIs-268.01, 303.01, 1888.01, 1925.01, 2728.01 & 3320.01 are exomoon candidates, based on an analysis of their transit timing. Here, we perform an independent investigation, which is framed in terms of three questions: 1) Are there significant excess TTVs? 2) Is there a significant periodic TTV? 3) Is there evidence for a non-zero moon mass? We applied rigorous statistical methods to these questions alongside a re-analysis of the Kepler photometry and find that none of the KOIs satisfy these three tests. Specifically, KOIs-268.01 & 3220.01 pass none of the tests and KOIs-303.01, 1888.01 & 1925.01 pass a single test each. Only KOI-2728.01 satisfies two, but fails the cross-validation test for predictions. Further, detailed photodynamical modeling reveals that KOI-2728.01 favours a negative radius moon (as does KOI-268.01). We also note that we find a significant photoeccentric for KOI-1925.01 indicating an eccentric orbit of e>(0.62+/-0.06). For comparison, we applied the same tests to Kepler-1625b, which reveals that 1) and 3) are passed, but 2) cannot be checked with the cross-validation method used here, due to the limited number of available epochs. In conclusion, we find no compelling evidence for exomoons amongst the six KOIs. Despite this, we're able to derive exomoon mass upper limits versus semi-major axis, with KOI-3220.01 leading to particularly impressive constraints of Ms/Mp < 0.4% [2 sigma] at a similar relative semi-major to that of the Earth-Moon.

preprint2020arXiv

An Objective Bayesian Analysis of Life's Early Start and Our Late Arrival

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

preprint2020arXiv

Detection of the phase curve and occultation of WASP-100b with TESS

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

preprint2020arXiv

Impossible moons -- Transit timing effects that cannot be due to an exomoon

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

preprint2020arXiv

Kepler's Dark Worlds: a Low Albedo for an Ensemble of Neptunian and Terran Exoplanets

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

preprint2020arXiv

Loose Ends for the Exomoon Candidate Host Kepler-1625b

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

preprint2016arXiv

Probabilistic inference of basic stellar parameters: application to flickering stars

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

preprint2016arXiv

Revised Masses and Densities of the Planets around Kepler-10

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

preprint2015arXiv

A disintegrating minor planet transiting a white dwarf

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

preprint2015arXiv

A novel method for identifying exoplanetary rings

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

preprint2014arXiv

Photo-Dynamical Analysis of Three Kepler Objects of Interest with Significant Transit Timing Variations

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

preprint2013arXiv

Dynamical Effects on the Habitable Zone for Earth-like Exomoons

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

preprint2013arXiv

KOI-142, the King of Transit Variations, is a Pair of Planets near the 2:1 Resonance

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

preprint2012arXiv

CfA4: Light Curves for 94 Type Ia Supernovae

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