{"id":"9291169b-7130-495c-9c34-327e92e862d8","arxiv_id":"1909.02554","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A cosmological Brans-Dicke model with a fixed cosmological constant fits current data better than standard LambdaCDM and lowers the H0 and sigma8 tensions, while its effective dark energy mimics quintessence near z=0.","lead":"Inside the standard model of cosmology, two key numbers disagree: the local expansion rate versus the cosmic microwave background value, and a measure of how clumpy matter is. This paper asks whether Brans-Dicke gravity, where the strength of gravity slowly changes over time, can make those numbers agree, and it finds a fit that noticeably reduces the disagreement.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The fitted |ω_BD|≈300 is excluded by Cassini-scale bounds unless a screening mechanism exists, but the model contains no scalar potential or other mechanism to provide one.","rationale":"I agree with the reader's weakest_assumption. The paper's own Section 3 text is the clearest statement of the load-bearing condition: the authors knowingly discard local bounds. The fitted parameter is close to GR but not close enough for solar-system tests; the only escape is screening. However, the model as written is standard massless Brans-Dicke (plus Λ), which has no screening: chameleon screening requires a scalar potential, Vainshtein requires derivative self-interactions, and symmetron requires a symmetry-breaking potential. None appears in the action. The authors cite 'other contexts' but do not import any such mechanism. This is a physical-viability concern rather than a mathematical inconsistency; it is appropriate to keep the verdict conditional unless the authors can exhibit a screening mechanism that leaves the cosmological predictions unchanged. Secondary issues (Eq. (10) quoted without derivation; R19 included in the fitted datasets) are noted in the paper and in the reader's rationale; they are less decisive because Eq. (10) is a standard subhorizon result and the R19 inclusion is transparent, whereas the screening assumption is both quantitatively excluded by existing data and structurally unsupported. Therefore no change to the reader's CONDITIONAL verdict.","tokens_in":9987,"tokens_out":10577,"duration_ms":126039,"concrete_test":"Compute the PPN parameter γ=(1+ϵ_BD)/(1+2ϵ_BD) at the posterior mean ϵ_BD=−0.003 (and at the 1σ upper and lower values) and compare |γ−1| with the Cassini bound 2.3×10^−5. Since the action in Section 2 has V(ψ)=0, also check whether any term in the Lagrangian can give the scalar an environment-dependent mass; the absence of such a term means the screening assumption has no dynamical support. If |γ−1| exceeds the bound and no screening term exists, the fitted parameter region is excluded.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on the fitted value ϵ_BD ≈ −0.003 (|ω_BD| ≈ 300) being physically admissible. Section 3 explicitly sets aside solar-system bounds ('we take the wider perspective that the BD theory, when applied to the cosmological level, is not restricted by the bounds obtained in the astrophysical neighborhood') and appeals to screening, citing Avilez & Skordis 2014. But the action used in Section 2 contains no potential V(ψ) and no nonlinear self-interactions; the BD scalar is massless and long-range in this theory. In the static weak-field limit the same model predicts a PPN parameter γ=(1+ϵ_BD)/(1+2ϵ_BD); at ϵ_BD=−0.003 this gives γ−1≈3×10^−3, more than two orders of magnitude above the Cassini bound |γ−1|<2.3×10^−5. No chameleon, Vainshtein, or symmetron mechanism is present to make ω environment-dependent. Thus the fitted region is ruled out by local gravitational tests unless the model is supplemented with a screening sector that is not specified. This is not an internal inconsistency, but it is the most load-bearing assumption: if local bounds apply, the cosmological preference is irrelevant to the viability of the theory.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript analyzes Brans-Dicke gravity with a bare cosmological constant and cold dark matter (BD-ΛCDM) against a combination of SnIa, H(z), BAO, RSD, weak lensing, and the full Planck 2015 CMB likelihood. The authors use two data sets, DS1 and DS2, with the public codes CLASS and MontePython, fitting eight parameters including the BD parameters ε_BD and φ_ini. They report that BD-ΛCDM is preferred over GR-based ΛCDM by ΔDIC ≈ 8.3–9.9 and ΔAIC ≈ 7.7–9.9, with H0 ≈ 71–72 km/s/Mpc, σ8 ≈ 0.80, and an effective dark-energy equation of state w_eff(0) ≈ −0.95 to −0.96, which they interpret as quintessence-like behavior at ≳3σ. The paper concludes that the H0 and σ8 tensions of ΛCDM essentially disappear in this context. The central load-bearing elements are the quoted matter-perturbation equation, the screening assumption that bypasses solar-system bounds on ω_BD, and the interpretation of the fitted BD field as an effective quintessence component.","tokens_in":10214,"tokens_out":5397,"duration_ms":62327,"significance":"If the central claim holds, BD-ΛCDM would be a minimal and economical modification of the concordance model that resolves two major tensions without introducing a new dark-energy fluid. The paper's strengths are its use of standard public pipelines and the full Planck 2015 likelihood, the explicit reporting of marginalized constraints and information-criterion differences, and a clear statement of the screening assumption. These features make the statistical comparison reproducible in principle. However, the quoted growth equation is not derived, the screening mechanism is not specified, and the quintessence signal is a reparameterization of the fitted BD field rather than an independent observable. As a result, the confidence that can currently be placed in the model-selection preference is limited, and the interpretation of the 3σ w_eff result needs to be framed more carefully.","major_comments":[{"comment":"The linear matter-perturbation equation is quoted without derivation, with the details deferred to a future 'more complete presentation.' Equation (10) is the basis for the reported σ8, RSD, WL, and bispectrum fits, so the central model-selection result depends on it. Please provide a derivation or a public implementation that demonstrates how the modified CLASS equations reproduce Eq. (10) in the subhorizon limit, including the treatment of the BD scalar-field perturbations and the (2+4ε)/(2+3ε) coupling factor. As written, the growth sector cannot be independently verified.","section":"§3, Eq. (10)"},{"comment":"The fitted value ε_BD ≈ −0.003 corresponds to |ω_BD| ≈ 300 and to a PPN parameter γ−1 ≈ 3×10^−3, which is more than two orders of magnitude above the Cassini bound. The manuscript explicitly sets aside local bounds and appeals to screening, but the action in §2 contains no potential for ψ, no nonlinear self-interactions, and no chameleon, Vainshtein, or symmetron mechanism that would make ω_BD environment-dependent. The manuscript should either specify a concrete screening mechanism that reconciles the cosmological and local values, or present the analysis as conditional on this assumption and quantify how severely the Cassini bound would restrict the fitted region if it applies.","section":"§3, screening of solar-system bounds"},{"comment":"The claim that BD-ΛCDM 'mimics quintessence at more than 3σ' is a restatement of the fitted BD field, not an independent constraint: w_eff is constructed from the same fitted φ that enters the background and growth equations, so the 3σ deviation from −1 is a derived parameter rather than a new observable. In addition, the functions f1 and f2 in Eq. (8) are not given explicitly, and the statement that they are 'numerically negligible' is asserted without a bound. Please provide explicit expressions and a quantitative assessment of the dropped terms near z = 0, and describe the reported significance as a constraint on the derived w_eff rather than as an independent detection of quintessence.","section":"§2.2, Eq. (8) and §4.2, w_eff(0)"},{"comment":"The numerical solution fixes φ'(z_ini) = 0 at z_ini = 10^14. This is a nontrivial theoretical input: the BD scalar is dynamical, and imposing a zero initial velocity at a finite initial epoch is a prior on the solution space. The paper reports no test of the sensitivity of H0, σ8, or ΔDIC/ΔAIC to the choice of z_ini and to the initial derivative. Please add a convergence/sensitivity test (for example, varying z_ini by an order of magnitude and allowing a non-zero initial derivative) to show that the model-selection result is not driven by this boundary condition.","section":"§4.2, initial conditions"}],"minor_comments":[{"comment":"There is a typo, 'Hearafter', that should read 'Hereafter'.","section":"§2.1, before Eq. (4)"},{"comment":"The sentence 'The two functions f1,2 need not be specified here' is unsatisfactory in a Letter whose central claim is based on w_eff; at minimum, the explicit expressions should be given in an appendix or in a companion public file.","section":"§2.2, Eq. (8)"},{"comment":"The compressed notation for the 1σ, 2σ, and 3σ intervals in the Table 1 caption is difficult to parse; please present the 68% intervals in the main table and move the extended intervals to a supplementary table.","section":"§4.2 and Table 1"},{"comment":"The Hildebrandt et al. entry is given only as an arXiv number; please add the journal reference if published.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of the journal and addresses a timely question. The main concern for the editor is the verification status of the perturbation implementation: the quoted growth equation is central to the σ8 and LSS constraints, and deferring its derivation to 'elsewhere' is not acceptable for a paper claiming model selection by ΔDIC/ΔAIC. I would require a derivation, an appendix, or a code-release note before reconsidering. The self-citation pattern is heavy but not disqualifying for a Letter extending prior work by the same group."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a competent, well-documented full-data fit of BD-LCDM, and it does soften the H0 and sigma8 tensions with a moderate DIC/AIC preference over LCDM. But the model's fitted |omega_BD| about 300 is only viable if solar-system bounds are switched off, and the paper has no screening mechanism to justify that switch. Treat it as a phenomenological cosmology paper, not as a resolution of the tensions.\n\nWhat is actually new: earlier papers from this group used approximate power-law solutions; this one solves the background and perturbation equations exactly and runs them against the full Planck 2015 likelihood together with SnIa, BAO, RSD, WL, and R19. That is a real step. The perturbation equation (quoted from Boisseau et al.) reduces to LCDM in the GR limit, the MCMC and model-selection tools are public and standard, and the residual H0 tension (1.1-1.8 sigma) and sigma8 about 0.80 are reported honestly.\n\nSoft spots in proportion. The biggest is explicitly acknowledged in Section 3: the authors 'take the wider perspective' that cosmological BD is not restricted by astrophysical-neighborhood bounds, citing screening. But the action they use has no scalar potential and no self-interactions; the BD field is massless. In the static weak-field limit the same theory gives gamma-1 about 3e-3 at the fitted epsilon_BD, two orders of magnitude above Cassini. So the screened escape route is an assumption, not a property of the model they wrote down. A referee should ask for the screening sector or a citation to a concrete mechanism in this exact model.\n\nSecond, the 'quintessence at more than 3 sigma' claim is a reparametrization of the fitted BD field parameters, not an independent detection. The significance inherits from the same chain that fixed epsilon_BD and phi_ini. It is a nice interpretive frame, not a new observable.\n\nThird, the perturbation equation is quoted without derivation, with details deferred to a future paper. Acceptable in a letter, but a serious referee should want the calculation or a fuller reference.\n\nFourth, because R19 is in the dataset, the reduced H0 tension is partly by construction. The model is flexible enough to pull H0 up while keeping sigma8 low; that is the real result, stated without overclaiming.\n\nWho this is for: modified-gravity and tension people. I would bring it to our reading group, and I would cite the exact numerical treatment if I worked in BD cosmology. The screening issue does not kill the paper as a statistical/physical analysis, but it prevents using the result as evidence for BD gravity at cosmological scales without additional physics. Deserves a serious referee; I would send it out, with the screening point pressed.","headline":"A technically solid BD-LCDM fit that softens H0 and sigma8 tensions, but its viability leans on an unproven screening assumption that lets |omega_BD|~300 evade Cassini bounds; deserves a real referee.","tokens_in":10791,"tokens_out":3925,"would_cite":true,"duration_ms":45051,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["98.80.-k","04.50.Kd"],"model":"deepseek-v4-flash","headline":"A Brans-Dicke gravity model with a bare cosmological constant and no extra dark energy fluid fits the combined cosmological data better than GR-based ΛCDM and dissolves the H0 and σ8 tensions.","keywords":["Brans-Dicke gravity","cosmological constant","Hubble tension","sigma8 tension","effective dark energy equation of state","quintessence","modified gravity","cosmological parameter estimation"],"falsifier":"A cosmological-scale constraint that forced $|\\omega_{\\rm BD}|\\gtrsim10^4$ at all epochs, or a null measurement of $\\dot G/G$ at the level of $10^{-14}\\,\\mathrm{yr}^{-1}$, would exclude the fitted $\\epsilon_{\\rm BD}\\simeq -0.003$ and remove the claimed tension relief.","tokens_in":9714,"feed_emoji":"🌌","tokens_out":7200,"duration_ms":72923,"temperature":0.7,"pith_summary":"The paper argues that the standard cosmological model's two most persistent headaches — the higher local value of the Hubble constant and the surplus structure-growth amplitude — both relax when gravity is described by Brans-Dicke theory with a bare cosmological constant and cold dark matter (BD-ΛCDM). Fitting supernova, H(z), BAO, redshift-space-distortion, weak-lensing, and full Planck 2015 CMB data, the model is preferred over the general-relativistic ΛCDM by standard information criteria. The best fit raises H0 from about 68.6 to 71–72 km/s/Mpc, cutting the tension with local distance-ladder measurements to about 1–2σ, while σ8 stays near 0.80. From the GR viewpoint the extra degrees of freedom look like a dynamical dark energy component whose equation of state is quintessence-like ($w_{\\rm eff}(0)\\simeq -0.96$) at above 3σ significance near the present epoch. If right, no new dark energy fluid is needed: a mildly time-varying gravitational coupling and the cosmological constant suffice to smooth out both tensions.","feed_headline":"Brans-Dicke gravity dissolves Hubble and sigma-8 tensions","feed_subtitle":"A bare cosmological constant plus a mildly varying G fits the data and raises H0 toward local measurements.","key_machinery":"The carrying object is the dimensionless Brans-Dicke field $\\phi(t)=G_N/G(t)$, with $\\epsilon_{\\rm BD}=1/\\omega_{\\rm BD}$ measuring the departure from GR. Rewriting the BD field equations as Friedmann equations produces an effective BD-fluid energy density and pressure, whose combined equation of state with the vacuum term is $w_{\\rm eff}(t)$; near $z=0$ this takes the approximate form $w_{\\rm eff}\\simeq -1 + \\frac{\\Delta\\phi\\,\\Omega_{m0}}{1-\\Omega_{m0}}(1+z)^3$. On the perturbation side, the linear matter-density contrast obeys a modified growth equation in which the effective gravitational constant is $(G_N/\\phi)(2+4\\epsilon_{\\rm BD})/(2+3\\epsilon_{\\rm BD})$. These two pieces — background equation of state and growth equation — let a single parameter $\\epsilon_{\\rm BD}$ shift both $H_0$ and $\\sigma_8$ in the observed directions.","core_discovery":"The central claim is that BD-ΛCDM — Brans-Dicke gravity with a constant vacuum term and no interaction between the scalar field and matter — provides a better fit to the combined cosmological dataset than the GR-based ΛCDM, with the H0 and σ8 tensions essentially disappearing. In the best fit, $\\epsilon_{\\rm BD}\\simeq -0.003$ (i.e. $\\omega_{\\rm BD}\\simeq -300$), $\\phi(0)\\simeq 0.90$, so G is slightly larger and growing at present, $\\dot G/G\\simeq +3\\times10^{-13}\\,\\mathrm{yr}^{-1}$. The model yields $H_0\\simeq 71$–$72\\,\\mathrm{km\\,s^{-1}\\,Mpc^{-1}}$, in line with the local distance ladder at the 1–2σ level, and $\\sigma_8\\simeq 0.80$, in between Planck and weak-lensing values. In the GR frame the BD field contributes an effective dark-energy fluid with a time-varying equation of state that near $z=0$ takes $w_{\\rm eff}\\simeq -0.96$, mimicking quintessence at more than 3σ. The DIC and AIC differences between the models lie between 5 and 10, which the authors read as strong evidence for BD-ΛCDM.","pith_inferences":["The paper leaves implicit that the same scalar field can in principle absorb both tensions at once, which suggests a common physical origin — a mildly running gravitational coupling — rather than two independent systematics.","If the model is right, future BAO and weak-lensing surveys should see the effective dark energy equation of state deviate from $-1$ at low redshift in the direction of quintessence, with the deviation growing roughly as $(1+z)^3$.","A null detection of the predicted present-day growth of $G$, $\\dot G/G\\simeq 3\\times10^{-13}\\,\\mathrm{yr}^{-1}$, with upcoming lunar-laser-ranging or atomic-clock experiments would put the model under direct observational pressure.","The screening assumption could be tested cross-model: if the cosmological Brans-Dicke parameter really is as small as $|\\omega_{\\rm BD}|\\simeq 300$, comparable tension relief should appear in other screened modified-gravity theories fitted to the same data."],"forward_implications":["The tension with the local distance-ladder value of $H_0$ drops from about 4.4σ in ΛCDM to roughly $1.8\\sigma$ (DS1) and $1.1\\sigma$ (DS2).","The mass-fluctuation amplitude stays at $\\sigma_8\\simeq 0.80$, sitting between the higher Planck value and the lower cosmic-shear value and thereby easing the $\\sigma_8$ tension.","Near the present epoch the effective dark energy equation of state is $w_{\\rm eff}\\simeq -0.96$, a quintessence-like signal at $\\gtrsim 3\\sigma$, which the authors describe as a smoking gun of the underlying BD dynamics.","The model is favored over ΛCDM by both DIC and AIC with differences between 5 and 10, which standard model-selection practice treats as strong evidence.","The fitted value of the Brans-Dicke parameter, $|\\omega_{\\rm BD}|\\simeq 300$, keeps the model close enough to GR to be phenomenological while predicting a mild, positive present-day time variation of $G$."],"supporting_citations":[{"why":"Defines the scalar-tensor gravity theory whose cosmological consequences are tested here.","marker":"Brans & Dicke (1961)"},{"why":"Supplies the linear matter perturbation equation used to model structure formation in BD gravity.","marker":"Boisseau et al. (2000)"},{"why":"Provides the local distance-ladder measurement of H0 (R19) that creates the tension the model is designed to relieve.","marker":"Riess et al. (2019)"},{"why":"Provides the full Planck 2015 TT+lowP+lensing likelihood used as the CMB dataset and the σ8 baseline.","marker":"Planck Collaboration XIII (2015)"},{"why":"Provides the weak-lensing data that anchor the low-σ8 end of the structure-growth tension.","marker":"Hildebrandt et al. (2018)"},{"why":"Supports the interpretation that the Brans-Dicke parameter at cosmological scales may evade solar-system bounds via screening.","marker":"Avilez & Skordis (2014)"},{"why":"Provides the earlier dataset compilation and fitting setup that DS1 and DS2 extend.","marker":"Solà, Gómez-Valent & de Cruz Pérez (2019)"},{"why":"Supplies the Pantheon+MCT supernova sample used in the DS1 fit.","marker":"Scolnic et al. (2018)"}],"fun_headline_variants":["Brans-Dicke gravity with Λ dissolves Hubble and σ8 tensions","A bare Λ plus Brans-Dicke scalar removes H0 and σ8 tensions","Mildly varying G resolves the H0 and σ8 tensions","Brans-Dicke scalar with Λ smooths out Hubble and σ8 tensions","When gravity varies slightly, H0 and σ8 align with data"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole fit rests on assuming the Solar System bound on the Brans-Dicke parameter, $|\\omega_{\\rm BD}|>40000$, does not apply at cosmological scales, allowing $|\\omega_{\\rm BD}|\\simeq 300$ through screening.","fun_headline_variants_meta":{"raw":{"variants":["Brans-Dicke gravity with Λ dissolves Hubble and σ8 tensions","A bare Λ plus Brans-Dicke scalar removes H0 and σ8 tensions","Mildly varying G resolves the H0 and σ8 tensions","Brans-Dicke scalar with Λ smooths out Hubble and σ8 tensions","When gravity varies slightly, H0 and σ8 align with data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001206,"raw_usage":{"total_tokens":5002,"prompt_tokens":1015,"completion_tokens":3987,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":631,"completion_tokens_details":{"reasoning_tokens":3891}},"tokens_in":631,"tokens_out":3987,"duration_ms":29477,"temperature":1.0,"reasoning_tokens":3891,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:47:01.838607+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A cosmological-scale constraint that forced $|\\omega_{\\rm BD}|\\gtrsim10^4$ at all epochs, or a null measurement of $\\dot G/G$ at the level of $10^{-14}\\,\\mathrm{yr}^{-1}$, would exclude the fitted $\\epsilon_{\\rm BD}\\simeq -0.003$ and remove the claimed tension relief.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the scalar-tensor gravity theory whose cosmological consequences are tested here."},{"cited_title":"G., et al","cited_arxiv_id":null,"evidence_quote":"Provides the local distance-ladder measurement of H0 (R19) that creates the tension the model is designed to relieve."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the full Planck 2015 TT+lowP+lensing likelihood used as the CMB dataset and the σ8 baseline."},{"cited_title":"M., et al","cited_arxiv_id":null,"evidence_quote":"Supplies the Pantheon+MCT supernova sample used in the DS1 fit."}],"review_version":1}