Pith. sign in

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 →

arxiv 2502.00325 v1 pith:TWXLBYXI submitted 2025-02-01 hep-lat hep-phnucl-th

classification hep-lathep-phnucl-th PACS 12.38.Gc14.20.Dh
keywords latticeQCDcolour-LorentzforceSiversasymmetrytransverseimpactparameterprotonstructuretwist-threeoperatorsformfactors
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper uses lattice QCD to compute the transverse spatial distribution of the colour-Lorentz force that a struck quark feels inside a proton, and argues that this force map underlies the Sivers single-spin asymmetry observed in polarised deep-inelastic scattering. The authors report a central, spin-independent attractive force that saturates at roughly 3 GeV/fm near 0.15 fm, nearly three times the phenomenological QCD string tension, together with a dipole-like spin-dependent distribution. If established, the result turns a twist-three matrix element of QCD into a concrete, visual account of confinement and of the final-state interactions responsible for the Sivers asymmetry.

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.

Watch

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

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

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)
  1. [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.
  2. [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.
  3. [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)
  1. [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.
  2. [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.
  3. [Abstract and Section 1] The abstract contains a typo: 'complimentary picture' should be 'complementary picture'.
  4. [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

0 steps flagged · score 0.0 of 10

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 4 free parameters · 4 assumptions · 0 invented entities

The central claim rests on standard QCD operator analysis and lattice methodology, plus one ad hoc factorization assumption (Eq. 12) and a multipole ansatz for the t-dependence. No new physical entities are introduced.

free parameters (4)
  • Dipole mass Λ1 and amplitude Φ1(0) for the Φ1 form factor (up and down)
    Fitted to lattice data in Eq. (9); controls the spatial extent and normalization of the unpolarized force distribution.
  • Dipole mass Λ3 and amplitude Φ3(0) for the Φ3 form factor (up quark)
    Fitted to lattice data; near-origin behavior of the spin-dependent force is sensitive to this model (Fig. 4b).
  • Multipole order n of the form factor ansatz = n = 2, 3, 4
    Varied in Eq. (13) to assess model dependence; results for F1 are stable while F3 diverges near the origin.
  • Quark density dipole parameters
    The density ρ(b⊥) is obtained from dipole fits to electromagnetic form factors following Refs [28,29]; used in Eq. (12) to divide out density weighting.
assumptions (4)
  • domain assumption Twist-three matrix element equals average colour-Lorentz force
    Identifies d2 and its off-forward generalization with a colour-Lorentz force via the chromodynamic lensing picture [1,18].
  • domain assumption Zero-skewness limit Δ+ = 0 gives probabilistic density interpretation
    Invoked in Section 2 after Eq. (4) following Burkardt [19].
  • ad hoc to paper Density-weighted force factorises as F^i = ρ F^i
    Eq. (12) in Section 4; necessary to extract unweighted force magnitudes, not derived from QCD.
  • ad hoc to paper Dipole or multipole ansatz for form factors
    Eqs. (9) and (13) used to model t-dependence; no theoretical derivation is given.

how reviews work

0 comments
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 reproduced from arXiv: 2502.00325 by the authors.

Figure 1
Figure 1. Form factors computed at our finest lattice spacing with corresponding dipole fits. the dipole model from Equation (9), however the down quark fit needs to be reconstructed from fits to the isovector (𝑢 − 𝑑) and isoscalar (𝑢 + 𝑑) combinations. By taking a 2D Fourier transform of the form factors, we obtain a visualisation of the colour￾Lorentz force distributions in the transverse plane, F 𝑗 𝑠 ′𝑠 = ∫ 𝑑 2Δ⊥ (2𝜋) 2 𝑒 … view at source ↗
Figure 2
Figure 2. Distribution of the colour-Lorentz force acting on an unpolarised up quark in an unpolarised proton, indicated by the vector field. Vectors have unit length and their magnitude is indicated by the colour bar. The vector field is superimposed on the density distribution for the up quark in an unpolarised proton. The second two terms in Equation (11) are proportional to the Dirac bilinear 𝑢(𝑝 ′ , 𝑠′ )𝑖𝜎+𝜇𝑢(𝑝, 𝑠), and … view at source ↗
Figure 3
Figure 3. Distribution of the colour-Lorentz force acting on an unpolarised up quark in a proton polarised in the 𝑥ˆ-direction, superimposed on the density distribution for the up quark in the polarised proton. Image conventions are identical to [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Density-removed force profiles. 2D Fourier transform of these form factors revealed distributions of large, local forces acting on the struck quark. To connect with phenomenology, these forces were interpreted in the context of the Sivers asymmetry, which provided a co…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

31 extracted references · 15 canonical work pages

  1. [2]

    Crawford, K.U

    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

  2. [1]

    Transverse Force on Quarks in DIS

    M. Burkardt,Transverse force on quarks in deep-inelastic scattering, Phys. Rev. D88 (2013) 114502 [1510.03112]

  3. [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]

  4. [4]

    Anselmino, M

    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]

  5. [5]

    COMPASScollaboration,COMPASS Results on Transverse Single-Spin Asymmetries,Czech. J. Phys.56(2006) F33 [hep-ex/0702002]

  6. [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

  7. [7]

    J.C.Collins, Fragmentationoftransverselypolarizedquarksprobedintransversemomentum distributions, Nucl. Phys. B396 (1993) 161 [hep-ph/9208213]

  8. [8]

    Collins,Leading twist single transverse-spin asymmetries: Drell-Yan and deep inelastic scattering,Phys

    J.C. Collins,Leading twist single transverse-spin asymmetries: Drell-Yan and deep inelastic scattering,Phys. Lett. B536 (2002) 43 [hep-ph/0204004]

Show all 31 references
  1. [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]

  2. [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

  3. [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

  4. [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]

  5. [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]

  6. [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...

  7. [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]

  8. [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

  9. [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

  10. [18]

    Aslan, M

    F.P. Aslan, M. Burkardt and M. Schlegel,Transverse force tomography, Phys. Rev. D100 (2019) 096021 [1904.03494]

  11. [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]

  12. [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]

  13. [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]

  14. [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

  15. [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

  16. [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]

  17. [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]

  18. [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

  19. [27]

    J. Wang, D. Fu and Y. Dong,A systematical study of the nucleon form factors with the pion cloud effect, 2410.14953

  20. [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

  21. [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]

  22. [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]

  23. [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

Pith tools

Reviewed August 9, 2026 · model on record in the stance chip above.