The Probe of Inflation and Cosmic Origins (PICO) is a proposed probe-scale space mission consisting of an imaging polarimeter operating in frequency bands between 20 and 800 GHz. We describe the science achievable by PICO, which has sensitivity equivalent to more than 3300 Planck missions, the technical implementation, the schedule and cost.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

The Probe of Inflation and Cosmic Origins (PICO) is an imaging polarimeter that will scan the sky for 5 years in 21 frequency bands spread between 21 and 799 GHz. It will produce full-sky surveys of intensity and polarization with a final combined-map noise level of 0.87 $\mu$K arcmin for the required specifications, equivalent to 3300 Planck missions, and with our current best-estimate would have a noise level of 0.61 $\mu$K arcmin (6400 Planck missions). PICO will either determine the energy scale of inflation by detecting the tensor to scalar ratio at a level $r=5\times 10^{-4}~(5\sigma)$, or will rule out with more than $5\sigma$ all inflation models for which the characteristic scale in the potential is the Planck scale. With LSST's data it could rule out all models of slow-roll inflation. PICO will detect the sum of neutrino masses at $>4\sigma$, constrain the effective number of light particle species with $\Delta N_{\rm eff}<0.06~(2\sigma)$, and elucidate processes affecting the evolution of cosmic structures by measuring the optical depth to reionization with errors limited by cosmic variance and by constraining the evolution of the amplitude of linear fluctuations $\sigma_{8}(z)$ with sub-percent accuracy. Cross-correlating PICO's map of the thermal Sunyaev-Zeldovich effect with LSST's gold sample of galaxies will precisely trace the evolution of thermal pressure with $z$. PICO's maps of the Milky Way will be used to determine the make up of galactic dust and the role of magnetic fields in star formation efficiency. With 21 full sky legacy maps in intensity and polarization, which cannot be obtained in any other way, the mission will enrich many areas of astrophysics. PICO is the only single-platform instrument with the combination of sensitivity, angular resolution, frequency bands, and control of systematic effects that can deliver this compelling, timely, and broad science.

The Probe of Inflation and Cosmic Origins (PICO) is a probe-class mission concept currently under study by NASA. PICO will probe the physics of the Big Bang and the energy scale of inflation, constrain the sum of neutrino masses, measure the growth of structures in the universe, and constrain its reionization history by making full sky maps of the cosmic microwave background with sensitivity 80 times higher than the Planck space mission. With bands at 21-799 GHz and arcmin resolution at the highest frequencies, PICO will make polarization maps of Galactic synchrotron and dust emission to observe the role of magnetic fields in Milky Way's evolution and star formation. We discuss PICO's optical system, focal plane, and give current best case noise estimates. The optical design is a two-reflector optimized open-Dragone design with a cold aperture stop. It gives a diffraction limited field of view (DLFOV) with throughput of 910 square cm sr at 21 GHz. The large 82 square degree DLFOV hosts 12,996 transition edge sensor bolometers distributed in 21 frequency bands and maintained at 0.1 K. We use focal plane technologies that are currently implemented on operating CMB instruments including three-color multi-chroic pixels and multiplexed readouts. To our knowledge, this is the first use of an open-Dragone design for mm-wave astrophysical observations, and the only monolithic CMB instrument to have such a broad frequency coverage. With current best case estimate polarization depth of 0.65 microK(CMB}-arcmin over the entire sky, PICO is the most sensitive CMB instrument designed to date.

The Probe of Inflation and Cosmic Origins (PICO) is a NASA-funded study of a Probe-class mission concept. The top-level science objectives are to probe the physics of the Big Bang by measuring or constraining the energy scale of inflation, probe fundamental physics by measuring the number of light particles in the Universe and the sum of neutrino masses, to measure the reionization history of the Universe, and to understand the mechanisms driving the cosmic star formation history, and the physics of the galactic magnetic field. PICO would have multiple frequency bands between 21 and 799 GHz, and would survey the entire sky, producing maps of the polarization of the cosmic microwave background radiation, of galactic dust, of synchrotron radiation, and of various populations of point sources. Several instrument configurations, optical systems, cooling architectures, and detector and readout technologies have been and continue to be considered in the development of the mission concept. We will present a snapshot of the baseline mission concept currently under development.

We study the possibility for constraining the topology of the Universe by means of the matched circles statistic applied to polarised cosmic microwave background (CMB) anisotropy maps. The advantages of using the CMB polarisation maps in studies of the topology over simply analysing the temperature data as has been done to-date are clearly demonstrated. We test our algorithm to search for pairs of matched circles on simulated CMB maps for a universe with the topology of 3-torus. It is found that the noise levels of both Planck and next generation CMB experiments data are no longer prohibitive and should be low enough to enable the use of the polarisation maps for such studies. For such experiments the minimum radius of the back-to-back matched circles which can be detected are determined. We also showed that the polarisation generated after reionisation does not have an impact on detectability of the matched circles.

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 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.

As well as primary fluctuations, CMB temperature maps contain a wealth of additional information in the form of secondary anisotropies. Secondary effects that can be identified with individual objects, such as the thermal and kinetic Sunyaev-Zel'dovich (SZ) effects due to galaxy clusters, are difficult to unambiguously disentangle from foreground contamination and the primary CMB however. We develop a Bayesian formalism for rigorously characterising anisotropies that are localised on the sky, taking the TSZ and KSZ effects as an example. Using a Gibbs sampling scheme, we are able to efficiently sample from the joint posterior distribution for a multi-component model of the sky with many thousands of correlated physical parameters. The posterior can then be exactly marginalised to estimate properties of the secondary anisotropies, fully taking into account degeneracies with the other signals in the CMB map. We show that this method is computationally tractable using a simple implementation based on the existing Commander component separation code, and also discuss how other types of secondary anisotropy can be accommodated within our framework.

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.

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}$.

We describe a simple but efficient method for deriving a consistent set of monopole and dipole corrections for multi-frequency sky map data sets, allowing robust parametric component separation with the same data set. The computational core of this method is linear regression between pairs of frequency maps, often called "T-T plots". Individual contributions from monopole and dipole terms are determined by performing the regression locally in patches on the sky, while the degeneracy between different frequencies is lifted when ever the dominant foreground component exhibits a significant spatial spectral index variation. Based on this method, we present two different, but each internally consistent, sets of monopole and dipole coefficients for the 9-year WMAP, Planck 2013, SFD 100 um, Haslam 408 MHz and Reich & Reich 1420 MHz maps. The two sets have been derived with different analysis assumptions and data selection, and provides an estimate of residual systematic uncertainties. In general, our values are in good agreement with previously published results. Among the most notable results are a relative dipole between the WMAP and Planck experiments of 10-15 uK (depending on frequency), an estimate of the 408 MHz map monopole of 8.9 +- 1.3 K, and a non-zero dipole in the 1420 MHz map of $0.15 +- 0.03 K pointing towards Galactic coordinates (l,b) = (308,-36) +- 14 degrees. These values represent the sum of any instrumental and data processing offsets, as well as any Galactic or extra-Galactic component that is spectrally uniform over the full sky.

We study hemispherical power asymmetry (Eriksen et. al. 2004, Hansen et. al. 2004, Hansen et. al. 2009) in the WMAP 9-year data. We analyse the combined V- and W-band sky maps, after application of the KQ85 mask, and find that the asymmetry is statistically significant at the 3.4sigma confidence level for ell=2-600, where the data is signal dominated, with a preferred asymmetry direction (l,b)=(227,-27). Individual asymmetry axes estimated from six independent multipole ranges are all consistent with this direction. Subsequently, we estimate cosmological parameters on different parts of the sky and show that the parameters A_s, n_s and Omega_b are the most sensitive to this power asymmetry. In particular, for the two opposite hemispheres aligned with the preferred asymmetry axis, we find n_s = 0.959 \pm 0.022 and n_s = 0.989 \pm 0.024, respectively.

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)

We impose constraints on the topology of the Universe determined from a search for matched circles in the cosmic microwave background (CMB) temperature anisotropy patterns of the 7-year WMAP data. We pay special attention to the sensitivity of the method to residual foreground contamination of the sky maps. A search for pairs of matched back-to-back circles in the higher resolution WMAP W-band map allows tighter constraints to be imposed on topology. Our results rule out universes with topologies that predict pairs of such circles with radii larger than approximately 10 degrees. This places a lower bound on the size of the fundamental domain for a flat universe of about 27.9 Gpc. We study also the possibility for constraining the topology of the Universe by means of the matched circles statistic applied to polarised CMB anisotropy maps. The advantages of using the CMB polarisation maps in studies of the topology over simply analysing the temperature data as has been done to-date are clearly demonstrated. It is found that the noise levels of both Planck and next generation CMB experiments data are no longer prohibitive and should be low enough to enable the use of the polarisation maps for such studies.

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.

We analyse WMAP 7-year temperature data, jointly modeling the cosmic microwave background (CMB) and Galactic foreground emission. We use the Commander code based on Gibbs sampling. Thus, from the WMAP7 data, we derive simultaneously the CMB and Galactic components on scales larger than 1deg with sensitivity improved relative to previous work. We conduct a detailed study of the low-frequency foreground with particular focus on the "microwave haze" emission around the Galactic center. We demonstrate improved performance in quantifying the diffuse galactic emission when Haslam 408MHz data are included together with WMAP7, and the spinning and thermal dust emission is modeled jointly. We also address the question of whether the hypothetical galactic haze can be explained by a spatial variation of the synchrotron spectral index. The excess of emission around the Galactic center appears stable with respect to variations of the foreground model that we study. Our results demonstrate that the new galactic foreground component - the microwave haze - is indeed present.

We demonstrate the feasibility to generate surrogates by Fourier-based methods for an incomplete data set. This is performed for the case of a CMB analysis, where astrophysical foreground emission, mainly present in the Galactic plane, is a major challenge. The shuffling of the Fourier phases for generating surrogates is now enabled by transforming the spherical harmonics into a new set of basis functions that are orthonormal on the cut sky. The results show that non-Gaussianities and hemispherical asymmetries in the CMB as identified in several former investigations, can still be detected even when the complete Galactic plane (|b| < 30{\deg}) is removed. We conclude that the Galactic plane cannot be the dominant source for these anomalies. The results point towards a violation of statistical isotropy.

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.

WMAP data when combined with ancillary data on free-free, synchrotron and dust allow an improved understanding of the spectrum of emission from each of these components. Here we examine the sky variation at intermediate and high latitudes using a cross-correlation technique. In particular, we compare the observed emission in several large partitions of the sky plus 33 selected sky regions to three "standard" templates. The regions are selected using a criterion based on the morphology of these template maps. The synchrotron emission shows evidence of steepening between GHz frequencies and the \emph{WMAP} bands. There are indications of spectral index variations across the sky but the current data are not precise enough to accurately quantify this from region-to-region. The emission correlated with the $H_{\alpha}$ template shows clear evidence of deviation from a free-free spectrum. The emission can be decomposed into a contribution from both free-free and spinning dust in the warm ionised medium of the Galaxy. The derived free-free emissivity corresponds to a mean electron temperature of $\sim 6000$ K, although the value depends critically on the impact of dust absorption on the $H_{\alpha}$ intensity. The WIM spinning dust emission has a peak emission in intensity in the range 40--50 GHz. The anomalous microwave emission associated with dust is detected at high significance in most of the 33 fields studied. The anomalous emission correlates well with the Finkbeiner et al. (1999) model 8 predictions (FDS8) at 94 GHz, and is well described globally by a power-law emission model with an effective spectral index between 20 and 60 GHz of $\beta \approx -2.7$. It is clear that attempts to explain the emission by spinning dust models require multiple components, which presumably relates to a complex mix of emission regions along a given line-of-sight.

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.

In the recent years, non-Gaussianity and statistical isotropy of the Cosmic Microwave Background (CMB) was investigated with various statistical measures, first and foremost by means of the measurements of the WMAP satellite. In this Review, we focus on the analyses that were accomplished with a measure of local type, the so-called Scaling Index Method (SIM). The SIM is able to detect structural characteristics of a given data set, and has proven to be highly valuable in CMB analysis. It was used for comparing the data set with simulations as well as surrogates, which are full sky maps generated by randomisation of previously selected features of the original map. During these investigations, strong evidence for non-Gaussianities as well as asymmetries and local features could be detected. In combination with the surrogates approach, the SIM detected the highest significances for non-Gaussianity to date.

We present a model-independent investigation of the WMAP data with respect to scale- dependent non-Gaussianities (NGs) by employing the method of constrained randomization. For generating so-called surrogate maps a shuffling scheme is applied to the Fourier phases of the original data, which allows to test for the presence of higher order correlations (HOCs) on well-defined scales. Using scaling indices as test statistics we find highly significant signatures for non-Gaussianities when considering all scales. We test for NGs in the bands l = [2,20], l = [20,60], l = [60,120] and l = [120,300]. We find highly significant signatures for non-Gaussianities and ecliptic hemispherical asymmetries for l = [2, 20]. We also obtain highly significant deviations from Gaussianity for the band l = [120,300]. The result for the full l-range can be interpreted as a superposition of the signatures found in the bands l = [2, 20] and l = [120, 300]. We find remarkably similar results when analyzing different ILC-like maps. We perform a set of tests to investigate if the detected anomalies can be explained by systematics. While no test can convincingly rule out the intrinsic nature of the anomalies for the low l case, the ILC map making procedure and/or residual noise in the maps can also lead to NGs at small scales. Our investigations prove that there are phase correlations in the WMAP data of the CMB. In the absence of an explanation in terms of Galactic foregrounds or known systematic artefacts, the signatures at low l must so far be taken to be cosmological at high significance. These findings strongly disagree with predictions of isotropic cosmologies with single field slow roll inflation. The task is now to elucidate the origin of the phase correlations and to understand the physical processes leading to these scale-dependent non-Gaussianities - if systematics as cause for them must be ruled out.

The Planck High Frequency Instrument (HFI) is designed to measure the temperature and polarization anisotropies of the Cosmic Microwave Background and galactic foregrounds in six wide bands centered at 100, 143, 217, 353, 545 and 857 GHz at an angular resolution of 10' (100 GHz), 7' (143 GHz), and 5' (217 GHz and higher). HFI has been operating flawlessly since launch on 14 May 2009. The bolometers cooled to 100 mK as planned. The settings of the readout electronics, such as the bolometer bias current, that optimize HFI's noise performance on orbit are nearly the same as the ones chosen during ground testing. Observations of Mars, Jupiter, and Saturn verified both the optical system and the time response of the detection chains. The optical beams are close to predictions from physical optics modeling. The time response of the detection chains is close to pre-launch measurements. The detectors suffer from an unexpected high flux of cosmic rays related to low solar activity. Due to the redundancy of Planck's observations strategy, the removal of a few percent of data contaminated by glitches does not affect significantly the sensitivity. The cosmic rays heat up significantly the bolometer plate and the modulation on periods of days to months of the heat load creates a common drift of all bolometer signals which do not affect the scientific capabilities. Only the high energy cosmic rays showers induce inhomogeneous heating which is a probable source of low frequency noise.

