Showing 1–50 of 228 results for author: Agol, E

The Implications of Thermal Hydrodynamic Atmospheric Escape on the TRAPPIST-1 Planets

Authors: Megan T. Gialluca, Rory Barnes, Victoria S. Meadows, Rodolfo Garcia, Jessica Birky, Eric Agol

Abstract: JWST observations of the 7-planet TRAPPIST-1 system will provide an excellent opportunity to test outcomes of stellar-driven evolution of terrestrial planetary atmospheres, including atmospheric escape, ocean loss and abiotic oxygen production. While most previous studies use a single luminosity evolution for the host star, we incorporate observational uncertainties in stellar mass, luminosity evo… ▽ More JWST observations of the 7-planet TRAPPIST-1 system will provide an excellent opportunity to test outcomes of stellar-driven evolution of terrestrial planetary atmospheres, including atmospheric escape, ocean loss and abiotic oxygen production. While most previous studies use a single luminosity evolution for the host star, we incorporate observational uncertainties in stellar mass, luminosity evolution, system age, and planetary parameters to statistically explore the plausible range of planetary atmospheric escape outcomes. We present probabilistic distributions of total water loss and oxygen production as a function of initial water content, for planets with initially pure water atmospheres and no interior-atmosphere exchange. We find that the interior planets are desiccated for initial water contents below 50 Earth oceans. For TRAPPIST-1e, f, g, and h, we report maximum water loss ranges of 8.0$^{+1.3}_{-0.9}$, 4.8$^{+0.6}_{-0.4}$, 3.4$^{+0.3}_{-0.3}$, and 0.8$^{+0.2}_{-0.1}$ Earth oceans, respectively, with corresponding maximum oxygen retention of 1290$^{+75}_{-75}$, 800$^{+40}_{-40}$, 560$^{+30}_{-25}$, and 90$^{+10}_{-10}$ bars. We explore statistical constraints on initial water content imposed by current water content, which could inform evolutionary history and planet formation. If TRAPPIST-1b is airless while TRAPPIST-1c possesses a tenuous oxygen atmosphere, as initial JWST observations suggest, then our models predict an initial surface water content of 8.2$^{+1.5}_{-1.0}$ Earth oceans for these worlds, leading to the outer planets retaining $>$1.5 Earth oceans after entering the habitable zone. Even if TRAPPIST-1c is airless, surface water on the outer planets would not be precluded. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 33 pages including appendix. Accepted for Publication in the Planetary Science Journal

Multiple Patchy Cloud Layers in the Planetary Mass Object SIMP0136+0933

Authors: Allison M. McCarthy, Philip S. Muirhead, Patrick Tamburo, Johanna M. Vos, Caroline V. Morley, Jacqueline Faherty, Daniella C. Bardalez Gagliuffi, Eric Agol, Christopher Theissen

Abstract: Multi-wavelength photometry of brown dwarfs and planetary-mass objects provides insight into their atmospheres and cloud layers. We present near-simultaneous $J-$ and $K_s-$band multi-wavelength observations of the highly variable T2.5 planetary-mass object, SIMP J013656.5+093347. We reanalyze observations acquired over a single night in 2015 using a recently developed data reduction pipeline. For… ▽ More Multi-wavelength photometry of brown dwarfs and planetary-mass objects provides insight into their atmospheres and cloud layers. We present near-simultaneous $J-$ and $K_s-$band multi-wavelength observations of the highly variable T2.5 planetary-mass object, SIMP J013656.5+093347. We reanalyze observations acquired over a single night in 2015 using a recently developed data reduction pipeline. For the first time, we detect a phase shift between $J-$ and $K_s-$band light curves, which we measure to be $39.9^{\circ +3.6}_{ -1.1}$. Previously, phase shifts between near-infrared and mid-infrared observations of this object were detected and attributed to probing different depths of the atmosphere, and thus different cloud layers. Using the Sonora Bobcat models, we expand on this idea to show that at least two different patchy cloud layers must be present to explain the measured phase shift. Our results are generally consistent with recent atmospheric retrievals of this object and other similar L/T transition objects. △ Less

Submitted 26 February, 2024; v1 submitted 22 February, 2024; originally announced February 2024.

Comments: 15 pages, 8 figures, Accepted for publication in ApJ

A roadmap to the efficient and robust characterization of temperate terrestrial planet atmospheres with JWST

Authors: TRAPPIST-1 JWST Community Initiative, :, Julien de Wit, René Doyon, Benjamin V. Rackham, Olivia Lim, Elsa Ducrot, Laura Kreidberg, Björn Benneke, Ignasi Ribas, David Berardo, Prajwal Niraula, Aishwarya Iyer, Alexander Shapiro, Nadiia Kostogryz, Veronika Witzke, Michaël Gillon, Eric Agol, Victoria Meadows, Adam J. Burgasser, James E. Owen, Jonathan J. Fortney, Franck Selsis, Aaron Bello-Arufe, Emeline Bolmont , et al. (57 additional authors not shown)

Abstract: Ultra-cool dwarf stars are abundant, long-lived, and uniquely suited to enable the atmospheric study of transiting terrestrial companions with JWST. Amongst them, the most prominent is the M8.5V star TRAPPIST-1 and its seven planets, which have been the favored targets of eight JWST Cycle 1 programs. While Cycle 1 observations have started to yield preliminary insights into the planets, they have… ▽ More Ultra-cool dwarf stars are abundant, long-lived, and uniquely suited to enable the atmospheric study of transiting terrestrial companions with JWST. Amongst them, the most prominent is the M8.5V star TRAPPIST-1 and its seven planets, which have been the favored targets of eight JWST Cycle 1 programs. While Cycle 1 observations have started to yield preliminary insights into the planets, they have also revealed that their atmospheric exploration requires a better understanding of their host star. Here, we propose a roadmap to characterize the TRAPPIST-1 system -- and others like it -- in an efficient and robust manner. We notably recommend that -- although more challenging to schedule -- multi-transit windows be prioritized to constrain stellar heterogeneities and gather up to 2$\times$ more transits per JWST hour spent. We conclude that in such systems planets cannot be studied in isolation by small programs, thus large-scale community-supported programs should be supported to enable the efficient and robust exploration of terrestrial exoplanets in the JWST era. △ Less

Submitted 24 October, 2023; v1 submitted 24 October, 2023; originally announced October 2023.

Potential Atmospheric Compositions of TRAPPIST-1 c constrained by JWST/MIRI Observations at 15 $μみゅー$m

Authors: Andrew P. Lincowski, Victoria S. Meadows, Sebastian Zieba, Laura Kreidberg, Caroline Morley, Michaël Gillon, Franck Selsis, Eric Agol, Emeline Bolmont, Elsa Ducrot, Renyu Hu, Daniel D. B. Koll, Xintong Lyu, Avi Mandell, Gabrielle Suissa, Patrick Tamburo

Abstract: The first JWST observations of TRAPPIST-1 c showed a secondary eclipse depth of 421+/-94 ppm at 15 um, which is consistent with a bare rock surface or a thin, O2-dominated, low CO2 atmosphere (Zieba et al. 2023). Here, we further explore potential atmospheres for TRAPPIST-1 c by comparing the observed secondary eclipse depth to synthetic spectra of a broader range of plausible environments. To sel… ▽ More The first JWST observations of TRAPPIST-1 c showed a secondary eclipse depth of 421+/-94 ppm at 15 um, which is consistent with a bare rock surface or a thin, O2-dominated, low CO2 atmosphere (Zieba et al. 2023). Here, we further explore potential atmospheres for TRAPPIST-1 c by comparing the observed secondary eclipse depth to synthetic spectra of a broader range of plausible environments. To self-consistently incorporate the impact of photochemistry and atmospheric composition on atmospheric thermal structure and predicted eclipse depth, we use a two-column climate model coupled to a photochemical model, and simulate O2-dominated, Venus-like, and steam atmospheres. We find that a broader suite of plausible atmospheric compositions are also consistent with the data. For lower pressure atmospheres (0.1 bar), our O2-CO2 atmospheres produce eclipse depths within 1$σしぐま$ of the data, consistent with the modeling results of Zieba et al. (2023). However, for higher-pressure atmospheres, our models produce different temperature-pressure profiles and are less pessimistic, with 1-10 bar O2, 100 ppm CO2 models within 2.0-2.2$σしぐま$ of the measured secondary eclipse depth, and up to 0.5% CO2 within 2.9$σしぐま$. Venus-like atmospheres are still unlikely. For thin O2 atmospheres of 0.1 bar with a low abundance of CO2 ($\sim$100 ppm), up to 10% water vapor can be present and still provide an eclipse depth within 1$σしぐま$ of the data. We compared the TRAPPIST-1 c data to modeled steam atmospheres of $\leq$ 3 bar, which are 1.7-1.8$σしぐま$ from the data and not conclusively ruled out. More data will be required to discriminate between possible atmospheres, or to more definitively support the bare rock hypothesis. △ Less

Submitted 10 August, 2023; originally announced August 2023.

Comments: 15 pages, accepted to APJL

No thick carbon dioxide atmosphere on the rocky exoplanet TRAPPIST-1 c

Authors: Sebastian Zieba, Laura Kreidberg, Elsa Ducrot, Michaël Gillon, Caroline Morley, Laura Schaefer, Patrick Tamburo, Daniel D. B. Koll, Xintong Lyu, Lorena Acuña, Eric Agol, Aishwarya R. Iyer, Renyu Hu, Andrew P. Lincowski, Victoria S. Meadows, Franck Selsis, Emeline Bolmont, Avi M. Mandell, Gabrielle Suissa

Abstract: Seven rocky planets orbit the nearby dwarf star TRAPPIST-1, providing a unique opportunity to search for atmospheres on small planets outside the Solar System (Gillon et al., 2017). Thanks to the recent launch of JWST, possible atmospheric constituents such as carbon dioxide (CO2) are now detectable (Morley et al., 2017, Lincowski et al., 2018}. Recent JWST observations of the innermost planet TRA… ▽ More Seven rocky planets orbit the nearby dwarf star TRAPPIST-1, providing a unique opportunity to search for atmospheres on small planets outside the Solar System (Gillon et al., 2017). Thanks to the recent launch of JWST, possible atmospheric constituents such as carbon dioxide (CO2) are now detectable (Morley et al., 2017, Lincowski et al., 2018}. Recent JWST observations of the innermost planet TRAPPIST-1 b showed that it is most probably a bare rock without any CO2 in its atmosphere (Greene et al., 2023). Here we report the detection of thermal emission from the dayside of TRAPPIST-1 c with the Mid-Infrared Instrument (MIRI) on JWST at 15 micron. We measure a planet-to-star flux ratio of fp/fs = 421 +/- 94 parts per million (ppm) which corresponds to an inferred dayside brightness temperature of 380 +/- 31 K. This high dayside temperature disfavours a thick, CO2-rich atmosphere on the planet. The data rule out cloud-free O2/CO2 mixtures with surface pressures ranging from 10 bar (with 10 ppm CO2) to 0.1 bar (pure CO2). A Venus-analogue atmosphere with sulfuric acid clouds is also disfavoured at 2.6 sigma confidence. Thinner atmospheres or bare-rock surfaces are consistent with our measured planet-to-star flux ratio. The absence of a thick, CO2-rich atmosphere on TRAPPIST-1 c suggests a relatively volatile-poor formation history, with less than 9.5 +7.5 -2.3 Earth oceans of water. If all planets in the system formed in the same way, this would indicate a limited reservoir of volatiles for the potentially habitable planets in the system. △ Less

Submitted 16 June, 2023; originally announced June 2023.

Comments: Published in Nature on June 19th. 2023, 10 figures, 4 tables

The Perkins INfrared Exosatellite Survey (PINES) II. Transit Candidates and Implications for Planet Occurrence around L and T Dwarfs

Authors: Patrick Tamburo, Philip S. Muirhead, Allison M. McCarthy, Murdock Hart, Johanna M. Vos, Eric Agol, Christopher Theissen, David Gracia, Daniella C. Bardalez Gagliuffi, Jacqueline Faherty

