Search

We present a Bayesian inference on the neutral hydrogen fraction of the intergalactic medium (IGM), $\overline{x}_{\small HI}$, at z ∼ 6–8 using the properties of Lyman break galaxies (LBGs) during the epoch of reionization. We use large samples of LBG candidates at 5.5 ≤ z ≤ 8.2 with spectroscopy from Keck/DEIMOS and Keck/MOSFIRE. For each galaxy, we incorporate either the Lyman-α (Lyα) equivalent width (EW) for detections or the EW limit spectrum for non-detections to parametrize the EW distribution at various ultraviolet brightnesses for a given redshift. Using our reference sample of galaxy candidates from the ionized universe at z ∼ 6.0, we are able to infer $\overline{x}_{\small HI}$ at two redshifts: z ∼ 6.7 and z ∼ 7.6. This work includes intrinsically faint, gravitationally lensed galaxies at z ∼ 6.0 in order to constrain the intrinsic faint-end Lyα EW distribution and provide a comparable population of galaxies to counterparts in our sample that are at higher redshift. The inclusion of faint galaxy candidates, in addition to a more sophisticated modelling framework, allows us to better isolate effects of the interstellar medium and circumgalactic medium on the observed Lyα distribution from those of the IGM. We infer an upper limit of $\overline{x}_{\small HI}$ ≤ 0.25 (0.44) at z = 6.7 ± 0.2 and a neutral fraction of $\overline{x}_{\small HI}$ = $0.83^{+0.08}_{-0.11}$ (0.83$^{+0.11}_{-0.21}$) at z = 7.6 ± 0.6, both within 68 per cent (95 per cent) uncertainty, results that favour a moderately late and fairly rapid reionization.

We report the discovery of a double-peaked Lyman-$\alpha$ (Ly$\alpha$) emitter (LAE) at $z=3.2177\pm0.0001$ in VLT/MUSE data. The galaxy is strongly lensed by the galaxy cluster RXC~J0018.5+1626 recently observed in the RELICS survey, and the double-peaked Ly$\alpha$ emission is clearly detected in the two counter images in the MUSE field-of-view. We measure a relatively high Ly$\alpha$ rest-frame equivalent width (EW) of $\mathrm{EW}_{\mathrm{Ly}\alpha,0}=(63\pm2)\,\mathring{\mathrm{A}}$. Additional near-infrared (NIR) spectroscopy allows us to measure the H$\beta$, [OIII]$\lambda4959\,\mathring{\mathrm{A}}$ and [OIII]$\lambda5007\,\mathring{\mathrm{A}}$ emission lines, which show moderate rest-frame EWs of the order of a few $\sim10-100\,\mathring{\mathrm{A}}$, an [OIII]$\lambda5007\,\mathring{\mathrm{A}}$/H$\beta$ ratio of $4.8\pm0.7$, and a lower limit on the [OIII]/[OII] ratio of $>5.6$. The galaxy has very blue UV-continuum slopes of $\beta_{\mathrm{FUV}}=-2.23\pm0.06$ and $\beta_{\mathrm{NUV}}=-3.0\pm0.2$, and is magnified by factors $\mu\sim7-10$ in each of the two images, thus enabling a view into a low-mass ($M_{\star}\simeq10^{7.5}\,\mathrm{M}_{\odot}$) high-redshift galaxy analog. Notably, the blue peak of the Ly$\alpha$ profile is significantly stronger than the red peak, which suggests an inflow of matter and possibly very low HI column densities in its circumgalactic gas. Combined with the high lensing magnification and image multiplicity, these properties make this galaxy a prime candidate for follow-up observations to search for LyC emission and constrain the LyC photon escape fraction., Comment: Updated to match published version

