
We discuss the potential of a next generation spaceborne Cosmic Microwave
Background (CMB) experiment for studies of extragalactic sources. Our analysis
has particular bearing on the definition of the future space project, CORE,
that has been submitted in response to ESA's call for a Mediumsize mission
opportunity as the successor of the Planck satellite. Even though the effective
telescope size will be somewhat smaller than that of Planck, CORE will have a
considerably better angular resolution at its highest frequencies, since, in
contrast with Planck, it will be diffraction limited at all frequencies. The
improved resolution implies a considerable decrease of the source confusion,
i.e. substantially fainter detection limits. In particular, CORE will detect
thousands of strongly lensed highz galaxies distributed over the full sky. The
extreme brightness of these galaxies will make it possible to study them, via
followup observations, in extraordinary detail. Also, the CORE resolution
matches the typical sizes of highz galaxy protoclusters much better than the
Planck resolution, resulting in a much higher detection efficiency; these
objects will be caught in an evolutionary phase beyond the reach of surveys in
other wavebands. Furthermore, CORE will provide unique information on the
evolution of the star formation in virialized groups and clusters of galaxies
up to the highest possible redshifts. Finally, thanks to its very high
sensitivity, CORE will detect the polarized emission of thousands of radio
sources and, for the first time, of dusty galaxies, at mm and submm
wavelengths, respectively.

We examine the cosmological constraints that can be achieved with a galaxy
cluster survey with the future CORE space mission. Using realistic simulations
of the millimeter sky, produced with the latest version of the Planck Sky
Model, we characterize the CORE cluster catalogues as a function of the main
mission performance parameters. We pay particular attention to telescope size,
key to improved angular resolution, and discuss the comparison and the
complementarity of CORE with ambitious future groundbased CMB experiments that
could be deployed in the next decade. A possible CORE mission concept with a
150 cm diameter primary mirror can detect of the order of 50,000 clusters
through the thermal SunyaevZeldovich effect (SZE). The total yield increases
(decreases) by 25% when increasing (decreasing) the mirror diameter by 30 cm.
The 150 cm telescope configuration will detect the most massive clusters
($>10^{14}\, M_\odot$) at redshift $z>1.5$ over the whole sky, although the
exact number above this redshift is tied to the uncertain evolution of the
cluster SZE fluxmass relation; assuming selfsimilar evolution, CORE will
detect $\sim 500$ clusters at redshift $z>1.5$. This changes to 800 (200) when
increasing (decreasing) the mirror size by 30 cm. CORE will be able to measure
individual cluster halo masses through lensing of the cosmic microwave
background anisotropies with a 1$\sigma$ sensitivity of $4\times10^{14}
M_\odot$, for a 120 cm aperture telescope, and $10^{14} M_\odot$ for a 180 cm
one. [abridged]

We present 279 galaxy cluster candidates at $z > 1.3$ selected from the 94
deg$^{2}$ Spitzer South Pole Telescope Deep Field (SSDF) survey. We use a
simple algorithm to select candidate highredshift clusters of galaxies based
on Spitzer/IRAC midinfrared data combined with shallow allsky optical data.
We identify distant cluster candidates in SSDF adopting an overdensity
threshold that results in a high purity (80%) cluster sample based on tests in
the Spitzer Deep, WideField Survey of the Bo\"otes field. Our simple algorithm
detects all three $1.4 < z \leq 1.75$ Xray detected clusters in the Bo\"otes
field. The uniqueness of the SSDF survey resides not just in its area, one of
the largest contiguous extragalactic fields observed with Spitzer, but also in
its deep, multiwavelength coverage by the South Pole Telescope (SPT),
Herschel/SPIRE and XMMNewton. This rich dataset will allow direct or stacked
measurements of SunyaevZel'dovich effect decrements or Xray masses for many
of the SSDF clusters presented here, and enable systematic study of the most
distant clusters on an unprecedented scale. We measure the angular correlation
function of our sample and find that these candidates show strong clustering.
Employing the COSMOS/UltraVista photometric catalog in order to infer the
redshift distribution of our cluster selection, we find that these clusters
have a comoving number density $n_c = (0.7^{+6.3}_{0.6}) \times 10^{7} h^{3}
\mathrm{Mpc}^{3}$ and a spatial clustering correlation scale length $r_0 = (32
\pm 7) h^{1} \rm{Mpc}$. Assuming our sample is comprised of dark matter halos
above a characteristic minimum mass, $M_{{\rm min}}$, we derive that at $z=1.5$
these clusters reside in halos larger than $M_{{\rm min}} = 1.5^{+0.9}_{0.7}
\times 10^{14} h^{1} M_{\odot}$. (abridged)

