Search

We present a clustering comparison of 12 galaxy formation models [including semi-analytic models (SAMs) and halo occupation distribution (HOD) models] all run on halo catalogues and merger trees extracted from a single _ cold dark matter N-body simulation. We compare the results of the measurements of the mean halo occupation numbers, the radial distribution of galaxies in haloes and the two-point correlation functions (2PCF). We also study the implications of the different treatments of orphan (galaxies not assigned to any dark matter subhalo) and non-orphan galaxies in these measurements. Our main result is that the galaxy formation models generally agree in their clustering predictions but they disagree significantly between HOD and SAMs for the orphan satellites. Although there is a very good agreement between the models on the 2PCF of central galaxies, the scatter between the models when orphan satellites are included can be larger than a factor of 2 for scales smaller than 1 h?1 Mpc. We also show that galaxy formation models that do not include orphan satellite galaxies have a significantly lower 2PCF on small scales, consistent with previous studies. Finally, we show that the 2PCF of orphan satellites is remarkably different between SAMs and HOD models. Orphan satellites in SAMs present a higher clustering than in HOD models because they tend to occupy more massive haloes. We conclude that orphan satellites have an important role on galaxy clustering and they are the main cause of the differences in the clustering between HOD models and SAMs.

We present a method that extends the capabilities of the PINpointing Orbit-Crossing Collapsed HIerarchical Objects (PINOCCHIO) code, allowing it to generate accurate dark matter halo mock catalogues in cosmological models where the linear growth factor and the growth rate depend on scale. Such cosmologies comprise, among others, models with massive neutrinos and some classes of modified gravity theories. We validate the code by comparing the halo properties from PINOCCHIO against N-body simulations, focusing on cosmologies with massive neutrinos: $\nu\Lambda$CDM. We analyse the halo mass function, halo two-point correlation function, halo power spectrum and the moments of the halo density field, showing that PINOCCHIO reproduces the results from simulations with the same level of precision as the original code ($\sim5-10\%$). We demonstrate that the abundance of halos in cosmologies with massless and massive neutrinos from PINOCCHIO matches very well the outcome of simulations, and point out that PINOCCHIO can reproduce the $\Omega_\nu-\sigma_8$ degeneracy that affects the halo mass function. We show that the clustering properties of the halos from PINOCCHIO matches accurately those from simulations both in real and redshift-space, in the latter case up to $k=0.3~h~{\rm Mpc}^{-1}$. We finally point out that the first moments of the halo density field from simulations are precisely reproduced by PINOCCHIO. We emphasize that the computational time required by PINOCCHIO to generate mock halo catalogues is orders of magnitude lower than the one needed for N-body simulations. This makes this tool ideal for applications like covariance matrix studies within the standard $\Lambda$CDM model but also in cosmologies with massive neutrinos or some modified gravity theories., Comment: 25 pages, 7 figures

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.

We present results of cosmological simulations of disk galaxies carried out with the GADGET-3 TreePM+SPH code, where star formation and stellar feedback are described using our MUlti Phase Particle Integrator (MUPPI) model. This description is based on simple multi-phase model of the interstellar medium at unresolved scales, where mass and energy flows among the components are explicitly followed by solving a system of ordinary differential equations. Thermal energy from SNe is injected into the local hot phase, so as to avoid that it is promptly radiated away. A kinetic feedback prescription generates the massive outflows needed to avoid the over-production of stars. We use two sets of zoomed-in initial conditions of isolated cosmological halos with masses (2-3) * 10^{12} Msun, both available at several resolution levels. In all cases we obtain spiral galaxies with small bulge-over-total stellar mass ratios (B/T \approx 0.2), extended stellar and gas disks, flat rotation curves and realistic values of stellar masses. Gas profiles are relatively flat, molecular gas is found to dominate at the centre of galaxies, with star formation rates following the observed Schmidt-Kennicutt relation. Stars kinematically belonging to the bulge form early, while disk stars show a clear inside-out formation pattern and mostly form after redshift z=2. However, the baryon conversion efficiencies in our simulations differ from the relation given by Moster et al. (2010) at a 3 sigma level, thus indicating that our stellar disks are still too massive for the Dark Matter halo in which they reside. Results are found to be remarkably stable against resolution. This further demonstrates the feasibility of carrying out simulations producing a realistic population of galaxies within representative cosmological volumes, at a relatively modest resolution., 19 pages, 21 figures, MNRAS accepted

