---
title: Twisted Mass Clover-Improved Fermions
url: https://www.emergentmind.com/topics/twisted-mass-clover-improved-fermions
type: topic
---

# Twisted Mass Clover-Improved Fermions

Searching arXiv for recent and foundational papers on twisted-mass clover-improved fermions.
Twisted mass clover-improved fermions are a Wilson-type lattice fermion formulation in which the twisted-mass action is supplemented by a Sheikholeslami–Wohlert clover term. In the literature, this formulation appears in both two-flavor and \(N_f=2+1+1\) realizations, typically tuned to maximal twist so that physical observables are automatically \(\mathcal O(a)\)-improved, while the clover term is used to suppress residual \(\mathcal O(a^2)\) cutoff effects, especially twisted-mass-induced isospin breaking in the pion sector and the associated instabilities near the physical point [1411.6842, 1712.09579, 1807.00495, 2104.06747].

## 1. Action, field content, and basis structure

The formulation combines the Wilson kernel, a twisted mass term, and a clover term. For \(N_f=2+1+1\) simulations, the light and heavy sectors are commonly written as separate twisted-basis doublets \(\chi_\ell=(u,d)^T\) and \(\chi_h=(s,c)^T\). One explicit form used in recent spectroscopy work is [2309.04401]
\[
\begin{split}
S_{tm}^\ell & = \sum_x \bar{\chi}_\ell(x)\left[D_W[U] + \frac{i}{4}c_{SW}\sigma_{\mu\nu}\mathcal{F}^{\mu\nu} + m_\ell + i\mu_\ell\gamma_5\tau^3\right]\chi_\ell(x),\\
S_{tm}^h & = \sum_x \bar{\chi}_h(x)\left[D_W[U] + \frac{i}{4}c_{SW}\sigma_{\mu\nu}\mathcal{F}^{\mu\nu} + m_h - \mu_\delta\tau^1 + i\mu_\sigma\gamma_5\tau^3\right]\chi_h(x),
\end{split}
\]
with the massless Wilson–Dirac operator
\[
D_W[U] = \frac{1}{2}\gamma_\mu\left(\nabla_\mu + \nabla_\mu^\star\right) - \frac{ar}{2}\nabla_\mu\nabla_\mu^\star.
\]

For \(N_f=2\) simulations at the physical point, the light-doublet action is often written as [1704.02647, 1507.05068]
\[
S_F[\chi,\bar\chi,U] = a^4 \sum_x \bar\chi(x)\left( D_W[U] + m_0 + i\mu_l\gamma_5\tau^3 -\frac{1}{4} c_{\rm SW}\sigma^{\mu\nu}\mathcal F^{\mu\nu}[U] \right)\chi(x).
\]
Representative studies use \(c_{\rm SW}=1.57551\) in this two-flavor setting [1704.02647, 1507.05068].

In ETMC-style \(N_f=2+1+1\) production runs, the heavy nondegenerate doublet can also be written in pseudofermion form through
\[
D_{1+1}(\mu_\sigma,\mu_\delta) = D_W(\kappa,c_{SW}) \otimes 1 + i \mu_\sigma (\gamma_5 \otimes \tau_3) - \mu_\delta (1 \otimes \tau_1),
\]
or equivalently as a flavor-block matrix, making explicit the average heavy twisted mass \(\mu_\sigma\) and the strange–charm splitting parameter \(\mu_\delta\) [1712.09579].

A recurrent feature of the literature is that some methodological papers write the action explicitly, whereas several application papers identify the discretization family but defer the full light and heavy actions, clover coefficient, and tuning details to external ensemble references. This is stated explicitly in disconnected-nucleon, strange-form-factor, and meson-form-factor proceedings, which use the formulation operationally rather than rederive it [1511.00433, 2605.02049, 2112.03953].

## 2. Maximal twist, automatic improvement, and the role of the clover term

The defining improvement mechanism is maximal twist. In the twisted-mass framework, the untwisted bare mass is tuned to criticality so that the PCAC mass vanishes,
\[
m_{\rm PCAC}\to 0,
\]
which yields automatic \(\mathcal O(a)\) improvement for physical observables [1411.6842, 1904.10013, 1712.09579, 2104.06747]. In practical ETMC tuning, this is implemented by choosing the critical hopping parameter or equivalently \(m_0=m_{\rm cr}\), with residual mistuning monitored through ratios such as \(Z_A m_{\rm PCAC}/\mu_\ell\) [1712.09579, 1807.00495, 2001.09116].

