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 examine the scaling relations between the mass of a supermassive black hole (SMBH) and its host galaxy properties at $1.2<z<1.7$ using both observational data and simulations. Recent measurements of 32 X-ray-selected broad-line Active Galactic Nucleus (AGNs) are compared with two independent state-of-the-art efforts, including the hydrodynamic simulation MassiveBlackII (MBII) and a semi-analytic model (SAM). After applying an observational selection function to the simulations, we find that both MBII and SAM agree well with the data, in terms of the central distribution. However, the dispersion in the mass ratio between black hole mass and stellar mass is significantly more consistent with the MBII prediction ($\sim0.3~$dex), than with the SAM ($\sim0.7~$dex), even when accounting for observational uncertainties. Hence, our observations can distinguish between the different recipes adopted in the models. The mass relations in the MBII are highly dependent on AGN feedback while the relations in the SAM are more sensitive to galaxy merger events triggering nuclear activity. Moreover, the intrinsic scatter in the mass ratio of our high-$z$ sample is comparable to that observed in the local sample, all but ruling out the proposed scenario the correlations are purely stochastic in nature arising from some sort of cosmic central limit theorem. Our results support the hypothesis of AGN feedback being responsible for a causal link between the SMBH and its host galaxy, resulting in a tight correlation between their respective masses.

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized.

The mass and structural assembly of galaxies is a matter of intense debate. Current theoretical models predict the existence of a linear relationship between galaxy size ($R_e$) and the host dark matter halo virial radius ($R_h$).\\ By making use of semi-empirical models compared to the size distributions of central galaxies from the Sloan Digital Sky Survey, we provide robust constraints on the normalization and scatter of the $R_e-R_h$ relation. We explore the parameter space of models in which the $R_e-R_h$ relation is mediated by either the spin parameter or the concentration of the host halo, or a simple constant the nature of which is in principle unknown. We find that the data require extremely tight relations for both early-type and late-type galaxies (ETGs,LTGs), especially for more massive galaxies. These constraints challenge models based solely on angular momentum conservation, which predict significantly wider distributions of galaxy sizes and no trend with stellar mass, if taken at face value. We discuss physically-motivated alterations to the original models that bring the predictions into better agreement with the data. We argue that the measured tight size distributions of SDSS disk galaxies can be reproduced by semi-empirical models in which the $R_e-R_h$ connection is mediated by the \emph{stellar} specific angular momenta $j_{star}.$ We find that current cosmological models of galaxy formation broadly agree with our constraints for LTGs, and justify the strong link between $R_e$ and $j_{star}$ that we propose, however the tightness of the $R_e-R_h$ relation found in such ab-initio theoretical models for ETGs is in tension with our semi-empirical findings.

Observational systematics complicate comparisons with theoretical models limiting understanding of galaxy evolution. In particular, different empirical determinations of the stellar mass function imply distinct mappings between the galaxy and halo masses, leading to diverse galaxy evolutionary tracks. Using our state-of-the-art STatistical sEmi-Empirical modeL, STEEL, we show fully self-consistent models capable of generating galaxy growth histories that simultaneously and closely agree with the latest data on satellite richness and star-formation rates at multiple redshifts and environments. Central galaxy histories are generated using the central halo mass tracks from state-of-the-art statistical dark matter accretion histories coupled to abundance matching routines. We show that too flat high-mass slopes in the input stellar-mass-halo-mass relations as predicted by previous works, imply non-physical stellar mass growth histories weaker than those implied by satellite accretion alone. Our best-fit models reproduce the satellite distributions at the largest masses and highest redshifts probed, the latest data on star formation rates and its bi-modality in the local Universe, and the correct fraction of ellipticals. Our results are important to predict robust and self-consistent stellar-mass-halo-mass relations and to generate reliable galaxy mock catalogues for the next generations of extra-galactic surveys such as Euclid and LSST.

A supermassive black hole has been found at the centre of nearly every galaxy observed with sufficient sensitivity. The masses of these black holes are observed to increase with either the total mass or the mean (random) velocity of the stars in their host galaxies. The origin of these correlations remains elusive. Observational systematics and biases severely limit our knowledge of the local demography of supermassive black holes thus preventing accurate model comparisons and progress in this field. Here we show that the large-scale spatial distribution of local active galactic nuclei (AGN), believed to be accreting supermassive black holes, can constrain the shape and normalization of the black hole-stellar mass relation thus bypassing resolution-related observational biases. In turn, our results can set more stringent constraints on the so-called "radiative efficiency", a fundamental parameter describing the inner physics of supermassive black holes that is closely linked to their spin, geometry, and ability to release energy. The mean value of the radiative efficiency can be estimated by comparing the average total luminous output of AGN with the relic mass density locked up in quiescent supermassive black holes at galaxy centres today. For currently accepted values of the AGN obscured fractions and bolometric corrections, our newest estimates of the local supermassive black hole mass density favour mean radiative efficiencies of ~10-20\%, suggesting that the vast majority of supermassive black holes are spinning moderately to rapidly. With large-scale AGN surveys coming online, our novel methodology will enable even tighter constraints on the fundamental parameters that regulate the growth of supermassive black holes.

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.

The study of the space density of bright AGNs at $z>4$ has been subject to extensive effort given its importance for the estimate of the cosmological ionizing emissivity and growth of supermassive black holes. In this context we have recently derived high space densities of AGNs at $z\sim 4$ and $-25<M_{1450}<-23$ in the COSMOS field from a spectroscopically complete sample. In the present paper we attempt to extend the knowledge of the AGN space density at fainter magnitudes ($-22.5<M_{1450}<-18.5$) in the $4<z<6.1$ redshift interval by means of a multiwavelength sample of galaxies in the CANDELS GOODS-South, GOODS-North and EGS fields. We use an updated criterion to extract faint AGNs from a population of NIR (rest-frame UV) selected galaxies at photometric $z>4$ showing X-ray detection in deep Chandra images available for the three CANDELS fields. We have collected a photometric sample of 32 AGN candidates in the selected redshift interval, six of which having spectroscopic redshifts. Including our COSMOS sample as well as other bright QSO samples allows a first guess on the shape of the UV luminosity function at $z\sim 4.5$. The resulting emissivity and photoionization rate appear consistent with that derived from the photoionization level of the intergalactic medium at $z\sim 4.5$. An extrapolation to $z\sim 5.6$ suggests an important AGN contribution to the IGM ionization if there are no significant changes in the shape of the UV luminosity function.

Recent advances in observations have provided a wealth of measurements of the expansions of outflows in galactic discs out to large radii in a variety of galactic hosts. To provide an updated baseline for the interpretation of such data, and to assess to what extent the present status of the modeling is consistent with the existing observations, we provide a compact two-dimensional description for the expansion of AGN-driven shocks in realistic galactic discs with exponential gas density profiles in a disc geometry. We derive solutions for the outflow expansion and the mass outflow rates in different directions with respect to the plane of the disc. These are expressed in terms of the global properties of the host galaxy and of the central AGN to allow for an easy and direct comparison with existing observations in a variety of galactic hosts with measured properties, and out to distances $\sim 10$ kpc from the centre. The results are compared with a state-of-the-art compilation of observed outflows in 19 galaxies with different measured gas and dynamical mass, allowing for a detailed, one-by-one comparison with the model predictions. The agreement we obtain for a wide range of host galaxy gas mass and AGN bolometric luminosity provides a quantitative systematic test for the modeling of AGN-driven outflows in galactic discs. We also consider a larger sample of galaxies with no reliable measurements of the gas and dynamical mass. In this case we perform a comparison of the model predictions for different bins of AGN luminosities assuming values for the gas mass and dynamical mass derived from scaling relations. The encouraging, quantitative agreement of the model predictions with a wide set of existing observations constitutes a baseline for the interpretation of forthcoming data, and for a more detailed treatment of AGN feedback in galaxy formation models.

The investigation of distant galaxy formation and evolution is a powerful tool to constrain dark matter scenarios, supporting and in some cases surpassing other astrophysical and experimental probes. The recent completion of the Hubble Frontier Fields (HFF) programme combining ultra-deep Hubble Space Telescope observations and the magnification power of gravitational lensing produced by foreground galaxy clusters has enabled the detection of the faintest primordial galaxies ever studied. Here we show how the number density of such primordial galaxies allows to constrain a variety of DM models alternative to CDM. In particular, it provides stringent limits on the mass of thermal WDM candidates, on the parameter space of sterile neutrino production models, and on other DM scenarios featuring particles in the keV mass range which is also supported by recent detections of a 3.5keV X-ray line. These constraints are robust and independent of the baryonic physics modeling of galaxy formation and evolution. Fuzzy DM (ultralight DM particles) results strongly disfavored.

In the local Universe, the existence of very young galaxies (VYGs), having formed at least half their stellar mass in the last 1 Gyr, is debated. We predict the present-day fraction of VYGs among central galaxies as a function of galaxy stellar mass. For this, we apply to high mass resolution Monte-Carlo halo merger trees (MCHMTs) three (one) analytical models of galaxy formation, where the ratio of stellar to halo mass (mass growth rate) is a function of halo mass and redshift. Galaxy merging is delayed until orbital decay by dynamical friction. With starbursts associated with halo mergers, our models predict typically one percent of VYGs up to galaxy masses of $10^{10}$ M$_\odot$, falling rapidly at higher masses, and VYGs are usually associated with recent major mergers of their haloes. Without these starbursts, two of the models have VYG fractions reduced by 1 or 2 dex at low or intermediate stellar masses, and VYGs are rarely associated with major halo mergers. In comparison, the state-of-the-art semi-analytical model (SAM) of Henriques et al. produces only 0.01% of VYGs at intermediate masses. Finally, the Menci et al. SAM run on MCMHTs with Warm Dark Matter cosmology generates 10 times more VYGs at masses below $10^8$ M$_\odot$ than when run with Cold Dark Matter. The wide range in these VYG fractions illustrates the usefulness of VYGs to constrain both galaxy formation and cosmological models.

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.

We investigate galaxy formation in models with dark matter (DM) constituted by sterile neutrinos. Given their large parameter space, defined by the combinations of sterile neutrino mass $m_{\nu}$ and mixing parameter $\sin^2(2\theta)$ with active neutrinos, we focus on models with $m_{\nu}=7$ keV, consistent with the tentative 3.5 keV line detected in several X-ray spectra of clusters and galaxies. We consider i) two resonant production models with $\sin^2(2\theta)=5\,10^{-11}$ and $\sin^2(2\theta)=2\,10^{-10}$, to cover the range of mixing parameter consistent with the 3.5 keV line; ii) two scalar-decay models, representative of the two possible cases characterizing such a scenario: a freeze-in and a freeze-out case. We also consider thermal Warm Dark Matter with particle mass $m_X=3$ keV. Using a semi-analytic model, we compare the predictions for the different DM scenarios with a wide set of observables. We find that comparing the predicted evolution of the stellar mass function, the abundance of satellites of Milky Way-like galaxies, and the global star formation history of galaxies with observations does not allow to disentangle the effects of the baryonic physics from those related to the different DM models. On the other hand, the distribution of the stellar-to-halo mass ratios, the abundance of faint galaxies in the UV luminosity function at $z\gtrsim 6$, and the specific star formation and age distribution of local, low-mass galaxies constitute potential probes for the considered DM scenarios. We discuss how next observations with upcoming facilities will enable to rule out or to strongly support DM models based on sterile neutrinos.

We have performed a high sensitivity observation of the UFO/BAL quasar APM 08279+5255 at z=3.912 with NOEMA at 3.2 mm, aimed at detecting fast moving molecular gas. We report the detection of blueshifted CO(4-3) with maximum velocity (v95\%) of $-1340$ km s$^{-1}$, with respect to the systemic peak emission, and a luminosity of $L' = 9.9\times 10^9 ~\mu^{-1}$ K km s$^{-1}$ pc$^{-2}$ (where $\mu$ is the lensing magnification factor). We discuss various scenarios for the nature of this emission, and conclude that this is the first detection of fast molecular gas at redshift $>3$. We derive a mass flow rate of molecular gas in the range $\rm \dot M=3-7.4\times 10^3$ M$_\odot$/yr, and momentum boost $\dot P_{OF} / \dot P_{AGN} \sim 2-6$, therefore consistent with a momentum conserving flow. For the largest $\dot P_{OF}$ the scaling is also consistent with a energy conserving flow with an efficiency of $\sim$10-20\%. The present data can hardly discriminate between the two expansion modes. The mass loading factor of the molecular outflow $\eta=\dot M_{OF}/SFR$ is $>>1$. We also detect a molecular emission line at a frequency of 94.83 GHz, corresponding to a rest frame frequency of 465.8 GHz, which we tentatively identified with the cation molecule $\rm N_2H^+$(5-4), which would be the first detection of this species at high redshift. We discuss the alternative possibility that this emission is due to a CO emission line from the, so far undetected, lens galaxy. Further observations of additional transitions of the same species with NOEMA can discriminate between the two scenarios.

We derive the contribution to the extragalactic gamma-ray background (EGB) from AGN winds and star-forming galaxies by including a physical model for the gamma-ray emission produced by relativistic protons accelerated by AGN-driven and supernova-driven shocks into a state-of-the-art semi-analytic model of galaxy formation. This is based on galaxy interactions as triggers of AGN accretion and starburst activity and on expanding blast wave as the mechanism to communicate outwards the energy injected into the interstellar medium by the active nucleus. We compare the model predictions with the latest measurement of the EGB spectrum performed by the Fermi-LAT in the range between 100 MeV and 820 GeV. We find that AGN winds can provide ~35$\pm$15% of the observed EGB in the energy interval E_{\gamma}=0.1-1 GeV, for ~73$\pm$15% at E_{\gamma}=1-10 GeV, and for ~60$\pm$20% at E_{\gamma}>10 GeV. The AGN wind contribution to the EGB is predicted to be larger by a factor of 3-5 than that provided by star-forming galaxies (quiescent plus starburst) in the hierarchical clustering scenario. The cumulative gamma-ray emission from AGN winds and blazars can account for the amplitude and spectral shape of the EGB, assuming the standard acceleration theory, and AGN wind parameters that agree with observations. We also compare the model prediction for the cumulative neutrino background from AGN winds with the most recent IceCube data. We find that for AGN winds with accelerated proton spectral index p=2.2-2.3, and taking into account internal absorption of gamma-rays, the Fermi-LAT and IceCube data could be reproduced simultaneously.

We report the discovery of 11 very faint (r< 23), low surface brightness ({\mu}_r< 27 mag/arcsec^2) dwarf galaxies in one deep field in the Virgo cluster, obtained by the prime focus cameras (LBC) at the Large Binocular Telescope (LBT). These extend our previous sample to reach a total number of 27 galaxies in a field of just of 0.17 deg^2 located at a median distance of 390 kpc from the cluster center. Their association with the Virgo cluster is supported by their separate position in the central surface brightness - total magnitude plane with respect to the background galaxies of similar total magnitude. For a significant fraction (26\%) of the sample the association to the cluster is confirmed by spectroscopic follow-up. We show that the mere abundance of satellite galaxies corresponding to our observed number in the target field provides extremely tight constraints on Dark Matter models with suppressed power spectrum compared to the Cold Dark Matter case, independently of the galaxy luminosity distribution. In particular, requiring the observed number of satellite galaxies not to exceed the predicted abundance of Dark Matter sub-halos yields a limit m_X >3 keV at 1-{\sigma} and m_X > 2.3 keV at 2-{\sigma} confidence level for the mass of thermal Warm Dark Matter particles. Such a limit is competitive with other limits set by the abundance of ultra-faint satellite galaxies in the Milky Way, is completely independent of baryon physics involved in galaxy formation, and has the potentiality for appreciable improvements with next observations. We extend our analysis to Dark Matter models based on sterile neutrinos, showing that our observations set tight constraints on the combination of sterile neutrino mass m_{\nu} and mixing parameter sin^2(2{\theta}). We discuss the robustness of our results with respect to systematics.

