Signatures of Exploding Supermassive PopIII Stars at High Redshift in JWST, EUCLID and Roman Space Telescope

Cosmic voids, the underdense regions of the cosmic web, are widely used to constrain cosmology. Voids contain few, isolated galaxies, presumably expected to be less evolved and preserving memory of the pristine Universe. We use the cosmological hydrodynamical simulation Horizon-AGN coupled to the void finder vide to investigate properties of galaxies in voids at z = 0. We find that, closer to void centres, low-mass galaxies are more common than their massive counterparts. At a fixed dark matter halo mass, they have smaller stellar masses than in denser regions. The star formation rate of void galaxies diminishes when approaching void centres, but their specific star formation rate slightly increases, suggesting that void galaxies form stars more efficiently with respect to their stellar mass. We find that this cannot only be attributed to the prevalence of low-mass galaxies. The inner region of voids also predominantly hosts low-mass black holes (BHs). However, the BH mass-to-galaxy mass ratios resemble those of the whole simulation at z = 0. Our results suggest that even if the growth channels in cosmic voids are different from those in denser environments, voids grow their galaxies and BHs in a similar way. While a large fraction of the BHs have low Eddington ratios, we find that $\text{$\sim$} 20{{\ \rm per\ cent}}$ could be observed as active galactic nuclei with $\log _{10} L_{\rm 2\!-\!10 \, keV}=41.5\!-\!42.5 \, \rm erg\, s^{-1}$. These results pave the way to future work with larger next-generation hydro-simulations, aiming to confirm our findings and prepare the application on data from upcoming large surveys such as Prime Focus Spectrograph, Euclid, and Wide Field Infrared Survey Telescope.

Context. In the supermassive black hole (SMBH)-galaxy coevolution scenario, heavily obscured active galactic nuclei (AGN) represent a fundamental phase of SMBH growth during which most of the BH mass is accreted and the scaling relations with the host galaxy are set. Obscured nuclei are thought to constitute a major fraction of the whole AGN population, but their statistics and evolution across cosmic time are still highly uncertain. Therefore, it is pivotal to identify new ways to detect this vast and hidden population of growing SMBHs. A promising way to select heavily obscured AGN is through radio emission, which is largely unaffected by obscuration and can be used as a proxy for nuclear activity. Aims. In this work, we study the AGN radio detection effectiveness in the major deep extragalactic surveys, considering different AGN obscuration levels, redshift, and AGN bolometric luminosities. We particularly focus on comparing their radio and X-ray detectability, making predictions for present and future radio surveys. Methods. We extrapolated the predictions of the AGN population synthesis model of the cosmic X-ray background (CXB) to the radio band, by deriving the 1.4 GHz luminosity functions of unobscured (i.e., with hydrogen column densities log NH < 22), obscured (22 < log NH < 24), and Compton-thick (CTK, log NH > 24) AGN. We then used these functions to forecast the number of detectable AGN based on the area, flux limit, and completeness of a given radio survey and compare it with the AGN number resulting from X-ray predictions. Results. When applied to deep extragalactic fields covered both by radio and X-ray observations, we show that, while X-ray selection is generally more effective in detecting unobscured AGN, the surface density of CTK AGN radio detected is on average ten times larger than the X-ray one, and even greater at high redshifts, considering the current surveys and facilities. Our results suggest that thousands of CTK AGN are already present in current radio catalogs, but most of them escaped any detection in the corresponding X-ray observations. We also present expectations for the number of AGN to be detected by the Square Kilometer Array Observatory (SKAO) in its future deep and wide radio continuum surveys, finding that it will be able to detect more than 2000 AGN at z > 6 and tens of them at z > 10, more than half of which are expected to be CTK.

