Next Article in Journal
History of Astronomy in Australia: Big-Impact Astronomy from World War II until the Lunar Landing (1945–1969)
Next Article in Special Issue
Optical Characteristics of Electromagnetic Radiation, Emitted by Particles or Stars Moving near Supermassive Black Hole
Previous Article in Journal
Neutrino Charge in a Magnetized Media
Previous Article in Special Issue
Leptonic and Hadronic Radiative Processes in Supermassive-Black-Hole Jets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Big and Young Supermassive Black Holes in the Early Universe

INAF—Osservatorio Astronomico di Brera, via E. Bianchi 46, I-23807 Merate, LC, Italy
Submission received: 29 January 2021 / Revised: 24 March 2021 / Accepted: 25 March 2021 / Published: 30 March 2021
(This article belongs to the Special Issue Jet Physics of Accreting Super Massive Black Holes)

Abstract

:
Blazars are active galactic nuclei characterized by relativistic jets launched in the vicinity of the central engine (i.e., a supermassive black hole), which are oriented close to our line of sight. Their peculiar orientation makes them very efficient tracers of the overall jetted population, and due to their brightness they can be visible up to very high redshifts. A deep knowledge of these objects can provide fundamental clues to the models of formation and growth of the first supermassive black holes, but the search for them in the early Universe must be careful and follow a systematic approach. The discovery in the last ∼15 years of extremely massive blazars at very high redshifts ( M BH > 10 9 M ,   z > 4 ) revolutionized our perception of their earliest evolution: there seem to be different formation epochs for extremely massive black holes hosted in jetted ( z 4 ) and non-jetted ( z 2.5 ) systems. This is not easy to explain, since one would expect jetted sources to accrete less efficiently. Small differences in the population are also derived from the search for such high-z sources. We will go through the open questions in order to understand where the common knowledge stands and which steps must be taken to better understand the formation and common evolution of supermassive black holes and jets in the early Universe.

1. Introduction

Active galactic nuclei (AGNs) are accreting supermassive black holes located in the center of massive galaxies. They are characterized by strong persistent emission covering the whole electromagnetic spectrum. In particular, optical-UV radiation is produced by the release of gravitational energy from the matter accreting towards the central black hole.
About 10% of the overall AGN population show collimated relativistic jets launched from the vicinity of the central engine [1]. When reaching the intra-cluster medium, these structures are slowed down and their release of energy creates extended lobes filled with strong magnetic fields, which allow them to emit isotropically in the radio band through synchrotron cooling. For this reason, an AGN showing strong radio emission is generally thought to be hosting a jet. A quantitative indication of jet presence was introduced by Fanaroff and Riley [2] in their effort to classify the position of bright radio regions of extragalactic sources: the radio-loudness R = F r / F opt measures the dominance of the radio over optical monochromatic emissions (i.e., the presence of a jet brighter than the accretion emission). Being emitted from extended lobes, radio emission is visible independently of the jet orientation with respect to an observer’s line of sight. Collimated jets instead emit at all wavelengths, and are strongly boosted due to relativistic beaming along their propagation direction. Their emission can thus dominate the whole AGN spectrum when aligned along an observer’s line of sight, in which case the sources are classified as blazars. The collimated radio emission from the jet is easily much brighter than the isotropic lobe radio emission. The radio-loudness will thus reach extreme values of few thousands in the case of a well-aligned jet and will be a first indicator of alignment.
Blazars’ broad-band emission is characterized by two humps: synchrotron emission is responsible for the one peaking around sub-millimeter to soft X-ray frequencies [3], while the inverse Compton on electrons or hadrons produces the high-energy peak that falls between X-rays and TeV γ -rays. This peculiar emission profile makes blazars easily detectable at very high frequencies. In fact, quasars detected in the γ -rays are generally blazars: AGN all-sky catalogs such as the one compiled by Fermi/Large Area Telescope (LAT) are mainly composed of such sources (see Section 2.1).

2. Finding High-Redshift Blazars

