In 1974 Georgi and Politzer showed that an asymptotically free non-Abelian gauge theory of the strong interactions - now QCD - reproduces Bjorken scaling up to computable logarithmic corrections: the moments of the deep-inelastic electroproduction structure functions are suppressed by inverse powers of ln⁡Q2\ln Q^2 for every spin-nn moment with n≠2n\neq 2, only the energy-momentum-tensor (n=2n=2) moment survives as Q2→∞Q^2\to\infty, and the limiting structure functions are singular distributions, 2xF1(x)=F2(x)=a δ(x)2xF_1(x)=F_2(x)=a\,\delta(x), with aa fixed by group theory (a=4/25=0.16a=4/25=0.16 for SU(3) with three quark triplets, the decimal value quoted in the paper's abstract). The paper itself acknowledged what this perturbative analysis cannot deliver: an unknown multiplicative constant from the strong-coupling end of the renormalization-group integral; unknown hadronic matrix elements multiplying every operator with nonzero anomalous dimension; an unsettled mass-generation mechanism for the gauge bosons (whether the Higgs mechanism could give them mass without destroying ultraviolet stability); lowest-order-only results with unexamined gauge and renormalization-scheme questions; and an unknown physical realization of the singular small-xx behavior. Later literature closed the purely perturbative items - anomalous dimensions and coefficient functions are now known to three loops, the scheme and gauge questions are understood, and the Higgs-mechanism question was answered and then mooted by confinement - while factorization theorems recast the nonperturbative items into the modern parton-distribution-function framework. What remains open is exactly the completion the 1974 analysis lacked: determining from QCD itself the nonperturbative inputs - the xx-dependent parton distributions, the leading higher-twist (twist-4) matrix elements, and the small-xx saturation dynamics that physically replaces the formal δ(x)\delta(x) limit - so that the structure functions at finite Q2Q^2 are predicted from first principles rather than fitted to the data they are meant to describe.

Background

Deep-inelastic scattering (DIS) is the scattering of a lepton, such as an electron, off a hadron target with large squared four-momentum transfer q2<0q^2<0 (one writes Q2=−q2Q^2=-q^2). In terms of the Bjorken variable x=Q2/(2p⋅q)x=Q^2/(2p\cdot q), with pp the target hadron's four-momentum - xx is the fraction of the target's longitudinal momentum carried, in the infinite-momentum frame, by the struck constituent - the inclusive cross section is parametrized by two structure functions, F1(x,Q2)F_1(x,Q^2) and F2(x,Q2)F_2(x,Q^2). Experiments at SLAC in the late 1960s found approximate Bjorken scaling, i.e. F1F_1 and F2F_2 depending nearly on xx alone, suggesting pointlike constituents (partons) inside the hadron.

Wilson's operator product expansion (OPE) - the short-distance expansion of a product of two local currents into a series of local operators with coefficient functions - gives the field-theoretic handle on DIS. Christ, Hasslacher, and Mueller (1972) showed that moments of the structure functions (∫01dx xn−2F2\int_0^1 dx\, x^{n-2}F_2 and analogous integrals) equal perturbatively calculable coefficient functions times hadronic matrix elements of local operators, with Q2Q^2-dependence governed by the anomalous dimensions of those operators - the deviations from canonical scaling produced by renormalization, computable from the renormalization-group (Callan-Symanzik) equations. The operators dominating the Bjorken limit are those of lowest twist (twist = mass dimension minus spin); the twist-2 operators are the symmetric tensors built from quark bilinears and from gluon field strengths.

In 1973, Gross and Wilczek and, independently, Politzer discovered asymptotic freedom: in a non-Abelian gauge theory the running coupling decreases logarithmically at short distances, so perturbation theory becomes increasingly accurate in the ultraviolet. Georgi and Politzer (1974), whose paper is the source of the present problem, applied the Christ-Hasslacher-Mueller technique to electroproduction in such a theory. Their one-loop results: (i) the anomalous-dimension matrix of the twist-2 operators has eigenvalues that are positive and increasing with nn, with a single zero eigenvalue at n=2n=2 corresponding to the conserved energy-momentum tensor; (ii) hence every moment with n≠2n\neq 2 vanishes as an inverse power of ln⁡Q2\ln Q^2 as Q2→∞Q^2\to\infty, so scaling holds up to calculable logarithms; (iii) only the n=2n=2 moment survives, so the limiting structure functions are singular distributions, 2xF1(x)=F2(x)=a δ(x)2xF_1(x)=F_2(x)=a\,\delta(x), with aa fixed by group-theoretic constants and the fermion charges (a=4/25=0.16a=4/25=0.16 for SU(3) with three quark triplets of charges 2/32/3, −1/3-1/3, −1/3-1/3, the decimal value quoted in the paper's abstract); and (iv) for n=2n=2 the required matrix element is that of the energy-momentum tensor, so the Q2→∞Q^2\to\infty limit is completely determined. The theory in question is what is now called QCD, the SU(3) gauge theory of quarks and gluons.

The same paper stated the limits of the calculation. It is lowest nontrivial order in the coupling gg, and the gauge-independence and coupling-choice arguments invoked hold only at that order. The vector bosons (gluons) must acquire mass somehow, and it was 'not known' whether scalar (Higgs) mass generation could do this without destroying the ultraviolet stability (asymptotic freedom) of the theory, so a nonperturbative mass mechanism was simply assumed, with fermion bare masses set to zero. The perturbative evaluation of the renormalization-group integral ∫0tdt′ γn(gˉ(g,t′))\int_0^t dt'\,\gamma_n(\bar g(g,t')) crosses a strong-coupling region (the neighborhood of the reference scale, the 'small-tt' region' in the paper's words), which introduces an unknown multiplicative constant into the asymptotic coefficient functions - a defect the authors set aside only because those terms are anyway multiplied by unknown matrix elements. The limiting structure functions exist only 'in the sense of distribution theory', so 'true scaling is never realized', and the physical realization of the implied singular small-xx behavior in a realistic theory was left unsettled. Finally, the matrix elements multiplying operators with nonzero anomalous dimensions are unknown nonperturbative inputs, so a complete quantitative prediction of the finite-Q2Q^2 structure functions requires nonperturbative information.

