Engineering Papers⌕ Search

Engineering topics

Eric Agol

Publications and source records attributed to Eric Agol.

15 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

TIC 172900988: A Transiting Circumbinary Planet Detected in One Sector of TESS Data

We report the first discovery of a transiting circumbinary planet detected from a single sector of Transiting Exoplanet Survey Satellite (TESS) data. During Sector 21, the planet TIC 172900988b transited the primary star and then five days later it transited the secondary star. The binary is itself eclipsing, with a period P ≈ 19.7 days and an eccentricity e ≈ 0.45. Archival data from ASAS-SN, Evryscope, KELT, and SuperWASP reveal a prominent apsidal motion of the binary orbit, caused by the dynamical interactions between the binary and the planet. A comprehensive photodynamical analysis of the TESS, archival and follow-up data yields stellar masses and radii of M(1) = 1.2384 ± 0.0007 Mꙩ and R(1) = 1.3827 ± 0.0016 Rꙩ for the primary and M(2) = 1.2019 ± 0.0007 Mꙩ and R(2) = 1.3124 ± 0.0012 Rꙩ for the secondary. The radius of the planet is R(3) = 11.25 ± 0.44 Rꚛ (1.004 ± 0.039R(Jup)). The planet’s mass and orbital properties are not uniquely determined—there are six solutions with nearly equal likelihood. Specifically, we find that the planet’s mass is in the range of 824 ≲ M3 ≲ 981 Mꚛ (2.65 ≲ M3 ≲ 3.09M(Jup)), its orbital period could be 188.8, 190.4, 194.0, 199.0, 200.4, or 204.1 days, and the eccentricity is between 0.02 and 0.09. At V = 10.141 mag, the system is accessible for high-resolution spectroscopic observations, e.g., the Rossiter–McLaughlin effect and transit spectroscopy.

Veselin B. Kostov↗
NASA NTRS2021

K2-138 g: Spitzer Spots a Sixth Planet for the Citizen Science System

K2 greatly extended Kepler's ability to find new planets, but it was typically limited to identifying transiting planets with orbital periods below 40 days. While analyzing K2 data through the Exoplanet Explorers project, citizen scientists helped discover one super-Earth and four sub-Neptune sized planets in the relatively bright (V = 12.21, K = 10.3) K2-138 system, all which orbit near 3:2 mean-motion resonances. The K2 light curve showed two additional transit events consistent with a sixth planet. Using Spitzer photometry, we validate the sixth planet's orbital period of 41.966 ± 0.006 days and measure a radius of 3.44 (+0.32,-.031)Rꚛ, solidifying K2-138 as the K2 system with the most currently known planets. There is a sizeable gap between the outer two planets, since the fifth planet in the system, K2-138 f, orbits at 12.76 days. We explore the possibility of additional nontransiting planets in the gap between f and g. Due to the relative brightness of the K2-138 host star, and the near resonance of the inner planets, K2-138 could be a key benchmark system for both radial velocity and transit-timing variation mass measurements, and indeed radial velocity masses for the inner four planets have already been obtained. With its five sub-Neptunes and one super-Earth, the K2-138 system provides a unique test bed for comparative atmospheric studies of warm to temperate planets of similar size, dynamical studies of near-resonant planets, and models of planet formation and migration.

Kevin K. Hardegree-Ullman↗
NASA NTRS2021

Refining the Transit Timing and Photometric Analysis of TRAPPIST-1: Masses, Radii, Densities, Dynamics, and Ephemerides

We have collected transit times for the TRAPPIST-1 system with the Spitzer Space Telescope over four years. We add to these ground-based, HST, and K2 transit-time measurements, and revisit an N-body dynamical analysis of the seven-planet system using our complete set of times from which we refine the mass ratios of the planets to the star. We next carry out a photodynamical analysis of the Spitzer light curves to derive the density of the host starand the planet densities. We find that all seven planets’ densities may be described with a single rocky mass–radius relation which is depleted in iron relative to Earth, with Fe 21 wt% versus 32 wt% for Earth, and otherwise Earth-like in composition. Alternatively, the planets may have an Earth-like composition but enhanced in light elements,such as a surface water layer or a core-free structure with oxidized iron in the mantle. We measure planet masses toa precision of 3%–5%, equivalent to a radial-velocity(RV)precision of 2.5 cm s−1, or two orders of magnitude more precise than current RV capabilities. We find the eccentricities of the planets are very small, the orbits are extremely coplanar, and the system is stable on 10 Myr timescales. We find evidence of infrequent timing outliers, which we cannot explain with an eighth planet; we instead account for the outliers using a robust likelihood function. We forecast JWST timing observations and speculate on possible implications of the planet densities forthe formation, migration, and evolution of the planet system.

