Source author record

Ruth Angus

Ruth Angus 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

24works
4topics
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

24 published item(s)

preprint2026arXiv

Quantifying the Contamination of TESS Ecliptic-Plane Light Curves by Minor Planets

Though missions devoted to time series photometry focus primarily on targets far beyond the solar system, their observations can be contaminated by foreground minor planets, especially near the ecliptic plane where solar system objects are most prevalent. Crucially, depending on one's choice of data reduction/background estimation algorithm, these objects can induce both apparent brightening and/or dimming events in processed light curves. To quantify the impact of these objects on archived TESS light curves, we used N-body integrations of all currently known minor planets to postdict all 600,000+ of their interactions with stars selected for high-cadence observations during TESS ecliptic plane sectors. We then created mock images of these moving sources and performed simple aperture photometry using the same target and background apertures used in SPOC processing. Our resulting 10,000+ target-specific light curves, which faithfully model the time-dependent positions and magnitudes of the actual solar system objects that approached each target, reveal that $>95\%$ of high-cadence ecliptic plane targets experience a minor planet crossing within 1 TESS pixel of the source. Additionally, 50% of all $T>13$ mag targets experience at least one instantaneous moment where the contaminating flux from minor planets exceeds 1% of the target flux. We discuss these population-level results and others, and highlight several case studies of bright flybys.

preprint2022arXiv

Further Evidence of Modified Spin-down in Sun-like Stars: Pileups in the Temperature-Period Distribution

We combine stellar surface rotation periods determined from NASA's Kepler mission with spectroscopic temperatures to demonstrate the existence of pileups at the long-period and short-period edges of the temperature-period distribution for main-sequence stars with temperatures exceeding $\sim 5500$K. The long-period pileup is well-described by a curve of constant Rossby number, with a critical value of $\mathrm{Ro_{crit}} \lesssim 2$. The long-period pileup was predicted by van Saders et al. (2019) as a consequence of weakened magnetic braking, in which wind-driven angular momentum losses cease once stars reach a critical Rossby number. Stars in the long-period pileup are found to have a wide range of ages ($\sim 2-6$Gyr), meaning that, along the pileup, rotation period is strongly predictive of a star's surface temperature but weakly predictive of its age. The short-period pileup, which is also well-described by a curve of constant Rossby number, is not a prediction of the weakened magnetic braking hypothesis but may instead be related to a phase of slowed surface spin-down due to core-envelope coupling. The same mechanism was proposed by Curtis et al. (2020) to explain the overlapping rotation sequences of low-mass members of differently aged open clusters. The relative dearth of stars with intermediate rotation periods between the short- and long-period pileups is also well-described by a curve of constant Rossby number, which aligns with the period gap initially discovered by McQuillan et al. (2013a) in M-type stars. These observations provide further support for the hypothesis that the period gap is due to stellar astrophysics, rather than a non-uniform star-formation history in the Kepler field.

preprint2022arXiv

TESS Giants Transiting Giants II: The hottest Jupiters orbiting evolved stars

Giant planets on short-period orbits are predicted to be inflated and eventually engulfed by their host stars. However, the detailed timescales and stages of these processes are not well known. Here we present the discovery of three hot Jupiters (P $<$ 10 d) orbiting evolved, intermediate-mass stars ($M_\star$ $\approx$ 1.5 M$_\odot$, 2 R$_\odot$ $<$ $R_\star < $ 5 R$_\odot$). By combining \tess photometry with ground-based photometry and radial velocity measurements, we report masses and radii for these three planets between 0.4 and 1.8 M$_\mathrm{J}$ and 0.8 and 1.8 R$_\mathrm{J}$. \planet has the shortest period (P=\period) of any planet discovered around a red giant star to date. Both \planettwo and \planetthree appear to be inflated, but \planet does not show any sign of inflation. The large radii and relatively low masses of \planettwo and \planetthree place them among the lowest density hot Jupiters currently known, while \planet is conversely one of the highest. All three planets have orbital eccentricities below 0.2. The large spread in radii for these systems implies that planet inflation has a complex dependence on planet mass, radius, incident flux, and orbital properties. We predict that \planet has the shortest orbital decay timescale of any planet currently known, but do not detect any orbital decay in this system. Transmission spectroscopy of \planettwo would provide a favorable opportunity for the detection of water, carbon dioxide and carbon monoxide features in the atmosphere of a planet orbiting an evolved star, and could yield new information about planet formation and atmospheric evolution.

