We investigate the accuracy of weak lensing simulations by comparing the results of five independently developed lensing simulation codes run on the same input $N$-body simulation. Our comparison focuses on the lensing convergence maps produced by the codes, and in particular on the corresponding PDFs, power spectra and peak counts. We find that the convergence power spectra of the lensing codes agree to $\lesssim 2\%$ out to scales $\ell \approx 4000$. For lensing peak counts, the agreement is better than $5\%$ for peaks with signal-to-noise $\lesssim 6$. We also discuss the systematic errors due to the Born approximation, line-of-sight discretization, particle noise and smoothing. The lensing codes tested deal in markedly different ways with these effects, but they nonetheless display a satisfactory level of agreement. Our results thus suggest that systematic errors due to the operation of existing lensing codes should be small. Moreover their impact on the convergence power spectra for a lensing simulation can be predicted given its numerical details, which may then serve as a validation test.

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.

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.

The forthcoming generation of galaxy redshift surveys will sample the large-scale structure of the Universe over unprecedented volumes with high-density tracers. This advancement will make robust measurements of three-point clustering statistics possible. In preparation for this improvement, we investigate how several methodological choices can influence inferences based on the bispectrum about galaxy bias and shot noise. We first measure the real-space bispectrum of dark-matter haloes extracted from 298 N-body simulations covering a volume of approximately $1000 h^{-3} \mathrm{Gpc}^3$. We then fit a series of theoretical models based on tree-level perturbation theory to the numerical data. To achieve this, we estimate the covariance matrix of the measurement errors by using 10,000 mock catalogues generated with the Pinocchio code. We study how the model constraints are influenced by the binning strategy for the bispectrum configurations and by the form of the likelihood function. We also use Bayesian model-selection techniques to single out the optimal theoretical description of our data. We find that a three-parameter bias model combined with Poissonian shot noise is necessary to model the halo bispectrum up to scales of $k_\mathrm{max}\lesssim 0.08 h \mathrm{Mpc}^{-1}$, although fitting formulae that relate the bias parameters can be helpful to reduce the freedom of the model without compromising accuracy. Our data clearly disfavour local Eulerian and local Lagrangian bias models and do not require corrections to Poissonian shot noise. We anticipate that model-selection diagnostics will be particularly useful to extend the analysis to smaller scales as, in this case, the number of model parameters will grow significantly large.

We present simulations of galaxy formation, based on the GADGET-3 code, in which a sub-resolution model for star formation and stellar feedback is interfaced with a new model for AGN feedback. Our sub-resolution model describes a multiphase ISM, accounting for hot and cold gas within the same resolution element: we exploit this feature to investigate the impact of coupling AGN feedback energy to the different phases of the ISM over cosmic time. Our fiducial model considers that AGN feedback energy coupling is driven by the covering factors of the hot and cold phases. We perform a suite of cosmological hydrodynamical simulations of disc galaxies ($M_{\rm halo, \, DM} \simeq 2 \cdot 10^{12}$ M$_{\odot}$, at $z=0$), to investigate: $(i)$ the effect of different ways of coupling AGN feedback energy to the multiphase ISM; $(ii)$ the impact of different prescriptions for gas accretion (i.e. only cold gas, both cold and hot gas, with the additional possibility of limiting gas accretion from cold gas with high angular momentum); $(iii)$ how different models of gas accretion and coupling of AGN feedback energy affect the coevolution of supermassive BHs and their host galaxy. We find that at least a share of the AGN feedback energy has to couple with the diffuse gas, in order to avoid an excessive growth of the BH mass. When the BH only accretes cold gas, it experiences a growth that is faster than in the case in which both cold and hot gas are accreted. If the accretion of cold gas with high angular momentum is reduced, the BH mass growth is delayed, the BH mass at $z=0$ is reduced by up to an order of magnitude, and the BH is prevented from accreting below $z \lesssim 2$, when the galaxy disc forms.

We present and test a new numerical method to determine second-order Lagrangian displacement fields in the context of modified gravity (MG) theories. We start from the extension of Lagrangian Perturbation Theory to a class of MG models that can be described by a parametrized Poisson equation, with the introduction of a scale-dependent function. We exploit fast Fourier transforms to compute the full source term of the differential equation for the second-order Lagrangian displacement field. We compare its mean to the source term computed for specific configurations for which a k-dependent solution can be found numerically. We choose the configuration that best matches the full source term, thus obtaining an approximate factorization of the second-order displacement field as the space term valid for standard gravity times a k-dependent, second-order growth factor $D_2(k,t)$. This approximation is used to compute second order displacements for particles. The method is tested against N-body simulations run with standard and $f(R)$ gravity: we rely on the results of a friends-of-friends code run on the N-body snapshots to assign particles to halos, then compute the halo power spectrum. We find very consistent results for the two gravity theories: second-order LPT (2LPT) allows to recover the halo power spectrum of N-body simulations within $\sim 10\%$ precision to $k\sim 0.2-0.4\ h\ {\rm Mpc}^{-1}$, as well as halo positions, with an error that is a fraction of the inter-particle distance. We show that, when considering the same level of non-linearity in the density field, the performance of 2LPT with MG is the same (within $1\%$) as the one obtained for the standard $\Lambda$CDM model with General Relativity. When implemented in a computer code, this formulation of 2LPT can quickly generate dark matter distributions with $f(R)$ gravity, and can easily be extended to other MG theories.

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.

We present a blind method to determine the properties of a foreground contamination, given by a visibility mask, that affects a deep galaxy survey. Angular cross correlations of density fields in different redshift bins are expected to vanish (apart from a contribution due to lensing), but are sensitive to the presence of a foreground that modulates the flux limit across the sky. After formalizing the expected effect of a foreground mask on the measured galaxy density, under a linear, luminosity-dependent bias model for galaxies, we construct two estimators that single out the mask contribution if a sufficient number of independent redshift bins is available. These estimators are combined to give a reconstruction of the mask. We use Milky-Way reddening as a prototype for the mask. Using a set of 20 large mock catalogs covering $1/4$-th of the sky and number-matched to $H\alpha$ emitters to mimic an Euclid-like sample, we demonstrate that our method can reconstruct the mask and its angular clustering at scales $\ell<100$, beyond which the cosmological signal becomes dominant. The uncertainty of this reconstruction is quantified to be $1/3$-rd of the sample variance of the signal. Such a reconstruction requires knowledge of the average and square average of the mask, but we show that it is possible to recover this information either from external models or internally from the data. It also relies on knowledge of how the impact of the foreground changes with redshift (due to the extinction curve in our case), but this can be tightly constrained by cross correlations of different redshift bins. The strong points of this blind reconstruction technique lies in the ability to find "unknown unknowns" that affect a survey, and in the facility to quantify, using sets of mock catalogs, how its uncertainty propagates to clustering measurements. [Abridged]

We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constrains related to growth rate of structure, Alcock-Paczynski distortions and biasing. We consider seven methods in three broad categories: algorithms that solve for halo density evolution deterministically using Lagrangian trajectories (ICE-COLA, Pinocchio and PeakPatch), methods that rely on halo assignment schemes onto dark-matter overdensities calibrated with a target N-body run (Halogen, Patchy) and two standard assumptions about the full density PDF (Gaussian and Lognormal). We benchmark their performance against a set of three hundred N-body simulations, running similar sets of approximate simulations with matched initial conditions, for each method. We find that most methods reproduce the monopole to within $5\%$, while residuals for the quadrupole are sometimes larger and scale dependent. The variance of the multipoles is typically reproduced within $10\%$. Overall, we find that covariances built from approximate simulations yield errors on model parameters within $10\%$ of those from the N-body based covariance.

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.

We calculate cosmic distributions in space and time of the formation sites of the first, "Pop III.1" stars, exploring a model in which these are the progenitors of all supermassive black holes (SMBHs), seen in the centers of most large galaxies. Pop III.1 stars are defined to form from primordial composition gas in dark matter minihalos with $\sim10^6\:M_\odot$ that are isolated from neighboring astrophysical sources by a given isolation distance, $d_{\rm{iso}}$. We assume Pop III.1 sources are seeds of SMBHs, based on protostellar support by dark matter annihilation heating that allows them to accrete a large fraction of their minihalo gas, i.e., $\sim10^5\:M_\odot$. Exploring $d_{\rm{iso}}$ from $10 - 100\:\rm{kpc}$ (proper distances), we predict the redshift evolution of Pop III.1 source and SMBH remnant number densities. The local, $z=0$ density of SMBHs constrains $d_{\rm{iso}}\lesssim 100\:\rm{kpc}$ (i.e., $3\:\rm{Mpc}$ comoving distance at $z\simeq30$). In our simulated ($\sim60\:\rm{Mpc}$)$^3$ comoving volume, Pop III.1 stars start forming just after $z=40$. Their formation is largely complete by $z\simeq25$ to $20$ for $d_{\rm{iso}}=100$ to $50\:\rm{kpc}$. We follow source evolution to $z=10$, by which point most SMBHs reside in halos with $\gtrsim10^8\:M_\odot$. Over this period, there is relatively limited merging of SMBHs for these values of $d_{\rm{iso}}$. We also predict SMBH clustering properties at $z=10$: feedback suppression of neighboring sources leads to relatively flat angular correlation functions.

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.

We characterize the properties of the intergalactic medium (IGM) around a sample of galaxies extracted from state-of-the-art hydrodynamical simulations of structure formation in a cosmological volume of 25 Mpc comoving at $z\sim 2$. The simulations are based on two different subresolution schemes for star formation and supernova feedback: the MUlti-Phase Particle Integrator (MUPPI) scheme and the Effective Model. We develop a quantitative and probabilistic analysis based on the apparent optical depth method of the properties of the absorbers as a function of impact parameter from their nearby galaxies: in such a way we probe different environments from circumgalactic medium to low-density filaments. Absorbers' properties are then compared with a spectroscopic observational data set obtained from high-resolution quasar spectra. Our main focus is on the N$_{\rm{ CIV}}$-N$_{\rm { HI}}$ relation around simulated galaxies: the results obtained with MUPPI and the Effective model are remarkably similar, with small differences only confined to regions at impact parameters $b = [1-3] \times r_{\rm {vir}}$. Using $\mbox{C IV}$ as a tracer of the metallicity, we obtain evidence that the observed metal absorption systems have the highest probability to be confined in a region of 150-400 kpc around galaxies. Near-filament environments have instead metallicities too low to be probed by present-day telescopes, but could be probed by future spectroscopical studies. Finally we compute $\mbox{C IV}$ covering fractions which are in agreement with observational data.

We investigate the impact of galactic outflow modelling on the formation and evolution of a disc galaxy, by performing a suite of cosmological simulations with zoomed-in initial conditions of a Milky Way-sized halo. We verify how sensitive the general properties of the simulated galaxy are to the way in which stellar feedback triggered outflows are implemented, keeping initial conditions, simulation code and star formation (SF) model all fixed. We present simulations that are based on a version of the GADGET3 code where our sub-resolution model is coupled with an advanced implementation of Smoothed Particle Hydrodynamics that ensures a more accurate fluid sampling and an improved description of gas mixing and hydrodynamical instabilities. We quantify the strong interplay between the adopted hydrodynamic scheme and the sub-resolution model describing SF and feedback. We consider four different galactic outflow models, including the one introduced by Dalla Vecchia and Schaye (2012) and a scheme that is inspired by the Springel and Hernquist (2003) model. We find that the sub-resolution prescriptions adopted to generate galactic outflows are the main shaping factor of the stellar disc component at low redshift. The key requirement that a feedback model must have to be successful in producing a disc-dominated galaxy is the ability to regulate the high-redshift SF (responsible for the formation of the bulge component), the cosmological infall of gas from the large-scale environment, and gas fall-back within the galactic radius at low redshift, in order to avoid a too high SF rate at $z=0$.

