Threshold resummation of Semi-Inclusive Deep-Inelastic Scattering

This paper derives threshold resummation expressions for semi-inclusive deep-inelastic scattering (SIDIS) by leveraging crossing symmetry with the Drell-Yan process, explicitly determining resummation coefficients up to next-to-next-to-leading log accuracy for both double-soft and single-soft limits while confirming consistency with recent fixed-order and next-to-leading power results.

The Gist

Imagine you are trying to predict the outcome of a very high-speed collision between two tiny particles, like smashing two grains of sand together at nearly the speed of light. In the world of particle physics, this is called Semi-Inclusive Deep-Inelastic Scattering (SIDIS).

Scientists use complex math to predict what happens after the crash. Usually, they calculate the "average" outcome. But sometimes, things get weird. When the collision happens right at the very edge of possibility—like a car braking just inches from a wall—the math breaks down. The numbers get huge and messy because of "soft radiation" (tiny, low-energy particles flying off) and "collinear radiation" (particles flying off in a straight line with the original ones).

This paper is like a new, super-precise GPS for these edge-case collisions. Here is how the authors fixed the problem, explained simply:

1. The Two "Knobs" of the Collision

In this experiment, there are two main dials or "knobs" that control the outcome:

• Knob A (The Incoming Part): How much of the original particle's energy was used? (Physicists call this x).
• Knob B (The Outgoing Part): How much of the energy ended up in the new particle we are watching? (Physicists call this z).

Usually, scientists study what happens when both knobs are turned all the way to the limit (the "Double Soft" limit). It's like a car coming to a complete stop.

2. The New Discovery: The "One-Knob" Limit

The authors realized that in the real world (and especially at the future Electron-Ion Collider), we often care about situations where only one knob is at the limit, while the other is somewhere in the middle.

• Scenario 1: The incoming energy is maxed out, but the outgoing particle is just "okay."
• Scenario 2: The outgoing particle is maxed out, but the incoming energy is just "okay."

The paper says: "Hey, we can predict these 'one-knob' scenarios just as well as the 'two-knob' ones!"

3. The Magic Trick: Crossing the Street

The authors didn't start from scratch. They used a clever trick called Crossing Symmetry.

• Imagine you have a map of a city (Drell-Yan process, a different type of particle collision).
• Now, imagine you flip the map over and look at it from the other side.
• The paper shows that the rules for the "one-knob" limit in our SIDIS experiment are mathematically identical to the rules for a specific limit in that other experiment.
• It's like realizing that the recipe for baking a cake is the same as the recipe for making a pie, just with the ingredients swapped around. They took the known "recipe" for the other experiment and adapted it for SIDIS.

4. The "Resummation" (Fixing the Messy Math)

When you get close to the limit, the math produces infinite loops of errors (like a feedback loop in a microphone).

• The Problem: Standard calculations say, "Add up the small errors." But near the limit, the errors are so big that adding them up one by one never works.
• The Solution: The authors developed a method called Resummation. Instead of adding errors one by one, they grouped them all together into a single, smooth curve.
• The Analogy: Imagine trying to count every single grain of sand on a beach. It's impossible. Instead, you measure the volume of the beach and calculate the total sand. That's what they did with the math—they stopped counting individual tiny particles and calculated the "flow" of the radiation instead.

5. Why This Matters

This isn't just abstract math. The Electron-Ion Collider (EIC) is being built to smash particles together to understand the "glue" (strong force) that holds the universe together.

• To get accurate data from the EIC, scientists need to know exactly what happens at these extreme limits.
• This paper provides the "instruction manual" for those extreme limits.
• It confirms that the theory works perfectly, even when you only push one of the dials to the limit.

Summary

Think of this paper as a specialized weather forecast.

• Standard forecasts tell you what happens on a normal day.
• This paper tells you exactly what happens during a hurricane (the "threshold" limit).
• Even better, it tells you what happens when the wind is blowing hard from only one direction (the "single soft" limit), which previous forecasts couldn't predict accurately.

By doing this, the authors have given future physicists the tools to interpret the most extreme and exciting data coming from the next generation of particle accelerators.

Technical Summary

1. Problem Statement

Semi-Inclusive Deep-Inelastic Scattering (SIDIS) is a critical process for understanding the internal structure of hadrons and is a primary physics goal for the upcoming Electron-Ion Collider (EIC). While fixed-order calculations for SIDIS have recently been completed up to Next-to-Next-to-Leading Order (NNLO) [3–5], these perturbative expansions suffer from large logarithmic corrections near the kinematic threshold.

In SIDIS, the cross-section depends on two scaling variables:

• Bjorken variable ($x$): The fraction of the incoming hadron's momentum carried by the struck parton.
• Fragmentation variable ($z$): The fraction of the final-state parton's momentum carried by the detected hadron.

Near the threshold ($x \to 1$ and/or $z \to 1$), large logarithms of the form $\ln(1-x)$ and $\ln(1-z)$ appear at every order in the strong coupling constant $\alpha_s$, potentially spoiling the convergence of the perturbative series. These logarithms must be resummed to all orders to achieve precision predictions.

While threshold resummation for the Drell-Yan (DY) process (related to SIDIS by crossing symmetry) is well-established, including cases where only one scaling variable approaches the threshold (single-soft limits), a systematic derivation of SIDIS resummation in these asymmetric single-soft limits was lacking. Specifically, the behavior when $x \to 1$ at fixed $z$, or $z \to 1$ at fixed $x$, required a dedicated formalism.

2. Methodology

The authors employ a direct QCD approach based on renormalization-group (RG) arguments, building upon their previous work on Drell-Yan rapidity distributions [16]. The methodology involves:

