收藏切换
Evidence for Quark Confinement in the Proton
收藏切换
PDF
Xiangdong Ji1, *, Gerald A. Miller2, Chen Yang3
Research. Vol 9 Article ID 1369
Less
收藏切换
Research. Vol 9 Article ID 1369
• RESEARCH ARTICLE •
Evidence for Quark Confinement in the Proton
Full
Xiangdong Ji1, *, Gerald A. Miller2, Chen Yang3
Affiliations
  • 1Tsung-Dao Lee Institute and School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 201210, China.
  • 2Department of Physics, University of Washington, Seattle, WA 98195-1560, USA.
  • 3Maryland Center for Fundamental Physics, Department of Physics, University of Maryland, College Park, MD 20742, USA.
About Author:

† These authors contributed equally to this work.

Published: 2026-08-20 doi: 10.34133/research.1369
Outline
收藏切换

The strong interaction is the fundamental force that holds quarks and the gluon force carriers together to form protons and neutrons and also binds the atomic nucleus. The theory governing quark–gluon interactions is quantum chromodynamics (QCD). A wide variety of experimental data teaches us that quarks and gluons cannot be observed in isolation, a phenomenon known as confinement that is unique to QCD. But no one has used QCD to mathematically prove confinement. Here, we show how to define and measure the force on quarks in the proton using available experimental data. Direct evidence for confinement is obtained because the force is found to be attractive and constant for a wide range of quark positions. This work guides future experimental efforts at future electron-ion colliders aimed at obtaining a rigorous quantitative understanding of confinement and the origin of nuclear mass.

