{"id":"7e3517cc-440e-41d4-aca0-e31e7f053e29","arxiv_id":"2509.01713","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Simulations of Fornax show a dark-matter-only fifth force stronger than about 40-60% of gravity would strip its stars, leaving β less than about 0.2 as consistent with observations.","lead":"This paper tests whether dark matter feels an extra force that normal matter does not, using simulations of the Fornax dwarf galaxy. If the force is too strong, Fornax's stars would be torn away, so the force must be weak or absent.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Initial conditions are tuned to present-day Fornax for each β, so the exclusion bound may be an artifact of that tuning rather than the fifth force.","rationale":"The reader's weakest_assumption focused on the MW potential and orbit reconstruction, which are indeed important systematics. However, I find a more fundamental issue in the simulation setup: the initial conditions are β-dependent and tuned to reproduce the present-day Fornax in isolation. This makes the exclusion bound contingent on the unverified assumption that Fornax's pre-infall structure was essentially the same as its present-day structure. If real Fornax started with a more extended DM halo or a larger stellar mass, the stripping for a given β would be weaker, potentially allowing larger β. The paper does not explore this, and the abstract/text threshold discrepancy (0.6 vs 0.4) further underscores that the bound is not pinned down. This concern does not invalidate the paper's mechanism or the qualitative trend, but it does mean the quantitative exclusion should be treated as conditional. The concrete test of using uniform initial conditions across β would settle whether the bound is a real fifth-force effect or an artifact of the fitting procedure. I therefore recommend keeping the reader's CONDITIONAL verdict, with the added condition that the initial-condition sensitivity be addressed.","tokens_in":10784,"tokens_out":13860,"duration_ms":159211,"concrete_test":"Re-run the β=0.4 and β=0.6 simulations starting from the β=0 (gravity-only) relaxed initial conditions (aDM=5.1 kpc, γ=1, same stellar profile and masses) instead of the β-tuned initial conditions. Compare the final bound stellar mass, surface-brightness, and line-of-sight velocity-dispersion profiles with the observed Fornax data. If the β=0.4 and β=0.6 models still show the claimed severe stripping, the bound is robust; if they remain broadly consistent with observations, the exclusion is an artifact of the initial-condition tuning. Additionally, produce and inspect the β=0.4 run's mass-loss curve and profiles, since the paper's stated threshold (β≳0.4) is explicitly 'not presented within the figures'.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central exclusion claim (β ≳ 0.4 in the text, β ≳ 0.6 in the abstract) depends critically on the initial conditions of the Fornax model. For each β, the initial DM halo is relaxed in isolation with the fifth force and then re-fit to reproduce Fornax's observed structural properties; Table II shows aDM decreasing from 5.1 kpc at β=0 to 2.7 kpc at β=1. This means the simulations do not hold the pre-infall Fornax model fixed while varying β; instead, each β gets a different initial model tuned to match the present-day galaxy. For high β, the more concentrated DM halo actually increases the initial stellar binding, so the dramatic stripping seen at β≥0.6 occurs despite a more tightly bound starting point. The paper never justifies why Fornax's pre-infall structure should be identical to its present-day structure (modulo the isolated relaxation), nor does it test the sensitivity of the bound to the assumed initial stellar mass or DM concentration. A forward model with a single, plausible pre-infall Fornax (e.g., the β=0 initial conditions) evolved under different β could yield a different threshold. Consequently, the exclusion bound is not robust to reasonable variations in the initial conditions, undermining the central claim.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":11147,"tokens_out":6696,"duration_ms":82597,"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":[{"comment":"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.","section":"§II.B, Table II"},{"comment":"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.","section":"§III, Fig. 4"},{"comment":"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.","section":"Abstract vs §IV"},{"comment":"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.","section":"Appendix C"}],"minor_comments":[{"comment":"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 β.","section":"§II.C"},{"comment":"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.","section":"Table IV"},{"comment":"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.","section":"Fig. 4"},{"comment":"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 β.","section":"§II.A.2"}],"recommendation":"major_revision","confidential_remarks":"The paper is in scope for the journal and the qualitative effect is interesting, but the current version overstates the constraint. The main technical issue is the re-fitting of initial conditions for each β: this must be addressed with a forward-model sensitivity study before the exclusion is credible. The abstract-versus-text threshold mismatch (β≳0.6 vs β≳0.4) and the lack of quantitative profile fits also need attention. If the authors can show that a fixed pre-infall model reproduces the same threshold, the paper would be a solid contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: This is a solid extension of the Keselman/Nusser Sgr idea to a cleaner system. Fornax is a better target than Sgr because its orbit is less eccentric and it lacks a stream, so you can use the internal structure directly. The authors do real N-body simulations with a live MW halo, a fifth force on DM, and an orbit reconstructed from Gaia. The qualitative result — that β ~ 0.5–1 rips the stars out of the DM halo — is physically plausible and consistent with earlier work.\n\nWhat's new: a Fornax-specific limit on β. They scan β, find that β≳0.6 (abstract) or β≳0.4 (Section IV) leaves a stellar remnant that doesn't match the observed surface brightness or velocity dispersion, while β≲0.2 looks fine. The appendices give an analytic tidal-stripping model that reproduces the qualitative trend. That's useful.\n\nThe soft spots are real. First, the abstract and the text disagree on the exclusion threshold: 0.6 vs 0.4. That's a red flag for a paper whose main number is a bound. Second, the comparison with observations is just overplotting simulated profiles on data; there's no likelihood or chi-square, so \"fails to match\" is a judgment call. Third, and most importantly, the initial conditions are re-fit for each β so that the isolated, relaxed Fornax matches its present-day properties. That means the β=0.6 run starts with a significantly more concentrated DM halo (aDM=3.3 kpc vs 5.1 kpc for β=0), which should make the stars harder to strip, yet they still get stripped. That argues for the fifth force being real, but it doesn't test the more natural question: what would a single pre-infall Fornax do under different β? The bound could shift if the progenitor was more or less extended. They never explore that, and the analytic check in Appendix C uses the β=0.2 DM distribution as input, so it's a consistency check, not an independent derivation. Fourth, the code link in the text isn't actually a URL.\n\nNone of these are fatal — the central mechanism holds up and the constraint is plausibly in the right ballpark — but they make the quantitative bound softer than the abstract claims. A referee should ask for a forward-model run with fixed initial conditions, a statistical comparison with the data, and a reconciliation of the abstract/text numbers.\n\nWho's this for: dark-sector phenomenologists and people working on WEP tests in dwarf galaxies. It deserves peer review, but it needs revision.","headline":"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.","tokens_in":11556,"tokens_out":3696,"would_cite":true,"duration_ms":42156,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["dark sector fifth force","weak equivalence principle breaking","dwarf spheroidal galaxy","tidal stripping","Fornax","N-body simulations","surface brightness","velocity dispersion"],"falsifier":"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.","tokens_in":10744,"feed_emoji":"🔭","tokens_out":9679,"duration_ms":95059,"temperature":0.7,"pith_summary":"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.","feed_headline":"A dark-only force above 40% of gravity would have destroyed Fornax","feed_subtitle":"Simulations show a dark-only force near half of gravity would strip the dwarf's stars; below 20 percent it survives.","key_machinery":"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","core_discovery":"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","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"supplies the self-bound stellar core picture for a satellite decoupled from its dark halo, and parameters used in the adopted Milky Way model.","marker":"[9]"},{"why":"introduced fifth-force-induced stream asymmetry in Sagittarius as a constraint, the approach this paper contrasts with Fornax.","marker":"[14]"},{"why":"provides the phase-space measurements of Fornax used to reconstruct its Galactic orbit.","marker":"[18]"},{"why":"gives the stellar and dark-matter mass-model constraints used to construct the Fornax initial conditions.","marker":"[19]"},{"why":"supplies independent Fornax mass-model constraints used in the initial setup and comparison.","marker":"[20]"},{"why":"provides Milky Way potential parameters adopted for the host galaxy model.","marker":"[25]"},{"why":"also contributes to the adopted Milky Way model parameters.","marker":"[26]"},{"why":"provides the observed surface-brightness profile against which simulated stellar profiles are compared.","marker":"[32]"},{"why":"provides the observed line-of-sight velocity-dispersion profile used to assess consistency.","marker":"[33]"}],"fun_headline_variants":["Dark force above 40% of gravity would strip Fornax's stars","Fornax rules out dark-sector fifth force stronger than 0.4","Simulations: dark force must be under 0.4 to keep Fornax intact","Fornax constrains dark-only force to below 40% of Newtonian"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Dark force above 40% of gravity would strip Fornax's stars","Fornax rules out dark-sector fifth force stronger than 0.4","Simulations: dark force must be under 0.4 to keep Fornax intact","Fornax constrains dark-only force to below 40% of Newtonian"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000286,"raw_usage":{"total_tokens":1509,"prompt_tokens":722,"completion_tokens":787,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":466,"completion_tokens_details":{"reasoning_tokens":700}},"tokens_in":466,"tokens_out":787,"duration_ms":8232,"temperature":1.0,"reasoning_tokens":700,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T12:15:01.155957+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Colladay and V","cited_arxiv_id":null,"evidence_quote":"supplies the self-bound stellar core picture for a satellite decoupled from its dark halo, and parameters used in the adopted Milky Way model."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"introduced fifth-force-induced stream asymmetry in Sagittarius as a constraint, the approach this paper contrasts with Fornax."},{"cited_title":"Mateo, Annual Review of Astronomy and Astrophysics 36, 435 (1998)","cited_arxiv_id":null,"evidence_quote":"provides the phase-space measurements of Fornax used to reconstruct its Galactic orbit."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"gives the stellar and dark-matter mass-model constraints used to construct the Fornax initial conditions."},{"cited_title":"K., Battaglia, G., Pawlowski, M","cited_arxiv_id":null,"evidence_quote":"supplies independent Fornax mass-model constraints used in the initial setup and comparison."},{"cited_title":"Dehnen, Monthly Notices of the Royal Astronomical Society 265, 250 (1993)","cited_arxiv_id":null,"evidence_quote":"provides Milky Way potential parameters adopted for the host galaxy model."},{"cited_title":"Ruggiero, galstep, https://github.com/ruggiero/ galstep (2014)","cited_arxiv_id":null,"evidence_quote":"also contributes to the adopted Milky Way model parameters."},{"cited_title":"Croker, Computer Physics Communications207, 478 (2016)","cited_arxiv_id":null,"evidence_quote":"provides the observed surface-brightness profile against which simulated stellar profiles are compared."},{"cited_title":"com/giulianoiorio/Gadget2-static_potential (2015)","cited_arxiv_id":null,"evidence_quote":"provides the observed line-of-sight velocity-dispersion profile used to assess consistency."}],"review_version":1}