Active galactic nuclei (AGN) forming in the early universe are thought to be the primary source of hard ionizing photons contributing to the reionization of intergalactic helium. However, the number density and spectral properties of high-redshift AGN remain largely unconstrained. In this work, we make use of physically-informed models calibrated with a wide variety of available observations to provide estimates for the role of AGN throughout the Epoch of Reionization. We present AGN luminosity functions in various bands between z = 2 to 7 predicted by the well-established Santa Cruz semi-analytic model, which includes modelling of black hole accretion and AGN feedback. We then combine the predicted AGN populations with a physical spectral model for self-consistent estimates of ionizing photon production rates, which depend on the mass and accretion rate of the accreting supermassive black hole. We then couple the predicted comoving ionizing emissivity with an analytic model to compute the subsequent reionization history of intergalactic helium and hydrogen. This work demonstrates the potential of coupling physically motivated analytic or semi-analytic techniques to capture multi-scale physical processes across a vast range of scales (here, from AGN accretion disks to cosmological scales). Our physical model predicts an intrinsic ionizing photon budget well above many of the estimates in the literature, meaning that helium reionization can comfortably be accomplished even with a relatively low escape fraction. We also make predictions for the AGN populations that are expected to be detected in future \emph{James Webb Space Telescope} surveys., Comment: 28 pages, 23 figures, accepted for publication in MNRAS

The transmission of Lyman-{\alpha} (Ly{\alpha}) in the spectra of distant quasars depends on the density, temperature, and ionization state of the intergalactic medium (IGM). Therefore, high-redshift (z > 5) Ly{\alpha} forests could be invaluable in studying the late stages of the epoch of reionization (EoR), as well as properties of the sources that drive it. Indeed, high-quality quasar spectra have now firmly established the existence of large-scale opacity fluctuations at z > 5, whose physical origins are still debated. Here we introduce a Bayesian framework capable of constraining the EoR and galaxy properties by forward-modelling the high-z Ly{\alpha} forest. Using priors from galaxy and CMB observations, we demonstrate that the final overlap stages of the EoR (when >95% of the volume was ionized) should occur at z < 5.6, in order to reproduce the large-scale opacity fluctuations seen in forest spectra. However, it is the combination of patchy reionization and the inhomogeneous UV background that produces the longest Gunn-Peterson troughs. Ly{\alpha} forest observations tighten existing constraints on the characteristic ionizing escape fraction of galaxies, with the combined observations suggesting f_{\rm esc} \approx 7^4_3%, and disfavoring a strong evolution with the galaxy's halo (or stellar) mass., Comment: 19 pages, 11 figures, submitted to MNRAS, comments are welcome

ALMA observations of $z>6$ galaxies have revealed abnormally high [OIII]$_{\rm 88\mu m}$/[CII]$_{\rm 158\mu m}$ ratios and [CII]$_{\rm 158\mu m}$ deficits compared to local galaxies. Numerous solutions have been proposed including differences in C and O abundance ratios, observational bias, and differences in ISM properties. In order to elucidate the underlying physics that drives this high-redshift phenomenon, we employ SPHINX$^{20}$, a state-of-the-art, cosmological radiation-hydrodynamics simulation, that resolves detailed ISM properties of thousands of galaxies in the epoch of reionization. We find that the observed $z>6$ [OIII]-SFR and [CII]-SFR relations can only be reproduced when the C/O abundance ratio is $\sim8\times$ lower than Solar and the total metal production is $\sim4\times$ higher than that of a Kroupa IMF. This implies that high-redshift galaxies are potentially primarily enriched by low-metallicity core-collapse supernovae with a more top-heavy IMF. As AGB stars and type-Ia supernova begin to contribute to the galaxy metallicity, both the [CII]-SFR and [CII] luminosity functions are predicted to converge to observed values at $z\sim4.5$. While we demonstrate that ionisation parameter, LyC escape fraction, ISM gas density, and CMB attenuation all drive galaxies towards higher [OIII]/[CII], observed values at $z>6$ can only be reproduced with substantially lower C/O abundances compared to Solar. The combination of [CII] and [OIII] can be used to predict the values of ionisation parameter, ISM gas density, and LyC escape fraction and we provide estimates of these quantities for nine observed $z>6$ galaxies. Finally, we demonstrate that [OI]$_{\rm 63\mu m}$ can be used as a replacement for [CII] in high-redshift galaxies where [CII] is unobserved and argue that more observation time should be used to target [OI] at $z>6$. (Abridged)

