Source author record

Vivien Parmentier

Vivien Parmentier appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

26 published item(s)

preprint2024arXiv

Modeling the day-night temperature variations of ultra-hot Jupiters: confronting non-grey general circulation models and observations

Ultra-hot Jupiters (UHJs) are natural laboratories to study extreme physics in planetary atmospheres and their rich observational data sets are yet to be confronted with models with varying complexities at a population level. In this work, we update the general circulation model of Tan & Komacek (2019) to include a non-grey radiative transfer scheme and apply it to simulate the realistic thermal structures, phase-dependent spectra, and wavelength-dependent phase curves of UHJs. We performed grids of models over a large range of equilibrium temperatures and rotation periods for varying assumptions, showing that the fractional day-night brightness temperature differences remain almost constant or slightly increase with increasing equilibrium temperature from the visible to mid-infrared wavelengths. This differs from previous work primarily due to the increasing planetary rotation rate with increasing equilibrium temperature for fixed host star type. Radiative effects of varying atmospheric compositions become more significant in dayside brightness temperature in longer wavelengths. Data-model comparisons of dayside brightness temperatures and phase curve amplitudes as a function of equilibrium temperature are in broad agreement. Observations show a large scatter compared to models even with a range of different assumptions, indicating significantly varying intrinsic properties in the hot Jupiter population. Our cloud-free models generally struggle to match all observations for individual targets with a single set of parameter choices, indicating the need for extra processes for understanding the heat transport of UHJs.

preprint2022arXiv

3D radiative-transfer for exoplanet atmospheres. gCMCRT: a GPU accelerated MCRT code

Radiative-transfer (RT) is a key component for investigating atmospheres of planetary bodies. With the 3D nature of exoplanet atmospheres being important in giving rise to their observable properties, accurate and fast 3D methods are required to be developed to meet future multi-dimensional and temporal data sets. We develop an open source GPU RT code, gCMCRT, a Monte Carlo RT forward model for general use in planetary atmosphere RT problems. We aim to automate the post-processing pipeline, starting from direct global circulation model (GCM) output to synthetic spectra. We develop albedo, emission and transmission spectra modes for 3D and 1D input structures. We include capability to use correlated-k and high-resolution opacity tables, the latter of which can be Doppler shifted inside the model. We post-process results from several GCM groups including ExoRad, SPARC/MITgcm THOR, UK Met Office UM, Exo-FMS and the Rauscher model. Users can therefore take advantage of desktop and HPC GPU computing solutions. gCMCRT is well suited for post-processing large GCM model grids produced by members of the community and for high-resolution 3D investigations.

preprint2022arXiv

A New Analysis of 8 Spitzer Phase Curves and Hot Jupiter Population Trends: Qatar-1b, Qatar-2b, WASP-52b, WASP-34b, and WASP-140b

With over 30 phase curves observed during the warm Spitzer mission, the complete data set provides a wealth of information relating to trends and three-dimensional properties of hot Jupiter atmospheres. In this work we present a comparative study of seven new Spitzer phase curves for four planets with equilibrium temperatures of T$_{eq}\sim$ 1300K: Qatar-2b, WASP-52b, WASP-34b, and WASP-140b, as well as the reanalysis of the 4.5 $\micron$ Qatar-1b phase curve due to the similar equilibrium temperature. In total, five 4.5 $\micron$ phase curves and three 3.6 $\micron$ phase curves are analyzed here with a uniform approach. Using these new results, in combination with literature values for the entire population of published Spitzer phase curves of hot Jupiters, we present evidence for a linear trend of increasing hot spot offset with increasing orbital period, as well as observational evidence for two classes of planets in apparent redistribution vs. equilibrium temperature parameter space, and tentative evidence for a dependence of hot spot offset on planetary surface gravity in our $\sim$ 1300 K sample. We do not find trends in apparent heat redistribution with orbital period or gravity. Non-uniformity in literature Spitzer data analysis techniques precludes a definitive determination of the sources or lack of trends.

preprint2022arXiv

Confirmation of Water Absorption in the Thermal Emission Spectrum of the Hot Jupiter WASP-77Ab with HST/WFC3