preprint2022arXiv

TESS Giants Transiting Giants III: An eccentric warm Jupiter supports a period-eccentricity relation for giant planets transiting evolved stars

The fate of planets around rapidly evolving stars is not well understood. Previous studies have suggested that relative to the main sequence population, planets transiting evolved stars ($P$ $<$ 100 d) tend to have more eccentric orbits. Here we present the discovery of TOI-4582 b, a 0.94 $\pm$ 0.12 R$_\mathrm{J}$, 0.53 $\pm$ 0.05 M$_\mathrm{J}$ planet orbiting an intermediate-mass subgiant star every 31.034 days. We find that this planet is also on a significantly eccentric orbit ($e$ = 0.51 $\pm$ 0.05). We then compare the population of planets found transiting evolved (log$g$ $<$ 3.8) stars to the population of planets transiting main sequence stars. We find that the rate at which median orbital eccentricity grows with period is significantly higher for evolved star systems than for otherwise similar main sequence systems, particularly for systems with only one planet detected. In general, we observe that mean planet eccentricity $<e>$ = $a$ + $b$log$_{10}$($P$) for the evolved population with a single transiting planet where $a$ = (-0.18 $\pm$ 0.08) and $b$ = (0.38 $\pm$ 0.06), significantly distinct from the main sequence planetary system population. This trend is seen even after controlling for stellar mass and metallicity. These systems do not appear to represent a steady evolution pathway from eccentric, long-period planetary orbits to circular, short period orbits, as orbital model comparisons suggest inspiral timescales are uncorrelated with orbital separation or eccentricity. Characterization of additional evolved planetary systems will distinguish effects of stellar evolution from those of stellar mass and composition.

preprint2022arXiv

The 3D Galactocentric velocities of Kepler stars: marginalizing over missing RVs

Precise Gaia measurements of positions, parallaxes, and proper motions provide an opportunity to calculate 3D positions and 2D velocities (i.e. 5D phase-space) of Milky Way stars. Where available, spectroscopic radial velocity (RV) measurements provide full 6D phase-space information, however there are now and will remain many stars without RV measurements. Without an RV it is not possible to directly calculate 3D stellar velocities, however one can infer 3D stellar velocities by marginalizing over the missing RV dimension. In this paper, we infer the 3D velocities of stars in the Kepler field in Cartesian Galactocentric coordinates (vx, vy, vz). We directly calculate velocities for around a quarter of all Kepler targets, using RV measurements available from the Gaia, LAMOST and APOGEE spectroscopic surveys. Using the velocity distributions of these stars as our prior, we infer velocities for the remaining three-quarters of the sample by marginalizing over the RV dimension. The median uncertainties on our inferred vx, vy, and vz velocities are around 4, 18, and 4 km/s, respectively. We provide 3D velocities for a total of 148,590 stars in the Kepler field. These 3D velocities could enable kinematic age-dating, Milky Way stellar population studies, and other scientific studies using the benchmark sample of well-studied Kepler stars. Although the methodology used here is broadly applicable to targets across the sky, our prior is specifically constructed from and for the Kepler field. Care should be taken to use a suitable prior when extending this method to other parts of the Galaxy.

preprint2022arXiv

The unpopular Package: a Data-driven Approach to De-trend TESS Full Frame Image Light Curves

