Pith. sign in

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 →

arxiv 2502.05274 v1 pith:LLKK3J6O submitted 2025-02-07 hep-lat hep-exhep-phnucl-th

classification hep-lathep-exhep-phnucl-th PACS 12.38.Gc14.20.Dh
keywords nucleonaxialchargetensorscalarsigma-piNstrangesigmatermcharmlatticeQCDcontinuumlimittwistedmassfermions
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper aims to establish reliable continuum-limit values for the nucleon axial, tensor, and scalar charges and for the pion-nucleon, strange, and charm $\sigma$-terms directly from lattice QCD. Its specific claim is that using three ensembles with all quark masses tuned near their physical values, and combining a two-state excited-state analysis with Akaike Information Criterion model averaging, gives continuum limits without any chiral extrapolation. If the claim holds, the quoted numbers, including $g_A^{u-d}=1.250(24)$ and $\sigma_{\pi N}=41.9(8.1)$ MeV, are direct QCD-based predictions rather than extrapolations from heavier pion masses. The same strategy can be transferred to other nucleon matrix elements.

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.

Watch

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

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

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)
  1. [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.
  2. [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.
  3. [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)
  1. [Abstract] The abstract contains the typo 'the the latter'; it should read 'the latter'.
  2. [Section 3] The sentence 'The extracted values show almost no dependence on on t_low^s' contains a duplicated 'on'.
  3. [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.
  4. [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.
  5. [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

0 steps flagged · score 0.0 of 10

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 2 free parameters · 5 assumptions · 0 invented entities

The calculation is a first-principles numerical evaluation, but the quoted values depend on several modeling choices: the form of the continuum extrapolation, the excited-state fit model, the renormalization scheme conversion, and the scale setting from prior ETMC work. No physically new entity is introduced. The most ad hoc element is the constant extrapolation for sigma_c, which is not justified with a systematic comparison.

free parameters (2)
  • Slope of the linear-in-a^2 continuum extrapolation (per observable) = not quoted
    The linear extrapolation model fits a slope to three lattice-spacing points and receives a 30% AIC weight for g_A; its value and uncertainty are not reported but influence the continuum values.
  • Two-state fit parameters (excited-state energy gaps and amplitudes, per ensemble) = not quoted
    Excited-state removal relies on fitted energy differences and amplitudes over varied t_s windows; AIC weights over these fits affect every quoted charge.
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).
    Invoked in Section 2 and used to justify the linear-in-a^2 continuum extrapolation in Section 3.
  • 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.
    This is the main systematic-error model described in Section 3 and Fig. 2; it is a standard but not guaranteed method in lattice QCD.
  • 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.
    Stated in Section 2 with reference [16]; any error in the scheme conversion propagates into all quoted charges.
  • 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.
    Used in Section 1 and Table 1; the ensemble pion masses are 136-141 MeV, close to but not exactly the physical pion mass.
  • ad hoc to paper A constant extrapolation is sufficient for the charm sigma-term sigma_c.
    Section 3, Table 3: a single constant fit is used for sigma_c without an AIC model average or a stated justification.

how reviews work

0 comments
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

Figures reproduced from arXiv: 2502.05274 by the authors.

Figure 1
Figure 1. We present the ratio and fit results for all ensembles, for the isovector axial charge. In the legend, we give the symbols used to denote the various values of 𝑡𝑠/𝑎. The top row shows the analysis for 𝑔 𝑢−𝑑 𝐴 for the B64 ensemble, the middle row for C80 and the bottom for D96. The first column shows results on the ratio versus 𝑡𝑖𝑛𝑠 − 𝑡𝑠/2. The horizontal bands are the model averaged values. In the second column, we … view at source ↗
Figure 2
Figure 2. In the left panel, we show the continuum limit of the nucleon isovector axial charge (open symbol and band) extrapolated using the B64, C80 and D96 ensembles (filled symbols). The extrapolation is the result of a model average, which combines linear and constant fits as explained in the text. In the right panel, we show the weights for each type of fit, namely the linear extrapolation is represented by a blue circle… view at source ↗
Figure 3
Figure 3. Comparison of the results of this work with other lattice QCD results, for the isovector axial, scalar and tensor charges. Our results are shown with the red square and red error band. The blue triangles show previous ETMC results, for 𝑔 𝑢−𝑑 𝐴 [20] and for 𝑔 𝑢−𝑑 𝑇 [21], while Ref. [22] gives results on all three isovector charges including 𝑔 𝑢−𝑑 𝑠 for the B64 ensemble. Open symbols represent results without a contin… view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Ratio and preliminary fit results for the currently available data of the E112 ensemble, for the isovector charges. The notation is the same as in [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

32 extracted references · 6 canonical work pages

  1. [16]

    Alexandrou, S

    C. Alexandrou, S. Bacchio, J. Finkenrath, C. Iona, G. Koutsou, Y. Li et al.,Nucleon charges and𝜎-terms in lattice QCD, 2412.01535

  2. [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]

  3. [2]

    COHERENTcollaboration,Observation of Coherent Elastic Neutrino-Nucleus Scattering, Science 357 (2017) 1123 [1708.01294]

  4. [3]

    High Energy Phys

    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

  5. [4]

    Marrodán Undagoitia and L

    T. Marrodán Undagoitia and L. Rauch,Dark matter direct-detection experiments,J. Phys. G 43 (2016) 013001 [1509.08767]

  6. [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]

  7. [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]

  8. [7]

    Alexandrou et al.,Simulating twisted mass fermions at physical light, strange and charm quark masses, Phys

    C. Alexandrou et al.,Simulating twisted mass fermions at physical light, strange and charm quark masses, Phys. Rev.D98 (2018) 054518 [1807.00495]

Show all 32 references
  1. [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]

  2. [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

  3. [10]

    Alpha collaboration, Lattice QCD with a chirally twisted mass term,JHEP 08 (2001) 058 [hep-lat/0101001]

  4. [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]

  5. [12]

    Sheikholeslami and R

    B. Sheikholeslami and R. Wohlert,Improved Continuum Limit Lattice Action for QCD with Wilson Fermions, Nucl. Phys.B259(1985) 572

  6. [13]

    Extended Twisted Mass collaboration, Quark masses using twisted-mass fermion gauge ensembles,Phys. Rev. D104 (2021) 074515 [2104.13408]

  7. [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]

  8. [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

  9. [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]

  10. [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]

  11. [19]

    B.Märkischetal., MeasurementoftheWeakAxial-VectorCouplingConstantintheDecayof Free Neutrons Using a Pulsed Cold Neutron Beam, Phys. Rev. Lett.122 (2019) 242501 [1812.04666]

  12. [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]

  13. [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]

  14. [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]

  15. [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]

  16. [24]

    Precision Neutron Decay Matrix Elements (PNDME) collaboration, Nucleon isovector axial form factors, Phys. Rev. D109 (2024) 014503 [2305.11330]

  17. [25]

    RQCD collaboration, Octet baryon isovector charges from Nf=2+1 lattice QCD,Phys. Rev. D 108 (2023) 034512 [2305.04717]

  18. [26]

    QCDSF/UKQCD/CSSMcollaboration, Constraining beyond the standard model nucleon isovector charges, Phys. Rev. D108 (2023) 094511 [2304.02866]

  19. [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]

  20. [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]

  21. [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]

  22. [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

  23. [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

  24. [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

Pith tools

Reviewed August 8, 2026 · model on record in the stance chip above.