Showing 201–250 of 305 results for author: Banday, A J

Planck Early Results XVIII: The power spectrum of cosmic infrared background anisotropies

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. Balbi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, E. Battaner, K. Benabed, A. Benoît, J. -P. Bernard, M. Bersanelli, R. Bhatia, K. Blagrave, J. J. Bock, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, M. Bucher , et al. (180 additional authors not shown)

Abstract: Using Planck maps of six regions of low Galactic dust emission with a total area of about 140 square degrees, we determine the angular power spectra of cosmic infrared background (CIB) anisotropies from multipole l = 200 to l = 2000 at 217, 353, 545 and 857 GHz. We use 21-cm observations of HI as a tracer of thermal dust emission to reduce the already low level of Galactic dust emission and use th… ▽ More Using Planck maps of six regions of low Galactic dust emission with a total area of about 140 square degrees, we determine the angular power spectra of cosmic infrared background (CIB) anisotropies from multipole l = 200 to l = 2000 at 217, 353, 545 and 857 GHz. We use 21-cm observations of HI as a tracer of thermal dust emission to reduce the already low level of Galactic dust emission and use the 143 GHz Planck maps in these fields to clean out cosmic microwave background anisotropies. Both of these cleaning processes are necessary to avoid significant contamination of the CIB signal. We measure correlated CIB structure across frequencies. As expected, the correlation decreases with increasing frequency separation, because the contribution of high-redshift galaxies to CIB anisotropies increases with wavelengths. We find no significant difference between the frequency spectrum of the CIB anisotropies and the CIB mean, with Delta I/I=15% from 217 to 857 GHz. In terms of clustering properties, the Planck data alone rule out the linear scale- and redshift-independent bias model. Non-linear corrections are significant. Consequently, we develop an alternative model that couples a dusty galaxy, parametric evolution model with a simple halo-model approach. It provides an excellent fit to the measured anisotropy angular power spectra and suggests that a different halo occupation distribution is required at each frequency, which is consistent with our expectation that each frequency is dominated by contributions from different redshifts. In our best-fit model, half of the anisotropy power at l=2000 comes from redshifts z<0.8 at 857 GHz and z<1.5 at 545 GHz, while about 90% come from redshifts z>2 at 353 and 217 GHz, respectively. △ Less

Submitted 17 August, 2011; v1 submitted 11 January, 2011; originally announced January 2011.

Comments: Final A&A version. Few typos corrected

Journal ref: Astronomy & Astrophysics, 2011, Volume 536, id.A18

Planck early results: Cluster Sunyaev-Zeldovich optical scaling relations

Authors: Planck Collaboration, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. Balbi, A. J. Banday, R. B. Barreiro, M. Bartelmann, J. G. Bartlett, E. Battaner, K. Benabed, A. Benoît, J. -P. Bernard, M. Bersanelli, R. Bhatia, J. J. Bock, A. Bonaldi, J. R. Bond, J. Borrill, F. R. Bouchet, M. L. Brown, M. Bucher, C. Burigana , et al. (180 additional authors not shown)

Abstract: We present the Sunyaev-Zeldovich (SZ) signal-to-richness scaling relation (Y500-N200) for the MaxBCG cluster catalogue. Employing a multi-frequency matched filter on the Planck sky maps, we measure the SZ signal for each cluster by adapting the filter according to weak-lensing calibrated mass-richness relations (N200-M500). We bin our individual measurements and detect the SZ signal down to the lo… ▽ More We present the Sunyaev-Zeldovich (SZ) signal-to-richness scaling relation (Y500-N200) for the MaxBCG cluster catalogue. Employing a multi-frequency matched filter on the Planck sky maps, we measure the SZ signal for each cluster by adapting the filter according to weak-lensing calibrated mass-richness relations (N200-M500). We bin our individual measurements and detect the SZ signal down to the lowest richness systems (N200=10) with high significance, achieving a detection of the SZ signal in systems with mass as low as M500~5e13 Msolar. The observed Y500-N200 relation is well modeled by a power law over the full richness range. It has a lower normalisation at given N200 than predicted based on X-ray models and published mass-richness relations. An X-ray subsample, however, does conform to the predicted scaling, and model predictions do reproduce the relation between our measured bin-average SZ signal and measured bin-average X-ray luminosities. At fixed richness, we find an intrinsic dispersion in the Y500-N200 relation of 60% rising to of order 100% at low richness. Thanks to its all-sky coverage, Planck provides observations for more than 13,000 MaxBCG clusters and an unprecedented SZ/optical data set, extending the list of known cluster scaling laws to include SZ-optical properties. The data set offers essential clues for models of galaxy formation. Moreover, the lower normalisation of the SZ-mass relation implied by the observed SZ-richness scaling has important consequences for cluster physics and cosmological studies with SZ clusters. △ Less

Submitted 11 January, 2011; originally announced January 2011.

Comments: Also available from this URL: http://www.rssd.esa.int

Planck Early Results XI: Calibration of the local galaxy cluster Sunyaev-Zeldovich scaling relations

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. Balbi, A. J. Banday, R. B. Barreiro, M. Bartelmann, J. G. Bartlett, E. Battaner, K. Benabed, A. Benoît, J. -P. Bernard, M. Bersanelli, R. Bhatia, J. J. Bock, A. Bonaldi, J. R. Bond, J. Borrill, F. R. Bouchet, H. Bourdin, M. L. Brown , et al. (185 additional authors not shown)

