Diquark Spectator Model Overview

• Diquark spectator model is an effective quark–diquark framework that approximates a hadron as an active quark and a spectator diquark, available in scalar and axial-vector variants.
• It employs phenomenological vertex regulators and form factors to control ultraviolet sensitivity and finite-momentum effects in calculations of TMDs, GPDs, and fragmentation observables.
• The model provides detailed insights into angular momentum decompositions, rescattering-induced phases, and quasi-distribution analyses relevant for lattice QCD studies.

Searching arXiv for recent and foundational papers on the diquark spectator model and closely related applications. The diquark spectator model is an effective quark–diquark description in which a hadron is approximated as an active quark plus a spectator diquark remnant. In the literature, this framework appears both in a minimal scalar diquark realization and in broader spectator constructions that include axial-vector or vector diquarks, light-cone wave functions, and phenomenological vertex regulators. It has been used to study standard and quasi parton distributions, transverse-momentum-dependent distributions, generalized parton distributions, generalized and gravitational form factors, fragmentation functions, and angular-momentum observables for the nucleon, hyperons, and spin-$\tfrac32$ baryons (Bhattacharya et al., 2018, Gamberg et al., 2014, Liu et al., 2014).

1. Core construction and spectator approximation

In its simplest field-theoretic form, the scalar diquark spectator model is a relativistic Yukawa-type theory containing a spin-$\tfrac12$ hadron field $\Psi$, a quark field $\psi$, and a scalar diquark field $\varphi$,

${\cal L}_{\rm SDM} = \bar{\Psi} \big(i \slashed{\partial}-M\big)\Psi + \bar{\psi} \big(i \slashed{\partial}-m_q\big)\psi + \frac{1}{2}\big(\partial_\mu\varphi\,\partial^\mu\varphi-m_s^2\varphi^2\big) + g\big(\bar{\Psi}\psi\varphi+\bar{\psi}\Psi\varphi\big),$

with the stability condition $M<m_s+m_q$ (Bhattacharya et al., 2018). In this form, the model deliberately retains only the lowest nontrivial ${\cal O}(g^2)$ contribution, neglecting virtual diagrams that contribute only at $x=1$ (Bhattacharya et al., 2018).

The defining spectator approximation consists of inserting a complete set of intermediate states into the quark correlator and truncating it to a single on-shell diquark state. In nucleon applications this state is taken to be either a scalar diquark $X=s$ or an axial-vector diquark $\tfrac12$0, so that the hadron effectively fluctuates as $\tfrac12$1 (Gamberg et al., 2014, Gamberg et al., 2015). In spectator language, the active quark is the line coupled to the probe, while the diquark carries the remainder of the hadron momentum and quantum numbers.

A closely related Abelian implementation is the QED-like scalar diquark model used for angular-momentum decomposition. There the proton is modeled as a quark plus scalar diquark with a point-like quark–nucleon–diquark vertex of coupling $\tfrac12$2, and photon couplings to the quark and spectator denoted $\tfrac12$3 and $\tfrac12$4. A key assumption is that the target is electromagnetically neutral, emulating color neutrality, even though $\tfrac12$5 and $\tfrac12$6 are kept as independent parameters in the displayed formulas (Amor-Quiroz et al., 2019).

2. Model variants, vertices, and spin structure

The diquark spectator model is not a single unique ansatz but a family of related constructions that differ mainly in the spectator spin assignment, vertex structure, and ultraviolet regulator. A common covariant choice uses a scalar vertex $\tfrac12$7 and an axial-vector vertex

$\tfrac12$8

with a dipole form factor

$\tfrac12$9

for $\Psi$0 (Tan et al., 2022). This regulator suppresses large virtualities and renders transverse-momentum integrals finite.

Other applications use different nonperturbative dressings. In $\Psi$1-hyperon fragmentation, the hyperon–quark–diquark vertex is taken as

$\Psi$2

with $\Psi$3, and the ultraviolet region is controlled by a Gaussian form factor

$\Psi$4

