Pith. sign in

REVIEW 3 major objections 5 minor 1 cited by

Isovector axial and pseudoscalar form factors from twisted mass lattice QCD at the physical point

T0 review · 3 major / 5 minor · reviewed 2026-08-08 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read This paper reports the first continuum-limit extraction of the nucleon's isovector axial and pseudoscalar form factors using only simulations at the physical pion mass, yielding $g_A = 1.245(28)(14)$ and verifying PCAC and pion pole…

desk verdict A competent proceedings that re-renders the same collaboration's PRD results; the PCAC/PPD claims are partly fit-enforced and the novelty text overstates what is new. read the letter →

arxiv 2502.07583 v1 pith:CQKT3U6B submitted 2025-02-11 hep-lat

Constantia Alexandrou , Simone Bacchio , Martha Constantinou , Jacob Finkenrath , Roberto Frezzotti , Bartosz Kostrzewa , Giannis Koutsou , Gregoris Spanoudes
show 1 more author
Carsten Urbach
This is my paper · ORCID
classification hep-lat
keywords latticeQCDnucleonaxialformfactortwistedmassfermionsphysicalpioncontinuumlimitpoledominancePCACcharge
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 continuum-limit values for the nucleon's isovector axial, induced pseudoscalar, and pseudoscalar form factors directly from lattice QCD at the physical pion mass, without any chiral extrapolation. Using three twisted-mass fermion ensembles with lattice spacings from 0.08 fm to 0.057 fm, the authors extract the axial charge, axial radius, induced pseudoscalar coupling at the muon capture point, and the pion-nucleon coupling constant. They further show that the PCAC relation and pion pole dominance are restored in the continuum limit. If correct, this removes a major model-dependent systematic from the weak-interaction form factors needed for neutrino scattering predictions.

What carries the argument

The argument is carried by three $N_f=2+1+1$ twisted-mass fermion ensembles tuned to the physical point, with lattice spacings $a=0.07957$, $0.06821$, and $0.05692$ fm and approximately equal physical volumes, so that a continuum extrapolation linear in $a^2$ can be taken at fixed physics. Nucleon matrix elements are extracted from two- and three-point correlation functions using two- and three-state spectral fits over source-sink separations from 0.5 to 1.5 fm, with model averaging based on the Akaike Information Criterion. Renormalization is performed with fully non-perturbative Ward identities, and disconnected quark loops are neglected because they are $\mathcal{O}(a^2)$ for isovector quantities in the twisted-mass formulation and vanish in the continuum. The $Q^2$ dependence is described by dipole and z-expansion parameterizations whose parameters carry the linear $a^2$ cut-off correction.

What would settle it

A fourth ensemble at $a=0.049$ fm, which the paper reports as ongoing, would test the cut-off assumption: if a four-point linear $a^2$ extrapolation moves $g_A$, $\langle r_A^2\rangle$, $g_P^*$, or $g_{\pi NN}$ by more than the quoted errors, the linearity assumption or the error budget would be contradicted. For the excited-state assumption, extending the source-sink separations beyond 1.5 fm or including four-state fits would show whether the two-versus-three-state difference captures the full contamination.

Watch

Extended reading notes

Core claim

The central discovery is a set of continuum-limit values at the physical point: $g_A = 1.245(28)(14)$, $\langle r_A^2\rangle = 0.339(48)(06)$ $\mathrm{fm}^2$, $g_P^* = 8.99(39)(49)$, and $g_{\pi NN} = 13.25(67)(69)$, where the first error is statistical and the second is the excited-state systematic. In addition, the ratios $r_{\mathrm{PCAC}}$ and $r_{\mathrm{PPD}}$ deviate sizably from unity at finite lattice spacing but become consistent with 1 in the continuum limit, demonstrating that both the partially conserved axial-vector current relation and pion pole dominance hold without chiral extrapolation. The paper claims this is the first calculation to take the continuum limit using solely physical-pion-mass simulations.

