The detection of signals hidden in noise is one of the oldest and common problems in astronomy. Various solutions have been proposed in the past such as the parametric approaches based on the least-squares fit of theoretical templates or the non-parametric techniques as the phase-folding method. Most of them, however, are suited only for signals with specific time evolution. For generic signals the spectral approach based on the periodogram is potentially the most effective. In astronomy the main problem in working with the periodogram is that often the sampling of the signals is irregular. This complicates its efficient computation (the fast Fourier transform cannot be directly used) but overall the determination of its statistical characteristics. The Lomb-Scargle periodogram (LSP) provides a solution to this last important issue, but its main drawback is the assumption of a very specific model of the datawhich is not correct for most of the practical applications. These issues are not always considered in literature with theoretical and practical consequences of no easy solution. Moreover, apart from pathological samplings, it is common believe that the LSP and the classical periodogram (CP) usually provide almost identical results. In general, this is true but here it is shown that there are situations where the LSP is less effective than the CP in the detection of signals in noise. There are no compelling reasons, therefore, to use the LSP instead of the CP which is directly connected to the correlation function of the observed signal with the sinusoidal functions at the various frequencies of interest.

Fixing the relationship among a set of experimental quantities is a fundamental issue in many scientific disciplines. In the two-dimensional case, the classical approach is to compute the linear correlation coefficient from a scatterplot. This method, however, implicitly assumes a linear relationship between the variables. Such assumption is not always correct. With the use of the partial correlation coefficients, an extension to the multi-dimensional case is possible. However, the problem of the assumed mutual linear relationship among the variables still remains. A relatively recent approach which permits to avoid this problem is modeling the joint probability density function (PDF) of the data with copulas. These are functions which contain all the information on the relationship between two random variables. Although in principle this approach can work also with multi-dimensional data, theoretical as well computational difficulties often limit its use to the two-dimensional case. In this paper, we consider an approach, based on so-called vine copulas, which overcomes this limitation and at the same time is amenable to a theoretical treatment and feasible from the computational point of view. We apply this method to published data on the near-IR and far-IR luminosities and atomic and molecular masses of the Herschel Reference Sample. We determine the relationship among the luminosities and gas masses and show that the far-IR luminosity can be considered as the key parameter which relates all the other three galaxy properties. Once removed from the 4D relation, the residual relation among the other three is negligible. This may be interpreted as that the correlation between the gas masses and near-IR luminosity is driven by the far-IR luminosity, likely by the star-formation activity of the galaxy.

We measure the 850-$\mu$m source densities of 46 candidate protoclusters selected from the Planck High-z catalogue (PHz) and the Planck Catalogue of Compact Sources (PCCS) that were followed up with Herschel-SPIRE and SCUBA-2. This paper aims to search for overdensities of 850-$\mu$m sources in order to select the fields that are most likely to be genuine protoclusters. Of the 46 candidate protoclusters, 25 have significant overdensities ($>$5 times the field counts), 11 have intermediate overdensities (3--5 times the field counts) and 10 have no overdensity ($<$3 times the field counts) of 850-$\mu$m sources. We find that the enhanced number densities are unlikely to be the result of sample variance. Compared with the number counts of another sample selected from Planck's compact source catalogues, this [PHz+PCCS]-selected sample has a higher fraction of candidate protoclusters with significant overdensities, though both samples show overdensities of 850-$\mu$m sources above intermediate level. Based on the estimated star-formation rate densities (SFRDs), we suggest that both samples can efficiently select protoclusters with starbursting galaxies near the redshift at which the global field SFRD peaks ($2 < z < 3$). Based on the confirmation of overdensities found here, future follow-up observations on other PHz targets may greatly increase the number of genuine DSFG-rich clusters/protoclusters.

We present a sub-kpc resolved study of the interstellar medium properties in SDP.81, a z=3.042 strongly gravitationally lensed dusty star-forming galaxy, based on high-resolution, multi-band ALMA observations of the FIR continuum, CO ladder and the [CII] line. Using a visibility-plane lens modelling code, we achieve a median source-plane resolution of ~200 pc. We use photon-dominated region (PDR) models to infer the physical conditions - far-UV field strength, density, and PDR surface temperature - of the star-forming gas on 200-pc scales, finding a FUV field strength of ~10^3-10^4 G0, gas density of ~10^5 cm^-3 and cloud surface temperatures up to 1500 K, similar to those in the Orion Trapezium region. The [CII] emission is significantly more extended than that FIR continuum: ~50 per cent of [CII] emission arises outside the FIR-bright region. The resolved [CII]/FIR ratio varies by almost 2 dex across the source, down to ~2x10^-4 in the star-forming clumps. The observed [CII]/FIR deficit trend is consistent with thermal saturation of the C+ fine-structure level occupancy at high gas temperatures. We make the source-plane reconstructions of all emission lines and continuum data publicly available.

We present the first detection of diffuse dust in the intra-cluster medium of the Virgo cluster out to $\sim$0.4 virial radii, and study the radial variation of its properties on a radial scale of the virial radius. Analysing near-UV - $i$ colours for a sample of $\sim12000$ background galaxies with redshifts $0.02 < z < 0.8$, we find significant colour reddening and relate it to variation in $E(B-V)$ values. The $E(B-V)$ mean profile shows a dust component characterised by an average reddening $E(B-V)\sim0.042 \pm 0.004$ mag within 1.5 degrees ($\sim0.3\, r_{vir}$) from the cluster centre. Assuming a Large Magellanic Cloud extinction law, we derive an average visual extinction $A_{V} = 0.14\pm 0.01$ for a total dust mass, $M_{d} = 2.5\pm0.2\times10^{9}M_{\odot}$, hence a dust-to-gas mass ratio $M_{d}/M_{g} = 3.0\pm 0.3 \times 10^{-4}$. Based on the upper limits on the flux density $\mathrm{I_{250\mu m} = 0.1\, MJy sr^{-1}} $ derived from $Herschel$ data, we estimate an upper limit for the dust temperature of $T_{d} \sim 10\, K$. However, similar densities can be obtained with dust at higher temperatures with lower emissivities. The Virgo cluster has diffuse dust in its intra-cluster medium characterised by different physical properties as those characterising the Milky Way dust. The diffuse dust in Virgo is transported into the cluster space through similar phenomena (stripping) as those building up the optical intra-cluster light, and it constitutes an additional cooling agent of the cluster gas.

A prominent jet-driven outflow of CO(2-1) molecular gas is found along the kinematic minor axis of the Seyfert 2 galaxy ESO 420-G13, at a distance of $340-600\, \rm{pc}$ from the nucleus. The wind morphology resembles a characteristic funnel shape, formed by a highly collimated filamentary emission at the base, likely tracing the jet propagation through a tenuous medium, until a bifurcation point at $440\, \rm{pc}$ where the jet hits a dense molecular core and shatters, dispersing the molecular gas into several clumps and filaments within the expansion cone. We also trace the jet in ionised gas within the inner $\lesssim 340\, \rm{pc}$ using the [NeII]$_{\rm 12.8 \mu m}$ line emission, where the molecular gas follows a circular rotation pattern. The wind outflow carries a mass of $\sim 8 \times 10^6\, \rm{M_\odot}$ at an average wind projected speed of $\sim 160\, \rm{km\,s^{-1}}$, which implies a mass outflow rate of $\sim 14\, \rm{M_\odot\,yr^{-1}}$. Based on the structure of the outflow and the budget of energy and momentum, we discard radiation pressure from the active nucleus, star formation, and supernovae as possible launching mechanisms. ESO 420-G13 is the second case after NGC 1377 where the presence of a previously unknown jet is revealed due to its interaction with the interstellar medium, suggesting that unknown jets in feeble radio nuclei might be more common than expected. Two possible jet-cloud configurations are discussed to explain the presence of an outflow at such distance from the AGN. The outflowing gas will likely not escape, thus a delay in the star formation rather than quenching is expected from this interaction, while the feedback effect would be confined within the central few hundred parsecs of the galaxy.

Dusty star-forming galaxies (DSFGs) detected at $z > 4$ provide important examples of the first generations of massive galaxies. However, few examples with spectroscopic confirmation are currently known, with Hershel struggling to detect significant numbers of $z > 6$ DSFGs. NGP6_D1 is a bright 850 $\mu m$ source (12.3 $\pm$ 2.5 mJy) with no counterparts at shorter wavelengths (a SPIRE dropout). Interferometric observations confirm it is a single source, with no evidence for any optical or NIR emission, or nearby likely foreground lensing sources. No $>3\sigma$ detected lines are seen in both LMT RSR and IRAM 30m EMIR spectra of NGP6_D1 across 32 $GHz$ of bandwidth despite reaching detection limits of $\sim 1 mJy/500 km~s^{-1}$, so the redshift remains unknown. Template fitting suggests that NGP6_D1 is most likely between $z = 5.8$ and 8.3. SED analysis finds that NGP6_D1 is a ULIRG, with a dust mass $\sim 10^8$ - $10^9$ $M_{\odot}$ and a SFR of $\sim$ 500 $M_{\odot}~yr^{-1}$. We place upper limits on the gas mass of NGP6_D1 of $M_{H2}$ $ < (1.1~\pm~3.5) \times 10^{11}$ $M_{\odot}$, consistent with a gas-to-dust ratio of $\sim$ 100 - 1000. We discuss the nature of NGP6_D1 in the context of the broader submm population, and find that comparable SPIRE dropouts account for $\sim$ 20% of all SCUBA-2 detected sources, but with a similar flux density distribution to the general population.

We present SCUBA-2 850-$\mu$m observations of 13 candidate starbursting protoclusters selected using Planck and Herschel data. The cumulative number counts of the 850-$\mu$m sources in 9/13 of these candidate protoclusters show significant overdensities compared to the field, with the probability $<$10$^{-2}$ assuming the sources are randomly distributed in the sky. Using the 250-, 350-, 500- and 850-$\mu$m flux densities, we estimate the photometric redshifts of individual SCUBA-2 sources by fitting spectral energy distribution (SED) templates with an MCMC method. The photometric redshift distribution, peaking at $2<z<3$, is consistent with that of known $z>2$ protoclusters and the peak of the cosmic star-formation rate density (SFRD). We find that the 850-$\mu$m sources in our candidate protoclusters have infrared luminosities of $L_{\mathrm{IR}}\gtrsim$10$^{12}L_{\odot}$ and star-formation rates of SFR=(500-1,500)$M_{\odot}$yr$^{-1}$. By comparing with results in the literature considering only Herschel photometry, we conclude that our 13 candidate protoclusters can be categorised into four groups: six of them being high-redshift starbursting protoclusters, one being a lower-redshift cluster/protocluster, three being protoclusters that contain lensed DSFG(s) or are rich in 850-$\mu$m sources, and three regions without significant Herschel or SCUBA-2 source overdensities. The total SFRs of the candidate protoclusters are found to be comparable or higher than those of known protoclusters, suggesting our sample contains some of the most extreme protocluster population. We infer that cross-matching Planck and Herschel data is a robust method for selecting candidate protoclusters with overdensities of 850-$\mu$m sources.

We report Atacama Large Millimetre/submillimeter Array (ALMA) Cycle-5 follow-up observations of two candidate [CII] emitters at z ~ 6 in the ALMA deep field in SSA22 (ADF22). The candidates were detected blindly in a Cycle-2 ALMA survey covering ~ 5 square arcmins, with a single tuning, along with two CO lines associated with galaxies at lower redshifts. Various tests suggested at least one of the two > 6-sigma [CII] candidates should be robust (Hayatsu et al. 2017). Nevertheless, our new, deeper observations recover neither candidate, demonstrating a higher contamination rate than expected. The cause of the spurious detections is under investigation but at present it remains unclarified.

The matched filter (MF) represents one of the main tools to detect signals from known sources embedded in the noise. In the Gaussian case the noise is assumed to be the realization of a Gaussian random field (GRF). The most important property of the MF, the maximization of the probability of detection subject to a constant probability of false detection or false alarm (PFA), makes it one of the most popular techniques. However, the MF technique relies upon the a priori knowledge of the number and the position of the searched signals in the GRF which usually are not available. A typical way out is to assume that the position of a signal coincides with one of the peaks in the matched filtered data. A detection is claimed when the probability that a given peak is due only to the noise (i.e. the PFA) is smaller than a prefixed threshold. In this case the probability density function (PDF) of the amplitudes has to be used for the computation of the PFA, which is different from the Gaussian. Moreover, the probability that a detection is false depends on the number of peaks present in the filtered GRF, the greater the number of peaks in a GRF, the higher the probability of peaks due to the noise that exceed the detection threshold. If not taken into account, the PFA can be severely underestimated. Many solutions proposed to this problem are non-parametric hence not able to exploit all the available information. This limitation has been overcome by means of two efficient parametric approaches, one based on the PDF of the peak amplitudes of a smooth and isotropic GRF whereas the other uses the Gumbel distribution (the asymptotic PDF of the corresponding extreme). Simulations and ALMA maps show that, although the two methods produce almost identical results, the first is more flexible and allows us to check the reliability of the detection procedure.

We present the science case for mapping several thousand galaxy (proto)clusters at z=1-10 with a large aperture single dish sub-mm facility, producing a high-redshift counterpart to local large surveys of rich clusters like the well-studied Abell catalogue. Principal goals of a large survey of distant clusters are the evolution of galaxy clusters over cosmic time and the impact of environment on the evolution and formation of galaxies. To make a big leap forward in this emerging research field, the community would benefit from a large-format, wide-band, direct-detection spectrometer (e.g., based on MKID technology), covering a wide field of ~1 square degree and a frequency coverage from 70 to 700 GHz.

The cycling of baryons in and out of galaxies is what ultimately drives galaxy formation and evolution. The circumgalactic medium (CGM) represents the interface between the interstellar medium and the cosmic web, hence its properties are directly shaped by the baryon cycle. Although traditionally the CGM is thought to consist of warm and hot gas, recent breakthroughs are presenting a new scenario according to which an important fraction of its mass may reside in the cold atomic and molecular phase. This would represent fuel that is readily available for star formation, with crucial implications for feeding and feedback processes in galaxies. However, such cold CGM, especially in local galaxies where its projected size on sky is expected to be of several arcminutes, cannot be imaged by ALMA due to interferometric spatial scale filtering of large-scale structures. We show that the only way to probe the multiphase CGM including its coldest component is through a large (e.g. 50-m) single dish (sub-)mm telescope.

The Sunyaev-Zeldovich (SZ) effect was first predicted nearly five decades ago, but has only recently become a mature tool for performing high resolution studies of the warm and hot ionized gas in and between galaxies, groups, and clusters. Galaxy groups and clusters are powerful probes of cosmology, and they also serve as hosts for roughly half of the galaxies in the Universe. In this white paper, we outline the advances in our understanding of thermodynamic and kinematic properties of the warm-hot universe that can come in the next decade through spatially and spectrally resolved measurements of the SZ effects. Many of these advances will be enabled through new/upcoming millimeter/submillimeter (mm/submm) instrumentation on existing facilities, but truly transformative advances will require construction of new facilities with larger fields of view and broad spectral coverage of the mm/submm bands.

Why do we write this note? It is erroneous to pretend to extract physical information from the experimental light curves (time series) of astrophysical systems by means of linear stochastic differential equations (LSDE). In general, the time evolution of these systems is governed by a set of nonlinear differential equations. Hence, the LSDEs are not suitable to model their dynamics. In spite of this, recently the LSDEs have been proposed as tools for the analysis of AGN light curves. Their use in this context seems to be dictated by their simplicity rather than by a real physical argument. We stress in this note that the correct approach to the analysis of signals coming from systems with nonlinear dynamics is to tackle the problem using methodologies in well defined physical contexts.

We have observed the giant elliptical galaxy M87 during the Virgo Environmental Survey Tracing Galaxy Evolution (VESTIGE), a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster carried out with MegaCam at the Canada French Hawaii Telescope (CFHT). The deep narrow-band image confirmed the presence of a filament of ionised gas extending up to ~ 3 kpc in the north-western direction and ~ 8 kpc to the south-east, with a couple of plumes of ionised gas, the weakest of which, at ~ 18 kpc from the nucleus, was previously unknown. The analysis of deep optical images taken from the NGVS survey confirms that this gas filament is associated with dust seen in absorption which is now detected up to ~ 2.4 kpc from the nucleus. We have also analysed the physical and kinematical properties of the ionised gas filament using deep IFU MUSE data covering the central 4.8 x 4.8 kpc^2 of the galaxy. The spectroscopic data confirms a perturbed kinematics of the ionised gas, with differences in velocity of ~ 700-800 km s^-1 on scales of < 1 kpc. The analysis of 2D diagnostic diagrams and the observed relationship between the shock-sensitive [OI]/Halpha line ratio and the velocity dispersion of the gas suggest that the gas is shock-ionised.

We present measurements of the angular correlation function of sub-millimeter (sub-mm) galaxies (SMGs) identified in four out of the five fields of the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) - GAMA-9h, GAMA-12h, GAMA-15h and NGP - with flux densities $S_{250\mu m}$>30 mJy at 250 {\mu}m. We show that galaxies selected at this wavelength trace the underlying matter distribution differently at low and high redshifts. We study the evolution of the clustering finding that at low redshifts sub-mm galaxies exhibit clustering strengths of $r_0$ $\sim$ 2 - 3 $h^{-1}$ Mpc, below z < 0.3. At high redshifts, on the other hand, we find that sub-mm galaxies are more strongly clustered with correlation lengths $r_0$ = 8.1 $\pm$ 0.5, 8.8 $\pm$ 0.8 and 13.9 $\pm$ 3.9 $h^{-1}$Mpc at z = 1 - 2, 2 - 3 and 3 - 5, respectively. We show that sub-mm galaxies across the redshift range 1 < z < 5, typically reside in dark-matter halos of mass of the order of ~ $10^{12.5}$ - $10^{13.0}$ $h^{-1} \, M_{\odot}$ and are consistent with being the progenitors of local massive elliptical galaxies that we see in the local Universe.

