REVIEW 3 major objections 4 minor 31 references
Transverse force distributions in the proton from lattice QCD
T0 review · 3 major / 4 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read Lattice QCD maps the colour-Lorentz force on a struck quark inside the proton, finding a central attractive force that plateaus near 3 GeV/fm.
desk verdict A proceedings summary of a real lattice calculation; the ~3 GeV/fm force magnitude is a model-dependent estimate, not a direct lattice result. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the off-forward matrix element of the twist-three light-cone operator $\bar{q}(0)\gamma^+ i g G^{+j}(0) q(0)$, decomposed into three form factors by Eq. (4); the operator's value is that it can be read, through the chromodynamic lensing picture, as the colour-Lorentz force $F^y = -(E_c + v \times B_c)^y$ acting on the quark at the instant the virtual photon is absorbed. The Fourier transform of these form factors with respect to the transverse momentum transfer $\Delta_\perp$, at zero skewness, converts the distribution into transverse impact parameter space $b_\perp$, giving a spatial image of the force. The extraction then relies on two modelling steps: a dipole (or higher-multipole) ansatz for the $t$-dependence of each form factor, and the factorisation assumption of Eq. (12), $F^i_{s',s}(b_\perp) = \rho_{s',s}(b_\perp) F^i_{s',s}(b_\perp)$, used to remove the quark-density weighting and quote unweighted force magnitudes.
What would settle it
A second lattice calculation that obtains the unweighted force without the factorisation of Eq. (12), for instance by computing the quark density and the force from the same off-forward matrix elements rather than from a dipole fit, would either reproduce the 3 GeV/fm plateau or reveal it as an artifact of the density division.
Extended reading notes
Core claim
Working with $N_f = 2 + 1$ dynamical fermions at the SU(3) symmetric point across three lattice spacings, the authors compute off-forward proton matrix elements of the twist-three operator $\bar{q}(0)\gamma^+ i g G^{+j}(0) q(0)$, whose Cartesian components mix chromo-electric and chromo-magnetic fields and therefore carry the QCD analogue of a Lorentz force. They parameterise these matrix elements in terms of three form factors $\Phi_1(t)$, $\Phi_2(t)$ and $\Phi_3(t)$, and take two-dimensional Fourier transforms at zero skewness to obtain the force distribution in transverse impact parameter space. The $\Phi_1$ term yields a universally attractive, spin-independent force, while the $\Phi_2$ and $\Phi_3$ terms produce a spin-dependent field that resembles a magnetic dipole. After dividing out the quark density, the unpolarised force rises from zero at the origin to a plateau of about 3 GeV/fm at $b \sim 0.15$ fm, remaining roughly constant at larger separations; the authors emphasise that this local force is nearly three times the phenomenological QCD string tension. The up quark has strong signals in all form factors, whereas the down quark's spin-dependent contribution is suppressed, consistent with a quark-diquark picture.
Load-bearing premise
The quoted force magnitudes, including the 3 GeV/fm plateau, rest on the assumption that the density-weighted force separates cleanly into a quark density times an unweighted force, and that the density itself is well described by a dipole fit to electromagnetic form factors; if that separation fails, the reported numbers are not the true local force on the quark.
Editorial extensions
If this is right
- The force map provides a real-space account of the Sivers asymmetry: the largest force pulls the struck up quark back toward the proton core in the region of highest quark density, in the direction opposite to the observed asymmetry of final states.
- The spin-independent plateau near 3 GeV/fm establishes a local confining force scale inside the proton that is markedly larger than the average force scale set by the QCD string tension.
- The suppression of the down quark's spin-dependent form factor supports a quark-diquark picture of the proton, where the down quark sits in a scalar diquark and contributes weakly to spin-dependent observables.
- The dipole-like spin-dependent force field predicts a repulsive force at large peripheral distances on the $+\hat{y}$ side, although the low quark density there makes a direct experimental signature unlikely.
- Extending the momentum range of the form factors would tighten the model dependence near the origin, where the spin-dependent force magnitude currently diverges.
Reading between the lines
- If the factorisation assumption holds, the 3 GeV/fm plateau implies that the force on a struck quark is concentrated in a core region rather than being uniformly spread, suggesting the confinement force field inside the proton is stronger than the flux-tube average.
- Applying the same imaging procedure to the neutron or to strange quarks would test whether the colour-Lorentz force distribution is flavour-blind at short distances and would connect the Sivers asymmetry of different hadrons to the same underlying force field.
- The dipole-like spin-dependent force pattern invites an interpretation in terms of an effective spin-orbit interaction; checking whether the integral of the force field reproduces the known second moment $d_2$ would provide a sum-rule-style consistency test of the imaging method.
- A future experiment measuring the Sivers asymmetry with high precision could, through the chromodynamic lensing relation, be confronted with the lattice force map; a mismatch would indicate that the twist-three force interpretation needs additional QCD corrections.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This proceedings contribution reports a lattice QCD calculation of transverse colour-Lorentz force distributions in the proton. Using Nf=2+1 dynamical clover fermions at the SU(3) symmetric point, the authors compute matrix elements of off-forward twist-three operators, extract three form factors Φ1, Φ2, Φ3, and take two-dimensional Fourier transforms to obtain force distributions in impact parameter space. The central results are: a spin-independent attractive central force that rises to about 3 GeV/fm near b≈0.15 fm and stays roughly constant, and a spin-dependent force pattern resembling a magnetic dipole. These distributions are interpreted in the chromodynamic lensing picture of the Sivers asymmetry. Only results at the finest lattice spacing (a=0.052 fm) are shown; the abstract claims three lattice spacings, with full details deferred to a companion paper (Ref. [2]).
Significance. If the unweighted-force interpretation holds, this work provides a new, intuitive, lattice-QCD-based picture of confinement and of the forces underlying the Sivers asymmetry. The technical machinery—RI'-MOM renormalization, operator mixing, two-state fits, and multipole model variation—is standard and appears to be applied carefully. The paper is transparent about its main assumption (Eq. (12)) and about the near-origin model dependence of the spin-dependent force. However, the headline numerical claim (3 GeV/fm, about three times the string tension) is a density-divided estimate that depends on an untested factorization ansatz and on a dipole model for the quark density; the robust lattice output is the density-weighted distribution. As a result, the significance is real but conditional on the validity and further testing of Eq. (12).
major comments (3)
- [Section 4, Eq. (12)] The central numerical claim of the paper—the unweighted colour-Lorentz force reaching about 3 GeV/fm, nearly three times the string tension—is obtained from Eq. (12), which states F^i_{s',s}(b⊥)=ρ_{s',s}(b⊥)F^i_{s',s}(b⊥) and is introduced with the words 'By assuming'. The left-hand side is the density-weighted quantity actually obtained from the lattice form factors, while the right-hand side is the unweighted force. This factorization is not derived from QCD and is not tested by the multipole variation in Eq. (13), which only changes the parametrization of the numerator. The abstract's phrasing 'We determine a central, spin-independent confining force, as well as spin-dependent force distributions with local forces larger than the QCD string tension' goes beyond what the lattice data alone establish. Please qualify the abstract and conclusions so that the density-weighted distributions are presented as the robust result and the unweighted magnitudes as model-dependent estimates under Eq. (12), or provide a dedicated test of the factorization.
- [Abstract and Sections 3–4] The abstract states that the calculation uses 'three lattice spacings', but the results shown in Figures 1–4 are exclusively from the finest lattice spacing, a=0.052 fm. No continuum extrapolation, scale-setting uncertainty, or systematic error budget is presented in this manuscript; the text instead refers to Ref. [2] for the full lattice details. For a stand-alone submission whose abstract makes a determination claim, this is a load-bearing gap. Please either include a summary of the multi-spacing results and systematic uncertainties, or explicitly state in the abstract and introduction that the results are preliminary, single-spacing demonstrations taken from the companion paper.
- [Section 4, Fig. 4b] The spin-dependent force magnitude along the x-axis diverges as b→0 and is strongly model-dependent between the dipole, tripole, and quadrupole parametrizations; the text itself notes that further Φ3 data at larger momentum are needed. The abstract's unqualified statement about 'spin-dependent force distributions with local forces larger than the QCD string tension' is therefore not supported by a stable quantitative result in the region where the force is largest. The claim should be restricted to the model-independent region around 0.25 fm, or the near-origin divergence should be explicitly acknowledged as a limitation in the abstract and conclusions.
minor comments (4)
- [Section 4, Eq. (12)] The notation in Eq. (12) is confusing because the same symbol F is used for the density-weighted lattice result and the unweighted force, making the equation look like a tautology; please use distinct symbols (e.g., F̃ and F) for the two quantities.
- [Figure 2 and surrounding text] The text states that 'the vectors shown in the figure are weighted by the corresponding quark density', while the figure caption says 'Vectors have unit length and their magnitude is indicated by the colour bar'; these statements are inconsistent and should be reconciled.
- [Abstract and Section 1] The abstract contains a typo: 'complimentary picture' should be 'complementary picture'.
- [Section 4, Eq. (10)] In Eq. (10), the exponent e^{-i b·Δ⊥} appears to be missing bold vector notation for b (the text uses b⊥ elsewhere); please make the vector notation uniform.
Circularity Check
No significant circularity: lattice matrix elements plus Fourier transforms drive the force distributions, and the unweighted force is an explicitly flagged factorization estimate rather than a fitted input recycled as a prediction.
full rationale
The central derivation is self-contained. The paper computes off-forward matrix elements of a twist-three operator on dynamical ensembles, extracts form factors Phi1, Phi2, Phi3 after RI'-MOM renormalization, fits them with multipole ansaetze, and obtains the transverse density-weighted colour-Lorentz force by 2D Fourier transform (Eqs. (10)-(11)). This is the lattice output, not an input. The headline unweighted force profiles (Figure 4) additionally require dividing by a quark density from electromagnetic form factors under the factorization assumption stated at Eq. (12) and introduced with 'By assuming...'; the text explicitly calls this an 'estimate', not a direct determination. Thus the 3 GeV/fm value is a model-dependent estimate, and the paper's own caveats (down-quark Phi3 reconstruction, divergence of the polarized force near the origin, need for larger momentum) acknowledge the limitation. The only self-referential items, Ref. [2] for details of the same calculation and Ref. [15] for a previously computed renormalization constant, are standard and do not smuggle in the target conclusion. The factorization assumption is a correctness/model risk, not a circularity: the unweighted force is not an input used to define the lattice observable. No circular step can be exhibited.
Assumptions & free parameters
free parameters (4)
- Dipole mass Λ1 and amplitude Φ1(0) for the Φ1 form factor (up and down)
- Dipole mass Λ3 and amplitude Φ3(0) for the Φ3 form factor (up quark)
- Multipole order n of the form factor ansatz =
n = 2, 3, 4
- Quark density dipole parameters
assumptions (4)
- domain assumption Twist-three matrix element equals average colour-Lorentz force
- domain assumption Zero-skewness limit Δ+ = 0 gives probabilistic density interpretation
- ad hoc to paper Density-weighted force factorises as F^i = ρ F^i
- ad hoc to paper Dipole or multipole ansatz for form factors
Cite this review
Pith. "Pith review of Transverse force distributions in the proton from lattice QCD." pith.science (2026). https://pith.science/paper/TWXLBYXI
@misc{pith2026250200325,
author = {Pith},
title = {Pith review of: Transverse force distributions in the proton from lattice QCD},
year = {2026},
howpublished = {\url{https://pith.science/paper/TWXLBYXI}},
note = {Machine review of arXiv:2502.00325}
}
abstract
Single-spin asymmetries observed in polarised deep-inelastic scattering are important probes of hadron structure. The Sivers asymmetry provides information about the transverse momentum of the struck quark and can be related to final-state interactions. Understanding these asymmetries at the quark level has been the subject of much interest in QCD phenomenology. In this contribution, we present a lattice QCD calculation of the transverse spatial distribution of a colour-Lorentz force acting on the struck quark in a proton. Our lattice calculations employ $N_f = 2 + 1$ flavours of dynamical fermions at the SU(3) symmetric point across three lattice spacings. We determine a central, spin-independent confining force, as well as spin-dependent force distributions with local forces larger than the QCD string tension. These distributions offer a new, complimentary picture that underlies the Sivers asymmetry in transversely polarised deep-inelastic scattering.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[2]
J.A. Crawford, K.U. Can, R. Horsley, P.E.L. Rakow, G. Schierholz, H. Stüben et al., Transverse force distributions in the proton from lattice QCD, 2408.03621
-
[1]
Transverse Force on Quarks in DIS
M. Burkardt,Transverse force on quarks in deep-inelastic scattering, Phys. Rev. D88 (2013) 114502 [1510.03112]
work page Pith review arXiv 2013
-
[3]
HERMES measurements of Collins and Sivers asymmetries from a transversely polarised hydrogen target
HERMES collaboration,HERMES measurements of Collins and Sivers asymmetries from a transversely polarised hydrogen target, in15th International Workshop on Deep-Inelastic Scattering and Related Subjects, pp. 579–582, 6, 2007, DOI [0706.2242]
work page Pith review arXiv 2007
-
[4]
M. Anselmino, M. Boglione, U. D’Alesio, A. Kotzinian, S. Melis, F. Murgia et al.,Sivers Effect for Pion and Kaon Production in Semi-Inclusive Deep Inelastic Scattering, Eur. Phys. J. A39(2009) 89 [0805.2677]
arXiv 2009
-
[5]
COMPASScollaboration,COMPASS Results on Transverse Single-Spin Asymmetries,Czech. J. Phys.56(2006) F33 [hep-ex/0702002]
work page Pith review arXiv 2006
-
[6]
Sivers,Hard-scattering scaling laws for single-spin production asymmetries,Phys
D. Sivers,Hard-scattering scaling laws for single-spin production asymmetries,Phys. Rev. D 43 (1991) 261
work page 1991
-
[7]
J.C.Collins, Fragmentationoftransverselypolarizedquarksprobedintransversemomentum distributions, Nucl. Phys. B396 (1993) 161 [hep-ph/9208213]
arXiv 1993
-
[8]
J.C. Collins,Leading twist single transverse-spin asymmetries: Drell-Yan and deep inelastic scattering,Phys. Lett. B536 (2002) 43 [hep-ph/0204004]
arXiv 2002
Show all 31 references
-
[9]
Burkardt,Chromodynamic lensing and transverse single spin asymmetries, Nucl
M. Burkardt,Chromodynamic lensing and transverse single spin asymmetries, Nucl. Phys. A 735 (2004) 185 [hep-ph/0302144]
2004 arXiv
-
[10]
Jaffe,𝑔2-The nucleon’s other spin-dependent structure function, Comments Nucl
R.L. Jaffe,𝑔2-The nucleon’s other spin-dependent structure function, Comments Nucl. Part. Phys. 19(1990) 239
1990
-
[11]
Jaffe and X.-D
R.L. Jaffe and X.-D. Ji,Studies of the transverse spin dependent structure function𝑔2(𝑥,𝑄 2), Phys. Rev. D43(1991) 724
1991
-
[12]
Blumlein and N
J. Blumlein and N. Kochelev,On the twist -2 and twist - three contributions to the spin dependent electroweak structure functions, Nucl. Phys. B498 (1997) 285 [hep-ph/9612318]
1997 arXiv
-
[13]
Göckeler, R
M. Göckeler, R. Horsley, W. Kurzinger, H. Oelrich, D. Pleiter, P.E.L. Rakow et al.,A lattice calculation of the nucleon’s spin dependent structure function𝑔2 revisited, Phys. Rev. D63 (2001) 074506 [hep-lat/0011091]
2001 arXiv
-
[14]
Göckeler, R
M. Göckeler, R. Horsley, D. Pleiter, P.E.L. Rakow, A. Schafer, G. Schierholz et al., Investigation of the second moment of the nucleon’s𝑔1 and𝑔2 structure functions in two-flavor lattice QCD,Phys. Rev. D72 (2005) 054507 [hep-lat/0506017]. 9 Transverse force distributions in th...
2005 arXiv
-
[15]
Bürger, T
S. Bürger, T. Wurm, M. Löffler, M. Göckeler, G. Bali, S. Collins et al.,Lattice results for the longitudinal spin structure and color forces on quarks in a nucleon, Phys. Rev. D105 (2022) 054504 [2111.08306]
2022 arXiv
-
[16]
Wandzura and F
S. Wandzura and F. Wilczek,Sum rules for spin dependent electroproduction: test of relativistic constituent quarks,Phys. Lett. B72(1977) 195
1977
-
[17]
Cortes, B
J.L. Cortes, B. Pire and J.P. Ralston,Measuring the transverse polarization of quarks in the proton, Z. Phys. C55(1992) 409
1992
-
[18]
Aslan, M
F.P. Aslan, M. Burkardt and M. Schlegel,Transverse force tomography, Phys. Rev. D100 (2019) 096021 [1904.03494]
2019 arXiv
-
[19]
Burkardt,Impact parameter dependent parton distributions and off forward parton distributions for𝜁→0,Phys
M. Burkardt,Impact parameter dependent parton distributions and off forward parton distributions for𝜁→0,Phys. Rev. D62(2000) 071503 [hep-ph/0005108]
2000 arXiv
-
[20]
Bietenholz, V
W. Bietenholz, V. Bornyakov, M. Göckeler, R. Horsley, W.G. Lockhart, Y. Nakamura et al., Flavour blindness and patterns of flavour symmetry breaking in lattice simulations of up, down and strange quarks, Phys. Rev. D84(2011) 054509 [1102.5300]
2011 arXiv
-
[21]
Cundy, M
N. Cundy, M. Göckeler, R. Horsley, T. Kaltenbrunner, A.D. Kennedy, Y. Nakamura et al., Non-perturbative improvement of stout-smeared three flavour clover fermions, Phys. Rev. D 79(2009) 094507 [0901.3302]
2009 arXiv
-
[22]
Baake, B
M. Baake, B. Gemunden and R. Odingen,Structure and representations of the symmetry group of the four-dimensional cube, J. Math. Phys.23(1982) 944
1982
-
[23]
Shuryak and A.I
E.V. Shuryak and A.I. Vainshtein,Theory of power corrections to deep inelastic scattering in quantum chromodynamics: (ii).𝑞−4 effects; polarized target, Nucl. Phys. B201 (1982) 141
1982
-
[24]
Martinelli, C
G. Martinelli, C. Pittori, C.T. Sachrajda, M. Testa and A. Vladikas,A general method for nonperturbative renormalization of lattice operators, Nucl. Phys. B445 (1995) 81 [hep-lat/9411010]
1995 arXiv
-
[25]
Göckeler, R
M. Göckeler, R. Horsley, H. Oelrich, H. Perlt, D. Petters, P.E.L. Rakow et al., Nonperturbative renormalization of composite operators in lattice QCD,Nucl. Phys. B544 (1999) 699 [hep-lat/9807044]
1999 arXiv
-
[26]
Cloët and G.A
I.C. Cloët and G.A. Miller,Nucleon form factors and spin content in a quark-diquark model with a pion cloud,Phys. Rev. C86(2012) 015208
2012
-
[27]
J. Wang, D. Fu and Y. Dong,A systematical study of the nucleon form factors with the pion cloud effect, 2410.14953
-
[28]
Göckeler, P
M. Göckeler, P. Hägler, R. Horsley, Y. Nakamura, D. Pleiter, P.E.L. Rakow et al.,Transverse spin structure of the nucleon from lattice QCD simulations,Phys. Rev. Lett.98 (2007) 222001 [hep-lat/0612032]. 10 Transverse force distributions in the proton from lattice QCD J. A. Crawford
2007 arXiv
-
[29]
Diehl and P
M. Diehl and P. Hägler,Spin densities in the transverse plane and generalized transversity distributions, Eur. Phys. J. C44(2005) 87 [hep-ph/0504175]
2005 arXiv
-
[30]
T.R. Haar, Y. Nakamura and H. Stüben,An update on the BQCD Hybrid Monte Carlo program, EPJ Web Conf.175 (2018) 14011 [1711.03836]
2018 arXiv
-
[31]
Edwards and B
R.G. Edwards and B. Joó,The chroma software system for lattice QCD, Nucl. Phys. B140 (2005) 832 [hep-lat/0409003]. 11
2005 arXiv
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.