Pith. sign in

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 →

arxiv 2509.01713 v1 pith:RH27CIWX submitted 2025-09-01 astro-ph.CO

classification astro-ph.CO
keywords darksectorfifthforceweakequivalenceprinciplebreakingdwarfspheroidalgalaxytidalstrippingFornaxN-bodysimulationssurfacebrightnessvelocitydispersion
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper sets out to test whether a long-range fifth force that couples only to dark matter—a form of weak-equivalence-principle breaking—is compatible with the Milky Way dwarf satellite Fornax. Using N-body simulations of a Fornax-like galaxy orbiting a live Milky Way halo, it shows that such a force with strength β greater than about 0.6 times gravity rips most of the stars out of the dark-matter halo, leaving a bound stellar remnant too small to match Fornax's surface-brightness and velocity-dispersion profiles. The authors therefore exclude β ≳ 0.4 under their adopted Milky Way potential and reconstructed orbit, while β ≲ 0.2 remains broadly consistent with observations. The result matters because it turns an otherwise hard-to-probe dark-sector coupling into a testable prediction using a single well-studied nearby dwarf galaxy.

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.

Watch

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

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

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)
  1. [§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.
  2. [§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.
  3. [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.
  4. [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)
  1. [§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 β.
  2. [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.
  3. [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.
  4. [§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

1 steps flagged · score 1.0 of 10

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.

  1. 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 3 free parameters · 5 assumptions · 0 invented entities

No new particles or forces are introduced; the fifth force is a known hypothesis. The central result depends on the fitted Dehnen scale radii of the initial Fornax model and on the assumed MW potential and orbit.

free parameters (3)
  • β (fifth-force amplitude) = scanned over 0, 0.2, 0.4, 0.6, 1.0
    The dark-sector force amplitude, the target of the constraint, is scanned rather than fitted to data.
  • 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
    Chosen to reproduce Fornax's observed structure; changes with β because the relaxed equilibrium depends on the fifth force.
  • Fornax stellar Dehnen scale radius (a*) = 0.4 kpc (β=0), 0.3 kpc (others)
    Chosen to match Fornax's observed stellar half-light radius.
assumptions (5)
  • domain assumption The fifth force acts only between DM particles and is F5 = β Fgravity (with rs much larger than the system)
    Motivated by string theory scalar forces [17]; this form is assumed throughout (Eq. 1-2).
  • domain assumption The MW gravitational field is described by the adopted disk, bulge, and Dehnen halo parameters (Table I)
    The host potential sets the tidal field; uncertainty here directly propagates to the β bound.
  • domain assumption Fornax and the MW halo are initially in dynamical equilibrium and can be represented by Dehnen profiles
    Initial conditions are generated by relaxing isolated profiles (Section II.B).
  • domain assumption The orbit of Fornax is reconstructed from Gaia astrometry and the adopted MW potential via backward integration
    Table III gives the 6D phase-space; the resulting pericenter ~90 kpc controls the tidal stripping.
  • domain assumption Standard Newtonian dynamics for stellar particles; the fifth force affects only DM particles
    This is the WEP-breaking setup; if stars also felt a fifth force, the signal would differ.

how reviews work

0 comments
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

Figures reproduced from arXiv: 2509.01713 by the authors.

Figure 1
Figure 1. FIG. 1. Snapshots projected onto the YZ-Galactic plane [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Contours of the effective potential [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Time evolution of the DM (solid line) and stellar [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Profiles of stellar surface density (left) and line-of-sight velocity dispersion (right) as functions of projected distance [PITH_FULL_IMAGE:figures/full_fig_p006_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

39 extracted references · 28 canonical work pages

  1. [1]

    The stellar bulge potential is taken as a spherical Hern- quist [22] in the form of, Φbulge =−GMb r +c

    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. [2]

    (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 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. [3]

    P. J. E. Peebles and B. Ratra, Rev. Mod. Phys.75, 559 (2003)

  4. [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]

  5. [5]

    Aghanim, Y

    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]

  6. [6]

    B. Carr, F. Kühnel, and M. Sandstad, Phys. Rev. D94, 083504 (2016), arXiv:1607.06077 [astro-ph.CO]

  7. [7]

    Essig, J

    R. Essig, J. A. Jaros, W. Wester, P. H. Adrian, S. Andreas, T. Averett, O. Baker, B. Batell, M. Battaglieri,et al., Dark sectors and new, light, weakly-coupled particles (2013), arXiv:1311.0029 [hep-ph]

  8. [8]

    DESI Collaboration, arXiv e-prints , arXiv:2503.14738 (2025), arXiv:2503.14738 [astro-ph.CO]

Show all 39 references
  1. [9]

    Colladay and V

    D. Colladay and V. A. Kostelecký, Physical Review D58, 116002 (1998)

  2. [10]

    Jacobson and D

    T. Jacobson and D. Mattingly, Phys. Rev. D64, 024028 (2001), arXiv:gr-qc/0007031 [gr-qc]

  3. [11]

    J. A. Keselman, A. Nusser, and P. J. E. Peebles, Phys. Rev. D 80, 063517 (2009)

  4. [12]

    J. A. Frieman and B. A. Gradwohl, Science260, 1441 (1993)

  5. [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

  6. [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

  7. [15]

    Nusser, S

    A. Nusser, S. S. Gubser, and P. J. E. Peebles, Phys. Rev. D 71, 083505 (2005)

  8. [16]

    Kesden and M

    M. Kesden and M. Kamionkowski, Phys. Rev. D 74, 083007 (2006)

  9. [17]

    Kesden, Phys

    M. Kesden, Phys. Rev. D 80, 083530 (2009), arXiv:0903.4458 [astro-ph.CO]

  10. [18]

    Mateo, Annual Review of Astronomy and Astrophysics 36, 435 (1998)

    M. Mateo, Annual Review of Astronomy and Astrophysics 36, 435 (1998)

  11. [19]

    S. S. Gubser and P. J. E. Peebles, Phys. Rev. D70, 123510 (2004)

  12. [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)

  13. [21]

    M. G. Walker and J. Peñarrubia, The Astrophysical Jour- nal 742, 20 (2011), arXiv:1108.2404 [astro-ph.CO]

  14. [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

  15. [23]

    Miyamoto and R

    M. Miyamoto and R. Nagai, Publications of the Astro- nomical Society of Japan27, 533 (1975)

  16. [24]

    Hernquist, The Astrophysical Journal356, 359 (1990)

    L. Hernquist, The Astrophysical Journal356, 359 (1990)

  17. [25]

    Dehnen, Monthly Notices of the Royal Astronomical Society 265, 250 (1993)

    W. Dehnen, Monthly Notices of the Royal Astronomical Society 265, 250 (1993)

  18. [26]

    Ruggiero, galstep, https://github.com/ruggiero/ galstep (2014)

    R. Ruggiero, galstep, https://github.com/ruggiero/ galstep (2014)

  19. [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

  20. [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

  21. [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

  22. [30]

    A. W. McConnachie, The Astronomical Journal144, 4 (2012)

  23. [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)

  24. [32]

    Croker, Computer Physics Communications207, 478 (2016)

    K. Croker, Computer Physics Communications207, 478 (2016)

  25. [33]

    com/giulianoiorio/Gadget2-static_potential (2015)

    Giuliano, Static potential extension, https://github. com/giulianoiorio/Gadget2-static_potential (2015)

  26. [34]

    M. G. Coleman, G. S. D. Costa, J. Bland-Hawthorn, and K. C. Freeman, The Astronomical Journal 129, 1443 (2005)

  27. [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]

  28. [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]

  29. [37]

    D. Blas, M. M. Ivanov, and S. Sibiryakov, JCAP2012 (10), 057, arXiv:1209.0464 [astro-ph.CO]

  30. [38]

    Bettoni, A

    D. Bettoni, A. Nusser, D. Blas, and S. Sibiryakov, JCAP 2017 (5), 024, arXiv:1702.07726 [astro-ph.CO]

  31. [39]

    King, Astronomical Journal67, 471 (1962)

    I. King, Astronomical Journal67, 471 (1962)

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.