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.

We present the analysis of the hot-Jupiter KELT-7 b using transmission and emission spectroscopy from the Hubble Space Telescope (HST), both taken with the Wide Field Camera 3 (WFC3). Our study uncovers a rich transmission spectrum which is consistent with a cloud-free atmosphere and suggests the presence of H2O and H-. In contrast, the extracted emission spectrum does not contain strong absorption features and, although it is not consistent with a simple blackbody, it can be explained by a varying temperature-pressure profile, collision induced absorption (CIA) and H-. KELT-7 b had also been studied with other space-based instruments and we explore the effects of introducing these additional datasets. Further observations with Hubble, or the next generation of space-based telescopes, are needed to allow for the optical opacity source in transmission to be confirmed and for molecular features to be disentangled in emission.

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).

This paper presents the atmospheric characterisation of three large, gaseous planets: WASP-127 b, WASP-79 b and WASP-62 b. We analysed spectroscopic data obtained with the G141 grism (1.088 - 1.68 $\mu$m) of the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST) using the Iraclis pipeline and the TauREx3 retrieval code, both of which are publicly available. For WASP-127 b, which is the least dense planet discovered so far and is located in the short-period Neptune desert, our retrieval results found strong water absorption corresponding to an abundance of log(H$_2$O) = -2.71$^{+0.78}_{-1.05}$, and absorption compatible with an iron hydride abundance of log(FeH)=$-5.25^{+0.88}_{-1.10}$, with an extended cloudy atmosphere. We also detected water vapour in the atmospheres of WASP-79 b and WASP-62 b, with best-fit models indicating the presence of iron hydride, too. We used the Atmospheric Detectability Index (ADI) as well as Bayesian log evidence to quantify the strength of the detection and compared our results to the hot Jupiter population study by Tsiaras et al. 2018. While all the planets studied here are suitable targets for characterisation with upcoming facilities such as the James Webb Space Telescope (JWST) and Ariel, WASP-127 b is of particular interest due to its low density, and a thorough atmospheric study would develop our understanding of planet formation and migration.

We analyse the transmission and emission spectra of the ultra-hot Jupiter WASP-76 b, observed with the G141 grism of the Hubble Space Telescope's Wide Field Camera 3 (WFC3). We reduce and fit the raw data for each observation using the open-source software Iraclis before performing a fully Bayesian retrieval using the publicly available analysis suite TauRex 3. Previous studies of the WFC3 transmission spectra of WASP-76 b found hints of titanium oxide (TiO) and vanadium oxide (VO) or non-grey clouds. Accounting for a fainter stellar companion to WASP-76, we reanalyse this data and show that removing the effects of this background star changes the slope of the spectrum, resulting in these visible absorbers no longer being detected, eliminating the need for a non-grey cloud model to adequately fit the data but maintaining the strong water feature previously seen. However, our analysis of the emission spectrum suggests the presence of TiO and an atmospheric thermal inversion, along with a significant amount of water. Given the brightness of the host star and the size of the atmospheric features, WASP-76 b is an excellent target for further characterisation with HST, or with future facilities, to better understand the nature of its atmosphere, to confirm the presence of TiO and to search for other optical absorbers.

The distant ice giants of the Solar System, Uranus and Neptune, have only been visited by one space mission, Voyager 2. The current knowledge on their composition remains very limited despite some recent advances. A better characterization of their composition is however essential to constrain their formation and evolution, as a significant fraction of their mass is made of heavy elements, contrary to the gas giants Jupiter and Saturn. An in situ probe like Galileo would provide us with invaluable direct ground-truth composition measurements. However, some of the condensibles will remain out of the grasp of a shallow probe. While additional constraints could be obtained from a complementary orbiter, thermochemistry and diffusion modeling can further help us to increase the science return of an in situ probe.

Most planets currently amenable to transit spectroscopy are close enough to their host star to exhibit a relatively strong day to night temperature gradient. For hot planets, this leads to cause a chemical composition dichotomy between the two hemispheres. In the extreme case of ultra hot jupiters, some species, such as molecular hydrogen and water, are strongly dissociated on the day-side while others, such as carbon monoxide, are not. However, most current retrieval algorithm rely on 1D forward models that are unable to model this effect. We thus investigate how the 3D structure of the atmosphere biases the abundances retrieved using commonly used algorithms. We study the case of Wasp-121b as a prototypical ultra hot Jupiter. We use the simulations of this planet performed with the SPARC/MIT global climate model (GCM) and generate transmission spectra that fully account for the 3D structure of the atmosphere with Pytmopsh3R. These spectra are then analyzed using the \taurex retrieval code. We find that such ultra hot jupiter's transmission spectra exhibit muted H$_2$O features that originate in the night-side where the temperature, hence the scale-height, is smaller than on the day-side. However, the spectral features of molecules present on the day-side are boosted by both its high temperature and low mean molecular weight. As a result, the retrieved parameters are strongly biased compared to the ground truth. In particular the [CO]/[H$_2$O] is overestimated by one to three orders of magnitude. This must be kept in mind when using such retrieval analysis to infer the C/O ratio of a planet's atmosphere. We also discuss whether indicators can allow us to infer the 3D structure of an observed atmosphere. Finally we show that HST/WFC3 transmission data of Wasp-121b are compatible with the day-night thermal and compositional dichotomy predicted by models.

Robustly modeling the inner edge of the habitable zone is essential for determining the most promising potentially habitable exoplanets for atmospheric characterization. Global climate models (GCMs) have become the standard tool for calculating this boundary, but divergent results have emerged among the various GCMs. In this study, we perform an inter-comparison of standard GCMs used in the field on a rapidly rotating planet receiving a G-star spectral energy distribution and on a tidally locked planet receiving an M-star spectral energy distribution. Experiments both with and without clouds are examined. We find relatively small difference (within 8 K) in global-mean surface temperature simulation among the models in the G-star case with clouds. In contrast, the global-mean surface temperature simulation in the M-star case is highly divergent (20-30 K). Moreover, even differences in the simulated surface temperature when clouds are turned off are significant. These differences are caused by differences in cloud simulation and/or radiative transfer, as well as complex interactions between atmospheric dynamics and these two processes. For example we find that an increase in atmospheric absorption of shortwave radiation can lead to higher relative humidity at high altitudes globally and, therefore, a significant decrease in planetary radiation emitted to space. This study emphasizes the importance of basing conclusions about planetary climate on simulations from a variety of GCMs and motivates the eventual comparison of GCM results with terrestrial exoplanet observations to improve their performance.

Context: The anomalously large radii of hot Jupiters has long been a mystery. However, by combining both theoretical arguments and 2D models, a recent study has suggested that the vertical advection of potential temperature leads to an adiabatic temperature profile in the deep atmosphere hotter than the profile obtained with standard 1D models. Aims: In order to confirm the viability of that scenario, we extend this investigation to three dimensional, time-dependent, models. Methods: We use a 3D GCM, DYNAMICO to perform a series of calculations designed to explore the formation and structure of the driving atmospheric circulations, and detail how it responds to changes in both upper and deep atmospheric forcing. Results: In agreement with the previous, 2D, study, we find that a hot adiabat is the natural outcome of the long-term evolution of the deep atmosphere. Integration times of order $1500$ years are needed for that adiabat to emerge from an isothermal atmosphere, explaining why it has not been found in previous hot Jupiter studies. Models initialised from a hotter deep atmosphere tend to evolve faster toward the same final state. We also find that the deep adiabat is stable against low-levels of deep heating and cooling, as long as the Newtonian cooling time-scale is longer than $\sim 3000$ years at $200$ bar. Conclusions: We conclude that the steady-state vertical advection of potential temperature by deep atmospheric circulations constitutes a robust mechanism to explain hot Jupiter inflated radii. We suggest that future studies of hot Jupiters are evolved for a longer time than currently done, and, when possible, include models initialised with a hot deep adiabat. We stress that this mechanism stems from the advection of entropy by irradiation induced mass flows and does not require (finely tuned) dissipative process, in contrast with most previously suggested scenarios.

