No Arabic abstract
The highly neutral inter-galactic medium (IGM) during the Epoch of Reionization (EoR) is expected to suppress Ly$alpha$ emission with damping-wing absorption, causing nearly no Ly$alpha$ detection from star-forming galaxies at $z{sim}8$. However, spectroscopic observations of the 4 brightest galaxies (${rm H}_{160}{sim}25$ mag) at these redshifts do reveal prominent Ly$alpha$ line, suggesting locally ionised IGM. In this paper, we explore the Ly$alpha$ IGM transmission and environment of bright galaxies during the EoR using the Meraxes semi-analytic model. We find brighter galaxies to be less affected by damping-wing absorption as they are effective at ionizing surrounding neutral hydrogen. Specifically, the brightest sources (${rm H}_{160}{lesssim}25.5$ mag) lie in the largest ionized regions in our simulation, and have low attenuation of their Ly$alpha$ from the IGM (optical depth ${<}1$). Fainter galaxies (25.5 mag${<}{rm H}_{160}{<}27.5$ mag) have transmission that depends on UV luminosity, leading to a lower incidence of Ly$alpha$ detection at fainter magnitudes. This luminosity-dependent attenuation explains why Ly$alpha$ has only been observed in the brightest galaxies at $z{sim}8$. Follow-up observations have revealed counterparts in the vicinity of these confirmed $z{sim}8$ Ly$alpha$ emitters. The environments of our modelled analogues agree with these observations in the number of nearby galaxies, which is a good indicator of whether Ly$alpha$ can be detected among fainter galaxies. At the current observational limit, galaxies with ${ge}2$--5 neighbours within $2{times}2$ are ${sim}2$--3 times more likely to show Ly$alpha$ emission. JWST will discover an order of magnitude more neighbours, revealing ${gtrsim}50$ galaxies in the largest ionizing bubbles and facilitating direct study of reionization morphology.
We investigate how the hydrostatic suppression of baryonic accretion affects the growth rate of dark matter halos during the Epoch of Reionization. By comparing halo properties in a simplistic hydrodynamic simulation in which gas only cools adiabatically, with its collisionless equivalent, we find that halo growth is slowed as hydrostatic forces prevent gas from collapsing. In our simulations, at the high redshifts relevant for reionization (between ${sim}6$ and ${sim}11$), halos that host dwarf galaxies ($lesssim 10^{9} mathrm{M_odot}$) can be reduced by up to a factor of 2 in mass due to the hydrostatic pressure of baryons. Consequently, the inclusion of baryonic effects reduces the amplitude of the low mass tail of the halo mass function by factors of 2 to 4. In addition, we find that the fraction of baryons in dark matter halos hosting dwarf galaxies at high redshift never exceeds ${sim}90%$ of the cosmic baryon fraction. When implementing baryonic processes, including cooling, star formation, supernova feedback and reionization, the suppression effects become more significant with further reductions of ${sim}30%$ to 60%. Although convergence tests suggest that the suppression may become weaker in higher resolution simulations, this suppressed growth will be important for semi-analytic models of galaxy formation, in which the halo mass inherited from an underlying N-body simulation directly determines galaxy properties. Based on the adiabatic simulation, we provide tables to account for these effects in N-body simulations, and present a modification of the halo mass function along with explanatory analytic calculations.
Motivated by recent measurements of the number density of faint AGN at high redshift, we investigate the contribution of quasars to reionization by tracking the growth of central supermassive black holes in an update of the Meraxes semi-analytic model. The model is calibrated against the observed stellar mass function at $zsim0.6-7$, the black hole mass function at $zlesssim0.5$, the global ionizing emissivity at $zsim2-5$ and the Thomson scattering optical depth. The model reproduces a Magorrian relation in agreement with observations at $z<0.5$ and predicts a decreasing black hole mass towards higher redshifts at fixed total stellar mass. With the implementation of an opening angle of 80 deg for quasar radiation, corresponding to an observable fraction of ${sim}23.4$ per cent due to obscuration by dust, the model is able to reproduce the observed quasar luminosity function at $zsim0.6-6$. The stellar light from galaxies hosting faint AGN contributes a significant or dominant fraction of the UV flux. At high redshift, the model is consistent with the bright end quasar luminosity function and suggests that the recent faint $zsim4$ AGN sample compiled by Giallongo et al. (2015) includes a significant fraction of stellar light. Direct application of this luminosity function to the calculation of AGN ionizing emissivity consequently overestimates the number of ionizing photons produced by quasars by a factor of 3 at $zsim6$. We conclude that quasars are unlikely to make a significant contribution to reionization.
We study dwarf galaxy formation at high redshift ($zge5$) using a suite of high- resolution, cosmological hydrodynamic simulations and a semi-analytic model (SAM). We focus on gas accretion, cooling and star formation in this work by isolating the relevant process from reionization and supernova feedback, which will be further discussed in a companion paper. We apply the SAM to halo merger trees constructed from a collisionless N-body simulation sharing identical initial conditions to the hydrodynamic suite, and calibrate the free parameters against the stellar mass function predicted by the hydrodynamic simulations at z = 5. By making comparisons of the star formation history and gas components calculated by the two modelling techniques, we find that semi-analytic prescriptions that are commonly adopted in the literature of low-redshift galaxy formation do not accurately represent dwarf galaxy properties in the hydrodynamic simulation at earlier times. We propose 3 modifications to SAMs that will provide more accurate high-redshift simulations. These include 1) the halo mass and baryon fraction which are overestimated by collisionless N-body simulations; 2) the star formation efficiency which follows a different cosmic evolutionary path from the hydrodynamic simulation; and 3) the cooling rate which is not well defined for dwarf galaxies at high redshift. Accurate semi-analytic modelling of dwarf galaxy formation informed by detailed hydrodynamical modelling will facilitate reliable semi-analytic predictions over the large volumes needed for the study of reionization.
Recent detections of Lyman alpha (Ly$alpha$) emission from $z>7.5$ galaxies were somewhat unexpected given a dearth of previous non-detections in this era when the intergalactic medium (IGM) is still highly neutral. But these detections were from UV bright galaxies, which preferentially live in overdensities which reionize early, and have significantly Doppler-shifted Ly$alpha$ line profiles emerging from their interstellar media (ISM), making them less affected by the global IGM state. Using a combination of reionization simulations and empirical ISM models we show, as a result of these two effects, UV bright galaxies in overdensities have $>2times$ higher transmission through the $zsim7$ IGM than typical field galaxies, and this boosted transmission is enhanced as the neutral fraction increases. The boosted transmission is not sufficient to explain the observed high Ly$alpha$ fraction of $M_mathrm{UV} lesssim -22$ galaxies (Stark et al. 2017), suggesting Ly$alpha$ emitted by these galaxies must be stronger than expected due to enhanced production and/or selection effects. Despite the bias of UV bright galaxies to reside in overdensities we show Ly$alpha$ observations of such galaxies can accurately measure the global neutral hydrogen fraction, particularly when Ly$alpha$ from UV faint galaxies is extinguished, making them ideal candidates for spectroscopic follow-up into the cosmic Dark Ages.
The epoch of reionization (6 < z < 10) marks the period in our universe when the first large galaxies grew to fruition, and began to affect the universe around them. Massive stars, and potentially accreting supermassive black holes, filled the universe with ionizing radiation, burning off the haze of neutral gas that had filled the intergalactic medium (IGM) since recombination (z~1000). The evolution of this process constrains key properties of these earliest luminous sources, thus observationally constraining reionization is a key science goal for the next decade. The measurement of Lyman-alpha emission from photometrically-identified galaxies is a highly constraining probe of reionization, as a neutral IGM will resonantly scatter these photons, reducing detectability. While significant work has been done with 8-10m telescopes, these observations require extremely large telescopes (ELTs); the flux limits available from todays 10m class telescopes are sufficient for only the brightest known galaxies (m < 26). Ultra-deep surveys with the Giant Magellan Telescope (GMT) and Thirty Meter Telescope (TMT) will be capable of detecting Lyman-alpha emission from galaxies 2-3 magnitudes fainter than todays deepest surveys. Wide-field fiber spectroscopy on the GMT combined with narrow-field AO-assisted slit spectroscopy on the TMT will be able to probe the expected size of ionized bubbles throughout the epoch of reionization, following up degree scale deep imaging surveys with the Wide Field Infrared Space Telescope. These data will provide the first resolved Lyman-alpha-based maps of the ionized intergalactic medium throughout the epoch of reionization, constraining models of both the temporal and spatial evolution of this phase change.