Open Access
Issue
A&A
Volume 691, November 2024
Article Number A30
Number of page(s) 19
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202451035
Published online 29 October 2024

© The Authors 2024

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model.

Open Access funding provided by Max Planck Society.

1. Introduction

The study of accreting super-massive black holes (SMBHs) in the first gigayear after the Big Bang reshapes our understanding of black hole formation and growth, with potentially profound implications on the evolution of their host galaxies (e.g., Habouzit et al. 2022). With the advent of deep infrared (IR) space-based imaging, as well as multi-object and slitless spectroscopy with the James Webb Space Telescope (JWST), the redshift frontier of accreting SMBHs has been pushed far into cosmic dawn, to z > 8 (e.g., Harikane et al. 2023; Larson et al. 2023; Bunker et al. 2023; Maiolino et al. 2024). Combining low-resolution prism spectroscopy with JWST/NIRSpec and with the high-spatial-resolution of the Chandra X-ray space telescope, Goulding et al. (2023) have recently discovered a (lensed) SMBH-powered active galactic nucleus (AGN) at z = 10.1, corresponding to a lookback-time of ∼13.2 Gyr. The most extreme SMBHs arguably reside in the cores of quasars, luminous persistent panchromatic beacons, that are radiating up to Lbol ∼ 1048 erg s−1 at cosmic times reaching back deep into the epoch of reionization. While the less massive SMBHs discovered by JWST may grant us access to crisp snapshots of earlier evolutionary stages of AGN and their hosts, the mere existence of matured quasars and their central ∼109M SMBHs in the first gigayear of the universe (up to z ∼ 7.5, e.g., Bañados et al. 2018; Wang et al. 2021a; Yang et al. 2021) poses a puzzling progenitor problem: How did the seeds of black holes evolve to such massive systems in the cosmic timespan of a few hundred million years? Even under the unlikely assumption of sustained Eddington-limited accretion, well-understood stellar seeding models based on the collapse of early Pop III stars struggle to explain the formation of the most massive SMBHs found in the cores of z > 6.5 quasars (for reviews see Inayoshi et al. 2020; Fan et al. 2023). Due to the steep decline of the AGN population with luminosity at high redshift, quasars with Lbol ∼ 1048 erg s−1 are rare in the first gigayear (Jiang et al. 2016; Vito et al. 2018; Kulkarni et al. 2019; Matsuoka et al. 2023; Schindler et al. 2023). The low space-density (∼1 Gpc−3) of these objects implies that wide-field surveys, sampling large cosmological volumes, are needed for their search. The optical and near-IR (NIR) Euclid Wide Survey will enable the discovery of numerous z > 7.5 quasars (Euclid Collaboration 2019) in the coming years.

At slightly lower redshifts, the most luminous quasars at z ∼ 6 impose less stringent constraints on black hole seed models but set a strict lower limit on the black hole population in the early universe and their accretion history (Vito et al. 2018; Wolf et al. 2021, 2023; Barlow-Hall et al. 2023). These rare AGN populating the extreme end of the luminosity and mass scale constitute powerful demographic probes (e.g., Belladitta et al. 2020). So far, about ∼475 reionisation-era quasars have been discovered in wide optical and IR surveys (Fan et al. 2023, and references therein). Selection methods are based on the blue-band dropout photometry imprinted by neutral hydrogen absorption by the intergalactic medium (IGM, Madau & Rees 2000). Quasar searches in imaging data are challenging as their optical and IR colours are similar to high-surface-density Galactic contaminants: late-type stars (M-, L- and T-dwarfs, e.g., Caballero et al. 2008). M-dwarfs, the most numerous stars in the galaxy (Chabrier 2001), emit in the X-rays due to coronal activity and flares (e.g., Briggs & Pye 2004; Rutledge et al. 2000; Stelzer et al. 2013; Magaudda et al. 2022). However, at z > 5.7, the bulk of galactic colour-contaminants has spectral type SpT M7 or later (Yang et al. 2017). These objects, dubbed ultra-cool dwarfs (UCDs), are intrinsically X-ray faint in their quiescent state. Stelzer et al. (2022) have investigated the X-ray emission of UCDs in the first pass all-sky survey (eRASS11, Merloni et al. 2024) of the extended Roentgen Survey with an Imaging Telescope Array (eROSITA, Predehl et al. 2021) onboard the Spektrum-Roentgen-Gamma (SRG) mission (Sunyaev et al. 2021). Out of 4754 spectroscopically confirmed UCDs (14 914 candidates) in their sample, 21 (37) are associated with an eRASS1 source, that is, less than 0.3% of the input catalogue. Stelzer et al. (2022) show that the eRASS1-detected UCDs have large X-ray to bolometric luminosity ratios (log LX/Lbol ∼ −3), indicative of flaring activity. They conclude that the eRASS1 survey is only sensitive to flaring or the closest quiescent UCDs (< 20 pc). We note that a volume complete X-ray survey of nearby M-dwarfs (< 10 pc) based on eRASS:4 was presented by Caramazza et al. (2023).

Thus, the detection of X-rays in eROSITA can significantly reduce the degeneracy between UCDs and high-redshift quasars using the latter’s ubiquitous, strong coronal and jet X-ray emission. In this work, we present a pilot z ∼ 6 quasar discovery program in eRASS1 through eRASS:4. Wolf et al. (2021) have already demonstrated the sensitivity of eROSITA to the X-ray emission of z ∼ 6 quasars: we expect eROSITA to detect ∼100 of the most X-ray luminous quasars (log LX > 45.5), at the end of the mission (eRASS:8). In this luminosity-redshift regime, the strong X-ray emission arises from a range of high-energetic physical effects including Doppler boosting of non-thermal jet emission (e.g., Belladitta et al. 2020), inverse Comptonization of cosmic microwave background photons (iC-CMB, e.g., Ighina et al. 2019; Medvedev et al. 2020) and rapid accretion (e.g., Wolf et al. 2023; Zappacosta et al. 2023). eROSITA thus probes a diversity of X-ray boosting mechanisms whose contributions are yet to be disentangled (e.g., Ighina et al. 2021). Our selection pipeline combines optical, IR and X-ray information. Similar panchromatic searches, spanning multiple orders of magnitudes in frequency, have recently been successfully performed by Gloudemans et al. (2022), Ighina et al. (2023). In these works, the authors based their selection on the detection of jets at radio frequencies with the LOw-Frequency ARray (LOFAR, van Haarlem et al. 2013) and the Australian square kilometre array pathfinder (ASKAP, Hotan et al. 2021) instead of accretion and jet signatures in X-rays. We report the discovery of 5 new z > 5.6 quasars in the joint Dark Energy Survey data release 2 (DES DR22, Dark Energy Survey Collaboration 2016) footprint, two of which are robustly detected in the eRASS catalogues. The newly discovered eRASS quasars are rare and highly energetic phenomenons at high redshift, which we showcase in this work: an X-ray variable broad absorption line (BAL) quasar and a blazar, a quasar for which the jet axis closely aligns with the line of sight.

This paper is organised as follows. In Section 2 we present the optical, IR, and X-ray based quasar selection pipeline. Discoveries through spectroscopic confirmation are presented in Section 3. In Section 4, we present the properties of the two newly discovered eRASS quasars. We discuss the implication of these discoveries and present our conclusions in Section 5. Throughout this work, magnitudes are given in the AB reference system. Unless specified otherwise, uncertainties are given at the 1σ level. Distance calculations assume a flat ΛCDM cosmology (Planck Collaboration VI 2020, H0 = 67.66 ± 0.42 km/(Mpc s), Ωm = 0.3111 ± 0.0056, ΩΛ = 0.6889 ± 0.0056, ΩK = 0.001 ± 0.002).

2. Panchromatic selection pipeline

2.1. Candidate selection in optical/IR

We have selected z > 5.6 quasar candidates using optical data from DES. In addition, NIR and mid-IR (MIR) photometry from the VISTA Hemisphere Survey (VHS, McMahon et al. 2013) Data Release 5 and CatWISE2020 (Marocco et al. 2021) complement the DES data in our selection pipeline. The combination of optical and infrared broad-band photometry covers a significant portion of the z > 5.6 quasar restframe UV/optical SED, and its efficiency for quasar discovery has previously been demonstrated (e.g., Reed et al. 2019).

2.1.1. Colour pre-selection

The main DES source catalogue contains close to 700 million sources. To increase the computational efficiency of our selection pipeline, we devised a set of loose colour and magnitude cuts to reduce the initial DES sample to photometric i-band dropouts. We imposed the following cuts on the grizY Kron photometry (mag_auto) of the main DES source catalogue:

( mag_auto_i mag_auto_z > 0.8 & mag_auto_z mag_auto_y < 0.12 & mag_auto_r > 22.5 ) OR ( mag_auto_i mag_auto_z > 2.2 & magerr_auto_r > 0.36 ) . $$ \begin{aligned}&(\mathtt {mag\_auto\_i} -\mathtt {mag\_auto\_z} > 0.8 \\& \& \,\mathtt {mag\_auto\_z} -\mathtt {mag\_auto\_y} < 0.12 \\& \& \,\mathtt {mag\_auto\_r} > 22.5) \\&\mathtt {OR} \\&(\mathtt {mag\_auto\_i} -\mathtt {mag\_auto\_z} > 2.2 \\& \& \mathtt {magerr\_auto\_r} > 0.36). \end{aligned} $$

In addition, we always imposed signal-to-noise thresholds on the z-band photometry (S/N > 10) and on the Y-band photometry (S/N > 5). We further required no detection in the g band (S/N < 3). These colour cuts are summarized in Fig. 1. The colour cuts are more relaxed than the ones used in previous z ∼ 6 quasar searches based on colour criteria (e.g., Fan et al. 2001; Mortlock et al. 2011; Reed et al. 2015, 2019; Bañados et al. 2023; Matsuoka et al. 2022), as the aim here is merely to pre-select dropout objects. We queried a sub-sample of 511 654 candidates.

thumbnail Fig. 1.

Colour-selection diagram. The hashed and grey regions correspond to the colour space from which we pre-selected high-redshift quasar candidates for the present pilot survey. Gaussian Kernel density estimates of DES detected M-dwarfs (Best et al. 2018) and L/T-dwarfs (Liu et al. 2002; Hawley et al. 2002; Cruz et al. 2003, 2007; Chiu et al. 2008; Reid et al. 2008; Schmidt et al. 2010; Burningham et al. 2010, 2013; Kirkpatrick et al. 2010, 2011; Wahhaj et al. 2011; Day-Jones et al. 2013; Marocco et al. 2013, 2014, 2015; Best et al. 2015; Cardoso et al. 2015; Tinney et al. 2018) are displayed as shaded green and shaded blue areas, respectively. DES detected, spectroscopically confirmed z > 5.6 quasars compiled in Fan et al. (2023) and additional sources from Yang et al. (2024) are shown as salmon red circles. 1D colour distributions of the samples from the literature are shown in the marginal histograms. Magenta stars indicate new quasars discovered in this work. Four of these discoveries are located at the edge of the grey selection box, close to the contamination region of brown dwarfs.

2.1.2. SED template fitting

As the main candidate filtering step, we opted for spectral energy distribution (SED) template fitting and photometric redshift computation. Again, the DES magnitudes in the bands g, r, i, z and Y are derived using the Kron aperture radius (mag_auto). In addition to the DES photometry, we obtained magnitudes in the NIR and MIR wavebands by astrometrically cross-matching the main DES catalogue with the VHS and CatWISE2020 source catalogues within < 1″. By construction, we only considered sources simultaneously detected in DES, VHS and CatWISE2020, reducing our sample to 220 221 sources. Combining the relatively deep DES optical photometry with these shallower IR surveys biases our selection towards very red, potentially dusty sources. Using the photometric code Le PHARE (v2.2, Arnouts & Ilbert 2011), we performed a χ2-based template fitting analysis on the collected photometry of the colour pre-selected sample (see Section 2.1.1). All magnitudes are directly read from the catalogues. For the VHS bands Y, J, H and Ks, we used the default point source aperture corrected magnitudes (aperture 3, 2″ in diameter). Finally, for CatWISE2020 MIR magnitude measurements, we selected profile-fit photometry corrected for proper motion in the bands W1 and W2 (w1mpropm and w2mpropm).

Le PHARE was supplemented with a custom template library for AGN and galaxies already used for the eROSITA Final Equatorial-Depth Survey by Salvato et al. (2022). The original AGN templates used in this work are from Polletta et al. (2007), Salvato et al. (2009), Brown et al. (2019). We also used mixed AGN-galaxy templates from Salvato et al. (2009) and Ananna et al. (2017), based on inactive galaxy templates from Noll et al. (2004). Due to their steep red spectra, UCDs of the spectral classes M, L and T are known to be the most prominent contaminants in photometric searches of neutral-IGM absorbed, high-redshift quasars, as both classes of objects appear as dropouts in optical bands. The spectral slope in the IR bands can help disentangle the typical quasar power-law accretion disk emission from the Rayleigh-Jeans tails of brown dwarf spectra (for a review and quasar searches and contamination issues see Fan et al. 2023). In addition to stellar templates from Pickles (1998) and Bohlin et al. (1995), we accounted for the spectral diversity of these Galactic colour-contaminants by including the stellar templates provided by Chabrier et al. (2000), Knapp et al. (2004), Golimowski et al. (2004), Chiu et al. (2006), Bayo et al. (2011).

We tested our SED fitting analysis on all spectroscopically confirmed quasars at z > 5.6 listed by Fan et al. (2023), with detections within < 1″ in DES, VHS and CatWISE2020. This sample contains 35 quasars with spectroscopic redshifts 5.61 < z < 6.89 (initial discoveries by Mortlock et al. 2009; Bañados et al. 2016, 2023; Reed et al. 2017, 2019; Yang et al. 2017, 2021; Chehade et al. 2018; Decarli et al. 2018; Shen et al. 2019; Eilers et al. 2020; Venemans et al. 2020; Matsuoka et al. 2022). The Le PHARE photo-z calculating sub-routine zphota returns the χ2 values for the best-fit AGN and stellar templates. Out of the 35 currently confirmed DES, VHS and CatWISE2020 z > 5.6 quasars, 34 are assigned a photometric redshift that follows the outlier rejection criterion ∣zspec − zphot ∣ /(1 + zspec) < 0.15 within its 1σ uncertainties. (e.g., Hildebrandt et al. 2010; Salvato et al. 2022) of which 32 have been fitted with AGN models with reduced χbest2 smaller than the reduced χstar2 values of the stellar models, that is, the AGN model is favoured. The comparison of our photometric estimates of the redshift of these quasars and their actual spectroscopic redshifts is presented in Fig. 2.