Abstract: We describe a new transit detection algorithm designed to detect single transit events in discontinuous Perkins INfrared Exosatellite Survey (PINES) observations of L and T dwarfs. We use this algorithm to search for transits in 131 PINES light curves and identify two transit candidates: 2MASS J18212815+1414010 (2MASS J1821+1414) and 2MASS J08350622+1953050 (2MASS J0835+1953). We disfavor 2MASS J1… ▽ More We describe a new transit detection algorithm designed to detect single transit events in discontinuous Perkins INfrared Exosatellite Survey (PINES) observations of L and T dwarfs. We use this algorithm to search for transits in 131 PINES light curves and identify two transit candidates: 2MASS J18212815+1414010 (2MASS J1821+1414) and 2MASS J08350622+1953050 (2MASS J0835+1953). We disfavor 2MASS J1821+1414 as a genuine transit candidate due to the known variability properties of the source. We cannot rule out the planetary nature of 2MASS J0835+1953's candidate event and perform follow-up observations in an attempt to recover a second transit. A repeat event has yet to be observed, but these observations suggest that target variability is an unlikely cause of the candidate transit. We perform a Markov chain Monte Carlo simulation of the light curve and estimate a planet radius ranging from $4.2^{+3.5}_{-1.6}R_\oplus$ to $5.8^{+4.8}_{-2.1}R_\oplus$, depending on the host's age. Finally, we perform an injection and recovery simulation on our light curve sample. We inject planets into our data using measured M dwarf planet occurrence rates and attempt to recover them using our transit search algorithm. Our detection rates suggest that, assuming M dwarf planet occurrence rates, we should have roughly a 1$\%$ chance of detecting a candidate that could cause the transit depth we observe for 2MASS J0835+1953. If 2MASS J0835+1953 b is confirmed, it would suggest an enhancement in the occurrence of short-period planets around L and T dwarfs in comparison to M dwarfs, which would challenge predictions from planet formation models. △ Less

Submitted 10 October, 2022; originally announced October 2022.

Comments: 23 pages, 15 figures, accepted to AJ

Analytic Light Curve for Mutual Transits of Two Bodies Across a Limb-darkened Star

Authors: Tyler A. Gordon, Eric Agol

Abstract: We present a solution for the light curve of two bodies mutually transiting a star with polynomial limb darkening. The term "mutual transit" in this work refers to a transit of the star during which overlap occurs between the two transiting bodies. These could be an exoplanet with an exomoon companion, two exoplanets, an eclipsing binary and a planet, or two stars eclipsing a third in a triple sta… ▽ More We present a solution for the light curve of two bodies mutually transiting a star with polynomial limb darkening. The term "mutual transit" in this work refers to a transit of the star during which overlap occurs between the two transiting bodies. These could be an exoplanet with an exomoon companion, two exoplanets, an eclipsing binary and a planet, or two stars eclipsing a third in a triple star system. We include analytic derivatives of the light curve with respect to the positions and radii of both bodies. We provide code that implements a photodynamical model for a mutual transit. We include two dynamical models, one for hierarchical systems in which a secondary body orbits a larger primary (e.g. an exomoon system) and a second for confocal systems in which two bodies independently orbit a central mass (e.g. two planets in widely separated orbits). Our code is fast enough to enable inference with MCMC algorithms, and the inclusion of derivatives allows for the use of gradient-based inference methods such as Hamiltonian Monte Carlo. While applicable to a variety of systems, this work was undertaken primarily with exomoons in mind. It is our hope that making this code publicly available will reduce barriers for the community to assess the detectability of exomoons, conduct searches for exomoons, and attempt to validate existing exomoon candidates. We also anticipate that our code will be useful for studies of planet-planet transits in exoplanetary systems, transits of circumbinary planets, and eclipses in triple-star systems. △ Less

Submitted 13 July, 2022; originally announced July 2022.

Comments: 29 pages, 18 figures. Code available at https://github.com/tagordon/gefera

Llamaradas Estelares: Modeling the Morphology of White-Light Flares

Authors: Guadalupe Tovar Mendoza, James R. A. Davenport, Eric Agol, James A. G. Jackman, Suzanne L. Hawley

Abstract: Stellar variability is a limiting factor for planet detection and characterization, particularly around active M-type stars. Here we revisit one of the most active stars from the Kepler mission, the M4 star GJ 1243, and use a sample of 414 flare events from 11 months of 1-minute cadence light curves to study the empirical morphology of white-light stellar flares. We use a Gaussian process detrendi… ▽ More Stellar variability is a limiting factor for planet detection and characterization, particularly around active M-type stars. Here we revisit one of the most active stars from the Kepler mission, the M4 star GJ 1243, and use a sample of 414 flare events from 11 months of 1-minute cadence light curves to study the empirical morphology of white-light stellar flares. We use a Gaussian process detrending technique to account for the underlying starspots. We present an improved analytic, continuous flare template that is generated by stacking the flares onto a scaled time and amplitude and uses a Markov Chain Monte Carlo analysis to fit the model. Our model is defined using classical flare events, but can also be used to model complex, multi-peaked flare events. We demonstrate the utility of our model using TESS data at the 10-minute, 2-minute and 20-second cadence modes. Our new flare model code is made publicly available on GitHub. △ Less

Submitted 11 May, 2022; originally announced May 2022.

Comments: 20 pages, 10 figures, accepted for publication in AJ. New flare model is available online https://github.com/lupitatovar/Llamaradas-Estelares

The Perkins INfrared Exosatellite Survey (PINES) I. Survey Overview, Reduction Pipeline, and Early Results

Authors: Patrick Tamburo, Philip S. Muirhead, Allison M. McCarthy, Murdock Hart, David Gracia, Johanna M. Vos, Daniella C. Bardalez Gagliuffi, Jacqueline Faherty, Christopher Theissen, Eric Agol, Julie N. Skinner, Sheila Sagear

Abstract: We describe the Perkins INfrared Exosatellite Survey (PINES), a near-infrared photometric search for short-period transiting planets and moons around a sample of 393 spectroscopically confirmed L- and T-type dwarfs. PINES is performed with Boston University's 1.8 m Perkins Telescope Observatory, located on Anderson Mesa, Arizona. We discuss the observational strategy of the survey, which was desig… ▽ More We describe the Perkins INfrared Exosatellite Survey (PINES), a near-infrared photometric search for short-period transiting planets and moons around a sample of 393 spectroscopically confirmed L- and T-type dwarfs. PINES is performed with Boston University's 1.8 m Perkins Telescope Observatory, located on Anderson Mesa, Arizona. We discuss the observational strategy of the survey, which was designed to optimize the number of expected transit detections, and describe custom automated observing procedures for performing PINES observations. We detail the steps of the $\texttt{PINES Analysis Toolkit}$ ($\texttt{PAT}$), software that is used to create light curves from PINES images. We assess the impact of second-order extinction due to changing precipitable water vapor on our observations and find that the magnitude of this effect is minimized in Mauna Kea Observatories $\textit{J}$-band. We demonstrate the validity of $\texttt{PAT}$ through the recovery of a transit of WASP-2 b and known variable brown dwarfs, and use it to identify a new variable L/T transition object: the T2 dwarf WISE J045746.08-020719.2. We report on the measured photometric precision of the survey and use it to estimate our transit detection sensitivity. We find that for our median brightness targets, assuming contributions from white noise only, we are sensitive to the detection of 2.5 $R_\oplus$ planets and larger. PINES will test whether the increase in sub-Neptune-sized planet occurrence with decreasing host mass continues into the L and T dwarf regime. △ Less

Submitted 21 April, 2022; v1 submitted 5 January, 2022; originally announced January 2022.

Comments: 25 pages, 15 figures, accepted to AJ

An upper limit on late accretion and water delivery in the Trappist-1 exoplanet system

Authors: Sean N. Raymond, Andre Izidoro, Emeline Bolmont, Caroline Dorn, Franck Selsis, Martin Turbet, Eric Agol, Patrick Barth, Ludmila Carone, Rajdeep Dasgupta, Michael Gillon, Simon L. Grimm

Abstract: The Trappist-1 system contains seven roughly Earth-sized planets locked in a multi-resonant orbital configuration, which has enabled precise measurements of the planets' masses and constrained their compositions. Here we use the system's fragile orbital structure to place robust upper limits on the planets' bombardment histories. We use N-body simulations to show how perturbations from additional… ▽ More The Trappist-1 system contains seven roughly Earth-sized planets locked in a multi-resonant orbital configuration, which has enabled precise measurements of the planets' masses and constrained their compositions. Here we use the system's fragile orbital structure to place robust upper limits on the planets' bombardment histories. We use N-body simulations to show how perturbations from additional objects can break the multi-resonant configuration by either triggering dynamical instability or simply removing the planets from resonance. The planets cannot have interacted with more than ${\sim 5\%}$ of an Earth mass (${M_\oplus}$) in planetesimals -- or a single rogue planet more massive than Earth's Moon -- without disrupting their resonant orbital structure. This implies an upper limit of ${10^{-4}}$ to ${10^{-2} M_\oplus}$ of late accretion on each planet since the dispersal of the system's gaseous disk. This is comparable to or less than the late accretion on Earth after the Moon-forming impact, and demonstrates that the Trappist-1 planets' growth was complete in just a few million years, roughly an order of magnitude faster than Earth's. Our results imply that any large water reservoirs on the Trappist-1 planets must have been incorporated during their formation in the gaseous disk. △ Less

Submitted 26 November, 2021; originally announced November 2021.

Comments: Published in Nature Astronomy (Nov 25, 2021). This is the authors' version including Methods and Supplementary Info

Thermal Phase Curves of XO-3b: an Eccentric Hot Jupiter at the Deuterium Burning Limit

Authors: Lisa Dang, Taylor J. Bell, Nicolas B. Cowan, Daniel Thorngren, Tiffany Kataria, Heather A. Knutson, Nikole K. Lewis, Keivan G. Stassun, Jonathan J. Fortney, Eric Agol, Gregory P. Laughlin, Adam Burrows, Karen A. Collins, Drake Deming, Diana Jovmir, Jonathan Langton, Sara Rastegar, Adam P. Showman

Abstract: We report \textit{Spitzer} full-orbit phase observations of the eccentric hot Jupiter XO-3b at 3.6 and 4.5 $μみゅー$m. Our new eclipse depth measurements of $1770 \pm 180$ ppm at 3.6 $μみゅー$m and $1610 \pm 70$ ppm at 4.5 $μみゅー$m show no evidence of the previously reported dayside temperature inversion. We also empirically derive the mass and radius of XO-3b and its host star using Gaia DR3's parallax measureme… ▽ More We report \textit{Spitzer} full-orbit phase observations of the eccentric hot Jupiter XO-3b at 3.6 and 4.5 $μみゅー$m. Our new eclipse depth measurements of $1770 \pm 180$ ppm at 3.6 $μみゅー$m and $1610 \pm 70$ ppm at 4.5 $μみゅー$m show no evidence of the previously reported dayside temperature inversion. We also empirically derive the mass and radius of XO-3b and its host star using Gaia DR3's parallax measurement and find a planetary mass $M_p=11.79 \pm 0.98 ~M_{\rm{Jup}}$ and radius $R_p=1.295 \pm 0.066 ~R_{\rm{Jup}}$. We compare our \textit{Spitzer} observations with multiple atmospheric models to constrain the radiative and advective properties of XO-3b. While the decorrelated 4.5 $μみゅー$m observations are pristine, the 3.6 $μみゅー$m phase curve remains polluted with detector systematics due to larger amplitude intrapixel sensitivity variations in this channel. We focus our analysis on the more reliable 4.5 $μみゅー$m phase curve and fit an energy balance model with solid body rotation to estimate the zonal wind speed and the pressure of the bottom of the mixed layer. Our energy balance model fit suggests an eastward equatorial wind speed of $3.13 ^{+0.26} _{-0.83}$ km/s, an atmospheric mixed layer down to $2.40 ^{+0.92} _{-0.16}$ bar, and Bond albedo of $0.106 ^{+0.008} _{-0.106}$. We assume that the wind speed and mixed layer depth are constant throughout the orbit. We compare our observations with a 1D planet-averaged model predictions at apoapse and periapse and 3D general circulation model (GCM) predictions for XO-3b. We also investigate the inflated radius of XO-3b and find that it would require an unusually large amount of internal heating to explain the observed planetary radius. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 23 pages, 14 figures, accepted for publication in AJ

TRAPPIST-1: Dynamical analysis of the transit-timing variations and origin of the resonant chain

Authors: Jean Teyssandier, Anne-Sophie Libert, Eric Agol

