skip to main content

Title: The radio galaxy population in the simba simulations
ABSTRACT We examine the 1.4 GHz radio luminosities of galaxies arising from star formation and active galactic nuclei (AGNs) within the state-of-the-art cosmological hydrodynamic simulation Simba. Simba grows black holes via gravitational torque limited accretion from cold gas and Bondi accretion from hot gas, and employs AGN feedback including jets at low Eddington ratios. We define a population of radio loud AGNs (RLAGNs) based on the presence of ongoing jet feedback. Within RLAGN, we define high and low excitation radio galaxies (HERGs and LERGs) based on their dominant mode of black hole accretion: torque limited accretion representing feeding from a cold disc, or Bondi representing advection-dominated accretion from a hot medium. Simba predicts good agreement with the observed radio luminosity function (RLF) and its evolution, overall as well as separately for HERGs and LERGs. Quiescent galaxies with AGN-dominated radio flux dominate the RLF at $\gtrsim 10^{22-23}$ W Hz−1, while star formation dominates at lower radio powers. Overall, RLAGNs have higher black hole accretion rates and lower star formation rates than non-RLAGN at a given stellar mass or velocity dispersion, but have similar black hole masses. Simba predicts an LERG number density of 8.53 Mpc−3, ∼10× higher than for HERGs, broadly as observed. While LERGs dominate among more » most massive galaxies with the largest black holes and HERGs dominate at high specific star formation rates, they otherwise largely populate similar-sized dark matter haloes and have similar host galaxy properties. Simba thus predicts that deeper radio surveys will reveal an increasing overlap between the host galaxy demographics of HERGs and LERGs. « less
; ; ;
Award ID(s):
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
3492 to 3509
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT The sensitivity of X-ray facilities and our ability to detect fainter active galactic nuclei (AGNs) will increase with the upcoming Athena mission and the AXIS and Lynx concept missions, thus improving our understanding of supermassive black holes (BHs) in a luminosity regime that can be dominated by X-ray binaries. We analyse the population of faint AGNs ($L_{\rm x, 2{-}10 \, keV}\leqslant 10^{42}\, \rm erg\,s^{ -1}$) in the Illustris, TNG100, EAGLE, and SIMBA cosmological simulations, and find that the properties of their host galaxies vary from one simulation to another. In Illustris and EAGLE, faint AGNs are powered by low-massmore »BHs located in low-mass star-forming galaxies. In TNG100 and SIMBA, they are mostly associated with more massive BHs in quenched massive galaxies. We model the X-ray binary (XRB) populations of the simulated galaxies, and find that AGNs often dominate the galaxy AGN + XRB hard X-ray luminosity at z > 2, while XRBs dominate in some simulations at z < 2. Whether the AGN or XRB emission dominates in star-forming and quenched galaxies depends on the simulations. These differences in simulations can be used to discriminate between galaxy formation models with future high-resolution X-ray observations. We compare the luminosity of simulated faint AGN host galaxies to observations of stacked galaxies from Chandra. Our comparison indicates that the simulations post-processed with our X-ray modelling tend to overestimate the AGN + XRB X-ray luminosity; luminosity that can be strongly affected by AGN obscuration. Some simulations reveal clear AGN trends as a function of stellar mass (e.g. galaxy luminosity drop in massive galaxies), which are not apparent in the observations.« less
  2. ABSTRACT Supermassive black holes (SMBHs) that reside at the centres of galaxies can inject vast amounts of energy into the surrounding gas and are thought to be a viable mechanism to quench star formation in massive galaxies. Here, we study the $10^{9-12.5}\, \mathrm{M_\odot }$ stellar mass central galaxy population of the IllustrisTNG simulation, specifically the TNG100 and TNG300 volumes at z = 0, and show how the three components – SMBH, galaxy, and circumgalactic medium (CGM) – are interconnected in their evolution. We find that gas entropy is a sensitive diagnostic of feedback injection. In particular, we demonstrate how themore »onset of the low-accretion black hole (BH) feedback mode, realized in the IllustrisTNG model as a kinetic, BH-driven wind, leads not only to star formation quenching at stellar masses $\gtrsim 10^{10.5}\, \mathrm{M_\odot }$ but also to a change in thermodynamic properties of the (non-star-forming) gas, both within the galaxy and beyond. The IllustrisTNG kinetic feedback from SMBHs increases the average gas entropy, within the galaxy and in the CGM, lengthening typical gas cooling times from $10\!-\!100\, \mathrm{Myr}$ to $1\!-\!10\, \mathrm{Gyr}$, effectively ceasing ongoing star formation and inhibiting radiative cooling and future gas accretion. In practice, the same active galactic nucleus (AGN) feedback channel is simultaneously ‘ejective’ and ‘preventative’ and leaves an imprint on the temperature, density, entropy, and cooling times also in the outer reaches of the gas halo, up to distances of several hundred kiloparsecs. In the IllustrisTNG model, a long-lasting quenching state can occur for a heterogeneous CGM, whereby the hot and dilute CGM gas of quiescent galaxies contains regions of low-entropy gas with short cooling times.« less
  3. ABSTRACT Active galactic nuclei (AGNs) feedback is responsible for maintaining plasma in global thermal balance in extended haloes of elliptical galaxies and galaxy clusters. Local thermal instability in the hot gas leads to the formation of precipitating cold gas clouds that feed the central supermassive black holes, thus heating the hot gas and maintaining global thermal equilibrium. We perform 3D magnetohydrodynamical (MHD) simulations of self-regulated AGNs feedback in a Perseus-like galaxy cluster with the aim of understanding the impact of the feedback physics on the turbulence properties of the hot and cold phases of the intracluster medium (ICM). We findmore »that, in general, the cold phase velocity structure function (VSF) is steeper than the prediction from Kolmogorov’s theory. We attribute the physical origin of the steeper slope of the cold phase VSF to the driving of turbulent motions primarily by the gravitational acceleration acting on the ballistic clouds. We demonstrate that, in the pure hydrodynamical case, the precipitating cold filaments may be the dominant agent driving turbulence in the hot ICM. The arguments in favour of this hypothesis are that: (i) the cold phase mass dominates over hot gas mass in the inner cool core; (ii) hot and cold gas velocities are spatially correlated; (iii) both the cold and hot phase velocity distributions are radially biased. We show that, in the MHD case, the turbulence in the ambient hot medium (excluding the jet cone regions) can also be driven by the AGN jets. The driving is then facilitated by enhanced coupling due to magnetic fields of the ambient gas and the AGN jets. In the MHD case, turbulence may thus be driven by a combination of AGN jet stirring and filament motions. We conclude that future observations, including those from high spatial and spectral resolution X-ray missions, may help to constrain self-regulated AGN feedback by quantifying the multitemperature VSF in the ICM.« less
  4. Abstract Accreting black holes can drive fast and energetic nuclear winds that may be an important feedback mechanism associated with active galactic nuclei (AGN). In this paper, we implement a scheme for capturing feedback from these fast nuclear winds and examine their impact in simulations of isolated disk galaxies. Stellar feedback is modeled using the FIRE physics and produces a realistic multiphase interstellar medium (ISM). We find that AGN winds drive the formation of a low-density, high-temperature central gas cavity that is broadly consistent with analytic model expectations. The effects of AGN feedback on the host galaxy are a strongmore »function of the wind kinetic power and momentum. Low and moderate luminosity AGN do not have a significant effect on their host galaxy: the AGN winds inefficiently couple to the ambient ISM and instead a significant fraction of their energy vents in the polar direction. For such massive black holes, accretion near the Eddington limit can have a dramatic impact on the host galaxy ISM: if AGN wind feedback acts for ≳ 20 − 30 Myr, the inner ∼1 − 10 kpc of the ISM is disrupted and the global galaxy star formation rate is significantly reduced. We quantify the properties of the resulting galaxy-scale outflows and find that the radial momentum in the outflow is boosted by a factor ∼2 − 3 relative to that initially supplied in the AGN wind for strong feedback scenarios, decreasing below unity for less energetic winds. In contrast to observations, however, the outflows are primarily hot, with very little atomic or molecular gas. We conjecture that merging galaxies and high-redshift galaxies, which have more turbulent and thicker disks and very different nuclear gas geometries, may be even more disrupted by AGN winds than found in our simulations.« less
  5. Abstract We present the results of an analysis of Wide-field Infrared Survey Explorer (WISE) observations of the full 2500 deg 2 South Pole Telescope (SPT)-Sunyaev–Zel’dovich cluster sample. We describe a process for identifying active galactic nuclei (AGN) in brightest cluster galaxies (BCGs) based on WISE mid-IR color and redshift. Applying this technique to the BCGs of the SPT-SZ sample, we calculate the AGN-hosting BCG fraction, which is defined as the fraction of BCGs hosting bright central AGNs over all possible BCGs. Assuming an evolving single-burst stellar population model, we find statistically significant evidence (>99.9%) for a mid-IR excess at highmore »redshift compared to low redshift, suggesting that the fraction of AGN-hosting BCGs increases with redshift over the range of 0 < z < 1.3. The best-fit redshift trend of the AGN-hosting BCG fraction has the form (1 + z ) 4.1±1.0 . These results are consistent with previous studies in galaxy clusters as well as as in field galaxies. One way to explain this result is that member galaxies at high redshift tend to have more cold gas. While BCGs in nearby galaxy clusters grow mostly by dry mergers with cluster members, leading to no increase in AGN activity, BCGs at high redshift could primarily merge with gas-rich satellites, providing fuel for feeding AGNs. If this observed increase in AGN activity is linked to gas-rich mergers rather than ICM cooling, we would expect to see an increase in scatter in the P cav versus L cool relation at z > 1. Last, this work confirms that the runaway cooling phase, as predicted by the classical cooling-flow model, in the Phoenix cluster is extremely rare and most BCGs have low (relative to Eddington) black hole accretion rates.« less