thumbnail Fig. 2.

Le PHARE test run on previously spectroscopically confirmed z > 5.6 quasars. Upper panel: The estimated photometric redshift zphot is compared to the true spectroscopic value zspec. The grey shaded area indicates the region within which Δz < 0.15 × (1 + zspec) where Δz = ∣zspec − zphot∣. Lower panel: The ratio of χ2 of the stellar and the AGN or galaxy template are shown in the lower panel. In only two cases, the stellar solution is preferred.

We also ran the SED fitting routine on the DES, VHS and CatWISE2020 photometry of known UCDs (Liu et al. 2002; Hawley et al. 2002; Cruz et al. 2003, 2007; Chiu et al. 2008; Reid et al. 2008; Schmidt et al. 2010; Burningham et al. 2010; Kirkpatrick et al. 2010, 2011; Wahhaj et al. 2011; Burningham et al. 2013; Day-Jones et al. 2013; Marocco et al. 2013, 2014, 2015; Best et al. 2015; Cardoso et al. 2015; Tinney et al. 2018). Out of 51 late-type M, L and T dwarfs, 47 are correctly classified as stellar objects by Le PHARE.

Computing photometric redshifts on the colour pre-selected sample of 220 221 sources with matches in VHS and CatWISE2020, we obtained a list of 2605 zphot > 5.6 candidates. Their colour-redshift distribution is shown in Fig. 3.

thumbnail Fig. 3.

Optical colour-plane of selected of 2605 zphot > 5.6 quasar candidates from our colour-cuts and SED fitting. The candidates (circles) are colour-coded according to their photometric redshift. We show redshift-colour tracks for Le PHARE templates: a typical quasar (pl_QSO_DR2_029_t0, an SDSS-DR5 composite modified by Salvato et al. 2009) and a Seyfert II template (Sey2_template_norm, Polletta et al. 2007), with increasing extinction values (0.0, 0.1, 0.2) displayed by increasing transparency. Squares indicate selected zphot > 5.6 candidates that are already confirmed in literature (references in Fan et al. 2023). The five quasars discovered in this work are shown as magenta stars.

2.1.3. Sample cleaning

We matched our quasar candidate sample of 2605 sources to compilations of known high-redshift quasars and UCDs. We had blindly selected 27 previously confirmed quasars at 5.7 < z < 6.8 from Mortlock et al. (2009), Jiang et al. (2008), Bañados et al. (2016, 2023), Reed et al. (2017, 2019), Decarli et al. (2018), Chehade et al. (2018), Eilers et al. (2020), Venemans et al. (2020), Matsuoka et al. (2022), 77 % of all currently confirmed quasars in the joint DES, VHS and CatWISE2020 surveys. Additionally, our sample contained 7 known M, L and T dwarfs from Liu et al. (2002), Tinney et al. (2018) and Best et al. (2018). These sources are discarded. We also removed one known low-redshift galaxy from Colless et al. (2001) from our sample.

We further extracted forced photometry on optical DES images of the quasar candidates to identify artefacts and problematic blue-band PSF-matched photometry (e.g., Bañados et al. 2016, 2023). Circular apertures of radius 3″ were centred on the optical coordinates of the candidates. We discarded candidates for which we measured a significant aperture flux in g-band, with errors on g-band aperture magnitudes ap_magerr_g > 0.36 when a positive g-band flux is measured. Given the size of the source extraction region compared to the size of the DES point-spread-function (PSF < 1.11″, Abbott et al. 2021), this step effectively removes photometric dropouts that lie close to (or in close pairs with) the g-band detected stars and lower-redshift galaxies. Example images of discarded candidates are presented in Appendix A. Following this step and a careful visual inspection of each candidate, we conserved a sample of 2098 quasar candidates.

2.2. eROSITA down-selection

We searched for X-ray luminous quasars in eRASS (Merloni et al. 2024), and we therefore only accounted for quasar candidates contained in the intersection region of the optical/IR and eRASS footprints (1604 candidates). The combined area of our survey is measured with a coarse multi-order coverage using the python package pymoc. It spans 3507 deg2. The survey footprint and the on-sky distribution of quasars discovered in this work are shown in Fig. 4.

thumbnail Fig. 4.

The intersection of the DES, VHS, CatWISE2020 and eRASS footprints (blue shaded area) covers a sky surface of ∼3500 deg2. The underlying eRASS1 source population is shown in grey. The quasars discovered in this work are marked as pink stars.

We verified whether our candidates were X-ray detected in any of the eRASS scans (eRASS1-4) or the deepest, complete cumulative survey eRASS:4. The eROSITA data analyzed in this work were processed with the standard processing pipeline (Brunner et al. 2022, processing version 020). We performed forced photometry on the eRASS1-4 and eRASS:4 maps at the coordinates of the quasar candidates following a method presented by Arcodia et al. (2024). In each eRASS scan and the deepest available cumulative image (up to eRASS:4), we extracted counts within a circular source aperture of radius 30″, corresponding the half-energy width (HEW) of the eROSITA PSF in the most sensitive energy range (0.2−2.3 keV, Dennerl et al. 2012; Predehl et al. 2021; Brunner et al. 2022; Merloni et al. 2024). We note, however, that the HEW of the PSF increases to 34.4″ in the 2.3−5.0 keV range. Additionally, background photons are extracted in an annulus of standard radii 4 times and 12 times the source region radius. Significantly detected sources in the background region are masked. The binomial no-source probability Pb (e.g., Weisskopf et al. 2007; Vito et al. 2019) is given by:

P b ( i s ) = i = s s + b ( s + b ) ! i ! ( i s b ) ! ( 1 1 + r ) i ( 2 + r 1 + r ) s + b i . $$ \begin{aligned} P_b(i\ge s) = \sum _{i=s}^{s+b} \frac{(s+b)!}{i!(i-s-b)!} \left(\frac{1}{1+r} \right)^i\left(\frac{2+r}{1+r}\right)^{s+b-i}. \end{aligned} $$(1)

Here s and b are counts within the source and background region, while r is given by the ratio of areas of these two extraction regions. We set the threshold of Pb to be less than 0.7% in the wavelength range 0.2−2.3 keV, where eROSITA is most sensitive. According to the eRASS1 simulations of Seppi et al. (2022), this corresponds to a maximal fraction of ∼25% of spurious X-ray detections of background fluctuations. We note that this threshold is looser than the one applied by Arcodia et al. (2024, Pb < 0.0003), which is expected to correspond to a spurious fraction of 1%. However, assuming typical X-ray to optical flux ratios for Type 1 AGN (e.g., Steffen et al. 2006) z > 5.6 quasars are expected to lie right at the sensitivity limit of the survey (Wolf et al. 2021) and we therefore applied a more generous threshold. The fraction of spurious background fluctuation detections in eRASS, fbkg, predicted by Seppi et al. (2022) is shown as a function of Pb in Fig. 5.

thumbnail Fig. 5.

Chance alignment. Upper panel: The cumulative distribution of expected random (unrelated) X-ray detections from forced photometry using a circular aperture of radius 30″ around 1604 optical quasar candidate positions. The black curve is the estimate obtained from Eq. (2). for a single eRASS survey. The black dashed curve shows an empirical estimate of the same number based on a random sample of positions in the surveyed field. The brown curve shows the number of expected contaminants obtained from a random sample in eRASS:4. Blue vertical lines indicate the Pb values of the best detection of quasars discovered in this work in any eRASS and in eRASS:4 (if available). Lower panel: Fraction of background sources fbkg expected to be detected in eRASS1 from a simulation by Seppi et al. (2022).

Due to the large extraction radius of the rigid source region aperture, a caveat of our detection method is the non-negligible possibility of chance alignment, that is, the astrometric overlap of the extraction region centred on the optical coordinates of a quasar candidate and an unrelated, neighbouring eRASS source. The issue of such spurious matches in cross-survey catalogue matching has led to the development of probabilistic approaches (e.g., Sutherland & Saunders 1992; Brusa et al. 2007; Budavári & Szalay 2008; Salvato et al. 2018, 2022). We first estimated the probability of a fully spurious association between a quasar candidate and a single-eRASS (or an eRASS:4) X-ray source based on the eROSITA catalogues. The eRASS1 source catalogue contains 1 277 477 sources. Here, we have also included sources from the supplementary catalogue, with a detection likelihood threshold L = −ln(P)≥5 in the 0.2−2.3 keV band (Merloni et al. 2024). Here, P is the probability of a source arising from a random background fluctuation (Brunner et al. 2022). It is directly comparable to the binomial no-source probability Pb of Eq. (1). We count 302 161 eRASS1 sources in the common DES, VHS, CatWISE2020 and eRASS1 footprint. Following Pineau et al. (2017), assuming that the quasar candidates and the eRASS1 sources are completely unrelated samples and using circular positional uncertainties, the total number of fully random associations is given by:

n spurious = n eRASS ( P b ) n opt Ω eRASS ( P b ) Ω opt Ω , $$ \begin{aligned} n_{\rm spurious} = n_{\rm eRASS}(P_b)n_{\rm opt}\frac{\Omega _{\rm eRASS}(P_b)\Omega \mathrm{opt}}{\Omega }, \end{aligned} $$(2)

where ΩeRASS(Pb) = πσeRASS(Pb)2 and Ωopt = πσopt2. σeRASS and σopt are, respectively, the fixed source extraction radius in the eRASS map and the quasar candidates. Ω is the area of the overlap region of the eRASS and DES DR2 footprints. neRASS and nopt are the number of eRASS1 sources and candidates in this overlap region. The number of eRASS1 sources considered in the surveyed footprint, as well as typical positional uncertainties, are functions of the detection likelihood (DET_LIKE in the eRASS catalogues), or equivalently the binomial no-source probability Pb. We therefore binned the full eRASS1 catalogue in DET_LIKE and computed median uncertainties in each bin from the RADEC_ERR = σ RA 2 + σ Dec 2 $ \sqrt{\sigma_{\texttt{RA}}^2 + \sigma_{\texttt{Dec}}^2} $ values. As RADEC_ERR is the quadratic sum of the 1-dimensional errors in right-ascension and declination, and we assumed circular errors for our calculation (σRA = σDec), we considered error radii σeRASS = RADEC_ERR/ 2 $ \sqrt{2} $ for the eRASS1 sources (Pineau et al. 2017). nspurious can then be calculated for each DET_LIKE bin. The resulting curve of expected spurious associations is shown in Fig. 5. At Pb ∼ 10−3, less than ten candidates out of 1604 are expected to be randomly overlapping within 30″ with an unrelated eRASS1 source.

In addition, we obtained an empirical estimate of chance alignments by generating 1604 random positions in the surveyed region of the sky for eRASS1 and eRASS:4. We extracted source and background counts following the procedure detailed earlier in this section and computed Pb for each of these random positions. We obtained a distribution displayed in Fig. 5.

The final eRASS z > 5.6 quasar candidate sample contains a total of 75 sources in 3507 deg2 of the sky. Table 1 summarizes the number of detections in each eRASS. To each quasar candidate, we associated a probability of being randomly associated within 30″ to a background fluctuation as Pspurious = nspurious(Pb)/nT, where nT is the total number of quasar candidates. This value was used to prioritize targets for spectroscopic follow-up (lower values were favoured).

Table 1.

Number of detected quasar candidates in each eRASS.

thumbnail Fig. 6.

Spectra of newly discovered eRASS quasars (left panel) and eRASS:4 (0.2–8 keV) postage stamps (right panel). The cutouts are 2′×2′ in size and are centred on the optical position of the quasars. The spectrum of J202040−621509 shows prominent broad absorption features. The continuum appears strongly absorbed over the entire spectral range bluewards of Lyman α. It appears in the eROSITA images as two distinct sources. This morphology is further investigated in Section 4.3.1 with a Chandra follow-up pointing. There, we argue that the eRASS source is related to the quasar (Fig. 11). The morphology J020916−562650 is more compact.

3. Spectroscopic candidate confirmation

From our 75 candidates, we obtained optical spectroscopy for two objects so far. They were selected based on observability from the Las Campanas Observatory. In addition, as this pipeline was developed and calibrated, seven more targets were submitted for spectroscopic follow-up. All of them were selected as significant X-ray sources in eRASS data processed with an earlier version of the eROSITA reduction pipeline. In the latest processing version (020), they are not detected anymore according to the criterion defined in Section 2.2. This is likely due to random fluctuations and to the fact that these objects lie at the sensitivity limit of the survey.

The targets were observed with the LDSS3-C spectrograph mounted on the Magellan-Clay telescope at the Las Campanas Observatory (Chile). The observations took place on July 22, 2021, as well as on 2021 October 12 and 26. Long-slit spectra with slit-width 1″ were obtained using the VPH-Red grism, which spans the wavelength range 6000−10 500 Å at a dispersion of about 1.16 Å/px. Total exposures varied between 1800 s and 3600 s. A summary of these observations is given in Table 2. Overscan subtractions, flat-field corrections, wavelength, and flux calibrations were performed with IRAF.

Table 2.

Summary of the 5 newly discovered quasars.

Both of the eRASS-detected quasars candidates are confirmed as high-redshift quasars3. Their LDSS3-C discovery spectra and eRASS:4 cutouts are presented in Fig. 6. In addition, 3 of the remaining 7 targets observed with LDSS3-C in the development phase of the pipeline are also confirmed as high-redshift quasars. The spectra of these objects are presented separately in Fig. 7. The quasars are identified by the clear, abrupt, and near-total absorption of their blue continua, attributable to the Lyman α absorption by the increasing density of reionization-era neutral hydrogen clouds. As the unabsorbed continuum and most emission lines of the quasars are shifted out of the sampled spectral range, redshifts can only be measured from the Lyman absorption break. The measurement precision is affected by the low signal-to-noise of some of the spectra. Absorption effects redwards of Lyman α such as IGM damping wings (Rybicki & Lightman 1979; Miralda-Escudé 1998) can also severely affect Lyman α redshift measurements. We fitted a step function with the python package lmfit (StepModel() + ConstantModel()) to the quasar spectra. We estimated the redshift and its error from the break centre by bootstrapping from the flux errors. We obtained values between z ∼ 5.6 − 6.1. The remaining 4 targets all have spectra corresponding to brown dwarf templates (the coordinates of these sources are presented in Appendix B). We further extrapolated the restframe absolute magnitude at 1450 Å, M1450, from the reddened y-band magnitudes in the DES DR2 catalogues using the spectral index aν = −0.44 (Vanden Berk et al. 2001). We note that because of the simultaneous DES, VHS, and CatWISE2020 selection and the respective depth of these surveys, the unabsorbed continuum slopes of the selected quasars are biased towards the red values. The choice of the extrapolation from the y-band mitigates the uncertainty due to the expected scatter in the spectral index, as its central wavelength is relatively close to 1450 Å ×(1 + z). M1450 measurements are listed in Table 2.

