David S. Amundsen

Accuracy tests of radiation schemes used in hot Jupiter global circulation models

The treatment of radiation transport in global circulation models (GCMs) is crucial for correctly describing Earth and exoplanet atmospheric dynamics processes. The two-stream approximation and correlated-k method are currently state-of-the-art approximations applied in both Earth and hot Jupiter GCM radiation schemes to facilitate the rapid calculation of fluxes and heating rates. Their accuracy have been tested extensively for Earth-like conditions, but verification of the methods’ applicability to hot Jupiter-like conditions is lacking in the literature. We are adapting the UK Met O ce GCM, the Unified Model (UM), for the study of hot Jupiters, and present in this work the adaptation of the Edwards-Slingo radiation scheme based on the two-stream approximation and the correlated-k method. We discuss the calculation of absorption coe cients from high-temperature line lists and highlight the large uncertainty in the pressure-broadened line widths. We compare fluxes and heating rates obtained with our adapted scheme to more accurate discrete ordinate (DO) line-by-line (LbL) calculations ignoring scattering e ects. We find that, in most cases, errors stay below 10% for both heating rates and fluxes using 10 k-coe cients in each band and a di usivity factor D = 1:66. The two-stream approximation and the correlated-k method both contribute non-negligibly to the total error. We also find that using band-averaged absorption coe cients, which have previously been used in radiative-hydrodynamical simulations of a hot Jupiter, may yield errors of 100%, and should thus be used with caution.

Transition from Static to Kinetic Friction: Insights from a 2D Model

We describe a 2D spring-block model for the transition from static to kinetic friction at an elastic-slider-rigid-substrate interface obeying a minimalistic friction law (Amontons-Coulomb). By using realistic boundary conditions, a number of previously unexplained experimental results on precursory microslip fronts are successfully reproduced. From the analysis of the interfacial stresses, we derive a prediction for the evolution of the precursor length as a function of the applied loads, as well as an approximate relationship between microscopic and macroscopic friction coefficients. We show that the stress buildup due to both elastic loading and microslip-related relaxations depends only weakly on the underlying shear crack propagation dynamics. Conversely, crack speed depends strongly on both the instantaneous stresses and the friction coefficients, through a nontrivial scaling parameter.

The unified model, a fully-compressible, non-hydrostatic, deep atmosphere global circulation model, applied to hot Jupiters - ENDGame for a HD 209458b test case

We are adapting the global circulation model (GCM) of the UK Met Office, the so-called unified model (UM), for the study of hot Jupiters. In this work we demonstrate the successful adaptation of the most sophisticated dynamical core, the component of the GCM which solves the equations of motion for the atmosphere, available within the UM, ENDGame (Even Newer Dynamics for General atmospheric modelling of the environment). Within the same numerical scheme ENDGame supports solution to the dynamical equations under varying degrees of simplification. We present results from a simple, shallow (in atmospheric domain) hot Jupiter model (SHJ), and a more realistic (with a deeper atmosphere) HD 209458b test case. For both test cases we find that the large-scale, time-averaged (over the 1200 days prescribed test period), dynamical state of the atmosphere is relatively insensitive to the level of simplification of the dynamical equations. However, problems exist when attempting to reproduce the results for these test cases derived from other models. For the SHJ case the lower (and upper) boundary intersects the dominant dynamical features of the atmosphere meaning the results are heavily dependent on the boundary conditions. For the HD 209458b test case, when using the more complete dynamical models, the atmosphere is still clearly evolving after 1200 days, and in a transient state. Solving the complete (deep atmosphere and non-hydrostatic) dynamical equations allows exchange between the vertical and horizontal momentum of the atmosphere, via Coriolis and metric terms. Subsequently, interaction between the upper atmosphere and the deeper more slowly evolving (radiatively inactive) atmosphere significantly alters the results, and acts over timescales longer than 1200 days.

A uniform analysis of HD 209458b Spitzer/IRAC light curves with Gaussian process models

We present an analysis of Spitzer/IRAC primary transit and secondary eclipse lightcurves measured for HD 209458b, using Gaussian process models to marginalise over the intrapixel sensitivity variations in the 3:6 m and 4:5 m channels and the ramp eect in the 5 :8 m and 8:0 m channels. The main advantage of this approach is that we can account for a broad range of degeneracies between the planet signal and systematics without actually having to specify a deterministic functional form for the latter. Our results do not conrm a previous claim of water absorption in transmission. Instead, our results are more consistent with a featureless transmission spectrum, possibly due to a cloud deck obscuring molecular absorption bands. For the emission data, our values are not consistent with the thermal inversion in the dayside atmosphere that was originally inferred from these data. Instead, we agree with another re-analysis of these same data, which concluded a non-inverted atmosphere provides a better t. We

Detection of H2O and Evidence for TiO VO in an Ultra Hot Exoplanet Atmosphere.

