Nucleon Helicity Parton Distribution Function in the Continuum Limit with Self-Renormalization

Jack Holligan,Huey-Wen Lin

DOI: https://doi.org/10.1016/j.physletb.2024.138731

2024-05-28

Abstract:We present the first lattice calculation of the nucleon isovector helicity parton distribution function (PDF) in the framework of large-momentum effective theory (LaMET) that uses the hybrid scheme with self-renormalization. We use ensembles generated by the MILC collaboration at lattice spacings $a=\{0.1207,0.0888,0.0582\}$ fm, with $N_f=2+1+1$ flavors of highly improved staggered quarks at sea pion mass of $M_{\pi}\approx 315$ MeV. We use clover-improved action for our valence quarks with nucleon boost momentum $P_z\approx 1.75$ GeV and high-statistics measurements for the LaMET matrix elements. We perform an extrapolation to the continuum limit and improve the handling of systematic errors using renormalization-group resummation (RGR) and leading-renormalon resummation (LRR). Our final nucleon helicity PDF is renormalized in the $\overline{\text{MS}}$ scheme at energy scale $\mu=2.0$ GeV. We compare our results with and without the two systematic improvements of RGR and LRR at each lattice spacing as well as the continuum limit, and we see that the application of RGR and LRR greatly reduces the systematic errors across the whole $x$ range. Our continuum results with both RGR and LRR show a small positive antiquark region for the nucleon helicity PDF as well as a change of as much as a factor of two in the central values compared to results with neither RGR or LRR. By contrast, the application of RGR and LRR only changes the central values by about 5\% in the quark region. We compare our lattice results with the global fits by the JAM, NNPDF and DSSV collaborations, and we observe some tension between our results.

High Energy Physics - Lattice

What problem does this paper attempt to address?

### Problems the paper attempts to solve

This paper aims to solve the problem of accurate calculation of the isovector helicity parton distribution function (PDF) of nucleons (protons and neutrons). Specifically, using the hybrid scheme and self - renormalization method in the framework of the large - momentum effective theory (LaMET), the paper calculates the isovector helicity parton distribution function of nucleons for the first time in the continuous limit.

### Main research contents

1. **Calculation methods**:
- Use the Hybrid Scheme with Self - Renormalization (HSR) to deal with the linear divergence and renormalon divergence in the matrix elements.
- Apply the Renormalization Group Resummation (RGR) and Leading Renormalon Resummation (LRR) methods to reduce systematic errors.
- Calculate at different lattice spacings and improve the accuracy of the results by extrapolating to the continuous limit.

2. **Experimental settings**:
- Use the ensembles with lattice spacings \(a=\{0.1207, 0.0888, 0.0582\}\) fm generated by the MILC collaboration, which contain \(N_f = 2 + 1+1\) flavors of highly improved staggered quarks, and the mass of the sea quarks is \(M_\pi\approx315\) MeV.
- Use valence quarks with Clover - improved action, and the nucleon momentum boost \(P_z\approx1.75\) GeV.

3. **Result analysis**:
- Compare the results at different lattice spacings and in the continuous limit before and after applying the RGR and LRR methods.
- It is found that applying the RGR and LRR methods significantly reduces systematic errors, especially in the anti - quark region.
- Compare with the global fitting results (such as those of the JAM, NNPDF and DSSV collaborations), and some inconsistencies are found.

### Key formulas

1. **Definition of matrix elements**:
\[
h_B(z, a)=\langle P_z|\bar{\psi}\left(-\frac{z}{2}\right)\Gamma W\left(-\frac{z}{2},\frac{z}{2}\right)\psi\left(\frac{z}{2}\right)|P_z\rangle
\]
where \(\Gamma = \gamma_z\gamma_5\) is the Dirac structure and \(W\) is the Wilson line connecting two space - time coordinates.

2. **Fitting form of the self - renormalization factor**:
\[
\ln\left(\frac{1}{Z(z, a)}\right)_{\text{NLO}}=k\frac{z}{a}\ln\left(\frac{a\Lambda_{\text{QCD}}}{\mu}\right)+g_1(z)+f_1(z)a+\frac{3C_F}{4\pi\beta_0}\ln\left(\frac{\ln\left(\frac{1}{a\Lambda_{\text{QCD}}}\right)}{\ln\left(\frac{\mu}{\Lambda_{\text{QCD}}}\right)}\right)+\ln\left(1 + d\ln\left(\frac{a\Lambda_{\text{QCD}}}{\mu}\right)\right)
\]

3. **Self - renormalization factor in the NLR case**:
\[
\ln\left(\frac{1}{Z(z, a)}\right)_{\text{NLR}}=\Delta I + k\frac{z}{a}\ln\left(\frac{a\Lambda_{\text{QCD}}}{\mu}\right)