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 xH 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 and (posterior medians and 68% credible intervals) when marginalized over quasar lifetimes of 103 ≤ tq ≤ 108 yr.
We decompose the Lyman- (Ly ) forest of an extensive sample of 75 high signal-to-noise ratio and high-resolution quasar spectra into a collection of Voigt profiles. Absorbers located near caustics in ...the peculiar velocity field have the smallest Doppler parameters, resulting in a low-b cutoff in the b- distribution set primarily by the thermal state of the intergalactic medium (IGM). We fit this cutoff as a function of redshift over the range , which allows us to measure the evolution of the IGM temperature-density ( ) relation parameters T0 and γ. We calibrate our measurements against mock Ly forest data generated using 26 hydrodynamic simulations with different thermal histories from the THERMAL suite, also encompassing different values of the IGM pressure smoothing scale. We adopt a forward-modeling approach and self-consistently apply the same algorithms to both data and simulations, propagating both statistical and modeling uncertainties via Monte Carlo. The redshift evolution of T0 (γ) shows a suggestive peak (dip) at z = 2.9 (z = 3). Our measured evolution of T0 and γ is generally in good agreement with previous determinations in the literature. Both the peak in the evolution of T0 at z = 2.8, as well as the high temperatures that we observe at 2.4 < z < 3.4, strongly suggest that a significant episode of heating occurred after the end of H i reionization, which was most likely the cosmic reionization of He ii.
We study quasar proximity zones in the redshift range 5.77 ≤ z ≤ 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 ∝ ( 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 105 years.
The He ii transverse proximity effect-enhanced He ii Ly 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 et al., we present a detailed model of the He ii 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 He ii ultraviolet background model, with the added effect of the radiation from a single bright foreground quasar. We vary the age tage and obscured sky fractions obsc of the foreground quasar, and explore the resulting effect on the He ii transverse proximity effect signal. Fluctuations in intergalactic medium density and the ultraviolet background, as well as the unknown orientation of the foreground quasar, result in a large variance of the He ii Ly transmission along the background sightline. We develop a fully Bayesian statistical formalism to compare far-ultraviolet He ii Ly transmission spectra of the background quasars to our models, and extract joint constraints on tage and obsc for the six Schmidt et al. foreground quasars with the highest implied He ii photoionization rates. Our analysis suggests a bimodal distribution of quasar emission properties, whereby one foreground quasar, associated with a strong He ii transmission spike, is relatively old (22 Myr) and unobscured ( obsc < 35 % ), whereas three others are either younger than 10 Myr or highly obscured ( obsc > 70 % ).
We study the formation of dark matter halos in the concordance Lambda CDM model over a wide range of redshifts, from z = 20 to the present. Our primary focus is the halo mass function, a key probe of ...cosmology. By performing a large suite of nested-box N-body simulations with careful convergence and error controls (60 simulations with box sizes from 4 to 256 h super(-1) Mpc), we determine the mass function and its evolution with excellent statistical and systematic errors, reaching a few percent over most of the considered redshift and mass range. Across the studied redshifts, the halo mass is probed over 6 orders of magnitude (10 super(7)-10 super(13.5) h super(-1) M unk). Historically, there has been considerable variation in the high-redshift mass function as obtained by different groups. We have made a concerted effort to identify and correct possible systematic errors in computing the mass function at high redshift and to explain the discrepancies between some of the previous results. We discuss convergence criteria for the required force resolution, simulation box size, halo mass range, initial and final redshifts, and time stepping. Because of conservative cuts on the mass range probed by individual boxes, our results are relatively insensitive to simulation volume, the remaining sensitivity being consistent with extended Press-Schechter theory. Previously obtained mass function fits near z = 0, when scaled by linear theory, are in good agreement with our results at all redshifts, although a mild redshift dependence consistent with that found by Reed et al. may exist at low redshifts. Overall, our results are consistent with a "universal" form for the mass function at high redshifts.
Many galaxy formation models predict alignments between galaxy spin and the cosmic web (i.e., directions of filaments and sheets), leading to an intrinsic alignment between galaxies that creates a ...systematic error in weak-lensing measurements. These effects are often predicted to be stronger at high redshifts (z 1) that are inaccessible to massive galaxy surveys on foreseeable instrumentation, but IGM tomography of the Ly forest from closely spaced quasars and galaxies is starting to measure the z ∼ 2-3 cosmic web with requisite fidelity. Using mock surveys from hydrodynamical simulations, we examine the utility of this technique, in conjunction with coeval galaxy samples, to measure alignment between galaxies and the cosmic web at z ∼ 2.5. We show that IGM tomography surveys with 5 h−1 Mpc sightline spacing can accurately recover the eigenvectors of the tidal tensor, which we use to define the directions of the cosmic web. For galaxy spins and shapes, we use a model parameterized by the alignment strength, , with respect to the tidal tensor eigenvectors from the underlying density field, and also consider observational effects such as errors in the galaxy position angle, inclination, and redshift. Measurements using the upcoming ∼1 deg2 CLAMATO tomographic survey and 600 coeval zCOSMOS-Deep galaxies should place 3 limits on extreme alignment models with , but much larger surveys encompassing >10,000 galaxies, such as Subaru PFS, will be required to constrain models with . These measurements will constrain models of galaxy-cosmic web alignment and test tidal torque theory at z ∼ 2, improving our understanding of the physics of intrinsic alignments.
Abstract Dense gas in minihalos with masses of 10 6 −10 8 M ⊙ can shield themselves from reionization for ∼100 Myr after being exposed to the UV background. These self-shielded systems, often ...unresolved in cosmological simulations, can introduce strong absorption in quasar spectra. This paper is the first systematic study on the impact of these systems on the Ly α forest. We first derive the H i column density profile of photoevaporating minihalos by conducting 1D radiation–hydrodynamics simulations. We utilize these results to estimate the Ly α opacity from minihalos in a large-scale simulation that cannot resolve self-shielding. When the ionization rate of the background radiation is 0.03 × 10 −12 s −1 , as expected near the end of reionization at z ∼ 5.5, we find that the incidence rate of damped Ly α absorbers increases by a factor of ∼2−4 compared to at z = 4.5. The Ly α flux is, on average, suppressed by ∼3% of its mean due to minihalos. The absorption features enhance the 1D power spectrum up to ∼5% at k ∼ 0.1 h Mpc −1 (or 10 −3 km −1 s), which is comparable to the enhancement caused by inhomogeneous reionization. The flux is particularly suppressed in the vicinity of large halos along the line-of-sight direction at separations of up to 10 h −1 Mpc at r ⊥ ≲ 2 h −1 Mpc. However, these effects become much smaller for higher ionizing rates (≳0.3 × 10 −12 s −1 ) expected in the post-reionization Universe. Our findings highlight the need to consider minihalo absorption when interpreting the Ly α forest at z ≳ 5.5. Moreover, the sensitivity of these quantities to the ionizing background intensity can be exploited to constrain the intensity itself.
The ultraviolet background (UVB) emitted by quasars and galaxies governs the ionization and thermal state of the intergalactic medium (IGM), regulates the formation of high-redshift galaxies, and is ...thus a key quantity for modeling cosmic reionization. The vast majority of cosmological hydrodynamical simulations implement the UVB via a set of spatially uniform photoionization and photoheating rates derived from UVB synthesis models. We show that simulations using canonical UVB rates reionize and, perhaps more importantly, spuriously heat the IGM, much earlier ( ) than they should. This problem arises because at , where observational constraints are nonexistent, the UVB amplitude is far too high. We introduce a new methodology to remedy this issue, and we generate self-consistent photoionization and photoheating rates to model any chosen reionization history. Following this approach, we run a suite of hydrodynamical simulations of different reionization scenarios and explore the impact of the timing of reionization and its concomitant heat injection on the thermal state of the IGM. We present a comprehensive study of the pressure smoothing scale of IGM gas, illustrating its dependence on the details of both hydrogen and helium reionization, and argue that it plays a fundamental role in interpreting Ly forest statistics and the thermal evolution of the IGM. The premature IGM heating we have uncovered implies that previous work has likely dramatically overestimated the impact of photoionization feedback on galaxy formation, which sets the minimum halo mass able to form stars at high redshifts. We make our new UVB photoionization and photoheating rates publicly available for use in future simulations.