(Yang et al., 2017, Wang et al., 2018). In proton twist-3 GPD calculations, an exponential form factor is used at the nucleon–quark–diquark vertex specifically to regularize the transverse-momentum integration and to drive all distributions to zero as $\Psi$5 and $\Psi$6 (Tan et al., 2024). In a light-cone spectator-diquark model for nucleon generalized form factors, the momentum dependence is encoded through the Brodsky–Huang–Lepage ansatz

$\Psi$7

for scalar and axial-vector spectators (Liu et al., 2014).

The treatment of the spin-1 spectator also varies. Some calculations include both transverse and longitudinal polarization states of the axial-vector diquark (Kaur et al., 29 Jan 2026), whereas the proton twist-3 GPD study keeps only the physical light-cone transverse polarization states through a specific polarization tensor (Tan et al., 2024). In still other quasi-TMD applications the simplification $\Psi$8 is adopted for the axial-vector propagator, with the explicit caveat that unphysical polarizations are then included (Tan et al., 2022). In light-cone spectator-diquark models, relativistic spin effects are incorporated by Melosh–Wigner rotation for both the quark and the spin-1 diquark, so that instant-form spin states are converted into light-cone helicity amplitudes before overlap calculations are performed (Liu et al., 2014).

3. TMDs, rescattering, and fragmentation observables

A major use of spectator-diquark models is the analytic calculation of transverse-momentum-dependent distributions. In a scalar-plus-axial spectator framework with dipolar vertices, the model yields explicit expressions for the T-even leading-twist TMDs $\Psi$9, $\psi$0, $\psi$1, $\psi$2, $\psi$3, and $\psi$4, and generates the Sivers and Boer–Mulders functions by expanding the gauge link to first order, i.e. through one-gluon exchange between the struck quark and the spectator (Bacchetta et al., 2010). This same framework was then used to construct weighted SIDIS and Drell–Yan asymmetries, with the transverse-momentum weights collapsing TMD convolutions into products of moments (Bacchetta et al., 2010).

In a spin-1 spectator realization of SIDIS, the proton is treated as a two-body light-front state consisting of a spin-$\psi$5 quark and a spin-1 diquark. There the single-spin asymmetry arises from the interference of a real tree amplitude with a one-loop amplitude containing a final-state gluon exchange; the required imaginary phase follows from an on-shell pole,

$\psi$6

and the resulting model gives explicit nonzero Sivers and Boer–Mulders functions through helicity-dependent light-front wave-function overlaps (Kumar et al., 2015). In that calculation the two functions are proportional, differing by a factor of $\psi$7 in the displayed formulas, and both exhibit a node controlled by the factor $\psi$8, with sign reversal near $\psi$9 for the chosen parameters (Kumar et al., 2015).

The same spectator logic has been extended to twist-3. In a scalar diquark model for the quark–gluon–quark correlator, eight time-reversal-odd interaction-dependent twist-3 TMDs are extracted from the Dirac decomposition of $\varphi$0. A dipole form factor,

$\varphi$1

is essential: with a point-like coupling, only $\varphi$2 and $\varphi$3 remain finite, while the other six distributions diverge (Lu et al., 2012). The same model reproduces the ETQS–Sivers identity

$\varphi$4

but does not satisfy the exact QCD equations-of-motion relations between interaction-dependent and interaction-independent twist-3 TMDs, which is explicitly attributed to the lowest-order truncation (Lu et al., 2012).

Fragmentation observables are treated analogously by replacing the initial baryon with a final hadron and a spectator remnant. For $\varphi$5 production, the spectator model with scalar and axial-vector diquarks fits $\varphi$6 to the DSV parametrization at $\varphi$7, using

$\varphi$8