Eric Agol↗
NASA NTRS2020

Multiple Transits during a Single Conjunction: Identifying Transiting Circumbinary Planetary Candidates from TESS

We present results of a study on identifying circumbinary planet candidates that produce multiple transits during one conjunction with eclipsing binary systems. The occurrence of these transits enables us to estimate the candidates’ orbital periods, which is crucial as the periods of the currently known transiting circumbinary planets are significantly longer than the typical observational baseline of the Transiting Exoplanet Survey Satellite (TESS). Combined with the derived radii, it also provides valuable information needed for follow-up observations and subsequent confirmation of a large number of circumbinary planet candidates from TESS. Motivated by the discovery of the 1108 day circumbinary planet Kepler-1647, we show the application of this technique to four of Kepler’s circumbinary planets that produce such transits. Our results indicate that in systems where the circumbinary planet is on a low-eccentricity orbit, the estimated planetary orbital period is within <10%–20% of the true value. This estimate is derived from photometric observations spanning less than 5% of the planet’s period, demonstrating the strong capability of the technique. Capitalizing on the current and future eclipsing binaries monitored by NASA’s TESS mission, we estimate that hundreds of circumbinary planet candidates producing multiple transits during one conjunction will be detected in the TESS data. Such a large sample will enable statistical understanding of the population of planets orbiting binary stars and shed new light on their formation and evolution.

Veselin Kostov↗
NASA NTRS2016

EVEREST: Pixel Level Decorrelation of K2 Light Curves

We present EPIC Variability Extraction and Removal for Exoplanet Science Targets (EVEREST), an open-source pipeline for removing instrumental noise from K2 light curves. EVEREST employs a variant of pixel level decorrelation to remove systematics introduced by the spacecraft’s pointing error and a Gaussian process to capture astrophysical variability. We apply EVEREST to all K2 targets in campaigns 0–7, yielding light curves with precision comparable to that of the original Kepler mission for stars brighter than K(sub p) ≈ 13, and within a factor of two of the Kepler precision for fainter targets. We perform cross-validation and transit injection and recovery tests to validate the pipeline and compare our light curves to the other de-trended light curves available for download at the MAST High Level Science Products archive. We find that EVEREST achieves the highest average precision of any of these pipelines for unsaturated K2 stars. The improved precision of these light curves will aid in exoplanet detection and characterization, investigations of stellar variability, asteroseismology, and other photometric studies. The EVEREST pipeline can also easily be applied to future surveys, such as the TESS mission, to correct for instrumental systematics and enable the detection of low signal-to-noise transiting exoplanets. The EVEREST light curves and the source code used to generate them are freely available online.

Catalogs↗
NASA NTRS2016

Transit Timing Variations for Planets Near Eccentricity-Type Mean Motion Resonances

We derive the transit timing variations (TTVs) of two planets near a second-order mean motion resonance (MMR) on nearly circular orbits. We show that the TTVs of each planet are given by sinusoids with a frequency of jn(sub 2) - ( j - 2)n(sub 1), where j ≽ 3 is an integer characterizing the resonance and n(sub 2) and n(sub 1) are the mean motions of the outer and inner planets, respectively. The amplitude of the TTV depends on the mass of the perturbing planet, relative to the mass of the star, and on both the eccentricities and longitudes of pericenter of each planet. The TTVs of the two planets are approximated anti-correlated, with phases of ϕ and ≈ϕ + π, where the phase ϕ also depends on the eccentricities and longitudes of pericenter. Therefore, the TTVs caused by proximity to a second-order MMR do not in general uniquely determine both planet masses, eccentricities, and pericenters. This is completely analogous to the case of TTVs induced by two planets near a first-order MMR. We explore how other TTV signals, such as the short-period synodic TTV or a first-order resonant TTV, in combination with the second-order resonant TTV, can break degeneracies. Finally, we derive approximate formulae for the TTVs of planets near any order eccentricity-type MMR; this shows that the same basic sinusoidal TTV structure holds for all eccentricity-type resonances. Our general formula reduces to previously derived results near first-order MMRs.

Celestial mechanics↗
NASA NTRS2016

Secure Mass Measurements From Transit Timing: 10 Kepler Exoplanets Between 3 and 8M⊕ with Diverse Densities and Incident Fluxes

