TRAPPIST-1 is a fantastic nearby (~39.14 light years) planetary system made of at least seven transiting terrestrial-size, terrestrial-mass planets all receiving a moderate amount of irradiation. To date, this is the most observationally favourable system of potentially habitable planets. Since the announcement of the discovery of TRAPPIST-1 planets in 2016, a growing number of techniques and approaches have been used and proposed to reveal its true nature. Here we have compiled a state-of-the-art overview of all the observational and theoretical constraints that have been obtained so far using these techniques and approaches. The goal is to get a better understanding of whether or not TRAPPIST-1 planets can have atmospheres, and if so, what they are made of. For this, we surveyed the literature on TRAPPIST-1 about topics as broad as irradiation environment, orbital architecture, transit observations, density measurements, stellar contamination, and numerical climate and escape models. Each of these topics adds a brick to our understanding of the likely atmospheres of the seven planets. We show that (i) HST transit observations, (ii) density measurements, (iii) atmospheric escape modelling, and (iv) gas accretion modelling altogether offer solid evidence against the presence of H2-dominated atmospheres around TRAPPIST-1 planets. This means they likely have either (i) a high molecular weight atmosphere or (ii) no atmosphere at all. There are several key challenges ahead to characterize the bulk compositions of the atmospheres (if present) of TRAPPIST-1 planets. The main one so far is characterizing and correcting for the effects of stellar contamination. Fortunately, a new wave of observations with the James Webb Space Telescope and near-infrared high-resolution ground-based spectrographs on very large telescopes will bring significant advances in the coming decade.

Mass-radius relationships for water-rich rocky planets are usually calculated assuming most water is present in condensed (either liquid or solid) form. Planet density estimates are then compared to these mass-radius relationships, even when these planets are more irradiated than the runaway greenhouse irradiation limit (around 1.1~times the insolation at Earth for planets orbiting a Sun-like star), for which water has been shown to be unstable in condensed form and would instead form a thick H2O-dominated atmosphere. Here we use the LMD Generic numerical climate model to derive new mass-radius relationships appropriate for water-rich rocky planets that are more irradiated than the runaway greenhouse irradiation limit, meaning planets endowed with a steam, water-dominated atmosphere. For a given water-to-rock mass ratio, these new mass-radius relationships lead to planet bulk densities much lower than calculated when water is assumed to be in condensed form. In other words, using traditional mass-radius relationships for planets that are more irradiated than the runaway greenhouse irradiation limit tends to dramatically overestimate -- possibly by several orders of magnitude -- their bulk water content. In particular, this result applies to TRAPPIST-1 b, c, and d, which can accommodate a water mass fraction of at most 2, 0.3 and 0.08 %, respectively, assuming planetary core with a terrestrial composition. In addition, we show that significant changes of mass-radius relationships (between planets less and more irradiated than the runaway greenhouse limit) can be used to remove bulk composition degeneracies in multiplanetary systems such as TRAPPIST-1. Finally, we provide an empirical formula for the H2O steam atmosphere thickness which can be used to construct mass-radius relationships for any water-rich, rocky planet more irradiated than the runaway greenhouse irradiation threshold.

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 $\mu$m. We also analyzed 29 occultations (secondary eclipses) of planet b and eight occultations of planet c observed at 4.5 $\mu$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 $\mu$m, and can only set 3$\sigma$ upper limits on their dayside brightness temperatures (611K for b 586K for c).

In an effort to measure the Rossiter-McLaughlin effect for the TRAPPIST-1 system, we performed high-resolution spectroscopy during transits of planets e, f, and b. The spectra were obtained with the InfraRed Doppler spectrograph on the Subaru 8.2-m telescope, and were supplemented with simultaneous photometry obtained with a 1-m telescope of the Las Cumbres Observatory Global Telescope. By analyzing the anomalous radial velocities, we found the projected stellar obliquity to be $\lambda=1\pm 28$ degrees under the assumption that the three planets have coplanar orbits, although we caution that the radial-velocity data show correlated noise of unknown origin. We also sought evidence for the expected deformations of the stellar absorption lines, and thereby detected the "Doppler shadow" of planet b with a false alarm probability of $1.7\,\%$. The joint analysis of the observed residual cross-correlation map including the three transits gave $\lambda=19_{-15}^{+13}$ degrees. These results indicate that the the TRAPPIST-1 star is not strongly misaligned with the common orbital plane of the planets, although further observations are encouraged to verify this conclusion.

The distribution of hot Jupiters, for which star-planet interactions can be significant, questions the evolution of exosystems. We aim to follow the orbital evolution of a planet along the rotational and structural evolution of the host star by taking into account the coupled effects of tidal and magnetic torques from ab initio prescriptions. It allows us to better understand the evolution of star-planet systems and to explain some properties of the distribution of observed close-in planets. To this end we use a numerical model of a coplanar circular star-planet system taking into account stellar structural changes, wind braking and star-planet interactions, called ESPEM (Benbakoura et al. 2019). We find that depending on the initial configuration of the system, magnetic effects can dominate tidal effects during the various phases of the evolution, leading to an important migration of the planet and to significant changes on the rotational evolution of the star. Both kinds of interactions thus have to be taken into account to predict the evolution of compact star-planet systems.

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.

The Transiting Exoplanet Survey Satellite, TESS, is currently carrying out an all-sky search for small planets transiting bright stars. In the first year of the TESS survey, steady progress was made in achieving the mission's primary science goal of establishing bulk densities for 50 planets smaller than Neptune. During that year, TESS's observations were focused on the southern ecliptic hemisphere, resulting in the discovery of three mini-Neptunes orbiting the star TOI-125, a V=11.0 K0 dwarf. We present intensive HARPS radial velocity observations, yielding precise mass measurements for TOI-125b, TOI-125c and TOI-125d. TOI-125b has an orbital period of 4.65 days, a radius of $2.726 \pm 0.075 ~\mathrm{R_{\rm E}}$, a mass of $ 9.50 \pm 0.88 ~\mathrm{M_{\rm E}}$ and is near the 2:1 mean motion resonance with TOI-125c at 9.15 days. TOI-125c has a similar radius of $2.759 \pm 0.10 ~\mathrm{R_{\rm E}}$ and a mass of $ 6.63 \pm 0.99 ~\mathrm{M_{\rm E}}$, being the puffiest of the three planets. TOI-125d, has an orbital period of 19.98 days and a radius of $2.93 \pm 0.17~\mathrm{R_{\rm E}}$ and mass $13.6 \pm 1.2 ~\mathrm{M_{\rm E}}$. For TOI-125b and TOI-125d we find unusual high eccentricities of $0.19\pm 0.04$ and $0.17^{+0.08}_{-0.06}$, respectively. Our analysis also provides upper mass limits for the two low-SNR planet candidates in the system; for TOI-125.04 ($R_P=1.36 ~\mathrm{R_{\rm E}}$, $P=$0.53 days) we find a $2\sigma$ upper mass limit of $1.6~\mathrm{M_{\rm E}}$, whereas TOI-125.05 ( $R_P=4.2^{+2.4}_{-1.4} ~\mathrm{R_{\rm E}}$, $P=$ 13.28 days) is unlikely a viable planet candidate with upper mass limit $2.7~\mathrm{M_{\rm E}}$. We discuss the internal structure of the three confirmed planets, as well as dynamical stability and system architecture for this intriguing exoplanet system.