Abstract: We present precise Sunyaev-Zeldovich (SZ) effect measurements in the direction of 62 nearby galaxy clusters (z <0.5) detected at high signal-to-noise in the first Planck all-sky dataset. The sample spans approximately a decade in total mass, 10^14 < M_500 < 10^15, where M_500 is the mass corresponding to a total density contrast of 500. Combining these high quality Planck measurements with deep XM… ▽ More We present precise Sunyaev-Zeldovich (SZ) effect measurements in the direction of 62 nearby galaxy clusters (z <0.5) detected at high signal-to-noise in the first Planck all-sky dataset. The sample spans approximately a decade in total mass, 10^14 < M_500 < 10^15, where M_500 is the mass corresponding to a total density contrast of 500. Combining these high quality Planck measurements with deep XMM-Newton X-ray data, we investigate the relations between D_A^2 Y_500, the integrated Compton parameter due to the SZ effect, and the X-ray-derived gas mass M_g,500, temperature T_X, luminosity L_X, SZ signal analogue Y_X,500 = M_g,500 * T_X, and total mass M_500. After correction for the effect of selection bias on the scaling relations, we find results that are in excellent agreement with both X-ray predictions and recently-published ground-based data derived from smaller samples. The present data yield an exceptionally robust, high-quality local reference, and illustrate Planck's unique capabilities for all-sky statistical studies of galaxy clusters. △ Less

Submitted 15 September, 2011; v1 submitted 11 January, 2011; originally announced January 2011.

Comments: 14 pages, 6 figures, to appear in A&A; final version corrects values for A1413 and A3921 in Table 1

Journal ref: Astron.Astrophys. 536 (2011) A11

Planck early results. IX. XMM-Newton follow-up for validation of Planck cluster candidates

Authors: Planck Collaboration, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. Balbi, A. J. Banday, R. B. Barreiro, M. Bartelmann, J. G. Bartlett, E. Battaner, K. Benabed, A. Benoît, J. -P. Bernard, M. Bersanelli, R. Bhatia, J. J. Bock, A. Bonaldi, J. R. Bond, J. Borrill, F. R. Bouchet, M. L. Brown, M. Bucher, C. Burigana , et al. (175 additional authors not shown)

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

Submitted 27 September, 2011; v1 submitted 11 January, 2011; originally announced January 2011.

Comments: 21 pages, 16 figures (some low resolution), final version accepted for publication in A&A

Planck Early Results VIII: The all-sky Early Sunyaev-Zeldovich cluster sample

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. Balbi, A. J. Banday, R. B. Barreiro, M. Bartelmann, J. G. Bartlett, E. Battaner, R. Battye, K. Benabed, A. Benoît, J. -P. Bernard, M. Bersanelli, R. Bhatia, J. J. Bock, A. Bonaldi, J. R. Bond, J. Borrill, F. R. Bouchet, M. L. Brown , et al. (213 additional authors not shown)

Abstract: We present the first all-sky sample of galaxy clusters detected blindly by the Planck satellite through the Sunyaev-Zeldovich (SZ) effect from its six highest frequencies. This early SZ (ESZ) sample is comprised of 189 candidates, which have a high signal-to-noise ratio ranging from 6 to 29. Its high reliability (purity above 95%) is further ensured by an extensive validation process based on Plan… ▽ More We present the first all-sky sample of galaxy clusters detected blindly by the Planck satellite through the Sunyaev-Zeldovich (SZ) effect from its six highest frequencies. This early SZ (ESZ) sample is comprised of 189 candidates, which have a high signal-to-noise ratio ranging from 6 to 29. Its high reliability (purity above 95%) is further ensured by an extensive validation process based on Planck internal quality assessments and by external cross-identification and follow-up observations. Planck provides the first measured SZ signal for about 80% of the 169 previously-known ESZ clusters. Planck furthermore releases 30 new cluster candidates, amongst which 20 meet the ESZ signal-to-noise selection criterion. At the submission date, twelve of the 20 ESZ candidates were confirmed as new clusters, with eleven confirmed using XMM-Newton snapshot observations, most of them with disturbed morphologies and low luminosities. The ESZ clusters are mostly at moderate redshifts (86% with z below 0.3) and span more than a decade in mass, up to the rarest and most massive clusters with masses above 10^15 Msol. △ Less

Submitted 7 July, 2011; v1 submitted 11 January, 2011; originally announced January 2011.

Comments: Accepted for publication in a special issue of A&A, 26 pages, this paper is part of a package of papers describing first results of the Planck mission

Planck Early Results. II. The thermal performance of Planck

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, M. Baker, A. Balbi, A. J. Banday, R. B. Barreiro, E. Battaner, K. Benabed, A. Benoit, J. P. Bernard, M. Bersanelli, P. Bhandari, R. Bhatia, J. J. Bock, A. Bonaldi, J. R. Bond, J. Borders, J. Borrill, F. R. Bouchet, B. Bowman , et al. (203 additional authors not shown)

Abstract: 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 satisf… ▽ More 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. △ Less

Submitted 2 January, 2012; v1 submitted 11 January, 2011; originally announced January 2011.

Comments: Accepted by A&A. Also available at http://www.rssd.esa.int

Report number: Planck Collaboration 2011b. / Planck2011-1.3

Planck Early Results: The Planck mission

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, M. Baker, A. Balbi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, E. Battaner, K. Benabed, K. Bennett, A. Benoît, J. -P. Bernard, M. Bersanelli, R. Bhatia, J. J. Bock, A. Bonaldi, J. R. Bond, J. Borrill, F. R. Bouchet, T. Bradshaw , et al. (250 additional authors not shown)

Abstract: The European Space Agency's Planck satellite was launched on 14 May 2009, and has been surveying the sky stably and continuously since 13 August 2009. Its performance is well in line with expectations, and it will continue to gather scientific data until the end of its cryogenic lifetime. We give an overview of the history of Planck in its first year of operations, and describe some of the key per… ▽ More The European Space Agency's Planck satellite was launched on 14 May 2009, and has been surveying the sky stably and continuously since 13 August 2009. Its performance is well in line with expectations, and it will continue to gather scientific data until the end of its cryogenic lifetime. We give an overview of the history of Planck in its first year of operations, and describe some of the key performance aspects of the satellite. This paper is part of a package submitted in conjunction with Planck's Early Release Compact Source Catalogue, the first data product based on Planck to be released publicly. The package describes the scientific performance of the Planck payload, and presents results on a variety of astrophysical topics related to the sources included in the Catalogue, as well as selected topics on diffuse emission. △ Less