Load-bearing premise

The results assume that the two- and three-state fits over source-sink separations up to about 1.5 fm remove all significant excited-state contamination, and that cut-off effects are linear in $a^2$ over the three lattice spacings.

Editorial extensions

If this is right

  • The method demonstrates that physical-point continuum limits for nucleon form factors are achievable with current twisted-mass ensembles, eliminating chiral extrapolation as a source of uncertainty.
  • The restoration of PCAC and pion pole dominance in the continuum limit confirms that the pion-pole behaviour of $G_P$ and $G_5$ is correctly reproduced once cut-off effects are removed.
  • The continuum values of $g_P^*$ and $g_{\pi NN}$ provide direct lattice inputs for neutrino-nucleus scattering and muon-capture phenomenology.

Reading between the lines

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

  • If the ongoing fourth ensemble at $a=0.049$ fm confirms the linear $a^2$ behaviour, the same three-ensemble physical-point strategy could be applied to other nucleon matrix elements, such as scalar and tensor charges, with the same reduction in systematic uncertainty.
  • The roughly one-standard-deviation gap between the lattice $g_A = 1.245(28)(14)$ and the experimental value near 1.276 suggests that a few percent of uncertainty remains to be explained; adding longer source-sink separations and fourth-state fits would directly test whether excited states are responsible.
  • Since pion pole dominance is restored in the continuum, the induced pseudoscalar coupling at the muon capture point is almost entirely fixed by $g_A$ and the axial radius through the pion pole; a more precise $g_A$ would therefore sharpen $g_P^*$ without a dedicated calculation.
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. This proceedings paper reports a lattice QCD determination of the isovector nucleon axial form factor G_A, the induced pseudoscalar form factor G_P, and the pseudoscalar form factor G_5 using three N_f=2+1+1 twisted-mass ensembles at the physical pion mass with lattice spacings 0.08, 0.068, and 0.057 fm. Excited-state contamination is treated with two- and three-state fits over sink-source separations from about 0.5 fm to 1.5 fm, followed by AIC model averaging. The form-factor shapes are fitted with dipole and z-expansion parameterizations, with a linear a^2 continuum extrapolation. The authors test the PCAC and pion-pole-dominance relations and quote continuum-limit values g_A = 1.245(28)(14), <r_A^2> = 0.339(48)(6) fm^2, g_P^* = 8.99(39)(49), and g_piNN = 13.25(67)(69).

Significance. If the results stand, this is one of the first, and by the authors' claim the first, continuum-limit determinations of these quantities obtained entirely at the physical pion mass, avoiding chiral extrapolation. The paper includes clear strengths: three physical-point ensembles with approximately equal volumes, nonperturbative renormalization, a broad set of sink-source separations, multiple parameterizations, and explicit two-state versus three-state comparisons. The physics outputs are directly relevant to neutrino-nucleus scattering modeling and to the muon-capture prediction for g_P^*. The main limitation is that the PCAC/PPD verification is partly entangled with a fit constraint that assumes PPD, and the continuum extrapolation rests on three lattice spacings with a linear a^2 ansatz; these issues need to be addressed or made explicit before the central claims can be fully accepted.