We develop the XFaster Cosmic Microwave Background (CMB) temperature and polarization anisotropy power spectrum and likelihood technique for the Planck CMB satellite mission. We give an overview of this estimator and its current implementation and present the results of applying this algorithm to simulated Planck data. We show that it can accurately extract the power spectrum of Planck data for the high-l multipoles range. We compare the XFaster approximation for the likelihood to other high-l likelihood approximations such as Gaussian and Offset Lognormal and a low-l pixel-based likelihood. We show that the XFaster likelihood is not only accurate at high-l, but also performs well at moderately low multipoles. We also present results for cosmological parameter Markov Chain Monte Carlo estimation with the XFaster likelihood. As long as the low-l polarization and temperature power are properly accounted for, e.g., by adding an adequate low-l likelihood ingredient, the input parameters are recovered to a high level of accuracy.

Accurate measurement of the cosmic microwave background (CMB) anisotropy requires precise knowledge of the instrument beam. We explore how well the Planck beams will be determined from observations of planets, developing techniques that are also appropriate for other experiments. We simulate planet observations with a Planck-like scanning strategy, telescope beams, noise, and detector properties. Then we employ both parametric and non-parametric techniques, reconstructing beams directly from the time-ordered data. With a faithful parameterization of the beam shape, we can constrain certain detector properties, such as the time constants of the detectors, to high precision. Alternatively, we decompose the beam using an orthogonal basis. For both techniques, we characterize the errors in the beam reconstruction with Monte Carlo realizations. For a simplified scanning strategy, we study the impact on estimation of the CMB power spectrum. Finally, we explore the consequences for measuring cosmological parameters, focusing on the spectral index of primordial scalar perturbations, n_s. The quality of the power spectrum measurement will be significantly influenced by the optical modeling of the telescope. In our most conservative case, using no information about the optics except the measurement of planets, we find that a single transit of Jupiter across the focal plane will measure the beam window functions to better than 0.3% for the channels at 100-217 GHz that are the most sensitive to the CMB. Constraining the beam with optical modeling can lead to much higher quality reconstruction. Depending on the optical modeling, the beam errors may be a significant contribution to the measurement systematics for n_s.

We introduce NeedATool (Needlet Analysis Tool), a software for data analysis based on needlets, a wavelet rendition which is powerful for the analysis of fields defined on a sphere. Needlets have been applied successfully to the treatment of astrophysical and cosmological observations, and in particular to the analysis of cosmic microwave background (CMB) data. Usually, such analyses are performed in real space as well as in its dual domain, the harmonic one. Both spaces have advantages and disadvantages: for example, in pixel space it is easier to deal with partial sky coverage and experimental noise; in harmonic domain, beam treatment and comparison with theoretical predictions are more effective. During the last decade, however, wavelets have emerged as a useful tool for CMB data analysis, since they allow to combine most of the advantages of the two spaces, one of the main reasons being their sharp localisation. In this paper, we outline the analytical properties of needlets and discuss the main features of the numerical code, which should be a valuable addition to the CMB analyst's toolbox.

We assess the strengths and weaknesses of several likelihood formalisms, including the XFaster likelihood. We compare the performance of the XFaster likelihood to that of the Offset Lognormal Bandpower likelihood on simulated data for the Planck satellite. Parameters estimated with these two likelihoods are in good agreement. The advantages of the XFaster likelihood can therefore be realized without compromising performance.

We discuss the effects of inhomogeneous sky-coverage on CMB lens reconstruction, focusing on application to the recently launched Planck satellite. We discuss the mean-field which is induced by noise inhomogeneities, as well as three approaches to lens reconstruction in this context: an optimal maximum-likelihood approach which is computationally expensive to evaluate, and two suboptimal approaches which are less intensive. The first of these is only sub-optimal at the five per-cent level for Planck, and the second prevents biasing due to uncertainties in the noise model.

We present skeleton studies of non-Gaussianity in the CMB temperature anisotropy observed in the WMAP5 data. The local skeleton is traced on the 2D sphere by cubic spline interpolation which leads to more accurate estimation of the intersection positions between the skeleton and the secondary pixels than conventional linear interpolation. We demonstrate that the skeleton-based estimator of non-Gaussianity of the local type (f_NL) - the departure of the length distribution from the corresponding Gaussian expectation - yields an unbiased and sufficiently converged f_NL-likelihood. We analyse the skeleton statistics in the WMAP5 combined V- and W-band data outside the Galactic base-mask determined from the KQ75 sky-coverage. The results are consistent with Gaussian simulations of the the best-fitting cosmological model, but deviate from the previous results determined using the WMAP1 data. We show that it is unlikely that the improved skeleton tracing method, the omission of Q-band data, the modification of the foreground-template fitting method or the absence of 6 extended regions in the new mask contribute to such a deviation. However, the application of the Kp0 base-mask in data processing does improve the consistency with the WMAP1 results. The f_NL-likelihoods of the data are estimated at 9 different smoothing levels. It is unexpected that the best-fit values show positive correlation with the smoothing scales. Further investigation argues against a point-source or goodness-of-fit explanation but finds that about 30% of either Gaussian or f_NL samples having better goodness-of-fit than the WMAP5 show a similar correlation. We present the estimate f_NL=47.3+/-34.9 (1sigma error) determined from the first four smoothing angles and f_NL=76.8+/-43.1 for the combination of all nine. The former result may be overestimated at the 0.21sigma-level because of point sources.

Probing Gaussianity represents one of the key questions in modern cosmology, because it allows to discriminate between different models of inflation. We test for large-scale non-Gaussianities in the cosmic microwave background (CMB) in a model-independent way. To this end, so-called first and second order surrogates are generated by first shuffling the Fourier phases belonging to the scales not of interest and then shuffling the remaining phases for the length scales under study. Using scaling indices as test statistics we find highly significant signatures for both non-Gaussianities and asymmetries on large scales for the WMAP data of the CMB. We find remarkably similar results when analyzing different ILC-maps based on the WMAP five and seven year data. Such features being independent from the map-making procedure would disfavor the fundamental principle of isotropy as well as canonical single-field slow-roll inflation - unless there is some undiscovered systematic error in the collection or reduction of the CMB data or yet unknown foreground contributions.

We study the alignments of the low multipoles of CMB anisotropies with specific directions in the sky (i.e. the dipole, the north Ecliptic pole, the north Galactic pole and the north Super Galactic pole). Performing $10^5$ random extractions we have found that: 1) separately quadrupole and octupole are mildly orthogonal to the dipole but when they are considered together, in analogy to \cite{Copi2006}, we find an unlikely orthogonality at the level of 0.8% C.L.; 2) the multipole vectors associated to $\ell=4$ are unlikely aligned with the dipole at $99.1 %$ C.L.; 3) the multipole vectors associated to $\ell=5$ are mildly orthogonal to the dipole but when we consider only maps that show exactly the same correlation among the multipoles as in the observed WMAP 5yr ILC, these multipole vectors are unlikely orthogonal to the dipole at $99.7 %$ C.L..

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.

