Review of lattice results concerning low energy particle physics

The latest version of the complete review as of December 2017 is accessible here. It contains updated sections as follows:
$V_{ud}$ and $V_{us}$: updated November 2016
 Lowenergy constants: updated July 2017
 Kaon mixing: updated December 2016
$B$meson mixing parameters: updated July/November 2017
The original complete 2015/2016 review is still accessible from EPJC. The separate sections can be downloaded as separate pdffiles following the links in the table of contents below.
The latest figures can be downloaded in eps, pdf and png format, together with a bibfile containing the bibtexentries for the calculations which contribute to the FLAG averages and estimates. The downloads are available via the menu in the sidebar.
The 2013/2014 review is accessible here or from EPJC.
Contents
Introduction
The introduction with the updated summary tables can be downloaded here.
We review lattice results related to pion, kaon, $D$ and $B$meson physics with the aim of making them easily accessible to the particle physics community. More specifically, we report on the determination of the lightquark masses, the form factor $f_+(0)$, arising in the semileptonic $K$→$\pi$ transition at zero momentum transfer, as well as the decay constant ratio $f_K/f_\pi$ and its consequences for the CKM matrix elements $V_{us}$ and $V_{ud}$. Furthermore, we describe the results obtained on the lattice for some of the lowenergy constants of SU(2)$_L$×SU(2)$_R$ and SU(3)$_L$×SU(3)$_R$ Chiral Perturbation Theory. We review the determination of the $B_K$ parameter of neutral kaon mixing as well as the additional four $B$ parameters that arise in theories of physics beyond the Standard Model. The latter quantities are an addition compared to the previous review. For the heavyquark sector, we provide results for $m_c$ and $m_b$ (also new compared to the previous review), as well as those for $D$ and $B$meson decay constants, form factors, and mixing parameters. These are the heavyquark quantities most relevant for the determination of CKM matrix elements and the global CKM unitaritytriangle fit. Finally, we review the status of lattice determinations of the strong coupling constant $\alpha_s$.
Flavour physics provides an important opportunity for exploring the limits of the Standard Model of particle physics and for constraining possible extensions that go beyond it. As the LHC explores a new energy frontier and as experiments continue to extend the precision frontier, the importance of flavour physics will grow, both in terms of searches for signatures of new physics through precision measurements and in terms of attempts to construct the theoretical framework behind direct discoveries of new particles. A major theoretical limitation consists in the precision with which stronginteraction effects can be quantified. Largescale numerical simulations of lattice QCD allow for the computation of these effects from first principles. The scope of the Flavour Lattice Averaging Group (FLAG) is to review the current status of lattice results for a variety of physical quantities in lowenergy physics. Set up in November 2007 it comprises experts in Lattice Field Theory, Chiral Perturbation Theory and Standard Model phenomenology. Our aim is to provide an answer to the frequently posed question “What is currently the best lattice value for a particular quantity?” in a way that is readily accessible to nonlatticeexperts. This is generally not an easy question to answer; different collaborations use different lattice actions (discretizations of QCD) with a variety of lattice spacings and volumes, and with a range of masses for the $u$− and $d$−quarks. Not only are the systematic errors different, but also the methodology used to estimate these uncertainties varies between collaborations. In the present work we summarize the main features of each of the calculations and provide a framework for judging and combining the different results. Sometimes it is a single result that provides the “best” value; more often it is a combination of results from different collaborations. Indeed, the consistency of values obtained using different formulations adds significantly to our confidence in the results.
The first two editions of the FLAG review were published in 2011 ^{[1]} and 2014 ^{[2]}. The second edition reviewed results related to both light ($u$, $d$ and $s$), and heavy ($c$ and $b$) flavours. The quantities related to pion and kaon physics were lightquark masses, the form factor $f_+(0)$ arising in semileptonic $K$→$\pi$ transitions (evaluated at zero momentum transfer), the decay constants $f_K$ and $f_\pi$, and the $B_K$ parameter from neutral kaon mixing. Their implications for the CKM matrix elements $V_{us}$ and $V_{ud}$ were also discussed. Furthermore, results were reported for some of the lowenergy constants of SU(2)$_L$×SU(2)$_R$ and SU(3)$_L$×SU(3)$_R$ Chiral Perturbation Theory. The quantities related to $D$ and $B$meson physics that were reviewed were the $B$ and $D$meson decay constants, form factors, and mixing parameters. These are the heavylight quantities most relevant to the determination of CKM matrix elements and the global CKM unitaritytriangle fit. Last but not least, the current status of lattice results on the QCD coupling $\alpha_s$ was reviewed.
In the present paper we provide updated results for all the abovementioned quantities, but also extend the scope of the review in two ways. First, we now present results for the charm and bottom quark masses, in addition to those of the three lightest quarks. Second, we review results obtained for the kaon mixing matrix elements of new operators that arise in theories of physics beyond the Standard Model. Our main results are collected in Tabs. 1 and 2.
Table 1 

