Researcher profile

Jérémy Leconte

Jérémy Leconte contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
27works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2021arXiv

Formation and dynamics of water clouds on temperate sub-Neptunes: The example of K2-18b

Hubble (HST) spectroscopic transit observations of the temperate sub-Neptune K2-18b were interpreted as the presence of water vapour with potential water clouds. 1D modelling studies also predict the formation of water clouds at some conditions. However, such models cannot predict the cloud cover, driven by atmospheric dynamics and thermal contrasts, and thus their real impact on spectra. The main goal of this study is to understand the formation, distribution and observational consequences of water clouds on K2-18b and other temperate sub-Neptunes. We simulated the atmospheric dynamics, water cloud formation and spectra of K2-18b for H2-dominated atmosphere using a 3D GCM. We analysed the impact of atmospheric composition (with metallicity from 1*solar to 1000*solar), concentration of cloud condensation nuclei and planetary rotation rate. Assuming that K2-18b has a synchronous rotation, we show that the atmospheric circulation in the upper atmosphere essentially corresponds to a symmetric day-to-night circulation. This regime preferentially leads to cloud formation at the substellar point or at the terminator. Clouds form for metallicity >100*solar with relatively large particles. For 100-300*solar metallicity, the cloud fraction at the terminators is small with a limited impact on transit spectra. For 1000*solar metallicity, very thick clouds form at the terminator. The cloud distribution appears very sensitive to the concentration of CCN and to the planetary rotation rate. Fitting HST transit data with our simulated spectra suggests a metallicity of ~100-300*solar. In addition, we found that the cloud fraction at the terminator can be highly variable, leading to a potential variability in transit spectra. This effect could be common on cloudy exoplanets and could be detectable with multiple transit observations.

preprint2021arXiv

Toward a multidimensional analysis of transmission spectroscopy. Part I: Computation of transmission spectra using a 1D, 2D, or 3D atmosphere structure

Considering the relatively high precision that will be reached by future observatories, it has recently become clear that one dimensional (1D) atmospheric models, in which the atmospheric temperature and composition of a planet are considered to vary only in the vertical, will be unable to represent exoplanetary transmission spectra with a sufficient accuracy. This is particularly true for warm to (ultra-) hot exoplanets because the atmosphere is unable to redistribute all the energy deposited on the dayside, creating a strong thermal and often compositional dichotomy on the planet. This situation is exacerbated by transmission spectroscopy, which probes the terminator region. This is the most heterogeneous region of the atmosphere. However, if being able to compute transmission spectra from 3D atmospheric structures (from a global climate model, e.g.) is necessary to predict realistic observables, it is too computationally expensive to be used in a data inversion framework. For this reason, there is a need for a medium-complexity 2D approach that captures the most salient features of the 3D model in a sufficiently fast implementation. With this in mind, we present a new open-source documented version of Pytmosph3R that handles the computation of transmission spectra for atmospheres with up to three spatial dimensions and can account for time variability. Taking the example of an ultra hot Jupiter, we illustrate how the changing orientation of the planet during the transit can allow us to probe the horizontal variations in the atmosphere. We further implement our algorithm in TauREx to allow the community to perform 2D retrievals. We describe our extensive cross-validation benchmarks and discuss the accuracy and numerical performance of each model.

preprint2020arXiv

A review of possible planetary atmospheres in the TRAPPIST-1 system

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.

preprint2020arXiv

Revised mass-radius relationships for water-rich rocky planets more irradiated than the runaway greenhouse limit

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.

preprint2020arXiv

Spectral binning of precomputed correlated-k coefficients

