pith. sign in

arxiv: 2604.24855 · v1 · submitted 2026-04-27 · 🌌 astro-ph.GA

Multi-component, axisymmetric dynamical models of dSphs based on distribution functions: inferences on dark matter and intermediate-mass black holes in Draco and Ursa Minor

R. Pascale , G. Battaglia , J.M. Arroyo-Polonio , E. Vasiliev , C. Nipoti , G.F. Thomas This is my paper

Pith reviewed 2026-05-08 02:17 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords dwarf spheroidal galaxiesdark matter density profilesaxisymmetric dynamical modelsDracoUrsa Minordistribution functionsintermediate-mass black holes
0
0 comments X

The pith

Axisymmetric models indicate a cuspy dark matter profile in Draco and a cored profile in Ursa Minor, with no evidence for intermediate-mass black holes.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops the first multi-component axisymmetric dynamical models of dwarf spheroidal galaxies using distribution functions. It applies the models to Draco and Ursa Minor by treating the stars as two chemically and kinematically distinct populations that trace a spherical dark matter potential, which may also contain a central black hole. Fits to Gaia astrometry and spectroscopic velocities and metallicities show that Draco hosts a centrally cusped dark matter density while Ursa Minor is better described by a shallower core. The same models place strict upper limits on any central black hole masses and demonstrate that spherical approximations tend to make dark matter profiles appear more cuspy than axisymmetric ones do. A reader would care because these two galaxies supply some of the cleanest local tests of dark matter behavior on small scales.

Core claim

The authors construct axisymmetric, multi-population distribution-function models in which chemo-dynamically distinct stellar components orbit in a spherical potential generated by a dominant dark matter halo plus an optional central intermediate-mass black hole. When fitted to discrete Gaia proper motions, line-of-sight velocities and metallicities, the models recover an inner dark matter density slope of 0.98 in Draco and 0.37 in Ursa Minor. Both galaxies are better described by two stellar populations than by one, and the data yield no support for black holes, with 95-percent upper mass limits of 10^5.2 solar masses in Draco and 10^3.33 solar masses in Ursa Minor.

What carries the argument

Axisymmetric distribution functions for two chemo-dynamically distinct stellar populations embedded in a spherical dark matter potential that may include a central black hole.

Load-bearing premise

The chemo-dynamically distinct stellar populations are assumed to be in dynamical equilibrium inside the spherical gravitational potential of the dark matter halo and any central black hole.

What would settle it

A measurement of the central stellar density slope or velocity dispersion profile in Draco that is inconsistent with an inner logarithmic slope near 1 would falsify the reported dark matter cusp.

Figures

Figures reproduced from arXiv: 2604.24855 by C. Nipoti, E. Vasiliev, G. Battaglia, G.F. Thomas, J.M. Arroyo-Polonio, R. Pascale.

Figure 1
Figure 1. Figure 1: Marginalised posterior distributions of the metallicity-related model parameters. The left set of panels refers to Draco, the view at source ↗
Figure 2
Figure 2. Figure 2: Structural properties of Draco (top) and Ursa Minor (bottom) inferred from the fit to the G26 sample. Left, middle and right view at source ↗
Figure 3
Figure 3. Figure 3: Kinematic properties of Draco inferred from the fit to the G26 sample. The top three rows show projected velocity-dispersion view at source ↗
Figure 4
Figure 4. Figure 4: Distribution of stars in the metallicity–line-of-sight velocity plane for Draco (top panels) and Ursa Minor (bottom panels). view at source ↗
Figure 5
Figure 5. Figure 5: Same as Fig view at source ↗
Figure 6
Figure 6. Figure 6: Difference between mean line-of-sight and systemic velocities (∆v) as a function of the azimuthal angle, ϕ, shown for the total (purple, top), MR (red, middle), and MP (blue, bottom) populations of Draco. Left and right columns correspond to the G26 and W23 datasets, respectively. Uncertainties on the mean values, which we show with bands, are derived from the Fisher information matrix. The vertical dashed… view at source ↗
Figure 7
Figure 7. Figure 7: Same as Fig view at source ↗
Figure 8
Figure 8. Figure 8: DM halo properties. Left panel: median DM halo density profile (solid curves) from the chemo-dynamical models of Sec view at source ↗
Figure 9
Figure 9. Figure 9: Marginalised one- and two-dimensional posterior distributions of the DM halo parameters of Draco (left) and Ursa Minor view at source ↗
read the original abstract

