We carry out a comprehensive Bayesian correlation analysis between hot halos and direct masses of supermassive black holes (SMBHs), by retrieving the X-ray plasma properties (temperature, luminosity, density, pressure, masses) over galactic to cluster scales for 85 diverse systems. We find new key scalings, with the tightest relation being the M • − T x , followed by M • − L x . The tighter scatter (down to 0.2 dex) and stronger correlation coefficient of all the X-ray halo scalings compared with the optical counterparts (as the M • − σ e ) suggest that plasma halos play a more central role than stars in tracing and growing SMBHs (especially those that are ultramassive). Moreover, M • correlates better with the gas mass than dark matter mass. We show the important role of the environment, morphology, and relic galaxies/coronae, as well as the main departures from virialization/self-similarity via the optical/X-ray fundamental planes. We test the three major channels for SMBH growth: hot/Bondilike models have inconsistent anti-correlation with X-ray halos and too low feeding; cosmological simulations find SMBH mergers as sub-dominant over most of the cosmic time and too rare to induce a central-limit-theorem effect; the scalings are consistent with chaotic cold accretion (CCA), the rain of matter condensing out of the turbulent X-ray halos that sustains a long-term self-regulated feedback loop. The new correlations are major observational constraints for models of SMBH feeding/feedback in galaxies, groups, and clusters (e.g., to test cosmological hydrodynamical simulations), and enable the study of SMBHs not only through X-rays, but also via the Sunyaev-Zel'dovich effect (Compton parameter), lensing (total masses), and cosmology (gas fractions).

Aims. We studied the star formation rate (SFR) in cosmological hydrodynamical simulations of galaxy (proto-)clusters in the redshift range 0 <  z <  4, comparing them to recent observational studies; we also investigated the effect of varying the parameters of the star formation model on galaxy properties such as SFR, star-formation efficiency, and gas fraction. Methods. We analyse a set of zoom-in cosmological hydrodynamical simulations centred on 12 clusters. The simulations are carried out with the GADGET-3 Tree-PM smoothed-particle hydro-dynamics code which includes various subgrid models to treat unresolved baryonic physics, including AGN feedback. Results. Simulations do not reproduce the high values of SFR observed within protocluster cores, where the values of SFR are underpredicted by a factor ≳4 both at z ∼ 2 and z ∼ 4. The difference arises as simulations are unable to reproduce the observed starburst population and is greater at z ∼ 2 because simulations underpredict the normalisation of the main sequence (MS) of star forming galaxies (i.e. the correlation between stellar mass and SFR) by a factor of ∼3. As the low normalisation of the MS seems to be driven by an underestimated gas fraction, it remains unclear whether numerical simulations miss starburst galaxies due to overly underpredicted gas fractions or overly low star formation efficiencies. Our results are stable against varying several parameters of the star formation subgrid model and do not depend on the details of AGN feedback. Conclusions. The subgrid model for star formation, introduced to reproduce the self-regulated evolution of quiescent galaxies, is not suitable to describe violent events like high-redshift starbursts. We find that this conclusion holds, independently of the parameter choice for the star formation and AGN models. The increasing number of multi-wavelength high-redshift observations will help to improve the current star formation model, which is needed to fully recover the observed star formation history of galaxy clusters.

We report a detailed CO(1−0) survey of a galaxy protocluster field at z = 2.16, based on 475 h of observations with the Australia Telescope Compact Array. We constructed a large mosaic of 13 individual pointings, covering an area of 21 arcmin2 and ±6500 km s−1 range in velocity. We obtained a robust sample of 46 CO(1−0) detections spanning z = 2.09 − 2.22, constituting the largest sample of molecular gas measurements in protoclusters to date. The CO emitters show an overdensity at z = 2.12 − 2.21, suggesting a galaxy super-protocluster or a protocluster connected to large-scale filaments of ∼120 cMpc in size. We find that 90% of CO emitters have distances >0.′5−4′ to the center galaxy, indicating that small area surveys would miss the majority of gas reservoirs in similar structures. Half of the CO emitters have velocities larger than escape velocities, which appears gravitationally unbound to the cluster core. These unbound sources are barely found within the R200 radius around the center, which is consistent with a picture in which the cluster core is collapsed while outer regions are still in formation. Compared to other protoclusters, this structure contains a relatively higher number of CO emitters with relatively narrow line widths and high luminosities, indicating galaxy mergers. We used these CO emitters to place the first constraint on the CO luminosity function and molecular gas density in an overdense environment. The amplitude of the CO luminosity function is 1.6 ± 0.5 orders of magnitude higher than that observed for field galaxy samples at z ∼ 2, and one order of magnitude higher than predictions for galaxy protoclusters from semi-analytical SHARK models. We derive a high molecular gas density of 0.6 − 1.3 × 109 M⊙ cMpc−3 for this structure, which is consistent with predictions for cold gas density of massive structures from hydro-dynamical DIANOGA simulations.