thumbnail Fig. 7.

Spectra of three newly discovered quasars selected during the calibration of the pipeline. These objects are not robustly detected in eROSITA. The red dashed line indicates the location of the Lyman α break.

4. Two X-ray and radio luminous quasars

In this section, we focus on the properties of the new eROSITA quasars J202040−621509 and J020916−562650. Both are radio-bright but are intrinsically very different sources. We identify them respectively as a high-redshift jetted broad absorption line quasar (BAL QSO) and a high-redshift blazar.

4.1. eROSITA photometry

We derived the X-ray properties of the two newly discovered and significantly eROSITA-detected quasars from the source and background spectra measured in the eRASS maps at their optical coordinates. The source and background regions are defined in Section 2.2 and the extraction bands are chosen to be broad (0.5−7 keV). Fluxes were derived in single eRASS scans, as well as in the cumulative eRASS:4 images. Extracted counts are converted to fluxes using a redshifted power-law model with fixed photon index (XSPEC model zpowerlw4, Γ = 1.9) and Galactic absorbing column density NH from HI4PI Collaboration (2016) with the interstellar absorption model tbabs (Wilms et al. 2000). We opted for a power-law slope Γ = 1.9 typically observed in lower-redshift AGN populations (e.g. Liu et al. 2022). We note, however, that steeper value of the photon index, potentially related to high accretion rates, have been reported in luminous X-ray spectra of high-redshift quasars (Γ ∼ 2.2 − 2.4, Vito et al. 2019; Wang et al. 2021b; Li et al. 2021; Zappacosta et al. 2023). The model was scaled using the Bayesian X-ray spectral analysis code BXA (Buchner et al. 2014). The spectra are rebinned with the ftgrouppha function of the FTOOLS package5, ensuring at least one bin in each count. The data are analysed using the Cash statistic (w-stat, Cash 1979; Wachter et al. 1979). We obtained median values, as well as the 15.9th and 84.1th percentiles6 of the posterior flux chains. For the single epoch photometry, we report the flux in the eRASS for which Pb reaches its minimum. We only consider epochs with Pb < 0.007. Quasar J202040−621509 is only detected in eRASS1 (Pb = 4.4 × 10−4), while J020916−56265 is detected in eRASS2, eRASS3, and eRASS4 (with a minimum of Pb = 2.2 × 10−7). The cumulative photometry is computed from the stacked survey eRASS:4 (the deepest, homogeneous, and complete cumulative eRASS survey). We summarize the X-ray properties of the discovered quasars in Table 3. We derived restframe hard X-ray luminosities L2 − 10 keV assuming the above AGN spectral model. Extrapolating from the dereddened y-band DES magnitudes assuming the spectral slope of α = −0.44 (Vanden Berk et al. 2001), we computed an estimate of the restframe monochromatic luminosity at 2500 Å, L2500. We computed the two-point spectral index αOX connecting the optical and X-ray part of the quasar SED as:

α OX = 0.384 × log ( L 2 k e V / L 2500 ) , $$ \begin{aligned} \alpha _{\mathrm{OX} } = 0.384 \times \mathrm{log} \, (L_{2\,\mathrm {keV}}/L_{2500}), \end{aligned} $$(3)

Table 3.

eROSITA X-ray properties from aperture photometry.

where L2500 is the luminosity at 2 keV.

In Fig. 8, left panel, we present the L 2 10 keV $ L_{2{-}10\,\rm keV} $–redshift plane X-ray detected quasars in a redshift range 5.5 < z < 7.5 and a luminosity range 44 < log L2 − 10 keV/(erg/s) < 46.5 (Nanni et al. 2017; Vito et al. 2019; Pons et al. 2020; Medvedev et al. 2020; Khorunzhev et al. 2021; Wang et al. 2021b; Wolf et al. 2021, 2023; Zappacosta et al. 2023). The average eROSITA sensitivity in units of fractional sensitive area is shown as a colour map (values for a typical equatorial eROSITA field). It was computed by forward-modelling the sensitivity of the instrument to a fiducial AGN, described here by an absorbed power-law with Γ = 2 and Galactic absorption, redshifted to z = 6 (see Wolf et al. 2021, 2023). The quasars J202040−621509 and J020916−562650 are among the most X-ray luminous quasars known to date. The spectral index αOX is known to anti-correlate steeply with rest-frame UV luminosity L2500 for Type 1 AGN (e.g., Steffen et al. 2006; Just et al. 2007; Lusso & Risaliti 2016). We show the relation measured by Nanni et al. (2017) for a sample of z > 5.5 quasars. In the right panel, we compare the αOXL2500 of the newly discovered quasars with the sample X-ray detected ones from the literature. The αOX values of J202040−621509, J020916−562650, and J230341−542730 significantly deviate from the measured relation. Their location indicates that they are X-ray over-luminous.

thumbnail Fig. 8.

X-ray emission of the new eRASS quasars. Left panel: X-ray luminosity redshift distribution of X-ray luminous quasars from literature (Nanni et al. 2017; Vito et al. 2019; Pons et al. 2020; Wang et al. 2021b; Medvedev et al. 2020; Belladitta et al. 2020; Wolf et al. 2021, 2023; Zappacosta et al. 2023) spanning a redshift range 5.5 < z < 7.5 and a luminosity range 44 < log L2 − 10 keV/(erg/s) < 46.5. The two newly discovered quasars (displayed as pink stars). The colour gradient shows the sensitivity to a fiducial absorbed power-law with Galactic absorption of the final cumulative eRASS in the equatorial region (in fractional of sensitive area). The quasars discovered in this work lie at the luminous end of the quasar population in the early universe. We show the luminosity derived from the Chandra follow-up observation of J202040−621509. Between the eRASS1 and the Chandra observation, its luminosity has decreased by an order of magnitude. Right panel: αOXL2500 distribution of the same sample of sources. The 1σ confidence interval of relation derived by Nanni et al. (2017) is shown by red dotted line. The hexagonal pattern shows a lower redshift AGN sample by Lusso & Risaliti (2016). The three quasars X-ray detected and newly discovered are over-luminous in the X-ray wavebands. Following its dimming observed in the recent Chandra observation, J202040–621509 is perfectly consistent αOXL2500 relation.

4.2. eRASS J020916–562650: A blazar at z = 5.6

4.2.1. Radio properties

J020916−562650 is clearly detected at 0.887 and 1.37 GHz in the Rapid ASKAP Continuum Survey (RACS, McConnell et al. 2020; Hale et al. 2021), at 0.887 in the ASKAP Variables and Slow Transients survey (VAST, Murphy et al. 2013) and at 0.943 GHz in the Evolutionary Map of the Universe survey (EMU, Norris et al. 2021). To quantify the flux densities at these frequencies, we performed a single Gaussian fit using the task IMFIT of the Common Astronomy Software Applications package (CASA, McMullin et al. 2007) on the images retrieved from the CASDA7 website. They are reported in Table C.2. By assuming < a single power-law for the continuum radio emission (Sν ∝ να) we estimated a radio spectral index (αr, see Fig. C.1) of 0.14 ± 0.10, which classifies it as a Flat Spectrum Radio Quasar (FSRQ, Urry & Padovani 1995). The two flux densities estimated from the VAST images show a hint of variability on a time scale of ∼2 months (in the observed frame, and thus, ∼9 days in the rest frame). Then, we computed the values of radio-loudness (R), which quantify how powerful is the radio (i.e., synchrotron) emission with respect to the optical/UV (i.e., accretion disk) one. We used the definition of Kellermann et al. (1989): R = S 5 GHz / S 4400 Å $ R = {S\!}_{\mathrm{5\,GHz}}/{S\!}_{\mathrm{4400\,\AA}} $. We derived the flux density at 5 GHz rest frame by using the available observed flux densities (0.87, 0.94 and 1.3 GHz) and the estimated radio spectral index. The value of the flux density at 4400 Å rest frame was derived from the z-band magnitude assuming the optical spectral index of αν = 0.44 (Vanden Berk et al. 2001). We obtained a mean R of 750 ± 170 (Log(R) = 2.87). The high value of radio loudness, the flat radio spectral index and tentative variability indicate that this source can be classified as a blazar, a radio-loud AGN with the relativistic jet pointed towards the Earth (e.g., Urry & Padovani 1995).

4.2.2. X-ray spectral analysis

To further investigate the potential blazar nature of the source, we attempted a fit of the eRASS:4 spectrum of J020916−562650 by modelling it once again by a power-law with Galactic absorption (NH = 2.18 × 1020 cm−2, HI4PI Collaboration 2016): zpowerlw × tbabs. The photon index Γ is left free to vary in the fit. We stress that the low-counting statistics (∼25 net counts in the source region) significantly limit our ability to obtain a well-constrained X-ray spectral shape for this source. The analysis is performed with BXA (Buchner et al. 2014). We define wide and uninformative priors for both the normalization and the slope of the power-law: Γ = 0.5 − 5 (uniform prior) and N = 1−6 − 1 (Jeffrey’s prior). The eRASS:4 source and background spectra are extracted in regions as defined in Section 2.2. We have used the srctool function of the software package eSASS (user version eSASSusers_240410) to obtain these spectra. We used the redshift of the source measured in Section 3. The posterior distributions of the photon index and the hard X-ray luminosity L2 − 10 keV we obtained are presented in Fig. 9. Due to the poor counting statistics, the photon index is not well constrained: Γ = 1 . 46 0.42 + 0.44 $ \Gamma = 1.46^{+0.44}_{-0.42} $. Flatter values than typical AGN are, however, preferred (e.g., Nandra & Pounds 1994; Reeves et al. 1997; Ishibashi & Courvoisier 2010). This flatter photon index value is consistent with values typically observed in the restframe hard X-ray spectra of lower redshift blazars (e.g., Langejahn et al. 2020). Deeper X-ray data will be required to confirm the flatness of the X-ray spectrum.

thumbnail Fig. 9.

Posteriors on the photon index Γ and the X-ray luminosity L2 − 10 keV from the X-ray spectral analysis performed on the eRASS:4 spectrum of J020916−562650.

The luminosity we obtained with this photon index is log L 2 10 keV / [ erg s 1 ] = 45 . 89 0.12 + 0.11 $ L_\mathrm{2{-}10\,keV}/[\mathrm{erg\,s^{-1}}] = 45.89^{+0.11}_{-0.12} $. From Eq. (3) and this luminosity estimate, we obtained a spectral index α OX = 1 . 15 0.01 + 0.01 $ \alpha_{\mathrm{OX}}=-1.15^{+0.01}_{-0.01} $. This relatively flat value, given the quasar’s optical monochromatic luminosity, indicates that the quasar is X-ray over-luminous with respect to typical Type 1 AGN. Such a flat value is typical for blazars (e.g., Donato et al. 2001).

4.2.3. A blazar SED

In the previous section, we were able to constrain Γ and αOX values for J020916−562650 that are typical for blazars. We further have shown that the radio morphology, spectral shape, and variability of the associated radio source are consistent with the hypothesis that the source is jetted and that its jet aligns with the line of sight. Fig. 10 shows the broadband rest frame SED of the source, with data from radio to the X-rays. The observed data are compared with the analytic approximation of a typical FSRQ SED as developed by Ghisellini et al. (2017), combined with the thermal accretion emission from a standard geometrically thin, optically thick α-disk (Shakura & Sunyaev 1973), with relative X-ray Corona and dusty IR torus. Such thermal emission clearly dominates the rest-frame NIR-optical-UV energy range in this case. Thus the DES, VHS and CATWISE2020 photometric points have been fitted with a Shakura & Sunyaev (1973) model, following a method already introduced by Sbarrato et al. (2012), Calderone et al. (2013), Ghisellini et al. (2015). Belladitta et al. (2019, 2022). Following this approach, we can conclude that the central engine is likely fast accreting. Its emission is consistent with a log MBH/M = 8.60 black hole, accreting at the Eddington limit. The radio and X-ray data points of J020916−562650 are instead clearly dominated by the emission from a relativistic jet, aligned closely to our line of sight. Even with a sparse sampling of the broad-band SED, the source shows the typical double hump SED profile of blazars: the hard and intense X-ray luminosity can only be justified by an inverse Compton emission, being strongly overluminous with respect to the expected Corona emission. The significant radio luminosity further confirms the blazar interpretation, tracing an intense synchrotron emission.

thumbnail Fig. 10.

The broadband rest frame SED of J020916−562650. A range of typical FSRQ SEDs (Lγ ∼ 1046 − 48 erg/s, grey shaded region) and a more tailored blazar analytical model (solid green line) are overplotted to the data. The radio and X-ray emission of the source are consistent with the typical blazar double hump SED. Specifically, the X-ray emission significantly overestimates the expected X-ray coronal emission, extrapolated by the source accretion emission that dominates the optical/UV restframe photometry (yellow dots). The thermal emission is modelled here with a standard accretion disk emission, IR torus, and X-ray Corona powered by a black hole with the mass of log MBH/M = 8.60 (blue solid and dashed lines, the latter showing the unabsorbed accretion disk SED).

4.2.4. AGN demographics from a single blazar discovery