The clover term is not introduced to obtain \(\mathcal O(a)\) improvement; that property already follows from maximal twist. Its purpose is to reduce the remaining \(\mathcal O(a^2)\) artifacts, most notably twisted-mass flavor breaking and the charged–neutral pion mass splitting. This point is stated directly in several papers: twisted mass at maximal twist retains automatic improvement, but the clover term suppresses the large \(\mathcal O(a^2)\) flavor-breaking artifacts that otherwise obstruct simulations close to the physical point [1411.6842, 1507.05068, 1712.09579, 1807.00495].

The quantitative evidence is formulation-defining. A dedicated \(N_f=2\) study reported a charged–neutral pion mass splitting of about \(20(20)\,\mathrm{MeV}\), including the disconnected contribution, and described this as roughly a factor-of-five reduction relative to earlier non-clover twisted-mass simulations on coarse lattices [1411.6842]. A subsequent physical-point \(N_f=2\) study found the full pion splitting compatible with zero within errors, with the charged–neutral difference at the physical point no larger than about \(13\,\mathrm{MeV}\), and the connected neutral-pion splitting about a factor of three smaller than in earlier twisted-mass simulations without clover [1507.05068]. In the first physical-point \(N_f=2+1+1\) clover-improved ensemble, the neutral–charged pion splitting was reported to be reduced by about a factor of six relative to earlier non-clover \(N_f=2+1+1\) simulations at similar lattice spacing [1807.00495].

This reduction has direct dynamical significance. Several papers connect the neutral-pion problem to the Sharpe–Singleton scenario and emphasize that the clover term is crucial for stable physical-point simulations at lattice spacings of order \(0.08\)–\(0.10\,\mathrm{fm}\) [1311.4522, 1712.09579, 2001.09116, 2104.06747]. A plausible implication is that twisted mass and clover improvement should be viewed as complementary rather than redundant: maximal twist controls odd-\(a\) effects, while the clover term reduces the size of the even-\(a\) artifacts that remain.

## 3. Realizations in \(N_f=2\) and \(N_f=2+1+1\) simulations

The formulation has been deployed in two principal dynamical-flavor settings. The two-flavor program established the viability of physical-point simulations with lattice spacing \(a\simeq 0.09\,\mathrm{fm}\), typically on \(48^3\times 96\) lattices with pion masses near \(130\,\mathrm{MeV}\) [1311.4522, 1507.05068, 1704.02647, 1807.11203]. These studies use degenerate light sea quarks and often add strange and charm only in the valence sector through Osterwalder–Seiler fermions, creating a mixed-action setup for heavy-flavor observables while keeping the light sea at maximal twist [1507.05068, 1704.02647].

The \(N_f=2+1+1\) program extends the formulation to fully dynamical light, strange, and charm quarks. A key milestone was the production of the first physical-point \(N_f=2+1+1\) clover-improved twisted-mass ensemble cB211.072.64 with
\[
128\times 64^3,\qquad \beta=1.778,\qquad a\mu_\ell=0.00072,\qquad a\mu_\sigma=0.1246864,\qquad a\mu_\delta=0.1315052,\qquad \kappa=0.1394265,\qquad c_{SW}=1.69,
\]
and lattice spacing \(a=0.08029(41)\,\mathrm{fm}\) [1807.00495]. Later work expanded this into a three-spacing and then four-spacing physical-point program with lattice spacings around \(0.080\), \(0.069\), \(0.057\), and \(0.049\,\mathrm{fm}\), spatial extents around \(5.1\)–\(5.5\,\mathrm{fm}\), and \(m_\pi L>3.6\), enabling continuum extrapolations directly at the physical pion mass [2309.04401, 2605.02049].

A characteristic feature of these \(N_f=2+1+1\) studies is the use of Osterwalder–Seiler valence strange and charm quarks on top of the dynamical twisted-mass sea when baryon or meson spectroscopy requires clean heavy-flavor interpolation. This is done to avoid unwanted \(s\)–\(c\) mixing through cutoff effects while preserving automatic \(\mathcal O(a)\) improvement [1704.02647, 2001.09116, 2104.06747, 2309.04401].

