BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

James S. Bolton

James S. Bolton appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-ph.GA

Catalog footprint

What is connected

32works
2topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2022arXiv

Limits on non-canonical heating and turbulence in the intergalactic medium from the low redshift Lyman-alpha forest

We examine the column density distribution function (CDDF) and Doppler parameter distribution from hydrodynamical simulations and Cosmic Origins Spectrograph (COS) observations of the Lyman-alpha forest at redshift $0\leq z\leq 0.2$. Allowing for a factor of two uncertainty in the metagalactic HI photoionisation rate, our hydrodynamical simulations are in good agreement ($1$-$1.5σ$) with the shape and amplitude of the observed CDDF at HI column densities $10^{13.3}\rm\,cm^{-2}\leq N_{\rm HI}\leq 10^{14.5}\rm\,cm^{-2}$. However, the Doppler widths of the simulated lines remain too narrow with respect to the COS data. We argue that invoking AGN feedback does not resolve this discrepancy. We also disfavour enhanced photoheating rates as a potential solution, as this requires an unphysically hard UV background spectrum. If instead appealing to a non-canonical source of heating, an additional specific heat injection of $u \lesssim 6.9\rm\,eV\,m_{\rm p}^{-1}$ is required at $z\lesssim 2.5$ for gas that has $N_{\rm HI}\simeq 10^{13.5}\rm\,cm^{-2}$ by $z=0.1$. Alternatively, there may be an unresolved line of sight turbulent velocity component of $v_{\rm turb}\lesssim 8.5\rm\,km\,s^{-1}(N_{\rm HI}/10^{13.5}\rm\,cm^{-2})^{0.21}$ for the coldest gas in the diffuse IGM.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The effect of inhomogeneous reionisation on the Lyman-$α$ forest power spectrum at redshift $z>4$: implications for thermal parameter recovery

We use the Sherwood-Relics suite of hybrid hydrodynamical and radiative transfer simulations to model the effect of inhomogeneous reionisation on the 1D power spectrum of the \Lya forest transmitted flux at redshifts $4.2\leq z \leq 5$. Relative to models that assume a homogeneous UV background, reionisation suppresses the power spectrum at small scales, $k \sim 0.1\rm\,km^{-1}\,s$ by $\sim 10$ per cent because of spatial variations in the thermal broadening kernel and the divergent peculiar velocity field associated with over-pressurised intergalactic gas. On larger scales, $k<0.03\rm\,km^{-1}\,s$, the power spectrum is instead enhanced by $10$--$50$ per cent by large scale spatial variations in the neutral hydrogen fraction. The effect of inhomogeneous reionisation must therefore be accounted for in analyses of forthcoming high precision measurements. We provide a correction for the Ly$α$ forest power spectrum at $4.1\leq z \leq 5.4$ that can be easily applied within other parameter inference frameworks using similar reionisation models. We perform a Bayesian analysis of mock data to assess the extent of systematic biases that may arise in measurements of the intergalactic medium if ignoring this correction. At the scales probed by current high resolution Ly$α$ forest data at $z>4$, $0.006 \rm \,km^{-1}\,s\leq k \leq 0.2 \rm\, km^{-1}\,s$, we find inhomogeneous reionisation does not introduce any significant bias in thermal parameter recovery for the current measurement uncertainties of $\sim 10$ per cent. However, for $5$ per cent uncertainties, $\sim 1σ$ shifts between the estimated and true parameters occur.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Probing the thermal state of the intergalactic medium at $z>5$ with the transmission spikes in high-resolution Ly$α$ forest spectra

We compare a sample of five high-resolution, high S/N Ly$α$ forest spectra of bright $6<z \lesssim 6.5$ QSOs aimed at spectrally resolving the last remaining transmission spikes at $z>5$ with those obtained from mock absorption spectra from the Sherwood and Sherwood-Relics suites of hydrodynamical simulations of the intergalactic medium (IGM). We use a profile fitting procedure for the inverted transmitted flux, $1-F$, similar to the widely used Voigt profile fitting of the transmitted flux $F$ at lower redshifts, to characterise the transmission spikes that probe predominately underdense regions of the IGM. We are able to reproduce the width and height distributions of the transmission spikes, both with optically thin simulations of the post-reionization Universe using a homogeneous UV background and full radiative transfer simulations of a late reionization model. We find that the width of the fitted components of the simulated transmission spikes is very sensitive to the instantaneous temperature of the reionized IGM. The internal structures of the spikes are more prominant in low temeperature models of the IGM. The width distribution of the observed transmission spikes, which require high spectral resolution ($\leq $ 8 km/s) to be resolved, is reproduced for optically thin simulations with a temperature at mean density of $T_0= (11000 \pm 1600,10500\pm 2100,12000 \pm 2200)$ K at $z= (5.4,5.6,5.8)$. This is weakly dependent on the slope of the temperature-density relation, which is favoured to be moderately steeper than isothermal. In the inhomogeneous, late reionization, full radiative transfer simulations where islands of neutral hydrogen persist to $z\sim5.3$, the width distribution of the observed transmission spikes is consistent with the range of $T_0$ caused by spatial fluctuations in the temperature-density relation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Inferring the IGM thermal history during reionisation with the Lyman-$α$ forest power spectrum at redshift $z \simeq 5$