Abstract: We analyze solutions drawn from the recently published posterior distribution of the TRAPPIST-1 system, which consists of seven Earth-size planets appearing to be in a resonant chain around a red dwarf. We show that all the planets are simultaneously in two-planet and three-planet resonances, apart from the innermost pair for which the two-planet resonant angles circulate. By means of a frequency… ▽ More We analyze solutions drawn from the recently published posterior distribution of the TRAPPIST-1 system, which consists of seven Earth-size planets appearing to be in a resonant chain around a red dwarf. We show that all the planets are simultaneously in two-planet and three-planet resonances, apart from the innermost pair for which the two-planet resonant angles circulate. By means of a frequency analysis, we highlight that the transit-timing variation (TTV) signals possess a series of common periods varying from days to decades, which are also present in the variations of the dynamical variables of the system. Shorter periods (e.g., the TTVs characteristic timescale of 1.3 yr) are associated with two-planet mean-motion resonances, while longer periods arise from three-planet resonances. By use of $N$-body simulations with migration forces, we explore the origin of the resonant chain of TRAPPIST-1 and find that for particular disc conditions, a chain of resonances - similar to the observed one - can be formed which accurately reproduces the observed TTVs. Our analysis suggests that while the 4-yr collected data of observations hold key information on the two-planet resonant dynamics, further monitoring of TRAPPIST-1 will soon provide signatures of three-body resonances, in particular the 3.3 and 5.1 yr periodicities expected for the current best-fit solution. dditional observations would help to assess whether the innermost pair of planets is indeed resonant (its proximity to the 8:5 resonance being challenging to explain), and therefore give additional constraints on formation scenarios. △ Less

Submitted 10 February, 2022; v1 submitted 7 October, 2021; originally announced October 2021.

Comments: 11 pages, 14 figures, accepted for publication in A&A

Journal ref: A&A 658, A170 (2022)

A differentiable N-body code for transit timing and dynamical modeling. I. Algorithm and derivatives

Authors: Eric Agol, David M. Hernandez, Zachary Langford

Abstract: When fitting N-body models to astronomical data - including transit times, radial velocity, and astrometric positions at observed times - the derivatives of the model outputs with respect to the initial conditions can help with model optimization and posterior sampling. Here we describe a general-purpose symplectic integrator for arbitrary orbital architectures, including those with close encounte… ▽ More When fitting N-body models to astronomical data - including transit times, radial velocity, and astrometric positions at observed times - the derivatives of the model outputs with respect to the initial conditions can help with model optimization and posterior sampling. Here we describe a general-purpose symplectic integrator for arbitrary orbital architectures, including those with close encounters, which we have recast to maintain numerical stability and precision for small step sizes. We compute the derivatives of the N-body coordinates and velocities as a function of time with respect to the initial conditions and masses by propagating the Jacobian along with the N-body integration. For the first time we obtain the derivatives of the transit times with respect to the initial conditions and masses using the chain rule, which is quicker and more accurate than using finite differences or automatic differentiation. We implement this algorithm in an open source package, NbodyGradient.jl, written in the Julia language, which has been used in the optimization and error analysis of transit-timing variations in the TRAPPIST-1 system. We present tests of the accuracy and precision of the code, and show that it compares favorably in speed to other integrators which are written in C. △ Less

Submitted 3 June, 2021; originally announced June 2021.

Comments: 24 pages; 7 figures. Revised for submission to MNRAS. Links in equations are provided to code at the repository https://github.com/ericagol/NbodyGradient.jl Some documentation and example scripts for running the code are provided in the repository

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

Authors: Veselin B. Kostov, Brian P. Powell, Jerome A. Orosz, William F. Welsh, William Cochran, Karen A. Collins, Michael Endl, Coel Hellier, David W. Latham, Phillip MacQueen, Joshua Pepper, Billy Quarles, Lalitha Sairam, Guillermo Torres, Robert F. Wilson, Serge Bergeron, Pat Boyce, Allyson Bieryla, Robert Buchheim, Caleb Ben Christiansen, David R. Ciardi, Kevin I. Collins, Dennis M. Conti, Scott Dixon, Pere Guerra , et al. (64 additional authors not shown)

Abstract: We report the first discovery of a transiting circumbinary planet detected from a single sector of TESS data. During Sector 21, the planet TIC 172900988b transited the primary star and then 5 days later it transited the secondary star. The binary is itself eclipsing, with a period of P = 19.7 days and an eccentricity of e = 0.45. Archival data from ASAS-SN, Evryscope, KELT, and SuperWASP reveal a… ▽ More We report the first discovery of a transiting circumbinary planet detected from a single sector of TESS data. During Sector 21, the planet TIC 172900988b transited the primary star and then 5 days later it transited the secondary star. The binary is itself eclipsing, with a period of P = 19.7 days and an eccentricity of 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 M1 = 1.2384 +/- 0.0007 MSun and R1 = 1.3827 +/- 0.0016 RSun for the primary and M2 = 1.2019 +/- 0.0007 MSun and R2 = 1.3124 +/- 0.0012 RSun for the secondary. The radius of the planet is R3 = 11.25 +/- 0.44 REarth (1.004 +/- 0.039 RJup). 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 MEarth (2.65 < M3 < 3.09 MJup), 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 a V = 10.141 mag, the system is accessible for high-resolution spectroscopic observations, e.g. Rossiter-McLaughlin effect and transit spectroscopy. △ Less

Submitted 27 August, 2021; v1 submitted 18 May, 2021; originally announced May 2021.

Comments: 57 pages, 30 figures, 25 tables; Accepted AJ

exoplanet: Gradient-based probabilistic inference for exoplanet data & other astronomical time series

Authors: Daniel Foreman-Mackey, Rodrigo Luger, Eric Agol, Thomas Barclay, Luke G. Bouma, Timothy D. Brandt, Ian Czekala, Trevor J. David, Jiayin Dong, Emily A. Gilbert, Tyler A. Gordon, Christina Hedges, Daniel R. Hey, Brett M. Morris, Adrian M. Price-Whelan, Arjun B. Savel

Abstract: "exoplanet" is a toolkit for probabilistic modeling of astronomical time series data, with a focus on observations of exoplanets, using PyMC3 (Salvatier et al., 2016). PyMC3 is a flexible and high-performance model-building language and inference engine that scales well to problems with a large number of parameters. "exoplanet" extends PyMC3's modeling language to support many of the custom functi… ▽ More "exoplanet" is a toolkit for probabilistic modeling of astronomical time series data, with a focus on observations of exoplanets, using PyMC3 (Salvatier et al., 2016). PyMC3 is a flexible and high-performance model-building language and inference engine that scales well to problems with a large number of parameters. "exoplanet" extends PyMC3's modeling language to support many of the custom functions and probability distributions required when fitting exoplanet datasets or other astronomical time series. While it has been used for other applications, such as the study of stellar variability, the primary purpose of "exoplanet" is the characterization of exoplanets or multiple star systems using time-series photometry, astrometry, and/or radial velocity. In particular, the typical use case would be to use one or more of these datasets to place constraints on the physical and orbital parameters of the system, such as planet mass or orbital period, while simultaneously taking into account the effects of stellar variability. △ Less

Submitted 24 June, 2021; v1 submitted 5 May, 2021; originally announced May 2021.

Comments: Published in the Journal of Open Source Software. Comments (still) welcome. Software available at https://docs.exoplanet.codes

Significant improvement in planetary system simulations from statistical averaging

Authors: David M. Hernandez, Eric Agol, Matthew J. Holman, Sam Hadden

Abstract: Symplectic integrators are widely used for the study of planetary dynamics and other $N$-body problems. In a study of the outer Solar system, we demonstrate that individual symplectic integrations can yield biased errors in the semi-major axes and possibly other orbital elements. The bias is resolved by studying an ensemble of initial conditions of the outer Solar system. Such statistical sampling… ▽ More Symplectic integrators are widely used for the study of planetary dynamics and other $N$-body problems. In a study of the outer Solar system, we demonstrate that individual symplectic integrations can yield biased errors in the semi-major axes and possibly other orbital elements. The bias is resolved by studying an ensemble of initial conditions of the outer Solar system. Such statistical sampling could significantly improve measurement of planetary system properties like their secular frequencies. We also compared the distributions of action-like variables between high and low accuracy integrations; traditional statistical metrics are unable to distinguish the distribution functions. △ Less

Submitted 5 April, 2021; originally announced April 2021.

Comments: 4 pages, 1 figure. Accepted by RNAAS

Journal ref: 2021 Res. Notes AAS 5 77

Analytic Light Curves in Reflected Light: Phase Curves, Occultations, and Non-Lambertian Scattering for Spherical Planets and Moons

Authors: Rodrigo Luger, Eric Agol, Fran Bartolić, Daniel Foreman-Mackey

Abstract: We derive efficient, closed form, differentiable, and numerically stable solutions for the flux measured from a spherical planet or moon seen in reflected light, either in or out of occultation. Our expressions apply to the computation of scattered light phase curves of exoplanets, secondary eclipse light curves in the optical, or future measurements of planet-moon and planet-planet occultations,… ▽ More We derive efficient, closed form, differentiable, and numerically stable solutions for the flux measured from a spherical planet or moon seen in reflected light, either in or out of occultation. Our expressions apply to the computation of scattered light phase curves of exoplanets, secondary eclipse light curves in the optical, or future measurements of planet-moon and planet-planet occultations, as well as to photometry of solar system bodies. We derive our solutions for Lambertian bodies illuminated by a point source, but extend them to model illumination sources of finite angular size and rough surfaces with phase-dependent scattering. Our algorithm is implemented in Python within the open-source starry mapping framework and is designed with efficient gradient-based inference in mind. The algorithm is 4-5 orders of magnitude faster than direct numerical evaluation methods and about 10 orders of magnitude more precise. We show how the techniques developed here may one day lead to the construction of two-dimensional maps of terrestrial planet surfaces, potentially enabling the detection of continents and oceans on exoplanets in the habitable zone. △ Less

Submitted 10 March, 2021; originally announced March 2021.

Comments: 61 pages, 17 figures. To be submitted to AAS journals

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

Authors: Kevin K. Hardegree-Ullman, Jessie L. Christiansen, David R. Ciardi, Ian J. M. Crossfield, Courtney D. Dressing, John H. Livingston, Kathryn Volk, Eric Agol, Thomas Barclay, Geert Barentsen, Björn Benneke, Varoujan Gorjian, Martti H. Kristiansen

Abstract: $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… ▽ More $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\pm0.006$ days and measure a radius of $3.44^{+0.32}_{-0.31}\,R_{\oplus}$, 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 non-transiting 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. △ Less

Submitted 17 February, 2021; originally announced February 2021.

Comments: 17 pages, 10 figures, 2 tables, 1 fabulous new planet. Accepted for publication in AJ

Stellar Rotation in the K2 Sample: Evidence for Modified Spindown

Authors: Tyler A. Gordon, James R. A. Davenport, Ruth Angus, Daniel Foreman-Mackey, Eric Agol, Kevin R. Covey, Marcel Agüeros, David Kipping

Abstract: We analyze light curves of 284,834 unique K2 targets using a Gaussian process model with a quasi-periodic kernel function. By crossmatching K2 stars to observations from Gaia Data Release 2, we have identified 69,627 likely main-sequence stars. From these we select a subsample of 8,977 stars on the main-sequence with highly precise rotation period measurements. With this sample we recover the gap… ▽ More We analyze light curves of 284,834 unique K2 targets using a Gaussian process model with a quasi-periodic kernel function. By crossmatching K2 stars to observations from Gaia Data Release 2, we have identified 69,627 likely main-sequence stars. From these we select a subsample of 8,977 stars on the main-sequence with highly precise rotation period measurements. With this sample we recover the gap in the rotation period-color diagram first reported by \cite{McQuillan2013}. While the gap was tentatively detected in \cite{Reinhold2020}, this work represents the first robust detection of the gap in K2 data for field stars. This is significant because K2 observed along many lines of sight at wide angular separation, in contrast to Kepler's single line of sight. Together with recent results for rotation in open clusters, we interpret this gap as evidence for a departure from the $t^{-1/2}$ Skumanich spin down law, rather than an indication of a bimodal star formation history. We provide maximum likelihood estimates and uncertainties for all parameters of the quasi-periodic light curve model for each of the 284,834 stars in our sample. △ Less

Submitted 1 April, 2021; v1 submitted 19 January, 2021; originally announced January 2021.

