QUIJOTE (Q-U-I JOint TEnerife) is an experiment designed to achieve CMB B-mode polarization detection and sensitive enough to detect a primordial gravitational-wave component if the B-mode amplitude is larger than r = 0.05. It consists in two telescopes and three instruments observing in the frequency range 10-42 GHz installed at the Teide Observatory in the Canary Islands, Spain. The observing strategy includes three raster scan deep integration fields for cosmology, a nominal wide survey covering the Northen Sky and specific raster scan deep integration observations in regions of specific interest. The main goals of the project are presented and the first scientific results obtained with the first instrument are reviewed.

We present the implications for cosmic inflation of the Planck measurements of the cosmic microwave background (CMB) anisotropies in both temperature and polarization based on the full Planck survey. The Planck full mission temperature data and a first release of polarization data on large angular scales measure the spectral index of curvature perturbations to be $n_\mathrm{s} = 0.968 \pm 0.006$ and tightly constrain its scale dependence to $d n_s/d \ln k =-0.003 \pm 0.007$ when combined with the Planck lensing likelihood. When the high-$\ell$ polarization data is included, the results are consistent and uncertainties are reduced. The upper bound on the tensor-to-scalar ratio is $r_{0.002} < 0.11$ (95% CL), consistent with the B-mode polarization constraint $r< 0.12$ (95% CL) obtained from a joint BICEP2/Keck Array and Planck analysis. These results imply that $V(\phi) \propto \phi^2$ and natural inflation are now disfavoured compared to models predicting a smaller tensor-to-scalar ratio, such as $R^2$ inflation. Three independent methods reconstructing the primordial power spectrum are investigated. The Planck data are consistent with adiabatic primordial perturbations. We investigate inflationary models producing an anisotropic modulation of the primordial curvature power spectrum as well as generalized models of inflation not governed by a scalar field with a canonical kinetic term. The 2015 results are consistent with the 2013 analysis based on the nominal mission data.

We discuss the effects on the CMB, CIB, and thermal SZ effect due to the peculiar motion of an observer with respect to the CMB rest frame, which induces boosting effects. We investigate the scientific perspectives opened by future CMB space missions, focussing on the CORE proposal. The improvements in sensitivity offered by a mission like CORE, together with its high resolution over a wide frequency range, will provide a more accurate estimate of the CMB dipole. The extension of boosting effects to polarization and cross-correlations will enable a more robust determination of purely velocity-driven effects that are not degenerate with the intrinsic CMB dipole, allowing us to achieve a S/N ratio of 13; this improves on the Planck detection and essentially equals that of an ideal cosmic-variance-limited experiment up to a multipole l of 2000. Precise inter-frequency calibration will offer the opportunity to constrain or even detect CMB spectral distortions, particularly from the cosmological reionization, because of the frequency dependence of the dipole spectrum, without resorting to precise absolute calibration. The expected improvement with respect to COBE-FIRAS in the recovery of distortion parameters (in principle, a factor of several hundred for an ideal experiment with the CORE configuration) ranges from a factor of several up to about 50, depending on the quality of foreground removal and relative calibration. Even for 1% accuracy in both foreground removal and relative calibration at an angular scale of 1 deg, we find that dipole analyses for a mission like CORE will be able to improve the recovery of the CIB spectrum amplitude by a factor of 17 in comparison with current results based on FIRAS. In addition to the scientific potential of a mission like CORE for these analyses, synergies with other planned and ongoing projects are also discussed.

We present an analysis of the main systematic effects that could impact the measurement of CMB polarization with the proposed CORE space mission. We employ timeline-to-map simulations to verify that the CORE instrumental set-up and scanning strategy allow us to measure sky polarization to a level of accuracy adequate to the mission science goals. We also show how the CORE observations can be processed to mitigate the level of contamination by potentially worrying systematics, including intensity-to-polarization leakage due to bandpass mismatch, asymmetric main beams, pointing errors and correlated noise. We use analysis techniques that are well validated on data from current missions such as Planck to demonstrate how the residual contamination of the measurements by these effects can be brought to a level low enough not to hamper the scientific capability of the mission, nor significantly increase the overall error budget. We also present a prototype of the CORE photometric calibration pipeline, based on that used for Planck, and discuss its robustness to systematics, showing how CORE can achieve its calibration requirements. While a fine-grained assessment of the impact of systematics requires a level of knowledge of the system that can only be achieved in a future study phase, the analysis presented here strongly suggests that the main areas of concern for the CORE mission can be addressed using existing knowledge, techniques and algorithms.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

We describe a space-borne, multi-band, multi-beam polarimeter aiming at a precise and accurate measurement of the polarization of the Cosmic Microwave Background. The instrument is optimized to be compatible with the strict budget requirements of a medium-size space mission within the Cosmic Vision Programme of the European Space Agency. The instrument has no moving parts, and uses arrays of diffraction-limited Kinetic Inductance Detectors to cover the frequency range from 60 GHz to 600 GHz in 19 wide bands, in the focal plane of a 1.2 m aperture telescope cooled at 40 K, allowing for an accurate extraction of the CMB signal from polarized foreground emission. The projected CMB polarization survey sensitivity of this instrument, after foregrounds removal, is 1.7 {\mu}K$\cdot$arcmin. The design is robust enough to allow, if needed, a downscoped version of the instrument covering the 100 GHz to 600 GHz range with a 0.8 m aperture telescope cooled at 85 K, with a projected CMB polarization survey sensitivity of 3.2 {\mu}K$\cdot$arcmin.

We discuss the potential of a next generation space-borne 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 Medium-size 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 high-z galaxies distributed over the full sky. The extreme brightness of these galaxies will make it possible to study them, via follow-up observations, in extraordinary detail. Also, the CORE resolution matches the typical sizes of high-z galaxy proto-clusters 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 sub-mm wavelengths, respectively.

We forecast the main cosmological parameter constraints achievable with the CORE space mission which is dedicated to mapping the polarisation of the Cosmic Microwave Background (CMB). CORE was recently submitted in response to ESA's fifth call for medium-sized mission proposals (M5). Here we report the results from our pre-submission study of the impact of various instrumental options, in particular the telescope size and sensitivity level, and review the great, transformative potential of the mission as proposed. Specifically, we assess the impact on a broad range of fundamental parameters of our Universe as a function of the expected CMB characteristics, with other papers in the series focusing on controlling astrophysical and instrumental residual systematics. In this paper, we assume that only a few central CORE frequency channels are usable for our purpose, all others being devoted to the cleaning of astrophysical contaminants. On the theoretical side, we assume LCDM as our general framework and quantify the improvement provided by CORE over the current constraints from the Planck 2015 release. We also study the joint sensitivity of CORE and of future Baryon Acoustic Oscillation and Large Scale Structure experiments like DESI and Euclid. Specific constraints on the physics of inflation are presented in another paper of the series. In addition to the six parameters of the base LCDM, which describe the matter content of a spatially flat universe with adiabatic and scalar primordial fluctuations from inflation, we derive the precision achievable on parameters like those describing curvature, neutrino physics, extra light relics, primordial helium abundance, dark matter annihilation, recombination physics, variation of fundamental constants, dark energy, modified gravity, reionization and cosmic birefringence. (ABRIDGED)

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 ground-based 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 Sunyaev-Zeldovich 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 flux-mass relation; assuming self-similar 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 results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter inflationary LCDM cosmology. From the Planck temperature and lensing data, for this cosmology we find a Hubble constant, H0= (67.8 +/- 0.9) km/s/Mpc, a matter density parameter Omega_m = 0.308 +/- 0.012 and a scalar spectral index with n_s = 0.968 +/- 0.006. (We quote 68% errors on measured parameters and 95% limits on other parameters.) Combined with Planck temperature and lensing data, Planck LFI polarization measurements lead to a reionization optical depth of tau = 0.066 +/- 0.016. Combining Planck with other astrophysical data we find N_ eff = 3.15 +/- 0.23 for the effective number of relativistic degrees of freedom and the sum of neutrino masses is constrained to < 0.23 eV. Spatial curvature is found to be |Omega_K| < 0.005. For LCDM we find a limit on the tensor-to-scalar ratio of r <0.11 consistent with the B-mode constraints from an analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP data leads to a tighter constraint of r < 0.09. We find no evidence for isocurvature perturbations or cosmic defects. The equation of state of dark energy is constrained to w = -1.006 +/- 0.045. Standard big bang nucleosynthesis predictions for the Planck LCDM cosmology are in excellent agreement with observations. We investigate annihilating dark matter and deviations from standard recombination, finding no evidence for new physics. The Planck results for base LCDM are in agreement with BAO data and with the JLA SNe sample. However the amplitude of the fluctuations is found to be higher than inferred from rich cluster counts and weak gravitational lensing. Apart from these tensions, the base LCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.

(abridged) We discuss the Galactic foreground emission between 20 and 100GHz based on observations by Planck/WMAP. The Commander component-separation tool has been used to separate the various astrophysical processes in total intensity. Comparison with RRL templates verifies the recovery of the free-free emission along the Galactic plane. Comparison of the high-latitude Halpha emission with our free-free map shows residuals that correlate with dust optical depth, consistent with a fraction (~30%) of Halpha having been scattered by high-latitude dust. We highlight a number of diffuse spinning dust morphological features at high latitude. There is substantial spatial variation in the spinning dust spectrum, with the emission peak ranging from below 20GHz to more than 50GHz. There is a strong tendency for the spinning dust component near many prominent HII regions to have a higher peak frequency, suggesting that this increase in peak frequency is associated with dust in the photodissociation regions around the nebulae. The emissivity of spinning dust in these diffuse regions is of the same order as previous detections in the literature. Over the entire sky, the commander solution finds more anomalous microwave emission than the WMAP component maps, at the expense of synchrotron and free-free emission. This can be explained by the difficulty in separating multiple broadband components with a limited number of frequency maps. Future surveys (5-20GHz), will greatly improve the separation by constraining the synchrotron spectrum. We combine Planck/WMAP data to make the highest S/N ratio maps yet of the intensity of the all-sky polarized synchrotron emission at frequencies above a few GHz. Most of the high-latitude polarized emission is associated with distinct large-scale loops and spurs, and we re-discuss their structure...

We present all-sky modelling of the high resolution Planck, IRAS, and WISE infrared (IR) observations using the physical dust model presented by Draine and Li in 2007 (DL). We study the performance and results of this model, and discuss implications for future dust modelling. The present work extends the DL dust modelling carried out on nearby galaxies using Herschel and Spitzer data to Galactic dust emission. We employ the DL dust model to generate maps of the dust mass surface density, the optical extinction Av, and the starlight intensity parametrized by Umin. The DL model reproduces the observed spectral energy distribution (SED) satisfactorily over most of the sky, with small deviations in the inner Galactic disk and in low ecliptic latitude areas. We compare the DL optical extinction Av for the diffuse interstellar medium with optical estimates for 2 10^5 quasi-stellar objects (QSOs) observed in the Sloan digital sky survey. The DL Av estimates are larger than those determined towards QSOs by a factor of about 2, which depends on Umin. The DL fitting parameter Umin, effectively determined by the wavelength where the SED peaks, appears to trace variations in the far-IR opacity of the dust grains per unit Av, and not only in the starlight intensity. To circumvent the model deficiency, we propose an empirical renormalization of the DL Av estimate, dependent of Umin, which compensates for the systematic differences found with QSO observations. This renormalization also brings into agreement the DL Av estimates with those derived for molecular clouds from the near-IR colours of stars in the 2 micron all sky survey. The DL model and the QSOs data are used to compress the spectral information in the Planck and IRAS observations for the diffuse ISM to a family of 20 SEDs normalized per Av, parameterized by Umin, which may be used to test and empirically calibrate dust models.

The Planck design and scanning strategy provide many levels of redundancy that can be exploited to provide tests of internal consistency. One of the most important is the comparison of the 70GHz and 100GHz channels. Based on different instrument technologies, with feeds located differently in the focal plane, analysed independently by different teams using different software, and near the minimum of diffuse foreground emission, these channels are in effect two different experiments. The 143GHz channel has the lowest noise level on Planck, and is near the minimum of unresolved foreground emission. In this paper, we analyse the level of consistency achieved in the 2013 Planck data. We concentrate on comparisons between the 70/100/143GHz channel maps and power spectra, particularly over the angular scales of the first and second acoustic peaks, on maps masked for diffuse Galactic emission and for strong unresolved sources. Difference maps covering angular scales from 8deg-15arcmin are consistent with noise, and show no evidence of cosmic microwave background structure. Including small but important corrections for unresolved-source residuals, we demonstrate agreement between 70 and 100GHz power spectra averaged over 70<l<390 at the 0.8% level, and agreement between 143 and 100GHz power spectra of 0.4% over the same l range. These values are within and consistent with the overall uncertainties in calibration given in the Planck 2013 results. We also present results based on the 2013 likelihood analysis showing consistency at the 0.35% between the 100/143/217GHz power spectra. We analyse calibration procedures and beams to determine what fraction of these differences can be accounted for by known approximations or systematic errors that could be controlled even better in the future, reducing uncertainties still further. Several possible small improvements are described...(abridged)

The Planck High Frequency Instrument (HFI) has observed the full sky at six frequencies (100, 143, 217, 353, 545, and 857 GHz) in intensity and at four frequencies in linear polarization (100, 143, 217, and 353 GHz). In order to obtain sky maps, the time-ordered information (TOI) containing the detector and pointing samples must be processed and the angular response must be assessed. The full mission TOI is included in the Planck 2015 release. This paper describes the HFI TOI and beam processing for the 2015 release. HFI calibration and map-making are described in a companion paper. The main pipeline has been modified since the last release (2013 nominal mission in intensity only), by including a correction for the non-linearity of the warm readout and by improving the model of the bolometer time response. The beam processing is an essential tool that derives the angular response used in all the Planck science papers and we report an improvement in the effective beam window function uncertainty of more than a factor 10 relative to the 2013 release. Noise correlations introduced by pipeline filtering function are assessed using dedicated simulations. Angular cross-power spectra using datasets that are decorrelated in time are immune to the main systematic effects.

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

