{"id":"de1b0482-cf42-4989-a7d2-b3bc8ba6051a","arxiv_id":"2506.06388","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":6,"one_line_summary":"A log-periodic deceleration parameter in f(R,G)=R+alpha R^2+beta e^{gamma G} gravity, fit to CC, BAO and Pantheon+SHOES data, gives H0 around 71.7 to 72.8 km/s/Mpc and q0 around -0.5, but the f(R,G) parameters are chosen, not constrained.","lead":"This paper fits a new wiggly formula for how fast the universe's expansion is changing, q(z)=q0+q1 sin[log(1+z)], inside a modified gravity model, using supernova, galaxy, and clock data. It reports a Hubble constant near 72 km/s/Mpc and present acceleration, but the gravity model's extra coefficients are chosen by hand, not measured from the data.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The physical conclusions are not consequences of the fitted model: α,β,γ are fixed by hand in Eq. (10), and the entropy derivation assumes 0<ω<1 while the fitted ω0≈−0.49, so the claim of a viable f(R,G) description is unsupported.","rationale":"The reader's weakest-assumption diagnosis is correct and the strongest load-bearing concern is exactly that the f(R,G) parameters are unconstrained. The H(z) fit itself is a reasonable kinematic reconstruction with three parameters, and the reported H0 and q0 values are broadly consistent with late-time constraints. But the paper's abstract and conclusion go beyond the kinematic fit: they assert that the f(R,G) model, with specific α, β, γ values, gives viable density, pressure, equation of state, energy conditions, and thermodynamics. Since the fitted data never touch α, β, or γ, the derived curves are not validated by the MCMC analysis. An independent check would be to vary α, β, and γ and see whether the qualitative claims survive; if they do not, the central physical conclusions collapse to a particular choice of hand-picked parameters. The entropy section adds a separable internal inconsistency: the derivation of T and S explicitly assumes 0<ω<1, whereas the fitted model has ω0≈−0.49 and negative ω over a broad redshift range. This means the thermodynamic 'consistency' claim is not supported even internally. Both issues point in the same direction: the paper provides a plausible kinematic description, but it does not demonstrate that the specific f(R,G) gravity model is observationally viable. Therefore the reader's REJECT verdict is unchanged.","tokens_in":17283,"tokens_out":2904,"duration_ms":38458,"concrete_test":"Take the best-fit H0, q0, q1 from the CC+BAO+Pantheon+SHOES row and rescan a coarse log-spaced grid over α, β, γ ∈ [10^-2, 10], recomputing Eq. (25), the NEC/DEC/SEC expressions, and Eq. (44). If ω0 or the sign of any energy condition changes across the grid, the plotted physical behavior is an artifact of the chosen (0.5, 0.6, 0.7), not a robust model prediction.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The kinematic fit in Eq. (13) is a genuine constraint on H0, q0, and q1 and is independent of α, β, and γ. However, all of the paper's derived physical quantities—ρ(z) in Eq. (23), p(z) in Eq. (24), ω(z) in Eq. (25), the energy conditions in Eqs. (26)-(28), and the thermodynamic quantities in Eqs. (44)-(45)—depend on α=0.5, β=0.6, γ=0.7, which are selected by hand in the Figure 4 caption and never constrained by the MCMC analysis. Different admissible values change ρ, p, and ω and can flip the signs of the energy conditions, so those plots are not predictions of the fitted model but illustrations of one arbitrary parameter slice. In addition, the thermodynamic derivation in Section 8 assumes a barotropic fluid with 0<ω<1 in Eqs. (42)-(43) to obtain T ∝ ρ^(ω/(1+ω)). The paper's own fits give ω0≈−0.49 and ω(z) negative over much of the plotted range, so Eqs. (44)-(45) are not applicable to the model's own solutions without further justification. Thus the central claim that this f(R,G) model offers a viable and observationally consistent description of cosmic acceleration rests on untested parameter choices and an internally inconsistent thermodynamic assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies late-time cosmology in the modified gravity theory f(R,G)=R+αR²+βe^{γG}, using a log-periodic deceleration parameter q(z)=q0+q1 sin[log(1+z)]. The Hubble parameter H(z) is derived analytically in Eq. (13), and the parameters H0, q0, q1 are constrained through MCMC using CC, BAO, and Pantheon+SHOES data. The best-fit parameters are then used to reconstruct the energy density, pressure, equation-of-state parameter, energy conditions, statefinder parameters, age of the Universe, and thermodynamic quantities. The authors claim the model offers a viable and observationally consistent description of cosmic acceleration, reporting H0≈71.7–72.8 km/s/Mpc, q0≈−0.48 to −0.52, a transition redshift from 0.879 to 0.744, and ω0≈−0.49.","tokens_in":17606,"tokens_out":5519,"duration_ms":64211,"significance":"The analytic reconstruction from q(z) to H(z) in Eq. (13) is clean and standard, and the use of multiple contemporary datasets is appropriate in principle. If the derived quantities were genuine consequences of a constrained f(R,G) model, the paper would be a useful addition to the modified-gravity phenomenology literature. However, the study's central claim is not supported: the f(R,G) coefficients α, β, γ are fixed by hand, the chi-square analyses ignore important covariance information, and the thermodynamic derivation relies on an assumption contradicted by the model's own fits. These issues are load-bearing for the paper's physical conclusions.","major_comments":[{"comment":"The f(R,G) coefficients α=0.5, β=0.6, γ=0.7 are fixed by hand, as stated in the Figure 4 caption, and are never constrained by the MCMC analysis. Consequently, all derived quantities in Sections 4–8—ρ(z) in Eq. (23), p(z) in Eq. (24), ω(z) in Eq. (25), the energy conditions in Eqs. (26)–(28), and the thermodynamic quantities in Eqs. (44)–(45)—depend on an untested and unjustified choice. Different admissible values of α, β, γ can change the sign of the energy conditions and alter ω(z) entirely. The paper therefore does not test the f(R,G) model; it illustrates one arbitrary parameter slice, so the claim that this model offers a viable and observationally consistent description is unsupported.","section":"§4.2, Eq. (10), Fig. 4"},{"comment":"The chi-square functions for BAO and Pantheon+SHOES use only diagonal uncertainties and ignore the covariance matrices. The DESI BAO measurements and the Pantheon+SHOES supernova sample have significant off-diagonal correlations that are routinely included in the published likelihoods. Omitting these covariances can bias the best-fit parameters and, more importantly, underestimate the uncertainties. Since the paper reports H0 with a relative uncertainty of about 1.7% and makes comparisons with other H0 measurements, this statistical treatment is not reliable and the reported parameter ranges cannot be taken at face value.","section":"§3.1.2–3.1.3, Eqs. (18) and (21)"},{"comment":"The thermodynamic derivation assumes a barotropic fluid with 0<ω<1 to obtain T ∝ ρ^{ω/(1+ω)} and S ∝ ρ^{1/(1+ω)}. However, the model's own fits give ω0≈−0.49 and ω(z) negative over much of the plotted redshift range. Thus Eqs. (44)–(45) are not applicable to the model's solutions without further justification. The temperature and entropy curves in Figure 9 and the conclusion that the generalized second law is satisfied are therefore not consequences of the fitted model; they rest on an internally inconsistent assumption.","section":"§8, Eqs. (42)–(44)"},{"comment":"The expression for the statefinder parameter s in Eq. (32) is inconsistent with its definition in Eq. (30). From Eq. (30), s=(r−1)/[3(q−1/2)], with r given by Eq. (31) as r=2q²+q+q1 cos[log(1+z)]. Substituting gives s=(2q²+q+q1 cos−1)/[3(q−1/2)], but Eq. (32) omits the '−1' in the numerator. The reported values {r0,s0}=(0.866,0.046) and the trajectory in Figure 7 are therefore computed with an incorrect formula and need to be re-evaluated.","section":"§6, Eqs. (31)–(32)"},{"comment":"The paper does not propagate the MCMC uncertainties to derived quantities such as q0, z_tr, ω0, r0, s0, and the age t0. Values like the transition redshift shift from 0.879 (CC) to 0.744 (full data) and the reported age range 13.01–13.59 Gyr are presented without error bars, so it is impossible to assess whether these differences are statistically significant. The plots in Figures 3–9 show only best-fit curves, not confidence regions.","section":"§§4–7"}],"minor_comments":[{"comment":"The dataset name is spelled inconsistently: 'Pantheon+SHOES' in some places and 'Pantheon+SH0ES' in others (e.g., Figures 3–9). Please standardise.","section":"§3.1.3"},{"comment":"The abstract states q0=−0.484 to −0.517 while the text reports q0=−0.484 (CC), −0.495 (CC+BAO), −0.517 (full data); please clarify whether these are best-fit values or 1σ ranges, and report the corresponding uncertainties consistently.","section":"§4.1"},{"comment":"The expressions for ρ(z) and p(z) in Eqs. (23)–(24) are extremely long and a typo in the text ('presuure' in the Figure 4 caption) suggests careful proofreading is needed; a symbolic derivation file or a more compact form would help the reader verify these equations.","section":"§4.2"},{"comment":"In Section 3.1.1, the sentence about the inability to constrain q(z) directly is unclear: 'due to the limited availability and high degree of uncertainty in observational data for q(z), we are unable to properly constrain the deceleration parameter forms.' Since the paper does constrain q0 and q1 through H(z) data, please rephrase to explain what is meant.","section":"§3.1.1"},{"comment":"The relation in Eq. (40) defines entropy as s=(ρ+p)V/T; the paper later uses S=(ρ+p)/T as entropy density. Please define the volume and clarify the distinction between total entropy and entropy density in the equations.","section":"§8"}],"recommendation":"reject","confidential_remarks":"The kinematic reconstruction from q(z) to H(z) is a legitimate exercise, but the paper's physical conclusions about f(R,G) gravity are not supported because α, β, γ are arbitrary, the statistical treatment ignores known covariance matrices, and the thermodynamic section assumes 0<ω<1 while the fits give negative ω. The statefinder formula error adds to the list of load-bearing problems. Even with substantial revision, the central claim would need a fundamentally different analysis in which α, β, γ are constrained or otherwise rigorously justified. As it stands, the paper does not meet the standards for publication in this journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The kinematic core of this paper is more solid than the abstract promises. The choice q(z)=q0+q1 sin[log(1+z)] is a clean, well-behaved parametrization, the integral to H(z) is straightforward, and the MCMC constraints on H0, q0, q1 from CC, BAO and Pantheon+SHOES are sensible and broadly consistent with late-time cosmology. That is a real, modest contribution to the kinematic reconstruction literature, and as far as I can tell the specific functional form is new.\n\nThe problem is that the paper sells itself as a test of f(R,G) gravity, and that part does not hold up. The three parameters α, β, γ in Eq. (10) are set by hand (α=0.5, β=0.6, γ=0.7) and are never constrained. Every derived quantity—ρ(z), p(z), the equation of state, the energy conditions, the temperature and entropy plots—depends on that arbitrary choice. With different α, β, γ those curves change and the energy conditions can flip sign. So they are not predictions of the fitted model; they are illustrations of one parameter slice. The paper cannot claim that this f(R,G) model is observationally viable on that basis.\n\nThe thermodynamics section is worse than just unmotivated. Section 8 assumes a barotropic fluid with 0<ω<1 to derive T ∝ ρ^{ω/(1+ω)} and S ∝ ρ^{1/(1+ω)}. The paper's own fits give ω0 around −0.49, with negative ω over most of the redshift range. The assumption contradicts the model's own solutions, so Eqs. (44)–(45) are inapplicable and the claim that the generalized second law holds is unsupported. That is an internal inconsistency with the paper's own equations, not a debatable modeling choice.\n\nThere are also smaller issues. The χ² for BAO and Pantheon+SHOES ignores covariance matrices, which likely understates the errors and makes the reported 1.7% uncertainty on H0 look overconfident. And the 'predictions' for statefinder, age, and transition redshift are just analytic restatements of the assumed q(z) form with fitted parameters; they are consistency checks of the parametrization, not independent successes.\n\nWhat the paper does well: it is clearly written around the kinematic fit, cites relevant work, and the log-periodic form is simple enough that the reconstruction could be reused. If the authors reframed the paper as a kinematic reconstruction, removed or constrained the f(R,G) parameters, and rewrote the thermodynamics with a consistent ω, I would take it seriously. As submitted, the central claim is not supported.\n\nWho should read it: people interested in new q(z) parametrizations and late-time kinematic fits. It is not a reliable constraint on modified gravity. I would not desk-reject it outright; the kinematic part deserves referee time, but the referee should demand major revision.","headline":"A clean new q(z) parametrization buried under an unsupported f(R,G) claim and a thermodynamics section that contradicts the paper's own equation of state.","tokens_in":18150,"tokens_out":3918,"would_cite":false,"duration_ms":45080,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["98.80.-k","95.36.+x","04.50.Kd"],"model":"deepseek-v4-flash","headline":"A log-periodic deceleration parameter in f(R,G) gravity fits cosmic-chronometer, BAO, and supernova data with H0 = 71.7–72.8 km/s/Mpc, q0 ≈ −0.5, and a transition redshift near 0.74–0.88.","keywords":["f(R,G) gravity","deceleration parameter","log-periodic parametrization","MCMC constraints","cosmological parameters","Hubble constant","energy conditions","thermodynamic analysis"],"falsifier":"Re-run the fit with $\\alpha$, $\\beta$, $\\gamma$ as free parameters: if the best-fit values move well away from $(0.5,0.6,0.7)$ or the resulting $H_0,q_0,q_1$ shift outside the reported ranges, the claimed viability of this $f(R,G)$ model would not hold up. A second check is to measure $q(z)$ at $z\\gtrsim2$ with better precision; the log-periodic ansatz predicts a specific oscillatory pattern that a monotonic expansion history would rule out.","tokens_in":17020,"feed_emoji":"🌌","tokens_out":11012,"duration_ms":112050,"temperature":0.7,"pith_summary":"The paper sets out to show that the observed late-time cosmic acceleration can be described by a modified gravity model $f(R,G)=R+\\alpha R^2+\\beta e^{\\gamma G}$ whose deceleration parameter oscillates logarithmically, $q(z)=q_0+q_1 \\sin[\\log(1+z)]$. From this ansatz it reconstructs the Hubble rate $H(z)$ and fits the three free parameters $H_0$, $q_0$, $q_1$ to cosmic-chronometer, BAO, and 1701-point supernova data. The fit returns $H_0$ between 71.7 and 72.8 km/s/Mpc, a present deceleration parameter $q_0$ between $-0.484$ and $-0.517$, and a transition from deceleration to acceleration at $z_{tr}$ between 0.879 and 0.744. If the reconstruction is right, the model also yields a radiation-era equation of state $\\omega\\approx 0.33$, a present $\\omega_0\\approx -0.49$, satisfies NEC and DEC while violating SEC at late times, and produces a cosmic age of 13.01–13.59 Gyr. The point of the exercise is that a single compact, model-independent parametrization of $q(z)$ can tie the expansion history to a specific $f(R,G)$ gravity and pass a battery of physical-consistency checks.","feed_headline":"Log-periodic expansion model pushes H0 toward 72 and q0 toward -0.5","feed_subtitle":"The same fit gives a transition redshift near 0.74 and a universe about 13 billion years old.","key_machinery":"The central object is the deceleration parameter $q(z)=q_0+q_1\\sin[\\log(1+z)]$, whose logarithmic sine makes the expansion oscillatory on the log-redshift scale while remaining finite as $z\\to-1$. Its closed-form integral, Eq. (13), fixes the kinematics; the model action $f(R,G)=R+\\alpha R^2+\\beta e^{\\gamma G}$ then turns that kinematics into dynamics through the modified Friedmann equations (8)–(9). Everything downstream—density, pressure, equation of state, energy conditions, statefinder pair, cosmic age, temperature, and entropy density—is a direct function of Eq. (13) with the fitted $H_0,q_0,q_1$ and the hand-set $\\alpha,\\beta,\\gamma$.","core_discovery":"On the paper's own terms, the central discovery is the reconstruction of the late-time Hubble expansion from the log-periodic deceleration ansatz: substituting $q(z)=q_0+q_1\\sin[\\log(1+z)]$ into the integral $H(z)=H_0\\exp[\\int_0^z (1+q(z'))/(1+z') dz']$ gives the closed form $H(z)=H_0(1+z)^{1+q_0}\\exp[q_1(1-\\cos[\\log(1+z)])]$ (Eq. 13). With $\\alpha$, $\\beta$, $\\gamma$ held at the chosen values, this $H(z)$ feeds the $f(R,G)$ field equations to produce analytic expressions for $\\rho(z)$ and $p(z)$, from which the equation of state, energy conditions, age integral, and thermodynamic quantities are derived. The paper takes the MCMC-constrained $H_0$, $q_0$, $q_1$, the negative $q_0$, the $z_{tr}$ shift, the $\\omega(z\\gg1)\\approx0.33$ limit, and the satisfaction/violation pattern of energy conditions as evidence that the model offers a viable and observationally consistent description of cosmic acceleration.","pith_inferences":["A natural next step the authors do not take is to constrain $\\alpha$, $\\beta$, $\\gamma$ with the same data; the reported $H_0$ and $q_0$ ranges could shift, and the exact $\\rho(z)$, $p(z)$ curves would carry parameter uncertainties rather than fixed values.","Because the deceleration parameter oscillates on the $\\log(1+z)$ scale, the model predicts small recurring acceleration–deceleration episodes at high redshift; higher-precision cosmic-chronometer or BAO measurements at $z\\gtrsim2$ could look for this signature.","The radiation-era limit $\\omega(z\\gg1)\\approx0.33$ is derived with the hand-chosen coefficients, so it is not a parameter-free prediction of the model; a fuller treatment would show how early-time behavior depends on $\\alpha$, $\\beta$, $\\gamma$.","The statefinder trajectory approaching the $\\Lambda$CDM point from the quintessence side suggests this parametrization could be used as a template for distinguishing modified gravity from a cosmological constant, but that use goes beyond what the paper claims."],"forward_implications":["Within the reconstruction, the Hubble constant lands at $H_0 = 71.7$–$72.8$ km/s/Mpc, closer to local distance-ladder measurements than to CMB-based values, so the model offers a route toward easing the Hubble tension.","The transition redshift moves from 0.879 with cosmic-chronometer data alone to 0.744 when BAO and supernova data are added, meaning the model predicts a later, sharper onset of acceleration as more data are included.","The equation-of-state parameter evolves from $\\omega\\approx0.33$ at high redshift to $\\omega_0\\approx -0.48$ to $-0.49$, so the same model covers radiation-dominated early behavior and quintessence-like late behavior.","The energy conditions come out as NEC and DEC satisfied with SEC violated at late times, which is the pattern expected for an accelerating fluid; the derived age of 13.01–13.59 Gyr and increasing total entropy keep the model within standard cosmological and thermodynamic bounds."],"supporting_citations":[{"why":"Type Ia supernova evidence for late-time accelerated expansion; supplies the phenomenon the model must reproduce.","marker":"[1]"},{"why":"Source of the $f(R,G)$ action and gravitational field equations on which the reconstruction is built.","marker":"[19]"},{"why":"Supplies the 31 cosmic-chronometer $H(z)$ measurements used in the CC likelihood.","marker":"[53]"},{"why":"Supplies the baryon acoustic oscillation measurements used for the 26 BAO points.","marker":"[56]"},{"why":"Supplies the method for computing the sound horizon entering the BAO distance ratios.","marker":"[57]"},{"why":"Supplies the 1701-light-curve supernova sample used in the distance-modulus likelihood.","marker":"[61]"},{"why":"Supplies the local-distance $H_0 = 73.4 \\pm 1.04$ km/s/Mpc baseline that the paper compares against as a late-time measurement.","marker":"[62]"},{"why":"Supplies the CMB-based $H_0 = 67.4 \\pm 0.5$ km/s/Mpc value whose tension with late-time measurements motivates the model.","marker":"[63]"},{"why":"Comparison f(R) model result with $H_0\\approx69.4$, $q_0\\approx-0.53$ used to benchmark the parameter discussion.","marker":"[64]"},{"why":"Recent modified-gravity result with $H_0\\approx72.5$, $q_0\\approx-0.48$ used to claim consistency with other frameworks.","marker":"[69]"}],"fun_headline_variants":["Log-periodic model sets H0 near 72","Oscillating q pins H0 to 72 and q0 to -0.5","f(R,G) gravity with log-periodic q predicts H0=72","Log-periodic deceleration yields H0 ~72, q0 ~ -0.5","Log-periodic q model predicts H0=72 and age 13 Gyr"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model's physical predictions all assume $\\alpha=0.5$, $\\beta=0.6$, $\\gamma=0.7$, which are chosen by hand in the caption of Figure 4 and never varied in the statistical fit; if those coefficients differ, the reported density, pressure, energy-condition, and thermodynamic behavior would change.","fun_headline_variants_meta":{"raw":{"variants":["Log-periodic model sets H0 near 72","Oscillating q pins H0 to 72 and q0 to -0.5","f(R,G) gravity with log-periodic q predicts H0=72","Log-periodic deceleration yields H0 ~72, q0 ~ -0.5","Log-periodic q model predicts H0=72 and age 13 Gyr"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001538,"raw_usage":{"total_tokens":6252,"prompt_tokens":1145,"completion_tokens":5107,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":761,"completion_tokens_details":{"reasoning_tokens":5003}},"tokens_in":761,"tokens_out":5107,"duration_ms":40933,"temperature":1.0,"reasoning_tokens":5003,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T10:33:11.414473+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the fit with $\\alpha$, $\\beta$, $\\gamma$ as free parameters: if the best-fit values move well away from $(0.5,0.6,0.7)$ or the resulting $H_0,q_0,q_1$ shift outside the reported ranges, the claimed viability of this $f(R,G)$ model would not hold up. A second check is to measure $q(z)$ at $z\\gtrsim2$ with better precision; the log-periodic ansatz predicts a specific oscillatory pattern that a monotonic expansion history would rule out.","supporting_citations":[{"cited_title":"Riess, et al., Astron","cited_arxiv_id":null,"evidence_quote":"Type Ia supernova evidence for late-time accelerated expansion; supplies the phenomenon the model must reproduce."},{"cited_title":"De Laurentis, M","cited_arxiv_id":null,"evidence_quote":"Source of the $f(R,G)$ action and gravitational field equations on which the reconstruction is built."},{"cited_title":"Di Valentino et al., Astropart","cited_arxiv_id":null,"evidence_quote":"Supplies the local-distance $H_0 = 73.4 \\pm 1.04$ km/s/Mpc baseline that the paper compares against as a late-time measurement."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Comparison f(R) model result with $H_0\\approx69.4$, $q_0\\approx-0.53$ used to benchmark the parameter discussion."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Recent modified-gravity result with $H_0\\approx72.5$, $q_0\\approx-0.48$ used to claim consistency with other frameworks."}],"review_version":1}