The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz., 15 pages, 15 figures, 1 table, Accepted to MNRAS

We compute the bispectra of the 21cm signal during the Epoch of Reionization for three different reionization scenarios that are based on a dark matter N-body simulation combined with a self-consistent, semi-numerical model of galaxy evolution and reionization. Our reionization scenarios differ in their trends of ionizing escape fractions ($f_\mathrm{esc}$) with the underlying galaxy properties and cover the physically plausible range, i.e. $f_\mathrm{esc}$ effectively decreasing, being constant, or increasing with halo mass. We find the 21cm bispectrum to be sensitive to the resulting ionization topologies that significantly differ in their size distribution of ionized and neutral regions throughout reionization. From squeezed to stretched triangles, the 21cm bispectra features a change of sign from negative to positive values, with ionized and neutral regions representing below-average and above-average concentrations contributing negatively and positively, respectively. The position of the change of sign provides a tracer of the size distribution of the ionized and neutral regions, and allows us to identify three major regimes that the 21cm bispectrum undergoes during reionization. In particular the regime during the early stages of reionization, where the 21cm bispectrum tracks the peak of the size distribution of the ionized regions, provides exciting prospects for pinning down reionization with the forthcoming Square Kilometre Array., Comment: 16 pages, 7 figures. Accepted for publication in MNRAS

A late end to reionisation at redshift $z\simeq 5.3$ is consistent with observed spatial variations in the Ly$\alpha$ forest transmission and the deficit of Ly$\alpha$ emitting galaxies around extended Ly$\alpha$ absorption troughs at $z=5.5$. In this model, large islands of neutral hydrogen should persist in the diffuse intergalactic medium (IGM) until $z\simeq 6$. We use a novel, hybrid approach that combines high resolution cosmological hydrodynamical simulations with radiative transfer to predict the incidence of strong 21 cm forest absorbers with optical depths $\tau_{21}>10^{-2}$ from the diffuse IGM in these late reionisation models. We include the effect of redshift space distortions on the simulated 21 cm forest spectra, and treat the highly uncertain heating of the pre-reionisation IGM by soft X-rays as a free parameter. For a model with only modest IGM pre-heating, such that average gas kinetic temperatures in the diffuse IGM remain below $T_{\rm K}\simeq 10^{2} \rm\, K$, we find that strong 21 cm forest absorption lines should persist until $z=6$. For a sample of $\sim 10$ sufficiently radio loud background sources, a null-detection of 21 cm forest absorbers at $z\simeq 6$ with SKA1-low or possibly LOFAR should provide an informative lower limit on the still largely unconstrained soft X-ray background at high redshift and the temperature of the pre-reionisation IGM.

We model the transmission of the Lyman-alpha line through the circum- and intergalactic media around dark matter haloes expected to host Lyman-alpha emitters (LAEs) at z > 5.7, using the high-dynamic-range Sherwood simulations. We find very different CGM environments around more massive haloes (~10^11 M_sun) compared to less massive haloes (~10^9 M_sun) at these redshifts, which can contribute to a different evolution of the Lyman-alpha transmission from LAEs within these haloes. Additionally we confirm that part of the differential evolution could result from bright LAEs being more likely to reside in larger ionized regions. We conclude that a combination of the CGM environment and the IGM ionization structure is likely to be responsible for the differential evolution of the bright and faint ends of the LAE luminosity function at z > 6. More generally, we confirm the suggestion that the self-shielded neutral gas in the outskirts of the host halo can strongly attenuate the Lyman-alpha emission from high redshift galaxies. We find that this has a stronger effect on the more massive haloes hosting brighter LAEs. The faint-end of the LAE luminosity function is thus a more reliable probe of the average ionization state of the IGM. Comparing our model for LAEs with a range of observational data we find that the favoured reionization histories are our previously advocated `Late' and `Very Late' reionization histories, in which reionization finishes rather rapidly at around z ~ 6., 26 pages, 16 figures, 2 appendices; accepted for publication by MNRAS

