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.

It is a widespread opinion that hydrogen reionization is mainly driven by primeval star-forming galaxies, with a minor role of high-z active galactic nuclei. Recent observations, however, challenge this notion, indicating a number of issues related to a galaxy-driven reionization scenario. We provide here an updated assessment of the space density of relatively faint (M1450~-22.5) AGNs at zspec~5.5 in order to improve the estimate of the photo-ionization rate contribution from accreting super massive black holes. Exploiting deep UV rest-frame ground-based spectra collected at the Very Large Telescope on the CANDELS/GOODS-South field and deep Chandra X-ray images in the CANDELS/GOODS-North and EGS areas, we find two relatively bright (M1450~-22.5) AGNs at zspec~5.5. We derive an AGN space density of Phi=1.29x10^-6 cMpc^-3 at z~5.5 and M1450~-22.5 by simply dividing their observed number by the cosmological volume in the range 5.0<z<6.1. Our estimate does not consider corrections for incompleteness, therefore it represents a lower limit, although uncertainties due to cosmic variance can still be significant. This value supports a high space density of AGNs at z>5, in contrast with previous claims mostly based on standard color selection, possibly affected by significant incompleteness. Our estimate for the AGN photo-ionization rate at z~5.5 is in agreement with the observed values at similar redshifts, which are needed to keep the intergalactic medium highly ionized. Upcoming JWST and giant ground based telescopes observations will improve the study of high-z AGNs and their contribution to the reionization of the Universe.

We present the results of a new study investigating the relationship between observed Ly$\alpha$ equivalent width ($W_{\lambda}$(Ly$\alpha$)) and the metallicity of the ionizing stellar population ($Z_{\star}$) for a sample of $768$ star-forming galaxies at $3 \leq z \leq 5$ drawn from the VANDELS survey. Dividing our sample into quartiles of rest-frame $W_{\lambda}$(Ly$\alpha$) across the range $-58 \unicode{xC5} \lesssim$ $W_{\lambda}$(Ly$\alpha$) $\lesssim 110 \unicode{xC5}$ we determine $Z_{\star}$ from full spectral fitting of composite far-ultraviolet (FUV) spectra and find a clear anti-correlation between $W_{\lambda}$(Ly$\alpha$) and $Z_{\star}$. Our results indicate that $Z_{\star}$ decreases by a factor $\gtrsim 3$ between the lowest $W_{\lambda}$(Ly$\alpha$) quartile ($\langle$$W_{\lambda}$(Ly$\alpha$)$\rangle=-18\unicode{xC5}$) and the highest $W_{\lambda}$(Ly$\alpha$) quartile ($\langle$$W_{\lambda}$(Ly$\alpha$)$\rangle=24\unicode{xC5}$). Similarly, galaxies typically defined as Lyman Alpha Emitters (LAEs; $W_{\lambda}$(Ly$\alpha$) $>20\unicode{xC5}$) are, on average, metal poor with respect to the non-LAE galaxy population ($W_{\lambda}$(Ly$\alpha$) $\leq20\unicode{xC5}$) with $Z_{\star}$$_{\rm{non-LAE}}\gtrsim 2 \times$ $Z_{\star}$$_{\rm{LAE}}$. Finally, based on the best-fitting stellar models, we estimate that the increasing strength of the stellar ionizing spectrum towards lower $Z_{\star}$ is responsible for $\simeq 15-25\%$ of the observed variation in $W_{\lambda}$(Ly$\alpha$) across our sample, with the remaining contribution ($\simeq 75-85\%$) being due to a decrease in the HI/dust covering fractions in low $Z_{\star}$ galaxies.

Reproducing the observed quenched fraction of satellite galaxies has been a long-standing issue for galaxy formation models. We modify the treatment of environmental effects in our state-of-the-art GAlaxy Evolution and Assembly (GAEA) semi-analytic model to improve our modelling of satellite galaxies. Specifically, we implement gradual stripping of hot gas, ram-pressure stripping of cold gas, and an updated algorithm to account for angular momentum exchanges between the gaseous and stellar disc components of model galaxies. Our updated model is in good agreement with the observed quenched fractions for central and satellite galaxies, and their dependencies on stellar mass and halo mass. We also find consistency between model predictions and observational estimates of quenching time for satellite galaxies; HI, H$_2$ fractions of central galaxies; and deficiencies of HI, H$_2$, SFR of galaxies in cluster haloes. Analysing results from our updated model, we find that the dominant quenching mechanisms are hot-gas stripping for low-mass galaxies, and AGN feedback for massive galaxies. The Ram-pressure stripping of cold gas has a limited impact on the predicted quenched fractions, nonetheless, is needed to reproduce the observed HI deficiencies.

Aims. Our aim is to estimate the intergalactic medium transmission towards UV-selected star-forming galaxies at redshift 4 and above and study the effect of the dust attenuation on these measurements. Methods. The ultra-violet spectrum of high redshift galaxies is a combination of their intrinsic emission and the effect of the Inter-Galactic medium (IGM) absorption along their line of sight. Using data coming from the unprecedented deep spectroscopy from the VANDELS ESO public survey carried out with the VIMOS instrument we compute both the dust extinction and the mean transmission of the IGM as well as its scatter from a set of 281 galaxies at z>3.87. Because of a degeneracy between the dust content of the galaxy and the IGM, we first estimate the stellar dust extinction parameter E(B-V) and study the result as a function of the dust prescription. Using these measurements as constraint for the spectral fit we estimate the IGM transmission Tr(Lyalpha). Both photometric and spectroscopic SED fitting are done using the SPectroscopy And photometRy fiTting tool for Astronomical aNalysis (SPARTAN) that is able to fit the spectral continuum of the galaxies as well as photometric data. Results. Using the classical Calzetti's attenuation law we find that E(B-V) goes from 0.11 at z=3.99 to 0.08 at z=5.15. These results are in very good agreement with previous measurements from the literature. We estimate the IGM transmission and find that the transmission is decreasing with increasing redshift from Tr(Lyalpha)=0.53 at z=3.99 to 0.28 at z=5.15. We also find a large standard deviation around the average transmission that is more than 0.1 at every redshift. Our results are in very good agreement with both previous measurements from AGN studies and with theoretical models.

One of the major challenges in observational cosmology is related to the redshift evolution of the average hydrogen ionization in the Universe. In order to probe the ionization level of the IGM and the ionization capabilities of bright QSOs at z=4, we have selected a sample of 2508 QSOs from the SDSS survey (DR14) at 3.6<z<4.6 and -29.0<M1450<-26.0. Starting from UV/optical rest-frame spectra of the whole QSO sample, we estimate the escape fraction and free path individually for each of them. We calculate the Lyman Continuum (LyC) escape fraction as the flux ratio blue-ward (~900 A rest-frame) and red-ward (~930 A rest-frame) of the Lyman limit. We then obtain the probability distribution function (PDF) of the individual free paths of the QSOs in the sample and study its evolution in luminosity and redshift, comparing our results with the ones found in literature so far. We find a lower limit to the mean LyC escape fraction of 0.49, in agreement with the values obtained for both brighter and fainter sources at the same redshift. We show that the free paths of ionizing photons are characterized by a skewed distribution function peaked at low values, with an average of ~49-59 proper Mpc at z~4, after excluding possible associated absorbers. This value is larger than the one obtained at the same redshift by many authors in the literature using different techniques. Finally, we also find that the redshift evolution of this parameter results to be possibly milder than previously thought. Our new determination of the mean free path at z~4 implies that previous estimates of the HI photo-ionization rate Gamma_HI available in the literature should be corrected by a factor of 1.2-1.7. These results have important implications if extrapolated at the epoch of reionization.

In this work, we analyze the role of AGN feedback in quenching star formation for massive, central galaxies in the local Universe. In particular, we compare the prediction of two semi-analytic models (L-GALAXIES and SAGE) featuring different schemes for AGN feedback, with the SDSS DR7 taking advantage of a novel technique for identifying central galaxies in an observational dataset. This enables us to study the correlation between the model passive fractions, which is predicted to be suppressed by feedback from an AGN, and the observed passive fractions in an observationally motivated parameter space. While the passive fractions for observed central galaxies show a good correlation with stellar mass and bulge mass, passive fractions in L-GALAXIES correlate with the halo and black hole mass. For SAGE, the passive fraction correlate with the bulge mass as well. Among the two models, SAGE has a smaller scatter in the black hole - bulge mass (M_BH - M_Bulge) relation and a slope that agrees better with the most recent observations at z \sim 0. Despite the more realistic prescription of radio mode feedback in SAGE, there are still tensions left with the observed passive fractions and the distribution of quenched galaxies. These tensions may be due to the treatment of galaxies living in non-resolved substructures and the resulting higher merger rates that could bring cold gas which is available for star formation.

