Context. As recently demonstrated, high-z submillimetre galaxies (SMGs) are the perfect background sample for tracing the mass density profiles of galaxies and clusters (baryonic and dark matter) and their time-evolution through gravitational lensing. Their magnification bias, a weak gravitational lensing effect, is a powerful tool for constraining the free parameters of a halo occupation distribution (HOD) model and potentially also some of the main cosmological parameters. Aims. The aim of this work is to test the capability of the magnification bias produced on high-z SMGs as a cosmological probe. We exploit cross-correlation data to constrain not only astrophysical parameters ($M_{min}$, $M_1$, and $\alpha$), but also some of the cosmological ones ($\Omega_m$, $\sigma_8$, and $H_0$) for this proof of concept. Methods. The measured cross-correlation function between a foreground sample of GAMA galaxies with spectroscopic redshifts in the range 0.2 < z < 0.8 and a background sample of H-ATLAS galaxies with photometric redshifts >1.2 is modelled using the traditional halo model description that depends on HOD and cosmological parameters. These parameters are then estimated by performing a Markov chain Monte Carlo analysis using different sets of priors to test the robustness of the results and to study the performance of this novel observable with the current set of data Results. With our current results, $\Omega_m$ and $H_0$ cannot be well constrained. However, we can set a lower limit of >0.24 at 95\% confidence level (CL) on $\Omega_m$ and we see a slight trend towards $H_0>70$ values. For our constraints on $\sigma_8$ we obtain only a tentative peak around 0.75, but an interesting upper limit of $\sigma_8\lesssim 1$ at 95\% CL. We also study the possibility to derive better constraints by imposing more restrictive priors on the astrophysical parameters.

Using the EMIR instrument on the IRAM 30m telescope, we conducted a spectroscopic redshift search of seven z$_{\rm phot}$ $\sim$ 4 sub-millimetre bright galaxies selected from the Herschel Bright Sources (HerBS) sample with fluxes at 500 $\mu$m greater than 80 mJy. For four sources, we obtained spectroscopic redshifts between 3.4 < z < 4.1 through the detection of multiple CO-spectral lines with J $\leq$ 3. Later, we detected low-J transitions for two of these sources with the GBT including the CO(1-0) transition. For the remaining three sources, more data are needed to determine the spectroscopic redshift unambiguously. The measured CO luminosities and line widths suggest that all these sources are gravitationally lensed. These observations demonstrate that the 2 mm window is indispensable to confirm robust spectroscopic redshifts for z < 4 sources. Finally, we present an efficient graphical method to correctly identify spectroscopic redshifts.

We generalize the analytic solutions presented in Pantoni et al. (2019) by including a simple yet effective description of wind recycling and galactic fountains, with the aim of self-consistently investigating the spatially-averaged time evolution of the gas, stellar, metal, and dust content in disc-dominated late-type galaxies (LTGs). Our analytic solutions, when supplemented with specific prescriptions for parameter setting and with halo accretion rates from $N-$body simulations, can be exploited to reproduce the main statistical relationships followed by local LTGs; these involve, as a function of the stellar mass, the star formation efficiency, the gas mass fraction, the gas/stellar metallicity, the dust mass, the star formation rate, the specific angular momentum, and the overall mass/metal budget. Our analytic solutions allow to easily disentangle the diverse role of the main physical processes ruling galaxy formation in LTGs; in particular, we highlight the crucial relevance of wind recycling and galactic fountains in efficiently refurnishing the gas mass, extending the star-formation timescale, and boosting the metal enrichment in gas and stars. All in all, our analytic solutions constitute a transparent, handy, and fast tool that can provide a basis for improving the (subgrid) physical recipes presently implemented in more sophisticated semi-analytic models and numerical simulations, and can offer a benchmark for interpreting and forecasting current and future spatially-averaged observations of local and higher redshift LTGs.

We investigate the shape of the Rotation Curves (RCs) of $z\sim 1$ star-forming galaxies and compare them with the local star-forming galaxies. For this purpose, we have used $409$ galaxies from the K-band Multi-Object Spectrograph (KMOS) for Redshift One Spectroscopic Survey (KROSS). This sample covers the redshift range $0.57\leq z \leq 1.04$, effective radii $0.69 \leq R_e \left[kpc \right] \leq 7.73$, absolute H-band magnitude $-24.46 \leq M_H \leq -18.85$ with median stellar mass $log\left(M_* \ \left[M_\odot \right]\right)=9.95$ and median total star-formation rate $log\left(SFR_{tot} \ \left[M_\odot \ yr^{-1} \right]\right)=1.49$. Using $^{3D}$BAROLO (Barolo), we extract $H\alpha$ kinematic maps and corresponding Rotation Curves (RCs). The main advantage of Barolo is that it incorporates the beam smearing in the 3D observational space, which provide us with the intrinsic rotation velocity even in the low spatial resolution data. Using Asymmetric Drift Correction (ADC), we have corrected the RCs for the pressure gradient effect, which seems to be a more dominant effect than beam smearing in high-z galaxies. Nearly all objects ($0.1 <v/ \sigma< 15$) are affected by the pressure gradient, and we noticed that ADC improves the rotation velocity of these systems by $\sim 10-87\%$. Only a combination of the three techniques (3D-kinematic modelling + 3D-Beamsmearing correction + ADC ) yields the intrinsic RC of an individual galaxy. Further, we present the co-added RCs constructed out of 237 high-quality objects to obtain intrinsic RC shapes out to 6.4 $\times$ disk scale length. We do not see any change in the shape of RCs with respect to the local star-forming disk-type galaxies. In contrast, we do find a significant evolution in the stellar-disk length ($R_D$) of the galaxies. Therefore, we conclude stellar disk evolves over cosmic time while total mass distribution stays constant.

The search for primordial gravitational waves through the $B$-mode polarization pattern in the CMB is one of the major goals of current and future CMB experiments. Besides foregrounds, a potential hurdle in this search is the anisotropic secondary $B$-mode polarization generated by the scattering of CMB photons off free electrons produced during patchy cosmological reionization. Robust predictions of these secondary anisotropies are challenging because of uncertainties in the reionization history. In this paper, we revise estimates of the reionization-induced $B$-mode signal by incorporating recent advances in the understanding of reionization through observations of the Lyman-$\alpha$ forest. To derive these $B$-mode estimates, we use high-dynamic-range radiative transfer simulations of reionization that are calibrated to the Ly$\alpha$ data. These simulations are also consistent with a variety of other high-redshift observations. We find that around multipoles $\ell\approx 100$, reionization induces $B$-mode power with $\ell(\ell+1)C_\ell^{BB}/2\pi\approx 4\times 10^{-6}\,\mu$K$^2$. This secondary signal is thus at the level of the primordial signal with the tensor-to-scalar ratio $r<10^{-4}$, and can increase by a factor of $\sim 50$ if reionization is sourced by highly clustered sources residing in haloes with mass of $\sim 10^{11}$ M$_\odot$. Our findings suggest that the contribution of patchy reionization to the search for primordial gravitational waves is unlikely to be a concern for currently planned CMB experiments.

The cross-correlation between fluctuations in the electron scattering optical depth $\tau_{\rm es}$ as probed by future Cosmic Microwave Background (CMB) experiments, and fluctuations in the 21cm differential brightness temperature $\Delta T_{\rm 21cm}$ as probed by ground-based radio interferometers, will trace the reionization history of the Universe. In particular, the $\tau_{\rm es}-$21cm cross-correlation should yield a determination of the characteristic bubble size distribution and ionization fraction as a function of redshift. When assuming that the cross-correlation signal is limited by instrumental noise rather than by foregrounds, we estimate its potential detectability by upcoming experiments. Specifically, the combination of HERA and Simons Observatory, CMB-S4 and PICO should yield a signal-to-noise ratio around 3 - 6, while and the exploitation of the SKA should increase it to 10-20. Finally, we have discussed how such levels of detectability can be affected when (simply modeled) 21cm foregrounds are present. For the most promising PICO$\times$SKA configuration, an efficiency of foreground removal to a level of $7\times 10^{-4}$ is needed to achieve a $5\sigma$ detection of the cross-correlation signal; in addition, safe avoidance of foreground contamination in the line-of-sight Fourier modes above $0.03 \,h\,\rm Mpc^{-1}$ would guarantee a detection significance around $3\sigma$.

We present a computational framework for "painting" galaxies on top of the Dark Matter Halo/Sub-Halo hierarchy obtained from N-body simulations. The method we use is based on the sub-halo clustering and abundance matching (SCAM) scheme which requires observations of the 1- and 2-point statistics of the target (observed) population we want to reproduce. This method is particularly tailored for high redshift studies and thereby relies on the observed high-redshift galaxy luminosity functions and correlation properties. The core functionalities are written in c++ and exploit Object Oriented Programming, with a wide use of polymorphism, to achieve flexibility and high computational efficiency. In order to have an easily accessible interface, all the libraries are wrapped in python and provided with an extensive documentation. We validate our results and provide a simple and quantitative application to reionization, with an investigation of physical quantities related to the galaxy population, ionization fraction and bubble size distribution.

We propose a new mechanism for the growth of supermassive black hole (BH) seeds in the star-forming progenitors of local early-type galaxies (ETGs) at $z\gtrsim 1$. This envisages the migration and merging of stellar compact remnants (neutron stars and stellar-mass BHs) via gaseous dynamical friction toward the central high-density regions of such galaxies. We show that, under reasonable assumptions and initial conditions, the process can build up central BH masses of order $10^4-10^6\, M_\odot$ within some $10^7$ yr, so effectively providing heavy seeds before standard disk (Eddington-like) accretion takes over to become the dominant process for further BH growth. Remarkably, such a mechanism may provide an explanation, alternative to super-Eddington accretion rates, for the buildup of billion solar masses BHs in quasar hosts at $z\gtrsim 7$, when the age of the Universe $\lesssim 0.8$ Gyr constitutes a demanding constraint; moreover, in more common ETG progenitors at redshift $z\sim 2-6$ it can concur with disk accretion to build such large BH masses even at moderate Eddington ratios $\lesssim 0.3$ within the short star-formation duration $\lesssim$ Gyr of these systems. Finally, we investigate the perspectives to detect the merger events between the migrating stellar remnants and the accumulating central supermassive BH via gravitational wave emission with future ground and space-based detectors such as the Einstein Telescope (ET) and the Laser Interferometer Space Antenna (LISA).

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.

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

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.

Magnification bias is a gravitational lensing effect that is normally overlooked because it is considered sub-optimal in comparison with the lensing shear. Thanks to the demonstrated optimal characteristics of the sub-millimetre galaxies (SMGs) for lensing analysis, in this work we were able to measure the magnification bias produced by a sample of QSOs acting as lenses, $0.2<z<1.0$, on the SMGs observed by Herschel at $1.2<z<4.0$. Two different methodologies were successfully applied: the traditional cross-correlation function approach and the Davis-Peebles estimator through stacking technique. The second one was found to be more robust for analysing the strong lensing regime ($<20-30$ arcsec in our case) and provides the possibility to take into account the positional errors of the sources in our samples. From the halo modelling of the cross-correlation function, the halo mass where the QSOs acting as lenses are located was estimated to be greater than $\log_{10}{(M_{min}/M_\odot)} > 13.6_{-0.4}^{+0.9}$, also confirmed by the mass density profile analysis ($M_{200c}\sim 10^{14} M_\odot$). These mass values indicate that we are observing the lensing effect of a cluster size halo signposted by the QSOs, as in previous studies of the magnification bias. Moreover, we were able to estimate the lensing convergence, $\kappa(\theta)$, for our magnification bias measurements down to a few kpcs. The derived mass density profile is in good agreement with a Navarro-Frank-White (NFW) profile. We also attempt an estimation of the halo mass and the concentration parameters, obtaining $M_{NFW}=1.0^{+0.4}_{-0.2}\times10^{14} M_\odot$ and $C=3.5_{-0.3}^{+0.5}$. This concentration value is rather low and it would indicate that the cluster halos around these QSOs are unrelaxed. However, higher concentration values still provides a compatible fit to the data.

We present a set of new analytic solutions aimed at self-consistently describing the spatially-averaged time evolution of the gas, stellar, metal, and dust content in an individual starforming galaxy hosted within a dark halo of given mass and formation redshift. Then, as an application, we show that our solutions, when coupled to specific prescriptions for parameter setting (inspired by in-situ galaxy-black hole coevolution scenarios) and merger rates (based on numerical simulations), can be exploited to reproduce the main statistical relationships followed by early-type galaxies and by their high-redshift starforming progenitors. Our analytic solutions allow to easily disentangle the diverse role of the main physical processes regulating galaxy formation, to quickly explore the related parameter space, and to make transparent predictions on spatially-averaged quantities. As such, our analytic solutions may provide a basis for improving the (subgrid) physical recipes presently implemented in theoretical approaches and numerical simulations, and can offer a benchmark for interpreting and forecasting current and future broadband observations of high-redshift starforming galaxies.

We investigate the merging rates of compact binaries in galaxies, and the related detection rate of gravitational wave (GW) events with AdvLIGO/Virgo and with the Einstein Telescope. To this purpose, we rely on three basic ingredients: (i) the redshift-dependent galaxy statistics provided by the latest determination of the star formation rate functions from UV+far-IR/(sub)millimeter/radio data; (ii) star formation and chemical enrichment histories for individual galaxies, modeled on the basis of observations; (iii) compact remnant mass distribution and prescriptions for merging of compact binaries from stellar evolution simulations. We present results for the intrinsic birthrate of compact remnants, the merging rates of compact binaries, GW detection rates and GW counts, attempting to differentiate the outcomes among BH-BH, NS-NS, and BH-NS mergers, and to estimate their occurrence in disk and spheroidal host galaxies. We compare our approach with the one based on cosmic SFR density and cosmic metallicity, exploited by many literature studies; the merging rates from the two approaches are in agreement within the overall astrophysical uncertainties. We also investigate the effects of galaxy-scale strong gravitational lensing of GW in enhancing the rate of detectable events toward high-redshift. Finally, we discuss the contribution of undetected GW emission from compact binary mergers to the stochastic background.

We investigate a sample of 152 dusty sources at 1:5 < z < 2:5 to understand the connection of enhanced Star-Formation-Rate (SFR) and Black-Hole-Accretion-Rate (BHAR). The sources are Herschel-selected, having stellar masses M*>10^10Msun and SFR (100-1000Msun/yr) elevated(>4?) above the star-forming "main sequence", classifying them as Starbursts (SB). Through a multiwavelength fitting approach (including a dusty torus component), we divided the sample into active SBs (dominated by an AGN emission, SBs-AGN, ? 23% of the sample) and purely star-forming SBs (SBs-SFR). We visually inspected their HST/UV-restframe maps: SBs-SFR are generally irregular and composite systems; ? 50% of SBs-AGN are instead dominated by regular compact morphologies. We then found archival ALMA continuum counterparts for 33 galaxies (12 SBs-AGN and 21 SBs-SFR). For these sources we computed dust masses, and, with standard assumptions, we also guessed total molecular gas-masses. SBs turn to be gas rich systems (fgas = Mgas=Mgas/(Mgas+M*) \sim 20%-70%), and the gas fractions of the two SB classes are very similar (fgas = 43 +/-4% and fgas = 42+/- 2%). Our results show that SBs are consistent with a mixture of: 1) highly star-forming merging systems (dominating the SBs-SFR), and 2) primordial galaxies, rapidly growing their M* together with their Black Hole (mainly the more compact SBs-AGN). Anyway, feedback effects have not reduced their fgas yet. Indeed, SBs at z = 2, with relatively low bolometric AGN luminosities in the range 10^44 < Lbol(AGN) < 10^46 erg/s (compared to bright optical and X-ray quasars), are still relatively far from the epoch when the AGN feedback will quench the SFR in the host and will substantially depress the gas fractions.

We present the mean spectral energy distribution (SED) of a sample of optically selected quasars (QSOs) at redshifts of $1 \le z \le 5$. To derive it, we exploit photometric information from SDSS, UKIDSS, and WISE surveys in combination with a stacking analysis of \textit{Herschel}, \textit{AKARI}, and \textit{Planck} maps at the location of the QSOs. The near-UV and optical parts of the reconstructed mean rest-frame SED are similar to those found in other studies. However, the SED shows an excess at 1-2 $\mu$m (when compared to the aforementioned SEDs normalized in the near-UV) and a prominent bump around 4-6 $\mu$m, followed by a decrease out to $\sim 20 \,\mu$m and a subsequent far-IR bump. From the fitted SEDs we estimate the average active galactic nuclei (AGN) luminosity $L_{\rm AGN}$ and star formation rate (SFR) as function of cosmic time, finding typical $L_{\rm AGN} \sim 10^{46} - 10^{47}$ erg/s and SFR $\sim 50 - 1000\, M_{\odot}/$yr. We develop mid-IR based criteria to split the QSO sample, finding that these allow us to move along the average relationship in the SFR vs. $L_{\rm AGN}$ diagram toward increasing AGN luminosities. When interpreted in the context of the in-situ coevolution scenario presented by Lapi et al. 2014, our results suggest that the detection in the far-IR band is an effective criterion to select objects where the star formation is on the verge of being affected by energy/momentum feedback from the central AGN.