Feedback from accreting SMBHs is often identified as the main mechanism responsible for regulating star-formation in AGN host galaxies. However, the relationships between AGN activity, radiation, winds, and star-formation are complex and still far from being understood. We study scaling relations between AGN properties, host galaxy properties and AGN winds. We then evaluate the wind mean impact on the global star-formation history, taking into account the short AGN duty cycle with respect to that of star-formation. We first collect AGN wind observations for 94 AGN with detected massive winds at sub-pc to kpc spatial scales. We then fold AGN wind scaling relations with AGN luminosity functions, to evaluate the average AGN wind mass-loading factor as a function of cosmic time. We find strong correlations between the AGN molecular and ionised wind mass outflow rates and the AGN bolometric luminosity. The power law scaling is steeper for ionised winds (slope 1.29+/-0.38) than for molecular winds (0.76+/-0.06), meaning that the two rates converge at high bolometric luminosities. The molecular gas depletion timescale and the molecular gas fraction of galaxies hosting powerful AGN winds are 3-10 times shorter and smaller than those of main-sequence galaxies with similar SFR, stellar mass and redshift. These findings suggest that, at high AGN bolometric luminosity, the reduced molecular gas fraction may be due to the destruction of molecules by the wind, leading to a larger fraction of gas in the atomic ionised phase. The AGN wind mass-loading factor $\eta=\dot M_{OF}/SFR$ is systematically higher than that of starburst driven winds. Our analysis shows that AGN winds are, on average, powerful enough to clean galaxies from their molecular gas only in massive systems at z<=2, i.e. a strong form of co-evolution between SMBHs and galaxies appears to break down for the least massive galaxies.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

We show that the measured abundance of ultra-faint lensed galaxies at $z\approx 6$ in the Hubble Frontier Fields (HFF) provides stringent constraints on the parameter space of i) Dark Matter models based on keV sterile neutrinos; ii) the "fuzzy" wavelike Dark Matter models, based on Bose-Einstein condensate of ultra-light particles. For the case of the sterile neutrinos, we consider two production mechanisms: resonant production through the mixing with active neutrinos and the decay of scalar particles. For the former model, we derive constraints for the combination of sterile neutrino mass $m_{\nu}$ and mixing parameter $\sin^2(2\theta)$ which provide the tightest lower bounds on the mixing angle (and hence on the lepton asymmetry) derived so far by methods independent of baryonic physics. For the latter we compute the allowed combinations of the scalar mass, its coupling to the Higgs field, and the Yukawa coupling of the scalar to the sterile neutrinos. We compare our results to independent, existing astrophysical bounds on sterile neutrinos in the same mass range. For the case of "fuzzy" Dark Matter, we show that the observed number density $\approx 1/{\rm Mpc}^3$ of high-redshift galaxies in the HFF sets a lower limit $m_\psi\geq 8\cdot 10^{-22}$ eV (at 3-$\sigma$ confidence level) on the particle mass, a result that strongly disfavors wavelike bosonic Dark Matter as a viable model for structure formation. We discuss the impact on our results of uncertainties due to systematics in the selection of highly magnified, faint galaxies at high redshifts.

We compute the non-thermal emissions produced by relativistic particles accelerated by the AGN-driven shocks in NGC 1068, and we compare the model predictions with the observed gamma-ray and radio spectra . The former is contributed by pion decay, inverse Compton scattering, and bremsstrahlung, while the latter is produced by synchrotron radiation. We derive the gamma-ray and radio emissions by assuming the standard acceleration theory, and we discuss how our results compare with those corresponding to other commonly assumed sources of gamma-ray and radio emissions, like Supernova remnants (SNR) or AGN jets. We find that the AGN-driven shocks observed in the circumnuclear molecular disk of such a galaxy provide a contribution to the gamma-ray emission comparable to that provided by the starburst activity when standard particle acceleration efficiencies are assumed, while they can yield the whole gamma-ray emission only when the parameters describing the acceleration efficiency and the proton coupling with the molecular gas are tuned to values larger than those assumed in standard, SNR-driven shocks. We discuss the range of acceleration efficiencies (for protons and electrons) and of proton calorimetric fractions required to account for the observed gamma-ray emission in the AGN outflow model. We further compare the neutrino flux expected in our model with constraints from current experiments, and we provide predictions for the detections by the upcoming KM3NeT neutrino telescope. This analysis strongly motivates observations of NGC 1068 at >TeV energies with current and future Cherenkov telescopes in order to gain insight into the nature of the gamma-rays source.

We show that the recently measured UV luminosity functions of ultra-faint lensed galaxies at z= 6 in the Hubble Frontier Fields provide an unprecedented probe for the mass m_X of the Warm Dark Matter candidates independent of baryonic physics. Comparing the measured abundance of the faintest galaxies with the maximum number density of dark matter halos in WDM cosmologies sets a robust limit m_X> 2.9 keV for the mass of thermal relic WDM particles at a 1-sigma confidence level, m_X> 2.4 keV at 2-sigma, and m_X> 2.1 keV at 3-sigma. These constitute the tightest constraints on WDM particle mass derived to date independently of the baryonic physics involved in galaxy formation. We discuss the impact of our results on the production mechanism of sterile neutrinos. In particular, if sterile neutrinos are responsible for the 3.5 keV line reported in observations of X-ray clusters, our results firmly rule out the Dodelson-Widrow production mechanism, and yield m_{sterile}> 6.1 keV for sterile neutrinos produced via the Shi-Fuller mechanism.

Recent observations indicate that the mass of Supermassive Black Holes (SMBHs) correlate differently with different galaxy stellar components. Comparing such observations with the results of "ab initio" galaxy formation models can provide insight on the mechanisms leading to the growth of SMBHs. Here we use a state-of-the-art semi-analytic model of galaxy formation to investigate the correlation of the different galaxy stellar components with the mass of the central SMBH. The stellar mass in the disc, in the bulge, and in the pseudo-bulge of galaxies is related to quiescent star formation, to galaxy interactions, and to the loss of angular momentum following disc instabilities, respectively. Consistently with recent findings, we find that while the predicted bulge masses are tightly correlated with the SMBH masses, the correlation between the latter and the galactic discs shows a much larger scatter, in particular when bulgeless galaxies are considered. In addition, we obtain that the predicted masses of pseudo-bulges shows little or no-correlation with the masses of SMBHs. We track the histories of merging, star formation, and SMBH accretion to investigate the physical processes at the origin of such findings within the context of cosmological models of galaxy formation. Finally, we discuss the effects of variations of our assumed fiducial model on the results.

It is widely believed that the large discrepancy between the observed number of satellite galaxies and the predicted number of dark subhalos can be resolved via a variety of baryonic effects which suppress star formation in low mass halos.Supporting this hypothesis, numerous high resolution simulations with star formation, and associated feedback have been shown to reproduce the satellite luminosity function around Milky Way-mass simulated galaxies at redshift zero. However, a more stringent test of these models is their ability to simultaneously match the satellite luminosity functions of a range of host halo masses and redshifts. In this work we measure the luminosity function of faint (sub-Small Magellanic Cloud luminosity) satellites around hosts with stellar masses 10.5$<\log_{10}$M$_*$/M$_\odot<11.5$ to an unprecedented redshift of 1.5. This new measurement of the satellite luminosity function provides powerful new constraining power; we compare these results with predictions from four different simulations and show that although the models perform similarly over-all, no one model reproduces the satellite luminosity function reliably at all redshifts and host stellar masses. This result highlights the continued need for improvement in understanding the fundamental physics that governs satellite galaxy evolution.

We compute the mass function of galactic dark matter halos for different values of the Warm Dark Matter (WDM) particle mass m_X and compare it with the abundance of ultra-faint galaxies derived from the deepest UV luminosity function available so far at redshift z~2. The magnitude limit M_UV=-13 reached by such observations allows us to probe the WDM mass functions down to scales close to or smaller than the half-mass mode mass scale ~10^9 M_sun. This allowed for an efficient discrimination among predictions for different m_X which turn out to be independent of the star formation efficiency adopted to associate the observed UV luminosities of galaxies to the corresponding dark matter masses. Adopting a conservative approach to take into account the existing theoretical uncertainties in the galaxy halo mass function, we derive a robust limit m_X>1.8 keV for the mass of thermal relic WDM particles when comparing with the measured abundance of the faintest galaxies, while m_X>1.5 keV is obtained when we compare with the Schechter fit to the observed luminosity function. The corresponding lower limit for sterile neutrinos depends on the modeling of the production mechanism; for instance m_sterile > 4 keV holds for the Shi-Fuller mechanism. We discuss the impact of observational uncertainties on the above bound on m_X. As a baseline for comparison with forthcoming observations from the HST Frontier Field, we provide predictions for the abundance of faint galaxies with M_UV=-13 for different values of m_X and of the star formation efficiency, valid up to z~4.

The triggering mechanisms for Active Galactic Nuclei (AGN) are still debated. Some of the most popular ones include galaxy interactions (IT) and disk instabilities (DI). Using an advanced semi analytic model (SAM) of galaxy formation, coupled to accurate halo occupation distribution modeling, we investigate the imprint left by each separate triggering process on the clustering strength of AGN at small and large scales. Our main results are as follows: i) DIs, irrespective of their exact implementation in the SAM, tend to fall short in triggering AGN activity in galaxies at the center of halos with $M_h>10^{13.5} h^{-1}M_{\odot}$. On the contrary, the IT scenario predicts abundance of active, central galaxies that generally agrees well with observations at every halo mass. ii) The relative number of satellite AGN in DIs at intermediate-to-low luminosities is always significantly higher than in IT models, especially in groups and clusters. The low AGN satellite fraction predicted for the IT scenario might suggest that different feeding modes could simultaneously contribute to the triggering of satellite AGN. iii) Both scenarios are quite degenerate in matching large-scale clustering measurements, suggesting that the sole average bias might not be an effective observational constraint. iv) Our analysis suggests the presence of both a mild luminosity and a more consistent redshift dependence in the AGN clustering, with AGN inhabiting progressively less massive dark matter halos as the redshift increases. We also discuss the impact of different observational selection cuts in measuring AGN clustering, including possible discrepancies between optical and X-ray surveys.

We have discovered 11 ultra-faint ($r\lesssim 22.1$) low surface brightness (LSB, central surface brightness $23\lesssim \mu_r\lesssim 26$) dwarf galaxy candidates in one deep Virgo field of just $576$ arcmin$^2$ obtained by the Large Binocular Camera (LBC) at the Large Binocular Telescope (LBT). Their association with the Virgo cluster is supported by their distinct position in the central surface brightness - total magnitude plane with respect to the background galaxies of similar total magnitude. They have typical absolute magnitudes and scale sizes, if at the distance of Virgo, in the range $-13\lesssim M_r\lesssim -9$ and $250\lesssim r_s\lesssim 850$ pc, respectively. Their colors are consistent with a gradually declining star formation history with a specific star formation rate of the order of $10^{-11}$ yr$^{-1}$, i.e. 10 times lower than that of main sequence star forming galaxies. They are older than the cluster formation age and appear regular in morphology. They represent the faintest extremes of the population of low luminosity LSB dwarfs that has been recently detected in wider surveys of the Virgo cluster. Thanks to the depth of our observations we are able to extend the Virgo luminosity function down to $M_r\sim -9.3$ (corresponding to total masses $M\sim 10^7$ M$_{\odot}$), finding an average faint-end slope $\alpha\simeq -1.4$. This relatively steep slope puts interesting constraints on the nature of the Dark Matter and in particular on warm Dark Matter (WDM) often invoked to solve the overprediction of the dwarf number density by the standard CDM scenario. We derive a lower limit on the WDM particle mass $>1.5$ keV.

We present the best sensitivity and angular resolution maps of the molecular disk and outflow of Mrk 231, as traced by CO observations obtained with IRAM/PdBI, and we analyze archival Chandra and NuSTAR observations. We constrain the physical properties of both the molecular disk and outflow, the presence of a highly-ionized ultra-fast nuclear wind, and their connection. The molecular outflow has a size of ~1 kpc, and extends in all directions around the nucleus, being more prominent along the south-west to north-east direction, suggesting a wide-angle biconical geometry. The maximum projected velocity of the outflow is nearly constant out to ~1 kpc, thus implying that the density of the outflowing material decreases from the nucleus outwards as $r^{-2}$. This suggests that either a large part of the gas leaves the flow during its expansion or that the bulk of the outflow has not yet reached out to ~1 kpc, thus implying a limit on its age of ~1 Myr. We find $\dot M_{OF}=[ 500-1000]~ M_{\odot}~yr^{-1}$ and $\dot E_{kin,OF}=[7-10]\times 10^{43}$ erg s$^{-1}$. Remarkably, our analysis of the X-ray data reveals a nuclear ultra-fast outflow (UFO) with velocity -20000 km s$^{-1}$, $\dot M_{UFO}=[0.3- 2.1] ~M_\odot yr^{-1}$, and momentum load $\dot P_{UFO}/\dot P_{rad}=[0.2-1.6]$.We find $\dot E_{kin,UFO}\sim \dot E_{kin,OF}$ as predicted for outflows undergoing an energy conserving expansion. This suggests that most of the UFO kinetic energy is transferred to mechanical energy of the kpc-scale outflow, strongly supporting that the energy released during accretion of matter onto super-massive black holes is the ultimate driver of giant massive outflows. We estimate a momentum boost $\dot P_{OF}/\dot P_{UFO}\approx [30-60]$. The ratios $\dot E_{kin, UFO}/L_{bol,AGN} =[ 1-5]\%$ and $\dot E_{kin,OF}/L_{bol,AGN} = [1-3]\%$ agree with the requirements of the most popular models of AGN feedback.

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.

Using an advanced semi analytic model (SAM) for galaxy formation, we have investigated the statistical effects of assuming two different mechanisms for triggering AGN activity on the properties of AGN host galaxies. We have considered a first accretion mode where AGN activity is triggered by disk instabilities (DI) in isolated galaxies, and a second feeding mode where such an activity is triggered by galaxy mergers and fly-by events (interactions, IT). We obtained the following results:i) for hosts with $M_* \lesssim 10^{11} M_{\bigodot}$, both DI and IT modes are able to account for the observed AGN hosts stellar mass function; for more massive hosts, the DI scenario predicts a lower space density than the IT model, lying below the observational estimates for z>0.8.ii) The analysis of the color-magnitude diagram (CMD) of AGN hosts for redshift z < 1.5 can provide a good observational test to effectively discriminate between the DI and IT mode, since DIs are expected to yield AGN host galaxy colors skewed towards bluer colors, while in the IT scenario the majority of hosts are expected to reside in the red sequence.iii) While both IT and DI scenarios can account for AGN triggered in main sequence or starburst galaxies, DIs fail in triggering AGN activity in passive galaxies.iv) The two modes are characterized by a different duration of the AGN phase, with DIs lasting even on time scales $\sim $ Gyr, much longer with respect to the IT scenario.v) The scatter of the $SFR-L_{bol}$ relation could represent another crucial diagnostics to discriminate between the two triggering modes, since the DI scenario predicts an appreciably lower scatter of the relation than the IT scenario. vi) Disk instabilities are not able to account for the observed fraction of AGN in groups for z < 1 and clusters for z < 0.7, while the IT scenario provides a good match to observational data.