The peculiar orientation of blazars makes them perfect tracers for the whole family of jetted AGNs. For each blazar observed and studied, one can assume the presence of many analogous jetted AGNs with their jets oriented in other directions with the same features (black hole mass, accretion rate, jet power, …). Knowing the viewing angle ( θ v ) and the jet beaming angle ( θ b ) of a blazar allows us to estimate the number of jetted quasars composing its parent population ( N parent ). If an object is observed under a viewing angle smaller than the beaming angle, in fact, one can estimate the number of objects that have the same properties but with jets pointing elsewhere.
We can assume that the emitting region is relativistically moving with an overall Lorentz factor Γ = 1 / 1 β 2 , and emitting isotropically in its own rest frame. According to an outside observer the region appears to emit half of its photons into a cone with semi–aperture θ b , where sin θ b = 1 / Γ [4]. The emitting region is highly relativistic, and therefore moves with β = v / c 1 . This implies a very small value for θ b , and therefore:
sin θ b θ b 1 Γ
We therefore define the beaming angle of a relativistic jet as θ b = 1 / Γ .
We can then derive the ratio between the number of objects observed with θ v 1 / Γ and the total number of analogous jetted AGNs, but randomly oriented. This estimate corresponds to the ratio f Ω between the solid angle subtended by the jet beaming angle Ω b and the all-sky surface Ω 4 π . Note that the jetted emission from a single AGN is formed by a jet (that we actually see in the case of a blazar) and a counterjet. Therefore, we must take into account both jet and counterjet in the calculation, and we will be able to express f Ω only as a function of the Lorentz factor Γ :
f Ω = 2 Ω b Ω 4 π = 1 β 1 β 2 2 = 1 2 Γ 2
That corresponds to the ratio between the number of objects seen at θ v 1 / Γ and the total number of jetted objects with jets oriented randomly, i.e., N parent = 1 / f Ω . Therefore, the number of analogous jetted AGNs oriented in the other directions is
N parent 2 Γ 2 = 338 Γ 13 2
where Γ = 13 is a typical value for a blazar [5]. At low redshift this is not particularly more efficient than compiling a complete catalog of misaligned AGNs to explore the parent population, since extensive studies and complete quasar catalogs are actually available. At high redshift, instead, a complete collection of jetted AGNs is not yet available, and therefore blazars can be the main actors in the systematic study of the jetted AGN population, thanks to their statistical relevancy.
In the 2000s, a few blazars were discovered serendipitously at z > 4 by observing their high-energy components. The most notable example is Q0906+6930 [6,7], the only high-redshift blazar ever detected by the Energetic Gamma–Ray Experiment Telescope (EGRET) and for 16 years the most distant blazar known, being located at z = 5.47 . Its blazar nature was further confirmed with very-long-baseline interferometry (VLBI), Very Large Array (VLA) and Chandra observations. The rest of the z > 4 blazars serendipitously discovered between 2000 and 2010 are consistently very radio-loud, and their blazar nature could be confirmed through their X-ray emission [8,9,10]. No z > 4 source other than Q0906+6930 has been serendipitously discovered in the GeV frequency range.

2.1. All-Sky Surveys

The systematic discovery of blazars is guaranteed at low redshift by high-energy instruments like Fermi/LAT. This instrument provides the largest γ –ray catalog of AGNs (the Third Catalog of AGNs detected by Fermi, 3LAC [11]), that mostly include blazars, along with very few radio galaxies and unclassified AGNs. 3LAC includes 414 Flat Spectrum Radio Quasars (FSRQs) and 604 BL Lacertae objects (BL Lacs), but they all have redshift z 3.1 .
This cut is not due only to the sensitivity limit, however, but the high redshift and the blazar sequence are mostly responsible for it. At high redshift, in fact, we expect to see only the most luminous objects, and the most luminous blazars have spectral energy distribution (SED) peaks at smaller frequencies [12,13]. This effect, along with the high redshift, shifts the high-energy hump of the SED below 100 MeV. The detection of this component by the Fermi/LAT therefore becomes unlikely. Instead, a powerful blazar can be more easily detected at frequencies just before the inverse Compton peak by hard X-ray instruments, like the Burst Alert Telescope (BAT) onboard the Swift satellite [14]. This is true already at z 3 . Unfortunately the Swift/BAT has a very high sensitivity limit, and therefore is extremely limited to objects up to z 4 . At redshift higher than 4, a systematic search at high energies is therefore impossible to perform. A different approach is needed.
There are only five blazars known at z > 3.1 detected by Fermi/LAT, which have been actively searched for in the all-sky data of the instrument [15]. The discovery of their γ -ray emission follows a systematic and specific approach, starting from known high-redshift radio-loud quasars and digging into ∼92 months of Fermi Pass 8 source class photons. A single source from this sample is located at z > 4 : NVSS J151002+570243 ( z = 4.3 ) is currently the farthest known γ -ray-emitting blazar, but it could not be discovered simply by the all-sky Fermi/LAT survey, because of the limited sensitivity.

2.2. Candidate Selection