Variability is a fundamental signature for active galactic nuclei (AGN) activity and serves as an unbiased indicator for rapid instability happening near the center of supermassive black holes (SMBHs). Previous studies showed that AGN variability does not have strong redshift evolution, and scales with their bolometric luminosity and BH mass, making it a powerful probe to identify low-mass, low-luminosity AGNs at high redshift. JWST has discovered a new population of high-redshift galaxies likely hosting moderate accreting BHs (>106 M ⊙)—the little red dots (LRDs; z ∼ 4–10). In this Letter, we study the variability of a sample of 22 LRDs with V-shaped spectral energy distributions in three JWST deep fields that also have reliable Hubble Space Telescope observations in closely paired filters at 1–2 μm (rest-frame UV), with the time difference between 6 and 11 yr. This LRD sample covers a redshift range of 3 < z < 8 with −21.3 < M UV < −18.4. Based on both photometry and imaging difference analyses, we find a mean magnitude difference of ∼0.15 ± 0.26 mag, with none of the LRDs showing photometric variability at 3σ significance. Extrapolation of Sloan Digital Sky Survey quasar variability predicts a magnitude change of order 0.3 mag for our LRD sample. This suggests an upper limit of about ∼30% AGN contribution to the total observed UV light in our sample of LRDs.

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.

JWST observations have opened a new chapter in supermassive black hole (SMBH) studies, stimulating discussion of two puzzles: the abundance of high-z SMBHs and the fraction of dual active galactic nuclei (AGNs). We argue that the answers to these puzzles may be linked to an interpretation of the data on the nanohertz gravitational waves (GWs) discovered by NANOGrav and other pulsar timing arrays as SMBH binaries whose evolution is driven by interactions with their environments down to O(0.1 pc) separations. We show that the stellar mass-black hole mass correlations found in JWST data and in low-ɀ inactive galaxies are similar, and present a global fit to these data, excluding low-ɀ AGNs. Matching the NANOGrav and dual-AGN data requires that binary evolution due to environmental effects at separations below O(1 kpc) be rapid on cosmological timescales. According to this interpretation, the SMBHs in low-ɀ AGNs are the tip of the iceberg of a local SMBH population in mainly inactive galaxies. This interpretation is consistent with the ‘little red dots’ observed with JWST being AGNs, and would favour the observability of GW signals from black hole binaries in LISA and decihertz GW detectors.

Context. The assembly and co-evolution of super-massive black holes (SMBHs) and their host galaxy stellar population is one of the key open questions in modern galaxy evolution. Observationally constraining this question is challenging. Important parameters of galaxies, such as the stellar mass (M⋆) and star formation rate (SFR), are inferred by modeling the spectral energy distribution (SED), with templates constructed on the basis of various assumptions on stellar evolution. In the case of galaxies triggering SMBH activity, the active galactic nucleus (AGN) contaminates the light of the host galaxy at all wavelengths, hampering inferences of host galaxy parameters. Underestimating the AGN contribution due to incomplete AGN templates results in a systematic overestimation of the stellar mass, biasing our understanding of AGN and galaxy co-evolution. This challenge has gained further attention with the advent of sensitive wide-area surveys with millions of newly detected luminous AGN, including those by eROSITA, Euclid, and LSST. Aims. We aim to robustly estimate the accuracy, bias, scatter, and uncertainty of AGN host galaxy parameters, including stellar masses, and improve these measurements relative to previously used techniques. Methods. This work makes two important contributions. Firstly, we present a new SED fitting code, GRAHSP, with an AGN model composed of a flexible power-law continuum with empirically determined broad and narrow lines and a FeII forest component, a flexible infrared torus that can reproduce the diverse dust temperature distributions, and appropriate attenuation on the galaxy and AGN light components. We verify that this model reproduces published X-ray to infrared SEDs of AGN to better than 20% accuracy. A fully Bayesian fit includes uncertainties in the model and the data, making the inference highly robust. The model is constrained with a fast nested sampling inference procedure supporting the many free model parameters. Secondly, we created a benchmark photometric data set where optically selected pure quasars are paired with non-AGN pure galaxies at the same redshift. Their photometry flux is summed into a hybrid (Chimera) object but with known galaxy and AGN properties. Based on this data-driven benchmark, true and retrieved stellar masses, SFR, and AGN luminosities can be compared, allowing for the evaluation and quantification of biases and uncertainties inherent in any given SED fitting methodology. Results. The Chimera benchmark, which we release with this paper, shows that previous codes systematically overestimate M⋆ and SFR by 0.5 dex with a wide scatter of 0.7 dex at AGN luminosities above 1044 erg s−1. In 20% of cases, the estimated error bars lie completely outside a 1 dex-wide band centreed around the true value, which we consider an outlier. In contrast, GRAHSP shows no measurable bias on M⋆ and SFR, with an outlier fraction of only about 5%. GRAHSP also estimates more realistic uncertainties. Conclusions. Unbiased characterization of galaxies hosting AGN enables characterization of the environmental conditions conducive to black hole growth, whether star formation is suppressed at high black hole activity, and identifying the mechanisms that prevent overluminous AGN relative to the host galaxy mass. It can also shed light on the long-standing questions of whether AGN obscuration is primarily an orientation effect or related to phases in galaxy evolution.

