{"id":"203a26a8-1625-4af6-97a0-5153161af521","arxiv_id":"1908.07312","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A combined analysis of RHIC pion and D-meson spin asymmetries yields first rough constraints on the two independent gluon Sivers functions in the CGI-GPM.","lead":"This paper uses RHIC data on single-spin asymmetries in pion and D-meson production to estimate the strength of two gluon Sivers functions, which describe how gluons are arranged inside a spinning proton. The result is a preliminary constraint on the process-dependent gluon Sivers effect within the color gauge invariant generalized parton model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The numeric constraints in Eq. (3.3) depend on an assumed x-shape and are not validated against both data sets, so the central claim is only as strong as that Ansatz.","rationale":"The reader's CONDITIONAL verdict is well matched to the actual status of the paper. The central claim is explicitly hedged as 'preliminary constraints' and depends on the CGI-GPM framework already developed in Refs. [5,6]. The load-bearing weakness is exactly the assumed x-dependence in Eqs. (3.1)-(3.2), combined with the absence of a proper fit or of an overlay of the final endpoint scenarios on both datasets. This is not an internal inconsistency: the model is well defined and the qualitative conclusion that the f-type gluon Sivers function is strongly reduced is plausible. However, the numeric ranges in Eq. (3.3) should not be cited as robust constraints until the shape dependence is tested and the endpoint scenarios are shown to reproduce the data. Since the reader already recommended a conditional acceptance with these caveats, the stress-test does not move the verdict.","tokens_in":6221,"tokens_out":4807,"duration_ms":54856,"concrete_test":"Recompute the CGI-GPM asymmetries for the two endpoint pairs of Eq. (3.3), namely (Ng(f),Ng(d))=(0.05,-0.15) and (-0.01,0.15), and overlay both the PHENIX pi0 pT spectrum (Ref. [8]) and the D->mu+- xF data (Ref. [9]) using the same code as in Figs. 1-2. If either endpoint curve falls outside the data uncertainties, the quoted ranges are not supported. As a secondary check, repeat the extraction with a modified x-Ansatz, e.g., comparing uniform Ng(x) with alpha=0.6, beta=0.6, and verify whether the allowed Ng(f) range shifts by more than the stated width; if so, the constraints are shape-dominated rather than data-dominated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The quantitative central claim, Eq. (3.3), is obtained by fixing the x-dependence of both gluon Sivers functions to the Ansatz (3.1)-(3.2), with rho=2/3 for pions and Ng(f,d)(x) taken constant for the D-meson analysis. The available PHENIX data cover only a few pT or xF bins, so they constrain an overall normalization for each channel, not the shape of the x-dependence. If the true x-dependence of f⊥g(f) or f⊥g(d) differs from the adopted form (for example, if alpha or beta are not the assumed values, or if the f- and d-type functions have different shapes), the quoted ranges +0.05>=Ng(f)>=-0.01 and -0.15<=Ng(d)<=+0.15 would change. The paper also does not display the final endpoint scenarios against both datasets: Fig. 1 right panel uses Ng(d)=-1, which is outside the final d-type range, while Fig. 2 shows only d-type contributions with Ng(f)=0.1. The statement that a 'very good description' of both pi0 and muon data is possible for the endpoint pairs in Eq. (3.3) is therefore asserted rather than demonstrated. The central claim of 'preliminary constraints' is valid only as an illustrative statement within the assumed parametrization, not as a measured property of the gluon Sivers functions.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This proceedings paper analyzes PHENIX single-spin asymmetries for p↑p → π0X at midrapidity and for muons from D-meson decays, as functions of pT and xF. In the color-gauge-invariant generalized parton model (CGI-GPM), two independent gluon Sivers functions, f-type and d-type, enter with modified hard functions. The authors adopt a factorized Gaussian ansatz for the Sivers x- and k⊥-dependence, saturate positivity for orientation, and then scan the normalizations Ng(f) and Ng(d). They conclude that the f-type gluon Sivers function is strongly reduced (+0.05 ≥ Ng(f) ≥ −0.01), while the d-type is small (−0.15 ≤ Ng(d) ≤ +0.15), and they compare these results with the simpler GPM.","tokens_in":6544,"tokens_out":6300,"duration_ms":62017,"significance":"The qualitative finding, if correct, is useful: it suggests that initial- and final-state interactions suppress the f-type gluon Sivers function relative to the GPM extraction, and it informs predictions for J/ψ and direct-photon processes at RHIC. The paper is transparent about its model assumptions, applies positivity constraints, and displays the relevant hard functions; the explicit GPM comparison is a valuable cross-check. The numerical ranges in Eq. (3.3) are, however, illustrative constraints within a fixed ansatz rather than a measured extraction, a limitation the authors acknowledge in the text.","major_comments":[{"comment":"The endpoint pairs in Eq. (3.3) are never actually shown against both data sets. The right panel of Fig. 1 uses Ng(d) = −1, which lies outside the final d-type range, and Fig. 2 shows only the d-type contribution. The statement that a very good description of both the μ± and π0 data is possible for the pairs in Eq. (3.3) is therefore asserted rather than demonstrated. Please add plots, or a table of χ² values, in which the final endpoint scenarios are overlaid on the PHENIX π0 and μ± data.","section":"Section 3, Eq. (3.3) and Figs. 1-2"},{"comment":"The constraints refer only to the normalizations Ng(f) and Ng(d) within the assumed x- and k⊥-shapes, with fixed α, β, and ρ, and with Ng(x) taken constant for the D-meson analysis. The PHENIX data cover only a few pT or xF bins, and the paper itself states that the number and precision of the available data points do not allow for a true fit. If the true x-dependence differs from the ansatz, the quoted ranges change. The abstract's phrase constraints on the two independent gluon Sivers functions should be explicitly qualified as constraints within the adopted parametrization.","section":"Section 3, Eqs. (3.1)-(3.3)"},{"comment":"No quantitative goodness-of-fit measure is given. The only uncertainty estimate is the ±20% band on Ng(f) in Fig. 1 at fixed endpoints, and no error bars are attached to the ranges in Eq. (3.3). It is therefore unclear whether Eq. (3.3) delimits an allowed region or merely a pair of manually chosen points. A χ² or profile-χ² scan over Ng(f) and Ng(d) would substantiate the word constraints.","section":"Section 3, Eq. (3.3)"}],"minor_comments":[{"comment":"The sentence with Ng(f) = 0.05 and Ng(f) = 0.15 appears to contain a typo; the second normalization should presumably be Ng(d) = 0.15 for the d-type curve.","section":"Section 3, text after Eq. (3.4)"},{"comment":"The abstract says for the first time, while the text is more cautious and calls the analysis a first attempt; the abstract should include the model-dependence qualification.","section":"Abstract and Section 3"},{"comment":"The figure captions do not fully specify all curves. Fig. 1 should state the exact values of Ng(f) and Ng(d) in the right panel, and Fig. 2 should clarify the mapping of the thin solid and dash-dotted lines to Ng(d) = +1 and Ng(d) = −1.","section":"Figs. 1 and 2"},{"comment":"The notation such as N( f )g is visually ambiguous; consistent use of Ng(f) would improve readability.","section":"Notation"}],"recommendation":"major_revision","confidential_remarks":"This is a conference proceedings, and the authors are candid that no true fit is possible. My main reason for major_revision is the gap between the abstract's constraints claim and what is actually displayed: the final endpoint scenarios are not compared with both data sets. If the editors accept a lower threshold for proceedings, a minor revision with the missing plots and qualifications might suffice; under the usual journal standard, the requested additions are necessary before the quantitative claim is supported."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a compact proceedings contribution from a group that has done the heavy lifting elsewhere. The genuinely new step is combining existing PHENIX pi0 and D-meson data within the CGI-GPM framework to get simultaneous ranges for the f-type and d-type gluon Sivers normalizations, Eq. (3.3). That is a legitimate extension of their own program, and the qualitative conclusion—f-type strongly suppressed, d-type small—looks consistent with the plotted curves. The GPM comparison provides a useful cross-check, and the authors are candid that the data 'do not allow for a true fit.' I give them credit for that honesty and for stating the preliminary nature of the result in the abstract.\n\nThe soft spots are real, though. The central numbers come from manual parameter scans, with no chi-square or error bars. More importantly, the quoted ranges depend on the assumed x-shape in Eqs. (3.1)–(3.2), including fixed exponents and rho = 2/3 for pions. If the true x-dependence differs, the ranges in (3.3) are not constraints on the actual functions. The stress-test note is right that the endpoint scenarios are not displayed against both data sets: Fig. 1 uses Ng(d) = -1, outside the final d-type range, and Fig. 2 shows only the d-type contribution with Ng(f) fixed at 0.1. So the claim that a 'very good description' of both pi0 and muon data is possible is asserted rather than demonstrated.\n\nThat said, these are soft spots in a proceedings-length paper, not fatal flaws. The framework is established, the external data are real, and the authors openly flag the limitations. If this were submitted as a proper research letter, I would send it to a referee—the topic is important for TMD spin physics and the analysis is a plausible first step. But the referee should ask for either a genuine fit or a prominent caveat that the numeric ranges are conditional on the adopted parametrization. I would not myself cite the numbers in (3.3) as extracted values without checking whether a later, more complete analysis has superseded them.","headline":"A short, honest proceedings paper that takes a small but real step toward constraining the two gluon Sivers functions; treat the numeric ranges as Ansatz-dependent illustrations, not measured constraints.","tokens_in":7108,"tokens_out":1651,"would_cite":false,"duration_ms":19884,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"RHIC pion and D-meson data put the first constraints on the two gluon Sivers functions","keywords":["gluon Sivers function","transverse single-spin asymmetry","transverse momentum dependent distributions","color gauge invariant generalized parton model","process dependence","RHIC","D meson production","neutral pion production"],"falsifier":"A measurement of single-spin asymmetries in $p^\\uparrow p \\to J/\\psi\\, X$ at RHIC, where the gluon-gluon fusion channel dominates and the $f$- and $d$-type modified hard functions weight differently than in pion and $D$-meson production, would settle the central claim: if the data demand an $f$-type normalization outside $+0.05$ to $-0.01$, the combined constraint from pion and $D$-meson data is contradicted.","tokens_in":5971,"feed_emoji":"🎯","tokens_out":8495,"duration_ms":81829,"temperature":0.7,"pith_summary":"The paper tries to establish that existing RHIC measurements of transverse single-spin asymmetries in neutral pion and D-meson production can already constrain the gluon Sivers function, the transverse-momentum distribution of gluons inside a transversely polarized proton. It works in a version of the parton model that includes initial- and final-state interactions, where two independent gluon Sivers functions enter, an f-type and a d-type. The combined data give, for the first time, preliminary ranges for their normalizations: $N_g^{(f)}$ between $+0.05$ and $-0.01$, and $N_g^{(d)}$ between $-0.15$ and $+0.15$, with the f-type strongly suppressed and a slight preference for a positive d-type. The paper also compares with the simpler generalized parton model and notes that current data cannot yet discriminate between the two frameworks.","feed_headline":"First RHIC constraints on two gluon Sivers functions","feed_subtitle":"Pion and D-meson spin asymmetries put the f-type coupling near zero and the d-type within ±0.15.","key_machinery":"The load-bearing object is the pair of modified hard-scattering functions $H^{\\mathrm{Inc}(f)}$ and $H^{\\mathrm{Inc}(d)}$, which replace the universal hard functions when an unpolarized gluon from the polarized proton undergoes initial- or final-state interactions with the remnant. They encode the process dependence: the $f$-type combination enters with different signs in channels such as $gq\\to gq$ versus $g\\bar q\\to g\\bar q$, and the $d$-type combination vanishes for $gg\\to gg$ and cancels between quark and antiquark channels, which is why the $d$-type contribution is naturally small. The second moving part is the Ansatz of Eqs. (3.1)-(3.2), which fixes the shape of both Sivers functions and leaves only the normalizations $N_g^{(f)}$ and $N_g^{(d)}$ free.","core_discovery":"The paper's central claim is that, at the level of the adopted parametrization, the RHIC pion and D-meson data put the first quantitative limits on the two independent gluon Sivers functions of the color gauge invariant generalized parton model. In that model, initial- and final-state interactions break the universality of the gluon Sivers function and generate two color-singlet combinations: an $f$-type, even under charge conjugation, and a $d$-type, odd under charge conjugation. The neutral-pion data alone would allow a large negative $d$-type contribution, but the muon-from-$D$ data limit $N_g^{(d)}$ to $|N_g^{(d)}|\\le 0.15$ with a slight preference for positive values; combining both data sets then forces $N_g^{(f)}$ into the narrow range $+0.05\\ge N_g^{(f)}\\ge -0.01$. The paper stresses that this is not a true fit and that the generalized parton model description still agrees with the data, so the main result is a first, model-dependent constraint rather than a measurement.","pith_inferences":["A natural next step would be to promote the fixed exponents in the Ansatz to free parameters and refit the combined data; this would test whether the narrow $N_g^{(f)}$ range is an artifact of the assumed $x$-shape or a stable feature.","Because the $d$-type function is odd under charge conjugation, processes that select the $gg$-initiated subprocess with no antiquark compensation, such as quarkonium production, should provide a sharper handle on the $d$-type function than the channels analyzed here.","The same modified hard-function machinery could be applied to future electron-ion collider data, where TMD measurements access gluon Sivers effects with different kinematic coverage and could break the remaining degeneracy between the $x$-shape and the normalization.","The quoted ranges come from a specific choice of $\\rho=2/3$; varying this transverse-momentum parameter would move the extracted normalizations, so the ranges should be read as conditional on that choice."],"forward_implications":["If the extracted ranges hold, the $f$-type gluon Sivers function is close to zero at RHIC kinematics, so future analyses of gluon-initiated single-spin asymmetries should use a small $f$-type normalization rather than the large values allowed by positivity alone.","The $d$-type function, although small, is not forced to vanish and is slightly preferred positive; processes sensitive to the $d$-type color combination should show small but nonzero asymmetries with a definite sign pattern.","Because the two functions enter different partonic channels with different signs, the gluon Sivers effect is genuinely process-dependent; no single universal function can describe both pion and open-heavy-flavor data.","The comparison with the generalized parton model shows that the simpler framework still fits the data, so additional channels are needed before process dependence itself is established."],"supporting_citations":[{"why":"Establishes the generalized parton model framework used as the baseline for computing single-spin asymmetries.","marker":"[1]"},{"why":"Defines the Sivers function and its azimuthal distribution inside a transversely polarized proton.","marker":"[2]"},{"why":"Derives the process-dependent quark Sivers function from initial- and final-state interactions, the basis of the CGI-GPM approach.","marker":"[3]"},{"why":"Extends the CGI-GPM to the gluon sector and provides the two independent gluon Sivers functions and their modified hard functions.","marker":"[5]"},{"why":"Provides the earlier GPM analysis of the gluon Sivers function and the parametrization used for the generalized parton model comparison.","marker":"[6]"},{"why":"Supplies the neutral pion single-spin asymmetry data at RHIC that drive the main constraint.","marker":"[8]"},{"why":"Supplies the muon-from-D decay single-spin asymmetry data used to constrain the d-type gluon Sivers function.","marker":"[9]"}],"fun_headline_variants":["First RHIC bounds on both gluon Sivers functions","RHIC pion and D data split gluon Sivers into two types","Gluon Sivers pair: RHIC sets first f-type and d-type limits","Initial and final state interactions yield two gluon Sivers at RHIC","Pion and D meson asymmetries constrain both gluon Sivers at RHIC"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper assumes the $x$- and transverse-momentum dependence of both gluon Sivers functions is fixed by the Ansatz of Eqs. (3.1)-(3.2), so the data constrain only the overall normalizations $N_g^{(f)}$ and $N_g^{(d)}$; if the true shape differs, the quoted ranges do not apply to the actual functions.","fun_headline_variants_meta":{"raw":{"variants":["First RHIC bounds on both gluon Sivers functions","RHIC pion and D data split gluon Sivers into two types","Gluon Sivers pair: RHIC sets first f-type and d-type limits","Initial and final state interactions yield two gluon Sivers at RHIC","Pion and D meson asymmetries constrain both gluon Sivers at RHIC"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000402,"raw_usage":{"total_tokens":2055,"prompt_tokens":865,"completion_tokens":1190,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":481,"completion_tokens_details":{"reasoning_tokens":1091}},"tokens_in":481,"tokens_out":1190,"duration_ms":11376,"temperature":1.0,"reasoning_tokens":1091,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:20:26.314866+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A measurement of single-spin asymmetries in $p^\\uparrow p \\to J/\\psi\\, X$ at RHIC, where the gluon-gluon fusion channel dominates and the $f$- and $d$-type modified hard functions weight differently than in pion and $D$-meson production, would settle the central claim: if the data demand an $f$-type normalization outside $+0.05$ to $-0.01$, the combined constraint from pion and $D$-meson data is contradicted.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Sivers function and its azimuthal distribution inside a transversely polarized proton."},{"cited_title":"Process dependent Sivers function and implications for single spin asymmetry in inclusive hadron production","cited_arxiv_id":"1009.1936","evidence_quote":"Derives the process-dependent quark Sivers function from initial- and final-state interactions, the basis of the CGI-GPM approach."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the neutral pion single-spin asymmetry data at RHIC that drive the main constraint."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the muon-from-D decay single-spin asymmetry data used to constrain the d-type gluon Sivers function."}],"review_version":1}