We present local extrema studies of two models that introduce a preferred direction into the observed cosmic microwave background (CMB) temperature field. In particular, we make a frequentist comparison of the one- and two-point statistics for the dipole modulation and ACW models with data from the five-year Wilkinson Microwave Anisotropy Probe (WMAP). This analysis is motivated by previously revealed anomalies in the WMAP data, and particularly the difference in the statistical nature of the temperature anisotropies when analysed in hemispherical partitions. The analysis of the one-point statistics indicates that the previously determined hemispherical variance difficulties can be apparently overcome by a dipole modulation field, but new inconsistencies arise if the mean and the l-dependence of the statistics are considered. The two-point correlation functions of the local extrema, the temperature pair product and the point-point spatial pair-count, demonstrate that the impact of such a modulation is to over-`asymmetrise' the temperature field on smaller scales than the wave-length of the dipole or quadrupole, and this is disfavored by the observed data.The results from the ACW model predictions, however, are consistent with the standard isotropic hypothesis. The two-point analysis confirms that the impact of this type of violation of isotropy on the temperature extrema statistics is relatively weak. From this work, we conclude that a model with more spatial structure than the dipole modulated or rotational-invariance breaking models are required to fully explain the observed large-scale anomalies in the WMAP data.

We develop a new Bayesian method for estimating white noise levels in CMB sky maps, and apply this algorithm to the 5-year WMAP data. We assume that the amplitude of the noise RMS is scaled by a constant value, alpha, relative to a pre-specified noise level. We then derive the corresponding conditional density, P(alpha | s, Cl, d), which is subsequently integrated into a general CMB Gibbs sampler. We first verify our code by analyzing simulated data sets, and then apply the framework to the WMAP data. For the foreground-reduced 5-year WMAP sky maps and the nominal noise levels initially provided in the 5-year data release, we find that the posterior means typically range between alpha=1.005 +- 0.001 and alpha=1.010 +- 0.001 depending on differencing assembly, indicating that the noise level of these maps are biased low by 0.5-1.0%. The same problem is not observed for the uncorrected WMAP sky maps. After the preprint version of this letter appeared on astro-ph., the WMAP team has corrected the values presented on their web page, noting that the initially provided values were in fact estimates from the 3-year data release, not from the 5-year estimates. However, internally in their 5-year analysis the correct noise values were used, and no cosmological results are therefore compromised by this error. Thus, our method has already been demonstrated in practice to be both useful and accurate.

Measurements of Cosmic Microwave Background (CMB) anisotropy have served as the best experimental probe of the early universe to date. The inflationary paradigm, inspired in part by the extreme isotropy of the CMB, is now a cornerstone in modern cosmology. Inflation has passed a series of rigorous experimental tests, but we still do not understand the physical mechanism or energy scale behind inflation. A general prediction of inflation and one that can provide certain insights into inflationary physics is a background of primordial gravitational waves. These perturbations leave a distinct signature in the CMB B-modes of polarization. The EPIC (Experimental Probe of Inflationary Cosmology) study team has investigated several CMB polarization mission concepts to carry out a definitive measurement of the inflationary B-mode polarization spectrum. In this report we study a mission with an aperture intermediate between the two missions discussed in our previous report. EPIC-IM's increased aperture allows access to a broader science case than the small EPIC-Low Cost mission. In addition to the search for inflationary gravitational waves, the increase aperture allows us to mine the scale polarization and lensing shear polarization signals down to cosmological limits, so that we extract virtually all the cosmological information available from the CMB. In addition, a modest number of channels operating at higher frequencies allows for an all-sky measurement of polarized Galactic dust, which will provide a rich dataset for Galactic science related to magnetic fields. Using a combination of a large sensitivity focal plane with a new optical design, and an efficient 4K mechanical cooler, EPIC-IM realizes higher sensitivity than EPIC-Comprehensive Science mission.

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

(Abridged)Motivated by the recent results of Hansen et al. (2008) concerning a noticeable hemispherical power asymmetry in the WMAP data on small angular scales, we revisit the dipole modulated signal model introduced by Gordon et al. (2005). This model assumes that the true CMB signal consists of a Gaussian isotropic random field modulated by a dipole, and is characterized by an overall modulation amplitude, A, and a preferred direction, p. Previous analyses of this model has been restricted to very low resolution due to computational cost. In this paper, we double the angular resolution, and compute the full corresponding posterior distribution for the 5-year WMAP data. The results from our analysis are the following: The best-fit modulation amplitude for l <= 64 and the ILC data with the WMAP KQ85 sky cut is A=0.072 +/- 0.022, non-zero at 3.3sigma, and the preferred direction points toward Galactic coordinates (l,b) = (224 degree, -22 degree) +/- 24 degree. The corresponding results for l <~ 40 from earlier analyses was A = 0.11 +/- 0.04 and (l,b) = (225 degree,-27 degree). The statistical significance of a non-zero amplitude thus increases from 2.8sigma to 3.3sigma when increasing l_max from 40 to 64, and all results are consistent to within 1sigma. Similarly, the Bayesian log-evidence difference with respect to the isotropic model increases from Delta ln E = 1.8 to Delta ln E = 2.6, ranking as "strong evidence" on the Jeffreys' scale. The raw best-fit log-likelihood difference increases from Delta ln L = 6.1 to Delta ln L = 7.3. Similar, and often slightly stronger, results are found for other data combinations. Thus, we find that the evidence for a dipole power distribution in the WMAP data increases with l in the 5-year WMAP data set, in agreement with the reports of Hansen et al. (2008).

A new and powerful probe of the origin and evolution of structures in the Universe has emerged and been actively developed over the last decade. In the coming decade, non-Gaussianity, i.e., the study of non-Gaussian contributions to the correlations of cosmological fluctuations, will become an important probe of both the early and the late Universe. Specifically, it will play a leading role in furthering our understanding of two fundamental aspects of cosmology and astrophysics: (i) the physics of the very early universe that created the primordial seeds for large-scale structures, and (ii) the subsequent growth of structures via gravitational instability and gas physics at later times. To date, observations of fluctuations in the Cosmic Microwave Background (CMB) and the Large-Scale Structure of the Universe (LSS) have focused largely on the Gaussian contribution as measured by the two-point correlations (or the power spectrum) of density fluctuations. However, an even greater amount of information is contained in non-Gaussianity and a large discovery space therefore still remains to be explored. Many observational probes can be used to measure non-Gaussianity, including CMB, LSS, gravitational lensing, Lyman-alpha forest, 21-cm fluctuations, and the abundance of rare objects such as clusters of galaxies and high-redshift galaxies. Not only does the study of non-Gaussianity maximize the science return from a plethora of present and future cosmological experiments and observations, but it also carries great potential for important discoveries in the coming decade.

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

We present an extensive frequentist analysis of the one-point statistics (number, mean, variance, skewness and kurtosis) and two-point correlation functions determined for the local extrema of the cosmic microwave background temperature field observed in five-years of Wilkinson Microwave Anisotropy Probe (WMAP) data. Application of a hypothesis test on the one-point statistics indicates a low-variance of hot and cold spots in all frequency bands of the WMAP data. The consistency of the observations with Gaussian simulations of the the best-fitting cosmological model is rejected at the 95% C.L. outside the WMAP KQ75 mask and the northern hemispheres in the Galactic and ecliptic coordinate frames. We demonstrate that it is unlikely that residual Galactic foreground emission contributes to the observed non-Gaussianities. However, the application of a high-pass filter that removes large angular scale power does improve the consistency with the best-fitting cosmological model. Two-point correlation functions of the local extrema are calculated for both the temperature pair product (T-T) and spatial pair-counting (P-P). The T-T observations demonstrate weak correlation on scales below 20 degree and lie completely below the lower 3-sigma confidence region once various temperature thresholds are applied to the extrema determined for the KQ75 mask and northern sky partitions. The P-P correlation structure corresponds to the clustering properties of the temperature extrema, and provides evidence that it is the large angular-scale structures, and some unusual properties thereof, that are intimately connected to the properties of the hot and cold-spots observed in the WMAP five-year data.

A well-tested and validated Gibbs sampling code, that performs component separation and CMB power spectrum estimation, was applied to the {\it WMAP} 5-yr data. Using a simple model consisting of CMB, noise, monopoles and dipoles, a ``per pixel'' low-frequency power-law (fitting for both amplitude and spectral index), and a thermal dust template with fixed spectral index, we found that the low-$\ell$ ($\ell < 50$) CMB power spectrum is in good agreement with the published {\it WMAP}5 results. Residual monopoles and dipoles were found to be small ($\lesssim 3 \mu$K) or negligible in the 5-yr data. We comprehensively tested the assumptions that were made about the foregrounds (e.g. dust spectral index, power-law spectral index prior, templates), and found that the CMB power spectrum was insensitive to these choices. We confirm the asymmetry of power between the north and south ecliptic hemispheres, which appears to be robust against foreground modeling. The map of low frequency spectral indices indicates a steeper spectrum on average ($\beta=-2.97\pm0.21$) relative to those found at low ($\sim$GHz) frequencies.

We present a model-independent method to test for scale-dependent non-Gaussianities in combination with scaling indices as test statistics. Therefore, surrogate data sets are generated, in which the power spectrum of the original data is preserved, while the higher order correlations are partly randomised by applying a scale-dependent shuffling procedure to the Fourier phases. We apply this method to the WMAP data of the cosmic microwave background (CMB) and find signatures for non-Gaussianities on large scales. Further tests are required to elucidate the origin of the detected anomalies.

How did the universe evolve? The fine angular scale (l>1000) temperature and polarization anisotropies in the CMB are a Rosetta stone for understanding the evolution of the universe. Through detailed measurements one may address everything from the physics of the birth of the universe to the history of star formation and the process by which galaxies formed. One may in addition track the evolution of the dark energy and discover the net neutrino mass. We are at the dawn of a new era in which hundreds of square degrees of sky can be mapped with arcminute resolution and sensitivities measured in microKelvin. Acquiring these data requires the use of special purpose telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and the South Pole Telescope (SPT). These new telescopes are outfitted with a new generation of custom mm-wave kilo-pixel arrays. Additional instruments are in the planning stages.

Modern cosmology has sharpened questions posed for millennia about the origin of our cosmic habitat. The age-old questions have been transformed into two pressing issues primed for attack in the coming decade: How did the Universe begin? and What physical laws govern the Universe at the highest energies? The clearest window onto these questions is the pattern of polarization in the Cosmic Microwave Background (CMB), which is uniquely sensitive to primordial gravity waves. A detection of the special pattern produced by gravity waves would be not only an unprecedented discovery, but also a direct probe of physics at the earliest observable instants of our Universe. Experiments which map CMB polarization over the coming decade will lead us on our first steps towards answering these age-old questions.

We repeat and extend the analysis of Eriksen et al 2004 and Hansen et al 2004 testing the isotropy of the Cosmic Microwave Background (CMB) fluctuations. We find that the hemispherical power asymmetry previously reported for the largest scales l=2-40 extend to much smaller scales. In fact, for the full multipole range l=2-600, significantly more power is found in the hemisphere centered at (theta=107 deg., phi=226 deg.) in galactic co-latitude and longitude than in the opposite hemisphere consistent with the previously detected direction of asymmetry for l=2-40. We adopt a model selection test where the direction and amplitude of asymmetry as well as the multipole range are free parameters. A model with an asymmetric distribution of power for l=2-600 is found to be preferred over the isotropic model at the 0.4% significance level taking into account the additional parameters required to describe it. A similar direction of asymmetry is found independently in all six subranges of 100 multipoles between l=2-600 and none of our 9800 isotropic simulated maps show a similarly consistent direction of asymmetry over such a large multipole range. No known systematic effects or foregrounds are found to be able to explain the asymmetry.

We present a new Monte Carlo Markov Chain algorithm for CMB analysis in the low signal-to-noise regime. This method builds on and complements the previously described CMB Gibbs sampler, and effectively solves the low signal-to-noise inefficiency problem of the direct Gibbs sampler. The new algorithm is a simple Metropolis-Hastings sampler with a general proposal rule for the power spectrum, C_l, followed by a particular deterministic rescaling operation of the sky signal. The acceptance probability for this joint move depends on the sky map only through the difference of chi-squared between the original and proposed sky sample, which is close to unity in the low signal-to-noise regime. The algorithm is completed by alternating this move with a standard Gibbs move. Together, these two proposals constitute a computationally efficient algorithm for mapping out the full joint CMB posterior, both in the high and low signal-to-noise regimes.

We describe and implement an exact, flexible, and computationally efficient algorithm for joint component separation and CMB power spectrum estimation, building on a Gibbs sampling framework. Two essential new features are 1) conditional sampling of foreground spectral parameters, and 2) joint sampling of all amplitude-type degrees of freedom (e.g., CMB, foreground pixel amplitudes, and global template amplitudes) given spectral parameters. Given a parametric model of the foreground signals, we estimate efficiently and accurately the exact joint foreground-CMB posterior distribution, and therefore all marginal distributions such as the CMB power spectrum or foreground spectral index posteriors. The main limitation of the current implementation is the requirement of identical beam responses at all frequencies, which restricts the analysis to the lowest resolution of a given experiment. We outline a future generalization to multi-resolution observations. To verify the method, we analyse simple models and compare the results to analytical predictions. We then analyze a realistic simulation with properties similar to the 3-yr WMAP data, downgraded to a common resolution of 3 degree FWHM. The results from the actual 3-yr WMAP temperature analysis are presented in a companion Letter.

Precision measurement of the scalar perturbation spectral index, n_s, from the Wilkinson Microwave Anisotropy Probe temperature angular power spectrum requires the subtraction of unresolved point source power. Here we reconsider this issue. First, we note a peculiarity in the WMAP temperature likelihood's response to the source correction: Cosmological parameters do not respond to increased source errors. An alternative and more direct method for treating this error term acts more sensibly, and also shifts n_s by ~0.3 sigma closer to unity. Second, we re-examine the source fit used to correct the power spectrum. This fit depends strongly on the galactic cut and the weighting of the map, indicating that either the source population or masking procedure is not isotropic. Jackknife tests appear inconsistent, causing us to assign large uncertainties to account for possible systematics. Third, we note that the WMAP team's spectrum was computed with two different weighting schemes: uniform weights transition to inverse noise variance weights at l = 500. The fit depends on such weighting schemes, so different corrections apply to each multipole range. For the Kp2 mask used in cosmological analysis, we prefer source corrections A = 0.012 +/- 0.005 muK^2 for uniform weighting and A = 0.015 +/- 0.005 muK^2 for N_obs weighting. Correcting WMAP's spectrum correspondingly, we compute cosmological parameters with our alternative likelihood, finding n_s = 0.970 +/- 0.017 and sigma_8 = 0.778 +/- 0.045 . This n_s is only 1.8 sigma from unity, compared to the ~2.6 sigma WMAP 3-year result. Finally, an anomalous feature in the source spectrum at l<200 remains, most strongly associated with W-band.

We present several diverse applications of the spherical fast convolution method suggested by Wandelt and Gorski (2001), which is useful for studies of telescope optical properties and for construction of shaped filters for analysis of all-sky data. We study sidelobe pickup in three design concepts for the CALISTO infrared telescope. The beam convolution allows for direct comparison of the performance of each telescope design. At 100 microns, the best of these designs limits sidelobe contamination from Galactic dust emission to < 0.1 MJy/sr for most of the sky with |b| > 25 deg. With the fast convolution method, we illustrate the impact of asymmetric primary beams on the recovery of the temperature power spectrum for the Planck microwave background mission. Finally, we use the fast convolution method to specify a class of orientable filters on the sphere, working through a pedagogical example.

Using a Gibbs sampling algorithm for joint CMB estimation and component separation, we compute the large-scale CMB and foreground posteriors of the 3-yr WMAP temperature data. Our parametric data model includes the cosmological CMB signal and instrumental noise, a single power law foreground component with free amplitude and spectral index for each pixel, a thermal dust template with a single free overall amplitude, and free monopoles and dipoles at each frequency. This simple model yields a surprisingly good fit to the data over the full frequency range from 23 to 94 GHz. We obtain a new estimate of the CMB sky signal and power spectrum, and a new foreground model, including a measurement of the effective spectral index over the high-latitude sky. A particularly significant result is the detection of a common spurious offset in all frequency bands of ~ -13muK, as well as a dipole in the V-band data. Correcting for these is essential when determining the effective spectral index of the foregrounds. We find that our new foreground model is in good agreement with template-based model presented by the WMAP team, but not with their MEM reconstruction. We believe the latter may be at least partially compromised by the residual offsets and dipoles in the data. Fortunately, the CMB power spectrum is not significantly affected by these issues, as our new spectrum is in excellent agreement with that published by the WMAP team. The corresponding cosmological parameters are also virtually unchanged.

We consider the issue of hemispherical power asymmetry in the three-year WMAP data, adopting a previously introduced modulation framework. Computing both frequentist probabilities and Bayesian evidences, we find that the model consisiting of an isotropic CMB sky modulated by a dipole field, gives a substantially better fit to the observations than the purely isotropic model, even when accounting for the larger prior volume. For the ILC map, the Bayesian log-evidence difference is ~1.8 in favour of the modulated model, and the raw improvement in maximum log-likelihood is 6.1. The best-fit modulation dipole axis points toward (l,b) = (225 deg,-27 deg), and the modulation amplitude is 0.114, in excellent agreement with the results from the first-year analyses. The frequentist probability of obtaining such a high modulation amplitude in an isotropic universe is ~1%. These results are not sensitive to data set or sky cut. Thus, the statistical evidence for a power asymmetry anomaly is both substantial and robust, although not decisive, for the currently available data. Increased sky coverage through better foreground handling and full-sky and high-sensitivity polarization maps may shed further light on this issue.

We apply a previously developed Gibbs sampling framework to the foreground corrected 3-yr WMAP polarization data and compute the power spectrum and residual foreground template amplitude posterior distributions. We first analyze the co-added Q- and V-band data, and compare our results to the likelihood code published by the WMAP team. We find good agreement, and thus verify the numerics and data processing steps of both approaches. However, we also analyze the Q- and V-bands separately, allowing for non-zero EB cross-correlations and including two individual foreground template amplitudes tracing synchrotron and dust emission. In these analyses, we find tentative evidence of systematics: The foreground tracers correlate with each of the Q- and V-band sky maps individually, although not with the co-added QV map; there is a noticeable negative EB cross-correlation at l <~ 16 in the V-band map; and finally, when relaxing the constraints on EB and BB, noticeable differences are observed between the marginalized band powers in the Q- and V-bands. Further studies of these features are imperative, given the importance of the low-l EE spectrum on the optical depth of reionization tau and the spectral index of scalar perturbations n_s.

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

We analyze the three-year WMAP temperature anisotropy data seeking to confirm the power spectrum and likelihoods published by the WMAP team. We apply five independent implementations of four algorithms to the power spectrum estimation and two implementations to the parameter estimation. Our single most important result is that we broadly confirm the WMAP power spectrum and analysis. Still, we do find two small but potentially important discrepancies: On large angular scales there is a small power excess in the WMAP spectrum (5-10% at l<~30) primarily due to likelihood approximation issues between 13 <= l <~30. On small angular scales there is a systematic difference between the V- and W-band spectra (few percent at l>~300). Recently, the latter discrepancy was explained by Huffenberger et al. (2006) in terms of over-subtraction of unresolved point sources. As far as the low-l bias is concerned, most parameters are affected by a few tenths of a sigma. The most important effect is seen in n_s. For the combination of WMAP, Acbar and BOOMERanG, the significance of n_s =/ 1 drops from ~2.7 sigma to ~2.3 sigma when correcting for this bias. We propose a few simple improvements to the low-l WMAP likelihood code, and introduce two important extensions to the Gibbs sampling method that allows for proper sampling of the low signal-to-noise regime. Finally, we make the products from the Gibbs sampling analysis publically available, thereby providing a fast and simple route to the exact likelihood without the need of expensive matrix inversions.

Earlier papers introduced a method of accurately estimating the angular cosmic microwave background (CMB) temperature power spectrum based on Gibbs sampling. Here we extend this framework to polarized data. All advantages of the Gibbs sampler still apply, and exact analysis of mega-pixel polarized data sets is thus feasible. These advantages may be even more important for polarization measurements than for temperature measurements. While approximate methods can alias power from the larger E-mode spectrum into the weaker B-mode spectrum, the Gibbs sampler (or equivalently, exact likelihood evaluations) allows for a statistically optimal separation of these modes in terms of power spectra. To demonstrate the method, we analyze two simulated data sets: 1) a hypothetical future CMBPol mission, with the focus on B-mode estimation; and 2) a Planck-like mission, to highlight the computational feasibility of the method.

We compare the performance of multiple codes written by different groups for making polarized maps from Planck-sized, all-sky cosmic microwave background (CMB) data. Three of the codes are based on a destriping algorithm; the other three are implementations of an optimal maximum-likelihood algorithm. Time-ordered data (TOD) were simulated using the Planck Level-S simulation pipeline. Several cases of temperature-only data were run to test that the codes could handle large datasets, and to explore effects such as the precision of the pointing data. Based on these preliminary results, TOD were generated for a set of four 217 GHz detectors (the minimum number required to produce I, Q, and U maps) under two different scanning strategies, with and without noise. Following correction of various problems revealed by the early simulation, all codes were able to handle the large data volume that Planck will produce. Differences in maps produced are small but noticeable; differences in computing resources are large.

Context. A specific example of Bianchi Type VIIh models, i.e. those including universal rotation (vorticity) and differential expansion (shear), has been shown in Jaffe et al. (2005) to correlate unexpectedly with the WMAP first-year data. Aims. We re-assess the signature of this model in the WMAP 3-year data. Methods. The cross-correlation methods are described in Jaffe et al. (2006a). We use the WMAP 3-year data release, including maps for individual years, and perform additional comparisons to assess the influence of both noise and residual foregrounds and eliminate potential non-cosmological sources for the correlation. Results. We confirm that the signal is detected in both the combined 3-year data and the individual yearly sky maps at a level consistent with our original analysis. The significance of the correlation is not affected by either noise or foreground residuals. Conclusions. The results of our previous study are unchanged.

WMAP data when combined with ancillary data on free-free, synchrotron and dust allow an improved understanding of the spectrum of emission from each of these components. Here we examine the sky variation at intermediate latitudes using a cross-correlation technique. In particular, we compare the observed emission in 15 selected sky regions to three ``standard'' templates. The free-free emission of the diffuse ionised gas is fitted by a well-known spectrum at K and Ka band, but the derived emissivity corresponds to a mean electron temperature of ~4000-5000K. This is inconsistent with estimates from galactic HII regions. The origin of the discrepancy is unclear. The anomalous emission associated with dust is clearly detected in most of the 15 fields studied; it correlates well with the Finkbeiner et al. model 8 predictions (FDS8) at 94 GHz, with an effective spectral index between 20 and 60GHz of -2.85. Furthermore, the emissivity varies by a factor of ~2 from cloud to cloud. A modestly improved fit to the anomalous dust at K-band is provided by modulating the template by an estimate of the dust colour temperature, specifically FDS8*T^n. We find a preferred value n~1.6. The synchrotron emission steepens between GHz frequencies and the WMAP bands. There are indications of spectral index variations across the sky but the current data are not precise enough to accurately quantify this from region to region. Our analysis of the WMAP data indicates strongly that the dust-correlated emission at the low WMAP frequencies has a spectrum which is compatible with spinning dust; we find no evidence for a synchrotron component correlated with dust (abridged).

The quality of CMB observations has improved dramatically in the last few years, and will continue to do so in the coming decade. Over a wide range of angular scales, the uncertainty due to instrumental noise is now small compared to the cosmic variance. One may claim with some justification that we have entered the era of precision CMB cosmology. However, some caution is still warranted: The errors due to residual foreground contamination in the CMB power spectrum and cosmological parameters remain largely unquantified, and the effect of these errors on important cosmological parameters such as the optical depth tau and spectral index n_s is not obvious. A major goal for current CMB analysis efforts must therefore be to develop methods that allows us to propagate such uncertainties from the raw data through to the final products. Here we review a recently proposed method that may be a first step towards that goal.

We explore methods of fitting templates to cosmic microwave background (CMB) data, and in particular demonstrate the application of the total convolver algorithm as a fast method of performing a search over all possible locations and orientations of the template relative to the sky. This analysis includes investigation of issues such as chance alignments and foreground residuals. We apply these methods to compare Bianchi models of type VII_h to WMAP first year data and confirm the basic result of our 2005 paper.

We present a new approach to foreground removal for Cosmic Microwave Background (CMB) maps. Rather than relying on prior knowledge about the foreground components, we first extract the necessary information about them directly from the microwave sky maps by taking differences of temperature maps at different frequencies. These difference maps, which we refer to as internal templates, consist only of linear combinations of galactic foregrounds and noise, with no CMB component. We obtain the foreground cleaned maps by fitting these internal templates to, and subsequently subtracting the appropriately scaled contributions of them from, the CMB dominated channels. The fitting operation is performed in wavelet space, making the analysis feasible at high resolution with only a minor loss of precision. Applying this procedure to the WMAP data, we obtain a power spectrum that matches the spectrum obtained by the WMAP team at the signal dominated scales. Finally, we have revisited previous claims about a north-south power asymmetry on large angular scales, and confirm that these remain unchanged with this completely different approach to foreground separation. This also holds when fitting the foreground contribution independently to the northern and southern hemisphere indicating that the asymmetry is unlikely to have its origin in different foreground properties of the hemispheres. This conclusion is further strengthened by the lack of any observed frequency dependence.

Motivated by the success of the Bianchi VII_h model in addressing many of the anomalies observed in the WMAP data (Jaffe et al.), we present calculations in real and in wavelet space of the Higher Criticism statistic of the Bianchi corrected Wilkinson Microwave Anisotropy Probe (WMAP) first year data. At the wavelet scale of 5 degrees the Higher Criticism of the WMAP map drops from a value above the 99% c.l. to a value below the 68% c.l. when corrected by the Bianchi template. An important property of the Higher Criticism statistic is its ability to locate the pixels that account for the deviation from Gaussianity. The analysis of the uncorrected WMAP data pointed to a cold spot in the southern hemisphere centered at (l,b)~ (209, -57). The Higher Criticism of the Bianchi corrected map indicates that this spot remains prominent, albeit at a level completely consistent with Gaussian statistics. Consequently, it is debatable how much emphasis should be placed on this residual feature, but we consider the effect of modestly increasing the scaling of the template. A factor of only 1.2 renders the spot indistinguishable from the background level, with no noticeable impact on the results published in Jaffe et al. for the low-l anomalies, large-scale power asymmetry or wavelet kurtosis. A trivial interpretation would be that the Bianchi template may require a small enhancement of power on scales corresponding to the wavelet scale of 5 degrees.

We generalize the predictions for the CMB anisotropy patterns arising in Bianchi type VIIh universes to include a dark energy component. We consider these models in light of the result of Jaffe et al. (2005a,b) in which a correlation was found on large angular scales between the WMAP data and the anisotropy structure in a low density Bianchi universe. We find that by including a term Omega_L > 0, the same best-fit anisotropy pattern is reproduced by several combinations of cosmological parameters. This sub-set of models can then be further constrained by current observations that limit the values of various cosmological parameters. In particular, we consider the so-called geometric degeneracy in these parameters imposed by the peak structure of the WMAP data itself. Apparently, despite the additional freedom allowed by the dark energy component, the modified Bianchi models are ruled out at high significance.

Motivated by the large-scale asymmetry observed in the cosmic microwave background sky, we consider a specific class of anisotropic cosmological models -- Bianchi type VII_h -- and compare them to the WMAP first-year data on large angular scales. Remarkably, we find evidence of a correlation which is ruled out as a chance alignment at the 3sigma level. The best fit Bianchi model corresponds to x=0.55, Omega_0=0.5, a rotation axis in the direction (l,b)=(222degr,-62degr), shear (sigma/H)_0=2.4e-10 and a right--handed vorticity (omega/H)_0=6.1e-10. Correcting for this component greatly reduces the significance of the large-scale power asymmetry, resolves several anomalies detected on large angular scales (ie. the low quadrupole amplitude and quadrupole/octopole planarity and alignment), and can account for a non--Gaussian "cold spot" on the sky. Despite the apparent inconsistency with the best-fit parameters required in inflationary models to account for the acoustic peaks, we consider the results sufficiently provocative to merit further consideration.

We apply the recently defined multipole vector framework to the frequency-specific first-year WMAP sky maps, estimating the low-l multipole coefficients from the high-latitude sky by means of a power equalization filter. While most previous analyses of this type have considered only heavily processed (and foreground-contaminated) full-sky maps, the present approach allows for greater control of residual foregrounds, and therefore potentially also for cosmologically important conclusions. The low-l spherical harmonics coefficients and corresponding multipole vectors are tabulated for easy reference. Using this formalism, we re-assess a set of earlier claims of both cosmological and non-cosmological low-l correlations based on multipole vectors. First, we show that the apparent l=3 and 8 correlation claimed by Copi et al. (2004) is present only in the heavily processed map produced by Tegmark et al. (2003), and must therefore be considered an artifact of that map. Second, the well-known quadrupole-octopole correlation is confirmed at the 99% significance level, and shown to be robust with respect to frequency and sky cut. Previous claims are thus supported by our analysis. Finally, the low-l alignment with respect to the ecliptic claimed by Schwarz et al. (2004) is nominally confirmed in this analysis, but also shown to be very dependent on severe a-posteriori choices. Indeed, we show that given the peculiar quadrupole-octopole arrangement, finding such a strong alignment with the ecliptic is not unusual.

We propose a solution to the CMB component separation problem based on standard parameter estimation techniques. We assume a parametric spectral model for each signal component, and fit the corresponding parameters pixel by pixel in a two-stage process. First we fit for the full parameter set (e.g., component amplitudes and spectral indices) in low-resolution and high signal-to-noise ratio maps using MCMC, obtaining both best-fit values for each parameter, and the associated uncertainty. The goodness-of-fit is evaluated by a chi^2 statistic. Then we fix all non-linear parameters at their low-resolution best-fit values, and solve analytically for high-resolution component amplitude maps. This likelihood approach has many advantages: The fitted model may be chosen freely, and the method is therefore completely general; all assumptions are transparent; no restrictions on spatial variations of foreground properties are imposed; the results may be rigorously monitored by goodness-of-fit tests; and, most importantly, we obtain reliable error estimates on all estimated quantities. We apply the method to simulated Planck and six-year WMAP data based on realistic models, and show that separation at the muK level is indeed possible in these cases. We also outline how the foreground uncertainties may be rigorously propagated through to the CMB power spectrum and cosmological parameters using a Gibbs sampling technique.

Following this astro-ph communication, we release a set of 10,000 WMAP ILC simulations, produced as described by Eriksen et al. (2004). We strongly encourage that an analysis of these simulations accompanies any scientific analysis of the observed ILC map, in order to assess the stability of the results. For examples of such analyses, we refer the reader to our original paper. In particular, Figure 1c) and d) show the bias and uncertainty induced by the ILC method, and Figure 7 shows a sobering examples of a low-l analysis; the reconstructed quadrupole-octopole alignment is plotted against the true input alignment, and the correlation is underwhelming. It can hardly be emphasized too strongly -- the heavily processed full-sky WMAP maps are not reliable for quantitative analysis, neither at large nor small scales, and they should only, at best, be used as supportive evidence. We hope that these simulations will reduce some of the confusion regarding these maps that currently exists in the community.

Recent analyses of the WMAP data have suggested that the low-order multipoles of the CMB anisotropy distribution show cosmologically interesting and unexpected morphologies and amplitudes. In this paper, we apply a power equalization (PE) filter to the high-latitude WMAP data in order to reconstruct these low-l multipoles free from the largest Galactic foreground modelling uncertainties in the Galactic plane. The characteristic spatial distributions of the modes of order l=2,3,4,5 have been determined as a function of frequency, sky coverage and two methods for foreground correction: using the template fitting results of Bennett et al. (2003) and by a simple linear projection scheme assuming the spectral dependence of the foreground components. Our analyses suggest that the K and Ka frequency bands remain significantly contaminated by residual foreground emission for the WMAP Kp2 mask. However, the l=3,4,5 multipole maps for the Q, V and W channels indicate that, after foreground cleaning, these multipoles are dominated by the CMB anisotropy component. We confirm the Tegmark et al. (2003) result that the octopole does indeed show structure in which its hot and cold spots are centred on a single plane in the sky, and show further that this is very stable with respect to the applied mask and foreground correction. The estimated quadrupole is much less stable showing non-negligible dependence on the Galactic foreground correction. Including these uncertainties is likely to weaken the statistical significance of the claimed alignment between the quadrupole and octopole. Nevertheless, these anisotropy patterns are also present for in COBE-DMR data, and are unlikely to be associated with instrumental systematic artifacts.

We present the first results from a Bayesian analysis of the WMAP first year data using a Gibbs sampling technique. Using two independent, parallel supercomputer codes we analyze the WMAP Q, V and W bands. The analysis results in a full probabilistic description of the information the WMAP data set contains about the power spectrum and the all-sky map of the cosmic microwave background anisotropies. We present the complete probability distributions for each C_l including any non-Gaussianities of the power spectrum likelihood. While we find good overall agreement with the previously published WMAP spectrum, our analysis uncovers discrepancies in the power spectrum estimates at low l multipoles. For example we claim the best-fit Lambda-CDM model is consistent with the C_2 inferred from our combined Q+V+W analysis with a 10% probability of an even larger theoretical C_2. Based on our exact analysis we can therefore attribute the "low quadrupole issue" to a statistical fluctuation.

We compute the two-, three- and four-point correlation functions from the Wilkinson Microwave Anisotropy Probe (WMAP) first-year data, and compare these to a Monte Carlo ensemble of 5000 realizations, based on the best-fit WMAP running-index spectrum of Gaussian fluctuations. The analysis is carried out in three steps, covering small (< 72'), intermediate (< 5 degrees) and large scales (up to 180 degrees). On the largest scales our results are consistent with the previously reported hemisphere power asymmetries: the northern ecliptic hemisphere is practically devoid of large scale fluctuations, while the southern hemisphere show relatively strong fluctuations. We also detect excess correlations in W-band difference maps as compared to the detailed noise simulations produced by the WMAP team, possibly indicative of unknown systematics. While unlikely to affect any temperature based results, this effect could potentially be important for the upcoming polarization data. On intermediate angular scales we find hints of a similar anisotropic distribution of power as seen on the very largest scales, but not to the same extent. In general the model is accepted on these scales. Finally, the same is also true on the smallest scales probed in this paper.

We study the Blackwell-Rao (BR) estimator of the probability distribution of the angular power spectrum, P(C_l|d), by applying it to samples of full-sky no-noise CMB maps generated via Gibbs sampling. We find the estimator, given a set of samples, to be very fast and also highly accurate, as determined by tests with simulated data. We also find that the number of samples required for convergence of the BR estimate rises rapidly with increasing l, at least at low l. Our existing sample chains are only long enough to achieve convergence at l less than about 40. In comparison with P(C_l|d) as reported by the WMAP team we find significant differences at these low l values. These differences lead to up to approximately 0.5 sigma shifts in the estimates of parameters in a 7-parameter Lambda CDM model with non-zero dn_s/dln(k). Fixing dn_s/dln(k)= 0 makes these shifts much less significant. Unlike existing analytic approximations, the BR estimator can be straightforwardly extended for the case of power spectra from correlated fields, such as temperature and polarization. We discuss challenges to extending the procedure to higher l and provide some solutions.

HEALPix -- the Hierarchical Equal Area iso-Latitude Pixelization -- is a versatile data structure with an associated library of computational algorithms and visualization software that supports fast scientific applications executable directly on very large volumes of astronomical data and large area surveys in the form of discretized spherical maps. Originally developed to address the data processing and analysis needs of the present generation of cosmic microwave background (CMB) experiments (e.g. BOOMERanG, WMAP), HEALPix can be expanded to meet many of the profound challenges that will arise in confrontation with the observational output of future missions and experiments, including e.g. Planck, Herschel, SAFIR, and the Beyond Einstein CMB polarization probe. In this paper we consider the requirements and constraints to be met in order to implement a sufficient framework for the efficient discretization and fast analysis/synthesis of functions defined on the sphere, and summarise how they are satisfied by HEALPix.

We study the Internal Linear Combination (ILC) method presented by the Wilkinson Microwave Anisotropy Probe (WMAP) science team, with the goal of determining whether it may be used for cosmological purposes, as a template-free alternative to existing foreground correction methods. We conclude that the method does have the potential to do just that, but great care must be taken both in implementation, and in a detailed understanding of limitations caused by residual foregrounds which can still affect cosmological results. As a first step we demonstrate how to compute the ILC weights both accurately and efficiently by means of Lagrange multipliers, and apply this method to the observed data to produce a new version of the ILC map. Next we describe how to generate Monte Carlo simulations of the ILC map, and find that these agree well with the observed map on angular scales up to l~200. Finally we make two comments to the on-going debates concerning the large-scale properties of the WMAP data. First, we note that the Galactic south-eastern quadrant is associated with notably different ILC weights than the other three quadrants, possibly indicating a foreground related anisotropy. Second, we study the properties of the quadrupole and octopole (amplitude, alignment and planarity), and confirm the previously reported anomalies. Even more interestingly, we find that the l=5 mode is spherically symmetric at about 3 sigma, and that the l=6 mode is planar at the 2 sigma level. However, we also assess the impact of residual foregrounds on these statistics, and find that the ILC map is not clean enough to allow for cosmological conclusions. Alternative methods must be developed to study these issues further.

The three Minkowski functionals and the recently defined length of the skeleton are estimated for the co-added first-year Wilkinson Microwave Anisotropy Probe (WMAP) data and compared with 5000 Monte Carlo simulations, based on Gaussian fluctuations with the a-priori best-fit running-index power spectrum and WMAP-like beam and noise properties. Several power spectrum-dependent quantities, such as the number of stationary points, the total length of the skeleton, and a spectral parameter, gamma, are also estimated. While the area and length Minkowski functionals and the length of the skeleton show no evidence for departures from the Gaussian hypothesis, the northern hemisphere genus has a chi^2 that is large at the 95% level for all scales. For the particular smoothing scale of 3.40 degrees FWHM it is larger than that found in 99.5% of the simulations. In addition, the WMAP genus for negative thresholds in the northern hemisphere has an amplitude that is larger than in the simulations with a significance of more than 3 sigma. On the smallest angular scales considered, the number of extrema in the WMAP data is high at the 3 sigma level. However, this can probably be attributed to the effect of point sources. Finally, the spectral parameter gamma is high at the 99% level in the northern Galactic hemisphere, while perfectly acceptable in the southern hemisphere. The results provide strong evidence for the presence of both non-Gaussian behavior and an unexpected power asymmetry between the northern and southern hemispheres in the WMAP data.

We revisit a recently introduced power spectrum estimation technique based on Gibbs sampling, with the goal of applying it to the high-resolution WMAP data. In order to facilitate this analysis, a number of sophistications have to be introduced, each of which is discussed in detail. We have implemented two independent versions of the algorithm to cross-check the computer codes, and to verify that a particular solution to any given problem does not affect the scientific results. We then apply these programs to simulated data with known properties at intermediate (N_side = 128) and high (N_side = 512) resolutions, to study effects such as incomplete sky coverage and white vs. correlated noise. From these simulations we also establish the Markov chain correlation length as a function of signal-to-noise ratio, and give a few comments on the properties of the correlation matrices involved. Parallelization issues are also discussed, with emphasis on real-world limitations imposed by current super-computer facilities. The scientific results from the analysis of the first-year WMAP data are presented in a companion letter.

Several recent papers have made claims about the detection of an asymmetric distribution of large scale power in the cosmic microwave background anisotropy as measured by the WMAP satellite. In this paper, we investigate how the estimates of particular cosmological parameters vary when inferred from power spectra computed separately on the northern and southern hemispheres in three frames of reference. We focus our investigation on the three specific parameters to which the analysis is most sensitive, the spectral index n, the amplitude of fluctuations A and the optical depth tau. When using a gaussian prior on the spectral index n centred at n=1 with a flat prior on the optical depth tau, the preferred value for the optical depth (derived in the reference frame of maximum asymmetry) in the north is tau=0 with upper limit tau<0.08, whereas in the south we find tau=0.24 +0.06/-0.07 (68% confidence level). The latter result is inconsistent with tau=0 at the 2 sigma level. The estimated optical depth of tau=0.17 on the (nearly) full sky found by the WMAP collaboration and confirmed independently here, could thus in large part originate in structure associated with the southern hemisphere. Furthermore, putting a prior on tau, we find values of the spectral index which are inconsistent between the opposite hemispheres. The exact values depend on the prior on tau. While our results should not be considered more than suggestive, the significance of the parameter differences in the two hemispheres being typically of order 2 sigma, if they are confirmed with the higher sensitivity WMAP 2-year data then it may be necessary to question the assumption of cosmological isotropy and the conceptual framework within which studies of the CMB anisotropy are made.

A flexible maximum-entropy component separation algorithm is presented that accommodates anisotropic noise, incomplete sky-coverage and uncertainties in the spectral parameters of foregrounds. The capabilities of the method are determined by first applying it to simulated spherical microwave data sets emulating the COBE-DMR, COBE-DIRBE and Haslam surveys. Using these simulations we find that is very difficult to determine unambiguously the spectral parameters of the galactic components for this data set due to their high level of noise. Nevertheless, we show that is possible to find a robust CMB reconstruction, especially at the high galactic latitude. The method is then applied to these real data sets to obtain reconstructions of the CMB component and galactic foreground emission over the whole sky. The best reconstructions are found for values of the spectral parameters: T_d=19 K, alpha_d=2, beta_ff=-0.19 and beta_syn=-0.8. The CMB map has been recovered with an estimated statistical error of \sim 22 muK on an angular scale of 7 degrees outside the galactic cut whereas the low galactic latitude region presents contamination from the foreground emissions.

We apply the Gabor transform methodology proposed in (Hansen et al. 2002, 2003) to the WMAP data in order to test the statistical properties of the CMB fluctuation field and specifically to evaluate the fundamental assumption of cosmological isotropy. We estimate the angular power spectrum locally in discs of various sizes positioned in different directions. We determine whether the spatial position of a few outliers in the angular power spectrum could suggest the presence of residual foregrounds or systematic effects. For multipoles close to the first peak, the most deviant local estimates from the best fit WMAP model are associated with a few particular areas close to the Galactic plane. Such deviations also include the ``dent'' in the spectrum just shortward of the first peak which was remarked upon by the WMAP team. Estimating the angular power spectrum excluding these areas gives a slightly higher first Doppler peak amplitude. Finally, we probe the isotropy of the largest angular scales by estimating the power spectrum on hemispheres and reconfirm strong indications of a north-south asymmetry previously reported by other authors. We find tentative evidence for residual foregrounds in the range l=2-4, which could be associated with the low measured quadrupole amplitudes and other anomalies on these angular scales. However, over the range l=5-40 the observed asymmetry is much harder to explain in terms of residual foregrounds and known systematic effects. Such asymmetric distributions of power on the sky provide a serious test for the cosmological principle of isotropy.

We report on the results from two independent but complementary statistical analyses of the WMAP first-year data, based on the power spectrum and N-point correlation functions. We focus on large and intermediate scales (larger than about 3 degrees) and compare the observed data against Monte Carlo ensembles with WMAP-like properties. In both analyses, we measure the amplitudes of the large-scale fluctuations on opposing hemispheres and study the ratio of the two amplitudes. The power-spectrum analysis shows that this ratio for WMAP, as measured along the axis of maximum asymmetry, is high at the 95%-99% level (depending on the particular multipole range included). The axis of maximum asymmetry of the WMAP data is weakly dependent on the multipole range under consideration but tends to lie close to the ecliptic axis. In the N-point correlation function analysis we focus on the northern and southern hemispheres defined in ecliptic coordinates, and we find that the ratio of the large-scale fluctuation amplitudes is high at the 98%-99% level. Furthermore, the results are stable with respect to choice of Galactic cut and also with respect to frequency band. A similar asymmetry is found in the COBE-DMR map, and the axis of maximum asymmetry is close to the one found in the WMAP data.

We develop, implement and test a set of algorithms for estimating N-point correlation functions from pixelized sky maps. These algorithms are slow, in the sense that they do not break the O(N_pix^N) barrier, and yet, they are fast enough for efficient analysis of data sets up to several hundred thousand pixels. The typical application of these methods is Monte Carlo analysis using several thousand realizations, and therefore we organize our programs so that the initialization cost is paid only once. The effective cost is then reduced to a few additions per pixel multiplet (pair, triplet etc.). Further, the algorithms waste no CPU time on computing undesired geometric configurations, and, finally, the computations are naturally divided into independent parts, allowing for trivial (i.e., optimal) parallelization.

With the advent of all-sky H-Alpha surveys it is possible to determine a reliable free-free template of the diffuse interstellar medium (Dickinson, Davies & Davis 2003) which can be used in conjunction with the synchrotron and dust templates to correct CMB observations for diffuse Galactic foregrounds. We have used the COBE-DMR data at 31.5, 53 and 90 GHz and employed cross- correlation techniques to re-evaluate the foreground contributions, particularly that due to dust which is known to be correlated with H-Alpha (and free-free) emission.The DMR microwave maps are found to contain, as well as the expected synchrotron and free-free components, a component tightly correlated to the 140 micron dust maps of DIRBE. At 31.5, 53 and 90 GHz this emission is 6.3 +/- 0.6, 2.4 +/- 0.4 and 2.2 +/- 0.4 microK/(MJy/sr) at 140 microns respectively. When corrected for the contribution from thermal dust following model 7 of Finkbeiner, Davis & Schlegel (1999), a strong anomalous dust emission component remains, which is well-fitted by a frequency spectrum of the form $\nu^{-\beta}$ where beta ~ 2.5 in the DMR frequency range; this is the dominant foreground at 31.5 GHz. The result implies the presence of an emission component with a dust-like morphology but a synchrotron-like spectrum. We discuss the possible origins of this component and compare it with the recent WMAP interpretation. The better knowledge of the individual foregrounds provided by the present study enables a larger area of the sky (|b|>15) to be used to re-appraise the CMB quadrupole normalisation, Qrms-PS, and the scalar perturbations spectral index, n. We find Qrms-PS=15.2^{+2.8}_{-2.3} with a power-law spectral index of n=1.2 +/- 0.2. These values are consistent with previous COBE-DMR analyses and the WMAP 1-year analysis.

We present an application of the fast Independent Component Analysis method to the COBE-DMR 4yr data. Although the signal-to-noise ratio in the COBE-DMR data is typically $\sim 1$, the approach is able to extract the CMB signal with high confidence when working at high galactic latitudes. The reconstructed CMB map shows the expected frequency scaling of the CMB. We fit the resulting CMB component for the rms quadrupole normalisation Qrms and primordial spectral index n and find results in excellent agreement with those derived from the minimum-noise combination of the 90 and 53 GHz DMR channels without galactic emission correction. Including additional channels (priors) such as the Haslam map of radio emission at 408 MHz and the DIRBE 140um map of galactic infra-red emission, the FastICA algorithm is able to both detect galactic foreground emission and separate it from the dominant CMB signal. Fitting the resulting CMB component for Qrms and n we find good agreement with the results from Gorski et al.(1996) in which the galactic emission has been taken into account by subtracting that part of the DMR signal observed to be correlated with these galactic template maps. We further investigate the ability of FastICA to evaluate the extent of foreground contamination in the COBE-DMR data. We include an all-sky Halpha survey (Dickinson, Davies & Davis 2003) to determine a reliable free-free template. In particular we find that, after subtraction of the thermal dust emission predicted by the Finkbeiner, Davis & Schlegel (1999) model 7, this component is the dominant foreground emission at 31.5 GHz. This indicates the presence of an anomalous dust correlated component which is well fitted by a power law spectral shape $\nu^{-\beta}$ with $\beta \sim 2.5$ in agreement with Banday et al. (2003).

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 calculate the two-, three- and (for the first time) four-point correlation functions of the COBE-DMR 4-year sky maps, and search for evidence of non-Gaussianity by comparing the data to Monte Carlo-simulations of the functions. The analysis is performed for the 53 and 90 GHz channels, and five linear combinations thereof. For each map, we simulate an ensemble of 10,000 Gaussian realisations based on an a priori best-fit scale-invariant cosmological power spectrum, the DMR beam pattern and instrument-specific noise properties. Each observed COBE-DMR map is compared to the ensemble using a simple chi^2 statistic, itself calibrated by simulations. In addition, under the assumption of Gaussian fluctuations, we find explicit expressions for the expected values of the four-point functions in terms of combinations of products of the two-point functions, then compare the observed four-point statistics to those predicted by the observed two-point function, using a redefined chi^2 statistic. Both tests accept the hypothesis that the DMR maps are consistent with Gaussian initial perturbations.

We extend the analysis of Gabor transforms on a Cosmic Microwave Background (CMB) temperature map (Hansen, Gorski and Hivon 2002) to polarisation. We study the temperature and polarisation power spectra on the cut sky, the so-called pseudo power spectra. The transformation kernels relating the full-sky polarisation power spectra and the polarisation pseudo power spectra are found to be similar to the kernel for the temperature power spectrum. This fact is used to construct a fast power spectrum estimation algorithm using the pseudo power spectrum of temperature and polarisation as data vectors in a maximum likelihood approach. Using the pseudo power spectra as input to the likelihood analysis solves the problem of having to invert huge matrices which makes the standard likelihood approach infeasible.

The Fourier transform of a dataset apodised with a window function is known as the Gabor transform. In this paper we extend the Gabor transform formalism to the sphere with the intention of applying it to CMB data analysis. The Gabor coefficients on the sphere known as the pseudo power spectrum is studied for windows of different size. By assuming that the pseudo power spectrum coefficients are Gaussian distributed, we formulate a likelihood ansatz using these as input parameters to estimate the full sky power spectrum from a patch on the sky. Since this likelihood can be calculated quickly without having to invert huge matrices, this allows for fast power spectrum estimation. By using the pseudo power spectrum from several patches on the sky together, the full sky power spectrum can be estimated from full-sky or nearly full-sky observations.

In this paper we pursue a new technique to search for evidence of a finite Universe, making use of a spherical mexican-hat wavelet decomposition of the microwave background fluctuations. Using the information provided by the wavelet coefficients at several scales we test whether compact orientable flat topologies are consistent with the COBE-DMR data. We consider topological sizes ranging from half to twice the horizon size. A scale-scale correlation test indicates that non-trivial topologies with appropriate topological sizes are as consistent with the COBE-DMR data as an infinite universe. Among the finite models the data seems to prefer a Universe which is about the size of the horizon for all but the hypertorus and the triple-twist torus. For the latter the wavelet technique does not seem a good discriminator of scales for the range of topological sizes considered here, while a hypertorus has a preferred size which is 80% of the horizon. This analysis allows us to find a best fit topological size for each model, although cosmic variance might limit our ability to distinguish some of the topologies.

We have analysed the angular power spectra of the Parkes radio continuum and polarisation survey of the Southern galactic plane at 2.4 GHz. We have found that in the multipole range l=40-250 the angular power spectrum of the polarised intensity is well described by a power-law spectrum with fitted spectral index alpha_L = 2.37 +- 0.21. In the same multipole range the angular power spectra of the E and B components of the polarised signal are significantly flatter, with fitted spectral indices respectively of alpha_E = 1.57 +- 0.12 and alpha_B = 1.45 +- 0.12. Temperature fluctuations in the E and B components are mostly determined by variations in polarisation angle. We have combined these results with other data from available radio surveys in order to produce a full-sky toy model of Galactic synchrotron intensity and linear polarisation at high frequencies (> 10 GHz). This can be used to study the feasibility of measuring the Cosmic Microwave Background polarisation with forthcoming experiments and satellite missions.

Removal of systematic effects is crucial in present and future CMB experiments mapping large fraction of the sky. Accurate CMB measurements ask for multi-feed array instruments observing the sky with a redundant scanning strategy covering the same region of the sky on different time scales and with different detectors for a better control of systematic effects. We investigate here the capability to suppress 1/f noise features in Time Ordered Data (TOD) by using the destriping technique described in Maino et al. (1999), under realistic assumptions for crossing condition between different scan circles and sky signal fluctuations on small angular scales. We consider as working case, Planck-LFI simulated observations with few arminutes pixel size convolved with LFI beam resolutions. We also extend the analysis to high values of the knee-frequency and found a residual additional noise rms ~31% larger than the pure white noise rms at fk=1 Hz which could be a critical issue in the extraction of the CMB angular power spectrum. Furthermore we verified that destriping quality does not significantly depend on the receiver sensitivity whereas it improves proportionally to the improvement of sampling rate. Therefore given a noise level, the higher the sampling rate, the better the destriping quality.

We propose an estimator for the trispectrum of a scalar random field on a sphere, discuss its geometrical and statistical properties, and outline its implementation. By estimating the trispectrum of the 4 year COBE-DMR data (in HEALPix pixelization) we find new evidence of a non-Gaussian signal associated with a known systematic effect. We find that by removing data from the sky maps for those periods of time perturbed by this effect, the amplitudes of the trispectrum coefficients become completely consistent with predictions for a Gaussian sky. These results reinforce the importance of statistical methods based in harmonic space for quantifying non-Gaussianity.

We present a new, fast, algorithm for the separation of astrophysical components superposed in maps of the sky, based on the fast Independent Component Analysis technique (FastICA). It allows to recover both the spatial pattern and the frequency scalings of the emissions from statistically independent astrophysical processes, present along the line-of-sight, from multi-frequency observations. We apply FastICA to simulated observations of the microwave sky with angular resolution and instrumental noise at the mean nominal levels for the Planck satellite, containing the most important known diffuse signals: the Cosmic Microwave Background (CMB), Galactic synchrotron, dust and free-free emissions. A method for calibrating the reconstructed maps of each component at each frequency has been devised. The spatial pattern of all the components have been recovered on all scales probed by the instrument. In particular, the CMB angular power spectra is recovered at the percent level up to $\ell_{max}\simeq 2000$. Frequency scalings and normalization have been recovered with better than percent precision for all the components at frequencies and in sky regions where their signal-to-noise ratio exceeds 1.5; the error increases at ten percent level for signal-to-noise ratios about 1. Runs have been performed on a Pentium III 600 MHz computer; FastICA typically took a time of the order of 10 minutes for all-sky simulations with 3.5 arcminutes pixel size. We conclude that FastICA is an extremly promising technique for analyzing the maps that will be obtained by the forthcoming high resolution CMB experiments.

Until now, the detailed interpretation of the observed microlensing events has suffered from the fact that the physical parameters of the phenomenon cannot be uniquely determined from the available astronomical measurements, i.e. the photometric lightcurves. The situation will change in the near-future with the availability of long-baseline, sensitive optical interferometers, which should be able to resolve the images of the lensed objects into their components. For this, it will be necessary to achieve a milliarcsecond resolution on sources with typical magnitudes K $\ga 12$. Indeed, brighter events have never been observed up to now by micro-lensing surveys. We discuss the possibilities opened by the use of long baseline interferometry in general, and in particular for one such facility, the ESO VLT Interferometer, which will attain the required performance. We discuss the expected accuracy and limiting magnitude of such measurements. On the basis of the database of the events detected by the OGLE experiment, we estimate the number of microlenses that could be available for measurements by the VLTI. We find that at least several tens of events could be observed each year. In conjunction with the photometric data, our ability to measure the angular separation between the microlensed images will enable a direct and unambiguous determination of both their masses and locations.

We measure the angular bispectrum of the cosmic microwave background (CMB) radiation anisotropy from the COBE Differential Microwave Radiometer (DMR) four-year sky maps. The angular bispectrum is the harmonic transform of the three-point correlation function, analogous to the angular power spectrum, the harmonic transform of the two-point correlation function. First, we study statistical properties of the bispectrum and the normalized bispectrum. We find the latter more useful for statistical analysis; the distribution of the normalized bispectrum is very much Gaussian, while the bare bispectrum distribution is highly non-Gaussian. Then, we measure 466 modes of the normalized bispectrum, all independent combinations of three-point configurations up to a maximum multipole of 20, the mode corresponding to the DMR beam size. By measuring 10 times as many modes as the sum of previous work, we test Gaussianity of the DMR maps. We compare the data with the simulated Gaussian realizations, finding no significant detection of the normalized bispectrum on the mode-by-mode basis. We also find that the previously reported detection of the normalized bispectrum is consistent with a statistical fluctuation. By fitting a theoretical prediction to the data for the primary CMB bispectrum, which is motivated by slow-roll inflation, we put a weak constraint on a parameter characterizing non-linearity in inflation. Simultaneously fitting the foreground bispectra estimated from interstellar dust and synchrotron template maps shows that neither dust nor synchrotron emissions significantly contribute to the bispectrum at high Galactic latitude. We conclude that the DMR map is consistent with Gaussianity.

We describe a fast and accurate method for estimation of the cosmic microwave background (CMB) anisotropy angular power spectrum --- Monte Carlo Apodised Spherical Transform EstimatoR. Originally devised for use in the interpretation of the Boomerang experimental data, MASTER is both a computationally efficient method suitable for use with the currently available CMB data sets (already large in size, despite covering small fractions of the sky, and affected by inhomogeneous and correlated noise), and a very promising application for the analysis of very large future CMB satellite mission products.

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.

We have analysed the Rhodes/HartRAO survey at 2326 MHz and derived the global angular power spectrum of Galactic continuum emission. In order to measure the angular power spectrum of the diffuse component, point sources were removed from the map by median filtering. A least-square fit to the angular power spectrum of the entire survey with a power law spectrum C_l proportional to l^{-alpha}, gives alpha = 2.43 +/- 0.01 for l = 2-100. The angular power spectrum of radio emission appears to steepen at high Galactic latitudes and for observed regions with |b| > 20 deg, the fitted spectral index is alpha = 2.92 +/- 0.07. We have extrapolated this result to 30 GHz (the lowest frequency channel of Planck) and estimate that no significant contribution to the sky temperature fluctuation is likely to come from synchrotron at degree-angular scales

Within the scope of the Planck IDIS (Integrated Data Information System) project we have started to develop the data model for time-ordered data and full-sky maps. The data model is part of the Data Management Component (DMC), a software system designed according to a three-tier architecture which allows complete separation between data storage and processing. The DMC is already being used for simulation activities and the modeling of some foreground components. We have ingested several Galactic surveys into the database and used the science data-access interface to process the data. The data structure for full-sky maps utilises the HEALPix tessellation of the sphere. We have been able to obtain consistent measures of the angular power spectrum of the Galactic radio continuum emission between 408 MHz and 2417 MHz.

We propose fast, exact and efficient algorithms for the convolution of two arbitrary functions on the sphere which speed up computations by a factor \order{\sqrt{N}} compared to present methods where $N$ is the number of pixels. No simplifying assumptions are made other than bandlimitation. This reduces typical computation times for convolving the full sky with the asymmetric beam pattern of a megapixel Cosmic Microwave Background (CMB) mission from months to minutes. Our methods enable realistic simulation and careful analysis of data from such missions, taking into account the effects of asymmetric ``point spread functions'' and far side lobes of the physical beam. While motivated by CMB studies, our methods are general and hence applicable to the convolution or filtering of any scalar field on the sphere with an arbitrary, asymmetric kernel. We show in an appendix that the same ideas can be applied to the inverse problems of map-making and beam reconstruction by similarly accelerating the {\em transpose convolution} which is needed for the iterative solution of the normal equations.

We discuss all-sky convolution of the instrument beam with the sky signal in polarimetry experiments, such as the Planck mission which will map the temperature anisotropy and polarization of the cosmic microwave background (CMB). To account properly for stray light (from e.g. the galaxy, sun, and planets) in the far side-lobes of such an experiment, it is necessary to perform the beam convolution over the full sky. We discuss this process in multipole space for an arbitrary beam response, fully including the effects of beam asymmetry and cross-polarization. The form of the convolution in multipole space is such that the Wandelt-Gorski fast technique for all-sky convolution of scalar signals (e.g. temperature) can be applied with little modification. We further show that for the special case of a pure co-polarized, axisymmetric beam the effect of the convolution can be described by spin-weighted window functions. In the limits of a small angle beam and large Legendre multipoles, the spin-weight 2 window function for the linear polarization reduces to the usual scalar window function used in previous analyses of beam effects in CMB polarimetry experiments. While we focus on the example of polarimetry experiments in the context of CMB studies, we emphasise that the formalism we develop is applicable to anisotropic filtering of arbitrary tensor fields on the sphere.

As the era of high precision cosmology approaches, the empirically determined power spectrum of the microwave background anisotropy $C_l$ will provide a crucial test for cosmological theories. We present an exact semi-analytic framework for the study of the ampling statistics of the $C_l$ resulting from observations with partial sky coverage and anisotropic noise distributions. This includes space-borne, air-borne and ground-based experiments. We apply this theory to demonstrate its power for constructing fast but unbiased approximate methods for the joint estimation of cosmological parameters. Further applications, such as a test for possible non-Gaussianity of the underlying theory and a ``poor man's power spectrum estimator'' are suggested. An appendix derives recursion relations for the efficient computation of the couplings between spherical harmonics on the cut sky.

We present an analysis of the Gaussianity of the 4-year COBE-DMR data (in HEALPix pixelisation) based on spherical wavelets. The skewness, kurtosis and scale-scale correlation spectra are computed from the detail wavelet coefficients at each scale. The sensitivity of the method to the orientation of the data is also taken into account. We find a single detection of non-Gaussianity at the $>99%$ confidence level in one of our statistics. We use Monte-Carlo simulations to assess the statistical significance of this detection and find that the probability of obtaining such a detection by chance for an underlying Gaussian field is as high as 0.69. Therefore, our analysis does not show evidence of non-Gaussianity in the COBE-DMR data.

We investigate the Gaussianity of the 4-year COBE-DMR data (in HEALPix pixelisation) using an analysis based on spherical Haar wavelets. We use all the pixels lying outside the Galactic cut and compute the skewness, kurtosis and scale-scale correlation spectra for the wavelet coefficients at each scale. We also take into account the sensitivity of the method to the orientation of the input signal. We find a detection of non-Gaussianity at $> 99$ per cent level in just one of our statistics. Taking into account the total number of statistics computed, we estimate that the probability of obtaining such a detection by chance for an underlying Gaussian field is 0.69. Therefore, we conclude that the spherical wavelet technique shows no strong evidence of non-Gaussianity in the COBE-DMR data.

We apply the Pseudo-C_l formalism to obtain an unbiased, approximate method for efficient simultaneous estimation of several cosmological parameters from large, almost full-sky cosmic microwave background data sets.

In this paper we pursue the origin of the non-Gaussianity determined by a bispectrum analysis of the COBE-DMR 4-year sky maps. The robustness of the statistic is demonstrated by the rebinning of the data into 12 coordinate systems. By computing the bispectrum statistic as a function of various data partitions - by channel, frequency, and time interval, we show that the observed non-Gaussian signal is driven by the 53 GHz data. This frequency dependence strongly rejects the hypothesis that the signal is cosmological in origin. A jack-knife analysis of the coadded 53 and 90 GHz sky maps reveals those sky pixels to which the bispectrum statistic is particularly sensitive. We find that by removing data from the 53 GHz sky maps for periods of time during which a known systematic effect perturbs the 31 GHz channels, the amplitudes of the bispectrum coefficients become completely consistent with that expected for a Gaussian sky. We conclude that the non-Gaussian signal detected by the normalised bispectrum statistic in the publicly available DMR sky maps is due to a systematic artifact. The impact of removing the affected data on estimates of the normalisation of simple models of cosmological anisotropy is negligible.

We propose a formalism for estimating the skewness and angular power spectrum of a general Cosmic Microwave Background data set. We use the Edgeworth Expansion to define a non-Gaussian likelihood function that takes into account the anisotropic nature of the noise and the incompleteness of the sky coverage. The formalism is then applied to estimate the skewness of the publicly available 4 year Cosmic Background Explorer (COBE) Differential Microwave Radiometer data. We find that the data is consistent with a Gaussian skewness, and with isotropy. Inclusion of non Gaussian degrees of freedom has essentially no effect on estimates of the power spectrum, if each $C_\ell$ is regarded as a separate parameter or if the angular power spectrum is parametrized in terms of an amplitude (Q) and spectral index (n). Fixing the value of the angular power spectrum at its maxiumum likelihood estimate, the best fit skewness is $S=6.5\pm6.0\times10^4(\muK)^3$; marginalizing over Q the estimate of the skewness is $S=6.5\pm8.4\times10^4(\muK)^3$ and marginalizing over n one has $S=6.5\pm8.5\times10^4(\muK)^3$.

We study the impact of the 1/f noise on the PLANCK Low Frequency Instrument (LFI) osbervations (Mandolesi et al 1998) and describe a simple method for removing striping effects from the maps for a number of different scanning stategies. A configuration with an angle between telescope optical axis and spin-axis just less than 90 degrees (namely 85 degress) shows good destriping efficiency for all receivers in the focal plane, with residual noise degradation < 1-2 %. In this configuration, the full sky coverage can be achieved for each channel separately with a 5 degrees spin-axis precession to maintain a constant solar aspect angle.

HEALPix is a Hierarchical, Equal Area, and iso-Latitude Pixelisation of the sphere designed to support efficiently - local operations on the pixel set, - a hierarchical tree structure for multi-resolution applications, and - the global Fast Spherical Harmonic transform. The HEALPix concept and the mathematical software based on it introduced in this primer meet the challenges which future high resolution and large volume CMB data sets including the MAP and Planck mission products will present.

We use Python I, II, and III cosmic microwave background anisotropy data to constrain cosmogonies. We account for the Python beamwidth and calibration uncertainties. We consider open and spatially-flat-Lambda cold dark matter cosmogonies, with nonrelativistic-mass density parameter Omega_0 in the range 0.1--1, baryonic-mass density parameter Omega_B in the range (0.005--0.029) h^{-2}, and age of the universe t_0 in the range (10--20) Gyr. Marginalizing over all parameters but Omega_0, the combined Python data favors an open (spatially-flat-Lambda) model with Omega_0 simeq 0.2 (0.1). At the 2 sigma confidence level model normalizations deduced from the combined Python data are mostly consistent with those drawn from the DMR, UCSB South Pole 1994, ARGO, MAX 4 and 5, White Dish, and SuZIE data sets.

I review our recent claim that there is evidence of non-Gaussianity in the 4 Year COBE DMR data. I describe the statistic we apply, the result we obtain and make a detailed list of the systematics we have analysed. I finish with a qualitative understanding of what it might be and its implications.

We review our recent claim that there is evidence of non-Gaussianity in the 4 Year COBE DMR data. We present some new results concerning the effect of the galactic cut upon the non-Gaussian signal. These findings imply a localization of the non-Gaussian signal on the Northern galactic hemisphere.

We use combinations of ten small-scale cosmic microwave background anisotropy data sets from the UCSB South Pole 1994, ARGO, MAX 4 and 5, White Dish and SuZIE experiments to constrain cosmogonies. We consider open and spatially-flat-Lambda cold dark matter cosmogonies. Marginalizing over all parameters but $\Omega_0$, the combined data favors an Omega_0 simeq 1 for both open and flat-Lambda models. Excluding the smallest angular scale SuZIE data, an Omega_0 simeq 0.3 (1.0) open (flat-Lambda) model is favored. Considering only multi-frequency data with error bars consistent with sample variance and noise considerations, an Omega_0 simeq 0.1 (1) open (flat-Lambda) model is favored. For both open and flat-$\Lambda$ models, after marginalizing over all the other parameters, a lower Omega_B h^2 (sim 0.005) or younger (t_0 sim 10 Gyr) universe is favored. However, the data do not rule out other values of Omega_0 in the flat-Lambda model and other values of Omega_B h^2 in both models. (abriged)

Multi-frequency, high resolution, full sky measurements of the anisotropy in both temperature and polarisation of the cosmic microwave background radiation are the goals of the satellite missions MAP (NASA) and Planck (ESA). The ultimate data products of these missions - multiple microwave sky maps, each of which will have to comprise more than 10^6 pixels in order to render the angular resolution of the instruments - will present serious challenges to those involved in the analysis and scientific exploitation of the results of both surveys. Some considerations of the relevant aspects of the mathematical structure of future CMB data sets are presented in this contribution. >>> for better on-screen rendition of the figures see <<< http://www.tac.dk/~healpix or http://www.mpa-garching.mpg.de/~cosmo/contributions.html

In a recent Letter we have shown how COBE-DMR maps may be used to disprove Gaussianity at a high confidence level. In this report we digress on a few issues closely related to this Letter. We present the general formalism for surveying non-Gaussianity employed. We present a few more tests for systematics. We wonder about the theoretical implications of our result.

Much attention has been given to the problem of estimating cosmological parameters from the $C_l$ measured by future experiments. Many of the approaches which are being used either invoke poorly controlled approximations or are computationally expensive. We derive exact results as well as fast and highly accurate approximations for mapping a theoretical model onto the observed power spectrum coefficients and computing their statistical properties. These results obtain from an analytic framework which applies for any azimuthally symmetric sky coverage regardless of the fraction of the sky observed by the experiment.

As the era of high precision cosmology approaches, the empirically determined power spectrum of the microwave background anisotropy, $C_l$, will provide a crucial test for cosmological theories. We present a unified semi-analytic framework for the study of the statistical properties of the $C_l$ coefficients computed from the results of balloon, ground based, and satellite experiments. An illustrative application shows that commonly used approximations {\it bias} the estimation of the baryon parameter $\Omega_b$ at the 1% level even for a satellite capturing as much as $\sim 70$% of the sky.

We use data from the ARGO cosmic microwave background (CMB) anisotropy experiment to constrain cosmogonies. We account for the ARGO beamwidth and calibration uncertainties, and marginalize over the offset removed from the data. Our derived amplitudes of the CMB anisotropy detected by the ARGO experiment are smaller than those derived previously. We consider open and spatially-flat-Lambda cold dark matter cosmogonies, with clustered-mass density parameter Omega_0 in the range 0.1-1, baryonic-mass density parameter Omega_B in the range (0.005-0.029)h^{-2}, and age of the universe t_0 in the range (10--20) Gyr. Marginalizing over all parameters but Omega_0, the ARGO data favors an open (spatially-flat-Lambda) model with Omega_0= 0.23 (0.1). However, these numerical values are model dependent. At the 2 sigma confidence level model normalizations deduced from the ARGO data are consistent with those drawn from the UCSB South Pole 1994, MAX 4+5, White Dish, and SuZIE data sets. The ARGO open model normalizations are also consistent with those deduced from the DMR data. However, for most spatially-flat-Lambda models the DMR normalizations are more than 2 sigma above the ARGO ones.

We report the development of numerical tools for the topological analysis of sub--degree resolution, all--sky maps. Software to be released in the HEALFAST (V0.9) package defines neighbour relationships for the HEALPIX tessellation of the sphere. We apply this routine to a fast extrema search which scales strictly linearly in the number of pixels, $N_{p}$. We also present a highly efficient algorithm for simulating the gradient vector and curvature tensor fields ``on--the--fly'' with the temperature map, needing only of order $N_{p} \log_2 N_{p}$ more operations.

We present a novel approach to quantifying the morphology of Cosmic Microwave Background (CMB) anisotropy maps. As morphological descriptors, we use shape parameters known as Minkowski functionals. Using the mathematical framework provided by the theory of integral geometry on arbitrary curved supports, we point out the differences to their characterization and interpretation in the case of flat space. With restrictions of real data -- such as pixelization and incomplete sky coverage, to mention just a few -- in mind, we derive and test unbiased estimators for all Minkowski functionals. Various examples, among them the analysis of the four-year COBE DMR data, illustrate the application of our method.

The COBE satellite has provided the only comprehensive multi-frequency full-sky observations of the microwave sky available today. Assessment of the observations requires a detailed likelihood analysis to extract the maximum amount of information present in the noisy data. I present a specific method for estimating the CMB anisotropy power spectrum independent of any assumptions about the underlying cosmology, and then use standard image processing techniques to generate the most revealing corresponding maps of the signal. The consistency of the data at the available frequencies provides strong support to the assertion that we are being provided with our first glimpse of the last scattering surface.

(Shortened Abstract). We analyze a redshift sample of Abell/ACO clusters and compare them with numerical simulations based on the truncated Zel'dovich approximation (TZA), for a list of eleven dark matter (DM) models. For each model we run several realizations, on which we estimate cosmic variance effects. We analyse correlation statistics, the probability density function, and supercluster properties from percolation analysis. As a general result, we find that the distribution of galaxy clusters provides a constraint only on the shape of the power spectrum, but not on its amplitude: a shape parameter 0.18 < \Gamma < 0.25 and an effective spectral index at 20Mpc/h in the range [-1.1,-0.9] are required by the Abell/ACO data. In order to obtain complementary constraints on the spectrum amplitude, we consider the cluster abundance as estimated using the Press--Schechter approach, whose reliability is explicitly tested against N--body simulations. We conclude that, of the cosmological models considered here, the only viable models are either Cold+Hot DM ones with \Omega_\nu = [0.2-0.3], better if shared between two massive neutrinos, and flat low-density CDM models with \Omega_0 = [0.3-0.5].

We investigate the effect of a putative X-ray emitting halo surrounding the Local Group of galaxies, and specifically the possible temperature anisotropies induced in the COBE-DMR four-year sky maps by an associated Sunyaev-Zel'dovich effect. By fitting the isothermal spherical halo model proposed by Suto et.al. (1996) to the coadded four-year COBE-DMR 53 and 90 GHz sky maps in Galactic coordinates, we find no significant evidence of a contribution. We therefore reject the claim that such a halo can affect the estimation of the primordial spectral index and amplitude of density perturbations as inferred from the DMR data. We find that correlation with the DMR data imposes constraints on the plausible contribution of such an X-ray emitting halo to a distortion in the CMB spectrum (as specified by the Compton-y parameter), up to a value for R -- the ratio of the core radius of the isothermal halo gas distribution to the distance to the Local Group centroid -- of 0.68. For larger values of R, the recent cosmological upper limit derived by COBE-FIRAS provides stronger constraints on the model parameters. Over the entire parameter space for R, we find an upper limit to the inferred sky-RMS anisotropy signal of 14 microKelvin (95% c.l.), a negligible amount relative to the 35 microKelvin signal observed in the COBE-DMR data.

Cut-sky orthogonal mode analyses of the COBE-DMR 53 and 90 GHz sky maps are used to determine the normalization of a variety of open cosmogonical models based on the cold dark matter scenario. To constrain the allowed cosmological-parameter range for these open cosmogonies, the predictions of the DMR-normalized models are compared to various observational measures of cosmography and large-scale structure, viz.: the age of the universe; small-scale dynamical estimates of the clustered-mass density parameter (Omega_0); constraints on the Hubble parameter (h), the X-ray cluster baryonic-mass fraction (Omega_B/Omega_0), and the matter power spectrum shape parameter; estimates of the mass perturbation amplitude; and constraints on the large-scale peculiar velocity field. The open-bubble inflation model is consistent with current determinations of the 95% confidence level range of these observational constraints.

Fundamental information about the Universe is encoded in anisotropies of the Cosmic Microwave Background (CMB) radiation. To make full use of this information, an experiment must image the entire sky with the angular resolution, sensitivity, and spectral coverage necessary to reach the limits set by cosmic variance on angular scales >~10'. Recent progress in detector technology allows this to be achieved by a properly designed space mission that fits well within the scope of NASA's Medium-class Explorer program. An essential component of the mission design is an observing strategy that minimizes systematic effects due to instrumental offset drifts. The detector advances make possible a `spin chopping' approach that has significant technical and scientific advantages over the strategy used by COBE, which reconstructed an image of the sky via inversion of a large matrix of differential measurements. The advantages include increased angular resolution, increased sensitivity, and simplicity of instrumentation and spacecraft operations. For the parameters typical of experiments like the Primordial Structures Investigation (PSI) and the Far InfraRed Explorer (FIRE), we show that the spin-chopping strategy produces images of the sky and power spectra of CMB anisotropies that contain no significant systematic artifacts.