We use cosmological hydrodynamical simulations to assess the feasibility of constraining the thermal history of the intergalactic medium during reionisation with the Ly$α$ forest at $z\simeq5$. The integrated thermal history has a measureable impact on the transmitted flux power spectrum that can be isolated from Doppler broadening at this redshift. We parameterise this using the cumulative energy per proton, $u_0$, deposited into a gas parcel at the mean background density, a quantity that is tightly linked with the gas density power spectrum in the simulations. We construct mock observations of the line of sight Ly$α$ forest power spectrum and use a Markov Chain Monte Carlo approach to recover $u_{0}$ at redshifts $5 \leq z \leq 12$. A statistical uncertainty of $\sim 20$ per cent is expected (at 68 per cent confidence) at $z\simeq 5$ using high resolution spectra with a total redshift path length of $Δz=4$ and a typical signal-to-noise ratio of $\rm S/N=15$ per pixel. Estimates for the expected systematic uncertainties are comparable, such that existing data should enable a measurement of $u_0$ to within $\sim 30$ per cent. This translates to distinguishing between reionisation scenarios with similar instantaneous temperatures at $z\simeq 5$, but with an energy deposited per proton that differs by $2$-$3\, \rm eV$ over the redshift interval $5\leq z \leq 12$. For an initial temperature of $T\sim 10^{4}\rm\,K$ following reionisation, this corresponds to the difference between early ($z_{\rm re}=12$) and late ($z_{\rm re}=7$) reionisation in our models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Testing the effect of galactic feedback on the IGM at z ~ 6 with metal-line absorbers

We present models of low- and high-ionization metal-line absorbers (O I, C II, C IV and Mg II) during the end of the reionization epoch, at z ~ 6. Using four cosmological hydrodynamical simulations with different feedback schemes (including the Illustris and Sherwood simulations) and two different choices of hydro-solver, we investigate how the overall incidence rate and equivalent width distribution of metal-line absorbers varies with the galactic wind prescription. We find that the O I and C II absorbers are reasonably insensitive to the feedback scheme. All models, however, struggle to reproduce the observations of C IV and Mg II, which are probing down to lower overdensities than O I and C II at z ~ 6, suggesting that the metals in the simulations are not being transported out into the IGM efficiently enough. The situation is improved but not resolved if we choose a harder (but still reasonable) and/or (locally) increased UV background at z ~ 6.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Sherwood simulation suite: overview and data comparisons with the Lyman-alpha forest at redshifts 2 < z < 5

We introduce a new set of large scale, high resolution hydrodynamical simulations of the intergalactic medium: the Sherwood simulation suite. These are performed in volumes 10^3-160^3 h^-3 comoving Mpc^{3}, span almost four orders of magnitude in mass resolution with up to 17.2 billion particles, and employ a variety of physics variations including warm dark matter and galactic outflows. We undertake a detailed comparison of the simulations to high resolution, high signal-to-noise observations of the Lyman-alpha forest over the redshift range 2 < z < 5. The simulations are in very good agreement with the observational data, lending further support to the paradigm that the Lyman-alpha forest is a natural consequence of the web-like distribution of matter arising in LCDM cosmological models. Only a small number of minor discrepancies remain with respect to the observational data. Saturated Lyman-alpha absorption lines with column densities N_HI > 10^14.5 cm^-2 at 2 < z < 2.5 are underpredicted in the models. An uncertain correction for continuum placement bias is required to match the distribution and power spectrum of the transmitted flux, particularly at z > 4. Finally, the temperature of intergalactic gas in the simulations may be slightly too low at z=2.7 and a flatter temperature-density relation is required at z=2.4, consistent with the expected effects of non-equilibrium ionisation during He-II reionisation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Constraining the temperature-density relation of the intergalactic medium with the Lyman-$α$ and $β$ forests

The post-reionization thermal state of the intergalactic medium is characterized by a power-law relationship between temperature and density, with a slope determined by the parameter $γ$. We describe a new method to measure $γ$ using the ratio of flux curvature in the Lyman-$α$ and $β$ forests. At a given redshift, this curvature ratio incorporates information from the different gas densities traced by Lyman-$α$ and $β$ absorption. It is relatively simple and fast to compute and appears robust against several observational uncertainties. We apply this technique to a sample of 27 high-resolution quasar spectra from the Very Large Telescope. While promising statistical errors on $γ$ appear to be achievable with these spectra, to reach its full potential, the dependence of the curvature ratio on the thermal state of the gas in the foreground Lyman-$α$ forest will require further, detailed forward modelling.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Evidence of patchy hydrogen reionization from an extreme Ly$α$ trough below redshift six

We report the discovery of an extremely long ($\sim$110 Mpc/$h$) and dark ($τ_{\rm eff} \gtrsim 7$) Ly$α$ trough extending down to $z \simeq 5.5$ towards the $z_{\rm em} \simeq 6.0$ quasar ULAS J0148+0600. We use these new data in combination with Ly$α$ forest measurements from 42 quasars at $4.5 \le z_{\rm em} \le 6.4$ to conduct an updated analysis of the line-of-sight variance in the intergalactic Ly$α$ opacity over $4 \le z \le 6$. We find that the scatter in transmission among lines of sight near $z \sim 6$ significantly exceeds theoretical expectations for either a uniform ultraviolet background (UVB) or simple fluctuating UVB models in which the mean free path to ionizing photons is spatially invariant. The data, particularly near $z \simeq 5.6$-5.8, instead require fluctuations in the volume-weighted hydrogen neutral fraction that are a factor of 3 or more beyond those expected from density variations alone. We argue that these fluctuations are most likely driven by large-scale variations in the mean free path, consistent with expectations for the final stages of inhomogeneous hydrogen reionization. Even by $z \simeq 5.6$, however, a large fraction of the data are consistent with a uniform UVB, and by $z \sim 5$ the data are fully consistent with opacity fluctuations arising solely from the density field. This suggests that while reionization may be ongoing at $z \sim 6$, it has fully completed by $z \sim 5$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lyman-alpha emitters gone missing: evidence for late reionization?