Serendipitous findings of high-redshift blazars are clearly not efficient for obtaining a complete census of these sources. For this reason, a systematic approach has been followed in the last years. As already pointed out, very bright blazars have been observed to have their signature humps shifted towards lower frequencies [12,13]. This effect, described by the so-called blazar sequence, only affects the jet emission. All the thermal radiation originated in the accretion system is thus left visible in the rest-frame UV-optical frequency range. This feature has allowed different research groups to select reliable blazar candidates. A systematic approach is possible if a complete quasar catalog is used as a starting point for searching for blazar candidates. This is feasible since only FSRQs can be observed at such high redshift: BL Lacs (if present) would not have reliable redshift estimates. Specifically in the optical range, therefore, we can be sure that the thermal emission from the accretion disc dominates the continuum emission, along with the emission lines from broad- and narrow-line regions that are visible in that wavelength range.
In order to select good high-redshift blazar candidates, one should limit the search to quasars with z > 4 if starting from an optical spectroscopic quasar catalog. However, it is possible to select from scratch new high-z quasars, as was done in order to discover the farthest blazar currently known, PSO J0309+27 [16]. These authors selected bright and compact i-dropout sources from the Pan-STARSS survey, with a relatively strong radio-counterpart. Their selection from the near-IR-optical survey is the standard procedure for finding new z > 5.7 quasars, followed by various groups focusing on the matter. Their selection in radio, instead, mirrors the search for extreme radio-loudness; at z > 4 the radio flux is bright enough to allow cuts towards extreme radio-loudness values (namely >100), while at higher redshift a strong radio flux is already a reliable indication of the presence of a jet.
This more systematic approach has currently led to the discovery and classification of 15 blazars at z > 4 , from well defined, flux-limited samples, that allow for the derivation of strong constraints on the overall jetted population (see the following Sections).

3. Discussion

The discovery of the high-redshift blazar population is still an open process, but the samples currently available can already tell us something about the massive black holes accreting in the early Universe. All results need to be considered as a lower limit on the population, since the complete census is still in the process of being completed.

3.1. How Many Jetted AGNs at High-z?

We have already demonstrated in detail how blazars are a statistically significant representation of the jetted population, because of their orientation. This peculiarity can be used to compare the jetted population with their non-jetted counterparts. We focus on the most massive sources, i.e., M BH > 10 9 M , to explore the evolution of the most extreme supermassive black holes across the cosmic time. Moreover, most quasars observed at very high-redshift with a measured mass host such a massive black hole. Henceforth, we can draw interesting conclusions on the AGN population that we are able to observe in the early Universe.
Blazars selected from complete optical quasar surveys are generally limited to a specific area of the sky. Those systematically selected at z > 4 from the Sloan Digital Sky Survey (SDSS) and Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) surveys, in particular, are limited to the solid angle Ω SDSS + FIRST 8770 deg 2 . It is therefore reasonable to assume that the amount of blazars we find in that specific sky region is a fraction Ω SDSS + FIRST / 4 π of the total number of blazars with those features distributed in the whole sky. Knowing the size of the selection area of the blazar candidates allows their discovery to be more statistically relevant than serendipitously discovered sources. The observation of one of them traces the presence of N parent jetted AGNs, where
N parent 2 Γ 2 4 π Ω SDSS + FIRST = 2 Γ 2 41253 deg 2 8770 deg 2 .
It is therefore possible to calculate how many jetted AGNs with same black hole mass and accretion rate can be inferred. This can thus be used efficiently to study the distribution of early extremely massive black holes, providing constraints for models of supermassive black hole formation and evolution.
Figure 1 shows the comoving number density of extremely massive black holes (i.e., with M BH > 10 9 M ) hosted by radio-quiet (blue line, derived from [17]) and radio-loud AGNs, i.e., non-jetted and jetted sources. On the other hand, the comoving number density of radio-loud quasars is derived from Fermi/LAT [18] and Swift/BAT [19] blazar luminosity functions, assuming that every blazar traces the existence of 2 Γ 2 jetted quasars [20]. The blazar surveys of Swift/BAT and Fermi/LAT are limited at z < 4 , thus at z > 4 there is no complete survey of blazars that provides a reliable number density. For z > 4 the density was previously assumed to decrease exponentially, as per the corresponding one for radio-quiet objects. Nevertheless, a hint of different density distributions between jetted and non-jetted objects was already visible in [20,21]. While radio-quiet objects seem to be more numerous at redshift z 2–2.5, the distribution of jetted quasars suggests that they do not show such a peak, but more likely a plateau that extends at least up to z 4 , or maybe a second peak at z 4 . Those papers did not reach a firm conclusion about the real distribution of jetted AGNs at z > 4 , lacking a significant number of objects in that redshift frame. In other words, they could not disentangle between a plateau and a high-redshift peak.
The number densities derived by classifying systematically selected blazar candidates are instead very significant. Current observations in fact clearly push towards an interesting conclusion: the density of extremely massive black holes hosted in jetted systems peaks at least around z 4 , while the non-jetted systems peak at z 2–2.5. In other words, the number of heavy black holes hosted in jetted systems increases with cosmic time, and peaks at z > 4 . This suggests two different epochs of supermassive black hole formation. The most striking feature of this result is that black holes that grow while developing a jet seem to be born earlier, or to grow faster, than black holes with same masses that do not show jets.

3.2. Early Quasar Evolution

