We present an analysis of the pairwise velocity statistics from a suite of cosmological N-body simulations describing the "Running Friedmann-Lema\^itre-Robertson-Walker" (R-FLRW) cosmological model. This model is based on quantum field theory in a curved space-time and extends {\Lambda}CDM with a time-evolving vacuum energy density. To enforce local conservation of matter a time-evolving gravitational coupling is also included. Our results constitute the first study of velocities in the R-FLRW cosmology, and we also compare with other dark energy simulations suites, repeating the same analysis. We find a strong degeneracy between the pairwise velocity and {\sigma}_8 at z=0 for almost all scenarios considered, which remains even when we look back to epochs as early as z=2. We also investigate various Coupled Dark Energy models, some of which show minimal degeneracy, and reveal interesting deviations from {\Lambda}CDM which could be readily exploited by future cosmological observations to test and further constrain our understanding of dark energy.

We present a suite of cosmological N-body simulations describing the "Running Friedmann-Lema{\"i}tre-Robertson-Walker" (R-FLRW) cosmological model. This model is based on quantum field theory in a curved space-time and extends {\Lambda}CDM with a time-evolving vacuum density, {\Lambda}(z), and time-evolving gravitational Newton's coupling, G(z). In this paper we review the model and introduce the necessary analytical treatment needed to adapt a reference N-body code. Our resulting simulations represent the first realisation of the full growth history of structure in the R-FLRW cosmology into the non-linear regime, and our normalisation choice makes them fully consistent with the latest cosmic microwave background data. The post-processing data products also allow, for the first time, an analysis of the properties of the halo and sub-halo populations. We explore the degeneracies of many statistical observables and discuss the steps needed to break them. Furthermore, we provide a quantitative description of the deviations of R-FLRW from {\Lambda}CDM, which could be readily exploited by future cosmological observations to test and further constrain the model.

This paper describes a new publicly available codebase for modelling galaxy formation in a cosmological context, the "Semi-Analytic Galaxy Evolution" model, or SAGE for short. SAGE is a significant update to that used in Croton et al. (2006) and has been rebuilt to be modular and customisable. The model will run on any N-body simulation whose trees are organised 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: 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.

With the installation of the Cosmic Origins Spectrograph on the Hubble Space Telescope, measurements of the metal content of the low redshift intergalactic medium (IGM) are now available. Using a new grid-based model for diffuse gas coupled to the SAGE semi-analytic model of galaxy formation, we examine the impact of supernova feedback on the pollution of the IGM. We consider different assumptions for the reheating and ejection of gas by supernovae and their dependence on galaxy circular velocity and gas surface density. Where metals are present, we find the most likely metallicity to be $-1.5 < $log$_{10}$(Z/Z$_{\odot}$)$< -1.0$ at $z = 0$, consistent with both observations and more sophisticated hydrodynamic simulations. Our model predicts that the regions of the IGM with the highest metallicities will be near galaxies with M$_{\star} \sim 10^{10.5}h^{-1}$M$_{\odot}$ and in environments of densities $\sim 10 \times$ the mean. We also find that 90% of IGM metals at $z = 0$ are ejected by galaxies with stellar masses less than $10^{10.33}h^{-1}$M$_{\odot}$.