REVIEW 5 minor 24 references
The spin-two/scalar response ratio in quarkonium-nucleon scattering factorizes into a simple multiplicative rescaling, cleanly separating quarkonium probe physics from nucleon scalar and gravitational structure.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-01 14:13 UTC pith:NQFSQO3S
load-bearing objection A careful, honest RG-consistent operator benchmark for quarkonium–nucleon scattering, but the headline factorization is an algebraic consequence of the coefficient parametrization, not a dynamical prediction.
Chromoelectric and chromomagnetic matching to scalar and spin-two nucleon structure
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is Eq. (38): R_{2/0}^Phi(t; rho) = [(1+rho)/(1-rho)] R_{2/0}^Phi(t; 0). Within a gluon-only local dimension-four matching condition, the scalar Wilson coefficient is proportional to (1-rho) and the spin-two coefficient to (1+rho), and because both RG evolution maps are linear, the chromomagnetic dependence factors out as a multiplicative rescaling. This identity is independent of the absolute quarkonium coupling, the target matrix elements, and the common scale, and holds in the threshold-aligned symmetric kinematics where the quarkonium velocity is parallel to the average nucleon momentum.
What carries the argument
The construction rests on two operator identities: the QCD trace identity expressing the gluonic scalar operator in terms of the RG-invariant total trace and the quark-mass operator, and the leading-logarithmic singlet mixing matrix for quark and gluon twist-two operators. The trace identity fixes a correlated quark-mass coefficient away from the matching scale; the mixing matrix induces a quark spin-two coefficient. Together they make the full scalar coefficient vector proportional to (1-rho) and the full spin-two vector proportional to (1+rho). The aligned velocity projection v^mu = P^mu/sqrt(P^2) reduces the traceless spin-two matrix element to the combination (3/4)M A_i(t) + (t/(16M))(3B
Load-bearing premise
The entire factorization rests on the gluon-only matching condition at the quarkonium scale: no independent quark scalar or spin-two operators, no higher-dimensional or nonlocal effects, and no heavy-flavor threshold crossovers within the leading-log evolution.
What would settle it
A lattice or experimental determination of the quarkonium-nucleon amplitude with both chromoelectric and chromomagnetic polarizabilities that shows a deviation from the (1+rho)/(1-rho) rescaling for fixed t, or the appearance of a non-zero quark spin-two Wilson coefficient at the matching scale, would falsify Eq. (38).
If this is right
- The ratio formula provides a controlled benchmark for quarkonium-nucleon elastic scattering without conflating probe matching with nucleon structure.
- The factorization identifies the combination 3B_i(t)-D_i(t) as the spin-two observable in aligned kinematics, distinct from the forward A_i(0) alone.
- Lattice or phenomenological extraction of the chromomagnetic-to-chromoelectric ratio rho can be translated directly into predictions for the spin-two/scalar response ratio.
- An independent chromomagnetic polarizability is shown to be necessary and sufficient to span the general CP-even, spin-independent, local two-gluon interaction at dimension four.
- The off-diagonal helicity representation for nonzero t is shown to be merely a kinematic effect of antiparallel momentum axes, not a new spin-flip structure.
Where Pith is reading between the lines
- If Eq. (38) survives tests, the same multiplicative rescaling should hold for other hadronic targets in the aligned-velocity limit, suggesting a universality of the rho dependence across targets.
- Measuring the ratio for two different quarkonium states with different rho values could isolate R(t;0) without knowing the absolute normalization, providing a practical extraction strategy.
- Beyond the local truncation, nonlocal or derivative gluon operators would introduce t-dependent corrections that break the simple multiplicative form; the formula's failure at larger |t| could map the domain of the local OPE.
- The same linearity argument extends to higher Mellin moments, where the mixing matrix has multiple non-zero eigenvalues, yielding a matrix-valued factorization.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper constructs an RG-consistent matching of compact heavy quarkonium to nucleon scalar and spin-two operators. At the matching scale the interaction is expressed in the chromoelectric/chromomagnetic basis with coefficients C_2 = α_E(1+ρ), C_S = -α_E(1-ρ). The scalar gluon operator is rewritten through the QCD trace identity into an invariant basis (Eqs. 13-19), while the spin-two sector is evolved at leading-logarithmic order in fixed-flavor MS using the N=2 singlet anomalous dimension matrix (Eqs. 22-26). The authors then derive the aligned-kinematics nucleon projection containing A_i(t) and 3B_i(t)-D_i(t), including an off-diagonal helicity representation in the Breit frame (App. B). Under the explicitly stated gluon-only, local dimension-four matching condition, the ratio of the spin-two to scalar reduced amplitudes factorizes as R_{2/0}(ρ) = [(1+ρ)/(1-ρ)] R_{2/0}(0) (Eq. 38). The paper states its scoping assumptions and consistency limits in Sec. IV.
Significance. If correct, the result provides a clean separation between quarkonium-dependent matching coefficients and target-dependent scalar/gravitational form factors, and it gives a controlled operator benchmark for lattice or phenomenological extractions of ρ. The derivations are internally coherent: the trace-identity transformations, the LL singlet evolution, the canonical-spin/helicity reduction, and the forward-limit checks (including the LMS normalization) are all explicit. The main caveat is that Eq. (38) is largely an algebraic consequence of the coefficient parametrization together with the stated gluon-only dimension-four assumptions; it is a consistency relation rather than a new dynamical prediction. The paper is honest about this, which is a strength rather than a weakness.
minor comments (5)
- [Sec. II, after Eq. (5)] The statement that evanescent finite terms do not modify the singlet anomalous-dimension matrix at leading-logarithmic accuracy is asserted but not justified. A short explanation or a reference would make the paper more self-contained.
- [Sec. III, Eq. (37)] The ratio is called 'signed' but the definition does not discuss the possibility that the denominator vanishes at isolated t values. Please state explicitly that Eq. (38) holds only where the scalar denominator is nonzero, as already partially noted.
- [Sec. IV, after Eq. (47)] The translation to the LMS forward matrix element is presented in one sentence. Showing the intermediate coefficient identification (e.g., how C_Θ is expressed in terms of c_E and c_B before substituting) would make the normalization check easier to verify.
- [Eq. (24) and App. A] The β0 coefficient uses T_F; please define T_F=1/2 and state explicitly that the fixed-n_f convention is used throughout the LL evolution, including any change in β0 when n_f changes.
- [References] References [18] and [23] are arXiv preprints; if the journal requires published versions, please update or note their status.
Circularity Check
No significant circularity: Eq. (38) is a derived algebraic corollary, not a fitted prediction; the only self-citation is illustrative and non-load-bearing.
full rationale
The paper's central claim, Eq. (38), is a derived consequence of the explicitly stated coefficient parametrization and matching assumptions, not a circular restatement. Equations (9), (14), (15), (25), and (26) fix the scalar and spin-two coefficient vectors to be proportional to (1−ρΦ) and (1+ρΦ), respectively; because the evolution is linear and ρΦ-independent, the ratio in Eq. (37) factorizes as written. RΦ2/0 enters only after being defined, and ρΦ is left as an undetermined input, so no fitted parameter is being renamed as a prediction. The only self-citation, Ref. [18], supplies an illustrative forward-range input for Fig. 1 and is explicitly not load-bearing; the paper states that Eq. (38) is 'not a model for ρΦ'. The trace-identity scalar transport, spin-two DGLAP evolution, and kinematic projections are derived from external standard results or explicit appendices, and the stated limitations (independent quark coefficients, higher-dimensional/nonlocal operators, J/ψ validity, fixed-nf truncation) are acknowledged in Section IV rather than hidden. The derivation is self-contained and no circular step reduces the central result to its own inputs.
Axiom & Free-Parameter Ledger
free parameters (3)
- rho_Phi = alpha_B^Phi / alpha_E^Phi
- mu_Phi (quarkonium matching scale)
- R_2/0^Phi(0;0) (pure-electric forward ratio) =
0.10--0.15 (illustrative)
axioms (6)
- domain assumption QCD trace identity Theta = (beta/(2g))[G^2] + (1+gamma_m)O_m, all orders in a mass-independent scheme (Eq. 11).
- domain assumption Leading-logarithmic singlet DGLAP evolution for the second Mellin moment with mixing matrix M2 (Eq. 22).
- ad hoc to paper Gluon-only matching at mu_Phi: quark scalar and quark spin-two Wilson coefficients vanish (Eqs. 13, 21).
- domain assumption Local dimension-four multipole expansion for compact quarkonium: radius smaller than gluon wavelength and suppressed color-octet gap (Section IV).
- domain assumption Aligned symmetric kinematics v_mu = P_mu/sqrt(P^2), v.Delta = 0 (Eq. 31).
- domain assumption Fixed-flavor MS scheme with trace subtraction before d->4 and minimal subtraction; evanescent terms do not affect LL anomalous dimensions (Section II).
read the original abstract
Compact heavy quarkonium couples through the multipole interaction to scalar and spin-two gluonic operators. At leading chromoelectric order the corresponding matching coefficients satisfy $C_2^\Phi=-C_S^\Phi$; an independent chromomagnetic polarizability lifts this relation within the general CP-even, spin-independent, local two-gluon interaction at dimension four and zero derivative order. We construct an RG-consistent realization in a fixed $MS$ convention. The QCD trace identity converts the gluon-only scalar matching condition into an invariant basis and fixes the correlated quark-mass coefficient required when the interaction is re-expressed in the scale-dependent basis away from the matching scale, whereas leading-logarithmic singlet evolution induces a quark spin-two coefficient. In threshold-aligned symmetric kinematics, the canonical-spin non-flip projection contains $A_i(t)$ and the combination $3B_i(t)-D_i(t)$. An explicit Breit-frame calculation relates this projection to an off-diagonal helicity representation for nonzero spacelike $t$; the off-diagonal form is kinematic rather than an additional dynamical spin flip. Linearity of the scalar and spin-two evolution factorizes the chromomagnetic dependence of their ratio as $R_{2/0}^{\Phi}(t;\rho_\Phi)=[(1+\rho_\Phi)/(1-\rho_\Phi)]R_{2/0}^{\Phi}(t;0)$ within the gluon-only dimension-four matching setup. The result separates state-dependent quarkonium matching from scalar and gravitational nucleon structure and states explicitly the assumptions under which this factorization holds.
Figures
Reference graph
Works this paper leans on
-
[1]
D. E. Kharzeev, Mass radius of the proton, Phys. Rev. D104, 054015 (2021), arXiv:2102.00110 [hep-ph]
Pith/arXiv arXiv 2021
-
[2]
This off-diagonal helicity matrix follows from comparing helicities defined along antiparallel momentum axes and does not represent an additional dynamical spin-flip structure
Thus the equal- helicity amplitudes vanish, whereas the two nonzero opposite-helicity amplitudes have a phase-convention-dependent relative sign. This off-diagonal helicity matrix follows from comparing helicities defined along antiparallel momentum axes and does not represent an additional dynamical spin-flip structure. Equation (32) is the canonical-spi...
-
[3]
M. E. Peskin, Short distance analysis for heavy quark systems. 1. diagrammatics, Nucl. Phys. B156, 365 (1979). 10
1979
-
[4]
Bhanot and M
G. Bhanot and M. E. Peskin, Short distance analysis for heavy quark systems. 2. applications, Nucl. Phys. B156, 391 (1979)
1979
-
[5]
M. B. Voloshin and V. I. Zakharov, Measuring QCD anomalies in hadronic transitions between onium states, Phys. Rev. Lett.45, 688 (1980)
1980
-
[6]
M. E. Luke, A. V. Manohar, and M. J. Savage, A QCD calculation of the interaction of quarkonium with nuclei, Phys. Lett. B288, 355 (1992), hep-ph/9204219
Pith/arXiv arXiv 1992
-
[7]
J.-W. Chen and M. J. Savage, Hadronic and electromagnetic interactions of quarkonia, Phys. Rev. D57, 2837 (1998), hep-ph/9710338
Pith/arXiv arXiv 1998
-
[8]
N. Brambilla, A. Pineda, J. Soto, and A. Vairo, Potential NRQCD: An effective theory for heavy quarkonium, Nucl. Phys. B566, 275 (2000), hep-ph/9907240
Pith/arXiv arXiv 2000
-
[9]
N. Brambilla, A. Pineda, J. Soto, and A. Vairo, Effective field theories for heavy quarkonium, Rev. Mod. Phys.77, 1423 (2005), hep-ph/0410047
Pith/arXiv arXiv 2005
-
[10]
N. Brambilla, G. Krein, J. Tarrus Castella, and A. Vairo, Long-range properties of1Sbottomonium states, Phys. Rev. D 93, 054002 (2016), arXiv:1510.05895 [hep-ph]
Pith/arXiv arXiv 2016
-
[11]
M. V. Polyakov and P. Schweitzer, Determination of the chromoelectric polarizability from lattice data, Phys. Rev. D98, 034030 (2018), arXiv:1801.08984 [hep-ph]
Pith/arXiv arXiv 2018
-
[12]
Y. Hatta, A. Rajan, and K. Tanaka, Quark and gluon contributions to the QCD trace anomaly, JHEP12(12), 008, arXiv:1810.05116 [hep-ph]
-
[13]
K. Tanaka, Three-loop formula for quark and gluon contributions to the QCD trace anomaly, JHEP01(01), 120, arXiv:1811.07879 [hep-ph]
-
[14]
T. Ahmed, L. Chen, and M. Czakon, A note on quark and gluon energy-momentum tensors, JHEP01, 077, arXiv:2208.01441 [hep-ph]
-
[15]
Ji, Gauge-invariant decomposition of nucleon spin, Phys
X.-D. Ji, Gauge-invariant decomposition of nucleon spin, Phys. Rev. Lett.78, 610 (1997), hep-ph/9603249
Pith/arXiv arXiv 1997
-
[16]
Ji, Deeply virtual compton scattering, Phys
X.-D. Ji, Deeply virtual compton scattering, Phys. Rev. D55, 7114 (1997), hep-ph/9609381
Pith/arXiv arXiv 1997
-
[17]
M. V. Polyakov and P. Schweitzer, Forces inside hadrons: Pressure, surface tension, mechanical radius, and all that, Int. J. Mod. Phys. A33, 1830025 (2018), arXiv:1805.06596 [hep-ph]
Pith/arXiv arXiv 2018
-
[18]
P. E. Shanahan and W. Detmold, Gluon gravitational form factors of the nucleon and the pion from lattice QCD, Phys. Rev. D99, 014511 (2019), arXiv:1810.04626 [hep-lat]
Pith/arXiv arXiv 2019
-
[19]
A. I. Syamtomov, Scalar and spin-two energy-momentum-tensor structure in near-threshold charmonium probes of the proton, arXiv e-prints (2026), arXiv:2606.08835 [hep-ph]
Pith/arXiv arXiv 2026
-
[20]
J. C. Collins, A. Duncan, and S. D. Joglekar, Trace and dilatation anomalies in gauge theories, Phys. Rev. D16, 438 (1977)
1977
-
[21]
Altarelli and G
G. Altarelli and G. Parisi, Asymptotic freedom in parton language, Nucl. Phys. B126, 298 (1977)
1977
-
[22]
A. Vogt, S.-O. Moch, and J. A. M. Vermaseren, The three-loop splitting functions in QCD: The singlet case, Nucl. Phys. B691, 129 (2004), hep-ph/0404111
Pith/arXiv arXiv 2004
-
[23]
Y. Guo, X. Ji, and F. Yuan, Proton’s gluon GPDs at large skewness and gravitational form factors from near-threshold heavy-quarkonium photoproduction, Phys. Rev. D109, 014014 (2024), arXiv:2308.13006 [hep-ph]
Pith/arXiv arXiv 2024
-
[24]
R. Ejima, D. Fujii, and M. Kawaguchi, Charmonium-nucleon femtoscopy as a possible probe of the nucleon gravitational form factor, arXiv e-prints (2026), arXiv:2607.00650 [hep-ph]
Pith/arXiv arXiv 2026
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.