{"id":"2e428646-2625-4f0a-a1d7-6004f07d94a0","arxiv_id":"2512.13798","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Spatially varying scalar fields rescale neutrino masses and arrival times, and neutrino data bound these rescaling parameters for Symmetron and Chameleon models.","lead":"This paper derives new constraints on scalar–tensor gravity theories by asking how a spatially varying scalar field would change the mass of neutrinos traveling through Earth and from supernovae. It matters because neutrinos can cross dense matter and cosmic distances, so they could reveal fifth-force effects other experiments cannot see.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Local-density approximation inconsistent with stated thin-shell regime; atmospheric-neutrino exclusions in Figs. 1–2 may be overestimated without a full Klein–Gordon profile.","rationale":"The reader identifies the local-density approximation as the weakest assumption; I agree and sharpen the concern. The paper's own validity condition (thin shell) points to the regime where the local-density approximation is least reliable: in a screened object, the scalar field is nearly constant in the interior, not following the local density. The proposed test—solving the Klein–Gordon equation for a representative point—would settle whether the local-density approximation is the controlling systematic. If the full profile yields the same bounds, the central claim stands; if not, the claim must be restricted or the analysis redone. This is a correctness risk, not a matter of consensus. Secondary issues include the missing statistical details of the IceCube DeepCore analysis and the inconsistency between 'eight years' in the main text and '9.3 years' in the supplemental caption, but these are reproducibility concerns; the physics approximation is more fundamental. The verdict remains CONDITIONAL: accept only after the profile dependence is checked and the IceCube analysis is made reproducible.","tokens_in":16735,"tokens_out":17258,"duration_ms":143524,"concrete_test":"Pick a representative point inside the claimed Symmetron exclusion (e.g., μ=1 meV, M_s=10^4 GeV, λ chosen so the effective coupling α_s,0 is O(0.1)). Solve the static, spherically symmetric Klein–Gordon equation for φ(r) in the Earth using the PREM density profile and the Symmetron potential, with the boundary condition φ→φ_vac at large radius. Compute A(φ(r)) and then the oscillation probability using Eq. (7) with this full profile. Compare the resulting L/E event distribution (and the derived 95% C.L. contour) with the local-density result of Fig. 5. If the full-profile distribution differs meaningfully from the local-density prediction, the exclusion contour in Figs. 1–2 is not robust. A simpler consistency check: verify whether |A(r)-A(ρ(r))|/A exceeds 10% at any point along a typical atmospheric neutrino chord; if yes, the stated thin-shell validity condition is not sufficient for t","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim is that atmospheric neutrino oscillations (IceCube DeepCore) and SN1987A time delays exclude previously unexplored scalar–tensor parameter space (Symmetron M_s∼10^4 GeV). The oscillation calculation uses Eq. (7)/(58), where the phase is ∫ A^2 H dr, with A(φ) replaced by A(ρ(r)) via Eq. (6). The authors acknowledge this and state the constraints apply when the scalar forms a thin shell around Earth (Compton wavelength between atomic and planetary scales). The problem is that in the thin-shell regime the scalar field is screened inside Earth: φ is pinned to its interior minimum and does not track the local PREM density. For the Symmetron, for ρ>ρ_crit the field sits at φ=0 so A=1; for ρ<ρ_crit and ρ/ρ_crit≪1 (which holds for M_s=10^4 GeV, μ=meV), A≈1+constant, again essentially independent of ρ. Thus the computed position-dependent distortion of the oscillation length (Fig. 4) does not correspond to the thin-shell field profile. The local-density approximation describes an unscreened, density-following field; the thin-shell condition is precisely where that approximation fails. Therefore the excluded regions in Fig. 1–2 are not established for the quoted validity regime unless a full Klein–Gordon profile calculation reproduces the assumed A(r). The SN time-delay bound is less affected because it depends on A_∞, though it also uses a single average density rather than a profile integral; this is subdominant.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives new constraints on scalar–tensor (ST) theories from neutrino physics. It argues that a spatial background scalar field φ, through a conformal factor A(φ), makes neutrino masses and matter potentials position-dependent. The central observables are the flavour-oscillation phase for atmospheric neutrinos crossing the Earth (Eqs. (7)/(58)) and the supernova time delay (Eq. (11)). Using an effective power-law parametrisation A(ρ) (Eq. (6)) and IceCube DeepCore data, the authors exclude part of the Symmetron parameter space around M_s ~ 10^4 GeV (Fig. 2), while SN1987A time delays give subdominant bounds. The Supplemental Material derives the conserved quantity K_ST and the phase in a frame-invariant way.","tokens_in":17002,"tokens_out":10802,"duration_ms":90700,"significance":"If correct, the bounds would provide a genuinely new probe of scalar–tensor theories: neutrinos traverse dense regions (Earth) and long baselines (supernova), complementing laboratory and solar-system tests. The derivation of the oscillation phase from K_ST and the frame-invariance argument in the Supplemental Material are physically coherent and a clear strength. The paper also performs a real data analysis (IceCube DeepCore) rather than a purely illustrative estimate, and the SN1987A time-delay bound is straightforwardly falsifiable. However, the central exclusion regions depend on an approximation that the paper itself flags, and the Symmetron mapping contains an internal inconsistency. These issues must be resolved before the bounds can be considered established.","major_comments":[{"comment":"The paper states that A(φ) was assumed to follow the local density exactly, and that the constraints apply when the scalar forms a thin shell around Earth and its core (10^-10 < m_φ/eV < 10^3). This is internally inconsistent. In the thin-shell regime the field is pinned to its density-dependent minimum in the interior and varies only in a thin surface layer, so A(r) is piecewise approximately constant along an atmospheric neutrino trajectory, not the smooth density-following function used in Eq. (6) and the phase integral Eq. (58). For the Symmetron, A_s=1 for ρ>ρ_crit and A_s≈1+const for ρ<ρ_crit (Eq. (14)), so the density-dependent distortion shown in Fig. 4 and the resulting IceCube exclusions in Figs. 1–2 do not follow from the stated validity regime. A full Klein–Gordon profile calculation is needed to support the quoted bounds.","section":"SYMMETRON AND CHAMELEON, after Eq. (15)"},{"comment":"The Symmetron is matched by setting α_{s,0}=μ^2/(2λM_s^2) with index n=0. However, the effective parametrization Eq. (6) is A≈1+α_n[(ρ/ρ0)^n − (ρ0/ρ0)^n]; for n=0 the bracket vanishes identically, so the density-dependent term is zero. The Symmetron conformal factor in the unscreened regime (Eq. (14)) is linear in ρ, A_s≈1+μ^2(1−ρ/ρ_crit)/(2λM_s^2), which should be matched with n=1 after normalizing at Earth’s surface. As written, the mapping cannot generate the claimed Symmetron exclusion region in Fig. 2(a). Please clarify or correct the index.","section":"Eqs. (6) and (15)"},{"comment":"The derivation of the effective density accessible via seismic measurements relies on the estimate B∝m_e^4. The authors themselves note that no general analytical formula for B exists. Since the matter potential V_CC and the oscillation phase depend on the inferred density, a different scaling of B with m_e would alter the constraints. The sensitivity of the IceCube bound to this assumption should be quantified, or the derivation should be made more rigorous.","section":"Section 'ANALYSIS', Eq. (9)"}],"minor_comments":[{"comment":"The main text states 'We analysed eight years of IceCube DeepCore data [79]', while the Supplemental Material and Fig. 5 refer to 9.3 years of data and Ref. [93]. The dataset used for the main bound should be specified consistently.","section":"Main text vs Supplemental Material"},{"comment":"The notation ρ0(x)/g/cm^3 is confusing. Clarify that densities are in units of g/cm^3 and that ρ0 is the Earth-surface density, a constant.","section":"Eq. (6)"},{"comment":"Several equations contain garbled symbols (e.g., '⌟roo⟪...') in the typeset version; these should be cleaned up.","section":"Eqs. (28) and (39)"},{"comment":"The word 'unicum' is non-standard; consider 'unique' or 'singular'.","section":"Introduction"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about the local-density approximation is valid and lands: the paper explicitly restricts its validity to the thin-shell regime, where the assumed density-following A(ρ) is precisely the wrong profile. Combined with the n=0 Symmetron matching issue, the central atmospheric-neutrino exclusion is not currently supported. The time-delay bound is less affected. I believe these issues are fixable within the scope of a revision, hence major_revision rather than reject."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick read on arXiv:2512.13798. The paper makes a real, if modest, contribution: it writes down the scalar-tensor corrections to neutrino oscillation phases and propagation times in a frame-consistent way, then translates existing IceCube DeepCore and SN1987A data into exclusions on the Symmetron and Chameleon parameters. The headline new result is an exclusion around M_s ~ 10^4 GeV for the Symmetron from atmospheric neutrinos; I don't know an existing bound there.\n\nWhat's well done: the derivation via the conserved quantity K_ST = g_00 p_0 is clean, and the frame-invariance argument for the observables is a useful correction to sloppier treatments in the literature. The density-mapping argument that seismological densities correspond to the Jordan frame (B ∝ m_e^4) is heuristic, but they flag it as an estimate. They are also transparent about the main limitation: A(φ) is assumed to follow the local density rather than solving the Klein-Gordon profile, and they state the thin-shell condition.\n\nNow the stress-test: it claims that in the thin-shell regime the field is pinned and does not track PREM density, so the position-dependent A(r) used for the atmospheric analysis is invalid. That objection does not hold up for the parameters they target. For M_s ~ 10^4 GeV and μ ~ meV, Eq.(16) gives ρ_crit ~ 20 g/cm^3, so Earth's densities are a fair fraction of ρ_crit and A varies appreciably with ρ. The Compton wavelength for μ ~ meV is sub-mm, well below Earth's density gradients, so the local-density approximation is justified in the stated regime. The stress-test's claim that ρ/ρ_crit ≪ 1 is numerically wrong in that region. The paper's own caveat is still real: the exclusion should be restricted to the validity window 10^-10 eV < m_φ < 10^3 eV, and the authors should check that the contours in Figs. 1-2 lie inside it.\n\nWhere the paper is soft: the IceCube analysis is deferred to the Supplemental Material, and the numbers are inconsistent—main text says eight years of DeepCore data, the supplement says 9.3 years in the figure caption. No likelihood or systematic description is given; the reader cannot reproduce the exclusion. The SN1987A bound rests on one average interstellar density, which they note is conservative, but the sub-dominance claim is fine. The mapping of the effective parameters to the Symmetron (α_s,0) is written in a way that might confuse; a referee should ask for the explicit identification.\n\nBottom line: this is a solid phenomenology paper within its stated assumptions. It deserves a serious referee, but acceptance should be conditional on the IceCube analysis being made public and the validity-region check. The stress-test's main concern is a misfire, but it points to the need to clarify the regime. Worth a reading-group session if you work on screening mechanisms or neutrino phenomenology; I'd cite it as the current neutrino-derived Symmetron constraint.","headline":"Solid, moderately novel neutrino-frame derivation that yields a genuine Symmetron exclusion; the stress-test's main objection is a misfire, but the deferred IceCube analysis and the paper's own thin-shell caveat are the real soft points.","tokens_in":17607,"tokens_out":9203,"would_cite":true,"duration_ms":79088,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.50.Kd","14.60.Pq"],"model":"deepseek-v4-flash","headline":"Neutrino data exclude a previously untested region of the Symmetron scalar-tensor model.","keywords":["scalar-tensor gravity","neutrino oscillations","mass-varying neutrinos","symmetron","chameleon","supernova time delay","atmospheric neutrinos","screening mechanisms"],"falsifier":"Solve the full Klein–Gordon equation for the Symmetron profile inside the Earth at the excluded parameters (M_s~10^4 GeV, μ~meV). If the field does not track the local density and the thin-shell condition fails, the A(r)^2 rescaling disappears and the excluded region evaporates; if the profile confirms the tracking assumption, the exclusion stands.","tokens_in":16521,"feed_emoji":"⚛️","tokens_out":7296,"duration_ms":61167,"temperature":0.7,"pith_summary":"Scalar–tensor theories of gravity predict that a background scalar field shifts Standard Model masses with the local density. Because neutrinos can cross dense matter and travel galactic distances, they are uniquely placed to feel this effect. This paper derives how a density-dependent rescaling A(φ) changes neutrino oscillation phases and time-of-flight, then uses atmospheric-neutrino oscillation data and the 1987 supernova neutrino burst to set new limits. It finds that Earth-crossing oscillations exclude a region of the Symmetron model around M_s ≈ 10^4 GeV, while current supernova timing is weaker but a future nearby supernova could reach inside the theory's valid region. If correct, neutrino observatories become a new, complementary probe of screened modified gravity.","feed_headline":"Atmospheric neutrinos exclude a new symmetron gravity region","feed_subtitle":"Mass-changing scalar fields leave a measurable imprint on Earth-crossing oscillations and supernova arrival times.","key_machinery":"The central object is the coupling/conformal factor A(φ) that rescales all Jordan-frame masses, expressed in an effective parametrization A(φ)=1+α_n[(ρ/ρ0)^n − 1]. The load-bearing identity is that both the neutrino oscillation phase ∫A(r)^2 H(r) dr and the supernova time delay ∫A(r)^2 dr are integrals of the square of this factor along the trajectory, together with the frame-invariant conserved quantity K_ST = g00 p0. This uniform A(r)^2 scaling is what lets two different observables — flavour oscillations in Earth and arrival-time delays from a supernova — constrain the same parameter combinations, and why adiabatic propagation in the Sun or a supernova shows no effect while non-adiabatic","core_discovery":"The paper's central claim is that the same background scalar field that mediates screened fifth forces also rescales particle masses by A(φ), and that neutrinos — because they pass through dense media and travel astrophysical distances — can observe the resulting spacetime variations. In the Jordan frame, the neutrino oscillation phase acquires an overall factor A(r)^2, rescaling the effective oscillation length but leaving mixing angles intact under adiabatic propagation; the supernova time delay is likewise multiplied by the path integral of A(r)^2. The authors show, from atmospheric-neutrino oscillation data with the effective parametrization A(φ)=1+α_n[(ρ/ρ0)^n − 1], that there is no pre","pith_inferences":["The same mass-rescaling mechanism would affect cosmic-ray propagation or dark-matter annihilation signals in dense regions; neutrino telescopes might therefore be used to map density profiles of dark-matter halos.","If a future measurement of the atmospheric L/E spectrum finds the first oscillation minimum shifted exactly as predicted by the A(r)^2 factor, it would confirm the density-tracking assumption; conversely, a null result would weaken the claimed exclusion.","The argument that the seismically inferred density is the Jordan-frame one hinges on the bulk modulus scaling as the fourth power of the electron mass; a different scaling would alter the calibration of ρ0 and shift the excluded region.","The constraints depend on the dark-matter halo model; if the Milky Way's halo is clumpy rather than smooth, the time-delay bound could strengthen considerably."],"forward_implications":["Atmospheric neutrino oscillation experiments can place leading constraints on quadratic scalar-field couplings (Symmetron) that laboratory fifth-force searches cannot easily reach.","A future Galactic supernova could lower the bound on the effective mass rescaling A_∞ to ~0.06 eV, reaching inside the theory's EFT-valid region.","The same A(r)^2 scaling applies to any dense medium; longer-baseline or high-statistics neutrino measurements through the Earth's core should sharpen the excluded region.","The parametrization covers any screening model with a density-power-law A(φ); constraints can be mapped to steeper potentials such as φ^{2n}.","Neutrino time delays become a probe of the integrated dark-matter density along the line of sight, not just of the theory's coupling."],"fun_headline_variants":["Neutrino oscillations constrain scalar-tensor gravity","Scalar fields tested by Earth-crossing neutrinos","Symmetron region excluded by atmospheric neutrinos","Neutrino time delays probe mass-changing scalars"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The constraints assume the scalar field instantly follows the local matter density (A(φ)=A(ρ(r))), so they hold only when the scalar forms a thin shell around Earth and its core; a full treatment would solve the Klein–Gordon equation for the field profile.","fun_headline_variants_meta":{"raw":{"variants":["Neutrino oscillations constrain scalar-tensor gravity","Scalar fields tested by Earth-crossing neutrinos","Symmetron region excluded by atmospheric neutrinos","Neutrino time delays probe mass-changing scalars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000191,"raw_usage":{"total_tokens":1131,"prompt_tokens":649,"completion_tokens":482,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":393,"completion_tokens_details":{"reasoning_tokens":421}},"tokens_in":393,"tokens_out":482,"duration_ms":5097,"temperature":1.0,"reasoning_tokens":421,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T16:20:11.614690+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the full Klein–Gordon equation for the Symmetron profile inside the Earth at the excluded parameters (M_s~10^4 GeV, μ~meV). If the field does not track the local density and the thin-shell condition fails, the A(r)^2 rescaling disappears and the excluded region evaporates; if the profile confirms the tracking assumption, the exclusion stands.","supporting_citations":[],"review_version":1}