The process that quenched star formation in galaxies at intermediate and high redshift is still the subject of considerable debate. One way to investigate this puzzling issue is to study the number density of quiescent galaxies at z~2, and its dependence on mass. Here we present the results of a new study based on very deep Ks-band imaging (with the HAWK-I instrument on the VLT) of two HST CANDELS fields (the UKIDSS Ultra-deep survey (UDS) field and GOODS-South). The new HAWK-I data (taken as part of the HUGS VLT Large Program) reach detection limits of Ks>26 (AB mag). We select a sample of passively-evolving galaxies in the redshift range 1.4<z<2.5. Thanks to the depth and large area coverage of our imaging, we have been able to extend the selection of quiescent galaxies a magnitude fainter than previous analyses. Through extensive simulations we demonstrate, for the first time, that the observed turn-over in the number of quiescent galaxies at K>22 is real. This has enabled us to establish unambiguously that the number counts of quiescent galaxies at z~2 flatten and slightly decline at magnitudes fainter than Ks~22(AB mag.). We show that this trend corresponds to a stellar mass threshold $M_*10^{10.8}\,{\rm M_{\odot}}$ below which the mechanism that halts the star formation in high-redshift galaxies seems to be inefficient. Finally we compare the observed pBzK number counts with those of quiescent galaxies extracted from four different semi-analytic models. We find that none of the models provides a statistically acceptable description of the number density of quiescent galaxies at these redshifts. We conclude that the mass function of quiescent galaxies as a function of redshift continues to present a key and demanding challenge for proposed models of galaxy formation and evolution.

Using a semi analytic model for galaxy formation we investigate the effects of Black Hole accretion triggered by disk instabilities (DI) in isolated galaxies on the evolution of AGN. Specifically, we took on, developed and expanded the Hopkins & Quataert (2011) model for the mass inflow following disk perturbations, and compare the corresponding evolution of the AGN population with that arising in a scenario where galaxy interactions trigger AGN (IT mode). We extended and developed the DI model by including different disk surface density profiles, to study the maximal contribution of DI to the evolution of the AGN population. We obtained the following results: i) for luminosities corresponding to $M_{1450}\gtrsim -26$ the DI mode can provide the BH accretion needed to match the observed AGN luminosity functions up to $z \approx 4.5$; in such a luminosity range and redshift, it can compete with the IT scenario as the main driver of cosmological evolution of AGN; ii) The DI scenario cannot provide the observed abundance of high-luminosity QSO with $M_{1450}\lesssim -26$ AGN, as well as the abundance of high-redhshift $z \approx 4.5$ QSOs with $M_{1450}\lesssim -24$, while the IT scenario provides an acceptable match up to $z \approx 6$, as found in our earliest works; iii) The dispersion of the distributions of Eddington ratio for low- and intermediate-luminosity AGN (bolometric $L_{AGN}$ = $10^{43}$ - $10^{45}$ erg/s) is predicted to be much smaller in the DI scenario compared to the IT mode; iv) The above conclusions are robust with respect to the explored variants of the Hopkins & Quataert (2011) model. We discuss the physical origin of our findings, and how it is possible to pin down the dominant fueling mechanism in the low-intermediate luminosity range $M_{1450}\gtrsim -26$ where both the DI and the IT modes are viable candidates as drivers for the AGN evolution.

We present near-infrared spectra, obtained with SINFONI and XShooter observations at ESO VLT, of nine dusty, red QSOs at 1.2<z<2.6. The sources are hard X-ray detected, characterized by cold absorption (N_H>10^{21} - 10^{22} cm^{-2}) and show a broad Ha component in the NIR spectra. We complement this sample with twelve additional sources taken from the literature with similar properties resulting in a total sample of 21 X-ray obscured, intermediate type (1.8-1.9), dusty reddened QSOs. From the broad Ha line we have computed the BH masses through the virial formula and derived Eddington ratios. Moreover, from optical/IR multi-component SED fitting we have derived the stellar mass of their host galaxies and their SFRs. We find that most of the sources in our sample are hosted in starburst and main sequence star-forming galaxies with Eddington ratios lambda>0.1. We find a strong trend with the BH mass i.e. less massive objects are scattered below and above the local relation while the most massive ones are mainly located above it. We also studied the evolution of these sources on the M_BH-M_star plane compared to a sample of optically blue type--1 QSOs and we find that obscured red QSOs show a ratio of M_BH to M_star that increases with redshift which is consistent with or slightly lower than what has been found for blue QSOs. These sources may represent the blow-out phase at the end of the rapid BH growth and immediately preceding the classical blue QSOs typically sampled in optical surveys. They in fact show evidence of outflows in the ionized gas component, but their BH has already fully formed.

We investigate various galaxy population properties of the massive X-ray luminous galaxy cluster XDCP J0044.0-2033 at z=1.58, which constitutes the most extreme matter density peak at this redshift currently known. We analyze deep VLT/HAWK-I NIR data in the J- and Ks-bands, complemented by Subaru imaging in i and V, Spitzer observations at 4.5 micron, and new spectroscopic observations with VLT/FORS2. We detect a cluster-associated excess population of about 90 galaxies, which follows a centrally peaked, compact NFW galaxy surface density profile with a concentration of c200~10. Based on the Spitzer 4.5 micron imaging data, we measure a stellar mass fraction of fstar,500=(3.3+-1.4)% consistent with local values. The total J- and Ks-band galaxy luminosity functions of the core region yield characteristic magnitudes J* and Ks* consistent with expectations from simple z_f=3 burst models. However, a detailed look at the morphologies and color distributions of the spectroscopically confirmed members reveals that the most massive galaxies are undergoing a very active mass assembly epoch through merging processes. Consequently, the bright end of the cluster red-sequence is not in place, while at intermediate magnitudes [Ks*,Ks*+1.6] a red-locus population is present, which is then sharply truncated at magnitudes fainter than Ks*+1.6. The dominant cluster core population comprises post-quenched galaxies transitioning towards the red-sequence at intermediate magnitudes, while additionally a significant blue cloud population of faint star-forming galaxies is present even in the densest central regions. Our observations lend support to the scenario in which the dominant effect of the dense z~1.6 cluster environment is an accelerated mass assembly timescale through merging activity that is responsible for driving core galaxies across the mass quenching threshold of log(Mstar/Msun)~10.4.

By means of a semi-analytic model of galaxy formation, we show how the local observed relation between age and galactic stellar mass is affected by assuming a DM power spectrum with a small-scale cutoff. We compare results obtained by means of both a Lambda-cold dark matter (LambdaCDM) and a Lambda-warm dark matter (LambdaWDM) power spectrum - suppressed with respect to the LambdaCDM at scales below ~ 1 Mpc. We show that, within a LWDM cosmology with a thermal relic particle mass of 0.75 keV, both the mass-weighted and the luminosity-weighted age-mass relations are steeper than those obtained within a LambdaCDM universe, in better agreement with the observed relations. Moreover, both the observed differential and cumulative age distributions are better reproduced within a LambdaWDM cosmology. In such a scenario, star formation appears globally delayed with respect to the LambdaCDM, in particular in low-mass galaxies. The difficulty of obtaining a full agreement between model results and observations is to be ascribed to our present poor understanding of baryonic physics.

We have evaluated the diffuse intracluster light (ICL) in the central core of the galaxy cluster CL0024+17 at $z\sim 0.4$ observed with the prime focus camera (LBC) at LBT. The measure required an accurate removal of the galaxies light within $\sim 200$ kpc from the center. The residual background intensity has then been integrated in circular apertures to derive the average ICL intensity profile. The latter shows an approximate exponential decline as expected from theoretical cold dark matter models. The radial profile of the ICL over the galaxies intensity ratio (ICL fraction) is increasing with decreasing radius but near the cluster center it starts to bend and then decreases where the overlap of the halos of the brightest cluster galaxies becomes dominant. Theoretical expectations in a simplified CDM scenario show that the ICL fraction profile can be estimated from the stripped over galaxy stellar mass ratio in the cluster. It is possible to show that the latter quantity is almost independent of the properties of the individual host galaxies but mainly depends on the average cluster properties. The predicted ICL fraction profile is thus very sensitive to the assumed CDM profile, total mass and concentration parameter of the cluster. Adopting values very similar to those derived from the most recent lensing analysis in CL0024+17 we find a good agreement with the observed ICL fraction profile. The galaxy counts in the cluster core have then been compared with that derived from composite cluster samples in larger volumes, up to the clusters virial radius. The galaxy counts in the CL0024+17 core appear flatter and the amount of bending respect to the average cluster galaxy counts imply a loss of total emissivity in broad agreement with the measured ICL fraction.

We use deep Herschel PACS and SPIRE observations in GOODSS, GOODSN and COSMOS to estimate the average dust mass (Mdust) of galaxies on a redshift-stellar mass (Mstar)-SFR grid. We study the scaling relations between Mdust, Mstar and SFR at z<=2.5. No clear evolution of Mdust is observed at fixed SFR and Mstar. We find a tight correlation between SFR and Mdust, likely a consequence of the Schmidt-Kennicutt (S-K) law. The Mstar-Mdust correlation observed by previous works flattens or sometimes disappears when fixing the SFR. Most of it likely derives from the combination of the Mdust-SFR and Mstar-SFR correlations. We then investigate the gas content as inferred by converting Mdust by assuming that the dust/gas ratio scales linearly with the gas metallicity. All galaxies in the sample follow, within uncertainties, the same SFR-Mgas relation (integrated S-K law), which broadly agrees with CO-based results for the bulk of the population, despite the completely different approaches. The majority of galaxies at z~2 form stars with an efficiency (SFE=SFR/Mgas) ~5 times higher than at z~0. It is not clear what fraction of such variation is an intrinsic redshift evolution and what fraction arises from selection effects. The gas fraction (fgas) decreases with Mstar and increases with SFR, and does not evolve with z at fixed Mstar and SFR. We explain these trends by introducing a universal relation between fgas, Mstar and SFR, non-evolving out to z~2.5. Galaxies move across this relation as their gas content evolves in time. We use the 3D fundamental fgas-Mstar-SFR relation and the redshift evolution of the Main Sequence to estimate the evolution of fgas in the average population of galaxies as a function of z and Mstar, and we find evidence a downsizing scenario.

We estimate the fraction of AGNs hosted in starburst galaxies (f_bursty) as a function of the AGN luminosity predicted under the assumption that starburst events and AGN activity are triggered by galaxy interactions during their merging histories. The latter are described through Monte Carlo realizations, and are connected to star formation and BH accretion using a semi-analytic model of galaxy formation in a cosmological framework. The predicted fraction f_bursty increases steeply with AGN luminosity from <0.2 at L_X < 10^44 erg/s to >0.9 at L_X > 10^45 erg/s over a wide redshift interval from z=0 to z=6. We compare the model predictions with new measurements of f_bursty from a sample of X-ray selected AGNs in the XMM-COSMOS field at 0.3< z< 2, and from a sample of QSOs (L_X > 10^45 erg/s) in the redshift range 2< z< 6.5. We find preliminary indications that under conservative assumptions half of the QSO hosts are starburst galaxies. This result provide motivation for future systematic studies of the stellar properties of high luminosity AGN hosts in order to constrain AGN triggering mechanisms.

The standard cosmological model based on cold dark matter (CDM) predicts a large number of subhalos for each galaxy-size halo. It is well known that matching the subhalos to the observed properties of luminous satellites of galaxies in the local universe poses a significant challenge to our understanding of the astrophysics of galaxy formation. We show that the cosmic evolution and host mass dependence of the luminosity function of satellites provides a powerful new diagnostic to disentangle astrophysical effects from variations in the underlying dark matter mass function. We illustrate this by comparing the results of recent observations of satellites out to $z=0.8$ based on Hubble Space Telescope images with the predictions of three different sets of state-of-the art semi-analytic models with underlying CDM power spectra and one semi-analytic model with an underlying Warm Dark Matter (WDM) power spectrum. We find that even though CDM models provide a reasonable fit to the local luminosity function of satellites around galaxies comparable or slightly larger than the Milky Way, they do not reproduce the data as well for different redshift and host galaxy stellar mass. This tension indicates that further improvements are likely to be needed in the description of star formation if the models are to be reconciled with the data. The WDM model matches the observed mass dependence and redshift evolution of satellite galaxies more closely than any of the CDM models, indicating that a modification of the underlying power spectrum may offer an alternative solution to this tension. We conclude by presenting predictions for the color magnitude relation of satellite galaxies to demonstrate how future observations will be able to further distinguish between these models and help constrain baryonic and non-baryonic physics.

Recent measurements of the abundance of AGN with low-luminosities (L_X< 10^44 erg/s in the 2-10 keV energy band) at high redshifts z>4 provide a serious challenge for Cold Dark Matter (CDM) models based on interaction-driven fueling of AGN. Using a semi-analytic model of galaxy formation we investigate how such observations fit in a Warm Dark Matter (WDM) scenario of galaxy formation, and compare the results with those obtained in the standard CDM scenario with different efficiencies for the stellar feedback. Taking on our previous exploration of galaxy formation in WDM cosmology, we assume as a reference case a spectrum which is suppressed - compared to the standard CDM case - below a cut-off scale ~ 0.2$ Mpc corresponding (for thermal relic WDM particles) to a mass m_X=0.75 keV. We run our fiducial semi-analytic model with such a WDM spectrum to derive AGN luminosity functions from z~6 to the present over a wide range of luminosities (10^43< L_X/erg/s < 10^46 in the 2-10 keV X-ray band), to compare with recent observations and with the results in the CDM case. When compared with the standard CDM case, the luminosity distributions we obtain assuming a WDM spectrum are characterized by a similar behaviour at low redshift, and by a flatter slope at faint magnitudes for z>3, which provide an excellent fit to present observations. We discuss how such a result compares with CDM models with maximized feedback efficiency, and how future deep AGN surveys will allow for a better discrimination between feedback and cosmological effects on the evolution of AGN in interaction-driven models for AGN fueling.

