Cosmology faces a paradox: the standard model of cold dark matter made of heavy particles (WIMPs) works beautifully on cluster scales, but stumbles in galaxies—cores too dense, too many satellites. The alternative is fuzzy dark matter (FDM), composed of ultralight bosons with a mass around 10⁻²² eV that form a Bose-Einstein condensate. Their de Broglie wavelength—about a hundred parsecs—is comparable to the size of halos, so quantum superposition and wave-particle duality create density ripples everywhere. These ripples should leave a signature in gravitational lensing, but until now they have been modeled simplistically—as a random field. New computational power has made it possible, for the first time, to evolve the wave function “honestly” and predict lensed images.
The researchers ran 1000 independent numerical simulations of FDM halos with a mass of 4×10¹¹ M⊙, solving the Schrödinger equation and the Einstein field equations (in the Newtonian approximation) on a 512³ grid in a 40 kpc cube over 3.3 billion years. For each three-dimensional density cloud, they developed a method to find the optimal orientation: the halo’s major axis was rotated along 103 directions, the surface density was projected, and then 10⁴ source positions were scanned, minimizing the mismatch in the relative geometry of the four images using a Procrustes metric. Lensing was computed with the lenstronomy package, which uses gravity in the thin-lens approximation. Separately, they constructed 1000 realizations of a Gaussian random field based on the best-fit smoothed NFW profile—the standard approach for FDM lensing.
For a particle mass of 10⁻²² eV, the wave simulations yield a median position anomaly A (the rms deviation from the observed geometry averaged over four images) of 12 mas for the first radio source and 6 mas for the second. This is noticeably better than the predictions of the smooth NFW model (A≈12–18 mas) and the Gaussian random field approximation (A≈16–7 mas). In individual realizations, the anomaly drops below 3 mas—within the 3σ errors of radio interferometry. For lighter particles (10⁻²³ eV), the anomaly jumps sharply to ~50 mas, linked to a larger wavelength and bigger fluctuations. Crucially, the simulations were not tuned to match the observations: the halo evolved from a random initial state with no control over the final parameters.
Matching the data without any tuning is a strong argument in favor of FDM as an explanation for anomalies in lensed quasars, which serve as natural laboratories for gravitational lensing—a phenomenon driven by spacetime curvature. The work also shows that the Gaussian random field approximation overestimates the fluctuation amplitude in the critical region around the critical curve where images form. Therefore, for sensitive tests of the nature of dark matter, full wave modeling is necessary, not statistical surrogates.
Future high-angular-resolution observations (for example, with the upcoming SKA observatory or next-generation VLBI) will measure the positions of lensed images to better than 1 mas. Comparing these with a library of simulations for various FDM masses and halo profiles will allow a direct constraint on the dark matter particle mass. Moreover, the method will be extended to flux and time-delay anomalies—two more classic tests of strong gravitational lensing.
The results will impact cosmology, galaxy astrophysics, and particle physics by providing a new, independent method to probe the nature of dark matter that bypasses the limitations of direct detection.
The next steps are to include the baryonic component of the lensing galaxy, which can further perturb the potential, and to systematically analyze flux and time-delay anomalies within the wave-FDM framework.
This work touches on the fundamental problem of the nature of dark matter and the small-scale structure issues in ΛCDM—in particular, the missing satellites problem and the core-cusp problem. The success of FDM in reproducing lensed images without tuning adds weight to ultralight candidates and, at the same time, raises the question of the necessary precision of numerical methods in cosmology.
This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.
What the work does
This study presents direct numerical simulations of wave dark matter (FDM) with a full solution of the Schrödinger equation, and computes lensing on realistic halos without parameter tuning. The wave model is found to describe the observed anomalies in quasar image positions better than a smooth NFW profile or a Gaussian field approximation.
Where its strength is
The strength of this work lies in the first application of full wave simulations to specific gravitational lens systems with quantitative comparison to data, demonstrating the potential of FDM as a solution to the small-scale problems without parameter tuning.
Where it could go next
Worth trying
Why this may matter
The agreement with data without parameter tuning strengthens the case for wave dark matter as a physically motivated alternative, and the proposed lensing method could become a standard tool for testing dark matter models on galactic scales.
Nearby works in our archive
The closest work in the archive sits at 0.69 by meaning — the usual level for our corpus.
Quasar HS 0810+2554 is a unique target: ultra-precise radio measurements with e-MERLIN and VLBI revealed not one, but two compact radio sources, each of which is lensed into four images. That makes eight independent positions with an accuracy of ~1 mas.