---
title: Asqtad-Improved Staggered Fermions
url: https://www.emergentmind.com/topics/asqtad-improved-staggered-fermions
type: topic
---

# Asqtad-Improved Staggered Fermions

Asqtad-improved staggered fermions are a class of lattice fermion discretizations employed in large-scale lattice QCD simulations to achieve full tree-level $\mathcal O(a^2)$ Symanzik improvement and strong suppression of taste-exchange interactions. The Asqtad program combines the staggered formalism—retaining a single component $\chi(x)$ per site representing four continuum tastes per flavor—with a set of gauge-link improvements, including Fat7 smearing, Lepage double-staple corrections, the Naik three-link dispersion correction, and tadpole improvement. Used extensively by the MILC Collaboration, Asqtad actions push lattice artifacts and taste-symmetry breaking down to the few-percent level at moderate lattice spacings while preserving computational efficiency.

## 1. Formal Structure of the Asqtad Action

The Asqtad-improved staggered action is formulated as
\[
S_F = a^4\sum_x \bar\chi(x)\Bigl[D^\text{Asqtad}(m)\Bigr]\chi(x)
\]
where the Dirac operator $D^\text{Asqtad}$ includes the following terms, all tadpole-improved:
- The one-link “thin” staggered derivative (naive term).
- Fat7 seven-link smeared derivative, suppressing high-momentum gluon exchange.
- Lepage five-link planar staple subtraction to remove $\mathcal O(a^2)$ taste breaking arising from Fat7 smearing.
- Naik three-link straight “hop” term, which corrects the fermion dispersion relation to $\mathcal O(a^2)$ [1204.3866, 1012.4484].

Explicitly, in the one-component notation:
\[
D_\mu^\text{Asqtad}\chi(x) = \eta_\mu(x)\left[
\rho_1 U_\mu^\text{fat7}(x)\chi(x+\hat\mu)
+ \rho_3 U_\mu^\text{lep}(x)\chi(x+\hat\mu)
+ \rho_N U_\mu^\text{Naik}(x)\chi(x+3\hat\mu)
-\text{(hermitian conj.)}
\right]
\]
with $\rho_N \approx -1/24$ and tadpole improvement factor $u_0$ applied to all link terms. $U_\mu^\text{fat7}$ and $U_\mu^\text{lep}$ implement multi-link staple averages and planar double-staple subtractions, respectively, while $U_\mu^\text{Naik}$ is a straight product of three links.

The gauge action, consistently, is the one-loop tree-level Symanzik-improved plaquette plus rectangle (“Lüscher-Weisz”) form, tadpole-improved:
\[
S_G = \beta \sum_x\left[ c_0 \sum_{\mu<\nu} (1-\tfrac{1}{3}\mathrm{Re\,tr\,}P_{\mu\nu}) + c_1\sum_{\mu\neq\nu}(1-\tfrac{1}{3}\mathrm{Re\,tr\,}R_{\mu\nu}) \right]
\]
with $c_1=-1/12$, $c_0=5/3$, $u_0$ measured as the average plaquette.

## 2. Motivation and Design of Improvement Terms

Each Asqtad improvement is chosen to target dominant sources of taste-symmetry breaking and discretization error:
- **Fat7 smearing** reduces coupling to gluons with $k\sim\pi/a$ that induce taste mixing, leading to suppression of tree-level taste-exchange interactions.
- **Lepage correction** subtracts residual $\mathcal O(a^2)$ taste-breaking vertices that fat7 over-smears into the action at high momentum.
- **Naik term** cancels the tree-level $\mathcal O(a^2)$ error in the free-fermion dispersion, restoring the correct energy-momentum relation $E^2 = p^2 + m^2+\mathcal O(a^4)$.
- **Tadpole improvement** absorbs large perturbative corrections by dividing all links by $u_0$, ensuring perturbative renormalization matching factors remain small and stable [1204.3866, 1012.4484, 1501.01991, 1112.4384].

Collectively, this design ensures no tree-level $\mathcal O(a^2)$ taste violations, achieves near-ideal dispersion, and pushes taste breaking down to $\mathcal O(\alpha_s a^2)$ at one loop.

## 3. Numerical Implementation and Performance

Asqtad improvement is implemented in large dynamical QCD simulations, notably by the MILC collaboration:
- **Ensembles**: 2+1 flavor QCD with $a\approx0.15,0.12,0.09,0.06$, and $0.045$ fm, volumes from $20^3\times64$ up to $48^3\times144$.
- **Quark masses**: $m_l/m_s$ ratios down to $1/20$, with physical strange mass matched to experimental meson masses via chiral extrapolation.
- **Action coefficients**: $\rho_N \rightarrow -1/24$, Lepage coefficient $-1/24$, fat7 staple weights fixed at tree-level, explicit tadpole improvement throughout [1204.3866, 1012.4484, 1501.01991, 1112.4384].

