Measuring the thermal and ionization state of the low-$z$ IGM using likelihood free inference

Abstract

We present a new approach to measure the power-law temperature density relationship $T=T_0(\rho /\bar{\rho }{)}^{\gamma -1}$ and the UV background photoionization rate ${\Gamma }_{\mathrm{HI}}$ of the IGM based on the Voigt profile decomposition of the Ly$\alpha$ forest into a set of discrete absorption lines with Doppler parameter $b$ and the neutral hydrogen column density $N_{\mathrm{HI}}$. Previous work demonstrated that the shape of the $b$-$N_{\mathrm{HI}}$ distribution is sensitive to the IGM thermal parameters $T_0$ and $\gamma$, whereas our new inference algorithm also takes into account the normalization of the distribution, i.e. the line-density d$N$/d$z$, and we demonstrate that precise constraints can also be obtained on ${\Gamma }_{\mathrm{HI}}$. We use density-estimation likelihood-free inference (DELFI) to emulate the dependence of the $b$-$N_{\mathrm{HI}}$ distribution on IGM parameters trained on an ensemble of 624 Nyx hydrodynamical simulations at $z=0.1$, which we combine with a Gaussian process emulator of the normalization. To demonstrate the efficacy of this approach, we generate hundreds of realizations of realistic mock HST/COS datasets, each comprising 34 quasar sightlines, and forward model the noise and resolution to match the real data. We use this large ensemble of mocks to extensively test our inference and empirically demonstrate that our posterior distributions are robust. Our analysis shows that by applying our new approach to existing Ly$\alpha$ forest spectra at $z\simeq 0.1$, one can measure the thermal and ionization state of the IGM with very high precision (${\sigma }_{\log T_0}\sim 0.08$ dex, ${\sigma }_{\gamma }\sim 0.06$, and ${\sigma }_{\log {\Gamma }_{\mathrm{HI}}}\sim 0.07$ dex).

20 pages, 17 figures, accepted for publication in MNRAS

References in corpus (23)

• A Cosmological Framework for the Co-Evolution of Quasars, Supermassive Black Holes, and Elliptical Galaxies: I. Galaxy Mergers & Quasar Activity
• Simulating galaxy formation with black hole driven thermal and kinetic feedback
• Introducing the Illustris Project: the evolution of galaxy populations across cosmic time
• Cosmic Reionization and Early Star-Forming Galaxies: A Joint Analysis of New Constraints from Planck and Hubble Space Telescope
• A unified model for AGN feedback in cosmological simulations of structure formation
• Model-independent evidence in favor of an end to reionization by z~6
• HeII Reionization and its Effect on the IGM
• Evolution of the Intergalactic Opacity: Implications for the Ionizing Background, Cosmic Star Formation, and Quasar Activity
• Fast likelihood-free cosmology with neural density estimators and active learning
• Possible evidence for an inverted temperature-density relation in the intergalactic medium from the flux distribution of the Lyman-alpha forest
• Voigt Profile Fitting to Quasar Absorption Lines: An Analytic Approximation to the Voigt-Hjerting Function
• The thermal history of the intergalactic medium down to redshift z=1.5: a new curvature measurement
• The End of Helium Reionization at z~2.7 Inferred from Cosmic Variance in HST/COS HeII Lyman Alpha Absorption Spectra
• IGM Constraints from the SDSS-III/BOSS DR9 Ly-alpha Forest Flux Probability Distribution Function
• A New Precision Measurement of the Small-Scale Line-of-Sight Power Spectrum of the Ly$\alpha$ Forest
• VoIgt profile Parameter Estimation Routine (VIPER): H I photoionization rate at z<0.5
• A new Measurement of the Intergalactic Temperature at $z\sim 2.55-2.95$
• The evolution of the helium-ionizing background at z ~ 2-3
• The broadening of Lyman-alpha forest absorption lines
• The Effect Of AGN Heating On The Low-Redshift Lyα Forest
• The effect of stellar and AGN feedback on the low redshift Lyman-$\alpha$ forest in the Sherwood simulation suite
• A Novel Statistical Method for Measuring the Temperature-Density Relation in the IGM Using the $b$-$N_{\text{HI}}$ Distribution of Absorbers in the Ly$\alpha$ Forest
• Improving IGM temperature constraints using wavelet analysis on high-redshift quasars