pith. sign in

arxiv: 2606.30078 · v1 · pith:AL64G5W6new · submitted 2026-06-29 · 🌌 astro-ph.GA · astro-ph.CO

Strong Stellar Diffusion from Wave DM Cosmological Simulation and Potential Unified Origin for dSphs, UFGs, and UDGs

A. Pozo , T. Broadhurst , J. Zhang , M. Oguri , K. Umetsu , R. Emami , L. Hernquist , P. Mocz
show 1 more author
M. Vogelsberger
This is my paper

Pith reviewed 2026-06-30 05:11 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.CO
keywords wave dark matterψDMstellar diffusiondwarf spheroidal galaxiesultra-diffuse galaxiessolitonSersic indexcosmological simulations
0
0 comments X

The pith

Wave dark matter simulations show stars diffuse outward in halos via soliton random walk, producing Gaussian profiles that expand with the square root of time.

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

The paper uses cosmological simulations of wave dark matter to demonstrate that stars undergo random-walk diffusion driven by intrinsic wave perturbations, spreading from dense cores into extended halos over cosmic time. This produces locally Gaussian stellar distributions with Sersic index 0.5 whose half-light radii grow proportionally to the square root of time scaled by the inverse square root of the dark matter particle mass. Diffusion strength increases with halo mass, so more massive systems become progressively more diffuse. The result interprets the observed sequence from compact dwarf spheroidals to ultra-diffuse galaxies as a simple age progression rather than requiring separate formation channels.

Core claim

Our ψDM simulations predict that stars diffuse throughout dark matter halos over the Hubble time through a random walk driven by the wave perturbations intrinsic to ψDM. The resulting stellar distribution locally follows a Gaussian profile (Sersic index n=0.5), expanding as R_{1/2}(t)≃(ℏ/m_ψ)^{0.5}√t, in good agreement with the core--halo profiles of typical ψDM dwarf spheroidal galaxies. The strength of this diffusion depends on halo mass and the corresponding soliton, naturally producing progressively more diffuse stellar systems in more massive halos. The observed continuity from faint dwarfs and compact dwarf spheroidals to ultra-diffuse galaxies can therefore be interpreted as an age se

What carries the argument

The random walk of the central soliton induced by wave perturbations intrinsic to ψDM, which scatters stars outward from the dense core into the halo.

If this is right

  • Stellar half-light radii follow R_{1/2}(t) ≃ (ℏ/m_ψ)^{0.5} √t and match core-halo structures in dwarf spheroidals.
  • Later-forming systems remain compact and bright while earlier ones become more diffuse, unifying dSphs, UFGs, and UDGs without distinct channels.
  • Stellar scattering creates extended envelopes around Local Group dwarfs and places globular clusters at larger radii.
  • Diffuse stellar halos observed by Euclid are direct signatures of wave dark matter granular dynamics across a wide mass range.

Where Pith is reading between the lines

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

  • Stellar population ages in ultra-diffuse versus compact dwarfs should correlate with their current sizes if the mechanism is at work.
  • The same diffusion process could be checked in higher-mass halos to test whether wave dark matter continues to shape galaxy structure beyond dwarfs.
  • Varying the particle mass m_ψ in simulations would predict measurable changes in the present-day size distribution of galaxies.
  • Globular cluster radial distributions in dwarfs could provide an independent tracer of the cumulative diffusion time.

Load-bearing premise

Stellar diffusion is produced solely by the soliton’s random walk from wave perturbations, with its strength set by halo mass and soliton properties so that observed galaxy sizes form a continuous age sequence.

What would settle it

A ψDM simulation or set of observed dwarf galaxies whose stellar half-light radii show no sqrt(t) growth, lack Gaussian profiles, or fail to increase in diffuseness with halo mass.

Figures

Figures reproduced from arXiv: 2606.30078 by A. Pozo, J. Zhang, K. Umetsu, L. Hernquist, M. Oguri, M. Vogelsberger, P. Mocz, R. Emami, T. Broadhurst.