Secondary eclipse observations of hot Jupiters can reveal both their compositions and thermal structures. Previous observations have shown a diversity of hot Jupiter eclipse spectra, including absorption features, emission features, and featureless blackbody-like spectra. We present a secondary eclipse spectrum of the hot Jupiter WASP-77Ab observed between $1-5$ $μ$m with the Hubble Space Telescope (HST) and the Spitzer Space Telescope. The HST observations show signs of water absorption indicative of a non-inverted thermal structure. We fit the data with both a one-dimensional free retrieval and a grid of one-dimensional self-consistent forward models to confirm this non-inverted structure. The free retrieval places a $3σ$ lower limit on the atmospheric water abundance of $\log(n_\mathrm{H_2O})>-4.78$ and can not constrain the CO abundance. The grid fit produces a slightly super-stellar metallicity and constrains the carbon-to-oxygen ratio to less than or equal to the solar value. We also compare our data to recent high-resolution observations of WASP-77Ab taken with the Gemini-South/IGRINS spectrograph and find that our observations are consistent with the best-fit model to the high-resolution data. However, the metallicity derived from the IGRINS data is significantly lower than that derived from our self-consistent model fit. We find that this difference may be due to disequilibrium chemistry, and the varying results between the models applied here demonstrate the difficulty of constraining disequilibrium chemistry with low-resolution, low wavelength coverage data alone. Future work to combine observations from IGRINS, HST, and JWST will improve our estimate of the atmospheric composition of WASP-77Ab.

preprint2022arXiv

Giant Planets from the Inside-Out

Giant planets acquire gas, ices and rocks during the early formation stages of planetary systems and thus inform us on the formation process itself. Proceeding from inside out, examining the connections between the deep interiors and the observable atmospheres, linking detailed measurements on giant planets in the solar system to the wealth of data on brown dwarfs and giant exoplanets, we aim to provide global constraints on interiors structure and composition for models of the formation of these planets. New developments after the Juno and Cassini missions point to both Jupiter and Saturn having strong compositional gradients and stable regions from the atmosphere to the deep interior. This is also the case of Uranus and Neptune, based on available, limited data on these planets. Giant exoplanets and brown dwarfs provide us with new opportunities to link atmospheric abundances to bulk, interior abundances \rev{and to link these abundances and isotopic ratios to formation scenarios. Analysing the wealth of data becoming available} will require new models accounting for the complexity of the planetary interiors and atmospheres

preprint2022arXiv

Spatially-resolving the terminator: Variation of Fe, temperature and winds in WASP-76 b across planetary limbs and orbital phase

Exoplanet atmospheres are inherently three-dimensional systems in which thermal/chemical variation and winds can strongly influence spectra. Recently, the ultra-hot Jupiter WASP-76 b has shown evidence for condensation and asymmetric Fe absorption with time. However, it is currently unclear whether these asymmetries are driven by chemical or thermal differences between the two limbs, as precise constraints on variation in these have remained elusive due to the challenges of modelling these dynamics in a Bayesian framework. To address this we develop a new model, HyDRA-2D, capable of simultaneously retrieving morning and evening terminators with day-night winds. We explore variations in Fe, temperature profile, winds and opacity deck with limb and orbital phase using VLT/ESPRESSO observations of WASP-76 b. We find Fe is more prominent on the evening for the last quarter of the transit, with $\log(X_\mathrm{Fe}) = {-4.03}^{+0.28}_{-0.31}$, but the morning shows a lower abundance with a wider uncertainty, $\log(X_\mathrm{Fe}) = {-4.59}^{+0.85}_{-1.0}$, driven by degeneracy with the opacity deck and the stronger evening signal. We constrain 0.1 mbar temperatures ranging from $2950^{+111}_{-156}$ K to $2615^{+266}_{-275}$ K, with a trend of higher temperatures for the more irradiated atmospheric regions. We also constrain a day-night wind speed of $9.8^{+1.2}_{-1.1}$ km/s for the last quarter, higher than $5.9^{+1.5}_{-1.1}$ km/s for the first, in line with general circulation models. We find our new spatially- and phase-resolved treatment is statistically favoured by 4.9$σ$ over traditional 1D-retrievals, and thus demonstrate the power of such modelling for robust constraints with current and future facilities.

preprint2021arXiv

HST/WFC3 Complete Phase-resolved Spectroscopy of White Dwarf-Brown Dwarf Binaries WD 0137 and EPIC 2122