Planck's all sky surveys at 30857 GHz provide an unprecedented opportunity
to follow the radio spectra of a large sample of extragalactic sources to
frequencies 220 times higher than allowed by past, large area, groundbased
surveys. We combine the results of the Planck Early Release Compact Source
Catalog (ERCSC) with quasisimultaneous groundbased observations, as well as
archival data, at frequencies below or overlapping Planck frequency bands, to
validate the astrometry and photometry of the ERCSC radio sources and study the
spectral features shown in this new frequency window opened by Planck. The
ERCSC source positions and flux density scales are found to be consistent with
the groundbased observations. We present and discuss the spectral energy
distributions (SEDs) of a sample of "extreme" radio sources to illustrate the
richness of the ERCSC for the study of extragalactic radio sources. Variability
is found to play a role in the unusual spectral features of some of these
sources.

We present the XMMNewton followup for confirmation of Planck cluster
candidates. Twentyfive candidates have been observed to date using snapshot
(~10 ksec) exposures, ten as part of a pilot programme to sample a low range of
signaltonoise ratios (4<S/N<6), and a further 15 in a programme to observe a
sample of S/N>5 candidates. The sensitivity and spatial resolution of
XMMNewton allows unambiguous discrimination between clusters and false
candidates. The 4 false candidates have S/N <= 4.1. A total of 21 candidates
are confirmed as extended Xray sources. Seventeen are single clusters, the
majority of which are found to have highly irregular and disturbed morphologies
(about ~70%). The remaining four sources are multiple systems, including the
unexpected discovery of a supercluster at z=0.45. For 20 sources we are able to
derive a redshift estimate from the Xray Fe K line (albeit of variable
quality). The new clusters span the redshift range 0.09 <= z <= 0.54, with a
median redshift of z~0.37. A first determination is made of their Xray
properties including the characteristic size, which is used to improve the
estimate of the SZ Compton parameter, Y_SZ. The followup validation programme
has helped to optimise the Planck candidate selection process. It has also
provided a preview of the Xray properties of these newlydiscovered clusters,
allowing comparison with their SZ properties, and to the Xray and SZ
properties of known clusters observed in the Planck survey. Our results suggest
that Planck may have started to reveal a nonnegligible population of massive
dynamically perturbed objects that is underrepresented in Xray surveys.
However, despite their particular properties, these new clusters appear to
follow the Y_SZY_X relation established for Xray selected objects, where Y_X
is the product of the gas mass and temperature.

Aims: Develop and validate tools to estimate residual noise covariance in
Planck frequency maps. Quantify signal error effects and compare different
techniques to produce lowresolution maps.
Methods: We derive analytical estimates of covariance of the residual noise
contained in lowresolution maps produced using a number of mapmaking
approaches. We test these analytical predictions using Monte Carlo simulations
and their impact on angular power spectrum estimation. We use simulations to
quantify the level of signal errors incurred in different resolution
downgrading schemes considered in this work.
Results: We find an excellent agreement between the optimal residual noise
covariance matrices and Monte Carlo noise maps. For destriping mapmakers, the
extent of agreement is dictated by the knee frequency of the correlated noise
component and the chosen baseline offset length. The significance of signal
striping is shown to be insignificant when properly dealt with. In map
resolution downgrading, we find that a carefully selected window function is
required to reduce aliasing to the subpercent level at multipoles, ell >
2Nside, where Nside is the HEALPix resolution parameter. We show that
sufficient characterization of the residual noise is unavoidable if one is to
draw reliable contraints on large scale anisotropy.
Conclusions: We have described how to compute the lowresolution maps, with a
controlled sky signal level, and a reliable estimate of covariance of the
residual noise. We have also presented a method to smooth the residual noise
covariance matrices to describe the noise correlations in smoothed, bandwidth
limited maps.

