No Arabic abstract
Strong gravitational lensing offers a wealth of astrophysical information on the background source it affects, provided the lensed source can be reconstructed as if it was seen in the absence of lensing. In the present work, we illustrate how sparse optimisation can address the problem. As a first step towards a full free-form lens modelling technique, we consider linear inversion of the lensed source under sparse regularisation and joint deblending from the lens light profile. The method is based on morphological component analysis, assuming a known mass model. We show with numerical experiments that representing the lens and source light using an undecimated wavelet basis allows us to reconstruct the source and to separate it from the foreground lens at the same time. Both the source and lens light have a non-analytic form, allowing for the flexibility needed in the inversion to represent arbitrarily small and complex luminous structures in the lens and source. in addition, sparse regularisation avoids over-fitting the data and does not require the use of any adaptive mesh or pixel grid. As a consequence, our reconstructed sources can be represented on a grid of very small pixels. Sparse regularisation in the wavelet domain also allows for automated computation of the regularisation parameter, thus minimising the impact of arbitrary choice of initial parameters. Our inversion technique for a fixed mass distribution can be incorporated in future lens modelling technique iterating over the lens mass parameters. The python package corresponding to the algorithms described in this article can be downloaded via the github platform at https://github.com/herjy/SLIT.
We describe a new method for analyzing gravitational lens images, for the case where the source light distribution is pixelized. The method is suitable for high resolution, high S/N data of a multiply-imaged extended source. For a given mass distribution, we show that the step of inverting the image to obtain the deconvolved pixelized source light distribution, and the uncertainties, is a linear one. This means that the only parameters of the non-linear problem are those required to model the mass distribution. This greatly simplifies the search for a min.-chi^2 fit to the data and speeds up the inversion. The method is extended in a straightforward way to include linear regularization. We apply the method to simulated Einstein ring images and demonstrate the effectiveness of the inversion for both the unregularized and regularized cases.
In this article we study the well-known strong lensing system SDSS J1004+4112. Not only does it host a large-separation lensed quasar with measured time-delay information, but several other lensed galaxies have been identified as well. A previously developed strong lens inversion procedure that is designed to handle a wide variety of constraints, is applied to this lensing system and compared to results reported in other works. Without the inclusion of a tentative central image of one of the galaxies as a constraint, we find that the model recovered by the other constraints indeed predicts an image at that location. An inversion which includes the central image provides tighter constraints on the shape of the central part of the mass map. The resulting model also predicts a central image of a second galaxy where indeed an object is visible in the available ACS images. We find masses of 2.5x10^13 M_O and 6.1x10^13 M_O within a radius of 60 kpc and 110 kpc respectively, confirming the results from other authors. The resulting mass map is compatible with an elliptical generalization of a projected NFW profile, with r_s = 58_{-13}^{+21} arcsec and c_vir = 3.91 +/- 0.74. The orientation of the elliptical NFW profile follows closely the orientation of the central cluster galaxy and the overall distribution of cluster members.
Using the combined resolving power of the Hubble Space Telescope and gravitational lensing, we resolve star-forming structures in a z~2.5 galaxy on scales much smaller than the usual kiloparsec diffraction limit of HST. SGAS J111020.0+645950.8 is a clumpy, star forming galaxy lensed by the galaxy cluster SDSS J1110+6459 at z = 0.659, with a total magnification ~30x across the entire arc. We use a hybrid parametric/non-parametric strong lensing mass model to compute the deflection and magnification of this giant arc, reconstruct the light distribution of the lensed galaxy in the source plane, and resolve the star formation into two dozen clumps. We develop a forward-modeling technique to model each clump in the source plane. We ray trace the model to the image plane, convolve with the instrumental point spread function (PSF), and compare with the GALFIT model of the clumps in the image plane, which decomposes clump structure from more extended emission. This technique has the advantage, over ray tracing, by accounting for the asymmetric lensing shear of the galaxy in the image plane and the instrument PSF. At this resolution, we can begin to study star formation on a clump-by-clump basis, toward the goal of understanding feedback mechanisms and the buildup of exponential disks at high redshift.
We present $Hubble Space Telescope$ ($HST$) imaging and grism spectroscopy of a strongly lensed LIRG at $z=0.816$, SGAS 143845.1$+$145407, and use the magnification boost of gravitational lensing to study the distribution of star formation throughout this galaxy. Based on the $HST$ imaging data, we create a lens model for this system; we compute the mass distribution and magnification map of the $z=0.237$ foreground lens. We find that the magnification of the lensed galaxy ranges between $2$ and $10$, with a total magnification (measured over all the images of the source) of $mu=11.8^{+4.6}_{-2.4}$. We find that the total projected mass density within $sim34$ kpc of the brightest cluster galaxy is $6.0^{+0.3}_{-0.7}times10^{12},M_{odot}$. Using the lens model we create a source reconstruction for SGAS 143845.1$+$145407, which paired with a faint detection of H$alpha$ in the grism spectroscopy, allows us to finally comment directly on the distribution of star formation in a $zsim1$ LIRG. We find widespread star formation across this galaxy, in agreement with the current understanding of these objects. However, we note a deficit of H$alpha$ emission in the nucleus of SGAS 143845.1$+$145407, likely due to dust extinction.
A fraction of light scalar dark matter, especially axions, may organize into Bose-Einstein condensates, gravitationally bound clumps, boson stars, and be present in large number in galactic halos today. We compute the expected number of gravitational microlensing events of clumps composed of the ordinary QCD axion and axion-like-particles and derive microlensing constraints from the EROS-2 survey and the Subaru Hyper Suprime-Cam observation. We perform a detailed lensing calculation, including the finite lens and source size effects in our analysis. We constrain the axion mass in terms of the fraction of dark matter collapsed into clumps, the individual clump densities, and the axion self-coupling. We also consider and constrain clumps composed of a generic scalar dark matter candidate with repulsive self-interactions. Our analysis opens up a new window for the potential discovery of dark matter.