We investigate the prevalence of galactic-scale outflows in post-starburst (PSB) galaxies at high redshift ($1 < z < 1.4$), using the deep optical spectra available in the UKIDSS Ultra Deep Survey (UDS). We use a sample of $\sim40$ spectroscopically confirmed PSBs, recently identified in the UDS field, and perform a stacking analysis in order to analyse the structure of strong interstellar absorption features such as Mg ii ($\lambda2800$ Ang.). We find that for massive ($M_* > 10^{10}\rm\,M_{\odot}$) PSBs at $z > 1$, there is clear evidence for a strong blue-shifted component to the Mg ii absorption feature, indicative of high-velocity outflows ($v_{\rm out}\sim1150\pm160\rm\,km\,s^{-1}$) in the interstellar medium. We conclude that such outflows are typical in massive PSBs at this epoch, and potentially represent the residual signature of a feedback process that quenched these galaxies. Using full spectral fitting, we also obtain a typical stellar velocity dispersion $\sigma_*$ for these PSBs of $\sim200\rm\,km\,s^{-1}$, which confirms they are intrinsically massive in nature (dynamical mass $M_{\rm d}\sim10^{11}\rm\,M_{\odot}$). Given that these high-$z$ PSBs are also exceptionally compact ($r_{\rm e}\sim1$--$2\rm\,kpc$) and spheroidal (Sersic index $n\sim3$), we propose that the outflowing winds may have been launched during a recent compaction event (e.g. major merger or disc collapse) that triggered either a centralised starburst or active galactic nuclei (AGN) activity. Finally, we find no evidence for AGN signatures in the optical spectra of these PSBs, suggesting they were either quenched by stellar feedback from the starburst itself, or that if AGN feedback is responsible, the AGN episode that triggered quenching does not linger into the post-starburst phase.

We extend our previous work focused at $z\sim0$, studying the redshift evolution of galaxy dynamical properties using the state-of-the-art semi-analytic model GAEA: we show that the predicted size-mass relation for disky/star forming and quiescent galaxies is in good agreement with observational estimates, up to $z\sim2$. Bulge dominated galaxies have sizes that are offset low with respect to observational estimates, mainly due to our implementation of disk instability at high redshift. At large masses, both quiescent and bulge dominated galaxies have sizes smaller than observed. We interpret this as a consequence of our most massive galaxies having larger gas masses than observed, and therefore being more affected by dissipation. We argue that a proper treatment of quasar driven winds is needed to alleviate this problem. Our model compact galaxies have number densities in agreement with observational estimates and they form most of their stars in small and low angular momentum high-$z$ halos. GAEA predicts that a significant fraction of compact galaxies forming at high-$z$ is bound to merge with larger structures at lower redshifts: therefore they are not the progenitors of normal-size passive galaxies at $z=0$. Our model also predicts a stellar-halo size relation that is in good agreement with observational estimates. The ratio between stellar size and halo size is proportional to the halo spin and does not depend on stellar mass but for the most massive galaxies, where AGN feedback leads to a significant decrease of the retention factor (from about 80 per cent to 20 per cent).

We present the results of a study utilising ultra-deep, rest-frame UV, spectroscopy to quantify the relationship between stellar mass and stellar metallicity for 681 star-forming galaxies at $2.5<z<5.0$ ($\langle z \rangle = 3.5 \pm 0.6$) drawn from the VANDELS survey. Via a comparison with high-resolution stellar population models, we determine stellar metallicities for a set of composite spectra formed from subsamples selected by mass and redshift. Across the stellar mass range $8.5 < \mathrm{log}(\langle M_{\ast} \rangle/\rm{M}_{\odot}) < 10.2$ we find a strong correlation between stellar metallicity and stellar mass, with stellar metallicity monotonically increasing from $Z_{\ast}/\mathrm{Z}_{\odot} < 0.09$ at $\langle M_{\ast} \rangle = 3.2 \times 10^{8} \rm{M}_{\odot}$ to $Z_{\ast}/Z_{\odot} = 0.27$ at $\langle M_{\ast} \rangle = 1.7 \times 10^{10} \rm{M}_{\odot}$. In contrast, at a given stellar mass, we find no evidence for significant metallicity evolution across the redshift range of our sample. However, comparing our results to the $z=0$ stellar mass-metallicity relation, we find that the $\langle z \rangle = 3.5$ relation is consistent with being shifted to lower metallicities by $\simeq 0.6$ dex. Contrasting our derived stellar metallicities with estimates of gas-phase metallicities at similar redshifts, we find evidence for enhanced $\rm{O}/\rm{Fe}$ ratios of the order (O/Fe) $\gtrsim 1.8$ $\times$ (O/Fe)$_{\odot}$. Finally, by comparing our results to simulation predictions, we find that the $\langle z \rangle = 3.5$ stellar mass-metallicity relation is consistent with current predictions for how outflow strength scales with galaxy mass. This conclusion is supported by an analysis of analytic models, and suggests that the mass loading parameter ($\eta=\dot{M}_{\mathrm{outflow}}/M_{\ast}$) scales as $\eta \propto M_{\ast}^{\beta}$ with $\beta \simeq -0.4$.

We investigate the physical properties of a sample of 52 Ly{\alpha} emitting galaxies in the VANDELS survey, with particular focus on the role of kinematics and neutral hydrogen column density in the escape and spatial distribution of Ly{\alpha} photons. For these galaxies, we derive different physical properties and characterize the Ly{\alpha} emission in terms of kinematics, EW, FWHM and spatial extension and then estimate the velocity of the neutral outflowing gas. We reproduce some of the well known correlations between Ly{\alpha} EW and stellar mass, dust extinction and UV \b{eta} slope, in the sense that the emission line appears brighter in lower mass, less dusty and bluer galaxies. We do not find any correlation with the SED-derived star formation rate, while we find that galaxies with brighter Ly{\alpha} tend to be more compact both in UV and in Ly{\alpha}. Our data reveal a new interesting correlation between the Ly{\alpha} velocity and the offset of the inter-stellar absorption lines with respect to the systemic redshift, in the sense that galaxies with larger inter-stellar medium (ISM) out-flow velocities show smaller Ly{\alpha} velocity shifts. We interpret this relation in the context of the shell-model scenario, where the velocity of the ISM and the HI column density contribute together in determining the Ly{\alpha} kinematics. In support to our interpretation, we observe that galaxies with high HI column densities have much more extended Ly{\alpha} spatial profiles, a sign of increased scattering. However, we do not find any evidence that the HI column density is related to any other physical properties of the galaxies, although this might be due in part to the limited range of parameters that our sample spans.

We study the shape and evolution of the star formation main sequence in three independently developed semi-analytic models of galaxy formation. We focus, in particular, on the characterization of the model galaxies that are significantly above the main sequence, and that can be identified with galaxies classified as `starburst' in recent observational work. We find that, in all three models considered, star formation triggered by merger events (both minor and major) contribute to only a very small fraction of the cosmic density of star formation. While mergers are associated to bursts of star formation in all models, galaxies that experienced recent merger events are not necessarily significantly above the main sequence. On the other hand, `starburst galaxies' are not necessarily associated with merger episodes, especially at the low-mass end. Galaxies that experienced recent mergers can have relatively low levels of star formation when/if the merger is gas-poor, and galaxies with no recent merger can experience episodes of starbursts due to a combination of large amount of cold gas available from cooling/accretion events and/or small disk radii which increases the cold gas surface density.

A wealth of observations recently challenged the notion of a universal stellar initial mass function (IMF) by showing evidences in favour of a variability of this statistical indicator as a function of galaxy properties. I present predictions from the semi-analytic model GAEA (GAlaxy Evolution and Assembly), which features (a) a detailed treatment of chemical enrichment, (b) an improved stellar feedback scheme, and (c) implements theoretical prescriptions for IMF variations. Our variable IMF realizations predict intrinsic stellar masses and mass-to-light ratios larger than those estimated from synthetic photometry assuming a universal IMF. This provides a self-consistent interpretation for the observed mismatch between photometrically inferred stellar masses of local early-type galaxies and those derived by dynamical and spectroscopic studies. At higher redshifts, the assumption of a variable IMF has a deep impact on our ability to reconstruct the evolution of the galaxy stellar mass function and the star formation history of galaxies.