The Planck satellite will observe the full sky at nine frequencies from 30 to
857 GHz. The goal of this paper is to examine the effects of four realistic
instrument systematics in the 30 GHz frequency maps: nonaxiallysymmetric
beams, sample integration, sorption cooler noise, and pointing errors. We
simulated one year long observations of four 30 GHz detectors. The simulated
timestreams contained CMB, foreground components (both galactic and
extragalactic), instrument noise (correlated and white), and the four
instrument systematic effects. We made maps from the timelines and examined the
magnitudes of the systematics effects in the maps and their angular power
spectra. We also compared the maps of different mapmaking codes to see how they
performed. We used five mapmaking codes (two destripers and three optimal
codes). None of our mapmaking codes makes an attempt to deconvolve the beam
from its output map. Therefore all our maps had similar smoothing due to beams
and sample integration. Temperature to polarization crosscoupling due to beam
mismatch causes a detectable bias in the TE spectrum of the CMB map. The
effects of cooler noise and pointing errors did not appear to be major concerns
for the 30 GHz channel. The only essential difference found so far between
mapmaking codes that affects accuracy (in terms of residual RMS) is baseline
length. All optimal codes give essentially indistinguishable results. A
destriper gives the same result as the optimal codes when the baseline is set
short enough. For longer baselines destripers require less computing resources
but deliver a noisier map.

This paper is one of a series describing the performance and accuracy of
mapmaking codes as assessed by the Planck CTP working group. We compare the
performance of multiple codes written by different groups for making polarized
maps from Plancksized, allsky cosmic microwave background (CMB) data. Three
of the codes are based on destriping algorithm, whereas the other three are
implementations of a maximumlikelihood algorithm. Previous papers in the
series described simulations at 100 GHz (Poutanen et al. 2006) and 217 GHz
(Ashdown et al. 2006). In this paper we make maps (temperature and
polarisation) from the simulated oneyear observations of four 30 GHz detectors
of Planck Low Frequency Instrument (LFI). We used Planck Level S simulation
pipeline to produce the observed timeordereddata streams (TOD). Our previous
studies considered polarisation observations for the CMB only. For this paper
we increased the realism of the simulations and included polarized galactic
foregrounds to our sky model. Our simulated TODs comprised of dipole, CMB,
diffuse galactic emissions, extragalactic radio sources, and detector noise.
The strong subpixel signal gradients arising from the foreground signals couple
to the output map through the mapmaking and cause an error (signal error) in
the maps. Destriping codes have smaller signal error than the
maximumlikelihood codes. We examined a number of schemes to reduce this error.
On the other hand, the maximumlikelihood mapmaking codes can produce maps
with lower residual noise than destriping codes.

We perform a Fisher matrix analysis to quantify cosmological constraints
obtainable from a 2dimensional SunyaevZel'dovich (SZ) cluster catalog using
the counts and the angular correlation function. Three kinds of SZ survey are
considered: the almost allsky Planck survey and two deeper groundbased
surveys, one with 10% sky coverage, the other one with a coverage of 250 square
degrees. With the counts and angular function, and adding the constraint from
the local Xray cluster temperature function, joint 10% to 30% errors (1 sigma)
are achievable on the cosmological parameter pair (sigma_8, Omega_m) in the
flat concordance model. Constraints from a 2D distribution remain relatively
robust to uncertainties in possible cluster gas evolution for the case of
Planck. Alternatively, we examine constraints on cluster gas physics when
assuming priors on the cosmological parameters (e.g., from cosmic microwave
background anisotropies and SNIa data), finding a poor ability to constrain gas
evolution with the 2dimensional catalog. From just the SZ counts and angular
correlation function we obtain, however, a constraint on the product between
the presentday cluster gas mass fraction and the normalization of the
masstemperature relation, T_*, with a precision of 15%. This is particularly
interesting because it would be based on a very large catalog and is
independent of any Xray data.

