pith. sign in

arxiv: 2606.27072 · v1 · pith:D2XKXXP3new · submitted 2026-06-25 · ✦ hep-lat

Performance of Low Mode Averaging on Twisted-Mass Fermion Ensembles at the physical pion mass point

Constantia Alexandrou , Simone Bacchio , Antonio Evangelista , Roberto Frezzotti , Francesca Margari , Francesco Sanfilippo , Christian Schneider This is my paper

Pith reviewed 2026-06-26 01:35 UTC · model grok-4.3

classification ✦ hep-lat
keywords low-mode averagingtwisted-mass fermionsphysical pion masschiral condensateBanks-Casher relationhadronic vacuum polarisationlattice QCDcorrelation functions
0
0 comments X

The pith

Low-mode averaging reduces statistical noise and cost for meson and baryon correlation functions on twisted-mass ensembles at the physical pion mass.

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

The paper tests low-mode averaging on twisted-mass fermion lattices with quark masses close to physical values. It compares exact low-mode treatment against a multigrid approximation for meson two-point functions, including the vector correlator needed for the muon anomalous magnetic moment, and applies the exact method to baryon two- and three-point functions. Clear reductions in noise appear at large Euclidean times, accompanied by lower computational cost. The same ensembles are used to compute the eigenvalue density of the massless Wilson operator and to extract the renormalized chiral condensate through the Banks-Casher relation, giving 269.5(4.5) MeV at 2 GeV in the MSbar scheme. The pion-mass dependence of this condensate also yields the scale-independent low-energy constant h-bar_1 equal to 5.2(1.1).

Core claim

On twisted-mass ensembles at the physical pion mass, low-mode averaging applied to the Dirac operator yields quantifiable reductions in variance for light-quark meson and baryon two- and three-point functions at large Euclidean times, with the exact low-mode implementation outperforming the multigrid approximation for baryons while both versions perform well for mesons; the same data produce a renormalized chiral condensate of 269.5(4.5) MeV via the Banks-Casher relation.

What carries the argument

Low-mode averaging, which isolates the contribution of the lowest eigenmodes of the Dirac operator for exact evaluation while approximating the high-mode remainder to reduce variance in correlation functions.

Load-bearing premise

The low-mode spectrum extracted from the massless Wilson operator on these ensembles remains representative after the chiral extrapolation used to obtain the renormalized condensate at 2 GeV.

What would settle it

Repeating the noise-reduction measurements on an independent twisted-mass ensemble at the same physical pion mass and finding that the reported variance reductions are absent within errors would falsify the performance claims.

Figures

Figures reproduced from arXiv: 2606.27072 by Antonio Evangelista, Christian Schneider, Constantia Alexandrou, Francesca Margari, Francesco Sanfilippo, Roberto Frezzotti, Simone Bacchio.

