Next Article in Journal
Effects of Two Quantum Correction Parameters on Chaotic Dynamics of Particles near Renormalized Group Improved Schwarzschild Black Holes
Previous Article in Journal
Fitting the Crab Supernova with a Gamma-Ray Burst
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Surveying the Onset and Evolution of Supermassive Black Holes at High-z with AXIS

1
Department of Physics, University of Miami, Coral Gables, FL 33124, USA
2
Department of Physics, University of Maryland Baltimore County, 1000 Hilltop Cir, Baltimore, MD 21250, USA
3
Department of Physics and Astronomy, Clemson University, Kinard Lab of Physics, Clemson, SC 29634, USA
4
Dipartimento di Fisica e Astronomia (DIFA), Università di Bologna, Via Gobetti 93/2, 40129 Bologna, Italy
5
INAF-Osservatorio di Astrofisica e Scienza dello Spazio, Via Piero Gobetti 93/3, 40129 Bologna, Italy
6
Center for Astrophysics, Harvard & Smithsonian, 60 Garden St., Cambridge, MA 02138, USA
7
Black Hole Initiative, Harvard University, 20 Garden St., Cambridge, MA 02138, USA
8
Zentrum für Astronomie der Universität Heidelberg, ITA, Albert-Ueberle-Str. 2, 69120 Heidelberg, Germany
9
Max-Planck-Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany
10
Leibniz-Institut fur Astrophysik Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany
11
Eureka Scientific, 2452 Delmer Street Suite 100, Oakland, CA 94602-3017, USA
12
Department of Astronomy and Joint Space-Science Institute, University of Maryland, College Park, MD 20742, USA
*
Author to whom correspondence should be addressed.
Universe 2024, 10(7), 276; https://doi.org/10.3390/universe10070276
Submission received: 17 April 2024 / Revised: 29 May 2024 / Accepted: 21 June 2024 / Published: 25 June 2024
(This article belongs to the Section Galaxies and Clusters)

Abstract

:
The nature and origin of supermassive black holes (SMBHs) remain an open matter of debate within the scientific community. While various theoretical scenarios have been proposed, each with specific observational signatures, the lack of sufficiently sensitive X-ray observations hinders the progress of observational tests. In this white paper, we present how AXIS will contribute to solving this issue. With an angular resolution of 1.5″ on-axis and minimal off-axis degradation, we designed a deep survey capable of reaching flux limits in the [0.5–2] keV range of approximately 2 × 10−18 erg s−1 cm−2 over an area of 0.13 deg2 in approximately 7 million seconds (7 Ms). Furthermore, we planned an intermediate depth survey covering approximately 2 deg2 and reaching flux limits of about 2 × 10−17 erg s−1 cm−2 in order to detect a significant number of SMBHs with X-ray luminosities (LX) of approximately 1042 erg s−1 up to z∼10. These observations will enable AXIS to detect SMBHs with masses smaller than 105 M, assuming Eddington-limited accretion and a typical bolometric correction for Type II AGN. AXIS will provide valuable information on the seeding and population synthesis models of SMBHs, allowing for more accurate constraints on their initial mass function (IMF) and accretion history from z∼0–10. To accomplish this, AXIS will leverage the unique synergy of survey telescopes such as the JWST, Roman, Euclid, Vera Rubin Telescope, and the new generation of 30 m class telescopes. These instruments will provide optical identification and redshift measurements, while AXIS will discover the smoking gun of nuclear activity, particularly in the case of highly obscured AGN or peculiar UV spectra as predicted and recently observed by the JWST in the early Universe.

1. Introduction

The origin of the very first SMBH in the Universe remains an open puzzle in astrophysics, attracting significant theoretical attention for the past two decades. The scarcity of observational data to address this question has added to its intrigue. Furthermore, the detection of numerous quasars powered by SMBHs that accrete matter at both high ( z > 6 ) and low redshifts ( z < 6 ) has underscored the importance of understanding their origins. It is predicted that the remnants of the first stars would have given rise to a population of relatively light black hole (BH) seeds during the early cosmic epochs. However, the key question remains: what is the actual mass of these remnants?
The reason for this puzzle is the detection with, e.g., UKIRT, DESI, and PanSTARRS of SMBHs with masses exceeding 109 M at redshifts (z) around 7 (see, e.g., [1,2,3]) and preliminary Chandra and JWST evidence of 107 M SMBHs up to z∼10 [4,5]. These observations have raised an intriguing paradox. The conventional mechanisms for stellar remnants to grow into such massive SMBHs within the limited available time seem inadequate. This inconsistency has spurred the investigation of alternative explanations, including the presence of heavy seeds like Direct Collapse Black Holes (DCBHs, see, e.g., [6,7]), remnants of relatively massive but lighter Population III stars (POPIII), or from BH mergers in compact clusters, or the collapse of very massive stars formed through mergers of stars [8].
DCBHs are hypothesized to originate directly from the collapse of primordial gas clouds, bypassing the conventional process of star formation. These unique entities, formed with substantial masses from the outset, offer a compelling explanation for the swift growth of SMBHs observed during the early stages of cosmic history.
Likewise, remnants of POPIII stars, which are believed to have formed in environments lacking heavy elements, may retain significant masses as they undergo stellar evolution. Although these remnants are lighter than DCBHs, they still have the potential to act as seeds for the subsequent growth into supermassive entities.
In recent times, the notion of Primordial Black Holes (see, e.g., [9,10,11]) has regained prominence as another plausible solution to the enigma of the formation of SMBHs. PBHs are black holes that might have emerged directly from the extreme density fluctuations present in the early Universe right after inflation. If PBHs exist, they could offer a plentiful supply of massive seeds, facilitating the rapid growth of SMBHs within the observable timespan.
The resurrection of the PBH concept, along with the consideration of DCBHs and remnants of POPIII stars, underscores the capacity for alternative mechanisms to explain the early formation and growth of SMBHs. These ideas challenge the traditional paradigm and call for further theoretical investigations, observational constraints, and potential detection methods to gain deeper insights into the enigma surrounding the origins of SMBHs in the cosmos.
In Figure 1, we present the growth of two hypothetical SMBHs that are seeded either from a “light” M∼10–100 M seed or from a “heavy” M∼104–6 M seed, in comparison with a collection of quasar detections at z > 6. The plot is compared with the sensitivity limit of a 7 Ms AXIS observation (see below). Assuming an indicative seed collapse at z = 25, the heavy seeds are capable of producing the observed high-redshift SMBH masses without requiring prolonged periods of super-Eddington accretion. In contrast, if SMBHs originated solely from light seeds, an extended phase of sustained super-Eddington growth would be necessary. However, in the case of PBHs, these assumptions are relaxed due to the considerably longer time available for black hole growth or, as suggested by certain models [10,12], the possibility of PBHs forming with very large masses at the e + e annihilation phase transition, a few seconds after the Big Bang.
The presence of X-ray emissions is widely considered to be a definitive indicator of the existence of an actively growing supermassive black hole (SMBH) in a galaxy, observed as AGN. Current, highly sensitive X-ray telescopes such as Chandra and XMM-Newton have played a crucial role in advancing our understanding of the growth and evolution of AGN up to redshifts of z∼4–5. Surveys conducted with Chandra have achieved deep [0.5–2] keV flux measurements on the order of 10−17 erg s−1 cm−2, although limited to relatively small fields spanning a few square arcminutes. This limitation is due to the rapid degradation of the Point Spread Function (PSF) with increasing off-axis angles and the presence of strong vignetting. On the other hand, XMM-Newton offers a more stable but larger PSF but has an angular resolution typically exceeding 5″ on-axis, limiting its sensitivity due to source confusion. Moreover, Chandra’s performance degradation has prevented it from reaching much deeper fluxes than what it has already obtained in the CDFS [16] with a reasonable time investment; therefore, the discovery space in the limiting flux–area space has stalled over the last decade.
With a proposed on-axis angular resolution of 1.5″ and an averaged PSF Half Energy Width (HEW) of 1.6″, coupled with an effective area of 4200 cm2 at 1 keV, 830 cm2 at 6 keV, and a 24′ diameter active field of view (i.e., ∼0.13 deg2), AXIS [17] is poised to be a cutting-edge “X-ray survey machine” for the next decade. Benefitting from the low background provided by its low Earth orbit, AXIS will surpass Chandra by achieving depths that are one order of magnitude greater over a high-angular-resolution area that is at least 10 times larger. The net result is an effective grasp (the number of detected sources per unit area) of the faintest fluxes that are at least two (three) orders of magnitude better than that of Chandra (XMM-Newton). As of the last release of the Chandra Source Catalog, (CSC, Martinez Galarza [18]), over 25 years of operation, Chandra covered about 700 deg2 of sky and detected about 400,000 sources with only a few hundreds at the faintest fluxes over of a few square arcminutes of survey area [16]. A similar number of sources was detected, over 1150 deg2, at around one order of magnitude shallower flux limit than Chandra, by XMM-Newton as of the latest release of the 4XMM Catalog [19]. In Figure 2, we present a simulated deep survey conducted by AXIS, allowing for a comparison of its imaging capabilities with those of Chandra.
In this paper, we present how the proposed AXIS probe will revolutionize the field of high-z X-ray surveys and the quest to understand the nature and growth of SMBHs.

2. The Role of X-ray Surveys in Understanding the Nature and the Evolution of SMBHs

X-ray surveys have brought about a revolutionary understanding of the evolution of SMBHs. It is now evident that in the redshift range of approximately 0 < z < 5, we observe a phenomenon known as “downsizing”. This implies that the highest density and volume emissivity of the most rapidly growing, luminous, and massive SMBHs occur earlier in cosmic time (around z = 3–4) compared to the less luminous, slowly growing, or less massive active SMBHs that dominate the Universe at z = 0–1.

2.1. Obscuration as a Challenge to Understand SMBH Growth

