Search

Context. Photochemistry is a key process driving planetary atmospheres away from local thermodynamic equilibrium. Recent observations of the H$_2$ dominated atmospheres of hot gas giants have detected SO$_2$ as one of the major products of this process. Aims. We investigate which chemical pathways lead to the formation of SO$_2$ in an atmosphere, and we investigate which part of the flux from the host star is necessary to initiate SO$_2$ production. Methods. We use the publicly available S-N-C-H-O photochemical network in the VULCAN chemical kinetics code to compute the disequilibrium chemistry of an exoplanetary atmosphere. Results. We find that there are two distinct chemical pathways that lead to the formation of SO$_2$. The formation of SO$_2$ at higher pressures is initiated by stellar flux >200 nm, whereas the formation of SO$_2$ at lower pressures is initiated by stellar flux <200 nm. In deeper layers of the atmosphere, OH is provided by the hydrogen abstraction of H$_2$O, and sulfur is provided by the photodissociation of SH and S$_2$, which leads to a positive feedback cycle that liberates sulfur from the stable H$_2$S molecule. In higher layers of the atmosphere, OH is provided by the photodissociation of H$_2$O, and sulfur can be liberated from H$_2$S by either photodissociation of SH and S$_2$, or by the hydrogen abstraction of SH. Conclusions. We conclude that the stellar flux in the 200-350 nm wavelength range as well as the ratio of NUV/UV radiation are important parameters determining the observability of SO$_2$. In addition we find that there is a diversity of chemical pathways to the formation of SO$_2$. This is crucial for the interpretation of SO$_2$ detections and derived elemental abundance ratios and overall metallicities., Comment: Accepted for publication in Astronomy & Astrophysics

We present a transmission spectrum of the misaligned hot Jupiter WASP-15b from 2.8--5.2 microns observed with JWST's NIRSpec/G395H grating. Our high signal to noise data, which has negligible red noise, reveals significant absorption by H$_2$O ($4.2\sigma$) and CO$_2$ ($8.9\sigma$). From independent data reduction and atmospheric retrieval approaches, we infer that WASP-15b's atmospheric metallicity is super-solar ($\gtrsim 15\times$ solar) and its C/O is consistent with solar, that together imply planetesimal accretion. Our GCM simulations for WASP-15b suggest that the C/O we measure at the limb is likely representative of the entire photosphere due to the mostly uniform spatial distribution of H$_2$O, CO$_2$ and CO. We additionally see evidence for absorption by SO$_2$ and absorption at 4.9$\mu$m, for which the current leading candidate is OCS, albeit with several caveats. If confirmed, this would be the first detection of OCS in an exoplanet atmosphere and point towards complex photochemistry of sulphur-bearing species in the upper atmosphere. These are the first observations from the BOWIE-ALIGN survey which is using JWST's NIRSpec/G395H instrument to compare the atmospheric compositions of aligned/low-obliquity and misaligned/high-obliquity hot Jupiters around F stars above the Kraft break. The goal of our survey is to determine whether the atmospheric composition differs across two populations of planets that have likely undergone different migration histories (disc versus disc-free) as evidenced by their obliquities (aligned versus misaligned)., Comment: 24 pages, 23 figures, 6 tables. Submitted to MNRAS

The characterization of young planets (< 300 Myr) is pivotal for understanding planet formation and evolution. We present the 3-5$\mu$m transmission spectrum of the 17 Myr, Jupiter-size ($R$ $\sim$10$R_{\oplus}$) planet, HIP 67522 b, observed with JWST/NIRSpec/G395H. To check for spot contamination, we obtain a simultaneous $g$-band transit with SOAR. The spectrum exhibits absorption features 30-50% deeper than the overall depth, far larger than expected from an equivalent mature planet, and suggests that HIP 67522 b's mass is $<$20 $M_{\oplus}$ irrespective of cloud cover and stellar contamination. A Bayesian retrieval analysis returns a mass constraint of $13.8\pm1.0M_{\oplus}$. This challenges the previous classification of HIP 67522 b as a hot Jupiter and instead, positions it as a precursor to the more common sub-Neptunes. With a density of $<$0.10g/cm$^{3}$, HIP 67522 b is one of the lowest density planets known. We find strong absorption from H$_{2}$O and CO$_{2}$ ($\ge7\sigma$), a modest detection of CO (3.5$\sigma$), and weak detections of H$_2$S and SO$_2$ ($\simeq2\sigma$). Comparisons with radiative-convective equilibrium models suggest supersolar atmospheric metallicities and solar-to-subsolar C/O ratios, with photochemistry further constraining the inferred atmospheric metallicity to 3$\times$10 Solar due to the amplitude of the SO$_2$ feature. These results point to the formation of HIP 67522 b beyond the water snowline, where its envelope was polluted by icy pebbles and planetesimals. The planet is likely experiencing substantial mass loss (0.01-0.03 M$_{\oplus}$ Myr$^{-1}$), sufficient for envelope destruction within a Gyr. This highlights the dramatic evolution occurring within the first 100 Myr of its existence., Comment: Accepted for publication in The Astronomical Journal; 32 pages, 18 figures, 7 tables

Due to the detection of phosphine PH3 in the Solar System gas giants Jupiter and Saturn, PH3 has long been suggested to be detectable in exosolar substellar atmospheres too. However, to date, a direct detection of phosphine has proven to be elusive in exoplanet atmosphere surveys. We construct an updated phosphorus-hydrogen-oxygen (PHO) photochemical network suitable for simulation of gas giant hydrogen-dominated atmospheres. Using this network, we examine PHO photochemistry in hot Jupiter and warm Neptune exoplanet atmospheres at Solar and enriched metallicities. Our results show for HD 189733b-like hot Jupiters that HOPO, PO and P2 are typically the dominant P carriers at pressures important for transit and emission spectra, rather than PH3. For GJ1214b-like warm Neptune atmospheres our results suggest that at Solar metallicity PH3 is dominant in the absence of photochemistry, but is generally not in high abundance for all other chemical environments. At 10 and 100 times Solar, small oxygenated phosphorus molecules such as HOPO and PO dominate for both thermochemical and photochemical simulations. The network is able to reproduce well the observed PH3 abundances on Jupiter and Saturn. Despite progress in improving the accuracy of the PHO network, large portions of the reaction rate data remain with approximate, uncertain or missing values, which could change the conclusions of the current study significantly. Improving understanding of the kinetics of phosphorus-bearing chemical reactions will be a key undertaking for astronomers aiming to detect phosphine and other phosphorus species in both rocky and gaseous exoplanetary atmospheres in the near future., Comment: Submitted to ApJ (12 July 2024) - Accepted ApJ (20 Oct 2024)

Temperate exoplanets between the sizes of Earth and Neptune, known as "sub-Neptunes", have emerged as intriguing targets for astrobiology. It is unknown whether these planets resemble Earth-like terrestrial worlds with a habitable surface, Neptune-like giant planets with deep atmospheres and no habitable surface, or something exotic in between. Recent JWST transmission spectroscopy observations of the canonical sub-Neptune K2-18 b revealed ~1% CH4, ~1% CO2, and a non-detection of CO in the atmosphere. While previous studies have proposed that the observed atmospheric composition could help constrain the lower atmosphere conditions and determine the interior structure of sub-Neptunes like K2-18 b, the possible interactions between the atmosphere and a hot, supercritical water ocean at its base remain unexplored. In this work, we investigate whether a global supercritical water ocean, resembling a planetary-scale hydrothermal system, can explain these observations on K2-18 b-like sub-Neptunes through equilibrium aqueous geochemical calculations. We find that the observed atmospheric CH4/CO2 ratio implies a minimum ocean temperature of ~715 K, whereas the corresponding CO/CO2 ratio allows ocean temperatures up to ~1060 K. These results indicate that a global supercritical water ocean on K2-18 b is plausible. While life cannot survive in this ocean, this work represents the first step towards understanding how a global supercritical water ocean may influence observable atmospheric characteristics on volatile-rich sub-Neptunes. Future observations with better constrained NH3 and CO mixing ratios could further help distinguish between possible interior compositions of K2-18 b., Comment: 15 pages, 5 figures, 1 table

