Black Hole Accretion Rate Research Articles

X-ray AGN in Boötes: the lack of growth of the most massive black holes since z = 4

ABSTRACT Supermassive black holes (BHs) are known to efficiently grow through gas accretion, but even sustained and intense mass build-up through this mechanism struggles to explain the assembly of the most massive BHs observed in the local Universe. Using the Chandra Deep-Wide Field Survey (CDFWS) in the Boötes field, we measure BH–galaxy assembly in massive galaxies ($M_\star \gtrsim 10^{10}\,\rm M_\odot$) through the active galactic nucleus (AGN) fraction and specific black hole accretion rate (sBHAR) distribution as a function of redshift and stellar mass. We determine stellar masses and star formation rates for a parent sample of optically selected galaxies as well as those with X-ray detections indicating the presence of an AGN through spectral energy distribution (SED) fitting. We derive a redshift-dependent mass completeness limit and extract X-ray information for every galaxy as to provide a comprehensive picture of the AGN population in massive galaxies. While X-ray AGN samples are dominated by moderately massive host galaxies of $M_{\star } \geqslant 10^{10}\rm {\rm M}_{\odot }$, we do not find a strong stellar mass dependence in AGN fraction (to limits in sBHAR), indicating a bias towards massive galaxies in the observed samples. We derive BH–galaxy growth tracks over time, which reveal that while most BH mass has been accumulated since $z=4$ for lower mass BHs, the assembly of the most massive BHs is more complex, with little to no relative mass gain since $z=4$, implying that rapid and intense growth episodes prior to $z=4$ were necessary to form these massive BHs.

X-ray reverberation modelling of the continuum, optical/UV time-lags in quasars

Extensive, multi-wavelength monitoring campaigns of nearby and higher redshift active galactic nuclei (AGN) have shown that the UV/optical variations are well correlated with time delays which increase with increasing wavelength. Such behaviour is expected in the context of the X-ray thermal reverberation of the accretion disc in AGN. Our main objective is to use time-lag measurements of luminous AGN and fit them with sophisticated X-ray reverberation time-lags models. In this way we can investigate whether X-ray reverberation can indeed explain the observed continuum time lags, and whether time-lag measurements can be used to measure physical parameters such as the X-ray corona height and the spin of the black hole (BH) in these systems. We use archival time-lag measurements for quasars from different surveys, and we compute their rest frame, mean time-lags spectrum. We fit the data with analytical X-ray reverberation models, using $ statistics, and fitting for both maximal and non spinning BHs, for various colour correction values and X-ray corona heights. We found that X-ray reverberation can explain very well the observed time lags, assuming the measured BH mass, accretion rate and X-ray luminosity of the quasars in the sample. The model agrees well with the data both for non-rotating and maximally rotating BHs, as long as the corona height is larger than $ 40$ gravitational radii. This is in agreement with previous results which showed that X-ray reverberation can also explain the disc radius in micro-lensed quasars, for the same corona heights. The corona height we measure depends on the model assumption of a perfectly flat disc. More realistic disc models may result in lower heights for the X-ray corona.

X-ray view of dissipative warm corona in active galactic nuclei

Context. Detections of the X-ray spectra of active galactic nuclei (AGNs) typically reveal a noticeable excess of soft X-rays beyond the extrapolation of the power-law trend observed between 2 and 10 keV. However, the cause of this surplus remains unclear. In the scenario of soft Comptonization, observations suggest a warm corona temperature, ranging from 0.1 to 1 keV, and an optical depth of approximately 10 to 30. Furthermore, according to radiative constraints derived from spectral analyses employing Comptonization models, it has been suggested that most of the accretion power is released within the warm corona. At the same time, the disk beneath it is largely non-dissipative, mainly emitting the reprocessed radiation from the corona. Aims. We tested the dissipative warm corona model using the radiative transfer code TITAN-NOAR on a sample of 82 XMM-Newton EPIC-pn observations of 21 AGNs. Based on a spectral modeling of the X-ray data, we aim to estimate the total amount of internal heating inside the warm corona on top of the accretion disk. Methods. By modeling the 0.3–10 keV EPIC-pn spectra with the TITAN-NOAR model component, we estimated the internal heating and optical depth of the warm corona and checked their correlations with such global parameters as the hot corona spectral index, black hole mass, and accretion rate. From the model normalization, we computed the radial extent of the warm corona on top of the cold accretion disk. Results. Our model infers the presence of dissipative warm corona, with optical depths distributed across the range of ∼6–30 and a total internal heating in the range of ∼1–29 × 10−23 erg s−1 cm3. We do not detect any variation between these properties and global properties, such as the black hole mass and accretion rate. The extent of the warm corona is spread across a broad range, from 7–408 gravitational radii, and we find that the warm corona is more extended for higher accretion rates. Conclusions. Soft excess emission is ubiquitous across a wide mass range and accretion rate in AGNs. We confirm that the warm corona responsible for producing the soft X-ray excess is highly dissipative with greater optical depths being associated with lower internal heating and vice versa. The cold standard accretion disk regulates the extent of the warm corona.