Supermassive primordial stars forming in atomically-cooled halos at z∼15−20 are currently thought to be the progenitors of the earliest quasars in the Universe. In this picture, the star evolves under accretion rates of 0.1−1 M⊙ yr−1 until the general relativistic instability triggers its collapse to a black hole at masses of ∼105 M⊙. However, the ability of the accretion flow to sustain such high rates depends crucially on the photospheric properties of the accreting star, because its ionising radiation could reduce or even halt accretion. Here we present new models of supermassive Population III protostars accreting at rates 0.001−10 M⊙ yr−1, computed with the GENEVA stellar evolution code including general relativistic corrections to the internal structure. We use the polytropic stability criterion to estimate the mass at which the collapse occurs, which has been shown to give a lower limit of the actual mass at collapse in recent hydrodynamic simulations. We find that at accretion rates higher than 0.001 M⊙ yr−1 the stars evolve as red, cool supergiants with surface temperatures below 104 K towards masses >105M⊙, and become blue and hot, with surface temperatures above 105 K,only for rates ≲0.001 M⊙ yr−1. Compared to previous studies, our results extend the range of masses and accretion rates at which the ionising feedback remains weak, reinforcing the case for direct collapse as the origin of the first quasars.

We introduce simulations aimed at assessing how well weak gravitational lensing of 21cm radiation from the Epoch of Reionization ($z \sim 8$) can be measured by an SKA-like radio telescope. A simulation pipeline has been implemented to study the performance of lensing reconstruction techniques. We show how well the lensing signal can be reconstructed using the three-dimensional quadratic lensing estimator in Fourier space assuming different survey strategies. The numerical code introduced in this work is capable of dealing with issues that can not be treated analytically such as the discreteness of visibility measurements and the inclusion of a realistic model for the antennae distribution. This paves the way for future numerical studies implementing more realistic reionization models, foreground subtraction schemes, and testing the performance of lensing estimators that take into account the non-Gaussian distribution of HI after reionization. If multiple frequency channels covering $z \sim 7-11.6$ are combined, Phase 1 of SKA-Low should be able to obtain good quality images of the lensing potential with a total resolution of $\sim 1.6$ arcmin. The SKA-Low Phase 2 should be capable of providing images with high-fidelity even using data from $z\sim 7.7 - 8.3$. We perform tests aimed at evaluating the numerical implementation of the mapping reconstruction. We also discuss the possibility of measuring an accurate lensing power spectrum. Combining data from $z \sim 7-11.6$ using the SKA2-Low telescope model, we find constraints comparable to sample variance in the range $L, 23 pages, 18 figures. Accepted for pubblication in MNRAS

We present new spectroscopy of the z = 3.62 gravitationally lensed quasar B1422+117 from the Gemini North GMOS integral field spectrograph. We observe significant differential magnifications between the broad emission lines and the continuum, as well as across the velocity structure of the Lyman-a line. We take advantage of this differential microlensing to algebraically decompose the quasar spectrum into the absorbed broad emission line and absorbed continuum components.We use the latter to derive the intrinsic Ly α forest absorption spectrum. The proximity effect is clearly detected, with a proximity zone edge of 8.6-17.3Mpc from the quasar, implying (perhaps intermittent) activity over at least 28Myr. The Ly α line profile exhibits a blue excess that is inconsistent with a symmetric fit to the unabsorbed red side. This has important implications for the use of this fitting technique in estimating the absorbed blue Ly α wings of Gunn-Peterson trough quasars.