The Cassini mission unveiled the intense and diverse activity in Saturn's atmosphere: banded jets, waves, vortices, equatorial oscillations. To set the path towards a better understanding of those phenomena, we performed high-resolution multi-annual numerical simulations of Saturn's atmospheric dynamics. We built a new Global Climate Model [GCM] for Saturn, named the Saturn DYNAMICO GCM, by combining a radiative-seasonal model tailored for Saturn to a hydrodynamical solver based on an icosahedral grid suitable for massively-parallel architectures. The impact of numerical dissipation, and the conservation of angular momentum, are examined in the model before a reference simulation employing the Saturn DYNAMICO GCM with a $1/2^{\circ}$ latitude-longitude resolution is considered for analysis. Mid-latitude banded jets showing similarity with observations are reproduced by our model. Those jets are accelerated and maintained by eddy momentum transfers to the mean flow, with the magnitude of momentum fluxes compliant with the observed values. The eddy activity is not regularly distributed with time, but appears as bursts; both barotropic and baroclinic instabilities could play a role in the eddy activity. The steady-state latitude of occurrence of jets is controlled by poleward migration during the spin-up of our model. At the equator, a weakly-superrotating tropospheric jet and vertically-stacked alternating stratospheric jets are obtained in our GCM simulations. The model produces Yanai (Rossby-gravity), Rossby and Kelvin waves at the equator, as well as extratropical Rossby waves, and large-scale vortices in polar regions. Challenges remain to reproduce Saturn's powerful superrotating jet and hexagon-shaped circumpolar jet in the troposphere, and downward-propagating equatorial oscillation in the stratosphere.

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.

Thermal atmospheric tides have a strong impact on the rotation of terrestrial planets. They can lock these planets into an asynchronous rotation state of equilibrium. We aim at characterizing the dependence of the tidal torque resulting from the semidiurnal thermal tide on the tidal frequency, the planet orbital radius, and the atmospheric surface pressure. The tidal torque is computed from full 3D simulations of the atmospheric climate and mean flows using a generic version of the LMDZ general circulation model (GCM) in the case of a nitrogen-dominated atmosphere. Numerical results are discussed with the help of an updated linear analytical framework. Power scaling laws governing the evolution of the torque with the planet orbital radius and surface pressure are derived. The tidal torque exhibits i) a thermal peak in the vicinity of synchronization, ii) a resonant peak associated with the excitation of the Lamb mode in the high frequency range, and iii) well defined frequency slopes outside these resonances. These features are well explained by our linear theory. Whatever the star-planet distance and surface pressure, the torque frequency spectrum -- when rescaled with the relevant power laws -- always presents the same behaviour. This allows us to provide a single and easily usable empirical formula describing the atmospheric tidal torque over the whole parameter space. With such a formula, the effect of the atmospheric tidal torque can be implemented in evolutionary models of the rotational dynamics of a planet in a computationally efficient, and yet relatively accurate way.

Transmission spectroscopy provides us with information on the atmospheric properties at the limb, which is often intuitively assumed to be a narrow annulus around the planet. Consequently, studies have focused on the effect of atmospheric horizontal heterogeneities along the limb. Here we demonstrate that the region probed in transmission -- the limb -- extends significantly toward the day and night sides of the planet. We show that the strong day-night thermal and compositional gradients expected on synchronous exoplanets create sufficient heterogeneities across the limb to result in important systematic effects on the spectrum and bias its interpretation. To quantify these effects, we developed a 3D radiative transfer model able to generate transmission spectra of atmospheres based on 3D atmospheric structures. We first produce synthetic JWST observations from a simulation of GJ 1214b and demonstrate the necessity of a real 3D approach to model data for such precise observatories. Second, we investigate how day-night temperature gradients cause a systematic bias in retrieval analysis performed with 1D forward models. For that purpose we synthesize a large set of forward spectra for prototypical HD 209458 b and GJ 1214 b type planets varying the temperatures of the day and night sides as well as the width of the transition region. We then perform typical retrieval analyses and compare the retrieved parameters to the ground truth of the input model. This study reveals systematic biases on the retrieved temperature (found to be higher than the terminator temperature) and abundances. This is due to the fact that the hotter dayside is more extended vertically and screens the nightside. These biases are difficult to detect as the 1D profiles used in the retrieval are found to provide an excellent match to the observed spectra. This needs to be kept in mind when interpreting real data.