We characterize the Cosmic Infrared Background (CIB)-lensing bispectrum which is one of the contributions to the three-point functions of Cosmic Microwave Background (CMB) maps in harmonic space. We show that the CIB-lensing bispectrum has a considerable strength and that it can be detected with high significance in the Planck high-frequency maps. We also present forecasts of the contamination on different shapes of the primordial non-Gaussianity $f_{nl}$ parameter produced by the CIB-lensing bispectrum and by the extragalactic point sources bispectrum in the Planck high-resolution CMB anisotropy maps. The local, equilateral and orthogonal shapes are considered for 'raw' single-frequency (i.e., without applying any component separation technique) and foreground-reduced Planck temperature maps. The CIB-lensing correlation seems to mainly affect orthogonal shapes of the bispectrum - with {\Delta}$f_{nl}$ = -21 and -88 for the 143 and 217 GHz bands respectively - while point sources mostly impact equilateral shapes, with {\Delta}$f_{nl}$ = 160, 54 and 60 at 100, 143 and 217 GHz. However, the results indicate that these contaminants do not induce any relevant bias on Planck $f_{nl}$ estimates when foreground-reduced maps are considered: using SEVEM for the component separation we obtain {\Delta}$f_{nl}$ = 10.5 due to the CIB-lensing and {\Delta}$f_{nl}$ = 30.4 due to point sources, corresponding to 0.30{\sigma} and 0.45{\sigma} in terms of the Planck 2013 $f_{nl}$ uncertainty. The component separation technique is, in fact, able to partially clean the extragalactic source contamination and the bias is reduced for all the shapes. We have further developed single- and multiple-frequency estimators based on the Komatsu, Spergel & Wandelt (2005) formalism that can be implemented to efficiently detect this signal.

Shortened abstract: Observations of the nearby Chamaeleon clouds in gamma rays with the Fermi Large Area Telescope and in thermal dust emission with Planck and IRAS have been used with the HI and CO radio data to (i) map the gas column densities in the different phases and at the dark neutral medium (DNM) transition between the HI-bright and CO-bright media; (ii) constrain the CO-to-$H_2$ conversion factor, $X_{CO}$; (iii) probe the dust properties per gas nucleon in each gas phase and spatially across the clouds. We have separated clouds in velocity in HI and CO emission and modelled the 0.4-100 GeV intensity, the dust optical depth at 353 GHz, the thermal radiance of the large grains, and an estimate of the dust extinction empirically corrected for the starlight intensity, $A_{VQ}$. The gamma-ray emissivity spectra confirm that the GeV-TeV cosmic rays uniformly permeate all gas phases up to the CO cores. The dust and cosmic rays reveal large amounts of DNM gas, with comparable spatial distributions and twice as much mass as in the CO-bright clouds. We give constraints on the HI-DNM-CO transitions and CO-dark $H_2$ fractions for separate clouds. The corrected extinction provides the best fit to the total gas traced by the gamma rays, but we find evidence for a rise in $A_{VQ}/N_H$ and a steep rise in opacity, with increasing $N_H$ and $H_2$ fraction, and with decreasing dust temperature. We observe less variations for the specific power of the grains, except for a decline by half in the CO cores. This combined information suggests grain evolution. The gamma rays and dust radiance yield consistent $X_{CO}$ estimates near $0.7\times10^{20}$ cm$^{-2}$ (K km/s)$^{-1}$. The other dust tracers yield biased values because of the grain opacity rise in the CO clouds. These results also confirm a factor of 2 difference between $X_{CO}$ estimates at pc and kpc scales.

We analyse the implications of the Planck data for cosmic inflation. The Planck nominal mission temperature anisotropy measurements, combined with the WMAP large-angle polarization, constrain the scalar spectral index to $n_s = 0.9603 \pm 0.0073$, ruling out exact scale invariance at over 5 $\sigma$. Planck establishes an upper bound on the tensor-to-scalar ratio of r < 0.11 (95% CL). The Planck data thus shrink the space of allowed standard inflationary models, preferring potentials with V" < 0. Exponential potential models, the simplest hybrid inflationary models, and monomial potential models of degree n > 2 do not provide a good fit to the data. Planck does not find statistically significant running of the scalar spectral index, obtaining $d n_s/d ln k = -0.0134 \pm 0.0090$. Several analyses dropping the slow-roll approximation are carried out, including detailed model comparison and inflationary potential reconstruction. We also investigate whether the primordial power spectrum contains any features. We find that models with a parameterized oscillatory feature improve the fit $\chi^2$ by ~ 10; however, Bayesian evidence does not prefer these models. We constrain several single-field inflation models with generalized Lagrangians by combining power spectrum data with bounds on $f_\mathrm{NL}$ measured by Planck. The fractional primordial contribution of CDM isocurvature modes in the curvaton and axion scenarios has upper bounds of 0.25% or 3.9% (95% CL), respectively. In models with arbitrarily correlated CDM or neutrino isocurvature modes, an anticorrelation can improve $\chi^2$ by approximatively 4 as a result of slightly lowering the theoretical prediction for the $\ell<40$ multipoles relative to the higher multipoles. Nonetheless, the data are consistent with adiabatic initial conditions.

We have constructed the first all-sky map of the thermal Sunyaev-Zeldovich (tSZ) effect by applying specifically tailored component separation algorithms to the 100 to 857 GHz frequency channel maps from the Planck survey. These maps show an obvious galaxy cluster tSZ signal that is well matched with blindly detected clusters in the Planck SZ catalogue. To characterize the signal in the tSZ map we have computed its angular power spectrum. At large angular scales ($\ell < 60$), the major foreground contaminant is the diffuse thermal dust emission. At small angular scales ($\ell > 500$) the clustered Cosmic Infrared Background (CIB) and residual point sources are the major contaminants. These foregrounds are carefully modelled and subtracted. We measure the tSZ power spectrum in angular scales, $0.17^{\circ} \lesssim \theta \lesssim 3.0^{\circ}$, that were previously unexplored. The measured tSZ power spectrum is consistent with that expected from the Planck catalogue of SZ sources, with additional clear evidence of signal from unresolved clusters and, potentially, diffuse warm baryons. We use the tSZ power spectrum to obtain the following cosmological constraints: $\sigma_8(\Omega_{\mathrm{m}}/0.28)^{3.2/8.1}=0.784 \pm 0.016 (68% C.L.). Marginalized band-powers of the Planck tSZ power spectrum and the best-fit model are given. The non-Gaussianity of the Compton parameter map is further characterized by computing its 1D probability distribution function and its bispectrum. These are used to place additional independent constraints on $\sigma_{8}$.

Using \Planck\ data combined with the Meta Catalogue of X-ray detected Clusters of galaxies (MCXC), we address the study of peculiar motions by searching for evidence of the kinetic Sunyaev-Zeldovich effect (kSZ). By implementing various filters designed to extract the kSZ generated at the positions of the clusters, we obtain consistent constraints on the radial peculiar velocity average, root mean square (rms), and local bulk flow amplitude at different depths. For the whole cluster sample of average redshift 0.18, the measured average radial peculiar velocity with respect to the cosmic microwave background (CMB) radiation at that redshift, i.e., the kSZ monopole, amounts to $72 \pm 60$ km s$^{-1}$. This constitutes less than 1% of the relative Hubble velocity of the cluster sample with respect to our local CMB frame. While the linear $\Lambda$CDM prediction for the typical cluster radial velocity rms at $z=0.15$ is close to 230km s$^{-1}$, the upper limit imposed by \Planck\ data on the cluster subsample corresponds to 800 km s$^{-1}$ at 95% confidence level, i.e., about three times higher. \Planck\ data also set strong constraints on the local bulk flow in volumes centred on the Local Group. There is no detection of bulk flow as measured in any comoving sphere extending to the maximum redshift covered by the cluster sample. A blind search for bulk flows in this sample has an upper limit of 254 km s$^{-1}$ (95% confidence level) dominated by CMB confusion and instrumental noise, indicating that the Universe is largely homogeneous on Gpc scales. In this context, in conjunction with supernova observations, \Planck\ is able to rule out a large class of inhomogeneous void models as alternatives to dark energy or modified gravity. The \Planck\ constraints on peculiar velocities and bulk flows are thus consistent with the $\Lambda$CDM scenario.

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in May 2013 as a large-class mission for investigating within the framework of the ESA Cosmic Vision program a set of important scientific questions that require high resolution, high sensitivity, full-sky observations of the sky emission at wavelengths ranging from millimeter-wave to the far-infrared. PRISM's main objective is to explore the distant universe, probing cosmic history from very early times until now as well as the structures, distribution of matter, and velocity flows throughout our Hubble volume. PRISM will survey the full sky in a large number of frequency bands in both intensity and polarization and will measure the absolute spectrum of sky emission more than three orders of magnitude better than COBE FIRAS. The aim of this Extended White Paper is to provide a more detailed overview of the highlights of the new science that will be made possible by PRISM

In the present paper we study the radio sources in the NRAO VLA SKY Survey (NVSS) analysing its power spectrum and galaxy distribution. There is a discrepancy between the theoretical models in the literature and the data at large scales. A new model for NVSS is proposed combining the power spectrum data from NVSS and the galaxy distribution from the Combined EIS-NVSS Survey Of Radio Sources (CENSORS). Taken into account these two data sets the differences in the power spectrum at large scales are reduced, but there is still some tension between the two data sets. Different models are compared using Bayesian evidence. A model for the galaxy distribution based on a gamma function provides a higher evidence against other models proposed in the literature. In addition, the effect of primordial non-Gaussianity has been also considered. The 2-sigma constraints on the local non-Gaussian f_NL parameter are -43 < f_NL < 142.

Rotational transition lines of CO play a major role in molecular radio astronomy and in particular in the study of star formation and the Galactic structure. Although a wealth of data exists in the Galactic plane and some well-known molecular clouds, there is no available CO high sensitivity all-sky survey to date. Such all-sky surveys can be constructed using the \Planck\ HFI data because the three lowest CO rotational transition lines at 115, 230 and 345 GHz significantly contribute to the signal of the 100, 217 and 353 GHz HFI channels respectively. Two different component separation methods are used to extract the CO maps from Planck HFI data. The maps obtained are then compared to one another and to existing external CO surveys. From these quality checks the best CO maps in terms of signal to noise and/or residual foreground contamination are selected. Three sets of velocity-integrated CO emission maps are produced: Type 1 maps of the CO (1-0), (2-1), and (3-2) rotational transitions with low foreground contamination but moderate signal-to-noise ratio; Type 2 maps for the (1-0) and (2-1) transitions with a better signal-to-noise ratio; and one Type 3 map, a line composite map with the best signal-to-noise ratio in order to locate the faintest molecular regions. The maps are described in detail. They are shown to be fully compatible with previous surveys of parts of the Galactic Plane and also of fainter regions out of the Galactic plane. The Planck HFI velocity-integrated CO maps for the (1-0), (2-1), and (3-2) rotational transitions provide an unprecedented all-sky CO view of the Galaxy. These maps are also of great interest to monitor potential CO contamination on CMB \Planck\ studies.

We present the Planck likelihood, a complete statistical description of the two-point correlation function of the CMB temperature fluctuations. We use this likelihood to derive the Planck CMB power spectrum over three decades in l, covering 2 <= l <= 2500. The main source of error at l <= 1500 is cosmic variance. Uncertainties in small-scale foreground modelling and instrumental noise dominate the error budget at higher l's. For l < 50, our likelihood exploits all Planck frequency channels from 30 to 353 GHz through a physically motivated Bayesian component separation technique. At l >= 50, we employ a correlated Gaussian likelihood approximation based on angular cross-spectra derived from the 100, 143 and 217 GHz channels. We validate our likelihood through an extensive suite of consistency tests, and assess the impact of residual foreground and instrumental uncertainties on cosmological parameters. We find good internal agreement among the high-l cross-spectra with residuals of a few uK^2 at l <= 1000. We compare our results with foreground-cleaned CMB maps, and with cross-spectra derived from the 70 GHz Planck map, and find broad agreement in terms of spectrum residuals and cosmological parameters. The best-fit LCDM cosmology is in excellent agreement with preliminary Planck polarisation spectra. The standard LCDM cosmology is well constrained by Planck by l <= 1500. For example, we report a 5.4 sigma deviation from n_s /= 1. Considering various extensions beyond the standard model, we find no indication of significant departures from the LCDM framework. Finally, we report a tension between the best-fit LCDM model and the low-l spectrum in the form of a power deficit of 5-10% at l <~ 40, significant at 2.5-3 sigma. We do not elaborate further on its cosmological implications, but note that this is our most puzzling finding in an otherwise remarkably consistent dataset. (Abridged)

In this paper we present forecasts of the contamination on different shapes of the primordial non-Gaussianity fnl parameter -- detectable on future Cosmic Microwave Background (CMB) high--resolution anisotropy maps -- produced by unresolved extragalactic point sources at frequencies of cosmological interest (45--375 GHz). We consider two scenarios: an ideal (noiseless) mission and a possible future space-borne satellite, with instrumental characteristics similar to the ones proposed for the Cosmic Origins Explorer (COrE). The local, equilateral, orthogonal and flat shapes are considered in both temperature (intensity) and polarized emission data. The angular power spectrum and bispectrum of extragalactic point sources are estimated by state-of-the-art models of source number counts. The impact of all the most relevant (far--IR and radio selected) source populations on these shapes at COrE frequencies is studied. The results of this analysis show that unresolved extragalactic point sources should not induce a very relevant non-Gaussian signal in the frequency range 100--200 GHz, thus not preventing a correct estimate of the CMB primordial fnl parameter. Polarization information allows one to significantly reduce the error--bars in the fnl parameter and the bias induced by unresolved sources and, hence, to widen the range of frequencies for fnl studies. On the contrary, at nu < 100 GHz or nu > 225 GHz, important non-Gaussian deviations in CMB anisotropy maps are expected due to unresolved extragalactic sources.

We present a map of the Cosmic Microwave Background (CMB) anisotropies induced by the late Integrated Sachs Wolfe effect. The map is constructed by combining the information of the WMAP 7-yr CMB data and the NRAO VLA Sky Survey (NVSS) through a linear filter. This combination improves the quality of the map that would be obtained using information only from the Large Scale Structure data. In order to apply the filter, a given cosmological model needs to be assumed. In particular, we consider the standard LCDM model. As a test of consistency, we show that the reconstructed map is in agreemet with the assumed model, which is also favoured against a scenario where no correlation between the CMB and NVSS catalogue is considered.

