No Arabic abstract
We show the Eddington ratio distributions of supermassive black holes at a wide redshift range (0 < z < 8) obtained with a semi-analytic model of galaxy formation. The distribution is broadly consistent with observational estimates at low redshift. We find that the growth rate of black holes at higher redshift is more likely to exceed the Eddington limit because the typical gas fraction of the host galaxies is higher at higher redshift. We also find that the super- Eddington growth is more common for less massive supermassive black holes, supporting an idea that supermassive black holes have been formed via super-Eddington accretion. These results indicate the slowing down of cosmic growth of supermassive black holes: the growth of supermassive black holes with a higher Eddington ratio peaks at higher redshift. We also show the effect of the sample selection on the shape of the Eddington ratio distribution functions and find that shallower observations will miss active galactic nuclei with not only the smaller but also higher Eddington ratios.
One of the main themes in extragalactic astronomy for the next decade will be the evolution of galaxies over cosmic time. Many future observatories, including JWST, ALMA, GMT, TMT and E-ELT will intensively observe starlight over a broad redshift range, out to the dawn of the modern Universe when the first galaxies formed. It has, however, become clear that the properties and evolution of galaxies are intimately linked to the growth of their central black holes. Understanding the formation of galaxies, and their subsequent evolution, will therefore be incomplete without similarly intensive observations of the accretion light from supermassive black holes (SMBH) in galactic nuclei. To make further progress, we need to chart the formation of typical SMBH at z>6, and their subsequent growth over cosmic time, which is most effectively achieved with X-ray observations. Recent technological developments in X-ray optics and instrumentation now bring this within our grasp, enabling capabilities fully matched to those expected from flagship observatories at longer wavelengths.
We investigate the rapid growth phase of supermassive black holes (BHs) within the hydrodynamical cosmological eagle simulation. This non-linear phase of BH growth occurs within $sim$$L_{*}$ galaxies, embedded between two regulatory states of the galaxy host: in sub $L_{*}$ galaxies efficient stellar feedback regulates the gas inflow onto the galaxy and significantly reduces the growth of the central BH, while in galaxies more massive than $L_{*}$ efficient AGN feedback regulates the gas inflow onto the galaxy and curbs further non-linear BH growth. We find evolving critical galaxy and halo mass scales at which rapid BH growth begins. Galaxies in the low-redshift Universe transition into the rapid BH growth phase in haloes that are approximately an order of magnitude more massive than their high-redshift counterparts (M{200} $approx 10^{12.4}$~Msol at $z approx 0$ decreasing to M{200} $approx 10^{11.2}$~Msol at $z approx 6$). Instead, BHs enter the rapid growth phase at a fixed critical halo virial temperature ($T_{mathrm{vir}} approx 10^{5.6}$~K). We additionally show that major galaxy--galaxy interactions ($mu geq frac{1}{4}$, where $mu$ is the stellar mass ratio) play a substantial role in triggering the rapid growth phase of BHs in the low-redshift Universe, whilst potentially having a lower influence at high redshift. Approximately 40% of BHs that initiate the rapid BH growth phase at $z approx 0$ do so within $pm 0.5$ dynamical times of a major galaxy--galaxy merger, a fourfold increase above what is expected from the background merger rate. We find that minor mergers ($frac{1}{10} leq mu < frac{1}{4}$) have a substantially lower influence in triggering the rapid growth phase at all epochs.
The next generation of electromagnetic and gravitational wave observatories will open unprecedented windows to the birth of the first supermassive black holes. This has the potential to reveal their origin and growth in the first billion years, as well as the signatures of their formation history in the local Universe. With this in mind, we outline three key focus areas which will shape research in the next decade and beyond: (1) What were the seeds of the first quasars; how did some reach a billion solar masses before z$sim7$? (2) How does black hole growth change over cosmic time, and how did the early growth of black holes shape their host galaxies? What can we learn from intermediate mass black holes (IMBHs) and dwarf galaxies today? (3) Can we unravel the physics of black hole accretion, understanding both inflows and outflows (jets and winds) in the context of the theory of general relativity? Is it valid to use these insights to scale between stellar and supermassive BHs, i.e., is black hole accretion really scale invariant? In the following, we identify opportunities for the Canadian astronomical community to play a leading role in addressing these issues, in particular by leveraging our strong involvement in the Event Horizon Telescope, the {it James Webb Space Telescope} (JWST), Euclid, the Maunakea Spectroscopic Explorer (MSE), the Thirty Meter Telescope (TMT), the Square Kilometer Array (SKA), the Cosmological Advanced Survey Telescope for Optical and ultraviolet Research (CASTOR), and more. We also discuss synergies with future space-based gravitational wave (LISA) and X-ray (e.g., Athena, Lynx) observatories, as well as the necessity for collaboration with the stellar and galactic evolution communities to build a complete picture of the birth of supermassive black holes, and their growth and their influence over the history of the Universe.
We use global three dimensional radiation magneto-hydrodynamical simulations to study accretion disks onto a $5times 10^8M_{odot}$ black hole with accretion rates varying from $sim 250L_{Edd}/c^2$ to $1500 L_{Edd}/c^2$. We form the disks with torus centered at $50-80$ gravitational radii with self-consistent turbulence initially generated by the magneto-rotational instability. We study cases with and without net vertical magnetic flux. The inner regions of all disks have radiation pressure $sim 10^4-10^6$ times the gas pressure. Non-axisymmetric density waves that steepen into spiral shocks form as gas flows towards the black hole. In simulations without net vertical magnetic flux, Reynolds stress generated by the spiral shocks are the dominant mechanism to transfer angular momentum. Maxwell stress from MRI turbulence can be larger than the Reynolds stress only when net vertical magnetic flux is sufficiently large. Outflows are formed with speed $sim 0.1-0.4c$. When the accretion rate is smaller than $sim 500 L_{Edd}/c^2$, outflows start around $10$ gravitational radii and the radiative efficiency is $sim 5%-7%$ with both magnetic field configurations. With accretion rate reaching $1500 L_{Edd}/c^2$, most of the funnel region close to the rotation axis becomes optically thick and the outflow only develops beyond $50$ gravitational radii. The radiative efficiency is reduced to $1%$. We always find the kinetic energy luminosity associated with the outflow is only $sim 15%-30%$ of the radiative luminosity. The mass flux lost in the outflow is $sim 15%-50%$ of the net mass accretion rates. We discuss implications of our simulation results on the observational properties of these disks.
Electron-positron pair creation near sub-Eddington accretion rate black holes is believed to be dominated by the Breit-Wheeler process (photon-photon collisions). The interacting high energy photons are produced when unscreened electric fields accelerate leptons either in coherent, macroscopic gaps or in incoherent structures embedded in the turbulent plasma flow. The latter type of acceleration results in a drizzle of pair production sourced by photons from the background radiation field whose energies are near the pair-production threshold. In this work, we use radiation GRMHD simulations to extend an earlier study of pair drizzle by Moscibrodzka et al. We focus on low-magnetization (SANE) accretion onto supermassive Kerr black holes and consider radiation due to synchrotron, bremsstrahlung, and Compton upscattering processes. We confirm that pair drizzle in M87 is sufficient to keep the magnetospheric charge density orders of magnitude above the Goldreich-Julian density. We also find that pair production peaks along the jet-disk boundary.