Pith. sign in

REVIEW 3 major objections 6 minor 1 cited by

Studies of nucleon isovector structure with the PACS10 superfine lattice

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

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

Pith's one-line read The paper claims that a time-derivative correlator method removes the leading pion-nucleon contamination from the induced pseudoscalar form factor, and that on the PACS10 superfine lattice the nucleon axial-vector, induced pseudoscalar…

desk verdict New superfine PACS10 data are a genuine step, but the g*_P and g_piNN results rest on an unshown 'complete removal' of piN contamination, so treat this as a progress report, not a definitive calculation. read the letter →

arxiv 2411.16784 v3 pith:PRUXGASA submitted 2024-11-25 hep-lat

Ryutaro Tsuji , Yasumichi Aoki , Ken-Ichi Ishikawa , Yoshinobu Kuramashi , Shoichi Sasaki , Kohei Sato , Eigo Shintani , Hiromasa Watanabe
show 1 more author
Takeshi Yamazaki
This is my paper · ORCID
classification hep-lat PACS 12.38.Gc14.20.Dh
keywords latticeQCDnucleonaxial-vectorcouplinginducedpseudoscalarpion-nucleonPACS10ensemblesexcited-statecontaminationtime-derivativecorrelatormethodphysicalpoint
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

This paper is a lattice QCD determination of three couplings that control the weak axial structure of the nucleon: the axial-vector coupling $g_A$, the induced pseudoscalar coupling $g_P^*$, and the pion-nucleon coupling $g_{\pi NN}$. The authors use the PACS10 ensembles at three lattice spacings (0.09, 0.06, and 0.04 fm) on a volume with large spatial extent of about 10 fm and physical quark masses, and they report that all three couplings are consistent across spacings and agree with experiment and with other lattice QCD results. The central technical claim is that a newly proposed time-derivative correlator method completely removes the leading pion-nucleon excited-state contamination from the induced pseudoscalar form factor $F_P$, which has long blocked accurate values of $g_P^*$ and $g_{\pi NN}$. The results for the superfine lattice are explicitly preliminary, and no continuum extrapolation is attempted; instead the paper argues that the absence of strong lattice-spacing dependence makes the three-spacing comparison a controlled start toward one. If the claim holds, precision lattice values of these couplings could feed neutrino oscillation experiments and the neutron lifetime puzzle.

What carries the argument

The load-bearing object is the time-derivative correlator method for $F_P$ (Ref. [21], in preparation), a correlator-based technique the paper says removes the leading pion-nucleon excited-state contribution once the exponential smearing parameters are tuned; the paper does not derive it, so the method itself carries the $g_P^*$ and $g_{\pi NN}$ results. Around it stand the PACS10 ensembles with the six stout-smeared $O(a)$-improved Wilson-clover quark action and Iwasaki gauge action at $\beta=1.82$, $2.00$, $2.20$; all-mode averaging for statistics; the Schr\"odinger functional renormalization factor $Z_A$; and the $z$-expansion parametrization of $(q^2+m_\pi^2)F_P(q^2)$, which factors out the pion pole before the couplings are read off.

What would settle it

Apply the time-derivative correlator method to $F_P$ on the superfine ensemble at several source-sink separations beyond the current $t_{\mathrm{sep}}/a = 20, 29$ and check whether $g_P^*$ and $g_{\pi NN}$ stay constant; a visible drift with $t_{\mathrm{sep}}$, or an analysis with a variational basis that disagrees with the new-method result, would show that the leading $\pi N$ contribution was not completely removed.

Watch

Extended reading notes

Core claim

