Quasar Population Research Articles

Discovery of Six Low-z, High-luminosity, and High-mass Quasars

Abstract This study presents the discovery of six low-redshift quasars, identified through spectral observations conducted with the RTT-150 telescope. These quasars, with redshifts ranging from 0.3 to 0.6, were selected incorporating their placement in the color–color diagrams and detection by ROSAT 2RXS, focusing on those resembling the spectral energy distribution of a quasar. Our analysis includes detailed modeling of their continuum and emission line properties, which allowed us to estimate their bolometric luminosities, central black hole masses, and Eddington ratios. The findings reveal that these quasars exhibit exceptionally high luminosities and high masses for their redshift values, making them rare examples in the low-redshift quasar population. The results contribute to our understanding of quasar characteristics and their role in probing the structure and evolution of the nearby universe. This research underscores the importance of continued exploration of low-redshift quasars to enhance our knowledge of cosmic evolution and the dynamics of supermassive black holes in the early universe.

A unified model for the clustering of quasars and galaxies at z ≈ 6

ABSTRACT Recent observations from the EIGER JWST program have measured for the first time the quasar–galaxy cross-correlation function at $z\approx 6$. The autocorrelation function of faint $z\approx 6$ quasars was also recently estimated. These measurements provide key insights into the properties of quasars and galaxies at high redshift and their relation with the host dark matter haloes. In this work, we interpret these data building upon an empirical quasar population model that has been applied successfully to quasar clustering and demographic measurements at $z\approx 2\!-\!4$. We use a new, large-volume N-body simulation with more than a trillion particles, FLAMINGO-10k, to model quasars and galaxies simultaneously. We successfully reproduce observations of $z\approx 6$ quasars and galaxies (i.e. their clustering properties and luminosity functions), and infer key quantities such as their luminosity–halo mass relation, the mass function of their host haloes, and their duty cycle/occupation fraction. Our key findings are (i) quasars reside on average in $\approx 10^{12.5}\, {\rm M}_{\odot }$ haloes (corresponding to $\approx 5\sigma$ fluctuations in the initial conditions of the linear density field), but the distribution of host halo masses is quite broad; (ii) the duty cycle of (UV-bright) quasar activity is relatively low ($\approx 1~{{\ \rm per\ cent}}$); (iii) galaxies (that are bright in [O iii]) live in much smaller haloes ($\approx 10^{10.9}\, {\rm M}_{\odot }$) and have a larger duty cycle (occupation fraction) of $\approx 13~{{\ \rm per\ cent}}$. Finally, we focus on the inferred properties of quasars and present a homogeneous analysis of their evolution with redshift. The picture that emerges reveals a strong evolution of the host halo mass and duty cycle of quasars at $z\approx 2\!-\!6$, and calls for new investigations of the role of quasar activity across cosmic time.

Red Quasars: Estimation of SMBH Spin, Mass, and Accretion Disk Inclination Angle

ABSTRACTWe estimated values of spin, mass, and inclination angle for sample of 42 red quasars. Our estimations show that two objects: F2MS J1113+1244 and F2MS J1434+0935 with the highest Eddington ratios may have geometrically thick disk. Six objects: SDSS J0036‐0113, S82X 0040+0058, S82X 0118+0018, S82X 0303‐0115, FBQS J1227+3214, S82X 2328‐0028 may have “retrograde” rotation. Analysis of estimated spin values shows that red quasar population may contain Seyfert galaxies and NLS1.

The SRG/eROSITA All-Sky Survey: X-ray beacons at late cosmic dawn

