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 all ...that relevant for inhomogeneously irradiated planets, or when the water content is limited (land planets). Based on results from our 3D global climate model, we present general features of the climate and large-scale circulation on close-in terrestrial planets. 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, changing the redistribution properties of the atmosphere. Using analytical and numerical arguments, we also demonstrate the presence of systematic biases among mean surface temperatures and among temperature profiles predicted from either 1D or 3D simulations. After including a complete modeling of the water cycle, we further demonstrate that two stable climate regimes can exist for land planets closer than the inner edge of the classical habitable zone. One is the classical runaway state where all the water is vaporized, 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 on the night side or near the poles, highlighting the dynamical nature of the runaway greenhouse effect. We also present synthetic spectra showing the observable signature of these two states. Taking the example of two prototype planets in this regime, namely Gl 581 c and HD 85512 b, we argue that depending on the rate of water delivery and atmospheric escape during the life of these planets, they could accumulate a significant amount of water ice at their surface. If such a thick ice cap is present, various physical mechanisms observed on Earth (e.g., gravity driven ice flows, 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 exoplanet characterization observatories, such as EChO and JWST.
Context.
The atmospheric composition of exoplanets with masses between 2 and 10
M
⊕
is poorly understood. In that regard, the sub-Neptune K2-18b, which is subject to Earth-like stellar irradiation, ...offers a valuable opportunity for the characterisation of such atmospheres. Previous analyses of its transmission spectrum from the
Kepler
,
Hubble
(HST), and
Spitzer
space telescopes data using both retrieval algorithms and forward-modelling suggest the presence of H
2
O and an H
2
–He atmosphere, but have not detected other gases, such as CH
4
.
Aims.
We present simulations of the atmosphere of K2-18 b using Exo-REM, our self-consistent 1D radiative-equilibrium model, using a large grid of atmospheric parameters to infer constraints on its chemical composition.
Methods.
We compared the transmission spectra computed by our model with the above-mentioned data (0.4–5
μ
m), assuming an H
2
–He dominated atmosphere. We investigated the effects of irradiation, eddy diffusion coefficient, internal temperature, clouds, C/O ratio, and metallicity on the atmospheric structure and transit spectrum.
Results.
We show that our simulations favour atmospheric metallicities between 40 and 500 times solar and indicate, in some cases, the formation of H
2
O-ice clouds, but not liquid H
2
O clouds. We also confirm the findings of our previous study, which showed that CH
4
absorption features nominally dominate the transmission spectrum in the HST spectral range. We compare our results with results from retrieval algorithms and find that the H
2
O-dominated spectrum interpretation is either due to the omission of CH
4
absorptions or a strong overfitting of the data. Finally, we investigated different scenarios that would allow for a CH
4
-depleted atmosphere. We were able to fit the data to those scenarios, finding, however, that it is very unlikely for K2-18b to have a high internal temperature. A low C/O ratio (≈0.01–0.1) allows for H
2
O to dominate the transmission spectrum and can fit the data but so far, this set-up lacks a physical explanation. Simulations with a C/O ratio <0.01 are not able to fit the data satisfactorily.
ABSTRACT 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 on average. This circulation must be present on a large range of synchronously rotating exoplanets with a 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 for solar metallicity to for the 100× solar metallicity. These values should favor an efficient formation of photochemical haze in the upper atmosphere of GJ1214b.
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 – actually extends significantly towards 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 that 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 apply this tool to a simulation of the atmosphere of GJ 1214 b to produce synthetic JWST observations and show that producing a spectrum using only atmospheric columns at the terminator results in errors greater than expected noise. This demonstrates the necessity for a real 3D approach to model data for such precise observatories. Secondly, we investigate how day–night temperature gradients cause a systematic bias in retrieval analysis performed with 1D forward models. For that purpose we synthesise 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 – a result of the non-linear properties of atmospheric transmission. These biases will be difficult to detect as the 1D profiles used in the retrieval procedure are found to provide an excellent match to the observed spectra based on standard fitting criteria. This must be kept in mind when interpreting current and future data.
Context. The observation of planets in their formation stage is a crucial but very challenging step in understanding when, how, and where planets form. PDS 70 is a young pre-main sequence star ...surrounded by a transition disk, in the gap of which a planetary-mass companion has recently been discovered. This discovery represents the first robust direct detection of such a young planet, possibly still at the stage of formation. Aims. We aim to characterize the orbital and atmospheric properties of PDS 70 b, which was first identified on May 2015 in the course of the SHINE survey with SPHERE, the extreme adaptive-optics instrument at the VLT. Methods. We obtained new deep SPHERE/IRDIS imaging and SPHERE/IFS spectroscopic observations of PDS 70 b. The astrometric baseline now covers 6 yr, which allowed us to perform an orbital analysis. For the first time, we present spectrophotometry of the young planet which covers almost the entire near-infrared range (0.96–3.8 μm). We use different atmospheric models covering a large parameter space in temperature, log g, chemical composition, and cloud properties to characterize the properties of the atmosphere of PDS 70 b. Results. PDS 70 b is most likely orbiting the star on a circular and disk coplanar orbit at ~22 au inside the gap of the disk. We find a range of models that can describe the spectrophotometric data reasonably well in the temperature range 1000–1600 K and log g no larger than 3.5 dex. The planet radius covers a relatively large range between 1.4 and 3.7 RJ with the larger radii being higher than expected from planet evolution models for the age of the planet of 5.4 Myr. Conclusions. This study provides a comprehensive data set on the orbital motion of PDS 70 b, indicating a circular orbit and a motion coplanar with the disk. The first detailed spectral energy distribution of PDS 70 b indicates a temperature typical of young giant planets. The detailed atmospheric analysis indicates that a circumplanetary disk may contribute to the total planetflux.
Context. Hubble
Space Telescope (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 in K2-18b’s atmosphere in some conditions. However, such models cannot predict the cloud cover, which is driven by atmospheric dynamics and thermal contrasts, and thus neither can they predict the real impact of clouds on spectra.
Aims.
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.
Methods.
We simulated the atmospheric dynamics, water cloud formation, and spectra of K2-18b for a H
2
-dominated atmosphere using a 3D general circulation model. We analysed the impact of atmospheric composition (with metallicity from 1× solar to 1000× solar), concentration of cloud condensation nuclei, and planetary rotation rate.
Results.
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 with very efficient heat redistribution. This regime preferentially leads to cloud formation at the sub-stellar point or at the terminator. Clouds form at metallicity ≥100× solar with relatively large particles (radius = 30–450
μ
m). At 100–300× solar metallicity, the cloud fraction at the terminators is small with a limited impact on transit spectra. At 1000× solar metallicity, very thick clouds form at the terminator, greatly flattening the transit spectrum. The cloud distribution appears very sensitive to the concentration of cloud condensation nuclei and to the planetary rotation rate, although the impact on transit spectra is modest in the near-infrared. Fitting HST transit data with our simulated spectra suggests a metallicity of ~100–300× solar, which is consistent with the mass-metallicity trend of giant planets in the Solar System. In addition, we found that the cloud fraction at the terminator can be highly variable in some conditions, leading to a potential variability in transit spectra that is correlated with spectral windows. This effect could be common on cloudy exoplanets and could be detectable with multiple transit observations. Finally, the complex cloud dynamics revealed in this study highlight the inherent 3D nature of clouds shaped by couplings between microphysics, radiation, and atmospheric circulation.
We developed a simple, physical, and self-consistent cloud model for brown dwarfs and young giant exoplanets. We compared different parametrizations for the cloud particle size, by fixing either ...particle radii or the mixing efficiency (parameter fsed), or by estimating particle radii from simple microphysics. The cloud scheme with simple microphysics appears to be the best parametrization by successfully reproducing the observed photometry and spectra of brown dwarfs and young giant exoplanets. In particular, it reproduces the L-T transition, due to the condensation of silicate and iron clouds below the visible/near-IR photosphere. It also reproduces the reddening observed for low-gravity objects, due to an increase of cloud optical depth for low gravity. In addition, we found that the cloud greenhouse effect shifts chemical equilibrium, increasing the abundances of species stable at high temperature. This effect should significantly contribute to the strong variation of methane abundance at the L-T transition and to the methane depletion observed on young exoplanets. Finally, we predict the existence of a continuum of brown dwarfs and exoplanets for absolute J magnitude = 15-18 and color = 0-3, due to the evolution of the L-T transition with gravity. This self-consistent model therefore provides a general framework to understand the effects of clouds and appears well-suited for atmospheric retrievals.
Context.
The TRAPPIST-1 planetary system is favourable for transmission spectroscopy and offers the unique opportunity to study rocky planets with possibly non-primary envelopes. We present here the ...transmission spectrum of the seventh planet of the TRAPPIST-1 system, TRAPPIST-1 h (
R
P
= 0.752
R
⊕
,
T
eq
= 173 K) using
Hubble
Space Telescope (HST), Wide Field Camera 3 Grism 141 (WFC3/G141) data.
Aims.
Our purpose is to reduce the HST observations of the seventh planet of the TRAPPIST-1 system and, by testing a simple atmospheric hypothesis, to put a new constraint on the composition and the nature of the planet.
Methods.
First we extracted and corrected the raw data to obtain a transmission spectrum in the near-infrared (NIR) band (1.1–1.7 μm). TRAPPIST-1 is a cold M-dwarf and its activity could affect the transmission spectrum. We corrected for stellar modulations using three different stellar contamination models; while some fit the data better, they are statistically not significant and the conclusion remains unchanged concerning the presence or lack thereof of an atmosphere. Finally, using a Bayesian atmospheric retrieval code, we put new constraints on the atmosphere composition of TRAPPIST-1h.
Results.
According to the retrieval analysis, there is no evidence of molecular absorption in the NIR spectrum. This suggests the presence of a high cloud deck or a layer of photochemical hazes in either a primary atmosphere or a secondary atmosphere dominated by heavy species such as nitrogen. This result could even be the consequence of the lack of an atmosphere as the spectrum is better fitted using a flat line. Variations in the transit depth around 1.3 μm are likely due to remaining scattering noise and the results do not improve while changing the spectral resolution. TRAPPIST-1 h has probably lost its atmosphere or possesses a layer of clouds and hazes blocking the NIR signal. We cannot yet distinguish between a primary cloudy or a secondary clear envelope using HST/WFC3 data; however, in most cases with more than 3
σ
confidence, we can reject the hypothesis of a clear atmosphere dominated by hydrogen and helium. By testing the forced secondary atmospheric scenario, we find that a CO-rich atmosphere (i.e. with a volume mixing ratio of 0.2) is one of the best fits to the spectrum with a Bayes factor of 1.01, corresponding to a 2.1
σ
detection.
Context. PDS 70 is a young (5.4 Myr), nearby (~113 pc) star hosting a known transition disk with a large gap. Recent observations with SPHERE and NACO in the near-infrared (NIR) allowed us to detect ...a planetary mass companion, PDS 70 b, within the disk cavity. Moreover, observations in Hα with MagAO and MUSE revealed emission associated to PDS 70 b and to another new companion candidate, PDS 70 c, at a larger separation from the star. PDS 70 is the only multiple planetary system at its formation stage detected so far through direct imaging. Aims. Our aim is to confirm the discovery of the second planet PDS 70 c using SPHERE at VLT, to further characterize its physical properties, and search for additional point sources in this young planetary system. Methods. We re-analyzed archival SPHERE NIR observations and obtained new data in Y, J, H and K spectral bands for a total of four different epochs. The data were reduced using the data reduction and handling pipeline and the SPHERE data center. We then applied custom routines (e.g., ANDROMEDA and PACO) to subtract the starlight. Results. We re-detect both PDS 70 b and c and confirm that PDS 70 c is gravitationally bound to the star. We estimate this second planet to be less massive than 5 MJup and with a Teff around 900 K. Also, it has a low gravity with logg between 3.0 and 3.5 dex. In addition, a third object has been identified at short separation (~0.12′′) from the star and gravitationally bound to the star. Its spectrum is however very blue, meaning that we are probably seeing stellar light reflected by dust and our analysis seems to demonstrate that it is a feature of the inner disk. We cannot however completely exclude the possibility that it is a planetary mass object enshrouded by a dust envelope. In this latter case, its mass should be of the order of a few tens of M⊕. Moreover, we propose a possible structure for the planetary system based on our data, and find that this structure cannot be stable on a long timescale.
Context.
Clouds are ubiquitous in exoplanet atmospheres and they represent a challenge for the model interpretation of their spectra. When generating a large number of model spectra, complex cloud ...models often prove too costly numerically, whereas more efficient models may be overly simplified.
Aims.
We aim to constrain the atmospheric properties of the directly imaged planet HR 8799e with a free retrieval approach.
Methods.
We used our radiative transfer code petitRADTRANS for generating the spectra, which we coupled to the PyMultiNest tool. We added the effect of multiple scattering which is important for treating clouds. Two cloud model parameterizations are tested: the first incorporates the mixing and settling of condensates, the second simply parameterizes the functional form of the opacity.
Results.
In mock retrievals, using an inadequate cloud model may result in atmospheres that are more isothermal and less cloudy than the input. Applying our framework on observations of HR 8799e made with the GPI, SPHERE, and GRAVITY, we find a cloudy atmosphere governed by disequilibrium chemistry, confirming previous analyses. We retrieve that C/O = 0.60
−0.08
+0.07
. Other models have not yet produced a well constrained C/O value for this planet. The retrieved C/O values of both cloud models are consistent, while leading to different atmospheric structures: either cloudy or more isothermal and less cloudy. Fitting the observations with the self-consistent Exo-REM model leads to comparable results, without constraining C/O.
Conclusions.
With data from the most sensitive instruments, retrieval analyses of directly imaged planets are possible. The inferred C/O ratio of HR 8799e is independent of the cloud model and thus appears to be a robust. This C/O is consistent with stellar, which could indicate that the HR 8799e formed outside the CO
2
or CO iceline. As it is the innermost planet of the system, this constraint could apply to all HR 8799 planets.