We study the signal from patchy reionization in view of the future high accuracy polarization measurements of the Cosmic Microwave Background (CMB). We implement an extraction procedure of the patchy reionization signal analogous to CMB lensing. We evaluate the signal to noise ratio (SNR) for the future Stage IV (S4) CMB experiment. The signal has a broad peak centered on the degree angular scales, with a long tail at higher multipoles. The CMB S4 experiment can effectively constrain the properties of reionization by measuring the signal on degree scales. The signal amplitude depends on the properties of the structure determining the reionization morphology. We describe bubbles having radii distributed log-normally. The expected S/N is sensitive to the mean bubble radius: $\bar{R}=5$ Mpc implies $S/N \approx 4$, $\bar{R}=10$ Mpc implies $S/N \approx 20$. The spread of the radii distribution strongly affects the integrated SNR, that changes by a factor of $10^2$ when $\sigma_{lnr}$ goes from $\ln 2$ to $\ln3$. Future CMB experiments will thus place important constraints on the physics of reionization.

We build templates of rotation curves as a function of the $I-$band luminosity via the mass modeling (by the sum of a thin exponential disk and a cored halo profile) of suitably normalized, stacked data from wide samples of local spiral galaxies. We then exploit such templates to determine fundamental stellar and halo properties for a sample of about $550$ local disk-dominated galaxies with high-quality measurements of the optical radius $R_{\rm opt}$ and of the corresponding rotation velocity $V_{\rm opt}$. Specifically, we determine the stellar $M_\star$ and halo $M_{\rm H}$ masses, the halo size $R_{\rm H}$ and velocity scale ${V_{\rm H}}$, and the specific angular momenta of the stellar $j_\star$ and dark matter $j_{\rm H}$ components. We derive global scaling relationships involving such stellar and halo properties both for the individual galaxies in our sample and for their mean within bins; the latter are found to be in pleasing agreement with previous determinations by independent methods (e.g., abundance matching techniques, weak lensing observations, and individual rotation curve modeling). Remarkably, the size of our sample and the robustness of our statistical approach allow us to attain an unprecedented level of precision over an extended range of mass and velocity scales, with $1\sigma$ dispersion around the mean relationships of less than $0.1$ dex. We thus set new standard local relationships that must be reproduced by detailed physical models, that offer a basis for improving the sub-grid recipes in numerical simulations, that provide a benchmark to gauge independent observations and check for systematics, and that constitute a basic step toward the future exploitation of the spiral galaxy population as a cosmological probe.

We use a 200 $h^{-1}Mpc$ a side N-body simulation to study the mass accretion history (MAH) of dark matter halos to be accreted by larger halos, which we call infall halos. We define a quantity $a_{\rm nf}\equiv (1+z_{\rm f})/(1+z_{\rm peak})$ to characterize the MAH of infall halos, where $z_{\rm peak}$ and $z_{\rm f}$ are the accretion and formation redshifts, respectively. We find that, at given $z_{\rm peak}$, their MAH is bimodal. Infall halos are dominated by a young population at high redshift and by an old population at low redshift. For the young population, the $a_{\rm nf}$ distribution is narrow and peaks at about $1.2$, independent of $z_{\rm peak}$, while for the old population, the peak position and width of the $a_{\rm nf}$ distribution both increases with decreasing $z_{\rm peak}$ and are both larger than those of the young population. This bimodal distribution is found to be closely connected to the two phases in the MAHs of halos. While members of the young population are still in the fast accretion phase at $z_{\rm peak}$, those of the old population have already entered the slow accretion phase at $z_{\rm peak}$. This bimodal distribution is not found for the whole halo population, nor is it seen in halo merger trees generated with the extended Press-Schechter formalism. The infall halo population at $z_{\rm peak}$ are, on average, younger than the whole halo population of similar masses identified at the same redshift. We discuss the implications of our findings in connection to the bimodal color distribution of observed galaxies and to the link between central and satellite galaxies.

[ABRIDGED] We aim to provide a holistic view on the typical size and kinematic evolution of massive early-type galaxies (ETGs), that encompasses their high-$z$ star-forming progenitors, their high-$z$ quiescent counterparts, and their configurations in the local Universe. Our investigation covers the main processes playing a relevant role in the cosmic evolution of ETGs. Specifically, their early fast evolution comprises: biased collapse of the low angular momentum gaseous baryons located in the inner regions of the host dark matter halo; cooling, fragmentation, and infall of the gas down to the radius set by the centrifugal barrier; further rapid compaction via clump/gas migration toward the galaxy center, where strong heavily dust-enshrouded star-formation takes place and most of the stellar mass is accumulated; ejection of substantial gas amount from the inner regions by feedback processes, which causes a dramatic puffing up of the stellar component. In the late slow evolution, passive aging of stellar populations and mass additions by dry merger events occur. We describe these processes relying on prescriptions inspired by basic physical arguments and by numerical simulations, to derive new analytical estimates of the relevant sizes, timescales, and kinematic properties for individual galaxies along their evolution. Then we obtain quantitative results as a function of galaxy mass and redshift, and compare them to recent observational constraints on half-light size $R_e$, on the ratio $v/\sigma$ between rotation velocity and velocity dispersion (for gas and stars) and on the specific angular momentum $j_\star$ of the stellar component; we find good consistency with the available multi-band data in average values and dispersion, both for local ETGs and for their $z\sim 1-2$ star-forming and quiescent progenitors.

We combine the latest datasets obtained with different surveys to study the frequency dependence of polarized emission coming from Extragalactic Radio Sources (ERS). We consider data over a very wide frequency range starting from $1.4$ GHz up to $217$ GHz. This range is particularly interesting since it overlaps the frequencies of the current and forthcoming Cosmic Microwave Background (\cmb) experiments. Current data suggest that at high radio frequencies, ($ \nu \geq 20$ GHz) the fractional polarization of ERS does not depend on the total flux density. Conversely, recent datasets indicate a moderate increase of polarization fraction as a function of frequency, physically motivated by the fact that Faraday depolarization is expected to be less relevant at high radio-frequencies. We compute ERS number counts using updated models based on recent data, and we forecast the contribution of unresolved ERS in CMB polarization spectra. Given the expected sensitivities and the observational patch sizes of forthcoming \cmb\ experiments about $\sim 200 $ ( up to $\sim 2000 $ ) polarized ERS are expected to be detected. Finally, we assess that polarized ERS can contaminate the cosmological B-mode polarization if the tensor-to-scalar ratio is $r< 0.05$ and they have to be robustly controlled to de-lens \cmb\ B-modes at the arcminute angular scales.

Clumping and turbulence are expected to affect the matter accreted onto the outskirts of galaxy clusters. To determine their impact on the thermodynamic properties of Abell 2142 we perform an analysis of the X-ray temperature data from XMM-Newton via our SuperModel, a state-of-the-art tool for investigating the astrophysics of the intracluster medium already tested on many individual clusters (since Cavaliere et al. 2009). Using the gas density profile corrected for clumpiness derived by Tchernin et al. (2016), we find evidence for the presence of a nonthermal pressure component required to sustain gravity in the cluster outskirts of Abell 2142, that amounts to about 30% of the total pressure at the virial radius. The presence of the nonthermal component implies the gas fraction to be consistent with the universal value at the virial radius and the electron thermal pressure profile to be in good agreement with that inferred from the SZ data. Our results indicate that the presence of gas clumping and of a nonthermal pressure component are both necessary to recover the observed physical properties in the cluster outskirts. Moreover, we stress that an alternative method often exploited in the literature (included Abell 2142) to determine the temperature profile k_BT = P_e/n_e basing on a combination of the Sunyaev-Zel'dovich (SZ) pressure P_e and of the X-ray electron density n_e does not allow to highlight the presence of nonthermal pressure support in the cluster outskirts.

In this work we measure and study the cross-correlation signal between a foreground sample of GAMA galaxies with spectroscopic redshifts in the range $0.2<z<0.8$, and a background sample of H-ATLAS galaxies with photometric redshifts $\gtrsim1.2$. It constitutes a substantial improvement over the cross-correlation measurements made by Gonzalez-Nuevo et al. (2014) with updated catalogues and wider area (with $S/N\gtrsim 5$ below 10' and reaching $S/N\sim 20$ below 30"). The better statistics allow us to split the sample in different redshift bins and to perform a tomographic analysis (with $S/N\gtrsim 3$ below 10 arcmin and reaching $S/N\sim 15$ below 30"). Moreover, we implement a halo model to extract astrophysical information about the background galaxies and the deflectors that are producing the lensing link between the foreground (lenses) and background (sources) samples. In the case of the sources, we find typical mass values in agreement with previous studies: a minimum halo mass to host a central galaxy, $M_{min}\sim 10^{12.26} M_\odot$, and a pivot halo mass to have at least one sub-halo satellite, $M_1\sim 10^{12.84} M_\odot$. However, the lenses are massive galaxies or even galaxy groups/clusters, with minimum mass of $M_{min}^{lens}\sim 10^{13.06} M_\odot$. Above a mass of $M_1^{lens}\sim 10^{14.57} M_\odot$ they contain at least one additional satellite galaxy which contributes to the lensing effect. The tomographic analysis shows that, while $M_1^{lens}$ is almost redshift independent, there is a clear evolution of increase $M_{min}^{lens}$ with redshift in agreement with theoretical estimations. Finally, the halo modeling allows us to identify a strong lensing contribution to the cross-correlation for angular scales below 30". This interpretation is supported by the results of basic but effective simulations.

Nuclear activity and star formation play relevant roles in the early stages of galaxy formation. We aim at identifying them in high redshift galaxies by exploiting high-resolution and sensitivity X-ray and mm data to confirm their presence and relative role in contributing to the galaxy SEDs and energy budget. We present the data, model and analysis in the X-ray and mm bands for two strongly lensed galaxies, SDP.9 and SDP.11, selected in the Herschel-ATLAS catalogues as having an excess emission in the mid-IR regime at z>1.5, suggesting nuclear activity in the early stages of galaxy formation. We observed both of them in X-ray with Chandra and analyzed the high-resolution mm data available in the ALMA Science Archive for SDP9, and, by combining the information available, we reconstructed the source morphology. Both the targets were detected in the X-ray, strongly indicating the presence of highly obscured nuclear activity. High resolution ALMA observations for SDP9 in continuum and CO(6-5) spectral line allowed us to estimate the lensed galaxy redshift to a better accuracy than pre-ALMA estimates and to model the emission of the optical, mm, and X-ray band emission for this galaxy. We demonstrated that the X-ray emission is generated in the nuclear environment and it strongly support the presence of nuclear activity in this object. Hence, we identified weak nuclear activity associated with high-z galaxies with large star formation rates, useful to extend the investigation of the relationship between star formation and nuclear activity to two intrinsically less luminous, high-z star forming galaxies than was possible so far. Given our results only for two objects, they solely cannot constrain the evolutionary models, but provide us with interesting hints and set an observational path towards addressing the role of star formation and nuclear activity in forming galaxies.

We discuss the effects on the CMB, CIB, and thermal SZ effect due to the peculiar motion of an observer with respect to the CMB rest frame, which induces boosting effects. We investigate the scientific perspectives opened by future CMB space missions, focussing on the CORE proposal. The improvements in sensitivity offered by a mission like CORE, together with its high resolution over a wide frequency range, will provide a more accurate estimate of the CMB dipole. The extension of boosting effects to polarization and cross-correlations will enable a more robust determination of purely velocity-driven effects that are not degenerate with the intrinsic CMB dipole, allowing us to achieve a S/N ratio of 13; this improves on the Planck detection and essentially equals that of an ideal cosmic-variance-limited experiment up to a multipole l of 2000. Precise inter-frequency calibration will offer the opportunity to constrain or even detect CMB spectral distortions, particularly from the cosmological reionization, because of the frequency dependence of the dipole spectrum, without resorting to precise absolute calibration. The expected improvement with respect to COBE-FIRAS in the recovery of distortion parameters (in principle, a factor of several hundred for an ideal experiment with the CORE configuration) ranges from a factor of several up to about 50, depending on the quality of foreground removal and relative calibration. Even for 1% accuracy in both foreground removal and relative calibration at an angular scale of 1 deg, we find that dipole analyses for a mission like CORE will be able to improve the recovery of the CIB spectrum amplitude by a factor of 17 in comparison with current results based on FIRAS. In addition to the scientific potential of a mission like CORE for these analyses, synergies with other planned and ongoing projects are also discussed.

The continuity equation is developed for the stellar mass content of galaxies, and exploited to derive the stellar mass function of active and quiescent galaxies over the redshift range $z\sim 0-8$. The continuity equation requires two specific inputs gauged on observations: (i) the star formation rate functions determined on the basis of the latest UV+far-IR/sub-mm/radio measurements; (ii) average star-formation histories for individual galaxies, with different prescriptions for discs and spheroids. The continuity equation also includes a source term taking into account (dry) mergers, based on recent numerical simulations and consistent with observations. The stellar mass function derived from the continuity equation is coupled with the halo mass function and with the SFR functions to derive the star formation efficiency and the main sequence of star-forming galaxies via the abundance matching technique. A remarkable agreement of the resulting stellar mass function for active and quiescent galaxies, of the galaxy main sequence and of the star-formation efficiency with current observations is found; the comparison with data also allows to robustly constrain the characteristic timescales for star formation and quiescence of massive galaxies, the star formation history of their progenitors, and the amount of stellar mass added by in-situ star formation vs. that contributed by external merger events. The continuity equation is shown to yield quantitative outcomes that must be complied by detailed physical models, that can provide a basis to improve the (sub-grid) physical recipes implemented in theoretical approaches and numerical simulations, and that can offer a benchmark for forecasts on future observations with multi-band coverage, as it will become routinely achievable in the era of JWST.

Follow-up observations at high-angular resolution of submillimeter galaxies showed that the single-dish sources are comprised of a blend of several galaxies. Consequently, number counts derived from low and high angular resolution observations are in disagreement. This demonstrates the importance of resolution effects and the need to have realistic simulations to explore them. We built a new 2deg^2 simulation of the extragalactic sky from the far-infrared to the submillimeter. It is based on an updated version of the two star-formation mode galaxy evolution model. Using global galaxy properties, we use the abundance matching technique to populate a dark-matter lightcone and thus simulate the clustering. We produce maps and extract the sources, and show that the limited angular resolution of single-dish instruments have a strong impact on (sub)millimeter continuum observations. Taking into account these resolution effects, we are reproducing a large set of observables, including number counts, redshift distributions, and cosmic infrared background power spectra. Our simulation describes consistently the number counts from single-dish telescopes and interferometers. In particular, at 350 and 500 um, we find that number counts measured by Herschel between 5 and 50 mJy are biased towards high values by a factor 2, and that redshift distributions are biased towards low z. We also show that the clustering has an important impact on the Herschel pixel histogram used to derive number counts from P(D) analysis. Finally, we demonstrate that the large number density of red Herschel sources found in observations but not in models could be an observational artifact caused by the combination of noise, resolution effects, and the steepness of color and flux density distributions. Our simulation (SIDES) is available at http://cesam.lam.fr/sides

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

We investigate the origin, the shape, the scatter, and the cosmic evolution in the observed relationship between specific angular momentum $j_\star$ and the stellar mass $M_\star$ in early-type (ETGs) and late-type galaxies (LTGs). Specifically, we exploit the observed star-formation efficiency and chemical abundance to infer the fraction $f_{\rm inf}$ of baryons that infall toward the central regions of galaxies where star formation can occur. We find $f_{\rm inf}\approx 1$ for LTGs and $\approx 0.4$ for ETGs with an uncertainty of about $0.25$ dex, consistent with a biased collapse. By comparing with the locally observed $j_\star$ vs. $M_\star$ relations for LTGs and ETGs we estimate the fraction $f_j$ of the initial specific angular momentum associated to the infalling gas that is retained in the stellar component: for LTGs we find $f_j\approx 1.11^{+0.75}_{-0.44}$, in line with the classic disc formation picture; for ETGs we infer $f_j\approx 0.64^{+0.20}_{-0.16}$, that can be traced back to a $z<1$ evolution via dry mergers. We also show that the observed scatter in the $j_{\star}$ vs. $M_{\star}$ relation for both galaxy types is mainly contributed by the intrinsic dispersion in the spin parameters of the host dark matter halo. The biased collapse plus mergers scenario implies that the specific angular momentum in the stellar components of ETG progenitors at $z\sim 2$ is already close to the local values, in pleasing agreement with observations. All in all, we argue such a behavior to be imprinted by nature and not nurtured substantially by the environment.

Observational investigations of the abundance of massive precursors of local galaxy clusters ("proto-clusters") allow us to test the growth of density perturbations, to constrain cosmological parameters that control it, to test the theory of non-linear collapse and how the galaxy formation takes place in dense environments. The Planck collaboration has recently published a catalogue of >~ 2000 cold extra-galactic sub-millimeter sources, i.e. with colours indicative of z >~ 2, almost all of which appear to be over-densities of star-forming galaxies. They are thus considered as proto-cluster candidates. Their number densities (or their flux densities) are far in excess of expectations from the standard scenario for the evolution of large-scale structure. Simulations based on a physically motivated galaxy evolution model show that essentially all cold peaks brighter than S_{545GHz} = 500 mJy found in Planck maps after having removed the Galactic dust emission can be interpreted as positive Poisson fluctuations of the number of high-z dusty proto-clusters within the same Planck beam, rather then being individual clumps of physically bound galaxies. This conclusion does not change if an empirical fit to the luminosity function of dusty galaxies is used instead of the physical model. The simulations accurately reproduce the statistic of the Planck detections and yield distributions of sizes and ellipticities in qualitative agreement with observations. The redshift distribution of the brightest proto-clusters contributing to the cold peaks has a broad maximum at 1.5 <~ z <~ 3. Therefore follow-up of Planck proto-cluster candidates will provide key information on the high-z evolution of large scale structure.

We investigate the astrophysics of radio-emitting star-forming galaxies and ac- tive galactic nuclei (AGNs), and elucidate their statistical properties in the radio band including luminosity functions, redshift distributions, and number counts at sub-mJy flux levels, that will be crucially probed by next-generation radio continuum surveys. Specifically, we exploit the model-independent approach by Mancuso et al. (2016a,b) to compute the star formation rate functions, the AGN duty cycles and the conditional probability of a star-forming galaxy to host an AGN with given bolometric luminosity. Coupling these ingredients with the radio emission properties associated to star formation and nuclear activity, we compute relevant statistics at different radio frequencies, and disentangle the relative con- tribution of star-forming galaxies and AGNs in different radio luminosity, radio flux, and redshift ranges. Finally, we highlight that radio-emitting star-forming galaxies and AGNs are expected to host supermassive black holes accreting with different Eddington ratio distributions, and to occupy different loci in the galaxy main sequence diagrams. These specific predictions are consistent with current datasets, but need to be tested with larger statistics via future radio data with multi-band coverage on wide areas, as it will become routinely achievable with the advent of the SKA and its precursors.

The assessment of the relationship between radio continuum luminosity and star formation rate (SFR) is of crucial importance to make reliable predictions for the forthcoming ultra-deep radio surveys and to allow a full exploitation of their results to measure the cosmic star formation history. We have addressed this issue by matching recent accurate determinations of the SFR function up to high redshifts with literature estimates of the 1.4 GHz luminosity functions of star forming galaxies (SFGs). This was done considering two options, proposed in the literature, for the relationship between the synchrotron emission ($L_{\rm synch}$), that dominates at 1.4 GHz, and the SFR: a linear relation with a decline of the $L_{\rm synch}$/SFR ratio at low luminosities or a mildly non-linear relation at all luminosities. In both cases we get good agreement with the observed radio luminosity functions but, in the non-linear case, the deviation from linearity must be small. The luminosity function data are consistent with a moderate increase of the $L_{\rm synch}$/SFR ratio with increasing redshift, indicated by other data sets, although a constant ratio cannot be ruled out. A stronger indication of such increase is provided by recent deep 1.4 GHz counts, down to $\mu$Jy levels. This is in contradiction with models predicting a decrease of that ratio due to inverse Compton cooling of relativistic electrons at high redshifts. Synchrotron losses appear to dominate up to $z\simeq 5$. We have also updated the Massardi et al. (2010) evolutionary model for radio loud AGNs.

By means of the updated PARSEC database of evolutionary tracks of massive stars, we compute the integrated stellar light, the ionizing photon budget and the supernova rates of young simple stellar populations (SSPs), for different metallicities and IMF upper mass limits. Using CLOUDY we compute and include in the SSP spectra the neb- ular emission contribution. We also revisit the thermal and non-thermal radio emission contribution from young stars. Using GRASIL we can thus predict the panchromatic spectrum and the main recombination lines of any type of star-forming galaxy, including the effects of dust absorption and re-emission. We check the new models against the spectral energy distributions (SEDs) of selected well-observed nearby galaxies. From the best-fit models we obtain a consistent set of star formation rate (SFR) calibrations at wavelengths ranging from ultraviolet (UV) to radio. We also provide analytical calibrations that take into account the dependence on metallcity and IMF upper mass limit of the SSPs. We show that the latter limit can be well constrained by combining information from the observed far infrared, 24 {\mu}m, 33 GHz and H{\alpha} luminosities. Another interesting property derived from the fits is that, while in a normal galaxy the attenuation in the lines is significantly higher than that in the nearby continuum, in individual star bursting regions they are similar, supporting the notion that this effect is due to an age selective extinction. Since in these conditions the Balmer decrement method may not be accurate, we provide relations to estimate the attenuation from the observed 24 {\mu}m or 33 GHz fluxes. These relations can be useful for the analysis of young high redshift galaxies.

We compare the set of local galaxies having dynamically measured black holes with a large, unbiased sample of galaxies extracted from the Sloan Digital Sky Survey. We confirm earlier work showing that the majority of black hole hosts have significantly higher velocity dispersions sigma than local galaxies of similar stellar mass. We use Monte-Carlo simulations to illustrate the effect on black hole scaling relations if this bias arises from the requirement that the black hole sphere of influence must be resolved to measure black hole masses with spatially resolved kinematics. We find that this selection effect artificially increases the normalization of the Mbh-sigma relation by a factor of at least ~3; the bias for the Mbh-Mstar relation is even larger. Our Monte Carlo simulations and analysis of the residuals from scaling relations both indicate that sigma is more fundamental than Mstar or effective radius. In particular, the Mbh-Mstar relation is mostly a consequence of the Mbh-sigma and sigma-Mstar relations, and is heavily biased by up to a factor of 50 at small masses. This helps resolve the discrepancy between dynamically-based black hole-galaxy scaling relations versus those of active galaxies. Our simulations also disfavour broad distributions of black hole masses at fixed sigma. Correcting for this bias suggests that the calibration factor used to estimate black hole masses in active galaxies should be reduced to values of fvir~1. Black hole mass densities should also be proportionally smaller, perhaps implying significantly higher radiative efficiencies/black hole spins. Reducing black hole masses also reduces the gravitational wave signal expected from black hole mergers.

We provide an holistic view of galaxy evolution at high redshift z>4, that incorporates the constraints from various astrophysical/cosmological probes, including the estimate of the cosmic SFR density from UV/IR surveys and long GRB rates, the cosmic reionization history after the latest Planck measurements, and the missing satellites issue. We achieve this goal in a model-independent way by exploiting the SFR functions derived by Mancuso et al. (2016) on the basis of an educated extrapolation of the latest UV/far-IR data from HST/Herschel, and already tested against a number of independent observables. Our SFR functions integrated down to an UV magnitude limit M_UV<-13 (or SFR limit around 10^-2 M_sun/yr) produces a cosmic SFR density in excellent agreement with recent determinations from IR surveys and, taking into account a metallicity ceiling Z<Z_sun/2, with the estimates from long GRB rates. They also yield a cosmic reionization history consistent with that implied by the recent measurements of the Planck mission on the electron scattering optical depth tau_es~0.058; remarkably, this result is obtained under a conceivable assumption regarding the average value f_esc~0.1 of the escape fraction for ionizing photons. We demonstrate via the abundance matching technique that the above constraints concurrently imply galaxy formation to become inefficient within dark matter halos of mass below a few 10^8 M_sun; pleasingly, such a limit is also required not to run into the missing satellite issue. Finally, we predict a downturn of the galaxy luminosity function faintward of M_UV<-12, and stress that its detailed shape, as plausibly probed in the next future by the JWST, will be extremely informative on the astrophysics of galaxy formation in small halos, or even on the microscopic nature of the dark matter.

We present a sample of 80 candidate strongly lensed galaxies with flux density above 100mJy at 500{\mu}m extracted from the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS), over an area of 600 square degrees. Available imaging and spectroscopic data allow us to confirm the strong lensing in 20 cases and to reject it in one case. For other 8 objects the lensing scenario is strongly supported by the presence of two sources along the same line of sight with distinct photometric redshifts. The remaining objects await more follow-up observations to confirm their nature. The lenses and the background sources have median redshifts z_L = 0.6 and z_S = 2.5, respectively, and are observed out to z_L = 1.2 and z_S = 4.2. We measure the number counts of candidate lensed galaxies at 500{\mu}m and compare them with theoretical predictions, finding a good agreement for a maximum magnification of the background sources in the range 10-20. These values are consistent with the magnification factors derived from the lens modelling of individual systems. The catalogue presented here provides sub- mm bright targets for follow-up observations aimed at exploiting gravitational lensing to study with un-precedented details the morphological and dynamical properties of dusty star forming regions in z >~ 1.5 galaxies.

