In a multiwavelength survey of 13 quasars at 5.8 ≲ z ≲ 6.5, which were preselected to be potentially young, we find five objects with extremely small proximity zone sizes that may imply UV-luminous quasar lifetimes of ≲100,000 yr. Proximity zones are regions of enhanced transmitted flux in the vicinity of quasars that are sensitive to the quasars' lifetimes because the intergalactic gas has a finite response time to their radiation. We combine submillimeter observations from the Atacama Large Millimetre Array and the NOrthern Extended Millimeter Array, as well as deep optical and near-infrared spectra from the medium-resolution spectrograph on the Very Large Telescope and on the Keck telescopes, in order to identify and characterize these new young quasars, which provide valuable clues about the accretion behavior of supermassive black holes in the early universe and pose challenges on current black hole formation models to explain the rapid formation of billion-solar-mass black holes. We measure the quasars' systemic redshifts, black hole masses, Eddington ratios, emission-line luminosities, and star formation rates of their host galaxies. Combined with previous results, we estimate the fraction of young objects within the high-redshift quasar population at large to be 5% ≲ f young ≲ 10%. One of the young objects, PSO J158-14, shows a very bright dust continuum flux (F cont = 3.46 0.02 mJy), indicating a highly starbursting host galaxy with a star formation rate of approximately 1420 M o˙ yr-1.

During reionization, neutral hydrogen in the intergalactic medium (IGM) imprints a damping wing absorption feature on the spectrum of high-redshift quasars. A detection of this signature provides compelling evidence for a significantly neutral universe, and enables measurements of the hydrogen neutral fraction x H i(z) at that epoch. Obtaining reliable quantitative constraints from this technique, however, is challenging due to stochasticity induced by the patchy inside-out topology of reionization, degeneracies with quasar lifetime, and the unknown unabsorbed quasar spectrum close to rest-frame Lyα. We combine a large-volume semi-numerical simulation of reionization topology with 1D radiative transfer through high-resolution hydrodynamical simulations of the high-redshift universe to construct models of quasar transmission spectra during reionization. Our state-of-the-art approach captures the distribution of damping wing strengths in biased quasar halos that should have reionized earlier, as well as the erosion of neutral gas in the quasar environment caused by its own ionizing radiation. Combining this detailed model with our new technique for predicting the quasar continuum and its associated uncertainty, we introduce a Bayesian statistical method to jointly constrain the neutral fraction of the universe and the quasar lifetime from individual quasar spectra. We apply this methodology to the spectra of the two quasars with the highest redshifts known, ULAS J1120+0641 and ULAS J1342+0928, and measure volume-averaged neutral fractions {xH 1}(Z = 7.09) = 0.48+0.26-0.26 and {xH 1}(Z = 7.54) = 0.60+0.20-0.23(posterior medians and 68% credible intervals) when marginalized over quasar lifetimes of 103 ≤ t q ≤ 108 yr.

