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.

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.

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.

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.

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.

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.

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.

Motivated by the work of Cooper & Showman, we revisit the chemical relaxation method, which seeks to enhance the computational efficiency of chemical-kinetics calculations by replacing the chemical network with a handful of independent source/sink terms. Chemical relaxation solves the evolution of the system and can treat disequilibrium chemistry, as the source/sink terms are driven towards chemical equilibrium on a prescribed chemical timescale, but it has surprisingly never been validated. First, we generalize the treatment by forgoing the use of a single chemical timescale, instead developing a pathway analysis tool that allows us to identify the rate-limiting reaction as a function of temperature and pressure. For the interconversion between methane and carbon monoxide and between ammonia, and molecular nitrogen, we identify the key rate-limiting reactions for conditions relevant to currently characterizable exo-atmospheres (500-3000 K, 0.1 mbar to 1 kbar). Second, we extend chemical relaxation to include carbon dioxide and water. Third, we examine the role of metallicity and carbon-to-oxygen ratio in chemical relaxation. Fourth, we apply our pathway analysis tool to diagnose the differences between our chemical network and that of Moses and Venot. Finally, we validate the chemical relaxation method against full chemical kinetics calculations in one dimension. For WASP-18b-, HD 189733b- and GJ 1214-b-like atmospheres, we show that chemical relaxation is mostly accurate to within an order of magnitude, a factor of 2 and $\sim 10\%$, respectively. The level of accuracy attained allows for the chemical relaxation method to be included in three-dimensional general circulation models.

We analyze the emission spectrum of the hot Jupiter WASP-12b using our HELIOS-R retrieval code and HELIOS-K opacity calculator. When interpreting Hubble and Spitzer data, the retrieval outcomes are found to be prior-dominated. When the prior distributions of the molecular abundances are assumed to be log-uniform, the volume mixing ratio of HCN is found to be implausibly high. A VULCAN chemical kinetics model of WASP-12b suggests that chemical equilibrium is a reasonable assumption even when atmospheric mixing is implausibly rigorous. Guided by (exo)planet formation theory, we set Gaussian priors on the elemental abundances of carbon, oxygen and nitrogen with the Gaussian peaks being centered on the measured C/H, O/H and N/H values of the star. By enforcing chemical equilibrium, we find substellar O/H and stellar to slightly superstellar C/H for the dayside atmosphere of WASP-12b. The superstellar carbon-to-oxygen ratio is just above unity, regardless of whether clouds are included in the retrieval analysis, consistent with Madhusudhan et al. (2011). Furthermore, whether a temperature inversion exists in the atmosphere depends on one's assumption for the Gaussian width of the priors. Our retrieved posterior distributions are consistent with the formation of WASP-12b in a solar-composition protoplanetary disk, beyond the water iceline, via gravitational instability or pebble accretion (without core erosion) and migration inwards to its present orbital location via a disk-free mechanism, and are inconsistent with both in-situ formation and core accretion with disk migration, as predicted by Madhusudhan et al. (2017). We predict that the interpretation of James Webb Space Telescope WASP-12b data will not be prior-dominated.

