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

Abstract

We present a new method for determining the thermal state of the intergalactic medium based on Voigt profile decomposition of the Ly$\alpha$ forest. The distribution of Doppler parameter and column density ($b$-$N_{\text{HI}}$ distribution) is sensitive to the temperature density relation $T=T_0(\rho /{\rho }_0{)}^{\gamma -1}$, and previous work has inferred $T_0$ and $\gamma$ by fitting its low-$b$ cutoff. This approach discards the majority of available data, and is susceptible to systematics related to cutoff determination. We present a method that exploits all information encoded in the $b$-$N_{\text{HI}}$ distribution by modeling its entire shape. We apply kernel density estimation to discrete absorption lines to generate model probability density functions, then use principal component decomposition to create an emulator which can be evaluated anywhere in thermal parameter space. We introduce a Bayesian likelihood based on these models enabling parameter inference via Markov chain Monte Carlo. The method's robustness is tested by applying it to a large grid of thermal history simulations. By conducting 160 mock measurements we establish that our approach delivers unbiased estimates and valid uncertainties for a 2D $(T_0,\gamma )$ measurement. Furthermore, we conduct a pilot study applying this methodology to real observational data at $z=2$. Using 200 absorbers, equivalent in pathlength to a single Ly$\alpha$ forest spectrum, we measure $\log T_0=4.092_{-0.055}^{+0.050}$ and $\gamma =1.49_{-0.074}^{+0.073}$ in excellent agreement with cutoff fitting determinations using the same data. Our method is far more sensitive than cutoff fitting, enabling measurements of $\log T_0$ and $\gamma$ with precision on $\log T_0$ ($\gamma$) nearly two (three) times higher for current dataset sizes.

Accepted for publication in ApJ, 23 pages, 15 figures

References in corpus (21)

• Cosmic Reionization and Early Star-Forming Galaxies: A Joint Analysis of New Constraints from Planck and Hubble Space Telescope
• 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
• A Direct Precision Measurement of the Intergalactic Lyman-alpha Opacity at 2<z<4.2
• Possible evidence for an inverted temperature-density relation in the intergalactic medium from the flux distribution of the Lyman-alpha forest
• Cosmic Calibration: Constraints from the Matter Power Spectrum and the Cosmic Microwave Background
• The thermal history of the intergalactic medium down to redshift z=1.5: a new curvature measurement
• Self-Consistent Modeling of Reionization in Cosmological Hydrodynamical Simulations
• An unbiased measurement of the UV background and its evolution via the proximity effect in quasar spectra
• The End of Helium Reionization at z~2.7 Inferred from Cosmic Variance in HST/COS HeII Lyman Alpha Absorption Spectra
• Cosmic Calibration
• IGM Constraints from the SDSS-III/BOSS DR9 Ly-alpha Forest Flux Probability Distribution Function
• The Second Data Release of the KODIAQ Survey
• 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
• AGN-driven helium reionization and the incidence of extended HeIII regions at redshift z>3
• The effect of stellar and AGN feedback on the low redshift Lyman-$\alpha$ forest in the Sherwood simulation suite