The 2-point temperature correlation function is evaluated from the 4-year COBE DMR microwave anisotropy maps. We examine the 2-point function, which is the Legendre transform of the angular power spectrum, and show that the data are statistically consistent from channel to channel and frequency to frequency. The most likely quadrupole normalization is computed for a scale-invariant power-law spectrum of CMB anisotropy, using a variety of data combinations. For a given data set, the normalization inferred from the 2-point data is consistent with that inferred by other methods. The smallest and largest normalization deduced from any data combination are 16.4 and 19.6 uK respectively, with a value ~18 uK generally preferred.

We employ a pixel-based likelihood technique to estimate the angular power spectrum of the COBE Differential Microwave Radiometer (DMR) 4-year sky maps. The spectrum is consistent with a scale-invariant power-law form with a normalization, expressed in terms of the expected quadrupole anisotropy, of Q_{rms-PS|n=1} = 18 +/- 1.4 uK, and a best-fit spectral index of 1.2 +/- 0.3. The normalization is somewhat smaller than we concluded from the 2-year data, mainly due to additional Galactic modeling. We extend the analysis to investigate the extent to which the ``small" quadrupole observed in our sky is statistically consistent with a power-law spectrum. The most likely quadrupole amplitude is somewhat dependent on the details of Galactic foreground subtraction and data selection, ranging between 7 and 10 uK, but in no case is there compelling evidence that the quadrupole is too small to be consistent with a power-law spectrum. We conclude with a likelihood analysis of the band power amplitude in each of four spectral bands between l = 2 and 40, and find no evidence for deviations from a simple power-law spectrum.