Summary of the main results of this review, grouped in terms of $N_f$, the number of dynamical quark flavours in lattice simulations. Quark masses and the quark condensate are given in the $\overline{\text{MS}}$ scheme at running scale $\mu = 2$ GeV or as indicated; the other quantities listed are specified in the quoted sections. For each result we provide the list of references that entered the FLAG average or estimate in the bibfile for download. We recommend to consult the detailed tables and figures in the relevant section for more significant information and for explanations on the source of the quoted errors. 

Quantity 
Sec. 
$N_f=2+1+1$ 
$N_f=2+1$ 
$N_f=2$ 
Refs. 
$m_s$ [MeV] 
93.9(1.1) 
92.0(2.1) 
101(3) 

$m_{ud}$ [MeV] 
3.70(17) 
3.373(80) 
3.6(2) 

$m_s/m_{ud}$ 
27.30(34) 
27.43(31) 
27.3(9) 

$m_u$ [MeV] 
2.36(24) 
2.16(9)(7) 
2.40(23) 
^{[3]} 

$m_d$ [MeV] 
5.03(26) 
4.68(14)(7) 
4.80(23) 
^{[3]} 

$m_u/m_d$ 
0.470(56) 
0.46(2)(2) 
0.50(4) 
^{[3]} 

$\overline{m}_c$(3 GeV) [GeV] 
0.996(25) 
0.987(6) 
1.03(4) 

$m_c/m_s$ 
11.70(6) 
11.82(16) 
11.74(35) 

$\overline{m}_b(\overline{m}_b)$ [GeV] 
4.190(21) 
4.164(23) 
4.256(81) 

$f_+(0)$ 
0.9706(27) 
0.9677(27) 
0.9560(57)(62) 

$f_{K^\pm}/f_{\pi^\pm}$ 
1.193(3) 
1.192(4) 
1.205(6)(17) 

$f_{\pi^\pm}$ [MeV] 

130.2(8) 


$f_{K^\pm}$ [MeV] 
155.7(3) 
155.7(7) 
157.5(2.4) 

$\Sigma^{1/3}$ [MeV] 
280(8)(15) 
272(5) 
266(10) 

$F_\pi/F$ [MeV] 
1.077(2)(2) 
1.062(7) 
1.073(15) 

$\bar{\ell}_3$ [MeV] 
3.53(5)(26) 
3.07(64) 
3.41(82) 

$\bar{\ell}_4$ [MeV] 
4.73(2)(10) 
4.02(45) 
4.40(28) 

$\bar{\ell}_6$ [MeV] 


15.1(1.2) 

$\hat B_K$ 
0.717(18)(16) 
0.7625(97) 
0.727(22)(12) 

$B_2$ 
0.46(1)(3) 
0.502(14) 
0.47(2)(1) 

$B_3$ 
0.79(2)(4) 
0.766(32) 
0.78(4)(2) 

$B_4$ 
0.78(2)(4) 
0.926(19) 
0.76(2)(2) 

$B_5$ 
0.49(3)(3) 
0.720(38) 
0.58(2)(2) 
Table 2 

Summary of the main results of this review, grouped in terms of $N_f$, the number of dynamical quark flavours in lattice simulations. The quantities listed are specified in the quoted sections. For each result we provide the list of references that entered the FLAG average or estimate in the bibfile for download. We recommend to consult the detailed tables and figures in the relevant section for more significant information and for explanations on the source of the quoted errors. 

Quantity 
Sec. 
$N_f=2+1+1$ 
$N_f=2+1$ 
$N_f=2$ 
Refs. 
$f_D$ [MeV] 
212.15(1.45) 
209.2(3.3) 
208(7) 

