Draft version January 18, 2021 Typeset using LATEX twocolumn style in AASTeX63 Refining the transit timing and photometric analysis of TRAPPIST-1: Masses, radii, densities, dynamics, and ephemerides. Eric Agol ,1 Caroline Dorn ,2 Simon L. Grimm ,3 Martin Turbet ,4 Elsa Ducrot ,5 Laetitia Delrez ,6, 4, 5 Michaël Gillon ,5 Brice-Olivier Demory ,3 Artem Burdanov ,7 Khalid Barkaoui ,8, 5 Zouhair Benkhaldoun ,8 Emeline Bolmont ,4 Adam Burgasser ,9 Sean Carey ,10 Julien de Wit ,7 Daniel Fabrycky ,11 Daniel Foreman-Mackey ,12 Jonas Haldemann ,13 David M. Hernandez ,14 James Ingalls ,10 Emmanuel Jehin ,6 Zachary Langford ,1 Jérémy Leconte ,15 Susan M. Lederer ,16 Rodrigo Luger ,12 Renu Malhotra ,17 Victoria S. Meadows ,1 Brett M. Morris ,3 Francisco J. Pozuelos ,6, 5 Didier Queloz ,18 Sean N. Raymond ,15 Franck Selsis ,15 Marko Sestovic ,3 Amaury H.M.J. Triaud ,19 and Valerie Van Grootel 6 1Astronomy Department and Virtual Planetary Laboratory, University of Washington, Seattle, WA 98195 USA 2University of Zurich, Institute of Computational Sciences, Winterthurerstrasse 190, CH-8057, Zurich, Switzerland 3Center for Space and Habitability, University of Bern, Gesellschaftsstrasse 6, CH-3012, Bern, Switzerland 4Observatoire de Genève, Université de Genève, 51 Chemin des Maillettes, CH-1290 Sauverny, Switzerland 5Astrobiology Research Unit, Université de Liège, Allée du 6 Août 19C, B-4000 Liège, Belgium 6Space Sciences, Technologies and Astrophysics Research (STAR) Institute, Université de Liège, Allée du 6 Août 19C, B-4000 Liège, Belgium 7Department of Earth, Atmospheric and Planetary Science, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA 8Oukaimeden Observatory, High Energy Physics and Astrophysics Laboratory, Cadi Ayyad University, Marrakech, Morocco 9Center for Astrophysics and Space Science, University of California San Diego, La Jolla, CA, 92093, USA 10IPAC, California Institute of Technology, 1200 E California Boulevard, Mail Code 314-6, Pasadena, California 91125, USA 11Department of Astronomy & Astrophysics, University of Chicago, 5640 S. Ellis Ave., Chicago, IL 60637 USA 12Center for Computational Astrophysics, Flatiron Institute, Simons Foundation, 162 5th Ave., New York, NY 10010, USA 13University of Bern, Gesellschaftsstrasse 6, CH-3012, Bern, Switzerland 14Harvard-Smithsonian Center for Astrophysics, 60 Garden St., MS 51, Cambridge, MA 02138, USA 15Laboratoire d’astrophysique de Bordeaux, Univ. Bordeaux, CNRS, B18N, Allée Geoffroy Saint-Hilaire, F-33615 Pessac, France 16NASA Johnson Space Center, 2101 NASA Pkwy., Houston TX, 77058, USA 17Lunar and Planetary Laboratory, The University of Arizona, Tucson, AZ 85721 USA 18Cavendish Laboratory, JJ Thomson Avenue, Cambridge, CB3 0H3, UK 19School of Physics & Astronomy, University of Birmingham, Edgbaston, Birmimgham B15 2TT, UK (Received June 26, 2020; Revised July 10, 2020; Accepted January 18, 2021) Submitted to PSJ ABSTRACT We have collected transit times for the TRAPPIST-1 system with the Spitzer Space Telescope over arXiv:2010.01074v2 [astro-ph.EP] 14 Jan 2021 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 Corresponding author: Eric Agol [email protected] 2 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. Keywords: infrared: planetary systems — planets and satellites: terrestrial planets — planets and satellites: compositions — planets and satellites: fundamental parameters 1. INTRODUCTION although for any given planet there is still a degeneracy The TRAPPIST-1 planetary system took the exo- between a larger core-mass fraction and a volatile enve- planet community by surprise thanks to the high mul- lope (Dorn et al. 2018). In a multi-planet system, the tiplicity of small transiting planets orbiting a very-low- bulk density as a function of planet orbital distance may be used to partly break the compositional degeneracy by mass star (∼ 0:09M ; Gillon et al. 2016, 2017; Luger et al. 2017b; Van Grootel et al. 2018). These unex- assuming a common refractory composition and a water pected nature stems from the fact that this system was composition which increases with orbital distance (Un- found in a survey of only 50 nearby ultracool dwarf stars terborn et al. 2018; Lichtenberg et al. 2019). (Jehin et al. 2011; Gillon et al. 2013), suggesting either The TRAPPIST-1 system was initially found with a a high-frequency of such systems around the latest of ground-based pilot survey using a 60-cm telescope, re- the M-dwarfs (He et al. 2016), or perhaps this discovery vealing two short-period transiting planets, and two ad- was fortuitous (Sestovic & Demory 2020; Sagear et al. ditional orphan transits (Gillon et al. 2016; Burdanov 2020). The proximity of the host star (∼12pc) makes et al. 2018). Subsequent ground-based study of the sys- it brighter in the infrared (K = 10:3) than most ul- tem revealed several additional orphan transits, lead- ing to an incomplete picture of the number of planets tracool dwarfs. Its small size (∼ 0:12R ) means that its planets’ masses and radii are large relative to those and the architecture of the system. A 20-day observa- of the star, which enables precise characterization of the tion campaign with the Spitzer Space Telescope (Werner planets’ properties. The system provides the first oppor- et al. 2004) resolved the confusion, revealing the peri- tunity for a detailed study of potentially rocky, Earth- ods of six of the seven transiting planets (Gillon et al. sized exoplanets with incident fluxes spanning the range 2017), while only a single transit observed of the out- of the terrestrial planets in our Solar System. As such, ermost planet left its orbit in question. A subsequent it has galvanized the exoplanet community to study this observation campaign of the system with the K2 mis- system in detail, both observationally and theoretically, sion included four additional transits of the outer planet, identifying its period, and revealed a series of generalized and has fueled hopes that atmospheric signatures (or 1 even biosignatures) might be detected with the James three-body Laplace relations (GLRs) between adjacent Webb Space Telescope (Barstow & Irwin 2016; Morley triplets of planets (Luger et al. 2017b). Additional ob- et al. 2017; Batalha et al. 2018; Krissansen-Totton et al. servations with Spitzer continued to monitor the transit 2018; Wunderlich et al. 2019; Fauchez et al. 2019; Lustig- times of the seven planets at higher precision than af- Yaeger et al. 2019). forded by ground-based observations. An initial analy- More conservatively, the system provides a poten- sis of the Spitzer data to determine planetary radii and tial laboratory for comparative planetology of terrestrial masses was presented in Delrez et al.(2018a) and Grimm planets, and may provide insight and constraints on the et al.(2018). In total, Spitzer observed TRAPPIST-1 formation and evolution of terrestrial planets around for more than 1075 hrs (nearly 45 days), and the result- the lowest-mass stars. In particular, transiting multi- ing time-series photometry includes 188 transits (Ducrot planet systems afford an opportunity to constrain the et al. 2020). In this paper we complement and extend interior compositions of exoplanets. Sizes from transit depths combined with masses from transit-timing varia- 1 −1 −1 −1 This refers to the condition pP1 − (p + q)P2 + qP3 ≈ 0, tions yield the densities of the planets (e.g. Agol & Fab- which is a generalization of the Laplace resonance with p = 1 rycky 2017). In the case of rocky planets with a thin at- and q = 2 (Papaloizou 2014). mosphere, the bulk density can constrain the core-mass fraction and/or Mg/Fe mass-ratio (Valencia et al. 2007), T1 masses, radii, densities 3 the analysis of Ducrot et al.(2020) to include a transit- dynamical state of all of the planets needs to be solved timing and photodynamic analysis of the system. together and their masses and radii are measured rela- Although the planets in the TRAPPIST-1 system tive to the star, so model comparisons with individual have short orbital periods, ranging from 1.5 to 19 days, planets are not independent. the dynamical interactions accumulate gradually with In this paper we revisit a transit-timing and photomet- time, which requires longer-timescale monitoring to ac- ric analysis with the completed Spitzer program using curately constrain the orbital model. The GLRs also the more extensive transit dataset we now have in hand. cause adjacent pairs of planets to reside near mean- The goal of this program is to provide a more precise −1 −1 motion resonances, for which jPi ≈ (j + k)Pi+1 for understanding of the masses, radii, and densities of the integers j and k for the ith and i + 1th planets.