In a multi-wavelength survey of $13$ quasars at $5.8\lesssim z\lesssim6.5$, that were pre-selected to be potentially young, we find five objects with extremely small proximity zone sizes that may imply UV-luminous quasar lifetimes of $\lesssim 100,000$ years. Proximity zones are regions of enhanced transmitted flux in the vicinity of the quasars that are sensitive to the quasars' lifetimes because the intergalactic gas has a finite response time to their radiation. We combine sub-mm observations from the Atacama Large Millimetre Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA), as well as deep optical and near-infrared spectra from medium-resolution spectrograph on the Very Large Telescope (VLT) 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 (SMBHs) 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\%\lesssim f_{\rm young}\lesssim 10\%$. One of the young objects, PSO J158-14, shows a very bright dust continuum flux ($F_{\rm cont}=3.46\pm 0.02\,\rm mJy$), indicating a highly star-bursting host galaxy with a star formation rate of approximately $1420\,M_{\odot}\,\rm yr^{-1}$.

We present deep (9 hours) Gemini-N/GNIRS near-infrared spectroscopic observations of ULAS J1342+0928, a luminous quasar at z=7.54. Various broad emission lines were detected, as well as the underlying continuum and iron forests over the rest-frame wavelength 970-2930A. There is a clear trend that higher-ionization emission lines show larger blueshifts with CIV1549 exhibiting 5510^{+240}_{-110} km s-1 blueshift with respect to the systematic redshift from the far-infrared [CII] 158um emission line. Those high ionization lines have wide profiles with FWHM more than 10000 km s-1. A modest blueshift of 340^{+110}_{-80} km s-1 is also seen in MgII, the lowest ionization line identified in the spectrum. The updated MgII-based black hole mass of M_BH=9.1_{-1.3}^{+1.4} x 10^8 M_sun and the Eddington ratio of L_bol/L_Edd=1.1_{-0.2}^{+0.2} confirm that ULAS J1342+0928 is powered by a massive and actively accreting black hole. There is no significant difference in the emission line ratios such as SiIV/CIV and AlIII/CIV when compared to lower-redshift quasars in a similar luminosity range, which suggests early metal pollution of the broad-line-region clouds. This trend also holds for the FeII/MgII line ratio, known as a cosmic clock that traces the iron enrichment in the early universe. Different iron templates and continuum fitting ranges were used to explore how the FeII/MgII measurement changes as a function of spectral modeling. Quasars at even higher redshift or at fainter luminosity range (L_bol<10^46 erg s-1) are needed to probe the sites of early metal enrichment and a corresponding change in the FeII/MgII ratio.