This downsizing trend is probably influenced by a combination of factors. According to the unified model of AGN, the accretion process is frequently concealed by a cold absorber, such as the torus or interstellar material. The intense radiation from luminous quasars can rapidly dissociate molecules, allowing radiation to escape rapidly. Conversely, if an SMBH accretes at a low rate or has a lower mass, the process is prolonged, and the line of sight remains obstructed for an extended period, delaying the clearing of the absorber.
The observed downsizing phenomenon could be attributed to the presence of cold absorbers, the ionizing radiation from luminous quasars, and the rate of accretion and mass of the SMBHs themselves. Understanding these mechanisms is crucial to unraveling the intricate processes that govern the growth and evolution of SMBH populations throughout cosmic history.
Indeed, recent studies have revealed an intriguing trend with regard to obscured quasars, indicating an increasing fraction of them with redshift. In particular, most (over 75%) of the quasars at z > 3–4 exhibit significant obscuration due to absorbing material with column densities N H > 1023 (e.g., [20]). These findings highlight the prevalence of heavily obscured sources at higher redshifts, which are very difficult to detect with restframe optical or UV observations.

2.2. High Angular Resolution X-ray Surveys to Unveil the Elusive Nature of High-z AGN

Preliminary results from early observations with the James Webb Space Telescope (JWST) [21] have provided valuable insights into the nature of obscured sources. Specifically, at z = 5, the number of Type II sources is estimated to be at least three times higher than previously inferred from observations conducted by the Chandra X-ray Observatory [22]. This discrepancy can be attributed to two primary factors.
First, the volume covered by high-resolution Chandra observations is relatively limited, effectively constrained to the inner few on-axis arcminutes of the survey. Consequently, these observations do not capture the rarer lower volume density and more luminous sources, potentially leading to underestimated numbers.
Secondly, these highly obscured sources are predicted to be inherently faint in X-ray emission, with [0.5–2] keV flux levels below the detection threshold of approximately 1–5 × 10−17 erg s−1 cm−2. As a result, they may go undetected in the Chandra deep X-ray surveys, further contributing to the underestimation of their abundance.
Taken together, these factors highlight the necessity of combining observations from different wavelengths, such as infrared (e.g., the JWST) and X-rays (e.g., Chandra or future missions like AXIS), to obtain a comprehensive understanding of obscured sources. By leveraging the capabilities of both telescopes, we can overcome the limitations imposed by the small volume probed by Chandra and the faintness of highly obscured sources, thus providing a more accurate and complete picture of the population of obscured quasars across cosmic history.
At the same time, distinguishing highly obscured AGN from normal star-forming galaxies based solely on Near-Mid Infrared observations is indeed challenging. This difficulty arises due to the presence of highly degenerate features shared by both star formation processes (such as dust and narrow spectral lines) and obscured AGN. In particular, for obscured AGN, this similarity in features makes it extremely challenging to untangle the true nature of a given source using only infrared observations.
In a recent paper, Lyu et al. [23] showed that MIRI can, in principle, disentangle AGN from SFG. However, at z > 4, its sensitivity to obscuration features significantly decreases as the rest-frame wavelengths sampled are optical-NIR.
In their sample of MIRI-selected AGN, 80% are not X-ray detected. This can be due to a large number of Compton-thick sources or the intrinsic X-ray faintness of high-z AGN at the faint end of the luminosity function. On the other hand, Vito et al. [24] found no difference in the spectral shape of high-z QSOs, while [25] found steeper-than-average X-ray spectra. Finally, Maiolino et al. [26] presented a population of the JWST selected. X-ray weak sources where heavy X-ray absorption by clouds with a large (Compton-thick) column density and low dust content, such as the Broad Line Region (BLR) clouds, can explain the X-ray weakness. As for the extreme redshifts that AXIS plans to sample, Natarajan et al. [27] predicts extremely bright X-ray/optical ratios for the early phase of SMBH growth at z∼9.
This suggests that, in order to understand the nature of accretion at high-z and confirm the expected large fraction of Compton-thick AGN (see, e.g., [20,28,29,30]), further studies in the X-ray band, paired with other wavelengths, are necessary.
To overcome this challenge and gain a more accurate understanding, a multiwavelength approach is crucial. Combining observations from different spectral ranges, such as X-ray, optical, and infrared, allows for a more comprehensive analysis. X-ray surveys, such as those conducted by Chandra or future missions like AXIS, provide valuable information about the presence of AGN activity through the detection of X-ray emissions associated with the accretion onto SMBHs.
Furthermore, optical observations can provide insight into the presence of broad emission lines and other AGN signatures. By integrating data from multiple wavelengths with Near-Mid Infrared observations, a more thorough examination of sources can be achieved, leading to more accurate classifications and differentiation between highly obscured AGN and normal star-forming galaxies.
Indeed, a commonly employed technique for selecting AGN involves the use of rest-frame (observed NIR at high-z) optical spectroscopy, typically employing the BPT diagram. However, early results from the JWST indicate that this approach may not be as effective at high redshifts, as the spectra of AGN at these cosmic epochs differ from those commonly observed in the local Universe. For example, a recent study by Maiolino et al. [5] presents the case of AGN at z∼10.6, where an unusually bright Ne [IV] line is observed. This emission line is atypical for AGN [31], highlighting the distinct nature of high-redshift AGN spectra. Additionally, a recent survey of high-redshift broad-line AGN using the JWST has revealed a lack of [NII] emission [32]. As a consequence, these AGN fall within the parameter space characterized by high [OIII]/H β ratios and low [NII]/H α ratios, similar to galaxies at z≳ 4 without broad-line emission. This finding suggests that the low metallicity of their host galaxies prevents their differentiation from regular star-forming galaxies using the BPT diagram.
However, it is worth noting that at higher luminosities, this phenomenon seems to be less pronounced. Hence, low-mass SMBHs (on the order of 106–7 M) growing in low-metallicity, high-column-density galaxies necessitate confirmation in the X-ray band to be definitively labeled as AGN. X-ray observations, such as those conducted by AXIS, can play a crucial role in confirming the presence of AGN activity in these sources.

2.3. The Power of Combining AXIS and JWST Surveys

The distinctive spectroscopic features observed in high-redshift AGN, coupled with their occurrence in low-metallicity environments, pose challenges for their identification using traditional methods such as the BPT diagram. Complementary observations in the X-ray band are essential for providing definitive confirmation of AGN activity in these specific scenarios.
Indeed, the combination of the JWST and AXIS is of utmost importance in comprehending the nature of galaxies hosting the first AGN, as well as the coevolution of SMBHs and galaxies.
Soft X-ray surveys play a crucial role in this endeavor as obscured AGN, including Compton-thick sources, exhibit peak emissivity in the 10–20 keV energy range. At redshifts z > 6, this peak is shifted into the [0.5–2] keV energy range, precisely where focusing X-ray telescopes like AXIS possess their maximum effective area, as depicted in Figure 3.
This advantageous alignment allows AXIS to detect sources that the JWST may not be able to identify as AGN because of the relatively unbiased nature of X-ray selection. By conducting soft X-ray surveys, AXIS can identify and characterize obscured AGN even in cases where they may be missed or misclassified based on other observational criteria. This unbiased selection is essential for a comprehensive understanding of the AGN population and its connection to the early Universe.
The synergistic combination of the JWST and an X-ray telescope like AXIS enables a holistic view of galaxy evolution, shedding light on the intricate interplay between SMBHs and their host galaxies. By harnessing the power of multiwavelength observations, we can explore the coevolutionary processes that shape these systems throughout cosmic history.
To study SMBHs at redshifts z = 6–10, it is crucial to have an observatory capable of detecting sources with intrinsic X-ray luminosities LX > 1042 erg s−1, and log N(H) > 23 as most of the accretion in this epoch is expected to be obscured (see, e.g., [33]). The X-ray luminosities of highly obscured, Eddington-limited, M ≳ 7 × 105 M sources correspond to fluxes of approximately 1–4 × 10−18 erg s−1 cm−2 in the [0.5–2] keV energy band. Such sensitivity levels are necessary to probe the presence and properties of rapidly growing BHs in the high-reshift Universe.
AXIS, with its designed capabilities and sensitivity, is precisely tailored to achieve these required limits. Combined with the survey strategy outlined below, AXIS has the potential to serendipitously detect hundreds of SMBHs at redshifts z = 6–10 and even extend the search to z∼12 in the case of pointed observations targeting preselected objects.
Figure 3. Simulation of a 7 Ms observation of a Compton-thick AGN with an intrinsic luminosity of LX = 2 × 1044 erg s−1 at redshifts around 6 and 8. AXIS will detect hundreds of counts, enabling the measurement of column density (NH) at a significance level of 5 σ . The continuous lines represent the best fits obtained using the XSPEC MyTorus model. This level of statistical accuracy is comparable to that achieved by NuSTAR, that unambiguously provided Compton-thick AGN line-of-sight column density measurements the local Universe [34].
Figure 3. Simulation of a 7 Ms observation of a Compton-thick AGN with an intrinsic luminosity of LX = 2 × 1044 erg s−1 at redshifts around 6 and 8. AXIS will detect hundreds of counts, enabling the measurement of column density (NH) at a significance level of 5 σ . The continuous lines represent the best fits obtained using the XSPEC MyTorus model. This level of statistical accuracy is comparable to that achieved by NuSTAR, that unambiguously provided Compton-thick AGN line-of-sight column density measurements the local Universe [34].
Universe 10 00276 g003

2.4. SMBH and Host Galaxy Coevolution over Cosmic Time with AXIS