The Spektrum-Roentgen-Gamma (SRG)/extended Roentgen Survey with an Imaging Telescope Array (eROSITA) All-Sky Survey (eRASS) is expected to contain $ quasars that emitted their light when the universe was less than a billion years old, that is, at $z>5.6$. By selection, these quasars populate the bright end of the active galactic nuclei (AGN) X-ray luminosity function, and their space density offers a powerful demographic diagnostic of the parent super-massive black hole (SMBH) population. Of the $ 400$ quasars that have been discovered at $z>5.6$ to date, less than 15<!PCT!> have been X-ray detected. We present a pilot survey to uncover the elusive X-ray luminous end of the distant quasar population. We have designed a quasar selection pipeline based on optical, infrared and X-ray imaging data from DES DR2, VHS DR5, CatWISE2020 and the eRASS (up to its four-pass cumulative version, eRASS:4). The core selection method relies on SED template fitting. We performed optical follow-up spectroscopy with the Magellan /LDSS3 instrument for the redshift confirmation of a subset of candidates. We have further obtained a deeper X-ray image of one of our candidates with Chandra ACIS-S. We report the discovery of five new quasars in the redshift range $5.6 < z < 6.1$. Two of these quasars are detected in eRASS and are, therefore, X-ray ultra-luminous by selection. We also report the detection of these quasars at radio frequencies. The first one is a broad absorption line quasar, which shows significant, order-of-magnitude X-ray dimming over 3.5 years, corresponding to six months in the quasar rest frame. The second X-ray detected quasar is a jetted source with compact morphology. We show that a blazar configuration is likely for this source, making it one of the most distant blazars known to date. With our pilot study, we demonstrate the power of eROSITA as a discovery machine for luminous quasars in the epoch of reionization. The X-ray emission of the two eROSITA detected quasars are likely to be driven by different high-energetic emission mechanisms, a diversity which we will further explore in a future systematic full-hemisphere survey.

The NuSTAR Serendipitous Survey: The 80 Month Catalog and Source Properties of the High-energy Emitting Active Galactic Nucleus and Quasar Population

We present a catalog of hard X-ray serendipitous sources detected in the first 80 months of observations by the Nuclear Spectroscopic Telescope Array (NuSTAR). The NuSTAR serendipitous survey 80 month (NSS80) catalog has an unprecedented ∼62 Ms of effective exposure time over 894 unique fields (a factor of 3 increase over the 40 month catalog, NSS40), with an areal coverage of ∼36 deg2, larger than all NuSTAR extragalactic surveys. NSS80 provides 1274 hard X-ray sources in the 3−24 keV band (822 new detections compared to the previous NSS40). Approximately 76% of the NuSTAR sources have lower-energy (<10 keV) X-ray counterparts from Chandra, XMM-Newton, and Swift-XRT. We have undertaken an extensive campaign of ground-based spectroscopic follow-up to obtain new source redshifts and classifications for 427 sources. Combining these with existing archival spectroscopy provides redshifts for 550 NSS80 sources, of which 547 are classified. The sample is primarily composed of active galactic nuclei (AGNs), detected over a large range in redshift (z = 0.012–3.43), but also includes 58 spectroscopically confirmed Galactic sources. In addition, five AGN–galaxy pairs, one dual AGN system, one BL Lac candidate, and a hotspot of 4C 74.26 (radio quasar) have been identified. The median rest-frame 10−40 keV luminosity and redshift of NSS80 are 〈L 10−40 keV〉 = 1.2 × 1044 erg s−1 and 〈z〉 = 0.56. We investigate the optical properties and construct composite optical spectra to search for subtle signatures not present in the individual spectra, finding an excess of redder BL AGNs compared to optical quasar surveys, predominantly due to the presence of the host galaxy and, at least in part, due to dust obscuration.

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

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

O iii] λ5007 emissions in extremely red quasars (ERQs) are compact

ABSTRACT ‘Extremely red quasars’ (ERQs) are a non-radio-selected, intrinsically luminous population of quasars at cosmic noon selected by their extremely red colour from rest-frame UV to mid-IR. ERQs are uniquely associated with exceptionally broad and blueshifted [O iii] $\lambda$5007 emission reaching speeds $\gt $6000 km s$^{-1}$. We obtained adaptive optics integral-field spectroscopic observations using Keck/OSIRIS and Gemini/NIFS of a sample of 10 ERQs with bolometric luminosities (10$^{47.0}$–10$^{47.9}$) erg s$^{-1}$ at $z\sim$ (2.3–3.0). The goal is to measure the sizes and spatially resolved kinematics of the [O iii]-emitting regions. We study the surface brightness maps and aperture-extracted spectra and model the point-spread functions. We identify signs of merger activities in the continuum emissions. We identify physically distinct [O iii] kinematic components that are bimodal and respectively trace ERQ-driven outflows of velocity dispersion $\gtrsim$250 km s$^{-1}$ and dynamically quiescent interstellar media. We find that the ERQ-driven ionized outflows are typically at $\sim$1 kpc scales whereas the quiescent ionized gas extend to a few kpc. Compared to normal quasars the extremely fast ERQ-driven [O iii] outflows tend to be more compact, supporting the notion that ERQs are in a young stage of quasar/galaxy evolution and represent unique physical conditions beyond orientation differences with normal quasar populations. The kinematically quiescent [O iii] emissions in ERQs tend to be spatially resolved but less extended than in normal quasars, which can be explained by global and patchy dust obscuration. The hint of ionization cones suggests some of the obscuration can be partially explained by a patchy torus.

