---
title: Light-Front Quark-Scalar-Diquark Model
url: https://www.emergentmind.com/topics/light-front-quark-scalar-diquark-model
type: topic
---

# Light-Front Quark-Scalar-Diquark Model

Searching arXiv for recent and foundational papers on the light-front quark–scalar-diquark model.
arXiv search query: "light-front quark scalar diquark model proton TMD GPD GTMD"
The light-front quark-scalar-diquark model is a two-body valence description of the proton in which an active quark is accompanied by a spectator scalar diquark, and the nonperturbative dynamics are encoded in light-front wave functions (LFWFs). In the cited literature, this framework appears both in soft-wall AdS/QCD-inspired form and in spectator formulations with a point-like quark-diquark vertex, and it is used to compute PDFs, TMDs, GPDs, GTMDs, Wigner distributions, gravitational form factors, and higher-twist correlators [1608.07776, 1702.04557, 1810.09099].

## 1. Effective degrees of freedom and model definition

In its minimal form, the model truncates the proton Fock expansion to a single active quark plus a scalar diquark spectator. For the scalar channel, the diquark has spin \(0\) and helicity \(\lambda_S=0\), while the proton retains helicity \(J^z=\pm \tfrac12\) and the quark helicity takes \(\lambda_q=\pm \tfrac12\) [1702.04557, 1810.09099]. The basic state content is therefore
\[
|q^{\lambda_q}(xP^+,p_\perp),\,S((1-x)P^+,-p_\perp)\rangle,
\]
with \(x\) the quark light-cone momentum fraction and \(p_\perp\) its transverse momentum [2307.16411].

A broader spin-flavor decomposition often embeds the scalar channel into a quark-diquark basis containing both scalar and axial-vector spectators. In the SU(4) basis used for nucleon modeling,
\[
|P;\pm\rangle
=
C_S\,|u\,S^0\rangle^\pm
+
C_V\,|u\,A^0\rangle^\pm
+
C_{VV}\,|d\,A^1\rangle^\pm,
\]
with the scalar sector isolated by retaining only \(|u\,S^0\rangle^\pm\) or, in a pure scalar-diquark treatment, by setting the scalar weight \(C_S=1\) [1608.07776, 2307.16411].

The model is also used in a spectator realization in which the proton is treated as a quark of mass \(m\) coupled to a scalar diquark of mass \(M_s\) through a point-like vertex of strength \(g_s\). In that formulation, no explicit Lagrangian is required for practical calculations because the dynamics are carried by the LFWFs themselves, although an effective Yukawa-type interaction \(g_S\,\bar\Psi_N\,\phi_S\,\psi + \mathrm{h.c.}\) is a standard underlying picture [1810.09099, 2310.03592].

## 2. Light-front state expansion and helicity structure

For the scalar-diquark sector, the proton state is expanded as
\[
|P;\pm\rangle_S
=
\sum_{\lambda_q=\pm \frac12}
\int [dx\,d^2p_\perp]\,
\psi^\pm_{\lambda_q}(x,p_\perp)\,
|\lambda_q; xP^+,p_\perp;\lambda_S=0\rangle,
\]
with
\[
[dx\,d^2p_\perp]
\equiv
\frac{dx\,d^2p_\perp}{2(2\pi)^3\sqrt{x(1-x)}}.
\]
In the soft-wall AdS/QCD-based construction, the four nonzero helicity amplitudes are [1702.04557, 1608.07776]
\[
\psi^+_{+}(x,p_\perp)=N_S\,\phi_1(x,p_\perp),\qquad
\psi^+_{-}(x,p_\perp)=-\,N_S\,\frac{p^1+i p^2}{xM}\,\phi_2(x,p_\perp),
\]
\[
\psi^-_{+}(x,p_\perp)=N_S\,\frac{p^1-i p^2}{xM}\,\phi_2(x,p_\perp),\qquad
\psi^-_{-}(x,p_\perp)=N_S\,\phi_1(x,p_\perp).
\]

A different but closely related spectator parameterization writes the scalar-channel spinor components as [1810.09099, 2307.16411]
\[
\psi^+_{+}(x,p_\perp)=\frac{m+xM}{x}\,\varphi(x,p_\perp),\qquad
\psi^+_{-}(x,p_\perp)=-\,\frac{p_x+i p_y}{x}\,\varphi(x,p_\perp),
\]
with \(\psi^-_{\pm}(x,p_\perp)=\bigl[\psi^+_{\mp}(x,p_\perp)\bigr]^*\). In both implementations, the spin dependence is carried by explicit transverse-momentum prefactors and the radial dependence is isolated in one or two scalar profile functions.