Submitted 16 June, 2011; v1 submitted 11 January, 2011; originally announced January 2011.

Comments: This is part of a package of Planck papers labelled in their titles as "Planck Early Results". The whole package can also be downloaded from http://www.rssd.esa.int/Planck. This paper was accepted by Astronomy & Astrophysics on 31 May 2011

Planck early results XIV: ERCSC validation and extreme radio sources

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, E. Angelakis, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. Balbi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, E. Battaner, K. Benabed, A. Benoît, J. P. Bernard, M. Bersanelli, R. Bhatia, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, M. Bucher, C. Burigana , et al. (184 additional authors not shown)

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

Submitted 16 December, 2011; v1 submitted 10 January, 2011; originally announced January 2011.

Comments: Revised to accommodate referee's comments. Changes made to sections 5.1 and 5.4. No change in figures or scientific conclusions

Journal ref: 2011 A&A 536

Constraints on the topology of the Universe derived from the 7-year WMAP data

Authors: P. Bielewicz, A. J. Banday

Abstract: We impose constraints on the topology of the Universe determined from a search for matched circles in the 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, and show that for a full sky estimate of the CMB signal (the ILC map) such residuals introduce a non-negligible effect on the… ▽ More We impose constraints on the topology of the Universe determined from a search for matched circles in the 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, and show that for a full sky estimate of the CMB signal (the ILC map) such residuals introduce a non-negligible effect on the statistics of matched circles. In order to reduce this effect, we perform the analysis on maps for which the most contaminated regions have been removed. 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 α_min \approx 10 degrees. This places a lower bound on the size of the fundamental domain for a flat universe of about 27.9 Gpc. This bound is close to the upper limit on the size of Universe possible to detect by the method of matched circles, i.e. the diameter of the observable Universe is 28.3 Gpc. △ Less

Submitted 6 February, 2018; v1 submitted 16 December, 2010; originally announced December 2010.

Comments: 8 pages, 3 figures, published in MNRAS, this research was supported by the Agence Nationale de Recherche (ANR-08-CEXC-0002-01)

Journal ref: MNRAS, 412 (2011), 2104

Scale-dependent non-Gaussianities in the WMAP data as identified by using surrogates and scaling indices

Authors: C. Raeth, A. J. Banday, G. Rossmanith, H. Modest, R. Suetterlin, K. M. Gorski, J. Delabrouille, G. E. Morfill

Abstract: 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 i… ▽ More 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. △ Less

Submitted 27 April, 2011; v1 submitted 14 December, 2010; originally announced December 2010.

Comments: 11 pages, 6 figures, matches MNRAS accepted version

The primordial non-Gaussianity of local type (f_NL) in the WMAP 5-year data: the length distribution of CMB skeleton

Authors: Zhen Hou, A. J. Banday, Krzysztof M. Gorski, Franz Elsner, Benjamin D. Wandelt

Abstract: 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-Gau… ▽ More 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. △ Less

Submitted 1 June, 2010; v1 submitted 30 May, 2010; originally announced May 2010.

Comments: 17 pages, 14 figures, 5 tables, accepted for publication in MNRAS

Probing non-Gaussianities on Large Scales in WMAP5 and WMAP7 Data using Surrogates

Authors: C. Raeth, G. Rossmanith, G. Morfill, A. J. Banday, K. M. Gorski

Abstract: 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… ▽ More 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. △ Less

Submitted 14 May, 2010; originally announced May 2010.

Comments: 4 pages, 3 figures, to appear in the Proceedings of Moriond Cosmology 2010

Joint 3D modelling of the polarized Galactic synchrotron and thermal dust foreground diffuse emission

Authors: L. Fauvet, J. F. Macías-Pérez, J. Aumont, F. X. Désert, T. R. Jaffe, A. J. Banday, M. Tristram, A. H. Waelkens, D. Santos

