Differentiable Modeling of Planet and Substellar Atmosphere: High-Resolution Emission, Transmission, and Reflection Spectroscopy with ExoJAX2

Abstract

Modeling based on differentiable programming holds great promise for astronomy, enabling advanced techniques such as gradient-based posterior sampling and optimization. This paradigm motivated us to develop ExoJAX (Kawahara et al. 2022), the first auto-differentiable spectrum model of exoplanets and brown dwarfs. ExoJAX directly calculates cross-sections as functions of temperature and pressure to minimize interpolation errors in high-dispersion spectra, although initial work focused on narrowband emission spectroscopy. Here, we introduce a fast, memory-efficient opacity algorithm and differentiable radiative transfer for emission, transmission, and reflection spectroscopy. In the era of data-rich JWST observations, retrieval analyses are often forced to bin high-resolution spectra due to computational bottlenecks. The new algorithm efficiently handles native-resolution data, preserving the full information content and dynamic range. The advances proposed in this paper enable broader applications, demonstrated by retrievals of GL229 B's high-dispersion emission, WASP-39 b's JWST mid-resolution transmission at original resolution (R $\sim$ 2,700), and Jupiter's reflection spectrum. We derive a C/O ratio for GL229 B consistent with its host star, constrain WASP-39 b's radial velocity from molecular line structures, and infer Jupiter's metallicity in line with previous estimates.

36 pages, 15 figures, accepted by ApJ. ExoJAX2 (version 2.0) has been released, available at https://github.com/HajimeKawahara/exojax

References in corpus (50)

• MCMC using Hamiltonian dynamics
• Precipitating Condensation Clouds in Substellar Atmospheres
• The chemical make-up of the Sun: A 2020 vision
• The Near-Infrared Spectrograph (NIRSpec) on the James Webb Space Telescope I. Overview of the instrument and its capabilities
• ExoMol molecular line lists XXX: a complete high-accuracy line list for water
• A Temperature and Abundance Retrieval Method for Exoplanet Atmospheres
• petitRADTRANS: a Python radiative transfer package for exoplanet characterization and retrieval
• The ExoMol database: molecular line lists for exoplanet and other hot atmospheres
• A Systematic Retrieval Analysis of Secondary Eclipse Spectra I: A Comparison of Atmospheric Retrieval Techniques
• ExoMol line lists III: An improved hot rotation-vibration line list for HCN and HNC
• Early Release Science of the Exoplanet WASP-39b with JWST NIRSpec G395H
• ExoMol molecular line lists XXXV: a rotation-vibration line list for hot ammonia
• Improved HCN/HNC linelist, model atmospheres and synthetic spectra for WZ Cas
• ExoMol molecular line lists - XVI: The rotation-vibration spectrum of hot H${}_2$S
• An accurate, extensive, and practical line list of methane for the HITEMP database
• Retrieving the Aerosol Complex Refractive Index using PyMieScatt: A Mie Computational Package with Visualization Capabilities
• ExoMol molecular line lists - XIV: The rotation-vibration spectrum of hot SO${}_2$
• In-orbit Performance of the Near-Infrared Spectrograph NIRSpec on the James Webb Space Telescope
• Optimal Estimation Retrievals of the Atmospheric Structure and Composition of HD 189733b from Secondary Eclipse Spectroscopy
• High-Resolution Spectroscopic Detection of TiO and Stratosphere in the Day-side of WASP-33b
• HELIOS-Retrieval: An Open-source, Nested Sampling Atmospheric Retrieval Code, Application to the HR 8799 Exoplanets and Inferred Constraints for Planet Formation
• ExoMol molecular line lists -- XXXVII: spectra of acetylene
• The Transiting Exoplanet Community Early Release Science Program for JWST
• Transiting Exoplanet Studies and Community Targets for JWST's Early Release Science Program
• Helios-r2 -- A new Bayesian, open-source retrieval model for brown dwarfs and exoplanet atmospheres
• The GAPS Programme with HARPS-N at TNG XVI: Measurement of the Rossiter-McLaughlin effect of the transiting planetary systems HAT-P-3, HAT-P-12, HAT-P-22, WASP-39 and WASP-60
• Vaex: Big Data exploration in the era of Gaia
• First Detection of Hydroxyl Radical Emission from an Exoplanet Atmosphere: High-dispersion Characterization of WASP-33b using Subaru/IRD
• Marvel analysis of the measured high-resolution rovibrational spectra of H2S
• JAX-COSMO: An End-to-End Differentiable and GPU Accelerated Cosmology Library
• Neural posterior estimation for exoplanetary atmospheric retrieval
• The Influence of H2O Pressure Broadening in High Metallicity Exoplanet Atmospheres
• Exploring and Validating Exoplanet Atmospheric Retrievals with Solar System Analog Observations
• CosmoPower-JAX: high-dimensional Bayesian inference with differentiable cosmological emulators
• A Catalogue of Exoplanet Atmospheric Retrieval Codes
• GIGA-Lens: Fast Bayesian Inference for Strong Gravitational Lens Modeling
• Auto-Differentiable Spectrum Model for High-Dispersion Characterization of Exoplanets and Brown Dwarfs
• Linearized Flux Evolution (LiFE): A Technique for Rapidly Adapting Fluxes from Full-Physics Radiative Transfer Models
• The cool brown dwarf Gliese 229 B is a close binary
• Information in the Reflected Light Spectra of Widely Separated Giant Exoplanets
• A retrieval challenge exercise for the Ariel mission
• Origin and extent of the opacity challenge for atmospheric retrievals of WASP-39 b
• Tidal Effects on the Radial Velocities of V723 Mon: Additional Evidence for a Dark $3M_{\odot }$ Companion
• The future of cosmological likelihood-based inference: accelerated high-dimensional parameter estimation and model comparison
• A spectroscopic thermometer: individual vibrational band spectroscopy with the example of OH in the atmosphere of WASP-33b
• Kernel Phase and Coronagraphy with Automatic Differentiation
• An Interpretable Machine Learning Framework for Modeling High-Resolution Spectroscopic Data
• Constraining the formation of WASP-39b using JWST transit spectroscopy
• Discovery of the Binarity of Gliese 229B, and Constraints on the System's Properties
• Probing Two-dimensional Asymmetries of an Exoplanet Atmosphere from Chromatic Transit Variation

Cited by in corpus (1)

• PyIRD: A Python-Based Data Reduction Pipeline for Subaru/IRD and REACH