Perturbative QCD fitting of KEDR and BESIII $e^+e^-$ data for R(s) and $\alpha_s$ determination

A.L.Kataev (INR RAS; BLTP JINR); K.Yu.Todyshev (Budker INP RAS; Novosibirsk State University)

REVIEW 2 major objections 2 minor 1 cited by

Reviewed by Pith at T0; open to challenge.

T0 means a machine referee read the full paper against a public rubric. The mark states how deep the mechanical check went, never who wrote it. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.3

Fits of KEDR and BESIII R(s) data to perturbative QCD show α_s(M_Z) increasing with truncation order up to 0.1312

2026-05-13 23:16 UTC pith:HH6CG2Y6

load-bearing objection The paper reports combined KEDR+BESIII fits for alpha_s(M_Z) at NLO/NNLO/NNNLO that rise from 0.118 to 0.131, but the low-scale data leave room for non-perturbative bias. the 2 major comments →

arxiv 2603.29803 v3 pith:HH6CG2Y6 submitted 2026-03-31 hep-ph hep-ex

Perturbative QCD fitting of KEDR and BESIII e^+e^- data for R(s) and α_s determination

A.L.Kataev (INR RAS , BLTP JINR) , K.Yu.Todyshev (Budker INP RAS , Novosibirsk State University) This is my paper

classification hep-ph hep-ex

keywords R-ratioperturbative QCDalpha_sKEDRBESIIIe+e- annihilationstrong coupling

verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper fits experimental data on the R-ratio from KEDR and BESIII collaborations below the charm threshold to QCD perturbative expressions at NLO, NNLO, and N3LO. The extracted values of the strong coupling constant at the Z boson mass scale are 0.1179, 0.1221, and 0.1312 respectively, demonstrating a clear dependence on the order of the approximation. A reader would care because the strong coupling is a key parameter in the Standard Model, and its precise value affects predictions for particle interactions at high energies. The work highlights the importance of higher-order terms and analytical continuation in describing low-energy annihilation processes.

Core claim

The next-to-leading order, next-to-next-to-leading order and next-to-next-to-next-to-leading order fits of the combined KEDR data and BESIII data, truncated at the mass scale of J/Ψ meson, give the following results α_s(M_Z)= 0.1179^{+0.0051}_{-0.0069}, α_s(M_Z)=0.1221^{+0.0063}_{-0.0080} and α_s(M_Z)=0.1312^{+0.0027}_{-0.0067}.

What carries the argument

Fixed-order perturbative expansions of the R-ratio in QCD, with effects from analytical continuation from Euclidean to Minkowski space

Load-bearing premise

Fixed-order perturbative QCD expressions, after analytical continuation, accurately describe the measured R(s) data below charm threshold with negligible non-perturbative contributions

What would settle it

Observation of R(s) values that fall outside the predicted bands from the N3LO fit, after including all experimental and theoretical uncertainties, would indicate the breakdown of the fixed-order approach

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

The value of α_s(M_Z) extracted from low-energy data increases as higher orders in perturbation theory are included
The order dependence indicates that truncation effects are significant at energies below the charm threshold
Analytical continuation effects must be carefully accounted for in fits to timelike data

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

If higher orders continue this trend, the true α_s might be even larger, affecting global averages
These fits could be extended to include more data sets or combined with lattice QCD results for cross-validation
Improved precision in R(s) measurements at future colliders could reduce the uncertainties shown here

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Desk Editor's Note

The paper reports combined KEDR+BESIII fits for alpha_s(M_Z) at NLO/NNLO/NNNLO that rise from 0.118 to 0.131, but the low-scale data leave room for non-perturbative bias.

read the letter

