---
title: 'HISQ: Highly Improved Staggered Quark Action'
url: https://www.emergentmind.com/topics/hisq
type: topic
---

# HISQ: Highly Improved Staggered Quark Action

Searching arXiv for recent HISQ papers and foundational references.
arxiv_search(query="HISQ lattice QCD highly improved staggered quark", max_results=10, sort_by="relevance")
arxiv_search(query="HISQ lattice QCD highly improved staggered quark", max_results=10, sort_by="relevance") is not available in this environment.
HISQ, the **Highly Improved Staggered Quark** action, is a staggered-fermion discretization in lattice QCD designed to reduce taste-symmetry breaking and suppress leading cutoff effects. In the formulations summarized across scaling, spectroscopy, renormalization, and algorithmic studies, HISQ combines two levels of link smearing, reunitarization, and a Naik correction so that leading discretization errors are pushed down to \(O(\alpha_S a^2, a^4)\), while taste-breaking splittings are reduced to an additional small \(O(\alpha_S^2 a^2)\) level [1012.1265]. Numerical comparisons with the earlier asqtad action consistently report substantially smaller lattice artifacts at the same lattice spacings and quark masses, and subsequent work has deployed HISQ in flavor physics, spectroscopy, scale setting, topological observables, flux-tube studies, nonperturbative renormalization, and large-scale solver development [1004.0342].

## 1. Definition and design objectives

The HISQ action was introduced to address two persistent limitations of staggered fermions. The first is **taste-symmetry breaking**, which causes the sixteen staggered pion tastes to become nondegenerate at finite lattice spacing. The second is the presence of leading \(O(a^2)\) errors in the fermion dispersion relation, which are especially important for charm quarks. The HISQ design goal is therefore twofold: to make the pion tastes nearly degenerate even at finite \(a\), and to improve the fermion action sufficiently that charm can be treated relativistically on fine enough lattices [1004.0342].

In the modern lattice-QCD usage summarized in the gradient-flow and flux-tube studies, HISQ is not an isolated fermion kernel but part of a broader lattice action. One couples the HISQ Dirac operator to a **tree-level Symanzik-improved gauge action**, so that
\[
S_{\rm HISQ} \;=\; S_{\rm gauge}^{\rm (Sym)}\;+\;\sum_{f=u,d,s,c} \bar\chi_f\,\bigl(D^{\rm HISQ}[U^{\rm (HISQ)}]+m_f\bigr)\,\chi_f\,.
\]
For the HISQ/tree setup used in \((2+1)\)-flavor flux-tube simulations, the gauge part is the tree-level improved Symanzik action with coefficients \(c_0=5/3\) and \(c_1=-1/12\) at tree level [1411.0068; 1701.03371].

A recurrent comparison point is the older **asqtad** action. Across the scaling studies and phenomenological applications, HISQ is described as a further \(O(a^2)\)-improved version of staggered fermions that reduces taste splittings by roughly a factor of two to three relative to asqtad, and in some summaries by roughly an order of magnitude depending on the observable and convention being emphasized [1212.0613; 1511.06547]. This suggests that the phrase “HISQ improvement” is best understood not as one isolated modification, but as a coordinated suppression of taste exchange, dispersion errors, and generic cutoff effects.

## 2. Fermion operator and link construction

The HISQ Dirac operator is built in stages. In the scaling study of Bazavov et al., the massless part is written as
\[
(/D_{\rm HISQ})_{x,y}
= \sum_{\mu}\bigl[\,\delta_{x+\hat\mu,y}\,X_\mu(x)\;-\;\delta_{x-\hat\mu,y}\,X_\mu^\dagger(x-\hat\mu)\bigr]
\;+\;(1+\epsilon_N)\sum_{\mu}\bigl[\,\delta_{x+3\hat\mu,y}\,W_\mu(x)\,W_\mu(x+\hat\mu)\,W_\mu(x+2\hat\mu)\;-\;\mathrm{h.c.}\bigr]\,,
\]
where \(X_\mu(x)\) is the twice-smeared link entering the one-hop term and \(W_\mu(x)\) is the once-smeared, unitarized link entering the Naik term [1004.0342].

The construction proceeds through four steps: **Fat7 smearing** of the original gauge links \(U\to V\); **projection** \(V\to W\) onto \(U(3)\) by polar decomposition or SVD; **Asqtad-style smearing** \(W\to X\); and addition of a **Naik term** with a mass-dependent coefficient \(\epsilon_N\) for the charm quark [1004.0342]. The same sequence is summarized elsewhere as Fat7 smearing, reunitarization, Lepage correction, and Naik improvement [1211.0751; 1703.05552].