We discuss the results of the relationships between the K-band and stellar mass, far-infrared luminosities, star formation rate, dust and gas masses of nearby galaxies computing the bivariate K-band Luminosity Function (BLF) and bivariate K-band Mass Function (BMF) of the Herschel Reference Survey (HRS), a volume-limited sample with full wavelength coverage. We derive the BLFs and BMFs from the K-band and stellar mass, far-infrared luminosities, star formation rate, dust and gas masses cumulative distributions using a copula method which is outlined in detail. The use of the bivariate computed taking into account the upper limits allows us to derive on a more solid statistical ground the relationship between the observed physical quantities. The analysis shows that the behaviour of the morphological (optically selected) subsamples is quite different. A statistically meaningful result can be obtained over the whole HRS sample only from the relationship between the K-band and the stellar mass, while for the remaining physical quantities (dust and gas masses, far-IR luminosity and star formation rate), the analysis is distinct for late-type (LT) and early-type galaxies (ETG). However, the number of ETGs is small to perform a robust statistical analysis, and in most of the case results are discussed only for the LTG subsample. The Luminosity and Mass Functions (LFs, MFs) of LTGs are generally dependent on the K-band and the various dependencies are discussed in detail. We are able to derive the corresponding LFs and MFs and compare them with those computed with other samples. Our statistical analysis allows us to characterise the HRS, that, although non homogeneously selected and partially biased towards low IR luminosities, may be considered as representative of the local LT galaxy population.

The matched filter (MF) is widely used to detect signals hidden within the noise. If the noise is Gaussian, its performances are well-known and describable in an elegant analytical form. The treatment of non-Gaussian noises is often cumbersome as in most cases there is no analytical framework. This is true also for Poisson noise which, especially in the low-number count regime, presents the additional difficulty to be discrete. For this reason in the past methods have been proposed based on heuristic or semi-heuristic arguments. Recently, an analytical form of the MF has been introduced but the computation of the probability of false detection or false alarm (PFA) is based on numerical simulations. To overcome this inefficient and time consuming approach we propose here an effective method to compute the PFA based on the saddle point approximation (SA). We provide the theoretical framework and support our findings by means of numerical simulations. We discuss also the limitations of the MF in practical applications.

Atomic Carbon can be an efficient tracer of the molecular gas mass, and when combined to the detection of high-J and low-J CO lines it yields also a sensitive probe of the power sources in the molecular gas of high redshift galaxies. The recently installed SEPIA5 receiver at the focus of the APEX telescope has opened up a new window at frequencies 159 - 211 GHz allowing the exploration of the Atomic Carbon in high-z galaxies, at previously inaccessible frequencies from the ground. We have targeted three gravitationally lensed galaxies at redshift of about 3 and conducted a comparative study of the observed high-J CO/CI ~ratios with well-studied nearby galaxies. Atomic Carbon (CI(2-1)) was detected in one of the three targets and marginally in a second, while in all three targets the $J=7\to6$ CO line is detected. The CO(7-6)/CI(2-1), CO(7-6)/CO(1-0) line ratios and the CO(7-6)/(far-IR continuum) luminosity ratio are compared to those of nearby objects. A large excitation status in the ISM of these high-z objects is seen, unless differential lensing unevenly boosts the CO line fluxes from the warm and dense gas more than the CO(1-0), CI(2-1), tracing a more widely distributed cold gas phase. We provide estimates of total molecular gas masses derived from the atomic Carbon and the Carbon monoxide CO(1-0), which within the uncertainties turn out to be equal.

Starburst galaxies at the peak of cosmic star formation are among the most extreme starforming engines in the universe, producing stars over ~100 Myr. The star formation rates of these galaxies, which exceed 100 $M_\odot$ per year, require large reservoirs of cold molecular gas to be delivered to their cores, despite strong feedback from stars or active galactic nuclei. Starburst galaxies are therefore ideal targets to unravel the critical interplay between this feedback and the growth of a galaxy. The methylidyne cation, CH$^+$, is a most useful molecule for such studies because it cannot form in cold gas without supra-thermal energy input, so its presence highlights dissipation of mechanical energy or strong UV irradiation. Here, we report the detection of CH$^+$(J=1-0) emission and absorption lines in the spectra of six lensed starburst galaxies at redshifts z~2.5. This line has such a high critical density for excitation that it is emitted only in very dense ($>10^5$ cm$^{-3}$) gas, and is absorbed in low-density gas. We find that the CH$^+$ emission lines, which are broader than 1000 km s$^{-1}$, originate in dense shock waves powered by hot galactic winds. The CH$^+$ absorption lines reveal highly turbulent reservoirs of cool ($T\sim 100$K), low-density gas, extending far outside (>10 kpc) the starburst cores (radii <1 kpc). We show that the galactic winds sustain turbulence in the 10 kpc-scale environments of the starburst cores, processing these environments into multi-phase, gravitationally bound reservoirs. However, the mass outflow rates are found to be insufficient to balance the star formation rates. Another mass input is therefore required for these reservoirs, which could be provided by on-going mergers or cold stream accretion. Our results suggest that galactic feedback, coupled jointly to turbulence and gravity, extends the starburst phase instead of quenching it.

The matched filter (MF) is one of the most popular and reliable techniques to the detect signals of known structure and amplitude smaller than the level of the contaminating noise. Under the assumption of stationary Gaussian noise, MF maximizes the probability of detection subject to a constant probability of false detection or false alarm (PFA). This property relies upon a priori knowledge of the position of the searched signals, which is usually not available. Recently, it has been shown that when applied in its standard form, MF may severely underestimate the PFA. As a consequence the statistical significance of features that belong to noise is overestimated and the resulting detections are actually spurious. For this reason, an alternative method of computing the PFA has been proposed that is based on the probability density function (PDF) of the peaks of an isotropic Gaussian random field. In this paper we further develop this method. In particular, we discuss the statistical meaning of the PFA and show that, although useful as a preliminary step in a detection procedure, it is not able to quantify the actual reliability of a specific detection. For this reason, a new quantity is introduced called the specific probability of false alarm (SPFA), which is able to carry out this computation. We show how this method works in targeted simulations and apply it to a few interferometric maps taken with the Atacama Large Millimeter/submillimeter Array (ALMA) and the Australia Telescope Compact Array (ATCA). We select a few potential new point sources and assign an accurate detection reliability to these sources.

We present a coherent multi-band modelling of the CO Spectral Energy Distribution of the local Seyfert Galaxy NGC7130 to assess the impact of the AGN activity on the molecular gas. We take advantage of all the available data from X-ray to the sub-mm, including ALMA data. The high-resolution (~0.2") ALMA CO(6-5) data constrain the spatial extension of the CO emission down to ~70 pc scale. From the analysis of the archival CHANDRA and NuSTAR data, we infer the presence of a buried, Compton-thick AGN of moderate luminosity, L_2-10keV ~ 1.6x10^{43} ergs-1. We explore photodissociation and X-ray-dominated regions (PDRs and XDRs) models to reproduce the CO emission. We find that PDRs can reproduce the CO lines up to J~6, however, the higher rotational ladder requires the presence of a separate source of excitation. We consider X-ray heating by the AGN as a source of excitation, and find that it can reproduce the observed CO Spectral Energy Distribution. By adopting a composite PDR+XDR model, we derive molecular cloud properties. Our study clearly indicates the capabilities offered by current-generation of instruments to shed light on the properties of nearby galaxies adopting state-of-the art physical modelling.

We present a coherent database of spectroscopic observations of far-IR fine-structure lines from the Herschel/PACS archive for a sample of 170 local AGN, plus a comparison sample of 20 starburst galaxies and 43 dwarf galaxies. Published Spitzer/IRS and Herschel/SPIRE line fluxes are included to extend our database to the full 10-600 $\mu m$ spectral range. The observations are compared to a set of CLOUDY photoionisation models to estimate the above physical quantities through different diagnostic diagrams. We confirm the presence of a stratification of gas density in the emission regions of the galaxies, which increases with the ionisation potential of the emission lines. The new [OIV]25.9$\mu m$/[OIII]88$\mu m$ vs [NeIII]15.6$\mu m$/[NeII]12.8$\mu m$ diagram is proposed as the best diagnostic to separate: $i)$ AGN activity from any kind of star formation; and $ii)$ low-metallicity dwarf galaxies from starburst galaxies. Current stellar atmosphere models fail to reproduce the observed [OIV]25.9$\mu m$/[OIII]88$\mu m$ ratios, which are much higher when compared to the predicted values. Finally, the ([NeIII]15.6$\mu m$ + [NeII]12.8$\mu m$)/([SIV]10.5$\mu m$ + [SIII]18.7$\mu m$) ratio is proposed as a promising metallicity tracer to be used in obscured objects, where optical lines fail to accurately measure the metallicity. The diagnostic power of mid- to far-infrared spectroscopy shown here for local galaxies will be of crucial importance to study galaxy evolution during the dust-obscured phase at the peak of the star formation and black-hole accretion activity ($1 < z < 4$). This study will be addressed by future deep spectroscopic surveys with present and forthcoming facilities such as JWST, ALMA, and SPICA.

We present new estimates of AGN accretion and star-formation luminosity in galaxies obtained for the local 12-$\mu$m sample of Seyfert galaxies (12MGS), by performing a detailed broad-band spectral energy distribution (SED) decomposition including the emission of stars, dust heated by star formation and a possible AGN dusty torus. Thanks to the availability of data from the X-rays to the sub-millimetre, we constrain and test the contribution of the stellar, AGN and star-formation components to the SEDs. The availability of Spitzer-IRS low resolution mid-infrared (mid-IR) spectra is crucial to constrain the dusty torus component at its peak wavelengths. The results of SED-fitting are also tested against the available information in other bands: the reconstructed AGN bolometric luminosity is compared to those derived from X-rays and from the high excitation IR lines tracing AGN activity like [Ne V] and [O IV]. The IR luminosity due to star-formation (SF) and the intrinsic AGN bolometric luminosity are shown to be strongly related to the IR line luminosity. Variations of these relations with different AGN fractions are investigated, showing that the relation dispersions are mainly due to different AGN relative contribution within the galaxy. Extrapolating these local relations between line and SF or AGN luminosities to higher redshifts, by means of recent Herschel galaxy evolution results, we then obtain mid- and far-IR line luminosity functions useful to estimate how many star-forming galaxies and AGN we expect to detect in the different lines at different redshifts and luminosities with future IR facilities (e.g., JWST, SPICA).

The detection reliability of weak signals is a critical issue in many astronomical contexts and may have severe consequences for determining number counts and luminosity functions, but also for optimising the use of telescope time in follow-up observations. Because of its optimal properties, one of the most popular and widely-used detection technique is the matched filter (MF). This is a linear filter designed to maximise the detectability of a signal of known structure that is buried in additive Gaussian random noise. In this work we show that in the very common situation where the number and position of the searched signals within a data sequence (e.g. an emission line in a spectrum) or an image (e.g. a point-source in an interferometric map) are unknown, this technique, when applied in its standard form, may severely underestimate the probability of false detection. This is because the correct use of the MF relies upon a-priori knowledge of the position of the signal of interest. In the absence of this information, the statistical significance of features that are actually noise is overestimated and detections claimed that are actually spurious. For this reason, we present an alternative method of computing the probability of false detection that is based on the probability density function (PDF) of the peaks of a random field. It is able to provide a correct estimate of the probability of false detection for the one-, two- and three-dimensional case. We apply this technique to a real two-dimensional interferometric map obtained with ALMA.

In recent years there has been a paradigm shift from centralised to geographically distributed resources. Individual entities are no longer able to host or afford the necessary expertise in-house, and, as a consequence, society increasingly relies on widespread collaborations. Although such collaborations are now the norm for scientific projects, more technical structures providing support to a distributed scientific community without direct financial or other material benefits are scarce. The network of European ALMA Regional Centre (ARC) nodes is an example of such an internationally distributed user support network. It is an organised effort to provide the European ALMA user community with uniform expert support to enable optimal usage and scientific output of the ALMA facility. The network model for the European ARC nodes is described in terms of its organisation, communication strategies and user support.

We investigate the star formation properties of ~800 sources detected in one of the deepest radio surveys at 1.4 GHz. Our sample spans a wide redshift range (~0.1 - 4) and about four orders of magnitude in star formation rate (SFR). It includes both star forming galaxies (SFGs) and active galactic nuclei (AGNs), further divided into radio-quiet and radio-loud objects. We compare the SFR derived from the far infrared luminosity, as traced by Herschel, with the SFR computed from their radio emission. We find that the radio power is a good SFR tracer not only for pure SFGs but also in the host galaxies of RQ AGNs, with no significant deviation with redshift or specific SFR. Moreover, we quantify the contribution of the starburst activity in the SFGs population and the occurrence of AGNs in sources with different level of star formation. Finally we discuss the possibility of using deep radio survey as a tool to study the cosmic star formation history.

We present a sub-100 pc-scale analysis of the CO molecular gas emission and kinematics of the gravitational lens system SDP.81 at redshift 3.042 using Atacama Large Millimetre/submillimetre Array (ALMA) science verification data and a visibility-plane lens reconstruction technique. We find clear evidence for an excitation dependent structure in the unlensed molecular gas distribution, with emission in CO (5-4) being significantly more diffuse and structured than in CO (8-7). The intrinsic line luminosity ratio is r_8-7/5-4 = 0.30 +/- 0.04, which is consistent with other low-excitation starbursts at z ~ 3. An analysis of the velocity fields shows evidence for a star-forming disk with multiple velocity components that is consistent with a merger/post-coalescence merger scenario, and a dynamical mass of M(< 1.56 kpc) = 1.6 +/- 0.6 x 10^10 M_sol . Source reconstructions from ALMA and the Hubble Space Telescope show that the stellar component is offset from the molecular gas and dust components. Together with Karl G. Jansky Very Large Array CO (1-0) data, they provide corroborative evidence for a complex ~2 kpc-scale starburst that is embedded within a larger ~15 kpc structure.

We present a sub-50 pc-scale analysis of the gravitational lens system SDP.81 at redshift 3.042 using Atacama Large Millimetre/submillimetre Array (ALMA) science verification data. We model both the mass distribution of the gravitational lensing galaxy and the pixelated surface brightness distribution of the background source using a novel Bayesian technique that fits the data directly in visibility space. We find the 1 and 1.3 mm dust emission to be magnified by a factor of u_tot = 17.6+/-0.4, giving an intrinsic total star-formation rate of 315+/-60 M_sol/yr and a dust mass of 6.4+/-1.5*10^8 M_sol. The reconstructed dust emission is found to be non-uniform, but composed of multiple regions that are heated by both diffuse and strongly clumped star-formation. The highest surface brightness region is a ~1.9*0.7 kpc disk-like structure, whose small extent is consistent with a potential size-bias in gravitationally lensed starbursts. Although surrounded by extended star formation, with a density of 20-30+/-10 M_sol/yr/kpc^2, the disk contains three compact regions with densities that peak between 120-190+/-20 M_sol/yr/kpc^2. Such star-formation rate densities are below what is expected for Eddington-limited star-formation by a radiation pressure supported starburst. There is also a tentative variation in the spectral slope of the different star-forming regions, which is likely due to a change in the dust temperature and/or opacity across the source.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations from the 2014 Long Baseline Campaign in dust continuum and spectral line emission from the HL Tau region. The continuum images at wavelengths of 2.9, 1.3, and 0.87 mm have unprecedented angular resolutions of 0.075 arcseconds (10 AU) to 0.025 arcseconds (3.5 AU), revealing an astonishing level of detail in the circumstellar disk surrounding the young solar analogue HL Tau, with a pattern of bright and dark rings observed at all wavelengths. By fitting ellipses to the most distinct rings, we measure precise values for the disk inclination (46.72pm0.05 degrees) and position angle (+138.02pm0.07 degrees). We obtain a high-fidelity image of the 1.0 mm spectral index ($\alpha$), which ranges from $\alpha\sim2.0$ in the optically-thick central peak and two brightest rings, increasing to 2.3-3.0 in the dark rings. The dark rings are not devoid of emission, we estimate a grain emissivity index of 0.8 for the innermost dark ring and lower for subsequent dark rings, consistent with some degree of grain growth and evolution. Additional clues that the rings arise from planet formation include an increase in their central offsets with radius and the presence of numerous orbital resonances. At a resolution of 35 AU, we resolve the molecular component of the disk in HCO+ (1-0) which exhibits a pattern over LSR velocities from 2-12 km/s consistent with Keplerian motion around a ~1.3 solar mass star, although complicated by absorption at low blue-shifted velocities. We also serendipitously detect and resolve the nearby protostars XZ Tau (A/B) and LkHa358 at 2.9 mm.

We present Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm continuum images of the asteroid 3 Juno obtained with an angular resolution of 0.042 arcseconds (60 km at 1.97 AU). The data were obtained over a single 4.4 hr interval, which covers 60% of the 7.2 hr rotation period, approximately centered on local transit. A sequence of ten consecutive images reveals continuous changes in the asteroid's profile and apparent shape, in good agreement with the sky projection of the three-dimensional model of the Database of Asteroid Models from Inversion Techniques. We measure a geometric mean diameter of 259pm4 km, in good agreement with past estimates from a variety of techniques and wavelengths. Due to the viewing angle and inclination of the rotational pole, the southern hemisphere dominates all of the images. The median peak brightness temperature is 215pm13 K, while the median over the whole surface is 197pm15 K. With the unprecedented resolution of ALMA, we find that the brightness temperature varies across the surface with higher values correlated to the subsolar point and afternoon areas, and lower values beyond the evening terminator. The dominance of the subsolar point is accentuated in the final four images, suggesting a reduction in the thermal inertia of the regolith at the corresponding longitudes, which are possibly correlated to the location of the putative large impact crater. These results demonstrate ALMA's potential to resolve thermal emission from the surface of main belt asteroids, and to measure accurately their position, geometric shape, rotational period, and soil characteristics.