Oceanic tides are a major source of tidal dissipation. They are a key actor for the orbital and rotational evolution of planetary systems, and contribute to the heating of icy satellites hosting a subsurface ocean. Oceanic tides are characterized by a highly frequency-resonant behavior, which is mainly due to the propagation of surface gravity waves in the case of thin oceans, and internal waves when they are deeper. In this work, we derive self-consistent ab initio expressions of the oceanic tidal torque as a function of the key physical parameters of the system (the ocean depth, the Brunt-V\"ais\"al\"a stratification frequency, the rotation rate, the tidal frequency, the Rayleigh friction). These solutions include the coupled mechanisms of internal and surface gravito-inertial waves, which allows us to study the case of planets hosting deep oceans and offer interesting prospects for the coupling between subsurface oceans and ice shells in the case of icy satellites.

Hot Jupiters are submitted to an intense stellar heating. The resulting thermal tides can torque their atmospheres into asynchronous rotation, while these planets are usually assumed to be locked into spin-orbit synchronization with their host star. Particularly, the thermal atmospheric torque can be greatly enhanced by the dynamical component of the tidal response, that is the component associated with the propagation of internal waves. Owing to the involved complex dynamics, semi-analytical approaches are crucial to understand the physical mechanisms that are responsible for the frequency-resonant behavior of thermal tides, and quantify the atmospheric tidal torque. In this work, we revisit the early works by Arras \& Socrates (2010) and present an improved modeling of thermal tides taking into account rotation and radiative cooling. Using this new modeling, we compute analytically the atmospheric tidal response of hot Jupiters and show that resonances associated with low-frequency internal gravity waves are able to drive asynchronous zonal flows in the range 1-30 days.

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.

An accurate estimate of the inner edge of the habitable zone is critical for determining which exoplanets are potentially habitable and for designing future telescopes to observe them. Here, we explore differences in estimating the inner edge among seven one-dimensional (1D) radiative transfer models: two line-by-line codes (SMART and LBLRTM) as well as five band codes (CAM3, CAM4_Wolf, LMDG, SBDART, and AM2) that are currently being used in global climate models. We compare radiative fluxes and spectra in clear-sky conditions around G- and M-stars, with fixed moist adiabatic profiles for surface temperatures from 250 to 360 K. We find that divergences among the models arise mainly from large uncertainties in water vapor absorption in the window region (10 um) and in the region between 0.2 and 1.5 um. Differences in outgoing longwave radiation increase with surface temperature and reach 10-20 Wm^-2; differences in shortwave reach up to 60 Wm^-2, especially at the surface and in the troposphere, and are larger for an M-dwarf spectrum than a solar spectrum. Differences between the two line-by-line models are significant, although smaller than among the band models. Our results imply that the uncertainty in estimating the insolation threshold of the inner edge (the runaway greenhouse limit) due only to clear-sky radiative transfer is ~10% of modern Earth's solar constant (i.e., ~34 Wm^-2 in global mean) among band models and ~3% between the two line-by-line models. These comparisons show that future work is needed focusing on improving water vapor absorption coefficients in both shortwave and longwave, as well as on increasing the resolution of stellar spectra in broadband models.

A large fraction of known rocky exoplanets are expected to have been spun-down to a state of synchronous rotation, including temperate ones. Studies about the atmospheric and surface processes occurring on such planets thus assume that the day/night sides are fixed with respect to the surface over geological timescales. Here we show that this should not be the case for many synchronous exoplanets. This is due to True Polar Wander (TPW), a well known process occurring on Earth and in the Solar System that can reorient a planet without changing the orientation of its rotational angular momentum with respect to an inertial reference frame. As on Earth, convection in the mantle of rocky exoplanets should continuously distort their inertia tensor, causing reorientation. Moreover, we show that this reorientation is made very efficient by the slower rotation rate of synchronous planets. This is due to the weakness of their combined rotational/tidal bulge---the main stabilizing factor limiting TPW. Stabilization by an elastic lithosphere is also shown to be inefficient. We thus expect the axes of smallest and largest moment of inertia to change continuously over time but to remain closely aligned with the star-planet and orbital axes, respectively.

The James Webb Space Telescope (JWST) presents the opportunity to transform our understanding of planets and the origins of life by revealing the atmospheric compositions, structures, and dynamics of transiting exoplanets in unprecedented detail. However, the high-precision, time-series observations required for such investigations have unique technical challenges, and prior experience with other facilities indicates that there will be a steep learning curve when JWST becomes operational. In this paper we describe the science objectives and detailed plans of the Transiting Exoplanet Community Early Release Science (ERS) Program, which is a recently approved program for JWST observations early in Cycle 1. The goal of this project, for which the obtained data will have no exclusive access period, is to accelerate the acquisition and diffusion of technical expertise for transiting exoplanet observations with JWST, while also providing a compelling set of representative datasets that will enable immediate scientific breakthroughs. The Transiting Exoplanet Community ERS Program will exercise the time-series modes of all four JWST instruments that have been identified as the consensus highest priorities, observe the full suite of transiting planet characterization geometries (transits, eclipses, and phase curves), and target planets with host stars that span an illustrative range of brightnesses. The observations in this program were defined through an inclusive and transparent process that had participation from JWST instrument experts and international leaders in transiting exoplanet studies. Community engagement in the project will be centered on a two-phase Data Challenge that culminates with the delivery of planetary spectra, time-series instrument performance reports, and open-source data analysis toolkits in time to inform the agenda for Cycle 2 of the JWST mission.

The study of extrasolar planets and of the Solar System provides complementary pieces of the mosaic represented by the process of planetary formation. Exoplanets are essential to fully grasp the huge diversity of outcomes that planetary formation and the subsequent evolution of the planetary systems can produce. The orbital and basic physical data we currently possess for the bulk of the exoplanetary population, however, do not provide enough information to break the intrinsic degeneracy of their histories, as different evolutionary tracks can result in the same final configurations. The lessons learned from the Solar System indicate us that the solution to this problem lies in the information contained in the composition of planets. The goal of the Atmospheric Remote-Sensing Infrared Exoplanet Large-survey (ARIEL), one of the three candidates as ESA M4 space mission, is to observe a large and diversified population of transiting planets around a range of host star types to collect information on their atmospheric composition. ARIEL will focus on warm and hot planets to take advantage of their well-mixed atmospheres, which should show minimal condensation and sequestration of high-Z materials and thus reveal their bulk composition across all main cosmochemical elements. In this work we will review the most outstanding open questions concerning the way planets form and the mechanisms that contribute to create habitable environments that the compositional information gathered by ARIEL will allow to tackle

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.

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%.

Oceanic tides are a major source of tidal dissipation. They drive the evolution of planetary systems and the rotational dynamics of planets. However, 2D models commonly used for the Earth cannot be applied to extrasolar telluric planets hosting potentially deep oceans because they ignore the three-dimensional effects related to the ocean vertical structure. Our goal is to investigate in a consistant way the importance of the contribution of internal gravity waves in the oceanic tidal response and to propose a modeling allowing to treat a wide range of cases from shallow to deep oceans. A 3D ab initio model is developed to study the dynamics of a global planetary ocean. This model takes into account compressibility, stratification and sphericity terms, which are usually ignored in 2D approaches. An analytic solution is computed and used to study the dependence of the tidal response on the tidal frequency and on the ocean depth and stratification. In the 2D asymptotic limit, we recover the frequency-resonant behaviour due to surface inertial-gravity waves identified by early studies. As the ocean depth and Brunt-V\"ais\"al\"a frequency increase, the contribution of internal gravity waves grows in importance and the tidal response become three-dimensional. In the case of deep oceans, the stable stratification induces resonances that can increase the tidal dissipation rate by several orders of magnitude. It is thus able to affect significantly the evolution time scale of the planetary rotation.

As brown dwarfs and young giant planets cool down, they are known to experience various chemical transitions --- for example from CO rich L-dwarfs to methane rich T-dwarfs. Those chemical transitions are accompanied by spectral transitions whose sharpness cannot be explained by chemistry alone. In a series of articles, Tremblin et al. proposed that some of the yet unexplained features associated to these transitions could be explained by a reduction of the thermal gradient near the photosphere. To explain, in turn, this more isothermal profile, they invoke the presence of an instability analogous to fingering convection -- compositional convection -- triggered by the change in mean molecular weight of the gas due to the chemical transitions mentioned above. In this short note, we use existing arguments to demonstrate that any turbulent transport, if present, would in fact increase the thermal gradient. This misinterpretation comes from the fact that turbulence mixes/homogenizes entropy (potential temperature) instead of temperature. So, while increasing transport, turbulence in an initially stratified atmosphere actually carries energy downward, whether it is due to fingering or any other type of compositional convection. These processes therefore cannot explain the features observed along the aforementioned transitions by reducing the thermal gradient in the atmosphere of substellar objects. Understanding the microphysical and dynamical properties of clouds at these transitions thus probably remains our best way forward.

Thermal tides can torque the atmosphere of hot Jupiters into asynchronous rotation, while these planets are usually assumed to be locked into spin-orbit synchronization with their host star. In this work, our goal is to characterize the tidal response of a rotating hot Jupiter to the tidal semidiurnal thermal forcing of its host star, by identifying the structure of tidal waves responsible for variation of mass distribution, their dependence on the tidal frequency and their ability to generate strong zonal flows. We develop an ab initio global modeling that generalizes the early approach of Arras & Socrates (2010) to rotating and non-adiabatic planets. We derive analytically the torque exerted on the body and the associated timescales of evolution, as well as the equilibrium tidal response of the atmosphere in the zero-frequency limit. Finally, we integrate numerically the equations of thermal tides for three cases including dissipation and rotation step by step. The resonances associated with tidally generated gravito-inertial waves amplify significantly the resulting tidal torque in the range 1-30 days. This torque can drive globally the atmosphere into asynchronous rotation, as its sign depends on the tidal frequency. The resonant behaviour of the tidal response is enhanced by rotation, which couples the forcing to several Hough modes in the general case, while the radiative cooling tends to regularize it and diminish its amplitude.

The ARIEL (Atmospheric Remote-sensing Exoplanet Large-survey) mission concept is one of the three M4 mission candidates selected by the European Space Agency (ESA) for a Phase A study, competing for a launch in 2026. ARIEL has been designed to study the physical and chemical properties of a large and diverse sample of exoplanets, and through those understand how planets form and evolve in our galaxy. Here we describe the assumptions made to estimate an optimal sample of exoplanets - including both the already known exoplanets and the "expected" ones yet to be discovered - observable by ARIEL and define a realistic mission scenario. To achieve the mission objectives, the sample should include gaseous and rocky planets with a range of temperatures around stars of different spectral type and metallicity. The current ARIEL design enables the observation of ~1000 planets, covering a broad range of planetary and stellar parameters, during its four year mission lifetime. This nominal list of planets is expected to evolve over the years depending on the new exoplanet discoveries.

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}$.

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.

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.