Birth of Rapidly Spinning, Overmassive Black Holes in the Early Universe

The James Webb Space Telescope (JWST) has unveiled numerous massive black holes (BHs) in faint, broad-line active galactic nuclei (AGNs). The discovery highlights the presence of dust-reddened AGN populations, referred to as “little red dots (LRDs),” more abundant than X-ray-selected AGNs, which are less influenced by obscuration. This finding indicates that the cosmic growth rate of BHs within this population does not decrease but rather increases at higher redshifts beyond z ∼ 6. The BH accretion rate density deduced from their luminosity function is remarkably higher than that from other AGN surveys in X-ray and infrared bands. To align the cumulative mass density accreted to BHs with the observed BH mass density at z ≃ 4–5, as derived from the integration of the BH mass function, the radiative efficiency must be doubled from the canonical 10% value, achieving significance beyond the >3σ confidence level. This suggests the presence of rapid spins with 96% of the maximum limit among these BHs under the thin-disk approximation, maintained by prolonged mass accretion instead of chaotic accretion with randomly oriented inflows. Moreover, we derive an upper bound for the stellar mass of galaxies hosting these LRDs, ensuring consistency with galaxy formation in the standard cosmological model, where the host stellar mass is limited by the available baryonic reservoir. Our analysis gives a lower bound for the BH-to-galaxy mass ratio that exceeds the typical value known in the nearby universe and aligns with that for JWST-detected unobscured AGNs. Accordingly, we propose a hypothesis that the dense, dust-rich environments within LRDs facilitate the emergence of rapidly spinning and overmassive BH populations during the epoch of reionization. This scenario predicts a potential association between relativistic jets and other high-energy phenomena with overmassive BHs in the early universe.

The obsidian model: three regimes of black hole feedback

ABSTRACT In theoretical models of galaxy evolution, black hole feedback is a necessary ingredient in order to explain the observed exponential decline in number density of massive galaxies. Most contemporary black hole feedback models in cosmological simulations rely on a constant radiative efficiency (usually $\eta \sim 0.1$) at all black hole accretion rates. We present the obsidian subgrid model, a synthesis model for the spin-dependent radiative efficiencies of three physical accretion rate regimes, i.e. $\eta = \eta (j, \dot{M}_\mathrm{acc})$, for use in large-volume cosmological simulations. The three regimes include: an advection-dominated accretion flow ($\dot{M}_\mathrm{acc}\lt 0.03\, \dot{M}_\mathrm{Edd}$), a quasar-like mode ($0.03 \lt \dot{M}_\mathrm{acc}/ \dot{M}_\mathrm{Edd}\lt 0.3$), and a slim disc mode ($\dot{M}_\mathrm{acc}\gt 0.3\, \dot{M}_\mathrm{Edd}$). Additionally, we include a large-scale powerful jet at low accretion rates. The black hole feedback model we present is a kinetic model that prescribes mass loadings but could be used in thermal models directly using the radiative efficiency. We implement the obsidian model into the simba galaxy evolution model to determine if it is possible to reproduce galaxy populations successfully, and provide a first calibration for further study. Using a $2\times 1024^3$ particle cosmological simulation in a $(150\, \mathrm{cMpc})^3$ volume, we found that the model is successful in reproducing the galaxy stellar mass function, black hole mass–stellar mass relationship, and stellar mass–halo mass relationship. Moving forward, this model opens new avenues for exploration of the impact of black hole feedback on galactic environments.

Hybrid GRMHD and force-free simulations of black hole accretion

