The X-ray-to-UV relation of active galactic nuclei (AGNs), commonly parametrized via the monochromatic luminosities at 2500 & Aring; and 2 keV, reflects the energetic interplay between the accretion disc and the X-ray-emitting corona, and is key for under-standing accretion physics. Previous studies suggest that disc-dominated emission becomes more prominent with increasing opti-cal luminosity. However, the redshift evolution of this relation remains debated, and a dependence on Eddington ratio, predicted by accretion flow models, is still observationally unconstrained. We revisit this relation using a large nearly all-sky sample by combin-ing the Sloan Digital Sky Survey Quasar Catalog Data Release 16 quasar (QSO) catalogue with X-ray data from XMM-Newton and the SRG/EROSITA All-Sky Survey Data Release 1, yielding 136 745 QSOs at redshifts 0.5 <= z < 3 We introduce a hierarchical Bayesian framework that treats X-ray detections and upper limits uniformly, enabling robust inference from both parametric and non-parametric models. We confirm a tight sublinear log L-x(2 keV)-log L-v (2500 A) correlation, but with a normalization at the lower end of previous estimates. Contrary to most literature results, we detect a mild but systematic redshift evolution: the relation flattens and its intrinsic scatter decreases at higher redshift. This trend is consistent with disc emission increasingly dominated by scattering and enhanced energy transfer to the X-ray corona, potentially indicating redshift evolution in the X-ray bolometric correction. We find no significant dependence on Eddington ratio, in tension with recent accretion flow models.

An important diagnostic of the inner structure of accretion flows onto supermassive black holes are the stochastic flux variations at X-ray wavelengths. Despite its significance, a systematic characterization of the statistical properties of the X-ray variability to the highest Eddington ratios and most massive black holes is still lacking. In this paper, we address this issue using SRG/eROSITA (Spectrum-Roentgen-Gamma/extended ROentgen Survey with an Imaging Telescope Array) five-epoch light curves to characterize the mean X-ray variability of optically selected Sloan Digital Sky Survey quasars extending to black holes masses of 10(10) solar and accretion rates close to the Eddington limit. The adopted variability statistic is the ensemble normalized excess variance, which is measured using a novel hierarchical Bayesian model (EBEXVAR) tailored to the Poisson nature of the X-ray light curves. We find a clear anticorrelation of the ensemble variability with black hole mass, extending previous results to time-scales of months. This can be interpreted as evidence for an X-ray corona size and/or physical conditions that scale with black holes mass. We also find an unexpected increase of the ensemble normalized excess variance close to the Eddington limit, which is contrary to the predictions of empirical variability models. This result suggests an additional variability component for fast growing black holes that may be related to systematic variations of the hot corona size with Eddington ratio or shielding of the hot corona by an inner puffed-up disc and/or outflows.

