2025/08/12 by Anirban Chowdhary, Chowdhary, Anirban, Suchetana Chatterjee +1 · 1 citation
Engineering · Physics and Astronomy · #Astronomy and Astrophysical Research #Astrophysics of Galaxies (astro-ph.GA) #Cosmology and Nongalactic Astrophysics (astro-ph.CO) #FOS: Physical sciences #Galaxies: Formation, Evolution, Phenomena #Space Technology and Applications
paper · pdf · doi:10.48550/arxiv.2508.08851
openalex publication_date 2025/08/12 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
We use cosmological hydrodynamic simulations (IllustrisTNG and SIMBA) to explore the redshift, luminosity, and black hole mass dependence of the quasar halo occupation distribution (HOD). In both simulations, we find that the quasar activity is quenched at a characteristic halo mass (∼ 1013 M\odot) scale affecting the nature of its occupation distribution function. We note that the quenching is more pronounced at low redshifts for quasars selected through a luminosity threshold. We show that a very significant bias (a factor of ∼ 10-50 in the central occupation and ∼ 10-70% in the satellite occupation fraction) is introduced in the reconstruction of quasar host halo mass distributions from the observed two-point-correlation function, if the HOD modeling does not account for the quenching effect in the central occupation function. While there is strong suppression of the occupation fraction of central quasars, the satellite occupation still follows a power-law like behavior. Our results show that the global satellite fraction of quasars increases monotonically from high to low redshifts, with 20-40 % of the quasars being satellite at intermediate redshifts, consistent with previous clustering based estimates. In addition, our study reveals that while the occupation function of quasars depends on redshift, luminosity, and feedback modes, there is hardly any evolution in the supermassive black hole (SMBH; mass-selected sample) occupation. The SMBH HOD over the entire parameter space is well-modeled by a power-law and a step function similar to what has been found for galaxies and low-luminosity active galactic nuclei.