Volonteri et al. (2011) and Ghisellini et al. (2010) developed the geometric argument that for each blazar, that is, a line-of-sight-aligned jet (viewing angle below the critical angle < 1/ΓLorentz, where ΓLorentz the bulk Lorentz factor), there must exist 2 × ΓLorentz2 parent jets in quasars pointed in any other direction. Following the above argument, Belladitta et al. (2020) derived a lower limit on the space density of AGN similar in the same redshift-magnitude bin than PSO J030947.49+271757.31, the most distant blazar known to date. Assuming a typical value for the bulk Lorentz factor (ΓLorentz = 10, e.g., Ghisellini et al. 2014) and considering a redshift interval z = 5.6 − 6.5 over the surveyed area of 3507 deg2, we obtained a number density of quasars at least as bright as J020916−562650 (M1450 < 25.76) of 6 . 89 5.70 15.86 $ 6.89^{15.86}_{-5.70} $ Gpc−3. The Poisson confidence limits are estimated according to Gehrels (1986). This lower limit on the space density of the parent jetted quasar population is consistent with Belladitta et al. (2019) at similar absolute magnitudes ( 1 . 10 0.91 + 2.53 $ 1.10^{+2.53}_{-0.91} $ Gpc−3 at M1450 ≤ 25.1 and within z = 5.5 − 6.5). Belladitta et al. (2019) compared their result to the integrated quasar luminosity function by Matsuoka et al. (2018) ( 3 . 64 0.03 + 0.09 $ 3.64^{+0.09}_{-0.03} $ Gpc−3 at M1450 < −25.1), concluding that 30% of AGN at z = 5.5 − 6.5 are radio loud.

Integrating the quasar luminosity function of Schindler et al. (2023) over the same redshift range down to M1450 < 25.7 yields a space density ∼1.32 Gpc−3. The density of luminous jetted quasars derived from the discovery of J020916−562650 is consistent with the previous density measurements at these redshifts and luminosities by Jiang et al. (2016), Matsuoka et al. (2018), Schindler et al. (2023). We note, however, that the argumentation above requires that all of the quasars we accounted for are similarly radio-loud and jetted. This suggests that either most of the jets of the quasars discovered so far in this redshift range have not yet been detected in radio surveys or that a largely hidden population of potentially obscured quasars have eluded dedicated surveys, as the radio-loud fraction of quasars is not expected to exceed ∼30% at these redshifts (Bañados et al. 2015; Belladitta et al. 2019). A similar conclusion was reached by Caccianiga et al. (2024), who compared the number of radio detected quasars to the number of blazars. A possible source of error in our estimate is the value we have assumed for the bulk Lorentz factor or the ΓLorentz2 dependence on the number density of parent jets. Lister et al. (2019) argue that the simple geometric argument we have made use of fails to account for biases arising from flux-limited sampling. These authors show through Monte Carlo simulations that the evolution of the fraction of jetted quasars to blazars (brighter than 1.5 Jy at 15 GHz) as a function of Γ is shallower than the assumed relation ∝ΓLorentz2 evolution at ΓLorentz > 15. We note, however, that Lister et al. (2019) a derive steep relative evolution of blazar jets and their parent population at ΓLorentz < 15, resulting in Nparent ∼ 103 × N> 1.5 Jy. Due to the large surveyed area and the rarity of AGN, the Poisson counting uncertainties fully dominate cosmic variance (e.g., Trenti & Stiavelli 2008; Dai et al. 2015). Our results underline the necessity of systematic, multi-wavelength searches for high-redshift blazars to better constrain AGN demographics in the poorly sampled luminous regime.

4.3. eRASS J202040–621509: X-ray dimming in a radio-bright z = 5.7 BAL quasar?

4.3.1. Confirming the X-ray detection with Chandra ACIS-S

To increase counting statistics and constrain the X-ray spectral shape of J202040−621509, we obtained an additional Chandra ACIS-S pointing. J202040−621509 was observed with Chandra ACIS-S on September 23, 2023 as Cycle 23 GTO target (PI: Predehl, Observer: Wolf, obsid: 26040). The total exposure time was 29.58 ks. We analysed the data with the X-ray ray data analysis software ciao 4.15. We reduced the data with the latest pipeline and calibration database version. To detect sources in the image, we ran the Mexican-hat wavelet algorithm wavdetect on the obtained broad band images (0.5−7 keV). We used the default detection parameters (sighthresh = 10−6 and bkgsigthresh = 0.001). A source is significantly detected within 0.3″ of the quasar with 9.9 ± 3.2 net counts. Given the astrometry, the X-ray source is unequivocally associated with the quasar. We refer to this source as ch1 henceforth. Two more sources, which we label ch2 and ch3, are detected at a lower significance level within 30″ of the quasar. They are located to the north and northeast of the quasar and do not overlap with the eROSITA X-ray centroids to the south. The eRASS:4 and Chandra extraction images are overlayed in Fig. 11 (left panel). From the eRASS:4 image, we counted 7 photons in the source region, making the source position and morphology highly uncertain. An eRASS1 catalogue source, 1eRASS J202039.8−621525, is located within 17.71″ of the optical coordinates of the quasar (green triangle). To help us determine the source of the eROSITA signal, we performed a cross-match between the eRASS1 catalogue and the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Survey (Dey et al. 2019) with the probabilistic cross-matching algorithm NWAY (Salvato et al. 2018). Our method combines optical/IR photometric priors derived from a truth sample of X-ray sources from the Fourth XMM-Newton Serendipitous catalogue (4XMM, Webb et al. 2020) with the astrometric configuration (distances, counterpart candidate surface density and positional errors). The procedure is presented in detail by Salvato et al. (2022). In our setup, all DESI Legacy Survey sources within a radius of 1’ around an eRASS1 source are considered as potential candidates. 151 such DESI Legacy Survey candidates are found. We show a DES DR2 z-band image showing the respective coordinates of the centroid of eRASS1 J202040−621509, ch1, ch2, ch3 and the DESI Legacy Survey neighbouring candidates in Fig. 11 (central panel). NWAY outputs the matching statistics p_any and p_i, which are respectively the probability that any of the candidates is matched to the X-ray source and the relative probability among candidates. In the case of 1eRASS J202039.8−621525 (DETUID:em01_306153_020_ML00171_002_c010, DET_LIKE_0 = 9.95), we obtained p_any = 0.012. The best matching source is a DESI Legacy Survey source matched within 0.16″ with the quasar DES coordinates (p_i = 0.97). We note that closer DESI Legacy Survey counterpart candidates are down-weighted by the custom photometric priors we employed, as their colours do not correspond to that of typical X-ray sources. Preliminary calibrations using a random test sample show that a p_any = 0.012 roughly corresponds to a purity of ∼45% at a completeness of ∼100% in the retrieved counterpart samples. From this exercise, we cannot conclude with certainty that 1eRASS J202039.8−621525 is associated with the quasar, but we can exclude the other DESI Legacy Survey sources in the field as the source of the X-ray emission to a high degree of certainty (p_i < 0.02). Many of the DESI Legacy Survey sources appear to be duplicate our spurious. The offset Chandra source ch2 appears matched to a catalogue DESI Legacy Survey source, while ch3 has no obvious counterpart. Performing forced photometry on background-subtracted DES images within circular source apertures of radii equal to the semi-major axis of the Chandra error ellipses (respectively 2.5″ and 1.5″). We obtained a S/N < 3 in all DES bands for both sources, suggesting that these two Chandra sources are not matched in DES.

thumbnail Fig. 11.

Astrometric configuration of J202040−621509. All images are centred on the optical position of the new quasar J202040−621509 and have sizes 100″ × 100″. Left panel: Broadband (0.5−7 keV) Chandra image. Chandra sources detected with 30″ (black, dashed circle) of the quasar coordinates with wavdetect are marked in red. The centroid of 1eRASS J202039.8−621525 catalogue detection is marked by a green triangle. The contours were obtained from the smoothed eRASS:4 0.5−7 keV image. Central panel: DES DR2 z-band image. Circular coloured markers denote counterpart candidates for 1eRASS J202039.8−621525 from the DESI Legacy Survey catalogue. They are colour-coded according to their relative probability p_i of being the best counterpart. The DESI Legacy Survey source with the highest p_i is spatially co-incident with the quasar coordinates. Right panel: CatWISE 2020 W1 image. Among the Chandra and eROSITA X-ray detections, only the quasar is unambiguously detected.

We extracted fluxes for the three Chandra sources using the srcflux routine of the ciao software package. For all three sources, we assumed a fiducial power-law model zpowerlw with Galactic absorption and fixed photon index (Γ = 1.9). We report the broad-band photometry for these three detections in Table 4. ch1 is the most significant detection among the three sources near the quasar. We note that no counts are registered in the ultrasoft band (0.2−0.5 keV) and the soft band (0.5−1 keV) for ch1. For this source, we measured 5.85 ± 2.45 net counts in the medium band (1−2 keV) and 4 ± 2 counts in the hard band (2−7 keV). This hard photometry can be indicative of a significant level of obscuration, as no counts can be measured in the restframe 1.34−6.7 keV energy of the quasar SED.

Table 4.

Chandra ACIS-S photometry for three sources detected within 30″ of the quasar.

In summary, we have presented evidence that supports that 1eRASS J202039.8−621525 is associated with the quasar emission, as (A) the Chandra sources ch2 and ch3 are marginally detected and do not overlap with the eRASS contours. (B) These sources are also not significantly detected in DES and CatWISE2020. (C) Among 151 potential DESI Legacy Survey candidates within 1′ of 1eRASS J202039.8−621525, the candidate with the highest relative probability p_i of being associated to the X-ray source is matched at the sub-arcsecond level with quasar coordinates. (D) No other Chandra source is detected to the southwest of the quasar, close to eROSITA events. While we cannot exclude the possibility of a transient phenomenon in a foreground galaxy, we conclude that it is likely that the bright eRASS1 source is associated with the quasar.

4.3.2. Long-term X-ray variability

As presented in Table 3, J202040−621509 is only detected in eRASS1, but not in any subsequent scan.

The eROSITA image presented in Fig. 11 displays two distinct lobes, which extend to the south and southwest of the coordinates of the quasar. One of these lobes coincides astrometrically with the quasar. In the Chandra image, three sources were detected with the wavdetect algorithm. One of these sources is matched at the sub-arcsecond level to the quasar’s coordinates. The two other sources are detected at a lower significance level to the north and the northwest of the quasar. Their location does not match the eRASS contours, and they are associated with DES sources within 1″. Assuming that the eRASS1 photometry is solely associated with the quasar, the eRASS and Chandra photometry presented respectively in Sections 4.1 and 4.3.1 imply a decrease in the measured 0.5−7 keV flux of the source of about an order of magnitude from F 0.5 7 keV = 7 . 82 4.47 + 11.06 × 10 14 $ F_{0.5{-}7\,\mathrm{keV}} = 7.82^{+11.06}_{-4.47} \times 10^{-14} $ erg s−1 cm−2 to F 0.5 7 keV = 5 . 83 2.63 + 3.70 × 10 15 $ F_{0.5{-}7\,\mathrm{keV}} = 5.83^{+3.70}_{-2.63} \times 10^{-15} $ erg s−1 cm−2. We note that the flux decrease is only significant at the ∼2σ level due to the low number of net source counts in eRASS1. The two observations are separated by 3.5 years in the observed frame, corresponding to 6 months in the quasar restframe.

Fig. 12 shows the flux-decrease in the broad 0.5−7 keV band. The BXA flux posteriors for the eRASS1 and Chandra, their median and 15.9th and 84.1th percentiles are shown. We have further derived L2 − 10 keV and αOX from the this Chandra flux posteriors. They are respectively marked in the left and right panels of Fig. 8.

thumbnail Fig. 12.

Distribution of BXAF0.5 − 7 keV fluxes for the eRASS 1 (in blue) and Chandra observations (red) of J202040−621509. The vertical black lines show the locations of the 15.9-th and 84.1-th percentiles of the distributions. The green and pink regions span the 2.3-th and 97.7-th percentiles of the distributions. The lower panel shows the time separation in days between the two observations. We measured a decrease in the source flux between the two epochs significant at the ∼2σ level.

4.3.3. Fast outflows

The discovery spectrum of J202040−621509 shows multi-trough, broad absorption lines (BALs). We first identified the following clear BALs in the quasar spectrum: NV, SiIV and CIV, which all show similar double-trough shapes (e.g., Korista et al. 1993). In quasar spectra, BALs are imprinted by extreme outflows. Such outflows are observed in about 12% to 20% of the quasar population (e.g., Knigge et al. 2008; Scaringi et al. 2009) and their velocity widths range can reach up to ∼30 000 km s−1 (so-called ultra-fast outflows, Rogerson et al. 2016). BALs are defined to have velocity widths > 2000 km s−1 and reach maximal depths of at least 10% below the quasar continuum level. This definition was introduced by Weymann et al. (1991) to trace outflows while avoiding contamination by host and extra-galactic absorbing systems. As the NV trough is affected by the neutral IGM absorbtion shortwards of Lyman α, we focused our analysis on the outflow velocity structure of CIV and SiIV. A measure of the trough widths, also called balnicity index (BI) was defined by Weymann et al. (1991) as:

BI = 3000 25 000 [ 1 f ( v ) 0.9 ] C d v , $$ \begin{aligned} \mathrm{BI} = \int _ {3000}^{25\,000} \left[1 - \frac{f(v)}{0.9}\right] C \, dv, \end{aligned} $$(4)

where f(v) denotes the continuum-normalized flux, C is a constant equal to 1 for velocities > 2000 km s−1 with respect to the trough edge and 0 elsewhere. The measurement of the continuum is challenging due to the numerous BALs. We approximated the continuum by using the principal component decomposition of z ∼ 3 SDSS quasar continua presented by Pâris et al. (2011). The spectrum is strongly attenuated toward longer wavelengths, which is not accounted for by the PCA model. We fitted an additional power-law in the restframe spectral regions 1300−1350 Å and 1400−1600 Å to the spectrum normalized by the PCA-reprojected model. Finally, the spectrum of the source was normalised by the combined PCA and power-law models. We accounted for the flux uncertainties by bootstrapping from the spectrum RMS. From the normalized spectrum we measured BICIV = 3394 ± 197 km s−1 and BISiIV = 2935 ± 110 km s−1. This BAL strength confirms the classification of J202040−621509 as BAL quasar.

