Luminosity Function Of Active Galactic Nuclei Research Articles

What Are the Pillars of Reionization? Revising the AGN Luminosity Function at z ∼ 5

In the past, high-z active galactic nuclei (AGNs) were given a minor role as possible drivers of reionization, despite initial evidence in favor of their large space densities at low luminosities by Chandra and the Hubble Space Telescope. Recent observations from JWST are finding relatively large numbers of faint AGNs at z > 4, convincingly confirming these early results. We present a sample of z ∼ 5 AGNs, both from wide, shallow ground-based surveys and from deep, pencil-beam observations from JWST, allowing us to estimate their space densities with unprecedented accuracy. The bright end (M 1450 < −26) of the z ∼ 5 AGN luminosity function is well constrained, with a rather steep slope. The faint end (M 1450 ≥ −22) indicates a high space density, the scatter is significant, and the knee (M 1450 ∼ −24) is mostly undetermined. Comparisons with state-of-the-art models find reasonable agreement with the observed AGN luminosity function at z = 5, while the predicted space density evolution at higher redshifts appears to be too fast with respect to observational constraints. Given the large variance at the faint end, we consider different options in fitting the luminosity functions and deriving the ionizing emissivity. Even in the most conservative scenario, the photoionization rate produced by z ∼ 5 AGNs is consistent with the ultraviolet background measurements. A slow evolution of the space density of faint AGNs is observed, indicating that active SMBHs are probably producing large amounts of ionizing photons at z > 6, well into the Epoch of Reionization. This is an important indication that high-z AGNs could be major contributors to the reionization of the Universe.

MUSE view of PDS 456: Kiloparsec-scale wind, extended ionized gas, and close environment

PDS 456 is the most luminous (Lbol ∼ 1047 erg s−1) radio-quiet quasar at z &lt; 0.3 and can be regarded as a local counterpart of the powerful quasars shining at Cosmic Noon. It hosts a strong nuclear X-ray ultra-fast (∼0.3c) outflow, and a massive and clumpy CO (3–2) molecular outflow extending up to ∼5 kpc from the nucleus. We analyzed the first MUSE Wide Field Mode (WFM) and Adaptive-Optics Narrow Field Mode (AO-NFM) optical integral field spectroscopic observations of PDS456. The AO-NFM observations provide an unprecedented spatial resolution, reaching up to ∼280 pc. Our findings reveal a complex circumgalactic medium around PDS 456, extending to a maximum projected size of ≈46 kpc. This includes a reservoir of gas with a mass of ∼107 − 108 M⊙, along with eight companion galaxies and a multi-phase outflow. WFM and NFM MUSE data reveal an outflow on a large scale (≈12 kpc from the quasar) in [O III], and on smaller scales (within 3 kpc) with higher resolution (about 280 pc) in Hα, respectively. The [O III] outflow mass rate is 2.3 ± 0.2 M⊙ yr−1 which is significantly lower than those typically found in other luminous quasars. Remarkably, the Hα outflow shows a similar scale, morphology, and kinematics to the CO (3–2) molecular outflow, with the latter dominating in terms of kinetic energy and mass outflow rate by two and one orders of magnitude, respectively. Our results therefore indicate that mergers, powerful active galactic nucleus (AGN) activity, and feedback through AGN-driven winds collectively contribute to shaping the host galaxy evolution of PDS 456, and likely that of similar objects at the brightest end of the AGN luminosity function across all redshifts. Moreover, the finding that the momentum boost of the total outflow deviates from the expected energy-conserving expansion for large-scale outflows highlights the need of novel AGN-driven outflow models to comprehensively interpret these phenomena.

A Candidate Supermassive Black Hole in a Gravitationally Lensed Galaxy at Z ≈ 10

