https://doi.org/10.1140/epjc/s10052-026-15428-2
Regular Article - Theoretical Physics
The Milky Way and M31 rotation curves in Yukawa gravity: phenomenology and Bayesian analysis
1
Department of Physics, National University of Singapore, 117551, Singapore, Singapore
2
Research Center for Computing, National Research and Innovation Agency (BRIN), 40173, Bandung, Indonesia
3
Departemen Fisika, FMIPA, Universitas Indonesia, 16424, Depok, Indonesia
a
This email address is being protected from spambots. You need JavaScript enabled to view it.
Received:
7
August
2025
Accepted:
8
February
2026
Published online:
27
February
2026
Abstract
Yukawa gravity provides a generalized framework for modeling gravity modification. We investigate the rotation curve profiles of spiral galaxies under Yukawa-like theories governed by the coupling strength 
and the interaction range 
. We develop a unified analytical and numerical framework to calculate rotational velocities under Yukawa gravity, which includes contributions from all major galactic components: stellar bulge, disk, dark matter (DM) halo, and central supermassive black hole. The calculations show that and strongly influence velocity distributions by shifting peaks, creating double-peaked structures, or enhancing dark matter dominance in the bulge or disk. To assess observational implications, we perform Bayesian analyses using data from the Milky Way (MW) and Andromeda (M31), which offer complementary characteristics: MW provides precise velocity profiles across multiple scales, while M31 includes broader morphological constraints. We examine four scenarios: Yukawa gravity without dark matter, dark matter with non-trivial coupling, fully modified gravity, and standard Newtonian gravity. Results show that MW models with 
kpc yield high Bayes factors but risk overfitting, as dark matter mimics baryonic kinematics, while M31’s photometric priors from conjugate observations mitigate this, yielding robust parameter estimates. However, in M31, Bayes factors favor Newtonian gravity, suggesting that current data lack the precision to resolve more complex models. This finding highlights two key needs: (i) realistic, physically or empirically informed priors to avoid biased constraints, and (ii) high-precision data with independent photometry to guard against overfitting. Our framework offers a scalable approach for testing gravity with large galactic rotation curve datasets.
© The Author(s) 2026
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Funded by SCOAP3.
