-
Temperature-chemistry coupling in the evolution of gas giant atmospheres driven by stellar flares
Authors:
Harrison Nicholls,
Eric Hébrard,
Olivia Venot,
Benjamin Drummond,
Elise Evans
Abstract:
The effect of enhanced UV irradiation associated with stellar flares on the atmospheric composition and temperature of gas giant exoplanets was investigated. This was done using a 1D radiative-convective-chemical model with self-consistent feedback between the temperature and the non-equilibrium chemistry.
It was found that flare-driven changes to chemical composition and temperature give rise t…
▽ More
The effect of enhanced UV irradiation associated with stellar flares on the atmospheric composition and temperature of gas giant exoplanets was investigated. This was done using a 1D radiative-convective-chemical model with self-consistent feedback between the temperature and the non-equilibrium chemistry.
It was found that flare-driven changes to chemical composition and temperature give rise to prolonged trends in evolution across a broad range of pressure levels and species. Allowing feedback between chemistry and temperature plays an important role in establishing the quiescent structure of these atmospheres, and determines their evolution due to flares. It was found that cooler planets are more susceptible to flares than warmer ones, seeing larger changes in composition and temperature, and that temperature-chemistry feedback modifies their evolution.
Long-term exposure to flares changes the transmission spectra of gas giant atmospheres; these changes differed when the temperature structure was allowed to evolve self-consistently with the chemistry. Changes in spectral features due to the effects of flares on these atmospheres can be associated with changes in composition. The effects of flares on the atmospheres of sufficiently cool planets will impact observations made with JWST. It is necessary to use self-consistent models of temperature and chemistry in order to accurately capture the effects of flares on features in the transmission spectra of cooler gas giants, but this depends heavily on the radiation environment of the planet.
△ Less
Submitted 6 June, 2023;
originally announced June 2023.
-
A modern-day Mars climate in the Met Office Unified Model: dry simulations
Authors:
Danny McCulloch,
Denis E. Sergeev,
Nathan Mayne,
Matthew Bate,
James Manners,
Ian Boutle,
Benjamin Drummond,
Kristzian Kohary
Abstract:
We present results from the Met Office Unified Model (UM), a world-leading climate and weather model, adapted to simulate a dry Martian climate. We detail the adaptation of the basic parameterisations and analyse results from two simulations, one with radiatively active mineral dust and one with radiatively inactive dust. These simulations demonstrate how the radiative effects of dust act to accel…
▽ More
We present results from the Met Office Unified Model (UM), a world-leading climate and weather model, adapted to simulate a dry Martian climate. We detail the adaptation of the basic parameterisations and analyse results from two simulations, one with radiatively active mineral dust and one with radiatively inactive dust. These simulations demonstrate how the radiative effects of dust act to accelerate the winds and create a mid-altitude isothermal layer during the dusty season. We validate our model through comparison with an established Mars model, the Laboratoire de Météorologie Dynamique planetary climate model (PCM), finding good agreement in the seasonal wind and temperature profiles but with discrepancies in the predicted dust mass mixing ratio and conditions at the poles. This study validates the use of the UM for a Martian atmosphere, highlighting how the adaptation of an Earth general circulation model (GCM) can be beneficial for existing Mars GCMs and provides insight into the next steps in our development of a new Mars climate model.
△ Less
Submitted 1 February, 2023;
originally announced February 2023.
-
Observability of signatures of transport-induced chemistry in clear atmospheres of hot gas giant exoplanets
Authors:
Maria Zamyatina,
Éric Hébrard,
Benjamin Drummond,
Nathan J. Mayne,
James Manners,
Duncan A. Christie,
Pascal Tremblin,
David K. Sing,
Krisztian Kohary
Abstract:
Transport-induced quenching, i.e., the homogenisation of chemical abundances by atmospheric advection, is thought to occur in the atmospheres of hot gas giant exoplanets. While some numerical modelling of this process exists, the three-dimensional nature of transport-induced chemistry is underexplored. Here we present results of 3D cloud- and haze-free simulations of the atmospheres of HAT-P-11b,…
▽ More
Transport-induced quenching, i.e., the homogenisation of chemical abundances by atmospheric advection, is thought to occur in the atmospheres of hot gas giant exoplanets. While some numerical modelling of this process exists, the three-dimensional nature of transport-induced chemistry is underexplored. Here we present results of 3D cloud- and haze-free simulations of the atmospheres of HAT-P-11b, HD 189733b, HD 209458b, and WASP-17b including coupled hydrodynamics, radiative transfer and chemistry. Our simulations were performed with two chemical schemes: a chemical kinetics scheme, which is capable of capturing transport-induced quenching, and a simpler, more widely used chemical equilibrium scheme. We find that transport-induced quenching is predicted to occur in atmospheres of all planets in our sample; however, the extent to which it affects their synthetic spectra and phase curves varies from planet to planet. This implies that there is a "sweet spot" for the observability of signatures of transport-induced quenching, which is controlled by the interplay between the dynamics and chemistry.
△ Less
Submitted 21 November, 2022; v1 submitted 16 November, 2022;
originally announced November 2022.
-
3D modelling of the impact of stellar activity on tidally locked terrestrial exoplanets: atmospheric composition and habitability
Authors:
Robert J. Ridgway,
Maria Zamyatina,
Nathan J. Mayne,
James Manners,
F. Hugo Lambert,
Marrick Braam,
Benjamin Drummond,
Éric Hébrard,
Paul I. Palmer,
Krisztian Kohary
Abstract:
Stellar flares present challenges to the potential habitability of terrestrial planets orbiting M dwarf stars through inducing changes in the atmospheric composition and irradiating the planet's surface in large amounts of ultraviolet light. To examine their impact, we have coupled a general circulation model with a photochemical kinetics scheme to examine the response and changes of an Earth-like…
▽ More
Stellar flares present challenges to the potential habitability of terrestrial planets orbiting M dwarf stars through inducing changes in the atmospheric composition and irradiating the planet's surface in large amounts of ultraviolet light. To examine their impact, we have coupled a general circulation model with a photochemical kinetics scheme to examine the response and changes of an Earth-like atmosphere to stellar flares and coronal mass ejections. We find that stellar flares increase the amount of ozone in the atmosphere by a factor of 20 compared to a quiescent star. We find that coronal mass ejections abiotically generate significant levels of potential bio-signatures such as N$_2$O. The changes in atmospheric composition cause a moderate decrease in the amount of ultraviolet light that reaches the planets surface, suggesting that while flares are potentially harmful to life, the changes in the atmosphere due to a stellar flare act to reduce the impact of the next stellar flare.
△ Less
Submitted 25 October, 2022; v1 submitted 24 October, 2022;
originally announced October 2022.
-
The Impact of Mixing Treatments on Cloud Modelling in 3D Simulations of Hot Jupiters
Authors:
D. A. Christie,
N. J. Mayne,
S. Lines,
V. Parmentier,
J. Manners,
I. Boutle,
B. Drummond,
T. Mikal-Evans,
D. K. Sing,
K. Kohary
Abstract:
We present results of 3D hydrodynamical simulations of HD209458b including a coupled, radiatively-active cloud model ({\sc EddySed}). We investigate the role of the mixing by replacing the default convective treatment used in previous works with a more physically relevant mixing treatment ($K_{zz}$) based on global circulation. We find that uncertainty in the efficiency of sedimentation through th…
▽ More
We present results of 3D hydrodynamical simulations of HD209458b including a coupled, radiatively-active cloud model ({\sc EddySed}). We investigate the role of the mixing by replacing the default convective treatment used in previous works with a more physically relevant mixing treatment ($K_{zz}$) based on global circulation. We find that uncertainty in the efficiency of sedimentation through the sedimentation factor $f_\mathrm{sed}$ plays a larger role in shaping cloud thickness and its radiative feedback on the local gas temperatures -- e.g. hot spot shift and day-to-night side temperature gradient -- than the switch in mixing treatment. We demonstrate using our new mixing treatments that simulations with cloud scales which are a fraction of the pressure scale height improve agreement with the observed transmission spectra, the emission spectra, and the Spitzer 4.5 $\mathrm{μm}$ phase curve, although our models are still unable to reproduce the optical and UV transmission spectra. We also find that the inclusion of cloud increases the transit asymmetry in the optical between the east and west limbs, although the difference remains small ($\lesssim 1\%$).
△ Less
Submitted 12 July, 2021;
originally announced July 2021.
-
Ground-Based Transmission Spectroscopy with VLT FORS2: Evidence for faculae and clouds in the optical spectrum of the warm Saturn WASP-110b
Authors:
Nikolay Nikolov,
Gracjan Maciejewski,
Savvas Constantinou,
Nikku Madhusudhan,
Jonathan J. Fortney,
Barry Smalley,
Aarynn L. Carter,
Ernst J. W. de Mooij,
Benjamin Drummond,
Neale P. Gibson,
Christiane Helling,
Nathan J. Mayne,
Thomas Mikal-Evans,
David K. Sing,
Jamie Wilson
Abstract:
We present a ground-based optical transmission spectrum for the warm Saturn-mass exoplanet WASP-110b from two transit observations made with the FOcal Reducer and Spectrograph (FORS2) on the Very Large Telescope (VLT). The spectrum covers the wavelength range from 4000 to 8333Å, which is binned in 46 transit depths measured to an averaged precision of 220 parts per million (ppm) over an averaged 8…
▽ More
We present a ground-based optical transmission spectrum for the warm Saturn-mass exoplanet WASP-110b from two transit observations made with the FOcal Reducer and Spectrograph (FORS2) on the Very Large Telescope (VLT). The spectrum covers the wavelength range from 4000 to 8333Å, which is binned in 46 transit depths measured to an averaged precision of 220 parts per million (ppm) over an averaged 80Å~bin for a Vmag=12.8 star. The measured transit depths are unaffected by a dilution from a close A-type field dwarf, which was fully resolved. The overall main characteristic of the transmission spectrum is an increasing radius with wavelength and a lack of the theoretically predicted pressure-broadened sodium and potassium absorption features for a cloud-free atmosphere. We analyze archival high-resolution optical spectroscopy and find evidence for low to moderate activity of the host star, which we take into account in the atmospheric retrieval analysis. Using the AURA retrieval code, we find that the observed transmission spectrum can be best explained by a combination of unocculted stellar faculae and a cloud deck. Transmission spectra of cloud-free and hazy atmospheres are rejected at a high confidence. With a possible cloud deck at its terminator, WASP-110b joins the increasing population of irradiated hot-Jupiter exoplanets with cloudy atmospheres observed in transmission.
△ Less
Submitted 13 May, 2021;
originally announced May 2021.
-
Pseudo-2D Modelling of Heat Redistribution Through H$_2$ Thermal Dissociation/Recombination: Consequences for Ultra-Hot Jupiters
Authors:
Alexander Roth,
Benjamin Drummond,
Eric Hébrard,
Pascal Tremblin,
Jayesh Goyal,
Nathan Mayne
Abstract:
Thermal dissociation and recombination of molecular hydrogen, H_2, in the atmospheres of ultra-hot Jupiters (UHJs) has been shown to play an important role in global heat redistribution. This, in turn, significantly impacts their planetary emission, yet only limited investigations on the atmospheric effects have so far been conducted. Here we investigate the heat redistribution caused by this diss…
▽ More
Thermal dissociation and recombination of molecular hydrogen, H_2, in the atmospheres of ultra-hot Jupiters (UHJs) has been shown to play an important role in global heat redistribution. This, in turn, significantly impacts their planetary emission, yet only limited investigations on the atmospheric effects have so far been conducted. Here we investigate the heat redistribution caused by this dissociation/recombination reaction, alongside feedback mechanisms between the atmospheric chemistry and radiative transfer, for a planetary and stellar configuration typical of UHJs. To do this, we have developed a time-dependent pseudo-2D model, including a treatment of time-independent equilibrium chemical effects. As a result of the reaction heat redistribution, we find temperature changes of up to $\sim$400 K in the atmosphere. When TiO and VO are additionally considered as opacity sources, these changes in temperature increase to over $\sim$800 K in some areas. This heat redistribution is found to significantly shift the region of peak atmospheric temperature, or hotspot, towards the evening terminator in both cases. The impact of varying the longitudinal wind speed on the reaction heat distribution is also investigated. When excluding TiO/VO, increased wind speeds are shown to increase the impact of the reaction heat redistribution up to a threshold wind speed. When including TiO/VO there is no apparent wind speed threshold, due to thermal stabilisation by these species. We also construct pseudo-2D phase curves from our model, and highlight both significant spectral flux damping and increased phase offset caused by the reaction heat redistribution.
△ Less
Submitted 29 April, 2021;
originally announced April 2021.
-
A Library of Self-Consistent Simulated Exoplanet Atmospheres
Authors:
Jayesh M. Goyal,
Nathan Mayne,
Benjamin Drummond,
David K. Sing,
Eric Hébrard,
Nikole Lewis,
Pascal Tremblin,
Mark W. Phillips,
Thomas Mikal-Evans,
Hannah R. Wakeford
Abstract:
We present a publicly available library of model atmospheres with radiative-convective equilibrium Pressure-Temperature ($P$-$T$) profiles fully consistent with equilibrium chemical abundances, and the corresponding emission and transmission spectrum with R$\sim$5000 at 0.2 $μ$m decreasing to R$\sim$35 at 30 $μ$m, for 89 hot Jupiter exoplanets, for four re-circulation factors, six metallicities an…
▽ More
We present a publicly available library of model atmospheres with radiative-convective equilibrium Pressure-Temperature ($P$-$T$) profiles fully consistent with equilibrium chemical abundances, and the corresponding emission and transmission spectrum with R$\sim$5000 at 0.2 $μ$m decreasing to R$\sim$35 at 30 $μ$m, for 89 hot Jupiter exoplanets, for four re-circulation factors, six metallicities and six C/O ratios. We find the choice of condensation process (local/rainout) alters the $P$-$T$ profile and thereby the spectrum substantially, potentially detectable by JWST. We find H$^-$ opacity can contribute to form a strong temperature inversion in ultra-hot Jupiters for C/O ratios $\geq$ 1 and can make transmission spectra features flat in the optical, alongside altering the entire emission spectra. We highlight how adopting different model choices such as thermal ionisation, opacities, line-wing profiles and the methodology of varying the C/O ratio, effects the $P$-$T$ structure and the spectrum. We show the role of Fe opacity to form primary/secondary inversion in the atmosphere. We use WASP-17b and WASP-121b as test cases to demonstrate the effect of grid parameters across their full range, while highlighting some important findings, concerning the overall atmospheric structure, chemical transition regimes and their observables. Finally, we apply this library to the current transmission and emission spectra observations of WASP-121b, which shows H$_2$O and tentative evidence for VO at the limb, and H$_2$O emission feature indicative of inversion on the dayside, with very low energy redistribution, thereby demonstrating the applicability of library for planning and interpreting observations of transmission and emission spectrum.
△ Less
Submitted 4 August, 2020;
originally announced August 2020.
-
Ground-Based Transmission Spectroscopy with FORS2: A featureless optical transmission spectrum and detection of H$_2$O for the ultra-hot Jupiter WASP-103b
Authors:
J. Wilson,
N. P. Gibson,
N. Nikolov,
S. Constantinou,
N. Madhusudhan,
J. Goyal,
J. K. Barstow,
A. L. Carter,
E. J. W. de Mooij,
B. Drummond,
T. Mikal-Evans,
C. Helling,
N. J. Mayne,
D. K. Sing
Abstract:
We report ground-based transmission spectroscopy of the highly irradiated and ultra-short period hot-Jupiter WASP-103b covering the wavelength range $\approx$ 400-600 nm using the FORS2 instrument on the Very Large Telescope. The light curves show significant time-correlated noise which is mainly invariant in wavelength and which we model using a Gaussian process. The precision of our transmission…
▽ More
We report ground-based transmission spectroscopy of the highly irradiated and ultra-short period hot-Jupiter WASP-103b covering the wavelength range $\approx$ 400-600 nm using the FORS2 instrument on the Very Large Telescope. The light curves show significant time-correlated noise which is mainly invariant in wavelength and which we model using a Gaussian process. The precision of our transmission spectrum is improved by applying a common-mode correction derived from the white light curve, reaching typical uncertainties in transit depth of $\approx$ 2x10$^{-4}$ in wavelength bins of 15 nm. After correction for flux contamination from a blended companion star, our observations reveal a featureless spectrum across the full range of the FORS2 observations and we are unable to confirm the Na absorption previously inferred using Gemini/GMOS or the strong Rayleigh scattering observed using broad-band light curves. We performed a Bayesian atmospheric retrieval on the full optical-infrared transmission spectrum using the additional data from Gemini/GMOS, HST/WFC3 and Spitzer observations and recover evidence for H$_2$O absorption at the 4.0$σ$ level. However, our observations are not able to completely rule out the presence of Na, which is found at 2.0$σ$ in our retrievals. This may in part be explained by patchy/inhomogeneous clouds or hazes damping any absorption features in our FORS2 spectrum, but an inherently small scale height also makes this feature challenging to probe from the ground. Our results nonetheless demonstrate the continuing potential of ground-based observations for investigating exoplanet atmospheres and emphasise the need for the application of consistent and robust statistical techniques to low-resolution spectra in the presence of instrumental systematics.
△ Less
Submitted 27 July, 2020;
originally announced July 2020.
-
A new set of atmosphere and evolution models for cool T-Y brown dwarfs and giant exoplanets
Authors:
Mark W. Phillips,
Pascal Tremblin,
Isabelle Baraffe,
Gilles Chabrier,
Nicole F. Allard,
Fernand Spiegelman,
Jayesh M. Goyal,
Ben Drummond,
Eric Hebrard
Abstract:
We present a new set of solar metallicity atmosphere and evolutionary models for very cool brown dwarfs and self-luminous giant exoplanets, which we term ATMO 2020. Atmosphere models are generated with our state-of-the-art 1D radiative-convective equilibrium code ATMO, and are used as surface boundary conditions to calculate the interior structure and evolution of…
▽ More
We present a new set of solar metallicity atmosphere and evolutionary models for very cool brown dwarfs and self-luminous giant exoplanets, which we term ATMO 2020. Atmosphere models are generated with our state-of-the-art 1D radiative-convective equilibrium code ATMO, and are used as surface boundary conditions to calculate the interior structure and evolution of $0.001-0.075\,\mathrm{M_{\odot}}$ objects. Our models include several key improvements to the input physics used in previous models available in the literature. Most notably, the use of a new H-He equation of state including ab initio quantum molecular dynamics calculations has raised the mass by $\sim1-2\%$ at the stellar-substellar boundary and has altered the cooling tracks around the hydrogen and deuterium burning minimum masses. A second key improvement concerns updated molecular opacities in our atmosphere model ATMO, which now contains significantly more line transitions required to accurately capture the opacity in these hot atmospheres. This leads to warmer atmospheric temperature structures, further changing the cooling curves and predicted emission spectra of substellar objects. We present significant improvement for the treatment of the collisionally broadened potassium resonance doublet, and highlight the importance of these lines in shaping the red-optical and near-infrared spectrum of brown dwarfs. We generate three different grids of model simulations, one using equilibrium chemistry and two using non-equilibrium chemistry due to vertical mixing, all three computed self-consistently with the pressure-temperature structure of the atmosphere. We show the impact of vertical mixing on emission spectra and in colour-magnitude diagrams, highlighting how the $3.5-5.5\,\mathrm{μm}$ flux window can be used to calibrate vertical mixing in cool T-Y spectral type objects.
△ Less
Submitted 30 March, 2020;
originally announced March 2020.
-
Implications of three-dimensional chemical transport in hot Jupiter atmospheres: results from a consistently coupled chemistry-radiation-hydrodynamics model
Authors:
Benjamin Drummond,
Eric Hebrard,
Nathan J. Mayne,
Olivia Venot,
Robert J. Ridgway,
Quentin Changeat,
Shang-min Tsai,
James Manners,
Pascal Tremblin,
Nathan Luke Abraham,
David Sing,
Krisztian Kohary
Abstract:
We present results from a set of simulations using a fully coupled three-dimensional (3D) chemistry-radiation-hydrodynamics model and investigate the effect of transport of chemical species by the large-scale atmospheric flow in hot Jupiter atmospheres. We couple a flexible chemical kinetics scheme to the Met Office Unified Model which enables the study of the interaction of chemistry, radiative t…
▽ More
We present results from a set of simulations using a fully coupled three-dimensional (3D) chemistry-radiation-hydrodynamics model and investigate the effect of transport of chemical species by the large-scale atmospheric flow in hot Jupiter atmospheres. We couple a flexible chemical kinetics scheme to the Met Office Unified Model which enables the study of the interaction of chemistry, radiative transfer and fluid dynamics. We use a newly-released "reduced" chemical network comprising 30 chemical species that has been specifically developed for application in 3D atmosphere models. We simulate the atmospheres of the well-studied hot Jupiters HD~209458b and HD~189733b which both have dayside--nightside temperature contrasts of several hundred Kelvin and superrotating equatorial jets. We find qualitatively quite different chemical structures between the two planets, particularly for methane (CH$_4$), when advection of chemical species is included. Our results show that consideration of 3D chemical transport is vital in understanding the chemical composition of hot Jupiter atmospheres. 3D mixing leads to significant changes in the abundances of absorbing gas-phase species compared with what would be expected by assuming local chemical equilibrium, or from models including 1D - and even 2D - chemical mixing. We find that CH$_4$, carbon dioxide (CO$_2$) and ammonia (NH$_3$) are particularly interesting as 3D mixing of these species leads to prominent signatures of out-of-equilibrium chemistry in the transmission and emission spectra, detectable with near-future instruments.
△ Less
Submitted 14 April, 2020; v1 submitted 30 January, 2020;
originally announced January 2020.
-
Idealised simulations of the deep atmosphere of hot jupiters: Deep, hot, adiabats as a robust solution to the radius inflation problem
Authors:
F. Sainsbury-Martinez,
P. Wang,
S. Fromang,
P. Tremblin,
T. Dubos,
Y. Meurdesoif,
A. Spiga,
J. Leconte,
I. Baraffe,
G. Chabrier,
N. Mayne,
B. Drummond,
F. Debras
Abstract:
Context: The anomalously large radii of hot Jupiters has long been a mystery. However, by combining both theoretical arguments and 2D models, a recent study has suggested that the vertical advection of potential temperature leads to an adiabatic temperature profile in the deep atmosphere hotter than the profile obtained with standard 1D models. Aims: In order to confirm the viability of that scena…
▽ More
Context: The anomalously large radii of hot Jupiters has long been a mystery. However, by combining both theoretical arguments and 2D models, a recent study has suggested that the vertical advection of potential temperature leads to an adiabatic temperature profile in the deep atmosphere hotter than the profile obtained with standard 1D models. Aims: In order to confirm the viability of that scenario, we extend this investigation to three dimensional, time-dependent, models. Methods: We use a 3D GCM, DYNAMICO to perform a series of calculations designed to explore the formation and structure of the driving atmospheric circulations, and detail how it responds to changes in both upper and deep atmospheric forcing. Results: In agreement with the previous, 2D, study, we find that a hot adiabat is the natural outcome of the long-term evolution of the deep atmosphere. Integration times of order $1500$ years are needed for that adiabat to emerge from an isothermal atmosphere, explaining why it has not been found in previous hot Jupiter studies. Models initialised from a hotter deep atmosphere tend to evolve faster toward the same final state. We also find that the deep adiabat is stable against low-levels of deep heating and cooling, as long as the Newtonian cooling time-scale is longer than $\sim 3000$ years at $200$ bar. Conclusions: We conclude that the steady-state vertical advection of potential temperature by deep atmospheric circulations constitutes a robust mechanism to explain hot Jupiter inflated radii. We suggest that future studies of hot Jupiters are evolved for a longer time than currently done, and, when possible, include models initialised with a hot deep adiabat. We stress that this mechanism stems from the advection of entropy by irradiation induced mass flows and does not require (finely tuned) dissipative process, in contrast with most previously suggested scenarios.
△ Less
Submitted 15 November, 2019;
originally announced November 2019.
-
Overcast on Osiris: 3D radiative-hydrodynamical simulations of a cloudy hot Jupiter using the parameterised, phase-equilibrium cloud formation code EddySed
Authors:
S. Lines,
N. J. Mayne,
J. Manners,
I. A. Boutle,
B. Drummond,
T. Mikal-Evans,
K. Kohary,
D. K. Sing
Abstract:
We present results from 3D radiative-hydrodynamical simulations of HD 209458b with a fully coupled treatment of clouds using the EddySed code, critically, including cloud radiative feedback via absorption and scattering. We demonstrate that the thermal and optical structure of the simulated atmosphere is markedly different, for the majority of our simulations, when including cloud radiative effect…
▽ More
We present results from 3D radiative-hydrodynamical simulations of HD 209458b with a fully coupled treatment of clouds using the EddySed code, critically, including cloud radiative feedback via absorption and scattering. We demonstrate that the thermal and optical structure of the simulated atmosphere is markedly different, for the majority of our simulations, when including cloud radiative effects, suggesting this important mechanism can not be neglected. Additionally, we further demonstrate that the cloud structure is sensitive to not only the cloud sedimentation efficiency (termed $f_{\textrm{sed}}$ in EddySed), but also the temperature-pressure profile of the deeper atmosphere. We briefly discuss the large difference between the resolved cloud structures of this work, adopting a phase-equilibrium and parameterised cloud model, and our previous work incorporating a cloud microphysical model, although a fairer comparison where, for example, the same list of constituent condensates is included in both treatments, is reserved for a future work. Our results underline the importance of further study into the potential condensate size distributions and vertical structures, as both strongly influence the radiative impact of clouds on the atmosphere. Finally, we present synthetic observations from our simulations reporting an improved match, over our previous cloud-free simulations, to the observed transmission, HST WFC3 emission and 4.5 $μ$m Spitzer phase curve of HD 209458b. Additionally, we find all our cloudy simulations have an apparent albedo consistent with observations.
△ Less
Submitted 27 June, 2019;
originally announced June 2019.
-
An emission spectrum for WASP-121b measured across the 0.8-1.1 micron wavelength range using the Hubble Space Telescope
Authors:
Thomas Mikal-Evans,
David K. Sing,
Jayesh M. Goyal,
Benjamin Drummond,
Aarynn L. Carter,
Gregory W. Henry,
Hannah R. Wakeford,
Nikole K. Lewis,
Mark S. Marley,
Pascal Tremblin,
Nikolay Nikolov,
Tiffany Kataria,
Drake Deming,
Gilda E. Ballester
Abstract:
WASP-121b is a transiting gas giant exoplanet orbiting close to its Roche limit, with an inflated radius nearly double that of Jupiter and a dayside temperature comparable to a late M dwarf photosphere. Secondary eclipse observations covering the 1.1-1.6 micron wavelength range have revealed an atmospheric thermal inversion on the dayside hemisphere, likely caused by high altitude absorption at op…
▽ More
WASP-121b is a transiting gas giant exoplanet orbiting close to its Roche limit, with an inflated radius nearly double that of Jupiter and a dayside temperature comparable to a late M dwarf photosphere. Secondary eclipse observations covering the 1.1-1.6 micron wavelength range have revealed an atmospheric thermal inversion on the dayside hemisphere, likely caused by high altitude absorption at optical wavelengths. Here we present secondary eclipse observations made with the Hubble Space Telescope Wide Field Camera 3 spectrograph that extend the wavelength coverage from 1.1 micron down to 0.8 micron. To determine the atmospheric properties from the measured eclipse spectrum, we performed a retrieval analysis assuming chemical equilibrium, with the effects of thermal dissociation and ionization included. Our best-fit model provides a good fit to the data with reduced chi^2=1.04. The data diverge from a blackbody spectrum and instead exhibit emission due to H- shortward of 1.1 micron. The best-fit model does not reproduce a previously reported bump in the spectrum at 1.25 micron, possibly indicating this feature is a statistical fluctuation in the data rather than a VO emission band as had been tentatively suggested. We estimate an atmospheric metallicity of [M/H]=1.09(-0.69,+0.57), and fit for the carbon and oxygen abundances separately, obtaining [C/H]=-0.29(-0.48,+0.61) and [O/H]=0.18(-0.60,+0.64). The corresponding carbon-to-oxygen ratio is C/O=0.49(-0.37,+0.65), which encompasses the solar value of 0.54, but has a large uncertainty.
△ Less
Submitted 16 July, 2019; v1 submitted 14 June, 2019;
originally announced June 2019.
-
The carbon-to-oxygen ratio: implications for the spectra of hydrogen-dominated exoplanet atmospheres
Authors:
Benjamin Drummond,
Aarynn L. Carter,
Eric Hébrard,
Nathan J. Mayne,
David K. Sing,
Thomas M. Evans,
Jayesh Goyal
Abstract:
We present results from one-dimensional atmospheric simulations investigating the effect of varying the carbon-to-oxygen (C/O) ratio on the thermal structure, chemical composition and transmission and emission spectra, for irradiated hydrogen-dominated atmospheres. We find that each of these properties of the atmosphere are strongly dependent on the individual abundances (relative to hydrogen) of…
▽ More
We present results from one-dimensional atmospheric simulations investigating the effect of varying the carbon-to-oxygen (C/O) ratio on the thermal structure, chemical composition and transmission and emission spectra, for irradiated hydrogen-dominated atmospheres. We find that each of these properties of the atmosphere are strongly dependent on the individual abundances (relative to hydrogen) of carbon and oxygen. We confirm previous findings that different chemical equilibrium compositions result from different sets of element abundances but with the same C/O ratio. We investigate the effect of this difference in composition on the thermal structure and simulated spectra. We also simulate observations using the PandExo tool and show that these differences are observationally significant with current (i.e. Hubble Space Telescope) and future (i.e. James Webb Space Telescope) instruments. We conclude that it is important to consider the full set of individual element abundances, with respect to hydrogen, rather than the ratios of only two elements, such as the C/O ratio, particularly when comparing model predictions with observed transmission and emission spectra.
△ Less
Submitted 26 March, 2019;
originally announced March 2019.
-
A reduced chemical scheme for modelling warm to hot hydrogen-dominated atmospheres
Authors:
Olivia Venot,
Roda Bounaceur,
Michel Dobrijevic,
Eric Hébrard,
Thibault Cavalié,
Pascal Tremblin,
Benjamin Drummond,
Benjamin Charnay
Abstract:
Three dimensional models that account for chemistry are useful tools to predict the chemical composition of (exo)planet and brown dwarf atmospheres and interpret observations of future telescopes. Recent Juno observations of the NH3 tropospheric distribution in Jupiter also indicate that 3D chemical modelling may be necessary to constrain the deep composition of the giant planets of the Solar Syst…
▽ More
Three dimensional models that account for chemistry are useful tools to predict the chemical composition of (exo)planet and brown dwarf atmospheres and interpret observations of future telescopes. Recent Juno observations of the NH3 tropospheric distribution in Jupiter also indicate that 3D chemical modelling may be necessary to constrain the deep composition of the giant planets of the Solar System. However, due to the high computational cost of chemistry calculations, 3D chemical modelling has so far been limited. A solution to include chemical kinetics in a 3D model is to use a reduced chemical scheme. In this view, we have developed a reduced scheme containing 30 species and 181 reversible reactions, from the full chemical scheme of Venot et al. 2012. This scheme reproduces accurately the vertical profiles of the observable species (H2O, CH4, CO, CO2, NH3, and HCN) and has a large range of validity. It can be used to study all kind of warm atmospheres, except hot C-rich ones, which contains a high amount of C2H2. It can be used in 1D models, for fast computations, but also in 3D models for hot giant (exo)planet and brown dwarf atmospheres.
△ Less
Submitted 13 February, 2019;
originally announced February 2019.
-
Thermo-compositional diabatic convection in the atmospheres of brown dwarfs and in Earth's atmosphere and oceans
Authors:
P. Tremblin,
T. Padioleau,
M. Phillips,
G. Chabrier,
I. Baraffe,
S. Fromang,
E. Audit,
H. Bloch,
A. J. Burgasser,
B. Drummond,
M. Gonzalez,
P. Kestener,
S. Kokh,
P. -O. Lagage,
M. Stauffert
Abstract:
By generalizing the theory of convection to any type of thermal and compositional source terms (diabatic processes), we show that thermohaline convection in Earth oceans, fingering convection in stellar atmospheres, and moist convection in Earth atmosphere are deriving from the same general diabatic convective instability. We show also that "radiative convection" triggered by CO/CH4 transition wit…
▽ More
By generalizing the theory of convection to any type of thermal and compositional source terms (diabatic processes), we show that thermohaline convection in Earth oceans, fingering convection in stellar atmospheres, and moist convection in Earth atmosphere are deriving from the same general diabatic convective instability. We show also that "radiative convection" triggered by CO/CH4 transition with radiative transfer in the atmospheres of brown dwarfs is analog to moist and thermohaline convection. We derive a generalization of the mixing length theory to include the effect of source terms in 1D codes. We show that CO/CH4 radiative convection could significantly reduce the temperature gradient in the atmospheres of brown dwarfs similarly to moist convection in Earth atmosphere thus possibly explaining the reddening in brown-dwarf spectra. By using idealized two-dimensional hydrodynamic simulations in the Ledoux unstable regime, we show that compositional source terms can indeed provoke a reduction of the temperature gradient. The L/T transition could be explained by a bifurcation between the adiabatic and diabatic convective transports and could be seen as a giant cooling crisis: an analog of the boiling crisis in liquid/steam-water convective flows. This mechanism with other chemical transitions could be present in many giant and earth-like exoplanets. The study of the impact of different parameters (effective temperature, compositional changes) on CO/CH4 radiative convection and the analogy with Earth moist and thermohaline convection is opening the possibility to use brown dwarfs to better understand some aspects of the physics at play in the climate of our own planet.
△ Less
Submitted 10 February, 2019;
originally announced February 2019.
-
The Limits of the Primitive Equations of Dynamics for Warm, Slowly Rotating Small Neptunes and Super Earths
Authors:
N. J. Mayne,
B. Drummond,
F. Debras,
E. Jaupart,
J. Manners,
I. A. Boutle,
I. Baraffe,
K. Kohary
Abstract:
We present significant differences in the simulated atmospheric flow for warm, tidally-locked small Neptunes and super Earths (based on a nominal GJ 1214b) when solving the simplified, and commonly used, primitive dynamical equations or the full Navier-Stokes equations. The dominant prograde, superrotating zonal jet is markedly different between the simulations which are performed using practicall…
▽ More
We present significant differences in the simulated atmospheric flow for warm, tidally-locked small Neptunes and super Earths (based on a nominal GJ 1214b) when solving the simplified, and commonly used, primitive dynamical equations or the full Navier-Stokes equations. The dominant prograde, superrotating zonal jet is markedly different between the simulations which are performed using practically identical numerical setups, within the same model. The differences arise due to the breakdown of the so-called `shallow-fluid' and traditional approximations, which worsens when rotation rates are slowed, and day-night temperature contrasts are increased. The changes in the zonal advection between simulations solving the full and simplified equations, give rise to significant differences in the atmospheric redistribution of heat, altering the position of the hottest part of the atmosphere and temperature contrast between the day and night sides. The implications for the atmospheric chemistry and, therefore, observations need to be studied with a model including a more detailed treatment of the radiative transfer and chemistry. Small Neptunes and super Earths are extremely abundant and important, potentially bridging the structural properties (mass, radius, composition) of terrestrial and gas giant planets. Our results indicate care is required when interpreting the output of models solving the primitive equations of motion for such planets.
△ Less
Submitted 6 December, 2018;
originally announced December 2018.
-
Fully scalable forward model grid of exoplanet transmission spectra
Authors:
Jayesh M. Goyal,
Hannah R. Wakeford,
Nathan J. Mayne,
Nikole K. Lewis,
Benjamin Drummond,
David K. Sing
Abstract:
Simulated exoplanet transmission spectra are critical for planning and interpretation of observations and to explore the sensitivity of spectral features to atmospheric thermochemical processes. We present a publicly available generic model grid of planetary transmission spectra, scalable to a wide range of H$_2$/He dominated atmospheres. The grid is computed using the 1D/2D atmosphere model ATMO…
▽ More
Simulated exoplanet transmission spectra are critical for planning and interpretation of observations and to explore the sensitivity of spectral features to atmospheric thermochemical processes. We present a publicly available generic model grid of planetary transmission spectra, scalable to a wide range of H$_2$/He dominated atmospheres. The grid is computed using the 1D/2D atmosphere model ATMO for two different chemical scenarios, first considering local condensation only, secondly considering global condensation and removal of species from the atmospheric column (rainout). The entire grid consists of 56,320 model simulations across 22 equilibrium temperatures (400 - 2600 K), four planetary gravities (5 - 50 ms$^{-2}$), five atmospheric metallicities (1x - 200x), four C/O ratios (0.35 - 1.0), four scattering haze parameters, four uniform cloud parameters, and two chemical scenarios. We derive scaling equations which can be used with this grid, for a wide range of planet-star combinations. We validate this grid by comparing it with other model transmission spectra available in the literature. We highlight some of the important findings, such as the rise of SO$_2$ features at 100x solar metallicity, differences in spectral features at high C/O ratios between two condensation approaches, the importance of VO features without TiO to constrain the limb temperature and features of TiO/VO both, to constrain the condensation processes. Finally, this generic grid can be used to plan future observations using the HST, VLT, JWST and various other telescopes. The fine variation of parameters in the grid also allows it to be incorporated in a retrieval framework, with various machine learning techniques.
△ Less
Submitted 30 October, 2018;
originally announced October 2018.
-
An optical transmission spectrum for the ultra-hot Jupiter WASP-121b measured with the Hubble Space Telescope
Authors:
Thomas M. Evans,
David K. Sing,
Jayesh Goyal,
Nikolay Nikolov,
Mark S. Marley,
Kevin Zahnle,
Gregory W. Henry,
Joanna K. Barstow,
Munazza K. Alam,
Jorge Sanz-Forcada,
Tiffany Kataria,
Nikole K. Lewis,
Panayotis Lavvas,
Gilda E. Ballester,
Lotfi Ben-Jaffel,
Sarah D. Blumenthal,
Vincent Bourrier,
Benjamin Drummond,
Antonio Garcia Munoz,
Mercedes Lopez-Morales,
Pascal Tremblin,
David Ehrenreich,
Hannah R. Wakeford,
Lars A. Buchhave,
Alain Lecavelier des Etangs
, et al. (2 additional authors not shown)
Abstract:
We present an atmospheric transmission spectrum for the ultra-hot Jupiter WASP-121b, measured using the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST). Across the 0.47-1 micron wavelength range, the data imply an atmospheric opacity comparable to - and in some spectroscopic channels exceeding - that previously measured at near-infrared wavelengths (1.15-1.65 m…
▽ More
We present an atmospheric transmission spectrum for the ultra-hot Jupiter WASP-121b, measured using the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST). Across the 0.47-1 micron wavelength range, the data imply an atmospheric opacity comparable to - and in some spectroscopic channels exceeding - that previously measured at near-infrared wavelengths (1.15-1.65 micron). Wavelength-dependent variations in the opacity rule out a gray cloud deck at a confidence level of 3.8-sigma and may instead be explained by VO spectral bands. We find a cloud-free model assuming chemical equilibrium for a temperature of 1500K and metal enrichment of 10-30x solar matches these data well. Using a free-chemistry retrieval analysis, we estimate a VO abundance of -6.6(-0.3,+0.2) dex. We find no evidence for TiO and place a 3-sigma upper limit of -7.9 dex on its abundance, suggesting TiO may have condensed from the gas phase at the day-night limb. The opacity rises steeply at the shortest wavelengths, increasing by approximately five pressure scale heights from 0.47 to 0.3 micron in wavelength. If this feature is caused by Rayleigh scattering due to uniformly-distributed aerosols, it would imply an unphysically high temperature of 6810+/-1530K. One alternative explanation for the short-wavelength rise is absorption due to SH (mercapto radical), which has been predicted as an important product of non-equilibrium chemistry in hot Jupiter atmospheres. Irrespective of the identity of the NUV absorber, it likely captures a significant amount of incident stellar radiation at low pressures, thus playing a significant role in the overall energy budget, thermal structure, and circulation of the atmosphere.
△ Less
Submitted 25 October, 2018;
originally announced October 2018.
-
The 3D thermal, dynamical and chemical structure of the atmosphere of HD 189733b: implications of wind-driven chemistry for the emission phase curve
Authors:
Benjamin Drummond,
Nathan J. Mayne,
James Manners,
Isabelle Baraffe,
Jayesh Goyal,
Pascal Tremblin,
David K. Sing,
Krisztian Kohary
Abstract:
In this paper we present three-dimensional atmospheric simulations of the hot Jupiter HD~189733b under two different scenarios: local chemical equilibrium and including advection of the chemistry by the resolved wind. Our model consistently couples the treatment of dynamics, radiative transfer and chemistry, completing the feedback cycle between these three important processes. The effect of wind-…
▽ More
In this paper we present three-dimensional atmospheric simulations of the hot Jupiter HD~189733b under two different scenarios: local chemical equilibrium and including advection of the chemistry by the resolved wind. Our model consistently couples the treatment of dynamics, radiative transfer and chemistry, completing the feedback cycle between these three important processes. The effect of wind--driven advection on the chemical composition is qualitatively similar to our previous results for the warmer atmosphere of HD~209458b, found using the same model. However, we find more significant alterations to both the thermal and dynamical structure for the cooler atmosphere of HD~189733b, with changes in both the temperature and wind velocities reaching $\sim10\%$. We also present the contribution function, diagnosed from our simulations, and show that wind--driven chemistry has a significant impact on its three--dimensional structure, particularly for regions where methane is an important absorber. Finally, we present emission phase curves from our simulations and show the significant effect of wind--driven chemistry on the thermal emission, particularly within the 3.6 \textmu m Spitzer/IRAC channel.
△ Less
Submitted 23 October, 2018;
originally announced October 2018.
-
Exonephology: Transmission spectra from a 3D simulated cloudy atmosphere of HD209458b
Authors:
S. Lines,
J. Manners,
N. J. Mayne,
J. Goyal,
A. L. Carter,
I. A. Boutle,
E. K. H. Lee,
Ch. Helling,
B. Drummond,
D. M. Acreman,
D. K. Sing
Abstract:
We present high resolution transmission spectra, calculated directly from a 3D radiative-hydrodynamics simulation that includes kinetic cloud formation, for HD209458b. We find that the high opacity of our vertically extensive cloud deck, composed of a large number density of sub-micron particles, flattens the transmission spectrum and obscures spectral features identified in observed data. We use…
▽ More
We present high resolution transmission spectra, calculated directly from a 3D radiative-hydrodynamics simulation that includes kinetic cloud formation, for HD209458b. We find that the high opacity of our vertically extensive cloud deck, composed of a large number density of sub-micron particles, flattens the transmission spectrum and obscures spectral features identified in observed data. We use the PandExo simulator to explore features of our HD209458b spectrum which may be detectable with the James Webb Space Telescope (JWST). We determine that an 8 - 12 micron absorption feature attributed to the mixed-composition, predominantly silicate cloud particles is a viable marker for the presence of cloud. Further calculations explore, and trends are identified with, variations in cloud opacity, composition heterogeneity and artificially scaled gravitational settling on the transmission spectrum. Principally, by varying the upper extent of our cloud decks, rainout is identified to be a key process for the dynamical atmospheres of hot-Jupiters and shown to dramatically alter the resulting spectrum. Our synthetic transmission spectra, obtained from the most complete, forward atmosphere simulations to--date, allow us to explore the model's ability to conform with observations. Such comparisons can provide insight into the physical processes either missing, or requiring improvement.
△ Less
Submitted 17 August, 2018;
originally announced August 2018.
-
An absolute sodium abundance for a cloud-free 'hot Saturn' exoplanet
Authors:
Nikolay Nikolov,
David K. Sing,
Jonathan J. Fortney,
Jayesh M. Goyal,
Benjamin Drummond,
Tom M. Evans,
Neale P. Gibson,
Ernst J. W. De Mooij,
Zafar Rustamkulov,
Hannah R. Wakeford,
Barry Smalley,
Adam J. Burgasser,
Coel Hellier,
Christiane Helling,
Nathan J. Mayne,
Nikku Madhusudhan,
Tiffany Kataria,
Josef Baines,
Aarynn L. Carter,
Gilda E. Ballester,
Joanna K. Barstow,
Jack McCleery,
Jessica J. Spake
Abstract:
Broad absorption signatures from alkali metals, such as the sodium (Na I) and potassium (K I) resonance doublets, have long been predicted in the optical atmospheric spectra of cloud-free irradiated gas-giant exoplanets1,2,3. However, observations have only revealed the narrow cores of these features rather than the full pressure-broadened profiles4-6. Cloud and haze opacity at the day-night plane…
▽ More
Broad absorption signatures from alkali metals, such as the sodium (Na I) and potassium (K I) resonance doublets, have long been predicted in the optical atmospheric spectra of cloud-free irradiated gas-giant exoplanets1,2,3. However, observations have only revealed the narrow cores of these features rather than the full pressure-broadened profiles4-6. Cloud and haze opacity at the day-night planetary terminator are considered responsible for obscuring the absorption-line wings, which hinders constraints on absolute atmospheric abundances7-9. Here we present an optical transmission spectrum for the 'hot-Saturn' WASP-96b obtained with the Very Large Telescope, which exhibits the complete pressure-broadened profile of the sodium absorption feature. The spectrum is in excellent agreement with cloud-free, solar-abundance models assuming chemical equilibrium. We are able to measure a precise, absolute sodium abundance of logε_Na=6.9+0.6-0.4, and use it as a proxy to the planet's atmospheric metallicity relative to the solar value (Z_p/Z_\star=2.3+8.9/--1.7). This result is consistent with the mass-metallicity trend observed for solar-system planets and exoplanets10-12.
△ Less
Submitted 15 June, 2018;
originally announced June 2018.
-
The Transiting Exoplanet Community Early Release Science Program for JWST
Authors:
Jacob L. Bean,
Kevin B. Stevenson,
Natalie M. Batalha,
Zachory Berta-Thompson,
Laura Kreidberg,
Nicolas Crouzet,
Björn Benneke,
Michael R. Line,
David K. Sing,
Hannah R. Wakeford,
Heather A. Knutson,
Eliza M. -R. Kempton,
Jean-Michel Désert,
Ian Crossfield,
Natasha E. Batalha,
Julien de Wit,
Vivien Parmentier,
Joseph Harrington,
Julianne I. Moses,
Mercedes Lopez-Morales,
Munazza K. Alam,
Jasmina Blecic,
Giovanni Bruno,
Aarynn L. Carter,
John W. Chapman
, et al. (77 additional authors not shown)
Abstract:
The James Webb Space Telescope (JWST) presents the opportunity to transform our understanding of planets and the origins of life by revealing the atmospheric compositions, structures, and dynamics of transiting exoplanets in unprecedented detail. However, the high-precision, time-series observations required for such investigations have unique technical challenges, and prior experience with other…
▽ More
The James Webb Space Telescope (JWST) presents the opportunity to transform our understanding of planets and the origins of life by revealing the atmospheric compositions, structures, and dynamics of transiting exoplanets in unprecedented detail. However, the high-precision, time-series observations required for such investigations have unique technical challenges, and prior experience with other facilities indicates that there will be a steep learning curve when JWST becomes operational. In this paper we describe the science objectives and detailed plans of the Transiting Exoplanet Community Early Release Science (ERS) Program, which is a recently approved program for JWST observations early in Cycle 1. The goal of this project, for which the obtained data will have no exclusive access period, is to accelerate the acquisition and diffusion of technical expertise for transiting exoplanet observations with JWST, while also providing a compelling set of representative datasets that will enable immediate scientific breakthroughs. The Transiting Exoplanet Community ERS Program will exercise the time-series modes of all four JWST instruments that have been identified as the consensus highest priorities, observe the full suite of transiting planet characterization geometries (transits, eclipses, and phase curves), and target planets with host stars that span an illustrative range of brightnesses. The observations in this program were defined through an inclusive and transparent process that had participation from JWST instrument experts and international leaders in transiting exoplanet studies. Community engagement in the project will be centered on a two-phase Data Challenge that culminates with the delivery of planetary spectra, time-series instrument performance reports, and open-source data analysis toolkits in time to inform the agenda for Cycle 2 of the JWST mission.
△ Less
Submitted 3 September, 2018; v1 submitted 13 March, 2018;
originally announced March 2018.
-
Simulating the cloudy atmospheres of HD 209458 b and HD 189733 b with the 3D Met Office Unified Model
Authors:
S. Lines,
N. J. Mayne,
Ian A. Boutle,
James Manners,
Graham K. H. Lee,
Ch. Helling,
Benjamin Drummond,
David S. Amundsen,
Jayesh Goyal,
David M. Acreman,
Pascal Tremblin,
Max Kerslake
Abstract:
To understand and compare the 3D atmospheric structure of HD 209458 b and HD 189733 b, focusing on the formation and distribution of cloud particles, as well as their feedback on the dynamics and thermal profile. We couple the 3D Met Office Unified Model (UM), including detailed treatments of atmospheric radiative transfer and dynamics, to a kinetic cloud formation scheme. The resulting model self…
▽ More
To understand and compare the 3D atmospheric structure of HD 209458 b and HD 189733 b, focusing on the formation and distribution of cloud particles, as well as their feedback on the dynamics and thermal profile. We couple the 3D Met Office Unified Model (UM), including detailed treatments of atmospheric radiative transfer and dynamics, to a kinetic cloud formation scheme. The resulting model self--consistently solves for the formation of condensation seeds, surface growth and evaporation, gravitational settling and advection, cloud radiative feedback via absorption and, crucially, scattering. Fluxes directly obtained from the UM are used to produce synthetic SEDs and phase curves. Our simulations show extensive cloud formation in both planets, however, cooler temperatures in the HD 189733 b result in higher cloud particle number densities. Sub-micron particles are suspended by vertical flows leading to extensive upper-atmosphere cloud cover. A combination of meridional advection and efficient cloud formation in cooler high latitude regions, result in enhanced cloud coverage for latitudes > 30 degrees and leads to a zonally banded structure for all our simulations. The cloud bands extend around the entire planet(s), as the temperatures, even on the day side, remain below the condensation temperature of silicates and oxides. Therefore, our simulated optical phase curve for HD 209458 b shows no `offset', in contrast to observations. Efficient scattering by cloud results an atmospheric cooling of up to 250K, and an advection-driven fluctuating cloud opacity causes temporal variability in the thermal emission. The inclusion of this fundamental cloud-atmosphere radiative feedback leads to significant differences with approaches neglecting these physical elements and suggests both a note of caution of interpretations neglecting such cloud feedback and scattering, and merits further study.
△ Less
Submitted 1 March, 2018;
originally announced March 2018.
-
Observable signatures of wind--driven chemistry with a fully consistent three dimensional radiative hydrodynamics model of HD 209458b
Authors:
Benjamin Drummond,
N. J. Mayne,
James Manners,
Aarynn L. Carter,
Ian A. Boutle,
Isabelle Baraffe,
Eric Hebrard,
Pascal Tremblin,
David K. Sing,
David S. Amundsen,
Dave Acreman
Abstract:
We present a study of the effect of wind-driven advection on the chemical composition of hot Jupiter atmospheres using a fully-consistent 3D hydrodynamics, chemistry and radiative transfer code, the Met Office Unified Model (UM). Chemical modelling of exoplanet atmospheres has primarily been restricted to 1D models that cannot account for 3D dynamical processes. In this work we couple a chemical r…
▽ More
We present a study of the effect of wind-driven advection on the chemical composition of hot Jupiter atmospheres using a fully-consistent 3D hydrodynamics, chemistry and radiative transfer code, the Met Office Unified Model (UM). Chemical modelling of exoplanet atmospheres has primarily been restricted to 1D models that cannot account for 3D dynamical processes. In this work we couple a chemical relaxation scheme to the UM to account for the chemical interconversion of methane and carbon monoxide. This is done consistently with the radiative transfer meaning that departures from chemical equilibrium are included in the heating rates (and emission) and hence complete the feedback between the dynamics, thermal structure and chemical composition. In this letter we simulate the well studied atmosphere of HD~209458b. We find that the combined effect of horizontal and vertical advection leads to an increase in the methane abundance by several orders of magnitude; directly opposite to the trend found in previous works. Our results demonstrate the need to include 3D effects when considering the chemistry of hot Jupiter atmospheres. We calculate transmission and emission spectra, as well as the emission phase curve, from our simulations. We conclude that gas-phase non-equilibrium chemistry is unlikely to explain the model-observation discrepancy in the 4.5\,{\textmu m} {\it Spitzer}/IRAC channel. However, we highlight other spectral regions, observable with the James Webb Space Telescope, where signatures of wind-driven chemistry are more prominant.
△ Less
Submitted 26 February, 2018;
originally announced February 2018.
-
The effect of metallicity on the atmospheres of exoplanets with fully coupled 3D hydrodynamics, equilibrium chemistry, and radiative transfer
Authors:
Benjamin Drummond,
N. J. Mayne,
Isabelle Baraffe,
Pascal Tremblin,
James Manners,
David S. Amundsen,
Jayesh Goyal,
Dave Acreman
Abstract:
In this work we have performed a series of simulations of the atmosphere of GJ~1214b assuming different metallicities using the Met Office Unified Model (UM). The UM is a general circulation model (GCM) that solves the deep, non-hydrostatic equations of motion and uses a flexible and accurate radiative transfer scheme, based on the two-stream and correlated-$k$ approximations, to calculate the hea…
▽ More
In this work we have performed a series of simulations of the atmosphere of GJ~1214b assuming different metallicities using the Met Office Unified Model (UM). The UM is a general circulation model (GCM) that solves the deep, non-hydrostatic equations of motion and uses a flexible and accurate radiative transfer scheme, based on the two-stream and correlated-$k$ approximations, to calculate the heating rates. In this work we consistently couple a well-tested Gibbs energy minimisation scheme to solve for the chemical equilibrium abundances locally in each grid cell for a general set of elemental abundances, further improving the flexibility and accuracy of the model. As the metallicity of the atmosphere is increased we find significant changes in the dynamical and thermal structure, with subsequent implications for the simulated phase curve. The trends that we find are qualitatively consistent with previous works, though with quantitative differences. We investigate in detail the effect of increasing the metallicity by splitting the mechanism into constituents, involving the mean molecular weight, the heat capacity and the opacities. We find the opacity effect to be the dominant mechanism in altering the circulation and thermal structure. This result highlights the importance of accurately computing the opacities and radiative transfer in 3D GCMs.
△ Less
Submitted 3 January, 2018;
originally announced January 2018.
-
The Complete transmission spectrum of WASP-39b with a precise water constraint
Authors:
Hannah R. Wakeford,
David K. Sing,
Drake Deming,
Nikole K. Lewis,
Jayesh Goyal,
Tom J. Wilson,
Joanna Barstow,
Tiffany Kataria,
Benjamin Drummond,
Thomas M. Evans,
Aarynn L. Carter,
Nikolay Nikolov,
Heather A. Knutson,
Gilda E. Ballester,
Avi M. Mandell
Abstract:
WASP-39b is a hot Saturn-mass exoplanet with a predicted clear atmosphere based on observations in the optical and infrared. Here we complete the transmission spectrum of the atmosphere with observations in the near-infrared (NIR) over three water absorption features with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) G102 (0.8-1.1 microns) and G141 (1.1-1.7 microns) spectroscopic gri…
▽ More
WASP-39b is a hot Saturn-mass exoplanet with a predicted clear atmosphere based on observations in the optical and infrared. Here we complete the transmission spectrum of the atmosphere with observations in the near-infrared (NIR) over three water absorption features with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) G102 (0.8-1.1 microns) and G141 (1.1-1.7 microns) spectroscopic grisms. We measure the predicted high amplitude H2O feature centered at 1.4 microns, and the smaller amplitude features at 0.95 and 1.2 microns, with a maximum water absorption amplitude of 2.4 planetary scale heights. We incorporate these new NIR measurements into previously published observational measurements to complete the transmission spectrum from 0.3-5 microns. From these observed water features, combined with features in the optical and IR, we retrieve a well constrained temperature Teq = 1030(+30,-20) K, and atmospheric metallicity 151 (+48,-46)x solar which is relatively high with respect to the currently established mass-metallicity trends. This new measurement in the Saturn-mass range hints at further diversity in the planet formation process relative to our solar system giants.
△ Less
Submitted 28 November, 2017;
originally announced November 2017.
-
A better characterization of the chemical composition of exoplanets atmospheres with ARIEL
Authors:
Olivia Venot,
Benjamin Drummond,
Yamila Miguel,
Ingo P. Waldmann,
Enzo Pascale,
Tiziano Zingales
Abstract:
Since the discovery of the first extrasolar planet more than twenty years ago, we have discovered nearly four thousand planets orbiting stars other than the Sun. Current observational instruments (on board the Hubble Space Telescope, Spitzer, and on ground-based facilities) have allowed the scientific community to obtain important information on the physical and chemical properties of these planet…
▽ More
Since the discovery of the first extrasolar planet more than twenty years ago, we have discovered nearly four thousand planets orbiting stars other than the Sun. Current observational instruments (on board the Hubble Space Telescope, Spitzer, and on ground-based facilities) have allowed the scientific community to obtain important information on the physical and chemical properties of these planets. However, for a more in-depth characterisation of these worlds, more powerful telescopes are needed. Thanks to the high sensitivity of their instruments, the next generation of space observatories (e.g. JWST, ARIEL) will provide observations of unprecedented quality, allowing us to extract far more information than what was previously possible. Such high quality observations will provide constraints on theoretical models of exoplanet atmospheres and lead to a greater understanding of their physics and chemistry. Important modelling efforts have been carried out during the past few years, showing that numerous parameters and processes (such as the elemental abundances, temperature, mixing, etc.) are likely to affect the atmospheric composition of exoplanets and subsequently the observable spectra. In this manuscript, we review the different parameters that can influence the molecular composition of exoplanet atmospheres. We show that the high-precision of ARIEL observations will improve our view and characterisation of exoplanet atmospheres. We also consider future developments that are necessary to improve atmospheric models, driven by the need to interpret the available observations.
△ Less
Submitted 2 May, 2018; v1 submitted 22 November, 2017;
originally announced November 2017.
-
A library of ATMO forward model transmission spectra for hot Jupiter exoplanets
Authors:
Jayesh M. Goyal,
Nathan Mayne,
David K. Sing,
Benjamin Drummond,
Pascal Tremblin,
David S. Amundsen,
Thomas Evans,
Aarynn L. Carter,
Jessica Spake,
Isabelle Baraffe,
Nikolay Nikolov,
James Manners,
Gilles Chabrier,
Eric Hebrard
Abstract:
We present a grid of forward model transmission spectra, adopting an isothermal temperature-pressure profile, alongside corresponding equilibrium chemical abundances for 117 observationally significant hot exoplanets (Equilibrium Temperatures of 547-2710 K). This model grid has been developed using a 1D radiative-convective-chemical equilibrium model termed ATMO, with up-to-date high temperature o…
▽ More
We present a grid of forward model transmission spectra, adopting an isothermal temperature-pressure profile, alongside corresponding equilibrium chemical abundances for 117 observationally significant hot exoplanets (Equilibrium Temperatures of 547-2710 K). This model grid has been developed using a 1D radiative-convective-chemical equilibrium model termed ATMO, with up-to-date high temperature opacities. We present an interpretation of observations of ten exoplanets, including best fit parameters and $χ^{2}$ maps. In agreement with previous works, we find a continuum from clear to hazy/cloudy atmospheres for this sample of hot Jupiters. The data for all the 10 planets are consistent with sub-solar to solar C/O ratio, 0.005 to 10 times solar metallicity and water rather than a methane dominated infrared spectra. We then explore the range of simulated atmospheric spectra for different exoplanets, based on characteristics such as temperature, metallicity, C/O-ratio, haziness and cloudiness. We find a transition value for the metallicity between 10 and 50 times solar, which leads to substantial changes in the transmission spectra. We also find a transition value of C/O ratio, from water to carbon species dominated infrared spectra, as found by previous works, revealing a temperature dependence of this transition point ranging from $\sim$0.56 to $\sim$1-1.3 for equilibrium temperatures from $\sim$900 to $\sim$2600 K. We highlight the potential of the spectral features of HCN and C$_2$H$_2$ to constrain the metallicities and C/O ratios of planets, using JWST observations. Finally, our entire grid ($\sim$460,000 simulations) is publicly available and can be used directly with the JWST simulator PandExo for planning observations.
△ Less
Submitted 29 July, 2020; v1 submitted 27 October, 2017;
originally announced October 2017.
-
Cloudless atmospheres for young low-gravity substellar objects
Authors:
P. Tremblin,
G. Chabrier,
I. Baraffe,
Michael. C. Liu,
E. A. Magnier,
P. -O. Lagage,
C. Alves de Oliveira,
A. J. Burgasser,
D. S. Amundsen,
B. Drummond
Abstract:
Atmospheric modeling of low-gravity (VL-G) young brown dwarfs remains a challenge. The presence of very thick clouds has been suggested because of their extremely red near-infrared (NIR) spectra, but no cloud models provide a good fit to the data with a radius compatible with evolutionary models for these objects. We show that cloudless atmospheres assuming a temperature gradient reduction caused…
▽ More
Atmospheric modeling of low-gravity (VL-G) young brown dwarfs remains a challenge. The presence of very thick clouds has been suggested because of their extremely red near-infrared (NIR) spectra, but no cloud models provide a good fit to the data with a radius compatible with evolutionary models for these objects. We show that cloudless atmospheres assuming a temperature gradient reduction caused by fingering convection provides a very good model to match the observed VL-G NIR spectra. The sequence of extremely red colors in the NIR for atmospheres with effective temperature from ~2000 K down to ~1200 K is very well reproduced with predicted radii typical of young low-gravity objects. Future observations with NIRSPEC and MIRI on the James Webb Space Telescope (JWST) will provide more constrains in the mid-infrared, helping to confirm/refute whether or not the NIR reddening is caused by fingering convection. We suggest that the presence/absence of clouds will be directly determined by the silicate absorption features that can be observed with MIRI. JWST will therefore be able to better characterize the atmosphere of these hot young brown dwarfs and their low-gravity exoplanet analogues.
△ Less
Submitted 7 October, 2017;
originally announced October 2017.
-
Advection of potential temperature in the atmosphere of irradiated exoplanets: a robust mechanism to explain radius inflation
Authors:
P. Tremblin,
G. Chabrier,
N. J. Mayne,
D. S. Amundsen,
I. Baraffe,
F. Debras,
B. Drummond,
J. Manners,
S. Fromang
Abstract:
The anomalously large radii of strongly irradiated exoplanets have remained a major puzzle in astronomy. Based on a 2D steady state atmospheric circulation model, the validity of which is assessed by comparison to 3D calculations, we reveal a new mechanism, namely the advection of the potential temperature due to mass and longitudinal momentum conservation, a process occuring in the Earth's atmosp…
▽ More
The anomalously large radii of strongly irradiated exoplanets have remained a major puzzle in astronomy. Based on a 2D steady state atmospheric circulation model, the validity of which is assessed by comparison to 3D calculations, we reveal a new mechanism, namely the advection of the potential temperature due to mass and longitudinal momentum conservation, a process occuring in the Earth's atmosphere or oceans. At depth, the vanishing heating flux forces the atmospheric structure to converge to a hotter adiabat than the one obtained with 1D calculations, implying a larger radius for the planet. Not only do the calculations reproduce the observed radius of HD209458b, but also the observed correlation between radius inflation and irradiation for transiting planets. Vertical advection of potential temperature induced by non uniform atmospheric heating thus provides a robust mechanism explaining the inflated radii of irradiated hot Jupiters.
△ Less
Submitted 18 April, 2017;
originally announced April 2017.
-
Exploring the climate of Proxima B with the Met Office Unified Model
Authors:
Ian A. Boutle,
Nathan J. Mayne,
Benjamin Drummond,
James Manners,
Jayesh Goyal,
F. Hugo Lambert,
David M. Acreman,
Paul D. Earnshaw
Abstract:
We present results of simulations of the climate of the newly discovered planet Proxima Centauri B, performed using the Met Office Unified Model (UM). We examine the responses of both an `Earth-like' atmosphere and simplified nitrogen and trace carbon dioxide atmosphere to the radiation likely received by Proxima Centauri B. Additionally, we explore the effects of orbital eccentricity on the plane…
▽ More
We present results of simulations of the climate of the newly discovered planet Proxima Centauri B, performed using the Met Office Unified Model (UM). We examine the responses of both an `Earth-like' atmosphere and simplified nitrogen and trace carbon dioxide atmosphere to the radiation likely received by Proxima Centauri B. Additionally, we explore the effects of orbital eccentricity on the planetary conditions using a range of eccentricities guided by the observational constraints. Overall, our results are in agreement with previous studies in suggesting Proxima Centauri B may well have surface temperatures conducive to the presence of liquid water. Moreover, we have expanded the parameter regime over which the planet may support liquid water to higher values of eccentricity (>= 0.1) and lower incident fluxes (881.7 Wm-2) than previous work. This increased parameter space arises because of the low sensitivity of the planet to changes in stellar flux, a consequence of the stellar spectrum and orbital configuration. However, we also find interesting differences from previous simulations, such as cooler mean surface temperatures for the tidally-locked case. Finally, we have produced high resolution planetary emission and reflectance spectra, and highlight signatures of gases vital to the evolution of complex life on Earth (oxygen, ozone and carbon dioxide).
△ Less
Submitted 24 August, 2018; v1 submitted 27 February, 2017;
originally announced February 2017.
-
The UK Met Office GCM with a sophisticated radiation scheme applied to the hot Jupiter HD 209458b
Authors:
David S. Amundsen,
Nathan J. Mayne,
Isabelle Baraffe,
James Manners,
Pascal Tremblin,
Benjamin Drummond,
Chris Smith,
David M. Acreman,
Derek Homeier
Abstract:
To study the complexity of hot Jupiter atmospheres revealed by observations of increasing quality, we have adapted the UK Met Office global circulation model (GCM), the Unified Model (UM), to these exoplanets. The UM solves the full 3D Navier-Stokes equations with a height-varying gravity, avoiding the simplifications used in most GCMs currently applied to exoplanets. In this work we present the c…
▽ More
To study the complexity of hot Jupiter atmospheres revealed by observations of increasing quality, we have adapted the UK Met Office global circulation model (GCM), the Unified Model (UM), to these exoplanets. The UM solves the full 3D Navier-Stokes equations with a height-varying gravity, avoiding the simplifications used in most GCMs currently applied to exoplanets. In this work we present the coupling of the UM dynamical core to an accurate radiation scheme based on the two-stream approximation and correlated-k method with state-of-the-art opacities from ExoMol. Our first application of this model is devoted to the extensively studied hot Jupiter HD 209458b. We have derived synthetic emission spectra and phase curves, and compare them to both previous models also based on state-of-the-art radiative transfer, and to observations. We find a reasonable a agreement between observations and both our day side emission and hot spot offset, however, our night side emission is too large. Overall our results are qualitatively similar to those found by Showman et al. (2009) with the SPARC/MITgcm, however, we note several quantitative differences: Our simulations show significant variation in the position of the hottest part of the atmosphere with pressure, as expected from simple timescale arguments, and in contrast to the "vertical coherency" found by Showman et al. (2009). We also see significant quantitative differences in calculated synthetic observations. Our comparisons strengthen the need for detailed intercomparisons of dynamical cores, radiation schemes and post-processing tools to understand these differences. This effort is necessary in order to make robust conclusions about these atmospheres based on GCM results.
△ Less
Submitted 9 September, 2016; v1 submitted 30 August, 2016;
originally announced August 2016.
-
The Effects of Consistent Chemical Kinetics Calculations on the Pressure-Temperature Profiles and Emission Spectra of Hot Jupiters
Authors:
Benjamin Drummond,
Pascal Tremblin,
Isabelle Baraffe,
David S. Amundsen,
Nathan J. Mayne,
Olivia Venot,
Jayesh Goyal
Abstract:
In this work we investigate the impact of calculating non-equilibrium chemical abundances consistently with the temperature structure for the atmospheres of highly-irradiated, close-in gas giant exoplanets. Chemical kinetics models have been widely used in the literature to investigate the chemical compositions of hot Jupiter atmospheres which are expected to be driven away from chemical equilibri…
▽ More
In this work we investigate the impact of calculating non-equilibrium chemical abundances consistently with the temperature structure for the atmospheres of highly-irradiated, close-in gas giant exoplanets. Chemical kinetics models have been widely used in the literature to investigate the chemical compositions of hot Jupiter atmospheres which are expected to be driven away from chemical equilibrium via processes such as vertical mixing and photochemistry. All of these models have so far used pressure--temperature (P-T) profiles as fixed model input. This results in a decoupling of the chemistry from the radiative and thermal properties of the atmosphere, despite the fact that in nature they are intricately linked. We use a one-dimensional radiative-convective equilibrium model, ATMO, which includes a sophisticated chemistry scheme to calculate P-T profiles which are fully consistent with non-equilibrium chemical abundances, including vertical mixing and photochemistry. Our primary conclusion is that, in cases of strong chemical disequilibrium, consistent calculations can lead to differences in the P-T profile of up to 100 K compared to the P-T profile derived assuming chemical equilibrium. This temperature change can, in turn, have important consequences for the chemical abundances themselves as well as for the simulated emission spectra. In particular, we find that performing the chemical kinetics calculation consistently can reduce the overall impact of non-equilibrium chemistry on the observable emission spectrum of hot Jupiters. Simulated observations derived from non-consistent models could thus yield the wrong interpretation. We show that this behaviour is due to the non-consistent models violating the energy budget balance of the atmosphere.
△ Less
Submitted 14 July, 2016;
originally announced July 2016.
-
Detection of H2O and evidence for TiO/VO in an ultra hot exoplanet atmosphere
Authors:
Thomas M. Evans,
David K. Sing,
Hannah R. Wakeford,
Nikolay Nikolov,
Gilda E. Ballester,
Benjamin Drummond,
Tiffany Kataria,
Neale P. Gibson,
David S. Amundsen,
Jessica Spake
Abstract:
We present a primary transit observation for the ultra hot (Teq~2400K) gas giant expolanet WASP-121b, made using the Hubble Space Telescope Wide Field Camera 3 in spectroscopic mode across the 1.12-1.64 micron wavelength range. The 1.4 micron water absorption band is detected at high confidence (5.4 sigma) in the planetary atmosphere. We also reanalyze ground-based photometric lightcurves taken in…
▽ More
We present a primary transit observation for the ultra hot (Teq~2400K) gas giant expolanet WASP-121b, made using the Hubble Space Telescope Wide Field Camera 3 in spectroscopic mode across the 1.12-1.64 micron wavelength range. The 1.4 micron water absorption band is detected at high confidence (5.4 sigma) in the planetary atmosphere. We also reanalyze ground-based photometric lightcurves taken in the B, r', and z' filters. Significantly deeper transits are measured in these optical bandpasses relative to the near-infrared wavelengths. We conclude that scattering by high-altitude haze alone is unlikely to account for this difference, and instead interpret it as evidence for titanium oxide and vanadium oxide absorption. Enhanced opacity is also inferred across the 1.12-1.3 micron wavelength range, possibly due to iron hydride absorption. If confirmed, WASP-121b will be the first exoplanet with titanium oxide, vanadium oxide, and iron hydride detected in transmission. The latter are important species in M/L dwarfs, and their presence is likely to have a significant effect on the overall physics and chemistry of the atmosphere, including the production of a strong thermal inversion.
△ Less
Submitted 8 April, 2016;
originally announced April 2016.
-
Cloudless atmospheres for L/T dwarfs and extra-solar giant planets
Authors:
P. Tremblin,
D. S. Amundsen,
G. Chabrier,
I. Baraffe,
B. Drummond,
S. Hinkley,
P. Mourier,
O. Venot
Abstract:
The admitted, conventional scenario to explain the complex spectral evolution of brown dwarfs (BD) since their first detections twenty years ago, has always been the key role played by micron-size condensates, called "dust" or "clouds", in their atmosphere. This scenario, however, faces major problems, in particular the J-band brightening and the resurgence of FeH absorption at the L to T transiti…
▽ More
The admitted, conventional scenario to explain the complex spectral evolution of brown dwarfs (BD) since their first detections twenty years ago, has always been the key role played by micron-size condensates, called "dust" or "clouds", in their atmosphere. This scenario, however, faces major problems, in particular the J-band brightening and the resurgence of FeH absorption at the L to T transition, and a physical first-principle understanding of this transition is lacking. In this paper, we propose a new, completely different explanation for BD and extrasolar giant planet (EGP) spectral evolution, without the need to invoke clouds. We show that, due to the slowness of the CO/CH4 and N2/NH3 chemical reactions, brown dwarf (L and T, respectively) and EGP atmospheres are subject to a thermo-chemical instability similar in nature to the fingering or chemical convective instability present in Earth oceans and at the Earth core/mantle boundary. The induced small-scale turbulent energy transport reduces the temperature gradient in the atmosphere, explaining the observed increase in near infrared J-H and J-K colors of L dwarfs and hot EGPs, while a warming up of the deep atmosphere along the L to T transition, as the CO/CH4 instability vanishes, naturally solves the two aforementioned puzzles, and provides a physical explanation of the L to T transition. This new picture leads to a drastic revision of our understanding of BD and EGP atmospheres and their evolution.
△ Less
Submitted 14 January, 2016;
originally announced January 2016.
-
Fingering convection and cloudless models for cool brown dwarf atmospheres
Authors:
P. Tremblin,
D. S. Amundsen,
P. Mourier,
I. Baraffe,
G. Chabrier,
B. Drummond,
D. Homeier,
O. Venot
Abstract:
This work aims to improve the current understanding of the atmospheres of brown dwarfs, especially cold ones with spectral type T and Y, whose modeling is a current challenge. Silicate and iron clouds are believed to disappear at the photosphere at the L/T transition, but cloudless models fail to reproduce correctly the spectra of T dwarfs, advocating for the addition of more physics, e.g. other t…
▽ More
This work aims to improve the current understanding of the atmospheres of brown dwarfs, especially cold ones with spectral type T and Y, whose modeling is a current challenge. Silicate and iron clouds are believed to disappear at the photosphere at the L/T transition, but cloudless models fail to reproduce correctly the spectra of T dwarfs, advocating for the addition of more physics, e.g. other types of clouds or internal energy transport mechanisms. We use a one-dimensional (1D) radiative/convective equilibrium code ATMO to investigate this issue. This code includes both equilibrium and out-of-equilibrium chemistry and solves consistently the PT structure. Included opacity sources are H2-H2, H2-He, H2O, CO, CO2, CH4, NH3, K, Na, and TiO, VO if they are present in the atmosphere. We show that the spectra of Y dwarfs can be accurately reproduced with a cloudless model if vertical mixing and NH3 quenching are taken into account. T dwarf spectra still have some reddening in e.g. J - H compared to cloudless models. This reddening can be reproduced by slightly reducing the temperature gradient in the atmosphere. We propose that this reduction of the stabilizing temperature gradient in these layers, leading to cooler structures, is due to the onset of fingering convection, triggered by the destabilizing impact of condensation of very thin dust.
△ Less
Submitted 13 April, 2015;
originally announced April 2015.