On the PACS10 ensembles, the authors compute the nucleon two- and three-point functions with exponentially smeared interpolating operators and all-mode averaging, extract $g_A$ from the axial form factor at zero momentum transfer, and extract $g_P^*$ and $g_{\pi NN}$ from $F_P(q^2)$ at seven momentum transfers using a $z$-expansion of $(q^2+m_\pi^2)F_P(q^2)$ that factors out the pion pole. Their main stated result is that, with careful tuning of the smearing parameters, the new time-derivative correlator method (Ref. [21], in preparation) 'completely removes the leading $\pi N$ contribution' from the $F_P$ analysis, so that the two pseudoscalar couplings can be determined accurately. At the three lattice spacings the nucleon dispersion relation deviates from the continuum relation by at most 1.1%, and the three couplings show no strong dependence on the lattice spacing; $g_A$, $g_P^*$, and $g_{\pi NN}$ agree with the experimental values and with other lattice QCD results within the quoted uncertainties. The superfine data are still preliminary, so the paper stops short of a continuum extrapolation.

Load-bearing premise

The load-bearing premise is that the unpublished time-derivative correlator method removes the entire leading pion-nucleon excited-state contamination from the $F_P$ data, rather than merely reducing it, so that the residual time dependence is negligible.

Editorial extensions

If this is right

  • A controlled continuum extrapolation of $g_A$, $g_P^*$, and $g_{\pi NN}$ becomes a statistics-limited exercise on the existing coarse, fine, and superfine ensembles rather than a search for new analysis methods.
  • The two least well-known axial couplings, $g_P^*$ and $g_{\pi NN}$, can be produced at physical quark mass with the leading excited-state contamination under control, removing a known source of lattice systematic error.
  • Reliable values of these couplings can be used as inputs to neutrino-nucleus cross-section calculations for long-baseline oscillation experiments, where the axial form factor is a significant uncertainty.
  • The $g_A$ result at all three spacings serves as a benchmark that the same ensembles and analysis reproduce a well-measured experimental quantity, supporting the credibility of the less well-known pseudoscalar results.
  • Combining the three spacings with proper systematic errors should allow a continuum limit with controlled discretization uncertainty, since the observed spacing dependence is already smaller than current statistical errors.

Reading between the lines

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

  • If the time-derivative method really removes all leading $\pi N$ contamination, the same idea could be applied to other nucleon matrix elements (scalar, tensor, vector) whose plateaus suffer from excited-state contamination; the paper does not discuss this extension.
  • Because the decisive method is cited as 'in preparation' and is not validated here, the stability of the $g_P^*$ and $g_{\pi NN}$ results cannot be independently checked until Ref. [21] appears; that single reference is the fastest possible falsifier.
  • The apparent absence of lattice-spacing dependence between 0.09 and 0.04 fm suggests the remaining discretization errors are below the statistical noise; if so, adding statistics on the superfine lattice will matter more than further action improvement.
  • A simple testable consequence is that the new method should make the $F_P$ plateau strictly independent of the source-sink separation; any residual $t_{\mathrm{sep}}$ drift would indicate the leading contamination is suppressed rather than completely removed.
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 / 6 minor

Summary. This manuscript reports lattice QCD calculations of the isovector nucleon axial-vector coupling g_A, the induced pseudoscalar coupling g*_P, and the pion-nucleon coupling g_piNN using PACS10 ensembles at three lattice spacings (0.09, 0.06, 0.04 fm) with a physical pion mass and a spatial extent of about 10 fm. The axial form factor is extracted with a plateau method, while the induced pseudoscalar form factor F_P is analyzed with a z-expansion including the pion-pole factor, using a new time-derivative correlator method (Ref. [21], in preparation) that is claimed to completely remove the leading piN excited-state contamination. Results for the three couplings are compared with experimental values and with other lattice QCD determinations. The superfine results are explicitly labeled preliminary and are shown with statistical errors only; no continuum extrapolation is performed.

Significance. If the new F_P method performs as claimed, the paper would demonstrate that the long-standing piN excited-state contamination problem for the induced pseudoscalar form factor can be overcome at the physical point, and the observed consistency across three lattice spacings would support the prospect of a controlled continuum extrapolation of g_A, g*_P, and g_piNN. The work has notable strengths: physical-point ensembles with large volume, the use of AMA, nonperturbative O(a) improvement, Schrödinger functional renormalization, and comparison with multiple independent lattice and experimental determinations. The central caveat is that the decisive F_P technique is not documented or validated in this manuscript, so the two pseudoscalar couplings currently cannot be independently checked by the reader. The paper is internally consistent and appropriately hedged about the preliminary nature of the superfine results.