Context. Active galactic nuclei (AGN) stand as extreme X-ray emitters where disk-corona interplay shapes their spectral energy distribution. The soft X-ray excess, a unique feature of AGN in the 0.5 - 2.0 keV, encodes critical information on the "warm corona" structure bridging the disk and hot corona. However, the systematic evolution of this feature with fundamental accretion parameters in large AGN samples - particularly those studied through the spectral stacking technique - remains observationally unconstrained. Aims. The eROSITA All-Sky Survey (eRASS:5) provides an unprecedented sample to statistically map AGN spectral properties. We present a multiwavelength investigation of how the average AGN X-ray spectra evolve with accretion parameters (alpha(ox), L-UV, lambda(Edd), M-BH), and we explore the disk-corona connection by further combining stacked UV data. Methods. We have developed Xstack, a novel X-ray spectral stacking code that consistently stacks rest-frame pulse invariant (PI) spectra and associated responses using optimized response weighting to preserve spectral shapes. With Xstack, we stacked 17 929 AGNs ("spec-z" sample, total exposure similar to 23 Ms) with similar X-ray loudness, alpha(ox), and UV luminosity, L-UV, and 4159 AGNs ("BH-mass" sample, similar to 3 Ms) with similar Eddington ratios, lambda(Edd), and black hole masses, M-BH. We analyzed the resulting stacked X-ray spectra with a phenomenological model for both samples. We further fit the stacked optical-UV X-ray SED with the physical AGNSED model on a 3 x 3 M-BH - lambda(Edd) grid. Results. We observed that the soft excess strength rises strongly with increasing alpha(ox) and lambda(Edd) binning (by a factor of five), while the hard X-ray spectral shape remains largely unchanged, consistent with the interpretation that soft excess is primarily driven by the warm corona rather than reflection. The trends are weaker with L-UV binning and reversed for M-BH binning. The analysis of the optical-UV X-ray SEDs with AGNSED revealed that the warm corona radius (in units of R-g) generally increases with lambda(Edd) and decreases with M-BH, or equivalently the disk-to-warm-corona transition consistently occurs near similar to 1 x 10(4) K. The hot corona contracts with lambda(Edd), and the radius remains independent of M-BH, aligning with disk evaporation predictions. Conclusions. The soft excess is likely warm-corona dominated, with the disk-to-warm-corona transition potentially linked to hydrogen ionization instability at similar to 1 x 10(4) K, which is consistent with previous work utilizing eFEDS-HSC stacked data. Our work highlights the power of spectral stacking for revealing the AGN disk-corona connection.

Context. During its calibration and performance verification phase, the eROSITA instrument aboard the Spectrum-RG satellite performed a uniform wide-area X-ray survey of approximately 140 deg(2), known as the eROSITA Final Equatorial Depth Survey (eFEDS). Aims. The primary aim of eFEDS is to demonstrate the scientific performance to be expected at the end of the eight-pass eROSITA all-sky survey. This survey will provide the first focussed image of the whole sky in the hard X-ray (>2 keV) bandpass. The expected source population in this energy range is thus of great interest, particularly for AGN studies. Methods. We used a 2.3-5 keV selection to construct a sample of 246 point-like hard X-ray sources for further study and characterisation. These sources are classified as either extragalactic (similar to 90%) or Galactic (similar to 10%), with the former consisting overwhelmingly of AGN and the latter active stars. We concentrated our further analysis on the extragalactic AGN sample, describing their X-ray and multi-wavelength properties and comparing them to the eFEDS main AGN sample selected in the softer 0.2-2.3 keV band. Results. The eROSITA hard band selects a subsample of sources that is a factor of more than ten brighter than the eFEDS main sample. The AGN within the hard population reach up to z = 3.2 but on the whole, they are relatively nearby, with median z = 0.34 compared to z = 0.94 for the main sample. The hard survey probes typical luminosities in the range log L-X = 43-46. The X-ray spectral analysis shows significant intrinsic absorption (with log N-H > 21) in similar to 20% of the sources, with a hard X-ray power law continuum with mean < Gamma >= 1.83 +/- 0.04, which is typical of AGN, but slightly harder than the soft-selected eROSITA sample. Around 10% of the hard sample show a significant 'soft excess' component. The sampled black hole mass distribution in the eFEDS broad-line AGN population is consistent with that of the deeper COSMOS survey that probes a higher redshift population. On the other hand, the Eddington ratios appear systematically lower, which is consistent with the idea that the decline in SMBH activity since z similar to 1 is due to a reduction in the typical accretion rate, rather than a shift towards activity in lower-mass black holes. Conclusions. The eFEDS hard sample provides a preview of what can be expected from the eRASS final survey in terms of data quality. This pilot survey indicates the power of eROSITA to shed new light on the demographics and evolution of AGN, and the potential for discovery of new and rare populations.

