{"id":"6c4dfac4-b790-4b39-afca-5ae6c3108ec1","arxiv_id":"2505.05172","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Gluon spin-orbit coupling in a fast proton generates a cos(phi) two-particle azimuthal correlation, offering a direct experimental probe of the proton's double helicity parton distributions.","lead":"A new theory paper proposes that the observed 'directed flow' sideways asymmetry in proton-proton and proton-nucleus collisions can come from a quantum spin-orbit effect in gluons, providing a new window into the proton's internal structure. The proposed signal has a distinctive momentum and nuclear-size dependence, so it can be searched for with existing forward detectors at the LHC and RHIC.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The v1 signal (18) rests on the imported small-x relation (17); if its coefficient or kinematic range fails at the k and Delta values of Figs. 1-2, the central mechanism loses support.","rationale":"The reader's weakest-assumption analysis identifies Eq. (17) as the fragile imported input, and my stress-test reaches the same conclusion. The paper is transparent about modeling choices for F_Delta q Delta q and frames the predictions as upper-limit estimates, which is good practice; those uncertainties affect the magnitude of v1 but not the existence of the mechanism. The truly load-bearing link is the relation between the gluon spin-orbit correlation C_g and the dipole S-matrix at small x. If that relation's coefficient, sign-independent magnitude, or range of validity is different from what is imported, the predicted cosine correlation either changes size or vanishes. The paper's self-citation of [8,21] is not by itself a flaw, but the absence of an independent check or an explicit error estimate in the relevant k and Delta range leaves the central claim conditional. I therefore keep the reader's CONDITIONAL verdict unchanged. A focused re-derivation and numerical point check of Eq. (17) would settle whether the concern is real or whether the mechanism is as robust as claimed.","tokens_in":11811,"tokens_out":13720,"duration_ms":142649,"concrete_test":"Re-derive Eq. (17) from the operator definitions (14)-(15) and the small-x helicity GTMD formalism of Ref. [8], keeping all terms of order k^2/M^2 and Delta^2/M^2. Then numerically evaluate the ratio R = (x C_g/M^2) / (-N_c S/(8 pi^4 alpha_s)) for x' = 10^-4 to 10^-2, k = 0.5 to 3 GeV, and Delta = 0.1 to 1 GeV using the GBW and MVe dipole amplitudes of Figs. 1-2. If |R-1| exceeds 30% anywhere in this kinematic box, Eq. (18) and the v1 predictions (20)-(25) require revision; if |R-1| stays small, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central derivation collapses unless Eq. (17), x C_g(x,k,Delta)/M^2 ≈ -N_c/(8 pi^4 alpha_s) S(x,k,Delta), is accurate at the actual kinematics of the numerical predictions. This relation is imported from prior self-cited work [8,21] and is explicitly labeled 'a crucial use'. The paper only quotes the condition x << 1, but the application involves k ~ 0.5-3 GeV and Delta ~ 0.1-1 GeV, for which neglected corrections of order k^2/M^2, Delta^2/M^2, and subleading x behavior could be of order unity. Because C_g appears squared in Eq. (18), the sign of (17) is immaterial, but the magnitude and functional form are fully load-bearing: a 30% change in the coefficient changes v1 by 30%, and if the relation is simply not valid away from k=Delta=0, the predicted cos(phi1-phi2) correlation need not exist at all. The paper does not provide an independent derivation of (17), nor does it quantify its theoretical uncertainty in the relevant regime. The model dependence of F_Delta q Delta q is honestly acknowledged, but that is secondary: even with perfectly known DPDFs, the mechanism stands or falls on (17).","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a new mechanism for directed flow (a cos(phi1-phi2) two-particle azimuthal correlation) in forward proton-proton and proton-nucleus collisions. The mechanism is double parton scattering in which the square of the helicity-dependent quark-gluon amplitude (1) is combined with the small-x gluon spin-orbit relation (17). The resulting cross section (18) is proportional to the proton's double helicity parton distribution F_Delta q Delta q times the square of the dipole S-matrix, and the ratio of (18) to the unpolarized DPS cross section gives v1 ~ (1/|P|) sqrt(|F_Delta q Delta q|/F_qq) with an A^{-1/3} target dependence. Numerical estimates are presented for LHC and RHIC kinematics using the GBW and MVe dipole models, with the maximally negative double-helicity PDF ansatz (26).","tokens_in":12099,"tokens_out":17041,"duration_ms":163471,"significance":"If the central relation (17) is valid in the kinematics used, the paper identifies a genuinely new, experimentally testable source of directed flow in cold QCD collisions, with distinctive 1/|P| and A^{-1/3} scalings that are falsifiable with LHCb forward data. The derivation is transparent about its main approximations: the color-octet DPDF is neglected in (10), the double-gluon distribution is factorized in (9), and the Gaussian and factorized-Delta models are introduced in (22) and (23). The paper also correctly emphasizes that the observable can probe the poorly known double helicity DPDF, whose sign can in principle be read off from the sign of v1. The central weakness is the complete dependence of the numerical predictions on the imported small-x spin-orbit relation (17), whose quantitative validity at the k and Delta values of Figures 1 and 2 is not established in the manuscript.","major_comments":[{"comment":"The v1 signal is built on the imported relation x C_g(x,k,Delta)/M^2 ~ -N_c/(8 pi^4 alpha_s) S(x,k,Delta). The paper states only 'valid when x << 1', but the numerical applications involve x' ~ 10^-2 to 10^-4, k ~ 0.5-3 GeV and |Delta| ~ 0.1-1 GeV. Since C_g enters squared in Eq. (18), a moderate error in the coefficient shifts v1 by the same relative amount, and a breakdown of the relation away from k=Delta=0 would remove the predicted cos(phi1-phi2) correlation entirely. The authors should provide an independent derivation of (17), or at least a quantitative estimate of the O(k^2/M^2), O(Delta^2/M^2) and subleading-x corrections, with an explicit statement of the range of (k, Delta, x') over which the relation is accurate for the kinematics of Figures 1 and 2.","section":"II, Eqs. (17)-(18)"},{"comment":"The Gaussian integration as written does not reproduce the quoted coefficient. With F_Delta q Delta q proportional to exp(-Delta^2/(2 sigma^2)) and S(x',k,Delta) proportional to exp(-R_A^2 Delta^2/4), the Delta-integrals in Eq. (20) yield the inverse of beta = 1/(2 sigma^2) + R_A^2/2, not 2/(sigma^2 + R_A^2). The correct analytic expression for v1 is v1 ~ 1/(2|P|) sqrt(2 sigma^2/(1 + sigma^2 R_A^2)), which for sigma = 0.26 GeV and R_A = 1.2 fm (A = 1) gives approximately 0.098/|P|; the formula displayed in Eq. (25) gives a different value. The equations should be corrected, and the unit convention relating sigma (GeV) and R_A (fm) should be stated explicitly, because all numerical results in Figures 1 and 2 are generated from Eq. (24).","section":"III, Eqs. (24)-(25)"}],"minor_comments":[{"comment":"The plotted magnitudes in Figures 1 and 2 are upper limits obtained by saturating the positivity bound with r = -1. The text does mention this, but the abstract and introduction should more prominently state that the central quantitative claim is conditional on this model input, given that the literature quoted in the same paragraph reports both negative and positive values of r.","section":"III, Eq. (26)"},{"comment":"The ellipsis in Eq. (15) hides the subleading terms in the expansion of the helicity GTMD; a brief statement about the expected size of these terms would help the reader judge the validity of the subsequent substitution.","section":"II, Eq. (15)"},{"comment":"The factor of 1/2 for identical partons is stated after Eq. (6) but is not explicitly carried through Eq. (18) or the ratio (20); the authors should state whether both the numerator and denominator of the v1 ratio include the identical-particle symmetry factor, since the cancellation in the ratio is not automatic.","section":"II, Eq. (6) and Eq. (18)"},{"comment":"The comparison of v1 magnitudes between LHC and RHIC is made with different yref and Pref values; the text should note that these are different kinematic reference choices and that an apples-to-apples comparison across energies would require identical reference selections.","section":"III, Figures 1 and 2"},{"comment":"The sentence 'Fig. 2 show the results' should read 'Fig. 2 shows the results'; there is also a missing comma after 'from [23]' in the sentence introducing Eq. (26).","section":"III, text after Eq. (24)"}],"recommendation":"major_revision","confidential_remarks":"The paper's central claim is conditional on the validity of Eq. (17), which is imported from the authors' previous work. I recommend that the editor select a referee familiar with Refs. [8] and [21] who can independently assess the range of validity of that relation, since the numerical predictions of the present paper inherit its full quantitative content. The manuscript is within the scope of the journal and the idea is novel."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The central idea is new: in the hybrid factorization/DPS framework, squaring the helicity-dependent piece of the quark-gluon scattering amplitude and using the small-x relation (17) generates a cos(phi1-phi2) correlation in forward pp/pA collisions, proportional to the projectile's double helicity PDF F_Delta q Delta q. That gives a concrete, possibly measurable forward v1 with a 1/|P| falloff and A^{-1/3} dependence. The derivation is transparent about its approximations—octet DPDFs neglected, double gluon distribution saturated by single nucleon insertions, Gaussian DPDF models—and the authors are honest that the numerics are an upper-limit estimate using r = -1. The analytic estimate (25) is a useful guide. The mechanism is genuinely different from the odderon mechanism of ref. [21].\n\nThe main soft spot is the imported relation (17). It is load-bearing: it enters squared, so the magnitude of v1 scales with its coefficient, and if it fails at the k ~ 0.5-3 GeV, Delta ~ 0.1-1 GeV used in the plots, the predicted correlation could be much weaker or absent. The paper cites [8,21] but does not give an independent check or estimate of corrections in that kinematic regime. That is not a fatal defect, but it should be addressed in revision. The factorization (9) is also an approximation, and the numerical results depend on the Gaussian widths; the paper acknowledges this, though a systematic uncertainty band would strengthen it. On the other hand, the shape of v1 (1/|P|, A^{-1/3}) follows from the structure of the formulas and is not circular, and the model dependence of F_Delta q Delta q is openly stated.\n\nOverall, this is a solid, honest paper. The observable is promising, especially for LHCb forward data. I would engage with it: it deserves peer review, and the referee should push on Eq. (17) and on presenting the predictions as conditional on its validity. I'd take it to the reading group.","headline":"A plausible new mechanism for forward directed flow from gluon spin-orbit coupling and double helicity PDFs; the quantitative claim leans on one imported relation that the authors should be asked to back up.","tokens_in":12671,"tokens_out":2415,"would_cite":true,"duration_ms":26025,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.-t","13.85.-t"],"model":"deepseek-v4-flash","headline":"A new source of directed flow in forward proton-proton and proton-nucleus collisions arises from parton spin-orbit coupling and directly probes the proton's double helicity PDFs.","keywords":["directed flow","double parton scattering","double helicity PDF","gluon spin-orbit coupling","small-x QCD","color glass condensate","forward rapidity","azimuthal correlations"],"falsifier":"A decisive check is to measure forward rapidity $v_1$ in $pp$ collisions after standard nonflow subtraction and compare its $|P|$ dependence with the predicted $1/|P|$ behavior: observing $v_1\\propto|P|$ in the low-$|P|$ region, or a sharp peak near $|P|\\sim0.2$ GeV of the odderon type, would rule out this mechanism as the dominant source.","tokens_in":11534,"feed_emoji":"⚛️","tokens_out":11317,"duration_ms":104848,"temperature":0.7,"pith_summary":"The paper claims that forward-rapidity proton-proton and proton-nucleus collisions contain a previously overlooked source of the two-particle azimuthal correlation called directed flow ($v_1$). The source is the square of the helicity-dependent term in the elementary quark-gluon scattering amplitude, which becomes sizable because gluons at small $x$ have a strong spin-orbit correlation. If the derivation is right, measuring this correlation gives a direct handle on the proton's double helicity parton distribution $F_{\\Delta q\\Delta q}$ and on the small-$x$ gluon spin-orbit coupling, quantities that are otherwise hard to isolate. The predicted signal is a $v_1$ that falls roughly as $1/|P|$, which distinguishes it from hydrodynamic and odderon-based explanations.","feed_headline":"Spin-orbit coupling can create directed flow in pp","feed_subtitle":"Forward two-particle correlations would expose the proton's double-helicity PDFs.","key_machinery":"The machinery is the double parton scattering cross section built from the squared helicity-dependent quark-gluon amplitude, together with a single imported small-$x$ identity. The identity, Eq. (17), says that the gluon spin-orbit correlation $C_g(x,k,\\Delta)$, scaled by the nucleon mass squared, is proportional to the dipole $S$-matrix $S(x,k,\\Delta)$ with a known coefficient. The double helicity PDF $F_{\\Delta q\\Delta q}(x_1,x_2,\\Delta)$ plays the role of the projectile-side weight; after a Gaussian model for its $\\Delta$ dependence and a factorized ansatz for the target's $S$-matrix, the $\\Delta$ integral collapses and leaves a compact expression for $v_1$ that depends mainly on the ratio $|F_{\\Delta q\\Delta q}|/F_{qq}$ and not on the detailed form of $S$.","core_discovery":"The central discovery is that a $\\cos(\\phi_1-\\phi_2)$ two-particle correlation, interpreted as directed flow, appears in unpolarized $pp$ and $pA$ collisions at forward rapidity from the square of the helicity-flip term of the quark-gluon amplitude (Eq. 1). Combined with the small-$x$ relation $x C_g(x,k,\\Delta)/M^2 \\approx -N_c/(8\\pi^4 \\alpha_s)\\,S(x,k,\\Delta)$ (Eq. 17), the corresponding double parton scattering cross section (Eq. 18) becomes proportional to $F_{\\Delta q\\Delta q}(x_1,x_2,\\Delta) S(x'_1,k_1,-\\Delta) S(x'_2,k_2,\\Delta)$. The sign of the gluon spin-orbit coupling is lost because the distribution is squared, but its magnitude survives in the size of the correlation. The ratio of this contribution to the unpolarized double parton cross section yields $v_1 \\simeq (1/|P|)\\sqrt{|F_{\\Delta q\\Delta q}|/F_{qq}}$ in the no-fragmentation limit, so forward $v_1$ is a direct probe of the double helicity PDF.","pith_inferences":["If the same mechanism is evaluated with a quasi-real photon as the projectile, the correlation should appear in photoproduction-like final states, offering an independent test at future lepton-hadron colliders.","Because the sign of the gluon spin-orbit coupling is lost when the amplitude is squared, measuring this observable alone cannot determine whether gluon helicity and orbital angular momentum are aligned or anti-aligned; a polarized-beam version would be needed to recover that sign.","The neglected color-octet double PDF could contaminate the $S^2$ form with a color-quadrupole term; a model or lattice estimate of that octet distribution would show how safe the large-$N_c$ truncation is."],"forward_implications":["Forward two-particle correlation data at high-energy hadron colliders already exist, so the predicted $v_1$ can be searched for in reanalyzed samples without new instrumentation.","If confirmed, the sign of $v_1$ in a chosen channel such as $uu\\to\\pi^+\\pi^+$ fixes the sign of $F_{\\Delta u\\Delta u}$ and helps settle conflicting model predictions.","The predicted $1/|P|$ falloff and the $A^{-1/3}$ nuclear-size dependence separate this mechanism from hydrodynamic ($v_1\\propto|P|$) and odderon (low-$|P|$ peak) interpretations.","At the partonic level, $v_1$ is nearly independent of the dipole $S$-matrix; the only sizable model dependence enters through fragmentation functions.","In the maximally negative scenario $F_{\\Delta u\\Delta u}=-F_{uu}$, the numerical estimates put $|v_1|$ at the percent level for $|P|\\sim1$ GeV at both collider energies considered, which is in the measurable range."],"supporting_citations":[{"why":"Supplies the small-x relation $x C_g/M^2 \\approx -N_c/(8\\pi^4\\alpha_s) S$ (Eq. 17) that converts the gluon spin-orbit correlation into the dipole S-matrix.","marker":"[8]"},{"why":"Co-cited for Eq. (17) and considered as the odderon mechanism for the same $\\cos(\\Delta\\phi)$ correlation that this paper contrasts with its own source.","marker":"[21]"},{"why":"Provides the definitions of double parton distributions $F_{qq}$ and $F_{\\Delta q\\Delta q}$ and the color-octet variant used in the derivation.","marker":"[24]"},{"why":"Supplies the Gaussian width $\\sigma=0.26$ GeV and one of the model estimates for the ratio $F_{\\Delta u\\Delta u}/F_{uu}$ used in the numerics.","marker":"[25]"},{"why":"Provides the hybrid factorization formula for single-particle production that the double parton scattering derivation extends.","marker":"[29]"},{"why":"Gives the double parton scattering cross section with two dipole S-matrices that Eq. (6) builds on.","marker":"[31]"},{"why":"Derives the azimuthal-correlation machinery for double parton scattering with $S(r,b)$ anisotropies, which the paper adapts to the spin-orbit case.","marker":"[33]"},{"why":"Provides the small-x relation between the gluon generalized TMD and the dipole S-matrix used in Eq. (9) and the final cross section.","marker":"[34]"},{"why":"Supplies the dipole model used to compute the numerical $v_1$ curves in the figures.","marker":"[48]"},{"why":"Supplies the alternative dipole model used for the second set of numerical predictions.","marker":"[49]"}],"fun_headline_variants":["Spin-orbit coupling seeds directed flow in pp and pA","Directed flow from gluon spin-orbit coupling","New mechanism for directed flow: spin-orbit coupling","Spin-orbit coupling yields forward directed flow","Probing double-helicity PDFs with directed flow"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the small-$x$ relation (17) imported from earlier work is correct in sign, magnitude, and kinematic range at $x'\\sim10^{-2}$ to $10^{-4}$ and $|\\Delta|\\sim0.1$ to $1$ GeV, together with the factorization of the target's double-gluon distribution into single distributions (Eq. 9).","fun_headline_variants_meta":{"raw":{"variants":["Spin-orbit coupling seeds directed flow in pp and pA","Directed flow from gluon spin-orbit coupling","New mechanism for directed flow: spin-orbit coupling","Spin-orbit coupling yields forward directed flow","Probing double-helicity PDFs with directed flow"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00044,"raw_usage":{"total_tokens":2191,"prompt_tokens":865,"completion_tokens":1326,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":481,"completion_tokens_details":{"reasoning_tokens":1250}},"tokens_in":481,"tokens_out":1326,"duration_ms":9354,"temperature":1.0,"reasoning_tokens":1250,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:11:09.857837+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive check is to measure forward rapidity $v_1$ in $pp$ collisions after standard nonflow subtraction and compare its $|P|$ dependence with the predicted $1/|P|$ behavior: observing $v_1\\propto|P|$ in the low-$|P|$ region, or a sharp peak near $|P|\\sim0.2$ GeV of the odderon type, would rule out this mechanism as the dominant source.","supporting_citations":[{"cited_title":"Directed flow from parton spin-orbit coupling in $pp$ and $pA$ collisions","cited_arxiv_id":"2505.05172","evidence_quote":"Co-cited for Eq. (17) and considered as the odderon mechanism for the same $\\cos(\\Delta\\phi)$ correlation that this paper contrasts with its own source."},{"cited_title":"Angular correlations in the double Drell-Yan process","cited_arxiv_id":"1210.5434","evidence_quote":"Provides the hybrid factorization formula for single-particle production that the double parton scattering derivation extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the small-x relation between the gluon generalized TMD and the dipole S-matrix used in Eq. (9) and the final cross section."}],"review_version":1}