Many super-Earths are on very short orbits around their host star and, therefore, more likely to be tidally locked. Because this locking can lead to a strong contrast between the dayside and nightside surface temperatures, these super-Earths could exhibit mantle convection patterns and tectonics that could differ significantly from those observed in the present-day solar system. The presence of an atmosphere, however, would allow transport of heat from the dayside towards the nightside and thereby reduce the surface temperature contrast between the two hemispheres. On rocky planets, atmospheric and geodynamic regimes are closely linked, which directly connects the question of atmospheric thickness to the potential interior dynamics of the planet. Here, we study the interior dynamics of super-Earth GJ 486b ($R=1.34$ $R_{\oplus}$, $M=3.0$ $M_{\oplus}$, T$_\mathrm{eq}\approx700$ K), which is one of the most suitable M-dwarf super-Earth candidates for retaining an atmosphere produced by degassing from the mantle and magma ocean. We investigate how the geodynamic regime of GJ 486b is influenced by different surface temperature contrasts by varying possible atmospheric circulation regimes. We also investigate how the strength of the lithosphere affects the convection pattern. We find that hemispheric tectonics, the surface expression of degree-1 convection with downwellings forming on one hemisphere and upwelling material rising on the opposite hemisphere, is a consequence of the strong lithosphere rather than surface temperature contrast. Anchored hemispheric tectonics, where downwellings und upwellings have a preferred (day/night) hemisphere, is favoured for strong temperature contrasts between the dayside and nightside and higher surface temperatures., Comment: 36 pages, 13 figures, resubmitted to JGR: Planets

Transmission spectroscopy has been a workhorse technique over the past two decades to constrain the physical and chemical properties of exoplanet atmospheres. One of its classical key assumptions is that the portion of the atmosphere it probes -- the terminator region -- is homogeneous. Several works in the past decade, however, have put this into question for highly irradiated, hot ($T_{eq}\gtrsim 1000$ K) gas giant exoplanets both empirically and via 3-dimensional modelling. While models predict clear differences between the evening (day-to-night) and morning (night-to-day) terminators, direct morning/evening transmission spectra in a wide wavelength range has not been reported for an exoplanet to date. Under the assumption of precise and accurate orbital parameters on WASP-39 b, here we report the detection of inhomogeneous terminators on the exoplanet WASP-39 b, which allows us to retrieve its morning and evening transmission spectra in the near-infrared ($2-5\ \mu$m) using JWST. We observe larger transit depths in the evening which are, on average, $405 \pm 88$ ppm larger than the morning ones, also having qualitatively larger features than the morning spectrum. The spectra are best explained by models in which the evening terminator is hotter than the morning terminator by $177^{+65}_{-57}$ K with both terminators having C/O ratios consistent with solar. General circulation models (GCMs) predict temperature differences broadly consistent with the above value and point towards a cloudy morning terminator and a clearer evening terminator., Comment: Published in Nature at https://www.nature.com/articles/s41586-024-07768-4. All code to produce plots (with data) can be found at https://github.com/nespinoza/wasp39-terminators

Understanding the surface temperature and interior structure of cold-to-temperate sub-Neptunes is critical for assessing their habitability, yet direct observations are challenging. In this study, we investigate the impact of water condensation on the atmospheric compositions of sub-Neptunes, focusing on the implications for JWST spectroscopic observations. By modeling the atmospheric photochemistry of two canonical sub-Neptunes, K2-18 b and LHS 1140 b, both with and without water condensation and with and without thick atmospheres, we demonstrate that water condensation can significantly affect the predicted atmospheric compositions. This effect is driven by oxygen depletion from the condensation of water vapor and primarily manifests as an increase in the C/O ratio within the photochemically active regions of the atmosphere. This change in composition particularly affects planets with thin H2-dominated atmospheres, leading to a transition in dominant nitrogen and carbon carriers from N2 and oxygen-rich species like CO/CO2 towards heavier hydrocarbons and nitriles. While our models do not fully account for the loss mechanisms of these higher-order species, such molecules can go on to form more refractory molecules or hazes. Planets with thin H2-rich atmospheres undergoing significant water condensation are thus likely to exhibit very hazy atmospheres. The relatively flat JWST spectra observed for LHS 1140 b could be consistent with such a scenario, suggesting a shallow surface with extensive water condensation or a high atmospheric C/O ratio. Conversely, the JWST observations of K2-18 b are better aligned with a volatile-rich mini-Neptune with a thick atmosphere., Comment: 21 pages, 7 figures

The recent inference of sulphur dioxide (SO$_2$) in the atmosphere of the hot ($\sim$1100 K), Saturn-mass exoplanet WASP-39b from near-infrared JWST observations suggests that photochemistry is a key process in high temperature exoplanet atmospheres. This is due to the low ($<$1 ppb) abundance of SO$_2$ under thermochemical equilibrium, compared to that produced from the photochemistry of H$_2$O and H$_2$S (1-10 ppm). However, the SO$_2$ inference was made from a single, small molecular feature in the transmission spectrum of WASP-39b at 4.05 $\mu$m, and therefore the detection of other SO$_2$ absorption bands at different wavelengths is needed to better constrain the SO$_2$ abundance. Here we report the detection of SO$_2$ spectral features at 7.7 and 8.5 $\mu$m in the 5-12 $\mu$m transmission spectrum of WASP-39b measured by the JWST Mid-Infrared Instrument (MIRI) Low Resolution Spectrometer (LRS). Our observations suggest an abundance of SO$_2$ of 0.5-25 ppm (1$\sigma$ range), consistent with previous findings. In addition to SO$_2$, we find broad water vapour absorption features, as well as an unexplained decrease in the transit depth at wavelengths longer than 10 $\mu$m. Fitting the spectrum with a grid of atmospheric forward models, we derive an atmospheric heavy element content (metallicity) for WASP-39b of $\sim$7.1-8.0 $\times$ solar and demonstrate that photochemistry shapes the spectra of WASP-39b across a broad wavelength range., Comment: Published in Nature

Magnetic fields connect an array of planetary processes, from atmospheric escape to interior convection. Despite their importance, exoplanet magnetic fields are largely unconstrained by both theory and observation. In this Letter, we propose a novel method for constraining the B field strength of hot gas giants: comparing the velocities of heavy ions and neutral gas with high-resolution spectroscopy. The core concept of this method is that ions are directly deflected by magnetic fields. While neutrals are also affected by B fields via friction with field-accelerated ions, ionic gas should be more strongly coupled to the underlying magnetic field than bulk neutral flow. Hence, measuring the difference between the two velocities yields rough constraints on the B field, provided an estimate of the stellar UV flux is known. We demonstrate that heavy ions are particularly well suited for this technique, because they are less likely to be entrained in complex hydrodynamic outflows than their lighter counterparts. We perform a proof-of-concept calculation with Ba II, an ion whose velocity has been repeatedly measured at high confidence with high-resolution spectroscopy. Our work shows that a 10G magnetic field would produce ~ km/s ion--neutral velocity differences at a microbar, whereas a 50G magnetic field would produce ~20km/s velocity difference. With new leverage on magnetic fields, we will be able to investigate magnetic field generation in the extreme edge cases of hot gas giants, with wide-ranging consequences for planetary interior structure, dynamo theory, and habitability., Comment: Accepted at ApJ Letters. 12 pages; 4 figures; 2 tables