We use our state-of-the-art Galaxy Evolution and Assembly (GAEA) semi-analytic model to study how and on which time-scales star formation is suppressed in satellite galaxies. Our fiducial stellar feedback model, implementing strong stellar driven outflows, reproduces relatively well the variations of passive fractions as a function of galaxy stellar mass and halo mass measured in the local Universe, as well as the `quenching' time-scales inferred from the data. We show that the same level of agreement can be obtained by using an alternative stellar feedback scheme featuring lower ejection rates at high redshift, and modifying the treatment for hot gas stripping. This scheme over-predicts the number densities of low to intermediate mass galaxies. In addition, a good agreement with the observed passive fractions can be obtained only by assuming that cooling can continue on satellites, at the rate predicted considering halo properties at infall, even after their parent dark matter substructure is stripped below the resolution of the simulation. For our fiducial model, the better agreement with the observed passive fractions can be ascribed to: (i) a larger cold gas fraction of satellites at the time of accretion, and (ii) a lower rate of gas reheating by supernovae explosions and stellar winds with respect to previous versions of our model. Our results suggest that the abundance of passive galaxies with stellar mass larger than ~10^10 Msun is primarily determined by the self-regulation between star formation and stellar feedback, with environmental processes playing a more marginal role.

We use a state-of-the-art semi-analytic model to study the size and the specific angular momentum of galaxies. Our model includes a specific treatment for the angular momentum exchange between different galactic components. Disk scale radii are estimated from the angular momentum of the gaseous/stellar disk, while bulge sizes are estimated assuming energy conservation. The predicted size--mass and angular momentum--mass relations are in fair agreement with observational measurements in the local Universe, provided a treatment for gas dissipation during major mergers is included. Our treatment for disk instability leads to unrealistically small radii of bulges formed through this channel, and predicts an offset between the size--mass relations of central and satellite early-type galaxies, that is not observed. The model reproduces the observed dependence of the size--mass relation on morphology, and predicts a strong correlation between specific angular momentum and cold gas content. This correlation is a natural consequence of galaxy evolution: gas-rich galaxies reside in smaller halos, and form stars gradually until present day, while gas-poor ones reside in massive halos, that formed most of their stars at early epochs, when the angular momentum of their parent halos is low. The dynamical and structural properties of galaxies can be strongly affected by a different treatment for stellar feedback, as this would modify their star formation history. A higher angular momentum for gas accreted through rapid mode does not affect significantly the properties of massive galaxies today, but has a more important effect on low-mass galaxies at higher redshift.

We present a comparison of the observed evolving galaxy stellar mass functions with the predictions of eight semi-analytic models and one halo occupation distribution model. While most models are able to fit the data at low redshift, some of them struggle to simultaneously fit observations at high redshift. We separate the galaxies into 'passive' and 'star-forming' classes and find that several of the models produce too many low-mass star-forming galaxies at high redshift compared to observations, in some cases by nearly a factor of 10 in the redshift range $2.5 < z < 3.0$. We also find important differences in the implied mass of the dark matter haloes the galaxies inhabit, by comparing with halo masses inferred from observations. Galaxies at high redshift in the models are in lower mass haloes than suggested by observations, and the star formation efficiency in low-mass haloes is higher than observed. We conclude that many of the models require a physical prescription that acts to dissociate the growth of low-mass galaxies from the growth of their dark matter haloes at high redshift.

In this paper, we present a new derivation of the shape and evolution of the integrated galaxy-wide initial mass function (IGIMF), incorporating explicitly the effects of cosmic rays (CRs) as regulators of the chemical and thermal state of the gas in the dense cores of molecular clouds. We predict the shape of the IGIMF as a function of star formation rate (SFR) and CR density, and show that it can be significantly different with respect to local estimates. In particular, we focus on the physical conditions corresponding to IGIMF shapes that are simultaneously shallower at high-mass end and steeper at the low-mass end than a Kroupa IMF. These solutions can explain both the levels of $\alpha$-enrichment and the excess of low-mass stars as a function of stellar mass, observed for local spheroidal galaxies. As a preliminary test of our scenario, we use idealized star formation histories to estimate the mean IMF shape for galaxies of different $z=0$ stellar mass. We show that the fraction of low-mass stars as a function of galaxy stellar mass predicted by these mean IMFs agrees with the values derived from high-resolution spectroscopic surveys.

VANDELS is a uniquely-deep spectroscopic survey of high-redshift galaxies with the VIMOS spectrograph on ESO's Very Large Telescope (VLT). The survey has obtained ultra-deep optical (0.48 < lambda < 1.0 micron) spectroscopy of ~2100 galaxies within the redshift interval 1.0 < z < 7.0, over a total area of ~0.2 sq. degrees centred on the CANDELS UDS and CDFS fields. Based on accurate photometric redshift pre-selection, 85% of the galaxies targeted by VANDELS were selected to be at z>=3. Exploiting the red sensitivity of the refurbished VIMOS spectrograph, the fundamental aim of the survey is to provide the high signal-to-noise ratio spectra necessary to measure key physical properties such as stellar population ages, masses, metallicities and outflow velocities from detailed absorption-line studies. Using integration times calculated to produce an approximately constant signal-to-noise ratio (20 < t_int < 80 hours), the VANDELS survey targeted: a) bright star-forming galaxies at 2.4 < z < 5.5, b) massive quiescent galaxies at 1.0 < z < 2.5, c) fainter star-forming galaxies at 3.0 < z < 7.0 and d) X-ray/Spitzer-selected active galactic nuclei and Herschel-detected galaxies. By targeting two extragalactic survey fields with superb multi-wavelength imaging data, VANDELS will produce a unique legacy data set for exploring the physics underpinning high-redshift galaxy evolution. In this paper we provide an overview of the VANDELS survey designed to support the science exploitation of the first ESO public data release, focusing on the scientific motivation, survey design and target selection.

This paper describes the observations and the first data release (DR1) of the ESO public spectroscopic survey "VANDELS, a deep VIMOS survey of the CANDELS CDFS and UDS fields". VANDELS' main targets are star-forming galaxies at 2.4<z<5.5 and massive passive galaxies at 1<z<2.5. By adopting a strategy of ultra-long exposure times, from 20 to 80 hours per source, VANDELS is designed to be the deepest ever spectroscopic survey of the high-redshift Universe. Exploiting the red sensitivity of the VIMOS spectrograph, the survey has obtained ultra-deep spectra covering the wavelength 4800-10000 A with sufficient signal-to-noise to investigate the astrophysics of high-redshift galaxy evolution via detailed absorption line studies. The VANDELS-DR1 is the release of all spectra obtained during the first season of observations and includes data for galaxies for which the total (or half of the total) scheduled integration time was completed. The release contains 879 individual objects with a measured redshift and includes fully wavelength and flux-calibrated 1D spectra, the associated error spectra, sky spectra and wavelength-calibrated 2D spectra. We also provide a catalog with the essential galaxy parameters, including spectroscopic redshifts and redshift quality flags. In this paper we present the survey layout and observations, the data reduction and redshift measurement procedure and the general properties of the VANDELS-DR1 sample. We also discuss the spectroscopic redshift distribution, the accuracy of the photometric redshifts and we provide some examples of data products. All VANDELS-DR1 data are publicly available and can be retrieved from the ESO archive. Two further data releases are foreseen in the next 2 years with a final release scheduled for June 2020 which will include improved re-reduction of the entire spectroscopic data set. (abridged)

Finding the sources responsible for the hydrogen reionization is one of the most pressing issues in cosmology. Bright QSOs are known to ionize their surrounding neighborhood, but they are too few to ensure the required HI ionizing background. A significant contribution by faint AGNs, however, could solve the problem, as recently advocated on the basis of a relatively large space density of faint active nuclei at z>4. We have carried out an exploratory spectroscopic program to measure the HI ionizing emission of 16 faint AGNs spanning a broad U-I color interval, with I~21-23 and 3.6<z<4.2. These AGNs are three magnitudes fainter than the typical SDSS QSOs (M1450<~-26) which are known to ionize their surrounding IGM at z>~4. The LyC escape fraction has been detected with S/N ratio of ~10-120 and is between 44 and 100% for all the observed faint AGNs, with a mean value of 74% at 3.6<z<4.2 and -25.1<M1450<-23.3, in agreement with the value found in the literature for much brighter QSOs (M1450<~-26) at the same redshifts. The LyC escape fraction of our faint AGNs does not show any dependence on the absolute luminosities or on the observed U-I colors. Assuming that the LyC escape fraction remains close to ~75% down to M1450~-18, we find that the AGN population can provide between 16 and 73% (depending on the adopted luminosity function) of the whole ionizing UV background at z~4, measured through the Lyman forest. This contribution increases to 25-100% if other determinations of the ionizing UV background are adopted. Extrapolating these results to z~5-7, there are possible indications that bright QSOs and faint AGNs can provide a significant contribution to the reionization of the Universe, if their space density is high at M1450~-23.

Recent studies proposed that cosmic rays (CR) are a key ingredient in setting the conditions for star formation, thanks to their ability to alter the thermal and chemical state of dense gas in the UV-shielded cores of molecular clouds. In this paper, we explore their role as regulators of the stellar initial mass function (IMF) variations, using the semi-analytic model for GAlaxy Evolution and Assembly (GAEA). The new model confirms our previous results obtained using the integrated galaxy-wide IMF (IGIMF) theory: both variable IMF models reproduce the observed increase of $\alpha$-enhancement as a function of stellar mass and the measured $z=0$ excess of dynamical mass-to-light ratios with respect to photometric estimates assuming a universal IMF. We focus here on the mismatch between the photometrically-derived ($M^{\rm app}_{\star}$) and intrinsic ($M_{\star}$) stellar masses, by analysing in detail the evolution of model galaxies with different values of $M_{\star}/M^{\rm app}_{\star}$. We find that galaxies with small deviations (i.e. formally consistent with a universal IMF hypothesis) are characterized by more extended star formation histories and live in less massive haloes with respect to the bulk of the galaxy population. While the IGIMF theory does not change significantly the mean evolution of model galaxies with respect to the reference model, a CR-regulated IMF implies shorter star formation histories and higher peaks of star formation for objects more massive than $10^{10.5} M_\odot$. However, we also show that it is difficult to unveil this behaviour from observations, as the key physical quantities are typically derived assuming a universal IMF.

We present a comparison of nine galaxy formation models, eight semi-analytical and one halo occupation distribution model, run on the same underlying cold dark matter simulation (cosmological box of co-moving width 125$h^{-1}$ Mpc, with a dark-matter particle mass of $1.24\times 10^9 h^{-1}$ Msun) and the same merger trees. While their free parameters have been calibrated to the same observational data sets using two approaches, they nevertheless retain some 'memory' of any previous calibration that served as the starting point (especially for the manually-tuned models). For the first calibration, models reproduce the observed z = 0 galaxy stellar mass function (SMF) within 3-{\sigma}. The second calibration extended the observational data to include the z = 2 SMF alongside the z~0 star formation rate function, cold gas mass and the black hole-bulge mass relation. Encapsulating the observed evolution of the SMF from z = 2 to z = 0 is found to be very hard within the context of the physics currently included in the models. We finally use our calibrated models to study the evolution of the stellar-to-halo mass (SHM) ratio. For all models we find that the peak value of the SHM relation decreases with redshift. However, the trends seen for the evolution of the peak position as well as the mean scatter in the SHM relation are rather weak and strongly model dependent. Both the calibration data sets and model results are publicly available.

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 update our recently published model for GAlaxy Evolution and Assembly (GAEA), to include a self-consistent treatment of the partition of cold gas in atomic and molecular hydrogen. Our model provides significant improvements with respect to previous ones used for similar studies. In particular, GAEA (i) includes a sophisticated chemical enrichment scheme accounting for non-instantaneous recycling of gas, metals, and energy; (ii) reproduces the measured evolution of the galaxy stellar mass function; (iii) reasonably reproduces the observed correlation between galaxy stellar mass and gas metallicity at different redshifts. These are important prerequisites for models considering a metallicity dependent efficiency of molecular gas formation. We also update our model for disk sizes and show that model predictions are in nice agreement with observational estimates for the gas, stellar and star forming disks at different cosmic epochs. We analyse the influence of different star formation laws including empirical relations based on the hydro-static pressure of the disk, analytic models, and prescriptions derived from detailed hydro-dynamical simulations. We find that modifying the star formation law does not affect significantly the global properties of model galaxies, neither their distributions. The only quantity showing significant deviations in different models is the cosmic molecular-to-atomic hydrogen ratio, particularly at high redshift. Unfortunately, however, this quantity also depends strongly on the modelling adopted for additional physical processes. Useful constraints on the physical processes regulating star formation can be obtained focusing on low mass galaxies and/or at higher redshift. In this case, self-regulation has not yet washed out differences imprinted at early time.

We study galaxy mass assembly and cosmic star formation rate (SFR) at high-redshift (z$\gt$4), by comparing data from multiwavelength surveys with predictions from the GAlaxy Evolution and Assembly (GAEA) model. GAEA implements a stellar feedback scheme partially based on cosmological hydrodynamical simulations, that features strong stellar driven outflows and mass-dependent timescale for the re-accretion of ejected gas. In previous work, we have shown that this scheme is able to correctly reproduce the evolution of the galaxy stellar mass function (GSMF) up to $z\sim3$. We contrast model predictions with both rest-frame Ultra-Violet (UV) and optical luminosity functions (LF), which are mostly sensible to the SFR and stellar mass, respectively. We show that GAEA is able to reproduce the shape and redshift evolution of both sets of LFs. We study the impact of dust on the predicted LFs and we find that the required level of dust attenuation is in qualitative agreement with recent estimates based on the UV continuum slope. The consistency between data and model predictions holds for the redshift evolution of the physical quantities well beyond the redshift range considered for the calibration of the original model. In particular, we show that GAEA is able to recover the evolution of the GSMF up to z$\sim$7 and the cosmic SFR density up to z$\sim$10.

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.

Star forming galaxies have long been considered the dominant sources of the cosmic ultraviolet background radiation at early epochs. However, observing and characterizing the galaxy population with significant ionizing emission has proven to be challenging. In particular, the fraction of ionizing radiation that escapes the local environment to the intergalactic medium is poorly known. We investigate the relation between the escape fraction and galaxy luminosity. We combine the deep ultraviolet observations of Hubble Ultra Deep Field (UVUDF) with the deep Multi Unit Spectroscopic Explorer (MUSE) observations of the same field, collecting a sample of 165 faint star forming galaxies in the $3 < z < 4$ redshift range with deep rest-frame observations of the Lyman continuum. In our sample, we do not find any galaxy with significant emission of LyC radiation. We bin the galaxies in various redshift and brightness intervals and stack their images. From stacked images we estimate the relative escape fraction upper limits as a function of the luminosity. Thanks to the depth of the sample we measure meaningful 1$\sigma$ upper limits of $f_{esc,rel} < 0.07, 0.2$ and 0.6 at $L \sim L_{\rm z=3}^{*}, 0.5L_{\rm z=3}^{*}$ and $0.1L_{\rm z=3}^{*}$, respectively. We use our estimates and theoretical predictions from the literature to study a possible dependence of the escape fraction on galaxy luminosity by modelling the ionizing background with different prescriptions of $f_{\rm esc} (M_{\rm UV})$. We show that the understanding of the luminosity dependence hinges on the ability to constrain the escape fraction down to $M_{\rm UV} \sim -18$ mag in the future.

We take advantage of our recently published model for Galaxy Evolution and Assembly (GAEA) to study the origin of the observed correlation between [$\alpha$/Fe] and galaxy stellar mass. In particular, we analyse the role of radio mode AGN feedback, that recent work has identified as a crucial ingredient to reproduce observations. In GAEA, this process introduces the observed trend of star formation histories extending over shorter time-scales for more massive galaxies, but does not provide a sufficient condition to reproduce the observed $\alpha$ enhancements of massive galaxies. In the framework of our model, this is possible only assuming that any residual star formation is truncated for galaxies more massive than $10^{10.5}\,{\rm M}_{\odot}$. This results, however, in even shorter star formation time-scales for the most massive galaxies, that translate in total stellar metallicities significantly lower than observed. Our results demonstrate that (i) trends of [$\alpha$/Fe] ratios cannot be simply converted in relative time-scale indicators, and (ii) AGN feedback cannot explain alone the positive correlation between [$\alpha$/Fe] and galaxy mass/velocity dispersion. Reproducing simultaneously the mass-metallicity relation and the $\alpha$ enhancements observed poses a challenge for hierarchical models, unless more exotic solutions are adopted such as metal-rich winds or a variable IMF.

In this work, we study the basic statistical properties of HI-selected galaxies extracted from six different semi-analytic models, all run on the same cosmological N-body simulation. One model includes an explicit treatment for the partition of cold gas into atomic and molecular hydrogen. All models considered agree nicely with the measured HI mass function in the local Universe, with the measured scaling relations between HI and galaxy stellar mass, and with the predicted 2-point correlation function for HI rich galaxies. One exception is given by one model that predicts very little HI associated with galaxies in haloes above 10^12 Msun: we argue this is due to a too efficient radio-mode feedback for central galaxies, and to a combination of efficient stellar feedback and instantaneous stripping of hot gas for satellites. We demonstrate that treatment of satellite galaxies introduces large uncertainties at low HI masses. While models assuming non instantaneous stripping of hot gas tend to form satellite galaxies with HI masses slightly smaller than those of centrals with the same stellar mass, instantaneous gas stripping does not translate necessarily in lower HI masses. In fact, the adopted stellar feedback and star formation affect the satellites too. We analyze the relation between HI content and spin of simulated haloes: low spin haloes tend to host HI poor galaxies while high spin haloes are populated by galaxies in a wide range of HI mass. In our simulations, this is due to a correlation between the initial gas disk size and the halo spin.

In this work, we investigate the implications of the Integrated Galaxy-wide stellar Initial Mass Function (IGIMF) approach in the framework of the semi-analytic model GAEA (GAlaxy Evolution and Assembly), which features a detailed treatment of chemical enrichment and stellar feedback. The IGIMF provides an analytic description of the dependence of the stellar IMF shape on the rate of star formation in galaxies. We find that our model with a universal IMF predicts a rather flat [$\alpha$/Fe]-stellar mass relation. The model assuming the IGIMF, instead, is able to reproduce the observed increase of $\alpha$-enhancement with stellar mass, in agreement with previous studies. This is mainly due to the fact that massive galaxies are characterized by larger star formation rates at high-redshift, leading to stronger $\alpha$-enhancement with respect to low-mass galaxies. At the same time, the IGIMF hypothesis does not affect significantly the trend for shorter star formation timescales for more massive galaxies. We argue that in the IGIMF scenario the [$\alpha$/Fe] ratios are good tracers of the highest star formation events. The final stellar masses and mass-to-light-ratio of our model massive galaxies are larger than those estimated from the synthetic photometry assuming a universal IMF, providing a self-consistent interpretation of similar recent results, based on dynamical analysis of local early type galaxies.

In this work, we investigate the abundance and distribution of metals in the intergalactic medium (IGM) at $\langle z \rangle \simeq 2.8$ through the analysis of an ultra-high signal-to-noise ratio UVES spectrum of the quasar HE0940-1050. In the CIV forest, our deep spectrum is sensitive at $3\,\sigma$ to lines with column density down to $\log N_{\rm CIV} \simeq 11.4$ and in 60 per cent of the considered redshift range down to $\simeq11.1$. In our sample, all HI lines with $\log N_{\rm HI} \ge 14.8$ show an associated CIV absorption. In the range $14.0 \le \log N_{\rm HI} <14.8$, 43 per cent of HI lines has an associated CIV absorption. At $\log N_{\rm HI} < 14.0$, the detection rates drop to $<10$ per cent, possibly due to our sensitivity limits and not to an actual variation of the gas abundance properties. In the range $\log N_{\rm HI} \ge 14$, we observe a fraction of HI lines with detected CIV a factor of 2 larger than the fraction of HI lines lying in the circum-galactic medium (CGM) of relatively bright Lyman-break galaxies hosted by dark matter haloes with $\langle M\rangle \sim10^{12}$ M$_{\odot}$. The comparison of our results with the output of a grid of photoionization models and of two cosmological simulations implies that the volume filling factor of the IGM gas enriched to a metallicity $\log Z/Z_{\odot} \ge -3$ should be of the order of $\sim 10-13$ percent. In conclusion, our results favour a scenario in which metals are found also outside the CGM of bright star-forming galaxies, possibly due to pollution by lower mass objects and/or to an early enrichment by the first sources.

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.

One major problem of current theoretical models of galaxy formation is given by their inability to reproduce the apparently `anti-hierarchical' evolution of galaxy assembly: massive galaxies appear to be in place since $z\sim 3$, while a significant increase of the number densities of low mass galaxies is measured with decreasing redshift. In this work, we perform a systematic analysis of the influence of different stellar feedback schemes, carried out in the framework of GAEA, a new semi-analytic model of galaxy formation. It includes a self-consistent treatment for the timings of gas, metal and energy recycling, and for the chemical yields. We show this to be crucial to use observational measurements of the metallicity as independent and powerful constraints for the adopted feedback schemes. The observed trends can be reproduced in the framework of either a strong ejective or preventive feedback model. In the former case, the gas ejection rate must decrease significantly with cosmic time (as suggested by parametrizations of the cosmological `FIRE' simulations). Irrespective of the feedback scheme used, our successful models always imply that up to 60-70 per cent of the baryons reside in an `ejected' reservoir and are unavailable for cooling at high redshift. The same schemes predict physical properties of model galaxies (e.g. gas content, colour, age, and metallicity) that are in much better agreement with observational data than our fiducial model. The overall fraction of passive galaxies is found to be primarily determined by internal physical processes, with environment playing a secondary role.

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 study the effect of mergers on the morphology of galaxies by means of the simulated merger tree approach first proposed by Moster et al. This method combines N-body cosmological simulations and semi-analytic techniques to extract realistic initial conditions for galaxy mergers. These are then evolved using high resolution hydrodynamical simulations, which include dark matter, stars, cold gas in the disc and hot gas in the halo. We show that the satellite mass accretion is not as effective as previously thought, as there is substantial stellar stripping before the final merger. The fraction of stellar disc mass transferred to the bulge is quite low, even in the case of a major merger, mainly due to the dispersion of part of the stellar disc mass into the halo. We confirm the findings of Hopkins et al., that a gas rich disc is able to survive major mergers more efficiently. The enhanced star formation associated with the merger is not localised to the bulge of galaxy, but a substantial fraction takes place in the disc too. The inclusion of the hot gas reservoir in the galaxy model contributes to reducing the efficiency of bulge formation. Overall, our findings suggest that mergers are not as efficient as previously thought in transforming discs into bulges. This possibly alleviates some of the tensions between observations of bulgeless galaxies and the hierarchical scenario for structure formation.

The distribution of galaxy morphological types is a key test for models of galaxy formation and evolution, providing strong constraints on the relative contribution of different physical processes responsible for the growth of the spheroidal components. In this paper, we make use of a suite of semi-analytic models to study the efficiency of galaxy mergers in disrupting galaxy discs and building galaxy bulges. In particular, we compare standard prescriptions usually adopted in semi-analytic models, with new prescriptions proposed by Kannan et al., based on results from high-resolution hydrodynamical simulations, and we show that these new implementations reduce the efficiency of bulge formation through mergers. In addition, we compare our model results with a variety of observational measurements of the fraction of spheroid-dominated galaxies as a function of stellar and halo mass, showing that the present uncertainties in the data represent an important limitation to our understanding of spheroid formation. Our results indicate that the main tension between theoretical models and observations does not stem from the survival of purely disc structures (i.e. bulgeless galaxies), rather from the distribution of galaxies of different morphological types, as a function of their stellar mass.

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.

Among the possible alternatives to the standard cosmological model ($\Lambda$CDM), coupled Dark Energy models postulate that Dark Energy (DE), seen as a dynamical scalar field, may interact with Dark Matter (DM), giving rise to a "fifth-force", felt by DM particles only. In this paper, we study the impact of these cosmologies on the statistical properties of galaxy populations by combining high-resolution numerical simulations with semi-analytic models (SAM) of galaxy formation and evolution. New features have been implemented in the reference SAM in order to have it run self-consistently and calibrated on these cosmological simulations. They include an appropriate modification of the mass temperature relation and of the baryon fraction in DM haloes, due to the different virial scalings and to the gravitational bias, respectively. Our results show that the predictions of our coupled-DE SAM do not differ significantly from theoretical predictions obtained with standard SAMs applied to a reference $\Lambda$CDM simulation, implying that the statistical properties of galaxies provide only a weak probe for these alternative cosmological models. On the other hand, we show that both galaxy bias and the galaxy pairwise velocity distribution are sensitive to coupled DE models: this implies that these probes might be successfully applied to disentangle among quintessence, $f(R)$-Gravity and coupled DE models.

The constraints on neutrino masses led to the revision of their cosmological role, since the existence of a cosmological neutrino background is a clear prediction of the standard cosmological model. In this paper, we study the impact of such background on the spatial distribution of both Dark Matter (DM) and galaxies, by coupling $N$-body numerical simulations with semi-analytic models (SAMs) of galaxy formation. Cosmological simulations including massive neutrinos predict a slower evolution of DM perturbations with respect to the $\Lambda$CDM runs with the same initial conditions and a suppression on the matter power spectrum on small and intermediate scales, thus impacting on the predicted properties of galaxy populations. We explicitly show that most of these deviations are driven by the different $\sigma_8$ predicted for cosmologies including a massive neutrino background. We conclude that independent estimates of $\sigma_8$ are needed, in order to unambiguously characterise the effect of this background on the growth of structures. Galaxy properties alone are a weak tracer of deviations with respect to the $\Lambda$CDM run, but their combination with the overall matter distribution at all scales allows to disentangle between different cosmological models. Moreover, these deviations go on opposite direction with respect to competing models like modified gravity, thus weakening any detectable cosmological signal. Given the ubiquitous presence of a neutrino background, these effects have to be taken into account in future missions aimed at constraining the properties of the "Dark" components of the Universe.

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.

Deviations from a universal, MW-like, Stellar Initial Mass Function (IMF) have been reported for distant galaxies, although the physical reason behind the observed variations is still matter of ongoing debate. In this paper, we present an exploratory study to assess the impact of the proposed IMF evolution on the statistical galaxy properties, as predicted by the Semi-Analytical model of galaxy formation and evolution MORGANA. In particular, we test different dependencies for the IMF shape, as a function of both model galaxy properties (such as star formation rate, velocity dispersion or stellar mass) and environment, and compare the predicted stellar mass functions and star formation rate functions with reference runs at fixed IMF. In most cases, MORGANA predictions show deviations of the order of a few tenths of dex with respect to a run assuming an Universal Kroupa IMF. Among the proposed IMF variations, an increasing Top-Heavy IMF at increasing star formation rates has the largest impact on predicted galaxy properties, while most of the models assuming an increasing Bottom-Heavy IMF at higher masses/velocity dispersion lead to variations in galaxy properties that are of the same order as the uncertainty on the mass and star formation rate determination. By comparing the predicted galaxy stellar mass functions, we conclude that the study of the high-mass end can provide useful constraints to disentangle models assuming an increasing Top-Heavy IMF in high star forming or Bottom-Heavy IMF in massive systems.

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.

We study the observed cosmic ionizing background as a constraint on the nature of the sources responsible for the reionization of the Universe. In earlier work, we showed that extrapolations of the Ultra-Violet Luminosity Function (LF) of Lyman Break Galaxies (LBGs) at fixed Lyman continuum photon escape fraction are not able to reproduce the redshift evolution of this background. Here, we employ extrapolations of the high-z LFs to describe the contribution of LBGs to the ionizing photon rate, taking into account the smoothing of the baryonic perturbations, due to the background itself (i.e. the filtering mass), as well as a possible sharp increase of the escape fraction in dwarf galaxies. Under the hypothesis of a dominant contribution of LBGs to cosmic reionization, our results suggest that sources fainter than the current observational limits should be characterised by escape fractions of the order of ~0.1-0.3 (larger than the current estimates for bright galaxies) to account for a z>6 reionization and the measured evolution of cosmic ionizing background, at the same time. The contribution to the background from quasars turns out to be relevant at z<3. Overall, our results support the case for dedicated observations of faint galaxies in the rest-frame UV, in order to better determine their physical properties. Observed escape fractions outside our proposed range bear relevant consequences on the nature of the astrophysical sources responsible for cosmic reionization and/or its buildup process.

Modifications of the equations of general relativity at large distances offer one possibility to explain the observed properties of our Universe without invoking a cosmological constant. Numerous proposals for such modified gravity cosmologies exist, but often their consequences for structure formation in the non-linear sector are not yet accurately known. In this work, we employ high-resolution numerical simulations of f(R)-gravity models coupled with a semi-analytic model (SAM) for galaxy formation to obtain detailed predictions for the evolution of galaxy properties. The f(R)-gravity models imply the existence of a `fifth-force', which is however locally suppressed, preserving the successes of general relativity on solar system scales. We show that dark matter haloes in f(R)-gravity models are characterized by a modified virial scaling with respect to the LCDM scenario, reflecting a higher dark matter velocity dispersion at a given mass. This effect is taken into account in the SAM by an appropriate modification of the mass--temperature relation. We find that the statistical properties predicted for galaxies (such as the stellar mass function and the cosmic star formation rate) in f(R)-gravity show generally only very small differences relative to LCDM, smaller than the dispersion between the results of different SAM models, which can be viewed as a measure of their systematic uncertainty. We also demonstrate that galaxy bias is not able to disentangle between f(R)-gravity and the standard cosmological scenario. However, f(R)-gravity imprints modifications in the linear growth rate of cosmic structures at large scale, which can be recovered from the statistical properties of large galaxy samples.

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 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 study the impact of early dark energy (EDE) cosmologies on galaxy properties by coupling high-resolution numerical simulations with semi-analytic modeling (SAM) of galaxy formation and evolution. EDE models are characterized by a non-vanishing high-redshift contribution of dark energy, producing an earlier growth of structures and a modification of large-scale structure evolution. They can be viewed as typical representatives of non-standard dark energy models in which only the expansion history is modified, and hence the impact on galaxy formation is indirect. We show that in EDE cosmologies the predicted space density of galaxies is enhanced at all scales with respect to the standard LCDM scenario, and the corresponding cosmic star formation history and stellar mass density is increased at high-redshift. We compare these results with a set of theoretical predictions obtained with alternative SAMs applied to our reference LCDM simulation, yielding a rough measure of the systematic uncertainty of the models. We find that the modifications in galaxy properties induced by EDE cosmologies are of the same order of magnitude as intra-SAM variations for a standard LCDM realization (unless rather extreme EDE models are considered), suggesting that is difficult to use such predictions alone to disentangle between different cosmological scenarios. However, when independent information on the underlying properties of host dark matter haloes is included, the SAM predictions on galaxy bias may provide important clues on the expansion history and the equation-of-state evolution.

