REVIEW 3 major objections 3 minor 69 references
Sivers Asymmetry in Photoproduction of $J/\psi$ and Jet at the EIC
T0 review · 3 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Almost back-to-back $J/\psi$-jet photoproduction at the EIC is argued to give a gluon-dominated Sivers asymmetry of a few percent, making the channel a practical probe of the gluon Sivers function.
desk verdict A transparent, useful pheno estimate of a few-percent gluon Sivers asymmetry for J/psi+jet at the EIC, with a real caveat about the unspecified f/d mixture in color-octet channels. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the weighted Sivers asymmetry $A_N^{\sin(\phi_q)}$, the ratio of the $\sin\phi_q$-modulated difference of cross sections for transversely polarized protons to the unpolarized cross section, where $\phi_q$ is the azimuthal angle of the total transverse momentum $q_\perp=P_{\Psi\perp}+P_{j\perp}$ of the $J/\psi$-jet pair. The argument is carried by the generalized parton model, which assumes TMD factorization and retains the initial parton transverse momentum in the hard-scattering amplitudes, applied to the subprocesses $\gamma g \to J/\psi g$ and $\gamma q \to J/\psi q$. Quarkonium production is treated with NRQCD, with color-singlet $^3S_1^{(1)}$ and color-octet $^3S_1^{(8)}$, $^1S_0^{(8)}$, $^3P_J^{(8)}$ channels; the gluon Sivers function enters through a Gaussian model whose $x$-dependent normalization is fixed by available fits (DGLAP) or by averaging the known quark Sivers normalizations (TMD-a/TMD-b). The TMD-evolution treatment adds standard $b_\perp$-space evolution with a $b_*$ prescription and nonperturbative Gaussian factors, and this difference in treatment is what produces the sign and size differences between the two schemes.
What would settle it
A measurement of the $\sin\phi_q$-weighted asymmetry in $e p^\uparrow \to J/\psi+\mathrm{jet}+X$ at the EIC at $\sqrt{s}=45$ GeV with $K_\perp=3$ GeV, $z=0.3$, and $0<q_\perp<1$ GeV would settle the claim: the prediction is a positive few-percent gluon-dominated asymmetry under DGLAP evolution and a negative few-percent asymmetry under the exploratory TMD evolution, with quark-initiated contributions suppressed. Observing a very different size, sign, or $q_\perp$ dependence, or a significant quark contribution, would falsify the proposed interpretation.
Extended reading notes
Core claim
The central claim is that the back-to-back $J/\psi+\mathrm{jet}$ photoproduction channel at the EIC is a viable gluon Sivers probe because the asymmetry is both sizable and gluon-dominated. For $K_\perp=3$ GeV, $z=0.3$, and $0<q_\perp<1$ GeV, the weighted asymmetry reaches roughly 3%, 1%, and 6% at $\sqrt{s}=45$ GeV for the SIDIS1, SIDIS2, and SIDIS3 DGLAP-evolved gluon Sivers parametrizations, drops to about 2% at $\sqrt{s}=100$ GeV, and reaches up to 8% (4% at $\sqrt{s}=100$ GeV) in the TMD-evolved TMD-b scheme with a negative sign. The largest contributions come from the $^3S_1^{(1)}$ color-singlet and $^1S_0^{(8)}$ color-octet NRQCD states, and the result is insensitive to which published long-distance matrix-element set is used. The paper also reports that the sign of the asymmetry flips from positive to negative when TMD evolution is applied, a feature it attributes to the different gluon Sivers parametrizations used in the two schemes.
Load-bearing premise
The load-bearing premise is that a factorization scheme in which the initial parton's transverse momentum is kept in the hard scattering, called the generalized parton model here, is valid for heavy-quarkonium-plus-jet photoproduction; the paper states this assumption but does not prove it, and it does not specify which gluon Sivers component enters the color-octet diagrams.
Editorial extensions
If this is right
- If the central claim is right, the EIC can measure the gluon Sivers function through $J/\psi+\mathrm{jet}$ photoproduction without requiring the $J/\psi$ to be produced in the forward region.
- Quark- and antiquark-initiated subprocesses contribute negligibly in the considered kinematics, so the measured asymmetry can be interpreted as essentially a gluon Sivers signal.
- Both $^3S_1^{(1)}$ and $^1S_0^{(8)}$ NRQCD channels must be kept; their combined contribution, and the insensitivity to the long-distance matrix-element set, make the prediction stable across model choices.
- The sign of the asymmetry distinguishes the evolution treatments: positive and a few percent under DGLAP-evolved fits, negative and up to about 8% under the exploratory TMD-evolved parametrization.
- Higher $K_\perp$ or larger $z$ suppresses the gluon channel because the parton momentum fraction grows quadratically with $K_\perp$, so the optimal measurement region is $K_\perp\sim M_{J/\psi}$ and $z\sim0.3$.
Reading between the lines
- Because the sign flip between DGLAP and TMD-evolved treatments is driven by the unknown gluon Sivers parametrization rather than by the hard scattering itself, a precise measurement of both sign and $q_\perp$ shape could discriminate between evolution schemes and constrain nonperturbative TMD parameters.
- A dedicated gluon Sivers fit using this channel, rather than rescaling quark Sivers parameters, could test the transverse-momentum sum rule with gluon data: the TMD-a proxy satisfies the sum rule to about 1%, while TMD-b violates it by about 19%.
- The numerical estimates omit resolved-photon contributions, $\psi(2S)$ feed-down (about 15%), and $\chi_c$ decays (about 1%); a full experimental analysis would likely need these corrections to reach the few-percent precision claimed.
- The paper leaves open which combination of f-type and d-type gluon Sivers functions enters the color-octet diagrams; disentangling these process-dependent components is a natural follow-up.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript calculates the transverse-spin (Sivers) asymmetry for the quasi-real photoproduction process e p^↑ → J/ψ + jet + X at the future Electron-Ion Collider, using the generalized parton model (GPM) together with NRQCD for the J/ψ production mechanism. Color-singlet and color-octet intermediate states (3S1^1, 3S1^8, 1S0^8, 3P_J^8) are included, with matrix elements taken from the authors' earlier gluon-channel paper and from new quark-channel expressions in Appendix B. The asymmetry is evaluated both with DGLAP-evolved Sivers parametrizations (SIDIS1, SIDIS2, SIDIS3) and with an exploratory TMD-evolution scheme in which the gluon Sivers function is built from u/d quark Sivers parameters. The authors find a gluon-dominated asymmetry of a few percent for DGLAP sets, up to about 6% at √s=45 GeV, and up to about 8% in the TMD-evolved TMD-b scheme. They conclude that almost back-to-back J/ψ+jet photoproduction is a promising probe of the gluon Sivers function.
Significance. The observable is well motivated and complementary to existing and proposed gluon-Sivers probes: the additional jet makes the kinematics more differential and avoids restricting J/ψ to the forward region. The manuscript provides explicit analytic matrix elements, a clear tabulation of input parameters, and a maximized-asymmetry study that isolates the role of each NRQCD channel. If the calculation were fully justified, the predicted asymmetry of several percent would be a useful benchmark for EIC planning. The main caveats are that the prediction rests on TMD factorization in the GPM for color-octet quarkonium production, which is not proven for this process, and that the color (f-type versus d-type) process dependence of the gluon Sivers function is not implemented for the color-octet channels that are shown to be numerically important.
major comments (3)
- [Section II, Eq. (9) and Fig. 2] The paper introduces in Section I the fact that the gluon Sivers function for any process is generally a combination of two independent color structures, conventionally called f-type and d-type, but the calculation in Eq. (9) inserts the same ΔN f_{g/p↑}(x_g,q⊥) for every partonic channel. This is not a harmless simplification for the NRQCD color-octet channels 3S1^8, 1S0^8, and 3P_J^8, whose amplitudes in Appendix B and Ref. [29] involve color-octet operators; the relevant gluon Sivers function for such channels is generally a process-dependent linear combination of f- and d-type functions, and that combination is never identified. Since Fig. 2 shows that the 1S0^8 channel is one of the two dominant contributions, the quoted few-percent asymmetries are not actually predictions for the f-type Sivers function unless the f/d combination is specified and modeled. The sign and magnitude of the color-octet contributions could shift, and this directly affects the central claim that the process is a clean probe of the gluon Sivers function.
- [Section III, Eqs. (28) and Fig. 4] The TMD-evolved results do not use a gluon Sivers function extracted from gluon-sensitive data. Instead, Eqs. (28) construct a gluon Sivers function from u/d quark parameters via Ng=(Nu+Nd)/2 (TMD-a) or Ng=Nd (TMD-b), and the text notes that the TMD-b choice violates the Burkardt sum rule by about 19%. The up-to-8% negative asymmetry in Fig. 4 and the statement in the Conclusion that the asymmetry becomes negative when TMD evolution is incorporated are therefore driven by an exploratory, quark-inspired model rather than by a quantitative gluon fit. This is an interesting illustration, but the manuscript should label the TMD-evolved numbers as model estimates in the abstract and conclusion, and should not present the sign flip as a robust prediction.
- [Section II, after Eq. (2) and Section IV] The calculation assumes TMD factorization in the GPM, as the authors state explicitly after Eq. (2), but no factorization proof is given for ep → J/ψ + jet + X with NRQCD color-octet contributions, and Ref. [39] is cited only for the status in pp collisions. Since the entire numerical estimate and the conclusion that the process can determine the gluon Sivers function rely on this assumption, the manuscript should include a critical discussion of its range of validity for the color-octet channels in the q⊥ ≪ K⊥ kinematics. If the GPM is intended as a phenomenological model rather than a proven factorization, this should be stated more prominently, and the consequences for the interpretation of the asymmetry should be assessed.
minor comments (3)
- [Throughout] There are several typographical errors that should be corrected, including 'inreaction' in Section II, 'avavilable' and 'wthin' in Section III, and the axis labels in Figs. 2–4 that are rendered as 'q¦ HGeVL' and 'AN sin IfqM'.
- [Section II] The sentence following Eq. (2) that introduces the Weizsäcker-Williams distribution states Qmax² = 1 GeV², but it is not explained how this choice interacts with the quasi-real photon condition Q² ≈ 0; a brief justification would improve clarity.
- [Section IV] The text says the range 0 ≤ q⊥ ≤ 1 GeV is considered to satisfy |q⊥| ≪ |K⊥| for K⊥ = 3 GeV; since the ratio reaches 1/3, the authors should either justify that this is sufficiently small for the back-to-back approximation or restrict the plotted range further.
Circularity Check
No circularity: the asymmetry is a model convolution of externally fitted Sivers functions with fixed NRQCD hard matrix elements.
full rationale
The derivation chain is self-contained as a phenomenological estimate. The Sivers asymmetry is built from Eqs. (2)-(10) and (12)-(17) as a convolution of independently fitted TMD inputs (gluon Sivers functions from Refs. [7,8] fitted to RHIC data, quark Sivers functions from Ref. [13] fitted to SIDIS data) with NRQCD hard-scattering amplitudes, evaluated at fixed EIC kinematics. Nothing in the calculation is fitted to the J/psi+jet asymmetry itself, and the numerator and denominator share the same parametrizations, so the reported few-percent asymmetries are not forced by construction. The only author-overlapping input is Ref. [29], cited for the gamma g -> J/psi g matrix elements; those are fixed perturbative NRQCD amplitudes, not outputs of the present paper, and the paper does not invoke any uniqueness theorem or fitted parameter to make its central prediction. The paper explicitly flags its limitations in Section II ('We have assumed TMD factorization in the GPM model') and Section III ('as so far a fit for the gluon Sivers function is not available in the TMD evolution approach... we use an exploratory approach'), but these are caveats about model validity rather than circular reductions. No quoted equation reduces to its own input by definition, so the appropriate finding is no significant circularity.
Assumptions & free parameters
free parameters (6)
- Na (gluon Sivers normalization) =
0.65 (SIDIS1), 0.05 (SIDIS2), 0.25 (SIDIS3)
- Sivers x-shape parameters alpha, beta =
alpha=2.8,0.8,0.6; beta=2.8,1.4,0.6
- Sivers Gaussian width parameter rho =
0.687, 0.576, 0.1
- Unpolarized TMD width <q_T^2> =
0.25 GeV^2 (SIDIS1/2), 1.0 GeV^2 (SIDIS3)
- TMD evolution nonperturbative parameters =
g2=0.16 GeV^2, bmax=1.5 GeV^-1, Q0=sqrt(2.4) GeV, <p_s^2>=0.282 GeV^2
- TMD-a/TMD-b gluon Sivers construction =
Nu=0.106, Nd=-0.163, alpha, beta, rho=0.38
assumptions (5)
- domain assumption TMD factorization (GPM) is valid for ep -> J/psi + jet photoproduction with NRQCD color octet contributions.
- domain assumption NRQCD factorization applies and the LDMEs of Ref. [64] are reliable.
- domain assumption The gluon Sivers function entering the asymmetry is the f-type distribution as parametrized in Refs. [7,8]; the d-type component is ignored.
- domain assumption The Weizsacker-Williams equivalent photon approximation describes the quasi-real photon flux.
- domain assumption Resolved photon, feed-down, and fragmentation contributions can be neglected in the chosen kinematics.
Cite this review
Pith. "Pith review of Sivers Asymmetry in Photoproduction of $J/\psi$ and Jet at the EIC." pith.science (2026). https://pith.science/paper/RJBY3JSV
@misc{pith2026190803698,
author = {Pith},
title = {Pith review of: Sivers Asymmetry in Photoproduction of $J/\psi$ and Jet at the EIC},
year = {2026},
howpublished = {\url{https://pith.science/paper/RJBY3JSV}},
note = {Machine review of arXiv:1908.03698}
}
abstract
We calculate the Sivers asymmetry in the photoproduction of almost back-to-back $J/\psi$-jet pair in the process $ep^\uparrow \to J/\psi+\mathrm{jet}+X$, which will be possible at the future planned electron-ion collider (EIC). We use the framework of generalized parton model (GPM), and NRQCD for calculating the $J/\psi$ production rate. We include contributions from both color singlet and color octate states in the asymmetry. We obtain sizable Sivers asymmetry that can be promising to determine the gluon Sivers function. We also investigate the effect of TMD evolution on the asymmetry.
Figures
Reference graph
Works this paper leans on
- [29]
-
[39]
R. M. Godbole, A. Misra, A. Mukherjee, and V. S. Rawoot, Phys. Rev. D85, 094013 (2012), 1201.1066
arXiv 2012
-
[1]
D. W. Sivers, Phys. Rev. D41, 83 (1990)
1990
- [2]
- [3]
-
[4]
A. Airapetian et al. (HERMES), Phys. Rev. Lett. 94, 012002 (2005), hep-ex/0408013
arXiv 2005
- [5]
-
[6]
M. Anselmino, M. Boglione, U. D’Alesio, F. Murgia, and A. Prokudin, JHEP 04, 046 (2017), 1612.06413
arXiv 2017
Show all 69 references
- [7]
-
[8]
D’Alesio, C
U. D’Alesio, C. Flore, F. Murgia, C. Pisano, and P. Taels, Phys. Rev. D99, 036013 (2019), 1811.02970
2019 arXiv
-
[10]
S. M. Aybat and T. C. Rogers, Phys. Rev. D83, 114042 (2011), 1101.5057
2011 arXiv
- [11]
-
[12]
M. G. Echevarria, I. Scimemi, and A. Vladimirov, JHEP 09, 004 (2016), 1604.07869
2016 arXiv
-
[13]
M. G. Echevarria, A. Idilbi, Z.-B. Kang, and I. Vitev, Phys. Rev. D89, 074013 (2014), 1401.5078
2014 arXiv
- [14]
-
[15]
S. J. Brodsky, F. Fleuret, C. Hadjidakis, and J. P. Lansberg, Phys. Rept. 522, 239 (2013), 1202.6585
2013 arXiv
-
[16]
Kikoa, M
D. Kikoa, M. G. Echevarria, C. Hadjidakis, J.-P. Lansberg, C. Lorc, L. Massacrier, C. M. Quintans, A. Signori, and B. Trzeciak, Few Body Syst. 58, 139 (2017), 1702.01546
2017 arXiv
-
[17]
Trzeciak, C
B. Trzeciak, C. Da Silva, E. G. Ferreiro, C. Hadjidakis, D. Kikola, J. P. Lansberg, L. Massacrier, J. Seixas, A. Uras, and Z. Yang, Few Body Syst. 58, 148 (2017), 1703.03726
2017 arXiv
-
[18]
P. J. Mulders and J. Rodrigues, Phys. Rev. D63, 094021 (2001), hep-ph/0009343
2001 arXiv
-
[19]
gives a constraint on the GSF. However, in order to implement the contstraint from the Burkardt sum rule contribution from all quark favours need to be included, and sea quark Sivers function is still not well constrained. Assuming contributions only from u andd quarks and glu...
-
[20]
Anselmino, M
M. Anselmino, M. Boglione, U. D’Alesio, A. Kotzinian, S. Melis, F. Murgia, A. Prokudin, and C. Turk, Eur. Phys. J. A39, 89 (2009), 0805.2677
2009 arXiv
- [21]
-
[22]
D. Boer, P. J. Mulders, and F. Pijlman, Nucl. Phys. B667, 201 (2003), hep-ph/0303034
2003 arXiv
-
[23]
The GSF for any process, in general can be written in terms of two independent functions, one of them has an operator structure that is C-even, and the other, C-odd
as well as COMPASS [24] seem to favour the sign change, however more data are needed [6] . The GSF for any process, in general can be written in terms of two independent functions, one of them has an operator structure that is C-even, and the other, C-odd. In the literature, t...
1908 arXiv
-
[24]
S. J. Brodsky, D. S. Hwang, and I. Schmidt, Nucl. Phys. B642, 344 (2002), hep-ph/0206259
2002 arXiv
-
[25]
Adamczyk et al
L. Adamczyk et al. (STAR), Phys. Rev. Lett. 116, 132301 (2016), 1511.06003
2016 arXiv
-
[26]
Aghasyan et al
M. Aghasyan et al. (COMPASS), Phys. Rev. Lett. 119, 112002 (2017), 1704.00488
2017 arXiv
-
[27]
C. J. Bomhof and P. J. Mulders, JHEP 02, 029 (2007), hep-ph/0609206
2007 arXiv
-
[28]
M. G. A. Buffing, A. Mukherjee, and P. J. Mulders, Phys. Rev. D88, 054027 (2013), 1306.5897
2013 arXiv
-
[30]
Lansberg, C
J.-P. Lansberg, C. Pisano, and M. Schlegel, Nucl. Phys. B920, 192 (2017), 1702.00305
2017 arXiv
-
[31]
Rajesh, R
S. Rajesh, R. Kishore, and A. Mukherjee, Phys. Rev. D98, 014007 (2018), 1802.10359
2018 arXiv
- [32]
-
[33]
D’Alesio, F
U. D’Alesio, F. Murgia, C. Pisano, and P. Taels, Phys. Rev. D96, 036011 (2017), 1705.04169
2017 arXiv
- [34]
- [35]
-
[36]
Matouek (COMPASS), J
J. Matouek (COMPASS), J. Phys. Conf. Ser. 678, 012050 (2016)
2016
- [37]
-
[38]
D’Alesio, C
U. D’Alesio, C. Flore, and F. Murgia, Phys. Rev. D95, 094002 (2017), 1701.01148
2017 arXiv
-
[40]
R. M. Godbole, A. Misra, A. Mukherjee, and V. S. Rawoot, Phys. Rev. D88, 014029 (2013), 1304.2584
2013 arXiv
-
[41]
M. G. Echevarria (2019), 1907.06494
2019 arXiv
-
[42]
C. E. Carlson and R. Suaya, Phys. Rev. D14, 3115 (1976)
1976
-
[43]
E. L. Berger and D. L. Jones, Phys. Rev. D23, 1521 (1981)
1981
-
[44]
Baier and R
R. Baier and R. Ruckl, Phys. Lett. 102B, 364 (1981)
1981
-
[45]
Baier and R
R. Baier and R. Ruckl, Nucl. Phys. B201, 1 (1982)
1982
-
[46]
Braaten and S
E. Braaten and S. Fleming, Phys. Rev. Lett. 74, 3327 (1995), hep-ph/9411365
1995 arXiv
-
[47]
P. L. Cho and A. K. Leibovich, Phys. Rev. D53, 150 (1996), hep-ph/9505329
1996 arXiv
-
[48]
P. L. Cho and A. K. Leibovich, Phys. Rev. D53, 6203 (1996), hep-ph/9511315
1996 arXiv
-
[49]
G. P. Lepage, L. Magnea, C. Nakhleh, U. Magnea, and K. Hornbostel, Phys. Rev. D46, 4052 (1992), hep-lat/9205007
1992 arXiv
-
[50]
Abe et al
F. Abe et al. (CDF), Phys. Rev. Lett. 79, 572 (1997)
1997
- [51]
- [52]
-
[53]
F. D. Aaron et al. (H1), Eur. Phys. J. C68, 401 (2010), 1002.0234
2010 arXiv
- [54]
- [55]
-
[56]
Frixione, M
S. Frixione, M. L. Mangano, P. Nason, and G. Ridolfi, Phys. Lett. B319, 339 (1993), hep-ph/9310350
1993 arXiv
-
[57]
D. Boer, P. J. Mulders, C. Pisano, and J. Zhou, JHEP 08, 001 (2016), 1605.07934
2016 arXiv
-
[58]
Bacchetta, U
A. Bacchetta, U. D’Alesio, M. Diehl, and C. A. Miller, Phys. Rev. D70, 117504 (2004), hep-ph/0410050
2004 arXiv
-
[59]
Bacchetta, C
A. Bacchetta, C. Bomhof, U. D’Alesio, P. J. Mulders, and F. Murgia, Phys. Rev. Lett. 99, 212002 (2007), hep-ph/0703153
2007 arXiv
-
[60]
Collins, Foundations of perturbative QCD (Cambridge University Press, 2013)
J. Collins, Foundations of perturbative QCD (Cambridge University Press, 2013)
2013
-
[61]
S. M. Aybat, J. C. Collins, J.-W. Qiu, and T. C. Rogers, Phys. Rev. D85, 034043 (2012), 1110.6428
2012 arXiv
-
[62]
Ji, J.-W
X. Ji, J.-W. Qiu, W. Vogelsang, and F. Yuan, Phys. Rev. Lett. 97, 082002 (2006), hep-ph/0602239
2006 arXiv
-
[63]
Kouvaris, J.-W
C. Kouvaris, J.-W. Qiu, W. Vogelsang, and F. Yuan, Phys. Rev. D74, 114013 (2006), hep-ph/0609238
2006 arXiv
- [64]
- [65]
-
[66]
Chao, Y.-Q
K.-T. Chao, Y.-Q. Ma, H.-S. Shao, K. Wang, and Y.-J. Zhang, Phys. Rev. Lett. 108, 242004 (2012), 1201.2675
2012 arXiv
-
[67]
Li and L.-s
Y.-d. Li and L.-s. Liu, Commun. Theor. Phys. 29, 99 (1998)
1998
-
[68]
Butenschoen and B
M. Butenschoen and B. A. Kniehl, Phys. Rev. Lett. 104, 072001 (2010), 0909.2798
2010 arXiv
-
[69]
Artoisenet, J
P. Artoisenet, J. M. Campbell, F. Maltoni, and F. Tramontano, Phys. Rev. Lett. 102, 142001 (2009), 0901.4352
2009 arXiv
-
[70]
Buckley, J
A. Buckley, J. Ferrando, S. Lloyd, K. Nordstrm, B. Page, M. Rfenacht, M. Schnherr, and G. Watt, Eur. Phys. J. C75, 132 (2015), 1412.7420
2015 arXiv
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.