skip to main content


Title: Probing the z ≳ 6 quasars in a universe with IllustrisTNG physics: impact of gas-based black hole seeding models
ABSTRACT

We explore implications of a range of black hole (BH) seeding prescriptions on the formation of the brightest $z$ ≳ 6 quasars in cosmological hydrodynamic simulations. The underlying galaxy formation model is the same as in the IllustrisTNG simulations. Using constrained initial conditions, we study the growth of BHs in rare overdense regions (forming $\gtrsim 10^{12}\, {\rm M}_{\odot }\,h^{-1}$ haloes by $z$ = 7) using a  (9 Mpc h−1)3 simulated volume. BH growth is maximal within haloes that are compact and have a low tidal field. For these haloes, we consider an array of gas-based seeding prescriptions wherein $M_{\mathrm{seed}}=10^4\!-\!10^6\, {\rm M}_{\odot }\,h^{-1}$ seeds are inserted in haloes above critical thresholds for halo mass and dense, metal-poor gas mass (defined as $\tilde{M}_{\mathrm{h}}$ and $\tilde{M}_{\mathrm{sf,mp}}$, respectively, in units of Mseed). We find that a seed model with $\tilde{M}_{\mathrm{sf,mp}}=5$ and $\tilde{M}_{\mathrm{h}}=3000$ successfully produces a $z$ ∼ 6 quasar with $\sim 10^9\, {\rm M}_{\odot }$ mass and ∼1047 erg s−1 luminosity. BH mergers play a crucial role at $z$ ≳ 9, causing an early boost in BH mass at a time when accretion-driven BH growth is negligible. With more stringent seeding conditions (e.g. $\tilde{M}_{\mathrm{sf,mp}}=1000$), the relative paucity of BH seeds results in a much lower merger rate. In this case, $z$ ≳ 6 quasars can only be formed if we enhance the maximum allowed BH accretion rates (by factors ≳10) compared to the accretion model used in IllustrisTNG. This can be achieved either by allowing for super-Eddington accretion, or by reducing the radiative efficiency. Our results demonstrate that progenitors of $z$ ∼ 6 quasars have distinct BH merger histories for different seeding models, which will be distinguishable with Laser Interferometer Space Antenna observations.

 
more » « less
Award ID(s):
1909933 2008490
NSF-PAR ID:
10370035
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
516
Issue:
1
ISSN:
0035-8711
Page Range / eLocation ID:
p. 138-157
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    Direct collapse black holes (BHs) are promising candidates for producing massive z ≳ 6 quasars, but their formation requires fine-tuned conditions. In this work, we use cosmological zoom simulations to study systematically the impact of requiring: (1) low gas angular momentum (spin), and (2) a minimum incident Lyman–Werner (LW) flux in order to form BH seeds. We probe the formation of seeds (with initial masses of $M_{\rm seed} \sim 10^4\!-\!10^6\, \mathrm{M}_{\odot }\, h^{-1})$ in haloes with a total mass >3000 × Mseed and a dense, metal-poor gas mass >5 × Mseed. Within this framework, we find that the seed-forming haloes have a prior history of star formation and metal enrichment, but they also contain pockets of dense, metal-poor gas. When seeding is further restricted to haloes with low gas spins, the number of seeds formed is suppressed by factors of ∼6 compared to the baseline model, regardless of the seed mass. Seed formation is much more strongly impacted if the dense, metal-poor gas is required to have a critical LW flux (Jcrit). Even for Jcrit values as low as 50J21, no $8\times 10^{5}~\mathrm{M}_{\odot }\, h^{-1}$ seeds are formed. While lower mass ($1.25\times 10^{4},1\times 10^{5}~\mathrm{M}_{\odot }\, h^{-1}$) seeds do form, they are strongly suppressed (by factors of ∼10–100) compared to the baseline model at gas mass resolutions of $\sim 10^4~\mathrm{M}_{\odot }\, h^{-1}$ (with even stronger suppression at higher resolutions). As a result, BH merger rates are also similarly suppressed. Since early BH growth is dominated by mergers in our models, none of the seeds are able to grow to the supermassive regime ($\gtrsim 10^6~\mathrm{M}_{\odot }\, h^{-1}$) by z = 7. Our results hint that producing the bulk of the z ≳ 6 supermassive BH population may require alternate seeding scenarios that do not depend on the LW flux, early BH growth dominated by rapid or super-Eddington accretion, or a combination of these possibilities.

     
    more » « less
  2. null (Ed.)
    Abstract Deciphering the formation of supermassive black holes (SMBHs) is a key science goal for upcoming observational facilities. In many theoretical channels proposed so far, the seed formation depends crucially on local gas conditions. We systematically characterize the impact of a range of gas-based black hole seeding prescriptions on SMBH populations using cosmological simulations. Seeds of mass Mseed ∼ 103–106 M⊙ h−1 are placed in haloes that exceed critical thresholds for star-forming, metal-poor gas mass and halo mass (defined as $\tilde{M}_{\mathrm{sf,mp}}$ and $\tilde{M}_{\mathrm{h}}$, respectively, in units of Mseed). We quantify the impact of these parameters on the properties of z ≥ 7 SMBHs. Lower seed masses produce higher black hole merger rates (by factors of ∼10 and ∼1000 at z ∼ 7 and z ∼ 15, respectively). For fixed seed mass, we find that $\tilde{M}_{\mathrm{h}}$ has the strongest impact on the black hole population at high redshift (z ≳ 15, where a factor of 10 increase in $\tilde{M}_{\mathrm{h}}$ suppresses merger rates by ≳ 100). At lower redshift (z ≲ 15), we find that $\tilde{M}_{\mathrm{sf,mp}}$ has a larger impact on the black hole population. Increasing $\tilde{M}_{\mathrm{sf,mp}}$ from 5–150 suppresses the merger rates by factors of ∼8 at z ∼ 7–15. This suggests that the seeding criteria explored here could leave distinct imprints on LISA merger rates. In contrast, AGN luminosity functions are much less sensitive to seeding criteria, varying by factors ≲ 2 − 3 within our models. Such variations will be challenging to probe even with future sensitive instruments such as Lynx or JWST. Our study provides a useful benchmark for development of seed models for large-volume cosmological simulations. 
    more » « less
  3. ABSTRACT

    We investigate the black hole mass function at z ∼ 5 using XQz5, our recent sample of the most luminous quasars between the redshifts 4.5 < z < 5.3. We include 72 quasars with black hole masses estimated from velocity-broadened emission-line measurements and single-epoch virial prescriptions in the footprint of a highly complete parent survey. The sample mean Eddington ratio and standard deviation is log λ ≈ −0.20 ± 0.24. The completeness-corrected mass function is modelled as a double power law, and we constrain its evolution across redshift assuming accretion-dominated mass growth. We estimate the evolution of the mass function from z = 5–4, presenting joint constraints on accretion properties through a measured dimensionless e-folding parameter, kef ≡ 〈λ〉U(1 − ϵ)/ϵ = 1.79 ± 0.06, where 〈λ〉 is the mean Eddington ratio, U is the duty cycle, and ϵ is the radiative efficiency. If these supermassive black holes were to form from seeds smaller than $10^8\, {\rm M}_{\odot }$, the growth rate must have been considerably faster at z ≫ 5 than observed from z = 5–4. A growth rate exceeding 3 × the observed rate would reduce the initial heavy seed mass to $10^{5-6}\, {\rm M}_{\odot }$, aligning with supermassive star and/or direct collapse seed masses. Stellar mass ($10^2\, {\rm M}_{\odot }$) black hole seeds would require ≳4.5 × the observed growth rate at z ≫ 5 to reproduce the measured active black hole mass function. A possible pathway to produce the most extreme quasars is radiatively inefficient accretion flow, suggesting black holes with low angular momentum or photon trapping in supercritically accreting thick discs.

     
    more » « less
  4. 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
  5. ABSTRACT Accretion discs around supermassive black holes are promising sites for stellar mass black hole mergers detectable with LIGO. Here we present the results of Monte Carlo simulations of black hole mergers within 1-d AGN disc models. For the spin distribution in the disc bulk, key findings are: (1) The distribution of χeff is naturally centred around $\tilde{\chi }_{\rm eff} \approx 0.0$, (2) the width of the χeff distribution is narrow for low natal spins. For the mass distribution in the disc bulk, key findings are: (3) mass ratios $\tilde{q} \sim 0.5\!-\!0.7$, (4) the maximum merger mass in the bulk is $\sim 100\!-\!200\, \mathrm{M}_{\odot }$, (5) $\sim 1{{\ \rm per\ cent}}$ of bulk mergers involve BH $\gt 50\, \mathrm{M}_{\odot }$ with (6) $\simeq 80{{\ \rm per\ cent}}$ of bulk mergers are pairs of first generation BH. Additionally, mergers at a migration trap grow an IMBH with typical merger mass ratios $\tilde{q}\sim 0.1$. Ongoing LIGO non-detections of black holes $\gt 10^{2}\, \mathrm{M}_{\odot }$ puts strong limits on the presence of migration traps in AGN discs (and therefore AGN disc density and structure) as well as median AGN disc lifetime. The highest merger rate occurs for this channel if AGN discs are relatively short-lived (≤1 Myr) so multiple AGN episodes can happen per Galactic nucleus in a Hubble time. 
    more » « less