The majority of observed pixels on the Transiting Exoplanet Survey Satellite (TESS) are delivered in the form of full frame images (FFI). However, the FFIs contain systematic effects such as pointing jitter and scattered light from the Earth and Moon that must be removed before downstream analysis. We present unpopular, an open-source Python package to de-trend TESS FFI light curves based on the causal pixel model method. Under the assumption that shared flux variations across multiple distant pixels are likely to be systematics, unpopular removes these common (i.e., popular) trends by modeling the systematics in a given pixel's light curve as a linear combination of light curves from many other distant pixels. To prevent overfitting we employ ridge regression and a train-and-test framework where the data points being de-trended are separated from those used to obtain the model coefficients. We also allow for simultaneous fitting with a polynomial model to capture any long-term astrophysical trends. We validate our method by de-trending different sources (e.g., supernova, tidal disruption event, exoplanet-hosting star, fast rotating star) and comparing our light curves to those obtained by other pipelines when appropriate. We also show that unpopular is able to preserve sector-length astrophysical signals, allowing for the extraction of multi-sector light curves from the FFI data. The unpopular source code and tutorials are freely available online.

preprint2022arXiv

wdwarfdate: A Python Package to Derive Bayesian Ages of White Dwarfs

White dwarfs have been successfully used as cosmochronometers in the literature, however their reach has been limited in comparison to their potential. We present wdwarfdate, a publicly available Python package to derive the Bayesian age of a white dwarf, based on its effective temperature (Teff) and surface gravity (logg). We make this software easy to use with the goal of transforming the usage of white dwarfs as cosmochronometers into an accessible tool. The code estimates the mass and cooling age of the white dwarf, as well as the mass and main-sequence age of the progenitor star, allowing for a determination of the total age of the object. We test the reliability of the method by estimating the parameters of white dwarfs from previous studies, and find agreement with the literature within measurement errors. By analyzing the limitation of the code we find a typical uncertainty of 10% on the total age when both input parameters have uncertainties of 1%, and an uncertainty of 25% on the total age when Teff has an uncertainty of 10% and logg of 1%. Furthermore, wdwarfdate assumes single star evolution, and can be applied to calculate the total age of a white dwarf with parameters in the range 1,500<Teff<90,000 K and 7.9<logg<9.3. Finally, the code assumes a uniform mixture of C/O in the core and single star evolution, which is reliable in the range of white dwarf masses 0.45-1.1 Msun (7.73<logg<8.8).

preprint2021arXiv

$\textit{TESS}$ Giants Transiting Giants I: A Non-inflated Hot Jupiter Orbiting a Massive Subgiant

While the population of confirmed exoplanets continues to grow, the sample of confirmed transiting planets around evolved stars is still limited. We present the discovery and confirmation of a hot Jupiter orbiting TOI-2184 (TIC 176956893), a massive evolved subgiant ($M_\star= 1.53 \pm 0.12 M_\odot$, $R_\star= 2.90 \pm 0.14 R_\odot$) in the $\textit{TESS}$ Southern Continuous Viewing Zone. The planet was flagged as a false positive by the $\textit{TESS}$ Quick-Look Pipeline due to periodic systematics introducing a spurious depth difference between even and odd transits. Using a new pipeline to remove background scattered light in $\textit{TESS}$ Full Frame Image (FFI) data, we combine space-based $\textit{TESS}$ photometry, ground-based photometry, and ground-based radial velocity measurements to report a planet radius of $R_p= 1.017 \pm 0.051 R_J$ and mass of $M_p= 0.65 \pm 0.16 M_J$. For a planet so close to its star, the mass and radius of TOI-2184b are unusually well matched to those of Jupiter. We find that the radius of TOI-2184b is smaller than theoretically predicted based on its mass and incident flux, providing a valuable new constraint on the timescale of post-main-sequence planet inflation. The discovery of TOI-2184b demonstrates the feasibility of detecting planets around faint ($\textit{TESS}$ magnitude $>12$) post-main sequence stars and suggests that many more similar systems are waiting to be detected in the $\textit{TESS}$ FFIs, whose confirmation may elucidate the final stages of planetary system evolution.

preprint2021arXiv

A decade of radial-velocity monitoring of Vega and new limits on the presence of planets