The formulation is also used in near-physical and heavier-than-physical ensembles for structure observables. Examples include \(m_\pi\simeq 260\)–\(265\,\mathrm{MeV}\) ensembles at \(a\simeq 0.093\)–\(0.095\,\mathrm{fm}\) for gluon and meson form factors, and physical-point large-volume ensembles at \(64^3\times128\) and larger for nucleon structure [2310.01389, 2112.03953, 1904.10013].

## 4. Algorithmic consequences and stochastic estimators

Twisted mass clover-improved fermions are not only a discretization choice; they also shape the estimator technology used in demanding disconnected and boosted-hadron calculations. The most distinctive formulation-specific ingredient is the one-end trick. In a physical-point disconnected-nucleon study, the twisted-mass flavor identities were written as [1511.00433]
\[
\Gamma\left[G_u - G_d\right] = -2i\mu a\,\Gamma G_d\gamma_5G_u,
\]
leading to the loop estimators
\[
\mathrm{Tr}\left[\Gamma(G_u-G_d)\right] = -\frac{2i\mu a}{N}\sum_{j=1}^N \left\langle s_j \Gamma\gamma_5 s_j\right\rangle
\]
and
\[
\mathrm{Tr}\left[\Gamma(G_u+G_d)\right] = \frac{2}{N}\sum_{j=1}^N \left\langle s_j \gamma_5\Gamma\gamma_5 D_W s_j\right\rangle.
\]
These identities permit volume-summed stochastic contractions and substantial variance reduction. The same literature emphasizes that the resulting all-site support gives the fermion loop at all insertion times simultaneously [1511.00433, 2605.02049].

Physical-point disconnected studies combine this twisted-mass structure with solver acceleration. For light disconnected loops, exact deflation of the Hermitian operator
\[
Q=\gamma_5 M
\]
is used, typically with Lanczos eigensolvers; for strange and charm loops, truncated-solver methods are employed where low-precision and high-precision solves remain sufficiently correlated [1511.00433, 1807.11203]. One physical-point \(N_f=2\) study reports 600 eigenpairs and 2250 stochastic sources per configuration for light loops, with the inversion cost becoming negligible after deflation, while strange loops use 63 high-precision and 1024 low-precision sources [1511.00433].

In \(N_f=2+1+1\) strange-form-factor calculations, the generalized one-end trick is combined with spin-color dilution and hierarchical probing using a four-dimensional coloring of distance eight [2605.02049]. In the gluon-PDF program, the same clover-improved twisted-mass ensemble supports very high statistics and a simultaneous treatment of gluon–quark singlet mixing; the quark-singlet part uses low-mode deflation, hierarchical probing, the one-end trick, and full spin/color dilution [2310.01389].

The ensemble-generation side is equally modern. Physical-mass \(N_f=2+1+1\) production uses HMC for the light sector, RHMC for the heavy sector, Hasenbusch mass preconditioning, DD-\(\alpha\)AMG multigrid, and nested minimal-norm integrators, with reported acceptance rates around \(77\%\) and healthy topological tunneling [1712.09579]. This computational infrastructure is part of the practical identity of the formulation in current ETMC-style simulations.

## 5. Spectroscopy, decay constants, and structure observables

The formulation has been used across a broad physics program. In low-lying baryon spectroscopy, a two-flavor physical-point calculation on a \(48^3\times96\) lattice found hyperon and charmed-baryon masses in agreement with experiment and reported that the clover term suppresses isospin symmetry breaking relative to earlier \(N_f=2+1+1\) twisted-mass results [1704.02647]. A later three-spacing \(N_f=2+1+1\) calculation carried the baryon spectrum to the continuum limit directly at the physical point, finding that most isospin splittings are already consistent with zero at finite lattice spacing and that the remaining ones vanish in the continuum limit [2309.04401].

