2024MNRAS.527.2639P

SDSS J1004+4112 is a well-studied gravitational lens with a recently measured time delay between its first and fourth arriving quasar images. Using this new constraint, we present updated free-form lens reconstructions using the lens inversion method GRALE, which only uses multiple image and time delay data as inputs. In addition, we obtain hybrid lens reconstructions by including a model of the brightest cluster galaxy as a Sersic lens. For both reconstructions, we use two sets of images as input: one with all identified images, and the other a revised set leaving out images that have been potentially misidentified. We also develop a source position optimization Markov chain Monte Carlo (MCMC) routine, performed on completed GRALE runs, that allows each model to better match observed image positions and time delays. All the reconstructions produce similar mass distributions, with the hybrid models finding a steeper profile in the centre. Similarly, all the mass distributions are fitted by the Navarro-Frenk-White profile, finding results consistent with previous parametric reconstructions and those derived from Chandra X-ray observations. We identify an ∼5 × 10¹¹ M_⊙ substructure apparently unaffiliated with any cluster member galaxy and present in all our models, and study its reality. Using our free-form and hybrid models, we predict a central quasar image time delay of ∼2980 ± 270 and ∼3280 ± 215 d, respectively. A potential future measurement of this time delay will, while being an observational challenge, further constrain the steepness of the central density profile.