F. K. Hansen

HEALPIX : A framework for high-resolution discretization and fast analysis of data distributed on the sphere

HEALPix—the Hierarchical Equal Area isoLatitude Pixelization—is a versatile structure for the pixelization of data on the sphere. An associated library of computational algorithms and visualization software supports fast scientific applications executable directly on discretized spherical maps generated from very large volumes of astronomical data. Originally developed to address the data processing and analysis needs of the present generation of cosmic microwave background 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 inflation probe. In this paper we consider the requirements and implementation constraints on a framework that simultaneously enables an efficient discretization with associated hierarchical indexation and fast analysis/synthesis of functions defined on the sphere. We demonstrate how these are explicitly satisfied by HEALPix.

Master of the cosmic microwave background anisotropy power spectrum: a fast method for statistical analysis of large and complex cosmic microwave background data sets

We describe a fast and accurate method for estimation of the cosmic microwave background (CMB) anisotropy angular power spectrum—Monte Carlo Apodized Spherical Transform Estimator (MASTER). 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.

Asymmetries in the Cosmic Microwave Background Anisotropy Field

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 report on the results from two independent but complementary statistical analyses of the Wilkinson Microwave Anisotropy Probe (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°) 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 Differential Microwave Radiometer (DMR) map, and the axis of maximum asymmetry is close to the one found in the WMAP data.

Testing the cosmological principle of isotropy: Local power spectrum estimates of the WMAP data

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. In particular, we apply the transform with several apodisation scales, thus allowing the determination of the positional dependence of the angular power spectrum with either high spatial localisation or high angular resolution (ie. narrow bins in multipole space). Practically, this implies that we estimate the angular power spectrum locally in discs of various sizes positioned in different directions: small discs allow the greatest sensitivity to positional dependence, whereas larger discs allow greater sensitivity to variations over different angular scales. In addition, 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. Indeed, there is a remarkable lack of power in a region associated with the north ecliptic pole. With the greater fidelity in l-space allowed by this larger sky coverage, 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 (eg. planarity and alignment). However, over the range l = 5 40 the observed asymmetry is much harder to explain in terms of residual foregrounds and known systematic effects. By reorienting the coordinate axes, we partition the sky into different hemispheres and search for the reference frame which maximises the asymmetric distribution of power. The north pole for this coordinate frame is found to intersect the sphere at (80 ◦ ,57 ◦ ) in Galactic co-latitude and longitude over almost the entire multipole range l = 5 40. Furthermore, the strong negative outlier at l = 21 and the strong positive outlier at l = 39 as determined from the global power spectrum by the WMAP team, are found to be associated with the northern and southern hemispheres respectively (in this frame of maximum asymmetry). Thus, these two outliers follow the general tendency of the multipoles l = 5 40 to be of systematically lower amplitude in the north and higher in the south. Such asymmetric distributions of power on the sky provide a serious test for the cosmological principle of isotropy.

HEMISPHERICAL POWER ASYMMETRY IN THE THIRD-YEAR WILKINSON MICROWAVE ANISOTROPY PROBE SKY MAPS

We consider the issue of hemispherical power asymmetry in the third-year WMAP data, adopting a previously introduced modulation framework. Computing both frequentist probabilities and Bayesian evidences, we find that the model consisting 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 favor 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°, -27°), 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.

Component separation methods for the PLANCK mission

Context. The PLANCK satellite will map the full sky at nine frequencies from 30 to 857 GHz. The CMB intensity and polarization that are its prime targets are contaminated by foreground emission. Aims. The goal of this paper is to compare proposed methods for separating CMB from foregrounds based on their different spectral and spatial characteristics, and to separate the foregrounds into “components” with different physical origins (Galactic synchrotron, free-free and dust emissions; extra-galactic and far-IR point sources; Sunyaev-Zeldovich effect, etc.) Methods. A component separation challenge has been organised, based on a set of realistically complex simulations of sky emission. Several methods including those based on internal template subtraction, maximum entropy method, parametric method, spatial and harmonic cross correlation methods, and independent component analysis have been tested. Results. Different methods proved to be effective in cleaning the CMB maps of foreground contamination, in reconstructing maps of diffuse Galactic emissions, and in detecting point sources and thermal Sunyaev-Zeldovich signals. The power spectrum of the residuals is, on the largest scales, four orders of magnitude lower than the input Galaxy power spectrum at the foreground minimum. The CMB power spectrum was accurately recovered up to the sixth acoustic peak. The point source detection limit reaches 100 mJy, and about 2300 clusters are detected via the thermal SZ effect on two thirds of the sky. We have found that no single method performs best for all scientific objectives. Conclusions. We foresee that the final component separation pipeline for PLANCK will involve a combination of methods and iterations between processing steps targeted at different objectives such as diffuse component separation, spectral estimation, and compact source extraction.

POWER ASYMMETRY IN COSMIC MICROWAVE BACKGROUND FLUCTUATIONS FROM FULL SKY TO SUB-DEGREE SCALES: IS THE UNIVERSE ISOTROPIC?

We repeat and extend the analysis of Eriksen et al. and Hansen et al., testing the isotropy of the cosmic microwave background fluctuations. We find that the hemispherical power asymmetry previously reported for the largest scales l = 2-40 extends to much smaller scales. In fact, for the full multipole range l = 2-600, significantly more power is found in the hemisphere centered at (θ = 107° ± 10°, ∅ = 226° ± 10°) 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. 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.

Asymmetries in the Local Curvature of the Wilkinson Microwave Anisotropy Probe Data

We use the local curvature to investigate the possible existence of non-Gaussianity/asymmetry in the Wilkinson Microwave Anisotropy Probe data. Considering the full sky, we find results that are consistent with the Gaussian assumption. However, strong non-Gaussian features emerge when considering the northern and southern Galactic hemisphere separately, particularly on scales between 1° and 5°. Quite interestingly, the maximum non-Gaussianity is found for hemispheres centered near the ecliptic poles, which might suggest the presence of some systematic effect. The direction of the asymmetry seems consistent with the findings by Eriksen et al.We use the local curvature to investigate the possible existence of non-Gaussianity/asymmetry in the WMAP data. Considering the full sky we find results which are consistent with the Gaussian assumption. However, strong non-Gaussian features emerge when considering the northern and southern galactic hemisphere separately, particularly on scales between 1 and 5 degrees. Quite interestingly, the maximum non-Gaussianity is found for hemispheres centered near the ecliptic poles, which might suggest the presence of some systematic effect. The direction of the asymmetry seems consistent with the findings by Eriksen et al. 2004. Subject headings: (cosmology:) cosmic microwave background — cosmology: observations — methods: data analysis — methods: statistical

Cosmic microwave background component separation by parameter estimation

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.We propose a method for CMB component separation based on standard Bayesian 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 approximated by a χ2 statistic. Then we fix all nonlinear 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 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 satellite and 6 year WMAP data based on realistic models and show that separation at the microkelvin 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.

Asymmetries in the CMB anisotropy field

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 report on the results from two independent but complementary statistical analyses of the Wilkinson Microwave Anisotropy Probe (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°) 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 Differential Microwave Radiometer (DMR) map, and the axis of maximum asymmetry is close to the one found in the WMAP data.