The large BALs are detached with respect to the quasar redshift; the troughs of all species are shifted shortwards of the expected corresponding quasar emission lines. We measured detached and terminal velocities from the start and the end of the absorption troughs. From the SiIV BAL we measured outflow velocities vSiIV = 0.002c − 0.020c and vCIV = 0.003c − 0.025c. We note that the optical spectrum and the X-ray data were not obtained simultaneously. The LDSS3 observations took place closer to the eRASS1 observation (∼1 year in the observed frame) than the Chandra ACIS-S pointing (∼2 years), closer to the high-luminosity state discussed in Section 4.3.2.

4.4. Radio properties

J202040−621509 is detected at radio frequencies, at 0.887 and 1.37 GHz in VAST and at 0.943 GHz in the EMU Pilot survey (Norris et al. 2021). We retrieved the corresponding radio images from CASDA and we performed a single Gaussian fit using the task IMFIT within the CASA software to quantify the flux densities (see Table C.1). The images are shown in Fig. C.2 (top two rows) in Appendix C. Several VAST images are available, corresponding to different observing dates. We extracted the flux density from each of them. By assuming a single power-law for the continuum radio emission, we obtained a radio spectral index (αr) of −0.07 ± 0.15 (see Fig. C.1) that makes J202040−621509 FSRQ, similar to J020916−562650. The VAST measurements are consistent with each other, indicating that the source does not appear to be variable on time scales of a few months (observed frame). Then we computed the values of radio-loudness (R) by following the same procedure as for J020916−562650. We obtained a mean R of 10 ± 4 which places the source at the boundary between objects classified as radio-quiet or radio-loud (Urry & Padovani 1995). This could indicate the presence of a relativistic jet (e.g., Sbarrato et al. 2021).

4.5. Obscuring fast accretion disc winds at high redshift and other variability scenarios

AGN shows substantial variability in flux and spectral shape at all frequencies (for a review, see Lira 2021). This variability arises from accretion physics (e.g., Nandra et al. 1997), broad line region kinematics (e.g., Cherepashchuk & Lyutyi 1973; Boksenberg et al. 1978; Blandford & McKee 1982; Peterson 1993) and jets (e.g., Albert et al. 2007). In X-rays, AGN has been observed to be variable on timescales of hours to days (e.g., Winkler & White 1975; Marshall et al. 1981; McHardy et al. 2004; Ponti et al. 2012). Little redshift evolution in the X-ray variabity of AGN, as traced by the variance-frequency diagram, has been observed from z = 0 to z = 3 (Paolillo et al. 2023). Interpreted as reverberation effect, that is, reprocessing of accretion disk photons in the electron corona (Haardt & Maraschi 1991), rapid X-ray variability can be used to study the structure of the innermost region of the AGN at a couple of gravitational radii away from the black hole (e.g., McHardy 2013). Coronal variability usually shows maximum amplitudes of a factor ∼2 (e.g., Gibson & Brandt 2012). However, more extreme variability, unrelated to coronal fluctuations, has been observed on longer timescales (e.g., Timlin et al. 2020). For instance, (LaMassa et al. 2015) reported the first discovery of a changing look quasar in which the AGN flux dropped by a factor of 6 over timescales of years while not showing absorption signatures. As presented in Section 4.3.3, the UV spectrum of J202040−621509, in contrast, is that of a typical BAL quasar. The X-ray variability in quasars containing BALs in their restframe UV spectra is thought to arise from absorption effects. More precisely, variations in the column density of intervening, obscuring clouds are expected to arise from large-scale outflows (e.g., Brandt et al. 2000; Green et al. 2001; Gallagher et al. 2002, 2006; Gibson et al. 2009; Saez et al. 2012). An often invoked physical model for the observed extreme outflow velocities is that of line-driven winds launched from the accretion disk, as high-accretion rates drive the UV radiation pressure of the disk up. High-density regions, located at the inner edge of this flow can shield the coronal X-ray emission through photo-electric absorption (NH > 1022 cm−2, Murray et al. 1995).

A possible scenario is that the observed X-ray dimming of J202040−621509 over one order of magnitude in X-ray flux may partly be explained by such obscuring winds. While too few counts have been collected in the Chandra observation to extract and fit an X-ray spectrum, we only collected counts at energies > 6.7 keV in the quasar rest frame. This hard photometry can, for example, be indicative of heavy X-ray absorption, affecting lower-energetic photons.

Extreme X-ray variability in high-redshift quasars has been reported previously (for example, at z = 5.41, Shemmer et al. 2005). Within the first gigayear, only two objects were known to display such a behaviour (Nanni et al. 2018; Vito et al. 2022). Nanni et al. (2018) reported a hardening of the spectrum and a flux decrease of factor 2.5 in the quasar SDSS J1030+0524 (z = 6.31) between two XMM-Newton observations spaced by 14 years. They proposed intrinsic variation of accretion rate as a potential driver of the observed variability. Vito et al. (2022) presented the surprising non-detection in a 100ks XMM-Newton observations of the quasar CFHQS J164121+375520 (z = 6.047), which had been detected 115 quasar-restframe days before with Chandra, implying a flux-dimming of factor > 7. The authors discuss intervening obscuring material on sub-pc scales, such as disk-shape variability in super-Eddington sources. Vito et al. (2022) also list a possible model by Giustini & Proga (2019), in which rapidly accreting black holes with ≥ 0.25, generate the favourable conditions for their optically-thick, thin disk emission to dominate the output of a cooled corona. These disks are expected to launch strong line-driven winds, imprinting BAL features in the spectrum of the quasars.

While the latter model yields a promising scenario for the phenomenology observed in J202040−621509, we cannot exclude any of the above-mentioned sources of variability as the counting statistics of our X-ray observations are poor and we are missing simultaneous UV photometry to simultaneously track the corona and disk emission. A deeper, simultaneous X-ray and UV monitoring campaign will be required to shed light on this seemingly extreme X-ray dimming event.

5. Conclusions

In this work, we have presented a pilot survey designed to uncover new X-ray luminous quasars at z > 5.6, combining eROSITA data with optical and IR photometry from DES, VHS and CatWISE2020. We have discovered 5 quasars out of 9 sources that were selected for spectroscopic confirmation with LDSS3-C. Two of these sources are securely detected in eRASS. These sources are also radio-detected. J020916−562650 is formally radio-loud, while J202040−621509 lies right at the threshold delimiting radio-loud from radio-quiet sources (R = 10). The reionisation-era quasars uncovered by our pilot study shed light on powerful X-ray emission mechanisms in the early universe and give us access to demographic diagnostics for the z > 5.6 quasar population.

J020916−562650 (z = 5.6) shows X-ray and radio properties of a blazar, making it one of the most distant AGN of this type known to date (Belladitta et al. 2020; An et al. 2023; Caccianiga et al. 2024). In their independent analysis, Ighina et al. (2024) also come to the conclusion that this source is a blazar. Recently, the discovery of a blazar at z = 7 was reported by Banados et al. (2024), setting a new redshift frontier for this type of source. From this detection and using simple geometric arguments, we infered a space density of similarly luminous, jetted quasars of 6 . 89 5.70 15.86 Gpc 3 $ 6.89^{15.86}_{-5.70}\,\mathrm{Gpc^{-3}} $. This value is consistent with the space density of all quasars, jetted or not, as measured by the most recent quasar luminosity functions at this redshift and brightness. Our measurement implies two potential scenarios: either all quasars discovered so far at z ∼ 6 and M1450 < −25 are jetted (and most of these jets have not been detected yet) or a substantial fraction of jetted quasars in this redshift-magnitude regime have eluded quasar searches so far, perhaps due to obscuration. Severe systematics that can affect our measurement of quasar space density are discussed in this paper.

J202040−621509, (z = 5.7) is an R = 10 BAL quasar. Through four consecutive eRASS observations and an additional Chandra pointing, we provided substantial evidence for an extreme X-ray dimming event (a factor ≥13 in X-ray flux). Potential dimming scenarios are nuclear obscuration or disk-dominated quasar SEDs in strongly outflowing quasars with rapidly accreting black holes. We conclude, however, that increased X-ray counting statistics and simultaneous UV monitoring will be required to draw clearer conclusions on this dimming event.

An eRASS hemisphere reionization quasar search will be based on this pilot study. With the advent of the next generation of wide-field optical and IR imaging surveys such as Euclid Wide Survey (Euclid Collaboration 2019, 2024), and surveys of the Rubin Observatory (Ivezić et al. 2019) and the Nancy Grace Roman Space Telescope, we will be able to explore new synergies with the eRASS, enabling the discovery of the most X-ray luminous quasars in the early universe.


1

In “eRASSn”, n refers to the number of all-sky scans, while “eRASS:n” corresponds to the cumulative survey up to the n-th scan.

2

In this work, DES always refers to the second data release of the survey.

3

We note that shortly after the submission of the present paper, the independent discovery of quasar J020916−562650 was reported by Ighina et al. (2024). The quasar J202040−621509 was also presented as quasar candidate by Ighina et al. (2023). However, our spectroscopic confirmation of this quasar also pre-dates this publication.

4

XSPEC 12.13.1.

6

These percentiles correspond to the μ ± 1σ locations of the normal distribution.

7

The CSIRO ASKAP Science Data Archive, https://data.csiro.au/domain/casda

Acknowledgments

