REVIEW 2 major objections 5 minor 1 cited by
A damping term in ultralight dark matter reduces dynamical friction enough to resolve the Fornax timing problem for globular cluster GC3.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-04 00:34 UTC pith:TI642QBS
load-bearing objection The damped friction formula is real and reusable; the Fornax "resolution" is a model-dependent claim resting on an unexamined damping term in the noninteracting limit. the 2 major comments →
Damping of dynamical friction force in self-interacting ultralight dark matter and Fornax timing problem
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The paper's central claim is that the phenomenological damping term in the generalized Gross-Pitaevskii–Poisson equation, with strength ξ = 2kBT/ℏ, essentially decreases the magnitude of the dynamical friction force acting on stars, globular clusters, or dwarf galaxies moving in an ultralight dark matter halo. Concretely, the authors derive an analytic expression for the tangential friction force on a circularly orbiting point mass (Eq. 16) and show that the damping introduces an exponential suppression e^{-ξτ/2} in the force integral, reducing the force compared to the undamped case. Applying this to the Fornax dwarf galaxy, they find that the infall time of globular cluster GC3 can reach t
What carries the argument
The key mechanism is the damping term in the generalized Gross-Pitaevskii equation (Eq. 1), with coefficient ξ = 2kBT/ℏ set by the halo temperature, identified with the stellar velocity dispersion through kBT = mσv²/2. In the linearized density-perturbation equation (Eq. 4) this term appears as ξ∂tα, shifting the response poles (Eq. 11) and producing an e^{-ξτ/2} factor in the friction-force integral (Eq. 12). The paper combines this with two analytical soliton density profiles—a Gaussian for weak or no self-interaction and a Thomas-Fermi profile for strong repulsive self-interaction—matched to an isothermal NFW envelope, to compute the force and the inspiral trajectory of GC3.
Load-bearing premise
The result hinges on the phenomenological damping term in the generalized Gross-Pitaevskii equation, with its strength fixed by equating kBT to mσv²/2 using the observed stellar velocity dispersion; if ultralight dark matter is not actually dissipative in this way, the reduced dynamical friction—and the proposed resolution of the Fornax timing problem—disappears. It also assumes the homogeneous-medium circular-orbit force applies locally during the inspiral.
What would settle it
An N-body or wave simulation of a self-interacting ULDM halo that includes the same damping coefficient and finds that GC3-like clusters still sink in a few Gyr—or an observational determination that the Fornax globular clusters' orbital decay times are much shorter than claimed—would falsify the proposed resolution. More directly, a measurement or bound on the ULDM temperature/dissipation that rules out ξ = 2kBT/ℏ would remove the effect.
If this is right
- The damping-induced reduction of dynamical friction applies not only to GC3 but to any compact object orbiting inside an ULDM halo, potentially altering predicted merger rates and orbital decays in dwarf galaxies.
- Fornax's timing problem can be resolved in nearly noninteracting (fuzzy) ULDM with m ≈ 3.09e-22 eV, provided GC3 formed at a galactocentric distance of about 1.5 kpc rather than much closer to the center.
- In strongly self-interacting ULDM, the timing problem is resolved for m ≳ 2.7e-21 eV at r0 = 1.5 kpc, and even at r0 = 1 kpc if the boson mass is sufficiently large.
- The analytic formula (Eq. 16) provides a ready tool for estimating dynamical friction in other self-interacting ULDM systems where the damping term is relevant.
- Because the damping term also changes the force's radial dependence, the predicted rotation curves and infall times could be used to distinguish damped ULDM from ordinary CDM cores.
Where Pith is reading between the lines
- If the damping term is real and as strong as assumed, similar suppression should affect the orbital decay of satellite galaxies and black-hole binaries in ULDM halos; testing these systems could independently confirm or bound the dissipation strength.
- The value of ξ is fixed by observational velocity dispersion via a thermodynamic relation; this is a phenomenological input rather than a derivation from a fundamental ULDM Lagrangian, so the central result hinges on that identification being valid.
- The homogeneous-medium, circular-orbit treatment likely underestimates the damping effect for eccentric orbits; extending the calculation to eccentric trajectories might either strengthen or weaken the claimed resolution, a testable modification.
- One could search for the damped wake's signature—reduced amplitude and altered phase of the density perturbation behind the perturber—in future high-resolution ULDM simulations or gravitational lensing observations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the dynamical friction force acting on globular cluster GC3 in the Fornax dwarf spheroidal, assuming an ultralight bosonic dark matter halo described by a generalized Gross-Pitaevskii-Poisson equation that includes a temperature-dependent damping term. Starting from linear-response theory, the authors derive an analytic expression for the dynamical friction force for a point mass on a circular orbit in a homogeneous ULDM medium in the presence of damping (Eq. 16). Using observationally motivated density profiles (Gaussian soliton core for weak/noninteracting ULDM, Thomas-Fermi core for strong repulsive self-interaction, matched to an NFW envelope), they compute the tangential force and integrate the orbital decay equations for GC3. With the damping term, they find that infall times can reach 12 Gyr for noninteracting ULDM with m≈3.09e-22 eV if the initial radius is 1.5 kpc, and for strongly interacting ULDM with m≳2.7e-21 eV at the same initial radius, resolving the Fornax timing problem in those cases.
Significance. Should the damping mechanism be physically valid, the paper offers an analytic and transparent route to resolving the Fornax timing problem within ULDM models, and it makes falsifiable predictions in the (m, r0) plane. The derivation of Eq. (16) is internally consistent, and the numerical parameter scan is clearly presented. However, the central claim rests on a phenomenological dissipation term whose microscopic justification is not established in the noninteracting limit, and the local application of a homogeneous-medium circular-orbit force involves uncontrolled approximations. The paper is therefore more a proof-of-principle exploration than a definitive resolution, but its identified parameter regions are concrete and testable.
major comments (2)
- [Sec. 2, Eqs. (1)-(2); Sec. 5, Eq. (46)] The damping term with strength ξ=2k_BT/ℏ is adopted from a finite-temperature dissipative BEC model (Ref. [14]) and is not derived from a fundamental action for ULDM. In the noninteracting limit a_s=0, which is the paper's headline Gaussian-regime solution (Sec. 6, Fig. 4b, m=3.09e-22 eV), there is no microscopic justification for this dissipation; one expects no thermal coupling for an ideal scalar field. Since the ξ=0 calculation already gives the known short infall times, the claimed resolution is entirely contingent on this unvetted term. Please provide a derivation or at least a quantitative estimate of ξ in the a_s→0 limit, and discuss why the same ξ applies inside the soliton core, where the equilibrium model takes k_BT effects to be negligible.
- [Sec. 3, Eq. (16); Appendix, Eq. (A.5); Sec. 6, Eq. (47)] The derived force assumes a circular orbit in a homogeneous medium of constant density ρ_DM. It is then applied at each instantaneous radius of a spiraling orbit in the strongly inhomogeneous core-envelope profile of Sec. 4, and the radial component F_r is dropped with the only justification that it is 'much smaller' than the gravitational term. No quantitative estimate is given. This is a load-bearing approximation for the 12-Gyr infall claims. Please test it, e.g., by comparing with a local expansion in ρ'(r)/ρ(r), by including F_r, or by a direct simulation in a spherically symmetric halo.
minor comments (5)
- [Sec. 7, first paragraph] The quoted lower bound 'mlow ≈ 3.09·10^-21 eV' contradicts Eq. (27) and Sec. 6, where the noninteracting mass is 3.09·10^-22 eV. This typo is in a key conclusion and should be corrected.
- [Sec. 5, Eq. (45)] The notation is inconsistent: the dimensionless wavevector should appear as K in the denominator (K^2 \tilde{D}(K)), not as k. Please harmonize the notation.
- [Fig. 3] Panels a) and b) use different y-axis scales (0.6 vs 1.5), which exaggerates the visual suppression by the damping term. Consider a common scale or explicit annotation.
- [Sec. 6, Eq. (47)] The initial angular momentum l(0) is not specified. The integration of Eq. (47) requires an initial tangential velocity; please state explicitly whether the orbit is initially circular, l = m_GC r_0 v_c(r_0).
- [Throughout] Typos: 'Euristic' should be 'Heuristic' (Sec. 1); 'ultra light' should be 'ultralight' (Abstract, Sec. 1). References [7] and [13] appear to refer to the same paper with inconsistent journal volume/year; please reconcile.
Circularity Check
No circular derivation: the claimed infall times are computed from externally fixed halo/damping inputs, not used to fit the 12-Gyr target.
full rationale
The derivation chain is self-contained in the relevant sense. The Fornax halo density is fixed by external observations (rho_obs at r_h, M = 5.8e7 M_sun; Sec. 4), and the damping coefficient xi = 2 k_B T / hbar is set from the observed stellar velocity dispersion sigma_v = 8 km/s via k_B T = m sigma_v^2 / 2 (Eqs. 2 and 46). The infall times in Sec. 6 are computed from the dynamical friction force in Eqs. (44)-(47), and the 12-Gyr GC3 age is not used to calibrate any parameter; the scan over m and r0 is openly exploratory. The only self-citations, Refs. [22,23], concern the Plummer-sphere generalization, which is not used in the point-perturber calculation, so they are contextual rather than load-bearing. The damping term is adopted from the external Chavanis model [14]; its microscopic validity in the as=0 limit is a physics assumption and a correctness risk, but it is a model input, not a circular restatement of the target result. No equation in the paper reduces by construction to the claimed Fornax resolution.
Axiom & Free-Parameter Ledger
free parameters (4)
- ULDM boson mass m =
scanned over 3.09×10^-22 to 3×10^-21 eV
- Initial GC3 orbital radius r0 =
1.0 or 1.5 kpc
- Soliton scale radii (R=581 pc for Gaussian, R_TF=1 kpc for TF) =
581 pc / 1 kpc
- NFW envelope parameters ρ_e and r_e =
ρ_e=2.59×10^3 M_sun/pc^3, r_e≈7.5 pc (Gaussian); ρ_e=2.06×10^3 M_sun/pc^3, r_e=8.4 pc (TF)
axioms (6)
- domain assumption Generalized Gross-Pitaevskii-Poisson equation with logarithmic nonlinearity and damping term (Eq. 1, last term) governs ULDM halo dynamics
- domain assumption Einstein/fluctuation-dissipation relation ξ = 2k_B T/ℏ (Eq. 2)
- domain assumption Halo temperature equated to observed stellar velocity dispersion: k_B T = m σ_v²/2 (Eq. 46)
- domain assumption Linearized density-perturbation equation (Eq. 4) with damping entering as +ξ∂_t α, solved in the homogeneous-medium approximation and applied locally to the inhomogeneous Fornax halo
- domain assumption Soliton density profiles: Gaussian ansatz (Eq. 25) and Thomas-Fermi profile (Eq. 20), and NFW envelope stitched by continuity (Secs. 4.1–4.4) describe the Fornax halo
- ad hoc to paper Circular-orbit dynamical-friction force is used at each instantaneous radius during the inspiraling orbit, and the radial component F_r is dropped
read the original abstract
The dynamics of globular clusters in the Fornax dwarf galaxy pose a challenge for the standard cold dark matter and can be used to test other models of dark matter. We study this dynamics in the context of an ultralight bosonic dark matter model, accounting for the damping term in a generalized Gross-Pitaevskii equation. Employing analytic formulas for the dynamical friction force, the infall time and evolution of globular clusters are compared in the cases with and without the damping term. It is argued that the damping term plays an important role in the Fornax timing problem in ultralight dark matter (ULDM) models. We found that the ULDM model with repulsive self-interaction can solve the Fornax timing problem in the absence of or with very small self-interaction, even if the initial position of the globular cluster is not far from the center of the galaxy. Still, the problem is resolved for strongly interacting repulsive ULDM, even for the most pressing case of globular cluster GC3, if its starting position exceeds 1.5 kpc.
Figures
Forward citations
Cited by 1 Pith paper
-
Vortex State of Ultralight Dark Matter and the Fornax Timing Problem
A vortex state of ultralight dark matter suppresses dynamical friction for co-rotating globular clusters, potentially resolving the Fornax timing problem.
Reference graph
Works this paper leans on
-
[1]
Gratton, A
R. Gratton, A. Bragaglia, E. Carretta, V. D’Orazi, S. Lucatello a nd A. Sollima, What is a globular cluster? an observational perspective , The Astronony and Astrophysics Review 27 (2019) 8
2019
-
[2]
Tremaine, J
S. Tremaine, J. Ostriker and S. Spitzer, The formation of the nuclei of galaxies. I. M31., The Astrophysical Journal 196 (1975) 407
1975
-
[3]
Pace, M.G
A.B. Pace, M.G. Walker, S.E. Koposov, N. Caldwell, M. Mateo, E.W. Ols zewski et al., Spectroscopic confirmation of the sixth globular cluster in the Fornax dwarf spheroidal galaxy , The Astrophysical Journal 923 (2021) 77
2021
-
[4]
Tremaine, The formation of the nuclei of galaxies
S. Tremaine, The formation of the nuclei of galaxies. II - The local group , The Astrophysical Journal 203 (1976) 345
1976
-
[5]
K. Oh, D. Lin and H. Richer, Globular clusters in the Fornax dwarf spheroidal galaxy, The Astrophysical Journal 531 (2000) 727
2000
-
[6]
Meadows, J
N. Meadows, J. Navarro, I. Santos-Santos, A. Benitez-Llamb ay and C. Frenk, Cusp or core? Revisiting the globular cluster timing problem in F ornax, Mon. Not. Roy. Astron. Soc. 491 (2020) 3336
2020
-
[7]
N. Bar, S. Danieli and K. Blum, Dynamical friction in globular cluster-rich ultra-diffuse galaxies: The case of NGC5846-UDG1 , The Astrophysical Journal Letters 932 (2021) L10
2021
-
[8]
Sanchez-Salcedo, J
F. Sanchez-Salcedo, J. Reyes-Iturbide and X. Hernandez, An extensive study of dynamical friction in dwarf galaxies: the role of stars, dar k matter, halo profiles and MOND , Monthly Notices of the Royal Astronomical Society 370 (2006) 1829
2006
-
[9]
D. Cole, W. Dehnen, J. Read and M. Wilkinson, The mass distribution of the Fornax dSph: constraints from its globular cluster distrib ution, Monthly Notices of the Royal Astronomical Society 426 (2012) 1829
2012
-
[10]
Goerdt, B
T. Goerdt, B. Moore, J. Read, J. Stadel and M. Zemp, Does the Fornax dwarf spheroidal have a central cusp or core? , Monthly Notices of the Royal Astronomical Society 368 (2006) 1073
2006
-
[11]
Boldrini, R
P. Boldrini, R. Mohayaee and J. Silk, Fornax globular cluster distributions: implications for the cusp-core problem , Monthly Notices of the Royal Astronomical Society 485 (2019) 2546
2019
-
[12]
Ferreira, Ultra-light dark matter , Astron
E.G.M. Ferreira, Ultra-light dark matter , Astron. Astrophys. Rev. 29 (2021) 7 [2005.03254]
Pith/arXiv arXiv 2021
-
[13]
N. Bar, S. Danieli and K. Blum, Dynamical friction in globular cluster-rich ultra-diffuse galaxies: the case of ngc5846-udg1 , The Astrophysical Journal 932 (2022) L10
2022
-
[14]
Chavanis, Predictive model of BEC dark matter halos with a solitonic co re and an isothermal atmosphere , Physical Review D 100 (2019) 083022
P.-H. Chavanis, Predictive model of BEC dark matter halos with a solitonic co re and an isothermal atmosphere , Physical Review D 100 (2019) 083022 . 16
2019
-
[15]
H.-Y. Schive, M.-H. Liao, T.-P. Woo, S.-K. Wong, T. Chiueh, T. Bro adhurst et al., Understanding the Core-Halo Relation of Quantum Wave Dark M atter from 3D Simulations, Phys. Rev. Lett. 113 (2014) 261302 [1407.7762]
Pith/arXiv arXiv 2014
-
[16]
L. Hui, J.P. Ostriker, S. Tremaine and E. Witten, Ultralight scalars as cosmological dark matter , Phys. Rev. D 95 (2017) 043541 [1610.08297]
Pith/arXiv arXiv 2017
-
[17]
L. Lancaster, C. Giovanetti, P. Mocz, Y. Kahn, M. Lisanti and D.N. Spergel, Dynamical Friction in a Fuzzy Dark Matter Universe , JCAP 01 (2020) 001 [1909.06381]
Pith/arXiv arXiv 2020
-
[18]
T.J.L. de Boer and M. Fraser, Four and one more: The formation history and total mass of globular clusters in the Fornax dSph , Astronomy & Astrophysics 590 (2016) A35 [1510.05642]
Pith/arXiv arXiv 2016
-
[19]
N. Bar, D. Blas, K. Blum and H. Kim, Assessing the Fornax globular cluster timing problem in different models of dark matter , Phys. Rev. D 104 (2021) 043021 [2102.11522]
Pith/arXiv arXiv 2021
-
[20]
Desjacques, A
V. Desjacques, A. Nusser and R. Buehler, Analytic solution to the dynamical friction acting on circularly moving perturbers , The Astrophysical Journal 928 (2022) 64
2022
-
[21]
R. Buehler and V. Desjacques, Dynamical friction in fuzzy dark matter: Circular orbits, Phys. Rev. D 107 (2023) 023516 [2207.13740]
Pith/arXiv arXiv 2023
-
[22]
V.M. Gorkavenko, A.I. Yakimenko, A.O. Zaporozhchenko and E.V . Gorbar, Dynamical friction in ultralight dark matter: Plummer sphe re perspective, Phys. Scripta 100 (2025) 075039 [2412.15428]
Pith/arXiv arXiv 2025
-
[23]
O.V. Barabash, T.V. Gorkavenko, V.M. Gorkavenko, O.M. Teslyk , N.S. Yakovenko, A.O. Zaporozhchenko et al., Analytic calculation of dynamical friction for Plummer sphere in ultralight dark matter , 2504.06448
-
[24]
L. Berezhiani, G. Cintia, V. De Luca and J. Khoury, Dynamical friction in dark matter superfluids: The evolution of black hole binaries , JCAP 06 (2024) 024 [2311.07672]
Pith/arXiv arXiv 2024
-
[25]
Jardel and K
J.R. Jardel and K. Gebhardt, The dark matter density profile of the Fornax dwarf , The Astrophysical Journal 746 (2012) 89
2012
-
[26]
McConnachie, The observed properties of dwarf galaxies in and around the local group, The Astronomical Journal 144 (2012) 4
A.W. McConnachie, The observed properties of dwarf galaxies in and around the local group, The Astronomical Journal 144 (2012) 4
2012
-
[27]
P.-H. Chavanis, Mass-radius relation of Newtonian self-gravitating Bose- Einstein condensates with short-range interactions: I. Analytical results, Phys. Rev. D 84 (2011) 043531 [1103.2050]
Pith/arXiv arXiv 2011
-
[28]
Chavanis, Self-gravitating Bose-Einstein condensates , Quantum Aspects of Black Holes (2015) 151
P.-H. Chavanis, Self-gravitating Bose-Einstein condensates , Quantum Aspects of Black Holes (2015) 151 . 17
2015
-
[29]
Chavanis, Mass-radius relation of self-gravitating Bose-Einstein c ondensates with a central black hole , The European Physical Journal Plus 134 (2019) 352
P.-H. Chavanis, Mass-radius relation of self-gravitating Bose-Einstein c ondensates with a central black hole , The European Physical Journal Plus 134 (2019) 352
2019
-
[30]
J.F. Navarro, C.S. Frenk and S.D.M. White, A Universal density profile from hierarchical clustering, Astrophys. J. 490 (1997) 493 [astro-ph/9611107]
Pith/arXiv arXiv 1997
-
[31]
Binney and S
J. Binney and S. Tremaine, Galactic dynamics , vol. 13, Princeton university press (2011)
2011
-
[32]
M.G. Walker, M. Mateo, E.W. Olszewski, O.Y. Gnedin, X. Wang, B. Se n et al., Velocity Dispersion Profiles of Seven Dwarf Spheroidal Gala xies, Astrophys. J. Lett. 667 (2007) L53 [0708.0010]. 18
Pith/arXiv arXiv 2007
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.