In this paper we discuss the contribution of different astrophysical sources to the ionization of neutral hydrogen at different redshifts. We critically revise the arguments in favour/against a substantial contribution of Active Galactic Nuclei (AGNs) and/or Lyman Break Galaxies (LBGs) to the reionization of the Universe at z>5. We consider extrapolations of the high-z QSO and LBG luminosity functions and their redshift evolution as well as indirect constraints on the space density of lower luminosity Active Galactic Nuclei based on the galaxy stellar mass function. Since the hypothesis of a reionization due to LBGs alone requires a significant contribution of faint dwarf galaxies and a LyC photon escape fraction (f_esc) of the order of ~20%, in tension with present observational constraints, we examine under which hypothesis AGNs and LBGs may provide a combined relevant contribution to the reionization. We show that a relatively steep faint-end of the AGN luminosity function, consistent with present constraints, provides a relevant (although sub-dominant) contribution, thus allowing us to recover the required ionizing photon rates with f_esc~5% up to z~7. At higher redshifts, we test the case for a luminosity-dependent f_esc scenario and we conclude that, if the observed LBGs are indeed characterized by very low f_esc, values of the order of f_esc~70% are needed for objects below our detection threshold, for this galaxy population to provide a substantial contribution to reionization. Clearly, the study of the properties of faint sources (both AGNs and LBGs) is crucial.

