CosmoLensNRE
Cosmology Inference from Strong Gravitational Lensing using Neural Ratio Estimation
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Cosmology Inference from Strong Gravitational Lensing using Neural Ratio Estimation
Cosmological observations and galaxy dynamics seem to imply that five out of six parts in mass of all matter in the universe is composed of dark matter, which is not accounted for by the Standard Model of particles. The cold dark matter (CDM) paradigm has been extremely successful at describing observations on large-cosmological scale. However, many different dark matter candidates have the same observable effects as CDM on large-length scales. One possible avenue to distinguish between these different models is to look on much smaller-length scales, where usually dark matter models become distinguishable. This research developed the theoretical framework and statistical tools needed to map the detailed distribution of dark matter on subgalactic scales using strong gravitational lensing. A second goal of this research was to build new statistical techniques that efficiently exploit the full power of cosmic-structure data from next-generation surveys. Galaxy clustering on large scales provides significant cosmological information through the power spectrum. Additional information can be gained from higher-order statistics. This research provided new ways to understand the initial conditions of the universe from current and upcoming cosmological surveys.
Hierarchical neural posterior estimation has its place Hierarchical Bayesian Modeling (HBM) combined with MCMC algorithms has been shown to provide more robust and accurate inference for real-world phenomena in which nature takes a nested form. However, MCMC-based inference can be computationally expensive, and its performance often suffers for complex posterior geometries. These costs are especially pertinent for HBM. Studies have recently demonstrated the potential for a flexible, expressive, and amortized hierarchical neural posterior estimator (HNPE) built on Normalizing Flows. These studies have mostly been performed on simple datasets, or they focus on a single parameter from each level of the hierarchy. A systematic study analyzing how both hierarchical methods compare for more complex and realistic datasets is necessary before applying HNPE for scientific measurements. Here, we re-explore the theory behind HNPE and conduct comparative numerical experiments of HNPE and MCMC-based HBM methods on real and synthetic data, including strong gravitational lensing simulations. In particular, we use a suite of diagnostics to show trade-offs in terms of accuracy, precision, time to train or sample, reproducibility, and the need for expert domain knowledge. Especially for higher dimensional and complex posteriors, HNPE is expected to drastically improve on time for inference, accuracy, and precision with an upfront training time cost.
Strong gravitational lenses are a rare and instructive type of astronomical object. Identification has long relied on serendipity, but different strategies—such as mixed spectroscopy of multiple galaxies along the line of sight, machine-learning algorithms, and citizen science—have been employed to identify these objects as new imaging surveys become available. We report on the comparison between spectroscopic, machine-learning, and citizen-science identification of galaxy–galaxy lens candidates from independently constructed lens catalogs in the common survey area of the equatorial fields of the Galaxy and Mass Assembly survey. In these, we have the opportunity to compare high completeness spectroscopic identifications against high-fidelity imaging from the Kilo Degree Survey used for both machine-learning and citizen-science lens searches. We find that the three methods—spectroscopy, machine learning, and citizen science—identify 47, 47, and 13 candidates, respectively, in the 180 square degrees surveyed. These identifications barely overlap, with only two identified by both citizen science and machine learning. We have traced this discrepancy to inherent differences in the selection functions of each of the three methods, either within their parent samples (i.e., citizen science focuses on low redshift) or inherent to the method (i.e., machine learning is limited by its training sample and prefers well-separated features, while spectroscopy requires sufficient flux from lensed features to lie within the fiber). These differences manifest as separate samples in estimated Einstein radius, lens stellar mass, and lens redshift. The combined sample implies a lens candidate sky density of ∼0.59 deg{sup −2} and can inform the construction of a training set spanning a wider mass–redshift space. A combined approach and refinement of automated searches would result in a more complete sample of galaxy–galaxy lens candidates for future surveys.
Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov and conservation of the Walker–Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images.
Light passing near a black hole can follow multiple paths from an emission source to an observer due to strong gravitational lensing. Photons following different paths take different amounts of time to reach the observer, which produces an echo signature in the image. The characteristic echo delay is determined primarily by the mass of the black hole, but it is also influenced by the black hole spin and inclination to the observer. In the Kerr geometry, echo images are demagnified, rotated, and sheared copies of the direct image and lie within a restricted region of the image. Echo images have exponentially suppressed flux, and temporal correlations within the flow make it challenging to directly detect light echoes from the total light curve. In this Letter, we propose a novel method to search for light echoes by correlating the total light curve with the interferometric signal at high spatial frequencies, which is a proxy for indirect emission. We explore the viability of our method using numerical general relativistic magnetohydrodynamic simulations of a near-face-on accretion system scaled to M87-like parameters. We demonstrate that our method can be used to directly infer the echo delay period in simulated data. An echo detection would be clear evidence that we have captured photons that have circled the black hole, and a high-fidelity echo measurement would provide an independent measure of fundamental black hole parameters. Our results suggest that detecting echoes may be achievable through interferometric observations with a modest space-based very long baseline interferometry mission.
We present near-infrared spectroscopy and imaging of the compact steep-spectrum radio source 3C 318 which shows it to be a quasar at redshift z = 1.574 (the z = 0.752 value previously reported is incorrect). 3C 318 is an IRAS, ISO and SCUBA source so its new redshift makes it the most intrinsically luminous far-infrared (FIR) source in the 3C catalogue (there is no evidence of strong gravitational lensing effects). Its bolometric luminosity greatly exceeds the 10(exp 13) solar luminosity level above which an object is said to be hyperluminous. Its spectral energy distribution (SED) requires that the quasar heats the dust responsible for the FIR flux, as is believed to be the case in other hyperluminous galaxies, and contributes (at the greater than 10% level) to the heating of the CIA dust responsible for the sub-mm emission. We cannot determine whether a starburst makes an important contribution to the heating of the coolest dust, so evidence for a high star-formation rate is circumstantial being based on the high dust, and hence gas, C-1 mass required by its sub-mm detection. We show that the current sub-mm and FIR data available for the highest-redshift radio galaxies are consistent with SEDs similar to that of 3C 318. This indicates that at least some of this population may be detected in the sub-mm because of dust heated by the quasar nucleus, and that interpreting sub-mm detection as evidence for very high (approx. less than 1000 solar mass/yr) star-formation rates may not always be valid. We show that the 3C318 quasar is slightly reddened (A(sub v) approx. = 0.5), the most likely cause of which is SMC-type dust in the host galaxy. If very distant radio galaxies are reddened in a similar way then we show that only slightly greater amounts of dust could obscure the quasars in these sources. We speculate that the low fraction of quasars amongst the very high redshift (z approx. greater than 3) objects in low-frequency radio-selected samples is the result of such obscuration. The highest-z objects might be preferentially obscured because like 3C318 they are inevitably observed very shortly after the jet-triggering event, or because their host galaxies are richer in dust and gas at earlier cosmic epochs, or because of some combination of these two effects.
The Sloan Lens ACS (SLACS) Survey is an efficient Hubble Space Telescope (HST) Snapshot imaging survey for new galaxy-scale strong gravitational lenses. The targeted lens candidates are selected spectroscopically from the Sloan Digital Sky Survey (SDSS) database of galaxy spectra for having multiple nebular emission lines at a redshift significantly higher than that of the SDSS target galaxy. The SLACS survey is optimized to detect bright early-type lens galaxies with faint lensed sources in order to increase the sample of known gravitational lenses suitable for detailed lensing, photometric, and dynamical modeling. In this paper, the first in a series on the current results of our HST Cycle 13 imaging survey, we present a catalog of 19 newly discovered gravitational lenses, along with nine other observed candidate systems that are either possible lenses, nonlenses, or nondetections. The survey efficiency is thus >=68%. We also present Gemini 8 m and Magellan 6.5 m integral-field spectroscopic data for nine of the SLACS targets, which further support the lensing interpretation. A new method for the effective subtraction of foreground galaxy images to reveal faint background features is presented. We show that the SLACS lens galaxies have colors and ellipticities typical of the spectroscopic parent sample from which they are drawn (SDSS luminous red galaxies and quiescent MAIN sample galaxies), but are somewhat brighter and more centrally concentrated. Several explanations for the latter bias are suggested. The SLACS survey provides the first statistically significant and homogeneously selected sample of bright early-type lens galaxies, furnishing a powerful probe of the structure of early-type galaxies within the half-light radius. The high confirmation rate of lenses in the SLACS survey suggests consideration of spectroscopic lens discovery as an explicit science goal of future spectroscopic galaxy surveys.
This thesis addresses the intersection of dark matter and galaxy formation physics by modeling the faintest galaxies in the Universe in a cosmological context. In the first part of this work, I present an empirical model for the connection between faint galaxies and small dark matter halos. Combining this model with tailored cosmological simulations of Milky Way analogs and state-of-the-art satellite galaxy observations yields new insights into dwarf galaxy formation and hierarchical structure formation. Next, I show that these same dwarf galaxy observations place stringent constraints on microphysical dark matter properties including its warmth, self-interactions, Standard Model interactions, minimum particle mass, and formation epoch. Finally, I describe a unified framework that combines dark matter constraints from dwarf galaxies and strong gravitational lensing, paving the way for next-generation surveys of small-scale cosmic structure to deliver unprecedented insights into dark matter physics.
We report a new galaxy cluster catalog, selected using the thermal Sunyaev-Zel'dovich (SZ) effect, from 5 years of observations of the SPT-3G Main field. Drawn from arcminute-resolution data with white noise levels of 3.2, 2.5, and 8.9 $μ$K-arcmin at 95, 150, and 220 GHz, respectively, the sample consists of 8,892 cluster candidates detected above significance $ξ=4$, with an expected purity of $>82\%$ (4,480 at $ξ\ge5$ with purity $>99\%$). Using optical and infrared data we have confirmed 7,190 candidates as clusters. The sample spans a mass range $7.9 \times 10^{13}M_\odot/h_{70}< M_\textrm{500c} < 1.6 \times 10^{15}M_\odot/h_{70}$ with a median mass of $1.65 \times 10^{14}M_\odot/h_{70}$, and a redshift range of $0.0371$ and 271 at $z>1.5$. Compared to previous SZ cluster samples from South Pole Telescope and Atacama Cosmology Telescope data, the SPT-3G sample is highly consistent in mass and redshift but is significantly deeper, with per-cluster detection signal-to-noise 2-4 times higher and a cluster density of 4.5 confirmed clusters/deg$^2$. We cross match with eRASS1 cluster and point source catalogs, finding 1,279 and 1,319 matches, respectively. The SPT and eROSITA cluster mass estimates are in relatively good agreement. We perform a series of validation checks using both internal data splits and comparisons to external samples. These tests show increasing correlated (dusty) emission with redshift, with a $\sim17\times$ larger 220 GHz temperature increment for clusters at $z\sim1.5$ than $z\sim0.25$, but only weak evidence for correlated synchrotron emission. Finally, a number of clusters are flagged as candidate strong gravitational lenses.
Here, we search for signatures of gravitational lensing in the gravitational-wave signals from compact binary coalescences detected by Advanced Laser Interferometer Gravitational-wave Observatory (LIGO) and Advanced Virgo during O3a, the first half of their third observing run. We study: (1) the expected rate of lensing at current detector sensitivity and the implications of a non-observation of strong lensing or a stochastic gravitational-wave background on the merger-rate density at high redshift; (2) how the interpretation of individual high-mass events would change if they were found to be lensed; (3) the possibility of multiple images due to strong lensing by galaxies or galaxy clusters; and (4) possible wave-optics effects due to point-mass microlenses. Several pairs of signals in the multiple-image analysis show similar parameters and, in this sense, are nominally consistent with the strong lensing hypothesis. However, taking into account population priors, selection effects, and the prior odds against lensing, these events do not provide sufficient evidence for lensing. Overall, we find no compelling evidence for lensing in the observed gravitational-wave signals from any of these analyses.
The goal of this Letter is to develop a machine learning model to analyze the main gravitational lens and detect dark substructure (subhalos) within simulated images of strongly lensed galaxies. Using the technique of image segmentation, we turn the task of identifying subhalos into a classification problem, where we label each pixel in an image as coming from the main lens, a subhalo within a binned mass range, or neither. Our network is only trained on images with a single smooth lens and either zero or one subhalo near the Einstein ring. On an independent test set with lenses with large ellipticities, quadrupole and octopole moments, and for source apparent magnitudes between 17–25, the area of the main lens is recovered accurately. On average, only 1.3% of the true area is missed and 1.2% of the true area is added to another part of the lens. In addition, subhalos as light as 10 8.5 M ⊙ can be detected if they lie in bright pixels along the Einstein ring. Furthermore, the model is able to generalize to new contexts it has not been trained on, such as locating multiple subhalos with varying masses or more than one large smooth lens.
We present a strong lensing analysis and reconstruct the mass distribution of SPT-CL J0356−5337, a galaxy cluster at redshift z = 1.034. Our model supersedes previous models by making use of new multiband Hubble Space Telescope data and Multi-Unit Spectroscopic Explorer (MUSE) spectroscopy. We identify two additional lensed galaxies to inform a more well-constrained model using 12 sets of multiple images in five separate lensed sources. The three previously known sources were spectroscopically confirmed by G. Mahler et al. at redshifts of z = 2.363, z = 2.364, and z = 3.048. We measured the spectroscopic redshifts of two of the newly discovered arcs using MUSE data, at z = 3.0205 and z = 5.3288. We increase the number of cluster member galaxies by a factor of 3 compared to previous work. We also report the detection of extended Lyα emission from several background galaxies. We measure the total projected mass density of the two major subcluster components, one dominated by the brightest cluster galaxy and the other by a compact group of luminous red galaxies. We find ${M}_{{\rm{B}}{\rm{C}}{\rm{G}}}(\lt 80\,\rm{kpc})=3.9{3}_{-0.14}^{+0.21}\times 1{0}^{13}$ M⊙ and ${M}_{{\rm{LRG}}}(\lt 80\,\rm{kpc})=2.9{2}_{-0.23}^{+0.16}\times 1{0}^{13}$ M⊙, yielding a mass ratio of $1.3{5}_{-0.08}^{+0.16}$ . The strong lensing constraints offer a robust estimate of the projected mass density regardless of modeling assumptions; allowing more substructure in this line of sight does not change the results or conclusions. Our results corroborate the conclusion that SPT-CL J0356−5337 is dominated by two mass components and is likely undergoing a major merger on the plane of the sky.
We present COOL J1323+0343, an early-type galaxy at z = 1.0153 ± 0.0006, strongly lensed by a cluster of galaxies at z = 0.353 ± 0.001. This object was originally imaged by DECaLS and noted as a gravitational lens by COOL-LAMPS, a collaboration initiated to find strong-lensing systems in recent public optical imaging data, and confirmed with follow-up data. With ground-based grzH imaging and optical spectroscopy from the Las Campanas Observatory and the Nordic Optical Telescope, we derive a stellar mass, metallicity, and star formation history from stellar-population synthesis modeling. The lens modeling implies a total magnification, summed over the three images in the arc, of μ ~ 113. The stellar mass in the source plane is M * ~ 10.64 M ⊙ and the 1σ upper limit on the star formation rate (SFR) in the source plane is SFR ~ 3.75 × 10 -2 M ⊙ yr -1 (log sSFR = -12.1 yr -1 ) in the youngest two age bins (0–100 Myr), closest to the epoch of observation. Our measurements place COOL J1323+0343 below the characteristic mass of the stellar mass function, making it an especially compelling target that could help clarify how intermediate-mass quiescent galaxies evolve. We reconstruct COOL J1323+0343 in the source plane and fit its light profile. This object is below the expected size evolution of an early-type galaxy at this mass with an effective radius r e ~ 0.5 kpc. This extraordinarily magnified and bright lensed early-type galaxy offers an exciting opportunity to study the morphology and star formation history of an intermediate-mass early-type galaxy in detail at z ~ 1.
We present the discovery of SN 2025mkn, a gravitationally lensed Type II supernova. First detected as a blue transient in Zwicky Transient Facility (ZTF), 0. "83 from a z = 0.42 elliptical galaxy, the follow-up SNIFS/UH2.2 m and LRIS/Keck spectra revealed absorption lines at z = 1.371. Later JWST NIRCam imaging shows that the bright transient is a close pair of point sources separated by 0. "07, and a 30 times fainter counterimage opposite the lens, for which NIRSpec reveals strong Hα emission also at z = 1.371. The lightcurves and spectra are consistent with the Type II supernova source being magnified ≳100 times, with ∼250 required to reconcile its luminosity with that of nearby events such as SN 2023ixf. Lens models are consistent with such high magnifications, and always show that the faint image arrived first (undetected in earlier ZTF imaging), consistent with the later spectral phase of this fainter image. A fourth image is also predicted and possibly detected in the NIRSpec data. Lightcurve-based time-delay measurements are not possible due to the first image being the faintest; however, the resolved NIRSpec spectra offer a future opportunity for time-delay cosmography through supernova phase measurements.
Supernovae (SNe) that have been multiply imaged by gravitational lensing are rare and powerful probes for cosmology. Each detection is an opportunity to develop the critical tools and methodologies needed as the sample of lensed SNe increases by orders of magnitude with the upcoming Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope. The latest such discovery is of the quadruply imaged Type Ia SN 2022qmx (aka, "SN Zwicky") at z = 0.3544. SN Zwicky was discovered by the Zwicky Transient Facility in spatially unresolved data. Here we present follow-up Hubble Space Telescope observations of SN Zwicky, the first from the multicycle "LensWatch (www.lenswatch.org)" program. We measure photometry for each of the four images of SN Zwicky, which are resolved in three WFC3/UVIS filters (F475W, F625W, and F814W) but unresolved with WFC3/IR F160W, and present an analysis of the lensing system using a variety of independent lens modeling methods. We find consistency between lens-model-predicted time delays (≲1 day), and delays estimated with the single epoch of Hubble Space Telescope colors (≲3.5 days), including the uncertainty from chromatic microlensing (~1–1.5 days). Our lens models converge to an Einstein radius of ${\theta }_{{\rm{E}}}=({0.168}_{-0.005}^{+0.009})^{\prime\prime} $, the smallest yet seen in a lensed SN system. The "standard candle" nature of SN Zwicky provides magnification estimates independent of the lens modeling that are brighter than predicted by $\sim {1.7}_{-0.6}^{+0.8}$ mag and $\sim {0.9}_{-0.6}^{+0.8}$ mag for two of the four images, suggesting significant microlensing and/or additional substructure beyond the flexibility of our image-position mass models.
Abstract Gravitational lensing by galaxy clusters has emerged as a powerful tool to probe the standard cold dark matter (CDM) paradigm of structure formation in the Universe. Despite the remarkable explanatory power of CDM on large scales, tensions with observations on small scales have been reported. Recent studies find that the observational cross section for galaxy–galaxy strong lensing (GGSL) in clusters exceeds the CDM prediction by more than an order of magnitude and persists even after rigorous examination of various possible systematics. We investigate the impact of modifying the internal structure of cluster dark matter subhalos on GGSL and report that altering the inner density profile, given by r γ , to steeper slopes with γ > 2.5 can alleviate the GGSL discrepancy. This is steeper than slopes obtained with the inclusion of the contribution of baryons to the inner regions of these subhalos. Deviating from the γ ~ 1.0 cusps that CDM predicts, these steeper slopes could arise in models of self-interacting dark matter undergoing core collapse. Our results motivate additional study of subhalo core collapse in dense cluster environments.
We report the combined results of eight searches for strong gravitational lens systems in the full 5000 square degrees of Dark Energy Survey (DES) observations. The observations accumulated by the end of the third observing season fully covered the DES footprint in five filters (grizY), with an i-band limiting magnitude (at 10σ) of 23.44. In four searches, a list of potential candidates was identified using a color and magnitude selection from the object catalogs created from the first three observing seasons. Three other searches were conducted at the locations of previously identified galaxy clusters. Cutout images of potential candidates were then visually scanned using an object viewer. An additional set of candidates came from a data-quality check of a subset of the color–coadd tiles created from the full DES six-season data set. A short list of the most promising strong-lens candidates was then numerically ranked according to whether or not we judged them to be bona fide strong gravitational lens systems. These searches discovered a diverse set of 247 strong-lens candidate systems, of which 81 are identified for the first time. We provide the coordinates, magnitudes, and photometric properties of the lens and source objects, and an estimate of the Einstein radius for 81 new systems and 166 previously reported systems. This catalog will be of use for selecting interesting systems for detailed follow up, studies of galaxy cluster and group mass profiles, as well as a training/validation set for automated strong-lens searches.