A Double-Precision Map of Galaxy Clusters express
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.
Original abstract
as written by the authors arXiv:2606.30713
We compute the galaxy power spectrum at two-loop order in cosmological perturbation theory (effective field theory, EFT). We derive galaxy bias operators through the fifth order and obtain two-loop renormalization conditions for the their bias coefficients. We compute the two-loop integrals using a renormalization scheme consistent with the CLASS-PT code, allowing for an easy interface of our new computations with standard tools used in the one-loop galaxy power spectrum and bispectrum analyses. We also derive the relevant higher-derivative and stochastic contributions, and implement IR resummation using time-sliced perturbation theory. Having identified the redundant operators, we find that the two-loop galaxy power spectrum requires 21 additional EFT parameters per galaxy sample. We compare our computation with the galaxy-galaxy and galaxy-matter power spectra from the PT Challenge N-body simulation at $$z=0.61$$ and find a per mille-level agreement up to $$k=0.85~h$$Mpc$$^{-1}$$. We show that even with conservative priors on all EFT parameters, the two-loop model produces an unbiased measurement of the mass fluctuation amplitude $$\sigma_8$$ with three times narrower error-bars than the linear theory model. The improvement over the one-loop model is $$\simeq 40\%$$. This suggests significant gains in the two-loop EFT analyses of galaxy clustering and galaxy--lensing two-point functions ($$2\times2$$ pt) from CMB lensing maps and imaging surveys like Euclid, LSST, and Roman. In addition, our two-loop computation offers a probe of new physics scenarios that modify the shape of the matter power spectrum at wavenumbers $$(0.4-0.8)~h$$Mpc$$^{-1}$$ such as the presence of ultra-light axion dark matter sub-components with masses $$m_a\sim 10^{-24}$$ eV.
The two-loop calculation uses 21 additional parameters for each galaxy type — roughly the same number as the keys on a small musical instrument like a bandoneon. Each 'key' tunes the model to the real pattern of clustering.
In Carl Sagan's novel 'Contact,' the alien signal included a detailed map of the large-scale structure of the universe — two-loop calculations take a step toward the precision with which such messages could be decoded.