We present initial results of very high resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of the $z$=3.042 gravitationally lensed galaxy HATLAS J090311.6+003906 (SDP.81). These observations were carried out using a very extended configuration as part of Science Verification for the 2014 ALMA Long Baseline Campaign, with baselines of up to 15 km. We present continuum imaging at 151, 236 and 290 GHz, at unprecedented angular resolutions as fine as 23 milliarcseconds (mas), corresponding to an un-magnified spatial scale of ~180 pc at z=3.042. The ALMA images clearly show two main gravitational arc components of an Einstein ring, with emission tracing a radius of ~1.5". We also present imaging of CO(10-9), CO(8-7), CO(5-4) and H2O line emission. The CO emission, at an angular resolution of ~170 mas, is found to broadly trace the gravitational arc structures but with differing morphologies between the CO transitions and compared to the dust continuum. Our detection of H2O line emission, using only the shortest baselines, provides the most resolved detection to date of thermal H2O emission in an extragalactic source. The ALMA continuum and spectral line fluxes are consistent with previous Plateau de Bure Interferometer and Submillimeter Array observations despite the impressive increase in angular resolution. Finally, we detect weak unresolved continuum emission from a position that is spatially coincident with the center of the lens, with a spectral index that is consistent with emission from the core of the foreground lensing galaxy.

We discuss the potential of a next generation space-borne CMB experiment for studies of extragalactic sources with reference to COrE+, a project submitted to ESA in response to the M4 call. We consider three possible options for the telescope size: 1m, 1.5m and 2m (although the last option is probably impractical, given the M4 boundary conditions). The proposed instrument will be far more sensitive than Planck and will have a diffraction-limited angular resolution. These properties imply that even the 1m telescope option will perform substantially better than Planck for studies of extragalactic sources. The source detection limits as a function of frequency have been estimated by means of realistic simulations. The most significant improvements over Planck results are presented for each option. COrE+ will provide much larger samples of truly local star-forming galaxies, making possible analyses of the properties of galaxies (luminosity functions, dust mass functions, star formation rate functions, dust temperature distributions, etc.) across the Hubble sequence. Even more interestingly, COrE+ will detect, at |b|> 30 deg, thousands of strongly gravitationally lensed galaxies. Such large samples are of extraordinary astrophysical and cosmological value in many fields. Moreover, COrE+ high frequency maps will be optimally suited to pick up proto-clusters of dusty galaxies, i.e. to investigate the evolution of large scale structure at larger redshifts than can be reached by other means. Thanks to its high sensitivity COrE+ will also yield a spectacular advance in the blind detection of extragalactic sources in polarization. This will open a new window for studies of radio source polarization and of the global properties of magnetic fields in star forming galaxies and of their relationships with SFRs.

We compare various star formation rate (SFR) indicators for star-forming galaxies at $1.4<z<2.5$ in the COSMOS field. The main focus is on the SFRs from the far-IR (PACS-Herschel data) with those from the ultraviolet, for galaxies selected according to the BzK criterion. FIR-selected samples lead to a vastly different slope of the SFR-stellar mass ($M_*$) relation, compared to that of the dominant main sequence population as measured from the UV, since the FIR selection picks predominantly only a minority of outliers. However, there is overall agreement between the main sequences derived with the two SFR indicators, when stacking on the PACS maps the BzK-selected galaxies. The resulting logarithmic slope of the SFR-{$M_*$} relation is $\sim0.8-0.9$, in agreement with that derived from the dust-corrected UV-luminosity. Exploiting deeper 24$\mu$m-Spitzer data we have characterized a sub-sample of galaxies with reddening and SFRs poorly constrained, as they are very faint in the $B$ band. The combination of Herschel with Spitzer data have allowed us to largely break the age/reddening degeneracy for these intriguing sources, by distinguishing whether a galaxy is very red in B-z because of being heavily dust reddened, or whether because star formation has been (or is being) quenched. Finally, we have compared our SFR(UV) to the SFRs derived by stacking the radio data and to those derived from the H$\alpha$ luminosity of a sample of star-forming galaxies at $1.4<z<1.7$. The two sets of SFRs are broadly consistent as they are with the SFRs derived from the UV and by stacking the corresponding PACS data in various mass bins.

We study the relationship between the K-band and the sub-millimetre (submm) emissions of nearby galaxies by computing the bivariate K-band-submm luminosity function (BLF) of the Herschel Reference Survey (HRS), a volume-limited sample observed in submm with Hersche/SPIRE. We derive the BLF from the K-band and submm cumulative distributions using a copula method. Using the BLF allows us to derive the relationship between the luminosities on more solid statistical ground. The analysis shows that over the whole HRS sample, no statistically meaningful conclusion can be derived for any relationship between the K-band and the submm luminosity. However, a very tight relationship between these luminosities is highlighted, by restricting our analysis to late-type galaxies. The luminosity function of late-type galaxies computed in the K-band and in the submm are dependent and the dependence is caused by the link, between the stellar mass and the cold dust mass, which has been already observed

The goal of this work is to compute the number density of far-IR selected galaxies in the comoving frame and along the past lightcone of observationally constrained LTB models and to compare those results with their standard model counterparts. We derived integral number densities and differential number densities using different cosmological distance definitions in the LTB models. Then, we computed selection functions and consistency functions for the luminosity functions in the combined fields of the Herschel/PACS evolutionary probe (PEP) survey in both standard and void cosmologies, from which we derived the observed values of the above-mentioned densities. We used the Kolmogorov-Smirnov statistics to study both the evolution of the consistency functions and its connection to the evolution of the comoving density of sources. Finally, we fitted the power-law behaviour of the densities along the observer's past lightcone. The analysis of the comoving number density shows that the increased flexibility of the LTB models is not enough to fit the observed redshift evolution of the number counts, if it is specialised to a recent best-fit giant void parametrisation. The results for the power-law fits of the densities along the observer's past lightcone show general agreement across both cosmological models studied here around a slope of -2.5 $\pm$ 0.1 for the integral number density on the luminosity-distance volumes. The differential number densities show much bigger slope discrepancies. We conclude that the differential number densities on the observer's past lightcone were still rendered dependent on the cosmological model by the flux limits of the PEP survey. In addition, we show that an intrinsic evolution of the sources must be assumed to fit the comoving number-density redshift evolution in the giant void parametrisation for the LTB models used in this work.

We present Karl G. Jansky Very Large Array (VLA) observations of 44 GHz continuum and CO J=2-1 line emission in BR1202-0725 at z=4.7 (a starburst galaxy and quasar pair) and BRI1335-0417 at z=4.4 (also hosting a quasar). With the full 8 GHz bandwidth capabilities of the upgraded VLA, we study the (rest-frame) 250 GHz thermal dust continuum emission for the first time along with the cold molecular gas traced by the Low-J CO line emission. The measured CO J=2-1 line luminosities of BR1202-0725 are L'(CO) = (8.7+/-0.8)x10^10 K km/s pc^2 and L'(CO) = (6.0+/-0.5)x10^10 K km/s pc^2 for the submm galaxy (SMG) and quasar, which are equal to previous measurements of the CO J=5-4 line luminosities implying thermalized line emission and we estimate a combined cold molecular gas mass of ~9x10^10 Msun. In BRI1335-0417 we measure L'(CO) = (7.3+/-0.6)x10^10 K km/s pc^2. We detect continuum emission in the SMG BR1202-0725 North (S(44GHz) = 51+/-6 microJy), while the quasar is detected with S(44GHz) = 24+/-6 microJy and in BRI1335-0417 we measure S(44GHz) = 40+/-7 microJy. Combining our continuum observations with previous data at (rest-frame) far-infrared and cm-wavelengths, we fit three component models in order to estimate the star-formation rates. This spectral energy distribution fitting suggests that the dominant contribution to the observed 44~GHz continuum is thermal dust emission, while either thermal free-free or synchrotron emission contributes less than 30%.

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

[Abridged] We study the evolution of the dust temperatures of galaxies in the SFR-M* plane up to z~2 using observations from the Herschel Space Observatory. Starting from a sample of galaxies with reliable star-formation rates (SFRs), stellar masses (M*) and redshift estimates, we grid the SFR-M* parameter space in several redshift ranges and estimate the mean Tdust of each SFR-M*-z bin. Dust temperatures are inferred using the stacked far-infrared flux densities of our SFR-M*-z bins. At all redshifts, Tdust increases with infrared luminosities (LIR), specific SFRs (SSFR; i.e., SFR/M*) and distances with respect to the main sequence (MS) of the SFR-M* plane (i.e., D_SSFR_MS=log[SSFR(galaxy)/SSFR_MS(M*,z)]). The Tdust-SSFR and Tdust-D_SSFR_MS correlations are statistically more significant than the Tdust-LIR one. While the slopes of these three correlations are redshift-independent, their normalizations evolve from z=0 and z~2. We convert these results into a recipe to derive Tdust from SFR, M* and z. The existence of a strong Tdust-D_SSFR_MS correlation provides us with information on the dust and gas content of galaxies. (i) The slope of the Tdust-D__SSFR_MS correlation can be explained by the increase of the star-formation efficiency (SFE; SFR/Mgas) with D_SSFR_MS as found locally by molecular gas studies. (ii) At fixed D_SSFR_MS, the constant Tdust observed in galaxies probing large ranges in SFR and M* can be explained by an increase or decrease of the number of star-forming regions with comparable SFE enclosed in them. (iii) At high redshift, the normalization towards hotter temperature of the Tdust-D_SSFR_MS correlation can be explained by the decrease of the metallicities of galaxies or by the increase of the SFE of MS galaxies. All these results support the hypothesis that the conditions prevailing in the star-forming regions of MS and far-above-MS galaxies are different.

This is the first paper of a series aiming at investigating galaxy formation and evolution in the giant-void class of the Lemaitre-Tolman-Bondi (LTB) models that best fits current cosmological observations. Here we investigate the Luminosity Function (LF) methodology, and how its estimates would be affected by a change on the cosmological model assumed in its computation. Are the current observational constraints on the allowed Cosmology enough to yield robust LF results? We use the far-infrared source catalogues built on the observations performed with the Herschel/PACS instrument, and selected as part of the PACS evolutionary probe (PEP) survey. Schechter profiles are obtained in redshift bins up to z approximately 4, assuming comoving volumes in both the standard model, that is, Friedmann-Lemaitre-Robertson-Walker metric with a perfect fluid energy-momentum tensor, and non-homogeneous LTB dust models, parametrized to fit the current combination of results stemming from the observations of supernovae Ia, the cosmic microwave background, and baryonic acoustic oscillations. We find that the luminosity functions computed assuming both the standard model and LTB void models show in general good agreement. However, the faint-end slope in the void models shows a significant departure from the standard model up to redshift 0.4. We demonstrate that this result is not artificially caused by the used LF estimator which turns out to be robust under the differences in matter-energy density profiles of the models. The differences found in the LF slopes at the faint end are due to variation in the luminosities of the sources, which depend on the geometrical part of the model. It follows that either the standard model is over-estimating the number density of faint sources or the void models are under-estimating it.

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

We report the PACS-100um/160um detections of a sample of 42 GALEX-selected and FIR-detected Lyman break galaxies (LBGs) at z ~ 1 located in the COSMOS field and analyze their ultra-violet (UV) to far-infrared (FIR) properties. The detection of these LBGs in the FIR indicates that they have a dust content high enough so that its emission can be directly detected. According to a spectral energy distribution (SED) fitting with stellar population templates to their UV-to-near-IR observed photometry, PACS-detected LBGs tend to be bigger, more massive, dustier, redder in the UV continuum, and UV-brighter than PACS-undetected LBGs. PACS-detected LBGs at z ~ 1 are mostly disk-like galaxies and are located over the green-valley and red sequence of the color-magnitude diagram of galaxies at their redshift. By using their UV and IR emission, we find that PACS-detected LBGs tend to be less dusty and have slightly higher total star-formation rates (SFRs) than other PACS-detected UV-selected galaxies within their same redshift range. As a consequence of the selection effect due to the depth of the FIR observations employed, all our PACS-detected LBGs are LIRGs. However, none of them are in the ULIRG regime, where the FIR observations are complete. The finding of ULIRGs-LBGs at higher redshifts suggests an evolution of the FIR emission of LBGs with cosmic time. In an IRX-$\beta$ diagram, PACS-detected LBGs at z ~ 1 tend to be located around the relation for local starburst similarly to other UV-selected PACS-detected galaxies at their same redshift. Consequently, the dust-correction factors obtained with their UV continuum slope allow to determine their total SFR, unlike at higher redshifts. However, the dust attenuation derived from UV to NIR SED fitting overestimates the total SFR for most of our PACS-detected LBGs in age-dependent way: the overestimation factor is higher in younger galaxies.

This paper deals with the detection problem of extragalactic point-sources in multi-frequency, microwave sky maps that will be obtainable in future cosmic microwave background radiation (CMB) experiments with instruments capable of very high spatial resolution. With spatial resolutions that can be of order of 0.1-1.0 arcsec or better, the extragalactic point-sources will appear isolated. The same holds also for the compact structures due to the Sunyaev-Zeldovich (SZ) effect (both thermal and kinetic). This situation is different from the maps obtainable with instruments as WMAP or PLANCK where, because of the smaller spatial resolution (approximately 5-30 arcmin), the point-sources and the compact structures due to the SZ effect form a uniform noisy background (the "confusion noise"). Hence, the point-source detection techniques developed in the past are based on the assumption that all the emissions that contribute to the microwave background can be modeled with homogeneous and isotropic (often Gaussian) random fields and make use of the corresponding spatial power-spectra. In the case of very high resolution observations such an assumption cannot be adopted since it still holds only for the CMB. Here, we propose an approach based on the assumption that the diffuse emissions that contribute to the microwave background can be locally approximated by two-dimensional low order polynomials. In particular, two sets of numerical techniques are presented containing two different algorithms each. The performance of the algorithms is tested with numerical experiments that mimic the physical scenario expected for high Galactic latitude observations with the Atacama Large Millimeter/Submillimeter Array (ALMA).

By making use of Herschel-PEP observations of the COSMOS and Extended Groth Strip fields, we have estimated the dependence of the clustering properties of FIR-selected sources on their 100um fluxes. Our analysis shows a tendency for the clustering strength to decrease with limiting fluxes: r0(S100um >8 mJy)~4.3 Mpc and r0(S100um >5 mJy)~5.8 Mpc. These values convert into minimum halo masses Mmin~10^{11.6} Msun for sources brighter than 8 mJy and Mmin~10^{12.4} Msun for S100um > 5 mJy galaxies. We show such an increase of the clustering strength to be due to an intervening population of z~2 sources, which are very strongly clustered and whose relative contribution, equal to about 10% of the total counts at S100um > 2 mJy, rapidly decreases for brighter flux cuts. By removing such a contribution, we find that z <~ 1 FIR galaxies have approximately the same clustering properties, irrespective of their flux level. The above results were then used to investigate the intrinsic dependence on cosmic epoch of the clustering strength of dusty star-forming galaxies between z~0 and z~2.5. In order to remove any bias in the selection process, the adopted sample only includes galaxies observed at the same rest-frame wavelength, lambda~60 um, which have comparable luminosities and therefore star-formation rates (SFR>~100 Msun/yr). Our analysis shows that the same amount of (intense) star forming activity takes place in extremely different environments at the different cosmological epochs. For z<~1 the hosts of such star forming systems are small, Mmin~10^{11} Msun, isolated galaxies. High (z~2) redshift star formation instead seems to uniquely take place in extremely massive/cluster-like halos, Mmin~10^{13.5} Msun, which are associated with the highest peaks of the density fluctuation field at those epochs. (abridged)

We present results from the deepest Herschel-PACS (Photodetector Array Camera and Spectrometer) far-infrared blank field extragalactic survey, obtained by combining observations of the GOODS (Great Observatories Origins Deep Survey) fields from the PACS Evolutionary Probe (PEP) and GOODS-Herschel key programmes. We describe data reduction and the construction of images and catalogues. In the deepest parts of the GOODS-S field, the catalogues reach 3-sigma depths of 0.9, 0.6 and 1.3 mJy at 70, 100 and 160 um, respectively, and resolve ~75% of the cosmic infrared background at 100um and 160um into individually detected sources. We use these data to estimate the PACS confusion noise, to derive the PACS number counts down to unprecedented depths and to determine the infrared luminosity function of galaxies down to LIR=10^11 Lsun at z~1 and LIR=10^12 Lsun at z~2, respectively. For the infrared luminosity function of galaxies, our deep Herschel far-infrared observations are fundamental because they provide more accurate infrared luminosity estimates than those previously obtained from mid-infrared observations. Maps and source catalogues (>3-sigma) are now publicly released. Combined with the large wealth of multi-wavelength data available for the GOODS fields, these data provide a powerful new tool for studying galaxy evolution over a broad range of redshifts.

(Abridged) We study relationships between the SFR and the nuclear properties of X-ray selected AGNs out to z=2.5, using far-IR data in three extragalactic deep fields as part of the PACS Evolutionary Probe (PEP) program. Guided by studies of intrinsic infra-red AGN SEDs, we show that the majority of the FIR emission in AGNs is produced by cold dust heated by star-formation. We uncover characteristic redshift-dependent trends between the mean FIR luminosity (L_fir) and accretion luminosity (L_agn) of AGNs. At low AGN luminosities, accretion and SFR are uncorrelated at all redshifts, consistent with a scenario where most low-luminosity AGNs are primarily fueled by secular processes in their host galaxies. At high AGN luminosities, a significant correlation is observed between L_fir and L_agn, but only among AGNs at low and moderate redshifts (z<1). We interpret this as a signature of the increasing importance of major-mergers in driving both the growth of super-massive black holes (SMBHs) and global star-formation in their hosts at high AGN luminosities. However, we also find that the enhancement of SFR in luminous AGNs weakens or disappears at high redshifts (z>1). This suggests that the role of mergers in SMBH-galaxy co-evolution is less important at these epochs. At all redshifts, we find essentially no relationship between L_fir and nuclear obscuration across five orders of magnitude in obscuring column density, suggesting that various different mechanisms are likely to be responsible for obscuring X-rays in active galaxies. We explain our results within a scenario in which two different modes of SMBH fueling operate among low- and high-luminosity AGNs. We postulate, guided by emerging knowledge about the properties of high redshift galaxies, that the dominant mode of accretion among high-luminosity AGNs evolves with redshift.