Much of the subsequent program is complete. The moment language was recast as evolution equations for parton distribution functions (PDFs - densities of partons carrying momentum fraction xx; Altarelli and Parisi, 1977) - the DGLAP equations, named after Dokshitzer, Gribov, Lipatov, Altarelli, and Parisi - whose perturbative Q2Q^2-evolution is governed by the same anomalous dimensions, now called splitting functions: computed to two loops by Floratos, Ross, and Sachrajda (1977), completed with the full next-to-leading-order coefficient functions and scheme analysis by Bardeen, Buras, Duke, and Muta (1978) and by Curci, Furmanski, and Petronzio (1980), reviewed by Buras (1980), and extended to three loops by Moch, Vermaseren, and Vogt (2004). Factorization theorems (Collins, Soper, and Sterman) make the separation of short-distance coefficients from long-distance matrix elements rigorous, and the renormalon analysis (Beneke, 1999; renormalons are the intrinsic ambiguities of the divergent perturbative series) shows how the residual nonperturbative ambiguity of the perturbative series reappears as power corrections - precisely the higher-twist contributions of order mp2/Q2m_p^2/Q^2 (mpm_p is the proton mass) that the 1974 analysis dropped (their operator structure was worked out by Jaffe and Soldate, 1982). In parton language the Georgi-Politzer result means that evolution drives the distributions toward small xx as Q2Q^2 grows; the quantitative small-xx regime is instead governed by the BFKL resummation of powers of αsln⁡(1/x)\alpha_s\ln(1/x) (with αs\alpha_s the strong coupling; Kuraev, Lipatov, and Fadin, 1977; Balitsky and Lipatov, 1978) and is tamed, at sufficiently small xx, by nonlinear saturation dynamics (Gribov, Levin, and Ryskin, 1983; and the color-glass-condensate framework, the effective theory of the saturated high-density gluon fields in fast hadrons and nuclei, reviewed by Gelis, Iancu, Jalilian-Marian, and Venugopalan, 2010). Data from HERA, the electron-proton collider operated at DESY (the combined H1 and ZEUS analysis), confirm the strong rise of F2F_2 at small xx. What has never been achieved is the first-principles determination of the nonperturbative inputs themselves: PDFs are extracted from experimental data by global fits, the twist-4 matrix elements are only partly constrained, and the lattice-QCD program (lattice QCD computes the theory by Monte Carlo evaluation of its discretized Euclidean space-time path integral) for computing xx-dependent PDFs from Euclidean correlators (initiated by Ji, 2013) is still maturing (Constantinou and Cichy, 2026).

Problem Statement

The open problem is the nonperturbative completion of the 1974 Georgi-Politzer electroproduction program. The paper's own analysis leaves the finite-q2q^2 structure functions undetermined: its footnote to the renormalization-group solution records that all terms with nonzero anomalous dimension enter the forward Compton amplitude 'multiplied by unknown matrix elements', its moment relations were derived 'Ignoring terms of order mp2/q2m_p^2/q^2', and the limiting structure functions exist only 'in the sense of distribution theory', so that 'true scaling is never realized'. Determine, from QCD alone:

