Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCMB power spectrum estimation using wavelets
135
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Observations of the Cosmic Microwave Background (CMB) provide increasingly accurate information about the structure of the Universe at the recombination epoch. Most of this information is encoded in the angular power spectrum of the CMB. The aim of this work is to propose a versatile and powerful method for spectral estimation on the sphere which can easily deal with non-stationarity, foregrounds and multiple experiments with various specifications. In this paper, we use needlets (wavelets) on the sphere to construct natural and efficient spectral estimators for partially observed and beamed CMB with non stationary noise. In the case of a single experiment, we compare this method with Pseudo-$C_ell$ methods. The performance of the needlet spectral estimators (NSE) compares very favorably to the best Pseudo--$C_ell$ estimators, over the whole multipole range. On simulations with a simple model (CMB + uncorrelated noise with known variance per pixel + mask), they perform uniformly better. Their distinctive ability to aggregate many different experiments, to control the propagation of errors and to produce a single wide-band error bars is highlighted. The needlet spectral estimator is a powerful, tunable tool which is very well suited to angular power spectrum estimation of spherical data such as incomplete and noisy CMB maps.
rate research
Read More
The cosmic microwave background (CMB) power spectrum is a powerful cosmological probe as it entails almost all the statistical information of the CMB perturbations. Having access to only one sky, the CMB power spectrum measured by our experiments is only a realization of the true underlying angular power spectrum. In this paper we aim to recover the true underlying CMB power spectrum from the one realization that we have without a need to know the cosmological parameters. The sparsity of the CMB power spectrum is first investigated in two dictionaries; Discrete Cosine Transform (DCT) and Wavelet Transform (WT). The CMB power spectrum can be recovered with only a few percentage of the coefficients in both of these dictionaries and hence is very compressible in these dictionaries. We study the performance of these dictionaries in smoothing a set of simulated power spectra. Based on this, we develop a technique that estimates the true underlying CMB power spectrum from data, i.e. without a need to know the cosmological parameters. This smooth estimated spectrum can be used to simulate CMB maps with similar properties to the true CMB simulations with the correct cosmological parameters. This allows us to make Monte Carlo simulations in a given project, without having to know the cosmological parameters. The developed IDL code, TOUSI, for Theoretical pOwer spectrUm using Sparse estImation, will be released with the next version of ISAP.
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.
We revisit the problem of exact CMB likelihood and power spectrum estimation with the goal of minimizing computational cost through linear compression. This idea was originally proposed for CMB purposes by Tegmark et al. (1997), and here we develop it into a fully working computational framework for large-scale polarization analysis, adopting WMAP as a worked example. We compare five different linear bases (pixel space, harmonic space, noise covariance eigenvectors, signal-to-noise covariance eigenvectors and signal-plus-noise covariance eigenvectors) in terms of compression efficiency, and find that the computationally most efficient basis is the signal-to-noise eigenvector basis, which is closely related to the Karhunen-Loeve and Principal Component transforms, in agreement with previous suggestions. For this basis, the information in 6836 unmasked WMAP sky map pixels can be compressed into a smaller set of 3102 modes, with a maximum error increase of any single multipole of 3.8% at $ellle32$, and a maximum shift in the mean values of a joint distribution of an amplitude--tilt model of 0.006$sigma$. This compression reduces the computational cost of a single likelihood evaluation by a factor of 5, from 38 to 7.5 CPU seconds, and it also results in a more robust likelihood by implicitly regularizing nearly degenerate modes. Finally, we use the same compression framework to formulate a numerically stable and computationally efficient variation of the Quadratic Maximum Likelihood implementation that requires less than 3 GB of memory and 2 CPU minutes per iteration for $ell le 32$, rendering low-$ell$ QML CMB power spectrum analysis fully tractable on a standard laptop.
We present two novel methods for the estimation of the angular power spectrum of cosmic microwave background (CMB) anisotropies. We assume an absolute CMB experiment with arbitrary asymmetric beams and arbitrary sky coverage. The methods differ from earlier ones in that the power spectrum is estimated directly from time-ordered data, without first compressing the data into a sky map, and they take into account the effect of asymmetric beams. In particular, they correct the beam-induced leakage from temperature to polarization. The methods are applicable to a case where part of the sky has been masked out to remove foreground contamination, leaving a pure CMB signal, but incomplete sky coverage. The first method (DQML) is derived as the optimal quadratic estimator, which simultaneously yields an unbiased spectrum estimate and minimizes its variance. We successfully apply it to multipoles up to $ell$=200. The second method is derived as a weak-signal approximation from the first one. It yields an unbiased estimate for the full multipole range, but relaxes the requirement of minimal variance. We validate the methods with simulations for the 70 GHz channel of {tt Planck} surveyor, and demonstrate that we are able to correct the beam effects in the $TT$, $EE$, $BB$, and $TE$ spectra up to multipole $ell$=1500. Together the two methods cover the complete multipole range with no gap in between.
We describe a class of observational strategies for probing the anisotropies in the cosmic microwave background (CMB) where the instrument scans on rings which can be combined into an n-torus, the {em ring torus}. This class has the remarkable property that it allows exact maximum likelihood power spectrum estimation in of order $N^2$ operations (if the size of the data set is $N$) under circumstances which would previously have made this analysis intractable: correlated receiver noise, arbitrary asymmetric beam shapes and far side lobes, non-uniform distribution of integration time on the sky and partial sky coverage. This ease of computation gives us an important theoretical tool for understanding the impact of instrumental effects on CMB observables and hence for the design and analysis of the CMB observations of the future. There are members of this class which closely approximate the MAP and Planck satellite missions. We present a numerical example where we apply our ring torus methods to a simulated data set from a CMB mission covering a 20 degree patch on the sky to compute the maximum likelihood estimate of the power spectrum $C_ell$ with unprecedented efficiency.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