Comments: 17 pages, 14 figures, submitted to ApJ, revised April 1, 2020

Refining the transit timing and photometric analysis of TRAPPIST-1: Masses, radii, densities, dynamics, and ephemerides

Authors: Eric Agol, Caroline Dorn, Simon L. Grimm, Martin Turbet, Elsa Ducrot, Laetitia Delrez, Michael Gillon, Brice-Olivier Demory, Artem Burdanov, Khalid Barkaoui, Zouhair Benkhaldoun, Emeline Bolmont, Adam Burgasser, Sean Carey, Julien de Wit, Daniel Fabrycky, Daniel Foreman-Mackey, Jonas Haldemann, David M. Hernandez, James Ingalls, Emmanuel Jehin, Zachary Langford, Jeremy Leconte, Susan M. Lederer, Rodrigo Luger , et al. (10 additional authors not shown)

Abstract: 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… ▽ More 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 star and 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 to a precision of 3-5%, equivalent to a radial-velocity (RV) precision of 2.5 cm/sec, 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 for the formation, migration and evolution of the planet system. △ Less

Submitted 14 January, 2021; v1 submitted 2 October, 2020; originally announced October 2020.

Comments: Final version to be published in the Planetary Sciences Journal. 56 pages, 30 figures. Data from the paper and a complete table of forecast JWST times may be found at https://github.com/ericagol/TRAPPIST1_Spitzer/

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

Authors: Veselin B. Kostov, William F. Welsh, Nader Haghighipour, Eric Agol, Daniel C. Fabrycky, Billy Quarles, Gongjie Li, Sean M. Mills, Laurance R. Doyle, Tsevi Mazeh, Jerome A. Orosz, David Martin, Brian Powell

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

Submitted 6 August, 2020; originally announced August 2020.

Comments: 16 pages, 6 figures, 2 tables, AJ accepted

Journal ref: AJ 160 174 (2020)

A Fast, 2D Gaussian Process Method Based on Celerite: Applications to Transiting Exoplanet Discovery and Characterization

Authors: Tyler Gordon, Eric Agol, Daniel Foreman-Mackey

Abstract: Gaussian processes (GPs) are commonly used as a model of stochastic variability in astrophysical time series. In particular, GPs are frequently employed to account for correlated stellar variability in planetary transit light curves. The efficient application of GPs to light curves containing thousands to tens of thousands of datapoints has been made possible by recent advances in GP methods, incl… ▽ More Gaussian processes (GPs) are commonly used as a model of stochastic variability in astrophysical time series. In particular, GPs are frequently employed to account for correlated stellar variability in planetary transit light curves. The efficient application of GPs to light curves containing thousands to tens of thousands of datapoints has been made possible by recent advances in GP methods, including the celerite method. Here we present an extension of the celerite method to two input dimensions, where, typically, the second dimension is small. This method scales linearly with the total number of datapoints when the noise in each large dimension is proportional to the same celerite kernel and only the amplitude of the correlated noise varies in the second dimension. We demonstrate the application of this method to the problem of measuring precise transit parameters from multiwavelength light curves and show that it has the potential to improve transit parameters measurements by orders of magnitude. Applications of this method include transit spectroscopy and exomoon detection, as well a broader set of astronomical problems. △ Less

Submitted 22 July, 2020; v1 submitted 11 July, 2020; originally announced July 2020.

Comments: Submitted to AAS Journals Updated 7/22/20 to include links to code

TRAPPIST-1: Global Results of the Spitzer Exploration Science Program {\it Red Worlds}

Authors: Elsa Ducrot, M. Gillon, L. Delrez, E. Agol, P. Rimmer, M. Turbet, M. N. Günther, B-O. Demory, A. H. M. J. Triaud, E. Bolmont, A. Burgasser, S. J. Carey, J. G. Ingalls, E. Jehin, J. Leconte, S. M. Lederer, D. Queloz, S. N. Raymond, F. Selsis, V. Van Grootel, J. de Wit

Abstract: With more than 1000 hours of observation from Feb 2016 to Oct 2019, the Spitzer Exploration Program Red Worlds (ID: 13067, 13175 and 14223) exclusively targeted TRAPPIST-1, a nearby (12pc) ultracool dwarf star orbited by seven transiting Earth-sized planets, all well-suited for a detailed atmospheric characterization with the upcoming JWST. In this paper, we present the global results of the proje… ▽ More With more than 1000 hours of observation from Feb 2016 to Oct 2019, the Spitzer Exploration Program Red Worlds (ID: 13067, 13175 and 14223) exclusively targeted TRAPPIST-1, a nearby (12pc) ultracool dwarf star orbited by seven transiting Earth-sized planets, all well-suited for a detailed atmospheric characterization with the upcoming JWST. In this paper, we present the global results of the project. We analyzed 88 new transits and combined them with 100 previously analyzed transits, for a total of 188 transits observed at 3.6 or 4.5 $μみゅー$m. We also analyzed 29 occultations (secondary eclipses) of planet b and eight occultations of planet c observed at 4.5 $μみゅー$m to constrain the brightness temperatures of their daysides. We identify several orphan transit-like structures in our Spitzer photometry, but all of them are of low significance. We do not confirm any new transiting planets. We estimate for TRAPPIST-1 transit depth measurements mean noise floors of $\sim$35 and 25 ppm in channels 1 and 2 of Spitzer/IRAC, respectively. most of this noise floor is of instrumental origins and due to the large inter-pixel inhomogeneity of IRAC InSb arrays, and that the much better interpixel homogeneity of JWST instruments should result in noise floors as low as 10ppm, which is low enough to enable the atmospheric characterization of the planets by transit transmission spectroscopy. We construct updated broadband transmission spectra for all seven planets which show consistent transit depths between the two Spitzer channels. We identify and model five distinct high energy flares in the whole dataset, and discuss our results in the context of habitability. Finally, we fail to detect occultation signals of planets b and c at 4.5 $μみゅー$m, and can only set 3$σしぐま$ upper limits on their dayside brightness temperatures (611K for b 586K for c). △ Less

Submitted 24 June, 2020; originally announced June 2020.

Comments: 50 pages, 21 figures. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 640, A112 (2020)

TESS unveils the optical phase curve of KELT-1b. Thermal emission and ellipsoidal variation from the brown dwarf companion, and activity from the star

Authors: C. von Essen, M. Mallonn, A. Piette, N. B. Cowan, N. Madhusudhan, E. Agol, V. Antoci, K. Poppenhaeger, K. G. Stassun, S. Khalafinejad, G. Tautvaišienė

Abstract: We present the detection and analysis of the phase curve of KELT-1b at optical wavelengths, analyzing data taken by the Transiting Exoplanet Survey Satellite (TESS). With a mass of ~27 M_J, KELT-1b is a low-mass brown dwarf. Due to the high mass and close proximity of its companion, the host star has a TESS light curve which shows clear ellipsoidal variations. We model the data with a six-componen… ▽ More We present the detection and analysis of the phase curve of KELT-1b at optical wavelengths, analyzing data taken by the Transiting Exoplanet Survey Satellite (TESS). With a mass of ~27 M_J, KELT-1b is a low-mass brown dwarf. Due to the high mass and close proximity of its companion, the host star has a TESS light curve which shows clear ellipsoidal variations. We model the data with a six-component model: secondary eclipse, phase curve accounting for reflected light and thermal emission, Doppler beaming, ellipsoidal variations, stellar activity and the primary transit. We determine the secondary eclipse depth in the TESS bandpass to be 304 +/- 75 parts-per-million (ppm), the most accurate eclipse depth determined so far for KELT-1b. We measure the amplitude of the phase curve to be 128 +/- 27 ppm, with a corresponding eastward offset between the region of maximum brightness and the substellar point of 19.2 +/- 9.6 degrees, in good agreement with Spitzer measurements. We determine day and night brightness temperatures of 3201 +/- 147 K and 1484 +/- 110 K, respectively, slightly higher than those from Spitzer 3.6 and 4.5 micrometer data. A one-dimensional self-consistent atmospheric model can explain the TESS and Spitzer day side brightness temperatures with thermal emission alone and no reflected light. The night side data can be explained by a model with an internal temperature of ~1100 K, which may be related to the inflated radius. The difference between the TESS and Spitzer brightness temperatures can be explained by stronger molecular opacity in the Spitzer bands. On the night side, this opacity is due primarily to CH4 and CO while on the day side it is due to H2-H2 and H2-He collision-induced absorption. △ Less

Submitted 15 March, 2021; v1 submitted 17 June, 2020; originally announced June 2020.

Comments: 18 pages, 12 figures

The Stellar Variability Noise Floor for Transiting Exoplanet Photometry with PLATO

Authors: Brett M. Morris, Monica G. Bobra, Eric Agol, Yu Jin Lee, Suzanne L. Hawley

Abstract: One of the main science motivations for the ESA PLAnetary Transit and Oscillations (PLATO) mission is to measure exoplanet transit radii with 3% precision. In addition to flares and starspots, stellar oscillations and granulation will enforce fundamental noise floors for transiting exoplanet radius measurements. We simulate light curves of Earth-sized exoplanets transiting continuum intensity imag… ▽ More One of the main science motivations for the ESA PLAnetary Transit and Oscillations (PLATO) mission is to measure exoplanet transit radii with 3% precision. In addition to flares and starspots, stellar oscillations and granulation will enforce fundamental noise floors for transiting exoplanet radius measurements. We simulate light curves of Earth-sized exoplanets transiting continuum intensity images of the Sun taken by the HMI instrument aboard SDO to investigate the uncertainties introduced on the exoplanet radius measurements by stellar granulation and oscillations. After modeling the solar variability with a Gaussian process, we find that the amplitude of solar oscillations and granulation is of order 100 ppm -- similar to the depth of an Earth transit -- and introduces a fractional uncertainty on the depth of transit of 0.73% assuming four transits are observed over the mission duration. However, when we translate the depth measurement into a radius measurement of the planet, we find a much larger radius uncertainty of 3.6%. This is due to a degeneracy between the transit radius ratio, the limb-darkening, and the impact parameter caused by the inability to constrain the transit impact parameter in the presence of stellar variability. We find that surface brightness inhomogeneity due to photospheric granulation contributes a lower limit of only 2 ppm to the photometry in-transit. The radius uncertainty due to granulation and oscillations, combined with the degeneracy with the transit impact parameter, accounts for a significant fraction of the error budget of the PLATO mission, before detector or observational noise is introduced to the light curve. If it is possible to constrain the impact parameter or to obtain follow-up observations at longer wavelengths where limb-darkening is less significant, this may enable higher precision radius measurements. △ Less

Submitted 19 February, 2020; originally announced February 2020.

Comments: 11 pages, 8 figures, accepted in MNRAS

The TRAPPIST-1 JWST Community Initiative

Authors: Michaël Gillon, Victoria Meadows, Eric Agol, Adam J. Burgasser, Drake Deming, René Doyon, Jonathan Fortney, Laura Kreidberg, James Owen, Franck Selsis, Julien de Wit, Jacob Lustig-Yaeger, Benjamin V. Rackham

Abstract: The upcoming launch of the James Webb Space Telescope (JWST) combined with the unique features of the TRAPPIST-1 planetary system should enable the young field of exoplanetology to enter into the realm of temperate Earth-sized worlds. Indeed, the proximity of the system (12pc) and the small size (0.12 Rsun) and luminosity (0.05 Lsun) of its host star should make the comparative atmospheric charact… ▽ More The upcoming launch of the James Webb Space Telescope (JWST) combined with the unique features of the TRAPPIST-1 planetary system should enable the young field of exoplanetology to enter into the realm of temperate Earth-sized worlds. Indeed, the proximity of the system (12pc) and the small size (0.12 Rsun) and luminosity (0.05 Lsun) of its host star should make the comparative atmospheric characterization of its seven transiting planets within reach of an ambitious JWST program. Given the limited lifetime of JWST, the ecliptic location of the star that limits its visibility to 100d per year, the large number of observational time required by this study, and the numerous observational and theoretical challenges awaiting it, its full success will critically depend on a large level of coordination between the involved teams and on the support of a large community. In this context, we present here a community initiative aiming to develop a well-defined sequential structure for the study of the system with JWST and to coordinate on every aspect of its preparation and implementation, both on the observational (e.g. study of the instrumental limitations, data analysis techniques, complementary space-based and ground-based observations) and theoretical levels (e.g. model developments and comparison, retrieval techniques, inferences). Depending on the outcome of the first phase of JWST observations of the planets, this initiative could become the seed of a major JWST Legacy Program devoted to the study of TRAPPIST-1. △ Less

