Counting Particles in Cosmic Halos: New Map Ties Dark Matter’s Mass to Galaxy Size

What if the invisible scaffolding of the universe is not a swarm of particles at all, but a single, gigantic quantum wave? That is the provocative premise behind a growing body of research into ultralight dark matter, and a new theoretical study has now mapped, with unprecedented precision, exactly when such cosmic condensates can hold themselves together. The work, published in The European Physical Journal C, provides one of the most systematic characterizations to date of the equilibrium configurations that self-gravitating Bose–Einstein condensates can form, and it delivers a sobering verdict on one of the most popular dark matter candidates.

The study, conducted by Francisco A. Guzmán, Elías Castellanos, and Jorge Mastache, tackles the Gross–Pitaevskii–Poisson (GPP) system, a pair of coupled equations that describes a Bose–Einstein condensate evolving under its own gravity. In this picture, the entire dark matter halo of a galaxy is represented by a single macroscopic wavefunction whose squared amplitude gives the mass density. The first equation is a nonlinear Schrödinger equation that includes the gravitational potential and a short-range repulsive interaction between the bosons, while the second is the Newtonian Poisson equation, which computes the gravitational field generated by the condensate’s own density. Solving these equations simultaneously is notoriously difficult because the problem is nonlinear and nonlocal: the gravitational potential depends on the wavefunction everywhere, not just at a single point.

What sets this analysis apart from much of the existing literature is a deliberate choice of bookkeeping. Many previous studies exploit mathematical scaling symmetries of the equations, solving everything in dimensionless units and only converting back to physical quantities at the end. The Mexican team instead kept three quantities explicit throughout: the boson mass, the scattering length that characterizes the strength of repulsive self-interaction between particles, and the total particle number N. That last parameter is crucial, because it directly fixes the total halo mass through the simple relation that the halo mass equals the boson mass multiplied by the particle count. By treating N as an independent control parameter, the researchers built a direct bridge between microscopic particle physics and macroscopic observables such as halo size, density profile, and rotation speed.

Technically, the team discretized the radial equations on a finite grid extending to 100 kiloparsecs and solved the resulting nonlinear eigenvalue problem with a Newton–Raphson iterative scheme, updating the wavefunction, the gravitational potential, and the chemical potential self-consistently until convergence to tolerances between 10 and 8 and 10 to the minus 10. To guard against artifacts of the starting guess, they initialized the solver with four different smooth, localized profiles: Gaussian, exponential, linear-exponential, and hyperbolic-secant shapes. A striking robustness emerged. Whenever different starting profiles converged to the same solution branch, the resulting stationary configurations were numerically indistinguishable, with relative differences below roughly one part in ten thousand for the wavefunction. The choice of initial guess affected only how often the iteration converged, not the physics of the answer.

Scanning a 40-by-40 grid of boson masses between 10 to the minus 23 and 10 to the minus 20 electronvolts and particle numbers between 10 to the 93 and 10 to the 99, the researchers classified every converged solution into three families. Ground states are nodeless, meaning the wavefunction never crosses zero, and they represent the lowest-energy equilibrium of the condensate. Excited states possess one or more radial nodes, producing oscillatory density shells around a central core, much like the excited orbitals of an atom. Unbound configurations fail a binding criterion, meaning the energy per particle exceeds the depth of the self-generated gravitational well, so no localized condensate can exist. The resulting phase diagrams show the ground-state branch occupying a well-defined diagonal band in the mass–particle-number plane: heavier bosons need fewer particles to bind, while lighter ones require vastly larger numbers.

One of the most consequential findings concerns how repulsive self-interactions reshape this map. As the scattering length increases, the ground-state branch systematically migrates toward larger boson masses and smaller particle numbers, because the repulsive pressure helps support the condensate against gravitational collapse. The excited-state and unbound regions, by contrast, shift far more slowly. This sensitivity matters for phenomenology: the location of the ground-state band, the region where realistic galactic halos live, is exquisitely dependent on the strength of the self-interaction, while the rest of the solution space is comparatively stubborn.