The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz.

We present an analysis of the evolution of the Lyman-series forest into the epoch of reionization using cosmological radiative transfer simulations in a scenario where reionization ends late. We explore models with different midpoints of reionization and gas temperatures. We find that once the simulations have been calibrated to match the mean flux of the observed Lyman-$\alpha$ forest at $4 < z < 6$, they also naturally reproduce the distribution of effective optical depths of the Lyman-$\beta$ forest in this redshift range. We note that the tail of the largest optical depths that is most challenging to match corresponds to the long absorption trough of ULAS J0148+0600, which we have previously shown to be rare in our simulations. We consider the evolution of the Lyman-series forest out to higher redshifts, and show that future observations of the Lyman-$\beta$ forest at $z>6$ will discriminate between different reionization histories. The evolution of the Lyman-$\alpha$ and Lyman-$\gamma$ forests are less promising as a tool for pushing studies of reionization to higher redshifts due to the stronger saturation and foreground contamination, respectively., Comment: 10 pages, 9 figures, 1 table. Accepted to MNRAS

We apply Bayesian statistics to perform model selection on different reionisation scenarios via the Multinest algorithm. Initially, we recover the results shown by 21CMMC for the parameter estimation of 21cmFAST models. We proceed to test several toy models of the Epoch of Reionisation (EoR) defined in contrasting morphology and scale. We find that LOFAR observations are unlikely to allow model selection even with long integration times. HERA would require 61 dipoles to perform the same analysis in 1080 hours, and becomes comparable to the SKA with 217 dipoles. We find the SKA requires only 324 hours of observation to conclusively distinguish between our models. Once model selection is achievable, an analysis of observational priors is performed finding that neutral fraction checks at specific redshifts add little to no inference. We show the difficulties in model selection at the level of distinguishing fiducial parameters within a model or distinguishing galaxies with a constant versus power law mass-to-light ratio. Finally, we explore the use of the Savage-Dickey density ratio to show the redundancy of the parameter Rmfp within 21cmFAST., Comment: 20 pages, 16 figures - Accepted for publication by MNRAS

We use the observed unresolved cosmic X-ray background (CXRB) in the 0.5-2 keV band and existing upper limits on the 21-cm power spectrum to constrain the high-redshift population of X-ray sources, focusing on their effect on the thermal history of the Universe and the cosmic 21-cm signal. Because the properties of these sources are poorly constrained, we consider hot gas, X-ray binaries and mini-quasars (i.e., sources with soft or hard X-ray spectra) as possible candidates. We find that (1) the soft-band CXRB sets an upper limit on the X-ray efficiency of sources that existed before the end of reionization, which is one-to-two orders of magnitude higher than typically assumed efficiencies, (2) hard sources are more effective in generating the CXRB than the soft ones, (3) the commonly-assumed limit of saturated heating is not valid during the first half of reionization in the case of hard sources, with any allowed value of X-ray efficiency, (4) the maximal allowed X-ray efficiency sets a lower limit on the depth of the absorption trough in the global 21-cm signal and an upper limit on the height of the emission peak, while in the 21-cm power spectrum it sets a minimum amplitude and frequency for the high-redshift peaks, and (5) the existing upper limit on the 21-cm power spectrum sets a lower limit on the X-ray efficiency for each model. When combined with the 21-cm global signal, the CXRB will be useful for breaking degeneracies and helping constrain the nature of high-redshift heating sources., Comment: Submitted to MNRAS, 13 pages, 4 figures, 2 tables

