Supermassive stars with masses of M* ≃ 104–105 M⊙ are invoked as possible seeds of high-redshift supermassive black holes, but it remains under debate whether their protostar indeed acquires sufficient mass via gas accretion overcoming radiative feedback. We investigate protostellar growth in dynamically heated atomic cooling haloes (ACHs) found in recent cosmological simulations, performing three-dimensional radiation hydrodynamical simulations that consider stellar evolution under variable mass accretion. We find that one of the ACHs feeds the central protostar at rates exceeding a critical value, above which the star evolves in a cool bloating phase and hardly produces ionizing photons. Consequently, the stellar mass reaches M* ≳ 104 M⊙ unimpeded by radiative feedback. In the other ACH, where the mass supply rate is lower, the star evolves almost as a hot main-sequence star, emitting intense ionizing radiation. Then, the stellar mass growth is terminated around 500 M⊙ by photoevaporation of the circumstellar disc. Our simulations provide a formula of the final stellar mass determined either by stellar feedback or their lifetime as a function of the mass supply rate from the parent cloud. Combining the results with the statistical properties of star-forming clouds in high-redshift quasar progenitor haloes, we construct a top-heavy mass distribution of primordial stars over M* ≃ 100–105 M⊙, approximately following a power-law spectrum of ${\propto} M_\ast ^{-1.3}$. Their black hole remnants would be further fed via the dense debris disc, powering ‘milliquasars’ with a bolometric luminosity of Lbol ≳ 1043 erg s−1.
- Award ID(s):
- 1743747
- NSF-PAR ID:
- 10188516
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 489
- Issue:
- 1
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- 802 to 819
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
ABSTRACT -
ABSTRACT The physical origin of the seeds of supermassive black holes (SMBHs), with postulated initial masses ranging from ∼105 M⊙ to as low as ∼102 M⊙, is currently unknown. Most existing cosmological hydrodynamic simulations adopt very simple, ad hoc prescriptions for BH seeding and seed at unphysically high masses ∼105–106 M⊙. In this work, we introduce a novel sub-grid BH seeding model for cosmological simulations that is directly calibrated to high-resolution zoom simulations that explicitly resolve ∼103 M⊙ seeds forming within haloes with pristine, dense gas. We trace the BH growth along galaxy merger trees until their descendants reach masses of ∼104 or 105 M⊙. The results are used to build a new stochastic seeding model that directly seeds these descendants in lower resolution versions of our zoom region. Remarkably, we find that by seeding the descendants simply based on total galaxy mass, redshift and an environmental richness parameter, we can reproduce the results of the detailed gas-based seeding model. The baryonic properties of the host galaxies are well reproduced by the mass-based seeding criterion. The redshift-dependence of the mass-based criterion captures the combined influence of halo growth, dense gas formation, and metal enrichment on the formation of ∼103 M⊙ seeds. The environment-based seeding criterion seeds the descendants in rich environments with higher numbers of neighbouring galaxies. This accounts for the impact of unresolved merger dominated growth of BHs, which produces faster growth of descendants in richer environments with more extensive BH merger history. Our new seed model will be useful for representing a variety of low-mass seeding channels within next-generation larger volume uniform cosmological simulations.
-
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.
-
null (Ed.)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
-
Supermassive stars (SMSs) with masses of 𝑀∗ ≃ 104–105 M⊙ are invoked as possible seeds of high-redshift supermassive black holes, but it remains under debate whether their protostar indeed acquires sufficient mass via gas accretion overcoming radiative feedback. We investigate protostellar growth in dynamically heated atomic-cooling haloes (ACHs) found in recent cosmological simulations, performing three-dimensional radiation hydrodynamical (RHD) simulations that consider stellar evolution under variable mass accretion. We find that one of the ACHs feeds the central protostar at rates exceeding a critical value, above which the star evolves in a cool bloating phase and hardly produces ionizing photons. Consequently, the stellar mass reaches 𝑀∗ 104 M⊙ unimpeded by radiative feedback. In the other ACH, where the mass supply rate is lower, the star spends most of its life as a hot main-sequence star, emitting intense ionizing radiation. Then, the stellar mass growth is terminated around 500 M⊙ by photoevaporation of the circumstellar disk. A series of our RHD simulations provide a formula of the final stellar mass determined either by stellar feedback or their lifetime as a function of the mass supply rate from the parent cloud in the absence of stellar radiation. Combining the results with the statistical properties of SMS-forming clouds in high-redshift quasar progenitor haloes, we construct a top-heavy mass distribution of primordial stars over 𝑀∗ ≃ 100–105 M⊙, approximately following a power-law spectrum of ∝ 𝑀−1.3 with a steeper decline at 𝑀 2 × 104 M . Their massive BH remnants would be ∗∗⊙ further fed via the dense debris disk, powering “milli-quasars" with a bolometric luminosity of 𝐿bol 1043 erg s−1.more » « less