Abstract The calculation of internal atmospheric (longwave) fluxes is a key component of any model of exoplanet atmospheres that requires radiative-transfer (RT) calculations. For atmospheres ...containing a strong scattering component such as cloud particles, most 1D multiple-scattering RT methods typically involve numerically expensive matrix inversions. This computational bottleneck is exacerbated when multitudes of RT calculations are required, such as in general circulation models (GCMs) and retrieval methods. In an effort to increase the speed of RT calculations without sacrificing too much accuracy, we investigate the applicability of approximate longwave scattering methods developed for the Earth science community to hot Jupiter atmospheres. We test the absorption approximation and variational iteration method (VIM) applied to typical cloudy hot Jupiter scenarios, using 64-stream DISORT calculations as reference solutions. We find the four-stream VIM variant is a highly promising method to explore for use in hot Jupiter GCM and retrieval modeling, and it shows excellent speed characteristics, with typical errors ∼1% for outgoing fluxes and within ∼50%, but with larger errors in the test case of a deep cloud layer, for vertical heating rates. Other methods explored in this study were found to typically produce similar error characteristics in vertical heating rates.
ABSTRACT
Transmission spectra contain a wealth of information about the atmospheres of transiting exoplanets. However, large thermal and chemical gradients along the line of sight can lead to biased ...inferences in atmospheric retrievals. In order to determine how far from the limb plane the atmosphere still impacts the transmission spectrum, we derive a new formula to estimate the opening angle of a planet. This is the angle subtended by the atmospheric region that contributes to the observation along the line of sight, as seen from the planet centre. We benchmark our formula with a 3D Monte Carlo radiative transfer code and we define an opening angle suitable for the interpretation of JWST observations, assuming a 10-ppm noise floor. We find that the opening angle is only a few degrees for planets cooler than ca. 500 Kelvins, while it can be as large as 25 degrees for (ultra-)hot Jupiters and 50 degrees for hot Neptunes. Compared to previous works, our more robust approach leads to smaller estimates for the opening angle across a wide range scale heights and planetary radii. Finally, we show that ultra-hot Jupiters have an opening angle that is smaller than the angle over which the planet rotates during the transit. This allows for time-resolved transmission spectroscopy observations that probe independent parts of the planetary limb during the first and second half of the transit.
ABSTRACT
Ultra-hot Jupiters are tidally locked gas giants with dayside temperatures high enough to dissociate hydrogen and other molecules. Their atmospheres are vastly non-uniform in terms of ...chemistry, temperature, and dynamics, and this makes their high-resolution transmission spectra and cross-correlation signal difficult to interpret. In this work, we use the SPARC/MITgcm global circulation model to simulate the atmosphere of the ultra-hot Jupiter WASP-76b under different conditions, such as atmospheric drag and the absence of TiO and VO. We then employ a 3D Monte Carlo radiative transfer code, hires-mcrt, to self-consistently model high-resolution transmission spectra with iron (Fe i) lines at different phases during the transit. To untangle the structure of the resulting cross-correlation map, we decompose the limb of the planet into four sectors, and we analyse each of their contributions separately. Our experiments demonstrate that the cross-correlation signal of an ultra-hot Jupiter is primarily driven by its temperature structure, rotation, and dynamics, while being less sensitive to the precise distribution of iron across the atmosphere. We also show that the previously published iron signal of WASP-76b can be reproduced by a model featuring iron condensation on the leading limb. Alternatively, the signal may be explained by a substantial temperature asymmetry between the trailing and leading limb, where iron condensation is not strictly required to match the data. Finally, we compute the Kp–Vsys maps of the simulated WASP-76b atmospheres, and we show that rotation and dynamics can lead to multiple peaks that are displaced from zero in the planetary rest frame.
ABSTRACT
Insight into the formation and global distribution of cloud particles in exoplanet atmospheres continues to be a key problem to tackle going into the JWST era. Understanding microphysical ...cloud processes and atmospheric feedback mechanisms in three-dimensional (3D) has proven to be a challenging prospect for exoplaneteers. In an effort to address the large computational burden of coupling these models in 3D simulations, we develop an open source, lightweight, and efficient microphysical cloud model for exoplanet atmospheres. ‘Mini-cloud’ is a microphysical based cloud model for exoplanet condensate clouds that can be coupled to contemporary general circulation models (GCMs) and other time-dependent simulations. We couple mini-cloud to the Exo-FMS GCM and use a prime JWST target, the hot Jupiter HAT-P-1b, as a test case for the cloud formation module. After 1000+ of days of integration with mini-cloud, our results show a complex 3D cloud structure with cloud properties relating closely the dynamical and temperature properties of the atmosphere. Current transit and emission spectra data are best fit with a reduced cloud particle number density compared to the nominal simulation, with our simulated JWST NIRISS SOSS spectra showing promising prospects for characterizing the atmosphere in detail. Overall, our study is another small step in first principles 3D exoplanet cloud formation microphysical modelling. We suggest that additional physics not included in the present model, such as coagulation, are required to reduce the number density of particles to appropriately observed levels.
Abstract
Inhomogeneous cloud formation and wavelength-dependent phenomena are expected to shape hot Jupiter atmospheres. We present a general circulation model with multiwavelength “picket fence” ...radiative transfer and radiatively active, temperature-dependent clouds, and compare the results to those of a double gray routine. The double gray method inherently fails to model polychromatic effects in hot Jupiter atmospheres, while picket fence captures these non-gray aspects and performs well compared to fully wavelength-dependent methods. We compare both methods with radiatively active clouds and cloud-free models, assessing the limitations of the double gray method. Although there are broad similarities, the picket fence models have larger dayside–nightside temperature differences, nonisothermal upper atmospheres, and multiwavelength effects in the presence of radiatively active clouds. We model the well-known hot Jupiters HD 189733 b and HD 209458 b. For the hotter HD 209458 b, the picket fence method prevents clouds from thermostating dayside temperatures, resulting in hotter upper atmospheres and the dissipation of dayside clouds. Differences in the temperature structures are then associated with nuanced differences in the circulation patterns and clouds. Models of the cooler HD 189733 b have global cloud coverage, regardless of the radiative transfer scheme, whereas there are larger differences in the models of HD 209458 b, particularly in the extent of the partial cloud coverage on its dayside. This results in minor changes to the thermal and reflected light phase curves of HD 189733 b, but more significant differences for the picket fence and double gray versions of HD 209458 b.
We present the transmission spectrum of HAT-P-12b through a joint analysis of data obtained from the Hubble Space Telescope Space Telescope Imaging Spectrograph and Wide Field Camera 3 and Spitzer, ...covering the wavelength range 0.3-5.0 m. We detect a muted water vapor absorption feature at 1.4 m attenuated by clouds, as well as a Rayleigh scattering slope in the optical indicative of small particles. We interpret the transmission spectrum using both the state-of-the-art atmospheric retrieval code SCARLET and the aerosol microphysics model CARMA. These models indicate that the atmosphere of HAT-P-12b is consistent with a broad range of metallicities between several tens to a few hundred times solar, a roughly solar C/O ratio, and moderately efficient vertical mixing. Cloud models that include condensate clouds do not readily generate the submicron particles necessary to reproduce the observed Rayleigh scattering slope, while models that incorporate photochemical hazes composed of soot or tholins are able to match the full transmission spectrum. From a complementary analysis of secondary eclipses by Spitzer, we obtain measured depths of 0.042% 0.013% and 0.045% 0.018% at 3.6 and 4.5 m, respectively, which are consistent with a blackbody temperature of 890+60−70 K and indicate efficient day-night heat recirculation. HAT-P-12b joins the growing number of well-characterized warm planets that underscore the importance of clouds and hazes in our understanding of exoplanet atmospheres.
ABSTRACT
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: semigrey, 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 semigrey 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 semigrey schemes.
ABSTRACT
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 idealized 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 and global spatial variations due to zonal flows. Chemical species are generally globally homogenized, showing variations of ±10 per cent or less, dependent on pressure level, and follow the dynamical structures present in the atmosphere. However, we find localized storm regions and eddies can show higher contrasts, up to ±100 per cent, in mixing ratio compared to the background global mean. This initial study represents another step in understanding the connection between 3D atmospheric flows in brown dwarfs and their rich chemical inventories.
ABSTRACT
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 three-dimensional (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 3D Exo-FMS general circulation model 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 on to 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.