In 1977 De Rújula, Georgi and Politzer showed that Bloom–Gilman local duality between the resonance region and the deep-inelastic scaling limit of electroproduction follows from the operator-product-expansion (OPE) moment predictions of QCD, and made this a quantitative test of asymptotic freedom at 1≲Q2≲161\lesssim Q^{2}\lesssim 16 GeV2^{2}, with a phenomenological fit of the proton's twist-two moments AnA_{n} and a parametrized, undetermined higher-twist sector BnkB_{nk}. The authors' own abstract acknowledged part of what the test still lacked: the existing R=σL/σTR=\sigma_{L}/\sigma_{T} data at large ξ\xi and Q2Q^{2} lay systematically above their twist-two prediction, the higher-twist sector BnkB_{nk} was parametrized, with the scale M0≈400M_{0}\approx 400 MeV obtained from resonance-region data, and the twist-two moments AnA_{n} came from a fit to scaling-region data; their companion paper ('Trouble with ξ\xi scaling?', Physical Review D 15 (1977) 2495) restated that the higher-twist corrections are negligible only for WW well above the proton mass. The subsequent literature identified the remaining gaps: a determined gluon distribution and controlled quark–gluon operator mixing below ξ≲0.3\xi\lesssim 0.3, higher-order perturbative input beyond the one-loop terms used in 1977, first-principles or independent determination of the nn-dependent twist-four and higher matrix elements Bnk(Q2)B_{nk}(Q^{2}), structure-function data at ξ≳0.8\xi\gtrsim 0.8, and a nonperturbative, first-principles determination of the proton matrix elements AnA_{n} and of the quark/gluon separation. Later literature has delivered the perturbative orders (next-to-next-to-leading order (NNLO) splitting functions, five-loop β\beta function), the phenomenological gluon distribution with full singlet mixing, precision separated structure functions and extensive duality verification (Jefferson Lab), and the first physical-continuum lattice determinations of the nucleon gluon distribution. What no later work supplies is the problem's open core: a first-principles determination of the nn-dependent higher-twist coefficients BnkB_{nk} and of the complete set of twist-two proton moments, enabling parameter-free local-duality predictions of the resonance-region structure functions, of RR at large ξ\xi, and of the elastic form factors through the threshold moment relations, tested against precision data.

Background

In inclusive electron–proton scattering, a virtual photon of four-momentum squared Q2Q^{2} and energy transfer ν\nu strikes the proton, and the cross section is written in terms of structure functions, conventionally F2(x,Q2)F_{2}(x,Q^{2}) (historically νW2\nu W_{2}) and F1(x,Q2)F_{1}(x,Q^{2}), where x=Q2/(2Mν)x=Q^{2}/(2M\nu) is the Bjorken variable and MM the proton mass. At large Q2Q^{2} these functions depend approximately only on xx ("scaling"), while asymptotic freedom — the decrease of the QCD coupling αs(Q2)\alpha_{s}(Q^{2}) at large Q2Q^{2} — predicts small logarithmic violations of scaling. The light-cone operator product expansion (OPE) organizes this: moments of the structure functions are sums of matrix elements of local quark and gluon operators of increasing "twist" (dimension minus spin), each multiplied by perturbatively computable coefficient functions. Twist-two operators give the leading logarithmic evolution (via their anomalous dimensions, which govern Q2Q^{2}-dependence), while twist-four and higher operators give power corrections suppressed by M02/Q2M_{0}^{2}/Q^{2} for some hadronic scale M0M_{0}, and enhanced by powers of the moment index nn. Because the proton mass is not negligible at moderate Q2Q^{2}, one works with the Nachtmann variable ξ=2x/(1+1+4x2M2/Q2)\xi = 2x/(1+\sqrt{1+4x^{2}M^{2}/Q^{2}}) (introduced by Nachtmann in 1973), whose moments ∫dξ ξnF\int d\xi\, \xi^{n}F project onto operators of definite spin and resum the kinematic target-mass corrections.

Two empirical phenomena in the moderate-Q2Q^{2} data motivated this framework. "Precocious scaling" is the observation that approximate scaling behavior sets in at unexpectedly low Q2Q^{2}. Bloom–Gilman duality (quark–hadron duality), observed by Bloom and Gilman in 1970–71, is the fact that the resonance-region structure function, averaged over an appropriate energy interval, approximately equals the smooth deep-inelastic scaling curve at the same ξ\xi; it is called "local" duality when the equality holds resonance by resonance rather than only after averaging over the whole resonance region. The ratio R=σL/σTR=\sigma_{L}/\sigma_{T} of longitudinal to transverse virtual-photon absorption cross sections isolates the longitudinal structure function FLF_{L}, which leading-twist (twist-two) QCD predicts to be nonzero but small at large Q2Q^{2}.