major comments (3)
  1. [Sec. 5, Sec. 6, Eq. (11)] Section 5 states: 'Since the pion pole dominance relation is satisfied at the continuum limit, we enforce the value of the pion-nucleon coupling constant g_piNN extracted from both form factors to be the same.' This makes the PPD verification in Sec. 6 and Fig. 4 partly tautological: the combined G_P/G_5 fit already imposes a common pole residue, so the near-pole behavior of the ratio rPPD in Eq. (11) is constrained by the fit rather than independently demonstrated. The PCAC ratio rPCAC also uses the same constrained G_P. To support the claimed restoration of PCAC/PPD, either repeat the analysis without the common-residue constraint and show that the two residues agree within errors, or explicitly state that the verification covers only the smooth parts of the form factors. This issue also bears on g_P^* and g_piNN, whose quoted values inherit the enforced pole residue.
  2. [Sec. 5, Sec. 7, Fig. 4] The continuum extrapolation is a linear function of a^2 over only three lattice spacings for every parameter, and the PCAC/PPD ratios show sizable cut-off effects that are then extrapolated to zero (Fig. 4). Three points do not allow a meaningful test of the linearity assumption or a robust assignment of truncation error. I ask for a quantitative statement of the cut-off systematic, for example a comparison with an a^4 term, a check of the slope using different pairs of ensembles, or a conservative error from the spread of the extrapolation. Without such a statement, the 'controlled continuum extrapolation' claim in Sec. 7 is stronger than the data support.
  3. [Sec. 5, Fig. 2] The final central values come from two-state fits, with the systematic error taken as the difference from three-state fits, but the three-state fits are stated to be stable only up to Q^2 ~ 0.5 GeV^2, while the z-expansion fits use two-state data up to 1 GeV^2. The excited-state contamination that the three-state fits are designed to remove is therefore unconstrained in the momentum region that contributes to the high-Q^2 shape of the form factors. Please show how the two-state versus three-state difference in the overlapping Q^2 region propagates to the continuum-limit g_A, <r_A^2>, g_P^*, and g_piNN, or restrict the central fits to Q^2 <= 0.5 GeV^2 where a direct three-state cross-check is available.
minor comments (5)
  1. [Sec. 4, Eq. (9)] The AIC weight in Eq. (9) uses N_dof as the penalty term; standard AIC uses the number of model parameters. Please clarify which information criterion from Refs. [14,15] is being used and how N_dof is defined in the presence of priors.
  2. [Sec. 3] The statement that neglected disconnected quark loop contributions are of order a^2 and vanish in the continuum limit should be supported by a one-line argument or an explicit reference; as written it can be misread as claiming that QCD disconnected contributions vanish in the continuum rather than only in the twisted-mass isovector combination considered here.
  3. [Fig. 2, caption] The caption for Fig. 2 should identify which bands represent the continuum limit and which represent finite lattice spacing, and what the red and yellow bands mean; the text says the red band comes from two-state fits and the yellow band is the systematic uncertainty, but the caption does not explain this.
  4. [Sec. 1] The claim to be 'the first calculation to use solely simulations at physical pion mass to take the continuum limit' would benefit from a qualifier such as 'to our knowledge' or a direct comparison with other physical-point efforts, to avoid an overclaim in a proceedings paper.
  5. [Throughout] Minor editorial issues: 'pseudo-scalar' and 'pseudoscalar' are used inconsistently, and the subscripts in Eq. (8) (A_i,j) would be clearer with a comma or by introducing a notation for the matrix element coefficients.

Circularity Check

1 steps flagged · score 5.0 of 10

PPD/PCAC validation is partly self-fulfilling: Sec. 5 enforces a common g_πNN by assuming PPD, so the Sec. 6/7 verification is not fully independent; g_A and r_A^2 remain independently extracted.

  1. self definitional [Sec. 5 (combined fit of G_P and G_5) and Sec. 6 (Eq. (11), Fig. 4)]
    "Since the pion pole dominance relation is satisfied at the continuum limit, we enforce the value of the pion-nucleon coupling constant gπNN extracted from both form factors to be the same."

    The rPPD ratio in Eq. (11) is presented as a test of PPD, but it uses G_P obtained from the Sec. 5 combined fit. That fit was constrained by explicitly assuming PPD: a single g_πNN pole residue was imposed on both G_P and G_5. Since equality of the two residues is one of the consequences of PPD together with PCAC, the fit removes a degree of freedom that the ratio would otherwise probe. The continuum-limit statement that PPD and PCAC are 'restored' is therefore partly a restatement of the assumption used to build the fit, not an independent check. The rPCAC ratio inherits the same constrained G_P. The primary g_A and <r_A^2> extraction is unaffected, and the g_P* and g_πNN numbers still depend on the lattice data, so the circularity is partial rather than total.

