M. Vaccari

Extragalactic optical-infrared background radiation, its time evolution and the cosmic photon-photon opacity

Context. The background radiation in the optical and the infrared cause energy loss in the propagation of high energy particles through space. In particular, TeV observations with Cherenkov telescopes of extragalactic sources are influenced by the opacity effects due to the interaction of the very high-energy source photons with the background light. Aims. With the aim of assessing with the best possible detail these opacity terms, we have modelled the extragalactic optical and infrared backgounds using available information on cosmic sources in the universe from far-UV to sub-millimeter wavelengths over a wide range of cosmic epochs. Methods. We have exploited the relevant cosmological survey data – including number counts, redshift distributions, luminosity functions – from ground-based observatories in the optical, near-IR, and sub-millimeter, as well as multi-wavelength information coming from the HST, ISO and Spitzer space telescopes. Additional constraints have been used from direct measurements or upper limits on the extragalactic backgrounds by dedicated missions (COBE). All data were fitted and interpolated with a multi-wavelength backward evolutionary model, allowing us to estimate the background photon density and its redshift evolution. From the redshift-dependent background spectrum, the photon-photon opacities for sources of high-energy emission at any redshifts were then computed. The same results can also be used to compute the optical depths for any kind of processes in the intergalactic space involving interactions with background photons (like scattering of cosmic-ray particles). Results. We have applied our photon-photon opacity estimates to the analysis of spectral data at TeV energies on a few BLAZARs of particular interest. The opacity-corrected TeV spectra are entirely consistent with standard photon-generation processes and show photon indices steeper than Γintrinsic = 1.6. Contrary to some previous claims, but in agreement with other reports, we find no evidence for any truly diffuse background components in addition to those from resolved sources. We have tested in particular the effects of a photon background originating at very high redshifts, as would be the emissions by a primeval population of Population III stars around z ∼ 10. We could not identify any opacity features in our studied BLAZAR spectra consistent with such an emission and place a stringent limit on such a diffuse photon intensity of ∼ 6n W/m 2 /sr between 1 and 4 μm. Conclusions. TeV observations of BLAZARs are consistent with background radiation contributed by resolved galaxies in the optical and IR, and exclude prominent additional components from very high-z unresolved sources.

The Herschel ATLAS

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 570 deg2 of the extragalactic sky, 4 times larger than all the other Herschel extragalactic surveys combined, in five far-infrared and submillimeter bands. We describe the survey, the complementary multiwavelength data sets that will be combined with the Herschel data, and the six major science programs we are undertaking. Using new models based on a previous submillimeter survey of galaxies, we present predictions of the properties of the ATLAS sources in other wave bands.

The SCUBA Half-Degree Extragalactic Survey - II. Submillimetre maps, catalogue and number counts

We present maps, source catalogue and number counts of the largest, most complete and unbiased extragalactic submillimetre survey: the 850-μm SCUBA Half-Degree Extragalactic Survey (SHADES). Using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), SHADES mapped two separate regions of sky: the Subaru/XMM–Newton Deep Field (SXDF) and the Lockman Hole East (LH). Encompassing 93 per cent of the overall acquired data (i.e. data taken up to 2004 February 1), these SCUBA maps cover 720 arcmin2 with a rms noise level of about 2 mJy and have uncovered >100 submillimetre galaxies. In order to ensure the utmost robustness of the resulting source catalogue, data reduction was independently carried out by four subgroups within the SHADES team, providing an unprecedented degree of reliability with respect to other SCUBA catalogues available from the literature. Individual source lists from the four groups were combined to produce a robust 120-object SHADES catalogue; an invaluable resource for follow-up campaigns aiming to study the properties of a complete and consistent sample of submillimetre galaxies. For the first time, we present deboosted flux densities for each submillimetre galaxy found in a large survey. Extensive simulations and tests were performed separately by each group in order to confirm the robustness of the source candidates and to evaluate the effects of false detections, completeness and flux density boosting. Corrections for these effects were then applied to the data to derive the submillimetre galaxy source counts. SHADES has a high enough number of detected sources that meaningful differential counts can be estimated, unlike most submillimetre surveys which have to consider integral counts. We present differential and integral source number counts and find that the differential counts are better fit with a broken power law or a Schechter function than with a single power law; the SHADES data alone significantly show that a break is required at several mJy, although the precise position of the break is not well constrained. We also find that a 850-μm survey complete down to 2 mJy would resolve 20–30 per cent of the far-infrared background into point sources.

The SCUBA HAlf Degree Extragalactic Survey – III. Identification of radio and mid-infrared counterparts to submillimetre galaxies

