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 (z ∼15) than they should. This problem arises because at z > 6, 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.

Cosmological hydrodynamic simulations can accurately predict the properties of the intergalactic medium (IGM), but only under the condition of retaining the high spatial resolution necessary to resolve density fluctuations in the IGM. This resolution constraint prohibits simulating large volumes, such as those probed by BOSS and future surveys, like DESI and 4MOST. To overcome this limitation, we present "Iteratively Matched Statistics" (IMS), a novel method to accurately model the Lyα forest with collisionless N-body simulations, where the relevant density fluctuations are unresolved. We use a small-box, high-resolution hydrodynamic simulation to obtain the probability distribution function (PDF) and the power spectrum of the real-space Lyα forest flux. These two statistics are iteratively mapped onto a pseudo-flux field of an N-body simulation, which we construct from the matter density. We demonstrate that our method can reproduce the PDF, line of sight and 3D power spectra of the Lyα forest with good accuracy (7%, 4%, and 7% respectively). We quantify the performance of the commonly used Gaussian smoothing technique and show that it has significantly lower accuracy (20%-80%), especially for N-body simulations with achievable mean inter-particle separations in large-volume simulations. In addition, we show that IMS produces reasonable and smooth spectra, making it a powerful tool for modeling the IGM in large cosmological volumes and for producing realistic "mock" skies for Lyα forest surveys.

Galaxy formation depends critically on the physical state of gas in the circumgalactic medium (CGM) and its interface with the intergalactic medium (IGM), determined by the complex interplay between inflow from the IGM and outflows from supernovae and/or AGN feedback. The average Lyα absorption profile around galactic halos represents a powerful tool to probe their gaseous environments. We compare predictions from Illustris and Nyx hydrodynamical simulations with the observed absorption around foreground quasars, damped Lyα systems, and Lyman-break galaxies. We show how large-scale BOSS and small-scale quasar pair measurements can be combined to precisely constrain the absorption profile over three decades in transverse distance . Far from galaxies, , the simulations converge to the same profile and provide a reasonable match to the observations. This asymptotic agreement arises because the ΛCDM model successfully describes the ambient IGM and represents a critical advantage of studying the mean absorption profile. However, significant differences between the simulations, and between simulations and observations, are present on scales , illustrating the challenges of accurately modeling and resolving galaxy formation physics. It is noteworthy that these differences are observed as far out as , indicating that the "sphere of influence" of galaxies could extend to approximately ∼7 times the halo virial radius. Current observations are very precise on these scales and can thus strongly discriminate between different galaxy formation models. We demonstrate that the Lyα absorption profile is primarily sensitive to the underlying temperature-density relationship of diffuse gas around galaxies, and argue that it thus provides a fundamental test of galaxy formation models.

We determine the thermal evolution of the intergalactic medium (IGM) over 3 Gyr of cosmic time 1.8 < z < 5.4 by comparing measurements of the Ly forest power spectrum to a suite of 70 hydrodynamical simulations. We conduct Bayesian inference of IGM thermal parameters using an end-to-end forward modeling framework whereby mock spectra generated from our simulation grid are used to build a custom emulator that interpolates the power spectrum between thermal grid points. The temperature at mean density T0 rises steadily from T0 ∼ 6000 K at z=5.4, peaks at 14,000 K for z3.4, and decreases at lower redshift, reaching T07000 K by z1.8. This evolution provides conclusive evidence for photoionization heating resulting from the reionization of He II, as well as the subsequent cooling of the IGM due to the expansion of the universe after all reionization events are complete. Our results are broadly consistent with previous measurements of thermal evolution based on a variety of approaches, but the sensitivity of the power spectrum, the combination of high-precision measurements of largescale modes (k 0.02 s km-1) from the Baryon Oscillation Spectroscopic Survey with our recent determination of the small-scale power, our large grid of models, and our careful statistical analysis allow us to break the wellknown degeneracy between the temperature at mean density T0 and the slope of the temperature-density relation that has plagued previous analyses. At the highest redshifts, z5, we infer lower temperatures than expected from the standard picture of IGM thermal evolution leaving little room for additional smoothing of the Ly forest by free streaming of warm dark matter.

We use N-body simulations to show that high-redshift galaxy counts provide an interesting constraint on the nature of darkmatter, specificallywarm darkmatter (WDM), owing to the lack of early structure formation these models. Our simulations include three WDM models with thermal-production masses of 0.8, 1.3 and 2.6 keV, as well as CDM. Assuming a relationship between dark halo mass and galaxy luminosity that is set by the observed luminosity function at bright magnitudes, we find that 0.8 keV WDM is disfavoured by direct galaxy counts in the Hubble Ultra Deep Field at >10σ. Similarly, 1.3 keV WDM is statistically inconsistent at 2.2σ. Future observations with JWST (and possibly HST via the Frontier Fields) could rule out 1.3 keV WDMat high significance, and may be sensitive to WDM masses greater than 2.6 keV. We also examine the ability of galaxies in these WDM models to reionize the universe, and find that 0.8 and 1.3 keV WDM produce optical depths to the cosmic microwave background (CMB) that are inconsistent at 68 per cent C.L. with current Planck results, even with extremely high ionizing radiation escape fractions, and 2.6 keV WDM requires an optimistic escape fraction to yield an optical depth consistent with Planck data. Although CMB optical depth calculations are model dependent, we find a strong challenge for stellar processes alone to reionize the universe in a 0.8 and 1.3 keV WDM cosmology. © 2014 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society.