Figure 1
Figure 1. Figure 1: FIG. 1. The normalised density of absolute eigenvalues [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Gain in the signal-to-noise ratio defined in Eq. (39), computed for the deflated two-point light-quark flavour non-singlet [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. On the C80 ensemble, the comparison of the two-point vector-vector correlator in isoQCD. The left panel shows the [PITH_FULL_IMAGE:figures/full_fig_p014_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. The gain in the signal-to-noise ratio in light-quark flavour non-singlet two-point correlation functions, as defined in [PITH_FULL_IMAGE:figures/full_fig_p015_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7 [PITH_FULL_IMAGE:figures/full_fig_p018_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Error reduction (top panels) and gain (bottom panels) in the statistical uncertainty of the correlation functions relative [PITH_FULL_IMAGE:figures/full_fig_p021_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. The measured (grey-blue dots) and predicted according to Eq. (67) (red solid line) runtime for the exact IR correlators [PITH_FULL_IMAGE:figures/full_fig_p023_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10 [PITH_FULL_IMAGE:figures/full_fig_p024_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11 [PITH_FULL_IMAGE:figures/full_fig_p025_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. The cost-gain [PITH_FULL_IMAGE:figures/full_fig_p027_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. The cost-gain [PITH_FULL_IMAGE:figures/full_fig_p027_13.png] view at source ↗
read the original abstract

We study the performance of low-mode averaging (LMA) on twisted-mass fermion ensembles at near-physical quark masses, assessing both its theoretical framework and practical cost-effectiveness in modern lattice QCD. In particular, we present a numerical study of light-quark meson and baryon observables. For mesons, we analyse two-point functions, including the vector-vector correlator relevant for the hadronic vacuum polarisation contribution to the muon anomalous magnetic moment, comparing two implementations of LMA: an exact approach based on explicit low modes and an approximate, high-statistics variant using multigrid techniques. For baryons, we restrict to the exact approach and study both two- and three-point functions, quantifying the resulting noise and cost reductions at large Euclidean times. In addition, we compute the eigenvalue density of the massless Wilson operator and determine the renormalised chiral condensate via the Banks-Casher relation, obtaining $\sqrt[3]{\Sigma_{\mathrm{R}}}=269.5(4.5)~\mathrm{MeV}$ for $N_f{=}2{+}1{+}1$ isospin-symmetric QCD at a scale $2~\mathrm{GeV}$ in the $\overline{\mathrm{MS}}$ scheme, with an uncertainty dominated by the chiral extrapolation. Additionally, from the pion-mass dependence of $\Sigma_{\mathrm{R}}$, we extract the scale-independent low-energy constant $\bar{h}_1=5.2(1.1)$.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript studies the performance of low-mode averaging (LMA) on twisted-mass fermion ensembles at near-physical quark masses. It reports quantitative noise and cost reductions for meson two-point functions (including the vector-vector correlator for the hadronic vacuum polarisation) using both exact low-mode and multigrid LMA implementations, and for baryon two- and three-point functions using the exact approach. In addition, the eigenvalue density of the massless Wilson operator is computed and the renormalised chiral condensate is extracted via the Banks-Casher relation, yielding √[3]{Σ_R} = 269.5(4.5) MeV at 2 GeV in the MSbar scheme (uncertainty dominated by chiral extrapolation) together with the low-energy constant ar{h}_1 = 5.2(1.1).

Significance. If the numerical demonstrations hold, the LMA results supply concrete, ensemble-specific evidence of efficiency gains for correlator computations at physical pion masses; this is directly relevant to reducing the cost of precision lattice QCD calculations such as those entering the muon anomalous magnetic moment. The Banks-Casher extraction adds an N_f=2+1+1 data point, though its quoted precision is limited by the extrapolation procedure.

major comments (1)
  1. [Banks-Casher analysis] Banks-Casher analysis (abstract and associated results): the statement that the uncertainty on √[3]{Σ_R} is dominated by the chiral extrapolation is load-bearing for the quoted central value and error, yet no details are supplied on the fit ansatz, the number or range of ensembles entering the extrapolation, or any volume-dependence checks; this prevents assessment of whether post-hoc choices affect the result.
minor comments (2)
  1. [Abstract] Clarify in the abstract and introduction whether the ensembles are exactly at the physical pion mass or only near-physical, and provide the precise m_π values used.
  2. [Results section on meson correlators] The comparison between exact and multigrid LMA for mesons would benefit from an explicit table of noise-reduction factors and CPU-cost ratios at the largest Euclidean times studied.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their positive assessment of the work and for the constructive comment on the Banks-Casher section. We address the point below and will revise the manuscript accordingly.

read point-by-point responses

  1. Referee: [Banks-Casher analysis] Banks-Casher analysis (abstract and associated results): the statement that the uncertainty on √[3]{Σ_R} is dominated by the chiral extrapolation is load-bearing for the quoted central value and error, yet no details are supplied on the fit ansatz, the number or range of ensembles entering the extrapolation, or any volume-dependence checks; this prevents assessment of whether post-hoc choices affect the result.

    Authors: We agree that the current manuscript does not supply sufficient detail on the chiral extrapolation procedure underlying the quoted uncertainty. In the revised version we will add an explicit description of the fit ansatz, list the ensembles (with their pion masses and spatial volumes) that enter the extrapolation, specify the fit range, and report the outcome of any volume-dependence checks that were performed. This will make the error budget fully transparent. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper's central results are direct numerical measurements of noise reduction and computational cost savings achieved by exact and multigrid LMA implementations on twisted-mass ensembles for meson and baryon correlators at large Euclidean times. These are empirical outcomes from explicit computations on the given ensembles, not predictions derived from fitted parameters or self-referential definitions. The Banks-Casher extraction of the condensate applies an external standard relation to the computed eigenvalue spectrum of the massless Wilson operator, followed by a chiral extrapolation whose dominant uncertainty is explicitly flagged; this does not reduce the main LMA claims to the input data by construction. No load-bearing self-citations, uniqueness theorems, or ansatze smuggled via prior work appear in the derivation chain. The work is self-contained as a performance benchmark study.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

Only the abstract is available; the ledger is therefore populated from statements that must be true for the reported numbers to hold. The chiral extrapolation and the applicability of Banks-Casher on finite-volume twisted-mass ensembles are the main unstated inputs.

free parameters (1)
  • chiral extrapolation parameters
    The uncertainty on the condensate is stated to be dominated by the chiral extrapolation, implying at least one fitted parameter or functional form chosen to reach the physical point.
axioms (2)
  • domain assumption Banks-Casher relation holds for the renormalized eigenvalue density of the massless Wilson operator on these ensembles
    Invoked to convert eigenvalue density into chiral condensate value.
  • domain assumption Twisted-mass fermion action at the simulated parameters produces ensembles representative of Nf=2+1+1 QCD at physical pion mass after extrapolation
    Required for the performance claims and the condensate result to apply to real-world QCD.

pith-pipeline@v0.9.1-grok · 5811 in / 1512 out tokens · 18611 ms · 2026-06-26T01:35:49.318761+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

49 extracted references · 25 linked inside Pith

  1. [1]

    In this case, the validity of Eq

    Quark loops Consider a generic quark loop of the form L(t, ⃗ p,Γ) = Tr T t q ΓS ,(A1) whereT t q projects onto a momentum boostqand time-slicetand Γ is an arbitrary Dirac matrix. In this case, the validity of Eq. (6) is immediate, but already illustrates the mechanism at work. Indeed, since for quark loops the observable depends linearly on the propagator...

  2. [2]

    Meson two-point correlation functions with stochastic sources Consider a generic meson two-point correlation function, computed using a backwards-running propagator via the γ5-hermiticity,S=γ 5S†γ54, defined as M(t1, t2,p 1,p 2,Γ 1,Γ 2) = Tr T t2 p2 γ5Γ2 S T t1 p1 Γ1γ5 S† .(A3) Then ifH η satisfies HηHη =H η , H † η =H η ,and [H η, T t1 p1 Γ1γ5] = 0,(A4) ...

  3. [3]

    (6) holds when point sources are employed

    Baryon two-point correlation functions with point sources The final example we consider is a baryon two-point correlation function, for which we show that Eq. (6) holds when point sources are employed. For simplicity, we consider a baryon correlation function of the form B(tx, ty,p,P,W) = X x,y ei(x−y)·p Pγγ ′ W abc αβ W a′b′c′ α′β′ S(x, y)αα′ aa′ S(x, y)...

  4. [4]

    Walker-Loud et al.,Lattice QCD Determination ofg A,PoSCD2018(2020) 020, [1912.08321]

    A. Walker-Loud et al.,Lattice QCD Determination ofg A,PoSCD2018(2020) 020, [1912.08321]. [2]RQCDcollaboration, G. S. Bali, S. Collins, S. Heybrock, M. L¨ offler, R. R¨ odl, W. S¨ oldner et al.,Octet baryon isovector charges from Nf=2+1 lattice QCD,Phys. Rev. D108(2023) 034512, [2305.04717]

    arXiv 2020

  5. [5]

    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]

    arXiv 2024

  6. [6]

    Alexandrou, S

    C. Alexandrou, S. Bacchio, J. Finkenrath, C. Iona, G. Koutsou, Y. Li et al.,Nucleon charges andσ-terms in lattice QCD, Phys. Rev. D111(2025) 054505, [2412.01535]. [5]Flavour Lattice A veraging Group (FLAG)collaboration, Y. Aoki et al.,FLAG review 2024,Phys. Rev. D113 (2026) 014508, [2411.04268]. [6]Precision Neutron Decay Matrix Elements (PNDME)collaborat...

    arXiv 2025

  7. [7]

    Djukanovic, G

    D. Djukanovic, G. von Hippel, J. Koponen, H. B. Meyer, K. Ottnad, T. Schulz et al.,Isovector axial form factor of the nucleon from lattice QCD,Phys. Rev. D106(2022) 074503, [2207.03440]. [8]Extended Twisted Masscollaboration, C. Alexandrou, S. Bacchio, M. Constantinou, J. Finkenrath, R. Frezzotti, B. Kostrzewa et al.,Nucleon axial and pseudoscalar form fa...

    arXiv 2022

  8. [8]

    Djukanovic, G

    D. Djukanovic, G. von Hippel, H. B. Meyer, K. Ottnad, M. Salg and H. Wittig,Electromagnetic form factors of the nucleon from Nf=2+1 lattice QCD,Phys. Rev. D109(2024) 094510, [2309.06590]. [10]PACScollaboration, R. Tsuji, Y. Aoki, K.-I. Ishikawa, Y. Kuramashi, S. Sasaki, K. Sato et al.,Nucleon form factors in Nf=2+1 lattice QCD at the physical point: Finit...

    arXiv 2024

  9. [9]

    Alexandrou, S

    C. Alexandrou, S. Bacchio, G. Koutsou, B. Prasad and G. Spanoudes,Proton and neutron electromagnetic form factors from lattice QCD in the continuum limit,2507.20910. [12]Extended Twisted Mass Collaboration (ETMC)collaboration, C. Alexandrou et al.,Probing the Energy-Smeared R Ratio Using Lattice QCD,Phys. Rev. Lett.130(2023) 241901, [2212.08467]. [13]Exte...

    arXiv 2023

  10. [10]

    De Santis et al.,Inclusive semileptonic decays of the Ds meson: A first-principles lattice QCD calculation,Phys

    A. De Santis et al.,Inclusive semileptonic decays of the Ds meson: A first-principles lattice QCD calculation,Phys. Rev. D112(2025) 054503, [2504.06063]

    arXiv 2025

  11. [11]

    Borsanyi et al.,Leading hadronic contribution to the muon magnetic moment from lattice QCD,Nature593(2021) 32 51–55, [2002.12347]

    S. Borsanyi et al.,Leading hadronic contribution to the muon magnetic moment from lattice QCD,Nature593(2021) 32 51–55, [2002.12347]

    arXiv 2021

  12. [12]

    Djukanovic, G

    D. Djukanovic, G. von Hippel, S. Kuberski, H. B. Meyer, N. Miller, K. Ottnad et al.,The hadronic vacuum polarization contribution to the muon g−2 at long distances,JHEP04(2025) 098, [2411.07969]

    arXiv 2025

  13. [13]

    Boccaletti et al.,Hybrid calculation of hadronic vacuum polarization in muon g−2 to 0.48%,Nature653(2026) 373–377, [2407.10913]

    A. Boccaletti et al.,Hybrid calculation of hadronic vacuum polarization in muon g−2 to 0.48%,Nature653(2026) 373–377, [2407.10913]. [20]Fermilab Lattice, HPQCD,, MILCcollaboration, A. Bazavov et al.,Hadronic Vacuum Polarization for the Muon g-2 from Lattice QCD: Long-Distance and Full Light-Quark Connected Contribution,Phys. Rev. Lett.135(2025) 011901, [2...

    Pith/arXiv arXiv 2026

  14. [14]

    H. Neff, N. Eicker, T. Lippert, J. W. Negele and K. Schilling,On the low fermionic eigenmode dominance in QCD on the lattice,Phys. Rev. D64(2001) 114509, [hep-lat/0106016]

    Pith/arXiv arXiv 2001

  15. [15]

    Giusti, P

    L. Giusti, P. Hernandez, M. Laine, P. Weisz and H. Wittig,Low-energy couplings of QCD from current correlators near the chiral limit,JHEP04(2004) 013, [hep-lat/0402002]

    Pith/arXiv arXiv 2004

  16. [16]

    T. A. DeGrand and S. Schaefer,Improving meson two point functions in lattice QCD,Comput. Phys. Commun.159 (2004) 185–191, [hep-lat/0401011]

    Pith/arXiv arXiv 2004

  17. [17]

    T. A. DeGrand and S. Schaefer,Chiral properties of two-flavor QCD in small volume and at large lattice spacing,Phys. Rev. D72(2005) 054503, [hep-lat/0506021]

    Pith/arXiv arXiv 2005

  18. [18]

    Giusti and S

    L. Giusti and S. Necco,Low-mode averaging for baryon correlation functions,PoSLA T2005(2006) 132, [hep-lat/0510011]

    Pith/arXiv arXiv 2006

  19. [19]

    G. S. Bali, S. Collins and A. Schafer,Effective noise reduction techniques for disconnected loops in Lattice QCD,Comput. Phys. Commun.181(2010) 1570–1583, [0910.3970]

    Pith/arXiv arXiv 2010

  20. [20]

    T. Blum, T. Izubuchi and E. Shintani,New class of variance-reduction techniques using lattice symmetries,Phys. Rev. D 88(2013) 094503, [1208.4349]

    Pith/arXiv arXiv 2013

  21. [21]

    Shintani, R

    E. Shintani, R. Arthur, T. Blum, T. Izubuchi, C. Jung and C. Lehner,Covariant approximation averaging,Phys. Rev. D 91(2015) 114511, [1402.0244]

    Pith/arXiv arXiv 2015

  22. [22]

    Luscher,Local coherence and deflation of the low quark modes in lattice QCD,JHEP07(2007) 081, [0706.2298]

    M. Luscher,Local coherence and deflation of the low quark modes in lattice QCD,JHEP07(2007) 081, [0706.2298]. [32]xQCDcollaboration, A. Li et al.,Overlap Valence on 2+1 Flavor Domain Wall Fermion Configurations with Deflation and Low-mode Substitution,Phys. Rev. D82(2010) 114501, [1005.5424]

    Pith/arXiv arXiv 2007

  23. [23]

    Luscher,Deflation acceleration of lattice QCD simulations,JHEP12(2007) 011, [0710.5417]

    M. Luscher,Deflation acceleration of lattice QCD simulations,JHEP12(2007) 011, [0710.5417]

    Pith/arXiv arXiv 2007

  24. [24]

    Brannick, R

    J. Brannick, R. C. Brower, M. A. Clark, J. C. Osborn and C. Rebbi,Adaptive Multigrid Algorithm for Lattice QCD, Phys. Rev. Lett.100(2008) 041601, [0707.4018]

    Pith/arXiv arXiv 2008

  25. [25]

    Babich, J

    R. Babich, J. Brannick, R. C. Brower, M. A. Clark, T. A. Manteuffel, S. F. McCormick et al.,Adaptive multigrid algorithm for the lattice Wilson-Dirac operator,Phys. Rev. Lett.105(2010) 201602, [1005.3043]

    Pith/arXiv arXiv 2010

  26. [26]

    Frommer, K

    A. Frommer, K. Kahl, S. Krieg, B. Leder and M. Rottmann,Adaptive Aggregation-Based Domain Decomposition Multigrid for the Lattice Wilson–Dirac Operator,SIAM J. Sci. Comput.36(2014) A1581–A1608, [1303.1377]

    Pith/arXiv arXiv 2014

  27. [27]

    Frommer, K

    A. Frommer, K. Kahl, S. Krieg, B. Leder and M. Rottmann,An adaptive aggregation based domain decomposition multilevel method for the lattice wilson dirac operator: multilevel results,1307.6101

    Pith/arXiv arXiv

  28. [28]

    Alexandrou, S

    C. Alexandrou, S. Bacchio, J. Finkenrath, A. Frommer, K. Kahl and M. Rottmann,Adaptive Aggregation-based Domain Decomposition Multigrid for Twisted Mass Fermions,Phys. Rev. D94(2016) 114509, [1610.02370]

    Pith/arXiv arXiv 2016

  29. [29]

    Alexandrou, S

    C. Alexandrou, S. Bacchio and J. Finkenrath,Multigrid approach in shifted linear systems for the non-degenerated twisted mass operator,Comput. Phys. Commun.236(2019) 51–64, [1805.09584]

    Pith/arXiv arXiv 2019

  30. [30]

    Brannick, A

    J. Brannick, A. Frommer, K. Kahl, B. Leder, M. Rottmann and A. Strebel,Multigrid Preconditioning for the Overlap Operator in Lattice QCD,Numer. Math.132(2016) 463–490, [1410.7170]

    Pith/arXiv arXiv 2016

  31. [31]

    R. C. Brower, M. A. Clark, D. Howarth and E. S. Weinberg,Multigrid for chiral lattice fermions: Domain wall,Phys. Rev. D102(2020) 094517, [2004.07732]

    arXiv 2020

  32. [32]

    R. C. Brower, M. A. Clark, A. Strelchenko and E. Weinberg,Multigrid algorithm for staggered lattice fermions,Phys. Rev. D97(2018) 114513, [1801.07823]

    Pith/arXiv arXiv 2018

  33. [33]

    Bacchio, C

    S. Bacchio, C. Alexandrou and J. Finkerath,Multigrid accelerated simulations for Twisted Mass fermions,EPJ Web Conf.175(2018) 02002, [1710.06198]

    Pith/arXiv arXiv 2018

  34. [34]

    Finkenrath, C

    J. Finkenrath, C. Alexandrou, S. Bacchio, P. Charalambous, P. Dimopoulos, R. Frezzotti et al.,Simulation of an ensemble ofN f = 2 + 1 + 1twisted mass clover-improved fermions at physical quark masses,EPJ Web Conf.175(2018) 02003, [1712.09579]

    Pith/arXiv arXiv 2018

  35. [35]

    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. D 98(2018) 054518, [1807.00495]

    Pith/arXiv arXiv 2018

  36. [36]

    Frommer, M

    A. Frommer, M. N. Khalil and G. Ramirez-Hidalgo,A Multilevel Approach to Variance Reduction in the Stochastic Estimation of the Trace of a Matrix,SIAM J. Sci. Comput.44(2022) A2536–A2556, [2108.11281]

    arXiv 2022

  37. [37]

    Gruber, T

    R. Gruber, T. Harris and M. Krstic Marinkovic,Multigrid low-mode averaging,Phys. Rev. D111(2025) 074508, [2412.06347]

    arXiv 2025

  38. [38]

    Frommer, J

    A. Frommer, J. Jimenez-Merchan, F. Knechtli, T. Korzec and G. Ramirez-Hidalgo,Using orthogonal projectors in multigrid multilevel Monte Carlo for trace estimation in lattice QCD,2509.11424. 33

    arXiv

  39. [39]

    Banks and A

    T. Banks and A. Casher,Chiral Symmetry Breaking in Confining Theories,Nucl. Phys. B169(1980) 103–125

    1980

  40. [40]

    Leutwyler and A

    H. Leutwyler and A. V. Smilga,Spectrum of Dirac operator and role of winding number in QCD,Phys. Rev. D46(1992) 5607–5632

    1992

  41. [41]

    Giusti and M

    L. Giusti and M. Luscher,Chiral symmetry breaking and the Banks-Casher relation in lattice QCD with Wilson quarks, JHEP03(2009) 013, [0812.3638]. [52]Alphacollaboration, R. Frezzotti, P. A. Grassi, S. Sint and P. Weisz,Lattice QCD with a chirally twisted mass term, JHEP08(2001) 058, [hep-lat/0101001]. [53]Extended Twisted Masscollaboration,Quark masses an...

    Pith/arXiv arXiv 2009

  42. [42]

    Detmold, D

    W. Detmold, D. J. Murphy, A. V. Pochinsky, M. J. Savage, P. E. Shanahan and M. L. Wagman,Sparsening algorithm for multihadron lattice QCD correlation functions,Phys. Rev. D104(2021) 034502, [1908.07050]

    arXiv 2021

  43. [43]

    Li, S.-C

    Y. Li, S.-C. Xia, X. Feng, L.-C. Jin and C. Liu,Field sparsening for the construction of the correlation functions in lattice QCD,Phys. Rev. D103(2021) 014514, [2009.01029]

    arXiv 2021

  44. [44]

    Christian and W

    S. Christian and W. Detmold,Aspects of propagator sparsening in lattice QCD,Phys. Rev. D111(2025) 114514, [2501.05404]. [57]Extended Twisted Masscollaboration, C. Alexandrou et al.,Quark masses using twisted-mass fermion gauge ensembles,Phys. Rev. D104(2021) 074515, [2104.13408]

    arXiv 2025

  45. [45]

    Bonanno, F

    C. Bonanno, F. D’Angelo and M. D’Elia,The chiral condensate of N f = 2 + 1 QCD from the spectrum of the staggered Dirac operator,JHEP11(2023) 013, [2308.01303]. [59]Flavour Lattice A veraging Group (FLAG)collaboration, Y. Aoki et al.,FLAG Review 2021,Eur. Phys. J. C82 (2022) 869, [2111.09849]

    arXiv 2023

  46. [46]

    Gasser and H

    J. Gasser and H. Leutwyler,Chiral Perturbation Theory to One Loop,Annals Phys.158(1984) 142

    1984

  47. [47]

    Amoros, J

    G. Amoros, J. Bijnens and P. Talavera,Two point functions at two loops in three flavor chiral perturbation theory,Nucl. Phys. B568(2000) 319–363, [hep-ph/9907264]

    Pith/arXiv arXiv 2000

  48. [48]

    Bijnens,Private communication, 2026

    J. Bijnens,Private communication, 2026. [63]QUDAcollaboration, M. A. Clark, R. Babich, K. Barros, R. C. Brower and C. Rebbi,Solving Lattice QCD systems of equations using mixed precision solvers on GPUs,Comput. Phys. Commun.181(2010) 1517–1528, [0911.3191]. [64]QUDAcollaboration, R. Babich, M. A. Clark, B. Joo, G. Shi, R. C. Brower and S. Gottlieb,Scaling...

    Pith/arXiv arXiv 2026

  49. [49]

    J¨ ulich Supercomputing Centre,JUWELS Cluster and Booster: Exascale Pathfinder with Modular Supercomputing Architecture at Juelich Supercomputing Centre,Journal of large-scale research facilities7(2021)

    2021