Baryonic assembly bias in X-ray-selected galaxy groups and clusters: insights from the Magneticum simulation

Galaxy groups and clusters trace the large-scale matter distribution, with their clustering usually interpreted mainly as a function of halo mass. Yet, at fixed mass, their baryonic properties retain information about halo growth, gas accretion, and feedback. The intrinsic scatter in X-ray luminosity and gas fraction suggests that X-ray-selected systems may not be a random subset of the halo population. If these observables correlate with halo assembly, they may trace secondary variations in halo bias. We test this using the Magneticum hydrodynamical simulation, measuring the clustering of systems selected by X-ray luminosity and gas fraction at fixed halo mass. We construct mass-matched subsamples by ranking halos in percentiles of X-ray luminosity and derive the linear halo-matter bias from the halo-matter cross-power spectrum. X-ray-bright halos are more strongly clustered than X-ray-faint halos at fixed mass. For the 84th-16th percentile split, we find $\Delta b_{\rm lin}=0.17\pm0.03$, corresponding to a $\sim17\%$ enhancement relative to the X-ray-faint sample. A 67th-33rd split gives a consistent signal, with $\Delta b_{\rm lin}=0.12\pm0.02$ and a $\sim12\%$ enhancement. The effect is strongest at group scales and negligible for cluster-size halos. Gas fraction shows an even stronger clustering dependence, with relative enhancements of $\sim39\%$ and $\sim26\%$ for the two percentile splits. This signal is present from $z\simeq2$, whereas X-ray luminosity becomes significant only at $z\simeq0.3$, once the gas thermodynamic state is more closely coupled to baryon retention. Matching halos by both mass and formation time reduces the large-scale bias difference to below $2\sigma$, indicating that formation time captures much of the signal. These results show that, in Magneticum, X-ray luminosity traces a baryonic manifestation of halo assembly bias beyond mass.