Mapping the Growth of Supermassive Black Holes as a Function of Galaxy Stellar Mass and Redshift

Fan Zou,Zhibo Yu,W. N. Brandt,Hyungsuk Tak,Guang Yang,Qingling Ni

DOI: https://doi.org/10.3847/1538-4357/ad27cc

2024-03-30

Abstract:The growth of supermassive black holes is strongly linked to their galaxies. It has been shown that the population mean black-hole accretion rate ($\overline{\mathrm{BHAR}}$) primarily correlates with the galaxy stellar mass ($M_\star$) and redshift for the general galaxy population. This work aims to provide the best measurements of $\overline{\mathrm{BHAR}}$ as a function of $M_\star$ and redshift over ranges of $10^{9.5}<M_\star<10^{12}~M_\odot$ and $z<4$. We compile an unprecedentedly large sample with eight thousand active galactic nuclei (AGNs) and 1.3 million normal galaxies from nine high-quality survey fields following a wedding-cake design. We further develop a semiparametric Bayesian method that can reasonably estimate $\overline{\mathrm{BHAR}}$ and the corresponding uncertainties, even for sparsely populated regions in the parameter space. $\overline{\mathrm{BHAR}}$ is constrained by X-ray surveys sampling the AGN accretion power and UV-to-infrared multi-wavelength surveys sampling the galaxy population. Our results can independently predict the X-ray luminosity function (XLF) from the galaxy stellar mass function (SMF), and the prediction is consistent with the observed XLF. We also try adding external constraints from the observed SMF and XLF. We further measure $\overline{\mathrm{BHAR}}$ for star-forming and quiescent galaxies and show that star-forming $\overline{\mathrm{BHAR}}$ is generally larger than or at least comparable to the quiescent $\overline{\mathrm{BHAR}}$.

Astrophysics of Galaxies

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: **How does the growth of super - massive black holes (SMBHs) depend on the properties of galaxies, especially the stellar mass and redshift of galaxies?**

Specifically, the goal of the paper is to provide the best measurements of the black hole accretion rate (BHAR) as a function of the stellar mass and redshift of galaxies in the range of stellar mass \(10^{9.5}<M_{\star}<10^{12}M_{\odot}\) and redshift \(z < 4\). The paper achieves this goal by compiling an unprecedentedly large sample, including 8,000 active galactic nuclei (AGNs) and 1.3 million normal galaxies, and developing a semi - parametric Bayesian method to reasonably estimate BHAR and its uncertainty.

### Main research contents:
1. **Data collection**:
- Compiled samples from nine high - quality survey fields, which cover a wide range from the nearby universe to high redshift.
- Used X - ray surveys to sample the accretion power of AGNs and used ultraviolet - to - infrared multi - band surveys to sample the galaxy population.

2. **Method development**:
- Developed a semi - parametric Bayesian method that can also reasonably estimate BHAR and its uncertainty in sparse regions of the parameter space.
- Considered the detection capabilities of different surveys through the X - ray detection probability function \(P_{\text{det}}(f_X)\), where \(f_X\) is the X - ray flux in the 2 - 10 keV energy band.

3. **Result analysis**:
- Independently predicted the X - ray luminosity function (XLF) and it is consistent with the observed XLF.
- Measured the BHAR of star - forming galaxies and quiescent galaxies and found that the BHAR of star - forming galaxies is usually greater than or at least comparable to that of quiescent galaxies.

### Key formulas:
- **Expression of BHAR**:
\[
\text{BHAR}(M_{\star}, z)=\int_{\log\lambda_{\text{min}}}^{+\infty}\frac{(1 - \epsilon)k_{\text{bol}}(M_{\star}\lambda)M_{\star}\lambda}{\epsilon c^{2}}p(\lambda|M_{\star}, z)d\log\lambda
\]
where \(\epsilon\) is the radiative efficiency of accretion, and \(k_{\text{bol}}(L_X)\) is the Boltzmann correction factor in the 2 - 10 keV energy band.

- **Detection probability function**:
\[
P_{\text{det}}(f_X)=\frac{1}{2}\left[\text{erf}(b(\log f_X - a)) + 1\right]
\]
where \(a\) and \(b\) are parameters that determine the shape of \(P_{\text{det}}(f_X)\).

- **Double - power - law model**:
\[
p(\lambda|M_{\star}, z)=
\begin{cases}
A\left(\frac{\lambda}{\lambda_c}\right)^{-\gamma_1},&\text{if }\lambda_{\text{min}}<\lambda\leq\lambda_c\\
A\left(\frac{\lambda}{\lambda_c}\right)^{-\gamma_2},&\text{if }\lambda > \lambda_c
\end{cases}
\]
where \(A,\lambda_c,\gamma_1,\gamma_2\) are parameters that depend on \(M_{\star}\) and \(z\).

### Conclusion:
This study provides the most accurate BHAR as a function of galaxies so far.