Distant quasars are unique tracers to study the formation of the earliest supermassive black holes (SMBHs) and the history of cosmic reionization. Despite extensive efforts, only two quasars have been found at z ≥ 7.5, due to a combination of their low spatial density and the high contamination rate in quasar selection. We report the discovery of a luminous quasar at z = 7.642, J0313−1806, the most distant quasar yet known. This quasar has a bolometric luminosity of 3.6 × 10¹³ L _⊙. Deep spectroscopic observations reveal a SMBH with a mass of (1.6 ± 0.4) × 10⁹ M _⊙ in this quasar. The existence of such a massive SMBH just ∼670 million years after the big bang challenges significantly theoretical models of SMBH growth. In addition, the quasar spectrum exhibits strong broad absorption line (BAL) features in C iv and Si iv, with a maximum velocity close to 20% of the speed of light. The relativistic BAL features, combined with a strongly blueshifted C iv emission line, indicate that there is a strong active galactic nucleus (AGN)-driven outflow in this system. Atacama Large Millimeter/submillimeter Array observations detect the dust continuum and [C ii] emission from the quasar host galaxy, yielding an accurate redshift of 7.6423 ± 0.0013 and suggesting that the quasar is hosted by an intensely star-forming galaxy, with a star formation rate of ∼200 M _⊙ yr⁻¹ and a dust mass of ∼7 × 10⁷ M _⊙. Follow-up observations of this reionization-era BAL quasar will provide a powerful probe of the effects of AGN feedback on the growth of the earliest massive galaxies.

We report the discovery of a luminous quasar, J1007+2115 at z = 7.515 (“Pōniuā‘ena”), from our wide-field reionization-era quasar survey. J1007+2115 is the second quasar now known at z > 7.5, deep into the reionization epoch. The quasar is powered by a (1.5 ± 0.2) × 10⁹ M _⊙ supermassive black hole (SMBH), based on its broad Mg ii emission-line profile from Gemini and Keck near-IR spectroscopy. The SMBH in J1007+2115 is twice as massive as that in quasar J1342+0928 at z = 7.54, the current quasar redshift record holder. The existence of such a massive SMBH just 700 million years after the Big Bang significantly challenges models of the earliest SMBH growth. Model assumptions of Eddington-limited accretion and a radiative efficiency of 0.1 require a seed black hole of ≳10⁴ M _⊙ at z = 30. This requirement suggests either a massive black hole seed as a result of direct collapse or earlier periods of rapid black hole growth with hyper-Eddington accretion and/or a low radiative efficiency. We measure the damping wing signature imprinted by neutral hydrogen absorption in the intergalactic medium (IGM) on J1007+2115's Lyα line profile, and find that it is weaker than that of J1342+0928 and two other z ≳ 7 quasars. We estimate an IGM volume-averaged neutral fraction . This range of values suggests a patchy reionization history toward different IGM sightlines. We detect the 158 μm [C ii] emission line in J1007+2115 with the Atacama Large Millimeter/submillimeter Array; this line centroid yields a systemic redshift of z = 7.5149 ± 0.0004 and indicates a star formation rate of ∼210 M _☉ yr⁻¹ in its host galaxy.