The angular power spectrum estimator developed by Peebles (1973) and Hauser & Peebles (1973) has been modified and applied to the 4 year maps produced by the COBE DMR. The power spectrum of the observed sky has been compared to the power spectra of a large number of simulated random skies produced with noise equal to the observed noise and primordial density fluctuation power spectra of power law form, with $P(k) \propto k^n$. The best fitting value of the spectral index in the range of spatial scales corresponding to spherical harmonic indices $3 \leq \ell \lesssim 30$ is an apparent spectral index $n_{app}$ = 1.13 (+0.3) (-0.4) which is consistent with the Harrison-Zel'dovich primordial spectral index $n_{pri} = 1$ The best fitting amplitude for $n_{app} = 1$ is $\langle Q_{RMS}^2\rangle^{0.5}$ = 18 uK.

We use the COBE Differential Microwave Radiometers (DMR) 4-year sky maps to model Galactic microwave emission at high latitudes (|b| > 20 deg). Cross-correlation of the DMR maps with Galactic template maps detects fluctuations in the high-latitude microwave sky brightness with the angular variation of the DIRBE far-infrared dust maps and a frequency dependence consistent with a superposition of dust and free-free emission. We find no significant correlations between the DMR maps and various synchrotron templates. On the largest angular scales (e.g., quadrupole), Galactic emission is comparable in amplitude to the anisotropy in the cosmic microwave background (CMB). The CMB quadrupole amplitude, after correction for Galactic emission, has amplitude $Q_{rms}$ = 10.7 uK with random uncertainty 3.6 uK and systematic uncertainty 7.1 uK from uncertainty in our knowledge of Galactic microwave emission.