full rationale

The paper's main numerical extraction is a standard lattice QCD analysis: three physical-point twisted-mass ensembles, excited-state fits via spectral decomposition truncated at N_st=2 or 3, AIC model averaging, z-expansion fits, and an a^2 continuum extrapolation. None of these steps reduces to the claimed results; g_A and <r_A^2> are determined from correlation-function data without assuming the target values. The self-citations to Refs. [9] and [12] provide the ensembles, statistics, and renormalization constants; they are prior data inputs rather than the conclusion, so they are not load-bearing circularity. The one genuine circular element is the PPD/PCAC validation: Sec. 5 states that 'Since the pion pole dominance relation is satisfied at the continuum limit, we enforce the value of the pion-nucleon coupling constant gπNN extracted from both form factors to be the same.' This imports the relation being tested into the fit, so the subsequent rPPD and rPCAC agreement in Fig. 4 and the Sec. 7 statement that 'both relations are satisfied' are partially self-fulfilling. However, rPPD still depends on G_A and on the regular part of G_P, so the agreement is not forced to unity purely by construction. On balance the central claim has substantial independent lattice content, with a material but partial circularity in the symmetry-relation verification; score 5.

Assumptions & free parameters 4 free parameters · 4 assumptions · 0 invented entities

The central results depend on the standard lattice QCD framework, nonperturbative renormalization from Ref. [12], and the assumption that two-/three-state fits control excited-state contamination. The z-expansion coefficients and a^2 slopes are fitted to the same data used to report physical quantities; they are not independent inputs. No new entities are introduced. Most assumptions are standard for this subfield, but the enforcement of a common g_piNN before reporting PPD means one chiral relation is partly imposed rather than purely tested.

free parameters (4)
  • z-expansion coefficients a_k = not quoted; fitted to lattice data at Q^2 up to 1 GeV^2
    The Q^2 shape of G_A (and pole-subtracted G_P/G_5) is determined by these fitted coefficients with Gaussian priors; they set the reported radius and couplings.
  • linear a^2 cut-off slope coefficients = not quoted
    The continuum extrapolation assumes each fit parameter varies linearly with a^2 over three spacings; these slopes are fit to the three ensembles.
  • dipole mass m_A (dipole Ansatz alternative) = not quoted
    Used as a cross-check of the axial form factor Q^2 dependence; the radius follows from r^2 = 12/m^2.
  • excited-state energies and overlaps in multi-state fits = not quoted
    Nuisance parameters in the two- and three-state spectral fits; they affect the central values and the systematic difference between the two analysis choices.
assumptions (4)
  • domain assumption PCAC relation Eq. (4) is the correct chiral Ward identity for the lattice currents with nonperturbative renormalization.
    The paper uses PCAC as a consistency check; the relation is expected from QCD chiral symmetry, but its validity on the lattice relies on the renormalization scheme and current definitions.
  • domain assumption Pion pole dominance Eq. (5) holds near Q^2 = -m_pi^2; it is used to relate G_P and G_A and to constrain g_piNN.
    Invoked in Section 5 to justify a combined fit of G_P and G_5 with a common g_piNN, before being tested in Section 6.
  • domain assumption Disconnected quark loop contributions to isovector matrix elements are O(a^2) in twisted mass QCD and vanish in the continuum limit.
    Stated in Section 3 and justified by Ref. [13]; this justifies neglecting disconnected diagrams.
  • ad hoc to paper Gaussian priors on z-expansion coefficients suppress higher-order terms and ensure smooth convergence.
    This is an analysis model choice described in Section 5; it is not derived from QCD and affects the extracted radius and couplings.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Isovector axial and pseudoscalar form factors from twisted mass lattice QCD at the physical point." pith.science (2026). https://pith.science/paper/CQKT3U6B