The explicit Fat7 coefficients quoted in the scaling study are: 1-link \(c=1/8\), 3-staple \(c=1/16\), 5-staple \(c=1/64\), and 7-staple \(c=1/384\), with all rotations and reflections included. After projection, the second smearing repeats the same staple structure and adds a negative Lepage term with \(c=-1/8\), while the Naik term carries \(c=-(1+\epsilon_N)/24\) [1004.0342]. In the HISQ/tree implementation described for flux-tube calculations, tadpole improvement uses
\[
u_0 = \langle \mathrm{Re}\,\mathrm{Tr}\,P_{\mu\nu}\rangle^{1/4}
\]
from the plaquette, and reunitarization is carried out by polar decomposition, projecting first to \(U(3)\) and then to \(SU(3)\) [1701.03371].

For light and strange quarks, \(\epsilon_N=0\). For the charm quark, the tree-level coefficient is expanded as
\[
\epsilon_N
= -\tfrac{27}{40}(am_c)^2+\tfrac{327}{1120}(am_c)^4-\tfrac{15607}{268800}(am_c)^6-\tfrac{73697}{3942400}(am_c)^8\,,
\]
which is chosen so that heavy-quark dispersion relations are accurate through higher order in \(am_c\) [1004.0342]. In related summaries this same role is described as tuning \(c_{\rm Naik}\) or \((1+\varepsilon)\) so that the tree-level kinetic mass equals the pole mass, or so that the dispersion relation is improved through \(O((ap)^4)\) [1212.0613; 1211.6966].

## 3. Error structure, taste symmetry, and scaling

The central theoretical claim attached to HISQ is that the two-stage smearing plus reunitarization suppresses taste-exchange interactions much more effectively than earlier staggered formulations. In the dynamical-HISQ overview, the leading discretization errors are summarized as \(O(\alpha_S a^2, a^4)\), and the remaining taste-breaking splittings as \(O(\alpha_S^2 a^2)\) [1012.1265]. In the flux-tube summary, the tree-level Symanzik gauge action removes \(O(a^2)\) errors in the pure-gauge sector, while Fat7 smearing, projection, and the Naik term cancel the leading taste-breaking \(O(\alpha_s a^2)\) artifacts of staggered fermions; the sequence “Fat7 \(\to\) project \(\to\) Asq-tone” pushes taste splittings down to \(O(\alpha_s^2 a^2)\) or smaller [1701.03371].

The standard scaling relation for non-Goldstone pion tastes is
\[
M_\xi^2 - M_{\rm G}^2
\;=\; C_\xi\;\alpha_S^2\;a^2 \;+\; O(a^4)\,,
\]
or, in the form used in the scaling study,
\[
\Delta M^2\cdot r_1^2 = (M^2_{\rm non\text{-}G}-M^2_G)\cdot r_1^2\,.
\]
Numerically, HISQ reduces these taste splittings by \(\simeq 3\times\) relative to asqtad at the same lattice spacing, in excellent agreement with the expected \(a^2\alpha_s^2\) scaling [1004.0342; 1012.1265]. In the correlator-fit summary, the largest non-Goldstone pion at \(a\approx0.12\) fm is \(\sim70\) MeV above the Goldstone, compared to \(\sim200\) MeV with asqtad [1212.0613].

This improved taste symmetry is accompanied by smaller cutoff effects in hadronic quantities. For light-light and heavy-light masses, the HISQ points lie very close to the asqtad continuum curve, and “a HISQ result at \(a\) matches an asqtad result at \(\simeq 2a/3\)” [1004.0342]. The same pattern is reported for pseudoscalar decay constants and topological susceptibility, with the latter improved by about a factor \((2/3)^2\) in \(a\) [1004.0342]. The dynamical-HISQ program likewise reports flatter \(a^2\) dependence in \(r_1 M_\rho\) and \(r_1 M_N\), with slopes reduced by roughly a factor of two relative to asqtad [1012.1265].

A common misconception is that HISQ only improves taste splittings in the light sector. The heavy-quark tuning of the Naik term is an equally central part of the construction. This is why later work could extend the **heavy-HISQ** method to valence heavy masses ranging from the charm quark up to very nearly the physical \(b\)-quark on the finest lattices, while maintaining controlled discretization effects in joint continuum and heavy-mass fits [2502.06713].