and finds $\varphi$9 to be strange-quark dominated, while the naive T-odd ${\cal L}_{\rm SDM} = \bar{\Psi} \big(i \slashed{\partial}-M\big)\Psi + \bar{\psi} \big(i \slashed{\partial}-m_q\big)\psi + \frac{1}{2}\big(\partial_\mu\varphi\,\partial^\mu\varphi-m_s^2\varphi^2\big) + g\big(\bar{\Psi}\psi\varphi+\bar{\psi}\Psi\varphi\big),$0 is driven mainly by ${\cal L}_{\rm SDM} = \bar{\Psi} \big(i \slashed{\partial}-M\big)\Psi + \bar{\psi} \big(i \slashed{\partial}-m_q\big)\psi + \frac{1}{2}\big(\partial_\mu\varphi\,\partial^\mu\varphi-m_s^2\varphi^2\big) + g\big(\bar{\Psi}\psi\varphi+\bar{\psi}\Psi\varphi\big),$1 quarks and yields a negative transverse polarization ${\cal L}_{\rm SDM} = \bar{\Psi} \big(i \slashed{\partial}-M\big)\Psi + \bar{\psi} \big(i \slashed{\partial}-m_q\big)\psi + \frac{1}{2}\big(\partial_\mu\varphi\,\partial^\mu\varphi-m_s^2\varphi^2\big) + g\big(\bar{\Psi}\psi\varphi+\bar{\psi}\Psi\varphi\big),$2 in both SIDIS and single-inclusive ${\cal L}_{\rm SDM} = \bar{\Psi} \big(i \slashed{\partial}-M\big)\Psi + \bar{\psi} \big(i \slashed{\partial}-m_q\big)\psi + \frac{1}{2}\big(\partial_\mu\varphi\,\partial^\mu\varphi-m_s^2\varphi^2\big) + g\big(\bar{\Psi}\psi\varphi+\bar{\psi}\Psi\varphi\big),$3 annihilation (Yang et al., 2017). The ${\cal L}_{\rm SDM} = \bar{\Psi} \big(i \slashed{\partial}-M\big)\Psi + \bar{\psi} \big(i \slashed{\partial}-m_q\big)\psi + \frac{1}{2}\big(\partial_\mu\varphi\,\partial^\mu\varphi-m_s^2\varphi^2\big) + g\big(\bar{\Psi}\psi\varphi+\bar{\psi}\Psi\varphi\big),$4 Collins function ${\cal L}_{\rm SDM} = \bar{\Psi} \big(i \slashed{\partial}-M\big)\Psi + \bar{\psi} \big(i \slashed{\partial}-m_q\big)\psi + \frac{1}{2}\big(\partial_\mu\varphi\,\partial^\mu\varphi-m_s^2\varphi^2\big) + g\big(\bar{\Psi}\psi\varphi+\bar{\psi}\Psi\varphi\big),$5 has also been computed with scalar and vector spectators and the same Gaussian vertex structure; in that case one-gluon-loop rescattering and eikonal-line contributions generate the T-odd phase, QCD evolution is implemented for the first transverse moment, and the weighted ${\cal L}_{\rm SDM} = \bar{\Psi} \big(i \slashed{\partial}-M\big)\Psi + \bar{\psi} \big(i \slashed{\partial}-m_q\big)\psi + \frac{1}{2}\big(\partial_\mu\varphi\,\partial^\mu\varphi-m_s^2\varphi^2\big) + g\big(\bar{\Psi}\psi\varphi+\bar{\psi}\Psi\varphi\big),$6 asymmetry in ${\cal L}_{\rm SDM} = \bar{\Psi} \big(i \slashed{\partial}-M\big)\Psi + \bar{\psi} \big(i \slashed{\partial}-m_q\big)\psi + \frac{1}{2}\big(\partial_\mu\varphi\,\partial^\mu\varphi-m_s^2\varphi^2\big) + g\big(\bar{\Psi}\psi\varphi+\bar{\psi}\Psi\varphi\big),$7 is found to be positive, sizable, and sensitive to evolution at ${\cal L}_{\rm SDM} = \bar{\Psi} \big(i \slashed{\partial}-M\big)\Psi + \bar{\psi} \big(i \slashed{\partial}-m_q\big)\psi + \frac{1}{2}\big(\partial_\mu\varphi\,\partial^\mu\varphi-m_s^2\varphi^2\big) + g\big(\bar{\Psi}\psi\varphi+\bar{\psi}\Psi\varphi\big),$8 (Wang et al., 2018).

4. GPDs, generalized form factors, and baryon structure

