skip to main content

Title: Accelerated Growth of Seed Black Holes by Dust in the Early Universe
Abstract

We explore the effect of dust on the growth of seed black holes (BHs) in the early universe. Previous 1D radiation-hydrodynamic (RHD) simulations show that increased radiation pressure on dust further suppresses the accretion rate than the case for the chemically pristine gas. Using the Enzo+Moray code, we perform a suite of 3D RHD simulations of accreting BHs in a dusty interstellar medium (ISM). We use the modified Grackle cooling library to consider dust physics in its nonequilibrium chemistry. The BH goes through an early evolutionary phase, where ionizing BH radiation creates an oscillating Hiiregion as it cycles between accretion and feedback. As the simulations proceed, dense cold gas accumulates outside the ionized region where inflow from the neutral medium meets the outflow driven by radiation pressure. In the late phase, high-density gas streams develop and break the quasi-spherical symmetry of the ionized region, rapidly boosting the accretion rate. The late phase is characterized by the coexistence of strong ionized outflows and fueling high-density gas inflows. The mean accretion rate increases with metallicity reaching a peak atZ∼ 0.01–0.1Z, one order of magnitude higher than the one for pristine gas. However, as the metallicity approaches the solar abundance, the mean more » accretion rate drops as the radiation pressure becomes strong enough to drive out the high-density gas. Our results indicate that a dusty metal-poor ISM can accelerate the growth rate of BHs in the early universe, but can also stun its growth as the ISM is further enriched toward the solar abundance.