Determining the AGN content in structures of different mass/velocity dispersion and comparing them to higher mass/lower redshift analogs is important to understand how the AGN formation process is related to environmental properties. We use our well-tested cluster finding algorithm to identify structures in the GOODS North and South fields, exploiting the available spectroscopic redshifts and accurate photometric redshifts. We identify 9 structures in GOODS-south (presented in a previous paper) and 8 new structures in GOODS-north. We only consider structures where at least 2/3 of the members brighter than M_R=-20 have a spectroscopic redshift. For those group members that coincide with X-ray sources in the 4 and 2 Msec Chandra source catalogs respectively, we determine if the X-ray emission originates from AGN activity or it is related to the galaxies' star-formation activity. We find that the fraction of AGN with Log L_H > 42 erg/s in galaxies with M_R < -20 is on average 6.3+-1.3%, much higher than in lower redshift groups of similar mass and more than double the fraction found in massive clusters at a similarly high redshift. We then explore the spatial distribution of AGN in the structures and find that they preferentially populate the outer regions. The colors of AGN host galaxies in structures tend to be confined to the green valley, thus avoiding the blue cloud and, partially, also the red-sequence, contrary to what happens in the field. We finally compare our results to the predictions of two sets of semi analytic models to investigate the evolution of AGN and evaluate potential triggering and fueling mechanisms. The outcome of this comparison attests the importance of galaxy encounters, not necessarily leading to mergers, as an efficient AGN triggering mechanism. (abridged)

An increasing amount of observational evidence supports the notion that there are two modes of star formation: a quiescent mode in disk-like galaxies, and a starburst mode, which is generally interpreted as driven by merging. Using a semi-analytic model of galaxy formation, we derive the relative contribution to the cosmic star formation rate density of quiescently starforming and starburst galaxies, predicted under the assumption that starburst events are triggered by galaxy encounters (merging and fly-by kind) during their merging histories. We show that, within this framework, quiescently starforming galaxies dominate the cosmic star formation rate density at all redshifts. The contribution of the burst-dominated starforming galaxies increases with redshift, rising from <5% at low redshift (z<0.1) to ~20% at z>5. We estimated that the fraction of the final (z=0) galaxy stellar mass which is formed through the burst component of star formation is ~10% for 10^10 M_\odot<M_*<10^11.5 M_\odot. Starburst galaxies, selected according to their distance from the galaxy main sequence, account for ~10% of the star formation rate density in the redshift interval 1.5<z<2.5, i.e. at the cosmic peak of the star formation activity.

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.

Massive AGN-driven outflows are invoked by AGN-galaxy co-evolutionary models to suppress both star formation and black hole accretion. Massive molecular outflows have recently been revealed in some AGN hosts. However, the physical properties and structure of these AGN-driven molecular outflows are still poorly constrained. Here we present new IRAM PdBI observations of Mrk231, the closest quasar known, targeting both the CO(1-0) and CO(2-1) transitions. We detect broad wings in both transitions, tracing a massive molecular outflow with velocities up to 800 km/s. The wings are spatially resolved at high significance level (5-11 sigma), indicating that the molecular outflow extends on the kpc scale. The CO(2-1)/CO(1-0) ratio of the red broad wings is consistent with the ratio observed in the narrow core, while the blue broad wing is less excited than the core. The latter result suggests that quasar driven outflow models invoking shocks (which would predict higher gas excitation) are not appropriate to describe the bulk of the outflow in Mrk231. However, we note that within the central 700 pc the CO(2-1)/CO(1-0) ratio of the red wing is slightly, but significantly, higher than in the line core, suggesting that shocks may play a role in the central region. We also find that the average size of the outflow anticorrelates with the critical density of the transition used as a wind tracer. This indicates that, although diffuse and dense clumps coexist in the outflowing gas, dense outflowing clouds have shorter lifetime and that they evaporate into the diffuse component along the outflow or, more simply, that diffuse clouds are more efficiently accelerated to larger distances by radiation pressure.

We have evaluated the contribution of the AGN population to the ionization history of the Universe based on a semi-analytic model of galaxy formation and evolution in the CDM cosmological scenario. The model connects the growth of black holes and of the ensuing AGN activity to galaxy interactions. In the model we have included a self consistent physical description of the escape of ionizing UV photons; this is based on the blast-wave model for the AGN feedback we developed in a previous paper to explain the distribution of hydrogen column densities in AGNs of various redshifts and luminosities, due to absorption by the host galaxy gas. The model predicts UV luminosity functions for AGNs which are in good agreement with those derived from the observations especially at low and intermediate redshifts (z=3). At higher redshifts (z>5) the model tends to overestimate the data at faint luminosities. Critical biases both in the data and in the model are discussed to explain such apparent discrepancies. The predicted hydrogen photoionization rate as a function of redshift is found to be consistent with that derived from the observations. All that suggests to reconsider the role of the AGNs as the main driver of the ionization history of the Universe.

We investigate for the first time the effects of a Warm Dark Matter (WDM) power spectrum on the statistical properties of galaxies using a semi-analytic model of galaxy formation. The WDM spectrum we adopt as a reference case is suppressed - compared to the standard Cold Dark Matter (CDM) case - below a cut-off scale ~ 1 Mpc corresponding (for thermal relic WDM particles) to a mass m_X=0.75 keV. This ensures consistency with present bounds provided by the microwave background WMAP data and by the comparison of hydrodynamical N-body simulations with observed Lyman-{\alpha} forest. We run our fiducial semi-analytic model with such a WDM spectrum to derive galaxy luminosity functions (in B, UV, and K bands) and the stellar mass distributions over a wide range of cosmic epochs, to compare with recent observations and with the results in the CDM case. The predicted color distribution of galaxies in the WDM model is also checked against the data. When compared with the standard CDM case, the luminosity and stellar mass distributions we obtain assuming a WDM spectrum are characterized by: i) a flattening of the faint end slope and ii) a sharpening of the cutoff at the bright end for z \lesssim 0.8. We discuss how the former result is directly related to the smaller number of low-mass haloes collapsing in the WDM scenario, while the latter is related to the smaller number of satellite galaxies accumulating in massive haloes at low redshift, thus suppressing the accretion of small lumps on the central, massive galaxies. These results shows how a adopting a WDM power spectrum may contribute to solve two major problems of CDM galaxy formation scenarios, namely, the excess of predicted faint (low mass) galaxies at low and - most of all - high redshifts, and the excess of bright (massive) galaxies at low redshifts.

We plan to analyse dust extinction in Lyman Break Galaxies (LBGs) by introducing a new and more reliable approach to their selection and to the characterization of their distribution of UV slopes beta, using deep IR images from HST. We exploit deep WFC3 IR observations of the ERS and HUDF fields over GOODS-South, combined with HST-ACS optical data, to select z~4 LBGs through a new (B-V) vs. (V-H) colour diagram. The UV slope of the selected galaxies is robustly determined by a linear fit over their observed I, Z, Y, J magnitudes, coherently with the original definition of beta. The same fit is used to determine their rest-frame UV magnitudes M1600 through a simple interpolation. We estimate the effect of observational uncertainties with detailed simulations that we also exploit, under a parametric maximum-likelihood approach, to constrain the probability density function of UV slopes PDF(beta) as a function of rest-frame magnitude. We find 142 and 25 robust LBGs in the ERS and HUDF fields respectively, limiting our sample to S/N(H)>10 objects. Our newly defined criteria improve the selection of z~4 LBGs and allow us to exclude red interlopers at lower redshift, especially z~3-3.5 objects. We find that z~4 LBGs are characterized by blue UV slopes, suggesting a low dust extinction: all L<L* galaxies have an average UV slope <beta>\simeq-2.1, while brighter objects only are slightly redder (<beta> -1.9). We find an intrinsic dispersion ~ 0.3 for PDF(beta) at all magnitudes. The SFRD at z~4 corrected according to these estimates turns out to be lower than previously found: log(SFRD)\simeq-1.09 M_sun/yr/Mpc^3. Finally, we discuss how the UV slope of z~4 galaxies changes as a function of the dust-corrected UV magnitude (i.e. SFR) showing that most galaxies with a high SFR (> 80 M_sun/yr) are highly extincted objects. [Abridged]

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.

We compute the number density of massive Black Holes (BHs) at the centre of galaxies at z=6 in different Dynamical Dark Energy (DDE) cosmologies, and compare it with existing observational lower limits, to derive constraints on the evolution of the Dark Energy equation of state parameter w. Our approach only assumes the canonical scenario for structure formation from the collapse of overdense regions of the Dark Matter dominated primordial density field on progressively larger scales; the Black Hole accretion and merging rate have been maximized in the computation so as to obtain robust constraints on w and on its look-back time derivative w_a. Our results provide independent constraints complementary to those obtained by combining Supernovae, Cosmic Microwave Background and Baryonic Acoustic Oscillations; while the latter concern combinations of w_0 and w_a leaving the time evolution of the state parameter w_a highly unconstrained, the BH abundance mainly provide upper limits on w_a, only weakly depending on w_0. Combined with the existing constraints, our results significantly restrict the allowed region in DDE parameter space, ruling out DDE models not providing cosmic time and fast growth factor large enough to allow for the building up of the observed abundance of BHs; in particular, models with -1.2 \leq w_0 \leq -1 and positive redshift evolution w_a > 0.8 - completely consistent with previous constraints - are strongly disfavoured by our independent constraints from BH abundance. Such range of parameters corresponds to "Quintom" DDE models, with w crossing -1 starting from larger values.

We present detection and analysis of faint X-ray sources in the Chandra deep field south (CDFS) using the 4 Msec Chandra observation and adopting a new detection algorithm, based on a targeted search at the position of known high-z galaxies. This optimized technique results in the identification of 54 z>3 AGNs, 29 of which are new detections. Applying stringent completeness criteria, we derive AGN luminosity functions in the redshift bins 3-4, 4-5 and >5.8 and for 42.75<log L(2-10 keV)<44.5. We join this data with the luminous AGN luminosity functions from optical surveys and find that the evolution of the high-z, wide luminosity range luminosity function can be best modeled by pure luminosity evolution with L* decreasing from 6x10^44 ergs/s at z=3 to L*=2x10^44 ergs/s at z=6. We compare the high-z luminosity function with the prediction of theoretical models using galaxy interactions as AGN triggering mechanism. We find that these models are broadly able to reproduce the high-z AGN luminosity functions. A better agreement is found assuming a minimum Dark Matter halo mass for black hole formation and growth. We compare our AGN luminosity functions with galaxy mass functions to derive high-z AGN duty cycle using observed Eddington ratio distributions to derive black hole masses. We find that the duty cycle increases with galaxy stellar mass and redshift by a factor 10-30 from z=0.25 to z=4-5. We also report on the detection of a large fraction of highly obscured, Compton thick AGN at z>3 (18+17-10%). Their optical counterparts are not strongly reddened and we thus conclude that the size of the X-ray absorber is likely smaller than the dust sublimation radius. We finally report the discovery of a highly star-forming galaxy at z=3.47. If confirmed, this would be one of the farthest objects in which stellar sources are detected in X-rays.

[Abridged] We use the large public spectroscopic database available in the GOODS-South field to estimate the dynamical mass and the virialization status of cluster ClG 0332-2747 at z=0.734. Cluster members selected from their photometric redshift are used with spectroscopic ones to analyse the galaxy population of the cluster. In the newly released Chandra 4Ms observations we detect a faint extended X-ray emission associated to the cluster. Finally, we compare the optical and X-ray properties of ClG 0332-2747 with the predictions of a well tested semianalytic model. We estimate the velocity dispersion and the virial mass considering all 44 spectroscopic members, or 20 red-sequence members only. We obtain sigma_v=634 +/- 105 Km/s, M_200=3.07 ^{+1.57}_{-1.16}~10^{14} M_sun in the former case, and slightly lower values in the latter case. The cluster appears to have reached the virial equilibrium: it shows a perfectly Gaussian velocity distribution and no evidence for substructures. ClG 0332-2747 contains a high fraction of bright red galaxies, and is dominated by a very massive (1.1 x 10^{12} M_sun) old brightest cluster galaxy (BCG), suggesting that it formed at an early epoch. We detect a faint extended X-ray source centered on the BCG, with a total X-ray luminosity of L_X ~ 2 x 10^{42} erg s^-1 (0.1-2.4 keV). This L_X is lower by a factor of ~10-20 than expected according to the M-L_X relation. We provide a possible explanation of this discrepancy as due to the effects of AGN feedback on the ICM: the semianalytic model reproduces the M-L_X relation measured from "X-ray bright" clusters, and it predicts a high scatter at low masses due to heating and expulsion of the cluster gas. Interestingly, in the model clusters with an evolved galaxy population like ClG 0332-2747 present the largest scatter in X-ray luminosity.

The most massive elliptical galaxies apparently formed the fastest, because the ratio of alpha elements (such as oxygen) to iron is the smallest. In fact, iron is mainly produced from type Ia supernovae on a timescale of ~ 0.1-1 billion years, while the alpha elements come from massive stars on timescales of a few tens of million years (Matteucci 1994). Reproducing such a alpha/Fe correlation has long been a severe problem for cosmological theories of galaxy formation, which envisage massive galaxies to assemble gradually from smaller progenitors, and to be characterized by a star formation history too much extended towards late cosmic times. While it has recently become clear that feedback from Active Galactic Nuclei (AGN) activity play a role in the late quenching of star formation (e.g. Cattaneo et al. 2009), and that early star formation history in the galaxy progenitors affect the alpha/Fe ratio (Calura & Menci 2009), major mergers alone cannot enhance the star formation in the high-redshift progenitors to the levels required to match the steepness of the observed alpha/Fe correlation (Spolaor et al. 2010). Here we report that the inclusion of the effects of fly-by 'harassments', that trigger lower level starbursts, combined with the AGN quenching of the starburst activity, considerably enhances the capability to account for the observed alpha/Fe ratio in ellipticals within cosmological galaxy formation models . The critical difference between the earlier work and the present result is the effect of starbursts driven by fly-by encounters that would have been very common amongst the high-redshift progenitors of massive galaxies and which would have boosted star formation in the first 2 billion years after the Big Bang, combined with quenching of the burst activity within the first 3-4 Gyr.

(Abridged) We present here the second half of an ESO Large Programme, which exploits the unique combination of area and sensitivity provided in the near-IR by the camera Hawk-I at the VLT. We have obtained - 30 observing hours with Hawk-I in the Y-band of two high galactic latitude fields. We combined the Y-band data with deep J and K Hawk-I observations, and with FORS1/FORS2 U, B, V, R, I, and Z observations to select z-drop galaxies having Z - Y > 1, no optical detection and flat Y - J and Y - K colour terms. We detect 8 high-quality candidates in the magnitude range Y = 25.5 - 26.5 that we add to the z-drop candidates selected in two Hawk-I pointings over the GOODS-South field. We use this full sample of 15 objects found in -161 arcmin^2 of our survey to constrain the average physical properties and the evolution of the number density of z ~ 7 LBGs. A stacking analysis yields a best-fit SED with photometric redshift z= 6.85 +0.20 -0.15 and an E(B-V)=0.05 +0.15 -0.05. We compute a binned estimate of the z ~ 7 LF and explore the effects of photometric scatter and model uncertainties on the statistical constraints. After accounting for the expected incompleteness through MonteCarlo simulations, we strengthen our previous finding that a Schechter luminosity function constant from z=6 to z=7 is ruled out at a >99% confidence level, even including the effects of cosmic variance. For galaxies brighter than M_1500= -19.0, we derive a luminosity density rho_UV = 1.5^{+2.1}{-0.8} x 10^25 erg/s/Hz/Mpc^3, implying a decrease by a factor 3.5 from z=6 to z=6.8. We find that, under standard assumptions, the emission rate of ionizing photons coming from UV bright galaxies is lower by at least a factor of two than the value required for reionization. Finally, we exploit deep Hawk-I J and K band observations to derive an upper limit on the number density of M1500<~ -22.0 LBGs at z-8 (Y-dropouts).