ABSTRACT We present a new approach for stably evolving general relativistic magnetohydrodynamic (GRMHD) simulations in regions where the magnetization $\sigma =b^2/\rho c^2$ becomes large. GRMHD codes typically struggle to evolve plasma above $\sigma \approx 100$ in simulations of black hole accretion. To ensure stability, GRMHD codes will inject mass density artificially to the simulation as necessary to keep the magnetization below a ceiling value $\sigma _{\rm max}$. We propose an alternative approach where the simulation transitions to solving the equations of general relativistic force-free electrodynamics (GRFFE) above a magnetization $\sigma _{\rm trans}$. We augment the GRFFE equations in the highly magnetized region with approximate equations to evolve the decoupled field-parallel velocity and plasma energy density. Our hybrid scheme is explicit and easily added to the framework of standard-volume GRMHD codes. We present a variety of tests of our method, implemented in the GRMHD code koral, and we show results from a 3D hybrid GRMHD + GRFFE simulation of a magnetically arrested disc (MAD) around a spinning black hole. Our hybrid MAD simulation closely matches the average properties of a standard GRMHD MAD simulation with the same initial conditions in low magnetization regions, but it achieves a magnetization $\sigma \approx 10^6$ in the evacuated jet funnel. We present simulated horizon-scale images of both simulations at 230 GHz with the black hole mass and accretion rate matched to M87*. Images from the hybrid simulation are less affected by the choice of magnetization cut-off $\sigma _{\rm cut}$ imposed in radiative transfer than images from the standard GRMHD simulation.

Reconstruction of Cosmic Black Hole Growth and Mass Distribution from Quasar Luminosity Functions at z > 4: Implications for Faint and Low-mass Populations in JWST

The evolution of the quasar luminosity function (QLF) is fundamental to understanding the cosmic evolution of black holes (BHs) through their accretion phases. In the era of the James Webb Space Telescope (JWST), Euclid, and Nancy Grace Roman Space Telescope, their unprecedented detection sensitivity and wide survey area can unveil the low-luminosity quasar and low-mass BH population, and provide new insights into quasar host galaxies. We present a theoretical model describing BH growth from initial seeding at z ≳ 20 to ∼ 4, incorporating the duration of accretion episodes, the distribution of Eddington ratios, and the mass dependency of BH accretion rates. By constraining the model parameters with the observed QLFs at 4 ≤ z ≤ 6 across a wide UV luminosity range, we find that the high-redshift BH population grows rapidly at z ≳ 6, and decelerates the pace in subsequent epochs. Toward lower redshifts (z < 6), mass-dependent accretion inhibits the growth of high-mass BHs with M • > 108 M ⊙, leading to mass saturation at M • ≳ 1010 M ⊙. We predict the BH mass function down to M • ∼ 106 M ⊙ for both unobscured and obscured quasar populations at 4 ≤ z ≤ 11, offering a benchmark for future observational tests. Our model accounts for the presence of both bright and faint quasars at z > 4, including those discovered by JWST. Furthermore, our findings suggest two distinct pathways for the early assembly of the BH–galaxy mass correlation: the population with a BH-to-stellar-mass ratio near the local value of M •/M ⋆ ≃ 5 × 10−3 maintains proximity to the relation via moderate growth, while the population that begins to grow above the local relation becomes as overmassive as M •/M ⋆ ∼ 0.01–0.1 by z ∼ 6 via rapid mass accretion.

The Role of Active Galactic Nucleus Winds in Galaxy Formation: Connecting AGN Outflows at Low Redshifts to the Formation/Evolution of Their Host Galaxies