In this PhD thesis I present the implementation of a method to study the evolution of the large scale structure in a weak-field approximation of General Relativity. Our starting point is the assumption that we can reproduce the classical evolution of the universe with an N-body code in which small scales---below 1 Mpc---interactions are particle-to-particle in the Newtonian limit, while in the large scales these interaction are mediated through the relativistic gravitational field in the weak-field approximation. We have implemented this idea by first redesigning Gevolution code from a monolithic application into a relativistic Particle-Mesh (PM) library libgevolution. Then we have modified Gadget4, a Newtonian TreePM code, adding libgevolution as a plug-in replacement to its original Newtonian PM. This project goes by the name of GrGadget. The advantage of a combined Tree+PM approach with respect to a pure PM method lies in the fact that we can simulate huge cubic boxes, with a length of the order of the 2 Gpc/h, representing a portion of the visible universe and still be able to resolve structures where dark matter halos hosting galaxies are formed on scales below 100 kpc/h using reasonable computational resources. Furthermore one can scale the size of the box without necessarily needing to increase the memory requirements or sacrifying the small scale resolution. Running GrGadget we have also realized that failing to resolve the non-linear structures, for example in the case of a pure PM code, would underestimate the amplitude of the highest Fourier modes of the relativistic fields by more that a 30% factor. In this PhD thesis I present the implementation of a method to study the evolution of the large scale structure in a weak-field approximation of General Relativity. Our starting point is the assumption that we can reproduce the classical evolution of the universe with an N-body code in which small scales---below 1 Mpc---interactions are particle-to-particle in the Newtonian limit, while in the large scales these interaction are mediated through the relativistic gravitational field in the weak-field approximation. We have implemented this idea by first redesigning Gevolution code from a monolithic application into a relativistic Particle-Mesh (PM) library libgevolution. Then we have modified Gadget4, a Newtonian TreePM code, adding libgevolution as a plug-in replacement to its original Newtonian PM. This project goes by the name of GrGadget. The advantage of a combined Tree+PM approach with respect to a pure PM method lies in the fact that we can simulate huge cubic boxes, with a length of the order of the 2 Gpc/h, representing a portion of the visible universe and still be able to resolve structures where dark matter halos hosting galaxies are formed on scales below 100 kpc/h using reasonable computational resources. Furthermore one can scale the size of the box without necessarily needing to increase the memory requirements or sacrifying the small scale resolution. Running GrGadget we have also realized that failing to resolve the non-linear structures, for example in the case of a pure PM code, would underestimate the amplitude of the highest Fourier modes of the relativistic fields by more that a 30% factor.

Data analysis in cosmology requires reliable covariance matrices. Covariance matrices derived from numerical simulations often require a very large number of realizations to be accurate. When a theoretical model for the covariance matrix exists, the parameters of the model can often be fit with many fewer simulations. We write a likelihood-based method for performing such a fit. We demonstrate how a model covariance matrix can be tested by examining the appropriate $\chi^2$ distributions from simulations. We show that if model covariance has amplitude freedom, the expectation value of second moment of $\chi^2$ distribution with a wrong covariance matrix will always be larger than one using the true covariance matrix. By combining these steps together, we provide a way of producing reliable covariances without ever requiring running a large number of simulations. We demonstrate our method on two examples. First, we measure the two-point correlation function of halos from a large set of $10000$ mock halo catalogs. We build a model covariance with $2$ free parameters, which we fit using our procedure. The resulting best-fit model covariance obtained from just $100$ simulation realizations proves to be as reliable as the numerical covariance matrix built from the full $10000$ set. We also test our method on a setup where the covariance matrix is large by measuring the halo bispectrum for thousands of triangles for the same set of mocks. We build a block diagonal model covariance with $2$ free parameters as an improvement over the diagonal Gaussian covariance. Our model covariance passes the $\chi^2$ test only partially in this case, signaling that the model is insufficient even using free parameters, but significantly improves over the Gaussian one., Comment: Accepted for publication in JCAP. 24 pages, 8 figures