Taking advantage of the all-sky coverage and broad frequency range of the Planck satellite, we study the Sunyaev-Zeldovich (SZ) and pressure profiles of 62 nearby massive clusters detected at high significance in the 14-month nominal survey. Careful reconstruction of the SZ signal indicates that most clusters are individually detected at least out to R500. By stacking the radial profiles, we have statistically detected the radial SZ signal out to 3 x R500, i.e., at a density contrast of about 50-100, though the dispersion about the mean profile dominates the statistical errors across the whole radial range. Our measurement is fully consistent with previous Planck results on integrated SZ fluxes, further strengthening the agreement between SZ and X-ray measurements inside R500. Correcting for the effects of the Planck beam, we have calculated the corresponding pressure profiles. This new constraint from SZ measurements is consistent with the X-ray constraints from XMM-Newton in the region in which the profiles overlap (i.e., [0.1-1]R500), and is in fairly good agreement with theoretical predictions within the expected dispersion. At larger radii the average pressure profile is slightly flatter than most predictions from numerical simulations. Combining the SZ and X-ray observed profiles into a joint fit to a generalised pressure profile gives best-fit parameters [P0, c500, gamma, alpha, beta] = [6.41, 1.81, 0.31, 1.33, 4.13]. Using a reasonable hypothesis for the gas temperature in the cluster outskirts we reconstruct from our stacked pressure profile the gas mass fraction profile out to 3 x R500. Within the temperature driven uncertainties, our Planck constraints are compatible with the cosmic baryon fraction and expected gas fraction in halos.

In this work we present constraints on different shapes of primordial non-Gaussianity using the Wilkinson Microwave Anisotropy Probe (WMAP) 7-year data and the spherical Mexican hat wavelet fnl estimator including the linear term correction. In particular we focus on the local, equilateral and orthogonal shapes. We first analyse the main statistical properties of the wavelet estimator and show the conditions to reach optimality. We include the linear term correction in our estimators and compare the estimates with the values already published using only the cubic term. The estimators are tested with realistic WMAP simulations with anisotropic noise and the WMAP KQ75 sky cut. The inclusion of the linear term correction shows a negligible improvement (< 1 per cent) in the error-bar for any of the shapes considered. The results of this analysis show that, in the particular case of the wavelet estimator, the optimality for WMAP anisotropy levels is basically achieved with the mean subtraction and in practical terms there is no need of including a linear term once the mean has been subtracted. Our best estimates are now: local fnl = 39.0 +/ 21.4, equilateral fnl = -62.8 +/- 154.0 and orthogonal fnl = -159.8 +/- 115.1 (all cases 68 per cent CL). We have also computed the expected linear term correction for simulated Planck maps with anisotropic noise at 143 GHz following the Planck Sky Model and including a mask. The improvement achieved in this case for the local fnl error-bar is also negligible (0.4 per cent).

We use a Quadratic Maximum Likelihood (QML) method to estimate the angular power spectrum of the cross-correlation between cosmic microwave background and large scale structure maps as well as their individual auto-spectra. We describe our implementation of this method and demonstrate its accuracy on simulated maps. We apply this optimal estimator to WMAP 7-year and NRAO VLA Sky Survey (NVSS) data and explore the robustness of the angular power spectrum estimates obtained by the QML method. With the correction of the declination systematics in NVSS, we can safely use most of the information contained in this survey. We then make use of the angular power spectrum estimates obtained by the QML method to derive constraints on the dark energy critical density in a flat $\Lambda$CDM model by different likelihood prescriptions. When using just the cross-correlation between WMAP 7 year and NVSS maps with 1.8$^\circ$ resolution, the best-fit model has a cosmological constant of approximatively 70% of the total energy density, disfavouring an Einstein-de Sitter Universe at more than 2 $\sigma$ CL (confidence level).

The performance of the Planck instruments in space is enabled by their low operating temperatures, 20K for LFI and 0.1K for HFI, achieved through a combination of passive radiative cooling and three active mechanical coolers. The scientific requirement for very broad frequency coverage led to two detector technologies with widely different temperature and cooling needs. Active coolers could satisfy these needs; a helium cryostat, as used by previous cryogenic space missions (IRAS, COBE, ISO, Spitzer, AKARI), could not. Radiative cooling is provided by three V-groove radiators and a large telescope baffle. The active coolers are a hydrogen sorption cooler (<20K), a 4He Joule-Thomson cooler (4.7K), and a 3He-4He dilution cooler (1.4K and 0.1K). The flight system was at ambient temperature at launch and cooled in space to operating conditions. The HFI bolometer plate reached 93mK on 3 July 2009, 50 days after launch. The solar panel always faces the Sun, shadowing the rest of Planck, andoperates at a mean temperature of 384K. At the other end of the spacecraft, the telescope baffle operates at 42.3K and the telescope primary mirror operates at 35.9K. The temperatures of key parts of the instruments are stabilized by both active and passive methods. Temperature fluctuations are driven by changes in the distance from the Sun, sorption cooler cycling and fluctuations in gas-liquid flow, and fluctuations in cosmic ray flux on the dilution and bolometer plates. These fluctuations do not compromise the science data.

(abridged) Planck has observed the entire sky from 30 GHz to 857GHz. The observed foreground emission contains contributions from different phases of the interstellar medium (ISM). We have separated the observed Galactic emission into the different gaseous components (atomic, molecular and ionised) in each of a number of Galactocentric rings. Templates are created for various Galactocentric radii using velocity information from atomic (neutral hydrogen) and molecular (12CO) observations. The ionised template is assumed to be traced by free-free emission as observed by WMAP, while 408 MHz emission is used to trace the synchrotron component. Gas emission not traced by the above templates, namely "ark gas", as evidenced using Planck data, is included as an additional template, the first time such a component has been used in this way. These templates are then correlated with each of the Planck frequency bands, as well as other ancillary data. The emission per column density of the gas templates allows us to create distinct spectral energy distributions (SEDs) per Galactocentric ring and in each of the gaseous tracers from 1.4 GHz to 25 THz (12\mu m). Apart from the thermal dust and free-free emission, we have probed the Galaxy for anomalous (e.g., spinning) dust as well as synchrotron emission. We show that anomalous dust emission is present in the atomic, molecular and dark gas phases throughout the Galactic disk. The derived dust propeties associated with the dark gas phase are derived but do not allow us to reveal the nature of this phase. For all environments, the anomalous emission is consistent with rotation from polycyclic aromatic hydrocarbons (PAHs) and, according to our simple model, accounts for $(25\pm5)%$ (statistical) of the total emission at 30 GHz.

We present models to predict high frequency counts of extragalactic radio sources using physically grounded recipes to describe the complex spectral behaviour of blazars, that dominate the mm-wave counts at bright flux densities. We show that simple power-law spectra are ruled out by high-frequency (nu>100 GHz) data. These data also strongly constrain models featuring the spectral breaks predicted by classical physical models for the synchrotron emission produced in jets of blazars (Blandford & Konigl 1979; Konigl 1981). A model dealing with blazars as a single population is, at best, only marginally consistent with data coming from current surveys at high radio frequencies. Our most successful model assumes different distributions of break frequencies, nu_M, for BL Lacs and Flat-Spectrum Radio Quasars (FSRQs). The former objects have substantially higher values of nu_M, implying that the synchrotron emission comes from more compact regions; therefore, a substantial increase of the BL Lac fraction at high radio frequencies and at bright flux densities is predicted. Remarkably, our best model is able to give a very good fit to all the observed data on number counts and on distributions of spectral indices of extragalactic radio sources at frequencies above 5 and up to 220 GHz. Predictions for the forthcoming sub-mm blazar counts from Planck, at the highest HFI frequencies, and from Herschel surveys are also presented.

This paper presents the first results of comparison of Planck along with IRAS data with Green Bank Telescope 21-cm observations in 14 fields covering more than 800 square degrees at high Galactic latitude. Galactic dust emission for fields with average HI column density lower than 2 x 10^20 cm^-2 is well correlated with 21-cm emission. The residual emission in these fields, once the HI-correlated emission is removed, is consistent with the expected statistical properties of the cosmic infrared background fluctuations. Fields with larger column densities show significant excess dust emission compared to the HI column density. Regions of excess lie in organized structures that suggest the presence of hydrogen in molecular form, though they are not always correlated with CO emission. Dust emission from intermediate-velocity clouds is detected with high significance. Its spectral properties are consistent with, compared to the local ISM values, significantly hotter dust (T~20 K), lower sub-millimeter dust opacity, and a relative abundance of very small grains to large grains about four times higher. These results are compatible with expectations for clouds that are part of the Galactic fountain in which there is dust shattering and fragmentation. Correlated dust emission in HVCs is not detected; the average of the 99.9% confidence upper limits to the emissivity is 0.15 times the local ISM value at 857 and 3000 GHz, in accordance with gas phase evidence for lower metallicity and depletion in these clouds. Unexpected anti-correlated variations of the dust temperature and emission cross-section per H atom are identified in the local ISM and IVCs, a trend that continues into molecular environments. This suggests that dust growth through aggregation, seen in molecular clouds, is active much earlier in the cloud condensation and star formation processes.

We present constraints on the non-linear coupling parameter fnl with the Wilkinson Microwave Anisotropy Probe (WMAP) data. We use the method based on the spherical Mexican hat wavelet (SMHW) to measure the fnl parameter for three of the most interesting shapes of primordial non-Gaussianity: local, equilateral and orthogonal. Our results indicate that this parameter is compatible with a Gaussian distribution within the two sigma confidence level (CL) for the three shapes and the results are consistent with the values presented by the WMAP team. We have included in our analysis the impact on fnl due to contamination by unresolved point sources. The point sources add a positive contribution of Delta(fnl) = 2.5 \pm 3.0, Delta(fnl) = 37 \pm 18 and Delta(fnl) = 25 \pm 14 for the local, equilateral and orthogonal cases respectively. As mentioned by the WMAP team, the contribution of the point sources to the orthogonal and equilateral form is expected to be larger than the local one and thus it cannot be neglected in future constraints on these parameters. Taking into account this contamination, our best estimates for fnl are -16.0 \leq fnl \leq 76.0, -382 \leq fnl \leq 202 and -394 \leq fnl \leq 34 at 95% CL for the local, equilateral and orthogonal cases respectively. The three shapes are compatible with zero at 95% CL (2{\sigma}). Our conclusion is that the WMAP 7-year data are consistent with Gaussian primordial fluctuations within ~2{\sigma} CL. We stress however the importance of taking into account the unresolved point sources in the measurement of fnl in future works, especially when using more precise data sets such as the forthcoming Planck data.

COrE (Cosmic Origins Explorer) is a fourth-generation full-sky, microwave-band satellite recently proposed to ESA within Cosmic Vision 2015-2025. COrE will provide maps of the microwave sky in polarization and temperature in 15 frequency bands, ranging from 45 GHz to 795 GHz, with an angular resolution ranging from 23 arcmin (45 GHz) and 1.3 arcmin (795 GHz) and sensitivities roughly 10 to 30 times better than PLANCK (depending on the frequency channel). The COrE mission will lead to breakthrough science in a wide range of areas, ranging from primordial cosmology to galactic and extragalactic science. COrE is designed to detect the primordial gravitational waves generated during the epoch of cosmic inflation at more than $3\sigma $ for $r=(T/S)>=10^{-3}$. It will also measure the CMB gravitational lensing deflection power spectrum to the cosmic variance limit on all linear scales, allowing us to probe absolute neutrino masses better than laboratory experiments and down to plausible values suggested by the neutrino oscillation data. COrE will also search for primordial non-Gaussianity with significant improvements over Planck in its ability to constrain the shape (and amplitude) of non-Gaussianity. In the areas of galactic and extragalactic science, in its highest frequency channels COrE will provide maps of the galactic polarized dust emission allowing us to map the galactic magnetic field in areas of diffuse emission not otherwise accessible to probe the initial conditions for star formation. COrE will also map the galactic synchrotron emission thirty times better than PLANCK. This White Paper reviews the COrE science program, our simulations on foreground subtraction, and the proposed instrumental configuration.

We study the spherical Mexican hat wavelet (SMHW) as a detector of primordial non-Gaussianity of the local type on the Cosmic Microwave Background (CMB) anisotropies. For this purpose we define third order statistics based on the wavelet coefficient maps and the original map. We find the dependence of these statistics in terms of the non-linear coupling parameter fnl and the bispectrum of this type of non-Gaussianity. We compare the analytical values for these statistics with the results obtained with non-Gaussian simulations for an ideal full-sky CMB experiment without noise. We study the power of this method to detect fnl, i. e. the variance of this parameter, and compare it with the variance obtained from the primary bispectrum for the same experiment. Finally we apply our wavelet based estimator on WMAP-like maps with incomplete sky and inhomogeneous noise and compare with the optimal bispectrum estimator. The results show that the wavelet cubic statistics are as efficient as the bispectrum as optimal detectors of this type of primordial non-Gaussianity.

One of the most interesting explanations for the non-Gaussian Cold Spot (CS) detected in the WMAP data by Vielva et al. 2004, is that it arises from the interaction of the CMB radiation with a cosmic texture (Cruz et al. 2007b). In this case, a lack of polarization is expected in the region of the spot, as compared to the typical values associated to large fluctuations of a GIRF. In addition, other physical processes related to a non-linear evolution of the gravitational field could lead to a similar scenario. However, some of these alternative scenarios (e.g., a large void in the large scale structure) have been shown to be very unlikely. In this work we characterise the polarization properties of the Cold Spot under both hypotheses: a large Gaussian spot and an anomalous feature generated, for instance, by a cosmic texture. We propose a methodology to distinguish between them, and we discuss its discrimination power as a function of the instrumental noise level. In particular, we address the cases of current experiments, like WMAP and Planck, and others in development as QUIJOTE. We find that for an ideal experiment the Gaussian hypothesis could be rejected at a significance level better than 0.8%. While WMAP is far from providing useful information in this respect, we find that Planck will be able to reach a significance of around 7%; in addition, we show that the ground-based experiment QUIJOTE could provide a significance of around 1%. If these results are combined with the significance level found for the CS in temperature, the capability of QUIJOTE and Planck to reject the alternative hypothesis becomes 0.025% and 0.124%, respectively.

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.

The 30 and 44 GHz Back End Modules (BEM) for the Planck Low Frequency Instrument are broadband receivers (20% relative bandwidth) working at room temperature. The signals coming from the Front End Module are amplified, band pass filtered and finally converted to DC by a detector diode. Each receiver has two identical branches following the differential scheme of the Planck radiometers. The BEM design is based on MMIC Low Noise Amplifiers using GaAs P-HEMT devices, microstrip filters and Schottky diode detectors. Their manufacturing development has included elegant breadboard prototypes and finally qualification and flight model units. Electrical, mechanical and environmental tests were carried out for the characterization and verification of the manufactured BEMs. A description of the 30 and 44 GHz Back End Modules of Planck-LFI radiometers is given, with details of the tests done to determine their electrical and environmental performances. The electrical performances of the 30 and 44 GHz Back End Modules: frequency response, effective bandwidth, equivalent noise temperature, 1/f noise and linearity are presented.