De Rújula, Georgi and Politzer ("DGP", Annals of Physics 103, 315 (1977)) gave this its QCD "demythification": they showed that local duality is a consequence of the OPE moment predictions, fitted the leading-twist evolution (one-loop anomalous dimensions and coefficient functions) to SLAC (Stanford Linear Accelerator Center) electroproduction data for 1≲Q2≲161\lesssim Q^{2}\lesssim 16 GeV2^{2}, used resonance-region moments to extract the higher-twist scale M0≈400M_{0}\approx 400 MeV from terms of the form (nM02/Q2)kBnk(Q2)(nM_{0}^{2}/Q^{2})^{k}B_{nk}(Q^{2}), predicted R=σL/σTR=\sigma_{L}/\sigma_{T}, and derived "threshold duality" relations expressing the elastic electromagnetic form factors GEG_{E}, GMG_{M} as integrals of the scaling structure function over ξ\xi between the pion threshold and ξ=1\xi=1, from which they extracted the proton magnetic form factor GMpG_{M}^{p} from large-ξ\xi resonance data. The whole construction, however, rested on inputs left undetermined in the analysis — some flagged by the authors themselves, others identified by the subsequent literature: the coefficients BnkB_{nk} were parametrized rather than computed; gluons were neglected below ξ≲0.3\xi\lesssim 0.3 and quark–gluon operator mixing (the mixing, under renormalization, of quark operators with gluon operators of the same spin, which couples the gluon distribution into the moment analysis) was ignored; the perturbative input was truncated at one loop with approximate coefficient functions; the fitted quark distribution F(ξ)F(\xi) was unconstrained at ξ≳0.8\xi\gtrsim 0.8 where no data existed; the RR data at large xx and Q2Q^{2} sat systematically above the prediction; and the proton matrix elements AnA_{n} came from fits, not from first principles. The open problem is to remove these crutches — to supply, from first principles or from data independent of the resonance region being tested, every input of the DGP local-duality construction, and to test the resulting parameter-free predictions.

Problem Statement

Complete the quantitative test of electroproduction local duality and precocious scaling formulated by De Rújula, Georgi and Politzer (1977), by replacing every phenomenological input of their analysis with a determination from first principles or from independent data, and confronting the resulting parameter-free predictions with precision measurements. The common scientific objective is a falsifiable, assumption-controlled test of the OPE moment basis of local duality (resonance-region averages equal to leading-twist moments plus controlled higher twist) in the moderate-Q2Q^{2} region; the six components below are the jointly necessary elements of that single test:

  1. Existing published determinations of R=σL/σTR=\sigma_{L}/\sigma_{T} at large Nachtmann ξ\xi and Q2Q^{2} to test the twist-two longitudinal prediction (the 1977 SLAC data lay systematically above that prediction at x≥0.67x\geq 0.67, Q2≥8Q^{2}\geq 8 GeV2^{2}).
  2. A quantitative determination of the gluon distribution at small ξ\xi and a controlled inclusion of quark–gluon operator mixing, so the moment analysis extends below ξ≲0.3\xi\lesssim 0.3, where the 1977 analysis neglected gluons by assumption.
  3. Higher-order perturbative input — anomalous dimensions, coefficient functions, and the β\beta function beyond the one-loop terms used in 1977 — so the theory side is controlled as αs(Q2)\alpha_{s}(Q^{2}) grows at lower Q2Q^{2}.
  4. Determination of the twist-four and higher-twist operator matrix elements Bnk(Q2)B_{nk}(Q^{2}), including their magnitudes, their nn dependence, and the higher-twist scale M0M_{0} entering the terms (nM02/Q2)kBnk(Q2)(nM_{0}^{2}/Q^{2})^{k}B_{nk}(Q^{2}), from a nonperturbative, first-principles determination, so the power-suppressed terms are determined rather than parametrized.
  5. Existing published structure-function data at large ξ\xi (ξ≳0.8\xi\gtrsim 0.8) in the scaling region, to test the predicted tail of the scaling function F(ξ)F(\xi) and thereby verify or violate the elastic-form-factor predictions obtained from the threshold moment relations.
  6. A nonperturbative, first-principles determination of the proton matrix elements AnA_{n} and of the quark/gluon operator separation (the decomposition of each twist-two moment into its quark and gluon operator matrix elements, distinct from the gluon distribution and perturbative mixing of component 2), replacing phenomenological fits. The local-duality observables — resonance-region averages of the structure functions, R(ξ,Q2)R(\xi,Q^{2}) at large ξ\xi, and the elastic form factors through the threshold ξ\xi-moment relations — are predicted with quantified uncertainties from inputs obtained without using the data being predicted, and are either verified or quantitatively violated by those data. Each determination in components 4 and 6 carries a quantified uncertainty, and the prediction uncertainties are small enough to resolve the higher-twist correction being tested: they do not exceed the magnitude of the (nM02/Q2)kBnk(nM_{0}^{2}/Q^{2})^{k}B_{nk} term at the quoted Q2Q^{2}.