Accretion is the dominant contribution to the cosmic massive black hole (MBH) density in the Universe today. However, modelling accretion in cosmological simulations is challenging due to the dynamic range involved, as well as the theoretical uncertainties of the underlying mechanisms driving accretion from galactic to black hole horizon scales. We present a simple, flexible parametrisation for gas inflows onto MBHs aimed at managing this uncertainty in the context of large-volume cosmological simulations. This study has been carried out as part of the ‘Learning the Universe’ collaboration, with an aim to jointly infer the initial conditions and physical processes governing the evolution of the Universe, using a Bayesian forward-modelling approach. To allow for this forward-modelling approach, we updated the prescription for accretion with a two-parameter free-fall based inflow estimate that allows for a radius-dependent inflow rate and added a simple model for unresolved accretion disks. We used uniform resolution cosmological hydrodynamical simulations and the IllustrisTNG framework to study the MBH population and its dependence on the introduced model parameters. Once the parameters of the accretion formula were chosen, aimed at achieving a zero black hole mass density (BHMD) at a roughly similar redshift, the differences caused by details in the accretion formula are moderate in the supermassive black hole (SMBH) regime, indicating that it is difficult to distinguish between accretion mechanisms based on luminous active galactic nuclei (AGNs) powered by SMBHs. Applying the same models to intermediate-mass black holes (IMBHs) at high redshifts, however, reveals significantly different accretion rates in high-redshift, moderate-luminosity AGNs, as well as different frequencies and mass distributions of IMBH mergers for the same black hole formation model. This difference in the early growth history will also likely lead to an accretion model-dependent SMBH population in low-mass black hole formation scenarios.

Destruction of a Star during the Evolution of a Star + Supermassive Black Hole System