We provide a novel, unifying physical interpretation on the origin, the average shape, the scatter, and the cosmic evolution for the main sequences of starforming galaxies and active galactic nuclei at high redshift z $\gtrsim$ 1. We achieve this goal in a model-independent way by exploiting: (i) the redshift-dependent SFR functions based on the latest UV/far-IR data from HST/Herschel, and re- lated statistics of strong gravitationally lensed sources; (ii) deterministic evolutionary tracks for the history of star formation and black hole accretion, gauged on a wealth of multiwavelength observations including the observed Eddington ratio distribution. We further validate these ingredients by showing their consistency with the observed galaxy stellar mass functions and AGN bolometric luminosity functions at different redshifts via the continuity equation approach. Our analysis of the main sequence for high-redshift galaxies and AGNs highlights that the present data are consistently interpreted in terms of an in situ coevolution scenario for star formation and black hole accretion, envisaging these as local, time coordinated processes.

We present an improved and extended analysis of the cross-correlation between the map of the Cosmic Microwave Background (CMB) lensing potential derived from the \emph{Planck} mission data and the high-redshift galaxies detected by the \emph{Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS) in the photometric redshift range $z_{\rm ph} \ge 1.5$. We compare the results based on the 2013 and 2015 \textit{Planck} datasets, and investigate the impact of different selections of the H-ATLAS galaxy samples. Significant improvements over our previous analysis have been achieved thanks to the higher signal-to-noise ratio of the new CMB lensing map recently released by the \textit{Planck} collaboration. The effective galaxy bias parameter, $b$, for the full galaxy sample, derived from a joint analysis of the cross-power spectrum and of the galaxy auto-power spectrum is found to be $b = 3.54^{+0.15}_{-0.14}$. Furthermore, a first tomographic analysis of the cross-correlation signal is implemented, by splitting the galaxy sample into two redshift intervals: $1.5 \le z_{\rm ph} < 2.1$ and $z_{\rm ph}\ge 2.1$. A statistically significant signal was found for both bins, indicating a substantial increase with redshift of the bias parameter: $b=2.89\pm0.23$ for the lower and $b=4.75^{+0.24}_{-0.25}$ for the higher redshift bin. Consistently with our previous analysis we find that the amplitude of the cross correlation signal is a factor of $1.45^{+0.14}_{-0.13}$ higher than expected from the standard $\Lambda$CDM model for the assumed redshift distribution. The robustness of our results against possible systematic effects has been extensively discussed although the tension is mitigated by passing from 4 to 3$\sigma$.

Galaxy groups differ from clusters primarily by way of their lower masses, M~10^14 M_sun vs. M~10^15 M_sun. We discuss how mass affects the thermal state of the intracluster or the intragroup medium, specifically as to their entropy levels and radial profiles. We show that entropy is produced in both cases by the continuing inflow of intergalactic gas across the system boundary into the gravitational potential well. The inflow is highly supersonic in clusters, but weakly so in groups. The former condition implies strong accretion shocks with substantial conversion of a large inflow kinetic into thermal energy, whereas the latter condition implies less effective conversion of lower energies. These features produce a conspicuous difference in entropy deposition at the current boundary. Thereafter, adiabatic compression of the hot gas into the potential well converts such time histories into radial profiles throughout a cluster or a group. In addition, in both cases a location of the system at low z in the accelerating universe or in a poor environment will starve out the inflow and the entropy production, and produce flattening or even bending down of the outer profile. We analyze in detail the sharp evidence provided by the two groups ESO 3060170 and RXJ1159+5531 that have been recently observed in X rays out to their virial radii, and find a close and detailed match with our expectations.