We present an analysis of 1524 spectra of Vega spanning 10 years, in which we search for periodic radial velocity variations. A signal with a periodicity of 0.676 days and a semi-amplitude of ~10 m/s is consistent with the rotation period measured over much shorter time spans by previous spectroscopic and spectropolarimetric studies, confirming the presence of surface features on this A0 star. The timescale of evolution of these features can provide insight into the mechanism that sustains the weak magnetic fields in normal A type stars. Modeling the radial velocities with a Gaussian process using a quasi-periodic kernel suggests that the characteristic spot evolution timescale is ~180 days, though we cannot exclude the possibility that it is much longer. Such long timescales may indicate the presence of failed fossil magnetic fields on Vega. TESS data reveal Vega's photometric rotational modulation for the first time, with a total amplitude of only 10 ppm, and a comparison of the spectroscopic and photometric amplitudes suggest the surface features may be dominated by bright plages rather than dark spots. For the shortest orbital periods, transit and radial velocity injection recovery tests exclude the presence of transiting planets larger than 2 Earth radii and most non-transiting giant planets. At long periods, we combine our radial velocities with direct imaging from the literature to produce detection limits for Vegan planets and brown dwarfs out to distances of 15 au. Finally, we detect a candidate radial velocity signal with a period of 2.43 days and a semi-amplitude of 6 m/s. If caused by an orbiting companion, its minimum mass would be ~20 Earth masses; because of Vega's pole-on orientation, this would correspond to a Jovian planet if the orbit is aligned with the stellar spin. We discuss the prospects for confirmation of this candidate planet.

preprint2021arXiv

Gyro-Kinematic Ages for around 30,000 Kepler Stars

Estimating stellar ages is important for advancing our understanding of stellar and exoplanet evolution and investigating the history of the Milky Way. However, ages for low-mass stars are hard to infer as they evolve slowly on the main sequence. In addition, empirical dating methods are difficult to calibrate for low-mass stars as they are faint. In this work, we calculate ages for Kepler F, G, and crucially K and M dwarfs, using their rotation and kinematic properties. We apply the simple assumption that the velocity dispersion of stars increases over time and adopt an age--velocity--dispersion relation (AVR) to estimate average stellar ages for groupings of coeval stars. We calculate the vertical velocity dispersion of stars in bins of absolute magnitude, temperature, rotation period, and Rossby number and then convert velocity dispersion to kinematic age via an AVR. Using this method, we estimate gyro-kinematic ages for 29,949 Kepler stars with measured rotation periods. We are able to estimate ages for clusters and asteroseismic stars with an RMS of 1.22 Gyr and 0.26 Gyr respectively. With our Astraea machine learning algorithm, which predicts rotation periods, we suggest a new selection criterion (a weight of 0.15) to increase the size of the McQuillian et al. (2014) catalog of Kepler rotation periods by up to 25%. Using predicted rotation periods, we estimated gyro-kinematic ages for stars without measured rotation periods and found promising results by comparing 12 detailed age--element abundance trends with literature values.

preprint2020arXiv

Astraea: Predicting Long Rotation Periods with 27-Day Light Curves

The rotation periods of planet-hosting stars can be used for modeling and mitigating the impact of magnetic activity in radial velocity measurements and can help constrain the high-energy flux environment and space weather of planetary systems. Millions of stars and thousands of planet hosts are observed with the Transiting Exoplanet Survey Satellite (TESS). However, most will only be observed for 27 contiguous days in a year, making it difficult to measure rotation periods with traditional methods. This is especially problematic for field M dwarfs, which are ideal candidates for exoplanet searches, but which tend to have periods in excess of the 27-day observing baseline. We present a new tool, Astraea, for predicting long rotation periods from short-duration light curves combined with stellar parameters from Gaia DR2. Using Astraea, we can predict the rotation periods from Kepler 4-year light curves with 13% uncertainty overall (and a 9% uncertainty for periods > 30 days). By training on 27-day Kepler light curve segments, Astraea can predict rotation periods up to 150 days with 9% uncertainty (5% for periods > 30 days). After training this tool on these 27-day Kepler light curve segments, we applied \texttt{Astraea} to real TESS data. For the 195 stars that were observed by both Kepler and TESS, we were able to predict the rotation periods with 55% uncertainty despite the wild differences in systematics.