@misc{pith2026250207583,
  author       = {Pith},
  title        = {Pith review of: Isovector axial and pseudoscalar form factors from twisted mass lattice QCD at the physical point},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/CQKT3U6B}},
  note         = {Machine review of arXiv:2502.07583}
}
abstract

We present the isovector axial, induced pseudoscalar, and pseudoscalar form factors of the nucleon using three twisted-mass fermion ensembles with degenerate up- and down-, strange-, and charm-quarks with masses tuned to their physical values (physical point). The three ensembles have lattice spacing $a$=0.08, 0.068, and 0.057 fm and approximately equal physical volume allowing for the continuum limit to be taken at the physical point. Excited-state contributions to the matrix elements are evaluated using several sink-source separations from 0.5 fm to 1.5 fm and multistate fits. We check the partially conserved axial-vector current (PCAC) hypothesis and the pion pole dominance (PPD) and show that in the continuum limit both relations are satisfied. We provide results at the continuum limit for the isovector nucleon axial charge, axial radius, pion-nucleon coupling constant, and for the induced pseudoscalar form factor at the muon capture point.

Figures

Figures reproduced from arXiv: 2502.07583 by the authors.

Figure 1
Figure 1. The ratio of Eq. (10) for the ensemble cC211.060.80 and for three cases, namely that yielding 𝑔𝐴 (top), that yielding a linear combination of the axial and induced pseudoscalar form factor at the first non-zero momentum transfer (middle) and that yielding the pseudoscalar form factor at the first non-zero momentum transfer (bottom). In the left column we plot the ratio versus the insertion time, in the middle we plo… view at source ↗
Figure 2
Figure 2. Left: The axial form factor 𝐺𝐴(𝑄 2 ) determined from two-state fits (red band) with systematic uncertainty from the difference with three-state fits (yellow band). Results at finite lattice spacing shown in blue, orange and green bands. Right: Results for the axial charge 𝑔𝐴 and radius ⟨𝑟 2 𝐴 ⟩ obtained using different analysis approaches. convergence of the z-expansion at 𝑘max = 3 and verify stability with respect … view at source ↗
Figure 3
Figure 3. Left: The induced pseudoscalar form factor 𝐺𝑃 (𝑄 2 ) and right: the pseudoscalar form factor 𝐺˜ 5 (𝑄 2 ) at finite lattice spacing (blue, orange and green bands) and in the continuum limit (red band). Results are obtained using two-state fits and z-expansion of order 3. The inner panels show the region near the pion pole enlarged. 0.0 0.2 0.4 0.6 0.8 1.0 Q 2 [GeV 2 ] 0.4 0.6 0.8 1.0 1.2 rPCAC(Q 2 ) cB211.072.64 cC21… view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Left: Ratio testing the PCAC relation. Right: Ratio testing pion pole dominance. Results at finite lattice spacing shown in blue, orange and green bands and in the continuum limit with red band. The yellow band includes systematic uncertainties from excited states. In …

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. The Nucleon Axial Form Factor from Elementary Target Data

    hep-ex 2025-12 conditional novelty 6.0 of 10

    The nucleon axial form factor from hydrogen and lattice-QCD data falls more slowly with Q² than deuterium-based fits, indicating deuterium extractions are biased low.

Reference graph

Works this paper leans on

