{"id":"d3ab86d7-c5bd-44d6-8a80-6fc50f2a240d","arxiv_id":"2608.10955","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Gluonic nucleon energy correlators in the CGC reduce at eikonal accuracy to the unpolarized and linearly polarized gluon components, both given by the adjoint dipole S-matrix, and the linearly polarized one drives a saturation-sensitive cos 2phi asymmetry in DIS.","lead":"This paper derives the small-x gluonic nucleon energy correlators and fracture functions inside the Color Glass Condensate, finding that only two gluon components survive at eikonal order and that both are fixed by the adjoint dipole S-matrix. It uses the linearly polarized component to predict a cos 2phi azimuthal asymmetry in the target-fragmentation energy pattern of deep inelastic scattering, a proposed signature of gluon saturation at the future Electron-Ion Collider.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The predicted cos2phi suppression rests on an untested endpoint replacement of the z-convolution in Eq. (4.12); if the approximation is inaccurate, the central saturation signature is not established.","rationale":"I read the central claim as two-tier: the analytic CGC statement that only f_1^g and h_1^{t,g} survive and are fixed by the adjoint dipole S-matrix, and the numerical claim that their ratio generates a sizable, saturation-sensitive cos2phi asymmetry. The analytic tier is supported by nontrivial cross-checks: the dilute limit matches the collinear fracture functions in Appendix A, the large-theta limit of the NECs matches Ref. [108], and the jet fracture function agrees with Ref. [117] in the overlapping x<<1, p_g_perp^2 << k_g_perp^2 << Q^2 limit. The self-identified limitation in Sec. 5 that the full CGC-to-collinear TFR matching is not performed is honest but means the phenomenological tier relies on factorization assumptions. Within that tier, I agree with the Reader that Eq. (4.13) is the weakest link: it is a numerical shortcut with no error estimate, and the figures in which the paper's main result appears all use it. I do not elevate the hand-chosen nuclear Q_s0A range to the primary concern because it is an input uncertainty that could be propagated; the endpoint replacement is a structural approximation inside the predicted ratio and cannot be cured by broadening error bars. A secondary issue worth checking in the same numerical pass is that the exact normalization of Eq. (4.4) contains A Sigma_T + B Sigma_L, while Eq. (4.6) drops Sigma_L; if Sigma_L is numerically non-negligible at x_B=5e-3, the absolute asymmetry would also shift. Both issues are resolvable by one direct evaluation of the unsimplified formulas.","tokens_in":32059,"tokens_out":27421,"duration_ms":280121,"concrete_test":"Implement Eqs. (3.34), (4.8), and (4.12) with the same rcBK inputs and kinematics (x_B=5e-3, Q^2=25 GeV^2, sqrt(s)=89 GeV): compute h_1^{t,g}(z,theta) on a fine grid z in [x_B,1], perform the convolution integral_{x_B}^1 (dz/z)(x_B/z)^3 z h_1^{t,g}(z,theta), and compare it with (1/3) x_B h_1^{t,g}(x_B,theta) for theta in [0.1,0.9]. Repeat for proton and Au with Q_s0A^2 = 3, 4, 5 Q_s0^2. If the ratio of the exact to the approximate result deviates from 1 by more than about 30% anywhere in this range, the endpoint replacement is not reliable and the asymmetry figures need to be recomputed or shown with an error band.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The numerical asymmetry prediction depends on the step from Eq. (4.12) to Eq. (4.13), where the z-integral over T_F (x_B/z)^3 z h_1^{t,g}(z,theta) is replaced by (T_F/3) x_B h_1^{t,g}(x_B,theta). This is not a perturbative power-counting result; it is a small-x ansatz. In the saturation window theta Q <= Q_s, h_1^{t,g} is computed from Eq. (3.34) with F_{x_g}, the cutoff theta(x_g-xi), and epsilon_f^2=(x/xi)p_perp^2, so its z-dependence is not the simple 1/z form that would make the replacement harmless. No estimate of the truncation error is given, yet Eq. (4.14) and all of Fig. 12 inherit it. If the true convolution differs from the endpoint value by an O(1) factor, the predicted asymmetry magnitude and its theta-dependence, and therefore the claimed eA/ep suppression, are not established. The formal CGC derivation may still be correct; this is a numerical-phenomenology concern, but it is the load-bearing step for the paper's main observable.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes leading-twist gluonic nucleon energy correlators (NECs) and fracture functions in the Color Glass Condensate effective theory at eikonal accuracy, starting from operator definitions. It finds that for an unpolarized target only the unpolarized gluon NEC f_1^g and the linearly polarized gluon NEC h_1^{t,g} survive, and both are expressed in terms of the adjoint dipole S-matrix (Eqs. (3.30) and (3.34)). The dilute limit is matched to collinear splitting results (Sec. 3.2.3 and Appendix A), and the unpolarized gluon jet fracture function is matched to Ref. [117] (Sec. 3.2.4). The paper then proposes a cos 2φ azimuthal asymmetry of the TFR DIS energy pattern, given by Eq. (4.14), and presents numerical predictions for EIC kinematics using rcBK dipole evolution. The predicted asymmetry is negative, grows with θ, and shows a sizable nuclear suppression in e+A relative to e+p collisions.","tokens_in":32403,"tokens_out":11036,"duration_ms":107037,"significance":"If the results hold, this is a substantial step in extending CGC-based TFR physics to the gluon sector. The analytic core is credible: the derivation is based on operator definitions, the reduction to two unpolarized-target components is physically sensible, and the cross-checks in Secs. 3.2.3, 3.2.4, and Appendix A provide nontrivial support. The proposed cos 2φ asymmetry is well motivated and, unlike TMD-based dijet asymmetries, is formulated in collinear factorization, so it is not affected by Sudakov suppression. The analytic matrix elements contain no adjustable parameters; the numerical inputs are the rcBK dipole parameters from a HERA fit. The main weakness is that the quantitative predictions in Sec. 4 rest on an untested endpoint approximation for the z-convolution, which directly affects the magnitude, θ-dependence, and nuclear suppression of the headline observable.","major_comments":[{"comment":"The step from Eq. (4.12) to Eq. (4.13) replaces the z-convolution by the endpoint value x_B h_1^{t,g}(x_B,θ). This is presented as a small-x approximation, but no derivation or error estimate is given. The function h_1^{t,g}(z,θ) computed from Eq. (3.34) has a nontrivial z-dependence through ε_f^2=(x/ξ)p_⊥^2 and through the dipole S_{x_g}, and the kernel x_B^3/z^4 in the convolution does not automatically make the endpoint replacement exact. Since Eq. (4.14) and all of Fig. 12 inherit this approximation, the predicted magnitude of ⟨cos 2φ⟩ and its nuclear suppression R_eA are not quantitatively established until the full convolution is evaluated or a controlled estimate of the truncation error is provided.","section":"Sec. 4.2-4.3, Eqs. (4.12)-(4.14)"},{"comment":"The x_g prescription in the convolution is underspecified. The text states that the authors follow Ref. [103] and choose x_g=x_B, but Eq. (4.12) requires h_1^{t,g}(z,θ) for z>x_B, and Eq. (3.34) depends on x_g through θ(x_g-ξ) and through S_{x_g}. If one literally uses x_g=x_B for all z, the kinematic constraint in Eq. (3.34) becomes z-independent, which is inconsistent with the dilute-limit identification x_g=x/z used in Eq. (3.38). The authors should specify how x_g is assigned for z>x_B and show the sensitivity of the numerical results to this choice.","section":"Sec. 4.3, Eqs. (3.34), (3.38), and (4.12)"}],"minor_comments":[{"comment":"The Taylor expansion in Eq. (3.35) contains corrupted and unreadable symbols in the derivative terms; the expression needs to be cleaned up.","section":"Eq. (3.35)"},{"comment":"There appears to be an unbalanced bracket in the curly-brace expression around \"p_⊥·k_g⊥)]][ε_f^2+...\"; this should be corrected.","section":"Eq. (3.34)"},{"comment":"The integration measure in Eq. (2.8) is garbled: \"dP+ h⊥d2Ph⊥\" should presumably be dP_h^+ d^2P_h⊥, consistent with the other formulas.","section":"Eq. (2.8)"},{"comment":"The section title contains a typo: \"F racture function\" should be \"Fracture function\".","section":"Sec. 3.2.4 title"},{"comment":"The paper correctly notes in Sec. 5 that a direct derivation of the TFR factorization from the full CGC cross section remains to be done; it would be helpful to state this caveat also in Sec. 4.2, where the factorization formula is used for the saturation-sensitive region.","section":"Sec. 5"}],"recommendation":"major_revision","confidential_remarks":"The formal CGC derivation appears sound and well cross-checked; my main concern is the untested numerical approximation in Eq. (4.13), which is load-bearing for the paper's central phenomenological claim. The requested revision is straightforward: evaluate the full convolution or provide a controlled comparison, and specify the x_g prescription for z>x_B. I do not see a basis for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The analytic core is the real contribution. This is the first systematic calculation of gluonic NECs and fracture functions in the CGC at eikonal accuracy, starting from operator definitions. The reduction of the eight leading-twist gluon components to just two for an unpolarized target — the unpolarized f_1^g and the linearly polarized h_1^{t,g}, both fixed by the adjoint dipole S-matrix — is new and clean. The cross-checks are the right ones: the dilute/collinear limit reproduces the known 1/theta^2 behavior and splitting kernels, and Sec. 3.2.4 shows agreement with the unpolarized gluon jet fracture function of Caucal–Salazar in the overlapping limit. That gives me confidence the formalism is correct.\n\nThe observable is attractive: a Sudakov-free, single-energy-flow cos 2phi asymmetry in the TFR, governed by the ratio h_1^{t,g}/f_1^q. Since the gluon NEC sits on the adjoint dipole and the quark NEC on the fundamental, the ratio should be saturation-sensitive. Good idea.\n\nThe soft spot is exactly where the stress-test note lands. Eq. (4.13) replaces the z-convolution in Eq. (4.12) with the endpoint value, claiming small-x dominance, but no error estimate is given. In the saturation window the z-dependence of h_1^{t,g} is not literally 1/z, so the replacement could be off by an O(1) factor, shifting the magnitude and possibly the theta-dependence of the predicted asymmetry and its eA/ep suppression. The qualitative direction of the suppression is likely robust, but the paper should say so and back it with a check — for example, comparing the full convolution against the endpoint term for a model z-dependence. The hand-picked range 3–5 Q_s0^2 for the nuclear saturation scale is also a choice; Fig. 12 shows a band, but it is not propagated as an uncertainty.\n\nThese are numerical-phenomenology issues, not flaws in the derivation. The analytic results stand. I would send it to a serious referee and ask for the convolution check and some error estimate on the small-x approximation. With those added, this becomes a strong paper for the EIC small-x program.","headline":"A solid analytic CGC calculation of gluonic NECs and fracture functions with a genuinely new observable, but the numerics rest on an untested small-x convolution approximation.","tokens_in":32924,"tokens_out":3556,"would_cite":true,"duration_ms":31842,"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":"This paper shows that gluonic energy correlators reduce to two dipole-determined components and that their cos 2φ asymmetry is a saturation-sensitive observable at the EIC.","keywords":["nucleon energy correlators","fracture functions","Color Glass Condensate","small-x gluon saturation","linearly polarized gluons","target fragmentation region","cos 2φ azimuthal asymmetry","Electron-Ion Collider"],"falsifier":"A direct evaluation of the $z$-convolution in Eq. (4.12) using the computed $h_1^{t,g}(z,\\theta)$ from Eq. (3.34), without the small-$x$ substitution, would show whether the predicted $\\cos 2\\phi$ asymmetry and its nuclear suppression survive; alternatively, at the EIC, measuring $\\langle\\cos 2\\phi\\rangle$ in $e{+}p$ and $e{+}Au$ at $x_B = 5\\times10^{-3}$, $Q^2=25\\,\\mathrm{GeV}^2$ and finding $R_{eA} \\ge 1$ in the small-$\\theta$ region would contradict the central prediction.","tokens_in":31878,"feed_emoji":"⚛️","tokens_out":11714,"duration_ms":93530,"temperature":0.7,"pith_summary":"The paper derives gluonic nucleon energy correlators (NECs) and fracture functions from their operator definitions within the Color Glass Condensate effective theory, a high-energy description of dense gluonic matter, at eikonal accuracy. It finds that for an unpolarized target only two gluon components survive at this order—the unpolarized $f_1^g$ and the linearly polarized $h_1^{t,g}$—and that both are completely fixed by the adjoint dipole $S$-matrix. It then shows that $h_1^{t,g}$ produces a $\\cos 2\\phi$ azimuthal modulation of the DIS energy pattern in the target fragmentation region, with magnitude set by the ratio of $h_1^{t,g}$ to the unpolarized quark NEC $f_1^q$. Because this ratio compares an adjoint-dipole quantity to a fundamental-dipole quantity, the paper argues the asymmetry is especially sensitive to the saturation scale, and it predicts a sizable nuclear suppression in $e{+}A$ relative to $e{+}p$ collisions at the future Electron-Ion Collider.","feed_headline":"Nuclear suppression of a gluon asymmetry exposes saturation at the EIC","feed_subtitle":"In the DIS target fragmentation region, the cos 2φ asymmetry is set by a gluon-to-quark NEC ratio, making the saturation scale the lever.","key_machinery":"The carrying object is the adjoint dipole $S$-matrix $S_{x_g}(\\mathbf{r}_\\perp,\\mathbf{b}_\\perp)=\\frac{1}{N_c^2-1}\\langle \\mathrm{Tr}[U(\\mathbf{b}_\\perp+\\mathbf{r}_\\perp/2)U^\\dagger(\\mathbf{b}_\\perp-\\mathbf{r}_\\perp/2)]\\rangle_{x_g}$ and its Fourier transform, the gluon dipole distribution $F_{x_g}(\\mathbf{k}_{g\\perp})$. The reality and $C$-evenness of this $S$-matrix is the mechanism that eliminates all spin-dependent gluon components at eikonal order, leaving $f_1^g$ and $h_1^{t,g}$. The derivation uses CGC shockwave amplitudes for a real gluon emitted before and after the target, whose squares and interference form the three diagrams of Fig. 5; only the interference diagram feeds $h_1^{t,g}$. The fracture functions are then converted into NECs by the inclusive energy sum rule, and the observable asymmetry is built from the ratio in Eq. (4.14), which isolates the dipole-size dependence.","core_discovery":"At eikonal accuracy the gluon dipole $S$-matrix is real and $C$-even, which forces all target-spin-dependent gluon components of NECs and fracture functions to vanish; only the two components of an unpolarized target, $f_1^g$ and $h_1^{t,g}$, survive. Both are expressed through the adjoint dipole gluon distribution $F_{x_g}(k_{g\\perp})$, the Fourier transform of the adjoint dipole $S$-matrix, with $h_1^{t,g}$ receiving contributions only from the interference of the two CGC emission amplitudes. The inclusive energy sum rule converts the resulting fracture functions into the NECs of Eqs. (3.30) and (3.34), and in the dilute large-$\\theta$ limit both recover the collinear gluon-splitting results with $1/\\theta^2$ scaling. The phenomenological claim is Eq. (4.14): at leading order in $\\alpha_s$ and small $x_B$, the $\\cos 2\\phi$ asymmetry of the TFR energy pattern is $-\\frac{\\alpha_s}{2\\pi}\\frac{T_F}{3}\\frac{B(y)}{A(y)}$ times the ratio of $\\sum_{q,\\bar q} e_q^2 x_B h_1^{t,g}(x_B,\\theta)$ to $\\sum_{q,\\bar q} e_q^2 x_B f_1^q(x_B,\\theta)$. The numerical study at $x_B=5\\times10^{-3}$, $Q^2=25\\,\\mathrm{GeV}^2$ finds this asymmetry negative, growing in magnitude with $\\theta$, and strongly suppressed in $e{+}\\mathrm{Au}$ relative to $e{+}p$.","pith_inferences":["Beyond the paper: the ratio form of the asymmetry means its nuclear suppression is largely insensitive to absolute normalization uncertainties in the dipole models, so even a first EIC measurement with moderate statistics could distinguish saturated from unsaturated dipole inputs.","Beyond the paper: the accuracy of the predicted $\\theta$-dependence rests on the untested replacement of the $z$-convolution in Eq. (4.13); numerically evaluating the exact convolution using Eq. (3.34) would be a decisive check of the quantitative claims.","Beyond the paper: the same adjoint-to-fundamental-dipole ratio structure should control $\\cos 2\\phi$ energy-pattern asymmetries in forward proton-nucleus collisions, where the target is gluon-dominated; the formalism developed here appears directly adaptable to that case.","Beyond the paper: since the vanishing of target-spin gluon NECs follows from the reality of the eikonal adjoint $S$-matrix, subeikonal corrections that introduce imaginary $C$-odd structures are the natural place to find the first nonzero spin-dependent gluon NECs, a direction the paper itself flags."],"forward_implications":["Across the $\\theta$ range, $f_1^g$ shows almost no nuclear suppression while $h_1^{t,g}$ and $f_1^q$ are suppressed, so the asymmetry's suppression is a cleaner saturation signal than the gluon NEC alone.","Because the asymmetry is a single-inclusive calorimetric energy flow, it can be measured without jet reconstruction or particle identification, and it avoids Sudakov suppression in contrast to TMD dijet asymmetries.","The predicted $R_{eA}<1$, most pronounced at small $\\theta$ where $\\theta Q \\lesssim Q_s$, gives a concrete EIC signature for the onset of gluon saturation.","In the dilute limit both gluon NECs match the collinear splitting kernels with $1/\\theta^2$ scaling, connecting the CGC result to the established collinear framework.","Together with the quark-sector NEC results, the gluon-sector calculation completes a unified eikonal-level tomography framework for the target fragmentation region at small $x$."],"supporting_citations":[{"why":"Supplies the unpolarized quark NEC $f_1^q$ from the fundamental dipole and the $x_g=x_B$ choice used in the numerical setup.","marker":"[103]"},{"why":"Provides the quark NEC $f_1^q$ used in comparisons and the odderon contrast showing why the gluon sector differs.","marker":"[104]"},{"why":"Gives the operator definitions, eight-component parametrizations, projections, the energy sum rule, and the structure-function factorization used in Sec. 4.","marker":"[79]"},{"why":"Previous eikonal result for the unpolarized gluon jet fracture function that this calculation reproduces in the overlapping limit and uses as a comparison point.","marker":"[117]"},{"why":"Introduced nucleon energy correlators and the TFR factorization formula quoted in Eq. (4.7).","marker":"[77]"},{"why":"Gluonic fracture-function decomposition and the hard kernel for $\\Sigma_{UU}^{\\cos 2\\phi}$ used in Eq. (4.12).","marker":"[76]"},{"why":"TFR factorization proof that underwrites the collinear factorization framework used for the energy pattern.","marker":"[73]"},{"why":"Supplies the rcBK evolution input and parameters ($\\gamma$, $\\Lambda$, $Q_{s0}^2$, $\\sigma_0$) used for the numerical predictions.","marker":"[131]"},{"why":"Supplies the relation between nuclear and proton transverse areas, Eq. (4.17), used in the $e{+}A$ predictions.","marker":"[65]"},{"why":"Collinear splitting-kernel results for $f_1^g$ and $h_1^{t,g}$ that the dilute large-$\\theta$ limit must reproduce.","marker":"[108]"}],"fun_headline_variants":["Gluon asymmetry ratio exposes saturation at the EIC","Cos 2φ asymmetry probes gluon saturation in TFR","Nuclear suppression of gluon asymmetry signals saturation","Gluon NEC ratio reveals saturation scale in DIS"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The numerical prediction for the asymmetry and its nuclear suppression depends on the approximation in Eq. (4.13) that replaces the integral over $z$ of $z h_1^{t,g}(z,\\theta)$ by $x_B h_1^{t,g}(x_B,\\theta)$, and on the hand-picked nuclear saturation scale range $3 Q_{s0}^2 < Q_{s0A}^2 < 5 Q_{s0}^2$; if these are not accurate, the magnitude and $\\theta$-dependence of the predicted $\\cos 2\\phi$ asymmetry could shift substantially.","fun_headline_variants_meta":{"raw":{"variants":["Gluon asymmetry ratio exposes saturation at the EIC","Cos 2φ asymmetry probes gluon saturation in TFR","Nuclear suppression of gluon asymmetry signals saturation","Gluon NEC ratio reveals saturation scale in DIS"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000892,"raw_usage":{"total_tokens":3944,"prompt_tokens":1137,"completion_tokens":2807,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":753,"completion_tokens_details":{"reasoning_tokens":2743}},"tokens_in":753,"tokens_out":2807,"duration_ms":18199,"temperature":1.0,"reasoning_tokens":2743,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:31:37.716274+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct evaluation of the $z$-convolution in Eq. (4.12) using the computed $h_1^{t,g}(z,\\theta)$ from Eq. (3.34), without the small-$x$ substitution, would show whether the predicted $\\cos 2\\phi$ asymmetry and its nuclear suppression survive; alternatively, at the EIC, measuring $\\langle\\cos 2\\phi\\rangle$ in $e{+}p$ and $e{+}Au$ at $x_B = 5\\times10^{-3}$, $Q^2=25\\,\\mathrm{GeV}^2$ and finding $R_{eA} \\ge 1$ in the small-$\\theta$ region would contradict the central prediction.","supporting_citations":[{"cited_title":"Gluonic contributions to semi-inclusive DIS in the target fragmentation region","cited_arxiv_id":"2402.15112","evidence_quote":"Gluonic fracture-function decomposition and the hard kernel for $\\Sigma_{UU}^{\\cos 2\\phi}$ used in Eq. (4.12)."}],"review_version":1}