We give the calibration and scientific performance parameters of the Planck Low Frequency Instrument (LFI) measured during the ground cryogenic test campaign. These parameters characterise the instrument response and constitute our best pre-launch knowledge of the LFI scientific performance. The LFI shows excellent $1/f$ stability and rejection of instrumental systematic effects; measured noise performance shows that LFI is the most sensitive instrument of its kind. The set of measured calibration parameters will be updated during flight operations through the end of the mission.

In this paper we present the Low Frequency Instrument (LFI), designed and developed as part of the Planck space mission, the ESA program dedicated to precision imaging of the cosmic microwave background (CMB). Planck-LFI will observe the full sky in intensity and polarisation in three frequency bands centred at 30, 44 and 70 GHz, while higher frequencies (100-850 GHz) will be covered by the HFI instrument. The LFI is an array of microwave radiometers based on state-of-the-art Indium Phosphide cryogenic HEMT amplifiers implemented in a differential system using blackbody loads as reference signals. The front-end is cooled to 20K for optimal sensitivity and the reference loads are cooled to 4K to minimise low frequency noise. We provide an overview of the LFI, discuss the leading scientific requirements and describe the design solutions adopted for the various hardware subsystems. The main drivers of the radiometric, optical and thermal design are discussed, including the stringent requirements on sensitivity, stability, and rejection of systematic effects. Further details on the key instrument units and the results of ground calibration are provided in a set of companion papers.

This paper provides an overview of the Low Frequency Instrument (LFI) programme within the ESA Planck mission. The LFI instrument has been developed to produce high precision maps of the microwave sky at frequencies in the range 27-77 GHz, below the peak of the cosmic microwave background (CMB) radiation spectrum. The scientific goals are described, ranging from fundamental cosmology to Galactic and extragalactic astrophysics. The instrument design and development are outlined, together with the model philosophy and testing strategy. The instrument is presented in the context of the Planck mission. The LFI approach to ground and inflight calibration is described. We also describe the LFI ground segment. We present the results of a number of tests demonstrating the capability of the LFI data processing centre (DPC) to properly reduce and analyse LFI flight data, from telemetry information to calibrated and cleaned time ordered data, sky maps at each frequency (in temperature and polarization), component emission maps (CMB and diffuse foregrounds), catalogs for various classes of sources (the Early Release Compact Source Catalogue and the Final Compact Source Catalogue). The organization of the LFI consortium is briefly presented as well as the role of the core team in data analysis and scientific exploitation. All tests carried out on the LFI flight model demonstrate the excellent performance of the instrument and its various subunits. The data analysis pipeline has been tested and its main steps verified. In the first three months after launch, the commissioning, calibration, performance, and verification phases will be completed, after which Planck will begin its operational life, in which LFI will have an integral part.

We present new constraints on the non-linear coupling parameter fnl with the Wilkinson Microwave Anisotropy Probe (WMAP) data. We use an updated method based on the spherical Mexican hat wavelet (SMHW) which provides improved constraints on the fnl parameter. This paper is a continuation of a previous work by Curto et al. where several third order statistics based on the SMHW were considered. In this paper, we use all the possible third order statistics computed from the wavelet coefficient maps evaluated at 12 angular scales. The scales are logarithmically distributed from 6.9 arcmin to 500 arcmin. Our analysis indicates that fnl is constrained to -18 < fnl < +80 at 95% confidence level (CL) for the combined V+W WMAP map. This value has been corrected by the presence of undetected point sources, which adds a positive contribution of Delta_fnl = 6 +- 5. Our result excludes at ~99% CL the best-fitting value fnl=87 reported by Yadav & Wandelt. We have also constrained fnl for the Q, V and W frequency bands separately, finding compatibility with zero at 95 % CL for the Q and V bands but not for the W band. We have performed some further tests to understand the cause of this deviation which indicate that systematics associated to the W radiometers could be responsible for this result. Finally we have performed a Galactic North-South analysis for fnl. We have not found any asymmetry, i.e. the best-fitting fnl for the northern pixels is compatible with the best-fitting fnl for the southern pixels.

The WMAP cold spot was found by applying spherical wavelets to the first year WMAP data. An excess of kurtosis of the wavelet coefficient was observed at angular scales of around 5 degrees. This excess was shown to be inconsistent with Gaussian simulations with a p-value of around 1%. A cold spot centered at (b = -57, l = 209) was shown to be the main cause of this deviation. Several hypotheses were raised to explain the origin of the cold spot. After performing a Bayesian template fit a collapsing cosmic texture was found to be the most probable hypothesis explaining the spot. Here we review the properties of the cold spot and the possible explanations.

We present a Gaussianity analysis of the WMAP 5-year Cosmic Microwave Background (CMB) temperature anisotropy data maps. We use several third order estimators based on the spherical Mexican hat wavelet. We impose constraints on the local non-linear coupling parameter fnl using well motivated non-Gaussian simulations. We analyse the WMAP maps at resolution of 6.9 arcmin for the Q, V, and W frequency bands. We use the KQ75 mask recommended by the WMAP team which masks out 28% of the sky. The wavelet coefficients are evaluated at 10 different scales from 6.9 to 150 arcmin. With these coefficients we compute the third order estimators which are used to perform a chi-squared analysis. The chi-squared statistic is used to test the Gaussianity of the WMAP data as well as to constrain the fnl parameter. Our results indicate that the WMAP data are compatible with the Gaussian simulations, and the fnl parameter is constrained to -8 < fnl < +111 at 95% CL for the combined V+W map. This value has been corrected for the presence of undetected point sources, which add a positive contribution of Delta_fnl = 3+-5 in the V+W map. Our results are very similar to those obtained by Komatsu et al (2008) using the bispectrum.

We present the current status of the QUIJOTE (Q-U-I JOint TEnerife) CMB Experiment, a new instrument which will start operations early 2009 at Teide Observatory, with the aim of characterizing the polarization of the CMB and other processes of galactic and extragalactic emission in the frequency range 10-30 GHz and at large angular scales. QUIJOTE will be a valuable complement at low frequencies for the PLANCK mission, and will have the required sensitivity to detect a primordial gravitational-wave component if the tensor-to-scalar ratio is larger than r=0.05.

The non-Gaussian cold spot found in the WMAP data has created controversy about its origin. Here we calculate the Bayesian posterior probability ratios for three different models that could explain the cold spot. A recent work claimed that the Spot could be caused by a cosmic texture, while other papers suggest that it could be due to the gravitational effect produced by an anomalously large void. Also the Sunyaev-Zeldovich effect caused by a cluster is taken into account as a possible origin. We perform a template fitting on a 20 degrees radius patch centered at Galactic coordinates (b = -57circ, l = 209) and calculate the posterior probability ratios for the void and Sunyaev-Zeldovich models, comparing the results to those obtained with textures. Taking realistic priors for the parameters, the texture interpretation is favored, while the void and Sunyaev-Zeldovich hypotheses are discarded. The temperature decrement produced by voids or clusters is negligible considering realistic values for the parameters.

The extraction of a signal from some observational data sets that contain different contaminant emissions, often at a greater level than the signal itself, is a common problem in Astrophysics and Cosmology. The signal can be recovered, for instance, using a simple Wiener filter. However, in certain cases, additional information may also be available, such as a second observation which correlates to a certain level with the sought signal. In order to improve the quality of the reconstruction, it would be useful to include as well this additional information. Under these circumstances, we have constructed a linear filter, the linear covariance-based filter, that extracts the signal from the data but takes also into account the correlation with the second observation. To illustrate the performance of the method, we present a simple application to reconstruct the so-called Integrated Sachs-Wolfe effect from simulated observations of the Cosmic Microwave Background and of catalogues of galaxies.

We present a Gaussianity analysis of the Archeops Cosmic Microwave Background (CMB) temperature anisotropy data maps at high resolution to constrain the non-linear coupling parameter fnl characterising well motivated non-Gaussian CMB models. We used the data collected by the most sensitive Archeops bolometer at 143 GHz. The instrumental noise was carefully characterised for this bolometer, and for another Archeops bolometer at 143 GHz used for comparison. Angular scales from 27 arcmin to 1.8 degrees and a large fraction of the sky, 21 %, covering both hemispheres (avoiding pixels with Galactic latitude |b| < 15 degrees) were considered. The three Minkowski functionals on the sphere evaluated at different thresholds were used to construct a chi-squared statistics for both the Gaussian and the non-Gaussian CMB models. The Archeops maps were produced with the Mirage optimal map-making procedure from processed time ordered data. The analysis is based on simulations of signal (Gaussian and non-Gaussian fnl CMB models) and noise which were processed in the time domain using the Archeops pipeline and projected on the sky using the Mirage optimal map-making procedure. The Archeops data were found to be compatible with Gaussianity after removal of highly noisy pixels at high resolution. The non-linear coupling parameter was constrained to -400 < fnl < 590 at 68% CL and -920 < fnl < 1075 at 95% CL, for realistic non-Gaussian CMB simulations.

In this chapter we review the present status of the Gaussianity studies of the CMB anisotropies, including physical effects producing non-Gaussianity, methods to test it and observational constraints.

We have estimated the CMB variance from the three-year WMAP data, finding a value which is significantly lower than the one expected from Gaussian simulations using the WMAP best-fit cosmological model, at a significance level of 98.7 per cent. This result is even more prominent if we consider only the north ecliptic hemisphere (99.8 per cent). Different analyses have been performed in order to identify a possible origin for this anomaly. In particular we have studied the behaviour of single radiometers and single year data as well as the effect of residual foregrounds and 1/f noise, finding that none of these possibilities can explain the low value of the variance. We have also tested the effect of varying the cosmological parameters, finding that the estimated CMB variance tends to favour higher values of $n_s$ than the one of the WMAP best-fit model. In addition, we have also tested the consistency between the estimated CMB variance and the actual measured CMB power spectrum of the WMAP data, finding a strong discrepancy. A possible interpretation of this result could be a deviation from Gaussianity and/or isotropy of the CMB.

The decomposition of a signal on the sphere with the steerable wavelet constructed from the second Gaussian derivative gives access to the orientation, signed-intensity, and elongation of the signal's local features. In the present work, the non-Gaussianity of the WMAP temperature data of the cosmic microwave background (CMB) is analyzed in terms of the first four moments of the statistically isotropic random fields associated with these local morphological measures, at wavelet scales corresponding to angular sizes between 27.5 arcminutes and 30 degrees on the celestial sphere. While no detection is made neither in the orientation analysis nor in the elongation analysis, a strong detection is made in the excess kurtosis of the signed-intensity of the WMAP data. The non-Gaussianity is observed with a significance level below 0.5% at a wavelet scale corresponding to an angular size around 10 degrees, and confirmed at neighbour scales. This supports a previous detection of an excess of kurtosis in the wavelet coefficient of the WMAP data with the axisymmetric Mexican hat wavelet (Vielva et al. 2004). Instrumental noise and foreground emissions are not likely to be at the origin of the excess of kurtosis. Large-scale modulations of the CMB related to some unknown systematics are rejected as possible origins of the detection. The observed non-Gaussianity may therefore probably be imputed to the CMB itself, thereby questioning the basic inflationary scenario upon which the present concordance cosmological model relies. Taking the CMB temperature angular power spectrum of the concordance cosmological model at face value, further analysis also suggests that this non-Gaussianity is not confined to the directions on the celestial sphere with an anomalous signed-intensity.

