Linearized iterative least-squares (LIL): A parameter fitting algorithm for component separation in multifrequency CMB experiments such as Planck

Abstract

We present an efficient algorithm for the least squares parameter fitting optimized for component separation in multi-frequency CMB experiments. We sidestep some of the problems associated with non-linear optimization by taking advantage of the quasi-linear nature of the foreground model. We demonstrate our algorithm, linearized iterative least-squares (LIL), on the publicly available Planck sky model FFP6 simulations and compare our result with the other algorithms. We work at full Planck resolution and show that degrading the resolution of all channels to that of the lowest frequency channel is not necessary. Finally we present results for the publicly available Planck data. Our algorithm is extremely fast, fitting 6 parameters to 7 lowest Planck channels at full resolution (50 million pixels) in less than 160 CPU-minutes (or few minutes running in parallel on few tens of cores). LIL is therefore easily scalable to future experiments which may have even higher resolution and more frequency channels. We also naturally propagate the uncertainties in different parameters due to noise in the maps as well as degeneracies between the parameters to the final errors on the parameters using Fisher matrix. One indirect application of LIL could be a front-end for Bayesian parameter fitting to find the maximum of the likelihood to be used as the starting point for the Gibbs sampling. We show for rare components, such as the carbon-monoxide emission, present in small fraction of sky, the optimal approach should combine parameter fitting with model selection. LIL may also be useful in other astrophysical applications which satisfy the quasi-linearity criteria.

Accepted version

References in corpus (40)

• Maps of Dust IR Emission for Use in Estimation of Reddening and CMBR Foregrounds
• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Final Maps and Results
• The Milky Way in Molecular Clouds: A New Complete CO Survey
• Planck 2013 results. I. Overview of products and scientific results
• The Cosmic Microwave Background Spectrum from the Full COBE/FIRAS Data Set
• BICEP2 I: Detection Of B-mode Polarization at Degree Angular Scales
• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Foreground Emission
• MASTER of the CMB Anisotropy Power Spectrum: A Fast Method for Statistical Analysis of Large and Complex CMB Data Sets
• The COBE Diffuse Infrared Background Experiment Search for the Cosmic Infrared Background: I. Limits and Detections
• The Primordial Inflation Explorer (PIXIE): A Nulling Polarimeter for Cosmic Microwave Background Observations
• Extrapolation of Galactic Dust Emission at 100 Microns to CMBR Frequencies Using FIRAS
• Electric Dipole Radiation from Spinning Dust Grains
• Diffuse Galactic Emission from Spinning Dust Grains
• Joint Bayesian component separation and CMB power spectrum estimation
• A full sky, low foreground, high resolution CMB map from WMAP
• Planck 2013 results. XII. Component separation
• Component separation methods for the Planck mission
• The pre-launch Planck Sky Model: a model of sky emission at submillimetre to centimetre wavelengths
• PRISM (Polarized Radiation Imaging and Spectroscopy Mission): An Extended White Paper
• Magnetic Dipole Microwave Emission from Dust Grains
• Planck 2013 results. IX. HFI spectral response
• Foreground separation methods for satellite observations of the cosmic microwave background
• MILCA, a Modified Internal Linear Combination Algorithm to extract astrophysical emissions from multi-frequency sky maps
• All-sky astrophysical component separation with Fast Independent Component Analysis (FastICA)
• CMB component separation by parameter estimation
• Planck 2013 results. VII. HFI time response and beams
• Planck 2013 results. XXI. All-sky Compton parameter power spectrum and high-order statistics
• Multi-resolution internal template cleaning: An application to the Wilkinson Microwave Anisotropy Probe 7-yr polarization data
• The contribution of the Unresolved Extragalactic Radio Sources to the Brightness Temperature of the sky
• Planck 2013 results. V. LFI calibration
• All-sky component separation for the Planck mission
• Planck 2013 results. IV. Low Frequency Instrument beams and window functions
• WMAP 3yr data with the CCA: anomalous emission and impact of component separation on the CMB power spectrum
• Foreground Subtraction of Cosmic Microwave Background Maps using WI-FIT (Wavelet based hIgh resolution Fitting of Internal Templates)
• H.E.S.S. Observations of the Crab during its March 2013 GeV Gamma-Ray Flare
• CMB power spectrum estimation and map reconstruction with the Expectation - Maximization algorithm
• Separation of correlated astrophysical sources using multiple-lag data covariance matrices
• Component separation with flexible models. Application to the separation of astrophysical emissions
• Constraints on $\mu$-distortion fluctuations and primordial non-Gaussianity from Planck data

Cited by in corpus (10)

• Thermal SZ fluctuations in the ICM: probing turbulence and thermodynamics in Coma cluster with $Planck$
• CMB spectral distortions revisited: a new take on $\mu$ distortions and primordial non-Gaussianities from FIRAS data
• Probing Cosmology and Cluster Astrophysics with Multi-Wavelength Surveys I. Correlation Statistics
• An alternative validation strategy for the Planck cluster catalog and $y$-distortion maps
• Limits on the fluctuating part of $y$-type distortion monopole from Planck and SPT results
• Constraints on non-resonant photon-axion conversion from the Planck satellite data
• Improved Galactic Foreground Removal for B-Modes Detection with Clustering Methods
• Spectral imaging of galaxy clusters with Planck
• Detection of WHIM in the Planck data using Stack First approach
• Data driven foreground clustering approach to component separation in multifrequency CMB experiments: A new Planck CMB map