With the discovery of a planetary system around TRAPPIST-1, there has been a surge of interest in ultracool dwarfs as potential planet hosts. Planetary systems around ultracool dwarfs represent our best chance of characterising temperate rocky-planet atmospheres with JWST. However, TRAPPIST-1 remains the only known system of its kind, and the occurrence rate of planets around ultracool dwarfs is still poorly constrained. We seek to perform a complete transit search on the ultracool dwarfs observed by NASA's K2 mission to constrain the occurrence rate of planets around these stars. We filter and characterise the sample of ultracool dwarfs observed by K2, by fitting their spectral energy distributions and using parallaxes from Gaia. We build an automatic pipeline to perform photometry, detrend the lightcurves, and search for transits. Using extensive injection-recovery tests of our pipeline, we compute the detection sensitivity of our search, and thus the completeness of our sample. We infer the planetary occurrence rates within a hierarchical Bayesian model (HBM) to treat uncertain planetary parameters. With the occurrence rate parametrised by a step-wise function, we show a convenient way to directly marginalise over the second level of our HBM (the planet parameters). Our method is applicable generally and can greatly speed up inference with larger catalogues of detected planets. We detect one planet in our sample of 702 ultracool dwarfs: a previously-validated mini-Neptune. We thus infer a mini-Neptune (2 - 4Re) occurrence rate of {\eta} = 0.20+0.16-0.11 within orbital periods of 1 - 20 days. For super-Earths (1 - 2Re) and ice/gas giants (4 - 6Re) within 1 - 20 days, we place 95% credible intervals of {\eta} < 1.14 and {\eta} < 0.29, respectively. If TRAPPIST-1-like systems were ubiquitous, we would have had a 96% chance of finding at least one.

Despite the existence of co-orbital bodies in the solar system, and the prediction of the formation of co-orbital planets by planetary system formation models, no co-orbital exoplanets (also called trojans) have been detected thus far. Here we study the signature of co-orbital exoplanets in transit surveys when two planet candidates in the system orbit the star with similar periods. Such pair of candidates could be discarded as false positives because they are not Hill-stable. However, horseshoe or long libration period tadpole co-orbital configurations can explain such period similarity. This degeneracy can be solved by considering the Transit Timing Variations (TTVs) of each planet. We then focus on the three planet candidates system TOI-178: the two outer candidates of that system have similar orbital period and had an angular separation near $\pi/3$ during the TESS observation of sector 2. Based on the announced orbits, the long-term stability of the system requires the two close-period planets to be co-orbitals. Our independent detrending and transit search recover and slightly favour the three orbits close to a 3:2:2 resonant chain found by the TESS pipeline, although we cannot exclude an alias that would put the system close to a 4:3:2 configuration. We then analyse in more detail the co-orbital scenario. We show that despite the influence of an inner planet just outside the 2:3 mean-motion resonance, this potential co-orbital system can be stable on the Giga-year time-scale for a variety of planetary masses, either on a trojan or a horseshoe orbit. We predict that large TTVs should arise in such configuration with a period of several hundred days. We then show how the mass of each planet can be retrieved from these TTVs.

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 $\mu$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.

As of today, hundreds of hot Jupiters have been found, yet the inflated radii of a large fraction of them remain unexplained. It is still unclear whether a single inflation mechanism is enough to explain the entire distribution of radii, or whether a combination of them is needed. We seek to understand the relationship of hot Jupiter radii with stellar irradiation and mass. We also aim to find the intrinsic physical scatter in their radii caused by unobservable parameters. By constructing a hierarchical Bayesian model, we infer the probabilistic relation between planet radius, mass and incident flux for a sample of 286 gas giants. We separately incorporate the observational uncertainties of the data and the intrinsic physical scatter in the population. This allows us to treat the intrinsic physical scatter in radii (due to latent parameters such as the heavy element fraction) as a parameter to be inferred. We find that the planetary mass plays a key role in the inflation extent, with planets in the range between 0.37 - 0.98 MJ showing the most inflated radii. At higher masses, the radius response to incident flux begins to decrease. Below a threshold of 0.37 +- 0.03 MJ we find that giant exoplanets as a population are unable to maintain inflated radii greater than 1.4 RJ but instead exhibit smaller sizes as the incident flux is increased beyond 10^6 Wm-2. We also find that below 1 MJ, there is a cutoff point at high incident flux beyond which we find no more inflated planets, and that this cutoff point decreases as the mass decreases. At incident fluxes higher than 1.6*10^6 Wm -2 and in a mass range 0.37 - 0.98 MJ, we find no evidence for a population of non-inflated hot Jupiters. Our study demonstrates the importance of population-level analysis to grasp the underlying properties of exoplanets.

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 (TTV). 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%.

The TRAPPIST-1 system is the first transiting planet system found orbiting an ultra-cool dwarf star. At least seven planets similar to Earth in radius and in mass were previously found to transit this host star. Subsequently, TRAPPIST-1 was observed as part of the K2 mission and, with these new data, we report the measurement of an 18.77 d orbital period for the outermost planet, TRAPPIST-1h, which was unconstrained until now. This value matches our theoretical expectations based on Laplace relations and places TRAPPIST-1h as the seventh member of a complex chain, with three-body resonances linking every member. We find that TRAPPIST-1h has a radius of 0.727 Earth radii and an equilibrium temperature of 173 K. We have also measured the rotational period of the star at 3.3 d and detected a number of flares consistent with a low-activity, middle-aged, late M dwarf.