We introduce meraxes, a new, purpose-built semi-analytic galaxy formation model designed for studying galaxy growth during reionization. meraxes is the first model of its type to include a temporally ...and spatially coupled treatment of reionization and is built upon a custom (100 Mpc)3
N-body simulation with high temporal and mass resolution, allowing us to resolve the galaxy and star formation physics relevant to early galaxy formation. Our fiducial model with supernova feedback reproduces the observed optical depth to electron scattering and evolution of the galaxy stellar mass function between z = 5 and 7, predicting that a broad range of halo masses contribute to reionization. Using a constant escape fraction and global recombination rate, our model is unable to simultaneously match the observed ionizing emissivity at z ≲ 6. However, the use of an evolving escape fraction of 0.05–0.1 at z ∼ 6, increasing towards higher redshift, is able to satisfy these three constraints. We also demonstrate that photoionization suppression of low-mass galaxy formation during reionization has only a small effect on the ionization history of the intergalactic medium. This lack of ‘self-regulation’ arises due to the already efficient quenching of star formation by supernova feedback. It is only in models with gas supply-limited star formation that reionization feedback is effective at regulating galaxy growth. We similarly find that reionization has only a small effect on the stellar mass function, with no observationally detectable imprint at M
* > 107.5 M⊙. However, patchy reionization has significant effects on individual galaxy masses, with variations of factors of 2–3 at z = 5 that correlate with environment.
Abstract
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 z ∼ 0.6–7, the black hole mass function at z ≲ 0.5, the global ionizing emissivity at z ∼ 2–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 ∼23.4 per cent due to obscuration by dust, the model is able to reproduce the observed quasar luminosity function at z ∼ 0.6–6. The stellar light from galaxies hosting faint active galactic nucleus (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 z ∼ 4 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 z ∼ 6. We conclude that quasars are unlikely to make a significant contribution to reionization.
We present the Dark-ages Reionization and Galaxy formation Observables from Numerical Simulations (DRAGONS) programme and Tiamat, the collisionless N-body simulation programme upon which DRAGONS is ...built. The primary trait distinguishing Tiamat from other large simulation programme is its density of outputs at high redshift (100 from z = 35 to z = 5; roughly one every 10 Myr) enabling the construction of very accurate merger trees at an epoch when galaxy formation is rapid and mergers extremely frequent. We find that the friends-of-friends halo mass function agrees well with the prediction of Watson et al. at high masses, but deviates at low masses, perhaps due to our use of a different halo finder or perhaps indicating a break from 'universal' behaviour. We then analyse the dynamical evolution of galaxies during the Epoch of Reionization finding that only a small fraction (~20 per cent) of galactic haloes are relaxed. We illustrate this using standard relaxation metrics to establish two dynamical recovery time-scales: (i) haloes need ~1.5 dynamical times following formation, and (ii) ~2 dynamical times following a major (3:1) or minor (10:1) merger to be relaxed. This is remarkably consistent across a wide mass range. Lastly, we use a phase-space halo finder to illustrate that major mergers drive long-lived massive phase-space structures which take many dynamical times to dissipate. This can yield significant differences in the inferred mass build-up of galactic haloes and we suggest that care must be taken to ensure a physically meaningful match between the galaxy formation physics of semi-analytic models and the halo finders supplying their input.
Low-frequency observatories are currently being constructed with the goal of detecting redshifted 21-cm emission from the epoch of reionization. These observatories will also be able to detect ...intensity fluctuations in the cumulative 21-cm emission after reionization, from hydrogen in unresolved damped Lyα absorbers (such as gas-rich galaxies) down to a redshift z∼ 3.5. The inferred power spectrum of 21-cm fluctuations at all redshifts will show acoustic oscillations, whose comoving scale can be used as a standard ruler to infer the evolution of the equation of state for the dark energy. We find that the first generation of low-frequency experiments (such as MWA or LOFAR) will be able to constrain the acoustic scale to within a few per cent in a redshift window just prior to the end of the reionization era, provided that foregrounds can be removed over frequency bandpasses of ≳8 MHz. This sensitivity to the acoustic scale is comparable to the best current measurements from galaxy redshift surveys, but at much higher redshifts. Future extensions of the first-generation experiments (involving an order of magnitude increase in the antennae number of the MWA) could reach sensitivities below 1 per cent in several redshift windows and could be used to study the dark energy in the unexplored redshift regime of 3.5 ≲z≲ 12. Moreover, new experiments with antennae designed to operate at higher frequencies would allow precision measurements (≲1 per cent) of the acoustic peak to be made at more moderate redshifts (1.5 ≲z≲ 3.5), where they would be competitive with ambitious spectroscopic galaxy surveys covering more than 1000 deg2. Together with other data sets, observations of 21-cm fluctuations will allow full coverage of the acoustic scale from the present time out to z∼ 12.
In this paper, we present calculations of the UV luminosity function (LF) from the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulations project, which combines N-body, ...semi-analytic and seminumerical modelling designed to study galaxy formation during the Epoch of Reionization. Using galaxy formation physics including supernova feedback, the model naturally reproduces the UV LFs for high-redshift star-forming galaxies from z ∼ 5 through to z ∼ 10. We investigate the luminosity–star formation rate (SFR) relation, finding that variable SFR histories of galaxies result in a scatter around the median relation of 0.1–0.3 dex depending on UV luminosity. We find close agreement between the model and observationally derived SFR functions. We use our calculated luminosities to investigate the LF below current detection limits, and the ionizing photon budget for reionization. We predict that the slope of the UV LF remains steep below current detection limits and becomes flat at M
UV ≳ −14. We find that 48 (17) per cent of the total UV flux at z ∼ 6 (10) has been detected above an observational limit of M
UV ∼ −17, and that galaxies fainter than M
UV ∼ −17 are the main source of ionizing photons for reionization. We investigate the luminosity–stellar mass relation, and find a correlation for galaxies with M
UV < −14 that has the form
$M_{\ast } \propto 10^{-0.47M_\mathrm{UV}}$
, in good agreement with observations, but which flattens for fainter galaxies. We determine the luminosity–halo mass relation to be
$M_\mathrm{vir} \propto 10^{-0.35M_\mathrm{UV}}$
, finding that galaxies with M
UV = −20 reside in host dark matter haloes of 1011.0±0.1 M⊙ at z ∼ 6, and that this mass decreases towards high redshift.
We investigate high-redshift galaxy sizes using a semi-analytic model constructed for the Dark-ages Reionization And Galaxy-formation Observables from Numerical Simulation project. Our fiducial ...model, including strong feedback from supernovae and photoionization background, accurately reproduces the evolution of the stellar mass function and ultraviolet (UV) luminosity function. Using this model, we study the size-luminosity relation of galaxies and find that the effective radius scales with UV luminosity as R sub( e) ... L super( 0.25) at z ~ 5-9. We show that recently discovered very luminous galaxies at z ... 7 and 11 lie on our predicted size-luminosity relations. We find that a significant fraction of galaxies at z > 8 will not be resolved by James Webb Space Telescope, but Giant Magellan Telescope will have the ability to resolve all galaxies in haloes above the atomic cooling limit. We show that our fiducial model successfully reproduces the redshift evolution of average galaxy sizes at z > 5. We also explore galaxy sizes in models without supernova feedback. The no-supernova feedback models produce galaxy sizes that are smaller than observations. We therefore confirm that supernova feedback plays an important role in determining the size-luminosity relation of galaxies and its redshift evolution during reionization. (ProQuest: ... denotes formulae/symbols omitted.)
We use the Dark-ages, Reionization And Galaxy formation Observables from Numerical Simulations (DRAGONS) framework to investigate the effect of galaxy formation physics on the morphology and ...statistics of ionized hydrogen (H ii) regions during the Epoch of Reioinization (EoR). DRAGONS self-consistently couples a semi-analytic galaxy formation model with the inhomogeneous ionizing UV background, and can therefore be used to study the dependence of morphology and statistics of reionization on feedback phenomena of the ionizing source galaxy population. Changes in galaxy formation physics modify the sizes of H ii regions and the amplitude and shape of 21-cm power spectra. Of the galaxy physics investigated, we find that supernova feedback plays the most important role in reionization, with H ii regions up to ≈20 per cent smaller and a fractional difference in the amplitude of power spectra of up to ≈17 per cent at fixed ionized fraction in the absence of this feedback. We compare our galaxy formation-based reionization models with past calculations that assume constant stellar-to-halo mass ratios and find that with the correct choice of minimum halo mass, such models can mimic the predicted reionization morphology. Reionization morphology at fixed neutral fraction is therefore not uniquely determined by the details of galaxy formation, but is sensitive to the mass of the haloes hosting the bulk of the ionizing sources. Simple EoR parametrizations are therefore accurate predictors of reionization statistics. However, a complete understanding of reionization using future 21-cm observations will require interpretation with realistic galaxy formation models, in combination with other observations.
Abstract
Massive quiescent galaxies (MQGs) are thought to have formed stars rapidly at early times followed by a long period of quiescence. The recent discovery of a MQG, ZF-COSMOS-20115 at z ∼ 4, ...only 1.5 Gyr after the big bang, places new constraints on galaxy growth and the role of feedback in early star formation. Spectroscopic follow-up confirmed ZF-COSMOS-20115 as a MQG at z = 3.717 with an estimated stellar mass of ∼1011 M⊙, showing no evidence of recent star formation. We use the Meraxes semi-analytic model to investigate how ZF-COSMOS-20115 analogues build stellar mass, and why they become quiescent. We identify three analogue galaxies with similar properties to ZF-COSMOS-20115. We find that ZF-COSMOS-20115 is likely hosted by a massive halo with virial mass of ∼1013 M⊙, having been through significant mergers at early times. These merger events drove intense growth of the nucleus, which later prevented cooling and quenched star formation. Therefore, ZF-COSMOS-20115 is unlikely to have experienced strong or extended star formation events at z < 3.7. We find that the analogues host the most massive black holes in our simulation and were luminous quasars at z ∼ 5, indicating that ZF-COSMOS-20115 and other MQGs may be the descendants of high-redshift quasars. In addition, the model suggests that ZF-COSMOS-20115 formed in a region of intergalactic medium that was reionized early.
We use high-resolution N-body simulations to study the concentration and spin parameters of dark matter haloes in the mass range 108 M⊙ h
−1 < M < 1011 M⊙ h
−1 and redshifts 5 < z < 10, corresponding ...to the haloes of galaxies thought to be responsible for reionization. We build a subsample of equilibrium haloes and contrast their properties to the full population that also includes unrelaxed systems. Concentrations are calculated by fitting both NFW and Einasto profiles to the spherically averaged density profiles of individual haloes. After removing haloes that are out of equilibrium, we find a z > 5 concentration–mass (c(M)) relation that is almost flat and well described by a simple power law for both NFW and Einasto fits. The intrinsic scatter around the mean relation is Δc
vir ∼ 1 (or 20 per cent) at z = 5. We also find that the analytic model proposed by Ludlow et al. reproduces the mass and redshift dependence of halo concentrations. Our best-fitting Einasto shape parameter, α, depends on peak height, ν, in a manner that is accurately described by α = 0.0070ν2 + 0.1839. The distribution of the spin parameter, λ, has a weak dependence on equilibrium state; λ peaks at roughly ∼0.033 for our relaxed sample, and at ∼0.04 for the full population. The spin–virial mass relation has a mild negative correlation at high redshift.
Abstract
The direct detection of regions of ionized hydrogen (H ii) has been suggested as a promising probe of cosmic reionization. Observing the redshifted 21-cm signal of hydrogen from the epoch of ...reionization (EoR) is a key scientific driver behind new-generation, low-frequency radio interferometers. We investigate the feasibility of combining low-frequency observations with the Square Kilometre Array and near infra-red survey data of the Wide-Field Infrared Survey Telescope to detect cosmic reionization by imaging H ii bubbles surrounding massive galaxies during the cosmic dawn. While individual bubbles will be too small to be detected, we find that by stacking redshifted 21-cm spectra centred on known galaxies, it will be possible to directly detect the EoR at z ∼ 9–12, and to place qualitative constraints on the evolution of the spin temperature of the intergalactic medium (IGM) at z ≳ 9. In particular, given a detection of ionized bubbles using this technique, it is possible to determine if the IGM surrounding them is typically in absorption or emission. Determining the globally averaged neutral fraction of the IGM using this method will prove more difficult due to degeneracy with the average size of H ii regions.