A serious concern for semi-analytical galaxy formation models, aiming to simulate multi-wavelength surveys and to thoroughly explore the model parameter space, is the extremely time consuming numerical solution of the radiative transfer of stellar radiation through dusty media. To overcome this problem, we have implemented an artificial neural network algorithm in the radiative transfer code GRASIL, in order to significantly speed up the computation of the infrared SED. The ANN we have implemented is of general use, in that its input neurons are defined as those quantities effectively determining the shape of the IR SED. Therefore, the training of the ANN can be performed with any model and then applied to other models. We made a blind test to check the algorithm, by applying a net trained with a standard chemical evolution model (i.e. CHE_EVO) to a mock catalogue extracted from the SAM MORGANA, and compared galaxy counts and evolution of the luminosity functions in several near-IR to sub-mm bands, and also the spectral differences for a large subset of randomly extracted models. The ANN is able to excellently approximate the full computation, but with a gain in CPU time by $\sim 2$ orders of magnitude. It is only advisable that the training covers reasonably well the range of values of the input neurons in the application. Indeed in the sub-mm at high redshift, a tiny fraction of models with some sensible input neurons out of the range of the trained net cause wrong answer by the ANN. These are extreme starbursting models with high optical depths, favorably selected by sub-mm observations, and difficult to predict a priori.