preprint2020arXiv

Data-driven derivation of stellar properties from photometric time series data using convolutional neural networks

Stellar variability is driven by a multitude of internal physical processes that depend on fundamental stellar properties. These properties are our bridge to reconciling stellar observations with stellar physics, and for understanding the distribution of stellar populations within the context of galaxy formation. Numerous ongoing and upcoming missions are charting brightness fluctuations of stars over time, which encode information about physical processes such as rotation period, evolutionary state (such as effective temperature and surface gravity), and mass (via asteroseismic parameters). Here, we explore how well we can predict these stellar properties, across different evolutionary states, using only photometric time series data. To do this, we implement a convolutional neural network, and with data-driven modeling we predict stellar properties from light curves of various baselines and cadences. Based on a single quarter of \textit{Kepler} data, we recover stellar properties, including surface gravity for red giant stars (with an uncertainty of $\lesssim$ 0.06 dex), and rotation period for main sequence stars (with an uncertainty of $\lesssim$ 5.2 days, and unbiased from $\approx$5 to 40 days). Shortening the \textit{Kepler} data to a 27-day TESS-like baseline, we recover stellar properties with a small decrease in precision, $\sim$0.07 dex for log $g$ and $\sim$5.5 days for $P_{\rm rot}$, unbiased from $\approx$5 to 35 days. Our flexible data-driven approach leverages the full information content of the data, requires minimal feature engineering, and can be generalized to other surveys and datasets. This has the potential to provide stellar property estimates for many millions of stars in current and future surveys.

preprint2020arXiv

Exploring the evolution of stellar rotation using Galactic kinematics

The rotational evolution of cool dwarfs is poorly constrained after around 1-2 Gyr due to a lack of precise ages and rotation periods for old main-sequence stars. In this work we use velocity dispersion as an age proxy to reveal the temperature-dependent rotational evolution of low-mass Kepler dwarfs, and demonstrate that kinematic ages could be a useful tool for calibrating gyrochronology in the future. We find that a linear gyrochronology model, calibrated to fit the period-Teff relationship of the Praesepe cluster, does not apply to stars older than around 1 Gyr. Although late-K dwarfs spin more slowly than early-K dwarfs when they are young, at old ages we find that late-K dwarfs rotate at the same rate or faster than early-K dwarfs of the same age. This result agrees qualitatively with semi-empirical models that vary the rate of surface-to-core angular momentum transport as a function of time and mass. It also aligns with recent observations of stars in the NGC 6811 cluster, which indicate that the surface rotation rates of K dwarfs go through an epoch of inhibited evolution. We find that the oldest Kepler stars with measured rotation periods are late-K and early-M dwarfs, indicating that these stars maintain spotted surfaces and stay magnetically active longer than more massive stars. Finally, based on their kinematics, we confirm that many rapidly rotating GKM dwarfs are likely to be synchronized binaries.

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.

preprint2020arXiv

WASP-107b's density is even lower: a case study for the physics of planetary gas envelope accretion and orbital migration

With a mass in the Neptune regime and a radius of Jupiter, WASP-107b presents a challenge to planet formation theories. Meanwhile, the planet's low surface gravity and the star's brightness also make it one of the most favorable targets for atmospheric characterization. Here, we present the results of an extensive 4-year Keck/HIRES radial-velocity (RV) follow-up program of the WASP-107 system and provide a detailed study of the physics governing the accretion of its gas envelope. We reveal that WASP-107b's mass is only 1.8 Neptune masses ($M_b = 30.5 \pm 1.7$ $M_\oplus$). The resulting extraordinarily low density suggests that WASP-107b has a H/He envelope mass fraction of $> 85$% unless it is substantially inflated. The corresponding core mass of $<4.6$ $M_\oplus$ at 3$σ$ is significantly lower than what is traditionally assumed to be necessary to trigger massive gas envelope accretion. We demonstrate that this large gas-to-core mass ratio most plausibly results from the onset of accretion at $\gtrsim 1$ AU onto a low-opacity, dust-free atmosphere and subsequent migration to the present-day $a_b = 0.0566 \pm 0.0017$ AU. Beyond WASP-107b, we also detect a second more massive planet ($M_c \sin i = 0.36 \pm 0.04$ $M_{J}$) on a wide eccentric orbit ($e_c = 0.28 \pm 0.07$) which may have influenced the orbital migration and spin-orbit misalignment of WASP-107b. Overall, our new RV observations and envelope accretion modeling provide crucial insights into the intriguing nature of WASP-107b and the system's formation history. Looking ahead, WASP-107b will be a keystone planet to understand the physics of gas envelope accretion.

