We perform a phenomenological study of helicity-dependent parton distribution functions (PDFs) using small-<math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>x</mi></math> helicity evolution equations, incorporating for the first time single-inclusive jet production data in polarized proton-proton (<math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>p</mi><mi>p</mi></math>) scattering at parton momentum fractions <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>x</mi><mo>&lt;</mo><mn>0.1</mn></math>. We also simultaneously include double-longitudinal spin asymmetries in inclusive and semi-inclusive deep-inelastic scattering probing <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>x</mi><mo>&lt;</mo><mn>0.1</mn></math>. Employing the polarized small-<math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>x</mi></math> pure-glue calculation of <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>p</mi><mi>p</mi><mo stretchy="false">→</mo><mi>g</mi><mi>X</mi></math> for the jet production cross section, we modify the large-<math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>c</mi></mrow></msub><mo>&amp;</mo><msub><mrow><mi>N</mi></mrow><mrow><mi>f</mi></mrow></msub></mrow></math> KPS-CTT evolution equations by setting <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msub><mi>N</mi><mi>f</mi></msub><mo>=</mo><mn>0</mn></math> to replicate the large-<math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msub><mi>N</mi><mi>c</mi></msub></math> (pure-glue) limit, while retaining external quark flavors for the spinor field operators. We find that the <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>p</mi><mi>p</mi></math> data have a considerable impact on the helicity PDFs at small <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>x</mi></math>, reducing their uncertainties and leading to a total quark and gluon helicity in the proton for <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>x</mi><mo>&lt;</mo><mn>0.1</mn></math> of <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mo>−</mo><mn>0.04</mn><mo>±</mo><mn>0.23</mn></math>. Combining our analysis with a recent JAM helicity PDF analysis of the world polarized data, which includes <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>x</mi><mo>&gt;</mo><mn>0.1</mn></math>, we find a total quark and gluon helicity contribution for <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>x</mi><mo>&gt;</mo><msup><mn>10</mn><mrow><mo>−</mo><mn>7</mn></mrow></msup></math> of between 0.02 and 0.51.

United States

Daniel Adamiak

Nicholas Baldonado

Yuri V. Kovchegov

Ming Li

W. Melnitchouk

Daniel Pitonyak

Nobuo Sato

Matthew D. Sievert

Andrey Tarasov

Thailand

Yossathorn Tawabutr

First study of polarized proton-proton scattering with small-<math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>x</mi></math> helicity evolution

We calculate the perturbative corrections in the strong coupling to the unpolarized quark transverse-momentum dependent parton distribution function (TMDPDF) operator within a background-field framework, extending the approach of Mukherjee et al. [<a href="http://dx.doi.org/10.1103/PhysRevD.109.034035">Phys. Rev. D 109, 034035 (2024)</a>]. We focus on ensuring gauge invariance, identifying two key components needed: a gauge-invariant TMDPDF operator with a transverse gauge link at spatial infinity, and accounting of the equations of motion (EoM) of the background fields. We go beyond next-to-leading order in strong coupling expansion, considering multiple interactions with the background field at all orders of strong coupling. By examining the interplay between quark and gluon contributions, we show that spurious singularities, proportional to EoM, can be misinterpreted as genuine divergences in quantum chromodynamics factorization unless properly identified and removed.

Swagato Mukherjee

Vladimir V. Skokov

Shaswat Tiwari

Perturbative corrections to quark TMDPDFs in the background-field method: Gauge invariance, equations of motion, and multiple interactions