(abridged) We describe XMMNewton Guaranteed Time observations of a sample of
eight high redshift (0.45<z<0.62) clusters. The goal of these observations was
to measure the luminosity and the temperature of the clusters to a precision of
\~10%, leading to constraints on the possible evolution of the
luminositytemperature relation, and ultimately on the values of the matter
density, Omega_M and, to a lesser extent, the cosmological constant Omega_L.
The clusters were drawn from the SHARC and 160 Square Degree (160SD) ROSAT
surveys. Here we describe our data analysis techniques and present, for the
first time with XMMNewton,LxTx relation. For each of the eight clusters in
the sample, we have measured total bolometric luminosities, performed
betamodel fits to the radial surface profiles and made spectral fits to a
single temperature isothermal model. We describe data analysis techniques that
pay particular attention to background mitigation. Characterizing the LxTx
relation as Lx = L_{6} (T/6keV)^{alpha},we find L_{6}=16.8 +7.6/5.2 10^{44}
erg/s and alpha=2.7 +/0.4 for a EdS H=50 cosmology at a typical redshift z
=0.55. Comparing with the low redshift study by Markevitch, assuming LT to
evolve as (1+z)^A, we find A=0.68 +/0.26 for the same cosmology and A=1.52
+0.26/0.27 for a concordance cosmology. We conclude that there is now evidence
from both XMMNewton and Chandra for an evolutionary trend in the LT relation.
Our observations lend support to the robustness and completeness of the SHARC
and 160SD surveys.

The evolution with redshift of the temperatureluminosity relation of Xray
galaxy clusters is a key ingredient to break degeneracies in the interpretation
of Xray clusters redshift number counts. We therefore take advantage of the
recent measurements of the temperatureluminosity relation of distant clusters
observed with XMMNewton and Chandra satellites to examine theoretical number
counts expected for different available Xrays cluster samples, namely the
RDCS, EMSS, SHARC, 160deg^2 and the MACS at redshift greater than 0.3. We
derive these counts without any adjustment, using models previously normalized
to the local temperature distribution function and to the highz (z = 0.33)
TDF. We find that these models having Omega_M in the range [0.851.] predict
counts in remarkable agreement with the observed counts in the different
samples. We illustrate that this conclusion is weakly sensitive to the various
ingredients of the modeling. Therefore number counts provide a robust evidence
of an evolving population. A realistic flat low density model (Omega_M = 0.3),
normalized to the local abundance of clusters is found to overproduce cluster
abundance at high redshift (above z = 0.5) by nearly an order of magnitude.
This result is in conflict with the popular concordance model. The conflict
could indicate a deviation from the expected scaling of the MT relation with
redshift.

We discuss the angular correlation function of SunyaevZel'dovich
(SZ)detected galaxy clusters as a cosmological probe. As a projection of the
realspace cluster correlation function, the angular function samples the
underlying SZ catalog redshift distribution. It offers a way to study cosmology
and cluster evolution directly with the twodimensional catalog, even before
extensive followup observations, thereby facilitating the immediate scientific
return from SZ surveys. As a simple illustration of the information content of
the angular function, we examine its dependence on the parameter pair Om_m,
sigma_8 in flat cosmologies. We discuss sources of modeling uncertainty and
consider application to the future Planck SZ catalog, showing how these two
parameters and the normalization of the SZ fluxmass relation can be
simultaneously found when the local Xray cluster abundance constraint is
included.

Application of a Goodnessoffit (GOF) statistic is an essential element of
parameter estimation. We discuss the computation of GOF when estimating
parameters from anisotropy measurements of the cosmic microwave background
(CMB), and we propose two GOF statistics to be used when employing approximate
bandpower likelihood functions. They are based on an approximate form for the
distribution of bandpower estimators that requires only minimal experimental
information to construct. Monte Carlo simulations of CMB experiments show that
the proposed form describes the true distributions quite well. We apply these
GOF statistics to current CMB anisotropy data and discuss the results.

We present the results of a combined analysis of cosmic microwave background
(CMB) and Xray galaxy clusters baryon fraction to deduce constraints over 6 in
flationnary cosmological parameters. Such a combination is necessary for breaki
ng degeneracies inherent to the CMB.

I discuss galaxy cluster surveys based on the SunyaevZel'dovich effect and
their relevance for cosmological studies. The unique aspects of cluster
selection by this method are emphasized and certain issues of surveying are
addressed. Finally, I briefly present prospects for upcoming surveys.

We describe a fast method to simulate in a semianalytical way consistent
maps of the thermal, kinetic and polarised SunyaevZel'dovich effect, featuring
both cluster spatial correlations and large scale velocity flows.