We used Early Release Science (ERS) observations taken with the Wide Field Camera 3 (WFC3) in the GOODS-S field to study the galaxy stellar mass function (GSMF) at 0.6<=z<4.5. Deep WFC3 near-IR data (for Y as faint as 27.3, J and H as faint as 27.4 AB mag at 5 sigma), as well as deep Ks (as faint as 25.5 at 5 sigma) Hawk-I band data, provide an exquisite data set with which determine in an unprecedented way the low-mass end of the GSMF, allowing an accurate probe of masses as low as M~7.6 10^9 Msun at z~3. Although the area used is relatively small (~33 arcmin^2), we found generally good agreement with previous studies on the entire mass range. Our results show that the slope of the faint-end increases with redshift, from alpha=-1.44+/-0.03 at z~0.8 to alpha=-1.86+/-0.16 at z~3, although indications exist that it does not steepen further between z~3 and z~4. This result is insensitive to any uncertainty in the M* parameter. The steepness of the GSMF faint-end solves the well-known disagreement between the stellar mass density (SMD) and the integrated star formation history at z>2. However, we confirm the that there appears to be an excess of integrated star formation with respect to the SMD at z<2, by a factor of ~2-3. Our comparison of the observations with theoretical predictions shows that the models forecast a greater abundance of low mass galaxies, at least up to z~3, as well as a dearth of massive galaxies at z~4 with respect to the data, and that the predicted SMD is generally overestimated at z<~2.