Determining an accurate position for a submillimetre (submm) galaxy (SMG) is the crucial step that enables us to move from the basic properties of an SMG sample – source counts and 2D clustering – to an assessment of their detailed, multiwavelength properties, their contribution to the history of cosmic star formation and their links with present-day galaxy populations. In this paper, we identify robust radio and/or infrared (IR) counterparts, and hence accurate positions, for over two-thirds of the SCUBA HAlf-Degree Extragalactic Survey (SHADES) Source Catalogue, presenting optical, 24-μm and radio images of each SMG. Observed trends in identification rate have given no strong rationale for pruning the sample. Uncertainties in submm position are found to be consistent with theoretical expectations, with no evidence for significant additional sources of error. Employing the submm/radio redshift indicator, via a parametrization appropriate for radio-identified SMGs with spectroscopic redshifts, yields a median redshift of 2.8 for the radio-identified subset of SHADES, somewhat higher than the median spectroscopic redshift. We present a diagnostic colour–colour plot, exploiting Spitzer photometry, in which we identify regions commensurate with SMGs at very high redshift. Finally, we find that significantly more SMGs have multiple robust counterparts than would be expected by chance, indicative of physical associations. These multiple systems are most common amongst the brightest SMGs and are typically separated by 2–6 arcsec, ~15–50/ sin i kpc at z∼ 2, consistent with early bursts seen in merger simulations.

HerMES: The SPIRE confusion limit

We report on the sensitivity of SPIRE photometers on the Herschel Space Observatory. Specifically, we measure the confusion noise from observations taken during the Science Demonstration Phase of the Herschel Multi-tiered Extragalactic Survey. Confusion noise is defined to be the spatial variation of the sky intensity in the limit of infinite integration time, and is found to be consistent among the different fields in our survey at the level of 5.8, 6.3 and 6.8 mJy/beam at 250, 350 and 500 microns, respectively. These results, together with the measured instrument noise, may be used to estimate the integration time required for confusion-limited maps, and provide a noise estimate for maps obtained by SPIRE.

HerMES : SPIRE galaxy number counts at 250, 350, and 500 μm

Emission at far-infrared wavelengths makes up a significant fraction of the total light detected from galaxies over the age of Universe. Herschel provides an opportunity for studying galaxies at the peak wavelength of their emission. Our aim is to provide a benchmark for models of galaxy population evolution and to test pre-existing models of galaxies. With the Herschel Multi-tiered Extra-galactic survey, HerMES, we have observed a number of fields of different areas and sensitivity using the SPIRE instrument on Herschel. We have determined the number counts of galaxies down to ~20 mJy. Our constraints from directly counting galaxies are consistent with, though more precise than, estimates from the BLAST fluctuation analysis. We have found a steep rise in the Euclidean normalised counts <100 mJy. We have directly resolved ~15% of the infrared extra-galactic background at the wavelength near where it peaks.

The Herschel Reference Survey

The Herschel Reference Survey is a Herschel guaranteed time key project and will be a benchmark study of dust in the nearby universe. The survey will complement a number of other Herschel key projects including large cosmological surveys that trace dust in the distant universe. We will use Herschel to produce images of a statistically-complete sample of 323 galaxies at 250, 350, and 500 μm. The sample is volume-limited, containing sources with distances between 15 and 25 Mpc and flux limits in the K band to minimize the selection effects associated with dust and with young high-mass stars and to introduce a selection in stellar mass. The sample spans the whole range of morphological types (ellipticals to late-type spirals) and environments (from the field to the center of the Virgo Cluster) and as such will be useful for other purposes than our own. We plan to use the survey to investigate (i) the dust content of galaxies as a function of Hubble type, stellar mass, and environment; (ii) the connection between the dust content and composition and the other phases of the interstellar medium; and (iii) the origin and evolution of dust in galaxies. In this article, we describe the goals of the survey, the details of the sample and some of the auxiliary observing programs that we have started to collect complementary data. We also use the available multifrequency data to carry out an analysis of the statistical properties of the sample.

Photometric redshifts in the SWIRE Survey