Xiangdong Ji, Gerald A. Miller, Chen Yang. Evidence for Quark Confinement in the Proton[J]. Research, 2026 , 9 (8) : 1369 - . DOI: 10.34133/research.1369
Rapid progress in understanding the role of the quark and gluon energy-momentum tensor (EMT) and the associated momentum conservation laws [1–4], along with the newly found ability to measure its properties [5–8], has recently occurred. But the link between these properties and the underlying mystery of quark confinement has not been established.
There are several important reasons for this. The difficulties can be understood by comparing the present situation with the history of the most prominent forces in our daily lives—gravity and electromagnetism. The inverse square laws of gravity and electricity were both discovered in table-top laboratory experiments by Cavendish and Coulomb using a torsion balance [9,10], but it is impossible to place quarks on a torsion balance and directly read the forces between them. This is because of the effects of confinement [11–13] that prevent the observation of a free quark. Furthermore, in the small-sized world of the proton, the uncertainty principle precludes finding quarks at a definite position. To make matters even more difficult, the quarks within the proton move at speeds very close to those of light. Einstein’s theory of relativity must be respected.
Thus, the fundamental strong forces acting between colored quarks have never been accessed by direct experimental measurement, but instead are postulated within a theory known as quantum chromodynamics (QCD) that operates in the sub-atomic quantum domain. The interactions between quarks separated by small distances, namely, asymptotic freedom [14,15], have been studied using high-energy scattering, but at larger distances, our knowledge is largely limited to models of hadron spectroscopy [16–18] and models of static potentials for massive quarks [19]. Finding a way to measure the forces on quarks that is directly connected to experimental measurements is therefore necessary to study inter-quark forces at large separations without using models. If this is accomplished, a test of the fundamental theory and understanding of the QCD confinement mechanism would be obtained.
In the following, we show how to overcome the various roadblocks mentioned above. The first step is to understand the proper definition of forces. This is done using the QCD EMT, the strong-interaction version of the well-known Maxwell stress tensor used to understand electromagnetic forces. The nucleon matrix element of the EMT is defined in terms of kinematic factors and relativistically invariant form factors with the specific forms chosen so that the conservation of energy and momentum is respected [1]. Then, one needs to obtain the form factors from various experimental measurements [20–22] and fundamental numerical calculations known as lattice QCD [11]. This information allows a force density to be obtained. The force itself is obtained by dividing its density by the density of quarks [23], with the latter also obtained from experimental measurements. All this is not yet enough. It is necessary to obtain an unambiguous way to define the position of a quark that respects the theory of relativity and the uncertainty principle. This is done by making evaluations in an infinite momentum frame (IMF) [24–26]. The net result is a well-defined expression for the force on quarks in a proton that can be evaluated using the relevant input information available.
We shall extend a method to extract the force on quarks that has been well-tested in nonrelativistic, quantum mechanical systems [27]. Force densities are obtained from the divergence of the EMT [28], which receives contributions from quarks, gluons, and an anomalous term [2]. For a stable system, the divergence of the total EMT vanishes, but the individual terms have nonvanishing divergences, giving rise to force components. Our interest here is in the forces on the quarks. The quark (q) part of the EMT is traceless and is given by
Tqμν=12ψ¯iD↔ (μγν)ψ,
The divergence of this EMT represents the local momentum change of the quarks. Analogous to Newton’s law, it thus yields [3,4] the density of the forces acting on the quarks:
Fqi≡∂μTqμi=gψ¯γμFμiψ
=gψ¯γ0F0iψ+gψ¯γjFjiψ
=gρcE→c+j→c×B→ci
where c is the color index, i represents a spatial variable, g is the strong coupling constant, ρcj→c are the quark color-charge and color-current densities, and E→cB→c are the chromo-electromagnetic fields. The right-hand side of Eq. (4) is the gluon version of the Lorentz force law, or the color-Lorentz force [3,4], familiar from electromagnetism.
Nucleon momentum densities and flows are measured through high-energy experiments as the matrix elements of the EMT between eigenstates of momentum, parametrized in terms of kinematic terms and EMT form factors, expressed as A,B,C,C¯ [1]. These are experimentally accessible through exclusive processes such as deeply virtual Compton scattering (DVCS) and deep exclusive meson production (DVMP) [1,5–8,20–22,29–33]. The form factors that describe the EMT can also be computed using first-principles lattice-QCD (LQCD) calculations [34].
The force density, Fqj, as shown in Eq. (4), is obtained from the term in the quark part of the EMT with a nonvanishing divergence; these are denoted as nonconserved (nc). The conserved terms have vanishing divergence and thus give no contribution to the force density. The relevant term is the C¯q form factor in the matrix elements [1],
P′TqμνPnc=U¯P′gμνMC¯qq2UP=U¯P′−14gμνMGs,qq2UP,
where U and U¯ are proton momentum-spin wave functions, M is the proton mass, qμ=P′μ−Pμ is the momentum transfer, and C¯qq2=−Gs,qq2/4 is the quark scalar form factor obtained from the constraint that the quark EMT is traceless [35],
Gs,qq2=Aqq2+Bqq2q24M2−Cqq23q2M2,
with q2=P′−P2. The form factor Bqq2 is relatively small and is therefore neglected in the following calculations. See the Supplementary Materials for details on the matrix elements and EMT form factors.
Experimental data and LQCD only measure the matrix elements of the EMT and force density Fqj between momentum states, as shown in Eq. (5), but such states do not have a well-defined position. Furthermore, the different momenta of the initial and final states lead to different internal dynamics and structures due to relativity, while a true density is defined by a matrix element between the same states.
To handle this, the proton state vector, for both the initial and final states, is taken to be an infinite superposition of states with fixed infinite momentum in one direction (z) and any momentum transverse to z: We use an infinite momentum frame [24–26]. In this frame, the transverse momentum of the system factorizes out of the wave function so that states of any momentum transverse to z have the same internal dynamics. This superposition yields a state of well-defined transverse position that can be taken to be the origin so that transverse positions r→⊥ relative to the origin are defined with the associated spatial distributions in this transverse plane [36–38], as shown in Fig. 1. The nonconserved part of the quark EMT in the transverse plane is then
Tqμνncr→⊥=−14gμνMGs,qr→⊥
which is given by the scalar form factor, and its spatial density Gs,qr→⊥ is obtained by Fourier–Bessel transforming Gs,q−q⊥2. Details are provided in the Supplementary Materials.
In the 4-dimensional divergence appearing in Eq. (4), 2 of these derivatives ∂+ and ∂− vanish when matrix elements are taken between the proton state vector. This is because the energy transfer vanishes because P+→∞, and the longitudinal momentum is fixed at an infinite value. On the other hand, the nonconserved part is proportional to gμν as in Eq. (7). Therefore, since g+i=g−i=0, only the transverse derivative survives, as shown in Eq. (S20). Thus, we are left with only the transverse parts of the divergence. This means that the force density is simply
Fqir→⊥=∂μTqμir→⊥=M4∇⊥iGs,qr→⊥
See Eqs. (S23) and (S24) for an explicit mathematical derivation.
The force on the quarks is the force density, Eq. (8), divided by the quark density, ρqr⊥. The function ρqr⊥ can be obtained, using the infinite momentum frame discussed above, by taking the 2-dimensional Fourier transform of the sum of the proton and neutron Dirac form factors 3F1pq2+F1nq2 [37], where the factor 3 reflects the number of valence quarks. Therefore, the color-Lorentz force on quarks is given by
F→qr→⊥≡F→qr→⊥ρqr→⊥=M/4∇→⊥Gs,qr→⊥3F1pr→⊥+F1nr→⊥,
where the spatial densities are obtained by Fourier–Bessel transformations of the corresponding form factors. See the Supplementary Materials for details. Due to rotational symmetry in the transverse plane, the force lies exclusively along the radial direction.
Equation (9) shows how to precisely determine the force once all the form factors are known for a large range of values of the momentum transfer. The accuracy of our force determination depends on the current incomplete nature of this knowledge.
We perform a comprehensive assessment of the state-of-the-art global phenomenological analyses, including both LQCD and experimental inputs to extract the quark scalar form factor [5,6,32–34]; the results are summarized in Fig. 2. The LQCD calculations of the EMT form factors (using pion mass mπ=169 MeV, and lattice spacing a = 0.091 fm) are shown as data points with purple error bars [34]. The Aqq2-form factor, which appears as part of Gs,qq2, is rescaled so that Aq0 of LQCD is consistent with the Coordinated Theoretical-Experimental Project on QCD (CTEQ) extraction [39]. The GYZ fit (orange curve with its uncertainty band) is a Bayesian analysis driven primarily by lattice constraints, with only weak sensitivity to near-threshold J/ψ photo-production measurements [32], as shown in figure 3 of Ref. [32]. The experimental information is mainly reflected by the black curve and band, obtained by combining the global extraction of generalized parton distributions (GUMP fit) [33] and the dispersive analysis of DVCS data (BEG fit) [5]. Due to the limited kinematic range available for current experiments, the inferred scalar form factor is poorly constrained at large values of momentum transfer, ∣q2∣. See the Supplementary Materials for numerical details and error analyses. Note that the current GUMP fit does not incorporate an error estimate for Aqq2, and the uncertainties of Gs,qq2 and the force below solely come from the BEG fit of Cqq2.
The force acting on the proton’s quarks in the IMF, using Eq. (9), is displayed in Fig. 3. The denominator of Eq. (9) is obtained from the Dirac proton and neutron form factors of Ref. [40]. The forces extracted from the GYZ fit and GUMP and BEG fits, together with the 1σ deviations, are shown as the orange and black curves and associated shaded areas. For comparison, we also show the constant confining force inferred from several relativistic quark models as the blue hatched band [41–44].
The experimental results from GUMP and BEG fits provide strong evidence for a net confining force acting on quarks. At intermediate distances, specifically r⊥=1.0∼1.4 fm, the total force acting on the quarks extracted from the experimental data is consistent with the QCD linear potential within uncertainties. Averaging over this interval yields F¯q=−0.38292 GeV/fm, where the uncertainty corresponds to a 1σ deviation. With more data availability in the future, a rigorous statistical analysis of this constant force behavior can be studied. On the other hand, the GYZ fit, dominated by the LQCD data, produces a mean value F¯q=−0.217−0.015+0.016 GeV/fm with uncertainties 1σ, for the range r⊥=0.7∼1.2 fm. The GYZ fit exhibits relatively small uncertainties, obtained from a Monte Carlo replica method. Note that these uncertainties may be underestimated as the systematic errors in the lattice calculations are not thoroughly considered. These values are also compatible with relativistic quark models [41–44]. The negative value of the force denotes an attractive potential on the quark because the force is the negative gradient of the potential.
Extending the curves to larger values of r⊥ than shown in Fig. 3 would involve a substantial increase in the uncertainties, particularly for experimental extraction. This illustrates the challenges in measuring the EMT form factors. These can only be accessed through processes sensitive to generalized parton distributions [1,21]. Measuring these processes generally requires data from high-energy and high-luminosity facilities with detectors capable of precise and complete mapping of scattering events. Therefore, data from the future electron-ion colliders [45,46] are crucially needed to improve the precision of EMT form factors and thereby achieve a rigorous quantitative understanding of color confinement.
The approximately constant strong force on the component quarks observed far away from the proton center that we find is consistent with the known QCD linear potential acting between heavy quarks. This confinement force is about half the QCD string tension σ≈0.9GeV/fm modeled from heavy quarkonium. This reduction is expected because the motion of the light quarks that make up the proton spreads out the postulated flux tube between the quarks [19] as is also indicated by the hopping expansion [47]. A lower string tension for light quarks is also consistent with the findings in light mesons [48].
We note that the GYZ fit, based mainly on LQCD results, yields a weaker confining force than the one inferred from the experimental data. This difference is plausibly attributable to various lattice artifacts such as excited state contamination, the large pion mass, and the finite lattice spacings used, as Aq0 shows a significant deviation from the data.
In summary, we show how to directly determine the strong force on quarks within hadrons using momentum current distributions and associated form factors in Eq. (9) accessible through high-energy experiments and lattice QCD calculations. By applying our approach to the proton in the infinite-momentum frame, we obtain the strong-interaction force acting on quarks and strong evidence for quark confinement. The total confinement force is consistent with the linear potential nature of QCD, although with significant experimental uncertainties. The future electron-ion colliders and other experiments will significantly improve this extraction by extending the kinematic coverage of the DVCS process up to Q2 ∼ 103 GeV2 and expanding the momentum transfer ∣q2∣ down to 0.02 GeV2 and Bjorken-x range down to 10−4. This will significantly reduce the uncertainties in the extracted quark scalar form factor Gs,q in the large-r⊥ region, enabling a high-precision tomography of the quark confinement force.
Details of the materials and methods used in this article are available in the Supplementary Materials.
1.
Ji XD.Gauge-invariant decomposition of nucleon spin. Phys Rev Lett. 1997;78(4):610–613.
2.
Ji X, Yang C.Momentum flow mechanisms and color-Lorentz forces on quarks in the nucleon. Research. 2026;9:1155.
3.
Tanaka K.Operator relations for gravitational form factors of a spin-0 hadron. Phys Rev D. 2018;98(3):Article 034009.
4.
Polyakov MV, Son HD.Nucleon gravitational form factors from instantons: Forces between quark and gluon subsystems. J High Energy Phys. 2018;2018(9):156.
5.
Burkert VD, Elouadrhiri L, Girod FX.The pressure distribution inside the proton. Nature. 2018;557(7705):396–399.
6.
Kumerički K.Measurability of pressure inside the proton. Nature. 2019;570(7759):E1–E2.
7.
Duran B, Meziani ZE, Joosten S, Jones MK, Prasad S, Peng C, Armstrong W, Atac H, Chudakov E, Bhatt H, et al.Determining the gluonic gravitational form factors of the proton. Nature. 2023;615(7954):813–816.
8.
Adhikari S, Afzal F, Akondi CS, Albrecht M, Amaryan M, Arroyave V, Asaturyan A, Austregesilo A, Baldwin Z, Barbosa F, et al.Measurement of the J/ψ photoproduction cross section over the full near-threshold kinematic region. Phys Rev C. 2023;108(2):Article 025201.
9.
Coulomb C. Premier-[troisième] mémoire sur l’electricité et le magnétisme. Nineteenth Century Collections Online (NCCO): Science, Technology, and Medicine: 1780-1925. Paris: Académie Royale des sciences; 1785.
10.
Cavendish H. The electrical researches of the honourable Henry Cavendish, FRS. Cambridge: University Press; 1879.
11.
Wilson KG.Confinement of quarks. Phys Rev D. 1974;10(8):2445–2459.
12.
Polyakov AM.Quark confinement and topology of gauge groups. Nucl Phys B. 1977;120(3):429–458.
13.
Greensite J.The confinement problem in lattice gauge theory. Prog Part Nucl Phys. 2003;51(1):1–83.
14.
Gross DJ, Wilczek F.Ultraviolet behavior of nonabelian gauge theories. Phys Rev Lett. 1973;30(26):1343–1346.
15.
Politzer HD.Reliable perturbative results for strong interactions? Phys Rev Lett. 1973;30(26):1346–1349.
16.
Eichten E, Gottfried K, Kinoshita T, Lane KD, Yan TM.Charmonium: The model. Phys Rev D. 1978;17:3090.
17.
Godfrey S, Isgur N.Mesons in a relativized quark model with chromodynamics. Phys Rev D. 1985;32(1):189–231.
18.
Brambilla N, Krämer M, Mussa R, Vairo A, Bali G, Bodwin GT, Braaten E, Eichten E, Eidelman S, Godfrey S, et al. Heavy quarkonium physics. Geneva: CERN; 2005.
19.
Bali GS.QCD forces and heavy quark bound states. Phys Rep. 2001;343(1–2):1–136.
20.
Ji XD.Deeply virtual Compton scattering. Phys Rev D. 1997;55(11):7114–7125.
21.
Müller D, Robaschik D, Geyer B, Dittes F, Horejsi J.Wave functions, evolution equations and evolution kernels from light ray operators of QCD. Fortsch Phys. 1994;42(2):101–141.
22.
Radyushkin AV.Asymmetric gluon distributions and hard diffractive electroproduction. Phys Lett B. 1996;385:333–342.
23.
Hofstadter R, McAllister RW.Electron scattering from the proton. Phys Rev. 1955;98(1):217–218.
24.
Dirac PAM.Forms of relativistic dynamics. Rev Mod Phys. 1949;21(3):392–399.
25.
Weinberg S.Dynamics at infinite momentum. Phys Rev. 1966;150(4):1313–1318.
26.
Kogut JB, Soper DE.Quantum electrodynamics in the infinite momentum frame. Phys Rev D. 1970;1(10):2901–2913.
27.
Ji X, Yang C.A journey of seeking pressure and forces in the nucleon. Nucl Phys B. 2026;1024:Article 117342.
28.
Landau LD, Lifshitz EM, Kosevich AM, Pitaevskii LP. Theory of elasticity. 3rd ed. Vol 7. Oxford: Butterworth-Heinemann; 1986.
29.
Collins JC, Frankfurt L, Strikman M.Factorization for hard exclusive electroproduction of mesons in QCD. Phys Rev D. 1997;56(5):2982–3006.
30.
Diehl M.Generalized parton distributions. Phys Rep. 2003;388(2–4):41–277.
31.
Belitsky AV, Radyushkin AV.Unraveling hadron structure with generalized parton distributions. Phys Rep. 2005;418(1–6):1–387.
32.
Guo Y, Yuan F, Zhao W.Bayesian inferring nucleon’s gravitation form factors via near-threshold J/ψ photoproduction. Phys Rev Lett. 2025;135(11):Article 111902.
33.
Guo Y, Aslan FP, Ji X, Santiago MG.First global extraction of generalized Parton distributions from experiment and lattice data with next-to-leading-order accuracy. Phys Rev Lett. 2025;135(26):Article 261903.
34.
Hackett DC, Pefkou DA, Shanahan PE.Gravitational form factors of the proton from lattice QCD. Phys Rev Lett. 2024;132(25):Article 251904.
35.
Ji X.Proton mass decomposition: Naturalness and interpretations. Front Phys. 2021;16:64601.
36.
Burkardt M.Impact parameter space interpretation for generalized parton distributions. Int J Mod Phys A. 2003;18(02):173–207.
37.
Miller GA.Charge density of the neutron. Phys Rev Lett. 2007;99:Article 112001.
38.
Miller GA.Impossibility of obtaining time-independent, three-dimensional, spherically symmetric densities of confined systems of relativistically moving constituents. Phys Rev C. 2025;112(4):Article 045204.
39.
Hou TJ, Gao J, Hobbs TJ, Xie K, Dulat S, Guzzi M, Huston J, Nadolsky P, Pumplin J, Schmidt C, et al.New CTEQ global analysis of quantum chromodynamics with high-precision data from the LHC. Phys Rev D. 2021;103(1):Article 014013.
40.
Kelly JJ.Simple parametrization of nucleon form factors. Phys Rev C. 2004;70(6):Article 068202.
41.
Leal Ferreira P, Helayel JA, Zagury N.A linear potential model for quark confinement. Nuovo Cim A. 1980;55(2):215.
42.
Eich E, Rein D, Rodenberg R.Static baryon properties in a relativistic quark model with linearly confining potential. Z Phys. 1985;28:225.
43.
Dziembowski Z, Fabre de la Ripelle M, Miller GA.Nonperturbative gluons and pseudoscalar mesons in baryon spectroscopy. Phys Rev C. 1996;53(5):R2038–R2042.
44.
Glozman LY, Plessas W, Varga K, Wagenbrunn RF.Unified description of light and strange baryon spectra. Phys Rev D. 1998;58(9):Article 094030.
45.
Khalek RA, Accardi A, Adam J, Adamiak D, Akers W, Albaladejo M, Al-Bataineh A, Alexeev MG, Ameli F, Antonioli P, et al.Science requirements and detector concepts for the electron-ion collider: EIC yellow report. Nucl Phys A. 2022;1026:Article 122447.
46.
Anderle DP, Bertone V, Cao X, Chang L, Chang N, Chen G, Chen X, Chen Z, Cui Z, Dai L, et al.Electron-ion collider in China. Front Phys. 2021;16(6):64701.
47.
Rothe HJ. Lattice gauge theories: An introduction. 4th ed. Vol. 82. Singapore: World Scientific Publishing Co. Pte. Ltd.; 2012.
48.
Kawanai T, Sasaki S.Interquark potential with finite quark mass from lattice QCD. Phys Rev Lett. 2011;107(9):Article 091601.
Year 2026 volume 9 Issue 8
PDF
102
57
Cite this Article
BibTeX
Article Info
doi: 10.34133/research.1369
  • Receive Date:2026-05-14
  • Online Date:2026-10-16
  • Published:2026-08-20
