Source author record

A. Lapi

A. Lapi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

70works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

70 published item(s)

preprint2022arXiv

Astroparticle Constraints from Cosmic Reionization and Primordial Galaxy Formation

We derive astroparticle constraints in different dark matter scenarios alternative to cold dark matter (CDM): thermal relic warm dark matter, WDM; fuzzy dark matter, $ψ$DM; self-interacting dark matter, SIDM; sterile neutrino dark matter, $ν$DM. Our framework is based on updated determinations of the high-redshift UV luminosity functions for primordial galaxies out to redshift $z\sim 10$, on redshift-dependent halo mass functions in the above DM scenarios from numerical simulations, and on robust constraints on the reionization history of the Universe from recent astrophysical and cosmological datasets. First, we build up an empirical model of cosmic reionization characterized by two parameters, namely the escape fraction $f_{\rm esc}$ of ionizing photons from primordial galaxies, and the limiting UV magnitude $M_{\rm UV}^{\rm lim}$ down to which the extrapolated UV luminosity functions are steeply increasing. Second, we perform standard abundance matching of the UV luminosity function and the halo mass function, obtaining a relationship between UV luminosity and halo mass whose shape depends on an astroparticle quantity $X$ specific of each DM scenario (e.g., WDM particle mass); we exploit such a relation to introduce in the analysis a constraint from primordial galaxy formation, in terms of the threshold halo mass above which primordial galaxies can efficiently form stars. Third, we implement a sequential updating Bayesian MCMC technique to perform joint inference on the three parameters $f_{\rm esc}$, $M_{\rm UV}^{\rm lim}$, $X$, and to compare the outcomes of different DM scenarios on the reionization history. Finally, we highlight the relevance of our astroparticle estimates in predicting the behavior of the high-redshift UV luminosity function at faint, yet unexplored magnitudes, that may be tested with the advent of the James Webb Space Telescope.

preprint2022arXiv

Selecting a complete sample of blazars in sub-millimetre catalogues

The \textit{Herschel} Astrophysical Terahertz Large Area Survey (H-ATLAS), that has covered about 642 sq. deg. in 5 bands from 100 to 500 $μ\rm m$, allows a blind flux-limited selection of blazars at sub-mm wavelengths. However, blazars constitute a tiny fraction of H-ATLAS sources and therefore identifying them is not a trivial task. Using the data on known blazars detected by the H-ATLAS we have defined a locus for 500\,$μ$m selected blazars and exploited it to select blazar candidates in the H-ATLAS fields. Candidates and known blazars in the H-ATLAS equatorial and South Galactic Pole fields were followed up with the Australia Telescope Compact Array (ATCA) or with the Karl G. Jansky Very Large Array (VLA), and matched with existing radio- and mm-catalogues to reconstruct the spectral behaviour over at least 6 orders of magnitude in frequency. We identified a selection approach that, combining the information in the sub-mm and radio domains, efficiently singles out genuine blazars. In this way, we identified a sample of 39 blazars brighter than $S_{500μ\rm m} = 35\,$mJy in the H-ATLAS fields. Tests made cross-matching the H-ATLAS catalogues with large catalogues of blazar candidates indicate that the sample is complete. The derived counts are compared with model predictions finding good consistency with the C2Ex model and with estimates based on ALMA data.

preprint2022arXiv

The Black Hole Mass Function Across Cosmic Times II. Heavy Seeds and (Super)Massive Black Holes

