{"id":"a15a4123-7d6e-46c5-be5e-368601ff8ba5","arxiv_id":"1908.03698","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Photoproduction of J/psi+jet at the EIC shows a gluon-dominated Sivers asymmetry of up to 6-8% in the models studied, making it a promising probe of the gluon Sivers function.","lead":"The authors predict the Sivers asymmetry for the photoproduction of back-to-back J/psi and jet pairs at the future Electron-Ion Collider. They find several-percent asymmetries that could help measure the poorly known gluon Sivers function.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Color-octet channels are evaluated with a single f-type gluon Sivers function, but the process-dependent f/d combination for those channels is never specified; this can shift or flip the few-percent asymmetry.","rationale":"The reader’s weakest-assumption statement already names the missing f/d specification for color-octet diagrams, and my concern sharpens that: it is not merely a caveat about TMD factorization, but an ambiguity about which gluon Sivers function is actually being probed at leading order in the GPM. The paper explicitly acknowledges the f/d distinction but never applies it to the CO matrix elements, while its own Figure 2 shows the 1S0^8 CO channel is one of the two largest contributions. This makes the handful-of-percent asymmetry and its sign underdetermined unless an additional model choice is made. The TMD-evolved sign flip is also driven by a quark-based proxy for the gluon Sivers function, so it cannot support a quantitative claim. These issues do not invalidate the paper as a phenomenological estimate: they mean the result is conditional on a stated (but unstated) process-dependence assumption, which is exactly the reader’s CONDITIONAL verdict. I therefore keep the verdict unchanged, while noting that the central claim should be framed as ‘sensitive to a process-dependent combination of the gluon Sivers functions’ until the f/d combination is specified and tested.","tokens_in":14127,"tokens_out":5880,"duration_ms":58377,"concrete_test":"At a fixed kinematic point (sqrt(s)=45 GeV, K_perp=3 GeV, z=0.3, q_perp=0.5 GeV), recompute the numerator of the asymmetry (Eq. 9) for the 1S0^8 and 3S1^8 channels, first with the f-type gluon Sivers function used in the paper, then with the process-dependent combination appropriate to the gamma-g -> J/psi-g color flow, using a model for the d-type component (e.g., the model in Ref. [28] or a Burkardt-sum-rule-consistent ansatz). If the color-octet contribution to A_N^{sin(phi_q)} shifts by more than about 1-2 percentage points or changes sign, the conclusion that this process cleanly determines the gluon Sivers function needs to be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the J/psi+jet asymmetry is a clean probe of the gluon Sivers function relies on inserting the same Delta N f_{g/p up arrow}(x_g, q_perp) in Eq. (9) for every partonic channel. The Introduction correctly notes that the gluon Sivers function has two independent color structures, f-type and d-type, but the calculation never specifies which combination enters the color-octet amplitudes in Appendix B (3S1^8, 1S0^8, 3P_J^8). For color-octet quarkonium production, the relevant gluon Sivers function is a process-dependent linear combination of f- and d-type functions, not the f-type function fitted in inclusive hadron production. Since Figure 2 shows the 1S0^8 color-octet channel is one of the two dominant contributions, the quoted asymmetry is not actually a prediction for “the” gluon Sivers function unless that combination is identified and modeled. In addition, the TMD-evolved results (Figure 4) use a gluon Sivers function built from u/d quark parameters via Eq. (28), not from a gluon fit, so the up-to-8% negative asymmetry is exploratory rather than a quantitative prediction.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":14397,"tokens_out":5744,"duration_ms":63316,"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":[{"comment":"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":"Section II, Eq. (9) and Fig. 2"},{"comment":"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":"Section III, Eqs. (28) and Fig. 4"},{"comment":"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.","section":"Section II, after Eq. (2) and Section IV"}],"minor_comments":[{"comment":"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":"Throughout"},{"comment":"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":"Section II"},{"comment":"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.","section":"Section IV"}],"recommendation":"major_revision","confidential_remarks":"The decisive issue for me is the f-type versus d-type color structure of the gluon Sivers function in the color-octet NRQCD channels. The authors are aware of the two color structures but do not implement or even discuss their process dependence for the CO amplitudes that make a large contribution to the asymmetry. This is fixable in a revision, either by computing the relevant f/d combination, by restricting the quantitative claim to the color-singlet channel, or by framing the color-octet result as a model uncertainty. The TMD-evolution section is clearly labeled as exploratory, but its sign-flip statement in the conclusion should be softened. The paper is otherwise a standard GPM/NRQCD estimate and could be publishable after these points are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is a straightforward, readable pheno paper, not a derivation. It extends the same group's earlier single-J/psi photoproduction calculation to the back-to-back J/psi+jet final state, adds both color-singlet and color-octet NRQCD contributions, and shows what TMD evolution does. The headline numbers are a Sivers asymmetry of a few percent, up to ~6% with DGLAP-evolved GSF fits and up to ~8% in an explicitly exploratory TMD-evolved scheme. The gluon channel dominates, so the process is indeed a plausible gluon-Sivers observable for the EIC.\n\nWhat the paper does well: the formalism is standard GPM+NRQCD and the presentation is clean. The quark-channel matrix elements are spelled out in Appendix B. The authors use published GSF fits where they exist and say plainly when they are making an exploratory construction (the TMD-a/TMD-b model from u and d quark Sivers parameters). They also state that TMD factorization is assumed, not proven. The maximized asymmetry plot is useful because it separates CS and CO contributions, which is exactly where the main caveat lives.\n\nThe soft spots, in proportion: the assumption of TMD factorization for color-octet J/psi production is a real limitation, but it is admitted and is standard for GPM estimates. The more specific problem is the f/d structure. The paper notes in the introduction that the gluon Sivers function is not one function; it is a process-dependent combination of f-type and d-type structures. Yet the calculation uses one f-type parametrization in Eq. (9) for every subprocess, including the color-octet channels. Figure 2 shows 1S0^8 is one of the dominant CO channels. Unless the authors specify which linear combination of f- and d-type functions enters each CO amplitude, the quoted asymmetries are not unambiguously predictions for 'the' gluon Sivers function. This could shift or even flip the sign of the few-percent effect. It does not destroy the paper as an estimate, but it should be fixed or at least discussed in a revision.\n\nAlso, the TMD-evolved sign change from positive to negative is driven by the u/d-based model choice for the GSF, not by a fitted gluon Sivers function. The authors call this exploratory, so it is a caveat rather than an error. The fixed kinematic choice (z=0.3, K_T=3 GeV) and the asserted independence from the LDME set are minor; the latter is asserted rather than demonstrated, but the conclusion is unlikely to change qualitatively.\n\nBottom line: a serious referee should engage with this. It is honest, transparent, and gives a concrete prediction that the EIC can eventually test. The f/d issue for the color-octet contributions needs attention, and the exploratory TMD-evolved numbers should be labeled as model studies. I would put this forward for peer review as a solid phenomenological contribution.","headline":"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.","tokens_in":14970,"tokens_out":3403,"would_cite":true,"duration_ms":35120,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["Sivers asymmetry","gluon Sivers function","transverse momentum dependent distributions","J/psi photoproduction","electron-ion collider","NRQCD","TMD evolution","back-to-back jets"],"falsifier":"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.","tokens_in":13895,"feed_emoji":"🎯","tokens_out":15506,"duration_ms":138494,"temperature":0.7,"pith_summary":"This paper argues that photoproduction of a nearly back-to-back $J/\\psi$-jet pair at the future Electron-Ion Collider (EIC) provides a clean, gluon-dominated window into the gluon Sivers function, the transverse-momentum-dependent distribution describing how gluons are distributed asymmetrically inside a transversely polarized proton. The authors calculate the $\\sin\\phi_q$-weighted single-spin asymmetry $A_N^{\\sin(\\phi_q)}$ for $ep^\\uparrow \\to J/\\psi+\\mathrm{jet}+X$ using the generalized parton model (GPM), which keeps the initial parton's intrinsic transverse momentum in the hard scattering, together with NRQCD for the quarkonium production rate. They include both color-singlet and color-octet mechanisms and find a few-percent asymmetry: about 3%, 1%, and 6% at $\\sqrt{s}=45$ GeV for three DGLAP-evolved gluon Sivers fits, and up to 8% in an exploratory TMD-evolved parametrization, while quark-initiated contributions are negligible. If this holds, the EIC could determine the gluon Sivers function in a kinematic region where the $J/\\psi$ does not have to be forward, complementing existing probes.","feed_headline":"EIC J/psi-jet pairs expose gluon Sivers asymmetry","feed_subtitle":"The predicted few-percent spin asymmetry is gluon-dominated, so the channel could pin down the gluon Sivers function.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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$."],"supporting_citations":[{"why":"Supplies the NRQCD matrix elements for the $\\gamma g\\to J/\\psi g$ subprocess that dominate the asymmetry.","marker":"[29]"},{"why":"Provides the DGLAP-evolved gluon Sivers parametrizations labeled SIDIS1 and SIDIS2 used in the numerical estimates.","marker":"[7]"},{"why":"Provides the re-fit gluon Sivers parametrization labeled SIDIS3 and the saturation choice $N_a=1$, $\\rho=2/3$.","marker":"[8]"},{"why":"Provides the TMD-evolution formalism and the u/d quark Sivers parametrizations used to construct the exploratory gluon parametrizations TMD-a and TMD-b.","marker":"[13]"},{"why":"Gives the equivalent-photon flux of the electron used for quasi-real photoproduction.","marker":"[54]"},{"why":"Defines the convention for the Sivers function and the azimuthal weight used in the asymmetry definition.","marker":"[56]"},{"why":"Supplies the color-singlet and color-octet long-distance matrix elements for $J/\\psi$ production.","marker":"[64]"},{"why":"Supplies the collinear parton distribution functions used as input to the Gaussian TMD and Sivers parametrizations.","marker":"[68]"}],"fun_headline_variants":["J/psi-jet photoproduction exposes gluon Sivers","Gluon Sivers from J/psi-jet at EIC","Sizable Sivers asymmetry in J/psi-jet channel","EIC J/psi-jet: gluon Sivers sign flip","Back-to-back J/psi-jet measures gluon Sivers"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["J/psi-jet photoproduction exposes gluon Sivers","Gluon Sivers from J/psi-jet at EIC","Sizable Sivers asymmetry in J/psi-jet channel","EIC J/psi-jet: gluon Sivers sign flip","Back-to-back J/psi-jet measures gluon Sivers"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000691,"raw_usage":{"total_tokens":3124,"prompt_tokens":935,"completion_tokens":2189,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":551,"completion_tokens_details":{"reasoning_tokens":2097}},"tokens_in":551,"tokens_out":2189,"duration_ms":17298,"temperature":1.0,"reasoning_tokens":2097,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:08:14.702058+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Measurements of inelastic J/psi and psi' photoproduction at HERA","cited_arxiv_id":"hep-ex/0211011","evidence_quote":"Gives the equivalent-photon flux of the electron used for quasi-real photoproduction."}],"review_version":1}