The abundance of highredshift galaxy clusters depends sensitively on the
matter density $\OmM$ and, to a lesser extent, on the cosmological constant
$\Lambda$. Measurements of this abundance therefore constrain these fundamental
cosmological parameters, and in a manner independent and complementary to other
methods, such as observations of the cosmic microwave background and distance
measurements. Cluster abundance is best measured by the Xray temperature
function, as opposed to luminosity, because temperature and mass are tightly
correlated, as demonstrated by numerical simulations. Taking advantage of the
sensitivity of XMMNewton, our Guaranteed Time program aims at measuring the
temperature of the highest redshift (z>0.4) SHARC clusters, with the ultimate
goal of constraining both $\OmM$ and $\Lambda$.

Most parameter constraints obtained from cosmic microwave background (CMB)
anisotropy data are based on power estimates and rely on approximate likelihood
functions; computational difficulties generally preclude an exact analysis
based on pixel values. With the specific goal of testing this kind of approach,
we have performed a complete (unapproximated) likelihood analysis combining
the COBE, Saskatoon and MAX data sets. We examine in detail the ability of
certain approximate techniques based on bandpower estimates to recover the
full likelihood constraints. The traditional $\chi^2$method does not always
find the same bestfit model as the likelihood analysis (a bias), due mainly to
the false assumption of Gaussian likelihoods that makes the method overly
sensitive to data outliers. Although an improvement, other approaches employing
nonGaussian flatband likelihoods do not always faithfully reproduce the
complete likelihood constraints either; not even when using the exact flatband
likelihood curves. We trace this to the neglect of spectral information by
simple flat bandpower estimates. A straightforward extension incorporating a
local effective slope (of the power spectrum, $C_l$) provides a faithful
representation of the likelihood surfaces without significantly increasing
computing cost. Finally, we also demonstrate that the bestfit model to this
particular data set is a {\em good fit}, or that the observations are
consistent with Gaussian sky fluctuations, according to our statistic.

Using an approximate likelihood method adapted to bandpower estimates, we
analyze the ensemble of first generation cosmic microwave background anisotropy
experiments to deduce constraints over a sixdimensional parameter space
describing Inflationgenerated adiabatic, scalar fluctuations. The basic
preferences of simple Inflation scenarios are consistent with the data set:
flat geometries $(\OmT \equiv 1\Omk \sim 1)$ and a scaleinvariant primeval
spectrum ($n\sim 1$) are favored. Models with significant negative curvature
($\OmT < 0.7$) are eliminated, while constraints on postive curvature are less
stringent. Degeneracies among the parameters prevent independent determinations
of the matter density $\OmM$ and the cosmological constant $\Lambda$, and the
Hubble constant $\Ho$ remains relatively unconstrained. We also find that the
height of the first Doppler peak relative to the amplitude suggested by data at
larger $l$ indicates a high baryon content ($\Omb h^2$), almost independently
of the other parameters. Besides the overall qualitative advance expected of
the next generation experiments, their improved dipole calibrations will be
particularly useful for constraining the peak height. Our analysis includes a
{\em GoodnessofFit} statistic applicable to power estimates and which
indicates that the maximum likelihood model provides an acceptable fit to the
data set.

Bandpower estimates of cosmic microwave background fluctuations are now
routinely used to place constraints on cosmological parameters. For this to be
done in a rigorous fashion, the full likelihood function of bandpower
estimates must be employed. Even for Gaussian theories, this likelihood
function is not itself Gaussian, for the simple reason that bandpowers measure
the {\em variance} of the random sky fluctuations. In the context of Gaussian
sky fluctuations, we use an ideal situation to motivate a general form for the
full likelihood function from a given experiment. This form contains only two
free parameters, which can be determined if the 68% and 95% confidence
intervals of the true likelihood function are known. The ansatz works
remarkably well when compared to the complete likelihood function for a number
of experiments. For application of this kind of approach, we suggest that in
the future both 68% and 95% (and perhaps also the 99.7%) confidence intervals
be given when reporting experimental results.