This is the second paper in a series aimed at modeling the black hole (BH) mass function, from the stellar to the (super)massive regime. In the present work we focus on (super)massive BHs and provide an ab-initio computation of their mass function across cosmic times. We consider two main mechanisms to grow the central BH, that are expected to cooperate in the high-redshift star-forming progenitors of local massive galaxies. The first is the gaseous dynamical friction process, that can cause the migration toward the nuclear regions of stellar-mass BHs originated during the intense bursts of star formation in the gas-rich host progenitor galaxy, and the buildup of a central heavy BH seed $M_\bullet\sim 10^{3-5}\, M_\odot$ within short timescales $\lesssim$ some $10^7$ yr. The second mechanism is the standard Eddington-type gas disk accretion onto the heavy BH seed, through which the central BH can become (super)massive $M_\bullet\sim 10^{6-10}\, M_\odot$ within the typical star-formation duration $\lesssim 1$ Gyr of the host. We validate our semi-empirical approach by reproducing the observed redshift-dependent bolometric AGN luminosity functions and Eddington ratio distributions, and the relationship between the star-formation and the bolometric luminosity of the accreting central BH. We then derive the relic (super)massive BH mass function at different redshifts via a generalized continuity equation approach, and compare it with present observational estimates. Finally, we reconstruct the overall BH mass function from the stellar to the (super)massive regime, over more than ten orders of magnitudes in BH mass.

preprint2022arXiv

The far-infrared/radio correlation for a sample of strongly lensed dusty star-forming galaxies detected by Herschel

We investigate the radio-far infrared (FIR) correlation for a sample of $28$ bright high-redshift ($1 \lesssim z \lesssim 4$) star-forming galaxies selected in the FIR from the Herschel-ATLAS fields as candidates to be strongly gravitationally lensed. The radio information comes either from high sensitivity dedicated ATCA observations at $2.1$ GHz or from cross-matches with the FIRST survey at $1.4$ GHz. By taking advantage of source brightness possibly enhanced by lensing magnification, we identify a weak evolution with redshift out to $z\lesssim 4$ of the FIR-to-radio luminosity ratio $q_{\rm FIR}$. We also find that the $q_{\rm FIR}$ parameter as a function of the radio power $L_{1.4\,\rm GHz}$ displays a clear decreasing trend, similarly to what is observed for optically/radio selected lensed quasars found in literature, yet covering a complementary region in the $q_{\rm FIR}-L_{1.4\,\rm GHz}$ diagram. We interpret such a behavior in the framework of an in-situ galaxy formation scenario, as a result of the transition from an early dust-obscured star-forming phase (mainly pinpointed by our FIR selection) to a late radio-loud quasar phase (preferentially sampled by the optical/radio selection).

preprint2022arXiv

The impact of very massive stars on the chemical evolution of extremely metal-poor galaxies

In recent observations of extremely metal-poor low-mass starburst galaxies, almost solar Fe/O ratios are reported, despite N/O ratios consistent with the low metallicity. We investigate if the peculiar Fe/O ratios can be a distinctive signature of an early enrichment produced by very massive objects dying as Pair-Instability Supernovae (PISN). We run chemical evolution models with yields that account for the contribution by PISN. We use both the recent non-rotating stellar yields from Goswami et al. 2021, and new yields from rotating very massive stars calculated on purpose in this work. We also search for the best initial mass function (IMF) that is able to reproduce the observations. We can reproduce the observations by adopting a bi-modal IMF and by including an initial burst of rotating very massive stars. Only with a burst of very massive stars can we reproduce the almost solar Fe/O ratios at the estimated young ages. We also confirm that rotation is absolutely needed to concomitantly reproduce the observed N/O ratios. These results stress the importance of very massive stars in galactic chemical evolution studies and strongly support a top-heavy initial mass function in the very early evolutionary stages of metal poor starburst galaxies.

preprint2022arXiv

Tomography-based observational measurements of the halo mass function via the submillimeter magnification bias