Context. Large-scale environment is one of the main physical drivers of galaxy evolution. The densest regions at high redshifts (i.e. z > 2 protoclusters) are gas-rich regions characterised by high star formation activity. The same physical properties that enhance star formation in protoclusters are also thought to boost the growth of supermassive black holes (SMBHs), most likely in heavily obscured conditions. Aims. We aim to test this scenario by probing the active galactic nucleus (AGN) content of SPT2349–56: a massive, gas-rich, and highly star-forming protocluster core at z = 4.3 discovered as an overdensity of dusty star-forming galaxies (DSFGs). We compare our results with data on the field environment and other protoclusters. Methods. We observed SPT2349–56 with Chandra (200 ks) and searched for X-ray emission from the known galaxy members. We also performed a spectral energy distribution fitting procedure to derive the physical properties of the discovered AGNs. Results. In the X-ray band, we detected two protocluster members: C1 and C6, corresponding to an AGN fraction among DSFGs in the structure of ≈10%. This value is consistent with other protoclusters at z = 2 − 4, but higher than the AGN incidence among DSFGs in the field environment. Both AGNs are heavily obscured sources, hosted in star-forming galaxies with ≈3 × 1010 M⊙ stellar masses. We estimate that the intergalactic medium in the host galaxies contributes to a significant fraction (or even entirely) to the nuclear obscuration. In particular, C1 is a highly luminous (LX = 2 × 1045 erg s−1) and Compton-thick (NH = 2 × 1024 cm−2) AGN, likely powered by a MBH > 6 × 108 M⊙ SMBH, assuming Eddington-limited accretion. Its high accretion rate suggests that it is in the phase of efficient growth that is generally required to explain the presence of extremely massive SMBHs in the centres of local galaxy clusters. Considering SPT2349–56 and DRC, a similar protocuster at z = 4, and under different assumptions on their volumes, we find that gas-rich protocluster cores at z ≈ 4 enhance the triggering of luminous (logLX/erg s−1 = 45 − 46) AGNs by three to five orders of magnitude with respect to the predictions from the AGN X-ray luminosity function at a similar redshift in the field environment. We note that this result is not solely driven by the overdensity of the galaxy population in the structures. Conclusions. Our results indicate that gas-rich protoclusters at high redshift boost the growth of SMBHs, which will likely impact the subsequent evolution of the structures. Therefore, they stand as key science targets to obtain a complete understanding of the relation between the environment and galaxy evolution. Dedicated investigations of similar protoclusters are required to definitively confirm this conclusion with a higher statistical significance.

This review examines the connection between X-ray-selected Active Galactic Nuclei (AGN) and their host galaxies, focusing on how X-ray observations provide insights into AGN structure and clustering. AGNs, powered by supermassive black holes, are key drivers of galaxy evolution, and X-ray data play a critical role in studying these energetic phenomena. The unified model of AGNs, which attributes differences between type 1 (unobscured) and type 2 (obscured) AGNs to orientation effects, is discussed. However, variations in clustering between these two types challenge this model, suggesting additional factors influence their evolution. Detecting AGN clusters in the X-ray band remains difficult due to observational biases and limitations, but such studies are vital for understanding how AGNs form and interact within large-scale structures. Host galaxy properties, including luminosity, stellar mass, and star formation rate, are analyzed for their impact on AGN clustering. Research indicates that AGN luminosity is strongly linked to the mass of the dark matter halos surrounding their host galaxies. This relationship may vary depending on the triggering mechanism of the AGN, such as galaxy mergers or internal instabilities. Differences in AGN clustering patterns provide insights into the diverse pathways through which AGNs are activated. AGN feedback, which describes how AGNs influence star formation in their host galaxies, is another key focus. Observations suggest that at higher redshifts, brighter AGNs tend to enhance star formation rates, showing a complex interplay between AGN activity and galaxy growth. By synthesizing recent observational results, this review highlights the central role of AGNs in shaping galaxies and their environments. It provides a deeper understanding of how AGNs interact with their host galaxies and larger cosmic structures, offering valuable insights into the processes driving galaxy evolution over cosmic time.

Understanding how Active Galactic Nuclei (AGN) evolve through cosmic time allows us to probe the physical processes that control their evolution. We use an updated model for the evolution of masses and spins of supermassive black holes (SMBHs), coupled to the latest version of the semi-analytical model of galaxy formation GALFORM using the Planck cosmology and a high resolution Millennium style dark matter simulation to make predictions for AGN and SMBH properties for $0 < z < 6$. We compare the model to the observed black hole mass function and the SMBH versus galaxy bulge mass relation at $z=0$, and compare the predicted bolometric, hard X-ray, soft X-ray and optical AGN luminosity functions to observations at $z < 6$, and find that the model is in good agreement with the observations. The model predicts that at $z<2$ and $L_{\mathrm{bol}} < 10^{43} \mathrm{ergs^{-1}}$, the AGN luminosity function is dominated by objects accreting in an Advection Dominated Accretion Flow (ADAF) disc state, while at higher redshifts and higher luminosities the dominant contribution is from objects accreting via a thin disc or at super-Eddington rates. The model also predicts that the AGN luminosity function at $z<3$ and $L_{\mathrm{bol}} < 10^{44} \mathrm{ergs^{-1}}$ is dominated by the contribution from AGN fuelled by quiescent hot halo accretion, while at higher luminosities and higher redshifts, the AGN luminosity function is dominated by the contribution from AGN fuelled by starbursts triggered by disc instabilities. We employ this model to predict the evolution of SMBH masses, Eddington ratios, and spins, finding that the median SMBH spin evolves very little for $0<z<6$.