We exploit the continuity equation approach and the `main sequence' star-formation timescales to show that the observed high abundance of galaxies with stellar masses > a few 10^10 M_sun at redshift z>4 implies the existence of a galaxy population featuring large star formation rates (SFRs) > 10^2 M_sun/yr in heavily dust-obscured conditions. These galaxies constitute the high-redshift counterparts of the dusty star-forming population already surveyed for z<3 in the far-IR band by the Herschel space observatory. We work out specific predictions for the evolution of the corresponding stellar mass and SFR functions out to z~10, elucidating that the number density at z<8 for SFRs >30 M_sun/yr cannot be estimated relying on the UV luminosity function alone, even when standard corrections for dust extinction based on the UV slope are applied. We compute the number counts and redshift distributions (including galaxy-scale gravitational lensing) of this galaxy population, and show that current data from AzTEC-LABOCA, SCUBA-2 and ALMA-SPT surveys are already digging into it. We substantiate how an observational strategy based on a color preselection in the far-IR or (sub-)mm band with Herschel and SCUBA-2, supplemented by photometric data via on-source observations with ALMA, can allow to reconstruct the bright end of the SFR functions out to z~8. In parallel, such a challenging task can be managed by exploiting current UV surveys in combination with (sub-)mm observations by ALMA and NIKA2 and/or radio observations by SKA and its precursors.

We review aspects of Cosmic Microwave Background (CMB) spectral distortions which do not appear to have been fully explored in the literature. In particular, implications of recent evidences of heating of the intergalactic medium (IGM) by feedback from active galactic nuclei are investigated. Taking also into account the IGM heating associated to structure formation, we argue that values of the y parameter of several*10^(-6), i.e. a factor of a few below the COBE/FIRAS upper limit, are to be expected. The Compton scattering by the re-ionized plasma also re-processes primordial distortions, adding a y-type contribution. Hence no pure Bose-Einstein-like distortions are to be expected. An assessment of Galactic and extragalactic foregrounds, taking into account the latest results from the Planck satellite as well as the contributions from the strong CII and CO lines from star-forming galaxies, demonstrates that a foreground subtraction accurate enough to fully exploit the PIXIE sensitivity will be extremely challenging. Motivated by this fact we also discuss methods to detect spectral distortions not requiring absolute measurements and show that accurate determinations of the frequency spectrum of the CMB dipole amplitude may substantially improve over COBE/FIRAS limits on distortion parameters. Such improvements may be at reach of next generation CMB anisotropy experiments. The estimated amplitude of the Cosmic Infrared Background (CIB) dipole might be detectable by careful analyses of Planck maps at the highest frequencies. Thus Planck might provide interesting constraints on the CIB intensity, currently known with a ~30% uncertainty.

Dark matter constitutes the great majority of the matter content in the Universe, but its microscopic nature remains an intriguing mystery, with profound implications for particle physics, astrophysics and cosmology. Here we shed light on the longstanding issue of whether the dark matter is warm or cold by combining the measurements of the galaxy luminosity functions out to high redshifts z~10 from the Hubble Space Telescope with the recent cosmological data on the reionization history of the Universe from the Planck mission. We derive robust and tight bounds on the mass of warm dark matter particle, finding that the current data require it to be in the narrow range between 2 and 3 keV. In addition, we show that a mass not exceeding 3 keV is also concurrently indicated by astrophysical constraints related to the local number of satellites in Milky Way-sized galaxies, though it is in marginal tension with analysis of the Lyman-alpha forest. For warm dark matter masses above 3 keV as well as for cold dark matter, to satisfy the Planck constraints on the optical depth and not to run into the satellite problem would require invoking astrophysical processes that inhibit galaxy formation in halos with mass M_H< a few 10^8 M_sun, corresponding to a limiting UV magnitude M_UV~-11. Anyway, we predict a downturn of the galaxy luminosity function at z~8 faintward of M_UV~-12, and stress that its detailed shape is extremely informative both on particle physics and on the astrophysics of galaxy formation in small halos. These expectations will be tested via the Hubble Frontier Fields and with the advent of the James Webb Space Telescope, which will enable probing the very faint end of the galaxy luminosity function out to z~8-10.

We have worked out predictions for the radio counts of star-forming galaxies down to nJy levels, along with redshift distributions down to the detection limits of the phase 1 Square Kilometer Array MID telescope (SKA1-MID) and of its precursors. Such predictions were obtained by coupling epoch dependent star formation rate (SFR) functions with relations between SFR and radio (synchrotron and free-free) emission. The SFR functions were derived taking into account both the dust obscured and the unobscured star-formation, by combining far-infrared (FIR), ultra-violet (UV) and H_alpha luminosity functions up to high redshifts. We have also revisited the South Pole Telescope (SPT) counts of dusty galaxies at 95\,GHz performing a detailed analysis of the Spectral Energy Distributions (SEDs). Our results show that the deepest SKA1-MID surveys will detect high-z galaxies with SFRs two orders of magnitude lower compared to Herschel surveys. The highest redshift tails of the distributions at the detection limits of planned SKA1-MID surveys comprise a substantial fraction of strongly lensed galaxies. We predict that a survey down to 0.25 microJy at 1.4 GHz will detect about 1200 strongly lensed galaxies per square degree, at redshifts of up to 10. For about 30% of them the SKA1-MID will detect at least 2 images. The SKA1-MID will thus provide a comprehensive view of the star formation history throughout the re-ionization epoch, unaffected by dust extinction. We have also provided specific predictions for the EMU/ASKAP and MIGHTEE/MeerKAT surveys.

[abridged] We investigate the coevolution of galaxies and hosted supermassive black holes throughout the history of the Universe by a statistical approach based on the continuity equation and the abundance matching technique. Specifically, we present analytical solutions of the continuity equation without source term to reconstruct the supermassive black hole (BH) mass function from the AGN luminosity functions. Such an approach includes physically-motivated AGN lightcurves tested on independent datasets, which describe the evolution of the Eddington ratio and radiative efficiency from slim- to thin-disc conditions. We nicely reproduce the local estimates of the BH mass function, the AGN duty cycle as a function of mass and redshift, along with the Eddington ratio function and the fraction of galaxies with given stellar mass hosting an AGN with given Eddington ratio. We exploit the same approach to reconstruct the observed stellar mass function at different redshift from the UV and far-IR luminosity functions associated to star formation in galaxies. These results imply that the buildup of stars and BHs in galaxies occurs via in-situ processes, with dry mergers playing a marginal role at least for stellar masses < 3 10^11 M_sun and BH masses < 10^9 M_sun, where the statistical data are more secure and less biased by systematic errors. In addition, we develop an improved abundance matching technique to link the stellar and BH content of galaxies to the gravitationally dominant dark matter component. The resulting relationships constitute a testbed for galaxy evolution models, highlighting the complementary role of stellar and AGN feedback in the star formation process. Finally, the clustering properties of BHs and galaxies are found to be in full agreement with current observations, so further validating our results from the continuity equation.

We present predictions for number counts and redshift distributions of galaxies detectable in continuum and in emission lines with the Mid-infrared (MIR) Instrument (SMI) proposed for the Space Infrared Telescope for Cosmology and Astrophysics (SPICA). We have considered 24 MIR fine-structure lines, four Polycyclic Aromatic Hydrocarbon (PAH) bands (at 6.2, 7.7, 8.6 and 11.3$\mu$m) and two silicate bands (in emission and in absorption) at 9.7$\mu$m and 18.0$\mu$m. Six of these lines are primarily associated with Active Galactic Nuclei (AGNs), the others with star formation. A survey with the SMI spectrometers of 1 hour integration per field-of-view (FoV) over an area of $1\,\hbox{deg}^2$ will yield $5\,\sigma$ detections of $\simeq 140$ AGN lines and of $\simeq 5.2\times10^{4}$ star-forming galaxies, $\simeq 1.6\times10^{4}$ of which will be detected in at least two lines. The combination of a shallow ($20.0\,\hbox{deg}^{2}$, $1.4\times10^{-1}$ h integration per FoV) and a deep survey ($6.9\times10^{-3}\,\hbox{deg}^{2}$, $635$ h integration time), with the SMI camera, for a total of $\sim$1000 h, will accurately determine the MIR number counts of galaxies and of AGNs over five orders of magnitude in flux density, reaching values more than one order of magnitude fainter than the deepest Spitzer $24\,\mu$m surveys. This will allow us to determine the cosmic star formation rate (SFR) function down to SFRs more than 100 times fainter than reached by the Herschel Observatory.

We present the first measurement of the correlation between the map of the cosmic microwave background (CMB) lensing potential derived from the \emph{Planck} nominal mission data and $z\gtrsim 1.5$ galaxies detected by the \emph{Herschel}-ATLAS (H-ATLAS) survey covering about $600\,\hbox{deg}^2$, i.e. about 1.4\% of the sky. We reject the hypothesis that there is no correlation between CMB lensing and galaxy detection at a $20\,\sigma$ significance, checking the result by performing a number of null tests. The significance of the detection of the theoretically expected cross-correlation signal is found to be $10\,\sigma$. The galaxy bias parameter, $b$, derived from a joint analysis of the cross-power spectrum and of the auto-power spectrum of the galaxy density contrast is found to be $b=2.80^{+0.12}_{-0.11}$, consistent with earlier estimates for H-ATLAS galaxies at similar redshifts. On the other hand, the amplitude of the cross-correlation is found to be a factor $1.62 \pm 0.16$ higher than expected from the standard model and also found by cross-correlation analyses with other tracers of the large-scale structure. The enhancement due to lensing magnification can account for only a fraction of the excess cross-correlation signal. We suggest that part of it may be due to an incomplete removal of the contamination of the CIB, that includes the H-ATLAS sources we are cross-correlating with. In any case, the highly significant detection reported here using a catalog covering only 1.4\% of the sky demonstrates the potential of CMB lensing correlations with submillimeter surveys.

We investigate the impact that warm dark matter (WDM) has in terms of 21cm intensity mapping in the post-reionization Universe at z = 3 - 5. We perform hydrodynamic simulations for 5 different models: cold dark matter and WDM with 1,2,3,4 keV (thermal relic) mass and assign the neutral hydrogen a-posteriori using two different methods that both reproduce observations in terms of column density distribution function of neutral hydrogen systems. Contrary to naive expectations, the suppression of power present in the linear and non-linear matter power spectra, results in an increase of power in terms of neutral hydrogen and 21cm power spectra. This is due to the fact that there is a lack of small mass halos in WDM models with respect to cold dark matter: in order to distribute a total amount of neutral hydrogen within the two cosmological models, a larger quantity has to be placed in the most massive halos, that are more biased compared to the cold dark matter cosmology. We quantify this effect and address significance for the telescope SKA1-LOW, including a realistic noise modeling. The results indicate that we will be able to rule out a 4 keV WDM model with 5000 hours of observations at z > 3, with a statistical significance of > 3 sigma, while a smaller mass of 3 keV, comparable to present day constraints, can be ruled out at more than 2 sigma confidence level with 1000 hours of observations at z > 5.

We have combined determinations of the epoch-dependent star formation rate (SFR) function with relationships between SFR and radio (synchrotron and free-free) emission to work out detailed predictions for the counts and the redshift distributions of star-forming galaxies detected by planned Square Kilometer Array (SKA) surveys. The evolving SFR function comes from recent models fitting the far-infrared (FIR) to millimeter-wave luminosity functions and the ultraviolet (UV) luminosity functions up to z=10, extended to take into account additional UV survey data. We used very deep 1.4 GHz number counts from the literature to check the relationship between SFR and synchrotron emission, and the 95 GHz South Pole Telescope (SPT) counts of dusty galaxies to test the relationship between SFR and free-free emission. We show that the SKA will allow us to investigate the SFRs of galaxies down to few Msun/yr up to z=10, thus extending by more than two orders of magnitude the high-z SFR functions derived from Herschel surveys. SKA1-MID surveys, down to microJy levels, will detect hundreds of strongly lensed galaxies per square degree; a substantial fraction of them will show at least two images above the detection limits.

We present an analysis of high-quality X-ray data out to the virial radius for the two galaxy clusters Abell 1246 and GMBCG J255.34805+64.23661 (J255) by means of our entropy-based SuperModel. For Abell 1246 we find that the spherically-averaged entropy profile of the intracluster medium (ICM) progressively flattens outwards, and that a nonthermal pressure component amounting to ~20% of the total is required to support hydrostatic equilibrium in the outskirts; there we also estimate a modest value C~1.6 of the ICM clumping factor. These findings agree with previous analyses on other cool-core, relaxed clusters, and lend further support to the picture by Lapi et al. (2010) that relates the entropy flattening, the development of nonthermal pressure component, and the azimuthal variation of ICM properties to weakening boundary shocks. In this scenario clusters are born in a high-entropy state throughout, and are expected to develop on similar timescales a low entropy state both at the center due to cooling, and in the outskirts due to weakening shocks. However, the analysis of J255 testifies how such a typical evolutionary course can be interrupted or even reversed by merging especially at intermediate redshift, as predicted by Cavaliere et al. (2011b). In fact, a merger has rejuvenated the ICM of this cluster at z~0.45 by reestablishing a high entropy state in the outskirts, while leaving intact or erasing only partially the low-entropy, cool core at the center.

In galaxy clusters the thermal Sunyaev-Zel'dovich (SZ) effect from the hot intracluster medium (ICM) provides a direct, self-contained measure of the pressure integrated over crossing lines of sight, that is intrinsically independent of redshift and well suited for evolutionary studies. We show in detail how the size of the effect and its pattern on the sky plane are directly related to the entropy levels in the ICM, and how they characterize the cluster cores and outskirts independently. We find that at redshifts z<0.3 the signals to be expected in the cores considerably exceed those detected at 10' resolution with the Planck satellite. We propose that at 1' resolutions as implemented on recent ground instrumentation for mapping features in individual clusters, the average patterns of the SZ signals can provide a direct and effective way to find and count cool, low-entropy cores and hot, high-entropy outskirts out to z~2. Such counts will tell the timing and the mode of the processes that drive the evolution of the ICM from the distant to the local cluster population.

We present new estimates of redshift-dependent luminosity functions of IR lines detectable by SPICA/SAFARI and excited both by star formation and by AGN activity. The new estimates improve over previous work by using updated evolutionary models and dealing in a self consistent way with emission of galaxies as a whole, including both the starburst and the AGN component. New relationships between line and AGN bolometric luminosity have been derived and those between line and IR luminosities of the starburst component have been updated. These ingredients were used to work out predictions for the source counts in 11 mid/far-IR emission lines partially or entirely excited by AGN activity. We find that the statistics of the emission line detection of galaxies as a whole is mainly determined by the star formation rate, because of the rarity of bright AGNs. We also find that the slope of the line integral number counts is flatter than 2 implying that the number of detections at fixed observing time increases more by extending the survey area than by going deeper. We thus propose a wide spectroscopic survey of 1 hour integration per field-of-view over an area of 5 deg$^{2}$ to detect (at 5$\sigma$) $\sim$760 AGNs in [OIV]25.89$\mu$m $-$ the brightest AGN mid-infrared line $-$ out to $z\sim2$. Pointed observations of strongly lensed or hyper-luminous galaxies previously detected by large area surveys such as those by Herschel and by the SPT can provide key information on the galaxy-AGN co-evolution out to higher redshifts.

We derive approximated, yet very accurate analytical expressions for the abundance and clustering properties of dark matter halos in the excursion set peak framework; the latter relies on the standard excursion set approach, but also includes the effects of a realistic filtering of the density field, a mass-dependent threshold for collapse, and the prescription from peak theory that halos tend to form around density maxima. We find that our approximations work excellently for diverse power spectra, collapse thresholds and density filters. Moreover, when adopting a cold dark matter power spectra, a tophat filtering and a mass-dependent collapse threshold (supplemented with conceivable scatter), our approximated halo mass function and halo bias represent very well the outcomes of cosmological $N-$body simulations.

We report a highly significant ($>10\sigma$) spatial correlation between galaxies with $S_{350\mu\rm m}\ge 30\,$mJy detected in the equatorial fields of the \textsl{Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS) with estimated redshifts $\gtrsim 1.5$, and SDSS or GAMA galaxies at $0.2\le z\le 0.6$. The significance of the cross-correlation is much higher than those reported so far for samples with non-overlapping redshift distributions selected in other wavebands. Extensive, realistic simulations of clustered sub-mm galaxies amplified by foreground structures confirm that the cross-correlation is explained by weak gravitational lensing ($\mu<2$). The simulations also show that the measured amplitude and range of angular scales of the signal are larger than can be accounted for by galaxy-galaxy weak lensing. However, for scales $\lesssim 2\,$arcmin, the signal can be reproduced if SDSS/GAMA galaxies act as signposts of galaxy groups/clusters with halo masses in the range $10^{13.2}$--$10^{14.5} M_{\odot}$. The signal detected on larger scales appears to reflect the clustering of such halos.

[Abridged] We present a physical model for the evolution of the ultraviolet (UV) luminosity function (LF) of high-z galaxies taking into account in a self-consistent way their chemical evolution and the associated evolution of dust extinction. The model yields good fits of the UV and Lyman-alpha LFs at z>~2. The weak evolution of both LFs between z=2 and z=6 is explained as the combined effect of the negative evolution of the halo mass function, of the increase with redshift of the star formation efficiency, and of dust extinction. The slope of the faint end of the UV LF is found to steepen with increasing redshift, implying that low luminosity galaxies increasingly dominate the contribution to the UV background at higher and higher redshifts. The observed range of UV luminosities at high-z implies a minimum halo mass capable of hosting active star formation M_crit <~ 10^9.8 M_odot, consistent with the constraints from hydrodynamical simulations. From fits of Lyman-alpha LFs plus data on the luminosity dependence of extinction and from the measured ratios of non-ionizing UV to Lyman-continuum flux density for samples of z=~3 Lyman break galaxies and Lyman-alpha emitters, we derive a simple relationship between the escape fraction of ionizing photons and the star formation rate, impling larger escape fraction for less massive galaxies. Galaxies already represented in the UV LF (M_UV <~ -18) can keep the universe fully ionized up to z=~6, consistent with (uncertain) data pointing to a rapid drop of the ionization degree above z~6. On the other side, the electron scattering optical depth, tau_es, inferred from CMB experiments favor an ionization degree close to unity up to z=~9-10. Consistency with CMB data can be achieved if M_crit =~ 10^8.5 M_odot, implying that the UV LFs extend to M_UV =~ -13, although the corresponding tau_es is still on the low side of CMB-based estimates.

Several physical processes and formation events are expected in cluster outskirts, a vast region up to now essentially not covered by observations. The recent Suzaku (X-ray) and Planck (Sunayev-Zeldovich effect) observations out to the virial radius have highlighted in these peripheral regions a rather sharp decline of the intracluster gas temperature, an entropy flattening in contrast with the theoretically expected power law increase, the break of the hydrostatic equilibrium even in some relaxed clusters, a derived gas mass fraction above the cosmic value measured from several CMB experiments, and a total X-ray mass lower than the weak lensing mass determinations. Here we present the analysis of four clusters (A1795, A2029, A2204 and A133) with the SuperModel that includes a nonthermal pressure component due to turbulence to sustain the hydrostatic equilibrium also in the cluster outskirts. In such way we obtain a correct determination of the total X-ray mass and of the gas mass fraction; this in turn allows to determine the level of the gas clumping that can affect the shape of the entropy profiles reported by the Suzaku observations. Our conclusion is that the role of the gas clumping is very marginal and that the observed entropy flattening is due to the rapid decrement of the temperature in the cluster outskirts caused by non gravitational effects. Moreover, we show that the X-ray/SZ joint analysis from ROSAT and Planck data, as performed in some recent investigations, is inadequate to discriminate between a power law increase and a flattening of the entropy.

We exploit the recent, wide samples of far-infrared (FIR) selected galaxies followed-up in X rays and of X-ray/optically selected active galactic nuclei (AGNs) followed-up in the FIR band, along with the classic data on AGN and stellar luminosity functions at high redshift z>1.5, to probe different stages in the coevolution of supermassive black holes (BHs) and host galaxies. The results of our analysis indicate the following scenario: (i) the star formation in the host galaxy proceeds within a heavily dust-enshrouded medium at an almost constant rate over a timescale ~0.5-1 Gyr, and then abruptly declines due to quasar feedback; over the same timescale, (ii) part of the interstellar medium loses angular momentum, reaches the circum-nuclear regions at a rate proportional to the star formation and is temporarily stored into a massive reservoir/proto-torus wherefrom it can be promptly accreted; (iii) the BH grows by accretion in a self-regulated regime with radiative power that can slightly exceed the Eddington limit L/L_Edd< 4, particularly at the highest redshifts; (iv) for massive BHs the ensuing energy feedback at its maximum exceeds the stellar one and removes the interstellar gas, thus stopping the star formation and the fueling of the reservoir; (v) afterwards, if the latter has retained enough gas, a phase of supply-limited accretion follows exponentially declining with a timescale of about 2 e-folding times. We show that the ratio of the FIR luminosity of the host galaxy to the bolometric luminosity of the AGN maps the various stages of the above sequence. Finally, we discuss how the detailed properties and the specific evolution of the reservoir can be investigated via coordinated, high-resolution observations of starforming, strongly-lensed galaxies in the (sub-)mm band with ALMA and in the X-ray band with Chandra and the next generation X-ray instruments.

We present the first set of maps and band-merged catalog from the Herschel Stripe 82 Survey (HerS). Observations at 250, 350, and 500 micron were taken with the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. HerS covers 79 deg$^2$ along the SDSS Stripe 82 to a depth of 13.0, 12.9, and 14.8 mJy beam$^{-1}$ (including confusion) at 250, 350, and 500 micron, respectively. HerS was designed to measure correlations with external tracers of the dark matter density field --- either point-like (i.e., galaxies selected from radio to X-ray) or extended (i.e., clusters and gravitational lensing) --- in order to measure the bias and redshift distribution of intensities of infrared-emitting dusty star-forming galaxies and AGN. By locating HeRS in Stripe 82, we maximize the overlap with available and upcoming cosmological surveys. The band-merged catalog contains 3.3x10$^4$ sources detected at a significance of >3 $\sigma$ (including confusion noise). The maps and catalog are available at http://www.astro.caltech.edu/hers/

While continuum imaging data at far-infrared to sub-millimeter wavelengths have provided tight constraints on the population properties of dusty star forming galaxies up to high redshifts, future space missions like the Space Infra-Red Telescope for Cosmology and Astrophysics (SPICA) and ground based facilities like the Cerro Chajnantor Atacama Telescope (CCAT) will allow detailed investigations of their physical properties via their mid-/far-infrared line emission. We present updated predictions for the number counts and the redshift distributions of star forming galaxies spectroscopically detectable by these future missions. These predictions exploit a recent upgrade of evolutionary models, that include the effect of strong gravitational lensing, in the light of the most recent Herschel and South Pole Telescope data. Moreover the relations between line and continuum infrared luminosity are re-assessed, considering also differences among source populations, with the support of extensive simulations that take into account dust obscuration. The derived line luminosity functions are found to be highly sensitive to the spread of the line to continuum luminosity ratios. Estimates of the expected numbers of detections per spectral line by SPICA/SAFARI and by CCAT surveys for different integration times per field of view at fixed total observing time are presented. Comparing with the earlier estimates by Spinoglio et al. (2012) we find, in the case of SPICA/SAFARI, differences within a factor of two in most cases, but occasionally much larger. More substantial differences are found for CCAT.

We report on deep near-infrared observations obtained with the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST) of the first five confirmed gravitational lensing events discovered by the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS). We succeed in disentangling the background galaxy from the lens to gain separate photometry of the two components. The HST data allow us to significantly improve on previous constraints of the mass in stars of the lensed galaxy and to perform accurate lens modelling of these systems, as described in the accompanying paper by Dye et al. We fit the spectral energy distributions of the background sources from near-IR to millimetre wavelengths and use the magnification factors estimated by Dye et al. to derive the intrinsic properties of the lensed galaxies. We find these galaxies to have star-formation rates of approximately 400 to 2000 M_sol/yr, with approximately (6-25)x10^10 M_sol of their baryonic mass already turned into stars. At these rates of star formation, all remaining molecular gas will be exhausted in less than 100 Myr, reaching a final mass in stars of a few 10^11 M_sol. These galaxies are thus proto-ellipticals caught during their major episode of star formation, and observed at the peak epoch z=1.5-3 of the cosmic star formation history of the Universe.

This paper presents a compilation of clustering results taken from the literature for galaxies with highly enhanced (SFR [30-10^3] Msun/yr) star formation activity observed in the redshift range z=[0-3]. We show that, irrespective of the selection technique and only very mildly depending on the star forming rate, the clustering lengths of these objects present a sharp increase of about a factor 3 between z~1 and z~2, going from values of ~5 Mpc to about 15 Mpc and higher. This behaviour is reflected in the trend of the masses of the dark matter hosts of star-forming galaxies which increase from ~10^11.5 Msun to ~10^13.5 Msun between z~1 and z~2. Our analysis shows that galaxies which actively form stars at high redshifts are not the same population of sources we observe in the more local universe. In fact, vigorous star formation in the early universe is hosted by very massive structures, while for z~1 a comparable activity is encountered in much smaller systems, consistent with the down-sizing scenario. The available clustering data can hardly be reconciled with merging as the main trigger for intense star formation activity at high redshifts. We further argue that, after a characteristic time-scale of ~1 Gyr, massive star-forming galaxies at z>~2 evolve into z<~1.5 passive galaxies with large (Mstellar=[10^11 - 10^12] Msun) stellar masses.

Strong gravitational lenses are now being routinely discovered in wide-field surveys at (sub)millimeter wavelengths. We present Submillimeter Array (SMA) high-spatial resolution imaging and Gemini-South and Multiple Mirror Telescope optical spectroscopy of strong lens candidates discovered in the two widest extragalactic surveys conducted by the Herschel Space Observatory: the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) and the Herschel Multi-tiered Extragalactic Survey (HerMES). From a sample of 30 Herschel sources with S_500>100 mJy, 21 are strongly lensed (multiply imaged), 4 are moderately lensed (singly imaged), and the remainder require additional data to determine their lensing status. We apply a visibility-plane lens modeling technique to the SMA data to recover information about the masses of the lenses as well as the intrinsic (i.e., unlensed) sizes (r_half) and far-infrared luminosities (L_FIR) of the lensed submillimeter galaxies (SMGs). The sample of lenses comprises primarily isolated massive galaxies, but includes some groups and clusters as well. Several of the lenses are located at z_lens>0.7, a redshift regime that is inaccessible to lens searches based on Sloan Digital Sky Survey spectroscopy. The lensed SMGs are amplified by factors that are significantly below statistical model predictions given the 500um flux densities of our sample. We speculate that this may reflect a deficiency in our understanding of the intrinsic sizes and luminosities of the brightest SMGs. The lensed SMGs span nearly one decade in L_FIR (median L_FIR=7.9x10^12 L_sun) and two decades in FIR luminosity surface density (median Sigma_FIR=6.0x10^11 L_sun kpc^-2). The strong lenses in this sample and others identified via (sub-)mm surveys will provide a wealth of information regarding the astrophysics of galaxy formation and evolution over a wide range in redshift.

We have selected a complete sample of flat-spectrum radio quasars (FSRQs) from the WMAP 7-yr catalog within the SDSS area, all with measured redshift, and have compared the black hole mass estimates based on fitting a standard accretion disk model to the `blue bump' with those obtained from the commonly used single epoch virial method. The sample comprises 79 objects with a flux density limit of 1 Jy at 23 GHz, 54 of which (68%) have a clearly detected `blue bump'. Thirty-four of the latter have, in the literature, black hole mass estimates obtained with the virial method. The mass estimates obtained from the two methods are well correlated. If the calibration factor of the virial relation is set to $f=4.5$, well within the range of recent estimates, the mean logarithmic ratio of the two mass estimates is equal to zero with a dispersion close to the estimated uncertainty of the virial method. The fact that the two independent methods agree so closely in spite of the potentially large uncertainties associated with each lends strong support to both of them. The distribution of black-hole masses for the 54 FSRQs in our sample with a well detected blue bump has a median value of $7.4\times 10^{8}\,M_\odot$. It declines at the low mass end, consistent with other indications that radio loud AGNs are generally associated with the most massive black holes, although the decline may be, at least partly, due to the source selection. The distribution drops above $\log(M_\bullet/M_\odot) = 9.4$, implying that ultra-massive black holes associated with FSRQs must be rare.