Aims. The main goal of this paper is to derive observational constraints on the halo mass fuction (HMF) by performing a tomographic analysis of the magnification bias signal on a sample of background submillimeter galaxies. The results can then be compared with those from a non-tomographic study. Methods. We measure the cross-correlation function between a sample of foreground GAMA galaxies with spectroscopic redshifts in the range $0.1 < z < 0.8$ (and divided up into four bins) and a sample of background submillimeter galaxies from H-ATLAS with photometric redshifts in the range $1.2 < z < 4.0$. We model the weak lensing signal within the halo model formalism and carry out a Markov chain Monte Carlo algorithm to obtain the posterior distribution of all HMF parameters, which we assume to follow the Sheth and Tormen (ST) three-parameter and two-parameter fits. Results. While the observational constraints on the HMF from the non-tomographic analysis are not stringent, there is a remarkable improvement in terms of uncertainty reduction when tomography is adopted. Moreover, with respect to the traditional ST triple of values from numerical simulations, the results from the three-parameter fit predict a higher number density of halos at masses below $10^{12}M_{\odot}/h$ at 95% credibility. The two-parameter fit yields even more restricting results, with a larger number density of halos below $10^{13}M_{\odot}/h$ and a lower one above $10^{14}M_{\odot}/h$, this time at more than 3$σ$ credibility. Our results are therefore in disagreement with the standard N-body values for the ST fit at 2$σ$ and 3$σ$, respectively.

preprint2022arXiv

Warm molecular and ionized gas kinematics in the Type-2 quasar J0945+1737

We analyze Near-Infrared Integral Field Spectrograph (NIFS) observations of the type-2 quasar (QSO2) SDSS J094521.33+173753.2 to investigate its warm molecular and ionized gas kinematics. This QSO2 has a bolometric luminosity of 10$^{45.7}$ erg s$^{-1}$ and a redshift of z = 0.128. The K-band spectra provided by NIFS cover a range of 1.99-2.40 $μ$m where low-ionization (Pa$α$ and Br$δ$), high ionization ([S XI]$λ$1.920 $μ$m and [Si~VI]$λ$1.963 $μ$m) and warm molecular lines (from H$_2$ 1-0S(5) to 1-0S(1)) are detected, allowing us to study the multi-phase gas kinematics. Our analysis reveals gas in ordinary rotation in all the emission lines detected and also outflowing gas in the case of the low- and high-ionization emission lines. In the case of the nuclear spectrum, which corresponds to a circular aperture of 0.3\arcsec~(686 pc) in diameter, the warm molecular lines can be characterized using a single Gaussian component of full width at half maximum (FWHM)= 350-400 km s$^{-1}$, while Pa$α$, Br$δ$, and [Si~VI] are best fitted with two blue-shifted Gaussian components of FWHM$\sim$800 and 1700 km s$^{-1}$, in addition to a narrow component of $\sim$300 km s$^{-1}$. We interpret the blue-shifted broad components as outflowing gas, which reaches the highest velocities, of up to $-$840 km s$^{-1}$, in the south-east direction (PA$\sim$125$^{\circ}$), extending up to a distance of $\sim$3.4 kpc from the nucleus. The ionized outflow has a maximum mass outflow rate of $\dot{\text{M}}_{\text{out, max}}$=42-51 M$_\odot$ yr$^{-1}$, and its kinetic power represents 0.1$\%$ of the quasar bolometric luminosity.

preprint2021arXiv

A direct and robust method to observationally constrain the halo mass function via the submillimeter magnification bias: Proof of concept

Aims. The main purpose of this work is to provide a method to derive tabulated observational constraints on the halo mass function (HMF) by studying the magnification bias effect on high-redshift submillimeter galaxies. Under the assumption of universality, we parametrize the HMF according to two traditional models, namely the Sheth and Tormen (ST) and Tinker fits and assess their performance in explaining the measured data within the Λ cold dark matter (ΛCDM) model. We also study the potential influence of the halo occupation distribution (HOD) parameters in this analysis and discuss two important aspects regarding the HMF parametrization. Methods. We measure the cross-correlation function between a foreground sample of GAMA galaxies with redshifts in the range $0.2<z<0.8$ and a background sample of H-ATLAS galaxies with redshifts in the range $1.2<z<4.0$ and carry out an MCMC algorithm to check this observable against its mathematical prediction within the halo model formalism. Results. If all HMF parameters are assumed to be positive, the ST fit only seems to fully explain the measurements by forcing the mean number of satellite galaxies in a halo to increase substantially from its prior mean value. The Tinker fit, on the other hand, provides a robust description of the data without relevant changes in the HOD parameters, but with some dependence on the prior range of two of its parameters. When the normalization condition for the HMF is dropped and we allow negative values of the $p_1$ parameter in the ST fit, all the involved parameters are better determined, unlike the previous models, thus deriving the most general HMF constraints. While all cases are in agreement with the traditional fits within the uncertainties, the last one hints at a slightly higher number of halos at intermediate and high masses, raising the important point of the allowed parameter range.