We report the discovery of a luminous quasar, J1007+2115 at $z=7.515$ ("P\={o}niu\={a}'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\pm0.2)\times10^9$ $M_{\odot}$ supermassive black hole (SMBH), based on its broad MgII 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 $\gtrsim 10^{4}$ $M_{\odot}$ 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$\alpha$ line profile, and find that it is weaker than that of J1342+0928 and two other $z\gtrsim7$ quasars. We estimate an IGM volume-averaged neutral fraction $\langle x\rm_{HI}\rangle=0.39^{+0.22}_{-0.13}$. This range of values suggests a patchy reionization history toward different IGM sightlines. We detect the 158 $\mu$m [C II] emission line in J1007+2115 with ALMA; this line centroid yields a systemic redshift of $z=7.5149\pm0.0004$ and indicates a star formation rate of $\sim210$ $M_{\odot}$ yr$^{-1}$ in its host galaxy.

Reconstruction techniques for intrinsic quasar continua are crucial for the precision study of Lyman-$\alpha$ (Ly-$\alpha$) and Lyman-$\beta$ (Ly-$\beta$) transmission at $z>5.0$, where the $\lambda<1215 A$ emission of quasars is nearly completely absorbed. While the number and quality of spectroscopic observations has become theoretically sufficient to quantify Ly-$\alpha$ transmission at $5.0<z<6.0$ to better than $1\%$, the biases and uncertainties arising from predicting the unabsorbed continuum are not known to the same level. In this paper, we systematically evaluate eight reconstruction techniques on a unified testing sample of $2.7<z<3.5$ quasars drawn from eBOSS. The methods include power-law extrapolation, stacking of neighbours, and six variants of Principal Component Analysis (PCA) using direct projection, fitting of components, or neural networks to perform weight mapping. We find that power-law reconstructions and the PCA with fewest components and smallest training sample display the largest biases in the Ly-$\alpha$ forest ($-9.58\%/+8.22\%$ respectively). Power-law extrapolations have larger scatters than previously assumed of $+13.1\%/-13.2\%$ over Ly-$\alpha$ and $+19.9\%/-20.1\%$ over Ly-$\beta$. We present two new PCAs which achieve the best current accuracies of $9\%$ for Ly-$\alpha$ and $17\%$ for Ly-$\beta$. We apply the eight techniques after accounting for wavelength-dependent biases and scatter to a sample $19$ quasars at $z>5.7$ with IR X-Shooter spectroscopy, obtaining well-characterised measurements for the mean flux transmission at $4.7<z<6.3$. Our results demonstrate the importance of testing and, when relevant, training, continuum reconstruction techniques in a systematic way.

PypeIt is a Python package for semi-automated reduction of astronomical, spectroscopic data. Its algorithms build on decades-long development of previous data reduction pipelines by the developers (Bernstein, Burles, & Prochaska, 2015; Bochanski et al., 2009). The reduction procedure -- including a complete list of the input parameters and available functionality -- is provided as online documentation hosted by Read the Docs, which is regularly updated. (https://pypeit.readthedocs.io/en/latest/). Release v1.0.3 serves the following spectrographs: Gemini/GNIRS, Gemini/GMOS, Gemini/FLAMINGOS 2, Lick/Kast, Magellan/MagE, Magellan/Fire, MDM/OSMOS, Keck/DEIMOS (600ZD, 830G, 1200G), Keck/LRIS, Keck/MOSFIRE (J and Y gratings tested), Keck/NIRES, Keck/NIRSPEC (low-dispersion), LBT/Luci-I, Luci-II, LBT/MODS (beta), NOT/ALFOSC (grism4), VLT/X-Shooter (VIS, NIR), VLT/FORS2 (300I, 300V), WHT/ISIS.

Observations of the early Universe suggest that reionization was complete by $z\sim6$, however, the exact history of this process is still unknown. One method for measuring the evolution of the neutral fraction throughout this epoch is via observing the Ly$\alpha$ damping wings of high-redshift quasars. In order to constrain the neutral fraction from quasar observations, one needs an accurate model of the quasar spectrum around Ly$\alpha$, after the spectrum has been processed by its host galaxy but before it is altered by absorption and damping in the intervening IGM. In this paper, we present a novel machine learning approach, using artificial neural networks, to reconstruct quasar continua around Ly$\alpha$. Our QSANNdRA algorithm improves the error in this reconstruction compared to the state-of-the-art PCA-based model in the literature by 14.2% on average, and provides an improvement of 6.1% on average when compared to an extension thereof. In comparison with the extended PCA model, QSANNdRA further achieves an improvement of 22.1% and 16.8% when evaluated on low-redshift quasars most similar to the two high-redshift quasars under consideration, ULAS J1120+0641 at $z=7.0851$ and ULAS J1342+0928 at $z=7.5413$, respectively. Using our more accurate reconstructions of these two $z>7$ quasars, we estimate the neutral fraction of the IGM using a homogeneous reionization model and find $\bar{x}_\mathrm{HI} = 0.25^{+0.05}_{-0.05}$ at $z=7.0851$ and $\bar{x}_\mathrm{HI} = 0.60^{+0.11}_{-0.11}$ at $z=7.5413$. Our results are consistent with the literature and favour a rapid end to reionization.

Luminous $z\ge7$ 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$\alpha$ 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 towards 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/NIRES and Gemini/GMOS. By using the MgII single-epoch virial method, we find that it hosts a $\rm (1.39\pm0.16) \times10^{9} ~M_\odot$ SMBH accreting at an Eddington ratio of $\lambda_{\rm Edd}=0.7\pm0.1$, consistent with the values seen in other luminous $z\sim 7$ quasars. Furthermore, the Ly$\alpha$ 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 $\langle x_{\rm HI} \rangle = 0.70^{+0.20}_{-0.23} (^{+0.28}_{-0.48})$ within 68% (95%) confidence intervals when marginalizing over quasar lifetimes of $10^3\le t_{\rm Q}\le10^8$ yr. This is the highest IGM neutral fraction yet measured using reionization-era quasar spectra.

We identify a strong Ly-alpha damping wing profile in the spectrum of the quasar P183+05 at z=6.4386. Given the detection of several narrow metal absorption lines at z=6.40392, the most likely explanation for the absorption profile is that it is due to a damped Ly-alpha system. However, in order to match the data a contribution of an intergalactic medium 5-38% neutral or additional weaker absorbers near the quasar is also required. The absorption system presented here is the most distant damped Ly-alpha system currently known. We estimate an HI column density ($10^{20.68\pm0.25}\,$cm$^{-2}$), metallicity ([O/H]$=-2.92\pm 0.32$), and relative chemical abundances of a system consistent with a low-mass galaxy during the first Gyr of the universe. This object is among the most metal-poor damped Ly-alpha systems known and, even though it is observed only ~850 Myr after the big bang, its relative abundances do not show signatures of chemical enrichment by Population III stars.

The discovery of quasars few hundred megayears after the Big Bang represents a major challenge to our understanding of black holes and galaxy formation and evolution. Their luminosity is produced by extreme gas accretion onto black holes, which already reached masses of 10$^9$ M$_\odot$ by z ~ 6. Simultaneously, their host galaxies form hundreds of stars per year, using up gas in the process. To understand which environments are able to sustain the rapid formation of these extreme sources we started a VLT/MUSE effort aimed at characterizing the surroundings of a sample of 5.7 < z < 6.6 quasars dubbed: the Reionization Epoch QUasar InvEstigation with MUSE (REQUIEM) survey. We here present results of our searches for extended Ly-Alpha halos around the first 31 targets observed as part of this program. Reaching 5-sigma surface brightness limits of 0.1-1.1 x 10$^{-17}$ erg/s/cm$^2$/arcsec$^2$ over a 1 arcsec$^2$ aperture, we were able to unveil the presence of 12 Ly-Alpha nebulae, 8 of which are newly discovered. The detected nebulae show a variety of emission properties and morphologies with luminosities ranging from 8 x 10$^{42}$ to 2 x 10$^{44}$ erg/s, FWHMs between 300 and 1700 km/s, sizes < 30 pkpc, and redshifts consistent with those of the quasar host galaxies. As the first statistical and homogeneous investigation of the circum-galactic medium of massive galaxies at the end of the reionization epoch, the REQUIEM survey enables the study of the evolution of the cool gas surrounding quasars in the first 3 Gyr of the Universe. A comparison with the extended Ly-Alpha emission observed around bright (M$_{1450}$ < -25 mag) quasars at intermediate redshift indicates little variations on the properties of the cool gas from z ~ 6 to z ~ 3 followed by a decline in the average surface brightness down to z ~ 2.

The supermassive black holes (SMBHs) observed at the centers of all massive galaxies are believed to have grown via luminous accretion during quasar phases in the distant past. The fraction of inflowing rest mass energy emitted as light, the radiative efficiency, has been inferred to be 10\%, in agreement with expectations from thin disk accretion models. But the existence of billion solar-mass SMBHs powering quasars at $z > 7$ challenges this picture: provided they respect the Eddington limit, there is not enough time to grow $z>7$ SMBHs from stellar remnant seeds unless the radiative efficiency is below 10\%. Here we show that one can constrain the radiative efficiencies of the most distant quasars known using foreground neutral intergalactic gas as a cosmological-scale ionizing photon counter. From the Ly$\alpha$ absorption profiles of ULAS J1120+0641 ($z=7.09$) and ULAS J1342+0928 ($z=7.54$), we determine posterior median radiative efficiencies of 0.08\% and 0.1\%, respectively, and the combination of the two measurements rule out the canonical 10\% efficiency at 99.8\% credibility after marginalizing over the unknown obscured fraction. This low radiative efficiency implies rapid mass accretion for the earliest SMBHs, greatly easing the tension between the age of the Universe and the SMBH masses. However, our measured efficiency may instead reflect nearly complete obscuration by dusty gas in the quasar host galaxies over the vast majority of their SMBH growth. Assuming 10\% efficiency during unobscured phases, we find that the obscured fraction would be $>82\%$ at 95\% credibility, and imply a $25.7^{+49.6}_{-16.5}$ times larger obscured than unobscured luminous quasar population at $z>7$.

We measure the intergalactic medium (IGM) opacity in the Ly$\alpha$ as well as in the Ly$\beta$ forest along $19$ quasar sightlines between $5.5\lesssim z_{\rm abs}\lesssim 6.1$, probing the end stages of the reionization epoch. Owing to its lower oscillator strength the Ly$\beta$ transition is sensitive to different gas temperatures and densities than Ly$\alpha$, providing additional constraints on the ionization and thermal state of the IGM. A comparison of our measurements to different inhomogeneous reionization models, derived from post-processing the Nyx cosmological hydrodynamical simulation to include spatial fluctuations in the ultraviolet background (UVB) or the gas temperature field, as well as to a uniform reionization model with varying thermal states of the IGM, leads to two primary conclusions: First, we find that including the effects of spectral noise is key for a proper data to model comparison. Noise effectively reduces the sensitivity to high opacity regions, and thus even stronger spatial inhomogeneities are required to match the observed scatter in the observations than previously inferred. Second, we find that models which come close to reproducing the distribution of Ly$\alpha$ effective optical depths nevertheless underpredict the Ly$\beta$ opacity at the same spatial locations. The origin of this disagreement is not entirely clear but models with an inversion in the temperature-density relation of the IGM just after reionization is completed match our measurements best, although they still do not fully capture the observations at $z\gtrsim 5.8$.

In this work we present a new hybrid method to simulate the thermal effects of the reionization in cosmological hydrodynamical simulations. The method improves upon the standard approach used in simulations of the intergalactic medium (IGM) and galaxy formation without a significant increase of the computational cost allowing for efficient exploration of the parameter space. The method uses a small set of phenomenological input parameters and combines a semi-numerical reionization model to solve for the topology of reionization and an approximate model of how reionization heats the IGM, with the massively parallel \texttt{Nyx} hydrodynamics code, specifically designed to solve for the structure of diffuse IGM gas. We have produced several large-scale high resolution cosmological hydrodynamical simulations ($2048^3$, $L_{\rm box} = 40$ Mpc/h) with different instantaneous and inhomogeneous HI reionization models that use this new methodology. We study the IGM thermal properties of these models and find that large scale temperature fluctuations extend well beyond the end of reionization. Analyzing the 1D flux power spectrum of these models, we find up to $\sim 50\%$ differences in the large scale properties (low modes, $k\lesssim0.01$ s/km) of the post-reionization power spectrum due to the thermal fluctuations. We show that these differences could allow one to distinguish between different reionization scenarios already with existing Ly$\alpha$ forest measurements. Finally, we explore the differences in the small-scale cutoff of the power spectrum and we find that, for the same heat input, models show very good agreement provided that the reionization redshift of the instantaneous reionization model happens at the midpoint of the inhomogeneous model.

The larger-than-expected scatter in the opacity of the Ly$\alpha$ forest suggests that the metagalactic ionizing background is strongly fluctuating at $z > 5.5$. Models for ionizing background fluctuations predict a strong positive bias on large scales, so the environments of massive $>10^{12}\,{\rm M}_\odot$ dark matter halos, e.g. $z\sim6$ quasar hosts, would be ideal laboratories to constrain the sources of ionizing photons. While the quasars themselves should overwhelm any plausible ionizing photon contribution from neighboring galaxies, proximate damped Ly$\alpha$ absorbers (DLAs) have recently been discovered in the foreground of $z\sim6$ quasars, and the Ly$\alpha$ forest in the shadow of these DLAs could probe the local ionization environment. Using Gpc$^3$ simulations of $z=6$ ionizing background fluctuations, we show that while the Ly$\alpha$ forest signal from ionization bias around a quasar host halo should be strong, it is likely suppressed by the associated intergalactic matter overdensity. We also show that the quasar itself may still overwhelm the clustering signal via a "ghost" of the proximity effect from the quasar radiation causing a large-scale bias in the ionizing photon mean free path. This ghost proximity effect is sensitive to the lifetime and geometry of quasar emission, potentially unlocking a new avenue for constraining these fundamental quasar properties. Finally, we present observations of a $z\sim6$ quasar with a proximate DLA which shows a strong excess in Ly$\alpha$ forest transmission at the predicted location of the ghost proximity effect.

We present measurements of the large-scale (~40 comoving Mpc) effective optical depth of HeII Ly$\alpha$ absorption, $\tau_{\rm eff}$, at 2.54<z<3.86 toward 16 HeII-transparent quasars observed with the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST), to characterize the ionization state of helium in the intergalactic medium (IGM). We provide the first statistical sample of $\tau_{\rm eff}$ measurements in six signal-to-noise ratio >3 HeII sightlines at z>3.5, and study the redshift evolution and sightline-to-sightline variance of $\tau_{\rm eff}$ in 24 HeII sightlines. We confirm an increase of the median $\tau_{\rm eff}$ from ~2 at z=2.7 to >5 at z>3, and a scatter in $\tau_{\rm eff}$ that increases with redshift. The z>3.5 HeII absorption is predominantly saturated, but isolated narrow ($\Delta v<650$ km/s) transmission spikes indicate patches of reionized helium. We compare our measurements to predictions for a range of UV background models applied to outputs of a large-volume (146 comoving Mpc)$^3$ hydrodynamical simulation by forward-modeling our sample's quality and size. At z>2.74 the variance in $\tau_{\rm eff}$ significantly exceeds expectations for a spatially uniform UV background, but is consistent with a fluctuating radiation field sourced by variations in the quasar number density and the mean free path in the post-reionization IGM. We develop a method to infer the approximate median HeII photoionization rate $\Gamma_{\rm HeII}$ of a fluctuating UV background from the median $\tau_{\rm eff}$, finding a factor ~5 decrease in $\Gamma_{\rm HeII}$ between z~2.6 and z~3.1. At z~3.1 a $\Gamma_{\rm HeII}=\left[9.1^{+1.1}_{-1.2}\,\mathrm{(stat.)}\,^{+2.4}_{-3.4}\,\mathrm{(sys.)}\right]\times10^{-16}\,\mathrm{s}^{-1}$ corresponds to a median HeII fraction of ~2.5%, indicating that our data probe the tail end of HeII reionization.

During reionization, the intergalactic medium is heated impulsively by supersonic ionization fronts (I-fronts). The peak gas temperatures behind the I-fronts, $T_\mathrm{reion}$, are a key uncertainty in models of the thermal history after reionization. Here we use high-resolution radiative transfer simulations to study the parameter space of $T_\mathrm{reion}$. We show that $T_\mathrm{reion}$ is only mildly sensitive to the spectrum of incident radiation over most of the parameter space, with temperatures set primarily by I-front speeds. We also explore what current models of reionization predict for $T_\mathrm{reion}$ by measuring I-front speeds in cosmological radiative transfer simulations. We find that the post-I-front temperatures evolve toward hotter values as reionization progresses. Temperatures of $T_\mathrm{reion} = 17,000-22,000$ K are typical during the first half of reionization, but $T_\mathrm{reion} = 25,000 - 30,000$ K may be achieved near the end of this process if I-front speeds reach $\sim10^4$ km/s as found in our simulations. Shorter reionization epochs lead to hotter $T_\mathrm{reion}$. We discuss implications for $z>5$ Ly$\alpha$ forest observations, which potentially include sight lines through hot, recently reionized patches of the Universe. Interpolation tables from our parameter space study are made publicly available, along with a simple fit for the dependence of $T_\mathrm{reion}$ on the I-front speed.

Since the discovery of $z\sim 6$ quasars two decades ago, studies of their Ly$\alpha$-transparent proximity zones have largely focused on their utility as a probe of cosmic reionization. But even when in a highly ionized intergalactic medium, these zones provide a rich laboratory for determining the timescales that govern quasar activity and the concomitant growth of their supermassive black holes. In this work, we use a suite of 1D radiative transfer simulations of quasar proximity zones to explore their time-dependent behaviour for activity timescales from $\sim10^3$ to $10^8$ years. The sizes of the simulated proximity zones, as quantified by the distance at which the smoothed Ly$\alpha$ transmission drops below 10% (denoted $R_p$), are in excellent agreement with observations, with the exception of a handful of particularly small zones that have been attributed to extremely short $\lesssim 10^4$ year lifetimes. We develop a physically motivated semi-analytic model of proximity zones which captures the bulk of their equilibrium and non-equilibrium behaviour, and use this model to investigate how quasar variability on $\lesssim10^5$ year timescales is imprinted on the distribution of observed proximity zone sizes. We show that large variations in the ionizing luminosity of quasars on timescales of $\lesssim10^4$ years are disfavored based on the good agreement between the observed distribution of $R_p$ and our model prediction based on "lightbulb" (i.e. steady constant emission) light curves.

The early phases of galaxy formation constitute one of the most exciting frontiers in astrophysics. It is during this era that the first luminous sources reionize the intergalactic medium - the moment when structure formation affects every baryon in the Universe. Here we argue that we will obtain a complete picture of this era by combining observations of galaxies with direct measurements of the reionization process: the former will provide a detailed understanding of bright sources, while the latter will constrain the (substantial) faint source population. We further describe how optimizing the comparison of these two measurements requires near-infrared galaxy surveys covering large volumes and retaining redshift information and also improvements in 21-cm analysis, moving those experiments into the imaging regime.

The epoch of reionization, when photons from early galaxies ionized the intergalactic medium about a billion years after the Big Bang, is the last major phase transition in the Universe's history. Measuring the characteristics of the transition is important for understanding early galaxies and the cosmic web and for modeling dwarf galaxies in the later Universe. But such measurements require probes of the intergalactic medium itself. Here we describe how the 21-cm line of neutral hydrogen provides a powerful probe of the reionization process and therefore important constraints on both the galaxies and intergalactic absorbers at that time. While existing experiments will make precise statistical measurements over the next decade, we argue that improved 21-cm analysis techniques - allowing imaging of the neutral gas itself - as well as improved theoretical models, are crucial for testing our understanding of this important era.

The Dark Ages are the period between the last scattering of the cosmic microwave background and the appearance of the first luminous sources, spanning approximately 1100 < z < 30. The only known way to measure fluctuations in this era is through the 21-cm line of neutral hydrogen. Such observations have enormous potential for cosmology, because they span a large volume while the fluctuations remain linear even on small scales. Observations of 21-cm fluctuations during this era can therefore constrain fundamental aspects of our Universe, including inflation and any exotic physics of dark matter. While the observational challenges to these low-frequency 21-cm observations are enormous, especially from the terrestrial environment, they represent an important goal for cosmology.

Absorption signatures in the spectra of QSOs are one of our most powerful tools for studying galactic and intergalactic environments at high redshifts. With the discovery of QSOs out to z > 7, QSO absorption lines are now tracing the end stages of reionization on multiple fronts using the hydrogen Lyman-$\alpha$ forest and heavy element absorbers. Next-generation QSO absorption line studies with large optical/IR telescopes will reveal in detail how the first galaxies emerged form the cosmic web, transformed their circum- and inter-galactic environments, and completed the last major phase transition of the Universe. These efforts will complement other upcoming studies of reionization, such as those with JWST, ALMA, and redshifted 21cm experiments.

The cosmic dawn and epoch of reionization mark the time period in the universe when stars, galaxies, and blackhole seeds first formed and the intergalactic medium changed from neutral to an ionized one. Despite substantial progress with multi-wavelength observations, astrophysical process during this time period remain some of the least understood with large uncertainties on our existing models of galaxy, blackhole, and structure formation. This white paper outlines the current state of knowledge and anticipated scientific outcomes with ground and space-based astronomical facilities in the 2020s. We then propose a number of scientific goals and objectives for new facilities in late 2020s to mid 2030s that will lead to definitive measurements of key astrophysical processes in the epoch of reionization and cosmic dawn.

Despite extensive efforts, only two quasars have been found at $z>7$ to date due to a combination of low spatial density and high contamination from more ubiquitous Galactic cool dwarfs in quasar selection. This limits our current knowledge of the super-massive black hole (SMBH) growth mechanism and reionization history. In this letter, we report the discovery of a luminous quasar at $z=7.021$, DELS J003836.10$-$152723.6 (hereafter J0038$-$1527), selected using photometric data from DESI Legacy imaging Survey (DELS), Pan-STARRS1 (PS1) imaging Survey, as well as Wide-field Infrared Survey Explore ($WISE$) mid-infrared all-sky survey. With an absolute magnitude of $M_{1450}$=$-$27.1 and bolometric luminosity of $L_{\rm Bol}$=5.6$\times$10$^{13}$ $L_\odot$, J0038$-$1527 is the most luminous quasar known at $z>7$. Deep optical to near infrared spectroscopic observations suggest that J0038-1527 hosts a 1.3 billion solar mass BH accreting at the Eddington limit, with an Eddington ratio of 1.25$\pm$0.19. The CIV broad emission line of J0038$-$1527 is blue-shifted by more than 3000 km s$^{-1}$ to the systemic redshift. More detailed investigations of the high quality spectra reveal three extremely high velocity CIV broad absorption lines (BALs) with velocity from 0.08 to 0.14 times the speed of light and total balnicity index of more than 5000 km s$^{-1}$, suggesting the presence of relativistic outflows. J0038$-$1527 is the first quasar found at the epoch of reionization (EoR) with such strong outflows and provides a unique laboratory to investigate AGN feedback on the formation and growth of the most massive galaxies in the early universe.

This is the third paper in a series aims at finding reionzation-era quasars with the combination of DESI Legacy imaging Surveys (DELS) and near-infrared imaging surveys, such as the UKIRT Hemisphere Survey (UHS), as well as the Wide-field Infrared Survey Explore ($WISE$) mid-infrared survey. In this paper, we describe the updated quasar candidate selection procedure, report the discovery of 16 quasars at $6.4\lesssim z \lesssim6.9$ from area of $\sim$13,020 deg$^2$, and present the quasar luminosity function (QLF) at $z\sim6.7$. The measured QLF follows $\Phi(L_{1450})\propto L_{1450}^{-2.35}$ in the magnitude range $27.6<M_{1450}<-25.5$. We determine the quasar comoving spatial density at $\langle z \rangle$=6.7 and $M_{1450}<-26.0$ to be $\rm 0.39\pm0.11 Gpc^{-3}$ and find that the exponential density evolution parameter to be $k=-0.78\pm0.18$ from $z\sim6$ to $z\sim6.7$, corresponding to a rapid decline by a factor of $\sim 6$ per unit redshift towards earlier epoch, a rate significantly faster than that at $z\sim 3- 5$. The cosmic time between $z\sim6$ and $z\sim6.7$ is only 121 Myrs. The quasar density declined by a factor of more than three within such short time requires that SMBHs must grow rapidly or they are less radiatively efficient at higher redshifts. We measured quasar comoving emissivity at $z\sim6.7$ which indicate that high redshift quasars are highly unlikely to make a significant contribution to hydrogen reionization. The broad absorption line (BAL) quasar fraction at $z\gtrsim6.5$ is measured to be $\gtrsim$22%. In addition, we also report the discovery of additional five quasars at $z\sim6$ in the appendix.

We publicly release a new sample of $34$ medium resolution quasar spectra at $5.77\leq z_{\rm em}\leq6.54$ observed with the Echellette Spectrograph and Imager (ESI) on the Keck telescope. This quasar sample represents an ideal laboratory to study the intergalactic medium (IGM) during the end stages of the epoch of reionization, and constrain the timing and morphology of the phase transition. For a subset of $23$ of our highest signal-to-noise ratio spectra (S/N$>7$, per $10\,{\rm km\,s^{-1}}$ pixel), we present a new measurement of the Lyman-$\alpha$ (Ly$\alpha$) forest opacity spanning the redshift range $4.8\lesssim z\lesssim6.3$. We carefully eliminate spectral regions that could be causing biases in our measurements due to additional transmitted flux in the proximity zone of the quasars, or extra absorption caused by strong intervening absorption systems along the line of sight. We compare the observed evolution of the IGM opacity with redshift to predictions from a hydrodynamical simulation with uniform ultraviolet background (UVB) radiation, as well as two semi-numerical patchy reionization models, one with a fluctuating UVB and another with a fluctuating temperature field. Our measurements show a steep rise in opacity at $z\gtrsim5.0$ and an increased scatter and thus support the picture of a spatially inhomogeneous reionization process, consistent with previous work. However, we measure significantly higher optical depths at $5.3\lesssim z\lesssim5.7$ than previous studies, which reduces the contrast between the highest opacity Gunn-Peterson troughs and the average opacity trend of the IGM, which may relieve some of the previously noted tension between these measurements and reionization models.

Quasars are the most luminous non-transient objects known and as a result they enable studies of the Universe at the earliest cosmic epochs. Despite extensive efforts, however, the quasar ULAS J1120+0641 at z=7.09 has remained the only one known at z>7 for more than half a decade. Here we report observations of the quasar ULAS J134208.10+092838.61 (hereafter J1342+0928) at redshift z=7.54. This quasar has a bolometric luminosity of 4e13 times the luminosity of the Sun and a black hole mass of 8e8 solar masses. The existence of this supermassive black hole when the Universe was only 690 million years old---just five percent of its current age---reinforces models of early black-hole growth that allow black holes with initial masses of more than about 1e4 solar masses or episodic hyper-Eddington accretion. We see strong evidence of absorption of the spectrum of the quasar redwards of the Lyman alpha emission line (the Gunn-Peterson damping wing), as would be expected if a significant amount (more than 10 per cent) of the hydrogen in the intergalactic medium surrounding J1342+0928 is neutral. We derive a significant fraction of neutral hydrogen, although the exact fraction depends on the modelling. However, even in our most conservative analysis we find a fraction of more than 0.33 (0.11) at 68 per cent (95 per cent) probability, indicating that we are probing well within the reionization epoch of the Universe.

The quasar lifetime $t_{\rm\,Q}$ is one of the most fundamental quantities for understanding quasar evolution and the growth of supermassive black holes (SMBHs), but remains uncertain by several orders of magnitude. In a recent study we uncovered a population of very young quasars ($t_{\rm Q}\lesssim10^4-10^5$ yr), based on the sizes of their proximity zones, which are regions of enhanced Ly$\alpha$ forest transmission near the quasar resulting from its own ionizing radiation. The presence of such young objects poses significant challenges to models of SMBH formation, which already struggle to explain the existence of SMBHs at such early cosmic epochs. We conduct the first comprehensive spectroscopic study of the youngest quasar known, $\rm SDSS\,J1335+3533$ at $z=5.9012$, whose lifetime is $t_{\rm Q}<10^4$ yr ($95\%$ confidence). A careful search of our deep optical and near-infrared spectra for HI and metal absorption lines allows us to convincingly exclude that its small proximity zone results from an associated absorption system rather than a short lifetime. We use the MgII emission line to measure its black hole mass $M_{\rm BH}=(4.09\pm0.58)\times 10^9 M_{\odot}$, implying an Eddington ratio of $0.30\pm0.04$ -- comparable to other co-eval quasars. We similarly find that the relationship between its black hole mass and dynamical mass are consistent with other $z\sim6$ quasars. The only possible anomaly associated with youth are its weak emission lines, but larger samples are needed to shed light on a potential causal connection. We discuss the implications of short lifetimes for various SMBH growth scenarios, and argue that future observations of young quasars with JWST could distinguish between them.

The observed scatter in intergalactic Lyman-$\alpha$ opacity at $z \lesssim 6$ requires large-scale fluctuations in the neutral fraction of the intergalactic medium (IGM) after the expected end of reionization. Post-reionization models that explain this scatter invoke fluctuations in either the ionizing ultraviolet background (UVB) or IGM temperature. These models make very different predictions, however, for the relationship between Lyman-$\alpha$ opacity and local density. Here we test these models using Lyman-$\alpha$ emitting galaxies (LAEs) to trace the density field surrounding the longest and most opaque known Lyman-$\alpha$ trough at $z < 6$. Using deep Subaru Hyper Suprime-Cam narrow-band imaging, we find a highly significant deficit of $z \simeq 5.7$ LAEs within 20 Mpc/$h$ of the trough. The results are consistent with a model in which the scatter in Lyman-$\alpha$ opacity near $z \sim 6$ is driven by large-scale UVB fluctuations, and disfavor a scenario in which the scatter is primarily driven by variations in IGM temperature. UVB fluctuations at this epoch present a boundary condition for reionization models, and may help shed light on the nature of the ionizing sources.

Recent observations have shown that the scatter in opacities among coeval segments of the Lyman-alpha forest increases rapidly at z > 5. In this paper, we assess whether the large scatter can be explained by fluctuations in the ionizing background in the post-reionization intergalactic medium. We find that matching the observed scatter at z ~ 5.5 requires a short spatially averaged mean free path of < 15 comoving Mpc/h, a factor of > 3 shorter than direct measurements at z ~ 5.2. We argue that such rapid evolution in the mean free path is difficult to reconcile with our measurements of the global H I photoionization rate, which stay approximately constant over the interval z ~ 4.8 - 5.5. However, we also show that measurements of the mean free path at z > 5 are likely biased towards higher values by the quasar proximity effect. This bias can reconcile the short values of the mean free path that are required to explain the large scatter in opacities. We discuss the implications of this scenario for cosmological reionization. Finally, we investigate whether other statistics applied to the z > 5 Lyman-alpha forest can shed light on the origin of the scatter. Compared to a model with a uniform ionizing background, models that successfully account for the scatter lead to enhanced power in the line-of-sight flux power spectrum on scales k < 0.1 h/Mpc. We find tentative evidence for this enhancement in observations of the high-redshift Lyman-alpha forest.

The central radio source in M87 provides the best opportunity to study jet formation because it has a large angular size for the gravitational radius of the black hole and has a bright jet that is well resolved by VLBI observations. We present intensive monitoring observations from 2007 and 2008, plus roughly annual observations that span 17 years, all made with the the Very Long Baseline Array at 43 GHz with a resolution of about 30 by 60 Rs. Our high-dynamic-range images clearly show the wide-opening-angle structure and the counter-jet. The jet and counter-jet are nearly symmetric in the inner 1.5 milli-arcseconds (mas; 0.12 pc in projection) with both being edge brightened. Both show deviations from parabolic shape in the form of an initial rapid expansion and subsequent contraction followed by further rapid expansion and, beyond the visible counter-jet, subsequent collimation. Proper motions and counter-jet/jet intensity ratios both indicate acceleration from apparent speeds of $\lesssim 0.5c$ to $\gtrsim 2c$ in the inner about 2 mas and suggest a helical flow. The jet displays a sideways shift with an approximately 8 to 10 year quasi-periodicity. The shift propagates outwards non-ballistically and significantly more slowly than the flow speed revealed by the fastest moving components. Polarization data show a systematic structure with magnetic field vectors that suggest a toroidal field close to the core.

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_{\rm HI}(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$\alpha$. 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 highest redshift quasars known, ULAS J1120+0641 and ULAS J1342+0928, and measured volume-averaged neutral fractions $\langle x_{\rm HI} \rangle(z=7.09)=0.48^{+0.26}_{-0.26}$ and $\langle x_{\rm HI} \rangle(z=7.54)=0.60^{+0.20}_{-0.23}$ (posterior medians and 68% credible intervals) when marginalized over quasar lifetimes of $10^3 \leq t_{\rm q} \leq 10^8$ years.

Measuring the proximity effect and the damping wing of intergalactic neutral hydrogen in quasar spectra during the epoch of reionization requires an estimate of the intrinsic continuum at rest-frame wavelengths $\lambda_{\rm rest}\sim1200$-$1260$ {\AA}. In contrast to previous works which used composite spectra with matched spectral properties or explored correlations between parameters of broad emission lines, we opted for a non-parametric predictive approach based on principal component analysis (PCA) to predict the intrinsic spectrum from the spectral properties at redder (i.e. unabsorbed) wavelengths. We decomposed a sample of $12764$ spectra of $z\sim2$-$2.5$ quasars from SDSS/BOSS into 10 red-side ($1280$ {\AA} $<\lambda_{\rm rest}<2900$ {\AA}) and 6 blue-side ($1180$ {\AA} $<\lambda_{\rm rest}<1280$ {\AA}) PCA basis spectra, and constructed a projection matrix to predict the blue-side coefficients from a fit to the red-side spectrum. We found that our method predicts the blue-side continuum with $\sim6$-$12\%$ precision and $\lesssim1\%$ bias by testing on the full training set sample. We then computed predictions for the blue-side continua of the two quasars currently known at $z>7$: ULAS J1120+0641 ($z=7.09$) and ULAS J1342+0928 ($z=7.54$). Both of these quasars are known to exhibit extreme emission line properties, so we individually calibrated the precision of the continuum predictions from similar quasars in the training set. We find that both $z>7$ quasars, and in particular ULAS J1342+0928, show signs of damping wing-like absorption at wavelengths redward of Ly$\alpha$.

The HeII transverse proximity effect - enhanced HeII Ly{\alpha} transmission in a background sightline caused by the ionizing radiation of a foreground quasar - offers a unique opportunity to probe the emission properties of quasars, in particular the emission geometry (obscuration, beaming) and the quasar lifetime. Building on the foreground quasar survey published in Schmidt+2017, we present a detailed model of the HeII transverse proximity effect, specifically designed to include light travel time effects, finite quasar ages, and quasar obscuration. We post-process outputs from a cosmological hydrodynamical simulation with a fluctuating HeII UV background model, plus the added effect of the radiation from a single bright foreground quasar. We vary the age $t_\mathrm{age}$ and obscured sky fractions $\Omega_\mathrm{obsc}$ of the foreground quasar, and explore the resulting effect on the HeII transverse proximity effect signal. Fluctuations in IGM density and the UV background, as well as the unknown orientation of the foreground quasar, result in a large variance of the HeII Ly{\alpha} transmission along the background sightline. We develop a fully Bayesian statistical formalism to compare far UV HeII Ly{\alpha} transmission spectra of the background quasars to our models, and extract joint constraints on $t_\mathrm{age}$ and $\Omega_\mathrm{obsc}$ for the six Schmidt+2017 foreground quasars with the highest implied HeII photoionization rates. Our analysis suggests a bimodal distribution of quasar emission properties, whereby one foreground quasar, associated with a strong HeII transmission spike, is relatively old $(22\,\mathrm{Myr})$ and unobscured $\Omega_\mathrm{obsc}<35\%$, whereas three others are either younger than $(10\,\mathrm{Myr})$ or highly obscured $(\Omega_\mathrm{obsc}>70\%)$.

Recent observations have discovered long (up to ~110 Mpc/h), opaque Gunn-Peterson troughs in the z ~ 5.5 Lyman-alpha forest, which are challenging to explain with conventional models of the post-reionization intergalactic medium. Here we demonstrate that observations of the galaxy populations in the vicinity of the deepest troughs can distinguish two competing models for these features: deep voids where the ionizing background is weak due to fluctuations in the mean free path of ionizing photons would show a deficit of galaxies, while residual temperature variations from extended, inhomogeneous reionization would show an overdensity of galaxies. We use large (~550 Mpc/h) semi-numerical simulations of these competing explanations to predict the galaxy populations in the largest of the known troughs at z ~ 5.7. We quantify the strong correlation of Lyman-alpha effective optical depth and galaxy surface density in both models and estimate the degree to which realistic surveys can measure such a correlation. While a spectroscopic galaxy survey is ideal, we also show that a relatively inexpensive narrowband survey of Lyman-alpha-emitting galaxies is ~90% likely to distinguish between the competing models.

We present a new method to quantify the value of the escape fraction of ionizing photons, and the existence of ultra-faint galaxies clustered around brighter objects during the epoch of cosmic reionization, using the diffuse Ly$\alpha$, continuum and H$\alpha$ emission observed around galaxies at $z\sim6$. We model the surface brightness profiles of the diffuse halos considering the fluorescent emission powered by ionizing photons escaping from the central galaxies, and the nebular emission from satellite star-forming sources, by extending the formalisms developed in Mas-Ribas & Dijkstra (2016) and Mas-Ribas et al. (2017). The comparison between our predicted profiles and Ly$\alpha$ observations at $z=5.7$ and $z=6.6$ favors a low ionizing escape fraction, $f_{\rm esc}^{\rm ion}\sim5\%$, for galaxies in the range $-19\gtrsim M_{\rm UV} \gtrsim -21.5$. However, uncertainties and possible systematics in the observations do not allow for firm conclusions. We predict H$\alpha$ and rest-frame visible continuum observations with JWST, and show that JWST will be able to detect extended (a few tens of kpc) fluorescent H$\alpha$ emission powered by ionizing photons escaping from a bright, $L\gtrsim 5L^*$, galaxy. Such observations can differentiate fluorescent emission from nebular emission by satellite sources. We discuss how observations and stacking of several objects may provide unique constraints on the escape fraction for faint galaxies and/or the abundance of ultra-faint radiation sources.

We study quasar proximity zones in the redshift range $5.77 \leq z \leq 6.54$ by homogeneously analyzing $34$ medium resolution spectra, encompassing both archival and newly obtained data, and exploiting recently updated systemic redshift and magnitude measurements. Whereas previous studies found strong evolution of proximity zone sizes with redshift, and argued that this provides evidence for a rapidly evolving intergalactic medium (IGM) neutral fraction during reionization, we measure a much shallower trend $\propto(1+z)^{-1.44}$. We compare our measured proximity zone sizes to predictions from hydrodynamical simulations post-processed with one-dimensional radiative transfer, and find good agreement between observations and theory irrespective of the ionization state of the ambient IGM. This insensitivity to IGM ionization state has been previously noted, and results from the fact that the definition of proximity zone size as the first drop of the smoothed quasar spectrum below the $10\%$ flux transmission level probes locations where the ionizing radiation from the quasar is an order of magnitude larger than the expected ultraviolet ionizing background that sets the neutral fraction of the IGM. Our analysis also uncovered three objects with exceptionally small proximity zones (two have $R_p < 1$proper Mpc), which constitute outliers from the observed distribution and are challenging to explain with our radiative transfer simulations. We consider various explanations for their origin, such as strong absorption line systems associated with the quasar or patchy reionization, but find that the most compelling scenario is that these quasars have been shining for $\lesssim 10^5$yr.

Large variations in the effective optical depth of the He II Ly$\alpha$ forest have been observed at $z\gtrsim2.7$, but the physical nature of these variations is uncertain: either the Universe is still undergoing the process of He II reionization, or the Universe is highly ionized but the He II-ionizing background fluctuates significantly on large scales. In an effort to build upon our understanding of the latter scenario, we present a novel model for the evolution of ionizing background fluctuations. Previous models have assumed the mean free path of ionizing photons to be spatially uniform, ignoring the dependence of that scale on the local ionization state of the intergalactic medium (IGM). This assumption is reasonable when the mean free path is large compared to the average distance between the primary sources of He II-ionizing photons, $\gtrsim L_\star$ quasars. However, when this is no longer the case, the background fluctuations become more severe, and an accurate description of the average propagation of ionizing photons through the IGM requires additionally accounting for the fluctuations in opacity. We demonstrate the importance of this effect by constructing 3D semi-analytic models of the helium ionizing background from $z=2.5$-$3.5$ that explicitly include a spatially varying mean free path of ionizing photons. The resulting distribution of effective optical depths at large scales in the He II Ly$\alpha$ forest is very similar to the latest observations with HST/COS at $2.5 \lesssim z \lesssim 3.5$.

The amplitude of the ionizing background that pervades the intergalactic medium (IGM) at the end of the epoch of reionization provides a valuable constraint on the emissivity of the sources which reionized the Universe. While measurements of the ionizing background at lower redshifts rely on a simulation-calibrated mapping between the photoionization rate and the mean transmission of the Ly$\alpha$ forest, at $z\gtrsim6$ the IGM becomes increasingly opaque, and transmission arises solely in narrow spikes separated by saturated Gunn-Peterson troughs. In this regime, the traditional approach of measuring the average transmission over large $\sim 50$ Mpc$/h$ regions is less sensitive and sub-optimal. Additionally, the five times smaller oscillator strength of the Ly$\beta$ transition implies the Ly$\beta$ forest is considerably more transparent at $z\gtrsim6$, even in the presence of contamination by foreground $z\sim 5$ Ly$\alpha$ forest absorption. In this work we present a novel statistical approach to analyze the joint distribution of transmission spikes in the co-spatial $z\sim 6$ Ly$\alpha$ and Ly$\beta$ forests. Our method relies on Approximate Bayesian Computation (ABC), which circumvents the necessity of computing the intractable likelihood function describing the highly correlated Ly$\alpha$ and Ly$\beta$ transmission. We apply ABC to mock data generated from a large-volume hydrodynamical simulation combined with a state-of-the-art model of ionizing background fluctuations in the post-reionization IGM, and show that it is sensitive to higher IGM neutral hydrogen fractions than previous techniques. As a proof of concept, we apply this methodology to a real spectrum of a $z=6.54$ quasar and measure the ionizing background from $5.4\leq z \leq 6.4$ along this sightline with $\sim0.2$ dex statistical uncertainties.

M87 is the best source in which to study a jet at high resolution in gravitational units because it has a very high mass black hole and is nearby. The angular size of the black hole is second only to Sgr A*, which does not have a strong jet. The jet structure is edge brightened with a wide opening angle base and a weak counterjet. We have roughly annual observations for 17 years plus intensive monitoring at three week intervals for a year and five day intervals for 2.5 months made with the Very Long Baseline Array (VLBA) at 43 GHz. The inner jet shows very complex dynamics, with apparent motions both along and across the jet. Speeds from zero to over 2c are seen, with acceleration observed over the first 3 milli-arcseconds. The counterjet decreases in brightness much more rapidly than the main jet, as is expected from relativistic beaming in an accelerating jet oriented near the line-of-sight. Details of the structure and dynamics are discussed. The roughly annual observations show side-to-side motion of the whole jet with a characteristic time scale of about 9 years.

The ionizing background of cosmic hydrogen is an important probe of the sources and absorbers of ionizing radiation in the post-reionization universe. Previous studies show that the ionization rate should be very sensitive to changes in the source population: as the emissivity rises, absorbers shrink in size, increasing the ionizing mean free path and, hence, the ionizing background. By contrast, observations of the ionizing background find a very flat evolution from z~2-5, before falling precipitously at z~6. We resolve this puzzling discrepancy by pointing out that, at z~2-5, optically thick absorbers are associated with the same collapsed halos that host ionizing sources. Thus, an increasing abundance of galaxies is compensated for by a corresponding increase in the absorber population, which moderates the instability in the ionizing background. However, by z~5-6, gas outside of halos dominates the absorption, the coupling between sources and absorbers is lost, and the ionizing background evolves rapidly. Our halo based model reproduces observations of the ionizing background, its flatness and sudden decline, as well as the redshift evolution of the ionizing mean free path. Our work suggests that, through much of their history, both star formation and photoelectric opacity in the universe track halo growth.

Extremely large opaque troughs in the Lyman-alpha forest have been interpreted as a sign of an extended reionization process below z~6. Such features are impossible to reproduce with simple models of the intergalactic ionizing background that assume a uniform mean free path of ionizing photons. We build a self-consistent model of the ionizing background that includes fluctuations in the mean free path due to the varying strength of the ionizing background and large-scale density field. The dominant effect is the suppression of the ionizing background in large-scale voids due to "self-shielding" by an enhanced number of optically thick absorbers. Our model results in a distribution of 50 Mpc/h Lyman-alpha forest effective optical depths that significantly improves agreement with the observations at z~5.6. Extrapolation to z~5.4 and z~5.8 appears promising, but matching the mean background evolution requires evolution in the absorber population beyond the scope of the present model. We also demonstrate the need for extremely large volumes (>400 Mpc on a side) to accurately determine the incidence of rare large-scale features in the Lyman-alpha forest.

The conditions within the ionization front of a quasar during reionization (T ~ 30,000 K, neutral hydrogen fraction ~ 0.5) are ideal for producing Lyman-alpha emission via collisional excitation of hydrogen atoms. Observations of this emission, which could subtend >10 arcmin$^2$ on the sky, would definitively demonstrate the presence of a neutral intergalactic medium at the observed epoch, placing valuable constraints on the progress of reionization. We find that the expected Lyman-alpha surface brightness is significantly weaker than previously determined and may be impossible to observe with current and near-future instruments. Past work calculated the Lyman-alpha emission from a quasar ionization front in a homogeneous medium with a clumping factor approximation to account for inhomogeneities. We find using 1D radiative transfer calculations that this approximation overestimates the emission by a factor of >3. Our calculations model the propagation of ionizing photons and compute the Lyman-alpha emission from quasar ionization fronts on sightlines from a hydrodynamic cosmological simulation at z = 7.1. To better understand the physical properties of the emission, we also develop an analytic model that accurately describes the results of the full radiative transfer calculation.

Interpretation of He II Ly{\alpha} absorption spectra after the epoch of He II reionization requires knowledge of the He II ionizing background. While past work has modelled the evolution of the average background, the standard cosmological radiative transfer technique assumes a uniform radiation field despite the discrete nature of the (rare) bright quasars that dominate the background. We implement a cosmological radiative transfer model that includes the most recent constraints on the ionizing spectra and luminosity function of quasars and the distribution of IGM absorbers. We also estimate, for the first time, the effects of fluctuations on the evolving continuum opacity in two ways: by incorporating the complete distribution of ionizing background amplitudes into the standard approach, and by explicitly treating the quasars as discrete -- but isolated -- sources. Our model results in a He II ionization rate that evolves steeply with redshift, increasing by a factor ~2 from z=3.0 to z=2.5. This causes rapid evolution in the mean He II Ly{\alpha} optical depth -- as recently observed -- without appealing to the reionization of He II. The observed behaviour could instead result from rapid evolution in the mean free path of ionizing photons as the helium in higher H I column density absorbers becomes fully ionized.

Abridged: The giant radio galaxy M 87 with its proximity, famous jet, and very massive black hole provides a unique opportunity to investigate the origin of very high energy (VHE; E>100 GeV) gamma-ray emission generated in relativistic outflows and the surroundings of super-massive black holes. M 87 has been established as a VHE gamma-ray emitter since 2006. The VHE gamma-ray emission displays strong variability on timescales as short as a day. In this paper, results from a joint VHE monitoring campaign on M 87 by the MAGIC and VERITAS instruments in 2010 are reported. During the campaign, a flare at VHE was detected triggering further observations at VHE (H.E.S.S.), X-rays (Chandra), and radio (43 GHz VLBA). The excellent sampling of the VHE gamma-ray light curve enables one to derive a precise temporal characterization of the flare: the single, isolated flare is well described by a two-sided exponential function with significantly different flux rise and decay times. While the overall variability pattern of the 2010 flare appears somewhat different from that of previous VHE flares in 2005 and 2008, they share very similar timescales (~day), peak fluxes (Phi(>0.35 TeV) ~= (1-3) x 10^-11 ph cm^-2 s^-1), and VHE spectra. 43 GHz VLBA radio observations of the inner jet regions indicate no enhanced flux in 2010 in contrast to observations in 2008, where an increase of the radio flux of the innermost core regions coincided with a VHE flare. On the other hand, Chandra X-ray observations taken ~3 days after the peak of the VHE gamma-ray emission reveal an enhanced flux from the core. The long-term (2001-2010) multi-wavelength light curve of M 87, spanning from radio to VHE and including data from HST, LT, VLA and EVN, is used to further investigate the origin of the VHE gamma-ray emission. No unique, common MWL signature of the three VHE flares has been identified.

The TeV radio galaxy M87 is the first radio galaxy detected in the TeV regime. The structure of its jet, which is not pointing towards our line of sight, is spatially resolved in X-ray (by Chandra), optical and radio observations. In 2008, the three main Atmospheric Cherenkov Telescope observatories VERITAS, MAGIC and H.E.S.S. coordinated their observations in a joint campaign from January to May with a total observation time of approx. 120 hours. In February, strong and rapid day-scale TeV flares were detected. VLBA monitoring observations during the same period showed that the 43 GHz radio flux density of the unresolved core began to rise at the time of the TeV flares and eventually reached levels above any previously seen with VLBI. New jet components appeared during the flare. The localization accuracy of the TeV instruments of many arcseconds, even for strong sources, is inadequate to constrain the origin of the emission in the inner jets of AGNs. For M87, with a 6 billion solar mass black hole and a distance of 16.7 Mpc, the VLBA resolution instead corresponds to 30 by 60 Schwarzschild radii. This is starting to resolve the jet collimation region. The temporal coincidence of the TeV and radio flares indicates that they are related and provides the first direct evidence that the TeV radiation from this source is produced within a few tens of RS of the radio core, thought to be coincident to within the VLBA resolution with the black hole.

M87, the central galaxy of the Virgo cluster, is the first radio galaxy detected in the TeV regime. The structure of its jet, which is not pointing toward the line of sight, is spatially resolved in X-ray (by Chandra), in optical and in radio observations. Time correlation between the TeV flux and emission at other wavelengths provides a unique opportunity to localize the very high energy gamma-ray emission process occurring in AGN. For 10 years, M87 has been monitored in the TeV band by atmospheric Cherenkov telescopes. In 2008, the three main atmospheric Cherenkov telescope observatories (H.E.S.S., MAGIC and VERITAS) coordinated their observations in a joint campaign from January to May with a total observation time of approx. 120 hours. The campaign largely overlapped with an intensive VLBA project monitoring the core of M87 at 43 GHz every 5 days. In February, high TeV activities with rapid flares have been detected. Contemporaneously, M87 was observed with high spatial resolution instruments in X-rays (Chandra). We discuss the results of the joint observation campaign in 2008.