While supermassive black holes (SMBHs) are widely observed in the nearby and distant Universe, their origin remains debated with two viable formation scenarios with light and heavy seeds. In the light seeding model, the seed of the first SMBHs form from the collapse of massive stars with masses of 10–100 M ⊙, while the heavy seeding model posits the formation of 104–5 M ⊙ seeds from direct collapse. The detection of SMBHs at redshifts z ≳ 10, edging closer to their formation epoch, provides critical observational discrimination between these scenarios. Here, we focus on the JWST-detected galaxy, GHZ 9, at z ≈ 10 that is lensed by the foreground cluster, A2744. Based on 2.1 Ms deep Chandra observations, we detect a candidate X-ray active galactic nucleus (AGN), which is spatially coincident with the high-redshift galaxy, GHZ 9. The SMBH candidate is inferred to have a bolometric luminosity of (1.0−0.4+0.5)×1046ergs−1 , which corresponds to a black hole (BH) mass of (8.0−3.2+3.7)×107M⊙ assuming Eddington-limited accretion. This extreme mass at such an early cosmic epoch suggests the heavy seed origin for this BH candidate. Based on the Chandra and JWST discoveries of extremely high-redshift quasars, we have constructed the first simple AGN luminosity function extending to z ≈ 10. Comparison of this luminosity function with theoretical models indicates an overabundant z ≈ 10 SMBH population, consistent with a higher-than-expected seed formation efficiency.

Crossing the Rubicon of Reionization with z ∼ 5 QSOs

One of the key open questions in cosmology is the nature of sources that completed cosmological hydrogen reionization at z ∼ 5.2. High-z primeval galaxies have been long considered the main drivers of reionization, with a minor role played by high-z active galactic nuclei (AGN). However, in order to confirm this scenario, it is fundamental to measure the photoionization rate produced by active SMBHs close to the Epoch of Reionization. Given the pivotal role played by spectroscopically complete observations of high-z QSOs, in this paper we present the first results of the RUBICON (Reionizing the Universe with Bright Cosmological Nuclei) survey. It consists of a color-selected sample of bona fide z ∼ 5 QSO candidates from the Hyper Suprime-Cam Subaru Strategic Survey. Our QSO candidates have been validated both by photometric redshifts based on spectral energy distribution fitting and by spectroscopic redshifts, confirming that they lie at 4.5 < z spec < 5.2. A relatively high space density of QSOs (Φ ∼ 1.4 × 10−8 c Mpc−3) is thus confirmed at z ∼ 5 and M 1450 ∼ −27, consistent with a pure density evolution of the AGN luminosity function from z = 4 to z = 5, with a mild density evolution rate of 0.25 dex. This indicates that AGN could play a nonnegligible role in cosmic reionization. The Rubicon of reionization has been crossed.

Uchuu–ν2GC galaxies and AGN: cosmic variance forecasts of high-redshift AGN for JWST, Euclid, and LSST

ABSTRACT Measurements of the luminosity function of active galactic nuclei (AGN) at high redshift (z ≳ 6) are expected to suffer from field-to-field variance, including cosmic and Poisson variances. Future surveys, such as those from the Euclid telescope and JWST, will also be affected by field variance. We use the Uchuu simulation, a state-of-the-art cosmological N-body simulation with 2.1 trillion particles in a volume of 25.7 Gpc3, combined with a semi-analytic galaxy and AGN formation model, to generate the Uchuu–ν2GC catalogue, publicly available, that allows us to investigate the field-to-field variance of the luminosity function of AGN. With this Uchuu–ν2GC model, we quantify the cosmic variance as a function of survey area, AGN luminosity, and redshift. In general, cosmic variance decreases with increasing survey area and decreasing redshift. We find that at z ∼ 6 − 7, the cosmic variance depends weakly on AGN luminosity. This is because the typical mass of dark matter haloes in which AGN reside does not significantly depend on luminosity. Due to the rarity of AGN, Poisson variance dominates the total field-to-field variance, especially for bright AGN. We also examine the effect of parameters related to galaxy formation physics on the field variance. We discuss uncertainties present in the estimation of the faint-end of the AGN luminosity function from recent observations, and extend this to make predictions for the expected number of AGN and their variance for upcoming observations with Euclid, JWST, and the Legacy Survey of Space and Time (LSST).