(a) the hadronic matrix elements of the leading-twist (twist-2) operators entering the operator-product-expansion moments of the electroproduction structure functions - equivalently, the xx-dependent parton distribution functions together with their normalization, which fixes the overall constant that the 1974 renormalization-group solution could not determine;

(b) the leading higher-twist (twist-4, order mp2/Q2m_p^2/Q^2) contributions to the structure functions and their moments at finite Q2Q^2;

(c) the physical behavior at small Bjorken xx toward which the evolution drives the structure functions - the quantitative form of the regime in which BFKL-type growth is tamed by saturation, which is what physically replaces the formal distributional limit F2→a δ(x)F_2\to a\,\delta(x) in a realistic theory with nonzero quark masses and nonperturbative dynamics;

and thereby predict the electroproduction structure functions F1(x,Q2)F_1(x,Q^2) and F2(x,Q2)F_2(x,Q^2) (or equivalently their moments) at finite Q2Q^2 from first principles rather than from inputs fitted to the data being described. Items (a), (b), and (c) are jointly required components of this completion, not independent alternatives: the completion combines them into a prediction of F1F_1 and F2F_2 (or equivalently their moments) at finite Q2Q^2, with the nonperturbative inputs derived from QCD itself rather than fitted to the electroproduction data being described. The purely perturbative items of the 1974 list (higher-order anomalous dimensions and coefficient functions, gauge- and scheme-dependence questions) and the formal Higgs-mechanism/ultraviolet-stability question are already closed by later literature and are not part of this problem; the mass-generation mechanism itself is settled by confinement in QCD.