[abridged] We study the evolution of the Star Formation Rate Function (SFRF) of massive galaxies over the 0.4<z<1.8 redshift range and its implications for our understanding of the physical processes responsible for galaxy evolution. We use multiwavelength observations included in the GOODS-MUSIC catalogue, which provides a suitable coverage of the spectral region from 0.3 to 24 micron and either spectroscopic or photometric redshifts for each object. Individual SFRs have been obtained by combining UV and 24 micron observations, when the latter were available. For all other sources an "SED fitting" SFR estimate has been considered. We then define a stellar mass limited sample, complete in the Mstar>1.e10 Msun range and determine the SFRF using the 1/Vmax algorithm. We define simulated galaxy catalogues based on three different semi-analytical models of galaxy formation and evolution. We show that the theoretical SFRFs are well described by a double power law functional form and its redshift evolution is approximated with high accuracy by a pure evolution of the typical SFR. We find good agreement between model predictions and the high-SFR end of the SFRF, when the observational errors on the SFR are taken into account. However, the observational SFRF is characterised by a double peaked structure, which is absent in its theoretical counterparts. At z>1.0 the observed SFRF shows a relevant density evolution, which is not reproduced by SAMs, due to the well known overprediction of intermediate mass galaxies at z~2. The agreement at the low-SFR end is poor: all models overpredict the space density of SFR~1 Msun/yr and no model reproduces the double peaked shape of the observational SFRF. If confirmed by deeper IR observations, this discrepancy will provide a key constraint on theoretical modelling of star formation and stellar feedback.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

We study the correlation between the morphological mix of cluster galaxies and the assembly history of the parent cluster by taking advantage of two independently developed semi-analytic models for galaxy formation and evolution. In our models, both the number of cluster members and that of elliptical members increase as a function of cluster mass, in such a way that the resulting elliptical fractions are approximately independent of cluster mass. The population of cluster ellipticals exhibit a marked bimodal distribution as a function of galaxy stellar mass, with a dip at masses $\sim 10^{10}\,{\rm M}_{\odot}$. In the framework of our models, this bimodality originates from the combination of a strongly decreasing number of galaxies with increasing stellar mass, and a correspondingly increasing probability of experiencing major mergers. We show that the correlation between the measured elliptical fraction and the assembly history of the parent cluster is weak, and that it becomes stronger in models that adopt longer galaxy merger times. We argue that this results from the combined effect of a decreasing bulge production due to a reduced number of mergers, and an increasing survival probability of pre-existing ellipticals, with the latter process being more important than the former.

It is generally assumed that the central galaxy in a dark matter halo, that is, the galaxy with the lowest specific potential energy, is also the brightest halo galaxy (BHG), and that it resides at rest at the centre of the dark matter potential well. This central galaxy paradigm (CGP) is an essential assumption made in various fields of astronomical research. In this paper we test the validity of the CGP using a large galaxy group catalogue constructed from the Sloan Digital Sky Survey. For each group we compute two statistics, ${\cal R}$ and ${\cal S}$, which quantify the offsets of the line-of-sight velocities and projected positions of brightest group galaxies relative to the other group members. By comparing the cumulative distributions of $|{\cal R}|$ and $|{\cal S}|$ to those obtained from detailed mock group catalogues, we rule out the null-hypothesis that the CGP is correct. Rather, the data indicate that in a non-zero fraction $f_{\rm BNC}(M)$ of all haloes of mass $M$ the BHG is not the central galaxy, but instead, a satellite galaxy. In particular, we find that $f_{\rm BNC}$ increases from $\sim 0.25$ in low mass haloes ($10^{12} h^{-1} {\rm M_{\odot}} \leq M \lsim 2 \times 10^{13} h^{-1}{\rm M_{\odot}}$) to $\sim 0.4$ in massive haloes ($M \gsim 5 \times 10^{13} h^{-1} {\rm M_{\odot}}$). We show that these values of $f_{\rm BNC}$ are uncomfortably high compared to predictions from halo occupation statistics and from semi-analytical models of galaxy formation. We end by discussing various implications of a non-zero $f_{\rm BNC}(M)$, with an emphasis on the halo masses inferred from satellite kinematics.

Interstellar dust grains are responsible for modifying the spectral energy distribution (SED) of galaxies, both absorbing starlight at UV and optical wavelengths and converting this energy into thermal emission in the infrared. The detailed description of these phenomena is of fundamental importance in order to compare the predictions of theoretical models of galaxy formation and evolution with the most recent observations in the infrared region. In this paper we compare the results of GRASIL, a code explicitly solving for the equation of radiative transfer in a dusty medium, with the predictions of a variety of IR template libraries. We employ star formation history samples extracted from the semi-analytical galaxy formation model MORGANA to create libraries of synthetic SEDs from the near- to the far-infrared. We consider model predictions at different redshift ranges to explore any possible influence in the shape and normalization of the SEDs due to the expected evolution of the galaxy properties. We compute the total absorbed starlight predicted by GRASIL at optical wavelengths to statistically compare the synthetic SEDs with the selected IR templates. We show that synthetic SEDs at a given total infrared luminosity are predicted to be systematically different at different redshift and for different properties of the underlying model galaxy. However, we determine spectral regions where the agreement between the results of radiative transfer and IR templates is good in a statistical sense (i.e. in terms of the luminosity functions). Moreover, we highlight some potentially relevant discrepancies between the different approaches, both in the region dominated by PAH emission and at sub-mm wavelengths. These results determine potentially critical issues in the infrared luminosity functions as predicted by semi-analytical models coupled with different IR flux estimators.

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.

AGN feedback is believed to play an important role in shaping a variety of observed galaxy properties, as well as the evolution of their stellar masses and star formation rates. In particular, in the current theoretical paradigm of galaxy formation, AGN feedback is believed to play a crucial role in regulating the levels of activity in galaxies, in relatively massive halos at low redshift. Only in recent years, however, has detailed statistical information on the dependence of galaxy activity on stellar mass, parent halo mass and hierarchy has become available. In this paper, we compare the fractions of galaxies belonging to different activity classes (star-forming, AGN and radio active) with predictions from four different and independently developed semi-analytical models. We adopt empirical relations to convert physical properties into observables (H_alpha emission lines, OIII line strength and radio power). We demonstrate that all models used in this study reproduce the overall distributions of galaxies belonging to different activity classes as a function of stellar mass and halo mass: star forming galaxies and the strongest radio sources are preferentially associated with low-mass and high-mass galaxies/halos respectively. However, model predictions differ from observational measurements in a number of ways. All models used in our study predict that almost every >1.e12 Msun dark matter halo and/or >1.e11 Msun galaxy should host a bright radio source, while only a small fraction of galaxies belong to this class in the data. In addition, radio brightness is expected to depend strongly on the mass of the parent halo mass in the models, while strong and weak radio galaxies are found in similar environments in data. Our results highlight that the distribution of AGN as a function of stellar mass provides one of the most promising discriminants between different gas accretion schemes.