Context. AGNs are strong X-ray emitters shaped by disk-corona interactions. The soft excess (0.5-2.0 keV) reveals key information about the "warm corona" bridging the disk and hot corona. Yet, how this feature evolves with accretion properties remains poorly constrained, especially in large samples using spectral stacking. Aims. The eROSITA All-Sky Survey (eRASS:5) provides an unprecedented sample. We investigate how the average AGN X-ray spectra evolve with accretion parameters, and explore disk-corona connection by further combining stacked UV data. Methods. We developed Xstack, a novel tool that stacks rest-frame X-ray spectra and responses while preserving spectral shape through optimized weighting. We stack 17929 AGNs ("spec-z" sample, 23 Ms) with similar X-ray loudness alpha_ox, UV luminosity L_UV, and 4159 AGNs ("BH-mass" sample, 3 Ms) with similar Eddington ratio lambda_Edd and black hole mass M_BH. The resulting stacked X-ray spectra are analyzed with a phenomenological model. We further fit the stacked optical-UV-Xray SED with AGNSED model. Results. Soft excess strengthens strongly with alpha_ox and lambda_Edd (~5), while the hard X-ray spectral shape remains largely unchanged, supporting that soft excess is dominated by warm corona rather than reflection. AGNSED modeling reveals that warm corona radius (R_g units) generally increases with lambda_Edd and decreases with M_BH, or equivalently the disk-to-warm-corona transition consistently occurs near 1e4 K. The hot corona contracts with lambda_Edd and is unaffected by M_BH, aligning with disk evaporation predictions. Conclusions. The soft excess likely originates from a warm corona, with the disk to warm corona transition tied to hydrogen ionization near 1e4 K - supporting earlier eFEDS-HSC stacking results (Hagen et al. 2024). This study shows the strength of spectral stacking in probing AGN disk-corona physics.

With its first All-Sky Survey (eRASS1), the extended ROentgen Survey with an Imaging Telescope Array (eROSITA) on board the Spectrum-Roentgen-Gamma (SRG) mission has offered an unprecedented, comprehensive view of the variable X-ray sky. Featuring enhanced sensitivity, broader energy coverage, and improved resolution compared to prior surveys, the eRASS1 Data Release 1 (DR1) catalogue underwent a variability analysis, and in this paper, we performed an advanced variability analysis focusing on a substantial subset of 128 669 sources, all exhibiting a net count exceeding ten. We performed multiple variability tests, utilising conventional normalised excess variance (NEV), maximum amplitude variability (AMP), and Bayesian excess variance methods (bexvar). The analysis focused on binned light curves; specifically, employing one eroday (a great circle scan with a duration of 4 hours) binning of the German part of the first eROSITA all-sky survey (eROSITA-DE) data, i.e., the source sample covers only half of the sky. Within the 128 669 DR1 sources with light curves, our research pinpointed 808 light curves that show hints of variability according to the AMP test, and 298 according to the NEV test. However, after applying suitable thresholds, 90 (123) sources were found to be significantly variable according to the AMP (NEV) tests. In addition, 1342 sources are considered variable according to the Bayesian test bexvar. The total number of unique sources is 1709, and they form the catalogue of variable sources released with this paper. We cross-matched with existing X-ray catalogues and identified 258, 318, 598, and 120 sources in 4XMM DR13, 2SXPS, 2RXS, and CSC2.1, respectively. Only 27 sources overlap across all catalogues, while 882 are new X-ray detections from eROSITA DR1. About 70% are coronal stars, 5% are Quasi-Stellar Objects, and 1.6% are normal galaxies. We further subclassified 18 sources as LMXBs, 11 as HMXBs, and 14 as bright stars. In this paper, we analyse the variability of eRASS1 sources on a timescale of only a few days. To study the physics of variable sources, we need more deeply pointed observations with other X-ray missions or at least the final depth of the eRASS: 8 observations. The timescale of the eRASS1 observations is not representative of the timescales of the expected upcoming eRASS catalogues. A substantial 52% of the eRASS1 variable sources were first discovered with eROSITA. The DR1 variability catalogue is excellent for follow-up observations with telescopes such as XMM-Newton , Chandra , or Swift .