Ultra-hot Jupiters are among the best targets for atmospheric characterization at high spectral resolution. Resolving their transmission spectra as a function of orbital phase offers a unique window into the 3D nature of these objects. In this work, we present three transits of the ultra-hot Jupiter WASP-121b observed with Gemini-S/IGRINS. For the first time, we measure the phase-dependent absorption signals of CO and H$_{\text{2}}$O in the atmosphere of an exoplanet, and we find that they are different. While the blueshift of CO increases during the transit, the absorption lines of H$_{\text{2}}$O become less blueshifted with phase, and even show a redshift in the second half of the transit. These measurements reveal the distinct spatial distributions of both molecules across the atmospheres of ultra-hot Jupiters. Also, we find that the H$_{\text{2}}$O signal is absent in the first quarter of the transit, potentially hinting at cloud formation on the evening terminator of WASP-121b. To further interpret the absorption trails of CO and H$_{\text{2}}$O, as well as the Doppler shifts of Fe previously measured with VLT/ESPRESSO, we compare the data to simulated transits of WASP-121b. To this end, we post-processes the outputs of global circulation models with a 3D Monte-Carlo radiative transfer code. Our analysis shows that the atmosphere of WASP-121b is subject to atmospheric drag, as previously suggested by small hotspot offsets inferred from phase-curve observations. Our study highlights the importance of phase-resolved spectroscopy in unravelling the complex atmospheric structure of ultra-hot Jupiters and sets the stage for further investigations into their chemistry and dynamics., Comment: 24 pages, 16 figures, accepted for publication in PASP

Spectroscopic phase curves of hot Jupiters measure their emission spectra at multiple orbital phases, thus enabling detailed characterisation of their atmospheres. Precise constraints on the atmospheric composition of these exoplanets offer insights into their formation and evolution. We analyse four phase-resolved emission spectra of the hot Jupiter WASP-43b, generated from a phase curve observed with the MIRI/LRS onboard the JWST, to retrieve its atmospheric properties. Using a parametric 2D temperature model and assuming a chemically homogeneous atmosphere within the observed pressure region, we simultaneously fit the four spectra to constrain the abundances of atmospheric constituents, thereby yielding more precise constraints than previous work that analysed each spectrum independently. Our analysis reveals statistically significant evidence of NH3 (4$\sigma$) in a hot Jupiter's emission spectra for the first time, along with evidence of H2O (6.5$\sigma$), CO (3.1$\sigma$), and a non-detection of CH4. With our abundance constraints, we tentatively estimate the metallicity of WASP-43b at 0.6-6.5$\times$solar and its C/O ratio at 0.6-0.9. Our findings offer vital insights into the atmospheric conditions and formation history of WASP-43b by simultaneously constraining the abundances of carbon, oxygen, and nitrogen-bearing species., Comment: 17 pages, 14 figures, accepted for publication in MNRAS. Comments welcome!

The signal from a transiting planet can be diluted by astrophysical contamination. In the case of circumstellar debris disks, this contamination could start in the mid-infrared and vary as a function of wavelength, which would then change the observed transmission spectrum for any planet in the system. The MIRI/LRS WASP-39b transmission spectrum shows an unexplained dip starting at $\sim$10 $\mu$m that could be caused by astrophysical contamination. The spectral energy distribution displays excess flux at similar levels to that which are needed to create the dip in the transmission spectrum. In this article, we show that this dip is consistent with the presence of a bright circumstellar debris disk, at a distance of $>$2 au. We discuss how a circumstellar debris disk like that could affect the atmosphere of WASP-39b. We also show that even faint debris disks can be a source of contamination in MIRI exoplanet spectra., Comment: accepted to ApJL