Detection of a signal hidden by noise within a time series is an important problem in many astronomical searches, i.e. for light curves containing the contributions of periodic/semi-periodic components due to rotating objects and all other astrophysical time-dependent phenomena. One of the most popular tools for use in such studies is the "periodogram", whose use in an astronomical context is often not trivial. The "optimal" statistical properties of the periodogram are lost in the case of irregular sampling of signals, which is a common situation in astronomical experiments. Parts of these properties are recovered by the "Lomb-Scargle" (LS) technique, but at the price of theoretical difficulties, that can make its use unclear, and of algorithms that require the development of dedicated software if a fast implementation is necessary. Such problems would be irrelevant if the LS periodogram could be used to significantly improve the results obtained by approximated but simpler techniques. In this work we show that in many astronomical applications simpler techniques provide results similar to those obtainable with the LS periodogram. The meaning of the "Nyquist frequency" is also discussed in the case of irregular sampling.

We propose two detection techniques that take advantage of a small sky area approximation and are based on modifications of the "internal linear combination" (ILC) method, an approach widely used in Cosmology for the separation of the various components that contribute to the microwave background. The main advantage of the proposed approach, especially in handling multi-frequency maps of the same region, is that it does not require the "a priori" knowledge of the spatial power-spectrum of either the CMB and/or the Galactic foreground. Hence, it is more robust, easier and more intuitive to use. The performance of the proposed algorithms is tested with numerical experiments that mimic the physical scenario expected for high Galactic latitude observations with the Atacama Large Millimeter/submillimeter Array (ALMA).

We use Herschel-PACS far-infrared data, combined with previous multi-band information and mid-IR spectra, to properly account for the presence of an active nucleus and constrain its energetic contribution in luminous infrared (IR) sources at z\sim2. The sample is composed of 24 sources in the GOODS-South field, with typical IR luminosity of 10^{12} Lo. Data from the 4 Ms Chandra X-ray imaging in this field are also used to identify and characterize AGN emission. We reproduce the observed spectral energy distribution (SED), decomposed into a host-galaxy and an AGN component. A smooth-torus model for circum-nuclear dust is used to account for the direct and re-processed contribution from the AGN. We confirm that galaxies with typical L_{8-1000um}\sim10^{12}Lo at z\sim2 are powered predominantly by star-formation. An AGN component is present in nine objects (\sim35% of the sample) at the 3sigma confidence level, but its contribution to the 8-1000 um emission accounts for only \sim5% of the energy budget. The AGN contribution rises to \sim23% over the 5-30 um range (in agreement with Spitzer IRS results) and to \sim60% over the narrow 2-6 um range. The presence of an AGN is confirmed by X-ray data for 3 (out of nine) sources, with X-ray spectral analysis indicating the presence of significant absorption, i.e. NH\sim10^{23}-10^{24} cm^{-2}. An additional source shows indications of obscured AGN emission from X-ray data. The comparison between the mid-IR--derived X-ray luminosities and those obtained from X-ray data suggests that obscuration is likely present also in the remaining six sources that harbour an AGN according to the SED-fitting analysis.

The analysis of the physical properties of low-redshift Ly$\alpha$ emitters (LAEs) can provide clues in the study of their high-redshift analogues. At $z \sim 0.3$, LAEs are bright enough to be detected over almost the entire electromagnetic spectrum and it is possible to carry out a more precise and complete study than at higher redshifts. In this study, we examine the UV and IR emission, dust attenuation, SFR and morphology of a sample of 23 GALEX-discovered star-forming (SF) LAEs at $z \sim 0.3$ with direct UV (GALEX), optical (ACS) and FIR (PACS and MIPS) data. Using the same UV and IR limiting luminosities, we find that LAEs at $z\sim 0.3$ tend to be less dusty, have slightly higher total SFRs, have bluer UV continuum slopes, and are much smaller than other galaxies that do not exhibit Ly$\alpha$ emission in their spectrum (non-LAEs). These results suggest that at $z \sim 0.3$ Ly$\alpha$ photons tend to escape from small galaxies with low dust attenuation. Regarding their morphology, LAEs belong to Irr/merger classes, unlike non-LAEs. Size and morphology represent the most noticeable difference between LAEs and non-LAEs at $z \sim 0.3$. Furthermore, the comparison of our results with those obtained at higher redshifts indicates that either the Ly$\alpha$ technique picks up different kind of galaxies at different redshifts or that the physical properties of LAEs are evolving with redshift.