It is widely believed that the carbonate-silicate cycle is the main agent to trigger deglaciations by CO$_2$ greenhouse warming on Earth and on Earth-like planets when they get in frozen state. Here we use a 3D Global Climate Model to simulate the ability of frozen planets to escape from glaciation by accumulating enough gaseous CO$_2$. We find that Earth-like planets orbiting a Sun-like star may never be able to escape from glaciation if their orbital distance is greater than $\sim$ 1.27 AU (Flux $<$ 847 W m$^{-2}$ or 62$\%$ of the Solar constant), because CO$_2$ would condense at the poles forming permanent CO$_2$ ice caps. This limits the amount of CO$_2$ in the atmosphere and thus its greenhouse effect. The amount of CO$_2$ that can be trapped in the polar caps depends on the efficiency of CO$_2$ ice to flow laterally as well as its graviational stability relative to subsurface water ice. The flow of CO$_2$ ice from poles to equator is mostly controlled by the bottom temperature, and hence by the internal heat flux. We find that a frozen Earth-like planet located at 1.30 AU of a Sun-like star could store as much as 1.5/4.5/15 bars of dry ice at the poles, for internal heat fluxes of 100/30/10 mW m$^{-2}$. But these amounts are lower limits. For planets with a significant water ice cover, we show that CO$_2$ ice deposits should be gravitationnally unstable. They get buried beneath the water ice cover in geologically short timescales of $\sim$~10$^4$ yrs, mainly controlled by the viscosity of water ice. CO$_2$ would be permanently sequestered underneath the water ice cover, in the form of CO$_2$ liquids, CO$_2$ clathrate hydrates and/or dissolved in subglacial water reservoirs (if any). This would considerably increase the amount of CO$_2$ trapped and further reduce the probability of deglaciation.

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.

Thermochemical models have been used in the past to constrain the deep oxygen abundance in the gas and ice giant planets from tropospheric CO spectroscopic measurements. Knowing the oxygen abundance of these planets is a key to better understand their formation. These models have widely used dry and/or moist adiabats to extrapolate temperatures from the measured values in the upper troposphere down to the level where the thermochemical equilibrium between H$_2$O and CO is established. The mean molecular mass gradient produced by the condensation of H$_2$O stabilizes the atmosphere against convection and results in a vertical thermal profile and H$_2$O distribution that departs significantly from previous estimates. We revisit O/H estimates using an atmospheric structure that accounts for the inhibition of the convection by condensation. We use a thermochemical network and the latest observations of CO in Uranus and Neptune to calculate the internal oxygen enrichment required to satisfy both these new estimates of the thermal profile and the observations. We also present the current limitations of such modeling.

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.

In an atmosphere, a cloud condensation region is characterized by a strong vertical gradient in the abundance of the related condensing species. On Earth, the ensuing gradient of mean molecular weight has relatively few dynamical consequences because N$_2$ is heavier than water vapor, so that only the release of latent heat significantly impacts convection. On the contrary, in an hydrogen dominated atmosphere (e.g. giant planets), all condensing species are significantly heavier than the background gas. This can stabilize the atmosphere against convection near a cloud deck if the enrichment in the given species exceeds a critical threshold. This raises two questions. What is transporting energy in such a stabilized layer, and how affected can the thermal profile of giant planets be? To answer these questions, we first carry out a linear analysis of the convective and double-diffusive instabilities in a condensable medium showing that an efficient condensation can suppress double-diffusive convection. This suggests that a stable radiative layer can form near a cloud condensation level, leading to an increase in the temperature of the deep adiabat. Then, we investigate the impact of the condensation of the most abundant species---water---with a steady-state atmosphere model. Compared to standard models, the temperature increase can reach several hundred degrees at the quenching depth of key chemical tracers. Overall, this effect could have many implications for our understanding of the dynamical and chemical state of the atmosphere of giant planets, for their future observations (with Juno for example), and for their internal evolution.

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.

Stellar activity and convection-related surface structures might cause bias in planet detection and characterization that use these transits. Surface convection simulations help to quantify the granulation signal. We used realistic three-dimensional radiative hydrodynamical simulations from the Stagger grid and synthetic images computed with the radiative transfer code Optim3D to model the transits of three prototype planets: a hot Jupiter, a hot Neptune, and a terrestrial planet. We computed intensity maps from RHD simulations of the Sun and a K-dwarf star at different wavelength bands from optical to far-infrared. We modeled the transit using synthetic stellar-disk images and emulated the temporal variation of the granulation intensity. We identified two types of granulation noise that act simultaneously during the planet transit: (i) the intrinsic change in the granulation pattern with timescales smaller than the usual planet transit, and (ii) the fact that the transiting planet occults isolated regions of the photosphere that differ in local surface brightness. Our modeling approach shows that the granulation pattern has a non-negligible effect on the light curve depth during the transit, and, consequentially on the determination of the radius of the planet transiting. The granulation noise appears to be correlated among the different wavelength ranges either in the visible or in the infrared regions. The granulation has to be considered as an intrinsic uncertainty (as a result of stellar variability) on the precise measurements of exoplanet transits of planets. The full characterization of the granulation is essential for determining the degree of uncertainty on the planet parameters. In this context, the use of 3D RHD simulations is important to measure the convection-related fluctuations.

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.

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.

We report on spectroscopic observations of Saturn's stratosphere in July 2011 with the Texas Echelon Cross Echelle Spectrograph (TEXES) mounted on the NASA InfraRed Telescope Facility (IRTF). The observations, targeting several lines of the CH$_4$ $\nu_4$ band and the H$_2$ S(1) quadrupolar line, were designed to determine how Saturn's stratospheric thermal structure was disturbed by the 2010 Great White Spot. A study of Cassini Composite Infrared Spectrometer (CIRS) spectra had already shown the presence of a large stratospheric disturbance centered at a pressure of 2~hPa, nicknamed the beacon B0, and a tail of warm air at lower pressures (Fletcher et al. 2012. Icarus 221, 560--586). Our observations confirm that the beacon B0 vertical structure determined by CIRS, with a maximum temperature of $180\pm1$K at 2~hPa, is overlain by a temperature decrease up to the 0.2-hPa pressure level. Our retrieved maximum temperature of $180\pm1$K is colder than that derived by CIRS ($200\pm1$K), a difference that may be quantitatively explained by terrestrial atmospheric smearing. We propose a scenario for the formation of the beacon based on the saturation of gravity waves emitted by the GWS. Our observations also reveal that the tail is a planet-encircling disturbance in Saturn's upper stratosphere, oscillating between 0.2 and 0.02~hPa, showing a distinct wavenumber-2 pattern. We propose that this pattern in the upper stratosphere is either the signature of thermal tides generated by the presence of the warm beacon in the mid-stratosphere, or the signature of Rossby wave activity.

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.

