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 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.  more » « less
Award ID(s):
2009687
NSF-PAR ID:
10237287
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
503
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
3492 to 3509
Format(s):
Medium: X
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-mass 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. 
    more » « less
  2. ABSTRACT

    Radio-loud active galactic nuclei (RLAGNs) are a unique AGN population and were thought to be preferentially associated with supermassive black holes (SMBHs) at low accretion rates. They could impact the host galaxy evolution by expelling cold gas through the jet-mode feedback. In this work, we studied CO(6−5) line emission and continuum emission in a high-redshift radio galaxy, MRC 0152−209, at z = 1.92 using ALMA (Atacama Large Millimeter/submillimeter Array) up to a 0.024″ resolution (corresponding to ∼200 pc at z = 1.92). This system is a starburst major merger comprising two galaxies: the north-west (NW) galaxy hosting the RLAGN with jet kinetic power Ljet ≳ 2 × 1046  erg s−1 and the other galaxy to the south-east (SE). Based on the spectral energy distribution fitting for the entire system (NW+SE galaxies), we find an AGN bolometric luminosity LAGN, bol ∼ 3 × 1046  erg s−1 with a lower limit of ∼0.9 × 1046  erg s−1 for the RLAGN. We estimate the black hole mass through MBH–M⋆ scaling relations and find an Eddington ratio of λEdd ∼ 0.07–4 conservatively by adopting the lower limit of LAGN, bol and considering the dispersion of the scaling relation. These results suggest that the RLAGN is radiatively efficient and the powerful jets could be launched from a super-Eddington accretion disc. ALMA Cycle 6 observations further reveal a massive (${M}_\mathrm{H_2}=(1.1-2.3)\times 10^9\ \rm M_\odot$), compact (∼500 pc), and monopolar molecular outflow perpendicular to the jet axis. The corresponding mass outflow rate ($1200^{+300}_{-300}-2600^{+600}_{-600}\ \mathrm{M_\odot }\ \rm yr^{-1}$) is comparable with the star formation rate of at least $\sim 2100\ \mathrm{M_\odot }\ \rm yr^{-1}$. Depending on the outflowing molecular gas mass, the outflow kinetic power/LAGN, bol ratio of ∼0.008–0.02, and momentum boost factor of ∼3–24 agree with a radiative-mode AGN feedback scenario. On the other hand, the jets can also drive the molecular outflow within its lifetime of ∼2 × 105 yr without additional energy supply from AGN radiation. The jet-mode feedback is then capable of removing all cold gas from the host galaxy through the long-term, episodic launching of jets. Our study reveals a unique object where starburst activity, powerful jets, and rapid BH growth co-exist, which may represent a fundamental stage of AGN-host galaxy co-evolution.

     
    more » « less
  3. ABSTRACT

    We present Trinity, a flexible empirical model that self-consistently infers the statistical connection between dark matter haloes, galaxies, and supermassive black holes (SMBHs). Trinity is constrained by galaxy observables from 0 < z < 10 [galaxies’ stellar mass functions, specific and cosmic star formation rates (SFRs), quenched fractions, and UV luminosity functions] and SMBH observables from 0 < z < 6.5 (quasar luminosity functions, quasar probability distribution functions, active black hole mass functions, local SMBH mass–bulge mass relations, and the observed SMBH mass distributions of high-redshift bright quasars). The model includes full treatment of observational systematics [e.g. active galactic nucleus (AGN) obscuration and errors in stellar masses]. From these data, Trinity infers the average SMBH mass, SMBH accretion rate, merger rate, and Eddington ratio distribution as functions of halo mass, galaxy stellar mass, and redshift. Key findings include: (1) the normalization and the slope of the SMBH mass–bulge mass relation increases mildly from z = 0 to z = 10; (2) The best-fitting AGN radiative+kinetic efficiency is ∼0.05–0.06, but can be in the range ∼0.035–0.07 with alternative input assumptions; (3) AGNs show downsizing, i.e. the Eddington ratios of more massive SMBHs start to decrease earlier than those of lower mass objects; (4) The average ratio between average SMBH accretion rate and SFR is ∼10−3 for low-mass galaxies, which are primarily star-forming. This ratio increases to ∼10−1 for the most massive haloes below z ∼ 1, where star formation is quenched but SMBHs continue to accrete.

     
    more » « less
  4. ABSTRACT

    We study gas inflows on to supermassive black holes using hydrodynamics simulations of isolated galaxies and idealized galaxy mergers with an explicit, multiphase interstellar medium (ISM). Our simulations use the recently developed ISM and stellar evolution model called Stars and MUltiphase Gas in GaLaxiEs (SMUGGLE). We implement a novel super-Lagrangian refinement scheme that increases the gas mass resolution in the immediate neighbourhood of the black holes (BHs) to accurately resolve gas accretion. We do not include black hole feedback in our simulations. We find that the complex and turbulent nature of the SMUGGLE ISM leads to highly variable BH accretion. BH growth in SMUGGLE converges at gas mass resolutions ≲3 × 103 M⊙. We show that the low resolution simulations combined with the super-Lagrangian refinement scheme are able to produce central gas dynamics and BH accretion rates very similar to that of the uniform high resolution simulations. We further explore BH fueling by simulating galaxy mergers. The interaction between the galaxies causes an inflow of gas towards the galactic centres and results in elevated and bursty star formation. The peak gas densities near the BHs increase by orders of magnitude resulting in enhanced accretion. Our results support the idea that galaxy mergers can trigger AGN activity, although the instantaneous accretion rate depends strongly on the local ISM. We also show that the level of merger-induced enhancement of BH fueling predicted by the SMUGGLE model is much smaller compared to the predictions by simulations using an effective equation of state model of the ISM.

     
    more » « less
  5. null (Ed.)
    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 the 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. 
    more » « less