In the standard scenario for galaxy evolution young star-forming galaxies transform into red bulge-dominated spheroids, where star formation has been quenched. To explain such a transformation, a strong negative feedback generated by accretion onto a super-massive black hole is often invoked. The depletion of gas resulting from quasar-driven outflows should eventually stop star-formation in the host galaxy and lead the black hole to "suicide" for starvation. Direct observational evidence for a major quasar feedback onto the host galaxy is still missing, since outflows previously observed in quasars are generally associated with the ionized component of the gas, which only accounts for a minor fraction of the total gas content, and typically occurring in the central regions. We used the IRAM PdB Interferometer to observe CO(1-0) in Mrk 231, the closest QSO known. Thanks to the wide band we detect broad wings of the CO line, with velocities up to 750 km/s and spatially resolved on the kpc scale. Such broad CO wings trace a giant molecular outflow of about 700 MSun/yr, far larger than the ongoing star-formation rate (~200 MSun/yr) observed in the host galaxy. This wind will totally expel the cold gas reservoir in Mrk 231 in about 1e7 yrs, therefore halting the star-formation activity on the same timescale. The inferred kinetic energy in the molecular outflow is ~1.2e44 erg/s, corresponding to a few percent of the AGN bolometric luminosity, very close to the fraction expected by models ascribing quasar feedback to highly supersonic shocks generated by radiatively accelerated nuclear winds. The direct observational evidence for quasar feedback reported here provides solid support to the scenarios ascribing the observed properties of local massive galaxies to quasar-induced large scale winds.

We derive the growth of SMBHs relative to the stellar content of their host galaxy predicted under the assumption of BH accretion triggered by galaxy encounters occurring during their merging histories. We show that, within this framework, the ratio Gamma=(M_BH/M_*)(z)/(M_BH/M_*)(z=0) between the Black Hole mass and the galactic stellar mass (normalized to the local value) depends on both BH mass and redshift. While the average value and the spread of Gamma(z) increase with z, such an effect is larger for massive BHs, reaching values Gamma=5 for massive Black Holes (M>10^9 M_{\odot}) at z>4, in agreement with recent observations of high-redshift QSOs; this is due to the effectiveness of interactions in triggering BH accretion in high-density environments at high redshifts. To test such a model against observations, we worked out specific predictions for sub-samples of the simulated galaxies corresponding to the different observational samples for which measurements of Gamma have been obtained. We found that for Broad Line AGNs at 1<z<2 values of Gamma=2 are expected, with a mild trend toward larger value for increasing BH mass. Instead, when we select from our Monte Carlo simulations only extremely gas rich, rapidly star forming galaxies at 2<z<3, we find low values 0.3<Gamma<1.5, consistent with recent observational findings on samples of sub-mm galaxies; in the framework of our model, these objects end up at z=0 in low-to-intermediate mass BHs (M<10^9 M_{\odot}), and they do not represent typical paths leading to local massive galaxies. The latter have formed preferentially through paths passing above the local M_*-M_BH relation. We discuss how the global picture emerging from the model is consistent with a downsizing scenario, where massive BHs accrete a larger fraction of their final mass at high redshifts z>4.

We perform a deep search for galaxies in the redshift range 6.5<z<7.5, to measure the evolution of the number density of luminous galaxies in this redshift range and derive useful constraints on the evolution of their Luminosity Function. We present here the first results of an ESO Large Program, that exploits the unique combination of area and sensitivity provided in the near-IR by the camera Hawk-I at the VLT. We have obtained two Hawk-I pointings on the GOODS South field for a total of 32 observing hours, covering ~90 arcmin2. The images reach Y=26.7 mags for the two fields. We have used public ACS images in the z band to select z-dropout galaxies with the colour criteria Z-Y>1, Y-J<1.5 and Y-K<2. The other public data in the UBVRIJHK bands are used to reject possible low redshift interlopers. The output has been compared with extensive Monte Carlo simulations to quantify the observational effects of our selection criteria as well as the effects of photometric errors. We detect 7 high quality candidates in the magnitude range Y=25.5-26.7. This interval samples the critical range for M* at z>6 (M_1500 ~- 19.5 to -21.5). After accounting for the expected incompleteness, we rule out at a 99% confidence level a Luminosity Function constant from z=6 to z=7, even including the effects of cosmic variance. For galaxies brighter than M_1500=-19.0 we derive a luminosity density rho_UV=1.5^{+2.0}_{-0.9} 10^25 erg/s/Hz/Mpc3, implying a decrease by a factor 3.5 from z=6 to z~6.8. On the basis of our findings, we make predictions for the surface densities expected in future surveys surveys, based on ULTRA-VISTA, HST-WFC3 or JWST-NIRCam, evaluating the best observational strategy to maximise their impact.

We investigate the chemical properties of local galaxies within a cosmological framework in the hierarchical picture of galaxy formation. To this aim, we use a hierarchical semi-analytic model which includes the contribution from (i) low and intermediate mass stars, (ii) type Ia Supernovae (SNe) and (iii) massive stars. - Abridged - We compare our predictions with available observations in the Milky Way (MW), in local dwarf galaxies and in local ellipticals. For Milky-Way-like galaxies, we can successfully reproduce the [O-Fe] vs [Fe/H] relation observed in disc stars and the stellar metallicity distribution (SMD). For dwarf galaxies, the stellar metallicity vs mass relation is reproduced by assuming that a substantial fraction of the heavy elements is lost through metal-enhanced outflows and a type Ia SN realization probability lower than the one of MW-like galaxies. - Abridged - In ellipticals, the observations indicate higher [alpha/Fe] values in larger galaxies. - Abridged - Our results computed with a standard Salpeter initial mass function (IMF) indicate a flat [alpha/Fe]-sigma relation. However, we suggest a possible solution to this problem and show how, by assuming a star formation-dependent IMF with a slope x=1.35 in systems with star formation rates < 100 M_sun/yr and slightly flatter (i.e. with x=1) in object with stronger star formation, the observed correlation between [alpha/Fe] and sigma can be accounted for on a large velocity dispersion range. Fundamental roles are played also by interaction-triggered starbursts and AGN.

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.

Aims: The aim of this work is to collect a complete, mass--selected sample of galaxies with very low specific star formation rate, for a comparison with the prediction of recent theoretical models. Method: We use the 24/K flux ratio, complemented by the SED fitting to the full 0.35-8.0 mum spectral distribution, to select quiescent galaxies from z~0.4 to z~4 in the GOODS--MUSIC sample. Our observational selection can be translated into thresholds on the specific star formation rate SFR/M_*, that can be used to compare with the theoretical predictions. Results: We find that, in the framework of the well known global decline of the quiescent fraction with redshift, a non-negligible fraction ~15-20% of massive galaxies with very low specific star formation rate exists up to z~4, including a tail of "Red&Dead" galaxies with SFR/M_*<10^{-11}/yr. Recent theoretical models vary to a large extent in the prediction of the fraction of galaxies with very low specific star formation rates, and are unable to provide a global match to our data.

By means of chemical evolution models of different morphological types, we study the mass-metallicity (MZ) relation and its evolution with redshift. Our aim is to understand the role of galaxies of different morphological types in the MZ relation at various redshift. One major result is that at high redshift, the majority of the galaxies falling on the MZ plot are apparently proto-ellipticals. Finally, we show some preliminary results of a study of the MZ relation in a framework of hierarchical galaxy formation.

The aim of this paper is to identify and study the properties and galactic content of groups and clusters in the GOODS-South field up to z\sim2.5, and to analyse the physical properties of galaxies as a continuous function of environmental density up to high redshift. We use the deep (z850\sim26), multi-wavelength GOODS-MUSIC catalogue, which has a 15% of spectroscopic redshifts and accurate photometric redshifts for the remaining fraction. On these data, we apply a (2+1)D algorithm, previously developed by our group, that provides an adaptive estimate of the 3D density field. We support our analysis with simulations to evaluate the purity and the completeness of the cluster catalogue produced by our algorithm. We find several high density peaks embedded in larger structures in the redshift range 0.4-2.5. From the analysis of their physical properties (mass profile, M200, \sigmav, LX, U-B vs. B diagram), we derive that most of them are groups of galaxies, while two are poor clusters with masses of few times 10^14 Mo. For these two clusters we find, from the Chandra 2Ms data, an X-ray emission significantly lower than expected from their optical properties, suggesting that the two clusters are either not virialised or gas poor. We also analyse the dependance on environment of galaxy colours, luminosities, stellar masses, ages and star formations. We find that galaxies in high density regions are, on average, more luminous and massive than field galaxies up to z\sim 2. The fraction of red galaxies increases with luminosity and with density up to z\sim 1.2. At higher z this dependance on density disappears. The variation of galaxy properties as a function of redshift and density suggests that a significant change occurs at z\sim 1.5-2. (abridged)

We discuss the uncertainties and the systematic effects that exist in the estimates of the stellar masses of high redshift galaxies, using broad band photometry, and how they affect the deduced galaxy stellar mass function. We use at this purpose the latest version of the GOODS-MUSIC catalog. In particular, we discuss the impact of different synthetic models, of the assumed initial mass function and of the selection band. Using Charlot & Bruzual 2007 and Maraston 2005 models we find masses lower than those obtained from Bruzual & Charlot 2003 models. In addition, we find a slight trend as a function of the mass itself comparing these two mass determinations with that from Bruzual & Charlot 2003 models. As consequence, the derived galaxy stellar mass functions show diverse shapes, and their slope depends on the assumed models. Despite these differences, the overall results and scenario remains unchanged. The masses obtained with the assumption of the Chabrier initial mass function are in average 0.24 dex lower than those from the Salpeter assumption, at all redshifts, causing a shift of galaxy stellar mass function of the same amount. Finally, using a 4.5 um-selected sample instead of a Ks-selected one, we add a new population of highly absorbed, dusty galaxies at z\simeq 2-3 of relatively low masses, yielding stronger constraints on the slope of the galaxy stellar mass function at lower masses.

(abridged) We take advantage of the deep Chandra and Spitzer coverage of a large area (more than 10 times the area covered by the Chandra deep fields, CDFs in the COSMOS field, to extend the search of highly obscured, Compton-thick active nuclei to higher luminosity. These sources have low surface density and large samples can be provided only through large area surveys, like the COSMOS survey. We analyze the X-ray properties of COSMOS MIPS sources with 24$\mu$m fluxes higher than 550$\mu$Jy. For the MIPS sources not directly detected in the Chandra images we produce stacked images in soft and hard X-rays bands. To estimate the fraction of Compton-thick AGN in the MIPS source population we compare the observed stacked count rates and hardness ratios to those predicted by detailed Monte Carlo simulations including both obscured AGN and star-forming galaxies. The density of lower luminosity Compton-thick AGN (logL(2-10keV)=43.5-44) at z=0.7--1.2 is $(3.7\pm1.1) \times10^{-5}$ Mpc$^{-3}$, corresponding to $\sim67%$ of that of X-ray selected AGN. The comparison between the fraction of infrared selected, Compton thick AGN to the X-ray selected, unobscured and moderately obscured AGN at high and low luminosity suggests that Compton-thick AGN follow a luminosity dependence similar to that discovered for Compton-thin AGN, becoming relatively rarer at high luminosities. We estimate that the fraction of AGN (unobscured, moderately obscured and Compton thick) to the total MIPS source population is $49\pm10%$, a value significantly higher than that previously estimated at similar 24$\mu$m fluxes. We discuss how our findings can constrain AGN feedback models.

We compare the results from a semi-analytic model of galaxy formation with spectro-photometric observations of distant galaxy clusters observed in the range 0.8< z< 1.3. We investigate the properties of their red sequence (RS) galaxies and compare them with those of the field at the same redshift. In our model we find that i) a well-defined, narrow RS is obtained already by z= 1.2; this is found to be more populated than the field RS, analogously to what observed and predicted at z=0; ii) the predicted U-V rest-frame colors and scatter of the cluster RS at z=1.2 have average values of 1 and 0.15 respectively, with a cluster-to-cluster variance of 0.2 and 0.06, respectively. The scatter of the RS of cluster galaxies is around 5 times smaller than the corresponding field value; iii) when the RS galaxies are considered, the mass growth histories of field and cluster galaxies at z=1.2 are similar, with 90 % of the stellar mass of RS galaxies at z=1.2 already formed at cosmic times t=2.5 Gyr, and 50 % at t=1 Gyr; v) the predicted distribution of stellar ages of RS galaxies at z=1.2 peaks at 3.7 Gyr for both cluster and field populations; however, for the latter the distribution is significantly skewed toward lower ages. When compared with observations, the above findings show an overall consistency, although the average value 0.07 of the observed cluster RS scatter (U-V colors) at z=1.2 is smaller than the corresponding model central value. We discuss the physical origin and the significance of the above results in the framework of cosmological galaxy formation.

We compute the effect of the galactic absorption on AGN emission in a cosmological context by including a physical model for AGN feeding and feedback in a semi-analytic model of galaxy formation. This is based on galaxy interactions as triggers for AGN accretion, and on expanding blast waves as a mechanism to propagate outwards the AGN energy injected into the interstellar medium at the center of galaxies. We first test our model against the observed number density of AGNs with different intrinsic luminosity as a function of redshift. The model yields a ''downsizing'' behavior in close agreement with the observed one for z<2. At higher redshifts, the model predicts an overall abundance of AGNs (including Compton-thick sources) larger than the observed Compton-thin sources by a factor around 2 for z>2 and L_X < 10^{44} erg/s. Thus, we expect that at such luminosities and redshifts about 1/2 of the total AGN population is contributed by Compton-thick sources. We then investigate the dependence of the absorbing column density N_H associated to cold galactic gas (and responsible for the Compton-thin component of the overall obscuration) on the AGN luminosity and redshift. We find that the absorbed fraction of AGNs with N_H>10^{22} cm^{-2} decreases with luminosity for z<1; in addition, the total (integrated over luminosity) absorbed fraction increases with redshift up to z around 2, and saturates to the value around 0.8 at higher redshifts. Finally, we predict the luminosity dependence of the absorbed fraction of AGNs with L_X< 3 10^{44} erg/s to weaken with increasing redshift. We compare our results with recent observations, and discuss their implications in the context of cosmological models of galaxy formation.

We derive stellar masses, ages and star formation histories of massive early-type galaxies in the z=1.237 RDCS1252.9-2927 cluster and compare them with those measured in a similarly mass-selected sample of field contemporaries drawn from the GOODS South Field. Robust estimates of these parameters are obtained by comparing a large grid of composite stellar population models with 8-9 band photometry in the rest-frame NUV, optical and IR, thus sampling the entire relevant domain of emission of the different stellar populations. Additionally, we present new, deep $U$-band photometry of both fields, giving access to the critical FUV rest-frame, in order to constrain empirically the dependence on the environment of the most recent star formation processes. We find that early-type galaxies, both in the cluster and in the field, show analogous optical morphologies, follow comparable mass vs. size relation, have congruent average surface stellar mass densities and lie on the same Kormendy relation. We also that a fraction of early-type galaxies in the field employ longer timescales, $\tau$, to assemble their mass than their cluster contemporaries. Hence we conclude that, while the formation epoch of early-type only depends on their mass, the environment does regulate the timescales of their star formation histories. Our deep $U$-band imaging strongly supports this conclusions. It shows that cluster galaxies are at least 0.5 mag fainter than their field contemporaries of similar mass and optical-to-infrared colors, implying that the last episode of star formation must have happened more recently in the field than in the cluster.