We present a new implementation of the GAlaxy Evolution and Assembly (GAEA) semi-analytic model, that features an improved modelling of the process of cold gas accretion onto supermassive black hole (SMBHs), derived from both analytic arguments and high-resolution simulations. We consider different scenarios for the loss of angular momentum required for the available cold gas to be accreted onto the central SMBHs, and we compare different combinations of triggering mechanisms, including galaxy mergers and disc instabilities in star forming discs. We compare our predictions with the luminosity function (LF) observed for Active Galactic Nuclei (AGN) and we confirm that a non-instantaneous accretion timescale (either in the form of a low-angular momentum reservoir or as an assumed light curve evolution) is needed in order to reproduce the measured evolution of the AGN-LF and the so-called AGN-downsizing trend. Moreover, we also study the impact of AGN feedback, in the form of AGN-driven outflows, on the SF properties of model galaxies, using prescriptions derived both from empirical studies or from numerical experiments. We show that AGN-driven outflows are effective in suppressing the residual star formation rate in massive galaxies ($> 10^{11} M_\odot$) without changing their overall assembly history. These winds also affect the SFR of lower mass galaxies, resulting in a too large fraction of passive galaxies at $< 10^{10} M_\odot$. Finally, we study the Eddington ratio distribution as a function of SMBH mass, showing that only objects more massive than $10^8 M_\odot$ are already in a self-regulated state as inferred from observations., 19 pages, 12 figures, 1 table, replaced with MNRAS accepted version

The generation of simulated convergence maps is of key importance in fully exploiting weak lensing by Large Scale Structure (LSS) from which cosmological parameters can be derived. In this paper we present an extension of the PINOCCHIO code which produces catalogues of dark matter haloes so that it is capable of simulating weak lensing by LSS. Like WL-MOKA, the method starts with a random realisation of cosmological initial conditions, creates a halo catalogue and projects it onto the past-light-cone, and paints in haloes assuming parametric models for the mass density distribution within them. Large scale modes that are not accounted for by the haloes are constructed using linear theory. We discuss the systematic errors affecting the convergence power spectra when Lagrangian Perturbation Theory at increasing order is used to displace the haloes within PINOCCHIO, and how they depend on the grid resolution. Our approximate method is shown to be very fast when compared to full ray-tracing simulations from an N-Body run and able to recover the weak lensing signal, at different redshifts, with a few percent accuracy. It also allows for quickly constructing weak lensing covariance matrices, complementing PINOCCHIO's ability of generating the cluster mass function and galaxy clustering covariances and thus paving the way for calculating cross covariances between the different probes. This work advances these approximate methods as tools for simulating and analysing surveys data for cosmological purposes., Comment: Replaced to match the accepted version for publication in MNRAS. Added: comparison of diagonal and off-diagonal terms of the covariance matrices; cosmological constraints. Acceptance date: 2020-06-02

We compare the measurements of the bispectrum and the estimate of its covariance obtained from a set of different methods for the efficient generation of approximate dark matter halo catalogs to the same quantities obtained from full N-body simulations. To this purpose we employ a large set of three-hundred realisations of the same cosmology for each method, run with matching initial conditions in order to reduce the contribution of cosmic variance to the comparison. In addition, we compare how the error on cosmological parameters such as linear and nonlinear bias parameters depends on the approximate method used for the determination of the bispectrum variance. As general result, most methods provide errors within 10% of the errors estimated from N-body simulations. Exceptions are those methods requiring calibration of the clustering amplitude but restrict this to two-point statistics. Finally we test how our results are affected by being limited to a few hundreds measurements from N-body simulation, and therefore to the bispectrum variance, by comparing with a larger set of several thousands realisations performed with one approximate method., Additional results with respect to v1, new section and new figures added. 25 pages, 1 table 18 figures