From the converged solutions, the team extracted empirical scaling laws for the characteristic halo radius, defined as the sphere enclosing 99 percent of the mass. In the weakly interacting regime, the fitted relation gives a radius proportional to the boson mass raised to the power of minus 2.996 and the particle number to the power of minus 0.998, values astonishingly close to the analytic prediction of minus 3 and minus 1 for non-interacting solitons. Rewritten in terms of halo mass, this recovers the famous mass–radius relation of Schrödinger–Poisson solitons, in which the radius scales inversely with the boson mass squared and the halo mass, providing a clean consistency check of the numerics. When self-interactions dominate, the scaling shifts toward the Thomas–Fermi expectation, with the radius growing as the square root of the scattering length and falling as the boson mass to the power of minus 1.545. Yet the fitted exponents reveal that the solutions occupy an intermediate regime where quantum pressure, self-interaction pressure, and gravity all contribute simultaneously, deviating measurably from the asymptotic Thomas–Fermi limit. A three-parameter fit confirms a clear hierarchy: the boson mass exerts the strongest influence on halo size, the particle number a moderate one, and the scattering length a softer but non-negligible effect.

The astrophysical payoff comes when these equilibrium configurations are confronted with real galaxies. Using ground-state solutions with a boson mass of 10 to the minus 22 electronvolts, the team modeled the rotation curves of two dwarf galaxies from the SPARC catalog, KK98-251 and UGC01281, adding the baryonic contribution from gas and stars to the condensate’s contribution computed from the enclosed mass profile. The results were impressive. For KK98-251, a halo of about 896 million solar masses with a characteristic radius of 5.05 kiloparsecs reproduced the observed kinematics, while UGC01281 required roughly 3.58 billion solar masses spread over 5.88 kiloparsecs. Remarkably, in both cases the solitonic core alone accounted for the dark matter, with no need for the extended interference-supported envelope that cosmological simulations of fuzzy dark matter typically predict around the central condensate.

Here, however, the story collides with cosmological constraints. Observations of the Lyman-alpha forest, the thicket of absorption features imprinted by intergalactic hydrogen on the light of distant quasars, place stringent lower bounds on the boson mass when ultralight bosons constitute all of the dark matter. Depending on the dataset and modeling assumptions, masses below roughly 10 to the minus 21 up to 2 times 10 to the minus 20 electronvolts are disfavored, excluding the canonical fuzzy dark matter scale of 10 to the minus 22 electronvolts. The new study shows that repulsive self-interactions can indeed rescue the equilibrium solutions: for scattering lengths around 10 to the minus 69 meters, ground-state branches extend into the mass range of roughly 2 times 10 to the minus 20 to 10 to the minus 18 electronvolts, comfortably compatible with the Lyman-alpha bounds. But there is a catch. Those high-mass, self-interacting configurations generate circular velocities far too small to match the dwarf-galaxy rotation curves that the lighter bosons fit so well. Within the class of stationary halos examined, the tension between the masses preferred by galactic dynamics and those demanded by the Lyman-alpha forest remains unresolved.

The broader significance of the work lies in its framework rather than any single verdict. By keeping the particle number explicit and mapping the full solution space, the authors have produced phase diagrams and scaling relations that connect fundamental boson properties directly to observable galactic quantities, offering a controlled laboratory for testing Bose–Einstein condensate dark matter models. The authors note that future work should probe the dynamical stability of the ground-state and excited branches, extend the analysis to halos with outer envelopes, and confront a much larger sample of rotation curves. For now, the message is twofold: quantum-wave halos remain a viable and elegant description of dwarf galaxies at the canonical mass scale, but escaping the tightening noose of Lyman-alpha constraints through self-interactions may cost the model the very galaxies it explains best. The universe, it seems, is still auditing every proposal for its missing matter.

Subject of Research: Equilibrium solutions of the Gross–Pitaevskii–Poisson system as a model for self-gravitating Bose–Einstein condensate dark matter halos

Article Title: Equilibrium halo solutions of the Gross–Pitaevskii–Poisson system: the role of the particle number

Article References: Guzmán, F. A., Castellanos, E., & Mastache, J. (2026). Equilibrium halo solutions of the Gross–Pitaevskii–Poisson system: the role of the particle number. The European Physical Journal C, 86(9), Article 1103. https://doi.org/10.1140/epjc/s10052-026-16373-w

Image Credits: AI Generated

DOI: 10.1140/epjc/s10052-026-16373-w

Keywords: dark matter, Bose–Einstein condensate, Gross–Pitaevskii–Poisson system, ultralight bosons, fuzzy dark matter, galaxy rotation curves, solitonic core, Thomas–Fermi limit, Lyman-alpha forest, self-interaction, halo scaling relations, computational astrophysics