major comments (3)
  1. [Section 4.2 and Summary (Ref. [21])] The claim that the new time-derivative correlator method (Ref. [21], 'in preparation') completely removes the leading piN contribution from the F_P data is the load-bearing step for the g*_P and g_piNN results. The manuscript gives no derivation, algorithmic description, or validation of this method, and no estimate of the residual contamination after the removal. Since F_P(q^2) is the sole input to the z-expansion that determines g*_P and g_piNN, any residual piN contamination biases these couplings directly. The authors should either include the method and demonstrate its effectiveness (for example on synthetic correlators with known piN contamination, or by a comparison with an independent extraction technique), or replace 'completely removes' with a quantitatively bounded residual-estimate statement. As written, the two pseudoscalar couplings cannot be independently assessed by the reader.
  2. [Section 3 (renormalization)] The renormalization constants Z_A and Z_V are said to be determined by the Schrödinger functional method, with a pointer to Appendix E of Ref. [20]. It is not stated whether the Z_A value used for the new superfine (beta=2.20) ensemble was computed directly at beta=2.20, interpolated from lower beta, or assumed unchanged, nor how the associated uncertainty is propagated. Because the reported g_A is proportional to Z_A, this point must be documented before the g_A agreement with experiment can be evaluated. Please report the Z_A values and their scale/beta dependence explicitly.
  3. [Section 4.2 and Summary (lattice-spacing dependence)] The statement that 'no strong dependence on the lattice spacing is observed' for the three couplings is not supported by a controlled continuum extrapolation. The superfine results in Fig. 2 carry only statistical errors, are shown at two separate tsep values without a combined analysis, and no correlated fit in the lattice spacing is performed. The text is mostly careful in calling the results preliminary, but the Summary phrase 'there is no strong dependence on the lattice spacing in our results' should be qualified to say that this is a qualitative consistency check within current large uncertainties, not an established continuum behavior.
minor comments (6)
  1. [Abstract] The phrase 'Combining the results obtained from the all of our coarse, fine and superfine lattices' is ungrammatical and overstates the analysis, since no combined fit is actually presented; please rephrase to refer to a comparison of results from the three ensembles.
  2. [Figure 2 caption] The caption says 'Two open triangle symbols represent the preliminary results... with tsep = 0.8 fm (the upper value...) and tsep = 1.2 fm (the lower value...)', but it is not immediately clear in the figure which symbol corresponds to which tsep; please add explicit labels to the figure.
  3. [Section 4.2] The z-expansion fit used for F_P is not specified in terms of truncation order, number of coefficients, priors, or treatment of correlations; a sentence giving these details or an explicit pointer to the exact implementation in Ref. [20] would improve reproducibility.
  4. [Section 2] Since the new time-derivative correlator method is central to the paper, please add a footnote or short appendix outlining its key idea even if the full derivation is deferred to Ref. [21].
  5. [Section 5] The phrase 'Needless to say' is informal for a journal article and should be replaced with a more neutral formulation.
  6. [References] Reference [9] appears to cite a preprint from 2016; please check whether a published version is available and update the citation if so.

Circularity Check

1 steps flagged · score 4.0 of 10

The g*_P and g_piNN results rest on the in-preparation self-citation [21] for complete piN removal, while g_A and the z-expansion analysis are otherwise independently checked.

  1. self citation load bearing [Section 4.2, supported by Section 2 and Ref. [21]]
    "We emphasize here that the excited-state contamination has been a long-standing obstacle to the accurate calculation of the two couplings associated with the induced pseudoscalar form factor F_P. However, in our calculations where careful tuning of the smearing parameters for the nucleon ground state was performed, we have succeeded in completely removing the leading piN contribution from the analysis of the F_P data by our newly proposed method [21]."

    The central numerical results for g*_P and g_piNN are obtained exclusively from F_P(q^2) through the z-expansion (Eq. (1) and Section 4.2). The only support offered for the decisive premise that the F_P data are free of the leading piN excited-state contamination is Ref. [21], an in-preparation paper by the same PACS authors. Section 2 likewise states only that 'our new method using time-derivative correlator is employed for F_P [21]' without derivation or internal validation. The paper does not exhibit the method, show that the time-derivative procedure eliminates the entire leading piN term, or cross-check it on synthetic data. Thus the accuracy claim for F_P reduces to an unverified self-citation rather than an independently substantiated input.