The combination of AXIS’s sensitivity to faint X-ray sources and its survey strategy will significantly advance our understanding of the early Universe, shedding light on the formation and evolution of SMBHs during these cosmic epochs. By enabling the detection of a large number of SMBHs at high redshifts, AXIS will contribute to our knowledge of the coevolution of galaxies and their central SMBH and its early onset, providing crucial insights into the early stages of cosmic structure formation.
Studying the growth of SMBHs from the earliest epochs of the Universe with redshifts greater than 6, through to the peak of star formation at redshifts 1–2, is crucial for unraveling the intertwined evolution of SMBHs and galaxies. During this period, most stars in the Universe formed, making it a pivotal phase in cosmic history.
Recent advances in our understanding of active galactic nuclei (AGN) have revealed that they are not constant sources of emission. Instead, they exhibit variability in their Eddington ratios on relatively short timescales, typically less than 105–6 years [35]. As a result, it is essential to explore the full range of accretion rates and luminosities across the entire cosmic timeline in order to comprehensively characterize the AGN population.
Additionally, based on new JWST data of galactic systems in the redshift range z = 4–7, with both SMBH and stellar mass measurements available, [36] recently found that SMBH masses are 10–100 times higher than what the host’s stellar mass would suggest, based on the local scaling relations. This study indicates a significant number of yet-to-be-detected SMBHs in current JWST surveys. They are still undetected because (i) they may be inactive (i.e., duty cycles 1–10%), (ii) the host overshines the AGN, or (iii) the AGN is obscured and not immediately recognizable by line diagnostics. Next-generation X-ray missions, such as AXIS, could then play a fundamental role in uncovering this, thus far, hidden population of sources. However, we caution the reader that this topic is a matter of fierce debate in the literature (see, e.g., [37]), particularly while methods for measuring SMBH masses at high-z are being refined.
This approach has been successfully implemented by [38] at redshift z∼0, and future studies extending this methodology to higher redshifts will be crucial for obtaining a complete understanding of SMBH growth and its coevolution with host galaxies. By investigating the interaction between the growth of large-scale structures and the accretion rates on SMBHs, profound advances can be made in our understanding of the AGN–galaxy coevolution (e.g., [39,40]).
By examining the intricate relationship between SMBH growth, star formation, and structure formation across different cosmic epochs, we can gain transformative insights into the processes that govern the coevolution of AGN and their host galaxies. This comprehensive understanding will contribute to our broader knowledge of galaxy evolution and the formation of the cosmic structures that we observe today.

3. AXIS Pillar 1: Determining the Nature of SMBH Seeds with AXIS at z = 6–10 X-ray Luminosity Function (XLF)

The growth of SMBHs by accretion and mergers gradually erases seeding information. However, the closer to the seeding epoch that we can probe AGN, the closer we get to their initial masses and host galaxy relationships. The epoch of z 7 is a special time: only at these extreme redshifts do quasars require seeds to have grown continuously near the Eddington limit to acquire their observed masses (e.g., [3,4]). This places interesting joint constraints on the astrophysics of SMBH seeding and accretion, which grow more and more informative as redshift increases.
Using a semi-analytical model (SAM) that integrated the SMBH assembly from z = 20 to z = 0 , Ricarte and Natarajan [33] simultaneously studied gravitational wave signals, luminosity functions, and BH occupation fractions (the fraction of galaxies that host massive BHs) to search for signatures that could be used to discriminate seeding models. This model contrasted rare and “heavy” seeds from direct collapse with abundant and low-mass light seeds from the remnants of Population III stars. Using empirical relations, these seeds grow at a rate capped at the Eddington limit with a duty cycle modulated by major halo mergers. Interestingly, the luminosity functions of these two cases are identical for z 8 and down to L X 1042 erg s−1, where AGN behavior was determined almost entirely by SMBH–galaxy coevolution. However, the two cases diverge at higher redshift, nearer to the seeding epoch, where the Eddington limit leads to very different high-to-moderate luminosity functions. We note that there are three major factors that determine the luminosity function at this epoch: BH masses, occupation fractions, and accretion rates. In combination with other facilities, AXIS will enable us to disentangle these effects.
In Figure 4, we plot the source counts predicted from [33] SAM (heavy seeds in red, light seeds in blue). We add in orange a compilation of LF obtained from the large-scale cosmological hydrodynamical simulations of the field [41]. The scatter at fixed Lx reflects the impact of different seeding, BH growth, and feedback modeling in these simulations1. The points with error bars represent the combination of the proposed “Deep” and “Intermediate” survey strategies, assuming the SAM points as true. SAM bolometric luminosity functions have been converted to 2–10 keV assuming bolometric corrections of Duras et al. [48] for Type II AGN and have been additionally suppressed by an ad hoc factor of 10, assuming an extreme Compton-thick fraction. Despite all being calibrated to eventually yield the correct local scaling relations between SMBH mass and galaxy mass, the cosmological simulations span a wide range, implying that additional constraints at this epoch would help constrain the astrophysics of this epoch. As discussed above, the two SAM curves begin to diverge at z 8 due to the assumed Eddington limit. Since heavy seeds were initialized at higher masses, they are able to produce AGN with higher luminosity at high redshifts than their light seed counterparts. A measurement of the evolution of the luminosity function during this epoch would place interesting joint constraints on seeding and accretion. The main goal of AXIS will be to unambiguously detect SMBH at z > 7 to derive the XLF to inform and constrain seeding models. This will only be possible by performing an X-ray survey with AXIS in the region where NIR measurements with the JWST and/or Roman exist or follow-ups are possible to properly characterize the host galaxies.
So far, extrapolations from population synthesis models differ by up to two orders of magnitude from theoretical models, but even the deepest Chandra surveys are not able to discriminate among them at z 6 . In fact, regardless of the assumed model, < < 1 AGN is expected to be detected with Chandra at z > 6 with flux < 10 17 erg s−1 cm−2, since only 10 arcmin 2 in the 7 Ms Chandra Deep Field South are sensitive to such emission level.

Accretion Mode and Signatures of the First SMBHs

To date, we detected more than 300 quasars at z > 5.9 , with tens of them already characterized by mass estimates [14]. Several confirmed detections at z > 7 , including some recent striking observations of candidate AGN by the JWST at z∼6–12 [5,15,49], provide us with a test bed for our models of black hole growth.
Seeds, light and heavy, are predicted to form in the redshift range z 20–30 [13]. Because the detection redshifts of quasars are increasing, the time between seed formation and SMBH detection is decreasing (see, e.g., [14,50]). The decreasing time between seeding and detection shrinks the parameter space of seed properties [51]: initial mass, average Eddington ratio, duty cycle, and radiative efficiency. Consequently, detecting SMBHs at a very high redshift allows us to pinpoint the properties of black hole seeds in the most efficient way. For example, detecting an SMBH with mass ∼ 10 10 M by z 9 would exclude the entire parameter space available for light seeds and dramatically reduce the one for heavy seeds [51]. The detection of a more moderate case (e.g., GOODS-North-z11 GNz11, if confirmed, a 10 7 M at z 10 ) would still provide important constraints on the seeding mechanisms, especially if the (instantaneous) Eddington ratio is constrained. For example, Pacucci and Loeb [51] shows that, in this case, the marginal distributions for the seed mass for eight quasars already detected at z > 7 can already be constrained down to two orders of magnitude in mass. The modeling of some of these quasars already shows a preference for heavy-seeding scenarios (see, e.g., [4,52]).
Detecting accreting SMBHs at z > 7–8 is one of the main goals of AXIS. With its surveys, AXIS will be able to detect low-luminosity AGN up to the highest redshifts probed spectroscopically with the JWST (see Figure 1). By detecting X-ray emissions from these sources, a first estimate of their black hole mass will be available, assuming accretion at the Eddington rate. This first analysis can be complemented by follow-up observations with the JWST, either with photometry or, for sources that are sufficiently bright, with spectroscopy. This additional information will constrain the mass and the redshift of the source. Hence, improved modeling with those as priors will further pinpoint the initial seed mass and the average growth rate (i.e., sub-Eddington vs. super-Eddington) to a precision that is unreachable today with current facilities.
In the electromagnetic realm, signatures of heavy vs. light seeding can be grouped into two categories: seed-related and host-related. Seed-related signatures are generated from the black hole itself. Early studies (see, e.g., [33,53,54]) show that heavy seeding will lead to steep, red infrared spectra and a bell-shaped X-ray emission peaked at ∼1 keV rest-frame. While current X-ray facilities, such as Chandra, are able to detect extremely massive, unobscured seeds, next-generation facilities are urgently needed to detect broader, fainter, and obscured populations. With a peak emission, at ∼1 keV, 1–2 orders of magnitude fainter, obscured populations of seeds will remain unobserved until next-generation instruments are operative [55]. AXIS will thus play a crucial role in detecting the Compton-thick population of heavy black hole seeds. The obscured fraction of seeds is expected to represent the norm at such high redshifts [56].
Host-related signatures, such as the black-hole-to-stellar-mass ratio, are essential in discriminating between light- and heavy-seeding models. In the heavy seed case, the ratio is expected to remain of order unity for several Myrs after the seeding process (see, e.g., [57]). On the contrary, in the light seeding case, the ratio is expected to be much lower and progressively stabilize around the local ratio of ∼10−4 over some Gyrs. Assuming typical effective radii nowadays detected for z 10 hosts, the JWST will be able to detect their infrared light and constrain their stellar mass [58,59]. As X-ray observations are crucial both for studying the Spectral Energy Distribution of the X-ray emission and for estimating the mass of the SMBH, AXIS will become the facility of reference to investigate both populations of black hole seeds.
The X-ray spectra of AGN are also a direct tool for studying the accretion mechanism and the properties of the accretion discs. In the redshift range 0 < z < 6 , Vito et al. [24] reports a lack of substantial evolution of the inner accretion-disk and hot-corona structure in QSOs. Surprisingly, recent results from [25] show that a sample of very luminous broad-line AGN z > 6 presents an unusually steep X-ray spectrum with spectral index Γ ∼2.4 compared to average Γ 1.9 observed at z < 6 . This points to a different physics of the disc and the corona of the AGN at high-z which, combined with the peculiar requirements of fast growth of the SMBH at high-z, makes the study of the physics of accretion at that epoch even more intriguing. AXIS will be able to study the spectra of hundreds of AGN at z > 6, allowing the study of accretion physics over a much wider range of luminosities and accretion rates. Another suggestion of [25], which cannot be currently tested because of the low count statistics, is that these sources could instead present an unusual extreme low-energy cutoff and a typical Γ = 1.9 spectral index. AXIS with its very large effective area will be able to produce much better spectra than XMM-Newton at high-z and further test the accretion physics of AGN in the early Universe.