We have returned to our previous Bianchi VII_h analysis in light of the Cruz et al. 2007 suggestion that the cold spot observed near the southern Galactic pole may be a remnant temperature perturbation of a cosmic texture. In Bridges et al. 2006b we found two favoured left handed Bianchi VII_h templates with restricted prior probabilities so that the template was centred close to the cold spot. Using WMAP data `corrected' for the texture fit we have now reexamined both models to assess any changes to these conclusions. We find that both models are left almost entirely unconstrained by the data and consequently exhibit significantly reduced Bayesian evidences. Both models are now disfavoured by the data. This result reinforces our previous suggestion that the cold spot was likely to be driving any Bianchi VII_h detection.

The Cosmic Microwave Background provides our most ancient image of the Universe and our best tool for studying its early evolution. Theories of high energy physics predict the formation of various types of topological defects in the very early universe, including cosmic texture which would generate hot and cold spots in the Cosmic Microwave Background. We show through a Bayesian statistical analysis that the most prominent, 5 degree radius cold spot observed in all-sky images, which is otherwise hard to explain, is compatible with having being caused by a texture. From this model, we constrain the fundamental symmetry breaking energy scale to be phi_0 ~ 8.7 x 10^(15) GeV. If confirmed, this detection of a cosmic defect will probe physics at energies exceeding any conceivable terrestrial experiment.

Significant alignment and signed-intensity anomalies of local features of the cosmic microwave background (CMB) are detected on the three-year WMAP data, through a decomposition of the signal with steerable wavelets on the sphere. Firstly, an alignment analysis identifies two mean preferred planes in the sky, both with normal axes close to the CMB dipole axis. The first plane is defined by the directions toward which local CMB features are anomalously aligned. A mean preferred axis is also identified in this plane, located very close to the ecliptic poles axis. The second plane is defined by the directions anomalously avoided by local CMB features. This alignment anomaly provides further insight on recent results (Wiaux et al. 2006). Secondly, a signed-intensity analysis identifies three mean preferred directions in the southern galactic hemisphere with anomalously high or low temperature of local CMB features: a cold spot essentially identified with a known cold spot (Vielva et al. 2004), a second cold spot lying very close to the southern end of the CMB dipole axis, and a hot spot lying close to the southern end of the ecliptic poles axis. In both analyses, the anomalies are observed at wavelet scales corresponding to angular sizes around 10 degress on the celestial sphere, with global significance levels around 1%. Further investigation reveals that the alignment and signed-intensity anomalies are only very partially related. Instrumental noise, foreground emissions, as well as some form of other systematics, are strongly rejected as possible origins of the detections. An explanation might still be envisaged in terms of a global violation of the isotropy of the Universe, inducing an intrinsic statistical anisotropy of the CMB.

A Gaussianity analysis using a goodness-of-fit test and the Minkowski functionals on the sphere has been performed to study the measured Archeops Cosmic Microwave Background (CMB) temperature anisotropy data for a 143 GHz Archeops bolometer. We consider large angular scales, greater than 1.8 degrees, and a large fraction of the North Galactic hemisphere, around 16%, with a galactic latitude b > 15 degrees. The considered goodness-of-fit test, first proposed by Rayner & Best (1989), has been applied to the data after a signal-to-noise decomposition. The three Minkowski functionals on the sphere have been used to construct a chi-square statistic using different thresholds. The first method has been calibrated using simulations of Archeops data containing the CMB signal and instrumental noise in order to check its asymptotic convergence. Two kind of maps produced with two different map-making techniques (coaddition and Mirage) have been analysed. Archeops maps for both Mirage and coaddition map-making, have been found to be compatible with Gaussianity. From these results we can exclude a dust and atmospheric contamination larger than 7.8% (90% CL). Also the non-linear coupling parameter f_{nl} can be constrained to be -800 < f_{nl} < 1100 at the 95% CL and on angular scales of 1.8 degrees. For comparison, the same method has been applied to data from the NASA WMAP satellite in the same region of sky. The 1-year and 3-year releases have been used. Results are compatible with those obtained with Archeops, implying in particular an upper limit for f_{nl} on degree angular scales.

The cosmic microwave background (CMB) is a relic radiation of the Big Bang and as such it contains a wealth of cosmological information. Statistical analyses of the CMB, in conjunction with other cosmological observables, represent some of the most powerful techniques available to cosmologists for placing strong constraints on the cosmological parameters that describe the origin, content and evolution of the Universe. The last decade has witnessed the introduction of wavelet analyses in cosmology and, in particular, their application to the CMB. We review here spherical wavelet analyses of the CMB that test the standard cosmological concordance model. The assumption that the temperature anisotropies of the CMB are a realisation of a statistically isotropic Gaussian random field on the sphere is questioned. Deviations from both statistical isotropy and Gaussianity are detected in the reviewed works, suggesting more exotic cosmological models may be required to explain our Universe. We also review spherical wavelet analyses that independently provide evidence for dark energy, an exotic component of our Universe of which we know very little currently. The effectiveness of accounting correctly for the geometry of the sphere in the wavelet analysis of full-sky CMB data is demonstrated by the highly significant detections of physical processes and effects that are made in these reviewed works.

Using a directional spherical wavelet analysis we detect the integrated Sachs-Wolfe (ISW) effect, indicated by a positive correlation between the first-year Wilkinson Microwave Anisotropy Probe (WMAP) and NRAO VLA Sky Survey (NVSS) data. Detections are made using both a directional extension of the spherical Mexican hat wavelet and the spherical butterfly wavelet. We examine the possibility of foreground contamination and systematics in the WMAP data and conclude that these factors are not responsible for the signal that we detect. The wavelet analysis inherently enables us to localise on the sky those regions that contribute most strongly to the correlation. On removing these localised regions the correlation that we detect is reduced in significance, as expected, but it is not eliminated, suggesting that these regions are not the sole source of correlation between the data. This finding is consistent with predictions made using the ISW effect, where one would expect weak correlations over the entire sky. In a flat universe the detection of the ISW effect provides direct and independent evidence for dark energy. We use our detection to constrain dark energy parameters by deriving a theoretical prediction for the directional wavelet covariance statistic for a given cosmological model. Comparing these predictions with the data we place constraints on the equation-of-state parameter $w$ and the vacuum energy density $\Omega_\Lambda$. We also consider the case of a pure cosmological constant, i.e. $w=-1$. For this case we rule out a zero cosmological constant at greater than the 99.9% significance level. All parameter estimates that we obtain are consistent with the standand cosmological concordance model values.

The study of the Gaussianity of the cosmic microwave background (CMB) radiation is a key topic to understand the process of structure formation in the Universe. In this paper, we review a very useful tool to perform this type of analysis, the Rayner & Best smooth tests of goodness of fit. We describe how the method has been adapted for its application to imaging and interferometric observations of the CMB and comment on some recent and future applications of this technique to CMB data.

Most of the cosmological information extracted from the CMB has been obtained through the power spectrum, however there is much more to be learnt from the statistical distribution of the temperature random field. We review some recent developments in the study of the Cosmic Microwave Background (CMB) anisotropies and present a description of the novel tools developed to analyse the properties of the CMB anisotropies beyond the power spectrum.

The non-Gaussian cold spot detected in wavelet space in the WMAP 1-year data, is detected again in the coadded WMAP 3-year data at the same position (b = -57, l = 209) and size in the sky (around 10 degrees). The present analysis is based on several statistical methods: kurtosis, maximum absolute temperature, number of pixels below a given threshold, volume and Higher Criticism. All these methods detect deviations from Gaussianity in the 3--year data set at a slightly higher confidence level than in the WMAP 1-year data. These small differences are mainly due to the new foreground reduction technique and not to the reduction of the noise level, which is negligible at the scale of the spot. In order to avoid "a posteriori" analyses, we recalculate for the WMAP 3-year data the significance of the deviation in the kurtosis. The skewness and kurtosis tests were the first tests performed with wavelets for the WMAP data. We obtain that the probability of finding an at least as high deviation in Gaussian simulations is 1.85%. The frequency dependence of the spot is shown to be extremely flat. Galactic foreground emissions are not likely to be responsible for the detected deviation from Gaussianity.

This paper reviews the application of a novel methodology for analysing the isotropy of the universe by probing the alignment of local structures in the CMB. The strength of the proposed methodology relies on the steerable wavelet filtering of the CMB signal. One the one hand, the filter steerability renders the computation of the local orientation of the CMB features affordable in terms of computation time. On the other hand, the scale-space nature of the wavelet filtering allows to explore the alignment of the local structures at different scales, probing possible different phenomena. We present the WMAP first-year data analysis recently performed by the same authors (Wiaux et al.), where an extremely significant anisotropy was found. In particular, a preferred plane was detected, having a normal direction with a northern end position close to the northern end of the CMB dipole axis. In addition, a most preferred direction was found in that plane, with a northern end direction very close to the north ecliptic pole. This result synthesised for the first time previously reported anomalies identified in the direction of the dipole and the ecliptic poles axes. In a forthcoming paper (Vielva et al.), we have extended our analysis to the study of individual frequency maps finding first indications for discarding foregrounds as the origin of the anomaly. We have also tested that the preferred orientations are defined by structures homogeneously distributed in the sky, rather than from localised regions. We have also analysed the WMAP 3-year data, finding the same anomaly pattern, although at a slightly lower significance level.

We study the power of several scalar quantities constructed on the sphere (presented in Monteserin et al. 2005) to detect non-Gaussianity on the temperature distribution of the cosmic microwave background (CMB). The test has been performed using non-Gaussian CMB simulations with injected skewness or kurtosis generated through the Edgeworth expansion. We have also taken into account in the analysis the effect of anisotropic noise and the presence of a Galactic mask. We find that the best scalars to detect an excess of skewness in the simulations are the derivative of the gradient, the fractional isotropy, the Laplacian and the shape index. For the kurtosis case, the fractional anisotropy, the Laplacian and the determinant are the quantities that perform better.

Using a directional spherical wavelet analysis we detect the integrated Sachs-Wolfe (ISW) effect, indicated by a positive correlation between the Wilkinson Microwave Anisotropy Probe (WMAP) and NRAO VLA Sky Survey (NVSS) data, at the 3.9 sigma level. In a flat universe the detection of the ISW effect provides direct and independent evidence for dark energy. Moreover, we use our detection to constrain the dark energy density Omega_Lambda. We obtain estimates for Omega_Lambda consistent with other analysis techniques and data sets and rule out a zero cosmological constant at greater than 99% significance.

(Abridged) We have used the Rayner & Best (1989) smooth tests of goodness-of-fit to study the Gaussianity of the Very Small Array (VSA) data. Out of the 41 published VSA individual pointings dedicated to cosmological observations, 37 are found to be consistent with Gaussianity, whereas four pointings show deviations from Gaussianity. In two of them, these deviations can be explained as residual systematic effects of a few visibility points which, when corrected, have a negligible impact on the angular power spectrum. The non-Gaussianity found in the other two (adjacent) pointings seems to be associated to a local deviation of the power spectrum of these fields with respect to the common power spectrum of the complete data set, at angular scales of the third acoustic peak (l = 700-900). No evidence of residual systematics is found in this case, and unsubstracted point sources are not a plausible explanation either. If those visibilities are removed, a cosmological analysis based on this new VSA power spectrum alone shows no differences in the parameter constraints with respect to our published results, except for the physical baryon density, which decreases by 10 percent. Finally, the method has been also used to analyse the VSA observations in the Corona Borealis supercluster region (Genova-Santos et al. 2005), which show a strong decrement which cannot be explained as primordial CMB. Our method finds a clear deviation (99.82%) with respect to Gaussianity in the second-order moment of the distribution, and which can not be explained as systematic effects. A detailed study shows that the non-Gaussianity is produced in scales of l~500, and that this deviation is intrinsic to the data (in the sense that can not be explained in terms of a Gaussian field with a different power spectrum).

We present a study of the effect of component separation on the recovered cosmic microwave background (CMB) temperature distribution, considering Gaussian and non-Gaussian input CMB maps. First, we extract the CMB component from simulated Planck data (in small patches of the sky) using the maximum entropy method (MEM), Wiener filter (WF) and a method based on the subtraction of foreground templates plus a linear combination of frequency channels (LCFC). We then apply a wavelet-based method to study the Gaussianity of the recovered CMB and compare it with the same analysis for the input map. When the original CMB map is Gaussian (and assuming that point sources have been removed), we find that none of the methods introduce non-Gaussianity (NG) in the CMB reconstruction. On the contrary, if the input CMB map is non-Gaussian, all the studied methods produce a reconstructed CMB with lower detections of NG than the original map. This effect is mainly due to the presence of instrumental noise. In this case, MEM tends to produce slightly higher non-Gaussian detections in the reconstructed map than WF whereas the detections are lower for the LCFC. We have also studied the effect of point sources in the MEM reconstruction. If no attempt to remove point sources is performed, they clearly contaminate the CMB reconstruction, introducing spurious NG. When the brightest point sources are removed from the data using the Mexican Hat Wavelet, the Gaussian character of the CMB is preserved. However, when analysing larger regions of the sky, the variance of our estimators will be appreciably reduced and, in this case, we expect the point source residuals to introduce spurious NG in the CMB. Thus, a careful subtraction (or masking) of point source emission is crucial when studying the Gaussianity of the CMB.

We question the global universe isotropy by probing the alignment of local structures in the cosmic microwave background (CMB) radiation. The original method proposed relies on a steerable wavelet decomposition of the CMB signal on the sphere. The analysis of the first-year WMAP data identifies a mean preferred plane with a normal direction close to the CMB dipole axis, and a mean preferred direction in this plane, very close to the ecliptic poles axis. Previous statistical anisotropy results are thereby synthesized, but further analyses are still required to establish their origin.

The non--Gaussian cold spot in the 1-year WMAP data, described in Vielva et al. and Cruz et al., is analysed in detail in the present paper. First of all, we perform a more rigorous calculation of the significance of the non-zero kurtosis detected in WMAP maps by Vielva et al. in wavelet space, mainly generated by the Spot. We confirm the robustness of that detection, since the probability of obtaining this deviation by chance is 0.69%. Afterwards, the morphology of the Spot is studied by applying Spherical Mexican Hat Wavelets with different ellipticities. The shape of the Spot is found to be almost circular. Finally, we discuss if the observed non-Gaussianity in wavelet space can arise from bad subtracted foreground residues in the WMAP maps. We show that the flat frequency dependence of the Spot cannot be explained by a thermal Sunyaev-Zeldovich effect. Based on our present knowledge of Galactic foreground emissions, we conclude that the significance of our detection is not affected by Galactic residues in the region of the Spot. Considering different Galactic foreground estimates, the probability of finding such a big cold spot in Gaussian simulations is always below 1%.

A method to compute several scalar quantities of Cosmic Microwave Background maps on the sphere is presented. We consider here four type of scalars: the Hessian matrix scalars, the distortion scalars, the gradient related scalars and the curvature scalars. Such quantities are obtained directly from the spherical harmonic coefficients (alm) of the map. We also study the probability density function of these quantities for the case of a homogeneous and isotropic Gaussian field, which are functions of the power spectrum of the initial field. From these scalars it is posible to construct a new set of scalars which are independent of the power spectrum of the field. We test our results using simulations and find a good agreement between the theoretical probability density functions and those obtained from simulations. Therefore, these quantities are proposed to investigate the presence of non-Gaussian features in CMB maps. Finally, we show how to compute the scalars in presence of anisotropic noise and realistic masks.

We report the first detection of the ISW effect in wavelet space, at scales in the sky around 7 degrees with a significance of around 3.3 sigma, by cross-correlating the WMAP first-year data and the NRAO VLA Sky Survey (NVSS). In addition, we present a detailed comparison among the capabilities of three different techniques for two different objectives: to detect the ISW and to put constraints in the nature of the dark energy. The three studied techniques are: the cross-angular power spectrum (CAPS, harmonic space), the correlation function (CCF, real space) and the covariance of the Spherical Mexican Hat Wavelet (SMHW) coefficients (CSMHW, wavelet space). We prove that the CSMHW is expected to provide a higher detection of the ISW effect for a certain scale. This prediction has been corroborated by the analysis of the data. The SMHW analysis shows that the cross-correlation signal is caused neither by systematic effects nor foreground contamination. However, by taking into account the information encoded in all the multipoles/scales/angles, the CAPS provides slightly better constraints than the SMHW in the cosmological parameters that define the nature of the dark energy. The limits provided by the CCF are wider than for the other two methods. Two different cases have been studied: 1) a flat Lambda-CDM universe and 2) a flat universe with an equation of state parameter different from -1. In the first case, the CAPS provides (for a bias value of b = 1.6) 0.59 < Lambda density < 0.84 (at 1 sigma CL). Moreover, the CAPS rejects the range Lambda density < 0.1 at 3.5 sigma, which is the highest detection of the dark energy reported up to date. In the second case, the CAPS gives 0.50 < dark energy density < 0.82 and -1.16 < w < 0.43 (at 1 sigma CL).

The standard inflationary model presents a simple scenario within which the homogeneity, isotropy and flatness of the universe appear as natural outcomes and, in addition, fluctuations in the energy density are originated during the inflationary phase. These seminal density fluctuations give rise to fluctuations in the temperature of the Cosmic Microwave Background (CMB) at the decoupling surface. Afterward, the CMB photons propagate almost freely, with slight gravitational interactions with the evolving gravitational field present in the large scale structure (LSS) of the matter distribution and a low scattering rate with free electrons after the universe becomes reionized. These secondary effects slightly change the shape of the intensity and polarization angular power spectra (APS) of the radiation. The APS contain very valuable information on the parameters characterizing the background model of the universe and those parametrising the power spectra of both matter density perturbations and gravitational waves. In the last few years data from sensitive experiments have allowed a good determination of the shape of the APS, providing for the first time a model of the universe very close to spatially flat. In particular the WMAP first year data, together with other CMB data at higher resolution and other cosmological data sets, have made possible to determine the cosmological parameters with a precision of a few percent. The most striking aspect of the derived model of the universe is the unknown nature of most of its energy contents. This and other open problems in cosmology represent exciting challenges for the CMB community. The future ESA Planck mission will undoubtely shed some light on these remaining questions.

We investigate which practical constraints are imposed by foregrounds to the detection of the B-mode polarization generated by gravitational waves in the case of experiments of the type currently being planned. Because the B-mode signal is probably dominated by foregrounds at all frequencies, the detection of the cosmological component depends drastically on our ability for removing foregrounds. We provide an analytical expression to estimate the level of the residual polarization for Galactic foregrounds, according to the method employed for their subtraction. We interpret this result in terms of the lower limit of the tensor-to-scalar ratio r that allows to disentangle the cosmological B-mode polarization from the foregrounds contribution. Polarized emission from extragalactic radio sources and gravitational lensing is also taken into account. As a first approach, we consider the ideal limit of an instrumental noise--free experiment: for a full--sky coverage and a degree resolution, we obtain a limit of r~10^(-4). This value can be improved by high--resolution experiments and, in principle, no clear fundamental limit on the detectability of gravitational waves polarization is found. Our analysis is also applied to planned or hypothetical future polarization experiments, taking into account expected noise levels.

We adapt the smooth tests of goodness of fit developed by Rayner and Best to the study of the non-Gaussianity of interferometric observations of the cosmic microwave background (CMB). The interferometric measurements (visibilities) are transformed into signal-to-noise eigenmodes, and then the method is applied directly in Fourier space. This transformation allows us to perform the analysis in different subsets of eigenmodes according to their signal-to-noise level. The method can also deal with non-uniform or incomplete coverage of the UV plane. We explore here two possibilities: we analyse either the real and imaginary parts of the complex visibilities (Gaussianly distributed under the Gaussianity hypothesis) or their phases (uniformly distributed under the Gaussianity hypothesis). The power of the method in discriminating between Gaussian and non-Gaussian distributions is studied by using several kinds of non-Gaussian simulations. On the one hand, we introduce a certain degree of non-Gaussianity directly into the Fourier space using the Edgeworth expansion, and afterwards the desired correlation is introduced. On the other hand, we consider interferometric observations of a map with topological defects (cosmic strings). To these previous non-Gaussian simulations we add different noise levels and quantify the required signal-to-noise ratio necessary to achieve a detection of these non-Gaussian features. Finally, we have also studied the ability of the method to constrain the so-called nonlinear coupling constant f_{NL} using chi^2 simulations.The whole method is illustrated here by application to simulated data from the Very Small Array interferometer.

An extremely cold and big spot in the WMAP 1-year data is analyzed. Our work is a continuation of a previous paper (Vielva et al. 2004) where non-Gaussianity was detected, with a method based on the Spherical Mexican Hat Wavelet (SMHW) technique. We study the spots at different thresholds on the SMHW coefficient maps, considering six estimators, namely number of maxima, number of minima, number of hot and cold spots, and number of pixels of the spots. At SMHW scales around 4 degrees (10 degrees on the sky), the data deviate from Gaussianity. The analysis is performed on all sky, the northern and southern hemispheres, and on four regions covering all the sky. A cold spot at (b = -57, l = 209) is found to be the source of this non-Gaussian signature. We compare the spots of our data with 10000 Gaussian simulations, and conclude that only around 0.2% of them present such a cold spot. Excluding this spot, the remaining map is compatible with Gaussianity and even the excess of kurtosis in Vielva et al. 2004, is found to be due exclusively to this spot. Finally, we study whether the spot causing the observed deviation from Gaussianity could be generated by systematics or foregrounds. None of them seem to be responsible for the non-Gaussian detection.

A non-Gaussian detection in the WMAP 1-year data is reported. The detection has been found in the combined Q-V-W map proposed by the WMAP team (Komatsu et al. 2003) after applying a wavelet technique based on the Spherical Mexican Hat Wavelet (SMHW). The skewness and the kurtosis of the SMHW coefficients are calculated at different scales. A non-Gaussian signal is detected at scales of the SMHW around 4 deg (size in the sky of around 10 deg). The right tail probability of the detection is approx. 0.4%. In addition, a study of Gaussianity is performed in each hemisphere. The northern hemisphere is compatible with Gaussianity, whereas the southern one deviates from Gaussianity with a right tail probability of approx. 0.1%. Systematics, foregrounds and uncertainties in the estimation of the cosmological parameters are carefully studied in order to identify the possible source of non-Gaussianity. The detected deviation from Gaussianity is not found to be caused by systematic effects: 1) each one of the Q, V and W receivers shows the same non-Gaussianity pattern, and 2) several combinations of the different receivers at each frequency band do not show this non-Gaussian pattern. Similarly, galactic foregrounds show a negligible contribution to the non-Gaussian detection: non-Gaussianity is detected in all the WMAP maps and no frequency dependence is observed. Moreover, the expected foreground contribution to the combined WMAP map was added to CMB Gaussian simulations showing a behaviour compatible with the Gaussian model. Influence of uncertainties in the CMB power spectrum estimation are also quantified. Hence, possible intrinsic temperature fluctuations (like secondary anisotropies and primordial features) can not be rejected as the source of this non-Gaussian detection.

We present a method to simulate the polarization properties of extragalactic radio sources at microwave frequencies. Polarization measurements of nearly 2x10^6 sources at 1.4 GHz are provided by the NVSS survey. Using this catalogue and the GB6 survey, we study the distribution of the polarization degree of both steep- and flat-spectrum sources. We find that the polarization degree is anti-correlated with the flux density for the former population, while no correlation is detected for the latter. The available high-frequency data are exploited to determine the frequency dependence of the distribution of polarization degrees. Using such information and the evolutionary model by Toffolatti et al. (1998), we estimate the polarization power spectrum of extragalactic radio sources at \geq 30 GHz and their contamination of CMB polarization maps. Two distinct methods to compute point-source polarization spectra are presented, extending and improving the one generally used in previous analyses. While extragalactic radio sources can significantly contaminate the CMB E-mode power spectrum only at low frequencies (<30 GHz), they can severely constrain the detectability of the CMB B-mode up to ~100 GHz.

We review two powerful methods to test the Gaussianity of the cosmic microwave background (CMB): one based on the distribution of spherical wavelet coefficients and the other on smooth tests of goodness-of-fit. The spherical wavelet families proposed to analyse the CMB are the Haar and the Mexican Hat ones. The latter is preferred for detecting non-Gaussian homogeneous and isotropic primordial models containing some amount of skewness or kurtosis. Smooth tests of goodness-of-fit have recently been introduced in the field showing some interesting properties. We will discuss the smooth tests of goodness-of-fit developed by Rayner and Best for the univariate as well as for the multivariate analysis.

In this work, goodness-of-fit tests are adapted and applied to CMB maps to detect possible non-Gaussianity. We use Shapiro-Francia test and two Smooth goodness-of-fit tests: one developed by Rayner and Best and another one developed by Thomas and Pierce. The Smooth tests test small and smooth deviations of a prefixed probability function (in our case this is the univariate Gaussian). Also, the Rayner and Best test informs us of the kind of non-Gaussianity we have: excess of skewness, of kurtosis, and so on. These tests are optimal when the data are independent. We simulate and analyse non-Gaussian signals in order to study the power of these tests. These non-Gaussian simulations are constructed using the Edgeworth expansion, and assuming pixel-to-pixel independence. As an application, we test the Gaussianity of the MAXIMA data. Results indicate that the MAXIMA data are compatible with Gaussianity. Finally, the values of the skewness and kurtosis of MAXIMA data are constrained by |S| \le 0.035 and |K| \le 0.036 at the 99% confidence level.

The LFI (Low Frequency Instrument) on board the ESA Planck satellite is constituted by an array of radiometric detectors actively cooled at 20 K in the 30-70 GHz frequency range in the focal plane of the Planck telescope. In this paper we present an overview of the LFI instrument, with a particular focus on the radiometer design. The adopted pseudo-correlation scheme uses a software balancing technique (with a tunable parameter called gain modulation factor) which is effective in reducing the radiometer susceptibility to amplifier instabilities also in presence of small non-idealities in the radiometric chain components, provided that the gain modulation factor is estimated with an accuracy of the order of 0.2%. These results have been recently confirmed by experimental laboratory measurements conducted on the LFI prototype radiometers at 30, 70 and 100 GHz.

We apply the iterative Expectation-Maximization algorithm (EM) to estimate the power spectrum of the CMB from multifrequency microwave maps. In addition, we are also able to provide a reconstruction of the CMB map. By assuming that the combined emission of the foregrounds plus the instrumental noise is Gaussian distributed in Fourier space, we have simplified the EM procedure finding an analytical expression for the maximization step. By using the simplified expression the CPU time can be greatly reduced. We test the stability of our power spectrum estimator with realistic simulations of Planck data, including point sources and allowing for spatial variation of the frequency dependence of the Galactic emissions. Without prior information about any of the components, our new estimator can recover the CMB power spectrum up to scales l~1500 with less than 10% error. This result is significantly improved if the brightest point sources are removed before applying our estimator. In this way, the CMB power spectrum can be recovered up to l~1700 with 10% error and up to l~2100 with 50% error. This result is very close to the one that would be obtained in the ideal case of only CMB plus white noise, for which all our assumptions are satisfied. Moreover, the EM algorithm also provides an straightforward mechanism to reconstruct the CMB map. The recovered cosmological signal shows a high degree of correlation (r = 0.98) with the input map and low residuals.

Goodness-of-Fit tests, including Smooth ones, are introduced and applied to detect non-Gaussianity in Cosmic Microwave Background simulations. We study the power of three different tests: the Shapiro-Francia test (1972), the uncategorised smooth test developed by Rayner and Best(1990) and the Neyman's Smooth Goodness-of-fit test for composite hypotheses (Thomas and Pierce 1979). The Smooth Goodness-of-Fit tests are designed to be sensitive to the presence of ``smooth'' deviations from a given distribution. We study the power of these tests based on the discrimination between Gaussian and non-Gaussian simulations. Non-Gaussian cases are simulated using the Edgeworth expansion and assuming pixel-to-pixel independence. Results show these tests behave similarly and are more powerful than tests directly based on cumulants of order 3, 4, 5 and 6. We have applied these tests to the released MAXIMA data. The applied tests are built to be powerful against detecting deviations from univariate Gaussianity. The Cholesky matrix corresponding to signal (based on an assumed cosmological model) plus noise is used to decorrelate the observations previous to the analysis. Results indicate that the MAXIMA data are compatible with Gaussianity.

We present an estimation of the point source (PS) catalogue that could be extracted from the forthcoming ESA Planck mission data. We have applied the Spherical Mexican Hat Wavelet (SMHW) to simulated all-sky maps that include CMB, Galactic emission (thermal dust, free-free and synchrotron), thermal Sunyaev-Zel'dovich effect and PS emission, as well as instrumental white noise. This work is an extension of the one presented in Vielva et al. (2001a). We have developed an algorithm focused on a fast local optimal scale determination, that is crucial to achieve a PS catalogue with a large number of detections and a low flux limit. An important effort has been also done to reduce the CPU time processor for spherical harmonic transformation, in order to perform the PS detection in a reasonable time. The presented algorithm is able to provide a PS catalogue above fluxes: 0.48 Jy (857 GHz), 0.49 Jy (545 GHz), 0.18 Jy (353 GHz), 0.12 Jy (217 GHz), 0.13 Jy (143 GHz), 0.16 Jy (100 GHz HFI), 0.19 Jy (100 GHz LFI), 0.24 Jy (70 GHz), 0.25 Jy (44 GHz) and 0.23 Jy (30 GHz). We detect around 27700 PS at the highest frequency Planck channel and 2900 at the 30 GHz one. The completeness level are: 70% (857 GHz), 75% (545 GHz), 70% (353 GHz), 80% (217 GHz), 90% (143 GHz), 85% (100 GHz HFI), 80% (100 GHz LFI), 80% (70 GHz), 85% (44 GHz) and 80% (30 GHz). In addition, we can find several PS at different channels, allowing the study of the spectral behaviour and the physical processes acting on them. We also present the basic procedure to apply the method in maps convolved with asymmetric beams. The algorithm takes ~ 72 hours for the most CPU time demanding channel (857 GHz) in a Compaq HPC320 (Alpha EV68 1 GHz processor) and requires 4 GB of RAM memory.

This paper considers filters (the Mexican hat wavelet, the matched and the scale-adaptive filters) that optimize the detection/separation of point sources on a background. We make a one-dimensional treatment, we assume that the sources have a Gaussian profile, i. e. $\tau (x) = e^{- x^2/2R^2}$, and a background modelled by an homogeneous and isotropic Gaussian random field, characterised by a power spectrum $P(q)\propto q^{-\gamma}, \gamma \geq 0$. Local peak detection is used after filtering. Then, the Neyman-Pearson criterion is used to define the confidence level for detections and a comparison of filters is done based on the number of spurious and true detections. We have performed numerical simulations to test theoretical ideas and conclude that the results of the simulations agree with the analytical results.

Nonlinearity introduced in slow-roll inflation will produce weakly non-Gaussian CMB temperature fluctuations. We have simulated non-Gaussian large scale CMB maps (including COBE-DMR constraints) introducing an additional quadratic term in the gravitational potential.The amount of nonlinearity being controlled by the so called nonlinear coupling parameter f_{nl}. An analysis based on the Spherical Mexican Hat wavelet was applied to these and to the COBE-DMR maps. Skewness values obtained at several scales were combined into a Fisher discriminant. Comparison of the Fisher discriminant distributions obtained for different nonlinear coupling parameters with the COBE-DMR values, sets a constraint of abs(f_{nl})<1100 at the 68% confidence level. This new constraint being tighter than the one previously obtained by using the bispectrum by Komatsu et al. (2002).

We show how one may expect a significant number of SZ detections in future Planck data without any of the typical assumptions needed in present component separation methods, such as about the power spectrum or the frequency dependence of any of the components, circular symmetry or a typical scale for the clusters. We reduce the background by subtracting an estimate of the point sources, dust and CMB. The final SZE map is estimated in Fourier space. The catalogue of returned clusters is complete above flux 200 mJy (353 GHz) while the lowest flux reached by our method is about 70 mJy (353 GHz). We predict a large number of detections (about 9000) in 4/5 of the sky. This large number of SZ detections will allow a robust and consistent analysis of the evolution of the cluster population with redshift and will have important implications for determining the best cosmological model.

We investigate the performance of spherical wavelets in discriminating between standard inflationary models (Gaussian) and non-Gaussian models. For the later we consider small perturbations of the Gaussian model in which an artificially specified skewness or kurtosis is introduced through the Edgeworth expansion. By combining all the information present in all the wavelet scales with the Fisher discriminant, we find that the spherical Mexican Hat wavelets are clearly superior to the spherical Haar wavelets. The former can detect levels of the skewness and kurtosis of ~1% for 33' resolution, an order of magnitude smaller than the later. Also, as expected, both wavelets are better for discriminating between the models than the direct consideration of moments of the temperature maps. The introduction of instrumental white noise in the maps, S/N=1, does not change the main results of this paper.

The CMB Planck satellite (ESA) will provide a catalogue of nearly 30000 galaxy clusters (Sunyaev-Zel'dovich effect) between redshift 0 and approx. 2-3 (depending on the cosmological model). This hughe catalogue of clusters will allow detailed studies of the evolution of the cluster population with redshift which will constrain the cosmological model. The identification of the high redshift clusters should be done with a large telescope (like GTC) since Planck data alone will not provide any clue about the redshift of the clusters. In this work we show that an optical follow up of only 300 clusters (randomly selected from the Planck catalogue) is needed in order to distinguish models with or without a cosmological constant.

The problem of detecting Sunyaev-Zel'dovich (SZ) clusters in multifrequency CMB observations is investigated using a number of filtering techniques. A multifilter approach is introduced, which optimizes the detection of SZ clusters on microwave maps. An alternative method is also investigated, in which maps at different frequencies are combined in an optimal manner so that existing filtering techniques can be applied to the single combined map. The SZ profiles are approximated by the circularly-symmetric template $\tau (x) = [1 +(x/r_c)^2]^{-\lambda}$, with $\lambda \simeq \tfrac{1}{2}$ and $x\equiv |\vec{x}|$, where the core radius $r_c$ and the overall amplitude of the effect are not fixed a priori, but are determined from the data. The background emission is modelled by a homogeneous and isotropic random field, characterized by a cross-power spectrum $P_{\nu_1 \nu_2}(q)$ with $q\equiv |\vec{q}|$. The filtering methods are illustrated by application to simulated Planck observations of a $12.8^\circ \times 12.8^\circ$ patch of sky in 10 frequency channels. Our simulations suggest that the Planck instrument should detect $\approx 10000$ SZ clusters in 2/3 of the sky. Moreover, we find the catalogue to be complete for fluxes $S > 170$ mJy at 300 GHz.

This paper presents scale-adaptive filters that optimize the detection/separation of compact sources on a background. We assume that the sources have a multiquadric profile, i. e. $\tau (x) = {[1 + {(x/r_c)}^2]}^{-\lambda}, \lambda \geq {1/2}, x\equiv |\vec{x}|$, and a background modeled by a homogeneous and isotropic random field, characterized by a power spectrum $P(q)\propto q^{-\gamma}, \gamma \geq 0, q\equiv |\vec{q}|$. We make an n-dimensional treatment but consider two interesting astrophysical applications related to clusters of galaxies (Sunyaev-Zel'dovich effect and X-ray emission).

Wavelet-related techniques have proven useful in the processing and analysis of one and two dimensional data sets (spectra in the former case, images in the latter). In this work we apply adaptive filters, introduced in a previous work (Sanz et al. 2001), to optimize the detection and extraction of point sources from a one-dimensional array of time-ordered data such as the one that will be produced by the future 30 GHz LFI28 channel of the ESA Planck mission. At a $4\sigma$ detection level 224 sources over a flux of 0.88 Jy are detected with a mean relative error (in absolute value) of 21% and a systematic bias of -7.7%. The position of the sources in the sky is determined with errors inferior to the size of the pixel. The catalogue of detected sources is complete at fluxes $\geq$ 4.3 Jy. The number of spurious detections is less than a 10% of the true detections. We compared the results with the ones obtained by filtering with a Gaussian filter and a Mexican Hat Wavelet of width equal to the scale of the sources. The adaptive filter outperforms the other filters in all considered cases. We conclude that optimal adaptive filters are well suited to detect and extract sources with a given profile embedded in a background of known statistical properties. In the Planck case, they could be useful to obtain a real-time preliminary catalogue of extragalactic sources, which would have a great scientific interest, e. g. for follow-up observations.

This paper introduces the use of pseudo-filters that optimize the detection/extraction of sources on a background. We assume as a first approach that such sources are described by a spherical (central) profile and that the background is represented by a homogeneous and isotropic random field. We make an n-dimensional treatment making emphasis in astrophysical applications for spectra, images and volumes, for the cases of exponential and Gaussian source profiles and scale-free power spectra to represent the background.

We present a maximum-entropy method (MEM) and `Mexican Hat' wavelet (MHW) joint analysis to recover the different components of the microwave sky from simulated observations by the ESA Planck satellite in a small patch of the sky. This combined method allows one to improve the CMB, Sunyaev-Zel'dovich and Galactic foregrounds separation achieved by the MEM technique alone. In particular, the reconstructed CMB map is free from any bright point source contamination. The joint analysis also produces point source catalogues at each Planck frequency which are more complete and accurate than those obtained by each method on its own. The results are especially improved at high frequencies where infrared galaxies dominate the point source contribution. Although this joint technique has been performed on simulated Planck data, it could be easily applied to other multifrequency CMB experiments, such as the forthcoming NASA MAP satellite or the recently performed Boomerang and MAXIMA experiments.