• Crossing Symmetry: Utilizing the relation between SIDIS and Drell-Yan. The SIDIS process $H(P_1) + G(q) \to h(P_2) + X$ is the crossed version of the DY process $Q_1 + Q_2 \to G + X$.
• Kinematic Analysis of Soft Limits: The authors rigorously define three distinct soft limits:
  1. Double Soft Limit: Both $x \to 1$ and $z \to 1$. This corresponds to the total available energy approaching the threshold, suppressing all radiation.
  2. Single Soft Limit ($x$-soft): $x \to 1$ with $z$ fixed. Physically, this implies the recoil system is collinear to the incoming parton.
  3. Single Soft Limit ($z$-soft): $z \to 1$ with $x$ fixed. Physically, this implies the recoil system is collinear to the outgoing (fragmenting) parton.
• Phase Space Factorization: By analyzing the phase space measure in $d = 4-2\epsilon$ dimensions, they demonstrate that in all soft limits, the dependence on the scaling variables factorizes into a single "soft scale" $\Lambda^2$.
  • Double soft: $\Lambda^2 \propto Q^2(1-x)(1-z)$.
  • Single soft: $\Lambda^2 \propto Q^2(1-x)$ or $Q^2(1-z)$.
• Mellin Space Resummation: The calculation is performed in Mellin space (conjugate to $x$ and $z$ via variables $N$ and $M$). In this space, the convolution of parton distributions, coefficient functions, and fragmentation functions becomes a simple product. The large logarithms transform into powers of $\ln N$ and $\ln M$.
• Matching to Fixed Order: The resummed expressions are expanded to $\mathcal{O}(\alpha_s^2)$ and matched against the known NNLO fixed-order results [3, 5] to explicitly determine the unknown resummation coefficients.

3. Key Contributions

The paper makes several specific theoretical advancements:

1. Derivation of Asymmetric Resummation: The authors derive the general resummation formalism for SIDIS in the single-soft limits. They show that in these limits, the radiation is purely collinear to either the incoming or the outgoing parton, distinct from the double-soft limit where radiation is soft.
2. Explicit Coefficient Determination: They calculate the resummation coefficients explicitly for the nonsinglet quark channel up to Next-to-Next-to-Leading Logarithmic (NNLL) accuracy.
   • They determine the functions $A_q$ (cusp anomalous dimension), $D_{qq}$ (soft non-collinear function), and $C_{c,qq}$ (hard coefficient).
   • Crucially, they derive the dependence of the $D$ and $C$ coefficients on the non-soft variable (e.g., $D(x\text{-soft})$ depends on $M$).
3. Consistency Checks:
   • They verify that the double-soft limit of their single-soft results reproduces the known inclusive resummation results.
   • They confirm that their results reproduce recent Next-to-Leading Power (NLP) results found in the literature [12].
   • They demonstrate that the single-soft resummation exponents coincide at NLL (due to the equality of spacelike and timelike splitting functions at LO) but differ at NNLL due to NLO corrections.
4. Formalism for Evolution Eigenstates: The paper addresses the technical requirement that resummation in the single-soft limit must be performed in terms of evolution eigenstates (which mix at higher orders), providing a framework for future singlet-channel extensions.

4. Key Results

The resummed coefficient function $C(N, M)$ takes the form:
$$C(N, M) = g_0(\alpha_s) \exp\left[ \frac{1}{\alpha_s} g_1(\lambda) + g_2(\lambda) + \alpha_s g_3(\lambda) + \dots \right]$$
where $\lambda = \beta_0 \alpha_s L$ and $L$ is the large logarithm ($\ln \bar{N}$ or $\ln \bar{M}$).

• Double Soft Limit: The result matches the standard inclusive resummation where $D^{(1)}_{qq} = 0$ and $D^{(2)}_{qq}$ is a constant.
• Single Soft Limits ($x$-soft and $z$-soft):
  • The coefficient $D^{(1)}$ is no longer zero but equals the subtracted leading-order anomalous dimension $\tilde{\gamma}^{(0)}(M)$ (or $\tilde{\gamma}^{(0)}(N)$). This reflects the fact that collinear radiation from the non-soft leg contributes at NLL.
  • The coefficient $D^{(2)}$ includes contributions from the NLO anomalous dimensions and specific functions $F(M)$ or $H(N)$ that account for the difference between the soft scale and the hard scale.
  • The constant term $g_0$ acquires a dependence on the non-soft variable, ensuring consistency with the fixed-order expansion.
• NLP Connection: The authors show that the expansion of their single-soft coefficients in powers of $1/N$ or $1/M$ naturally reproduces the Next-to-Leading Power corrections previously derived for the double-soft limit.

5. Significance

• Precision Physics for the EIC: This work provides the necessary theoretical tools to make high-precision predictions for SIDIS at the future Electron-Ion Collider, particularly in kinematic regions where threshold effects are dominant (large $x$ and/or large $z$).
• Validation of Crossing Symmetry: It provides a non-trivial, explicit verification that the crossing symmetry between SIDIS and Drell-Yan holds not just for fixed-order calculations but also for the complex structure of threshold resummation in asymmetric limits.
• Foundation for Full Phenomenology: While currently restricted to the nonsinglet channel, this paper establishes the "proof of concept" for asymmetric resummation. It paves the way for including singlet contributions (gluons and sea quarks), which are essential for describing data at small $x$ or specific $z$ ranges.
• Approximate Higher Orders: The derived resummation formulas can be used to construct approximate NNLO and N$^3$LO corrections for SIDIS, improving the accuracy of global QCD fits and parton distribution function (PDF) extractions.

In summary, this paper successfully extends the state-of-the-art in QCD resummation to the semi-inclusive case, resolving the behavior of large logarithms when only one of the two kinematic variables approaches the threshold, and providing explicit coefficients ready for phenomenological application.