The total rest-frame UV luminosity function from 3 &amp;lt; z &amp;lt; 5: a simultaneous study of AGN and galaxies from −28 &amp;lt; MUV &amp;lt; −16

ABSTRACT We present measurements of the rest-frame ultraviolet luminosity function (LF) at redshifts z = 3, z = 4, and z = 5, using 96894, 38655, and 7571 sources, respectively, to map the transition between active galactic nuclei (AGN) and galaxy-dominated ultraviolet emission shortly after the epoch of reionization (EoR). Sources are selected using a comprehensive photometric redshift approach, using 10 $\, {\rm deg}^2$ of deep extragalactic legacy fields covered by both HSC and VISTA. The use of template fitting spanning a wavelength range of 0.3–2.4 μm achieves 80–90 per cent completeness, much higher than the classical colour–colour cut methodology. The measured LF encompasses −26 &amp;lt; MUV &amp;lt; −19.25. This is further extended to −28.5 &amp;lt; MUV &amp;lt; −16 using complementary results from other studies, allowing for the simultaneous fitting of the combined AGN and galaxy LF. We find that there are fewer UV luminous galaxies (MUV &amp;lt; −22) at z ∼ 3 than z ∼ 4, indicative of an onset of widespread quenching alongside dust obscuration, and that the evolution of the AGN LF is very rapid, with their number density rising by around two orders of magnitude from 3 &amp;lt; z &amp;lt; 6. It remains difficult to determine if a double power law functional form is preferred over the Schechter function to describe the galaxy UV LF. Estimating the hydrogen ionizing photon budget from our UV LFs, we find that AGN can contribute to, but cannot solely maintain, the reionization of the Universe at z = 3–5. However, the rapidly evolving AGN LF strongly disfavours a significant contribution within the EoR.

Constraints From Dwarf Galaxies on Black Hole Seeding and Growth Models With Current and Future Surveys

Dwarf galaxies are promising test beds for constraining models of supermassive and intermediate-mass black holes (MBHs) via their BH occupation fraction (BHOF). Disentangling seeding from the confounding effects of mass assembly over a Hubble time is a challenging problem that we tackle in this study with a suite of semianalytical models (SAMs). We show how the measured BHOF depends on the lowest BH mass or active galactic nucleus (AGN) luminosity achieved by a survey. To tell seeding models apart, we need to detect or model all AGNs brighter than 1037 erg s−1 in galaxies of M * ∼ 108−10 M ⊙. Shallower surveys, like eRASS, cannot distinguish between seed models even with the compensation of a much larger survey volume. We show that the AMUSE survey, with its inference of the MBH population underlying the observed AGNs, strongly favors heavy seed models, growing with either a power-law Eddington ratio distribution function or one in which BH accretion is tied to the star formation rate (i.e., the AGN-main sequence, AGN-MS, model). These two growth channels can then be distinguished by the AGN luminosity function at >1040 erg s−1, with the AGN-MS model requiring more accretion than observed at z ∼ 0. Thus, current X-ray observations favor heavy seeds whose Eddington ratios follow a power-law distribution. The different models also predict different radio scaling relations, which we quantify using the fundamental plane of BH activity. We close with recommendations for the design of upcoming multiwavelength campaigns that can optimally detect MBHs in dwarf galaxies.

Constraints on the X-ray luminosity function of AGN at z = 5.7–6.4 with the Extragalactic Serendipitous Swift Survey

