We present a detailed comparison of the substructure properties of a single Milky Way sized dark matter halo from the Aquarius suite at five different resolutions, as identified by a variety of different (sub)halo finders for simulations of cosmic structure formation. These finders span a wide range of techniques and methodologies to extract and quantify substructures within a larger non‐homogeneous background density (e.g. a host halo). This includes real‐space‐, phase‐space‐, velocity‐space‐ and time‐space‐based finders, as well as finders employing a Voronoi tessellation, Friends‐of‐Friends techniques or refined meshes as the starting point for locating substructure. A common post‐processing pipeline was used to uniformly analyse the particle lists provided by each finder. We extract quantitative and comparable measures for the subhaloes, primarily focusing on mass and the peak of the rotation curve for this particular study. We find that all of the finders agree extremely well in the presence and location of substructure and even for properties relating to the inner part of the subhalo (e.g. the maximum value of the rotation curve). For properties that rely on particles near the outer edge of the subhalo the agreement is at around the 20 per cent level. We find that the basic properties (mass and maximum circular velocity) of a subhalo can be reliably recovered if the subhalo contains more than 100 particles although its presence can be reliably inferred for a lower particle number limit of 20. We finally note that the logarithmic slope of the subhalo cumulative number count is remarkably consistent and <1 for all the finders that reached high resolution. If correct, this would indicate that the larger and more massive, respectively, substructures are the most dynamically interesting and that higher levels of the (sub)subhalo hierarchy become progressively less important.

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and V max function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -to a minor extent -the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

We present results from the semi-analytic model of galaxy formation sag applied on the MultiDark simulation MDPL2. sag features an updated supernova (SN) feedback scheme and a robust modelling of the environmental effects on satellite galaxies. This incorporates a gradual starvation of the hot gas halo driven by the action of ram pressure stripping (RPS), that can affect the cold gas disc, and tidal stripping (TS), which can act on all baryonic components. Galaxy orbits of orphan satellites are integrated providing adequate positions and velocities for the estimation of RPS and TS. The star formation history and stellar mass assembly of galaxies are sensitive to the redshift dependence implemented in the SN feedback model. We discuss a variant of our model that allows to reconcile the predicted star formation rate density at z 3 with the observed one, at the expense of an excess in the faint end of the stellar mass function at z = 2. The fractions of passive galaxies as a function of stellar mass, halo mass and the halo-centric distances are consistent with observational measurements. The model also reproduces the evolution of the main sequence of star forming central and satellite galaxies. The similarity between them is a result of the gradual starvation of the hot gas halo suffered by satellites, in which RPS plays a dominant role. RPS of the cold gas does not affect the fraction of quenched satellites but it contributes to reach the right atomic hydrogen gas content for more massive satellites (M 10 10 M ).

We present the public release of the MULTIDARK-GALAXIES: three distinct galaxy catalogues derived from one of the Planck cosmology MULTIDARK simulations (i.e. MDPL2, with a volume of (1 h −1 Gpc) 3 and mass resolution of 1.5 × 10 9 h −1 M ) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data for a selection of fundamental properties of galaxies like stellar mass function, star formation rate, cold gas fractions, and metallicities -noting that they sometimes perform differently reflecting model designs and calibrations. We have further selected galaxy subsamples of the catalogues by number densities in stellar mass, cold gas mass, and star formation rate in order to study the clustering statistics of galaxies. We show that despite different treatment of orphan galaxies, i.e. galaxies that lost their dark-matter host halo due to the finite mass resolution of the N -body simulation or tidal stripping, the clustering signal is comparable, and reproduces the observations in all three models -in particular when selecting samples based upon stellar mass. Our catalogues provide a powerful tool to study galaxy formation within a volume comparable to those probed by on-going and future photometric and redshift surveys. All model data consisting of a range of galaxy propertiesincluding broad-band SDSS magnitudes -are publicly available.