preprint2021arXiv

Galaxy cluster mass density profile derived using the submillimetre galaxies magnification bias

In this work we want to study the average mass density profile of tens to hundreds of clusters of galaxies acting as lenses that produce a magnification bias on the SMGs, and to estimate their associated masses and concentrations for different richness ranges. The background sample is composed of SMGs observed by Herschel with 1.2<z<4.0 (mean redshift at ~2.3) while the foreground sample is made up of galaxy clusters extracted from the SDSS III with photometric redshifts of 0.05< z< 0.8 (mean redshift at ~0.38). Measurements are obtained by stacking the SMG--cluster pairs to estimate the cross-correlation function using the Davis-Peebles estimator. This methodology allows us to derive the mass density profile for a wide range of angular scales, ~2-250 arcsec or ~10-1300 kpc for z=0.38, with a high radial resolution. We find that It is impossible to fit the data with a single mass density profile at all scales. As for the outer part, the estimated average masses increase from $M_{200c}=5.8$ to $51.5\times 10^{13} M_\odot$ and the concentration parameter from C=0.74 to 1.74. In the small-scale regions, the obtained average masses fluctuate around $M_{200c}=3-4 \times 10^{13}M_\odot$ with average C~4. The total average masses are in perfect agreement with the M-R relationship estimated from the cluster catalogue. While the estimated average C values of the central galactic halos are in agreement with traditional M-C relationships, we find low concentrations for the outer part. Moreover, C decrease for lower R values, probably indicating that the group of galaxies cannot be considered to be relaxed systems. Finally, we notice a systematic lack of signal at the transition between the dominance of the cluster halo and the central galactic halo (~100 kpc). This feature is also present in previous studies using different catalogues and/or methodologies, but is never discussed.

preprint2020arXiv

Cosmological constraints with the sub-millimetre galaxies Magnification Bias after large scale bias corrections

The study of the magnification bias produced on high-redshift sub-millimetre galaxies by foreground galaxies through the analysis of the cross-correlation function was recently demonstrated as an interesting independent alternative to the weak-lensing shear as a cosmological probe. In the case of the proposed observable, most of the cosmological constraints mainly depend on the largest angular separation measurements. Therefore, we aim to study and correct the main large-scale biases that affect foreground and background galaxy samples to produce a robust estimation of the cross-correlation function. Then we analyse the corrected signal to derive updated cosmological constraintsWe measured the large-scale, bias-corrected cross-correlation functions using a background sample of H-ATLAS galaxies with photometric redshifts > 1.2 and two different foreground samples (GAMA galaxies with spectroscopic redshifts or SDSS galaxies with photometric ones, both in the range 0.2 < z < 0.8). These measurements are modelled using the traditional halo model description that depends on both halo occupation distribution and cosmological parameters. We then estimated these parameters by performing a Markov chain Monte Carlo under multiple scenarios to study the performance of this observable and how to improve its results. After the large-scale bias corrections, we obtain only minor improvements with respect to the previous magnification bias results, mainly confirming their conclusions: a lower bound on $Ω_m > 0.22$ at $95\%$ C.L. and an upper bound $σ_8 < 0.97$ at $95\%$ C.L. (results from the $z_{spec}$ sample). However, by combining both foreground samples into a simplified tomographic analysis, we were able to obtain interesting constraints on the $Ω_m$-$σ_8$ plane as follows: $Ω_m= 0.50_{- 0.20}^{+ 0.14}$ and $σ_8= 0.75_{- 0.10}^{+ 0.07}$ at 68\% CL.