Submitted 13 February, 2020; v1 submitted 12 February, 2020; originally announced February 2020.

Comments: Comments and suggestions are welcome and can be addressed to "t1jwstci_board@u.washington.edu". 23 pages, 8 figures

On the impact of tides on the transit-timing fits to the TRAPPIST-1 system

Authors: Emeline Bolmont, Brice-Olivier Demory, Sergi Blanco-Cuaresma, Eric Agol, Simon L. Grimm, Pierre Auclair-Desrotour, Franck Selsis, Adrien Leleu

Abstract: Transit Timing Variations, or TTVs, can be a very efficient way of constraining masses and eccentricities of multi-planet systems. Recent measurements of the TTVs of TRAPPIST-1 led to an estimate of the masses of the planets, enabling an estimate of their densities. A recent TTV analysis using data obtained in the past two years yields a 34% and 13% increase in mass for TRAPPIST-1b and c, respecti… ▽ More Transit Timing Variations, or TTVs, can be a very efficient way of constraining masses and eccentricities of multi-planet systems. Recent measurements of the TTVs of TRAPPIST-1 led to an estimate of the masses of the planets, enabling an estimate of their densities. A recent TTV analysis using data obtained in the past two years yields a 34% and 13% increase in mass for TRAPPIST-1b and c, respectively. In most studies to date, a Newtonian N-body model is used to fit the masses of the planets, while sometimes general relativity is accounted for. Using the Posidonius N-body code, we show that in the case of the TRAPPIST-1 system, non-Newtonian effects might be also relevant to correctly model the dynamics of the system and the resulting TTVs. In particular, using standard values of the tidal Love number $k_2$ (accounting for the tidal deformation) and the fluid Love number $k_{2f}$ (accounting for the rotational flattening) leads to differences in the TTVs of TRAPPIST-1b and c similar to the differences caused by general relativity. We also show that relaxing the values of tidal Love number $k_2$ and the fluid Love number $k_{2f}$ can lead to TTVs which differ by as much as a few 10~s on a $3-4$-year timescale, which is a potentially observable level. The high values of the Love numbers needed to reach observable levels for the TTVs could be achieved for planets with a liquid ocean, which, if detected, might then be interpreted as a sign that TRAPPIST-1b and TRAPPIST-1c could have a liquid magma ocean. For TRAPPIST-1 and similar systems, the models to fit the TTVs should potentially account for general relativity, for the tidal deformation of the planets, for the rotational deformation of the planets and, to a lesser extent, for the rotational deformation of the star, which would add up to 7x2+1 = 15 additional free parameters in the case of TRAPPIST-1. △ Less

Submitted 5 February, 2020; originally announced February 2020.

Comments: 7 pages, 3 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 635, A117 (2020)

An Automated Method to Detect Transiting Circumbinary Planets

Authors: Diana Windemuth, Eric Agol, Josh Carter, Eric B. Ford, Nader Haghighipour, Jerome A. Orosz, William F. Welsh

Abstract: To date a dozen transiting "Tatooines" or circumbinary planets (CBPs) have been discovered, by eye, in the data from the Kepler mission; by contrast, thousands of confirmed circumstellar planets orbiting around single stars have been detected using automated algorithms. Automated detection of CBPs is challenging because their transits are strongly aperiodic with irregular profiles. Here, we descri… ▽ More To date a dozen transiting "Tatooines" or circumbinary planets (CBPs) have been discovered, by eye, in the data from the Kepler mission; by contrast, thousands of confirmed circumstellar planets orbiting around single stars have been detected using automated algorithms. Automated detection of CBPs is challenging because their transits are strongly aperiodic with irregular profiles. Here, we describe an efficient and automated technique for detecting circumbinary planets that transit their binary hosts in Kepler light curves. Our method accounts for large transit timing and duration variations (TTVs and TDVs), induced by binary reflex motion, in two ways: 1) We directly correct for large-scale TTVs and TDVs in the light curves by using Keplerian models to approximate binary and CBP orbits; and 2) We allow additional aperiodicities on the corrected light curves by employing the Quasi-periodic Automated Transit Search algorithm (QATS). We demonstrate that our method dramatically improves detection significance using simulated data and two previously identified CBP systems, Kepler-35 and Kepler-64. △ Less

Submitted 16 September, 2019; originally announced September 2019.

Comments: 14 pages, 8 figures. Accepted to MNRAS

The Discovery of the Long-Period, Eccentric Planet Kepler-88 d and System Characterization with Radial Velocities and Photodynamical Analysis

Authors: Lauren M. Weiss, Daniel C. Fabrycky, Eric Agol, Sean M. Mills, Andrew W. Howard, Howard Isaacson, Erik A. Petigura, Benjamin J. Fulton, Lea Hirsch, Evan Sinukoff

Abstract: We present the discovery of Kepler-88 d ($P_d = 1403\pm14$ days, $M\mathrm{sin}i_d = 965\pm44\,M_\oplus = 3.04\pm0.14\,M_J$, $e_d = 0.41\pm0.03$) based on six years of radial velocity (RV) follow-up from the W. M. Keck Observatory HIRES spectrograph. Kepler-88 has two previously identified planets. Kepler-88 b (KOI-142.01) transits in the NASA \Kepler\ photometry and has very large transit timing… ▽ More We present the discovery of Kepler-88 d ($P_d = 1403\pm14$ days, $M\mathrm{sin}i_d = 965\pm44\,M_\oplus = 3.04\pm0.14\,M_J$, $e_d = 0.41\pm0.03$) based on six years of radial velocity (RV) follow-up from the W. M. Keck Observatory HIRES spectrograph. Kepler-88 has two previously identified planets. Kepler-88 b (KOI-142.01) transits in the NASA \Kepler\ photometry and has very large transit timing variations. \citet{Nesvorny2013} perfomed a dynamical analysis of the TTVs to uniquely identify the orbital period and mass of the perturbing planet (Kepler-88 c), which was later was confirmed with RVs from the Observatoire de Haute-Provence (OHP, Barros et al. 2014). To fully explore the architecture of this system, we performed photodynamical modeling on the \Kepler\ photometry combined with the RVs from Keck and OHP and stellar parameters from spectroscopy and Gaia. Planet d is not detectable in the photometry, and long-baseline RVs are needed to ascertain its presence. A photodynamical model simultaneously optimized to fit the RVs and \Kepler\ photometry yields the most precise planet masses and orbital properties yet for b and c: $P_b = 10.91647\pm0.00014\,\mathrm{days}$, $M_b=9.5\pm1.2\,M_\oplus$, $P_c=22.2649\pm0.0007\,\mathrm{days}$, and $M_c=214.1\pm5.3\,M_\oplus$. The photodynamical solution also finds that planets b and c have low eccentricites and low mutual inclination, are apsidally anti-aligned, and have conjunctions on the same hemisphere of the star. Continued RV follow-up of systems with small planets will improve our understanding of the link between inner planetary system architectures and giant planets. △ Less

Submitted 29 April, 2020; v1 submitted 5 September, 2019; originally announced September 2019.

Comments: 25 pages, 14 figures, published in AJ

Analytic Planetary Transit Light Curves and Derivatives for Stars with Polynomial Limb Darkening

Authors: Eric Agol, Rodrigo Luger, Daniel Foreman-Mackey

Abstract: We derive analytic, closed-form solutions for the light curve of a planet transiting a star with a limb darkening profile which is a polynomial function of the stellar elevation, up to arbitrary integer order. We provide improved analytic expressions for the uniform, linear, and quadratic limb-darkened cases, as well as novel expressions for higher order integer powers of limb darkening. The formu… ▽ More We derive analytic, closed-form solutions for the light curve of a planet transiting a star with a limb darkening profile which is a polynomial function of the stellar elevation, up to arbitrary integer order. We provide improved analytic expressions for the uniform, linear, and quadratic limb-darkened cases, as well as novel expressions for higher order integer powers of limb darkening. The formulae are crafted to be numerically stable over the expected range of usage. We additionally present analytic formulae for the partial derivatives of instantaneous flux with respect to the radius ratio, impact parameter, and limb darkening coefficients. These expressions are rapid to evaluate, and compare quite favorably in speed and accuracy to existing transit light curve codes. We also use these expressions to numerically compute the first partial derivatives of exposure-time averaged transit light curves with respect to all model parameters. An additional application is modeling eclipsing binary or eclipsing multiple star systems in cases where the stars may be treated as spherically symmetric. We provide code which implements these formulae in C++, Python, IDL, and Julia, with tests and examples of usage. △ Less

Submitted 11 October, 2019; v1 submitted 8 August, 2019; originally announced August 2019.

Comments: 57 pages, 18 figures. Accepted to AJ

Modeling Kepler Eclipsing Binaries: Homogeneous Inference of Orbital & Stellar Properties

Authors: Diana Windemuth, Eric Agol, Aleezah Ali, Flavien Kiefer

Abstract: We report on the properties of eclipsing binaries from the Kepler mission with a newly developed photometric modeling code, which uses the light curve, spectral energy distribution of each binary, and stellar evolution models to infer stellar masses without the need for radial velocity measurements. We present solutions and posteriors to orbital and stellar parameters for 728 systems, forming the… ▽ More We report on the properties of eclipsing binaries from the Kepler mission with a newly developed photometric modeling code, which uses the light curve, spectral energy distribution of each binary, and stellar evolution models to infer stellar masses without the need for radial velocity measurements. We present solutions and posteriors to orbital and stellar parameters for 728 systems, forming the largest homogeneous catalogue of full Kepler binary parameter estimates to date. Using comparisons to published radial velocity measurements, we demonstrate that the inferred properties (e.g., masses) are reliable for well-detached main-sequence binaries, which make up the majority of our sample. The fidelity of our inferred parameters degrades for a subset of systems not well described by input isochrones, such as short-period binaries that have undergone interactions, or binaries with post-main sequence components. Additionally, we identify 35 new systems which show evidence of eclipse timing variations, perhaps from apsidal motion due to binary tides or tertiary companions. We plan to subsequently use these models to search for and constrain the presence of circumbinary planets in Kepler eclipsing binary systems. △ Less

Submitted 31 July, 2019; originally announced August 2019.

Comments: 36 pages, 16 figures; accepted 2019 July 30 to MNRAS

Detection of Hundreds of New Planet Candidates and Eclipsing Binaries in K2 Campaigns 0-8

Authors: Ethan Kruse, Eric Agol, Rodrigo Luger, Daniel Foreman-Mackey

Abstract: We implement a search for exoplanets in campaigns zero through eight (C0-8) of the K2 extension of the Kepler spacecraft. We apply a modified version of the QATS planet search algorithm to K2 light curves produced by the EVEREST pipeline, carrying out the C0-8 search on $1.5 \times 10^5$ target stars with magnitudes in the range of Kp = 9-15. We detect 818 transiting planet candidates, of which 37… ▽ More We implement a search for exoplanets in campaigns zero through eight (C0-8) of the K2 extension of the Kepler spacecraft. We apply a modified version of the QATS planet search algorithm to K2 light curves produced by the EVEREST pipeline, carrying out the C0-8 search on $1.5 \times 10^5$ target stars with magnitudes in the range of Kp = 9-15. We detect 818 transiting planet candidates, of which 374 were undiscovered by prior searches, with {64,15,5,2,1} in {2,3,4,5,6}-planet multi-planet candidate systems, respectively. Of the new planets detected, 100 orbit M dwarfs, including one which is potentially rocky and in the habitable zone. 154 of our candidates reciprocally transit with our Solar System: they are geometrically aligned to see at least one Solar System planet transit. We find candidates which display transit timing variations and dozens of candidates on both period extremes with single transits or ultra-short periods. We point to evidence that our candidates display similar patterns in frequency and size-period relation as confirmed planets, such as tentative evidence for the radius gap. Confirmation of these planet candidates with follow-up studies will increase the number of K2 planets by up to 50%, and characterization of their host stars will improve statistical studies of planet properties. Our sample includes many planets orbiting bright stars amenable for radial velocity follow-up and future characterization with JWST. We also list the 579 eclipsing binary systems detected as part of this search. △ Less