## 4. Ensembles, scale setting, and lines of constant physics

HISQ simulations are typically organized along a **line of constant physics**. In the flux-tube study with \((2+1)\) flavors, this means choosing \((m_l,m_s)\) at each gauge coupling \(\beta\) so that hadron masses remain at their physical values, with \(m_s\) fixed to its physical value and \(m_l=m_s/20\), implying \(M_\pi \simeq 160\) MeV [1701.03371]. The stated purpose is to guarantee comparability of flux-tube profiles measured at different lattice spacings and to ensure that quark masses track the physical point as \(a\to0\) [1701.03371].

The early dynamical-HISQ scaling studies used four flavors of dynamical HISQ quarks at lattice spacings approximately 0.15, 0.12, and 0.09 fm, with \(m_l=0.2\,m_s\), \(m_s\approx m_s^{\rm phys}\), and \(m_c\approx m_c^{\rm phys}\) [1004.0342]. Later MILC/Fermilab HISQ ensembles extended to \(2+1+1\) flavors, lattice spacings from 0.15 fm down to 0.06 fm, and light sea-quark masses from \(\approx0.2\,m_s\) down to the physical point [1212.0613]. Gradient-flow analyses on these ensembles explicitly covered four lattice spacings from \(a\approx0.15\) fm down to \(0.06\) fm and both physical and unphysical quark masses [1411.0068].

Several scale-setting prescriptions recur in the HISQ literature. One is the static-potential quantity \(r_1\), defined through
\[
r_1^2 F(r_1)=-1\,,
\]
with \(F(r)=-\partial V/\partial r\) [1004.0342]. Another is the decay constant of a fictitious “unmixed \(s\bar s\)” pseudoscalar, \(f_{ss}\), obtained by tuning the valence mass until the pseudoscalar has the physical ratio \(M_{ss}/f_{ss}\) [1004.0342]. A later review of dynamical-HISQ simulations gives representative lattice spacings from both \(r_1\) and \(f_{ss}\), and emphasizes that differences \(a(f_{ss})-a(r_1)\) vanish as \(a^2\to0\) [1012.1265].

Gradient-flow scales provide a third approach. Using Symanzik flow and the cloverleaf definition of \(\langle E\rangle\), the scales \(t_0\) and \(w_0\) are defined by
\[
t^2\,\langle E(t)\rangle\Big|_{t=t_0}=0.3,\qquad
t\,\frac{d}{dt}\Bigl[t^2\langle E(t)\rangle\Bigr]\Big|_{t=w_0^2}=0.3\,.
\]
On MILC \(N_f=2+1+1\) HISQ ensembles, the reported preliminary physical values are \(\sqrt{t_0}=0.1422(7)\) fm and \(w_0=0.1732(10)\) fm [1411.0068]. This suggests that HISQ calculations have supported a transition from purely potential-based scale setting to a mixed ecosystem of potential, pseudoscalar, and flow-based standards.

## 5. Phenomenology and nonperturbative applications

HISQ has been used extensively in precision flavor physics. In two-point correlator analyses on nineteen \(2+1+1\)-flavor HISQ ensembles, simultaneous fits with the full covariance matrix and Gaussian constraints were used to extract light, strange, and charmed pseudoscalar observables, supporting preliminary results
\[
f_D = 209(4)_{\rm stat}\,(^{+3}_{-1})_{\rm syst}\ \mathrm{MeV},\qquad
f_{D_s} = 248(3)_{\rm stat}\,(^{+2}_{-1})_{\rm syst}\ \mathrm{MeV},
\]
together with \(f_{D_s}/f_D = 1.188(25)\) and quark-mass ratios \(m_s/m_l=27.2(1.5)\), \(m_c/m_s=11.8(2)\) [1212.0613].

In kaon semileptonic decays, HISQ valence quarks were first used on asqtad \(N_f=2+1\) ensembles and then on full HISQ \(N_f=2+1+1\) ensembles including physical light-quark masses. The partially quenched staggered-chiral analysis of \(f_+^{K\pi}(0)\) reported projected total uncertainties of \(0.35\%\)–\(0.50\%\) on the asqtad-sea data and anticipated \(0.2\%\)–\(0.3\%\) on full HISQ data because sea-strange tuning and dynamical charm are correct by construction and taste splittings are another factor of two smaller in HISQ sea [1211.0751]. In semileptonic \(B\) and \(B_s\) work with NRQCD heavy quarks and HISQ light valence quarks, the HISQ component was specifically used to suppress taste-breaking and reduce renormalization uncertainties in form-factor calculations over a range of momentum transfer [1210.6992; 1211.6966].