Normalization is imposed sector by sector. A typical condition is
\[
\sum_{\lambda_q}\int dx\,d^2p_\perp\,|\psi^+_{\lambda_q}(x,p_\perp)|^2=1,
\]
or its equivalent with the light-front measure written explicitly [1702.04557, 1608.07776]. In the AdS/QCD-based nucleon model, the free two-body invariant mass is
\[
M_0^2(x,k_\perp)=\frac{k_\perp^2+m_q^2}{x}+\frac{k_\perp^2+m_S^2}{1-x},
\]
with \(m_q=0\) in the massless-quark version of the model [1608.07776].

## 3. Radial wave-function ansätze and parameterizations

The best-known implementation uses soft-wall AdS/QCD-inspired profile functions. In momentum space,
\[
\phi_i(x,p_\perp)
=
\frac{4\pi}{\kappa}\,
\sqrt{\frac{\ln(1/x)}{1-x}}\,
x^{a_i}(1-x)^{b_i}
\exp\!\left[
-\frac{\delta\,p_\perp^2}{2\kappa^2}\,
\frac{\ln(1/x)}{(1-x)^2}
\right],
\]
with \(\kappa=0.4\,\mathrm{GeV}\); the strict AdS/QCD limit corresponds to \(a_i=b_i=0\) and \(\delta=1\) [1702.04557, 1608.07776]. The parameters \((a_i,b_i,\delta)\) are fitted either to unpolarized PDFs or to nucleon form-factor data, depending on the analysis [1702.04557, 1608.07776].

A spectator alternative employs a single scalar profile
\[
\varphi(x,p_\perp)
=
-\,\frac{g_s}{\sqrt{1-x}}\;
\frac{x(1-x)}{p_\perp^2+\bigl[xM_s^2+(1-x)m^2-x(1-x)M^2\bigr]},
\]
together with constituent masses and couplings fitted to form factors [1810.09099]. Higher-twist analyses also use an AdS/QCD-motivated Gaussian form,
\[
\phi(x,p_\perp)=
N_s\,x^{a_s}(1-x)^{b_s}
\exp\!\left\{
-\frac{p_\perp^2+m^2(1-x)+M_s^2x-x(1-x)M^2}{2\kappa^2x(1-x)}
\right\},
\]
with \(N_s\) fixed by \(\int dx\,d^2p_\perp\,|\phi|^2=1\) [2307.16411].

| Implementation | Radial input | Representative parameters |
|---|---|---|
| Soft-wall AdS/QCD nucleon model | Two profiles \(\phi_{1,2}\) with Gaussian \(p_\perp\)-dependence and powers \(x^{a_i}(1-x)^{b_i}\) | \(\kappa=0.4\,\mathrm{GeV}\); at \(\mu_0=0.313\,\mathrm{GeV}\): for \(u\), \(a_1=0.280\), \(b_1=0.1716\), \(a_2=0.840\), \(b_2=0.2284\), \(\delta=1.0\), \(N_S=2.0191\) [1608.07776, 2301.09536] |
| TMD-focused AdS/QCD scalar sector | Same functional form, fitted to PDFs at a different initial scale | \(\kappa=0.4\,\mathrm{GeV}\); at \(\mu_0=0.8\,\mathrm{GeV}\): \(a_1^u\simeq0.15\), \(b_1^u\simeq0.70\), \(a_2^u\simeq0.28\), \(b_2^u\simeq0.65\), \(\delta^u\simeq1.00\) [1702.04557] |
| Spectator scalar-diquark model | Single propagator-like profile \(\varphi\) | \(m=0.33\,\mathrm{GeV}\), \(M_s=0.822\pm0.053\,\mathrm{GeV}\), \(g_s=0.847\pm0.111\) [1810.09099] |

A more dynamical realization arises in the light-front Faddeev-Bethe-Salpeter approach, where a contact two-quark interaction produces a scalar-diquark subamplitude \(\mathcal F(M_{12}^2)\), the proton valence wave function is obtained by symmetrizing a bachelor-quark vertex \(\Gamma(x,k_\perp)\), and model parameters are tuned to reproduce the proton mass and Dirac form factor [2211.10959]. This suggests a route from phenomenological quark-diquark LFWFs to a bound-state construction in which the scalar diquark is not merely an external spectator parameter.

