Cosmology

Abstract:

We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models' is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here.

Abstract:

To better understand the impact of supernova (SN) explosions on the evolution of galaxies, we perform a suite of high-resolution (12 pc), zoom-in cosmological simulations of a Milky Way-like galaxy at z = 3 with adaptive mesh refinement. We find that SN explosions can efficiently regulate star formation, leading to the stellar mass and metallicity consistent with the observed mass–metallicity relation and stellar mass–halo mass relation at z ~ 3. This is achieved by making three important changes to the classical feedback scheme: (i) the different phases of SN blast waves are modelled directly by injecting radial momentum expected at each stage, (ii) the realistic time delay of SNe is required to disperse very dense gas before a runaway collapse sets in, and (iii) a non-uniform density distribution of the interstellar medium (ISM) is taken into account below the computational grid scale for the cell in which an SN explodes. The simulated galaxy with the SN feedback model shows strong outflows, which carry approximately 10 times larger mass than star formation rate, as well as smoothly rising circular velocity. Although the metallicity of the outflow depends sensitively on the feedback model used, we find that the accretion rate and metallicity of the cold flow around the virial radius is impervious to SN feedback. Our results suggest that understanding the structure of the turbulent ISM may be crucial to assess the role of SN and other feedback processes in galaxy formation theory.

Abstract:

Galaxies and the dark matter haloes that host them are not spherically symmetric, yet spherical symmetry is a helpful simplifying approximation for idealized calculations and analysis of observational data. The assumption leads to an exact conservation of angular momentum for every particle, making the dynamics unrealistic. But how much does that inaccuracy matter in practice for analyses of stellar distribution functions, collisionless relaxation, or dark matter core-creation? We provide a general answer to this question for a wide class of aspherical systems; specifically, we consider distribution functions that are 'maximally stable', i.e. that do not evolve at first order when external potentials (which arise from baryons, large-scale tidal fields or infalling substructure) are applied. We show that a spherically symmetric analysis of such systems gives rise to the false conclusion that the density of particles in phase space is ergodic (a function of energy alone). Using this idea we are able to demonstrate that: (a) observational analyses that falsely assume spherical symmetry are made more accurate by imposing a strong prior preference for near-isotropic velocity dispersions in the centre of spheroids; (b) numerical simulations that use an idealized spherically symmetric setup can yield misleading results and should be avoided where possible; and (c) triaxial dark matter haloes (formed in collisionless cosmological simulations) nearly attain our maximally stable limit, but their evolution freezes out before reaching it.