We use ultra-deep ultraviolet VLT/VIMOS intermediate-band and VLT/FORS1 narrow-band imaging in the GOODS Southern field to derive limits on the distribution of the escape fraction (f_esc) of ionizing radiation for L >~ L*(z=3) Lyman Break Galaxies (LBGs) at redshift 3.4--4.5. Only one LBG, at redshift z=3.795, is detected in its Lyman continuum (LyC; S/N~5.5), the highest redshift galaxy currently known with a direct detection. Its ultraviolet morphology is quite compact (R_eff=0.8, kpc physical). Three out of seven AGN are also detected in their LyC, including one at redshift z=3.951 and z850 = 26.1. From stacked data (LBGs) we set an upper limit to the average f_esc in the range 5%--20%, depending on the how the data are selected (e.g., by magnitude and/or redshift). We undertake extensive Monte Carlo simulations that take into account intergalactic attenuation, stellar population synthesis models, dust extinction and photometric noise in order to explore the moments of the distribution of the escaping radiation. Various distributions (exponential, log-normal and Gaussian) are explored. We find that the median f_esc is lower than ~6% with an 84% percentile limit not larger than 20%. If this result remains valid for fainter LBGs down to current observational limits, then the LBG population might be not sufficient to account for the entire photoionization budget at the redshifts considered here, with the exact details dependent upon the assumed ionizing background and QSO contribution thereto. It is possible that f_esc depends on the UV luminosity of the galaxies, with fainter galaxies having higher f_esc, and estimates of f_esc from a sample of faint LBG from the HUDF (i775<28.5) are in broad quantitative agreement with such a scenario.

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 test the luminosity function of Milky Way satellites as a constraint for the nature of Dark Matter particles. We perform dissipationless high-resolution N-body simulations of the evolution of Galaxy-sized halo in the standard Cold Dark Matter (CDM) model and in four Warm Dark Matter (WDM) scenarios, with a different choice for the WDM particle mass (m_w). We then combine the results of the numerical simulations with semi-analytic models for galaxy formation, to infer the properties of the satellite population. Quite surprisingly we find that even WDM models with relatively low m_w values (2-5 keV) are able to reproduce the observed abundance of ultra faint (Mv<-9) dwarf galaxies, as well as the observed relation between Luminosity and mass within 300 pc. Our results suggest a lower limit of 1 keV for thermal warm dark matter, in broad agreement with previous results from other astrophysical observations like Lyman-alpha forest and gravitational lensing.

[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]

Using a large SDSS galaxy group catalogue, we study how the stellar ages and metallicities of central and satellite galaxies depend on stellar mass and halo mass. We find that satellites are older and metal-richer than centrals of the same stellar mass. In addition, the slopes of the age-stellar mass and metallicity-stellar mass relations are found to become shallower in denser environments. This is due to the fact that the average age and metallicity of low mass satellite galaxies increase with the mass of the halo in which they reside. A comparison with the semi-analytical model of Wang et al. (2008) shows that it succesfully reproduces the fact that satellites are older than centrals of the same stellar mass and that the age difference increases with the halo mass of the satellite. This is a consequence of strangulation, which leaves the stellar populations of satellites to evolve passively, while the prolonged star formation activity of centrals keeps their average ages younger. The resulting age offset is larger in more massive environments because their satellites were accreted earlier. The model fails, however, in reproducing the halo mass dependence of the metallicities of low mass satellites, yields metallicity-stellar mass and age-stellar mass relations that are too shallow, and predicts that satellite galaxies have the same metallicities as centrals of the same stellar mass, in disagreement with the data. We argue that these discrepancies are likely to indicate the need to (i) modify the recipes of both supernova feedback and AGN feedback, (ii) use a more realistic description of strangulation, and (iii) include a proper treatment of the tidal stripping, heating and destruction of satellite 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.

The star formation rate (SFR) and black hole accretion rate (BHAR) functions are measured to be proportional to each other at z < ~3. This close correspondence between SF and BHA would naturally yield a BH mass-galaxy mass correlation, whereas a BH mass-bulge mass correlation is observed. To explore this apparent contradiction we study the SF in spheroid-dominated galaxies between z=1 and the present day. We use 903 galaxies from the COMBO-17 survey with M* >2x10^10M_sun, ultraviolet and infrared-derived SFRs from Spitzer and GALEX, and morphologies from GEMS HST/ACS imaging. Using stacking techniques, we find that <25% of all SF occurs in spheroid-dominated galaxies (Sersic index n>2.5), while the BHAR that we would expect if the global scalings held is three times higher. This rules out the simplest picture of co-evolution, in which SF and BHA trace each other at all times. These results could be explained if SF and BHA occur in the same events, but offset in time, for example at different stages of a merger event. However, one would then expect to see the corresponding star formation activity in early-stage mergers, in conflict with observations. We conclude that the major episodes of SF and BHA occur in different events, with the bulk of SF happening in isolated disks and most BHA occurring in major mergers. The apparent global co-evolution results from the regulation of the BH growth by the potential well of the galactic spheroid, which includes a major contribution from disrupted disk stars.

We study the star formation and the mass assembly process of 0.3<=z<2.5 galaxies using their IR emission from MIPS 24um band. We used an updated version of the GOODS-MUSIC catalog, extended by the addition of mid-IR fluxes. We compared two different estimators of the Star Formation Rate: the total infrared emission derived from 24um, estimated using both synthetic and empirical IR templates, and the multiwavelength fit to the full galaxy SED. For both estimates, we computed the SFR Density and the Specific SFR. The two SFR tracers are roughly consistent, given the uncertainties involved. However, they show a systematic trend, IR-based estimates exceeding the fit-based ones as the SFR increases. We show that: a) at z>0.3, the SFR is well correlated with stellar mass, and this relationship seems to steepen with redshift (using IR-based SFRs); b) the contribution to the global SFRD by massive galaxies increases with redshift up to ~2.5, more rapidly than for galaxies of lower mass, but appears to flatten at higher z; c) despite this increase, the most important contributors to the SFRD at any z are galaxies of about, or immediately lower than, the characteristic stellar mass; d) at z~2, massive galaxies are actively star-forming, with a median SFR 300 Msun/yr. During this epoch, they assemble a substantial part of their final stellar mass; e) the SSFR shows a clear bimodal distribution. The analysis of the SFRD and the SSFR seems to support the downsizing scenario, according to which high mass galaxies have formed their stars earlier and faster than their low mass counterparts. A comparison with theoretical models indicates that they follow the global increase in the SSFR with redshift and predict the existence of quiescent galaxies even at z>1.5, but they systematically underpredict the average SSFR.

Deep multicolour surveys are the main tool to explore the formation and evolution of the faint galaxies which are beyond the spectroscopic limit with the present technology. The photometric properties of these faint galaxies are usually compared with current renditions of semianalytical models to provide constraints on the fundamental physical processes involved in galaxy formation and evolution, namely the mass assembly and the star formation. Galaxy counts over large sky areas in the near-UV band are important because they are difficult to obtain given the low efficiency of near-UV instrumentation, even at 8m class telescopes. A large instrumental field of view helps in minimizing the biases due to the cosmic variance. We have obtained deep images in the 360nm U band provided by the blue channel of the Large Binocular Camera at the prime focus of the Large Binocular Telescope. We have derived over an area of ~0.4 sq. deg. the galaxy number counts down to U=27 in the Vega system (corresponding to U=27.86 in the AB system) at a completeness level of 30% reaching the faintest current limit for this wavelength and sky area. The shape of the galaxy counts in the U band can be described by a double power-law, the bright side being consistent with the shape of shallower surveys of comparable or greater areas. The slope bends over significantly at U>23.5 ensuring the convergence of the contribution by star forming galaxies to the EBL in the near-UV band to a value which is more than 70% of the most recent upper limits derived for this band. We have jointly compared our near-UV and K band counts collected from the literature with few selected hierarchical CDM models emphasizing critical issues in the physical description of the galaxy formation and evolution.

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.

The QSO luminosity function at z>5 provides strong constraints on models of joint evolution of QSO and their hosts. However, these observations are challenging because the low space densities of these objects necessitate surveying of large areas, in order to obtain statistically meaningful samples, while at the same time cosmological redshifting and dimming means that rather deep Near Infrared (NIR) imaging must be carried out. Several upcoming and proposed facilities with wide-field NIR imaging capabilities will open up this new region of parameter space. In this paper we present predictions for the expected number counts of z>5 QSOs, based on simple empirical and semi-empirical models of QSO evolution, as a function of redshift, depth and surveyed area. We compute the evolution of observed-frame QSO magnitudes and colors in a representative photometric system covering the wavelength range 550nm<lambda<1800nm, and combine this information with different estimates for the evolution of the QSO luminosity function. We conclude that planned ground-based surveys such as Pan-STARRS and VISTA should be able to detect a large number of luminous QSOs up to z<7.5, but that space-based missions such as EUCLID (formerly SPACE/DUNE) or SNAP are probably required in order to obtain substantial samples at higher redshift. We also use our models to predict the expected number counts for future X-ray space missions (such as XEUS and Constellation-X), and show that because of their small field-of-view, these telescopes are unlikely to discover significant numbers of AGN at very high redshift. However, X-ray follow-up of objects detected at longer wavelength will be an important means of confirming their identity as AGN and constraining obscuration.

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.

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 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.

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]

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.