We perform a suite of cosmological hydrodynamical simulations of disc galaxies, with zoomed-in initial conditions leading to the formation of a halo of mass $M_{\rm halo, \, DM} \simeq 2 \cdot 10^{12}$ M$_{\odot}$ at redshift $z=0$. These simulations aim at investigating the chemical evolution and the distribution of metals in a disc galaxy, and at quantifying the effect of $(i)$ the assumed IMF, $(ii)$ the adopted stellar yields, and $(iii)$ the impact of binary systems originating SNe Ia on the process of chemical enrichment. We consider either a Kroupa et al. (1993) or a more top-heavy Kroupa (2001) IMF, two sets of stellar yields and different values for the fraction of binary systems suitable to give rise to SNe Ia. We investigate stellar ages, SN rates, stellar and gas metallicity gradients, and stellar $\alpha$-enhancement in simulations, and compare predictions with observations. We find that a Kroupa et al. (1993) IMF has to be preferred when modelling late-type galaxies in the local universe. On the other hand, the comparison of stellar metallicity profiles and $\alpha$-enhancement trends with observations of Milky Way stars shows a better agreement when a Kroupa (2001) IMF is assumed. Comparing the predicted SN rates and stellar $\alpha$-enhancement with observations supports a value for the fraction of binary systems producing SNe Ia of $0.03$, at least for late-type galaxies and for the considered IMFs. Adopted stellar yields are crucial in regulating cooling and star formation, and in determining patterns of chemical enrichment for stars, especially for those located in the galaxy bulge., Comment: Accepted for publications in MNRAS. Comments welcome

This paper is the first in a set that analyses the covariance matrices of clustering statistics obtained from several approximate methods for gravitational structure formation. We focus here on the covariance matrices of anisotropic two-point correlation function measurements. Our comparison includes seven approximate methods, which can be divided into three categories: predictive methods that follow the evolution of the linear density field deterministically (ICE-COLA, Peak Patch, and Pinocchio), methods that require a calibration with N-body simulations (Patchy and Halogen), and simpler recipes based on assumptions regarding the shape of the probability distribution function (PDF) of density fluctuations (log-normal and Gaussian density fields). We analyse the impact of using covariance estimates obtained from these approximate methods on cosmological analyses of galaxy clustering measurements, using as a reference the covariances inferred from a set of full N-body simulations. We find that all approximate methods can accurately recover the mean parameter values inferred using the N-body covariances. The obtained parameter uncertainties typically agree with the corresponding N-body results within 5% for our lower mass threshold, and 10% for our higher mass threshold. Furthermore, we find that the constraints for some methods can differ by up to 20% depending on whether the halo samples used to define the covariance matrices are defined by matching the mass, number density, or clustering amplitude of the parent N-body samples. The results of our configuration-space analysis indicate that most approximate methods provide similar results, with no single method clearly outperforming the others., Comment: 23 pages, 11 figures. Replaced to match accepted MNRAS version. Included Kullback-Leibler divergence