The presence of a jet in AGNs is commonly linked to high values of black hole spin. The black hole rotational energy, in fact, is a possible source of power to launch the jet [23,24], while the angular momentum of the black hole itself easily identifies its launching direction. This does not facilitate fast accretion, according to common knowledge. Rapidly spinning black holes, in fact, are generally associated to high accretion-based radiative efficiencies, since the innermost stable orbit of the accretion disc approaches the gravitational radius. Thorne [25] showed that the maximum spin that a black hole can reach through accretion is a = 0.998 (instead of a = 1 ), since the hole slows down by capturing accretion disc radiation, which exerts a counteracting torque. The accretion efficiency that corresponds to this spin value is η = 0.3 , larger than the efficiency in the case of a non-rotating black hole ( η = 0.08 ). A larger efficiency corresponds to a lower overall accretion rate (since M · = L d / η c 2 ). For a fixed luminosity, then, a black hole accretes less matter and grows slower. A spinning black hole accreting at the Eddington rate (i.e., at the highest rate allowed for a standard disc), if we assume that it doubles its mass on a Salpeter [26] time, would need 3.1 Gyr to grow from a seed of 100 M to 10 9 M [21]. Hence, such massive black holes would not be visible at z > 2.1 , while we see from Figure 1 that their preferential formation epoch seems to be around z 4 .
Some options for a faster accretion in the presence of a jet are explored in the literature [21,27]. A fascinating toy model suggests how the jet could facilitate a fast accretion, contrary to what is generally thought. The available energy, in fact, is not all radiated away, but it may contribute in amplifying the magnetic field and thus launching the jet. We can define two different contributions to the global efficiency η : one related to the disc luminosity (i.e., the radiative efficiency η d ), and one implied in generating the magnetic field that helps in launching the jet η j :
η = η d + η j .
In this case, to produce a fixed luminosity, we need a larger accretion rate, since part of the gravitational energy is spent for some non radiative processes necessary to form the jet. In other words, a spinning black hole can accrete faster if some of its accretion efficiency is spent for the jet ( η j > 0 and η = 0.3 ). Considering this, the accretion is faster, but black holes with M BH > 10 9 M are still hard to form before z 4–5.
The blazars found at z 5 (∼1Gyr) all have M BH 10 9 M . If the black hole is maximally spinning ( η = 0.3 ), its mere presence requires η j 0.2 and η d 0.1 , as can be seen in Figure 2, if we assume a black hole seed of 100 M . Assuming that a more massive seed helps, we are still forced to consider a maximally spinning black hole that does not convert its entire gravitational energy into radiation, but also uses this energy in generating a magnetic field to launch the jet. This is reasonable and (probably) necessary for high-redshift objects, but is not required to justify the presence of heavy black holes in jetted systems at low redshift. At z 2 , in fact, more time is available to form such massive black holes. Note also that a fast accretion requires a large amount of cold gas in the surroundings of the black hole, to continuously fuel the accretion. This is a reasonable picture at high redshift, but the black hole environment at lower redshift is generally poorer. This means that other less extreme and more plausible processes can drive supermassive black hole formation and growth. This is the reason why the presence of such massive objects is not a major problem, while their presence at z > 4 is harder to justify.
At this stage, the toy model here described is only a plausible idea to justify such a relevant population of heavy black holes in jetted systems, and therefore the presence of these sources at high redshift still remains an interesting open issue. Another straightforward option is that these supermassive black holes are accreting in a super-Eddington regime, but there is no observation-based indication that this could be the case for the currently known sources.

4. Open Issues

4.1. Missing Jetted Sources

