High Resolution Grid-based Simulations of the Warm-Hot Intergalactic Medium

We present high-resolution cosmological hydrodynamic simulations of the Warm-Hot Intergalactic Medium (WHIM) using the GPU-optimized grid-based hydrodynamic code Kratos. Employing a uniform 4096 cubed grid in a (100 h-1 Mpc) cubed comoving volume, we achieve a spatial resolution of 24.5 h-1 kpc, sufficient to resolve the Jeans scale of gas at T 10 4 K and n H 10-3 to 10-2 cm-3. This calculation ranks among the largest grid-based cosmological hydrodynamic simulations performed to date. We find that 23.4% of cosmic baryons reside in the WHIM phase (T = 10 5 - 10 7 K) at z = 0, significantly below the 40-50% found in earlier, lower-resolution simulations. Through a suite of lower-resolution simulations, we demonstrate that spatial resolution plays a pivotal role in determining the WHIM fraction: resolving gas near its Jeans scale allows it to reach higher densities, where enhanced radiative cooling transfers a substantial fraction of baryons out of the WHIM temperature range. The remaining WHIM resides predominantly in filaments and accretion-shock structures in the vicinity of halos, where hierarchical structure formation and ongoing hydrodynamic accretion provide continued shock heating. Synthetic observations of Lyman-alpha, O VI, O VII, and O VIII emission reveal distinct morphological and kinematic signatures, with O VI tracing filament-halo interfaces and the X-ray lines probing hotter gas associated with massive halos. These predictions underscore the importance of current and future missions such as XRISM, ATHENA, and HUBS for mapping WHIM thermodynamics and kinematics.