We search the high-latitude portion of the COBE Differential Microwave Radiometers (DMR) 4-year sky maps for evidence of a non-Gaussian temperature distribution in the cosmic microwave background. The genus, 3-point correlation function, and 2-point correlation function of temperature maxima and minima are all in excellent agreement with the hypothesis that the CMB anisotropy on angular scales of 7 degrees or larger represents a random-phase Gaussian field. A likelihood comparison of the DMR sky maps to a set of random-phase non-Gaussian toy models selects the exact Gaussian model as most likely. Monte Carlo simulations show that the 2-point correlation of the peaks and valleys in the maps provides the greatest discrimination among the class of models tested.

Fourier analysis and power spectrum estimation of the cosmic microwave background anisotropy on an incompletely sampled sky developed by Gorski (1994) has been applied to the high-latitude portion of the 4-year COBE DMR 31.5, 53 and 90 GHz sky maps. Likelihood analysis using newly constructed Galaxy cuts (extended beyond |b| = 20deg to excise the known foreground emission) and simultaneously correcting for the faint high latitude galactic foreground emission is conducted on the DMR sky maps pixelized in both ecliptic and galactic coordinates. The Bayesian power spectrum estimation from the foreground corrected 4-year COBE DMR data renders n ~ 1.2 +/- 0.3, and Q_{rms-PS} ~ 15.3^{+3.7}_{-2.8} microK (projections of the two-parameter likelihood). These results are consistent with the Harrison-Zel'dovich n=1 model of amplitude Q_{rms-PS} ~ 18 microK detected with significance exceeding 14sigma (dQ/Q < 0.07). (A small power spectrum amplitude drop below the published 2-year results is predominantly due to the application of the new, extended Galaxy cuts.)