With the major increase in the volume of the spectroscopic line lists needed to perform accurate radiative transfer calculations, disseminating accurate radiative data has become almost as much a challenge as computing it. Considering that many planetary science applications are only looking for heating rates or mid-to-low resolution spectra, any approach enabling such computations in an accurate and flexible way at a fraction of the computing and storage costs is highly valuable. For many of these reasons, the correlated-k approach has become very popular. Its major weakness has been the lack of ways to adapt the spectral grid/resolution of precomputed k-coefficients, making it difficult to distribute a generic database suited for many different applications. Currently, most users still need to have access to a line-by-line transfer code with the relevant line lists or high-resolution cross sections to compute k-coefficient tables at the desired resolution. In this work, we demonstrate that precomputed k-coefficients can be binned to a lower spectral resolution without any additional assumptions, and show how this can be done in practice. We then show that this binning procedure does not introduce any significant loss in accuracy. Along the way, we quantify how such an approach compares very favorably with the sampled cross section approach. This opens up a new avenue to deliver accurate radiative transfer data by providing mid-resolution k-coefficient tables to users who can later tailor those tables to their needs on the fly. To help with this final step, we briefly present Exo_k, an open-access, open-source Python library designed to handle, tailor, and use many different formats of k-coefficient and cross-section tables in an easy and computationally efficient way.

preprint2020arXiv

Strong biases in retrieved atmospheric composition caused by strong day-night chemical heterogeneities

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.

preprint2020arXiv

The deep composition of Uranus and Neptune from in situ exploration and thermochemical modeling

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.

preprint2016arXiv

Habitability of planets on eccentric orbits: the limits of the mean flux approximation

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.

preprint2016arXiv

Shear-driven instabilities and shocks in the atmospheres of hot Jupiters

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.

preprint2016arXiv

Stratospheric aftermath of the 2010 Storm on Saturn as observed by the TEXES instrument. I. Temperature structure

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$ $ν_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.

preprint2015arXiv

3D modeling of GJ1214b's atmosphere: formation of inhomogeneous high clouds and observational implications

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.

preprint2015arXiv

3D modeling of GJ1214b's atmosphere: vertical mixing driven by an anti-Hadley circulation

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.

preprint2015arXiv

Asynchronous rotation of Earth-mass planets in the habitable zone of lower-mass stars

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.

preprint2014arXiv

Tidal evolution in multiple planet systems: application to Kepler-62 and Kepler-186

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.

preprint2013arXiv

3D climate modeling of close-in land planets: Circulation patterns, climate moist bistability and habitability

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.

preprint2013arXiv

Exploring the faint young Sun problem and the possible climates of the Archean Earth with a 3-D GCM

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°C and 20°C). Such amounts of greenhouse gases remain consistent with the geological data. Removing continents produces a warming lower than +4°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°C. In our model, ice-free water belts can be maintained up to 25°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.

preprint2013arXiv

Increased insolation threshold for runaway greenhouse processes on Earth like planets

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.

preprint2013arXiv

Layered convection as the origin of Saturn's luminosity anomaly

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.

preprint2013arXiv

Searching for Terrestrial Planets Orbiting in the Habitable Zone of Ultra-Cool Stars and Brown Dwarfs

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.

preprint2012arXiv

A new vision on giant planet interiors: the impact of double diffusive convection

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.

preprint2011arXiv

Distorted, non-spherical transiting planets: impact on the transit depth and on the radius determination

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&#39;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.

preprint2011arXiv

Habitability of extrasolar planets and tidal spin evolution

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.

preprint2011arXiv

Tidal evolution of planets around brown dwarfs

The tidal evolution of planets orbiting brown dwarfs (BDs) presents an interesting case study because BDs&#39; 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&#39;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&#39;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.

preprint2011arXiv

Tidal obliquity evolution of potentially habitable planets

Stellar insolation has been used as the main constraint on a planet&#39;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&#39;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 &#39;tilt erosion time&#39;, 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&#39;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.

preprint2010arXiv

Is tidal heating sufficient to explain bloated exoplanets? Consistent calculations accounting for finite initial eccentricity

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&#39;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&#39;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&#39;s tidal bulges, or inefficient convection in the planet&#39;s interior must be invoked, together with tidal dissipation, to provide all the pieces of the abnormally large exoplanet puzzle.

preprint2010arXiv

The Lyot Project Direct Imaging Survey of Substellar Companions: Statistical Analysis and Information from Nondetections

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

preprint2010arXiv

The radius anomaly in the planet/brown dwarf overlapping mass regime

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ébrard 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.