4. Assembly of Black Holes with Large-Scale Structures as Traced by High-z AGN Clustering

A complementary way to study the growth of SMBHs is by connecting them to their larger-scale dark matter environments, probed by the spatial clustering of AGN. Together with the luminosity function, AGN clustering constrains the relationship between black holes and their host dark matter halos and infers any preferred cosmic environments for AGN activity. Since the growth and evolution of halos are well understood, clustering can be used to statistically connect SMBHs to their progenitors in the earlier Universe, which would infer how certain populations evolve over time.
The area and depth of the intermediate AXIS survey (with redshift measurements from future multiwavelength facilities; see Section 5.2.2) would allow the uncertainties of AGN clustering at z 3 to be at the 10% level or better on both the large and small scales (0.1–10 Mpc/h). This is an order of magnitude better than the NewAthena ESA concept mission, and unfeasible with current instruments. Cross-correlations with galaxies at intermediate redshifts ( z 1.5 ) are predicted to provide uncertainties less than 5%. The unprecedented clustering statistics at these epochs (especially on 1-halo scales; <1 Mpc/h) would allow for measurements in several bins of AGN luminosity, which would provide powerful constraints on AGN fueling and feedback mechanisms at the peak of SMBH growth (e.g., determining whether or not the majority of AGN are triggered by galaxy mergers, as well as constraining the environmental dependence of the distribution of Eddington ratios).
Finally, the JWST will dramatically improve the characterization of the close environments around high-redshift ( z > 6 ) AGN via observations of galaxy fields around AXIS AGN of varying luminosities. This will determine whether or not high-z AGN reside in overdensities for a wide range of accretion rates and mass scales. Although theory predicts that the most massive black holes at high redshift reside in the progenitors of today’s massive clusters, observations have so far led to conflicting conclusions, most likely due to the limited HST sensitivity. The JWST will detect high-z galaxies much more efficiently, and together with the wide range of AGN luminosities from the AXIS surveys, this long-standing problem will be solved.
Moreover, the clustering statistics of over- (or under-)massive SMBHs at high-z will shed light on their seeding mechanisms, as those that originated from DCBHs are predicted to be significantly more clustered than SMBHs formed via other mechanisms. This is because the Lyman–Werner radiation from one DMBH would likely trigger additional DCBHs within the same region, resulting in highly clustered AGN statistics in the early Universe. Therefore, the clustering amplitudes of AXIS AGN will provide powerful additional constraints on such SMBH formation models.

5. AXIS Survey of the First AGN

Surveys with the AXIS probe were first simulated and discussed in [60], to which we refer for a more detailed explanation of the simulations. Both the Marchesi et al. [60] simulations and the ones used to obtain the numbers reported in this paper were performed using the Monte Carlo code simulation of X-ray telescopes (hereafter SIXTE [61]). This software allows one to simulate an observation with an X-ray telescope by creating a photon list, which includes the arrival time, energy, and position of each photon based on the simulated telescope setup (i.e., effective area, field of view, PSF, vignetting, read-out properties, and redistribution matrix). In doing so, SIXTE generates an impact list that contains the energy and arrival time of each photon, as well as its position on the detector. The final event file is obtained from this intermediate list, reprocessed to take into account the simulated detector read-out properties and redistribution matrix file.
The input catalog we used for our simulations is based on the [28] AGN population synthesis model. AGN were simulated down to a 0.5–2 keV luminosity L0.5–2 = 1040 erg s−1 and up to redshift z = 3. The mocks were checked to ensure the correct reproduction of the trends with luminosity, redshift, and column density of the AGN densities as a function of luminosity, redshift, and column density and are in close agreement with the observational evidence at all fluxes sampled by current X-ray surveys. In the high-redshift regime (i.e., at z > 3, where the AGN space density starts declining), we use a separate mock catalog, built from the [22] z > 3 AGN luminosity function, which nicely describes the observational evidence from the deepest X-ray surveys currently available up to redshift z 5 (e.g., [62]). Finally, neither AGN nor host galaxy clustering are included in the mock generation. The catalogs we used are available online at http://cxb.oas.inaf.it/mock.html (accessed on 15 July 2023) in FITS format and are ready to be used within SIXTE.

5.1. The AXIS Wedding Cake Survey

Following a standard approach extensively used in past X-ray surveys, AXIS plans to use a so-called “Wedding cake” strategy (Figure 5) to perform its extragalactic surveys as follows:
  • A deep, pencil-beam survey; that is, a 7 Ms observation of a single AXIS pointing (∼0.13 deg2).
  • An intermediate-area and intermediate-depth survey, which would cover 2 deg2 with a uniform, 360 ks exposure, for an overall time request of 6 Ms.
  • Finally, while no wide-area survey is currently planned, it would be possible to obtain an AXIS Serendipitous field by combining Guest Observer observations. Assuming 20 Ms of GO (Guest Observer) non-galactic plane time, with a median of 50 ks per pointing, it would be possible to cover 50 deg2 with a sensitivity ∼10−16 erg s−1 cm−2.
We report in Table 1 a summary of the properties of these surveys: the number of detections and flux limits were obtained from end-to-end simulations with the SIXTE tool described in the previous section. Note that in this table, we report only the expected number of sources at z < 6, as the model used to generate the mocks employs a simple power-law extrapolation of z∼3–4 XLF. In the AXIS Deep Field, the 3�� σ upper boundary of [22] XLF predicts that one should detect ∼25 AGN in the redshift range z = [6–7]; in the Intermediate Field, instead, in the same redshift range, one should detect ∼80 sources. In the same redshift range, the SAM heavy seed scenario predicts ∼120 detections in the Deep Field and ∼240 in the Intermediate Field, while the SAM light seed scenario predicts 120 detections in the Deep Field and ∼230 in the Intermediate one. At higher redshifts, as shown in Figure 4, the two SAM models diverge more significantly, and the expected number of AXIS detections is much larger. For example, at z = 9, the SAM heavy seed scenario predicts 10–12 detections in the two AXIS fields, while the SAM light seed scenario predicts only 1–2 detections.

5.2. Synergies