The spherical Mexican Hat wavelet is introduced in this paper, with the aim of testing the Gaussianity of the Cosmic Microwave Background temperature fluctuations. Using the information given by the wavelet coefficients at several scales, we have performed several statistical tests on the COBE-DMR maps to search for evidence of non-Gaussianity. Skewness, kurtosis, scale-scale correlations (for two and three scales) and Kolmogorov-Smirnov tests indicate that the COBE-DMR data are consistent with a Gaussian distribution. We have extended the analysis to compare temperature values provided by COBE-DMR data with distributions (obtained from Gaussian simulations) at each pixel and at each scale. The number of pixels with temperature values outside the 95% and the 99% is consistent with that obtained from Gaussian simulations, at all scales. Moreover, the extrema values for COBE-DMR data (maximum and minimum temperatures in the map) are also consistent with those obtained from Gaussian simulations.

In this work we show how well can the cosmological parameters be constrained using galaxy cluster data. We also show how in the process of fitting the model which attempts to describe the data it is possible to get some information about the cluster scaling relations T-M and L_x-T. Among other conclusions we found that only low density universes (Omega approx. 0.3 with or without cosmological constant) are compatible with recent data sets. These constraints will be improved with future SZ data.

We have used recent X-ray and optical data in order to impose some constraints on the cosmology and cluster scaling relations. Generically two kind of hypotheses define our model. First we consider that the cluster population is well described by the standard Press-Schechter (PS) formalism, and second, these clusters are supposed to follow scaling relations with mass: Temperature-Mass (T-M) and X-ray Luminosity-Mass (L_x - M). As a difference with many other authors we do not assume specific scaling relations to model cluster properties such as the usual $T-M$ virial relation or one observational determination of the $L_x-T$ relation. Instead we consider general free parameter scaling relations. With the previous model (PS plus scalings) we fit our free parameters to several X-ray and optical data with the advantage over many other works that we consider all the data sets at the same time. This prevents us from being inconsistent with some of the available observations. Among other interesting conclusions, we find that only low-density universes are compatible with all the data considered and that the degeneracy between $\Omega_m$ and $\sigma_8$ is broken. Also we obtain interesting limits on the parameters characterizing the scaling relations.

