Researcher profile

Liliya L. R. Williams

Liliya L. R. Williams contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

9 published item(s)

preprint2025arXiv

Transient star B/R ratio and star formation in $z\gtrsim 1$ lensed galaxies

The extreme magnification from galaxy clusters and microlenses therein allows the detection of individual, luminous stars in lensed galaxies as transient events, and hence provides a valuable window into the high mass stellar population in $z>1$ galaxies. As these bright stars can only be formed at specific ages, the relative abundance of transient events at blue (B) and red (R) optical wavelengths ($B/R$ ratio) can provide insights into the recent star formation history of galaxies that are not well constrained by their spectral energy distributions (SEDs). Here, we forward model the transient detection rates in an idealized mock scenario to find that the $B/R$ ratio of strongly lensed $z>1$ galaxies decreases quickly with increasing age. This ratio has moderate sensitivity to metallicity and comparatively low sensitivity to dust attenuation, with no significant dependency on the stellar initial mass function. Fitting model stellar populations to either the SED or $B/R$ ratio alone of ``Warhol'' arc ($z = 0.94$), we find that neither a simple single starburst nor a more complex star formation can simultaneously reproduce both constraints. We then demonstrate that a best-fit model constrained by both the B/R ratio and SED requires a star-formation rate that has varied quite dramatically over the past $\sim$50 Myr, for which the total stellar mass formed over this time is a factor of 10 (with $2-3σ$ significance) different from the best-fit models to the SED alone. Our work shows that the transient $B/R$ ratio can be used as an additional powerful constraint on the recent star formation history of higher-redshift galaxies in future works that are strongly lensed by galaxy clusters.

preprint2021arXiv

The Hubble constant from eight time-delay galaxy lenses

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

The BUFFALO HST Survey

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.

preprint2020arXiv

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

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.

preprint2019arXiv

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

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.

preprint2019arXiv

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

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.