Leading-Twist Quark Wigner Distributions

• Leading-twist quark Wigner distributions are light-front phase-space maps that combine GTMD, TMD, and GPD data to reveal quark spin and orbital dynamics.
• They use transverse Fourier transforms to translate momentum transfer into impact-parameter space, enabling a comprehensive nucleon tomography.
• These distributions provide a framework to extract quark orbital angular momentum and analyze spin-spin and spin-orbit correlations through multipole decompositions.

Leading-twist quark Wigner distributions are light-front phase-space distributions of the nucleon that depend on the longitudinal momentum fraction $x$, the quark transverse momentum $\boldsymbol{k}_\perp$, and the transverse impact parameter $\boldsymbol{b}_\perp$. They are obtained as transverse Fourier transforms of generalized transverse-momentum dependent distributions (GTMDs), and thereby combine in a single picture the information separately contained in transverse-momentum dependent distributions (TMDs) and generalized parton distributions (GPDs). At leading twist, they form a complete set of real-valued but non-positive-definite quasi-distributions, with 16 independent polarization structures for a spin-$1/2$ target. Their main theoretical significance is that they expose spin-orbit and spin-spin correlations inaccessible in pure GPD or TMD projections, and they provide a direct phase-space route to quark orbital angular momentum (OAM) (Lorce et al., 2011, Lorce et al., 2012).

1. Operator definition and light-front kinematics

The standard construction starts from the off-forward quark correlator at fixed light-front time $z^+=0$,