We present one of the largest uniform optical spectroscopic surveys of X-ray selected sources to date that were observed as a pilot study for the Black Hole Mapper (BHM) survey. The BHM program of the Sloan Digital Sky Survey (SDSS)-V is designed to provide optical spectra for hundreds of thousands of X-ray selected sources from the SRG/eROSITA all-sky survey. This significantly improves our ability to classify and characterise the physical properties of large statistical populations of X-ray emitting objects. Our sample consists of 13 079 sources in the eROSITA eFEDS performance verification field, 12 011 of which provide reliable redshifts from 0 less than or similar to z <= 5.8. The vast majority of these objects were detected as point-like sources (X-ray flux limit F0.5 - 2 keV greater than or similar to 6.5 x 10(-15) erg/s/cm(2)) and were observed for about 20 years with fibre-fed SDSS spectrographs. After including all available redshift information for the eFEDS sources from the dedicated SDSS-V plate programme and archival data, we visually inspected the SDSS optical spectra to verify the reliability of these redshift measurements and the performance of the SDSS pipeline. The visual inspection allowed us to recover reliable redshifts (for 99% of the spectra with a signal-to-noise ratio of > 2) and to assign classes to the sources, and we confirm that the vast majority of our sample consists of active galactic nuclei (AGNs). Only similar to 3% of the eFEDS/SDSS sources are Galactic objects. We analysed the completeness and purity of the spectroscopic redshift catalogue, in which the spectroscopic completeness increases from 48% (full sample) to 81% for a cleaner, brighter (r(AB) < 21.38) sample that we defined by considering a high X-ray detection likelihood, a reliable counterpart association, and an optimal sky coverage. We also show the diversity of the optical spectra of the X-ray selected AGNs and provide spectral stacks with a high signal-to-noise ratio in various sub-samples with different redshift and optical broad-band colours. Our AGN sample contains optical spectra of (broad-line) quasars, narrow-line galaxies, and optically passive galaxies. It is considerably diverse in its colours and in its levels of nuclear obscuration.

Aims. The interplay between the accretion of supermassive black holes (SMBHs) and the stellar mass growth of the host galaxies is still a matter of hot debate. The accretion of the central SMBHs is expected to release energy under the form of active galactic nuclei. This energy is believed to impact the star formation activity and contribute to the quenching of the host galaxies. Here, we address this key unsolved issue with our cosmological semi-empirical model DECODE (Discrete statistical sEmi-empiriCal mODEl). Methods. In DECODE , we grow galaxies with their star formation rate linked to halo accretion rate distributions via abundance matching. SMBHs are evolved following the stellar mass growth of their host galaxies by assigning an accretion rate at each redshift from the empirical Eddington ratio distributions and duty cycles. We tested the assumption that galaxies permanently quench when their central SMBHs approach the limit imposed by the observed M BH − σ ★ relation, as a proxy of SMBH disruptive feedback. Results. We find that simply imposing the M BH − σ ★ condition is sufficient to generate a fraction of quenched galaxies consistent with current data, including the newest ones from Euclid. In addition, our minimal data-driven model also predicts SMBH scaling relations consistent in slope and normalisation with those that have been observed, and an M BH − M ★ relation weakly evolving with redshift. The model also naturally generates SMBH accretion rates peaking within 1 Gyr of their host star formation histories. Interestingly, we note that all the main predictions on galaxy quenched fractions and SMBH growth histories and scaling relations are degenerate with those expected in a halo quenching model. Conclusions. The comprehensive data-driven model presented in this work represents an invaluable tool to investigate SMBH demography across time and environments in an accurate, physically motivated manner, ideally suited to rapidly exploring the implications from large surveys, such as Euclid and Rubin-LSST.

