{"id":"db05a101-8a6f-47eb-bb8c-021f95487b40","arxiv_id":"2504.15769","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"The paper derives spin-decomposed Wilson coefficients for the twist-2 gluon condensate and uses them to argue that charmonium spin fractions shift by a few percent near T_c, with the total spin unchanged by definition.","lead":"This paper uses QCD sum rules to estimate how the spin, orbital, and gluon pieces of charmonium states change just above the QCD critical temperature. The author reports small compensating shifts, but the headline cancellation is an automatic consequence of how the spin fractions are defined.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed thermal spin shifts rest entirely on the new twist-2 Wilson coefficients in Eqs. (7)-(12), yet no derivation or independent cross-check is provided, leaving the paper's only non-tautological numerical content unverified.","rationale":"The reader's REJECT is appropriate and my analysis does not move it. I agree with the reader that the total-spin-preservation sentence is an identity following from Eqs. (16)-(18), and that the new coefficients are presented without derivation. I would go further: the condensate-only temperature dependence (Sec. IV) is a standard assumption in thermal QCDSR and is not the weakest link. The weakest link is that the only quantitative predictions — the signs and few-percent sizes of the individual spin-component shifts in Figs. 2 and 3 — are linear in the newly introduced C_{4,2} coefficients, and those coefficients are asserted, not shown. This is not an internal inconsistency that I can demonstrate from the text, but it is a load-bearing gap: even a single wrong sign in one of the J_n combinations would reverse the claimed trend. The paper honestly notes the omission of O(α_s) spin decomposition and the empirical Borel-window choices in Sec. V, which I credit; those are limitations, not the core problem. The decisive check is an independent OPE calculation of C_{4,2}, or a comparison of the sum of the component coefficients with the known total twist-2 coefficient. Until that is done, the component slopes are unsupported and the paper should not be accepted as is.","tokens_in":7586,"tokens_out":8290,"duration_ms":76603,"concrete_test":"Independently compute the OPE of the right-circular component Π_{12}(ω)/ω in Eq. (5) from the three-point correlator in Eq. (1), using standard background-field or diagrammatic techniques at leading order in α_s, and compare the coefficient of G_2 with Eqs. (7)-(12) for both channels. A decisive sub-check is to sum Eqs. (7)-(9) and (10)-(12) and compare with the known total dimension-4 twist-2 coefficient for the vector and axial two-point functions in a thermal medium; a mismatch would prove the spin decomposition inconsistent. If the coefficients match, the concern is resolved and the component trends stand.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central numerical content of the paper is the few-percent rise of s_q and compensating fall of l_q between 1.00 and 1.04 T_c in Figs. 2 and 3. These slopes are generated entirely by the new dimension-4 twist-2 Wilson coefficients C_{4,2}^{S_q}, C_{4,2}^{L_q}, and C_{4,2}^{J_g} in Eqs. (7)-(12), multiplied by the temperature dependence of G_2(T). The coefficients are introduced in Sec. II as 'newly computed' and then displayed, but no derivation, diagram set, or independent check is given; Appendix A lists only the older G_0 coefficients. Because the spin fractions are defined by Eqs. (16)-(18) as ratios of components of the same total, s_q + l_q + j_g = 1 is an identity, so the advertised 'total spin preservation' is not a testable prediction. The only falsifiable statements are the individual component slopes, and those are linear in the unverified coefficients. An error in any J_n combination in Eqs. (7)-(12) would change or reverse the predicted trend while remaining invisible in the stated cancellation. One also needs the operator identity Eq. (4), asserted as a confirmation rather than proven, for the decomposition itself to be meaningful. Thus, unless the coefficients (or Eq. (4)) are independently established, the paper's numerical conclusion is unsupported.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper extends a previous QCD sum-rule analysis of the spin decomposition of spin-1 charmonium (Ref. [20]) to finite temperature near T_c. It introduces a new contribution from the dimension-4 twist-2 gluon condensate G_2 to the rotating-frame two-point function of vector and axial-vector currents, with the relevant Wilson coefficients displayed in Eqs. (7)-(12). From Borel-transformed sum rules, the authors define quark spin, quark orbital, and gluon angular-momentum fractions and compute their temperature dependence between 1.00 T_c and 1.04 T_c. The reported behavior is that the quark spin fraction rises by a few percent, the quark orbital fraction falls by a comparable amount, the gluon fraction is nearly unchanged, and the total is preserved.","tokens_in":7910,"tokens_out":6289,"duration_ms":57334,"significance":"If the new coefficients and the operator identity in Eq. (4) were independently established, the paper would provide a first finite-temperature QCD sum-rule estimate of the spin decomposition of charmonium near T_c, complementing recent lattice studies. The presentation is clear and the Borel-window treatment follows standard practice; the authors are also explicit about the neglect of O(alpha_s) spin decomposition. However, the central advertised conclusion that the total spin is preserved is a normalization identity following from the definitions, and the nontrivial numerical content is carried entirely by unverified Wilson coefficients. The paper is therefore not acceptable in its present form.","major_comments":[{"comment":"The claim that thermal changes cancel and the total spin is preserved is tautological: because s_q, l_q, and j_g are each defined as a spin-decomposed component of the OPE divided by the same total M(M^2,s_0), the sum s_q + l_q + j_g equals 1 identically for every Borel mass and temperature. The abstract's statement that 'these thermal changes cancel each other, ensuring that the total spin is preserved even at finite temperature' is therefore not a dynamical prediction. The paper should reframe the individual component shifts as the actual predictions and explicitly state that the sum is a normalization constraint.","section":"Sec. III, Eqs. (16)-(18); Abstract and Sec. IV"},{"comment":"The numerical temperature dependence of all three spin components is generated by the newly introduced Wilson coefficients C_{4,2}^{S_q}, C_{4,2}^{L_q}, and C_{4,2}^{J_g}. These coefficients are displayed in Eqs. (7)-(12) with the statement that they are 'newly computed', but no derivation, no list of contributing diagrams, and no consistency check is provided anywhere in the paper; Appendix A contains only the older G_0 coefficients. Since an error in any one of the J_n combinations would change or reverse the predicted slopes in Figs. 2 and 3 while leaving the advertised 'total spin preservation' unaffected, these coefficients are load-bearing. A derivation, or at least a check against known limits, must be added before the numerical conclusions can be assessed.","section":"Sec. II, Eqs. (7)-(12)"},{"comment":"The identity in Eq. (4), which relates the rotating-frame two-point function to the inertial-frame two-point function and is used to justify the spin-orbital decomposition, is asserted without proof ('we confirm'). This identity is needed for the physical meaning of M_{S_q}, M_{L_q}, and M_{J_g}; without a derivation or a reference containing one, the interpretation of the fractions in Eqs. (16)-(18) is not established.","section":"Sec. II, Eq. (4)"},{"comment":"The numerical analysis assumes that all temperature dependence is encoded in the gluon condensates G_0(T) and G_2(T), while the Wilson coefficients remain temperature-independent. This is a standard assumption when T is well below the OPE scale [25], but the calculation is performed very close to T_c, where that condition is less secure. Because the predicted shifts are only a few percent and are linear in G_2(T), the paper should quantify the sensitivity to this assumption, for example by comparing with the size of possible T-dependent Wilson-coefficient corrections or by varying the condensate input within the lattice uncertainties.","section":"Sec. IV"}],"minor_comments":[{"comment":"The vertical axes are labeled '(%)', but the j_g panels use inconsistent tick labels such as '0.' and '1.'; a uniform notation would improve readability.","section":"Figs. 2 and 3"},{"comment":"The notation '1 ++' with a space is typographically inconsistent; please use '1^{++}' consistently.","section":"Throughout the text"},{"comment":"Table I reports j_g = -1.5 +/- 0.7% for the axial vector channel, yet the text describes the gluon contribution as 'nearly unchanged' without discussing the negative sign; a negative gluon angular-momentum fraction deserves a physical comment.","section":"Table I and Sec. IV"},{"comment":"The stated uncertainties reflect only the spread over the Borel window and not the uncertainties in m_c, alpha_s, s_0, or the condensate values; the paper mentions this in passing, but it should be stated directly with Table I.","section":"Sec. IV and Table I"},{"comment":"The Borel window and continuum threshold are taken from an inertial-frame OPE that includes O(alpha_s) corrections, while the spin decomposition omits those corrections; this inconsistency is acknowledged but should be flagged more prominently near Eq. (15) to avoid overstating the precision of the extracted spin fractions.","section":"Sec. III and Appendix B"}],"recommendation":"major_revision","confidential_remarks":"I am recommending major revision rather than rejection because the missing derivations and the normalization issue are addressable in a revision: the authors can supply a derivation of Eqs. (7)-(12), add independent consistency checks, and reframe the cancellation claim as a normalization identity. If the derivation cannot be provided, rejection would be appropriate. I also note that the paper's main advertised conclusion is tautological, so even if the coefficients are correct, the significance should be reframed around the individual spin-component shifts rather than the preservation of the total spin."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline result—that thermal changes in quark spin and orbital angular momentum cancel so total spin is preserved—is true by construction. Eqs. (16)–(18) define the three components as fractions of the same total M, so their sum is identically 1 at every temperature and every Borel mass. That part of the abstract and conclusion should be reframed; it is not a physical prediction.\n\nWhat is genuinely new is the set of spin-decomposed Wilson coefficients for the dimension-4 twist-2 gluon operator in Eqs. (7)–(12). I believe these are not in the cited literature, and if they are correct they give the first QCD sum rule estimate of how the spin decomposition of vector and axial-vector charmonium shifts near Tc. The paper inherits the standard machinery from Ref. [20], the vacuum numbers look consistent with that earlier work, and the author is honest about the limitation that O(alpha_s) corrections are not spin-decomposed.\n\nThe soft spots are real and concentrated in the one place that matters. The numerical slopes in Figs. 2 and 3 are generated entirely by the new twist-2 coefficients multiplied by the temperature dependence of G2(T). Those coefficients are displayed without derivation, no diagram set is given, and no independent check against known limits or alternative computation is offered. An error in any Jn combination would reverse or change the predicted trends while leaving the advertised cancellation intact. The assumption that all temperature dependence sits in the gluon condensates is standard for T well below the OPE scale, but near 1.0–1.04 Tc it deserves more scrutiny than a citation. Finally, the error bars shown are only standard deviations over the Borel window; they ignore uncertainties in the condensates, threshold, and quark mass, so the apparent significance of the few-percent shifts is overstated.\n\nNone of this makes the paper worthless. The calculation is compact, the decomposition is clean, and a referee who can verify Eqs. (7)–(12) will either confirm a useful result or catch an error. I would send this to peer review with a request for a derivation or a detailed cross-check of the new coefficients, for a reframed conclusion that presents the individual component slopes as the claim, and for propagated uncertainties. As it stands, I would not cite it in my own work until the coefficients are independently confirmed.","headline":"The total spin preservation claim is a normalization artifact, but the new twist-2 Wilson coefficients are a genuine calculation that deserves referee time.","tokens_in":8437,"tokens_out":2124,"would_cite":false,"duration_ms":19885,"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":"For both the vector (1−−) and axial-vector (1++) charmonium states, heating toward the critical temperature raises the quark-spin share by a few percent, lowers the quark orbital share by about the same amount, and leaves the gluon share…","keywords":["charmonium","QCD sum rules","spin decomposition","finite temperature","twist-2 gluon condensate","quark orbital angular momentum","gluon angular momentum","heavy quarkonium"],"falsifier":"A lattice QCD computation of the gauge-invariant spin decomposition of the vector or axial-vector charmonium at T/Tc = 1.00 and 1.04 that finds no few-percent rise in the quark-spin fraction, or a gluon angular momentum change larger than about one percentage point, would falsify the prediction. A shorter check is to include the O(αs) corrections with their spin decomposition in the same sum rules: if the shifts reverse sign or stop cancelling, the central claim fails.","tokens_in":7388,"feed_emoji":"🌀","tokens_out":13425,"duration_ms":110071,"temperature":0.7,"pith_summary":"This paper asks what happens to the spin budget of spin-1 charmonium states—the vector 1−− (J/ψ-like) and axial-vector 1++ partners—as the system is heated toward the quark-gluon plasma transition. Using QCD sum rules, the author computes the dimension-4 twist-2 gluon operator's contribution to the two-point function of heavy vector and axial currents in a rotating frame, the piece that carries the medium's anisotropy. The central finding is a compensation: between T/Tc = 1.00 and 1.04, the quark spin fraction rises by a few percent, the quark orbital angular momentum fraction falls by a comparable amount, and the gluon angular momentum fraction barely changes. Because these thermal shifts cancel, the total spin of each charmonium state is preserved at finite temperature. The result matters because it turns the spin structure of quarkonia into a concrete, parameter-light observable of the hot medium, with the entire thermal response entering through the temperature-dependent gluon condensates.","feed_headline":"Near Tc, charmonium trades orbital spin for quark spin","feed_subtitle":"The total spin stays fixed even as its pieces redistribute near the quark-gluon plasma transition.","key_machinery":"The load-bearing object is the spin-decomposed two-point function in a rotating frame, $\\Pi^{\\mu\\nu}(q)=i\\int d^4x\\,d^4y\\,e^{iqx}\\langle 0|T\\{j^\\mu(x) M_z(y) j^\\nu(0)\\}|0\\rangle$, with $M_z$ the z-component of the total angular momentum density of quarks and gluons. The gauge-invariant decomposition $M_z = \\frac{1}{2}\\bar\\psi\\gamma^\\perp\\gamma_5\\psi + \\psi^\\dagger(\\vec x \\times (-i\\vec D))\\psi + \\vec x \\times (\\vec E \\times \\vec B)$ separates quark spin, quark orbital angular momentum, and gluon total angular momentum, and the identity of Eq. (4) connects the rotating-frame correlator to the inertial-frame correlator through spin and orbital derivative operators, proving total angular momentum conservation at the OPE level. The new ingredient is the dimension-4 twist-2 gluon condensate $G_2$, defined by $\\langle 0|\\frac{\\alpha_s}{\\pi} G_{\\mu\\alpha}G_{\\nu}^{\\ \\alpha}|0\\rangle = g_{\\mu\\nu}G_0 + (u_\\mu u_\\nu - \\frac{1}{4}g_{\\mu\\nu})G_2$, whose Wilson coefficients $C_{4,2}$ in the vector and axial channels are computed here for the first time. A Borel transform with continuum subtraction converts this OPE into sum rules for the individual fractions $s_q$, $\\ell_q$, $j_g$, whose temperature dependence is carried by $G_0(T)$ and $G_2(T)$.","core_discovery":"The paper's discovery is that the spin-decomposed operator product expansion for spin-1 charmonium can be carried to dimension 4 in a medium by adding the twist-2 gluon condensate $G_2$, with analytic Wilson coefficients for the quark-spin, quark-orbital, and gluon parts given in Eqs. (7)--(12). Using Borel-transformed QCD sum rules with these coefficients and lattice inputs for $G_0(T)$ and $G_2(T)$, the paper finds that in both the $1^{--}$ and $1^{++}$ channels, as $T/T_c$ rises from 1.00 to 1.04, the quark-spin share increases by a few percent, the quark orbital angular momentum share decreases by about the same amount, and the gluon share stays essentially constant. The three changes sum to zero, so the total spin is conserved at finite temperature. The author reads this as orbital angular momentum being transferred into quark spin as the heavy-quark binding weakens, and expects the qualitative trend to survive the missing $O(\\alpha_s)$ corrections because those corrections are nearly temperature-independent below the OPE separation scale.","pith_inferences":["Beyond the paper, the same compensation mechanism predicts that bottomonium spin-1 states will show much smaller thermal spin reshuffling than charmonium, since condensate contributions are suppressed by powers of $1/m_b$; a parallel sum-rule calculation would quantify the effect.","Beyond the paper, the newly computed twist-2 Wilson coefficients allow the sum rules to be inverted: lattice-QCD measurements of charmonium spin fractions near $T_c$ could be used to extract the temperature dependence of the twist-2 gluon condensate $G_2$, which is otherwise hard to isolate.","Beyond the paper, the near-exact cancellation suggests that total quarkonium spin is a more robust observable near the transition than any single component, so measurements of quarkonium polarization in heavy-ion collisions could see little net change even while the microscopic spin budget shifts.","Beyond the paper, recomputing the same OPE at nonzero three-momentum or in a moving medium would test whether the compensation persists outside the rest frame; the conservation identity of Eq. (4) suggests it should, but the Wilson coefficients would need to be re-derived."],"forward_implications":["Across T/Tc = 1.00–1.04, both the vector and axial-vector channels show the same pattern: quark spin grows by a few percent, quark orbital angular momentum shrinks by a comparable amount, and the gluon share is essentially flat.","The individual thermal shifts cancel, so the total spin sum $s_q+\\ell_q+j_g$ stays at unity at every temperature considered, meaning total angular momentum conservation holds in the finite-temperature sum rule.","The new $C_{4,2}$ Wilson coefficients complete the leading-order spin-decomposed OPE for spin-1 quarkonia in a medium; future sum-rule studies of quarkonium at finite temperature can reuse them directly.","Because the Wilson coefficients are assumed temperature-independent, the magnitude of the predicted shifts is tied to the condensate changes $G_0(T)$ and $G_2(T)$; different lattice condensate inputs would rescale the effect without changing its compensating structure."],"supporting_citations":[{"why":"Previous vacuum spin decomposition of spin-1 quarkonia in a rotating frame; supplies the method, the scalar-gluon condensate Wilson coefficients, and the vacuum spin fractions that this paper extends to finite temperature.","marker":"[20]"},{"why":"Same rotating-frame QCD sum-rule approach applied to the proton spin decomposition; establishes the decomposition methodology for other hadrons.","marker":"[21]"},{"why":"Finite-temperature charmonium QCD sum rules; source of the charm quark mass, strong coupling, Borel windows, thresholds, and the temperature-dependent condensate inputs.","marker":"[17]"},{"why":"Justifies the assumption that Wilson coefficients are temperature-independent below the OPE separation scale, making the gluon condensates the sole carriers of temperature dependence.","marker":"[25]"},{"why":"Pure SU(3) lattice gauge thermodynamics; provides the numerical values of the gluon condensates G0(T) and G2(T) used in the finite-temperature analysis.","marker":"[26]"},{"why":"Gauge-invariant decomposition of hadron spin into quark spin, quark orbital angular momentum, and gluon angular momentum; defines the Mz decomposition used here.","marker":"[23]"},{"why":"Standard QCD sum-rule review whose O(αs)-corrected OPE determines the Borel window and continuum threshold adopted in the spin decomposition.","marker":"[15]"}],"fun_headline_variants":["Charmonium spin splits: orbital gives way to quark spin near Tc","At Tc, charmonium's orbital spin flows into quark spin","Quark spin wins, orbital loses: charmonium near Tc","Total spin intact as charmonium pieces reshuffle near Tc","Spin conservation holds in charmonium as T rises to Tc"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that temperature enters the calculation only through the strengths of the gluon background fields, $G_0(T)$ and $G_2(T)$, while the computed coefficients multiplying those fields are unchanged by temperature; if those coefficients shift measurably near the critical temperature, the predicted few-percent compensation between quark spin and orbital angular momentum would not follow.","fun_headline_variants_meta":{"raw":{"variants":["Charmonium spin splits: orbital gives way to quark spin near Tc","At Tc, charmonium's orbital spin flows into quark spin","Quark spin wins, orbital loses: charmonium near Tc","Total spin intact as charmonium pieces reshuffle near Tc","Spin conservation holds in charmonium as T rises to Tc"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000778,"raw_usage":{"total_tokens":3402,"prompt_tokens":871,"completion_tokens":2531,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":487,"completion_tokens_details":{"reasoning_tokens":2441}},"tokens_in":487,"tokens_out":2531,"duration_ms":15044,"temperature":1.0,"reasoning_tokens":2441,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:18:20.516334+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A lattice QCD computation of the gauge-invariant spin decomposition of the vector or axial-vector charmonium at T/Tc = 1.00 and 1.04 that finds no few-percent rise in the quark-spin fraction, or a gluon angular momentum change larger than about one percentage point, would falsify the prediction. A shorter check is to include the O(αs) corrections with their spin decomposition in the same sum rules: if the shifts reverse sign or stop cancelling, the central claim fails.","supporting_citations":[{"cited_title":"Spin-1 quarkonia in a rotating frame and their spin contents","cited_arxiv_id":"2212.14570","evidence_quote":"Previous vacuum spin decomposition of spin-1 quarkonia in a rotating frame; supplies the method, the scalar-gluon condensate Wilson coefficients, and the vacuum spin fractions that this paper extends to finite temperature."},{"cited_title":"Proton Spin Decomposition via QCD Sum Rules","cited_arxiv_id":"2503.20225","evidence_quote":"Same rotating-frame QCD sum-rule approach applied to the proton spin decomposition; establishes the decomposition methodology for other hadrons."},{"cited_title":"Heavy quarkonium correlators at finite temperature: QCD sum rule approach","cited_arxiv_id":"0908.2856","evidence_quote":"Finite-temperature charmonium QCD sum rules; source of the charm quark mass, strong coupling, Borel windows, thresholds, and the temperature-dependent condensate inputs."},{"cited_title":"Hatsuda, Y","cited_arxiv_id":null,"evidence_quote":"Justifies the assumption that Wilson coefficients are temperature-independent below the OPE separation scale, making the gluon condensates the sole carriers of temperature dependence."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Standard QCD sum-rule review whose O(αs)-corrected OPE determines the Borel window and continuum threshold adopted in the spin decomposition."}],"review_version":1}