General circulation models of the atmosphere of hot Jupiters have shown the existence of a supersonic eastward equatorial jet. In this paper, we investigate the effects of compressibility on the atmospheric dynamics by solving the standard Euler equations. This is done by means of a series of simulations performed in the framework of the equatorial beta-plane approximation using the finite volume shock-capturing code RAMSES. At low resolution, we recover the classical results described in the literature: we find a strong and steady supersonic equatorial jet of a few km/s that displays no signature of shocks. We next show that the jet zonal velocity depends significantly on the grid meridional resolution. When that resolution is fine enough to properly resolve the jet, the latter is subject to a Kelvin-Helmholtz instability. The jet zonal mean velocity displays regular oscillations with a typical timescale of few days and a significant amplitude of about 15% of the jet velocity. We also find compelling evidence for the development of a vertical shear instability at pressure levels of a few bars. It seems to be responsible for an increased downward kinetic energy flux, significantly affecting the temperature of the deep atmosphere, and appears to act as a form of drag on the equatorial jet. This instability also creates velocity fluctuations that propagate upward and steepen into weak shocks at pressure levels of a few mbars. We conclude that hot Jupiter equatorial jets are potentially unstable to both a barotropic Kelvin-Helmholtz instability and a vertical shear instability. Upon confirmation using more realistic models, both instabilities could result in significant time variability of the atmospheric winds, may provide a small scale dissipation mechanism in the flow, and might have consequences for the internal evolution of hot Jupiters.

Clouds have been shown to be present in many exoplanetary atmospheres. Cloud formation modeling predicts considerable inhomogeneities of cloud cover, consistent with optical phase curve observations. However, optical phase curves cannot resolve some existing degeneracies between cloud location and cloud optical properties. We present a conceptually simple technique to detect inhomogeneous cloud cover on exoplanets. Such an inhomogeneous cloud cover produces an asymmetric primary transit of the planet in front of the host star. Asymmetric transits produce characteristic residuals compared to a standard symmetric model. Furthermore, bisector spans can be used to determine asymmetries in the transit light curve. We apply a model of asymmetric transits to the light curves of HAT-P-7b, Kepler-7b and HD209458b and search for possible cloud signatures. The nearly uninterrupted Kepler photometry is particularly well-suited for this method since it allows for a very high time resolution. We do not find any statistically sound cloud signature in the data of the considered planets. For HAT-P-7b, a tentative detection of an asymmetric cloud cover is found, consistent with analysis of the optical phase curve. Based on Bayesian probability arguments, a symmetric model with an offset in the transit ephemeris remains, however, the most viable model. Still, this work demonstrates that for suitable targets, namely low-gravity planets around bright stars, the method can be used to constrain cloud cover characteristics and is thus a helpful additional tool to study exoplanetary atmospheres.

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.

The warm sub-Neptune GJ1214b has a featureless transit spectrum which may be due to the presence of high and thick clouds or haze. Here, we simulate the atmosphere of GJ1214b with a 3D General Circulation Model for cloudy hydrogen-dominated atmospheres, including cloud radiative effects. We show that the atmospheric circulation is strong enough to transport micrometric cloud particles to the upper atmosphere and generally leads to a minimum of cloud at the equator. By scattering stellar light, clouds increase the planetary albedo to 0.4-0.6 and cool the atmosphere below 1 mbar. However, the heating by ZnS clouds leads to the formation of a stratospheric thermal inversion above 10 mbar, with temperatures potentially high enough on the dayside to evaporate KCl clouds. We show that flat transit spectra consistent with HST observations are possible if cloud particle radii are around 0.5 micron, and that such clouds should be optically thin at wavelengths > 3 microns. Using simulated cloudy atmospheres that fit the observed spectra we generate transit, emission and reflection spectra and phase curves for GJ1214b. We show that a stratospheric thermal inversion would be readily accessible in near and mid-infrared atmospheric spectral windows. We find that the amplitude of the thermal phase curves is strongly dependent on metallicity, but only slightly impacted by clouds. Our results suggest that primary and secondary eclipses and phase curves observed by the James Webb Space Telescope in the near to mid-infrared should provide strong constraints on the nature of GJ1214b's atmosphere and clouds.

GJ1214b is a warm sub-Neptune transiting in front of a nearby M dwarf star. Recent observations indicate the presence of high and thick clouds or haze whose presence requires strong atmospheric mixing. In order to understand the transport and distribution of such clouds/haze, we study the atmospheric circulation and the vertical mixing of GJ1214b with a 3D General Circulation Model for cloud-free hydrogen-dominated atmospheres (metallicity of 1, 10 and 100 times the solar value) and for a water-dominated atmosphere. We analyze the effect of the atmospheric metallicity on the thermal structure and zonal winds. We also analyze the zonal mean meridional circulation and show that it corresponds to an anti-Hadley circulation in most of the atmosphere with upwelling at mid-latitude and downwelling at the equator in average. This circulation must be present on a large range of synchronously rotating exoplanets with strong impact on cloud formation and distribution. Using simple tracers, we show that vertical winds on GJ1214b can be strong enough to loft micrometric particles and that the anti-Hadley circulation leads to a minimum of tracers at the equator. We find that the strength of the vertical mixing increases with metallicity. We derive 1D equivalent eddy diffusion coefficients and find simple parametrizations from Kzz=7x10^2xP_{bar}^{-0.4} m^2/s for solar metallicity to Kzz=3x10^3xP_{bar}^{-0.4} m^2/s for the 100xsolar metallicity. These values should favor an efficient formation of photochemical haze in the upper atmosphere of GJ1214b.

Planets in the habitable zone of lower-mass stars are often assumed to be in a state of tidally synchronized rotation, which would considerably affect their putative habitability. Although thermal tides cause Venus to rotate retrogradely, simple scaling arguments tend to attribute this peculiarity to the massive Venusian atmosphere. Using a global climate model, we show that even a relatively thin atmosphere can drive terrestrial planets' rotation away from synchronicity. We derive a more realistic atmospheric tide model that predicts four asynchronous equilibrium spin states, two being stable, when the amplitude of the thermal tide exceeds a threshold that is met for habitable Earth-like planets with a 1-bar atmosphere around stars more massive than 0.5-0.7Msun. Thus, many recently discovered terrestrial planets could exhibit asynchronous spin-orbit rotation, even with a thin atmosphere.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

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.

We present the radial-velocity follow-up of two Kepler planetary transiting candidates (KOI-189 and KOI-686) carried out with the SOPHIE spectrograph at the Observatoire de Haute Provence. These data promptly discard these objects as viable planet candidates and show that the transiting objects are in the regime of very low-mass stars, where a strong discrepancy between observations and models persists for the mass and radius parameters. By combining the SOPHIE spectra with the Kepler light curve and photometric measurements found in the literature, we obtain a full characterization of the transiting companions, their orbits, and their host stars. The two companions are in significantly eccentric orbits with relatively long periods (30 days and 52.5 days), which makes them suitable objects for a comparison with theoretical models, since the effects invoked to understand the discrepancy with observations are weaker for these orbital distances. KOI-189 B has a mass M = 0.0745 +/- 0.0033 Msun and a radius R = 0.1025 +/- 0.0024 Rsun. The density of KOI-189 B is significantly lower than expected from theoretical models for a system of its age. We explore possible explanations for this difference. KOI-189 B is the smallest hydrogen-burning star with such a precise determination of its fundamental parameters. KOI-686 B is larger and more massive (M = 0.0915 +/- 0.0043 Msun; R = 0.1201 +/- 0.0033 Rsun), and its position in the mass-radius diagram agrees well with theoretical expectations.