$$W^{[\Gamma]q}_{\Lambda'\Lambda}(P,x,\boldsymbol{k}_\perp,\Delta,n) =\frac{1}{2}\int\frac{dz^-\,d^2\boldsymbol{z}_\perp}{(2\pi)^3}\, e^{ixP^+z^- - i\boldsymbol{k}_\perp\cdot\boldsymbol{z}_\perp} \,\langle p',\Lambda'|\,\bar\psi\!\left(-\tfrac{z}{2}\right)\Gamma\,\mathcal W\,\psi\!\left(\tfrac{z}{2}\right)\,|p,\Lambda\rangle\Big|_{z^+=0},$$

with $P=(p'+p)/2$, $\Delta=p'-p$, and $x=k^+/P^+$. The Wilson line $\mathcal W$ ensures color gauge invariance and is specified by a path choice, typically straight or staple-like (Lorce et al., 2012).

The five-dimensional Wigner distribution is then defined by Fourier transforming the transverse momentum transfer to impact-parameter space,

$\boldsymbol{k}_\perp$0

In the infinite-momentum frame, or equivalently the Drell–Yan–West frame with $\boldsymbol{k}_\perp$1, these objects admit a particularly clean interpretation because transverse boosts form a Galilean subgroup, allowing transversely localized target states without additional relativistic ambiguities (Lorce et al., 2011).

The leading-twist quark bilinears are

$\boldsymbol{k}_\perp$2

corresponding respectively to unpolarized, longitudinally polarized, and transversely polarized quark probes. Although GTMDs are generally complex, their two-dimensional Fourier transforms are real, consistent with a phase-space interpretation. The resulting Wigner distributions are nevertheless quasi-probabilistic rather than genuine probabilities, because the underlying position and momentum operators do not commute (Lorce et al., 2012).

2. Leading-twist classification and polarization content

For a spin-$\boldsymbol{k}_\perp$3 target, leading twist contains 16 independent Wigner distributions covering all target/quark polarization combinations. In the longitudinal sector, a convenient decomposition is

$\boldsymbol{k}_\perp$4

with $\boldsymbol{k}_\perp$5 denoting target and quark longitudinal polarizations (Lorce et al., 2011).

In terms of GTMDs, the longitudinal channels are

$\boldsymbol{k}_\perp$6

This makes explicit that $\boldsymbol{k}_\perp$7 is controlled by $\boldsymbol{k}_\perp$8 and $\boldsymbol{k}_\perp$9 by $\boldsymbol{b}_\perp$0 (Lorce et al., 2011).

A recurring feature of the literature is that notation is not fully uniform. Some papers use $\boldsymbol{b}_\perp$1 for unpolarized quarks in a longitudinally polarized target, while others note that some literature uses $\boldsymbol{b}_\perp$2 for the same physical configuration; the operator definition removes the ambiguity (Pasquini et al., 2011).

Beyond the longitudinal sector, the full leading-twist set can be organized by basic multipoles in transverse momentum and transverse position space. The complete T-even and T-odd multipole decomposition shows that each polarization channel carries a definite spin-spin or spin-orbit content: monopoles encode unpolarized or helicity densities, dipoles encode net transverse or azimuthal flows, and quadrupoles encode transversity or pretzelous-type structures (Lorcé et al., 2015). This multipole language is especially useful because it simultaneously clarifies which structures survive in TMD or GPD projections and which require the full phase-space description.

3. Relations to GTMDs, TMDs, and GPDs

Wigner distributions are the impact-parameter representation of GTMDs, and GTMDs are their purely momentum-space counterparts. The standard reduction relations are

$\boldsymbol{b}_\perp$3

and

$\boldsymbol{b}_\perp$4

Thus, integrating over $\boldsymbol{b}_\perp$5 yields the forward limit of GTMDs, i.e. TMDs, while integrating over $\boldsymbol{b}_\perp$6 yields impact-parameter GPDs (Lorce et al., 2012).

Some GTMDs are genuine “mother distributions” in the sense that both their TMD and GPD limits survive. In particular,

$\boldsymbol{b}_\perp$7

while their $\boldsymbol{b}_\perp$8 integrals generate $\boldsymbol{b}_\perp$9 and $1/2$0. By contrast, integrating over either $1/2$1 or $1/2$2 kills $1/2$3 and $1/2$4, showing that $1/2$5 and $1/2$6 encode correlations absent from any leading-twist TMD or GPD alone (Lorce et al., 2011).

The symmetry constraints are correspondingly strong. Hermiticity enforces the reality of the Wigner distributions, parity fixes the allowed tensor structures, and time reversal interacts with the gauge-link geometry. With a straight lightlike Wilson line, naively T-odd components vanish; with a staple-like Wilson line, naively time-reversal odd structures are allowed and encode effects such as the Sivers and Boer–Mulders functions (Lorce et al., 2012).

A central conceptual limitation concerns observability. The leading-twist GTMD combination associated with the unpolarized-quark/longitudinally polarized-proton configuration is the $1/2$7 term,

$1/2$8

which is proportional to $1/2$9. Parity constraints imply that this structure is parity-odd and nonzero only for imaginary values of the helicity amplitudes. In coplanar tree-level kinematics it therefore cannot contribute at leading twist to deeply virtual Compton scattering (DVCS) without final-state interactions; moreover, integrating $z^+=0$0 over $z^+=0$1 gives zero, so it drops out of GPDs (Courtoy et al., 2013). This is the main reason why OAM-sensitive observables in DVCS appear at twist three rather than directly through leading-twist $z^+=0$2.

4. Orbital angular momentum and competing definitions

The most prominent application of leading-twist quark Wigner distributions is the extraction of the longitudinal quark OAM from the phase-space average of $z^+=0$3,

$z^+=0$4

In impact-parameter language, this projects precisely the $z^+=0$5 structure (Lorce et al., 2012).

For a straight gauge link one obtains the compact relation

$z^+=0$6

The path choice of the Wilson line determines the OAM notion. A straight lightlike Wilson line yields kinetic OAM in the sense of Ji’s decomposition, whereas a staple-like Wilson line yields canonical OAM in the sense of Jaffe–Manohar (Lorce et al., 2012). Earlier model analyses also emphasized that in the absence of explicit gauge-field degrees of freedom these distinctions can collapse for the flavor-summed total even while flavor-separated values differ (Lorce et al., 2011).

The Ji sum rule gives the kinetic OAM through twist-2 GPDs,

$z^+=0$7

The same quantity also appears in twist-three form. Using the twist-three parametrization of Meissner–Metz–Schlegel, Lorcé and Pasquini derived

$z^+=0$8

which is equivalent to

$z^+=0$9

These relations connect longitudinal OAM to the Penttinen–Polyakov–Shuvaev–Strikman sum rule and clarify why twist-three GPDs are central to observability (Lorce et al., 2012).

A frequently discussed proxy is the pretzelosity-based expression

$$W^{[\Gamma]q}_{\Lambda'\Lambda}(P,x,\boldsymbol{k}_\perp,\Delta,n) =\frac{1}{2}\int\frac{dz^-\,d^2\boldsymbol{z}_\perp}{(2\pi)^3}\, e^{ixP^+z^- - i\boldsymbol{k}_\perp\cdot\boldsymbol{z}_\perp} \,\langle p',\Lambda'|\,\bar\psi\!\left(-\tfrac{z}{2}\right)\Gamma\,\mathcal W\,\psi\!\left(\tfrac{z}{2}\right)\,|p,\Lambda\rangle\Big|_{z^+=0},$$0

but the literature explicitly stresses that no rigorous expression of OAM in terms of leading-twist TMDs is known (Lorce et al., 2012).

The observability analysis sharpens this distinction further. In DVCS, the OAM-relevant helicity structure is accessed through twist-three GTMDs and GPDs, especially $$W^{[\Gamma]q}_{\Lambda'\Lambda}(P,x,\boldsymbol{k}_\perp,\Delta,n) =\frac{1}{2}\int\frac{dz^-\,d^2\boldsymbol{z}_\perp}{(2\pi)^3}\, e^{ixP^+z^- - i\boldsymbol{k}_\perp\cdot\boldsymbol{z}_\perp} \,\langle p',\Lambda'|\,\bar\psi\!\left(-\tfrac{z}{2}\right)\Gamma\,\mathcal W\,\psi\!\left(\tfrac{z}{2}\right)\,|p,\Lambda\rangle\Big|_{z^+=0},$$1 or $$W^{[\Gamma]q}_{\Lambda'\Lambda}(P,x,\boldsymbol{k}_\perp,\Delta,n) =\frac{1}{2}\int\frac{dz^-\,d^2\boldsymbol{z}_\perp}{(2\pi)^3}\, e^{ixP^+z^- - i\boldsymbol{k}_\perp\cdot\boldsymbol{z}_\perp} \,\langle p',\Lambda'|\,\bar\psi\!\left(-\tfrac{z}{2}\right)\Gamma\,\mathcal W\,\psi\!\left(\tfrac{z}{2}\right)\,|p,\Lambda\rangle\Big|_{z^+=0},$$2, rather than through leading-twist $$W^{[\Gamma]q}_{\Lambda'\Lambda}(P,x,\boldsymbol{k}_\perp,\Delta,n) =\frac{1}{2}\int\frac{dz^-\,d^2\boldsymbol{z}_\perp}{(2\pi)^3}\, e^{ixP^+z^- - i\boldsymbol{k}_\perp\cdot\boldsymbol{z}_\perp} \,\langle p',\Lambda'|\,\bar\psi\!\left(-\tfrac{z}{2}\right)\Gamma\,\mathcal W\,\psi\!\left(\tfrac{z}{2}\right)\,|p,\Lambda\rangle\Big|_{z^+=0},$$3. In the Wandzura–Wilczek approximation,

$$W^{[\Gamma]q}_{\Lambda'\Lambda}(P,x,\boldsymbol{k}_\perp,\Delta,n) =\frac{1}{2}\int\frac{dz^-\,d^2\boldsymbol{z}_\perp}{(2\pi)^3}\, e^{ixP^+z^- - i\boldsymbol{k}_\perp\cdot\boldsymbol{z}_\perp} \,\langle p',\Lambda'|\,\bar\psi\!\left(-\tfrac{z}{2}\right)\Gamma\,\mathcal W\,\psi\!\left(\tfrac{z}{2}\right)\,|p,\Lambda\rangle\Big|_{z^+=0},$$4

while beyond Wandzura–Wilczek the difference between kinetic and canonical OAM is attributed to final-state interactions and gauge-link structure (Courtoy et al., 2013).

5. Model calculations and tomographic patterns

The earliest detailed pictures were obtained in valence-quark light-front models. In the light-front constituent quark model (LFCQM), the $$W^{[\Gamma]q}_{\Lambda'\Lambda}(P,x,\boldsymbol{k}_\perp,\Delta,n) =\frac{1}{2}\int\frac{dz^-\,d^2\boldsymbol{z}_\perp}{(2\pi)^3}\, e^{ixP^+z^- - i\boldsymbol{k}_\perp\cdot\boldsymbol{z}_\perp} \,\langle p',\Lambda'|\,\bar\psi\!\left(-\tfrac{z}{2}\right)\Gamma\,\mathcal W\,\psi\!\left(\tfrac{z}{2}\right)\,|p,\Lambda\rangle\Big|_{z^+=0},$$5-integrated distribution of unpolarized quarks in a longitudinally polarized nucleon exhibits a clear dipole pattern in impact-parameter space at fixed transverse momentum, for example at $$W^{[\Gamma]q}_{\Lambda'\Lambda}(P,x,\boldsymbol{k}_\perp,\Delta,n) =\frac{1}{2}\int\frac{dz^-\,d^2\boldsymbol{z}_\perp}{(2\pi)^3}\, e^{ixP^+z^- - i\boldsymbol{k}_\perp\cdot\boldsymbol{z}_\perp} \,\langle p',\Lambda'|\,\bar\psi\!\left(-\tfrac{z}{2}\right)\Gamma\,\mathcal W\,\psi\!\left(\tfrac{z}{2}\right)\,|p,\Lambda\rangle\Big|_{z^+=0},$$6 with $$W^{[\Gamma]q}_{\Lambda'\Lambda}(P,x,\boldsymbol{k}_\perp,\Delta,n) =\frac{1}{2}\int\frac{dz^-\,d^2\boldsymbol{z}_\perp}{(2\pi)^3}\, e^{ixP^+z^- - i\boldsymbol{k}_\perp\cdot\boldsymbol{z}_\perp} \,\langle p',\Lambda'|\,\bar\psi\!\left(-\tfrac{z}{2}\right)\Gamma\,\mathcal W\,\psi\!\left(\tfrac{z}{2}\right)\,|p,\Lambda\rangle\Big|_{z^+=0},$$7. For $$W^{[\Gamma]q}_{\Lambda'\Lambda}(P,x,\boldsymbol{k}_\perp,\Delta,n) =\frac{1}{2}\int\frac{dz^-\,d^2\boldsymbol{z}_\perp}{(2\pi)^3}\, e^{ixP^+z^- - i\boldsymbol{k}_\perp\cdot\boldsymbol{z}_\perp} \,\langle p',\Lambda'|\,\bar\psi\!\left(-\tfrac{z}{2}\right)\Gamma\,\mathcal W\,\psi\!\left(\tfrac{z}{2}\right)\,|p,\Lambda\rangle\Big|_{z^+=0},$$8 quarks the dipole indicates OAM aligned with the proton helicity, whereas for $$W^{[\Gamma]q}_{\Lambda'\Lambda}(P,x,\boldsymbol{k}_\perp,\Delta,n) =\frac{1}{2}\int\frac{dz^-\,d^2\boldsymbol{z}_\perp}{(2\pi)^3}\, e^{ixP^+z^- - i\boldsymbol{k}_\perp\cdot\boldsymbol{z}_\perp} \,\langle p',\Lambda'|\,\bar\psi\!\left(-\tfrac{z}{2}\right)\Gamma\,\mathcal W\,\psi\!\left(\tfrac{z}{2}\right)\,|p,\Lambda\rangle\Big|_{z^+=0},$$9 quarks the distortion tends to favor OAM anti-aligned with the proton helicity. The average transverse momentum field,

$P=(p'+p)/2$0

shows a vortex-like flow, and for $P=(p'+p)/2$1 quarks the sign of the OAM density changes at about $P=(p'+p)/2$2 fm from the transverse center. A similar qualitative picture is obtained in the light-front chiral quark-soliton model (LF$P=(p'+p)/2$3QSM) (Lorce et al., 2012).

The same papers compared three OAM extractions in the valence sector: kinetic OAM from Ji’s sum rule, canonical OAM from the Wigner integral, and the pretzelosity proxy. The reported values are:

Model	Definition	Values
LFCQM	Kinetic (Ji)	$P=(p'+p)/2$4
LFCQM	Canonical (Jaffe–Manohar)	$P=(p'+p)/2$5
LFCQM	TMD-based proxy	$P=(p'+p)/2$6
LF$P=(p'+p)/2$7QSM	Kinetic (Ji)	$P=(p'+p)/2$8
LF$P=(p'+p)/2$9QSM	Canonical (Jaffe–Manohar)	$\Delta=p'-p$0
LF$\Delta=p'-p$1QSM	TMD-based proxy	$\Delta=p'-p$2

In both models, the total quark OAM coincides across definitions, whereas flavor-separated values do not. In the LFCQM both $\Delta=p'-p$3 and $\Delta=p'-p$4 kinetic OAM are positive, while canonical and TMD-based extractions give a negative $\Delta=p'-p$5-quark OAM. In the LF$\Delta=p'-p$6QSM, kinetic OAM gives a small negative $\Delta=p'-p$7 and positive $\Delta=p'-p$8, closer in sign to lattice indications for the isovector combination (Lorce et al., 2012).

Models with explicit gluonic degrees of freedom alter this pattern. In the dressed-quark calculation of Mukherjee, Nair, and Ojha, performed in light-front gauge with the gauge link set to unity, $\Delta=p'-p$9 peaks at the origin and $x=k^+/P^+$0 exhibits a clear dipole structure. In that model both kinetic and canonical OAM are negative, and their magnitudes differ once gluons are included, in contrast to quark-only models (Mukherjee et al., 2014). Related perturbative dressed-quark analyses found a negative longitudinal spin-orbit correlation, with $x=k^+/P^+$1 and consequently $x=k^+/P^+$2 (Mukherjee, 2015).

The model landscape is broader. Light-cone spectator models compute all 16 leading-twist distributions for $x=k^+/P^+$3 and $x=k^+/P^+$4 quarks and emphasize strong directional dependence in the transverse-spin sectors, while soft-wall AdS/QCD-inspired quark-diquark models report circular $x=k^+/P^+$5 and $x=k^+/P^+$6 distributions, dipolar $x=k^+/P^+$7 and $x=k^+/P^+$8, quadrupolar pretzelous structures, and a Soffer-bound-type inequality for Wigner distributions (Liu et al., 2015, Chakrabarti et al., 2017).

6. Multipoles, skewness, and longitudinal extensions

A major later development was the complete multipole decomposition of the transverse phase space. This framework classifies every leading-twist Wigner distribution into basic multipoles in $x=k^+/P^+$9 and $\mathcal W$0, separately in the T-even and T-odd sectors. It provides a transparent interpretation of the spin-spin and spin-orbit correlations encoded in each channel, and it shows how impact-parameter GPD multipoles and TMD multipoles arise as projections of the more general phase-space patterns (Lorcé et al., 2015).

Nonzero skewness introduces an additional layer of structure. In a light-front dressed-quark model, increasing $\mathcal W$1 produces characteristic distortions in the spatial and momentum correlations, including dipole and quadrupole patterns, asymmetries, and localization effects. These features were attributed to spin-orbit correlations and quantum interference between light-front wave-function components with differing orbital angular momentum. The forward-limit symmetry progressively deforms as $\mathcal W$2 increases, and mixed-space ridges sharpen, indicating stronger localization in phase space (Jana et al., 24 Jul 2025).

A complementary line of work introduces a boost-invariant longitudinal coordinate, either $\mathcal W$3 or $\mathcal W$4, as the Fourier conjugate of the skewness. In this formulation, the Fourier transform of $\mathcal W$5-dependent GTMDs defines Wigner distributions in longitudinal position space, which exhibit diffraction patterns analogous to optical single-slit interference. This extends the usual five-dimensional picture to a six-dimensional light-front phase space $\mathcal W$6 and retains the same 16 independent leading-twist polarization structures (Maji et al., 2022, Yang et al., 10 May 2025).

These generalized constructions sharpen two points. First, the standard five-dimensional distributions are not the only possible light-front phase-space representation; they are the $\mathcal W$7 sector of a broader off-forward geometry. Second, while Wigner distributions themselves are not directly measurable, their $\mathcal W$8 dependence and twist-three projections create experimentally relevant observables. In DVCS, the twist-three effective axial Compton form factor enters the $\mathcal W$9 modulation of the longitudinal target-spin asymmetry,

$\boldsymbol{k}_\perp$00

with the coefficient $\boldsymbol{k}_\perp$01 controlled by $\boldsymbol{k}_\perp$02 and therefore sensitive to $\boldsymbol{k}_\perp$03 and OAM-associated GTMD/Wigner content. HERMES reported a sizable $\boldsymbol{k}_\perp$04 modulation, enabling an extraction of $\boldsymbol{k}_\perp$05 within the Wandzura–Wilczek approximation (Courtoy et al., 2013).

Taken together, these developments establish leading-twist quark Wigner distributions as the central phase-space language for nucleon tomography. They unify GPD and TMD information, isolate the GTMD structures that carry spin-orbit content, and provide the most direct leading-twist description of how transverse position and transverse momentum correlations encode quark orbital motion. At the same time, the literature makes clear that gauge-link geometry, parity constraints, and twist-three dynamics are not technical details but defining ingredients of the subject.