Cluster strong lensing with hierarchical inference: Formalism, functional tests, and public code release

Context. Lensing by galaxy clusters is a versatile probe of cosmology and extragalactic astrophysics, but the accuracy of some of its predictions is limited by the simplified models adopted to reduce the (otherwise intractable) number of degrees of freedom.Aims. We aim to explore cluster lensing models in which the parameters of all cluster member galaxies are free to vary around some common scaling relations with non-zero scatter and deviate significantly from these relations if, and only if, the data require this.Methods. We devised a Bayesian hierarchical inference framework that enables the determination of all lensing parameters and the scaling relation hyperparameters, including intrinsic scatter, from lensing constraints and (if given) stellar kinematic measurements. We achieve this through BAYESLENS, a purpose-built wrapper around common parametric lensing codes that can sample the full posterior on parameters and hyperparameters; we release BAYESLENS with this paper.Results. We ran functional tests of our code against simple mock cluster lensing datasets with realistic uncertainties. The parameters and hyperparameters are recovered within their 68% credibility ranges and the positions of all the "observed" multiple images are accurately reproduced by the BAYELENS best-fit model, without over-fitting.Conclusions. We have shown that an accurate description of cluster member galaxies is attainable, despite a large number of degrees of freedom, through fast and tractable inference. This extends beyond the most updated cluster lensing models. The precise impact on studies of cosmography, galaxy evolution, and high-redshift galaxy populations can then be quantified on real galaxy clusters. While other sources of systematics exist and may be significant in real clusters, our results show that the contribution of intrinsic scatter in cluster member populations can now be controlled.