From Sub-eikonal DIS to Quark Distributions and their High-Energy Evolution

Abstract: Relating the high-energy dipole description of deep-inelastic scattering to the standard light-ray operator formulation at finite Bjorken $x_B$ is essential for connecting the small-$x$ framework to the usual partonic description. I demonstrate that this connection already emerges at the first sub-eikonal order. At the differential level, the first sub-eikonal correction is governed by a quark TMD-like light-ray operator. In the inclusive limit, after complete phase-space integration, it reconstructs the standard nonlocal quark and helicity distributions at nonzero $x_B$. I then show independently that the same inclusive operator content follows from the high-energy limit of the leading-twist non-local operator product expansion, thereby establishing an explicit operator-level bridge between the shock-wave formalism and the non-local light-cone expansion. I further discuss the high-energy evolution of the corresponding operators at $x_B=0$. Rewriting the evolution equations in terms of dipole-type operator combinations, I identify an operator basis whose bilocal building blocks vanish in the zero-dipole-size limit, making the small-dipole behavior and the leading-logarithmic structure manifest. In the double-logarithmic approximation the evolution equations admit the usual mixed longitudinal-transverse Bessel-type solution when the transverse phase space is treated independently. When the transverse phase space is instead constrained by longitudinal ordering, the second logarithm is converted into a logarithm of energy, and in the symmetric double-logarithmic regime one recovers the fixed-coupling Kirschner-Lipatov exponent with the full finite-$N_c$ color factor $C_F$.

• The paper presents the first sub-eikonal corrections in DIS, establishing an operator bridge between shock-wave methods and light-ray OPE formalisms.
• It derives explicit operator structures for quark and helicity distributions and demonstrates matching with mixed and pure double-logarithmic evolution regimes.
• The work provides a framework for extending small-x resummation to incorporate finite-xB effects, crucial for precision QCD analyses at modern colliders.

Sub-eikonal Corrections in Deep-Inelastic Scattering and the Operator Bridge to Quark Distributions

Overview

This work systematically investigates the formal and practical connections between the high-energy dipole (shock-wave/Wilson line) description of deep-inelastic scattering (DIS) at small Bjorken $x_B$ and the light-ray nonlocal operator product expansion (OPE) underlying standard parton distribution functions at finite $x_B$. By computing first sub-eikonal corrections in the small-$x$ framework, the analysis establishes how TMD-like and helicity quark operators emerge from the Wilson-line formalism and maps their evolution, revealing both mixed and pure double-logarithmic regimes. The resulting formalism explicitly demonstrates how partonic light-cone quark and helicity distributions are reconstructed from the high-energy operator content when going beyond the strict eikonal (leading-power) limit, and clarifies the operator basis and evolution equations controlling dynamics at $x_B=0$.

Operator-Level Connection: From Dipoles to Light-Ray Distributions

In the high-energy (small-$x_B$) limit, DIS is commonly described in terms of Wilson lines and dipole operators, successfully capturing gluon saturation and resumming large high-energy logarithms (BFKL/BK/JIMWLK equations). However, at finite $x_B$, the standard language involves nonlocal quark and helicity operators along the light cone, as in DGLAP or light-ray OPE.

This work calculates the first sub-eikonal corrections to the high-energy dipole limit using the shock-wave formalism with quark propagators partially in and out of the background field. A crucial observation is that these corrections, at the differential level, are governed by light-ray operators similar to quark TMDs and, after full phase-space integration, reconstruct the standard longitudinally nonlocal (light-cone) quark and helicity distributions at nonzero $x_B$. The explicit operator structures are constructed as

$$Q_1(x_\perp, x_B) \sim \int dx^+ dy^+\, e^{ix_B P^-(x^+ - y^+)}\, \bar\psi(y^+,x_\perp)[y^+,x^+] n_1\,\psi(x^+,x_\perp)$$

with analogous expressions for the helicity operator $Q_5$.

An essential technical result is the explicit noncommutativity between the strict small-$x_B$ limit and full phase-space integration. Projecting too early (setting $x_B$0 in the integrand) yields only naive collinear PDFs at $x_B$1, while correct inclusive results at finite $x_B$2 are obtained by completing the integration first, with the high-energy limit eikonalizing the transverse kernel while retaining nontrivial longitudinal $x_B$3-dependence.