We limit the possible contributions from non-cosmological sources to the COBE-DMR four-year sky maps. The DMR data are cross-correlated with maps of rich clusters, extragalactic IRAS sources, HEAO-1 A-2 X-ray emission and 5 GHz radio sources using a Fourier space technique. There is no evidence of significant contamination by such sources at an rms level of ~8 uK (95% confidence level at 7 degree resolution) in the most sensitive 53 GHz sky map. This level is consistent with previous limits set by analysis of earlier DMR data and by simple extrapolations from existing source models. We place a 95% confidence level rms limit on the Comptonization parameter averaged over the high-latitude sky of $\delta y\ <\ 1\ \times\ 10^{-6}$. Extragalactic sources have an insignificant effect on CMB power spectrum parameterizations determined from the DMR data.

The sky-RMS is the simplest model-independent characterization of a cosmological anisotropy signal. The RMS temperature fluctuations determined from the $COBE$-DMR four-year sky maps are frequency independent, consistent with the hypothesis that they are cosmological in origin, with a typical amplitude at 7 degrees of ~ 35 \pm\ 2 uK and at 10 degrees of ~29 \pm\ 1 uK. A joint analysis of the 7 deg and 10 deg "cross"-RMS derived from the data in both Galactic and Ecliptic coordinates is used to determine the RMS quadrupole normalization, $Q_{rms-PS}$, for a scale-invariant Harrison-Zel'dovich power law model. The difference in the inferred normalizations either including or excluding the quadrupole is a consequence of the sensitivity of the method to the Galaxy contaminated observed sky quadrupole. Whilst there are variations depending on the data selection, all results are consistent with an inferred \qrms\ normalization of ~18 \pm\ 2 uK.