We combine high resolution hydrodynamical simulations with an intermediate resolution, dark matter only simulation and an analytical model for the growth of ionized regions to estimate the large scale distribution and redshift evolution of the visibility of Lyman-alpha emission in 6<=z<=8 galaxies. The inhomogeneous distribution of neutral hydrogen during the reionization process results in significant fluctuations in the Lyman-alpha transmissivity on large scales. The transmissivity depends not only on the ionized fraction of the intergalactic medium by volume and the amplitude of the local ionizing background, but is also rather sensitive to the evolution of the relative velocity shift of the Lyman-alpha emission line due to resonant scattering. We reproduce a decline in the space density of Lyman-alpha emitting galaxies as rapid as observed with a rather rapidly evolving neutral fraction between z=6-8, and a typical Lyman-alpha line velocity offset of 100 km/s redward of systemic at z=6 which decreases toward higher redshift. The new (02/2015) Planck results indicate such a recent end to reionization is no longer disfavoured by constraints from the cosmic microwave background.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Reionization and high-redshift galaxies: the view from quasar absorption lines

Determining when and how the first galaxies reionized the intergalactic medium (IGM) promises to shed light on both the nature of the first objects and the cosmic history of baryons. Towards this goal, quasar absorption lines play a unique role by probing the properties of diffuse gas on galactic and intergalactic scales. In this review we examine the multiple ways in which absorption lines trace the connection between galaxies and the IGM near the reionization epoch. We first describe how the Ly$α$ forest is used to determine the intensity of the ionizing ultraviolet background and the global ionizing emissivity budget. Critically, these measurements reflect the escaping ionizing radiation from all galaxies, including those too faint to detect directly. We then discuss insights from metal absorption lines into reionization-era galaxies and their surroundings. Current observations suggest a buildup of metals in the circumgalactic environments of galaxies over $z \sim 6$ to 5, although changes in ionization will also affect the evolution of metal line properties. A substantial fraction of metal absorbers at these redshifts may trace relatively low-mass galaxies. Finally, we review constraints from the Ly$α$ forest and quasar near zones on the timing of reionization. Along with other probes of the high-redshift Universe, absorption line data are consistent with a relatively late end to reionization ($5.5 \lesssim z \lesssim 7$); however the constraints are still fairly weak. Significant progress is expected to come through improved analysis techniques, increases in the number of known high-redshift quasars from optical and infrared sky surveys, large gains in sensitivity from next-generation observing facilities, and synergies with other probes of the reionization era.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The impact of temperature fluctuations on the large-scale clustering of the Ly$α$ forest

We develop a semi-analytic method for assessing the impact of the large-scale IGM temperature fluctuations expected following He${\rm\,{\scriptstyle II}}$ reionization on three-dimensional clustering measurements of the Ly$α$ forest. Our methodology builds upon the existing large volume, mock Ly$α$ forest survey simulations presented by Greig et al. by including a prescription for a spatially inhomogeneous ionizing background, temperature fluctuations induced by patchy He${\rm\,{\scriptstyle II}}$ photoheating and the clustering of quasars. This approach enables us to achieve a dynamic range within our semi-analytic model substantially larger than currently feasible with computationally expensive, fully numerical simulations. The results agree well with existing numerical simulations, with large-scale temperature fluctuations introducing a scale-dependent increase in the spherically averaged 3D Ly$α$ forest power spectrum of up to 20-30 per cent at wavenumbers $k\sim0.02$ Mpc$^{-1}$. Although these large-scale thermal fluctuations will not substantially impact upon the recovery of the baryon acoustic oscillation scale from existing and forthcoming dark energy spectroscopic surveys, any complete forward modelling of the broad-band term in the Ly$α$ correlation function will none the less require their inclusion.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The photoheating of the intergalactic medium in synthesis models of the UV background

We compare cosmological hydrodynamical simulations combined with the homogeneous metagalactic UV background (UVB) of Haardt & Madau (2012) (HM2012) to observations of the Lyman-alpha forest that are sensitive to the thermal and ionization state of the intergalactic medium (IGM). The transition from optically thick to thin photoheating predicted by the simple one-zone, radiative transfer model implemented by HM2012 predicts a thermal history that is in remarkably good agreement with the observed rise of the IGM temperature at z~3 if we account for the expected evolution of the volume filling factor of HeIII. Our simulations indicate that there may be, however, some tension between the observed peak in the temperature evolution and the rather slow evolution of the HeII opacities suggested by recent Hubble Space Telescope/COS measurements. The HM2012 UVB also underpredicts the metagalactic hydrogen photoionization rate required by our simulations to match the observed opacity of the forest at z>4 and z<2.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gas around galaxy haloes: methodology comparisons using hydrodynamical simulations of the intergalactic medium

We perform cosmological simulations of the intergalactic medium (IGM) at redshift z ~ 3 using the numerical gravity-hydrodynamics codes GADGET-3 and Enzo for the purpose of modelling the gaseous environments of galaxies. We identify haloes in the simulations using three different algorithms. Different rank orderings of the haloes by mass result, introducing a limiting factor in identifying haloes with observed galaxies. We also compare the physical properties of the gas between the two codes, focussing primarily on the gas outside the virial radius, motivated by recent HI absorption measurements of the gas around z ~ 2 - 3 galaxies. The internal dispersion velocities of the gas in the haloes have converged for a box size of 30 comoving Mpc, but the centre-of-mass peculiar velocities of the haloes have not up to a box size of 60 comoving Mpc. The density and temperature of the gas within the instantaneous turn-around radii of the haloes are adequately captured for box sizes 30 Mpc on a side, but the results are highly sensitive to the treatment of unresolved, rapidly cooling gas, with the gas mass fraction within the virial radius severely depleted by star formation in the GADGET-3 simulations. Convergence of the gas peculiar velocity field on large scales requires a box size of at least 60 Mpc. Outside the turn-around radius, the physical state of the gas agrees to 30 percent or better both with box size and between simulation methods. We conclude that generic IGM simulations make accurate predictions for the intergalactic gas properties beyond the halo turn-around radii, but the gas properties on smaller scales are highly dependent on star formation and feedback implementations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

