Target-Mass Corrections in the OPE Sum-Rule Approach to Quarkonium-Nucleon Interactions with Global-Fit PDFs: an $x$-Resolved Analysis

This paper presents an $x$-resolved analysis of target-mass corrections in quarkonium-nucleon interactions using modern global-fit PDFs, revealing how the interplay between kinematic weights and the specific $x$-distribution of gluon moments across different PDF sets determines the magnitude of these corrections in sum-rule predictions.

The Gist

Imagine you are trying to predict how a tiny, heavy marble (a quarkonium, like a J/ψ particle) bounces off a large, soft pillow (a nucleon, like a proton or neutron).

In the world of particle physics, we usually treat the pillow as if it were weightless and perfectly still. We assume the marble is so small and the collision so fast that the pillow's own weight doesn't matter. This is the "standard recipe" physicists have used for decades.

However, this paper argues that for heavy marbles, the pillow's weight does matter. The authors are revisiting this recipe using the most modern, high-precision maps of how the pillow is made (called PDFs or Parton Distribution Functions) to see exactly how the pillow's weight changes the bounce.

Here is a breakdown of what they did, using simple analogies:

1. The Old Recipe vs. The New Map

• The Old Recipe: In the 1990s, physicists used a rough sketch of the pillow's internal structure. They calculated the bounce but ignored the fact that the pillow has mass. They assumed the pillow was just a cloud of weightless dust.
• The New Map: Today, we have incredibly detailed GPS maps of the pillow (called ABMP16, MSHT20, CT18, and NNPDF4.0). These maps tell us exactly where the "dust" (gluons) is concentrated: is it in the center? Is it on the edges? Is it sparse or dense?
• The Goal: The authors didn't just want to swap the old sketch for the new map. They wanted to trace the entire journey of the calculation to see exactly where the "weight of the pillow" changes the outcome.

2. The "Weight" Problem (Target-Mass Corrections)

Imagine you are pushing a shopping cart.

• Scenario A: The cart is empty (massless). You push it, and it accelerates instantly.
• Scenario B: The cart is full of bricks (massive). You push it, and it resists.

In particle physics, the "Target-Mass Correction" (TMC) is the math that accounts for the fact that the pillow (the nucleon) isn't weightless.

• The Twist: Usually, when we add weight to a calculation, we just add a small "correction factor" at the very end.
• The Discovery: This paper shows that the weight doesn't just add a small factor at the end. It changes the shape of the calculation all along the way. It acts like a filter that suppresses certain parts of the pillow's internal structure more than others.

3. The "x-Resolved" Analysis (Zooming In)

This is the most creative part of the paper. Instead of looking at the whole pillow at once, the authors sliced the pillow into different zones based on how much "momentum" (energy) the dust inside carries. Let's call these zones:

• Small-x: The fluffy, light dust on the very outside of the pillow.
• Large-x: The heavy, dense bricks near the core.

They asked: "When we account for the pillow's weight, which zone gets suppressed the most?"

The Finding:
The "weight correction" acts like a heavy hand pressing down on the large-x (heavy brick) zones.

• If your map of the pillow says the heavy bricks are concentrated in the center (Large-x), the "weight correction" will crush the predicted bounce significantly.
• If your map says the heavy bricks are spread out, the effect is different.
• Key Insight: The final result depends on both the physics of the weight (the universal rule) AND the specific shape of the map (the PDF) you are using.

4. The "Direct Convolution" (Avoiding the Shortcut)

In the past, to calculate the bounce, physicists used a "shortcut" (a mathematical guess called a parametric ansatz). They would measure a few points and draw a smooth curve through them.

• The Problem: With modern, high-precision maps, this shortcut breaks. It creates a "spike" in the prediction that looks nothing like reality (like predicting the marble bounces 100 times higher than it should).
• The Solution: The authors stopped using the shortcut. Instead, they used a Direct Convolution.
  • Analogy: Instead of guessing the shape of the hill, they walked every single step of the path, measuring the ground at every inch. They let the math flow naturally from the "dust" distribution to the final "bounce" without forcing it into a pre-made shape.

5. The Results: What Happened?

When they ran the numbers with the new maps and the "no-shortcut" method:

1. Near the Threshold (The Slow Bounce): When the marble hits the pillow slowly (near the minimum energy needed to make a bounce), the "weight correction" is huge. It suppresses the bounce probability by about 40%.
   • Why? Because at low speeds, the heavy bricks (Large-x) dominate the interaction, and the weight correction crushes them.
