A&A 438, 475-480 (2005)
DOI: 10.1051/0004-6361:20042158
L. Toffolatti1 - M. Negrello2 - J. González-Nuevo1 - G. De Zotti2,3 - L. Silva4 - G. L. Granato2,3 - F. Argüeso5
1 - Departamento de Física, Universidad de Oviedo, c.
Calvo Sotelo s/n, 33007 Oviedo, Spain
2 -
International School for Advanced Studies,
SISSA/ISAS, Via Beirut 2-4, 34014 Trieste, Italy
3 -
INAF - Osservatorio Astronomico di Padova, Vicolo
dell'Osservatorio 5, 35122 Padova, Italy
4 -
INAF - Osservatorio Astronomico di Trieste, Via Tiepolo 11,
34131 Trieste, Italy
5 -
Departamento de Matemáticas, Universidad de Oviedo, c.
Calvo Sotelo s/n, 33007 Oviedo, Spain
Received 11 October 2004 / Accepted 21 April 2005
Abstract
The possible contributions of the various classes of
extragalactic sources (including, in addition to the canonical
radio sources, GHz Peaked Spectrum sources, advection-dominated
sources, starburst galaxies, high-redshift proto-spheroidal
galaxies) to the arcminute scale fluctuations measured by the CBI,
BIMA, and ACBAR experiments are discussed. At 30 GHz, fluctuations
due to radio sources undetected by ancillary low-frequency surveys
may be higher than estimated by the CBI and BIMA groups.
High-redshift dusty galaxies, whose fluctuations may be strongly
enhanced by the effect of clustering, could contribute to the BIMA
excess signal, and dominate at 150 GHz (the ACBAR frequency).
Moreover, in the present data situation, the dust emission of
these high-redshift sources set an unavoidable limit to the
detection of primordial CMB anisotropies at high multipoles, even
at frequencies as low as 30 GHz. It is concluded that the
possibility that the excess power at high multipoles is dominated
by unsubtracted extragalactic sources cannot be ruled out. On the
other hand, there is room for a contribution from the
Sunyaev-Zeldovich effect within clusters of galaxies, with a
density fluctuation amplitude parameter
consistent with
the values preferred by current data.
Key words: cosmic microwave background - galaxies: general - radio continuum: galaxies
In the past few years different experiments (BIMA: Dawson et al.
2002; CBI: Mason et al. 2003; Readhead et al. 2004; and ACBAR: Kuo
et al. 2004), aimed at measuring the anisotropies of the cosmic
microwave background (CMB) on arcmin angular scales, have detected
signals at multipoles
in excess of the expected
primordial CMB anisotropies. The origin of this excess signal, in
the range
K, is not well understood
yet, although several possibilities have been discussed in the
literature.
All experimental groups argue that it cannot be due to
point-source contamination. If so, the most likely candidate is
the thermal Sunyaev-Zeldovich (SZ) effect, which is expected to
dominate CMB anisotropies on angular scales of a few arcminutes
(Gnedin & Jaffe 2001). However, an interpretation on terms of SZ
effects from clusters of galaxies (Bond et al. 2005; Komatsu &
Seljak 2002) or to the inhomogeneous plasma distribution during
the formation of large scale structure (Zhang et al. 2002)
requires values of
(the rms density fluctuation on a
scale of 8 h-1 Mpc) significantly higher than indicated by
current data. SZ effects associated with the formation and the
early evolutionary phases of massive spheroidal galaxies could
account for the BIMA signal, although some parameters need to be
stretched to their boundary values (De Zotti et al. 2004).
Alternative interpretations advocate non-standard inflationary
models (Cooray & Melchiorri 2002; Griffiths et al. 2003).
In this paper we revisit the contributions of extragalactic point sources to the power spectrum on arcminute scales at the relevant frequencies, including the possible role of faint sources, with flux densities too weak to be filtered out, and the effect of clustering. The outline of the paper is as follows. In Sect. 2 we describe the different source populations which give the dominant contributions to number counts at cm and mm wavelengths. In Sect. 3 we present our estimates of arcminute-scale CMB anisotropies due to extragalactic sources, while in Sect. 4 we summarize our main conclusions.
A flat CDM cosmology with
has been used throughout the paper.
The estimated contributions of the various populations of extragalactic sources to the counts at 30 GHz (the frequency of BIMA and CBI experiments) and at 150 GHz (ACBAR experiment), obtained from the model of De Zotti et al. (2005), which updates the model by Toffolatti et al. (1998), are shown in Fig. 1. In addition to the canonical flat- and steep-spectrum radio sources, the model takes into account star-forming galaxies with their complex spectra including both radio (synchrotron plus free-free) and dust emission, and the source populations characterized by spectra peaking at high radio frequencies, such as extreme GHz Peaked Spectrum (GPS) sources and accretion flows on almost inactive supermassive black-holes in early type galaxies (ADAF/ADIOS sources).
Because of their inverted low-frequency spectra, GPS and ADAF/ADIOS sources are potentially worrisome. However, GPS sources are rare and the analysis made by De Zotti et al. (2000) implies that they likely have very flat counts and therefore are minor contributors to small scale fluctuations. Furthermore, the repeated multifrequency measurements by Tinti et al. (2005) have demonstrated that most GPS candidates identified with quasars in the sample of Dallacasa et al. (2000) are in fact flaring blazars, so that the surface densities of bona-fide GPS sources is probably substantially lower than estimated by De Zotti et al. (2000), a conclusion further supported by an examination, carried out by De Zotti et al. (2005), of GPS candidates in the WMAP sample (Bennett et al. 2003).
ADAF/ADIOS sources are far more numerous, but have a low radio power. The estimate by De Zotti et al. (2005) of their counts is well below that by Perna & Di Matteo (2000), whose results are probably affected by a numerical error, and implies that also these sources do not contribute significantly to the fluctuations measured by the BIMA and CBI experiments. On the other hand, Pierpaoli & Perna (2004; model A) pointed out that if the standard ADAF model (Narayan & Yi 1994) is used, these sources could make up to 40-50% of the BIMA and CBI excesses. We note, however, that the standard ADAF scenario faces a number of serious difficulties, some of which are summarized in Sect. 4.3 of De Zotti et al. (2005), suggesting that the radio emission is suppressed by massive outflows. It is therefore likely that the results of model A by Pierpaoli & Perna (2004) should be regarded as, probably generous, upper limits.
As for starburst galaxies, the slope of their differential counts
can exceed 3, if these objects have to account for the very steep
ISOCAM m counts below a few mJy, as implied by recent
analyses (Gruppioni et al. 2003; Franceschini et al. 2003; Pozzi
et al. 2004; Silva et al. 2005). In this case, their main
contribution to small scale fluctuations comes from weak sources,
at
Jy levels, far fainter than those removed from CBI and
BIMA maps. On the other hand, the counts of such sources are
tightly constrained by
Jy counts at 1.4 GHz (Richards 2000),
5 GHz (Fomalont et al. 1991), and 8.4 GHz (Fomalont et al. 2002).
Taking such constraints into account and applying an average
spectral index
,
appropriate for this class of
sources, we find that they can only provide a minor contribution
to the excess power detected by BIMA and CBI:
at
,
by using the nominal
detection limit,
Jy, for BIMA and
at
(
mJy) for CBI,
respectively. Their contribution is negligible at the ACBAR
frequency.
An additional contribution is expected from dusty proto-spheroidal
galaxies, which may account for galaxies selected by SCUBA and
MAMBO surveys (Granato et al. 2001, 2004), whose counts at
m and 1.2 mm appear to fall down very rapidly at flux
densities above several mJy (Scott et al. 2002; Borys et al. 2003;
Greve et al. 2004). The spectral energy distribution of nearby
dusty galaxies is dominated by synchrotron plus free-free emission
at
-3 mm (Bressan et al. 2002), while at shorter
wavelengths dust emission, rapidly raising with frequency (
), takes over. Since these sources are at typical
redshifts >2 (Chapman et al. 2003), dust emission can
significantly contribute to the counts even at 30 GHz. When dust
emission comes in, the counts, already steep because of the effect
of the strong cosmological evolution, are boosted by the large
negative K-correction.
![]() |
Figure 1:
Differential counts,
![]() |
Open with DEXTER |
The poor knowledge of the millimeter emission of these sources,
however, makes estimates of their contributions to the 30 GHz
counts quite uncertain. The two short-dashed lines in Fig. 1 show
the counts we obtain using the physical evolutionary model by
Granato et al. (2004) but with two choices for the spectral energy
distribution (SED). The lower (thicker) line refers to the SED
produced by the code GRASIL (originally described by Silva et al.
1998). An excess emission by a factor 2 at
mm was however detected in several Galactic clouds, combining
Archeops with WMAP and DIRBE data (Bernard et al. 2003; Dupac et al. 2003), and in NGC 1569 (Galliano et al. 2003). The origin of
the excess is still not understood. Possibilities discussed in the
literature are that the grain sizes or composition change in dense
environments or that there is an intrinsic dependence of the dust
emissivity index on temperature (Dupac et al. 2004). If the excess
is due to very cold grains (Reach et al. 1995; Galliano et al.
2003) it cannot be present in the high-z proto-spheroids. But if
it is a general property of the SED of dusty galaxies, the
predicted counts of dusty proto-spheroids are given by the upper
(thin) short-dashed curve.
The strategy of the CBI group (Mason et al. 2003; Readhead et al.
2004) to remove the point source contamination comprises pointed
31 GHz observations with the OVRO 40 m telescope of all NVSS
sources with 1.4 GHz flux density 6 mJy, and direct counts
at 31 GHz using the CBI deep and mosaic maps. Although the
threshold of OVRO observations is 6 mJy, the survey is
99% complete only at
mJy. The limiting flux
density ranges from 6 to 12 mJy in the deep CBI maps, and from 18
to 25 mJy in the mosaic maps. Subtraction of OVRO detected sources
removes two-thirds of the observed power level.
Furthermore, they have adopted the constraint matrix approach to
remove from their dataset all NVSS sources with flux densities
greater than 3.4 mJy at 1.4 GHz (Readhead et al. 2004), and have
estimated the contribution to fluctuations due to sources below
the NVSS cutoff using the observed OVRO-NVSS distribution of
spectral indices and adopting a rather shallow power-law slope for
the counts (
); for comparison, Richards
(2000) finds
for
mJy (see also Windhorst et al. 1993).
At the flux-density levels relevant for the CBI experiment, apart
from SZ effects (see Sect. 1), the dominant contribution to
fluctuations due to extragalactic sources is expected to come from
the classical steep- and flat-spectrum radio sources. However, an
accurate determination of the 30 GHz fluctuations due to sources
with
mJy is very difficult because of the
effect of sources with inverted spectra (in fact, Readhead et al.
2004 report the detection of a source, NGC 1068, with
mJy, not removed by the constraint matrix), and
of variability. The difficulty is illustrated by the results of
high-frequency surveys. Ricci et al. (2004) found that the 18 GHz
flux densities of extragalactic sources detected by the ATCA pilot
survey are not significantly correlated with the SUMSS flux
densities at 0.84 GHz (i.e. at a frequency not far from that of
the NVSS survey). Waldram et al. (2003) also reported a large
spread (about a factor of 10) of the 15 to 1.4 GHz flux density
ratios of sources detected in their 9C survey at 15 GHz, although
the flux densities at the two frequencies are correlated. They
also noted that pointed 15 GHz observations of the NVSS sources
with
mJy in the area covered by their survey
would have detected 434 sources above the 9C survey limit of 25 mJy
but would have missed 31 sources having
mJy
but
S1.4 < 25 mJy. The distribution of spectral indices has
a systematic drift towards flatter values with decreasing
low-frequency flux density down to
mJy
(Windhorst et al. 1993), so that the fraction of sources with
inverted spectra is expected to be higher at the fainter flux
density levels of interest here.
High frequency surveys emphasize flat-spectrum sources. The dominant flat-spectrum population are blazars, that are highly variable on timescales of years and whose variability amplitude increases with frequency (Impey & Neugebauer 1988; Ciaramella et al. 2004). The monitoring campaigns at 22, 37 and 87 GHz by the Metsähovi group (Teräsranta et al. 1998) have shown that intensity variations by factors of several are common at these frequencies, so that a substantial fraction of such sources may have had, at the moment of the CBI observations, 30 GHz fluxes higher than 3.4 mJy, even by a considerable factor. Variability can indeed account, to a large extent, for the lack of a correlation between the ATCA 18 GHz and the SUMSS 0.84 GHz flux densities (Ricci et al. 2004) and for the large spread of the 15 to 1.4 GHz flux density ratios (Waldram et al. 2003).
To appraise residual fluctuations at 30 GHz due to unsubtracted
sources we have adopted the analytical description of the counts
below a few mJy by Richards (2000;
with
and
)
and computed the Poisson fluctuations at 30 GHz of those
sources with
S1.4 < 3.4 mJy, assuming a Gaussian
distribution of spectral indices with mean
(Fomalont et al. 1991; Windhorst et al. 1993) and two values of
the dispersion (
or 0.4), based on the width of the
distribution of
of Waldram et al. (2003,
their Fig. 9). A comparison with the 30 GHz counts yielded by the
De Zotti et al. (2005) model, which takes into account the
available information from high-frequency surveys, shows that, in
the 30 GHz flux density range relevant to estimate fluctuations,
the counts extrapolated using the upper values of A and
are somewhat too high if
;
more consistent counts
(only slightly above the model predictions) are obtained with the
central values A=8.3,
.
To bound the plausible
range of residual fluctuations we have therefore considered the
cases A=8.3,
,
(upper dotted line in
the left-hand panel of Fig. 2) and A=7.9,
,
(lower dotted line).
![]() |
Figure 2:
Angular power spectrum (
![]() |
Open with DEXTER |
The additional contribution to CMB fluctuations given by
correlated positions in the sky of canonical steep- and
flat-spectrum radio sources has been recently analyzed by
González-Nuevo et al. (2005). Their outcomes indicate that the
extra power due to the clustering of radio sources cannot, by
itself, explain the excess signal detected by CBI and BIMA. Using
the
estimated by Blake & Wall (2002) from sources in
the NVSS survey down to
mJy - which can represent a
realistic approximation to the clustering properties of faint
undetected sources in the CBI fields - they found that clustered
radio sources at
mJy can give an extra
power
-4
K, which has to be summed up -
in quadrature - to the Poisson term,
-22
K.
The dominance of Poisson over clustering fluctuations even
at faint fluxes is due to the strong dilution of the clustering
signal of extragalactic radio sources by the broadness of their
luminosity function and of their redshift distribution (Dunlop &
Peacock 1990; Toffolatti et al. 1998; Negrello et al. 2004).
To remove the point source contamination, the BIMA group (Dawson
et al. 2002) have carried out a VLA survey at 4.8 GHz of their
fields. These observations reached a rms flux of Jy
for a 9' FWHM region with center
coinciding with the center of the corresponding BIMA field.
Sources with flux density
within 8' of the
pointing center were projected out. On the other hand, point
sources with flux densities
Jy, lying at an
angular distance
from the BIMA field center, cannot be
detected (and removed) by VLA observations if
,
where
is
the VLA response function, assumed Gaussian.
Therefore, we have estimated the fluctuations in the BIMA fields due to sources fainter than
Jy
,
where
is the angular distance from the pointing direction, by adopting the
number counts at 4.8 GHz of Fomalont et al. (1991),
arcmin-2, extrapolated to 30 GHz with a mean spectral index
(
),
appropriate for the relevant flux-density range (Fomalont et al. 1991; Windhorst et al. 1993). The corresponding
power spectrum is shown in Fig. 2 (right-hand panel), where the shaded area reflects the range of values
corresponding to the uncertainties in the counts of Fomalont et al. (1991) and in the dust emission spectrum of
high redshift proto-spheroidal galaxies (see below).
The fluctuations due to forming spheroidal galaxies, not represented in the 4.8 GHz counts, get comparable
contributions from both the Poisson and the clustering term, while the latter term turns out to be small,
compared to the former, for the other classes of sources relevant here. Adopting the standard expression for the
two-point correlation function,
with a constant comoving clustering length
r0=8.3 h-1 Mpc, h being the Hubble constant in units of
(see Negrello et al. 2004), we find, at
,
a Poisson contribution of
K and a clustering contribution of
K. Clearly these
contributions, to be summed in quadrature to the contribution discussed above, have a minor effect. If, however,
these sources show the mm excess mentioned in Sect. 2, their contribution to fluctuations would be approximately
doubled (see Fig. 2; right-hand panel), and, summed in quadrature to the above estimate of the contribution of
radio sources, could account for the reported excess signal.
The ACBAR measurements in the 150 GHz band reported by Kuo et al.
(2004) up to multipoles
are consistent with the
primordial CMB power spectrum predicted by standard cosmological
models. On the other hand, the measured power in the highest
multipole bin is also consistent with the excess detected by the
CBI experiment. Fluctuations due to extragalactic radio sources in
the ACBAR band are quite small (of a few
K at
). Quoting results by Blain et al. (1998), Kuo et al.
(2004) conclude that dusty proto-galaxies are also not expected to
contribute significantly to the observed signal. However, the
Blain et al. (1998) estimate actually refers to Poisson
fluctuations only, while it has been pointed out by several
authors (Scott & White 1999; Haiman & Knox 2000; Magliocchetti
et al. 2001; Perrotta et al. 2003; Negrello et al. 2004) that the
strong positional correlation of SCUBA galaxies, indicated by
observational data and theoretical arguments, implies that the
main contribution on ACBAR scales comes from source clustering. On
the other hand, the rest-frame mm-wave emission of these sources
is poorly known. If dusty proto-spheroidal galaxies have the
excess emission (compared to model expectations) discussed in
Sect. 3.2 (Fig. 3, right-hand panel), the power measured in the
highest multipole bin is easily explained in terms of fluctuations
due to undetected high-z galaxies. Anyway, dusty proto-spheroids
are always giving a relevant contribution to the measured power,
even without excess emission (Fig. 3, left-hand panel).
![]() |
Figure 3:
Angular power spectrum (
![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
Contamination from extragalactic point sources appears to be a likely candidate to account for a large fraction, perhaps for most of the excess power on arcminute scales detected by the CBI and BIMA experiments. The fluctuations due to radio sources undetected by the ancillary low-frequency surveys may in fact be higher than estimated by the CBI and BIMA groups. On the other hand, we argue that extreme GHz Peaked Spectrum sources and advection dominated sources, potentially worrisome because of their spectra peaking at high microwave/mm-wave wavelengths, should provide only a minor contribution to the CBI and BIMA signals.
Although the very steep m counts of starburst galaxies
below a few mJy imply a very strong cosmological evolution, the
radio surveys down to
Jy levels constrain the contribution of
their radio emission to fluctuations at 30 GHz to be
relatively small. On the other hand, the dust emission at
rest-frame mm wavelengths from star-forming galaxies at
high-redshifts, such as those detected by SCUBA and MAMBO surveys,
can be redshifted down to 30 GHz. Using the spectral energy
distributions given by the Granato et al. (2004) model, we find
that these sources yield fluctuations of a few to several
K
on arcminute scales. Their rest-frame spectral energy distribution
at mm wavelengths, however, is poorly known, and may well be
higher than implied by the adopted model.
Moreover, observational indications and the theoretical arguments
converge in suggesting that they are highly clustered (see
Negrello et al. 2004, and references therein), so that their
fluctuations may be strongly super-Poissonian for multipoles
(remember that the clustering-to-Poisson ratio
increases with the angular scale, i.e. with decreasing multipole
number, De Zotti et al. 1996). Clustering fluctuations of the
high-z galaxies detected by (sub)-mm SCUBA and MAMBO surveys may
indeed dominate the contamination by extragalactic sources of the
signal measured by the ACBAR experiment at 150 GHz.
Because the dust emission spectrum rises very steeply with
frequency, lower frequency surveys cannot be used to remove their
effect from 30 GHz maps. In the present data situation, they
therefore set an unavoidable limit to the determination of
the primordial CMB angular power spectrum at high multipoles, even
at frequencies as low as 30 GHz.
We stress that the present results are fully compatible with the
estimated contributions of Sunyaev-Zeldovich effects in clusters
of galaxies to the arcminute scale anisotropies. In fact, while an
interpretation of the full CBI and BIMA signals in terms of SZ
fluctuations would require a density fluctuation amplitude
(measured by the parameter )
at or above the limit
allowed by current data (Bond et al. 2005), our analysis leaves
room for an SZ contribution corresponding to the
values
favoured by analyses of CMB, cosmic shear, and large scale
structure data (Spergel et al. 2003; Pierpaoli et al. 2003; Van Waerbecke et al. 2005).
New interesting constraints on the CMB angular power spectrum up to
at 34 GHz should be provided
in the near future by the VSA experiment (see, e.g., Rebolo et al. 2004). The reduced noise level of the new
configuration and an effective cleaning of deep fields down to
5 mJy - by dedicated observations with
the Ryle Telescope at 15 GHz - will shed new light on the nature of the excess at high multipole and on the
point source populations mainly contributing to the number counts at
a few mJy. Moreover, Planck HFI data as well as the forthcoming surveys by the Herschel telescope - at frequencies where the
emission due to cold dust grains is the dominant one - shall be unique in determining much better the
cosmological evolution, the emission and the clustering properties of high-redshift dusty galaxies.
Acknowledgements
We are grateful to K. Dawson and to G. Holder for very stimulating comments and clarifications on the point source subtraction for the BIMA experiment, and to the referee whose comments greatly helped improving the paper. LT, JGN and FA thank the Spanish MEC (Ministerio de Educación y Ciencia) for partial financial support (project ESP2004-07067-C03-01). J.G.N. acknowledges a FPU fellowship and an "Ayuda'' for Short Research Periods of the Spanish Ministry of Education (MEC). J.G.N. also thanks the SISSA-ISAS, International School for Advanced Studies (Trieste, Italy), where his share of this work was completed, for the warm hospitality.