IGM Constraints from the SDSS-III/BOSS DR9 Ly-alpha Forest Flux Probability Distribution Function

The Ly$α$ forest transmission probability distribution function (PDF) is an established probe of the intergalactic medium (IGM) astrophysics, especially the temperature-density relationship of the IGM. We measure the transmission PDF from 3393 Baryon Oscillations Spectroscopic Survey (BOSS) quasars from SDSS Data Release 9, and compare with mock spectra that include careful modeling of the noise, continuum, and astrophysical uncertainties. The BOSS transmission PDFs, measured at $\langle z \rangle = [2.3,2.6,3.0]$, are compared with PDFs created from mock spectra drawn from a suite of hydrodynamical simulations that sample the IGM temperature-density relationship, $γ$, and temperature at mean-density, $T_0$, where $T(Δ) = T_0 Δ^{γ-1}$. We find that a significant population of partial Lyman-limit systems with a column-density distribution slope of $β_\mathrm{pLLS} \sim -2$ are required to explain the data at the low-transmission end of transmission PDF, while uncertainties in the mean Ly$α$ forest transmission affect the high-transmission end. After modelling the LLSs and marginalizing over mean-transmission uncertainties, we find that $γ=1.6$ best describes the data over our entire redshift range, although constraints on $T_0$ are affected by systematic uncertainties. Within our model framework, isothermal or inverted temperature-density relationships ($γ\leq 1$) are disfavored at a significance of over 4$σ$, although this could be somewhat weakened by cosmological and astrophysical uncertainties that we did not model.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Suite of Hydrodynamical Simulations for the Lyman-Alpha Forest with Massive Neutrinos

The signature left in quasar spectra by the presence of neutral hydrogen in the Universe allows one to constrain the sum of the neutrino masses with improved sensitivity, with respect to laboratory experiments, and may shed a new light on the neutrino mass hierarchy and on the absolute mass scale of neutrinos. Constraints on cosmological parameters and on the dark energy equation of state can also be derived, from a joint parameter estimation procedure. However, this requires a detailed modeling of the line-of-sight power spectrum of the transmitted flux in the Lyman-Alpha (LyA) forest on scales ranging from a few to hundreds of Mpcs, which in turns demands the inclusion and careful treatment of cosmological neutrinos. To this end, we present here a suite of state-of-the-art hydrodynamical simulations with cold dark matter, baryons and massive neutrinos, specifically targeted for modeling the low-density regions of the IGM as probed by the LyA forest at high-redshift. The simulations span volumes ranging from (25 Mpc/h)^3 to (100 Mpc/h)^3, and are made using either 3 X 192^3~21 millions or 3 X 768^3~1.4 billion particles. The resolution of the various runs can be further enhanced, so that we can reach the equivalent of 3 X 3072^3~87 billion particles in a (100 Mpc/h)^3 box size. The chosen cosmological parameters are compatible with the latest Planck (2013) results, although we also explore the effect of slight variations in the main cosmological and astrophysical parameters. We adopt a particle-type implementation of massive neutrinos, and consider three degenerate species having masses M_nu =0.1, 0.2, 0.3, 0.4 and 0.8 eV, respectively. We improve on previous studies in several ways, in particular with updated routines for IGM radiative cooling and heating processes, and initial conditions based on 2LPT rather than the Zeldovich approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The thermal history of the intergalactic medium down to redshift z=1.5: a new curvature measurement

According to the photo-heating model of the intergalactic medium (IGM), He II reionization is expected to affect its thermal evolution. Evidence for additional energy injection into the IGM has been found at $3\lesssim z\lesssim4$, though the evidence for the subsequent fall-off below $z\sim2.8$ is weaker and depends on the slope of the temperature--density relation, $γ$. Here we present, for the first time, an extension of the IGM temperature measurements down to the atmospheric cut-off of the H I Lyman-$α$ forest at $z\simeq1.5$. Applying the curvature method on a sample of 60 UVES spectra we investigated the thermal history of the IGM at $z<3$ with precision comparable to the higher redshift results. We find that the temperature of the cosmic gas traced by the Ly-$α$ forest [$T(\barΔ)]$ increases for increasing overdensity from $T(\barΔ)\sim 22670$ K to 33740 K in the redshift range $z\sim2.8-1.6$. Under the assumption of two reasonable values for $γ$, the temperature at the mean density ($T_{0}$) shows a tendency to flatten at $z\lesssim 2.8$. In the case of $γ\sim1.5$, our results are consistent with previous ones which indicate a falling $T_{0}$ for redshifts $z\lesssim2.8$. Finally, our $T(\barΔ)$ values show reasonable agreement with moderate blazar heating models.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A consistent determination of the temperature of the intergalactic medium at redshift z=2.4