References

  1. H. Georgi and H. D. Politzer, 'Electroproduction scaling in an asymptotically free theory of strong interactions', Phys. Rev. D 9, 416 (1974), DOI: 10.1103/PhysRevD.9.416, https://doi.org/10.1103/PhysRevD.9.416
  2. K. G. Wilson, 'Non-Lagrangian Models of Current Algebra', Phys. Rev. 179, 1499 (1969), DOI: 10.1103/PhysRev.179.1499, https://doi.org/10.1103/PhysRev.179.1499
  3. N. Christ, B. Hasslacher, and A. H. Mueller, 'Light-Cone Behavior of Perturbation Theory', Phys. Rev. D 6, 3543 (1972), DOI: 10.1103/PhysRevD.6.3543, https://doi.org/10.1103/PhysRevD.6.3543
  4. D. J. Gross and F. Wilczek, 'Ultraviolet Behavior of Non-Abelian Gauge Theories', Phys. Rev. Lett. 30, 1343 (1973), DOI: 10.1103/PhysRevLett.30.1343, https://doi.org/10.1103/PhysRevLett.30.1343
  5. H. D. Politzer, 'Reliable perturbative results for strong interactions?', Phys. Rev. Lett. 30, 1346 (1973), DOI: 10.1103/PhysRevLett.30.1346, https://doi.org/10.1103/PhysRevLett.30.1346
  6. D. J. Gross and F. Wilczek, 'Asymptotically free gauge theories. I', Phys. Rev. D 8, 3633 (1973), DOI: 10.1103/PhysRevD.8.3633, https://doi.org/10.1103/PhysRevD.8.3633
  7. E. S. Fradkin and O. K. Kalashnikov, 'Asymptotically free SU(3) model of strong interactions', Phys. Lett. B 59, 159 (1975), DOI: 10.1016/0370-2693(75)90692-9, https://doi.org/10.1016/0370-2693(75)90692-9
  8. A. Salam and J. A. Strathdee, 'Higgs couplings and asymptotic freedom', Phys. Rev. D 18, 4713 (1978), DOI: 10.1103/PhysRevD.18.4713, https://doi.org/10.1103/PhysRevD.18.4713
  9. E. H. Fradkin and S. H. Shenker, 'Phase diagrams of lattice gauge theories with Higgs fields', Phys. Rev. D 19, 3682 (1979), DOI: 10.1103/PhysRevD.19.3682, https://doi.org/10.1103/PhysRevD.19.3682
  10. G. Altarelli and G. Parisi, 'Asymptotic freedom in parton language', Nucl. Phys. B 126, 298 (1977), DOI: 10.1016/0550-3213(77)90384-4, https://doi.org/10.1016/0550-3213(77)90384-4
  11. E. G. Floratos, D. A. Ross, and C. T. Sachrajda, 'Higher-order effects in asymptotically free gauge theories: the anomalous dimensions of Wilson operators', Nucl. Phys. B 129, 66 (1977) [Erratum: Nucl. Phys. B 139, 545 (1978)], DOI: 10.1016/0550-3213(77)90020-7, https://doi.org/10.1016/0550-3213(77)90020-7
  12. W. A. Bardeen, A. J. Buras, D. W. Duke, and T. Muta, 'Deep-inelastic scattering beyond the leading order in asymptotically free gauge theories', Phys. Rev. D 18, 3998 (1978), DOI: 10.1103/PhysRevD.18.3998, https://doi.org/10.1103/PhysRevD.18.3998
  13. G. Curci, W. Furmanski, and R. Petronzio, 'Evolution of parton densities beyond leading order: the nonsinglet case', Nucl. Phys. B 175, 27 (1980), DOI: 10.1016/0550-3213(80)90003-6, https://doi.org/10.1016/0550-3213(80)90003-6
  14. A. J. Buras, 'Asymptotic freedom in deep inelastic processes in the leading order and beyond', Rev. Mod. Phys. 52, 199 (1980), DOI: 10.1103/RevModPhys.52.199, https://doi.org/10.1103/RevModPhys.52.199
  15. S. Moch, J. A. M. Vermaseren, and A. Vogt, 'The three-loop splitting functions in QCD: the non-singlet case', Nucl. Phys. B 688, 101 (2004), arXiv:hep-ph/0403192, DOI: 10.1016/j.nuclphysb.2004.03.030, https://arxiv.org/abs/hep-ph/0403192
  16. A. Vogt, S. Moch, and J. A. M. Vermaseren, 'The three-loop splitting functions in QCD: the singlet case', Nucl. Phys. B 691, 129 (2004), arXiv:hep-ph/0404111, DOI: 10.1016/j.nuclphysb.2004.04.024, https://arxiv.org/abs/hep-ph/0404111
  17. J. C. Collins, D. E. Soper, and G. Sterman, 'Factorization of hard processes in QCD', Adv. Ser. Direct. High Energy Phys. 5, 1 (1989), arXiv:hep-ph/0409313, DOI: 10.1142/9789814503266_0001, https://arxiv.org/abs/hep-ph/0409313
  18. M. Beneke, 'Renormalons', Phys. Rep. 317, 1 (1999), arXiv:hep-ph/9807443, DOI: 10.1016/S0370-1573(98)00130-6, https://arxiv.org/abs/hep-ph/9807443
  19. R. L. Jaffe and M. Soldate, 'Twist-4 effects in electroproduction: Canonical operators and coefficient functions', Phys. Rev. D 26, 49 (1982), DOI: 10.1103/PhysRevD.26.49, https://doi.org/10.1103/PhysRevD.26.49
  20. E. A. Kuraev, L. N. Lipatov, and V. S. Fadin, 'The Pomeranchuk singularity in nonabelian gauge theories', Sov. Phys. JETP 45, 199 (1977) [Zh. Eksp. Teor. Fiz. 72, 377 (1977)], INSPIRE literature record 123995 (no DOI exists for this pre-digital journal), https://inspirehep.net/literature/123995
  21. I. I. Balitsky and L. N. Lipatov, 'The Pomeranchuk singularity in quantum chromodynamics', Sov. J. Nucl. Phys. 28, 822 (1978) [Yad. Fiz. 28, 1597 (1978)], INSPIRE literature record 137229 (no DOI exists for this pre-digital journal), https://inspirehep.net/literature/137229
  22. L. V. Gribov, E. M. Levin, and M. G. Ryskin, 'Semihard processes in QCD', Phys. Rep. 100, 1 (1983), DOI: 10.1016/0370-1573(83)90022-4, https://doi.org/10.1016/0370-1573(83)90022-4
  23. F. Gelis, E. Iancu, J. Jalilian-Marian, and R. Venugopalan, 'The color glass condensate', Annu. Rev. Nucl. Part. Sci. 60, 463 (2010), arXiv:1002.0333, DOI: 10.1146/annurev.nucl.010909.083629, https://arxiv.org/abs/1002.0333
  24. H1 and ZEUS Collaborations, 'Combination of measurements of inclusive deep inelastic e±p scattering cross sections and QCD analysis of HERA data', Eur. Phys. J. C 75, 580 (2015), arXiv:1506.06042, DOI: 10.1140/epjc/s10052-015-3710-4, https://arxiv.org/abs/1506.06042
  25. X. Ji, 'Parton physics on a Euclidean lattice', Phys. Rev. Lett. 110, 262002 (2013), arXiv:1305.1539, DOI: 10.1103/PhysRevLett.110.262002, https://arxiv.org/abs/1305.1539
  26. M. Constantinou and K. Cichy, 'Parton distribution functions from lattice QCD', arXiv:2608.00251 (2026), https://arxiv.org/abs/2608.00251