Full-Sky Weak Lensing Simulation with 70 Billion Particles

Share this post on:


Authors: R. Teyssier, S. Pires,  ... , J.-L. Starck et al.
Journal: A&A
Year: 2009
Download: ADS | arXiv


We have performed a 70 billion dark-matter particles N-body simulation in a 2 h1 Gpc periodic box, using the concordance, cosmological model as favored by the latest WMAP3 results. We have computed a full-sky convergence map with a resolution of Δθ0.74 arcmin2, spanning 4 orders of magnitude in angular dynamical range. Using various high-order statistics on a realistic cut sky, we have characterized the transition from the linear to the nonlinear regime at 1000 and shown that realistic galactic masking affects high-order moments only below <200. Each domain (Gaussian and non-Gaussian) spans 2 decades in angular scale. This map is therefore an ideal tool for testing map-making algorithms on the sphere. As a first step in addressing the full map reconstruction problem, we have benchmarked in this paper two denoising methods: 1) Wiener filtering applied to the Spherical Harmonics decomposition of the map and 2) a new method, called MRLens, based on the modification of the Maximum Entropy Method on a Wavelet decomposition. While the latter is optimal on large spatial scales, where the signal is Gaussian, MRLens outperforms the Wiener method on small spatial scales, where the signal is highly non-Gaussian. The simulated full-sky convergence map is freely available to the community to help the development of new map-making algorithms dedicated to the next generation of weak-lensing surveys.

Share this post on:

Author: Samuel Farrens

I have been a postdoctoral researcher at CEA Saclay since October 2015. I am currently working on the DEDALE project and the Euclid mission with Jean-Luc Starck.

My background is in optical detection of clusters of galaxies and photometric redshift estimation. I am now branching out into the field of PSF estimation using sparse signal processing techniques.

View all posts by Samuel Farrens >