One approach to extracting the global 21-cm signal from total-power measurements at low radio frequencies is to parametrize the different contributions to the data and then fit for these parameters. We examine parametrizations of the 21-cm signal itself, and propose one based on modelling the Ly α background, intergalactic medium temperature and hydrogen ionized fraction using tanh functions. This captures the shape of the signal from a physical modelling code better than an earlier parametrization based on interpolating between maxima and minima of the signal, and imposes a greater level of physical plausibility. This allows less biased constraints on the turning points of the signal, even though these are not explicitly fit for. Biases can also be alleviated by discarding information which is less robustly described by the parametrization, for example by ignoring detailed shape information coming from the covariances between turning points or from the high-frequency parts of the signal, or by marginalizing over the high-frequency parts of the signal by fitting a more complex foreground model. The fits are sufficiently accurate to be usable for experiments gathering 1000 h of data, though in this case it may be important to choose observing windows which do not include the brightest areas of the foregrounds. Our assumption of pointed, single-antenna observations and very broad-band fitting makes these results particularly applicable to experiments such as the Dark Ages Radio Explorer, which would study the global 21-cm signal from the clean environment of a low lunar orbit, taking data from the far side.

We present a search for metal absorption line systems at the highest redshifts to date using a deep (30h) VLT/X-Shooter spectrum of the z = 7.084 quasi-stellar object (QSO) ULAS J1120+0641. We detect seven intervening systems at z > 5.5, with the highest-redshift system being a C IV absorber at z = 6.51. We find tentative evidence that the mass density of C IV remains flat or declines with redshift at z < 6, while the number density of C II systems remains relatively flat over 5 < z < 7. These trends are broadly consistent with models of chemical enrichment by star formation-driven winds that include a softening of the ultraviolet background towards higher redshifts. We find a larger number of weak ( W_rest < 0.3A ) Mg II systems over 5.9 < z < 7.0 than predicted by a power-law fit to the number density of stronger systems. This is consistent with trends in the number density of weak Mg II systems at z = 2.5, and suggests that the mechanisms that create these absorbers are already in place at z = 7. Finally, we investigate the associated narrow Si IV, C IV, and N V absorbers located near the QSO redshift, and find that at least one component shows evidence of partial covering of the continuum source., 18 pages, 23 figures, accepted for publication in MNRAS