Dwarf spheroidal galaxies (dSphs) are prime laboratories for studying dark matter (DM) and the black hole demographics in the low-mass regime. These systems are also often flattened; nevertheless most studies rely on spherical models, potentially affecting dynamical inferences. We introduce the first multi-component, axisymmetric dynamical models of dSphs based on distribution functions and apply them to the Milky Way dSphs Draco and Ursa Minor. The stellar distribution is described by chemo-dynamically distinct axisymmetric populations tracing a spherical potential generated by a dominant DM halo and a central intermediate-mass BH (IMBH). The models are fitted to discrete stellar data from a Gaia-based astrometric sample and two spectroscopic datasets providing line-of-sight velocities and metallicities, testing robustness across samples. We compare the DM properties under different modelling assumptions, including flattened one- and spherical two-component models. Both galaxies are better described by two stellar populations: a metal-rich, kinematically colder and concentrated component, and a more extended metal-poor one with hotter kinematics. We detect weak rotation, dynamically unimportant and ignored in the models. We measure a cuspy DM density profile in Draco ($\gamma=0.98_{-0.26}^{+0.28}$), and a more cored distribution ($\gamma=0.37_{-0.24}^{+0.31}$) for Ursa Minor. The DM halo of Draco remains stable across all models and datasets, making it the most robustly determined in the Local Group and highly relevant for indirect DM searches. We show that modelling flattened systems with spherical models can bias the DM inner slope towards cuspier values, while we find no degeneracy between inner halo density and inclination. We find no evidence for IMBHs and place upper limits on their masses, $\log M_{\rm BH}[M_{\odot}] < 5.2$ for Draco and $< 3.33$ for Ursa Minor (95% confidence).

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 3 minor

Summary. The paper introduces the first multi-component axisymmetric distribution-function (DF) models for dwarf spheroidal galaxies and applies them to Draco and Ursa Minor. Chemo-dynamically distinct axisymmetric stellar populations are modeled as tracers in a spherical potential generated by a dominant DM halo (with free inner logarithmic slope γ) plus an optional central IMBH. The models are fitted to discrete Gaia astrometric data plus two spectroscopic datasets, with robustness tests across samples and comparisons to flattened one-component and spherical two-component models. The main results are a cuspy DM profile in Draco (γ = 0.98_{-0.26}^{+0.28}) that remains stable across all variants, a more cored profile in Ursa Minor (γ = 0.37_{-0.24}^{+0.31}), weak rotation that is dynamically unimportant, and no evidence for IMBHs with 95% upper limits log M_BH < 5.2 (Draco) and < 3.33 (Ursa Minor).

Significance. If the central results hold, the work supplies some of the most robust DM inner-slope constraints available for any Local Group dwarf, particularly the stable cuspy profile in Draco that is directly relevant to indirect DM searches. The methodological advance of axisymmetric DF models for flattened dSphs, together with the explicit demonstration that spherical models bias γ cuspier, constitutes a clear step forward. Credit is due for the systematic robustness tests across independent datasets and the direct comparison of modeling assumptions (flattened one-component vs. spherical two-component).

major comments (1)
  1. [Abstract and §3] Abstract and §3 (model description): The stellar DFs are axisymmetric while the total gravitational potential is forced to be strictly spherical (DM halo + optional IMBH). Although the manuscript states that spherical models bias the recovered inner slope γ toward cuspier values, no quantitative test of the residual bias introduced by retaining a spherical potential inside an otherwise axisymmetric framework is presented, nor is an orbit-integration or non-spherical-potential check described. Because the headline γ values (0.98 for Draco, 0.37 for Ursa Minor) and their uncertainties are the central claim, this assumption is load-bearing and requires either a direct quantification of the systematic or an explicit justification that the bias is negligible compared with the quoted errors.
