{"id":"116b45d7-b06f-49ad-abc0-8a5f276c3cf8","arxiv_id":"1909.02552","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"Helioseismic frequencies of 19 selected solar modes could constrain the DHOST fifth-force coupling Y to about 10^-3, far tighter than previous stellar bounds.","lead":"This paper proposes using the Sun's acoustic oscillations, measured by helioseismology, to test for a fifth force predicted by modified gravity theories. It estimates that current solar data could constrain the fifth-force strength about two orders of magnitude more tightly than white dwarf observations.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central Y constraint rests on assuming f'_theory ≈ f'_pol with a single scalar ξ in [-0.22, 0.22]; no modified-gravity solar model validates this, so Eq. (7) is only as strong as the polytropic proxy.","rationale":"The reader's weakest-assumption analysis identifies exactly the same load-bearing point: the unknown true frequency response f'_theory is replaced by the polytropic slope f'_pol times a nuisance ξ marginalized over a range estimated from the same polytropic and WKB approximations. My reading of the full text confirms that this is the decisive fragility. The paper contains no modified-gravity evolutionary solar model, no full non-radial pulsation calculation in modified gravity, and no independent calibration of the ξ prior; the quoted interval in Eq. (7) is governed by f'_pol and the chosen ξ range. The stated 'marginalisation range has no practical effect' is true only within the assumed prior, so it does not address the possibility that the proxy itself is wrong in shape or normalization. The concern is not internal inconsistency: the paper is transparent about the approximations and explicitly labels the result an order-of-magnitude estimate. For that reason, the appropriate verdict remains CONDITIONAL rather than REJECT: the proposal is plausible and the method is worth pursuing, but the central number should be treated as an illustration pending validation with realistic modified-gravity solar models. The single concrete check that would settle the concern is to build such models and measure f'_theory directly, comparing it mode-by-mode with f'_pol and then recomputing the constraint.","tokens_in":13580,"tokens_out":5041,"duration_ms":57790,"concrete_test":"Construct MESA solar models with the modified hydrostatic equilibrium (Eq. 3) for at least Y = ±0.01, ±0.03, and ±0.1, matching the same calibration as the Y = 0 reference model, and compute non-radial l = 5–35, n = 6–11 eigenfrequencies with GYRE, ideally without the Cowling approximation. Measure f'_theory(n,l) by numerical differentiation and compare with f'_pol(n,l). If all ratios lie within [0.78, 1.22] and the mode dependence is consistent with a single ξ, recompute Eq. (7); if any ratio falls outside that range, or the recomputed 2σ interval moves by more than a factor of about two, the headline constraint is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The result quoted as Eq. (7) is obtained from a likelihood in which the predicted frequencies are f_theory(Y) = f_theory(0) + f'_pol(n,l)(1+ξ)Y, with ξ = (f'_pol - f'_theory)/f'_pol treated as a single nuisance parameter and marginalized over |ξ| ≤ 0.22. The width of the final Y interval is therefore set almost entirely by f'_pol(n,l), the frequency slope of a polytropic solar model under Eq. (3). The authors do not compute f'_theory with a realistic modified-gravity solar model; they instead estimate |ξ| ≤ 5.2% from a WKB treatment of f̄_theory, and then conservatively widen to 22% by analogy with the 17.5–21.5% frequency offset of the polytrope. This does not test the load-bearing premise. The solar response to the fifth force depends on how the modified hydrostatic equilibrium reshapes density, pressure, and sound speed in a model with real microphysics; the polytrope can misrepresent that response in magnitude, in l/n dependence, or even in sign for some modes. In addition, Eq. (5) assumes a single scalar ξ for all 19 modes; if f'_theory/f'_pol varies mode-by-mode, marginalizing over one scalar does not correctly propagate the unknown. The paper is transparent about this, explicitly stating 'f'_theory is unknown', and the authors call Eq. (7) an order-of-magnitude estimate, so the proposal survives as a proof-of-principle. But the specific 2σ interval is not a measurement of Y unless the proxy is validated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes helioseismology as a new probe of fifth forces in DHOST scalar-tensor theories, modeled through a modified hydrostatic equilibrium equation. Using a polytropic solar model with a linear-in-Y frequency correction, and selecting 19 GONG modes that already agree with a standard-gravity evolutionary model at the 1σ level, the authors derive a 2σ constraint on the fifth-force coupling Y of -1.8e-3 ≤ Y ≤ 1.2e-3 (Eq. 7), claiming an improvement of more than two orders of magnitude over existing white-dwarf bounds. They also translate this into constraints on the cosmological parameters α_H and β_1. The paper explicitly acknowledges that the true frequency response of the Sun to the fifth force is unknown and is treated as a nuisance parameter ξ, so the quoted interval is presented as an order-of-magnitude estimate of the method's constraining power.","tokens_in":13979,"tokens_out":2503,"duration_ms":25977,"significance":"If the central inference were validated, the proposed method would open a genuinely new observational window on modified gravity at stellar scales, with potential constraints on DHOST parameters orders of magnitude tighter than current ones. The paper is transparent about its assumptions, uses public helioseismic data, and combines established codes (MESA, GYRE) with a clear statistical framework. The key strengths are the identification of a concrete observable (acoustic mode frequencies) that responds to the fifth force, and a plausible order-of-magnitude estimate of the achievable sensitivity. However, the scientific value of the specific quoted interval depends entirely on the credibility of the polytropic proxy for the true solar response, which is not established in the manuscript.","major_comments":[{"comment":"The central constraint Eq. (7) rests on replacing the unknown true frequency derivative f'_theory with f'_pol(1+ξ) and marginalizing over |ξ| ≤ 0.22. The paper never computes f'_theory from a solar model with modified gravity; the WKB estimate for ξ is derived from the same polytropic framework and the conservative range is justified by analogy with the polytrope's frequency offset. This does not test the load-bearing premise that the polytrope captures the solar response. If the true response differs in magnitude, sign, or mode dependence, the quoted Y interval is not a valid measurement. A realistic modified-gravity solar model (or at least a mode-dependent systematic treatment) is needed to support Eq. (7).","section":"Modelling and computation, Eq. (5)"},{"comment":"The 19 modes used for the likelihood are selected because they satisfy |f_theory(Y=0) - f_obs| < σ_obs under the standard-gravity evolutionary model, and the paper states that including a broader set of modes causes tension with Newtonian gravity at 2σ. This selection means the analysis is performed on modes that are already consistent with the null hypothesis, so the narrowness of the resulting interval partly reflects the selection criterion rather than the fifth-force sensitivity. The paper does not model the selection effect or its impact on the posterior, making it difficult to interpret Eq. (7) as a rigorous statistical bound.","section":"Constraining power, data selection"},{"comment":"The nuisance parameter ξ is assumed to be a single scalar common to all modes, but the paper's own Fig. 1 shows that f'_pol varies with degree and overtone, and the true response f'_theory could in principle vary mode-by-mode in a different way. Marginalizing over one scalar ξ does not propagate the unknown theoretical error when the ratio f'_theory/f'_pol is not constant across the mode set. The analysis should either allow for mode-dependent systematics or justify why a single scalar suffices.","section":"Eq. (5) and Fig. 1"}],"minor_comments":[{"comment":"In the sentence 'Marginalising over |ξ| ≤22' the value 22 should read 0.22; as written it suggests a marginalization range orders of magnitude larger than intended.","section":"Constraining power, text after Eq. (7)"},{"comment":"The polytropic index is given as n_pol = 3.069 in one place and n_pol = 3.068 in another; the discrepancy should be resolved.","section":"Modelling and computation, Eq. (4) context"},{"comment":"The caption states 'Continuous (dashed) curves correspond to weaker (stronger) gravity with Y > 0 (Y < 0)', but the text earlier says Y > 0 tends to weaken gravity; the caption is consistent, but the wording could be clarified to avoid confusion about the sign convention.","section":"Fig. 1 caption"},{"comment":"Reference [30] has an empty title field; the entry should include the full GONG data description.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is a proof-of-principle with a clearly stated but unvalidated proxy for the solar frequency response. The title and abstract overstate the result as an achieved precision constraint when Eq. (7) is explicitly an order-of-magnitude estimate under an unverified assumption. I would encourage the authors to either compute f'_theory with a realistic modified-gravity solar model or reframe the paper as a sensitivity forecast without quoting a formal confidence interval. The data-selection issue also needs a more careful treatment before the result can be taken as a measurement."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing you should know: this paper proposes helioseismology as a probe of fifth forces in DHOST theories, and that is genuinely new. It does not overclaim: the authors explicitly call Eq. (7) an order-of-magnitude estimate. That framing is accurate, because the central number rests on an unvalidated proxy.\n\nWhat's good: the polytropic scaling analysis is clean, the use of GYRE with public solar models is reproducible in principle, and the treatment of the unknown frequency response via a nuisance parameter ξ is honest. The bound on Y is two orders of magnitude tighter than white-dwarf constraints, and the mapping to (α_H, β_1) is clearly presented. The paper is well written and transparent about its data selection and approximations.\n\nThe soft spot is exactly where the stress-test points: f'_theory is never computed with a modified-gravity solar model. The entire constraint comes from assuming the polytropic slope f'_pol is a good proxy, with a single scalar ξ marginalized over [-0.22, 0.22]. If the real fifth-force response differs mode-by-mode in sign or magnitude, the 2σ interval is not a measurement. The authors acknowledge this, but that does not make the number load-bearing. Also, the 19 modes are selected because they already fit standard gravity within 1σ; the paper honestly reports that including more modes creates tension, which is a red flag for systematics.\n\nStill, this is a proof-of-principle, not a fake result. The idea is plausible and the next step — building a modified-gravity solar model — is clearly identified. I would send this to a serious referee because the community needs this proposal on the record, with the caveats attached. The paper should be published with the framing that Eq. (7) is an illustrative sensitivity estimate, not a final constraint.\n\nIf you work on modified gravity or helioseismology, it's worth a read. I'd bring it to reading group, and I'd probably cite it as the first helioseismic fifth-force constraint, with a caveat.","headline":"A genuinely new idea with an honest but fragile central number: helioseismology as a fifth-force probe deserves serious consideration, but Eq. (7) is a sensitivity estimate, not a measurement.","tokens_in":14447,"tokens_out":1842,"would_cite":true,"duration_ms":19622,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Solar pulsations can tighten fifth-force bounds by two orders of magnitude.","keywords":["helioseismology","modified gravity","fifth force","DHOST scalar-tensor theories","solar oscillations","stellar pulsations","dark energy","Vainshtein screening"],"falsifier":"Build a modified-gravity solar evolution model with the same input physics as the reference model but nonzero $Y$, compute the eigenfrequencies without the Cowling approximation, and compare the derivative $f'_{\\rm theory}$ with $f'_{\\rm pol}$; if the implied $|\\xi|$ exceeds 0.22 for the selected modes, the interval does not follow. A simpler check is to fit all available GONG modes without the $1\\sigma$ selection criterion: if the resulting constraint shifts far outside $[-1.8\\times10^{-3},1.2\\times10^{-3}]$ or disappears, the mode-selection step is driving the claimed sensitivity.","tokens_in":13361,"feed_emoji":"🌞","tokens_out":6645,"duration_ms":68973,"temperature":0.7,"pith_summary":"This paper argues that the Sun's acoustic oscillation spectrum carries a measurable imprint of a fifth force produced by the most general scalar-tensor theories of dark energy, the DHOST theories. Within these theories, the force adds a term proportional to the coupling $Y$ to the stellar hydrostatic equilibrium, slightly reshaping the interior sound-speed profile and thereby shifting the predicted pulsation frequencies. Comparing a selection of 19 GONG solar frequencies with frequencies computed from a polytropic solar model plus a standard-gravity reference model, the paper finds $-1.8\\times10^{-3}\\le Y\\le1.2\\times10^{-3}$ at $2\\sigma$ after marginalising over a nuisance parameter for the unknown solar response. If this holds, it would tighten existing white-dwarf bounds on $Y$ by more than two orders of magnitude and translate into $10^{-3}$-level constraints on the DHOST parameters $\\alpha_H$ and $\\beta_1$. The paper presents helioseismology as a new, local, high-precision test of gravity at astrophysical scales.","feed_headline":"Solar pulsations tighten fifth-force bounds 100-fold","feed_subtitle":"A 19-mode fit to solar oscillation frequencies yields a 2σ interval from -0.0018 to 0.0012, far sharper than white-dwarf limits.","key_machinery":"The argument turns on a linear expansion of the acoustic frequency around $Y=0$: $f_{\\rm theory}(Y)=f_{\\rm pol}(Y)+f'_{\\rm pol}\\,\\xi\\,Y$, where $f_{\\rm pol}$ is computed in a polytropic model with index $n_{\\rm pol}=3.069$ and $\\xi\\equiv(f'_{\\rm pol}-f'_{\\rm theory})/f'_{\\rm pol}$ at $Y=0$ is an unknown nuisance parameter measuring how the true solar response differs from the polytropic proxy. The WKB asymptotic formula $f=(n+l/2+\\alpha)\\bar f$, with $\\bar f$ the inverse round-trip acoustic travel time, supplies an upper bound $|\\xi|\\le5.2\\%$ from the derivative of $\\bar f$ under the fifth force; a deliberately conservative choice widens this to $|\\xi|\\le0.22$. The polytropic frequencies come from solving the non-radial adiabatic pulsation equations under the Cowling approximation, with a zeroth-order microphysics correction calibrated to a standard-gravity solar model.","core_discovery":"The central claim is that helioseismic observations can constrain the fifth-force strength $Y$ in DHOST scalar-tensor theories to $10^{-3}$ accuracy, far beyond current astrophysical bounds. The fifth force modifies the hydrostatic equilibrium through the term $(G Y/4)(d^2m/dr^2)\\rho$, altering the sound-speed profile and hence the acoustic eigenfrequencies; the paper computes these frequencies for polytropic solar models, calibrates the polytrope against a standard-gravity solar evolution model, and compares 19 well-fitting modes with GONG data. Using a linear correction with a nuisance parameter $\\xi$ for the unknown frequency response, it obtains $-1.8\\times10^{-3}\\le Y\\le1.2\\times10^{-3}$ at $2\\sigma$. The paper frames this interval as an order-of-magnitude illustration of helioseismology's constraining power rather than as a complete helioseismic inversion, since a broader set of modes produces tension that would need improved modelling to resolve.","pith_inferences":["A full helioseismic inversion using modified-gravity solar evolution models could convert the paper's mode-selection dependence into an all-mode consistency test, confirming the $10^{-3}$ scale or exposing a background-modelling bias that mimics a fifth force.","The same frequency-shift technique could be applied to other well-observed stars through asteroseismology, where different internal structures might amplify or suppress the fifth-force signature and thereby provide independent checks on the solar result.","Because the quoted interval barely changes when the nuisance range is widened from $|\\xi|\\le5.2\\%$ to $|\\xi|\\le0.22$, the statistical constraint is robust to the unknown response within that range; the decisive uncertainty is whether the true response lies inside that range at all."],"forward_implications":["Existing white-dwarf bounds on the fifth-force coupling ($Y>-0.48$ and $Y<0.18$) would be replaced by an interval roughly two orders of magnitude tighter, making local stellar tests competitive with cosmological probes.","The DHOST parameters $\\alpha_H$ and $\\beta_1$ would be constrained at the $10^{-3}$ level, significantly sharper than the current combined bounds from pulsar and white-dwarf observations.","The constraint relies on selecting only modes whose standard-gravity predictions agree with observations to within $1\\sigma$; including all modes produces tension, indicating that a full helioseismic inversion is needed to separate genuine fifth-force effects from background-modelling artifacts.","The result gives a concrete numerical target for future modified-gravity solar models: non-Cowling pulsation calculations and evolutionary solar models with nonzero $Y$ should either reproduce the interval or reveal the systematic offset."],"supporting_citations":[{"why":"Supplies the mapping from the stellar fifth-force coupling $Y$ to the cosmological parameters $\\alpha_H$ and $\\beta_1$, together with the earlier constraints this paper compares against.","marker":"[14]"},{"why":"Provides the lower bound $Y>-0.48$ from relativistic stars in beyond-Horndeski theories, the baseline the new interval improves on.","marker":"[18]"},{"why":"Provides the white-dwarf upper bound $Y<0.18$ that serves as the comparison point for the claimed improvement.","marker":"[19]"},{"why":"Supplies the WKB asymptotic frequency formula used to estimate the nuisance parameter $\\xi$.","marker":"[22]"},{"why":"Supplies the linear non-radial pulsation equations and boundary conditions used for the eigenfrequency computation.","marker":"[24]"},{"why":"Provides the open-source oscillation solver used to compute the polytropic and evolutionary-model eigenfrequencies.","marker":"[28]"},{"why":"Supplies the observed solar frequencies used in the likelihood.","marker":"[30]"},{"why":"Provides the standard-gravity solar evolution code used to build the reference model and calibrate the polytrope.","marker":"[31]"}],"fun_headline_variants":["Helioseismology pins fifth force to parts per thousand","Helioseismic modes shrink fifth-force bounds 100-fold","Sun's oscillations improve fifth-force constraints 100x","Solar pulsations tighten fifth-force bounds to 10^-3"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the way the real Sun's oscillation frequencies would respond to a fifth force is close to the response of the simplified polytropic model, within the range allowed for the nuisance parameter; if the true response differs more than that, the quoted bound is not a valid measurement of $Y$.","fun_headline_variants_meta":{"raw":{"variants":["Helioseismology pins fifth force to parts per thousand","Helioseismic modes shrink fifth-force bounds 100-fold","Sun's oscillations improve fifth-force constraints 100x","Solar pulsations tighten fifth-force bounds to 10^-3"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000528,"raw_usage":{"total_tokens":2545,"prompt_tokens":939,"completion_tokens":1606,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":555,"completion_tokens_details":{"reasoning_tokens":1537}},"tokens_in":555,"tokens_out":1606,"duration_ms":13571,"temperature":1.0,"reasoning_tokens":1537,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:46:48.601278+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Build a modified-gravity solar evolution model with the same input physics as the reference model but nonzero $Y$, compute the eigenfrequencies without the Cowling approximation, and compare the derivative $f'_{\\rm theory}$ with $f'_{\\rm pol}$; if the implied $|\\xi|$ exceeds 0.22 for the selected modes, the interval does not follow. A simpler check is to fit all available GONG modes without the $1\\sigma$ selection criterion: if the resulting constraint shifts far outside $[-1.8\\times10^{-3},1.2\\times10^{-3}]$ or disappears, the mode-selection step is driving the claimed sensitivity.","supporting_citations":[{"cited_title":"Linear adiabatic stellar pulsation.,","cited_arxiv_id":null,"evidence_quote":"Supplies the WKB asymptotic frequency formula used to estimate the nuisance parameter $\\xi$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the linear non-radial pulsation equations and boundary conditions used for the eigenfrequency computation."},{"cited_title":"GONG is managed by the National Solar Observatory, which is operated by AURA, Inc","cited_arxiv_id":null,"evidence_quote":"Supplies the observed solar frequencies used in the likelihood."}],"review_version":1}