16 extracted references · 6 canonical work pages · cited by 1 Pith paper

  1. [9]

    Alexandrou, S

    Extended Twisted Mass collaboration, C. Alexandrou, S. Bacchio, M. Constantinou, J. Finkenrath, R. Frezzotti et al.,Nucleon axial and pseudoscalar form factors using twisted-mass fermion ensembles at the physical point,Phys. Rev. D109(2024) 034503 [2309.05774]

    arXiv 2024

  2. [1]

    Cai et al.,Measurement of the axial vector form factor from antineutrino–proton scattering, Nature 614 (2023) 48

    MINERvAcollaboration, T. Cai et al.,Measurement of the axial vector form factor from antineutrino–proton scattering, Nature 614 (2023) 48

    work page 2023

  3. [2]

    UCNAcollaboration, M. A. P. Brown et al.,New result for the neutron𝛽-asymmetry parameter𝐴0 from UCNA, Phys. Rev. C97 (2018) 035505 [1712.00884]

    arXiv 2018

  4. [3]

    Darius et al.,Measurement of the Electron-Antineutrino Angular Correlation in Neutron 𝛽 Decay, Phys

    G. Darius et al.,Measurement of the Electron-Antineutrino Angular Correlation in Neutron 𝛽 Decay, Phys. Rev. Lett.119 (2017) 042502

    work page 2017

  5. [4]

    UCNAcollaboration, M. P. Mendenhall et al.,Precision measurement of the neutron 𝛽-decay asymmetry,Phys. Rev. C87(2013) 032501 [1210.7048]

    arXiv 2013

  6. [5]

    D. Mund, B. Maerkisch, M. Deissenroth, J. Krempel, M. Schumann et al.,Determination of the Weak Axial Vector Coupling from a Measurement of the Beta-Asymmetry Parameter A in Neutron Beta Decay, Phys. Rev. Lett.110 (2013) 172502 [1204.0013]

    work page Pith review arXiv 2013

  7. [6]

    Liu, J.-M

    K.-F. Liu, J.-M. Wu, S.-J. Dong and W. Wilcox,Lattice calculation of nucleon axial form-factors, Nucl. Phys. B Proc. Suppl.20(1991) 467

    work page 1991

  8. [7]

    K. F. Liu, S. J. Dong, T. Draper, J. M. Wu and W. Wilcox,Nucleon axial form-factor from lattice QCD, Phys. Rev. D49(1994) 4755 [hep-lat/9305025]

    work page Pith review arXiv 1994

Show all 16 references
  1. [8]

    Alexandrou, G

    C. Alexandrou, G. Koutsou, T. Leontiou, J. W. Negele and A. Tsapalis,Axial Nucleon and Nucleon to Delta form factors and the Goldberger-Treiman Relations from Lattice QCD, Phys. Rev. D76(2007) 094511 [0706.3011], [Erratum:Phys. Rev. D80 (2009) 099901]

    2007 arXiv

  2. [10]

    Finkenrath et al.,Twisted mass gauge ensembles at physical values of the light, strange and charm quark masses, PoS LATTICE2021(2022) 284 [2201.02551]

    J. Finkenrath et al.,Twisted mass gauge ensembles at physical values of the light, strange and charm quark masses, PoS LATTICE2021(2022) 284 [2201.02551]

    2022 arXiv

  3. [11]

    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]

    2018 arXiv

  4. [12]

    Extended Twisted Mass collaboration, C. Alexandrou et al.,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

  5. [13]

    Frezzotti and G

    R. Frezzotti and G. C. Rossi,Chirally improving Wilson fermions. 1. O(a) improvement, JHEP 08 (2004) 007 [hep-lat/0306014]

    2004 arXiv

  6. [14]

    W. I. Jay and E. T. Neil,Bayesian model averaging for analysis of lattice field theory results, Phys. Rev. D103 (2021) 114502 [2008.01069]. 9 Isovector axial FFs Giannis Koutsou

    2021 arXiv

  7. [15]

    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]

    2024 arXiv

  8. [16]

    Alexandrou, S

    C. Alexandrou, S. Bacchio, C. Iona, G. Koutsou, Y. Li et al.,Nucleon axial, tensor, and scalar charges and𝜎-terms in lattice QCD, PoS LATTICE2024(2025) 316. 10

    2025

Pith tools

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