AXIS is proposed to fly during an unprecedented epoch for extragalactic surveys. In the early 2030s, NASA and ESA will boast an exceptional fleet of survey satellites, while powerful ground-based telescopes will observe the first light of the Universe. The JWST, having completed its first decade of operation, will have already amassed a wealth of near-infrared to mid-infrared photometric and spectroscopic data from various deep survey fields. Roman and Euclid, covering a significant portion of the extragalactic sky in the visible and near-infrared bands, will contribute to the comprehensive view of the Universe. Meanwhile, Legacy Survey of Space and Time (LSST)–Rubin [63] will observe the extragalactic sky with high cadence, capturing transient and variable phenomena, while the ELT (Extremely Large Telescope, https://www.eso.org/public/teles-instr/elt/, accessed on 15 July 2023) will make its debut as a cutting-edge 30-meter class telescope. Additionally, the SKA (Square Kilometre Array, [64]) will provide a groundbreaking perspective of the radio sky.
Within this vibrant observational landscape, AXIS will serve as the ideal X-ray counterpart, completing the panchromatic view of the extragalactic Universe. By harnessing its capabilities in X-ray observations, AXIS will play a vital role in unequivocally detecting accreting objects, such as AGN and other X-ray-emitting sources. Its synergistic partnership with the aforementioned observatories and telescopes will allow for a comprehensive understanding of the properties, dynamics, and coevolution of objects across the electromagnetic spectrum. This collaborative effort will significantly advance our understanding of the extragalactic Universe and provide a holistic view of the diverse phenomena that shape our cosmic landscape.

5.2.1. Optical Identifications

While X-ray emission serves as a clear indicator of accretion in galaxies, mere detection is not sufficient to obtain informative insights. However, due to the high angular resolution of AXIS, we will have the ability to achieve optical identification with follow-up observations for more than 95% of the detected sources, as estimated by the COSMOS-Legacy survey [65]. To accomplish this, we will rely on NIR source catalogs [66] from observatories such as HST3, Subaru [67], the JWST4, Roman5, Euclid6, and potentially the ELT.
Optical identification involves more than just positional cross-matching, as obtaining precise source centroids is crucial. It also requires leveraging prior knowledge of the colors and magnitude distribution of potential counterparts. To facilitate these complex optical identifications, various tools have been developed, including LYR and N-WAY. These tools use well-tested Bayesian approaches to efficiently perform the intricate process of optical identification [20,68,69].
By combining the X-ray detections from AXIS with accurate optical identifications, we can establish robust associations between the X-ray sources and their optical counterparts. This multiwavelength approach not only provides valuable information about the nature of the accreting objects but also enables detailed studies of their properties, including redshifts, spectral features, and host galaxy characteristics. The combination of X-ray and optical data, supported by the advanced identification techniques mentioned above, will significantly enhance our understanding of the connection between accretion processes and the larger astrophysical context.

5.2.2. Photometric and Spectroscopic Redshifts

To accurately measure the luminosity function of high-redshift AGN, obtaining their redshift information is of fundamental importance. The survey fields targeted by AXIS will be strategically chosen in highly investigated areas of the sky, including regions such as COSMOS, Chandra Deep Fields, the JWST Deep Fields, and areas covered by either Roman or Euclid or both.
The counterparts of high-redshift AGN are expected to be extremely faint, with magnitudes on the order of m A B 3.6 μ m 26–29 (e.g., [70]). Preliminary results from the JWST indicate that the NIR magnitudes of the first AGN candidates identified at redshifts z∼7–10 fall within the range of 25–26 [5]. By cross-matching these sources with AXIS data, they can be easily identified, especially considering that the typical flux limit of deep JWST surveys is around 29 (AB mag). In cases where optical identifications are available, spectroscopic redshifts can be obtained either from previous JWST measurements or through follow-up observations using JWST-NIRSPEC for fainter sources. Spectroscopic campaigns will also be initiated using 10m and 30m class telescopes such as Keck MOSFIRE7, Subaru PFS8, or the ELT. Additionally, grism spectroscopic campaigns conducted by Roman and Euclid will provide valuable sources of redshifts for the brightest sources with redshifts z < 6–8.
Spectra play a crucial role in determining the masses of black holes, which is essential for understanding the distributions of Eddington rates and tracking the evolution of the MBH σ relation. By obtaining spectroscopic information, AXIS, in conjunction with other observatories and telescopes, will contribute to advancing our knowledge of the masses of black holes in high-redshift AGN. This, in turn, will enable detailed investigations into the accretion processes and the relationships between black hole masses and the properties of their host galaxies.
To assemble an even larger sample of redshift measurements, photometric redshifts (Photo-z) will be relied upon. Photo-z has become a well-established technique to determine the redshifts of AGN [30,71]. The abundant amount of ancillary photometric data available from NIR surveys will enable precise measurements of Photo-z. Using this photometric information, we can obtain high-quality estimates of redshifts for almost all sources detected in the survey.
The Photo-z measurements can be effectively employed in determining the XLF. The sources can be distributed across the redshift range based on their associated probability density functions, allowing for a comprehensive determination of the XLF (e.g., [30]). This approach enables us to derive statistical information about the AGN population and its evolution over cosmic time.
As a valuable by-product of this work, we will obtain measurements of the Spectral Energy Distribution (SED) for all sources. The SEDs provide crucial insights into the host properties of AGN, including stellar mass, gas mass, extinction, and indirectly the black hole mass and Eddington rates. By analyzing the SEDs, we can gain a deeper understanding of the physical characteristics and compositions of the host galaxies, shedding light on the intricate interplay between the central supermassive black hole and its environment.
Through the utilization of photometric redshifts and the comprehensive analysis of SEDs, AXIS will significantly contribute to the study of AGN populations, their properties, and their coevolution with host galaxies. This approach will provide valuable information to advance our understanding of the growth and dynamics of SMBHs across cosmic time.

5.3. Stacking Analysis of JWST Detected High-z Galaxies

The strategic stacking of X-ray images at locations of known galaxies offers a remarkable increase in sensitivity. This approach provides a powerful means to investigate the average properties of a population that may otherwise be challenging to study individually. An illustrative example of the efficacy of this technique was demonstrated by [62], who stacked X-ray signals from HST-selected massive galaxies at redshifts of 4–5 using the 7 Ms CDFS (Chandra Deep Field South). Using this method, they achieved effective exposure times of approximately 109 s, exceeding Chandra’s capabilities, and successfully detected the cumulative X-ray signal up to redshifts around five.
Considering a short AXIS exposure of 50 ks in regions covered by the JWST, assuming a magnitude limit of 28th at 7 μm, it is estimated that the expected source density would be on the order of 20,000 deg−2. If 1% of these sources are AGN, statistically significant fluxes on the order of 5 × 10−19 erg s−1 cm−2 can be achieved. At a redshift of approximately four, this corresponds to luminosities of around 1041 erg s−1. In the Deep Fields, AXIS will reach even deeper flux levels, although the scaling with exposure time is not linear because of the transition from being photon-limited to being background-limited at such depths. This capability will enable AXIS to study populations of very faint sources detected by the JWST, extending investigations to redshifts well beyond 10.
By employing the stacking technique and leveraging the synergistic observations with the JWST, AXIS will significantly contribute to our understanding of the high-redshift Universe and the properties of faint AGN populations. This combined effort will push the boundaries of knowledge, enabling investigations into the early stages of galaxy and black hole formation and evolution.

6. Conclusions

In this white paper, we thoroughly explored the necessity and potential of AXIS in investigating the early accretion of SMBHs in the Universe. Specifically:
  • AXIS will measure the XLF of AGN up to redshifts around 10, providing crucial insights into their evolution over cosmic time with its deep and intermediate surveys in combination with ancillary information from the JWST, Roman, Euclid, LSST–Rubin, and the ELT. With the XLF, AXIS will inform and constrain models of SMBH seeding, potentially helping to distinguish between scenarios where they grow from light or heavy seeds.
  • AXIS will significantly expand the current flux–area parameter space of X-ray surveys, covering a range of 0.13–2.5 deg2 and reaching much fainter fluxes. This, coupled with the constant PSF across the field of view, will result in a leap forward in survey grasp, achieving a two-order-of-magnitude improvement compared to Chandra.
  • AXIS will play a crucial role in determining the nature of high-redshift galaxies observed by the JWST that exhibit characteristics typical of star-forming galaxies (SFGs), aiding in their classification as AGN in particular in the case of highly obscured and low-metallicity sources.
  • The excellent angular resolution of AXIS will facilitate the clear identification of X-ray sources with NIR detections, enabling comprehensive studies of their Spectral Energy Distributions (SEDs) and spectral features.
  • With access to a wealth of ancillary data, AXIS will enable the study of the onset and evolution of galaxy–SMBH coevolution.
  • The high source surface density observed by AXIS will allow for the characterization of the coevolution of early AGN with their environments through clustering studies. It will also enable investigations into the environment of individual high-redshift quasars.
  • Leveraging the low background and exceptional angular resolution, stacking AXIS images at the positions of very faint JWST sources will facilitate the study and constraint of the AGN population, potentially extending to redshifts beyond 10 and probing low-luminosity (black hole mass) sources at lower redshifts.
This white paper, along with the entire series, emphasizes the crucial and unparalleled importance of adopting a high-angular-resolution X-ray telescope such as AXIS as a key instrument for studying the Universe in the early 2030s. AXIS, a remarkable X-ray survey machine building upon the legacy of Chandra, will surpass the achievements of its predecessor and push the boundaries of black hole exploration into new and uncharted territories.

Author Contributions

All the authors contributed equally to the manuscript that was assembled and written by N.C. All authors have read and agreed to the published version of the manuscript.

Funding

This project was supported by the University of Miami.

Data Availability Statement

This work is based on existing data available on public archives.

Acknowledgments

We kindly acknowledge the entire AXIS team for their outstanding scientific and technical work. This work is the result of several months of discussion in the AXIS-AGN SWG. N.C. thanks the University of Miami for its support during the proposal and writing phase. This white paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS white papers can be found at the AXIS website with a mission overview here.

Conflicts of Interest

The authors declare no conflicts of interest.

Notes

1
Values are collected in Habouzit et al. in prep. from Horizon-AGN [42], Illustris [43], EAGLE [44], TNG100 and TNG300 [45], SIMBA [46], and Astrid [47].
2
The wide Serendipitous is not an official AXIS program but it is an estimate of the extragalatic survey sensitivity obtained by grouping all the expected GO observations over the 5 year of mission lifetime.
3
https://www.stsci.edu/hst, accessed on 15 July 2023
4
https://www.stsci.edu/JWST, accessed on 15 July 2023
5
https://roman.gsfc.nasa.gov, accessed on 15 July 2023
6
7
8
https://pfs.ipmu.jp, accessed on 15 July 2023

References

  1. 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]
  2. Bañados, E.; Venemans, B.P.; Morganson, E.; Decarli, R.; Walter, F.; Chambers, K.C.; Rix, H.W.; Farina, E.P.; Fan, X.; Jiang, L.; et al. Discovery of Eight z∼6 Quasars from Pan-STARRS1. Astron. J. 2014, 148, 14. [Google Scholar] [CrossRef]
  3. Wang, F.; Yang, J.; Fan, X.; Hennawi, J.F.; Barth, A.J.; Banados, E.; Bian, F.; Boutsia, K.; Connor, T.; Davies, F.B.; et al. A Luminous Quasar at Redshift 7.642. Astrophys. J. Lett. 2021, 907, L1. [Google Scholar] [CrossRef]
  4. Bogdan, A.; Goulding, A.; Natarajan, P.; Kovacs, O.; Tremblay, G.; Chadayammuri, U.; Volonteri, M.; Kraft, R.; Forman, W.; Jones, C.; et al. Detection of an X-ray quasar in a gravitationally-lensed z = 10.3 galaxy suggests that early supermassive black holes originate from heavy seeds. arXiv 2023, arXiv:2305.15458. [Google Scholar] [CrossRef]
  5. Maiolino, R.; Scholtz, J.; Witstok, J.; Carniani, S.; D’Eugenio, F.; de Graaff, A.; Uebler, H.; Tacchella, S.; Curtis-Lake, E.; Arribas, S.; et al. A small and vigorous black hole in the early Universe. arXiv 2023, arXiv:2305.12492. [Google Scholar] [CrossRef] [PubMed]
  6. Lodato, G.; Natarajan, P. Supermassive black hole formation during the assembly of pre-galactic discs. Mon. Not. R. Astron. Soc. 2006, 371, 1813–1823. [Google Scholar] [CrossRef]
  7. Begelman, M.C.; Volonteri, M.; Rees, M.J. Formation of supermassive black holes by direct collapse in pre-galactic haloes. Mon. Not. R. Astron. Soc. 2006, 370, 289–298. [Google Scholar] [CrossRef]
  8. Volonteri, M. Formation of supermassive black holes. Astron. Astrophys. Rev. 2010, 18, 279–315. [Google Scholar] [CrossRef]
  9. Hawking, S. Gravitationally collapsed objects of very low mass. Mon. Not. R. Astron. Soc. 1971, 152, 75. [Google Scholar] [CrossRef]
  10. Carr, B.; Kohri, K.; Sendouda, Y.; Yokoyama, J. Constraints on primordial black holes. Rep. Prog. Phys. 2021, 84, 116902. [Google Scholar] [CrossRef]
  11. Cappelluti, N.; Hasinger, G.; Natarajan, P. Exploring the High-redshift PBH-ΛCDM Universe: Early Black Hole Seeding, the First Stars and Cosmic Radiation Backgrounds. Astrophys. J. 2022, 926, 205. [Google Scholar] [CrossRef]
  12. García-Bellido, J. Massive Primordial Black Holes as Dark Matter and their detection with Gravitational Waves. J. Phys. Conf. Ser. 2017, 840, 012032. [Google Scholar] [CrossRef]
  13. Barkana, R.; Loeb, A. In the beginning: The first sources of light and the reionization of the universe. Phys. Rep. 2001, 349, 125–238. [Google Scholar] [CrossRef]
  14. Fan, X.; Banados, E.; Simcoe, R.A. Quasars and the Intergalactic Medium at Cosmic Dawn. arXiv 2022, arXiv:2212.06907. [Google Scholar] [CrossRef]
  15. Larson, R.L.; Finkelstein, S.L.; Kocevski, D.D.; Hutchison, T.A.; Trump, J.R.; Arrabal Haro, P.; Bromm, V.; Cleri, N.J.; Dickinson, M.; Fujimoto, S.; et al. A CEERS Discovery of an Accreting Supermassive Black Hole 570 Myr after the Big Bang: Identifying a Progenitor of Massive z > 6 Quasars. arXiv 2023, arXiv:2303.08918. [Google Scholar] [CrossRef]
  16. Luo, B.; Brandt, W.N.; Xue, Y.Q.; Lehmer, B.; Alexander, D.M.; Bauer, F.E.; Vito, F.; Yang, G.; Basu-Zych, A.R.; Comastri, A.; et al. The Chandra Deep Field-South Survey: 7 Ms Source Catalogs. Astrophys. J. 2017, 228, 2. [Google Scholar] [CrossRef]
  17. Reynolds, C.S.; Kara, E.A.; Mushotzky, R.F.; Ptak, A.; Koss, M.J.; Williams, B.J.; Allen, S.W.; Bauer, F.E.; Bautz, M.; Bodaghee, A.; et al. Overview of the Advanced X-ray Imaging Satellite (AXIS). arXiv 2023, arXiv:2311.00780. [Google Scholar] [CrossRef]
  18. Martinez Galarza, J. The Chandra Source Catalog version 2.1: New Avenues for Discovery in X-ray Datasets. AAS High Energy Astrophys. Div. 2023, 20, 404.01. [Google Scholar]
  19. Webb, N.A.; Coriat, M.; Traulsen, I.; Ballet, J.; Motch, C.; Carrera, F.J.; Koliopanos, F.; Authier, J.; de la Calle, I.; Ceballos, M.T.; et al. The XMM-Newton serendipitous survey. IX. The fourth XMM-Newton serendipitous source catalogue. Astron. Astrophys. 2020, 641, A136. [Google Scholar] [CrossRef]
  20. Peca, A.; Cappelluti, N.; Urry, C.M.; LaMassa, S.; Marchesi, S.; Ananna, T.T.; Baloković, M.; Sanders, D.; Auge, C.; Treister, E.; et al. On the Cosmic Evolution of AGN Obscuration and the X-Ray Luminosity Function: XMM-Newton and Chandra Spectral Analysis of the 31.3 deg2 Stripe 82X. Astrophys. J. 2023, 943, 162. [Google Scholar] [CrossRef]
  21. Yang, G.; Caputi, K.I.; Papovich, C.; Arrabal Haro, P.; Bagley, M.B.; Behroozi, P.; Bell, E.F.; Bisigello, L.; Buat, V.; Burgarella, D.; et al. CEERS Key Paper. VI. JWST/MIRI Uncovers a Large Population of Obscured AGN at High Redshifts. Astrophys. J. 2023, 950, L5. [Google Scholar] [CrossRef]
  22. Vito, F.; Gilli, R.; Vignali, C.; Comastri, A.; Brusa, M.; Cappelluti, N.; Iwasawa, K. The hard X-ray luminosity function of high-redshift (3 < z ≲ 5) active galactic nuclei. Mon. Not. R. Astron. Soc. 2014, 445, 3557–3574. [Google Scholar] [CrossRef]
  23. Lyu, J.; Alberts, S.; Rieke, G.H.; Shivaei, I.; Pérez-González, P.G.; Sun, F.; Hainline, K.N.; Baum, S.; Bonaventura, N.; Bunker, A.J.; et al. Active Galactic Nuclei Selection and Demographics: A New Age with JWST/MIRI. Astrophys. J. 2024, 966, 229. [Google Scholar] [CrossRef]
  24. Vito, F.; Brandt, W.N.; Bauer, F.E.; Calura, F.; Gilli, R.; Luo, B.; Shemmer, O.; Vignali, C.; Zamorani, G.; Brusa, M.; et al. The X-ray properties of z > 6 quasars: No evident evolution of accretion physics in the first Gyr of the Universe. Astron. Astrophys. 2019, 630, A118. [Google Scholar] [CrossRef]
  25. Zappacosta, L.; Piconcelli, E.; Fiore, F.; Saccheo, I.; Valiante, R.; Vignali, C.; Vito, F.; Volonteri, M.; Bischetti, M.; Comastri, A.; et al. HYPerluminous quasars at the Epoch of ReionizatION (HYPERION). A new regime for the X-ray nuclear properties of the first quasars. arXiv 2023, arXiv:2305.02347. [Google Scholar] [CrossRef]
  26. Maiolino, R.; Risaliti, G.; Signorini, M.; Trefoloni, B.; Juodzbalis, I.; Scholtz, J.; Uebler, H.; D’Eugenio, F.; Carniani, S.; Fabian, A.; et al. JWST meets Chandra: A large population of Compton thick, feedback-free, and X-ray weak AGN, with a sprinkle of SNe. arXiv 2024, arXiv:2405.00504. [Google Scholar] [CrossRef]
  27. Natarajan, P.; Pacucci, F.; Ferrara, A.; Agarwal, B.; Ricarte, A.; Zackrisson, E.; Cappelluti, N. Unveiling the First Black Holes With JWST:Multi-wavelength Spectral Predictions. Astrophys. J. 2017, 838, 117. [Google Scholar] [CrossRef]
  28. Gilli, R.; Comastri, A.; Hasinger, G. The synthesis of the cosmic X-ray background in the Chandra and XMM-Newton era. Astron. Astrophys. 2007, 463, 79–96. [Google Scholar] [CrossRef]
  29. Treister, E.; Urry, C.M.; Virani, S. The Space Density of Compton-Thick Active Galactic Nucleus and the X-Ray Background. Astrophys. J. 2009, 696, 110–120. [Google Scholar] [CrossRef]
  30. Ananna, T.T.; Salvato, M.; LaMassa, S.; Urry, C.M.; Cappelluti, N.; Cardamone, C.; Civano, F.; Farrah, D.; Gilfanov, M.; Glikman, E.; et al. AGN Populations in Large-volume X-Ray Surveys: Photometric Redshifts and Population Types Found in the Stripe 82X Survey. Astrophys. J. 2017, 850, 66. [Google Scholar] [CrossRef]
  31. Terao, K.; Nagao, T.; Onishi, K.; Matsuoka, K.; Akiyama, M.; Matsuoka, Y.; Yamashita, T. Multiline Assessment of Narrow-line Regions in z 3 Radio Galaxies. Astrophys. J. 2022, 929, 51. [Google Scholar] [CrossRef]
  32. Harikane, Y.; Zhang, Y.; Nakajima, K.; Ouchi, M.; Isobe, Y.; Ono, Y.; Hatano, S.; Xu, Y.; Umeda, H. JWST/NIRSpec First Census of Broad-Line AGNs at z = 4–7: Detection of 10 Faint AGNs with M_BH~106–107 M_sun and Their Host Galaxy Properties. arXiv 2023, arXiv:2303.11946. [Google Scholar] [CrossRef]
  33. Ricarte, A.; Natarajan, P. The observational signatures of supermassive black hole seeds. Mon. Not. R. Astron. Soc. 2018, 481, 3278–3292. [Google Scholar] [CrossRef]
  34. Baloković, M.; Comastri, A.; Harrison, F.A.; Alexander, D.M.; Ballantyne, D.R.; Bauer, F.E.; Boggs, S.E.; Brandt, W.N.; Brightman, M.; Christensen, F.E.; et al. The NuSTAR View of Nearby Compton-thick Active Galactic Nuclei: The Cases of NGC 424, NGC 1320, and IC 2560. Astrophys. J. 2014, 794, 111. [Google Scholar] [CrossRef]
  35. Hickox, R.C.; Mullaney, J.R.; Alexander, D.M.; Chen, C.T.J.; Civano, F.M.; Goulding, A.D.; Hainline, K.N. Black Hole Variability and the Star Formation-Active Galactic Nucleus Connection: Do All Star-forming Galaxies Host an Active Galactic Nucleus? Astrophys. J. 2014, 782, 9. [Google Scholar] [CrossRef]
  36. Pacucci, F.; Nguyen, B.; Carniani, S.; Maiolino, R.; Fan, X. JWST CEERS and JADES Active Galaxies at z = 4–7 Violate the Local M-M🟉 Relation at >3σ: Implications for Low-mass Black Holes and Seeding Models. Astrophys. J. 2023, 957, L3. [Google Scholar] [CrossRef]
  37. Li, J.; Silverman, J.D.; Shen, Y.; Volonteri, M.; Jahnke, K.; Zhuang, M.Y.; Scoggins, M.T.; Ding, X.; Harikane, Y.; Onoue, M.; et al. Tip of the iceberg: Overmassive black holes at 4 < z < 7 found by JWST are not inconsistent with the local MBH-M🟉 relation. arXiv 2024, arXiv:2403.00074. [Google Scholar] [CrossRef]
  38. Ananna, T.T.; Urry, C.M.; Ricci, C.; Natarajan, P.; Hickox, R.C.; Trakhtenbrot, B.; Treister, E.; Weigel, A.K.; Ueda, Y.; Koss, M.J.; et al. Probing the Structure and Evolution of BASS Active Galactic Nuclei through Eddington Ratios. Astrophys. J. 2022, 939, L13. [Google Scholar] [CrossRef]
  39. Powell, M.C.; Cappelluti, N.; Urry, C.M.; Koss, M.; Finoguenov, A.; Ricci, C.; Trakhtenbrot, B.; Allevato, V.; Ajello, M.; Oh, K.; et al. The Swift/BAT AGN Spectroscopic Survey. IX. The Clustering Environments of an Unbiased Sample of Local AGNs. Astrophys. J. 2018, 858, 110. [Google Scholar] [CrossRef]
  40. Powell, M.C.; Urry, C.M.; Cappelluti, N.; Johnson, J.T.; LaMassa, S.M.; Ananna, T.T.; Kollmann, K.E. The Clustering of X-Ray Luminous Quasars. Astrophys. J. 2020, 891, 41. [Google Scholar] [CrossRef]
  41. Habouzit, M.; Somerville, R.S.; Li, Y.; Genel, S.; Aird, J.; Anglés-Alcázar, D.; Davé, R.; Georgiev, I.Y.; McAlpine, S.; Rosas-Guevara, Y.; et al. Supermassive black holes in cosmological simulations - II: The AGN population and predictions for upcoming X-ray missions. Mon. Not. R. Astron. Soc. 2022, 509, 3015–3042. [Google Scholar] [CrossRef]
  42. Dubois, Y.; Peirani, S.; Pichon, C.; Devriendt, J.; Gavazzi, R.; Welker, C.; Volonteri, M. The HORIZON-AGN simulation: Morphological diversity of galaxies promoted by AGN feedback. Mon. Not. R. Astron. Soc. 2016, 463, 3948–3964. [Google Scholar] [CrossRef]
  43. Genel, S.; Vogelsberger, M.; Springel, V.; Sijacki, D.; Nelson, D.; Snyder, G.; Rodriguez-Gomez, V.; Torrey, P.; Hernquist, L. Introducing the Illustris project: The evolution of galaxy populations across cosmic time. Mon. Not. R. Astron. Soc. 2014, 445, 175–200. [Google Scholar] [CrossRef]
  44. Schaye, J.; Crain, R.A.; Bower, R.G.; Furlong, M.; Schaller, M.; Theuns, T.; Dalla Vecchia, C.; Frenk, C.S.; McCarthy, I.G.; Helly, J.C.; et al. The EAGLE project: Simulating the evolution and assembly of galaxies and their environments. Mon. Not. R. Astron. Soc. 2015, 446, 521–554. [Google Scholar] [CrossRef]
  45. Pillepich, A.; Springel, V.; Nelson, D.; Genel, S.; Naiman, J.; Pakmor, R.; Hernquist, L.; Torrey, P.; Vogelsberger, M.; Weinberger, R.; et al. Simulating galaxy formation with the IllustrisTNG model. Mon. Not. R. Astron. Soc. 2018, 473, 4077–4106. [Google Scholar] [CrossRef]
  46. Davé, R.; Anglés-Alcázar, D.; Narayanan, D.; Li, Q.; Rafieferantsoa, M.H.; Appleby, S. SIMBA: Cosmological simulations with black hole growth and feedback. Mon. Not. R. Astron. Soc. 2019, 486, 2827–2849. [Google Scholar] [CrossRef]
  47. Ni, Y.; Di Matteo, T.; Bird, S.; Croft, R.; Feng, Y.; Chen, N.; Tremmel, M.; DeGraf, C.; Li, Y. The ASTRID simulation: The evolution of supermassive black holes. Mon. Not. R. Astron. Soc. 2022, 513, 670–692. [Google Scholar] [CrossRef]
  48. Duras, F.; Bongiorno, A.; Ricci, F.; Piconcelli, E.; Shankar, F.; Lusso, E.; Bianchi, S.; Fiore, F.; Maiolino, R.; Marconi, A.; et al. Universal bolometric corrections for active galactic nuclei over seven luminosity decades. Astron. Astrophys. 2020, 636, A73. [Google Scholar] [CrossRef]
  49. Juodžbalis, I.; Conselice, C.J.; Singh, M.; Adams, N.; Ormerod, K.; Harvey, T.; Austin, D.; Volonteri, M.; Cohen, S.H.; Jansen, R.A.; et al. EPOCHS VII: Discovery of high redshift (6.5 < z < 12) AGN candidates in JWST ERO and PEARLS data. arXiv 2023, arXiv:2307.07535. [Google Scholar]
  50. Euclid Collaboration; Barnett, R.; Warren, S.J.; Mortlock, D.J.; Cuby, J.G.; Conselice, C.; Hewett, P.C.; Willott, C.J.; Auricchio, N.; Balaguera-Antolínez, A.; et al. Euclid preparation. V. Predicted yield of redshift 7 < z < 9 quasars from the wide survey. Astron. Astrophys. 2019, 631, A85. [Google Scholar] [CrossRef]
  51. Pacucci, F.; Loeb, A. The search for the farthest quasar: Consequences for black hole growth and seed models. Mon. Not. R. Astron. Soc. 2022, 509, 1885–1891. [Google Scholar] [CrossRef]
  52. Natarajan, P.; Pacucci, F.; Ricarte, A.; Bogdán, Á.; Goulding, A.D.; Cappelluti, N. First Detection of an Overmassive Black Hole Galaxy UHZ1: Evidence for Heavy Black Hole Seed Formation from Direct Collapse. Astrophys. J. 2024, 960, L1. [Google Scholar] [CrossRef]
  53. Pacucci, F.; Ferrara, A.; Volonteri, M.; Dubus, G. Shining in the dark: The spectral evolution of the first black holes. Mon. Not. R. Astron. Soc. 2015, 454, 3771–3777. [Google Scholar] [CrossRef]
  54. Valiante, R.; Schneider, R.; Zappacosta, L.; Graziani, L.; Pezzulli, E.; Volonteri, M. Chasing the observational signatures of seed black holes at z > 7: Candidate observability. Mon. Not. R. Astron. Soc. 2018, 476, 407–420. [Google Scholar] [CrossRef]
  55. Pacucci, F.; Baldassare, V.; Cappelluti, N.; Fan, X.; Ferrara, A.; Haiman, Z.; Natarajan, P.; Ozel, F.; Schneider, R.; Tremblay, G.R.; et al. Detecting the Birth of Supermassive Black Holes Formed from Heavy Seeds. Bull. Am. Astron. Soc. 2019, 51, 117. [Google Scholar] [CrossRef]
  56. Pacucci, F.; Natarajan, P.; Volonteri, M.; Cappelluti, N.; Urry, C.M. Conditions for Optimal Growth of Black Hole Seeds. Astrophys. J. 2017, 850, L42. [Google Scholar] [CrossRef]
  57. Visbal, E.; Haiman, Z. Identifying Direct Collapse Black Hole Seeds through Their Small Host Galaxies. Astrophys. J. 2018, 865, L9. [Google Scholar] [CrossRef]
  58. Naidu, R.P.; Oesch, P.A.; van Dokkum, P.; Nelson, E.J.; Suess, K.A.; Brammer, G.; Whitaker, K.E.; Illingworth, G.; Bouwens, R.; Tacchella, S.; et al. Two Remarkably Luminous Galaxy Candidates at z ≈ 10-12 Revealed by JWST. Astrophys. J. 2022, 940, L14. [Google Scholar] [CrossRef]
  59. Adams, N.J.; Conselice, C.J.; Ferreira, L.; Austin, D.; Trussler, J.A.A.; Juodžbalis, I.; Wilkins, S.M.; Caruana, J.; Dayal, P.; Verma, A.; et al. Discovery and properties of ultra-high redshift galaxies (9 < z < 12) in the JWST ERO SMACS 0723 Field. Mon. Not. R. Astron. Soc. 2023, 518, 4755–4766. [Google Scholar] [CrossRef]
  60. Marchesi, S.; Gilli, R.; Lanzuisi, G.; Dauser, T.; Ettori, S.; Vito, F.; Cappelluti, N.; Comastri, A.; Mushotzky, R.; Ptak, A.; et al. Mock catalogs for the extragalactic X-ray sky: Simulating AGN surveys with ATHENA and with the AXIS probe. Astron. Astrophys. 2020, 642, A184. [Google Scholar] [CrossRef]
  61. Dauser, T.; Falkner, S.; Lorenz, M.; Kirsch, C.; Peille, P.; Cucchetti, E.; Schmid, C.; Brand, T.; Oertel, M.; Smith, R.; et al. SIXTE: A generic X-ray instrument simulation toolkit. Astron. Astrophys. 2019, 630, A66. [Google Scholar] [CrossRef]
  62. Vito, F.; Brandt, W.N.; Yang, G.; Gilli, R.; Luo, B.; Vignali, C.; Xue, Y.Q.; Comastri, A.; Koekemoer, A.M.; Lehmer, B.D.; et al. High-redshift AGN in the Chandra Deep Fields: The obscured fraction and space density of the sub-L🟉 population. Mon. Not. R. Astron. Soc. 2018, 473, 2378–2406. [Google Scholar] [CrossRef]
  63. Ivezić, Ž.; Kahn, S.M.; Tyson, J.A.; Abel, B.; Acosta, E.; Allsman, R.; Alonso, D.; AlSayyad, Y.; Anderson, S.F.; Andrew, J.; et al. LSST: From Science Drivers to Reference Design and Anticipated Data Products. Astrophys. J. 2019, 873, 111. [Google Scholar] [CrossRef]
  64. Bourke, T.; Braun, R.; Bonaldi, A.; Garcia-Miro, C.; Keane, E.; Wagg, J.; SKAO Science Team. Expected Science Performance of the Square Kilometre Array Phase 1 (SKA1). Am. Astron. Soc. Meet. Abstr. 2018, 231, 152.07. [Google Scholar]
  65. Civano, F.; Marchesi, S.; Comastri, A.; Urry, M.C.; Elvis, M.; Cappelluti, N.; Puccetti, S.; Brusa, M.; Zamorani, G.; Hasinger, G.; et al. The Chandra Cosmos Legacy Survey: Overview and Point Source Catalog. Astrophys. J. 2016, 819, 62. [Google Scholar] [CrossRef]
  66. Ochsenbein, F.; Bauer, P.; Marcout, J. The VizieR database of astronomical catalogues. Astron. Astrophys. 2000, 143, 23–32. [Google Scholar] [CrossRef]
  67. Iye, M.; Karoji, H.; Ando, H.; Kaifu, N.; Kodaira, K.; Aoki, K.; Aoki, W.; Chikada, Y.; Doi, Y.; Ebizuka, N.; et al. Current Performance and On-Going Improvements of the 8.2 m Subaru Telescope. Publ. Astron. Soc. Jpn. 2004, 56, 381–397. [Google Scholar] [CrossRef]
  68. Salvato, M.; Buchner, J.; Budavári, T.; Dwelly, T.; Merloni, A.; Brusa, M.; Rau, A.; Fotopoulou, S.; Nandra, K. Finding counterparts for all-sky X-ray surveys with NWAY: A Bayesian algorithm for cross-matching multiple catalogues. Mon. Not. R. Astron. Soc. 2018, 473, 4937–4955. [Google Scholar] [CrossRef]
  69. Buchner, J.; Salvato, M.; BudavaÌri, T.; Fotopoulou, S. nway: Bayesian cross-matching of astronomical catalogs. arXiv 2021, arXiv:2102.014. [Google Scholar]
  70. Finkelstein, S.L.; Bagley, M.B.; Ferguson, H.C.; Wilkins, S.M.; Kartaltepe, J.S.; Papovich, C.; Yung, L.Y.A.; Arrabal Haro, P.; Behroozi, P.; Dickinson, M.; et al. CEERS Key Paper. I. An Early Look into the First 500 Myr of Galaxy Formation with JWST. Astrophys. J. 2023, 946, L13. [Google Scholar] [CrossRef]
  71. Salvato, M.; Ilbert, O.; Hasinger, G.; Rau, A.; Civano, F.; Zamorani, G.; Brusa, M.; Elvis, M.; Vignali, C.; Aussel, H.; et al. Dissecting Photometric Redshift for Active Galactic Nucleus Using XMM- and Chandra-COSMOS Samples. Astrophys. J. 2011, 742, 61. [Google Scholar] [CrossRef]