We study the evolution of the galaxy population up to z\sim3 as a function of its colour properties. In particular, luminosity functions and luminosity densities have been derived as a function of redshift for the blue/late and red/early populations. We use data from the GOODS-MUSIC catalogue which have typical magnitude limits z<26 and Ks<23.5 for most of the sample. About 8% of the galaxies have spectroscopic redshifts; the remaining have well calibrated photometric redshifts derived from the extremely wide multi-wavelength coverage in 14 bands (from the U band to the Spitzer 8 \mu m band). We have derived a catalogue of galaxies complete in rest-frame B-band, which has been divided in two subsamples according to their rest-frame U-V colour (or derived specific star formation rate, SSFR) properties. We confirm a bimodality in the U-V colour and SSFR of the galaxy sample up to z\sim 3. This bimodality is used to compute the LFs of the blue/late and red/early subsamples. The LFs of the blue/late and total samples are well represented by steep Schechter functions evolving in luminosity with increasing redshifts. The volume density of the LFs of the red/early populations decreases with increasing redshift. The shape of the red/early LFs shows an excess of faint red dwarfs with respect to the extrapolation of a flat Schechter function and can be represented by the sum of two Schechter functions. Our model for galaxy formation in the hierarchical clustering scenario, which also includes external feedback due to a diffuse UV background, shows a general broad agreement with the LFs of both populations, the larger discrepancies being present at the faint end for the red population. Hints on the nature of the red dwarf population are given on the basis of their stellar mass and spatial distributions.

We present a 900 sec, wide-field U image of the inner region of the Andromeda galaxy obtained during the commissioning of the blue channel of the Large Binocular Camera mounted on the prime focus of the Large Binocular Telescope. Relative photometry and absolute astrometry of individual sources in the image was obtained along with morphological parameters aimed at discriminating between stars and extended sources, e.g. globular clusters. The image unveils the near-ultraviolet view of the inner ring of star formation recently discovered in the infrared by the Spitzer Space Telescope and shows in great detail the fine structure of the dust lanes associated with the galaxy inner spiral arms. The capabilities of the blue channel of the Large Binocular Camera at the Large Binocular Telescope (LBC-Blue) are probed by direct comparison with ultraviolet GALEX observations of the same region in M31. We discovered 6 new candidate stellar clusters in this high-background region of M31. We also recovered 62 bona-fide globulars and 62 previously known candidates from the Revised Bologna Catalogue of the M31 globular clusters, and firmly established the extended nature of 19 of them.

By their star content, the galaxies split out into a red and a blue population; their color index peaked around u-r=2.5 or u-r=1, respectively, quantifies the ratio of the blue stars newly formed from cold galactic gas, to the redder ones left over by past generations. On the other hand, upon accreting substantial gas amounts the central massive black holes energize active galactic nuclei (AGNs); here we investigate whether these show a similar, and possibly related, bimodal partition as for current accretion activity relative to the past. To this aim we use an updated semianalytic model; based on Monte Carlo simulations, this follows with a large statistics the galaxy assemblage, the star generations and the black hole accretions in the cosmological framework over the redshift span from z=10 to z=0. We test our simulations for yielding in close detail the observed split of galaxies into a red, early and a blue, late population. We find that the black hole accretion activities likewise give rise to two source populations: early, bright quasars and later, dimmer AGNs. We predict for their Eddington parameter $\lambda_E$ -- the ratio of the current to the past black hole accretions -- a bimodal distribution; the two branches sit now under $\lambda_E \approx 0.01$ (mainly contributed by low-luminosity AGNs) and around $\lambda_E \approx 0.3-1$. These not only mark out the two populations of AGNs, but also will turn out to correlate strongly with the red or blue color of their host galaxies.

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

We investigate the effect of AGN feedback associated to the bright QSO phase onto the color distribution of galaxies from z=0 up to z=4. To this aim, we insert a blast-wave model of AGN feedback in our semi-analytic model of galaxy formation, which includes the growth of supermassive black holes and the AGN activity triggered by interactions of the host galaxies. The AGN feedback is directly related to the impulsive, luminous quasar phase. We test our model by checking the consistency of its results against i) the QSO luminosity functions from z=0 to z=4; ii) the observed local relation between the black hole mass m_{BH} and the mass of the host galaxy. At low redshift the inclusion of AGN feedback enhances the number of red bright galaxies, so that the color distribution of M_r<-22 objects is entirely dominated by red (u-r>1.5) galaxies; at 0.5< z< 2 it yields rest-frame U-V color distribution in agreement with existing observations. In the range z=1.5-2.5, we find that 31% of galaxies contribute to the EROs population with m_K<20 (Vega system); at such a magnitude, the model yields an EROs surface density of 6.3 10^3 /deg^2 matching existing data; about 80 % of the above density is contributed by passive galaxies with old stellar populations. Extending our analysis to z=4, the model matches the observed surface density 1.5 10^3 /deg^2 of DRGs at m_K=20; such a population is predicted to be dominated by galaxies with old stellar populations for z>2.5.

We use the GOODS-MUSIC sample, a catalog of ~3000 Ks-selected galaxies based on VLT and HST observation of the GOODS-South field with extended multi-wavelength coverage (from 0.3 to 8 micron) and accurate estimates of the photometric redshifts to select 179 DRGs with J-Ks>1.3 in an area of 135 sq. arcmin. We first show that the J-Ks>1.3 criterion selects a rather heterogeneous sample of galaxies, going from the targeted high-redshift luminous evolved systems, to a significant fraction of lower redshift (1<z<2) and less luminous dusty starbursts. These low-redshift DRGs are significantly less clustered than higher-z DRGs. With the aid of extreme and simplified theoretical models of clustering evolution we show that it is unlikely that the two samples are drawn from the same population observed at two different stages of evolution. High-z DRGs likely represent the progenitors of the more massive and more luminous galaxies in the local Universe and might mark the regions that will later evolve into structures of intermediate mass, like groups or small galaxy clusters. Low-z DRGs, on the other hand, will likely evolve into slightly less massive field galaxies.

We show how the observed bimodality in the color distribution of galaxies can be explained in the framework of the hierarchical clustering picture in terms of the interplay between the properties of the merging histories and the feedback/star-formation processes in the progenitors of local galaxies. Using a semi-analytic model of hierarchical galaxy formation, we compute the color distributions of galaxies with different luminosities and compare them with the observations. Our fiducial model matches the fundamental properties of the observed distributions, namely: 1) the distribution of objects brighter than M_r = -18 is clearly bimodal, with a fraction of red objects increasing with luminosity; 2) for objects brighter than M_r = -21 the color distribution is dominated by red objects with color u-r = 2.2-2.4; 3) the spread on the distribution of the red population is smaller than that of the blue population; 4) the fraction of red galaxies is larger in denser environments, even for low-luminosity objects; 5) the bimodality in the distribution persists up to z = 1.5. We discuss the role of the different physical processes included in the model in producing the above results.

The K20 survey is a near infrared-selected, deep (Ks < 20) redshift survey targeting galaxies in two independent regions of the sky, the CDFS and the q0055-2659 field. The total Ks-selected sample includes 545 objects. Optical spectra for 525 of them have been obtained, providing 501 spectroscopic identifications (including 12 type-1 AGN and 45 stars). This paper describes the final K20 spectroscopic catalogue, along with the technique used to determine redshifts, measure the spectral features and characterize the spectra. The classification of the galaxy spectra has been performed according to a simple parametric recipe that uses the equivalent widths of the two main emission lines (O[II] and Halpha+[N II]) and two continuum indices (the 4000A break index, and a near-UV color index, C(28-39)). We defined three main spectroscopic classes: red early-type galaxies, blue emission-line galaxies and the intermediate galaxies, which show emission lines but a red continuum. More than 95% of the examined galaxies is included in one of these spectral types and a composite spectrum is built for each of the three galaxy classes. The blue emission-line and the early-type galaxies have been divided in redshift bins, and the corresponding composite spectra have been built, in order to investigate the evolution of the spectral properties of the K20 galaxies with redshift.

We investigate how the hierarchical merging of dark matter halos, the radiative cooling of baryons, and the energy feedback from supernovae and active galactic nuclei or quasars combine to govern the amount and the thermal state of the hot plasma pervading groups and clusters of galaxies. We show that, by itself, supernova preheating of the external gas flowing into clusters falls short of explaining the observed X-ray scaling relations of the plasma luminosity L_X or the plasma entropy K versus the X-ray temperature T. To account for the scaling laws from rich to poor clusters takes preheating enhanced by the energy input from active galactic nuclei. In groups, on the other hand, the internal impacts of powerful quasars going off in member galaxies can blow some plasma out of the structure. So they depress L_X and raise K to the observed average levels; meanwhile, the sporadic nature of such impulsive events generates the intrinsic component of the wide scatter apparent in the data. The same quasar feedback gives rise in groups to entropy profiles as steep as observed, a feature hard to explain with simple preheating schemes. Finally, we argue a close connection of the L_X-T or the K-T relation with the M_{\bullet}-\sigma correlation between the host velocity dispersion and the masses of the black holes, relics of the quasar activity.

We have explored the redshift evolution of the luminosity function of red and blue galaxies up to $z=3.5$. This was possible joining a deep I band composite galaxy sample, which includes the spectroscopic K20 sample and the HDFs samples, with the deep $H_{AB}=26$ and $K_{AB}=25$ samples derived from the deep NIR images of the Hubble Deep Fields North and South, respectively. About 30% of the sample has spectroscopic redshifts and the remaining fraction well-calibrated photometric redshifts. This allowed to select and measure galaxies in the rest-frame blue magnitude up to $z\sim 3$ and to derive the redshift evolution of the B-band luminosity function of galaxies separated by their rest-frame $U-V$ color or specific (i.e. per unit mass) star-formation rate. The class separation was derived from passive evolutionary tracks or from their observed bimodal distributions. Both distributions appear bimodal at least up to $z\sim 2$ and the locus of red/early galaxies is clearly identified up to these high redshifts. Both luminosity and density evolutions are needed to describe the cosmological behaviour of the red/early and blue/late populations. The density evolution is greater for the early population with a decrease by one order of magnitude at $z\sim 2-3$ with respect to the value at $z\sim 0.4$. The luminosity densities of the early and late type galaxies with $M_B<-20.6$ appear to have a bifurcation at $z>1$. Indeed while star-forming galaxies slightly increase or keep constant their luminosity density, "early" galaxies decrease in their luminosity density by a factor $\sim 5-6$ from $z\sim 0.4$ to $z\sim 2.5-3$. A comparison with one of the latest versions of the hierarchical CDM models shows a broad agreement with the observed number and luminosity density evolutions of both populations.

We study the evolution of the cosmic star formation by computing the luminosity density (LD) in the UV, B, J, and K bands, and the stellar mass density (MD) of galaxies in two reference models of galaxy evolution: the pure luminosity evolution (PLE) model developed by Calura & Matteucci (2003) and the semi-analytical model (SAM) of hierarchical galaxy formation by Menci et al. (2002). The former includes a detailed description of the chemical evolution of galaxies of different morphological types with no density evolution; the latter includes the merging histories of the galactic DM haloes, as predicted by the hierarchical clustering scenario, but it does not contain morphological classification nor chemical evolution. We find that at z< 1.5 both models are consistent with the available data on the LD of galaxies in all the considered bands. At high z, the LDs predicted in the PLE model show a peak due to the formation of ellipticals, whereas the SAM predicts a gradual decrease of the star formation and of the LD for z> 2.5. At such redshifts the PLE predictions tend to overestimate the present data in the B band whereas the SAM tends to underestimate the observed UV LD. As for the stellar MD, the PLE picture predicts that nearly 50% and 85% of the present stellar mass are in place at z=4 and z=1, respectively. According to the SAM, 50% and 60% of the present stellar mass are in place at z=1.2 and z=1, respectively. Both predictions fit the observed MD up to z=1. At z>1, the PLE model and the SAM tend to overestimate and underestimate the observed values, respectively. We discuss the origin of the above model results, and the role of observational uncertainties (such as dust extinction) in comparing models with observations.

We present a detailed analysis of the stellar mass content of galaxies up to z=2.5 in the K20 galaxy sample, that has a 92% spectroscopic completeness and a complete $UBVRIzJK_s$ multicolor coverage. We find that the M/L ratio decreases with redshift: in particular, the average M/L ratio of early type galaxies decreases with $z$, with a scatter that is indicative of a range of star--formation time-scales and redshift of formation. More important, the typical M/L of massive early type galaxies is larger than that of less massive ones, suggesting that their stellar population formed at higher z. The final K20 galaxy sample spans a range of stellar masses from M*=10^9Msun to M*=10^12Msun, with massive galaxies ($M*>10^11Msun) detected up to z~2. We compute the Galaxy Stellar Mass Function at various z, of which we observe only a mild evolution (i.e. by 20-30%) up to z~1. At z>1, the evolution of the GSMF appears to be much faster: at z~2, about 35% of the present day stellar mass in objects with M*~10^11Msun appear to have assembled. We also detect a change in the physical nature of the most massive galaxies, since at z>1 a population of massive star--forming galaxies progressively appears. We finally analyze our results in the framework of Lambda-CDM hierarchical models. First, we show that the large number of massive galaxies detected at high z does not violate any fundamental Lambda-CDM constraint based on the number of massive DM halos. Then, we compare our results with the predictions of renditions of both semianalytic and hydro-dynamical models, that range from severe underestimates to slight overestimates of the observed mass density at z<~2. We discuss how the differences among these models are due to the different implementation of the main physical processes. (Abridged)

We have incorporated the description of the X-ray properties of Active Galactic Nuclei (AGNs) into a semi-analytic model of galaxy formation, adopting physically motivated scaling laws for accretion triggered by galaxy encounters. Our model reproduces the level of the cosmic X-ray background at 30 keV; we predict that the largest contribution (around 2/3) comes from sources with intermediate X-ray luminosity 10^{43.5}< L_X/erg/s <10^{44.5}, with 50 % of the total specific intensity produced at z<2. The predicted number density of luminous X-ray AGNs (L_X>10^{44.5} erg/s in the 2-10 keV band) peaks at z around 2 with a decline of around 3 dex to z=0; for the low luminosity sources (10^{43}<L_X/erg/s <10^{44}) it has a broaderand less pronounced maximum around z 1.5. The comparison with the data shows a generally good agreement. The model predictions slightly exceed the observed number of low-luminosity AGNs at z around 1.5, with the discrepancy progressively extending to intermediate-luminosity objects at higher redshifts; we discuss possible origins for the mismatch. Finally, we predict the source counts and the flux distribution at different redshifts in the hard (20-100 keV) X-ray band for the sources contributing to the X- ray background.