A quasar-galaxy merger at z ∼ 6.2: Black hole mass and quasar properties from the NIRSpec spectrum

We present JWST/NIRSpec integral field data of the quasar PJ308-21 at z = 6.2342. As shown by previous ALMA and HST imaging, the quasar has two companion sources, interacting with the quasar host galaxy. The high-resolution G395H/290LP NIRSpec spectrum covers the 2.87 − 5.27 μm wavelength range and shows the rest-frame optical emission of the quasar with exquisite quality (signal-to-noise ratio ∼100 − 400 per spectral element). Based on the Hβ line from the broad line region, we obtain an estimate of the black hole mass MBH, Hβ ∼ 2.7 × 109 M⊙. This value is within a factor ≲1.5 of the Hα-based black hole mass from the same spectrum (MBH, Hα ∼ 1.93 × 109 M⊙) and is consistent with a previous estimate relying on the Mg IIλ2799 line (MBH, MgII ∼ 2.65 × 109 M⊙). All these MBH estimates are within the ∼0.5 dex intrinsic scatter of the adopted mass calibrations. The high Eddington ratio of PJ308-21 λEdd, Hβ ∼ 0.67 (λEdd, Hα ∼ 0.96) is in line with the overall quasar population at z ≳ 6. The relative strengths of the [O III], Fe II, and Hβ lines are consistent with the empirical “Eigenvector 1” correlations as observed for low redshift quasars. We find evidence for blueshifted [O III] λ5007 emission with a velocity offset Δv[O III] = −1922 ± 39 km s−1 from the systemic velocity and a full width at half maximum (FWHM) FWHM([O III]) = 2776−74+75 km s−1. This may be the signature of outflowing gas from the nuclear region, despite the true values of Δv[O III] and FWHM([O III]) likely being more uncertain due to the blending with Hβ and Fe II lines. Our study demonstrates the unique capabilities of NIRSpec in capturing quasar spectra at cosmic dawn and studying their properties in unprecedented detail.

CatNorth: An Improved Gaia DR3 Quasar Candidate Catalog with Pan-STARRS1 and CatWISE

A complete and pure sample of quasars with accurate redshifts is crucial for quasar studies and cosmology. In this paper, we present CatNorth, an improved Gaia Data Release 3 (Gaia DR3) quasar candidate catalog with more than 1.5 million sources in the 3π sky built with data from Gaia, Pan-STARRS1, and CatWISE2020. The XGBoost algorithm is used to reclassify the original Gaia DR3 quasar candidates as stars, galaxies, and quasars. To construct training/validation data sets for the classification, we carefully built two different master stellar samples in addition to the spectroscopic galaxy and quasar samples. An ensemble classification model is obtained by averaging two XGBoost classifiers trained with different master stellar samples. Using a probability threshold of p QSO_mean > 0.95 in our ensemble classification model and an additional cut on the logarithmic probability density of zero proper motion, we retrieved 1,545,514 reliable quasar candidates from the parent Gaia DR3 quasar candidate catalog. We provide photometric redshifts for all candidates with an ensemble regression model. For a subset of 89,100 candidates, accurate spectroscopic redshifts are estimated with the convolutional neural network from the Gaia BP/RP spectra. The CatNorth catalog has a high purity of ∼90%, while maintaining high completeness, which is an ideal sample to understand the quasar population and its statistical properties. The CatNorth catalog is used as the main source of input catalog for the Large Sky Area Multi-Object Fiber Spectroscopic Telescope phase III quasar survey, which is expected to build a highly complete sample of bright quasars with i < 19.5.

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.

Absence of radio-bright dominance in a near-infrared selected sample of red quasars

