Gluon-Proton Spin Correlations in QCD

• Gluon–proton spin correlations are defined by the interplay between gluon helicity and orbital angular momentum, fundamental to the proton spin sum rule in QCD.
• Experimental studies using longitudinal double-spin and transverse single-spin asymmetries at RHIC, LHC, and EIC quantify these correlations across various kinematic regimes.
• Global QCD analyses and advanced nonperturbative models combine DIS, SIDIS, and top decay measurements to provide quantitative constraints on gluon spin contributions.

Gluon–Proton Spin Correlations

Gluon–proton spin correlations encode the fundamental interplay between gluon spin degrees of freedom and the total angular momentum of the proton in quantum chromodynamics (QCD). This subject is central to resolving the proton spin decomposition problem, quantifying how gluon helicity and orbital angular momentum contribute to the proton’s total spin in both longitudinally and transversely polarized regimes. The field combines experimental measurements—especially at facilities such as RHIC and the LHC—with sophisticated global QCD analyses, nonperturbative modeling, and operator-based field theory. This article provides a comprehensive technical overview of the operational definitions, measurement strategies, phenomenological status, and theoretical frameworks underpinning gluon–proton spin correlations.

1. Operator Definitions and Spin Decomposition

The polarized gluon parton distribution function (PDF), denoted $\Delta g(x, Q^2) = g_+(x, Q^2) - g_-(x, Q^2)$, quantifies the net helicity of gluons with momentum fraction $x$ in a proton of given longitudinal spin. The correlator is defined through light-cone field strength operators: $$\Delta g(x, Q^2) = \frac{1}{2P^+}\int\frac{dz^-}{2\pi}e^{ixP^+z^-}\langle P, S_L| G^{+ \mu}(0)\,\tilde{G}_{\mu}^{~+}(z^-) |P, S_L\rangle~,$$ with $$\tilde{G}^{\mu\nu} \equiv \frac{1}{2}\epsilon^{\mu\nu\alpha\beta}G_{\alpha\beta}$$ and $S_L=\pm1$ the proton helicity.

The proton spin sum rule is formulated as: $$\frac{1}{2} = \frac{1}{2}\Delta\Sigma(Q^2) + \Delta G(Q^2) + L_q(Q^2) + L_g(Q^2)~,$$ where $\Delta\Sigma$ is the total quark (and antiquark) helicity contribution, $\Delta G(Q^2) = \int_0^1 dx\, \Delta g(x, Q^2)$ the gluon helicity, and $L_q$, $L_g$ the quark and gluon orbital angular momentum (OAM). The decomposition can be formulated in both the Ji (gauge-invariant, kinetic) and Jaffe–Manohar (canonical) schemes (Lin, 2024, Tan et al., 2023, Bhattacharya et al., 2022).

For transverse spin–momentum correlations, three-gluon (“trigluon”) twist–3 correlators and transverse-momentum-dependent (TMD) functions such as the gluon Sivers function $x$0 enter in field-theory definitions, capturing PT–odd correlations between the proton's spin and intrinsic gluon kinematics (Lewis, 2020, Acharya et al., 2021, Abdulameer et al., 2022).

2. Experimental Probes of Gluon Spin Correlations

2.1 Longitudinal Double–Spin Asymmetries at RHIC

The principal observable constraining the gluon helicity in the proton is the longitudinal double–spin asymmetry: $x$1 where $x$2 ($x$3) is the cross section for like (opposite) proton helicities. In leading-twist collinear factorization,

$x$4

RHIC measurements of inclusive jet, dijet, and pion-photon production at $x$5–$x$6 sample gluon $x$7 in $x$8, and have yielded nonzero $x$9, establishing a positive $$\Delta g(x, Q^2) = \frac{1}{2P^+}\int\frac{dz^-}{2\pi}e^{ixP^+z^-}\langle P, S_L| G^{+ \mu}(0)\,\tilde{G}_{\mu}^{~+}(z^-) |P, S_L\rangle~,$$0 for $$\Delta g(x, Q^2) = \frac{1}{2P^+}\int\frac{dz^-}{2\pi}e^{ixP^+z^-}\langle P, S_L| G^{+ \mu}(0)\,\tilde{G}_{\mu}^{~+}(z^-) |P, S_L\rangle~,$$1 (Surrow, 2013, Li, 2014, Walker, 2011, Lin, 2024).