We extend our prior results on the worldline computation of the axial vector-vector-vector (AVV) triangle anomaly in polarized deeply inelastic scattering (DIS) to the finite mass case by computing in addition the pseudoscalar-vector-vector (PVV) triangle graph. For the well-studied QED case, we show explicitly how the off-forward AVV pole exactly cancels an identical PVV pole. We then demonstrate the dramatic difference in QCD due to the chiral condensate, which qualitatively modifies anomalous Ward identities. As in the massless case, the anomaly pole in QCD is canceled by the dynamics of a primordial isosinglet pseudoscalar <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mover accent="true"><mi>η</mi><mo stretchy="false">¯</mo></mover></math>-meson, whose Wess-Zumino-Witten coupling to the topological charge density shifts the pole to the physical <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msup><mi>η</mi><mo>′</mo></msup></math> mass, with the finite quark mass contribution differing by <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>O</mi><mo stretchy="false">(</mo><mn>10</mn><mo>%</mo><mo stretchy="false">)</mo></math> from the Witten-Veneziano formula. We obtain a compact analytic expression for the finite mass corrections to Shore and Veneziano’s result that the proton’s net quark helicity <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi mathvariant="normal">Δ</mi><mi mathvariant="normal">Σ</mi><mo>∝</mo><msqrt><mrow><msubsup><mi>χ</mi><mi>QCD</mi><mo>′</mo></msubsup><msub><mo stretchy="false">|</mo><mrow><mi>m</mi><mo>=</mo><mn>0</mn></mrow></msub><mo stretchy="false">(</mo><mn>0</mn><mo stretchy="false">)</mo></mrow></msqrt></math>, the forward slope of the topological susceptibility in the chiral limit, and show they are of the order of a few percentages. Our prior prediction that the polarized DIS structure function <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><msub><mi>g</mi><mn>1</mn></msub></math> is quenched by sphaleronlike topological transitions at small <math xmlns="http://www.w3.org/1998/Math/MathML" display="inline"><mi>x</mi></math> is unaffected by quark mass effects. Our results illustrate how worldline computations of anomalous processes, in synergy with lattice computations and nonet chiral perturbation theory, can uncover novel nonperturbative features of QCD at the Electron-Ion Collider.

Raju Venugopalan

Role of the chiral anomaly in polarized deeply inelastic scattering. III. Wess-Zumino-Witten contributions and chiral Ward identities for finite quark mass

Finland

Florian Cougoulic

Yuri Kovchegov

Erratum to: Quark and gluon helicity evolution at small x: revised and updated

This paper introduces a transverse-momentum dependent (TMD) factorization scheme designed to unify both large and small Bjorken-x regimes. We compute the next-to-leading order quantum chromodynamics (QCD) corrections to the gluon TMD operator for an unpolarized hadron within this proposed scheme. This leads to the emergence of a new TMD evolution, incorporating those in transverse momentum, rapidity, and Bjorken-x. When matched to the collinear factorization scheme, our factorization scheme faithfully reproduces the well-established Dokshitzer-Gribov-Lipatov-Altarelli-Parisi and Collins-Soper-Sterman (CSS) evolutions. Conversely, matching with high-energy factorization not only yields the Balitsky-Fadin-Kuraev-Lipatov evolution but also reveals distinctive signatures of CSS logarithms. The development of this novel TMD factorization scheme, capable of seamlessly reconciling disparate Bjorken-x regimes and faithfully reproducing established QCD evolution equations, has the potential to significantly advance our comprehension of high-energy processes and three-dimensional parton structures of hadrons.

Unified description of DGLAP, CSS, and BFKL evolution: TMD factorization bridging large and small <math><mrow><mi>x</mi></mrow></math>

We analyze the world polarized deep-inelastic scattering (DIS) and semi-inclusive DIS (SIDIS) data at low values of <math><mi>x</mi><mo>&lt;</mo><mn>0.1</mn></math>, using small-<math><mi>x</mi></math> evolution equations for the flavor singlet and nonsinglet helicity parton distribution functions (hPDFs), which resum all powers of both <math><msub><mi>α</mi><mi>s</mi></msub><msup><mi>ln</mi><mn>2</mn></msup><mo>(</mo><mn>1</mn><mo>/</mo><mi>x</mi><mo>)</mo></math> and <math><msub><mi>α</mi><mi>s</mi></msub><mi>ln</mi><mo>(</mo><mn>1</mn><mo>/</mo><mi>x</mi><mo>)</mo><mi>ln</mi><mo>(</mo><msup><mi>Q</mi><mn>2</mn></msup><mo>/</mo><msubsup><mi>Q</mi><mn>0</mn><mn>2</mn></msubsup><mo>)</mo></math>, with <math><msub><mi>α</mi><mi>s</mi></msub></math> being the strong coupling constant. The hPDFs for quarks, antiquarks, and gluons are extracted and evolved to lower values of <math><mi>x</mi></math> to make predictions for the future Electron-Ion Collider (EIC). We improve on our earlier work by employing the more realistic large-<math><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi>c</mi></mrow></msub><mo>&amp;</mo><msub><mrow><mi>N</mi></mrow><mrow><mi>f</mi></mrow></msub></mrow></math> limit of the revised small-<math><mi>x</mi></math> helicity evolution, and by incorporating running coupling corrections along with SIDIS data into the fit. We find an anticorrelation between the signs of the gluon and <math><mi>C</mi></math>-even quark hPDFs, as well as the <math><msub><mi>g</mi><mn>1</mn></msub></math> structure function. While the existing low-<math><mi>x</mi></math> polarized DIS and SIDIS data are insufficient to constrain the initial conditions for the polarized dipole amplitudes in the helicity evolution equations, future EIC data will allow more precise predictions for hPDFs and the <math><msub><mi>g</mi><mn>1</mn></msub></math> structure function for <math><mi>x</mi></math> values beyond those probed at the EIC. Using the obtained hPDFs, we discuss the contributions to the proton spin from quark and gluon spins at small <math><mi>x</mi></math>.

