[Paper Review] Nonlinear Weak Lensing reconstruction for Galaxy Clusters

Yuan Shi, Li Cui|arXiv (Cornell University)|Jan 6, 2026

Advanced Fluorescence Microscopy Techniques0 citations

TL;DR

The paper develops a nonlinear weak-lensing mass reconstruction framework for galaxy clusters that handles masked regions and uses a model-based initial guess to improve convergence in the nonlinear regime, validated with simulated data.

ABSTRACT

We present a numerical investigation of nonlinear cluster lens reconstruction using weak lensing mass mapping. Recent advances in imaging and shear estimation have pushed reliable reduced shear measurements closer to cluster cores, making mass reconstruction accessible in the nonlinear regime. However, the Kaiser-Squires based algorithm becomes unstable in cluster cores, where convergence $κ$ significantly deviates from zero and the linear approximation breaks down. To address this limitation, we develop a reconstruction framework with two key modifications: applying smooth masks to these regions and using a model-derived analytical solution as the initial guess, rather than assuming $κ= 0$. We validate our framework using simulated cluster lensing data with known mass distributions, incorporating realistic masks that arise from limitations in reduced shear measurements. We show that in the absence of shape noise, our framework yields high-fidelity mass reconstruction in regions of large reduced shear, with the best-performing method achieving residuals below $0.02 σ$ in the unmasked regions. This pushes mass reconstruction to higher accuracy in the nonlinear regime.

Motivation & Objective

• Motivate accurate mass reconstruction in cluster cores where κ is non-negligible and the KS approach fails.
• Develop a framework that mitigates masking biases by using a model-derived initial guess and a smooth mask.
• Quantify performance of KS and AKRA-based reconstructions under realistic masks and reduced-shear data limitations.
• Demonstrate improved reconstruction fidelity in nonlinear regimes using simulated cluster data with known mass distributions.

Proposed method

• Review KS and AKRA formulations for κ and γ, including reduced shear g and its relation to κ.
• Introduce a model-based initial guess for κ, using a Singular Isothermal Sphere (SIS) to initialize κ^(0).
• Replace binary masks with a smooth masking function to reduce numerical instabilities.
• Formulate the masked shear as a linear system γ^m = A κ + n with A encoding lensing response and mask convolution.
• Solve with AKRA estimator: κ̂ = (A^T N^{-1} A + R)^{-1} A^T N^{-1} γ^m, iterating with g^m(1−κ^(i−1)).
• Evaluate five configurations combining KS/AKRA, K1/K2, A1/A2/A3, with varying initial guesses and mask forms.

Experimental results

Research questions

• RQ1 How do masks from unreliable reduced-shear measurements bias nonlinear mass reconstruction in cluster cores?
• RQ2 Can a model-based initial κ and smooth masking improve convergence and reduce residuals in nonlinear weak-lensing maps?
• RQ3 What is the impact of iteration count on residuals and bias in KS and AKRA-based reconstructions under realistic masking?
• RQ4 What are the limits of mass-mapping accuracy as the amount of available reduced-shear data is varied?

Key findings

• AKRA-based methods dramatically reduce bias in unmasked regions compared to KS-based methods (mean residuals ~0.003 vs ~0.12 in toy model).
• A model-based initial guess (SIS) improves stability and reduces variance over iterations compared to κ^(0)=0.
• Smooth masking provides greater numerical stability and lowers residuals at convergence than binary masking, especially in nonlinear regimes.
• Convergence can occur with low residuals even at high masking thresholds (g_th up to 0.8), though edge effects and iteration stopping influence bias.
• K2 (KS with improved initialization) and A2 (AKRA with model-based initialization) offer favorable trade-offs between efficiency and accuracy.
• A3 (AKRA with optimized settings) achieves the highest stability and accuracy, particularly when reduced-shear data are strongly nonlinear.