The quark and gluon gravitational form factors (GFFs) of the nucleon encode how mass, momentum, and mechanical pressure are spatially distributed among quarks and gluons. A PDF-grounded extraction (Selyugin and Teryaev, 2024) found a dipole-like quark GFF but a tripole-like (n≈3n\approx3) gluon GFF with a gluon gravitational radius close to the proton electromagnetic radius, conjecturally caused by an extra gluonic propagator in quark–gluon Fock components of the nucleon wave function, and reliable only for momentum transfers ∣t∣≲2 GeV2|t|\lesssim2~\mathrm{GeV}^2. The open problem is to determine the tt-dependence of Aq(t)A_q(t) and Ag(t)A_g(t) — the precise multipole power and scale of the gluon GFF, its behavior at larger ∣t∣|t|, and the QCD mechanism behind the differing quark and gluon shapes — with controlled systematics from lattice QCD, near-threshold quarkonium photoproduction, and other reactions. Later work (2023–2025) leaves all three elements open: analyses still fix the multipole power by assumption, lattice results stop at ∣t∣=2 GeV2|t|=2~\mathrm{GeV}^2, radius determinations conflict, and rival frameworks predict either stronger or identical gluon versus quark tt-dependence.

Background

Generalized parton distributions (GPDs) are functions F(x,ξ,t)F(x,\xi,t) that describe how partons (quarks and gluons) with longitudinal momentum fraction xx are distributed in transverse space in a hadron, with tt the square of the momentum transfer and ξ\xi the skewness (longitudinal momentum asymmetry between initial and final hadron). Their second Mellin moments in xx at ξ=0\xi=0 give the matrix elements of the QCD energy–momentum tensor (the conserved current coupling to gravity), parameterized by the gravitational form factors (GFFs) A(t)A(t), B(t)B(t), and D(t)D(t) (X. Ji, Phys. Rev. D 55, 7114 (1997)). A(0)A(0) obeys the momentum sum rule — quark plus gluon contributions Aq(0)+Ag(0)=1A_q(0)+A_g(0)=1 — so the slope of Aq,g(t)A_{q,g}(t) at t=0t=0 defines quark and gluon 'gravitational radii' of the mass–energy distribution, while D(t)D(t) encodes internal pressure and shear forces (the 'mechanical' structure). Because gravity is far too weak to probe hadrons directly, GFFs must be accessed indirectly. The prime experimental route is near-threshold heavy-quarkonium photoproduction, e.g. γp→J/ψ p\gamma p\to J/\psi\, p: at threshold the process is gluon-dominated, so the cross section constrains the gluon GFFs; this was first exploited in the Jefferson Lab J/psi-007 extraction (Duran et al., Nature 615, 813 (2023)), which parameterized the gluon GFFs Ag(t)A_g(t) and C(t)=−D(t)/(2A(t))C(t)=-D(t)/(2A(t)) in a tripole (multipole power n=3n=3) form, and in the updated analysis of Guo, Ji, Liu, and Yang (Phys. Rev. D 108, 034003 (2023)). On the theory side, lattice QCD (discrete-space-time numerical simulation of the QCD path integral) can in principle determine the flavor-separated GFFs from first principles. A convenient phenomenological shape is the multipole ansatz A(t)=[1−t/L2]−nA(t)=[1-t/L^2]^{-n}: n=1n=1 monopole, n=2n=2 dipole, n=3n=3 tripole; larger nn means a harder (more compact) large-∣t∣|t| falloff. Perturbative quark-counting rules (Brodsky–Farrar/Matveev–Muradyan–Tavkhelidze scaling) predict power-law falloffs of form factors at very large momentum transfer and motivate connecting observed powers to the number of propagators in the underlying partonic process. Selyugin and Teryaev (arXiv:2311.12988, Phys. At. Nucl. 87, 518 (2024)) built GPD tt-profiles anchored to parton distribution functions (PDFs, the forward limit xx-distributions measured in deep-inelastic scattering) and validated against the 912-point electromagnetic form factor database; taking the limit of zero skewness they obtained Aq(t)A_q(t) in dipole form with L2=1.6 GeV2L^2=1.6~\mathrm{GeV}^2, while the gluon GFF Ag(t)A_g(t), propagated from several modern gluon PDFs, was best described by a multipole with power nn closest to 3 and L2≈0.9 GeV2L^2\approx0.9~\mathrm{GeV}^2, with a gluon gravitational radius of 0.850.85–0.87 fm0.87~\mathrm{fm}, approximately equal to the proton electromagnetic charge radius. They conjectured that the stronger gluonic tt-dependence arises from an additional gluonic propagator associated with quark–gluon Fock components (higher particle-number components) of the nucleon wave function, in line with counting rules, and stated that their profile is validated only for ∣t∣≲2 GeV2|t|\lesssim2~\mathrm{GeV}^2 and that 'various experimental, theoretical and model tests of this picture are highly desirable'.

Problem Statement