Our current knowledge of star formation and accretion luminosity at high-redshift (z>3-4), as well as the possible connections between them, relies mostly on observations in the rest-frame ultraviolet (UV), which are strongly affected by dust obscuration. Due to the lack of sensitivity of past and current infrared (IR) instrumentation, so far it has not been possible to get a glimpse into the early phases of the dust-obscured Universe. Among the next generation of IR observatories, SPICA, observing in the 12-350 micron range, will be the only facility that can enable us to make the required leap forward in understanding the obscured star-formation rate and black-hole accretion rate densities (SFRD and BHARD, respectively) with respect to what Spitzer and Herschel achieved in the mid- and far-IR at z<3. In particular, SPICA will have the unique ability to trace the evolution of the obscured SFRD and BHARD over cosmic time, from the peak of their activity back to the reionisation epoch (i.e., 3<z<6-7), where its predecessors had severe limitations. Here we discuss the potential of both deep and shallow photometric surveys performed with the SPICA mid-IR instrument (SMI), enabled by the very low level of impact of dust obscuration in a band centred at 34 micron. These unique unbiased photometric surveys that SPICA will perform will be followed up by observations both with the SPICA spectrometers and with other facilities at shorter and longer wavelengths, with the aim to fully characterise the evolution of AGNs and star-forming galaxies after re-ionisation.

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 {\Lambda}CDM 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 2-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 Mpc/h. 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 test to quantify how well some approximate methods, designed to reproduce the mildly non-linear evolution of perturbations, are able to reproduce the clustering of DM halos once the grouping of particles into halos is defined and kept fixed. The following methods have been considered: Lagrangian Perturbation Theory (LPT) up to third order, Truncated LPT, Augmented LPT, MUSCLE and COLA. The test runs as follows: halos are defined by applying a friends-of-friends (FoF) halo finder to the output of an N-body simulation. The approximate methods are then applied to the same initial conditions of the simulation, producing for all particles displacements from their starting position and velocities. The position and velocity of each halo are computed by averaging over the particles that belong to that halo, according to the FoF halo finder. This procedure allows us to perform a well-posed test of how clustering of the matter density and halo density fields are recovered, without asking to the approximate method an accurate reconstruction of halos. We have considered the results at $z=0,0.5,1$, and we have analysed power spectrum in real and redshift space, object-by-object difference in position and velocity, density Probability Distribution Function (PDF) and its moments, phase difference of Fourier modes. We find that higher LPT orders are generally able to better reproduce the clustering of halos, while little or no improvement is found for the matter density field when going to 2LPT and 3LPT. Augmentation provides some improvement when coupled with 2LPT, while its effect is limited when coupled with 3LPT. Little improvement is brought by MUSCLE with respect to Augmentation. The more expensive particle-mesh code COLA outperforms all LPT methods [abridged]

The forecasted accuracy of upcoming surveys of large-scale structure cannot be achieved without a proper quantification of the error induced by foreground removal (or other systematics like 0-point photometry offset). Because these errors are highly correlated on the sky, their influence is expected to be especially important at very large scales. In this work we quantify how the uncertainty in the visibility mask of a survey influences the measured power spectrum of a sample of tracers of the density field and its covariance matrix. We start from a very large set of 10,000 catalogs of dark matter halos in periodic cosmological boxes, produced with the PINOCCHIO approximate method. To make an analytic approach feasible, we assume luminosity-independent halo bias and an idealized geometry for the visibility mask. We find that the power spectrum of these biased tracers can be expressed as the sum of a cosmological term, a mask term and a term involving their convolution. The mask and convolution terms scale like $P\propto l^2\sigma_A^2$, where $\sigma_A^2$ is the variance of the uncertainty on the visibility mask. With $l=30-100$ Mpc$/h$ and $\sigma_A=5-20$\%, the mask term can be significant at $k\sim0.01-0.1\ h/$Mpc, and the convolution term can amount to $\sim 1-10$\% of the total. For the power spectrum covariance, the coupling of the convolution term with the other two gives rise to several mixed terms, that we quantify by difference using the mock catalogs. These are found to be of the same order of the mask covariance, and to introduce non-diagonal terms at large scales. Then, the power spectrum covariance matrix cannot be expressed as the sum of a cosmological and of a mask term. Our results lie down the theoretical bases to quantify the impact that uncertainties in the mask calibration have on the derivation of cosmological constraints from large spectroscopic surveys. [Abridged]

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.

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.

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.

We use a sample of 1669 QSOs ($r<20.15$, $3.6<z<4.0$) from the BOSS survey to study the intrinsic shape of their continuum and the Lyman continuum photon escape fraction (f$_{esc}$), estimated as the ratio between the observed flux and the expected intrinsic flux (corrected for the intergalactic medium absorption) in the wavelength range 865-885 \AA\ rest-frame. Modelling the intrinsic QSO continuum shape with a power-law, $F_{\lambda}\propto\lambda^{-\gamma}$, we find a median $\gamma=1.30$ (with a dispersion of $0.38$, no dependence on the redshift and a mild intrinsic luminosity dependence) and a mean f$_{esc}=0.75$ (independent of the QSO luminosity and/or redshift). The f$_{esc}$ distribution shows a peak around zero and a long tail of higher values, with a resulting dispersion of $0.7$. If we assume for the QSO continuum a double power-law shape (also compatible with the data) with a break located at $\lambda_{\rm br}=1000$ \AA\ and a softening $\Delta\gamma=0.72 $ at wavelengths shorter than $\lambda_{\rm br}$, the mean f$_{esc}$ rises to $=0.82$. Combining our $\gamma$ and f$_{esc}$ estimates with the observed evolution of the AGN luminosity function (LF) we compute the AGN contribution to the UV ionizing background (UVB) as a function of redshift. AGN brighter than one tenth of the characteristic luminosity of the LF are able to produce most of it up $z\sim 3$, if the present sample is representative of their properties. At higher redshifts a contribution of the galaxy population is required. Assuming an escape fraction of Lyman continuum photons from galaxies between $5.5$ and $7.6\%$, independent of the galaxy luminosity and/or redshift, a remarkably good fit to the observational UVB data up to $z\sim 6$ is obtained. At lower redshift the extrapolation of our empirical estimate agrees well with recent UVB observations, dispelling the so-called Photon Underproduction Crisis.

We implement novel numerical models of AGN feedback in the SPH code GADGET-3, where the energy from a supermassive black hole (BH) is coupled to the surrounding gas in the kinetic form. Gas particles lying inside a bi-conical volume around the BH are imparted a one-time velocity (10,000 km/s) increment. We perform hydrodynamical simulations of isolated cluster (total mass 10^14 /h M_sun), which is initially evolved to form a dense cool core, having central T<10^6 K. A BH resides at the cluster center, and ejects energy. The feedback-driven fast wind undergoes shock with the slower-moving gas, which causes the imparted kinetic energy to be thermalized. Bipolar bubble-like outflows form propagating radially outward to a distance of a few 100 kpc. The radial profiles of median gas properties are influenced by BH feedback in the inner regions (r<20-50 kpc). BH kinetic feedback, with a large value of the feedback efficiency, depletes the inner cool gas and reduces the hot gas content, such that the initial cool core of the cluster is heated up within a time 1.9 Gyr, whereby the core median temperature rises to above 10^7 K, and the central entropy flattens. Our implementation of BH thermal feedback (using the same efficiency as kinetic), within the star-formation model, cannot do this heating, where the cool core remains. The inclusion of cold gas accretion in the simulations produces naturally a duty cycle of the AGN with a periodicity of 100 Myr.

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.

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.

We present a comparison of major methodologies of fast generating mock halo or galaxy catalogues. The comparison is done for two-point and the three-point clustering statistics. The reference catalogues are drawn from the BigMultiDark N-body simulation. Both friend-of-friends (including distinct halos only) and spherical overdensity (including distinct halos and subhalos) catalogs 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 quasilinear 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% level) at small scales, i.e., r<25 Mpc/h or k>0.15 h/Mpc. Full N-body solutions will remain indispensable to produce reference catalogues. Nevertheless, we have demonstrated that the far more efficient approximate solvers can reach a few percent accuracy in terms of clustering statistics at the scales interesting for the large-scale structure analysis after calibration with a few reference N-body calculations. 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.

We present a direct comparison between the observed star formation rate functions (SFRF) and the state-of-the-art predictions of semi-analytic models (SAM) of galaxy formation and evolution. We use the PACS Evolutionary Probe Survey (PEP) and Herschel Multi-tiered Extragalactic Survey (HerMES) data-sets in the COSMOS and GOODS-South fields, combined with broad-band photometry from UV to sub-mm, to obtain total (IR+UV) instantaneous star formation rates (SFRs) for individual Herschel galaxies up to z~4, subtracted of possible active galactic nucleus (AGN) contamination. The comparison with model predictions shows that SAMs broadly reproduce the observed SFRFs up to z~2, when the observational errors on the SFR are taken into account. However, all the models seem to under-predict the bright-end of the SFRF at z>2. The cause of this underprediction could lie in an improper modelling of several model ingredients, like too strong (AGN or stellar) feedback in the brighter objects or too low fall-back of gas, caused by weak feedback and outflows at earlier epochs.

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.