## 4. Overlap representations for PDFs and TMDs

The scalar-diquark model is particularly suited to overlap calculations. For collinear PDFs, the unpolarized, helicity, and transversity distributions are obtained by integrating LFWF overlaps over transverse momentum. In the scalar channel these take the form [1801.01718, 1608.07776]
\[
f_1^q(x)=
\int\frac{d^2k_\perp}{16\pi^3}
\left[
|\phi_1|^2+\frac{k_\perp^2}{x^2M_N^2}|\phi_2|^2
\right],
\]
\[
g_1^q(x)=
\int\frac{d^2k_\perp}{16\pi^3}
\left[
|\phi_1|^2-\frac{k_\perp^2}{x^2M_N^2}|\phi_2|^2
\right],
\]
and \(h_1^q(x)\) from the corresponding chiral-odd overlap. In the nucleon model based on soft-wall AdS/QCD LFWFs, the helicity and transversity distributions were reported to agree with phenomenological fits, axial and tensor charges were reported to agree with experimental or phenomenological extractions, and the Soffer bound was verified numerically at all scales studied [1608.07776, 1801.01718].

At leading twist and for T-even TMDs, the scalar sector yields explicit overlap expressions. Representative examples are [1702.04557]
\[
f_1^s(x,p_\perp^2)=\frac{1}{16\pi^3}
\left(
|\psi^+_{+}|^2+|\psi^+_{-}|^2
\right),
\]
\[
g_{1L}^s(x,p_\perp^2)=\frac{1}{16\pi^3}
\left(
|\psi^+_{+}|^2-|\psi^+_{-}|^2
\right),
\]
\[
h_1^s(x,p_\perp^2)=\frac{1}{16\pi^3}
\left[
\psi^+_{+}\psi^{-\,\dagger}_{-}
+
\psi^+_{-}\psi^{-\,\dagger}_{+}
\right].
\]
The remaining leading-twist T-even functions \(g_{1T}\), \(h_{1L}^{\perp}\), \(h_{1T}\), and \(h_{1T}^{\perp}\) follow from off-diagonal helicity overlaps [1702.04557].

A notable numerical result is the approximate \(x\)-\(p_\perp^2\) factorization,
\[
f_1^s(x,p_\perp^2)
\simeq
f_1^s(x)\,
\frac{e^{-p_\perp^2/\langle p_\perp^2\rangle}}{\pi\,\langle p_\perp^2\rangle},
\]
with
\[
\langle p_\perp^2\rangle
=
\frac{\int dx\,d^2p_\perp\,p_\perp^2\,f_1^s(x,p_\perp^2)}
{\int dx\,d^2p_\perp\,f_1^s(x,p_\perp^2)}.
\]
The AdS/QCD scalar-diquark model does not factorize analytically into a product of pure \(x\)- and \(p_\perp\)-functions, but the Gaussian ansatz was found to be realized to very good approximation [1702.04557].

## 5. Generalized distributions, spatial imaging, and orbital structure