We build a sample of 298 spectroscopically-confirmed galaxies at redshift z~2, selected in the z-band from the GOODS-MUSIC catalog. By exploiting the rest frame 8 um luminosity as a proxy of the star formation rate (SFR) we check the accuracy of the standard SED-fitting technique, finding it is not accurate enough to provide reliable estimates of the galaxy physical parameters. We then develop a new SED-fitting method that includes the IR luminosity as a prior and a generalized Calzetti law with a variable RV . Then we exploit such a new method to re-analyze our galaxy sample, and to robustly determine SFRs, stellar masses and ages. We find that there is a general trend of increasing attenuation with the SFR. Moreover, we find that the SFRs range between a few to 1000 solar mass per year, the masses from one billion to 400 billion solar masses, while the ages from a few tens of Myr to more than 1 Gyr. We discuss how individual age easurements of highly attenuated objects indicate that dust must form within a few tens of Myr and be copious already at ~100 Myr. In addition, we find that low luminous galaxies harbor, on average, significantly older stellar populations and are also less massive than brighter ones; we discuss how these findings and the well known 'downsizing' scenario are consistent in a framework where less massive galaxies form first, but their star formation lasts longer. Finally, we find that the near-IR attenuation is not scarce for luminous objects, contrary to what is customarily assumed; we discuss how this affects the interpretation of the observed mass-to-light ratios.

[Abridged] Since 1971 observations in X rays of thousands galaxy clusters have uncovered huge amounts of hot baryons filling up the deep gravitational potential wells provided by dark matter (DM) halos with sizes of millions light-years and masses of some 10^15 M_sun. At temperatures T~10^8 K and with average densities of n~1 particle per liter, such baryons add up to some 10^14 M_sun. With the neutralizing electrons, they constitute the best proton-electron plasma in the Universe (Intra Cluster Plasma, ICP). A key physical feature of the ICP is constituted by its good local Thermal equilibrium, and by its overall hydrostatic condition in the DM wells, modulated by entropy. The latter is set up in the cluster center by the initial halo collapse, and is progressively added at the outgrowing cluster boundary by standing shocks in the supersonic flow of intergalactic gas into the DM wells. We review these entropy-based models and discuss their outcomes and predictions concerning the ICP observables in X rays and in microwaves. The results provide a baseline for disentangling a number of additional and intriguing physical processes superposed to the general equilibrium. We cover: the central entropy erosion produced by radiative cooling vs. the intermittent energy inputs mainly due to active galactic nuclei and mergers; outer turbulent support linked with weakening shocks and decreasing inflow through the virial boundary, causing reduced entropy production; the development from high to low entropy levels throughout a typical cluster; perturbations of the equilibrium up to outright disruption due to deep impacts of infalling galaxy groups or collisions with comparable companion clusters; relativistic energy distributions of electrons accelerated during such events, producing extended radio emission by synchrotron radiation, and contributing to non-thermal pressure support for the ICP.

Upon its completion the Herschel ATLAS (H-ATLAS) will be the largest submillimetre survey to date, detecting close to half-a-million sources. It will only be possible to measure spectroscopic redshifts for a small fraction of these sources. However, if the rest-frame spectral energy distribution (SED) of a typical H-ATLAS source is known, this SED and the observed Herschel fluxes can be used to estimate the redshifts of the H-ATLAS sources without spectroscopic redshifts. In this paper, we use a subset of 40 H-ATLAS sources with previously measured redshifts in the range 0.5<z<4.2 to derive a suitable average template for high redshift H-ATLAS sources. We find that a template with two dust components T_c = 23.9 K, T_h = 46.9 K and ratio of mass of cold dust to mass of warm dust of 30.1) provides a good fit to the rest-frame fluxes of the sources in our calibration sample. We use a jackknife technique to estimate the accuracy of the redshifts estimated with this template, finding a root mean square of Delta z/(1+z) = 0.26. For sources for which there is prior information that they lie at z > 1 we estimate that the rms of Delta z/(1+z) = 0.12. We have used this template to estimate the redshift distribution for the sources detected in the H-ATLAS equatorial fields, finding a bimodal distribution with a mean redshift of 1.2, 1.9 and 2.5 for 250, 350 and 500 um selected sources respectively. \end{abstract}

We exploit the excursion set approach in integral formulation to derive novel, accurate analytic approximations of the unconditional and conditional first crossing distributions, for random walks with uncorrelated steps and general shapes of the moving barrier; we find the corresponding approximations of the unconditional and conditional halo mass functions for Cold Dark Matter power spectra to represent very well the outcomes of state-of-the-art cosmological N-body simulations. In addition, we apply these results to derive, and confront with simulations, other quantities of interest in halo statistics, including the rates of halo formation and creation, the average halo growth history, and the halo bias. Finally, we discuss how our approach and main results change when considering random walks with correlated instead of uncorrelated steps, and Warm instead of Cold Dark Matter power spectra.

The total mass derived from X-ray emission is biased low in a large number of clusters when compared with the mass estimated via strong and weak lensing. Suzaku and Chandra observations out to the virial radius report in several relaxed clusters steep temperature gradients that on assuming pure thermal hydrostatic equilibrium imply an unphysically decreasing mass profile. Moreover, the gas mass fraction appears to be inconsistent with the cosmic value measured from the CMB. Such findings can be interpreted as an evidence for an additional nonthermal pressure in the outskirts of these clusters. This nonthermal component may be due to turbulence stirred by residual bulk motions of extragalactic gas infalling into the cluster. Here we present a SuperModel analysis of Abell 1835 observed by Chandra out to the virial radius. The SuperModel formalism can include in the equilibrium a nonthermal component whose level and distribution are derived imposing that the gas mass fraction f_{gas} equals the cosmic value at the virial radius. Including such a nonthermal component, we reconstruct from X rays an increasing mass profile consistent with the hydrostatic equilibrium also in the cluster outskirts and in agreement at the virial boundary with the weak lensing value. The increasing f_{gas} profile confirms that the baryons are not missing but located at the cluster outskirts.

