REVIEW 3 major objections 4 minor 1 cited by
Determination of the pseudoscalar decay constant from SU(2) with two fundamental flavors
T0 review · 3 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read For SU(2) gauge theory with two fundamental flavors, the paper reports a continuum-extrapolated pseudoscalar decay constant $w_0 f_{\rm PS}=0.1436(19)$ at reference mass $(w_0 M_{\rm PS})^2=1.09(2)$, obtained with non-perturbatively…
desk verdict First continuum estimate of f_PS in SU(2) with two fundamental flavors, but the central value rests on only three lattice spacings and the paper itself admits a finer point is needed. 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 that carries the result is the mixed-action combination of two non-perturbative O(a)-improvement mechanisms. The sea sector uses Wilson fermions with non-perturbatively tuned exponential clover improvement, defined through the Dirac operator in Eq. (4), which suppresses discretization effects and improves stability. The valence sector adds a chirally rotated twisted-mass term and is tuned to maximal twist, meaning the PCAC mass is set to zero and the valence pseudoscalar mass is matched to the sea, so that the pseudoscalar decay constant is automatically O(a)-improved and, crucially, renormalization-free: $f_{\rm PS}^R = 2\mu_0 \langle 0|P|\pi\rangle / M_{\rm PS}^2$ depends only on the bare twisted mass, the pseudoscalar mass, and the matrix element. This removes a renormalization-constant uncertainty, allowing the precise scale setting in units of $w_0$ that leads to the continuum value.
What would settle it
Generate one additional ensemble at a finer lattice spacing, for example $\beta = 2.4$ or a smaller value of $w_0/a$, matched to the same reference mass $(w_0 M_{\rm PS})^2 = 1.09(2)$, and check whether its $w_0 f_{\rm PS}$ falls on the O($a^{2}$) continuum extrapolation through the three existing points within errors. If a residual O(a) term is present, the new point will deviate from the quoted $0.1436(19)$ by more than the quoted uncertainty.
Extended reading notes
Core claim
On its own terms, the paper claims that the continuum limit of the pseudoscalar decay constant can already be taken with high precision and small discretization effects for SU(2) with two fundamental flavors. At the reference mass $(w_0 M_{\rm PS})^2 = 1.09(2)$, the continuum-extrapolated value is $w_0 f_{\rm PS} = 0.1436(19)$, with the three available lattice spacings behaving in an O($a^{2}$)-compatible way; the result from a linear fit is stated to be only in slight tension with this continuum value. This is the first study to quantify the size of discretization effects with exponential clover improvement for the SU(2) gauge group, finding effects below 10% in the explored lattice-spacing range, compared with about 30% in the earlier unimproved-Wilson study. The authors are careful to call the result preliminary: one more ensemble at a finer lattice spacing is required to certify that the O(a) improvement is exact enough that only O($a^{2}$) corrections remain.
Load-bearing premise
The continuum extrapolation assumes the non-perturbative O(a) improvement is exact enough that the remaining lattice-spacing dependence is O($a^{2}$), an assumption the paper itself flags as not yet confirmed because only three lattice spacings are available.
Editorial extensions
If this is right
- The quoted continuum value gives a scale-setting anchor for the SU(2) composite-Higgs theory at a fixed reference mass, removing one systematic uncertainty from predictions of the spectrum.
- The renormalization-free formula at maximal twist means future ensembles can improve the decay constant without requiring a separately computed renormalization constant, reducing a major source of error.
- With the same action and analysis, adding one finer lattice spacing is enough to turn the preliminary continuum value into a certified O(a^2)-extrapolated result.
- A future chiral extrapolation to $(w_0 M_{\rm PS})^2 \to 0$ can convert this fixed-mass value into the chiral-limit decay constant needed for composite-Higgs parameter constraints.
Reading between the lines
- The same mixed-action, maximal-twist strategy should transfer directly to other composite-Higgs candidate gauge theories, since the renormalization-free property of $f_{\rm PS}$ and the automatic O(a) improvement do not depend on the specific gauge group.
- If a finer-spacing ensemble confirms O(a^2) scaling, the residual slight tension in the current linear fit would most naturally indicate a small next-order O(a^4) contribution, and the quoted $0.1436(19)$ would likely move by less than the present error.
- With GPU acceleration now available for the simulation code, producing the additional chiral ensembles and finer spacings described as future work is computationally feasible, so a percent-level continuum result for the full chiral-limit decay constant appears within reach.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports a preliminary lattice determination of the pseudoscalar decay constant f_PS in SU(2) gauge theory with two fundamental Dirac flavors, a candidate composite-Higgs theory. Ensembles are generated with non-perturbatively O(a)-improved (exponential clover) Wilson fermions in the sea and measured with maximally twisted valence quarks, so that f_PS is obtained from a renormalization-free expression involving the bare twisted mass, the pseudoscalar mass, and a matrix element. Using ensembles at beta = 2.15, 2.2, and 2.3, the authors interpolate af_PS and w0/a to a common reference mass (w0 MPS)^2 = 1.09(2) and take a continuum limit, quoting w0 f_PS = 0.1436(19). They note that an additional finer lattice spacing is needed to confirm that O(a) effects are eliminated, and they frame the work as preliminary.
Significance. Assuming the continuum extrapolation is confirmed with an additional finer lattice spacing, this would be the first continuum determination of the pseudoscalar decay constant for SU(2) with two fundamental flavors from an O(a)-improved mixed-action setup, and it would provide a useful benchmark for composite-Higgs model building. The paper has a genuine methodological strength: Eq. (7) avoids renormalization constants at maximal twist, so no fitted parameter enters the final formula, and the statistical errors are small. The reported reduction of discretization effects relative to the earlier unimproved study (below 10% versus roughly 30%) is a valuable demonstration of the exponential clover and mixed-action strategy. The main caveat is that the central continuum number is not yet established beyond a statistical-only extrapolation with three lattice spacings.
major comments (3)
- [Section 4 (continuum extrapolation) and concluding paragraph] The quoted continuum value w0 f_PS = 0.1436(19) at reference mass (w0 MPS)^2 = 1.09(2) is obtained from three lattice spacings only (beta = 2.15, 2.2, 2.3), with one interpolated point per beta. The paper itself states that one more point with finer lattice spacing is needed to confirm that O(a) effects are eliminated; with three points it is not possible to discriminate between the expected O(a^2) scaling and residual O(a) contamination. The error quoted is statistical only. This is load-bearing because the continuum value is the primary result; the authors should quantify a discretization systematic (for example, by comparing O(a) and O(a^2) fits, adding a term linear in a, or quoting a conservative uncertainty) or clearly label the number as preliminary rather than as the continuum result.
- [Section 4] The 'linear fit' whose result is in slight tension with 0.1436(19) is not defined: it is not stated whether the fit is linear in a or in a^2, what its continuum intercept is, or what its chi^2 per degree of freedom is. If the fit is linear in a, the 'slight tension' is exactly the symptom of the residual O(a) effect that the authors say they cannot exclude; if it is linear in a^2, the comparison is not informative about O(a) contamination. Please report the fit form, the fitted continuum value, and the fit quality so that the reader can judge the robustness of the central number.
- [Section 2.3, Eqs. (7)-(8)] The renormalization-free formula assumes exact maximal twist (m_PCAC = 0). The manuscript does not report the residual PCAC masses or the size of the correction from Eq. (8), even though any mistuning enters the decay constant directly. Since the final precision is about 1.3%, the tuning uncertainty should be quantified or the Z_A correction should be applied; otherwise the quoted error does not include a potentially relevant systematic effect.
minor comments (4)
- [Table 1] The row for beta = 2.3 reads '364'; this should presumably be '36^4'. Please correct this typographical issue.
- [Section 4] The interpolation to the reference mass is described only verbally; the functional form (for example, a polynomial in M_PS^2), the fit ranges, and the chi^2 values are not given. Reporting these details would allow the reader to assess the interpolation error and the compatibility of the beta = 2.2 and beta = 2.15 points.
- [Figure 4] The continuum extrapolation plot shows points but no fit curve or error band; adding the fit line and the continuum value with its error would make the claimed O(a^2) behavior visible and would help the reader evaluate the 'slight tension' mentioned in the text.
- [Section 2.2 and Fig. 1] The c_sw tuning curve is said to be published in companion papers [8,9]; please state explicitly which data or figures are new in this proceedings and which are reproduced, so that the non-perturbative improvement claim does not rest on unpublished details.
Circularity Check
No circularity: f_PS is a direct twisted-mass matrix-element measurement; continuum extrapolation uses measured inputs, and self-citations are auxiliary tuning references.
full rationale
The paper's central quantity, the pseudoscalar decay constant, is obtained from Eq. (7), f^R_PS = 2 mu0 <0|P|pi>/M_PS^2, where mu0 is a bare twisted-mass parameter, <0|P|pi> is a measured matrix element, and M_PS is the measured pseudoscalar mass. No parameter in this formula is fitted to the f_PS values, so the result does not reduce by construction to an input. The continuum limit is taken by interpolating the measured w0 f_PS and a w0 at each beta to a common reference mass and extrapolating in (a/w0)^2; this is a standard scaling fit to measured data, not a renamed prediction. The paper's own caveat that one more point with finer lattice spacing is needed 'to confirm that the O(a) effects are eliminated' is an extrapolation and systematic-risk statement, not a circularity. The only self-references are the non-perturbative c_sw tuning curve described as 'also published as part of the work in [8, 9]' and the mixed-action correction formula (8) attributed to [17], whose authors overlap with this paper; neither is used to define the quoted w0 fPS = 0.1436(19) value, and the correction formula is an analytic relation, not an empirical fit. Hence the central claim has independent numerical content and no significant circularity.
Assumptions & free parameters
free parameters (1)
- c_sw(β) =
not tabulated; shown in Fig. 1 and published in [8,9]
assumptions (5)
- domain assumption At maximal twist, the pseudoscalar decay constant needs no renormalization constant (Eq. 7).
- domain assumption Finite-volume effects are negligible for the ensembles used.
- domain assumption Exact O(a) improvement makes residual discretization effects of order O(a^2).
- domain assumption The Schrödinger-functional c_sw curve from companion papers [8,9] is valid for these ensembles.
- domain assumption SU(2) gauge theory with two fundamental flavors is a relevant composite Higgs candidate.
Cite this review
Pith. "Pith review of Determination of the pseudoscalar decay constant from SU(2) with two fundamental flavors." pith.science (2026). https://pith.science/paper/R2YV2QTD
@misc{pith2026241206471,
author = {Pith},
title = {Pith review of: Determination of the pseudoscalar decay constant from SU(2) with two fundamental flavors},
year = {2026},
howpublished = {\url{https://pith.science/paper/R2YV2QTD}},
note = {Machine review of arXiv:2412.06471}
}
read the original abstract
The SU(2) gauge group with two fundamental flavors is a candidate for a composite Higgs extension of the Standard Model. Central to Higgs phenomenology is a non-perturbative determination of observables of the theory, such as the decay constant of the pseudo-Nambu-Goldstone Bosons. We present preliminary results for the continuum limit of the pseudoscalar decay constant using a mixed-action setup, with non-perturbatively improved stabilized Wilson Fermions on the sea, and maximally twisted valence quarks. Pivotal to this study is the recent porting of our simulation suite HiRep to GPU architecture.
Figures
Figures from the paper (1 more)
Forward citations
Cited by 1 Pith paper
-
The singlet scalar state in a chiral ensemble in $SU(2)$ with two fundamental flavours
The first lattice measurement of the sigma mass with exponential clover fermions in SU(2) with two fundamental flavours finds the effective mass agrees with the two-pion threshold, 2 m_PS.
Reference graph
Works this paper leans on
-
[6]
SU(2) gauge theory with two fundamental flavors: A minimal template for model building.Phys
Rudy Arthur, Vincent Drach, Martin Hansen, Ari Hietanen, Claudio Pica, and Francesco Sannino. SU(2) gauge theory with two fundamental flavors: A minimal template for model building.Phys. Rev. D, 94(9):094507, 2016
work page 2016
-
[1]
FundamentalCompositeDynamics: A Review.Phys
GiacomoCacciapaglia,ClaudioPica,andFrancescoSannino. FundamentalCompositeDynamics: A Review.Phys. Rept., 877:1–70, 2020
work page 2020
-
[2]
Higher representations on the lattice: Numerical simulations
Luigi Del Debbio, Agostino Patella, and Claudio Pica. Higher representations on the lattice: Numerical simulations. SU(2) with adjoint fermions.Phys. Rev. D, 81:094503, 2010
work page 2010
-
[3]
GPU- accelerated Higher Representations of Wilson Fermions with HiRep.PoS, EuroPLEx2023:035, 2024
Sofie Martins, Erik Kjellgren, Emiliano Molinaro, Claudio Pica, and Antonio Rago. GPU- accelerated Higher Representations of Wilson Fermions with HiRep.PoS, EuroPLEx2023:035, 2024. 7 Determination of the pseudoscalar decay constant from𝑆𝑈(2) with two fundamental flavors Sofie Martins
work page 2024
-
[4]
Scaling SU(2) to 1000 GPUs using HiRep
Sofie Martins, Erik Kjellgren, Emiliano Molinaro, Claudio Pica, and Antonio Rago. Scaling SU(2) to 1000 GPUs using HiRep. In41st International Symposium on Lattice Field Theory, 11 2024
work page 2024
-
[5]
Speeding up the hybrid Monte Carlo algorithm for dynamical fermions
Martin Hasenbusch. Speeding up the hybrid Monte Carlo algorithm for dynamical fermions. Phys. Lett. B, 519:177–182, 2001
work page 2001
-
[7]
Master-field simulations of O(𝑎)-improved lattice QCD: Algorithms, stability and exactness.Comput
Anthony Francis, Patrick Fritzsch, Martin Lüscher, and Antonio Rago. Master-field simulations of O(𝑎)-improved lattice QCD: Algorithms, stability and exactness.Comput. Phys. Commun., 255:107355, 2020
work page 2020
-
[8]
2-flavour𝑆𝑈(2) gauge theory with exponential clover Wilson fermions.PoS, LATTICE2023:094, 2024
Laurence Sebastian Bowes, Vincent Drach, Patrick Fritzsch, Antonio Rago, and Fernando Romero-Lopez. 2-flavour𝑆𝑈(2) gauge theory with exponential clover Wilson fermions.PoS, LATTICE2023:094, 2024
work page 2024
Show all 18 references
-
[9]
The mass of the𝜎 in a chiral ensemble in𝑆𝑈(2) with two fundamental flavours
Laurence Sebastian Bowes, Vincent Drach, Patrick Fritzsch, Sofie Martins, Antonio Rago, and Fernando Romero-López. The mass of the𝜎 in a chiral ensemble in𝑆𝑈(2) with two fundamental flavours. 2024. LATTICE 2024, to appear
2024
-
[10]
Nonperturbative O(a) improvement of lattice QCD.Nucl
Martin Luscher, Stefan Sint, Rainer Sommer, Peter Weisz, and Ulli Wolff. Nonperturbative O(a) improvement of lattice QCD.Nucl. Phys. B, 491:323–343, 1997
1997
-
[11]
Simulations with different lattice Dirac operators for valence and sea quarks.Phys
Oliver Bar, Gautam Rupak, and Noam Shoresh. Simulations with different lattice Dirac operators for valence and sea quarks.Phys. Rev. D, 67:114505, 2003
2003
-
[12]
A Local formulation of lattice QCD without unphysical fermion zero modes.Nucl
Roberto Frezzotti, Pietro Antonio Grassi, Stefan Sint, and Peter Weisz. A Local formulation of lattice QCD without unphysical fermion zero modes.Nucl. Phys. B Proc. Suppl., 83:941–946, 2000
2000
-
[13]
Lattice QCD with a chirally twisted mass term.JHEP, 08:058, 2001
Roberto Frezzotti, Pietro Antonio Grassi, Stefan Sint, and Peter Weisz. Lattice QCD with a chirally twisted mass term.JHEP, 08:058, 2001
2001
-
[14]
O(a) improved twisted mass lattice QCD.JHEP, 07:048, 2001
Roberto Frezzotti, Stefan Sint, and Peter Weisz. O(a) improved twisted mass lattice QCD.JHEP, 07:048, 2001
2001
-
[15]
Lattice QCD with a chiral twist
Stefan Sint. Lattice QCD with a chiral twist. InWorkshop on Perspectives in Lattice QCD, 2 2007
2007
-
[16]
Twisted mass lattice QCD.Phys
Andrea Shindler. Twisted mass lattice QCD.Phys. Rept., 461:37–110, 2008
2008
-
[17]
Large𝑁𝑐 scaling of meson masses and decay constants.Eur
Pilar Hernández, Carlos Pena, and Fernando Romero-López. Large𝑁𝑐 scaling of meson masses and decay constants.Eur. Phys. J. C, 79(10):865, 2019. 8 Determination of the pseudoscalar decay constant from𝑆𝑈(2) with two fundamental flavors Sofie Martins
2019
-
[18]
Singlet channel scattering in a composite Higgs model on the lattice.Eur
Vincent Drach, Patrick Fritzsch, Antonio Rago, and Fernando Romero-López. Singlet channel scattering in a composite Higgs model on the lattice.Eur. Phys. J. C, 82(1):47, 2022. 9
2022
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.