Spectator-diquark models provide a compact framework for off-forward observables. In a light-cone spectator-diquark model with scalar and axial-vector spectators, plus Melosh–Wigner rotation for both the quark and the axial-vector diquark, the twist-two generalized form factors at ${\cal L}_{\rm SDM} = \bar{\Psi} \big(i \slashed{\partial}-M\big)\Psi + \bar{\psi} \big(i \slashed{\partial}-m_q\big)\psi + \frac{1}{2}\big(\partial_\mu\varphi\,\partial^\mu\varphi-m_s^2\varphi^2\big) + g\big(\bar{\Psi}\psi\varphi+\bar{\psi}\Psi\varphi\big),$9 are computed from overlaps of light-cone wave functions. The resulting $M<m_s+m_q$0- and $M<m_s+m_q$1-quark contributions to $M<m_s+m_q$2, $M<m_s+m_q$3, $M<m_s+m_q$4, $M<m_s+m_q$5, $M<m_s+m_q$6, and $M<m_s+m_q$7 are reported to be comparable with lattice QCD results (Liu et al., 2014).

At twist $M<m_s+m_q$8, a proton spectator model with scalar and axial-vector diquarks and an exponential vertex form factor has been used to calculate the chiral-even GPDs at nonzero skewness. The correlator is evaluated by contour integration in $M<m_s+m_q$9, with the contributing pole depending on whether ${\cal O}(g^2)$0 lies in the ERBL or DGLAP region. A central result is that almost all twist-3 GPDs are discontinuous at ${\cal O}(g^2)$1, with ${\cal O}(g^2)$2 and ${\cal O}(g^2)$3 as the explicit exceptions; the forward limit yields the twist-3 PDF ${\cal O}(g^2)$4, and the same calculation compares kinetic orbital angular momentum and spin-orbit correlations as extracted from twist-2 and twist-3 GPDs (Tan et al., 2024).

The framework has also been pushed beyond the nucleon. For strange and non-strange baryons, a quark–diquark spectator model with scalar and axial-vector channels computes gravitational form factors from the second Mellin moments of chiral-odd GPDs,

${\cal O}(g^2)$5

and related relations for ${\cal O}(g^2)$6 and ${\cal O}(g^2)$7 (Kaur et al., 29 Jan 2026). In the reported ${\cal O}(g^2)$8-quark sector, ${\cal O}(g^2)$9 for the proton and $x=1$0 for $x=1$1, while $x=1$2 for the proton; this is presented as a complementary sign pattern tied to flavor dependence in transversity-related GPD moments (Kaur et al., 29 Jan 2026).

For the spin-$x=1$3 $x=1$4, the spectator approach generates the full leading-twist set of eight unpolarized and eight polarized GPDs and verifies the required time-reversal symmetry under $x=1$5 (Fu et al., 2023). In the forward limit, the structure functions are obtained from $x=1$6, $x=1$7, $x=1$8, and $x=1$9, and after evolution from $X=s$0 to $X=s$1 the model predicts, for example, a decrease of the struck $X=s$2-quark longitudinal momentum fraction in $X=s$3 from about $X=s$4 to about $X=s$5 (Fu et al., 2023). A later extension to transversity GPDs introduces $X=s$6 independent spin-$X=s$7 tensor GPDs; with the simplifying choice $X=s$8, only $X=s$9, $\tfrac12$00, $\tfrac12$01, and $\tfrac12$02 survive in practice, and the model gives the $\tfrac12$03-quark tensor charge in $\tfrac12$04 as

$\tfrac12$05

with $\tfrac12$06 quoted for the isobar (Fu et al., 21 Aug 2025).

5. Quasi-distributions and finite-momentum extrapolation

A distinct branch of spectator-model work concerns quasi-distributions, i.e. equal-time spatial correlators intended as lattice-QCD proxies for light-cone observables. For standard collinear PDFs one replaces the light-cone bilocal operators by spatial correlators along $\tfrac12$07,

$\tfrac12$08

with analogous definitions for $\tfrac12$09 and $\tfrac12$10 (Gamberg et al., 2015). In the spectator model these quasi-PDFs have support on the entire real axis, reduce analytically to the standard distributions as $\tfrac12$11, and provide a controlled estimate of the finite-boost systematics relevant for lattice QCD (Gamberg et al., 2014, Gamberg et al., 2015).

