{"id":"5dad19c2-2df7-4960-9e01-c361fecc2be9","arxiv_id":"2607.25841","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"First full Boltzmann-code implementation of entropic dark energy, with MCMC constraints from CMB+BAO+SN, yields α≈1 and a fit statistically indistinguishable from ΛCDM.","lead":"Cosmologists built the first full simulation-code implementation of an \"entropic acceleration\" model, in which cosmic acceleration comes from horizon entropy gain rather than a cosmological constant, and fitted it to CMB, galaxy, and supernova data. It fits about as well as the standard model, but the key parameter is measured, not independently predicted.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim is conditional on an un-derived perturbation sector: the 'full CMB spectra' are those of a smooth effective fluid, and the c_s² scan does not bracket non-adiabatic or non-local δfμν effects.","rationale":"The reader's weakest assumption identifies the same load-bearing concern: the un-derived perturbed entropic source δfμν is approximated by a smooth effective fluid with c_s²=1 and PPF closure. This is indeed the point on which the central claim—that the paper provides the full CMB and matter power spectra of GREA—is least secure. The code itself appears transparent and the background validation, including the fσ8 reproduction and normalization test, gives genuine support to the implementation. But the perturbation sector is not derived from GREA, and the c_s² scan does not cover the non-adiabatic or non-local effects that δfμν could introduce. The paper is honest about this limitation, yet the abstract and conclusions overstate the status of the resulting spectra and constraints. Since the reader already recommends CONDITIONAL and my concern supports that assessment, no change to the verdict is needed. The concrete derivation-and-reanalysis test would settle whether the approximation is harmless or whether the constraints shift.","tokens_in":29915,"tokens_out":7850,"duration_ms":75267,"concrete_test":"Derive δfμν to linear order from the covariant bulk-viscosity formalism (Gagnon & Lesgourgues 2011), computing δd_H from the perturbed light cone and the resulting non-adiabatic pressure. Implement this exact source in the public CLASS_GREA code and re-run the CMB-SPA+DESI+PP and CMB-SPA+DESDovekie+DESI chains. If the inferred α or the low-ℓ TT/lensing spectra shift by more than the quoted 1σ uncertainties—or by more than the <1.6% bracket from Appendix B—the effective-fluid approximation is not adequate for the headline constraints.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III explicitly leaves the perturbed entropic source undetermined: 'Whether this non-local piece can be neglected (δζ=0) is the central open question of a complete derivation. We do not settle it here.' The implementation then replaces δfμν with a local imperfect fluid characterized by w(a), c_s²=1, and PPF closure. This is an assumption, not a consequence of GREA. Appendix B's robustness scan only varies c_s² over [0.01,1], which brackets the local clustering-to-smooth transition, but it does not sample the non-adiabatic pressure or the non-local light-cone response that a genuine δfμν could generate. The paper itself identifies the low-ℓ ISW and lensing amplitude as the channels where δfμν matters. Therefore the claimed full CMB, lensing, and matter power spectra—and the CMB-SPA constraints on α—are not uniquely determined by GREA; they are determined only by an unvalidated effective-fluid approximation. The conclusion that α≈1 is in 'excellent agreement with the theoretical prediction' additionally conflates a fitted parameter with an independent prediction, but the more load-bearing technical issue is the perturbation-sector assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports the first implementation of the General Relativistic Entropic Acceleration (GREA) model inside the CLASS Einstein–Boltzmann solver. The GREA background is integrated directly, while the entropic component is treated at linear order as an effective dark-energy fluid with the GREA equation of state, PPF closure, and sound speed c_s^2=1. The code is validated against semi-analytic growth and used to compute CMB TT/TE/EE, lensing, lensing-induced BB, matter power spectrum, and fσ8. The authors run MCMC analyses with COBAYA against CMB-SPA, DESI DR2 BAO, and several SN samples, report α≈1 in most combinations with |Δχ²|≲6 relative to ΛCDM, and interpret the results as confirming the GREA prediction α∼1 and a second phantom crossing at z≈2.","tokens_in":30268,"tokens_out":4855,"duration_ms":46392,"significance":"If the effective-fluid treatment is accepted as a faithful proxy for the GREA perturbation sector, this is a useful and reproducible technical step: it opens full primary-CMB and lensing likelihoods to a one-parameter, thermodynamically motivated alternative to ΛCDM. The public code, validation harness, and MCMC setup are concrete strengths, and the background-normalization subtlety is handled carefully. However, the significance is conditional, because the perturbed entropic source δfμν is not derived and the low-ℓ ISW and lensing predictions are not uniquely determined by GREA at this stage.","major_comments":[{"comment":"The paper explicitly states that the perturbed entropic sector δfμν is not derived: 'Whether this non-local piece can be neglected (δζ=0) is the central open question of a complete derivation. We do not settle it here.' The implementation then replaces δfμν by a local effective fluid with w(a), PPF closure, and c_s^2=1. This is an assumption, not a consequence of GREA. Consequently the abstract's claim of 'full CMB and matter power spectra' and the later statement that spectra are computed 'self-consistently' overstate what is uniquely predicted. The spectra are those of the effective-fluid ansatz. This is a load-bearing issue because the low-ℓ ISW and lensing channels are exactly the parts of the data that depend on this assumption; the Table II constraints are therefore conditional on the ansatz. I recommend explicitly reframing the title/abstract/conclusions as an implementation of th","section":"Section III (paragraph starting 'A fully consistent Einstein–Boltzmann treatment')"},{"comment":"The robustness scan for c_s^2∈[0.01,1] varies only the local clustering-to-smooth transition of a single adiabatic sound speed. It does not sample non-adiabatic pressure or the non-local light-cone response that a genuine δfμν could generate. The sentence in Section III that 'the plausible size of the effect is bracketed by the c_s^2 sensitivity analysis' is therefore too strong: the scan brackets only one family of local closures. A non-adiabatic or non-local perturbation could shift the low-ℓ ISW and lensing power by a different amount. The paper itself acknowledges this residual caveat, but the claimed 1.6% robustness should be stated as applying only to the local adiabatic sound-speed ambiguity.","section":"Appendix B and Section III (c_s^2 sensitivity)"},{"comment":"The conclusion that α≈1 is 'in excellent agreement with the theoretical prediction' conflates parameter estimation with confirmation of a prediction. The quantity √(−k)η0 is sampled with a flat prior and later mapped to α via Eq. (8); the posterior is a fitted parameter, not an independent test of a prediction unless the theory supplies a separate, prior expectation with uncertainty. The data being consistent with α=1 at the 1σ level is a meaningful result, but it should be phrased as consistency, not as the data 'singling out' a predicted value. Similarly, the second phantom crossing at z≈2 is beyond the last bin of the reconstruction shown in Fig. 9 and is evaluated at the best fit, so it is a model prediction not yet probed by the data; this should be stated more carefully.","section":"Eq. (8), Table II, Section VI"}],"minor_comments":[{"comment":"Typographical: 'aclassimplementation' should be 'a CLASS implementation' with spaces. Also, 'CMB-SP A' appears with an irregular space in several places in Section V; use 'CMB-SPA' consistently.","section":"Title / Abstract"},{"comment":"The caption notes that 'most of these configurations have two crossings of the phantom divide,' while the text elsewhere refers to 'a transient phantom crossing' in the singular. Please clarify how many crossings occur in the baseline α≈1 case and which crossing is being discussed.","section":"Figure 2 caption"},{"comment":"The sentence 'the largest is|lnB| ≃2.72' is grammatically incomplete; specify for which dataset combination this value occurs. It appears to be the CMB-SPA+PP+DESI column in Table II, but should be stated explicitly.","section":"Section V, near Eq. (13)"},{"comment":"When introducing the sampled variable √(−k)η0 = α D_H(0), it may help readers to state the implied prior range for α, since the flat prior [2.5,4.5] on √(−k)η0 translates into a non-uniform prior on α through the background mapping. This would also clarify the 'prediction' discussion in the conclusions.","section":"Section IV"},{"comment":"The top panel caption says the shaded band indicates the 68% confidence region of the GREA reconstruction, but the plotted curves are said to be evaluated at the best fit; please specify whether the band is from the posterior or from the parameter covariance.","section":"Figure 10 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper is honest about its main limitation, which I appreciate. My recommendation is driven by the gap between the abstract's 'full CMB and matter power spectra' language and the fact that the perturbation sector is an un-derived effective-fluid placeholder. This can be fixed by careful reframing and by softening the 'prediction confirmed' language. The technical implementation and MCMC pipeline appear solid, and the public release is a strength. I would not reject the paper, but the current central claims go beyond what the model as implemented actually predicts."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I've read D'Onofrio et al. on GREA in CLASS. The honest take: this is a technically competent, transparent implementation of a one-parameter dark-energy alternative, and the first to push GREA through a full Einstein-Boltzmann likelihood analysis. The code is public, the validation harness is sensible (reproducing the semi-analytic fσ8 is a good check), and the c_s^2 robustness scan is a genuine plus. The background normalization subtlety — H0 in the definition of τ is not the physical present-day rate — is handled correctly and explained clearly.\n\nWhat's actually new: full CMB temperature/polarization/lensing spectra for GREA, MCMC against CMB-SPA+DESI+SN, and constraints on α from full primary CMB rather than compressed distance priors. That is a real step past the previous background-only analyses [135,136]. The paper also flags its central limitation candidly: the perturbed entropic source δfμν is not derived, and the linear sector is modeled as an effective fluid with PPF closure and c_s^2=1. Appendix B shows that varying c_s^2 between 0.01 and 1 shifts all observables by less than 1.6%, so the main constraints on α are probably robust to that choice. Good.\n\nThe soft spots, in order of importance.\n\nFirst, the α≈1 story. α is not an independent prediction of the theory; it is a defined ratio, and the sampled parameter is sqrt(-k)η0 with a flat prior. Recovering α≈1 from the posterior and calling it 'excellent agreement with the theoretical prediction' is circular. The abstract and conclusions lean on that framing heavily. The authors should either derive a genuine prior or prediction for α, or simply report the constraints without the vindication language. This is a rhetorical issue, not a mathematical one, but it matters for how readers interpret the headline result.\n\nSecond, the perturbation sector. The 'full CMB spectra of GREA' are spectra of an effective fluid, not of the complete theory. The non-adiabatic pressure and non-local light-cone response of δfμν could, in principle, alter the low-ℓ ISW and lensing amplitude. The c_s^2 scan does not bracket those effects. The paper acknowledges this explicitly, and the impact on the main constraints is probably small, but the claim that these are 'self-consistent' GREA predictions is a stretch. Reframing the sections to say 'spectra of the effective-fluid approximation' would be more accurate.\n\nThird, the second phantom crossing at z≈2 is a model prediction but lies beyond the current data bins; the authors admit this, so it's not a fatal flaw, but it's not a test yet.\n\nThe paper deserves a serious referee. It is not a desk reject. The code and data are public, the implementation is careful, and the model is interesting enough to warrant scrutiny. My recommendation: send it to review, but require the authors to reframe the α prediction claim and to be explicit that the perturbation-level outputs are provisional pending a derivation of δfμν.\n\nI'd bring it to a reading group if we want to dissect 'prediction vs posterior' language. I'd cite it for the CLASS implementation and constraints, though not for the theory itself.","headline":"Solid first Boltzmann implementation of GREA, but the α≈1 'prediction' is a fitted parameter and the perturbation sector is an effective-fluid placeholder — still worth refereeing.","tokens_in":30775,"tokens_out":4058,"would_cite":true,"duration_ms":46462,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83F05","83C55","83-08"],"pacs":["98.80.-k","95.36.+x","04.20.-q"],"model":"deepseek-v4-flash","headline":"This paper reports the first Boltzmann-solver implementation of GREA, showing that the model's single parameter α is constrained to ≈1 by CMB, BAO, and supernova data, matching ΛCDM's fit within |Δχ²|≲6.","keywords":["General Relativistic Entropic Acceleration","horizon thermodynamics","dark energy","cosmological constant problem","Einstein-Boltzmann solver","CMB power spectra","parametrized post-Friedmann","phantom crossing"],"falsifier":"Measure the growth index γ(z) from redshift-space distortions (e.g., DESI, Euclid): GREA predicts a rising γ(z) with dγ/dz>0, while ΛCDM and most quintessence models predict a decreasing γ(z). If future data show dγ/dz<0, the model's growth prediction is falsified; alternatively, a detection of non-adiabatic pressure in the dark-energy fluid would invalidate the δζ=0 closure.","tokens_in":29764,"feed_emoji":"🌌","tokens_out":5328,"duration_ms":54624,"temperature":0.7,"pith_summary":"The paper pushes General Relativistic Entropic Acceleration (GREA) from a background-level idea into full perturbation theory, implementing it inside an Einstein–Boltzmann solver to compute CMB, lensing, and matter power spectra for the first time. It then confronts the model with the full primary-CMB, DESI BAO, and Type Ia supernova data via Markov-chain Monte Carlo. The central result is that GREA's single coupling α is measured to be about 1, exactly where the thermodynamic derivation places it, and that the one-parameter model fits the data as well as ΛCDM, with |Δχ²|≤6 across nine dataset combinations. If correct, this shows that a thermodynamically motivated, parameter-light alternative to the cosmological constant can survive the same precision tests as Λ, making acceleration a consequence of horizon entropy rather than vacuum energy.","feed_headline":"α ≈ 1: entropic-acceleration model matches ΛCDM on full data","feed_subtitle":"First Boltzmann-level test of the horizon-entropy model: its one parameter is pinned to the predicted value, with |Δχ²| ≲ 6.","key_machinery":"The key object is the entropic-force tensor f_μν and its homogeneous limit: an effective dark-energy density ρ_GREA ∝ sinh(2τ)/a², where τ is the dimensionless conformal time tied to the causal horizon. Its dynamics are fixed by a single parameter α, defined by α D_H(z=0) = √(-k) η0. At the perturbation level, the machinery is the effective-fluid description: the component is assigned the GREA equation of state w(a) and evolved as a non-clustering fluid (c_s²=1) under the parametrized-post-Friedmann closure, which handles the w=-1 crossing. The implementation also requires a subtle normalization rescaling so that the physical present-day Hubble rate, not the fiducial H0, is what the Boltzman","core_discovery":"The authors integrate the GREA background directly into a Boltzmann solver and evolve the entropic component as an effective dark-energy fluid with sound speed c_s² = 1, regulated by the parametrized-post-Friedmann scheme so perturbations remain regular through the phantom crossing. The resulting angular power spectra and growth functions are new. When fit to the CMB-SPA (Planck+ACT+SPT) likelihood, DESI DR2 BAO, and Pantheon+/DES Dovekie supernovae, the inferred coupling α—the ratio of spatial-curvature scale to the causal horizon today—clusters tightly around unity (α≈1.00–1.08 for the CMB-anchored combinations), in agreement with the model's parameter-free prediction. The fit is statistic","pith_inferences":["A full first-principles derivation of the perturbed entropic source δf_μν might reveal non-adiabatic pressure or non-local horizon terms; if those alter the low-ℓ ISW response by more than the ~1.6% sound-speed bracket, the reported constraints could shift, though likely within current error bars.","The model's thermodynamic origin hints that black-hole entropy growth (already discussed by the authors) could produce analogous acceleration signatures; testing GREA with gravitational-wave standard sirens could further distinguish it from scalar-field dark energy.","Since GREA is not nested in ΛCDM, the Δχ²≈0 is more telling than a Bayesian preference; a future comparison using Bayesian evidence with theory-motivated priors (fine-tuning penalty for Λ) would likely shift the model odds toward GREA."],"forward_implications":["If the central claim holds, GREA provides a single-parameter, thermodynamically grounded rival to ΛCDM that reproduces the expansion history, CMB anisotropies, and growth data at the same quality, eliminating the need for a finely tuned cosmological constant.","The distinctive second phantom crossing at z≈2 is a clean, falsifiable prediction that future higher-redshift dark-energy reconstructions (e.g., from DESI Lyman-α and Euclid) can test.","The model predicts enhanced growth at fixed primordial amplitude, raising σ8 relative to ΛCDM; a dedicated weak-lensing analysis within the GREA framework would determine whether this worsens or reframes the S8 tension.","The full Boltzmann-level implementation opens the way to constrain GREA with lensing, ISW, and growth-rate data from upcoming surveys, where the model's signatures are expected to be decisive.","The sound-speed sensitivity analysis (all observables shift by <1.6%) indicates current constraints are robust to the one undetermined ingredient, so the α≈1 measurement is stable under the effective-fluid assumptions."],"fun_headline_variants":["Entropic acceleration passes first Boltzmann test, α≈1","Horizon entropy predicts α≈1 and matches ΛCDM on all data","No cosmological constant: entropy growth accelerates universe, α≈1","CMB and BAO pin entropic model to α≈1, matching ΛCDM","Phantom crossing at z≈2 emerges from entropic gravity fit"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The perturbed entropic source δf_μν is not derived from first principles; the paper sets its perturbation to zero (δζ=0) and models the component as a smooth fluid with sound speed c_s²=1, so if the true perturbations carry non-adiabatic or non-local pieces, the computed ISW and lensing spectra—and therefore the reported constraints—are not uniquely determined by GREA alone.","fun_headline_variants_meta":{"raw":{"variants":["Entropic acceleration passes first Boltzmann test, α≈1","Horizon entropy predicts α≈1 and matches ΛCDM on all data","No cosmological constant: entropy growth accelerates universe, α≈1","CMB and BAO pin entropic model to α≈1, matching ΛCDM","Phantom crossing at z≈2 emerges from entropic gravity fit"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000204,"raw_usage":{"total_tokens":1254,"prompt_tokens":799,"completion_tokens":455,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":543,"completion_tokens_details":{"reasoning_tokens":360}},"tokens_in":543,"tokens_out":455,"duration_ms":5050,"temperature":1.0,"reasoning_tokens":360,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T01:19:45.631667+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the growth index γ(z) from redshift-space distortions (e.g., DESI, Euclid): GREA predicts a rising γ(z) with dγ/dz>0, while ΛCDM and most quintessence models predict a decreasing γ(z). If future data show dγ/dz<0, the model's growth prediction is falsified; alternatively, a detection of non-adiabatic pressure in the dark-energy fluid would invalidate the δζ=0 closure.","supporting_citations":[],"review_version":1}