Multi-component power spectra estimation method for multi-detector observations of the Cosmic Microwave Background

We present a new method for multi-component power spectra estimation in multi-frequency observations of the CMB. Our method is based on matching a model to the cross and auto power spectra of observed maps. All the component power spectra are estimated, as well as their mixing matrix. Noise power spectra are also estimated. The method has been applied to full-sky Planck simulations containing five astrophysical components and white noise. The beam smoothing effect is taken into account.Comment: 7 pages, 1 figure, proc. of the CMBnet workshop, 20-21 Feb. 2003 Oxford, UK. New Astronomy Reviews (eds. A. Melchiorri, J.I. Silk) in pres

Bayesian blind component separation for Cosmic Microwave Background observations

We present a technique for the blind separation of components in CMB data. The method uses a spectral EM algorithm which recovers simultaneously component templates, their emission law as a function of wavelength, and noise levels. We test the method on Planck HFI simulated observations featuring 3 astrophysical components.Comment: 15 pages, 5 figures, to appear in the Proceedings of the MAXENT 2001 international worksho

Independent Component analysis of the Cosmic Microwave Background

This paper presents an application of ICA to astronomical imaging. A first section describes the astrophysical context and motivates the use of source separation ideas. A second section describes our approach to the problem: the use of a noisy Gaussian stationary model. This technique uses spectral diversity and take explicitly into account contamination by additive noise. Preliminary and extremely encouraging results on realistic synthetic signals and on real data will be presented at the conferenc

A Method for Individual Source Brightness Estimation in Single- and Multi-band Data

We present a method of reliably extracting the flux of individual sources from sky maps in the presence of noise and a source population in which number counts are a steeply falling function of flux. The method is an extension of a standard Bayesian procedure in the millimeter/submillimeter literature. As in the standard method, the prior applied to source flux measurements is derived from an estimate of the source counts as a function of flux, dN/dS. The key feature of the new method is that it enables reliable extraction of properties of individual sources, which previous methods in the literature do not. We first present the method for extracting individual source fluxes from data in a single observing band, then we extend the method to multiple bands, including prior information about the spectral behavior of the source population(s). The multi-band estimation technique is particularly relevant for classifying individual sources into populations according to their spectral behavior. We find that proper treatment of the correlated prior information between observing bands is key to avoiding significant biases in estimations of multi-band fluxes and spectral behavior, biases which lead to significant numbers of misclassified sources. We test the single- and multi-band versions of the method using simulated observations with observing parameters similar to that of the South Pole Telescope data used in Vieira, et al. (2010).Comment: 11 emulateapj pages, 3 figures, revised to match published versio

Multi-Detector Multi-Component spectral matching and applications for CMB data analysis

We present a new method for analyzing multi--detector maps containing contributions from several components. Our method, based on matching the data to a model in the spectral domain, permits to estimate jointly the spatial power spectra of the components and of the noise, as well as the mixing coefficients. It is of particular relevance for the analysis of millimeter--wave maps containing a contribution from CMB anisotropies.Comment: 15 pages, 7 Postscript figures, submitted to MNRA

BLAST Observations of the South Ecliptic Pole field: Number Counts and Source Catalogs

We present results from a survey carried out by the Balloon-borne Large Aperture Submillimeter Telescope (BLAST) on a 9 deg^2 field near the South Ecliptic Pole at 250, 350 and 500 {\mu}m. The median 1{\sigma} depths of the maps are 36.0, 26.4 and 18.4 mJy, respectively. We apply a statistical method to estimate submillimeter galaxy number counts and find that they are in agreement with other measurements made with the same instrument and with the more recent results from Herschel/SPIRE. Thanks to the large field observed, the new measurements give additional constraints on the bright end of the counts. We identify 132, 89 and 61 sources with S/N>4 at 250, 350, 500 {\mu}m, respectively and provide a multi-wavelength combined catalog of 232 sources with a significance >4{\sigma} in at least one BLAST band. The new BLAST maps and catalogs are available publicly at http://blastexperiment.info.Comment: 25 pages, 6 figures, 4 tables, Accepted by ApJS. Maps and catalogs available at http://blastexperiment.info

BLAST: the Redshift Survey

The Balloon-borne Large Aperture Submillimeter Telescope (BLAST) has recently surveyed ~8.7 deg^2 centered on GOODS-South at 250, 350, and 500 microns. In Dye et al. (2009) we presented the catalogue of sources detected at 5-sigma in at least one band in this field and the probable counterparts to these sources in other wavebands. In this paper, we present the results of a redshift survey in which we succeeded in measuring redshifts for 82 of these counterparts. The spectra show that the BLAST counterparts are mostly star-forming galaxies but not extreme ones when compared to those found in the Sloan Digital Sky Survey. Roughly one quarter of the BLAST counterparts contain an active nucleus. We have used the spectroscopic redshifts to carry out a test of the ability of photometric redshift methods to estimate the redshifts of dusty galaxies, showing that the standard methods work well even when a galaxy contains a large amount of dust. We have also investigated the cases where there are two possible counterparts to the BLAST source, finding that in at least half of these there is evidence that the two galaxies are physically associated, either because they are interacting or because they are in the same large-scale structure. Finally, we have made the first direct measurements of the luminosity function in the three BLAST bands. We find strong evolution out to z=1, in the sense that there is a large increase in the space-density of the most luminous galaxies. We have also investigated the evolution of the dust-mass function, finding similar strong evolution in the space-density of the galaxies with the largest dust masses, showing that the luminosity evolution seen in many wavebands is associated with an increase in the reservoir of interstellar matter in galaxies.Comment: Accepted for publication in the Astrophysical Journal. Maps and associated results are available at http://blastexperiment.info

Impact of calibration errors on CMB component separation using FastICA and ILC

The separation of emissions from different astrophysical processes is an important step towards the understanding of observational data. This topic of component separation is of particular importance in the observation of the relic Cosmic Microwave Background Radiation, as performed by the WMAP satellite and the more recent Planck mission, launched May 14th, 2009 from Kourou and currently taking data. When performing any sort of component separation, some assumptions about the components must be used. One assumption that many techniques typically use is knowledge of the frequency scaling of one or more components. This assumption may be broken in the presence of calibration errors. Here we compare, in the context of imperfect calibration, the recovery of a clean map of emission of the Cosmic Microwave Background from observational data with two methods: FastICA (which makes no assumption of the frequency scaling of the components), and an `Internal Linear Combination' (ILC), which explicitly extracts a component with a given frequency scaling. We find that even in the presence of small calibration errors with a Planck-style mission, the ILC method can lead to inaccurate CMB reconstruction in the high signal-to-noise regime, because of partial cancellation of the CMB emission in the recovered map. While there is no indication that the failure of the ILC will translate to other foreground cleaning or component separation techniques, we propose that all methods which assume knowledge of the frequency scaling of one or more components be careful to estimate the effects of calibration errors.Comment: 13 pages, 5 figure

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

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.Comment: 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.pd