ABSTRACT X-ray luminosity functions (XLFs) of active galactic nuclei (AGNs) trace the growth and evolution of supermassive black hole populations across cosmic time. However, current XLF models are poorly constrained at redshifts of z &amp;gt; 6, with a lack of spectroscopic constraints at these high redshifts. In this work, we place limits on the bright-end of the XLF at z = 5.7–6.4 using high-redshift AGN identified within the Extragalactic Serendipitous Swift Survey (ExSeSS) catalogue. Within ExSeSS, we find one serendipitously X-ray detected z &amp;gt; 6 AGN, ATLAS J025.6821-33.4627, with an X-ray luminosity of $L_\mathrm{X}=8.47^{+3.40}_{-3.13}\times 10^{44}\mathrm{erg.s^{-1}}$ and z = 6.31 ± 0.03, making it the highest redshift, spectroscopically confirmed, serendipitously X-ray detected quasar known to date. We also calculate an upper limit on the space density at higher luminosities where no additional sources are found, enabling us to place constraints on the shape of the XLF. Our results are consistent with the rapid decline in the space densities of high-luminosity AGN towards high redshift as predicted by extrapolations of existing parametric models of the XLF. We also find that our X-ray based measurements are consistent with estimates of the bolometric quasar luminosity function based on UV measurements at z ≳ 6, although they require a large X-ray to bolometric correction factor at these high luminosities.

A new method to determine X-ray luminosity functions of AGN and their evolution with redshift

ABSTRACTAlmost all massive galaxies today are understood to contain supermassive black holes (SMBH) at their centres. SMBHs grew by accreting material from their surroundings, emitting X-rays as they did so. X-ray luminosity functions (XLFs) of active galactic nuclei (AGN) have been extensively studied in order to understand the AGN population’s cosmological properties and evolution. We present a new fixed rest-frame method to achieve a more accurate study of the AGN XLF evolution over cosmic time. Normally, XLFs are constructed in a fixed observer-frame energy band, which can be problematic because it probes different rest-frame energies at different redshifts. In the new method, we construct XLFs in the fixed rest-frame band instead, by varying the observed energy band with redshift. We target a rest-frame 2–8 keV band using XMM-Newton and HEAO 1 X-ray data, with seven observer-frame energy bands that vary with redshift for 0 &amp;lt; z &amp;lt; 3. We produce the XLFs using two techniques; one to construct a binned XLF, and one using a maximum likelihood (ML) fit, which makes use of the full unbinned source sample. We find that our ML best-fitting pure luminosity evolution results for both methods are consistent with each other, suggesting that performing XLF evolution studies with the high-redshift data limited to high-luminosity AGN is not very sensitive to the choice of fixed observer-frame or rest-frame energy band, which is consistent with our expectation that high-luminosity AGN typically show little absorption. We have demonstrated the viability of the new method in measuring the XLF evolution.

On the Coevolution of the AGN and Star-forming Galaxy Ultraviolet Luminosity Functions at 3 < z < 9

Studies of the rest-frame ultraviolet (UV) luminosity functions (LFs) typically treat star-forming galaxies and active galactic nuclei (AGNs) separately. However, modern ground-based surveys now probe volumes large enough to discover AGNs at depths sensitive enough for fainter galaxies, bridging these two populations. Using these observations as constraints, we present a methodology to empirically jointly model the evolution of the rest-UV LFs at z = 3–9. Our critical assumptions are that both populations have LFs well described by double power laws modified to allow for a flattening at the faint-end, and that all LF parameters evolve smoothly with redshift. This methodology provides a good fit to the observations and makes predictions to volume densities not yet observed, finding that the volume density of bright (M UV = −28) AGNs rises by ∼105 from z = 9 to z = 3, while bright (M UV = −21) star-forming galaxies rise by only ∼102 across the same epoch. The observed bright-end flattening of the z = 9 LF is unlikely to be due to AGN activity, and rather is due to a shallowing of the bright-end slope, implying a reduction of feedback in bright galaxies at early times. The intrinsic ionizing emissivity is dominated by star-forming galaxies at z > 3, even after applying a notional escape fraction. We find decent agreement between our AGN LFs and predictions based on different black hole seeding models, though all models underpredict the observed abundance of bright AGNs. We show that the wide-area surveys of the upcoming Euclid and Roman observatories should be capable of discovering AGNs to z ∼ 8.