Figure 1
Figure 1. Figure 1 [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Same as [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Radial migration properties of stars in the most massive simulated galaxy at z = 5.56 for ψDM (top row), “WDM” (middle row), and CDM (bottom row), shown in comoving units. The panels show: present-day radius versus birth radius, the distribution of radial changes ∆r = rnow −rbirth, the birth-radius distributions of inward and outward migrating stars, their cumulative distribution functions, and the statist… view at source ↗
Figure 4
Figure 4. Figure 4: Evolution of the oldest stellar population in the ψDM galaxy, considering only stars already formed at z ≥ 10.86, shown in comoving units. The top-left panel shows the evolution of the stellar density contours from z = 10.86 to z = 5.56, illustrating the progressive expansion and diffusion of the initial compact stellar core. The top-middle panel presents the evolution of the mean radial displacement, whil… view at source ↗
Figure 5
Figure 5. Figure 5: Evolution of the stellar population born within the central core of the ψDM galaxy at all redshifts between z = 10.86 and z = 5.56, shown in comoving units. The top-left panel shows the evolution of the stellar density contours, illustrating the gradual expansion of the stellar distribution over time despite the continuous replenishment of the central core by newly formed stars. The top-middle and top-righ… view at source ↗
Figure 6
Figure 6. Figure 6: Properties of the stellar population constituting the central core of the ψDM galaxy at z = 5.56, shown in comoving units. The top panels display the present-day stellar core in three different projections, illustrating the compact central structure surrounded by a diffuse stellar component. The bottom panels characterize the migration properties of these stars through the comparison of present-day and bir… view at source ↗
Figure 7
Figure 7. Figure 7: Orbital evolution of stars born in the central core of the ψDM galaxy at z ≥ 10.86, shown in comoving units. Each panel displays the stellar distribution at different redshifts, illustrating the progressive outward diffusion of stellar orbits over time. The black points mark the initial compact stellar core, while the coloured stellar distributions show how stars become increasingly extended and diffuse du… view at source ↗
Figure 8
Figure 8. Figure 8: Summary of the soliton-driven stellar diffusion scenario in ψDM halos. The upper panel illustrates how the stochastic random walk of the central soliton scatters stars from the inner core toward the outer halo, naturally producing increasingly diffuse stellar systems from dwarf galaxies to UDGs. The lower panel summarizes the numerical evidence from the simulations, showing the outward migration of stars, … view at source ↗
Figure 9
Figure 9. Figure 9: Predicted stellar expansion produced by soliton-driven scattering as a function of halo mass in the ψDM framework. Left panel: halo masses are estimated directly from the core radii reported in Pozo et al. (2024b), although these measurements become increasingly uncertain for the most compact systems. Middle panel: halo masses are inferred from the observed velocity dispersions using the boson masses prefe… view at source ↗
read the original abstract

Our $\psi$DM simulations predict that stars diffuse throughout dark matter halos over the Hubble time through a random walk driven by the wave perturbations intrinsic to $\psi$DM. The resulting stellar distribution locally follows a Gaussian profile (Sersic index $n=0.5$), as expected from the central limit theorem, expanding as $R_{1/2}(t)\simeq(\hbar/m_\psi)^{0.5}\sqrt{t}$, in good agreement with the core--halo profiles of typical $\psi$DM dwarf spheroidal galaxies. The strength of this diffusion depends on halo mass and the corresponding soliton, naturally producing progressively more diffuse stellar systems in more massive halos. The observed continuity from faint dwarfs and compact dwarf spheroidals to ultra-diffuse galaxies can therefore be interpreted as an age sequence, with later-forming dwarfs experiencing less diffusion and thus remaining smaller and brighter. Stellar scattering arises from the random walk of the soliton, gradually transporting stars from the dense central core into the outer halo, creating the extended stellar envelopes observed around Local Group dwarfs. Rather than being unique to Ultra-Diffuse Galaxies (UDGs), this wave-driven stellar diffusion may provide a unified mechanism explaining galaxy structure across a vast mass range, from ultra-faint dwarfs to the most massive UDGs, without requiring distinct formation channels or "failed galaxy" scenarios. The diffuse stellar halos and globular cluster distributions recently revealed by Euclid may therefore represent direct observational signatures of the granular dynamics of wave dark matter.

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 paper reports results from ψDM cosmological simulations indicating that stars diffuse throughout dark matter halos over a Hubble time via random walks driven by wave perturbations intrinsic to ψDM. This produces locally Gaussian stellar distributions (Sersic index n=0.5) expanding as R_{1/2}(t) ≃ (ℏ/m_ψ)^{0.5} √t, with diffusion strength depending on halo mass and soliton properties. The work interprets the observed continuity from dSphs and UFGs to UDGs as an age sequence arising naturally without distinct formation channels, and suggests this as a unified mechanism with potential signatures in Euclid data.

Significance. If the simulation results hold, the work provides a mechanism tying wave dark matter granular dynamics directly to stellar structure across dwarf to UDG scales, offering a parameter-light explanation for extended envelopes and a falsifiable age-sequence interpretation. The invocation of the central limit theorem for the n=0.5 profile is a clear strength, as is the absence of invented entities beyond the standard ψDM framework.

major comments (2)
  1. [Abstract] Abstract: the manuscript asserts that 'our ψDM simulations predict' the diffusion, scaling, and mass-dependent continuity, but supplies no simulation details, resolution checks, error estimates, or quantitative fit statistics, so the data-to-claim link cannot be verified.
  2. [Abstract] Abstract: the R_{1/2}(t) ≃ (ℏ/m_ψ)^{0.5} √t scaling follows directly from the random-walk assumption and central limit theorem as stated; m_ψ remains a free parameter whose value can be adjusted to match observed profiles, so the result is not an independent prediction of the simulations beyond the modeling assumption.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thoughtful review and constructive comments on our manuscript. We address each major comment below in detail. Revisions will be made to improve clarity and strengthen the link between our simulation results and the claims in the abstract.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the manuscript asserts that 'our ψDM simulations predict' the diffusion, scaling, and mass-dependent continuity, but supplies no simulation details, resolution checks, error estimates, or quantitative fit statistics, so the data-to-claim link cannot be verified.

    Authors: The full manuscript contains the requested details on simulation setup, resolution convergence tests, error estimates, and quantitative profile fits in the Methods and Results sections. We agree, however, that the abstract is too concise to convey this information adequately. We will revise the abstract to include a brief reference to the simulation methodology and key validation statistics, thereby strengthening the data-to-claim connection without altering the length substantially. revision: yes

  2. Referee: [Abstract] Abstract: the R_{1/2}(t) ≃ (ℏ/m_ψ)^{0.5} √t scaling follows directly from the random-walk assumption and central limit theorem as stated; m_ψ remains a free parameter whose value can be adjusted to match observed profiles, so the result is not an independent prediction of the simulations beyond the modeling assumption.

    Authors: The scaling is indeed a direct consequence of the random-walk model and central limit theorem, as we state. The simulations are used to confirm that stellar diffusion in ψDM halos proceeds according to this random-walk behavior, with the mass dependence emerging naturally from the soliton properties realized in the runs. While m_ψ is a parameter (constrained elsewhere), the simulations provide concrete, mass-dependent outcomes for the resulting stellar profiles that go beyond the pure analytic assumption. We will revise the abstract and discussion to more explicitly separate the analytic scaling from the simulation-validated results. revision: partial

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained from simulations

full rationale

The paper presents the stellar diffusion, Gaussian (n=0.5) profile, and R_{1/2}(t) scaling as outputs of ψDM simulations driven by intrinsic wave perturbations and soliton random walks. The text explicitly ties the Gaussian form to the central limit theorem and the time scaling to the expected diffusive random-walk behavior, without reducing either to a fitted input or self-citation. No load-bearing self-citations, uniqueness theorems, or ansatzes imported from prior author work appear in the provided manuscript text. The mass-dependent diffusion strength is asserted to emerge from the simulated soliton properties rather than being imposed by construction. This satisfies the criteria for an independent, simulation-backed result.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on the ψDM framework, the random-walk assumption for stellar motion, and the particle mass m_ψ as the scale-setting parameter; no independent evidence for the diffusion mechanism is supplied beyond the simulation itself.

free parameters (1)
  • m_ψ
    Mass of the wave dark matter particle that controls the diffusion length scale in the R_{1/2}(t) relation.
axioms (2)
  • domain assumption Wave perturbations intrinsic to ψDM drive a random walk of stars inside halos
    Invoked as the physical driver of diffusion in the abstract.
  • standard math Central limit theorem produces Gaussian stellar profiles from the random walk
    Used to justify the Sersic n=0.5 result.

pith-pipeline@v0.9.1-grok · 5849 in / 1578 out tokens · 46748 ms · 2026-06-30T05:11:31.375121+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

50 extracted references

  1. [1]

    Amorisco, N. C. & Loeb, A. 2016, MNRAS, 459, L51

    2016

  2. [2]

    N., Smoot, G

    Broadhurst, T., de Martino, I., Luu, H. N., Smoot, G. F., & Tye, S. H. H. 2020, Phys. Rev. D., 101

    2020

  3. [3]

    K., Kereš, D., Wetzel, A., et al

    Chan, T. K., Kereš, D., Wetzel, A., et al. 2018, MNRAS, 478, 906

    2018

  4. [4]

    P., et al

    Chiti, A., Frebel, A., Ji, A. P., et al. 2023, AJ, 165, 55

    2023

  5. [5]

    D., et al

    Chiti, A., Frebel, A., Simon, J. D., et al. 2021, Nature Astronomy, 5, 392

    2021

  6. [6]

    2016, Astro

    Collaboration, P. 2016, Astro. & Astrophys., 594, 63

    2016

  7. [7]

    Collins, M. L. M., Read, J. I., Ibata, R. A., et al. 2021, ArXiv

    2021

  8. [8]

    Collins, M. L. M., Tollerud, E. J., Rich, R., et al. 2019, Mon. Not. R. Astron. Soc., 491, 3496

    2019

  9. [9]

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

    Conn, B. C., Jerjen, H., Kim, D., & Schirmer, M. 2018, ApJ, 852, 68 de Blok, W. J. G. 2010, Advances in Astronomy, 2010, 789293 Di Cintio, A., Brook, C. B., Dutton, A. A., et al. 2017, MNRAS, 466, L1

    2018

  10. [10]

    S., et al

    Drlica-Wagner, A., Bechtol, K., Rykoff, E. S., et al. 2015, ApJ, 813, 109

    2015

  11. [11]

    S., Ferré-Mateu, A., Forbes, D

    Gannon, J. S., Ferré-Mateu, A., Forbes, D. A., et al. 2024, MNRAS, 531, 1856

    2024

  12. [12]

    S., Forbes, D

    Gannon, J. S., Forbes, D. A., Romanowsky, A. J., et al. 2022, MNRAS, 510, 946

    2022

  13. [13]

    2000, Phys

    Hu, W., Barkana, R., & Gruzinov, A. 2000, Phys. Rev. Lett., 85, 1158

    2000

  14. [14]

    P., Tremaine, S., & Witten, E

    Hui, L., Ostriker, J. P., Tremaine, S., & Witten, E. 2017, Phys. Rev. D, 95, 043541

    2017

  15. [15]

    R., Romanowsky, A

    Janssens, S. R., Romanowsky, A. J., Abraham, R., et al. 2022, MNRAS, 517, 858

    2022

  16. [16]

    D., & Cohen, J

    Kirby, E, N., Simon, J. D., & Cohen, J. G. 2015, Astrophys. J., 810, 12

    2015

  17. [17]

    V ., & Valenzuela, O

    Klypin, A., Kravtsov, A. V ., & Valenzuela, O. Prada, F. 1999, Astrophys. J., 522, 82

    1999

  18. [18]

    Li, X., Hui, L., & Yavetz, T. D. 2021, Phys. Rev. D, 103, 023508

    2021

  19. [19]

    S., et al

    Liao, S., Gao, L., Frenk, C. S., et al. 2019, MNRAS, 490, 5182

    2019

  20. [20]

    2018, Phys

    Lin, S.-C., Schive, H.-Y ., Wong, S.-K., & Chiueh, T. 2018, Phys. Rev. D, 97, 103523

    2018

  21. [21]

    & Ishak, M

    Lin, W. & Ishak, M. 2016, Journal of Cosmo. and Astroparticle Phy

    2016

  22. [22]

    1927, Zeitschrift fur Physik, 40, 322

    Madelung, E. 1927, Zeitschrift fur Physik, 40, 322

    1927

  23. [23]

    Marsh, D. J. E. & Silk, J. 2014, Mon. Not. R. Astron. Soc., 437, 2652

    2014

  24. [24]

    F., Nidever, D

    Martin, N. F., Nidever, D. L., Besla, G., et al. 2015, ApJ, 804, L5

    2015

  25. [25]

    & Springel, V

    May, S. & Springel, V . 2021, Mon. Not. R. Astron. Soc., Advanced access

    2021

  26. [26]

    2019, Phys

    Mocz, P., Fialkov, A., V ogelsberger, M., et al. 2019, Phys. Rev. Lett., 123, 141301

    2019

  27. [27]

    H., et al

    Mocz, P., V ogelsberger, M., Robles, V . H., et al. 2017, Mon. Not. R. Astron. Soc., 471, 4559

    2017

  28. [28]

    Mocz, P. et al. 2020, Mon. Not. R. Astron. Soc., 494, 2027

    2020

  29. [29]

    1994, Nature, 370, 629

    Moore, B. 1994, Nature, 370, 629

    1994

  30. [30]

    Naoz, S., Yoshida, N., & Gnedin, N. Y . 2012, ApJ, 747, 128

    2012

  31. [31]

    F., Frenk, C

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

    1996

  32. [32]

    & Loeb, A

    Padmanabhan, H. & Loeb, A. 2021, Astron. Astrophys., 646, 4

    2021

  33. [33]

    2024, ApJ, 963, 54

    Pandya, V ., Zhang, H., Huertas-Company, M., et al. 2024, ApJ, 963, 54

    2024

  34. [34]

    2020, Mon

    Pozo, A., Broadhurst, T., de Martino, I., et al. 2020, Mon. Not. R. Astron. Soc., 504, 2868

    2020

  35. [35]

    2021, MNRAS, 504, 2868

    Pozo, A., Broadhurst, T., de Martino, I., et al. 2021, MNRAS, 504, 2868

    2021

  36. [36]

    2022, MNRAS, 515, 2624

    Pozo, A., Broadhurst, T., Emami, R., & Smoot, G. 2022, MNRAS, 515, 2624

    2022

  37. [37]

    H., Bullock, J

    Robles, V . H., Bullock, J. S., & Boylan-Kolchin, M. 2019, MNRAS, 483, 289

    2019

  38. [38]

    E., Tempel, E., Puzia, T

    Rong, Y ., Mancera Piña, P. E., Tempel, E., Puzia, T. H., & De Rijcke, S. 2020, MNRAS, 498, L72

    2020

  39. [39]

    & Loeb, A

    Safarzadeh, M. & Loeb, A. 2021, Arxiv

    2021

  40. [40]

    S., Okamoto, S., Yagi, M., et al

    Sato, K. S., Okamoto, S., Yagi, M., et al. 2025, ApJ, 993, L7

    2025

  41. [41]

    2020, Phys

    Schive, H.-Y ., Chiueh, T., & Broadhurst, T. 2020, Phys. Rev. Lett., 124, 201301

    2020

  42. [42]

    2016, ApJ, 818, 89

    Schive, H.-Y ., Chiueh, T., Broadhurst, T., & Huang, K.-W. 2016, ApJ, 818, 89

    2016

  43. [43]

    C., & Engels, J

    Schwabe, B., Niemeyer, J. C., & Engels, J. F. 2016, Phys. Rev. D., 94

    2016

  44. [44]

    2020, in American Astronomical So- ciety Meeting Abstracts, V ol

    Tremmel, M., Wright, A., Brooks, A., et al. 2020, in American Astronomical So- ciety Meeting Abstracts, V ol. 235, American Astronomical Society Meeting Abstracts #235, 316.02 van Dokkum, P., Danieli, S., Cohen, Y ., et al. 2018, Nature, 555, 629

    2020

  45. [45]

    C., & Schwabe, B

    Veltmaat, J., Niemeyer, J. C., & Schwabe, B. 2018, Phys. Rev. D., 98

    2018

  46. [46]

    J., Brodie, J., et al

    Villaume, A., Romanowsky, A. J., Brodie, J., et al. 2022, ApJ, 924, 32

    2022

  47. [47]

    Widrow, L. M. & Kaiser, N. 1993, Astrophys. J., 416, L71

    1993

  48. [48]

    & Chiueh, T

    Woo, T.-P. & Chiueh, T. 2009, ApJ, 697, 850

    2009

  49. [49]

    2023, ApJS, 267, 27

    Zaritsky, D., Donnerstein, R., Dey, A., et al. 2023, ApJS, 267, 27

    2023

  50. [50]

    2022, ApJS, 261, 11 Article number, page xiv of xiv

    Zaritsky, D., Donnerstein, R., Karunakaran, A., et al. 2022, ApJS, 261, 11 Article number, page xiv of xiv

    2022