Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMeasuring galaxy cluster masses with CMB lensing using a Maximum Likelihood estimator: Statistical and systematic error budgets for future experiments
74
0
0.0
(
0
)
Added by
Srinivasan Raghunathan
Publication date
2017
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
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.
We present and describe im3shape, a new publicly available galaxy shape measurement code for weak gravitational lensing shear. im3shape performs a maximum likelihood fit of a bulge-plus-disc galaxy model to noisy images, incorporating an applied point spread function. We detail challenges faced and choices made in its design and implementation, and then discuss various limitations that affect this and other maximum likelihood methods. We assess the bias arising from fitting an incorrect galaxy model using simple noise-free images and find that it should not be a concern for current cosmic shear surveys. We test im3shape on the GREAT08 Challenge image simulations, and meet the requirements for upcoming cosmic shear surveys in the case that the simulations are encompassed by the fitted model, using a simple correction for image noise bias. For the fiducial branch of GREAT08 we obtain a negligible additive shear bias and sub-two percent level multiplicative bias, which is suitable for analysis of current surveys. We fall short of the sub-percent level requirement for upcoming surveys, which we attribute to a combination of noise bias and the mis-match between our galaxy model and the model used in the GREAT08 simulations. We meet the requirements for current surveys across all branches of GREAT08, except those with small or high noise galaxies, which we would cut from our analysis. Using the GREAT08 metric we we obtain a score of Q=717 for the usable branches, relative to the goal of Q=1000 for future experiments. The code is freely available from https://bitbucket.org/joezuntz/im3shape
We introduce the Galaxy Intensity Mapping cross-COrrelation estimator (GIMCO), which is a new tomographic estimator for the gravitational lensing potential, based on a combination of intensity mapping (IM) and galaxy number counts. The estimator can be written schematically as IM$(z_f)times$galaxy$(z_b)$ $-$ galaxy$(z_f)times$IM$(z_b)$ for a pair of distinct redshifts $(z_f,z_b)$; this combination allows to greatly reduce the contamination by density-density correlations, thus isolating the lensing signal. As an estimator constructed only from cross-correlations, it is additionally less susceptible to systematic effects. We show that the new estimator strongly suppresses cosmic variance and consequently improves the signal-to-noise ratio (SNR) for the detection of lensing, especially on linear scales and intermediate redshifts. %This makes it particularly valuable for future studies of dark energy and modified gravity. For cosmic variance dominated surveys, the SNR of our estimator is a factor 30 larger than the SNR obtained from the correlation of galaxy number counts only. Shot noise and interferometer noise reduce the SNR. For the specific example of the Dark Energy Survey (DES) cross-correlated with the Hydrogen Intensity mapping and Real time Analysis eXperiment (HIRAX), the SNR is around 4, whereas for Euclid cross-correlated with HIRAX it reaches 52. This corresponds to an improvement of a factor 4-5 compared to the SNR from DES alone. For Euclid cross-correlated with HIRAX the improvement with respect to Euclid alone strongly depends on the redshift. We find that the improvement is particularly important for redshifts below 1.6, where it reaches a factor of 5. This makes our estimator especially valuable to test dark energy and modified gravity, that are expected to leave an impact at low and intermediate redshifts.
Gravitational lensing of the CMB is a valuable cosmological signal that correlates to tracers of large-scale structure and acts as a important source of confusion for primordial $B$-mode polarization. State-of-the-art lensing reconstruction analyses use quadratic estimators, which are easily applicable to data. However, these estimators are known to be suboptimal, in particular for polarization, and large improvements are expected to be possible for high signal-to-noise polarization experiments. We develop a method and numerical code, $rm{LensIt}$, that is able to find efficiently the most probable lensing map, introducing no significant approximations to the lensed CMB likelihood, and applicable to beamed and masked data with inhomogeneous noise. It works by iteratively reconstructing the primordial unlensed CMB using a deflection estimate and its inverse, and removing residual lensing from these maps with quadratic estimator techniques. Roughly linear computational cost is maintained due to fast convergence of iterative searches, combined with the local nature of lensing. The method achieves the maximal improvement in signal to noise expected from analytical considerations on the unmasked parts of the sky. Delensing with this optimal map leads to forecast tensor-to-scalar ratio parameter errors improved by a factor $simeq 2 $ compared to the quadratic estimator in a CMB stage IV configuration.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