Super-Earths with solid surfaces such CoRot-7b and Kepler-10b are expected to be extremely hot. It has been suggested that they must be atmosphere-free and that a lava ocean is present on their hot dayside. Here, we use several dedicated thermal models to study how observations with NIRSPEC on the JWST could further confirm and constrain, or reject the atmosphere-free lava ocean planet model for very hot super-Earths. Using CoRoT-7b as a working case, we explore the consequences on the phase-curve of a non tidal-locked rotation, with/without an atmosphere, and for different values of the albedo. We simulate future observations of the reflected light and thermal emission from CoRoT-7b with NIRSPEC-JWST and look for detectable signatures, such as time lag, as well as the possibility to retrieve the latitudinal surface temperature distribution. We demonstrate that we should be able to constrain several parameters after observations of two orbits (42 h) with a broad range of wavelengths: i)The Bond albedo is retrieved to within +/- 0.03 in most cases. ii) The lag effect allows us to retrieve the rotation period within 3 hours for a planet, with rotation half the orbital period. iii) Any spin period shorter than a limit in the range 30 - 800 h, depending on the thickness of the thermal layer in the soil, would be detected. iv) The presence of a thick gray atmosphere with a pressure of one bar, and a specific opacity higher than 1.0E-5 m-2 kg-1 is detectable. v) With spectra up to 4.5 mu, the latitudinal temperature profile can be retrieved to within 30 K with a risk of a totally wrong solution in 5 % of the cases. We conclude that it should thus be possible to distinguish the reference situation from other cases. Surface temperature map and the albedo brings important constraints on the nature or the physical state of the soil of hot super-Earths.

Because the solar luminosity increases over geological timescales, Earth climate is expected to warm, increasing water evaporation which, in turn, enhances the atmospheric greenhouse effect. Above a certain critical insolation, this destabilizing greenhouse feedback can "runaway" until all the oceans are evaporated. Through increases in stratospheric humidity, warming may also cause oceans to escape to space before the runaway greenhouse occurs. The critical insolation thresholds for these processes, however, remain uncertain because they have so far been evaluated with unidimensional models that cannot account for the dynamical and cloud feedback effects that are key stabilizing features of Earth's climate. Here we use a 3D global climate model to show that the threshold for the runaway greenhouse is about 375 W/m$^2$, significantly higher than previously thought. Our model is specifically developed to quantify the climate response of Earth-like planets to increased insolation in hot and extremely moist atmospheres. In contrast with previous studies, we find that clouds have a destabilizing feedback on the long term warming. However, subsident, unsaturated regions created by the Hadley circulation have a stabilizing effect that is strong enough to defer the runaway greenhouse limit to higher insolation than inferred from 1D models. Furthermore, because of wavelength-dependent radiative effects, the stratosphere remains cold and dry enough to hamper atmospheric water escape, even at large fluxes. This has strong implications for Venus early water history and extends the size of the habitable zone around other stars.

What kind of environment may exist on terrestrial planets around other stars? In spite of the lack of direct observations, it may not be premature to speculate on exoplanetary climates, for instance to optimize future telescopic observations, or to assess the probability of habitable worlds. To first order, climate primarily depends on 1) The atmospheric composition and the volatile inventory; 2) The incident stellar flux; 3) The tidal evolution of the planetary spin, which can notably lock a planet with a permanent night side. The atmospheric composition and mass depends on complex processes which are difficult to model: origins of volatile, atmospheric escape, geochemistry, photochemistry. We discuss physical constraints which can help us to speculate on the possible type of atmosphere, depending on the planet size, its final distance for its star and the star type. Assuming that the atmosphere is known, the possible climates can be explored using Global Climate Models analogous to the ones developed to simulate the Earth as well as the other telluric atmospheres in the solar system. Our experience with Mars, Titan and Venus suggests that realistic climate simulators can be developed by combining components like a "dynamical core", a radiative transfer solver, a parametrisation of subgrid-scale turbulence and convection, a thermal ground model, and a volatile phase change code. On this basis, we can aspire to build reliable climate predictors for exoplanets. However, whatever the accuracy of the models, predicting the actual climate regime on a specific planet will remain challenging because climate systems are affected by strong positive destabilizing feedbacks (such as runaway glaciations and runaway greenhouse effect). They can drive planets with very similar forcing and volatile inventory to completely different states.

Different solutions have been proposed to solve the "faint young Sun problem", defined by the fact that the Earth was not fully frozen during the Archean despite the fainter Sun. Most previous studies were performed with simple 1-D radiative convective models and did not account well for the clouds and ice-albedo feedback or the atmospheric and oceanic transport of energy. We apply a global climate model (GCM) to test the different solutions to the faint young Sun problem. We explore the effect of greenhouse gases (CO2 and CH4), atmospheric pressure, cloud droplet size, land distribution, and Earth's rotation rate. We show that neglecting organic haze, 100 mbar of CO2 with 2 mbar of CH4 at 3.8 Ga and 10 mbar of CO2 with 2 mbar of CH4 at 2.5 Ga allow a temperate climate (mean surface temperature between 10{\deg}C and 20{\deg}C). Such amounts of greenhouse gases remain consistent with the geological data. Removing continents produces a warming lower than +4{\deg}C. The effect of rotation rate is even more limited. Larger droplets (radii of 17 microns versus 12 microns) and a doubling of the atmospheric pressure produce a similar warming of around +7{\deg}C. In our model, ice-free water belts can be maintained up to 25{\deg}N/S with less than 1 mbar of CO2 and no methane. An interesting cloud feedback appears above cold oceans, stopping the glaciation. Such a resistance against full glaciation tends to strongly mitigate the faint young Sun problem.

Several studies have already considered the influence of tides on the evolution of systems composed of a star and a close-in companion to tentatively explain different observations such as the spin-up of some stars with hot Jupiters, the radius anomaly of short orbital period planets and the synchronization or quasi-synchronization of the stellar spin in some extreme cases. However, the nature of the mechanism responsible for the tidal dissipation in such systems remains uncertain. In this paper, we claim that the so-called elliptical instability may play a major role in these systems, explaining some systematic features present in the observations. This hydrodynamic instability, arising in rotating flows with elliptical streamlines, is suspected to be present in both planet and star of such systems, which are elliptically deformed by tides. The presence and the influence of the elliptical instability in gaseous bodies, such as stars or hot Jupiters, are most of the time neglected. In this paper, using numerical simulations and theoretical arguments, we consider several features associated to the elliptical instability in hot-Jupiter systems. In particular, the use of ad hoc boundary conditions makes it possible to estimate the amplitude of the elliptical instability in gaseous bodies. We also consider the influence of compressibility on the elliptical instability, and compare the results to the incompressible case. We demonstrate the ability for the elliptical instability to grow in the presence of differential rotation, with a possible synchronized latitude, provided that the tidal deformation and/or the rotation rate of the fluid are large enough. Moreover, the amplitude of the instability for a centrally-condensed mass of fluid is of the same order of magnitude as for an incompressible fluid for a given distance to the threshold of the instability. Finally, we show that the assumption of the elliptical instability being the main tidal dissipation process in eccentric inflated hot Jupiters and misaligned stars is consistent with current data.

We propose to use Kepler in 2-wheel mode to conduct a detailed search for Earth-sized planets orbiting ultra-cool stars and brown dwarfs (spectral types from M7 to L3). This population of objects presents several advantages for exoplanet surveys. First, ultra-cool stars and brown dwarfs are small and thus result in favorable planet-to-star area ratios. Second, because of their low effective temperature, the inner edge of their habitable zone is extremely close (2 to 3 days only). Third, our targets are bright at infrared wavelengths, which will enable detailed follow-up studies. Our program therefore represents a unique opportunity to find a transiting Earth-size exoplanet for which atmospheric features (including biosignatures) could be detected with near-to-come facilities such as JWST. Such exoplanet has not been discovered yet. Kepler in 2-wheel mode provides the required stability and photometric precision to make this survey successful. Our initial target sample includes 60 ultra-cool stars and brown dwarfs from which we expect to detect at least one transiting planet. We propose to monitor each source for 4 days, resulting in a total program duration of ~240 days.

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.

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.

