{"id":"e0be0d22-5096-42ad-bf4f-1e83952ec625","arxiv_id":"2606.02666","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":2.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":1,"one_line_summary":"Exponential f(R) gravity tested on combined CMB, DESI-DR2 and supernova data yields higher S8 and slightly lower H0 than LambdaCDM, alleviating S8 tension by up to 1.2 sigma but not H0 tension.","lead":"The paper constrains parameters of an exponential f(R) gravity model using CMB, DESI-DR2, supernova and other datasets, reporting a moderate reduction in S8 tension but no relief for H0 tension. A generalist might read it to see how one modified-gravity alternative performs against current cosmological data tensions.","discovery_kind":"unclear","skeptic_critique":{"model":"grok-4.3","headline":"Scalar perturbation implementation in Boltzmann solver is the unverified link for the reported S8 shift","rationale":"The reader's weakest_assumption correctly isolates the single assumption whose failure would invalidate the S8 improvement claim. With the full text now available, no other internal inconsistency (e.g., dataset handling or parameter counting) appears more load-bearing than this implementation step.","tokens_in":1758,"tokens_out":292,"duration_ms":14033,"concrete_test":"Extract the linear perturbation equations and initial conditions stated in the methods section; independently recompute the growth factor D(a) and the effective gravitational constant for the best-fit b value using the standard f(R) formulas; if the resulting S8 differs by >0.5σ from the reported posterior, the solver implementation is the source of the shift.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The headline result (higher S8, ~1.2σ tension relief with late-time data) requires that the extra scalar degree of freedom in the exponential f(R) model produces the correct modified growth and slip relations inside the Boltzmann solver. The analysis combines PPS+BBN+CC+DESI-DR2+CMB likelihoods; any mismatch between the analytic perturbation equations (effective G, anisotropic stress, or initial conditions for the scalaron) and the numerical implementation would directly alter the inferred S8 posterior without reflecting the model's actual cosmology.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The paper considers an exponential f(R) gravity model as an extension to GR for late-time acceleration. By combining PPS, BBN, CC, DESI-DR2, and CMB datasets, it constrains the free parameter b and standard cosmological parameters. The analysis finds slightly lower H0 values than in ΛCDM, indicating no significant relaxation of the H0 tension, but systematically higher S8 values leading to a moderate alleviation of the S8 tension by up to ~1.2σ with late-time datasets. The model is concluded to provide a consistent improvement in the description of large-scale structure formation without resolving all tensions simultaneously.","tokens_in":1877,"tokens_out":491,"duration_ms":28844,"significance":"If the scalar perturbation implementation in the Boltzmann solver is accurate, this work demonstrates that exponential f(R) gravity can offer a partial alleviation of the S8 tension through modified growth of structure, while being consistent with the combined datasets. The inclusion of DESI-DR2 data is a positive aspect for the constraints. However, the improvement is obtained by fitting to the same data defining the tension, and the lack of explicit perturbation equations reduces the ability to assess the result independently. The paper does not provide machine-checked proofs or reproducible code.","major_comments":[{"comment":"Abstract: The abstract states the datasets and the direction of the S8 shift but supplies no equations, error budgets, or description of how the f(R) perturbation equations were solved, so the numerical support for the 1.2-sigma claim cannot be verified.","section":"Abstract"},{"comment":"Section describing the perturbation implementation (likely §3 or equivalent): The headline result of higher S8 (~1.2σ tension relief) requires that the extra scalar degree of freedom produces the correct modified growth and slip relations inside the Boltzmann solver. No details are provided on the analytic perturbation equations (effective G, anisotropic stress, or initial conditions for the scalaron), which is load-bearing for the reported S8 posterior.","section":"Section describing the perturbation implementation (likely §3 or equivalent)"}],"minor_comments":[{"comment":"The abstract could be expanded to include the best-fit value and uncertainty on the parameter b for context.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the detailed and constructive report. The comments highlight the need for greater transparency on the perturbation implementation, which we address below by committing to specific revisions. We agree that additional details will strengthen the manuscript's reproducibility and verifiability.","responses":[{"response":"We agree that the abstract is too concise and omits key methodological information. In the revised manuscript we will expand the abstract to include a brief statement on the solution of the f(R) perturbation equations within the Boltzmann solver and the resulting impact on the S8 posterior.","revision_made":"yes","referee_comment":"[Abstract] Abstract: The abstract states the datasets and the direction of the S8 shift but supplies no equations, error budgets, or description of how the f(R) perturbation equations were solved, so the numerical support for the 1.2-sigma claim cannot be verified."},{"response":"The referee correctly identifies that explicit analytic expressions for the perturbation equations are absent from the current text. We will add a new subsection (or expand the existing section on the model) that derives and presents the effective gravitational constant, the anisotropic stress, the slip relation, and the initial conditions for the scalaron, together with their numerical implementation. This will allow independent verification of the modified growth responsible for the reported S8 shift.","revision_made":"yes","referee_comment":"[Section describing the perturbation implementation (likely §3 or equivalent)] Section describing the perturbation implementation (likely §3 or equivalent): The headline result of higher S8 (~1.2σ tension relief) requires that the extra scalar degree of freedom produces the correct modified growth and slip relations inside the Boltzmann solver. No details are provided on the analytic perturbation equations (effective G, anisotropic stress, or initial conditions for the scalaron), which is load-bearing for the reported S8 posterior."}],"tokens_in":1477,"tokens_out":408,"duration_ms":16202,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main takeaway is that this exponential f(R) model, when fit to CMB plus DESI-DR2 and the other listed sets, produces higher S8 values that reduce the tension with late-time structure data by about 1.2 sigma while H0 stays low and offers no relief on that side.\n\nThe incremental step is the use of DESI-DR2 BAO measurements with this particular f(R) form. Earlier papers studied the same exponential function, but the new data combination lets them update the constraints on the single extra parameter b and show the resulting shifts in the standard parameters.\n\nThe paper does a clear job of reporting the direction and size of those shifts and of stating that the model improves the large-scale structure description without fixing both tensions at once.\n\nThe soft spot is the scalar perturbation implementation. The S8 result depends on the Boltzmann solver correctly handling the extra degree of freedom, including the modified growth, slip, and initial conditions for the scalaron. The abstract gives no equations or solver checks, so any mismatch between the analytic expectations for this f(R) and the numerical code would move the posterior without reflecting the actual model. The stress-test note on this point stands; the full text needs to show explicit verification that the numerics match the theory.\n\nThis work is for cosmologists who compile constraints on modified gravity models using current datasets. Someone updating a review of f(R) fits to DESI data would find the numbers useful for comparison.\n\nIt deserves peer review. The dataset combination is timely and the claim is modest, so referees can focus on confirming the perturbation implementation without the paper overreaching.","headline":"Exponential f(R) raises S8 enough to cut tension by 1.2 sigma with DESI-DR2 but leaves H0 untouched and the Boltzmann implementation unverified.","tokens_in":2379,"tokens_out":412,"would_cite":false,"duration_ms":21613,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"An exponential f(R) model raises S8 enough to ease the structure tension by up to 1.2 sigma while leaving the Hubble tension untouched.","keywords":["f(R) gravity","cosmological tensions","S8 tension","H0 tension","large-scale structure","DESI-DR2","CMB constraints"],"falsifier":"Repeating the full MCMC analysis with an independent Boltzmann code that evolves the same f(R) scalar degree of freedom and finding that the S8 posterior no longer shifts upward by the amount reported here would falsify the claimed improvement.","tokens_in":2655,"feed_emoji":"🌌","tokens_out":737,"duration_ms":18467,"temperature":0.7,"pith_summary":"The paper tests whether a specific exponential modification of gravity can ease the mismatch between early-universe and late-universe observations. When the model is fitted to CMB, DESI-DR2, supernova, and other data, the extra scalar degree of freedom produces slightly lower H0 values than Lambda CDM, so the Hubble tension stays roughly the same. At the same time the model systematically increases the predicted S8, bringing the amplitude of matter fluctuations into better agreement with large-scale structure measurements. The net result is a modest improvement in the description of structure formation without simultaneously fixing both tensions.","feed_headline":"Exponential f(R) model lifts S8 and cuts structure tension by 1.2 sigma","feed_subtitle":"CMB plus DESI-DR2 data show higher matter clustering amplitude that better matches observations while H0 stays low.","key_machinery":"The exponential f(R) function whose extra scalar degree of freedom is evolved inside the Boltzmann solver and jointly constrained by the combined likelihood of PPS, BBN, CC, DESI-DR2 and CMB data.","core_discovery":"The exponential f(R) gravity model yields H0 values slightly lower than those obtained in Lambda CDM, indicating no significant relaxation of the H0 tension. In contrast, the model predicts systematically higher values of S8, leading to a moderate alleviation of the S8 tension by up to approximately 1.2 sigma when late-time datasets are included. Overall the results show that the model does not resolve all cosmological tensions simultaneously yet provides a consistent improvement in the description of large-scale structure formation.","pith_inferences":["Future surveys that tighten the S8 measurement from below could quickly rule out or confirm the size of the upward shift reported here.","The same f(R) scalar degree of freedom might be tested independently through its effect on the growth rate f sigma8 at intermediate redshifts.","If the implementation of the scalar mode in the solver is the dominant source of the S8 shift, repeating the analysis with a different code would be a direct check."],"forward_implications":["The model leaves the Hubble tension at its Lambda CDM level or slightly worse.","Inclusion of late-time data produces an S8 value that reduces the tension with weak-lensing surveys by up to 1.2 sigma.","The improvement is tied specifically to the higher S8 prediction rather than to any change in the background expansion history.","The model remains statistically consistent with all datasets used but does not outperform Lambda CDM on every tension at once."],"fun_headline_variants":["Exponential f(R) raises S8 easing tension 1.2 sigma","f(R) model yields higher S8 cutting tension 1.2 sigma","Exponential f(R) gravity increases S8 easing 1.2 sigma tension","Higher S8 from exponential f(R) eases tension 1.2 sigma"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The extra scalar degree of freedom is correctly implemented in the Boltzmann solver and the combined datasets carry no unaccounted systematics that would change the reported S8 shift.","fun_headline_variants_meta":{"raw":{"variants":["Exponential f(R) raises S8 easing tension 1.2 sigma","f(R) model yields higher S8 cutting tension 1.2 sigma","Exponential f(R) gravity increases S8 easing 1.2 sigma tension","Higher S8 from exponential f(R) eases tension 1.2 sigma"]},"model":"grok-4.3","cost_usd":0.00897,"raw_usage":{"total_tokens":4053,"prompt_tokens":715,"num_sources_used":0,"completion_tokens":74,"cost_in_usd_ticks":89699500,"prompt_tokens_details":{"text_tokens":715,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":3264,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":715,"tokens_out":74,"duration_ms":22924,"temperature":1.0,"reasoning_tokens":3264,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-28T13:34:31.711668+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"Repeating the full MCMC analysis with an independent Boltzmann code that evolves the same f(R) scalar degree of freedom and finding that the S8 posterior no longer shifts upward by the amount reported here would falsify the claimed improvement.","supporting_citations":[],"review_version":1}