== Glossary == <> The material below is from the first edition of the review. The updated glossary from the present edition can be downloaded as a [[attachment:Media/FLAG_glossary.pdf|pdf-file|&do=get]]. === Gauge actions === The simplest and most widely used discretisation of the Yang-Mills part of the QCD action is the Wilson plaquette action <>: $S_\text{G} = \beta\sum_{x} \sum_{\mu<\nu}\Big(1-\frac{1}{3}\text{Re Tr}\,W_{\mu\nu}^{1\times1}(x)\Big),$ where the plaquette, $W_{\mu\nu}^{1\times1}(x)$, is the product of link variables around an elementary square of the lattice, i.e. $W_{\mu\nu}^{1\times1}(x) \equiv U_\mu(x)U_\nu(x+a\hat{\mu})U_\mu(x+a\hat{\nu})^{-1} U_\nu(x)^{-1}.$ This expression reproduces the Euclidean Yang-Mills action in the continuum up to corrections of order $a^2$. There is a general formalism, known as the ``Symanzik improvement programme`` <> <>, which is designed to cancel the leading lattice artefacts, such that observables have an accelerated rate of convergence to the continuum limit. The improvement programme is implemented by adding higher-dimensional operators, whose coefficients must be tuned appropriately in order to cancel the leading lattice artefacts. The effectiveness of this procedure depends largely on the method with which the coefficients are determined. The most widely applied methods (in ascending order of effectiveness) include perturbation theory, tadpole-improved (partially resummed) perturbation theory, renormalisation group methods, and the non-perturbative evaluation of improvement conditions. In the case of Yang-Mills theory, the simplest version of an improved lattice action is obtained by adding rectangular 1$\times$2 loops to the plaquette action, i.e. $S_\text{G}^\text{imp} = \beta\sum_{x}\left\{ c_0\sum_{\mu<\nu}\Big(1-\frac{1}{3}\text{Re Tr}\,W_{\mu\nu}^{1\times1}(x)\Big)+c_1\sum_{\mu,\nu} \Big(1-\frac{1}{3}\text{Re Tr}\,W_{\mu\nu}^{1\times2}(x)\Big) \right\},$ where the coefficients $c_0, c_1$ satisfy the normalisation condition $c_0+8c_1=1$. The ''Symanzik-improved'' <>, ''Iwasaki'' <>, and ''DBW2'' <><> actions are all defined through the equation above via particular choices for $c_0, c_1$. Details are listed in the following Table, together with the abbreviations used in the summary tables. <
> <
> ||<:20%#eeeeee>'''Abbrev.'''||<:15%#eeeeee>'''$c_1$'''||<:#eeeeee>'''Description''' || || Wilson || 0 || Wilson plaquette action || || tlSym || -1/12 || tree-level Symanzik-improved gauge action || || tadSym || variable || tadpole Symanzik-improved gauge action || || Iwasaki || −0.331 || Renormalisation group improved (“Iwasaki”) action|| || DBW2 || −1.4088 || Renormalisation group improved (“DBW2”) action || The leading lattice artefacts are $O(a^2)$ or better for all discretisations. ==== References ==== <> === Light-quark actions === If one attempts to discretise the quark action, one is faced with the fermion doubling problem: the naive lattice transcription produces a 16-fold degeneracy of the fermion spectrum. ==== Wilson fermions ==== Wilson's solution to the doubling problem is based on adding a dimension-5 operator which removes the doublers from the low-energy spectrum. The Wilson-Dirac operator for the massless case reads \cite{Wilson:1974sk} $D_\text{w} = \frac{1}{2}\gamma_\mu(\nabla_\mu+\nabla_\mu^*)+a\nabla_\mu^*\nabla_\mu,$ where $\nabla_\mu,\,\nabla_\mu^*$ denote lattice versions of the covariant derivative. Adding the Wilson term, $a\nabla_\mu^*\nabla_\mu$, results in an explicit breaking of chiral symmetry even in the massless theory. Furthermore, the leading order lattice artefacts are of order $a$. With the help of the Symanzik improvement programme, the leading artefacts can be cancelled by adding the so-called ''Clover'' or Sheikholeslami-Wohlert (SW) term. The resulting expression in the massless case reads $D_\text{sw} = D_\text{w}+\frac{ia}{4}\,c_\text{sw}\sigma_{\mu\nu}\widehat{F}_{\mu\nu},$ where $\sigma_{\mu\nu}=\frac{i}{2}[\gamma_\mu,\gamma_\nu]$, and $\widehat{F}_{\mu\nu}$ is a lattice transcription of the gluon field strength tensor $F_{\mu\nu}$. Provided that the coefficient $c_\text{sw}$ is suitably tuned, observables computed using $D_\text{sw}$ will approach the continuum limit with a rate proportional to~$a^2$. Chiral symmetry remains broken, though. The coefficient $c_\text{sw}$ can be determined perturbatively at tree-level (tree-level impr., $c_\text{sw} = 1$ or tlSW in short), via a mean field approach <> (mean-field impr. or mfSW) or via a non-perturbative approach <> (non-perturbativley impr. or npSW). Finally, we mention ''twisted mass QCD'' as a method which was originally designed to address another problem of Wilson's discretisation: the Wilson-Dirac operator is not protected against the occurrence of unphysical zero modes, which manifest themselves as ''exceptional'' configurations. They occur with a certain frequency in numerical simulations with Wilson quarks and can lead to strong statistical fluctuations. The problem can be cured by introducing a so-called ''chirally twisted'' mass term, after which the fermionic part of the QCD action in the continuum assumes the form <> $S_\text{F}^\text{tm;cont} = \int d^4{x}\, \bar \psi(x)(\gamma_\mu D_\mu + m + i\mu_\text{q}\gamma_5\tau^3)\psi(x).$ Here, $\mu_\text{q}$ is the twisted mass parameter, and $\tau^3$ is a Pauli matrix. The standard action in the continuum can be recovered via a global chiral field rotation. The lattice action of twisted mass QCD (tmWil) for $N_f=2$ flavours is defined as $S_\text{F}^\text{tm}[U,\bar \psi,\psi] = a^4\sum_{x\in\Lambda_\text{E}}\bar \psi(x)(D_\text{w}+m_0+i\mu_\text{q}\gamma_5\tau^3)\psi(x).$ Although this formulation breaks physical parity and flavour symmetries, is has a number of advantages over standard Wilson fermions. In particular, the presence of the twisted mass parameter $\mu_\text{q}$ protects the discretised theory against unphysical zero modes. Another attractive feature of twisted mass lattice QCD is the fact that the leading lattice artefacts are of order $a^2$ without the need to add the Sheikholeslami-Wohlert term <>. Although the problem of explicit chiral symmetry breaking remains, the twisted formulation is particularly useful to circumvent some of the problems that are encountered in connection with the renormalization of local operators on the lattice, such as those required to determine $B_\text{K}$. ==== Staggered fermions ==== An alternative procedure to deal with the doubling problem is based on so-called ''staggered'' or Kogut-Susskind (KS) fermions <>. Here the degeneracy is only lifted partially, from 16 down to~4. It has become customary to refer to these residual doublers as ''tastes'' in order to distinguish them from physical flavours. At order~$a^2$ different tastes can interact via gluon exchange, thereby generating large lattice artefacts. The improvement programme can be used to suppress taste-changing interactions, leading to ''improved staggered fermions'', with the so-called ''Asqtad'' action as one of its widely used versions <>. The standard procedure to remove the residual doubling of staggered quarks (''four tastes per flavour'') is to take fractional powers of the quark determinant in the QCD functional integral. This is usually referred to as the ''fourth root trick''. The validity of this procedure has not been rigorously proven so far. In fact, it has been questioned by several authors, and the issue is still hotly debated (for both sides of the argument see the reviews in refs. <> <> <> <> <>. ==== Ginsparg-Wilson fermions ==== Fermionic lattice actions, which do not suffer from the doubling problem whilst preserving chiral symmetry go under the name of ''Ginsparg-Wilson fermions''. In the continuum the massless Dirac operator anti-commutes with $\gamma_5$. At non-zero lattice spacing chiral symmetry can be realised even if this condition is relaxed according to <> <> $\left\{D,\gamma_5\right\} = aD\gamma_5 D,$ which is now known as the Ginsparg-Wilson relation <>. A lattice Dirac operator which satisfies the GW relation can be constructed in several ways. The ''domain wall'' construction proceeds by introducing a fifth dimension of length $N_5$ and coupling the fermions to a mass defect (i.e. a negative mass term) <>. The five-dimensional action can be constructed such that modes of opposite chirality are trapped at the four dimensional boundaries in the limit of an infinite extent of the extra dimension <>. In any real simulation, though, one has to work with a finite value of $N_5$, so that the decoupling of chiral modes is not exact. This leads to a residual breaking of chiral symmetry, which, however, is exponentially suppressed. A doubler-free, (approximately) chirally symmetric quark action can thus be realised at the expense of simulating a five-dimensional theory. The so-called ''overlap'' or Neuberger-Dirac operator can be derived from the domain wall formulation <>. It acts in four space-time dimensions and is, in its simplest form, defined by $D_{\rm N} = \frac{1}{\bar a} \left( 1-\frac{A}{\sqrt{A^\dagger{A}}}\right), \quad A=1+s-aD_{\rm w},\quad \bar a=\frac{a}{1+s},$ where $D_{\rm w}$ is the massless Wilson-Dirac operator, and $|s|<1$ is a tunable parameter. The overlap operator $D_\text{N}$ removes all doublers from the spectrum, and can easily be shown to satisfy the Ginsparg-Wilson relation. The occurrence of an inverse square root in $D_\text{N}$ renders the application of $D_\text{N}$ in a computer program potentially very costly, since it must be implemented using, for instance, a polynomial approximation. The third example of an operator which satisfies the Ginsparg-Wilson relation is the so-called fixed-point action <> <>. This construction proceeds via a renormalisation group approach. A related formalism are the so-called ''chirally improved'' fermions <>. ==== Smearing ==== A simple modification which can help improve the action as well as the computational performance is the use of smeared gauge fields in the covariant derivatives of the fermionic action. Any smearing procedure is acceptable as long as it consists of only adding irrelevant (local) operators. Moreover, it can be combined with any discretisation of the quark action. The ''Asqtad'' staggered quark action mentioned above <> is an example which makes use of so-called ''Asqtad'' smeared (or ''fat'') links. Another example is the use of n-HYP smeared <> or n-stout smeared <> gauge links in the tree-level clover improved discretisation of the quark action, denoted by ''n-HYP tlSW'' and ''n-stout tlSW'' in the following. In the table below we summarise the most widely used discretisations of the quark action and their main properties together with the abbreviations used in the summary tables. Note that in order to maintain the leading lattice artefacts of the actions as given in the table in non-spectral observables (like operator matrix elements) the corresponding non-spectral operators need to be improved as well. ||<:#eeeeee>'''Abbrev.'''||<:#eeeeee>'''Discretisation'''||<:#eeeeee>'''Leading lattice artefacts'''||<:#eeeeee>'''Chiral symmetry'''||<:#eeeeee>'''Remarks'''|| || Wilson || Wilson || $O(a)$ || broken || -- || || tmWil ||Twisted Mass Wilson || $O(a^2)$ at maximal twist || broken || flavour symmetry breaking: $(M_\text{PS}^{0})^2-(M_\text{PS}^\pm)^2\sim O(a^2)$ || || tlSW || Sheikholeslami-Wohlert || $O(g^2 a)$ || broken || tree-level improved, $c_\text{sw}=1$ || || n-HYP tlSW || Sheikholeslami-Wohlert || $O(g^2 a)$ || broken || tree-level improved, $c_\text{sw}=1$, n-HYP smeared gauge links || || n-stout tlSW || Sheikholeslami-Wohlert || $O(g^2 a)$ || broken || tree-level improved, $c_\text{sw}=1$, n-stout smeared gauge links || || mfSW || Sheikholeslami-Wohlert || $O(g^2 a)$ || broken || mean-field improved || || npSW || Sheikholeslami-Wohlert || $O(a^2)$ || broken || non-perturbatively improved || || KS || Staggered || $O(a^2)$ || U(1)$\times$U(1) subgroup unbroken || rooting for $N_f<4$ || || Asqtad || Staggered || $O(a^2)$ || U(1)$\times$U(1) subgroup unbroken || rooting for $N_f<4$, Asqtad smeared gauge links || || DW || Domain Wall || asymptotically $O(a^2)$ || remnant breaking exponentially suppressed || exact chiral symmetry and $O(a)$ improvement only in the limit $L_s\rightarrow \infty$ || || overlap || Neuberger || $O(a^2)$ || exact || -- || ==== References ==== <> === Matching and running === ||<:#eeeeee>'''Abbrev.''' ||<:#eeeeee>'''Description''' || || RI || Regularisation-independent momentum subtraction scheme || || SF || Schroedinger functional scheme || || PT1L || matching/running computed in perturbation theory at one loop || || PT2L || matching/running computed in perturbation theory at two loop || === Chiral extrapolation === Symanzik’s framework can be combined with Chiral Perturbation Theory. The well-known terms occurring in the chiral effective Lagrangian are then supplemented by contributions proportional to powers of the lattice spacing a. The additional terms are constrained by the symmetries of the lattice action and therefore depend on the specific choice of the discretization. The resulting effective theory can be used to analyze the a-dependence of the various quantities of interest – provided the quark masses and the momenta considered are in the range where the truncated chiral perturbation series yields an adequate approximation. Understanding the dependence on the lattice spacing is of central importance for a controlled extrapolation to the continuum limit. For staggered fermions, this program has first been carried out for a single staggered flavor (a single staggered field) <> at $O(a^2)$. In the following, this effective theory is denoted by SχPT. It was later generalized to an arbitrary number of flavours <> <>, and to next-to-leading order <>. The corresponding theory is commonly called Rooted Staggered chiral perturbation theory and is denoted by RSχPT. For Wilson fermions, the effective theory has been developed in <> <> <> and is called WχPT, while the theory for Wilson twisted mass fermions <> <> <> is termed tmWχPT. Another important approach is to consider theories in which the valence and sea quark masses are chosen to be different. These theories are called partially quenched. The acronym for the corresponding chiral effective theory is PQχPT <> <> <> <>. Finally, one can also consider theories where the fermion discretizations used for the sea and the valence quarks are different. The effective chiral theories for these “mixed action” theories are referred to as MAχPT <><><>. ==== References ==== <> === Summary of simulated lattice actions === Summary of simulated lattice actions with $N_f = 2+1$ and $N_f = 2+1+1$ quark flavours. ||<:30%#eeeeee tablestyle="width: 70%;">||<:#eeeeee>'''Ref.'''||<:#eeeeee>'''$N_f$'''||<:#eeeeee>'''gauge action'''||<:#eeeeee>'''quark action'''|| || Aubin 08, 09 ||<:> [[http://arXiv.org/abs/0810.4328|1]], [[http://arXiv.org/abs/0905.3947|2]] ||<:> 2+1 ||<:> tadSym ||<:> Asqtad ^A)^ || || Blum 10 ||<:> [[http://arXiv.org/abs/1006.1311|3]] ||<:> 2+1 ||<:> Iwasaki ||<:> DW || || BMW 10A-C ||<:> [[http://arXiv.org/abs/1011.2403|4]], [[http://arXiv.org/abs/1011.2711|5]], [[http://arXiv.org/abs/1011.4189|6]] ||<:> 2+1 ||<:> tlSym ||<:> 2-HEX tlSW || || BMW 10 ||<:> [[http://arXiv.org/abs/1001.4692|7]] ||<:> 2+1 ||<:> tlSym ||<:> 6-stout tlSW || || CP-PACS/JLQCD 07 ||<:> [[http://arXiv.org/abs/0704.1937|8]] ||<:> 2+1 ||<:> Iwasaki ||<:> npSW || || ETM 10, 10E ||<:> [[http://arXiv.org/abs/1004.5284|9]], [[http://arXiv.org/abs/1012.0200|10]] ||<:> 2+1+1 ||<:> Iwasaki ||<:> tmWil || || HPQCD/UKQCD 06 ||<:> [[http://arXiv.org/abs/hep-lat/0603023|11]] ||<:> 2+1 ||<:> tadSym ||<:> Asqtad || || HPQCD/UKQCD 07 ||<:> [[http://arXiv.org/abs/0706.1726|12]] ||<:> 2+1 ||<:> tadSym ||<:> Asqtad ^B)^ || || HPQCD/MILC/UKQCD 04 ||<:> [[http://arXiv.org/abs/hep-lat/0405022|13]] ||<:> 2+1 ||<:> tadSym ||<:> Asqtad || || JLQCD 09 ||<:> [[http://arXiv.org/abs/0911.5555|14]] ||<:> 2+1 ||<:> Iwasaki ||<:> overlap || || JLQCD/TWQCD 08B, 09A ||<:> [[http://arXiv.org/abs/0810.0085|15]], [[http://arXiv.org/abs/0910.5532|16]] ||<:> 2+1 ||<:> Iwasaki ||<:> overlap || || JLQCD/TWQCD 10 ||<:> [[http://arXiv.org/abs/1012.4052|17]] ||<:> 2+1, 3 ||<:> Iwasaki ||<:> overlap || || LHP 04 ||<:> [[http://arXiv.org/abs/hep-lat/0411028|18]] ||<:> 2+1 ||<:> tadSym ||<:> Asqtad ^A)^ || || MILC 04, 07, 09, 09A ||<:> [[http://arXiv.org/abs/0903.3598|19]], [[http://arXiv.org/abs/hep-lat/0407028|20]] , [[http://arXiv.org/abs/0710.1118|21]] ,[[http://arXiv.org/abs/0910.2966|22]] ||<:> 2+1 ||<:> tadSym ||<:> Asqtad || || NPLQCD 06 ||<:> [[http://arXiv.org/abs/hep-lat/0606023|23]] ||<:> 2+1 ||<:> tadSym ||<:> Asqtad ^A)^ || || PACS-CS 08, 08B, 09, 10 ||<:> [[http://arXiv.org/abs/0911.2561|24]], [[http://arXiv.org/abs/0807.1661|25]], [[http://arXiv.org/abs/0811.2630|26]], [[http://arXiv.org/abs/1006.1164|27]] ||<:> 2+1 ||<:> Iwasaki ||<:> npSW || || RBC/UKQCD 07, 08, 08A, 10, 10A-B ||<:> [[http://arXiv.org/abs/0804.0473|28]], [[http://arXiv.org/abs/1004.0886|29]], [[http://arXiv.org/abs/0804.3971|30]], [[http://arXiv.org/abs/hep-ph/0702042|31]], [[http://arXiv.org/abs/1011.0892|32]], [[http://arXiv.org/abs/1012.4178|33]] ||<:> 2+1 ||<:> Iwasaki ||<:> DW || || SWME 10, 11 ||<:> [[http://arXiv.org/abs/1008.5179|34]],[[http://arXiv.org/abs/1101.2685|35]] ||<:> 2+1 ||<:> tadSym ||<:> Asqtad ^C)^ || || TWQCD 08 ||<:> [[http://arXiv.org/abs/0810.3406|36]] ||<:> 2+1 ||<:> Iwasaki ||<:> DW || ^A)^ The calculation uses domain wall fermions in the valence quark sector. ^B)^ The calculation uses HISQ staggered fermions in the valence quark sector. ^C)^ The calculation uses HYP smeared improved staggered fermions in the valence quark sector. <
> <
> <
> Summary of simulated lattice actions with $N_f = 2$ quark flavours. ||<:30%#eeeeee tablestyle="width: 70%;"> ||<:15%#eeeeee>'''Ref.''' ||<:#eeeeee>'''$N_f$''' ||<:#eeeeee>'''gauge action''' ||<:#eeeeee>'''quark action''' || || ALPHA 05 ||<:> [[http://arXiv.org/abs/hep-lat/0507035|37]] ||<:> 2 ||<:> Wilson ||<:> npSW || || Bernardoni 10 ||<:> [[http://arXiv.org/abs/1008.1870|38]] ||<:> 2 ||<:> Wilson ||<:> npSW ^A)^ || || CERN-TOV 06 ||<:> [[http://arXiv.org/abs/hep-lat/0610059|39]] ||<:> 2 ||<:> Wilson ||<:> Wilson/npSW || || CERN 08 ||<:> [[http://arXiv.org/abs/0812.3638|40]] ||<:> 2 ||<:> Wilson ||<:> Wilson/npSW || || CP-PACS 01 ||<:> [[http://arXiv.org/abs/hep-lat/0105015|41]] ||<:> 2 ||<:> Iwasaki ||<:> mfSW || || ETM 07, 07A, 08, 09, 09A-D, 10B, D ||<:> [[http://arXiv.org/abs/0709.4574|42]], [[http://arXiv.org/abs/hep-lat/0701012|43]], [[http://arXiv.org/abs/0812.4042|44]], [[http://arXiv.org/abs/0904.0954|45]], [[http://arXiv.org/abs/0906.4728|46]], [[http://arXiv.org/abs/0911.5061|47]], [[http://arXiv.org/abs/0911.1931|48]], [[http://arXiv.org/abs/0910.4838|49]], [[http://arXiv.org/abs/1012.3573|50]], [[http://arXiv.org/abs/1010.3659|51]] ||<:> 2 ||<:> tlSym ||<:> tmWil || || HHS 08 ||<:> [[http://arXiv.org/abs/0806.4586|52]] ||<:> 2 ||<:> tadSym ||<:> n-HYP tlSW || || JLQCD 08 ||<:> [[http://arXiv.org/abs/0801.4186|53]] ||<:> 2 ||<:> Iwasaki ||<:> overlap || || JLQCD 02, 05 ||<:> [[http://arXiv.org/abs/0806.0894|54]], [[http://arXiv.org/abs/hep-lat/0510068|55]] ||<:> 2 ||<:> Wilson ||<:> npSW || || JLQCD/TWQCD 07, 08A ||<:> [[http://arXiv.org/abs/hep-lat/0212039|56]], [[http://arXiv.org/abs/0711.4965|57]] ||<:> 2 ||<:> Iwasaki ||<:> overlap || || QCDSF 07 ||<:> [[http://arXiv.org/abs/0710.2100|58]] ||<:> 2 ||<:> Wilson ||<:> npSW || || QCDSF/UKQCD 04, 06, 06A, 07 ||<:> [[http://arXiv.org/abs/hep-lat/0601004|59]], [[http://arXiv.org/abs/hep-ph/0409312|60]], [[http://www.physik.uni-regensburg.de/lat07/hevea/schierholz.pdf|61]], [[http://arXiv.org/abs/hep-lat/0608021|62]] ||<:> 2||<:> Wilson ||<:> npSW || || RBC 07 ||<:> [[http://arXiv.org/abs/0708.0484|63]] ||<:> 2 ||<:> DBW2 ||<:> DW || || RBC 04, 06 ||<:> [[http://arXiv.org/abs/hep-ph/0607162|64]], [[http://arXiv.org/abs/hep-lat/0411006|65]] ||<:> 2 ||<:> DBW2 ||<:> DW || || SPQcdR 05 ||<:> [[http://arXiv.org/abs/hep-lat/0510014|66]] ||<:> 2 ||<:> Wilson ||<:> Wilson || || UKQCD 04, 07 ||<:> [[http://arXiv.org/abs/0710.5136|67]], [[http://arXiv.org/abs/hep-lat/0406013|68]] ||<:> 2 ||<:> Wilson ||<:> npSW || ^A)^ The calculation uses overlap fermions in the valence quark sector.