As they keep cooling and contracting, Solar System giant planets radiate more energy than they receive from the Sun. Applying the first and second principles of thermodynamics, one can determine their cooling rate, luminosity, and temperature at a given age. Measurements of Saturn's infrared intrinsic luminosity, however, reveal that this planet is significantly brighter than predicted for its age. This excess luminosity is usually attributed to the immiscibility of helium in the hydrogen-rich envelope, leading to "rains" of helium-rich droplets. Existing evolution calculations, however, suggest that the energy released by this sedimentation process may not be sufficient to resolve the puzzle. Here, we demonstrate using planetary evolution models that the presence of layered convection in Saturn's interior, generated, like in some parts of Earth oceans, by the presence of a compositional gradient, significantly reduces its cooling. It can explain the planet's present luminosity for a wide range of configurations without invoking any additional source of energy. This suggests a revision of the conventional homogeneous adiabatic interior paradigm for giant planets, and questions our ability to assess their heavy element content. This reinforces the possibility for layered convection to help explaining the anomalously large observed radii of extrasolar giant planets.

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.

The inner edge of the classical habitable zone is often defined by the critical flux needed to trigger the runaway greenhouse instability. This 1D notion of a critical flux, however, may not be so relevant for inhomogeneously irradiated planets, or when the water content is limited (land planets). Here, based on results from our 3D global climate model, we find that the circulation pattern can shift from super-rotation to stellar/anti stellar circulation when the equatorial Rossby deformation radius significantly exceeds the planetary radius. Using analytical and numerical arguments, we also demonstrate the presence of systematic biases between mean surface temperatures or temperature profiles predicted from either 1D or 3D simulations. Including a complete modeling of the water cycle, we further demonstrate that for land planets closer than the inner edge of the classical habitable zone, two stable climate regimes can exist. One is the classical runaway state, and the other is a collapsed state where water is captured in permanent cold traps. We identify this "moist" bistability as the result of a competition between the greenhouse effect of water vapor and its condensation. We also present synthetic spectra showing the observable signature of these two states. Taking the example of two prototype planets in this regime, namely Gl581c and HD85512b, we argue that they could accumulate a significant amount of water ice at their surface. If such a thick ice cap is present, gravity driven ice flows and geothermal flux should come into play to produce long-lived liquid water at the edge and/or bottom of the ice cap. Consequently, the habitability of planets at smaller orbital distance than the inner edge of the classical habitable zone cannot be ruled out. Transiting planets in this regime represent promising targets for upcoming observatories like EChO and JWST.

On the basis of geological evidence, it is often stated that the early martian climate was warm enough for liquid water to flow on the surface thanks to the greenhouse effect of a thick atmosphere. We present 3D global climate simulations of the early martian climate performed assuming a faint young sun and a CO2 atmosphere with pressure between 0.1 and 7 bars. The model includes a detailed radiative transfer model using revised CO2 gas collision induced absorption properties, and a parameterisation of the CO2 ice cloud microphysical and radiative properties. A wide range of possible climates is explored by using various values of obliquities, orbital parameters, cloud microphysic parameters, atmospheric dust loading, and surface properties. Unlike on present day Mars, for pressures higher than a fraction of a bar, surface temperatures vary with altitude because of the adiabatic cooling and warming of the atmosphere when it moves vertically. In most simulations, CO2 ice clouds cover a major part of the planet but greenhouse effect does not exceed +15 K. We find that a CO2 atmosphere could not have raised the annual mean temperature above 0{\deg}C anywhere on the planet. The collapse of the atmosphere into permanent CO2 ice caps is predicted for pressures higher than 3 bar, or conversely at pressure lower than one bar if the obliquity is low enough. Summertime diurnal mean surface temperatures above 0{\deg}C (a condition which could have allowed rivers to form) are predicted for obliquity larger than 40{\deg} at high latitudes but not in locations where most valley networks are observed. In the absence of other warming mechanisms, our climate model results are thus consistent with a cold early Mars scenario in which non climatic mechanisms must occur to explain the evidence for liquid water. In a companion paper by Wordsworth et al., we simulate the hydrological cycle on such a planet.

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.

While conventional interior models for Jupiter and Saturn are based on the simplistic assumption of a solid core surrounded by a homogeneous gaseous envelope, we derive new models with an inhomogeneous distribution of heavy elements, i.e. a gradient of composition, within these planets. Such a compositional stratification hampers large scale convection which turns into double-diffusive convection, yielding an inner thermal profile which departs from the traditionally assumed adiabatic interior, affecting these planet heat content and cooling history. To address this problem, we develop an analytical approach of layered double-diffusive convection and apply this formalism to Solar System gaseous giant planet interiors. These models satisfy all observational constraints and yield a metal enrichment for our gaseous giants up to 30 to 60% larger than previously thought. The models also constrain the size of the convective layers within the planets. As the heavy elements tend to be redistributed within the gaseous envelope, the models predict smaller than usual central cores inside Saturn and Jupiter, with possibly no core for this latter. These models open a new window and raise new challenges on our understanding of the internal structure of giant (solar and extrasolar) planets, in particular on the determination of their heavy material content, a key diagnostic for planet formation theories.

We quantify the systematic impact of the non-spherical shape of transiting planets and brown dwarfs, due to tidal forces and rotation, on the observed transit depth. Such a departure from sphericity leads to a bias in the derivation of the transit radius from the light curve and affects the comparison with planet structure and evolution models which assume spherical symmetry. As the tidally deformed planet projects its smallest cross section area during the transit, the measured effective radius is smaller than the one of the unperturbed spherical planet. This effect can be corrected by calculating the theoretical shape of the observed planet. We derive simple analytical expressions for the ellipsoidal shape of a fluid object (star or planet) accounting for both tidal and rotational deformations and calibratre it with fully numerical evolution models in the 0.3Mjup-75Mjup mass range. Our calculations yield a 20% effect on the transit depth, i.e. a 10% decrease of the measured radius, for the extreme case of a 1Mjup planet orbiting a Sun-like star at 0.01AU. For the closest planets detected so far (< 0.05 AU), the effect on the radius is of the order of 1 to 10%, by no means a negligible effect, enhancing the puzzling problem of the anomalously large bloated planets. These corrections must thus be taken into account for a correct determination of the radius from the transit light curve. Our analytical expressions can be easily used to calculate these corrections, due to the non-spherical shape of the planet, on the observed transit depth and thus to derive the planet's real equilibrium radius. They can also be used to model ellipsoidal variations of the stellar flux now detected in the CoRoT and Kepler light curves. We also derive directly usable analytical expressions for the moment of inertia, oblateness and Love number (k_2) of a fluid planet as a function of its mass.

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.

Stellar radiation has conservatively been used as the key constraint to planetary habitability. We review here the effects of tides, exerted by the host star on the planet, on the evolution of the planetary spin. Tides initially drive the rotation period and the orientation of the rotation axis into an equilibrium state but do not necessarily lead to synchronous rotation. As tides also circularize the orbit, eventually the rotation period does equal the orbital period and one hemisphere will be permanently irradiated by the star. Furthermore, the rotational axis will become perpendicular to the orbit, i.e. the planetary surface will not experience seasonal variations of the insolation. We illustrate here how tides alter the spins of planets in the traditional habitable zone. As an example, we show that, neglecting perturbations due to other companions, the Super-Earth Gl581d performs two rotations per orbit and that any primordial obliquity has been eroded.