Figure 1 clearly shows another peculiar behavior of high-redshift jetted sources with respect to non-jetted ones and their lower redshift counterparts. Not only do radio-loud sources seem to be differently distributed, but at z 4.5 they seem almost to outnumber the radio-quiet ones, or at least to have comparable number densities. This is strikingly different from the ratio of the two populations at low redshift: it is calculated that jetted AGNs are ∼10% of the overall AGN population, certainly not half, or its entirety.
A discrepancy in the population ratio is already evident when considering only the jetted population. The systematic search for blazar candidates, as in [28], leads to a blazar-to-radio-loud ratio that is very different from the lower redshift one; at least 8 blazars (3 other sources will be classified as blazars in an upcoming paper, making the blazar-to-radio ratio even more extreme) have already been classified out of the sample of the SDSS and FIRST spectroscopic quasar catalog, where there were only 53 radio-loud sources. This translates to a ratio of ∼15% blazars in a sample of radio-loud quasars at z 4 , clearly much higher than what was expected from the alignment-based calculations, i.e., less than 1 % [29]. Might misaligned jetted sources be missing, or not visible?
Clearly, the discrepancy could be due to different jet features at very high redshift, which would twist the ratio between aligned and misaligned jetted sources from the standard 2 Γ 2 . This is rather hard to justify, though, since the broad-band SED profile of z > 4 blazars does not show major differences when compared to lower redshift ones. A known effect is the quenching of extended radio emissions in favor of up-scattering of CMB photons by the lobe electron population happening at z > 3 [30]. The CMB energy density at such high redshifts becomes comparable with the magnetic energy density in lobe regions of jetted AGNs, offering an alternative and preferential interaction to electrons that would normally emit by synchrotron cooling. This effect is widely confirmed by the observation of extended X-ray emissions visible along jet directions in high-z quasars [31] and easily justifies the loss of widely misaligned blazar counterparts at z > 4 .
Numbers suggest instead that even the mildly misaligned sources are missing. This discrepancy can be fixed by taking into consideration the larger IR luminosity observed in high-z quasars, which has been observed in z > 7 quasars and partially obscured z > 4 ones [32,33]; high-redshift AGNs might be over-obscured by “bubbles” of dust, possibly typical in the early Universe [29]. Such bubbles would be able to completely obscure the nuclear emission, up to a threshold luminosity that would eventually eject the dust because of its strong radiation pressure. According to this scenario, only the jet would be able to pierce through the bubble, opening a small aperture that would allow direct observation of the thermal emission only when the jet was closely aligned to the line-of-sight. This would prevent slightly misaligned jetted AGNs from being visible as standard quasars in optical-UV photometric and spectroscopic catalogs during the over-obscured phase.
These features reconcile the results in the frame of jetted sources, but do not help in understanding why radio-loud (jetted) quasars are as or more numerous than radio-quiet (non-jetted) quasars. Might the presence of a jet actually help in fast accretion?

4.2. Where Does the Jet Point?

Many sources classified as blazars on the basis of X-ray observations from the z 4 SDSS and FIRST samples were also observed with Very Large Array (VLA) or very-long-baseline interferometry (VLBI) high-resolution configurations. While in most cases the two approaches yield consistent results [34], some of these sources do not show the radio features expected from blazar-like AGNs [35]. As an example, some show extended radio-emitting regions that suggest larger viewing angles, or a strong core-dominance is lacking as in the case of PSO J0309+27 [36].
Different possibilities could be invoked to handle this interesting inconsistency. It is possible that low- and high-redshift jets in AGNs emit with different processes, which still have to be completely explored. A more simple explanation relies on a different geometry of the jet: X-ray and radio emissions are emitted at different distances from the central engine, and while at low-z the relativistic jet is clearly able to propagate along the same direction for the whole time, at very high-z this might change. Different densities of the medium where the jet propagates might make its path difficult, and ultimately bend the jet. This would explain the different viewing angles derived from X-ray and radio emissions.
This option clearly needs to be further explored, and it will constitute an interesting field to be investigated in the upcoming years, with more sensitive facilities and more detailed modeling of the early formation of supermassive black holes and relativistic jets. A tight collaboration between X-ray, radio and theoretical astronomers will be crucial, being clearly necessary for a completely holistic approach to solving this puzzle.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AGNActive galactic nucleus
BL LacBL Lacertae object
Fermi/LATThe Large Area Telescope onboard the Fermi Gamma-Ray Space Telescope
FIRSTFaint Images of the Radio Sky at Twenty-Centimeters
FRSQFlat spectrum radio quasars
SDSSSloan Digital Sky Survey
Swift/BATBurst Alert Telescope onboard the Niel Gehrels Swift Observatory
VLAVery Large Array
VLBIVery-long-baseline interferometry