Brown dwarfs in close-in orbits around white dwarfs offer an excellent opportunity to investigate properties of fast-rotating, tidally-locked, and highly-irradiated atmospheres. We present Hubble Space Telescope Wide Field Camera 3 G141 phase-resolved observations of two brown dwarf-white dwarf binaries: WD 0137-349 and EPIC 212235321. Their 1.1 to 1.7 $μ$m phase curves demonstrate rotational modulations with semi-amplitudes of $5.27\pm0.02$% and $29.1\pm0.1$%; both can be well fit by multi-order Fourier series models. The high-order Fourier components have the same phase as the first order and are likely caused by hot spots located at the substellar points, suggesting inefficient day/night heat transfer. Both brown dwarfs' phase-resolved spectra can be accurately represented by linear combinations of their day- and night-side spectra. Fitting the irradiated brown dwarf model grids to the day-side spectra require a filling factor of ~50%, further supporting a hot spot dominating the emission of the day-sides. The night-side spectrum of WD 0137-349B is reasonably well fit by non-irradiated substellar models and the one of EPIC 212235321B can be approximated by a Planck function. We find strong spectral variations in the brown dwarfs' day/night flux and brightness temperature contrasts, which highlights the limitations of band-integrated measurements in probing heat transfer in irradiated objects. On the color-magnitude diagram, WD 0137-349B evolves along a cloudless model track connecting the early-L and mid-T spectral types, demonstrating that clouds and disequilibrium chemistry have a negligible effect on this object. A full interpretation of these high-quality phase-resolved spectra calls for new models that couple atmospheric circulation and radiative transfer under high-irradiation conditions.

preprint2021arXiv

Towards multi-dimensional analysis of transmission spectroscopy. Part II: Day-night induced biases in retrievals from hot to ultra-hot Jupiters

Hot Jupiters are very good targets for transmission spectroscopy analysis. Their atmospheres have a large scale height implying a high signal to noise ratio. As these planets orbit close to their stars, they often present strong thermal and chemical hetereogeneities between the day and the night side of their atmosphere. For the hottest ones, the thermal dissociation of several species occurs in their atmospheres which leads to a stronger chemical dichotomy between the two hemispheres. It has already been shown that the current retrieval algorithms, which are using 1D forward models, find biased molecular abundances in ultra hot Jupiters. Here, we quantify the effective temperature domain over which these biases are present. We use a set of 12 simulations of typical Hot Jupiters from Teq = 1000 K to Teq = 2100 K performed with the Substellar and Planetary Atmospheric Radiation and Circulation global climate model and generate transmission spectra that fully account for the 3D structure of the atmosphere with Pytmosph3R. These spectra are then analyzed using the 1D TauREx retrieval code. We find that for JWST-like data, accounting for non-isothermal vertical temperature profiles is required over the whole temperature range. We further find that 1D retrieval codes start to estimate wrong parameter values for planets with equilibrium temperatures greater than 1400 K if there are absorbers in the visible able to create a hot stratosphere. The high temperatures at low pressures indeed entail a thermal dissociation of species which creates a strong chemical day-night dichotomy. As a by-product, we demonstrate that when using synthetic observations to assess the detectability of a given feature or process using a Bayesian framework, it is valid to use non-randomized input data, as long as the anticipated observational uncertainties are correctly taken into account.

preprint2020arXiv

A transition between the hot and the ultra-hot Jupiter atmospheres

[Abridged] A key hypothesis in the field of exoplanet atmospheres is the trend of atmospheric thermal structure with planetary equilibrium temperature. We explore this trend and report here the first statistical detection of a transition in the near-infrared (NIR) atmospheric emission between hot and ultra-hot Jupiters. We measure this transition using secondary eclipse observations and interpret this phenomenon as changes in atmospheric properties, and more specifically in terms of transition from non-inverted to inverted thermal profiles. We examine a sample of 78 hot Jupiters with secondary eclipse measurements at 3.6 μm and 4.5 μm measured with Spitzer Infrared Array Camera (IRAC). We measure the deviation of the data from the blackbody, which we define as the difference between the observed 4.5 μm eclipse depth and that expected at this wavelength based on the brightness temperature measured at 3.6 μm. We study how the deviation between 3.6 and 4.5 μm changes with theoretical predictions with equilibrium temperature and incoming stellar irradiation. We reveal a clear transition in the observed emission spectra of the hot Jupiter population at 1660 +/- 100 K in the zero albedo, full redistribution equilibrium temperature. We find the hotter exoplanets have even hotter daysides at 4.5 μm compared to 3.6 μm, which manifests as an exponential increase in the emitted power of the planets with stellar insolation. We propose that the measured transition is a result of seeing carbon monoxide in emission due to the formation of temperature inversions in the atmospheres of the hottest planets. These thermal inversions could be caused by the presence of atomic and molecular species with high opacities in the optical and/or the lack of cooling species. We find that the population of hot Jupiters statistically disfavors high C/O planets (C/O>= 0.85).

preprint2020arXiv

Evidence for H$_{2}$ Dissociation and Recombination Heat Transport in the Atmosphere of KELT-9b