Stellar insolation has been used as the main constraint on a planet's habitability. However, as more Earth-like planets are discovered around low-mass stars (LMSs), a re-examination of the role of tides on the habitability of exoplanets has begun. Those studies have yet to consider the misalignment between a planet's rotational axis and the orbital plane normal, i.e. the planetary obliquity. We apply two equilibrium tide theories to compute the obliquity evolution of terrestrial planets orbiting in the habitable zones around LMSs. The time for the obliquity to decrease from an Earth-like obliquity of 23.5 deg to 5 deg, the 'tilt erosion time', is compared to the traditional insolation habitable zone (IHZ) as a function of semi-major axis, eccentricity, and stellar mass. We also compute tidal heating and equilibrium rotation caused by obliquity tides. The Super-Earth Gl581d and the planet candidate Gl581g are studied as examples for tidal processes. Earth-like obliquities of terrestrial planets in the IHZ around <0.25 solar mass (M_sun) stars are eroded in <0.1 Gyr. Only terrestrial planets orbiting stars with masses >0.9 M_sun experience tilt erosion times larger than 1 Gyr throughout the IHZ. Terrestrial planets in the IHZ of stars with masses <0.25 M_sun undergo significant tidal heating due to obliquity tides. The predictions of the two tidal models diverge significantly for e>0.3. In our two-body simulations, Gl581d's obliquity is eroded to 0 and its rotation period reached its equilibrium state of half its orbital period in <0.1 Gyr. Tidal surface heating on the putative Gl581g is <150 mW/m^2 as long as its eccentricity is <0.3. Obliquity tides modify the concept of the habitable zone. Tilt erosion of terrestrial planets orbiting LMSs should be included by atmospheric modelers. Tidal heating needs to be considered by geologists.

In this short review, we summarize our present understanding (and non-understanding) of exoplanet formation, structure and evolution, in the light of the most recent discoveries. Recent observations of transiting massive brown dwarfs seem to remarkably confirm the predicted theoretical mass-radius relationship in this domain. This mass-radius relationship provides, in some cases, a powerful diagnostic to distinguish planets from brown dwarfs of same mass, as for instance for Hat-P-20b. If confirmed, this latter observation shows that planet formation takes place up to at least 8 Jupiter masses. Conversely, observations of brown dwarfs down to a few Jupiter masses in young, low-extinction clusters strongly suggest an overlapping mass domain between (massive) planets and (low-mass) brown dwarfs, i.e. no mass edge between these two distinct (in terms of formation mechanism) populations. At last, the large fraction of heavy material inferred for many of the transiting planets confirms the core-accretion scenario as been the dominant one for planet formation.

The recent detection of the transit of very massive substellar companions (CoRoT-3b, Deleuil et al. 2008; CoRoT-15b, Bouchy et al. 2010; WASP-30b, Anderson et al. 2010; Hat-P-20b, Bakos et al. 2010) provides a strong constraint to planet and brown dwarf formation and migration mechanisms. Whether these objects are brown dwarfs originating from the gravitational collapse of a dense molecular cloud that, at the same time, gave birth to the more massive stellar companion, or whether they are planets that formed through core accretion of solids in the protoplanetary disk can not always been determined unambiguously and the mechanisms responsible for their short orbital distances are not yet fully understood. In this contribution, we examine the possibility to constrain the nature of a massive substellar object from the various observables provided by the combination of Radial Velocity and Photometry measurements (e.g. M_p, R_p, M_s, Age, a, e...). In a second part, developments in the modeling of tidal evolution at high eccentricity and inclination - as measured for HD 80 606 with e=0.9337 (Naef et al. 2001), XO-3 with a stellar obliquity >37.3+-3.7 deg (H\'ebrard et al. 2008; Winn et al. 2009) and several other exoplanets - are discussed along with their implication in the understanding of the radius anomaly problem of extrasolar giant planets.

The Lyot project used an optimized Lyot coronagraph with Extreme Adaptive Optics at the 3.63m Advanced Electro-Optical System telescope (AEOS) to observe 86 stars from 2004 to 2007. In this paper we give an overview of the survey results and a statistical analysis of the observed nondetections around 58 of our targets to place constraints on the population of substellar companions to nearby stars. The observations did not detect any companion in the substellar regime. Since null results can be as important as detections, we analyzed each observation to determine the characteristics of the companions that can be ruled out. For this purpose we use a Monte Carlo approach to produce artificial companions, and determine their detectability by comparison with the sensitivity curve for each star. All the non-detection results are combined using a Bayesian approach and we provide upper limits on the population of giant exoplanets and brown dwarfs for this sample of stars. Our nondetections confirm the rarity of brown dwarfs around solar-like stars and we constrain the frequency of massive substellar companions (M>40Mjup) at orbital separation between and 10 and 50 AU to be <20%.

In this paper, we present the consistent evolution of short-period exoplanets coupling the tidal and gravothermal evolution of the planet. Contrarily to previous similar studies, our calculations are based on the complete tidal evolution equations of the Hut model, valid at any order in eccentricity, obliquity and spin. We demonstrate, both analytically and numerically, that, except if the system was formed with a nearly circular orbit (e<0.2), solving consistently the complete tidal equations is mandatory to derive correct tidal evolution histories. We show that calculations based on tidal models truncated at second order in eccentricity, as done in all previous studies, lead to erroneous tidal evolutions. As a consequence, tidal energy dissipation rates are severely underestimated in all these calculations and the characteristic timescales for the various orbital parameters evolutions can be wrong by up to three orders in magnitude. Based on these complete, consistent calculations, we revisit the viability of the tidal heating hypothesis to explain the anomalously large radius of transiting giant planets. We show that, even though tidal dissipation does provide a substantial contribution to the planet's heat budget and can explain some of the moderately bloated hot-Jupiters, this mechanism can not explain alone the properties of the most inflated objects, including HD 209458b. Indeed, solving the complete tidal equations shows that enhanced tidal dissipation and thus orbit circularization occur too early during the planet's evolution to provide enough extra energy at the present epoch. In that case another mechanisms, such as stellar irradiation induced surface winds dissipating in the planet's tidal bulges, or inefficient convection in the planet's interior must be invoked, together with tidal dissipation, to provide all the pieces of the abnormally large exoplanet puzzle.

We present detailed structure and evolution calculations for the first transiting extrasolar planets discovered by the space-based CoRoT mission. Comparisons between theoretical and observed radii provide information on the internal composition of the CoRoT objects. We distinguish three different categories of planets emerging from these discoveries and from previous ground-based surveys: (i) planets explained by standard planetary models including irradiation, (ii) abnormally bloated planets and (iii) massive objects belonging to the overlapping mass regime between planets and brown dwarfs. For the second category, we show that tidal heating can explain the relevant CoRoT objects, providing non-zero eccentricities. We stress that the usual assumption of a quick circularization of the orbit by tides, as usually done in transit light curve analysis, is not justified a priori, as suggested recently by Levrard et al. (2009), and that eccentricity analysis should be carefully redone for some observations. Finally, special attention is devoted to CoRoT-3b and to the identification of its very nature: giant planet or brown dwarf ? The radius determination of this object confirms the theoretical mass-radius predictions for gaseous bodies in the substellar regime but, given the present observational uncertainties, does not allow an unambiguous identification of its very nature. This opens the avenue, however, to an observational identification of these two distinct astrophysical populations, brown dwarfs and giant planets, in their overlapping mass range, as done for the case of the 8 Jupiter-mass object Hat-P-2b. (abridged)

In this review, we summarize our present knowledge of the behaviour of the mass-radius relationship from solar-type stars down to terrestrial planets, across the regime of substellar objects, brown dwarfs and giant planets. Particular attention is paid to the identification of the main physical properties or mechanisms responsible for this behaviour. Indeed, understanding the mechanical structure of an object provides valuable information about its internal structure, composition and heat content as well as its formation history. Although the general description of these properties is reasonably well mastered, disagreement between theory and observation in certain cases points to some missing physics in our present modelling of at least some of these objects. The mass-radius relationship in the overlaping domain between giant planets and low-mass brown dwarfs is shown to represent a powerful diagnostic to distinguish between these two different populations and shows once again that the present IAU distinction between these two populations at a given mass has no valid foundation.