Statistical Strong Lensing V. Inference with Microlensed Quasars

Microlensing of strongly lensed quasars provides direct information on the stellar component of lens galaxies, and can help break the degeneracy between the stellar initial mass function (IMF) and the dark matter. However, a statistical framework to forward model microlensing data while accounting for strong lensing selection effects is still missing. We incorporated microlensing effects into the hierarchical model of Paper III, and investigated to what extent flux measurements of strongly lensed quasars can break the degeneracy between the stellar IMF and the dark matter distribution. We tested the method on a mock sample of 219 doubly imaged quasars, lensed by circularly symmetric lenses. When microlensing effects are modelled, the inference using image fluxes recovers the input foreground population parameters and provides tighter constraints than the inference without flux constraints. The improvement is clearest in the $μ_α-μ_5$ plane, where $μ_α$ describes the mean stellar population synthesis mismatch parameter $\logα_{\rm SPS}$ and $μ_5$ describes the mean projected dark matter mass within $5\,{\rm kpc}$. The 68% credible interval on $μ_α$ shrinks from 0.16 dex to 0.09 dex, with the degeneracy between $μ_α$ and $μ_5$ significantly reduced. By contrast, fitting the image fluxes without accounting for microlensing leads to large biases on the inferred parameters. This work is a first step towards an accurate modelling of microlensing data. Application to real lens samples will require dropping the circular symmetry assumption and adopting more realistic models for the microlensing magnification.