During my Ph.D, my research activity has been mainly focused on the statistical analysis of the large-scale distribution of galaxies. The accurate determination of the cosmological parameters has become one of the key activities in modern cosmology. Choosing between different models requires accurate theoretical predictions of the observables and precise modeling of their statistical uncertainties. In the first part of the thesis I focused my attention on the study of systematic effects, such as galactic extinction, on the clustering of biased tracers of the cosmic density field, e.g. galaxies, that will be of capital importance for future surveys. I have built an analytic model to account for the effects of a generic foreground and we have measured the 2-point function in Fourier space and computed its covariance matrix using a large set (10’000) of synthetic realizations, run with the approximated method PINOCCHIO. I have analyzed all the corrections to the power spectrum and its covariance due to the presence of this foreground, including its coupling with the cosmological signal. This analysis has been published in a JCAP paper this years (Colavincenzo et al. 2017). Then I focused on the Euclid “Covariance comparison projects”, aimed at testing if covariance matrices of clustering measures, obtained with various approximated methods, are unbiased; I want to quantify their impact on the errors on cosmological parameters. I am in charge of the bispectrum covariance matrix comparison (Colavincenzo et al. in prep.). In the second part of this thesis I have used the large set of simulated galaxy catalogs I have access to test models of the galaxy power spectrum covariance matrix. In particular, I aim to get an analytic expression for the non-linear galaxy power spectrum covariance matrix (colavincenzo et al. in prep.). These same catalogs can be used to perform a likelihood analysis on the galaxy bias parameters using the bispectrum covariance matrix. The joint analysis of power spectrum and bispectrum can be used to improve the constraints on cosmological and galaxy bias parameters.

We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z=0, due mainly to the different implementations of feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas-rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disk is discernible in most simulations, though its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observation. Models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disk component. Our conclusions hold at two different numerical resolutions. Some differences can also be traced to the numerical techniques: more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high-angular momentum stellar disks remains a challenge., 28 pages, 16 figures, MNRAS accepted. High resolution figures and additional images can be found at http://www.aip.de/People/cscannapieco/aquila/

We present the latest version of Pinocchio, a code that generates catalogues of DM haloes in an approximate but fast way with respect to an N-body simulation. This code version extends the computation of particle and halo displacements up to 3rd-order Lagrangian Perturbation Theory (LPT), in contrast with previous versions that used Zeldovich approximation (ZA). We run Pinocchio on the same initial configuration of a reference N-body simulation, so that the comparison extends to the object-by-object level. We consider haloes at redshifts 0 and 1, using different LPT orders either for halo construction - where displacements are needed to decide particle accretion onto a halo or halo merging - or to compute halo final positions. We compare the clustering properties of Pinocchio haloes with those from the simulation by computing the power spectrum and 2-point correlation function (2PCF) in real and redshift space (monopole and quadrupole), the bispectrum and the phase difference of halo distributions. We find that 2LPT and 3LPT give noticeable improvement. 3LPT provides the best agreement with N-body when it is used to displace haloes, while 2LPT gives better results for constructing haloes. At the highest orders, linear bias is typically recovered at a few per cent level. In Fourier space and using 3LPT for halo displacements, the halo power spectrum is recovered to within 10 per cent up to $k_{max}\sim0.5\ h/$Mpc. The results presented in this paper have interesting implications for the generation of large ensemble of mock surveys aimed at accurately compute covariance matrices for clustering statistics., 20 pages, 20 figures, submitted to MNRAS

Recent inspections of local available data suggest that the almost linear relation between the stellar mass of spheroids ($M_{\rm sph}$) and the mass of the super massive Black Holes (BHs) residing at their centres, shows a break below $M_{\rm sph} \sim 10^{10}\ {\rm M}_\odot$, with a steeper, about quadratic relation at smaller masses. We investigate the physical mechanisms responsible for the change in slope of this relation, by comparing data with the results of the semi-analytic model of galaxy formation MORGANA, which already predicted such a break in its original formulation. We find that the change of slope is mostly induced by effective stellar feedback in star-forming bulges. The shape of the relation is instead quite insensitive to other physical mechanisms connected to BH accretion such as disc instabilities, galaxy mergers, Active Galactic Nucleus (AGN) feedback, or even the exact modelling of accretion onto the BH, direct or through a reservoir of low angular momentum gas. Our results support a scenario where most stars form in the disc component of galaxies and are carried to bulges through mergers and disc instabilities, while accretion onto BHs is connected to star formation in the spheroidal component. Therefore, a model of stellar feedback that produces stronger outflows in star-forming bulges than in discs will naturally produce a break in the scaling relation. Our results point to a form of co-evolution especially at lower masses, below the putative break, mainly driven by stellar feedback rather than AGN feedback., MNRAS accepted, 10 pages, 6 figures, 1 table

