GASTLI: An open-source coupled interior-atmosphere model to unveil gas giant composition

Abstract

The metal mass fractions of gas giants are a powerful tool to constrain their formation mechanisms and evolution. The metal content is inferred by comparing mass and radius measurements with interior structure and evolution models. In the midst of the JWST, CHEOPS, TESS, and the forthcoming PLATO era, we are at the brink of obtaining unprecedented precision in radius, age and atmospheric metallicity measurements. To prepare for this wealth of data, we present the GAS gianT modeL for Interiors (GASTLI), an easy-to-use, publicly available Python package. The code is optimized to rapidly calculate mass-radius relations, and radius and luminosity thermal evolution curves for a variety of envelope compositions and core mass fractions. Its applicability spans planets with masses $17M_{\oplus }<M<6M_{Jup}$, and equilibrium temperatures $T_{eq}<1000$ K. The interior model is stratified in a core composed of water and rock, and an envelope constituted by H/He and metals (water). The interior is coupled to a grid of self-consistent, cloud-free atmospheric models to determine the atmospheric and boundary interior temperature, as well as the contribution of the atmosphere to the total radius. We successfully validate GASTLI by comparing it to previous work and data of the Solar System's gas giants and Neptune. We also test GASTLI on the Neptune-mass exoplanet HAT-P-26 b, finding a bulk metal mass fraction between 0.60-0.78 and a core mass of 8.5-14.4 $M_{\oplus }$. Finally, we explore the impact of different equations of state and assumptions, such as C/O ratio and transit pressure, in the estimation of bulk metal mass fraction. These differences between interior models entail a change in radius of up to 2.5% for Jupiter-mass planets, but more than 10\% for Neptune-mass. These are equivalent to variations in core mass fraction of 0.07, or 0.10 in envelope metal mass fraction.

18 pages, 9 figures. In review in Astronomy & Astrophysics

References in corpus (51)

• Modules for Experiments in Stellar Astrophysics (MESA)
• Modules for Experiments in Stellar Astrophysics (MESA): Pulsating Variable Stars, Rotation, Convective Boundaries, and Energy Conservation
• Growth Model Interpretation of Planet Size Distribution
• On the radiative equilibrium of irradiated planetary atmospheres
• On the Luminosity of Young Jupiters
• Comparing Jupiter interior structure models to Juno gravity measurements and the role of a dilute core
• Planetary population synthesis coupled with atmospheric escape: a statistical view of evaporation
• Can we constrain interior structure of rocky exoplanets from mass and radius measurements?
• Inflating Hot Jupiters With Ohmic Dissipation
• Ab initio Equation of State data for hydrogen, helium, and water and the internal structure of Jupiter
• Methane, Carbon Monoxide, and Ammonia in Brown Dwarfs and Self-Luminous Giant Planets
• HAT-P-26b: A Neptune-Mass Exoplanet with a Well Constrained Heavy Element Abundance
• Observing transiting planets with JWST -- Prime targets and their synthetic spectral observations
• A new equation of state for dense hydrogen-helium mixtures
• A Systematic Study of Departures from Chemical Equilibrium in the Atmospheres of Substellar Mass Objects
• The water abundance in Jupiter's equatorial zone
• The interior and atmosphere of the habitable-zone exoplanet K2-18b
• Beyond Equilibrium Temperature: How the Atmosphere/Interior Connection Affects the Onset of Methane, Ammonia, and Clouds in Warm Transiting Giant Planets
• A new equation of state for dense hydrogen-helium mixtures II: taking into account hydrogen-helium interactions
• Jupiter's formation and its primordial internal structure
• Toward the analysis of JWST exoplanet spectra: Identifying troublesome model parameters
• Constraints on Super-Earths Interiors from Stellar Abundances
• Idealised simulations of the deep atmosphere of hot jupiters: Deep, hot, adiabats as a robust solution to the radius inflation problem
• Convection and Mixing in Giant Planet Evolution
• Models of Saturn's Interior Constructed with Accelerated Concentric Maclaurin Spheroid Method
• Evidence for a Dichotomy in the Interior Structures of Jupiter and Saturn from Helium Phase Separation
• Jupiter's interior from Juno: Equation-of-state uncertainties and dilute core extent
• Characterisation of the hydrospheres of TRAPPIST-1 planets
• Effect of Non-Adiabatic Thermal Profiles on the Inferred Compositions of Uranus and Neptune
• The Metal-Rich Atmosphere of the Neptune HAT-P-26b
• Do Metal-Rich Stars Make Metal-Rich Planets? New Insights on Giant Planet Formation from Host Star Abundances
• Planet formation around M dwarfs via disc instability: Fragmentation conditions and protoplanet properties
• Detection of the tidal deformation of WASP-103b at $3\sigma$ with CHEOPS
• Latitudinal variation in the abundance of methane (CH4) above the clouds in Neptune's atmosphere from VLT/MUSE Narrow Field Mode Observations
• An estimate of the $k_2$ Love number of WASP-18Ab from its radial velocity measurements
• Possible Internal Structures and Compositions of Proxima Centauri b
• Accounting for non-ideal mixing effects in the hydrogen-helium equation of state
• Synthetic Evolution Tracks of Giant Planets
• MAGRATHEA: an open-source spherical symmetric planet interior structure code
• A subsolar oxygen abundance or a radiative region deep in Jupiter revealed by thermochemical modelling
• Super-adiabatic Temperature Gradient at Jupiter's Equatorial Zone and Implications for the Water Abundance
• No evidence for radius inflation in hot Jupiters from vertical advection of heat
• Latitudinal variation of methane mole fraction above clouds in Neptune's atmosphere from VLT/MUSE-NFM: Limb-darkening reanalysis
• A new method to measure the spectra of transiting exoplanet atmospheres using multi-object spectroscopy
• Water content trends in K2-138 and other low-mass multiplanetary systems
• On the Effects of Planetary Oblateness on Exoplanet Studies
• HST/STIS capability for Love number measurement of WASP-121b
• Evidence of apsidal motion and a possible co-moving companion star detected in the WASP-19 system
• Connecting gravity field, moment of inertia, and core properties in Jupiter through empirical structure models
• The effect of cloudy atmospheres on the thermal evolution of warm giant planets from an interior modelling perspective
• Constraining the oblateness of transiting planets with photometry and spectroscopy