« less
Authors:
; ;
Publication Date:
NSF-PAR ID:
10371097
Journal Name:
The Astrophysical Journal
Volume:
936
Issue:
2
Page Range or eLocation-ID:
Article No. 116
ISSN:
0004-637X
Publisher:
DOI PREFIX: 10.3847
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    Formation of supermassive black holes (BHs) remains a theoretical challenge. In many models, especially beginning from stellar relic ‘seeds,’ this requires sustained super-Eddington accretion. While studies have shown BHs can violate the Eddington limit on accretion disc scales given sufficient ‘fuelling’ from larger scales, what remains unclear is whether or not BHs can actually capture sufficient gas from their surrounding interstellar medium (ISM). We explore this in a suite of multiphysics high-resolution simulations of BH growth in magnetized, star-forming dense gas complexes including dynamical stellar feedback from radiation, stellar mass-loss, and supernovae, exploring populations of seeds with masses $\sim 1\!-\!10^{4}\, \mathrm{M}_{\odot }$. In this initial study, we neglect feedback from the BHs: so this sets a strong upper limit to the accretion rates seeds can sustain. We show that stellar feedback plays a key role. Complexes with gravitational pressure/surface density below $\sim 10^{3}\, \mathrm{M}_{\odot }\, {\rm pc^{-2}}$ are disrupted with low star formation efficiencies so provide poor environments for BH growth. But in denser cloud complexes, early stellar feedback does not rapidly destroy the clouds but does generate strong shocks and dense clumps, allowing $\sim 1{{\ \rm per\ cent}}$ of randomly initialized seeds to encounter a dense clump withmore »low relative velocity and produce runaway, hyper-Eddington accretion (growing by orders of magnitude). Remarkably, mass growth under these conditions is almost independent of initial BH mass, allowing rapid intermediate-mass black hole (IMBH) formation even for stellar-mass seeds. This defines a necessary (but perhaps not sufficient) set of criteria for runaway BH growth: we provide analytic estimates for the probability of runaway growth under different ISM conditions.

    « less
  2. ABSTRACT

    The early growth of black holes (BHs) in high-redshift galaxies is likely feedback regulated. While radiative feedback has been extensively studied, the role of mechanical feedback has received less scrutiny to date. Here, we use high-resolution parsec-scale hydrodynamical simulations to study jet propagation and its effect on 100 M⊙ BH accretion in the dense, low-metallicity gas expected in early protogalaxies. As the jet propagates, it shocks the surrounding gas forming a jet cocoon. The cocoon consists of a rapidly cooling cold phase at the interface with the background gas and an overpressured subsonic phase of reverse shock-heated gas filling the interior. We vary the background gas density and temperature, BH feedback efficiency, and the jet model. We found that the width of the jet cocoon roughly follows a scaling derived by assuming momentum conservation in the jet-propagation direction and energy conservation in the lateral directions. Depending on the assumed gas and jet properties, the cocoon either stays elongated to large radii or isotropizes before reaching the Bondi radius, forming a nearly spherical bubble. Lower jet velocities and higher background gas densities result in self-regulation to higher momentum fluxes and elongated cocoons. In all cases, the outward cocoon momentum flux balancesmore »the inward inflowing gas momentum flux near the Bondi radius, which ultimately regulates BH accretion. The time-averaged accretion rate always remains below the Bondi rate, and exceeds the Eddington rate only if the ambient medium is dense and cold, and/or the jet is weak (low velocity and mass loading).

    « less
  3. We study the long-term evolution of the global structure of axisymmetric accretion flows onto a black hole (BH) at rates substantially higher than the Eddington value (Mdot,Edd)performing two-dimensional hydrodynamical simulations with and without radiative diffusion. In the high-accretion optically-thick limit, where the radiation energy is efficiently trapped within the inflow, the accretion flow becomes adiabatic and comprises of turbulent gas in the equatorial region and strong bipolar outflows. As a result, the mass inflow rate decreases toward the center as Mdot,in∝r_p with p∼0.5−0.7 and a small fraction of the inflowing gas feeds the nuclear BH. Thus, super-Eddington accretion is sustained only when a larger amount of gas is supplied from larger radii at >100−1000 Mdot, Edd. The global structure of the flow settles down to a quasi-steady state in millions of the orbital timescale at the BH event horizon, which is >10−100 times longer than that addressed in previous (magneto-)RHD simulation studies. Energy transport via radiative diffusion accelerates the outflow near the poles in the inner region but does not change the overall properties of the accretion flow compared to the cases without diffusion. Based on our simulation results, we provide a mechanical feedback model for super-Eddington accreting BHs. Thismore »can be applied as a sub-grid model in large-scale cosmological simulations that do not sufficiently resolve galactic nuclei, and to the formation of the heaviest gravitational-wave sources via accretion in dense environments.« less
  4. Abstract The metallicity and gas density dependence of interstellar depletions, the dust-to-gas (D/G), and dust-to-metal (D/M) ratios have important implications for how accurately we can trace the chemical enrichment of the universe, either by using FIR dust emission as a tracer of the ISM or by using spectroscopy of damped Ly α systems to measure chemical abundances over a wide range of redshifts. We collect and compare large samples of depletion measurements in the Milky Way (MW), Large Magellanic Cloud (LMC) ( Z = 0.5 Z ⊙ ), and Small Magellanic Cloud (SMC) ( Z = 0.2 Z ⊙ ). The relations between the depletions of different elements do not strongly vary between the three galaxies, implying that abundance ratios should trace depletions accurately down to 20% solar metallicity. From the depletions, we derive D/G and D/M. The D/G increases with density, consistent with the more efficient accretion of gas-phase metals onto dust grains in the denser ISM. For log N (H) > 21 cm −2 , the depletion of metallicity tracers (S, Zn) exceeds −0.5 dex, even at 20% solar metallicity. The gas fraction of metals increases from the MW to the LMC (factor 3) and SMC (factor 6),more »compensating for the reduction in total heavy element abundances and resulting in those three galaxies having the same neutral gas-phase metallicities. The D/G derived from depletions are respective factors of 2 (LMC) and 5 (SMC) higher than the D/G derived from FIR, 21 cm, and CO emission, likely due to the combined uncertainties on the dust FIR opacity and on the depletion of carbon and oxygen.« less
  5. 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 themore »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.« less