The dissipation of the kinetic energy of wave-like tidal flows within the convective envelope of low-mass stars is one of the key physical mechanisms that shapes the orbital and rotational dynamics of short-period exoplanetary systems. Although low-mass stars are magnetically active objects, the question of how the star's magnetic field impacts large-scale tidal flows and the excitation, propagation and dissipation of tidal waves still remains open. Our goal is to investigate the impact of stellar magnetism on the forcing of tidal waves, and their propagation and dissipation in the convective envelope of low-mass stars as they evolve. We have estimated the amplitude of the magnetic contribution to the forcing and dissipation of tidally induced magneto-inertial waves throughout the structural and rotational evolution of low-mass stars (from M to F-type). For this purpose, we have used detailed grids of rotating stellar models computed with the stellar evolution code STAREVOL. The amplitude of dynamo-generated magnetic fields is estimated via physical scaling laws at the base and the top of the convective envelope. We find that the large-scale magnetic field of the star has little influence on the excitation of tidal waves in the case of nearly-circular orbits and coplanar hot-Jupiter planetary systems, but that it has a major impact on the way waves are dissipated. Our results therefore indicate that a full magneto-hydrodynamical treatment of the propagation and dissipation of tidal waves is needed to properly assess the impact of star-planet tidal interactions throughout the evolutionary history of low-mass stars hosting short-period massive planets.

Eccentricity tides generate a torque that can drive an ocean planet towards asynchronous rotation states of equilibrium when enhanced by resonances associated with the oceanic tidal modes. We investigate the impact of eccentricity tides on the rotation of rocky planets hosting a thin uniform ocean and orbiting cool dwarf stars such as TRAPPIST-1, with orbital periods ~1-10 days. Combining the linear theory of oceanic tides in the shallow water approximation with the Andrade model for the solid part of the planet, we develop a global model including the coupling effects of ocean loading, self-attraction, and deformation of the solid regions. We derive from this model analytic solutions for the tidal Love numbers and torque exerted on the planet. These solutions are used with realistic values of parameters provided by advanced models of the internal structure and tidal oscillations of solid bodies to explore the parameter space both analytically and numerically. Our model allows us to fully characterise the frequency-resonant tidal response of the planet, and particularly the features of resonances associated with the oceanic tidal modes (eigenfrequencies, resulting maxima of the tidal torque and Love numbers) as functions of the planet parameters (mass, radius, Andrade parameters, ocean depth and Rayleigh drag frequency). Resonances associated with the oceanic tide decrease the critical eccentricity beyond which asynchronous rotation states distinct from the usual spin-orbit resonances can exist. We provide an estimation and scaling laws for this critical eccentricity, which is found to be lowered by roughly one order of magnitude, switching from ~0.3 to ~0.06 in typical cases and to ~0.01 in extremal ones.

Planets similar to Earth - but slightly more irradiated - are expected to enter into a runaway greenhouse state, where all surface water rapidly evaporates, forming an optically thick H2O-dominated atmosphere. For Earth, this extreme climate transition is thought to occur for a ~6% increase only of the solar luminosity, though the exact limit at which the transition would occur is still a highly debated topic. In general, the runaway greenhouse is believed to be a fundamental process in the evolution of Earth-size, temperate planets. Using 1-D radiative-convective climate calculations accounting for thick, hot water vapour-dominated atmospheres, we evaluate the transit atmospheric thickness of a post-runaway greenhouse atmosphere, and find that it could possibly reach over a thousand kilometers (i.e., a few tens of % of Earth radius). This abrupt radius inflation - resulting from the runaway-greenhouse-induced transition - could be detected statistically by ongoing and upcoming space missions such as TESS, CHEOPS and PLATO (combined with precise radial velocity mass measurements with ground-based spectrographs such as ESPRESSO, CARMENES or SPIRou), or even in particular cases in multiplanetary systems such as TRAPPIST-1 (when masses and radii will be known with good enough precision). This result provides the community with an observational test of (1) the concept of runaway greenhouse, that defines the inner edge of the traditional Habitable Zone, and the exact limit of the runaway greenhouse transition. In particular, this could provide an empirical measurement of the irradiation at which Earth analogs transition from a temperate to a runaway greenhouse climate state. This astronomical measurement would make it possible to statistically estimate how close Earth is from the runaway greenhouse. (2) the presence (and statistical abundance) of water in temperate, Earth-size exoplanets.