Lebanon

Global analysis of polarized DIS and SIDIS data with improved small-<math><mi>x</mi></math> helicity evolution

In [1], we developed a first-quantized worldline formalism for all-order computations of amplitudes in QED. In particular, we demonstrated in this framework an all-order proof of the infrared safety of the Faddeev-Kulish (FK) <math><mi>S</mi></math>-matrix for virtual exchanges in the scattering of charged fermions. In this work, we extend the worldline formalism for both the Dyson and FK <math><mi>S</mi></math>-matrix to consider further the emission and absorption of arbitrary numbers of photons. We show how Low’s theorem follows in this framework and derive Weinberg’s theorem for the exponentiation of IR divergences. In particular, we extend our all-order proof of the IR safety of the FK <math><mi>S</mi></math>-matrix to both virtual exchanges and real photon emissions. We argue that the worldline approach leads to a modern Wilsonian interpretation of the IR safety of the FK <math><mi>S</mi></math>-matrix and provides a novel template for the treatment of IR divergences in real-time problems. Using Grassmannian integration methods, we derive a simple and powerful result for N-th rank vacuum polarization tensors. Applications of these methods will be discussed in follow-up work.

Spain

Xabier Feal

QED as a many-body theory of worldlines. II. All-order <math><mi>S</mi></math>-matrix formalism

We discuss a reformulation of QED in which matter and gauge fields are integrated out explicitly, resulting in a many-body Lorentz covariant theory of <math><mn>0</mn><mo>+</mo><mn>1</mn></math>-dimensional worldlines describing superpairs of spinning charges interacting through Lorentz forces. This provides a powerful, string inspired definition of amplitudes to all-loop orders. In particular, one obtains a more general formulation of Wilson loops and lines, with exponentiated dynamical fields and spin precession contributions, and worldline contour averages exactly defined through first quantized path integrals. We discuss in detail the attractive features of this formalism for high order perturbative computations. We show that worldline S-matrix elements, to all-loop orders in perturbation theory, can be constructed to be manifestly free of soft singularities, with infrared (IR) divergences captured and removed by endpoint photon exchanges at infinity that are equivalent to the soft coherent dressings of the Dyson S-matrix proposed by Faddeev and Kulish. We discuss these IR structures and make connections with soft theorems, the Abelian exponentiation of IR divergences and cusp anomalous dimensions. Follow-up papers will discuss the efficient computation of cusp anomalous dimensions and universal features of soft theorems.

QED as a many-body theory of worldlines: General formalism and infrared structure