minor comments (3)
  1. [§4.2 and Table 2] §4.2 and Table 2: The parameter table lists the DM scale radius and density normalization but does not explicitly show the derived γ posterior for each dataset/model variant; adding a dedicated row or supplementary table would improve traceability of the robustness claim.
  2. [Figure 5] Figure 5 caption: The inclination angle is stated to be unconstrained, yet the text claims “no degeneracy between inner halo density and inclination.” A brief quantitative statement (e.g., correlation coefficient or marginalised posterior width) would clarify this point.
  3. [Introduction] References: The comparison to earlier spherical Jeans or Schwarzschild models of Draco and Ursa Minor would benefit from citing the specific works whose γ values are being contrasted (currently referenced only generically in the introduction).

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their positive and constructive report, which recognizes the methodological advance and the robustness of the DM slope constraints. We address the single major comment below.

read point-by-point responses

  1. Referee: [Abstract and §3] Abstract and §3 (model description): The stellar DFs are axisymmetric while the total gravitational potential is forced to be strictly spherical (DM halo + optional IMBH). Although the manuscript states that spherical models bias the recovered inner slope γ toward cuspier values, no quantitative test of the residual bias introduced by retaining a spherical potential inside an otherwise axisymmetric framework is presented, nor is an orbit-integration or non-spherical-potential check described. Because the headline γ values (0.98 for Draco, 0.37 for Ursa Minor) and their uncertainties are the central claim, this assumption is load-bearing and requires either a direct quantification of the systematic or an explicit justification that the bias is negligible compared with the quoted errors.

    Authors: We agree that the spherical-potential assumption within an otherwise axisymmetric framework is an important approximation whose residual effect on γ merits explicit discussion. The manuscript already demonstrates that fully spherical (DF + potential) models bias γ cuspier relative to our axisymmetric-DF models; however, we did not perform a direct test with a non-spherical potential. A full axisymmetric or triaxial potential would require a substantial extension of the DF machinery and orbit library, which lies outside the scope of the present study. In the revised manuscript we add an explicit justification in §3: the DM halo dominates the potential by more than an order of magnitude, so any flattening induced by the stellar component is expected to be negligible; the observed stellar axis ratios (~0.7) and the lack of detectable rotation further support a nearly spherical total potential. We also state that the quoted uncertainties on γ are statistical only and do not formally include this systematic, although we argue it is sub-dominant to the reported errors. We are prepared to expand this discussion or add a caveat if the referee requests specific wording. revision: partial

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper constructs multi-component axisymmetric distribution-function models and fits their parameters (including DM inner slope γ and optional IMBH mass) directly to independent Gaia astrometric plus spectroscopic datasets for Draco and Ursa Minor. The central inferences on DM profiles emerge from likelihood maximization against these external observations; comparisons across spherical vs. flattened assumptions are performed as robustness tests rather than as self-referential inputs. No quoted step reduces a claimed prediction to a fitted quantity by construction, nor does any load-bearing premise rest on an unverified self-citation chain. The modeling framework therefore remains independent of its own outputs.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard assumptions in galactic dynamics and fitted parameters for the DM profile and BH mass.

free parameters (2)
  • DM inner slope gamma = 0.98 for Draco, 0.37 for Ursa Minor
    Fitted parameter for the inner density profile of the dark matter halo.
  • IMBH mass upper limit = <5.2 log solar masses for Draco
    Tested parameter with no detection, bounds placed from data fit.
axioms (2)
  • domain assumption The gravitational potential is spherical despite axisymmetric stellar distributions.
    Stated in the model description for the DM halo and IMBH.
  • domain assumption Stellar populations are in steady-state equilibrium.
    Fundamental for distribution function modeling of discrete stellar data.

pith-pipeline@v0.9.0 · 5700 in / 1363 out tokens · 69765 ms · 2026-05-08T02:17:24.618130+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

7 extracted references · 7 canonical work pages

  1. [1]

    Abe, S., Abhir, J., Abhishek, A., et al. 2024, J. Cosmology Astropart. Phys., 2024, 047 Aditya, K. & Mangalam, A. 2025, arXiv e-prints, arXiv:2512.14146 Akins, H. B., Casey, C. M., Lambrides, E., et al. 2025, ApJ, 991, 37 Albert, A., Alfaro, R., Alvarez, C., et al. 2020, Phys. Rev. D, 101, 103001 Aparicio, A., Carrera, R., & Martínez-Delgado, D. 2001, AJ,...

    work page arXiv 2024

  2. [2]

    Each MCMC run is post-processed to ensure that the retained samples are independent and representative of the converged posterior

    and the snooker update (ter Braak & Vrugt 2008), which provide robust performance across a wide range of posterior shapes. Each MCMC run is post-processed to ensure that the retained samples are independent and representative of the converged posterior. In particular: i) we remove an initial burn-in phase long enough to ensure model convergence; ii) we pe...

    work page 2008

  3. [3]

    2018; Battaglia et al

    We retained only stars with a full five-parameter as- trometric solution and high-quality astrometry, requiring astrometric_params_solved>31 and a renormalised unit weight errorruwe<1.4 (Lindegren et al. 2018; Battaglia et al. 2022). Also, to ensure reliable photometric measurements, we kept stars with non-nullG BP −G RP colour, and applied a quality filt...

    work page 2018

  4. [4]

    2019): G0 =G−2.664E(B−V), (G−G RP)0 =(G−G RP)−0.643E(B−V).(B.3) We do not consider errors on the estimates of theE(B−V), thus ∆(G−G RP)0 = ∆(G−G RP) and∆G 0 = ∆G

    and apply the following formulae to compute the de-reddenedG 0 band magnitude and (G−G RP)0 colour (see Sestito et al. 2019): G0 =G−2.664E(B−V), (G−G RP)0 =(G−G RP)−0.643E(B−V).(B.3) We do not consider errors on the estimates of theE(B−V), thus ∆(G−G RP)0 = ∆(G−G RP) and∆G 0 = ∆G. Finally, we limited the sample to stars withG 0-band mag- nitudes brighter ...

    work page 2019

  5. [5]

    (2015) and W23 catalogues, who pro- vide multi-epoch, medium-resolution spectra for thousands of targets observed with the M2FS and Hectochelle spectrographs

    The first chemo-kinematic dataset adopted here is constructed from the Walker et al. (2015) and W23 catalogues, who pro- vide multi-epoch, medium-resolution spectra for thousands of targets observed with the M2FS and Hectochelle spectrographs. These data include line-of-sight velocities, atmospheric param- eters, chemical abundances, and signal-to-noise (...

    work page 2015

  6. [6]

    Together, these six criteria define a baseline set of stable, well-behaved spectra suitable for identifying Draco and Ursa Minor kinematic members

    are excluded to avoid binary-induced mo- tions14. Together, these six criteria define a baseline set of stable, well-behaved spectra suitable for identifying Draco and Ursa Minor kinematic members. After establishing these base samples, we apply additional filters targeting the quality of the velocity and metallicity de- terminations. For line-of-sight ve...

    work page 2026

  7. [7]

    The quantity Reffdenotes the effective (or half-counts) radius, defined as the elliptical radius enclosing half of the total star counts

    represents a possible offset of the system’s centre with respect to the reference value, and the ellipticity is defined ase≡1−b/a, whereaandbare the semi-major and semi- minor axes of the isodensity contours, respectively. The quantity Reffdenotes the effective (or half-counts) radius, defined as the elliptical radius enclosing half of the total star coun...

    work page 2020