(abridged) We present a new determination of the local temperature function
of Xray clusters. We use a new sample comprising fifty clusters for which
temperature information is now available, making it the largest complete sample
of its kind. It is therefore expected to significantly improve the estimation
of the temperature distribution function of moderately hot clusters. We find
that the resulting temperature function is higher than previous estimations,
but agrees well with the temperature distribution function inferred from the
BCS and RASS luminosity function. We have used this sample to constrain the
amplitude of the matter fluctuations on cluster's scale of
$8\sqrt[3]{\Omega_0}^{1}h^{1}$Mpc, assuming a masstemperature relation based
on recent numerical simulations. We find $\sigma_8 = 0.6\pm 0.02$ for an
$\Omega_0 = 1$ model. Our sample provides an ideal reference at $z \sim 0$ to
use in the application of the cosmological test based on the evolution of Xray
cluster abundance (Oukbir & Blanchard 1992, 1997). Using Henry's sample, we
find that the abundance of clusters at $z = 0.33$ is significantly smaller, by
a factor larger than 2, which shows that the EMSS sample provides strong
evidence for evolution of the cluster abundance. A likelihood analysis leads to
a rather high value of the mean density parameter of the universe: $\Omega
=0.92 \pm 0.22$ (open case) and $\Omega =0.86 \pm 0.25$ (flat case), which is
consistent with a previous, independent estimation based on the full EMSS
sample by Sadat et al.(1998). Some systematic uncertainties which could alter
this result are briefly discussed.

The majority of present efforts to constrain cosmological parameters with
cosmic microwave background (CMB) anisotropy data employ approximate likelihood
functions, the time consuming nature of a complete analysis being a major
obstacle. We have performed a full (unapproximated) likelihood analysis on
several experiments that allows us to examine the various assumptions made in
these approximate methods and to evaluate their performance. Our results
indicate that care should be taken when using such approaches. With an improved
approximate method, we present some constraints on cosmological parameters
using the entire present CMB data set.

We present analyses of the twopoint correlation properties of the ESP galaxy
redshift survey. From the redshiftspace correlation function xi(s), we see
positive clustering out to separations ~50/h Mpc, with a smooth break on larger
scales and zerocrossing between 60 and 80/h Mpc. xi(s) is reasonably well
described by a shallow power law with \gamma~1.5 between 3 and 50/h Mpc, while
on smaller scales (0.22/h Mpc) it has a shallower slope (\gamma~ 1). We
examine the full effect of redshiftspace distortions through the
twodimensional correlation function xi(rp,pi), from which we project out the
realspace xi(r) below 10/h Mpc. This function is well described by a powerlaw
model (r/r_o)^{\gamma}, with r_o=4.15^{+0.20}_{0.21} h^{1} Mpc and
\gamma=1.67^{+0.07}_{0.09}. Comparison to other redshift surveys shows a
consistent picture in which clustering remains positive out to separations of
50/h Mpc or larger, in substantial agreement with the results obtained from
angular surveys like the APM and EDSGC. Also the shape of the twopoint
correlation function is remarkably unanimous among these data sets, in all
cases requiring more power above 5/h Mpc (a `shoulder'), than a simple
extrapolation of the canonical xi(r)=(r/5)^{1.8}. xi(s) for volumelimited
subsamples shows evidence of luminosity segregation only for the most luminous
sample with M_{b_J}\le 20.5. When redshiftspace distortions are removed
through projection of xi(rp,pi), however, a weak dependence on luminosity is
seen at small separations also at fainter magnitudes. This effect is masked in
redshift space, as the mean pairwise velocity dispersion experiences a parallel
increase, basically erasing the effect of the clustering growth on xi(s).

We discuss the problem of constraining cosmological parameters with cosmic
microwave background bandpower estimates. Because these latter are variances,
they do not have gaussian distribution functions and, hence, the standard
$\chi^2$approach is not strictly applicable. A general purpose approximation
to experimental bandpower likelihood functions is proposed, which requires
only limited experimental details. Comparison with the full likelihood function
calculated for several experiments shows that the approximation works well.

We discuss the constraints one can place on cosmological parameters using
current cosmic microwave background data. A standard $\chi^2$minimization
over bandpower estimates is first presented, followed by a discussion of the
more correct likelihood approach. We propose an approximation to the complete
likelihood function of an arbitrary experiment requiring only limited and
easily found information about the observations. Examination of both open
models  $(\Omega,h,Q,n)$  and flat models ($\Omega+\Omega_\Lambda=1$) 
$(\Omega,\Omega_b,h,Q,n)$  leaves one rather robust result: models with small
curvature are favored.