preprint2020arXiv

Cosmology with the submillimetre galaxies magnification bias: Proof of concept

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 $α$), but also some of the cosmological ones ($Ω_m$, $σ_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, $Ω_m$ and $H_0$ cannot be well constrained. However, we can set a lower limit of >0.24 at 95\% confidence level (CL) on $Ω_m$ and we see a slight trend towards $H_0>70$ values. For our constraints on $σ_8$ we obtain only a tentative peak around 0.75, but an interesting upper limit of $σ_8\lesssim 1$ at 95\% CL. We also study the possibility to derive better constraints by imposing more restrictive priors on the astrophysical parameters.

preprint2020arXiv

Growth of Supermassive Black Hole Seeds in ETG Star-Forming Progenitors: Multiple Merging of Stellar Compact Remnants via Gaseous Dynamical Friction and Gravitational Wave Emission

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

preprint2020arXiv

IRAM 30m-EMIR Redshift Search of z = 3-4 Lensed Dusty Starbursts selected from the HerBS sample

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

preprint2020arXiv

Knocking on giants' doors: I. The evolution of the dust-to-stellar mass ratio in distant dusty galaxies

The dust-to-stellar mass ratio ($M_{\rm dust}$/$M_{\rm \star}$) is a crucial yet poorly constrained quantity to understand the production mechanisms of dust, metals and stars in galaxy evolution. In this work we explore and interpret the nature of $M_{\rm dust}$/$M_{\rm \star}$ in 300 massive ($M_{\star}>10^{10}M_{\odot}$), dusty star-forming galaxies detected with ALMA up to $z\approx5$. We find that $M_{\rm dust}$/$M_{\rm \star}$ evolves with redshift, stellar mass, specific SFR and integrated dust size, differently for main sequence and starburst galaxies. In both galaxy populations $M_{\rm dust}$/$M_{\rm \star}$ rises until $z\sim2$ followed by a roughly flat trend towards higher redshifts. We show that the inverse relation between $M_{\rm dust}$/$M_{\rm \star}$ and $M_{\star}$ holds up to $z\approx5$ and can be interpreted as an evolutionary transition from early to late starburst phases. We demonstrate that $M_{\rm dust}$/$M_{\rm \star}$ in starbursts mirrors the increase in molecular gas fraction with redshift, and is enhanced in objects with the most compact dusty star-formation. The state-of-the-art cosmological simulation SIMBA broadly matches the evolution of $M_{\rm dust}$/$M_{\rm \star}$ in main sequence galaxies, but underestimates it in starbursts. The latter is found to be linked to lower gas-phase metallicities and longer dust growth timescales relative to data. Our data are well reproduced by analytical model that includes recipes for rapid metal enrichment, strongly suggesting that high $M_{\rm dust}$/$M_{\rm \star}$ is due to fast grain growth in metal enriched ISM. Our work highlights multifold benefits of using $M_{\rm dust}$/$M_{\rm \star}$ as a diagnostic tool for: (1) separating main sequence and starburst galaxies until $z\sim5$; (2) probing the evolutionary phases of dusty galaxies, and (3) refining the treatment of dust life cycle in simulations.

preprint2020arXiv

New Analytic Solutions for Galaxy Evolution II: Wind Recycling, Galactic Fountains and Late-Type Galaxies

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.

preprint2016arXiv

Another look to distortions of the CMB spectrum

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.

preprint2016arXiv

The Herschel-ATLAS: a sample of 500μm-selected lensed galaxies over 600 square degrees

We present a sample of 80 candidate strongly lensed galaxies with flux density above 100mJy at 500μ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μ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.

preprint2016arXiv

The Intragroup versus the Intracluster Medium

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.

preprint2016arXiv

The Quest for Dusty Star-forming Galaxies at High Redshift z>4

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.

preprint2015arXiv

Black Hole and Galaxy Coevolution from Continuity Equation and Abundance Matching

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

preprint2015arXiv

Cold or Warm? Constraining Dark Matter with Primeval Galaxies and Cosmic Reionization after Planck

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.

preprint2015arXiv

Cross-correlation between the CMB lensing potential measured by Planck and high-z sub-mm galaxies detected by the Herschel-ATLAS survey

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\,σ$ 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\,σ$. 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.

preprint2015arXiv

Predictions for surveys with the SPICA Mid-infrared Instrument

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$μ$m) and two silicate bands (in emission and in absorption) at 9.7$μ$m and 18.0$μ$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\,σ$ 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\,μ$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.

