skip to main content


Title: Why do black holes trace bulges (& central surface densities), instead of galaxies as a whole?
ABSTRACT

Previous studies of fueling black holes in galactic nuclei have argued (on scales ${\sim}0.01{-}1000\,$pc) accretion is dynamical with inflow rates $\dot{M}\sim \eta \, M_{\rm gas}/t_{\rm dyn}$ in terms of gas mass Mgas, dynamical time tdyn, and some η. But these models generally neglected expulsion of gas by stellar feedback, or considered extremely high densities where expulsion is inefficient. Studies of star formation, however, have shown on sub-kpc scales the expulsion efficiency fwind = Mejected/Mtotal scales with the gravitational acceleration as $(1-f_{\rm wind})/f_{\rm wind}\sim \bar{a}_{\rm grav}/\langle \dot{p}/m_{\ast }\rangle \sim \Sigma _{\rm eff}/\Sigma _{\rm crit}$ where $\bar{a}_{\rm grav}\equiv G\, M_{\rm tot}(\lt r)/r^{2}$ and $\langle \dot{p}/m_{\ast }\rangle$ is the momentum injection rate from young stars. Adopting this as the simplest correction for stellar feedback, $\eta \rightarrow \eta \, (1-f_{\rm wind})$, we show this provides a more accurate description of simulations with stellar feedback at low densities. This has immediate consequences, predicting the slope and normalization of the MBH − σ and MBH − Mbulge relation, LAGN −SFR relations, and explanations for outliers in compact Es. Most strikingly, because star formation simulations show expulsion is efficient (fwind ∼ 1) below total-mass surface density $M_{\rm tot}/\pi \, r^{2}\lt \Sigma _{\rm crit}\sim 3\times 10^{9}\, \mathrm{M}_{\odot }\, {\rm kpc^{-2}}$ (where $\Sigma _{\rm crit}=\langle \dot{p}/m_{\ast }\rangle /(\pi \, G)$), BH mass is predicted to specifically trace host galaxy properties above a critical surface brightness Σcrit (B-band $\mu _{\rm B}^{\rm crit}\sim 19\, {\rm mag\, arcsec^{-2}}$). This naturally explains why BH masses preferentially reflect bulge properties or central surface densities (e.g. $\Sigma _{1\, {\rm kpc}}$), not ‘total’ galaxy properties.

 
more » « less
Award ID(s):
2009687 2108230 1715216 1652522 2108318 2001905
NSF-PAR ID:
10361033
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
510
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
p. 630-638
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT It has been established for decades that rotation curves deviate from the Newtonian gravity expectation given baryons alone below a characteristic acceleration scale $g_{\dagger }\sim 10^{-8}\, \rm {cm\, s^{-2}}$, a scale promoted to a new fundamental constant in MOND. In recent years, theoretical and observational studies have shown that the star formation efficiency (SFE) of dense gas scales with surface density, SFE ∼ Σ/Σcrit with $\Sigma _{\rm crit} \sim \langle \dot{p}/m_{\ast }\rangle /(\pi \, G)\sim 1000\, \rm {M_{\odot }\, pc^{-2}}$ (where $\langle \dot{p}/m_{\ast }\rangle$ is the momentum flux output by stellar feedback per unit stellar mass in a young stellar population). We argue that the SFE, more generally, should scale with the local gravitational acceleration, i.e. that SFE ${\sim}g_{\rm tot}/g_{\rm crit}\equiv (G\, M_{\rm tot}/R^{2}) / \langle \dot{p}/m_{\ast }\rangle$, where Mtot is the total gravitating mass and $g_{\rm crit}=\langle \dot{p}/m_{\ast }\rangle = \pi \, G\, \Sigma _{\rm crit} \approx 10^{-8}\, \rm {cm\, s^{-2}} \approx \mathit{ g}_{\dagger }$. Hence, the observed g† may correspond to the characteristic acceleration scale above which stellar feedback cannot prevent efficient star formation, and baryons will eventually come to dominate. We further show how this may give rise to the observed acceleration scaling $g_{\rm obs}\sim (g_{\rm baryon}\, g_{\dagger })^{1/2}$ (where gbaryon is the acceleration due to baryons alone) and flat rotation curves. The derived characteristic acceleration g† can be expressed in terms of fundamental constants (gravitational constant, proton mass, and Thomson cross-section): $g_{\dagger }\sim 0.1\, G\, m_{\mathrm{ p}}/\sigma _{\rm T}$. 
    more » « less
  2. ABSTRACT

    Galaxy sizes correlate closely with the sizes of their parent dark matter haloes, suggesting a link between halo formation and galaxy growth. However, the precise nature of this relation and its scatter remains to be understood fully, especially for low-mass galaxies. We analyse the galaxy–halo size relation (GHSR) for low-mass ($M_\star \sim 10^{7-9}\, {\rm M}_\odot$) central galaxies over the past 12.5 billion years with the help of cosmological volume simulations (FIREbox) from the Feedback in Realistic Environments (FIRE) project. We find a nearly linear relationship between the half-stellar mass galaxy size R1/2 and the parent dark matter halo virial radius Rvir. This relation evolves only weakly since redshift z = 5: $R_{1/2}\, [{\rm kpc}] = (0.053\pm 0.002)(R_{\rm vir}/35\, {\rm kpc})^{0.934\pm 0.054}$, with a nearly constant scatter $\langle \sigma \rangle = 0.084\, [{\rm dex}]$. While this ratio is similar to what is expected from models where galaxy disc sizes are set by halo angular momentum, the low-mass galaxies in our sample are not angular momentum supported, with stellar rotational to circular velocity ratios vrot/vcirc ∼ 0.15. Introducing redshift as another parameter to the GHSR does not decrease the scatter. Furthermore, this scatter does not correlate with any of the halo properties we investigate – including spin and concentration – suggesting that baryonic processes and feedback physics are instead critical in setting the scatter in the GHSR. Given the relatively small scatter and the weak dependence of the GHSR on redshift and halo properties for these low-mass central galaxies, we propose using galaxy sizes as an independent method from stellar masses to infer halo masses.

     
    more » « less
  3. ABSTRACT We revisit the question of ‘hot mode’ versus ‘cold mode’ accretion on to galaxies using steady-state cooling flow solutions and idealized 3D hydrodynamic simulations. We demonstrate that for the hot accretion mode to exist, the cooling time is required to be longer than the free-fall time near the radius where the gas is rotationally supported, Rcirc, i.e. the existence of the hot mode depends on physical conditions at the galaxy scale rather than on physical conditions at the halo scale. When allowing for the depletion of the halo baryon fraction relative to the cosmic mean, the longer cooling times imply that a virialized gaseous halo may form in halo masses below the threshold of $\sim 10^{12}\, {\rm M_{\odot }}$ derived for baryon-complete haloes. We show that for any halo mass there is a maximum accretion rate for which the gas is virialized throughout the halo and can accrete via the hot mode of ${\dot{M}}_{\rm crit}\approx 0.7(v_{\rm c}/100\, \rm km\ s^{-1})^{5.4}(R_{\rm circ}/10\, {\rm kpc})(Z/\, {\rm Z_{\odot }})^{-0.9}\, {\rm M_{\odot }}\, {\rm yr}^{-1}$, where Z and vc are the metallicity and circular velocity measured at Rcirc. For accretion rates $\gtrsim {\dot{M}}_{\rm crit}$ the volume-filling gas phase can in principle be ‘transonic’ – virialized in the outer halo but cool and free-falling near the galaxy. We compare ${\dot{M}}_{\rm crit}$ to the average star formation rate (SFR) in haloes at 0 < z < 10 implied by the stellar-mass–halo-mass relation. For a plausible metallicity evolution with redshift, we find that ${\rm SFR}\lesssim {\dot{M}}_{\rm crit}$ at most masses and redshifts, suggesting that the SFR of galaxies could be primarily sustained by the hot mode in halo masses well below the classic threshold of $\sim 10^{12}\, {\rm M_{\odot }}$. 
    more » « less
  4. ABSTRACT We present improved results of the measurement of the correlation between galaxies and the intergalactic medium transmission at the end of reionization. We have gathered a sample of 13 spectroscopically confirmed Lyman-break galaxies (LBGs) and 21 Lyman-α emitters (LAEs) at angular separations 20 arcsec ≲ θ ≲ 10 arcmin (∼0.1–4 pMpc at z ∼ 6) from the sightlines to eight background z ≳ 6 quasars. We report for the first time the detection of an excess of Lyman-α transmission spikes at ∼10–60 cMpc from LAEs (3.2σ) and LBGs (1.9σ). We interpret the data with an improved model of the galaxy–Lyman-α transmission and two-point cross-correlations, which includes the enhanced photoionization due to clustered faint sources, enhanced gas densities around the central bright objects and spatial variations of the mean free path. The observed LAE(LBG)–Lyman-α transmission spike two-point cross-correlation function (2PCCF) constrains the luminosity-averaged escape fraction of all galaxies contributing to reionization to $\langle f_{\rm esc} \rangle _{M_{\rm UV}\lt -12} = 0.14_{-0.05}^{+0.28}\, (0.23_{-0.12}^{+0.46})$. We investigate if the 2PCCF measurement can determine whether bright or faint galaxies are the dominant contributors to reionization. Our results show that a contribution from faint galaxies ($M_{\rm UV} \gt -20 \, (2\sigma)$) is necessary to reproduce the observed 2PCCF and that reionization might be driven by different sub-populations around LBGs and LAEs at z ∼ 6. 
    more » « less
  5. null (Ed.)
    ABSTRACT The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations, focusing on the MBH − M⋆ relation and the BH mass function. All simulations predict tight MBH − M⋆ relations, and struggle to produce BHs of $M_{\rm BH}\leqslant 10^{7.5}\, \rm M_{\odot }$ in galaxies of $M_{\star }\sim 10^{10.5}\!-\!10^{11.5}\, \rm M_{\odot }$. While the time evolution of the mean MBH − M⋆ relation is mild ($\rm \Delta M_{\rm BH}\leqslant 1\, dex$ for 0 $\leqslant z \leqslant$ 5) for all the simulations, its linearity (shape) and normalization varies from simulation to simulation. The strength of SN feedback has a large impact on the linearity and time evolution for $M_{\star }\leqslant 10^{10.5}\, \rm M_{\odot }$. We find that the low-mass end is a good discriminant of the simulation models, and highlights the need for new observational constraints. At the high-mass end, strong AGN feedback can suppress the time evolution of the relation normalization. Compared with observations of the local Universe, we find an excess of BHs with $M_{\rm BH}\geqslant 10^{9}\, \rm M_{\odot }$ in most of the simulations. The BH mass function is dominated by efficiently accreting BHs ($\log _{10}\, f_{\rm Edd}\geqslant -2$) at high redshifts, and transitions progressively from the high-mass to the low-mass end to be governed by inactive BHs. The transition time and the contribution of active BHs are different among the simulations, and can be used to evaluate models against observations. 
    more » « less