{"id":"468227b3-c76c-4c47-9ecf-6f7b85ea5979","arxiv_id":"2510.17704","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A single relativistic scalar field with a Bekenstein-type kinetic term and a mass-generating term reproduces deep-MOND acceleration in the weak-field limit and predicts stronger-than-MOND light bending.","lead":"A single scalar-field theory of gravity is shown to reproduce MOND's low-acceleration dynamics while carrying real mass-energy, so it acts like dark matter in galaxy clusters. It predicts light bending different from ordinary MOND, but leaves cosmology untouched.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The transition function (7) is implemented backwards: as written, h=0 in the deep-MOND regime and h=1 in the Newton regime, so the Bekenstein/mass terms are off where eq. (40) is supposed to apply.","rationale":"The reader's weakest assumption (dropping p3 in (26)→(27)) is secondary. In the baryonic vacuum the p3 terms are O(a1 r) relative to the retained ∇(|∇σ|∇σ) term, so over 1–500 kpc their effect on the radial profile is at the 10^−4 level; the dimensional comparison to Λ is sloppy but not the main failure. The more load-bearing problem is that the transition function h, as explicitly defined in (7), is active in the wrong regime: with the published increasing g, h is 0 for small gradients and 1 for large gradients, exactly opposite to the stated assumption that the MOND terms are switched on only below a threshold. This directly invalidates eq. (40) as a consequence of the Lagrangian unless the definition is corrected. It is an internal inconsistency, checkable by direct evaluation, rather than a disagreement with existing consensus. I keep the CONDITIONAL verdict: the issue is local and plausibly fixable, but the manuscript as written does not support the central claim.","tokens_in":52052,"tokens_out":26259,"duration_ms":218140,"concrete_test":"Evaluate (7) at x=0.05 (deep-MOND side, a_N/a0≈0.1): (0.05−0.1)/9.9 < 0, so h=0, while the derivation of (40) requires h=1; at x=20, h=1 although this is the Newtonian side. Plot h over 0.01 ≤ x ≤ 100 to confirm the inversion. If confirmed, replace h by its complement (or redefine x) and recompute the radial-acceleration relation and rotation curves; the central claim cannot stand with the published definition.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (7) defines h(x;α,β)=g((x−α)/(β−α)) with g=f(x)/(f(x)+f(1−x)), f=e^{−1/x} for x>0 and f=0 for x≤0. This g is strictly increasing from 0 to 1, so h=0 for x≤α=0.1 and h=1 for x≥β=10. With x=|∇σ|/a1 (p.11) or y=a_N/a0 (p.25), the deep-MOND regime is x,y≪1; (7) gives h=0 there. Then the Lagrangian (8) reduces to L_V alone, and equations (26)–(27) — derived under 'h≡1 in the Milgrom regime' — are not the governing equations in the regime where MOND is claimed. The text later says 'agrees exactly with our (58) for h≡1 (i.e. for a_N≤α a0)' and treats h=1 on the deep-MOND side, but that is not the function defined in (7). This is a direct internal inconsistency in the gating of the central mechanism, not a matter of coefficient size.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a relativistic scalar-field extension of Einstein gravity, formulated in integrable Weyl geometry, with a non-minimally coupled scalar field. The Lagrangian contains a Bekenstein-type cubic kinetic term and a second-order mass-generating term, both supposed to be active only when the scalar-field gradient is spacelike and below a MOND-type threshold. The central claim is that, in the weak-field Einstein-gauge limit, the scalar field obeys the deep-MOND equation (40), so that free-fall trajectories are MONDian while light deflection is governed by the modified metric (55), with the scalar potential contributing twice as strongly to lensing as baryonic matter of equal potential. The paper also presents a numerical central-symmetric solution, a Kuzmin-disk analysis, a radial-acceleration relation, and a Coma-cluster estimate.","tokens_in":52380,"tokens_out":6983,"duration_ms":61345,"significance":"If the derivation were sound, the paper would be a noteworthy contribution: it offers a single-scalar-field relativistic MOND framework, with explicit variational calculations (Appendix 6.2), a transparent weak-field reduction, a concrete lensing prediction that differs from the usual MOND expectation, and an attempt to address cluster mass discrepancies. The honest discussion of the model's limitations and the detailed numerical checks are also strengths. However, the central result is compromised by an internal inconsistency in the transition function and by an unjustified smallness estimate for the cubic correction; as written, the deep-MOND equation is not actually the governing equation in the regime where the model claims to reproduce MOND.","major_comments":[{"comment":"The transition function is implemented backwards. With h(x;ᾱ,β̄)=g((x−ᾱ)/(β̄−ᾱ)) and g increasing from 0 to 1, h=0 for x≤ᾱ=0.1 and h=1 for x≥β̄=10. Since x=|∇σ|/a1, the deep-MOND regime is x≪1, so h=0 there. Thus the Lagrangian (8) reduces to L_V alone precisely where the Bekenstein and mass terms are supposed to act. Equations (26)–(27) are derived under 'h≡1 in the Milgrom regime', and §3.3 later states that (58) agrees with observation for 'h≡1 (i.e. for a_N≤ᾱ a0)', which contradicts the definition (7). This is not a coefficient-size issue: the gating of the central mechanism is inverted and must be corrected (e.g. by using 1−h or a decreasing transition function) before the derivation of (40) and (58) can be accepted.","section":"Eq. (7), (8), §3.3"},{"comment":"The discard of p3(x)=(x−2β^{−1}a1)x^2 is not justified. The stated bound |p3|≤c^{−3}a0^3≪Λ compares a quantity of dimension L^{−3} with Λ of dimension L^{−2} in the paper's own dimensional conventions; the comparison is dimensionally inconsistent. Moreover, even setting dimensions aside, the model's own transition function keeps the Bekenstein and mass terms active up to x=10, where p3∼900 a1^3, far from negligible. Since (27) is the basis for the deep-MOND equation (40) and for the radial-acceleration relation (58), the derivation needs either a correct estimate over the full active interval or an explicit restriction of the claimed regime.","section":"Eqs. (26)–(27), p. 15"},{"comment":"The paper's headline 'single scalar field' is qualified by several additional, independently chosen structures: a non-dynamical timelike unit vector field A^μ in the mass term, the smooth threshold function h (called 'at best metaphoric' on p. 11–12), and constants α=−4, β=2, c1=12√(a1M), a0, λ, and the transition boundaries. The deep-MOND equation (40) is therefore inherited by construction rather than predicted. This is acceptable for a phenomenological model, but the abstract and introduction should state the framework as a parametrized scalar-tensor model, not as a 'single scalar field' derivation of MOND.","section":"§2.1, Eqs. (5), (10), (11)"}],"minor_comments":[{"comment":"The title contains a typo: 'relativistc' should be 'relativistic'.","section":"Title"},{"comment":"The figure-2 caption and the surrounding text include a large, apparently accidental quotation from Hossenfelder and Mistele's paper, with headers such as '4 Comparison with Observation' and discussion of Verlinde-matching. This passage is not integrated into the present argument and should be removed or clearly set off as a quotation with explicit attribution.","section":"Fig. 2 and §3.3"},{"comment":"The upper smoothing function h̃(x;α̂,β̂)=1−g((x−α̂)/(β̂−α̂)) used for the halo cut-off has the opposite monotonicity from h in Eq. (7). This is presumably the intended direction for a suppression mechanism, but the two uses of 'h' with opposite behavior are confusing and should be reconciled notationally.","section":"§4.1"},{"comment":"The transition variable is x=|∇σ|/a1 in Eq. (7) but y=a_N/a0 in Eq. (57), with the same symbol h used in both places. Since the threshold conditions are different (gradient versus acceleration), the paper should explicitly state the mapping and any approximation used to convert one into the other.","section":"§3.3"}],"recommendation":"major_revision","confidential_remarks":"The central idea is potentially salvageable: if the transition function is corrected so that the Bekenstein/mass terms are active in the deep-MOND regime, and if the p3 smallness estimate is repaired or the active regime restricted, the paper could become a useful contribution. In its current form, however, the main result (40) does not follow from the stated Lagrangian. The inserted quotation in §3.3 also needs editorial cleanup. I recommend major revision rather than rejection, because the gating inversion appears to be a fixable technical error rather than a fundamentally flawed approach."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: there is a sign error in the gating function. Equation (7) defines h as increasing from 0 to 1, so h=0 for small |∇σ|/a1 and h=1 for large values. The authors intend the Bekenstein and mass terms to be active below the MOND threshold and inactive above. The Lagrangian in (8) multiplies those terms by h, so with (7) they vanish exactly in the deep-MOND regime where the paper derives equation (40). The text even says 'agrees exactly with our (58) for h≡1 (i.e., for a_N ≤ α a0)' — but for a_N ≤ α a0, h=0, not 1. This is not a coefficient-size issue; it is an internal contradiction in the central mechanism. The derivation of the Milgrom equation from the Lagrangian is therefore not valid as written.\n\nThe paper does a lot right. The Weyl-geometric framework is used carefully, the variational derivatives in the appendix are detailed and consistent, and the second-order mass term L_2φ is genuinely new: it gives the factor-2/-6 lensing modifications and the scalar field a dark-matter-like energy content. The survey of other relativistic MOND attempts is useful and honest. The author also concedes the main limitations: no fundamental theory, no cosmology, and the cluster estimate leaves a discrepancy up to 8.7 at small radii.\n\nThe soft spots beyond the sign error: the p3 drop at (26)-(27) is justified only for |∇σ|≤a1 with a dimensionally inconsistent estimate, while the transition interval extends to 10 a1. The screening mechanism is borrowed from superfluid dark matter and explicitly called 'at best metaphoric.' The Coma numbers are tuned to match baryonic masses and have no error bars.\n\nBottom line: the variational machinery may be salvageable, and the lensing predictions are worth testing once the gating is corrected. But as it stands, the central claim that this Lagrangian produces MOND dynamics does not follow. I'd send it for peer review — a good referee can pinpoint the gating fix — but the revised version must re-derive the weak-field limit with a corrected h.","headline":"The transition function in this otherwise careful paper is backwards, so the MOND terms are switched off precisely in the deep-MOND regime; the central derivation needs a fix before it can be trusted.","tokens_in":52941,"tokens_out":3114,"would_cite":false,"duration_ms":26568,"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 single relativistic scalar field, with two threshold-gated kinetic terms in integrable Weyl geometry, is claimed to reproduce deep-MOND dynamics in the weak-field limit and to modify gravitational light deflection by exactly twice the ord","keywords":["MOND","scalar field","Weyl geometry","deep MOND equation","gravitational lensing","radial acceleration relation","dark matter","modified gravity"],"falsifier":"For a spherically symmetric galaxy with a known baryonic mass distribution and a measured rotation curve in the deep-MOND regime, the model predicts that the scalar-field contribution to the gravitational lensing deflection angle is exactly twice the value expected from the phantom-matter distribution under the usual factor-2 rule. A measurement of the Einstein radius in a strong-lensing galaxy, at 10% precision, that agrees with the standard factor-2 rule would falsify this prediction. Alternatively, numerically integrating the full scalar field equation including p3 for a realistic mass dist","tokens_in":51747,"feed_emoji":"🌌","tokens_out":6342,"duration_ms":49807,"temperature":0.7,"pith_summary":"This paper argues that one scalar field can do double duty: its Bekenstein-type cubic kinetic term generates deep-MOND dynamics for test-particle motion, while a second-order mass-generating term gives the field a non-negligible energy-momentum tensor that acts like a form of dark matter and changes how light is deflected. In the weak-field (Newton-Milgrom) approximation, the total gravitational potential is the sum of the ordinary baryonic Newtonian potential and a scalar-field potential that satisfies the deep MOND equation ∇·(|∇Φ|∇Φ)=a0(4πG)ρ_m. Because the scalar-field energy tensor is traceless and carries pressure, the scalar potential enters the metric with a relative factor: in the centrally symmetric case, lensing from the scalar field is twice as strong as from ordinary matter with the same potential. The model thus offers a concrete bridge between the dark-matter and modified-gravity readings of the missing-mass problem, with distinctive, testable lensing predictions.","feed_headline":"One scalar field yields MOND plus dark-matter lensing","feed_subtitle":"The model's scalar potential follows the deep MOND equation and bends light twice as strongly as ordinary matter.","key_machinery":"The argument is carried by the scalar field φ with Weyl weight -1, written in the Einstein gauge as σ = -ln(φ/φ0). Its Lagrangian contains a Bekenstein-type cubic kinetic term L_φ3 ∝ φ^{-2}|Dφ|^3 and a second-order term L_2φ ∝ φ DλDλφ AλAλ, both gated by a transition function h(|∇σ|/a1) that suppresses them when the gradient is timelike or above the MOND threshold. The reduced scalar field equation, obtained by combining the trace of the Einstein equation with the ϕ-variation, simplifies to the relativistic Milgrom equation in the Milgrom regime. On the gravity side, the key identity is that the scalar field's contribution to the source in the Poisson equation includes a pressure term, givin","core_discovery":"In the Einstein gauge, the model's scalar field σ is related to the Newtonian potential by Φ^(φ)=c^2 σ. The paper derives the relativistic Milgrom equation ∇λ(|∇σ|∂^λσ) = -a1 β^{-1}(8πκ) tr T from the reduced scalar field equation, after dropping a cubic polynomial p3(x)=(x-2β^{-1}a1)x^2 as negligible in the Milgrom regime. In the flat-space weak-field limit this becomes the deep MOND equation (40). The scalar field's effective energy tensor, extracted from the Einstein equation, is traceless and dominated by second-order derivatives; its Newtonian mass equivalent is twice its energy density because of pressure. Combining the baryonic and scalar potentials, the perturbed metric is g = -(1+2Φ","pith_inferences":["If the neglected cubic polynomial p3 is retained in the scalar field equation, the deep MOND equation receives corrections near the upper transition boundary; these could be tested with high-precision rotation-curve data in the 1–10 a0 range.","The anisotropic lensing prediction for disk galaxies is a distinctive signature that could distinguish this model from other relativistic MOND approaches.","The paper's suggestion that the scalar field's energy content accounts for cluster missing mass implies a natural extension to cosmological structure formation, though the model is currently silent on early-universe dark matter."],"forward_implications":["Galactic rotation curves can be explained without particle dark matter: the scalar field produces MONDian accelerations and also contributes its own energy-momentum as an effective dark component.","Gravitational lensing in deep-MOND galaxies should be stronger than standard MOND's phantom-matter estimate by a factor of two in spherically symmetric cases, and anisotropic close to disk planes.","The model preserves exact Newton/Einstein behavior above the threshold (a_N > 10 a0), so solar-system tests are unaffected.","In galaxy clusters, the scalar-field halos of galaxies and hot gas add a Newtonian mass equivalent that lowers the missing mass ratio; for Coma the estimate reduces the discrepancy but may not fully close it."],"fun_headline_variants":["Scalar field bridges dark matter and MOND","One field explains MOND and dark-matter lensing","Relativistic scalar field mimics dark matter and MOND","Unified model: MOND dynamics plus lensing","Scalar field yields MOND and bends light like dark matter"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the cubic polynomial p3(x)=(x-2β^{-1}a1)x^2 is negligible throughout the Milgrom regime, so the relativistic Milgrom equation and the deep MOND equation follow from the full scalar field equation; if p3 is not small when the field gradient approaches the transition threshold, the model's MOND predictions acquire corrections.","fun_headline_variants_meta":{"raw":{"variants":["Scalar field bridges dark matter and MOND","One field explains MOND and dark-matter lensing","Relativistic scalar field mimics dark matter and MOND","Unified model: MOND dynamics plus lensing","Scalar field yields MOND and bends light like dark matter"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00013,"raw_usage":{"total_tokens":948,"prompt_tokens":713,"completion_tokens":235,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":457,"completion_tokens_details":{"reasoning_tokens":166}},"tokens_in":457,"tokens_out":235,"duration_ms":2557,"temperature":1.0,"reasoning_tokens":166,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T08:58:58.779551+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For a spherically symmetric galaxy with a known baryonic mass distribution and a measured rotation curve in the deep-MOND regime, the model predicts that the scalar-field contribution to the gravitational lensing deflection angle is exactly twice the value expected from the phantom-matter distribution under the usual factor-2 rule. A measurement of the Einstein radius in a strong-lensing galaxy, at 10% precision, that agrees with the standard factor-2 rule would falsify this prediction. Alternatively, numerically integrating the full scalar field equation including p3 for a realistic mass dist","supporting_citations":[],"review_version":1}