The Differential Microwave Radiometers (DMR) instrument aboard the Cosmic Background Explorer (COBE) has mapped the full microwave sky to mean sensitivity 26 uK per 7 deg field of view. The absolute calibration is determined to 0.7% with drifts smaller than 0.2% per year. We have analyzed both the raw differential data and the pixelized sky maps for evidence of contaminating sources such as solar system foregrounds, instrumental susceptibilities, and artifacts from data recovery and processing. Most systematic effects couple only weakly to the sky maps. The largest uncertainties in the maps result from the instrument susceptibility to the Earth's magnetic field, microwave emission from the Earth, and upper limits to potential effects at the spacecraft spin period. Systematic effects in the maps are small compared to either the noise or the celestial signal: the 95% confidence upper limit for the pixel-pixel rms from all identified systematics is less than 6 uK in the worst channel. A power spectrum analysis of the (A-B)/2 difference maps shows no evidence for additional undetected systematic effects.

The cosmic microwave background radiation provides unique constraints on cosmological models. In this Letter we present a summary of the spatial properties of the cosmic microwave background radiation based on the full 4 years of COBE DMR observations, as detailed in a set of companion Letters. The anisotropy is consistent with a scale-invariant power law model and Gaussian statistics. With full use of the multi-frequency 4-year DMR data, including our estimate of the effects of Galactic emission, we find a power-law spectral index of $n=1.2\pm 0.3$ and a quadrupole normalization $Q_{rms-PS}=15.3^{+3.8}_{-2.8}$ $\mu$K. For $n=1$ the best-fit normalization is $Q_{rms-PS}\vert_{n=1}=18\pm 1.6$ $\mu$K. These values are consistent with both our previous 1-year and 2-year results. The results include use of the $\ell=2$ quadrupole term; exclusion of this term gives consistent results, but with larger uncertainties. The 4-year sky maps, presented in this Letter, portray an accurate overall visual impression of the anisotropy since the signal-to-noise ratio is ~2 per 10 degree sky map patch. The improved signal-to-noise ratio of the 4-year maps also allows for improvements in Galactic modeling and limits on non-Gaussian statistics.

We compute the CMB anisotropy in an open inflation CDM model which is normalized to the two-year DMR sky map. Presently available CMB data is consistent with $\Omega_0 \sim 0.3 - 0.4$, but does not strongly disfavour either $\Omega_0 \sim 0.1$ or $\Omega_0 \sim 1$. Accumulating CMB data should strengthen the observational discriminative power, however, a better understanding of systematic effects will be needed before it will be possible to draw robust conclusions about model viability.

We investigate the consequences of flat, dark-matter dominated cosmogonies with hybrid isocurvature and adiabatic initial perturbations and with Harrison-Zel'dovich primordial spectrum normalised to the $COBE$-DMR two-year measurements. We show that whilst the $COBE$-DMR data alone shows no preference for a specific admixture of these modes, acceptable combinations are strongly constrained by other observational data. Nevertheless, in some cases a suitable mixture of these modes still may be used in an attempt to avoid some of the observed problems of purely adiabatic models. Specifically, we consider critical density, cold dark matter (CDM) and mixed dark matter (MDM) models.

We cross-correlate the COBE-DMR 2-year sky maps with spatial templates from long-wavelength radio surveys and the far-infrared COBE DIRBE maps. We place an upper limit on the spectral index of synchrotron radiation beta_{synch} < -2.9 between 408 MHz and 31.5 GHz. We obtain a statistically significant cross-correlation with the DIRBE maps whose dependence on the DMR frequencies indicates a superposition of dust and free-free emission. The high-latitude dust emission (|b| > 30 deg) is well fitted by a single dust component with temperature T = 18 (+3, -7) K and emissivity epsilon ~ nu^beta with beta = 1.9 (+3.0, -0.5). The free-free emission is spatially correlated with the dust on angular scales larger than the 7 degree DMR beam, with rms variations 5.3 +/- 1.8 uK at 53 GHz and angular power spectrum P ~ ell^{-3}. If this correlation persists to smaller angular scales, free-free emission should not be a significant contaminant to measurements of the cosmic microwave anisotropy at degree angular scales for frequencies above 20 GHz.

The two-year COBE-DMR 53 and 90 GHz sky maps, in both galactic and ecliptic coordinates, are used to determine the normalisation of inflationary universe models with a flat global geometry and adiabatic density perturbations. The appropriately normalised cold and mixed dark matter models and cosmological constant dominated, cold dark matter models, computed for a range of values of Omega_b and h, are then compared to various measures of structure in the universe. Critical density CDM models appear to be irreconcilable with observations on both large and small scales simultaneously, whereas MDM models provide a somewhat better fit to the data. Although the COBE-DMR data alone prefer a nearly critical value for the total density, low-density cosmological constant models with Omega_0 greater than or equal to 0.15 can not be rejected at a confidence level exceeding 95%. Such models may also provide a significantly better fit to the matter distribution data than critical density CDM.

We compute the three-point temperature correlation function of the COBE Differential Microwave Radiometer (DMR) two-year sky maps to search for evidence of non-Gaussian temperature fluctuations. We detect three-point correlations in our sky with a substantially higher signal-to-noise ratio than from the first year data. However, the magnitude of the signal is consistent with the level of cosmic variance expected from Gaussian fluctuations, even when the low order multipole moments, up to l = 9, are filtered from the data. These results do not strongly constrain most existing models of structure formation, but the absence of intrinsic three-point correlations on large angular scales is an important consistency test for such models.

The three papers summarize the methods employed by the COBE-DMR group to analyze the power spectrum of the two-year DMR anisotropy maps. Results for the slope and normalization of the power spectrum are presented and compared. The implications of these results for various dark matter-based structure formation models are also discussed.

The two-year COBE-DMR 53 and 90 GHz sky maps are used to determine the normalization of inflationary, cosmological constant dominated, cold dark matter models. Although the COBE-DMR data prefer a nearly critical value for the total density, low-density models with Omega_0 > 0.15 can not be rejected at a confidence level > 95%. We compare theoretical predictions for various large scale statistics, computed for a range of values of Omega_0 and h=0.5, or 0.8, to various measures of structure in the universe.

Fourier analysis on a cut sky of the two-year COBE DMR 53 and 90 GHz sky maps (in galactic coordinates) is used to determine the normalisation of an open inflation model based on the cold dark matter scenario. The normalised model is compared to measures of large-scale structure in the universe. Although the DMR data alone does not provide sufficient discriminative power to prefer a particular value of the mass density parameter, the open model appears to be reasonably consistent with observations when Omega_0 ~ 0.3 - 0.4 and merits further study.

The two-year COBE-DMR 53 and 90 GHz sky maps in both ecliptic and galactic coordinates, are used to determine the normalisation of inflationary, flat, dark matter dominated universe models. The appropriately normalized cold and mixed dark matter models, computed for a range of values of Omega_baryon and h, and several values of the hot to cold dark matter abundance ratio, are then compared to various measures of structure in the universe.

The frequency-independent RMS temperature fluctuations determined from the COBE-DMR two year sky maps are used to infer the parameter Q_{rms-PS}, which characterizes the normalization of power law models of primordial cosmological temperature anisotropy. In particular, a 'cross'-RMS statistic is used to determine Q_{rms-PS} for a forced fit to a scale-invariant Harrison-Zel'dovich (n = 1) spectral model. Using a joint analysis of the 7 degree and 10 degree RMS temperature derived from both the 53 and 90 GHz sky maps, we find Q_{rms-PS} = 17.0^{+2.5}_{-2.1} uK when the low quadrupole is included, and Q_{rms-PS} = 19.4^{+2.3}_{-2.1} uK excluding the quadrupole. These results are consistent with the n = 1 fits from more sensitive methods (e.g. power spectrum, correlation function). The effect of the low quadrupole derived from the COBE-DMR data on the inferred Q_{rms-PS} normalization is investigated. A bias to lower Q_{rms-PS} is found when the quadrupole is included. The higher normalization for a forced n = 1 fit is then favored by the cross-RMS technique. As initially pointed out in Wright et al. (1994a) and further discussed here, analytic formulae for the RMS sky temperature fluctuations will NOT provide the correct normalization amplitude.

The first two years of COBE DMR observations of the CMB anisotropy are analyzed and compared with our previously published first year results. The results are consistent, but the addition of the second year of data increases the precision and accuracy of the detected CMB temperature fluctuations. The two-year 53 GHz data are characterized by RMS temperature fluctuations of DT=44+/-7 uK at 7 degrees and DT=30.5+/-2.7 uK at 10 degrees angular resolution. The 53X90 GHz cross-correlation amplitude at zero lag is C(0)^{1/2}=36+/-5 uK (68%CL) for the unsmoothed 7 degree DMR data. A likelihood analysis of the cross correlation function, including the quadrupole anisotropy, gives a most likely quadrupole-normalized amplitude Q_{rms-PS}=12.4^{+5.2}_{-3.3} uK (68% CL) and a spectral index n=1.59^{+0.49}_{-0.55} for a power law model of initial density fluctuations, P(k)~k^n. With n fixed to 1.0 the most likely amplitude is 17.4 +/-1.5 uK (68% CL). Excluding the quadrupole anisotropy we find Q_{rms-PS}= 16.0^{+7.5}_{-5.2} uK (68% CL), n=1.21^{+0.60}_{-0.55}, and, with n fixed to 1.0 the most likely amplitude is 18.2+/-1.6 uK (68% CL). Monte Carlo simulations indicate that these derived estimates of n may be biased by ~+0.3 (with the observed low value of the quadrupole included in the analysis) and {}~+0.1 (with the quadrupole excluded). Thus the most likely bias-corrected estimate of n is between 1.1 and 1.3. Our best estimate of the dipole from the two-year DMR data is 3.363+/-0.024 mK towards Galactic coordinates (l,b)= (264.4+/-0.2 degrees, +48.1+/-0.4 degrees), and our best estimate of the RMS quadrupole amplitude in our sky is 6+/-3 uK.

The natural approach to a spectral analysis of data distributed on the sky employs spherical harmonic decomposition. A common problem encountered in practical astronomy is the lack of full sky coverage in the available data. For example, the removal of the galactic plane data from the COBE DMR sky maps compromises Fourier analysis of the cosmic microwave background (CMB) temperature distribution due to the loss of orthogonality of the spherical harmonics. An explicit method for constructing orthonormal functions on an incomplete (e.g. Galaxy-cut) sphere is presented. These functions should be used in the proper Fourier analysis of the COBE DMR sky maps to provide the correct input for the determination of the spectrum of primordial inhomogeneity. The results of such an analysis are presented in an accompanying Letter. A similar algebraic construction of appropriate functions can be devised for other astronomical applications.

A new technique of Fourier analysis on a cut sky (Gorski, 1994) has been applied to the two year COBE DMR sky maps. The Bayesian power spectrum estimation results are consistent with the Harrison-Zel'dovich n=1 model. The maximum likelihood estimates of the parameters of the power spectrum of primordial perturbations are n=1.22 (1.02) and Q_{rms-PS}=17 (20) uK including (excluding) the quadrupole anisotropy. The marginal likelihood function on n renders n=1.10 \pm 0.32 (0.87 \pm 0.36).

We compute the three-point temperature correlation function of the {\it COBE} Differential Microwave Radiometer (DMR) first-year sky maps to search for non-Gaussian temperature fluctuations. The level of fluctuations seen in the computed correlation function are too large to be attributable solely to instrument noise. However the fluctuations are consistent with the level expected to result from a superposition of instrument noise and sky signal arising from a Gaussian power law model of initial fluctuations, with a quadrupole normalized amplitude of 17 $\mu$K and a power law spectral index $n = 1$. We place limits on the amplitude of intrinsic three-point correlations with a variety of predicted functional forms.