Auto-Differentiable Spectrum Model for High-Dispersion Characterization of Exoplanets and Brown Dwarfs

Abstract

We present an auto-differentiable spectral modeling of exoplanets and brown dwarfs. This model enables a fully Bayesian inference of the high--dispersion data to fit the ab initio line-by-line spectral computation to the observed spectrum by combining it with the Hamiltonian Monte Carlo in recent probabilistic programming languages. An open source code, exojax, developed in this study, was written in Python using the GPU/TPU compatible package for automatic differentiation and accelerated linear algebra, JAX (Bradbury et al. 2018). We validated the model by comparing it with existing opacity calculators and a radiative transfer code and found reasonable agreements of the output. As a demonstration, we analyzed the high-dispersion spectrum of a nearby brown dwarf, Luhman 16 A and found that a model including water, carbon monoxide, and ${\mathrm{H}}_2/He$ collision induced absorption was well fitted to the observed spectrum ($R=10^5$ and 2.28-2.30 $\mu$m). As a result, we found that $T_0=1295_{-32}^{+35}$ K at 1 bar and C/O $=0.62\pm 0.03$, which is slightly higher than the solar value. This work demonstrates the potential of full Bayesian analysis of brown dwarfs and exoplanets as observed by high-dispersion spectrographs and also directly-imaged exoplanets as observed by high-dispersion coronagraphy.

arXiv comment29 pages, 18 figures, Accepted by ApJS. ExoJAX can now handle a large number of lines, thanks to a new rapid opacity calculator. Available at https://github.com/HajimeKawahara/exojax

References in corpus (79)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• The NumPy array: a structure for efficient numerical computation
• The chemical composition of the Sun
• The No-U-Turn Sampler: Adaptively Setting Path Lengths in Hamiltonian Monte Carlo
• Efficient, uninformative sampling of limb darkening coefficients for two-parameter laws
• The orbital motion, absolute mass, and high-altitude winds of exoplanet HD209458b
• The chemical make-up of the Sun: A 2020 vision
• ExoMol molecular line lists XXX: a complete high-accuracy line list for water
• petitRADTRANS: a Python radiative transfer package for exoplanet characterization and retrieval
• The ExoMol database: molecular line lists for exoplanet and other hot atmospheres
• A Conceptual Introduction to Hamiltonian Monte Carlo
• The signature of orbital motion from the dayside of the planet tau Bootis b
• The fast spin-rotation of a young extrasolar planet
• Detection of water absorption in the day side atmosphere of HD 189733 b using ground-based high-resolution spectroscopy at 3.2 microns
• A Systematic Retrieval Analysis of Secondary Eclipse Spectra I: A Comparison of Atmospheric Retrieval Techniques
• Nightside condensation of iron in an ultra-hot giant exoplanet
• Retrieving Temperatures and Abundances of Exoplanet Atmospheres with High-Resolution Cross-Correlation Spectroscopy
• Multi-Wavelength Constraints on the Day-Night Circulation Patterns of HD 189733b
• Combining high-dispersion spectroscopy (HDS) with high contrast imaging (HCI): Probing rocky planets around our nearest neighbors
• Atomic iron and titanium in the atmosphere of the exoplanet KELT-9b
• Composable Effects for Flexible and Accelerated Probabilistic Programming in NumPyro
• A spectral survey of an ultra-hot Jupiter: Detection of metals in the transmission spectrum of KELT-9 b
• Transmission spectral properties of clouds for hot Jupiter exoplanets
• exoplanet: Gradient-based probabilistic inference for exoplanet data & other astronomical time series
• Detection of Fe I in the atmosphere of the ultra-hot Jupiter WASP-121b, and a new likelihood-based approach for Doppler-resolved spectroscopy
• A global cloud map of the nearest known brown dwarf
• Discovery of a Binary Brown Dwarf at 2 Parsecs from the Sun
• Atmospheric characterization of the ultra-hot Jupiter MASCARA-2b/KELT-20b
• Five carbon- and nitrogen-bearing species in a hot giant planet's atmosphere
• High-Resolution Spectroscopic Detection of TiO and Stratosphere in the Day-side of WASP-33b
• Weighing The Non-Transiting Hot Jupiter Tau BOO b
• Dust in the Photospheric Environment: Unified Cloudy Models of M, L, and T Dwarfs
• First direct detection of an exoplanet by optical interferometry; Astrometry and K-band spectroscopy of HR8799 e
• HELIOS-Retrieval: An Open-source, Nested Sampling Atmospheric Retrieval Code, Application to the HR 8799 Exoplanets and Inferred Constraints for Planet Formation
• HELIOS-K 2.0 Opacity Calculator and Open-source Opacity Database for Exoplanetary Atmospheres
• Discovery of water at high spectral resolution in the atmosphere of 51 Peg b
• Observing Exoplanets with High Dispersion Coronagraphy. I. The scientific potential of current and next-generation large ground and space telescopes
• A framework to combine low- and high-resolution spectroscopy for the atmospheres of transiting exoplanets
• Where is the Water? Jupiter-like C/H ratio but strong H${}_2$O depletion found on $\tau$ Boötis b using SPIRou
• Ionized calcium in the atmospheres of two ultra-hot exoplanets WASP-33b and KELT-9b
• Theoretical transmission spectra of exoplanet atmospheres with hydrocarbon haze: Effect of creation, growth, and settling of haze particles. I. Model description and first results
• Observing Exoplanets with High-Dispersion Coronagraphy. II. Demonstration of an Active Single-Mode Fiber Injection Unit
• Medium-resolution integral-field spectroscopy for high-contrast exoplanet imaging: Molecule maps of the $\beta$ Pictoris system with SINFONI
• The slow spin of the young sub-stellar companion GQ Lupi b and its orbital configuration
• Detection of neutral atomic species in the ultra-hot jupiter WASP-121b
• Exoplanet atmospheres with GIANO II. Detection of molecular absorption in the dayside spectrum of HD 102195b
• Analytical Models of Exoplanetary Atmospheres. II. Radiative Transfer via the Two-stream Approximation
• Resolved Near-Infrared Spectroscopy of WISE J104915.57-531906.1AB: A Flux-Reversal Binary at the L dwarf/T dwarf Transition
• Evidence for multiple molecular species in the hot Jupiter HD 209458b
• Evidence against a strong thermal inversion in HD 209458 b from high-dispersion spectroscopy
• Detection of ionized calcium in the atmosphere of the ultra-hot Jupiter KELT-9b
• A temperature inversion with atomic iron in the ultra-hot dayside atmosphere of WASP-189b
• First Detection of Hydroxyl Radical Emission from an Exoplanet Atmosphere: High-dispersion Characterization of WASP-33b using Subaru/IRD
• Detection of Fe and evidence for TiO in the dayside emission spectrum of WASP-33b
• Cloud structure of the nearest brown dwarfs: Spectroscopic variability of Luhman 16AB from the Hubble Space Telescope
• Detection of Fe I and Fe II in the atmosphere of MASCARA-2b using a cross-correlation method
• High-resolution Transmission Spectroscopy of MASCARA-2 b with EXPRES
• Detection of Fe\,{\sc i} Emission in the Day-side Spectrum of WASP-33b
• HyDRA-H: Simultaneous Hybrid Retrieval of Exoplanetary Emission Spectra
• Neutral Cr and V in the atmosphere of ultra hot jupiter WASP-121 b
• The Influence of H2O Pressure Broadening in High Metallicity Exoplanet Atmospheres
• Detecting Exomoons Via Doppler Monitoring of Directly Imaged Exoplanets
• Direct characterization of young giant exoplanets at high spectral resolution by coupling SPHERE and CRIRES+
• Updated astrometry and masses of the LUH 16 brown dwarf binary
• CO Detected in CI Tau b: Hot Start Implied by Planet Mass and M${}_K$
• The Detectability and Constraints of Biosignature Gases in the Near & Mid-Infrared from Transit Transmission Spectroscopy
• Individual, Model-Independent Masses of the Closest Known Brown Dwarf Binary to the Sun
• High-Resolution Infrared Spectroscopy of the Brown Dwarf Epsilon Indi Ba
• Spectroscopic Coronagraphy for Planetary Radial Velocimetry of Exoplanets
• Neutral metals in the atmosphere of HD149026b
• Tidal Effects on the Radial Velocities of V723 Mon: Additional Evidence for a Dark $3M_{\odot }$ Companion
• Kernel Phase and Coronagraphy with Automatic Differentiation
• The Keck Planet Imager and Characterizer: Demonstrating advanced exoplanet characterization techniques for future extremely large telescopes
• The Spin Effect on Planetary Radial Velocimetry of Exoplanets
• Algorithm 916: computing the Faddeyeva and Voigt functions
• JAXNS: a high-performance nested sampling package based on JAX
• Characterizing Earth-like Planets Using a Combination of High-Dispersion Spectroscopy and High-Contrast Instruments: Doppler-shifted Water and Oxygen Lines
• Efficient Multi-Accuracy Computations of Complex Functions with Complex Arguments
• Occultation mapping of Io's surface in the near-infrared I: Inferring static maps