Context. The dichotomy between red and blue quasars is still an open question. It is debated whether red quasars are simply blue quasars that are observed at certain inclination angles or if they provide insight into a transitional phase in the evolution of quasars. Aims. We investigate the relation between quasar colors and radio-detected fraction because radio observations of quasars provide a powerful tool in distinguishing between quasar models. Methods. We present the eHAQ+GAIA23 sample, which contains quasars from the High A(V) Quasar (HAQ) Survey, the Extended High A(V) Quasar (eHAQ) Survey, and the Gaia quasar survey. All quasars in this sample have been found using a near-infrared color selection of target candidates that have otherwise been missed by the Sloan Digital Sky Survey (SDSS). We implemented a redshift-dependent color cut in g* − i* to select red quasars in the sample and divided them into redshift bins, while using a nearest-neighbors algorithm to control for luminosity and redshift differences between our red quasar sample and a selected blue sample from the SDSS. Within each bin, we cross-matched the quasars to the Faint Images of the Radio Sky at Twenty centimeters (FIRST) survey and determined the radio-detection fraction. Results. For redshifts 0.8 &lt; z ≤ 1.5, the red and blue quasars have a radio-detection fraction of 0.153−0.032+0.037 and 0.132−0.030+0.034, respectively. The red and blue quasars with redshifts 1.5 &lt; z ≤ 2.4 have radio-detection fractions of 0.059−0.016+0.019 and 0.060−0.016+0.019, respectively, and the red and blue quasars with redshifts z &gt; 2.4 have radio-detection fractions of 0.029−0.012+0.017 and 0.058−0.019+0.024, respectively. For the WISE color-selected red quasars, we find a radio-detection fraction of 0.160−0.034+0.038 for redshifts 0.8 &lt; z ≤ 1.5, 0.063−0.017+0.020 for redshifts 1.5 &lt; z ≤ 2.4, and 0.051−0.022+0.030 for redshifts z &gt; 2.4. In other words, we find similar radio-detection fractions for red and blue quasars within &lt; 1σ uncertainty, independent of redshift. This disagrees with what has been found in the literature for red quasars in SDSS. It should be noted that the fraction of broad absorption line (BAL) quasars in red SDSS quasars is about five times lower. BAL quasars have been observed to be more frequently radio quiet than other quasars, therefore the difference in BAL fractions could explain the difference in radio-detection fraction. Conclusions. The dusty torus of a quasar is transparent to radio emission. When we do not observe a difference between red and blue quasars, it leads us to argue that orientation is the main cause of quasar redness. Moreover, the observed higher proportion of BAL quasars in our dataset relative to the SDSS sample, along with the higher rate of radio detections, indicates an association of the redness of quasars and the inherent BAL fraction within the overall quasar population. This correlation suggests that the redness of quasars is intertwined with the inherent occurrence of BAL quasars within the entire population of quasars. In other words, the question why some quasars appear red or exhibit BAL characteristics might not be isolated; it could be directly related to the overall prevalence of BAL quasars in the quasar population. This finding highlights the need to explore the underlying factors contributing to both the redness and the frequency of BAL quasars, as they appear to be interconnected phenomena.

How does the radio enhancement of broad absorption line quasars relate to colour and accretion rate?

ABSTRACT The origin of radio emission in different populations of radio-quiet quasars is relatively unknown, but recent work has uncovered various drivers of increased radio-detection fraction. In this work, we pull together three known factors: optical colour (g − i), C iv distance (a proxy for L/LEdd), and whether or not the quasar contains broad absorption lines (BALQSOs) which signify an outflow. We use SDSS (Sloan Digital Sky Survey) DR14 spectra along with the LOFAR Two Metre Sky Survey Data Release 2 and find that each of these properties have an independent effect. BALQSOs are marginally more likely to be radio-detected than non-BALQSOs at similar colours and L/LEdd, moderate reddening significantly increases the radio-detection fraction and the radio detection increases with L/LEdd above a threshold for all populations. We test a widely used simple model for radio wind shock emission and calculate energetic efficiencies that would be required to reproduce the observed radio properties. We discuss interpretations of these results concerning radio-quiet quasars more generally. We suggest that radio emission in BALQSOs is connected to a different physical origin than the general quasar population since they show different radio properties independent of colour and C iv distance.

Extreme Variability Quasars in Their Various States. II. Spectral Variation Revealed with Multiepoch Spectra

We previously built a sample of 14,012 extremely variable quasars (EVQs) based on Sloan Digital Sky Survey (SDSS) and Pan-STARRS1 photometric observations. In this work we present the spectral fitting to their SDSS spectra and study the spectral variation in 1259 EVQs with multiepoch SDSS spectra (after prudently excluding spectra with potentially unreliable spectroscopic photometry). We find a clear “bluer-when-brighter” trend in EVQs, consistent with previous findings of normal quasars and active galactic nuclei. We detect significant intrinsic Baldwin effect (iBeff, i.e., smaller line equivalent width at higher continuum flux in individual active galactic nuclei) in the broad Mg ii and C iv lines of EVQs. Meanwhile, no systematical iBeff is found for the broad Hβ line, which could be attributed to strong host contamination at longer wavelengths. Remarkably, by comparing the iBeff slope of EVQs with archived changing-look quasars, we show that the changing-look quasars identified in the literature are most likely a biased (due to its definition) subpopulation of EVQs, rather than a distinct population of quasars. We also found no significant broad line breathing of Hβ, Mg ii, or C iv, suggesting the broad line breathing in quasars may disappear at longer timescales (∼3000 days).

