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 for every spin- moment with , only the energy-momentum-tensor () moment survives as , and the limiting structure functions are singular distributions, , with fixed by group theory ( 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- 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 -dependent parton distributions, the leading higher-twist (twist-4) matrix elements, and the small- saturation dynamics that physically replaces the formal limit - so that the structure functions at finite 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 (one writes ). In terms of the Bjorken variable , with the target hadron's four-momentum - 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, and . Experiments at SLAC in the late 1960s found approximate Bjorken scaling, i.e. and depending nearly on 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 ( and analogous integrals) equal perturbatively calculable coefficient functions times hadronic matrix elements of local operators, with -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 , with a single zero eigenvalue at corresponding to the conserved energy-momentum tensor; (ii) hence every moment with vanishes as an inverse power of as , so scaling holds up to calculable logarithms; (iii) only the moment survives, so the limiting structure functions are singular distributions, , with fixed by group-theoretic constants and the fermion charges ( for SU(3) with three quark triplets of charges , , , the decimal value quoted in the paper's abstract); and (iv) for the required matrix element is that of the energy-momentum tensor, so the 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 , 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 crosses a strong-coupling region (the neighborhood of the reference scale, the 'small-' 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- 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- 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 ; Altarelli and Parisi, 1977) - the DGLAP equations, named after Dokshitzer, Gribov, Lipatov, Altarelli, and Parisi - whose perturbative -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 ( 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 as grows; the quantitative small- regime is instead governed by the BFKL resummation of powers of (with the strong coupling; Kuraev, Lipatov, and Fadin, 1977; Balitsky and Lipatov, 1978) and is tamed, at sufficiently small , 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 at small . 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 -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- 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 ', 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 -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 ) contributions to the structure functions and their moments at finite ;
(c) the physical behavior at small Bjorken 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 in a realistic theory with nonzero quark masses and nonperturbative dynamics;
and thereby predict the electroproduction structure functions and (or equivalently their moments) at finite 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 and (or equivalently their moments) at finite , 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- 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
- M. Constantinou and K. Cichy, 'Parton distribution functions from lattice QCD', arXiv:2608.00251 (2026), https://arxiv.org/abs/2608.00251