We combine Planck HFI data at 857, 545, 353 & 217GHz with data from WISE, Spitzer, IRAS & Herschel to investigate the properties of a flux limited sample of local star-forming galaxies. A 545GHz flux density limit was chosen so that the sample is 80% complete at this frequency, giving a sample of 234 local galaxies. We investigate the dust emission and star formation properties of the sample via various models & calculate the local dust mass function. Although 1-component modified black bodies fit the dust emission longward of 80um very well (median beta=1.83) the degeneracy between dust temp & beta also means that the SEDs are very well described by a dust emissivity index fixed at beta=2 and 10<T<25 K. Although a second, warmer dust component is required to fit shorter wavelength data, & contributes ~1/3 of the total infrared emission, its mass is negligible. No evidence is found for a very cold (6-10 K) dust component. The temp of the cold dust component is strongly influenced by the ratio of the star formation rate to the total dust mass. This implies, contrary to what is often assumed, that a significant fraction of even the emission from ~20 K dust is powered by ongoing star formation, whether or not the dust itself is associated with star forming clouds or `cirrus'. There is statistical evidence of a free-free contribution to the 217GHz flux densities of <20%. We find a median dust-to-stellar mass ratio of 0.0046; & that this ratio is anti-correlated with galaxy mass. There is good correlation between dust mass & atomic gas mass (median M_d/M_HI = 0.022), suggesting that galaxies that have more dust have more interstellar medium in general. Our derived dust mass function implies a mean dust mass density of the local Universe (for dust within galaxies), of 7.0+-1.4 x 10^5 M_solar/Mpc, significantly greater than that found in the most recent estimate using Herschel data.

In this paper we present forecasts of the contamination on different shapes of the primordial non-Gaussianity fnl parameter -- detectable on future Cosmic Microwave Background (CMB) high--resolution anisotropy maps -- produced by unresolved extragalactic point sources at frequencies of cosmological interest (45--375 GHz). We consider two scenarios: an ideal (noiseless) mission and a possible future space-borne satellite, with instrumental characteristics similar to the ones proposed for the Cosmic Origins Explorer (COrE). The local, equilateral, orthogonal and flat shapes are considered in both temperature (intensity) and polarized emission data. The angular power spectrum and bispectrum of extragalactic point sources are estimated by state-of-the-art models of source number counts. The impact of all the most relevant (far--IR and radio selected) source populations on these shapes at COrE frequencies is studied. The results of this analysis show that unresolved extragalactic point sources should not induce a very relevant non-Gaussian signal in the frequency range 100--200 GHz, thus not preventing a correct estimate of the CMB primordial fnl parameter. Polarization information allows one to significantly reduce the error--bars in the fnl parameter and the bias induced by unresolved sources and, hence, to widen the range of frequencies for fnl studies. On the contrary, at nu < 100 GHz or nu > 225 GHz, important non-Gaussian deviations in CMB anisotropy maps are expected due to unresolved extragalactic sources.

[Abridged] We present a comprehensive investigation of the cosmological evolution of the luminosity function (LF) of galaxies and active galactic nuclei (AGN) in the infrared (IR). Based on the observed dichotomy in the ages of stellar populations of early-type galaxies on one side and late-type galaxies on the other, the model interprets the epoch-dependent LFs at z \geq 1.5 using a physical model for the evolution of proto-spheroidal galaxies and of the associated AGNs, while IR galaxies at z<1.5 are interpreted as being mostly late-type 'cold' (normal) and 'warm' (starburst) galaxies. As for proto-spheroids, in addition to the epoch-dependent LFs of stellar and AGN components separately, we have worked out the evolving LFs of these objects as a whole (stellar plus AGN component). The model provides a physical explanation for the observed positive evolution of both galaxies and AGNs up to z \simeq 2.5 and for the negative evolution at higher redshifts, for the sharp transition from Euclidean to extremely steep counts at (sub-)mm wavelengths, as well as the (sub-)mm counts of strongly lensed galaxies, that are hard to account for by alternative, physical or phenomenological, approaches. The evolution of late-type galaxies and of z<1.5 AGNs is described using a parametric phenomenological approach. The modeled AGN contributions to the counts and to the cosmic infrared background (CIB) are always subdominant with maximal at mid-IR wavelengths. The model provides a good fit to the multi-wavelength (from the mid-IR to millimeter waves) data on LFs at different redshifts and on number counts (both global and per redshift slices). A prediction of the present model is a systematic variation with wavelength of the populations dominating the counts and the contributions to the CIB intensity. The implied specific trend for cross-wavelength CIB power spectra is found to be in good agreement with the data.

The Planck collaboration has recently published precise and resolved measurements of the Sunyaev-Zel'dovich effect in Abell 1656 (the Coma cluster of galaxies), so directly gauging the electron pressure profile in the intracluster plasma. On the other hand, such a quantity may be also derived from combining the density and temperature provided by X-ray observations of the thermal bremsstrahlung radiation emitted by the plasma. We find a model-independent tension between the SZ and the X-ray pressure, with the SZ one being definitely lower by 15-20%. We propose that such a challenging tension can be resolved in terms of an additional, non-thermal support to the gravitational equilibrium of the intracluster plasma. This can be straightforwardly included in our Supermodel, so as to fit in detail the Planck SZ profile while being consistent with the X-ray observables. Possible origins of the nonthermal component include cosmic-ray protons, ongoing turbulence, and relativistic electrons; given the existing observational constraints on the first two options, here we focus on the third. For this to be effective, we find that the electron population must include not only an energetic tail accelerated to gamma> 10^3 responsible for the Coma radiohalo, but also many more, lower energy electrons. The electron acceleration is to be started by merging events similar to those which provided the very high central entropy of the thermal intracluster plasma in Coma.

The Planck Early Release Compact Source Catalog (ERCSC) has offered the first opportunity to accurately determine the luminosity function of dusty galaxies in the very local Universe (i.e. distances <~ 100 Mpc), at several (sub-)millimetre wavelengths, using blindly selected samples of low redshift sources, unaffected by cosmological evolution. This project, however, requires careful consideration of a variety of issues including the choice of the appropriate flux density measurement, the separation of dusty galaxies from radio sources and from Galactic sources, the correction for the CO emission, the effect of density inhomogeneities, and more. We present estimates of the local luminosity functions at 857 GHz (350 microns), 545 GHz (550 microns) and 353 GHz (850 microns) extending across the characteristic luminosity L_star, and a preliminary estimate over a limited luminosity range at 217 GHz (1382 microns). At 850 microns and for luminosities L >~ L_star our results agree with previous estimates, derived from the SCUBA Local Universe Galaxy Survey (SLUGS), but are higher than the latter at L <~ L_star. We also find good agreement with estimates at 350 and 500 microns based on preliminary Herschel survey data.

We present a solid relationship between the neutral hydrogen (HI) disk mass and the stellar disk mass of late-type galaxies in the local universe. This relationship is derived by comparing the stellar disk mass function from the Sloan Digital Sky Survey and the HI mass function from the HI Parkes All Sky Survey (HIPASS). We find that the HI mass in late-type galaxies tightly correlates with the stellar mass over three orders of magnitude in stellar disk mass. We cross-check our result with that obtained from a sample of HIPASS objects for which the stellar mass has been obtained by inner kinematics. In addition, we derive the HI versus halo mass relationship and the dependence of all the baryonic components in spirals on the host halo mass. These relationships bear the imprint of the processes ruling galaxy formation, and highlight the inefficiency of galaxies both in forming stars and in retaining their pristine HI gas.

We have worked out simple analytical formulae that accurately approximate the relationship between the position of the source with respect to the lens center and the amplification of the images, hence the lens cross section, for realistic lens profiles. We find that, for essentially the full range of parameters either observationally determined or yielded by numerical simulations, the combination of dark matter and star distribution can be very well described, for lens radii relevant to strong lensing, by a simple power-law whose slope is very weakly dependent on the parameters characterizing the global matter surface density profile and close to isothermal in agreement with direct estimates for individual lens galaxies. Our simple treatment allows an easy insight into the role of the different ingredients that determine the lens cross section and the distribution of gravitational amplifications. They also ease the reconstruction of the lens mass distribution from the observed images and, vice-versa, allow a fast application of ray-tracing techniques to model the effect of lensing on a variety of source structures. The maximum amplification depends primarily on the source size. Amplifications larger than ~20 are indicative of compact source sizes at high-z, in agreement with expectations if galaxies formed most of their stars during the dissipative collapse of cold gas. Our formalism has allowed us to reproduce the counts of strongly lensed galaxies found in the H-ATLAS SDP field. While our analysis is focussed on spherical lenses, we also discuss the effect of ellipticity and the case of late-type lenses (showing why they are much less common, even though late-type galaxies are more numerous). Furthermore we discuss the effect of a cluster halo surrounding the early-type lens and of a supermassive black hole at its center.

We present IRAM Plateau de Bure Interferometer observations of the 12CO(3-2) emission from two far-infrared luminous QSOs at z ~ 2.5 selected from the Herschel-ATLAS survey. These far-infrared bright QSOs were selected to have supermassive black holes (SMBHs) with masses similar to those thought to reside in sub-millimetre galaxies (SMGs) at z ~ 2.5; making them ideal candidates as systems in transition from an ultraluminous infrared galaxy phase to a sub-mm faint, unobscured, QSO. We detect 12CO(3-2) emission from both QSOs and we compare their baryonic, dynamical and SMBH masses to those of SMGs at the same epoch. We find that these far-infrared bright QSOs have similar dynamical but lower gas masses than SMGs. In particular we find that far-infrared bright QSOs have ~50+-23% less warm/dense gas than SMGs, which combined with previous results showing the QSOs lack the extended, cool reservoir of gas seen in SMGs, suggests that they are at a different evolutionary stage. This is consistent with the hypothesis that far-infrared bright QSOs represent a short (~1Myr) but ubiquitous phase in the transformation of dust obscured, gas-rich, starburst-dominated SMGs into unobscured, gas-poor, QSOs.

We identify near-infrared Ks band counterparts to Herschel-ATLAS sub-mm sources, using a preliminary object catalogue from the VISTA VIKING survey. The sub-mm sources are selected from the H-ATLAS Phase 1 catalogue of the GAMA 9h field, which includes all objects detected at 250, 350 or 500 um with the SPIRE instrument. We apply and discuss a likelihood ratio (LR) method for VIKING candidates within a search radius of 10" of the 22,000 SPIRE sources with a 5 sigma detection at 250 um. We find that 11,294(51%) of the SPIRE sources have a best VIKING counterpart with a reliability $R\ge 0.8$, and the false identification rate of these is estimated to be 4.2%. We expect to miss ~5% of true VIKING counterparts. There is evidence from Z-J and J-Ks colours that the reliable counterparts to SPIRE galaxies are marginally redder than the field population. We obtain photometric redshifts for ~68% of all (non-stellar) VIKING candidates with a median redshift of 0.405. Comparing to the results of the optical identifications supplied with the Phase I catalogue, we find that the use of medium-deep near-infrared data improves the identification rate of reliable counterparts from 36% to 51%.

We have investigated the Cosmic Infrared Background (CIB) anisotropies in the framework of the physical evolutionary model for proto-spheroidal galaxies by Granato et al. (2004). After having re-calibrated the cumulative flux function $dS/dz$ at $\lambda \ge 850\,\mu$m using the available determinations of the shot noise amplitude (the original model already correctly reproduces it at shorter wavelengths) the CIB power spectra at wavelengths from $250\,\mu$m to $2\,$mm measured by {\it Planck}, {\it Herschel}, SPT and ACT experiments have been fitted using the halo model with only 2 free parameters, the minimum halo mass and the power-law index of the mean occupation function of satellite galaxies. The best-fit {\it minimum} halo mass is $\log(M_{\rm min}/M_\odot) = 12.24 \pm 0.06$, higher than, but consistent within the errors, with the estimate by Amblard et al. (2011) and close to the estimate by Planck Collaboration (2011). The redshift evolution of the volume emissivity of galaxies yielded by the model is found to be consistent with that inferred from the data. The derived {\it effective} halo mass, $M_{\rm eff} \simeq 5\times 10^{12}\,M_\odot$, of $z\simeq 2$ sub-millimeter galaxies is close to that estimated for the most efficient star-formers at the same redshift. The effective bias factor and the comoving clustering radius at $z\simeq 2$ yielded by the model are substantially lower than those found for a model whereby the star formation is fueled by steady gas accretion, but substantially higher than those found for a merging-driven galaxy evolution with a top-heavy initial mass function.

We present the results of a cross-correlation of the Planck Early Release Compact Source Catalog (ERCSC) with the catalog of Herschel-ATLAS sources detected in the Phase 1 fields, covering 134.55 deg2. There are 28 ERCSC sources detected by Planck at 857 GHz in this area. As many as 16 of them are probably high Galactic latitude cirrus; 10 additional sources can be clearly identified as bright, low-z galaxies; one further source is resolved by Herschel as two relatively bright sources; and the last is resolved into an unusual condensation of low-flux, probably high-redshift point sources, around a strongly lensed Herschel-ATLAS source at z = 3.26. Our results demonstrate that the higher sensitivity and higher angular resolution H-ATLAS maps provide essential information for the interpretation of candidate sources extracted from Planck sub-mm maps.

We use the Herschel-ATLAS survey to conduct the first large-scale statistical study of the submm properties of optically selected galaxies. Using ~80,000 r-band selected galaxies from 126 deg^2 of the GAMA survey, we stack into submm imaging at 250, 350 and 500{\mu}m to gain unprecedented statistics on the dust emission from galaxies at z < 0.35. We find that low redshift galaxies account for 5% of the cosmic 250{\mu}m background (4% at 350{\mu}m; 3% at 500{\mu}m), of which approximately 60% comes from 'blue' and 20% from 'red' galaxies (rest-frame g - r). We compare the dust properties of different galaxy populations by dividing the sample into bins of optical luminosity, stellar mass, colour and redshift. In blue galaxies we find that dust temperature and luminosity correlate strongly with stellar mass at a fixed redshift, but red galaxies do not follow these correlations and overall have lower luminosities and temperatures. We make reasonable assumptions to account for the contaminating flux from lensing by red sequence galaxies and conclude that galaxies with different optical colours have fundamentally different dust emission properties. Results indicate that while blue galaxies are more luminous than red galaxies due to higher temperatures, the dust masses of the two samples are relatively similar. Dust mass is shown to correlate with stellar mass, although the dust/stellar mass ratio is much higher for low stellar mass galaxies, consistent with the lowest mass galaxies having the highest specific star formation rates. We stack the 250{\mu}m/NUV luminosity ratio, finding results consistent with greater obscuration of star formation at lower stellar mass and higher redshift. Submm luminosities and dust masses of all galaxies are shown to evolve strongly with redshift, indicating a fall in the amount of obscured star formation in ordinary galaxies over the last four billion years.

While the selection of strongly lensed galaxies with 500{\mu}m flux density S_500>100 mJy has proven to be rather straightforward (Negrello et al. 2010), for many applications it is important to analyze samples larger than the ones obtained when confining ourselves to such a bright limit. Moreover, only by probing to fainter flux densities is possible to exploit strong lensing to investigate the bulk of the high-z star-forming galaxy population. We describe HALOS (the Herschel-ATLAS Lensed Objects Selection), a method for efficiently selecting fainter candidate strongly lensed galaxies, reaching a surface density of ~1.5-2 deg^-2, i.e. a factor of about 4 to 6 higher than that at the 100 mJy flux limit. HALOS will allow the selection of up to ~1000 candidate strongly lensed galaxies (with amplifications \mu>2) over the full H-ATLAS survey area. Applying HALOS to the H-ATLAS Science Demonstration Phase field (~14.4 deg^2) we find 31 candidate strongly lensed galaxies, whose candidate lenses are identified in the VIKING near-infrared catalog. Using the available information on candidate sources and candidate lenses we tentatively estimate a ~72% purity of the sample. The redshift distribution of the candidate lensed sources is close to that reported for most previous surveys for lensed galaxies, while that of candidate lenses extends to substantially higher redshifts than found in the other surveys. The counts of candidate strongly lensed galaxies are also in good agreement with model predictions (Lapi et al. 2011). Even though a key ingredient of the method is the deep near-infrared VIKING photometry, we show that H-ATLAS data alone allow the selection of a similarly deep sample of candidate strongly lensed galaxies with an efficiency close to 50%; a slightly lower surface density (~1.45 deg^-2) can be reached with a ~70% efficiency.

The pressure profiles of the Intracluster Plasma in galaxy clusters show a wide variance when observed in X rays at low redshifts z<0.2. We find the profiles to follow two main patterns, featuring either a steep or a shallow shape throughout both core and outskirts. We trace these shapes back to a physical dichotomy of clusters into two classes, marked by either low entropy (LE) or high entropy (HE) throughout. From X-ray observations and Sunyaev-Zel'dovich stacked data at higher 0.2<z<0.4, we elicit evidence of an increasing abundance of HEs relative to LEs. We propose this to constitute a systematic trend toward high z; specifically, we predict the pressure profiles to converge into a truly universal HE-like template for z>0.5. We submit our physical templates and converging trend for further observational tests, in view of the current and upcoming measurements of individual, stacked, and integrated Sunyaev-Zel'dovich signals.

We investigate the structure of cold dark matter halos using advanced models of spherical collapse and accretion in an expanding Universe. These base on solving time-dependent equations for the moments of the phase-space distribution function in the fluid approximation; our approach includes non-radial random motions, and most importantly, an advanced treatment of both dynamical relaxation effects that takes place in the infalling matter: phase-mixing associated to shell crossing, and collective collisions related to physical clumpiness. We find self-similar solutions for the spherically-averaged profiles of mass density rho(r), pseudo phase-space density Q(r) and anisotropy parameter beta(r). These profiles agree with the outcomes of state-of-the-art N-body simulations in the radial range currently probed by the latter; at smaller radii, we provide specific predictions. In the perspective provided by our self-similar solutions we link the halo structure to its two-stage growth history, and propose the following picture. During the early fast collapse of the inner region dominated by a few merging clumps, efficient dynamical relaxation plays a key role in producing a closely universal mass density and pseudo phase-space density profiles; in particular, these are found to depend only weakly on the detailed shape of the initial perturbation and the related collapse times. The subsequent inside-out growth of the outer regions feeds on the slow accretion of many small clumps and diffuse matter; thus the outskirts are only mildly affected by dynamical relaxation but are more sensitive to asymmetries and cosmological variance.

We take up from a library of 12 galaxy clusters featuring extended X-ray observations of their Intra Cluster Plasma (ICP), analyzed with our entropy-based Supermodel. Its few intrinsic parameters - basically, the central level and the outer slope of the entropy profile - enable us to uniformly derive not only robust snapshots of the ICP thermal state, but also the 'concentration' parameter marking the age of the host dark matter halo. We test these profiles for consistency with numerical simulations and observations. We find the central and the outer entropy to correlate, so that these clusters split into two main classes defined on the basis of low (LE) or high (HE) entropy conditions prevailing throughout the ICP. We also find inverse correlations between the central/outer entropy and the halo concentration. We interpret these in terms of mapping the ICP progress on timescales around 5 Gyr toward higher concentrations, under the drive of the dark matter halo development. The progress proceeds from HEs to LEs, toward states of deeper entropy erosion by radiative cooling in the inner regions, and of decreasing outer entropy production as the accretion peters out. We propose these radial and time features to constitute a cluster Grand Design, that we use here to derive a number of predictions. For HE clusters we predict sustained outer temperature profiles. For LEs we expect the outer entropy ramp to bend over, hence the temperature decline to steepen at low z; this feature goes together with an increasing turbulent support, a condition that can be directly probed with the SZ effect.

Exploiting the Herschel-ATLAS Science Demonstration Phase (SDP) survey data, we have determined the luminosity functions (LFs) at rest-frame wavelengths of 100 and 250 micron and at several redshifts z>1, for bright sub-mm galaxies with star formation rates (SFR) >100 M_sun/yr. We find that the evolution of the comoving LF is strong up to z~2.5, and slows down at higher redshifts. From the LFs and the information on halo masses inferred from clustering analysis, we derived an average relation between SFR and halo mass (and its scatter). We also infer that the timescale of the main episode of dust-enshrouded star formation in massive halos (M_H>3*10^12 M_sun) amounts to ~7*10^8 yr. Given the SFRs, which are in the range 10^2-10^3 M_sun/yr, this timescale implies final stellar masses of order of 10^11-10^12 M_sun. The corresponding stellar mass function matches the observed mass function of passively evolving galaxies at z>1. The comparison of the statistics for sub-mm and UV selected galaxies suggests that the dust-free, UV bright phase, is >10^2 times shorter than the sub-mm bright phase, implying that the dust must form soon after the onset of star formation. Using a single reference Spectral Energy Distribution (SED; the one of the z~2.3 galaxy SMM J2135-0102), our simple physical model is able to reproduce not only the LFs at different redshifts > 1 but also the counts at wavelengths ranging from 250 micron to ~1 mm. Owing to the steepness of the counts and their relatively broad frequency range, this result suggests that the dispersion of sub-mm SEDs of z>1 galaxies around the reference one is rather small.

The Supermodel provides an accurate description of the thermal contribution by the hot intracluster plasma which is crucial for the analysis of the hard excess. In this paper the thermal emissivity in the Coma cluster is derived starting from the intracluster gas temperature and density profiles obtained by the Supermodel analysis of X-ray observables: the XMM-Newton temperature profile and the Rosat brightness distribution. The Supermodel analysis of the BeppoSAX/PDS hard X-ray spectrum confirms our previous results, namely an excess at the c.l. of ~4.8sigma and a nonthermal flux of 1.30+-0.40x 10^-11 erg cm^-2 s^-1 in the energy range 20-80 keV. A recent joint XMM-Newton/Suzaku analysis reports an upper limit of ~6x10^-12 erg cm^-2 s^-1 in the energy range 20-80 keV for the nonthermal flux with an average gas temperature of 8.45+-0.06 keV, and an excess of nonthermal radiation at a confidence level above 4sigma, without including systematic effects, for an average XMM-Newton temperature of 8.2 keV in the Suzaku/HXD-PIN FOV, in agreement with our earlier PDS analysis. Here we present a further evidence of the compatibility between the Suzaku and BeppoSAX spectra, obtained by our Supermodel analysis of the PDS data, when the smaller size of the HXD-PIN FOV and the two different average temperatures derived by XMM-Newton and by the joint XMM-Newton/Suzaku analysis are taken into account. The consistency of the PDS and HXD-PIN spectra reaffirms the presence of a nonthermal component in the hard X-ray spectrum of the Coma cluster. The Supermodel analysis of the PDS data reports an excess at c.l. above 4sigma also for the higher average temperature of 8.45 keV thanks to the PDS FOV considerably greater than the HXD-PIN FOV.

The Seminar "Dark Matter in Galaxies" was delivered, within the Dark Matter Awareness Week (1-8 December 2010) at 140 institutes in 46 countries and it was followed by 4200 people. A documentation of this worldwide initiative is at http://www.youtube.com/watch?v=AOBit8a-1Fw. A reference presentation, prepared by a coordinated pool of leading scientists in the field, was available to speakers. In response to feedbacks and suggestions, we upgraded it to a "Presentation Review" of which we provide here the Explanatory Notes, the link to the .pptx file, and some image of the slides. This Presentation Review is an innovative scientific product to meet the request of information about the phenomenology of the DM mystery at galactic scale. This Review concerns the mass discrepancy phenomenon detected in galaxies, usually accounted by postulating the presence of a non luminous non baryonic component. In the theoretical framework of Newtonian Gravity we recall the properties of Dark Matter halos as emerging from the state-of-the-art of numerical simulations performed in the current $\Lambda CDM$ scenario. Then, the simple but much-telling phenomenology of the distribution of dark and luminous matter in Spirals, Ellipticals, and dwarf Spheroidals is reported. We show that a coherent observational framework emerges from reliable data of different large samples of objects and it is obtained by different methods of investigation. We then highlight the impressive evidence that the distribution of dark and luminous matter are closely correlated and that have universal features. Hints on the cosmological role of this phenomenological scenario are then given. Finally, we discuss the constraints on the elusive nature of the dark particle that the actual distribution of DM around galaxies pose on its direct and indirect searches.

We submit that non thermalized support for the outer intracluster medium in relaxed galaxy clusters is provided by turbulence, driven by inflows of intergalactic gas across the virial accretion shocks. We expect this component to increase briskly during the cluster development for z<1/2, due to three factors. First, the accretion rates of gas and dark matter subside, when they feed on the outer wings of the initial perturbations in the accelerating Universe. Second, the infall speeds decrease across the progressively shallower gravitational potential at the shock position. Third, the shocks eventually weaken, and leave less thermal energy to feed the intracluster entropy, but relatively more bulk energy to drive turbulence into the outskirts. The overall outcome from these factors is physically modeled and analytically computed; thus we ascertain how these concur in setting the equilibrium of the outer intracluster medium, and predict how the observables in X rays and microwaves are affected, so as to probe the development of outer turbulence over wide cluster samples. By the same token, we quantify the resulting negative bias to be expected in the total mass evaluated from X-ray measurements.

We investigate the dynamical basis of the classic empirical models (specifically, Sersic-Einasto and generalized NFW) that are widely used to describe the distributions of collisionless matter in galaxies. We submit that such a basis is provided by our \alpha-profiles, shown to constitute solutions of the Jeans dynamical equilibrium with physical boundary conditions. We show how to set the parameters of the empirical in terms of the dynamical models; we find the empirical models, and specifically Sersic-Einasto, to constitute a simple and close approximation to the dynamical models. Finally, we discuss how these provide an useful baseline for assessing the impact of the small-scale dynamics that may modulate the density slope in the central galaxy regions.

[abridged] Massive, passively evolving galaxies at redshifts z>1 exhibit on the average physical sizes smaller by factors ~3 than local early type galaxies (ETGs) endowed with the same stellar mass. Small sizes are in fact expected on theoretical grounds, if dissipative collapse occurs. Recent results show that the size evolution at z<1 is limited to less than 40%, while most of the evolution occurs at z>1, where both compact and already extended galaxies are observed and the scatter in size is remarkably larger than locally. The presence at high z of a significant number of ETGs with the same size as their local counterparts as well as of ETGs with quite small size, points to a timescale to reach the new, expanded equilibrium configuration of less than the Hubble time. We demonstrate that the projected mass of compact, high-z galaxies and that of local ETGs within the *same physical radius*, the nominal half-luminosity radius of high-z ETGs, differ substantially, in that the high-z ETGs are on the average significantly denser. We propose that quasar activity, which peaks at z~2, can remove large amounts of gas from central galaxy regions on a timescale shorter than of the dynamical one, triggering a puffing up of the stellar component at constant stellar mass; in this case the size increase goes together with a decrease of the central mass. The size evolution is expected to parallel that of the quasars and the inverse hierarchy, or downsizing, seen in the quasar evolution is mirrored in the size evolution.

To investigate the poorly constrained sub-mm counts and spectral properties of blazars we searched for these in the Herschel-ATLAS (H-ATLAS) science demostration phase (SDP) survey catalog. We cross-matched 500$\mu$m sources brighter than 50 mJy with the FIRST radio catalogue. We found two blazars, both previously known. Our study is among the first blind blazar searches at sub-mm wavelengths, i.e., in the spectral regime where little is still known about the blazar SEDs, but where the synchrotron peak of the most luminous blazars is expected to occur. Our early results are consistent with educated extrapolations of lower frequency counts and question indications of substantial spectral curvature downwards and of spectral upturns at mm wavelengths. One of the two blazars is identified with a Fermi/LAT $\gamma$-ray source and a WMAP source. The physical parameters of the two blazars are briefly discussed.These observations demonstrate that the H-ATLAS survey will provide key information about the physics of blazars and their contribution to sub-mm counts.

In galaxy clusters the entropy distribution of the IntraCluster Plasma modulates the latter's equilibrium within the Dark Matter gravitational wells, as rendered by our Supermodel. We argue the entropy production at the boundary shocks to be reduced or terminated as the accretion rates of DM and intergalactic gas peter out; this behavior is enforced by the slowdown in the outskirt development at late times, when the Dark Energy dominates the cosmology while the outer wings of the initial perturbation drive the growth. In such conditions, we predict the ICP temperature profiles to steepen into the cluster outskirts. The detailed expectations from our simple formalism agree with the X-ray data concerning five clusters whose temperature profiles have been recently measured out to the virial radius. We predict steep temperature declines to prevail in clusters at low redshift, tempered only by rich environs including adjacent filamentary structures.

The SAGACE experiment consists of a mm/sub-mm telescope with a 3-m diameter primary mirror, coupled to a cryogenic multi-beam differential spectrometer. SAGACE explores the sky in the 100-760 GHz frequency range, using four diffraction-limited bolometer arrays. The instrument is designed to perform spectroscopic surveys of the Sunyaev-Zeldovich effects of thousands of galaxy clusters, of the spectral energy distribution of active galactic nuclei, and of the [CII] line of a thousand galaxies in the redshift desert. In 2008 a full phase-A study for a national small mission was completed and delivered to the Italian Space Agency (ASI). We have shown that taking advantage of the differential operation of the Fourier Transform Spectrometer, this ambitious instrument can operate from a Molniya orbit, and can be built and operated within the tight budget of a small mission.

We discuss the gamma-ray signal to be expected from dark matter (DM)annihilations at the Galactic Center. To describe the DM distribution in the Galactic halo we base on the Jeans equation for self-gravitating, anisotropic equilibria. In solving the Jeans equation, we adopt the specific correlation between the density \rho(r) and the velocity dispersion \sigma^2_r(r) expressed by the powerlaw behavior of the DM `entropy' K= \sigma_r^2/\rho^{2/3} ~ r^\alpha with \alpha ~ 1.25-1.3. Indicated (among others) by several recent N-body simulations, this correlation is privileged by the form of the radial pressure term in the Jeans equation, and yields a main body profile consistent with the classic self-similar development of DM halos. In addition, we require the Jeans solutions to satisfy regular boundary conditions both at the center (finite pressure, round gravitational potential) and in the outskirts (finite overall mass). With these building blocks we derive physical solutions, dubbed `\alpha-profiles'. We find the one with \alpha=1.25, suitable for the Galaxy halo, to be intrinsically flatter at the center relative to the empirical NFW formula, yet steeper than the empirical Einasto profile. So on scales of 10^{-1} deg it yields annihilation fluxes lower by a factor 5 than the former yet higher by a factor 10 than the latter; such fluxes will eventually fall within the reach of the Fermi satellite. We show the effectiveness of the \alpha-profile in relieving the astrophysical uncertainties related to the macroscopic DM distribution, and discuss its expected performance as a tool instrumental to interpret the upcoming gamma-ray data in terms of DM annihilation.

