A cosmological map is like a snapshot of ripples on water: from it you can reconstruct what stone was thrown. The stone is the properties of the universe — how much ordinary matter it contains, how much dark matter, how fast it is expanding. The ripples are the pattern of galaxies that telescopes measure. Until now, calculations for such maps stopped at the first loop of corrections, in so-called one-loop models. The new work pushes them to the second loop, the two-loop level, and for the first time provides a complete set of formulas for comparison with observations.

The authors derived how galaxy brightness relates to the invisible mass distribution, including effects up to fifth order in perturbation theory. The result is a set of 21 additional parameters for each galaxy type: these numbers act like corrections to camera lenses, removing distortions. The calculation is integrated with the popular CLASS-PT code so that other groups can use it. A test on a computer simulation of the universe showed agreement at the level of tenths of a percent down to scales of about 7.4 million light-years.

The practical payoff is direct: the Euclid telescope, the Vera Rubin Observatory, and the Roman Space Telescope are collecting data right now. A two-loop analysis will extract more information from them, especially at the intersection with cosmic microwave background maps. Furthermore, the method is sensitive to the shape of the power spectrum on scales of 7.4–15 million light-years. This window is important for the search for ultralight axion dark matter particles with masses around 10⁻²⁴ eV — candidates that leave a characteristic dip in the galaxy distribution.

The ripples on the water are no longer blurry. The second loop of corrections turns a coarse-grained photograph of the universe into a snapshot where every fold is visible.