Phase curve observations provide an opportunity to study the full energy budgets of exoplanets by quantifying the amount of heat redistributed from their daysides to their nightsides. Theories explaining the properties of phase curves for hot Jupiters have focused on the balance between radiation and dynamics as the primary parameter controlling heat redistribution. However, recent phase curves have shown deviations from the trends that emerge from this theory, which has led to work on additional processes that may affect hot Jupiter energy budgets. One such process, molecular hydrogen dissociation and recombination, can enhance energy redistribution on ultra-hot Jupiters with temperatures above $\sim2000$ K. In order to study the impact of H$_{2}$ dissociation on ultra-hot Jupiters, we present a phase curve of KELT-9b observed with the Spitzer Space Telescope at 4.5 $μ$m. KELT-9b is the hottest known transiting planet, with a 4.5-$μ$m dayside brightness temperature of $4566^{+140}_{-136}$ K and a nightside temperature of $2556^{+101}_{-97}$ K. We observe a phase curve amplitude of $0.609 \pm 0.020$ and a hot spot offset of $18.7^{+2.1}_{-2.3}$ degrees. The observed amplitude is too small to be explained by a simple balance between radiation and advection. General circulation models (GCMs) and an energy balance model that include the effects of H$_{2}$ dissociation and recombination provide a better match to the data. The GCMs, however, predict a maximum hot spot offset of $5$ degrees, which disagrees with our observations at $>5σ$ confidence. This discrepancy may be due to magnetic effects in the planet's highly ionized atmosphere.

preprint2020arXiv

Global Chemistry and Thermal Structure Models for the Hot Jupiter WASP-43b and Predictions for JWST

The James Webb Space Telescope (JWST) is expected to revolutionize the field of exoplanets. The broad wavelength coverage and the high sensitivity of its instruments will allow characterization of exoplanetary atmospheres with unprecedented precision. Following the Call for the Cycle 1 Early Release Science Program, the Transiting Exoplanet Community was awarded time to observe several targets, including WASP-43b. The atmosphere of this hot Jupiter has been intensively observed but still harbors some mysteries, especially concerning the day-night temperature gradient, the efficiency of the atmospheric circulation, and the presence of nightside clouds. We will constrain these properties by observing a full orbit of the planet and extracting its spectroscopic phase curve in the 5--12 $μ$m range with JWST/MIRI. To prepare for these observations, we performed an extensive modeling work with various codes: radiative transfer, chemical kinetics, cloud microphysics, global circulation models, JWST simulators, and spectral retrieval. Our JWST simulations show that we should achieve a precision of 210 ppm per 0.1 $μ$m spectral bin on average, which will allow us to measure the variations of the spectrum in longitude and measure the night-side emission spectrum for the first time. If the atmosphere of WASP-43b is clear, our observations will permit us to determine if its atmosphere has an equilibrium or disequilibrium chemical composition, providing eventually the first conclusive evidence of chemical quenching in a hot Jupiter atmosphere. If the atmosphere is cloudy, a careful retrieval analysis will allow us to identify the cloud composition.

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

Understanding the atmospheric properties and chemical composition of the ultra-hot Jupiter HAT-P-7b II. Mapping the effects of gas kinetics

The atmospheres of ultra-hot Jupiters are commonly considered to be at thermochemical equilibrium. We aim to provide disequilibrium chemistry maps for a global understanding of the chemistry in HAT-P-7b's atmosphere and assess the importance of disequilibrium chemistry on UHJs. We apply a hierarchical modelling approach utilising 97 1D atmospheric profiles from 3D GCM of HAT-P-7b. For each 1D profile, we evaluate our kinetic cloud formation model consistently with the local gas-phase composition in chemical equilibrium. We then evaluate quenching results from a zeroth-order approximation in comparison to a kinetic gas-phase approach. We find that the zeroth-order approach of estimating quenching points agrees well with the full gas-kinetic modeling results. Chemical disequilibrium has the greatest effect on the nightside and morning abundance of species such as H, H$_2$O, CH$_4$, CO$_2$, HCN, and all C$_n$H$_m$ molecules; heavier C$_n$H$_m$ molecules are more affected by disequilibrium processes. CO abundance, however, is affected only marginally. While dayside abundances also notably change, those around the evening terminator of HAT-P-7b are the least affected by disequilibrium processes. The latter finding may partially explain the consistency of observed transmission spectra of UHJs with atmospheres in thermochemical equilibrium. Photochemistry only negligibly affects molecular abundances and quenching levels. In general, the quenching points of HAT-P-7b's atmosphere are at much lower pressures in comparison to the cooler hot-jupiters. We propose several avenues to look for the effect of disequilibrium processes on UHJs that are, in general, based on abundance and opacity measurements at different local times. It remains a challenge to completely disentangle this from the chemical effects of clouds and that of a primordial non-solar abundance.

preprint2019arXiv

Transit Signatures of Inhomogeneous Clouds on Hot Jupiters: Insights From Microphysical Cloud Modeling

We determine the observability in transmission of inhomogeneous cloud cover on the limbs of hot Jupiters through post processing a general circulation model to include cloud distributions computed using a cloud microphysics model. We find that both the east and west limb often form clouds, but that the different properties of these clouds enhances the limb to limb differences compared to the clear case. Using JWST it should be possible to detect the presence of cloud inhomogeneities by comparing the shape of the transit lightcurve at multiple wavelengths because inhomogeneous clouds impart a characteristic, wavelength dependent signature. This method is statistically robust even with limited wavelength coverage, uncertainty on limb darkening coefficients, and imprecise transit times. We predict that the short wavelength slope varies strongly with temperature. The hot limb of the hottest planets form higher altitude clouds composed of smaller particles leading to a strong rayleigh slope. The near infrared spectral features of clouds are almost always detectable, even when no spectral slope is visible in the optical. In some of our models a spectral window between 5 and 9 microns can be used to probe through the clouds and detect chemical spectral features. Our cloud particle size distributions are not log-normal and differ from species to species. Using the area or mass weighted particle size significantly alters the relative strength of the cloud spectral features compared to using the predicted size distribution. Finally, the cloud content of a given planet is sensitive to a species' desorption energy and contact angle, two parameters that could be constrained experimentally in the future.

preprint2016arXiv

The Impact of Non-Uniform Thermal Structure on the Interpretation of Exoplanet Emission Spectra

The determination of atmospheric structure and molecular abundances of planetary atmospheres via spectroscopy involves direct comparisons between models and data. While varying in sophistication, most model-spectra comparisons fundamentally assume "1D" model physics. However, knowledge from general circulation models and of solar system planets suggests that planetary atmospheres are inherently "3D" in their structure and composition. We explore the potential biases resulting from standard "1D" assumptions within a Bayesian atmospheric retrieval framework. Specifically, we show how the assumption of a single 1-dimensional thermal profile can bias our interpretation of the thermal emission spectrum of a hot Jupiter atmosphere that is composed of two thermal profiles. We retrieve upon spectra of unresolved model planets as observed with a combination of $HST$ WFC3+$Spitzer$ IRAC as well as $JWST$ under varying differences in the two thermal profiles. For WFC3+IRAC, there is a significantly biased estimate of CH$_4$ abundance using a 1D model when the contrast is 80%. For $JWST$, two thermal profiles are required to adequately interpret the data and estimate the abundances when contrast is greater than 40%. We also apply this preliminary concept to the recent WFC3+IRAC phase curve data of the hot Jupiter WASP-43b. We see similar behavior as present in our simulated data: while the H$_2$O abundance determination is robust, CH$_4$ is artificially well-constrained to incorrect values under the 1D assumption. Our work demonstrates the need to evaluate model assumptions in order to extract meaningful constraints from atmospheric spectra and motivates exploration of optimal observational setups.

preprint2016arXiv

The Influence of Non-Uniform Cloud Cover on Transit Transmission Spectra

We model the impact of non-uniform cloud cover on transit transmission spectra. Patchy clouds exist in nearly every solar system atmosphere, brown dwarfs, and transiting exoplanets. Our major findings suggest that fractional cloud coverage can exactly mimic high mean molecular weight atmospheres and vice-versa over certain wavelength regions, in particular, over the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) bandpass (1.1-1.7 $μ$m). We also find that patchy cloud coverage exhibits a signature that is different from uniform global clouds. Furthermore, we explain analytically why the "patchy cloud-high mean molecular weight" degeneracy exists. We also explore the degeneracy of non-uniform cloud coverage in atmospheric retrievals on both synthetic and real planets. We find from retrievals on a synthetic solar composition hot Jupiter with patchy clouds and a cloud free high mean molecular weight warm Neptune, that both cloud free high mean molecular weight atmospheres and partially cloudy atmospheres can explain the data equally well. Another key find is that the HST WFC3 transit transmission spectra of two well observed objects, the hot Jupiter HD189733b and the warm Neptune HAT-P-11b, can be explained well by solar composition atmospheres with patchy clouds without the need to invoke high mean molecular weight or global clouds. The degeneracy between high molecular weight and solar composition partially cloudy atmospheres can be broken by observing the molecular Rayleigh scattering differences between the two. Furthermore, the signature of partially cloudy limbs also appears as a $\sim$100 ppm residual in the ingress and egress of the transit light curves, provided the transit timing is known to seconds.

preprint2016arXiv

Transitions in the cloud composition of hot Jupiters