In the pseudoscalar sector, physical-point \(N_f=2\) studies extracted quark-mass ratios and decay constants such as
\[
\mu_s/\mu_l=27.63(13),\qquad \mu_c/\mu_l=339.6(2.2),\qquad \mu_c/\mu_s=12.29(10),
\]
together with
\[
f_K=153.9(7.5)\,\mathrm{MeV},\qquad f_D=219(11)\,\mathrm{MeV},\qquad f_{D_s}=255(12)\,\mathrm{MeV},
\]
while explicitly noting the absence of continuum and finite-volume corrections in that early stage [1411.6842]. A later \(N_f=2\) physical-point paper quoted
\[
f_K/f_\pi = 1.1976(21)\left({}^{+06}_{-07}\right),
\]
and renormalized light, strange, and charm masses in the \(\overline{\rm MS}\) scheme at \(2\,\mathrm{GeV}\) [1507.05068]. In the mature \(N_f=2+1+1\) program, the ratio
\[
(f_K/f_\pi)^{\rm isoQCD}=1.1995(44)
\]
and the gradient-flow scales
\[
w_0=0.17383(63)\,\mathrm{fm},\qquad \sqrt{t_0}=0.14436(61)\,\mathrm{fm},\qquad t_0/w_0=0.11969(62)\,\mathrm{fm}
\]
were obtained from near-physical multi-spacing ensembles [2104.06747, 2111.14710].

Nucleon structure applications are especially prominent. Physical-point \(N_f=2\) studies used twisted-mass clover-improved ensembles for disconnected contributions to nucleon sigma terms, axial charges, tensor charges, and form factors, combining exact deflation, one-end-trick estimators, and large statistics [1511.00433, 1807.11203]. In strange electromagnetic form factors, four \(N_f=2+1+1\) physical-point clover-improved ensembles at \(a=0.080\), \(0.068\), \(0.057\), and \(0.049\,\mathrm{fm}\) were used to take the continuum limit directly at the physical pion mass [2605.02049].

The formulation also supports boosted-hadron and nonlocal-operator programs. A gluon-PDF study employed one \(N_f=2+1+1\) maximally twisted-mass clover-improved ensemble with Iwasaki gluons, \(a=0.0938(2)(3)\,\mathrm{fm}\), \(32^3\times64\), and \(m_\pi=260\,\mathrm{MeV}\), using the pseudo-distribution method and removing gluon–quark singlet mixing on the same ensemble [2310.01389]. Meson scalar, vector, and tensor form factors were likewise computed on a clover-improved \(N_f=2+1+1\) ensemble at \(m_\pi\approx265\,\mathrm{MeV}\) [2112.03953].

## 6. Scaling patterns, limitations, and conceptual boundaries

Because maximally twisted clover-improved Wilson fermions are expected to have leading cutoff effects of order \(a^2\), many continuum analyses are built around linear \(a^2\) ansätze. This is explicit in continuum extrapolations of strange form factors,
\[
c_k(a^2)=c_{k,0}+a^2 c_{k,2},
\]
in dipole and \(z\)-expansion fits [2605.02049], in baryon spectroscopy with coefficients and masses extrapolated linearly in \(a^2\) [2309.04401], and in scale-setting and decay-constant fits where \(a^2\) and \(a^2\xi\) terms suffice to describe the data [2104.06747, 2111.14710]. These studies report mild linear dependence in \(a^2\), consistent with the standard maximal-twist expectation.

At the same time, the literature is explicit about limitations. Several proceedings and application papers identify the discretization family but do not print the full action, twisted/physical-basis rotations, clover coefficient, or maximal-twist tuning condition, relying instead on earlier ETMC ensemble papers [1511.00433, 2605.02049, 1807.11203, 2112.03953]. Early physical-point studies were limited to one lattice spacing and often one volume, so finite-volume and discretization effects could only be inferred indirectly [1411.6842, 1507.05068, 1704.02647]. Mixed-action use of Osterwalder–Seiler valence strange and charm quarks is a standard expedient rather than a fully unitary formulation [1507.05068, 1704.02647, 2104.06747, 2309.04401].

A further conceptual boundary is that twisted mass and clover improvement should not be conflated. A recent Hamiltonian Schwinger-model study analyzes twisted-mass Wilson fermions without a clover term and confirms automatic \(\mathcal O(a)\) improvement at maximal twist while also exhibiting finite-\(a\) isospin breaking [2509.02329]. This clarifies that the twisted-mass mechanism and the clover suppression of residual artifacts are distinct ingredients. In the lattice-QCD literature surveyed here, “twisted mass clover-improved fermions” therefore denotes a compound strategy: maximal twist provides automatic improvement, and the clover term makes the physical-point formulation numerically stable and phenomenologically useful by reducing the dominant residual \(\mathcal O(a^2)\) distortions [1411.6842, 1712.09579, 1807.00495, 2104.06747].

Source: https://www.emergentmind.com/topics/twisted-mass-clover-improved-fermions