We present new measurements of the thermal state of the intergalactic medium (IGM) at $z\sim2.4$ derived from absorption line profiles in the Ly$α$ forest. We use a large set of high-resolution hydrodynamical simulations to calibrate the relationship between the temperature-density ($T$--$Δ$) relation in the IGM and the distribution of HI column densities, $N_{\rm HI}$, and velocity widths, $b_{\rm HI}$, of discrete Ly$α$ forest absorbers. This calibration is then applied to the measurement of the lower cut-off of the $b_{\rm HI}$--$N_{\rm HI}$ distribution recently presented by Rudie et al. (2012). We infer a power-law $T$--$Δ$ relation, $T=T_{0}Δ^{γ-1}$, with a temperature at mean density, $T_{0}=[1.00^{+0.32}_{-0.21}]\times10^{4}\rm\,K$ and slope $(γ-1)=0.54\pm0.11$. The slope is fully consistent with that advocated by the analysis of Rudie et al (2012); however, the temperature at mean density is lower by almost a factor of two, primarily due to an adjustment in the relationship between column density and physical density assumed by these authors. These new results are in excellent agreement with the recent temperature measurements of Becker et al. (2011), based on the curvature of the transmitted flux in the Ly$α$ forest. This suggests that the thermal state of the IGM at this redshift is reasonably well characterised over the range of densities probed by these methods.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

New Measurements of the Ionizing Ultraviolet Background over 2 < z < 5 and Implications for Hydrogen Reionization

We present new measurements of the intensity of the ionizing ultraviolet background and the global emissivity of ionizing photons over 2 < z < 5. Our results are based on a suite of updated measurements of physical properties of the high-redshift intergalactic medium (IGM), including gas temperatures and the opacity of the IGM to Ly-alpha and ionizing photons. Consistent with previous works, we find a relatively flat hydrogen photoionization rate over 2 < z < 5, although our measurements are roughly a factor of two higher than the 2008 values of Faucher-Giguere et al., due primarily to our lower gas temperatures. The ionizing emissivity we derive is also generally higher than other recent estimates due to a combination of lower gas temperatures, higher ionizing opacity, and an accounting of cosmological radiative transfer effects. We find evidence that the emissivity increases from z~3 to 5, reaching ~5 ionizing photons per atom per gigayear at z=4.75 for realistic galaxy spectra. We further find that galaxies must dominate the emissivity near 1 Ryd at z > 4, and possibly at all redshifts z > 2.4. Our results suggest that the globally-averaged ionizing "efficiency" of star-forming galaxies increases substantially with redshift over 3.2 < z < 4.75. This trend is consistent with the conclusion often drawn from reionization models that the ionizing efficiency of galaxies must be higher during reionization in order for galaxies to reionize the IGM by z=6. Our emissivity values at z~5 suggest that ionizing photons may have been a factor of two more abundant during the final stages of reionization than previously indicated. The evolution of the ionizing emissivity over 2 < z < 5 suggests, moreover, that the steep decline in the photoionization rate from z~5 to 6 may indicate a rapid evolution in the mean free path at z > 5.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Probing the metallicity and ionization state of the circumgalactic medium at z~6 and beyond with OI absorption

Low ionization metal absorption due to OI has been identified as an important probe of the physical state of the inter-/circumgalactic medium at the tail-end of reionization. We use here high-resolution hydrodynamic simulations to interpret the incidence rate of OI absorbers at z~6 as observed by Becker et al. (2011). We infer weak OI absorbers (EW > 0.1 A) to have typical HI column densities in the range of sub-DLAs, densities of 80 times the mean baryonic density and metallicities of about 1/500 th solar. This is similar to the metallicity inferred at similar overdensities at z~3, suggesting that the metal enrichment of the circumgalactic medium around low-mass galaxies has already progressed considerably by z~6. The apparently rapid evolution of the incidence rates for OI absorption over the redshift range 5<z<6 mirrors that of self-shielded Lyman-Limit systems at lower redshift and is mainly due to the rapid decrease of the meta-galactic photo-ionization rate at z>5. We predict the incidence rate of OI absorbers to continue to rise rapidly with increasing redshift as the IGM becomes more neutral. If the distribution of metals extends to lower density regions, OI absorbers will allow the metal enrichment of the increasingly neutral filamentary structures of the cosmic web to be probed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The one-dimensional Ly-alpha forest power spectrum from BOSS

We have developed two independent methods to measure the one-dimensional power spectrum of the transmitted flux in the Lyman-$α$ forest. The first method is based on a Fourier transform, and the second on a maximum likelihood estimator. The two methods are independent and have different systematic uncertainties. The determination of the noise level in the data spectra was subject to a novel treatment, because of its significant impact on the derived power spectrum. We applied the two methods to 13,821 quasar spectra from SDSS-III/BOSS DR9 selected from a larger sample of over 60,000 spectra on the basis of their high quality, large signal-to-noise ratio, and good spectral resolution. The power spectra measured using either approach are in good agreement over all twelve redshift bins from $<z> = 2.2$ to $<z> = 4.4$, and scales from 0.001 $\rm(km/s)^{-1}$ to $0.02 \rm(km/s)^{-1}$. We determine the methodological and instrumental systematic uncertainties of our measurements. We provide a preliminary cosmological interpretation of our measurements using available hydrodynamical simulations. The improvement in precision over previously published results from SDSS is a factor 2--3 for constraints on relevant cosmological parameters. For a $Λ$CDM model and using a constraint on $H_0$ that encompasses measurements based on the local distance ladder and on CMB anisotropies, we infer $σ_8 =0.83\pm0.03$ and $n_s= 0.97\pm0.02$ based on \ion{H}{i} absorption in the range $2.1<z<3.7$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the rapid demise of Lyman-alpha emitters at z>7 due to the increasing incidence of optically thick absorption systems