Using Sloan Digital Sky Survey (SDSS) spectra, we applied an automatic method to search for outflows (OFs) in three large samples of narrow-line active galactic nuclei (AGN) at low redshifts (z < 0.4), separated into three spectral activity classes: radio-loud galaxies (RGs), 15,793; radio-quiet Seyfert 2 AGN (Sy2), 18,585; and LINERs, 25,656. In general, the probability of detecting an OF decreases along the sequence Sy1→Sy2→LINER/RG and independently of the AGN class, the wind velocity, traced by W80, increases with the AGN luminosity. Moreover W80 is systematically higher in RGs or any of the other AGN classes when detected in radio. These results support the idea that there are two main modes of production of OF, the radiative mode dominant in radio-quiet AGN and the jet mode dominant in RGs, although both modes could also happen simultaneously at different levels. From the spectra and SDSS photometry, the characteristics of the AGN host galaxies and their supermassive black holes (SMBHs) were also retrieved using the stellar population synthesis code STARLIGHT. This revealed that, independently of the AGN spectral class, (1) galaxy hosts with OFs have systematically later morphological types and higher star formation rates (SFRs) than their counterparts without OF, (2) the AGN occupy different positions in the specific diagnostic diagram (specific black hole accretion rate (sBHAR) versus specific SFR), which suggests they follow different evolutionary paths congruent with the morphology of their galaxy hosts, and (3) they show no evidence of AGN quenching or triggering of star formation. These results are consistent with a scenario explaining the different AGN classes as consequences of different formation processes of galaxies: early-type galaxies (LINERs and RGs) formed bigger bulges and more massive SMBHs, exhausting their reservoir of gas more rapidly than late-type galaxies (Sy2 and Sy1), and thereby quenching their star formation and starving their SMBHs.

Ensemble power spectral density of SDSS quasars in UV/optical bands

Context. Quasar variability has proven to be a powerful tool to constrain the properties of their inner engine and the accretion process onto supermassive black holes. Correlations between UV/optical variability and physical properties have been long studied with a plethora of different approaches and time-domain surveys, although the detailed picture is not yet clear. Aims. We analysed archival data from the SDSS Stripe-82 region to study how the quasar power spectral density (PSD) depends on the black hole mass, bolometric luminosity, accretion rate, redshift, and rest-frame wavelength. We developed a model-independent analysis framework that could be easily applied to upcoming large surveys such as the Legacy Survey of Space and Time (LSST). Methods. We used light curves of 8042 spectroscopically confirmed quasars, observed in at least six yearly seasons in five filters ugriz. We split the sample into bins of similar physical properties containing at least 50 sources, and we measured the ensemble PSD in each of them. Results. We find that a simple power law is a good fit to the power spectra in the frequency range explored. Variability does not depend on the redshift at a fixed wavelength. Instead, both PSD amplitude and slope depend on the black hole mass, accretion rate, and rest-frame wavelength. We provide scaling relations to model the observed variability as a function of the physical properties, and discuss the possibility of a universal PSD shape for all quasars, where frequencies scale with the black hole mass, while normalization and slope(s) are fixed (at any given wavelength and accretion rate).

Modelling the galaxy radio continuum from star formation and active galactic nuclei in the Shark semi-analytic model

ABSTRACT We present a model of radio continuum emission associated with star formation (SF) and active galactic nuclei (AGNs) implemented in the Shark semi-analytic model of galaxy formation. SF emission includes free-free and synchrotron emission, which depend on the free-electron density and the rate of core-collapse supernovae with a minor contribution from supernova remnants, respectively. AGN emission is modelled based on the jet production rate, which depends on the black hole mass, accretion rate, and spin, and includes synchrotron self-absorption. Shark reproduces radio luminosity functions (RLFs) at $1.4\, \rm GHz$ and $150\, \rm MHz$ for 0 ≤ z ≤ 4, and scaling relations between radio luminosity, star formation rate, and infrared luminosity of galaxies in the local and distant universe in good agreement with observations. The model also reproduces observed number counts of radio sources from 150 MHz to 8.4 GHz to within a factor of 2 on average, though larger discrepancies are seen at the very bright fluxes at higher frequencies. We use this model to understand how the radio continuum emission from radio-quiet AGNs can affect the measured RLFs of galaxies. We find current methods to exclude AGNs from observational samples result in large fractions of radio-quiet AGNs contaminating the ‘star-forming galaxies’ selection and a brighter end to the resulting RLFs.We investigate how this affects the infrared-radio correlation (IRRC) and show that AGN contamination can lead to evolution of the IRRC with redshift. Without this contamination, our model predicts a redshift- and stellar mass-independent IRRC, except at the dwarf-galaxy regime.

The X-ray variability of active galactic nuclei: Power spectrum and variance analysis of the Swift/BAT light curves