We investigate the atomic (HI) and molecular (H_2) Hydrogen content of normal galaxies by combining observational studies linking galaxy stellar and gas budgets to their host dark matter (DM) properties, with a physically grounded galaxy formation model. This enables us to analyse empirical relationships between the virial, stellar, and gaseous masses of galaxies and explore their physical origins. Utilising a semi-analytic model (SAM) to study the evolution of baryonic material within evolving DM halos, we study the effects of baryonic infall and various star formation and feedback mechanisms on the properties of formed galaxies using the most up-to-date physical recipes. We find that in order to significantly improve agreement with observations of low-mass galaxies we must suppress the infall of baryonic material and exploit a two-phase interstellar medium (ISM), where the ratio of HI to H_2 is determined by the galactic disk structure. Modifying the standard Schmidt-Kennicutt star formation law, which acts upon the total cold gas in galaxy discs and includes a critical density threshold, and employing a star formation law which correlates with the H_2 gas mass results in a lower overall star formation rate. This in turn, allows us to simultaneously reproduce stellar, HI and H_2 mass functions of normal galaxies.

We use our model for the formation and evolution of galaxies within a two-phase galaxy formation scenario, showing that the high-redshift domain typically supports the growth of spheroidal systems, whereas at low redshifts the predominant baryonic growth mechanism is quiescent and may therefore support the growth of a disc structure. Under this framework we investigate the evolving galaxy population by comparing key observations at both low and high-redshifts, finding generally good agreement. By analysing the evolutionary properties of this model, we are able to recreate several features of the evolving galaxy population with redshift, naturally reproducing number counts of massive star-forming galaxies at high redshifts, along with the galaxy scaling relations, star formation rate density and evolution of the stellar mass function. Building upon these encouraging agreements, we make model predictions that can be tested by future observations. In particular, we present the expected evolution to z=2 of the super-massive black hole mass function, and we show that the gas fraction in galaxies should decrease with increasing redshift in a mass, with more and more evolution going to higher and higher masses. Also, the characteristic transition mass from disc to bulge dominated system should decrease with increasing redshift.

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 510 square degrees of the extragalactic sky, four times larger than all the other Herschel surveys combined, in five far-infrared and submillimetre bands. We describe the survey, the complementary multi-wavelength datasets that will be combined with the Herschel data, and the six major science programmes we are undertaking. Using new models based on a previous submillimetre survey of galaxies, we present predictions of the properties of the ATLAS sources in other wavebands.

[abridged] We present the analysis of the X-ray brightness and temperature profiles for six clusters belonging to both the Cool Core and Non Cool Core classes, in terms of the Supermodel (SM) developed by Cavaliere, Lapi & Fusco-Femiano (2009). Based on the gravitational wells set by the dark matter halos, the SM straightforwardly expresses the equilibrium of the IntraCluster Plasma (ICP) modulated by the entropy deposited at the boundary by standing shocks from gravitational accretion, and injected at the center by outgoing blastwaves from mergers or from outbursts of Active Galactic Nuclei. The cluster set analyzed here highlights not only how simply the SM represents the main dichotomy Cool vs. Non Cool Core clusters in terms of a few ICP parameters governing the radial entropy run, but also how accurately it fits even complex brightness and temperature profiles. For Cool Core clusters like A2199 and A2597, the SM with a low level of central entropy straightforwardly yields the characteristic peaked profile of the temperature marked by a decline toward the center, without requiring currently strong radiative cooling and high mass deposition rates. Non Cool Core clusters like A1656 require instead a central entropy floor of a substantial level, and some like A2256 and even more A644 feature structured temperature profiles that also call for a definite floor extension; in such conditions the SM accurately fits the observations, and suggests that in these clusters the ICP has been just remolded by a merger event, in the way of a remnant cool core. The SM also predicts that dark matter halos with high concentration should correlate with flatter entropy profiles and steeper brightness in the outskirts; this is indeed the case with A1689.

We propose and test a scenario for the assembly and evolution of luminous matter in galaxies which substantially differs from that adopted by other semianalytic models. As for the dark matter (DM), we follow the detailed evolution of halos within the canonical LCDM cosmology using standard Montecarlo methods. However, when overlaying prescriptions for baryon evolution, we take into account an effect pointed out in the past few years by a number of studies mostly based on intensive N-body simulations, namely that typical halo growth occurs in two phases: an early, fast collapse phase featuring several major merger events, followed by a late, quiescent accretion onto the halo outskirts. We propose that the two modes of halo growth drive two distinct modes for the evolution of baryonic matter, favoring the development of the spheroidal and disc components of galaxies, respectively. We test this idea using the semianalytic technique. Our galaxy formation model envisages an early coevolution of spheroids and the central supermassive black holes, already tested in our previous works, followed by a relatively quiescent growth of discs around the preformed spheroids. In this exploratory study, we couple our model to the spectrophotometric code GRASIL, and compare our results on several properties of the local galaxy population with observations, Finding an encouraging agreement.

[abridged] The equilibria of the intracluster plasma (ICP) and of the gravitationally dominant dark matter (DM) are governed by the hydrostatic and the Jeans equation. Jeans, with the DM `entropy' set to K ~ r^\alpha and \alpha ~ 1.25 - 1.3 applying from groups to rich clusters, yields our radial \alpha-profiles. In the ICP the entropy run k(r) is mainly shaped by shocks, as steadily set by supersonic accretion of gas at the cluster boundary, and intermittently driven from the center by merging events or by AGNs; the resulting equilibrium is described by the exact yet simple formalism constituting our ICP Supermodel. With a few parameters, this accurately represents the runs of density n(r) and temperature T(r) as required by recent X-ray data on surface brightness and spectroscopy for both cool core (CC) and non cool core (NCC) clusters; the former are marked by a middle temperature peak, whose location is predicted from rich clusters to groups. The Supermodel inversely links the inner runs of n(r) and T(r), and highlights their central scaling with entropy n_c ~ k_c^-1 and T_c ~ k_c^0.35, to yield radiative cooling times t_c ~ 0.3 (k_c/15 keV cm^2)^1.2 Gyr. We discuss the stability of the central values so focused either in CC and NCC clusters. From the Supermodel we derive as limiting cases the classic polytropic \beta-models, and the `mirror' model with T(r) ~ \sigma^2(r) suitable for NCC and CC clusters, respectively; these highlight how the ICP temperature T(r) tends to mirror the DM velocity dispersion \sigma^2(r) away from entropy injections. Finally, we discuss how the Supermodel connects information derived from X-ray and gravitational lensing observations.

We test with gravitational lensing data the dark matter (DM) halos embedding the luminous baryonic component of galaxy clusters; our benchmark is provided by their two-stage cosmogonical development that we compute with its variance, and by the related '\alpha-profiles' we derive. The latter solve the Jeans equation for the self-gravitating, anisotropic DM equilibria, and yield the radial runs of the density \rho(r) and the velocity dispersion \sigma_r^2(r) in terms of the DM 'entropy' K = \sigma_r^2/\rho^(2/3) ~ r^(\alpha) highlighted by recent N-body simulations; the former constrains the slope to the narrow range \alpha ~ 1.25 - 1.3. These physically based \alpha-profiles meet the overall requirements from gravitational lensing observations, being intrinsically flatter at the center and steeper in the outskirts relative to the empirical NFW formula. Specifically, we project them along the l.o.s. and compare with a recent extensive dataset from strong and weak lensing observations in and around the cluster A1689. We find an optimal fit at both small and large scales in terms of a halo constituted by an early body with \alpha ~ 1.25 and by recent extensive outskirts, that make up an overall mass 10^15 M_sun with a concentration parameter c ~ 10 consistent with the variance we compute in the \LambdaCDM cosmogony. The resulting structure corresponds to a potential well shallow in the outskirts as that inferred from the X rays radiated from the hot electrons and baryons constituting the intracluster plasma.