Spectroscopy and current normalization provide another major application. In the nonperturbative comparison of clover and HISQ strange-quark correlators, local vector and axial currents made from HISQ, clover, and mixed HISQ–clover fields were normalized nonperturbatively, and the physical \(\phi\)-meson quantities were reported as
\[
m_{\phi} = 1.023(6)\ \mathrm{GeV},\qquad
f_{\phi} = 0.238(3)\ \mathrm{GeV},
\]
in good agreement with experiment, using HISQ valence quarks on \(N_f=2+1+1\) ensembles with light sea masses down to the physical point [1703.05552]. The same study emphasizes that the remnant \(U(1)\) chiral symmetry of HISQ provides an absolutely normalized decay constant for the Goldstone-taste pseudoscalar through the PCAC relation [1703.05552].

Beyond hadron structure and flavor physics, HISQ has been used in gauge-field diagnostics and vacuum-structure studies. In the flux-tube calculation with \((2+1)\) HISQ fermions, Monte Carlo simulations using the HISQ/tree action in publicly available MILC code were used to investigate the transverse profile of the chromoelectric field generated by a quark-antiquark pair [1701.03371]. In nonperturbative renormalization, HISQ bilinears have been analyzed in RI/MOM and RI/SMOM schemes on \(n_f=2+1+1\) HISQ ensembles, where the SMOM setup was reported to suppress infrared contamination strongly relative to exceptional MOM kinematics [1511.06547].

## 6. Solvers, eigenmodes, and terminological scope

The computational cost of HISQ inversions has driven substantial algorithmic work. Standard propagator calculations typically use mixed-precision CG with even–odd preconditioning, but this suffers from critical slowing down as the light-quark mass approaches its physical value [2502.00152]. On a \(144^3\times288\), \(a=0.04\) fm ensemble at the physical light-quark mass, the baseline undeflated CG reported \(33\,003\) iterations and \(47.8\) s per solve at residual \(10^{-8}\); deflation with \(1024\) eigenvectors reduced this to \(5\,593\) iterations and \(8.72\) s, while \(2048\) eigenvectors reduced it to \(2\,909\) iterations and \(4.79\) s, a \(10.0\times\) speedup over CG [2502.00152].

Earlier numerical experiments already showed the same qualitative behavior. On a \(24^3\times48\), \(m_\pi=220\) MeV HISQ ensemble, exact deflation reduced the CG iteration count from \(120\) to \(70\) with \(50\) deflated modes and to \(40\) with \(100\) modes, corresponding to speed-ups of \(1.7\) and \(3.0\) [1710.07219]. The same study emphasized that eigenpairs must be converged to a residual at least an order of magnitude below the target solve residual; otherwise the CG residual levels off [1710.07219].

For eigenmode calculations, studies of the Hermitian positive-definite operator \(D^\dagger D\) with the HISQ Dirac operator compared Implicitly Restarted Lanczos, Thick-Restart Lanczos, and Block Lanczos as implemented in **Grid** and **QUDA**. Chebyshev polynomial filtering with \((\alpha,\beta)\approx(0.05,25)\) and degree \(p\approx100\) was reported as optimal for extracting \(O(10^2)\) modes, QUDA’s TRL outperformed non-restarted Lanczos in all tested \((n,p)\), and Block Lanczos with Split-Grid was nearly competitive with the unblocked code while reducing communication overhead [2201.03755]. This suggests that the practical history of HISQ is inseparable from the development of libraries such as **MILC**, **QUDA**, and **Grid**.

Finally, the acronym **HISQ** is not unique across arXiv. In lattice QCD it denotes **Highly Improved Staggered Quark** and refers to the fermion action discussed throughout this literature [1004.0342]. In a 2025 quantum-control architecture paper, however, **HISQ** denotes a universal instruction set defined as an extension of RV32I and paired with a booking-based synchronization protocol in distributed control hardware [2509.04798]. The coexistence of these two usages is terminological rather than conceptual. In the lattice-QCD literature, HISQ retains its established meaning: a staggered-quark discretization built from Fat7 smearing, reunitarization, Lepage correction, and Naik improvement, paired in practice with improved gauge actions and used for high-precision calculations across a broad range of observables [1012.1265; 1701.03371].

Source: https://www.emergentmind.com/topics/hisq