full rationale

The remainder of the paper is not circular. The axial-vector coupling g_A is extracted from F_A by the standard plateau method and compared with independent experimental and FLAG/lattice results. The renormalization constant Z_A is taken from the published Ref. [20] (Schrodinger functional determination), not from a fit to the target quantities. The z-expansion is a standard model-independent parameterization whose original reference is Hill and Paz [22], and no continuum-limit extrapolation is claimed. The statement that no strong lattice-spacing dependence is observed is a direct comparison of the three PACS10 ensembles, not a fitted parameter renamed as a prediction. Therefore the only significant circularity concern is the reliance on the in-preparation self-citation [21] for the complete removal of the leading piN contamination in F_P. Because that reliance is load-bearing for g*_P and g_piNN but the rest of the analysis has independent content, a score of 4 is appropriate rather than a higher score reserved for derivations that reduce to their inputs by construction.

Assumptions & free parameters 2 free parameters · 6 assumptions · 0 invented entities

The calculation contributes data (new superfine ensemble) but relies on the collaboration's previous papers for ensemble generation, renormalization, and the new F_P method, plus standard form-factor machinery. No new physical entities are introduced; the main ledger items are the tuned smearing parameters, the z-expansion fit coefficients, and the unshown time-derivative method.

free parameters (2)
  • Exponential smearing parameters (A, B) = A=1.2, B=0.16 (128^4); A=1.2, B=0.11 (160^4); A=1.2, B=0.07 (256^4)
    Chosen by hand per ensemble to maximize ground-state overlap of the nucleon interpolating operator; they affect the quality of the plateaus and thus the extracted couplings, but they are not fitted to the target couplings themselves.
  • Coefficients of the z-expansion for (q^2 + m_pi^2) F_P(q^2) = not given in the paper
    The induced pseudoscalar form factor at 7 momentum transfers is fitted with a z-expansion truncated at some order; the fitted coefficients determine g*_P and g_piNN. The values and truncation order are not reported in this proceedings paper.
assumptions (6)
  • domain assumption Lattice QCD with the six stout-smeared O(a) improved Wilson-clover quark action and Iwasaki gauge action provides a valid discretization of QCD, with lattice artifacts controlled by O(a) improvement and the continuum limit in principle reachable.
    The entire calculation rests on this discretization; the paper checks the dispersion relation to bound cut-off effects but does not take a continuum limit for the couplings.
  • standard math The asymptotic ratio in Eqs. (6)-(7) isolates the ground-state matrix element in the plateau region.
    Standard spectral decomposition of two- and three-point functions; the paper uses the plateau method for F_A and the new method for F_P.
  • standard math The Lorentz decomposition of the axial matrix element in Eq. (2) is complete for the nucleon axial-vector current.
    Relies on Lorentz covariance and parity; standard in the lattice and phenomenology literature.
  • domain assumption Renormalization factors Z_A and Z_V determined by the Schrodinger functional method in Ref. [20] apply to the superfine ensemble.
    The paper states Z_A and Z_V are determined by the SF method (Appendix E of Ref. [20]) and uses them at all three lattice spacings, including the new superfine one, without a re-determination shown here.
  • standard math The z-expansion with the pion-pole factor (q^2 + m_pi^2) is a model-independent parameterization that can be extrapolated to the pion pole for g_piNN.
    The z-expansion is a conformal-mapping parameterization; the pole factor removes the dominant pion-pole singularity. The reliability depends on the number of fit points and truncation order, which are not fully specified.
  • ad hoc to paper The new time-derivative correlator method (Ref. [21], in preparation) removes the leading piN excited-state contamination from F_P.
    This is the paper's central methodological premise for the g*_P and g_piNN results; it is not derived or validated in this paper, only cited as in preparation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Studies of nucleon isovector structure with the PACS10 superfine lattice." pith.science (2026). https://pith.science/paper/PRUXGASA

