ABSTRACT We present an open-source Thermochemical Equilibrium Abundances (TEA) code that calculates the abundances of gaseous molecular species. The code is based on the methodology of White et al. ...and Eriksson. It applies Gibbs free-energy minimization using an iterative, Lagrangian optimization scheme. Given elemental abundances, TEA calculates molecular abundances for a particular temperature and pressure or a list of temperature-pressure pairs. We tested the code against the method of Burrows & Sharp, the free thermochemical equilibrium code Chemical Equilibrium with Applications (CEA), and the example given by Burrows & Sharp. Using their thermodynamic data, TEA reproduces their final abundances, but with higher precision. We also applied the TEA abundance calculations to models of several hot-Jupiter exoplanets, producing expected results. TEA is written in Python in a modular format. There is a start guide, a user manual, and a code document in addition to this theory paper. TEA is available under a reproducible-research, open-source license via https://github.com/dzesmin/TEA.
The James Webb Space Telescope (JWST) is expected to revolutionize the field of exoplanets. The broad wavelength coverage and the high sensitivity of its instruments will allow characterization of ...exoplanetary atmospheres with unprecedented precision. Following the Call for the Cycle 1 Early Release Science Program, the Transiting Exoplanet Community was awarded time to observe several targets, including WASP-43b. The atmosphere of this hot Jupiter has been intensively observed but still harbors some mysteries, especially concerning the day-night temperature gradient, the efficiency of the atmospheric circulation, and the presence of nightside clouds. We will constrain these properties by observing a full orbit of the planet and extracting its spectroscopic phase curve in the 5-12 m range with JWST/MIRI. To prepare for these observations, we performed extensive modeling work with various codes: radiative transfer, chemical kinetics, cloud microphysics, global circulation models, JWST simulators, and spectral retrieval. Our JWST simulations show that we should achieve a precision of 210 ppm per 0.1 m spectral bin on average, which will allow us to measure the variations of the spectrum in longitude and measure the nightside emission spectrum for the first time. If the atmosphere of WASP-43b is clear, our observations will permit us to determine if its atmosphere has an equilibrium or disequilibrium chemical composition, eventually providing the first conclusive evidence of chemical quenching in a hot Jupiter atmosphere. If the atmosphere is cloudy, a careful retrieval analysis will allow us to identify the cloud composition.
The carbon-to-oxygen ratio (C/O) in a planet provides critical information about its primordial origins and subsequent evolution. A primordial C/O greater than 0.8 causes a carbide-dominated ...interior, as opposed to the silicate-dominated composition found on Earth; the atmosphere can also differ from those in the Solar System. The solar C/O is 0.54 (ref. 3). Here we report an analysis of dayside multi-wavelength photometry of the transiting hot-Jupiter WASP-12b (ref. 6) that reveals C/O ≥ 1 in its atmosphere. The atmosphere is abundant in CO. It is depleted in water vapour and enhanced in methane, each by more than two orders of magnitude compared to a solar-abundance chemical-equilibrium model at the expected temperatures. We also find that the extremely irradiated atmosphere (T > 2,500 K) of WASP-12b lacks a prominent thermal inversion (or stratosphere) and has very efficient day-night energy circulation. The absence of a strong thermal inversion is in stark contrast to theoretical predictions for the most highly irradiated hot-Jupiter atmospheres.
Celotno besedilo
Dostopno za:
DOBA, IJS, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Abstract We present eclipse maps of the two-dimensional thermal emission from the dayside of the hot-Jupiter WASP-43b, derived from an observation of a phase curve with the JWST MIRI/LRS instrument. ...The observed eclipse shapes deviate significantly from those expected for a planet emitting uniformly over its surface. We fit a map to this deviation, constructed from spherical harmonics up to order ℓ max = 2 , alongside the planetary, orbital, stellar, and systematic parameters. This yields a map with a meridionally averaged eastward hot-spot shift of (7.75 ± 0.36)°, with no significant degeneracy between the map and the additional parameters. We show the latitudinal and longitudinal contributions of the dayside emission structure to the eclipse shape, finding a latitudinal signal of ∼200 ppm and a longitudinal signal of ∼250 ppm. To investigate the sensitivity of the map to the method, we fix the parameters not used for mapping and derive an “eigenmap” fitted with an optimized number of orthogonal phase curves, which yields a similar map to the ℓ max = 2 map. We also fit a map up to ℓ max = 3 , which shows a smaller hot-spot shift, with a larger uncertainty. These maps are similar to those produced by atmospheric simulations. We conclude that there is a significant mapping signal which constrains the spherical harmonic components of our model up to ℓ max = 2 . Alternative mapping models may derive different structures with smaller-scale features; we suggest that further observations of WASP-43b and other planets will drive the development of more robust methods and more accurate maps.
Abstract 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/Low Resolution Spectrometer WASP-39b transmission spectrum shows an unexplained dip starting at ∼10 μ 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 Letter, 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.
Using the atmospheric structure from a 3D global radiation-hydrodynamic simulation of HD 189733b and the open-source Bayesian Atmospheric Radiative Transfer (BART) code, we investigate the difference ...between the secondary-eclipse temperature structure produced with a 3D simulation and the best-fit 1D retrieved model. Synthetic data are generated by integrating the 3D models over the Spitzer, the Hubble Space Telescope (HST), and the James Web Space Telescope (JWST) bandpasses, covering the wavelength range between 1 and 11 m where most spectroscopically active species have pronounced features. Using the data from different observing instruments, we present detailed comparisons between the temperature-pressure profiles recovered by BART and those from the 3D simulations. We calculate several averages of the 3D thermal structure and explore which particular thermal profile matches the retrieved temperature structure. We implement two temperature parameterizations that are commonly used in retrieval to investigate different thermal profile shapes. To assess which part of the thermal structure is best constrained by the data, we generate contribution functions for our theoretical model and each of our retrieved models. Our conclusions are strongly affected by the spectral resolution of the instruments included, their wavelength coverage, and the number of data points combined. We also see some limitations in each of the temperature parametrizations, as they are not able to fully match the complex curvatures that are usually produced in hydrodynamic simulations. The results show that our 1D retrieval is recovering a temperature and pressure profile that most closely matches the arithmetic average of the 3D thermal structure. When we use a higher resolution, more data points, and a parametrized temperature profile that allows more flexibility in the middle part of the atmosphere, we find a better match between the retrieved temperature and pressure profile and the arithmetic average. The Spitzer and HST simulated observations sample deep parts of the planetary atmosphere and provide fewer constraints on the temperature and pressure profile, while the JWST observations sample the middle part of the atmosphere, providing a good match with the middle and most complex part of the arithmetic average of the 3D temperature structure.
Abstract
We present a novel physically motivated, parametrized temperature model for phase-curve retrieval, able to self-consistently assess the variation in thermal structure in multidimensions. To ...develop this approach, we drew motivation from both full three-dimensional general circulation models and analytic formulations, accounting for the dominant dynamical feature of tidally locked planets, the planetary jet. Our formulation shows notable flexibility. It can generate planetary jets of various characteristics and redistribution efficiencies seen in the literature, including both standard eastward and unusual westward offset hotspots, as well as more exotic configurations for potential future observations. In our modeling scheme we utilize a tractable set of parameters efficient enough to enable future Bayesian analysis and, in addition to the resolved temperature structure, we return physical insights not yet derived from retrievals: the amplitude and the phase offset, and the location and the extent of the equatorial jet.
ABSTRACT
We present the open-source python Radiative-transfer in a Bayesian framework (pyrat bay) framework for exoplanet atmospheric modelling, spectral synthesis, and Bayesian retrieval. The ...modular design of the code allows the users to generate atmospheric 1D parametric models of the temperature, abundances (in thermochemical equilibrium or constant-with-altitude), and altitude profiles (in hydrostatic equilibrium); sample ExoMol and HITRAN line-by-line cross-sections with custom resolving power and line-wing cutoff values; compute emission or transmission spectra considering cross-sections from molecular line transitions, collision-induced absorption, Rayleigh scattering, grey clouds, and alkali resonance lines; and perform Markov chain Monte Carlo atmospheric retrievals for a given transit or eclipse data set. We benchmarked the pyrat bay framework by reproducing line-by-line sampling of ExoMol cross-sections, producing transmission and emission spectra consistent with petitradtrans models, accurately retrieving the atmospheric properties of simulated transmission and emission observations generated with taurex models, and closely reproducing aura retrieval analyses of the space-based transmission spectrum of HD 209458 b. Finally, we present a retrieval analysis of a population of transiting exoplanets, focusing on those observed in transmission with the HST WFC3/G141 grism. We found that this instrument alone can confidently identify when a data set shows H2O-absorption features; however, it cannot distinguish whether a muted H2O feature is caused by clouds, high atmospheric metallicity, or low H2O abundance. Our results are consistent with previous retrieval analyses. The pyrat bay code is available at PyPI (pip install pyratbay) and conda. The code is heavily documented (https://pyratbay.readthedocs.io) and tested to provide maximum accessibility to the community and long-term development stability.
ABSTRACT Time-correlated noise is a significant source of uncertainty when modeling exoplanet light-curve data. A correct assessment of correlated noise is fundamental to determine the true ...statistical significance of our findings. Here, we review three of the most widely used correlated-noise estimators in the exoplanet field, the time-averaging, residual-permutation, and wavelet-likelihood methods. We argue that the residual-permutation method is unsound in estimating the uncertainty of parameter estimates. We thus recommend to refrain from this method altogether. We characterize the behavior of the time averaging's rms-versus-bin-size curves at bin sizes similar to the total observation duration, which may lead to underestimated uncertainties. For the wavelet-likelihood method, we note errors in the published equations and provide a list of corrections. We further assess the performance of these techniques by injecting and retrieving eclipse signals into synthetic and real Spitzer light curves, analyzing the results in terms of the relative-accuracy and coverage-fraction statistics. Both the time-averaging and wavelet-likelihood methods significantly improve the estimate of the eclipse depth over a white-noise analysis (a Markov-chain Monte Carlo exploration assuming uncorrelated noise). However, the corrections are not perfect; when retrieving the eclipse depth from Spitzer data sets, these methods covered the true (injected) depth within the 68% credible region in only ∼45%-65% of the trials. Lastly, we present our open-source model-fitting tool, Multi-Core Markov-Chain Monte Carlo (MC3). This package uses Bayesian statistics to estimate the best-fitting values and the credible regions for the parameters for a (user-provided) model. MC3 is a Python/C code, available at https://github.com/pcubillos/MCcubed.
Heng & Tsai developed an analytic framework to obtain thermochemical-equilibrium abundances for , CO, , , , , , , and for a system with known temperature, pressure, and elemental abundances ...(hydrogen, carbon, nitrogen, and oxygen). However, the implementation of their approach can become numerically unstable under certain circumstances, leading to inaccurate solutions (e.g., C/O ≥ 1 atmospheres at low pressures). Building on their approach, we identified the conditions that prompt inaccurate solutions, and developed a new framework to avoid them, providing a reliable implementation for arbitrary values of temperature (200 to ∼2000 K), pressure ( to bar), and CNO abundances ( to ∼ × solar elemental abundances), for hydrogen-dominated atmospheres. The accuracy of our analytic framework is better than 10% for the more abundant species that have mixing fractions larger than , whereas the accuracy is better than 50% for the less abundant species. Additionally, we added the equilibrium-abundance calculation of atomic and molecular hydrogen into the system, and explored the physical limitations of this approach. Efficient and reliable tools, such as this one, are highly valuable for atmospheric Bayesian studies, which need to evaluate a large number of models. We implemented our analytic framework into the rate Python open-source package, available at https://github.com/pcubillos/rate.