In this work we present a new hybrid method to simulate the thermal effects of 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 in the computational cost, thereby allowing for efficient exploration of the parameter space. The method uses a small set of phenomenological input parameters, and combines a seminumerical reionization model to solve for the topology of reionization with an approximate model of how reionization heats the IGM, using the massively parallel Nyx hydrodynamics code which is specifically designed to solve for the structure of diffuse IGM gas. We have produced several medium-scale, high-resolution simulations (20483, Lbox = 40 Mpc h-1) with various instantaneous and inhomogeneous ${\rm H\,{\small I}}$ 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. By analysing the 1D flux power spectrum of these models, we find up to ${\sim } 50{{\rm per\cent}}$ differences in the large-scale properties (low modes, k ≲0.01 s km-1) of the post-reionization power spectrum as a result of the thermal fluctuations. We show that these differences could allow one to distinguish between different reionization scenarios with existing Lyα forest measurements. Finally, we explore the differences in the small-scale cut-off of the power spectrum, finding that, for the same heat input, models show very good agreement provided that the reionization redshift of the instantaneous reionization model occurs at the midpoint of the inhomogeneous model.

The intense radiation emitted by luminous quasars dramatically alters the ionization state of their surrounding IGM. This so-called proximity effect extends out to tens of Mpc, and manifests as large coherent regions of enhanced Ly (Ly) forest transmission in absorption spectra of background sightlines. Here, we present a novel method based on Ly forest tomography, which is capable of mapping these quasar "light echoes" in three dimensions. Using a dense grid (10-100) of faint () background galaxies as absorption probes, one can measure the ionization state of the IGM in the vicinity of a foreground quasar, yielding detailed information about the quasar's radiative history and emission geometry. An end-to-end analysis-combining cosmological hydrodynamical simulations post-processed with a quasar emission model, realistic estimates of galaxy number densities, and instrument + telescope throughput-is conducted to explore the feasibility of detecting quasar light echoes. We present a new, fully Bayesian statistical method that allows one to reconstruct quasar light echoes from thousands of individual low-S/N transmission measurements. Armed with this tool, we undertake an exhaustive parameter study and show that light echoes can be convincingly detected for luminous (M 1450 <-27.5 mag, corresponding to m 1450 < 18.4 mag at ) quasars at redshifts 3 < z QSO < 5, and that a relative precision better than 20% on the quasar age can be achieved for individual objects in the expected range of ages between 1 and 100 Myr. The observational requirements are relatively modest: Moderate-resolution (R 750), multiobject spectroscopy at a low signal-to-noise ratio (S/N > 5) is sufficient, requiring three-hour integrations using existing instruments on 8 m class telescopes.

We present a new measurement of the Lyα forest power spectrum at 1.8 < z < 3.4 using 74 Keck/HIRES and VLT/UVES high-resolution, high-signal-to-noise-ratio quasar spectra. We developed a custom pipeline to measure the power spectrum and its uncertainty, which fully accounts for finite resolution and noise and corrects for the bias induced by masking missing data, damped Lyα absorption systems, and metal absorption lines. Our measurement results in unprecedented precision on the small-scale modes k > 0.02 s km-1, inaccessible to previous SDSS/BOSS analyses. It is well known that these high-k modes are highly sensitive to the thermal state of the intergalactic medium, but contamination by narrow metal lines is a significant concern. We quantify the effect of metals on the small-scale power and find a modest effect on modes with k < 0.1 s km-1. As a result, by masking metals and restricting to k < 0.1 s km-1, their impact is completely mitigated. We present an end-to-end Bayesian forward-modeling framework whereby mock spectra with the same noise, resolution, and masking as our data are generated from Lyα forest simulations. These mock spectra are used to build a custom emulator, enabling us to interpolate between a sparse grid of models and perform Markov chain Monte Carlo fits. Our results agree well with BOSS on scales < 0.02 s km-1, where the measurements overlap. The combination of the percent-level low-k precision of BOSS with our 5%-15% high-k measurements results in a powerful new data set for precisely constraining the thermal history of the intergalactic medium, cosmological parameters, and the nature of dark matter. The power spectra and their covariance matrices are provided as electronic tables.

We show that the canonical oscillation-based (nonresonant) production of sterile neutrino dark matter is inconsistent at >99% confidence with observations of galaxies in the Local Group. We set lower limits on the nonresonant sterile neutrino mass of 2.5 keV (equivalent to 0.7 keV thermal mass) using phase-space densities derived for dwarf satellite galaxies of the Milky Way as well as limits of 8.8 keV (equivalent to 1.8 keV thermal mass) based on subhalo counts of N-body simulations of M 31 analogs. Combined with improved upper mass limits derived from significantly deeper x-ray data of M 31 with full consideration for background variations, we show that there remains little room for nonresonant production if sterile neutrinos are to explain 100% of the dark matter abundance. Resonant and nonoscillation sterile neutrino production remain viable mechanisms for generating sufficient dark matter sterile neutrinos. © 2014 American Physical Society.

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 T 0 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 T 0 (γ) shows a suggestive peak (dip) at z = 2.9 (z = 3). Our measured evolution of T 0 and γ is generally in good agreement with previous determinations in the literature. Both the peak in the evolution of T 0 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.