# Proton isovector helicity PDF at NNLO and the twist\-3 moment $\\tilde\{d\}\_2$ from lattice QCD at physical quark masses

**INSPIRE:** [3138451](https://inspirehep.net/literature/3138451)
**arXiv:** [2604\.00143](https://arxiv.org/abs/2604.00143)

**Authors:** Gao, Xiang, Hanlon, Andrew D\., Mukherjee, Swagato, Petreczky, Peter, Shu, Hai\-Tao, Yao, Fei, Zhang, Rui, Zhao, Yong

**Submitted:** 31 March 2026

**Subjects:**
- hep\-lat
- hep\-ph
- nucl\-ex
- nucl\-th
- Lattice

## Abstract

We present a lattice quantum chromodynamics calculation of the $x$\-dependent isovector quark helicity parton distribution function \(PDF\) of the proton in the large momentum effective theory \(LaMET\) framework\. Through operator product expansion \(OPE\) we also extract the $\\tilde\{d\}\_2$ moment of the twist\-3 PDF $g\_T\(x\)$ for the first time in the $\\overline\{\\rm MS\}$ scheme, which is proportional to the average color Lorentz force experienced by the quark in the proton\. This calculation is performed on a lattice of spacing $a$ \= 0\.076 fm at physical quark masses\. The quasi\-PDF matrix elements are measured in proton states boosted to momenta $P\_z\=\\\{0, 0\.25, 1\.02, 1\.53\\\}$ GeV\. We first extract the lowest few helicity PDF moments from the renormalization\-group \(RG\) invariant ratios of the matrix elements with OPE\. Combined with the matrix elements relevant for $g\_T\(x\)$, we obtain $\\tilde\{d\}\_2\^\{u\-d\}\(2\\ \{\\rm GeV\}\)\=0\.0024\(46\)$ at next\-to\-leading order in $\\overline\{\\rm MS\}$\. Then, the helicity quasi\-PDF matrix elements are renormalized in the hybrid scheme with linear renormalon resummation and Fourier transformed to the $x$\-space after an asymptotic extrapolation\. The quasi\-PDF is perturbatively matched to the $\\overline\{\\rm MS\}$ PDF with RG and threshold resummations at next\-to\-leading power and next\-to\-next\-to\-leading logarithmic accuracies\. After resummations, we determine the PDF in the region $x\\in\[0\.25,0\.75\]$ with controlled systematic uncertainties\. The end\-point regions are then parameterized, combined with the LaMET prediction at moderate $x$, and fitted to the short\-distance matrix elements in coordinate space\.