Cited by in corpus (15)

• The ESO SupJup Survey I: Chemical and isotopic characterisation of the late L-dwarf DENIS J0255-4700 with CRIRES${}^{+}$
• A Catalogue of Exoplanet Atmospheric Retrieval Codes
• The ESO SupJup Survey II: The 12C/13C ratios of three young brown dwarfs with CRIRES${}^{+}$
• The ESO SupJup Survey IV. Unveiling the carbon isotope ratio of GQ Lup B and its host star
• The JWST Weather Report: retrieving temperature variations, auroral heating, and static cloud coverage on SIMP-0136
• An Interpretable Machine Learning Framework for Modeling High-Resolution Spectroscopic Data
• Multi-modal atmospheric characterization of $\beta$ Pictoris b: Adding high-resolution continuum spectra from GRAVITY
• PICASO 4.0: Clouds and Photochemistry in Climate Models of Brown Dwarfs and Exoplanets
• The ESO SupJup Survey VIII. Chemical fingerprints of young L dwarf twins
• Radiative Transfer and Inversion codes for characterizing planetary atmospheres: an overview
• Probing Two-dimensional Asymmetries of an Exoplanet Atmosphere from Chromatic Transit Variation
• Differentiable Modeling of Planet and Substellar Atmosphere: High-Resolution Emission, Transmission, and Reflection Spectroscopy with ExoJAX2
• Analytic auto-differentiable $\Lambda$CDM cosmography
• A 1.48-2.48 $\mu$m R=28000 spectroscopic atlas of the L7.5 and T0.5 components of the nearest pair of brown dwarfs: Luhman 16AB
• Atmospheric retrieval of Subaru/IRD high-resolution spectrum of the archetype T-type brown dwarf Gl 229 B