Dwarf AGNs from variability for the origins of seeds (DAVOS): Intermediate-mass black hole demographics from optical synoptic surveys

ABSTRACT We present a phenomenological forward Monte Carlo model for forecasting the population of active galactic nuclei (AGNs) in dwarf galaxies observable via their optical variability. Our model accounts for expected changes in the spectral energy distribution of AGNs in the intermediate-mass black hole (IMBH) mass range and uses observational constraints on optical variability as a function of black hole (BH) mass to generate mock light curves. Adopting several different models for the BH occupation function, including one for off-nuclear IMBHs, we quantify differences in the predicted local AGN mass and luminosity functions in dwarf galaxies. As a result, we are able to model the fraction of variable AGNs as a function of important galaxy host properties, such as host galaxy stellar mass, in the presence of selection effects. We find that our adopted occupation fractions for the ‘heavy’ and ‘light’ initial BH seeding scenarios can be distinguished with variability at the 2–3σ level for galaxy host stellar masses below ∼108M⊙ with data from the upcoming Vera C. Rubin Observatory. We also demonstrate the prevalence of a selection bias whereby recovered IMBH masses fall, on average, above the predicted value from the local host galaxy–BH mass scaling relation with the strength of this bias dependent on the survey sensitivity. Our methodology can be used more broadly to calibrate AGN demographic studies in synoptic surveys. Finally, we show that a targeted ∼ hourly cadence program over a few nights with the Rubin Observatory can provide strong constraints on IMBH masses given their expected rapid variability time-scales.

Consistent Analysis of the AGN LF in X-Ray and MIR in the XMM-LSS Field

Abstract The luminosity function of active galactic nuclei (AGN) probes the history of supermassive black hole assembly and growth across cosmic time. To mitigate selection biases, we present a consistent analysis of the AGN luminosity functions (LFs) derived for both X-ray and mid-infrared (MIR) selected AGN in the XMM-Large Scale Structure field. There are 4268 AGN used to construct the MIR luminosity function (IRLF) and 3427 AGN used to construct the X-ray luminosity function (XLF), providing the largest census of the AGN population out to z = 4 in both bands with significant reduction in uncertainties. We are able for the first time to see the knee of the IRLF at z &gt; 2 and observe a flattening of the faint-end slope as redshift increases. The bolometric luminosity density, a proxy for the cosmic black hole accretion history, computed from our LFs, shows a peak at z ≈ 2.25, consistent with recent estimates of the peak in the star formation rate density (SFRD). However, at earlier epochs, the AGN luminosity density is flatter than the SFRD. If confirmed, this result suggests that the build up of black hole mass outpaces the growth of stellar mass in high-mass systems at z ≳ 2.5. This is consistent with observations of redshift z ∼ 6 quasars that lie above the local M − σ relationship. The luminosity density derived from the IRLF is higher than that from the XLF at all redshifts. This is consistent with the dominant role of obscured AGN activity in the cosmic growth of supermassive black holes.

The Space Density of Ultra-luminous QSOs at the End of Reionization Epoch by the QUBRICS Survey and the AGN Contribution to the Hydrogen Ionizing Background

