Revisiting the Extreme Clustering of $z \approx 4$ Quasars with Large Volume Cosmological Simulations

Abstract: Observations from wide-field quasar surveys indicate that the quasar auto-correlation length increases dramatically from $z\approx2.5$ to $z\approx4$. This large clustering amplitude at $z\approx4$ has proven hard to interpret theoretically, as it implies that quasars are hosted by the most massive dark matter halos residing in the most extreme environments at that redshift. In this work, we present a model that simultaneously reproduces both the observed quasar auto-correlation and quasar luminosity functions. The spatial distribution of halos and their relative abundance are obtained via a novel method that computes the halo mass and halo cross-correlation functions by combining multiple large-volume dark-matter-only cosmological simulations with different box sizes and resolutions. Armed with these halo properties, our model exploits the conditional luminosity function framework to describe the stochastic relationship between quasar luminosity, $L$, and halo mass, $M$. Assuming a simple power-law relation $L\propto M^\gamma$ with log-normal scatter, $\sigma$, we are able to reproduce observations at $z\sim 4$ and find that: (a) the quasar luminosity-halo mass relation is highly non-linear ($\gamma\gtrsim2$), with very little scatter ($\sigma \lesssim 0.3$ dex); (b) luminous quasars ($\log_{10} L/{\rm erg s}^{-1} \gtrsim 46.5-47$) are hosted by halos with mass $\log_{10} M/{\rm M}\odot\gtrsim 13-13.5$; and (c) the implied duty cycle for quasar activity approaches unity ($\varepsilon{\rm DC}\approx10-60\%$). We also consider observations at $z\approx2.5$ and find that the quasar luminosity-halo mass relation evolves significantly with cosmic time, implying a rapid change in quasar host halo masses and duty cycles, which in turn suggests concurrent evolution in black hole scaling relations and/or accretion efficiency.