This PhD Thesis is devoted to the accurate analysis of the physical properties of Active Galactic Nuclei (AGN) and the AGN/host-galaxy interplay. Due to the broad-band AGN emission (from radio to hard X-rays), a multi-wavelength approach is mandatory. Our research is carried out over the COSMOS field, within the context of the XMM-Newton wide-field survey. To date, the COSMOS field is a unique area for comprehensive multi-wavelength studies, allowing us to define a large and homogeneous sample of QSOs with a well-sampled spectral coverage and to keep selection effects under control. Moreover, the broad-band information contained in the COSMOS database is well-suited for a detailed analysis of AGN SEDs, bolometric luminosities and bolometric corrections. In order to investigate the nature of both obscured (Type-2) and unobscured (Type-1) AGN, the observational approach is complemented with a theoretical modelling of the AGN/galaxy co-evolution. The X-ray to optical properties of an X-ray selected Type-1 AGN sample are discussed in the first part. The relationship between X-ray and optical/UV luminosities, parametrized by the spectral index αox, provides a first indication about the nature of the central engine powering the AGN. Since a Type-1 AGN outshines the surrounding environment, it is extremely difficult to constrain the properties of its host-galaxy. Conversely, in Type-2 AGN the host-galaxy light is the dominant component of the optical/near-IR SEDs, severely affecting the recovery of the intrinsic AGN emission. Hence a multi-component SED-fitting code is developed to disentangle the emission of the stellar populationof the galaxy from that associated with mass accretion. Bolometric corrections, luminosities, stellar masses and star-formation rates, correlated with the morphology of Type-2 AGN hosts, are presented in the second part, while the final part concerns a physically-motivated model for the evolution of spheroidal galaxies with a central SMBH. The model is able to reproduce two important stages of galaxy evolution, namely the obscured cold-phase and the subsequent quiescent hot-phase.

Bolometric luminosity is an important quantity that tells us the radiative energy output of an active galactic nucleus (AGN). A common way to estimate bolometric luminosity is to use ultraviolet (UV) or optical luminosities as its proxies, but the UV- or optical-based-bolometric luminosity estimators can be easily affected by dust extinction. In this study, we present new methods for estimating bolometric luminosities using infrared (IR) hydrogen Paschen and Brackett line luminosities to alleviate the dust extinction effects. We show that there exist tight correlations between the bolometric luminosities and the IR hydrogen line luminosities, and present the IR hydrogen line-based-bolometric luminosity estimators. As an example, we apply the relation to dust obscured red AGNs, confirming previous results that red AGNs have higher Eddington rations than unobscured AGNs. The new bolometric luminosity estimator will be useful for studying obscured AGNs (e.g. red AGNs, Compton thick AGNs, and Type 2 AGNs), during the coming era of sensitive near-infrared (NIR) and mid-infrared (MIR) missions, such as the James Webb Space Telescope.

