2026/05/31 by Prasun Panthi, Md Shahrier Islam Arham
Physics and Astronomy · #astro-ph.CO #astro-ph.GA #astro-ph.IM
21 pages, 7 figures
arxiv created 2026/08/05 · arxiv updated 2026/08/06
Ultralight-dark-matter rotation-curve analyses often infer the particle mass after fixing a relation between the central soliton and its host halo. The resulting mass constraint is then conditional on a population-level relation that the data may not independently support. We present a reusable hierarchical Bayesian framework that instead infers the ultralight particle mass and the core-halo scaling exponent jointly from a galaxy population. The differentiable forward model combines a Schive-normalized soliton with a smoothly matched, regularized NFW envelope and allows galaxy-level halo and stellar parameters to be inferred together with the global dark-matter parameters. We apply the framework to 106 SPARC galaxies, including 26 systems with bulges, and sample the resulting 346-dimensional posterior with JAX/NumPyro NUTS. The free-scaling run has zero divergences and R≃ 1.000 for the global parameters. The posterior moves to the high-mass, weak-scaling boundary, with log10(mϕ/eV)=-19.20+0.12-0.11 and α=0.014+0.023-0.011. In this regime, the solitonic cores lie below the radial scales probed by the rotation curves, while the baryonic terms and outer NFW envelope carry the visible fits. The same boundary behaviour remains after removing UGC06787 and after widening the high-mass stellar-to-halo-mass prior. The selected SPARC sample therefore does not give an interior joint constraint on the particle mass and core-halo scaling relation within the adopted model. The framework makes this unresolved-core limit explicit instead of interpreting the prior boundary as an interior mass constraint. Its modular structure also allows the same analysis to be applied to synthetic recovery tests, expanded or higher-resolution rotation-curve samples, and alternative halo models.