Over a large range of equilibrium temperatures, clouds shape the transmission spectrum of hot Jupiter atmospheres, yet their composition remains unknown. Recent observations show that the Kepler lightcurves of some hot Jupiters are asymmetric: for the hottest planets, the lightcurve peaks before secondary eclipse, whereas for planets cooler than $\sim1900\rm\,K$, it peaks after secondary eclipse. We use the thermal structure from 3D global circulation models to determine the expected cloud distribution and Kepler lightcurves of hot Jupiters. We demonstrate that the change from an optical lightcurve dominated by thermal emission to one dominated by scattering (reflection) naturally explains the observed trend from negative to positive offset. For the cool planets the presence of an asymmetry in the Kepler lightcurve is a telltale sign of the cloud composition, because each cloud species can produce an offset only over a narrow range of effective temperatures. By comparing our models and the observations, we show that the cloud composition of hot Jupiters likely varies with equilibrium temperature. We suggest that a transition occurs between silicate and manganese sulfide clouds at a temperature near $1600\rm\,K$, analogous to the L/T transition on brown dwarfs. The cold trapping of cloud species below the photosphere naturally produces such a transition and predicts similar transitions for other condensates, including TiO. We predict that most hot Jupiters should have cloudy nightsides, that partial cloudiness should be common at the limb and that the dayside hot spot should often be cloud-free.

preprint2014arXiv

A non-grey analytical model for irradiated atmospheres. II: Analytical vs. numerical solutions

The recent discovery and characterization of the diversity of the atmospheres of exoplanets and brown dwarfs calls for the development of fast and accurate analytical models. We quantify the accuracy of the analytical solution derived in paper I for an irradiated, non-grey atmosphere by comparing it to a state-of-the-art radiative transfer model. Then, using a grid of numerical models, we calibrate the different coefficients of our analytical model for irradiated solar-composition atmospheres of giant exoplanets and brown dwarfs. We show that the so-called Eddington approximation used to solve the angular dependency of the radiation field leads to relative errors of up to 5% on the temperature profile. We show that for realistic non-grey planetary atmospheres, the presence of a convective zone that extends to optical depths smaller than unity can lead to changes in the radiative temperature profile on the order of 20% or more. When the convective zone is located at deeper levels (such as for strongly irradiated hot Jupiters), its effect on the radiative atmosphere is smaller. We show that the temperature inversion induced by a strong absorber in the optical, such as TiO or VO is mainly due to non-grey thermal effects reducing the ability of the upper atmosphere to cool down rather than an enhanced absorption of the stellar light as previously thought. Finally, we provide a functional form for the coefficients of our analytical model for solar-composition giant exoplanets and brown dwarfs. This leads to fully analytical pressure-temperature profiles for irradiated atmospheres with a relative accuracy better than 10% for gravities between 2.5m/s^2 and 250 m/s^2 and effective temperatures between 100 K and 3000 K. This is a great improvement over the commonly used Eddington boundary condition.

preprint2014arXiv

Evolution of Exoplanets and their Parent Stars

Studying exoplanets with their parent stars is crucial to understand their population, formation and history. We review some of the key questions regarding their evolution with particular emphasis on giant gaseous exoplanets orbiting close to solar-type stars. For masses above that of Saturn, transiting exoplanets have large radii indicative of the presence of a massive hydrogen-helium envelope. Theoretical models show that this envelope progressively cools and contracts with a rate of energy loss inversely proportional to the planetary age. The combined measurement of planetary mass, radius and a constraint on the (stellar) age enables a global determination of the amount of heavy elements present in the planet interior. The comparison with stellar metallicity shows a correlation between the two, indicating that accretion played a crucial role in the formation of planets. The dynamical evolution of exoplanets also depends on the properties of the central star. We show that the lack of massive giant planets and brown dwarfs in close orbit around G-dwarfs and their presence around F-dwarfs are probably tied to the different properties of dissipation in the stellar interiors. Both the evolution and the composition of stars and planets are intimately linked.

preprint2014arXiv

Planetary population synthesis coupled with atmospheric escape: a statistical view of evaporation

We apply hydrodynamic evaporation models to different synthetic planet populations that were obtained from a planet formation code based on a core-accretion paradigm. We investigated the evolution of the planet populations using several evaporation models, which are distinguished by the driving force of the escape flow (X-ray or EUV), the heating efficiency in energy-limited evaporation regimes, or both. Although the mass distribution of the planet populations is barely affected by evaporation, the radius distribution clearly shows a break at approximately 2 $R_{\oplus}$. We find that evaporation can lead to a bimodal distribution of planetary sizes (Owen & Wu 2013) and to an "evaporation valley" running diagonally downwards in the orbital distance - planetary radius plane, separating bare cores from low-mass planet that have kept some primordial H/He. Furthermore, this bimodal distribution is related to the initial characteristics of the planetary populations because low-mass planetary cores can only accrete small primordial H/He envelopes and their envelope masses are proportional to their core masses. We also find that the population-wide effect of evaporation is not sensitive to the heating efficiency of energy-limited description. However, in two extreme cases, namely without evaporation or with a 100\% heating efficiency in an evaporation model, the final size distributions show significant differences; these two scenarios can be ruled out from the size distribution of $Kepler$ candidates.