We present predictions of Spectral Energy Distributions (SEDs), from the UV to the FIR, of simulated galaxies at $z=0$. These were obtained by post-processing the results of an N-body+hydro simulation of a small cosmological volume, that uses the Multi-Phase Particle Integrator (MUPPI) for star formation and stellar feedback, with the GRASIL-3D radiative transfer code, that includes reprocessing of UV light by dust. Physical properties of galaxies resemble observed ones, though with some tension at small and large stellar masses. Comparing predicted SEDs of simulated galaxies with different samples of local galaxies, we find that these resemble observed ones, when normalised at 3.6 $\mu$m. A comparison with the Herschel Reference Survey shows that, when binning galaxies in Star Formation Rate (SFR), average SEDs are reproduced to within a factor of $\sim2$ even in normalization, while binning in stellar mass highlights the same tension that is present in the stellar mass -- SFR plane. We use our sample to investigate the correlation of IR luminosity in Spitzer and Herschel bands with several galaxy properties. SFR is the quantity that best correlates with IR light up to $160\ \mu$m, while at longer wavelengths better correlations are found with molecular mass and, at $500\ \mu$m, with dust mass. However, using the position of the FIR peak as a proxy for cold dust temperature, we assess that heating of cold dust is mostly determined by SFR, with stellar mass giving only a minor contribution. We finally show how our sample of simulated galaxies can be used as a guide to understand the physical properties and selection biases of observed samples., Comment: 18 pages, 11 figures, Submitted to MNRAS

Precision cosmology has recently triggered new attention on the topic of approximate methods for the clustering of matter on large scales, whose foundations date back to the period from late '60s to early '90s. Indeed, although the prospect of reaching sub-percent accuracy in the measurement of clustering poses a challenge even to full N-body simulations, an accurate estimation of the covariance matrix of clustering statistics, not to mention the sampling of parameter space, requires usage of a large number (hundreds in the most favourable cases) of simulated (mock) galaxy catalogs. Combination of few N-body simulations with a large number of realizations performed with approximate methods gives the most promising approach to solve these problems with a reasonable amount of resources. In this paper I review this topic, starting from the foundations of the methods, then going through the pioneering efforts of the '90s, and finally presenting the latest extensions and a few codes that are now being used in present-generation surveys and thoroughly tested to assess their performance in the context of future surveys., Revised version of the paper submitted to Galaxies, special issue: "Dark Matter: Large versus Small Scale Structures", ed. J. Gaite & A. Diaferio. 41 pages. Comments welcome