We compare the average star formation (SF) activity in X-ray selected AGN hosts with mass-matched control inactive galaxies,including star forming and quiescent sources, at 0.5<z<2.5. Recent observations carried out by PACS, the 60-210um Herschel photometric camera, in GOODS-S, GOODS-N and COSMOS allow us to unbiasedly estimate the far-IR luminosity, and hence the SF properties, of the two samples. Accurate AGN host stellar masses are measured by decomposing their total emission into the stellar and nuclear components. We find a higher average SF activity in AGN hosts with respect to non-AGNs. The level of SF enhancement is modest (~0.26dex at ~3sigma) at low X-ray luminosities (Lx<~10^43.5erg/s) and more pronounced (0.56dex at >10sigma) for bright AGNs. However, when comparing to star forming galaxies only, AGN hosts are broadly consistent with the locus of their `main sequence'. We investigate the relative far-IR luminosity distributions of active and inactive galaxies, and find a higher fraction of PACS detected, hence normal and highly star forming systems among AGN hosts. Although different interpretations are possible, we explain our findings as a consequence of a twofold AGN growth path: faint AGNs evolve through secular processes, with instantaneous AGN accretion not tightly linked to the current total SF in the host, while luminous AGNs co-evolve with their hosts through periods of enhanced AGN activity and SF, possibly through major mergers. While an increased SF with respect to non-AGNs of similar mass is expected in the latter, we interpret the modest SF offsets measured in low-Lx AGN hosts as either a) generated by non-synchronous accretion and SF histories in a merger scenario or b) due to possible connections between instantaneous SF and accretion that can be induced by smaller scale (non-major merger) mechanisms. Far-IR luminosity distributions favour the latter scenario.

In this work we analyze the physical properties of a sample of 56 spectroscopically selected star-forming (SF) Ly$\alpha$ emitting galaxies at 2.0$\lesssim$z$\lesssim$3.5 using both a spectral energy distribution (SED) fitting procedure from rest-frame UV to mid-IR and direct 160$\mu$m observations taken with the Photodetector Array Camera & Spectrometer (PACS) instrument onboard \emph{Herschel Space Observatory}. We define LAEs as those Ly$\alpha$ emitting galaxies whose rest-frame Ly$\alpha$ equivalent widths (Ly$\alpha$ EW$_{rest-frame}$) are above 20\AA, the typical threshold in narrow-band searches. Ly$\alpha$ emitting galaxies with Ly$\alpha$ EW$_{rest-frame}$ are called non-LAEs. As a result of an individual SED fitting for each object, we find that the studied sample of LAEs contains galaxies with ages mostly below 100Myr and a wide variety of dust attenuations, SFRs, and stellar masses. The heterogeneity in the physical properties is also seen in the morphology, ranging from bulge-like galaxies to highly clumpy systems. In this way, we find that LAEs at 2.0$\lesssim$z$\lesssim$3.5 are very diverse, and do not have a bimodal nature, as suggested in previous works. Furthermore, the main difference between LAEs and non-LAEs is their dust attenuation, because LAEs are not as dusty as non-LAEs. On the FIR side, four galaxies of the sample (two LAEs and two non-LAEs) have PACS-FIR counterparts. Their total IR luminosity place all of them in the ULIRG regime and are all dusty objects, with A$_{1200}$$\gtrsim$4mag. This is an indication from direct FIR measurements that dust and Ly$\alpha$ emission are not mutually exclusive. This population of red and dusty LAEs is not seen at z$\sim$0.3, suggesting an evolution with redshift of the IR nature of galaxies selected via their Ly$\alpha$ emission.

We study a sample of 61 submillimetre galaxies (SMGs) selected from ground-based surveys, with known spectroscopic redshifts and observed with Herschel as part of the PACS Evolutionary Probe (PEP) and the Herschel Multi-tiered Extragalactic Survey (HerMES) key programmes. We use the broad far-infrared wavelength coverage (100-600um) provided by the combination of PACS and SPIRE observations. Using a power-law temperature distribution model to derive infrared luminosities and dust temperatures, we measure a dust emissivity spectral index for SMGs of beta=2.0+/-0.2. Our results unveil the diversity of the SMG population. Some SMGs exhibit extreme infrared luminosities of ~10^13 Lsun and relatively warm dust components, while others are fainter (~10^12 Lsun) and are biased towards cold dust temperatures. The extreme infrared luminosities of some SMGs (LIR>10^12.7 Lsun, 26/61 systems) imply SFRs of >500Msun yr^-1. Such high SFRs are difficult to reconcile with a secular mode of star formation, and may instead correspond to a merger-driven stage in the evolution of these galaxies. Another observational argument in favour of this scenario is the presence of dust temperatures warmer than that of SMGs of lower luminosities (~40K as opposed to ~25K), consistent with observations of local ULIRGs triggered by major mergers and with results from hydrodynamic simulations of major mergers combined with radiative transfer calculations. Luminous SMGs are also offset from normal star-forming galaxies in the stellar mass-SFR plane, suggesting that they are undergoing starburst events with short duty cycles, compatible with the major merger scenario. On the other hand, a significant fraction of the low infrared luminosity SMGs have cold dust temperatures, are located close to the main sequence of star formation, and thus might be evolving through a secular mode of star formation. [abridged]

We combine Herschel-PACS data from the PEP program with Spitzer 24 um and 16 um photometry and ultra deep IRS mid-infrared spectra, to measure the mid- to far-infrared spectral energy distribution (SED) of 0.7<z<2.5 normal star forming galaxies around the main sequence (the redshift-dependent relation of star formation rate and stellar mass). Our deep data confirm from individual far-infrared detections that z~2 star formation rates are overestimated if based on 24 um fluxes and SED templates that are calibrated via local trends with luminosity. Galaxies with similar ratios of rest-frame nuLnu(8) to 8-1000 um infrared luminosity (LIR) tend to lie along lines of constant offset from the main sequence. We explore the relation between SED shape and offset in specific star formation rate (SSFR) from the redshift-dependent main sequence. Main sequence galaxies tend to have a similar nuLnu(8)/LIR regardless of LIR and redshift, up to z~2.5, and nuLnu(8)/LIR decreases with increasing offset above the main sequence in a consistent way at the studied redshifts. We provide a redshift-independent calibration of SED templates in the range of 8--60 um as a function of log(SSFR) offset from the main sequence. Redshift dependency enters only through the evolution of the main sequence with time. Ultra deep IRS spectra match these SED trends well and verify that they are mostly due to a change in ratio of PAH to LIR rather than continua of hidden AGN. Alternatively, we discuss the dependence of nuLnu(8)/LIR on LIR. Same nuLnu(8)/LIR is reached at increasingly higher LIR at higher redshift, with shifts relative to local by 0.5 and 0.8 dex in log(LIR) at redshifts z~1 and z~2. Corresponding SED template calibrations are provided for use if no stellar masses are in hand. For most of those z~2 star forming galaxies that also host an AGN, the mid-infrared is dominated by the star forming component.

Star formation in massive galaxies is quenched at some point during hierarchical mass assembly. To understand where and when the quenching processes takes place, we study the evolution of the total star formation rate per unit total halo mass (\Sigma(SFR/M)) in three different mass scales: low mass halos (field galaxies), groups, and clusters, up to a redshift ~1.6. We use deep far-infrared PACS data at 100 and 160 um to accurately estimate the total star formation rate of the Luminous Infrared Galaxy population of 9 clusters with mass ~10^{15} M_{\odot}, and 9 groups/poor clusters with mass ~ 5 x 10^{13} M_{\odot}. Estimates of the field \Sigma(SFR/M) are derived from the literature, by dividing the star formation rate density by the mean comoving matter density of the universe. The field \Sigma(SFR/M) increases with redshift up to z~1 and it is constant thereafter. The evolution of the \Sigma(SFR/M)-z relation in galaxy systems is much faster than in the field. Up to redshift z~0.2, the field has a higher \Sigma(SFR/M) than galaxy groups and galaxy clusters. At higher redshifts, galaxy groups and the field have similar \Sigma(SFR/M), while massive clusters have significantly lower \Sigma(SFR/M) than both groups and the field. There is a hint of a reversal of the SFR activity vs. environment at z~1.6, where the group \Sigma(SFR/M) lies above the field \Sigma(SFR/M)-z relation. We discuss possible interpretations of our results in terms of the processes of downsizing, and star-formation quenching.

One remaining open question regarding the physical properties of Ly$\alpha$ emitters (LAEs) is their dust content and its evolution with redshift. The variety of results is large and with those reported by now is difficult to establish clear relations between dust, other fundamental parameters of galaxies (star-formation rate, metallicity or age) and redshift. In this Letter, we report \emph{Herschel} PACS-100$\mu$m, PACS-160$\mu$m and \emph{Spitzer} MIPS-24$\mu$m detections of a sample of spectroscopically GALEX-selected LAEs at z$\sim$0.3 and $\sim$1.0. Five out of ten and one out of two LAEs are detected in, at least, one PACS band at z$\sim$0.3 and $\sim$1.0, respectively. These measurements have a great importance given that they allow us to quantify, for the first time, the dust content of LAEs from direct FIR observations. MIPS-24$\mu$m detections allow us to determine IR properties of the PACS-undetected LAEs. We obtain that mid-IR/FIR detected star-forming (SF) LAEs at z$\sim$0.3 have dust content within 0.75$\lesssim$ $A_{1200\AA}$ $\lesssim$2.0, with a median value of A$_{1200\textrm{\AA}}$$\sim$1.1. This range broadens out to 0.75$\lesssim$ $A_{1200\AA}$ $\lesssim$2.5 when considering those LAEs at z$\sim$1.0. Only one SF LAE is undetected both in MIPS-24$\mu$m and PACS, with $A_{1200\AA}$ $\lesssim$0.75. These results seem to be larger than those reported for high-redshift LAEs and, therefore, although an evolutionary trend is not clearly seen, it could point out that low-redshift LAEs are dustier than high-redshift ones. However, the diverse methods used could introduce a systematic offset in the results.

Two main modes of star formation are know to control the growth of galaxies: a relatively steady one in disk-like galaxies, defining a tight star formation rate (SFR)-stellar mass sequence, and a starburst mode in outliers to such a sequence which is generally interpreted as driven by merging. Such starburst galaxies are rare but have much higher SFRs, and it is of interest to establish the relative importance of these two modes. PACS/Herschel observations over the whole COSMOS and GOODS-South fields, in conjunction with previous optical/near-IR data, have allowed us to accurately quantify for the first time the relative contribution of the two modes to the global SFR density in the redshift interval 1.5<z<2.5, i.e., at the cosmic peak of the star formation activity. The logarithmic distributions of galaxy SFRs at fixed stellar mass are well described by Gaussians, with starburst galaxies representing only a relatively minor deviation that becomes apparent for SFRs more than 4 times higher than on the main sequence. Such starburst galaxies represent only 2% of mass-selected star forming galaxies and account for only 10% of the cosmic SFR density at z~2. Only when limited to SFR>1000M(sun)/yr, off-sequence sources significantly contribute to the SFR density (46+/-20%). We conclude that merger-driven starbursts play a relatively minor role for the formation of stars in galaxies, whereas they may represent a critical phase towards the quenching of star formation and morphological transformation in galaxies.

We compare multi-wavelength SFR indicators out to z~3 in GOODS-South. Our analysis uniquely combines U-to-8um photometry from FIREWORKS, MIPS 24um and PACS 70, 100, and 160um photometry from the PEP survey, and Ha spectroscopy from the SINS survey. We describe a set of conversions that lead to a continuity across SFR indicators. A luminosity-independent conversion from 24um to total infrared luminosity yields estimates of LIR that are in the median consistent with the LIR derived from PACS photometry, albeit with significant scatter. Dust correction methods perform well at low to intermediate levels of star formation. They fail to recover the total amount of star formation in systems with large SFR_IR/SFR_UV ratios, typically occuring at the highest SFRs (SFR_UV+IR \gtrsim 100 Msun/yr) and redshifts (z \gtrsim 2.5) probed. Finally, we confirm that Ha-based SFRs at 1.5<z<2.6 are consistent with SFR_SED and SFR_UV+IR provided extra attenuation towards HII regions is taken into account (Av,neb = Av,continuum / 0.44). With the cross-calibrated SFR indicators in hand, we perform a consistency check on the star formation histories inferred from SED modeling. We compare the observed SFR-M relations and mass functions at a range of redshifts to equivalents that are computed by evolving lower redshift galaxies backwards in time. We find evidence for underestimated stellar ages when no stringent constraints on formation epoch are applied. We demonstrate how resolved SED modeling, or alternatively deep UV data, may help to overcome this bias. The age bias is most severe for galaxies with young stellar populations, and reduces towards older systems. Finally, our analysis suggests that SFHs typically vary on timescales that are long (at least several 100 Myr) compared to the galaxies' dynamical time.

Deep far-infrared photometric surveys studying galaxy evolution and the nature of the cosmic infrared background are a key strength of the Herschel mission. We describe the scientific motivation for the PACS Evolutionary Probe (PEP) guaranteed time key program and its role in the complement of Herschel surveys, and the field selection which includes popular multiwavelength fields such as GOODS, COSMOS, Lockman Hole, ECDFS, EGS. We provide an account of the observing strategies and data reduction methods used. An overview of first science results illustrates the potential of PEP in providing calorimetric star formation rates for high redshift galaxy populations, thus testing and superseeding previous extrapolations from other wavelengths, and enabling a wide range of galaxy evolution studies.

The cosmic infrared background (CIB) includes roughly half of the energy radiated by all galaxies at all wavelengths across cosmic time, as observed at the present epoch. The PACS Evolutionary Probe (PEP) survey is exploited here to study the CIB and its redshift differential, at 70, 100 and 160 micron, where the background peaks. Combining PACS observations of the GOODS-S, GOODS-N, Lockman Hole and COSMOS areas, we define number counts spanning over more than two orders of magnitude in flux: from ~1 mJy to few hundreds mJy. Stacking of 24 micron sources and P(D) statistics extend the analysis down to ~0.2 mJy. Taking advantage of the wealth of ancillary data in PEP fields, differential number counts and CIB are studied up to z=5. Based on these counts, we discuss the effects of confusion on PACS blank field observations and provide confusion limits for the three bands considered. The total CIB surface brightness emitted above PEP 3 sigma flux limits is 4.52 +/- 1.18, 8.35 +/- 0.95 and 9.49 +/- 0.59 [nW/m2/sr] at 70, 100, and 160 micron, respectively. These values correspond to 58 +/- 7% and 74 +/- 5% of the COBE/DIRBE CIB direct measurements at 100 and 160 micron. Employing the P(D) analysis, these fractions increase to ~65% and ~89%. More than half of the resolved CIB was emitted at redshift z<=1. The 50%-light redshifts lie at z=0.58, 0.67 and 0.73 at the three PACS wavelengths. The distribution moves towards earlier epochs at longer wavelengths: while the 70 micron CIB is mainly produced by z<=1.0 objects, the contribution of z>1.0 sources reaches 50% at 160 micron. Most of the CIB resolved in the three PACS bands was emitted by galaxies with infrared luminosities in the range 1e11-1e12 L(sun).

ABRIDGED-This paper presents the first direct estimate of the 3D clustering properties of far-infrared sources up to z~3. This has been possible thanks to the Pacs Evolutionary Probe (PEP) survey of the GOODS South field performed with the PACS instrument onboard the Herschel Satellite. An analysis of the two-point correlation function over the whole redshift range spanned by the data reports for the correlation length, r_0~6.3 Mpc and r_0~6.7 Mpc, respectively at 100um and 160um, corresponding to dark matter halo masses M>~10^{12.4} M_sun. Objects at z~2 instead seem to be more strongly clustered, with r_0~19 Mpc and r_0~17 Mpc in the two considered PACS channels. This dramatic increase of the correlation length between z~1 and z~2 is connected with the presence of a wide, M>~10^{14} M_sun, filamentary structure which includes more than 50% of the sources detected at z~2. An investigation of the properties of such sources indicates the possibility for boosted star-forming activity in those which reside within the overdense environment with respect of more isolated galaxies found in the same redshift range. Lastly, we also present our results on the evolution of the relationship between luminous and dark matter in star-forming galaxies between z~1 and z~2. We find that the increase of (average) stellar mass in galaxies <M*> between z~1 and z~2 is about a factor 10 lower than that of the dark matter haloes hosting such objects (<M*>[z~1]/<M*>[z~2] ~ 0.4 vs M_{halo}[z~1]/M_{halo}[z~2] ~ 0.04). Our findings agree with the evolutionary picture of downsizing whereby massive galaxies at z~2 were more actively forming stars than their z~1 counterparts, while at the same time contained a lower fraction of their mass in the form of luminous matter.

In experimental microwave maps, point-sources can strongly affect the estimation of the power-spectrum and/or the test of Gaussianity of the Cosmic Microwave Background (CMB) component. As a consequence, their removal from the sky maps represents a critical step in the analysis of the CMB data. Before removing a source, however, it is necessary to detect it and source extraction consists of a delicate preliminary operation. In the literature, various techniques have been presented to detect point-sources in the sky maps. The most sophisticated ones exploit the multi-frequency nature of the observations that is typical of the CMB experiments. These techniques have "optimal" theoretical properties and, at least in principle, are capable of remarkable performances. Actually, they are rather difficult to use and this deteriorates the quality of the obtainable results. In this paper, we present a new technique, the "weighted matched filter" (WMF), that is quite simple to use and hence more robust in practical applications. Such technique shows particular efficiency in the detection of sources whose spectra have a slope different from zero. We apply this method to three Southern Hemisphere sky regions - each with an area of 400 square degrees - of the seven years Wilkinson Microwave Anisotropy Probe (WMAP) maps and compare the resulting sources with those of the two seven-year WMAP point-sources catalogues. In these selected regions we find seven additional sources not previously listed in WMAP catalogues and discuss their most likely identification and spectral properties.

We study the dust properties of galaxies in the redshift range 0.1<z<2.8 observed by the Herschel Space Observatory in the field of the Great Observatories Origins Deep Survey-North as part of PEP and HerMES key programmes. Infrared (IR) luminosity (L_IR) and dust temperature (T_dust) of galaxies are derived from the spectral energy distribution (SED) fit of the far-infrared (FIR) flux densities obtained with PACS and SPIRE instruments onboard Herschel. As a reference sample, we also obtain IR luminosities and dust temperatures of local galaxies at z<0.1 using AKARI and IRAS data in the field of the Sloan Digital Sky Survey. We compare the L_IR-T_dust relation between the two samples and find that: the median T_dust of Herschel-selected galaxies at z>0.5 with L_IR>5x10^{10} L_\odot, appears to be 2-5 K colder than that of AKARI-selected local galaxies with similar luminosities; and the dispersion in T_dust for high-z galaxies increases with L_IR due to the existence of cold galaxies that are not seen among local galaxies. We show that this large dispersion of the L_IR-T_dust relation can bridge the gap between local star-forming galaxies and high-z submillimeter galaxies (SMGs). We also find that three SMGs with very low T_dust (<20 K) covered in this study have close neighbouring sources with similar 24-\mum brightness, which could lead to an overestimation of FIR/(sub)millimeter fluxes of the SMGs.

Ly-alpha emitters (LAEs) are seen everywhere in the redshift domain from local to z~7. Far-infrared (FIR) counterparts of LAEs at different epochs could provide direct clues on dust content, extinction, and spectral energy distribution (SED) for these galaxies. We search for FIR counterparts of LAEs that are optically detected in the GOODS-North field at redshift z~2.2 using data from the Herschel Space Telescope with the Photodetector Array Camera and Spectrometer (PACS). The LAE candidates were isolated via color-magnitude diagram using the medium-band photometry from the ALHAMBRA Survey, ancillary data on GOODS-North, and stellar population models. According to the fitting of these spectral synthesis models and FIR/optical diagnostics, most of them seem to be obscured galaxies whose spectra are AGN-dominated. From the analysis of the optical data, we have observed a fraction of AGN or composite over source total number of ~0.75 in the LAE population at z~2.2, which is marginally consistent with the fraction previously observed at z=2.25 and even at low redshift (0.2<z<0.45), but significantly different from the one observed at redshift ~3, which could be compatible either with a scenario of rapid change in the AGN fraction between the epochs involved or with a non detection of obscured AGN in other z=2-3 LAE samples due to lack of deep FIR observations. We found three robust FIR (PACS) counterparts at z~2.2 in GOODS-North. This demonstrates the possibility of finding dust emission in LAEs even at higher redshifts.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We present first insights into the far-IR properties for a sample of IRAC and MIPS-24um detected Lyman Break Galaxies (LBGs) at z ~ 3, as derived from observations in the northern field of the Great Observatories Origins Survey (GOODS-N) carried out with the PACS instrument on board the Herschel Space Observatory. Although none of our galaxies are detected by Herschel, we employ a stacking technique to construct, for the first time, the average spectral energy distribution of infrared luminous LBGs from UV to radio wavelengths. We derive a median IR luminosity of L_{IR} = 1.6 x 10^12 Lo, placing the population in the class of ultra luminous infrared galaxies (ULIRGs). Complementing our study with existing multi-wavelength data, we put constraints on the dust temperature of the population and find that for their L_{IR}, MIPS-LBGs are warmer than submm-luminous galaxies while they fall in the locus of the L_{IR}-T_{d} relation of the local ULIRGs. This, along with estimates based on the average SED, explains the marginal detection of LBGs in current sub-mm surveys and suggests that these latter studies introduce a bias towards the detection of colder ULIRGs in the high-z universe, while missing high-z ULIRGS with warmer dust.

We present a search for Herschel-PACS counterparts of dust-obscured, high-redshift objects previously selected at submillimeter and millimeter wavelengths in the Great Observatories Origins Deep Survey North field. We detect 22 of 56 submillimeter galaxies (39%) with a SNR of >=3 at 100 micron down to 3.0 mJy, and/or at 160 micron down to 5.7 mJy. The fraction of SMGs seen at 160 micron is higher than that at 100 micron. About 50% of radio-identified SMGs are associated with PACS sources. We find a trend between the SCUBA/PACS flux ratio and redshift, suggesting that these flux ratios could be used as a coarse redshift indicator. PACS undetected submm/mm selected sources tend to lie at higher redshifts than the PACS detected ones. A total of 12 sources (21% of our SMG sample) remain unidentified and the fact that they are blank fields at Herschel-PACS and VLA 20 cm wavelength may imply higher redshifts for them than for the average SMG population (e.g., z>3-4). The Herschel-PACS imaging of these dust-obscured starbursts at high-redshifts suggests that their far-infrared spectral energy distributions have significantly different shapes than template libraries of local infrared galaxies.

The periodogram is a popular tool that tests whether a signal consists only of noise or if it also includes other components. The main issue of this method is to define a critical detection threshold that allows identification of a component other than noise, when a peak in the periodogram exceeds it. In the case of signals sampled on a regular time grid, determination of such a threshold is relatively simple. When the sampling is uneven, however, things are more complicated. The most popular solution in this case is to use the "Lomb-Scargle" periodogram, but this method can be used only when the noise is the realization of a zero-mean, white (i.e. flat-spectrum) random process. In this paper, we present a general formalism based on matrix algebra, which permits analysis of the statistical properties of a periodogram independently of the characteristics of noise (e.g. colored and/or non-stationary), as well as the characteristics of sampling.

We use deep observations taken with the Photodetector Array Camera and Spectrometer (PACS), on board the Herschel satellite as part of the PACS evolutionary probe (PEP) guaranteed project along with submm ground-based observations to measure the dust mass of a sample of high-z submillimeter galaxies (SMGs). We investigate their dust content relative to their stellar and gas masses, and compare them with local star-forming galaxies. High-z SMGs are dust rich, i.e. they have higher dust-to-stellar mass ratios compared to local spiral galaxies (by a factor of 30) and also compared to local ultraluminous infrared galaxies (ULIRGs, by a factor of 6). This indicates that the large masses of gas typically hosted in SMGs have already been highly enriched with metals and dust. Indeed, for those SMGs whose gas mass is measured, we infer dust-to-gas ratios similar or higher than local spirals and ULIRGs. However, similarly to other strongly star-forming galaxies in the local Universe and at high-z, SMGs are characterized by gas metalicities lower (by a factor of a few) than local spirals, as inferred from their optical nebular lines, which are generally ascribed to infall of metal-poor gas. This is in contrast with the large dust content inferred from the far-IR and submm data. In short, the metalicity inferred from the dust mass is much higher (by more than an order of magnitude) than that inferred from the optical nebular lines. We discuss the possible explanations of this discrepancy and the possible implications for the investigation of the metalicity evolution at high-z.

Gravitational lensing by massive galaxy clusters allows study of the population of intrinsically faint infrared galaxies that lie below the sensitivity and confusion limits of current infrared and submillimeter telescopes. We present ultra-deep PACS 100 and 160 microns observations toward the cluster lens Abell 2218, to penetrate the Herschel confusion limit. We derive source counts down to a flux density of 1 mJy at 100 microns and 2 mJy at 160 microns, aided by strong gravitational lensing. At these levels, source densities are 20 and 10 beams/source in the two bands, approaching source density confusion at 160 microns. The slope of the counts below the turnover of the Euclidean-normalized differential curve is constrained in both bands and is consistent with most of the recent backwards evolutionary models. By integrating number counts over the flux range accessed by Abell 2218 lensing (0.94-35 mJy at 100 microns and 1.47-35 mJy at 160 microns, we retrieve a cosmic infrared background (CIB) surface brightness of ~8.0 and ~9.9 nW m^-2 sr^-1, in the respective bands. These values correspond to 55% (+/- 24%) and 77% (+/- 31%) of DIRBE direct measurements. Combining Abell 2218 results with wider/shallower fields, these figures increase to 62% (+/- 25%) and 88% (+/- 32%) CIB total fractions, resolved at 100 and 160 microns, disregarding the high uncertainties of DIRBE absolute values.

The Photodetector Array Camera and Spectrometer (PACS) is one of the three science instruments on ESA's far infrared and submillimetre observatory. It employs two Ge:Ga photoconductor arrays (stressed and unstressed) with 16x25 pixels, each, and two filled silicon bolometer arrays with 16x32 and 32x64 pixels, respectively, to perform integral-field spectroscopy and imaging photometry in the 60-210\mu\ m wavelength regime. In photometry mode, it simultaneously images two bands, 60-85\mu\ m or 85-125\mu\m and 125-210\mu\ m, over a field of view of ~1.75'x3.5', with close to Nyquist beam sampling in each band. In spectroscopy mode, it images a field of 47"x47", resolved into 5x5 pixels, with an instantaneous spectral coverage of ~1500km/s and a spectral resolution of ~175km/s. We summarise the design of the instrument, describe observing modes, calibration, and data analysis methods, and present our current assessment of the in-orbit performance of the instrument based on the Performance Verification tests. PACS is fully operational, and the achieved performance is close to or better than the pre-launch predictions.

The star formation rate (SFR) is a key parameter in the study of galaxy evolution. The accuracy of SFR measurements at z~2 has been questioned following a disagreement between observations and theoretical models. The latter predict SFRs at this redshift that are typically a factor 4 or more lower than the measurements. We present star-formation rates based on calorimetric measurements of the far-infrared (FIR) luminosities for massive 1.5<z<2.5, normal star-forming galaxies (SFGs), which do not depend on extinction corrections and/or extrapolations of spectral energy distributions. The measurements are based on observations in GOODS-N with the Photodetector Array Camera & Spectrometer (PACS) onboard Herschel, as part of the PACS Evolutionary Probe (PEP) project, that resolve for the first time individual SFGs at these redshifts at FIR wavelengths. We compare FIR-based SFRs to the more commonly used 24 micron and UV SFRs. We find that SFRs from 24 micron alone are higher by a factor of ~4-7.5 than the true SFRs. This overestimation depends on luminosity: gradually increasing for log L(24um)>12.2 L_sun. The SFGs and AGNs tend to exhibit the same 24 micron excess. The UV SFRs are in closer agreement with the FIR-based SFRs. Using a Calzetti UV extinction correction results in a mean excess of up to 0.3 dex and a scatter of 0.35 dex from the FIR SFRs. The previous UV SFRs are thus confirmed and the mean excess, while narrowing the gap, is insufficient to explain the discrepancy between the observed SFRs and simulation predictions.

The Herschel Space Observatory enables us to accurately measure the bolometric output of starburst galaxies and active galactic nuclei (AGN) by directly sampling the peak of their far-infrared (IR) emission. Here we examine whether the spectral energy distribution (SED) and dust temperature of galaxies have strongly evolved since z~2.5. We use Herschel deep extragalactic surveys from 100 to 500um to compute total IR luminosities in galaxies down to the faintest levels, using PACS and SPIRE in the GOODS-North field (PEP and HerMES key programs). We show that measurements in the SPIRE bands can be used below the statistical confusion limit if information at higher spatial resolution is used to identify isolated galaxies whose flux is not boosted by bright neighbors. Below z~1.5, mid-IR extrapolations are correct for star-forming galaxies with a dispersion of only 40% (0.15dex), therefore similar to z~0 galaxies. This narrow distribution is puzzling when considering the range of physical processes that could have affected the SED of these galaxies. Extrapolations from only one of the 160um, 250um or 350um bands alone tend to overestimate the total IR luminosity. This may be explained by the lack of far-IR constraints around and above ~150um (rest-frame) on local templates. We also note that the dust temperature of luminous IR galaxies around z~1 is mildly colder by 10-15% than their local analogs and up to 20% for ULIRGs at z~1.6. Above z=1.5, distant galaxies are found to exhibit a substantially larger mid- over far-IR ratio, which could either result from stronger broad emission lines or warm dust continuum heated by a hidden AGN. Two thirds of the AGNs identified in the field with a measured redshift exhibit the same behavior as purely star-forming galaxies. Hence a large fraction of AGNs harbor star formation at very high SFR and in conditions similar to purely star-forming galaxies.

Sensitive Herschel far-infrared observations can break degeneracies that were inherent to previous studies of star formation in high-z AGN hosts. Combining PACS 100 and 160um observations of the GOODS-N field with 2Msec Chandra data, we detect ~20% of X-ray AGN individually at >3sig. The host far-infrared luminosity of AGN with L2-10~10^43erg/s increases with redshift by an order of magnitude from z=0 to z~1. In contrast, there is little dependence of far-infrared luminosity on AGN luminosity, for L2-10<~10^44erg/s AGN at z>~1. We do not find a dependence of far-infrared luminosity on X-ray obscuring column, for our sample which is dominated by L2-10<10^44erg/s AGN. In conjunction with properties of local and luminous high-z AGN, we interpret these results as reflecting the interplay between two paths of AGN/host coevolution. A correlation of AGN luminosity and host star formation is traced locally over a wide range of luminosities and also extends to luminous high z AGN. This correlation reflects an evolutionary connection, likely via merging. For lower AGN luminosities, star formation is similar to that in non-active massive galaxies and shows little dependence on AGN luminosity. The level of this secular, non-merger driven star formation increasingly dominates over the correlation at increasing redshift.

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

We exploit the deepest existing far-infrared (FIR) data obtained so far by Herschel at 100 and 160 um in the GOODS-N, as part of the PACS Evolutionary Probe (PEP) survey, to derive for the first time the evolution of the rest-frame 60-um, 90-um, and total IR luminosity functions (LFs) of galaxies and AGNs from z=0 to unprecedented high redshifts (z~2-3). The PEP LFs were computed using the 1/Vmax method. The FIR sources were classified by means of a detailed broad- band SED-fitting analysis and spectral characterisation. Based on the best-fit model results, k-correction and total IR (8-1000 um) luminosity were obtained for each source. LFs (monochromatic and total) were then derived for various IR populations separately in different redshift bins and compared to backward evolution model predictions. We detect strong evolution in the LF to at least z~2. Objects with SEDs similar to local spiral galaxies are the major contributors to the star formation density (SFD) at z< 0.3, then, as redshift increases, moderate SF galaxies - most likely containing a low-luminosity AGN - start dominating up to z ~= 1.5. At >1.5 the SFD is dominated by the contributions of starburst galaxies. In agreement with previous findings, the comoving IR LD derived from our data evolves approximately as (1 + z)^(3.8+/-0.3) up to z~1, there being some evidence of flattening up to z~2.

We use deep observations obtained with the Photodetector Array Camera and Spectrometer (PACS) onboard the Herschel space observatory to study the far-infrared (FIR) properties of submm and optically faint radio galaxies (SMGs and OFRGs). From literature we compiled a sample of 35 securely identified SMGs and nine OFRGs located in the GOODS-N and the A2218 fields. This sample is cross-matched with our PACS 100 um and 160 um multi-wavelength catalogs. About half of the galaxies in our sample are detected with PACS. The dust temperatures and the infrared luminosities of our galaxies are derived by fitting their PACS and SCUBA 850 um (only the upper limits for the OFRGs) flux densities with a single modified (beta=1.5) black body function. The median dust temperature of our SMG sample is T=36+/-8K while for our OFRG sample it is T=47+/-3K. For both samples, median dust temperatures derived from Herschel data agree well with previous estimates. In particular, Chapman et al. (2005) found a dust temperature of T=36+/-7K for a large sample of SMGs assuming the validity of the FIR/radio correlation. The agreement between our studies confirms that the local FIR/radio correlation effectively holds at high redshift even though we find <q>=2.17+/-0.19, a slightly lower value than that observed in local systems. The median IR luminosities of SMGs and OFRGs are 4.6*10^12 Lsun and 2.6*10^12 Lsun, respectively. We note that for both samples the IR luminosity estimates from the radio part of the spectral energy distribution are accurate, while estimates from the mid-IR are considerably (x3) more uncertain. Our observations confirm the remarkably high luminosities of SMGs and thus imply median star-formation rates of 960Msun yr^-1 for SMGs with S(850um)>5mJy and 460Msun yr^-1 for SMGs with S(850um)>2mJy, assuming a Chabrier IMF and no dominant AGN contribution to the far-infrared luminosity.

We exploit deep observations of the GOODS-N field taken with PACS, on board of Herschel, as part of the PEP guaranteed time, to study the link between star formation and stellar mass in galaxies to z~2. Starting from a stellar mass-selected sample of ~4500 galaxies with mag[4.5mu]<23 (AB), we identify ~350 objects with a PACS detection at 100 or 160mu and ~1500 with only Spitzer 24 mu counterpart. Stellar masses and total IR luminosities (LIR) are estimated by fitting the SEDs. Consistently with other Herschel results, we find that LIR based only on 24 mu data is overestimated by a median factor ~1.8 at z~2, whereas it is underestimated (with our approach) up to a factor ~1.6 at 0.5<z<1.0. We then exploit this calibration to correct LIR based on the MIPS fluxes. These results clearly show how Herschel is fundamental to constrain LIR, and hence the SFR, of high redshift galaxies. Using the galaxies detected with PACS (and/or MIPS), we investigate the existence and evolution of the relations between the star formation rate (SFR), the specific star formation rate (SSFR=SFR/mass) and the stellar mass. Moreover, in order to avoid selection effects, we also repeat this study through a stacking analysis on the PACS images to fully exploit the far-IR information also for the Herschel and Spitzer undetected subsamples. We find that the SSFR-mass relation steepens with redshift, being almost flat at z<1.0 and reaching a slope of alpha=-0.50^(+0.13)_(-0.16) at z~2, at odds with recent works based on radio-stacking analysis at the same redshift. The mean SSFR of galaxies increases with redshift, by a factor ~15 for massive M>10^(11) (M_sun) galaxies from z=0 to z=2, and seems to flatten at z>1.5 in this mass range. Moreover, the most massive galaxies have the lowest SSFR at any z, implying that they have formed their stars earlier and more rapidly than their low mass counterparts.

The constituents of the cosmic IR background (CIB) are studied at its peak wavelengths (100 and 160 um) by exploiting Herschel/PACS observations of the GOODS-N, Lockman Hole, and COSMOS fields in the PACS Evolutionary Probe (PEP) guaranteed-time survey. The GOODS-N data reach 3 sigma depths of ~3.0 mJy at 100 um and ~5.7 mJy at 160 um. At these levels, source densities are 40 and 18 beams/source, respectively, thus hitting the confusion limit at 160 um. Differential number counts extend from a few mJy up to 100-200 mJy, and are approximated as a double power law, with the break lying between 5 and 10 mJy. The available ancillary information allows us to split number counts into redshift bins. At z<=0.5 we isolate a class of luminous sources (L(IR)~1e11 Lsun), whose SEDs resemble late-spiral galaxies, peaking at ~130 um restframe and significantly colder than what is expected on the basis of pre-Herschel models. By integrating number counts over the whole covered flux range, we obtain a surface brightness of 6.36 +/- 1.67 and 6.58 +/-1.62 [nW m^-2 sr^-1] at 100 and 160 um, resolving ~45% and ~52% of the CIB, respectively. When stacking 24 um sources, the inferred CIB lies within 1.1 sigma and 0.5 sigma from direct measurements in the two bands, and fractions increase to 50% and 75%.Most of this resolved CIB fraction was radiated at z<=1.0, with 160 um sources found at higher redshift than 100 um ones.

We set out to determine the ratio, q(IR), of rest-frame 8-1000um flux, S(IR), to monochromatic radio flux, S(1.4GHz), for galaxies selected at far-IR and radio wavelengths, to search for signs that the ratio evolves with redshift, luminosity or dust temperature, and to identify any far-IR-bright outliers - useful laboratories for exploring why the far-IR/radio correlation is generally so tight when the prevailing theory suggests variations are almost inevitable. We use flux-limited 250-um and 1.4-GHz samples, obtained in GOODS-N using Herschel (HerMES; PEP) and the VLA. We determine bolometric IR output using ten bands spanning 24-1250um, exploiting data from PACS and SPIRE, as well as Spitzer, SCUBA, AzTEC and MAMBO. We also explore the properties of an L(IR)-matched sample, designed to reveal evolution of q(IR) with z, spanning log L(IR) = 11-12 L(sun) and z=0-2, by stacking into the radio and far-IR images. For 1.4-GHz-selected galaxies, we see tentative evidence of a break in the flux ratio, q(IR), at L(1.4GHz) ~ 10^22.7 W/Hz, where AGN are starting to dominate the radio power density, and of weaker correlations with z and T(d). From our 250-um-selected sample we identify a small number of far-IR-bright outliers, and see trends of q(IR) with L(1.4GHz), L(IR), T(d) and z, noting that some of these are inter-related. For our L(IR)-matched sample, there is no evidence that q(IR) changes significantly as we move back into the epoch of galaxy formation: we find q(IR) goes as (1+z)^gamma, where gamma = -0.04 +/- 0.03 at z=0-2; however, discounting the least reliable data at z<0.5 we find gamma = -0.26 +/- 0.07, modest evolution which may be related to the radio background seen by ARCADE2, perhaps driven by <10uJy radio activity amongst ordinary star-forming galaxies at z>1.

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.

Although the ``internal linear'' combination method (ILC) is a technique widely used for the separation of the Cosmic Microwave Background signal from the Galactic foregrounds, its characteristics are not yet well defined. This can lead to misleading conclusions about the actual potentialities and limits of such approach in real applications. Here we discuss briefly some facts about ILC that to our knowledge are not fully worked out in literature and yet have deep effects in the interpretation of the results.

Optically very faint R>25.5 sources detected by the Spitzer Space Telescope at 24um represent a very interesting population at redshift z [1.5-3]. They exhibit strong clustering properties, implying that they are hosted by very massive halos, and their mid-IR emission could be powered by either dust-enshrouded star-formation and/or by an obscured AGN. We report observations carried out with the MAMBO array at the IRAM 30m antenna on Pico Veleta of a candidate protocluster with five optically obscured sources selected from the 24um Spitzer sample of the First Look Survey. Interestingly, these sources appear to lie on a high density filament aligned with the two radio jets of an AGN. Four out of five of the observed sources were detected. We combine these measurements with optical, infrared and radio observations to probe}the nature of the candidate protocluster members. Our preliminary conclusions can be summarized as follows: the Spectral Energy Distributions of all sources include both AGN and starburst contributions; the AGN contribution to the bolometric luminosities ranges between 14 and 26% of the total. Such a contribution is enough for the AGN to dominate the emission at 5.8, 8 and 24um, while the stellar component, inferred from SED fitting, prevails at 1.25mm and at lambda<4.5um. The present analysis suggests a coherent interplay at high-z between extended radio activity and the development of filamentary large-scale structures.

The detection of point-sources in microwave maps is a critical step in the analysis of the Cosmic Microwave Background (CMB) data. If not properly removed, these sources have adverse effects on the estimation of the power-spectrum and/or the test of Gaussianity of the CMB component. The techniques to extract point sources from an observed image used up to now do not show a general consensus about their real performance and properties. Their characteristics have been studied essentially through numerical simulations based on semi-empirical models of the CMB and the Galactic foreground. Such models often have different levels of sophistication and/or are based on different physical assumptions (e.g. the number of Galactic components and level of the noise). Moreover, the application of a given technique to a set of data (either simulated or experimental) requires the tuning of one or more parameters that unavoidably is a subjective operation. Hence, a reliable comparison is difficult. What is missing is a statistical analysis of the properties of the proposed methodologies. This is the aim of the present paper. The statistical properties of the detection techniques in the context of two different criteria, i.e. the Neyman-Pearson criterion and the maximization of the signal-to-noise ratio, are analyzed through an analytical approach. Various examples, techniques and the conditions under which they work satisfactorily are shown together with their limits and implementation details. Exploiting some a priori information, it is possible to develop simple algorithms with performances similar to those of more sophisticated but complex techniques. In this respect, a detection algorithm, useful in future ground-based experiments as, for example the Atacama Large Millimeter/submillimeter Array, is presented.

Foreground contamination is the fundamental hindrance to the cosmic microwave background (CMB) signals and its separation from it represents a fundamental question in Cosmology. One of the most popular algorithm used to disentangle foregrounds from the CMB signals is the "internal linear combination" method (ILC). In its original version, this technique is applied directly to the observed maps. In recent literature, however, it is suggested that in the harmonic (Fourier) domain it is possible to obtain better results since a separation can be attempted where the various Fourier frequencies are given different weights. This is seen as a useful characteristic in the case of noisy data. Here, we argue that the benefits of using such an approach are overestimated. Better results can be obtained if a classic procedure is adopted where data are filtered before the separation is carried out.

The Atacama Large Millimetre/submillimetre Array (ALMA) is one of the largest and most complicated observatories ever built. Constructing and operating an observatory at high altitude (5000m) in a cost effective and safe manner, with minimal effect on the environment creates interesting challenges. Since the array will have to adapt quickly to prevailing weather conditions, ALMA will be operated exclusively in service mode. By the time of full science operations, the fundamental ALMA data product shall be calibrated, deconvolved data cubes and images, but raw data and data reduction software will be made available to users as well. User support is provided by the ALMA Regional Centres (ARCs) located in Europe, North America and Japan. These ARCs constitute the interface between the user community and the ALMA observatory in Chile. For European users the European ARC is being set up as a cluster of nodes located throughout Europe, with the main centre at the ESO Headquarters in Garching. The main centre serves as the access portal and in synergy with the distributed network of ARC nodes, the main aim of the ARC is to optimize the ALMA science output and to fully exploit this unique and powerful facility. The aim of this article is to introduce the process of proposing for observing time, subsequent execution of the observations, obtaining and processing of the data in the ALMA epoch. The complete end-to-end process of the ALMA data flow from the proposal submission to the data delivery is described.

The effects of the in-flight behaviour of the bolometer arrays of the Herschel/PACS instrument under impacts of Galactic cosmic rays are explored. This instrument is part of the ESA-Herschel payload, which will be launched at the end of 2008 and will operate at the Lagrangian L2 point of the Sun-Earth system. We find that the components external to the detectors (the spacecraft, the cryostat, the PACS box, collectively referred to as the `shield') are the major source of secondary events affecting the detector behaviour. The impacts deposit energy on the bolometer chips and influence the behaviour of nearby pixels. 25% of hits affect the adjacent pixels. The energy deposited raises the bolometer temperature by a factor ranging from 1 to 6 percent of the nominal value. We discuss the effects on the observations and compare simulations with laboratory tests.