$f_{D_s}$ [MeV] 
248.83(1.27) 
249.8(2.3) 
250(7) 

$f_{D_s}/f_D$ 
1.1716(32) 
1.187(12) 
1.20(2) 

$f_+^{D\pi}(0)$ [MeV] 

0.666(29) 


$f_+^{DK}(0)$ [MeV] 

0.747(19) 


$f_B$ [MeV] 
186(4) 
192.0(4.3) 
188(7) 

$f_{B_s}$ [MeV] 
224(5) 
228.4(3.7) 
227(7) 

$f_{B_s}/f_B$ 
1.205(7) 
1.201(16) 
1.206(23) 

$f_{B_d}\sqrt{\hat{B}_{B_d}}$ [MeV] 

225(9) 
216(10) 

$f_{B_s}\sqrt{\hat{B}_{B_s}}$ [MeV] 

274(8) 
262(10) 

$\hat{B}_{B_d}$ 

1.30(10) 
1.30(6) 

$\hat{B}_{B_s}$ 

1.35(6) 
1.32(5) 

$\xi$ 

1.206(17) 
1.225(31) 

$B_{B_s}/B_{B_d}$ 

1.032(38) 
1.007(21) 

$\alpha^{(5)}_{\overline{\text{MS} } }(M_Z)$ 
0.1182(12) 


$\Lambda^{(5)}_{\overline{\text{MS} } }$ [MeV] 
211(14) 

Our plan is to continue providing FLAG updates, in the form of a peer reviewed paper, roughly on a biennial basis. This effort is supplemented by our more frequently updated website here, where figures as well as pdffiles for the individual sections can be downloaded. The papers reviewed in the present edition have appeared before the closing date 30 November 2015. The section on leptonic and semileptonic kaon and pion decay and $Vud$ and $Vus$ has been updated in November 2016, while the section on kaon mixing has been updated in December 2016.
FLAG composition, guidelines and rules
FLAG strives to be representative of the lattice community, both in terms of the geographical location of its members and the lattice collaborations to which they belong. We aspire to provide the particlephysics community with a single source of reliable information on lattice results.
In order to work reliably and efficiently, we have adopted a formal structure and a set of rules by which all FLAG members abide. The collaboration presently consists of an Advisory Board (AB), an Editorial Board (EB), and seven Working Groups (WG). The role of the Advisory Board is that of general supervision and consultation. Its members may interfere at any point in the process of drafting the paper, expressing their opinion and offering advice. They also give their approval of the final version of the preprint before it is rendered public. The Editorial Board coordinates the activities of FLAG, sets priorities and intermediate deadlines, and takes care of the editorial work needed to amalgamate the sections written by the individual working groups into a uniform and coherent review. The working groups concentrate on writing up the review of the physical quantities for which they are responsible, which is subsequently circulated to the whole collaboration for critical evaluation.
The current list of FLAG members and their Working Group assignments is:
Advisory Board (AB) 
S. Aoki, C. Bernard, M. Golterman, H. Leutwyler, C. Sachrajda 

Editorial Board (EB) 
G. Colangelo, A. Jüttner, S. Hashimoto, S. Sharpe, T. Vladikas, U. Wenger 

Working Groups (WG) 

Quark masses 
L. Lellouch, T. Blum, V. Lubicz 

$V_{us}, V_{ud}$ 
S. Simula, P. Boyle,^{[4]} T. Kaneko 

LEC 
S. Dürr, H. Fukaya, U.M. Heller 

$B_K$ 
H. Wittig, P. Dimopoulos, R. Mawhinney 

$f_{B_{(s)}}$, $f_{D_{(s)}}$, $B_B$ 
M. Della Morte, Y. Aoki, D. Lin 

$B_{(s)}$, $D$ semileptonic and radiative decays 
E. Lunghi, D. Becirevic, S. Gottlieb, C. Pena 

$\alpha_s$ 
R. Sommer, R. Horsley, T. Onogi 
As some members of the WG on quark masses were faced with unexpected hindrances, S. Simula has kindly assisted in the completion of the relevant section during the final phases of its composition.
The most important FLAG guidelines and rules are the following:
 the composition of the AB reflects the main geographical areas in which lattice collaborations are active, with members from America, Asia/Oceania and Europe;
 the mandate of regular members is not limited in time, but we expect that a certain turnover will occur naturally;
 whenever a replacement becomes necessary this has to keep, and possibly improve, the balance in FLAG, so that different collaborations, from different geographical areas are represented;
