Inferring subhalo effective density slopes from strong lensing observations with neural likelihood-ratio estimation

The Big Picture

Imagine trying to understand the contents of a room by watching shadows on a wall. You can’t see the objects directly, only the distorted shapes they cast. At cosmic scales, galaxies bend light from even more distant galaxies into rings, arcs, and multiple images. This is strong gravitational lensing, and it’s one of the few tools we have for probing one of the universe’s deepest mysteries: dark matter.

Dark matter makes up about 84% of all matter in the universe, yet it emits no light, reflects none, and interacts with ordinary matter only through gravity. We know it exists from its gravitational pull, but figuring out its fundamental nature requires looking at the smallest scales, far below individual galaxies. Sub-galactic clumps called subhalos contain no stars, making them essentially invisible except through their gravitational fingerprints on passing light. Strong lensing can detect exactly that.

The problem: existing methods are either imprecise or grind to a halt computationally when analyzing more than one or two subhalos at a time. Researchers at Harvard and MIT have introduced a machine learning approach (a neural network trained to judge which dark matter model best fits an observed lensing image) that analyzes entire populations of subhalos simultaneously, pulling out statistical power previously out of reach.

Key Insight: By replacing traditional sampling with a trained neural network, this approach combines information from many subhalos across multiple lensing images at once, enabling population-level inferences about dark matter that were previously computationally infeasible.

How It Works

Most previous studies measured the subhalo mass function, a census of subhalo masses, but doing so requires assumptions about internal density profiles in regions not directly observed. That introduces systematic errors. The team instead focused on the effective density slope: how steeply a subhalo’s density falls off from its center, measured at the scale where lensing is most sensitive. This slope can be extracted directly from lensing data without assuming any particular profile shape. It’s cleaner and encodes complementary information about dark matter physics that pure mass measurements miss.

The forward model works like this:

1. Simulate a lensing system: Use lenstronomy to generate mock 100×100 pixel images at 0.08 arcseconds per pixel, including a background source galaxy, a foreground lens galaxy, and a population of dark matter subhalos drawn from some distribution.
2. Vary the population parameter: The key parameter is the mean effective density slope, a number that distinguishes cold dark matter (cuspy, concentrated profiles) from alternatives like warm or self-interacting dark matter (cored, less concentrated profiles).
3. Train a neural network: The network learns to estimate the likelihood ratio between two hypotheses: that an image came from a population with slope θ₁ versus slope θ₂.
4. Deploy across many images: After training, the network processes new lensing images rapidly, combining evidence from each without multiplying computational cost.

A lensing image contains enormous variation in individual subhalo details: exact masses, positions, individual slopes. None of that directly reveals the population-level pattern. The network’s key trick is that it bypasses measuring any of those individual properties. It learns to extract only the signal that distinguishes one dark matter population from another. This is simulation-based inference at its most powerful: train on realistic simulations, then deploy the trained network on new observations.

Why It Matters

The timing here is no coincidence. Upcoming sky surveys from the Euclid space telescope and the Rubin Observatory’s Legacy Survey of Space and Time are expected to increase known strong lensing systems from the current few hundred to potentially hundreds of thousands. Each system contains information about dark matter substructure. Without a computationally efficient analysis method, that information just sits there.

Traditional Markov chain Monte Carlo sampling works by randomly exploring all possible explanations until it homes in on the ones that best fit the data. It scales badly with both the number of subhalos per image and the number of images per dataset. The neural likelihood-ratio approach pays a one-time training cost, then processes each new image in a fraction of a second. That’s the difference between a method that handles a handful of systems and one that operates at survey scale.

Beyond the immediate application, the work shows how simulation-based inference can be adapted to population-level problems across astrophysics, a template for similarly hard inference challenges throughout cosmology.

Bottom Line: By combining a carefully chosen dark matter observable with neural likelihood-ratio estimation, this method opens up population-level dark matter science from strong gravitational lensing, and it grows more powerful with every new lensing system future surveys discover.