AIMS: The separation of foreground contamination from cosmic microwave background (CMB) observations is one of the most challenging and important problem of digital signal processing in Cosmology. In literature, various techniques have been presented, but no general consensus about their real performances and properties has been reached. This is due to the characteristics of these techniques that have been studied essentially through numerical simulations based on semi-empirical models of the CMB and the Galactic foregrounds. Such models often have different level of sophistication and/or are based on different physical assumptions (e.g., the number of the Galactic components and the level of the noise). Hence, a reliable comparison is difficult. What actually is missing is a statistical analysis of the properties of the proposed methodologies. Here, we consider the "Internal Linear Combination" method (ILC) which, among the separation techniques, requires the smallest number of "a priori" assumptions. This feature is of particular interest in the context of the CMB polarization measurements at small angular scales where the lack of knowledge of the polarized backgrounds represents a serious limit. METHODS: The statistical characteristics of ILC are examined through an analytical approach and the basic conditions are fixed in a way to work satisfactorily. RESULTS: ILC provides satisfactory results only under rather restrictive conditions. This is a critical fact to take into consideration in planning the future ground-based observations (e.g., with ALMA) where, contrary to the satellite experiments, there is the possibility to have a certain control of the experimental conditions.

We present Spitzer rest-frame mid-infrared spectroscopy of twelve z~2 mm-bright type 1 QSOs, selected from unlensed and lensed QSO samples and covering a range of AGN optical luminosities L_5100=10^45 to 10^47 erg/s. On top of the AGN continuum, we detect PAH emission from luminous star formation in nine objects individually as well as in the composite spectrum for the full sample. PAH luminosity and rest frame far-infrared luminosity correlate and extend the similar correlation for lower luminosity local QSOs. This provides strong evidence for intense star formation in the hosts of these mm-bright QSOs, sometimes exceeding 1000 Msun/yr and dominating their rest frame far-infrared emission. The PAH-based limit on star formation rates is lower for luminous z~2 QSOs that are not preselected for their mm emission. Partly dependent on systematic changes of the AGN dust covering factor and the dust spectral energy distribution of the AGN proper, the spectral energy distributions of mm-faint high-z QSOs may be AGN dominated out to rest frame far-infrared wavelengths. Towards the most luminous high-z QSOs, there is a flattening of the relation between star formation and AGN luminosity that is observed for lower redshift QSOs. No QSO in our sample has a PAH-measured star formation rate in excess of 3000 Msun/yr.

AIMS: One of the most challenging and important problem of digital signal processing in Cosmology is the separation of foreground contamination from cosmic microwave background (CMB). This problem becomes even more difficult in situations, as the CMB polarization observations, where the amount of available "a priori" information is limited. In this case, it is necessary to resort to the "blind separation" methods. One important member of this class is represented by the "Independent Components Analysis" (ICA). In its original formulation, this method has various interesting characteristics, but also some limits. One of the most serious is the difficulty to take into account any information available in advance. In particular, ICA is not able to exploit the fact that emission of CMB is the same at all the frequencies of observations. Here, we show how to deal with this question. The connection of the proposed methodology with the "Internal Linear Composition" (ILC) technique is also illustrated. METHODS: A modification of the classic ICA approach is presented and its characteristics are analyzed both analytically and by means of numerical experiments. RESULTS: The modified version of ICA appears to provide more stable results and of better quality.

This paper analyses the radio properties of a subsample of optically obscured (R>25.5) galaxies observed at 24um by the Spitzer Space Telescope within the First Look Survey. 96 F[24um]>0.35 mJy objects out of 510 are found to have a radio counterpart at 1.4 GHz, 610 MHz or at both frequencies respectively down to ~40uJy and ~200uJy. IRAC photometry sets the majority of them in the redshift interval z [1-3] and allows for a broad distinction between AGN-dominated galaxies (~47% of the radio-identified sample) and systems powered by intense star-formation (~13%), the remaining objects being impossible to classify. The percentage of radio identifications is a strong function of 24um flux. The radio number counts at both radio frequencies suggest that the physical process(es) responsible for radio activity in these objects have a common origin regardless of whether the source shows mid-IR emission compatible with being an obscured AGN or a star-forming galaxy. We also find that both candidate AGN and star-forming systems follow (although with a large scatter) the relationship between 1.4 GHz and 24um fluxes reported by Appleton et al. (2004) which identifies sources undergoing intense star formation activity. On the other hand, the inferred radio spectral indices alpha indicate that a large fraction of objects in our sample (~60% of all galaxies with estimated alpha) may belong to the population of Ultra Steep Spectrum (USS) Sources, typically 'frustrated' radio-loud AGN. We interpret our findings as a strong indication for concurrent AGN and star-forming activity, whereby the 1.4 GHz flux is of thermal origin, while that at 610 GHz mainly stems from the nuclear source.

We report the first detection of the 6.2micron and 7.7micron infrared `PAH' emission features in the spectrum of a high redshift QSO, from the Spitzer-IRS spectrum of the Cloverleaf lensed QSO (H1413+117, z~2.56). The ratio of PAH features and rest frame far-infrared emission is the same as in lower luminosity star forming ultraluminous infrared galaxies and in local PG QSOs, supporting a predominantly starburst nature of the Cloverleaf's huge far-infrared luminosity (5.4E12 Lsun, corrected for lensing). The Cloverleaf's period of dominant QSO activity (Lbol ~ 7E13 Lsun) is coincident with an intense (star formation rate ~1000 Msun/yr) and short (gas exhaustion time ~3E7yr) star forming event.

Massive tori with $\approx 10^{8-9}M_{\odot}$ are predicted to extend on $\sim $100 pc scale around the centre of elliptical galaxy progenitors by a model of a supermassive black hole (SMBH) growth coeval to the spheroidal population of the host galaxy. Direct detection of such massive tori would cast light on a key physical condition that allows the rapid growth of SMBHs and the appearance of QSOs at high redshift. For this reason, we examine the detectability of such structures at substantial redshift with the Atacama Large Millimeter Array (ALMA). We propose that submillimeter galaxies (SMGs) are the best targets to test our predictions. In order to assess the observational feasibility, we estimate the expected number counts of SMGs with massive tori and check the detectability with the ALMA instrument, the unique facility which can resolve the central region of high redshift objects. Our work shows that ALMA will be able to resolve and detect high-$J$ ($J >$ 4) CO emissions from $\sim$100 pc scale extended massive tori up to $z\approx2$. Observations of lensed SMGs will yield excellent spatial resolution, allowing even to resolve their massive tori at higher redshift. We discuss further the detectability of the HCN molecule, as a better tracer of the high density gas expected in such tori. The final goal of these kind of observations is to pinpoint possible physical mechanisms that storage in the very central galactic regions very large amount gas on timescale of several 10$^{8}$ yr.

The Cosmic Microwave Background (CMB) physics can be used to constrain the dark energy dynamics: B modes of the polarization of the diffuse CMB emission as well as the polarized signal towards clusters of galaxies are sensitive to the Hubble expansion rate and thus the dark energy abundance in the early stage of cosmic acceleration. The first effect is sourced by gravitational lensing of large scale cosmic structures, the second is due to scattering of the primary CMB temperature anisotropy quadrupole by free electrons in cluster plasma. We are investigating the capabilities of ALMA to detect these effects and constrain the high-redshift dark energy abundance through measurements of sub-arcminute CMB anisotropies.

We present results of simulations performed with the Geant4 software code of the effects of Galactic Cosmic Ray impacts on the photoconductor arrays of the PACS instrument. This instrument is part of the ESA-Herschel payload, which will be launched in late 2007 and will operate at the Lagrangian L2 point of the Sun-Earth system. Both the Satellite plus the cryostat (the shield) and the detector act as source of secondary events, affecting the detector performance. Secondary event rates originated within the detector and from the shield are of comparable intensity. The impacts deposit energy on each photoconductor pixel but do not affect the behaviour of nearby pixels. These latter are hit with a probability always lower than 7%. The energy deposited produces a spike which can be hundreds times larger than the noise. We then compare our simulations with proton irradiation tests carried out for one of the detector modules and follow the detector behaviour under 'real' conditions.

The arguments presented by Gonzalez-Nuevo et. al (2006) in favour of the Mexican Hat Wavelet Family (MHWF) are critically discussed here. These authors allege the optimal properties of this new class of filters in the detection of point sources embedded in a noise background but their claim is not based upon a solid mathematical foundation and proof.

The Kluzniak & Abramowicz model explains high frequency, double peak, "3:2" QPOs observed in neutron star and black hole sources in terms of a non-linear parametric resonance between radial and vertical epicyclic oscillations of an almost Keplerian accretion disk. The 3:2 ratio of epicyclic frequencies occurs only in strong gravity. Rebusco (2004) and Horak (2004) studied the model analytically: they proved that a small forcing may indeed excite the parametric 3:2 resonance, but they have not explained the physical nature of the forcing. Here we integrate their equations numerically, dropping the ad hoc forcing, and adding instead a stochastic term to mimic the action of the very complex processes that occur in disks as, for example, MRI turbulence. We demonstrate that the presence of the stochastic term triggers the resonance in epicyclic oscillations of nearly Keplerian disks, and influences their pattern.

In this note we re-propose the arguments presented in Vio (2005) examining the superiority of the 'bi-parametric scale adaptive filter' (BSAF) when dealing with source detection as claimed by (Lopez-Caniego 2005a) and confirmed in Lopez-Caniego (2005b). We suggest here that the dispute can be easily solved if these authors provide the community with a detailed derivation of a basic equation.

In this note we stress the necessity of a careful check of both the theoretical arguments and the numerical experiments used by Lopez-Caniego et al. (2005) to support the superior performances of the biparametric scale adaptive filter (BSAF) with respect to the classic matched filter (MF) in the detection of sources on a random Gaussian background.

In this draft we comment on the results concerning the performances of matched filters, scale adaptive filters and Mexican hat wavelet that recently appeared in literature in the context of point source detection in Cosmic Microwave Background maps. In particular, we show that, contrary to what has been claimed, the use of the matched filters still appear to be the most reliable and efficient method to disantangle point sources from the backgrounds, even when using detection criterion that, differently from the classic $n\sigma$ thresholding rule, takes into account not only the height of the peaks in the signal corresponding to the candidate sources but also their curvature.

In this second paper in a series dedicated to developing efficient numerical techniques for the deblurring Cosmic Microwave Background (CMB) maps, we consider the case of asymmetric point spread functions (PSF). Although conceptually this problem is not different from the symmetric case, there are important differences from the computational point of view because it is no longer possible to use some of the efficient numerical techniques that work with symmetric PSFs. We present procedures that permit the use of efficient techniques even when this condition is not met. In particular, two methods are considered: a procedure based on a Kronecker approximation technique that can be implemented with the numerical methods used with symmetric PSFs but that has the limitation of requiring only mildly asymmetric PSFs. The second is a variant of the classic Tikhonov technique that works even with very asymmetric PSFs but that requires discarding the edges of the maps. We provide details for efficient implementations of the algorithms. Their performance is tested on simulated CMB maps.

Averaged spectral energy distributions (SEDs) of active and starburst galaxies from the 12 micron sample in the Local Universe and Quasars, from an optically selected sample at a mean redshift <z>=0.7, are built from optical/near-IR/far-IR (IRAS & ISO) photometric observations. These SEDs are then used to predict at various redshifts the number of Seyfert type 1 and type 2, starburst, normal galaxies, and quasars, that will be detected by future Space Missions dedicated to far-infrared and submillimeter astronomy, like SIRTF and Herschel. These predictions are then compared with the expected capabilities and detection limits of future deep far-IR surveys. Possible ways to identify AGN candidates on far-IR colour-colour plots for follow-up observations are then explored.

Digital deblurring of images is an important problem that arises in multifrequency observations of the Cosmic Microwave Background (CMB) where, because of the width of the point spread functions (PSF), maps at different frequencies suffer a different loss of spatial resolution. Deblurring is useful for various reasons: first, it helps to restore high frequency components lost through the smoothing effect of the instrument's PSF; second, emissions at various frequencies observed with different resolutions can be better studied on a comparable resolution; third, some map-based component separation algorithms require maps with similar level of degradation. Because of computational efficiency, deblurring is usually done in the frequency domain. But this approach has some limitations as it requires spatial invariance of the PSF, stationarity of the noise, and is not flexible in the selection of more appropriate boundary conditions. Deblurring in real space is more flexible but usually not used because of its high computational cost. In this paper (the first in a series on the subject) we present new algorithms that allow the use of real space deblurring techniques even for very large images. In particular, we consider the use of Tikhonov deblurring of noisy maps with applications to {\it PLANCK}. We provide details for efficient implementations of the algorithms. Their performance is tested on Gaussian and non-Gaussian simulated CMB maps, and PSFs with both circular and elliptical symmetry. Matlab code is made available.

We present the ISO far-IR photometry of a complete sub-sample of optically selected bright quasars belonging to two complete surveys selected through multicolour (U,B,V,R,I) techniques. The ISOPHOT camera on board of the ISO Satellite was used to target these quasars at wavelengths of 7.3, 11.5, 60, 100 and 160 micron. Almost two thirds of the objects were detected at least in one ISOPHOT band. The detection rate is independent of the source redshift, very likely due to the negative K-correction of the far-IR thermal emission. More than a half of the optically selected QSOs show significant emission between 4 and 100 micron in the quasar rest-frame. These fluxes have a very likely thermal origin, although in a few objects an additional contribution from a non-thermal component is plausible in the long wavelength bands. In a colour-colour diagram these objects span a wide range of properties from AGN-dominated to ULIRG-like. The far-IR composite spectrum of the quasar population presents a broad far-IR bump between 10 and 30 micron and a sharp drop at wavelengths greater than 100 micron in the quasar restframe. The amount of energy emitted in the far-IR, is on average a few times larger than that emitted in the blue and the ratio L(FIR)/L(B) increases with the bolometric luminosity. Objects with fainter blue magnitudes have larger ratios between the far-IR (wavelengths > 60 micron) fluxes and the blue band flux, which is attributed to extinction by dust around the central source. No relation between the blue absolute magnitude and the dust colour temperature is seen, suggesting that the dominant source of FIR energy could be linked to a concurrent starburst rather than to gravitational energy produced by the central engine.

We provide theoretical procedures and practical recipes to simulate non-Gaussian correlated, homogeneous random fields with prescribed marginal distributions and cross-correlation structure, either in a N-dimensional Cartesian space or on the celestial sphere. We illustrate our methods using far-infrared maps obtained with the Infrared Space Observatory. However, the methodology presented here can be used in other astrophysical applications that require modeling correlated features in sky maps, for example, the simulation of multifrequency sky maps where backgrounds, sources and noise are correlated and can be modeled by random fields.

We study the clustering properties of the galaxy clusters detectable for the Planck satellite due to their thermal Sunyaev-Zel'dovich effect. We take the past light-cone effect and the redshift evolution of both the underlying dark matter correlation function and the cluster bias factor into account. A theoretical mass-temperature relation allows us to convert the sensitivity limit of a catalogue into a minimum mass for the dark matter haloes hosting the clusters. We confirm that the correlation length is an increasing function of the sensitivity limits defining the survey. Using the expected characteristics of the Planck cluster catalogue, which will be a quite large and unbiased sample, we predict the two-point correlation function and power spectrum for different cosmological models. We show that the wide redshift distribution of the Planck survey, will allow to constrain the cluster clustering properties up to z=1. The dependence of our results on the main cosmological parameters (the matter density parameter, the cosmological constant and the normalisation of the density power-spectrum) is extensively discussed. We find that the future Planck clustering data place only mild constraints on the cosmological parameters, because the results depend on the physical characteristics of the intracluster medium, like the baryon fraction and the mass-temperature relation. Once the cosmological model and the Hubble constant are determined, the clustering data will allow a determination of the baryon fraction with an accuracy of few per cent.

New far-infrared photometry with ISOPHOT, onboard the Infrared Space Observatory, is presented for 58 galaxies with homogeneous published data for another 32 galaxies all belonging to the 12 micron galaxy sample. In total 29 Seyfert 1's, 35 Seyfert 2's and 12 starburst galaxies, about half of the 12 micron active galaxy sample, plus 14 normal galaxies for comparison. The ISO and the IRAS data are used to define color-color diagrams and spectral energy distributions (SED). Thermal dust emission at two temperatures (one cold at 15-30K and one warm at 50-70K) can fit the 60-200 micron SED, with a dust emissivity law proportional to the inverse square of the wavelength. Seyfert 1's and Seyfert 2's are indistinguishable longward of 100 micron, while, as already seen by IRAS, the former have flatter SEDs shortward of 60 micron. A mild anti-correlation is found between the [200 - 100] color and the "60 micron excess". We infer that this is due to the fact that galaxies with a strong starburst component, and thus a strong 60 micron flux, have a steeper far-infrared turnover. In non-Seyfert galaxies, increasing the luminosity corresponds to increasing the star formation rate, that enhances the 25 and 60 micron emission. This shifts the peak emission from around 150 micron in the most quiescent spirals to shorter than 60 micron in the strongest starburst galaxies.

We present the results of sub-mm, mm (850 um, 450 um and 1250 um) and radio (1.4 and 4.8 GHz) continuum observations of a sample of 27 K-selected Extremely Red Objects, or EROs, (14 of which form a complete sample with K < 20 and I-K > 5) aimed at detecting dusty starbursts, deriving the fraction of UltraLuminous Infrared Galaxies (ULIGs) in ERO samples, and constraining their redshifts using the radio-FIR correlation. One ERO was tentatively detected at 1250 um and two were detected at 1.4 GHz, one of which has a less secure identification as an ERO counterpart. Limits on their redshifts and their star forming properties are derived and discussed. We stacked the observations of the undetected objects at 850 um, 1250 um and 4.8 GHz in order to search for possible statistical emission from the ERO population as a whole, but no significant detections were derived either for the whole sample or as a function of the average NIR colours. These results strongly suggest that the dominant population of EROs with K < 20 is not comprised of ULIGs like HR 10, but is probably made of radio-quiet ellipticals and weaker starburst galaxies with L < 10^{12} L_sun and SFR < 100 M_sun/yr.

Increasing observational evidence supports a picture for a close link among AGN phenomenum, star-formation processes and galaxy formation. Since the physical phenomena related to the onset of both AGN and star formation are very likely characterized by a strong far infrared (FIR) emission it becomes mandatory to investigate the FIR energy domain. However, because of their faintness in the FIR energy range, very little is known about the FIR properties of type 1 AGN, and type 2 AGN, which are supposed to be missed in classical searches, still lack definition and samples. Our aim is to investigate the FIR properties of optically selected type 1 AGN, compare them to Ultraluminous Infrared Galaxies (ULIRGs) and derive some general characteristics of the population. Expectations for future surveys are also presented.

We present far-infrared and millimetric observations of a sample of 3C objects. Millimetre data were taken at 1.25 mm with the IRAM 30m antenna feeding the MPIfR bolometer array. Mid-infrared (MIR) and far-infrared (FIR) photometry were carried out with the ISOCAM and ISOPHOT cameras on the ISO Satellite. Additional FIR IRAS observations are also included. We present the entire Spectral Energy Distributions (SEDs) from the UV to radio and discuss the emitting mechanisms. Two composite spectra, one for the radio galaxies and one for the radio quasars, are built in the object rest frame. While the SEDs of the two classes differ, they are indistinguishable in the MIR and FIR range where they probably arise as thermal emission from a dusty torus and a larger-scale (cooler) dust distribution in the host galaxy.

The IRAM Plateau de Bure mm interferometer and deep K-band imaging have been used to identify the brightest submm source detected in the Lockman field of the UK 8mJy SCUBA survey. The near infrared counterpart is an extended (20-30kpc), clumpy, and extremely red object. The spectral energy distribution suggests it to be a dusty star forming object at a redshift of about 3 (2-4). Its star formation rate and near-infrared properties are consistent with Lockman850.1 being a massive elliptical in formation.

In this paper we will consider the problem of the numerical simulation of non-Gaussian, scalar random fields with a prescribed correlation structure provided either by a theoretical model or computed on a set of observational data. Although, the numerical generation of a generic, non-Gaussian random field is a trivial operation, the task becomes tough when constraining the field with a prefixed correlation structure. At this regards, three numerical methods, useful for astronomical applications, are presented. The limits and capabilities of each method are discussed and the pseudo-codes describing the numerical implementation are provided for two of them.

We report the detection of a diffuse radio halo source in the hottest known cluster of galaxies 1E0657-56 (RXJ0658-5557). The radio halo has a morphology similar to the X-ray emission from the hot intracluster medium. The presence of a luminous radio halo in such a hot cluster is further evidence for a steep correlation between the radio halo power and the X-ray temperature. We favour models for the origin of radio halo sources involving a direct connection between the X-ray emitting thermal particles and the radio emitting relativistic particles.

CO J = 5 - 4 and J = 2 - 1 emission lines were detected towards the extremely red galaxy (ERG) HR10 (J164502+4626.4) at z=1.44. The CO intensities imply a molecular gas mass M(H2) of 1.6 x 10^{11} h^{-2}_50 solar masses, and, combined with the intensity of the dust continuum, a gas-to-dust mass ratio around 200-400 (assuming galactic values for the conversion factors). The peak of the CO lines are at the same redshift as the [O II]3727 line, but blue- shifted by 430 km s^{-1} from the Halpha line. These CO detections confirm the previous results that HR10 is a highly obscured object with a large thermal far-infrared luminosity and a high star-formation rate. The overall properties of HR10 (CO detection, L_FIR to L'_CO ratio, and FIR to radio flux ratio) clearly favour the hypothesis that its extreme characteristics are related to star-formation processes rather than to a hidden AGN.

We present near-infrared VLT-UT1+ISAAC spectroscopy of a sample of 9 extremely red galaxies (ERGs) with R-K>5 and K<19.0. Neither strong emission lines (F_lim<1-5 10^{-16} erg s^{-1}cm^{-2}) nor continuum breaks are detected. From near-infrared spectrophotometry, complemented with broad-band optical and near-IR photometry, we estimate ``spectro-photometric'' redshifts to be in the range of 0.8<z_sphot<1.8. We derive upper limits on the star formation rates in range of SFR<6-30h_50^{-2} M_solar yr^{-1}. Two of the observed ERGs are dusty starburst candidates because they require strong dust reddening to reproduce their global spectral energy distributions. The other ERGs are consistent with being dustless old passively evolved spheroidals at z>0.8. We discuss the general implications of our findings in relation with the problem of the formation of early type galaxies.

The IRAM interferometer has been used to detect the submm source HDF 850.1 found by Hughes et al. (1998) in the Hubble Deep Field. The flux density measured at 1.3mm is 2.2 mJy, in agreement with the flux density measured at the JCMT. The flux densities and upper limits measured at 3.4, 2.8, 1.3, 0.85, and 0.45 mm show that the emission is from dust. We suggest that the 1.3mm dust source is associated with the optical arc-like feature, 3-593.0, that has a photometric redshift of about 1.7. If HDF 850.1 is at this redshift and unlensed, its spectral energy distribution, combined with that of 3-593.0, matches closely that of the ultraluminous galaxy VII Zw 31. Another possibility is that the dust source may be gravitationally lensed by the elliptical galaxy 3-586.0 at a redshift of 1.

Preliminary results of a project aiming at unveiling the nature of the extremely red galaxies (ERGs) (objects with colours R-K>6 and I-K>5) found in deep optical-NIR surveys are presented. Very little is known about these objects, the critical issue being whether they are old ellipticals at z>1 or distant star-forming galaxies strongly reddened by dust extinction. We expect to shed light onto the unknown nature of these galaxies by completing our three-step project: (1) the construction of two very deep optical/NIR surveys to select ERGs, (2) subsequent VLT/NIR spectroscopy; (3) observations in the sub-mm region with SCUBA at the JCMT and with MPIfRbolo at the IRAM 30m antenna.

We test emission models of circum-nuclear dust torii around quasars, at low and high redshifts, by using a large collection of photometric data for an unbiased sample of 120 optically-selected objects with millimetric and sub-millimetric fluxes, including new unpublished data. Under the assumption that the dust is heated by a point-like source with a power-law primary spectrum, as defined by the observed optical-UV continuum, we infer the basic model parameters, such as dust masses, temperature distributions and torus sizes, by numerically solving the radiative transfer equation in the dust distribution. In addition to the substantiated statistics, an essential improvement over previous analyses comes from the use of optical-UV data to constrain the primary illuminating continuum, which is needed to estimate dust temperatures and sizes. Dependences of the best-fit parameters on luminosity and redshift are studied and the contribution of dust in the host galaxy to the observed fluxes is briefly mentioned. This analysis constrains the properties of the enriched interstellar medium in the galaxies hosting the quasars. The dust abundance does not display appreciable trends as a function of redshift, from z ~ 1 to almost 5, and shows that dust and metals are at least as, and often more, abundant at these early epochs than they are in local galactic counterparts. This evolutionary pattern is remarkably at variance with respect to what is expected for disk galaxies, like the Milky Way, slowly building metals during the whole Hubble time. It rather points in favour of a much more active phase of star-formation at early epochs, probably related to the formation of massive spheroidal galaxies.

We report simultaneous observations at 1.2 and 2 mm, with a double channel photometer on the SEST Telescope, of the X-ray cluster RXJ0658-5557 in search for the Sunyaev-Zel'dovich (S-Z). The S-Z data were analyzed using the relativistically correct expression for the Comptonization parameter and we find from the detected decrement (2.60 +/- 0.79) ~ 10^{-4}, which is consistent with that computed using the X-ray (ROSAT and ASCA) observations. The uncertainty includes contributions due to statistical uncertainty in the detection and systematics and calibration. The 1.2 {mm} channel data alone gives rise to a larger Comptonization parameter and this result is discussed in terms of contamination from foreground sources and/or dust in the cluster or from a possible systematic effect. We then make use of the combined analysis of the ROSAT and ASCA X-ray satellite observations to determine an isothermal model for the S-Z surface brightness. Since the cluster is asymmetrical and probably in a merging process, models are only approximate. The associated uncertainty can, however, be estimated by exploring a set of alternative models. We then find as the global uncertainty on the Comptonization parameter a factor of 1.3. Combining the S-Z and the X-ray measurements, we determine a value for the Hubble constant. The 2 mm data are consistent with H_0(q_0 = 1/2)= 53^{+38}_{-28} km/s Mpc^{-1}, where the uncertainty is dominated by the uncertainty in models of the X-ray plasma halo.

The nature of the extremely red galaxies is poorly known, the critical issue being whether they are old ellipticals at $z>1.5$ or distant star-forming galaxies strongly reddened by dust extinction. Here we review the little we know about them, and we present recent and preliminary results from our ongoing submillimeter -- millimeter observations.

Near-infrared surveys have revealed a substantial population of enigmatic faint galaxies with extremely red optical-to-near-infrared colours and with a sky surface density comparable to that of faint quasars. There are two scenarios for these extreme colours: (i) these distant galaxies have formed virtually all their stars at very high redshifts and, due to the absence of recently formed stars, the colours are extremely red and (ii) these distant galaxies contain large amounts of dust, severely reddening the rest-frame UV--optical spectrum. HR10 ($z = 1.44$) is considered the archetype of the extremely red galaxies. Here we report the detection of the continuum emission from HR10 at 850$\mu$m and at 1250$\mu$m, demonstrating that HR10 is a very dusty galaxy undergoing a major episode of star formation. Our result provides a clear example of a high-redshift galaxy where the star formation rate inferred from the ultraviolet luminosity would be underestimated by a factor up to 1000, and shows that great caution should be used to infer the global star formation history of the Universe from optical observations only.

Recent results have questioned the description of the QSO luminosity function in terms of a pure luminosity evolution and call for a luminosity dependent luminosity evolution. Measurements of the QSO clustering amplitude and evolution allow further distinguishing among the various physical scenarios proposed to interpret the QSO phenomenon. The general properties of the QSO population would arise naturally if quasars are short-lived events connected to a characteristic halo mass $\sim 5 \cdot 10^{12}$ M$_\odot$. This is the typical mass of groups of galaxies in which the interactions triggering the QSO activity preferentially take place.

The evolution of QSO clustering is investigated with a new sample of 388 QSOs with 0.3<z<2.2, B<20.5 and M_B<-23. Evidence is found for an increase of the clustering amplitude with increasing redshift. These measurements allow to further distinguish among the various physical scenarios proposed to interpret the QSO phenomenon. A single population model is inconsistent with the observations. The general properties of the QSO population would arise naturally if quasars are short-lived events connected to a characteristic halo mass ~5x10^{12} M_sun.

The evolution of quasar clustering is investigated with a new sample of 388 quasars with 0.3<z<=2.2, B<=20.5 and Mb<-23, selected over an area of 24.6 sq. deg. in the South Galactic Pole. Assuming a two-point correlation function of the form xi(r) = (r/r_o)^-1.8, we detect clustering with r_0=6.2 +/- 1.6 h^-1 comoving Mpc at an average redshift of z=1.3. We find a 2 sigma significant increase of the quasar clustering between z=0.95 and z=1.8, independent of the quasar absolute magnitude and inconsistent with recent evidence on the evolution of galaxy clustering. If other quasar samples are added (resulting in a total data-set of 737 quasars) the increase of the quasar clustering is still favoured although it becomes less significant. We find epsilon=-2.5. Evolutionary parameters epsilon>0.0 are excluded at a 0.3% probability level, to be compared with epsilon=0.8 found for galaxies. The observed clustering properties appear qualitatively consistent with a scenario of Omega=1 CDM in which a) the difference between the quasar and the galaxy clustering can be explained as a difference in the effective bias and redshift distributions, and b) the quasars, with a lifetime of t~10^8 yr, sparsely sample halos of mass greater than M_min~10^12-10^13 h^-1 M_sun. We discuss also the possibility that the observed change in the quasar clustering is due to an increase in the fraction of early-type galaxies as quasar hosts at high z.

We exploit observations at 1.25 mm with the ESO-SEST telescope of a southern galaxy sample, selected from the IRAS PSC and complete to $S_{60}=2 Jy$, to derive the FIR and mm luminosity functions and the conditional probability distributions of FIR and mm luminosity of galaxies. The reliability of these estimates is ensured by the good observed correlation of the far-infrared and mm emissions. This detailed knowledge of the millimetric properties of galaxies is used to simulate the extragalactic sub-mm sky (background intensity, small-scale anisotropy signals, and discrete source statistics) -- that will be investigated soon by a variety of ground-based and space observatories. We find, in particular, that the recent tentative detection of a sub-mm background would require, if confirmed, strong evolution with cosmic time of the galaxy long-wavelength emissivity. We finally emphasize the difficulty to test such evolution through observations from currently available millimetric sites on ground.

We use photometric data, from the optical to the mm for a large sample of optically selected radio-quiet quasars, at low and high redshifts, to test emission models from circum-nuclear dusty torii around them. Model parameters, such as dust mass, temperature distributions and torus sizes are inferred, under the assumption that the dust is heated only by the central nuclear source. Dependences of best-fit parameters on luminosity and redshift are studied and the contribution of dust in the host galaxy to the observed fluxes is briefly mentioned. Only tentative conclusions can be drawn from this yet large dataset and detailed modelling, because crucial spectral information is missing at wavelengths where the QSO dust emission is more prominent, i.e. at $\lambda \simeq 50$ to 600 $\mu m$. Poor information in this wavelength range affects our understanding not only of the high-redshift objects, but also of the low-redshift sample. We then expect fundamental contributions by the FIRST mission in this field, which is obviously related to the more general one of galaxy formation.

The long-term variability of a sample of 149 optically selected QSOs in the field of the Selected Area 94 has been studied in the R-band. The relations between variability and luminosity and between variability and redshift have been investigated by means of ``robust'' statistical estimators, allowing to disentangle the effects of the measurement errors. The results are compared with the corresponding properties of the variability in the B-band for the same sample. An anti-correlation between the R-band variability and the intrinsic luminosity is found, analogously to what is observed in the B-band. The amplitude of the R-band variability turns out to be smaller (of a factor $1.13 \pm 0.05$) than the B-band variability. The implications in terms of the black-hole, starburst and microlensing models are discussed.

An observing campaign was devoted to the search for the Sunyaev-Zeldovich (S-Z) effect towards X-ray ROSAT Clusters in the millimetric spectral domain. A double channel (1.2 and 2 {\it mm}) photometer was installed at the focus of the 15m Swedish ESO Submillimeter Telescope (SEST) in Chile in september 1994 and 1995 and observations of the targets S1077, A2744, S295 and RXJ0658-5557 were gathered. Detections were found for A2744 at 1 {\it mm} and in both channels (at 1.2 and 2 {\it mm}) towards RXJ0658-5557. For the first time there is evidence for the S-Z enhancement and both the latter and the decrement were detected on the same source. We discuss astrophysical and systematic effects which could give origin to these signals.

We have studied the FIR/{\it mm} spectrum of IR galaxies by combining IRAS photometry with new {\it mm} data on a complete southern IRAS galaxy sample. The observed spectra and a dust model emphasize a dicothomy in the galaxy population: half of the objects with a lot of warm dust are characterized by higher values of the bolometric (UV-FIR) luminosity, of the dust-to-gas mass ratio, of the dust optical depths and extinction, while those dominated by cold ({\it cirrus}) dust show opposite trends. From these data we derive the {\it mm} luminosity function of galaxies and estimate their contribution to the sub-{\it mm} background (BKG).

We present 1.25 mm continuum data for a southern galaxy sample selected from the IRAS PSC and complete to S_60=2 Jy. Two thirds of the galaxies have been detected and significant limits on the remaining objects have been set. We find, on a statistical basis, indications that the dust emission in these galaxies is somewhat more centrally concentrated than that of the optical light, possibly tracing a higher metal content in the inner galactic regions. This result also allows to estimate the aperture corrections to the millimetric data. The latter, together with IRAS photometric data, have been used to compare the broad-band FIR/mm spectra with a simple dust model. According to their far-IR/mm spectrum, the sample galaxies show a dichothomy: almost half of the objects, those displaying bright 25-60um fluxes ascribed to warm dust residing in 'starburst' regions, are characterized by higher values of the bolometric (optical + FIR) luminosity, of the dust-to-gas mass ratio, of the dust optical depths and of the overall extinction. A complementary class of objects dominated by cold dust ('cirrus') shows opposite trends. Because of the favourable observational setup, selection wavelength and completeness, we believe these data provide an exhaustive and unbiased view of dust properties in spiral galaxies.

We report on observations of the Sunyaev-Zeldovich effect towards X-ray ROSAT clusters taken with a double channel (1.2 and 2 mm) photometer installed at the focus of the 15m SEST antenna in Chile. This paper describes the first results obtained for the high-z clusters S1077, A2744 and S295. Marginal detections were found for A2744 and at 1 mm for S1077. We discuss these data in terms of contamination of sources along the line of sight and give a constraint on the amplitude of the kinematic effect.

The long-term variability of a sample of 486 optically selected QSOs in the fields of the SA94, SA57 and SGP has been studied. The relations of variability with luminosity and redshift have been investigated by means of statistical estimators that are ``robust'' and allow at the same time to eliminate the influence of the measurement errors. The analysis in the QSOs rest frame of both the ensemble structure function ($SF$) and individual variability indices show: 1) a negative correlation between variability and luminosity; 2) a significant positive correlation between variability and $z$; 3) such correlations may be equally well parameterized either with a model in which the timescale of the variability is fixed for all the QSOs ($\tau \sim 2.4$ yr), while the amplitude linearly increases with the absolute magnitude and $z$, or with a model in which the timescale of the variability linearly depends on the absolute magnitude and the amplitude is only a function of $z$. The same analysis carried out in the observer's frame shows: 1) a negative correlation between variability and luminosity; 2) the timescale of variability does not depend significantly either on the absolute magnitude or on $z$; 3) the ensemble structure function is well represented by a parameterization in which, with a fixed timescale of about 5.5 yr, the amplitude linearly increases with the absolute magnitude; 4) although the general behaviour of the $SF$ does not show a systematic variation of the timescale and/or amplitude with $z$, if we examine the average variability index for objects with $-25>M_B>-27$, we find that below $z=1$ QSOs are significantly less variable than at higher redshift. Implications in terms of the black-hole, starburst and microlensing models are briefly discussed.

This is the first paper in a series aimed at defining a statistically significant sample of QSOs in the range $ 15 < B < 18.75$ and $ 0.3 < z < 2.2$. The selection is carried out using direct plates obtained at the ESO and UK Schmidt Telescopes, scanned with the COSMOS facility and searched for objects with an ultraviolet excess. Follow-up spectroscopy, carried out at ESO La Silla, is used to classify each candidate. In this initial paper, we describe the scientific objectives of the survey; the selection and observing techniques used. We present the first sample of 285 QSOs ($M_B < -23$) in a 153 deg$^2$ area, covered by the six ``deep'' fields, intended to obtain significant statistics down $B \simeq 18.75$ with unprecedented photometric accuracy. From this database, QSO counts are determined in the magnitude range $ 17 < B < 18.75$.

We have investigated the spatial distribution of quasars and its relationship with redshift by using the two-point correlation function, the variance of cell counts and the conditional density as a function of redshift. By comparing our results with those found by other authors we conclude that quasars have a correlation function with amplitude similar to or greater than that of galaxies, but the issue of its evolution still needs to be settled. We discuss how this can be achieved by adding new observations and comparing real data with survey simulations.

Using IRAM 30m, SEST 15m, and Nancay radiotelescopes, we have gathered the 1mm continuum emission, spectra of the J=1--0 CO line and of the HI line at 21cm for two samples of IRAS galaxies. Using these data, we have estimated the HI masses, the H$_2$ masses and the dust masses. We find that for far-infrared selected galaxies: (1) the median value of M$_{H2}$/M$_{HI}$ is 0.5, meaning that the atomic phase dominates in these galaxies. The fraction of gas in molecular form increases with increasing FIR luminosity but does not show any obvious trend with other galaxy properties. (2) the $H2$ surface density (s.d.) is better correlated with the cold dust s.d. than the HI s.d., but the correlation of HI with dust is not negligible. Globally in these galaxies the cold dust emission is associated with both the molecular and atomic phases. (3) the FIR surface brightness increases as the third power of the S60um/S100um ratio, it shows a tight correlation with both the H$_2$ and dust surface densities and a weaker one with the HI surface density: a large part of the far-infrared emission of these galaxies originates in the molecular medium. (4) the gas-to-dust ratio ranges between 100 and 1000 and its average value is 230, close to the Galactic value. There is indeed a clear trend: this ratio decreases as the FIR surface density increases, a result which can be explained in the framework of an enhancement of metallicity in galaxy discs having a higher star formation rate.

We report on observations performed at 1.25 $mm$ of a southern galaxy sample, selected from the IRAS PSC and complete to $S_{60}=2 Jy$. We detected 18 sources and set significant limits on 10 further objects. We use these data to discuss the spatial distribution of cold dust, the broad-band far-IR/mm spectra, the overall amount of dust, the gas-to-dust mass ratio, the dust optical depth, and the overall extinction, for such a representative galaxy sample. These results are also supported by a successful comparison with values inferred for the Galaxy. Because of the favourable observational setup, selection wavelength and completeness, we believe these data provide an unbiased view of dust properties in spiral galaxies.

{}From three quasar samples with a total of 1038 objects in the redshift range $1.0 \div 2.2$ we measure the variance $\sigma^2$ of counts in cells of volume $V_u$. By a maximum likelihood analysis applied separately on these samples we obtain estimates of $\sigma^2(\ell)$, with $\ell \equiv V_u^{1/3}$. The analysis from a single catalog for $\ell = ~40~h^{-1}$ Mpc and from a suitable average over the three catalogs for $\ell = ~60,~80$ and $100~h^{-1}$ Mpc, gives $\sigma^2(\ell) = 0.46^{+0.27}_{-0.27}$, $0.18^{+0.14}_{-0.15}$, $0.05^{+0.14}_{-0.05}$ and $0.12^{+0.13}_{-0.12}$, respectively, where the $70\%$ confidence ranges account for both sampling errors and statistical fluctuations in the counts. This allows a comparison of QSO clustering on large scales with analogous data recently obtained both for optical and IRAS galaxies: QSOs seem to be more clustered than these galaxies by a biasing factor $b_{QSO}/b_{gal} \sim 1.4 - 2.3$.