Aims. We study the X-ray power spectrum of active galactic nuclei (AGN) in order to investigate whether Seyfert I and II power spectra are similar or not and whether AGN variability depends on the mass and accretion rate of black holes as well as to compare the power spectra of AGN with the power spectra of Galactic X-ray black hole binaries. Method. We used 14–195 keV band light curves from the 157-month Swift/BAT hard X-ray survey, and we computed the mean power spectrum and excess variance of AGN in narrow black hole mass and AGN luminosity bins. We fitted a power-law model to the AGN power spectra, and we investigated whether the power spectrum parameters and the excess variance depend on the black hole mass, luminosity, and accretion rate of AGN. Results. We found the Seyfert I and Seyfert II power spectra to be identical, in agreement with AGN unification models. The mean AGN X-ray power spectrum has the same power-law like shape, with a slope of −1 in all AGN irrespective of their luminosity and black hole mass. We did not detect any flattening to a slope of zero at frequencies as low as 10−9 Hz. We detected an anti-correlation between the power spectral density function (PSD) amplitude and the accretion rate, similar to what has been seen in the past in the 2–10 keV band. This implies that the variability amplitude in AGN decreases with an increasing accretion rate. The universal AGN power spectrum is consistent with the mean 2–9 keV band Cyg X-1 power spectrum in its soft state. We detected a small difference in amplitude, but this is probably due to the difference in energy. Conclusions. The mean low-frequency AGN X-ray power spectrum is consistent with the extension of the mean 0.01–25 Hz Cyg X-1 power spectrum in its soft state to lower frequencies. We cannot prove that the mean AGN PSD is analogous to the mean Cyg X-1 PSD in its soft state, as we do not know the location of the high-frequency break in the hard X-ray AGN PSDs. However, if this is the case, then the accretion disc in AGN probably extends to the radius of the innermost circular stable orbit (as is probably the case with the black hole binaries in their soft state). The X-ray corona will then be located on top, illuminating the disc and producing the X-ray reflection and disc reverberation phenomena commonly observed in these objects. Furthermore, the agreement between the PSD amplitude in AGN and the Cyg X-1 (either in the soft or the hard state) over many decades in frequency indicates that the X-ray variability process is probably the same in all accreting objects, irrespective of the mass of the compact object. We plan to investigate this issue further in the near future.

Active galactic nucleus X-ray luminosity function and absorption function in the Early Universe (3 ≤ z ≤ 6)

The X-ray luminosity function (XLF) of active galactic nuclei (AGN) offers a robust tool to study the evolution and the growth of the supermassive black-hole population over cosmic time. Owing to the limited area probed by X-ray surveys, optical surveys are routinely used to probe the accretion in the high-redshift Universe z ≥ 3. However, optical surveys may be incomplete because they are strongly affected by dust redenning. In this work we derive the XLF and its evolution at high redshifts (z ≥ 3) using a large sample of AGN selected in different fields with various areas and depths covering a wide range of luminosities. Additionally, we put the tightest yet constraints on the absorption function in this redshift regime. In particular, we used more than 600 soft X-ray selected (0.5 − 2 keV) high-z sources in the Chandra deep fields, the Chandra COSMOS Legacy survey, and the XMM-XXL northern field. We derived the X-ray spectral properties for all sources via spectral fitting, using a consistent technique and model. To model the parametric form of the XLF and the absorption function, we used a Bayesian methodology, allowing us to correctly propagate the uncertainties for the observed X-ray properties of our sources and also the absorption effects. The evolution of XLF is in agreement with a pure density evolution model similar to what is witnessed at optical wavelengths, although a luminosity-dependent density evolution model cannot be securely ruled out. A large fraction (∼60%) of our sources are absorbed by column densities of NH ≥ 1023 cm−2, while ∼17% of the sources are Compton-Thick. Our results favour a scenario where both the interstellar medium of the host and the AGN torus contribute to the obscuration. The derived black hole accretion rate density is roughly in agreement with the large-scale cosmological hydrodynamical simulations, if one takes into account the results that the X-ray AGN are hosted by massive galaxies, while it differs from that derived using JWST data. The latter could be due to the differences in the AGN and host-galaxy properties.

Impact of Active Galactic Nuclei Feedback on the Dynamics of Gas: A Review across Diverse Environments

