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 s⁻¹, 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.

Context. The TRAPPIST-1 system hosts seven Earth-sized, temperate exoplanets orbiting an ultra-cool dwarf star. As such, it represents a remarkable setting to study the formation and evolution of terrestrial planets that formed in the same protoplanetary disk. While the sizes of the TRAPPIST-1 planets are all known to better than 5% precision, their densities have significant uncertainties (between 28% and 95%) because of poor constraints on the planet’s masses. Aims. The goal of this paper is to improve our knowledge of the TRAPPIST-1 planetary masses and densities using transit-timing variations (TTVs). The complexity of the TTV inversion problem is known to be particularly acute in multi-planetary systems (convergence issues, degeneracies and size of the parameter space), especially for resonant chain systems such as TRAPPIST-1. Methods. To overcome these challenges, we have used a novel method that employs a genetic algorithm coupled to a full N-body integrator that we applied to a set of 284 individual transit timings. This approach enables us to efficiently explore the parameter space and to derive reliable masses and densities from TTVs for all seven planets. Results. Our new masses result in a five- to eight-fold improvement on the planetary density uncertainties, with precisions ranging from 5% to 12%. These updated values provide new insights into the bulk structure of the TRAPPIST-1 planets. We find that TRAPPIST-1 c and e likely have largely rocky interiors, while planets b, d, f, g, and h require envelopes of volatiles in the form of thick atmospheres, oceans, or ice, in most cases with water mass fractions less than 5%.

Context. KELT-9 b exemplifies a newly emerging class of short-period gaseous exoplanets that tend to orbit hot, early type stars – termed ultra-hot Jupiters. The severe stellar irradiation heats their atmospheres to temperatures of ~4000 K, similar to temperatures of photospheres of dwarf stars. Due to the absence of aerosols and complex molecular chemistry at such temperatures, these planets offer the potential of detailed chemical characterization through transit and day-side spectroscopy. Detailed studies of their chemical inventories may provide crucial constraints on their formation process(es) and evolution history. Aims. We aim to search the optical transmission spectrum of KELT-9 b for absorption lines by metals using the cross-correlation technique. Methods. We analysed two transit observations obtained with the HARPS-N spectrograph. We used an isothermal equilibrium chemistry model to predict the transmission spectrum for each of the neutral and singly ionized atoms with atomic numbers between three and 78. Of these, we identified the elements that are expected to have spectral lines in the visible wavelength range and used those as cross-correlation templates. Results. We detect (>5σ) absorption by Na I, Cr II, Sc II and Y II, and confirm previous detections of Mg I, Fe I, Fe II, and Ti II. In addition, we find evidence of Ca I, Cr I, Co I, and Sr II that will require further observations to verify. The detected absorption lines are significantly deeper than predicted by our model, suggesting that the material is transported to higher altitudes where the density is enhanced compared to a hydrostatic profile, and that the material is part of an extended or outflowing envelope. There appears to be no significant blue-shift of the absorption spectrum due to a net day-to-night side wind. In particular, the strong Fe II feature is shifted by 0.18 ± 0.27 km s⁻¹, consistent with zero. Using the orbital velocity of the planet we derive revised masses and radii of the star and the planet: M_* = 1.978 ± 0.023 M_⊙, R_* = 2.178 ± 0.011 R_⊙, m_p = 2.44 ± 0.70 M_J and R_p = 1.783 ± 0.009 R_J.