Figure 1. Evolution tracks for black hole seeds (light and heavy), assuming Eddington-limited growth. The heavy seeds range between 10 4 and 10 6 M in initial mass, while the light seeds range between 10 M and 100 M . Seeding occurs at a typical redshift of z∼25 [13]. The black points show the masses and redshifts of known z > 6 quasars [5,14,15]. The X-ray bolometric correction is assumed to be 10%. AXIS will be able to detect all black holes within the green-shaded area.
Figure 1. Evolution tracks for black hole seeds (light and heavy), assuming Eddington-limited growth. The heavy seeds range between 10 4 and 10 6 M in initial mass, while the light seeds range between 10 M and 100 M . Seeding occurs at a typical redshift of z∼25 [13]. The black points show the masses and redshifts of known z > 6 quasars [5,14,15]. The X-ray bolometric correction is assumed to be 10%. AXIS will be able to detect all black holes within the green-shaded area.
Universe 10 00276 g001
Figure 2. The large image on the left displays a simulated AXIS deep 7 Ms field. The two top right panels provide a zoomed-in view of the green area, showcasing a comparison of the on-axis performances between AXIS and Chandra. Within this area, an extended source was included in the simulation to demonstrate how AXIS can readily detect point sources embedded within diffuse emission. The central and bottom right panels depict the yellow and red areas, respectively. The performance of AXIS optics exhibits minimal degradation in angular resolution even at large off-axis angles, whereas Chandra’s PSF and effective area experience substantial degradation. AXIS’s performance remains largely unaffected by the off-axis angle.
Figure 2. The large image on the left displays a simulated AXIS deep 7 Ms field. The two top right panels provide a zoomed-in view of the green area, showcasing a comparison of the on-axis performances between AXIS and Chandra. Within this area, an extended source was included in the simulation to demonstrate how AXIS can readily detect point sources embedded within diffuse emission. The central and bottom right panels depict the yellow and red areas, respectively. The performance of AXIS optics exhibits minimal degradation in angular resolution even at large off-axis angles, whereas Chandra’s PSF and effective area experience substantial degradation. AXIS’s performance remains largely unaffected by the off-axis angle.
Universe 10 00276 g002
Figure 4. The expected source counts from both semi-analytical models (SAMs) and cosmological simulations are shown in the figure. The red and blue curves represent heavy- and light-seeding scenarios, respectively, from the SAM developed by Ricarte and Natarajan [33] (see text for more details). These curves diverge at redshifts above 9, highlighting the different predictions for these scenarios. The cosmological simulations compiled by Habouzit et al. [41] exhibit a wide range of luminosity functions, underscoring the need for additional high redshift constraints to better constrain the models. Additionally, at redshift 6, the 3 σ upper uncertainty on the number counts from the X-ray luminosity function by Vito et al. [22] is shown for reference. This uncertainty is extrapolated from observational evidence at redshift 4, providing additional context for the expected source counts at high redshifts.
Figure 4. The expected source counts from both semi-analytical models (SAMs) and cosmological simulations are shown in the figure. The red and blue curves represent heavy- and light-seeding scenarios, respectively, from the SAM developed by Ricarte and Natarajan [33] (see text for more details). These curves diverge at redshifts above 9, highlighting the different predictions for these scenarios. The cosmological simulations compiled by Habouzit et al. [41] exhibit a wide range of luminosity functions, underscoring the need for additional high redshift constraints to better constrain the models. Additionally, at redshift 6, the 3 σ upper uncertainty on the number counts from the X-ray luminosity function by Vito et al. [22] is shown for reference. This uncertainty is extrapolated from observational evidence at redshift 4, providing additional context for the expected source counts at high redshifts.
Universe 10 00276 g004
Figure 5. The figure displays the 0.5–2 keV area–flux curves for the AXIS deep, intermediate, and wide Serendipitous surveys2, as outlined in this paper. The plotted lines represent the sensitivity curves of the surveys in the 0.5–2 keV range, covering an area range of 20% to 80% of the survey coverage. The light blue and pink areas indicate the parameter space covered by the SAM heavy and light black hole (BH) seed models, respectively. These areas represent the predicted range of AGN detections in the redshift ranges of z = [8–9] and z = [9–10], assuming the detection of 10 AGN in each range. As discussed in the text, the AXIS deep and intermediate surveys have the potential to detect a sufficient number of sources at z 8–10, resulting in a significant reduction in the prediction range of the current SAM models. For comparison, the figure also includes a yellow area representing the parameter space covered by current X-ray surveys using Chandra, XMM-Newton, and eROSITA. This provides context for the coverage and capabilities of these existing surveys compared to the proposed AXIS surveys.
Figure 5. The figure displays the 0.5–2 keV area–flux curves for the AXIS deep, intermediate, and wide Serendipitous surveys2, as outlined in this paper. The plotted lines represent the sensitivity curves of the surveys in the 0.5–2 keV range, covering an area range of 20% to 80% of the survey coverage. The light blue and pink areas indicate the parameter space covered by the SAM heavy and light black hole (BH) seed models, respectively. These areas represent the predicted range of AGN detections in the redshift ranges of z = [8–9] and z = [9–10], assuming the detection of 10 AGN in each range. As discussed in the text, the AXIS deep and intermediate surveys have the potential to detect a sufficient number of sources at z 8–10, resulting in a significant reduction in the prediction range of the current SAM models. For comparison, the figure also includes a yellow area representing the parameter space covered by current X-ray surveys using Chandra, XMM-Newton, and eROSITA. This provides context for the coverage and capabilities of these existing surveys compared to the proposed AXIS surveys.
Universe 10 00276 g005
Table 1. Properties of three reference AXIS surveys simulated in this work. The flux limits are calculated at 1 %, 20 %, and 80 % of the covered field and correspond to the fluxes at which 1 %, 20 %, and 80 % of the input sources are detected. The number of detections is computed in the 0.5–7 keV band.
Table 1. Properties of three reference AXIS surveys simulated in this work. The flux limits are calculated at 1 %, 20 %, and 80 % of the covered field and correspond to the fluxes at which 1 %, 20 %, and 80 % of the input sources are detected. The number of detections is computed in the 0.5–7 keV band.
SurveyArea
deg2
Tile Exposure
ks
Total Exposure
Ms
Flux Limit (0.5–2 keV)
erg s−1 cm−2
Number of
AGN Detections
1 %20 %80 %
Deep0.13580074.5 × 10−191.9 × 10−184.3 × 10−182800
Intermediate236064.5 × 10−181.1 × 10−172.6 × 10−1721,000
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cappelluti, N.; Foord, A.; Marchesi, S.; Pacucci, F.; Ricarte, A.; Habouzit, M.; Vito, F.; Powell, M.; Koss, M.; Mushotzky, R. Surveying the Onset and Evolution of Supermassive Black Holes at High-z with AXIS. Universe 2024, 10, 276. https://doi.org/10.3390/universe10070276

AMA Style

Cappelluti N, Foord A, Marchesi S, Pacucci F, Ricarte A, Habouzit M, Vito F, Powell M, Koss M, Mushotzky R. Surveying the Onset and Evolution of Supermassive Black Holes at High-z with AXIS. Universe. 2024; 10(7):276. https://doi.org/10.3390/universe10070276

Chicago/Turabian Style

Cappelluti, Nico, Adi Foord, Stefano Marchesi, Fabio Pacucci, Angelo Ricarte, Melanie Habouzit, Fabio Vito, Meredith Powell, Michael Koss, and Richard Mushotzky. 2024. "Surveying the Onset and Evolution of Supermassive Black Holes at High-z with AXIS" Universe 10, no. 7: 276. https://doi.org/10.3390/universe10070276

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