This is an erratum for the publication Bolmont & Mathis 2016 (Celestial Mechanics and Dynamical Astronomy, 126, 275-296, https://doi.org/10.1007/s10569-016-9690-3). There was a small mistake for the spin integration of our code which we corrected and we take advantage of this erratum to investigate a bit further the influence of a planet on the spin of its host star.

The energy dissipation of wave-like tidal flows in the convective envelope of low-mass stars is one of the key physical mechanisms that shape the orbital and rotational dynamics of short-period planetary systems. Tidal flows, and the excitation, propagation, and dissipation of tidally-induced inertial waves can be modified by stellar magnetic fields (e.g., Wei 2016, 2018, Lin and Ogilvie 2018). It is thus important to assess for which stars, at which location of their internal structure, and at which phase of their evolution, one needs to take into account the effects of magnetic fields on tidal waves. Using scaling laws that provide the amplitude of dynamo-generated magnetic fields along the rotational evolution of these stars (e.g., Christensen et al. 2009, Brun et al. 2015), combined with detailed grids of stellar rotation models (e.g., Amard et al. 2016), we examine the influence of a magnetic field on tidal forcing and dissipation near the tachocline of solar-like stars. We show that full consideration of magnetic fields is required to compute tidal dissipation, but not necessarily for tidal forcing.

TRAPPIST-1 (Gillon et al. 2017) is an extremely compact planetary system: seven earth-sized planets orbit at distances lower than 0.07 AU around one of the smallest M-dwarf known in the close neighborhood of the Sun (with a mass of less than 0.09 $M_\odot$). With 3 planets within the classical habitable zone, this system represents an interesting observational target for future instruments such as the JWST (e.g. Barstow & Irwin 2016). As the planets are close-in, tidal interactions play a crucial role in the evolution of the system by controlling both orbital configurations and rotational states of the planets. For the closest planets, the associated tidal dissipation could have an influence on their internal evolution and potentially on their climate and habitability Turbet et al. (2018). Following (Tobie et al. 2005), we build multilayer models of the internal structure of the TRAPPIST-1 planets accounting for the mass and radius of Grimm et al. (2018), then we compute the tidal response and estimate the tidal heat flux of each planet as well as the profile of tidal heating with depth. Finally, we compare our results to the homogeneous model of Efroimsky (2012) and assess the impact heating rate on the thermal state of each layer of the planet.

The planetary system of TRAPPIST-1, discovered in 2016-2017, is a treasure-trove of information. Thanks to a combination of observational techniques, we have estimates of the radii and masses of the seven planets of this very exotic system. With three planets within the traditional Habitable Zone limits, it is one of the best constrained system of astrobiological interest. I will review here the theoretical constraints we can put on this system by trying to reconstruct its history: its atmospheric evolution which depends on the luminosity evolution of the dwarf star, and its tidal dynamical evolution. These constraints can then be used as hypotheses to assess the habitability of the outer planets of the system with a Global Climate Model.

The surface angular velocity evolution of low-mass stars is now globally understood and the main physical mechanisms involved in it are observationally quite constrained. Additionally, recent observations showed anomalies in the rotation period distribution of open clusters main sequence early K-type stars that cannot be reproduced by current angular momentum evolution model. In this work, we study the parameter space of star-planet system's configurations to investigate if including the tidal star-planet interaction in angular momentum evolution models could reproduce these rotation period distribution's anomalies. To study this effect, we use a parametric angular momentum evolution model that allows for core-envelope decoupling and angular momentum extraction by magnetized stellar wind that we coupled to an orbital evolution code where we take into account the torque due to the tides raised on the star by the planet. We explore different stellar and planetary configurations (stellar mass from 0.5 to 1.0 $\rm M_{\odot}$ and planetary mass from 10 $\rm M_{\oplus}$ to 13 $\rm M_{\rm jup}$) to study their effect on the planetary orbital and stellar rotational evolution. The stellar angular momentum is the most impacted by the star-planet interaction when the planet is engulfed during the early main sequence phase. Thus, if a close-in Jupiter mass planet is initially located around 50\% of the stellar corotation radius, a kink in the rotational period distribution opens around late and early K-type stars during the early main sequence phase. Tidal star-planet interactions can create a kink in the rotation period distribution of low-mass stars, which could possibly account for unexpected scatter seen in the rotational period distribution of young stellar clusters.

Seven temperate Earth-sized exoplanets readily amenable for atmospheric studies transit the nearby ultracool dwarf star TRAPPIST-1 (refs 1,2). Their atmospheric regime is unknown and could range from extended primordial hydrogen-dominated to depleted atmospheres (refs 3-6). Hydrogen in particular is a powerful greenhouse gas that may prevent the habitability of inner planets while enabling the habitability of outer ones (refs 6-8). An atmosphere largely dominated by hydrogen, if cloud-free, should yield prominent spectroscopic signatures in the near-infrared detectable during transits. Observations of the innermost planets have ruled out such signatures (ref 9). However, the outermost planets are more likely to have sustained such a Neptune-like atmosphere (refs 10,11). Here, we report observations for the four planets within or near the system's habitable zone, the circumstellar region where liquid water could exist on a planetary surface (refs 12-14). These planets do not exhibit prominent spectroscopic signatures at near-infrared wavelengths either, which rules out cloud-free hydrogen-dominated atmospheres for TRAPPIST-1 d, e and f, with significance of 8, 6 and 4 sigma, respectively. Such an atmosphere is instead not excluded for planet g. As high-altitude clouds and hazes are not expected in hydrogen-dominated atmospheres around planets with such insolation (refs 15,16), these observations further support their terrestrial and potentially habitable nature.

The recently detected TRAPPIST-1 planetary system, with its seven planets transiting a nearby ultracool dwarf star, offers the first opportunity to perform comparative exoplanetology of temperate Earth-sized worlds. To further advance our understanding of these planets' compositions, energy budgets, and dynamics, we are carrying out an intensive photometric monitoring campaign of their transits with the $\textit{Spitzer Space Telescope}$. In this context, we present 60 new transits of the TRAPPIST-1 planets observed with $\textit{Spitzer}$/IRAC in February and March 2017. We combine these observations with previously published $\textit{Spitzer}$ transit photometry and perform a global analysis of the resulting extensive dataset. This analysis refines the transit parameters and provides revised values for the planets' physical parameters, notably their radii, using updated properties for the star. As part of our study, we also measure precise transit timings that will be used in a companion paper to refine the planets' masses and compositions using the transit timing variations method. TRAPPIST-1 shows a very low level of low-frequency variability in the IRAC 4.5-$\mu$m band, with a photometric RMS of only 0.11$\%$ at a 123-s cadence. We do not detect any evidence of a (quasi-)periodic signal related to stellar rotation. We also analyze the transit light curves individually, to search for possible variations in the transit parameters of each planet due to stellar variability, and find that the $\textit{Spitzer}$ transits of the planets are mostly immune to the effects of stellar variations. These results are encouraging for forthcoming transmission spectroscopy observations of the TRAPPIST-1 planets with the $\textit{James Webb Space Telescope}$.

The very recent discovery of planets orbiting very low mass stars sheds light on these exotic objects. Planetary systems around low-mass stars and brown dwarfs are very different from our solar system: the planets are expected to be much closer than Mercury, in a layout that could resemble the system of Jupiter and its moons. The recent discoveries point in that direction with, for example, the system of Kepler-42 and especially the system of TRAPPIST-1 which has seven planets in a configuration very close to the moons of Jupiter. Low-mass stars and brown dwarfs are thought to be very common in our neighborhood and are thought to host many planetary systems. The planets orbiting in the habitable zone of brown dwarfs (and very low-mass stars) represent one of the next challenges of the following decades: they are the only planets of the habitable zone whose atmosphere we will be able to probe (e.g. with the JWST).

Planetary systems with several planets in compact orbital configurations such as TRAPPIST-1 are surely affected by tidal effects. Its study provides us with important insight about its evolution. We developed a second generation of a N-body code based on the tidal model used in Mercury-T, re-implementing and improving its functionalities using Rust as programming language (including a Python interface for easy use) and the WHFAST integrator. The new open source code ensures memory safety, reproducibility of numerical N-body experiments, it improves the spin integration compared to Mercury-T and allows to take into account a new prescription for the dissipation of tidal inertial waves in the convective envelope of stars. Posidonius is also suitable for binary system simulations with evolving stars.

TRAPPIST-1 planets are invaluable for the study of comparative planetary science outside our Solar System and possibly habitability. First, we derive from N-body simulations possible planetary evolution scenarios, and show that each of the planets are likely to be in synchronous rotation. We then use a 3-D Global Climate Model to explore the possible climates of cool planets of the TRAPPIST-1 system. In particular, we look at the conditions required for cool planets to prevent possible volatile species to be lost by permanent condensation, irreversible burying or photochemical destruction. We also explore the resilience of the same volatiles (when in condensed phase) to a runaway greenhouse process. We find that background atmospheres made of N2, CO or O2 are resistant to atmospheric collapse. However, it should be difficult for TRAPPIST-1 planets to accumulate significant greenhouse gases like CO2, CH4, or NH3. CO2 can easily condense on the nightside, forming glaciers that would flow toward the substellar region. A complete CO2 ice cover is possible on TRAPPIST-1g and h only, although CO2 ice deposits could be gravitationally unstable and get buried beneath the water ice shell in geologically short timescales. Given TRAPPIST-1 planets large EUV irradiation (at least 1000x Titan's flux), CH4 and NH3 should be photodissociated rapidly and thus be hard to accumulate in the atmosphere. Photochemical hazes could then sedimentate and form a surface layer of tholins. Regarding habitability, we confirm that few bars of CO2 would suffice to warm the surface of TRAPPIST-1f and g above the melting point of water. We also show that TRAPPIST-1e is a remarkable candidate for surface habitability. If the planet is today synchronous and abundant in water, then it should always sustain surface liquid water at least in the substellar region, whatever the atmosphere considered.

Planets in close-in orbits interact magnetically and tidally with their host stars. These interactions lead to a net torque that makes close-in planets migrate inward or outward depending on their orbital distance. We compare systematically the strength of magnetic and tidal torques for typical observed star-planet systems (T-Tauri & hot Jupiter, M dwarf & Earth-like planet, K star & hot Jupiter) based on state-of-the-art scaling-laws. We find that depending on the characteristics of the system, tidal or magnetic effects can dominate. For very close-in planets, we find that both torques can make a planet migrate on a timescale as small as 10 to 100 thousands of years. Both effects thus have to be taken into account when predicting the evolution of compact systems.

The ultracool dwarf star TRAPPIST-1 hosts seven Earth-size transiting planets, some of which could harbour liquid water on their surfaces. UV observations are essential to measure their high-energy irradiation, and to search for photodissociated water escaping from their putative atmospheres. Our new observations of TRAPPIST-1 Ly-$\alpha$ line during the transit of TRAPPIST-1c show an evolution of the star emission over three months, preventing us from assessing the presence of an extended hydrogen exosphere. Based on the current knowledge of the stellar irradiation, we investigated the likely history of water loss in the system. Planets b to d might still be in a runaway phase, and planets within the orbit of TRAPPIST-1g could have lost more than 20 Earth oceans after 8 Gyr of hydrodynamic escape. However, TRAPPIST-1e to h might have lost less than 3 Earth oceans if hydrodynamic escape stopped once they entered the habitable zone. We caution that these estimates remain limited by the large uncertainty on the planet masses. They likely represent upper limits on the actual water loss because our assumptions maximize the XUV-driven escape, while photodissociation in the upper atmospheres should be the limiting process. Late-stage outgassing could also have contributed significant amounts of water for the outer, more massive planets after they entered the habitable zone. While our results suggest that the outer planets are the best candidates to search for water with the JWST, they also highlight the need for theoretical studies and complementary observations in all wavelength domains to determine the nature of the TRAPPIST-1 planets, and their potential habitability.

Observations of hot Jupiter type exoplanets suggest that their orbital period distribution depends on the metallicity of their host star. We investigate here whether the impact of the stellar metallicity on the evolution of the tidal dissipation inside the convective envelope of rotating stars and its resulting effect on the planetary migration might be a possible explanation for this observed statistical trend. We use a frequency-averaged tidal dissipation formalism coupled to an orbital evolution code and to rotating stellar evolution models to estimate the effect of a change of stellar metallicity on the evolution of close-in planets. We consider here two different stellar masses: 0.4 and 1.0 $M_{\odot}$ evolving from the early pre-main sequence phase up to the red giant branch. We show that the metallicity of a star has a strong effect on the stellar parameters which in turn strongly influence the tidal dissipation in the convective region. While on the pre-main sequence the dissipation of a metal poor Sun-like star is higher than the dissipation of a metal rich Sun-like star, on the main sequence it is the opposite. However, for the $0.4~M_{\odot}$ star, the dependence of the dissipation with metallicity is much less visible. Using an orbital evolution model, we show that changing the metallicity leads to different orbital evolutions (e.g., planets migrate farther out from an initially fast rotating metal rich star). By using this model, we qualitatively reproduced the observational trends of the population of hot Jupiters with the metallicity of their host stars. However, more steps are needed to improve our model to try to quantitatively fit our results to the observations. Namely, we need to improve the treatment of the rotation evolution in the orbital evolution model and ultimately we need to consistently couple of the orbital model to the stellar evolution model.

Star-planet interactions must be taken into account in stellar models to understand the dynamical evolution of close-in planets. The dependence of the tidal interactions on the structural and rotational evolution of the star is of peculiar importance and should be correctly treated. We quantify how tidal dissipation in the convective envelope of rotating low-mass stars evolves from the pre-main sequence up to the red-giant branch depending on the initial stellar mass. We investigate the consequences of this evolution on planetary orbital evolution. We couple the tidal dissipation formalism described in Mathis (2015) to the stellar evolution code STAREVOL and apply it to rotating stars with masses between 0.3 and 1.4 M$_\odot$. In addition, we generalize the work of Bolmont & Mathis (2016) by following the orbital evolution of close-in planets using the new tidal dissipation predictions for advanced phases of stellar evolution. On the PMS the evolution of tidal dissipation is controlled by the evolution of the internal structure of the contracting star. On the MS it is strongly driven by the variation of surface rotation that is impacted by magnetized stellar winds braking. The main effect of taking into account the rotational evolution of the stars is to lower the tidal dissipation strength by about four orders of magnitude on the main-sequence, compared to a normalized dissipation rate that only takes into account structural changes. The evolution of the dissipation strongly depends on the evolution of the internal structure and rotation of the star. From the pre-main sequence up to the tip of the red-giant branch, it varies by several orders of magnitude, with strong consequences for the orbital evolution of close-in massive planets. These effects are the strongest during the pre-main sequence, implying that the planets are mainly sensitive to the star's early history.

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.

The discovery of Proxima b, a terrestrial temperate planet, presents the opportunity of studying a potentially habitable world in optimal conditions. A key aspect to model its habitability is to understand the radiation environment of the planet in the full spectral domain. We characterize the X-rays to mid-IR radiative properties of Proxima with the goal of providing the top-of-atmosphere fluxes on the planet. We also aim at constraining the fundamental properties of the star. We employ observations from a large number of facilities and make use of different methodologies to piece together the full spectral energy distribution of Proxima. In the high-energy domain, we pay particular attention to the contribution by rotational modulation, activity cycle, and flares so that the data provided are representative of the overall radiation dose received by the atmosphere of the planet. We present the full spectrum of Proxima covering 0.7 to 30000 nm. The integration of the data shows that the top-of-atmosphere average XUV irradiance on Proxima b is 0.293 W m^-2, i.e., nearly 60 times higher than Earth, and that the total irradiance is 877+/-44 W m^-2, or 64+/-3% of the solar constant but with a significantly redder spectrum. We also provide laws for the XUV evolution of Proxima corresponding to two scenarios. Regarding the fundamental properties of Proxima, we find M=0.120+/-0.003 Msun, R=0.146+/-0.007 Rsun, Teff=2980+/-80 K, and L=0.00151+/-0.00008 Lsun. In addition, our analysis reveals a ~20% excess in the 3-30 micron flux of the star that is best interpreted as arising from warm dust in the system. The data provided here should be useful to further investigate the current atmospheric properties of Proxima b as well as its past history, with the overall aim of firmly establishing the habitability of the planet.

One focus of modern astronomy is to detect temperate terrestrial exoplanets well-suited for atmospheric characterisation. A milestone was recently achieved with the detection of three Earth-sized planets transiting (i.e. passing in front of) a star just 8% the mass of the Sun 12 parsecs away. Indeed, the transiting configuration of these planets with the Jupiter-like size of their host star - named TRAPPIST-1 - makes possible in-depth studies of their atmospheric properties with current and future astronomical facilities. Here we report the results of an intensive photometric monitoring campaign of that star from the ground and with the Spitzer Space Telescope. Our observations reveal that at least seven planets with sizes and masses similar to the Earth revolve around TRAPPIST-1. The six inner planets form a near-resonant chain such that their orbital periods (1.51, 2.42, 4.04, 6.06, 9.21, 12.35 days) are near ratios of small integers. This architecture suggests that the planets formed farther from the star and migrated inward. The seven planets have equilibrium temperatures low enough to make possible liquid water on their surfaces.

The TRAPPIST-1 system offers the opportunity to characterize terrestrial, potentially habitable planets orbiting a nearby ultracool dwarf star. We performed a four-orbit reconnaissance with the Space Telescope Imaging Spectrograph onboard the Hubble Space Telescope to study the stellar emission at Lyman-$\alpha$, to assess the presence of hydrogen exospheres around the two inner planets, and to determine their UV irradiation. We detect the Lyman-$\alpha$ line of TRAPPIST-1, making it the coldest exoplanet host star for which this line has been measured. We reconstruct the intrinsic line profile, showing that it lacks broad wings and is much fainter than expected from the stellar X-ray emission. TRAPPIST-1 has a similar X-ray emission as Proxima Cen but a much lower Ly-$\alpha$ emission. This suggests that TRAPPIST-1 chromosphere is only moderately active compared to its transition region and corona. We estimated the atmospheric mass loss rates for all planets, and found that despite a moderate extreme UV emission the total XUV irradiation could be strong enough to strip the atmospheres of the inner planets in a few billions years. We detect marginal flux decreases at the times of TRAPPIST-1b and c transits, which might originate from stellar activity, but could also hint at the presence of extended hydrogen exospheres. Understanding the origin of these Lyman-$\alpha$ variations will be crucial in assessing the atmospheric stability and potential habitability of the TRAPPIST-1 planets.

The astrophysics community uses different tools for computational tasks such as complex systems simulations, radiative transfer calculations or big data. Programming languages like Fortran, C or C++ are commonly present in these tools and, generally, the language choice was made based on the need for performance. However, this comes at a cost: safety. For instance, a common source of error is the access to invalid memory regions, which produces random execution behaviors and affects the scientific interpretation of the results. In 2015, Mozilla Research released the first stable version of a new programming language named Rust. Many features make this new language attractive for the scientific community, it is open source and it guarantees memory safety while offering zero-cost abstraction. We explore the advantages and drawbacks of Rust for astrophysics by re-implementing the fundamental parts of Mercury-T, a Fortran code that simulates the dynamical and tidal evolution of multi-planet systems.

Since 1995, numerous close-in planets have been discovered around low-mass stars (M to A-type stars). These systems are susceptible to be tidally evolving, in particular the dissipation of the kinetic energy of tidal flows in the host star may modify its rotational evolution and also shape the orbital architecture of the surrounding planetary system. Recent theoretical studies have shown that the amplitude of the stellar dissipation can vary over several orders of magnitude as the star evolves, and that it also depends on the stellar mass and rotation. We present here one of the first studies of the dynamics of close-in planets orbiting low-mass stars (from $0.6~M_\odot$ to $1.2~M_\odot$) where we compute the simultaneous evolution of the star's structure, rotation and tidal dissipation in its external convective envelope. We demonstrate that tidal friction due to the stellar dynamical tide, i.e. tidal inertial waves (their restoring force is the Coriolis acceleration) excited in the convection zone, can be larger by several orders of magnitude than the one of the equilibrium tide currently used in celestial mechanics. This is particularly true during the Pre Main Sequence (PMS) phase and to a lesser extent during the Sub Giant (SG) phase. Numerical simulations show that only the high dissipation occurring during the PMS phase has a visible effect on the semi-major axis of close-in planets. We also investigate the effect of the metallicity of the star on the tidal evolution of planets. We find that the higher the metallicity of the star, the higher the dissipation and the larger the tidally-induced migration of the planet.

Ultracool dwarfs (UCD; $T_{\rm eff}<\sim3000~$K) cool to settle on the main sequence after $\sim$1 Gyr. For brown dwarfs, this cooling never stops. Their habitable zone (HZ) thus sweeps inward at least during the first Gyr of their lives. Assuming they possess water, planets found in the HZ of UCDs have experienced a runaway greenhouse phase too hot for liquid water prior to entering the HZ. It has been proposed that such planets are desiccated by this hot early phase and enter the HZ as dry worlds. Here we model the water loss during this pre-HZ hot phase taking into account recent upper limits on the XUV emission of UCDs and using 1D radiation-hydrodynamic simulations. We address the whole range of UCDs but also focus on the planets recently found around the $0.08~M_\odot$ dwarf TRAPPIST-1. Despite assumptions maximizing the FUV-photolysis of water and the XUV-driven escape of hydrogen, we find that planets can retain significant amounts of water in the HZ of UCDs, with a sweet spot in the $0.04$-$0.06~M_\odot$ range. We also studied the TRAPPIST-1 system using observed constraints on the XUV-flux. We find that TRAPPIST-1b and c may have lost as much as 15 Earth Oceans and planet d -- which might be inside the HZ -- may have lost less than 1 Earth Ocean. Depending on their initial water contents, they could have enough water to remain habitable. TRAPPIST-1 planets are key targets for atmospheric characterization and could provide strong constraints on the water erosion around UCDs.

Proxima b is a planet with a minimum mass of 1.3 MEarth orbiting within the habitable zone (HZ) of Proxima Centauri, a very low-mass, active star and the Sun's closest neighbor. Here we investigate a number of factors related to the potential habitability of Proxima b and its ability to maintain liquid water on its surface. We set the stage by estimating the current high-energy irradiance of the planet and show that the planet currently receives 30 times more EUV radiation than Earth and 250 times more X-rays. We compute the time evolution of the star's spectrum, which is essential for modeling the flux received over Proxima b's lifetime. We also show that Proxima b's obliquity is likely null and its spin is either synchronous or in a 3:2 spin-orbit resonance, depending on the planet's eccentricity and level of triaxiality. Next we consider the evolution of Proxima b's water inventory. We use our spectral energy distribution to compute the hydrogen loss from the planet with an improved energy-limited escape formalism. Despite the high level of stellar activity we find that Proxima b is likely to have lost less than an Earth ocean's worth of hydrogen before it reached the HZ 100-200 Myr after its formation. The largest uncertainty in our work is the initial water budget, which is not constrained by planet formation models. We conclude that Proxima b is a viable candidate habitable planet.

Radial velocity monitoring has found the signature of a $M \sin i = 1.3$~M$_\oplus$ planet located within the Habitable Zone (HZ) of Proxima Centauri \citep{Anglada16}. Despite a hotter past and an active host star the planet Proxima~b could have retained enough volatiles to sustain surface habitability \citep{Ribas2016}. Here we use a 3D Global Climate Model (GCM) to simulate Proxima b's atmosphere and water cycle for its two likely rotation modes (1:1 and 3:2 spin-orbit resonances) while varying the unconstrained surface water inventory and atmospheric greenhouse effect. We find that a broad range of atmospheric compositions allow surface liquid water. On a tidally-locked planet with sufficient surface water inventory, liquid water is always present, at least in the substellar region. With a non-synchronous rotation, this requires a minimum greenhouse warming ($\sim$10~mbar of CO$_2$ and 1~bar of N$_2$). If the planet is dryer, $\sim$0.5~bar/1.5~bars of CO$_2$ (respectively for asynchronous/synchronous rotation) suffice to prevent the trapping of any arbitrary small water inventory into polar/nightside ice caps. We produce reflection/emission spectra and phase curves for the simulated climates. We find that atmospheric characterization will be possible by direct imaging with forthcoming large telescopes. The angular separation of $7 \lambda/D$ at 1~$\mu$m (with the E-ELT) and a contrast of $\sim$10$^{-7}$ will enable high-resolution spectroscopy and the search for molecular signatures, including H$_2$O, O$_2$, and CO$_2$. The observation of thermal phase curves can be attempted with JWST, thanks to a contrast of $2\times10^{-5}$ at 10~$\mu$m. Proxima~b will also be an exceptional target for future IR interferometers. Within a decade it will be possible to image Proxima~b and possibly determine whether this exoplanet's surface is habitable.

Contrary to Earth, which has a small orbital eccentricity, some exoplanets discovered in the insolation habitable zone (HZ) have high orbital eccentricities (e.g., up to an eccentricity of $\sim0.97$ for HD~20782~b). This raises the question of the capacity of these planets to host surface liquid water. In order to assess the habitability of an eccentric planet, the mean flux approximation is often used. It states that a planet on an eccentric orbit is called habitable if it receives on average a flux compatible with the presence of surface liquid water. However, as the planets do experience important insolation variations over one orbit and even spend some time outside the HZ for high eccentricities, the question of their habitability might not be as straightforward. We performed a set of simulations using the Global Climate Model LMDz, exploring the limits of the mean flux approximation when varying the luminosity of the host star and the eccentricity of the planet. We computed the climate of tidally locked ocean covered planets with orbital eccentricity from 0 to 0.9 receiving a mean flux equal to Earth's, around stars of luminosity ranging from $L_\odot$ to $10^{-4}~L_\odot$. Using here a definition of habitability based on the presence of surface liquid water, we find that most of the planets considered can sustain surface liquid water on the dayside with an ice cap on the nightside. However, for high eccentricity and high luminosity, planets cannot sustain surface liquid water during the whole orbital period. They completely freeze at apoastron and when approaching periastron an ocean appears around the substellar point. We conclude that the higher the eccentricity and the higher the luminosity of the star, the less reliable the mean flux approximation.

Over the past decade, observations of giant exoplanets (Jupiter-size) have provided key insights into their atmospheres, but the properties of lower-mass exoplanets (sub-Neptune) remain largely unconstrained because of the challenges of observing small planets. Numerous efforts to observe the spectra of super-Earths (exoplanets with masses of one to ten times that of Earth) have so far revealed only featureless spectra. Here we report a longitudinal thermal brightness map of the nearby transiting super-Earth 55 Cancri e revealing highly asymmetric dayside thermal emission and a strong day-night temperature contrast. Dedicated space-based monitoring of the planet in the infrared revealed a modulation of the thermal flux as 55 Cancri e revolves around its star in a tidally locked configuration. These observations reveal a hot spot that is located 41 +- 12 degrees east of the substellar point (the point at which incident light from the star is perpendicular to the surface of the planet). From the orbital phase curve, we also constrain the nightside brightness temperature of the planet to 1380 +- 400 kelvin and the temperature of the warmest hemisphere (centred on the hot spot) to be about 1300 kelvin hotter (2700 +- 270 kelvin) at a wavelength of 4.5 microns, which indicates inefficient heat redistribution from the dayside to the nightside. Our observations are consistent with either an optically thick atmosphere with heat recirculation confined to the planetary dayside, or a planet devoid of atmosphere with low-viscosity magma flows at the surface.

Since twenty years, a large population of close-in planets orbiting various classes of low-mass stars (from M to A-type stars) has been discovered. In such systems, the dissipation of the kinetic energy of tidal flows in the host star may modify its rotational evolution and shape the orbital architecture of the surrounding planetary system. In this context, recent works demonstrated that the amplitude of this dissipation can vary over several orders of magnitude as a function of stellar mass, age and rotation. In addition, the stellar rotation evolution strongly impacts angular momentum exchanges within star-planet systems. Therefore, it is now necessary to take into account the structural and rotational evolution of stars when studying the orbital evolution of close-in planets. The presence of planets may also modify the rotation of the host stars and as a consequence their evolution, magnetic activity and mixing. In this work, we present the first study of the dynamics of close-in planets of various masses orbiting low-mass stars (from $0.6$ to $1.2~M_\star$) where we compute the simultaneous evolution of the star's structure, rotation and tidal dissipation in its external convective enveloppe. We find that tidal friction due to the stellar dynamical tide, i.e. tidal inertial waves excited in the convection zone, can be larger by several orders of magnitude than the one of the equilibrium tide, especially during the PMS phase. Consequently, the orbital migration of the planet is now more pronounced and depends more on the stellar mass, rotation and age. This would very weakly affect the planets in the habitable zone because they are located at orbital distances such that stellar tide-induced migration happens on very long timescales. We also demonstrate that the rotational evolution of host stars is only weakly affected by the presence of planets except for massive companions.

A large proportion of observed planetary systems contain several planets in a compact orbital configuration, and often harbor at least one close-in object. These systems are then most likely tidally evolving. We investigate how the effects of planet-planet interactions influence the tidal evolution of planets. We introduce for that purpose a new open-source addition to the Mercury N-body code, Mercury-T, which takes into account tides, general relativity and the effect of rotation-induced flattening in order to simulate the dynamical and tidal evolution of multi-planet systems. It uses a standard equilibrium tidal model, the constant time lag model. Besides, the evolution of the radius of several host bodies has been implemented (brown dwarfs, M-dwarfs of mass $0.1~M_\odot$, Sun-like stars, Jupiter). We validate the new code by comparing its output for one-planet systems to the secular equations results. We find that this code does respect the conservation of total angular momentum. We applied this new tool to the planetary system Kepler-62. We find that tides influence the stability of the system in some cases. We also show that while the four inner planets of the systems are likely to have slow rotation rates and small obliquities, the fifth planet could have a fast rotation rate and a high obliquity. This means that the two habitable zone planets of this system, Kepler-62e ad f are likely to have very different climate features, and this of course would influence their potential at hosting surface liquid water.

Our knowledge of planets' orbital dynamics, which was based on Solar System studies, has been challenged by the diversity of exoplanetary systems. Around cool and ultra cool dwarfs, the influence of tides on the orbital and spin evolution of planets can strongly affect their climate and their capacity to host surface liquid water. We illustrate the role of tides and dynamics with the extreme case of planets orbiting around brown dwarfs. In multiple planet systems, the eccentricity is excited by planet-planet interactions. Planets are therefore heated up from the inside by the tidally-induced friction. This process can heat a habitable zone planet to such a level that surface liquid water cannot exist. We also talk about the newly discovered potentially habitable Earth-sized planet Kepler-186f. Given the poorly estimated age of the system, the planet could still be evolving towards synchronization and have a high obliquity or be pseudo-synchronized with a zero obliquity. These two configurations would have a different effect on the climate of this planet.

A large number of observed exoplanets are part of multiple planet systems. Most of these systems are sufficiently close-in to be tidally evolving. In such systems, there is a competition between the excitation caused by planet-planet interactions and tidal damping. Using as an example two multiple planet systems, which host planets in the surface liquid water habitable zone (HZ): Kepler-62 and Kepler-186, we show the importance and effect of both planetary and stellar tides on the dynamical evolution of planets and on the climate of the HZ planets.

The Kepler-186 system consists of five planets orbiting an early-M dwarf. The planets have physical radii of 1.0-1.50 R$_\oplus$ and orbital periods of 4 to 130 days. The $1.1~$R$_\oplus$ Kepler-186f with a period of 130 days is of particular interest. Its insolation of roughly $0.32~S_\odot$places it within the liquid water habitable zone. We present a multi-faceted study of the Kepler-186 system. First, we show that the distribution of planet masses can be roughly reproduced if the planets accreted from a high-surface density disk presumably sculpted by an earlier phase of migration. However, our simulations predict the existence of 1-2 undetected planets between planets e and f. Next, we present a dynamical analysis of the system including the effect of tides. The timescale for tidal evolution is short enough that the four inner planets must have small obliquities and near-synchronous rotation rates. Tidal evolution of Kepler-186f is slow enough that its current spin state depends on a combination of its dissipation rate and the stellar age. Finally, we study the habitability of Kepler-186f with a 1-D climate model. The planet's surface temperature can be raised above 273 K with 0.5-5 bars of CO$_2$, depending on the amount of N$_2$ present. Kepler-186f represents a case study of an Earth-sized planet in the cooler regions of the habitable zone of a cool star.

The quest for Earth-like planets represents a major focus of current exoplanet research. While planets that are Earth-sized and smaller have been detected, these planets reside in orbits that are too close to their host star to allow liquid water on their surface. We present the detection of Kepler-186f, a 1.11+\-0.14 Earth radius planet that is the outermost of five planets - all roughly Earth-sized - that transit a 0.47+\-0.05 Rsun star. The intensity and spectrum of the star's radiation places Kepler-186f in the stellar habitable zone, implying that if Kepler-186f has an Earth-like atmosphere and H2O at its surface, then some of this H2O is likely to be in liquid form.

Short period (<50 days) low-mass (<10Mearth) exoplanets are abundant and the few of them whose radius and mass have been measured already reveal a diversity in composition. Some of these exoplanets are found on eccentric orbits and are subjected to strong tides affecting their rotation and resulting in significant tidal heating. Within this population, some planets are likely to be depleted in volatiles and have no atmosphere. We model the thermal emission of these "Super Mercuries" to study the signatures of rotation and tidal dissipation on their infrared light curve. We compute the time-dependent temperature map at the surface and in the subsurface of the planet and the resulting disk-integrated emission spectrum received by a distant observer for any observation geometry. We calculate the illumination of the planetary surface for any Keplerian orbit and rotation. We include the internal tidal heat flow, vertical heat diffusion in the subsurface and generate synthetic light curves. We show that the different rotation periods predicted by tidal models (spin-orbit resonances, pseudo-synchronization) produce different photometric signatures, which are observable provided that the thermal inertia of the surface is high, like that of solid or melted rocks (but not regolith). Tidal dissipation can also directly affect the light curves and make the inference of the rotation more difficult or easier depending on the existence of hot spots on the surface. Infrared light curve measurement with the James Webb Space Telescope and EChO can be used to infer exoplanets' rotation periods and dissipation rates and thus to test tidal models. This data will also constrain the nature of the (sub)surface by constraining the thermal inertia.

Given the very close proximity of their habitable zones, brown dwarfs represent high-value targets in the search for nearby transiting habitable planets that may be suitable for follow-up occultation spectroscopy. In this paper we develop search strategies to find habitable planets transiting brown dwarfs depending on their maximum habitable orbital period (PHZ out). Habitable planets with PHZ out shorter than the useful duration of a night (e.g. 8-10 hrs) can be screened with 100 percent completeness from a single location and in a single night (near-IR). More luminous brown dwarfs require continuous monitoring for longer duration, e.g. from space or from a longitude-distributed network (one test scheduling achieved - 3 telescopes, 13.5 contiguous hours). Using a simulated survey of the 21 closest known brown dwarfs (within 7 pc) we find that the probability of detecting at least one transiting habitable planet is between 4.5 +5.6-1.4 and 56 +31-13 percent, depending on our assumptions. We calculate that brown dwarfs within 5-10 pc are characterizable for potential biosignatures with a 6.5 m space telescope using approx. 1 percent of a 5-year mission's life-time spread over a contiguous segment only one fifth to one tenth of this duration.

Exoplanetary science has reached a historic moment. The James Webb Space Telescope will be capable of probing the atmospheres of rocky planets, and perhaps even search for biologically produced gases. However this is contingent on identifying suitable targets before the end of the mission. A race therefore, is on, to find transiting planets with the most favorable properties, in time for the launch. Here, we describe a realistic opportunity to discover extremely favorable targets - rocky planets transiting nearby brown dwarfs - using the Spitzer Space Telescope as a survey instrument. Harnessing the continuous time coverage and the exquisite precision of Spitzer in a 5,400 hour campaign monitoring nearby brown dwarfs, we will detect a handful of planetary systems with planets as small as Mars. The survey we envision is a logical extension of the immense progress that has been realized in the field of exoplanets and a natural outcome of the exploration of the solar neighborhood to map where the nearest habitable rocky planets are located (as advocated by the 2010 Decadal Survey). Our program represents an essential step towards the atmospheric characterization of terrestrial planets and carries the compelling promise of studying the concept of habitability beyond Earth-like conditions. In addition, our photometric monitoring will provide invaluable observations of a large sample of nearby brown dwarfs situated close to the M/L transition. This is why, we also advocate an immediate public release of the survey data, to guarantee rapid progress on the planet search and provide a treasure trove of data for brown dwarf science.

We use the super Earth 55 Cnc e as a case study to address an observable effect of tidal heating. We investigate whether planet-planet interactions can force the eccentricity of this planet to a level affecting the eclipse depth observed with Spitzer. Using the constant time lag tidal model, we first calculate the observed planet flux as a function of albedo and eccentricity, for different tidal dissipation constants and for two extreme cases: a planet with no heat redistribution and a planet with full heat redistribution. We derive the values of albedo and eccentricity that match the observed transit depth. We then perform N-body simulations of the planetary system including tides and General Relativity to follow the evolution of the eccentricity of planet e. We compare the range of eccentricities given by the simulations with the eccentricities required to alter the eclipse depth. We find that the eccentricity of planet e can be large enough to contribute at a measurable level to the thermal emission measured with Spitzer. This affects the constraints on the albedo of the planet, which can be as high as 0.9 (instead of 0.55 when ignoring tidal heating). We also derive a maximum value for the eccentricity of planet e directly from the eclipse depth: e<0.015 assuming Earth's dissipation constant. Transiting exoplanets in multiple planet systems are exceptional targets for testing tidal models because their tidal luminosity may be observable. Future multi-wavelengths observations of eclipse depth and phase curves (for instance with EChO and JWST) should allow us to better resolve the temperature map of these planets and break the degeneracy between albedo and tidal heating that remains for single band observations. In addition, an accurate determination of the eccentricity will make it possible to constrain the dissipation rate of the planet and to probe its internal structure.

We investigate how the evolution of the stellar spin rate affects, and is affected by, planets in close orbits, via star-planet tidal interactions. To do this, we used a standard equilibrium tidal model to compute the orbital evolution of single planets orbiting both Sun-like stars and 0.1 M\odot M-dwarfs. We tested two stellar spin evolution profiles, one with fast initial rotation (P=1.2 day) and one with slow initial rotation (P=8 day). We tested the effect of varying the stellar and planetary dissipation and the planet's mass and initial orbital radius. Conclusions: Tidal evolution allows to differentiate the early behaviors of extremely close-in planets orbiting either a rapidly rotating star or a slowly rotating star. The early spin-up of the star allows the close-in planets around fast rotators to survive the early evolution. For planets around M-dwarfs, surviving the early evolution means surviving on Gyr timescales whereas for Sun-like stars the spin-down brings about late mergers of Jupiter planets. In light of this study, we can say that differentiating between one spin evolution from another given the present position of planets can be very tricky. Unless we can observe some markers of former evolution it is nearly impossible to distinguish the two very different spin profiles, let alone intermediate spin profiles. Though some conclusions can still be drawn from statistical distributions of planets around fully convective M-dwarfs. However, if the tidal evolution brings about a merger late in its history it can also entail a noticeable acceleration of the star in late ages, so that it is possible to have old stars that spin rapidly. This raises the question of better constraining the age of stars.

The tidal evolution of planets orbiting brown dwarfs (BDs) presents an interesting case study because BDs' terrestrial planet forming region is located extremely close-in. In fact, the habitable zones of BDs range from roughly 0.001 to 0.03 AU and for the lowest-mass BDs are located interior to the Roche limit. In contrast with stars, BDs spin up as they age. Thus, the corotation distance moves inward. This has important implications for the tidal evolution of planets around BDs. We used a standard equilibrium tidal model to compute the orbital evolution of a large ensemble of planet-BD systems. We tested the effect of numerous parameters such as the initial semi-major axis and eccentricity, the rotation period of the BD, the masses of both the BD and planet, and the tidal dissipation factors. We find that all planets that form at or beyond the corotation distance and with initial eccentricities smaller than \sim 0.1 are repelled from the BD. Some planets initially interior to corotation can survive if their inward tidal evolution is slower than the BD's spin evolution, but most initially close-in planets fall onto the BD. We find that the most important parameter for the tidal evolution is the initial orbital distance with respect to the corotation distance. Some planets can survive in the habitable zone for Gyr timescales, although in many cases the habitable zone moves inward past the planet's orbit in just tens to hundreds of Myr. Surviving planets can have orbital periods of less than 10 days (as small as 10 hrs), so they could be observable by transit.