The Impact of the WHIM on the IGM Thermal State Determined from the Low-z Lyman-α Forest

At $z \lesssim 1$, shock heating caused by large-scale velocity flows and possibly violent feedback from galaxy formation, converts a significant fraction of the cool gas ($T\sim 10^4$ K) in the intergalactic medium (IGM) into warm-hot phase (WHIM) with $T >10^5$K, resulting in a significant deviation from the previously tight power-law IGM temperature-density relationship, $T=T_0 (\rho / {\bar{\rho}})^{\gamma -1}$. This study explores the impact of the WHIM on measurements of the low-$z$ IGM thermal state, $[T_0,\gamma]$, based on the $b$-$N_{H I}$ distribution of the Lyman-$\alpha$ forest. Exploiting a machine learning-enabled simulation-based inference method trained on Nyx hydrodynamical simulations, we demonstrate that [$T_0$, $\gamma$] can still be reliably measured from the $b$-$N_{H I}$ distribution at $z=0.1$, notwithstanding the substantial WHIM in the IGM. To investigate the effects of different feedback, we apply this inference methodology to mock spectra derived from the IllustrisTNG and Illustris simulations at $z=0.1$. The results suggest that the underlying $[T_0,\gamma]$ of both simulations can be recovered with biases as low as $|\Delta \log(T_0/\text{K})| \lesssim 0.05$ dex, $|\Delta \gamma | \lesssim 0.1$, smaller than the precision of a typical measurement. Given the large differences in the volume-weighted WHIM fractions between the three simulations (Illustris 38\%, IllustrisTNG 10\%, Nyx 4\%) we conclude that the $b$-$N_{H I}$ distribution is not sensitive to the WHIM under realistic conditions. Finally, we investigate the physical properties of the detectable Lyman-$\alpha$ absorbers, and discover that although their $T$ and $\Delta$ distributions remain mostly unaffected by feedback, they are correlated with the photoionization rate used in the simulation.