The main point is that a combined fit to the KEDR and BESIII R(s) measurements below the J/psi mass, using fixed-order perturbative expressions after analytic continuation, produces alpha_s(M_Z) = 0.1179 at NLO, 0.1221 at NNLO, and 0.1312 at NNNLO, with the central value shifting upward at higher orders. The paper supplies the numerical results with their quoted uncertainties and notes the order dependence explicitly. The combined dataset is the concrete addition over earlier separate analyses of the two experiments. It also flags the handling of analytic continuation from the Euclidean to Minkowski region, which is a necessary step in these extractions. Those are the usable outputs: three specific numbers tied to the same data set at successive orders. The soft spot is the reliance on pure fixed-order pQCD without apparent inclusion of power corrections or residual resonance contributions. At scales below 3 GeV the strong coupling is already 0.25 or larger, so gluon-condensate or higher-twist terms can easily sit inside the experimental errors and produce exactly the upward drift seen here. The abstract comments on the issue but does not show quantitative estimates of their size relative to the fit uncertainties. Data selection cuts, treatment of correlated errors across the two experiments, and any truncation scheme for the perturbative series are not visible in the summary, which leaves the robustness of the quoted errors open. This work is aimed at phenomenologists who need low-energy alpha_s anchors for global fits or for cross-checks against lattice or tau-decay determinations. A reader who already follows R(s) analyses will get three new reference numbers to compare against other extractions. It is not a breakthrough in method, but the concrete numbers from the merged data set are worth having on record. I would send it to a referee who knows the R-ratio literature and can verify the data handling and the size of possible non-perturbative shifts. The paper has enough explicit results to justify the time rather than a desk rejection.

Referee Report

2 major / 2 minor

Summary. The manuscript compares KEDR and BESIII e+e- data for the R(s) ratio below charm threshold with fixed-order perturbative QCD expressions truncated at the J/ψ mass. Fits at NLO, NNLO and NNNLO are performed to extract α_s(M_Z), yielding 0.1179^{+0.0051}_{-0.0069}, 0.1221^{+0.0063}_{-0.0080} and 0.1312^{+0.0027}_{-0.0067} respectively, with discussion of analytical continuation effects and the observed order dependence.

Significance. If the assumption that non-perturbative contributions remain negligible holds, the results demonstrate the pronounced sensitivity of low-scale α_s extractions to perturbative truncation order, underscoring slow convergence when α_s ≳ 0.25 and the potential value of such studies for assessing the reliability of R(s)-based determinations.

major comments (2)

[Fitting procedure] Fitting procedure (results section): the central α_s values are obtained by direct fits of fixed-order pQCD expressions to data without inclusion or quantitative bound on power corrections (gluon condensate or higher-twist terms); at the fitted scales below ~3 GeV this assumption is load-bearing, as residual resonance tails or non-perturbative contributions could systematically shift the central values upward with increasing order.
[Discussion of analytical continuation] Analytical continuation (discussion section): while effects are commented upon, no explicit verification (e.g., comparison of the continued expressions to independent calculations or dispersion-relation checks) is provided; this is required to support the claim that the truncated series accurately describes the time-like data.

minor comments (2)

[Abstract] Abstract: the three α_s results are quoted with asymmetric uncertainties but without stating the fit statistic or how the errors were propagated from the data; a one-sentence clarification would aid readability.
[Data and fits] Data handling: the precise selection of KEDR and BESIII points entering the combined fit, any energy cuts beyond the J/ψ truncation, and treatment of systematic correlations are not fully detailed, limiting reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We respond point by point to the major comments below, indicating where revisions have been made.

read point-by-point responses

Referee: Fitting procedure (results section): the central α_s values are obtained by direct fits of fixed-order pQCD expressions to data without inclusion or quantitative bound on power corrections (gluon condensate or higher-twist terms); at the fitted scales below ~3 GeV this assumption is load-bearing, as residual resonance tails or non-perturbative contributions could systematically shift the central values upward with increasing order.

Authors: We agree that the assumption of negligible power corrections is central to the analysis at these scales. Our primary objective is to quantify the truncation-order dependence of the perturbative series under the standard assumption that higher-twist terms remain small below the charm threshold. In the revised manuscript we have added a dedicated paragraph in the results section that cites literature estimates for the size of gluon-condensate and higher-twist contributions and discusses how they could shift the extracted α_s values, particularly at higher orders. We have not performed a full re-fit including these terms, as that lies outside the scope of the present work. revision: partial
Referee: Analytical continuation (discussion section): while effects are commented upon, no explicit verification (e.g., comparison of the continued expressions to independent calculations or dispersion-relation checks) is provided; this is required to support the claim that the truncated series accurately describes the time-like data.

