High-redshift Quasars Research Articles

Machine Learning–based Search of High-redshift Quasars

We present a machine learning search for high-redshift (5.0 < z < 6.5) quasars using the combined photometric data from the Dark Energy Spectroscopic Instrument (DESI) Imaging Legacy Surveys and the Wide-field Infrared Survey Explorer survey. We explore the imputation of missing values for high-redshift quasars, discuss the feature selections, compare different machine learning algorithms, and investigate the selections of class ensemble for the training sample, then we find that the random forest model is very effective in separating the high-redshift quasars from various contaminators. The 11 class random forest model can achieve a precision of 96.43% and a recall of 91.53% for high-redshift quasars for the test set. We demonstrate that the completeness of the high-redshift quasars can reach as high as 82.20%. The final catalog consists of 216,949 high-redshift quasar candidates with 476 high probable ones in the entire Legacy Surveys DR9 footprint, and we make the catalog publicly available. Using Multi Unit Spectroscopic Explorer (MUSE) and DESI early data release (EDR) public spectra, we find that 14 true high-redshift quasars (11 in the training sample) out of 21 candidates are correctly identified for MUSE, and 20 true high-redshift quasars (11 in the training sample) out of 21 candidates are correctly identified for DESI-EDR. Additionally, we estimate photometric redshift for the high-redshift quasar candidates using a random forest regression model with a high precision.

Reverberation Mapping of High-mass and High-redshift Quasars Using Gravitational Time Delays

Mass estimates of black holes (BHs) in the centers of active galactic nuclei (AGNs) often rely on the radius–luminosity relation. However, this relation, usually probed by reverberation mapping (RM), is poorly constrained in the high-luminosity and high-redshift ends due to the very long expected RM lag times. Multiply imaged AGNs may offer a unique opportunity to explore the radius–luminosity relation at these ends. In addition to comprising several magnified images enabling a more efficient light-curve sampling, the time delay between multiple images of strongly lensed quasars can also aid in making such RM measurements feasible on reasonable timescales: if the strong-lensing time delay is, for example, of the order of the expected RM time lag, changes in the emission lines in the leading image can be observed around the same time as the changes in the continuum in the trailing image. In this work we probe the typical time-delay distribution in galaxy-cluster lenses and estimate the number of both high-mass (∼109−1010 M ⊙) and high-redshift (z ≳ 4−12) quasars that are expected to be strongly lensed by clusters. We find that up to several tens of thousands of M BH ∼ 106–108 M ⊙ broad-line AGNs at z > 4 should be multiply imaged by galaxy clusters and detectable with JWST, hundreds with Euclid, and several thousand with the Roman Space Telescope, across the whole sky. These could supply an important calibration for the BH mass scaling in the early Universe.

Feedback-regulated seed black hole growth in star-forming molecular clouds and galactic nuclei

Context. The detection of supermassive black holes (SMBHs) in high-redshift luminous quasars may require a phase of rapid accretion, and as a precondition, substantial gas influx toward seed black holes (BHs) from kiloparsec or parsec scales. Our previous research demonstrated the plausibility of such gas supply for BH seeds within star-forming giant molecular clouds (GMCs) with high surface density (∼104 M⊙ pc−2), facilitating “hyper-Eddington” accretion via efficient feeding by dense clumps, which are driven by turbulence and stellar feedback. Aims. This article presents an investigation of the impacts of feedback from accreting BHs on this process, including radiation, mechanical jets, and highly relativistic cosmic rays. Methods. We ran a suite of numerical simulations to explore diverse parameter spaces of BH feedback, including the subgrid accretion model, feedback energy efficiency, mass loading factor, and initial metallicity. Results. Using radiative feedback models inferred from the slim disk, we find that hyper-Eddington accretion is still achievable, yielding BH bolometric luminosities of as high as 1041 − 1044 erg/s, depending on the GMC properties and specific feedback model assumed. We find that the maximum possible mass growth of seed BHs (ΔMmax BH) is regulated by the momentum-deposition rate from BH feedback, ṗfeedback/(ṀBHc), which leads to an analytic scaling that agrees well with simulations. This scenario predicts the rapid formation of ∼104 M⊙ intermediate-massive BHs (IMBHs) from stellar-mass BHs within ∼1 Myr. Furthermore, we examine the impacts of subgrid accretion models and how BH feedback may influence star formation within these cloud complexes.

