X-ray emission maps and scaling relations in IllustrisTNG and MillenniumTNG: Differences between cluster and group regimes

Abstract

X-ray observations are a primary probe of the intracluster medium, widely used to infer galaxy cluster masses and scaling relations. We present and validate a pipeline to generate X-ray emission maps of galaxy groups and clusters from large cosmological simulations, and use it to study the origin of deviations from self-similarity across the group-to-cluster transition. We apply this pipeline to the Illustris-TNG300 and MillenniumTNG simulations, constructing \mbox{X-ray} emission maps and spectra based on APEC cooling functions. For simulations that do not explicitly track individual chemical abundances, we introduce a metallicity-based prescription that accurately reproduces the full spectral emission. We derive the $L_{\mathrm{X}}$--$M_{500}$, $M_{\mathrm{gas}}$--$M_{500}$, and $T$--$M_{500}$ scaling relations over $10^{12.5} \leq M_{500} \leq 10^{15.5}\,\mathrm{M}_\odot$, compare them with observational data, and quantify the hydrostatic equilibrium and spectroscopic temperature biases through synthetic X-ray analyses. The simulated scaling relations are in good overall agreement with observations and are best described by broken power laws with a pivot at $M_{500}=10^{13.67}\,\mathrm{M}_\odot$. At high masses, the slopes are close to self-similar expectations; at lower masses the relations steepen significantly, reflecting the growing importance of AGN feedback. X-ray hydrostatic masses are systematically underestimated by $15\%$, independently of cluster mass. When spectroscopic effects are included, the bias becomes mass-dependent, ranging from $15\%$ at low masses to $21\%$ at high masses. The recovered X-ray luminosity, measured in the $0.15$--$1\,R_{500}$ aperture, is also mass-dependent: high-mass clusters are underestimated by $18\%$, while low-mass systems show discrepancies of up to $33\%$, driven by single-temperature spectral modelling of the gas outside the core.