pith. sign in

arxiv: 2603.00257 · v2 · submitted 2026-02-27 · 🌌 astro-ph.GA

A dynamical attractor in the evolution of dwarf spheroidal galaxies

Jorge Pe\~narrubia , Ethan O. Nadler This is my paper

Pith reviewed 2026-05-15 18:07 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords dwarf spheroidal galaxiesdynamical evolutiondark matter haloestidal strippingstellar orbitssubhaloesMilky Way satelliteshalo mass relation
0
0 comments X

The pith

Stellar orbits in dwarf spheroidal galaxies irreversibly expand toward a dynamical attractor with half-light radius matching the dark halo peak radius.

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

Controlled N-body experiments demonstrate that stellar orbits in dSphs embedded in haloes with many subhaloes gain energy from stochastic force fluctuations and expand irreversibly. They approach an attractor state where the stellar half-light radius equals the halo's radius of maximum circular velocity and the velocity dispersion is about half that peak velocity. This evolution occurs in isolation but is sped up by tidal stripping. Assuming Milky Way dSphs have reached this state, their inferred dark halo masses form tight sequences when plotted against half-light radius, with smaller galaxies following tracks expected for cuspy haloes and larger ones suggesting cored profiles. The results indicate that the observed variety in dSph structures arises mainly from this evolutionary process rather than initial star formation conditions.

Core claim

We use controlled N-body experiments to study the dynamical evolution of dwarf spheroidal galaxies embedded in dark-matter haloes containing a large population of dark subhaloes. Stellar orbits subject to stochastic force fluctuations irreversibly gain energy and expand toward a dynamical attractor characterized by a stellar half-light radius r_half ≈ r_max and a velocity dispersion σ ≈ 0.5 v_max, where v_max is the peak circular velocity of the host halo at radius r_max. This state is reached both in isolation and under tidal stripping, although tidal mass loss significantly accelerates the evolution. Assuming that the Milky Way dSphs have reached this state, the inferred halo masses of MW

What carries the argument

The dynamical attractor state reached when stellar half-light radius equals the dark halo's r_max and velocity dispersion is half v_max, driven by stochastic heating from subhalo encounters.

Load-bearing premise

That the Milky Way dwarf spheroidal galaxies have already reached the proposed dynamical attractor state.

What would settle it

If measurements show that the inferred dark matter halo masses of Milky Way dSphs do not form narrow sequences when plotted against their half-light radii, or if isolated dSphs are not found to be larger than satellites.

Figures

Figures reproduced from arXiv: 2603.00257 by Ethan O. Nadler, Jorge Pe\~narrubia.

Figure 1
Figure 1. Figure 1: Left panel: Stellar velocity dispersion 𝜎 as a function of half-light radius 𝑟half, derived from the virial theorem for stellar tracers with Plummer (solid) or Hernquist (dashed) profiles embedded in Dehnen (colored) and NFW (black) haloes with 𝑀ℎ = 109 𝑀⊙ and 𝑐ℎ = 2.26 kpc. The dependence on the stellar profile is weak in cuspy haloes and increases as the inner density slope approaches 𝛾 → 0. Middle panel… view at source ↗
Figure 2
Figure 2. Figure 2: Top panel: Time evolution of the bound subhalo mass fraction (black solid line) for a dSph with 𝑀ℎ = 109 𝑀⊙ and 𝑐ℎ = 2.26 kpc orbiting a MW￾like halo with 𝑀𝐺 = 1012 𝑀⊙ and 𝑐𝐺 = 21.5 kpc. The orbit has apocentre 𝑟apo = 120 kpc and pericentre 𝑟peri = 25 kpc (shown with a black dotted line in arbitrary units). The bound stellar mass fraction of smooth and clumpy halo models are shown in orange and blue, respe… view at source ↗
Figure 3
Figure 3. Figure 3: Left panel: Dwarf galaxy masses inferred from Equation (3), with (𝜂𝑟 , 𝜂𝜎 ) = (1, 0.5), as a function of half-light radius. Here, we adopt 𝛾 = 1, although 𝑀HA values are largely insensitive to this choice. Orange circles indicate galaxies with associated stellar streams, while green circles mark systems exhibiting signatures of tidal stripping (Pace et al. 2022). Magenta squares correspond to galaxies that… view at source ↗
read the original abstract

We use controlled $N$-body experiments to study the dynamical evolution of dwarf spheroidal galaxies (dSphs) embedded in dark-matter (DM) haloes containing a large population of dark subhaloes. We show that stellar orbits subject to stochastic force fluctuations irreversibly gain energy and expand toward a dynamical attractor characterized by a stellar half-light radius $r_{\rm half} \approx r_{\rm max}$ and a velocity dispersion $\sigma \approx 0.5\,v_{\rm max}$, where $v_{\rm max}$ is the peak circular velocity of the host halo at radius $r_{\rm max}$. This state is reached both in isolation and under tidal stripping, although tidal mass loss significantly accelerates the evolution. Assuming that the Milky Way (MW) dSphs have reached this state, we find that the inferred halo masses collapse onto narrow sequences as a function of $r_{\rm half}$. Under this assumption, MW satellites with $r_{\rm half} \lesssim 1\,\mathrm{kpc}$ follow the tidal tracks of cuspy haloes, while larger systems deviate in a manner consistent with cored DM profiles. Moreover, the mass--luminosity relation follows the slope expected from abundance matching, but with halo masses systematically lowered from their peak values at fixed luminosity. These results suggest that the structural diversity of dSphs is largely an evolutionary outcome driven by internal heating and tides, rather than by the conditions of star formation. This framework predicts that isolated, early-quenched dSphs should have systematically larger sizes than satellites, a prediction testable with upcoming surveys.

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

2 major / 2 minor

Summary. The paper uses controlled N-body experiments to demonstrate that stellar orbits in dwarf spheroidal galaxies (dSphs) embedded in dark-matter halos with abundant subhalos undergo irreversible energy gain from stochastic force fluctuations, driving expansion toward a dynamical attractor state with r_half ≈ r_max and σ ≈ 0.5 v_max. This attractor is reached both in isolation and under tidal stripping (with tides accelerating the process). Assuming all Milky Way dSphs have reached this state, the inferred halo masses collapse onto narrow sequences versus r_half: systems with r_half ≲ 1 kpc follow tidal tracks of cuspy halos, while larger ones deviate consistently with cored profiles. The mass-luminosity relation matches the slope from abundance matching but with systematically lower halo masses at fixed luminosity. The work concludes that dSph structural diversity is primarily an evolutionary outcome of internal heating and tides rather than star-formation conditions, and predicts that isolated early-quenched dSphs should be systematically larger than satellites.

Significance. If the central results hold, the paper provides a dynamical mechanism that could unify the observed range of dSph sizes and velocity dispersions, reduce tensions between cuspy and cored halo inferences, and shift emphasis from formation physics to post-infall evolution. The framework yields a clear, observationally testable prediction for field versus satellite dSph sizes that upcoming surveys can address. The N-body demonstration of an attractor reached under both isolated and tidal conditions is a concrete strength.

major comments (2)
  1. [Abstract and §3] Abstract and §3 (results on MW dSphs): The claim that inferred halo masses collapse onto narrow sequences as a function of r_half rests entirely on the assumption that every MW dSph has already reached the attractor; no independent verification is supplied, nor is any comparison given between the simulated heating timescale (set by subhalo density and mass function) and the known star-formation quenching times or orbital infall times of the MW satellites. If a subset of dSphs remain on the pre-attractor branch, both the reported mass sequences and the cuspy-versus-cored distinction would not hold.
  2. [§2] §2 (N-body setup): The manuscript provides no tabulated values for subhalo number density, mass function, particle number, softening lengths, or convergence tests that establish the robustness of the attractor location (r_half ≈ r_max, σ ≈ 0.5 v_max). Without these, it is impossible to judge whether the factor 0.5 is a robust outcome or sensitive to the specific subhalo population chosen.
minor comments (2)
  1. [Abstract] The abstract introduces r_max and v_max without a brief definition; a parenthetical reminder that these are the peak circular-velocity radius and value of the host halo would aid readability.
  2. [Figures] Figure captions should explicitly state the number of realizations and the range of subhalo parameters explored so that the attractor’s stability can be assessed at a glance.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment of the paper's significance and for the constructive major comments. We respond to each point below and indicate the revisions we will make.

read point-by-point responses

  1. Referee: [Abstract and §3] Abstract and §3 (results on MW dSphs): The claim that inferred halo masses collapse onto narrow sequences as a function of r_half rests entirely on the assumption that every MW dSph has already reached the attractor; no independent verification is supplied, nor is any comparison given between the simulated heating timescale (set by subhalo density and mass function) and the known star-formation quenching times or orbital infall times of the MW satellites. If a subset of dSphs remain on the pre-attractor branch, both the reported mass sequences and the cuspy-versus-cored distinction would not hold.

    Authors: We agree that the assumption that all observed MW dSphs have reached the attractor is central to the interpretation in §3 and that a direct comparison of timescales would strengthen the argument. The simulations demonstrate that the attractor is reached in both isolated and tidally stripped cases, with tides accelerating the process. In the revised manuscript we will add a dedicated paragraph in §3 that compares the simulated heating timescales (determined by the adopted subhalo population) to literature values for MW satellite infall times and quenching epochs. This addition will clarify the conditions under which the reported mass sequences and the cuspy/cored distinction are expected to hold. revision: yes

  2. Referee: [§2] §2 (N-body setup): The manuscript provides no tabulated values for subhalo number density, mass function, particle number, softening lengths, or convergence tests that establish the robustness of the attractor location (r_half ≈ r_max, σ ≈ 0.5 v_max). Without these, it is impossible to judge whether the factor 0.5 is a robust outcome or sensitive to the specific subhalo population chosen.

    Authors: We thank the referee for highlighting this presentational gap. While the relevant parameters are described in the text of §2, we will add a new table in the revised manuscript that explicitly lists the subhalo number density, mass function, particle numbers, softening lengths, and the results of the convergence tests performed. This will allow readers to assess directly the robustness of the attractor location, including the specific ratio σ ≈ 0.5 v_max. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper derives the dynamical attractor (r_half ≈ r_max and σ ≈ 0.5 v_max) from independent controlled N-body experiments, both in isolation and under tides. This result is not fitted to or defined by the MW dSph observations. The subsequent step of inferring halo masses and their collapse onto sequences as a function of r_half is explicitly conditional on the stated assumption that observed dSphs have reached the attractor state. No equations reduce by construction to prior inputs, no parameters are renamed as predictions, and no self-citation chains or uniqueness theorems are invoked to force the result. The mass sequences are an interpretive application of the simulation-derived relations to data, which is a standard non-circular step. The paper remains self-contained against external benchmarks such as the N-body runs.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The results rely on standard N-body simulation techniques without additional free parameters or invented entities beyond the approximate relations derived from the runs.

free parameters (1)
  • Approximate factor 0.5 in velocity dispersion relation
    The σ ≈ 0.5 v_max is presented as an approximate outcome of the simulations.
axioms (1)
  • standard math Collisionless dynamics in Newtonian gravity for stellar and dark matter particles.
    Core assumption of N-body simulations.

pith-pipeline@v0.9.0 · 5593 in / 1298 out tokens · 98188 ms · 2026-05-15T18:07:12.136850+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

23 extracted references · 23 canonical work pages

  1. [1]

    F., Fritz T

    Battaglia G., Taibi S., Thomas G. F., Fritz T. K., 2022, A&A, 657, A54. doi:10.1051/0004-6361/202141528 Boldrini P., 2021, Galax, 10,

    work page doi:10.1051/0004-6361/202141528 2022

  2. [2]

    doi:10.3390/galaxies10010005 Dehnen W., 1993, MNRAS, 265,

    work page doi:10.3390/galaxies10010005 1993

  3. [3]

    , year = 1993, month = nov, volume =

    doi:10.1093/mnras/265.1.250 Dutta Chowdhury D., van den Bosch F. C., van Dokkum P., Robles V. H., Schive H.-Y., Chiueh T., 2023, ApJ, 949,

    work page doi:10.1093/mnras/265.1.250 2023

  4. [5]

    G., 2018, MNRAS, 481,

    doi:10.1093/mnras/stad581 Errani R., Pe ˜narrubia J., Walker M. G., 2018, MNRAS, 481,

    work page doi:10.1093/mnras/stad581 2018

  5. [6]

    doi:10.1093/mnras/sty2505 Errani R., Pe ˜narrubia J., 2020, MNRAS, 491,

    work page doi:10.1093/mnras/sty2505 2020

  6. [7]

    F., 2021, MNRAS, 505,

    doi:10.1093/mnras/stz3349 Errani R., Navarro J. F., 2021, MNRAS, 505,

    work page doi:10.1093/mnras/stz3349 2021

  7. [8]

    F., Ibata R., Pe ˜narrubia J., 2022, MNRAS, 511,

    doi:10.1093/mnras/stab1215 Errani R., Navarro J. F., Ibata R., Pe ˜narrubia J., 2022, MNRAS, 511,

    work page doi:10.1093/mnras/stab1215 2022

  8. [9]

    G., Rozier S., Pe ˜narrubia J., Navarro J

    doi:10.1093/mnras/stac476 Errani R., Walker M. G., Rozier S., Pe ˜narrubia J., Navarro J. F., 2025, ApJ, 992,

    work page doi:10.1093/mnras/stac476 2025

  9. [10]

    I., Wyse R

    doi:10.3847/1538-4357/adfa27 Gilmore G., Wilkinson M. I., Wyse R. F. G., Kleyna J. T., Koch A., Evans N. W., Grebel E. K., 2007, ApJ, 663,

    work page doi:10.3847/1538-4357/adfa27 2007

  10. [11]

    doi:10.1086/518025 Hu C.-Y., Naab T., Glover S. C. O., Walch S., Clark P. C., 2017, MNRAS, 471,

    work page doi:10.1086/518025 2017

  11. [12]

    doi:10.1093/mnras/stx1773 Jethwa P., Erkal D., Belokurov V., 2018, MNRAS, 473,

    work page doi:10.1093/mnras/stx1773 2018

  12. [13]

    L., 1998, ARA&A, 36,

    doi:10.1093/mnras/stx2330 Mateo M. L., 1998, ARA&A, 36,

    work page doi:10.1093/mnras/stx2330 1998

  13. [14]

    J., Crnojevi ´c D., Drlica-Wagner A., Caldwell N., Guhathakurta P., Seth A

    doi:10.1146/annurev.astro.36.1.435 Mutlu-Pakdil B., Sand D. J., Crnojevi ´c D., Drlica-Wagner A., Caldwell N., Guhathakurta P., Seth A. C., et al., 2021, ApJ, 918,

    work page doi:10.1146/annurev.astro.36.1.435 2021

  14. [16]

    O., Drlica-Wagner A., Bechtol K., Mau S., Wechsler R

    doi:10.3847/1538-4357/ab94b0 Nadler E. O., Drlica-Wagner A., Bechtol K., Mau S., Wechsler R. H., Gluscevic V., Boddy K., et al., 2021, PhRvL, 126, 091101. doi:10.1103/PhysRevLett.126.091101 Pace A. B., Erkal D., Li T. S., 2022, ApJ, 940,

    work page doi:10.3847/1538-4357/ab94b0 2021

  15. [17]

    S., Bindi, V., Zhang, M., Feng, X., 2019, A Numerical Study of Cosmic Proton Modulation Using AMS-02 Observations , Astrophys

    doi:10.3847/1538- 4357/ac997b Pe˜narrubia J., McConnachie A. W., Navarro J. F., 2008a, ApJ, 672,

    work page doi:10.3847/1538-

  16. [18]

    F., McConnachie A

    doi:10.1086/521543 Pe˜narrubia J., Navarro J. F., McConnachie A. W., 2008b, ApJ, 673,

    work page doi:10.1086/521543

  17. [19]

    J., Walker M

    doi:10.1086/523686 Pe˜narrubia J., Benson A. J., Walker M. G., Gilmore G., McConnachie A. W., Mayer L., 2010, MNRAS, 406,

    work page doi:10.1086/523686 2010

  18. [20]

    doi:10.1111/j.1365- 2966.2010.16762.x Pe˜narrubia J., 2023, MNRAS, 519,

    work page doi:10.1111/j.1365- 2010

  19. [21]

    G., 2025, MNRAS.tmp

    doi:10.1093/mnras/stac3642 Pe˜narrubia J., Errani R., Vitral E., Walker M. G., 2025, MNRAS.tmp. doi:10.1093/mnras/staf1800 (P25) Read J. I., Mamon G. A., Vasiliev E., Watkins L. L., Walker M. G., Pe˜narrubia J., Wilkinson M., et al., 2021, MNRAS, 501,

    work page doi:10.1093/mnras/stac3642 2025

  20. [22]

    D., 2019, ARA&A, 57,

    doi:10.1093/mnras/staa3663 Simon J. D., 2019, ARA&A, 57,

    work page doi:10.1093/mnras/staa3663 2019

  21. [23]

    L.; Errani, R.; Pe ˜narrubia, J.; Walker, M

    doi:10.1146/annurev-astro-091918- 104453 Splawska, S. L.; Errani, R.; Pe ˜narrubia, J.; Walker, M. G., 2026; arXiv:2602.11273 Tulin S., Yu H.-B., 2018, PhR, 730,

    work page doi:10.1146/annurev-astro-091918- 2026

  22. [24]

    2018, Physics Reports, 730, 1, 10.1016/j.physrep.2017.11.004

    doi:10.1016/j.physrep.2017.11.004 Vitral E., Pe ˜narrubia J., Walker M. G., 2025a, MNRAS.tmp. doi:10.1093/mnras/staf2013 Vitral E., van der Marel R. P., Sohn S. T., Pe ˜narrubia J., Patel E., Watkins L. L., Libralato M., et al., 2025b, arXiv, arXiv:2508.20711. doi:10.48550/arXiv.2508.20711 Walker M. G., Mateo M., Olszewski E. W., Pe ˜narrubia J., Evans N....

    work page doi:10.1016/j.physrep.2017.11.004 2017

  23. [25]

    doi:10.1088/0004-637X/704/2/1274 MNRAS000, 000–000 (0000)

    work page doi:10.1088/0004-637x/704/2/1274