preprint2016arXiv

Maximizing Science in the Era of LSST, Stars Study Group Report: Rotation and Magnetic Activity in the Galactic Field Population and in Open Star Clusters

This is the stars chapter of the Kavli workshop report, which resulted from the community-based study of needed US OIR capabilities in the LSST era. The full report, which will include this chapter, is anticipated to be available in Fall 2016. See NOAO website (http://www.noao.edu/meetings/lsst-oir-study/) for more details.

preprint2016arXiv

Maximizing Science in the Era of LSST: A Community-Based Study of Needed US Capabilities

The Large Synoptic Survey Telescope (LSST) will be a discovery machine for the astronomy and physics communities, revealing astrophysical phenomena from the Solar System to the outer reaches of the observable Universe. While many discoveries will be made using LSST data alone, taking full scientific advantage of LSST will require ground-based optical-infrared (OIR) supporting capabilities, e.g., observing time on telescopes, instrumentation, computing resources, and other infrastructure. This community-based study identifies, from a science-driven perspective, capabilities that are needed to maximize LSST science. Expanding on the initial steps taken in the 2015 OIR System Report, the study takes a detailed, quantitative look at the capabilities needed to accomplish six representative LSST-enabled science programs that connect closely with scientific priorities from the 2010 decadal surveys. The study prioritizes the resources needed to accomplish the science programs and highlights ways that existing, planned, and future resources could be positioned to accomplish the science goals.

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.

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

Calibrating Gyrochronology using Kepler Asteroseismic targets

Among the available methods for dating stars, gyrochronology is a powerful one because it requires knowledge of only the star's mass and rotation period. Gyrochronology relations have previously been calibrated using young clusters, with the Sun providing the only age dependence, and are therefore poorly calibrated at late ages. We used rotation period measurements of 310 Kepler stars with asteroseismic ages, 50 stars from the Hyades and Coma Berenices clusters and 6 field stars (including the Sun) with precise age measurements to calibrate the gyrochronology relation, whilst fully accounting for measurement uncertainties in all observable quantities. We calibrated a relation of the form $P=A^n\times(B-V-c)^b$, where $P$ is rotation period in days, $A$ is age in Myr, $B$ and $V$ are magnitudes and $a$, $b$ and $n$ are the free parameters of our model. We found $a = 0.40^{+0.3}_{-0.05}$, $b = 0.31^{+0.05}_{-0.02}$ and $n = 0.55^{+0.02}_{-0.09}$. Markov Chain Monte Carlo methods were used to explore the posterior probability distribution functions of the gyrochronology parameters and we carefully checked the effects of leaving out parts of our sample, leading us to find that no single relation beween rotation period, colour and age can adequately describe all the subsets of our data. The Kepler asteroseismic stars, cluster stars and local field stars cannot all be described by the same gyrochronology relation. The Kepler asteroseismic stars may be subject to observational biases, however the clusters show unexpected deviations from the predicted behaviour, providing concerns for the overall reliability of gyrochronology as a dating method.

preprint2015arXiv

Systematics-insensitive periodic signal search with K2

From pulsating stars to transiting exoplanets, the search for periodic signals in K2 data, Kepler's 2-wheeled extension, is relevant to a long list of scientific goals. Systematics affecting K2 light curves due to the decreased spacecraft pointing precision inhibit the easy extraction of periodic signals from the data. We here develop a method for producing periodograms of K2 light curves that are insensitive to pointing-induced systematics; the Systematics-Insensitive Periodogram (SIP). Traditional sine-fitting periodograms use a generative model to find the frequency of a sinusoid that best describes the data. We extend this principle by including systematic trends, based on a set of 'Eigen light curves', following Foreman-Mackey et al. (2015), in our generative model as well as a sum of sine and cosine functions over a grid of frequencies. Using this method we are able to produce periodograms with vastly reduced systematic features. The quality of the resulting periodograms are such that we can recover acoustic oscillations in giant stars and measure stellar rotation periods without the need for any detrending. The algorithm is also applicable to the detection of other periodic phenomena such as variable stars, eclipsing binaries and short-period exoplanet candidates. The SIP code is available at https://github.com/RuthAngus/SIPK2.

preprint2014arXiv

Characterizing K2 Planet Discoveries: A super-Earth transiting the bright K-dwarf HIP 116454

We report the first planet discovery from the two-wheeled Kepler (K2) mission: HIP 116454 b. The host star HIP 116454 is a bright (V = 10.1, K = 8.0) K1-dwarf with high proper motion, and a parallax-based distance of 55.2 +/- 5.4 pc. Based on high-resolution optical spectroscopy, we find that the host star is metal-poor with [Fe/H] = -.16 +/- .18, and has a radius R = 0.716 +/- .0024 R_sun and mass M = .775 +/- .027 Msun. The star was observed by the Kepler spacecraft during its Two-Wheeled Concept Engineering Test in February 2014. During the 9 days of observations, K2 observed a single transit event. Using a new K2 photometric analysis technique we are able to correct small telescope drifts and recover the observed transit at high confidence, corresponding to a planetary radius of Rp = 2.53 +/- 0.18 Rearth. Radial velocity observations with the HARPS-N spectrograph reveal a 11.82 +/- 1.33 Mearth planet in a 9.1 day orbit, consistent with the transit depth, duration, and ephemeris. Follow-up photometric measurements from the MOST satellite confirm the transit observed in the K2 photometry and provide a refined ephemeris, making HIP 116454 b amenable for future follow-up observations of this latest addition to the growing population of transiting super-Earths around nearby, bright stars.

preprint2013arXiv

Maximizing Kepler science return per telemetered pixel: Detailed models of the focal plane in the two-wheel era

Kepler's immense photometric precision to date was maintained through satellite stability and precise pointing. In this white paper, we argue that image modeling--fitting the Kepler-downlinked raw pixel data--can vastly improve the precision of Kepler in pointing-degraded two-wheel mode. We argue that a non-trivial modeling effort may permit continuance of photometry at 10-ppm-level precision. We demonstrate some baby steps towards precise models in both data-driven (flexible) and physics-driven (interpretably parameterized) modes. We demonstrate that the expected drift or jitter in positions in the two-weel era will help with constraining calibration parameters. In particular, we show that we can infer the device flat-field at higher than pixel resolution; that is, we can infer pixel-to-pixel variations in intra-pixel sensitivity. These results are relevant to almost any scientific goal for the repurposed mission; image modeling ought to be a part of any two-wheel repurpose for the satellite. We make other recommendations for Kepler operations, but fundamentally advocate that the project stick with its core mission of finding and characterizing Earth analogs. [abridged]

preprint2013arXiv

Maximizing Kepler science return per telemetered pixel: Searching the habitable zones of the brightest stars

In today's mailing, Hogg et al. propose image modeling techniques to maintain 10-ppm-level precision photometry in Kepler data with only two working reaction wheels. While these results are relevant to many scientific goals for the repurposed mission, all modeling efforts so far have used a toy model of the Kepler telescope. Because the two-wheel performance of Kepler remains to be determined, we advocate for the consideration of an alternate strategy for a >1 year program that maximizes the science return from the "low-torque" fields across the ecliptic plane. Assuming we can reach the precision of the original Kepler mission, we expect to detect 800 new planet candidates in the first year of such a mission. Our proposed strategy has benefits for transit timing variation and transit duration variation studies, especially when considered in concert with the future TESS mission. We also expect to help address the first key science goal of Kepler: the frequency of planets in the habitable zone as a function of spectral type.