2.2 Transverse Single–Spin Asymmetries and Three-Gluon Correlators

Transverse spin–momentum correlations are probed via the single–spin asymmetry $$\Delta g(x, Q^2) = \frac{1}{2P^+}\int\frac{dz^-}{2\pi}e^{ixP^+z^-}\langle P, S_L| G^{+ \mu}(0)\,\tilde{G}_{\mu}^{~+}(z^-) |P, S_L\rangle~,$$2 for inclusive and semi-inclusive channels sensitive to gluon subprocesses. Processes such as open–heavy–flavor production and midrapidity direct photon production in $$\Delta g(x, Q^2) = \frac{1}{2P^+}\int\frac{dz^-}{2\pi}e^{ixP^+z^-}\langle P, S_L| G^{+ \mu}(0)\,\tilde{G}_{\mu}^{~+}(z^-) |P, S_L\rangle~,$$3 collisions are predominantly gluon-initiated and isolate twist–3 three-gluon correlators and the gluon Sivers function. The null results for $$\Delta g(x, Q^2) = \frac{1}{2P^+}\int\frac{dz^-}{2\pi}e^{ixP^+z^-}\langle P, S_L| G^{+ \mu}(0)\,\tilde{G}_{\mu}^{~+}(z^-) |P, S_L\rangle~,$$4 at RHIC, with sensitivities down to $$\Delta g(x, Q^2) = \frac{1}{2P^+}\int\frac{dz^-}{2\pi}e^{ixP^+z^-}\langle P, S_L| G^{+ \mu}(0)\,\tilde{G}_{\mu}^{~+}(z^-) |P, S_L\rangle~,$$5, place strong constraints on gluon trigluon matrix elements and the first $$\Delta g(x, Q^2) = \frac{1}{2P^+}\int\frac{dz^-}{2\pi}e^{ixP^+z^-}\langle P, S_L| G^{+ \mu}(0)\,\tilde{G}_{\mu}^{~+}(z^-) |P, S_L\rangle~,$$6 moment of $$\Delta g(x, Q^2) = \frac{1}{2P^+}\int\frac{dz^-}{2\pi}e^{ixP^+z^-}\langle P, S_L| G^{+ \mu}(0)\,\tilde{G}_{\mu}^{~+}(z^-) |P, S_L\rangle~,$$7 (Abdulameer et al., 2022, Lewis, 2020, Acharya et al., 2021).

2.3 Heavy–Flavor and LHC Spin Correlations

At the LHC, top–antitop pair production is dominated by $$\Delta g(x, Q^2) = \frac{1}{2P^+}\int\frac{dz^-}{2\pi}e^{ixP^+z^-}\langle P, S_L| G^{+ \mu}(0)\,\tilde{G}_{\mu}^{~+}(z^-) |P, S_L\rangle~,$$8, and spin correlation observables in top decays provide access to $$\Delta g(x, Q^2) = \frac{1}{2P^+}\int\frac{dz^-}{2\pi}e^{ixP^+z^-}\langle P, S_L| G^{+ \mu}(0)\,\tilde{G}_{\mu}^{~+}(z^-) |P, S_L\rangle~,$$9 at high scales. The extraction relies on measuring double–spin and single–spin asymmetries in dilepton, lepton+jets, and fully hadronic channels, mapping out gluon polarizations in the relevant $$\tilde{G}^{\mu\nu} \equiv \frac{1}{2}\epsilon^{\mu\nu\alpha\beta}G_{\alpha\beta}$$0 and $$\tilde{G}^{\mu\nu} \equiv \frac{1}{2}\epsilon^{\mu\nu\alpha\beta}G_{\alpha\beta}$$1 domains (Goldstein et al., 2017).

2.4 Electron–Ion Collider and Exclusive/Dijet Channels

Exclusive dijet production in deep inelastic scattering (ep→e' jj p') at the planned EIC enables direct and differential measurement of the gluon orbital angular momentum via azimuthal double–spin or single–spin asymmetries, with the relevant modulations linked to CFFs of the gluon GTMDs and OAM distributions (Bhattacharya et al., 2022, Ji et al., 2016).

3. Global QCD Analysis and Gluon Helicity Extraction

