{"id":"fe99b778-b8ba-4379-a076-69acaaa4649a","arxiv_id":"2502.03190","paper_version":2,"verdict":"REJECT","confidence":"MODERATE","novelty_score":2.0,"correctness_risk":"high","formal_verification":"none","parameter_count":6,"one_line_summary":"A fit of Hu-Sawicki f(R) gravity with a fitted neutrino interaction to CMB, lensing, BAO, CC, and Pantheon data returns H0 near 70 km/s/Mpc, partially easing but not resolving the Hubble and S8 tensions.","lead":"This paper fits a modified gravity model, f(R) gravity with a neutrino interaction, to five cosmological datasets and reports that the model pulls the inferred Hubble constant from 67.4 to about 70 km/s/Mpc, between the Planck and local values. The claimed easing of the Hubble and S8 tensions is a consequence of fitting extra parameters, including an interaction strength, rather than an independent prediction.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The neutrino-coupled model treats Γ both as a derived Hu-Sawicki quantity (Eq. 23) and as an independent free parameter, so the fitted H0 constraints may not describe the claimed theory.","rationale":"The reader identified the missing implementation link as the weakest assumption: without code or a mapping from the phase-space equations to the CLASS likelihood, the numerical constraints cannot be attributed to the stated theory. I agree that this is a serious concern, but I find a more specific and more decisive internal inconsistency: Γ is simultaneously defined as a derived quantity (Eq. 23) and fitted as an independent free parameter (Section VII). This is not merely a missing reproducibility detail; it is a contradiction between the theoretical model and the analysis performed. If Γ is derived from Hu-Sawicki f(R), then it cannot be freely varied, and the reported posterior on Γ has no meaning; if it is free, then Eqs. (20)–(23) do not define the model that was fitted. Either way, the central claim—that neutrino coupling within perturbed f(R) gravity alleviates the Hubble tension—is not supported by the presented constraints. Because this concern reinforces rather than changes the reader's rejection, the verdict remains REJECT, expressed as UNCHANGED. The concrete test isolates the contradiction: evaluating Eq. (23) with the reported best-fit parameters and background will show whether Γ(N) can plausibly be a constant equal to 0.64; if not, the free-Γ analysis is a different model. This is one specific, feasible check that would settle whether the concern lands.","tokens_in":22784,"tokens_out":4752,"duration_ms":44500,"concrete_test":"Compute Γ(N) from Eq. (23) using the best-fit Hu-Sawicki parameters (e.g., c1 = 1.15 × 10^-3, c2 = 6.55 × 10^-5, n = 4) and the best-fit background H(N) reported in Tables IV/VI, then plot Γ(N) over the fitted redshift range. If Γ(N) varies by more than its own mean or is inconsistent with the reported constant Γ = 0.64 ± 0.14 at 68% confidence, the free-constant treatment is invalid and the quoted H0 constraints do not represent the derived model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline claim—that neutrino-coupled perturbed f(R) gravity reduces the Planck-2018/R22 tension to 0.74σ/0.95σ—depends entirely on the fitted H0 and S8 values in Tables V and VII. Those fits treat Γ, the neutrino–gravity coupling, as an independent constant (Section VII: Γ = 0.64 ± 0.14) and place it in the autonomous equation for η8 (Eq. 34). But the theory section defines Γ as u^μ∇_μ f_R (Eq. 20) and then states it is 'explicitly written in terms of the Hu-Sawicki model' as Γ = d f_R/dN (Eq. 23), a function of the Ricci scalar and background dynamics—not a free parameter. A derived quantity cannot be varied independently; if Γ is fixed by (c1, c2, n) and H(N), then the posterior Γ = 0.64 ± 0.14 is meaningless, and the H0 constraints are not fits of the stated model. If Γ is instead meant as a free interaction coefficient, then Eqs. (20)–(23) do not define it and the connection to Hu-Sawicki f(R) is broken. In addition, no modified Boltzmann solver, parameter file, or mapping between the η_i phase-space variables (24)–(39) and CLASS perturbation variables is provided, so the numerical results cannot be checked against the claimed equations. Either the free-Γ fit is a different model, or the derived model cannot produce the quoted constraints; both possibilities remove the support for the central claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies two extensions of ΛCDM, perturbed f(R) gravity and perturbed f(R) gravity coupled to neutrinos, fitting them to CMB, BAO, cosmic chronometer, Pantheon, and lensing data with MontePython/CLASS. The central claim is that the coupled model shifts H0 toward the local R22 value, reducing the tension with Planck 2018 to 0.74σ and with R22 to 0.95σ, and also reduces the S8 tension relative to Planck, KiDS-1000, and DES-Y3. Parameter tables for H0, S8, Γ, and neutrino mass limits are reported for several dataset combinations.","tokens_in":23234,"tokens_out":4058,"duration_ms":37935,"significance":"If the numerical results were reliable, the conclusion that a neutrino coupling in perturbed f(R) gravity can simultaneously ease the Hubble and S8 tensions would be of interest to the modified-gravity and cosmological-tensions communities. The paper assembles a broad set of public datasets and presents parameter tables in a transparent format. However, the central claim is currently undermined by an internal inconsistency in the treatment of Γ, by the absence of any demonstration that the claimed equations are what the likelihood code actually solves, and by inconsistencies between the quoted tension values and the tabulated errors. No machine-checked proofs, modified Boltzmann code, chain files, or parameter files are provided, so the numerical constraints cannot be independently verified. For these reasons the significance of the result, as presented, is not established.","major_comments":[{"comment":"The parameter Γ is defined in Eq. (20) as Γ = u^μ∇_μ f_R and in Eq. (23) is stated to be explicitly written in terms of the Hu-Sawicki model as Γ = d f_R/dN, i.e., a derived function of the Ricci scalar and background dynamics. In contrast, Section VII treats Γ as an independent free constant and reports constraints such as Γ = 0.64 ± 0.14. A derived quantity cannot be varied independently of the Hu-Sawicki parameters (c1, c2, n) and the background H(N); either the reported fits are not fits of the stated model, or Eq. (23) is not part of the fitted model and the connection to Hu-Sawicki f(R) is broken. This ambiguity directly affects the headline claim that neutrino physics in f(R) gravity alleviates the Hubble tension.","section":"§IV and §VII, Eqs. (23), (34)"},{"comment":"The manuscript states that the likelihood analysis is performed with MontePython and CLASS (refs. [111–113]), but it never shows that the phase-space system (24)–(39), the Hu-Sawicki form (22), or the neutrino interaction (20)–(21) are implemented in the Boltzmann solver. Standard CLASS does not contain this autonomous system, and no modified code, parameter files, or mapping between the η_i variables and CLASS perturbation variables is given. Without such a mapping, the quoted H0 and S8 constraints cannot be checked against the claimed theory, so the numerical results do not currently constrain the model presented in Sections II–IV.","section":"§VI and §V, Eqs. (24)–(39)"},{"comment":"The tension values quoted in Table VII are not consistent with the quoted errors. For the CMB+All row, H0 = 70.46 ± 2.01 km/s/Mpc gives a tension with Planck 2018 (67.4 ± 0.5) of (70.46 − 67.4)/sqrt(2.01^2 + 0.5^2) ≈ 1.48σ, not 1.42σ, and with R22 (73.5 ± 1.04) of (73.5 − 70.46)/sqrt(2.01^2 + 1.04^2) ≈ 1.34σ, not 1.12σ. Similar discrepancies appear in the CMB+Lensing row (0.77σ versus 0.74σ for Planck; 1.12σ versus 0.95σ for R22). Since the abstract's claim of reduction to 0.74σ and 0.95σ is based on these numbers, the headline quantitative result is not supported by the tabulated data.","section":"Table VII and abstract"},{"comment":"The neutrino-mass section introduces η9 in Eq. (53) but the autonomous system (24)–(39) defines only η1 through η8, and no evolution equation for η9 is provided. In addition, Eq. (55) states that substituting ρν into 'Equation (1)' yields the expression for η9, but the relevant definition is Eq. (53), not the action Eq. (1). This makes the connection between the reported Σm_ν constraints and the coupled f(R)+neutrino system unclear and should be corrected before the neutrino-mass bounds can be assessed.","section":"§VII, Eqs. (53)–(55)"}],"minor_comments":[{"comment":"The text refers to 'Table VII' and 'Table IX' in places that do not match the printed table captions; for example, the Hubble-tension table for the coupled model is Table VII, while the S8 tables are Tables VIII and IX. The cross-references in Sections VIII and IX need to be harmonized with the actual table numbers.","section":"Table numbering and cross-references"},{"comment":"The action (17) includes L_m and L_int, and the text mentions L_ν in the same sentence; it should be clarified whether L_ν is part of L_m or a separate term, since the subsequent field equation (18) contains only T^(m) and T^(int).","section":"Section III, Eq. (17)"},{"comment":"The background equations contain notation that is not defined at first use, such as the meaning of the prime in Eq. (3) and the role of c_s^2 before the dust-dominated assumption is introduced; defining these quantities in place would improve readability.","section":"Section II, Eqs. (3)–(5)"},{"comment":"Several BAO entries in Table I are listed without uncertainties (e.g., the rows at z = 0.38 and z = 0.51), which is inconsistent with their use in a likelihood; the source of each value and error should be stated.","section":"Section V, Table I"},{"comment":"Figure 1 is described as showing the comparison of Σm_ν and H0 (top) and Γ and H0 (low), but the axes and contours are not legible in the text version; a higher-resolution figure with axis labels and caption definitions would be needed.","section":"Section VII, Fig. 1"}],"recommendation":"reject","confidential_remarks":"The manuscript draws heavily on a series of previous papers by the same author (refs. [98–101]) for the theoretical setup and for comparisons, but the present text does not provide enough detail to verify that the numerical pipeline matches the stated equations. The derived-versus-free Γ inconsistency and the missing code-to-theory mapping are not cosmetic; they remove the evidential basis for the central claim. These issues would require substantial new derivation, code release, and refitting to address, which is beyond a routine revision. I therefore recommend rejection, although I would encourage the authors to resubmit a version that clearly specifies the model actually implemented in the likelihood code and provides reproducible chains."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nWhat you need to know: this is a re-run of Yarahmadi's existing neutrino-coupled Hu-Sawicki f(R) model with extra dataset combinations and an S8 analysis. That is the whole new content. The model is not new here; the central mechanism—Γ shifting H0 upward—is already in refs [98]–[101] according to the paper's own text.\n\nWhat it does well: it is honest that tensions remain; the tables make dataset dependencies clear; the BAO/CC/Pantheon compilation is standard and the priors are stated; the S8 comparison with Planck, KiDS-1000, and DES-Y3 is a useful addition. On its face, the coupled model pulls H0 to ~70 km/s/Mpc and lowers S8 to ~0.785, which is the direction needed.\n\nThe soft spots are serious. First, the theory defines Γ as u^μ∇_μ f_R and writes it explicitly as d f_R/dN, a background-derived function of the Hu-Sawicki parameters. Yet Section VII fits Γ as an independent constant, Γ=0.64±0.14. Either the fitted model is not the model defined by Eqs. (20)–(23), or the derived relation is being ignored. That is load-bearing, not cosmetics.\n\nSecond, the numerical results cannot be checked. The phase-space system (24)–(39) is not shown to be what the likelihood code solves. No modified CLASS or MGCLASS implementation, no parameter files, no chain files. The text cites MontePython and CLASS, but a standard CLASS run does not solve these equations. Without that mapping, the quoted H0 and S8 constraints do not demonstrably constrain the claimed theory.\n\nThird, the tension numbers are sloppy. In Table VII, CMB+Lensing gives H0=69.82±3.12; combined with R22=73.5±1.04 that is about 1.1σ, not 0.95σ. The 0.74σ Planck tension is about right, but the abstract and conclusion cherry-pick the most favorable lines while the CMB+All line shows 1.42σ/1.12σ. The S8 tensions also treat asymmetric errors as symmetric, which is minor but adds to the impression of post-hoc selection.\n\nFinally, the \"alleviation\" is achieved by fitting Γ, c1, c2, and Σmν to the same data used to define the tension. That is not circular in a fatal way—parameter fitting is how these models work—but it means the result is a parameter scan, not a prediction. The claim that neutrino physics plays a significant role rests on a fitted constant, and there is no same-data ΛCDM baseline or model-selection statistic to show the extra parameters buy anything.\n\nAll that said, the underlying idea—interacting neutrinos in f(R) gravity, constrained by CMB+BAO+CC+Pantheon+lensing—is not crazy, and the author cites his own prior work appropriately. The problem is execution and transparency.\n\nMy recommendation: as submitted, I would desk reject. If the author supplies the modified Boltzmann code, shows that the likelihood actually solves Eqs. (24)–(39), and either drops the free-Γ fit or changes the theory to justify it, the S8 part could be worth a proper review. But this version does not deserve referee time.","headline":"Reuses the author's prior neutrino-coupled f(R) model and adds S8 constraints, but the free-Γ fit contradicts the derived-Γ theory and no code connects the equations to the quoted likelihood results.","tokens_in":23726,"tokens_out":3885,"would_cite":false,"duration_ms":37707,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.50.Kd","98.80.-k","95.36.+x","14.60.Pq"],"model":"deepseek-v4-flash","headline":"Coupling neutrinos to perturbed Hu–Sawicki $f(R)$ gravity moves the fitted Hubble constant to about $70$ km/s/Mpc and lowers $S_8$, leaving residual tensions near one sigma.","keywords":["Hubble tension","f(R) gravity","neutrino coupling","Hu-Sawicki model","S8 tension","modified gravity","cosmological parameter estimation"],"falsifier":"Fix $\\Gamma=0$ and re-run the same likelihood pipeline on the same datasets; the recovered $H_0$ and $S_8$ should match the paper's uncoupled $f(R)$ tables, and if they instead reproduce the coupled-model values, the attribution of tension relief to neutrino coupling is an implementation artifact.","tokens_in":22528,"feed_emoji":"🌌","tokens_out":20585,"duration_ms":160233,"temperature":0.7,"pith_summary":"This paper tries to establish that adding a neutrino–gravity interaction to perturbed $f(R)$ gravity—in the Hu–Sawicki form—raises the inferred Hubble constant from the CMB-preferred value to about $70$ km/s/Mpc, easing the roughly $4\\sigma$ mismatch with local distance-ladder measurements, while also lowering the clustering amplitude $S_8$ toward the values preferred by weak-lensing surveys. It compares the same model with and without a single neutrino interaction parameter $\\Gamma$ on five dataset combinations built from CMB, lensing, BAO, cosmic chronometers, and supernova data. In the coupled model, the residual tension with the CMB-based $H_0$ can drop to $0.74\\sigma$ and the tension with the local value to $0.95\\sigma$, and the $S_8$ tension with weak-lensing priors falls below $1\\sigma$ for several combinations. The paper's conclusion is not that the tensions disappear, but that the neutrino-coupled model aligns the datasets better than $f(R)$ gravity alone.","feed_headline":"Neutrino-coupled f(R) gravity trims Hubble and S8 tensions","feed_subtitle":"Adding a neutrino interaction moves H0 to about 70 and trims residual tensions toward one sigma.","key_machinery":"The machinery is the Hu–Sawicki $f(R)$ function, $f(R) = -m^2 c_1 (R/m^2)^n / (c_2 (R/m^2)^n + 1)$ with $n=4$, together with the neutrino interaction $\\Gamma = u^\\mu\\nabla_\\mu f_R$ and $Q_\\nu=-\\Gamma\\rho_\\nu$. The paper rewrites the perturbed field equations as a first-order autonomous system in the eight phase-space variables $\\eta_1,\\dots,\\eta_8$ (normalized perturbations of the metric potentials, curvature, matter density, and neutrino density), evolving in $N=\\ln a$ according to equations (25)–(34). That system is what the numerical likelihood calculation integrates; the interaction parameter $\\Gamma$ and the total neutrino mass are the extra degrees of freedom that shift $H_0$ upward and $S_8$ downward relative to the uncoupled model.","core_discovery":"The central claim, on the paper's own terms, is that a neutrino interaction term in perturbed $f(R)$ gravity acts as a lever on both the expansion rate and the growth of structure. The paper defines the coupling as $Q_\\nu=-\\Gamma\\rho_\\nu$ with $\\Gamma=u^\\mu\\nabla_\\mu f_R$, treats $\\Gamma$ as a free constant alongside the neutrino mass sum, and finds that including it shifts $H_0$ upward and $S_8$ downward relative to the uncoupled model. Across all dataset combinations the coupled model returns $H_0$ between $69.82$ and $70.57$ km/s/Mpc, and for the CMB+lensing combination the residual tensions are $0.74\\sigma$ with the early-universe reference and $0.95\\sigma$ with the local distance ladder. For the full dataset the model gives $S_8=0.785\\pm0.046$, with tensions of $1.13\\sigma$ against the CMB reference and $0.70\\sigma$ and $0.23\\sigma$ against the two weak-lensing priors. The paper reads this as evidence that modified gravity combined with neutrino physics can adjust both the background expansion and matter clustering in a way the standard model cannot.","pith_inferences":["The paper does not report an information-criterion comparison; adding $\\Gamma$ and loosening the neutrino mass adds parameters, so the improved agreement might not survive a penalty for model complexity.","Treating $\\Gamma$ as a constant is the simplest parametrization; a redshift-dependent interaction would trace a different $H_0$–$S_8$ trajectory that future BAO and weak-lensing data could distinguish.","Because $\\Gamma$ and the neutrino mass sum are fitted simultaneously, a laboratory measurement that pins down the mass sum would sharpen the allowed range of $\\Gamma$.","The same phase-space system could be used to predict growth observables such as $f\\sigma_8$; comparing those predictions to redshift-space distortion measurements at low redshift would test the model independently of the $H_0$ fit."],"forward_implications":["If the model is right, a single mechanism can move $H_0$ toward late-time measurements without wrecking the CMB fit, so the Hubble tension no longer has to be blamed entirely on systematics in one of the two probes.","The same coupling suppresses $S_8$, addressing the growth-structure tension simultaneously; with the full dataset the residual tension with the stronger weak-lensing prior is $0.23\\sigma$.","The full dataset fixes the interaction parameter at $\\Gamma = 0.64 \\pm 0.14$, meaning the data prefer a nonzero energy exchange between neutrinos and the modified gravity sector.","The full dataset also tightens the neutrino mass sum to $\\sum m_\\nu < 0.119$ eV at 95% confidence, consistent with the oscillation lower bound and more restrictive than the CMB-only bound.","Since the uncoupled $f(R)$ model leaves tensions in the $1.1\\sigma$–$1.5\\sigma$ range, the comparison isolates the neutrino coupling as the ingredient doing most of the work."],"supporting_citations":[{"why":"It supplies the early-universe reference values of $H_0$ and $S_8$ against which the Hubble and growth tensions are computed.","marker":"[3]"},{"why":"It supplies the local distance-ladder $H_0$ that defines the late-time side of the Hubble tension.","marker":"[4, 5]"},{"why":"It provides the CMB temperature and polarization likelihood used in the fits.","marker":"[80]"},{"why":"It provides the CMB lensing likelihood and the reference $S_8$ value used for comparison.","marker":"[82]"},{"why":"It supplies a weak-lensing $S_8$ measurement used as a Gaussian prior and comparison target.","marker":"[21]"},{"why":"It supplies a second weak-lensing $S_8$ measurement used as a Gaussian prior and comparison target.","marker":"[22]"},{"why":"It is the source of the $f(R)$+neutrino action, the $\\Gamma$ coupling prescription, and the Hu–Sawicki expressions that define the model.","marker":"[101]"},{"why":"It is the MCMC sampler that produces the posterior parameter constraints.","marker":"[111, 112]"},{"why":"It is the Boltzmann solver that computes the theoretical spectra against which the likelihoods are evaluated.","marker":"[113]"},{"why":"It supplies the supernova sample included in the comprehensive dataset combination.","marker":"[78]"}],"fun_headline_variants":["Neutrino-coupled f(R) gravity cuts Hubble and S8 tension","Neutrino coupling shifts H0 and S8 in f(R) gravity","f(R) gravity with neutrinos eases H0-S8 discrepancy","Neutrino-modified f(R) trims Hubble and structure tensions"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the numerical code used in the likelihood analysis actually solves the paper's phase-space equations for the Hu–Sawicki model with the neutrino interaction term, since the paper supplies the equations but not the code or a detailed mapping from those equations to the chi-squared likelihoods.","fun_headline_variants_meta":{"raw":{"variants":["Neutrino-coupled f(R) gravity cuts Hubble and S8 tension","Neutrino coupling shifts H0 and S8 in f(R) gravity","f(R) gravity with neutrinos eases H0-S8 discrepancy","Neutrino-modified f(R) trims Hubble and structure tensions"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000752,"raw_usage":{"total_tokens":3402,"prompt_tokens":1058,"completion_tokens":2344,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":674,"completion_tokens_details":{"reasoning_tokens":2265}},"tokens_in":674,"tokens_out":2344,"duration_ms":16629,"temperature":1.0,"reasoning_tokens":2265,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T05:34:59.678636+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fix $\\Gamma=0$ and re-run the same likelihood pipeline on the same datasets; the recovered $H_0$ and $S_8$ should match the paper's uncoupled $f(R)$ tables, and if they instead reproduce the coupled-model values, the attribution of tension relief to neutrino coupling is an implementation artifact.","supporting_citations":[{"cited_title":"Yarahmadi et al","cited_arxiv_id":null,"evidence_quote":"It is the source of the $f(R)$+neutrino action, the $\\Gamma$ coupling prescription, and the Hu–Sawicki expressions that define the model."}],"review_version":1}