preprint2014arXiv

Entropy Flattening, Gas Clumping and Turbulence in Galaxy Clusters

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.

preprint2014arXiv

Exploring the relationship between black hole accretion and star formation with blind mid-/far-infrared spectroscopic surveys

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$σ$) $\sim$760 AGNs in [OIV]25.89$μ$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.

preprint2014arXiv

Herschel-ATLAS/GAMA: SDSS cross-correlation induced by weak lensing

We report a highly significant ($>10σ$) spatial correlation between galaxies with $S_{350μ\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 ($μ<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.

preprint2014arXiv

Probing the evolution of galaxy clusters with the SZ Effect

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.

preprint2014arXiv

Radio Observations of Star Forming Galaxies in the SKA era

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.

preprint2014arXiv

Statistics of Dark Matter Halos in the Excursion Set Peak Framework

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.

preprint2014arXiv

SuperModel Analysis of Abell 1246 and J255: on the Evolution of Galaxy Clusters from High to Low Entropy States

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.

preprint2014arXiv

The Coevolution of Supermassive Black Holes and Massive Galaxies at High Redshift

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.

preprint2014arXiv

The Herschel Stripe 82 Survey (HerS): Maps and Early Catalog

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 $σ$ (including confusion noise). The maps and catalog are available at http://www.astro.caltech.edu/hers/

preprint2013arXiv

Dust and star-formation properties of a complete sample of local galaxies drawn from the Planck Early Release Compact Source Catalogue

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.

preprint2013arXiv

Exploring the early dust-obscured phase of galaxy formation with blind mid-/far-IR spectroscopic surveys

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.

preprint2013arXiv

Forecasts on the contamination induced by unresolved point sources in primordial non-Gaussianity beyond Planck

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.

preprint2013arXiv

Gravitational Lens Models Based on Submillimeter Array Imaging of Herschel-selected Strongly Lensed Sub-millimeter Galaxies at z>1.5

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.

preprint2013arXiv

H-ATLAS: estimating redshifts of Herschel sources from sub-mm fluxes

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}

preprint2013arXiv

Herschel-ATLAS: deep HST/WFC3 imaging of strongly lensed submillimeter galaxies

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.

preprint2013arXiv

Statistics of Dark Matter Halos from the Excursion Set Approach

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.

preprint2013arXiv

The Astrophysics of the Intracluster Plasma

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

preprint2013arXiv

Turbulence in the SuperModel: Mass Reconstruction with Nonthermal Pressure for Abell 1835

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.

preprint2012arXiv

Clustering of sub-millimeter galaxies in a self-regulated baryon collapse model

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 $λ\ge 850\,μ$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\,μ$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.

preprint2012arXiv

Effective Models for Statistical Studies of Galaxy-Scale Gravitational Lensing

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.

preprint2012arXiv

Herschel-ATLAS: Planck sources in the Phase 1 fields

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.

preprint2012arXiv

Herschel-ATLAS: towards a sample of ~1000 strongly-lensed galaxies

While the selection of strongly lensed galaxies with 500μ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 μ>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.

preprint2012arXiv

Herschel-ATLAS: VISTA VIKING near-IR counterparts in the Phase 1 GAMA 9h data

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

preprint2012arXiv

Herschel-ATLAS/GAMA: a census of dust in optically selected galaxies from stacking at submillimetre wavelengths

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μ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μm background (4% at 350μm; 3% at 500μ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μ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.

preprint2012arXiv

The evolutionary connection between QSOs and SMGs: molecular gas in far-infrared luminous QSOs at z ~ 2.5

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.

preprint2011arXiv

A Grand Design for Galaxy Clusters: Connections and Predictions

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.

preprint2011arXiv

Herschel-ATLAS Galaxy Counts and High Redshift Luminosity Functions: The Formation of Massive Early Type Galaxies

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.

preprint2011arXiv

Self-Similar Dynamical Relaxation of Dark Matter Halos in an Expanding Universe

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.

preprint2011arXiv

The Intracluster Plasma: a Universal Pressure Profile?

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.

preprint2010arXiv

Cosmic Evolution of Size and Velocity Dispersion for Early Type Galaxies

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

preprint2010arXiv

Dark Matter Halos: The Dynamical Basis of Effective Empirical Models

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

preprint2010arXiv

Herschel-ATLAS: blazars in the SDP field

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$μ$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 $γ$-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.

preprint2010arXiv

Non Thermal Support for the Outer Intracluster Medium

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.

preprint2010arXiv

Probing the Astrophysics of Cluster Outskirts

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.

preprint2010arXiv

SAGACE: the Spectroscopic Active Galaxies And Clusters Explorer

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.

preprint2009arXiv

Dark Matter Equilibria in Galaxies and Galaxy Systems

[abridged] In the dark matter (DM) halos embedding galaxies and galaxy systems the `entropy' K = σ^2 / ρ^{2/3} (a quantity that combines the radial velocity dispersion σwith the density ρ) is found from intensive N-body simulations to follow a powerlaw run K ~ r^α throughout the halos' bulk, with αaround 1.25. Taking up from phenomenology just that α~ 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 α-profiles, marked by the entropy slope αitself, and by the maximal gravitational pull κ_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 αto within the narrow range 1.25-1.29 from galaxies to galaxy systems. Our range of α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 α-profiles with parameters derived from the Jeans equation. We finally stress how our findings and predictions as to αand κ_crit contribute to understand hitherto unsolved issues concerning the fundamental structure of DM halos.

preprint2009arXiv

Galaxy Clusters, a Novel Look at Diffuse Baryons Withstanding Dark Matter Gravity

[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^αand α~ 1.25 - 1.3 applying from groups to rich clusters, yields our radial α-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 β-models, and the `mirror' model with T(r) ~ σ^2(r) suitable for NCC and CC clusters, respectively; these highlight how the ICP temperature T(r) tends to mirror the DM velocity dispersion σ^2(r) away from entropy injections. Finally, we discuss how the Supermodel connects information derived from X-ray and gravitational lensing observations.

preprint2009arXiv

Gamma Rays from Annihilations at the Galactic Center in a Physical Dark Matter Distribution

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 ρ(r) and the velocity dispersion σ^2_r(r) expressed by the powerlaw behavior of the DM `entropy' K= σ_r^2/ρ^{2/3} ~ r^αwith α~ 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 `α-profiles'. We find the one with α=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 α-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.

preprint2009arXiv

Supermodel Analysis of Galaxy Clusters

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

preprint2009arXiv

The Herschel ATLAS

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.

preprint2009arXiv

Two phase galaxy formation: The Gas Content of Normal Galaxies

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.

preprint2009arXiv

Two-Phase Galaxy Formation

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.

preprint2007arXiv

The Baryonic and Dark Matter Properties of High Redshift Gravitationally Lensed Disk Galaxies

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.

preprint2007arXiv

The Role of the Dust in Primeval Galaxies: A Simple Physical Model for Lyman Break Galaxies and Lyman Alpha Emitters

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.

preprint2006arXiv

New Relationships between Galaxy Properties and Host Halo Mass, and the Role of Feedbacks in Galaxy Formation

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.

preprint2006arXiv

Quasar Luminosity Functions from Joint Evolution of Black Holes and Host Galaxies

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.

preprint2006arXiv

Searching for AGN-driven Shocks in Galaxy Clusters

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.