Global fits incorporating DIS, SIDIS, and RHIC $$\tilde{G}^{\mu\nu} \equiv \frac{1}{2}\epsilon^{\mu\nu\alpha\beta}G_{\alpha\beta}$$2 data—exemplified by DSSV, NNPDFpol, and JAM frameworks—simultaneously constrain quark and gluon helicity distributions (Lin, 2024, Li, 2014, Binder et al., 2011). Typical results at $$\tilde{G}^{\mu\nu} \equiv \frac{1}{2}\epsilon^{\mu\nu\alpha\beta}G_{\alpha\beta}$$3 find:

• Truncated integral: $$\tilde{G}^{\mu\nu} \equiv \frac{1}{2}\epsilon^{\mu\nu\alpha\beta}G_{\alpha\beta}$$4,
• Extending to lower $$\tilde{G}^{\mu\nu} \equiv \frac{1}{2}\epsilon^{\mu\nu\alpha\beta}G_{\alpha\beta}$$5 via dijet topologies: $$\tilde{G}^{\mu\nu} \equiv \frac{1}{2}\epsilon^{\mu\nu\alpha\beta}G_{\alpha\beta}$$6.

These results establish the gluon spin contribution as a major component of the proton's angular momentum, comparable in magnitude to the total quark plus antiquark spin $$\tilde{G}^{\mu\nu} \equiv \frac{1}{2}\epsilon^{\mu\nu\alpha\beta}G_{\alpha\beta}$$7 (Surrow, 2013, Lin, 2024, Li, 2014).

The uncertainties at $$\tilde{G}^{\mu\nu} \equiv \frac{1}{2}\epsilon^{\mu\nu\alpha\beta}G_{\alpha\beta}$$8 remain large, underscoring the need for next-generation measurements at the EIC to fully map $$\tilde{G}^{\mu\nu} \equiv \frac{1}{2}\epsilon^{\mu\nu\alpha\beta}G_{\alpha\beta}$$9 to small $S_L=\pm1$0 and complete the spin puzzle.

4. Nonperturbative Structure and Spin–Orbit Correlations

4.1 GPDs, TMDs, and Wigner Distributions

Theoretical modeling via light-front quark–gluon spectator models, AdS/QCD, and BLFQ Hamiltonians has yielded explicit expressions for gluon GPDs ($S_L=\pm1$1, $S_L=\pm1$2, $S_L=\pm1$3), TMDs, and GTMDs. The GTMD $S_L=\pm1$4 encodes the gluon OAM density as the $S_L=\pm1$5 moment, while the GPD $S_L=\pm1$6 is directly linked to OAM through Ji's sum rule: $S_L=\pm1$7 Phase-space (Wigner) distributions and GTMDs additionally reveal spin–orbit correlations, quantified by canonical OAM and spin–OAM correlators such as $S_L=\pm1$8 (Chakrabarti et al., 17 Sep 2025, Tan et al., 2023, Choudhary et al., 2024).

4.2 Quantum Entanglement Measures

Light-front Hamiltonian approaches (BLFQ) demonstrate that the inclusion of dynamical gluons in Fock sector expansions enhances spin–spin entanglement entropy among proton constituents, shifting the maximum entanglement to mid-$S_L=\pm1$9 and directly connecting entanglement entropy to measured helicity distributions via binary entropy relations. These quantum-informational diagnostics provide new avenues for probing multi-parton correlations in future experiments (Qian et al., 2024).

4.3 Spin–Orbit Dynamics and Model Uncertainties

Model calculations indicate that the gluon spin and OAM can be sizeable and often opposite in sign, with kinetic OAM $$\frac{1}{2} = \frac{1}{2}\Delta\Sigma(Q^2) + \Delta G(Q^2) + L_q(Q^2) + L_g(Q^2)~,$$0 typically negative for positive $$\frac{1}{2} = \frac{1}{2}\Delta\Sigma(Q^2) + \Delta G(Q^2) + L_q(Q^2) + L_g(Q^2)~,$$1 (Tan et al., 2023, Choudhary et al., 2024). These findings are robust across soft-wall AdS/QCD, holographic light-front models, and statistical parton frameworks (Bourrely et al., 2014, Gurjar et al., 2022), with extracted $$\frac{1}{2} = \frac{1}{2}\Delta\Sigma(Q^2) + \Delta G(Q^2) + L_q(Q^2) + L_g(Q^2)~,$$2 compatible with both lattice QCD and phenomenological fits.