Theoretical predictions and observational data indicate a class of sub-Neptune exoplanets may have water-rich interiors covered by hydrogen-dominated atmospheres. Provided suitable climate conditions, such planets could host surface liquid oceans. Motivated by recent JWST observations of K2-18 b, we self-consistently model the photochemistry and potential detectability of biogenic sulfur gases in the atmospheres of temperate sub-Neptune waterworlds for the first time. On Earth today, organic sulfur compounds produced by marine biota are rapidly destroyed by photochemical processes before they can accumulate to significant levels. Domagal-Goldman et al. (2011) suggest that detectable biogenic sulfur signatures could emerge in Archean-like atmospheres with higher biological production or low UV flux. In this study, we explore biogenic sulfur across a wide range of biological fluxes and stellar UV environments. Critically, the main photochemical sinks are absent on the nightside of tidally locked planets. To address this, we further perform experiments with a 3D GCM and a 2D photochemical model (VULCAN 2D (Tsai et al. 2024)) to simulate the global distribution of biogenic gases to investigate their terminator concentrations as seen via transmission spectroscopy. Our models indicate that biogenic sulfur gases can rise to potentially detectable levels on hydrogen-rich waterworlds, but only for enhanced global biosulfur flux ($\gtrsim$20 times modern Earth's flux). We find that it is challenging to identify DMS at 3.4 $\mu m$ where it strongly overlaps with CH$_4$, whereas it is more plausible to detect DMS and companion byproducts, ethylene (C$_2$H$_4$) and ethane (C$_2$H$_6$), in the mid-infrared between 9 and 13 $\mu m$., Comment: 9 pages, 4 figures, accepted for publication in ApJL

JWST recently measured the transmission spectrum of K2-18b, a habitable-zone sub-Neptune exoplanet, detecting CH$_4$ and CO$_2$ in its atmosphere. The discovery paper argued the data are best explained by a habitable "Hycean" world, consisting of a relatively thin H$_2$-dominated atmosphere overlying a liquid water ocean. Here, we use photochemical and climate models to simulate K2-18b as both a Hycean planet and a gas-rich mini-Neptune with no defined surface. We find that a lifeless Hycean world is hard to reconcile with the JWST observations because photochemistry only supports $< 1$ part-per-million CH$_4$ in such an atmosphere while the data suggest about $\sim 1\%$ of the gas is present. Sustaining %-level CH$_4$ on a Hycean K2-18b may require the presence of a methane-producing biosphere, similar to microbial life on Earth $\sim 3$ billion years ago. On the other hand, we predict that a gas-rich mini-Neptune with $100 \times$ solar metallicity should have 4% CH$_4$ and nearly 0.1% CO$_2$, which are compatible with the JWST data. The CH$_4$ and CO$_2$ are produced thermochemically in the deep atmosphere and mixed upward to the low pressures sensitive to transmission spectroscopy. The model predicts H$_2$O, NH$_3$ and CO abundances broadly consistent with the non-detections. Given the additional obstacles to maintaining a stable temperate climate on Hycean worlds due to H$_2$ escape and potential supercriticality at depth, we favor the mini-Neptune interpretation because of its relative simplicity and because it does not need a biosphere or other unknown source of methane to explain the data., Comment: Accepted for publication at ApJL

With JWST slated to gain high fidelity time dependent data on brown dwarf atmospheres, it is highly anticipated to do the same for directly imaged, sub-Jupiter exoplanets. With this new capability, the need for a full 3D understanding to explain spectral features and their time dependence is becoming a vital aspect for consideration. To examine the atmospheric properties of directly imaged sub-Jupiter exoplanets, we use the three dimensional Exo-FMS general circulation model (GCM) to simulate a metal enhanced generic young sub-Jupiter object. We couple Exo-FMS to a kinetic chemistry scheme, a tracer based cloud formation scheme and a spectral radiative-transfer model to take into account the chemical and cloud feedback on the atmospheric thermochemical and dynamical properties. Our results show a highly complex feedback between clouds and chemistry onto the 3D temperature structure of the atmosphere, bringing about latitudinal differences and inducing time-dependent stormy features at photospheric pressures. This suggests a strong connection and feedback between the spatial cloud coverage and chemical composition of the atmosphere, with the temperature changes and dynamical motions induced by cloud opacity and triggered convection feedback driving chemical species behaviour. In addition, we also produce synthetic latitude dependent and time dependent spectra of our model to investigate atmospheric variability and periodicity in commonly used photometric bands. Overall, our efforts put the included physics in 3D simulations of exoplanets on par with contemporary 1D radiative-convective equilibrium modelling., Comment: Accepted MNRAS (Feb 2024), 17 pages

The atmospheric dynamics of tidally-locked hot Jupiters is characterized by strong equatorial winds. Understanding the interaction between global circulation and chemistry is crucial in atmospheric studies and interpreting observations. Two-dimensional (2D) photochemical transport models shed light on how the atmospheric composition depends on circulation. In this paper, we introduce the 2D photochemical (horizontal and vertical) transport model, VULCAN 2D, which improves on the pseudo-2D approaches by allowing for non-uniform zonal winds. We extensively validate our VULCAN 2D with analytical solutions and benchmark comparisons. Applications to HD 189733 b and HD 209458 b reveal a transition in mixing regimes: horizontal transport predominates below $\sim$0.1 mbar while vertical mixing is more important at higher altitudes above 0.1 mbar. Motivated by the previously inferred carbon-rich atmosphere, we find that HD 209458 b with super-solar carbon-to-oxygen ratio (C/O) exhibits pronounced C$_2$H$_4$ absorption on the morning limb but not on the evening limb, owing to horizontal transport from the nightside. We discuss when a pseudo-2D approach is a valid assumption and its inherent limitations. Finally, we demonstrate the effect of horizontal transport in transmission observations and its impact on the morning-evening limb asymmetry with synthetic spectra, highlighting the need to consider global transport when interpreting exoplanet atmospheres., Comment: 19 pages, 20 figures, published in ApJ

The atmospheres of brown dwarfs have been long observed to exhibit a multitude of non-equilibrium chemical signatures and spectral variability across the L, T and Y spectral types. We aim to investigate the link between the large-scale 3D atmospheric dynamics and time-dependent chemistry in the brown dwarf regime, and to assess its impact on spectral variability. We couple the miniature kinetic chemistry module `mini-chem' to the Exo-FMS general circulation model (GCM). We then perform a series of idealised brown dwarf regime atmospheric models to investigate the dynamical 3D chemical structures produced by our simulations. The GCM output is post-processed using a 3D radiative-transfer model to investigate hemisphere-dependent spectral signatures and rotational variability. Our results show the expected strong non-equilibrium chemical behaviour brought on by vertical mixing as well as global spacial variations due to zonal flows. Chemical species are generally globally homogenised, showing variations of $\pm$10\% or less, dependent on pressure level, and follow the dynamical structures present in the atmosphere. However, we find localised storm regions and eddies can show higher contrasts, up to $\pm$100\%, in mixing ratio compared to the background global mean. This initial study represents another step in understanding the connection between three-dimensional atmospheric flows in brown dwarfs and their rich chemical inventories., Comment: MNRAS Accepted: 5 June 2023

JWST has recently detected the first robust photochemical product on an exoplanet: sulfur dioxide (SO$_2$) on WASP-39b (Rustamkulov et al. 2023; Alderson et al. 2023; Tsai et al. 2023b). The data from the NIRISS instrument also reveal signs of partial coverage of clouds (Feinstein et al. 2023). Most of the previous studies have focused on interpreting spectral data with 1D models. To explore how the chemical species and cloud particles are altered by global circulation, we applied a 2D photochemical model and a 2D microphysical cloud model separately to post-process the thermal and dynamical structures simulated by a 3D general circulation model (GCM) of WASP-39b. We found that SO$_2$ produced by photochemistry on the dayside can be transported to the nightside owing to the efficient replenishment of horizontal transport. The morning-evening limb differences in methane (CH$_4$) abundances predicted by the 1D models disappeared after horizontal transport is included. Similarly, the inclusion of horizontal transport also reduced the limb differences in SO$_2$. Our modeling results suggest that the fast zonal wind results in minimal or negligible limb asymmetry in composition. Based on the synthetic spectra generated by our 2D atmosphere simulations, we propose that observing SO$_2$ absorption in the emission spectra of WASP-39b at different phases may offer opportunities to probe the horizontal quenching process of photochemical products. We will focus on the gas-phase chemistry in this paper and leave the results regarding clouds in the subsequent paper as part of the series., Comment: 10 pages, 6 figures, published in ApJL

Hycean worlds are a proposed subset of sub-Neptune exoplanets with substantial water inventories, liquid surface oceans and extended hydrogen-dominated atmospheres that could be favourable for habitability. In this work, we aim to quantitatively define the inner edge of the Hycean habitable zone using a 1D radiative-convective model. As a limiting case, we model a dry hydrogen-helium envelope above a surface ocean. We find that 10 to 20 bars of atmosphere produces enough greenhouse effect to drive a liquid surface ocean supercritical when forced with current Earth-like instellation. Introducing water vapour into the atmosphere, we show the runaway greenhouse instellation limit is greatly reduced due to the presence of superadiabatic layers where convection is inhibited. This moves the inner edge of the habitable zone from $\approx$ 1 AU for a G-star to 1.6 AU (3.85 AU) for a Hycean world with a H$_2$-He inventory of 1 bar (10 bar). For an M-star, the inner edge is equivalently moved from 0.17 AU to 0.28 AU (0.54 AU). Our results suggest that most of the current Hycean world observational targets are not likely to sustain a liquid water ocean. We present an analytical framework for interpreting our results, finding that the maximum possible OLR scales approximately inversely with the dry mass inventory of the atmosphere. We discuss the possible limitations of our 1D modelling and recommend the use of 3D convection-resolving models to explore the robustness of superadiabatic layers., Comment: Submitted to ApJ, comments welcome

Atmospheric metal enrichment (i.e., elements heavier than helium, also called "metallicity") is a key diagnostic of the formation of giant planets. The giant planets of the solar system exhibit an inverse relationship between mass and both their bulk metallicities and atmospheric metallicities. Extrasolar giant planets also display an inverse relationship between mass and bulk metallicity. However, there is significant scatter in the relationship and it is not known how atmospheric metallicity correlates with either planet mass or bulk metallicity. Here we show that the Saturn-mass exoplanet HD 149026b has an atmospheric metallicity 59 - 276 times solar (at 1 $\sigma$), which is greater than Saturn's atmospheric metallicity of ~7.5 times solar at >4 $\sigma$ confidence. This result is based on modeling CO$_2$ and H$_2$O absorption features in the thermal emission spectrum of the planet measured by JWST. HD 149026b is the most metal-rich giant planet known, with an estimated bulk heavy element abundance of 66 $\pm$ 2% by mass. We find that the atmospheric metallicities of both HD 149026b and the solar system giant planets are more correlated with bulk metallicity than planet mass., Comment: Published online in Nature on March 27, 2023; a JWST thermal emission spectrum of a planet that definitely has an atmosphere ;)

The chemical inventory of hot Jupiter (HJ) exoplanets atmospheres continue to be observed by various ground and space based instruments in increasing detail and precision. It is expected that some HJs will exhibit strong non-equilibrium chemistry characteristics in their atmospheres, which might be inferred from spectral observations. We aim to model the three dimensional thermochemical non-equilibrium chemistry in the atmospheres of the HJs WASP-39b and HD 189733b. We couple a lightweight, reduced chemical network `mini-chem' that utilises net reaction rate tables to the Exo-FMS General Circulation Model (GCM). We perform GCM models of the exoplanets WASP-39b and HD 189733b as case studies of the coupled mini-chem scheme. The GCM results are then post-processed using the 3D radiative-transfer model gCMCRT to produce transmission and emission spectra to assess the impact of non-equilibrium chemistry on their observable properties. Both simulations show significant departures from chemical equilibrium (CE) due to the dynamical motions of the atmosphere. The spacial distribution of species generally follows closely the dynamical features of the atmosphere rather than the temperature field. Each molecular species exhibits a different quench level in the simulations, also dependent on the latitude of the planet. Major differences are seen in the transmission and emission spectral features between the CE and kinetic models. Our simulations indicate that considering the 3D kinetic chemical structures of HJ atmospheres has an important impact on physical interpretation of observational data. Drawing bulk atmospheric parameters from fitting feature strengths may lead to inaccurate interpretation of chemical conditions in the atmosphere of HJs. Our open source mini-chem module is simple to couple with contemporary HJ GCM models without substantially increasing required computational resources., Comment: A&A Accepted 17 Feb 2023

Photochemistry is a fundamental process of planetary atmospheres that regulates the atmospheric composition and stability. However, no unambiguous photochemical products have been detected in exoplanet atmospheres to date. Recent observations from the JWST Transiting Exoplanet Early Release Science Program found a spectral absorption feature at 4.05 $\mu$m arising from SO$_2$ in the atmosphere of WASP-39b. WASP-39b is a 1.27-Jupiter-radii, Saturn-mass (0.28 M$_J$) gas giant exoplanet orbiting a Sun-like star with an equilibrium temperature of $\sim$1100 K. The most plausible way of generating SO$_2$ in such an atmosphere is through photochemical processes. Here we show that the SO$_2$ distribution computed by a suite of photochemical models robustly explains the 4.05 $\mu$m spectral feature identified by JWST transmission observations with NIRSpec PRISM (2.7$\sigma$) and G395H (4.5$\sigma$). SO$_2$ is produced by successive oxidation of sulphur radicals freed when hydrogen sulphide (H$_2$S) is destroyed. The sensitivity of the SO$_2$ feature to the enrichment of the atmosphere by heavy elements (metallicity) suggests that it can be used as a tracer of atmospheric properties, with WASP-39b exhibiting an inferred metallicity of $\sim$10$\times$ solar. We further point out that SO$_2$ also shows observable features at ultraviolet and thermal infrared wavelengths not available from the existing observations., Comment: 39 pages, 14 figures, accepted to be published in Nature

Measuring the metallicity and carbon-to-oxygen (C/O) ratio in exoplanet atmospheres is a fundamental step towards constraining the dominant chemical processes at work and, if in equilibrium, revealing planet formation histories. Transmission spectroscopy provides the necessary means by constraining the abundances of oxygen- and carbon-bearing species; however, this requires broad wavelength coverage, moderate spectral resolution, and high precision that, together, are not achievable with previous observatories. Now that JWST has commenced science operations, we are able to observe exoplanets at previously uncharted wavelengths and spectral resolutions. Here we report time-series observations of the transiting exoplanet WASP-39b using JWST's Near InfraRed Camera (NIRCam). The long-wavelength spectroscopic and short-wavelength photometric light curves span 2.0 - 4.0 $\mu$m, exhibit minimal systematics, and reveal well-defined molecular absorption features in the planet's spectrum. Specifically, we detect gaseous H$_2$O in the atmosphere and place an upper limit on the abundance of CH$_4$. The otherwise prominent CO$_2$ feature at 2.8 $\mu$m is largely masked by H$_2$O. The best-fit chemical equilibrium models favour an atmospheric metallicity of 1-100$\times$ solar (i.e., an enrichment of elements heavier than helium relative to the Sun) and a sub-stellar carbon-to-oxygen (C/O) ratio. The inferred high metallicity and low C/O ratio may indicate significant accretion of solid materials during planet formation or disequilibrium processes in the upper atmosphere., Comment: 35 pages, 13 figures, 3 tables, Nature, accepted

Measuring the abundances of carbon and oxygen in exoplanet atmospheres is considered a crucial avenue for unlocking the formation and evolution of exoplanetary systems. Access to an exoplanet's chemical inventory requires high-precision observations, often inferred from individual molecular detections with low-resolution space-based and high-resolution ground-based facilities. Here we report the medium-resolution (R$\sim$600) transmission spectrum of an exoplanet atmosphere between 3-5 $\mu$m covering multiple absorption features for the Saturn-mass exoplanet WASP-39b, obtained with JWST NIRSpec G395H. Our observations achieve 1.46x photon precision, providing an average transit depth uncertainty of 221 ppm per spectroscopic bin, and present minimal impacts from systematic effects. We detect significant absorption from CO$_2$ (28.5$\sigma$) and H$_2$O (21.5$\sigma$), and identify SO$_2$ as the source of absorption at 4.1 $\mu$m (4.8$\sigma$). Best-fit atmospheric models range between 3 and 10x solar metallicity, with sub-solar to solar C/O ratios. These results, including the detection of SO$_2$, underscore the importance of characterising the chemistry in exoplanet atmospheres, and showcase NIRSpec G395H as an excellent mode for time series observations over this critical wavelength range., Comment: 44 pages, 11 figures, 3 tables. Resubmitted after revision to Nature

The atmospheres of ultra hot Jupiters (UHJs) are prime targets for the detection of molecules and atoms at both low and high spectral resolution. We study the atmospheres of the UHJs WASP-121b and WASP-189b by performing 3D general circulation models (GCMs) of these planets using high temperature correlated-k opacity schemes with ultra-violet (UV) absorbing species included. The GCM results are then post-processed at low and high spectral resolutions and compared to available data. The high resolution results are cross-correlated with molecular and atomic templates to produce mock molecular detections. Our GCM models produce similar temperature-pressure (T-p) structure trends to previous 1D radiative-convective equilibrium models of UHJs. Furthermore, the inclusion of UV opacities greatly shapes the thermal and dynamical properties of the high-altitude, low-pressure regions of the UHJ atmospheres, with sharp T-p inversions due to the absorption of UV light. This suggests that optical wavelength, high-resolution observations probe a dynamically distinct upper atmospheric region, rather than the deeper jet forming layers., Comment: MNRAS accepted: 8 August 2022, published 11 August 2022

The gas giant HD 80606 b has a highly eccentric orbit (e $\sim$ 0.93). The variation due to the rapid shift of stellar irradiation provides a unique opportunity to probe the physical and chemical timescales and to study the interplay between climate dynamics and atmospheric chemistry. In this work, we present integrated models to study the atmospheric responses and the underlying physical and chemical mechanisms of HD 80606 b. We first run three-dimensional general circulation models (GCMs) to establish the atmospheric thermal and dynamical structures for different atmospheric metallicities and internal heat. Based on the GCM output, we then adopted a 1D time-dependent photochemical model to investigate the compositional variation along the eccentric orbit. The transition of the circulation patterns of HD 80606 b matched the dynamics regimes in previous works. Our photochemical models show that efficient vertical mixing leads to deep quench levels of the major carbon and nitrogen species and the quenching behavior does not change throughout the eccentric orbit. Instead, photolysis is the main driver of the time-dependent chemistry. While CH$_4$ dominates over CO through most of the orbits, a transient state of [CO]/[CH$_4$}] $>$ 1 after periastron is confirmed for all metallicity and internal heat cases. The upcoming JWST Cycle 1 GO program will be able to track this real-time CH$_4$--CO conversion and infer the chemical timescale. Furthermore, sulfur species initiated by sudden heating and photochemical forcing exhibit both short-term and long-term cycles, opening an interesting avenue for detecting sulfur on exoplanets., Comment: 18 pages, 18 figures, published in MNRAS