The same overlap machinery extends to GPDs, GTMDs, and Wigner distributions. For chiral-odd GPDs, the scalar-diquark model provides overlap formulas for \(H_T\), \(E_T\), \(\widetilde H_T\), and \(\widetilde E_T\) at both zero and nonzero skewness, with active-quark momentum arguments shifted as
\[
x'=\frac{x-\zeta}{1-\zeta},\qquad
x''=\frac{x+\zeta}{1+\zeta},
\]
and similarly for \(k'_\perp\) and \(k''_\perp\) [1509.00598]. Their two-dimensional Fourier transforms define impact-parameter-space GPDs, and a one-dimensional Fourier transform in \(\zeta\) introduces a longitudinal impact coordinate \(\sigma=\tfrac12 b^-P^+\). In the scalar model, \(H_T\), \(E_T\), and \(\widetilde H_T\) show diffraction-like patterns in \(\sigma\)-space, whereas \(\widetilde E_T\) does not display a pronounced diffraction pattern because it is odd in \(\zeta\) [1509.00598].

Phase-space imaging is developed further through Wigner distributions and GTMDs. The twist-2 Wigner operator is
\[
W^{[\Gamma]}(\Delta_\perp,p_\perp,x;S)
=
\frac12
\int\frac{dz^-\,d^2z_\perp}{(2\pi)^3}
e^{i\,p\cdot z}
\langle P'';S|
\overline\psi(-\tfrac z2)\,\Gamma\,\mathcal W\,\psi(\tfrac z2)
|P';S\rangle_{z^+=0},
\]
and its Fourier transform in \(\Delta_\perp\) yields
\[
\rho^{[\Gamma]}(b_\perp,p_\perp,x;S)
=
\int\frac{d^2\Delta_\perp}{(2\pi)^2}
e^{-i\,\Delta_\perp\cdot b_\perp}
W^{[\Gamma]}(\Delta_\perp,p_\perp,x;S).
\]
Combining proton and quark polarizations generates 16 independent twist-2 Wigner distributions, while the GTMD decomposition contains 16 complex functions of
\[
(x,\zeta,k_\perp^2,k_\perp\!\cdot\!\Delta_\perp,\Delta_\perp^2),
\]
of which only 10 survive at \(\zeta=0\) [1810.09099]. In the scalar-diquark sector, \(\rho_{UU}\) and \(\rho_{LL}\) are circularly symmetric in both \(b_\perp\) and \(p_\perp\); \(\rho_{UL}\) and \(\rho_{LU}\) display dipolar patterns and a quadrupole in mixed space; and \(F_{1,4}\) at \(\zeta=0\) is directly related to quark orbital angular momentum [1810.09099].

Momentum-space densities also reveal transverse-shape information. In the scalar model,
\[
\rho_{UU}(p_\perp)\equiv f_1^{(1)}(p_\perp^2),\qquad
\rho_{TL}(p_\perp;S_\perp,\lambda)
\equiv
\frac12 f_1^{(1)}(p_\perp^2)
+
\frac{\lambda\,p_\perp\!\cdot\!S_\perp}{2M}\,
g_{1T}^{(1)}(p_\perp^2),
\]
with \(f^{(1)}(p_\perp^2)\equiv \int dx\,f(x,p_\perp^2)\). The transverse-shape relation
\[
\frac{\rho_{\mathrm{REL}\,T}(p_\perp,n)}{M\,\tilde f_1(p_\perp^2)}
=
1
+
\frac{\tilde h_1(p_\perp^2)}{\tilde f_1(p_\perp^2)}\cos\phi_n
+
\frac{p_\perp^2}{2M^2}
\frac{\tilde h_{1T}^{\perp}(p_\perp^2)}{\tilde f_1(p_\perp^2)}
\cos(2\phi-\phi_n)
\]
shows non-spherical, “pretzel” distortions driven by \(h_{1T}^{\perp}\) [1702.04557].

For energy-momentum structure, the scalar-diquark LFWFs generate gravitational form factors through overlaps with \(T_q^{++}\):
\[
A_q^{(S)}(Q^2)=
\int\frac{dx\,d^2k_\perp}{16\pi^3}\,x\,[\cdots],\qquad
B_q^{(S)}(Q^2)=
-\frac{2M}{q^1-iq^2}
\int\frac{dx\,d^2k_\perp}{16\pi^3}\,x\,[\cdots].
\]
Their two-dimensional Fourier transforms define the longitudinal momentum density \(\rho_q(b_\perp)\) and the angular momentum density \(J_q(b_\perp)\) in the transverse plane [1712.02110].

T-odd extensions require an explicit final-state interaction. In the holographic scalar-diquark model with one-gluon-exchange phase,
\[
\Phi(x,k_\perp)
=
\frac{e_1e_2}{8\pi}\,[k_\perp^2+B(x)]\,g(x,k_\perp^2),
\]
the T-odd GTMDs \(F_{1,2}^o\) and \(H_{1,1}^o\) emerge from the imaginary part of \(\psi^*\psi\), and in the forward limit \(\Delta_\perp\to0\) they reduce to the Sivers and Boer-Mulders functions [1902.07051].

## 6. Evolution, higher-twist sectors, and limitations

The scalar-diquark model is defined at a low hadronic scale and then evolved. For integrated PDFs, ordinary DGLAP evolution is applied, either directly or through a simulated DGLAP procedure in which the profile parameters \(a_i^\nu(\mu)\), \(b_i^\nu(\mu)\), and \(\delta^\nu(\mu)\) acquire scale dependence [1608.07776, 1801.01718]. For TMDs, a Collins-Soper-Sterman-type evolution is written in \(b_\perp\)-space:
\[
\tilde F(x,b;\mu)
=
\tilde F(x,b;\mu_0)\,
\exp\!\left\{
\ln(\mu/\mu_0)\,\widetilde K(b;\mu_0)
+
\int_{\mu_0}^{\mu}\frac{d\mu'}{\mu'}\,
\gamma_F[\mu',\mu^2/\mu'^2]
\right\},
\]
with
\[
\widetilde K(b;\mu)
=
\widetilde K(b_*;\mu_b)
+
\int_{\mu}^{\mu_b}\frac{d\mu'}{\mu'}\,\gamma_K(\mu')
-
g_K(b),
\]
followed by the Fourier transform back to momentum space [1702.04557].

Beyond leading twist, the model remains analytically tractable. For T-even twist-3 TMDs, the scalar sector yields explicit formulas for \(e\), \(f^\perp\), \(g_L^\perp\), \(g_T'\), \(g_T^\perp\), \(h_L\), \(h_T\), and \(h_T^\perp\). Because explicit quark-gluon interaction terms are absent and the gauge link is set to unity for T-even observables, Wandzura-Wilczek-type relations become exact in the model, including
\[
e=\frac{m}{M}f_1,\qquad
f^\perp=f_1
\]
for the corresponding twist-3 functions [2302.07165].

At twist 4, the scalar-diquark sector gives closed expressions for \(f_3\), \(g_{3L}\), \(h_{3L}^{\perp}\), \(g_{3T}\), \(h_{3T}\), and \(h_{3T}^{\perp}\), with model relations tying them to leading-twist TMDs. The simplest is
\[
x^2\,f_3^\nu(x,p_\perp^2)
=
\frac{p_\perp^2+m^2}{M^2}\,
f_1^\nu(x,p_\perp^2),
\]
and analogous formulas relate \(g_{3L}\), \(h_{3L}^{\perp}\), \(g_{3T}\), \(h_{3T}\), and \(h_{3T}^{\perp}\) to \(g_{1L}\), \(h_{1L}^\perp\), \(g_{1T}\), \(h_1\), \(h_{1T}\), and \(h_{1T}^\perp\) [2301.09536]. For the twist-4 T-even function \(h_3\), the scalar channel gives
\[
x^2 h_3^\nu(x,p_\perp^2)
=
\frac{1}{16\pi^3}C_S^2N_s^2
\left[
\frac{m}{M}|\varphi_1^\nu|^2
+
\frac{p_\perp^2}{M^2x}|\varphi_2^\nu|^2
\right]^2,
\]
with reported numerical averages
\[
\langle p_\perp\rangle_u\simeq0.26\,\mathrm{GeV},\quad
\langle p_\perp^2\rangle_u\simeq0.07\,\mathrm{GeV}^2,\quad
\langle p_\perp\rangle_d\simeq0.27\,\mathrm{GeV},\quad
\langle p_\perp^2\rangle_d\simeq0.08\,\mathrm{GeV}^2
\]
[2307.16411]. Twist-4 GTMDs have also been parameterized in the scalar sector, and their TMD limit reproduces the previously derived twist-4 T-even TMDs [2310.03592].

The main limitations are tied to the same simplifications that make the model analytic. Many calculations are valence-only and retain only the minimal quark-diquark Fock sector [1608.07776, 2211.10959]. In the chiral-odd GPD construction, the analysis is restricted to the DGLAP region \(x>\zeta\), and all valence GPDs vanish at \(x=\zeta\) because no explicit ERBL contribution is included [1509.00598]. T-odd observables are absent unless an explicit final-state interaction phase is introduced [1902.07051]. In the light-front Faddeev-Bethe-Salpeter approach, the evolved non-polarized valence PDFs suggest that higher Fock components beyond the valence one are needed for comparison with NNPDF4.0 [2211.10959]. A plausible implication is that the light-front quark-scalar-diquark model is most reliable as a controlled low-scale valence framework and as a laboratory for analytic relations among light-front distributions, rather than as a complete description of proton structure at all \(x\) and all twists.

Source: https://www.emergentmind.com/topics/light-front-quark-scalar-diquark-model