Raymond Pierrehumbert FRS

Abstract:

The earliest atmospheres of rocky planets originate from extensive volatile release during magma ocean epochs that occur during assembly of the planet. These establish the initial distribution of the major volatile elements between different chemical reservoirs that subsequently evolve via geological cycles. Current theoretical techniques are limited in exploring the anticipated range of compositional and thermal scenarios of early planetary evolution, even though these are of prime importance to aid astronomical inferences on the environmental context and geological history of extrasolar planets. Here, we present a coupled numerical framework that links an evolutionary, vertically‐resolved model of the planetary silicate mantle with a radiative‐convective model of the atmosphere. Using this method we investigate the early evolution of idealized Earth‐sized rocky planets with end‐member, clear‐sky atmospheres dominated by either H2, H2O, CO2, CH4, CO, O2, or N2. We find central metrics of early planetary evolution, such as energy gradient, sequence of mantle solidification, surface pressure, or vertical stratification of the atmosphere, to be intimately controlled by the dominant volatile and outgassing history of the planet. Thermal sequences fall into three general classes with increasing cooling timescale: CO, N2, and O2 with minimal effect, H2O, CO2, and CH4 with intermediate influence, and H2 with several orders of magnitude increase in solidification time and atmosphere vertical stratification. Our numerical experiments exemplify the capabilities of the presented modeling framework and link the interior and atmospheric evolution of rocky exoplanets with multi‐wavelength astronomical observations.

Abstract:

White dwarf–brown dwarf short-period binaries (Porb ≲ 2 h) are some of the most extreme irradiated atmospheric environments known. These systems offer an opportunity to explore theoretical and modelling efforts of irradiated atmospheres different to typical hot Jupiter systems. We aim to investigate the three-dimensional (3D) atmospheric structural and dynamical properties of the brown dwarf WD 0137−349B. We use the 3D global circulation model (GCM) Exo-Flexible Modelling System (FMS) with a dual-band grey radiative transfer scheme to model the atmosphere of WD 0137−349B. The results of the GCM model are post-processed using the 3D Monte Carlo radiative transfer model CMCRT. Our results suggest inefficient day–night energy transport and a large day–night temperature contrast for WD 0137−349B. Multiple flow patterns are present, shifting energy asymmetrically eastward or westward depending on their zonal direction and latitude. Regions of overturning are produced on the western terminator. We are able to reproduce the start of the system near-infrared (IR) emission excess at ≳1.95 μm as observed by the Gemini Near-Infrared Spectrograph (GNIRS) instrument. Our model overpredicts the IR phase curve fluxes by factors of ≈1–3, but generally fits the shape of the phase curves well. Chemical kinetic modelling using VULCAN suggests a highly ionized region at high altitudes can form on the dayside of the brown dwarf. We present a first attempt at simulating the atmosphere of a short-period white dwarf–brown dwarf binary in a 3D setting. Further studies into the radiative and photochemical heating from the ultraviolet irradiation are required to more accurately capture the energy balance inside the brown dwarf atmosphere. Cloud formation may also play an important role in shaping the emission spectra of the brown dwarf.