We infer dynamical masses in eight multiplanet systems using transit times measured from Keplerʼs complete dataset, including short-cadence data where available. Of the 18 dynamical masses that we infer, 10 pass multiple tests for robustness. These are in systems Kepler-26 (KOI-250), Kepler-29 (KOI-738), Kepler-60 (KOI-2086), Kepler-105 (KOI-115), and Kepler-307 (KOI-1576). Kepler -105 c has a radius of 1.3R⊕ and a density consistent with an Earth-like composition. Strong transit timing variation (TTV) signals were detected from additional planets, but their inferred masses were sensitive to outliers or consistent solutions could not be found with independently measured transit times, including planets orbiting Kepler-49 (KOI-248), Kepler-57 (KOI-1270), Kepler-105 (KOI-115), and Kepler-177 (KOI-523). Nonetheless, strong upper limits on the mass of Kepler-177 cimply an extremely low density of 0.1 g cm−3. In most cases, individual orbital eccentricities were poorly constrained owing to degeneracies in TTV inversion. For five planet pairs in our sample, strong secular interactions imply a moderate to high likelihood of apsidal alignment over a wide range of possible eccentricities. We also find solutions for the three planets known to orbit Kepler-60 in a Laplace-like resonance chain. However, nonlibrating solutions also match the transit timing data. For six systems, we calculate more precise stellar parameters than previously known, enabling useful constraints on planetary densities where we have secure mass measurements. Placing these exoplanets on the mass–radius diagram, we find that a wide range of densities is observed among sub-Neptune-mass planets and that the range in observed densities is anticorrelated with incident flux.

Stars individual↗
NASA NTRS2016

Secure Mass Measurements From Transit Timing: 10 Kepler Exoplanets Between 3 and 8 M⊕ With Diverse Densities and Incident Fluxes

We infer dynamical masses in eight multiplanet systems using transit times measured from Keplerʼs complete data set, including short-cadence data where available. Of the 18 dynamical masses that we infer, 10 pass multiple tests for robustness. These are in systems Kepler-26 (KOI-250), Kepler-29 (KOI-738), Kepler-60 (KOI-2086), Kepler-105 (KOI-115), and Kepler-307 (KOI-1576). Kepler-105 c has a radius of 1.3 R(sub ⨁) and a density consistent with an Earth-like composition. Strong transit timing variation (TTV) signals were detected from additional planets, but their inferred masses were sensitive to outliers or consistent solutions could not be found with independently measured transit times, including planets orbiting Kepler-49 (KOI-248), Kepler-57 (KOI-1270), Kepler-105 (KOI-115), and Kepler-177 (KOI-523). Nonetheless, strong upper limits on the mass of Kepler-177 c imply an extremely low density of ∽0.1 g cm(exp −3). In most cases, individual orbital eccentricities were poorly constrained owing to degeneracies in TTV inversion. For five planet pairs in our sample, strong secular interactions imply a moderate to high likelihood of apsidal alignment over a wide range of possible eccentricities. We also find solutions for the three planets known to orbit Kepler-60 in a Laplace-like resonance chain. However, nonlibrating solutions also match the transit timing data. For six systems, we calculate more precise stellar parameters than previously known, enabling useful constraints on planetary densities where we have secure mass measurements. Placing these exoplanets on the mass–radius diagram, we find that a wide range of densities is observed among sub-Neptune-mass planets and that the range in observed densities is anticorrelated with incident flux.

Individual stars↗
NASA NTRS2015

Spectral Eclipse Timing

We utilize multi-dimensional simulations of varying equatorial jet strength to predict wavelength-dependent variations in the eclipse times of gas-giant planets. A displaced hot spot introduces an asymmetry in the secondary eclipse light curve that manifests itself as a measured offset in the timing of the center of eclipse. A multiwavelength observation of secondary eclipse, one probing the timing of barycentric eclipse at short wavelengths and another probing at longer wavelengths, will reveal the longitudinal displacement of the hot spot and break the degeneracy between this effect and that associated with the asymmetry due to an eccentric orbit. The effect of time offsets was first explored in the IRAC wavebands by Williams et al. Here we improve upon their methodology, extend to a broad range of wavelengths, and demonstrate our technique on a series of multi-dimensional radiative-hydrodynamical simulations of HD 209458b with varying equatorial jet strength and hot-spot displacement. Simulations with the largest hot-spot displacement result in timing offsets of up to 100 s in the infrared. Though we utilize a particular radiative hydrodynamical model to demonstrate this effect, the technique is model independent. This technique should allow a much larger survey of hot-spot displacements with the James Webb Space Telescope than currently accessible with time-intensive phase curves, hopefully shedding light on the physical mechanisms associated with thermal energy advection in irradiated gas giants.

Hydrodynamics↗