ABSTRACT We present an improved study of the relation between supermassive black hole growth and their host galaxy properties in the local Universe (z < 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.

ABSTRACT In active galactic nuclei, the relationship between UV and X-ray luminosity is well studied (often characterised by αox) but often with heterogeneous samples. We have parametrized the intrinsic distribution of X-ray luminosity, LX, for the optically selected sample of SDSS quasars in the Stripe 82 and XXL fields across redshifts 0.5–3.5. We make use of the available XMM observations and a custom pipeline to produce Bayesian sensitivity curves that are used to derive the intrinsic X-ray distribution in a hierarchical Bayesian framework. We find that the X-ray luminosity distribution is well described by a Gaussian function in log10 LX space with a mean that is dependent on the monochromatic 2500 Å UV luminosity, L2500. We also observe some redshift dependence of the distribution. The mean of the LX distribution increases with redshift while the width decreases. This weak but significant redshift dependence leads to L2500–LX and L2500–αox relations that evolve with redshift, and we produce a redshift- and L2500-dependent αox equation. Neither black hole mass nor Eddington ratio appear to be potential drivers of the redshift evolution.

ABSTRACT We adopt a Bayesian X-ray spectral approach to investigate the accretion properties of unobscured ($20\lt \log (N_{\rm H}/{\rm cm}^{-2}\lt 22$) and obscured ($22\lt \log (N_{\rm H}/{\rm cm}^{-2}\lt 24$) active galactic nuclei (AGNs) to shed light on the orientation versus evolution scenarios for the origin of the obscuring material. For a sample of 3882 X-ray-selected AGN from the Chandra COSMOS Legacy, AEGIS, and CDFS extragalactic surveys, we constrain their stellar masses, $M_\star$, intrinsic X-ray luminosities, $L_{\rm X}$, obscuring column densities, $N_{\rm H}$, and specific accretion rates $\lambda \propto L_{\rm X}/M_\star$. By combining these observables within a Bayesian non-parametric approach, we infer, for the first time, the specific accretion rate distribution (SARD) of obscured and unobscured AGN to $z\approx 3$, i.e. the probability of a galaxy with mass $M_\star$ at redshift z hosting an AGN with column density $N_{\rm H}$ and specific accretion rate $\lambda$. Our findings indicate that (1) both obscured and unobscured SARDs share similar shapes, shifting towards higher accretion rates with redshift, (2) unobscured SARDs exhibit a systematic offset towards higher $\lambda$ compared to obscured SARD for all redshift intervals, (3) the obscured AGN fraction declines sharply at $\log \lambda _{\rm break} \sim -2$ for $z \lt 0.5$, but shifts to higher $\lambda$ values with increasing redshift, (4) the incidence of AGN within the theoretically unstable blow-out region of the $\lambda -N_{\rm H}$ plane increases with redshift. These observations provide compelling evidence for AGN ‘downsizing’ and radiation-regulated nuclear-scale obscuration with an increasing host galaxy contribution towards higher redshifts.

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 & ouml;tes field, we measure BH-galaxy assembly in massive galaxies (M star greater than or similar to 10(10)M(circle dot)) 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 >= 10(10)M(circle dot), 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.

The recent serendipitous discovery of a new population of short duration X-ray transients, thought to be associated with collisions of compact objects or stellar explosions in distant galaxies, has motivated efforts to build up statistical samples by mining X-ray telescope archives. Most searches to date however, do not fully exploit recent developments in the signal and imaging processing research domains to optimise searches for short X-ray flashes. This paper addresses this issue by presenting a new source detection pipeline, STATiX (Space and Time Algorithm for Transients in X-rays), which directly operates on 3-dimensional X-ray data cubes consisting of two spatial and one temporal dimension. The algorithm leverages wavelet transforms and the principles of sparsity to denoise X-ray observations and then detect source candidates on the denoised data cubes. The light curves of the detected sources are then characterised using the Bayesian blocks algorithm to identify flaring periods. We describe the implementation of STATiX in the case of XMM-Newton data, present extensive validation and performance verification tests based on simulations and also apply the pipeline to a small subset of seven XMM-Newton observations, which are known to contain transients sources. In addition to known flares in the selected fields we report a previously unknown short duration transient found by our algorithm that is likely associated with a flaring Galactic star. This discovery demonstrates the potential of applying STATiX to the full XMM-Newton archive.

Finding massive black holes (MBHs,MBH≈ 104–107M⊙) in the nuclei of low-mass galaxies $\left( {{M_*}\mathop {\mathop < \limits_ }\limits_ {{10}^{10}}{M_ \odot }} \right)$ is crucial to constrain seeding and growth of black holes over cosmic time, but it is particularly challenging due to their low accretion luminosities. Variability selection via long-term photometric ultraviolet, optical, or infrared (UVOIR) light curves has proved effective and identifies lower-Eddington ratios compared to broad and narrow optical spectral lines searches. In the inefficient accretion regime, X-ray and radio searches are effective, but they have been limited to small samples. Therefore, differences between selection techniques have remained uncertain. Here, we present the first large systematic investigation of the X-ray properties of a sample of known MBH candidates in dwarf galaxies. We extracted X-ray photometry and spectra of a sample of ~200 UVOIR variability-selected MBHs and significantly detected 17 of them in the deepest available SRG/eROSITA image, of which four are newly discovered X-ray sources and two are new secure MBHs. This implies that tens to hundreds of LSST MBHs will have SRG/eROSITA counterparts, depending on the seeding model adopted. Surprisingly, the stacked X-ray images of the many non-detected MBHs are incompatible with standard disk-corona relations, typical of active galactic nuclei, inferred from both the optical and radio fluxes. They are instead compatible with the X-ray emission predicted for normal galaxies. After careful consideration of potential biases, we identified that this X-ray weakness needs a physical origin. A possibility is that a canonical X-ray corona might be lacking in the majority of this population of UVOIR-variability selected low-mass galaxies or that unusual accretion modes and spectral energy distributions are in place for MBHs in dwarf galaxies. This result reveals the potential for severe biases in occupation fractions derived from data from only one waveband combined with SEDs and scaling relations of more massive black holes and galaxies.

The assembly and co-evolution of supermassive black holes (SMBH) and their host galaxy stellar population is a key open questions in galaxy evolution. Stellar mass ($M_\star$) and star formation rate (SFR), are inferred by modeling the spectral energy distribution (SED). For galaxies triggering SMBH activity, the active galactic nucleus (AGN) contaminates the light at all wavelengths, hampering the inference of galaxy parameters. Incomplete AGN templates can lead to systematic overestimates of the stellar mass, biasing our understanding of AGN-galaxy co-evolution. This challenge has gained further impetus with the advent of sensitive wide-area surveys with millions of luminous AGN, including by eROSITA, Euclid and LSST. We aim to estimate the accuracy and bias of AGN host galaxy parameters and improve upon existing techniques. This work makes two contributions: 1) a new SED fitting code, GRAHSP, with a flexible, empirically motivated AGN model including a power law continuum emission lines, a FeII forest and a flexible infrared torus. We verify that our model reproduces published X-ray to infrared SEDs of AGN to better than 20\% accuracy. A fully Bayesian fit with nested sampling includes uncertainties in the model and the data, making the inference highly robust. 2) we created a benchmark photometric dataset where pure quasars are merged with non-AGN pure galaxies into a hybrid (Chimera) object but with known galaxy and AGN properties. Comparing the true and retrieved $M_\star$, SFR and AGN luminosities shows that previous codes systematically over-estimate $M_\star$ and SFR by 0.5 dex with a wide scatter of 0.7 dex, at AGN luminosities above 10^44 erg/s. In contrast, GRAHSP shows no bias on $M_\star$ and SFR. GRAHSP also estimates more realistic uncertainties. GRAHSP enables characterization of the environmental conditions conducive to black hole growth. (abridged)

The eROSITA telescope array aboard the Spektrum Roentgen Gamma (SRG) satellite began surveying the sky in December 2019, with the aim of producing all-sky X-ray source lists and sky maps of an unprecedented depth. Here we present catalogues of both point-like and extended sources using the data acquired in the first six months of survey operations (eRASS1; completed June 2020) over the half sky whose proprietary data rights lie with the German eROSITA Consortium. We describe the observation process, the data analysis pipelines, and the characteristics of the X-ray sources. With nearly 930 000 entries detected in the most sensitive 0.2-2.3 keV energy range, the eRASS1 main catalogue presented here increases the number of known X-ray sources in the published literature by more than 60%, and provides a comprehensive inventory of all classes of X-ray celestial objects, covering a wide range of physical processes. A smaller catalogue of 5466 sources detected in the less sensitive but harder 2.3-5 keV band is the result of the first true imaging survey of the entire sky above 2 keV. We present methods to identify and flag potential spurious sources in the catalogues, which we applied for this work, and we tested and validated the astrometric accuracy via cross-comparison with other X-ray and multi-wavelength catalogues. We show that the number counts of X-ray sources in eRASS1 are consistent with those derived over narrower fields by past X-ray surveys of a similar depth, and we explore the number counts variation as a function of the location in the sky. Adopting a uniform all-sky flux limit (at 50% completeness) of F0.5-2 keV > 5 x 10(-14) erg s(-1) cm(-2), we estimate that the eROSITA all-sky survey resolves into individual sources about 20% of the cosmic X-ray background in the 1-2 keV range. The catalogues presented here form part of the first data release (DR1) of the SRG/eROSITA all-sky survey. Beyond the X-ray catalogues, DR1 contains all detected and calibrated event files, source products (light curves and spectra), and all-sky maps. Illustrative examples of these are provided.

Although flux variability is one of the defining properties of accretion flows onto supermassive black holes, its dependence on physical parameters such as the mass of the compact object and the Eddington ratio remain under discussion. In this paper we address this issue using the structure function statistic to measure the variability at X-ray wavelengths of a sample of optically selected QSOs with available black hole masses and Eddington ratios. We present a new Bayesian methodology for estimating the structure function tailored to the Poisson nature of the X-ray data. This is applied to 15,548 SDSS DRQ16 QSOs with repeat observations in the XMM-Newton archive and/or the SRG/eROSITA All Sky Survey. The X-ray structure function monotonically increases to time intervals of about 10-15 years, suggesting a preference for scenarios in which instabilities of the accretion disk contribute to the X-ray variability on long timescales. Additionally, there is evidence that the amplitude of the stochastic X-ray flux variations rises with decreasing black hole mass and Eddington ratio. This finding imposes stringent constraints on empirical models of Active Galactic Nuclei variability derived from local samples, emphasizing the significance of high-redshift population studies for comprehending the stochastic flux variations in active black holes.

Radio jets are present in a diverse sample of AGN. However, the mechanisms of jet powering are not fully understood, and it is yet unclear to what extent they obey mass-invariant scaling relations, similar to those found for the triggering and fuelling of X-ray selected AGN. We study the incidence of eROSITA/eFEDS X-ray and LOFAR radio AGN as a function of several stellar mass normalised AGN power indicators. A new sample of radio AGN from the LOFAR-eFEDS survey is defined and we publicly release this catalogue, including host galaxy counterparts from the Legacy Survey DR9, LOFAR radio morphologies and host galaxy properties from the complete, spectroscopic (z<0.4) GAMA09 survey. The fraction of GAMA09 galaxies hosting radio, X-ray and both radio and X-ray AGN are calculated as a function of the specific black hole kinetic ($\lambda_{\rm Jet}$) and radiative ($\lambda_{\rm Edd}$) power. The incidence of eFEDS X-ray AGN as a function of $\lambda_{\rm Edd}$ shows the same mass-invariance as found in past studies. Meanwhile, radio AGN, regardless of their morphology, are more likely to be hosted in more massive galaxies, at all $\lambda_{\rm Jet}$. Across the stellar mass range, the compact radio AGN incidence follows the same power-law distribution, showing that it is not only high mass galaxies that host high power radio AGN and vice versa. On the other hand, the incidence of compact and complex radio AGN is boosted at the highest jet powers, diverging from a simple power-law. Interestingly, this increased incidence cannot be explained by more powerful radio AGN lying in more dense environments which could naturally boost their radio luminosity. Overall, we show that statistical incidence studies are a powerful method to probe disk-jet coupling for different AGN accretion modes, although future work on a more reliable determination of jet power for diverse samples of radio AGN is needed.

Context. The soft X-ray instrument eROSITA on board the Spectrum-Roentgen-Gamma (SRG) observatory has successfully completed four of the eight planned all-sky surveys, detecting almost one million X-ray sources during the first survey (eRASS1). The catalog of this survey will be released as part of the first eROSITA data release (DR1). Aims. Based on X-ray aperture photometry, we provide flux upper limits for eRASS1 in several energy bands. We cover galactic longitudes between 180 degrees <= l <= 360 degrees (eROSITA-DE). These data are crucial for studying the X-ray properties of variable and transient objects, as well as non-detected sources in the eROSITA all-sky survey data. Methods. We performed aperture photometry on every pixel of the SRG/eROSITA standard pipeline data products for all available sky tiles in the single detection band (0.2-2.3 keV). Simultaneously, we performed the same analysis in the three-band detection at soft (0.2-0.6 keV), medium (0.6-2.3 keV), and hard (2.3-5.0 keV) energy bands. Based on the combination of products for the individual bands, we are also able to provide aperture photometry products and flux upper limits for the 0.2-5.0 keV energy band. The upper limits were calculated based on a Bayesian approach that utilizes detected counts and background within the circular aperture. Results. The final data products consist of tables with the aperture photometry products (detected counts, background counts, and exposure time), a close-neighbor flag, and the upper flux limit based on an absorbed power-law spectral model (Gamma = 2.0, N-H = 3 x 10(20) cm(-2)). The upper limits are calculated using the one-sided 3 sigma confidence interval (CL) of a normal distribution, representing CL = 99.87%. The aperture photometry products allow for an easy computation of upper limits at any other confidence interval and spectral model. These data are stored in a database with hierarchical indexes in order to offer a fast query option. Conclusions. We provide a detailed description of the process of retrieving SRG/eROSITA upper limits for a large set of input positions, as well as of the eROSITA data, the X-ray aperture photometry, the upper limit calculation, and the final data products. The eROSITA upper flux limits represent either the maximum flux of potential non-detections or the 3 sigma upper flux uncertainty of detected sources. We emphasize the importance of choosing the right spectral model that ought to match the spectral shape of the source of interest: the wrong spectral model can produce discrepancies of up to 30% in the final flux upper limit value. We also describe the architecture of the database and the web tool, which are designed to handle large queries on input positions.

ABSTRACT Environmental effects are believed to play an important yet poorly understood role in triggering accretion events onto the supermassive black holes (SMBHs) of galaxies (active galactic nuclei; AGNs). Massive clusters, which represent the densest structures in the Universe, provide an excellent laboratory to isolate environmental effects and study their impact on black hole growth. In this work, we critically review observational evidence for the preferential activation of SMBHs in the outskirts of galaxy clusters. We develop a semi-empirical model under the assumption that the incidence of AGN in galaxies is independent of environment. We demonstrate that the model is broadly consistent with recent observations on the AGN halo occupation at z = 0.2, although it may overpredict satellite AGN in massive haloes at that low redshift. We then use this model to interpret the projected radial distribution of X-ray sources around high redshift (z ≈ 1) massive ($\gt 5 \times 10^{14} \, M_\odot$) clusters, which show excess counts outside their virial radius. Such an excess naturally arises in our model as a result of sample variance. Up to 20 per cent of the simulated projected radial distributions show excess counts similar to the observations, which are however, because of background/foreground AGN and hence, not physically associated with the cluster. Our analysis emphasizes the importance of projection effects and shows that current observations of z ≈ 1 clusters remain inconclusive on the activation of SMBHs during infall.