Article Data
Affiliations
History
  • Received:2026-05-14
  • Revised:2026-07-02
  • Accepted:2026-07-02
Affiliations
    1Tsung-Dao Lee Institute and School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 201210, China.
    2Department of Physics, University of Washington, Seattle, WA 98195-1560, USA.
    3Maryland Center for Fundamental Physics, Department of Physics, University of Maryland, College Park, MD 20742, USA.

Corresponding:

* Address correspondence to:
References
Share
https://castjournals.cast.org.cn/joweb/research/EN/10.34133/research.1369
Share to
QR

Scan QR to access full text

Cite this article
BibTeX
Citations
表12种不同金属材料的力学参数
科
Family
属数
Number of
genus
种数
Number of
species
占总种数比例
Percentage of
total species (%)
属
Genus
种数
Number of
species
占总种数比例
Percentage of total
species (%)
鹅膏菌科Amanitaceae 2 11 5.26 鹅膏菌属 Amanita 10 4.78
小菇科 Mycenaceae 2 12 5.74 丝盖伞属 Inocybe 5 2.39
多孔菌科 Polyporaceae 8 14 6.70 蜡蘑属 Laccaria 5 2.39
红菇科 Russulaceae 3 23 11.00 小皮伞属 Marasmius 6 2.87
小菇属 Mycena 11 5.26
光柄菇属 Pluteus 5 2.39
红菇属 Russula 17 8.13
栓菌属 Trametes 5 2.39
关闭全屏
  • BibTeX
  • EndNote
  • RefWorks
  • TxT