This review examines the relationship between black hole activity and kinematic gas–star misalignment in brightest group galaxies (BGGs) with different merger rates. The formation history of galaxy groups is assessed through “age-dating” as an indicator of distinct major mergers involving the BGGs. BGGs within groups characterized by a higher frequency of major mergers are more likely to host active SMBHs. A consistent correlation is identified between the level of black hole activity, as indicated by the 1.4 GHz and 325 MHz radio emissions, and the degree of kinematic misalignment between the gas and stellar components in BGGs. In dynamically fossil groups, where black hole accretion rate is relatively (∼1 dex) lower due to the lack of recent (≤1 Gyr) major mergers, there is reduced (∼30%) misalignment between the gas and stellar components of BGGs compared to non-fossil groups. Additionally, this study reveals that BGGs in non-fossil groups show higher levels of star formation rate and increased occurrence of mergers, contributing to observed color differences. Exploring the properties and dynamics of the gas disk influenced by mechanical AGN feedback through hydrodynamic simulations suggests that AGN wind-induced effects further lead to the persistent gas misalignment in the disk around the supermassive black hole.

Supermassive Black Hole Winds in X-rays -- SUBWAYS. III. A population study on ultra-fast outflows

The detection of blueshifted absorption lines likely associated with ionized iron K-shell transitions in the X-ray spectra of many active galactic nuclei (AGNs) suggests the presence of a highly ionized gas outflowing with mildly relativistic velocities (0.03c - 0.6c) named ultra-fast outflow (UFO). Within the SUBWAYS project, we characterized these winds starting from a sample of 22 radio-quiet quasars at an intermediate redshift (0.1 leq z leq 0.4) and compared the results with similar studies in the literature on samples of local Seyfert galaxies (i.e., 42 radio-quiet AGNs observed with XMM-Newton at z leq 0.1) and high redshift radio-quiet quasars (i.e., 14 AGNs observed with XMM-Newton and Chandra at z geq 1.4). The scope of our work is a statistical study of UFO parameters and incidence considering the key physical properties of the sources, such as supermassive black hole (SMBH) mass, bolometric luminosity, accretion rates, and spectral energy distribution (SED) with the aim of gaining new insights into the UFO launching mechanisms. We find indications that highly luminous AGNs with a steeper X-ray/UV ratio, $ ox $, are more likely to host UFOs. The presence of UFOs is not significantly related to any other AGN property in our sample. These findings suggest that the UFO phenomenon may be transient. Focusing on AGNs with UFOs, other important findings from this work include: (1) faster UFOs have larger ionization parameters and column densities; (2) X-ray radiation plays a more crucial role in driving highly ionized winds compared to UV; (3) the correlation between outflow velocity and luminosity is significantly flatter than what is expected for radiatively driven winds; (4) more massive black holes experience higher wind mass losses, suppressing the accretion of matter onto the black hole; (5) the UFO launching radius is positively correlated with the Eddington ratio. Furthermore, our analysis suggests the involvement of multiple launching mechanisms, including radiation pressure and magneto-hydrodynamic processes, rather than pointing to a single, universally applicable mechanism.

Mapping the Growth of Supermassive Black Holes as a Function of Galaxy Stellar Mass and Redshift

The growth of supermassive black holes is strongly linked to their galaxies. It has been shown that the population mean black hole accretion rate ( BHAR¯ ) primarily correlates with the galaxy stellar mass (M ⋆) and redshift for the general galaxy population. This work aims to provide the best measurements of BHAR¯ as a function of M ⋆ and redshift over ranges of 109.5 < M ⋆ < 1012 M ⊙ and z < 4. We compile an unprecedentedly large sample with 8000 active galactic nuclei (AGNs) and 1.3 million normal galaxies from nine high-quality survey fields following a wedding cake design. We further develop a semiparametric Bayesian method that can reasonably estimate BHAR¯ and the corresponding uncertainties, even for sparsely populated regions in the parameter space. BHAR¯ is constrained by X-ray surveys sampling the AGN accretion power and UV-to-infrared multiwavelength surveys sampling the galaxy population. Our results can independently predict the X-ray luminosity function (XLF) from the galaxy stellar mass function (SMF), and the prediction is consistent with the observed XLF. We also try adding external constraints from the observed SMF and XLF. We further measure BHAR¯ for star-forming and quiescent galaxies and show that star-forming BHAR¯ is generally larger than or at least comparable to the quiescent BHAR¯ .