We present an open-source and validated chemical kinetics code for studying hot exoplanetary atmospheres, which we name VULCAN. It is constructed for gaseous chemistry from 500 to 2500 K using a reduced C-H-O chemical network with about 300 reactions. It uses eddy diffusion to mimic atmospheric dynamics and excludes photochemistry. We have provided a full description of the rate coefficients and thermodynamic data used. We validate VULCAN by reproducing chemical equilibrium and by comparing its output versus the disequilibrium-chemistry calculations of Moses et al. and Rimmer & Helling. It reproduces the models of HD 189733b and HD 209458b by Moses et al., which employ a network with nearly 1600 reactions. We also use VULCAN to examine the theoretical trends produced when the temperature-pressure profile and carbon-to-oxygen ratio are varied. Assisted by a sensitivity test designed to identify the key reactions responsible for producing a specific molecule, we revisit the quenching approximation and find that it is accurate for methane but breaks down for acetylene, because the disequilibrium abundance of acetylene is not directly determined by transport-induced quenching, but is rather indirectly controlled by the disequilibrium abundance of methane. Therefore, we suggest that the quenching approximation should be used with caution and must always be checked against a chemical kinetics calculation. A one-dimensional model atmosphere with 100 layers, computed using VULCAN, typically takes several minutes to complete. VULCAN is part of the Exoclimes Simulation Platform (ESP; http://www.exoclime.net/) and publicly available at http://github.com/exoclime/VULCAN .

We present the open-source radiative transfer code named HELIOS, which is constructed for studying exoplanetary atmospheres. In its initial version, the model atmospheres of HELIOS are one-dimensional and plane-parallel, and the equation of radiative transfer is solved in the two-stream approximation with non-isotropic scattering. A small set of the main infrared absorbers is employed, computed with the opacity calculator HELIOS-K and combined using a correlated-$k$ approximation. The molecular abundances originate from validated analytical formulae for equilibrium chemistry. We compare HELIOS with the work of Miller-Ricci & Fortney using a model of GJ 1214b, and perform several tests, where we find: model atmospheres with single-temperature layers struggle to converge to radiative equilibrium; $k$-distribution tables constructed with $\gtrsim 0.01$ cm$^{-1}$ resolution in the opacity function ($ \lesssim 10^3$ points per wavenumber bin) may result in errors $\gtrsim 1$-10 % in the synthetic spectra; and a diffusivity factor of 2 approximates well the exact radiative transfer solution in the limit of pure absorption. We construct "null-hypothesis" models (chemical equilibrium, radiative equilibrium and solar element abundances) for 6 hot Jupiters. We find that the dayside emission spectra of HD 189733b and WASP-43b are consistent with the null hypothesis, while it consistently under-predicts the observed fluxes of WASP-8b, WASP-12b, WASP-14b and WASP-33b. We demonstrate that our results are somewhat insensitive to the choice of stellar models (blackbody, Kurucz or PHOENIX) and metallicity, but are strongly affected by higher carbon-to-oxygen ratios. The code is publicly available as part of the Exoclimes Simulation Platform (ESP; exoclime.net).

We present novel, analytical, equilibrium-chemistry formulae for the abundances of molecules in hot exoplanetary atmospheres that include the carbon, oxygen and nitrogen networks. Our hydrogen-dominated solutions involve acetylene (C$_2$H$_2$), ammonia (NH$_3$), carbon dioxide (CO$_2$), carbon monoxide (CO), ethylene (C$_2$H$_4$), hydrogen cyanide (HCN), methane (CH$_4$), molecular nitrogen (N$_2$) and water (H$_2$O). By considering only the gas phase, we prove that the mixing ratio of carbon monoxide is governed by a decic equation (polynomial equation of degree 10). We validate our solutions against numerical calculations of equilibrium chemistry that perform Gibbs free energy minimization and demonstrate that they are accurate at the $\sim 1\%$ level for temperatures from 500 to 3000 K. In hydrogen-dominated atmospheres, the ratio of abundances of HCN to CH$_4$ is nearly constant across a wide range of carbon-to-oxygen ratios, which makes it a robust diagnostic of the metallicity in the gas phase. Our validated formulae allow for the convenient benchmarking of chemical kinetics codes and provide an efficient way of enforcing chemical equilibrium in atmospheric retrieval calculations.

We present a self-consistent formalism for computing and understanding the atmospheric chemistry of exoplanets from the viewpoint of an astrophysicist. Starting from the first law of thermodynamics, we demonstrate that the van't Hoff equation (which describes the equilibrium constant), Arrhenius equation (which describes the rate coefficients) and procedures associated with the Gibbs free energy (minimisation, rescaling) have a common physical and mathematical origin. We address an ambiguity associated with the equilibrium constant, which is used to relate the forward and reverse rate coefficients, and restate its two definitions. By necessity, one of the equilibrium constants must be dimensionless and equate to an exponential function involving the Gibbs free energy, while the other is a ratio of rate coefficients and must therefore possess physical units. We demonstrate that the Arrhenius equation takes on a functional form that is more general than previously stated without recourse to tagging on ad hoc functional forms. Finally, we derive analytical models of chemical systems, in equilibrium, with carbon, hydrogen and oxygen. We include acetylene and are able to reproduce several key trends, versus temperature and carbon-to-oxygen ratio, published in the literature. The rich variety of behavior that mixing ratios exhibit as a function of the carbon-to-oxygen ratio is merely the outcome of stoichiometric book-keeping and not the direct consequence of temperature or pressure variations.

Three-dimensional equatorial trapped waves excited by stellar isolation and the resulting equatorial superrotating jet in a vertical stratified atmosphere of a tidally-locked hot Jupiter are investigated. Taking the hot Jupiter HD 189733b as a fiducial example, we analytically solve a set of linear equations subject to stationary stellar heating with a uniform zonal-mean flow included. We also extract wave information in the final equilibrium state of the atmosphere from the radiative hydrodynamical simulation for HD 189733b by Dobbs-Dixon & Agol (2013). We find that the analytic wave solutions are able to qualitatively explain the three-dimensional simulation results. Studying the vertical structure of waves allows us to explore new wave features such as the westward tilt of wavefronts related to the Rossby-wave resonance as well as double gyres of dispersive Rossby waves. We also make an attempt to apply our linear wave analysis to explain some numerical features associated with the equatorial jet development seen in the GCM by Showman & Polvani (2011). During the spin-up phase of the equatorial jet, the acceleration of the jet as a result of the divergence of the wave momentum flux can be in principle boosted by the Rossby-wave resonance as the superrotating jet speed is close to the phase speed of the Rossby waves. However, we also find that as the jet speed increases, the Rossby-wave structure shifts eastward, while the Kelvin-wave structure remains approximately stationary, leading to the decline of the acceleration rate. Our analytic model of jet evolution implies that there exists only one stable equilibrium state of the atmosphere, possibly implying that the final state of the atmosphere is independent of initial conditions in the linear regime. Limitations and future improvements are also discussed.