preprint2014arXiv

Pseudo 2D chemical model of hot Jupiter atmospheres: application to HD 209458b and HD 189733b

We have developed a pseudo two-dimensional model of a planetary atmosphere, which takes into account thermochemical kinetics, photochemistry, vertical mixing, and horizontal transport, the latter being modeled as a uniform zonal wind. We have applied the model to the atmospheres of the hot Jupiters HD 209458b and HD 189733b. The adopted eddy diffusion coefficients are calculated by following the behaviour of passive tracers in three-dimensional general circulation models, which results in eddy values significantly below previous estimates. We find that the distribution of molecules with altitude and longitude in the atmospheres of these two hot Jupiters is complex because of the interplay of the various physical and chemical processes at work. Much of the distribution of molecules is driven by the strong zonal wind and the limited extent of vertical transport, resulting in an important homogenisation of the chemical composition with longitude. In general, molecular abundances are quenched horizontally to values typical of the hottest dayside regions, and thus the composition in the cooler nightside regions is highly contaminated by that of warmer dayside regions. As a consequence, the abundance of methane remains low, even below the predictions of previous one-dimensional models, which is likely to be in conflict with the high CH4 content inferred from observations of the dayside of HD 209458b. Another consequence of the important longitudinal homogenisation of the abundances is that the variability of the chemical composition has little effect on the way the emission spectrum is modified with phase and on the changes in the transmission spectrum from the transit ingress to the egress, these variations in the spectra being mainly due to changes in the temperature rather than in the composition between the different sides of the planet.

preprint2014arXiv

Unveiling the atmospheres of giant exoplanets with an EChO-class mission

More than a thousand exoplanets have been discovered over the last decade. Perhaps more excitingly, probing their atmospheres has become possible. With current data we have glimpsed the diversity of exoplanet atmospheres that will be revealed over the coming decade. However, numerous questions concerning their chemical composition, thermal structure, and atmospheric dynamics remain to be answered. More observations of higher quality are needed. In the next years, the selection of a space-based mission dedicated to the spectroscopic characterization of exoplanets would revolutionize our understanding of the physics of planetary atmospheres. Such a mission was proposed to the ESA cosmic vision program in 2014. Our paper is therefore based on the planned capabilities of the Exoplanet Characterization Observatory (EChO), but it should equally apply to any future mission with similar characteristics. With its large spectral coverage ($4-16\, \rm{μm}$), high spectral resolution ($Δλ/λ>300$ below $5\,\rm{μm}$ and $Δλ/λ>30$ above $5\,\rm{μm}$) and $1.5\rm{m}$ mirror, a future mission such as EChO will provide spectrally resolved transit lightcurves, secondary eclipses lightcurves, and full phase curves of numerous exoplanets with an unprecedented signal-to-noise ratio. In this paper, we review some of today's main scientific questions about gas giant exoplanets atmospheres, for which a future mission such as EChO will bring a decisive contribution.

preprint2013arXiv

3D mixing in hot Jupiter atmospheres. I. application to the day/night cold trap in HD 209458b

Hot Jupiters exhibit atmospheric temperatures ranging from hundreds to thousands of Kelvin. Because of their large day-night temperature differences, condensable species that are stable in the gas phase on the dayside, such as TiO and silicates, may condense and gravitationally settle on the nightside. Atmospheric circulation may counterbalance this tendency to gravitationally settle. This three dimensional (3D) mixing of chemical species has not previously been studied for hot Jupiters, yet it is crucial to assess the existence and distribution of TiO and silicates in the atmospheres of these planets. We perform 3D global circulation models of HD209458b including passive tracers that advect with the 3D flow, including a source/sink on the nightside to represent condensation and gravitational settling of haze particles. We show that global advection patterns produce strong vertical mixing that can keep condensable species lofted as long as they are trapped in particles of sizes of a few microns or less on the night side. We show that vertical mixing results not from small-scale convection but from the large-scale circulation driven by the day-night heating contrast. Although this vertical mixing is not diffusive in any rigorous sense, a comparison of our results with idealized diffusion models allows a rough estimate of the vertical diffusion coefficient. Kzz=5x10**4/Sqrt(Pbar) m2/s can be used in 1D models of HD 209458b. Moreover, our models exhibit strong spatial and temporal variability in the tracer concentration that could result in observable variations during transit or secondary eclipse measurements. Finally, we apply our model to the case of TiO in HD209458b and show that the day-night cold trap would deplete TiO if it condenses into particles bigger than a few microns on the planet's night side, making it unable to create the observed stratosphere of the planet.