X-ray stacking reveals average SMBH accretion properties of star-forming galaxies and their cosmic evolution over 4 ≲ z ≲ 7

ABSTRACT With an X-ray stacking analysis of $\simeq 12\, 000$ Lyman-break galaxies (LBGs) using the Chandra Legacy Survey image, we investigate average supermassive black hole (SMBH) accretion properties of star-forming galaxies (SFGs) at 4 ≲ z ≲ 7. Although no X-ray signal is detected in any stacked image, we obtain strong 3σ upper limits for the average black hole accretion rate (BHAR) as a function of star formation rate (SFR). At z ∼ 4 (5) where the stacked image is deeper, the 3σ BHAR upper limits per SFR are ∼1.5 (1.0) dex lower than the local black hole-to-stellar mass ratio, indicating that the SMBHs of SFGs in the inactive (BHAR $\lesssim 1 \, \mathrm{M}_\odot$ yr−1) phase are growing much more slowly than expected from simultaneous evolution. We obtain a similar result for BHAR per dark halo accretion rate. QSOs from the literature are found to have ∼1 dex higher SFRs and ≳ 2 dex higher BHARs than LBGs with the same dark halo mass. We also make a similar comparison for dusty starburst galaxies and quiescent galaxies from the literature. A duty-cycle corrected analysis shows that for a given dark halo, the SMBH mass increase in the QSO phase dominates over that in the much longer inactive phase. Finally, a comparison with the TNG300, TNG100, SIMBA100, and EAGLE100 simulations finds that they overshoot our BHAR upper limits by ≲ 1.5 dex, possibly implying that simulated SMBHs are too massive.

The properties of supermassive black holes and their host galaxies for type 1 and 2 active galactic nuclei in the eFEDS and COSMOS fields

In this study, our primary objective is to compare the properties of supermassive black holes (SMBHs) and their host galaxies between type 1 and type 2 active galactic nuclei (AGNs). In our analysis, we use X-ray detected sources in two fields, namely the eFEDS and the COSMOS-Legacy. To classify the X-ray sources, we performed a spectral energy distribution (SED) fitting analysis, using the CIGALE code. The robustness of our analysis was paramount so, to ensure this, we imposed stringent selection criteria. Thus, only sources with extensive photometric data across the optical, near- and mid-infrared part of the spectrum and reliable host galaxy properties and classifications were included. The final sample consists of 3312 AGNs, of which 3049 are classified as type 1 and 263 as type 2. The sources span a redshift range of 0.5 &lt; z &lt; 3.5 and encompass a wide range of X-ray luminosities, falling within 42 &lt; log,[LX,2−10 keV(erg s−1)] &lt; 46. Our results show that type 2 AGNs exhibit a tendency to inhabit more massive galaxies, by 0.2 − 0.3 dex (on a logarithmic scale), compared to type 1 sources. Type 2 AGNs also display, on average, lower specific black hole accretion rates, a proxy of the Eddington ratio, compared to type 1 AGNs. These differences persist across all redshifts and LX considered within our dataset. Moreover, our analysis uncovers that type 2 sources tend to have lower star formation rates compared to type 1 AGNs at z &lt; 1. This picture reverses at z &gt; 2 and log,[LX,2−10 keV(erg s−1)] &gt; 44. Similar patterns emerge when we categorize AGNs based on their X-ray obscuration levels (NH). However, in this case, the observed differences are pronounced only for low-to-intermediate LX AGNs and are also sensitive to the NH threshold applied for the AGN classification. These comprehensive findings enhance our understanding of the intricate relationships governing AGN types and their host galaxy properties across diverse cosmic epochs and luminosity regimes.

The link among X-ray spectral properties, AGN structure, and the host galaxy