Abstract Motivated by evidences favoring a rapid and late hydrogen reionization process completing at z ∼ 5.2–5.5 and mainly driven by rare and luminous sources, we have reassessed the estimate of the space density of ultra-luminous QSOs at z ∼ 5 in the framework of the QUBRICS survey. A ∼ 90% complete sample of 14 spectroscopically confirmed QSOs at M 1450 ≤ −28.3 and 4.5 ≤ z ≤ 5.0 has been derived in an area of 12,400 deg2, thanks to multiwavelength selection and Gaia astrometry. The space density of z ∼ 5 QSOs within −29.3 ≤ M 1450 ≤ −28.3 is three times higher than previous determinations. Our results suggest a steep bright-end slope for the QSO luminosity function at z ∼ 5 and a mild redshift evolution of the space density of ultrabright QSOs (M 1450 ∼ −28.5) at 3 &lt; z &lt; 5.5, in agreement with the redshift evolution of the much fainter active galactic nucleus (AGN) population at M 1450 ∼ −23. These findings are consistent with a pure density evolution for the AGN population at z &gt; 3. Adopting our z ∼ 4 QSO luminosity function and applying a mild density evolution in redshift, a photoionization rate of Γ HI = 0.46 − 0.09 + 0.17 × 10 − 12 s − 1 has been obtained at z = 4.75, assuming an escape fraction of ∼70% and a steep faint-end slope of the AGN luminosity function. The derived photoionization rate is ∼50–100% of the ionizing background measured at the end of the reionization epoch, suggesting that AGNs could play an important role in the cosmological reionization process.

The impact of black hole feedback on the UV luminosity and stellar mass assembly of high-redshift galaxies

ABSTRACT We employ the Delphi semi-analytical model to study the impact of black hole growth on high-redshift galaxies, both in terms of the observed ultraviolet (UV) luminosity and of the star formation rate. To do this, first, we assess the contribution of active galactic nuclei (AGN) to the total galaxy UV luminosity as a function of stellar mass and redshift. We find that for MUV ≲ −24 mag and z ≈ 5−6, the galaxies for which the black hole UV luminosity outshines the stellar UV emission become the majority, and we estimate their duty cycle. Secondly, we study the evolution of the AGN and stellar luminosity functions (LFs), finding that it is driven both by changes in their characteristic luminosities (i.e. evolution of the intrinsic brightness of galaxies) and in their normalizations (i.e. evolution of the number densities of galaxies), depending on the luminosity range considered. Finally, we follow the mass assembly history for three different halo mass bins, finding that the magnitude of AGN-driven outflows depends on the host halo mass. We show that AGN feedback is most effective when the energy emitted by the accreting black hole is approximately $1{{\ \rm per\ cent}}$ of the halo binding energy, and that this condition is met in galaxies in haloes with $M_\mathrm{h} \sim 10^{11.75} \, {\rm \rm M_\odot }$ at z = 4. In such cases, AGN feedback can drive outflows that are up to 100 times more energetic than supernovae-driven outflows, and the star formation rate is a factor of 3 lower than for galaxies of the same mass without black hole activity.

Diffuse flux of PeV neutrinos from centrifugally accelerated protons in active galactic nuclei

Evidence for high-energy astrophysical PeV neutrinos has been found in the IceCube experiment from an analysis with 7.5-year (2010–2017) data. Active galactic nuclei (AGN) are among the most prominent objects in the universe, and are widely speculated to be emitters of ultra-high-energy (UHE) cosmic rays with proton domination. Based on the standard two-step LLCD mechanism of particle acceleration, a transformation of energy occurs from AGN's central super-massive black hole (SMBH) rotation to high-energy protons. Protons can be accelerated up to energies and above, and might generate PeV neutrinos in the energy range through plausible hadronic interactions. The theoretically estimated revised extragalactic diffuse muon neutrino flux employing the “luminosity-dependent density evolution (LDDE)” model for the AGN luminosity function (LF) is found consistent with the IceCube level if only a fraction, of the total bolometric luminosity (BL) of AGN is being realizable to power the PeV neutrinos. In the Λ CDM cosmological framework with the LDDE-modeled LF and photon index distribution, about of the total BL is enough to power the IceCube neutrinos.

Supermassive black holes in cosmological simulations – II: the AGN population and predictions for upcoming X-ray missions