A variety of independent observational studies have now reported a significant decline in the fraction of Lyman-break galaxies which exhibit Ly-a emission over the redshift interval z=6-7. In combination with the strong damping wing extending redward of Ly-a in the spectrum of the bright z=7.085 quasar ULAS 1120+0641, this has strengthened suggestions that the hydrogen in the intergalactic medium (IGM) is still substantially neutral at z~7. Current theoretical models imply HI fractions as large as 40-90 per cent may be required to explain these data assuming there is no intrinsic evolution in the Ly-a emitter population. We propose that such large neutral fractions are not necessary. Based on a hydrodynamical simulation which reproduces the absorption spectra of high-redshift (z~6-7) quasars, we demonstrate that the opacity of the intervening IGM redward of rest-frame Ly-a can rise rapidly in average regions of the Universe simply because of the increasing incidence of absorption systems which are optically thick to Lyman continuum photons as the tail-end of reionisation is approached. Our simulations suggest these data do not require a large change in the IGM neutral fraction by several tens of per cent from z=6-7, but may instead be indicative of the rapid decrease in the typical mean free path for ionising photons expected during the final stages of reionisation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Thermal constraints on the reionisation of hydrogen by population-II stellar sources

Measurements of the intergalactic medium (IGM) temperature provide a potentially powerful constraint on the reionisation history due to the thermal imprint left by the photo-ionisation of neutral hydrogen. However, until recently IGM temperature measurements were limited to redshifts 2 < z < 4.8, restricting the ability of these data to probe the reionisation history at z > 6. In this work, we use recent measurements of the IGM temperature in the near-zones of seven quasars at z ~ 5.8 - 6.4, combined with a semi-numerical model for inhomogeneous reionisation, to establish new constraints on the redshift at which hydrogen reionisation completed. We calibrate the model to reproduce observational constraints on the electron scattering optical depth and the HI photo-ionisation rate, and compute the resulting spatially inhomogeneous temperature distribution at z ~ 6 for a variety of reionisation scenarios. Under standard assumptions for the ionising spectra of population-II sources, the near-zone temperature measurements constrain the redshift by which hydrogen reionisation was complete to be z > 7.9 (6.5) at 68 (95) per cent confidence. We conclude that future temperature measurements around other high redshift quasars will significantly increase the power of this technique, enabling these results to be tightened and generalised.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fast, large volume, GPU enabled simulations for the Ly-alpha forest: power spectrum forecasts for baryon acoustic oscillation experiments

High redshift measurements of the baryonic acoustic oscillation scale (BAO) from large Ly-alpha forest surveys represent the next frontier of dark energy studies. As part of this effort, efficient simulations of the BAO signature from the Ly-alpha forest will be required. We construct a model for producing fast, large volume simulations of the Ly-alpha forest for this purpose. Utilising a calibrated semi-analytic approach, we are able to run very large simulations in 1 Gpc^3 volumes which fully resolve the Jeans scale in less than a day on a desktop PC using a GPU enabled version of our code. The Ly-alpha forest spectra extracted from our semi-analytical simulations are in excellent agreement with those obtained from a fully hydrodynamical reference simulation. Furthermore, we find our simulated data are in broad agreement with observational measurements of the flux probability distribution and 1D flux power spectrum. We are able to correctly recover the input BAO scale from the 3D Ly-alpha flux power spectrum measured from our simulated data, and estimate that a BOSS-like 10^4 deg^2 survey with ~15 background sources per square degree and a signal-to-noise of ~5 per pixel should achieve a measurement of the BAO scale to within ~1.4 per cent. We also use our simulations to provide simple power-law expressions for estimating the fractional error on the BAO scale on varying the signal-to-noise and the number density of background sources. The speed and flexibility of our approach is well suited for exploring parameter space and the impact of observational and astrophysical systematics on the recovery of the BAO signature from forthcoming large scale spectroscopic surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

How neutral is the intergalactic medium surrounding the redshift z=7.085 quasar ULAS J1120+0641?

The quasar ULAS J1120+0641 at redshift z=7.085 has a highly ionised near zone which is smaller than those around quasars of similar luminosity at z~6. The spectrum also exhibits evidence for a damping wing extending redward of the systemic Lya redshift. We use radiative transfer simulations in a cosmological context to investigate the implications for the ionisation state of the inhomogeneous IGM surrounding this quasar. Our simulations show that the transmission profile is consistent with an IGM in the vicinity of the quasar with a volume averaged HI fraction of f_HI>0.1 and that ULAS J1120+0641 has been bright for 10^6--10^7 yr. The observed spectrum is also consistent with smaller IGM neutral fractions, f_HI ~ 10^-3--10-4, if a damped Lya system in an otherwise highly ionised IGM lies within 5 proper Mpc of the quasar. This is, however, predicted to occur in only ~5 per cent of our simulated sight-lines for a bright phase of 10^6--10^7 yr. Unless ULAS J1120+0641 grows during a previous optically obscured phase, the low age inferred for the quasar adds to the theoretical challenge of forming a 2x10^9 M_sol black hole at this high redshift.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Improved measurements of the intergalactic medium temperature around quasars: possible evidence for the initial stages of He-II reionisation at z~6