We present the SWIRE Photometric Redshift Catalogue 1 025 119 redshifts of unprecedented reliability and of accuracy comparable with or better than previous work. Our methodology is based on fixed galaxy and quasi-stellar object templates applied to data at 0.36–4.5 μm, and on a set of four infrared emission templates fitted to infrared excess data at 3.6–170 μm. The galaxy templates are initially empirical, but are given greater physical validity by fitting star formation histories to them, which also allows us to estimate stellar masses. The code involves two passes through the data, to try to optimize recognition of active galactic nucleus (AGN) dust tori. A few carefully justified priors are used and are the key to supression of outliers. Extinction, A_V , is allowed as a free parameter. The full reduced χ^2_ν (z) distribution is given for each source, so the full error distribution can be used, and aliases investigated. We use a set of 5982 spectroscopic redshifts, taken from the literature and from our own spectroscopic surveys, to analyse the performance of our method as a function of the number of photometric bands used in the solution and the reduced χ^2_ν . For seven photometric bands (5 optical + 3.6, 4.5 μm), the rms value of (z_(phot)−z_(spec)/(1 +z_(spec) is 3.5 per cent, and the percentage of catastrophic outliers [defined as >15 per cent error in (1 +z)], is ∼1 per cent. These rms values are comparable with the best achieved in other studies, and the outlier fraction is significantly better. The inclusion of the 3.6- and 4.5-μm IRAC bands is crucial in supression of outliers. We discuss the redshift distributions at 3.6 and 24 μm. In individual fields, structure in the redshift distribution corresponds to clusters which can be seen in the spectroscopic redshift distribution, so the photometric redshifts are a powerful tool for large-scale structure studies. 10 per cent of sources in the SWIRE photometric redshift catalogue have z > 2, and 4 per cent have z > 3, so this catalogue is a huge resource for high-redshift galaxies. A key parameter for understanding the evolutionary status of infrared galaxies is L_(ir)/L_(opt) . For cirrus galaxies this is a measure of the mean extinction in the interstellar medium of the galaxy. There is a population of ultraluminous galaxies with cool dust and we have shown SEDs for some of the reliable examples. For starbursts, we estimate the specific star formation rate, φ_*/M_* . Although the very highest values of this ratio tend to be associated with Arp220 starbursts, by no means all ultraluminous galaxies are. We discuss an interesting population of galaxies with elliptical-like spectral energy distributions in the optical and luminous starbursts in the infrared. For dust tori around type 1 AGN, L_(tor)/L_(opt) is a measure of the torus covering factor and we deduce a mean covering factor of 40 per cent. Our infrared templates also allow us to estimate dust masses for all galaxies with an infrared excess.

Mid- and far-infrared luminosity functions and galaxy evolution from multiwavelength Spitzer observations up to z ~ 2.5

Context. Studies of the infrared (IR) emission of cosmic sources have proven essential to constraining the evolutionary history of cosmic star formation and the gravitational accretion of nuclear black holes, because many of these events occur inside heavily dust-extinguished environments. Aims. The Spitzer Space Telescope has provided a large amount of data to constrain the nature and cosmological evolution of infrared source populations. In the present paper we exploit a large homogeneous dataset to derive a self-consistent picture of IR emission based on the time-dependent λ_(eff) = 24, 15, 12, and 8 μm monochromatic and bolometric IR luminosity functions (LF) over the full 0 < z < 2.5 redshift range. Methods. Our present analysis is based on a combination of data from deep Spitzer surveys of the VIMOS VLT Deep Survey (VVDS-SWIRE) and GOODS fields. To our limiting flux of S_(24) = 400 μJy, our sample derived from VVDS-SWIRE includes 1494 sources, and 666 and 904 sources brighter than S_(24) = 80 μJy are catalogued in GOODS-S and GOODS-N, respectively, for a total area of ~0.9 square degrees. Apart from a few galaxies, we obtain reliable optical identifications and redshifts for all these sources, providing a rich and robust dataset for our luminosity function determination. The final combined reliable sample includes 3029 sources, the fraction with photometric redshifts being 72% over all redshifts and almost all galaxies at z > 1.5. Based on the multiwavelength information available in these areas, we constrain the LFs at 8, 12, 15, and 24 μm. We also infer the total IR luminosities from our best-fit model of the observed SEDs of each source, and use this to derive the bolometric (8–1000 μm) LF and comoving volume emissivity to z ~ 2.5. Results. In the redshift interval 0 1. The mean redshift of the peak in the source number density shifts with luminosity: the brightest IR galaxies appear to form stars at earlier cosmic times (z > 1.5), while star formation in the less luminous galaxies continues until more recent epochs (z ~ 1 for L_(IR) < 10^(11)_☉), in overall agreement with similar analyses in the literature. Conclusions. Our results are indicative of a rapid increase in the galaxy IR comoving volume emissivity up to z ~ 1 and a constant average emissivity at z > 1. We also appear to measure a difference in the evolutionary rate of the source number densities as a function of luminosity, which is consistent with the downsizing evolutionary patterns reported for other samples of cosmic sources.

HerMES: Far infrared properties of known AGN in the HerMES fields

Nuclear and starburst activity are known to often occur concomitantly. Herschel-SPIRE provides sampling of the FIR SEDs of type 1 and type 2 AGN, allowing for the separation between the hot dust (torus) and cold dust (starburst) emission. We study large samples of spectroscopically confirmed type 1 and type 2 AGN lying within the Herschel Multi-tiered Extragalactic Survey (HerMES) fields observed during the science demonstration phase, aiming to understand their FIR colour distributions and constrain their starburst contributions. We find that one third of the spectroscopically confirmed AGN in the HerMES fields have 5-sigma detections at 250um, in agreement with previous (sub)mm AGN studies. Their combined Spitzer-MIPS and Herschel-SPIRE colours - specifically S(250)/S(70) vs. S(70)/S(24) - quite clearly separate them from the non-AGN, star-forming galaxy population, as their 24-um flux is dominated by the hot torus emission. However, their SPIRE colours alone do not differ from those of non-AGN galaxies. SED fitting shows that all those AGN need a starburst component to fully account for their FIR emission. For objects at z > 2, we find a correlation between the infrared luminosity attributed to the starburst component, L(SB), and the AGN accretion luminosity, L(acc), with L(SB) propto L(acc)^0.35. Type 2 AGN detected at 250um show on average higher L(SB) than type 1 objects but their number is still too low to establish whether this trend indicates stronger star-formation activity.