Abstract: We present for the first time a coherent model of the polarized Galactic synchrotron and thermal dust emissions which are the main diffuse foreground for the measurement of the polarized power spectra of the CMB fluctuations with the Planck satellite mission. We produce 3D models of the Galactic magnetic field including regular and turbulent components, and of the distribution of matter in the Gal… ▽ More We present for the first time a coherent model of the polarized Galactic synchrotron and thermal dust emissions which are the main diffuse foreground for the measurement of the polarized power spectra of the CMB fluctuations with the Planck satellite mission. We produce 3D models of the Galactic magnetic field including regular and turbulent components, and of the distribution of matter in the Galaxy, relativistic electrons and dust grains. By integrating along the line of sight we construct maps of the polarized Galactic synchrotron and thermal dust emission for each of these models and compare them to currently available data. We consider the 408 MHz all-sky continuum survey, the 23 GHz band of the Wilkinson Microwave Anisotropy Probe and the 353 GHz Archeops data.}{The best-fit parameters obtained are consistent with previous estimates in the literature based only on synchrotron emission and pulsar rotation measurements. They allows us to reproduce the large scale structures observed on the data. Poorly understood local Galactic structures and turbulence make difficult an accurate reconstruction of the observations in the Galactic plane. Finally, using the best-fit model we are able to estimate the expected polarized foreground contamination at the Planck frequency bands. For the CMB bands, 70, 100, 143 and 217 GHz, at high Galactic latitudes although the CMB signal dominates in general, a significant foreground contribution is expected at large angular scales. In particular, this contribution will dominate the CMB signal for the B modes expected from realistic models of a background of primordial gravitational waves. △ Less

Submitted 23 March, 2010; originally announced March 2010.

Journal ref: A&A, 2010

Planck pre-launch status: the Planck-LFI programme

Authors: N. Mandolesi, M. Bersanelli, R. C. Butler, E. Artal, C. Baccigalupi, A. Balbi, A. J. Banday, R. B. Barreiro, M. Bartelmann, K. Bennett, P. Bhandari, A. Bonaldi, J. Borrill, M. Bremer, C. Burigana, R. C. Bowman, P. Cabella, C. Cantalupo, B. Cappellini, T. Courvoisier, G. Crone, F. Cuttaia, L. Danese, O. D'Arcangelo, R. D. Davies , et al. (118 additional authors not shown)

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

Submitted 15 January, 2010; originally announced January 2010.

Comments: 25 pages, 16 figures. In press on Astronomy and Astrophysics

The two-and three-point correlation functions of the polarized five-year WMAP sky maps

Authors: E. Gjerløw, H. K. Eriksen, A. J. Banday, K. M. Górski, P. B. Lilje

Abstract: We present the two- and three-point real space correlation functions of the five-year WMAP sky maps, and compare the observed functions to simulated LCDM concordance model ensembles. In agreement with previously published results, we find that the temperature correlation functions are consistent with expectations. However, the pure polarization correlation functions are acceptable only for the 3… ▽ More We present the two- and three-point real space correlation functions of the five-year WMAP sky maps, and compare the observed functions to simulated LCDM concordance model ensembles. In agreement with previously published results, we find that the temperature correlation functions are consistent with expectations. However, the pure polarization correlation functions are acceptable only for the 33GHz band map; the 41, 61, and 94 GHz band correlation functions all exhibit significant large-scale excess structures. Further, these excess structures very closely match the correlation functions of the two (synchrotron and dust) foreground templates used to correct the WMAP data for galactic contamination, with a cross-correlation statistically significant at the 2sigma-3sigma confidence level. The correlation is slightly stronger with respect to the thermal dust template than with the synchrotron template. △ Less

Submitted 23 January, 2010; v1 submitted 30 October, 2009; originally announced October 2009.

Comments: 10 pages, 5 figures, published in ApJ. v2: New title, minor changes to appendix, and fixed some typos. v3: Matches version published in ApJ

Journal ref: Astrophysical Journal 710 (2010) 689-697

New insights into foreground analysis of the WMAP five-year data using FASTICA

Authors: M. Bottino, A. J. Banday, D. Maino

Abstract: In this paper, we present a foreground analysis of the WMAP 5-year data using the FASTICA algorithm, improving on the treatment of the WMAP 3-year data in Bottino et al 2008. We revisit the nature of the free-free spectrum with the emphasis on attempting to confirm or otherwise the spectral feature claimed in Dobbler et al 2008b and explained in terms of spinning dust emission in the warm ionise… ▽ More In this paper, we present a foreground analysis of the WMAP 5-year data using the FASTICA algorithm, improving on the treatment of the WMAP 3-year data in Bottino et al 2008. We revisit the nature of the free-free spectrum with the emphasis on attempting to confirm or otherwise the spectral feature claimed in Dobbler et al 2008b and explained in terms of spinning dust emission in the warm ionised medium. With the application of different Galactic cuts, the index is always flatter than the canonical value of 2.14 except for the Kp0 mask which is steeper. Irrespective of this, we can not confirm the presence of any feature in the free-free spectrum. We experiment with a more extensive approach to the cleaning of the data, introduced in connection with the iterative application of FASTICA. We confirm the presence of a residual foreground whose spatial distribution is concentrated along the Galactic plane, with pronounced emission near the Galactic center. This is consistent with the WMAP haze detected in Finkbeiner 2004. Finally, we attempted to perform the same analysis on full-sky maps. The code returns good results even for those regions where the cross-talk among the components is high. However, slightly better results in terms of the possibility of reconstructing a full-sky CMB map, are achieved with a simultaneous analysis of both the five WMAP maps and foreground templates. Nonetheless, some residuals are still present and detected in terms of an excess in the CMB power spectrum, on small angular scales. Therefore, a minimal mask for the brightest regions of the plane is necessary, and has been defined. △ Less

Submitted 27 October, 2009; v1 submitted 22 October, 2009; originally announced October 2009.

Comments: Accepted for publication in MNRAS, 25 pages, 17 figures, 4 tables. Version with full resolution figures available at: http://www.mpa-garching.mpg.de/~bottino/downloads/bottino_etal.pdf

Frequentist comparison of CMB local extrema statistics in the five-year WMAP data with two anisotropic cosmological models

Authors: Zhen Hou, A. J. Banday, K. M. Gorski, N. E. Groeneboom, H. K. Eriksen

Abstract: 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 re… ▽ More 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. △ Less

Submitted 18 October, 2009; originally announced October 2009.

Comments: 10 pages, 4 figures, 3 tables, accepted for publication in MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.401:2379-2387, 2010

Modelling the Galactic Magnetic Field on the Plane in 2D

Authors: T. R. Jaffe, J. P. Leahy, A. J. Banday, S. M. Leach, S. R. Lowe, A. Wilkinson

Abstract: We present a method for parametric modelling of the physical components of the Galaxy's magnetised interstellar medium, simulating the observables, and mapping out the likelihood space using a Markov Chain Monte-Carlo analysis. We then demonstrate it using total and polarised synchrotron emission data as well as rotation measures of extragalactic sources. With these three datasets, we define and… ▽ More We present a method for parametric modelling of the physical components of the Galaxy's magnetised interstellar medium, simulating the observables, and mapping out the likelihood space using a Markov Chain Monte-Carlo analysis. We then demonstrate it using total and polarised synchrotron emission data as well as rotation measures of extragalactic sources. With these three datasets, we define and study three components of the magnetic field: the large-scale coherent field, the small-scale isotropic random field, and the ordered field. In this first paper, we use only data along the Galactic plane and test a simple 2D logarithmic spiral model for the magnetic field that includes a compression and a shearing of the random component giving rise to an ordered component. We demonstrate with simulations that the method can indeed constrain multiple parameters yielding measures of, for example, the ratios of the magnetic field components. Though subject to uncertainties in thermal and cosmic ray electron densities and depending on our particular model parametrisation, our preliminary analysis shows that the coherent component is a small fraction of the total magnetic field and that an ordered component comparable in strength to the isotropic random component is required to explain the polarisation fraction of synchrotron emission. We outline further work to extend this type of analysis to study the magnetic spiral arm structure, the details of the turbulence as well as the 3D structure of the magnetic field. △ Less

Submitted 29 January, 2010; v1 submitted 23 July, 2009; originally announced July 2009.

Comments: 18 pages, 11 figures, updated to published MNRAS version

Non-Gaussian Signatures in the five-year WMAP data as identified with isotropic scaling indices

Authors: G. Rossmanith, C. Raeth, A. J. Banday, G. Morfill

Abstract: We continue the analysis of non-Gaussianities in the CMB by means of the scaling index method (SIM, Raeth, Schuecker & Banday 2007) by applying this method on the 5-year WMAP data. We compare each of the results with 1000 Monte Carlo simulations mimicing the Gaussian properties of the best fit $ΛCDM$-model. Based on the scaling indices, scale-dependent empirical probability distributions, moment… ▽ More We continue the analysis of non-Gaussianities in the CMB by means of the scaling index method (SIM, Raeth, Schuecker & Banday 2007) by applying this method on the 5-year WMAP data. We compare each of the results with 1000 Monte Carlo simulations mimicing the Gaussian properties of the best fit $ΛCDM$-model. Based on the scaling indices, scale-dependent empirical probability distributions, moments of these distributions and $χ^2$-combinations of them are calculated, obtaining similar results as in the former analysis of the 3-year data: We derive evidence for non-Gaussianity with a probability of up to 97.3% for the mean when regarding the KQ75-masked full sky and summing up over all considered length scales by means of a diagonal $χ^2$-statistics. Looking at only the northern or southern hemisphere, we obtain up to 98.5% or 96.6%, respectively. For the standard deviation, these results appear as 95.6% for the full sky (99.7% north, 89.4% south) and for a $χ^2$-combination of both measurements as 97.4% (99.1% north, 95.5% south). By performing an analysis of rotated hemispheres, we detect an obvious asymmetry in the data. In addition to these investigations, we present a method of filling the mask with Gaussian noise to eliminate boundary effects caused by the mask. With the help of this technique, we identify several local features on the map, of which the most significant one turns out to be the well-known cold spot. When excluding all these spots from the analysis, the deviation from Gaussianity increases, which shows that the discovered local anomalies are not the reason of the global detection of non-Gaussianity, but actually were damping the deviations on average. Our analyses per band and per year suggest, however, that it is very unlikely that the detected anomalies are due to foreground effects. △ Less

Submitted 21 July, 2009; v1 submitted 18 May, 2009; originally announced May 2009.

Comments: 14 pages, 10 figures; minor revisions; accepted for publication by MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.399:1921-1933,2009

The Hot and Cold Spots in Five-Year WMAP Data

Authors: Zhen Hou, A. J. Banday, K. M. Gorski

Abstract: 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-va… ▽ More 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. △ Less

Submitted 25 March, 2009; originally announced March 2009.

Comments: 15 pages, 6 figures, 8 tables, accepted for publication in MNRAS

Journal ref: Mon.Not.R.Astron.Soc. 396 (2009) 1273-1286

Bayesian component separation and CMB estimation for the 5-year WMAP temperature data

Authors: C. Dickinson, H. K. Eriksen, A. J. Banday, J. B. Jewell, K. M. Gorski, G. Huey, C. R. Lawrence, I. J. O'Dwyer, B. D. Wandelt

Abstract: 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-… ▽ More 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 μ$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 ($β=-2.97\pm0.21$) relative to those found at low ($\sim$GHz) frequencies. △ Less

Submitted 25 March, 2009; originally announced March 2009.

Comments: 18 pages, 19 figures, submitted to ApJ

Journal ref: Astrophys.J.705:1607-1623,2009

Increasing evidence for hemispherical power asymmetry in the five-year WMAP data

Authors: J. Hoftuft, H. K. Eriksen, A. J. Banday, K. M. Gorski, F. K. Hansen, P. B. Lilje

Abstract: (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 modulat… ▽ More (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). △ Less

Submitted 20 April, 2009; v1 submitted 6 March, 2009; originally announced March 2009.

Comments: 6 pages, 2 figures; added references and minor comments. Accepted for publication in ApJ

Journal ref: Astrophys.J.699:985-989,2009

Power Asymmetry in Cosmic Microwave Background Fluctuations from Full Sky to Sub-degree Scales: Is the Universe Isotropic?

Authors: F. K. Hansen, A. J. Banday, K. M. Gorski, H. K. Eriksen, P. B. Lilje

Abstract: 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 (thet… ▽ More 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. △ Less

Submitted 19 December, 2008; originally announced December 2008.

Comments: Submitted to ApJ

Journal ref: Astrophys.J.704:1448-1458,2009

A model-independent test for scale-dependent non-Gaussianities in the CMB

Authors: C. Raeth, G. Morfill, G. Rossmanith, A. J. Banday, K. M. Gorski

Abstract: 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 meth… ▽ More 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. △ Less

Submitted 18 March, 2009; v1 submitted 21 October, 2008; originally announced October 2008.

Comments: accepted for publication in PRL, minor revisions, results unchanged, l(cut)-dependency added

Journal ref: Phys. Rev. Lett., 102, 131301 (2009)

Foreground analysis of the WMAP three-year data with FASTICA

Authors: M. Bottino, A. J. Banday, D. Maino

Abstract: We present an analysis of the foreground emission present in the WMAP 3-year data as determined by the method of Independent Component Analysis. We derived coupling coefficients between the WMAP data and foreground templates which are then used to infer the spectral behaviour for three foreground components -- synchrotron, anomalous dust-correlated emission and free-free. For the first two compo… ▽ More We present an analysis of the foreground emission present in the WMAP 3-year data as determined by the method of Independent Component Analysis. We derived coupling coefficients between the WMAP data and foreground templates which are then used to infer the spectral behaviour for three foreground components -- synchrotron, anomalous dust-correlated emission and free-free. For the first two components, we find values consistent with previous results although slightly steeper. We confirm the inconsistency in the scaling between the Ha template and free-free emission at K- and Ka-bands where an electron temperature of ~ 4000 K is indicated. We also see evidence of significantly flatter spectral behaviour to higher frequencies than expected theoretically and previously noted by Dobler et al.(2008a), but only when analysing the Kp2 sky coverage. We further apply FASTICA 'iteratively', using data pre-cleaned using foreground templates scaled to the WMAP frequencies by coupling coefficients determined by a prior FASTICA analysis. This multi-frequency analysis allows us to determine the presence of residual foreground emission not traced by the templates. We confirm the existence of a component spatially distributed along the Galactic plane and particularly enhanced near the center (the 'WMAP haze'). This emission is less extended when using the WMAP K-Ka data as the synchrotron template confirming that it can be considered a better template for foreground cleaning of the WMAP data. However its use complicates the physical interpretation of the nature of the foreground emission and residuals. since it contains a mixture of several, physically distinct emission mechanisms. △ Less

Submitted 11 July, 2008; originally announced July 2008.

Comments: Accepted for publication in MNRAS, 24 pages, 14 figures

The scalar perturbation spectral index n_s: WMAP sensitivity to unresolved point sources

Authors: K. M. Huffenberger, H. K. Eriksen, F. K. Hansen, A. J. Banday, K. M. Gorski

Abstract: 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 sour… ▽ More 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. △ Less

Submitted 10 October, 2007; originally announced October 2007.

Comments: 9 pages, 10 figures, 3 tables. Submitted to ApJ

Joint Bayesian component separation and CMB power spectrum estimation

Authors: H. K. Eriksen, J. B. Jewell, C. Dickinson, A. J. Banday, K. M. Gorski, C. R. Lawrence

Abstract: 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 t… ▽ More 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. △ Less

Submitted 5 November, 2007; v1 submitted 7 September, 2007; originally announced September 2007.

Comments: 23 pages, 16 figures; version accepted for publication in ApJ -- only minor changes, all clarifications. More information about the WMAP3 analysis available at http://www.astro.uio.no/~hke under the Research tab

Journal ref: Astrophys.J.676:10-32,2008

The joint large-scale foreground-CMB posteriors of the 3-year WMAP data

Authors: H. K. Eriksen, C. Dickinson, J. B. Jewell, A. J. Banday, K. M. Gorski, C. R. Lawrence

Abstract: 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 fr… ▽ More 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. △ Less

Submitted 7 September, 2007; originally announced September 2007.

Comments: 5 pages, 4 figures, submitted to ApJL. Background data are available at http://www.astro.uio.no/~hke under the Research tab

Journal ref: Astrophys.J.672:L87-L90,2008

Bayesian analysis of the low-resolution polarized 3-year WMAP sky maps

Authors: H. K. Eriksen, Greg Huey, A. J. Banday, K. M. Gorski, J. B. Jewell, I. J. O'Dwyer, B. D. Wandelt

Abstract: 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… ▽ More 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. △ Less

Submitted 24 May, 2007; originally announced May 2007.

Comments: 5 pages, 4 figures, submitted to ApJL

Journal ref: Astrophys.J.665:L1-L4,2007

A Scaling Index Analysis of the WMAP three year data: Signatures of non-Gaussianities and Asymmetries in the CMB

Authors: C. Raeth, P. Schuecker, A. J. Banday

Abstract: Local scaling properties of the co-added foreground-cleaned three-year Wilkinson Microwave Anisotropy Probe (WMAP) data are estimated using weighted scaling indices. The scaling index method (SIM) is - for the first time - adapted and applied to the case of spherical symmetric spatial data. The results are compared with 1000 Monte Carlo simulations based on Gaussian fluctuations with a best fit… ▽ More Local scaling properties of the co-added foreground-cleaned three-year Wilkinson Microwave Anisotropy Probe (WMAP) data are estimated using weighted scaling indices. The scaling index method (SIM) is - for the first time - adapted and applied to the case of spherical symmetric spatial data. The results are compared with 1000 Monte Carlo simulations based on Gaussian fluctuations with a best fit $Λ$CDM power spectrum and WMAP-like beam and noise properties. Statistical quantities based on the scaling indices, namely the moments of the distribution and probability-based measures are determined. We find for most of the test statistics significant deviations from the Gaussian hypothesis. We find pronounced asymmetries, which can be interpreted as a global lack of structure in the northern hemisphere, which is consistent with previous findings. Furthermore, we detect a localized anomaly in the southern hemisphere, which gives rise to highly significant signature for non-Gaussianity in the spectrum of scaling indices. We identify this signature as the cold spot, which was also already detected in the first year WMAP data. Our results provide further evidence for both the presence of non-Gaussianities and asymmetries in the WMAP three-year data. More detailed bandand year-wise analyses are needed to elucidate the origin of the detected anomalies. In either case the scaling indices provide powerful nonlinear statistics to analyse CMB maps. △ Less

Submitted 6 February, 2007; originally announced February 2007.

Comments: submitted to MNRAS

Journal ref: MNRAS, 380, 466 (2007)

Hemispherical power asymmetry in the three-year Wilkinson Microwave Anisotropy Probe sky maps

Authors: H. K. Eriksen, A. J. Banday, K. M. Gorski, F. K. Hansen, P. B. Lilje

Abstract: 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… ▽ More 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. △ Less

Submitted 29 May, 2007; v1 submitted 4 January, 2007; originally announced January 2007.

Comments: 4 pages, 2 figures, accepted for publication in ApJL

Journal ref: Astrophys.J.660:L81-L84,2007

CMB signal in WMAP 3yr data with FastICA

Authors: D. Maino, S. Donzelli, A. J. Banday, F. Stivoli, C. Baccigalupi

Abstract: We present an application of the fast Independent Component Analysis (FastICA) to the WMAP 3yr data with the goal of extracting the CMB signal. We evaluate the confidence of our results by means of Monte Carlo simulations including CMB, foreground contaminations and instrumental noise specific of each WMAP frequency band. We perform a complete analysis involving all or a subset of the WMAP chann… ▽ More We present an application of the fast Independent Component Analysis (FastICA) to the WMAP 3yr data with the goal of extracting the CMB signal. We evaluate the confidence of our results by means of Monte Carlo simulations including CMB, foreground contaminations and instrumental noise specific of each WMAP frequency band. We perform a complete analysis involving all or a subset of the WMAP channels in order to select the optimal combination for CMB extraction, using the frequency scaling of the reconstructed component as a figure of merit. We found that the combination KQVW provides the best CMB frequency scaling, indicating that the low frequency foreground contamination in Q, V and W bands is better traced by the emission in the K band. The CMB angular power spectrum is recovered up to the degree scale, it is consistent within errors for all WMAP channel combination considered, and in close agreement with the WMAP 3yr results. We perform a statistical analysis of the recovered CMB pattern, and confirm the sky asymmetry reported in several previous works with independent techniques. △ Less

Submitted 8 September, 2006; originally announced September 2006.

Comments: 10 pages, 7 figures, submitted to MNRAS

Journal ref: MNRAS 374 (2007) 1207

A re-analysis of the three-year WMAP temperature power spectrum and likelihood

Authors: H. K. Eriksen, Greg Huey, R. Saha, F. K. Hansen, J. Dick, A. J. Banday, K. M. Gorski, P. Jain, J. B. Jewell, L. Knox, D. L. Larson, I. J. O'Dwyer, T. Souradeep, B. D. Wandelt

Abstract: 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… ▽ More 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. △ Less

Submitted 2 October, 2006; v1 submitted 5 June, 2006; originally announced June 2006.

Comments: 14 pages, 7 figures. Accepted for publication in ApJ. Numerical results unchanged, but interpretation sharpened: Likelihood approximation issues at l=13-30 far more important than potential foreground issues at l <= 12. Gibbs products (spectrum and sky samples, and "easy-to-use" likelihood module) available from http://www.astro.uio.no/~hke/ under "Research"

Journal ref: Astrophys.J.656:641-652,2007

Bianchi Type VII_h Models and the WMAP 3-year Data

Authors: T. R. Jaffe, A. J. Banday, H. K. Eriksen, K. M. Gorski, F. K. Hansen

Abstract: 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 t… ▽ More 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. △ Less

Submitted 2 June, 2006; originally announced June 2006.

Comments: 5 pages, 4 figures, submitted as A&A Research Note

Journal ref: Astron.Astrophys.460:393,2006

Bayesian foreground analysis with CMB data

Authors: H. K. Eriksen, C. Dickinson, C. R. Lawrence, C. Baccigalupi, A. J. Banday, K. M. Gorski, F. K. Hansen, E. Pierpaoli, M. D. Seiffert

Abstract: 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 er… ▽ More 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. △ Less

Submitted 7 April, 2006; originally announced April 2006.

Comments: 8 pages; for proceedings from CMB workshop at Irvine, March 2006

Journal ref: New Astron.Rev. 50 (2006) 861-867

Fast and Efficient Template Fitting of Deterministic Anisotropic Cosmological Models Applied to WMAP Data

Authors: T. R. Jaffe, A. J. Banday, H. K. Eriksen, K. M. Gorski, F. K. Hansen

Abstract: 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… ▽ More 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. △ Less

Submitted 31 March, 2006; originally announced March 2006.

Comments: 14 pages, 11 figure, to be published in ApJ

Journal ref: Astrophys.J. 643 (2006) 616-629

Foreground Subtraction of Cosmic Microwave Background Maps using WI-FIT (Wavelet based hIgh resolution Fitting of Internal Templates)

Authors: F. K. Hansen, A. J. Banday, H. K. Eriksen, K. M. Gorski, P. B. Lilje

Abstract: 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 onl… ▽ More 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. △ Less

Submitted 13 March, 2006; originally announced March 2006.

Comments: Submitted to ApJ

Journal ref: Astrophys.J.648:784-796,2006

No Higher Criticism of the Bianchi Corrected WMAP Data

Authors: L. Cayon, A. J. Banday, T. Jaffe, H. K. Eriksen, F. K. Hansen, K. M. Gorski, J. Jin

Abstract: 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… ▽ More 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. △ Less

Submitted 1 February, 2006; originally announced February 2006.

Comments: submitted to MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.369:598,2006

On the Viability of Bianchi Type VIIh Models with Dark Energy

Authors: T. R. Jaffe, S. Hervik, A. J. Banday, K. M. Gorski

Abstract: 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… ▽ More 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. △ Less

Submitted 16 December, 2005; originally announced December 2005.

Comments: submitted to ApJ

Journal ref: Astrophys.J.644:701-708,2006

Lyman-alpha forest-CMB cross-correlation and the search for the ionized baryons at high redshift

Authors: Rupert A. C. Croft, A. J. Banday, Lars Hernquist

Abstract: The intergalactic neutral hydrogen which is responsible for the Lyman alpha forest of quasar absorption is a tracer of much larger amounts of ionised hydrogen. The ionised component has yet to be detected directly, but is expected to scatter CMB photons via the Sunyaev-Zel'dovich (SZ) effect. We use hydrodynamic simulations of a LambdaCDM universe to create mock quasar spectra and CMB sky maps.… ▽ More The intergalactic neutral hydrogen which is responsible for the Lyman alpha forest of quasar absorption is a tracer of much larger amounts of ionised hydrogen. The ionised component has yet to be detected directly, but is expected to scatter CMB photons via the Sunyaev-Zel'dovich (SZ) effect. We use hydrodynamic simulations of a LambdaCDM universe to create mock quasar spectra and CMB sky maps. We find that the high-z Lya forest gas causes temperature fluctuations of the order of 1 muK rms in the CMB on arcmin scales. The kinetic and thermal SZ effects have a similar magnitude at z=3, with the thermal effect becoming relatively weaker as expected at higher z. The CMB signal associated with lines of sight having HI column densities > 10^18 cm^-2 is only marginally stronger than that for lower column densities. The strong dependence of rms temperature fluctuation on mean Lya absorbed flux, however, suggests that the CMB signal effectively arises in lower density material. We investigate the use of the cross-correlation of the Lya forest and the microwave background to detect the SZ effect at redshifts 2-4. In so doing we are able to set direct limits on the density of diffuse ionised intergalactic baryons. We carry out a preliminary comparison at a mean redshift z=3 of 3488 quasar spectra from SDSS Data Release 3 and the WMAP first year data. Assuming that the baryons are clustered as in a LambdaCDM cosmology, and have the same mean temperature, the cross-correlation yields a weak limit on the cosmic density of ionised baryons Omega_(b,I), which is Omega_(b,I) < 0.8 at 95% confidence. With data from upcoming CMB telescopes, we anticipate that a direct detection of the high redshift ionised IGM will soon be possible, providing an important consistency check on cosmological models. △ Less

Submitted 14 December, 2005; originally announced December 2005.

Comments: 14 pages, 10 figures, submitted to MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.369:1090-1102,2006

A determination of the Spectra of Galactic components observed by WMAP

Authors: R. D. Davies, C. Dickinson, A. J. Banday, T. R. Jaffe, K. M. Gorski, R. J. Davis

Abstract: 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 emissio… ▽ More 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). △ Less

Submitted 7 April, 2006; v1 submitted 14 November, 2005; originally announced November 2005.

Comments: 18 pages, 6 figures, revised version uses cross-correlation method rather than T-T method. Paper re-organised and sent back to referee

Journal ref: Mon.Not.Roy.Astron.Soc.370:1125-1139,2006

CMB component separation by parameter estimation

Authors: H. K. Eriksen, C. Dickinson, C. R. Lawrence, C. Baccigalupi, A. J. Banday, K. M. Gorski, F. K. Hansen, P. B. Lilje, E. Pierpaoli, M. D. Seiffert, K. M. Smith, K. Vanderlinde

Abstract: 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… ▽ More 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. △ Less

Submitted 11 August, 2005; originally announced August 2005.

Comments: 20 pages, 10 figures, submitted to ApJ. For a high-resolution version, see http://www.astro.uio.no/~hke/docs/eriksen_et_al_fgfit.ps

Journal ref: Astrophys.J.641:665-682,2006

Astro-ph communication: Simulations of the WMAP Internal Linear Combination sky map

Authors: H. K. Eriksen, A. J. Banday, K. M. Gorski, P. B. Lilje

Abstract: 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, Fig… ▽ More 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. △ Less

Submitted 8 August, 2005; originally announced August 2005.

Comments: 2 pages, submitted to astro-ph only

Multipole vector anomalies in the first-year WMAP data: a cut-sky analysis

Authors: P. Bielewicz, H. K. Eriksen, A. J. Banday, K. M. Gorski, P. B. Lilje

Abstract: 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… ▽ More 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. △ Less

Submitted 15 August, 2005; v1 submitted 8 July, 2005; originally announced July 2005.

Comments: 13 pages, 2 figures; corrected typos; added reference. Accepted for publication in ApJ

Journal ref: Astrophys.J.635:750-760,2005

Evidence of vorticity and shear at large angular scales in the WMAP data: a violation of cosmological isotropy?

Authors: T. R. Jaffe, A. J. Banday, H. K. Eriksen, K. M. Gorski, F. K. Hansen

Abstract: 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… ▽ More 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. △ Less

Submitted 2 September, 2005; v1 submitted 9 March, 2005; originally announced March 2005.

Comments: 4 pages, 3 figures; emulateapj.cls; ApJL accepted version plus fixed error in vorticity calculation (sqrt(2) off in Table 1, abstract, and conclusions); basic conclusions unchanged

Journal ref: Astrophys.J.629:L1-L4,2005

HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere

Authors: K. M. Gorski, E. Hivon, A. J. Banday, B. D. Wandelt, F. K. Hansen, M. Reinecke, M. Bartelman

Abstract: 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 proce… ▽ More 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. △ Less

Submitted 21 September, 2004; originally announced September 2004.

Journal ref: Astrophys.J.622:759-771,2005

The N-point correlation functions of the first-year Wilkinson Microwave Anisotropy Probe sky maps

Authors: H. K. Eriksen, A. J. Banday, K. M. Gorski, P. B. Lilje

Abstract: 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 18… ▽ More 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. △ Less

Submitted 6 December, 2004; v1 submitted 13 July, 2004; originally announced July 2004.

Comments: 14 pages, 9 figures, version to appear in ApJ

Journal ref: Astrophys.J. 622 (2005) 58-71

Power spectrum estimation from high-resolution maps by Gibbs sampling

Authors: H. K. Eriksen, I. J. O'Dwyer, J. B. Jewell, B. D. Wandelt, D. L. Larson, K. M. Gorski, S. Levin, A. J. Banday, P. B. Lilje

Abstract: 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 tha… ▽ More 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. △ Less

Submitted 23 August, 2004; v1 submitted 1 July, 2004; originally announced July 2004.

Comments: 15 pages, 11 figures, accepted for publication in ApJS

Journal ref: Astrophys.J.Suppl. 155 (2004) 227-241

Bayesian Power Spectrum Analysis of the First-Year WMAP data

Authors: I. J. O'Dwyer, H. K. Eriksen, B. D. Wandelt, J. B. Jewell, D. L. Larson, K. M. Gorski, A. J. Banday, S. Levin, P. B. Lilje

Abstract: 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 anisotrop… ▽ More 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. △ Less

Submitted 7 December, 2004; v1 submitted 1 July, 2004; originally announced July 2004.

Comments: 5 pages. 4 figures. For additional information and data see http://www.astro.uiuc.edu/~iodwyer/research#wmap

Journal ref: Astrophys.J. 617 (2004) L99-102

Cosmological parameters and the WMAP data revisited

Authors: F. K. Hansen, A. Balbi, A. J. Banday, K. M. Gorski

Abstract: 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 ref… ▽ More 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. △ Less

Submitted 9 June, 2004; originally announced June 2004.

Comments: submitted to MNRAS