{"id":"4949a4e1-a65b-4db2-9e48-887742647588","arxiv_id":"2506.15789","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Using a parametrized deceleration parameter, the authors constrain H(z)=H0(1+az+bz^2) with cosmic chronometer and Pantheon data, then compute dark energy and energy conditions in two f(R,L_m) models.","lead":"The paper fits a simple quadratic formula for the cosmic expansion rate to galaxy-age and supernova data, then plugs the fit into two modified gravity models to infer dark-energy-like properties. A generalist reading shows how an assumed expansion history, not the gravity theory itself, drives the reported results.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Model I's Friedmann equations (26)-(27) drop the time-derivative terms of F_R from Eqs. (14)-(15), invalidating all derived Model I density, pressure, EoS, and energy-condition results.","rationale":"The kinematic part of the paper is straightforward: the q(z) parametrization integrates to H(z)=H0(1+az+bz^2), and Model II's reduced equations (29)-(30) do follow from Eqs. (14)-(15). The reader's concern about the assumed q(z) form is a legitimate model-selection issue, but it does not identify an internal error. The most serious defect is in Model I's Section 5.1: the time-derivative of F_R is dropped without justification. This is a direct algebraic inconsistency, checkable by substitution, and it does not depend on external assumptions or priors. Since the derived rho, p, omega, and energy conditions for Model I are the paper's main evidence of physical viability within the f(R,L_m) framework, the omission undermines the central claim as written. Correcting the equations would require solving differential equations for rho and p, likely changing all Model I results; the current version therefore cannot be accepted as-is.","tokens_in":17335,"tokens_out":17390,"duration_ms":140725,"concrete_test":"Recompute Eq. (26) from Eq. (14) for f=R/2+(1+eta R)rho, keeping the 3H dot(F_R) term. Verify that the result contains -6H eta dot(rho). Then substitute the reported Eq. (31) for rho(z) into both versions using the joint best-fit parameters (H0=68.7, a=0.458, b=0.313) and eta=1.03; if the corrected equation has a residual of order 6H eta dot(rho) while the printed equation is exactly satisfied, the Model I derivation is confirmed inconsistent.","verdict_should_be":"REJECT","load_bearing_attack":"In Section 5.1, the authors derive Friedmann equations for Model I, f(R,L_m)=R/2+(1+eta R)L_m with L_m=rho. Starting from the paper's own Eq. (14), with F_R=1/2+eta rho, the 3H dot(F_R) term contributes 3H eta dot(rho). Substituting F - F_Lm L_m - F_R R = -eta rho R and R=6(dot H + 2H^2) gives the correct first Friedmann equation: 3H^2(6eta rho - 1) + 12eta dot(H) rho + rho - 6H eta dot(rho) = 0. The printed Eq. (26) omits the -6H eta dot(rho) term; similarly, Eq. (27) omits 2eta ddot(rho) + 6H eta dot(rho). Consequently, Eq. (31) for rho(z) is not a solution of the stated field equations; it solves an algebraic equation obtained by setting dot(rho)=0. All derived Model I pressure, EoS, and energy-condition results (Figs. 4-12 and Section 6) inherit this error. Because the central viability claim is explicitly stated for both models, Model I's invalid derivation removes a major part of the support for the conclusion.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies two f(R,L_m) gravity models, f(R,L_m)=R/2+(1+\\eta R)L_m and f(R,L_m)=R/2+L_m^\\eta, in a flat FLRW universe with a perfect fluid. It adopts the two-parameter deceleration parameter q(z) of Eq. (16), integrates it to H(z)=H_0(1+a z+b z^2), and fits (H_0,a,b) and the supernova magnitude offset M to 31 cosmic chronometer data points and to a joint CC+Pantheon sample using MCMC. It then derives the energy density, pressure, equation-of-state parameter, energy conditions, cosmographic parameters, and cosmic age for each model. The central claim is that the chosen q(z) parametrization inside f(R,L_m) gravity provides a viable and compelling account of the observed late-time acceleration.","tokens_in":17667,"tokens_out":14405,"duration_ms":129029,"significance":"If all results were correct, the paper would give a compact phenomenological description of late-time acceleration in modified gravity, with posterior constraints on H_0, the transition redshift, and derived dark-energy quantities. The strengths are the transparent parametrization, the use of standard CC and Pantheon datasets, and the fact that the Model II Friedmann equations, Eqs. (29)-(30), are derived correctly from the stated field equations. However, the paper's central viability claim currently rests on Model I results that do not follow from the field equations, on a hand-set coupling \\eta=1.03, and on derived quantities that are algebraic consequences of the assumed q(z) rather than dynamical predictions of f(R,L_m). The significance is therefore moderate and conditional on substantial revision.","major_comments":[{"comment":"The Friedmann equations for Model I omit the time-derivative terms of F_R. Substituting f=R/2+(1+\\eta R)L_m with L_m=\\rho into the paper's own Eq. (14) and using R=6(\\dot H+2H^2), I obtain 3H^2(1-6\\eta\\rho)-12\\eta\\rho\\dot H+6H\\eta\\dot\\rho-\\rho=0, whereas the printed Eq. (26) is equivalent to 3H^2(1-6\\eta\\rho)-12\\eta\\rho\\dot H-\\rho=0, i.e., the 6H\\eta\\dot\\rho term is missing. Similarly, Eq. (27) drops the 2\\eta\\ddot\\rho+6H\\eta\\dot\\rho terms. Consequently, Eq. (31), Eq. (32), the EoS parameter Eq. (35), and all Model I energy-condition plots in Section 6.4 do not solve the stated field equations; they solve the algebraic system obtained by setting \\dot\\rho=\\ddot\\rho=0. This invalidates a major part of the support for the Section 7 conclusion that both models are viable.","section":"Section 5.1, Eqs. (14)-(15), (26)-(27)"},{"comment":"The value \\eta=1.03 is introduced by hand in Section 6.2 with no observational constraint, no error propagation, and no sensitivity analysis. Since \\eta is a free parameter in both models and enters Model II as an exponent in Eqs. (29)-(30), the derived energy densities, pressures, EoS parameters, and energy-condition transition redshifts all depend strongly on this arbitrary choice. The paper should constrain \\eta jointly with H_0,a,b and M, marginalize over it, or at the minimum show how the derived quantities change with \\eta.","section":"Section 6.2 and Section 5"},{"comment":"The conclusion that the chosen parametrization 'provides a viable and compelling approach' to account for late-time acceleration is overstated relative to what is tested. The fitted quantities q_0, z_t, j_0, s_0, and the derived \\omega_0 are deterministic functions of the assumed q(z) family in Eq. (16) and the fitted parameters a and b; they do not test the f(R,L_m) dynamics. A comparison with \\LambdaCDM or with alternative q(z) parametrizations, together with a model-selection statistic, would be needed to support the claim that the model is compelling rather than merely consistent with the data.","section":"Sections 6.1, 6.5, and 7"}],"minor_comments":[{"comment":"The joint-data best fit is listed as a=0.458 in Table 1, but Section 7 prints a=-0.458 with a minus sign; this typo should be corrected.","section":"Section 7 and Table 1"},{"comment":"The text states that setting \\eta=0 recovers the conventional Friedmann equations of GR for Model II, but for f=R/2+L_m^\\eta the GR limit is \\eta=1, not \\eta=0; \\eta=0 gives R/2+1, which is GR with a cosmological constant rather than the standard Friedmann equations.","section":"Section 5.2"},{"comment":"The sentence saying that the energy-condition results are \"consistent with the fulfillment ... of the fundamental energy conditions (EC), namely the Null, Weak, Dominant and Strong EC\" is internally contradictory, since the same paragraph states that the Strong EC is violated; it should say all conditions except the SEC are satisfied.","section":"Section 6.4"},{"comment":"The derived quantities q_0, z_t, j_0, and s_0 are quoted without credible intervals, even though they come from the MCMC chains; the paper should propagate the posterior uncertainties to these quantities.","section":"Table 1"}],"recommendation":"major_revision","confidential_remarks":"The Model I derivation error in Section 5.1 is serious and affects one of the two models central to the paper's conclusion. The authors should either correct the Friedmann equations by retaining the time-derivative terms and re-deriving the physical quantities, or explicitly restrict the viability claim to Model II. The hand-set value \\eta=1.03 also needs to be justified or marginalized. With these changes the paper may become publishable, but in its current form the central claim is not sufficiently supported."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the stress-test note is correct and it is fatal for half of the paper. Starting from the paper's own Eq. (14) with f = R/2 + (1 + eta R)rho, you get F_R = 1/2 + eta rho and a non-zero 3H dot(F_R) = 3H eta dot(rho). The printed Eq. (26) drops that term, and Eq. (27) drops the corresponding second-order terms. So the density in Eq. (31), and everything downstream for Model I (pressure, EoS, energy conditions, Figs. 4-12), solves an algebraic equation with dot(rho)=0 rather than the actual field equations. That is not a stylistic slip; it is an internal contradiction with the authors' own starting equations. I checked it myself.\n\nWhat is actually good: Model II, with f = R/2 + L_m^eta, does follow correctly from Eqs. (14)-(15). The MCMC constraints on H0, a, b, M are reported clearly, the H0 values are plausible, and the integration from q(z) to H(z) is algebraically right. The citation pattern looks normal for this subfield.\n\nNow the soft spots, in proportion. Even where the algebra holds, the paper does not test f(R,L_m) gravity. The quadratic H(z) is assumed, so q0, zt, jerk, snap, and omega are algebraic functions of the fitted a and b; they cannot fail to produce acceleration. Calling that a prediction is circular. Also: eta = 1.03 is set by hand with no observational constraint, most derived quantities have no error bars, and there is no model comparison such as AIC or BIC against Lambda CDM. There is also a sign inconsistency: the joint-data value of a is 0.458 in Table 1 but appears as -0.458 in the conclusions. These are fixable, but together with the Model I error they make the central conclusion unsupported.\n\nWho is this for? Readers tracking f(R,L_m) phenomenology might want the valid Model II part after major revision. But in current form I would not spend referee time on it: the load-bearing error is demonstrable in-house, and the novelty is only a new combination of known elements. If the authors correct Model I, constrain eta, and reframe the fits as constraints rather than predictions, it could become a modest acceptable paper.","headline":"Model I's Friedmann equations contradict the paper's own field equations: the 3H dot(F_R) term is dropped, so the headline viability claim for Model I is void; the rest is a standard kinematic fit with a circular interpretation.","tokens_in":18185,"tokens_out":8600,"would_cite":false,"duration_ms":78877,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83F05","83D05"],"pacs":["98.80.-k","04.50.Kd"],"model":"deepseek-v4-flash","headline":"This paper claims that a two-parameter deceleration ansatz, integrated to a quadratic Hubble law inside $f(R,L_m)$ modified gravity, accounts for late-time acceleration with a transition at $z_t=0.655$.","keywords":["late-time cosmic acceleration","deceleration parameter","f(R,L_m) gravity","quadratic Hubble law","equation of state","energy conditions","cosmographic parameters","cosmic chronometer data"],"falsifier":"Fit $H(z)=H_0(1+az+bz^2)$ to the same 31 cosmic-chronometer and 1048 Pantheon points and examine the residuals as a function of redshift. If a nonparametric or higher-order reconstruction of $H(z)$ deviates from the best-fit quadratic by more than the quoted uncertainties, for example by requiring an inflection or a second transition in $q(z)$ away from $z\\approx 0.655$, the central claim fails. A direct calculation would be to redo the MCMC with a three-parameter or nonparametric $q(z)$ and compare model selection; the claim requires that the two-parameter quadratic match or beat those alternatives.","tokens_in":2070,"feed_emoji":"🌌","tokens_out":5625,"duration_ms":105728,"temperature":0.7,"pith_summary":"This paper tries to establish that a simple kinematic assumption inside $f(R,L_m)$ modified gravity can describe the observed late-time acceleration of the universe. The assumption is a two-parameter deceleration parameter that integrates to the quadratic Hubble law $H(z)=H_0(1+az+bz^2)$. Fitting this expansion history to cosmic-chronometer and Pantheon supernova data gives a present-day deceleration parameter $q_0\\approx -0.46$ (CC) or $-0.54$ (CC+Pantheon) and a transition from deceleration to acceleration at $z_t=0.655$. The same expansion history, inserted into two nonlinear $f(R,L_m)$ models, yields positive energy density, negative late-time pressure, an equation-of-state parameter in the quintessence range for one model, violation of the strong energy condition, and a cosmic age near 13 Gyr. If correct, the paper shows that this parametric family is a workable description of the acceleration epoch and that $f(R,L_m)$ gravity can reproduce its main features.","feed_headline":"Cosmic acceleration switch pinned at z = 0.655","feed_subtitle":"A two-parameter deceleration law in modified gravity fits chronometer and supernova data.","key_machinery":"The load-bearing object is the two-parameter deceleration ansatz, Eq. (16), a rational function of redshift engineered so that integrating $\\dot H=-(1+q)H^2$ yields exactly the quadratic Hubble law $H(z)=H_0(1+az+bz^2)$. This quadratic law is the entire expansion history: every derived quantity, including energy density, pressure, equation of state, energy conditions, jerk, snap, and cosmic age, is a function of that quadratic and of the Friedmann equations of the two chosen gravity models. The two $f(R,L_m)$ forms supply the modified Friedmann equations, Eqs. (26)-(27) and Eqs. (29)-(30), that convert the kinematic $H(z)$ into fluid properties and dark-energy behavior.","core_discovery":"The central discovery is that the ansatz $q(z)=-1+\\frac{(1+z)(a+2bz)}{1+az+bz^2}$ is exactly equivalent, through $H(z)=H_0\\,\\exp\\!\\left(\\int_0^z \\frac{1+q(x)}{1+x}\\,dx\\right)$, to the normalized Hubble parameter $H(z)/H_0=1+az+bz^2$. With median parameters from MCMC fits, $q(z)$ crosses zero at $z_t=0.655$ for both datasets, so the paper's kinematic core places the present epoch at $q_0<0$, a decelerated matter-dominated past with $q\\to\\tfrac12$, and a future de Sitter-like phase with $q\\to-1$. Feeding this $H(z)$ into the modified Friedmann equations of two $f(R,L_m)$ models, one with nonminimal curvature-matter coupling and one minimal $L_m^\\eta$ coupling, gives positive energy density, currently negative pressure, present equation-of-state values $\\omega_0\\approx -0.28$ and $-0.39$ for Model I and $-0.63$ and $-0.69$ for Model II, and violation of the strong energy condition. The paper therefore asserts that the chosen parametric deceleration form within $f(R,L_m)$ gravity is a viable account of the observed late-time cosmic acceleration.","pith_inferences":["Editorial inference: because $H(z)=H_0(1+az+bz^2)$ is fixed purely by the kinematic ansatz before any gravity theory is chosen, the fitted values of $q_0$ and $z_t$ would be unchanged in any modified gravity that adopts the same $q(z)$; the $f(R,L_m)$ analysis recasts that same expansion history in terms of modified fluid densities and pressures rather than independently testing the gravity theory","A testable extension the paper does not perform is to treat $\\eta$ as a free parameter in the MCMC fit instead of fixing $\\eta=1.03$, and to add other cosmological datasets, in order to see whether the gravity-model parameters and the $q(z)$ parameters remain mutually consistent.","The two-parameter quadratic $H(z)$ could be checked against nonparametric reconstructions: if $H(z)/H_0$ shows curvature beyond a quadratic, or if the deceleration parameter has a second transition, then the assumption underlying all derived dark-energy results would fail."],"forward_implications":["The fitted expansion history places the deceleration-to-acceleration transition at $z_t\\approx 0.655$ and the present deceleration parameter at $q_0\\approx -0.46$ (CC) or $-0.54$ (CC+Pantheon), so the model's kinematic core is compatible with late-time acceleration.","Model II, with minimal coupling $f(R,L_m)=R/2+L_m^\\eta$, gives a present equation-of-state parameter $\\omega_0\\approx -0.63$ to $-0.69$, in the quintessence region and closer to dark-energy behavior than Model I.","The strong energy condition is violated in both models, as required for accelerated expansion, while the null, weak, and dominant energy conditions remain satisfied.","The model returns a cosmic age $t_0\\approx 12.8$ to $13.0$ Gyr, and the joint dataset gives a jerk parameter $j_0\\approx 0.92$, close to the $\\Lambda$CDM value $j_0=1$, while the CC-only jerk deviates more."],"supporting_citations":[{"why":"Defines $f(R,L_m)$ gravity and gives the action and field equations used for both cosmological models.","marker":"[24]"},{"why":"Supplies the MCMC sampler used to obtain the median parameter values and posterior contours.","marker":"[43]"},{"why":"Provides cosmic-chronometer Hubble-parameter measurements used in the CC fit.","marker":"[44]"},{"why":"Provides the compiled cosmic-chronometer dataset used for the Hubble-parameter fit.","marker":"[45]"},{"why":"Gives the differential-age relation $H(z)=-\\frac{1}{1+z}\\frac{dz}{dt}$ that converts galaxy ages into Hubble data.","marker":"[48]"},{"why":"Supplies the Pantheon compilation of 1048 Type Ia supernovae used in the joint CC+Pantheon fit.","marker":"[51]"},{"why":"Supplies the generalized curvature-matter coupling framework from which Model II, $f(R,L_m)=R/2+L_m^\\eta$, is drawn.","marker":"[59]"},{"why":"Supports the choice $L_m=\\rho$ used to derive the modified Friedmann equations and the thermodynamic quantities.","marker":"[61]"}],"fun_headline_variants":["Deceleration to acceleration pivot at z=0.655 in f(R,L_m) gravity","Modified gravity with q(z) crossing zero fits cosmic acceleration","Parametric q(z) in f(R,L_m) theory explains late-time acceleration","Cosmic acceleration pivot at z=0.655 in modified gravity models","f(R,L_m) gravity with quadratic Hubble law matches cosmic data"],"cache_read_input_tokens":20224,"weakest_assumption_plain":"The load-bearing premise is that the true deceleration history is exactly the two-parameter rational form $q(z)=-1+\\frac{(1+z)(a+2bz)}{1+az+bz^2}$, equivalently $H(z)=H_0(1+az+bz^2)$; the paper also assumes $L_m=\\rho$ and fixes $\\eta=1.03$ by hand, and none of these choices is tested against alternative forms.","fun_headline_variants_meta":{"raw":{"variants":["Deceleration to acceleration pivot at z=0.655 in f(R,L_m) gravity","Modified gravity with q(z) crossing zero fits cosmic acceleration","Parametric q(z) in f(R,L_m) theory explains late-time acceleration","Cosmic acceleration pivot at z=0.655 in modified gravity models","f(R,L_m) gravity with quadratic Hubble law matches cosmic data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000697,"raw_usage":{"total_tokens":3193,"prompt_tokens":1032,"completion_tokens":2161,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":648,"completion_tokens_details":{"reasoning_tokens":2063}},"tokens_in":648,"tokens_out":2161,"duration_ms":14643,"temperature":1.0,"reasoning_tokens":2063,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:31:39.372349+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fit $H(z)=H_0(1+az+bz^2)$ to the same 31 cosmic-chronometer and 1048 Pantheon points and examine the residuals as a function of redshift. If a nonparametric or higher-order reconstruction of $H(z)$ deviates from the best-fit quadratic by more than the quoted uncertainties, for example by requiring an inflection or a second transition in $q(z)$ away from $z\\approx 0.655$, the central claim fails. A direct calculation would be to redo the MCMC with a three-parameter or nonparametric $q(z)$ and compare model selection; the claim requires that the two-parameter quadratic match or beat those alternatives.","supporting_citations":[],"review_version":2}