In this work, we compare the supermassive black hole (SMBH) and host galaxy properties of X-ray obscured and unobscured AGN. For that purpose, we used ∼35 000 X-ray detected AGN in the 4XMM-DR11 catalogue for which there are available measurements for their X-ray spectral parameters, such as the hydrogen column density, NH, and photon index, Γ, from the XMM2Athena Horizon 2020 European project. We constructed the spectral energy distributions (SEDs) of the sources, and we calculated the host galaxy properties via SED fitting analysis, utilising the CIGALE code. We applied strict photometric requirements and quality selection criteria to include only sources with robust X-ray and SED fitting measurements. Our sample consists of 1443 AGN. In the first part of our analysis, we used different NH thresholds (1023 cm−2 or 1022 cm−2) while also taking into account the uncertainties associated with the NH measurements in order to classify these sources as obscured and unobscured (or mildly obscured). We find that obscured AGN tend to live in more massive systems (by ∼0.1 dex) that have a lower star-formation rate, SFR, (by ∼0.25 dex) compared to their unobscured counterparts. However, only the difference in stellar mass, M*, appears statistically significant (&gt; 2σ). The results do not depend on the NH threshold used to classify AGN. The differences in M* and SFR are not statistically significant for luminous AGN (log (LX,2−10 KeV/erg s−1) &gt; 44). Our findings also show that unobscured AGN have, on average, higher specific black hole accretion rates, λsBHAR, compared to their obscured counterparts, a parameter which is often used as a proxy of the Eddington ratio. In the second part of our analysis, we cross-matched the 1443 X-ray AGN with the SDSS DR16 quasar catalogue of Wu and Shen to obtain information on the SMBH properties of our sources. This resulted in 271 type 1 AGN at z &lt; 1.9. Our findings show that type 1 AGN with increased NH (&gt; 1022 cm−2) tend to have higher black hole masses, MBH, compared to AGN with lower NH values at similar M*. The MBH/M* ratio remains consistent for NH values below 1022 cm−2, but it exhibits signs of increasing at higher NH values. Finally, we detected a correlation between Γ and Eddington ratios, but only for type 1 sources with NH &lt; 1022 cm−2.

Cold quasar investigation: comparing star formation rates to black hole growth

ABSTRACT Cold quasars are a rare population of luminous, unobscured quasars associated with host galaxies that have a high star formation rate. We aimed to study the host galaxies of 64 of these cold quasars in order to probe how the supermassive black holes and host galaxies were co-evolving. We compiled data from the XMM-XXL survey and cross-matched with the VHS, WISE, and HerMES surveys to obtain multiwavelength photometry spanning the X-ray to the infrared and including optical spectroscopy. From the data, we calculated the supermassive black hole’s mass using broad emission from the Mg ii and Hbeta lines. We compared this with the stellar mass of the entire galaxy and found that the black holes are significantly more massive than would be predicted by local relations, indicating that the majority of black hole growth precedes the bulk of the the stellar mass formation. In addition to this, we created a spectral energy distribution for each galaxy to calculate the star formation rate. We compared the star formation rate with the black hole accretion rate and find that the stellar mass is rapidly increasing at a relative rate faster than the black hole growth, supporting the picture where the black hole grows first.

Probing supermassive black hole growth and its dependence on stellar mass and star formation rate in low-redshift galaxies

ABSTRACT We present an improved study of the relation between supermassive black hole growth and their host galaxy properties in the local Universe (z &amp;lt; 0.33). To this end, we build an extensive sample combining spectroscopic measurements of star formation rate (SFR) and stellar mass from Sloan Digital Sky Survey, with specific Black Hole accretion rate (sBHAR, $\lambda _{\mathrm{sBHAR}} \propto L_{\rm X}/\mathcal {M}_{\ast }$) derived from the XMM–Newton Serendipitous Source Catalogue (3XMM–DR8) and the Chandra Source Catalogue (CSC2.0). We find that the sBHAR probability distribution for both star-forming and quiescent galaxies has a power-law shape peaking at log λsBHAR ∼ −3.5 and declining towards lower sBHAR in all stellar mass ranges. This finding confirms the decrease of active galactic nucleus (AGN) activity in the local Universe compared to higher redshifts. We observe a significant correlation between $\log \, \lambda _{\rm sBHAR}$ and $\log \, {\rm SFR}$ in almost all stellar mass ranges, but the relation is shallower compared to higher redshifts, indicating a reduced availability of accreting material in the local Universe. At the same time, the BHAR-to-SFR ratio for star-forming galaxies strongly correlates with stellar mass, supporting the scenario where both AGN activity and stellar formation primarily depend on the stellar mass via fuelling by a common gas reservoir. Conversely, this ratio remains constant for quiescent galaxies, possibly indicating the existence of the different physical mechanisms responsible for AGN fuelling or different accretion mode in quiescent galaxies.