ABSTRACT In large-scale hydrodynamical cosmological simulations, the fate of massive galaxies is mainly dictated by the modelling of feedback from active galactic nuclei (AGNs). The amount of energy released by AGN feedback is proportional to the mass that has been accreted on to the black holes (BHs), but the exact subgrid modelling of AGN feedback differs in all simulations. While modern simulations reliably produce populations of quiescent massive galaxies at z ≤ 2, it is also crucial to assess the similarities and differences of the responsible AGN populations. Here, we compare the AGN populations of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations. The AGN luminosity function (LF) varies significantly between simulations. Although in agreement with current observational constraints at z = 0, at higher redshift the agreement of the LFs deteriorates with most simulations producing too many AGNs of $L_{\rm x, 2\!-\!10 \, keV}\sim 10^{43\!-\!44}\, \rm erg\, s^{-1}$. AGN feedback in some simulations prevents the existence of any bright AGN with $L_{\rm x, 2\!-\!10 \, keV}\geqslant 10^{45}\rm \,erg\, s^{-1}$ (although this is sensitive to AGN variability), and leads to smaller fractions of AGN in massive galaxies than in the observations at z ≤ 2. We find that all the simulations fail at producing a number density of AGN in good agreement with observational constraints for both luminous ($L_{\rm x, 2\!-\!10 \, keV}\sim 10^\text{43-45}\, \rm erg\, s^{-1}$) and fainter ($L_{\rm x, 2\!-\!10 \, keV}\sim 10^\text{42-43}\, \rm erg\, s^{-1}$) AGNs and at both low and high redshifts. These differences can aid us in improving future BH and galaxy subgrid modelling in simulations. Upcoming X-ray missions (e.g. Athena, AXIS, and LynX) will bring faint AGNs to light and new powerful constraints. After accounting for AGN obscuration, we find that the predicted number density of detectable AGNs in future surveys spans at least one order of magnitude across the simulations, at any redshift.

Ionizing radiation from AGNs at z &amp;gt; 3.3 with the Subaru Hyper Suprime-Cam Survey and the CFHT Large Area U-band Deep Survey (CLAUDS)

ABSTRACT We use deep and wide imaging data from the CFHT Large Area U-band Deep Survey (CLAUDS) and the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) to constrain the ionizing radiation (Lyman Continuum, LyC) escape fraction from active galactic nuclei (AGNs) at z ∼ 3–4. For 94 AGNs with spectroscopic redshifts at 3.3 &amp;lt; z &amp;lt; 4.0, we use their U/i-band flux ratios to estimate LyC transmission of individual AGNs. The distribution of their LyC transmission shows values lower than the range of LyC transmission values for IGM of the same redshift range, which suggests that LyC escape fraction of AGNs at z &amp;gt; 3.3 is considerably lower than unity in most cases. We do not find any trend in LyC transmission values depending on their UV luminosities. Based on the photometry of stacked images we find the average flux ratio of LyC and non-ionizing UV photons escaping from the objects (fLyC/fUV)out = 0.182 ± 0.043 for AGNs at 3.3 &amp;lt; z &amp;lt; 3.6, which corresponds to LyC escape fraction fesc = 0.303 ± 0.072 if we assume a fiducial intrinsic SED of AGN. Based on the estimated LyC escape fraction and the UV luminosity function (UVLF) of AGNs, we argue that UV-selected AGNs’ contribution to the LyC emissivity at the epoch is minor, although the size of their contribution largely depends on the shape of the UVLF.

Modelling radio luminosity functions of radio-loud AGN by the cosmological evolution of supermassive black holes