We present the first discoveries from a survey of z ≳ 6 quasars using imaging data from the DECam Legacy Survey (DECaLS) in the optical, the UKIRT Deep Infrared Sky Survey (UKIDSS) and a preliminary version of the UKIRT Hemisphere Survey (UHS) in the near-IR, and ALLWISE in the mid-IR. DECaLS will image 9000 deg2 of sky down to z AB ∼ 23.0, and UKIDSS and UHS will map the northern sky at 0 < decl. < +607deg;, reaching J VEGA ∼ 19.6 (5-σ). The combination of these data sets allows us to discover quasars at redshift z ≳ 7 and to conduct a complete census of the faint quasar population at z ≳ 6. In this paper, we report on the selection method of our search, and on the initial discoveries of two new, faint z ≳ 6 quasars and one new z = 6.63 quasar in our pilot spectroscopic observations. The two new z ∼ 6 quasars are at z = 6.07 and z = 6.17 with absolute magnitudes at rest-frame wavelength 1450 Å being M 1450 = -25.83 and M 1450 = -25.76, respectively. These discoveries suggest that we can find quasars close to or fainter than the break magnitude of the Quasar Luminosity Function (QLF) at z ≳ 6. The new z = 6.63 quasar has an absolute magnitude of M 1450 = -25.95. This demonstrates the potential of using the combined DECaLS and UKIDSS/UHS data sets to find z ≳ 7 quasars. Extrapolating from previous QLF measurements, we predict that these combined data sets will yield ∼200 z ∼ 6 quasars to z AB < 21.5, ∼1000 z ∼ 6 quasars to z AB < 23, and ∼30 quasars at z > 6.5 to J VEGA < 19.5.

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 × 1013Le. Deep spectroscopic observations reveal a SMBH with a mass of (1.6 ± 0.4) × 109M⊙ 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 ∼200M⊙ yr-1 and a dust mass of ∼7 × 107M⊙. Follow-up observations of this reionizationera 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 ("Poniua'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) × 109 M o 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 104 M o 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-1 in its host galaxy.

Radio sources at the highest redshifts can provide unique information on the first massive galaxies and black holes, the densest primordial environments, and the epoch of reionization. The number of astronomical objects identified at z > 6 has increased dramatically over the last few years, but previously only three radio-loud (R 2500 = f ν,5 GHz/f ν,2500 Å > 10) sources had been reported at z > 6, with the most distant being a quasar at z = 6.18. Here we present the discovery and characterization of PSO J172.3556+18.7734, a radio-loud quasar at z = 6.823. This source has an Mg ii-based black hole mass of ∼3 × 108 M o˙ and is one of the fastest accreting quasars, consistent with super-Eddington accretion. The ionized region around the quasar is among the largest measured at these redshifts, implying an active phase longer than the average lifetime of the z ⪆ 6 quasar population. From archival data, there is evidence that its 1.4 GHz emission has decreased by a factor of two over the last two decades. The quasar's radio spectrum between 1.4 and 3.0 GHz is steep (α = -1.31). Assuming the measured radio slope and extrapolating to rest-frame 5 GHz, the quasar has a radio-loudness parameter R 2500 ∼ 90. A second steep radio source (α = -0.83) of comparable brightness to the quasar is only 23.″1 away (∼120 kpc at z = 6.82; projection probability <2%), but shows no optical or near-infrared counterpart. Further follow-up is required to establish whether these two sources are physically associated.

Luminous z ≥ 7 quasars provide direct probes of the evolution of supermassive black holes (SMBHs) and the intergalactic medium (IGM) during the epoch of reionization (EoR). The Lyα damping wing absorption imprinted by neutral hydrogen in the IGM can be detected in a single EoR quasar spectrum, allowing the measurement of the IGM neutral fraction toward that line of sight. However, damping wing features have only been detected in two z > 7 quasars in previous studies. In this paper, we present new high-quality optical and near-infrared spectroscopy of the z = 7.00 quasar DES J025216.64-050331.8 obtained with Keck/Near-Infrared Echellette Spectrometer and Gemini/GMOS. By using the Mg ii single-epoch virial method, we find that it hosts a SMBH accreting at an Eddington ratio of λ Edd = 0.7 0.1, consistent with the values seen in other luminous z ∼ 7 quasars. Furthermore, the Lyα region of the spectrum exhibits a strong damping wing absorption feature. The lack of associated metal absorption in the quasar spectrum indicates that this absorption is imprinted by a neutral IGM. Using a state-of-the-art model developed by Davies et al., we measure a volume-averaged neutral hydrogen fraction at z = 7 of within 68% (95%) confidence intervals when marginalizing over quasar lifetimes of. This is the highest IGM neutral fraction yet measured using reionization-era quasar spectra.

The Dark Energy Spectroscopic Instrument (DESI) survey will measure large-scale structures using quasars as direct tracers of dark matter in the redshift range 0.9 < z < 2.1 and using Lyα forests in quasar spectra at z > 2.1. We present several methods to select candidate quasars for DESI, using input photometric imaging in three optical bands (g, r, z) from the DESI Legacy Imaging Surveys and two infrared bands (W1, W2) from the Wide-field Infrared Survey Explorer. These methods were extensively tested during the Survey Validation of DESI. In this paper, we report on the results obtained with the different methods and present the selection we optimized for the DESI main survey. The final quasar target selection is based on a random forest algorithm and selects quasars in the magnitude range of 16.5 < r < 23. Visual selection of ultra-deep observations indicates that the main selection consists of 71% quasars, 16% galaxies, 6% stars, and 7% inconclusive spectra. Using the spectra based on this selection, we build an automated quasar catalog that achieves a fraction of true QSOs higher than 99% for a nominal effective exposure time of ∼1000 s. With a 310 deg−2 target density, the main selection allows DESI to select more than 200 deg−2 quasars (including 60 deg−2 quasars with z > 2.1), exceeding the project requirements by 20%. The redshift distribution of the selected quasars is in excellent agreement with quasar luminosity function predictions.

The DESI Legacy Imaging Surveys (http://legacysurvey.org/) are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image ≈14,000 deg2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12, and 22 μm) observed by the Wide-field Infrared Survey Explorer satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.