No Arabic abstract
We present an improved, hybrid CPU-GPU atmospheric retrieval code, Helios-r2, which is applicable to medium-resolution emission spectra of brown dwarfs, in preparation for precision atmospheric spectroscopy in the era of the James Webb Space Telescope. The model is available as open-source code on the Exoclimes Simulation Platform. We subject Helios-r2 to a battery of tests of varying difficulty. The simplest test involves a mock retrieval on a forward model generated using the same radiative transfer technique, the same implementation of opacities, and the same chemistry model. The least trivial test involves a mock retrieval on synthetic spectra from the Sonora model grid, which uses a different radiative transfer technique, a different implementation of opacities, and a different chemistry model. A calibration factor, which is included to capture uncertainties in the brown dwarf radius, distance to the brown dwarf and flux calibration of the spectrum, may compensate, sometimes erroneously, for discrepancies in modeling choices and implementation. We analyze spectra of the benchmark brown dwarf GJ 570 D and the binary brown dwarf companions in the Epsilon Indi system. The retrieved surface gravities are consistent with previous studies and/or values inferred from dynamical masses (for Epsilon Indi Ba and Bb only). There remains no clear criterion on how to reject unphysical values of the retrieved brown dwarf radii. The inferred radii and corresponding masses should be taken with great caution. The retrieved carbon-to-oxygen ratios and metallicity depend on whether chemical equilibrium is assumed.
Computing and using opacities is a key part of modeling and interpreting data of exoplanetary atmospheres. Since the underlying spectroscopic line lists are constantly expanding and currently include up to ~ 10^10 - 10^11 transition lines, the opacity calculator codes need to become more powerful. Here we present major upgrades to the HELIOS-K GPU-accelerated opacity calculator and describe the necessary steps to process large line lists within a reasonable amount of time. Besides performance improvements, we include more capabilities and present a toolbox for handling different atomic and molecular data sets: from downloading and pre-processing the data to performing the opacity calculations in a user-friendly way. HELIOS-K supports line lists from ExoMol, HITRAN, HITEMP, NIST, Kurucz and VALD3. By matching the resolution of 0.1 cm^-1 and cutting length of 25 cm^-1 used by the ExoCross code for timing performance (251 seconds excluding data read-in time), HELIOS-K can process the ExoMol BT2 water line list in 12.5 seconds. Using a resolution of 0.01 cm^-1, it takes 45 seconds - equivalent to about 10^7 lines per second. As a wavenumber resolution of 0.01 cm^-1 suffices for most exoplanetary atmosphere spectroscopic calculations, we adopt this resolution in calculating opacity functions for several hundred atomic and molecular species, and make them freely available on the open-access DACE database. For the opacity calculations of the database, we use a cutting length of 100 cm^-1 for molecules and no cutting length for atoms. Our opacities are available for downloading from https://dace.unige.ch/opacityDatabase and may be visualized using https://dace.unige.ch/opacity.
Aims: ARCiS, a novel code for the analysis of exoplanet transmission and emission spectra is presented. The aim of the modelling framework is to provide a tool able to link observations to physical models of exoplanet atmospheres. Methods: The modelling philosophy chosen in this paper is to use physical and chemical models to constrain certain parameters while keeping free the parts where our physical understanding is still more limited. This approach, in between full physical modelling and full parameterisation, allows us to use the processes we understand well and parameterise those less understood. A Bayesian retrieval framework is implemented and applied to the transit spectra of a set of 10 hot Jupiters. The code contains chemistry and cloud formation and has the option for self consistent temperature structure computations. Results: The code presented is fast and flexible enough to be used for retrieval and for target list simulations for e.g. JWST or the ESA Ariel missions. We present results for the retrieval of elemental abundance ratios using the physical retrieval framework and compare this to results obtained using a parameterised retrieval setup. Conclusions: We conclude that for most of the targets considered the current dataset is not constraining enough to reliably pin down the elemental abundance ratios. We find no significant correlations between different physical parameters. We confirm that planets in our sample with a strong slope in the optical transmission spectrum are the planets where we find cloud formation to be most active. Finally, we conclude that with ARCiS we have a computationally efficient tool to analyse exoplanet observations in the context of physical and chemical models.
The last few years has seen a dramatic increase in the number of exoplanets known and in the range of methods for characterising their atmospheric properties. At the same time, new discoveries of increasingly cooler brown dwarfs have pushed down their temperature range which now extends down to Y-dwarfs of <300 K. Modelling of these atmospheres has required the development of new techniques to deal with the molecular chemistry and clouds in these objects. The atmospheres of brown dwarfs are relatively well understood, but some problems remain, in particular the behavior of clouds at the L/T transition. Observational data for exoplanet atmosphere characterization is largely limited to giant exoplanets that are hot because they are near to their star (hot Jupiters) or because they are young and still cooling. For these planets there is good evidence for the presence of CO and H2O absorptions in the IR. Sodium absorption is observed in a number of objects. Reflected light measurements show that some giant exoplanets are very dark, indicating a cloud free atmosphere. However, there is also good evidence for clouds and haze in some other planets. It is also well established that some highly irradiated planets have inflated radii, though the mechanism for this inflation is not yet clear. Some other issues in the composition and structure of giant exoplanet atmospheres such as the occurence of inverted temperature structures, the presence or absence of CO2 and CH4, and the occurrence of high C/O ratios are still the subject of investigation and debate.
The study of the composition of brown dwarf atmospheres helped to understand their formation and evolution. Similarly, the study of exoplanet atmospheres is expected to constrain their formation and evolutionary states. We use results from 3D simulations, kinetic cloud formation and kinetic ion-neutral chemistry to investigate ionisation processes which will affect their atmosphere chemistry: The dayside of super-hot Jupiters is dominated by atomic hydrogen, and not H$_2$O. Such planetary atmospheres exhibit a substantial degree of thermal ionisation and clouds only form on the nightside where lightning leaves chemical tracers (e.g. HCN) for possibly long enough to be detectable. External radiation may cause exoplanets to be enshrouded in a shell of highly ionised, H$_3^+$-forming gas and a weather-driven aurora may emerge. Brown dwarfs enable us to study the role of electron beams for the emergence of an extrasolar, weather-system driven aurora-like chemistry, and the effect of strong magnetic fields on cold atmospheric gases. Electron beams trigger the formation of H$_3^+$ in the upper atmosphere of a brown dwarf (e.g. LSR-J1835) which may react with it to form hydronium, H$_3$O$^+$, as a longer lived chemical tracer. Brown dwarfs and super-hot gas giants may be excellent candidates to search for H$_3$O$^+$ as an H$_3^+$ product.
The cloud formation process starts with the formation of seed particles, after which, surface chemical reactions grow or erode the cloud particles. We investigate which materials may form cloud condensation seeds in the gas temperature and pressure regimes (T$_{rm gas}$ = 100-2000 K, p$_{rm gas}$ = 10$^{-8}$-100 bar) expected to occur in planetary and brown dwarf atmospheres. We apply modified classical nucleation theory which requires surface tensions and vapour pressure data for each solid species, which are taken from the literature. We calculate the seed formation rates of TiO$_{2}$[s] and SiO[s] and find that they efficiently nucleate at high temperatures of T$_{rm gas}$ = 1000-1750 K. Cr[s], KCl[s] and NaCl[s] are found to efficiently nucleate across an intermediate temperature range of T$_{rm gas}$ = 500-1000 K. We find CsCl[s] may serve as the seed particle for the water cloud layers in cool sub-stellar atmospheres. Four low temperature ice species, H$_{2}$O[s/l], NH$_{3}$[s], H$_{2}$S[s/l] and CH$_{4}$[s], nucleation rates (T$_{rm gas}$ = 100-250 K) are also investigated for the coolest sub-stellar/planetary atmospheres. Our results suggest a possibly, (T$_{rm gas}$, p$_{rm gas}$) distributed hierarchy of seed particle formation regimes throughout the sub-stellar and planetary atmospheric temperature-pressure space. In order to improve the accuracy of the nucleation rate calculation, further research into the small cluster thermochemical data for each cloud species is warranted. The validity of these seed particle scenarios will be tested by applying it to more complete cloud models in the future.