{"id":"582b3772-5e8c-43fd-81fe-85e06757308b","arxiv_id":"2412.03065","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"high","formal_verification":"none","parameter_count":7,"one_line_summary":"A two-variable Chebyshev polynomial reconstruction of f(Q,T) gravity is shown to be algebraically identical to a truncated Taylor series, then fitted to Pantheon+SH0ES data.","lead":"The paper fits a Taylor expansion of a modified gravity function f(Q,T) to supernova data using a Chebyshev-polynomial re-writing of the same series. The fit matches Lambda-CDM, but the claimed new Chebyshev method is actually a re-expression of a Taylor series, with parameters fitted rather than predicted.","discovery_kind":"incremental","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (23) sets T = +rho for dust, but under the stated (-,+,+,+) signature the trace is T = -rho; this sign error propagates into the reconstructed f(Q,T) coefficients via Eqs. (37)-(41).","rationale":"The paper's main contribution is a two-variable Chebyshev reconstruction of f(Q,T). For this to be correct, the relation between the fitted coefficients and the theory's variables must be consistent. The reader's identified sign error is exactly such an inconsistency: the trace T is defined by contraction, so its sign is fixed by the metric signature; using T=+rho for dust under (-,+,+,+) is factually wrong. Since all the reconstructed coefficients are derived from T and its derivatives, this is the most load-bearing concern. The alternative concern that the Chebyshev expansion reduces to the Taylor series is valid but is a novelty or statistical issue; even if the method is not new, the reconstructed function could in principle be correct if the algebra were consistent. The sign error, by contrast, directly invalidates the central reconstruction. I therefore agree with the reader that REJECT is appropriate. The paper does include an MCMC analysis and AIC/BIC comparison, which are useful, but they do not overcome the sign inconsistency because the fitted parameters are interpreted as f(Q,T) coefficients under the wrong trace convention.","tokens_in":14084,"tokens_out":12101,"duration_ms":108871,"concrete_test":"Evaluate T = g^{mu nu} T_{mu nu} for a dust fluid in metric (16) with signature (-,+,+,+). If the result is -3 H0^2 Omega_m0 / a^3, recompute T', T'', T''' with the corrected sign and re-derive Eqs. (37)-(41) and the resulting MCMC constraints in Table I. If the best-fit eta^(1) and eta^(2) move by more than their 1-sigma error bars, or if the chi^2 changes, the paper's reconstructed f(Q,T) is not the one claimed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In the FLRW metric (16) with the stated signature (-,+,+,+), the trace of the dust energy-momentum tensor is T = g^{mu nu} T_{mu nu} = -rho, with rho = 3 H0^2 Omega_m0 / a^3 in the paper's units. Equation (23) instead sets T = +rho. This is not a cosmetic convention issue: T, T', T'', and T''' in Eq. (23) enter the derivation of the present-day field equations and the cosmographic parameters in Section IV C. The coefficients eta^(1) and eta^(2) in Eqs. (37)-(41) are derivatives of eta(T) with respect to the trace; if the physical trace is -rho, these derivatives are evaluated at the wrong argument and their signs relative to rho are inverted. Consequently, the fitted coefficients in Table I do not describe the function f(Q,T) for the theory as defined by the paper. The central claim that a functional form of f(Q,T) has been reconstructed therefore fails for the stated theory; fixing the sign would change the reconstruction.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a cosmographic reconstruction of f(Q,T) gravity using two-variable Chebyshev polynomials. It assumes the minimally coupled form f(Q,T)=γ(Q)+η(T), expands both functions in Taylor series around present-day values, rewrites the expansion as Eq. (32), and uses the f(Q,T) Friedmann equations to express the deceleration, jerk, and snap parameters in terms of the Taylor coefficients (Eqs. (39)-(41)). It then constructs a luminosity-distance expression from the standard cosmographic Taylor series, re-expresses it in terms of Chebyshev polynomials, and performs an MCMC fit to Pantheon+SH0ES data with free parameters H0, γ(1), γ(3), γ(4), η(1), η(2), and Ωm0. The authors report a distance modulus that closely follows ΛCDM and quote AIC/BIC values. The central claims are that this yields a model-independent two-variable Chebyshev reconstruction of f(Q,T) and that the reconstructed theory is consistent with current supernova data.","tokens_in":14338,"tokens_out":9858,"duration_ms":100697,"significance":"If the advertised method worked, it would extend single-variable cosmography to a two-variable Chebyshev reconstruction of f(Q,T) and would provide a genuinely model-independent route to fixing the functional form. The manuscript is clearly organized, derives the standard f(Q,T) field equations, and presents a complete MCMC pipeline with posterior contours and a comparison with Pantheon+SH0ES. However, two load-bearing problems undermine the central claim. First, Eq. (23) sets the dust trace T to +ρ, whereas with the stated (-,+,+,+) signature the trace is -ρ; since η(1), η(2) and all higher trace derivatives enter Eqs. (37)-(41), the fitted coefficients in Table I do not describe f(Q,T) for the theory defined in the paper. Second, the 'Chebyshev reconstruction' in Eq. (32) is algebraically identical to the Taylor expansion (31), and the luminosity distance (35)-(36) is a polynomial basis change of the Taylor series (34), so the Chebyshev machinery does no new work.","major_comments":[{"comment":"The trace of the dust energy-momentum tensor for the FLRW metric (16) with signature (-,+,+,+) is T = g^{μν}T_{μν} = -ρ, where ρ = 3H0^2 Ωm0/a^3. Equation (23) instead sets T = +3H0^2 Ωm0/a^3. This is not a cosmetic convention choice: T, T', T'', and T''' from Eq. (23) enter the Taylor expansion (31), the present-day field equations (37)-(38), and the cosmographic relations (39)-(41) through η(T) and its derivatives. With T = -ρ, the derivatives η(1), η(2) are evaluated at an argument of opposite sign from ρ, so the fitted values in Table I are not the coefficients of f(Q,T) for the theory defined by Eq. (8). The authors must either correct the sign and redo the reconstruction or explicitly adopt and justify a nonstandard convention for T; as it stands, the central result is not a reconstruction of the stated f(Q,T) theory.","section":"Section III, Eq. (23)"},{"comment":"Equation (32) is not a Chebyshev expansion. The Chebyshev coefficients αi,j defined in Eq. (30) never appear in the derivation; Eq. (32) is simply the Taylor polynomial (31) rewritten in nested form, as can be verified by expanding (Q-Q0) and (T-T0). The phrase 'By incorporating (30) and (31) in the Chebyshev series (29)' therefore does not describe what is actually done. To support the central claim of a two-variable Chebyshev reconstruction, the authors would need to compute the Chebyshev coefficients from an appropriately normalized domain and show that the resulting series differs from the Taylor truncation. As it stands, the advertised method reduces to a Taylor expansion of f(Q,T) around (Q0,T0).","section":"Section IV A, Eq. (32)"},{"comment":"The luminosity distance used in the likelihood is the standard Taylor series (34) re-expressed in Chebyshev polynomials. Since Tn(z) are polynomials of degree n, the expression in Eq. (35) is exactly a quartic polynomial in z and is algebraically identical to the Taylor truncation (34). The fit therefore constrains the kinematic coefficients H0, q0, j0, s0, and the agreement with Pantheon+SH0ES and ΛCDM in Fig. 1 is a property of this polynomial distance modulus, not an independent test of f(Q,T). The subsequent use of Eqs. (39)-(41) to convert the fitted kinematic parameters into f(Q,T) coefficients is a consistency inversion: the output functional form is dictated by the input Taylor model. This is the sense in which the reconstruction is circular by construction, and the Chebyshev basis change adds no new information.","section":"Section V, Eqs. (34)-(36)"},{"comment":"The statistical comparison is misreported. The paper states that ΔAIC = 1.88 'indicates strong evidence in favor of the model,' but a difference of 1.88 is at best weak-to-moderate support, especially for a model with seven free parameters. More importantly, the authors also report ΔBIC = 14.43, which under standard criteria is strong evidence against the model, not a neutral 'slightly higher' value. The concluding claim that the model makes 'an excellent match' to ΛCDM is therefore not supported by the paper's own information-criterion results.","section":"Section V B and Section VI"}],"minor_comments":[{"comment":"The Chebyshev coefficient integrals in Eq. (30) are over the square [-1,1]^2, but Q = 6H^2 and T = 3H0^2Ωm0/a^3 are not normalized to this interval and have physical dimensions. An affine rescaling of Q and T is needed before the coefficients αi,j are defined; otherwise the integrals and the weight function are not meaningful as written.","section":"Section IV A, Eq. (30)"},{"comment":"The polynomial variable z in Eq. (35) exceeds 1 for redshifts up to 2.26 in the Pantheon+SH0ES sample. Although the finite polynomial identity (33) holds for all z, the Chebyshev convergence properties invoked in Section IV A apply to functions on [-1,1]. To exploit those properties, the authors should map z (or dL) to a bounded variable such as (2z - zmax)/zmax.","section":"Section V A, Eq. (35)"},{"comment":"The zero-order constants γ and η are omitted from the MCMC fit and from Table I. Since only the combination γ+η appears in Eq. (37), the individual values of γ(Q0) and η(T0) are not determined, and the full functional form f(Q,T) is not reconstructed. The paper should state whether these constants are fixed by a convention or simply left free.","section":"Section IV C and Table I"},{"comment":"The notation mixes T for the trace of the matter energy-momentum tensor and Tn for Chebyshev polynomials; this makes several equations, such as Eq. (29), unnecessarily confusing. Renaming the trace, for example to τ, would improve readability.","section":"Throughout"}],"recommendation":"reject","confidential_remarks":"The sign error in Eq. (23) and the fact that Eq. (32) is not a Chebyshev expansion are independently disqualifying for the central claims. The first is fixable in principle, but the second points to a deeper problem: the advertised new method of two-variable Chebyshev reconstruction is not actually implemented, and the luminosity-distance fit is a polynomial expansion of the standard cosmographic series. The agreement with data is therefore an expected consequence of the assumed kinematic expansion rather than evidence for the reconstructed f(Q,T). I do not see a viable path to acceptance without a substantial reformulation of the method and a reanalysis with the correct trace sign. I recommend rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nYou should know two things about this paper. The central claim—a model-independent, two-variable Chebyshev reconstruction of f(Q,T)—does not survive contact with its own equations. And there is a sign error in the trace T that propagates into the reconstructed coefficients. This version is not a reliable reconstruction.\n\nWhat the paper does well is concrete: it actually writes the f(Q,T) functional and the luminosity distance in a two-variable Chebyshev basis and then fits the resulting parameters to Pantheon+SH0ES with MCMC, reporting best-fit ranges and AIC/BIC. That is a real exercise, and the numbers in Table I are internally consistent with the equations as written.\n\nThe soft spots are serious. Equation (31) is a Taylor series; Eq. (32) is exactly that same series after substituting Chebyshev polynomials. The 'reconstruction' is a basis re-expression of an assumed Taylor form, not an independent method. Likewise, Eq. (35) with coefficients (36) is the Taylor d_L of Eq. (34) rewritten via z^n in terms of Chebyshev polynomials. So the 'excellent match' to data is a self-consistency check, not an independent confirmation.\n\nThe trace error is load-bearing. With the stated (−,+,+,+) signature, for pressureless dust the trace is T = −ρ. Equation (23) sets T = +ρ. That sign enters the field equations (37)–(38) and the parameter relations (39)–(41), so the fitted coefficients do not describe f(Q,T) for the stated theory. Fixing the sign changes all of the reconstructed ranges.\n\nThe statistical claims are also overstated. ΔBIC = 14.43 is substantial evidence against the model, and calling ΔAIC = 1.88 \"strong evidence\" for it misreads the standard criteria.\n\nIf the sign is fixed and the claims are reframed as a Taylor-based cosmographic fit in a Chebyshev basis, the paper might become a modest methodological note. As it stands, the central reconstruction claim is circular and the numerical results are for the wrong trace. I would not send this to peer review; I'd return it to the authors with these points.","headline":"A Taylor series in a Chebyshev basis plus a sign error in the trace does not make a model-independent reconstruction.","tokens_in":14896,"tokens_out":4264,"would_cite":false,"duration_ms":40766,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83D05","83F05"],"pacs":["04.50.Kd","98.80.-k"],"model":"deepseek-v4-flash","headline":"The paper proposes a two-variable Chebyshev-polynomial reconstruction of the f(Q,T) gravity Lagrangian and reports that its distance modulus matches Pantheon+SH0ES supernovae and ΛCDM.","keywords":["Chebyshev polynomial","f(Q,T) gravity","cosmography","luminosity distance","Pantheon+SH0ES","non-metricity","MCMC","distance modulus"],"falsifier":"Take the same f(Q,T) field equations but use the standard dust trace T = −ρ, redo the Markov-chain Monte Carlo fit, and compare the best-fit ranges to Table I; if the ranges shift by more than the quoted 1-σ errors, the reconstruction as stated does not describe the theory. A direct check is to insert the fitted γ(1), γ(3), γ(4), η(1), η(2), and Ωm0 into Eqs. (18)–(19) and ask whether the resulting H(z) reproduces the Pantheon+SH0ES distance moduli within the published covariance.","tokens_in":13871,"feed_emoji":"🔭","tokens_out":9722,"duration_ms":86097,"temperature":0.7,"pith_summary":"This paper sets out to reconstruct the functional form of f(Q,T) gravity, an extended symmetric teleparallel theory in which the Lagrangian depends on non-metricity Q and the trace T of the matter energy-momentum tensor, directly from distance observations rather than from a pre-chosen model. It uses first-kind Chebyshev polynomials in two variables to expand f(Q,T) under a minimally coupled ansatz γ(Q)+η(T), and converts the standard Taylor luminosity distance into a Chebyshev series whose coefficients are explicit functions of the present-time cosmographic parameters H0, q0, j0, s0. Fitting that distance modulus to the 1701 Pantheon+SH0ES supernovae gives 1-sigma ranges for six free derivatives, with H0 around 73 km/s/Mpc and Ωm0 around 0.304. The paper's central claim is that the resulting curve matches both the data and the ΛCDM distance modulus across the entire sampled redshift range. A reader should care because this is a model-independent path from supernova distances to the coupling structure of a modified gravity theory.","feed_headline":"Chebyshev method reconstructs f(Q,T) gravity from supernova data","feed_subtitle":"A two-variable Chebyshev expansion matches Pantheon+SH0ES distance moduli and ΛCDM across the full redshift range.","key_machinery":"The load-bearing objects are the first-kind Chebyshev polynomials T_n(x)=cos(n arccos x), with orthogonality on [-1,1] and recurrence T_{n+1}=2xT_n(x)-T_{n-1}(x). The paper uses them in two roles: to expand the two-variable function f(Q,T) as Σ α_{i,j} T_i(Q) T_j(T), and to convert the Taylor luminosity distance into a Chebyshev series dL(z) = (c/H0) Σ_{n=0}^{4} c_n T_n(z), where c_n are explicit rational functions of q0, j0, s0. The non-metricity scalar Q=$6H^{2}$ and the assumed matter trace T=$3H0^{2}$ Ωm0/$a^{3}$ feed the f(Q,T) Friedmann equations, whose present-time solution ties the cosmographic parameters to the Chebyshev coefficients and lets the Markov-chain Monte Carlo run constrain the six derivatives.","core_discovery":"The paper claims that the two-variable Chebyshev series reconstructs the f(Q,T) Lagrangian as f(Q,T) ≈ γ(Q)+η(T) with γ and η expanded to fourth order in (Q−Q0) and (T−T0), and that the luminosity distance can be written as dL(z) = (c/H0) Σ_{n=0}^{4} c_n T_n(z), where c0,...,c4 are closed-form rational functions of q0, j0, s0 (with α=1/192). Substituting the field equations at the present time yields expressions for q0, j0, and s0 in terms of γ(1), γ(3), γ(4), η(1), η(2), and Ωm0, so a Markov-chain Monte Carlo fit to the Pantheon+SH0ES distance moduli pins down those unknowns. The reported best-fit ranges are H0 = 73.$0^{{+1.0}}$_{-0.87}, Ωm0 = 0.$304^{{+0.042}}$_{-0.019}, γ(1) = 12528.74 ± 0.99, γ(3) = −1080.2 ± 1.0, γ(4) = −18.92 ± 0.98, η(1) = 0.$19^{{+0.41}}$_{-0.74}, and η(2) = 0.$88^{{+0.32}}$_{-0.55}. The paper interprets the resulting distance modulus as an excellent match to the 1701 data points and to ΛCDM, and reports ΔAIC = 1.88 as strong evidence in favor while ΔBIC = 14.43 offers no supportive evidence.","pith_inferences":["The same two-variable Chebyshev coefficient scheme can be transplanted to other coupled theories, such as f(R,T) or f(Q,L_m), by replacing the trace and the matter Lagrangian; the paper only demonstrates f(Q,T).","The sign of T is convention-dependent: redoing the chain with T = −ρ, the standard dust trace under the declared signature, is likely to shift the fitted derivative ranges even if the qualitative match to ΛCDM survives, so the quoted numbers should be read within that convention.","The method's redshift reach could be tested by applying the fitted distance modulus to high-redshift probes such as quasars, gamma-ray bursts, or fast radio bursts beyond the z = 2.26 supernova ceiling; the paper does not perform that test.","The paper fits a kinematic reconstruction, not a full cosmological model, so further work would be needed to show that the reconstructed f(Q,T) also predicts the observed growth of structure and cosmic microwave background anisotropies."],"forward_implications":["The fitted coefficients assemble into an explicit, data-anchored functional form for f(Q,T), something the paper argues single-variable cosmography cannot do for coupled theories.","The Chebyshev luminosity distance formula with the stated coefficients is a ready-made model-independent distance expression that can be reused with other distance catalogs or priors.","Because Chebyshev series converge exponentially for analytic functions, the reconstruction is intended to remain valid beyond the z < 1 range where Taylor cosmography breaks down.","The distance modulus matches ΛCDM over 0.001 ≤ z ≤ 2.2613, so the reconstructed theory is consistent with standard cosmology at the kinematic level.","The AIC/BIC split (1.88 versus 14.43) means the model gains strong support on fit quality but is penalized by its large number of parameters.","The reconstructed f(Q,T) is kinematic, fit to distances; connecting it to structure growth or perturbation theory would require additional constraints the paper does not address."],"supporting_citations":[{"why":"It introduces f(Q,T) gravity, whose action, field equations, and Friedmann equations form the basis of the reconstruction.","marker":"[31]"},{"why":"It establishes Chebyshev-polynomial cosmography for one variable and the z-to-Chebyshev conversion that the paper extends to two variables.","marker":"[49]"},{"why":"It provides the two-variable Chebyshev series and coefficient formulas used to expand f(Q,T).","marker":"[56]"},{"why":"It defines the cosmographic parameters H, q, j, s and their Taylor relations used throughout the derivation.","marker":"[54]"},{"why":"It supplies the Pantheon+ sample, the covariance matrix, and the distance-modulus calibration entering the chi-square.","marker":"[60]"},{"why":"It provides the SH0ES Cepheid distance measurements and Hubble-constant anchor that define the Pantheon+SH0ES combination.","marker":"[58]"}],"fun_headline_variants":["Chebyshev cosmography constrains f(Q,T) gravity","Two-variable Chebyshev reconstructs f(Q,T) gravity","Chebyshev expansion fits f(Q,T) to supernova data","New Chebyshev method recovers f(Q,T) gravity","Chebyshev series matches f(Q,T) to Pantheon+SH0ES"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the matter trace is T = $3H0^{2}$ Ωm0/$a^{3}$ with a positive sign under the paper's (−,+,+,+) metric signature; with the standard dust trace −ρ, every fitted coefficient in Eqs. (37)–(41) would change, so a sign mistake would collapse the reconstruction.","fun_headline_variants_meta":{"raw":{"variants":["Chebyshev cosmography constrains f(Q,T) gravity","Two-variable Chebyshev reconstructs f(Q,T) gravity","Chebyshev expansion fits f(Q,T) to supernova data","New Chebyshev method recovers f(Q,T) gravity","Chebyshev series matches f(Q,T) to Pantheon+SH0ES"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000171,"raw_usage":{"total_tokens":1312,"prompt_tokens":1025,"completion_tokens":287,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":641,"completion_tokens_details":{"reasoning_tokens":195}},"tokens_in":641,"tokens_out":287,"duration_ms":2894,"temperature":1.0,"reasoning_tokens":195,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T22:49:33.966450+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the same f(Q,T) field equations but use the standard dust trace T = −ρ, redo the Markov-chain Monte Carlo fit, and compare the best-fit ranges to Table I; if the ranges shift by more than the quoted 1-σ errors, the reconstruction as stated does not describe the theory. A direct check is to insert the fitted γ(1), γ(3), γ(4), η(1), η(2), and Ωm0 into Eqs. (18)–(19) and ask whether the resulting H(z) reproduces the Pantheon+SH0ES distance moduli within the published covariance.","supporting_citations":[{"cited_title":"In the literature, a plethora of works have been carried out on the astronomical and cosmological implications of this modified gravity [32–39]","cited_arxiv_id":null,"evidence_quote":"It introduces f(Q,T) gravity, whose action, field equations, and Friedmann equations form the basis of the reconstruction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It establishes Chebyshev-polynomial cosmography for one variable and the z-to-Chebyshev conversion that the paper extends to two variables."},{"cited_title":"On the Chebyshev approximation of a function with two variables","cited_arxiv_id":"1504.04693","evidence_quote":"It provides the two-variable Chebyshev series and coefficient formulas used to expand f(Q,T)."}],"review_version":1}