Feedback reshapes the baryon distribution within haloes, in halo outskirts, and beyond: the closure radius from dwarfs to massive clusters

We explore three sets of cosmological hydrodynamical simulations, IllustrisTNG (TNG50, TNG100, and TNG300), EAGLE, and SIMBA, to investigate the physical processes impacting the distribution of baryons in and around haloes across an unprecedented mass range of $10^8\lt M_{\rm 200c}/\rm {{\rm M}_{\odot }}\lt 10^{15}$, from the halo centre out to scales as large as $30\, R_{\rm 200c}$. We demonstrate that baryonic feedback mechanisms significantly redistribute gas, lowering the baryon fractions inside haloes while simultaneously accumulating this material outside the virial radius. To understand this large-scale baryonic redistribution and identify the dominant physical processes responsible, we examine several variants of TNG that selectively exclude stellar and active galactic nucleus (AGN) feedback, cooling, and radiation. We find that heating from the UV background in low-mass haloes ($10^{8}\le {M_{\rm 200c}}/\rm {{\rm M}_{\odot }}\lt 10^{10}$), stellar feedback in intermediate-mass haloes ($10^{10}\le {M_{\rm 200c}}/\rm {{\rm M}_{\odot }}\lt 10^{12}$), and AGN feedback in groups ($10^{12}\le {M_{\rm 200c}}/\rm {{\rm M}_{\odot }}\lt 10^{14}$) are the dominant processes. Galaxy clusters ($M_{\rm 200c}/\rm {{\rm M}_{\odot }}\gt 10^{14}$) are the least influenced by these processes on large scales. We define a halo mass-dependent characteristic scale, the closure radius R_c, within which all baryons associated with haloes are found. For groups and clusters, we introduce a universal relation between this scale and the halo baryon fraction: R_c/R_{200c, 500c} - 1 = β(z)(1 - f_b(< R_{200c, 500c})/f_{b, cosmic}), where $\beta (z)=\alpha \, (1+z)^\gamma$, and α and γ are free parameters fit using the simulations. Accordingly, we predict that all baryons associated with observed X-ray haloes can be found within R_c ~ 1.5-2.5R_200c. Our results can be used to constrain theoretical models, particularly the physics of supernova and AGN feedback, as well as their interplay with environmental processes, through comparison with current and future X-ray and Sunyaev-Zel'dovich (SZ) observations.