The authors would like to thank the referee for their prompt and thoughtful review. This work is based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. These observations are associated with program GO-14468. The authors are grateful to the WASP-121 discovery team for generously providing the ground-based photometric light curves. The research leading to these results received funding from the European Research Council under the European Union Seventh Framework Program (FP7/2007-2013) ERC grant agreement no. 336792. H.R.W. acknowledges support by an appointment to the NASA Postdoctoral Program at Goddard Space Flight Center, administered by ORAU and USRA through a contract with NASA. N.P.G. gratefully acknowledges support from the Royal Society in the form of a University Research Fellowship.

Fingering convection and cloudless models for cool brown dwarf atmospheres

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.

The mineral clouds on HD 209458b and HD 189733b

We highlight financial support of the European Community under the FP7 by the ERC starting grant 257431 and by an ERC advanced grant 247060. JK acknowledges the Rosen fellowship from the Brooklyn College New York, US.

An ultrahot gas-giant exoplanet with a stratosphere

Infrared radiation emitted from a planet contains information about the chemical composition and vertical temperature profile of its atmosphere. If upper layers are cooler than lower layers, molecular gases will produce absorption features in the planetary thermal spectrum. Conversely, if there is a stratosphere—where temperature increases with altitude—these molecular features will be observed in emission. It has been suggested that stratospheres could form in highly irradiated exoplanets, but the extent to which this occurs is unresolved both theoretically and observationally. A previous claim for the presence of a stratosphere remains open to question, owing to the challenges posed by the highly variable host star and the low spectral resolution of the measurements. Here we report a near-infrared thermal spectrum for the ultrahot gas giant WASP-121b, which has an equilibrium temperature of approximately 2,500 kelvin. Water is resolved in emission, providing a detection of an exoplanet stratosphere at 5σ confidence. These observations imply that a substantial fraction of incident stellar radiation is retained at high altitudes in the atmosphere, possibly by absorbing chemical species such as gaseous vanadium oxide and titanium oxide.

Treatment of overlapping gaseous absorption with the correlated-k method in hot Jupiter and brown dwarf atmosphere models

The correlated-k method is frequently used to speed up radiation calculations in both one-dimensional and three-dimensional atmosphere models. An inherent difficulty with this method is how to treat overlapping absorption, i.e. absorption by more than one gas in a given spectral region. We have evaluated the applicability of three different methods in hot Jupiter and brown dwarf atmosphere models, all of which have been previously applied within models in the literature: (i) Random overlap, both with and without resorting and rebinning, (ii) equivalent extinction and (iii) pre-mixing of opacities, where (i) and (ii) combine k-coefficients for different gases to obtain k-coefficients for a mixture of gases, while (iii) calculates k-coefficients for a given mixture from the corresponding mixed line-by-line opacities. We find that the random overlap method is the most accurate and flexible of these treatments, and is fast enough to be used in one-dimensional models with resorting and rebinning. In three-dimensional models such as GCMs it is too slow, however, and equivalent extinction can provide a speed-up of at least a factor of three with only a minor loss of accuracy while at the same time retaining the flexibility gained by combining k-coefficients computed for each gas individually. Pre-mixed opacities are significantly less flexible, and we also find that particular care must be taken when using this method in order to to adequately resolve steep variations in composition at important chemical equilibrium boundaries. We use the random overlap method with resorting and rebinning in our one-dimensional atmosphere model and equivalent extinction in our GCM, which allows us to e.g. consistently treat the feedback of non-equilibrium chemistry on the total opacity and therefore the calculated P-T profiles in our models.

HAT-P-26b: A Neptune-mass exoplanet with a well-constrained heavy element abundance

How much water is in that exoplanet? Thousands of exoplanets have been identified, but little is known about their atmospheres, especially for bodies smaller than Jupiter. The extent and composition of an atmosphere can provide evidence for how an exoplanet formed. Wakeford et al. used the Hubble and Spitzer space telescopes to measure the spectrum of the atmosphere around HAT-P-26b, a Neptune-sized exoplanet discovered in 2011. They detected signatures of water and clouds; this allowed them to constrain the atmospheric composition, which appears not to have been altered substantially since it formed. Science, this issue p. 628 The Neptune-sized exoplanet HAT-P-26b has an atmosphere containing water and clouds. A correlation between giant-planet mass and atmospheric heavy elemental abundance was first noted in the past century from observations of planets in our own Solar System and has served as a cornerstone of planet-formation theory. Using data from the Hubble and Spitzer Space Telescopes from 0.5 to 5 micrometers, we conducted a detailed atmospheric study of the transiting Neptune-mass exoplanet HAT-P-26b. We detected prominent H2O absorption bands with a maximum base-to-peak amplitude of 525 parts per million in the transmission spectrum. Using the water abundance as a proxy for metallicity, we measured HAT-P-26b’s atmospheric heavy element content (4.8−4.0+21.5 times solar). This likely indicates that HAT-P-26b’s atmosphere is primordial and obtained its gaseous envelope late in its disk lifetime, with little contamination from metal-rich planetesimals.