References

  1. Urry, C.M.; Padovani, P. Unified Schemes for Radio-Loud Active Galactic Nuclei. Publ. Astron. Soc. Pac. 1995, 107, 803. [Google Scholar] [CrossRef] [Green Version]
  2. Fanaroff, B.L.; Riley, J.M. The morphology of extragalactic radio sources of high and low luminosity. Mon. Not. R. Astron. Soc. 1974, 167, 31P–36P. [Google Scholar] [CrossRef] [Green Version]
  3. Padovani, P.; Alexander, D.M.; Assef, R.J.; De Marco, B.; Giommi, P.; Hickox, R.C.; Richards, G.T.; Smolčić, V.; Hatziminaoglou, E.; Mainieri, V.; et al. Active galactic nuclei: What’s in a name? Astron. Astrophys. Rev. 2017, 25, 2. [Google Scholar] [CrossRef] [Green Version]
  4. Rybicki, G.B.; Lightman, A.P. Radiative Processes in Astrophysics; John Wiley & Sons: Hoboken, NJ, USA, 1985. [Google Scholar]
  5. Ghisellini, G.; Tavecchio, F. Canonical high-power blazars. Mon. Not. R. Astron. Soc. 2009, 397, 985–1002. [Google Scholar] [CrossRef]
  6. Romani, R.W.; Sowards-Emmerd, D.; Greenhill, L.; Michelson, P. Q0906+6930: The Highest Redshift Blazar. Astrophys. J. Lett. 2004, 610, L9–L11. [Google Scholar] [CrossRef] [Green Version]
  7. Romani, R.W. The Spectral Energy Distribution of the High-z Blazar Q0906+6930. Astron. J. 2006, 132, 1959–1963. [Google Scholar] [CrossRef] [Green Version]
  8. Bassett, L.C.; Brandt, W.N.; Schneider, D.P.; Vignali, C.; Chartas, G.; Garmire, G.P. Chandra Observations of Radio-Loud Quasars at z>4: X-rays from the Radio Beacons of the Early Universe. Astron. J. 2004, 128, 523–533. [Google Scholar] [CrossRef]
  9. Shemmer, O.; Brandt, W.N.; Vignali, C.; Schneider, D.P.; Fan, X.; Richards, G.T.; Strauss, M.A. The X-ray Spectral Properties and Variability of Luminous High-Redshift Active Galactic Nuclei. Astrophys. J. 2005, 630, 729–739. [Google Scholar] [CrossRef]
  10. Young, M.; Elvis, M.; Risaliti, G. The Fifth Data Release Sloan Digital Sky Survey/XMM-Newton Quasar Survey. Astrophys. J. Suppl. 2009, 183, 17–32. [Google Scholar] [CrossRef] [Green Version]
  11. Ackermann, M.; Ajello, M.; Atwood, W.B.; Baldini, L.; Ballet, J.; Barbiellini, G.; Bastieri, D.; Becerra Gonzalez, J.; Bellazzini, R.; Bissaldi, E.; et al. The Third Catalog of Active Galactic Nuclei Detected by the Fermi Large Area Telescope. Astrophys. J. 2015, 810, 14. [Google Scholar] [CrossRef] [Green Version]
  12. Fossati, G.; Maraschi, L.; Celotti, A.; Comastri, A.; Ghisellini, G. A unifying view of the spectral energy distributions of blazars. Mon. Not. R. Astron. Soc. 1998, 299, 433–448. [Google Scholar] [CrossRef]
  13. Ghisellini, G.; Righi, C.; Costamante, L.; Tavecchio, F. The Fermi blazar sequence. Mon. Not. R. Astron. Soc. 2017, 469, 255–266. [Google Scholar] [CrossRef]
  14. Gehrels, N.; Chincarini, G.; Giommi, P.; Mason, K.O.; Nousek, J.A.; Wells, A.A.; White, N.E.; Barthelmy, S.D.; Burrows, D.N.; Cominsky, L.R.; et al. The Swift Gamma-Ray Burst Mission. Astrophys. J. 2004, 611, 1005–1020. [Google Scholar] [CrossRef] [Green Version]
  15. Ackermann, M.; Ajello, M.; Baldini, L.; Ballet, J.; Barbiellini, G.; Bastieri, D.; Becerra Gonzalez, J.; Bellazzini, R.; Bissaldi, E.; Blandford, R.D.; et al. Gamma-Ray Blazars within the First 2 Billion Years. Astrophys. J. Lett. 2017, 837, L5. [Google Scholar] [CrossRef]
  16. Belladitta, S.; Moretti, A.; Caccianiga, A.; Spingola, C.; Severgnini, P.; Della Ceca, R.; Ghisellini, G.; Dallacasa, D.; Sbarrato, T.; Cicone, C.; et al. The first blazar observed at z > 6. Astron. Astrophys. 2020, 635, L7. [Google Scholar] [CrossRef] [Green Version]
  17. Shen, X.; Hopkins, P.F.; Faucher-Giguère, C.A.; Alexander, D.M.; Richards, G.T.; Ross, N.P.; Hickox, R.C. The bolometric quasar luminosity function at z = 0–7. Mon. Not. R. Astron. Soc. 2020, 495, 3252–3275. [Google Scholar] [CrossRef]
  18. Ajello, M.; Shaw, M.S.; Romani, R.W.; Dermer, C.D.; Costamante, L.; King, O.G.; Max-Moerbeck, W.; Readhead, A.; Reimer, A.; Richards, J.L.; et al. The Luminosity Function of Fermi-detected Flat-spectrum Radio Quasars. Astrophys. J. 2012, 751, 108. [Google Scholar] [CrossRef] [Green Version]
  19. Ajello, M.; Costamante, L.; Sambruna, R.M.; Gehrels, N.; Chiang, J.; Rau, A.; Escala, A.; Greiner, J.; Tueller, J.; Wall, J.V.; et al. The Evolution of Swift/BAT Blazars and the Origin of the MeV Background. Astrophys. J. 2009, 699, 603–625. [Google Scholar] [CrossRef] [Green Version]
  20. Ghisellini, G.; Della Ceca, R.; Volonteri, M.; Ghirlanda, G.; Tavecchio, F.; Foschini, L.; Tagliaferri, G.; Haardt, F.; Pareschi, G.; Grindlay, J. Chasing the heaviest black holes of jetted active galactic nuclei. Mon. Not. R. Astron. Soc. 2010, 405, 387–400. [Google Scholar] [CrossRef] [Green Version]
  21. Ghisellini, G.; Haardt, F.; Della Ceca, R.; Volonteri, M.; Sbarrato, T. The role of relativistic jets in the heaviest and most active supermassive black holes at high redshift. Mon. Not. R. Astron. Soc. 2013, 432, 2818–2823. [Google Scholar] [CrossRef]
  22. Hopkins, P.F.; Richards, G.T.; Hernquist, L. An Observational Determination of the Bolometric Quasar Luminosity Function. Astrophys. J. 2007, 654, 731–753. [Google Scholar] [CrossRef]
  23. Blandford, R.D.; Znajek, R.L. Electromagnetic extraction of energy from Kerr black holes. Mon. Not. R. Astron. Soc. 1977, 179, 433–456. [Google Scholar] [CrossRef]
  24. Tchekhovskoy, A.; Narayan, R.; McKinney, J.C. Efficient generation of jets from magnetically arrested accretion on a rapidly spinning black hole. Mon. Not. R. Astron. Soc. 2011, 418, L79–L83. [Google Scholar] [CrossRef]
  25. Thorne, K.S. Disk-Accretion onto a Black Hole. II. Evolution of the Hole. Astrophys. J. 1974, 191, 507–520. [Google Scholar] [CrossRef]
  26. Salpeter, E.E. Accretion of Interstellar Matter by Massive Objects. Astrophys. J. 1964, 140, 796–800. [Google Scholar] [CrossRef]
  27. Mazzucchelli, C.; Bañados, E.; Venemans, B.P.; Decarli, R.; Farina, E.P.; Walter, F.; Eilers, A.C.; Rix, H.W.; Simcoe, R.; Stern, D.; et al. Physical Properties of 15 Quasars at z ≳ 6.5. Astrophys. J. 2017, 849, 91. [Google Scholar] [CrossRef] [Green Version]
  28. Sbarrato, T.; Ghisellini, G.; Nardini, M.; Tagliaferri, G.; Greiner, J.; Rau, A.; Schady, P. Blazar candidates beyond redshift 4 observed with GROND. Mon. Not. R. Astron. Soc. 2013, 433, 2182–2193. [Google Scholar] [CrossRef] [Green Version]
  29. Ghisellini, G.; Sbarrato, T. Dark bubbles around high-redshift radio-loud active galactic nucleus. Mon. Not. R. Astron. Soc. 2016, 461, L21–L25. [Google Scholar] [CrossRef] [Green Version]
  30. Ghisellini, G.; Celotti, A.; Tavecchio, F.; Haardt, F.; Sbarrato, T. Radio-loud active galactic nuclei at high redshifts and the cosmic microwave background. Mon. Not. R. Astron. Soc. 2014, 438, 2694–2700. [Google Scholar] [CrossRef] [Green Version]
  31. Yuan, W.; Fabian, A.C.; Celotti, A.; Jonker, P.G. Extended X-ray emission in the high-redshift quasar GB 1508+5714 at z = 4.3. Mon. Not. R. Astron. Soc. 2003, 346, L7–L10. [Google Scholar] [CrossRef]
  32. Mortlock, D.J.; Warren, S.J.; Venemans, B.P.; Patel, M.; Hewett, P.C.; McMahon, R.G.; Simpson, C.; Theuns, T.; Gonzáles-Solares, E.A.; Adamson, A.; et al. A luminous quasar at a redshift of z = 7.085. Nature 2011, 474, 616–619. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Fan, L.; Chen, W.; An, T.; Xie, F.G.; Han, Y.; Knudsen, K.K.; Yang, J. The Hyperluminous, Dust-obscured Quasar W2246-0526 at z = 4.6: Detection of Parsec-scale Radio Activity. Astrophys. J. Lett. 2020, 905, L32. [Google Scholar] [CrossRef]
  34. Frey, S.; Paragi, Z.; Fogasy, J.O.; Gurvits, L.I. The first estimate of radio jet proper motion at z > 5. Mon. Not. R. Astron. Soc. 2015, 446, 2921–2928. [Google Scholar] [CrossRef]
  35. Cao, H.M.; Frey, S.; Gabányi, K.É.; Paragi, Z.; Yang, J.; Cseh, D.; Hong, X.Y.; An, T. VLBI observations of four radio quasars at z > 4: Blazars or not? Mon. Not. R. Astron. Soc. 2017, 467, 950–960. [Google Scholar] [CrossRef] [Green Version]
  36. Spingola, C.; Dallacasa, D.; Belladitta, S.; Caccianiga, A.; Giroletti, M.; Moretti, A.; Orienti, M. Parsec-scale properties of the radio brightest jetted AGN at z > 6. Astron. Astrophys. 2020, 643, L12. [Google Scholar] [CrossRef]