2. At High Speeds: When the marble hits very fast, the effect disappears. The pillow acts almost like it's weightless again.
3. The Map Matters: Different maps (ABMP16 vs. NNPDF4.0) gave slightly different results, but they all agreed on the main point: The weight of the target matters a lot near the start of the collision.

Summary

This paper is like a mechanic who says, "We used to think the car's engine worked the same way whether the car was empty or full of passengers. But now that we have better sensors, we see that the passengers (the target mass) actually change how the gears turn, especially when you start driving slowly."

They didn't just update the numbers; they built a transparent, step-by-step model showing exactly where and why the passengers change the ride, proving that ignoring the target's mass leads to a wrong picture of how heavy particles interact with matter.

Technical Summary

1. Problem Statement

The interaction between heavy quarkonia (such as $J/\psi$ and $\Upsilon$) and nucleons is a critical probe of QCD dynamics, bridging perturbative and non-perturbative scales. This interaction is typically modeled using the Operator Product Expansion (OPE) and multipole expansion, where the short-distance dynamics are encoded in Wilson coefficients and the long-distance structure in nucleon matrix elements of gluonic operators.

Key Issues Addressed:

• Target-Mass Corrections (TMC): In original analyses, trace terms in nucleon matrix elements were often neglected (assuming $m_N^2/\epsilon_0^2 \ll 1$). However, for $J/\psi$, this ratio is not negligible ($\tau \approx 8.6$), necessitating a consistent treatment of the full traceless tensor structure of gluonic operators.
• Outdated PDF Inputs: Previous sum-rule analyses relied on parton distribution functions (PDFs) from the 1990s. Modern global fits (ABMP16, MSHT20, CT18, NNPDF4.0) incorporate vastly more data (DIS, collider, electroweak) and offer improved constraints on the gluon distribution across small, intermediate, and large-$x$ regions.
• Lack of Mechanistic Insight: Existing studies often treated TMC as a final correction to the cross-section or focused on phenomenological fits. There was a lack of analysis resolving how TMC effects propagate from the gluon distribution $g(x, Q)$ through Mellin moments to the final cross-section $\sigma_{\Phi N}(s)$.
• Pathology in Parametric Fits: Attempts to fit simple parametric ansatzes (e.g., $\sigma \propto (s-s_{th})^a$) to modern PDF moments resulted in unphysically large threshold exponents ($a \approx 17\text{--}20$) and cross-section shapes diverging from experimental data (HERA).

2. Methodology

The author employs a rigorous x-resolved analysis to trace the full chain of the calculation:
$$g(x, Q) \longrightarrow A_n(Q) \longrightarrow \text{Sum Rules} \longrightarrow \sigma_{\Phi N}(s)$$

A. Theoretical Framework

• OPE Amplitude: The forward $\Phi N$ scattering amplitude is expanded in terms of local twist-2 gluonic operators.
• Mellin Moments: Defined as $A_n(Q^2) = \int_0^1 dx \, x^{n-2} g(x, Q^2)$. The analysis restricts $n \ge 4$ to avoid infrared divergences present in the $n=2$ moment for modern PDFs with negative small-$x$ exponents.
• Target-Mass Corrections (TMC): Instead of neglecting trace terms, the author uses the full traceless tensor structure. This replaces the standard moments $A_n$ with TMC-weighted moments $\tilde{A}_n$:
  $$\tilde{A}_n(Q^2) = \int_0^1 dx \, x^{n-2} g(x, Q^2) T_n(x)$$
  where $T_n(x)$ is a generalized hypergeometric function ($_3F_2$) acting as a kinematic weight dependent on the target mass $m_N$ and binding energy $\epsilon_0$.

B. x-Resolved Moment Analysis

To understand the origin of TMC effects, the paper introduces:

• Moment Densities: $w_n(x) = x^{n-2}g(x)$ and $\tilde{w}_n(x) = x^{n-2}g(x)T_n(x)$.
• Log-x Profiles: Visualization of which $x$-regions dominate specific moments.
• Partial Moments: Decomposition of integrals into finite $x$-windows to determine how TMC suppresses contributions from small, intermediate, and large-$x$ regions differently for each moment order $n$.

C. Cross-Section Reconstruction