Correlations between the mass of a supermassive black hole (SMBH) and the properties of its host galaxy (e.g., total stellar mass M *, luminosity L host) suggest an evolutionary connection. A powerful test of a coevolution scenario is to measure the relations  BH –L host and  BH –M * at high redshift and compare with local estimates. For this purpose, we acquired Hubble Space Telescope (HST) imaging with WFC3 of 32 X-ray-selected broad-line (type 1) active galactic nuclei at 1.2 &lt; z &lt; 1.7 in deep survey fields. By applying state-of-the-art tools to decompose the HST images including available ACS data, we measured the host galaxy luminosity and stellar mass along with other properties through the two-dimensional model fitting. The black hole mass (  BH ) was determined using the broad Hα line, detected in the near-infrared with the Subaru Fiber Multi-Object Spectrograph, which potentially minimizes systematic effects using other indicators. We find that the observed ratio of  BH to total M * is 2.7× larger at z ∼ 1.5 than in the local universe, while the scatter is equivalent between the two epochs. A nonevolving mass ratio is consistent with the data at the 2σ–3σ confidence level when accounting for selection effects (estimated using two independent and complementary methods) and their uncertainties. The relationship between  BH and host galaxy total luminosity paints a similar picture. Therefore, our results cannot distinguish whether SMBHs and their total host stellar mass and luminosity proceed in lockstep or whether the growth of the former somewhat overshoots the latter, given the uncertainties. Based on a statistical estimate of the bulge-to-total mass fraction, the ratio  BH /M *,bulge is offset from the local value by a factor of ∼7, which is significant even accounting for selection effects. Taken together, these observations are consistent with a scenario in which stellar mass is subsequently transferred from an angular momentum–supported component of the galaxy to a pressure-supported one through secular processes or minor mergers at a faster rate than mass accretion onto the SMBH.

Context. The determination of the relative frequency of active galactic nuclei (AGN) versus other spectral classes, for example, HII region-like (HII), transition objects (TRAN), passive (PAS), and retired (RET), in a complete set of galaxies in the local Universe is of primary importance to discriminate the source of ionization in the nuclear region of galaxies (e.g., supermassive black holes vs. young and old stars). Aims. Here we aim to provide a spectroscopic characterization of the nuclei of galaxies belonging to the Herschel Reference Survey (HRS), a volume and magnitude limited sample representative of the local Universe, which has become a benchmark for local and high-z studies, for semianalytical models and cosmological simulations. The comparison between the nuclear spectral classification and the one determined on the global galactic scale provides information about how galaxy properties change from the nuclear to the outer regions. Moreover, the extrapolation of the global star formation (SF) properties from the SDSS fiber spectroscopy compared to the one computed by Hα photometry can be useful for testing the method based on aperture correction for determining the global star formation rate for local galaxies. Methods. By collecting the existing nuclear spectroscopy available from the literature, complemented with new observations obtained using the Loiano 1.52 m telescope, we analyze the 322 nuclear spectra of HRS galaxies; their integrated spectroscopy is available from the literature as well. Results. Using two diagnostic diagrams (the BPT and the WHAN) we provide a nuclear and an integrated spectral classification for the HRS galaxies. The BPT and the WHAN methods for nuclei consistently give a frequency of 53–64% HII, around 21–27% AGNs (including TRAN), and 15–20% of PAS (including RET), whereas for integrated spectra they give 69–84% HII, 4–11% of AGNs and 12–20% PAS. Solely among late-type galaxies (LTGs) do the nuclear percentages become 67–77% HII, 22–27% AGNs (including TRAN), and only 1–7% of PAS. For the integrated spectra these frequencies become: 80–85% HII, 9–11% AGNs and 4–9% PAS. Conclusions. We find that the fraction of HII region-like spectra is strongly anticorrelated with the stellar mass. On the contrary the frequency of AGNs increases significantly with stellar mass, such that at M* &gt; 1010.0 M⊙~ 66% of the LTGs are AGNs or TRAN. Moreover there is not a significant dependence of the frequency of AGNs as a function of environment: AGNs+TRAN above 109.0 M⊙ are consistent with ~30% irrespective of their membership to the Virgo cluster, suggesting that the AGNs population is not sensitive to the environment. Finally, extrapolation of the global SF properties from the nuclear spectroscopy including aperture corrections leads to underestimates with respect to values derived from direct integrated Hα photometry.