Independent Confirmation: Non-local OPE in the High-Energy Limit

The author independently derives the same operator structures by taking the high-energy limit of the leading-twist nonlocal OPE à la Balitsky-Braun, evaluating the time-ordered product of electromagnetic currents. In this limit, the straight gauge link in the light-ray quark operator reduces to the sum of semi-infinite Wilson lines characteristic of the shock-wave approach, and the matrix elements coincide (for the relevant Lorentz projections) with those constructed from sub-eikonal shock-wave diagrams. This cross-verification establishes an explicit operator-level bridge between the small-$x_B$4 dipole-based and finite-$x_B$5 operator-based formalisms, making transparent the connection between the languages of high-energy and standard collinear factorization.

High-Energy Evolution and Emergence of Double Logarithms

The high-energy evolution of the resulting operators at $x_B$6 is systematically developed. The evolution equations for the leading sub-eikonal quark operators $x_B$7 and $x_B$8 involve dipole-type (bilocal) combinations,

$x_B$9

Suppressing subleading local terms, in the double-logarithmic approximation, both polarized and unpolarized quark operators obey the same ladder-type evolution: $x$0 where $x$1 encodes rapidity/energy logarithm, and $x$2 is a transverse logarithm.

The regime of the second logarithm is clarified:

• Independent Integration: When the transverse phase is unconstrained, this yields the canonical mixed double-logarithmic (DLA) Bessel-type solution, with terms $x$3.
• Constrained Phase Space: When the transverse domain is tied to the longitudinal ordering (e.g., via $x$4), the second logarithm becomes a longitudinal energy logarithm, restoring the Kirschner-Lipatov exponent $x$5 (with full finite $x$6 factor $x$7).

This dichotomy demonstrates that the ladder expansion naturally reproduces both the usual DGLAP/BFKL resummation and, with appropriate kinematic constraint, the correct pure double-logarithmic $x$8 behavior. In the helicity (polarized) channel, the structural features persist, although possible differences beyond the strict ladder approximation are anticipated, consistent with the findings in polarized small-$x$9 literature (cf. Bartels-Ermolaev-Ryskin, Kovchegov et al.).

Implications and Future Developments

The explicit operator realization of sub-eikonal effects in DIS paves the way for systematic power-suppressed operator analyses at high energy, essential for precision matching between the low- and moderate-$x_B=0$0 regions. The equivalence between the high-energy operator evolution and the well-known DGLAP/OPE-based evolution clarifies that the transition from the dipole picture to light-ray operators is already realized at first sub-eikonal order.

This formalism provides a natural framework for extending small-$x_B=0$1 resummation—including polarization effects—directly at the operator level, leveraging the shock-wave approach and systematically incorporating subleading and finite-$x_B=0$2 effects. Further, the identification of the operator basis suitable for evolution in terms of dipole vanishing at zero size disentangles the genuinely small-$x_B=0$3 logarithms from local artifacts.

Key future avenues include:

• Derivation and solution of a closed evolution system for sub-eikonal operators, including all operator mixing (both singlet and non-singlet/quark-gluon).
• Inclusion of higher sub-eikonal orders to probe the relation between higher-twist and high-energy expansions.
• Quantitative phenomenology for polarized DIS structure functions in the small-$x_B=0$4 regime, particularly crucial for the Electron-Ion Collider kinematics.

Conclusion

The paper rigorously establishes that the first sub-eikonal quark correction in the high-energy expansion creates an explicit operator bridge between the small-$x_B=0$5 shock-wave formalism and the conventional finite-$x_B=0$6 light-cone OPE framework. The evolution equations for the associated operators naturally interpolate between mixed and pure double-logarithmic resummations, with the Kirschner-Lipatov exponent emerging precisely when the transverse phase space is kinematically constrained by longitudinal ordering. This operator-level connection provides a robust foundation for future high-energy QCD analyses at the sub-eikonal level and for the consistent treatment of the moderate-to-small-$x_B=0$7 transition in both unpolarized and polarized DIS.