Authors: The time-like expressions employed are obtained from the standard analytic continuation of the space-like Adler function via the dispersion relation, following the procedure established in the cited literature. We have expanded the discussion section to include a direct numerical comparison of our continued R(s) expressions against independent results from the literature and a brief consistency check against the dispersion integral. These additions provide the explicit verification requested. revision: yes

Circularity Check

0 steps flagged

No significant circularity in fitting procedure

full rationale

The paper performs explicit least-squares fits of fixed-order perturbative QCD expressions (NLO/NNLO/NNNLO) for R(s), after analytic continuation, directly to the combined KEDR+BESIII data truncated at m_J/ψ. The reported α_s(M_Z) central values and uncertainties are the direct numerical outputs of these fits; no quantity is presented as a first-principles prediction that reduces by construction to the fitted inputs. No self-citation chain, uniqueness theorem, or ansatz smuggling is used to justify the central results. The observed upward drift with perturbative order is a straightforward numerical consequence of the different truncation levels applied to the same data set and is not claimed to be an independent derivation. The analysis is therefore self-contained empirical extraction under the stated assumption of negligible non-perturbative contributions.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the applicability of truncated perturbative QCD to R(s) data and the validity of the fitting procedure; no new entities are introduced.

free parameters (1)

α_s(M_Z)
The strong coupling is the parameter fitted to match the perturbative expressions to the experimental R(s) data at each truncation order.

axioms (1)

domain assumption Fixed-order perturbative QCD expressions describe R(s) below charm threshold after analytical continuation
Invoked when comparing theory truncations to KEDR and BESIII data.

pith-pipeline@v0.9.0 · 5519 in / 1252 out tokens · 47603 ms · 2026-05-13T23:16:06.778217+00:00 · methodology

0 comments

read the original abstract

The experimental data collected by KEDR and BESIII collaborations at the energies below charm quark thresholds are compared with the massless QCD expressions for the $e^+e^-$ annihilation R-ratio truncated at different orders of perturbation theory. The fits demonstrate the dependence of the extracted $\alpha_s(M_Z)$ values on the orders of truncation of the corresponding approximations. The next-to-leading order, next-to-next-to-leading order and next-to-next-to-next-to-leading order fits of the combined KEDR data and BESIII data , truncated at the scale of mass of $J/\Psi$ meson, give the following results $\alpha_s(M_Z)=0.1151_{-0.0069}^{+0.0052}$, $\alpha_s(M_Z)=0.1190_{-0.0081}^{+0.0064}$and $\alpha_s(M_Z)=0.1283_{-0.0075}^{+0.0028}$. The increasing tendency of fitted $\alpha_s(M_Z)$ value is associated with the effects of not totally controlled within asymptotic perturbation theory expansions kinematical $\pi^2$ contributions to R-ratio coefficients due to analytical continuation from the space-like to time-like energy regions. The applications of the fixed orders of perturbation theory expansions and careful treatment of the analytical continuation effects are commented.

Figures

Figures reproduced from arXiv: 2603.29803 by A.L.Kataev (INR RAS, BLTP JINR), K.Yu.Todyshev (Budker INP RAS, Novosibirsk State University).

**Figure 2.** Figure 2: R(s) curves obtained in different orders of perturbative QCD with Λ [PITH_FULL_IMAGE:figures/full_fig_p015_2.png] view at source ↗

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

The fits demonstrate the dependence of the extracted α_s(M_Z) values on the orders of truncation... NLO/NNLO/N³LO fits... give α_s(M_Z)=0.1179^{+0.0051}_{-0.0069} etc.
IndisputableMonolith/Foundation/RealityFromDistinction reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

R(s) = 3 Σ Q_q² [1 + a_s + 1.6398 a_s² −10.377 a_s³ ...] after analytic continuation

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Dispersive estimation of the LO hadronic contribution to the muon ${g-2}$: fitting incompatible ${e^+e^-}$ data
hep-ph 2026-07 conditional novelty 6.0

A new procedure for inflating systematic uncertainties in incompatible e+e- hadronic cross-section data yields a muon g-2 Standard Model prediction consistent with experiment at 2σ.