[abridged] In the dark matter (DM) halos embedding galaxies and galaxy systems the `entropy' K = \sigma^2 / \rho^{2/3} (a quantity that combines the radial velocity dispersion \sigma with the density \rho) is found from intensive N-body simulations to follow a powerlaw run K ~ r^{\alpha} throughout the halos' bulk, with \alpha around 1.25. Taking up from phenomenology just that \alpha ~ const applies, we cut through the rich analytic contents of the Jeans equation describing the self-gravitating equilibria of the DM; we specifically focus on computing and discussing a set of novel physical solutions that we name \alpha-profiles, marked by the entropy slope \alpha itself, and by the maximal gravitational pull \kappa_crit required for a viable equilibrium to hold. We then use an advanced semianalytic description for the cosmological buildup of halos to constrain the values of \alpha to within the narrow range 1.25-1.29 from galaxies to galaxy systems. Our range of \alpha applies since the transition time that - both in oursemianalytic description and in state-of-the-art numerical simulations - separates two development stages: an early violent collapse that comprises a few major mergers and enforces dynamical mixing, followed by smoother mass addition through slow accretion. We also give an accurate analytic representation of the \alpha-profiles with parameters derived from the Jeans equation. We finally stress how our findings and predictions as to \alpha and \kappa_crit contribute to understand hitherto unsolved issues concerning the fundamental structure of DM halos.

Observations have evidenced that passively evolving massive galaxies at high redshift are much more compact than local galaxies with the same stellar mass. We argue that the observed strong evolution in size is directly related to the quasar feedback, which removes huge amounts of cold gas from the central regions in a Salpeter time, inducing an expansion of the stellar distribution. The new equilibrium configuration, with a size increased by a factor >~ 3, is attained after ~ 40 dynamical times, corresponding to ~ 2 Gyr. This means that massive galaxies observed at z >~ 1 will settle on the Fundamental Plane by z ~ 0.8-1. In less massive galaxies (M_star <~ 2 10^10 M_sun), the nuclear feedback is subdominant, and the mass loss is mainly due to stellar winds. In this case, the mass loss timescale is longer than the dynamical time and results in adiabatic expansion that may increase the effective radius by a factor of up to ~ 2 in 10 Gyr, although a growth by a factor of ~ 1.6 occurs within the first 0.5 Gyr. Since observations are focused on relatively old galaxies, with ages >~ 1 Gyr, the evolution for smaller galaxies is more difficult to perceive. Significant evolution of velocity dispersion is predicted for both small and large galaxies.

Motivated by the recent observational and theoretical evidence that long Gamma-Ray Bursts (GRBs) are likely associated with low metallicity, rapidly rotating massive stars, we examine the cosmological star formation rate (SFR) below a critical metallicity Z_crit Z_sun/10 - Z_sun/5, to estimate the event rate of high-redshift long GRB progenitors. To this purpose, we exploit a galaxy formation scenario already successfully tested on a wealth of observational data on (proto)spheroids, Lyman break galaxies, Lyman alpha emitters, submm galaxies, quasars, and local early-type galaxies. We find that the predicted rate of long GRBs amounts to about 300 events/yr/sr, of which about 30 per cent occur at z>~6. Correspondingly, the GRB number counts well agree with the bright SWIFT data, without the need for an intrinsic luminosity evolution. Moreover, the above framework enables us to predict properties of the GRB host galaxies. Most GRBs are associated with low mass galaxy halos M_H<~10^11 M_sun, and effectively trace the formation of small galaxies in such halos. The hosts are young, with age smaller than 5*10^7 yr, gas rich, but poorly extincted (A_V<~0.1) because of their chemical immaturity; this also implies high specific SFR and quite extreme alpha-enhancement. Only the minority of hosts residing in large halos with M_H>~10^12 M_sun have larger extinction (A_V~0.7-1), SFRs exceeding 100 M_sun/yr and can be detected at submm wavelengths. Most of the hosts have UV magnitudes in the range -20 <~M_1350<~ -16, and Lyman alpha luminosity in the range 2*10^40 <~L_Lya<~2*10^42 erg/s. GRB hosts are thus tracing the faint end of the luminosity function of Lyman break galaxies and Lyman alpha emitters.

From clusters to groups of galaxies, the powerful bremsstrahlung radiation L_X emitted in X rays by the intracluster plasma is observed to decline sharply with lowering virial temperatures T (i.e., at shallower depths of the gravitational wells) after a steep local L_X-T correlation; this implies increasing scarcity of diffuse baryons relative to dark matter, well under the cosmic fraction. We show how the widely debated issue concerning these `missing baryons' is solved in terms of the thermal and/or dynamical effects of the kinetic (at low redshifts z) and radiative (at high z) energy inputs from central active galactic nuclei, of which independent evidence is being observed. From these inputs we compute shape and z-evolution expected for L_X-T correlation which agree with the existing data, and provide a predictive pattern for future observations.

In the standard galaxy formation scenario plasma clouds with a high thermal energy content must exist at high redshifts since the proto-galactic gas is shock heated to the virial temperature, and extensive cooling, leading to efficient star formation, must await the collapse of massive halos (as indicated by the massive body of evidence, referred to as "downsizing"). Massive plasma clouds are potentially observable through the thermal and kinetic Sunyaev-Zel'dovich effects and their free-free emission. We find that the detection of substantial numbers of galaxy-scale thermal SZ signals is achievable by blind surveys with next generation radio telescope arrays such as EVLA, ALMA and SKA. This population is even detectable with the 10% SKA, and wide field of view options at high frequency on any of these arrays would greatly increase survey speed. An analysis of confusion effects and of the contamination by radio and dust emissions shows that the optimal frequencies are those in the range 10-35 GHz. Predictions for the redshift distributions of detected sources are also worked out.

We present a detailed study of the structural properties of four gravitationally lensed disk galaxies at z=1. Modelling the rotation curves on sub-kpc scales we derive the values for the disk mass, the reference dark matter density and core radius, and the angular momentum per unit mass. The derived models suggest that the rotation curve profile and amplitude are best fit with a dark matter component similar to those of local spiral galaxies. The stellar component also has a similar length scale, but with substantially smaller masses than similarly luminous disk galaxies in the local universe. Comparing the average dark matter density inside the optical radius we find that the disk galaxies at z=1 have larger densities (by up to a factor of 7) than similar disk galaxies in the local Universe. Furthermore, the angular momentum per unit mass versus reference velocity is well matched to the local relation, suggesting that the angular momentum of the disk remains constant between high redshifts and the present day. Though statistically limited, these observations point towards a spirals' formation scenario in which stellar disks are slowly grown by the accretion of angular momentum conserving material.

We explore the onset of star formation in the early Universe, exploiting the observations of high-redshift Lyman-break galaxies (LBGs) and Lyman alpha emitters (LAEs), in the framework of the galaxy formation scenario elaborated by Granato et al. (2004) already successfully tested against the wealth of data on later evolutionary stages. Complementing the model with a simple, physically plausible, recipe for the evolution of dust attenuation in metal poor galaxies we reproduce the luminosity functions (LFs) of LBGs and of LAEs at different redshifts. This recipe yields a much faster increase with galactic age of attenuation in more massive galaxies, endowed with higher star formation rates. These objects have therefore shorter lifetimes in the LAE and LBG phases, and are more easily detected in the dusty submillimeter bright phase. The short UV bright lifetimes of massive objects strongly mitigate the effect of the fast increase of the massive halo density with decreasing redshift, thus accounting for the weaker evolution of the LBG LF, compared to that of the halo mass function, and the even weaker evolution between z~6 and z~3 of the LAE LF. LAEs are on the average expected to be younger, with lower stellar masses, and associated to less massive halos than LBGs. Finally, we show that the intergalactic medium can be completely reionized at redshift z~6-7 by massive stars shining in protogalactic spheroids with halo masses from a few 10^10 to a few 10^11 M_sun, showing up as faint LBGs with magnitude in the range -17<M_1350<-20, without resorting to any special stellar initial mass function.

In the current LambdaCDM cosmological scenario, N-body simulations provide us with a Universal mass profile, and consequently a Universal equilibrium circular velocity of the virialized objects, as galaxies. In this paper we obtain, by combining kinematical data of their inner regions with global observational properties, the Universal Rotation Curve (URC) of disk galaxies and the corresponding mass distribution out to their virial radius. This curve extends the results of Paper I, concerning the inner luminous regions of Sb-Im spirals, out to the edge of the galaxy halos.

The ~800 optically unseen (R>25.5) 24mum-selected sources in the complete Spitzer First Look Survey sample (Fadda et al. 2006) with F[24mum]>0.35 mJy are found to be very strongly clustered. If, as indicated by several lines of circumstantial evidence, they are ultraluminous far-IR galaxies at z ~ [1.6-2.7], the amplitude of their spatial correlation function is very high. The associated comoving clustering length is estimated to be r_0=14.0_{-2.4}^{+2.1} Mpc, value which puts these sources amongst the most strongly clustered populations of our known universe. Their 8mum-24mum colours suggest that the AGN contribution dominates above F[24mum] ~ 0.8 mJy, consistent with earlier analyses. The properties of these objects (number counts, redshift distribution, clustering amplitude) are fully consistent with those of proto-spheroidal galaxies in the process of forming most of their stars and of growing their active nucleus, as described by the Granato et al. (2004) model. In particular, the inferred space density of such galaxies at z ~ 2 is much higher than what expected from most semi-analytic models. Matches of the observed projected correlation function w(\theta) with models derived within the so-called Halo Occupation Scenario show that these sources have to be hosted by haloes more massive than ~10^{13.4} M_\odot. This value is significantly higher than that for the typical galactic haloes hosting massive elliptical galaxies, suggesting a duration of the starburst phase of massive high-redshift dusty galaxies of T_B ~ 0.5 Gyr.

Shocks and blastwaves are expected to be driven driven into the intracluster medium filling galaxy groups and clusters by powerful outbursts of active galactic nuclei or quasars in the member galaxies; the first footprints of shock fronts have been tentatively traced out with X-ray imaging. We show how overpressures in the blasts behind the shock can prove the case and also provide specific marks of the nuclear activity: its strength, its current stage, and the nature of its prevailing output. We propose to detect these marks with the aimed pressure probe constituted by the resolved Sunyaev-Zel'dovich effect. We compute and discuss the outcomes to be expected in nearby and distant sources at different stages of their activity.

We show that our previously proposed anti-hierarchical baryon collapse scenario for the joint evolution of black holes and host galaxies predicts quasar luminosity functions at redshifts 1.5<z<6 and local properties in nice agreement with observations. In our model the quasar activity marks and originates the transition between an earlier phase of violent and heavily dust-enshrouded starburst activity promoting rapid black hole growth, and a later phase of almost passive evolution; the former is traced by the submillimeter-selected sources, while the latter accounts for the high number density of massive galaxies at substantial redshifts z>1.5, the population of Extremely Red Objects, and the properties of local ellipticals.

We propose that angular momentum transfer from the baryons to the Dark Matter (DM) during the early stages of galaxy formation can flatten the halo inner density profile and modify the halo dynamics. We compute the phase-space distribution function of DM halos, that corresponds to the density and anisotropy profiles obtained from N-body simulations in the concordance cosmology. We then describe an injection of angular momentum into the halo by modifying the distribution function, and show that the system evolves into a new equilibrium configuration; the latter features a constant central density and a tangentially-dominated anisotropy profile in the inner regions, while the structure is nearly unchanged beyond 10% of the virial radius. Then we propose a toy model to account for such a halo evolution, based on the angular momentum exchange due to dynamical friction; at the epoch of galaxy formation this is efficiently exerted by the DM onto the gas clouds spiralling down the potential well. The comparison between the angular momentum profile gained by the halo through dynamical friction and that provided by the perturbed distribution function reveals a surprising similarity, hinting at the reliability of the process.

We show that the ABC scenario we proposed for the co-evolution of spheroids and QSOs predicts accretion rates and masses of supermassive black holes in sub-mm galaxies in keeping with recent X-ray determinations. These masses are well below the local values, and those predicted by alternative models. The observed column densities may be mostly due to ISM in the galaxy. The contribution of the associated nuclear activity to the X-ray background is likely negligible, while they may contribute a sizeable fraction $\sim 10 %$ to hard-X cumulative counts at the faintest observed fluxes.

At variance with most semi-analytic models, in the Anti-hierarchical Baryon Collapse scenario (Granato et al. 2001, 2004) the main driver of the galaxy formation and evolution is not the merging sequence but are baryon processes. This approach emphasizes, still in the framework of the hierarchical clustering paradigm for dark matter halos, feedback processes from supernova explosions and from active nuclei, that tie together star formation in spheroidal galaxies and the growth of black holes at their centers. We review some recent results showing the remarkably successful predictive power of this scenario, which allows us to account for the evolution with cosmic time of a broad variety of properties of galaxies and active nuclei, which proved to be very challenging for competing models.

We present new relationships between halo masses M_h and several galaxy properties, including r*-band luminosities L_r, stellar (M_star) and baryonic masses, stellar velocity dispersions (sigma), and black hole masses (M_BH). Approximate analytic expressions are given. In the galaxy halo mass range 3*10^(10) M_sun <= M_h <= 3*10^(13) M_sun the M_h-L_r, M_star-M_h, and M_BH-M_h are well represented by a double power law, with a break at M_h_break~3*10^(11) M_sun, corresponding to a mass in stars M_star~ 1.2*10^(10) M_sun, to a r*-band luminosity L_r~ 5*10^9 L_sun, to a stellar velocity dispersion sigma ~ 88 km/s, and to a black hole mass M_BH~ 9*10^(6) M_sun. The sigma-M_h relation can be approximated by a single power law, though a double power law is a better representation. We interpret these relations in terms of the effect of feedback from supernovae and from the active nucleus on the interstellar medium. We argue that the break of the power laws occurs at a mass which marks the transition between the dominance of the stellar and the AGN feedback.

We compute the angular momentum, the spin parameter and the related distribution function for Dark Matter halos hosting a spiral galaxy. We base on scaling laws, inferred from observations, that link the properties of the galaxy to those of the host halo; we further assume that the Dark Matter has the same total specific angular momentum of the baryons. Our main results are: (i) we find that the gas component of the disk significantly contributes to the total angular momentum of the system; (ii) by adopting for the Dark Matter the observationally supported Burkert profile, we compute the total angular momentum of the disk and its correlation with the rotation velocity; (iii) we find that the distribution function of the spin parameter $\lambda$ peaks at a value of about 0.03, consistent with a no-major-merger scenario for the late evolution of spiral galaxies.

The temperature and density of the hot diffuse medium pervading galaxy groups and clusters combine into one significant quantity, the entropy. Here we express the entropy levels and profiles in model-independent forms by joining two observables, the X-ray luminosity and the change in the CMB intensity due to the Sunyaev-Zel'dovich (SZ) effect. Thus we present both global scaling relations for the entropy levels from clusters and groups, and a simple expression yielding the entropy profiles in individual clusters from resolved X-ray surface brightness and SZ spatial distributions. We propose that our approach provides two useful tools for comparing large data samples with models, in order to probe the processes that govern the thermal state of the hot intracluster medium. The feasibility of using such a diagnostic for the entropy is quantitatively assessed, based on current X-ray and upcoming SZ measurements.

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.

The Sunyaev-Zeldovich (SZ) effect has proven to be an extremely powerful tool to study the physical and evolutionary properties of rich clusters of galaxies. Upcoming SZ experiments, with their much improved sensitivity and angular resolution, will provide unique information also on phases of galaxy evolution characterized by the presence of large amounts of hot proto-galactic gas. We present a preliminary analysis of the SZ signals that can be expected at the collapse of the proto-galaxy when, according to the standard scenario, the gas is heated at its virial temperature, and during episodes of strong energy injections from the active nucleus. The contributions of such signals to excess power on arc-minute scales recently found by CBI and BIMA experiments are briefly discussed.

The energy fed by active galactic nuclei to the surrounding diffuse baryons changes their amount, temperature, and distribution; so in groups and in member galaxies it affects the X-ray luminosity and also the Sunyaev-Zel'dovich effect. Here we compute how the latter is enhanced by the transient blastwave driven by an active quasar, and is depressed when the equilibrium is recovered with a depleted density. We constrain such depressions and enhancements with the masses of relic black holes in galaxies and the X-ray luminosities in groups. We discuss how all these linked observables can tell the quasar contribution to the thermal history of the baryons pervading galaxies and groups.

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.