• Direct Convolution: Instead of fitting parametric ansatzes to a few moments, the cross-section is computed via a direct partonic convolution integral derived from the optical theorem:
  $$\sigma_{\Phi N}(\lambda) = C \int_{\epsilon_0/\lambda}^1 \frac{dx}{x} g(x, Q) K(x, \lambda)$$
  where the kernel $K$ encodes the correct threshold behavior $\sigma \propto (s-s_{th})^{3/2}$.
• TMC Implementation: For TMC, the kernel is modified by summing the series with the $T_n(x)$ weights. The series is truncated at $n_{max}=10$, which is shown to converge rapidly (error $< 10^{-5}$).
• PDF Inputs: Analytic approximations and tabulated values from ABMP16, MSHT20, CT18, and NNPDF4.0 are used at scales $Q = 10, 100, 1000$ GeV.

3. Key Contributions

1. Mechanistic Decomposition: The paper explicitly separates the TMC effect into three layers:
   • The universal local kinematic suppression encoded in $T_n(x)$.
   • The PDF-dependent redistribution of moment support across $x$-regions.
   • The final observable distortion of the cross-section.
2. Resolution of Parametric Pathology: By using direct convolution rather than fitting moments to a power-law ansatz, the study avoids the unphysical threshold exponents ($a \approx 17\text{--}20$) found in previous works, yielding a cross-section shape consistent with HERA and high-energy data.
3. x-Window Diagnostics: The analysis identifies that higher Mellin moments ($n \ge 6$) are dominated by the large-$x$ tail of the gluon distribution, where different global PDF sets diverge most significantly.
4. NNPDF4.0 Integration: The work addresses the specific challenges of using Monte Carlo-based NNPDF4.0 by supplementing analytic fits with direct tabulated LHAPDF inputs to ensure accuracy in the large-$x$ region.

4. Results

• TMC Suppression Magnitude:
  • The TMC effect is most pronounced near the production threshold ($\sqrt{s} \approx 4.04$ GeV).
  • The ratio of TMC-corrected to uncorrected cross-section, $R_{TMC}$, drops to 0.60–0.63 (a ~40% suppression) near threshold, regardless of the PDF set.
  • At higher energies, the integration region shifts to small $x$ where $T_n(x) \to 1$, causing $R_{TMC}$ to approach unity.
• Moment-Level Suppression ($Q=10$ GeV):
  • The integrated suppression ratio $R_n = \tilde{A}_n/A_n$ varies significantly with $n$:
    • $n=4$: $R_4 \approx 0.53\text{--}0.62$
    • $n=6$: $R_6 \approx 0.09\text{--}0.16$
    • $n=8$: $R_8 \approx 0.009\text{--}0.027$
  • Higher moments are suppressed much more strongly because their support is shifted to large $x$, where the kinematic weight $T_n(x)$ is smallest.
• PDF Dependence:
  • While the magnitude of suppression is universal near threshold, the absolute cross-section values and the shape of the distribution depend on the specific PDF set.
  • Differences in the large-$x$ gluon shape among ABMP16, MSHT20, CT18, and NNPDF4.0 lead to variations in the cross-section spread, which persists through the convolution to the observable level.
• Scale Stability: The TMC suppression ratios are largely stable when evolving from $Q=10$ GeV to $Q=100$ GeV, confirming that the effect is driven primarily by the kinematic weight rather than DGLAP evolution.

5. Significance

This work provides a transparent reanalysis of quarkonium-nucleon interactions that clarifies the role of modern PDF information.

• Theoretical Clarity: It demonstrates that TMC should not be viewed as a simple global suppression factor but as an $x$-dependent reweighting of twist-2 moments.
• Phenomenological Reliability: By avoiding parametric ansatzes and using direct convolution, the study produces cross-section predictions that are physically consistent with experimental data, resolving previous discrepancies in threshold behavior.
• Future Applications: The framework establishes a robust method for incorporating finite target-mass effects in OPE sum rules, which is essential for precise predictions in heavy quarkonium physics, particularly for upcoming experiments at the Electron-Ion Collider (EIC) and the LHC.

In summary, the paper moves beyond simple updates of PDF inputs to provide a deep, x-resolved understanding of how target-mass corrections propagate from the gluon distribution to the observable cross-section, identifying the large-$x$ gluon region as the critical domain for these corrections.