Submitted 13 September, 2019; v1 submitted 24 July, 2019; originally announced July 2019.

Comments: Updated to match published version in ApJS. 37 pages, 32 figures, two very long tables. Machine readable tables available in source

K2-146: Discovery of Planet c, Precise Masses from Transit Timing, and Observed Precession

Authors: Aaron Hamann, Benjamin T. Montet, Daniel C. Fabrycky, Eric Agol, Ethan Kruse

Abstract: K2-146 is a mid-M dwarf ($M_\star = 0.331 \pm 0.009 M_\odot$; $R_\star = 0.330 \pm 0.010 R_\odot$), observed in Campaigns 5, 16, and 18 of the K2 mission. In Campaign 5 data, a single planet was discovered with an orbital period of $2.6$~days and large transit timing variations due to an unknown perturber. Here we analyze data from Campaigns 16 and 18, detecting the transits of a second planet, c,… ▽ More K2-146 is a mid-M dwarf ($M_\star = 0.331 \pm 0.009 M_\odot$; $R_\star = 0.330 \pm 0.010 R_\odot$), observed in Campaigns 5, 16, and 18 of the K2 mission. In Campaign 5 data, a single planet was discovered with an orbital period of $2.6$~days and large transit timing variations due to an unknown perturber. Here we analyze data from Campaigns 16 and 18, detecting the transits of a second planet, c, with an orbital period of $4.0$~days, librating in a 3:2 resonance with planet b. Large, anti-correlated timing variations of both planets exist due to their resonant perturbations. The planets have a mutual inclination of $2.40^\circ\pm0.25^\circ$, which torqued planet c more closely into our line-of-sight. Planet c was grazing in Campaign 5 and thus missed in previous searches; in Campaigns 16 and 18 it is fully transiting, and its transit depth is three times larger. We improve the stellar properties using data from Gaia DR2, and using dynamical fits find that both planets are sub-Neptunes: their masses are $5.77\pm0.18$ and $7.50\pm0.23 M_{\oplus}$ and their radii are $2.04\pm0.06$ and $2.19\pm0.07$ R$_\oplus$, respectively. These mass constraints set the precision record for small exoplanets (a few gas giants have comparable relative precision). These planets lie in the photoevaporation valley when viewed in Radius-Period space, but due to the low-luminosity M-dwarf host star, they lie among the atmosphere-bearing planets when viewed in Radius-Irradiation space. This, along with their densities being 60%-80% that of Earth, suggests that they may both have retained a substantial gaseous envelope. △ Less

Submitted 24 July, 2019; originally announced July 2019.

Comments: 24 pages, 19 figures, accepted for publication in AJ

EXOFASTv2: A public, generalized, publication-quality exoplanet modeling code

Authors: Jason D. Eastman, Joseph E. Rodriguez, Eric Agol, Keivan G. Stassun, Thomas G. Beatty, Andrew Vanderburg, B. Scott Gaudi, Karen A. Collins, Rodrigo Luger

Abstract: We present the next generation public exoplanet fitting software, EXOFASTv2. It is capable of fitting an arbitrary number of planets, radial velocity data sets, astrometric data sets, and/or transits observed with any combination of wavelengths. We model the star simultaneously in the fit and provide several state-of-the-art ways to constrain its properties, including taking advantage of the now-u… ▽ More We present the next generation public exoplanet fitting software, EXOFASTv2. It is capable of fitting an arbitrary number of planets, radial velocity data sets, astrometric data sets, and/or transits observed with any combination of wavelengths. We model the star simultaneously in the fit and provide several state-of-the-art ways to constrain its properties, including taking advantage of the now-ubiquitous all-sky catalog photometry and Gaia parallaxes. EXOFASTv2 can model the star by itself, too. Multi-planet systems are modeled self-consistently with the same underlying stellar mass that defines their semi-major axes through Kepler's law and the planetary period. Transit timing, duration, and depth variations can be modeled with a simple command line option. We explain our methodology and rationale as well as provide an improved version of the core transit model that is both 25\% faster and more accurate. We highlight several potential pitfalls in exoplanet modeling, including the handling of eccentricity in transit-only fits, that the standard exoplanet convention for $ωおめが$ uses a left-handed coordinate system, contrary to most modern textbooks, how to avoid an important degeneracy when allowing negative companion masses, and a widely unappreciated, potential 10-minute ambiguity in the reported transit times. EXOFASTv2 is available at https://github.com/jdeast/EXOFASTv2 . The code is written in IDL, and includes an executable that can be run freely and legally without an IDL license or any knowledge of the language. Extensive documentation and tutorials are included in the distribution for a variety of example fits. Advanced amateurs and undergrads have successfully performed sophisticated global fits of complex planetary systems with EXOFASTv2. It is therefore a powerful tool for education and outreach as well as the broader professional community. △ Less

Submitted 22 July, 2019; originally announced July 2019.

Comments: 51 pages, 14 figures, 3 Tables. Submitted to PASP. Comments welcome

Kepler Object of Interest Network III. Kepler-82f: A new non-transiting $21 M_\bigoplus$ planet from photodynamical modelling

Authors: J. Freudenthal, C. von Essen, A. Ofir, S. ~Dreizler, E. Agol, S. Wedemeyer, B. M. Morris, A. C. Becker, H. J. Deeg, S. Hoyer, M. Mallonn, K. Poppenhaeger, E. Herrero, I. Ribas, P. Boumis, A. Liakos

Abstract: Context. The Kepler Object of Interest Network (KOINet) is a multi-site network of telescopes around the globe organised for follow-up observations of transiting planet candidate Kepler objects of interest (KOIs) with large transit timing variations (TTVs). The main goal of KOINet is the completion of their TTV curves as the Kepler telescope stopped observing the original Kepler field in 2013. A… ▽ More Context. The Kepler Object of Interest Network (KOINet) is a multi-site network of telescopes around the globe organised for follow-up observations of transiting planet candidate Kepler objects of interest (KOIs) with large transit timing variations (TTVs). The main goal of KOINet is the completion of their TTV curves as the Kepler telescope stopped observing the original Kepler field in 2013. Aims. We ensure a comprehensive characterisation of the investigated systems by analysing Kepler data combined with new ground-based transit data using a photodynamical model. This method is applied to the Kepler-82 system leading to its first dynamic analysis. Methods. In order to provide a coherent description of all observations simultaneously, we combine the numerical integration of the gravitational dynamics of a system over the time span of observations with a transit light curve model. To explore the model parameter space, this photodynamical model is coupled with a Markov chain Monte Carlo algorithm. Results. The Kepler-82b/c system shows sinusoidal TTVs due to their near 2:1 resonance dynamical interaction. An additional chopping effect in the TTVs of Kepler-82c hints to a further planet near the 3:2 or 3:1 resonance. We photodynamically analysed Kepler long- and short-cadence data and three new transit observations obtained by KOINet between 2014 and 2018. Our result reveals a non-transiting outer planet with a mass of $m_f=20.9\pm1.0\;M_\bigoplus$ near the 3:2 resonance to the outermost known planet, Kepler-82c. Furthermore, we determined the densities of planets b and c to the significantly more precise values $ρろー_b=0.98_{-0.14}^{+0.10}\;\text{g cm}^{-3}$ and $ρろー_c=0.494_{-0.077}^{+0.066}\;\text{g cm}^{-3}$. △ Less

Submitted 15 July, 2019; originally announced July 2019.

Comments: 18 pages, 9 figures. Accepted for publication in A&A

Journal ref: A&A 628, A108 (2019)

Stellar Properties of Active G and K Stars: Exploring the Connection Between Starspots and Chromospheric Activity

Authors: Brett M. Morris, Jason L. Curtis, Charli Sakari, Suzanne L. Hawley, Eric Agol

Abstract: We gathered high resolution spectra for an ensemble of 55 bright active and inactive stars using the ARC 3.5 m Telescope Echelle Spectrograph at Apache Point Observatory ($R\approx$31,500). We measured spectroscopic effective temperatures, surface gravities and metallicities for most stars in the sample with SME and MOOG. Our stellar property results are consistent with the photometric effective t… ▽ More We gathered high resolution spectra for an ensemble of 55 bright active and inactive stars using the ARC 3.5 m Telescope Echelle Spectrograph at Apache Point Observatory ($R\approx$31,500). We measured spectroscopic effective temperatures, surface gravities and metallicities for most stars in the sample with SME and MOOG. Our stellar property results are consistent with the photometric effective temperatures from the Gaia DR2 pipeline. We also measured their chromospheric $S$ and $\log R^\prime_\mathrm{HK}$ indices to classify the stars as active or inactive and study the connection between chromospheric activity and starspots. We then attempted to infer the starspot covering fractions on the active stars by modeling their spectra as a linear combination of hot and cool inactive stellar spectral templates. We find that it is critical to use precise colors of the stars to place stringent priors on the plausible spot covering fractions. The inferred spot covering fractions generally increase with the chromospheric activity indicator $\log R^\prime_\mathrm{HK}$, though we are largely insensitive to spot coverages smaller than $f_S \lesssim 20$%. We find a dearth of stars with small $\log R^\prime_\mathrm{HK}$ and significant spot coverages. △ Less

Submitted 30 June, 2019; originally announced July 2019.

Comments: Accepted for publication in AJ

APOGEE/Kepler Overlap Yields Orbital Solutions for a Variety of Eclipsing Binaries

Authors: Joni Marie Clark Cunningham, Meredith L. Rawls, Diana Windemuth, Aleezah Ali, Jason Jackiewicz, Eric Agol, Keivan G. Stassun

Abstract: Spectroscopic Eclipsing Binaries (SEBs) are fundamental benchmarks in stellar astrophysics and today are observed in breathtaking detail by missions like TESS, Kepler, and APOGEE. We develop a methodology for simultaneous analysis of high precision Kepler light curves and high resolution near-IR spectra from APOGEE and present orbital solutions and evolutionary histories for a subset of SEBs withi… ▽ More Spectroscopic Eclipsing Binaries (SEBs) are fundamental benchmarks in stellar astrophysics and today are observed in breathtaking detail by missions like TESS, Kepler, and APOGEE. We develop a methodology for simultaneous analysis of high precision Kepler light curves and high resolution near-IR spectra from APOGEE and present orbital solutions and evolutionary histories for a subset of SEBs within this overlap. Radial velocities extracted from APOGEE spectra using the Broadening Function technique are combined with Kepler light curves and to yield binary orbital solutions. The Broadening Function approach yields more precise radial velocities than the standard Cross-Correlation Function, which in turn yields more precise orbital parameters and enables the identification of tertiary stars. The orbital periods of these seven SEBs range from 4 to 40 days. Four of the systems (KIC 5285607, KIC 6864859, KIC 6778289, and KIC 4285087) are well-detached binaries. The remaining three systems have apparent tertiary companions, but each exhibits two eclipses along with at least one spectroscopically varying component (KIC 6449358, KIC 6131659, and KIC 6781535). Gaia distances are available for four targets which we use to estimate temperatures of both members of these SEBs. We explore evolutionary histories in H-R diagram space and estimate ages for this subset of our sample. Finally, we consider the implications for the formation pathways of close binary systems via interactions with tertiary companions. Our methodology combined with the era of big data and observation overlap opens up the possibility of discovering and analyzing large numbers of diverse SEBs, including those with high flux ratios and those in triple systems. △ Less

Submitted 28 June, 2019; originally announced July 2019.

Comments: 25 pages, 23 figures, 4 tables, accepted for publication in AJ

Kepler-62f: Kepler's First Small Planet in the Habitable Zone, but Is It Real?

Authors: William Borucki, Susan E. Thompson, Eric Agol, Christina Hedges

