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.

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.