JW acknowledges support by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC-2094 – 390783311. JW is grateful to Dr. Eduardo Bañados for guidance on the topic of high-redshift quasar searches. He would also like to thank Dr. Peter Predehl and Dr. Vadim Burwitz for helping with the Chandra observations presented in this project. T. Urrutia acknowledges funding from the ERC-AdG grant SPECMAP-CGM, GA 101020943. This work is based on data from eROSITA, the soft X-ray instrument aboard SRG, a joint Russian-German science mission supported by the Russian Space Agency (Roskosmos), in the interests of the Russian Academy of Sciences represented by its Space Research Institute (IKI), and the Deutsches Zentrum für Luft- und Raumfahrt (DLR). The SRG spacecraft was built by Lavochkin Association (NPOL) and its subcontractors, and is operated by NPOL with support from the Max Planck Institute for Extraterrestrial Physics (MPE). The development and construction of the eROSITA X-ray instrument was led by MPE, with contributions from the Dr. Karl Remeis Observatory Bamberg & ECAP (FAU Erlangen-Nuernberg), the University of Hamburg Observatory, the Leibniz Institute for Astrophysics Potsdam (AIP), and the Institute for Astronomy and Astrophysics of the University of Tübingen, with the support of DLR and the Max Planck Society. The Argelander Institute for Astronomy of the University of Bonn and the Ludwig Maximilians Universität Munich also participated in the science preparation for eROSITA. The eROSITA data shown here were processed using the eSASS/NRTA software system developed by the German eROSITA consortium. The scientific results reported in this article are based to a significant degree on observations made by the Chandra X-ray Observatory. his project used public archival data from the Dark Energy Survey (DES). Funding for the DES Projects has been provided by the U.S. Department of Energy, the U.S. National Science Foundation, the Ministry of Science and Education of Spain, the Science and Technology Facilities Council of the United Kingdom, the Higher Education Funding Council for England, the National Center for Supercomputing Applications at the University of Illinois at Urbana–Champaign, the Kavli Institute of Cosmological Physics at the University of Chicago, the Center for Cosmology and Astro-Particle Physics at the Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, Financiadora de Estudos e Projetos, Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Cientfico e Tecnológico and the Ministério da Ciência, Tecnologia e Inovação, the Deutsche Forschungsgemeinschaft and the Collaborating Institutions in the Dark Energy Survey. The Collaborating Institutions are Argonne National Laboratory, the University of California at Santa Cruz, the University of Cambridge, Centro de Investigaciones Enérgeticas, Medioambientales y Tecnológicas–Madrid, the University of Chicago, University College London, the DES-Brazil Consortium, the University of Edinburgh, the Eidgenössische Technische Hochschule (ETH) Zürich, Fermi National Accelerator Laboratory, the University of Illinois at Urbana-Champaign, the Institut de Ciències de l’Espai (IEEC/CSIC), the Institut de Física d’Altes Energies, Lawrence Berkeley National Laboratory, the Ludwig-Maximilians Universität München and the associated Excellence Cluster Universe, the University of Michigan, the National Optical Astronomy Observatory, the University of Nottingham, The Ohio State University, the OzDES Membership Consortium, the University of Pennsylvania, the University of Portsmouth, SLAC National Accelerator Laboratory, Stanford University, the University of Sussex, and Texas A&M University. Based in part on observations at Cerro Tololo Inter-American Observatory, National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. CatWISE uses data products from WISE, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory (JPL)/California Institute of Technology (Caltech), funded by the National Aeronautics and Space Administration (NASA), and from NEOWISE, which is a joint project of JPL/Caltech and the University of Arizona funded by NASA. CatWISE is led by JPL/Caltech, with funding from NASA’s Astrophysics Data Analysis Program. Based on observations obtained as part of the VISTA Hemisphere Survey, ESO Progam, 179.A-2010 (PI: McMahon) The Legacy Surveys consist of three individual and complementary projects: the Dark Energy Camera Legacy Survey (DECaLS; Proposal ID #2014B-0404; PIs: David Schlegel and Arjun Dey), the Beijing-Arizona Sky Survey (BASS; NOAO Prop. ID #2015A-0801; PIs: Zhou Xu and Xiaohui Fan), and the Mayall z-band Legacy Survey (MzLS; Prop. ID #2016A-0453; PI: Arjun Dey). DECaLS, BASS and MzLS together include data obtained, respectively, at the Blanco telescope, Cerro Tololo Inter-American Observatory, NSF’s NOIRLab; the Bok telescope, Steward Observatory, University of Arizona; and the Mayall telescope, Kitt Peak National Observatory, NOIRLab. Pipeline processing and analyses of the data were supported by NOIRLab and the Lawrence Berkeley National Laboratory (LBNL). The Legacy Surveys project is honored to be permitted to conduct astronomical research on Iolkam Du’ag (Kitt Peak), a mountain with particular significance to the Tohono O’odham Nation. NOIRLab is operated by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. LBNL is managed by the Regents of the University of California under contract to the U.S. Department of Energy. This project used data obtained with the Dark Energy Camera (DECam), which was constructed by the Dark Energy Survey (DES) collaboration. Funding for the DES Projects has been provided by the U.S. Department of Energy, the U.S. National Science Foundation, the Ministry of Science and Education of Spain, the Science and Technology Facilities Council of the United Kingdom, the Higher Education Funding Council for England, the National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign, the Kavli Institute of Cosmological Physics at the University of Chicago, Center for Cosmology and Astro-Particle Physics at the Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, Financiadora de Estudos e Projetos, Fundacao Carlos Chagas Filho de Amparo, Financiadora de Estudos e Projetos, Fundacao Carlos Chagas Filho de Amparo a Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Cientifico e Tecnologico and the Ministerio da Ciencia, Tecnologia e Inovacao, the Deutsche Forschungsgemeinschaft and the Collaborating Institutions in the Dark Energy Survey. The Collaborating Institutions are Argonne National Laboratory, the University of California at Santa Cruz, the University of Cambridge, Centro de Investigaciones Energeticas, Medioambientales y Tecnologicas-Madrid, the University of Chicago, University College London, the DES-Brazil Consortium, the University of Edinburgh, the Eidgenossische Technische Hochschule (ETH) Zurich, Fermi National Accelerator Laboratory, the University of Illinois at Urbana-Champaign, the Institut de Ciencies de l’Espai (IEEC/CSIC), the Institut de Fisica d’Altes Energies, Lawrence Berkeley National Laboratory, the Ludwig Maximilians Universitat Munchen and the associated Excellence Cluster Universe, the University of Michigan, NSF’s NOIRLab, the University of Nottingham, the Ohio State University, the University of Pennsylvania, the University of Portsmouth, SLAC National Accelerator Laboratory, Stanford University, the University of Sussex, and Texas A&M University. BASS is a key project of the Telescope Access Program (TAP), which has been funded by the National Astronomical Observatories of China, the Chinese Academy of Sciences (the Strategic Priority Research Program “The Emergence of Cosmological Structures” Grant # XDB09000000), and the Special Fund for Astronomy from the Ministry of Finance. The BASS is also supported by the External Cooperation Program of Chinese Academy of Sciences (Grant # 114A11KYSB20160057), and Chinese National Natural Science Foundation (Grant # 12120101003, # 11433005). The Legacy Survey team makes use of data products from the Near-Earth Object Wide-field Infrared Survey Explorer (NEOWISE), which is a project of the Jet Propulsion Laboratory/California Institute of Technology. NEOWISE is funded by the National Aeronautics and Space Administration. The Legacy Surveys imaging of the DESI footprint is supported by the Director, Office of Science, Office of High Energy Physics of the U.S. Department of Energy under Contract No. DE-AC02-05CH1123, by the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility under the same contract; and by the U.S. National Science Foundation, Division of Astronomical Sciences under Contract No. AST-0950945 to NOAO. This research made use of Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration 2013, 2018). This research has made use of software provided by the Chandra X-ray Center (CXC) in the application packages CIAO and Sherpa. This research made use of Astropy (http://www.astropy.org), a community-developed core Python package for Astronomy (Astropy Collaboration 2013, 2018). This research made use of Photutils, an Astropy package for detection and photometry of astronomical sources (Bradley et al. 2023). This research has made use of “Aladin sky atlas” developed at CDS, Strasbourg Observatory, France. In addition this research made use of BXA (https://johannesbuchner.github.io/BXA, Buchner et al. 2014), the corner package (Foreman-Mackey 2016), Pymoc, Aladin (Bonnarel et al. 2000), topcat (Taylor 2005), lmfit (Newville et al. 2016), seaborn, cmocean, matplotlib.

References

  1. Abbott, T. M. C., Adamów, M., Aguena, M., et al. 2021, ApJS, 255, 20 [NASA ADS] [CrossRef] [Google Scholar]
  2. Albert, J., Aliu, E., Anderhub, H., et al. 2007, ApJ, 669, 862 [Google Scholar]
  3. An, T., Wang, A., Liu, Y., et al. 2023, MNRAS, 519, 4047 [CrossRef] [Google Scholar]
  4. Ananna, T. T., Salvato, M., LaMassa, S., et al. 2017, ApJ, 850, 66 [Google Scholar]
  5. Arcodia, R., Merloni, A., Comparat, J., et al. 2024, A&A, 681, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Arnouts, S., & Ilbert, O. 2011, Astrophysics Source Code Library [record ascl:1108.009] [Google Scholar]
  7. Astropy Collaboration (Robitaille, T. P., et al.) 2013, A&A, 558, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Astropy Collaboration (Price-Whelan, A. M., et al.) 2018, AJ, 156, 123 [Google Scholar]
  9. Bañados, E., Venemans, B. P., Morganson, E., et al. 2015, ApJ, 804, 118 [Google Scholar]
  10. Bañados, E., Venemans, B. P., Decarli, R., et al. 2016, ApJS, 227, 11 [Google Scholar]
  11. Bañados, E., Venemans, B. P., Mazzucchelli, C., et al. 2018, Nature, 553, 473 [Google Scholar]
  12. Bañados, E., Schindler, J.-T., Venemans, B. P., et al. 2023, ApJS, 265, 29 [CrossRef] [Google Scholar]
  13. Banados, E., Momjian, E., Connor, T., et al. 2024, ArXiv e-prints [arXiv:2407.07236] [Google Scholar]
  14. Barlow-Hall, C. L., Delaney, J., Aird, J., et al. 2023, MNRAS, 519, 6055 [CrossRef] [Google Scholar]
  15. Bayo, A., Barrado, D., Stauffer, J., et al. 2011, A&A, 536, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Belladitta, S., Moretti, A., Caccianiga, A., et al. 2019, A&A, 629, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Belladitta, S., Moretti, A., Caccianiga, A., et al. 2020, A&A, 635, L7 [EDP Sciences] [Google Scholar]
  18. Belladitta, S., Caccianiga, A., Diana, A., et al. 2022, A&A, 660, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Best, W. M. J., Liu, M. C., Magnier, E. A., et al. 2015, ApJ, 814, 118 [CrossRef] [Google Scholar]
  20. Best, W. M. J., Magnier, E. A., Liu, M. C., et al. 2018, Am. Astron. Soc. Meet. Abstr., 231, 137.01 [NASA ADS] [Google Scholar]
  21. Blandford, R. D., & McKee, C. F. 1982, ApJ, 255, 419 [Google Scholar]
  22. Bohlin, R. C., Colina, L., & Finley, D. S. 1995, AJ, 110, 1316 [NASA ADS] [CrossRef] [Google Scholar]
  23. Boksenberg, A., Snijders, M. A. J., Wilson, R., et al. 1978, Nature, 275, 404 [NASA ADS] [CrossRef] [Google Scholar]
  24. Bonnarel, F., Fernique, P., Bienaymé, O., et al. 2000, A&AS, 143, 33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Bradley, L., Sipőcz, B., Robitaille, T., et al. 2023, https://doi.org/10.5281/zenodo.7946442 [Google Scholar]
  26. Brandt, W. N., Laor, A., & Wills, B. J. 2000, ApJ, 528, 637 [CrossRef] [Google Scholar]
  27. Briggs, K. R., & Pye, J. P. 2004, MNRAS, 353, 673 [NASA ADS] [CrossRef] [Google Scholar]
  28. Brown, M. J. I., Duncan, K. J., Landt, H., et al. 2019, MNRAS, 489, 3351 [Google Scholar]
  29. Brunner, H., Liu, T., Lamer, G., et al. 2022, A&A, 661, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Brusa, M., Zamorani, G., Comastri, A., et al. 2007, ApJS, 172, 353 [NASA ADS] [CrossRef] [Google Scholar]
  31. Buchner, J., Georgakakis, A., Nandra, K., et al. 2014, A&A, 564, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Budavári, T., & Szalay, A. S. 2008, ApJ, 679, 301 [CrossRef] [Google Scholar]
  33. Bunker, A. J., Saxena, A., Cameron, A. J., et al. 2023, A&A, 677, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Burningham, B., Pinfield, D. J., Lucas, P. W., et al. 2010, MNRAS, 406, 1885 [NASA ADS] [Google Scholar]
  35. Burningham, B., Cardoso, C. V., Smith, L., et al. 2013, MNRAS, 433, 457 [CrossRef] [Google Scholar]
  36. Caballero, J. A., Burgasser, A. J., & Klement, R. 2008, A&A, 488, 181 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Caccianiga, A., Ighina, L., Moretti, A., et al. 2024, A&A, 684, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Calderone, G., Ghisellini, G., Colpi, M., & Dotti, M. 2013, MNRAS, 431, 210 [NASA ADS] [CrossRef] [Google Scholar]
  39. Caramazza, M., Stelzer, B., Magaudda, E., et al. 2023, A&A, 676, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Cardoso, C. V., Burningham, B., Smart, R. L., et al. 2015, MNRAS, 450, 2486 [NASA ADS] [CrossRef] [Google Scholar]
  41. Cash, W. 1979, ApJ, 228, 939 [Google Scholar]
  42. Chabrier, G. 2001, ApJ, 554, 1274 [NASA ADS] [CrossRef] [Google Scholar]
  43. Chabrier, G., Baraffe, I., Allard, F., & Hauschildt, P. 2000, ApJ, 542, 464 [Google Scholar]
  44. Chehade, B., Carnall, A. C., Shanks, T., et al. 2018, MNRAS, 478, 1649 [Google Scholar]
  45. Cherepashchuk, A. M., & Lyutyi, V. M. 1973, Astrophys. Lett., 13, 165 [NASA ADS] [Google Scholar]
  46. Chiu, K., Fan, X., Leggett, S. K., et al. 2006, AJ, 131, 2722 [Google Scholar]
  47. Chiu, K., Liu, M. C., Jiang, L., et al. 2008, MNRAS, 385, L53 [CrossRef] [Google Scholar]
  48. Colless, M., Dalton, G., Maddox, S., et al. 2001, MNRAS, 328, 1039 [Google Scholar]
  49. Cruz, K. L., Reid, I. N., Liebert, J., Kirkpatrick, J. D., & Lowrance, P. J. 2003, AJ, 126, 2421 [Google Scholar]
  50. Cruz, K. L., Reid, I. N., Kirkpatrick, J. D., et al. 2007, AJ, 133, 439 [NASA ADS] [CrossRef] [Google Scholar]
  51. Dai, X., Griffin, R. D., Kochanek, C. S., Nugent, J. M., & Bregman, J. N. 2015, ApJS, 218, 8 [NASA ADS] [CrossRef] [Google Scholar]
  52. Dark Energy Survey Collaboration (Abbott, T., et al.) 2016, MNRAS, 460, 1270 [Google Scholar]
  53. Day-Jones, A. C., Marocco, F., Pinfield, D. J., et al. 2013, MNRAS, 430, 1171 [NASA ADS] [CrossRef] [Google Scholar]
  54. Decarli, R., Walter, F., Venemans, B. P., et al. 2018, ApJ, 854, 97 [Google Scholar]
  55. Dennerl, K., Burkert, W., Burwitz, V., et al. 2012, SPIE Conf. Ser., 8443, 844350 [Google Scholar]
  56. Dey, A., Schlegel, D. J., Lang, D., et al. 2019, AJ, 157, 168 [Google Scholar]
  57. Donato, D., Ghisellini, G., Tagliaferri, G., & Fossati, G. 2001, A&A, 375, 739 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Eilers, A.-C., Hennawi, J. F., Decarli, R., et al. 2020, ApJ, 900, 37 [Google Scholar]
  59. Euclid Collaboration (Barnett, R., et al.) 2019, A&A, 631, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  60. Euclid Collaboration (Mellier, Y., et al.) 2024, A&A, in press, https://doi.org/10.1051/0004-6361/202450810 [Google Scholar]
  61. Fan, X., Narayanan, V. K., Lupton, R. H., et al. 2001, AJ, 122, 2833 [Google Scholar]
  62. Fan, X., Bañados, E., & Simcoe, R. A. 2023, ARA&A, 61, 373 [NASA ADS] [CrossRef] [Google Scholar]
  63. Foreman-Mackey, D. 2016, J. Open Source Softw., 1, 24 [Google Scholar]
  64. Gallagher, S. C., Brandt, W. N., Chartas, G., Garmire, G. P., & Sambruna, R. M. 2002, ApJ, 569, 655 [NASA ADS] [CrossRef] [Google Scholar]
  65. Gallagher, S. C., Brandt, W. N., Chartas, G., et al. 2006, ApJ, 644, 709 [Google Scholar]
  66. Gehrels, N. 1986, ApJ, 303, 336 [Google Scholar]
  67. Ghisellini, G., Della Ceca, R., Volonteri, M., et al. 2010, MNRAS, 405, 387 [NASA ADS] [Google Scholar]
  68. Ghisellini, G., Sbarrato, T., Tagliaferri, G., et al. 2014, MNRAS, 440, L111 [CrossRef] [Google Scholar]
  69. Ghisellini, G., Haardt, F., Ciardi, B., et al. 2015, MNRAS, 452, 3457 [Google Scholar]
  70. Ghisellini, G., Righi, C., Costamante, L., & Tavecchio, F. 2017, MNRAS, 469, 255 [NASA ADS] [CrossRef] [Google Scholar]
  71. Gibson, R. R., & Brandt, W. N. 2012, ApJ, 746, 54 [NASA ADS] [CrossRef] [Google Scholar]
  72. Gibson, R. R., Brandt, W. N., Gallagher, S. C., & Schneider, D. P. 2009, ApJ, 696, 924 [NASA ADS] [CrossRef] [Google Scholar]
  73. Giustini, M., & Proga, D. 2019, A&A, 630, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Gloudemans, A. J., Duncan, K. J., Saxena, A., et al. 2022, A&A, 668, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Golimowski, D. A., Leggett, S. K., Marley, M. S., et al. 2004, AJ, 127, 3516 [Google Scholar]
  76. Goulding, A. D., Greene, J. E., Setton, D. J., et al. 2023, ApJ, 955, L24 [NASA ADS] [CrossRef] [Google Scholar]
  77. Green, P. J., Aldcroft, T. L., Mathur, S., Wilkes, B. J., & Elvis, M. 2001, ApJ, 558, 109 [CrossRef] [Google Scholar]
  78. Haardt, F., & Maraschi, L. 1991, ApJ, 380, L51 [Google Scholar]
  79. Habouzit, M., Onoue, M., Bañados, E., et al. 2022, MNRAS, 511, 3751 [NASA ADS] [CrossRef] [Google Scholar]
  80. Hale, C. L., McConnell, D., Thomson, A. J. M., et al. 2021, PASA, 38, e058 [NASA ADS] [CrossRef] [Google Scholar]
  81. Harikane, Y., Zhang, Y., Nakajima, K., et al. 2023, ApJ, 959, 39 [NASA ADS] [CrossRef] [Google Scholar]
  82. Hawley, S. L., Covey, K. R., Knapp, G. R., et al. 2002, AJ, 123, 3409 [NASA ADS] [CrossRef] [Google Scholar]
  83. HI4PI Collaboration (Ben Bekhti, N., et al.) 2016, A&A, 594, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  84. Hildebrandt, H., Arnouts, S., Capak, P., et al. 2010, A&A, 523, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  85. Hotan, A. W., Bunton, J. D., Chippendale, A. P., et al. 2021, PASA, 38, e009 [NASA ADS] [CrossRef] [Google Scholar]
  86. Ighina, L., Caccianiga, A., Moretti, A., et al. 2019, MNRAS, 489, 2732 [NASA ADS] [CrossRef] [Google Scholar]
  87. Ighina, L., Caccianiga, A., Moretti, A., et al. 2021, MNRAS, 505, 4120 [NASA ADS] [CrossRef] [Google Scholar]
  88. Ighina, L., Caccianiga, A., Moretti, A., et al. 2023, MNRAS, 519, 2060 [Google Scholar]
  89. Ighina, L., Caccianiga, A., Moretti, A., et al. 2024, A&A, 687, A242 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  90. Inayoshi, K., Visbal, E., & Haiman, Z. 2020, ARA&A, 58, 27 [NASA ADS] [CrossRef] [Google Scholar]
  91. Ishibashi, W., & Courvoisier, T. J. L. 2010, A&A, 512, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  92. Ivezić, Ž., Kahn, S. M., Tyson, J. A., et al. 2019, ApJ, 873, 111 [Google Scholar]
  93. Jiang, L., Fan, X., Annis, J., et al. 2008, AJ, 135, 1057 [NASA ADS] [CrossRef] [Google Scholar]
  94. Jiang, L., McGreer, I. D., Fan, X., et al. 2016, ApJ, 833, 222 [Google Scholar]
  95. Just, D. W., Brandt, W. N., Shemmer, O., et al. 2007, ApJ, 665, 1004 [Google Scholar]
  96. Kellermann, K. I., Sramek, R., Schmidt, M., et al. 1989, AJ, 98, 1195 [NASA ADS] [CrossRef] [Google Scholar]
  97. Khorunzhev, G. A., Meshcheryakov, A. V., Medvedev, P. S., et al. 2021, Astron. Lett., 47, 123 [NASA ADS] [CrossRef] [Google Scholar]
  98. Kirkpatrick, J. D., Looper, D. L., Burgasser, A. J., et al. 2010, ApJS, 190, 100 [Google Scholar]
  99. Kirkpatrick, J. D., Cushing, M. C., Gelino, C. R., et al. 2011, ApJS, 197, 19 [NASA ADS] [CrossRef] [Google Scholar]
  100. Knapp, G. R., Leggett, S. K., Fan, X., et al. 2004, AJ, 127, 3553 [Google Scholar]
  101. Knigge, C., Scaringi, S., Goad, M. R., & Cottis, C. E. 2008, MNRAS, 386, 1426 [NASA ADS] [CrossRef] [Google Scholar]
  102. Korista, K. T., Voit, G. M., Morris, S. L., & Weymann, R. J. 1993, ApJS, 88, 357 [NASA ADS] [CrossRef] [Google Scholar]
  103. Kulkarni, G., Worseck, G., & Hennawi, J. F. 2019, MNRAS, 488, 1035 [Google Scholar]
  104. LaMassa, S. M., Cales, S., Moran, E. C., et al. 2015, ApJ, 800, 144 [NASA ADS] [CrossRef] [Google Scholar]
  105. Langejahn, M., Kadler, M., Wilms, J., et al. 2020, A&A, 637, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  106. Larson, R. L., Finkelstein, S. L., Kocevski, D. D., et al. 2023, ApJ, 953, L29 [NASA ADS] [CrossRef] [Google Scholar]
  107. Li, J.-T., Wang, F., Yang, J., et al. 2021, ApJ, 906, 135 [NASA ADS] [CrossRef] [Google Scholar]
  108. Lira, P. 2021, in Nuclear Activity in Galaxies Across Cosmic Time, eds. M. Pović, P. Marziani, J. Masegosa, et al., 356, 101 [NASA ADS] [Google Scholar]
  109. Lister, M. L., Homan, D. C., Hovatta, T., et al. 2019, ApJ, 874, 43 [NASA ADS] [CrossRef] [Google Scholar]
  110. Liu, M. C., Wainscoat, R., Martín, E. L., Barris, B., & Tonry, J. 2002, ApJ, 568, L107 [NASA ADS] [CrossRef] [Google Scholar]
  111. Liu, T., Merloni, A., Comparat, J., et al. 2022, A&A, 661, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  112. Lusso, E., & Risaliti, G. 2016, ApJ, 819, 154 [Google Scholar]
  113. Madau, P., & Rees, M. J. 2000, ApJ, 542, L69 [Google Scholar]
  114. Magaudda, E., Stelzer, B., Raetz, S., et al. 2022, A&A, 661, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  115. Maiolino, R., Scholtz, J., Curtis-Lake, E., et al. 2024, A&A, in press, https://doi.org/10.1051/0004-6361/202347640 [Google Scholar]
  116. Marocco, F., Andrei, A. H., Smart, R. L., et al. 2013, AJ, 146, 161 [NASA ADS] [CrossRef] [Google Scholar]
  117. Marocco, F., Day-Jones, A. C., Lucas, P. W., et al. 2014, MNRAS, 439, 372 [NASA ADS] [CrossRef] [Google Scholar]
  118. Marocco, F., Jones, H. R. A., Day-Jones, A. C., et al. 2015, MNRAS, 449, 3651 [NASA ADS] [CrossRef] [Google Scholar]
  119. Marocco, F., Eisenhardt, P. R. M., Fowler, J. W., et al. 2021, ApJS, 253, 8 [Google Scholar]
  120. Marshall, N., Warwick, R. S., & Pounds, K. A. 1981, MNRAS, 194, 987 [NASA ADS] [CrossRef] [Google Scholar]
  121. Matsuoka, Y., Strauss, M. A., Kashikawa, N., et al. 2018, ApJ, 869, 150 [Google Scholar]
  122. Matsuoka, Y., Iwasawa, K., Onoue, M., et al. 2022, ApJS, 259, 18 [NASA ADS] [CrossRef] [Google Scholar]
  123. Matsuoka, Y., Onoue, M., Iwasawa, K., et al. 2023, ApJ, 949, L42 [NASA ADS] [CrossRef] [Google Scholar]
  124. McConnell, D., Hale, C. L., Lenc, E., et al. 2020, PASA, 37 [CrossRef] [Google Scholar]
  125. McHardy, I. M. 2013, MNRAS, 430, L49 [NASA ADS] [CrossRef] [Google Scholar]
  126. McHardy, I. M., Papadakis, I. E., Uttley, P., Page, M. J., & Mason, K. O. 2004, MNRAS, 348, 783 [NASA ADS] [CrossRef] [Google Scholar]
  127. McMahon, R. G., Banerji, M., Gonzalez, E., et al. 2013, The Messenger, 154, 35 [NASA ADS] [Google Scholar]
  128. McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, ASP Conf. Ser., 376, 127 [Google Scholar]
  129. Medvedev, P., Sazonov, S., Gilfanov, M., et al. 2020, MNRAS, 497, 1842 [Google Scholar]
  130. Merloni, A., Lamer, G., Liu, T., et al. 2024, A&A, 682, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  131. Miralda-Escudé, J. 1998, ApJ, 501, 15 [CrossRef] [Google Scholar]
  132. Mortlock, D. J., Patel, M., Warren, S. J., et al. 2009, A&A, 505, 97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  133. Mortlock, D. J., Warren, S. J., Venemans, B. P., et al. 2011, Nature, 474, 616 [Google Scholar]
  134. Murphy, T., Chatterjee, S., Kaplan, D. L., et al. 2013, PASA, 30, e006 [NASA ADS] [CrossRef] [Google Scholar]
  135. Murray, N., Chiang, J., Grossman, S. A., & Voit, G. M. 1995, ApJ, 451, 498 [NASA ADS] [CrossRef] [Google Scholar]
  136. Nandra, K., & Pounds, K. A. 1994, MNRAS, 268, 405 [Google Scholar]
  137. Nandra, K., George, I. M., Mushotzky, R. F., Turner, T. J., & Yaqoob, T. 1997, ApJ, 476, 70 [Google Scholar]
  138. Nanni, R., Vignali, C., Gilli, R., Moretti, A., & Brandt, W. N. 2017, A&A, 603, A128 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  139. Nanni, R., Gilli, R., Vignali, C., et al. 2018, A&A, 614, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  140. Newville, M., Stensitzki, T., Allen, D. B., et al. 2016, Astrophysics Source Code Library [record ascl:1606.014] [Google Scholar]
  141. Noll, S., Mehlert, D., Appenzeller, I., et al. 2004, A&A, 418, 885 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  142. Norris, R. P., Marvil, J., Collier, J. D., et al. 2021, PASA, 38 [Google Scholar]
  143. Paolillo, M., Papadakis, I. E., Brandt, W. N., et al. 2023, A&A, 673, A68 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  144. Pâris, I., Petitjean, P., Rollinde, E., et al. 2011, A&A, 530, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  145. Peterson, B. M. 1993, PASP, 105, 247 [NASA ADS] [CrossRef] [Google Scholar]
  146. Pickles, A. J. 1998, PASP, 110, 863 [NASA ADS] [CrossRef] [Google Scholar]
  147. Pineau, F. X., Derriere, S., Motch, C., et al. 2017, A&A, 597, A89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  148. Planck Collaboration VI. 2020, A&A, 641, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  149. Polletta, M., Tajer, M., Maraschi, L., et al. 2007, ApJ, 663, 81 [NASA ADS] [CrossRef] [Google Scholar]
  150. Pons, E., McMahon, R. G., Banerji, M., & Reed, S. L. 2020, MNRAS, 491, 3884 [Google Scholar]
  151. Ponti, G., Papadakis, I., Bianchi, S., et al. 2012, A&A, 542, A83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  152. Predehl, P., Andritschke, R., Arefiev, V., et al. 2021, A&A, 647, A1 [EDP Sciences] [Google Scholar]
  153. Reed, S. L., McMahon, R. G., Banerji, M., et al. 2015, MNRAS, 454, 3952 [Google Scholar]
  154. Reed, S. L., McMahon, R. G., Martini, P., et al. 2017, MNRAS, 468, 4702 [NASA ADS] [CrossRef] [Google Scholar]
  155. Reed, S. L., Banerji, M., Becker, G. D., et al. 2019, MNRAS, 487, 1874 [Google Scholar]
  156. Reeves, J. N., Turner, M. J. L., Ohashi, T., & Kii, T. 1997, MNRAS, 292, 468 [NASA ADS] [Google Scholar]
  157. Reid, I. N., Cruz, K. L., Kirkpatrick, J. D., et al. 2008, AJ, 136, 1290 [Google Scholar]
  158. Rogerson, J. A., Hall, P. B., Rodríguez Hidalgo, P., et al. 2016, MNRAS, 457, 405 [NASA ADS] [CrossRef] [Google Scholar]
  159. Rutledge, R. E., Basri, G., Martín, E. L., & Bildsten, L. 2000, ApJ, 538, L141 [NASA ADS] [CrossRef] [Google Scholar]
  160. Rybicki, G. B., & Lightman, A. P. 1979, Radiative Processes in Astrophysics (New York: Wiley) [Google Scholar]
  161. Saez, C., Brandt, W. N., Gallagher, S. C., Bauer, F. E., & Garmire, G. P. 2012, ApJ, 759, 42 [NASA ADS] [CrossRef] [Google Scholar]
  162. Salvato, M., Hasinger, G., Ilbert, O., et al. 2009, ApJ, 690, 1250 [CrossRef] [Google Scholar]
  163. Salvato, M., Buchner, J., Budavári, T., et al. 2018, MNRAS, 473, 4937 [Google Scholar]
  164. Salvato, M., Wolf, J., Dwelly, T., et al. 2022, A&A, 661, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  165. Sbarrato, T., Ghisellini, G., Maraschi, L., & Colpi, M. 2012, MNRAS, 421, 1764 [Google Scholar]
  166. Sbarrato, T., Ghisellini, G., Giovannini, G., & Giroletti, M. 2021, A&A, 655, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  167. Scaringi, S., Cottis, C. E., Knigge, C., & Goad, M. R. 2009, MNRAS, 399, 2231 [NASA ADS] [CrossRef] [Google Scholar]
  168. Schindler, J.-T., Bañados, E., Connor, T., et al. 2023, ApJ, 943, 67 [NASA ADS] [CrossRef] [Google Scholar]
  169. Schmidt, S. J., West, A. A., Hawley, S. L., & Pineda, J. S. 2010, AJ, 139, 1808 [NASA ADS] [CrossRef] [Google Scholar]
  170. Seppi, R., Comparat, J., Bulbul, E., et al. 2022, A&A, 665, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  171. Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 [NASA ADS] [Google Scholar]
  172. Shemmer, O., Brandt, W. N., Vignali, C., et al. 2005, ApJ, 630, 729 [Google Scholar]
  173. Shen, Y., Wu, J., Jiang, L., et al. 2019, ApJ, 873, 35 [Google Scholar]
  174. Steffen, A. T., Strateva, I., Brandt, W. N., et al. 2006, AJ, 131, 2826 [Google Scholar]
  175. Stelzer, B., Marino, A., Micela, G., López-Santiago, J., & Liefke, C. 2013, MNRAS, 431, 2063 [Google Scholar]
  176. Stelzer, B., Klutsch, A., Coffaro, M., Magaudda, E., & Salvato, M. 2022, A&A, 661, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  177. Sunyaev, R., Arefiev, V., Babyshkin, V., et al. 2021, A&A, 656, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  178. Sutherland, W., & Saunders, W. 1992, MNRAS, 259, 413 [Google Scholar]
  179. Taylor, M. B. 2005, ASP Conf. Ser., 347, 29 [Google Scholar]
  180. Timlin, J. D., Brandt, W. N., Zhu, S., et al. 2020, MNRAS, 498, 4033 [NASA ADS] [CrossRef] [Google Scholar]
  181. Tinney, C. G., Kirkpatrick, J. D., Faherty, J. K., et al. 2018, ApJS, 236, 28 [Google Scholar]
  182. Trenti, M., & Stiavelli, M. 2008, ApJ, 676, 767 [Google Scholar]
  183. Urry, C. M., & Padovani, P. 1995, PASP, 107, 803 [NASA ADS] [CrossRef] [Google Scholar]
  184. van Haarlem, M. P., Wise, M. W., Gunst, A. W., et al. 2013, A&A, 556, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  185. Vanden Berk, D. E., Richards, G. T., Bauer, A., et al. 2001, AJ, 122, 549 [Google Scholar]
  186. Venemans, B. P., Walter, F., Neeleman, M., et al. 2020, ApJ, 904, 130 [Google Scholar]
  187. Vito, F., Brandt, W. N., Yang, G., et al. 2018, MNRAS, 473, 2378 [Google Scholar]
  188. Vito, F., Brandt, W. N., Bauer, F. E., et al. 2019, A&A, 630, A118 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  189. Vito, F., Mignoli, M., Gilli, R., et al. 2022, A&A, 663, A159 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  190. Volonteri, M., Haardt, F., Ghisellini, G., & Della Ceca, R. 2011, MNRAS, 416, 216 [Google Scholar]
  191. Wachter, K., Leach, R., & Kellogg, E. 1979, ApJ, 230, 274 [NASA ADS] [CrossRef] [Google Scholar]
  192. Wahhaj, Z., Liu, M. C., Biller, B. A., et al. 2011, ApJ, 729, 139 [NASA ADS] [CrossRef] [Google Scholar]
  193. Wang, F., Yang, J., Fan, X., et al. 2021a, ApJ, 907, L1 [Google Scholar]
  194. Wang, F., Fan, X., Yang, J., et al. 2021b, ApJ, 908, 53 [Google Scholar]
  195. Webb, N. A., Coriat, M., Traulsen, I., et al. 2020, A&A, 641, A136 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  196. Weisskopf, M. C., Wu, K., Trimble, V., et al. 2007, ApJ, 657, 1026 [Google Scholar]
  197. Weymann, R. J., Morris, S. L., Foltz, C. B., & Hewett, P. C. 1991, ApJ, 373, 23 [NASA ADS] [CrossRef] [Google Scholar]
  198. Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 [Google Scholar]
  199. Winkler, P. F., & White, A. E. 1975, ApJ, 199, L139 [NASA ADS] [CrossRef] [Google Scholar]
  200. Wolf, J., Nandra, K., Salvato, M., et al. 2021, A&A, 647, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  201. Wolf, J., Nandra, K., Salvato, M., et al. 2023, A&A, 669, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  202. Yang, J., Fan, X., Wu, X.-B., et al. 2017, AJ, 153, 184 [NASA ADS] [CrossRef] [Google Scholar]
  203. Yang, J., Wang, F., Fan, X., et al. 2021, ApJ, 923, 262 [NASA ADS] [CrossRef] [Google Scholar]
  204. Yang, D.-M., Schindler, J.-T., Nanni, R., et al. 2024, MNRAS, 528, 2679 [CrossRef] [Google Scholar]
  205. Zappacosta, L., Piconcelli, E., Fiore, F., et al. 2023, A&A, 678, A201 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

Appendix A: Example DES images and forced photometry: Discarded quasar candidates

In Fig. A.1, we present example DES cutouts of candidates that were selected and candidates that were rejected after visual inspection and forced photometry. While the first example is a close pair, the second and the third are affected by artifacts such as bleed-trails and missing images.

thumbnail Fig. A.1.

Discarded candidates after visual inspection and forced photometry.

Appendix B: Contaminant spectra and coordinates

The coordinates of the contaminants observed with LDSS3 are listed in Table B.1.

Table B.1.

Coordinates of the contaminants that were spectroscopically confirmed during the calibration of the quasar discovery pipeline.

Appendix C: Radio flux densities and radio images

Two of the newly discovered quasars are detected at radio frequencies in several surveys carried out with the Australian SKA Pathfinder (ASKAP, see the text in Section 4 for all the details). Here, we list the available flux densities and radio images for these two objects (J2020−6215 and J0209−5626).

Table C.1.

Summary of the archival radio observations of J202040−621509.

Table C.2.

Summary of the archival radio observations of J020916−562650.

thumbnail Fig. C.1.

Observed radio spectral energy distribution of J202040−621509 (left and J020916−562650 (right from 0.887 to 1.37 GHz. The estimated flux densities are reported in the legend, together with the date of observation and the value of the best fit spectral index. The radio data of J020916−562650 at 0.887 GHz suggest a hint of variability. Simultaneous radio observations are crucial to reveal the true radio spectral energy distribution.

thumbnail Fig. C.2.

4′ × 4′ cutouts radio images of J2020−6215 (first and second row) and J0209−5626 (third and fourth row), ordered in increasing frequency. Both sources are clearly detected in all the images. The red cross represents the position of the optical source, and the synthesised beam of each observation is shown in the bottom left-hand corner.

All Tables

Table 1.

Number of detected quasar candidates in each eRASS.

Table 2.

Summary of the 5 newly discovered quasars.

Table 3.

eROSITA X-ray properties from aperture photometry.

Table 4.

Chandra ACIS-S photometry for three sources detected within 30″ of the quasar.

Table B.1.

Coordinates of the contaminants that were spectroscopically confirmed during the calibration of the quasar discovery pipeline.

Table C.1.

Summary of the archival radio observations of J202040−621509.

Table C.2.

Summary of the archival radio observations of J020916−562650.

All Figures

thumbnail Fig. 1.

Colour-selection diagram. The hashed and grey regions correspond to the colour space from which we pre-selected high-redshift quasar candidates for the present pilot survey. Gaussian Kernel density estimates of DES detected M-dwarfs (Best et al. 2018) and L/T-dwarfs (Liu et al. 2002; Hawley et al. 2002; Cruz et al. 2003, 2007; Chiu et al. 2008; Reid et al. 2008; Schmidt et al. 2010; Burningham et al. 2010, 2013; Kirkpatrick et al. 2010, 2011; Wahhaj et al. 2011; Day-Jones et al. 2013; Marocco et al. 2013, 2014, 2015; Best et al. 2015; Cardoso et al. 2015; Tinney et al. 2018) are displayed as shaded green and shaded blue areas, respectively. DES detected, spectroscopically confirmed z > 5.6 quasars compiled in Fan et al. (2023) and additional sources from Yang et al. (2024) are shown as salmon red circles. 1D colour distributions of the samples from the literature are shown in the marginal histograms. Magenta stars indicate new quasars discovered in this work. Four of these discoveries are located at the edge of the grey selection box, close to the contamination region of brown dwarfs.

In the text
thumbnail Fig. 2.

Le PHARE test run on previously spectroscopically confirmed z > 5.6 quasars. Upper panel: The estimated photometric redshift zphot is compared to the true spectroscopic value zspec. The grey shaded area indicates the region within which Δz < 0.15 × (1 + zspec) where Δz = ∣zspec − zphot∣. Lower panel: The ratio of χ2 of the stellar and the AGN or galaxy template are shown in the lower panel. In only two cases, the stellar solution is preferred.

In the text
thumbnail Fig. 3.

Optical colour-plane of selected of 2605 zphot > 5.6 quasar candidates from our colour-cuts and SED fitting. The candidates (circles) are colour-coded according to their photometric redshift. We show redshift-colour tracks for Le PHARE templates: a typical quasar (pl_QSO_DR2_029_t0, an SDSS-DR5 composite modified by Salvato et al. 2009) and a Seyfert II template (Sey2_template_norm, Polletta et al. 2007), with increasing extinction values (0.0, 0.1, 0.2) displayed by increasing transparency. Squares indicate selected zphot > 5.6 candidates that are already confirmed in literature (references in Fan et al. 2023). The five quasars discovered in this work are shown as magenta stars.

In the text
thumbnail Fig. 4.

The intersection of the DES, VHS, CatWISE2020 and eRASS footprints (blue shaded area) covers a sky surface of ∼3500 deg2. The underlying eRASS1 source population is shown in grey. The quasars discovered in this work are marked as pink stars.

In the text
thumbnail Fig. 5.

Chance alignment. Upper panel: The cumulative distribution of expected random (unrelated) X-ray detections from forced photometry using a circular aperture of radius 30″ around 1604 optical quasar candidate positions. The black curve is the estimate obtained from Eq. (2). for a single eRASS survey. The black dashed curve shows an empirical estimate of the same number based on a random sample of positions in the surveyed field. The brown curve shows the number of expected contaminants obtained from a random sample in eRASS:4. Blue vertical lines indicate the Pb values of the best detection of quasars discovered in this work in any eRASS and in eRASS:4 (if available). Lower panel: Fraction of background sources fbkg expected to be detected in eRASS1 from a simulation by Seppi et al. (2022).

In the text
thumbnail Fig. 6.

Spectra of newly discovered eRASS quasars (left panel) and eRASS:4 (0.2–8 keV) postage stamps (right panel). The cutouts are 2′×2′ in size and are centred on the optical position of the quasars. The spectrum of J202040−621509 shows prominent broad absorption features. The continuum appears strongly absorbed over the entire spectral range bluewards of Lyman α. It appears in the eROSITA images as two distinct sources. This morphology is further investigated in Section 4.3.1 with a Chandra follow-up pointing. There, we argue that the eRASS source is related to the quasar (Fig. 11). The morphology J020916−562650 is more compact.

In the text
thumbnail Fig. 7.

Spectra of three newly discovered quasars selected during the calibration of the pipeline. These objects are not robustly detected in eROSITA. The red dashed line indicates the location of the Lyman α break.

In the text
thumbnail Fig. 8.

X-ray emission of the new eRASS quasars. Left panel: X-ray luminosity redshift distribution of X-ray luminous quasars from literature (Nanni et al. 2017; Vito et al. 2019; Pons et al. 2020; Wang et al. 2021b; Medvedev et al. 2020; Belladitta et al. 2020; Wolf et al. 2021, 2023; Zappacosta et al. 2023) spanning a redshift range 5.5 < z < 7.5 and a luminosity range 44 < log L2 − 10 keV/(erg/s) < 46.5. The two newly discovered quasars (displayed as pink stars). The colour gradient shows the sensitivity to a fiducial absorbed power-law with Galactic absorption of the final cumulative eRASS in the equatorial region (in fractional of sensitive area). The quasars discovered in this work lie at the luminous end of the quasar population in the early universe. We show the luminosity derived from the Chandra follow-up observation of J202040−621509. Between the eRASS1 and the Chandra observation, its luminosity has decreased by an order of magnitude. Right panel: αOXL2500 distribution of the same sample of sources. The 1σ confidence interval of relation derived by Nanni et al. (2017) is shown by red dotted line. The hexagonal pattern shows a lower redshift AGN sample by Lusso & Risaliti (2016). The three quasars X-ray detected and newly discovered are over-luminous in the X-ray wavebands. Following its dimming observed in the recent Chandra observation, J202040–621509 is perfectly consistent αOXL2500 relation.

In the text
thumbnail Fig. 9.

Posteriors on the photon index Γ and the X-ray luminosity L2 − 10 keV from the X-ray spectral analysis performed on the eRASS:4 spectrum of J020916−562650.

In the text
thumbnail Fig. 10.

The broadband rest frame SED of J020916−562650. A range of typical FSRQ SEDs (Lγ ∼ 1046 − 48 erg/s, grey shaded region) and a more tailored blazar analytical model (solid green line) are overplotted to the data. The radio and X-ray emission of the source are consistent with the typical blazar double hump SED. Specifically, the X-ray emission significantly overestimates the expected X-ray coronal emission, extrapolated by the source accretion emission that dominates the optical/UV restframe photometry (yellow dots). The thermal emission is modelled here with a standard accretion disk emission, IR torus, and X-ray Corona powered by a black hole with the mass of log MBH/M = 8.60 (blue solid and dashed lines, the latter showing the unabsorbed accretion disk SED).

In the text
thumbnail Fig. 11.

Astrometric configuration of J202040−621509. All images are centred on the optical position of the new quasar J202040−621509 and have sizes 100″ × 100″. Left panel: Broadband (0.5−7 keV) Chandra image. Chandra sources detected with 30″ (black, dashed circle) of the quasar coordinates with wavdetect are marked in red. The centroid of 1eRASS J202039.8−621525 catalogue detection is marked by a green triangle. The contours were obtained from the smoothed eRASS:4 0.5−7 keV image. Central panel: DES DR2 z-band image. Circular coloured markers denote counterpart candidates for 1eRASS J202039.8−621525 from the DESI Legacy Survey catalogue. They are colour-coded according to their relative probability p_i of being the best counterpart. The DESI Legacy Survey source with the highest p_i is spatially co-incident with the quasar coordinates. Right panel: CatWISE 2020 W1 image. Among the Chandra and eROSITA X-ray detections, only the quasar is unambiguously detected.

In the text
thumbnail Fig. 12.

Distribution of BXAF0.5 − 7 keV fluxes for the eRASS 1 (in blue) and Chandra observations (red) of J202040−621509. The vertical black lines show the locations of the 15.9-th and 84.1-th percentiles of the distributions. The green and pink regions span the 2.3-th and 97.7-th percentiles of the distributions. The lower panel shows the time separation in days between the two observations. We measured a decrease in the source flux between the two epochs significant at the ∼2σ level.

In the text
thumbnail Fig. A.1.

Discarded candidates after visual inspection and forced photometry.

In the text
thumbnail Fig. C.1.

Observed radio spectral energy distribution of J202040−621509 (left and J020916−562650 (right from 0.887 to 1.37 GHz. The estimated flux densities are reported in the legend, together with the date of observation and the value of the best fit spectral index. The radio data of J020916−562650 at 0.887 GHz suggest a hint of variability. Simultaneous radio observations are crucial to reveal the true radio spectral energy distribution.

In the text
thumbnail Fig. C.2.

4′ × 4′ cutouts radio images of J2020−6215 (first and second row) and J0209−5626 (third and fourth row), ordered in increasing frequency. Both sources are clearly detected in all the images. The red cross represents the position of the optical source, and the synthesised beam of each observation is shown in the bottom left-hand corner.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.