Abstract: Kepler-62f is the first exoplanet small enough to plausibly have a rocky composition orbiting within the habitable zone (HZ) discovered by the Kepler Mission. The planet is 1.4 times the size of the Earth and has an orbital period of 267 days. At the time of its discovery, it had the longest period of any small planet in the habitable zone of a multi-planet system. Because of its long period, only… ▽ More Kepler-62f is the first exoplanet small enough to plausibly have a rocky composition orbiting within the habitable zone (HZ) discovered by the Kepler Mission. The planet is 1.4 times the size of the Earth and has an orbital period of 267 days. At the time of its discovery, it had the longest period of any small planet in the habitable zone of a multi-planet system. Because of its long period, only four transits were observed during Kepler's interval of observations. It was initially missed by the Kepler pipeline, but the first three transits were identified by an independent search by Eric Agol, and it was identified as a planet candidate in subsequent Kepler catalogs. However in the latest catalog of exoplanets (Thompson et al., 2018), it is labeled as a false positive. Recent exoplanet catalogues have evolved from subjective classification to automatic classifications of planet candidates by algorithms (such as `Robovetter'). While exceptionally useful for producing a uniform catalogue, these algorithms sometimes misclassify planet candidates as a false positive, as is the case of Kepler-62f. In particularly valuable cases, i.e., when a small planet has been found orbiting in the habitable zone (HZ), it is important to conduct comprehensive analyses of the data and classification protocols to provide the best estimate of the true status of the detection. In this paper we conduct such analyses and show that Kepler-62f is a true planet and not a false positive. The table of stellar and planet properties has been updated based on GAIA results. △ Less

Submitted 15 May, 2019; v1 submitted 14 May, 2019; originally announced May 2019.

Comments: Published in New Astronomy Reviews special issue on key Kepler discoveries. Published version available here: https://authors.elsevier.com/a/1Z2U34xWxHUnO2

Journal ref: New Astronomy Reviews, Volume 83, November 2018, Pages 28-36

Discovery and Characterization of Kepler-36b

Authors: Eric Agol, Joshua A. Carter

Abstract: We describe the circumstances that led to the discovery of Kepler-36b, and the subsequent characterization of its host planetary system. The Kepler-36 system is remarkable for its physical properties: the close separation of the planets, the contrasting densities of the planets despite their proximity, and the short chaotic timescale. Its discovery and characterization was also remarkable for the… ▽ More We describe the circumstances that led to the discovery of Kepler-36b, and the subsequent characterization of its host planetary system. The Kepler-36 system is remarkable for its physical properties: the close separation of the planets, the contrasting densities of the planets despite their proximity, and the short chaotic timescale. Its discovery and characterization was also remarkable for the novelty of the detection technique and for the precise characterization due to the large transit-timing variations caused by the close proximity of the planets, as well as the precise stellar parameters due to asteroseismology. This was the first multi-planet system whose transit data was processed using a fully consistent photometric-dynamical model, using population Markov Chain Monte Carlo techniques to precisely constrain system parameters. Amongst those parameters, the stellar density was found to be consistent with a complementary, concurrent asteroseismic analysis. In a first, the 3D orientation of the planets was constrained from the lack of transit-duration variations. The system yielded insights into the composition and evolution of short-period planet systems. The denser planet appears to have an Earth-like composition, with uncertainties comparable to the highest precision rocky exoplanet measurements, and the planet densities foreshadowed the rocky/gaseous boundary. The formation of this system remains a mystery, but should yield insights into the migration and evolution of compact exoplanet systems. △ Less

Submitted 13 May, 2019; originally announced May 2019.

Comments: 31 pages, 9 figures. In press with New Astronomy Reviews special issue on key discoveries with Kepler

Journal ref: New Astronomy Reviews, Volume 83, 2018, Pages 18-27, ISSN 1387-6473

The Mass of the White Dwarf Companion in the Self-Lensing Binary KOI-3278: Einstein vs. Newton

Authors: Daniel A. Yahalomi, Yossi Shvartzvald, Eric Agol, Avi Shporer, David W. Latham, Ethan Kruse, John M. Brewer, Lars A. Buchhave, Benjamin J. Fulton, Andrew W. Howard, Howard Isaacson, Erik A. Petigura, Samuel N. Quinn

Abstract: KOI-3278 is a self-lensing stellar binary consisting of a white-dwarf secondary orbiting a Sun-like primary star. Kruse and Agol (2014) noticed small periodic brightenings every 88.18 days in the Kepler photometry and interpreted these as the result of microlensing by a white dwarf with about 63$\%$ of the mass of the Sun. We obtained two sets of spectra for the primary that allowed us to derive t… ▽ More KOI-3278 is a self-lensing stellar binary consisting of a white-dwarf secondary orbiting a Sun-like primary star. Kruse and Agol (2014) noticed small periodic brightenings every 88.18 days in the Kepler photometry and interpreted these as the result of microlensing by a white dwarf with about 63$\%$ of the mass of the Sun. We obtained two sets of spectra for the primary that allowed us to derive three sets of spectroscopic estimates for its effective temperature, surface gravity, and metallicity for the first time. We used these values to update the Kruse and Agol (2014) Einsteinian microlensing model, resulting in a revised mass for the white dwarf of $0.539^{+0.022}_{-0.020} \, M_{\odot}$. The spectra also allowed us to determine radial velocities and derive orbital solutions, with good agreement between the two independent data sets. An independent Newtonian dynamical MCMC model of the combined velocities yielded a mass for the white dwarf of $0.5122^{+0.0057}_{-0.0058} \, M_{\odot}$. The nominal uncertainty for the Newtonian mass is about four times better than for the Einsteinian, $\pm 1.1\%$ vs. $\pm 4.1\%$ and the difference between the two mass determinations is $5.2 \%$. We then present a joint Einsteinian microlensing and Newtonian radial velocity model for KOI-3278, which yielded a mass for the white dwarf of $0.5250^{+0.0082}_{-0.0089} \, M_{\odot}$. This joint model does not rely on any white dwarf evolutionary models or assumptions on the white dwarf mass-radius relation. We discuss the benefits of a joint model of self-lensing binaries, and how future studies of these systems can provide insight into the mass-radius relation of white dwarfs. △ Less

Submitted 4 June, 2019; v1 submitted 24 April, 2019; originally announced April 2019.

Comments: ApJ Accepted; 22 Pages, 8 Figures, 6 Tables and 4 Supplementary Tables

Discovery of a Third Transiting Planet in the Kepler-47 Circumbinary System

Authors: Jerome A. Orosz, William F. Welsh, Nader Haghighipour, Billy Quarles, Donald R. Short, Sean M. Mills, Suman Satyal, Guillermo Torres, Eric Agol, Daniel C. Fabrycky, Daniel Jontof-Hutter, Gur Windmiller, Tobias W. A. Müller, Tobias C. Hinse, William D. Cochran, Michael Endl, Eric B. Ford, Tsevi Mazeh, Jack J. Lissauer