The main numerical conclusion is consistent across the two quasi-PDF studies: at intermediate momentum fractions, roughly $\tfrac12$12, quasi-PDFs are within about $\tfrac12$13 of the standard PDFs once $\tfrac12$14, whereas near $\tfrac12$15 much larger boosts, $\tfrac12$16, are required for satisfactory agreement (Gamberg et al., 2014, Gamberg et al., 2015). The breakdown at large $\tfrac12$17 is traced to the spectator on-shell kinematics, specifically to factors of $\tfrac12$18, and the same studies emphasize that the Soffer positivity bound should not be assumed for quasi-PDFs at finite $\tfrac12$19 (Gamberg et al., 2015).

Generalized quasi-distributions behave similarly. In the scalar diquark model, quasi-GPDs are defined from spatial correlators with either $\tfrac12$20 or $\tfrac12$21, remain continuous functions of $\tfrac12$22, and reduce analytically to the standard GPDs as $\tfrac12$23 (Bhattacharya et al., 2018). Finite-momentum mismatch is largest near $\tfrac12$24, $\tfrac12$25, and in narrow ERBL regions, and the moment analysis of the extended twist-2 study shows that the $\tfrac12$26 definitions are preferred for second moments related to quark total angular momentum, because the $\tfrac12$27 projection receives an additional higher-twist contamination from the $\tfrac12$28 form factor (Bhattacharya et al., 2019).

The quasi-program has also been extended to T-odd TMDs. In a spectator model with scalar and axial-vector diquarks, the quasi Sivers and quasi Boer–Mulders functions are defined by replacing the light-cone projectors $\tfrac12$29 and $\tfrac12$30 with $\tfrac12$31 and $\tfrac12$32. The resulting quasi-distributions reduce analytically to the standard TMDs when $\tfrac12$33, and their first transverse moments are reported to be fair approximations to the standard ones, again at the $\tfrac12$34 level, in the region $\tfrac12$35 when $\tfrac12$36 (Tan et al., 2022).

6. Angular momentum, spectator torque, and known limitations

One of the most detailed conceptual uses of the scalar diquark spectator model concerns the distinction between Ji and Jaffe–Manohar angular-momentum decompositions. The Ji relation for total angular momentum,

$\tfrac12$37

is contrasted with the canonical Jaffe–Manohar definition through the “potential” term

$\tfrac12$38

(Amor-Quiroz et al., 2019). In the scalar diquark model this difference vanishes at one loop,

$\tfrac12$39

but appears at two loops at order $\tfrac12$40, supporting the interpretation that the struck quark acquires a spectator-induced torque through the gauge-link structure and the same rescattering mechanism that produces the Sivers effect (Amor-Quiroz et al., 2019).

The spectator framework has also been used to test factorization itself. In a scalar diquark model for inclusive DIS at sub-asymptotic $\tfrac12$41, the hadronic tensor can be computed exactly and compared with a collinear-factorized approximation built with an inclusive jet function. This analysis argues that a mass-corrected scaling variable,

$\tfrac12$42

extends the region of validity of leading-twist collinear factorization beyond what is obtained with Bjorken $\tfrac12$43, but also shows an intrinsic breakdown near threshold due to exact four-momentum conservation, with failure beginning roughly when $\tfrac12$44 (Guerrero et al., 2020). The same study emphasizes that the decisive limitation is kinematic: in the exact model, $\tfrac12$45 is bounded by phase space, whereas in the factorized approximation the transverse integration extends to infinity (Guerrero et al., 2020).

Across applications, several limitations recur. These models are low-scale effective descriptions, not full nonperturbative QCD calculations; they commonly require dipole, Gaussian, exponential, or BHL-type regulators to suppress ultraviolet sensitivity; and lowest-order truncations can violate exact field-theoretic identities. Concrete examples include the failure of exact QCD equations-of-motion relations at twist $\tfrac12$46 (Lu et al., 2012), approximate or explicit violations of positivity bounds at large $\tfrac12$47 for T-odd $\tfrac12$48-fragmentation functions (Yang et al., 2017, Wang et al., 2018), and the non-guaranteed status of Soffer positivity for quasi-PDFs at finite $\tfrac12$49 (Gamberg et al., 2015). These features do not nullify the spectator approach; rather, they delimit its role as a controlled phenomenological laboratory in which quark–spectator dynamics, rescattering phases, and finite-momentum effects can be computed explicitly.