5. Constraints, Implications, and Future Prospects

5.1 Summary of Present Knowledge

A cohesive picture emerges: gluon helicity $$\frac{1}{2} = \frac{1}{2}\Delta\Sigma(Q^2) + \Delta G(Q^2) + L_q(Q^2) + L_g(Q^2)~,$$3 is established as a positive, non-negligible component of proton spin for $$\frac{1}{2} = \frac{1}{2}\Delta\Sigma(Q^2) + \Delta G(Q^2) + L_q(Q^2) + L_g(Q^2)~,$$4, with substantial theoretical and experimental evidence (Surrow, 2013, Li, 2014, Lin, 2024). The corresponding OAM $$\frac{1}{2} = \frac{1}{2}\Delta\Sigma(Q^2) + \Delta G(Q^2) + L_q(Q^2) + L_g(Q^2)~,$$5 is inferred to be negative and large enough to nearly balance $$\frac{1}{2} = \frac{1}{2}\Delta\Sigma(Q^2) + \Delta G(Q^2) + L_q(Q^2) + L_g(Q^2)~,$$6 in certain models, conforming to the proton spin sum rule across multiple theoretical frameworks (Zhu et al., 2015, Tan et al., 2023, Choudhary et al., 2024). Constraints on transverse spin–momentum correlations, including the gluon Sivers function and trigluon twist-3 correlators, are strong: current data indicate minimal transverse spin–momentum correlations for gluons at moderate $$\frac{1}{2} = \frac{1}{2}\Delta\Sigma(Q^2) + \Delta G(Q^2) + L_q(Q^2) + L_g(Q^2)~,$$7 (Lewis, 2020, Abdulameer et al., 2022, Acharya et al., 2021).

5.2 Outstanding Questions and the Role of the EIC

The central open issue is the behavior of $$\frac{1}{2} = \frac{1}{2}\Delta\Sigma(Q^2) + \Delta G(Q^2) + L_q(Q^2) + L_g(Q^2)~,$$8 for $$\frac{1}{2} = \frac{1}{2}\Delta\Sigma(Q^2) + \Delta G(Q^2) + L_q(Q^2) + L_g(Q^2)~,$$9 and the detailed mapping of $\Delta\Sigma$0 and spin–orbit correlations at small $\Delta\Sigma$1. The upcoming Electron–Ion Collider is designed to address these gaps using exclusive and semi–inclusive channels sensitive to gluon GTMDs and OAM distributions, exploiting azimuthal spin asymmetries and reconstructing Wigner phase-space distributions (Bhattacharya et al., 2022, Ji et al., 2016). Such measurements will enable direct experimental determination of the elusive gluon OAM and provide a definitive test of the QCD spin sum rule.

6. Tables: Key Experimental Results and Theoretical Constructs

Table 1. RHIC Constraints on Integrated Gluon Helicity

$\Delta\Sigma$2 Range	$\Delta\Sigma$3 Estimate	Source / Analysis
$\Delta\Sigma$4	$\Delta\Sigma$5	DSSV14 Global Fits (Li, 2014, Lin, 2024)
$\Delta\Sigma$6	$\Delta\Sigma$7	STAR (incl. forward dijets) (Lin, 2024)

Table 2. Theoretical Model Results for Gluon Angular Momentum

Model / Approach	$\Delta\Sigma$8	$\Delta\Sigma$9	Reference
Light-cone spectator	0.19	$\Delta G(Q^2) = \int_0^1 dx\, \Delta g(x, Q^2)$0	(Tan et al., 2023)
Light-front AdS/QCD	$\Delta G(Q^2) = \int_0^1 dx\, \Delta g(x, Q^2)$1	—	(Gurjar et al., 2022)
LF spectator (DIS2024)	0.48	$\Delta G(Q^2) = \int_0^1 dx\, \Delta g(x, Q^2)$2	(Choudhary et al., 2024)
BLFQ (spin entropy)	—	—	(Qian et al., 2024)

These results reinforce the central finding that gluon spin and OAM are dynamically entangled components of the proton's spin budget, with both experimental and theoretical efforts converging towards a detailed, quantitative understanding of gluon–proton spin correlations.