The radio-loud fraction of quasars at z &amp;gt; 6

ABSTRACT Quasars at redshifts z &amp;gt; 6 are an excellent probe of the formation and evolution of supermassive black holes in the early Universe. The population of radio-luminous quasars is of particular interest, as such quasars could potentially be used to study the neutral intergalactic medium during cosmic reionization via H i 21 cm absorption studies. However, the lack of deep radio observations of z &amp;gt; 6 quasars leaves the population poorly constrained, and suitable candidates for an H i 21 cm absorption study have yet to be found. In this work, we present Jansky Very Large Array (VLA) 1–2 GHz radio continuum observations of 138 quasars at redshifts 6.0 ≤ z &amp;lt; 7.6. We detect the radio continuum emission of the z = 6.1 quasar J1034−1425, with a 1.6 GHz flux density of $170\pm 36\, \mu$Jy. This quasar is radio-quiet with radio-loudness, R ≡ f5 GHz/fν, 4400 Å = 2.4 ± 0.5. In addition, we detect seven other quasars at z &amp;gt; 6, which have previously been characterized in the literature at these frequencies. Using the full sample, we estimate the radio-loud fraction to be $3.8^{+6.2}_{-2.4}\ \hbox{per cent}$, where the uncertainties are 95 per cent confidence intervals. This is lower than recent estimates of the radio-loud fraction in the literature, but is still marginally consistent with no redshift evolution of the radio-loud fraction. We explore the undetected quasar population by stacking their continuum images at their optical positions and obtain a median stacked flux density of 13.8 ± 3.9 µJy and luminosity of log L5 GHz/(W Hz−1) = 24.2 ± 0.1.

A benchmark for extreme conditions of the multiphase interstellar medium in the most luminous hot dust-obscured galaxy at z = 4.6

WISE J224607.6–052634.9 (W2246–0526) is a hot dust-obscured galaxy at z = 4.601, and the most luminous obscured quasar known to date. W2246–0526 harbors a heavily obscured supermassive black hole that is most likely accreting above the Eddington limit. We present observations with the Atacama Large Millimeter/submillimeter Array (ALMA) in seven bands, including band 10, of the brightest far-infrared (FIR) fine-structure emission lines of this galaxy: [OI]63 μm, [OIII]88 μm, [NII]122 μm, [OI]145 μm, [CII]158 μm, [NII]205 μm, [CI]370 μm, and [CI]609 μm. A comparison of the data to a large grid of CLOUDY radiative transfer models reveals that a high hydrogen density (nH ∼ 3 × 103 cm−3) and extinction (AV ∼ 300 mag), together with extreme ionization (log(U) = − 0.5) and a high X-ray to UV ratio (αox ≥ −0.8) are required to reproduce the observed nuclear line ratios. The values of αox and U are among the largest found in the literature and imply the existence of an X-ray-dominated region (XDR). In fact, this component explains the a priori very surprising non-detection of the [OIII]88 μm emission line, which is actually suppressed, instead of boosted, in XDR environments. Interestingly, the best-fitted model implies higher X-ray emission and lower CO content than what is detected observationally, suggesting the presence of a molecular gas component that should be further obscuring the X-ray emission over larger spatial scales than the central region that is being modeled. These results highlight the need for multiline infrared observations to characterize the multiphase gas in high redshift quasars and, in particular, W2246–0526 serves as an extreme benchmark for comparisons of interstellar medium conditions with other quasar populations at cosmic noon and beyond.

Accretion and obscuration in merger-dominated luminous red quasars