@misc{pith2026241116784,
  author       = {Pith},
  title        = {Pith review of: Studies of nucleon isovector structure with the PACS10 superfine lattice},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/PRUXGASA}},
  note         = {Machine review of arXiv:2411.16784}
}
abstract

We present the results for the nucleon axial-vector, induced pseudoscalar and pion-nucleon couplings obtained from 2+1 flavor lattice QCD at the physical point with a large spatial extent of about 10 fm. Our calculations are performed with the PACS10 gauge configurations generated by the PACS Collaboration with the six stout-smeared $O(a)$ improved Wilson-clover quark action and Iwasaki gauge action at $\beta$ = 1.82, 2.00 and 2.20 corresponding to lattice spacings of 0.09 fm (coarse), 0.06 fm (fine) and 0.04 fm (superfine), respectively. We first evaluate the value of the nucleon axial-vector coupling. In addition, the induced pseudoscalar and pion-nucleon couplings from the induced pseudoscalar form factor are also investigated. Combining the results obtained from the all of our coarse, fine and superfine lattices, we finally discuss the systematic uncertainties in our calculation based on the comparison with both of the experimental values and lattice QCD results provided by the other collaborations.

Figures

Figures reproduced from arXiv: 2411.16784 by the authors.

Figure 1
Figure 1. Check of the dispersion relation for the nucleon with an improved method described in Ref. [20]. A dashed line represents the relativistic continuum dispersion relation, while red and blue dotted lines are given by the linear fit of the coarse and fine data set, respectively. statistical uncertainties. The continuum-limit extrapolation requires a further calculation for the superfine lattice, and all of the systemat… view at source ↗
Figure 2
Figure 2. Summary for our current status, the experimental values and other lattice QCD results for the axial-vector coupling (left), induced pseudoscalar coupling (middle) and pion-nucleon coupling (right). The inner error bars represent the statistical errors, while the outer error bars are evaluated by both the statistical and systematic errors added in quadrature. Uncertainties stemming from the excited-state contaminatio… view at source ↗

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. A proposal for removing $\pi N$-state contamination from the nucleon induced pseudoscalar form factor in lattice QCD

    hep-lat 2025-01 conditional novelty 7.0 of 10

    A time-derivative subtraction of axial-vector correlators removes leading pion-nucleon contamination from the nucleon induced pseudoscalar form factor, yielding plateau values that match the pion-pole-dominance model.

Reference graph

Works this paper leans on

25 extracted references · 4 canonical work pages · cited by 1 Pith paper

  1. [21]

    Yamazaki (PACS) (2024), in preparation

    S.Sasaki,Y.Aoki,K.-I.Ishikawa,Y.Kuramashi,K.Sato,E.Shintani,R.Tsuji,H.Watanabe, and T. Yamazaki (PACS) (2024), in preparation

    work page 2024

  2. [20]

    Yamazaki (PACS), Phys

    R.Tsuji,Y.Aoki,K.-I.Ishikawa,Y.Kuramashi,S.Sasaki,K.Sato,E.Shintani,H.Watanabe, and T. Yamazaki (PACS), Phys. Rev. D109, 094505 (2024),2311.10345

    arXiv 2024

  3. [1]

    Abe et al

    K. Abe et al. (T2K), Nucl. Instrum. Meth. A659, 106 (2011),1106.1238

    arXiv 2011

  4. [2]

    Abe et al

    K. Abe et al. (Hyper-Kamiokande Proto-), PTEP2015, 053C02 (2015),1502.05199

    arXiv 2015

  5. [3]

    Abe et al

    K. Abe et al. (T2K), Nature580, 339 (2020), 1910.03887, [Erratum: Nature 583, E16 (2020)]

    arXiv 2020

  6. [4]

    M. A. Acero et al. (NOvA), Phys. Rev. Lett.123, 151803 (2019),1906.04907

    arXiv 2019

  7. [5]

    Abi et al

    B. Abi et al. (DUNE) (2020),2002.03005

    arXiv 2020

  8. [6]

    120,202002(2018), 1802.01804

    A.Czarnecki,W.J.Marciano,andA.Sirlin,Phys.Rev.Lett. 120,202002(2018), 1802.01804

    arXiv 2018