We report the results of near-infrared spectroscopic observations of 37 quasars in the redshift range 6.3 < z ≤ 7.64, including 32 quasars at z > 6.5, forming the largest quasar near-infrared spectral sample at this redshift. The spectra, taken with Keck, Gemini, VLT, and Magellan, allow investigations of central black hole mass and quasar rest-frame ultraviolet spectral properties. The black hole masses derived from the Mg ii emission lines are in the range (0.3–3.6) × 10⁹ M _⊙, which requires massive seed black holes with masses ≳10³–10⁴ M _⊙, assuming Eddington accretion since z = 30. The Eddington ratio distribution peaks at λ _Edd ∼ 0.8 and has a mean of 1.08, suggesting high accretion rates for these quasars. The C iv–Mg ii emission-line velocity differences in our sample show an increase of C iv blueshift toward higher redshift, but the evolutionary trend observed from this sample is weaker than the previous results from smaller samples at similar redshift. The Fe ii/Mg ii flux ratios derived for these quasars up to z = 7.6, compared with previous measurements at different redshifts, do not show any evidence of strong redshift evolution, suggesting metal-enriched environments in these quasars. Using this quasar sample, we create a quasar composite spectrum for z > 6.5 quasars and find no significant redshift evolution of quasar broad emission lines and continuum slope, except for a blueshift of the C iv line. Our sample yields a strong broad absorption line quasar fraction of ∼24%, higher than the fractions in lower-redshift quasar samples, although this could be affected by small sample statistics and selection effects.

This is the third paper in a series aimed at finding reionization-era quasars with the combination of DESI Legacy imaging Surveys (DELS), the Pan-STARRS1 (PS1) Survey, and near-infrared imaging surveys, such as the UKIRT Hemisphere Survey (UHS), as well as the Wide-field Infrared Survey Explorer (WISE) mid-infrared survey. In this paper, we describe the updated quasar candidate selection procedure, report the discovery of 16 quasars at 6.4 ≲ z ≲ 6.9 from an area of ∼13,020 deg², and present the quasar luminosity function (QLF) at z ∼ 6.7. The measured QLF follows in the magnitude range −27.6 < M ₁₄₅₀ < −25.5. We determine the quasar comoving spatial density at = 6.7 and M ₁₄₅₀ < −26.0 to be 0.39 ± 0.11 Gpc⁻³ and find the exponential density evolution parameter to be k = −0.78 ± 0.18 from z ∼ 6 to z ∼ 6.7, corresponding to a rapid decline by a factor of ∼6 per unit redshift toward earlier epochs. This indicates that the rapid decline of quasar spatial density at z > 5 that was found by previous works continues to z > 6, at a rate significantly faster than the average decline rate between z ∼ 3 and 5. We measured quasar comoving emissivity at z ∼ 6.7, which indicates that high-redshift quasars are highly unlikely to make a significant contribution to hydrogen reionization. The broad absorption line quasar fraction at z ≳ 6.5 is measured to be ≳22%. In addition, we also report the discovery of six additional quasars at z ∼ 6 in the Appendix.

We summarize and reanalyze observations bearing on missing galactic baryons, where we propose a consistent picture for halo gas in L ≳ L* galaxies. The hot X-ray-emitting halos are detected to 50–70 kpc, where typically M _hot(<50 kpc) ∼ 5 × 10⁹ M _⊙, and with density n ∝ r ^−3/2. When extrapolated to R ₂₀₀, the gas mass is comparable to the stellar mass, but about half of the baryons are still missing from the hot phase. If extrapolated to 1.7R ₂₀₀–3R ₂₀₀, the ratio of baryon to dark matter approaches the cosmic value. Significantly flatter density profiles are unlikely for R < 50 kpc, and they are disfavored but not ruled out for R > 50 kpc. For the Milky Way, the hot halo metallicity lies in the range 0.3–1 solar for R < 50 kpc. Planck measurements of the thermal Sunyaev–Zel’dovich (SZ) effect toward stacked luminous galaxies (primarily early type) indicate that most of their baryons are hot, are near the virial temperature, and extend beyond R ₂₀₀. This stacked SZ signal is nearly an order of magnitude larger than that inferred from the X-ray observations of individual (mostly spiral) galaxies with M _* > 10^11.3 M _⊙. This difference suggests that the hot halo properties are distinct for early- and late-type galaxies, possibly due to different evolutionary histories. For the cooler gas detected in UV absorption line studies, we argue that there are two absorption populations: extended halos, and disks extending to ∼50 kpc, containing most of this gas, and with masses a few times lower than the stellar masses. Such extended disks are also seen in 21 cm H i observations and in simulations.