ABSTRACT We present an analysis of the X-ray properties of 10 luminous, dust-reddened quasars from the FIRST-2MASS (F2M) survey based on new and archival Chandra observations. These systems are interpreted to be young, transitional objects predicted by merger-driven models of quasar/galaxy co-evolution. The sources have been well studied from the optical through mid-infrared, have Eddington ratios above 0.1, and possess high-resolution imaging, most of which shows disturbed morphologies indicative of a recent or ongoing merger. When combined with previous X-ray studies of five other F2M red quasars, we find that the sources, especially those hosted by mergers, have moderate-to-high column densities (NH ≃ 1022.5–23.5 cm−2) and Eddington ratios high enough to enable radiation pressure to blow out the obscuring material. We confirm previous findings that red quasars have dust-to-gas ratios that are significantly lower than the value for the Milky Way’s interstellar medium, especially when hosted by a merger. The dust-to-gas ratio for two red quasars that lack evidence for merging morphology is consistent with the Milky Way and they do not meet the radiative feedback conditions for blow-out. These findings support the picture of quasar/galaxy co-evolution in which a merger results in feeding of and feedback from an AGN. We compare the F2M red quasars to other obscured and reddened quasar populations in the literature, finding that, although morphological information is lacking, nearly all such samples meet blow-out conditions and exhibit outflow signatures suggestive of winds and feedback.

The intrinsic X-ray luminosity distribution of an optically selected SDSS quasar population

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

Constraining the Properties of Black Hole Seeds from the Farthest Quasars

Over 60 yr after the discovery of the first quasar, more than 275 such sources are identified in the epoch of reionization at z > 6. JWST is now exploring higher redshifts (z ≳ 8) and lower-mass (≲107 M ⊙) ranges. The discovery of progressively farther quasars is instrumental to constraining the properties of the first population of black holes (BHs), or BH seeds, formed at z ∼ 20–30. For the first time, we use Bayesian analysis of the most comprehensive catalog of quasars at z > 6 to constrain the distribution of BH seeds. We show that the mass distribution of BH seeds can be effectively described by combining a power law and a lognormal function tailored to the mass ranges associated with light and heavy seeds, assuming Eddington-limited growth and early seeding time. Our analysis reveals a power-law slope of and a lognormal mean of . The inferred values of the Eddington ratio, the duty cycle, and the mean radiative efficiency are , , and , respectively. Models that solely incorporate a power law or a lognormal distribution within the specific mass range corresponding to light and heavy seeds are statistically strongly disfavored, unlike models not restricted to this specific range. Our results suggest that including both components is necessary to comprehensively account for the masses of high-redshift quasars, and that both light and heavy seeds formed in the early Universe and grew to form the population of quasars we observe.

Constraining the AGN Luminosity Function from JWST with the X-Ray Background

We predict the X-ray background (XRB) expected from the population of quasars detected by the James Webb Space Telescope spectroscopic surveys over the redshift range z ∼ 4–7. We find that the measured UV emissivities imply a ∼10 times higher unresolved XRB than constrained by current experiments. We illustrate the difficulty of simultaneously matching the faint end of the quasar luminosity function and the XRB constraints. We discuss possible origins and consequences of this discrepancy.

Quasar Factor Analysis—An Unsupervised and Probabilistic Quasar Continuum Prediction Algorithm with Latent Factor Analysis

Since their first discovery, quasars have been essential probes of the distant Universe. However, due to our limited knowledge of its nature, predicting the intrinsic quasar continua has bottlenecked their usage. Existing methods of quasar continuum recovery often rely on a limited number of high-quality quasar spectra, which might not capture the full diversity of the quasar population. In this study, we propose an unsupervised probabilistic model, quasar factor analysis (QFA), which combines factor analysis with physical priors of the intergalactic medium to overcome these limitations. QFA captures the posterior distribution of quasar continua through generatively modeling quasar spectra. We demonstrate that QFA can achieve the state-of-the-art performance, ∼2% relative error, for continuum prediction in the Lyα forest region compared to previous methods. We further fit 90,678 2 < z < 3.5, signal-to-noise ratio >2 quasar spectra from Sloan Digital Sky Survey Data Release 16 and found that for ∼30% quasar spectra where the continua were ill-determined with previous methods, QFA yields visually more plausible continua. QFA also attains ≲1% error in the 1D Lyα power spectrum measurements at z ∼ 3 and ∼4% in z ∼ 2.4. In addition, QFA determines latent factors representing more physical motivation than principal component analysis. We investigate the evolution of the latent factors and report no significant redshift or luminosity dependency except for the Baldwin effect. The generative nature of QFA also enables outlier detection robustly; we showed that QFA is effective in selecting outlying quasar spectra, including damped Lyα systems and potential Type II quasar spectra.