An estimation of the number and amplitude (in flux) of the extragalactic point sources that will be observed by the Planck Mission is presented in this paper. The study is based on the Mexican Hat wavelet formalism introduced by Cayon et al. 2000. Simulations at Planck observing frequencies are analysed, taking into account all the possible cosmological, Galactic and Extragalactic emissions together with noise. With the technique used in this work the Planck Mission will produce a catalogue of extragalactic point sources above fluxes: 1.03 Jy (857 GHz), 0.53 Jy (545 GHz), 0.28 Jy (353 GHz), 0.24 Jy (217 GHz), 0.32 Jy (143 GHz), 0.41 Jy (100 GHz HFI), 0.34 Jy (100 GHz LFI), 0.57 Jy (70 GHz), 0.54 Jy (44 GHz) and 0.54 Jy (30 GHz), which are only slightly model dependent (see text). Amplitudes of these sources are estimated with errors below 15%. Moreover, we also provide a complete catalogue (for the point sources simulation analysed) with errors in the estimation of the amplitude below 10%. In addition we discuss the possibility of identifying different point source populations in the Planck catalogue by estimating their spectral indices.

We show how future measurements of the SZ effect (SZE) can be used to constrain the cosmological parameters. We combine the SZ information expected from the Planck full-sky survey, N(S), where no redshift information is included, with the N(z) obtained from an optically-identified SZ-selected survey covering less than 1% of the sky. We demonstrate how with a small subsample (approx. 300 clusters) of the whole SZ catalogue observed optically it is possible to drastically reduce the degeneracy among the cosmological parameters. We have studied the requirements for performing the optical follow-up and we show the feasibility of such a project. Finally we have compared the cluster expectations for Planck with those expected for Newton-XMM during their lifetimes. It is shown that, due to its larger sky coverage, Planck will detect a factor approx. 5 times more clusters than Newton-XMM and also with a larger redshift coverage.

We present a combined maximum-entropy method (MEM) and Mexican Hat wavelet (MHW) analysis in order to recover the different components of the microwave sky. We apply this technique to simulated observations by the ESA Planck satellite in small patches of the sky. In particular, the introduction of the MHW allows one to detect and subtract the brightest point sources present in the input data and therefore to improve the reconstructions of the CMB and foreground components achieved by MEM on its own. In addition, a point source catalogue at each Planck frequency is produced, which is more complete and accurate than those obtained by each technique independently.

We study two different methods to test Gaussianity in CMB maps. One of them is based on the partition function and the other on the morphology of hot and cold spots. The partition function contains information on all the moments and scales, being a useful quantity to compress the large data sets expected from future space missions like Planck. In particular, it contains much richer information than the one available through the radiation power spectrum. The second method utilizes morphological properties of hot and cold spots such as the eccentricity and number of spots in CMB maps. We study the performance of both methods in detecting non-Gaussian features in small scale CMB simulated maps as those which will be provided by the Planck mission.

In this work we show how galaxy clusters can be used to discriminate among different cosmological models. We have used available X-ray & optical cluster data to constrain the cosmological parameters as well as the cluster scaling relations, T-M and L-T. We also show the power of future SZE data to constrain even more these parameters.

We investigate the power of geometrical estimators on detecting non-Gaussianity in the cosmic microwave background. In particular the number, eccentricity and Gaussian curvature of excursion sets above (and below) a threshold are studied. We compare their different performance when applied to non-Gaussian simulated maps of small patches of the sky, which take into account the angular resolution and instrumental noise of the Planck satellite. These non-Gaussian simulations are obtained as perturbations of a Gaussian field in two different ways which introduce a small level of skewness or kurtosis in the distribution. A comparison with a classical estimator, the genus, is also shown. We find that the Gaussian curvature is the best of our estimators in all the considered cases. Therefore we propose the use of this quantity as a particularly useful test to look for non-Gaussianity in the CMB.

