-
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.
-
The lens SW05 J143454.4+522850: a fossil group at redshift 0.6?
Authors:
Philipp Denzel,
Onur Çatmabacak,
Jonathan P. Coles,
Claude Cornen,
Robert Feldmann,
Ignacio Ferreras,
Xanthe Gwyn Palmer,
Rafael Küng,
Dominik Leier,
Prasenjit Saha,
Aprajita Verma
Abstract:
Fossil groups are considered the end product of natural galaxy group evolution in which group members sink towards the centre of the gravitational potential due to dynamical friction, merging into a single, massive, and X-ray bright elliptical. Since gravitational lensing depends on the mass of a foreground object, its mass concentration, and distance to the observer, we can expect lensing effects…
▽ More
Fossil groups are considered the end product of natural galaxy group evolution in which group members sink towards the centre of the gravitational potential due to dynamical friction, merging into a single, massive, and X-ray bright elliptical. Since gravitational lensing depends on the mass of a foreground object, its mass concentration, and distance to the observer, we can expect lensing effects of such fossil groups to be particularly strong. This paper explores the exceptional system $\mathrm{J}143454.4+522850$. We combine gravitational lensing with stellar population-synthesis to separate the total mass of the lens into stars and dark matter. The enclosed mass profiles are contrasted with state-of-the-art galaxy formation simulations, to conclude that SW05 is likely a fossil group with a high stellar to dark matter mass fraction $0.027\pm0.003$ with respect to expectations from abundance matching $0.012\pm0.004$, indicative of a more efficient conversion of gas into stars in fossil groups.
△ Less
Submitted 7 April, 2021;
originally announced April 2021.
-
A new strategy for matching observed and simulated lensing galaxies
Authors:
Philipp Denzel,
Sampath Mukherjee,
Prasenjit Saha
Abstract:
The study of strong-lensing systems conventionally involves constructing a mass distribution that can reproduce the observed multiply-imaging properties. Such mass reconstructions are generically non-unique. Here, we present an alternative strategy: instead of modelling the mass distribution, we search cosmological galaxy-formation simulations for plausible matches. In this paper we test the idea…
▽ More
The study of strong-lensing systems conventionally involves constructing a mass distribution that can reproduce the observed multiply-imaging properties. Such mass reconstructions are generically non-unique. Here, we present an alternative strategy: instead of modelling the mass distribution, we search cosmological galaxy-formation simulations for plausible matches. In this paper we test the idea on seven well-studied lenses from the SLACS survey. For each of these, we first pre-select a few hundred galaxies from the EAGLE simulations, using the expected Einstein radius as an initial criterion. Then, for each of these pre-selected galaxies, we fit for the source light distribution, while using MCMC for the placement and orientation of the lensing galaxy, so as to reproduce the multiple images and arcs. The results indicate that the strategy is feasible, and even yields relative posterior probabilities of two different galaxy-formation scenarios, though these are not statistically significant yet. Extensions to other observables, such as kinematics and colours of the stellar population in the lensing galaxy, is straightforward in principle, though we have not attempted it yet. Scaling to arbitrarily large numbers of lenses also appears feasible. This will be especially relevant for upcoming wide-field surveys, through which the number of galaxy lenses will rise possibly a hundredfold, which will overwhelm conventional modelling methods.
△ Less
Submitted 19 February, 2021;
originally announced February 2021.
-
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.
-
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.
-
Lessons from a blind study of simulated lenses: image reconstructions do not always reproduce true convergence
Authors:
Philipp Denzel,
Sampath Mukherjee,
Jonathan P. Coles,
Prasenjit Saha
Abstract:
In the coming years, strong gravitational lens discoveries are expected to increase in frequency by two orders of magnitude. Lens-modelling techniques are being developed to prepare for the coming massive influx of new lens data, and blind tests of lens reconstruction with simulated data are needed for validation. In this paper we present a systematic blind study of a sample of 15 simulated strong…
▽ More
In the coming years, strong gravitational lens discoveries are expected to increase in frequency by two orders of magnitude. Lens-modelling techniques are being developed to prepare for the coming massive influx of new lens data, and blind tests of lens reconstruction with simulated data are needed for validation. In this paper we present a systematic blind study of a sample of 15 simulated strong gravitational lenses from the EAGLE suite of hydrodynamic simulations. We model these lenses with a free-form technique and evaluate reconstructed mass distributions using criteria based on shape, orientation, and lensed image reconstruction. Especially useful is a lensing analogue of the Roche potential in binary star systems, which we call the $\textit{lensing Roche potential}$. This we introduce in order to factor out the well-known problem of steepness or mass-sheet degeneracy. Einstein radii are on average well recovered with a relative error of ${\sim}5\%$ for quads and ${\sim}25\%$ for doubles; the position angle of ellipticity is on average also reproduced well up to $\pm10^{\circ}$, but the reconstructed mass maps tend to be too round and too shallow. It is also easy to reproduce the lensed images, but optimising on this criterion does not guarantee better reconstruction of the mass distribution.
△ Less
Submitted 17 June, 2020; v1 submitted 22 October, 2019;
originally announced October 2019.
-
Molecular dynamics simulations of bubble nucleation in dark matter detectors
Authors:
Philipp Denzel,
Jürg Diemand,
Raymond Angélil
Abstract:
Bubble chambers and droplet detectors used in dosimetry and dark matter particle search experiments use a superheated metastable liquid in which nuclear recoils trigger bubble nucleation. This process is described by the classical heat spike model of F. Seitz [Phys. Fluids (1958-1988) 1, 2 (1958)], which uses classical nucleation theory to estimate the amount and the localization of the deposited…
▽ More
Bubble chambers and droplet detectors used in dosimetry and dark matter particle search experiments use a superheated metastable liquid in which nuclear recoils trigger bubble nucleation. This process is described by the classical heat spike model of F. Seitz [Phys. Fluids (1958-1988) 1, 2 (1958)], which uses classical nucleation theory to estimate the amount and the localization of the deposited energy required for bubble formation. Here we report on direct molecular dynamics simulations of heat-spike-induced bubble formation. They allow us to test the nanoscale process described in the classical heat spike model. 40 simulations were performed, each containing about 20 million atoms, which interact by a truncated force-shifted Lennard-Jones potential. We find that the energy per length unit needed for bubble nucleation agrees quite well with theoretical predictions, but the allowed spike length and the required total energy are about twice as large as predicted. This could be explained by the rapid energy diffusion measured in the simulation: contrary to the assumption in the classical model, we observe significantly faster heat diffusion than the bubble formation time scale. Finally we examine α-particle tracks, which are much longer than those of neutrons and potential dark matter particles. Empirically, α events were recently found to result in louder acoustic signals than neutron events. This distinction is crucial for the background rejection in dark matter searches. We show that a large number of individual bubbles can form along an α track, which explains the observed larger acoustic amplitudes.
△ Less
Submitted 27 January, 2016;
originally announced January 2016.