REVIEW 4 major objections 4 minor 39 references
Probing the Dark Sector using the Fornax Satellite
T0 review · 4 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read A long-range dark-sector force stronger than roughly 40 percent of gravity would have stripped most of Fornax's stars, contradicting its observed surface brightness and velocity dispersion; forces below about 20 percent are consistent.
desk verdict Fornax gives a plausible new upper bound on a dark-sector fifth force, but the exclusion threshold is poorly pinned down and the abstract overstates it. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the fifth-force law φ(r) = -β G m² / r exp(-r/r_s) between dark-matter particles, with screening length r_s larger than the system so that the extra acceleration is F5 = β F_gravity. Only dark matter feels this force, so the satellite's dark halo is pulled more strongly than its stars, producing differential tidal stripping. The numerical machinery is a hybrid N-body code: stars move under ordinary gravity, dark-matter particles also feel the β-scaled force, and the Milky Way halo is live so that dynamical friction and the halo's β-dependent response are included; the disk and bulge are static analytic potentials. An analytic companion computes the effective potential i
What would settle it
A deep, wide-field search for diffuse stellar debris around Fornax would settle the stripping prediction: high-β runs eject most stars into a broad, asymmetric envelope, while β ≈ 0 leaves the stellar body compact. If no such envelope appears at the predicted surface brightness, the strong-stripping scenarios are excluded; if it is found, the β constraint tightens.
Extended reading notes
Core claim
The central discovery is a differential tidal effect: a fifth force between dark-matter particles amplifies the effective gravity felt by the dark halo but not by the stars, so the two components of a satellite are pulled apart. In the simulations, the dark halo's mass loss is nearly independent of β, but the stellar mass loss rises sharply with β: at β = 1 essentially all stars are ejected after one pericenter passage; at β = 0.6 about half survive the first passage and the rest are stripped on subsequent orbits; and at β = 0.2 stripping is mild. Comparing the simulated surface-brightness and line-of-sight velocity-dispersion profiles with Fornax data, the paper concludes that β ≳ 0.4 is ex
Load-bearing premise
The adopted Milky Way mass model and the reconstructed orbit of Fornax are accurate; if the Milky Way halo is lighter or Fornax's pericenter is more distant, tidal stripping weakens and larger β values could still match the observations.
Editorial extensions
If this is right
- Fornax alone places an upper bound β ≲ 0.4 (with β ≳ 0.6 clearly ruled out by the simulations) on a long-range dark-sector force whose screening length exceeds the Milky Way halo, under the adopted host model and orbit.
- Because Fornax is distant and on a mildly eccentric orbit, closer-in satellites experience stronger tides, so the same mechanism should produce even tighter limits, or clearer signatures, for other dwarf spheroidals.
- The stellar-to-dark-matter mass-loss ratio becomes a direct observable: high-β runs strip stars far more efficiently than dark matter, so measuring bound stellar mass relative to inferred halo mass tests the force.
- Tidal-stream asymmetry is a complementary diagnostic: a fifth force makes stellar streams asymmetric, while at very high β a surviving self-bound stellar core can leave a symmetric stream, as previously proposed for Sagittarius.
Reading between the lines
- The quoted β limits are tied to the adopted Milky Way mass model and reconstructed orbit; a lighter halo or a more distant pericenter would push the exclusion threshold upward, so the bounds should be read as model-dependent rather than universal.
- A joint analysis of several dwarf spheroidals, each weighted by its orbit and internal structure, could turn this single-galaxy bound into a population-level constraint and potentially reach β below 0.2.
- Deep imaging that searches for a faint, asymmetric stellar envelope around Fornax would provide a direct test: high-β evolution predicts such debris, while β ≈ 0 leaves the stellar body compact.
- The infinite-screening-length assumption is the simplest case; a fifth force with screening length comparable to the Milky Way size would change the orbital phase dependence and could alter the derived bound.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper investigates the observable consequences of a long-range fifth force acting only between dark matter (DM) particles, quantified by the parameter β relative to Newtonian gravity, for the Fornax dwarf spheroidal galaxy. Using GADGET2-based N-body simulations with a live Milky Way DM halo, static disk and bulge, and a Fornax-like two-component satellite, the authors evolve several β values and compare the resulting stellar mass, surface-brightness profile, and line-of-sight velocity-dispersion profile with observations. They report that β ≳ 0.6 strips and unbinds most of the stellar component, while β ≲ 0.2 remains broadly consistent with data. An analytic tidal-stripping model and a critical-β estimate are presented in appendices.
Significance. If the result holds, the paper provides a new astrophysical probe of dark-sector fifth forces using a relatively distant, low-eccentricity dSph, complementing earlier work on Sagittarius. The numerical setup is described in reasonable detail, includes dynamical friction through a live MW halo, and the qualitative trend—stronger β leading to enhanced stellar stripping—is physically transparent. The main value is in setting a concrete, falsifiable constraint on β from the survival of Fornax’s stellar component. However, as presented the constraint is not yet quantitative, and the central exclusion threshold is not robustly established because the initial conditions are re-fit for each β and the profile comparisons are visual.
major comments (4)
- [§II.B, Table II] The initial Fornax models are re-fit for each β: Table II shows aDM decreasing from 5.1 kpc at β=0 to 2.7 kpc at β=1.0, with a⋆ also changing. Thus the simulations do not hold a fixed pre-infall progenitor while varying the fifth-force strength; each β is assigned a different, more concentrated initial DM halo tuned to present-day Fornax. The central exclusion (§IV) therefore conflates the effect of β with the effect of the chosen initial conditions. A forward test with a single plausible pre-infall model (e.g., the β=0 relaxed profile) evolved under different β, plus a sensitivity scan over initial stellar mass and DM concentration, is needed before the bound can be considered robust.
- [§III, Fig. 4] The claim that a given β 'fails to match' or 'remains broadly consistent' with observations is based on visual inspection of the surface-brightness and velocity-dispersion profiles. No χ², likelihood, or other quantitative metric is computed, and the plotted observational data are not accompanied by an explicit treatment of uncertainties in the comparison. Since the paper presents a constraint on β, a statistical comparison—e.g., a profile likelihood over the photometric and kinematic data with a stated covariance model—is required to define the exclusion threshold and to translate the disagreement seen at β=0.4–0.6 into a confidence level.
- [Abstract vs §IV] The abstract states that β ≳ 0.6 strips and unbinds the bulk of the stars and that β ≲ 0.2 is broadly consistent, while §IV concludes with 'exclusion of β ≳ 0.4 (not presented within the figures)'. These thresholds are materially different. Table IV shows that at β=0.4 the bound stellar mass is 0.9×10^7 M⊙, about 37% of the initial 2.4×10^7 M⊙, and the central surface density is 0.4×10^7 M⊙ kpc^-2, which is not 'negligible'. The authors should reconcile the abstract, text, and tabulated values, and state explicitly which observable and criterion define the exclusion.
- [Appendix C] The analytic βcrit estimate is explicitly post-hoc: the DM distribution is taken from the β=0.2 simulation to obtain βcrit≃0.3, and from the β=1 simulation at t=4 Gyr to obtain βcrit≈0.9. As acknowledged, this is a consistency check rather than an independent derivation, but the wording in §IV ('as a complementary check ... reproduces the qualitative trends') risks overstating its role. It should be clearly labeled as non-predictive and should not be used to support the quantitative exclusion without a forward calculation from fixed initial conditions.
minor comments (4)
- [§II.C] The backward-time integration used to set the orbit is not described in enough detail. Specify which gravitational potential is used for the backward integration (static β=0 MW model, or the β-dependent live halo) and whether the resulting pericenter distance changes with β.
- [Table IV] The M⋆(rt) row has only four entries while the table has five β columns. If β=1 leaves no bound stars, state this explicitly or use a dash; otherwise the table is confusing.
- [Fig. 4] The right panel text says the simulated velocity-dispersion profiles are compared with the 'observed mean line-of-sight value'; please clarify whether the plotted black points are binned data with error bars and whether the normalization convention is the same as for the simulated profiles.
- [§II.A.2] Table I lists 'ainitial' and 'γinitial' for the Dehnen halo, but the text says the scale radius and slope are 'best-fit values after the relaxation'. Clarify whether the reported values are the initial or relaxed parameters, and whether the relaxed profile is refit for each β.
Circularity Check
Main β constraint is a forward N-body scan; only the Appendix C analytic check is a post-hoc consistency check using simulation output as input.
-
other
[Section IV (Discussion), paragraph reporting analytic check; Appendix C]
"Using the DM distribution from the β=0.2 simulation as input, the analytic estimate yields a critical βcrit≃0.3 above which stars and DM are expected to segregate (see appendix C), consistent with the absence of full segregation in that run."
The analytic estimate takes the DM distribution from the very β=0.2 simulation whose outcome it is used to corroborate, so the agreement is built from the run being explained rather than from independent input. It is a post-hoc consistency check, not a prediction. It does not feed back into the N-body exclusion bound, so the central claim retains independent content.
full rationale
The central exclusion bound (β≳0.4 in Section IV, β≳0.6 in the abstract) is obtained by direct N-body simulations in which β is a scanned input parameter; the simulated stellar mass, surface-brightness, and velocity-dispersion profiles are then compared with Fornax observations. No parameter is fitted to the final profiles to produce the bound, so the main claim does not reduce to its inputs by construction. The initial conditions are relaxed in isolation for each β and matched to Fornax's observed structural parameters (Table II); this makes the low-β consistency partly by construction, but it does not force the high-β disruption that drives the exclusion, and the sensitivity to assumed initial structure is an acknowledged modeling-systematic uncertainty rather than a circular step. The analytic βcrit estimate in Appendix C is explicitly a complementary check: it uses the DM distribution from the β=0.2 simulation as input and finds consistency with that same run, which is self-referential but not load-bearing. Self-citations ([9], [13], [36]) supply context and some MW parameters that are also supported by external references; no uniqueness theorem or ansatz is imported solely from the authors' prior work. The paper candidly lists limitations (MW potential, orbit reconstruction, assembly history, neglected physics), which are correctness risks, not circularity. Overall, no significant circularity; score 1 reflects only the minor post-hoc analytic check.
Assumptions & free parameters
free parameters (3)
- β (fifth-force amplitude) =
scanned over 0, 0.2, 0.4, 0.6, 1.0
- Fornax DM Dehnen scale radius (aDM) =
5.1, 4.0, 3.6, 3.3, 2.7 kpc for β=0,0.2,0.4,0.6,1.0
- Fornax stellar Dehnen scale radius (a*) =
0.4 kpc (β=0), 0.3 kpc (others)
assumptions (5)
- domain assumption The fifth force acts only between DM particles and is F5 = β Fgravity (with rs much larger than the system)
- domain assumption The MW gravitational field is described by the adopted disk, bulge, and Dehnen halo parameters (Table I)
- domain assumption Fornax and the MW halo are initially in dynamical equilibrium and can be represented by Dehnen profiles
- domain assumption The orbit of Fornax is reconstructed from Gaia astrometry and the adopted MW potential via backward integration
- domain assumption Standard Newtonian dynamics for stellar particles; the fifth force affects only DM particles
Cite this review
Pith. "Pith review of Probing the Dark Sector using the Fornax Satellite." pith.science (2026). https://pith.science/paper/RH27CIWX
@misc{pith2026250901713,
author = {Pith},
title = {Pith review of: Probing the Dark Sector using the Fornax Satellite},
year = {2026},
howpublished = {\url{https://pith.science/paper/RH27CIWX}},
note = {Machine review of arXiv:2509.01713}
}
abstract
A long-range force acting only in the dark sector can displace the stellar component of satellite galaxies relative to their dark matter (DM) halos, thereby \emph{breaking} the weak equivalence principle (WEP) between DM and baryons. We investigate observational signatures of such WEP breaking using $N$-body simulations of a Fornax-like Milky Way (MW) satellite, implementing a fifth force between DM particles of amplitude $\beta$ relative to Newtonian gravity (with a screening length much larger than the MW halo size). We find that $\beta \gtrsim 0.6$ strips and unbinds the bulk of the stellar component, leaving at most a negligible bound remnant that fails to match the observed stellar content, surface-brightness, and line-of-sight velocity-dispersion profiles of Fornax. By contrast, $\beta \lesssim 0.2$ yields only mild stellar stripping and remains broadly consistent with current photometric and kinematic constraints within our modeling assumptions.
Figures
Reference graph
Works this paper leans on
-
[1]
Static Potential Components A baryonic Miyamoto-Nagai [21] gravitational potential for the MW disk is adopted, Φdisk =− GMd√ r2 ⊥ + ( ad + √ z2 +d2)2 , (3) where z and r⊥ are cylindrical coordinates, withz per- pendicular to the disk. The stellar bulge potential is taken as a spherical Hern- quist [22] in the form of, Φbulge =−GMb r +c . (4) The parameter...
-
[2]
The MW Halo The MW halo is represented as an equilibriumN-body realization of a Dehnen density profile [23], ρhalo(r) = (3−γ)MDM 4π a rγ(r +a)4−γ . (5) Both the scale radiusa and γ, which control the slope of the density profile near the center, are given as the best-fit values after the relaxation (see table I). The initial conditions were generated with...
-
[3]
P. J. E. Peebles and B. Ratra, Rev. Mod. Phys.75, 559 (2003)
2003
-
[4]
A. G. Riess, A. V. Filippenko, P. Challis, A. Clocchiatti, A. Diercks, Garnavich,et al., The Astronomical Journal 116, 1009 (1998), arXiv:astro-ph/9805201 [astro-ph]
arXiv 1998
-
[5]
Planck Collaboration, N. Aghanim, Y. Akrami, M. Ash- down, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, Bartolo,et al., A & A641, A6 (2020), arXiv:1807.06209 [astro-ph.CO]
arXiv 2020
-
[6]
B. Carr, F. Kühnel, and M. Sandstad, Phys. Rev. D94, 083504 (2016), arXiv:1607.06077 [astro-ph.CO]
arXiv 2016
- [7]
-
[8]
DESI Collaboration, arXiv e-prints , arXiv:2503.14738 (2025), arXiv:2503.14738 [astro-ph.CO]
arXiv 2025
Show all 39 references
-
[9]
Colladay and V
D. Colladay and V. A. Kostelecký, Physical Review D58, 116002 (1998)
1998
-
[10]
Jacobson and D
T. Jacobson and D. Mattingly, Phys. Rev. D64, 024028 (2001), arXiv:gr-qc/0007031 [gr-qc]
2001 arXiv
-
[11]
J. A. Keselman, A. Nusser, and P. J. E. Peebles, Phys. Rev. D 80, 063517 (2009)
2009
-
[12]
J. A. Frieman and B. A. Gradwohl, Science260, 1441 (1993)
1993
-
[13]
W. A. Hellwing, S. R. Knollmann, and A. Knebe, Monthly Notices of the Royal Astronomical Society: Letters408, L104 (2010), https://academic.oup.com/mnrasl/article- pdf/408/1/L104/54672309/mnrasl_408_1_l104.pdf
2010
-
[14]
W. A. Hellwing, B. Li, C. S. Frenk, and S. Cole, Monthly Notices of the Royal Astronomical Society435, 2806 (2013), https://academic.oup.com/mnras/article- pdf/435/4/2806/3367715/stt1430.pdf
2013
-
[15]
Nusser, S
A. Nusser, S. S. Gubser, and P. J. E. Peebles, Phys. Rev. D 71, 083505 (2005)
2005
-
[16]
Kesden and M
M. Kesden and M. Kamionkowski, Phys. Rev. D 74, 083007 (2006)
2006
- [17]
-
[18]
Mateo, Annual Review of Astronomy and Astrophysics 36, 435 (1998)
M. Mateo, Annual Review of Astronomy and Astrophysics 36, 435 (1998)
1998
-
[19]
S. S. Gubser and P. J. E. Peebles, Phys. Rev. D70, 123510 (2004)
2004
-
[20]
K., Battaglia, G., Pawlowski, M
Fritz, T. K., Battaglia, G., Pawlowski, M. S., Kallivayalil, N., van der Marel, R., Sohn, S. T., Brook, C., and Besla, G., Astronomy & Astrophysics619, A103 (2018)
2018
-
[21]
M. G. Walker and J. Peñarrubia, The Astrophysical Jour- nal 742, 20 (2011), arXiv:1108.2404 [astro-ph.CO]
2011 arXiv
-
[22]
N. C. Amorisco, A. Agnello, and N. W. Evans, Monthly Notices of the Royal Astronomical Society: Letters429, L89 (2012), https://academic.oup.com/mnrasl/article- pdf/429/1/L89/54662740/mnrasl_429_1_l89.pdf
2012
-
[23]
Miyamoto and R
M. Miyamoto and R. Nagai, Publications of the Astro- nomical Society of Japan27, 533 (1975)
1975
-
[24]
Hernquist, The Astrophysical Journal356, 359 (1990)
L. Hernquist, The Astrophysical Journal356, 359 (1990)
1990
-
[25]
Dehnen, Monthly Notices of the Royal Astronomical Society 265, 250 (1993)
W. Dehnen, Monthly Notices of the Royal Astronomical Society 265, 250 (1993)
1993
-
[26]
Ruggiero, galstep, https://github.com/ruggiero/ galstep (2014)
R. Ruggiero, galstep, https://github.com/ruggiero/ galstep (2014)
2014
-
[27]
P. J. McMillan, Monthly Notices of the Royal Astronomical Society 414, 2446 (2011), https://academic.oup.com/mnras/article- pdf/414/3/2446/3525100/mnras0414-2446.pdf. 9
2011
-
[28]
Schönrich, J
R. Schönrich, J. Binney, and W. Dehnen, Monthly Notices of the Royal Astronomical Society 403, 1829 (2010), https://academic.oup.com/mnras/article- pdf/403/4/1829/18575828/mnras0403-1829.pdf
2010
-
[29]
Borukhovetskaya, R
A. Borukhovetskaya, R. Errani, J. F. Navarro, A. Fattahi, and I. Santos-Santos, Monthly No- tices of the Royal Astronomical Society 509, 5330 (2021), https://academic.oup.com/mnras/article- pdf/509/4/5330/41745368/stab2912.pdf
2021
-
[30]
A. W. McConnachie, The Astronomical Journal144, 4 (2012)
2012
-
[31]
Pietrzyński, M
G. Pietrzyński, M. Górski, W. Gieren, V. D. Ivanov, F. Bresolin, and R.-P. Kudritzki, The Astronomical Jour- nal 138, 459 (2009)
2009
-
[32]
Croker, Computer Physics Communications207, 478 (2016)
K. Croker, Computer Physics Communications207, 478 (2016)
2016
-
[33]
com/giulianoiorio/Gadget2-static_potential (2015)
Giuliano, Static potential extension, https://github. com/giulianoiorio/Gadget2-static_potential (2015)
2015
-
[34]
M. G. Coleman, G. S. D. Costa, J. Bland-Hawthorn, and K. C. Freeman, The Astronomical Journal 129, 1443 (2005)
2005
-
[35]
M. G. Walker, M. Mateo, E. W. Olszewski, O. Y. Gnedin, X. Wang, B. Sen, and M. Woodroofe, The Astrophysical Journal 667, L53 (2007), arXiv:0708.0010 [astro-ph]
2007 arXiv
-
[36]
M. G. Walker, M. Mateo, E. W. Olszewski, J. Peñarrubia, N. W. Evans, and G. Gilmore, The Astronomical Journal 704, 1274 (2009), arXiv:0906.0341 [astro-ph.CO]
2009 arXiv
-
[37]
D. Blas, M. M. Ivanov, and S. Sibiryakov, JCAP2012 (10), 057, arXiv:1209.0464 [astro-ph.CO]
-
[38]
Bettoni, A
D. Bettoni, A. Nusser, D. Blas, and S. Sibiryakov, JCAP 2017 (5), 024, arXiv:1702.07726 [astro-ph.CO]
2017 arXiv
-
[39]
King, Astronomical Journal67, 471 (1962)
I. King, Astronomical Journal67, 471 (1962)
1962
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.