In the effort to understand the link between the structure of galaxy clusters and their galaxy populations, we focus on MACSJ1206.2-0847 at z~0.44 and probe its substructure in the projected phase space through the spectrophotometric properties of a large number of galaxies from the CLASH-VLT survey. Our analysis is mainly based on an extensive spectroscopic dataset of 445 member galaxies, mostly acquired with VIMOS@VLT as part of our ESO Large Programme, sampling the cluster out to a radius ~2R200 (4 Mpc). We classify 412 galaxies as passive, with strong Hdelta absorption (red and blue galaxies, and with emission lines from weak to very strong. A number of tests for substructure detection are applied to analyze the galaxy distribution in the velocity space, in 2D space, and in 3D projected phase-space. Studied in its entirety, the cluster appears as a large-scale relaxed system with a few secondary, minor overdensities in 2D distribution. We detect no velocity gradients or evidence of deviations in local mean velocities. The main feature is the WNW-ESE elongation. The analysis of galaxy populations per spectral class highlights a more complex scenario. The passive galaxies and red strong Hdelta galaxies trace the cluster center and the WNW-ESE elongated structure. The red strong Hdelta galaxies also mark a secondary, dense peak ~2 Mpc at ESE. The emission line galaxies cluster in several loose structures, mostly outside R200. The observational scenario agrees with MACS J1206.2-0847 having WNW-ESE as the direction of the main cluster accretion, traced by passive galaxies and red strong Hdelta galaxies. The red strong Hdelta galaxies, interpreted as poststarburst galaxies, date a likely important event 1-2 Gyr before the epoch of observation. The emission line galaxies trace a secondary, ongoing infall where groups are accreted along several directions.

We present a quantification of the properties of bars in two N-body+SPH cosmological simulations of spiral galaxies, named GA and AqC. The initial conditions were obtained using the zoom-in technique and represent two dark matter (DM) halos of $2-3\times10^{12}\ {\rm M}_\odot$, available at two different resolutions. The resulting galaxies are presented in the companion paper of Murante et al. (2014). We find that the GA galaxy has a bar of length $8.8$ kpc, present at the two resolution levels even though with a slightly different strength. Classical bar signatures (e.g. pattern of streaming motions, high $m=2$ Fourier mode with roughly constant phase) are consistently found at both resolutions. Though a close encounter with a merging satellite at $z\sim0.6$ (mass ratio $1:50$) causes a strong, transient spiral pattern and some heating of the disk, we find that bar instability is due to secular process, caused by a low Toomre parameter $Q\lesssim1$ due to accumulation of mass in the disk. The AqC galaxy has a slightly different history: it suffers a similar tidal disturbance due to a merging satellite at $z\sim0.5$ but with a mass ratio of $1:32$, that triggers a bar in the high-resolution simulation, while at low resolution the merging is found to take place at a later time, so that both secular evolution and merging are plausible triggers for bar instability.

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.

We measure and quantify properties of galactic outflows and diffuse gas at $z \geq 1$ in cosmological hydrodynamical simulations. Our novel sub-resolution model, MUPPI, implements supernova feedback using fully local gas properties, where the wind velocity and mass loading are not given as input. We find the following trends at $z = 2$ by analysing central galaxies having a stellar mass higher than $10^{9} M_{\odot}$. The outflow velocity and mass outflow rate ($\dot{M}_{\rm out}$) exhibit positive correlations with galaxy mass and with the star formation rate (SFR). However, most of the relations present a large scatter. The outflow mass loading factor ($\eta$) is between $0.2 - 10$. The comparison Effective model generates a constant outflow velocity, and a negative correlation of $\eta$ with halo mass. The number fraction of galaxies where outflow is detected decreases at lower redshifts, but remains more than $80 \%$ over $z = 1 - 5$. High SF activity at $z \sim 2 - 4$ drives strong outflows, causing the positive and steep correlations of velocity and $\dot{M}_{\rm out}$ with SFR. The outflow velocity correlation with SFR becomes flatter at $z = 1$, and $\eta$ displays a negative correlation with halo mass in massive galaxies. Our study demonstrates that both the MUPPI and Effective models produce significant outflows at $\sim 1 / 10$ of the virial radius; at the same time shows that the properties of outflows generated can be different from the input speed and mass loading in the Effective model. Our MUPPI model, using local properties of gas in the sub-resolution recipe, is able to develop galactic outflows whose properties correlate with global galaxy properties, and consistent with observations.

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.

We compare predictions of cooled masses and cooling rates from three stripped-down Semi-Analytic Models (SAMs) of galaxy formation with the results of N-body+SPH simulations with gas particle mass of 3.9x10^6 Msun, where radiative cooling of a gas of primordial composition is implemented. We also run a simulation where cooling is switched on at redshift ~2, in order to test cooling models in a regime in which their approximations are expected to be valid. We confirm that cooling models implemented in SAMs are able to predict the amount of cooled mass at z=0 to within ~20 per cent. However, some relevant discrepancies are found. (i) When the contribution from poorly resolved halos is subtracted out, SAMs tend to under-predict by ~30 per cent the mass that cools in the infall-dominated regime. (ii) At large halo masses SAMs tend to over-predict cooling rates, though the numerical result may be affected by the use of SPH. (iii) As found in our previous work, cooling rates are found to be significantly affected by model details: simulations disfavour models with large cores and with quenching of cooling at major mergers. (iv) When cooling is switched on at z~2, cold gas accumulates very quickly in the simulated halos. This accumulation is reproduced by SAMs with varying degrees of accuracy.

We compare state-of-the-art semi-analytic models of galaxy formation as well as advanced sub-halo abundance matching models with a large sample of early-type galaxies from SDSS at z < 0.3. We focus our attention on the dependence of median sizes of central galaxies on host halo mass. The data do not show any difference in the structural properties of early-type galaxies with environment, at fixed stellar mass. All hierarchical models considered in this work instead tend to predict a moderate to strong environmental dependence, with the median size increasing by a factor of about 1.5-3 when moving from low to high mass host haloes. At face value the discrepancy with the data is highly significant, especially at the cluster scale, for haloes above log Mhalo > 14. The convolution with (correlated) observational errors reduces some of the tension. Despite the observational uncertainties, the data tend to disfavour hierarchical models characterized by a relevant contribution of disc instabilities to the formation of spheroids, strong gas dissipation in (major) mergers, short dynamical friction timescales, and very short quenching timescales in infalling satellites. We also discuss a variety of additional related issues, such as the slope and scatter in the local size-stellar mass relation, the fraction of gas in local early-type galaxies, and the general predictions on satellite galaxies.

[Abridged] In this paper, we carry out a detailed analysis of the performance of two different methods to identify the diffuse stellar light in cosmological hydrodynamical simulations of galaxy clusters. One method is based on a dynamical analysis of the stellar component. The second method is closer to techniques commonly employed in observational studies. Both the dynamical method and the method based on the surface brightness limit criterion are applied to the same set of hydrodynamical simulations for a large sample about 80 galaxy clusters. We find significant differences between the ICL and DSC fractions computed with the two corresponding methods, which amounts to about a factor of two for the AGN simulations, and a factor of four for the CSF set. We also find that the inclusion of AGN feedback boosts the DSC and ICL fractions by a factor of 1.5-2, respectively, while leaving the BCG+ICL and BCG+DSC mass fraction almost unchanged. The sum of the BCG and DSC mass stellar mass fraction is found to decrease from ~80 per cent in galaxy groups to ~60 per cent in rich clusters, thus in excess of what found from observational analysis. We identify the average surface brightness limits that yields the ICL fraction from the SBL method close to the DSC fraction from the dynamical method. These surface brightness limits turn out to be brighter in the CSF than in the AGN simulations. This is consistent with the finding that AGN feedback makes BCGs to be less massive and with shallower density profiles than in the CSF simulations. The BCG stellar component, as identified by both methods, are slightly older and more metal-rich than the stars in the diffuse component.

We explore the scale dependence of halo bias using real space cross-correlation measurements in N-body simulations and in Pinocchio, an algorithm based on Lagrangian Perturbation Theory. Recent work has shown how to interpret such real space measurements in terms of k-dependent bias in Fourier space, and how to remove the k-dependence to reconstruct the k-independent peak-background split halo bias parameters. We compare our reconstruction of the linear bias, which requires no free parameters, with previous estimates from N-body simulations which were obtained directly in Fourier space at large scales, and find very good agreement. Our reconstruction of the quadratic bias is similarly parameter-free, although in this case there are no previous Fourier space measurements to compare with. Our analysis of N-body simulations explicitly tests the predictions of the excursion set peaks (ESP) formalism of Paranjape et al. (2013) for the scale dependence of bias; we find that the ESP predictions accurately describe our measurements. In addition, our measurements in Pinocchio serve as a useful, successful consistency check between Pinocchio and N-body simulations that is not accessible to traditional measurements.

We exploit the deep and extended far infrared data sets (at 70, 100 and 160 um) of the Herschel GTO PACS Evolutionary Probe (PEP) Survey, in combination with the HERschel Multi tiered Extragalactic Survey (HerMES) data at 250, 350 and 500 um, to derive the evolution of the restframe 35 um, 60 um, 90 um, and total infrared (IR) luminosity functions (LFs) up to z~4. We detect very strong luminosity evolution for the total IR LF combined with a density evolution. In agreement with previous findings, the IR luminosity density increases steeply to z~1, then flattens between z~1 and z~3 to decrease at z greater than 3. Galaxies with different SEDs, masses and sSFRs evolve in very different ways and this large and deep statistical sample is the first one allowing us to separately study the different evolutionary behaviours of the individual IR populations contributing to the IR luminosity density. Galaxies occupying the well established SFR/stellar mass main sequence (MS) are found to dominate both the total IR LF and luminosity density at all redshifts, with the contribution from off MS sources (0.6 dex above MS) being nearly constant (~20% of the total IR luminosity density) and showing no significant signs of increase with increasing z over the whole 0.8<z<2.2 range. Sources with mass in the 10< log(M/Msun) <11 range are found to dominate the total IR LF, with more massive galaxies prevailing at the bright end of the high-z LF. A two-fold evolutionary scheme for IR galaxies is envisaged: on the one hand, a starburst-dominated phase in which the SMBH grows and is obscured by dust, is followed by an AGN dominated phase, then evolving toward a local elliptical. On the other hand, moderately starforming galaxies containing a low-luminosity AGN have various properties suggesting they are good candidates for systems in a transition phase preceding the formation of steady spiral galaxies.

We present in this paper the general real- and redshift-space clustering properties of galaxies as measured in the first data release of the VIPERS survey. VIPERS is a large redshift survey designed to probe the distant Universe and its large-scale structure at 0.5 < z < 1.2. We describe in this analysis the global properties of the sample and discuss the survey completeness and associated corrections. This sample allows us to measure the galaxy clustering with an unprecedented accuracy at these redshifts. From the redshift-space distortions observed in the galaxy clustering pattern we provide a first measurement of the growth rate of structure at z = 0.8: f\sigma_8 = 0.47 +/- 0.08. This is completely consistent with the predictions of standard cosmological models based on Einstein gravity, although this measurement alone does not discriminate between different gravity models.

Galaxies grow primarily via accretion-driven star formation in discs and merger-driven growth of bulges. These processes are implicit in semi-analytical models of galaxy formation, with bulge growth in particular relating directly to the hierarchical build-up of halos and their galaxies. In this paper, we consider several implementations of two semi-analytical models. Focusing on implementations in which bulges are formed during mergers only, we examine the fractions of elliptical galaxies and both passive and star-forming disk galaxies as functions of stellar and halo mass, for central and satellite systems. This is compared to an observational cross-matched SDSS+RC3 z ~ 0 sample of galaxies with accurate visual morphological classifications and M_{stellar} > 10^10.5 M_{sol}. The models qualitatively reproduce the observed increase of elliptical fraction with stellar mass, and with halo mass for central galaxies, supporting the idea that observed ellipticals form during major mergers. However, the overall elliptical fraction produced by the models is much too high compared with the z ~ 0 data. Since the "passive" -- i.e. non-star-forming -- fractions are approximately reproduced, and since the fraction which are star-forming disc galaxies is also reproduced, the problem is that the models overproduce ellipticals at the expense of passive S0 and spiral galaxies. Bulge-growth implementations (tuned to reproduce simulations) which allow the survival of residual discs in major mergers still destroy too much of the disc. Increasing the lifetime of satellites, or allowing significant disc regrowth around merger remnants, merely increases the fraction of star-forming disc galaxies. Instead, it seems necessary to reduce the mass ratios of merging galaxies, so that most mergers produce modest bulge growth in disc-galaxy remnants instead of ellipticals. [Abridged]

We present a new parallel implementation of the PINpointing Orbit Crossing-Collapsed HIerarchical Objects (PINOCCHIO) algorithm, a quick tool, based on Lagrangian Perturbation Theory, for the hierarchical build-up of Dark Matter halos in cosmological volumes. To assess its ability to predict halo correlations on large scales, we compare its results with those of an N-body simulation of a 3 Gpc/h box sampled with 2048^3 particles taken from the MICE suite, matching the same seeds for the initial conditions. Thanks to the FFTW libraries and to the relatively simple design, the code shows very good scaling properties. The CPU time required by PINOCCHIO is a tiny fraction (~1/2000) of that required by the MICE simulation. Varying some of PINOCCHIO numerical parameters allows one to produce a universal mass function that lies in the range allowed by published fits, although it underestimates the MICE mass function of FoF halos in the high mass tail. We compare the matter-halo and the halo-halo power spectra with those of the MICE simulation and find that these 2-point statistics are well recovered on large scales. In particular, when catalogs are matched in number density, agreement within ten per cent is achieved for the halo power spectrum. At scales k>0.1 h/Mpc, the inaccuracy of the Zel'dovich approximation in locating halo positions causes an underestimate of the power spectrum that can be modeled as a Gaussian factor with a damping scale of d=3 Mpc/h at z=0, decreasing at higher redshift. Finally, a remarkable match is obtained for the reduced halo bispectrum, showing a good description of nonlinear halo bias. Our results demonstrate the potential of PINOCCHIO as an accurate and flexible tool for generating large ensembles of mock galaxy surveys, with interesting applications for the analysis of large galaxy redshift surveys.

(Abridged) We explore new observationally-constrained sub-resolution models of galactic outflows and investigate their impact on the circumgalactic medium (CGM) over redshifts z = 2 - 4. We perform cosmological hydrodynamic simulations, including star formation, chemical enrichment, and four cases of SNe-driven outflows: no wind (NW), an energy-driven constant velocity wind (CW), a radially varying wind (RVWa) where the outflow velocity has a positive correlation with galactocentric distance (r), and a RVW with additional dependence on halo mass (RVWb). Overall, we find that the outflows expel metal-enriched gas away from galaxies, significantly quench star formation, and enrich the CGM. At z = 2, the radial profiles of gas properties around galaxy centers are most sensitive to the choice of the wind model for halo masses (10^9 - 10^11) M_sun. We infer that the RVWb model is similar to the NW case, except that it substantially enriches the CGM: the carbon metallicity (Z_C) is 10 times higher in RVWb than in NW at r > R_200; and the warm gas of 10^4 - 10^5 K and delta < 100 is enriched to 50 times higher than in NW. We also find that the impact of models CW and RVWa are similar, with the following differences. RVWa causes a greater suppression of star formation rate at z < 5, and has a higher fraction of low-density (delta < 10), warm-hot (10^4 - 10^6 K) gas than in CW. Outflows in CW produce a higher and earlier enrichment of some IGM phases than RVWa. By visual inspection, we note that the RVWa model shows more pronounced bipolar outflows and galactic disks. We present fitting formulae for [Z_C-delta] and [Z_C-r], also for the abundance of CIV as a function of r. We predict observational diagnostics to distinguish between different outflow scenarios: Z_C of the CGM gas at r = (30 - 300) kpc/h comoving, and CIV fraction of the inner gas at r < (4 - 5) kpc/h comoving.

We study the impact of star formation and stellar feedback prescriptions on galaxy properties predicted by means of "stripped-down" versions of independently developed semi-analytic models (SAMs). These include cooling, star formation, feedback from supernovae (SNe) and simplified prescriptions for galaxy merging, but no chemical evolution, disc instabilities or AGN feedback. We run these versions on identical samples of dark matter (DM) haloes extracted from high-resolution N-body simulations in order to perform both statistical analysis and object-by-object comparisons. We compare our results with previous work based on stripped-down versions of the same SAMs including only gas cooling, and show that all feedback models provide coherent modifications in the distribution of baryons between the various gas phases. In particular, we find that the predicted hot gas fractions are considerably increased by up to a factor of three, while the corresponding cold gas fractions are correspondingly decreased, and a significant amount of mass is ejected from the DM halo. Nonetheless, we also find relevant differences in the predicted properties of model galaxies among the three SAMs: these deviations are more relevant at mass scales comparable to that of our own Galaxy, and are reduced at larger masses, confirming the varying impact of stellar feedback at different mass scales. We also check the effect of enhanced star formation events (i.e. starbursts modes), defined in connection with galaxy mergers. We find that, in general, these episodes have a limited impact in the overall star formation histories of model galaxies, even in massive DM halos where merger-driven star formation has often been considered very important.

We address the fundamental question of matching the rest-frame K-band luminosity function (LF) of galaxies over the Hubble time using semi-analytic models, after modification of the stellar population modelling. We include the Maraston evolutionary synthesis models, that feature a higher contribution by the Thermally Pulsating - Asymptotic Giant Branch (TP-AGB) stellar phase, into three different semi-analytic models, namely the De Lucia and Blaizot version of the Munich model, MORGANA and the Menci model. We leave all other input physics and parameters unchanged. We find that the modification of the stellar population emission can solve the mismatch between models and the observed rest-frame K-band luminosity from the brightest galaxies derived from UKIDSS data at high redshift. For all explored semi-analytic models this holds at the redshifts - between 2 and 3 - where the discrepancy was recently pointed out. The reason for the success is that at these cosmic epochs the model galaxies have the right age (~1 Gyr) to contain a well-developed TP-AGB phase which makes them redder without the need of changing their mass or age. At the same time, the known overestimation of the faint end is enhanced in the K-band when including the TP-AGB contribution. At lower redshifts (z<2) some of the explored models deviate from the data. This is due to too short merging timescales and inefficient 'radio-mode' AGN feedback. Our results show that a strong evolution in mass predicted by hierarchical models is compatible with no evolution on the bright-end of the K-band LF from z=3 to the local universe. This means that, at high redshifts and contrary to what is commonly accepted, K-band emission is not necessarily a good tracer of galaxy mass.

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.

We present results from high--resolution cosmological hydrodynamical simulations of a Milky--Way-sized halo, aimed at studying the effect of feedback on the nature of gas accretion. Simulations include a model of inter-stellar medium and star formation, in which SN explosions provide effective thermal feedback. We distinguish between gas accretion onto the halo, which occurs when gas particles cross the halo virial radius, and gas accretion onto the central galaxy, which takes place when gas particles cross the inner one-tenth of the virial radius. Gas particles can be accreted through three different channels, depending on the maximum temperature value, $T_{\rm max}$, reached during the particles' past evolution: a cold channel for $T_{\rm max}<2.5 \times 10^5$ K, a hot one for $T>10^6$K, and a warm one for intermediate values of $T_{\rm max}$. We find that the warm channel is at least as important as the cold one for gas accretion onto the central galaxy. This result is at variance with previous findings that the cold mode dominates gas accretion at high redshift. We ascribe this difference to the different supernova feedback scheme implemented in our simulations. While results presented so far in the literature are based on uneffective SN thermal feedback schemes and/or the presence of a kinetic feedback, our simulations include only effective thermal feedback. We argue that observational detections of a warm accretion mode in the high--redshift circum-galactic medium would provide useful constraints on the nature of the feedback that regulates star formation in galaxies.

We study several versions of the Schmidt-Kennicutt (SK) relation obtained for isolated spiral galaxies in TreeSPH simulations run with the GADGET3 code including the novel MUlti-Phase Particle Integrator (MUPPI) algorithm for star formation and stellar feedback. [...] The standard SK relation between surface densities of cold (neutral+molecular) gas and star formation rate of simulated galaxies shows a steepening at low gas surface densities, starting from a knee whose position depends on disc gas fraction: for more gas-rich discs the steepening takes place at higher surface densities. Because gas fraction and metallicity are typically related, this environmental dependence mimics the predictions of models where the formation of H2 is modulated by metallicity. The cold gas surface density at which HI and molecular gas surface densities equate can range from ~10 up to 34 Msun/pc^2. As expected, the SK relation obtained using molecular gas shows much smaller variations among simulations. We find that disc pressure is not well represented by the classical external pressure of a disc in vertical hydrostatic equilibrium. Instead is well fit by the expression P_fit = Sigma_cold sigma_cold kappa / 6, where the three quantities on the right-hand side are cold gas surface density, vertical velocity dispersion and epicyclic frequency. When the "dynamical" SK relation, i.e. the relation that uses gas surface density divided by orbital time, is considered, we find that all of our simulations stay on the same relation. We interpret this as a manifestation of the equilibrium between energy injection and dissipation in stationary galaxy discs, when energetic feedback is effective and pressure is represented by the expression given above. These findings further support the idea that a realistic model of the structure of galaxy discs should take into account energy injection by SNe. [Abridged]

Dark matter halo merger trees are now routinely extracted from cosmological simulations of structure formation. These trees are frequently used as inputs to semi-analytic models of galaxy formation to provide the backbone within which galaxy formation takes place. By necessity, these merger trees are constructed from a finite set of discrete "snapshots" of the N-body simulation and so have a limited temporal resolution. To date, there has been little consideration of how this temporal resolution affects the properties of galaxies formed within these trees. In particular, the question of how many snapshots are needed to achieve convergence in galaxy properties has not be answered. Therefore, we study the convergence in the stellar and total baryonic masses of galaxies, distribution of merger times, stellar mass functions and star formation rates in the Galacticus model of galaxy formation as a function of the number of "snapshot" times used to represent dark matter halo merger trees. When utilizing snapshots between z=20 and z=0, we find that at least 128 snapshots are required to achieve convergence to within 5% for galaxy masses. This convergence is obtained for mean quantities averaged over large samples of galaxies - significant variance for individual galaxies remains even when using very large numbers of snapshots. We find only weak dependence of the rate of convergence on the distribution of snapshots in time - snapshots spaced uniformly in the expansion factor, uniformly in the logarithm of expansion factor or uniformly in the logarithm of critical overdensity for collapse work equally well in almost all cases. We provide input parameters to Galacticus which allow this type of convergence study to be tuned to other simulations and to be carried out for other galaxy properties.

We study the physical properties, formation histories, and environment of galaxies without a significant "classical" spheroidal component, as predicted by semi-analytical models of galaxy formation and evolution. This work is complementary to the analysis presented in De Lucia et al., (2011), where we focus on the relative contribution of various physical mechanisms responsible for bulge assembly in a Lambda-CDM cosmology. We find that the fraction of bulgeless galaxies is a strong decreasing function of stellar mass: they represent a negligible fraction of the galaxy population with M* > 10^12 Msun, but dominate at M* < 10^10 Msun. We find a clear dichotomy in this galaxy population, between central galaxies of low-mass dark matter haloes, and satellite galaxies in massive groups/clusters. We show that bulgeless galaxies are relatively young systems, that assemble most of their mass at low-redshift, but they can also host very old stellar populations. Since galaxy-galaxy mergers are assumed to lead to the formation of a spheroidal component, in our models these galaxies form preferentially in low-mass haloes that host a small number of satellites galaxies. We show that the adopted modelling for galaxy mergers represents a key ingredient in determining the actual number of bulgeless galaxies. Our results show that these galaxies are not a rare population in theoretical models: at z~0, galaxies with no classical bulge (but often including galaxies with the equivalent of pseudo-bulges) account for up to 14% of the galaxies with 10^11 < M*/Msun < 10^12.

We analyze predictions from two independently developed galaxy formation models to study the mechanisms, environments, and characteristic times of bulge formation in a LambdaCDM cosmogony. For each model, we test different prescriptions for bulge formation in order to quantify the relative importance of different channels. Our results show that the strong correlation between galaxy and halo mass for central galaxies, and the richer merger history of more massive systems naturally give rise to a strong correlation between galaxy mass and morphology, and between halo mass and morphological type of central galaxies. Large fractions of the bulge mass are acquired through major and minor mergers, but disc instability plays an important role, particularly for intermediate mass galaxies. We find that the modelling of disc instability events, as well as of the galaxy merger times, can affect significantly the timing of bulge formation, and the relative importance of different channels. Bulge dominated galaxies acquire their morphology through major mergers, but this can be modified by cooling of gas from the surrounding hot halo. We find that disc regrowth is a non negligible component of the evolution of bulge dominated galaxies, particularly for low to intermediate masses, and at high redshifts.

On May 27th 2010, the Italian astronomical community learned with concern that the National Institute for Astrophysics (INAF) was going to be suppressed, and that its employees were going to be transferred to the National Research Council (CNR). It was not clear if this applied to all employees (i.e. also to researchers hired on short-term contracts), and how this was going to happen in practice. In this letter, we give a brief historical overview of INAF and present a short chronicle of the few eventful days that followed. Starting from this example, we then comment on the current situation and prospects of astronomical research in Italy.

We use stripped-down versions of three semi-analytic galaxy formation models to study the influence of different assumptions about gas cooling and galaxy mergers. By running the three models on identical sets of merger trees extracted from high-resolution cosmological N-body simulations, we are able to perform both statistical analyses and halo-by-halo comparisons. Our study demonstrates that there is a good statistical agreement between the three models used here, when operating on the same merger trees, reflecting a general agreement in the underlying framework for semi-analytic models. We also show, however, that various assumptions that are commonly adopted to treat gas cooling and galaxy mergers can lead to significantly different results, at least in some regimes. In particular, we find that the different models adopted for gas cooling lead to similar results for mass scales comparable to that of our own Galaxy. Significant differences, however, arise at larger mass scales. These are largely (but not entirely) due to different treatments of the `rapid cooling' regime, and different assumptions about the hot gas distribution. At this mass regime, the predicted cooling rates can differ up to about one order of magnitude, with important implications on the relative weight that these models give to AGN feedback in order to counter-act excessive gas condensation in relatively massive haloes at low redshift. Different assumptions in the modelling of galaxy mergers can also result in significant differences in the timings of mergers, with important consequences for the formation and evolution of massive galaxies.

We present evidence that 'bona fide' disks and starburst systems occupy distinct regions in the gas mass versus star formation (SF) rate plane, both for the integrated quantities and for the respective surface densities. This result is based on CO observations of galaxy populations at low and high redshifts, and on the current consensus for the CO luminosity to gas mass conversion factors. The data suggest the existence of two different star formation regimes: a long-lasting mode for disks and a more rapid mode for starbursts, the latter probably occurring during major mergers or in dense nuclear SF regions. Both modes are observable over a large range of SF rates. The detection of CO emission from distant near-IR selected galaxies reveals such bimodal behavior for the first time, as they allow us to probe gas in disk galaxies with much higher SF rates than are seen locally. The different regimes can potentially be interpreted as the effect of a top-heavy IMF in starbursts. However, we favor a different physical origin related to the fraction of molecular gas in dense clouds. The IR luminosity to gas mass ratio (i.e., the SF efficiency) appears to be inversely proportional to the dynamical (rotation) timescale. Only when accounting for the dynamical timescale, a universal SF law is obtained, suggesting a direct link between global galaxy properties and the local SF rate.

We present a new multi-phase sub-resolution model for star formation and feedback in SPH numerical simulations of galaxy formation. Our model, called MUPPI (MUlti-Phase Particle Integrator), describes each gas particle as a multi-phase system, with cold and hot gas phases, coexisting in pressure equilibrium, and a stellar component. Cooling of the hot tenuous gas phase feeds the cold gas phase. Stars are formed out of molecular gas with a given efficiency, which scales with the dynamical time of the cold phase. Our prescription for star formation is not based on imposing the Schmidt-Kennicutt relation, which is instead naturally produced by MUPPI. Energy from supernova explosions is deposited partly into the hot phase of the gas particles, and partly to that of neighboring particles. Mass and energy flows among the different phases of each particle are described by a set of ordinary differential equations which we explicitly integrate for each gas particle, instead of relying on equilibrium solutions. This system of equations also includes the response of the multi-phase structure to energy changes associated to the thermodynamics of the gas. We apply our model to two isolated disk galaxy simulations and two spherical cooling flows. MUPPI is able to reproduce the Schmidt-Kennicutt relation for disc galaxies. It also reproduces the basic properties of the inter-stellar medium in disc galaxies, the surface densities of cold and molecular gas, of stars and of star formation rate, the vertical velocity dispersion of cold clouds and the flows connected to the galactic fountains. Quite remarkably, MUPPI also provides efficient stellar feedback without the need to include a scheme of kinetic energy feedback. [abridged]

[abridged] The properties of galaxies with the reddest observed R-K colors (Extremely Red Objects, EROs), including their apparent division into passive and obscured active objects with roughly similar number densities, are a known challenge for models of galaxy formation. We produce mock catalogues generated by interfacing the predictions of the semi-analytical MORGANA model for the evolution of galaxies in a Lambda-CDM cosmology with the spectro-photometric + radiative transfer code GRASIL and Infrared (IR) template library to show that the model correctly reproduces number counts, redshift distributions and active fractions of R-K>5 sources. We test the robustness of our results against different dust attenuations and, most importantly, against the inclusion of TP-AGB stars in Simple Stellar Populations used to generate galaxy spectra, and find that the inclusion of TP-AGBs has a relevant effect, in that it allows to increase by a large factor the number of very red active objects at all color cuts. We find that though the most passive and the most obscured active galaxies have a higher probability of being selected as EROs, many EROs have intermediate properties and the population does not show bimodality in specific star formation rate (SSFR). We predict that deep observations in the Far-IR, from 100 to 500 micron, are the most efficient way to constrain the SSFR of these objects; we give predictions for future Herschel observations. Finally, we test whether a simple evolutionary sequence for the formation of z=0 massive galaxies, going through a sub-mm-bright phase and then a ERO phase, are typical in this galaxy formation model. We find that this sequence holds for ~25 per cent of z=0 massive galaxies, while the model typically shows a more complex connection between sub-mm, ERO and massive galaxies. [abridged]

We study the luminosity function and the radial distribution of satellite galaxies within Milky Way sized haloes as predicted in Cold Dark Matter based models of galaxy formation, making use of numerical N-body techniques as well as three different semi-analytic model (SAMs) galaxy formation codes. We extract merger trees from very high-resolution dissipationless simulations of four Galaxy-sized DM haloes, and use these as common input for the semi-analytic models. We present a detailed comparison of our predictions with the observational data recently obtained on the Milky Way satellite luminosity function (LF). We find that semi-analytic models with rather standard astrophysical ingredients are able to reproduce the observed luminosity function over six orders of magnitude in luminosity, down to magnitudes as faint as M_V=-2. We also perform a comparison with the actual observed number of satellites as a function of luminosity, by applying the selection criteria of the SDSS survey to our simulations instead of correcting the observations for incompleteness. Using this approach we again find good agreement for both the luminosity and radial distributions of MW satellites. We investigate which physical processes in our models are responsible for shaping the predicted satellite LF, and find that tidal destruction, suppression of gas infall by a photo-ionizing background, and supernova feedback all make important contributions. We conclude that the number and luminosity of Milky Way satellites can be naturally accounted for within the (Lambda) Cold Dark Matter paradigm, and this should no longer be considered a problem.

It has recently been shown that galaxy formation models within the LambdaCDM cosmology predict that, compared to the observed population, small galaxies (with stellar masses < 10^{11} M_sun) form too early, are too passive since z ~ 3 and host too old stellar populations at z=0. We then expect an overproduction of small galaxies at z > 4 that should be visible as an excess of faint Lyman-break galaxies. To check whether this excess is present, we use the MORGANA galaxy formation model and GRASIL spectro-photometric + radiative transfer code to generate mock catalogues of deep fields observed with HST-ACS. We add observational noise and the effect of Lyman-alpha emission, and perform color-color selections to identify Lyman-break galaxies. The resulting mock candidates have plausible properties that closely resemble those of observed galaxies. We are able to reproduce the evolution of the bright tail of the luminosity function of Lyman-break galaxies (with a possible underestimate of the number of the brightest i-dropouts), but uncertainties and degeneracies in dust absorption parameters do not allow to give strong constraints to the model. Besides, our model shows a clear excess with respect to observations of faint Lyman-break galaxies, especially of z_{850} ~ 27 V-dropouts at z ~ 5. We quantify the properties of these "excess" galaxies and discuss the implications: these galaxies are hosted in dark matter halos with circular velocities in excess of 100 km s^{-1}, and their suppression may require a deep re-thinking of stellar feedback processes taking place in galaxy formation.

[abridged] It has been widely claimed that several lines of observational evidence point towards a "downsizing" (DS) of the process of galaxy formation over cosmic time. This behavior is sometimes termed "anti-hierarchical", and contrasted with the "bottom-up" assembly of the dark matter structures in Cold Dark Matter models. In this paper we address three different kinds of observational evidence that have been described as DS: the stellar mass assembly, star formation rate and the ages of the stellar populations in local galaxies. We compare a broad compilation of available data-sets with the predictions of three different semi-analytic models of galaxy formation within the Lambda-CDM framework. In the data, we see only weak evidence at best of DS in stellar mass and in star formation rate. We find that, when observational errors on stellar mass and SFR are taken into account, the models acceptably reproduce the evolution of massive galaxies, over the entire redshift range that we consider. However, lower mass galaxies are formed too early in the models and are too passive at late times. Thus, the models do not correctly reproduce the DS trend in stellar mass or the archaeological DS, while they qualitatively reproduce the mass-dependent evolution of the SFR. We demonstrate that these discrepancies are not solely due to a poor treatment of satellite galaxies but are mainly connected to the excessively efficient formation of central galaxies in high-redshift haloes with circular velocities ~100-200 km/s. [abridged]

We investigate the correlation of star formation quenching with internal galaxy properties and large scale environment (halo mass) in empirical data and theoretical models. We make use of the halo-based Group Catalog of Yang and collaborators, which is based on the Sloan Digital Sky Survey. Data from the Galaxy Evolution Explorer (GALEX) are also used to extract the recent star formation rate. In order to investigate the environmental effects, we examine the properties of "central" and "satellite" galaxies separately. For central galaxies, we are unable to conclude whether star formation quenching is primarily connected with halo mass or stellar mass, because these two quantities are themselves strongly correlated. For satellite galaxies, a nearly equally strong dependence on halo mass and stellar mass is seen. We make the same comparison for five different semi-analytic models based on three independently developed codes. We find that the models with AGN feedback reproduce reasonably well the dependence of the fraction of central red and passive galaxies on halo mass and stellar mass. However, for satellite galaxies, the same models badly overproduce the fraction of red/passive galaxies and do not reproduce the empirical trends with stellar mass or halo mass. This {\em satellite overquenching problem} is caused by the too-rapid stripping of the satellites' hot gas halos, which leads to rapid strangulation of star formation.

The treatment of dust attenuation is crucial in order to compare the predictions of galaxy formation models with multiwavelength observations. Most past studies have either used simple analytic prescriptions or else full radiative transfer (RT) calculations. Here, we couple star formation histories and morphologies predicted by the semi-analytic galaxy formation model MORGANA with RT calculations from the spectrophotometric and dust code GRASIL to create a library of galaxy SEDs from the UV/optical through the far Infrared, and compare the predictions of the RT calculations with analytic prescriptions. We consider a low and high redshift sample, as well as an additional library constructed with empirical, non-cosmological star formation histories and simple (pure bulge or disc) morphologies. Based on these libraries, we derive fitting formulae for the effective dust optical depth as a function of galaxy physical properties such as metallicity, gas mass, and radius. We show that such fitting formulae can predict the V-band optical depth with a scatter smaller than 0.4 dex for both the low and high redshift samples, but that there is a large galaxy-to-galaxy scatter in the shapes of attenuation curves, probably due to geometrical variations, which our simple recipe does not capture well. However, our new recipe provides a better approximation to the GRASIL results at optical wavelength than standard analytic prescriptions from the literature, particularly at high redshift.

We revisit the proximity effect produced by QSOs at redshifts 2.1-3.3 applying the FLO approach (Saitta et al. 2008) to a sample of ~6300 Ly-alpha lines fitted in 21 high resolution, high signal-to-noise spectra. This new technique allows to recover the hydrogen density field from the HI column densities of the lines in the Ly-alpha forest, on the basis of simple assumptions on the physical state of the gas. To minimize the systematic uncertainties that could affect the density recovering in the QSO vicinity, we carefully determined the redshifts of the QSOs in our sample and modelled in detail their spectra to compute the corresponding ionising fluxes. The mean density field obtained from the observed spectra shows a significant over-density in the region within 4 proper Mpc from the QSO position, confirming that QSOs are hosted in high density peaks. The absolute value of rho/<rho> for the peak is uncertain by a factor of ~3, depending on the assumed QSO spectral slope and the minimum HI column density detectable in the spectra. We do not confirm the presence of a significant over-density extending to separations of ~15 proper Mpc from the QSO, claimed in previous works at redshifts <z>=2.5 and 3.8. Our best guess for the UV background ionisation rate based on the IGM mean density recovered by FLO is Gamma_UVB ~ 10^{-12} s^{-1}. However, values of Gamma_UVB ~ 3x10^{-12} s^{-1} could be viable if an inverted temperature-density relation with index alpha=-0.5 is adopted.

[Abridged] We present FLO (From Lines to Over-densities), a new technique to reconstruct the hydrogen density field for the Lya forest lines observed in high resolution QSO spectra. The method is based on the hypothesis that the Lya lines arise in the low to intermediate density intergalactic gas and that the Jeans length is the typical size of the Lya absorbers. The reliability of FLO is tested against mock spectra obtained from cosmological simulations. The recovering algorithm gives satisfactory results in the range from the mean density to over-densities of ~30 and reproduces correctly the correlation function of the density field and the 1D power spectrum on scales between ~20 and 60 comoving Mpc. A sample of Lya forests from 22 high resolution QSO spectra is analysed, covering the redshift range 1.7<z<3.5. For each line of sight, we fit Voigt profiles to the lines of the Lya forest, providing the largest, homogeneous sample of fitted Lya lines ever studied. The line number density evolution with redshift follows a power-law relation: dn/dz=(166 +/- 4) [(1+z)/3.5]^{(2.8 +/- 0.2)} (1 sigma errors). The two-point correlation function of lines shows a signal up to separations of ~2 comoving Mpc; weak lines (log N(HI)<13.8) also show a significant clustering but on smaller scales (r<1.5 comoving Mpc). We estimate with FLO the hydrogen density field toward the 22 observed lines of sight. The redshift distribution of the average densities computed for each QSO is consistent with the cosmic mean hydrogen density in the analysed redshift range. The two-point correlation function and the 1D power spectrum of the delta field are estimated. The correlation function shows clustering signal up to ~4 comoving Mpc.

In order to study the process of cooling in dark-matter (DM) halos and assess how well simple models can represent it, we run a set of radiative SPH hydrodynamical simulations of isolated halos, with gas sitting initially in hydrostatic equilibrium within Navarro-Frenk-White (NFW) potential wells. [...] After having assessed the numerical stability of the simulations, we compare the resulting evolution of the cooled mass with the predictions of the classical cooling model of White & Frenk and of the cooling model proposed in the MORGANA code of galaxy formation. We find that the classical model predicts fractions of cooled mass which, after about two central cooling times, are about one order of magnitude smaller than those found in simulations. Although this difference decreases with time, after 8 central cooling times, when simulations are stopped, the difference still amounts to a factor of 2-3. We ascribe this difference to the lack of validity of the assumption that a mass shell takes one cooling time, as computed on the initial conditions, to cool to very low temperature. [...] The MORGANA model [...] better agrees with the cooled mass fraction found in the simulations, especially at early times, when the density profile of the cooling gas is shallow. With the addition of the simple assumption that the increase of the radius of the cooling region is counteracted by a shrinking at the sound speed, the MORGANA model is also able to reproduce for all simulations the evolution of the cooled mass fraction to within 20-50 per cent, thereby providing a substantial improvement with respect to the classical model. Finally, we provide a very simple fitting function which accurately reproduces the cooling flow for the first ~10 central cooling times. [Abridged]

In order to gain insight into the physical mechanisms leading to the formation of stars and their assembly in galaxies, we compare the predictions of the MOdel for the Rise of GAlaxies aNd Active nuclei (MORGANA) to the properties of K- and 850 micron-selected galaxies (such as number counts, redshift distributions and luminosity functions) by combining MORGANA with the spectrophotometric model GRASIL. We find that it is possible to reproduce the K- and 850 micron-band datasets at the same time and with a standard Salpeter IMF, and ascribe this success to our improved modeling of cooling in DM halos. We then predict that massively star-forming discs are common at z~2 and dominate the star-formation rate, but most of them merge with other galaxies within ~100 Myr. Our preferred model produces an overabundance of bright galaxies at z<1; this overabundance might be connected to the build-up of the diffuse stellar component in galaxy clusters, as suggested by Monaco et al. (2006), but a naive implementation of the mechanism suggested in that paper does not produce a sufficient slow-down of the evolution of these objects. Moreover, our model over-predicts the number of 10^{10}-10^{11} M_sun galaxies at z~1; this is a common behavior of theoretical models as shown by Fontana et al. (2006). These findings show that, while the overall build-up of the stellar mass is correctly reproduced by galaxy formation models, the ``downsizing'' trend of galaxies is not fully reproduced yet. This hints to some missing feedback mechanism in order to reproduce at the same time the formation of both the massive and the small galaxies.

We measure the luminosity function of QSOs in the redshift range 3.5 < z < 5.2 for the absolute magnitude interval -21 < M_{145} < -28. Suitable criteria are defined to select faint QSOs in the GOODS fields, checking their effectiveness and completeness in detail. The confirmed sample of faint QSOs is compared with a brighter one derived from the SDSS. Using a Monte-Carlo technique we estimate the properties of the luminosity function. Our results show that models based on pure density evolution show better agreement with observation than models based on pure luminosity evolution, even if a different break magnitude with respect to z ~ 2.1 is required at 3.5 < z < 5.2. According to our modeling a faint-end slope steeper than low-redshift observations is required to reproduce the data, moreover models with a steep bright-end slope score a higher probability than models with a bright-end flattening. Determining the faint-end of the luminosity function at these redshifts provides important constraints on models of the joint evolution of galaxies and AGNs.

We present some results of the new MORGANA model for the rise of galaxies and active nuclei, and show that the improved physical motivation of the description of star formation and feedback allows to get hints on the physical processes at play. We propose that the high level of turbulence in star-forming bulges is at the base of the observed downsizing of AGNs. In this framework it is also possible to reproduce the recently obtained evidence that most low-redshift accretion is powered by relatively massive, slowly accreting black holes. Besides, we notice that many galaxy formation models (including MORGANA) fail to reproduce a basic observable, namely the number density of 10^{11} Msun galaxies at z~1, as traced by the GOODS-MUSIC sample. This points to a possibly missing ingredient in the modeling of stellar feedback.

We present the MOdel for the Rise of GAlaxies aNd Active nuclei (MORGANA). Starting from the merger trees of dark matter halos and a model for the evolution of substructure within the halos, the complex physics of baryons is modeled with a set of state-of-the-art models that describe the mass, metal and energy flows between the various components and phases of a galaxy. The processes of shock-heating and cooling, star formation, feedback, galactic winds and super-winds, accretion onto BHs and AGN feedback are described by new models. In particular, the evolution of the halo gas explicitly follows the thermal and kinetic energies of the hot and cold phases, while star formation and feedback follow the results of the multi-phase model by Monaco (2004a). The increased level of sophistication allows to move from a phenomenological description of gas physics, based on simple scalings with the depth of the DM halo potential, toward a fully physically motivated one. The comparison of the predictions of MORGANA with a basic set of galactic data reveals from the one hand an overall rough agreement, and from the other hand highlights a number of well- or less-known problems: (i) producing the cutoff of the luminosity function requires to force the quenching of the late cooling flows by AGN feedback, (ii) the normalization of the Tully-Fisher relation of local spirals cannot be recovered unless the dark matter halos are assumed to have a very low concentration, (iii) the mass function of HI gas is not easily fitted at small masses, unless a similarly low concentration is assumed, (iv) there is an excess of small elliptical galaxies at z=0. These discrepancies, more than the points of agreement with data, give important clues on the missing ingredients of galaxy formation. (ABRIDGED)

Aims: In this work the luminosity function of QSOs is measured in the redshift range 3.5<z<5.2 for the absolute magnitude interval -21<M_{145}<-28. The determination of the faint end of the luminosity function at these redshifts provides important constraints on models of joint evolution of galaxies and AGNs. Methods: We have defined suitable criteria to select faint QSOs in the GOODS fields, checking in detail their effectiveness and completeness. Spectroscopic follow-up of the resulting QSO candidates has been carried out. The confirmed sample of faint QSOs is compared with a brighter one derived from the SDSS. We have used a Monte-Carlo technique to estimate the properties of the luminosity function, checking various parameterizations for its shape and evolution. Results: Models based on Pure Density Evolution show better agreement with observation than models based on Pure Luminosity Evolution. However a different break magnitude with respect to z~2.1 is required at 3.5<z<5.2. Models with a steeper faint end score a higher probability. We do not find any evidence for a flattening of the bright end at redshift z>3.5. Conclusions: The estimated space density evolution of QSOs indicates a suppression of the formation and/or feeding of Supermassive Black Holes at these redshifts. The QSO contribution to the UV background is insufficient to ionize the IGM at 3.5<z<5.2.

The high end of the stellar mass function of galaxies is observed to have little evolution since z~1. This represents a stringent constraint for merger--based models, aimed at explaining the evolution of the most massive galaxies in the concordance LambdaCDM cosmology. In this Letter we show that it is possible to remove the tension between the above observations and model predictions by allowing a fraction of stars to be scattered to the Diffuse Stellar Component (DSC) of galaxy clusters at each galaxy merger, as recently suggested by the analysis of N-body hydrodynamical simulations. To this purpose, we use the MORGANA model of galaxy formation in a minimal version, in which gas cooling and star formation are switched off after z=1. In this way, any predicted evolution of the galaxy stellar mass function is purely driven by mergers. We show that, even in this extreme case, the predicted degree of evolution of the high end of the stellar mass function is larger than that suggested by data. Assuming instead that a significant fraction, ~30 per cent, of stars are scattered in the DSC at each merger event, leads to a significant suppression of the predicted evolution, in better agreement with observational constraints, while providing a total amount of DSC in clusters which is consistent with recent observational determinations.

In order to constrain the physical processes that regulate and downsize the AGN population, the predictions of the MOdel for the Rise of GAlaxies aNd Active nuclei (MORGANA) are compared to luminosity functions (LFs) of AGNs in the optical, soft X-ray and hard X-ray bands, to the local BH-bulge mass relation, and to the observed X-ray number counts and background. We also give predictions on the accretion rate of AGNs in units of the Eddington rate and on the BH--bulge relation expected at high redshift. We find that it is possible to reproduce the downsizing of AGNs within the hierarchical LambdaCDM cosmogony, and that the most likely responsible for this downsizing is the stellar kinetic feedback that arises in star-forming bulges as a consequence of the high level of turbulence and leads to a massive removal of cold gas in small elliptical galaxies. At the same time, to obtain good fits to the number of bright quasars we need to require that quasar-triggered galactic winds self-limit the accretion onto BHs. In all cases, the predicted BH--bulge relation steepens considerably with respect to the observed one at bulge masses <10^{11} Msun; this problem is related to a known excess in the predicted number of small bulges, common to most similar models, so that the reproduction of the correct number of faint AGNs is done at the cost of underestimating their BH masses. This highlights an insufficient downsizing of elliptical galaxies, and hints for another feedback mechanism able to act on the compact discs that form and soon merge at high redshift. The results of this paper reinforce the need for direct investigations of the feedback mechanisms in active galaxies, that will be possible with the next generation of astronomical telescopes from sub-mm to X-rays.

The goal of this work is to measure the evolution of the Galaxy Stellar Mass Function and of the resulting Stellar Mass Density up to redshift ~4, in order to study the assembly of massive galaxies in the high redshift Universe. We have used the GOODS-MUSIC catalog, containing ~3000 Ks-selected galaxies with multi-wavelength coverage extending from the U band to the Spitzer 8 micron band, of which 27% have spectroscopic redshifts and the remaining fraction have accurate photometric redshifts. On this sample we have applied a standard fitting procedure to measure stellar masses. We compute the Galaxy Stellar Mass Function and the resulting Stellar Mass Density up to redshift ~4, taking into proper account the biases and incompleteness effects. Within the well known trend of global decline of the Stellar Mass Density with redshift, we show that the decline of the more massive galaxies may be described by an exponential timescale of ~6 Gyrs up to z~1.5, and proceeds much faster thereafter, with an exponential timescale of ~0.6 Gyrs. We also show that there is some evidence for a differential evolution of the Galaxy Stellar Mass Function, with low mass galaxies evolving faster than more massive ones up to z~1-1.5 and that the Galaxy Stellar Mass Function remains remarkably flat (i.e. with a slope close to the local one) up to z~1-1.3. The observed behaviour of the Galaxy Stellar Mass Function is consistent with a scenario where about 50% of present-day massive galaxies formed at a vigorous rate in the epoch between redshift 4 and 1.5, followed by a milder evolution until the present-day epoch.

We discuss the role of feedback from AGNs on the formation of spheroidal galaxies. The energy released by an accreting Black Hole (BH) may be injected into the ISM through blast waves arising directly from the central engine, radiation pressure or radiative heating. A scenario is described in which radiative heating perturbs the methabolism of a star-forming spheroid, leading to a critical stage where SNe form a cold expanding shell, pushed out of the galaxy by radiation pressure from the AGN. This mechanism can regulate the BH--bulge relation to the observed value. However, this relation may be simply imprinted by the mechanism responsible for the nearly complete loss of angular momentum of the gas that accretes onto the BH. Using a novel model of galaxy formation that includes AGNs, we show that models without self-regulation have problems in reproducing the correct slope of the AGN luminosity function, while models with winds give a much better fit; however, all these models are almost indistingishable as far as their predicted BH--bulge relation is concerned. Finally, we show that the downsizing of the faint AGNs is most likely due to kinetic feedback in star-forming bulges.

We investigate the effect of peculiar velocities on the redshift space distribution of z>~2 galaxies, and we focus in particular on Ly-alpha emitters. We generate catalogues of dark matter (DM) halos and identify emitters with halos of the same co-moving space density (M(Ly-alpha emitters) ~ 3x10^11 M_sun). We decompose the peculiar velocity field of halos into streaming, gradient and random components, and compute and analyse these as a function of scale. Streaming velocities are determined by fluctuations on very large scales, strongly affected by sample variance, but have a modest impact on the interpretation of observations. Gradient velocities are the most important as they distort structures in redshift space, changing the thickness and orientation of sheets and filaments. Random velocities are typically below or of the same order as the typical observational uncertainty on the redshift. We discuss the importance of these effects for the interpretation of data on the large-scale structure as traced by Ly-alpha emitters (or similar kinds of astrophysical high-redshift objects), focusing on the induced errors in the viewing angles of filaments. We compare our predictions of velocity patterns for Ly-alpha emitters to observations and find that redshift clumping of Ly-alpha emitters, as reported for instance in the fields of high-redshift radio galaxies, does not allow to infer whether an observed field is sampling an early galaxy overdensity.

The shining of quasars is a likely trigger of massive galatic winds, able to remove most ISM from a star-forming spheroid. However, the mechanism responsible for the deposition of energy into the ISM is still unclear. Starting from a model for feedback in galaxy formation with a two-phase medium (Monaco 2004a), we propose that the perturbation induced by radiative heating from a quasar on the ISM triggers a critical change of feedback regime. In the feedback model, SNRs expanding in the hot and pressurized phase of a star-forming spheroid tipically become pressure-confined before the hot interior gas is able to cool. Whenever the evaporation flow due to radiative heating of the quasar is significant with respect to the star-formation rate, the SNRs reach the point where their interior gas cools before being confined, forming a thick cold shell. We show that in this conditions the shells percolate into a super-shell of cold gas that sweeps the whole galaxy. Radiation pressure then pushes the shell out of the galaxy. This self-limiting mechanism leads to a correlation between black hole and bulge masses. The insertion of a motivated wind trigger criterion in a hierarchical galaxy formation model shows however that winds are not necessary to obtain a good black hole--bulge correlation. In absence of winds, good results are obtained if the mechanism responsible for the creation of a reservoir of low-angular momentum gas (able to accrete onto the black hole) deposits mass at a rate proportional to the star-formation rate. Using a novel galaxy formation model, we show under which conditions black hole masses are self-limited by the wind mechanism described above, and outline the possible observational consequences of this self-limitation. [Abridged]

Type II supernovae (SNe), probably the most important contributors to stellar feedback in galaxy formation, explode within the very dense star-forming clouds, where the injected energy is most easily radiated away. The efficiency of type II SNe in injecting energy into the interstellar medium (ISM) and in re-heating a fraction of the original star-forming cloud is estimated with the aid of a two-phase model for the ISM of the cloud. We argue that when SNe explode the star-forming cloud has already been partially destroyed by ionizing light and winds from massive stars. SN remnants (SNRs) will first cause the collapse of most of the cloud gas into cold fragments, until the diffuse hot phase has a low enough density to make further radiative losses negligible. This is completed in ~3 Myr, with a modest energy loss of ~5 per cent of the total budget. We compute that a fraction ranging from 5 to 30 per cent of the cloud is reheated to a high temperature (from 10^5 to 10^7 K); these numbers are very uncertain, due to the very complicated nature of the problem, and to the uncertain role of thermal evaporation. Small star-forming clouds, less massive than ~10^4 Msun, will be destroyed by a single SN. In all cases, a high fraction of the energy from type II SNe (~80 per cent for large clouds, smaller but still significant for small clouds) will be available for heating the ISM.

We present a new (semi-)analytic model for feedback in galaxy formation. The ISM is modeled as a two-phase medium in pressure equilibrium. The remnants of exploding type II SNe percolate into super-bubbles (SBs) that sweep the ISM, heating the hot phase (if the SB is adiabatic) or cooling it (in the snowplow stage, when the interior gas of the SB has cooled). The resulting feedback regimes occur in well-defined regions of the space defined by vertical scale-length and surface density of the structure. When SBs blow out in the adiabatic regime, the efficiency of SNe in heating the ISM is ~5 per cent, with \~80 per cent of the energy budget injected into the external halo, and the outcoming ISM is self-regulated to a state similar to that found in the Milky Way. Feedback is most efficient when SBs are pressure-confined in the adiabatic regime. In some significant regions of the parameter space confinement takes place in the snowplow stage; then the hot phase has a lower temperature and star formation is quicker. In some critical cases, the hot phase is strongly depleted and the cold phase percolates the whole volume, giving rise to a sudden burst of star formation. Strong galactic winds are predicted to happen only in critical cases. This model provides a starting point for constructing a realistic grid of feedback solutions to be used in galaxy formation codes. The predictive power of this model extends to many properties of the ISM, so that most parameters can be constrained by reproducing the main properties of the Milky Way. (Abridged)

The Great Observatories Origins Deep Survey provides significant constraints on the space density of less luminous QSOs at high redshift, which is particularly important to understand the interplay between the formation of galaxies and super-massive black holes and to measure the QSO contribution to the UV ionizing background. We present the results of a search for high-z QSOs, identified in the two GOODS fields on the basis of deep imaging in the optical (with HST) and X-ray (Chandra), and discuss the allowed space density of QSOs in the early universe.

We present a sample of 17 high-redshift (3.5<z<5.2) QSO candidates in the 320 sq.arcmin area of the Great Observatories Origins Deep Survey, selected in the magnitude range 22.45<z_{850}<25.25 using deep imaging with the Advanced Camera for Surveys onboard the Hubble Space Telescope and the Advanced CCD Imaging Spectrometer onboard the Chandra X-ray Observatory. On the basis of seven spectroscopic and ten photometric redshifts we estimate that the final sample will contain between two and four QSOs with 4<z<5.2. A dearth of high-redshift, moderate-luminosity (M_{145}=~-23) QSOs is observed with respect to predictions based on a) the extrapolation of the z~2.7 luminosity function (LF), according to a pure luminosity evolution calibrated by the results of the Sloan Digital Sky Survey; and b) a constant universal efficiency in the formation of super-massive black holes (SMBHs) in dark-matter halos. Evidence is gathered in favor of a density evolution of the LF at high redshift and of a suppression of the formation or feeding of SMBHs in low-mass halos.

We apply the ZTRACE algorithm to the optical NOG and infra-red PSCz galaxy catalogues to reconstruct the pattern of primordial fluctuations that have generated our local Universe. We check that the density fields traced by the two catalogues are well correlated, and consistent with a linear relation (either in $\delta$ or in log(1+$\delta$)) with relative bias (of NOG with respect to PSCz) b_rel = 1.1 $\pm$ 0.1. The relative bias relation is used to fill the optical zone of avoidance at |b| < 20$^\circ$ using the PSCz galaxy density field. We perform extensive testing on simulated galaxy catalogues to optimize the reconstruction. The quality of the reconstruction is predicted to be good at large scales, up to a limiting wavenumber k_lim $\simeq$ 0.4 h/Mpc beyond which all information is lost. We find that the improvement due to the denser sampling of the optical catalogue is compensated by the uncertainties connected to the larger zone of avoidance. The initial conditions reconstructed from the NOG catalogue are found (analogously to those from the PSCz) to be consistent with Gaussian paradigm. We use the reconstructions to produce sets of initial conditions ready to be used for constrained simulations of our local Universe.

We report the discovery of an unusual low--ionization broad absorption line quasar at z=1.776 which exhibits absorption lines from many excited states of Feii. This member of the rare class of 'FeLoBAL' quasars was serendipitously found in a mid-infrared (MIR) survey of distant clusters carried out with the ISOCAM camera. ISO J005645.1-273816 has a high MIR to UV luminosity ratio, suggesting a strong dust obscuration plus emission from very hot dust. This characteristic makes MIR surveys particularly efficient at detecting LoBAL quasars.

PINOCCHIO (PINpointing Orbit-Crossing Collapsed Hierarchical Objects) is a new algorithm for identifying dark matter halos in a given numerical realisation of the linear density field in a hierarchical universe (Monaco et al. 2001). It is shown that Lagrangian perturbation theory, and in particular its ellipsoidal truncation, is able to predict accurately the collapse, in the orbit-crossing sense, of generic mass elements. Some points that have undergone orbit crossing are assigned to the network of filaments and sheets that connects the halos; it is demonstrated that this network resembles closely that found in N-body simulations. The code generates a catalogue of dark matter halos with known mass, position, velocity, merging history and angular momentum. It is shown that the predictions of the code are very accurate when compared with the results of large N-body simulations that cover a range of cosmological models, box sizes and numerical resolutions. The mass function is recovered with an accuracy of better than 10 per cent in number density for halos with at least 30-50 particles. A similar accuracy is reached in the estimate of the correlation length r_0. The good agreement is still valid on the object-by-object level, with 70-100 per cent of the objects with more than 50 particles in the simulations also identified by our algorithm. For these objects the masses are recovered with an error of 20-40 per cent, and positions and velocities with a root mean square error of ~1-2 Mpc (0.5-2 grid lengths) and ~100 km/s, respectively. The recovery of the angular momentum of halos is considerably noisier and accuracy at the statistical level is achieved only by introducing free parameters. The algorithm requires negligible computer time as compared with performing a numerical $N$-body simulation.

We study the ability of PINOCCHIO (PINpointing Orbit-Crossing Collapsed HIerarchical Objects) to predict the merging histories of dark matter (DM) haloes, comparing the PINOCCHIO predictions with the results of two large N-body simulations run from the same set of initial conditions. We focus our attention on quantities most relevant to galaxy formation and large-scale structure studies. PINOCCHIO is able to predict the statistics of merger trees with a typical accuracy of 20 per cent. Its validity extends to higher-order moments of the distribution of progenitors. The agreement is valid also at the object-by-object level, with 70-90 per cent of the progenitors cleanly recognised when the parent halo is cleanly recognised itself. Predictions are presented also for quantities that are usually not reproduced by semi-analytic codes, such as the two-point correlation function of the progenitors of massive haloes and the distribution of initial orbital parameters of merging haloes. For the accuracy of the prediction and for the facility with which merger histories are produced, PINOCCHIO provides a means to generate catalogues of DM haloes which is extremely competitive to large-scale N-body simulations, making it a suitable tool for galaxy formation and large-scale structure studies.

We present a new algorithm (PINOCCHIO, PINpointing Orbit-Crossing Collapsed HIerarchical objects) to predict accurately the formation and evolution of individual dark matter haloes in a given realization of an initial linear density field. Compared with the halo population formed in a large (360^3 particles) collisionless simulation of a CDM universe, our method is able to predict to better than 10 per cent statistical quantities such as the mass function, two-point correlation function and progenitor mass function of the haloes. Masses of individual haloes are estimated accurately as well, with errors typically of order 30 per cent in the mass range well resolved by the numerical simulation. These results show that the hierarchical formation of dark matter haloes can be accurately predicted using local approximations to the dynamics when the correlations in the initial density field are properly taken into account. The approach allows one to automatically generate a large ensemble of accurate merging histories of haloes with complete knowledge of their spatial distribution. The construction of the full merger tree for a 256^3 realisation requires a few hours of CPU-time on a personal computer, orders of magnitude faster than the corresponding N-body simulation would take, and without needing any extensive post-processing. The technique can be efficiently used, for instance, for generating the input for galaxy formation modeling.

Direct and indirect observational evidence leads to the conclusion that high redshift QSOs did shine in the core of early type proto-galaxies during their main episode of star formation. Exploting this fact, we derive the rate of formation of this kind of stellar systems at high redshift by using the QSO Luminosity Function. The elemental proportions in elliptical galaxies, the descendents of the QSO hosts, suggest that the star formation was more rapid in more massive objects. We show that this is expected to occur in Dark Matter haloes, when the processes of cooling and heating is considered. This is also confirmed by comparing the observed sub-mm counts to those derived by coupling the formation rate and the star formation rate of the spheroidal galaxies with a detailed model for their SED evolution. In this scenario SCUBA galaxies and Lyman Break Galaxies are early type proto-galaxies forming the bulk of their stars before the onset of QSO activity.

The algorithm ZTRACE of Monaco & Efstathiou (1999) is applied to the IRAS PSCz catalogue to reconstruct the initial conditions of our local Universe with a resolution down to ~5 Mpc/h. The 1-point PDF of the reconstructed initial conditions is consistent with the assumptions that (i) IRAS galaxies trace mass on scales of ~5 Mpc/h, and (ii) the statistics of primordial density fluctuations is Gaussian. We use simulated PSCz catalogues, constructed from N-body simulations with Gaussian initial conditions, to show that local non-linear bias can cause the recovered initial PDF (assuming no bias) to be non-Gaussian. However, for plausible bias models, the distortions of the recovered PDF would be difficult to detect using the volume finely sampled by the PSCz catalogue. So, for Gaussian initial conditions, a range of bias models remain compatible with our PSCz reconstruction results.

We use rotation curve kinematics of 1000 spiral galaxies to investigate if these galaxies could host BH (Black Hole) remnants that once powered the QSO phenomenon. The rotation curves of both early and late-type spirals are used to place upper limits on the BH mass. We find that the central objects in spirals are 10-100 times less massive than those detected in ellipticals, suggesting that the contribution to the QSO phenomenon by the BH remnants hosted in these objects is negligible. The upper limits on the mass of the BHs resident in late type spirals are stringent: $M_{BH}\leq 10^{6}-10^{7} M_{\odot}$, indicating that only low luminosity activity possibly occurred in these objects.

We show that the mass function of black holes expected from the past quasar activity (both visible and obscured) is consistent with the number of dormant black holes found in the bulges of nearby galaxies. The joint formation of quasars and bulges is addressed by means of an analytical model for galaxy formation, based on the hierarchical clustering of cold dark matter halos. The model is able to reproduce the main statistical properties of both populations under the hypotheses that (i) star formation and quasar shining follow an anti-hierarchical order, and (ii) galaxy morphology and final black hole mass are determined by the same physical process.

Older and more recent pieces of observational evidence suggest a strong connection between QSOs and galaxies; in particular, the recently discovered correlation between black hole and galactic bulge masses suggests that QSO activity is directly connected to the formation of galactic bulges. The cosmological problem of QSO formation is analyzed in the framework of an analytical model for galaxy formation; for the first time a joint comparison with galaxy and QSO observables is performed. In this model it is assumed that the same physical variable which determines galaxy morphology is able to modulate the mass of the black hole responsible for QSO activity. Both halo spin and the occurence of a major merger are considered as candidates to this role. The predictions of the model are compared to available data for the type-dependent galaxy mass functions, the star-formation history of elliptical galaxies, the QSO luminosity function and its evolution (including the obscured objects contributing to the hard-X-ray background), the mass function of dormant black holes and the distribution of black-hole -- bulge mass ratios. A good agreement with observations is obtained if the halo spin modulates the efficiency of black-hole formation, and if the galactic halos at $z=0$ have shone in an inverted order with respect to the hierarchical one (i.e., stars and black holes in bigger galactic halos have formed before those in smaller ones). This inversion of hierarchical order for galaxy formation, which reconciles galaxy formation with QSO evolution, is consistent with many pieces of observational evidence.

The cosmological mass function problem is analyzed in full detail in the case of 1D gravity, with analytical, semi-analytical and numerical techniques. The extended Press & Schechter theory is improved by detailing the relation between smoothing radius and mass of the objects. This is done by introducing in the formalism the concept of a growth curve for the objects. The predictions of the extended Press & Schechter theory are compared to large N-body simulations of flat expanding 1D universes with scale-free power spectra of primordial perturbations. The collapsed objects in the simulations are located with a clump-finding algorithm designed to find regions that have undergone orbit crossing or that are in the multi-stream regime (these are different as an effect of the finite size of the multi-stream regions). It is found that the semi-analytical mass function theory, which has no free parameters, is able to recover the properties of collapsed objects both statistically and object by object. In particular, the predictions of regions in orbit crossing are optimized by the use of Gaussian filtering, while the use of sharp k-space filtering apparently allows to reproduce the larger multi-stream regions. The mass function theory does not reproduce well the clumps found with the standard friends-of-friends algorithm; however, the performance of this algorithm has not been thoroughly tested in the 1D cosmology. Our preliminary analyses of the 3D case confirms that the techniques developed in this paper are precious in understanding the cosmological mass function problem in 3D.

In order to map the galaxy density field on small scales in the local universe, we use the Nearby Optical Galaxy (NOG) sample, which is currently one of the largest, nearly complete, magnitude-limited ($B\leq$ 14 mag), all-sky sample of nearby optical galaxies ($\sim$ 6400 galaxies with cz< 5500 km/s). We have corrected the redshift-dependent distances of these galaxies for non-cosmological motions by means of peculiar velocity field models. Relying on group assignments and on total B magnitudes fully corrected for internal and Galactic extinctions, we determine the total and morphological-type specific luminosity functions for field and grouped galaxies using their locations in real distance space. The related determination of the selection function is meant to be an important step in recovering the galaxy density field on small scales from the NOG sample. Local galaxy density parameters will be used in statistical studies of environmental effects on galaxy properties.

We present and test a new method for the reconstruction of cosmological initial conditions from a full-sky galaxy catalogue. This method, called ZTRACE, is based on a self-consistent solution of the growing mode of gravitational instabilities according to the Zel'dovich approximation and higher order in Lagrangian perturbation theory. Given the evolved redshift-space density field, smoothed on some scale, ZTRACE finds via an iterative procedure, an approximation to the initial density field for any given set of cosmological parameters; real-space densities and peculiar velocities are also reconstructed. The method is tested by applying it to N-body simulations of an Einstein-de Sitter and an open cold dark matter universe. It is shown that errors in the estimate of the density contrast dominate the noise of the reconstruction. As a consequence, the reconstruction of real space density and peculiar velocity fields using non-linear algorithms is little improved over those based on linear theory. The use of a mass-preserving adaptive smoothing, equivalent to a smoothing in Lagrangian space, allows an unbiased (although noisy) reconstruction of initial conditions, as long as the (linearly extrapolated) density contrast does not exceed unity. The probability distribution function of the initial conditions is recovered to high precision, even for Gaussian smoothing scales of ~ 5 Mpc/h, except for the tail at delta >~ 1. This result is insensitive to the assumptions of the background cosmology.

We derive the mass function of the relic black holes and compared with that of the Massive Dark Objects in galaxies. Under the assumption that accretion onto massive BH's powers the Active Galactic Nuclei, the mass function of the BH responsibile for the past activity of QSO/AGN is computed. Our results support the scenario in which the QSO phase has exclusively occurred in every proto-elliptical.

In this paper we derive the galaxy luminosity function from the Nearby Optical Galaxy (NOG) sample, which is a nearly complete, magnitude-limited (B<14 mag), all-sky sample of nearby optical galaxies (~6400 galaxies with cz<5500 km/s). For this local sample, we use galaxy distance estimates based on different peculiar velocity models. Therefore, the derivation of the luminosity function is carried out using the locations of field and grouped galaxies in real distance space. The local field galaxy luminosity function in the B system is well described by a Schechter function. The exact values of the Schechter parameters slightly depend on the adopted peculiar velocity field models. The shape of the luminosity function of spiral galaxies does not differ significantly from that of E-S0 galaxies. On the other hand, the late-type spirals and irregulars have a very steeply rising luminosity function towards the faint end, whereas the ellipticals appreciably decrease in number towards low luminosities. The presence of galaxy systems in the NOG sample does not affect significantly the field galaxy luminosity function, since environmental effects on the total luminosity function appear to be marginal.

The Press & Schechter ``numerical recipe'' is briefly reviewed, together with the recently proposed dynamical mass function theory, in which the mass function is constructed by using the powerful Lagrangian perturbation theory. The dynamical mass function is found in good agreement with the recent N-body simulations of Governato et al. (1998), in the case of an Einstein-de Sitter Universe. The definition of collapse, the relation between mass and smoothing radius, and the definition of structure in 1D Universes are discussed. A detailed comparison of the dynamical mass function to simulations reveals that the orbit-crossed regions in the simulation are correctly reproduced, while the fragmentation of the collapsed medium into structures cannot be done in a univocal way. Finally, we try to answer the question: why the hell does the Press & Schechter work?

We use an all--sky, complete sample of nearby galaxies, extracted from the LEDA data base, to map the optical galaxy density field in the nearby universe. In order to determine this field, we correct the redshift--dependent distances by testing some peculiar velocity field models and we correct the galaxy number density for the incompletion of the galaxy sample at large distances through the derivation of the galaxy luminosity function. Local galaxy density parameters calculated for different smoothing scales are meant to be used in forthcoming statistical studies of environmental effects on galaxy properties.

In our recent paper (Salucci et al. 1998) we have investigated the mass distribution function of massive dark objects in galaxies, exploiting the available optical and radio observations. Under the assumption that massive black holes power active galactic nuclei, we have compared the mass functions of massive dark objects and black holes responsible for the observed activity. We have found that a scenario with a single short burst per active galactic nucleus is in a good agreement with the available data. Here we summarize and discuss the main points of our study.

By correcting the redshift--dependent distances for peculiar motions through a number of peculiar velocity field models, we recover the true distances of a wide, all-sky sample of nearby galaxies (~ 6400 galaxies with velocities cz<5500 km/s), which is complete up to the blue magnitude B=14 mag. Relying on catalogs of galaxy groups, we treat ~2700 objects as members of galaxy groups and the remaining objects as field galaxies. We model the peculiar velocity field using: i) a cluster dipole reconstruction scheme; ii) a multi--attractor model fitted to the Mark II and Mark III catalogs of galaxy peculiar velocities. According to Mark III data the Great Attractor has a smaller influence on local dynamics than previously believed, whereas the Perseus-Pisces and Shapley superclusters acquire a specific dynamical role. Remarkably, the Shapley structure, which is found to account for nearly half the peculiar motion of the Local Group, is placed by Mark III data closer to the zone of avoidance with respect to its optical position. Our multi--attractor model based on Mark III data favors a cosmological density parameter Omega ~ 0.5 (irrespective of a biasing factor of order unity). Differences among distance estimates are less pronounced in the ~ 2000 - 4000 km/s distance range than at larger or smaller distances. In the last regions these differences have a serious impact on the 3D maps of the galaxy distribution and on the local galaxy density --- on small scales.

This thesis aims to review the cosmological mass function problem, both from the theoretical and the observational point of view, and to present a new mass function theory, based on realistic approximations for the dynamics of gravitational collapse. Chapter 1 gives a general introduction on gravitational dynamics in cosmological models. Chapter 2 gives a complete review of the mass function theory. Chapters 3 and 4 present the ``dynamical'' mass function theory, based on truncated Lagrangian dynamics and on the excursion set approach. Chapter 5 reviews the observational state-of-the-art and the main applications of the mass function theories described before. Finally, Chapter 6 gives conclusions and future prospects.

The statistical tools needed to obtain a mass function from realistic collapse time estimates are presented. Collapse dynamics has been dealt with in paper I of this series by means of the powerful Lagrangian perturbation theory and the simple ellipsoidal collapse model. The basic quantity considered here is the inverse collapse time F; it is a non-linear functional of the initial potential, with a non-Gaussian distribution. In the case of sharp k-space smoothing, it is demonstrated that the fraction of collapsed mass can be determined by extending to the F process the diffusion formalism introduced by Bond et al. (1991). The problem is then reduced to a random walk with a moving absorbing barrier, and numerically solved; an accurate analytical fit, valid for small and moderate resolutions, is found. For Gaussian smoothing, the F trajectories are strongly correlated in resolution. In this case, an approximation proposed by Peacock & Heavens (1990) can be used to determine the mass functions. Gaussian smoothing is preferred, as it optimizes the performances of dynamical predictions and stabilizes the F trajectories. The relation between resolution and mass is treated at a heuristic level, and the consequences of this approximation are discussed. The resulting mass functions, compared to the classical Press & Schechter (1974) one, are shifted toward large masses (confirming the findings of Monaco 1995), and tend to give more intermediate-mass objects at the expense of small-mass objects. However, the small-mass part of the mass function, which depends on uncertain dynamics and is likely to be affected by uncertainties in the resolution--mass relation, is not considered a robust prediction of this theory.

A new theory for determining the mass function of cosmic structures is presented. It relies on a realistic treatment of collapse dynamics. Gravitational collapse is analyzed in the Lagrangian perturbative framework. Lagrangian perturbations provide an approximation of truncated type, i.e. small-scale structure is filtered out. The collapse time is suitably defined as the instant at which orbit crossing takes place. The convergence of the Lagrangian series in predicting the collapse time of a homogeneous ellipsoid is demonstrated; it is also shown that third-order calculations are necessary in predicting collapse. Then, the Lagrangian prediction, with a correction for quasi-spherical perturbations, can be used to determine the collapse time of a homogeneous ellipsoid in a fast and precise way. Furthermore, ellipsoidal collapse can be considered as a particular truncation of the Lagrangian series. Gaussian fields with scale-free power spectra are then considered. The Lagrangian series for the collapse time is found to converge when the collapse time is not large. In this case, ellipsoidal collapse gives a fast and accurate approximation of the collapse time; spherical collapse is found to poorly reproduce the collapse time, even in a statistical sense. Analytical fits of the distribution functions of the inverse collapse times, as predicted by the ellipsoid model and by third-order Lagrangian theory, are given. These will be necessary for a determination of the mass function, which will be given in paper II.

Relying on samples of disk galaxies for which a detailed photometric bulge/disk decomposition has been provided in the literature, we examine the dependence of the bulge-to-disk luminosity ratio (B/D) on the blue absolute luminosity and on the environmental density. In our statistical analysis of various B/D data sets we pay particular attention to disentangling the role played by the galaxy morphology--galaxy density relation. Besides, we focus our attention on nearby ($z<0.01$) galaxies, for which we can provide a three-dimensional characterization of the local galaxy density. We find that the observed tendency of galaxies to have greater B/D with increasing galaxy density simply reflects the average decline of B/D towards later morphological types together with the morphology--density relation. This relation tends to give rise also to a greater proportion of bright bulges in denser regions, because the decrease of B/D towards later types is mostly due to a dimming of the bulge rather than to a brightening of the disk. But when we remove the effect induced by the morphology--density relation, we detect no clear evidence of a dependence of B/D on galaxy density. Furthermore, B/D turns out to be substantially unrelated to the blue absolute magnitude of the galaxy. We briefly discuss to what extent our results (partially) disagree with previous claims.

Non-spherical dynamical approximations and models for the gravitational collapse are used to extend the well-known Press \& Schechter (PS) approach, in order to determine analytical expressions for the mass function of cosmic structures. The problem is rigorously set up by considering the intrinsic Lagrangian nature of the mass function. The Lagrangian equations of motion of a cold and irrotational fluid in single-stream regime show that the shear, which is non-locally determined by all the matter field, is the quantity which characterizes non-spherical perturbations. The Zel'dovich approximation, being a self-consistent first-order Lagrangian and local one, is used as a suitable guide to develop realistic estimates of the collapse time of a mass clump, starting from the local initial values of density and shear. Both Zel'dovich-based \an\ and models and the homogeneous ellipsoidal model predict that more large-mass objects are expected to form than the usual PS relation. In particular, the homogeneous ellipsoid model is consistent at large masses with a Press \& Schechter mass function with a lower value of the \dc\ parameter, in the range 1.4$\div$1.6. This gives a dynamical explanation of why lower \dc\ values have been found to fit the results of several N-body simulations. When more small-scale structure is present, highly non-linear dynamical effects can effectively slow down the collapse rate of a perturbation, increasing the effective value of \dc. This may have interesting consequences on the abundance of large-mass high-redshift objects.

Using an extensive sample of nearby galaxies (the Nearby Galaxies Catalog, by Tully), we investigate the environment of the galaxies hosting low-luminosity AGNs (Seyferts and LINERs). We define the local galaxy density, adopting a new correction for the incompleteness of the galaxy sample at large distances. We consider both a complete sample of bright and nearby AGNs, identified from the nuclear spectra obtained in available wide optical spectroscopic surveys, and a complete sample of nearby Seyferts. Basically, we compare the local galaxy density distributions of the AGNs with those of non-AGN samples, chosen in order to match the magnitude and morphological type distributions of the AGN samples. We find, only for the early-type spirals more luminous than $\sim M^*$, that both LINERs and Seyferts tend to reside in denser environments on all the scales tested, from tenths of Mpc to a few Mpc; moreover Seyferts show an enhanced small-scale density segregation with respect to LINERs. This gives support to the idea that AGNs can be stimulated by interactions. On larger scales, tens of Mpc, we find that the AGNs hosted in luminous early-type spirals show a tendency to stay near the center of the Local Supercluster. Finally we discuss the interpretations of our findings and their consequences for some possible scenarios of AGN formation and evolution and for the problem of how AGNs trace the large-scale structures.

We show that observational data for four Abell clusters of galaxies support the Gur\-za\-dyan-Kocharyan relation between the Hausdorff dimension and the dynamical properties of a galaxy system. The Hausdorff dimension is calculated using the two-point correlation function, while the dynamical parameters are estimated using available data and reasonable assumptions on the mass function of galaxies. This result can have essential consequences in the understanding of the dynamical mechanisms that determine the fractal distribution of galaxies.

We examine the 10\m\ emission of the central regions of 281 spiral galaxies, after having compiled all ground-based, small-aperture ($\sim$5") broad-band photometric observations at \l$\sim$10\m\ (N magnitudes) published in the literature. We evaluate the compactness of the $\sim$10\m\ emission of galaxy nuclei by comparing these small-beam measures with the large-beam IRAS 12\m\ fluxes. In the analysis of different subsets of objects, we apply survival analysis techniques in order to exploit the information contained in ``censored'' data (i.e., upper limits on the fluxes). Seyferts are found to contain the most powerful nuclear sources of mid-infrared emission, which in $\sim$1/3 of cases provide the bulk of the emission of the entire galaxy; thus, mid-infrared emission in the outer disc regions is not uncommon in Seyferts. The 10\m\ emission of Seyferts appears to be unrelated to their X-ray emission. HII region-like nuclei are stronger mid-infrared sources than normal nuclei and LINER nuclei (whose level of emission is not distinguishable from that of normal nuclei). Interacting objects have, on average, greater 10\m\ luminosities than non-interacting ones and exhibit more compact emission. Early-type spirals have stronger and more compact 10\m\ emission than late-type ones. Barred spirals are brighter at $\sim$10\m\ than unbarred systems, essentially because they more frequently contain HII region-like nuclei. The results of our detailed comparison between the behaviour of various categories of objects stress that the 10\m\ emission of spiral nuclei is closely linked to the (predominantly non-thermal synchrotron) radio

We have examined the effect of the environmental density on the arm classification of an extensive sample of spiral galaxies included in the Nearby Galaxy Catalog (Tully, 1988a). We have also explored the dependence of the arm class of a galaxy on other factors, such as its blue absolute magnitude and its disk-to-total mass ratio, inferred in the literature either from the gradient of a good galaxy rotation curve or from a photometric mass decomposition method. We have found that the arm class is strongly related to the absolute magnitude in the mid-type spirals (in the sense that grand design galaxies are, on average, more luminous than flocculent objects), whilst this relation is considerably weaker in the early and late types. In general the influence of the local density on the arm structure appears to be much weaker than that of the absolute magnitude. The local density acts essentially in strengthening the arm class--absolute magnitude relation for the mid types, whereas no environmental density effects are observed in the early and late types. Using the most recent estimates of the disk-to-total mass ratio, we do not confirm this ratio to be a significant factor which affects the arm class; nevertheless, owing to poor statistics and large uncertanties, the issue remains open. Neither a local density effect nor an unambiguous bar effect on the disk-to-total mass ratio is detectable; the latter finding may challenge some theoretical viewpoints on the formation of bar structures.

Using a variety of parameters of local galaxy density, we have examined the effects of the environmental density on the presence of bar structures in spiral galaxies of various morphological types. For an extensive sample of nearby galaxies, listed in the ``Nearby Galaxies Catalogue'' (Tully, 1988a), we have found that the spirals characterized by a high local density tend to be barred if they are early-type and early-type if they are barred (at the significance level of $\sim$3 sigma), confirming some earlier suggestion of low statistical significance. This fact, which is observed substantially in low-luminosity spirals, indicates that galaxy interactions can stimulate the formation of bars, primarily in early-type, low-luminosity spirals. This is in partial, qualitative agreement with the most recent relevant N-body simulations. On the other hand, no significant density segregation is observed between pure S-shaped (S(s)) spirals and spirals with inner rings (S(r)), which are often associated with bars. {\it Subject headings:} galaxies: general --- galaxies: structure --- galaxies: interactions --- galaxies: clustering