Bayesian analysis of Gaia DR3 reveals the Milky Way’s dark matter profile

The Milky Way may be finally surrendering one of its best-kept secrets. In a new Bayesian model comparison analysis published in Astrophysics and Space Science, Aryan Singh and Shantanu Desai of the Department of Physics at IIT Hyderabad have pitted seven different dark matter halo models against each other, along with three variants of modified Newtonian dynamics, to determine which mathematical description best matches the rotation curve of our Galaxy as measured by the European Space Agency’s Gaia mission. Their verdict, drawn from four independent compilations of Gaia DR3-based rotation curve data, is strikingly clear: the Milky Way’s dark matter halo is best described by the Einasto profile, a smooth, mathematically elegant density law that outperforms the widely used Navarro–Frenk–White (NFW) profile across most combinations of datasets and baryonic mass models.

The question of how dark matter is distributed in the Galactic halo is far more than an exercise in curve-fitting. The shape of the halo determines the local dark matter density, which in turn sets the expected signal rate for direct detection experiments searching for weakly interacting massive particles deep underground. It also governs how astronomers model the orbits of stellar streams, the dynamics of satellite galaxies, and the total mass of the Milky Way out to its farthest reaches. For decades, simulations of cold dark matter have suggested that halos should follow the NFW profile, whose density diverges as a power law toward the galactic center, producing a so-called “cuspy” core. Yet observations of many galaxies, particularly dwarf spirals, have long hinted that their dark matter is more centrally diffuse, or “cored,” a tension that ranks among the most persistent small-scale challenges to the standard cosmological paradigm.

Singh and Desai approached the problem with the full machinery of modern Bayesian statistics. Rather than simply asking which model fits the data best in a least-squares sense, they computed the Bayesian evidence for each model, a quantity that balances goodness of fit against model complexity and naturally penalizes models that add parameters without earning them. The computations were performed using DYNESTY, a dynamic nested sampling package that efficiently estimates both posterior distributions and evidences for high-dimensional models. The team combined three different baryonic mass models, describing the visible matter in the Galactic bulge, disk, and gas, with seven dark matter halo profiles, including the NFW profile, the Einasto profile, the Burkert profile, the Plummer profile, and several cored variants such as the pseudo-isothermal halo and the core-modified profile of Lazar and collaborators.

On the data side, the authors exploited four recent rotation curve compilations built on Gaia Data Release 3, the third and most precise catalog of positions, distances, and proper motions for nearly two billion stars delivered by the Gaia spacecraft. These include the circular velocity curve of the Milky Way derived from luminous red giant branch stars by Zhou and colleagues, the kinematic mapping of the Galactic disk to roughly 30 kiloparsecs by Wang and collaborators, and related datasets that trace the Galaxy’s rotation from about 5 to 25 kiloparsecs from the Galactic center. Each dataset encodes the same fundamental observable: the orbital speed of stars and gas as a function of distance from the Galactic center. Because the visible matter alone cannot explain why orbital speeds remain roughly flat, or in some recent analyses even decline in Keplerian fashion, at large radii, the shape of the rotation curve at those distances is a direct probe of the invisible halo.

The results carry implications that ripple across dark matter physics. In nearly every dataset-baryon combination the authors tested, the Einasto profile emerged as the preferred phenomenological description, beating the NFW profile by decisive margins in the Bayesian evidence. The Einasto profile, first introduced by Jaan Einasto in 1965 in a completely different context as a model for galactic light distributions, describes the halo density as a slowly varying power of the radius, falling more gradually near the center than NFW’s sharp cusp. Its victory here suggests that the Milky Way’s inner halo is smoother and less centrally concentrated than the classic cold dark matter prediction, at least as probed by Gaia-era kinematics.

Equally significant was the performance of the cored profiles. Models in which the dark matter density flattens to a finite central value, such as the Burkert profile originally proposed to explain dwarf galaxy halos, were systematically preferred over the cuspy NFW model. This finding aligns the Milky Way with a broader pattern seen in external galaxies and lends indirect support to scenarios in which baryonic feedback, the energetic outflows from star formation that push gas and dark matter around, or even self-interacting dark matter, reshapes the inner halo. The authors caution, however, that their conclusions apply within the adopted modeling framework and the specific Gaia-based rotation curve datasets used; the comparison is phenomenological rather than a falsification of any particular particle physics model.

The study also delivered a pointed verdict on modified gravity. Modified Newtonian Dynamics, or MOND, proposed by Mordehai Milgrom in 1983, dispenses with particle dark matter altogether by altering the law of gravity at accelerations below a characteristic scale of roughly 10^-10 meters per second squared. Singh and Desai tested MOND using three different interpolating functions, the mathematical bridges that connect the Newtonian and deep-MOND regimes, coupled to the same baryonic models. Within the implementations they considered, all three MOND variants provided poorer fits to the Gaia-based Milky Way rotation curve than the best dark matter profiles, and the Bayesian evidence firmly ranked them below the Einasto and cored halo models. While the authors are careful to note that this does not exhaust every possible MOND formulation, particularly fully relativistic extensions, it adds the Milky Way’s own kinematics to the list of arenas where modified gravity struggles to compete on equal statistical footing.

One perhaps counterintuitive outcome of the analysis is what did not matter: the baryonic models. The team’s three descriptions of the Galaxy’s luminous content, differing in how they treat the stellar disk, thick disk, and gas distribution, produced no decisive winner. None of the baryonic models was consistently favored over the others across the dark matter profiles and datasets. This suggests that, at the level of precision currently offered by the rotation curve data, systematic uncertainties in the baryonic component are not the dominant driver of model preference; the dark matter profile itself is doing the discriminating work. It also underscores how much room remains for improvement, since future data with tighter errors at large Galactic radii could well sharpen the sensitivity to the visible mass distribution.

The timing of this work is notable. Recent analyses of Gaia DR3 data, including the detection of a Keplerian decline in the Milky Way rotation curve beyond the solar circle by Jiao and collaborators, have reignited debate about the Galaxy’s total mass and halo structure, with some authors arguing that a lighter, faster-declining Milky Way carries cosmological consequences. Assessing the robustness of rotation curves inferred through the Jeans equations against cosmological simulations, as Koop and colleagues have done, remains an active concern, since systematics in distance estimates and stellar selection can bias the inferred velocities. By folding multiple independent rotation curve compilations into a single Bayesian comparison framework, Singh and Desai have provided a way to see which conclusions survive the choice of dataset, a robustness check that single-dataset analyses cannot offer.

The broader stakes extend beyond Galactic astronomy. Dark matter remains one of the most profound unsolved problems in physics, evidenced across galaxy rotation curves, gravitational lensing, and the cosmic microwave background, yet still escaping direct detection after decades of increasingly sensitive experiments. Pinning down the precise shape of the Milky Way’s halo matters for interpreting those experiments, because the local density and velocity distribution of dark matter set the expected event rates in detectors. A smoother, cored, or Einasto-like halo changes the translation between detector limits and particle properties. As Gaia continues to refine our map of the Galaxy and future surveys extend rotation measurements to ever larger radii, the Bayesian machinery demonstrated in this study offers a template: rather than assuming a halo model, let the data adjudicate among them. For now, the Milky Way’s dark matter appears to prefer the gentle mathematics of Einasto over the sharp cusp of NFW, and to leave little statistical room for gravity itself to do the dark matter’s work.