We present measurements of the intergalactic medium (IGM) temperature within ~5 proper Mpc of seven luminous quasars at z~6. The constraints are obtained from the Doppler widths of Lyman-alpha absorption lines in the quasar near-zones and build upon our previous measurement for the z=6.02 quasar SDSS J0818+1722. The expanded data set, combined with an improved treatment of systematic uncertainties, yields an average temperature at the mean density of log (T_0/K) = 4.21 \pm 0.03 (\pm^0.06_0.07) at 68 (95) per cent confidence for a flat prior distribution over 3.2 < log (T_0/K) < 4.8. In comparison, temperatures measured from the general IGM at z~5 are ~0.3 dex cooler, implying an additional source of heating around these quasars which is not yet present in the general IGM at slightly lower redshift. This heating is most likely due to the recent reionisation of He-II in vicinity of these quasars, which have hard, non-thermal ionising spectra. The elevated temperatures may therefore represent evidence for the earliest stages of He-II reionisation in the most biased regions of the high-redshift Universe. The temperature as a function of distance from the quasars is consistent with being constant, log(T_0/K)~4.2, with no evidence for a line-of-sight thermal proximity effect. However, the limited extent of the quasar near-zones prevents the detection of He-III regions larger than ~5 proper Mpc. Under the assumption the quasars have reionised the He-II in their vicinity, we infer the data are consistent with an average optically bright phase of duration in excess of 10^6.5 yr. These measurements represent the highest-redshift IGM temperature constraints to date, and thus provide a valuable data set for confronting models of H-I reionisation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The impact of spatial fluctuations in the ultra-violet background on intergalactic carbon and silicon

Spatial inhomogeneities in the spectral shape of the ultra-violet background (UVB) at the tail-end of HeII reionisation are thought to be the primary cause of the large fluctuations observed in the HeII to HI Ly-a forest optical depth ratio, tau_HeII/tau_HI, at z~2-3. These spectral hardness fluctuations will also influence the ionisation balance of intergalactic metals; we extract realistic quasar absorption spectra from a large hydrodynamical simulation to examine their impact on intergalactic SiIV and CIV absorbers. Using a variety of toy UVB models, we find that while the predicted spatial inhomogeneities in spectral hardness have a significant impact on tau_HeII/tau_HI, the longer mean free path for photons with frequencies above and below the HeII ionisation edge means these fluctuations have less effect on the SiIV and CIV ionisation balance. Furthermore, UVB models which produce the largest fluctuations in specific intensity at the HeII ionisation edge also have the softest ionising spectra, and thus result in photo-ionisation rates which are too low to produce significant fluctuations in the observed tau_SiIV/tau_CIV. Instead, we find spatial variations in the IGM metallicity will dominate any scatter in tau_SiIV/tau_CIV. Our results suggest that observational evidence for homogeneity in the observed tau_SiIV/tau_CIV distribution does not rule out the possibility of significant fluctuations in the UVB spectral shape at z~2-3. On the other hand, the scatter in metallicity inferred from observations of intergalactic CIV and SiIV absorption at z~2-3 using spatially uniform ionisation corrections is likely intrinsic, and therefore provides a valuable constraint on intergalactic metal enrichment scenarios at these redshifts.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A first direct measurement of the intergalactic medium temperature around a quasar at z=6

The thermal state of the intergalactic medium (IGM) provides an indirect probe of both the HI and HeII reionisation epochs. Current constraints on the IGM temperature from the Lya forest are restricted to the redshift range 2<z<4.5, limiting the ability to probe the thermal memory of HI reionisation toward higher redshift. In this work, we present the first direct measurement of the IGM temperature around a z=6 quasar by analysing the Doppler widths of Lya absorption lines in the proximity zone of SDSS J0818+1722. We use a high resolution (R= 40000) Keck/HIRES spectrum in combination with detailed numerical modelling to obtain the temperature at mean density, T_0=23600\pm^5000_6900K (\pm^9200_9300K) at 68 (95) per cent confidence assuming a prior probability 13500K<T_0<38500 K following HI and HeII reionisation. This enables us to place an upper limit on the redshift of HI reionisation, z_H, within 33 comoving Mpc of SDSS J0818+1722. If the quasar reionises the HeII in its vicinity, then in the limit of instantaneous reionisation we infer z_H<9.0 (11.0) at 68 (95) per cent confidence assuming photoheating is the dominant heat source and that HI reionisation is driven by ionising sources with soft spectra, typical of population II stars. If the HI and HeII in the IGM around SDSS J0818+1722 are instead reionised simultaneously by a population of massive metal-free stars, characterised by very hard ionising spectra, we obtain a tighter upper limit of z_H<8.4 (9.4). Initiating reionisation at higher redshifts produces temperatures which are too low with respect to our constraint unless the HI ionising sources or the quasar itself have spectra significantly harder than typically assumed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Detection of Extended He II Reionization in the Temperature Evolution of the Intergalactic Medium

We present new measurements of the temperature of the intergalactic medium (IGM) derived from the Lyman-alpha forest over 2.0 < z < 4.8. The small-scale structure in the forest of 61 high-resolution QSO spectra is quantified using a new statistic, the curvature, and the conversion to temperature calibrated using a suite of hydrodynamic simulations. At each redshift we focus on obtaining the temperature at an optimal overdensity probed by the Lyman-alpha forest, T(Delta), where the temperature is nearly a one-to-one function of the curvature regardless of the slope of the temperature-density relation. The median 2-sigma statistical uncertainty in these measurements is 8 per cent, though there may be comparable systematic errors due to the unknown amount of Jeans smoothing in the IGM. We use our T(Delta) results to infer the temperature at the mean density, T0. Even for a maximally steep temperature-density relation, T0 must increase from ~8000 K at z ~ 4.4 to >~12000 K at z ~ 2.8. This increase is not consistent with the monotonic decline in T0 expected in the absence of He II reionization. We therefore interpret the observed rise in temperature as evidence of He II reionization beginning at z >~ 4.4. The evolution of T0 is consistent with an end to He II reionization at z ~ 3, as suggested by opacity measurements of the He II Lyman-alpha forest, although the redshift at which T0 peaks will depend somewhat on the evolution of the temperature-density relation. These new temperature measurements suggest that the heat input due to the reionization of He II dominates the thermal balance of the IGM over an extended period with Delta_z >~ 1.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Measurements of the UV background at 4.6 < z < 6.4 using the quasar proximity effect