To address the problem concerning the early formation of stars in massive galaxies, we present the results of a semi-analytic model of galaxy formation which includes a physical description of starbursts triggered by galaxy interactions. These originate from the destabilization of cold galactic gas occurring in galaxy encounters, which in part feeds the accretion onto black holes powering quasars, and in part drives circumnuclear starsbursts at redshifts $z\approx 2-4$, preferentially in massive objects. This speeds up the formation of stars in massive galaxies at high redshifts without altering it in low mass galactic halos. Thus, at intermediate $z\approx 1.5-2$ we find that a considerable fraction of the stellar content of massive galaxies is already in place, at variance with the predictions of previous hierarchical models. The resulting high-$z$ star formation rate and B-band luminosity functions, and the luminosity and redshift distribution of galaxies in K-band at $z\lesssim 2$ are all in good agreement with the existing observations concerning the bright galaxy population.

We use a deep K(AB)<25 galaxy sample in the Hubble Deep Field South to trace the evolution of the cosmological stellar mass density from z~ 0.5 to z~3. We find clear evidence for a decrease of the average stellar mass density at high redshift, 2<z<3.2, that is 15^{+25}_{-5}% of the local value, two times higher than what observed in the Hubble Deep Field North. To take into account for the selection effects, we define a homogeneous subsample of galaxies with 10^{10}M_\odot \leq M_* \leq 10^{11}M_\odot: in this sample, the mass density at z>2 is 20^{+20}_{-5} % of the local value. In the mass--limited subsample at z>2, the fraction of passively fading galaxies is at most 25%, although they can contribute up to about 40% of the stellar mass density. On the other hand, star--forming galaxies at z>2 form stars with an average specific rate at least <\dot M/M_*> ~4 x10^{-10} yr$^{-1}$, 3 times higher than the z<~1 value. This implies that UV bright star--forming galaxies are substancial contributors to the rise of the stellar mass density with cosmic time. Although these results are globally consistent with $\Lambda$--CDM scenarios, the present rendition of semi analytic models fails to match the stellar mass density produced by more massive galaxies present at z>2.

We present an estimate of the cosmological evolution of the field galaxy luminosity function (LF) in the rest frame 4400 Angstrom B -band up to redshift z=3.5. To this purpose, we use a composite sample of 1541 I--selected galaxies selected down to I_(AB)=27.2 and 138 galaxies selected down to K_(AB)=25 from ground-based and HST multicolor surveys, most notably the new deep JHK images in the Hubble Deep Field South (HDF-S) taken with the ISAAC instrument at the ESO-VLT telescope. About 21% of the sample has spectroscopic redshifts, and the remaining fraction well calibrated photometric redshifts. The resulting blue LF shows little density evolution at the faint end with respect to the local values, while at the bright end (M_B(AB)<-20) a brightening increasing with redshift is apparent with respect to the local LF. Hierarchical CDM models overpredict the number of faint galaxies by about a factor 3 at z=1. At the bright end the predicted LFs are in reasonable agreement only at low and intermediate redshifts (z=1), but fail to reproduce the pronounced brightening observed in the high redshift (z=2-3) LF. This brightening could mark the epoch where a major star formation activity is present in the galaxy evolution.

We present the galaxy rest-frame near-IR Luminosity Function (LF) and its cosmic evolution to z=1.5 based on a spectroscopic survey of a magnitude limited sample of galaxies with Ks<20 (the K20 survey, Cimatti et al. 2002b). The LFs have been derived in the rest-frame J and Ks bands using 3 z bins (z_mean= 0.5, 1, 1.5) and compared to the local near-IR LF. The faint-end of the LFs is consistent with the local estimates, with no evidence for a change either in the slope or normalization up to z<1.3. Viceversa, the density of luminous galaxies (M_Ks-5logh_70<-25.5) is higher than locally at all z and relatively constant or mildly increasing with z within our sample. The data are consistent with a mild luminosity evolution both in the J and Ks-band up to z=1.5 (DeltaM_J=-0.69+-0.12 and DeltaM_K=-0.54+-0.12 at z=1). Moreover, we find that red and early-type galaxies dominate the bright-end of the LF, and that their number density shows at most a small decrease (<30%) up to z=1, thus suggesting that massive elliptical galaxies were already in place at z=1 and they should have formed their stars and assembled their mass at higher z. There appears to be a correlation of the optical/near-IR colors with near-IR luminosities, the most luminous/massive galaxies being red/old. We find a slow evolution with z of the near-IR comoving luminosity density to z=1.5. Finally, we show that hierarchical models overpredict significantly the density of low luminosity galaxies at z<=1 and underpredict the density of luminous galaxies at z>=1, whereas PLE models are more consistent with the data up to z=1.5. The GIF model (Kaufmann et al. 1999) shows a clear deficiency of red luminous galaxies at z=1 compared to our observations and predicts a decrease of luminous galaxies with z not observed in our sample.

We link the evolution of the galaxies in the hierarchical clustering scenario with the changing accretion rates of cold gas onto the central massive black holes that power the quasars. We base on galaxy interactions as main triggers of accretion; the related scaling laws are taken up from Cavaliere & Vittorini (2000), and grafted to a semi-analytic code for galaxy formation. As a result, at high $z$ the protogalaxies grow rapidly by hierarchical merging; meanwhile, much fresh gas is imported and also destabilized, so the holes are fueled at their full Eddington rates. At lower $z$ the galactic dynamical events are mostly encounters in hierarchically growing groups; now the refueling peters out, as the residual gas is exhausted while the destabilizing encounters dwindle. So, with no parameter tuning other than needed for stellar observables, our model uniquely produces at $z>3$ a rise, and at $z\lesssim 2.5 $ a decline of the bright quasar population as steep as observed. In addition, our results closely fit the observed luminosity functions of quasars, their space density at different magnitudes from $z\approx 5$ to $z\approx 0$, and the local $m_{BH}-\sigma$ relation.

We discuss Tree+SPH simulations of galaxy clusters and groups, aimed at studying the effect of cooling and non-gravitational heating on observable properties of the ICM. We simulate at high resolution four halos,with masses in the range (0.2-4)10^{14}M_sol. We discuss the effects of using different SPH implementations and show that high resolution is mandatory to correctly follow the cooling pattern of the ICM. All of our heating schemes which correctly reproduce the X-ray scaling properties of clusters and groups do not succeed in reducing the fraction of collapsed gas below a level of 20 (30) per cent at the cluster (group) scale. Finally, gas compression in cooling cluster regions causes an increase of the temperature and a steepening of the temperature profiles, independent of the presence of non-gravitational heating processes. This is inconsistent with recent observational evidence for a decrease of gas temperature towards the center of relaxed clusters. Provided these discrepancies persist even for a more refined modeling of energy feedback, they may indicate that some basic physical process is still missing in hydrodynamical simulations.

We present the results of a VLT and HST imaging survey aimed at the identification of 4.5<z<6 galaxies. In the VLT data, a set of broad and intermediate band filters has been used to select 13 high--$z$ candidates in a Z(AB)<25 catalog, over an area of about 30 arcmin2. Discrimination against lower redshift interlopers (mainly early--type galaxies at high redshift and cool Galactic stars) has been done combining morphological and spectral classification. This sample has been combined with a deeper I<27.2 sample obtained from the Hubble Deep Field campaigns. The VLT final sample consists of 13 high z candidates, 4 of which are identified with high confidence as z>4.5 galaxies. The resulting integral surface density of the Z<25 candidates at z>4.5z is in the range 0.13-0.44/arcmin^2 and that in the highest redshift bin 5<z<6 is between 0.07-0.13/arcmin^2. In the two HDFs, we identify at I<27.2 25 galaxies in the range 4.5<z<5 and 16 at 5<z<6$, corresponding to surface densities of 3.1 arcmin^-2 and 2 arcmin^-2, respectively. We show that the observed Z<25 UV luminosity density appear to drop by about one order of magnitude from z~3 to z~6. However, if we apply a threshold to obtain an absolute--magnitude limited sample, the UV luminosity density results to be roughly constant up to z~6. We finally show that recent semi-analytic hierarchical models for galaxy formation, while predicting a nearly constant total UV luminosity density up to z~6, under-predict the observed UV luminosity density at Z<25 and over-predicts the I<27.2 one. This behaviour can be understood in term of a poor match to the slope of the UV luminosity function.

Galaxy groups are quite underluminous in X-rays compared to clusters, so the intracluster medium has to be considerably underdense in the former. We consider this to be due to substantial energy fed back into the ICM when the baryons in the member galaxies condense into stars ending up in SNe, or accrete on to central supermassive black holes energizing AGNs. We compute the outflow and the blowout effects driven by the AGNs and the resulting, steep luminosity-temperature correlation L_X-T. We compare this with the SN contribution and with the X-ray data; the latter require the AGN energy to be coupled to the surrounding ICM at fractional levels around 5 10^{-2}. We link the L_X-T behavior with the parallel effects of the AGN feedback on the gas in the host galaxy; we find that these yield a correlation steep up to M_{\bullet}\propto \sigma^5 between the galactic velocity dispersions and the central BH masses.

We present the redshift distribution of a complete sample of 480 galaxies with Ks<20 distributed over two independent fields covering a total area of 52 arcmin^2. The redshift completeness is 87% and 98% respectively with spectroscopic and high-quality and tested photometric redshifts. The redshift distribution of field galaxies has a median redshift z_{med}\sim 0.80, with \sim 32% and \sim 9% of galaxies at z>1 and z>1.5 respectively. A ``blind'' comparison is made with the predictions of a set of the most recent LambdaCDM hierarchical merging and pure luminosity evolution (PLE) models. The hierarchical merging models overpredict and underpredict the number of galaxies at low-z and high-z respectively, whereas the PLE models match the median redshift and the low-z distribution, still being able to follow the high-z tail of N(z). We briefly discuss the implications of this comparison and the possible origins of the observed discrepancies. We make the redshift distribution publicly available.

The K20 survey is an ESO VLT optical and near-infrared spectroscopic survey aimed at obtaining spectral information and redshifts of a complete sample of about 550 objects to K_s\leq20.0 over two independent fields with a total area of 52 arcmin^2. In this paper we discuss the scientific motivation of such a survey, we describe the photometric and spectroscopic properties of the sample, and we release the $K_s$-band photometric catalog. Extensive simulations showed that the sample is photometrically highly complete to K_s=20. The observed galaxy counts and the R-K_s color distribution are consistent with literature results. We observed spectroscopically 94% of the sample, reaching a spectroscopic redshift identification completeness of 92% to K_s\leq20.0 for the observed targets, and of 87% for the whole sample (i.e. counting also the unobserved targets). Deep spectroscopy was complemented with multi-band deep imaging in order to derive tested and reliable photometric redshifts for the galaxies lacking spectroscopic redshifts. The results show a very good agreement between the spectroscopic and the photometric redshifts with <z_{spe}-z_{phot}>=0.01 and with a dispersion of \sigma_{\Delta z}=0.09. Using both the spectroscopic and the photometric redshifts, we reached an overall redshift completeness of about 98%. The size of the sample, the redshift completeness, the availability of high quality photometric redshifts and multicolor spectral energy distributions make the K20 survey database one of the most complete samples available to date for constraining the currently competing scenarios of galaxy formation and for a variety of other galaxy evolution studies.

We present a set of gas-dynamical simulations of galaxy groups and clusters aimed at exploring the effect of non-gravitational heating. We use GASOLINE, a parallel Tree+SPH code, to simulate the formation of four cosmic halos with temperature 0.5<T<8 keV. Non-gravitational heating is implemented in two different ways: (1) by imposing a minimum entropy floor at a given redshift, 1<z<5; (2) by gradually heating gas, proportionally to the SN rate expected from semi-analytical modeling of galaxy formation. Our main results are the following. (a) An extra heating energy of about 1 keV per gas particle is required to reproduce the observed Lx-T relation, independent of whether it is provided so as to create an entropy floor of 50-100 keV cm^2, or is modulated in redshift; our SN feedback recipe provides only 1/3 keV/part. (b) The M-T relation is almost unaffected by non-gravitational heating and follows the M T^{3/2} scaling, with a normalization ~40% higher than observed, independent of the heating scheme. The inclusion of cooling in a run of a small group has the effects of increasing T_ew by ~30%, possibly reconciling simulated and observed M-T relations, and of decreasing Lx by ~40%. In spite of the inclusion of SN feedback energy, almost 40% of the gas becomes cold, in excess of current observational estimates. (abridged)

We develop a semi-analytic model of hierarchical galaxy formation with an improved treatment of the evolution of galaxies inside dark matter haloes. We take into account not only dynamical friction processes building up the central dominant galaxy, but also binary aggregations of satellite galaxies inside a common halo described using the kinetic Smoluchowski equation. The description of gas cooling, star formation and evolution, and Supernova feedback follows the standard prescriptions widely used in semi-analytic modelling. We find that binary aggregations are effective in depleting the number of small/intermediate mass galaxies over the redshift range 1<z<3, thus flattening the slope of the luminosity function at the faint end. At z~0 the flattening occurs for -20<M_B<-18, but an upturn is obtained at the very faint end for M_B>-16. We compare our predicted luminosity functions with those obtained from deep multicolor surveys in the rest-frame B and UV bands for the redshift ranges 0<z<1 and 2.5<z<3.5, respectively. The comparison shows that the discrepancy of the predictions of other semi-analytic models with the observations is considerably reduced at z>1 and even more at z ~ 3 by the effect of binary aggregations. The predictions from our dynamical model are discussed and compared with the effects of complementary processes which may conspire in affecting the shape of the luminosity function.

We compare the 3D clustering of old passively-evolving and dusty star-forming z~1 EROs from the K20 survey. With detailed simulations of clustering, the comoving correlation length of dusty star-forming EROs is constrained to be less than r_0~2.5 Mpc/h. In contrast, the old EROs are much more positively correlated, with 5.5<~r_0/(Mpc/h)<~16, consistent with previous claims for z~1 field early-type galaxies based on analyses of ERO angular clustering. The low level of clustering of dusty star-forming EROs does not support these to be major mergers building up an elliptical galaxy, or typical counterparts of SCUBA sources, but it is instead consistent with the weak clustering of high redshift blue galaxies and of luminous local IRAS galaxies. Current hierarchical merging models can explain the large r_0 for z~1 field early-type galaxies, but fail in matching their high number density and overall old ages.

We present the results of VLT optical spectroscopy of a complete sample of 78 EROs with R-Ks\geq5 over a field of 52 arcmin^2. About 70% of the 45 EROs with Ks\leq19.2 have been spectroscopically identified with old passively evolving and dusty star-forming galaxies at 0.7<z<1.5. The two classes are about equally populated and for each of them we present and discuss the average spectrum. From the old ERO average spectrum and for Z=Z_{\odot} we derive a minimum age of \sim 3 Gyr, corresponding to a formation redshift of z_f \gtsima 2.4. PLE models with such formation redshifts well reproduce the density of old EROs (consistent with being passively evolving ellipticals), whereas the predictions of the current hierarchical merging models are lower than the observed densities by large factors (up to an order of magnitude). From the average spectrum of the star-forming EROs we estimate a substantial dust extinction with E(B-V) \gtsima 0.5. The star formation rates, corrected for the average reddening, suggest a significant contribution from EROs to the cosmic star-formation density at z \sim 1.