Carbon dioxide (CO2) is a key chemical species that is found in a wide range of planetary atmospheres. In the context of exoplanets, CO2 is an indicator of the metal enrichment (i.e., elements heavier than helium, also called "metallicity"), and thus formation processes of the primary atmospheres of hot gas giants. It is also one of the most promising species to detect in the secondary atmospheres of terrestrial exoplanets. Previous photometric measurements of transiting planets with the Spitzer Space Telescope have given hints of the presence of CO2 but have not yielded definitive detections due to the lack of unambiguous spectroscopic identification. Here we present the detection of CO2 in the atmosphere of the gas giant exoplanet WASP-39b from transmission spectroscopy observations obtained with JWST as part of the Early Release Science Program (ERS). The data used in this study span 3.0 to 5.5 {\mu}m in wavelength and show a prominent CO2 absorption feature at 4.3 {\mu}m (26{\sigma} significance). The overall spectrum is well matched by one-dimensional, 10x solar metallicity models that assume radiative-convective-thermochemical equilibrium and have moderate cloud opacity. These models predict that the atmosphere should have water, carbon monoxide, and hydrogen sulfide in addition to CO2, but little methane. Furthermore, we also tentatively detect a small absorption feature near 4.0 {\mu}m that is not reproduced by these models., Comment: 27 pages, 6 figures, Accepted for publication in Nature, data and models available at https://doi.10.5281/zenodo.6959427

