Are Models of Strong Gravitational Lensing by Clusters Converging or Diverging?

TL;DR

This study investigates whether increasing the number of images in gravitational lens models leads to convergence on a true mass distribution, finding no evidence of convergence or divergence with more data.

Contribution

The paper introduces novel comparison metrics and demonstrates that more images do not necessarily improve model similarity, challenging assumptions about model convergence.

Findings

Models with more images are not more similar to each other.

Model similarity does not increase with the number of input images.

Current models may be limited by degeneracies and unexploited priors.

Abstract

The increasingly large numbers of multiple images in cluster-scale gravitational lenses have allowed for tighter constraints on the mass distributions of these systems. Most lens models have progressed alongside this increase in image number. The general assumption is that these improvements would result in lens models converging to a common solution, suggesting that models are approaching the true mass distribution. To test whether or not this is occurring, we examine a sample of lens models of MACS J0416.1$-$2403 containing varying number of images as input. Splitting the sample into two bins (those including $<150$ and $>150$ images), we quantify the similarity of models in each bin using three comparison metrics, two of which are novel: Median Percent Difference, Frechet Distance, and Wasserstein Distance. In addition to quantifying similarity, the Frechet distance metric seems to also be an indicator of the mass sheet degeneracy. Each metric indicates that models with a greater number of input images are no more similar between one another than models with fewer input images. This suggests that lens models are neither converging nor diverging to a common solution for this system, regardless of method. With this result, we suggest that future models more carefully investigate lensing degeneracies and anomalous mass clumps (mass features significantly displaced from baryonic counterparts) to rigorously evaluate their model's validity. We also recommend further study into alternative, underutilized lens model priors (e.g. flux ratios) as an additional input constraint to image positions in hopes of breaking existing degeneracies.