preprint2013arXiv

A non-grey analytical model for irradiated atmospheres. I: Derivation

Context. Semi-grey atmospheric models (with one opacity for the visible and one opacity for the infrared) are useful to understand the global structure of irradiated atmospheres, their dynamics and the interior structure and evolution of planets, brown dwarfs and stars. But when compared to direct numerical radiative transfer calculations for irradiated exoplanets, these models systematically overestimate the temperatures at low optical depth, independently of the opacity parameters. We wish to understand why semi-grey models fail at low optical depths, and provide a more accurate approximation to the atmospheric structure by accounting for the variable opacity in the infrared. Our analytical irradiated non-grey model is found to provide a range of temperatures that is consistent with that obtained by numerical calculations. We find that even for slightly non-grey thermal opacities the temperature structure differs significantly from previous semi-grey models. For small values of beta (expected when lines are dominant), we find that the non-grey effects are confined to low-optical depths. However, for beta larger than 0.5 (appropriate in the presence of bands with a wavelength-dependence smaller or comparable with the width of the Planck function), we find that the temperature structure is affected even down to infrared optical depths unity and deeper as a result of the so-called blanketing effect. The expressions that we derive may be used to provide a proper functional form for algorithms that invert the atmospheric properties from spectral information. Because a full atmospheric structure can be calculated directly, these expressions should be useful for simulations of the dynamics of these atmospheres and of the thermal evolution of the planets. Finally, they should be used to test full radiative transfer models and improve their convergence.

preprint2013arXiv

Bulk Composition of GJ 1214b and other sub-Neptune exoplanets

GJ1214b stands out among the detected low-mass exoplanets, because it is, so far, the only one amenable to transmission spectroscopy. Up to date there is no consensus about the composition of its envelope although most studies suggest a high molecular weight atmosphere. In particular, it is unclear if hydrogen and helium are present or if the atmosphere is water dominated. Here, we present results on the composition of the envelope obtained by using an internal structure and evolutionary model to fit the mass and radius data. By examining all possible mixtures of water and H/He, with the corresponding opacities, we find that the bulk amount of H/He of GJ1214b is at most 7% by mass. In general, we find the radius of warm sub-Neptunes to be most sensitive to the amount of H/He. We note that all (Kepler-11b,c,d,f, Kepler-18b, Kepler-20b, 55Cnc-e, Kepler-36c and Kepler-68b) but two (Kepler-11e and Kepler-30b) of the discovered low-mass planets so far have less than 10% H/He. In fact, Kepler-11e and Kepler-30b have 10-18% and 5-15% bulk H/He. Conversely, little can be determined about the H2O or rocky content of sub-Neptune planets. We find that although a 100% water composition fits the data for GJ1214b, based on formation constraints the presence of heavier refractory material on this planet is expected, and hence, so is a component lighter than water required. A robust determination by transmission spectroscopy of the composition of the upper atmosphere of GJ1214b will help determine the extent of compositional segregation between the atmosphere and envelope.

preprint2013arXiv

Inference of Inhomogeneous Clouds in an Exoplanet Atmosphere

We present new visible and infrared observations of the hot Jupiter Kepler-7b to determine its atmospheric properties. Our analysis allows us to 1) refine Kepler-7b's relatively large geometric albedo of Ag=0.35+-0.02, 2) place upper limits on Kepler-7b thermal emission that remains undetected in both Spitzer bandpasses and 3) report a westward shift in the Kepler optical phase curve. We argue that Kepler-7b's visible flux cannot be due to thermal emission or Rayleigh scattering from H2 molecules. We therefore conclude that high altitude, optically reflective clouds located west from the substellar point are present in its atmosphere. We find that a silicate-based cloud composition is a possible candidate. Kepler-7b exhibits several properties that may make it particularly amenable to cloud formation in its upper atmosphere. These include a hot deep atmosphere that avoids a cloud cold trap, very low surface gravity to suppress cloud sedimentation, and a planetary equilibrium temperature in a range that allows for silicate clouds to potentially form in the visible atmosphere probed by Kepler. Our analysis does not only present evidence of optically thick clouds on Kepler-7b but also yields the first map of clouds in an exoplanet atmosphere.