M-dwarfs are thought to be hostile environments for exoplanets. Stellar events are very common on such stars. These events might cause the atmospheres of exoplanets to change significantly over time. It is not only the major stellar flare events that contribute to this disequilibrium, but the smaller flares might also affect the atmospheres in an accumulating manner. In this study, we aim to investigate the effects of time-dependent stellar activity on the atmospheres of known exoplanets. We simulate the chemistry of GJ876c, GJ581c, and GJ832c that go from H$_2$-dominated to N$_2$-dominated atmospheres using observed stellar spectra from the MUSCLES-collaboration. We make use of the chemical kinetics code VULCAN and implement a flaring routine that stochastically generates synthetic flares based on observed flare statistics. Using the radiative transfer code petitRADTrans we also simulate the evolution of emission and transmission spectra. We investigate the effect of recurring flares for a total of 11 days covering 515 flares. Results show a significant change in abundance for some relevant species such as H, OH, and CH$_4$, with factors going up to 3 orders of magnitude difference with respect to the preflare abundances. We find a maximum change of $\sim$12 ppm for CH$_4$ in transmission spectra on GJ876c. These changes in the spectra remain too small to observe. We also find that the change in abundance and spectra of the planets accumulate throughout time, causing permanent changes in the chemistry. We conclude this small but gradual change in chemistry arises due to the recurring flares., Comment: Accepted for publication in MNRAS

The elemental ratios of carbon, nitrogen, and oxygen in the atmospheres of hot Jupiters may hold clues to their formation locations in the protostellar disc. In this work, we adopt gas phase chemical abundances of C, N and O from several locations in a disc chemical kinetics model as sources for the envelope of the hot Jupiter HD 209458b and evolve the planet's atmospheric composition using a 1D chemical kinetics model, treating both vertical mixing and photochemistry. We consider two atmospheric pressure-temperature profiles, one with and one without a thermal inversion. From each of the resulting 32 atmospheric composition profiles, we find that the molecules CH4, NH3, HCN, and C2H2 are more prominent in the atmospheres computed using a realistic non-inverted P-T profile in comparison to a prior equilibrium chemistry based work which used an analytical P-T profile. We also compute the synthetic transmission and emission spectra for these atmospheres and find that many spectral features vary with the location in the disc where the planet's envelope was accreted. By comparing with the species detected using the latest high-resolution ground-based observations, our model suggests HD 209458b could have accreted most of its gas between the CO2 and CH4 icelines with a super solar C/O ratio from its protostellar disc, which in turn directly inherited its chemical abundances from the protostellar cloud. Finally, we simulate observing the planet with the James Webb Space Telescope (JWST) and show that differences in spectral signatures of key species can be recognized. Our study demonstrates the enormous importance of JWST in providing new insights into hot Jupiter's formation environments., Comment: Accepted for publication in ApJ, 28 Pages, 11 figures, 2 Tables (Appendix: 5 Figures, 1 Table)

Growing evidence has indicated that the global composition distribution plays an indisputable role in interpreting observational data. 3D general circulation models (GCMs) with a reliable treatment of chemistry and clouds are particularly crucial in preparing for the upcoming observations. In the effort of achieving 3D chemistry-climate modeling, the challenge mainly lies in the expensive computing power required for treating a large number of chemical species and reactions. Motivated by the need for a robust and computationally efficient chemical scheme, we devise a mini-chemical network with a minimal number of species and reactions for H$_2$-dominated atmospheres. We apply a novel technique to simplify the chemical network from a full kinetics model -- VULCAN by replacing a large number of intermediate reactions with net reactions. The number of chemical species is cut down from 67 to 12, with the major species of thermal and observational importance retained, including H$_2$O, CH$_4$, CO, CO$_2$, C$_2$H$_2$, NH$_3$, and HCN. The size of the total reactions is greatly reduced from $\sim$ 800 to 20. The mini-chemical scheme is validated by verifying the temporal evolution and benchmarking the predicted compositions in four exoplanet atmospheres (GJ 1214b, GJ 436b, HD 189733b, HD 209458b) against the full kinetics of VULCAN. It reproduces the chemical timescales and composition distributions of the full kinetics well within an order of magnitude for the major species in the pressure range of 1 bar -- 0.1 mbar across various metallicities and carbon-to-oxygen (C/O) ratios. The small scale of the mini-chemical scheme permits simple use and fast computation, which is optimal for implementation in a 3D GCM or a retrieval framework. We focus on the thermochemical kinetics of net reactions in this paper and address photochemistry in a follow-up paper., Comment: 9 pages, 5 figures, accepted for publication in A&A