Show all 25 references
  1. [7]

    R. L. Workman et al. (Particle Data Group), PTEP2022, 083C01 (2022)

    2022

  2. [8]

    V. A. Andreev et al. (MuCap), Phys. Rev. C91, 055502 (2015),1502.00913

    2015 arXiv

  3. [9]

    V. A. Babenko and N. M. Petrov (2016),1604.02912

    2016 arXiv

  4. [10]

    Limkaisang, K

    V. Limkaisang, K. Harada, J. Nagata, H. Yoshino, Y. Yoshino, M. Shoji, and M. Matsuda, Prog. Theor. Phys.105, 233 (2001)

    2001

  5. [11]

    Aoki et al

    Y. Aoki et al. (Flavour Lattice Averaging Group (FLAG)), Eur. Phys. J. C82, 869 (2022), 2111.09849, and all relevant references therein. 8 Studies of nucleon isovector structure with the PACS10 superfine lattice Ryutaro Tsuji

    2022 arXiv

  6. [12]

    G.S.Bali,S.Collins,S.Heybrock,M.Löffler,R.Rödl,W.Söldner,andS.Weishäupl(2023), 2305.04717

    2023 arXiv

  7. [13]

    D.Djukanovic, G.vonHippel, J.Koponen, H.B.Meyer, K.Ottnad, T.Schulz, andH.Wittig, Phys. Rev. D106, 074503 (2022),2207.03440

    2022 arXiv

  8. [14]

    Tsuji, N

    R. Tsuji, N. Tsukamoto, Y. Aoki, K.-I. Ishikawa, Y. Kuramashi, S. Sasaki, E. Shintani, and T. Yamazaki (PACS), Phys. Rev. D106, 094505 (2022),2207.11914

    2022 arXiv

  9. [15]

    S. Park, R. Gupta, B. Yoon, S. Mondal, T. Bhattacharya, Y.-C. Jang, B. Joó, and F. Winter (Nucleon Matrix Elements (NME)), Phys. Rev. D105, 054505 (2022),2103.05599

    2022 arXiv

  10. [16]

    Pohl et al., Nature466, 213 (2010)

    R. Pohl et al., Nature466, 213 (2010)

    2010

  11. [17]

    Tomalak, R

    O. Tomalak, R. Gupta, and T. Bhattacharya, Phys. Rev. D108, 074514 (2023),2307.14920

    2023 arXiv

  12. [18]

    R.Gupta, PoS LATTICE2023, 124(2024), 2401.16614, andallrelevantreferencestherein

    2024 arXiv

  13. [19]

    Shintani, K.-I

    E. Shintani, K.-I. Ishikawa, Y. Kuramashi, S. Sasaki, and T. Yamazaki, Phys. Rev. D99, 014510 (2019),1811.07292, [Erratum: Phys. Rev.D 102, 019902 (2020)]

    2019 arXiv

  14. [22]

    R. J. Hill and G. Paz, Phys. Rev. D82, 113005 (2010),1008.4619

    2010 arXiv

  15. [23]

    T. Blum, T. Izubuchi, and E. Shintani, Phys. Rev. D88, 094503 (2013),1208.4349

    2013 arXiv

  16. [24]

    E.Shintani,R.Arthur,T.Blum,T.Izubuchi,C.Jung,andC.Lehner,Phys.Rev.D 91,114511 (2015), 1402.0244

    2015 arXiv

  17. [25]

    von Hippel, T

    G. von Hippel, T. D. Rae, E. Shintani, and H. Wittig, Nucl. Phys. B914, 138 (2017), 1605.00564. 9

    2017 arXiv

Pith tools

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