REVIEW
The four-loop quark-to-quark splitting function in QCD
High Energy Physics - Phenomenology · Submitted 29 Sep 2026 · 67 references
Pith's verdict: The all-n four-loop pure-singlet splitting function is real and independently cross-checked; the only real gap is the paper never shows the moment count and ansatz dimension behind the harmonic-sum reconstruction.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
The paper computes the four-loop pure-singlet quark-to-quark splitting function of QCD analytically for all momentum fractions, completing the quark-to-quark entry of the N3LO evolution matrix and removing its approximation uncertainty.
The reading
What carries the argument
The load-bearing object is the four-loop off-shell quark operator matrix element $\langle q|O^B_q|q\rangle$ for the flavour-singlet twist-two operator $O_q(n)$, computed in dimensional regularization; its ultraviolet pole, after subtracting operator mixing with the gluon, gauge-variant operators, and lower-order counterterms, fixes the anomalous dimension $\gamma^{(3)}_{qq}$, and subtracting the known non-singlet part yields $\gamma^{(3)}_{\rm ps}$. Two mechanisms carry the argument: an auxiliary parameter $t$ converts the spin-$n$ powers in the operator into linear propagators so that the expansion in $t$ produces the individual Mellin moments, and a reconstruction ansatz — harmonic sums up to weight six with rational prefactors restricted to denominators $n$, $n\pm1$, $n\pm2$ — turns a finite set of exactly computed moments into a closed all-$n$ expression. The inverse Mellin transform into harmonic polylogarithms is what converts that expression into the momentum-fraction space where the splitting function and its small-$x$ limit are actually used.
Extended reading notes
Core claim
The central claim is that $\gamma^{(3)}_{\rm ps}(n)$, the Mellin-space anomalous dimension of the pure-singlet quark operator at four loops, is now known exactly for every integer moment $n$, and hence that $P^{(3)}_{\rm ps}(x)$ is known analytically for all momentum fractions. The result is obtained by computing the four-loop off-shell two-quark operator matrix element in dimensional regularization for a set of Mellin moments, subtracting the known non-singlet contribution, and reconstructing the all-$n$ dependence within the space of harmonic sums up to weight six with rational prefactors whose denominators involve only $n$, $n\pm1$, and $n\pm2$; an inverse Mellin transform then gives $P^{(3)}_{\rm ps}(x)$ in terms of harmonic polylogarithms and values of $\zeta$ up to weight six. The authors state that this expression agrees with all previously available partial results, including the moments through $n=22$, the analytic $\zeta_3$, $\zeta_4$, $\zeta_5$ coefficients, and the special $n=2$ structure in which the rational part and $S_{-2}(n)$ conspire to generate $\zeta_3$ contributions. The one reported disagreement is with the small-$x$ resummation of [18] at subleading colour: its $\zeta_2$ coefficient at $\log^4(x)$ is missing the $(d^{abc}_F)^2$ colour structure, a discrepancy whose structural origin the authors leave open.
Load-bearing premise
The all-$n$ formula is not established diagram by diagram but inferred by matching a finite set of exactly computed Mellin moments to an ansatz built from harmonic sums of weight at most six with rational prefactors whose only denominators are $n$, $n\pm1$, and $n\pm2$, so if the true four-loop anomalous dimension contained a structure outside that space the reconstructed expression could match every computed moment and still fail at other values of $n$.
Editorial extensions
If this is right
- Parton evolution at N3LO can now use the exact four-loop quark-to-quark splitting function at every momentum fraction instead of an approximation, removing the approximation uncertainty in this entry of the evolution matrix.
- The all-$n$ result confirms the fixed moments through $n=22$, the partial analytic $\zeta_3/\zeta_4/\zeta_5$ coefficients, and the special $n=2$ structure where the rational part and $S_{-2}(n)$ combine to produce $\zeta_3$ contributions.
- The exact $\gamma^{(3)}_{qq}$ supplies the quark contribution to the singlet trace $\gamma_{qq}+\gamma_{gg}$, whose reciprocity structure can help reconstruct the still-unknown four-loop gluon anomalous dimension $\gamma^{(3)}_{gg}$ from fixed moments.
- The small-$x$ expansion fixes the previously unknown next-to-leading-power coefficients $m^{(3)}_{-1,1}$ and $m^{(3)}_{-1,0}$, and the accompanying fit ($10^{-5}$ absolute, $10^{-8}$ relative) makes the result directly usable in evolution codes.
- The one reported disagreement — the missing $(d^{abc}_F)^2$ structure in the resummed $\zeta_2\log^4(x)$ coefficient of [18] — means resummation-based small-$x$ estimates for this channel are incomplete at subleading colour.
Reading between the lines
- If the reconstruction assumption is sound, the same moment-matching ansatz should yield exact all-$n$ forms for the remaining four-loop singlet entries ($P_{qg}$, $P_{gq}$, $P_{gg}$), and the success here is weak evidence that those quantities live in the same weight-six harmonic-sum space; applying the method to them is the natural next step.
- The missing $(d^{abc}_F)^2$ term in the resummation of [18] raises the possibility that other subleading-colour small-$x$ coefficients of the singlet evolution matrix carry analogous omissions; now that one entry is exact, those coefficients can be benchmarked against it.
- The conspicuous $n=2$ cancellation pattern (rational term $3/4$ combining with $S_{-2}(n)/(n-2)$ to produce $7/4 - (3/2)\zeta_3$) looks like the fingerprint of a structural constraint; deriving that constraint could supply the deeper structural explanation the authors call for and might predict which subleading-colour terms a corrected resummation must contain.
Original abstract
The scale evolution of parton distributions is governed by splitting functions. We compute the four-loop pure-singlet contribution to the quark-to-quark splitting function in perturbative QCD. Together with the known non-singlet contribution, our result completes the quark-to-quark splitting function at this order. We confirm previous partial results and obtain, for the first time, a fully analytic expression valid for all momentum fractions. We examine its behavior at small momentum fractions and provide precise numerical representations suitable for parton evolution.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
The scale evolution of parton distributions is governed by splitting functions. We compute the four-loop pure-singlet contribution to the quark-to-quark splitting function in perturbative QCD. Together with the known non-singlet contribution, our result completes the quark-to-quark splitting function at this order. We confirm previous partial results and obtain, for the first time, a fully analytic expression valid for all momentum fractions. We examine its behavior at small momentum fractions and provide precise numerical representations suitable for parton evolution.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | The scale evolution of parton distributions is governed by splitting functions. | background 1.00 | 0.03 | 0.02 | asserted 1.00 |
| s1 | We compute the four-loop pure-singlet contribution to the quark-to-quark splitting function in perturbative QCD. | contribution 0.99 | 0.97 | 0.35 | asserted 1.00 |
| s2 | Together with the known non-singlet contribution, our result completes the quark-to-quark splitting function at this order. | result 0.57 | 0.91 | 0.47 | asserted 1.00 |
| s3 | We confirm previous partial results and obtain, for the first time, a fully analytic expression valid for all momentum fractions. | result 0.81 | 0.96 | 0.44 | asserted 1.00 |
| s4 | We examine its behavior at small momentum fractions and provide precise numerical representations suitable for parton evolution. | contribution 0.70 | 0.95 | 0.09 | asserted 1.00 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Figures reproduced from arXiv:2609.37563 by the authors.
References
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.