We present measurements of the ionising ultraviolet background (UVB) at z ~ 5-6 using the quasar proximity effect. The fifteen quasars in our sample cover the range 4.6 < z_q < 6.4, enabling the first proximity effect measurements of the UVB at z > 5. The metagalactic hydrogen ionisation rate, Gamma_bkg, was determined by modelling the combined ionisation field from the quasar and the UVB in the proximity zone on a pixel-by-pixel basis. The optical depths in the spectra were corrected for the expected effect of the quasar until the mean flux in the proximity region equalled that in the average Ly-alpha forest, and from this we make a measurement of Gamma_bkg. A number of systematic effects were tested using synthetic spectra. Noise in the flux was found to be the largest source of bias at z ~ 5, while uncertainties in the mean transmitted Ly-alpha flux are responsible for the largest bias at z ~ 6. The impacts of large-scale overdensities and Lyman limit systems on Gamma_bkg were also investigated, but found to be small at z > 5. We find a decline in Gamma_bkg with redshift, from log(Gamma_bkg) = -12.15 $\pm$ 0.16 at z ~ 5 to log(Gamma_bkg) = -12.84 $\pm$ 0.18 at z ~ 6 (1 sigma errors). Compared to UVB measurements at lower redshifts, our measurements suggest a drop of a factor of five in the HI photoionisation rate between z ~ 4 and z ~ 6. The decline of Gamma_bkg appears to be gradual, and we find no evidence for a sudden change in the UVB at any redshift that would indicate a rapid change in the attenuation length of ionising photons. Combined with recent measurements of the evolution of the mean free path of ionising photons, our results imply decline in the emissivity of ionising photons by roughly a factor of two from z ~ 5 to 6, albeit with significant uncertainty due to the measurement errors in both Gamma_bkg and the mean free path.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Near-zone sizes and the rest frame extreme ultra-violet spectral index of the highest redshift quasars

The discovery of quasars with redshifts higher than six has prompted a great deal of discussion in the literature regarding the role of quasars, both as sources of reionization, and as probes of the ionization state of the IGM. However the extreme ultra-violet (EUV) spectral index cannot be measured directly for high redshift quasars owing to absorption at frequencies above the Lyman limit, and as a result, studies of the impact of quasars on the intergalactic medium during reionization must assume a spectral energy distribution in the extreme ultra-violet based on observations at low redshift, z<1. In this paper we use regions of high Ly-alpha transmission (near-zones) around the highest redshift quasars to measure the quasar EUV spectral index at z~6. We jointly fit the available observations for variation of near-zone size with both redshift and luminosity, and propose that the observed relation provides evidence for an EUV spectral index that varies with absolute magnitude in the high redshift quasar sample, becoming softer at higher luminosity. Using a large suite of detailed numerical simulations we find that the typical value of spectral index for a luminous quasar at z~6 is constrained to be alpha=1.3+/-0.4 for a specific luminosity of the form L\proptoν^{-alpha}. We find the scatter in spectral index among individual quasars to be in the range ~0.75-1.25. These values are in agreement with direct observations at low redshift, and indicate that there has been no significant evolution in the EUV spectral index of quasars over 90% of cosmic time.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Resolving the high redshift Lyman-alpha forest in smoothed particle hydrodynamics simulations

We use a large set of cosmological smoothed particle hydrodynamics (SPH) simulations to examine the effect of mass resolution and box size on synthetic Lya forest spectra at 2 \leq z \leq 5. The mass resolution requirements for the convergence of the mean Lya flux and flux power spectrum at z=5 are significantly stricter than at lower redshift. This is because transmission in the high redshift Lya forest is primarily due to underdense regions in the intergalactic medium (IGM), and these are less well resolved compared to the moderately overdense regions which dominate the Lya forest opacity at z~2-3. We further find that the gas density distribution in our simulations differs significantly from previous results in the literature at large overdensities (Δ>10). We conclude that studies of the Lya forest at z=5 using SPH simulations require a gas particle mass of M_gas \leq 2x10^5 M_sol/h, which is >8 times the value required at z=2. A box size of at least 40 Mpc/h is preferable at all redshifts.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The evolution of the Lyman-alpha forest effective optical depth following HeII reionisation

Three independent observational studies have now detected a narrow (Δz ~ 0.5) dip centred at z=3.2 in the otherwise smooth redshift evolution of the Lya forest effective optical depth. This feature has previously been interpreted as an indirect signature of rapid photo-heating in the IGM during the epoch of HeII reionisation. We examine this interpretation using a semi-analytic model of inhomogeneous HeII reionisation and high resolution hydrodynamical simulations of the Lya forest. We instead find that a rapid (Δz ~ 0.2) boost to the IGM temperature (ΔT ~ 10^4 K) beginning at z=3.4 produces a well understood and generic evolution in the Lya effective optical depth, where a sudden reduction in the opacity is followed by a gradual, monotonic recovery driven largely by adiabatic cooling in the low density IGM. This behaviour is inconsistent with the narrow feature in the observational data. If photo-heating during HeII reionisation is instead extended over several redshift units, as recent theoretical studies suggest, then the Lya opacity will evolve smoothly with redshift. We conclude that the sharp dip observed in the Lya forest effective optical depth is instead most likely due to a narrow peak in the hydrogen photo-ionisation rate around z=3.2, and suggest that it may arise from the modulation of either reprocessed radiation during HeII reionisation, or the opacity of Lyman limit systems.

0 citations0 reviewsNo code linkedOpen access