اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدCMB likelihood approximation by a Gaussianized Blackwell-Rao estimator
62
0
0.0
(
0
)
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
We introduce a new CMB temperature likelihood approximation called the Gaussianized Blackwell-Rao (GBR) estimator. This estimator is derived by transforming the observed marginal power spectrum distributions obtained by the CMB Gibbs sampler into standard univariate Gaussians, and then approximate their joint transformed distribution by a multivariate Gaussian. The method is exact for full-sky coverage and uniform noise, and an excellent approximation for sky cuts and scanning patterns relevant for modern satellite experiments such as WMAP and Planck. A single evaluation of this estimator between l=2 and 200 takes ~0.2 CPU milliseconds, while for comparison, a single pixel space likelihood evaluation between l=2 and 30 for a map with ~2500 pixels requires ~20 seconds. We apply this tool to the 5-year WMAP temperature data, and re-estimate the angular temperature power spectrum, $C_{ell}$, and likelihood, L(C_l), for l<=200, and derive new cosmological parameters for the standard six-parameter LambdaCDM model. Our spectrum is in excellent agreement with the official WMAP spectrum, but we find slight differences in the derived cosmological parameters. Most importantly, the spectral index of scalar perturbations is n_s=0.973 +/- 0.014, 1.9 sigma away from unity and 0.6 sigma higher than the official WMAP result, n_s = 0.965 +/- 0.014. This suggests that an exact likelihood treatment is required to higher ls than previously believed, reinforcing and extending our conclusions from the 3-year WMAP analysis. In that case, we found that the sub-optimal likelihood approximation adopted between l=12 and 30 by the WMAP team biased n_s low by 0.4 sigma, while here we find that the same approximation between l=30 and 200 introduces a bias of 0.6 sigma in n_s.
قيم البحث
اقرأ أيضاً
We propose an efficient Bayesian MCMC algorithm for estimating cosmological parameters from CMB data without use of likelihood approximations. It builds on a previously developed Gibbs sampling framework that allows for exploration of the joint CMB s
ky signal and power spectrum posterior, P(s,Cl|d), and addresses a long-standing problem of efficient parameter estimation simultaneously in high and low signal-to-noise regimes. To achieve this, our new algorithm introduces a joint Markov Chain move in which both the signal map and power spectrum are synchronously modified, by rescaling the map according to the proposed power spectrum before evaluating the Metropolis-Hastings accept probability. Such a move was already introduced by Jewell et al. (2009), who used it to explore low signal-to-noise posteriors. However, they also found that the same algorithm is inefficient in the high signal-to-noise regime, since a brute-force rescaling operation does not account for phase information. This problem is mitigated in the new algorithm by subtracting the Wiener filter mean field from the proposed map prior to rescaling, leaving high signal-to-noise information invariant in the joint step, and effectively only rescaling the low signal-to-noise component. To explore the full posterior, the new joint move is then interleaved with a standard conditional Gibbs sky map move. We apply our new algorithm to simplified simulations for which we can evaluate the exact posterior to study both its accuracy and performance, and find good agreement with the exact posterior; marginal means agree to less than 0.006 sigma, and standard deviations to better than 3%. The Markov Chain correlation length is of the same order of magnitude as those obtained by other standard samplers in the field.
We present ECLIPSE (Efficient Cmb poLarization and Intensity Power Spectra Estimator), an optimized implementation of the Quadratic Maximum Likelihood (QML) method for the estimation of the power spectra of the Cosmic Microwave Background (CMB). This
approach allows one to reduce significantly the computational costs associated to this technique, allowing to estimate the power spectra up to higher multipoles than previous implementations. In particular, for a resolution of $N_mathrm{side}=64$, $ell_{mathrm{max}}=192$ and a typical Galactic mask, the number of operations can be reduced by approximately a factor of 1000 in a full analysis including intensity and polarization with respect to an efficient direct implementation of the method. In addition, if one is interested in studying only polarization, it is possible to obtain the power spectra of the E and B modes with a further reduction of computational resources without degrading the results. We also show that for experiments observing a small fraction of the sky, the Fisher matrix becomes singular and, in this case, the standard QML can not be applied. To solve this problem, we have developed a binned version of the method that is unbiased and of minimum variance. We also test the robustness of the QML estimator when the assumed fiducial model differs from that of the sky and show the performance of an iterative approach. Finally, we present a comparison of the results obtained by QML and a pseudo-$C_{ell}$ estimator (NaMaster) for a next-generation satellite, showing that, as expected, QML produces significantly smaller errors at low multipoles. The ECLIPSE fast QML code developed in this work will be made publicly available.
The maximum likelihood estimator plays a fundamental role in statistics. However, for many models, the estimators do not have closed-form expressions. This limitation can be significant in situations where estimates and predictions need to be compute
d in real-time, such as in applications based on embedded technology, in which numerical methods can not be implemented. This paper provides a modification in the maximum likelihood estimator that allows us to obtain the estimators in closed-form expressions under some conditions. Under mild conditions, the estimator is invariant under one-to-one transformations, consistent, and has an asymptotic normal distribution. The proposed modified version of the maximum likelihood estimator is illustrated on the Gamma, Nakagami, and Beta distributions and compared with the standard maximum likelihood estimator.
We develop a Maximum Likelihood estimator (MLE) to measure the masses of galaxy clusters through the impact of gravitational lensing on the temperature and polarization anisotropies of the cosmic microwave background (CMB). We show that, at low noise
levels in temperature, this optimal estimator outperforms the standard quadratic estimator by a factor of two. For polarization, we show that the Stokes Q/U maps can be used instead of the traditional E- and B-mode maps without losing information. We test and quantify the bias in the recovered lensing mass for a comprehensive list of potential systematic errors. Using realistic simulations, we examine the cluster mass uncertainties from CMB-cluster lensing as a function of an experiments beam size and noise level. We predict the cluster mass uncertainties will be 3 - 6% for SPT-3G, AdvACT, and Simons Array experiments with 10,000 clusters and less than 1% for the CMB-S4 experiment with a sample containing 100,000 clusters. The mass constraints from CMB polarization are very sensitive to the experimental beam size and map noise level: for a factor of three reduction in either the beam size or noise level, the lensing signal-to-noise improves by roughly a factor of two.
A great deal of experimental effort is currently being devoted to the precise measurements of the cosmic microwave background (CMB) sky in temperature and polarisation. Satellites, balloon-borne, and ground-based experiments scrutinize the CMB sky at
multiple scales, and therefore enable to investigate not only the evolution of the early Universe, but also its late-time physics with unprecedented accuracy. The pipeline leading from time ordered data as collected by the instrument to the final product is highly structured. Moreover, it has also to provide accurate estimates of statistical and systematic uncertainties connected to the specific experiment. In this paper, we review likelihood approaches targeted to the analysis of the CMB signal at different scales, and to the estimation of key cosmological parameters. We consider methods that analyze the data in the spatial (i.e., pixel-based) or harmonic domain. We highlight the most relevant aspects of each approach and compare their performance.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
