This paper describes a new publicly available codebase for modeling galaxy formation in a cosmological context, the "Semi-Analytic Galaxy Evolution" model, or sage for short. super(5) sage is a ...significant update to the 2006 model of Croton et al. and has been rebuilt to be modular and customizable. The model will run on any N-body simulation whose trees are organized in a supported format and contain a minimum set of basic halo properties. In this work, we present the baryonic prescriptions implemented in sage to describe the formation and evolution of galaxies, and their calibration for three N-body simulations: Millennium, Bolshoi, and GiggleZ. Updated physics include the following: gas accretion, ejection due to feedback, and reincorporation via the galactic fountain; a new gas cooling-radio mode active galactic nucleus (AGN) heating cycle; AGN feedback in the quasar mode; a new treatment of gas in satellite galaxies; and galaxy mergers, disruption, and the build-up of intra-cluster stars. Throughout, we show the results of a common default parameterization on each simulation, with a focus on the local galaxy population.
The most luminous galaxies at high redshift are generally considered to be hosted in massive dark-matter halos of comparable number density, hence residing at the center of ...over-densities/protoclusters. We assess the validity of this assumption by investigating the clustering around the brightest galaxies populating the cosmic web at redshift z ∼ 8-9 through a combination of semi-analytic modeling and Monte Carlo simulations of mock Hubble Space Telescope WFC3 observations. The innovative aspect of our approach is the inclusion of a log-normal scatter parameter in the galaxy luminosity versus halo mass relation, extending the conditional luminosity function framework extensively used at low redshift to high z. Our analysis shows that the larger the value of , the less likely it is that the brightest source in a given volume is hosted in the most massive halo, and hence the weaker the overdensity of neighbors. We derive a minimum value of as a function of redshift by considering stochasticity in the halo assembly times, which affects galaxy ages and star formation rates in our modeling. We show that min(z) ∼ 0.15-0.3, with min increasing with redshift as a consequence of shorter halo assembly periods at higher redshifts. Current observations (mAB ∼ 27) of the environment of spectroscopically confirmed bright sources at z > 7.5 do not show strong evidence of clustering and are consistent with our modeling predictions for ≥ min. Deeper future observations reaching mAB ∼ 28.2-29 would have the opportunity to clearly quantify the clustering strength and hence to constrain , investigating the physical processes that drive star formation in the early universe.
ABSTRACT
We use the latest multiredshift (z = 6.5−8.7) upper limits on the 21-cm signal from the Murchison Widefield Array (MWA) to explore astrophysical models which are inconsistent with the data. ...We explore these limits in the context of reionization astrophysics by using 21CMMC to connect the disfavoured regions of parameter space to existing observational constraints on reionization such as high-z galaxy ultraviolet (UV) luminosity functions, the background UV photoionization rate, the intergalactic medium (IGM) neutral fraction, the electron scattering optical depth and the soft-band X-ray emissivity. We find the vast majority of disfavoured models to already be inconsistent with existing observational constraints. These can be broadly classified into two types of models: (i) ‘cold’ reionization and (ii) pure matter density fluctuations in a cold, neutral IGM (i.e. no reionization). Interestingly, a small subsample of models inconsistent with the MWA is consistent with the aforementioned constraints (excluding the X-ray emissivity). This implies that the current MWA limits are already providing unique information to disfavour models of reionization, albeit extremely weakly. We also provide the first limits on the soft-band X-ray emissivity from galaxies at high redshifts, finding 1σ lower limits of ϵX, 0.5−2 keV ≳ 1034.5 erg s−1 Mpc−3. Finally, we recover 95 per cent disfavoured limits on the IGM spin temperature of $\bar{T}_{\rm S}\lesssim$ 1.3, 1.4, 1.5, 1.8, 2.1, and 2.4 K at z = 6.5, 6.8, 7.1, 7.8, 8.2, and 8.7. With this, we infer the IGM must have undergone, at the very least, a small amount of X-ray heating. Note, the limits on ϵX, 0.5–2 keV and $\bar{T}_{\rm S}$ are conditional on the IGM neutral fraction.
We present the new semi-analytic model of galaxy evolution, Dark Sage, a heavily modified version of the publicly available sage code. The model is designed for detailed evolution of galactic discs. ...We evolve discs in a series of annuli with fixed specific angular momentum, which allows us to make predictions for the radial and angular-momentum structure of galaxies. Most physical processes, including all channels of star formation and associated feedback, are performed in these annuli. We present the surface density profiles of our model spiral galaxies, both as a function of radius and specific angular momentum, and find that the discs naturally build a pseudo-bulge-like component. Our main results are focused on predictions relating to the integrated mass-specific angular momentum relation of stellar discs. The model produces a distinct sequence between these properties in remarkable agreement with recent observational literature. We investigate the impact Toomre disc instabilities have on shaping this sequence and find they are crucial for regulating both the mass and spin of discs. Without instabilities, high-mass discs would be systematically deficient in specific angular momentum by a factor of ~2.5, with increased scatter. Instabilities also appear to drive the direction in which the mass-spin sequence of spiral galaxy discs evolves. With them, we find galaxies of fixed mass have higher specific angular momentum at later epochs.
We introduce a simple model to self-consistently connect the growth of galaxies to the formation history of their host dark matter haloes. Our model is defined by two simple functions: the 'baryonic ...growth function' which controls the rate at which new baryonic material is made available for star formation, and the 'physics function' which controls the efficiency with which this material is converted into stars. Using simple, phenomenologically motivated forms for both functions that depend only on a single halo property, we demonstrate the model's ability to reproduce the z = 0 red and blue stellar mass functions. Furthermore, by adding redshift as a second input variable to the physics function we show that the reproduction of the global stellar mass function out to z = 3 is improved. We conclude by discussing the general utility of our new model, highlighting its usefulness for creating mock galaxy samples which have a number of key advantages over those generated by other techniques.
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.
We introduce the Theoretical Astrophysical Observatory (TAO), an online virtual laboratory that houses mock observations of galaxy survey data. Such mocks have become an integral part of the modern ...analysis pipeline. However, building them requires expert knowledge of galaxy modeling and simulation techniques, significant investment in software development, and access to high performance computing. These requirements make it difficult for a small research team or individual to quickly build a mock catalog suited to their needs. To address this TAO offers access to multiple cosmological simulations and semi-analytic galaxy formation models from an intuitive and clean web interface. Results can be funnelled through science modules and sent to a dedicated supercomputer for further processing and manipulation. These modules include the ability to (1) construct custom observer light cones from the simulation data cubes; (2) generate the stellar emission from star formation histories, apply dust extinction, and compute absolute and/or apparent magnitudes; and (3) produce mock images of the sky. All of TAO's features can be accessed without any programming requirements. The modular nature of TAO opens it up for further expansion in the future.
Abstract
We introduce gbptrees: an algorithm for constructing merger trees from cosmological simulations, designed to identify and correct for pathological cases introduced by errors or ambiguities ...in the halo finding process. gbptrees is built upon a halo matching method utilizing pseudo-radial moments constructed from radially sorted particle ID lists (no other information is required) and a scheme for classifying merger tree pathologies from networks of matches made to-and-from haloes across snapshots ranging forward-and-backward in time. Focusing on subfind catalogues for this work, a sweep of parameters influencing our merger tree construction yields the optimal snapshot cadence and scanning range required for converged results. Pathologies proliferate when snapshots are spaced by ≲0.128 dynamical times; conveniently similar to that needed for convergence of semi-analytical modelling, as established by Benson et al. Total merger counts are converged at the level of ∼5 per cent for friends-of-friends (FoF) haloes of size n
p ≳ 75 across a factor of 512 in mass resolution, but substructure rates converge more slowly with mass resolution, reaching convergence of ∼10 per cent for n
p ≳ 100 and particle mass m
p ≲ 109 M⊙. We present analytic fits to FoF and substructure merger rates across nearly all observed galactic history (z ≤ 8.5). While we find good agreement with the results presented by Fakhouri et al. for FoF haloes, a slightly flatter dependence on merger ratio and increased major merger rates are found, reducing previously reported discrepancies with extended Press–Schechter estimates. When appropriately defined, substructure merger rates show a similar mass ratio dependence as FoF rates, but with stronger mass and redshift dependencies for their normalization.
ABSTRACT
The highly neutral intergalactic medium (IGM) during the epoch of reionization (EoR) is expected to suppress Ly α emission with damping-wing absorption, causing nearly no Ly α detection from ...star-forming galaxies at z ∼ 8. However, spectroscopic observations of the four brightest galaxies (H160 ∼ 25 mag) at these redshifts do reveal prominent Ly α line, suggesting locally ionized IGM. In this paper, we explore the Ly α 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 (H160 ≲ 25.5 mag) lie in the largest ionized regions in our simulation, and have low attenuation of their Ly α from the IGM (optical depth <1). Fainter galaxies (25.5 mag < H160 < 27.5 mag) have transmission that depends on UV luminosity, leading to a lower incidence of Ly α detection at fainter magnitudes. This luminosity-dependent attenuation explains why Ly α has only been observed in the brightest galaxies at z ∼ 8. Follow-up observations have revealed counterparts in the vicinity of these confirmed z ∼ 8 Ly α 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 α can be detected among fainter galaxies. At the current observational limit, galaxies with ≥2–5 neighbours within 2 arcmin × 2 arcmin are ∼2–3 times more likely to show Ly α emission. JWST will discover an order of magnitude more neighbours, revealing ≳50 galaxies in the largest ionizing bubbles and facilitating direct study of reionization morphology.
ABSTRACT
The fraction of ionizing photons that escape their host galaxies to ionize hydrogen in the intergalactic medium (IGM) is a critical parameter in analyses of the reionization era. In this ...paper, we use the meraxes semi-analytic galaxy formation model to infer the mean ionizing photon escape fraction and its dependence on galaxy properties through joint modelling of the observed high redshift galaxy population and existing constraints on the reionization history. Using a Bayesian framework, and under the assumption that escape fraction is primarily related to halo mass, we find that the joint constraints of the ultraviolet luminosity function, cosmic microwave background optical depth, and the Ly α forest require an escape fraction of $(18\pm 5)$ per cent for galaxies within haloes of M ≲ 109 M⊙ and $(5\pm 2)$ per cent for more massive haloes. In terms of galaxy properties, this transition in escape fraction occurs at stellar masses of M⋆ ∼ 107 M⊙, nearly independent of redshift. As a function of redshift, reionization is dominated by the smaller M⋆ ≲ 107 M⊙ galaxies with high escape fractions at z ≳ 6 and by the larger M⋆ ≳ 107 M⊙ galaxies with lower escape fractions at z ≲ 6. Galaxies with star formation rates of 10−2.5 M⊙yr−1 to 10−1.5 M⊙yr−1 provide the dominant source of ionizing photons throughout reionization. Our results are consistent with recent direct measurements of a $\sim 5~{{\ \rm per\ cent}}$ escape fraction from massive galaxies at the end of reionization and support the picture of low mass galaxies being the dominant sources of ionizing photons during reionization.