Performance metrics:
- Largest non-Goldstone pion mass-squared splitting drops from $\sim0.36\,\mathrm{GeV}^2$ (at $a\approx0.15\,\mathrm{fm}$) to $<0.01\,\mathrm{GeV}^2$ ($a\approx0.045\,\mathrm{fm}$).
- Residual taste violations on the finest lattices are $\sim(100-150\,\mathrm{MeV})^2$, a reduction by factors $\sim3-4$ from coarsest lattices.
- Light-quark discretization error in decay-constant ratios $<0.3\%$ for $f_{D_s}/f_D$, $f_{B_s}/f_B$, and $f_B/f_D$ [1501.01991].

## 4. Suppression of Taste-Symmetry Breaking

The central achievement of Asqtad improvement is the suppression of taste-symmetry violations. In naive staggered fermions, exchange of hard gluons produces sizable taste mixings, splitting the sixteen pion tastes into multiplets separated by $\mathcal O(\alpha_s a^2)$ in mass squared. Asqtad smearing and operator design remove all such splitting at tree level, leaving only one-loop residuals. Empirically, the ratio of taste-splitting between $a\approx0.12$ fm and $a\approx0.09$ fm lattices is measured as $0.39(1)$, matching the expected $O(a^2)$ scaling [1112.4384].

Table: Scaling of Taste-Splitting in Asqtad Ensembles

| Lattice Spacing $a$ (fm) | Max $\Delta M_\pi^2$ ($\mathrm{GeV}^2$) | Taste Splitting Reduction |
|--------------------------|-----------------------------------|--------------------------|
| 0.15                     | $\sim0.36$                        | Baseline                 |
| 0.09                     | $\sim0.04$                        | $\sim 9\times$ lower     |
| 0.045                    | $<0.01$                           | $>35\times$ lower        |

## 5. Perturbative and Nonperturbative Renormalization

Renormalization factors for bilinear operators in the Asqtad scheme have been computed both in one-loop perturbation theory and via nonperturbative RI-MOM methods:
- **One-loop matching**: Perturbative corrections for bilinear operators are reduced compared to naive staggered, with spread of matching coefficients in the scalar sector $8.9$ (Wilson glue) and $8.7$ (Symanzik glue), much smaller than $\sim42$ for naive sparse links [1004.4039].
- **Nonperturbative results**: NPR factors $Z_q^\text{RI-MOM}(2\,\text{GeV})=1.0764(44)$, $Z_m^\text{RI-MOM}(2\,\text{GeV})=1.246(15)$, scalar channel $Z_{S\otimes P}=1.079(18)$, vector $Z_{V\otimes V}=0.982(11)$, tensor $Z_{T\otimes T}=1.293(16)$ on MILC coarse ensembles [1310.4269].
- **Consistency**: Vector and tensor channels show good agreement between NPR and PT (few percent); scalar channels may depart by up to $0.15$–$0.20$ at $2\,\text{GeV}$, indicating the necessity of nonperturbative matching for precision work [1306.3881].

## 6. Applications in Lattice QCD and Quantum Simulation

Asqtad ensembles have underpinned key lattice calculations:
- **Hadron Spectrum**: Precision work in $B$- and $D$-meson decay constants, chiral condensates, and flavor-singlet masses ($\eta$, $\eta'$) [1501.01991, 1112.4384].
- **Finite-Temperature QCD**: Determination of crossover temperature $T_c=(164\pm6)$ MeV (continuum-extrapolated), with chiral condensate, susceptibility, and Polyakov loop pointing to smooth crossover; taste breaking produces a $5$–$10$ MeV shift in $T_c$ relative to HISQ [1012.4484].
- **Nucleon Matrix Elements**: The intrinsic strangeness $\langle N| \bar{s}s | N \rangle=0.637(55)_{\text{stat}}(74)_{\text{sys}}$ and preliminary charm results; the hybrid three-point/Feynman–Hellman method leverages the per-timeslice condensate structure of Asqtad [1204.3866].
- **Quantum Simulation**: The Asqtad-inspired Hamiltonian has been mapped to qubit circuits, demonstrating that vector-meson mass extrapolation is $30$–$50\%$ closer to the continuum for improved actions versus naive Kogut–Susskind, at the cost of increased gate counts scaling as $O(d)$ in dimension [2402.04317].

## 7. Comparison with Other Improved Staggered Actions and Future Prospects

The HISQ (“highly improved staggered quark”) action extends Asqtad improvement by introducing additional reunitarized smearing and further corrections to cancel one-loop $\mathcal O(a^2)$ taste-exchange at tree level. HISQ yields taste splittings $\sim$ one-third of Asqtad at fixed $a$ and enables physical quark masses at coarser lattices. However, at the time of key calculations, essential quantities for improved-hybrid methods (per-timeslice condensates, full propagators) were not available for HISQ ensembles [1204.3866].

A plausible implication is that future large-scale dynamical QCD projects seeking sub-percent accuracy on hadronic observables will migrate towards HISQ or related actions, while Asqtad remains a benchmark for systematic control of taste exchange and is well-matched for precision studies where link and operator improvement are critical.

---

**References:**  
[1204.3866], [1012.4484], [1501.01991], [1112.4384], [1004.4039], [1310.4269], [1306.3881], [2402.04317]

Source: https://www.emergentmind.com/topics/asqtad-improved-staggered-fermions