We have obtained spectra of 163 quasars at $z_\mathrm{em}>4.4$ with the Gemini Multi Object Spectrometers on the Gemini North and South telescopes, the largest publicly available sample of high-quality, low-resolution spectra at these redshifts. From this homogeneous data set, we generated stacked quasar spectra in three redshift intervals at $z\sim 5$. We have modelled the flux below the rest-frame Lyman limit ($\lambda_\mathrm{r}, Comment: 16 pages, 11 figures, plus 20 pages of online material; submitted to MNRAS

We investigate the origin of extragalactic continuum emission and its relation to the stellar population of a recently discovered peculiar z=3.344 Lyman alpha emitter. Based on an analysis of the broad-band colors and morphology we find further support for the idea that the underlying galaxy is being fed by a large-scale (L > 35 kpc) accretion stream. Archival HST images show small scale (~5 kpc) tentacular filaments converging near a hot-spot of star-formation, possibly fueled by gas falling in along the filaments. The spectral energy distribution of the tentacles is broadly compatible with either (1) non-ionizing rest-frame far-UV continuum emission from stars formed in an 60 million-year-old starburst; (2) nebular 2-photon-continuum radiation, arising from collisional excitation cooling, or (3) a recombination spectrum emitted by hydrogen fluorescing in response to ionizing radiation escaping from the galaxy. The latter possibility simultaneously accounts for the presence of asymmetric Lyman alpha emission from the large-scale gaseous filament and the nebular continuum in the smaller-scale tentacles as caused by the escape of ionizing radiation from the galaxy. Possible astrophysical explanations for the nature of the tentacles include: a galactic wind powered by the starburst; in-falling gas during cold accretion, or tails of interstellar medium dragged out of the galaxy by satellite halos that have plunged through the main halo. The possibility of detecting extragalactic 2-photon continuum emission in space-based, broad-band images suggests a tool for studying the gaseous environment of high redshift galaxies at much greater spatial detail than possible with Lyman alpha or other resonance line emission., 11 pages, 6 figures, accepted for publication by MNRAS

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., LCK acknowledges the support of an Isaac Newton Studentship, the Cambridge Trust and STFC. Support by the FP7 ERC Advanced Grant Emergence-320596 is gratefully acknowledged. EP acknowledges support from the Kavli Foundation. SB was supported by the National Aeronautics and Space Administration through Einstein Postdoctoral Fellowship Award Number PF5-160133 issued by the Chandra X-ray Observatory Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf of the National Aeronautics Space Administration under contract NAS8-03060. JSB acknowledges the support of a Royal Society University Research Fellowship. This work used the DIRAC Shared Memory Processing system at the University of Cambridge, operated by the COSMOS Project at the Department of Applied Mathematics and Theoretical Physics on behalf of the STFC DiRAC HPC Facility (www.dirac.ac.uk). This equipment was funded by BIS National E-infrastructure capital grant ST/J005673/1, STFC capital grant ST/H008586/1 and STFC DiRAC Operations grant ST/K00333X/1. This work also used the DiRAC Data Analytic system at the University of Cambridge, operated by the University of Cambridge High Performance Computing Service on behalf of the STFC DiRAC HPC Facility (www.dirac.ac.uk). This equipment was funded by BIS National E-infrastructure capital grant (ST/K001590/1), STFC capital grants ST/H008861/1 and ST/H00887X/1, and STFC DiRAC Operations grant ST/K00333X/1. DiRAC is part of the National E-Infrastructure. We acknowledge PRACE for awarding us access to the Curie supercomputer, based in France at the Tres Grand Centre de Calcul (TGCC), through the 8th regular call., This is the final version of the article. It first appeared from Oxford University Press via http://dx.doi.org/ 10.1093/mnras/stw1306

While the Lyman-$\alpha$ ($\mathrm{Ly}\alpha$) emission line serves as an important tool in the study of galaxies at $z\lesssim 6$, finding Ly$\alpha$ emitters (LAE) at significantly higher redshifts has been more challenging, probably because of the increasing neutrality of the intergalactic medium above $z\sim6$. Galaxies with extremely high rest-frame Ly$\alpha$ equivalent widths, EW(Ly$\alpha$) $\gtrsim 150$ \AA{}, at $z>6$ are good candidates for Ly$\alpha$ follow-up observations, and can stand out in multiband imaging surveys because of their unusual colors. We have conducted a photometric search for such objects in the Cluster Lensing And Supernova survey with Hubble (CLASH), and report here the identification of three likely gravitationally-lensed images of a single LAE candidate at $z\sim6.3$, behind the galaxy cluster Abell 2261($z = 0.225$). In the process, we also measured with Keck/MOSFIRE the first spectroscopic redshift of a multiply-imaged galaxy behind Abell 2261, at $z = 3.337$. This allows us to calibrate the lensing model, which in turn is used to study the properties of the candidate LAE. Population III galaxy spectral energy distribution (SED) model fits to the CLASH broadband photometry of the possible LAE provide a slightly better fit than Population I/II models. The best fitted model suggests intrinsic EW(Ly$\alpha$) $\approx 160$ \AA{} after absorption in the interstellar and intergalactic medium. Future spectroscopic observations will examine this prediction as well as shed more light on the morphology of this object, which indicates it may be a merger of two smaller galaxies., Comment: 13 pages, 6 figures

The quasar ULAS J1120+0641 at redshift z = 7.085 has a highly ionized 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 Lyα redshift. We use radiative transfer simulations in a cosmological context to investigate the implications for the ionization state of the inhomogeneous intergalactic medium (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 H I fraction of ? 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, ?, if a damped Lyα system in an otherwise highly ionized IGM lies within 5 proper Mpc of the quasar. This is, however, predicted to occur in only ˜5 per cent of our simulated sightlines 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 2 × 10^9 Msun black hole at this high redshift.