On the Physical Nature of Lyα Transmission Spikes in High-redshift Quasar Spectra

We investigate Lyman-alpha (Lyα) transmission spikes at 5.2 < z < 6.8 using synthetic quasar spectra from the “Cosmic Reionization on Computers” simulations. We focus on understanding the relationship between these spikes and the properties of the intergalactic medium (IGM). Disentangling the complex interplay between IGM physics and the influence of galaxies on the generation of these spikes presents a significant challenge. To address this, we employ Explainable Boosting machines, an interpretable machine learning algorithm, to quantify the relative impact of various IGM properties on the Lyα flux. Our findings reveal that gas density is the primary factor influencing absorption strength, followed by the intensity of background radiation and the temperature of the IGM. Ionizing radiation from local sources (i.e., galaxies) appears to have a minimal effect on Lyα flux. The simulations show that transmission spikes predominantly occur in regions of low gas density. Our results challenge recent observational studies suggesting the origin of these spikes in regions with enhanced radiation. We demonstrate that Lyα transmission spikes are largely a product of the large-scale structure, of which galaxies are biased tracers.

EIGER. VI. The Correlation Function, Host Halo Mass, and Duty Cycle of Luminous Quasars at z ≳ 6

We expect luminous (M 1450 ≲ −26.5) high-redshift quasars to trace the highest-density peaks in the early Universe. Here, we present observations of four z ≳ 6 quasar fields using JWST/NIRCam in the imaging and wide-field slitless spectroscopy mode and report a wide range in the number of detected [O iii]-emitting galaxies in the quasars’ environments, ranging between a density enhancement of δ ≈ 65 within a 2 cMpc radius—one of the largest protoclusters during the Epoch of Reionization discovered to date—to a density contrast consistent with zero, indicating the presence of a UV-luminous quasar in a region comparable to the average density of the Universe. By measuring the two-point cross-correlation function of quasars and their surrounding galaxies, as well as the galaxy autocorrelation function, we infer a correlation length of quasars at 〈z〉 = 6.25 of r0QQ=22.0−2.9+3.0cMpch−1 , while we obtain a correlation length of the [O iii]-emitting galaxies of r0GG=4.1±0.3cMpch−1 . By comparing the correlation functions to dark-matter-only simulations we estimate the minimum mass of the quasars’ host dark matter halos to be log10(Mhalo,min/M⊙)=12.43−0.15+0.13 (and log10(Mhalo,min[OIII]/M⊙)=10.56−0.03+0.05 for the [O iii] emitters), indicating that (a) luminous quasars do not necessarily reside within the most overdense regions in the early Universe, and that (b) the UV-luminous duty cycle of quasar activity at these redshifts is f duty ≪ 1. Such short quasar activity timescales challenge our understanding of early supermassive black hole growth and provide evidence for highly dust-obscured growth phases or episodic, radiatively inefficient accretion rates.

A lack of Lyman α emitters within 5 Mpc of a luminous quasar in an overdensity at z = 6.9: Potential evidence of negative quasar feedback at protocluster scales

High-redshift quasars are thought to live in the densest regions of space, which should be made evident by an overdensity of galaxies around them. However, campaigns to identify these overdensities by searching for Lyman-break galaxies (LBGs) and Lyman α emitters (LAEs) have had mixed results. These may be explained by either the small field of view of some of the experiments, the broad redshift ranges targeted by LBG searches, and the inherently high uncertainty of quasar redshifts estimated from ultraviolet emission lines, which makes it difficult to place the Ly-α emission line within a narrowband filter. Here, we present a 3 square degree search (∼1000 pMpc2) for LAEs around the z = 6.9 quasar VIK J2348–3054 using the Dark Energy CAMera (DECam) housed on the 4m Blanco telescope, finding 38 LAEs. The systemic redshift of VIK J2348–3054 is known from ALMA [CII] observations and places the Ly-α emission line of companions within the NB964 narrowband of DECam. This is the largest field-of-view LAE search around a z &gt; 6 quasar conducted to date. We find that this field is ∼ten times more overdense than Chandra Deep-Field South, observed previously with the same instrumental setup as well as several combined blank fields. This is strong evidence that VIK J2348–3054 resides in an overdensity of LAEs over several Mpc. Surprisingly, we find a lack of LAEs within 5 physical Mpc of the quasar and take this to most likely be evidence of quasar-suppressing star formation in its immediate vicinity. This result highlights the importance of performing overdensity searches over large areas to properly assess the density of those regions of the Universe.