Abstract: Of the nine confirmed transiting circumbinary planet systems, only Kepler-47 is known to contain more than one planet. Kepler-47 b (the "inner planet") has an orbital period of 49.5 days and a radius of about $3\,R_{\oplus}$. Kepler-47 c (the "outer planet") has an orbital period of 303.2 days and a radius of about $4.7\,R_{\oplus}$. Here we report the discovery of a third planet, Kepler-47 d (the… ▽ More Of the nine confirmed transiting circumbinary planet systems, only Kepler-47 is known to contain more than one planet. Kepler-47 b (the "inner planet") has an orbital period of 49.5 days and a radius of about $3\,R_{\oplus}$. Kepler-47 c (the "outer planet") has an orbital period of 303.2 days and a radius of about $4.7\,R_{\oplus}$. Here we report the discovery of a third planet, Kepler-47 d (the "middle planet"), which has an orbital period of 187.4 days and a radius of about $7\,R_{\oplus}$. The presence of the middle planet allows us to place much better constraints on the masses of all three planets, where the $1σしぐま$ ranges are less than $26\,M_{\oplus}$, between $7-43\,M_{\oplus}$, and between $2-5\,M_{\oplus}$ for the inner, middle, and outer planets, respectively. The middle and outer planets have low bulk densities, with $ρろー_{\rm middle} < 0.68$ g cm$^{-3}$ and $ρろー_{\rm outer} < 0.26$ g cm$^{-3}$ at the $1σしぐま$ level. The two outer planets are "tightly packed," assuming the nominal masses, meaning no other planet could stably orbit between them. All of the orbits have low eccentricities and are nearly coplanar, disfavoring violent scattering scenarios and suggesting gentle migration in the protoplanetary disk. △ Less

Submitted 15 April, 2019; originally announced April 2019.

Comments: 68 pages, 30 figures

Journal ref: AJ, 157, 174 (2019)

Near-resonance in a system of sub-Neptunes from TESS

Authors: Samuel N. Quinn, Juliette C. Becker, Joseph E. Rodriguez, Sam Hadden, Chelsea X. Huang, Timothy D. Morton, Fred Adams, David Armstrong, Jason D. Eastman, Jonathan Horner, Stephen R. Kane, Jack J. Lissauer, Joseph D. Twicken, Andrew Vanderburg, Rob Wittenmyer, George R. Ricker, Roland K. Vanderspek, David W. Latham, Sara Seager, Joshua N. Winn, Jon M. Jenkins, Eric Agol, Khalid Barkaoui, Charles A. Beichman, François Bouchy , et al. (56 additional authors not shown)

Abstract: We report the Transiting Exoplanet Survey Satellite ($TESS$) detection of a multi-planet system orbiting the $V=10.9$ K0 dwarf TOI 125. We find evidence for up to five planets, with varying confidence. Three high signal-to-noise transit signals correspond to sub-Neptune-sized planets ($2.76$, $2.79$, and $2.94\ R_{\oplus}$), and we statistically validate the planetary nature of the two inner plane… ▽ More We report the Transiting Exoplanet Survey Satellite ($TESS$) detection of a multi-planet system orbiting the $V=10.9$ K0 dwarf TOI 125. We find evidence for up to five planets, with varying confidence. Three high signal-to-noise transit signals correspond to sub-Neptune-sized planets ($2.76$, $2.79$, and $2.94\ R_{\oplus}$), and we statistically validate the planetary nature of the two inner planets ($P_b = 4.65$ days, $P_c = 9.15$ days). With only two transits observed, we report the outer object ($P_{.03} = 19.98$ days) as a high signal-to-noise ratio planet candidate. We also detect a candidate transiting super-Earth ($1.4\ R_{\oplus}$) with an orbital period of only $12.7$ hours and a candidate Neptune-sized planet ($4.2\ R_{\oplus}$) with a period of $13.28$ days, both at low signal-to-noise. This system is amenable to mass determination via radial velocities and transit timing variations, and provides an opportunity to study planets of similar size while controlling for age and environment. The ratio of orbital periods between TOI 125 b and c ($P_c/P_b = 1.97$) is slightly smaller than an exact 2:1 commensurability and is atypical of multiple planet systems from $Kepler$, which show a preference for period ratios just $wide$ of first-order period ratios. A dynamical analysis refines the allowed parameter space through stability arguments and suggests that, despite the nearly commensurate periods, the system is unlikely to be in resonance. △ Less

Submitted 25 January, 2019; originally announced January 2019.

Comments: Submitted to AAS Journals. 20 pages, 10 figures, 5 tables

The Solar Benchmark: Rotational Modulation of the Sun Reconstructed from Archival Sunspot Records

Authors: Brett M. Morris, James R. A. Davenport, Helen A. C. Giles, Leslie Hebb, Suzanne L. Hawley, Ruth Angus, Peter A. Gilman, Eric Agol

Abstract: We use archival daily spot coverage measurements from Howard et al. (1984) to study the rotational modulation of the Sun as though it were a distant star. A quasi-periodic Gaussian process measures the solar rotation period $P_\mathrm{rot} = 26.3 \pm 0.1$ days, and activity cycle period $P_\mathrm{cyc} = 10.7 \pm 0.3$ years. We attempt to search for evidence of differential rotation in variations… ▽ More We use archival daily spot coverage measurements from Howard et al. (1984) to study the rotational modulation of the Sun as though it were a distant star. A quasi-periodic Gaussian process measures the solar rotation period $P_\mathrm{rot} = 26.3 \pm 0.1$ days, and activity cycle period $P_\mathrm{cyc} = 10.7 \pm 0.3$ years. We attempt to search for evidence of differential rotation in variations of the apparent rotation period throughout the activity cycle and do not detect a clear signal of differential rotation, consistent with the null results of the hare-and-hounds exercise of Aigrain et al. (2015). The full reconstructed solar light curve is available online. △ Less

Submitted 14 January, 2019; originally announced January 2019.

Comments: Accepted for publication in MNRAS

A Higher Cadence Subsurvey Located in the Galactic Plane

Authors: Michael B. Lund, Keivan G. Stassun, Jay Farihi, Eric Agol, Markus Rabus, Avi Shporer, Keaton J. Bell

Abstract: Presently, the Galactic plane receives relatively few observations compared to most of the LSST footprint. While this may address static science, the plane will also represent the highest density of variable Galactic sources. The proper characterization of variability of these sources will benefit greatly from observations at a higher cadence. Presently, the Galactic plane receives relatively few observations compared to most of the LSST footprint. While this may address static science, the plane will also represent the highest density of variable Galactic sources. The proper characterization of variability of these sources will benefit greatly from observations at a higher cadence. △ Less

Submitted 7 December, 2018; originally announced December 2018.

Comments: 12 pages, 2 figures, submitted to Call for White Papers on LSST Cadence Optimization

STARRY: Analytic Occultation Light Curves

Authors: Rodrigo Luger, Eric Agol, Daniel Foreman-Mackey, David P. Fleming, Jacob Lustig-Yaeger, Russell Deitrick

Abstract: We derive analytic, closed form, numerically stable solutions for the total flux received from a spherical planet, moon or star during an occultation if the specific intensity map of the body is expressed as a sum of spherical harmonics. Our expressions are valid to arbitrary degree and may be computed recursively for speed. The formalism we develop here applies to the computation of stellar trans… ▽ More We derive analytic, closed form, numerically stable solutions for the total flux received from a spherical planet, moon or star during an occultation if the specific intensity map of the body is expressed as a sum of spherical harmonics. Our expressions are valid to arbitrary degree and may be computed recursively for speed. The formalism we develop here applies to the computation of stellar transit light curves, planetary secondary eclipse light curves, and planet-planet/planet-moon occultation light curves, as well as thermal (rotational) phase curves. In this paper we also introduce STARRY, an open-source package written in C++ and wrapped in Python that computes these light curves. The algorithm in STARRY is six orders of magnitude faster than direct numerical integration and several orders of magnitude more precise. STARRY also computes analytic derivatives of the light curves with respect to all input parameters for use in gradient-based optimization and inference, such as Hamiltonian Monte Carlo (HMC), allowing users to quickly and efficiently fit observed light curves to infer properties of a celestial body's surface map. △ Less

Submitted 15 October, 2018; originally announced October 2018.

Comments: 55 pages, 20 figures. Accepted to the Astronomical Journal. Check out the code at https://github.com/rodluger/starry

Are Starspots and Plages Co-Located on Active G and K Stars?

Authors: Brett M. Morris, Jason L. Curtis, Stephanie T. Douglas, Suzanne L. Hawley, Marcel A. Agüeros, Monica G. Bobra, Eric Agol

Abstract: We explore the connection between starspots and plages of three main-sequence stars by studying the chromospheric and photospheric activity over several rotation periods. We present simultaneous photometry and high-resolution ($R\sim 31,500$) spectroscopy of KIC 9652680, a young, superflare-producing G1 star with a rotation period of 1.4 days. Its Kepler light curve shows rotational modulation con… ▽ More We explore the connection between starspots and plages of three main-sequence stars by studying the chromospheric and photospheric activity over several rotation periods. We present simultaneous photometry and high-resolution ($R\sim 31,500$) spectroscopy of KIC 9652680, a young, superflare-producing G1 star with a rotation period of 1.4 days. Its Kepler light curve shows rotational modulation consistent with a bright hemisphere followed by a relatively dark hemisphere, generating photometric variability with a semi-amplitude of 4%. We find that KIC 9652680 is darkest when its $S$-index of Ca II H & K emission is at its maximum. We interpret this anti-correlation between flux and $S$ to indicate that dark starspots in the photosphere are co-located with the bright plages in the chromosphere, as they are on the Sun. Moving to lower masses and slower rotators, we present K2 observations with simultaneous spectroscopy of EPIC 211928486 (K5V) and EPIC 211966629 (K4V), two active stars in the 650 Myr-old open cluster Praesepe. The K2 photometry reveals that both stars have rotation periods of 11.7 days; while their flux varies by 1 and 2% respectively, their Ca II H & K $S$-indices seem to hold relatively constant as a function of rotational phase. This suggests that extended chromospheric networks of plages are not concentrated into regions of emission centered on the starspots that drive rotational modulation, unlike KIC 9652680. We also note that the Ca II emission of EPIC 211928486 dipped and recovered suddenly over the duration of one rotation, suggesting that the evolution timescale of plages may be of order the rotation period. △ Less

Submitted 12 September, 2018; originally announced September 2018.

Comments: Accepted for publication in AJ

Non-detection of Contamination by Stellar Activity in the Spitzer Transit Light Curves of TRAPPIST-1

Authors: Brett M. Morris, Eric Agol, Leslie Hebb, Suzanne L. Hawley, Michaël Gillon, Elsa Ducrot, Laetitia Delrez, James Ingalls, Brice-Olivier Demory

Abstract: We apply the transit light curve self-contamination technique of Morris et al. (2018) to search for the effect of stellar activity on the transits of the ultracool dwarf TRAPPIST-1 with 2018 Spitzer photometry. The self-contamination method fits the transit light curves of planets orbiting spotted stars, allowing the host star to be a source of contaminating positive or negative flux which influen… ▽ More We apply the transit light curve self-contamination technique of Morris et al. (2018) to search for the effect of stellar activity on the transits of the ultracool dwarf TRAPPIST-1 with 2018 Spitzer photometry. The self-contamination method fits the transit light curves of planets orbiting spotted stars, allowing the host star to be a source of contaminating positive or negative flux which influences the transit depths but not the ingress/egress durations. We find that none of the planets show statistically significant evidence for self-contamination by bright or dark regions of the stellar photosphere. However, we show that small-scale magnetic activity, analogous in size to the smallest sunspots, could still be lurking in the transit photometry undetected. △ Less

Submitted 8 August, 2018; originally announced August 2018.

Comments: Accepted for publication in ApJL

A second terrestrial planet orbiting the nearby M dwarf LHS 1140

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

Abstract: LHS 1140 is a nearby mid-M dwarf known to host a temperate rocky super-Earth (LHS 1140 b) on a 24.737-day orbit. Based on photometric observations by MEarth and Spitzer as well as Doppler spectroscopy from HARPS, we report the discovery of an additional transiting rocky companion (LHS 1140 c) with a mass of $1.81\pm0.39~{\rm M_{Earth}}$ and a radius of $1.282\pm0.024~{\rm R_{Earth}}$ on a tighter,… ▽ More LHS 1140 is a nearby mid-M dwarf known to host a temperate rocky super-Earth (LHS 1140 b) on a 24.737-day orbit. Based on photometric observations by MEarth and Spitzer as well as Doppler spectroscopy from HARPS, we report the discovery of an additional transiting rocky companion (LHS 1140 c) with a mass of $1.81\pm0.39~{\rm M_{Earth}}$ and a radius of $1.282\pm0.024~{\rm R_{Earth}}$ on a tighter, 3.77795-day orbit. We also obtain more precise estimates of the mass and radius of LHS 1140 b to be $6.98\pm0.89~{\rm M_{Earth}}$ and $1.727\pm0.032~{\rm R_{Earth}}$. The mean densities of planets b and c are $7.5\pm1.0~\rm{g/cm^3}$ and $4.7\pm1.1~\rm{g/cm^3}$, respectively, both consistent with the Earth's ratio of iron to magnesium silicate. The orbital eccentricities of LHS 1140 b and c are consistent with circular orbits and constrained to be below 0.06 and 0.31, respectively, with 90% confidence. Because the orbits of the two planets are co-planar and because we know from previous analyses of Kepler data that compact systems of small planets orbiting M dwarfs are commonplace, a search for more transiting planets in the LHS 1140 system could be fruitful. LHS 1140 c is one of the few known nearby terrestrial planets whose atmosphere could be studied with the upcoming James Webb Space Telescope. △ Less

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

Comments: 22 pages, 8 figures, accepted by AJ

Robust Transiting Exoplanet Radii in the Presence of Starspots from Ingress and Egress Durations

Authors: Brett M. Morris, Eric Agol, Leslie Hebb, Suzanne L. Hawley

Abstract: We typically measure the radii of transiting exoplanets from the transit depth, which is given by the ratio of cross-sectional areas of the planet and star. However, if a star has dark starspots (or bright regions) distributed throughout the transit chord, the transit depth will be biased towards smaller (larger) values, and thus the inferred planet radius will be smaller (larger) if unaccounted f… ▽ More We typically measure the radii of transiting exoplanets from the transit depth, which is given by the ratio of cross-sectional areas of the planet and star. However, if a star has dark starspots (or bright regions) distributed throughout the transit chord, the transit depth will be biased towards smaller (larger) values, and thus the inferred planet radius will be smaller (larger) if unaccounted for. We reparameterize the transit light curve to account for "self-contamination" by photospheric inhomogeneities by splitting the parameter $R_p/R_\star$ into two parameters: one for the radius ratio -- which controls the duration of ingress and egress -- and another which measures the possibly contaminated transit depth. We show that this is equivalent to the formulation for contamination by a second star (with positive or negative flux), and that it is sensitive to time-steady inhomogeneity of the stellar photosphere. We use synthetic light curves of spotted stars at high signal-to-noise to show that the radius recovered from measurement of the ingress/egress duration can recover the true radii of planets transiting spotted stars with axisymmetric spot distributions if the limb-darkening parameters are precisely known. We fit time-averaged high signal-to-noise transit light curves from Kepler and Spitzer of ten planets to measure the planet radii and search for evidence of spot distributions. We find that this sample has a range of measured depths and ingress durations which are self-consistent, providing no strong evidence for contamination by spots. However, there is suggestive evidence for occultation of starspots on Kepler-17, and that relatively bright regions are occulted by the planets of Kepler-412 and HD 80606. Future observations with the James Webb Space Telescope may enable this technique to yield accurate planetary radii in the presence of stellar inhomogeneities. △ Less

Submitted 12 July, 2018; originally announced July 2018.

Comments: Accepted for publication in AJ

The 0.8-4.5$μみゅー$m broadband transmission spectra of TRAPPIST-1 planets

Authors: E. Ducrot, M. Sestovic, B. M. Morris, M. Gillon, A. H. M. J. Triaud, J. de Wit, D. Thimmarayappa, E. Agol, Y. Almleaky, A. Burdanov, A. J. Burgasser, L. Delrez, B-O. Demory, E. Jehin, J. Leconte, J. McCormac, C. Murray, D. Queloz, F. Selsis, S. Thompson, V. Van Grootel

Abstract: The TRAPPIST-1 planetary system represents an exceptional opportunity for the atmospheric characterization of temperate terrestrial exoplanets with the upcoming James Webb Space Telescope (JWST). Assessing the potential impact of stellar contamination on the planets' transit transmission spectra is an essential precursor step to this characterization. Planetary transits themselves can be used to s… ▽ More The TRAPPIST-1 planetary system represents an exceptional opportunity for the atmospheric characterization of temperate terrestrial exoplanets with the upcoming James Webb Space Telescope (JWST). Assessing the potential impact of stellar contamination on the planets' transit transmission spectra is an essential precursor step to this characterization. Planetary transits themselves can be used to scan the stellar photosphere and to constrain its heterogeneity through transit depth variations in time and wavelength. In this context, we present our analysis of 169 transits observed in the optical from space with K2 and from the ground with the SPECULOOS and Liverpool telescopes. Combining our measured transit depths with literature results gathered in the mid/near-IR with Spitzer/IRAC and HST/WFC3, we construct the broadband transmission spectra of the TRAPPIST-1 planets over the 0.8-4.5 $μみゅー$m spectral range. While planets b, d, and f spectra show some structures at the 200-300ppm level, the four others are globally flat. Even if we cannot discard their instrumental origins, two scenarios seem to be favored by the data: a stellar photosphere dominated by a few high-latitude giant (cold) spots, or, alternatively, by a few small and hot (3500-4000K) faculae. In both cases, the stellar contamination of the transit transmission spectra is expected to be less dramatic than predicted in recent papers. Nevertheless, based on our results, stellar contamination can still be of comparable or greater order than planetary atmospheric signals at certain wavelengths. Understanding and correcting the effects of stellar heterogeneity therefore appears essential to prepare the exploration of TRAPPIST-1's with JWST. △ Less

Submitted 22 October, 2018; v1 submitted 3 July, 2018; originally announced July 2018.

Comments: Accepted for publication on AJ on 2018 August 31, published on 2018 October 22

Journal ref: The Astronomical Journal, Volume 156, Number 5, Published 2018 October 22