A description of the dynamics of a collisionless, self-gravitating fluid is developed and applied to follow the development of Large Scale Structures in the Universe. Such description takes on one of the assumptions of the Adhesion Approximation (AA) model, i.e., the introduction of an artificial viscosity $\nu$ in the Euler equation, but extends it to deeper non-linear stages, where the extrapolation of the linear relation $\psi=-\phi$ between the velocity and the gravitational potentials -- at the basis of both the Zel'dovich and the Adhesion models -- is no longer valid. This is achieved by expanding the relation bewteen the gravitational potential $\phi$ and the velocity potential $\psi$, in powers of the small viscosity $\nu$. When the results are compared to the outcomes of an N-body simulation, the proposed extension of the Adhesion approach provides a better description of small-scale, deeply non-linear regions, as it is quantitatively shown by the computation of some statistical indicators. At larger scales, the satisfactory description of the large scale texture and of the voids given by the canonical Adhesion Approximation is preserved in the extended model.

We present results from high-resolution Tree+SPH simulations of galaxy clusters and groups, aimed at studying the effect of non-gravitational heating on the entropy of the ICM. We simulate three systems, having emission-weighted temperature T=0.6,1 and 3 keV, with spatial resolution better than 1% of the virial radius. We consider the effect of different prescriptions for non-gravitational ICM heating, such as SN energy feedback, as predicted by semi-analytical models of galaxy formation, and two different minimum entropy floors, S_fl=50 and 100 keV cm^2, imposed at z=3. Simulations with only gravitational heating nicely reproduce predictions from self-similar ICM models, while extra heating is shown to break the self-similarity, by a degree which depends on total injected energy and on cluster mass. We use observational results on the excess entropy in central regions of galaxy systems, to constrain the amount of extra-heating required. We find that setting the entropy floor S_fl=50 keV cm^2, which corresponds to an extra heating energy of about 1 keV per particle, is able to reproduce the observed excess of ICM entropy.

We investigate how strongly and at what scales the Sunyaev- Zel'dovich effect reflects the shifting balance between the two processes that compete for governing the density and the thermodynamic state of the hot intra-cluster medium pervading clusters and groups of galaxies: the hierarchical clustering of the DM; the non-gravitational energy and momentum fed back into the ICM by the condensing baryons. We base on a SAM of galaxy formation and clustering to describe how the baryons are partitioned among the hot, the cool and the stellar phase; the partition shifts as the galaxies cluster hierarchically, and as the feedback by stellar winds and SN explosions follows the star formation. Their impact is amplified by the same large scale accretion shocks that thermalize the gravitational energy of gas falling into the growing potential wells. We compute the Compton parameter $y$, and find a relation of $y$ with the ICM temperature, the $y-T$ relation, which departs from the self-similar scaling and bends down at temperatures typical of galaxy groups. We model-independently relate this with the analogous behaviour of the L_x - T relation, and discuss to what extent our results are generic of the hierarchical models of galaxy formation and clustering.

The galaxy Luminosity Function (LF) has been estimated in the rest frame B luminosity at 0<z<1.25 and at 1700 {\AA} for 2.5<z<4.5 from deep multicolor surveys in the HDF-N, HDF-S, NTT-DF. The results have been compared with a recent version of galaxy formation models in the framework of hierarchical clustering in a flat Cold Dark Matter Universe with cosmological constant. The results show a general agreement for z<= 1, although the model LF has a steeper average slope at the faint end; at z~3 such feature results in an overprediction of the number of faint (I_{AB}~ 27) galaxies, while the agreement at the bright end becomes critically sensitive to the details of dust absorption at such redshifts. The discrepancies at the faint end show that a refined treatement of the physical processes involving smaller galaxies is to be pursued in the models, in terms of aggregation processes and/or stellar feedback heavily affecting the luminosity of the low luminosity objects. The implications of our results on the evolution of the cosmological star formation rate are discussed.

The intrinsic sizes of the field galaxies with I<26 in the Hubble and ESO-NTT Deep Fields are shown as a function of their redshifts and absolute magnitudes using photometric redshifts derived from the multicolor catalogs and are compared with the CDM predictions. Extending to lower luminosities and to higher z our previous analysis performed on the NTT field alone, we find that the distribution of the galaxy disk sizes at different cosmic epochs is within the range predicted by typical CDM models. However, the observed size distribution of faint (M_B>-19) galaxies is skewed with respect to the CDM predictions and an excess of small-size disks (R_d<2 kpc) is already present at z~ 0.5. The excess persists up to z~3 and involves brighter galaxies . Such an excess may be reduced if luminosity-dependent effects, like starburst activity in interacting galaxies, are included in the physical mechanisms governing the star formation history in CDM models.

The current models of early star and galaxy formation are based upon the hierarchical growth of dark matter halos, within which the baryons condense into stars after cooling down from a hot diffuse phase. The latter is replenished by infall of outer gas into the halo potential wells; this includes a fraction previously expelled and preheated, due to momentum and energy fed back by the SNe which follow the star formation. We identify such an implied hot phase with the medium known to radiate powerful X-rays in clusters and in groups of galaxies. We show that the amount of the hot component required by the current star formation models is enough to be observable out to redshifts $z \approx 1.5$ in forthcoming deep surveys from {\it Chandra} and {\it XMM}, especially in case the star formation rate is high at such and earlier $z$. These X-ray emissions constitute a necessary counterpart, and will provide a much wanted probe of the SF process itself (in particular, of the SN feedback), to parallel and complement the currently debated data from optical and IR observations of the young stars.

The sizes of the field galaxies with I<25 have been measured in the NTT Deep Field. Intrinsic sizes have been obtained after deconvolution of the PSF with a multigaussian method. The reliability of the method has been tested using both simulated data and HST observations of the same field. The distribution of the half light radii is peaked at r_{hl} 0.3 arcsec, in good agreement with that derived from HST images at the same magnitude. An approximate morphological classification has been obtained using the asymmetry and concentration parameters. The intrinsic sizes of the galaxies are shown as a function of their redshifts and absolute magnitudes using photometric redshifts derived from the multicolor catalog. While the brighter galaxies with morphological parameters typical of the normal spirals show a flat distribution in the range r_{d}=1-6 kpc, the fainter population at 0.4<z<0.8 dominates at small sizes. To explore the significance of this behaviour, an analytical rendition of the standard CDM model for the disc size evolution has been computed. The model showing the best fit to the local luminosity function and the Tully-Fisher relation is able to reproduce at intermediate redshifts a size distribution in general agreement with the observations, although it tends to underestimate the number of galaxies fainter than M_B~ -19 with disk sizes r_d~ 1-2 kpc.

A combination of ground-based (NTT and VLT) and HST (HDF-N and HDF-S) public imaging surveys have been used to collect a sample of 1712 I-selected and 319 $K\leq 21$ galaxies. Photometric redshifts have been obtained for all these galaxies. The results have been compared with the prediction of an analytic rendition of the current CDM hierarchical models for galaxy formation. We focus in particular on two observed quantities: the galaxy redshift distribution at K<21 and the evolution of the UV luminosity density. The derived photometric redshift distribution is in agreement with the hierarchical CDM prediction, with a fraction of only 5% of galaxies detected at z>2. This result strongly supports hierarchical scenarios where present-day massive galaxies are the result of merging processes. The observed UV luminosity density in the I-selected sample is confined within a factor of 4 over the whole range 0<z<4.5. CDM models in a critical Universe are not able to produce the density of UV photons that is observed at z>3. CDM models in $\Lambda$-dominated universe are in better agreement at 3<z<4.5, but predict a pronounced peak at z~1.5 and a drop by a factor of 8 from z=1.5 to z=4 that is not observed in the data. We conclude that improvements are required in the treatment of the physical processes directly related to the SFR, e.g. the starbust activity in merger processes and/or different feedback to the star formation activity.

We relate the star formation from cold baryons in virialized structures to the X-ray properties of the associated diffuse, hot baryonic component. Our computations use the standard ``semi-analytic'' models to describe i) the evolution of dark matter halos through merging after the hierarchical clustering, ii) the star formation governed by radiative cooling and by supernova feedback, iii) the hydro- and thermodynamics of the hot gas, rendered with our Punctuated Equilibria model. So we relate the X-ray observables concerning the intra-cluster medium to the thermal energy of the gas pre-heated and expelled by supernovae following star formation, and then accreted during the subsequent merging events. We show that at fluxes fainter than $F_X\approx 10^{-15}$ erg/cm$^2 $ s (well within the reach of next generation X-ray observatories) the X-ray counts of extended extragalactic sources (as well as the faint end of the luminosity function, the contribution to the soft X-ray background, and the $L_X-T$ correlation at the group scales) increase considerably when the star formation rate is enhanced for z>1 as indicated by growing optical/infrared evidence. Specifically, the counts in the range 0.5-2 keV are increased by factors $\sim 4$ when the the feedback is decreased and star formation is enhanced as to yield a flat shape of the star formation rate for 2<z<4.

Using current semi-analytic models of structure and star formation based on hierarchical clustering, we connect the cosmic star formation rate (SFR) to the X-ray emission properties of groups and clusters of galaxies. We show that when the baryons cool down from a hot phase to condense into stars with a SFR flat for z>2, then such hot phase yields deep X-ray source counts (at fluxes F=10^{-15} erg/cm^2 s in the energy band 0.5-2 keV) about 5 times larger than in the case of a SFR peaked at z ~ 1.5. We also discuss the effect of the SFR on the shape and evolution of the X-ray L-T relation.

We develop a model to describe the evolution of the intra-cluster X-ray emitting baryons, as they are included in the dark matter potential wells of galaxy clusters evolving through subsequent merging events in the framework of hierarchical clustering. The gas is assumed to re-adjusts to a new hydrostatic equilibrium after each merging event. Before merging it is gravitationally heated at the local virial temperature when bound in subclusters; at early $z$ the gas is preheated by supernova activity following star formation. In detail, we compute analytically the following steps: the dynamic histories of dark matter halos with their merging events; the associated infall of gas into a halo, with compressions and shocks estabilishing the conditions at the cluster boundary; the updated disposition of the gas in the potential well matching such conditions; the statistical convolution of the key quantities over the merging histories. The model predicts the density and surface brightness profiles with no free parameters; the so-called $\beta$ parameter is itself an outcome of the model, and the polytropic index $\gamma$ is internally constrained to a narrow range. We obtain declining temperature profiles, and profiles for the density and for the surface brightness shallower in groups compared with clusters. Our model groups also contain a lower baryonic fraction on average, but with a scatter considerably larger than at cluster scales. Various statistics are obtained analytically upon averaging over the merging histories. In particular, we predict in different cosmologies the statistical correlation $L-T$ of luminosity with temperature; similarly we derivethe correlation $R_X-T$ for the size of the X-ray emitting region. The intrinsic scatter in both correlations is also predicted.

The X-ray emission from clusters of galaxies is one of the best observational probe to investigate the distribution of dark matter at intermediate and high redshifts. Since the disposition of the intracluster plasma (ICP) responsible of the emission is crucial to link X-ray properties to the global properties of the dark matter halos, we propose a semi--analytical approach for the diffuse baryons. This comprises the following blocks: Monte Carlo merging histories to describe the dynamics of dark matter halos; the central hydrostatic disposition for the ICP; conditions of shock, or of closely adiabatic compression at the boundary with the external gas, preheated by stellar energy feedbacks. From our model we predict the $L-T$ correlation, consistent with the data as for shape and scatter.

To predict the X-ray observables associated to the diffuse baryons in clusters of galaxies, we develop a new physical approach to model such a hot intra-cluster plasma. Such approach is based on punctuated equilibria. and comprises the following blocks: Monte Carlo ``merging histories'' of dark matter potential wells; the central hydrostatic disposition for the ICP, reset to a new equilibrium after each merging episode; conditions of shock, or of closely adiabatic compression at the boundary with the external gas, preheated by stellar energy feedbacks. We predict the L-T relation, consistent with the data as for shape and scatter. This we combine with the mass distribution provided by the hierarchical clustering for different COBE-normalized CDM models, to predict the z-resolved luminosity functions, the source counts, the redshift distributions and contribution of the unresolved groups and clusters to the soft X-ray background. When compared with the recent ROSAT surveys, our results confirm that the critical cosmology with Standard CDM is ruled out by its overproduction of local clusters. On account of underproduction, instead, we rule out open cosmologies except for a narrow range around \Omega_o=0.5; even there, we find the consistency with the full data base to be hardly marginal. For \Omega_o=0.3 in flat geometry, we obtain acceptable fits. For the tilted CDM perturbation spectrum with high baryonic content in the critical universe, we obtain marginal consistency. Finally, we discuss the effective limitations of X-ray clusters and groups as cosmological signposts, and their brighter prospects toward the astrophysics of the ICP and the cosmogony of large, high-contrast structures.

We compute the evolution of the space-dependent mass distribution of galaxies in clusters due to binary aggregations by solving a space-dependent Smoluchowski equation. We derive the distribution of intergalactic distance for different ranges of mass (and of corresponding magnitude). We compare the results with the observed distributions, and find that the different degrees of luminosity segregation observed in clusters are well accounted for by our merging model. In addition, the presence of luminosity segregation is related to dynamical effects which also show up in different but connected observables, such as galaxy velocity profiles decreasing toward the center and X-ray measured beta-parameters smaller than 1. We predict both luminosity segregation and the observables above (being a product of binary aggregations) to be inversely correlated with the core radius and with the galaxy velocity dispersion; we discuss how the whole set of predictions compares with up-to-date observations.

We model the effects of shocks on the diffuse, X-ray emitting baryons in clusters of galaxies. Shocks separate the infalling from the inner gas nearly at equilibrium, and dominate the compression and the density gradients of the latter in the dark-matter potential of the cluster. We find that, independently of the detailed shape of the potential, the density gradient is steeper and the compression factor larger for the richer clusters. We show, considering the different merging histories, that in the hierarchical cosmogony the above effects lead, in X-rays, to a luminosity-temperature relation L \propto T^5 at the scale of groups which flattens down to L \propto T^{3} for rich clusters in accord with the observations, and then saturates toward L \propto T^2 for higher temperatures. From the merging histories we also compute statistical fluctuations of the L-T correlation.

We probe to what extent not only the counts of the faint blue galaxies and their redshift distribution but also their z-resolved luminosity functions can be explained in terms of binary merging (aggregations). We present a dynamical theory of such interactions. On this basis we find that ``minimal aggregations'' taking place within large scale structures and triggering starbusts yield rates and timing such as to explain the observations of the local, flat luminosity function and of its progressive rising and steepening for redshifts out to $z\approx 1$. Correspondingly, we predictfaint blue counts still rising up to $m_B\sim 28 -29$, redshift distributions shifting toward larger z with increasing $m_B$; in addition, we predict an upturn of the faint end of the luminosity function more pronounced in clusters than in the field. We propose that our picture provides the differential dynamics missing in the canonical hierarchical theories. This reconciles with the observations the steep luminosity functions predicted at high z by such theories.

We present a new approach to the statistics of the cosmic density field and to the mass distribution of high-contrast structures, based on the formalism of Cayley trees. Our approach includes in one random process both fluctuations and interactions of the density perturbations. We connect tree-related quantities, like the partition function or its generating function, to the mass distribution. The Press \& Schechter mass function and the Smoluchowski kinetic equation are naturally recovered as two limiting cases corresponding to independent Gaussian fluctuations, and to aggregation of high-contrast condensations, respectively. Numerical realizations of the complete random process on the tree yield an excess of large-mass objects relative to the Press \& Schechter function. When interactions are fully effective, a power-law distribution with logarithmic slope -2 is generated.