pith. sign in

arxiv: 2606.18854 · v1 · pith:H2AVNPG3new · submitted 2026-06-17 · ✦ hep-ph · hep-th· nucl-th

Diffractive Production of Heavy Quarkonia at the Electron Ion Collider

Mohammad Yousuf Jamal , Ajaharul Islam , Aritra Bandyopadhyay , Santosh K. Das This is my paper

Pith reviewed 2026-06-26 20:34 UTC · model grok-4.3

classification ✦ hep-ph hep-thnucl-th
keywords diffractive productionheavy quarkoniaElectron-Ion Collidertwo-gluon kernellight-cone wave functionscharmoniumbottomoniumtransverse momentum dependence
0
0 comments X

The pith

Retaining the full transverse-momentum dependence of the two-gluon kernel improves the description of radially excited heavy quarkonia in diffractive production at the EIC.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper calculates diffractive photo- and electroproduction of J/ψ, ψ(2S), and Υ(nS) states at EIC energies while keeping the complete transverse-momentum dependence in the hard two-gluon kernel instead of using the small-size color-dipole expansion. Quarkonia light-cone wave functions are constructed from Cornell-potential solutions of the Schrödinger equation and normalized to measured leptonic widths, then combined with a modern collinear gluon distribution. After benchmarking the approach against the full set of HERA charmonium cross-section ratio data, the calculation shows that the full ℓ_t treatment systematically improves agreement for radially excited states relative to the dipole limit. The authors supply bottomonium ratio predictions for EIC kinematics and identify the kinematic windows where the difference between the two treatments is largest.

Core claim

By retaining the full ℓ_t dependence of the hard two-gluon kernel in the production amplitude for S-wave heavy quarkonia, rather than expanding the impact-parameter Bessel kernel in the small-size color-dipole limit, and employing light-cone wave functions built from Cornell-potential Schrödinger solutions normalized to leptonic widths together with a modern collinear gluon distribution, the framework improves the description of radially excited states after validation on HERA charmonium data and yields consistent bottomonium cross-section ratio predictions in EIC kinematics, with the largest gains in identified kinematic windows.

What carries the argument

The full transverse-momentum (ℓ_t) resolved two-gluon kernel in the production amplitude, which avoids the Bessel expansion of the impact-parameter kernel used in the dipole limit and thereby captures improved overlap with the radially excited quarkonia wave functions.

If this is right

  • Cross-section ratios for radially excited states to ground states will be larger than those obtained in the dipole limit.
  • The difference between full ℓ_t and dipole treatments reaches its maximum in specific kinematic windows at EIC energies.
  • A consistent set of bottomonium production ratios is obtained once the framework is benchmarked on HERA charmonium data.
  • The same full-kernel approach applies to both photoproduction and electroproduction channels.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • EIC data on bottomonium ratios could be used to test whether the full ℓ_t treatment is required for accurate extractions of gluon distributions from diffractive processes.
  • The method could be applied to other vector-meson production channels where measurements of excited states are available at similar energies.
  • Re-analysis of existing HERA excited-state data with this kernel might reduce tensions between theory and experiment in those ratios.

Load-bearing premise

The light-cone wave functions from Cornell-potential Schrödinger solutions, normalized to leptonic widths and combined with a modern collinear gluon distribution, accurately represent the production dynamics when the full transverse-momentum dependence of the hard two-gluon kernel is retained.

What would settle it

EIC measurements of the ψ(2S) to J/ψ cross-section ratio in the kinematic windows where the paper predicts the largest difference that instead match the dipole-limit values rather than the full ℓ_t calculation would falsify the claim of systematic improvement.

Figures

Figures reproduced from arXiv: 2606.18854 by Ajaharul Islam, Aritra Bandyopadhyay, Mohammad Yousuf Jamal, Santosh K. Das.

Figure 1
Figure 1. Figure 1: FIG. 1: Schematic of exclusive diffractive vector-meson [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: (a) Cornell potential of Eq. (2) for [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Ratio [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: presents the bottomonium ratios σ(Υ(nS))/σ(Υ(1S)) for n = 2, 3, for which we are not aware of any existing experimental data. The heavier b mass reduces typical dipole sizes, so the node￾induced cancellations are milder than for charmonium, yet the full and dipole results still differ appreciably over a broad Q2 range, as seen in Fig. 4a. The full result is smoother and more monotonic, while the dipole lim… view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: (a) Photon LC wave function [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: The [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: shows the full z-dependent amplitude A(z; Q2 , W) [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
read the original abstract

We study diffractive photo- and electroproduction of the $S$-wave heavy quarkonia $J/\psi$, $\psi(2S)$, and $\Upsilon(nS)$ at energies relevant for the Electron-Ion Collider (EIC). The production amplitude is evaluated while retaining the full transverse-momentum ($\ell_t$) dependence of the hard two-gluon kernel, that is, without expanding the impact-parameter Bessel kernel as is done in the small-size color-dipole limit. The quarkonia light-cone wave functions are built from Cornell-potential solutions of the Schr\"odinger equation, normalized to the measured leptonic widths, and combined with a modern collinear gluon distribution. After benchmarking the framework against the full set of HERA charmonium cross-section ratio data, we provide a consistent set of bottomonium cross-section ratio predictions in EIC kinematics. We find that the full $\ell_t$-resolved treatment systematically improves the description of the radially excited states relative to the leading dipole limit, and we identify the kinematic windows where this difference is largest.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript calculates diffractive photo- and electroproduction of S-wave heavy quarkonia (J/ψ, ψ(2S), Υ(nS)) at EIC energies. The amplitude retains the full transverse-momentum dependence of the two-gluon kernel (no small-dipole Bessel expansion). Light-cone wave functions are obtained from Cornell-potential Schrödinger solutions, normalized to measured leptonic widths, and combined with a modern collinear gluon PDF. After benchmarking the framework against the complete HERA charmonium cross-section ratio data set, the authors supply bottomonium ratio predictions and state that the full ℓ_t treatment systematically improves the description of radially excited states relative to the dipole limit, with the largest differences in identifiable kinematic windows.

Significance. If the central claim holds, the work supplies a consistent set of EIC predictions that highlight the practical importance of the unexpanded kernel for excited-state ratios. The explicit benchmarking against the full HERA data set is a clear strength, as is the identification of kinematic windows where the difference is largest. The result is of direct relevance for EIC planning in the heavy-quarkonium sector.

major comments (2)
  1. [Abstract and benchmarking section] Abstract and benchmarking paragraph: the central claim that the full ℓ_t treatment improves the description of radially excited states rests on the Cornell wave functions accurately capturing the radial dependence at dipole sizes weighted by the unexpanded Bessel kernel. Normalization is performed solely to leptonic widths (primarily fixing |R(0)|); for nodal states such as ψ(2S) this leaves the shape at larger r less constrained. Because the same wave-function model is used both to generate the predictions and to benchmark against HERA data, the reported improvement could be an artifact of the potential choice rather than a robust feature of the ℓ_t treatment.
  2. [Method / wave-function section] Wave-function construction (method section): no sensitivity study is described that varies the potential parameters or substitutes alternative wave functions (e.g., from other quarkonium potentials or lattice inputs) while keeping the full kernel fixed. Such a test is required to establish that the reported improvement for excited states survives changes in the radial shape.
minor comments (2)
  1. [Numerical implementation] Notation: the precise numerical implementation of the full ℓ_t kernel (integration limits, regularization of the Bessel function at large arguments) should be stated explicitly so that the results can be reproduced.
  2. [Results figures] Figure clarity: the plots comparing full-kernel vs. dipole-limit ratios would benefit from explicit indication of the HERA data points used in the benchmark and the kinematic cuts applied.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We respond to the major comments point by point below.

read point-by-point responses

  1. Referee: [Abstract and benchmarking section] Abstract and benchmarking paragraph: the central claim that the full ℓ_t treatment improves the description of radially excited states rests on the Cornell wave functions accurately capturing the radial dependence at dipole sizes weighted by the unexpanded Bessel kernel. Normalization is performed solely to leptonic widths (primarily fixing |R(0)|); for nodal states such as ψ(2S) this leaves the shape at larger r less constrained. Because the same wave-function model is used both to generate the predictions and to benchmark against HERA data, the reported improvement could be an artifact of the potential choice rather than a robust feature of the ℓ_t treatment.

    Authors: We thank the referee for highlighting this subtlety. Importantly, the full ℓ_t calculation and the dipole-limit calculation are performed with identical Cornell wave functions (and the same gluon PDF). The reported improvement in the HERA charmonium ratios for excited states is therefore a direct consequence of the more accurate kernel treatment at the larger dipole sizes relevant to radially excited states, rather than a difference in wave-function modeling. The Cornell potential is a standard choice that reproduces the observed quarkonium spectrum, and normalization to leptonic widths follows the conventional procedure used in the literature. The HERA benchmarking, which includes excited-state ratios, provides empirical support for the framework in the kinematics of interest. Nevertheless, we agree that robustness to the radial shape merits explicit discussion; in the revised manuscript we will add a paragraph in the benchmarking section addressing the expected sensitivity based on existing comparisons of alternative quarkonium potentials in the literature. revision: partial

  2. Referee: [Method / wave-function section] Wave-function construction (method section): no sensitivity study is described that varies the potential parameters or substitutes alternative wave functions (e.g., from other quarkonium potentials or lattice inputs) while keeping the full kernel fixed. Such a test is required to establish that the reported improvement for excited states survives changes in the radial shape.

    Authors: We agree that a dedicated sensitivity study varying the potential or employing alternative wave functions (e.g., from other potentials or lattice QCD) would strengthen the claim of robustness. However, because the comparison between full ℓ_t and dipole limit uses exactly the same wave functions, the relative improvement for excited states is attributable to the kernel and is independent of the specific radial shape chosen. A comprehensive scan over multiple wave-function models would require substantial additional computational work outside the primary scope of demonstrating the effect of the unexpanded kernel. In the revised version we will expand the method section with a short discussion of this point, including references to prior studies that compare different light-cone wave functions for heavy quarkonia, while noting that the HERA data agreement constrains the present choice. revision: partial

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The derivation uses external inputs: Cornell-potential Schrödinger solutions normalized to measured leptonic widths (external data fixing |R(0)|), a modern collinear gluon PDF, and a numerical two-gluon kernel evaluated both in full ℓ_t form and in the dipole approximation. Benchmarking against the full HERA charmonium ratio dataset is a comparison to independent external measurements using the same fixed framework on both sides; no parameters are fitted inside the paper to the target ratios, and no self-citation chain or self-definitional loop is invoked to justify the central claim. The reported systematic improvement for radially excited states is therefore a numerical outcome of retaining the unexpanded Bessel kernel rather than a quantity forced by construction from the paper's own inputs.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The framework rests on established potential-model wave functions and literature gluon distributions; the only data-driven elements are the normalization constants taken from measured leptonic widths.

free parameters (1)
  • wave-function normalization constants
    Light-cone wave functions are normalized to measured leptonic widths as stated in the abstract.
axioms (2)
  • domain assumption Cornell potential solutions of the Schrödinger equation provide accurate light-cone wave functions for S-wave heavy quarkonia
    Used to build the wave functions for J/ψ, ψ(2S), and Υ(nS).
  • domain assumption A modern collinear gluon distribution is sufficient for the diffractive amplitude at the relevant kinematics
    Combined with the wave functions to evaluate the production amplitude.

pith-pipeline@v0.9.1-grok · 5731 in / 1547 out tokens · 39304 ms · 2026-06-26T20:34:28.448531+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

32 extracted references · 20 linked inside Pith

  1. [1]

    M. G. Ryskin, Z. Phys. C57, 89 (1993)

    1993

  2. [2]

    M. G. Ryskin, R. G. Roberts, A. D. Martin, and E. M. Levin, Z. Phys. C76, 231 (1997), arXiv:hep-ph/9511228

    Pith/arXiv arXiv 1997

  3. [3]

    Frankfurt, W

    L. Frankfurt, W. Koepf, and M. Strikman, Phys. Rev. D 57, 512 (1998), arXiv:hep-ph/9702216

    Pith/arXiv arXiv 1998

  4. [4]

    Chekanovet al.(ZEUS), Eur

    S. Chekanovet al.(ZEUS), Eur. Phys. J. C24, 345 (2002), arXiv:hep-ex/0201043

    Pith/arXiv arXiv 2002

  5. [5]

    Abdul Khaleket al., Nucl

    R. Abdul Khaleket al., Nucl. Phys. A1026, 122447 (2022), arXiv:2103.05419 [physics.ins-det]

    Pith/arXiv arXiv 2022

  6. [6]

    D. Boer, C. Pisano, and P. Taels, Phys. Rev. D103, 074012 (2021), arXiv:2102.00003 [hep-ph]

    arXiv 2021

  7. [7]

    Fanelliet al., Nucl

    C. Fanelliet al., Nucl. Instrum. Meth. A1047, 167748 (2023), arXiv:2205.09185 [physics.ins-det]

    arXiv 2023

  8. [8]

    A. D. Martin, C. Nockles, M. G. Ryskin, and T. Teubner, Phys. Lett. B662, 252 (2008), arXiv:0709.4406 [hep-ph]

    Pith/arXiv arXiv 2008

  9. [9]

    S. P. Jones, A. D. Martin, M. G. Ryskin, and T. Teubner, JHEP11, 085, arXiv:1307.7099 [hep-ph]

    Pith/arXiv arXiv

  10. [10]

    Lappi and H

    T. Lappi and H. Mantysaari, Phys. Rev. C87, 032201 (2013), arXiv:1301.4095 [hep-ph]

    Pith/arXiv arXiv 2013

  11. [11]

    Tu (STAR), PoSHardProbes2020, 100 (2021), arXiv:2009.04860 [nucl-ex]

    Z. Tu (STAR), PoSHardProbes2020, 100 (2021), arXiv:2009.04860 [nucl-ex]

    arXiv 2021

  12. [12]

    Nemchik, N

    J. Nemchik, N. N. Nikolaev, E. Predazzi, B. G. Zakharov, and V. R. Zoller, J. Exp. Theor. Phys.86, 1054 (1998), arXiv:hep-ph/9712469

    Pith/arXiv arXiv 1998

  13. [13]

    Wusthoff and A

    M. Wusthoff and A. D. Martin, J. Phys. G25, R309 (1999), arXiv:hep-ph/9909362

    Pith/arXiv arXiv 1999

  14. [14]

    Adloffet al.(H1), Phys

    C. Adloffet al.(H1), Phys. Lett. B541, 251 (2002), arXiv:hep-ex/0205107

    Pith/arXiv arXiv 2002

  15. [15]

    L. L. Frankfurt, M. F. McDermott, and M. Strikman, JHEP02, 002, arXiv:hep-ph/9812316

    Pith/arXiv arXiv

  16. [16]

    Hoyer and S

    P. Hoyer and S. Peigne, Phys. Rev. D61, 031501 (2000), arXiv:hep-ph/9909519

    Pith/arXiv arXiv 2000

  17. [17]

    J.Cepila, J.Nemchik, M.Krelina,andR.Pasechnik,Eur. Phys. J. C79, 495 (2019), arXiv:1901.02664 [hep-ph]

    Pith/arXiv arXiv 2019

  18. [18]

    M. A. Peredo and M. Hentschinski, Phys. Rev. D109, 014032 (2024), arXiv:2308.15430 [hep-ph]

    arXiv 2024

  19. [19]

    Suzuki, A

    K. Suzuki, A. Hayashigaki, K. Itakura, J. Alam, and T. Hatsuda, Phys. Rev. D62, 031501 (2000), arXiv:hep- ph/0005250

    arXiv 2000

  20. [20]

    Buchmuller and S

    W. Buchmuller and S. H. H. Tye, Phys. Rev. D24, 132 (1981)

    1981

  21. [21]

    Quigg and J

    C. Quigg and J. L. Rosner, Phys. Lett. B71, 153 (1977)

    1977

  22. [22]

    Abramowiczet al.(H1, ZEUS), Eur

    H. Abramowiczet al.(H1, ZEUS), Eur. Phys. J. C75, 580 (2015), arXiv:1506.06042 [hep-ex]

    Pith/arXiv arXiv 2015

  23. [23]

    Buckleyet al., Eur

    A. Buckleyet al., Eur. Phys. J. C75, 132 (2015)

    2015

  24. [24]

    Accardiet al., Eur

    A. Accardiet al., Eur. Phys. J. A52, 268 (2016), arXiv:1212.1701 [nucl-ex]

    Pith/arXiv arXiv 2016

  25. [25]

    M. B. Gay Ducati, F. Kopp, M. V. T. Machado, and S. Martins, Phys. Rev. D94, 094023 (2016), arXiv:1610.06647 [hep-ph]

    Pith/arXiv arXiv 2016

  26. [26]

    Abtet al.(ZEUS), JHEP12, 164, arXiv:2206.13343 [hep-ex]

    I. Abtet al.(ZEUS), JHEP12, 164, arXiv:2206.13343 [hep-ex]

    arXiv

  27. [27]

    Chekanovet al.(ZEUS), Nucl

    S. Chekanovet al.(ZEUS), Nucl. Phys. B695, 3 (2004), arXiv:hep-ex/0404008

    Pith/arXiv arXiv 2004

  28. [28]

    Aktaset al.(H1), Eur

    A. Aktaset al.(H1), Eur. Phys. J. C46, 585 (2006), arXiv:hep-ex/0510016

    Pith/arXiv arXiv 2006

  29. [29]

    Adloffet al.(H1), Eur

    C. Adloffet al.(H1), Eur. Phys. J. C3, 13 (1998), arXiv:hep-ex/9712020

    Pith/arXiv arXiv 1998

  30. [30]

    See Supplemental Material at [URL] for parameters, nu- merical details, the convergence analysis of the small-ℓtb expansion, and three-dimensional visualizationsof the in- termediate quantities

  31. [31]

    Alexaet al.(H1), Eur

    C. Alexaet al.(H1), Eur. Phys. J. C73, 2466 (2013), arXiv:1304.5162 [hep-ex]

    Pith/arXiv arXiv 2013

  32. [32]

    Navaset al.(Particle Data Group), Phys

    S. Navaset al.(Particle Data Group), Phys. Rev. D110, 030001 (2024). Supplemental Material This Supplemental Material collects the inputs, the nu- merical procedures, and the auxiliary results referenced in the main text. Section I lists the parameters and the diffractive-slope intercepts, Section II describes the numerical evaluation of the amplitude and...

    2024