We revisit the problem of small Bjorken-x evolution of the gluon and flavor-singlet quark helicity distributions in the shock wave (s-channel) formalism. Earlier works on the subject in the same framework [1–3] resulted in an evolution equation for the gluon field-strength F 12 and quark “axial current” <math> <mover> <mi>ψ</mi> <mo>¯</mo> </mover> <mi>γ</mi> </math> $$ \overline{\psi}\gamma $$ + γ 5 ψ operators (sandwiched between the appropriate light-cone Wilson lines) in the double-logarithmic approximation (summing powers of α s ln2(1/x) with α s the strong coupling constant). In this work, we observe that an important mixing of the above operators with another gluon operator, <math> <mmultiscripts> <mrow></mrow> <mprescripts></mprescripts> <mi>D</mi> <mrow> <mo>←</mo> <mi>i</mi> </mrow> </mmultiscripts> </math> $$ {}_D{}^{\leftarrow i} $$ D i , also sandwiched between the light-cone Wilson lines (with the repeated transverse index i = 1, 2 summed over), was missing in the previous works [1–3]. This operator has the physical meaning of the sub-eikonal (covariant) phase: its contribution to helicity evolution is shown to be proportional to another sub-eikonal operator, D i − <math> <mmultiscripts> <mrow></mrow> <mprescripts></mprescripts> <mi>D</mi> <mrow> <mo>←</mo> <mi>i</mi> </mrow> </mmultiscripts> </math> $$ {}_D{}^{\leftarrow i} $$ , which is related to the Jaffe-Manohar polarized gluon distribution [4]. In this work we include this new operator into small-x helicity evolution, and construct novel evolution equations mixing all three operators (D i − <math> <mmultiscripts> <mrow></mrow> <mprescripts></mprescripts> <mi>D</mi> <mrow> <mo>←</mo> <mi>i</mi> </mrow> </mmultiscripts> </math> $$ {}_D{}^{\leftarrow i} $$ , F 12, and <math> <mover> <mi>ψ</mi> <mo>¯</mo> </mover> <mi>γ</mi> </math> $$ \overline{\psi}\gamma $$ + γ 5 ψ), generalizing the results of [1–3]. We also construct closed double-logarithmic evolution equations in the large-N c and large-N c &amp;N f limits, with N c and N f the numbers of quark colors and flavors, respectively. Solving the large-N c equations numerically we obtain the following small-x asymptotics of the quark and gluon helicity distributions ∆Σ and ∆G, along with the g 1 structure function, <math> <mo>∆</mo> <mi>Σ</mi> <mfenced> <mi>x</mi> <msup> <mi>Q</mi> <mn>2</mn> </msup> </mfenced> <mo>∼</mo> <mo>∆</mo> <mi>G</mi> <mfenced> <mi>x</mi> <msup> <mi>Q</mi> <mn>2</mn> </msup> </mfenced> <mo>∼</mo> <msub> <mi>g</mi> <mn>1</mn> </msub> <mfenced> <mi>x</mi> <msup> <mi>Q</mi> <mn>2</mn> </msup> </mfenced> <mo>∼</mo> <msup> <mfenced> <mfrac> <mn>1</mn> <mi>x</mi> </mfrac> </mfenced> <mrow> <mn>3.66</mn> <msqrt> <mfrac> <mrow> <msub> <mi>α</mi> <mi>s</mi> </msub> <msub> <mi>N</mi> <mi>c</mi> </msub> </mrow> <mrow> <mn>2</mn> <mi>π</mi> </mrow> </mfrac> </msqrt> </mrow> </msup> </math> $$ \Delta \Sigma \left(x,{Q}^2\right)\sim \Delta G\left(x,{Q}^2\right)\sim {g}_1\left(x,{Q}^2\right)\sim {\left(\frac{1}{x}\right)}^{3.66\sqrt{\frac{\alpha_s{N}_c}{2\pi }}} $$

Quark and gluon helicity evolution at small x: revised and updated

We calculate single-logarithmic corrections to the small-x flavor-singlet helicity evolution equations derived recently [1–3] in the double-logarithmic approximation. The new single-logarithmic part of the evolution kernel sums up powers of α  s  ln(1/x), which are an important correction to the dominant powers of α  s  ln2(1/x) summed up by the double-logarithmic kernel from [1–3] at small values of Bjorken x and with α  s  the strong coupling constant. The single-logarithmic terms arise separately from either the longitudinal or transverse momentum integrals. Consequently, the evolution equations we derive employing the light-cone perturbation theory simultaneously include the small-x evolution kernel and the leading-order polarized DGLAP splitting functions. We further enhance the equations by calculating the running coupling corrections to the kernel.