Figure 1. Comoving number density of active supermassive black holes with M BH 10 9 M hosted in jetted (orange line) and non-jetted active galactic nuclei (AGNs) (blue line). The non-jetted density is derived from the luminosity function of radio-quiet quasars in the Sloan Digital Sky Survey (SDSS) [17], by assuming that all quasars with L opt > 10 47 erg / cm 2 / s host black holes with masses we are interested in (following the assumptions considered by Ghisellini et al. [20] while extracting the same quantities from Hopkins et al. [22]). As a comparison, the light blue dashed line is the comoving number density derived from Hopkins et al. [22]. The jetted-related curve is instead formed by the sum of luminosity functions from Fermi/Large Area Telescope (LAT) (yellow line and points) and Swift/Burst Alert Telescope (BAT) (red line and points), up to z 4 . Green hexagonal lower limits are derived from the systematic classification of blazar candidates from SDSS and Faint Images of the Radio Sky at Twenty-Centimeters (FIRST), based on X-ray observations and corrected by the sky area covered by the survey. The yellow diamond lower limit refers to PSO J0309+27, the farthest blazar currently known.
Figure 1. Comoving number density of active supermassive black holes with M BH 10 9 M hosted in jetted (orange line) and non-jetted active galactic nuclei (AGNs) (blue line). The non-jetted density is derived from the luminosity function of radio-quiet quasars in the Sloan Digital Sky Survey (SDSS) [17], by assuming that all quasars with L opt > 10 47 erg / cm 2 / s host black holes with masses we are interested in (following the assumptions considered by Ghisellini et al. [20] while extracting the same quantities from Hopkins et al. [22]). As a comparison, the light blue dashed line is the comoving number density derived from Hopkins et al. [22]. The jetted-related curve is instead formed by the sum of luminosity functions from Fermi/Large Area Telescope (LAT) (yellow line and points) and Swift/Burst Alert Telescope (BAT) (red line and points), up to z 4 . Green hexagonal lower limits are derived from the systematic classification of blazar candidates from SDSS and Faint Images of the Radio Sky at Twenty-Centimeters (FIRST), based on X-ray observations and corrected by the sky area covered by the survey. The yellow diamond lower limit refers to PSO J0309+27, the farthest blazar currently known.
Galaxies 09 00023 g001
Figure 2. Black hole mass evolution as a function of the age of the Universe. Red points are the known quasars with virial-based mass measures at z > 5.7 . The hexagons show the mass range of the highest-z blazar known [16]. The grey contours show the whole SDSS DR7 quasar sample, including the few tens of radio-loud sources (blue points). The solid lines show the Eddington-limited evolution based on the Salpeter timescale, starting from two different black hole seeds, and assuming that (i) the energy released during the accretion η coincides with the accretion disc radiative efficiency η d in the case of a Shakura–Sunyaev disk around a non-spinning black hole (yellow line), (ii) same but in the case of a maximally spinning black hole (green line), (iii) the black hole is maximally spinning and thus the energy release is η = 0.3 , but only a fraction of it is radiated through the accretion disc (light blue line).
Figure 2. Black hole mass evolution as a function of the age of the Universe. Red points are the known quasars with virial-based mass measures at z > 5.7 . The hexagons show the mass range of the highest-z blazar known [16]. The grey contours show the whole SDSS DR7 quasar sample, including the few tens of radio-loud sources (blue points). The solid lines show the Eddington-limited evolution based on the Salpeter timescale, starting from two different black hole seeds, and assuming that (i) the energy released during the accretion η coincides with the accretion disc radiative efficiency η d in the case of a Shakura–Sunyaev disk around a non-spinning black hole (yellow line), (ii) same but in the case of a maximally spinning black hole (green line), (iii) the black hole is maximally spinning and thus the energy release is η = 0.3 , but only a fraction of it is radiated through the accretion disc (light blue line).
Galaxies 09 00023 g002
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sbarrato, T. Big and Young Supermassive Black Holes in the Early Universe. Galaxies 2021, 9, 23. https://0-doi-org.brum.beds.ac.uk/10.3390/galaxies9020023

AMA Style

Sbarrato T. Big and Young Supermassive Black Holes in the Early Universe. Galaxies. 2021; 9(2):23. https://0-doi-org.brum.beds.ac.uk/10.3390/galaxies9020023

Chicago/Turabian Style

Sbarrato, Tullia. 2021. "Big and Young Supermassive Black Holes in the Early Universe" Galaxies 9, no. 2: 23. https://0-doi-org.brum.beds.ac.uk/10.3390/galaxies9020023

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop