pith. sign in

arxiv: 2604.06304 · v1 · submitted 2026-04-07 · 🌌 astro-ph.GA

Collisional Dynamics of Stars and Dark Matter in Ultra-Faint Galaxies

Rapha\"el Errani , Nicolas Esser , Jorge Pe\~narrubia , Matthew G. Walker This is my paper

Pith reviewed 2026-05-10 19:49 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords ultra-faint galaxiesdynamical frictiondark matter cuspsN-body simulationsstellar binariestidal strippinggalaxy evolutiondark matter cores
0
0 comments X

The pith

Dynamical friction between stars and subsolar-mass dark matter particles depletes central dark matter in ultra-faint galaxies, turning cusps into cores.

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

The paper employs controlled N-body simulations to examine the collisional exchange of energy between stars and dark matter in ultra-faint galaxies. Dynamical friction causes stars to transfer energy to lighter dark matter particles, driving dark matter outward from the center while stars sink inward. This process flattens dark matter density profiles from cusps to constant-density cores and proves most efficient in galaxies that are tidally limited and have low stellar velocity dispersion. The stellar component compacts into a dense central cluster that becomes baryon-dominated within its half-light radius while retaining a surrounding dark matter halo, and the same interactions promote the formation of stellar binaries. The contraction slows as central dark matter densities fall and binaries form, illustrating how the faintest galaxies can evolve toward structures that resemble globular clusters.

Core claim

Controlled N-body simulations reveal that dynamical friction between stars and subsolar-mass dark matter particles depletes dark matter from the centers of ultra-faint galaxies, transforming cusps into constant-density cores. This effect is strongest in tidally limited systems with low stellar velocity dispersion. High-mass stars sink to the center, reducing the dark matter to stellar mass ratio within the half-light radius and compacting the stars into a dense, baryon-dominated cluster surrounded by a dark matter halo. The collisional cooling also facilitates the formation of stellar binaries, with the contraction eventually slowing as central dark matter densities drop and binaries form.

What carries the argument

Dynamical friction between stars and subsolar-mass dark matter particles, which transfers orbital energy from stars to dark matter and thereby depletes dark matter from galactic centers while allowing stars to sink inward.

If this is right

  • The dynamical-to-stellar mass ratio within the stellar half-light radius decreases monotonically as high-mass stars sink toward the center.
  • The stellar population compacts into a dense, baryon-dominated cluster surrounded by a dark matter halo.
  • Such a cluster shares the chemical composition of an ultra-faint galaxy yet appears virtually dark matter-free within its half-light radius.
  • Collisional cooling with dark matter particles provides an efficient pathway for the formation of stellar binaries in the contracting cluster.
  • The contraction slows due to decreasing central dark matter densities and the formation of stellar binaries.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Similar collisional evolution could occur in other low-mass, dark-matter-dominated systems that experience tidal limitation.
  • Surveys searching for dense stellar clusters with ultra-faint-galaxy chemical signatures but little internal dark matter could identify objects shaped by this process.
  • Measurements of elevated binary fractions or flattened central dark matter profiles in ultra-faint galaxies would offer direct observational tests.
  • The dynamical similarity between the end states and globular clusters suggests that some objects currently classified as galaxies may be reclassified based on evolutionary history.

Load-bearing premise

Dark matter particles must have masses much smaller than the Sun for dynamical friction with stars to operate efficiently, and the galaxies must be tidally limited with low stellar velocity dispersion.

What would settle it

N-body simulations with dark matter particle masses at or above solar mass that show no cusp-to-core transformation in tidally limited ultra-faint galaxy models, or observations of ultra-faint galaxies that retain cuspy dark matter profiles despite low velocity dispersions and tidal limitation, would falsify the central claim.

Figures

Figures reproduced from arXiv: 2604.06304 by Jorge Pe\~narrubia, Matthew G. Walker, Nicolas Esser, Rapha\"el Errani.

Figure 1
Figure 1. Figure 1: Collisional energy exchange between stars and dark matter transforms a dark matter cusp into a constant-density core, as high-mass stars sink to the center of the ultra-faint galaxy. Top row: Simulation snapshots taken at 𝑡 = 0, 0.5, 1 and 2 Gyr with parameters as listed in [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: The initial dark matter cusp is transformed into a constant-density core. Left panel: Radial dark matter density profile 𝜌DM (𝑟 ) for the simulations with an initial dynamical-to-stellar mass ratio of Υdyn0 = 3 (orange), 10 (green) and 30 (purple). For comparison, the dark matter-only run is shown in black. The smaller Υdyn0 is, the larger the size of the constant-density core. To reduce noise, for each va… view at source ↗
Figure 4
Figure 4. Figure 4: Dynamical friction exerted by dark matter particles on stars drives the contraction of stellar half-light radii. Shown is the evolution of the half-light radii 𝑟h for the two stellar populations of low-mass and high-mass stars (yellow and blue lines, showing the median over all four random 𝑁-body realizations; shaded bands span the 16th–84th percentile of the underlying distribution) for the Υdyn0 = 10 mod… view at source ↗
Figure 6
Figure 6. Figure 6: As the population of high-mass stars contracts, the number of dynamically formed binaries 𝑁bin increases. Top panel: Shown is the binary fraction 𝑓bin ≡ 𝑁bin/𝑁★ for binaries consisting of two high-mass stars as a function of time. For the model with Υdyn0 = 10, the number of binaries with semi-major axes 𝑎bin ≤ 0.1 pc rises sharply around ∼2 Gyr. Shaded bands span the 16th–84th percentile of the underlying… view at source ↗
Figure 7
Figure 7. Figure 7: Evolution of the total energy 𝐸 (𝑡) − 𝐸 (0) and the kinetic energy 𝐾 (𝑡) − 𝐾 (0) for dark matter (gray) and high-mass stars (color coded by Υdyn). As the high-mass stars contract, heat is transferred to the dark matter, and both the stellar- and the dark matter systems have a positive heat capacity 𝜕𝐸/𝜕𝐾 > 0. Once the stars become self-gravitating (Υdyn ≈ 1), the stellar population heats up as it contracts… view at source ↗
read the original abstract

We use controlled N-body simulations to study the collisional exchange of energy between stars and dark matter in ultra-faint galaxies. We find that dynamical friction between stars and subsolar-mass dark matter particles results in the depletion of dark matter from the galaxies' centers, thereby transforming dark matter cusps into constant-density cores. The process is particularly effective in tidally limited galaxies with low stellar velocity dispersion. As high-mass stars sink toward the center of the dark matter halo, the dynamical-to-stellar mass ratio within the stellar half-light radius decreases monotonically. The stellar population of a dark matter-dominated galaxy is thereby compacted into a dense, baryon-dominated cluster, surrounded by a dark matter halo. Such a cluster would share the chemical composition of an ultra-faint galaxy, yet would be virtually dark matter-free within its half-light radius. We moreover find that the collisional cooling with dark matter particles provides an efficient pathway for the formation of stellar binaries in the contracting cluster. The contraction is eventually slowed down due to the decreasing central dark matter densities and the formation of stellar binaries. Our models highlight that the dynamical processes governing the faintest galaxies give rise to a rich phenomenology, blurring the line between the dynamics of globular clusters and galaxies.

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 / 0 minor

Summary. The manuscript presents controlled N-body simulations of collisional energy exchange between stars and dark matter in ultra-faint galaxies. It claims that dynamical friction between stars and subsolar-mass dark matter particles depletes dark matter from galaxy centers, converting cusps into constant-density cores. This effect is strongest in tidally limited systems with low stellar velocity dispersion, leading to monotonic decrease in the dynamical-to-stellar mass ratio within the half-light radius, compaction of stars into a dense baryon-dominated cluster surrounded by a dark matter halo, and efficient formation of stellar binaries via collisional cooling. The contraction slows due to declining central dark matter densities and binary formation, with the models highlighting a rich phenomenology that blurs globular clusters and galaxies.

Significance. If the results hold under the stated assumptions, the work provides a novel collisional pathway for cusp-to-core transformation in ultra-faint galaxies without relying on baryonic feedback. It demonstrates how such dynamics can produce virtually dark matter-free stellar clusters sharing ultra-faint galaxy chemical compositions, offering a mechanism to explain observed structures and the dynamical overlap between the faintest galaxies and globular clusters. The use of controlled N-body simulations enables direct, parameter-controlled exploration of these effects.

major comments (2)
  1. Abstract and simulation description: the central claim of cusp-to-core transformation via dynamical friction rests entirely on the assumption of subsolar-mass dark matter particles, yet no sensitivity tests to this mass choice, comparisons to standard CDM particle masses, or independent justification against ultra-faint galaxy constraints are provided; this assumption is load-bearing and not independently validated.
  2. Abstract: no details are given on numerical resolution (particle number, softening lengths, time-stepping, or convergence tests), nor is there quantitative validation of the reported energy exchange or friction timescales against analytic expectations for dynamical friction in the low-dispersion regime; this undermines assessment of whether the depletion and core formation are numerically robust.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed review of our manuscript. We address each major comment below and will revise the paper accordingly to strengthen the presentation of our results.

read point-by-point responses

  1. Referee: Abstract and simulation description: the central claim of cusp-to-core transformation via dynamical friction rests entirely on the assumption of subsolar-mass dark matter particles, yet no sensitivity tests to this mass choice, comparisons to standard CDM particle masses, or independent justification against ultra-faint galaxy constraints are provided; this assumption is load-bearing and not independently validated.

    Authors: We acknowledge that the subsolar-mass DM particle assumption is central to realizing the collisional regime in our controlled N-body experiments, where individual star-DM encounters enable efficient dynamical friction and the reported cusp-to-core evolution. This regime differs from the collisionless limit of standard CDM, in which simulation particles are typically far more massive. We will add a dedicated subsection (and associated figures) presenting sensitivity tests across a range of DM particle masses, including higher-mass cases that approach standard CDM resolutions, to quantify the mass threshold below which the core-formation effect operates. We will also expand the introduction and discussion to provide context for exploring this mass scale, including references to literature on possible DM properties or resolution considerations in ultra-faint systems, while noting that direct observational constraints on DM particle masses remain limited. revision: yes

  2. Referee: Abstract: no details are given on numerical resolution (particle number, softening lengths, time-stepping, or convergence tests), nor is there quantitative validation of the reported energy exchange or friction timescales against analytic expectations for dynamical friction in the low-dispersion regime; this undermines assessment of whether the depletion and core formation are numerically robust.

    Authors: We agree that the abstract should summarize the numerical setup and that direct comparison to analytic expectations is needed to establish robustness. Although the methods section describes the N-body integrator, particle numbers, softening, and time-stepping criteria, we will revise the abstract to include these key parameters. We will also add a new subsection that quantitatively validates the simulated energy-exchange rates and dynamical-friction timescales against the Chandrasekhar formula adapted to the low-dispersion regime, together with convergence tests that vary particle number and softening length. These additions will confirm that the central DM depletion and stellar compaction are not sensitive to numerical choices. revision: yes

Circularity Check

0 steps flagged

No significant circularity in simulation-based results

full rationale

The paper reports outcomes from controlled N-body simulations that integrate Newtonian equations of motion forward in time for stars and dark matter particles. No analytical derivation chain, fitted parameters, or self-referential definitions are invoked to produce the reported cusp-to-core transformation or cluster compaction; these emerge directly as numerical results under the chosen initial conditions and particle masses. The abstract and described setup contain no load-bearing self-citations, ansatzes smuggled via prior work, or renaming of known results as new predictions. The central claim is therefore self-contained as a direct consequence of the simulation methodology rather than equivalent to its inputs by construction.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim depends on the assumption of subsolar-mass dark matter particles and on the galaxies being tidally limited with low stellar velocity dispersion; no new entities are introduced.

free parameters (1)
  • dark matter particle mass = subsolar
    Set to subsolar values to enable efficient dynamical friction with stars.
axioms (2)
  • standard math Newtonian gravity and collisionless N-body dynamics govern star-dark matter interactions
    Basis for all reported simulation outcomes.
  • domain assumption The modeled galaxies are tidally limited and possess low stellar velocity dispersion
    Explicitly stated as the regime where the process is particularly effective.

pith-pipeline@v0.9.0 · 5532 in / 1350 out tokens · 42086 ms · 2026-05-10T19:49:55.371826+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

58 extracted references · 58 canonical work pages

  1. [1]

    C., & Evans, N

    Amorisco, N. C., & Evans, N. W. 2012, MNRAS, 419, 184

    work page 2012

  2. [2]

    W., & Diaferio, A

    Angus, G. W., & Diaferio, A. 2009, MNRAS, 396, 887

    work page 2009

  3. [3]

    Banik, U., & van den Bosch, F. C. 2022, ApJ, 926, 215

    work page 2022

  4. [4]

    1986, Nature, 324, 446

    Barnes, J., & Hut, P. 1986, Nature, 324, 446

    work page 1986

  5. [5]

    2022, MNRAS, 510, 3531

    Hilker, M. 2022, MNRAS, 510, 3531

    work page 2022

  6. [6]

    1987, Galactic dynamics (Princeton University Press)

    Binney, J., & Tremaine, S. 1987, Galactic dynamics (Princeton University Press)

    work page 1987

  7. [7]

    2022, MNRAS, 509, 5330

    Santos-Santos, I. 2022, MNRAS, 509, 5330

    work page 2022

  8. [8]

    S., Pace, A

    Cerny, W., Li, T. S., Pace, A. B., et al. 2026a, arXiv:2602.17652

    work page arXiv

  9. [9]

    2003, PASP, 115, 763

    Chabrier, G. 2003, PASP, 115, 763

    work page 2003

  10. [10]

    1941, ApJ, 94, 511 —

    Chandrasekhar, S. 1941, ApJ, 94, 511 —. 1943, ApJ, 97, 255

    work page 1941

  11. [11]

    R., Dehnen, W., Read, J

    Cole, D. R., Dehnen, W., Read, J. I., & Wilkinson, M. I. 2012, MNRAS, 426, 601

    work page 2012

  12. [12]

    C., Jerjen, H., Kim, D., & Schirmer, M

    Conn, B. C., Jerjen, H., Kim, D., & Schirmer, M. 2018, ApJ, 852, 68 Crnogorčević, M., & Linden, T. 2024, PhRvD, 109, 083018

    work page 2018

  13. [13]

    Delos, M. S. 2019, PhRvD, 100, 083529 Devlin,S.,Baumgardt,H.,&Sweet,S.M.2025,MNRAS,539,2485 Di Cintio, A., Brook, C. B., Dutton, A. A., et al. 2014, MNRAS, 441, 2986

    work page 2019

  14. [14]

    2007, ApJ, 657, 262

    Diemand, J., Kuhlen, M., & Madau, P. 2007, ApJ, 657, 262

    work page 2007

  15. [15]

    2001, ApJ, 560, 636

    El-Zant, A., Shlosman, I., & Hoffman, Y. 2001, ApJ, 560, 636

    work page 2001

  16. [16]

    Errani, R., & Navarro, J. F. 2021, MNRAS, 505, 18

    work page 2021

  17. [17]

    F., Ibata, R., & Peñarrubia, J

    Errani, R., Navarro, J. F., Ibata, R., & Peñarrubia, J. 2022, MNRAS, 511, 6001

    work page 2022

  18. [18]

    2020, MNRAS, 491, 4591

    Errani, R., & Peñarrubia, J. 2020, MNRAS, 491, 4591

    work page 2020

  19. [19]

    Errani, R., Peñarrubia, J., Laporte, C. F. P., & Gómez, F. A. 2017, MNRAS, 465, L59 Collisional Dynamics of Stars and Dark Matter 11

    work page 2017

  20. [20]

    Errani, R., Peñarrubia, J., & Walker, M. G. 2018, MNRAS, 481, 5073 —. 2025, ApJ, 993, 160

    work page 2018

  21. [21]

    2026, arXiv:2602.20245

    Esser, N. 2026, arXiv:2602.20245

    work page arXiv 2026

  22. [22]

    2022, arXiv:2201.09788

    Facchinetti, G., Stref, M., & Lavalle, J. 2022, arXiv:2201.09788

    work page arXiv 2022

  23. [23]

    Fellhauer, M., & Lin, D. N. C. 2007, MNRAS, 375, 604

    work page 2007

  24. [24]

    X., et al

    Geha, M., Pelliccia, D., Prochaska, J. X., et al. 2026, ApJ, 999, 140

    work page 2026

  25. [25]

    Giersz, M., & Heggie, D. C. 1994, MNRAS, 270, 298

    work page 1994

  26. [26]

    I., Stadel, J., & Zemp, M

    Goerdt, T., Moore, B., Read, J. I., Stadel, J., & Zemp, M. 2006, MNRAS, 368, 1073

    work page 2006

  27. [27]

    1998, MNRAS, 297, 517

    Hernandez, X., & Gilmore, G. 1998, MNRAS, 297, 517

    work page 1998

  28. [28]

    A., Lewis, G

    Ibata, R. A., Lewis, G. F., Irwin, M. J., & Quinn, T. 2002, MNRAS, 332, 915

    work page 2002

  29. [29]

    P., Simon, J

    Ji, A. P., Simon, J. D., Frebel, A., Venn, K. A., & Hansen, T. T. 2019, ApJ, 870, 83

    work page 2019

  30. [30]

    V., Spergel, D

    Johnston, K. V., Spergel, D. N., & Haydn, C. 2002, ApJ, 570, 656

    work page 2002

  31. [31]

    G., Peñarrubia, J., & Koposov, S

    Kervick, C., Walker, M. G., Peñarrubia, J., & Koposov, S. E. 2022, ApJ, 929, 77

    work page 2022

  32. [32]

    P., Mackey, D., & Da Costa, G

    Kim, D., Jerjen, H., Milone, A. P., Mackey, D., & Da Costa, G. S. 2015, ApJ, 803, 63

    work page 2015

  33. [33]

    Laporte, C. F. P., & Peñarrubia, J. 2015, MNRAS, 449, L90

    work page 2015

  34. [34]

    2010, A&A, 509, A46

    Longhitano, M., & Binggeli, B. 2010, A&A, 509, A46

    work page 2010

  35. [35]

    B., et al

    Mau, S., Cerny, W., Pace, A. B., et al. 2020, ApJ, 890, 136 Moliné, Á., Sánchez-Conde, M. A., Palomares-Ruiz, S., & Prada, F. 2017, MNRAS, 466, 4974

    work page 2020

  36. [36]

    F., Frenk, C

    Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 563

    work page 1996

  37. [37]

    F., Frenk, C

    Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 Oñorbe, J., Boylan-Kolchin, M., Bullock, J. S., et al. 2015, MNRAS, 454, 2092 Peñarrubia, J., Benson, A. J., Walker, M. G., et al. 2010, MNRAS, 406, 1290 Peñarrubia, J., Errani, R., Vitral, E., & Walker, M. G. 2025, MNRAS, 544, 2557 Peñarrubia, J., Ludlow, A. D., Chanamé, J., & Walker, M....

    work page 1997

  38. [38]

    I., Goerdt, T., Moore, B., et al

    Read, J. I., Goerdt, T., Moore, B., et al. 2006, MNRAS, 373, 1451

    work page 2006

  39. [39]

    D., & Loeb, A

    Safarzadeh, M., Simon, J. D., & Loeb, A. 2022, ApJ, 930, 54

    work page 2022

  40. [40]

    2025, PASP, 137, 104103

    Shariat, C., El-Badry, K., Gennaro, M., et al. 2025, PASP, 137, 104103

    work page 2025

  41. [41]

    D., Li, T

    Simon, J. D., Li, T. S., Drlica-Wagner, A., et al. 2017, ApJ, 838, 11

    work page 2017

  42. [42]

    D., Li, T

    Simon, J. D., Li, T. S., Ji, A. P., et al. 2024, ApJ, 976, 256

    work page 2024

  43. [43]

    Smith, S. E. T., Cerny, W., Hayes, C. R., et al. 2024, ApJ, 961, 92

    work page 2024

  44. [44]

    L., Errani, R., Peñarrubia, J., & Walker, M

    Splawska, S. L., Errani, R., Peñarrubia, J., & Walker, M. G. 2026, arXiv:2602.11273

    work page arXiv 2026

  45. [45]

    2005, MNRAS, 364, 1105

    Springel, V. 2005, MNRAS, 364, 1105

    work page 2005

  46. [46]

    2021, MNRAS, 506, 2871

    Springel, V., Pakmor, R., Zier, O., & Reinecke, M. 2021, MNRAS, 506, 2871

    work page 2021

  47. [47]

    Springel, V., Yoshida, N., & White, S. D. M. 2001, NewA, 6, 79

    work page 2001

  48. [48]

    2008, MNRAS, 391, 1685 Stücker, J., Ogiya, G., Angulo, R

    Springel, V., Wang, J., Vogelsberger, M., et al. 2008, MNRAS, 391, 1685 Stücker, J., Ogiya, G., Angulo, R. E., Aguirre-Santaella, A., & Sánchez-Conde, M. A. 2023, MNRAS, 521, 4432

    work page 2008

  49. [49]

    Y., Cerny, W., Pace, A

    Tan, C. Y., Cerny, W., Pace, A. B., et al. 2026, ApJ, 1000, 46

    work page 2026

  50. [50]

    R., & White, S

    Tormen, G., Bouchet, F. R., & White, S. D. M. 1997, MNRAS, 286, 865

    work page 1997

  51. [51]

    2025, arXiv:2509.15345 van den Bosch, F

    Untzaga, J., Mezcua, M., Bonoli, S., et al. 2025, arXiv:2509.15345 van den Bosch, F. C., Ogiya, G., Hahn, O., & Burkert, A. 2018, MNRAS, 474, 3043

    work page arXiv 2025

  52. [52]

    A., Irwin, M., Shetrone, M

    Venn, K. A., Irwin, M., Shetrone, M. D., et al. 2004, AJ, 128, 1177

    work page 2004

  53. [53]

    G., Combet, C., Hinton, J

    Walker, M. G., Combet, C., Hinton, J. A., Maurin, D., & Wilkinson, M. I. 2011, ApJL, 733, L46

    work page 2011

  54. [54]

    S., et al

    Wang, J., Bose, S., Frenk, C. S., et al. 2020, Nature, 585, 39

    work page 2020

  55. [55]

    White, S. D. M., & Rees, M. J. 1978, MNRAS, 183, 341

    work page 1978

  56. [56]

    D., Bullock, J

    Wolf, J., Martinez, G. D., Bullock, J. S., et al. 2010, MNRAS, 406, 1220

    work page 2010

  57. [57]

    R., et al

    Zaremba, D., Venn, K., Hayes, C. R., et al. 2025, ApJ, 987, 217

    work page 2025

  58. [58]

    S., et al

    Zheng, H., Bose, S., Frenk, C. S., et al. 2024, MNRAS, 528, 7300

    work page 2024