Abstract
We present an update of the open-source photochemical kinetics code VULCAN to include C–H–N–O–S networks and photochemistry. The 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., Venot et al., 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
zz
) are explored for GJ 436b, and we test updated planetary properties for 51 Eridani b with S
8
condensates. We find that 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.
ABSTRACT 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 (C2H2), ammonia (NH3), carbon dioxide (CO2), carbon monoxide (CO), ethylene (C2H4), hydrogen cyanide (HCN), methane (CH4), molecular nitrogen (N2), and water (H2O). By considering only the gas phase, we prove that the mixing ratio of carbon monoxide is governed by a decic equation (polynomial equation of 10 degrees). We validate our solutions against numerical calculations of equilibrium chemistry that perform Gibbs free energy minimization and demonstrate that they are accurate at the level for temperatures from 500 to 3000 K. In hydrogen-dominated atmospheres, the ratio of abundances of HCN to CH4 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.
Abstract
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 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.
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; exoclime.net) and publicly available at https://github.com/exoclime/VULCAN.
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, CO2, H2O, and CH4. In WASP-43b the formation of the equatorial jet has an important impact on the chemical distribution of the different species across the atmosphere. At low latitudes the chemistry is longitudinally quenched, except for CO2 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 on the chemical transport. We determine graphically the effect of disequilibrium on the 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 CH4 become non-negligible. In some spectral ranges the emission spectra can have more than 15% departure from the equilibrium solution.
The Peculiar Atmospheric Chemistry of KELT-9b Kitzmann, Daniel; Heng, Kevin; Rimmer, Paul B. ...
Astrophysical journal/The Astrophysical journal,
08/2018, Volume:
863, Issue:
2
Journal Article
Peer reviewed
Open access
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 kinetic 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 cloud-free, 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.
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. General circulation model 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.
We present new methodological features and physical ingredients included in the one-dimensional radiative transfer code HELIOS, improving the hemispheric two-stream formalism. We conduct a thorough ...intercomparison survey with several established forward models, including COOLTLUSTY and 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, and alkali atoms (and ionized hydrogen) playing an important role in the absorption of shortwave radiation (in very hot surroundings). Moreover, if these species are sufficiently abundant, they are likely to induce nonmonotonic T-P profiles. Second, without the stellar path length adjustment, the incoming stellar flux is significantly underestimated for zenith angles above 80°, 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.
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 toward 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 the 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 ∼10%, respectively. The level of accuracy attained allows for the chemical relaxation method to be included in three-dimensional general circulation models.
ABSTRACT
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 Cancri 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 and 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 (JWST) 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, and 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 and 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 $\rm{\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.