Showing 1–50 of 105 results for author: Williams, L L R

Elucidating Galaxy Population Properties Using a Model-Free Analysis of Quadruply Imaged Quasar Lenses From Large Surveys

Authors: John Miller Jr, Liliya L. R. Williams

Abstract: The population of strong lensing galaxies is a sub-set of intermediate-redshift massive galaxies, whose population-level properties are not yet well understood. In the near future, thousands of multiply imaged systems are expected to be discovered by wide-field surveys like Rubin Observatory's Legacy Survey of Space and Time (LSST) and Euclid. With the soon-to-be robust population of quadruply len… ▽ More The population of strong lensing galaxies is a sub-set of intermediate-redshift massive galaxies, whose population-level properties are not yet well understood. In the near future, thousands of multiply imaged systems are expected to be discovered by wide-field surveys like Rubin Observatory's Legacy Survey of Space and Time (LSST) and Euclid. With the soon-to-be robust population of quadruply lensed quasars, or quads, in mind, we introduce a novel technique to elucidate the empirical distribution of the galaxy population properties. Our re-imagining of the prevailing strong lensing analysis does not fit mass models to individual lenses, but instead starts with parametric models of many galaxy populations, which include generally ignored mass distribution complexities and exclude external shear for now. We construct many mock galaxy populations with different properties and obtain populations of quads from each of them. The mock `observed' population of quads is then compared to those from the mocks using a model-free analysis based on a 3D sub-space of directly observable quad image properties. The distance between two quad populations in the space of image properties is measured by a metric $η$, and the distance between their parent galaxy populations in the space of galaxy properties is measured by $ζ$. We find a well defined relation between $η$ and $ζ$. The discovered relation between the space of image properties and the space of galaxy properties allows for the observed galaxy population properties to be estimated from the properties of their quads, which will be conducted in a future paper. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: 18 pages and 13 figures

A high-resolution view of the source-plane magnification near cluster caustics in wave dark matter models

Authors: Jose M. Diego, Alfred Amruth, Jose M. Palencia, Tom Broadhurst, Sung Kei Li, Jeremy Lim, Rogier A. Windhorst, Adi Zitrin, Alexei V. Filippenko, Liliya L. R. Williams, Ashish K. Meena, Wenlei Chen, Patrick L. Kelly

Abstract: We present the highest resolution images to date of caustics formed by wave dark matter ($ψ$DM) fluctuations near the critical curves of cluster gravitational lenses. We describe the basic magnification features of $ψ$DM in the source plane at high macromodel magnification and discuss specific differences between the $ψ$DM and standard cold dark matter (CDM) models. The unique generation of demagn… ▽ More We present the highest resolution images to date of caustics formed by wave dark matter ($ψ$DM) fluctuations near the critical curves of cluster gravitational lenses. We describe the basic magnification features of $ψ$DM in the source plane at high macromodel magnification and discuss specific differences between the $ψ$DM and standard cold dark matter (CDM) models. The unique generation of demagnified counterimages formed outside the Einstein radius for $ψ$DM is highlighted. Substructure in CDM cannot generate such demagnified images of positive parity, thus providing a definitive way to distinguish $ψ$DM from CDM. Highly magnified background sources with sizes $r\approx 1pc$, or approximately a factor of ten smaller than the expected de Broglie wavelength of $ψ$DM, offer the best possibility of discriminating between $ψ$DM and CDM. These include objects such as very compact stellar clusters at high redshift that JWST is finding in abundance. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 8 pages

Dark Matter distinguished by skewed microlensing in the "Dragon Arc"

Authors: Tom Broadhurst, Sung Kei Li, Amruth Alfred, Jose M. Diego, Paloma Morilla, Patrick L. Kelly, Fengwu Sun, Masamune Oguri, Hayley Williams, Rogier Windhorst, Adi Zitrin, Katsuya T. Abe, Wenlei Chen, Yoshinobu Fudamoto, Hiroki Kawai, Jeremy Lim, Tao Liu, Ashish K. Meena, Jose M. Palencia, George F. Smoot, Liliya L. R. Williams

Abstract: Microlensed stars recently discovered by JWST & HST follow closely the winding critical curve of A370 along all sections of the ``Dragon Arc" traversed by the critical curve. These transients are fainter than $m_{AB}>26.5$, corresponding to the Asymptotic Giant Branch (AGB) and microlensed by diffuse cluster stars observed with $\simeq 18M_\odot/pc^2$, or about $\simeq 1$\% of the projected dark m… ▽ More Microlensed stars recently discovered by JWST & HST follow closely the winding critical curve of A370 along all sections of the ``Dragon Arc" traversed by the critical curve. These transients are fainter than $m_{AB}>26.5$, corresponding to the Asymptotic Giant Branch (AGB) and microlensed by diffuse cluster stars observed with $\simeq 18M_\odot/pc^2$, or about $\simeq 1$\% of the projected dark matter density. Most microlensed stars appear along the inner edge of the critical curve, following an asymmetric band of width $\simeq 4$kpc that is skewed by $-0.7\pm0.2$kpc. Some skewness is expected as the most magnified images should form along the inner edge of the critical curve with negative parity, but the predicted shift is small $\simeq -0.04$kpc and the band of predicted detections is narrow, $\simeq 1.4$kpc. Adding CDM-like dark halos of $10^{6-8}M_\odot$ broadens the band as desired but favours detections along the outer edge of the critical curve, in the wrong direction, where sub-halos generate local Einstein rings. Instead, the interference inherent to ``Wave Dark Matter" as a Bose-Einstein condensate ($ψ$DM) forms a symmetric band of critical curves that favours negative parity detections. A de Broglie wavelength of $\simeq 10$pc matches well the observed $4$kpc band of microlenses and predicts negative skewness $\simeq -0.6$kpc, similar to the data. The implied corresponding boson mass is $\simeq 10^{-22}$eV, in good agreement with estimates from dwarf galaxy cores when scaled by momentum. Further JWST imaging may reveal the pattern of critical curves by simply ``joining the dots" between microlensed stars, allowing wave corrugations of $ψ$DM to be distinguished from CDM sub-halos △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 12 pages, 5 figures

Flashlights: Microlensing vs Stellar Variability of Transients in the Star Clusters of the Dragon Arc

Authors: Sung Kei Li, Patrick L. Kelly, Jose M. Diego, Jeremy Lim, WenLei Chen, Amruth Alfred, Liliya L. R. Williams, Thomas J. Broadhurst, Ashish. K. Meena, Adi Zitrin, Alex Chow

Abstract: We study the nature of transient events detected in the "Dragon Arc", a star-forming galaxy at a redshift of $0.7251$ that is gravitationally lensed by the galaxy cluster Abell 370. In particular, we focus on a subset of ten transients that are identified as unresolved young star clusters in the deep broadband, F200LP, taken as part of the "Flashlights" Hubble Space Telescope program, showing flux… ▽ More We study the nature of transient events detected in the "Dragon Arc", a star-forming galaxy at a redshift of $0.7251$ that is gravitationally lensed by the galaxy cluster Abell 370. In particular, we focus on a subset of ten transients that are identified as unresolved young star clusters in the deep broadband, F200LP, taken as part of the "Flashlights" Hubble Space Telescope program, showing flux variations of $\sim 10-20\%$ over a period of about a year. Here we develop several methods to address whether stellar microlensing alone is capable of explaining the transients, or whether intrinsic stellar outbursts or variability are required to explain them. We first present a lens model that has new constraints in the Dragon Arc itself to understand the properties of the lensed young star clusters. Using our improved galaxy-cluster lens model, we simulate the effect of microlensing on the flux variation for unresolved stars within lensed young star clusters. We find good agreement between the observed and the expected detection rates of microlensing events by intracluster stars of young star clusters within $1σ$. However, we cannot fully exclude the possibility that a minority of these transients are caused by intrinsic stellar variability such as outbursts of Luminous Blue Variables (LBVs). With JWST observations taken recently or coming in the near future, the color information will be able to break the degeneracy and definitively test whether or not these lensed young star cluster transients are caused by stellar microlensing. △ Less

Submitted 12 April, 2024; originally announced April 2024.

Comments: 28 pages, 14 figures. To be submitted, comments welcomed

Imaging dark matter at the smallest scales with $z\approx1$ lensed stars

Authors: J. M. Diego, Sung Kei Li, Alfred Amruth, Ashish K. Meena, Tom J. Broadhurst, Patrick L. Kelly, Alexei V. Filippenko, Liliya L. R. Williams, Adi Zitrin, William E. Harris, Marta Reina-Campos, Carlo Giocoli, Liang Dai, Mitchell F. Struble, Tommaso Treu, Yoshinobu Fudamoto, Daniel Gilman, Anton M. Koekemoer, Jeremy Lim, J. M. Palencia, Fengwu Sun, Rogier A. Windhorst

Abstract: Observations of caustic-crossing galaxies at redshift $0.7<z<1$ show a wealth of transient events. Most of them are believed to be microlensing events of highly magnified stars. Earlier work predicted such events should be common near the critical curves (CCs) of galaxy clusters, but some are found relatively far away from these CCs. We consider the possibility that substructure on milliarcsecond… ▽ More Observations of caustic-crossing galaxies at redshift $0.7<z<1$ show a wealth of transient events. Most of them are believed to be microlensing events of highly magnified stars. Earlier work predicted such events should be common near the critical curves (CCs) of galaxy clusters, but some are found relatively far away from these CCs. We consider the possibility that substructure on milliarcsecond scales (few parsecs in the lens plane) is boosting the microlensing signal. We study the combined magnification from the macrolens, millilenses, and microlenses (3M-lensing). After considering realistic populations of millilenses and microlenses, we conclude that the enhanced microlensing rate around millilenses is not sufficient to explain the high fraction of observed events in the far region. Instead we find a that the shape of the luminosity function (LF) of the lensed stars combined with the amount of substructure in the lens plane determines the number of mcirolensing events found near and far from the CC. By measuring $β$ (the exponent of the LF), and the number density of microlensing events at each location, one can create a pseudoimage of the underlying distribution of mass on small scales. We identify two regimes: (i) positive imaging regime where $β>2$ and the number density of events is greater around substructureand the number density of events is greater around substructures, and (ii) negative imaging regime where $β<2$. We study the particular case of seven microlensing events found by HST in the Dragon arc (at z=0.725). We find that a population of supergiant stars with a steep LF with $β=2.55$ fits the distribution of these events. We identify a small region of high density of microlensing events, and interpret it as evidence of a possible invisible substructure, for which we derive a mass of $\sim 1.3 \times 10^8\,\Msun$ (within its Einstein radius). △ Less

Submitted 22 April, 2024; v1 submitted 11 April, 2024; originally announced April 2024.

Comments: 23 pages, 16 figures

Extrapolating the projected potential of gravitational lens models: property-preserving degeneracies

Authors: Jori Liesenborgs, Derek Perera, Liliya L. R. Williams

Abstract: While gravitational lens inversion holds great promise to reveal the structure of the light-deflecting mass distribution, both light and dark, the existence of various kinds of degeneracies implies that care must be taken when interpreting the resulting lens models. This article illustrates how thinking in terms of the projected potential helps to gain insight into these matters. Additionally it i… ▽ More While gravitational lens inversion holds great promise to reveal the structure of the light-deflecting mass distribution, both light and dark, the existence of various kinds of degeneracies implies that care must be taken when interpreting the resulting lens models. This article illustrates how thinking in terms of the projected potential helps to gain insight into these matters. Additionally it is shown explicitly how, when starting from a discretised version of the projected potential of one particular lens model, the technique of quadratic programming can be used to create a multitude of equivalent lens models that preserve all or a subset of lens properties. This method is applied to a number of scenarios, showing the lack of grasp on the mass outside the strong lensing region, revisiting mass redistribution in between images and applying this to a recent model of the SDSS J1004+4112 cluster, as well as illustrating the generalised mass sheet degeneracy and source-position transformation. In the case of J1004 we show that this mass redistribution did not succeed at completely eliminating a dark mass clump recovered by GRALE near one of the quasar images. △ Less

Submitted 27 February, 2024; originally announced February 2024.

Comments: 12 pages, 9 figures, MNRAS accepted

Essentials of strong gravitational lensing

Authors: Prasenjit Saha, Dominique Sluse, Jenny Wagner, Liliya L. R. Williams

Abstract: Of order one in 10^3 quasars and high-redshift galaxies appears in the sky as multiple images as a result of gravitational lensing by unrelated galaxies and clusters that happen to be in the foreground. While the basic phenomenon is a straightforward consequence of general relativity, there are many non-obvious consequences that make multiple-image lensing systems (aka strong gravitational lenses)… ▽ More Of order one in 10^3 quasars and high-redshift galaxies appears in the sky as multiple images as a result of gravitational lensing by unrelated galaxies and clusters that happen to be in the foreground. While the basic phenomenon is a straightforward consequence of general relativity, there are many non-obvious consequences that make multiple-image lensing systems (aka strong gravitational lenses) remarkable astrophysical probes in several different ways. This article is an introduction to the essential concepts and terminology in this area, emphasizing physical insight. The key construct is the Fermat potential or arrival-time surface: from it the standard lens equation, and the notions of image parities, magnification, critical curves, caustics, and degeneracies all follow. The advantages and limitations of the usual simplifying assumptions (geometrical optics, small angles, weak fields, thin lenses) are noted, and to the extent possible briefly, it is explained how to go beyond these. Some less well-known ideas are discussed at length: arguments using wavefronts show that much of the theory carries over unchanged to the regime of strong gravitational fields; saddle-point contours explain how even the most complicated image configurations are made up of just two ingredients. Orders of magnitude, and the question of why strong lensing is most common for objects at cosmological distance, are also discussed. The challenges of lens modeling, and diverse strategies developed to overcome them, are discussed in general terms, without many technical details. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: To appear in Space Science Reviews, Topical Collection "Strong Gravitational Lensing", eds. J. Wambsganss et al

Free-Form and Hybrid Lens Models for SDSS J1004+4112: Substructure and Central Image Time Delay Constraints

Authors: Derek Perera, Liliya L. R. Williams, Jori Liesenborgs, Agniva Ghosh, Prasenjit Saha

Abstract: 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 {\tt GRALE}, which only uses multiple image and time delay data as inputs. In addition, we obtain hybrid lens reconstructions by including a model of… ▽ More 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 {\tt 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 (BCG) 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 MCMC routine, performed on completed {\tt 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 center. Similarly, all the mass distributions are fit by the Navarro-Frenk-White (NFW) profile, finding results consistent with previous parametric reconstructions and those derived from Chandra X-ray observations. We identify a $\sim 5 \times 10^{11} M_{\odot}$ 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 $\sim 2980 \pm 270$ and $\sim 3280 \pm 215$ days, respectively. A potential future measurement of this time delay will, while being an observational challenge, further constrain the steepness of the central density profile. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 13 pages, 5 figures, MNRAS Accepted

What Multiple Images Say About the Large-Scale Mass Maps of Galaxy Clusters

Authors: Kekoa Lasko, Liliya L. R. Williams, Agniva Ghosh

Abstract: All lens modeling methods, simply-parametrized, hybrid, and free-form, use assumptions to reconstruct galaxy clusters with multiply imaged sources, though the nature of these assumptions (priors) can differ considerably between methods. This raises an important question in strong lens modeling: how much information about the mass model comes from the lensed images themselves, and how much is a con… ▽ More All lens modeling methods, simply-parametrized, hybrid, and free-form, use assumptions to reconstruct galaxy clusters with multiply imaged sources, though the nature of these assumptions (priors) can differ considerably between methods. This raises an important question in strong lens modeling: how much information about the mass model comes from the lensed images themselves, and how much is a consequence of model priors. One way to assess the relative contributions of the lensing data vs. model priors is to estimate global lens properties through images alone, without any prior assumptions about the mass distribution. This is our approach. We use 200 mock cluster lenses, half of which have substructures which vary from clumpy and compact to smooth and extended; a simulated cluster Ares; and real clusters Abell 1689 and RXJ1347.5-1145 to show that the center, ellipticity, and position angle can be estimated quite well, and nearly perfectly for weakly substructured clusters, implying that the recovery of these properties is largely driven by the images, not priors. However, the correlation between the true and image-estimated amount of substructure has a lot of scatter, suggesting that multiple images do not uniquely constrain substructure. Therefore in general, lens model priors have a stronger effect on smaller scales. Our analysis partly explains why reconstructions using different methodologies can produce qualitatively different mass maps on substructure scales. Our analysis is not meant to aide or replace lens inversion methods, but only to investigate what cluster properties are constrained with multiple images. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 16 Pages, 14, MNRAS

Journal ref: MNRAS, Volume 525, Issue 4 (2023)

Flashlights: Properties of Highly Magnified Images Near Cluster Critical Curves in the Presence of Dark Matter Subhalos

Authors: Liliya L. R. Williams, Patrick L. Kelly, Tommaso Treu, Alfred Amruth, Jose M. Diego, Sung Kei Li, Ashish K. Meena, Adi Zitrin, Thomas J. Broadhurst, Alexei V. Filippenko

Abstract: Dark matter subhalos with extended profiles and density cores, and globular stars clusters of mass $10^6-10^8 M_\odot$, that live near the critical curves in galaxy cluster lenses can potentially be detected through their lensing magnification of stars in background galaxies. In this work we study the effect such subhalos have on lensed images, and compare to the case of more well studied microlen… ▽ More Dark matter subhalos with extended profiles and density cores, and globular stars clusters of mass $10^6-10^8 M_\odot$, that live near the critical curves in galaxy cluster lenses can potentially be detected through their lensing magnification of stars in background galaxies. In this work we study the effect such subhalos have on lensed images, and compare to the case of more well studied microlensing by stars and black holes near critical curves. We find that the cluster density gradient and the extended mass distribution of subhalos are important in determining image properties. Both lead to an asymmetry between the image properties on the positive and negative parity sides of the cluster that is more pronounced than in the case of microlensing. For example, on the negative parity side, subhalos with cores larger than about $50\,$pc do not generate any images with magnification above $\sim 100$ outside of the immediate vicinity of the cluster critical curve. We discuss these factors using analytical and numerical analysis, and exploit them to identify observable signatures of subhalos: subhalos create pixel-to-pixel flux variations of $\gtrsim 0.1$ magnitudes, on the positive parity side of clusters. These pixels tend to cluster around (otherwise invisible) subhalos. Unlike in the case of microlensing, signatures of subhalo lensing can be found up to $1''$ away from the critical curves of massive clusters. △ Less

Submitted 19 December, 2023; v1 submitted 12 April, 2023; originally announced April 2023.

Comments: ApJ, accepted, 22 pages, 17 figures

Constraints on the Inner Regions of Lensing Galaxies from Central Images using a Recent AGN Offset Distribution

Authors: Derek Perera, Liliya L. R. Williams, Claudia Scarlata

Abstract: In gravitational lensing, central images in quads can serve as a powerful probe of the inner regions of lens galaxies. The presence of an offset central supermassive black hole (SMBH) has the potential to distort the time-delay surface in a way such that 3 central images form: a strongly de-magnified image near the SMBH, and two less de-magnified (and potentially observable) images at a central ma… ▽ More In gravitational lensing, central images in quads can serve as a powerful probe of the inner regions of lens galaxies. The presence of an offset central supermassive black hole (SMBH) has the potential to distort the time-delay surface in a way such that 3 central images form: a strongly de-magnified image near the SMBH, and two less de-magnified (and potentially observable) images at a central maximum and saddle point. Using a quad lens macro model, we simulate the constraints that could be placed on various lens galaxy parameters based on their central images probability of detection or non-detection. Informed by a recent low-redshift distribution of off-nucleus AGN, we utilize Bayesian inference to constrain the mean SMBH off-nucleus distance and galactic core radius for a sample of 6 quads. In general, we find that a detection of the central image in any quad would favor larger SMBH off-nucleus distances and galaxy core sizes. Assuming a linear relationship between core radii and velocity dispersion $r_c = bσ$, these results similarly imply strong constraints on $b$, where the likely case of a central image non-detection in each quad constraining $b$ to $3.11^{+2.72}_{-2.26} \times 10^{-4}$ kpc km$^{-1}$ s. Our results show that tight constraints on lens galaxy parameters can be made regardless of a detection or non-detection of a central image. Therefore, we recommend observational searches for the central image, possibly using our suggested novel detection technique in UV filters, to formalize stronger constraints on lens galaxy parameters. △ Less

Submitted 7 April, 2023; originally announced April 2023.

Comments: MNRAS, Accepted, 20 pages, 15 figures

Flashlights: More than A Dozen High-Significance Microlensing Events of Extremely Magnified Stars in Galaxies at Redshifts z=0.7-1.5

Authors: Patrick L. Kelly, Wenlei Chen, Amruth Alfred, Thomas J. Broadhurst, Jose M. Diego, Najmeh Emami, Alexei V. Filippenko, Allison Keen, Sung Kei Li, Jeremy Lim, Ashish K. Meena, Masamune Oguri, Claudia Scarlata, Tommaso Treu, Hayley Williams, Liliya L. R. Williams, Rui Zhou, Adi Zitrin, Ryan J. Foley, Saurabh W. Jha, Nick Kaiser, Vihang Mehta, Steven Rieck, Laura Salo, Nathan Smith , et al. (1 additional authors not shown)

Abstract: Once only accessible in nearby galaxies, we can now study individual stars across much of the observable universe aided by galaxy-cluster gravitational lenses. When a star, compact object, or multiple such objects in the foreground galaxy-cluster lens become aligned, they can magnify a background individual star, and the timescale of a magnification peak can limit its size to tens of AU. The numbe… ▽ More Once only accessible in nearby galaxies, we can now study individual stars across much of the observable universe aided by galaxy-cluster gravitational lenses. When a star, compact object, or multiple such objects in the foreground galaxy-cluster lens become aligned, they can magnify a background individual star, and the timescale of a magnification peak can limit its size to tens of AU. The number and frequency of microlensing events therefore opens a window into the population of stars and compact objects, as well as high-redshift stars. To assemble the first statistical sample of stars in order to constrain the initial mass function (IMF) of massive stars at redshift z=0.7-1.5, the abundance of primordial black holes in galaxy-cluster dark matter, and the IMF of the stars making up the intracluster light, we are carrying out a 192-orbit program with the Hubble Space Telescope called "Flashlights," which is now two-thirds complete owing to scheduling challenges. We use the ultrawide F200LP and F350LP long-pass WFC3 UVIS filters and conduct two 16-orbit visits separated by one year. Having an identical roll angle during both visits, while difficult to schedule, yields extremely clean subtraction. Here we report the discovery of more than a dozen bright microlensing events, including multiple examples in the famous "Dragon Arc" discovered in the 1980s, as well as the "Spocks" and "Warhol" arcs that have hosted already known supergiants. The ultradeep observer-frame ultraviolet-through-optical imaging is sensitive to hot stars, which will complement deep James Webb Space Telescope infrared imaging. We are also acquiring Large Binocular Telescope LUCI and Keck-I MOSFIRE near-infrared spectra of the highly magnified arcs to constrain their recent star-formation histories. △ Less

Submitted 4 November, 2022; originally announced November 2022.

Flashlights: An Off-Caustic Lensed Star at Redshift $z$ = 1.26 in Abell 370

Authors: Ashish Kumar Meena, Wenlei Chen, Adi Zitrin, Patrick L. Kelly, Miriam Golubchik, Rui Zhou, Amruth Alfred, Tom Broadhurst, Jose M. Diego, Masamune Oguri, Liliya L. R. Williams, Alexei V. Filippenko, Sung Kei Li

Abstract: We report the discovery of a transient seen in a strongly lensed arc at redshift $z_{\rm s}=1.2567$ in \emph{Hubble Space Telescope} imaging of the Abell 370 galaxy cluster. The transient is detected at $29.51\pm0.14$ AB mag in a WFC3/UVIS F200LP difference image made using observations from two different epochs, obtained in the framework of the \emph{Flashlights} program, and is also visible in t… ▽ More We report the discovery of a transient seen in a strongly lensed arc at redshift $z_{\rm s}=1.2567$ in \emph{Hubble Space Telescope} imaging of the Abell 370 galaxy cluster. The transient is detected at $29.51\pm0.14$ AB mag in a WFC3/UVIS F200LP difference image made using observations from two different epochs, obtained in the framework of the \emph{Flashlights} program, and is also visible in the F350LP band ($m_{\rm F350LP} \approx 30.53\pm0.76$ AB mag). The transient is observed on the negative-parity side of the critical curve at a distance of $\sim 0.6"$ from it, greater than previous examples of lensed stars. The large distance from the critical curve yields a significantly smaller macromagnification, but our simulations show that bright, O/B-type supergiants can reach sufficiently high magnifications to be seen at the observed position and magnitude. In addition, the observed transient image is a trailing image with an observer-frame time delay of $\sim+0.8$ days from its expected counterpart, so that any transient lasting for longer than that should have also been seen on the minima side and is thus excluded. This, together with the blue colour we measure for the transient ($m_{\rm F200LP} - m_{\rm F350LP} \approx [-0.3,-1.6]$ AB), rules out most other transient candidates such as (kilo)novae, for example, and makes a lensed star the prime candidate. Assuming the transient is indeed a lensed star as suggested, many more such events should be detected in the near future in cluster surveys with the \emph{Hubble Space Telescope} and \emph{James Webb Space Telescope}. △ Less

Submitted 5 April, 2023; v1 submitted 2 November, 2022; originally announced November 2022.

Comments: 9 pages. 5 figures. Accepted in MNRAS

Strong Lensing by Galaxies

Authors: A. J. Shajib, G. Vernardos, T. E. Collett, V. Motta, D. Sluse, L. L. R. Williams, P. Saha, S. Birrer, C. Spiniello, T. Treu

Abstract: Strong gravitational lensing at the galaxy scale is a valuable tool for various applications in astrophysics and cosmology. The primary uses of galaxy-scale lensing are to study elliptical galaxies' mass structure and evolution, constrain the stellar initial mass function, and measure cosmological parameters. Since the discovery of the first galaxy-scale lens in the 1980s, this field has made sign… ▽ More Strong gravitational lensing at the galaxy scale is a valuable tool for various applications in astrophysics and cosmology. The primary uses of galaxy-scale lensing are to study elliptical galaxies' mass structure and evolution, constrain the stellar initial mass function, and measure cosmological parameters. Since the discovery of the first galaxy-scale lens in the 1980s, this field has made significant advancements in data quality and modeling techniques. In this review, we describe the most common methods for modeling lensing observables, especially imaging data, as they are the most accessible and informative source of lensing observables. We then summarize the primary findings from the literature on the astrophysical and cosmological applications of galaxy-scale lenses. We also discuss the current limitations of the data and methodologies and provide an outlook on the expected improvements in both areas in the near future. △ Less

Submitted 1 May, 2024; v1 submitted 19 October, 2022; originally announced October 2022.

Comments: Submitted to Space Science Reviews, Topical Collection "Strong Gravitational Lensing", eds. J. Wambsganss et al. This version: updated after peer review

Statistical Mechanics of Collisionless Orbits. V. The approach to equilibrium for idealized self-gravitating systems

Authors: Liliya L. R. Williams, Jens Hjorth

Abstract: Self-gravitating Newtonian systems consisting of a very large number of particles have generally defied attempts to describe them using statistical mechanics. This is paradoxical since many astronomical systems, or simulations thereof, appear to have universal, equilibrium structures for which no physical basis exist. A decade ago we showed that extremizing the number of microstates with a given e… ▽ More Self-gravitating Newtonian systems consisting of a very large number of particles have generally defied attempts to describe them using statistical mechanics. This is paradoxical since many astronomical systems, or simulations thereof, appear to have universal, equilibrium structures for which no physical basis exist. A decade ago we showed that extremizing the number of microstates with a given energy per unit mass, under the constraints of conserved total energy and mass, leads to the maximum entropy state, $n(E) \propto \exp (-β(E-Φ_0))-1$, known as DARKexp. This differential energy distribution, and the resulting density structures, closely approximate those of dark-matter halos with central cusps, $ρ\sim r^{-1}$, and outer parts, $ρ\sim r^{-4}$. Here we define a non-equilibrium functional, $S_D$, which is maximized for DARKexp and increases monotonically during the evolution towards equilibrium of idealized collisionless systems of the Extended Spherical Infall Model. Systems that undergo more mixing more closely approach DARKexp. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Comments: 12 pages, 6 figures, 1 appendix, accepted to ApJ

An excursion into the core of the cluster lens Abell 1689

Authors: Agniva Ghosh, Dominic Adams, Liliya L. R. Williams, Jori Liesenborgs, Anahita Alavi, Claudia Scarlata

Abstract: Abell 1689 is a well studied cluster of galaxies and one of the largest gravitational lens systems ever observed. We have obtained a reconstruction of the cluster Abell 1689 using Grale, a free-form lens inversion method that relies exclusively on the multiple image data. Non-inclusion of any data related to cluster member galaxies ensures an unbiased measure of the mass distribution, which is the… ▽ More Abell 1689 is a well studied cluster of galaxies and one of the largest gravitational lens systems ever observed. We have obtained a reconstruction of the cluster Abell 1689 using Grale, a free-form lens inversion method that relies exclusively on the multiple image data. Non-inclusion of any data related to cluster member galaxies ensures an unbiased measure of the mass distribution, which is the most notable feature of free-form methods like Grale. We used two different sets of multiple image systems from the available strong lensing data - one containing only the secure systems (107 images), and the other containing all available systems, only excluding some very non-secure systems (151 images). For the very well-constrained central $\sim$100 kpc region of the cluster we made detailed comparison of the Grale reconstructed lensing mass and stellar mass retrieved by the Spectral Energy Distribution (SED) fitting software FAST++. We found a light-unaccompanied mass peak in this region, whose existence, while tentative, is favored by the distribution of nearby images that are local maxima in the Fermat potential. However, further tests, using different methodologies are needed to confirm the reality of this feature. If it shown to be real, this light-unaccompanied mass peak is consistent with dark matter self-interaction cross-section $σ\lesssim 1$cm$^2$/g, while being in tension with larger cross-sections. △ Less

Submitted 7 August, 2023; v1 submitted 17 June, 2022; originally announced June 2022.

Comments: 14 pages, 11 figures

Bridging the Gap Between Simply Parametrized and Free-Form Pixelated Models of Galaxy Lenses: The Case of WFI 2033-4723 Quad

Authors: Bernardo Barrera, Liliya L. R. Williams, Jonathan P. Coles, Philipp Denzel

Abstract: We study the radial and azimuthal mass distribution of the lensing galaxy in WFI2033-4723. Mindful of the fact that modeling results depend on modeling assumptions, we examine two very different recent models: simply parametrized (SP) models from the H0LiCOW collaboration, and pixelated free-form (FF) GLASS models. In addition, we fit our own models which are a compromise between the astrophysical… ▽ More We study the radial and azimuthal mass distribution of the lensing galaxy in WFI2033-4723. Mindful of the fact that modeling results depend on modeling assumptions, we examine two very different recent models: simply parametrized (SP) models from the H0LiCOW collaboration, and pixelated free-form (FF) GLASS models. In addition, we fit our own models which are a compromise between the astrophysical grounding of SP, and the flexibility of FF approaches. Our models consist of two offset parametric mass components, and generate many solutions, all fitting the quasar point image data. Among other results, we show that to reproduce point image properties the lensing mass must be lopsided, but the origin of this asymmetry can reside in the main lens plane or along the line of sight. We also show that there is a degeneracy between the slope of the density profile and the magnitude of external shear, and that the models from various modeling approaches are connected not by the mass sheet degeneracy, but by a more generalized transformation. Finally, we discuss interpretation degeneracy which afflicts all mass modeling: inability to correctly assign mass to the main lensing galaxy vs. nearby galaxies or line of sight structures. While this may not be a problem for the determination of $H_0$, interpretation degeneracy may become a major issue for the detailed study of galaxy structure. △ Less

Submitted 27 September, 2021; v1 submitted 9 August, 2021; originally announced August 2021.

Comments: 13 pages, 7 figures, submitted to OJA

Journal ref: The Open Journal of Astrophysics, Vol. 4, September 29, 2021

Further support for a trio of mass-to-light deviations in Abell 370: free-form Grale lens inversion using BUFFALO strong lensing data

Authors: Agniva Ghosh, Liliya L. R. Williams, Jori Liesenborgs, Ana Acebron, Mathilde Jauzac, Anton M. Koekemoer, Guillaume Mahler, Anna Niemiec, Charles Steinhardt, Andreas L. Faisst, David Lagattuta, Priyamvada Natarajan

Abstract: We use the Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) strong lensing image catalog of the merging galaxy cluster Abell 370 to obtain a mass model using the free-form lens inversion algorithm GRALE. The improvement of the strong lensing data quality results in a lens plane rms of only 0.45 arcsec, about a factor of two lower than that of our existing HFF v4 reconstruction.… ▽ More We use the Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) strong lensing image catalog of the merging galaxy cluster Abell 370 to obtain a mass model using the free-form lens inversion algorithm GRALE. The improvement of the strong lensing data quality results in a lens plane rms of only 0.45 arcsec, about a factor of two lower than that of our existing HFF v4 reconstruction. We attribute the improvement to spectroscopic data and use of the full reprocessed HST mosaics. In our reconstructed mass model, we found indications of three distinct mass features in Abell 370: (i) a $\sim 35$ kpc offset between the northern BCG and the nearest mass peak, (ii) a $\sim 100$ kpc mass concentration of roughly critical density $\sim 250$ kpc east of the main cluster, and (iii) a probable filament-like structure passing N-S through the cluster. While (i) is present in some form in most publicly available reconstructions spanning the range of modeling techniques: parametric, hybrid, and free-form, (ii) and (iii) are recovered by only about half of the reconstructions. We tested our hypothesis on the presence of the filament-like structure by creating a synthetic cluster - Irtysh IIIc - mocking the situation of a cluster with external mass. We also computed the source plane magnification distributions. Using them we estimated the probabilities of magnifications in the source plane, and scrutinized their redshift dependence. Finally, we explored the lensing effects of Abell 370 on the luminosity functions of sources at $z_s=9.0$, finding it consistent with published results. △ Less

Submitted 30 June, 2021; v1 submitted 23 April, 2021; originally announced April 2021.

Comments: 16 pages, 13 figures

Galaxy-lens determination of $H_0$: the effect of the ellipse+shear modeling assumption

Authors: Matthew R. Gomer, Liliya L. R. Williams

Abstract: Galaxy lenses are frequently modeled as an elliptical mass distribution with external shear and isothermal spheres to account for secondary and line-of-sight galaxies. There is statistical evidence that some fraction of observed quads are inconsistent with these assumptions, and require a dipole-like contribution to the mass with respect to the light. Simplifying assumptions about the shape of mas… ▽ More Galaxy lenses are frequently modeled as an elliptical mass distribution with external shear and isothermal spheres to account for secondary and line-of-sight galaxies. There is statistical evidence that some fraction of observed quads are inconsistent with these assumptions, and require a dipole-like contribution to the mass with respect to the light. Simplifying assumptions about the shape of mass distributions can lead to the incorrect recovery of parameters such as $H_0$. We create several tests of synthetic quad populations with different deviations from an elliptical shape, then fit them with an ellipse+shear model, and measure the recovered values of $H_0$. Kinematic constraints are not included. We perform two types of fittings -- one with a single point source and one with an array of sources emulating an extended source. We carry out two model-free comparisons between our mock quads and the observed population. One result of these comparisons is a statistical inconsistency not yet mentioned in the literature: the image distance ratios with respect to the lens center of observed quads appear to span a much wider range than those of synthetic or simulated quads. Bearing this discrepancy in mind, our mock populations can result in biases on $H_0$ $\sim10\%$. △ Less

Submitted 30 March, 2021; originally announced March 2021.

Comments: 15 pages, 5 figures; to be published in MNRAS

Exotic Image Formation in Strong Gravitational Lensing by Clusters of Galaxies -- II: Uncertainties

Authors: Ashish Kumar Meena, Agniva Ghosh, J. S. Bagla, L. L. R. Williams

Abstract: Due to the finite amount of observational data, the best-fit parameters corresponding to the reconstructed cluster mass have uncertainties. In turn, these uncertainties affect the inferences made from these mass models. Following our earlier work, we have studied the effect of such uncertainties on the singularity maps in simulated and actual galaxy clusters. The mass models for both simulated and… ▽ More Due to the finite amount of observational data, the best-fit parameters corresponding to the reconstructed cluster mass have uncertainties. In turn, these uncertainties affect the inferences made from these mass models. Following our earlier work, we have studied the effect of such uncertainties on the singularity maps in simulated and actual galaxy clusters. The mass models for both simulated and real clusters have been constructed using grale. The final best-fit mass models created using grale give the simplest singularity maps and a lower limit on the number of point singularities that a lens has to offer. The simple nature of these singularity maps also puts a lower limit on the number of three image (tangential and radial) arcs that a cluster lens has. Hence, we estimate the number of galaxy sources giving rise to the three image arcs, which can be observed with the James Webb Space Telescope (JWST). We find that we expect to observe at least 20-30 tangential and 5-10 radial three-image arcs in the Hubble Frontier Fields cluster lenses with the JWST. △ Less

Submitted 24 June, 2021; v1 submitted 25 March, 2021; originally announced March 2021.

Comments: 15 pages, 7 figures. Accepted in MNRAS

The Hubble constant from eight time-delay galaxy lenses

Authors: Philipp Denzel, Jonathan P. Coles, Prasenjit Saha, Liliya L. R. Williams

Abstract: We present a determination of the Hubble constant from the joint, free-form analysis of 8 strongly, quadruply lensing systems. In the concordance cosmology, we find $H_0 = 71.8^{+3.9}_{-3.3}\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{Mpc}^{-1}$ with a precision of $4.97\%$. This is in agreement with the latest measurements from Supernovae Type Ia and Planck observations of the cosmic microwave backgro… ▽ More We present a determination of the Hubble constant from the joint, free-form analysis of 8 strongly, quadruply lensing systems. In the concordance cosmology, we find $H_0 = 71.8^{+3.9}_{-3.3}\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{Mpc}^{-1}$ with a precision of $4.97\%$. This is in agreement with the latest measurements from Supernovae Type Ia and Planck observations of the cosmic microwave background. Our precision is lower compared to these and other recent time-delay cosmography determinations, because our modelling strategies reflect the systematic uncertainties of lensing degeneracies. We furthermore are able to find reasonable lensed image reconstructions by constraining to either value of $H_0$ from local and early Universe measurements. This leads us to conclude that current lensing constraints on $H_0$ are not strong enough to break the "Hubble tension" problem of cosmology. △ Less

Submitted 19 February, 2021; v1 submitted 28 July, 2020; originally announced July 2020.

Comments: 19 pages, 11 figures, published in MNRAS

Journal ref: MNRAS, Volume 501, Issue 1, February 2021, Pages 784-801

Two-component mass models of the lensing galaxy in the quadruply imaged supernova iPTF16geu

Authors: Liliya L. R. Williams, David Zegeye

Abstract: The first resolved, multiply imaged supernova Type Ia, iPTF16geu, was observed 4 years ago, five decades after such systems were first envisioned. Because of the unique properties of the source, these systems hold a lot of promise for the study of galaxy structure and cosmological parameters. However, this very first example presented modelers with a few puzzles. It was expected that to explain im… ▽ More The first resolved, multiply imaged supernova Type Ia, iPTF16geu, was observed 4 years ago, five decades after such systems were first envisioned. Because of the unique properties of the source, these systems hold a lot of promise for the study of galaxy structure and cosmological parameters. However, this very first example presented modelers with a few puzzles. It was expected that to explain image fluxes a contribution from microlensing by stars would be required, but to accommodate the magnitude of microlensing, the density slope of the elliptical power law lens model had to be quite shallow, $ρ_{2D} \propto r^{-0.7}$. Furthermore, the center of mass had to be displaced from that of observed light by ~0.1 kpc, and the position angle of light distribution was misaligned with that of mass by ~40 degrees. In this paper we present mass models that resolve the first two problems, and suggest a resolution of the third. Motivated by observations of local ellipticals, and some recent analysis of galaxy-scale lenses, our mass models consist of two offset (baryonic) mass components. The resulting mass distributions have a single centroid, but are lopsided, and have isodensity contours that are not purely elliptical and not self-similar with radius. For many of our models the microlensing requirements are modest, and the ring formed by the extended supernova host galaxy resembles the observed one. △ Less

Submitted 11 September, 2020; v1 submitted 16 June, 2020; originally announced June 2020.

Comments: 8 pages, 6 figures, accepted to The Open Journal of Astrophysics

Journal ref: The Open Journal of Astrophysics, Vol 3, September 15, 2020

Time Delay Lens Modelling Challenge

Authors: X. Ding, T. Treu, S. Birrer, G. C. -F. Chen, J. Coles, P. Denzel, M. Frigo A. Galan, P. J. Marshall, M. Millon, A. More, A. J. Shajib, D. Sluse, H. Tak, D. Xu, M. W. Auger, V. Bonvin, H. Chand, F. Courbin, G. Despali, C. D. Fassnacht, D. Gilman, S. Hilbert, S. R. Kumar, Y. -Y. Lin, J. W. Park , et al. (4 additional authors not shown)

Abstract: In recent years, breakthroughs in methods and data have enabled gravitational time delays to emerge as a very powerful tool to measure the Hubble constant $H_0$. However, published state-of-the-art analyses require of order 1 year of expert investigator time and up to a million hours of computing time per system. Furthermore, as precision improves, it is crucial to identify and mitigate systematic… ▽ More In recent years, breakthroughs in methods and data have enabled gravitational time delays to emerge as a very powerful tool to measure the Hubble constant $H_0$. However, published state-of-the-art analyses require of order 1 year of expert investigator time and up to a million hours of computing time per system. Furthermore, as precision improves, it is crucial to identify and mitigate systematic uncertainties. With this time delay lens modelling challenge we aim to assess the level of precision and accuracy of the modelling techniques that are currently fast enough to handle of order 50 lenses, via the blind analysis of simulated datasets. The results in Rung 1 and Rung 2 show that methods that use only the point source positions tend to have lower precision ($10 - 20\%$) while remaining accurate. In Rung 2, the methods that exploit the full information of the imaging and kinematic datasets can recover $H_0$ within the target accuracy ($ |A| < 2\%$) and precision ($< 6\%$ per system), even in the presence of poorly known point spread function and complex source morphology. A post-unblinding analysis of Rung 3 showed the numerical precision of the ray-traced cosmological simulations to be insufficient to test lens modelling methodology at the percent level, making the results difficult to interpret. A new challenge with improved simulations is needed to make further progress in the investigation of systematic uncertainties. For completeness, we present the Rung 3 results in an appendix, and use them to discuss various approaches to mitigating against similar subtle data generation effects in future blind challenges. △ Less

Submitted 17 February, 2021; v1 submitted 15 June, 2020; originally announced June 2020.

Comments: 23 pages, 12 figures, 6 tables, MNRAS accepted

Journal ref: MNRAS, (2021), 503, 1096

Free-form Grale lens inversion of galaxy clusters with up to 1000 multiple images

Authors: Agniva Ghosh, Liliya L. R. Williams, Jori Liesenborgs

Abstract: In the near future, ultra deep observations of galaxy clusters with HST or JWST will uncover $300-1000$ lensed multiple images, increasing the current count per cluster by up to an order of magnitude. This will further refine our view of clusters, leading to a more accurate and precise mapping of the total and dark matter distribution in clusters, and enabling a better understanding of background… ▽ More In the near future, ultra deep observations of galaxy clusters with HST or JWST will uncover $300-1000$ lensed multiple images, increasing the current count per cluster by up to an order of magnitude. This will further refine our view of clusters, leading to a more accurate and precise mapping of the total and dark matter distribution in clusters, and enabling a better understanding of background galaxy population and their luminosity functions. However, to effectively use that many images as input to lens inversion will require a re-evaluation of, and possibly upgrades to the existing methods. In this paper we scrutinize the performance of the free-form lens inversion method Grale in the regime of $150-1000$ input images, using synthetic massive galaxy clusters. Our results show that with an increasing number of input images, Grale produces improved reconstructed mass distributions, with the fraction of the lens plane recovered at better than $10\%$ accuracy increasing from $40-50\%$ for $\sim\!\!150$ images to $65\%$ for $\sim\!1000$ images. The reconstructed time delays imply a more precise measurement of $H_0$, with $\lesssim 1\%$ bias. While the fidelity of the reconstruction improves with the increasing number of multiple images used as model constraints, $\sim 150$ to $\sim 1000$, the lens plane rms deteriorates from $\sim 0.11''$ to $\sim 0.28''$. Since lens plane rms is not necessarily the best indicator of the quality of the mass reconstructions, looking for an alternative indicator is warranted. △ Less

Submitted 3 April, 2020; originally announced April 2020.

Comments: 18 pages, 12 figures, MNRAS, in press

Extended lens reconstructions with Grale: exploiting time domain, substructural and weak-lensing information

Authors: Jori Liesenborgs, Liliya L. R. Williams, Jenny Wagner, Sven De Rijcke

Abstract: The information about the mass density of galaxy clusters provided by the gravitational lens effect has inspired many inversion techniques. In this article, updates to the previously introduced method in Grale are described, and explored in a number of examples. The first looks into a different way of incorporating time delay information, not requiring the unknown source position. It is found that… ▽ More The information about the mass density of galaxy clusters provided by the gravitational lens effect has inspired many inversion techniques. In this article, updates to the previously introduced method in Grale are described, and explored in a number of examples. The first looks into a different way of incorporating time delay information, not requiring the unknown source position. It is found that this avoids a possible bias that leads to "over-focusing" the images, i.e. providing source position estimates that lie in a considerably smaller region than the true positions. The second is inspired by previous reconstructions of the cluster of galaxies MACS J1149.6+2223, where a multiply-imaged background galaxy contained a supernova, SN Refsdal, of which four additional images were produced by the presence of a smaller cluster galaxy. The inversion for the cluster as a whole, was not able to recover sufficient detail interior to this quad. We show how constraints on such different scales, from the entire cluster to a single member galaxy, can now be used, allowing such small scale substructures to be resolved. Finally, the addition of weak lensing information to this method is investigated. While this clearly helps recover the environment around the strong lensing region, the mass sheet degeneracy may make a full strong and weak inversion difficult, depending on the quality of the ellipticity information at hand. We encounter ring-like structure at the boundary of the two regimes, argued to be the result of combining strong and weak lensing constraints, possibly affected by degeneracies. △ Less

Submitted 23 March, 2020; originally announced March 2020.

Comments: 24 pages, 19 figures, accepted for publication in MNRAS

The BUFFALO HST Survey

Authors: Charles L. Steinhardt, Mathilde Jauzac, Ana Acebron, Hakim Atek, Peter Capak, Iary Davidzon, Dominique Eckert, David Harvey, Anton M. Koekemoer, Claudia D. P. Lagos, Guillaume Mahler, Mireia Montes, Anna Niemiec, Mario Nonino, P. A. Oesch, Johan Richard, Steven A. Rodney, Matthieu Schaller, Keren Sharon, Louis-Gregory Strolger, Joseph Allingham, Adam Amara, Yannick Bah'e, Celine Boehm, Sownak Bose , et al. (70 additional authors not shown)

Abstract: The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has no… ▽ More The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters. △ Less

Submitted 13 February, 2020; v1 submitted 27 January, 2020; originally announced January 2020.

Comments: Accepted ApJS; MAST archive will be live concurrent with publication

Image Flux Ratios of Gravitationally Lensed HS 0810+2554 with High Resolution Infrared Imaging

Authors: Terry Jay Jones, Liliya L. R. Williams, Steve Ertel, Philip M. Hinz, Amali Vaz, Shane Walsh, Ryan Webster

Abstract: We report near simultaneous imaging using LMIRCam on the LBTI of the quadruply imaged lensed quasar HS 0810+2554 at wavelengths of 2.16, 3.7 and $4.78~μ$m with a Full Width Half Max (FWHM) spatial resolution of $0^{\prime\prime}\!\!.13$, $0^{\prime\prime}\!\!.12$ and $0^{\prime\prime}\!\!.15$ respectively, comparable to HST optical imaging. In the $\rm{z} = 1.5$ rest frame of the quasar, the obser… ▽ More We report near simultaneous imaging using LMIRCam on the LBTI of the quadruply imaged lensed quasar HS 0810+2554 at wavelengths of 2.16, 3.7 and $4.78~μ$m with a Full Width Half Max (FWHM) spatial resolution of $0^{\prime\prime}\!\!.13$, $0^{\prime\prime}\!\!.12$ and $0^{\prime\prime}\!\!.15$ respectively, comparable to HST optical imaging. In the $\rm{z} = 1.5$ rest frame of the quasar, the observed wavelengths correspond to 0.86, 1.48, and $1.91~μ$m respectively. The two brightest images in the quad, A and B, are clearly resolved from each other with a separation of $0.187^{\prime\prime}$. The flux ratio of these two images (A/B) trends from 1.79 to 1.23 from 2.16 to $4.78~μ$m. The trend in flux ratio is consistent with the $2.16~μ$m flux originating from a small sized accretion disk in the quasar that experiences only microlensing. The excess flux above the contribution from the accretion disk at the two longer wavelengths originates from a larger sized region that experiences no microlensing. A simple model employing multiplicative factors for image B due to stellar microlensing $(m)$ and sub-structure millilensing $(M)$ is presented. The result is tightly constrained to the product $m\times M=1.79$. Given the observational errors, the 60\% probability contour for this product stretches from $m= 2.6$, $M = 0.69$ to $m= 1.79$, $M = 1.0$, where the later is consistent with microlensing only. △ Less

Submitted 23 October, 2019; originally announced October 2019.

Comments: accepted AJ

Model-independent and model-based local lensing properties of B0128+437 from resolved quasar images

Authors: Jenny Wagner, Liliya L. R. Williams

Abstract: The galaxy-scale gravitational lens B0128+437 generates a quadrupole-image configuration of a background quasar that shows milli-arcsecond-scale subcomponents in the multiple images observed with VLBI. As this multiple-image configuration including the subcomponents has eluded a parametric lens-model characterisation so far, we determine local lens properties at the positions of the multiple image… ▽ More The galaxy-scale gravitational lens B0128+437 generates a quadrupole-image configuration of a background quasar that shows milli-arcsecond-scale subcomponents in the multiple images observed with VLBI. As this multiple-image configuration including the subcomponents has eluded a parametric lens-model characterisation so far, we determine local lens properties at the positions of the multiple images with our model-independent approach. Using PixeLens, we also succeed in setting up a global free-form mass density reconstruction including all subcomponents as constraints. We compare the model-independent local lens properties with those obtained by PixeLens and those obtained by the parametric modelling algorithm Lensmodel. A comparison of all three approaches and a model-free analysis based on the relative polar angles of the multiple images corroborate the hypothesis that elliptically symmetric models are too simplistic to characterise the asymmetric mass density distribution of this lenticular or late-type galaxy. In addition, the model-independent approach efficiently determines local lens properties on the scale of the quasar subcomponents, which are computationally intensive to obtain by free-form model-based approaches. As only 40% of the small-scale subcomponent local lens properties overlap within the 1-$σ$ confidence bounds, mass density gradients on milli-arcsecond scales cannot be excluded. Hence, aiming at a global reconstruction of the deflecting mass density distribution, increasingly detailed observations require flexible free-form models that allow for density fluctuations on milli-arcsecond scale to replace parametric ones, especially for asymmetric lenses or lenses with localised inhomogeneities like B0128. △ Less

Submitted 3 September, 2019; originally announced September 2019.

Comments: 15 pages, submitted to Astronomy & Astrophysics, comments welcome

Journal ref: A&A 635, A86 (2020)

Free-form Grale reconstruction of Abell 2744: robustness of uncertainties against changes in lensing data

Authors: Kevin Sebesta, Liliya L. R. Williams, Jori Liesenborgs, Elinor Medezinski, Nobuhiro Okabe

Abstract: Abell 2744, a massive Hubble Frontier Fields merging galaxy cluster with many multiple images in the core has been the subject of many lens inversions using different methods. While most existing studies compare various inversion methods, we focus on a comparison of reconstructions that use different input lensing data. Since the quantity and quality of lensing data is constantly improving, it mak… ▽ More Abell 2744, a massive Hubble Frontier Fields merging galaxy cluster with many multiple images in the core has been the subject of many lens inversions using different methods. While most existing studies compare various inversion methods, we focus on a comparison of reconstructions that use different input lensing data. Since the quantity and quality of lensing data is constantly improving, it makes sense to ask if the estimated uncertainties are robust against changes in the data. We address this question using free-form Grale, which takes only image information as input, and nothing pertaining to cluster galaxies. We reconstruct Abell 2744 using two sets of strong lensing data from the Hubble Frontier Fields community. Our first and second reconstructions use 55 and 91 images, respectively, and only 10 of the 91 images have the same positions and redshifts as in the first reconstruction. Comparison of the two mass maps shows that Grale uncertainties are robust against these changes, as well as small modifications in the inversion routine. Additionally, applying the methods used in Sebesta et al. (2016) for MACS J0416, we conclude that, in a statistical sense, light follows mass in Abell 2744, with brighter galaxies clustering stronger with the recovered mass than the fainter ones. We also show that the faintest galaxies are anti-correlated with mass, which is likely the result of light contamination from bright galaxies, and lensing magnification bias acting on galaxies background to the cluster. △ Less

Submitted 22 July, 2019; originally announced July 2019.

Comments: 12 pages, 7 figure, MNRAS, in press

Galaxy-lens determination of $H_0$: constraining density slope in the context of the mass sheet degeneracy

Authors: Matthew R. Gomer, Liliya L. R. Williams

Abstract: Gravitational lensing offers a competitive method to measure $H_0$ with the goal of 1% precision. A major obstacle comes in the form of lensing degeneracies, such as the mass sheet degeneracy (MSD), which make it possible for a family of density profiles to reproduce the same lensing observables but return different values of $H_0$. The modeling process artificially selects one choice from this fa… ▽ More Gravitational lensing offers a competitive method to measure $H_0$ with the goal of 1% precision. A major obstacle comes in the form of lensing degeneracies, such as the mass sheet degeneracy (MSD), which make it possible for a family of density profiles to reproduce the same lensing observables but return different values of $H_0$. The modeling process artificially selects one choice from this family, potentially biasing $H_0$. The effect is more pronounced when the profile of a given lens is not perfectly described by the lens model, which will always be the case to some extent. To explore this, we quantify the bias and spread in $H_0$ by creating quads from two-component mass models and fitting them with a power-law ellipse+shear model. We find that the bias does not correspond to the estimate one would calculate by transforming the profile into a power law near the image radius. We also emulate the effect of including stellar kinematics by performing fits where the slope is constrained to the true value. Informing the fit using the true value near the image radius can introduce substantial bias (0-23% depending on the model). We confirm using Jeans arguments that kinematic constraints can result in a biased value of $H_0$, though the degree of bias depends on the region kinematic modeling probes in specific lenses. We conclude that lensing degeneracies manifest through commonplace modeling approaches in a more complicated way than is assumed in the literature. If stellar kinematics incorrectly break the MSD, their inclusion may introduce more bias than their omission. △ Less

Submitted 3 November, 2020; v1 submitted 19 July, 2019; originally announced July 2019.

Comments: 37 pages, 11 figures, accepted by JCAP

Journal ref: JCAP 11 (2020) 045

Power law pseudo phase-space density profiles of dark matter halos: fluke of physics?

Authors: Arpit Arora, Liliya L. R. Williams

Abstract: It has been known for nearly 20 years that the pseudo phase-space density profile of equilibrium simulated dark matter halos, $ρ(r)/σ^3(r)$, is well described by a power law over 3 decades in radius, even though both the density $ρ(r)$, and the velocity dispersion $σ(r)$ deviate significantly from power laws. The origin of this scale-free behavior is not understood. It could be an inherent propert… ▽ More It has been known for nearly 20 years that the pseudo phase-space density profile of equilibrium simulated dark matter halos, $ρ(r)/σ^3(r)$, is well described by a power law over 3 decades in radius, even though both the density $ρ(r)$, and the velocity dispersion $σ(r)$ deviate significantly from power laws. The origin of this scale-free behavior is not understood. It could be an inherent property of self-gravitating collisionless systems, or it could be a mere coincidence. To address the question we work with equilibrium halos, and more specifically, the second derivative of the Jeans equation, which, under the assumptions of (i) Einasto density profile, (ii) linear velocity anisotropy - density slope relation, and (iii) $ρ/σ^3\propto r^{-α}$, can be transformed from a differential equation to a cubic algebraic equation. Relations (i)-(iii) are all observed in numerical simulations, and are well parametrized by a total of 4 or 6 model parameters. We do not consider dynamical evolution of halos; instead, taking advantage of the fact that the algebraic Jeans equation for equilibrium halos puts relations (i)-(iii) on the same footing, we study the (approximate) solutions of this equation in the 4 and 6 dimensional spaces. We argue that the distribution of best solutions in these parameter spaces is inconsistent with $ρ/σ^3\propto r^{-α}$ being an fundamental property of gravitational evolution, and conclude that the scale-free nature of this quantity is likely to be a fluke. △ Less

Submitted 12 March, 2020; v1 submitted 7 April, 2019; originally announced April 2019.

Comments: ApJ, Accepted

Journal ref: ApJ 893 53 (2020)

Searching for Highly Magnified Stars at Cosmological Distances: Discovery of a Redshift 0.94 Blue Supergiant in Archival Images of the Galaxy Cluster MACS J0416.1-2403

Authors: Wenlei Chen, Patrick L. Kelly, Jose M. Diego, Masamune Oguri, Liliya L. R. Williams, Adi Zitrin, Tommaso L. Treu, Nathan Smith, Thomas J. Broadhurst, Nick Kaiser, Ryan J. Foley, Alexei V. Filippenko, Laura Salo, Jens Hjorth, Jonatan Selsing

Abstract: Individual highly magnified stars have been recently discovered at lookback times of more than half the age of the Universe, in lensed galaxies that straddle the critical curves of massive galaxy clusters. Having confirmed their detectability, it is now important to carry out systematic searches for them in order to establish their frequency, and in turn learn about the statistical properties of h… ▽ More Individual highly magnified stars have been recently discovered at lookback times of more than half the age of the Universe, in lensed galaxies that straddle the critical curves of massive galaxy clusters. Having confirmed their detectability, it is now important to carry out systematic searches for them in order to establish their frequency, and in turn learn about the statistical properties of high-redshift stars and of the granularity of matter in the foreground deflector. Here we report the discovery of a highly magnified star at redshift $z = 0.94$ in a strongly lensed arc behind a Hubble Frontier Field galaxy cluster, MACS J0416.1-2403, discovered as part of a systematic archival search. The bright transient (dubbed "Warhol") was discovered in Hubble Space Telescope data taken on 2014 September 15 and 16. This single image faded over a period of two weeks, and observations taken on 2014 September 1 show that the duration of the microlensing event was at most four weeks in total. The light curve may also exhibit slow changes over a period of years consistent with the level of microlensing expected from stars responsible for the intracluster light (ICL) of the cluster. Optical and infrared observations taken near peak brightness can be fit by a stellar spectrum with moderate host-galaxy extinction. A blue supergiant matches the measured spectral energy distribution near peak, implying a temporary magnification of at least several thousand. While the spectrum of an O-type star would also fit the transient's spectral energy distribution, extremely luminous O-type stars are much less common than blue supergiants. The short timescale of the event and the estimated effective temperature indicate that the lensed source is an extremely magnified star. △ Less

Submitted 14 February, 2019; originally announced February 2019.

Testing the uniqueness of gravitational lens mass models

Authors: Levi G. Walls, Liliya L. R. Williams

Abstract: The positions of images produced by the gravitational lensing of background sources provide unique insight in to galaxy-lens mass distribution. However, even quad images of extended sources are not able to fully characterize the central regions of the host galaxy. Most previous work has focused either on the radial density profile of the lenses or localized substructure clumps. Here, we concentrat… ▽ More The positions of images produced by the gravitational lensing of background sources provide unique insight in to galaxy-lens mass distribution. However, even quad images of extended sources are not able to fully characterize the central regions of the host galaxy. Most previous work has focused either on the radial density profile of the lenses or localized substructure clumps. Here, we concentrate on the azimuthal mass asymmetries near the image circle. The motivation for considering such mass inhomogeneities is that the transition between the central stellar dominated region and the outer dark matter dominated region, though well represented by a power law density profile, is unlikely to be featureless, and encodes information about the dynamical state and assembly history of galaxies. It also happens to roughly coincide with the Einstein radius. We ask if galaxies that have mass asymmetries beyond ellipticity can be modeled with simpler lenses, i.e., can complex mass distributions masquerade as simple elliptical+shear lenses? Our preliminary study indicates that for galaxies with elliptical stellar and dark matter distributions, but with no mass asymmetry, and an extended source filling the diamond caustic, an elliptical+shear lens model can reproduce the images well, thereby hiding the potential complexity of the actual mass distribution. For galaxies with non-zero mass asymmetry, the answer depends on the size and brightness distribution of the source, and its location within the diamond caustic. In roughly half of the cases we considered the mass asymmetries can easily evade detection. △ Less

Submitted 20 August, 2018; v1 submitted 16 August, 2018; originally announced August 2018.

Comments: 12 pages, 9 figures, MNRAS, accepted

The role of multiple images and model priors in measuring $H_0$ from supernova Refsdal in galaxy cluster MACS J1149.5+2223

Authors: Liliya L. R. Williams, Jori Liesenborgs

Abstract: Multiple image gravitational lensing systems with measured time delays provide a promising one-step method for determining $H_0$. MACS J1149, which lenses SN Refsdal into a quad S1-S4, and two other widely separated images, SX and SY, is a perfect candidate. If time delays are pinned down, the remaining uncertainty arises from the mass distribution in the lens. In MACS J1149, the mass in the relev… ▽ More Multiple image gravitational lensing systems with measured time delays provide a promising one-step method for determining $H_0$. MACS J1149, which lenses SN Refsdal into a quad S1-S4, and two other widely separated images, SX and SY, is a perfect candidate. If time delays are pinned down, the remaining uncertainty arises from the mass distribution in the lens. In MACS J1149, the mass in the relevant lens plane region can be constrained by (i) many multiple images, (ii) the mass of the galaxy splitting S1-S4 (which, we show, is correlated with $H_0$), (iii) magnification of SX (also correlated with $H_0$), and (iv) prior assumptions on the mass distribution. Our goal is not to estimate $H_0$, but to understand its error budget, i.e., estimate uncertainties associated with each of these constraints. Using multiple image positions alone, yields very large uncertainty, despite the fact that the position of SX is recovered to within $\!\le\!0.036$" (rms $\!\le\!0.36$") by GRALE lens inversion. Fixing the mass of the galaxy that splits S1-S4 reduces $1σ$ uncertainties to $\sim 23\%$, while fixing the magnification of SX yields $1σ$ uncertainties of $32\%$. We conclude that smaller uncertainties, of order few percent, are a consequence of imposing prior assumptions on the shapes of the galaxy and cluster mass distributions, which may or may not apply in a highly non-equilibrium environment of a merging cluster. We propose that if a measurement of $H_0$ is to be considered reliable, it must be supported by a wide range of lens inversion methods. △ Less

Submitted 8 November, 2018; v1 submitted 28 June, 2018; originally announced June 2018.

Comments: 14 pages, 11 figures, accepted to MNRAS

The impact of $Λ$CDM substructure and baryon-dark matter transition on the image positions of quad galaxy lenses

Authors: Matthew R. Gomer, Liliya L. R. Williams

Abstract: The positions of multiple images in galaxy lenses are related to the galaxy mass distribution. Smooth elliptical mass profiles were previously shown to be inadequate in reproducing the quad population. In this paper, we explore the deviations from such smooth elliptical mass distributions. Unlike most other work, we use a model-free approach based on the relative polar image angles of quads, and t… ▽ More The positions of multiple images in galaxy lenses are related to the galaxy mass distribution. Smooth elliptical mass profiles were previously shown to be inadequate in reproducing the quad population. In this paper, we explore the deviations from such smooth elliptical mass distributions. Unlike most other work, we use a model-free approach based on the relative polar image angles of quads, and their position in 3D space with respect to the Fundamental Surface of Quads. The FSQ is defined by quads produced by elliptical lenses. We have generated thousands of quads from synthetic populations of lenses with substructure consistent with $Λ$CDM simulations, and found that such perturbations are not sufficient to match the observed distribution of quads relative to the FSQ. The result is unchanged even when subhalo masses are increased by a factor of ten, and the most optimistic lensing selection bias is applied. We then produce quads from galaxies created using two components, representing baryons and dark matter. The transition from the mass being dominated by baryons in inner radii to being dominated by dark matter in outer radii can carry with it asymmetries, which would affect relative image angles. We run preliminary experiments using lenses with two elliptical mass components with nonidentical axis ratios and position angles, perturbations from ellipticity in the form of nonzero Fourier coefficients $a_4$ and $a_6$, and artificially offset ellipse centers as a proxy for asymmetry at image radii. We show that combination of these effects is a promising way of accounting for quad population properties. We conclude that the quad population provides a unique and sensitive tool for constraining detailed mass distribution in the centers of galaxies. △ Less

Submitted 21 February, 2018; v1 submitted 20 December, 2017; originally announced December 2017.

Comments: 18 pages, 15 figures, 2 tables

Journal ref: MNRAS 475 (2018) 1987

Evidence for the line of sight structure in the largest Hubble Frontier Field cluster, MACSJ0717.5+3745

Authors: Liliya L. R. Williams, Kevin Sebesta, Jori Liesenborgs

Abstract: MACS J0717 is the most massive and extended of the Hubble Frontier Field clusters. It is one of the more difficult clusters to model, and we argue that this is in part due to the line of sight structure (LoS) at redshifts beyond 2. We show that the Grale mass reconstruction based on sources at 3<z_s<4.1 has at least 10^{13}M_sun more mass than that based on nearby sources, z_s<2.6, and attribute t… ▽ More MACS J0717 is the most massive and extended of the Hubble Frontier Field clusters. It is one of the more difficult clusters to model, and we argue that this is in part due to the line of sight structure (LoS) at redshifts beyond 2. We show that the Grale mass reconstruction based on sources at 3<z_s<4.1 has at least 10^{13}M_sun more mass than that based on nearby sources, z_s<2.6, and attribute the excess mass to a putative LoS, which is at least 75" from the cluster center. Furthermore, the lens-model fitted z_s's of the recent Kawamata et al. reconstruction are biased systematically low compared to photometric z_s's, and the bias is a function of images' distance from the cluster center. We argue that these mimic the effect of LoS. We conclude that even in the presence of 100-200 images, lens-model adjusted source redshifts can conceal the presence of LoS, demonstrating the existence of degeneracies between z_s and (sub)structure. Also, a very good fit to image positions is not a sufficient condition for having a high fidelity mass map: Kawamata et al. obtain an rms of 0.52" for 173 images of 60 sources; our Grale reconstruction of the exact same data yields a somewhat different map, but similarly low rms, 0.62". In contrast, a Grale model that uses reasonable, but fixed z_s gives a worse rms of 1.28" for 44 sources with 126 images. Unaccounted for LoS can bias the mass map, affecting the magnification and luminosity function estimates of high redshift sources. △ Less

Submitted 1 August, 2018; v1 submitted 14 November, 2017; originally announced November 2017.

Comments: 14 pages, 8 figures, MNRAS, in press

Dynamics of merging: Post-merger mixing and relaxation of an Illustris galaxy

Authors: Anthony M. Young, Liliya L. R. Williams, Jens Hjorth

Abstract: During the merger of two galaxies, the resulting system undergoes violent relaxation and seeks stable equilibrium. However, the details of this evolution are not fully understood. Using Illustris simulation, we probe two physically related processes, mixing and relaxation. Though the two are driven by the same dynamics---global time-varying potential for the energy, and torques caused by asymmetri… ▽ More During the merger of two galaxies, the resulting system undergoes violent relaxation and seeks stable equilibrium. However, the details of this evolution are not fully understood. Using Illustris simulation, we probe two physically related processes, mixing and relaxation. Though the two are driven by the same dynamics---global time-varying potential for the energy, and torques caused by asymmetries for angular momentum---we measure them differently. We define mixing as the redistribution of energy and angular momentum between particles of the two merging galaxies. We assess the degree of mixing as the difference between the shapes of their N(E)s, and their N(L^2)s. We find that the difference is decreasing with time, indicating mixing. To measure relaxation, we compare N(E) of the newly merged system to N(E) of a theoretical prediction for relaxed collisionless systems, DARKexp, and witness the system becoming more relaxed, in the sense that N(E) approaches DARKexp N(E). Because the dynamics driving mixing and relaxation are the same, the timescale is similar for both. We measure two sequential timescales: a rapid, 1 Gyr phase after the initial merger, during which the difference in N(E) of the two merging halos decreases by ~80%, followed by a slow phase, when the difference decreases by ~50% over ~8.5 Gyrs. This is a direct measurement of the relaxation timescale. Our work also draws attention to the fact that when a galaxy has reached Jeans equilibrium it may not yet have reached a fully relaxed state given by DARKexp, in that it retains information about its past history. This manifests itself most strongly in stars being centrally concentrated. We argue that it is particularly difficult for stars, and other tightly bound particles, to mix because they have less time to be influenced by the fluctuating potential, even across multiple merger events. △ Less

Submitted 14 February, 2018; v1 submitted 31 August, 2017; originally announced September 2017.

Comments: accepted for publication in JCAP

Two Peculiar Fast Transients in a Strongly Lensed Host Galaxy

Authors: S. A. Rodney, I. Balestra, M. Bradac, G. Brammer, T. Broadhurst, G. B. Caminha, G. Chirivi, J. M. Diego, A. V. Filippenko, R. J. Foley, O. Graur, C. Grillo, S. Hemmati, J. Hjorth, A. Hoag, M. Jauzac, S. W. Jha, R. Kawamata, P. L. Kelly, C. McCully, B. Mobasher, A. Molino, M. Oguri, J. Richard, A. G. Riess , et al. (10 additional authors not shown)

Abstract: A massive galaxy cluster can serve as a magnifying glass for distant stellar populations, with strong gravitational lensing exposing details in the lensed background galaxies that would otherwise be undetectable. The MACS J0416.1-2403 cluster (hereafter MACS0416) is one of the most efficient lenses in the sky, and in 2014 it was observed with high-cadence imaging from the Hubble Space Telescope (H… ▽ More A massive galaxy cluster can serve as a magnifying glass for distant stellar populations, with strong gravitational lensing exposing details in the lensed background galaxies that would otherwise be undetectable. The MACS J0416.1-2403 cluster (hereafter MACS0416) is one of the most efficient lenses in the sky, and in 2014 it was observed with high-cadence imaging from the Hubble Space Telescope (HST). Here we describe two unusual transient events that appeared behind MACS0416 in a strongly lensed galaxy at redshift $z = 1.0054 \pm 0.0002$. These transients---designated HFF14Spo-NW and HFF14Spo-SE and collectively nicknamed "Spock"---were faster and fainter than any supernova (SN), but significantly more luminous than a classical nova. They reached peak luminosities of $\sim10^{41}$ erg s$^{-1}$ ($M_{\rm AB} < -14$ mag) in 5 rest-frame days, then faded below detectability in roughly the same time span. Models of the cluster lens suggest that these events may be spatially coincident at the source plane, but are most likely not temporally coincident. We find that HFF14Spo can be explained as a luminous blue variable (LBV), a recurrent nova (RN), or a pair of stellar microlensing events. To distinguish between these hypotheses will require a clarification of the positions of nearby critical curves, along with high-cadence monitoring of the field that could detect new transient episodes in the host galaxy. △ Less

Submitted 8 July, 2017; originally announced July 2017.

Comments: 14 pages (Main text + figures) + 40 pages (Methods, References and Supplementary Info)

The Frontier Fields Lens Modeling Comparison Project

Authors: M. Meneghetti, P. Natarajan, D. Coe, E. Contini, G. De Lucia, C. Giocoli, A. Acebron, S. Borgani, M. Bradac, J. M. Diego, A. Hoag, M. Ishigaki, T. L. Johnson, E. Jullo, R. Kawamata, D. Lam, M. Limousin, J. Liesenborgs, M. Oguri, K. Sebesta, K. Sharon, L. L. R. Williams, A. Zitrin

Abstract: Gravitational lensing by clusters of galaxies offers a powerful probe of their structure and mass distribution. Deriving a lens magnification map for a galaxy cluster is a classic inversion problem and many methods have been developed over the past two decades to solve it. Several research groups have developed techniques independently to map the predominantly dark matter distribution in cluster l… ▽ More Gravitational lensing by clusters of galaxies offers a powerful probe of their structure and mass distribution. Deriving a lens magnification map for a galaxy cluster is a classic inversion problem and many methods have been developed over the past two decades to solve it. Several research groups have developed techniques independently to map the predominantly dark matter distribution in cluster lenses. While these methods have all provided remarkably high precision mass maps, particularly with exquisite imaging data from the Hubble Space Telescope (HST), the reconstructions themselves have never been directly compared. In this paper, we report the results of comparing various independent lens modeling techniques employed by individual research groups in the community. Here we present for the first time a detailed and robust comparison of methodologies for fidelity, accuracy and precision. For this collaborative exercise, the lens modeling community was provided simulated cluster images -- of two clusters Ares and Hera -- that mimic the depth and resolution of the ongoing HST Frontier Fields. The results of the submitted reconstructions with the un-blinded true mass profile of these two clusters are presented here. Parametric, free-form and hybrid techniques have been deployed by the participating groups and we detail the strengths and trade-offs in accuracy and systematics that arise for each methodology. We note in conclusion that lensing reconstruction methods produce reliable mass distributions that enable the use of clusters as extremely valuable astrophysical laboratories and cosmological probes. △ Less

Submitted 14 June, 2016; originally announced June 2016.

Comments: 38 pages, 25 figures, submitted to MNRAS, version with full resolution images can be found at http://pico.bo.astro.it/~massimo/papers/FFsims.pdf

Lens Models Under the Microscope: Comparison of Hubble Frontier Field Cluster Magnification Maps

Authors: Jett Priewe, Liliya L. R. Williams, Jori Liesenborgs, Dan Coe, Steven A. Rodney

Abstract: Using the power of gravitational lensing magnification by massive galaxy clusters, the Hubble Frontier Fields provide deep views of six patches of the high redshift Universe. The combination of deep Hubble imaging and exceptional lensing strength has revealed the greatest numbers of multiply-imaged galaxies available to constrain models of cluster mass distributions. However, even with… ▽ More Using the power of gravitational lensing magnification by massive galaxy clusters, the Hubble Frontier Fields provide deep views of six patches of the high redshift Universe. The combination of deep Hubble imaging and exceptional lensing strength has revealed the greatest numbers of multiply-imaged galaxies available to constrain models of cluster mass distributions. However, even with $\mathcal{O}(100)$ images per cluster, the uncertainties associated with the reconstructions are not negligible. The goal of this paper is to show the diversity of model magnification predictions. We examine 7 and 9 mass models of Abell 2744 and MACS J0416, respectively, submitted to the Mikulski Archive for Space Telescopes for public distribution in September 2015. The dispersion between model predictions increases from 30% at common low magnifications ($μ\sim2$) to 70% at rare high magnifications ($μ\sim40$). MACS J0416 exhibits smaller dispersions than Abell 2744 for $2<μ<10$. We show that magnification maps based on different lens inversion techniques typically differ from each other by more than their quoted statistical errors. This suggests that some models underestimate the true uncertainties, which are primarily due to various lensing degeneracies. Though the exact mass sheet degeneracy is broken, its generalized counterpart is not broken at least in Abell 2744. Other, local degeneracies are also present in both clusters. Our comparison of models is complementary to the comparison of reconstructions of known synthetic mass distributions. By focusing on observed clusters, we can identify those that are best constrained, and therefore provide the clearest view of the distant Universe. △ Less

Submitted 28 October, 2016; v1 submitted 24 May, 2016; originally announced May 2016.

Comments: 17 pages, 18 figures, MNRAS, in press

Ubiquity of density slope oscillations in the central regions of galaxy and cluster-sized systems

Authors: Anthony M. Young, Liliya L. R. Williams, Jens Hjorth

Abstract: One usually thinks of a radial density profile as having a monotonically changing logarithmic slope, such as in NFW or Einasto profiles. However, in two different classes of commonly used systems, this is often not the case. These classes exhibit non-monotonic changes in their density profile slopes which we call oscillations for short. We analyze these two unrelated classes separately. Class 1 co… ▽ More One usually thinks of a radial density profile as having a monotonically changing logarithmic slope, such as in NFW or Einasto profiles. However, in two different classes of commonly used systems, this is often not the case. These classes exhibit non-monotonic changes in their density profile slopes which we call oscillations for short. We analyze these two unrelated classes separately. Class 1 consists of systems that have density oscillations and that are defined through their distribution function $f(E)$, or differential energy distribution $N(E)$, such as isothermal spheres, King profiles, or DARKexp, a theoretically derived model for relaxed collisionless systems. Systems defined through $f(E)$ or $N(E)$ generally have density slope oscillations. Class 1 system oscillations can be found at small, intermediate, or large radii but we focus on a limited set of Class 1 systems that have oscillations in the central regions, usually at $\log(r/r_{-2})\lesssim -2$, where $r_{-2}$ is the largest radius where $d\log(ρ)/d\log(r)=-2$. We show that the shape of their $N(E)$ can roughly predict the amplitude of oscillations. Class 2 systems which are a product of dynamical evolution, consist of observed and simulated galaxies and clusters, and pure dark matter halos. Oscillations in the density profile slope seem pervasive in the central regions of Class 2 systems. We argue that in these systems, slope oscillations are an indication that a system is not fully relaxed. We show that these oscillations can be reproduced by small modifications to $N(E)$ of DARKexp. These affect a small fraction of systems' mass and are confined to $\log(r/r_{-2})\lesssim 0$. The size of these modifications serves as a potential diagnostic for quantifying how far a system is from being relaxed. △ Less

Submitted 5 April, 2016; originally announced April 2016.

Comments: accepted by the Journal of Cosmology and Astroparticle Physics (JCAP)

Testing DARKexp against energy and density distributions of Millennium-II halos

Authors: Chris Nolting, Liliya L. R. Williams, Michael Boylan-Kolchin, Jens Hjorth

Abstract: We test the DARKexp model for relaxed, self-gravitating, collisionless systems against equilibrium dark matter halos from the Millennium-II simulation. While limited tests of DARKexp against simulations and observations have been carried out elsewhere, this is the first time the testing is done with a large sample of simulated halos spanning a factor of ~ 50 in mass, and using independent fits to… ▽ More We test the DARKexp model for relaxed, self-gravitating, collisionless systems against equilibrium dark matter halos from the Millennium-II simulation. While limited tests of DARKexp against simulations and observations have been carried out elsewhere, this is the first time the testing is done with a large sample of simulated halos spanning a factor of ~ 50 in mass, and using independent fits to density and energy distributions. We show that DARKexp, a one shape parameter family, provides very good fits to the shapes of density profiles, ρ(r), and differential energy distributions, N(E), of individual simulated halos. The best fit shape parameter $φ_{0}$ obtained from the two types of fits are correlated, though with scatter. Our most important conclusions come from ρ(r) and N(E) that have been averaged over many halos. These show that the bulk of the deviations between DARKexp and individual Millennium-II halos come from halo-to-halo fluctuations, likely driven by substructure, and other density perturbations. The average ρ(r) and N(E) are quite smooth and follow DARKexp very closely. The only deviation that remains after averaging is small, and located at most bound energies for N(E) and smallest radii for ρ(r). Since the deviation is confined to 3-4 smoothing lengths, and is larger for low mass halos, it is likely due to numerical resolution effects. △ Less

Submitted 6 September, 2016; v1 submitted 23 March, 2016; originally announced March 2016.

Comments: 24 pages, 12 figures, Accepted for publication in the Journal of Cosmology and Astroparticle Physics (JCAP)

Journal ref: JCAP09(2016)042

Non-universality of dark-matter halos: cusps, cores, and the central potential

Authors: Jens Hjorth, Liliya L. R. Williams, Radosław Wojtak, Michael McLaughlin

Abstract: Dark-matter halos grown in cosmological simulations appear to have central NFW-like density cusps with mean values of $d\logρ/d\log r \approx -1$, and some dispersion, which is generally parametrized by the varying index $α$ in the Einasto density profile fitting function. Non-universality in profile shapes is also seen in observed galaxy clusters and possibly dwarf galaxies. Here we show that non… ▽ More Dark-matter halos grown in cosmological simulations appear to have central NFW-like density cusps with mean values of $d\logρ/d\log r \approx -1$, and some dispersion, which is generally parametrized by the varying index $α$ in the Einasto density profile fitting function. Non-universality in profile shapes is also seen in observed galaxy clusters and possibly dwarf galaxies. Here we show that non-universality, at any given mass scale, is an intrinsic property of DARKexp, a theoretically derived model for collisionless self-gravitating systems. We demonstrate that DARKexp - which has only one shape parameter, $φ_0$ - fits the dispersion in profile shapes of massive simulated halos as well as observed clusters very well. DARKexp also allows for cored dark-matter profiles, such as those found for dwarf spheroidal galaxies. We provide approximate analytical relations between DARKexp $φ_0$, Einasto $α$, or the central logarithmic slope in the Dehnen-Tremaine analytical $γ$-models. The range in halo parameters reflects a substantial variation in the binding energies per unit mass of dark-matter halos. △ Less

Submitted 10 August, 2015; originally announced August 2015.

Comments: ApJ, in press, 10 pages, 7 figures

Journal ref: ApJ 811, 2 (2015)

Testing light-traces-mass in Hubble Frontier Fields Cluster MACS-J0416.1-2403

Authors: Kevin Sebesta, Liliya L. R. Williams, Irshad Mohammed, Prasenjit Saha, Jori Liesenborgs

Abstract: We reconstruct the projected mass distribution of a massive merging Hubble Frontier Fields cluster MACSJ0416 using the genetic algorithm based free-form technique called Grale. The reconstructions are constrained by 149 lensed images identified by Jauzac et al. using HFF data. No information about cluster galaxies or light is used, which makes our reconstruction unique in this regard. Using visual… ▽ More We reconstruct the projected mass distribution of a massive merging Hubble Frontier Fields cluster MACSJ0416 using the genetic algorithm based free-form technique called Grale. The reconstructions are constrained by 149 lensed images identified by Jauzac et al. using HFF data. No information about cluster galaxies or light is used, which makes our reconstruction unique in this regard. Using visual inspection of the maps, as well as galaxy-mass correlation functions we conclude that overall light does follow mass. Furthermore, the fact that brighter galaxies are more strongly clustered with mass is an important confirmation of the standard biasing scenario in galaxy clusters. On the smallest scales, approximately less than a few arcseconds, the resolution afforded by 149 images is still not sufficient to confirm or rule out galaxy-mass offsets of the kind observed in ACO 3827. We also compare the mass maps of MACSJ0416 obtained by three different groups: Grale, and two parametric Lenstool reconstructions from the CATS and Sharon/Johnson teams. Overall, the three agree well; one interesting discrepancy between Grale and Lenstool galaxy-mass correlation functions occurs on scales of tens of kpc and may suggest that cluster galaxies are more biased tracers of mass than parametric methods generally assume. △ Less

Submitted 15 July, 2016; v1 submitted 31 July, 2015; originally announced July 2015.

Comments: 10 pages, 5 figures, 2 tables; accepted for publication in MNRAS

Journal ref: MNRAS 461 (2016) 2126-2134

Model-free analysis of quadruply imaged gravitationally lensed systems and substructured galaxies

Authors: Addishiwot G. Woldesenbet, Liliya L. R. Williams

Abstract: Multiple image gravitational lens systems, and especially quads are invaluable in determining the amount and distribution of mass in galaxies. This is usually done by mass modeling using parametric or free-form methods. An alternative way of extracting information about lens mass distribution is to use lensing degeneracies and invariants. Where applicable, they allow one to make conclusions about… ▽ More Multiple image gravitational lens systems, and especially quads are invaluable in determining the amount and distribution of mass in galaxies. This is usually done by mass modeling using parametric or free-form methods. An alternative way of extracting information about lens mass distribution is to use lensing degeneracies and invariants. Where applicable, they allow one to make conclusions about whole classes of lenses without model fitting. Here, we use approximate, but observationally useful invariants formed by the three relative polar angles of quad images around the lens center to show that many smooth elliptical+shear lenses can reproduce the same set of quad image angles within observational error. This result allows us to show in a model-free way what the general class of smooth elliptical+shear lenses looks like in the three dimensional (3D) space of image relative angles, and that this distribution does not match that of the observed quads. We conclude that, even though smooth elliptical+shear lenses can reproduce individual quads, they cannot reproduce the quad population. What is likely needed is substructure, with clump masses larger than those responsible for flux ratio anomalies in quads, or luminous or dark nearby perturber galaxies. △ Less

Submitted 18 October, 2015; v1 submitted 20 July, 2015; originally announced July 2015.

Comments: MNRAS, in press; 18 pages, 10 figures, 1 table

Quantifying substructures in {\it Hubble Frontier Field} clusters: comparison with $ΛCDM$ simulations

Authors: Irshad Mohammed, Prasenjit Saha, Liliya L. R. Williams, Jori Liesenborgs, Kevin Sebesta

Abstract: The Hubble Frontier Fields (HFF) are six clusters of galaxies, all showing indications of recent mergers, which have recently been observed for lensed images. As such they are the natural laboratories to study the merging history of galaxy clusters. In this work, we explore the 2D power spectrum of the mass distribution $P_{\rm M}(k)$ as a measure of substructure. We compare $P_{\rm M}(k)$ of thes… ▽ More The Hubble Frontier Fields (HFF) are six clusters of galaxies, all showing indications of recent mergers, which have recently been observed for lensed images. As such they are the natural laboratories to study the merging history of galaxy clusters. In this work, we explore the 2D power spectrum of the mass distribution $P_{\rm M}(k)$ as a measure of substructure. We compare $P_{\rm M}(k)$ of these clusters (obtained using strong gravitational lensing) to that of $Λ$CDM simulated clusters of similar mass. To compute lensing $P_{\rm M}(k)$, we produced free-form lensing mass reconstructions of HFF clusters, without any light traces mass (LTM) assumption. The inferred power at small scales tends to be larger if (i)~the cluster is at lower redshift, and/or (ii)~there are deeper observations and hence more lensed images. In contrast, lens reconstructions assuming LTM show higher power at small scales even with fewer lensed images; it appears the small scale power in the LTM reconstructions is dominated by light information, rather than the lensing data. The average lensing derived $P_{\rm M}(k)$ shows lower power at small scales as compared to that of simulated clusters at redshift zero, both dark-matter only and hydrodynamical. The possible reasons are: (i)~the available strong lensing data are limited in their effective spatial resolution on the mass distribution, (ii)~HFF clusters have yet to build the small scale power they would have at $z\sim 0$, or (iii)~simulations are somehow overestimating the small scale power. △ Less

Submitted 23 March, 2016; v1 submitted 6 July, 2015; originally announced July 2015.

Comments: 13 pages, 10 figures, 1 table; Accepted for publication in MNRAS

Report number: FERMILAB-PUB-16-083-A

Illuminating a Dark Lens : A Type Ia Supernova Magnified by the Frontier Fields Galaxy Cluster Abell 2744

Authors: Steven A. Rodney, Brandon Patel, Daniel Scolnic, Ryan J. Foley, Alberto Molino, Gabriel Brammer, Mathilde Jauzac, Marusa Bradac, Dan Coe, Tom Broadhurst, Jose M. Diego, Or Graur, Jens Hjorth, Austin Hoag, Saurabh W. Jha, Traci L. Johnson, Patrick Kelly, Daniel Lam, Curtis McCully, Elinor Medezinski, Massimo Meneghetti, Julian Merten, Johan Richard, Adam Riess, Keren Sharon , et al. (5 additional authors not shown)

Abstract: SN HFF14Tom is a Type Ia Supernova (SN) discovered at z = 1.3457 +- 0.0001 behind the galaxy cluster Abell 2744 (z = 0.308). In a cosmology-independent analysis, we find that HFF14Tom is 0.77 +- 0.15 magnitudes brighter than unlensed Type Ia SNe at similar redshift, implying a lensing magnification of mu_obs = 2.03 +- 0.29. This observed magnification provides a rare opportunity for a direct empir… ▽ More SN HFF14Tom is a Type Ia Supernova (SN) discovered at z = 1.3457 +- 0.0001 behind the galaxy cluster Abell 2744 (z = 0.308). In a cosmology-independent analysis, we find that HFF14Tom is 0.77 +- 0.15 magnitudes brighter than unlensed Type Ia SNe at similar redshift, implying a lensing magnification of mu_obs = 2.03 +- 0.29. This observed magnification provides a rare opportunity for a direct empirical test of galaxy cluster lens models. Here we test 17 lens models, 13 of which were generated before the SN magnification was known, qualifying as pure "blind tests". The models are collectively fairly accurate: 8 of the models deliver median magnifications that are consistent with the measured mu to within 1-sigma. However, there is a subtle systematic bias: the significant disagreements all involve models overpredicting the magnification. We evaluate possible causes for this mild bias, and find no single physical or methodological explanation to account for it. We do find that model accuracy can be improved to some extent with stringent quality cuts on multiply-imaged systems, such as requiring that a large fraction have spectroscopic redshifts. In addition to testing model accuracies as we have done here, Type Ia SN magnifications could also be used as inputs for future lens models of Abell 2744 and other clusters, providing valuable constraints in regions where traditional strong- and weak-lensing information is unavailable. △ Less

Submitted 16 August, 2015; v1 submitted 22 May, 2015; originally announced May 2015.

Comments: 21 pages, 8 figures, accepted to ApJ

Mass-Galaxy offsets in Abell 3827, 2218 and 1689: intrinsic properties or line-of-sight substructures?

Authors: Irshad Mohammed, Jori Liesenborgs, Prasenjit Saha, Liliya L. R. Williams

Abstract: We have made mass maps of three strong-lensing clusters, Abell 3827, Abell 2218 and Abell 1689, in order to test for mass-light offsets. The technique used is GRALE, which enables lens reconstruction with minimal assumptions, and specifically with no information about the cluster light being given. In the first two of these clusters, we find local mass peaks in the central regions that are displac… ▽ More We have made mass maps of three strong-lensing clusters, Abell 3827, Abell 2218 and Abell 1689, in order to test for mass-light offsets. The technique used is GRALE, which enables lens reconstruction with minimal assumptions, and specifically with no information about the cluster light being given. In the first two of these clusters, we find local mass peaks in the central regions that are displaced from the nearby galaxies by a few to several kpc. These offsets {\em could\/} be due to line of sight structure unrelated to the clusters, but that is very unlikely, given the typical levels of chance line-of-sight coincidences in $ΛCDM$ simulations --- for Abell 3827 and Abell 2218 the offsets appear to be intrinsic. In the case of Abell 1689, we see no significant offsets in the central region, but we do detect a possible line of sight structure: it appears only when sources at $z\ga 3$ are used for reconstructing the mass. We discuss possible origins of the mass-galaxy offsets in Abell 3827 and Abell 2218: these include pure gravitational effects like dynamical friction, but also non-standard mechanisms like self-interacting dark-matter. △ Less

Submitted 17 February, 2014; originally announced February 2014.

Comments: 14 pages, 9 figures; Accepted for publication in MNRAS

Statistical mechanics of collisionless orbits. IV. Distribution of angular momentum

Authors: Liliya L. R. Williams, Jens Hjorth, Radoslaw Wojtak

Abstract: It has been shown in previous work that DARKexp, which is a theoretically derived, maximum entropy, one shape parameter model for isotropic collisionless systems, provides very good fits to simulated and observed dark-matter halos. Specifically, it fits the energy distribution, N(E), and the density profiles, including the central cusp. Here, we extend DARKexp N(E) to include the distribution in a… ▽ More It has been shown in previous work that DARKexp, which is a theoretically derived, maximum entropy, one shape parameter model for isotropic collisionless systems, provides very good fits to simulated and observed dark-matter halos. Specifically, it fits the energy distribution, N(E), and the density profiles, including the central cusp. Here, we extend DARKexp N(E) to include the distribution in angular momentum, L^2, for spherically symmetric systems. First, we argue, based on theoretical, semi-analytical, and simulation results, that while dark-matter halos are relaxed in energy, they are not nearly as relaxed in angular momentum, which precludes using maximum entropy to uniquely derive N(E,L^2). Instead, we require that when integrating N(E,L^2) over squared angular momenta one retrieves the DARKexp N(E). Starting with a general expression for N(E,L^2) we show how the distribution of particles in L^2 is related to the shape of the velocity distribution function, VDF, and velocity anisotropy profile, β(r). We then demonstrate that astrophysically realistic halos, as judged by the VDF shape and β(r), must have linear or convex distributions in L^2, for each separate energy bin. The distribution in energy of the most bound particles must be nearly flat, and become more tilted in favor of radial orbits for less bound particles. These results are consistent with numerical simulations and represent an important step towards deriving the full distribution function for spherically symmetric dark-matter halos. △ Less

Submitted 20 January, 2014; originally announced January 2014.

Comments: 17 pages, 10 figures, accepted to ApJ. The first three in this series are: Paper I: arXiv:1010.0265, Paper II: arXiv:1010.0266, Paper III: arXiv:1010.0267

Journal ref: Astrophysical Journal 783 (2014) 13

Collisionless dynamics in Globular Clusters

Authors: Liliya L. R. Williams, Eric I. Barnes, Jens Hjorth

Abstract: Since globular clusters (GCs) are old, low-N systems their dynamics is widely believed to be fully dominated by collisional two-body processes, and their surface brightness profiles are fit by King models. However, for many GCs, especially those with HST-resolved central regions, and `extra-tidal' features, King models provide poor fits. We suggest that this is partly because collisionless dynamic… ▽ More Since globular clusters (GCs) are old, low-N systems their dynamics is widely believed to be fully dominated by collisional two-body processes, and their surface brightness profiles are fit by King models. However, for many GCs, especially those with HST-resolved central regions, and `extra-tidal' features, King models provide poor fits. We suggest that this is partly because collisionless dynamics is also important and contribute to shaping the cluster properties. We show using time-scale and length-scale arguments that except for the very centers of clusters, collisionless dynamics should be more important than collisional. We then fit 38 GCs analyzed by Noyola and Gebhardt (2006) with (collisional) King and (collisionless) DARKexp models over the full available radial range, and find that the latter provide a better fit to 29 GCs; for six of these the fit is at least ~5x better in term of rms. DARKexp models are theoretically derived maximum entropy equilibrium states of self-gravitating collisionless systems and have already been shown to fit the results of dark matter N-body simulations. (We do not attempt fits with ad hoc fitting functions.) △ Less

Submitted 23 April, 2012; originally announced April 2012.

Comments: 27 pages, 18 figures; accepted to MNRAS

Journal ref: Mon. Not. Roy. Astron. Soc. 423 (2012) 3589-3600