Galaxy surveys targeting emission lines are characterising the evolution of star-forming galaxies, but there is still little theoretical progress in modelling their physical properties. We predict nebular emission from star-forming galaxies within a cosmological galaxy formation model. Emission lines are computed by combining the semi-analytical model SAG with the photoionisation code MAPPINGS-III. We characterise the interstellar medium (ISM) of galaxies by relating the ionisation parameter of gas in galaxies to their cold gas metallicity, obtaining a reasonable agreement with the observed Hα, [OII]λ 3727, [OIII]λ 5007 luminosity functions, and the the BPT diagram for local star-forming galaxies. The average ionisation parameter is found to increase towards low star-formation rates and high redshifts, consistent with recent observational results. The predicted link between different emission lines and their associated star-formation rates is studied by presenting scaling relations to relate them. Our model predicts that emission line galaxies have modest clustering bias, and thus reside in dark matter haloes of masses below M halo 10 12 [h −1 M ]. Finally, we exploit our modelling technique to predict galaxy number counts up to z ∼ 10 by targeting far-infrared (FIR) emission lines detectable with submillimetre facilities.

We make use of a semi-analytical model of galaxy formation to investigate the origin of the observed correlation between [α/Fe] abundance ratios and stellar mass in elliptical galaxies. We implement a new galaxy-wide stellar initial mass function (Top Heavy Integrated Galaxy Initial Mass Function, TH-IGIMF) in the semi-analytic model SAG and evaluate its impact on the chemical evolution of galaxies. The SFR-dependence of the slope of the TH-IGIMF is found to be key to reproducing the correct [α/Fe]-stellar mass relation. Massive galaxies reach higher [α/Fe] abundance ratios because they are characterized by more top-heavy IMFs as a result of their higher SFR. As a consequence of our analysis, the value of the minimum embedded star cluster mass and of the slope of the embedded cluster mass function, which are free parameters involved in the TH-IGIMF theory, are found to be as low as 5 M ⊙ and 2, respectively. A mild downsizing trend is present for galaxies generated assuming either a universal IMF or a variable TH-IGIMF. We find that, regardless of galaxy mass, older galaxies (with formation redshifts 2) are formed in shorter time-scales ( 2 Gyr), thus achieving larger [α/Fe] values. Hence, the time-scale of galaxy formation alone cannot explain the slope of the [α/Fe]-galaxy mass relation, but is responsible for the big dispersion of [α/Fe] abundance ratios at fixed stellar mass. We further test the hyphothesis of a TH-IGIMF in elliptical galaxies by looking into mass-to-light ratios, and luminosity functions. Models with a TH-IGIMF are also favoured by these constraints. In particular, mass-to-light ratios agree with observed values for massive galaxies while being overpredicted for less massive ones; this overprediction is present regardless of the IMF considered.

We present a study of a comparison of spin distributions of subhaloes found associated with a host halo. The subhaloes are found within two cosmological simulation families of Milky Way-like galaxies, namely the Aquarius and GHALO simulations. These two simulations use different gravity codes and cosmologies. We employ ten different substructure finders, which span a wide range of methodologies from simple overdensity in configuration space to full 6-d phase space analysis of particles. We subject the results to a common post-processing pipeline to analyse the results in a consistent manner, recovering the dimensionless spin parameter. We find that spin distribution is an excellent indicator of how well the removal of background particles (unbinding) has been carried out. We also find that the spin distribution decreases for substructure the nearer they are to the host halo's, and that the value of the spin parameter rises with enclosed mass towards the edge of the substructure. Finally subhaloes are less rotationally supported than field haloes, with the peak of the spin distribution having a lower spin parameter.

We present a new technique for modelling submillimetre galaxies (SMGs): the 'Count Matching' approach. Using light cones drawn from a semi-analytic model of galaxy formation, we choose physical galaxy properties given by the model as proxies for their submillimetre luminosities, assuming a monotonic relationship. As recent interferometric observations of the Extended Chandra Deep Field-South show that the brightest sources detected by single-dish telescopes are comprised by emission from multiple fainter sources, we assign the submillimetre fluxes so that the combined Large APEX BOlometer CAmera (LABOCA) plus bright-end Atacama Large Millimeter/submillimeter Array observed number counts for this field are reproduced. After turning the model catalogues given by the proxies into submillimetre maps, we perform a source extraction to include the effects of the observational process on the recovered counts and galaxy properties. We find that for all proxies, there are lines of sight giving counts consistent with those derived from LABOCA observations, even for input sources with randomized positions in the simulated map. Comparing the recovered redshift, stellar mass and host halo mass distributions for model SMGs with observational data, we find that the best among the proposed proxies is that in which the submillimetre luminosity increases monotonically with the product between dust mass and star formation rate (SFR). This proxy naturally reproduces a positive trend between SFR and bolometric IR luminosity. The majority of components of blended sources are spatially unassociated.