in all working groups the three members must belong to three different lattice collaborations;^{[5]}
 a paper is in general not reviewed (nor colourcoded, as described in the next section) by any of its authors;
 lattice collaborations not represented in FLAG will be consulted on the colour coding of their calculation;
 there are also internal rules regulating our work, such as voting procedures.
Citation policy
We draw attention to this particularly important point. As stated above, our aim is to make lattice QCD results easily accessible to nonlatticeexperts and we are well aware that it is likely that some readers will only consult the present paper and not the original lattice literature. It is very important that this paper be not the only one cited when our results are quoted. We strongly suggest that readers also cite the original sources. In order to facilitate this, in Tabs. 1 and 2, besides summarizing the main results of the present review, we also cite the original references from which they have been obtained. In addition, for each figure we make a bibtexfile available on this webpage which contains the bibtexentries of all the calculations contributing to the FLAG average or estimate. The bibliography at the end of this paper should also make it easy to cite additional papers. Indeed we hope that the bibliography will be one of the most widely used elements of the whole paper.
General issues
Several general issues concerning the present review are thoroughly discussed in Sec. 1.1 of our initial 2010 paper ^{[1]} and we encourage the reader to consult the relevant pages. In the remainder of the present subsection, we focus on a few important points. Though the discussion has been duly updated, it is essentially that of Sec. 1.2 of the 2013 review ^{[2]}.
The present review aims to achieve two distinct goals: first, to provide a description of the work done on the lattice concerning lowenergy particle physics; and, second, to draw conclusions on the basis of that work, summarizing the results obtained for the various quantities of physical interest.
The core of the information about the work done on the lattice is presented in the form of tables, which not only list the various results, but also describe the quality of the data that underlie them. We consider it important that this part of the review represents a generally accepted description of the work done. For this reason, we explicitly specify the quality requirements^{[6]} used and provide sufficient details in appendices so that the reader can verify the information given in the tables.
On the other hand, the conclusions drawn on the basis of the available lattice results are the responsibility of FLAG alone. Preferring to err on the side of caution, in several cases we draw conclusions that are more conservative than those resulting from a plain weighted average of the available lattice results. This cautious approach is usually adopted when the average is dominated by a single lattice result, or when only one lattice result is available for a given quantity. In such cases one does not have the same degree of confidence in results and errors as when there is agreement among several different calculations using different approaches. The reader should keep in mind that the degree of confidence cannot be quantified, and it is not reflected in the quoted errors.
Each discretization has its merits, but also its shortcomings. For most topics covered in this review we have an increasingly broad database, and for most quantities lattice calculations based on totally different discretizations are now available. This is illustrated by the dense population of the tables and figures in most parts of this review. Those calculations that do satisfy our quality criteria indeed lead to consistent results, confirming universality within the accuracy reached. In our opinion, the consistency between independent lattice results, obtained with different discretizations, methods, and simulation parameters, is an important test of lattice QCD, and observing such consistency also provides further evidence that systematic errors are fully under control.
In the sections dealing with heavy quarks and with $\alpha_s$, the situation is not the same. Since the $b$quark mass cannot be resolved with current lattice spacings, all lattice methods for treating $b$ quarks use effective field theory at some level. This introduces additional complications not present in the lightquark sector. An overview of the issues specific to heavyquark quantities is given in the introduction of Sec. 8. For $B$ and $D$ meson leptonic decay constants, there already exists a good number of different independent calculations that use different heavyquark methods, but there are only one or two independent calculations of semileptonic $B$ and $D$ meson form factors and $B$ meson mixing parameters. For $\alpha_s$, most lattice methods involve a range of scales that need to be resolved and controlling the systematic error over a large range of scales is more demanding. The issues specific to determinations of the strong coupling are summarized in Sec. 9.
Number of sea quarks in lattice simulations:
Lattice QCD simulations currently involve two, three or four flavours of dynamical quarks. Most simulations set the masses of the two lightest quarks to be equal, while the strange and charm quarks, if present, are heavier (and tuned to lie close to their respective physical values). Our notation for these simulations indicates which quarks are nondegenerate, e.g. $N_f =2+1$ if $m_u=m_d<m_s$ and $N_f =2+1+1$ if $m_u=m_d <m_s <m_c$. Calculations with $N_f = 2$, i.e. two degenerate dynamical flavours, often include strange valence quarks interacting with gluons, so that bound states with the quantum numbers of the kaons can be studied, albeit neglecting strange seaquark fluctuations. The quenched approximation ($N_f = 0$), in which sea quark contributions are omitted, has uncontrolled systematic errors and is no longer used in modern lattice simulations with relevance to phenomenology. Accordingly, we will review results obtained with $N_f =2$, $N_f =2+1$, and $N_f =2+1+1$, but omit earlier results with $N_f = 0$. The only exception concerns the QCD coupling constant $\alpha_s$. Since this observable does not require valence light quarks, it is theoretically well defined also in the $N_f = 0$ theory, which is simply pure gluodynamics. The $N_f$dependence of $\alpha_s$, or more precisely of the related quantity $r_0\Lambda_{MS}$, is a theoretical issue of considerable interest; here $r_0$ is a quantity with the dimension of length, which sets the physical scale, as discussed in Appendix A.2 in the Glossary. We stress, however, that only results with $N_f$ ≥ 3 are used to determine the physical value of αs at a high scale.
Lattice actions, simulation parameters and scale setting:
The remarkable progress in the precision of lattice calculations is due to improved algorithms, better computing resources and, last but not least, conceptual developments. Examples of the latter are improved actions that reduce lattice artifacts and actions that preserve chiral symmetry to very good approximation. A concise characterization of the various discretizations that underlie the results reported in the present review is given in Appendix A.1 in the Glossary.
Physical quantities are computed in lattice simulations in units of the lattice spacing so that they are dimensionless. For example, the pion decay constant that is obtained from a simulation is $f_\pi a$, where $a$ is the spacing between two neighboring lattice sites. To convert these results to physical units requires knowledge of the lattice spacing a at the fixed values of the bare QCD parameters (quark masses and gauge coupling) used in the simulation. This is achieved by requiring agreement between the lattice calculation and experimental measurement of a known quantity, which thus “sets the scale” of a given simulation. A few details on this procedure are provided in Appendix A.2 in the Glossary.
Renormalization and scheme dependence:
Several of the results covered by this review, such as quark masses, the gauge coupling, and Bparameters, are for quantities defined in a given renormalization scheme and at a specific renormalization scale. The schemes employed (e.g. regularizationindependent MOM schemes) are often chosen because of their specific merits when combined with the lattice regularization. For a brief discussion of their properties, see Appendix A.3 in the Glossary. The conversion of the results, obtained in these socalled intermediate schemes, to more familiar regularization schemes, such as the MSbarscheme, is done with the aid of perturbation theory. It must be stressed that the renormalization scales accessible in simulations are limited, because of the presence of an ultraviolet (UV) cutoff of ∼ $\pi/a$. To safely match to MSbar, a scheme defined in perturbation theory, Renormalization Group (RG) running to higher scales is performed, either perturbatively or nonperturbatively (the latter using finitesize scaling techniques).
Extrapolations:
Because of limited computing resources, lattice simulations are often performed at unphysically heavy pion masses, although results at the physical point have become increasingly common. Further, numerical simulations must be done at nonzero lattice spacing, and in a finite (fourdimensional) volume. In order to obtain physical results, lattice data are obtained at a sequence of pion masses and a sequence of lattice spacings, and then extrapolated to the physical pion mass and to the continuum limit. In principle, an extrapolation to infinite volume is also required. However, for most quantities discussed in this review, finitevolume effects are exponentially small in the linear extent of the lattice in units of the pion mass and, in practice, one often verifies volume independence by comparing results obtained on a few different physical volumes, holding other parameters equal. To control the associated systematic uncertainties, these extrapolations are guided by effective theories. For lightquark actions, the latticespacing dependence is described by Symanzik’s effective theory; for heavy quarks, this can be extended and/or supplemented by other effective theories such as HeavyQuark Effective Theory (HQET). The pionmass dependence can be parameterized with Chiral Perturbation Theory (χPT), which takes into account the NambuGoldstone nature of the lowest excitations that occur in the presence of light quarks. Similarly, one can use HeavyLight Meson Chiral Perturbation Theory (HMχPT) to extrapolate quantities involving mesons composed of one heavy ($b$ or $c$) and one light quark. One can combine Symanzik’s effective theory with χPT to simultaneously extrapolate to the physical pion mass and the continuum; in this case, the form of the effective theory depends on the discretization. See Appendix A.4 in the Glossary for a brief description of the different variants in use and some useful references. Finally, χPT can also be used to estimate the size of finitevolume effects measured in units of the inverse pion mass, thus providing information on the systematic error due to finitevolume effects in addition to that obtained by comparing simulations at different volumes.
Critical slowing down:
The lattice spacings reached in recent simulations go down to 0.05 fm or even smaller. In this regime, long autocorrelation times slow down the sampling of the configurations. Many groups check for autocorrelations in a number of observables, including the topological charge, for which a rapid growth of the autocorrelation time is observed with decreasing lattice spacing. This is often referred to as topological freezing. A solution to the problem consists in using open boundary conditions in time, instead of the more common antiperiodic ones. More recently two other approaches have been proposed, one based on a multiscale thermalization algorithm and another based on defining QCD on a nonorientable manifold. The problem is also touched upon in Sec. 9.2, where it is stressed that attention must be paid to this issue. While large scale simulations with open boundary conditions are already far advanced, unfortunately so far no results reviewed here have been obtained with any of the above methods. It is usually assumed that the continuum limit can be reached by extrapolation from the existing simulations and that potential systematic errors due to the long autocorrelation times have been adequately controlled.
Simulation algorithms and numerical errors:
Most of the modern latticeQCD simulations use exact algorithms, which do not produce any systematic errors when exact arithmetic is available. In reality, one uses numerical calculations at double (or in some cases even single) precision, and some errors are unavoidable. More importantly, the inversion of the Dirac operator is carried out iteratively and it is truncated once some accuracy is reached, which is another source of potential systematic error. In most cases, these errors have been confirmed to be much less than the statistical errors. In the following we assume that this source of error is negligible. Some of the most recent simulations use an inexact algorithm in order to speedup the computation, though it may produce systematic effects. Currently available tests indicate that errors from the use of inexact algorithms are under control.
References
G. Colangelo et al., Review of lattice results concerning low energy particle physics, Eur. Phys. J. C71 (2011) 1695, arXiv:1011.4408 (1 2)
S. Aoki et al., Review of lattice results concerning lowenergy particle physics, Eur. Phys. J. C74 (2014) 2890, arXiv:1310.8555. (1 2)
This is a FLAG estimate, based on $\chi$PT and the isospin averaged up and downquark mass $m_{ud}$. (1 2 3)
Peter Boyle had participated actively in the early stages of the current FLAG effort.
Unfortunately, due to other commitments, it was impossible for him to contribute until the end, and he decided to withdraw from the collaboration. (1)The WG on semileptonic D and B decays has currently four members, but only three of them belong to lattice collaborations. (1)
We also use terms like “quality criteria”, “rating”, “colour coding” etc. when referring to the classification of results, as described in Sec. 2. (1)
Quality criteria
The section on the quality criteria can be downloaded here.
Quark masses
The section on the quark masses can be downloaded here.
$\small{V_{ud}}$ and $\small{V_{us}}$
The section on $V_{ud}$ and $V_{us}$ updated in November 2016 can be downloaded here. The original section is still available here.
Lowenergy constants
The section on the Lowenergy constants updated in Juy 2017 can be downloaded here. The original section is still available here.
Kaon mixing
The section on the Kaon mixing updated in December 2016 can be downloaded here. The original section is still available here.
$\small{D}$meson decay constants and form factors
The section on the $D$meson decay constants and form factors can be downloaded here.
$\small{B}$meson decay constants, mixing parameters, and form factors
The section on the $B$meson decay constants, mixing parameters, and form factors updated in July and November 2017 can be downloaded here. The November 2017 update contains additional explanations on the error treatment for B meson bag parameters.
The strong coupling $\alpha_s$
The section on the strong coupling $\alpha_s$ can be downloaded here.
Glossary
The glossary can be downloaded here.
Notes
Notes to the various sections can be downloaded here.
This page has been accessed 86773 times.