The radio properties of z&gt;3.5 quasars: Are most high-redshift radio-loud active galactic nuclei obscured?

We explore the radio properties of powerful (rest-frame luminosity $ at 500 MHz) high-redshift ($z 3.5$) quasars. The aim of this study is to gain a better understanding of radio-loud sources at the epoch when they reach the highest space density. We selected 29 radio-loud quasars at low radio frequencies (76 MHz). Their radio spectra, covering the range from 76 MHz to 5 GHz, are generally well reproduced by a single power law. We created samples that were matched in radio luminosity at lower redshift (from $z to $z to investigate any spectral evolution. We find that the fraction of flat-spectrum radio quasars (FSRQs) strongly increases with redshift (from $ 8&lt;!PCT!&gt;$ at z=1.2 to $ 45&lt;!PCT!&gt;$ at z$&gt;$3.5). This effect is also observed in quasars with lower luminosities (down to $ The increase in the fraction of FSRQs with redshift corresponds to a decrease in the steep- spectrum radio quasars. This result can be explained when we assume that the beaming factor and the slope of the luminosity function do not change with redshift and if high-redshift radio-loud sources can be recognized as quasars only when they are seen at a small viewing angle ($ but most of them, about 90&lt;!PCT!&gt; are obscured in the UV and optical bands. We also found a trend for the size of radio sources to decrease with increasing redshift. Because projection effects are insufficient to cause this trend, we suggest that the large amount of gas causing the nuclear obscuration also hampers the growth of the more distant sources.

Chronicling the Reionization History at 6 ≲ z ≲ 7 with Emergent Quasar Damping Wings

The spectra of high-redshift (z ≳ 6) quasars contain valuable information on the progression of the Epoch of Reionization. At redshifts z < 6, the observed Lyman-series forest shows that the intergalactic medium is nearly ionized, while at z > 7 the observed quasar damping wings indicate high neutral gas fractions. However, there remains a gap in neutral gas fraction constraints at 6 ≲ z ≲ 7 where the Lyman-series forest becomes saturated but damping wings have yet to fully emerge. In this work, we use a sample of 18 quasar spectra at redshifts 6.0 < z < 7.1 to close this gap. We apply neural networks to reconstruct the quasars’ continuum emission around the partially absorbed Lyα line to normalize their spectra, and stack these continuum-normalized spectra in three redshift bins. To increase the robustness of our results, we compare the stacks to a grid of models from two hydrodynamical simulations, ATON and CROC, and we measure the volume-averaged neutral gas fraction, x¯HI , while jointly fitting for the mean quasar lifetime, t Q, for each stacked spectrum. We chronicle the evolution of neutral gas fraction using the ATON (CROC) models as follows: x¯HI=0.21−0.07+0.17 ( x¯HI=0.10<10−4+0.73 ) at 〈z〉 = 6.10, x¯HI=0.21−0.07+0.33 ( x¯HI=0.57−0.47+0.26 ) at 〈z〉 = 6.46, and x¯HI=0.37−0.17+0.17 ( x¯HI=0.57−0.21+0.26 ) at 〈z〉 = 6.87. At the same time, we constrain the average quasar lifetime to be t Q ≲ 7 Myr across all redshift bins, in good agreement with previous studies.

Explanation of Time Dilation of High Redshift Quasars, Surface Brightness, and Cosmic Microwave Background with the Stress Cosmology

The Big Bang theory is believed to be based on three problems to the tired light model. In this report, “time dilation of high redshift quasars” is first explained with the stress cosmology. A proceeding (delaying) speed of time is shown as a logarithm of changed energy. Second, “surface brightness” relates to “time dilation” and the combined luminosity per unit time. It decreases with time dilation. Third, according to the stress cosmology, the “cosmic microwave background” is explained with a relation between movement distance and decreasing energy quantity of discharged light. Thus, three problems can be explained with the stress cosmology being part of the tired light model. Therefore, there is no absolute proof of the Big Bang theory. Moreover, there is a fatal contradiction relating to the first law of thermodynamics in the Big Bang theory. The Big Bang theory required that the universe be a closed system according to the first law of thermodynamics. Nevertheless, the ekpyrotic universe theory is utilized to explain the Big Bang. The first law of thermodynamics indicates that our universe was an open system. The Big Bang theory is optional.

Intermediate-mass Black Hole Progenitors from Stellar Collisions in Dense Star Clusters

Very massive stars (VMSs) formed via a sequence of stellar collisions in dense star clusters have been proposed as the progenitors of massive black hole seeds. VMSs could indeed collapse to form intermediate-mass black holes, which would then grow by accretion to become the supermassive black holes observed at the centers of galaxies and powering high-redshift quasars. Previous studies have investigated how different cluster initial conditions affect the formation of a VMS, including mass segregation, stellar collisions, and binaries, among others. In this study, we investigate the growth of VMSs with a new grid of Cluster Monte Carlo star cluster simulations—the most expansive to date. The simulations span a wide range of initial conditions, varying the number of stars, cluster density, stellar initial mass function (IMF), and primordial binary fraction. We find a gradual shift in the mass of the most massive collision product across the parameter space; in particular, denser clusters born with top-heavy IMFs provide strong collisional regimes that form VMSs with masses easily exceeding 1000 M ⊙. Our results are used to derive a fitting formula that can predict the typical mass of a VMS formed as a function of the star cluster properties. Additionally, we study the stochasticity of this process and derive a statistical distribution for the mass of the VMS formed in one of our models, recomputing the model 50 times with different initial random seeds.

Probing New physics with high-redshift quasars: axions and non-standard cosmology

The Hubble diagram of quasars, as candidates to “standardizable” candles, has been used to measure the expansion history of the Universe at late times, up to very high redshifts (z ∼ 7). It has been shown that this history, as inferred from the quasar dataset, deviates at ≳ 3σ level from the concordance (ΛCDM) cosmology model preferred by the cosmic microwave background (CMB) and other datasets. In this article, we investigate whether new physics beyond ΛCDM (BΛCDM) or beyond the Standard Model (BSM) could make the quasar data consistent with the concordance model. We first show that an effective redshift-dependent relation between the quasar UV and X-ray luminosities, complementing previous phenomenological work in the literature, can potentially remedy the discrepancy. Such a redshift dependence can be realized in a BSM model with axion-photon conversion in the intergalactic medium (IGM), although the preferred parameter space is in tension with various other astrophysical constraints on axions, at a level depending on the specific assumptions made regarding the IGM magnetic field. We briefly discuss a variation of the axion model that could evade these astrophysical constraints. On the other hand, we show that models beyond ΛCDM such as one with a varying dark energy equation of state (wCDM) or the phenomenological cosmographic model with a polynomial expansion of the luminosity distance, cannot alleviate the tension. The code for our analysis, based on emcee [1] and corner.py [2], is publicly available at github.com/ChenSun-Phys/high_z_candles.

Large-scale Overdensity of Lyman Break Galaxies around the z = 6.3 Ultraluminous Quasar J0100 + 2802

We study the environment of the z = 6.33 ultraluminous quasar SDSS J010013.02+280225.8 (J0100) to understand its association with large-scale structure. Theoretical models propose high-redshift quasars as markers of galaxy overdensities residing in the most massive dark matter halos (DMHs) in the early Universe. J0100 is an ultraluminous quasar with the most massive black hole known at z ≳ 6, suggesting a high likelihood of residing in a massive DMH. We present wide-field (∼522 arcmin2) imaging in the r, i, and z bands from the Large Binocular Cameras on the Large Binocular Telescope, with Y- and J-band imaging from the Wide-field Infrared Camera on the Canada–France–Hawaii Telescope, centered on J0100. Applying color selections, we identify 23 objects as i-dropout Lyman break galaxy (LBG) candidates in the J0100 field. We use the deep photometric catalog in the 1.27 deg2 COSMOS field to calculate the density of LBGs in a blank field, and to estimate the selection completeness and purity. The observed surface density of LBG candidates in the J0100 field corresponds to a galaxy overdensity of δ = 4 (at 8.4σ). This large-scale overdensity suggests that the ∼22 arcmin2 overdensity found by Kashino et al. using JWST data extends out to much larger scales. We calculate the angular autocorrelation function of the candidates and find a positive correlation on ≲10′ scales as well as evidence of asymmetries in their spatial distribution, further suggesting the direct detection of large-scale structure in the field of the ultraluminous quasar J0100.

High-Redshift Quasars at z ≥ 3: Radio Variability and MPS/GPS Candidates

We present a study of the radio variability of bright, S1.4≥100 mJy, high-redshift quasars at z≥3 on timescales of up to 30–40 yrs. The study involved simultaneous RATAN-600 measurements at the frequencies of 2.3, 4.7, 8.2, 11.2, and 22.3 GHz in 2017–2020. In addition, data from the literature were used. We have found that the variability index, VS, which quantifies the normalized difference between the maximum and minimum flux density while accounting for measurement uncertainties, ranges from 0.02 to 0.96 for the quasars. Approximately half of the objects in the sample exhibit a variability index within the range from 0.25 to 0.50, which is comparable to that observed in blazars at lower redshifts. The distribution of VS at 22.3 GHz is significantly different from that at 2.3–11.2 GHz, which may be attributed to the fact that a compact AGN core dominates at the source’s rest frame frequencies greater than 45 GHz, leading to higher variability indices obtained at 22.3 GHz (the VS distribution peaks around 0.4) compared to the lower frequencies (the VS distribution at 2.3 and 4.7 GHz peaks around 0.1–0.2). Several source groups with distinctive variability characteristics were found using the cluster analysis of quasars. We propose seven new candidates for gigahertz-peaked spectrum (GPS) sources and five new megahertz-peaked spectrum (MPS) sources based on their spectrum shape and variability features. Only 6 out of the 23 sources previously reported as GPS demonstrate a low variability level typical of classical GPS sources (VS&lt;0.25) at 4.7–22.3 GHz. When excluding the highly variable peaked-spectrum blazars, we expect no more than 20% of the sources in the sample to be GPS candidates and no more than 10% to be MPS candidates.

The black hole masses of high-redshift QSOs

ABSTRACT Observations of high-redshift quasars frequently promote suggestions of large black hole masses, whose presence so early in cosmic time is not easily explicable. I consider the parallel with ultraluminous X-ray sources (ULXs) – now known to be stellar-mass black hole (and neutron star) binaries apparently radiating far above their Eddington luminosities LEdd. The true luminosity in ULXs is actually only of order LEdd, for stellar-mass accretors, but has a very anisotropic (‘beamed’) component, plus a near-isotropic component of similar luminosity but much lower specific intensity. Observers viewing ULXs from within the beam but assuming spherical symmetry deduce a luminosity ≫LEdd. These features appear because the accretors are fed mass at highly super-Eddington rates, most of it expelled in high-speed (v &amp;gt; 0.2c) outflows from the accretion disc. I show that in similarly beamed AGN, emission-line properties would be essentially the same as in unbeamed sources, but standard virial mass indicators unusable because velocity widths are dominated by the outflows, not bound motions about the black holes. In an ensemble of this kind the apparently most luminous systems are always the most distant, but have the lowest black hole masses. Interpreting observations of this ensemble without knowing that they are beamed leads instead to very high black hole mass estimates. The analogy with ULXs therefore suggests that high-redshift quasars might actually have central black hole masses which could have grown from stellar values within the lookback time. I consider how one might test these ideas observationally.

Primordial origin of supermassive black holes from axion bubbles

We study a modification of the primordial black hole (PBH) formation model from axion bubbles.We assume that the Peccei-Quinn scalar rolls down in the radial direction from a large field value to the potential minimum during inflation, which suppresses the axion fluctuations and weakens the clustering of PBHs on large scales.We find that the modified model can produce a sufficient number of PBHs that seed the supermassive black holes while avoiding the observational constraints from isocurvature perturbations and angular correlation of the high-redshift quasars.

EIGER. V. Characterizing the Host Galaxies of Luminous Quasars at z ≳ 6

Abstract We report JWST/NIRCam measurements of quasar host galaxy emissions and supermassive black hole (SMBH) masses for six quasars at 5.9 &lt; z &lt; 7.1 in the Emission-line galaxies and Intergalactic Gas in the Epoch of Reionization (EIGER) project. We obtain deep NIRCam imaging in the F115W, F200W, and F356W bands, as well as F356W grism spectroscopy of the quasars. We use bright unsaturated stars to construct models of the point-spread functions (PSFs) and estimate the errors of these PSFs. We then measure or constrain the fluxes and morphology of the quasar host galaxies by fitting the quasar images as a point source plus an exponential disk. We successfully detect the host galaxies of three quasars, which have host-to-quasar-flux ratios of ∼1%–5%. Spectral energy distribution fitting suggests that these quasar host galaxies have stellar masses of M * ≳ 1010 M ⊙. For quasars with host galaxy nondetections, we estimate the upper limits of their stellar masses. We use the grism spectra to measure the Hβ line profile and the continuum luminosity, then estimate the SMBH masses for the quasars. Our results indicate that the positive relation between SMBH masses and host galaxy stellar masses already exists at redshift z ≳ 6. The quasars in our sample show a high BH-to-stellar-mass ratio of M BH/M * ∼ 0.15, which is about ∼2 dex higher than local relations. We find that selection effects only contribute partially to the high M BH/M * ratios of high-redshift quasars. This result hints at a possible redshift evolution of the M BH–M * relation.

Constraining on the non-standard cosmological models combining the observations of high-redshift quasars and BAO

In this work, we studied four types of cosmological models with different mechanisms driving the accelerated expansion of the universe, include Braneworld models, Chaplygin Gas models, Emergent Dark Energy models, and cosmological torsion models. Considering that the dynamics of these models at low redshifts are very similar and difficult to distinguish, we used the latest and largest UV and X-ray measurements of quasars (QSOs) observations covering the range of redshift 0.009<z<7.5\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$0.009<z<7.5$$\\end{document}. However, the high intrinsic dispersion of this sample and the degeneracy between cosmological model parameters, we added 2D-BAO and 3D-BAO datasets to help us constrain the parameters of these cosmological models. Our results suggest that standard cold dark matter scenario may not be the best cosmological model preferred by the high-redshift observations. The Generalized Chaplygin Gas (GCG) and cosmological constant plus torsion (named Case II) models perform best by Akaike Information Criterion (AIC), but the Λ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Lambda $$\\end{document}CDM is the best cosmological model preferred by Bayesian Information Criterion (BIC). Our work also supports that the Phenomenologically Emergent Dark Energy and cosmological torsion models may alleviate the Hubble tension, the reported value of the Hubble constant obtained from QSO+BAO datasets combination lies between Planck 2018 observations and local measurements from the SH0ES collaboration, while other cosmological models all support that the Hubble constant tends to be closer to recent Planck 2018 results, but these model are penalized by information criterion.

The CluMPR galaxy cluster-finding algorithm and DESI legacy survey galaxy cluster catalogue

ABSTRACT Galaxy clusters enable unique opportunities to study cosmology, dark matter, galaxy evolution, and strongly lensed transients. We here present a new cluster-finding algorithm, CluMPR (Clusters from Masses and Photometric Redshifts), that exploits photometric redshifts (photo-z’s) as well as photometric stellar mass measurements. CluMPR uses a 2D binary search tree to search for overdensities of massive galaxies with similar redshifts on the sky and then probabilistically assigns cluster membership by accounting for photo-z uncertainties. We leverage the deep DESI Legacy Survey grzW1W2 imaging over one-third of the sky to create a catalogue of $\sim 300\, 000$ galaxy cluster candidates out to z = 1, including tabulations of member galaxies and estimates of each cluster’s total stellar mass. Compared to other methods, CluMPR is particularly effective at identifying clusters at the high end of the redshift range considered (z = 0.75–1), with minimal contamination from low-mass groups. These characteristics make it ideal for identifying strongly lensed high-redshift supernovae and quasars that are powerful probes of cosmology, dark matter, and stellar astrophysics. As an example application of this cluster catalogue, we present a catalogue of candidate wide-angle strongly lensed quasars in Appendix C. The nine best candidates identified from this sample include two known lensed quasar systems and a possible changing-look lensed QSO with SDSS spectroscopy. All code and catalogues produced in this work are publicly available (see Data Availability).

A Candidate Supermassive Black Hole in a Gravitationally Lensed Galaxy at Z ≈ 10

While supermassive black holes (SMBHs) are widely observed in the nearby and distant Universe, their origin remains debated with two viable formation scenarios with light and heavy seeds. In the light seeding model, the seed of the first SMBHs form from the collapse of massive stars with masses of 10–100 M ⊙, while the heavy seeding model posits the formation of 104–5 M ⊙ seeds from direct collapse. The detection of SMBHs at redshifts z ≳ 10, edging closer to their formation epoch, provides critical observational discrimination between these scenarios. Here, we focus on the JWST-detected galaxy, GHZ 9, at z ≈ 10 that is lensed by the foreground cluster, A2744. Based on 2.1 Ms deep Chandra observations, we detect a candidate X-ray active galactic nucleus (AGN), which is spatially coincident with the high-redshift galaxy, GHZ 9. The SMBH candidate is inferred to have a bolometric luminosity of (1.0−0.4+0.5)×1046ergs−1 , which corresponds to a black hole (BH) mass of (8.0−3.2+3.7)×107M⊙ assuming Eddington-limited accretion. This extreme mass at such an early cosmic epoch suggests the heavy seed origin for this BH candidate. Based on the Chandra and JWST discoveries of extremely high-redshift quasars, we have constructed the first simple AGN luminosity function extending to z ≈ 10. Comparison of this luminosity function with theoretical models indicates an overabundant z ≈ 10 SMBH population, consistent with a higher-than-expected seed formation efficiency.

Quasar Island – three new z ∼ 6 quasars, including a lensed candidate, identified with contrastive learning

ABSTRACT Of the hundreds of z ≳ 6 quasars discovered to date, only one is known to be gravitationally lensed, despite the high lensing optical depth expected at z ≳ 6. High-redshift quasars are typically identified in large-scale surveys by applying strict photometric selection criteria, in particular by imposing non-detections in bands blueward of the Lyman-α line. Such procedures by design prohibit the discovery of lensed quasars, as the lensing foreground galaxy would contaminate the photometry of the quasar. We present a novel quasar selection methodology, applying contrastive learning (an unsupervised machine learning technique) to Dark Energy Survey imaging data. We describe the use of this technique to train a neural network which isolates an ‘island’ of 11 sources, of which seven are known z ∼ 6 quasars. Of the remaining four, three are newly discovered quasars (J0109−5424, z = 6.07; J0122−4609, z = 5.99; J0603−3923, z = 5.94), as confirmed by follow-up and archival spectroscopy, implying a 91 per cent efficiency for our novel selection method; the final object on the island is a brown dwarf. In one case (J0109−5424), emission below the Lyman limit unambiguously indicates the presence of a foreground source, though high-resolution optical/near-infrared imaging is still needed to confirm the quasar’s lensed (multiply imaged) nature. Detection in the g band has led this quasar to escape selection by traditional colour cuts. Our findings demonstrate that machine learning techniques can thus play a key role in unveiling populations of quasars missed by traditional methods.