We use a sample of 1669 QSOs ($r, Comment: 8 pages, 9 figures, MNRAS accepted

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society ©: 2015 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved, We present a comparison of major methodologies of fast generating mock halo or galaxy catalogues. The comparison is done for two-point (power spectrum and two-point correlation function in real and redshift space), and the three-point clustering statistics (bispectrum and three-point correlation function). The reference catalogues are drawn from the BigMultiDark N-body simulation. Both friend-of-friends (including distinct haloes only) and spherical overdensity (including distinct haloes and subhalos) catalogues have been used with the typical number density of a large volume galaxy surveys. We demonstrate that a proper biasing model is essential for reproducing the power spectrum at quasi-linear and even smaller scales. With respect to various clustering statistics, a methodology based on perturbation theory and a realistic biasing model leads to very good agreement with N-body simulations. However, for the quadrupole of the correlation function or the power spectrum, only the method based on semi-N-body simulation could reach high accuracy (1 per cent level) at small scales, i.e. r < 25 h-1 Mpc or k > 0.15 h Mpc-1. Full N-body solutions will remain indispensable to produce reference catalogues. Nevertheless, we have demonstrated that the more efficient approximate solvers can reach a few per cent accuracy in terms of clustering statistics at the scales interesting for the large-scale structure analysis. This makes them useful for massive production aimed at covariance studies, to scan large parameter spaces, and to estimate uncertainties in data analysis techniques, such as baryon acoustic oscillation reconstruction, redshift distortion measurements, etc, The authors would like to express special thanks to the Instituto de Fisica Teorica (IFT-UAM/CSIC in Madrid) for its hospitality and support, via the Centro de Excelencia Severo Ochoa Program under Grant No. SEV-2012-0249, during the three-weekworkshop ‘nIFTy Cosmology’ where this work developed. We further acknowledge the financial support of the University of Western 2014 Australia Research Collaboration Award for ‘Fast Approximate Synthetic Universes for the SKA’, the ARC Centre of Excellence for All Sky Astrophysics (CAASTRO) grant number CE110001020, and the two ARC Discovery Projects DP130100117 and DP140100198. CC and FP were supported by the Spanish MICINNs Consolider- Ingenio 2010 Programme under grant MultiDark CSD2009-00064 and AYA2010-21231-C02-01 grant, and Spanish MINECOs Centro de Excelencia Severo Ochoa Programme under grant SEV-2012- 0249. CZ and CT acknowledges support from Tsinghua University, and 973 student programme No. 2013CB834906. CZ also thanks the support from the MultiDark summer student program to visit the Instituto de Física Teórica, (UAM/CSIC), Spain. AK is supported by the Ministerio de Economía y Competitividad (MINECO) in Spain through grant AYA2012-31101, as well as the Consolider- Ingenio 2010 Programme of the Spanish Ministerio de Ciencia e Innovación (MICINN) under grant MultiDark CSD2009-00064. He also acknowledges support from the Australian Research Council (ARC) grants DP130100117 and DP140100198. EM and PM have been supported by a FRA2012 grant of the University of Trieste, PRIN2010-2011 (J91J12000450001) from MIUR, and by Consorzio per la Fisica di Trieste. AI is supported by the JAE program grant from the Spanish National Science Council (CSIC). GY acknowledges support from the Spanish MINECO under research grants AYA2012-31101, FPA2012-34694, AYA2010-21231, Consolider Ingenio SyeC CSD2007-0050 and from Comunidad de Madrid under ASTROMADRID project (S2009/ESP-1496). VT and SG have been supported by the German Science Foundation under DFG grant GO563/23-1. FAM is supported by the Australian Research Council Centre of Excellence for All-Sky Astrophysics (CAASTRO) through project number CE110001020

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., La lista completa de autores que integran el documento puede consultarse en el archivo., Instituto de Astrofísica de La Plata, Facultad de Ciencias Astronómicas y Geofísicas

A well calibrated method to describe the environment of galaxies at all redshifts is essential for the study of structure formation. Such a calibration should include well understood correlations with halo mass, and the possibility to identify galaxies which dominate their potential well (centrals), and their satellites. Focusing on z = 1 and 2 we propose a method of environmental calibration which can be applied to the next generation of low to medium resolution spectroscopic surveys. Using an up-to-date semi-analytic model of galaxy formation, we measure the local density of galaxies in fixed apertures on different scales. There is a clear correlation of density with halo mass for satellite galaxies, while a significant population of low mass centrals is found at high densities in the neighbourhood of massive haloes. In this case the density simply traces the mass of the most massive halo within the aperture. To identify central and satellite galaxies, we apply an observationally motivated stellar mass rank method which is both highly pure and complete, especially in the more massive haloes where such a division is most meaningful. Finally we examine a test case for the recovery of environmental trends: the passive fraction of galaxies and its dependence on stellar and halo mass for centrals and satellites. With careful calibration, observationally defined quantities do a good job of recovering known trends in the model. This result stands even with reduced redshift accuracy, provided the sample is deep enough to preserve a wide dynamic range of density., Comment: 19 pages, 12 figures, accepted for publication in MNRAS