Planets smaller than Neptune and larger than Earth make up the majority of the discovered exoplanets. Those with H$_2$-rich atmospheres are prime targets for atmospheric characterization. The transition between the two main classes, super-Earths and sub-Neptunes, is not clearly understood as the rocky surface is likely not accessible to observations. Tracking several trace gases (specifically the loss of ammonia (NH$_3$) and hydrogen cyanide (HCN)) has been proposed as a proxy for the presence of a shallow surface. In this work, we revisit the proposed mechanism of nitrogen conversion in detail and find its timescale on the order of a million years. NH$_3$ exhibits dual paths converting to N$_2$ or HCN, depending on the UV radiation of the star and the stage of the system. In addition, methanol (CH$_3$OH) is identified as a robust and complementary proxy for a shallow surface. We follow the fiducial example of K2-18b with a 2D photochemical model (VULCAN) on an equatorial plane. We find a fairly uniform composition distribution below 0.1 mbar controlled by the dayside, as a result of slow chemical evolution. NH$_3$ and CH$_3$OH are concluded to be the most unambiguous proxies to infer surfaces on sub-Neptunes in the era of the James Webb Space Telescope (JWST)., Comment: Accepted for publication in ApJL

We present an update of the open-source photochemical kinetics code VULCAN (Tsai et al. 2017; https://github.com/exoclime/VULCAN) to include C-H-N-O-S networks and photochemistry. Additional new features are advection transport, condensation, various boundary conditions, and temperature-dependent UV cross-sections. First, we validate our photochemical model for hot Jupiter atmospheres by performing an intercomparison of HD 189733b models between Moses et al. (2011), Venot et al. (2012), and VULCAN, to diagnose possible sources of discrepancy. Second, we set up a model of Jupiter extending from the deep troposphere to upper stratosphere to verify the kinetics for low temperature. Our model reproduces hydrocarbons consistent with observations, and the condensation scheme successfully predicts the locations of water and ammonia ice clouds. We show that vertical advection can regulate the local ammonia distribution in the deep atmosphere. Third, we validate the model for oxidizing atmospheres by simulating Earth and find agreement with observations. Last, VULCAN is applied to four representative cases of extrasolar giant planets: WASP-33b, HD 189733b, GJ 436b, and 51 Eridani b. We look into the effects of the C/O ratio and chemistry of titanium/vanadium species for WASP-33b; we revisit HD 189733b for the effects of sulfur and carbon condensation; the effects of internal heating and vertical mixing ($K_{\textrm{zz}}$) are explored for GJ 436b; we test updated planetary properties for 51 Eridani b with S$_8$ condensates. We find sulfur can couple to carbon or nitrogen and impact other species such as hydrogen, methane, and ammonia. The observable features of the synthetic spectra and trends in the photochemical haze precursors are discussed for each case., Comment: 50 pages, 41 figures, accepted for publication in ApJ

Radiative-transfer (RT) is a fundamental part of modelling exoplanet atmospheres with general circulation models (GCMs). An accurate RT scheme is required for estimates of the atmospheric energy transport and for gaining physical insight from model spectra. We implement three RT schemes for Exo-FMS: semi-grey, non-grey `picket fence', and real gas with correlated-k. We benchmark the Exo-FMS GCM using these RT schemes to hot Jupiter simulation results from the literature. We perform a HD 209458b-like simulation with the three schemes and compare their results. These simulations are then post-processed to compare their observable differences. The semi-grey scheme results show qualitative agreement with previous studies in line with variations seen between GCM models. The real gas model reproduces well the temperature and dynamical structures from other studies. After post-processing our non-grey picket fence scheme compares very favourably with the real gas model, producing similar transmission spectra, emission spectra and phase curve behaviours. Exo-FMS is able to reliably reproduce the essential features of contemporary GCM models in the hot gas giant regime. Our results suggest the picket fence approach offers a simple way to improve upon RT realism beyond semi-grey schemes., Comment: MNRAS accepted 22 June 2021 - V2, typos fixed

Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths. It is the first mission dedicated to measuring the chemical composition and thermal structures of hundreds of transiting exoplanets, enabling planetary science far beyond the boundaries of the Solar System. The payload consists of an off-axis Cassegrain telescope (primary mirror 1100 mm x 730 mm ellipse) and two separate instruments (FGS and AIRS) covering simultaneously 0.5-7.8 micron spectral range. The satellite is best placed into an L2 orbit to maximise the thermal stability and the field of regard. The payload module is passively cooled via a series of V-Groove radiators; the detectors for the AIRS are the only items that require active cooling via an active Ne JT cooler. The Ariel payload is developed by a consortium of more than 50 institutes from 16 ESA countries, which include the UK, France, Italy, Belgium, Poland, Spain, Austria, Denmark, Ireland, Portugal, Czech Republic, Hungary, the Netherlands, Sweden, Norway, Estonia, and a NASA contribution., Comment: Ariel Definition Study Report, 147 pages. Reviewed by ESA Science Advisory Structure in November 2020. Original document available at: https://www.cosmos.esa.int/documents/1783156/3267291/Ariel_RedBook_Nov2020.pdf/

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 H$_2$, H$_2$O, CO$_2$, CH$_4$, CO, O$_2$, or N$_2$. 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, N$_2$, and O$_2$ with minimal effect, H$_2$O, CO$_2$, and CH$_4$ with intermediate influence, and H$_2$ 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., Comment: Accepted for publication in JGR Planets; 23 pages, 10 figures; summaries available at https://bit.ly/MO21_blog (blog) and http://bit.ly/MO21_video (video)

The atmospheric circulation of tidally locked planets is dominated by a superrotating eastward equatorial jet. We develop a predictive theory for the formation of this jet, proposing a mechanism in which the three-dimensional stationary waves induced by the day-night forcing gradient produce an equatorial acceleration. This is balanced in equilibrium by an interaction between the resulting jet and the vertical motion of the atmosphere. The three-dimensional structure of the zonal acceleration is vital to this mechanism. We demonstrate this mechanism in a hierarchy of models. We calculate the three-dimensional stationary waves induced by the forcing on these planets, and show the vertical structure of the zonal acceleration produced by these waves, which we use to suggest a mechanism for how the jet forms. GCM simulations are used to confirm the equilibrium state predicted by this mechanism, where the acceleration from these waves is balanced by an interaction between the zonal-mean vertical velocity and the jet. We derive a simple model of this using the "Weak Temperature Gradient" approximation, which gives an estimate of the jet speed on a terrestrial tidally locked planet. We conclude that the proposed mechanism is a good description of the formation of an equatorial jet on a terrestrial tidally locked planet, and should be useful for interpreting observations and simulations of these planets. The mechanism requires assumptions such as a large equatorial Rossby radius and weak acceleration due to transient waves, and a different mechanism may produce the equatorial jets on gaseous tidally locked planets., Comment: Accepted by ApJ

We explore the chemistry and observability of nitrogen dominated atmospheres for ultra-short-period super-Earths. We base the assumption, that super-Earths could have nitrogen filled atmospheres, on observations of 55 Cnc e that favour a scenario with a high-mean-molecular-weight atmosphere. We take Titan's elemental budget as our starting point and using chemical kinetics compute a large range of possible compositions for a hot super-Earth. We use analytical temperature profiles and explore a parameter space spanning orders of magnitude in C/O \& N/O ratios, while always keeping nitrogen the dominant component. We generate synthetic transmission and emission spectra and assess their potential observability with the future James Webb Space Telescope and ARIEL. Our results suggest that HCN is a strong indicator of a high C/O ratio, which is similar to what is found for H-dominated atmospheres. We find that these worlds are likely to possess C/O > 1.0, and that HCN, CN, CO should be the primary molecules to be searched for in thermal emission. For lower temperatures (T < 1500 K), we additionally find NH3 in high N/O ratio cases, and C2H4, CH4 in low N/O ratio cases to be strong absorbers. Depletion of hydrogen in such atmospheres would make CN, CO and NO exceptionally prominent molecules to look for in the 0.6 - 5.0 $\mu$m range. Our models show that the upcoming JWST and ARIEL missions will be able to distinguish atmospheric compositions of ultra-short period super-Earths with unprecedented confidence., Comment: Accepted for publication MNRAS

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., Comment: Accepted in A&A

Context: White dwarf - Brown dwarf short period binaries (P$_{\rm orb}$ $\lesssim$ 2 hours) 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. Aims: We aim to investigate the three dimensional atmospheric structural and dynamical properties of the Brown dwarf WD0137-349B. Methods: We use the three dimensional GCM model Exo-FMS, with a dual-band grey radiative-transfer scheme to model the atmosphere of WD0137-349B. The results of the GCM model are post-processed using the three dimensional Monte Carlo radiative-transfer model \textsc{cmcrt}. Results: Our results suggest inefficient day-night energy transport and a large day-night temperature contrast for WD0137-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-IR emission excess at $\gtrsim$ 1.95 $\mu$m as observed by the GNIRS instrument. Our model over predicts the IR phase curve fluxes by factors of $\approx$1-3, but generally fits the shape of the phase curves well. Chemical kinetic modelling using \textsc{vulcan} suggests a highly ionised region at high altitudes can form on the dayside of the Brown dwarf. Conclusions: 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 UV irradiation is 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., Comment: Version 2 - Accepted June 19 2020 (MNRAS)

THOR is the first open-source general circulation model (GCM) developed from scratch to study the atmospheres and climates of exoplanets, free from Earth- or Solar System-centric tunings. It solves the general non-hydrostatic Euler equations (instead of the primitive equations) on a sphere using the icosahedral grid. In the current study, we report major upgrades to THOR, building upon the work of Mendon\c{c}a et al. (2016). First, while the Horizontally Explicit Vertically Implicit (HEVI) integration scheme is the same as that described in Mendon\c{c}a et al. (2016), we provide a clearer description of the scheme and improved its implementation in the code. The differences in implementation between the hydrostatic shallow (HSS), quasi-hydrostatic deep (QHD) and non-hydrostatic deep (NHD) treatments are fully detailed. Second, standard physics modules are added: two-stream, double-gray radiative transfer and dry convective adjustment. Third, THOR is tested on additional benchmarks: tidally-locked Earth, deep hot Jupiter, acoustic wave, and gravity wave. Fourth, we report that differences between the hydrostatic and non-hydrostatic simulations are negligible in the Earth case, but pronounced in the hot Jupiter case. Finally, the effects of the so-called "sponge layer", a form of drag implemented in most GCMs to provide numerical stability, are examined. Overall, these upgrades have improved the flexibility, user-friendliness, and stability of THOR., Comment: 57 pages, 31 figures, revised, accepted for publication in ApJS

We present new methodological features and physical ingredients included in the 1D radiative transfer code HELIOS, improving the hemispheric two-stream formalism. We conduct a thorough intercomparison survey with several established forward models, including COOLTLUSTY, PHOENIX, and find satisfactory consistency with their results. Then, we explore the impact of (i) different groups of opacity sources, (ii) a stellar path length adjustment, and (iii) a scattering correction on self-consistently calculated atmospheric temperatures and planetary emission spectra. First, we observe that temperature-pressure (T-P) profiles are very sensitive to the opacities included, with metal oxides, hydrides, the alkali atoms (and ionized hydrogen) playing an important role for the absorption of shortwave radiation (in very hot surroundings). Moreover, if these species are sufficiently abundant, they are likely to induce non-monotonic T-P profiles. Second, without the stellar path length adjustment, the incoming stellar flux is significantly underestimated for zenith angles above 80{\deg}, which somewhat affects the upper atmospheric temperatures and the planetary emission. Third, the scattering correction improves the accuracy of the computation of the reflected stellar light by ~10%. We use HELIOS to calculate a grid of cloud-free atmospheres in radiative-convective equilibrium for self-luminous planets for a range of effective temperatures, surface gravities, metallicities, and C/O ratios, to be used by planetary evolution studies. Furthermore, we calculate dayside temperatures and secondary eclipse spectra for a sample of exoplanets for varying chemistry and heat redistribution. These results may be used to make predictions on the feasibility of atmospheric characterizations with future observations., Comment: 33 pages, 18 figures, accepted for publication in AJ

Spectral features in the observed spectra of exoplanets depend on the composition of their atmospheres. A good knowledge of the main atmospheric processes that drive the chemical distribution is therefore essential to interpret exoplanetary spectra. An atmosphere reaches chemical equilibrium if the rates of the forward and backward chemical reactions converge to the same value. However, there are atmospheric processes, such as atmospheric transport, that destabilize this equilibrium. In this work we study the changes in composition driven by a 3D wind field in WASP-43b using our Global Circulation Model, THOR. Our model uses validated temperature- and pressure-dependent chemical timescales that allow us to explore the disequilibrium chemistry of CO, CO$_2$, H$_2$O and CH$_4$. In WASP-43b the formation of the equatorial jet has an important impact in the chemical distribution of the different species across the atmosphere. At low latitudes the chemistry is longitudinally quenched, except for CO$_2$ at solar abundances. The polar vortexes have a distinct chemical distribution since these are regions with lower temperature and atmospheric mixing. Vertical and latitudinal mixing have a secondary impact in the chemical transport. We determine graphically the effect of disequilibrium on observed emission spectra. Our results do not show any significant differences in the emission spectra between the equilibrium and disequilibrium solutions for C/O = 0.5. However, if C/O is increased to 2.0, differences in the spectra due to the disequilibrium chemistry of CH$_4$ become non-negligible. In some spectral ranges the emission spectra can have more than 15$\%$ departures from the equilibrium solution., Comment: Accepted for publication in ApJ

The atmospheric temperatures of the ultra-hot Jupiter KELT-9b straddle the transition between gas giants and stars, and therefore between two traditionally distinct regimes of atmospheric chemistry. Previous theoretical studies assume the atmosphere of KELT-9b to be in chemical equilibrium. Despite the high ultraviolet flux from KELT-9, we show using photochemical kinetics calculations that the observable atmosphere of KELT-9b is predicted to be close to chemical equilibrium, which greatly simplifies any theoretical interpretation of its spectra. It also makes the atmosphere of KELT-9b, which is expected to be cloudfree, a tightly constrained chemical system that lends itself to a clean set of theoretical predictions. Due to the lower pressures probed in transmission (compared to emission) spectroscopy, we predict the abundance of water to vary by several orders of magnitude across the atmospheric limb depending on temperature, which makes water a sensitive thermometer. Carbon monoxide is predicted to be the dominant molecule under a wide range of scenarios, rendering it a robust diagnostic of the metallicity when analyzed in tandem with water. All of the other usual suspects (acetylene, ammonia, carbon dioxide, hydrogen cyanide, methane) are predicted to be subdominant at solar metallicity, while atomic oxygen, iron and magnesium are predicted to have relative abundances as high as 1 part in 10,000. Neutral atomic iron is predicted to be seen through a forest of optical and near-infrared lines, which makes KELT-9b suitable for high-resolution ground-based spectroscopy with HARPS-N or CARMENES. We summarize future observational prospects of characterizing the atmosphere of KELT-9b., Comment: Accepted by ApJ. 9 pages, 6 figures. Corrected minor errors in Figures 1a and 1b (some line styles were switched by accident), text and conclusions unchanged, these minor changes will be updated in final ApJ proof