References

  1. A. De Rújula, H. Georgi, H. D. Politzer, 'Demythification of electroproduction local duality and precocious scaling', Annals of Physics 103 (1977) 315–353. DOI: 10.1016/S0003-4916(97)90003-8. https://doi.org/10.1016/S0003-4916(97)90003-8
  2. E. D. Bloom, F. J. Gilman, 'Scaling and the behavior of nucleon resonances in inelastic electron-nucleon scattering', Physical Review D 4 (1971) 2901–2916. DOI: 10.1103/PhysRevD.4.2901. https://doi.org/10.1103/PhysRevD.4.2901
  3. O. Nachtmann, 'Positivity constraints for anomalous dimensions', Nuclear Physics B 63 (1973) 237–247 (introduces the Nachtmann variable ξ\xi and the finite-Q2Q^{2} projection onto operators of definite spin). DOI: 10.1016/0550-3213(73)90144-2. https://doi.org/10.1016/0550-3213(73)90144-2
  4. W. Melnitchouk, R. Ent, C. E. Keppel, 'Quark-hadron duality in electron scattering', Physics Reports 406 (2005) 127–301. arXiv:hep-ph/0501217. DOI: 10.1016/j.physrep.2004.10.004. https://arxiv.org/abs/hep-ph/0501217
  5. R. Ent, C. E. Keppel, I. Niculescu, 'Nucleon elastic form factors and local duality', Physical Review D 62 (2000) 073008. DOI: 10.1103/PhysRevD.62.073008. https://doi.org/10.1103/PhysRevD.62.073008
  6. S. Simula, 'Comment on "Nucleon elastic form factors and local duality"', Physical Review D 64 (2001) 038301. arXiv:hep-ph/0012193. DOI: 10.1103/PhysRevD.64.038301. https://arxiv.org/abs/hep-ph/0012193
  7. S. Simula, G. Ricco, M. Anghinolfi, M. Ripani, M. Taiuti, 'Bloom-Gilman duality of the nucleon structure function and the elastic peak contribution', in 'Few-Body Problems in Physics 98', Few Body Systems Supplement 10 (1999) 423–426. arXiv:hep-ph/9909547. DOI: 10.1007/978-3-7091-6798-4_79. https://doi.org/10.1007/978-3-7091-6798-4_79
  8. I. Niculescu et al. (Jefferson Lab Collaboration), 'Experimental verification of quark-hadron duality', Physical Review Letters 85 (2000) 1186–1189. DOI: 10.1103/PhysRevLett.85.1186. https://doi.org/10.1103/PhysRevLett.85.1186
  9. S. Liuti, R. Ent, C. E. Keppel, I. Niculescu, 'Perturbative QCD analysis of local duality in a fixed W2W^{2} framework', Physical Review Letters 89 (2002) 162001. arXiv:hep-ph/0111063. DOI: 10.1103/PhysRevLett.89.162001. https://arxiv.org/abs/hep-ph/0111063
  10. Y. Liang et al. (Jefferson Lab Hall C E94-110 Collaboration), 'Measurement of R=σL/σTR=\sigma_{L}/\sigma_{T} and the separated longitudinal and transverse structure functions in the nucleon-resonance region', Physical Review C 105 (2022) 065205. arXiv:nucl-ex/0410027. DOI: 10.1103/PhysRevC.105.065205. https://arxiv.org/abs/nucl-ex/0410027
  11. S. Moch, J. A. M. Vermaseren, A. Vogt, 'The three-loop splitting functions in QCD: the non-singlet case', Nuclear Physics B 688 (2004) 101–134. arXiv:hep-ph/0403192. DOI: 10.1016/j.nuclphysb.2004.03.030. https://arxiv.org/abs/hep-ph/0403192
  12. S. Capitani, M. Göckeler, R. Horsley, D. Petters, D. Pleiter, P. E. L. Rakow, G. Schierholz, 'Higher twist corrections to nucleon structure functions from lattice QCD', Nuclear Physics B Proceedings Supplements 79 (1999) 173–175. arXiv:hep-ph/9906320. DOI: 10.1016/S0920-5632(99)00666-0. https://arxiv.org/abs/hep-ph/9906320
  13. X. Ji, 'Parton physics on a Euclidean lattice', Physical Review Letters 110 (2013) 262002. arXiv:1305.1539. DOI: 10.1103/PhysRevLett.110.262002. https://arxiv.org/abs/1305.1539
  14. Z. Fan, W. Good, H.-W. Lin, 'Gluon parton distribution of the nucleon from (2+1+1)(2+1+1)-flavor lattice QCD in the physical-continuum limit', Physical Review D 108 (2023) 014508. arXiv:2210.09985. DOI: 10.1103/PhysRevD.108.014508. https://arxiv.org/abs/2210.09985
  15. J. F. Owens, A. Accardi, W. Melnitchouk, 'Global parton distributions with nuclear and finite-Q2Q^{2} corrections', Physical Review D 87 (2013) 094012. arXiv:1212.1702. DOI: 10.1103/PhysRevD.87.094012. https://arxiv.org/abs/1212.1702
  16. A. Courtoy, S. Liuti, 'Extraction of αs\alpha_{s} from deep inelastic scattering at large xx', Physics Letters B 726 (2013) 320–325. arXiv:1302.4439. DOI: 10.1016/j.physletb.2013.08.058. https://arxiv.org/abs/1302.4439
  17. M. J. Glatzmaier, S. Mantry, M. J. Ramsey-Musolf, 'Higher twist in electroproduction: flavor nonsinglet QCD evolution', Physical Review C 88 (2013) 025202. arXiv:1208.2998. DOI: 10.1103/PhysRevC.88.025202. https://arxiv.org/abs/1208.2998
  18. E. Moffat, W. Melnitchouk, T. C. Rogers, N. Sato, 'What are the low-QQ and large-xx boundaries of collinear QCD factorization theorems?', Physical Review D 95 (2017) 096008. arXiv:1702.03955. DOI: 10.1103/PhysRevD.95.096008. https://arxiv.org/abs/1702.03955
  19. A. B. Sheckler, G. A. Miller, 'Mystery of Bloom-Gilman duality: a light-front holographic QCD perspective', Physical Review D 103 (2021) 096018. arXiv:2101.00100. DOI: 10.1103/PhysRevD.103.096018. https://arxiv.org/abs/2101.00100
  20. V. Lagerquist, S. E. Kuhn, N. Sato, 'Detailed study of quark-hadron duality in spin structure functions of the proton and neutron', Physical Review C 107 (2023) 045201. arXiv:2205.01218. DOI: 10.1103/PhysRevC.107.045201. https://arxiv.org/abs/2205.01218
  21. M. Cerutti, A. Accardi, I. P. Fernando, S. Li, J. F. Owens, S. Park, 'Systematic uncertainties from higher-twist corrections in DIS at large xx', Physical Review D 111 (2025) 094013. arXiv:2501.06849. DOI: 10.1103/PhysRevD.111.094013. https://arxiv.org/abs/2501.06849
  22. A. Accardi, M. Cerutti, C. E. Keppel, S. Li, J. F. Owens, S. Park, P. Risse, 'CJ26 global QCD analysis with large-xx Jefferson Lab 6 and 12 GeV data', arXiv:2605.31424 (2026). https://arxiv.org/abs/2605.31424
  23. C. Cocuzza, W. Melnitchouk, N. Sato, A. W. Thomas, 'Stability of parton distributions at high xx: impact of nuclear and power corrections', arXiv:2605.00666 (2026). https://arxiv.org/abs/2605.00666
  24. S. I. Alekhin, S. A. Kulagin, R. Petti, 'On the interplay of nuclear and higher-twist corrections in nuclear structure functions', Proceedings of Science DIS2024 (2025) 067. arXiv:2410.21472. DOI: 10.22323/1.469.0067. https://doi.org/10.22323/1.469.0067
  25. Y. Wunderlich, A. Bulgakov, K. Joo, T.-S. H. Lee, V. I. Mokeev, 'Quark-hadron duality in inclusive electron-proton scattering at high Q2Q^{2}: structure functions and truncated moments from CLAS12', arXiv:2606.29690 (2026). https://arxiv.org/abs/2606.29690
  26. A. De Rújula, H. Georgi, H. D. Politzer, 'Trouble with ξ\xi scaling?', Physical Review D 15 (1977) 2495 (companion paper restating that the higher-twist corrections to the ξ\xi-scaling analysis are negligible only for WW well above the proton mass). DOI: 10.1103/PhysRevD.15.2495. https://doi.org/10.1103/PhysRevD.15.2495