Aims. We develop a formalism to model the luminosity functions (LFs) of radio-loud active galactic nuclei (AGN) at GHz frequencies by the cosmological evolution of the supermassive black hole (SMBH) hosted in their nuclei. The mass function and Eddington ratio distributions of SMBHs computed in a previous work published by one of the authors have been taken as the starting point for this analysis. Methods. Our approach is based on physical and phenomenological relations that allow us to statistically calculate the radio luminosity of AGN cores, corrected for beaming effects, by linking it with the SMBH at their centre, through the fundamental plane of black hole activity. Moreover, radio luminosity from extended jets and lobes is also computed through a power-law relationship that reflects the expected correlation between the inner radio core and the extended jets and lobes. By following a classification scheme well established in the field, radio-loud AGN are further divided into two classes, characterized by different accretion modes onto the central BH. If the Eddington ratio, λ, is ≤0.01 they are called low-kinetic (LK) mode AGN; if λ ≥ 0.01, they are called high-kinetic (HK) mode AGN, this critical value roughly corresponding to the transition between radiatively inefficient and efficient accretion flows. The few free parameters used in the present model are determined by fitting two different types of observational data sets: local (or low-redshift) LFs of radio-loud AGN at 1.4 GHz and differential number counts of extragalactic radio sources at 1.4 and 5 GHz. Results. Our present model fits well almost all published data on LFs of LK mode AGN and of the total AGN population up to redshifts z ≤ 1.5 and also in the full range of luminosities currently probed by data. On the other hand, it tends to underestimate some recent measures of the LF of HK mode AGN at low redshifts, but only at low radio luminosities. All in all, the good performance of our model in this redshift range is remarkable, considering that all the free parameters used but the fraction of HK mode AGN are redshift independent. The present model is also able to provide a very good fit to almost all data on number counts of radio-loud sources at 1.4 and 5 GHz.

The rapid transition from star formation to AGN-dominated rest-frame ultraviolet light at z ≃ 4

ABSTRACT With the advent of deep optical-to-near-infrared extragalactic imaging on the degree scale, samples of high-redshift sources are being selected that contain both bright star-forming (SF) galaxies and faint active galactic nuclei (AGN). In this study, we investigate the transition between SF- and AGN-dominated systems at z ≃ 4 in the rest-frame ultraviolet (UV). We find a rapid transition to AGN-dominated sources brightward of MUV ≃ −23.2. The effect is observed in the rest-frame UV morphology and size–luminosity relation, where extended clumpy systems become point-source-dominated, and also in the available spectra for the sample. These results allow us to derive the rest-frame UV luminosity function (LF) for the SF- and AGN-dominated subsamples. We find the SF-dominated LF is best fit with a double power law, with a lensed Schechter function being unable to explain the existence of extremely luminous SF galaxies at MUV ≃ −23.5. If we identify AGN-dominated sources according to a point-source morphology criterion, we recover the relatively flat faint-end slope of the AGN LF determined in previous studies. If we instead separate the LF according to the current spectroscopic AGN fraction, we find a steeper faint-end slope of α = −1.83 ± 0.11. Using a simple model to predict the rest-frame AGN LF from the z = 4 galaxy LF, we find that the increasing impact of host galaxy light on the measured morphology of faint AGN can explain our observations.

Accretion History of AGNs. III. Radiative Efficiency and AGN Contribution to Reionization

Abstract The cosmic history of supermassive black hole (SMBH) growth is important for understanding galaxy evolution, reionization, and the physics of accretion. Recent NuSTAR, Swift-BAT, and Chandra hard X-ray surveys have provided new constraints on the space density of heavily obscured active galactic nuclei (AGNs). Using the new X-ray luminosity function derived from these data, we here estimate the accretion efficiency of SMBHs and their contribution to reionization. We calculate the total ionizing radiation from AGNs as a function of redshift, based on the X radiation and distribution of obscuring column density, converted to ultraviolet (UV) wavelengths. Limiting the luminosity function to unobscured AGNs only, our results agree with current UV luminosity functions of unobscured AGNs. For realistic assumptions about the escape fraction, the contribution of all AGNs to cosmic reionization is ∼4 times lower than the galaxy contribution (23% at z ∼ 6). Our results also offer an observationally constrained prescription that can be used in simulations or models of galaxy evolution. To estimate the average efficiency with which SMBHs convert mass to light, we compare the total radiated energy, converted from X-ray light using a bolometric correction, with the most recent local black hole mass density. The most likely value, η ∼ 0.3–0.34, approaches the theoretical limit for a maximally rotating Kerr black hole, η = 0.42, implying that on average growing SMBHs are spinning rapidly.