REVIEW 3 major objections 5 minor 32 references
Nucleon axial, tensor, and scalar charges and $\sigma$-terms in lattice QCD
T0 review · 3 major / 5 minor · reviewed 2026-08-08 · deepseek-v4-flash
Pith's one-line read Using three ensembles with quarks at their physical masses, the paper claims continuum-limit nucleon charges and sigma-terms without chiral extrapolation, finding $g_A^{u-d}=1.250(24)$.
desk verdict Solid ETMC proceedings with the real new content in the preliminary E112 analysis; the continuum-limit values are shared with the companion paper and the charm sigma-term needs a stronger discretization systematic. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is three $N_f=2+1+1$ twisted-mass clover-improved ensembles with lattice spacings near 0.080, 0.068, and 0.057 fm and pion masses around 140 MeV, analyzed through two-state fits to two- and three-point correlation functions. Excited-state contamination is controlled by Akaike Information Criterion weighting over a wide fit-parameter space, and continuum-limit systematics are controlled by averaging a linear-in-$a^2$ extrapolation and two constant extrapolations, one omitting the coarsest ensemble. Nonperturbative renormalization is carried out in the RI'/MOM scheme with perturbative conversion to the $\overline{\mathrm{MS}}$ scheme at 2 GeV.
What would settle it
Add the fourth, finer ensemble at $a\approx 0.05$ fm with full statistics and complete disconnected contributions and repeat the continuum extrapolation; if the model-averaged values shift by more than the quoted uncertainties, the three-ensemble error model is inadequate. A simpler check is whether the linear and constant fits' intercepts agree within their errors.
Extended reading notes
Core claim
The central discovery asserted is that continuum-limit nucleon charges and $\sigma$-terms can be obtained using only physical-point ensembles, eliminating the chiral extrapolation that has dominated systematics in previous determinations. Concretely, the paper reports $g_A^{u-d}=1.250(24)$, compatible with experiment, $g_T^{u-d}=0.955(29)$, $g_S^{u-d}=1.08(31)$, $\sigma_{\pi N}=41.9(8.1)$ MeV, $\sigma_s=30(17)$ MeV, and $\sigma_c=82(29)$ MeV, together with single-flavor charges for up, down, strange, and charm. All values include valence and sea (disconnected) quark contributions and are renormalized nonperturbatively to the $\overline{\mathrm{MS}}$ scheme at 2 GeV.
Load-bearing premise
The results rest on the assumption that the true discretization error is bracketed by an AIC-weighted average of a linear-in-$a^2$ fit and two constant fits using only three lattice spacings; if the $a^2$ dependence is nonlinear at these spacings, or the coarsest ensemble lies outside the scaling region, the central values and errors would be biased.
Editorial extensions
If this is right
- The nucleon charges and sigma-terms can serve as direct QCD inputs for neutron decay, dark-matter direct detection, and searches for physics beyond the Standard Model.
- Agreement of $g_A^{u-d}$ with the measured value supports the claim that excited-state and continuum systematics are controlled at the physical point.
- The flavor-separated axial charges give the quark spin contribution to the nucleon spin, providing a target for future deep-inelastic-scattering measurements.
- The physical-point-only continuum strategy avoids chiral extrapolation and can be applied to other nucleon matrix elements, such as form factors and parton moments.
Reading between the lines
- The quoted uncertainty on the continuum limit rests on three lattice spacings; a full-statistics analysis of the finer $a\approx 0.05$ fm ensemble may shift or shrink the central values, especially for the sigma-terms.
- If confirmed, the sizable $\sigma_c=82(29)$ MeV implies a non-negligible charm-quark contribution to the nucleon mass, with consequences for dark-matter couplings and the proton mass decomposition.
- The success of this approach would suggest that the spread among published lattice values for these quantities is dominated by chiral extrapolation rather than by lattice artifacts.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analyzes three N_f=2+1+1 twisted-mass clover-improved ensembles (B64, C80, D96) with quark masses tuned near their physical values and extracts the nucleon axial, scalar, and tensor charges and the pion-nucleon, strange, and charm sigma-terms from two- and three-point correlation functions. The analysis uses two-state spectral fits, AIC model averaging over excited-state and continuum-extrapolation models, and nonperturbative RI'/MOM renormalization with conversion to MS at 2 GeV. The continuum limit is obtained from a model average of a linear-in-a^2 fit through all three ensembles, a constant fit through all three, and a constant fit omitting the coarsest B64 ensemble. A fourth finer ensemble E112 is used for preliminary isovector results only. The main outputs are g_A^{u-d}=1.250(24), g_S^{u-d}=1.08(31), g_T^{u-d}=0.955(29), sigma_piN=41.9(8.1) MeV, sigma_s=30(17) MeV, and sigma_c=82(29) MeV. The paper claims that this is the first continuum-limit determination of these quantities using only physical-point ensembles, avoiding chiral extrapolations.
Significance. If the quoted values are reliable, they provide a valuable QCD determination of nucleon charges and sigma-terms at the physical point, with the isovector axial charge agreeing with experiment and the tensor and scalar charges serving as inputs for searches of beyond-Standard-Model interactions. The use of physical-point ensembles avoids chiral extrapolation, the AIC framework is a modern and appropriate tool for systematic error evaluation, and the nonperturbative renormalization is a strength. However, the significance of the paper depends critically on the continuum extrapolation, which rests on only three lattice spacings, and on the treatment of the charm sigma-term, where a single constant extrapolation is used. These issues need to be addressed before the central claim can be considered fully supported.
major comments (3)
- [Section 3, Fig. 2] The continuum-limit claim is the central assertion of the paper, and it rests entirely on three lattice spacings (B64, C80, D96) spanning only a factor of about two in a^2. With these three points, the model set {linear in a^2, constant, constant omitting B64} cannot distinguish a genuine O(a^2) trend from scatter; indeed, the AIC weights for g_A (30% linear, 40% constant with B64, 30% constant without B64) show that the result is dominated by the constant fits. The quoted uncertainty on g_A and on the Table 2 entries therefore does not include a meaningful contribution from the unknown curvature of the discretization error. I ask the authors to add a quantitative stability test, for example a two-point linear extrapolation using C80 and D96, a comparison with the E112 point as a fourth spacing, or an explicit breakdown of the systematic error attributable to the choice of continuum-extrapolation model.
- [Section 3, Table 3] The quoted value sigma_c=82(29) MeV is obtained from a single constant extrapolation, as stated in Table 3. This is not consistent with the AIC-model-averaging strategy used for the other quantities and is difficult to justify: at the coarsest lattice spacing a=0.07957 fm, the dimensionless combination m_c a is of order 0.4, so an O((m_c a)^2) discretization correction is expected to be the dominant systematic effect. A constant model omits exactly this term, and the 29 MeV error bar therefore contains no contribution from the continuum-extrapolation systematic. Please either include a linear-in-a^2 model for sigma_c (possibly with a constrained prior), or provide a concrete argument that the charm scalar matrix element has negligible a^2 dependence on these ensembles, and add the resulting model uncertainty to the final value.
- [Section 4, Conclusions] The paper claims 'for the first time' a continuum-limit determination using only physical-point ensembles. This is a strong novelty claim. Given that the continuum limit is based on only three ensembles and that the fourth, finer E112 ensemble is presented only as preliminary and is not used in the extrapolation, the claim is stronger than the evidence shown. I recommend either including the E112 data at least as a cross-check for the isovector quantities, or softening the claim to explicitly state that the present result uses three physical-point ensembles and should be considered an intermediate step pending the finer ensemble.
minor comments (5)
- [Abstract] The abstract contains the typo 'the the latter'; it should read 'the latter'.
- [Section 3] The sentence 'The extracted values show almost no dependence on on t_low^s' contains a duplicated 'on'.
- [Fig. 2 caption] The caption refers to the 'C90 and D96 ensembles', but the ensemble is labeled C80 in Table 1; please correct this label.
- [Eqs. (2)-(4)] The notation for the three-point function is inconsistent: Eq. (3) uses C_mu, while Eq. (4) uses C_3pt. Please align the notation for clarity.
- [Tables 2 and 3] The results are quoted without a decomposition into statistical and systematic uncertainties. Since the paper's method is AIC model averaging, a statistical/systematic breakdown would make the size of the continuum-extrapolation systematic transparent and would help the reader assess the source of the quoted error bars.
Circularity Check
No circularity: the charges and sigma-terms are obtained from spectral fits to lattice correlation functions, with an AIC model average for excited states and continuum extrapolation; the only self-citation is a companion-paper reporting overlap, not a load-bearing derivation.
full rationale
The derivation chain is self-contained. Nucleon charges are extracted from two- and three-point correlation functions through the spectral decomposition in Eqs. (2)-(4); Eq. (4) is a standard large-time limit of a ratio, not an ansatz fitted to reproduce the quoted charges. Excited states are controlled by two-state fits with AIC model averaging, and the continuum limit uses a pre-specified set of linear-in-a^2 and constant fits with AIC weights; the data, not the target values, select the weights. The sigma-terms follow from Eq. (6) using separately computed scalar matrix elements (including disconnected loops) and quark masses, with the twisted-mass renormalization cancellation being a known property of the formulation. No equation reduces a claimed output to a fitted input by construction. The only self-reference is Ref. [16], cited for renormalization details and per-ensemble results; this is a companion-proceedings overlap, not a uniqueness theorem or an ansatz that carries the derivation. External benchmarks (g_A vs experiment and Fig. 3 comparisons with other collaborations) make the results falsifiable outside the paper's own fitted values. The short continuum lever arm and the single constant fit used for sigma_c are systematic-error concerns, not circularity.
Assumptions & free parameters
free parameters (2)
- Slope of the linear-in-a^2 continuum extrapolation (per observable) =
not quoted
- Two-state fit parameters (excited-state energy gaps and amplitudes, per ensemble) =
not quoted
assumptions (5)
- domain assumption Twisted-mass fermions at maximal twist with a clover term provide automatic O(a) improvement, so discretization errors begin at O(a^2).
- domain assumption The AIC model average over the chosen set of continuum extrapolations (linear, constant with B64, constant without B64) gives an unbiased estimate of the continuum limit and its systematic error.
- domain assumption RI'/MOM nonperturbative renormalization with perturbative conversion to the MS scheme at 2 GeV is valid for the axial, scalar, and tensor operators.
- domain assumption The quark masses and lattice spacings taken from Refs [8,9] are accurate enough that the physical-point condition holds and no chiral extrapolation is needed.
- ad hoc to paper A constant extrapolation is sufficient for the charm sigma-term sigma_c.
Cite this review
Pith. "Pith review of Nucleon axial, tensor, and scalar charges and $\sigma$-terms in lattice QCD." pith.science (2026). https://pith.science/paper/LLKK3J6O
@misc{pith2026250205274,
author = {Pith},
title = {Pith review of: Nucleon axial, tensor, and scalar charges and $\sigma$-terms in lattice QCD},
year = {2026},
howpublished = {\url{https://pith.science/paper/LLKK3J6O}},
note = {Machine review of arXiv:2502.05274}
}
abstract
We determine the nucleon axial, scalar and tensor charges at the continuum limit by analyzing three $N_f=2+1+1$ twisted mass fermion ensembles with all quark masses tuned to approximately their physical values. We include all contributions from valence and sea quarks. We use the Akaike Information Criterion to evaluate systematic errors due to excited states and the continuum extrapolation. For the nucleon isovector axial charge we find $g_A^{u-d}=1.250(24)$, in agreement with the experimental value. We compute the axial, tensor and scalar charges for each quark flavor. The axial charge provides crucial information on the intrinsic spin carried by quark in the nucleon and the the latter two provide input for experimental searches of physics beyond the standard model. Moreover, we extract the nucleon $\sigma$-terms and find $\sigma_{\pi N}=41.9(8.1)$ MeV, for the strange $\sigma_{s}=30(17)$ MeV and for the charm $\sigma_{c}=82(29)$ MeV. We also present preliminary results on the isovector quantities using a fourth ensemble at smaller lattice spacing.
Figures
Reference graph
Works this paper leans on
-
[16]
C. Alexandrou, S. Bacchio, J. Finkenrath, C. Iona, G. Koutsou, Y. Li et al.,Nucleon charges and𝜎-terms in lattice QCD, 2412.01535
-
[1]
General Neutrino Interactions at the DUNE Near Detector
I. Bischer and W. Rodejohann,General Neutrino Interactions at the DUNE Near Detector, Phys. Rev. D99(2019) 036006 [1810.02220]
work page Pith review arXiv 2019
-
[2]
COHERENTcollaboration,Observation of Coherent Elastic Neutrino-Nucleus Scattering, Science 357 (2017) 1123 [1708.01294]
arXiv 2017
-
[3]
A.Beda, V.Brudanin, V.Egorov, D.Medvedev, V.Pogosov, M.Shirchenkoetal.,Theresults of search for the neutrino magnetic moment in GEMMA experiment,Adv. High Energy Phys. 2012 (2012) 350150
work page 2012
-
[4]
T. Marrodán Undagoitia and L. Rauch,Dark matter direct-detection experiments,J. Phys. G 43 (2016) 013001 [1509.08767]
arXiv 2016
-
[5]
Search for Tensor Interactions in Kaon Decays at DA$\Phi$NE
M.V. Chizhov,Search for tensor interactions in kaon decays at DAPHNE,Phys. Lett. B381 (1996) 359 [hep-ph/9511287]
work page Pith review arXiv 1996
-
[6]
Dhargyal,New tensor interaction as the source of the observed CP asymmetry in 𝜏→𝐾𝑆𝜋𝜈𝜏,Springer Proc
L. Dhargyal,New tensor interaction as the source of the observed CP asymmetry in 𝜏→𝐾𝑆𝜋𝜈𝜏,Springer Proc. Phys.203 (2018) 329 [1610.06293]
arXiv 2018
-
[7]
C. Alexandrou et al.,Simulating twisted mass fermions at physical light, strange and charm quark masses, Phys. Rev.D98 (2018) 054518 [1807.00495]
arXiv 2018
Show all 32 references
-
[8]
Extended Twisted Mass collaboration, Lattice calculation of the short and intermediate time-distance hadronic vacuum polarization contributions to the muon magnetic moment using twisted-mass fermions, Phys. Rev. D107 (2023) 074506 [2206.15084]
2023 arXiv
-
[9]
Extended Twisted Mass Collaboration (ETMC) collaboration, Strange and charm quark contributions to the muon anomalous magnetic moment in lattice QCD with twisted-mass fermions, 2411.08852
-
[10]
Alpha collaboration, Lattice QCD with a chirally twisted mass term,JHEP 08 (2001) 058 [hep-lat/0101001]
2001 arXiv
-
[11]
Frezzotti and G.C
R. Frezzotti and G.C. Rossi,Chirally improving Wilson fermions. 1. O(a) improvement, JHEP 08 (2004) 007 [hep-lat/0306014]
2004 arXiv
-
[12]
Sheikholeslami and R
B. Sheikholeslami and R. Wohlert,Improved Continuum Limit Lattice Action for QCD with Wilson Fermions, Nucl. Phys.B259(1985) 572
1985
-
[13]
Extended Twisted Mass collaboration, Quark masses using twisted-mass fermion gauge ensembles,Phys. Rev. D104 (2021) 074515 [2104.13408]
2021 arXiv
-
[14]
Jay and E.T
W.I. Jay and E.T. Neil,Bayesian model averaging for analysis of lattice field theory results, Phys. Rev. D103 (2021) 114502 [2008.01069]
2021 arXiv
-
[15]
Neil and J.W
E.T. Neil and J.W. Sitison,Improved information criteria for Bayesian model averaging in lattice field theory,Phys. Rev. D109 (2024) 014510 [2208.14983]. 9 Nucleon axial, tensor, and scalar charges and𝜎-terms in lattice QCD Christos Iona
2024 arXiv
-
[17]
Bär,𝑁𝜋-state contamination in lattice calculations of the nucleon axial form factors, Phys
O. Bär,𝑁𝜋-state contamination in lattice calculations of the nucleon axial form factors, Phys. Rev. D99(2019) 054506 [1812.09191]
2019 arXiv
-
[18]
Bär,Nucleon-pion-state contribution in lattice calculations of the nucleon charges𝑔𝐴,𝑔𝑇 and𝑔𝑆, Phys
O. Bär,Nucleon-pion-state contribution in lattice calculations of the nucleon charges𝑔𝐴,𝑔𝑇 and𝑔𝑆, Phys. Rev. D94 (2016) 054505 [1606.09385]
2016 arXiv
-
[19]
B.Märkischetal., MeasurementoftheWeakAxial-VectorCouplingConstantintheDecayof Free Neutrons Using a Pulsed Cold Neutron Beam, Phys. Rev. Lett.122 (2019) 242501 [1812.04666]
2019 arXiv
-
[20]
Extended Twisted Mass collaboration, Nucleon axial and pseudoscalar form factors using twisted-mass fermion ensembles at the physical point,Phys. Rev. D109 (2024) 034503 [2309.05774]
2024 arXiv
-
[21]
Alexandrou et al.,Moments of the nucleon transverse quark spin densities using lattice QCD, Phys
C. Alexandrou et al.,Moments of the nucleon transverse quark spin densities using lattice QCD, Phys. Rev. D107 (2023) 054504 [2202.09871]
2023 arXiv
-
[22]
Alexandrou, S
C. Alexandrou, S. Bacchio, M. Constantinou, J. Finkenrath, K. Hadjiyiannakou, K. Jansen et al.,Nucleon axial, tensor, and scalar charges and𝜎-terms in lattice QCD,Phys. Rev. D 102 (2020) 054517 [1909.00485]
2020 arXiv
-
[23]
Djukanovic, G
D. Djukanovic, G. von Hippel, H.B. Meyer, K. Ottnad and H. Wittig,Improved analysis of isovector nucleon matrix elements with Nf=2+1 flavors of O(a) improved Wilson fermions, Phys. Rev. D109(2024) 074507 [2402.03024]
2024 arXiv
-
[24]
Precision Neutron Decay Matrix Elements (PNDME) collaboration, Nucleon isovector axial form factors, Phys. Rev. D109 (2024) 014503 [2305.11330]
2024 arXiv
-
[25]
RQCD collaboration, Octet baryon isovector charges from Nf=2+1 lattice QCD,Phys. Rev. D 108 (2023) 034512 [2305.04717]
2023 arXiv
-
[26]
QCDSF/UKQCD/CSSMcollaboration, Constraining beyond the standard model nucleon isovector charges, Phys. Rev. D108 (2023) 094511 [2304.02866]
2023 arXiv
-
[27]
Nucleon Matrix Elements (NME) collaboration, Precision nucleon charges and form factors using (2+1)-flavor lattice QCD,Phys. Rev. D105(2022) 054505 [2103.05599]
2022 arXiv
-
[28]
Walker-Loud et al.,Lattice QCD Determination of𝑔𝐴, PoS CD2018(2020) 020 [1912.08321]
A. Walker-Loud et al.,Lattice QCD Determination of𝑔𝐴, PoS CD2018(2020) 020 [1912.08321]
2020 arXiv
-
[29]
Liang, Y.-B
J. Liang, Y.-B. Yang, T. Draper, M. Gong and K.-F. Liu,Quark spins and Anomalous Ward Identity, Phys. Rev. D98(2018) 074505 [1806.08366]
2018 arXiv
-
[30]
𝜒QCDcollaboration, Nucleon isovector scalar charge from overlap fermions, Phys. Rev. D 104 (2021) 094503 [2103.12933]. 10 Nucleon axial, tensor, and scalar charges and𝜎-terms in lattice QCD Christos Iona
2021 arXiv
-
[31]
Jülich Supercomputing Centre,JUWELS: Modular Tier-0/1 Supercomputer at the Jülich Supercomputing Centre, Journal of large-scale research facilities5 (2019) A135
2019
-
[32]
Jülich Supercomputing Centre,JUWELS Cluster and Booster: Exascale Pathfinder with Modular Supercomputing Architecture at Juelich Supercomputing Centre, Journal of large-scale research facilities7 (2021) A183. 11
2021
Reviewed August 8, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.