Papers › A Strong Gravitational Lens Is Worth a Thousand Dark Matter Halos: Inference on...

A Strong Gravitational Lens Is Worth a Thousand Dark Matter Halos: Inference on Small-Scale Structure Using Sequential Methods

22 Apr 2024arXiv:2404.14487links table onlyarchive 2025-07-28

Sebastian Wagner-Carena, Jaehoon Lee, Jeffrey Pennington, Jelle Aalbers, Simon Birrer, Risa H. Wechsler

The archive published only this paper's code-link row. Authors, date and abstract are from arXiv's metadata (CC0), read from the Kaggle arXiv metadata snapshot of 2026-09-12 where its title matched the archive's; the title is the archive's.

Strong gravitational lenses are a singular probe of the universe's small-scale structure x2013 they are sensitive to the gravitational effects of low-mass (<10¹⁰ M_⊙) halos even without a luminous counterpart. Recent strong-lensing analyses of dark matter structure rely on simulation-based inference (SBI). Modern SBI methods, which leverage neural networks as density estimators, have shown promise in extracting the halo-population signal. However, it is unclear whether the constraining power of these models has been limited by the methodology or the information content of the data. In this study, we introduce an accelerator-optimized simulation pipeline that can generate lens images with realistic subhalo populations in a matter of milliseconds. Leveraging this simulator, we identify the main methodological limitation of our fiducial SBI analysis: training set size. We then adopt a sequential neural posterior estimation (SNPE) approach, allowing us to iteratively refine the distribution of simulated training images to better align with the observed data. Using only one-fifth as many mock Hubble Space Telescope (HST) images, SNPE matches the constraints on the low-mass halo population produced by our best non-sequential model. Our experiments suggest that an over three order-of-magnitude increase in training set size and GPU hours would be required to achieve an equivalent result without sequential methods. While the full potential of the existing strong lens sample remains to be explored, the notable improvement in constraining power enabled by our sequential approach highlights that the current constraints are limited primarily by methodology and not the data itself. Moreover, our results emphasize the need to treat training set generation and model optimization as interconnected stages of any cosmological analysis using simulation-based inference techniques.

PaperPDFCodeCode Syntology ran

In Syntology Open this paper in Syntology's Atlas, the map of the papers in Syntology's graph and their citations.

For agents, Syntology's MCP tool lists every function and class Syntology harvested from this paper and whether it ran (how to connect): get_harvested_code_for_paper(arxiv_id="2404.14487")

Code

Syntology Ran 13 of 13 code samples harvested from 1 repository linked to this paper; 0 have no recorded run. Of those that ran: 13 ran with no contract checked.

By repository: official repository: 13 samples from 1 repository, 13 ran. The run record, sample by sample. “Ran” means executed on a synthesized input, not that the code is correct or reproduces the paper.

swagnercarena/paltax officialmentioned in papermentioned on GitHubjaxApache-2.0 report

Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.

Code Syntology ran Syntology

13 samples harvested; 13 ran; 0 honoured the contract we drafted; 0 have no recorded run. Read from Syntology's graph 2026-09-24; that is when this build read the record, not when the samples ran.

13ran

Licence: 0 of the 13 samples are pointer only, meaning Syntology does not serve that copy's text. This page shows no code text for any sample; each one links to its file in the repository.

Harvested from swagnercarena/paltax. “Ran” means the sample executed on a synthesized input. It does not mean the output is correct, and nothing here reproduces the paper's results. “Honoured” and “violated” refer to a contract Syntology drafted from the code itself; “our draft was wrong” and “fixture could not drive it” are failures of Syntology's instrument, not of the code.

Each sample ends with its code_sha256, Syntology's identity for that exact code. An agent fetches the stored sample with Syntology's MCP tool get_code(code_sha256="…") (how to connect); click an identity to copy that call.

Repository labels, per sample. official repository: The archive marks this repository official for the paper. named in the paper: The archive records that the paper mentions this repository; it is not marked official. community (archive-listed): In the archive's code links for this paper, not marked official and not recorded as mentioned in the paper. found in paper text by Syntology: Syntology found this repository in the paper's own text; whether it is the authors' implementation is not asserted. community: Not in the archive's code links for this paper; a community repository Syntology harvested. Samples from a repository marked official are listed first. Licence labels name the repository's licence as recorded at harvest. “Pointer only” means Syntology does not serve that copy's text, for one of four reasons: no licence file was found; the licence was not identified; the licence is recorded as permissive but that copy's record is not marked cleared; or the licence is outside the permissive list Syntology serves text under (MIT, Apache-2.0, BSD and similar). Some licences outside that list permit redistribution, such as WTFPL, and GPL-3.0 under its conditions; they are simply not on the list. Hover a licence label for the reason. File links open the file on GitHub at the default branch, which may have changed since the harvest.

angular_diameter_distance swagnercarena/paltax/paltax/cosmology_utils.py official repository ran Apache-2.0 (permissive) · 46dfe4a99d641bb0 · report
comoving_distance swagnercarena/paltax/paltax/cosmology_utils.py official repository ran Apache-2.0 (permissive) · 23876623204ee870 · report
cored_nfw_draws swagnercarena/paltax/paltax/nfw_functions.py official repository ran Apache-2.0 (permissive) · aabe89bafcb5d03e · report
encode_constant swagnercarena/paltax/paltax/input_pipeline.py official repository ran Apache-2.0 (permissive) · b38e1a2c3f4f88cf · report
encode_normal swagnercarena/paltax/paltax/input_pipeline.py official repository ran Apache-2.0 (permissive) · efddc701eaed6cfa · report
encode_uniform swagnercarena/paltax/paltax/input_pipeline.py official repository ran Apache-2.0 (permissive) · cc194b3ddf6e1b17 · report
gaussian_product_analytical swagnercarena/paltax/paltax/hierarchical_inference.py official repository ran Apache-2.0 (permissive) · 0ec30f4c247665b0 · report
log_p_omega swagnercarena/paltax/paltax/hierarchical_inference.py official repository ran Apache-2.0 (permissive) · 87414ad5fa237b5f · report
nu_function swagnercarena/paltax/paltax/los.py official repository ran fingerprinted Apache-2.0 (permissive) · 5a8520efd4895d45 · report
power_law_draw swagnercarena/paltax/paltax/power_law.py official repository ran Apache-2.0 (permissive) · a4f6e55a048d7d39 · report
power_law_integrate swagnercarena/paltax/paltax/power_law.py official repository ran fingerprinted Apache-2.0 (permissive) · 255476f8621cb6d6 · report
psf_convolution swagnercarena/paltax/paltax/image_simulation.py official repository ran Apache-2.0 (permissive) · c3b548c9c454e1cb · report
suppressed_power_law_draw swagnercarena/paltax/paltax/power_law.py official repository ran Apache-2.0 (permissive) · ca9fea87be931ee0 · report

Results from the paper archive 2025-07-28

No leaderboard rows for this paper in the archive.

Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections