Towards noiseless gravitational lensing simulations

The microphysical properties of the dark matter (DM) particle can, in principle, be constrained by the properties and abundance of substructures in galaxy clusters, as measured through strong gravitational lensing. Unfortunately, there is a lack of accurate theoretical predictions for the lensing signal of these substructures, mainly because of the discreteness noise inherent to N-body simulations. Here, we present a method, dubbed as Recursive-TCM, that is able to provide lensing predictions with an arbitrarily low discreteness noise. This solution is based on a novel way of interpreting the results of N-body simulations, where particles simply trace the evolution and distortion of Lagrangian phase-space volume elements. We discuss the advantages and limitations of this method compared to the widely used density estimators based on cloud-in-cells and adaptive-kernel smoothing. Applying the new method to a cluster-sized DM halo simulated in warm and cold DM scenarios, we show how the expected differences in their substructure population translate into differences in convergence and magnification maps. We anticipate that our method will provide the high-precision theoretical predictions required to interpret and fully exploit strong gravitational lensing observations.