Gravitational lensing predicts an enhancement of the density of bright, distant QSOs around foreground galaxies. We measure this QSO-galaxy correlation w_qg for two complete samples of radio-loud quasars, the southern 1Jy and Half-Jansky samples. The existence of a positive correlation between z~1 quasars and z~0.15 galaxies is confirmed at a p=99.0% significance level (>99.9%) if previous measurements on the northern hemisphere are included). A comparison with the results obtained for incomplete quasar catalogs (e.g. the Veron-Cetty and Veron compilation) suggests the existence of an `identification bias', which spuriously increases the estimated amplitude of the quasar-galaxy correlation for incomplete samples. This effect may explain many of the surprisingly strong quasar-galaxy associations found in the literature. Nevertheless, the value of w_qg that we measure in our complete catalogs is still considerably higher than the predictions from weak lensing. We consider two effects which could help to explain this discrepancy: galactic dust extinction and strong lensing.

In this work we show that combining different cluster data sets is a powerful tool to constrain both, the cosmology and cluster properties. We assume a model with 9 parameters and fit them to 5 cluster data sets. From that fit, we conclude that only low density universes are compatible with the previous data and we also get some interesting conclusions on the rest of parameters.

We consider the 2D wavelet transform with two scales to study sky maps of temperature anisotropies in the cosmic microwave background radiation (CMB). We apply this technique to simulated maps of small sky patches of size 12.8 \times 12.8 square degrees and 1.5' \times 1.5' pixels. The relation to the standard approach, based on the cl's is established through the introduction of the scalogram. We consider temperature fluctuations derived from standard, open and flat-Lambda CDM models. We analyze CMB anisotropies maps plus uncorrelated Gaussian noise (uniform and non-uniform) at idfferent S/N levels. We explore in detail the denoising of such maps and compare the results with other techniques already proposed in the literature. Wavelet methods provide a good reconstruction of the image and power spectrum. Moreover, they are faster than previously proposed methods.

Analysis and denoising of Cosmic Microwave Background (CMB) maps are performed using wavelet multiresolution techniques. The method is tested on $12^{\circ}.8\times 12^{\circ}.8$ maps with resolution resembling the experimental one expected for future high resolution space observations. Semianalytic formulae of the variance of wavelet coefficients are given for the Haar and Mexican Hat wavelet bases. Results are presented for the standard Cold Dark Matter (CDM) model. Denoising of simulated maps is carried out by removal of wavelet coefficients dominated by instrumental noise. CMB maps with a signal-to-noise, $S/N \sim 1$, are denoised with an error improvement factor between 3 and 5. Moreover we have also tested how well the CMB temperature power spectrum is recovered after denoising. We are able to reconstruct the $C_{\ell}$'s up to $l\sim 1500$ with errors always below $20% $ in cases with $S/N \ge 1$.

We propose an alternative analysis of the microwave background temperature anisotropy maps that is based on the study of the roughness of natural surfaces. We apply it to large angle anisotropies, such as those measured by COBE-DMR. We show that for a large signal to noise experiment, the spectral index can be determined independently of the normalization. We then analyze the 4 yr COBE map and find for a flat $\Omega=1$ universe, that the best-fitting value for the spectral index is $n = 1.15^{+0.39}_{-0.34}$ and for the amplitude $Q_{rms-PS}= 14.1^{+3.9}_{-3.5}\mu K$. For $n=1$, the best-fitting normalization is $Q_{rms-PS}|_{n=1}= 16.2^{+1.4}_{-1.3}\mu K$.

We present an alternative method to analyse cosmic microwave background (CMB) maps. We base our analysis on the study of the partition function. This function is used to examine the CMB maps making use of the different information embedded at different scales and moments. Using the partition function in a likelihood analysis in two dimensions (Q_rms,n), we find the best-fitting model to the best data available at present the COBE--DMR 4 years data set. By means of this analysis we find a maximum in the likelihood function for n=1.8 (-0.65 +0.35) and Q_rms-PS = 10 (-2.5 +3) muK (95 % confidence level) in agreement with the results of other similar analyses (Smoot et al. 1994 (1 yr), Bennet et al. 1996 (4 yr)). Also making use of the partition function we perform a multifractal analysis and study the possible fractal nature of the CMB sky. We find that the measure used in the analysis is not a fractal. Finally, we use the partition function for testing the statistical distribution of the COBE--DMR data set. We conclude that no evidence of non-Gaussianity can be found by means of this method.

We imaged the X-ray error box of GRB 980703, beginning 22.5 hours after the gamma--ray event, in both the optical R and near-infrared H bands. A fading optical/IR object was detected within the X-ray error box, coincident with the variable radio source reported by Frail et al. (1998a), who also detected the optical transient independently of us. Further imagery revealed the GRB host galaxy, with R = 22.49 +/- 0.04 and H = 20.5 +/- 0.25, the brightest so far detected. When excluding its contribution to the total flux, both the R and H-band light curves are well-fit by a power-law decay with index alpha~1.4. Our data suggest an intrinsic column density in the host galaxy of ~ 3.5 x 10^21 cm^-2 which indicates the existence of a dense and rich-gas medium in which the GRB occurred, thus supporting the hypernova model scenarios.

We present a K survey of 29 fields covering approximately 90 arcmin^2 from the Medium Deep Survey (MDS) catalogue down to a completeness magnitude of K=18.0 (limiting magnitude K=19.0). The morphology obtained by the MDS team using high resolution images from HST/WFPC2 along with our NIR observations allow a Colour-Magnitude and Colour-Colour analysis that agrees in general with spectral evolution models (Bruzual & Charlot 1998) especially if a reasonable range of metallicities for the Simple Stellar Populations used (0.2 < Z/Zsolar < 2.5) is considered. However, a significant population of spheroids was found, which appear bluer than expected, confirming previous observations (Forbes et al. 1996, Koo et al. 1996, Glazebrook et al. 1998). This blueness might possibly signal the existence of non-negligible star formation in ellipticals and bulges at medium redshift. A number counts calculation for different morphological types show disks become the dominant population at faint magnitudes. The median redshift of the sample is z~0.2, from a photometric redshift estimation using V-K and I-K. A search for EROs in the survey field was also performed, with no detection of objects having I-K > 4.5, setting an upper limit to the number density of EROs at dn_{EROs}/d\Omega < 0.011 arcmin^{-2} (K < 18.0).

We present a list of candidates for high redshift late-type galaxies in the field around the z=3.79 quasar PC1643+4631A. Deep U,V and R imaging has been used to search for objects with a strong Lyman break between U and V, characteristic of galaxies with high hydrogen column densities at $z\sim 3$. A further study of the red objects detected by Hu & Ridgway (1994) has been done, allowing the temptative identification of many of them as low redshift ($z\sim 0.4$) elliptical galaxies.

We present a method, based on the correlation function of excursion sets above a given threshold, to test the Gaussianity of the CMB temperature fluctuations in the sky. In particular, this method can be applied to discriminate between standard inflationary scenarios and those producing non-Gaussianity such as topological defects. We have obtained the normalized correlation of excursion sets, including different levels of noise, for 2-point probability density functions constructed from the Gaussian, \chi_n^2 and Laplace 1-point probability density functions in two different ways. Considering subdegree angular scales, we find that this method can distinguish between different distributions even if the corresponding marginal probability density functions and/or the radiation power spectra are the same.

A thorough study of QSO-galaxy correlations has been done on a region close to the North Galactic Pole using a complete subsample of the optically selected CFHT/MMT QSO survey and the galaxy catalog of Odewahn and Aldering (1995). Although a positive correlation between bright QSOs and galaxies is expected because of the magnification bias effect, none is detected. On the contrary, there is a significant (>99.6%) anticorrelation between z<1.6 QSOs and red galaxies on rather large angular distances. This anticorrelation is much less pronounced for high redshift z>1.6 QSOs, which seems to exclude dust as a cause of the QSO underdensity. This result suggests that the selection process employed in the CFHT/MMT QSO survey is losing up to 50% of low redshift z<1.6 QSOs in regions of high galaxy density. The incompleteness in the whole z<1.6 QSO sample may reach 10% and have important consequences in the estimation of QSO evolution and the QSO autocorrelation function.

Correlations between the magnification or polarization of background sources, induced by gravitational lensing due to the large-scale structure, and the positions of foreground galaxies are investigated. We found that their amplitude is enhanced with respect to correlations for a single population. We analize the dependence of the correlations with the density parameter \Omega considering a nonlinear evolution of the matter power spectrum. The contribution of the linear evolution is negligible at scales below several arcminutes. Detailed results on the dependence of the correlations on the redshift of the foreground and background populations for different cosmological models are presented. The effect reaches its maximum amplitude for foreground populations with relatively small redshifts due to the fast increase of the nonlinear matter power spectrum at recent times.

K-band imaging of a sample of 21 radiogalaxies with z \sim 1.5 reveals the existence of a statistical association of foreground galaxies with the positions of the radiosources. The excess is detected within a 1' radius at a high significance level (>99.8%). K-band light is a good tracer of stellar mass, so this result indicates the existence of an association between foreground mass perturbations and background radiosources, as expected from the magnification bias effect, and confirms previous results obtained for other high-z radio-loud AGN.

We have studied the effect of gravitational lensing on the Cosmic Microwave Background (CMB) anisotropy in flat and open universes. We develop a formalism to calculate the changes on the radiation power spectrum induced by lensing in the Newtonian and synchronous-comoving gauges. The previously considered negligible contribution to the CMB radiation power spectrum of the anisotropic term of the lensing correlation is shown to be appreciable. However, considering the nonlinear evolution of the matter power spectrum produces only slight differences on the results based on linear evolution. The general conclusion for flat as well as open universes is that lensing slightly smoothes the radiation power spectrum. For a given range of multipoles the effect of lensing increases with Omega but for the same acoustic peak it decreases with $\Omega$. The maximum contribution of lensing to the radiation power spectrum for $l\leq 2000$ is $\sim$ 5% for $\Omega$ values in the range 0.1-1.

We present properties of the peaks (maxima) of the CMB anisotropies expected in flat and open CDM models. We obtain analytical expressions of several topological descriptors: mean number of maxima and the probability distribution of the gaussian curvature and the eccentricity of the peaks. These quantities are calculated as functions of the radiation power spectrum, assuming a gaussian distribution of temperature anisotropies. We present results for angular resolutions ranging from 5' to 20' (antenna FWHM), scales that are relevant for the MAP and COBRAS/SAMBA space missions and the ground-based interferometer experiments. Our analysis also includes the effects of noise. We find that the number of peaks can discriminate between standard CDM models, and that the gaussian curvature distribution provides a useful test for these various models, whereas the eccentricity distribution can not distinguish between them.

We have performed RIJHK_S imaging of the field around the z=2.67 quasar Q0836+113, which presents several metal line and a damped Ly\alpha absorption systems in its spectrum. The images reveal the existence of a red K_S=18.9 object \approx 11 arcsec from the quasar. On the basis of the empirical relationships between absorption radius and luminosity we conclude that this object may be the CIV absorber at z=1.82. This could be the first detection of a high redshift galaxy causing high-ionisation absorption. After carefully subtracting the QSO, we do not detect, up to a 3\sigma limiting magnitude for extended objects of K_S=20.8, the damped Ly\alpha absorber apparently detected as a Ly\alpha emitter at z=2.47. It is also suggested, that object ``SW'' from Wolfe et al. (1992) could be the galaxy responsible for the claimed MgII absorption at z=0.37.

We have studied the distribution of $B_J<20.5$ galaxies from the ROE/NRL COSMOS/UKST catalogue around two samples of $z>0.3$ QSOs with similar redshift distributions. The first sample is formed by 144 radio-loud QSOs from the Parkes Catalogue, and the other contains 167 optically selected QSOs extracted from the Large Bright Quasar Survey. It is found that there is a $\approx 99.0%$ significance level excess of COSMOS/UKST galaxies around the PKS QSOs, whereas there is a marginal defect of galaxies around the LBQS QSOs. When the distribution of galaxies around both samples is compared, we found that there is an overdensity of galaxies around the PKS sample with respect to the LBQS sample anticorrelated with the distance from the QSOs at a $99.7%$ significance level. Although this result apparently agrees with the predictions of the double magnification bias, it is difficult to explain by gravitational lensing effects alone; dust in the foreground galaxies and selection effects in the detection of LBQS QSOs should be taken into account. It has been established that the lines of sight to PKS flat-spectrum QSOs go through significatively higher foreground galaxy densities than the directions to LBQS quasars, what may be partially related with the reported reddening of PKS QSOs.

The major contribution to the anisotropy of the temperature of the Cosmic Microwave Background (CMB) radiation is believed to come from the interaction of linear density perturbations with the radiation previous to the decoupling time. Assuming a standard thermal history for the gas after recombination, only the gravitational field produced by the linear density perturbations present on a $\Omega\neq 1$ universe can generate anisotropies at low z (these anisotropies would manifest on large angular scales). However, secondary anisotropies are inevitably produced during the nonlinear evolution of matter at late times even in a universe with a standard thermal history. Two effects associated to this nonlinear phase can give rise to new anisotropies: the time-varying gravitational potential of nonlinear structures (Rees-Sciama RS effect) and the inverse Compton scattering of the microwave photons with hot electrons in clusters of galaxies (Sunyaev-Zeldovich SZ effect). These two effects can produce distinct imprints on the CMB temperature anisotropy. We discuss the amplitude of the anisotropies expected and the relevant angular scales in different cosmological scenarios. Future sensitive experiments will be able to probe the CMB anisotropies beyong the first order primary contribution.

A second--order perturbation theory approach is developed to calculate temperature anisotropies in the cosmic microwave background. Results are given for open universes and fluctuations corresponding to CDM models with either Harrison-Zeldovich (HZ) or Lyth-Stewart-Ratra-Peebles (LSRP) primordial energy--density fluctuation power spectrum. Our perturbation theory approach provides a distinctive multipole contribution as compared to the primary one, the amplitude of the effect being very dependent on normalization. For low--$\Omega$ models, the contribution of the secondary multipoles to the radiation power spectrum is negligible both for standard recombination and reionized scenarios, with the 2--year COBE--DMR normalization. For a flat universe this contribution is $\approx 0.1-10\%$ depending on the reionization history of the universe and on the normalization of the power spectrum.

We have made a likelihood statistical analysis of the angular correlations in the {\it COBE}-DMR two-year sky maps by Monte Carlo simulation of the temperature fluctuations. We assume an open universe and consider as primordial power spectrum the Harrison-Zeldovich one, $P(k)=Ak$. We find that the flatness of the universe is not implied by the data. The quadrupole normalization amplitude, $Q_{rms-PS}$, is related to the density parameter, $\Omega$, by $Q_{rms-PS} = 10.67 + 55.81 \Omega - 128.59 \Omega^2 + 81.26 \Omega^3\ \mu$K. We have determined the p.d.f. of $\Omega$ due to cosmic plus sampling (i.e. $20^\circ$ galactic cut) variance which generically shows a bimodal shape. The uncertainty as given by the r.m.s. is $\approx 0.35$, therefore to better constrain $\Omega$ experiments sensitive to higher multipoles ($l>20$) should be considered.

We have found a statistically significant (99.1 \%) excess of red ($O-E>2$) galaxies with photographic magnitudes $E<19.5$, $O< 21$ taken from the APM Sky Catalogue around $z \sim 1$ radiosources from the 1Jy catalogue. The amplitude, scale and dependence on galaxy colours of the observed overdensity are consistent with its being a result of the magnification bias caused by the weak gravitational lensing of large scale structures at redshift $z \approx 0.2-0.4$ and are hardly explained by other causes, as obscuration by dust.

We have carried out spectroscopic observations of the high redshift ($z\approx 3.4$) radio galaxy 0902+34 at intermediate resolution with the William Herschel Telescope. The dynamical spectral ranges covered are 4600-5480 \AA and 5920-7680 \AA with resolutions of 5.4 \AA and 9.5 \AA, respectively. We detect a continuum that is almost flat and resolve three emission lines: Ly$\alpha$, C IV $\lambda$1549 and He II $\lambda$1640, the last one previously undetected. The line ratios are similar to the typical values found for narrow-line high redshift radio galaxies. Line ratios observed in different regions of the galaxy seem to indicate the presence of strong ionization and/or dust density gradients. We have not detected any Ly$\alpha$ absorption at z=3.3968 (red wing of the Ly$\alpha$ emission line) as might be expected from the absorption found at 21 cm by other authors using the VLA and Arecibo antennas. We discuss possible models for the H I absorbing cloud.

A statistical analysis of texture on the {\it COBE}-DMR first year sky maps based on the genus and spot number is presented. A generalized $\chi^2$ statistic is defined in terms of ``observable'' quantities: the genus and spot density that would be measured by different cosmic observers. This strategy together with the use of Monte Carlo simulations of the temperature fluctuations, including all the relevant experimental parameters, represent the main difference with previous analyses. Based on the genus analysis we find a strong anticorrelation between the quadrupole amplitude $Q_{rms-PS}$ and the spectral index $n$ of the density fluctuation power spectrum at recombination of the form $Q_{rms-PS}= 22.2 \pm 1.7 - (4.7 \pm 1.3) \times n\ \mu$K for fixed $n$, consistent with previous works. The result obtained based on the spot density is consistent with this $Q_{rms-PS} (n)$ relation. In addition to the previous results we have determined, using Monte Carlo simulations, the minimum uncertainty due to cosmic variance for the determination of the spectral index with the genus analysis. This uncertainty is $\delta n\approx 0.2$.

We have taken deep $R$-band images of fields around five radiogalaxies: 0956+47, 1217+36, 3C256, 3C324 and 3C294 with $1<z<2$ . 0956+47 is found to show a double nucleus. Our data on 1217+36 suggest the revision of its classification as a radiogalaxy. We found a statistically significant excess of bright ($19.5<R<22$) galaxies on scales of 2 arcmin around the radiogalaxies (which have $R \approx 21.4$) in our sample. The excess has been determined empirically to be at $\gtrsim 99.5\%$ level. It is remarkable that this excess is not present for $22<R<23.75$ galaxies within the same area, suggesting that the excess is not physically associated to the galaxies but due to intervening groups and then related to gravitational lensing.