Determine the momentum-transfer dependence of the quark and gluon gravitational form factors Aq(t)A_q(t) and Ag(t)A_g(t) of the nucleon, so that the following linked elements are settled with controlled systematics: (i) the shape of the gluon GFF — whether it indeed follows a multipole form with power n≈3n\approx3 (tripole) and scale L2≈0.9 GeV2L^2\approx0.9~\mathrm{GeV}^2 while the quark GFF is dipole-like (n=2n=2, L2≈1.6 GeV2L^2\approx1.6~\mathrm{GeV}^2), i.e., a determination of the exponent nn and scale L2L^2 from data and/or lattice QCD with quantified rather than assumed values; (ii) the behavior of both form factors at momentum transfers beyond the region ∣t∣≲2 GeV2|t|\lesssim2~\mathrm{GeV}^2 where PDF-grounded GPD tt-profiles are validated and lattice data exist, including whether the dipole-versus-tripole distinction persists or is an artifact of the fitted window; and (iii) the microscopic QCD origin of the relative quark/gluon tt-dependence — in particular whether the stronger gluonic falloff is caused by the extra gluonic propagator associated with quark-gluon Fock components of the nucleon wave function, as suggested by quark-counting arguments — together with the resulting quark and gluon gravitational radii. Answers may come from any sound route — combined analyses of near-threshold quarkonium photoproduction and other exclusive reactions, lattice QCD with controlled continuum/chiral/renormalization systematics, holographic or Dyson-Schwinger frameworks, or direct QCD derivations — but must confront, not assume, the multipole shape and must state the ∣t∣|t| domain over which the claimed behavior is established.

References

  1. O. V. Selyugin and O. V. Teryaev, 'GPDs and gravitational form factors of nucleons (quark and gluon contributions)', Physics of Atomic Nuclei 87, 518-522 (2024), arXiv:2311.12988, doi:10.1134/S1063778824700492, https://doi.org/10.1134/S1063778824700492.
  2. X. Ji, 'Deeply virtual Compton scattering', Phys. Rev. D 55, 7114 (1997), doi:10.1103/PhysRevD.55.7114, https://doi.org/10.1103/PhysRevD.55.7114.
  3. B. Duran et al., 'Determining the gluonic gravitational form factors of the proton', Nature 615, 813-816 (2023), doi:10.1038/s41586-023-05730-4, https://doi.org/10.1038/s41586-023-05730-4.
  4. K. A. Mamo and I. Zahed, 'Nucleon mass radii and distribution: Holographic QCD, lattice QCD, and GlueX data', Phys. Rev. D 103, 094010 (2021), arXiv:2103.03186, doi:10.1103/PhysRevD.103.094010, https://arxiv.org/abs/2103.03186.
  5. Y. Guo, X. Ji, Y. Liu, and J. Yang, 'Updated analysis of near-threshold heavy quarkonium production for probe of proton's gluonic gravitational form factors', Phys. Rev. D 108, 034003 (2023), arXiv:2305.06992, doi:10.1103/PhysRevD.108.034003, https://arxiv.org/abs/2305.06992.
  6. D. C. Hackett, D. A. Pefkou, and P. E. Shanahan, 'Gravitational form factors of the proton from lattice QCD', Phys. Rev. Lett. 132, 251904 (2024), arXiv:2310.08484, doi:10.1103/PhysRevLett.132.251904, https://arxiv.org/abs/2310.08484.
  7. Y. Guo, F. Yuan, and W. Zhao, 'Bayesian inferring nucleon gravitational form factors via near-threshold J/psi photoproduction', Phys. Rev. Lett. 135, 111902 (2025), arXiv:2501.10532, doi:10.1103/3x7r-ythq, https://doi.org/10.1103/3x7r-ythq.
  8. X.-H. Cao, F.-K. Guo, Q.-Z. Li, and D.-L. Yao, 'Dispersive determination of nucleon gravitational form factors', Nature Communications 16, 6979 (2025), arXiv:2411.13398, doi:10.1038/s41467-025-62278-9, https://arxiv.org/abs/2411.13398.
  9. Z.-Q. Yao, Y.-Z. Xu, D. Binosi, Z.-F. Cui, M. Ding, K. Raya, C. D. Roberts, J. Rodriguez-Quintero, and S. M. Schmidt, 'Nucleon gravitational form factors', Eur. Phys. J. A 61, 92 (2025), arXiv:2409.15547, doi:10.1140/epja/s10050-025-01557-x, https://arxiv.org/abs/2409.15547.
  10. C. Lorcé, 'Electromagnetic and gravitational form factors of the nucleon', PoS(SPIN2023)010 (2024), arXiv:2402.00429, doi:10.22323/1.456.0010, https://arxiv.org/abs/2402.00429.
  11. W. Broniowski and E. Ruiz Arriola, 'Gravitational form factors and mechanical properties of the nucleon in a meson dominance approach', Phys. Rev. D 112, 054028 (2025), arXiv:2503.09297, doi:10.1103/ylml-mrlh, https://doi.org/10.1103/ylml-mrlh.