{"id":"0f41c621-eac4-4a86-9b6d-0ec69a16f70a","arxiv_id":"2412.02707","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Constraints from DESI, cosmic chronometers, gamma-ray bursts, and supernovae show the CEPN modified gravity model fits the cosmic expansion history with Omega_m0 around 0.32 and H0 around 72 to 74 km/s/Mpc, matching LambdaCDM within information criteria in most combinations.","lead":"This paper fits a modified gravity theory called Coupling Extended Proca-Nuevo gravity to recent cosmological datasets, including DESI, supernovae, and gamma-ray bursts. It reports that the model matches the expansion history about as well as the standard cosmological model, with a Hubble constant around 72 to 74 km/s/Mpc.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The CEPN background reduction to Eq. (27) is assumed, not derived; if the non-minimal L2/L3 terms contribute to the vector-field constraint, the fitted H(z) is not the theory's background.","rationale":"The reader's conditional verdict rests on the same deferred derivation, and I agree that this is the weakest load-bearing point. The statistical issues (missing covariances, no convergence diagnostics) are real but they weaken precision, not the existence of a fit; a wrong background equation would make the central compatibility claim vacuous. The concern is not that the cited papers are wrong—only that the manuscript's own argument treats the most consequential step as an imported result without verifying that the arbitrary functions α_n,d_n and non-minimal couplings indeed drop out at the background level. The proposed computation is straightforward and would resolve the issue. Therefore the conditional verdict is appropriate and no change is needed.","tokens_in":18409,"tokens_out":11184,"duration_ms":100540,"concrete_test":"Independently compute the full Euler-Lagrange equation for φ from the action S=∫d^4x√-g [M_Pl^2/2 R + L_EPN + L_M] with L_EPN as in eqs. (11)-(15), on the FLRW background (16)-(17), retaining all terms from L2 and L3. If the resulting equation is not α0,X + 3(α1,X+d1,X)Hφ/Λ^2 = 0 for arbitrary α2,d2,α3,d3, then eq. (20) is special-case, and the paper must either show the cancellation or re-derive eq. (27). A simpler numerical cross-check: substitute the best-fit Ω_m0 and H0 into the full background equations to see whether eq. (27) is actually satisfied.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's observable prediction is the Friedmann equation (27), obtained by eliminating the vector field with the constraint (20). But (20) is stated without derivation from the action (10)-(15). The displayed Lagrangian contains non-minimal terms proportional to R and G^{μν} in L2 and L3 (eqs. 14-15) with arbitrary functions α2,d2,α3,d3; varying those terms with respect to φ generically generates contributions proportional to R, G, and their derivatives, which are not visible in (20). If such contributions do not cancel, the relation ρ_EPN ∝ H^{-2/3} (eq. 26) and hence the fitted H(z) of eq. (27) is not the generic CEPN background. Because the likelihood in Section 4 compares only this H(z), any error or hidden fine-tuning in this reduction invalidates every reported constraint, including Ω_m0=0.317±0.011 and H0=73.89±0.19. The paper cites de Rham et al. 2022a and Anagnostopoulos & Saridakis 2024 but does not reproduce the computation, so the load-bearing step is unverified in this manuscript.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies Covariant Extended Proca-Nuevo (CEPN) gravity, a non-linear spin-1 extension of dRGT massive gravity, and fits its background cosmological solutions to DESI BAO, Cosmic Chronometer (CCh), Gamma-Ray Burst (GRB), and Pantheon+ Type Ia supernova data. The central phenomenological input is the Friedmann equation H^2 = H0^2 Ω_m0 (1+z)^3 + (1-Ω_m0) H0^{8/3} H^{-2/3} (Eq. 27), obtained after eliminating the vector field. The authors report constraints on Ω_m0 and H0 (and r_d for DESI), the present-day deceleration parameter q0, the transition redshift z_t, and information criteria comparing the model with ΛCDM. They conclude that the CEPN model is statistically compatible with the data and that the late-time acceleration is quintessence-like.","tokens_in":18655,"tokens_out":6451,"duration_ms":61194,"significance":"If the background reduction underlying Eq. (27) is valid, the paper provides a simple two-parameter dark-energy model that is competitive with ΛCDM on multiple modern datasets, and the reported H0 ≈ 73.9 km/s/Mpc for the full combination is relevant to the current Hubble tension. The paper has several strengths: it uses the recent DESI DR1 BAO sample, it compares the model against ΛCDM with AIC/BIC/DIC, it gives the explicit fitted Friedmann equation, and it presents derived cosmographic quantities. The fitted parameters are genuinely free and constrained by the data, so the central constraints are not an example of prediction-as-fit circularity. The main limitations are that the phenomenological density (26) is assumed from prior derivations rather than re-derived in the manuscript, and the likelihood ignores covariance information; these issues must be resolved before the observational claims can be accepted.","major_comments":[{"comment":"The vector-field constraint (20) is the load-bearing step, but it is stated without derivation from the action (10)-(15). The Lagrangian contains non-minimal couplings in L2 and L3 (Eqs. (14)-(15)) with arbitrary functions α2(X), d2(X), α3(X), d3(X); varying with respect to φ(t) generically produces terms proportional to R and G_00 (and their derivatives) that are not visible in (20). If such terms do not cancel, Eqs. (23)-(26) do not follow, and the fitted H(z) of Eq. (27) is not the generic CEPN background. The manuscript cites de Rham et al. (2022a) and Anagnostopoulos & Saridakis (2024) but does not reproduce the computation. Because every reported constraint depends on Eq. (27), please provide a self-contained derivation of (20) (or an explicit statement of the exact assumptions under which it holds), or alternatively demonstrate that the non-minimal terms' contributions vanish identically on the FLRW ansatz.","section":"Section 2.2, Eqs. (10)-(15) and (20)"},{"comment":"The likelihood is written as a sum of independent per-point chi-squared terms. This neglects the covariance matrices of the DESI BAO observables (D_M/r_d, D_H/r_d, D_V/r_d in Figure 3) and of the GRB distance moduli, which are correlated quantities. Using only diagonal errors can bias the best-fit values and under- or over-estimate the quoted uncertainties, and it directly affects the ΔAIC/ΔBIC values in Table 3. Please include the full covariance matrices for the BAO and GRB data, or state and justify the assumption that off-diagonal correlations are negligible.","section":"Section 4, Eq. (35)"},{"comment":"The MCMC analysis is not fully reproducible: the manuscript does not report the prior distributions used, the chain length, the burn-in period, acceptance rates, or any convergence diagnostics (e.g., Gelman-Rubin). These details are necessary to assess whether the 1σ intervals in Table 1 and Figure 1 are stable and to allow independent verification. Please add them.","section":"Section 4 and Table 1"}],"minor_comments":[{"comment":"The values of q_CEPN,0 and z_t are quoted without uncertainties, so the claim of 'up to 3σ confidence' is not supported for these derived quantities; please propagate the parameter errors.","section":"Table 2"},{"comment":"In the DESI row, Ω_m0 is reported as 0.298^{+0.034}_{-0.20}; the lower error bar is implausibly large compared with the upper one and is likely a typographical error. Please check.","section":"Table 1"},{"comment":"The legend in the first-column panels reads 'CDM model' where the red curve is the CEPN model; the caption should be corrected to avoid confusion.","section":"Figure 2"},{"comment":"The CCh data are said to be taken from Table 1 of Sudharani et al. (2024a); for reproducibility, please include the actual data points or the table in an appendix, or provide a public link.","section":"Section 3.2"},{"comment":"The abstract and conclusion state that the theory 'is shown to yield reliable, ghost-free cosmological solutions,' but this manuscript only fits the background equations; the ghost-free and stability claims are inherited from prior work and should be flagged as such rather than presented as results of this paper.","section":"Abstract and Section 6"},{"comment":"The histograms lack axis labels with units; please clarify that the horizontal axis is H0 in km/s/Mpc and label the vertical axis as counts or probability density.","section":"Figure 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is essentially an observational application of a background reduction taken from de Rham et al. (2022a) and Anagnostopoulos & Saridakis (2024). The key question for the journal is whether the unverified reduction (Eq. 20 to Eq. 26) is acceptable as a prior result. If the reduction is trusted, the covariance-matrix omission becomes the main statistical concern. I would encourage the editor to ask the authors to either present the constraint derivation in an appendix or explicitly delimit the regime in which Eq. (27) is the correct background; otherwise the reported constraints are conditional on an unexamined assumption. The self-citation for the CCh table is not a problem per se, but the data should be made publicly available."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a competent but incremental application of Coupled Extended Proca-Nuevo gravity to more datasets. The two things you should know: the novelty is real but modest—DESI DR1, GRBs, and a WiggleZ comparison on top of the existing CCh+SNeIa fit in Anagnostopoulos & Saridakis (2024). And the main risk isn't the statistics; it's that the background equation (27) is taken from earlier work without derivation, and the non-minimal couplings in the Lagrangian could, in principle, modify the constraint (20) that makes that equation valid.\n\nWhat the paper does well: the chi-square fits, contours, and information criteria are standard and reproducible in principle. They compare with the earlier CCh+SNeIa result and show consistency. The data selection is reasonable. The conclusion is appropriately cautious about the perturbative level being out of scope.\n\nSoft spots, in rough order of seriousness. First, the constraint (20) and the density (26) are asserted, not derived. The Lagrangian includes terms like alpha2(X) R and alpha3(X) K G, and varying those with respect to phi should generically produce R- and G-dependent terms in the vector equation. The paper just says the result reduces to (20). If that reduction is wrong—or requires delicate cancellations that aren't verified here—then every parameter constraint in Table 1 is modeling a theory that may not be the one in the action. This is exactly the kind of thing a referee should push on. It may be fine; de Rham et al. 2022a presumably did the work. But the authors need to show or cite the specific derivation.\n\nSecond, the statistics are a bit rougher than the tables imply. The DESI BAO and GRB data are treated with diagonal errors; no covariance matrices are used. There are no MCMC priors, chain lengths, or convergence diagnostics reported. q0 and z_t are quoted without uncertainties. None of these are fatal, but they limit how much weight you can put on the reported 1-sigma values.\n\nThird, the H0 tension claim is overstated. The full fit gives H0 = 73.89, but that's driven by the Cepheid-calibrated SNeIa in Pantheon+, not by DESI. The paper admits this later, but the abstract and conclusion frame DESI as the key to alleviating tension. That's misleading.\n\nOverall, the central result—that CEPN fits the combined data about as well as LambdaCDM—is plausible and consistent with the tables. The soft spots are addressable. I'd send it to a serious referee, mainly to check the theoretical reduction and to request proper covariance handling and MCMC details. It's not a breakthrough, but it's a solid data paper if the theory check passes.","headline":"Workmanlike extension of existing CEPN constraints to DESI/GRB data; the key theoretical reduction is cited, not derived, and the statistics are rougher than the tables suggest, but the central fit is plausible.","tokens_in":19212,"tokens_out":4459,"would_cite":false,"duration_ms":39855,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83F05","83D05"],"pacs":["98.80.-k","98.80.Es","04.50.Kd"],"model":"deepseek-v4-flash","headline":"A modified gravity model with a massive spin-1 field fits the latest cosmic expansion data at up to 3-sigma.","keywords":["Coupling Extended Proca-Nuevo gravity","massive spin-1 field","dark energy","Hubble tension","DESI BAO","cosmic chronometer","gamma-ray burst","Pantheon+ supernovae"],"falsifier":"High-precision measurements of the Hubble parameter in the redshift window $0.5 \\lesssim z \\lesssim 1.5$ could settle the claim, because that is where the $H^{-2/3}$ term departs most from $\\Lambda$CDM; if those measurements trace the $\\Lambda$CDM curve within about one percent while the CEPN best fit requires a several-percent deviation, the model's central claim fails. A second decisive check is to compute the sound horizon $r_d$ from the CEPN background at early times and compare it with the DESI-inferred value, since the paper fixes $r_d$ from the fit rather than deriving it from the model itself.","tokens_in":18148,"feed_emoji":"🌌","tokens_out":9863,"duration_ms":72000,"temperature":0.7,"pith_summary":"This paper claims that Covariant Extended Proca-Nuevo (CEPN) gravity, a non-linear theory of a massive spin-1 field extending dRGT massive gravity, can produce the late-time acceleration of the Universe without a cosmological constant. After reducing the cosmological background to a single algebraic Friedmann equation in the Hubble rate, the paper fits that equation to DESI BAO, cosmic chronometer, gamma-ray burst, and Pantheon+ supernova data and reports statistical compatibility at up to 3-$\\sigma$. The full combined fit yields $\\Omega_{m0} = 0.317 \\pm 0.011$, $H_0 = 73.89 \\pm 0.19$ km/s/Mpc, a present deceleration parameter around $-0.575$, and a transition redshift near $0.60$, which the authors read as quintessence-like behavior. The reason to care is that this is a concrete test of whether a theoretically motivated massive-vector-field theory can accommodate the expansion data and speak to the Hubble tension.","feed_headline":"Spin-1 gravity model fits cosmic expansion to 3-sigma","feed_subtitle":"CEPN gravity with DESI, GRB, chronometer, and supernova data yields H0 = 73.9 and quintessence-like acceleration.","key_machinery":"The load-bearing object is the modified Friedmann equation obtained after the vector-field constraint is solved and constants eliminated: $H^2 = H_0^2 \\Omega_{m0}(1+z)^3 + (1-\\Omega_{m0})H_0^{8/3}H^{-2/3}$. This single algebraic equation encodes the entire CEPN dark-energy sector for the background, and it is the expression the MCMC code feeds into the likelihood to produce theoretical Hubble parameters, distance moduli, and BAO distances. Its distinctive feature is the fractional-power term $H^{-2/3}$: it is subdominant at high redshift, reproduces $\\Lambda$CDM-like behavior today, and produces a small but detectable deviation at intermediate redshifts, which is exactly where the data in this paper have constraining power.","core_discovery":"On the paper's own terms, the central discovery is that the CEPN background Friedmann equation, with its distinctive $H^{-2/3}$ dark-energy term, is compatible with the current generation of expansion-rate and distance data. The fit to the DESI+GRBs+CCh+SNeIa combination gives $H_0 = 73.89 \\pm 0.19$ km/s/Mpc and $\\Omega_{m0} = 0.317^{+0.011}_{-0.012}$, with the deceleration parameter $q_0$ near $-0.575$ and the transition from deceleration to acceleration near $z_t \\approx 0.60$. The paper argues that the information-criteria comparison shows the CEPN model is statistically compatible with $\\Lambda$CDM for most dataset combinations, and that the current accelerated expansion is quintessential. It does not claim the model replaces $\\Lambda$CDM; it claims the model is a viable ghost-free alternative whose extra term becomes significant at intermediate redshifts.","pith_inferences":["A natural next step the paper leaves implicit is a joint fit that includes the CMB, since the model's early-time sound horizon is not fixed by the background fit; such a fit would either validate or expose the $H^{-2/3}$ term at recombination.","The fractional-power $H^{-2/3}$ dependence is an unusual prediction: it implies the dark-energy density grows relative to matter as the Universe expands into the far future, unlike a cosmological constant, a difference that upcoming deep surveys could in principle observe.","Because the paper fixes $r_d$ rather than computing it, the quoted $H_0$ constraint inherits the DESI sound-horizon calibration; a self-consistent computation would propagate that uncertainty into $H_0$."],"forward_implications":["The CEPN model provides a background-level dark-energy extension with no extra free parameters beyond $\\Omega_{m0}$ and $H_0$, so future expansion-rate surveys can directly test its distinctive $H^{-2/3}$ signature.","DESI BAO data alone put a tighter constraint on $H_0$ than the older WiggleZ sample, and the paper argues this stronger constraining power can help in addressing the Hubble tension.","When Pantheon+ Cepheid distances are included, the combined fit returns $H_0$ on the higher side, $73.89$ km/s/Mpc, so the model does not eliminate the discrepancy with the lower Planck value; it shifts toward the local measurement.","The information-criteria analysis shows the CEPN model and $\\Lambda$CDM are statistically indistinguishable for the DESI, DESI+GRBs, and DESI+GRBs+CCh combinations, meaning the CEPN term is a viable alternative rather than a disfavoured one.","The derived $q_0$ and $z_t$ values place the model in the quintessence camp, which is a concrete, falsifiable prediction about the equation of state of the dark-energy sector."],"supporting_citations":[{"why":"supplies the covariant CEPN Lagrangian and the background reduction that yields the effective dark-energy density.","marker":"de Rham et al. 2022a"},{"why":"provides the Friedmann equation with the $H^{-2/3}$ term and the CCh+SNeIa comparison that the paper regenerates.","marker":"Anagnostopoulos & Saridakis 2024"},{"why":"introduces Proca-Nuevo theory and its null-eigenvector constraint structure.","marker":"de Rham & Pozsgay 2020"},{"why":"the DESI BAO survey data used in the fit.","marker":"Adame et al. 2024a"},{"why":"the compiled table of DESI BAO measurements at the redshifts used.","marker":"Pourojaghi et al. 2024"},{"why":"the Pantheon+ supernova sample with Cepheid distances.","marker":"Scolnic et al. 2022"},{"why":"the cosmic chronometer Hubble-parameter measurements and their error treatment.","marker":"Moresco et al. 2020"},{"why":"the 162-long-GRB sample and its redshift distribution.","marker":"Demianski et al. 2017"},{"why":"the emcee MCMC sampler used for parameter estimation.","marker":"Foreman-Mackey et al. 2013"}],"fun_headline_variants":["Ghost-free spin-1 gravity fits expansion data to 3σ","CEPN gravity fits DESI, GRB, and SNe data","New gravity model yields H0=73.9, quintessence-like expansion","Spin-1 gravity passes cosmic expansion tests"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire fit relies on the prior derivation that the full CEPN Lagrangian, with its arbitrary functions, collapses to the simple $H^{-2/3}$ dark-energy density, a reduction the paper takes from earlier work and does not re-derive, so a failure or heavy fine-tuning in that reduction would invalidate all the reported constraints.","fun_headline_variants_meta":{"raw":{"variants":["Ghost-free spin-1 gravity fits expansion data to 3σ","CEPN gravity fits DESI, GRB, and SNe data","New gravity model yields H0=73.9, quintessence-like expansion","Spin-1 gravity passes cosmic expansion tests"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000901,"raw_usage":{"total_tokens":3867,"prompt_tokens":921,"completion_tokens":2946,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":537,"completion_tokens_details":{"reasoning_tokens":2872}},"tokens_in":537,"tokens_out":2946,"duration_ms":22904,"temperature":1.0,"reasoning_tokens":2872,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:16:48.932258+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"High-precision measurements of the Hubble parameter in the redshift window $0.5 \\lesssim z \\lesssim 1.5$ could settle the claim, because that is where the $H^{-2/3}$ term departs most from $\\Lambda$CDM; if those measurements trace the $\\Lambda$CDM curve within about one percent while the CEPN best fit requires a several-percent deviation, the model's central claim fails. A second decisive check is to compute the sound horizon $r_d$ from the CEPN background at early times and compare it with the DESI-inferred value, since the paper fixes $r_d$ from the fit rather than deriving it from the model itself.","supporting_citations":[{"cited_title":"K., & Saridakis, E","cited_arxiv_id":null,"evidence_quote":"provides the Friedmann equation with the $H^{-2/3}$ term and the CCh+SNeIa comparison that the paper regenerates."},{"cited_title":"2020, Phys","cited_arxiv_id":null,"evidence_quote":"introduces Proca-Nuevo theory and its null-eigenvector constraint structure."}],"review_version":1}