# Chiral behaviour of the pion decay constant

in Qcd

###### Abstract

As increased statistics and new ensembles with light pions have
become available within the CLS effort, we complete previous
work by inspecting the chiral behaviour of the pion decay constant.
We discuss the validity of Chiral Perturbation Theory (PT) and
examine the results concerning the pion decay constant and the ensuing
scale setting, the pion mass squared in units of the quark mass,
and the ratio of decay constants ; along
the way, the relevant low-energy constants of SU(2) PT are estimated.
All simulations were performed with two dynamical flavours of
nonperturbatively O(a)-improved Wilson fermions, on volumes with
, pion masses 192 MeV and lattice spacings
down to 0.048 fm. Our error analysis takes into account
the effect of slow modes on the autocorrelations.

DESY 13-213

Chiral behaviour of the pion decay constant

in QCD

for the ALPHA Collaboration

\abstract@cs

## 1 Introduction

The nonperturbative features of quantum chromodynamics (QCD) are more and more prominent as one approaches the low-energy regime; there, the most successful strategy to determine its properties is, to date, numerical simulation on the lattice, complemented by insight provided by chiral perturbation theory (PT) as to how the chiral limit is reached.

In the past years, new ensembles have been produced within the Coordinated Lattice Simulations (CLS) effort, closer and closer to both the physical point and the continuum limit. Here we report on the chiral behaviour of the pion decay constant (and related quantities) obtained from the latest set of CLS ensembles, and the subsequent determinations of the PT low-energy constants involved. The analysis presented here can be seen as a complement to Ref. [1] (to which we refer for most of the setup details), where analogous investigations were carried on, among other topics, concerning the kaon decay constant .

All ensembles were generated using -improved Wilson lattice action with two degenerate dynamical flavours, implementing either Domain-Decomposition [2] or Mass-Preconditioning [3]. Table 1 summarises the relevant information on the ensembles employed in the present analysis.

Two-point functions were measured using 10 to 20 stochastic sources per configuration, and subsequently used to extract meson- and PCAC-masses and decay constants, as detailed in [1], with a statistical uncertainty on the level of 1% or better. For renormalisation, the -factors come from one-loop perturbation theory [4], and the -factors are nonperturbatively determined [1].

## 2 Chiral analyses

In the data analysis, first an independent variable is built, parameterising the approach to the chiral and the physical points; in contrast to the one employed in [1], here we deal with the purely pionic (i.e. light-light flavours) variable

(2.0) |

appearing as expansion parameter in the PT formulae. and are the measured quantities at the finite quark masses; the physical point corresponds to .

Generally, fits are performed, simultaneously for all values of , to the functional form coming from PT, taking correlations among measurement into account. Our goal being the chiral and physical points, we inspect fits in the range with MeV; as it turns out, beyond NLO, a quadratic term has to be included for the fit to be stable in ; systematic uncertainties in the resulting parameters were estimated also by altering the fit functions.

( fm) | ( fm) | ( fm) | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|

MDU | MDU | MDU | |||||||||

A2 | 629 | 8000 | 120 | E4 | 580 | 2496 | 10 | N4 | 551 | 3752 | 4 |

A3 | 492 | 8032 | 120 | E5f | 436 | 16000 | 60 | N5 | 440 | 3808 | 4 |

A4 | 383 | 8096 | 120 | E5g | 436 | 16000 | 120 | N6 | 340 | 8040 | 40 |

A5 | 330 | 4004 | 160 | F6 | 311 | 4800 | 36 | O7 | 267 | 3920 | 20 |

B6 | 281 | 1272 | 50 | F7 | 266 | 9416 | 70 | ||||

G8 | 192 | 1114 | 20 |

Correlations are propagated down to the final quantities. Special care is taken for the integrated autocorrelation time, by applying the technique developed in [5]: the effects of slow modes in the transition matrix, to which an observable can couple, are estimated by attaching an tail, with coming from prior determinations, to the autocorrelation function where its statistical error is too large for direct determination. For all quantities, however, the largest contribution (about 30 to 60%) to the full uncertainty comes from the error on the renormalisation factor .

### 2.1 Analysis of

Here and in the following, we denote decay constants in lattice units with an uppercase symbol, as in , with taking the value MeV at the physical point. The outcome of this study, focused on how and to which extent can PT be trusted to describe the behaviour of , will then provide a scale-setting prescription.

The main fits are performed to the PT NLO functional [6], plus a quadratic term:

(2.0) |

where the three (one for each ) will yield the scale setting and encodes the NLO low-energy constant (LEC) :

(2.0) |

This functional form fits the data in a stable way for various (Fig. 1, left), and we take the results from the cut at 500 MeV as median value; systematic uncertainties come from comparing functional forms without and/or the chiral logarithms. The absence of -effects is illustrated by the collapse-plot of Fig. 1, right top, where data for each are rescaled to a common curve. An important remark is the following: the result we get for is such that, in the data range, the -term contributes by as much as 15% to the curve. This is reflected in Fig. 1, right bottom, that illustrates the behaviour of as a function of the pion-mass cut.

This procedure leads to the following determination of the lattice spacing , compatible with the one of [1] albeit slightly smaller in central value (the first error is statistical, the second systematic):

(2.0) |

This is also compatible with an update of [1], based on and including all new ensembles.

Another relevant quantity is the ratio between at the physical- and chiral-points, which – to NLO – encodes directly the LEC :

(2.0) |

the former falls within the world average by the 2013 FLAG review [7].

On a related note, we mention that a continuum-extrapolation of , done in a similar fashion, yields for the hadron parameter a slightly smaller (but still compatible) value than presented in [1]: fm.

### 2.2 Analysis of

We now turn to the ratio of the squared pion mass to the (renormalised) quark mass , expressed, as the chiral condensate below, in the scheme at a scale GeV. The PT NLO expression extends the Gell-Mann-Oakes-Renner (GMOR) relation [8] to:

(2.0) |

with the linear coefficient encoding the NLO low-energy constant .

One is primarily interested in extracting an estimate for the chiral condensate in the continuum limit: with this goal in mind, we encode the expected -scaling directly in the overall amplitude of the above formula, and look for the (third root of the) condensate in units of the physical pion decay constant. We also allow for the usual NNLO analytic term in the fit function:

(2.0) |

(the validity of the Ansatz of -scaling is corroborated by variants of the analysis in which first a -dependent amplitude is obtained, then the continuum limit is taken; the values of used here are those coming from the -based scale setting).

Eq. 2.2 describes well the data, which lie approximately flat in (Fig. 2, left), and the value of is stable against the usual variants of the fit procedure; , on the other hand, is more pion-mass-cut- and fit-function-dependent, leading to a large systematic error on . From the fit to the above equation, with MeV, we get the following values:

(2.0) |

The former result compares favorably to most recent two-flavours lattice determinations [9, 10]. We remark that a more first-principle lattice determination of , minimally relying on PT assumptions, is currently being carried on [11] based on the method proposed in [12] and using the same CLS configurations employed here; the same approach, in the context of twisted-mass fermions, is pursued in Ref. [10], with an outcome which overshoots slightly what found here – by about 1.5 standard deviations) – once recast as (there, this combination is used to overcome current discrepancies in the physical determination of ).

### 2.3 The ratio

We now consider a heavier valence quark and turn to study “kaon” properties, by moving along the trajectory (dubbed ‘strategy 1’ in [1]) defined by . Partially-quenched PT dictates the behaviour of the quantity [13], which in our range is almost linear in . Nevertheless, as usual, we allow for a term in the fit function:

(2.0) |

Thanks to the stabilising effect of the ensemble with the lightest pion (G8), the ratio at the physical point can be accurately determined, see Fig. 2, right; moreover, the presence of the quadratic term has an effect similar to that of Fig. 1, bottom right. Adding -terms as a check did not alter the picture at all. We then quote the result from the fit to the above equation with MeV, which agrees with most two-flavours lattice estimates, and the CKM matrix element obtained from the latter through the values in [14, 15]:

(2.0) |

## 3 Conclusions

This work presents an analysis of the chiral behaviour of pion-related quantities based on two-flavour lattice QCD. The agreement of the lattice spacing determination with the previous ones (from ) encourages us; compared to the kaon-based analysis, however, here -terms are necessary to stabilise the fits. The ubiquitous need for such terms makes one wonder to which extent is PT-behaviour observed: indeed, their rôle in the fits is to approximately cancel the curvature of the NLO chiral logarithms and restore a roughly linear dependence in a wide range of . A possible reading of this finding is that the chiral logarithms set in only beyond the pion masses well probed by our data. Similar observations, namely that the NNLO terms effectively cancel out the NLO logarithms, thus restoring a linear behaviour (e.g. for ), and that linearity sets in at relatively low pion masses already, have been reiterated in the detailed analysis of [16]. The continuum and chiral limits are carefully taken thanks to our ensembles covering a comfortable region in the plane (in particular, the lightest-pion ensemble G8 seems pivotal for stability).

To complete this preliminary investigation, a similar work on the -expansion for PT formulae (i.e. in terms of quark mass instead of pion mass squared), as well as full inclusion of NNLO chiral logarithms, are planned, with the goal of a better assessment of systematic uncertainties; moreover, comparative analyses of the same data, based on the Sommer scale [17] and on the recently introduced flow-time scale [18], are currently being undertaken [19]. A more complete and articulated account of the present analysis is deferred to a forthcoming paper.

The authors gratefully acknowledge access to HPC resources in the form of
a regular GCS/NIC project, a JUROPA/NIC project^{1}^{1}1See
http://www.fz-juelich.de/ias/jsc/EN/Expertise/Supercomputers/ComputingTime/Acknowledgements.html.
and through PRACE-2IP, receiving funding from the
European Community’s Seventh Framework Programme (FP7/2007-2013) under grant
agreement RI-283493.
This work is supported in part by the grants SFB/TR9
of the Deutsche Forschungsgemeinschaft.

## References

- [1] P. Fritzsch, F. Knechtli, B. Leder, M. Marinkovic, S. Schaefer, et al., The strange quark mass and Lambda parameter of two flavor QCD, Nucl.Phys. B865 (2012) 397–429, [arXiv:1205.5380].
- [2] M. Lüscher, Schwarz-preconditioned HMC algorithm for two-flavour lattice QCD, Comput.Phys.Commun. 165 (2005) 199–220, [hep-lat/0409106].
- [3] M. Marinkovic and S. Schaefer, Comparison of the mass preconditioned HMC and the DD-HMC algorithm for two-flavour QCD, PoS LATTICE2010 (2010) 031, [arXiv:1011.0911].
- [4] S. Sint and P. Weisz, Further one loop results in O(a) improved lattice QCD, Nucl.Phys.Proc.Suppl. 63 (1998) 856–858, [hep-lat/9709096].
- [5] ALPHA Collaboration, S. Schaefer, R. Sommer, and F. Virotta, Critical slowing down and error analysis in lattice QCD simulations, Nucl.Phys. B845 (2011) 93–119, [arXiv:1009.5228].
- [6] J. Gasser and H. Leutwyler, Chiral Perturbation Theory to One Loop, Annals Phys. 158 (1984) 142.
- [7] FLAG Working Group, “Review of lattice results concerning low energy particle physics.” http://itpwiki.unibe.ch/flag. (Made available as wiki prior to submission).
- [8] M. Gell-Mann, R. Oakes, and B. Renner, Behavior of current divergences under SU(3) x SU(3), Phys.Rev. 175 (1968) 2195–2199.
- [9] B. B. Brandt, A. Jüttner, and H. Wittig, The pion vector form factor from lattice QCD and NNLO chiral perturbation theory, arXiv:1306.2916.
- [10] K. Cichy, E. Garcia-Ramos, and K. Jansen, Chiral condensate from the twisted mass Dirac operator spectrum, JHEP 1310 (2013) 175, [arXiv:1303.1954].
- [11] G. P. Engel, L. Giusti, S. Lottini, and R. Sommer, Chiral condensate from the Banks-Casher relation, PoS LATTICE2013 (2013) 119, [arXiv:1309.4537].
- [12] L. Giusti and M. Luscher, Chiral symmetry breaking and the Banks-Casher relation in lattice QCD with Wilson quarks, JHEP 0903 (2009) 013, [arXiv:0812.3638].
- [13] S. R. Sharpe, Enhanced chiral logarithms in partially quenched QCD, Phys.Rev. D56 (1997) 7052–7058, [hep-lat/9707018].
- [14] M. Antonelli, V. Cirigliano, G. Isidori, F. Mescia, M. Moulson, et al., An Evaluation of Vus and precise tests of the Standard Model from world data on leptonic and semileptonic kaon decays, Eur.Phys.J. C69 (2010) 399–424, [arXiv:1005.2323].
- [15] J. Hardy and I. Towner, Superallowed 0+ —> 0+ nuclear beta decays: A New survey with precision tests of the conserved vector current hypothesis and the standard model, Phys.Rev. C79 (2009) 055502, [arXiv:0812.1202].
- [16] S. Dürr, Z. Fodor, C. Hoelbling, S. Krieg, T. Kurth, et al., Lattice QCD at the physical point meets SU(2) chiral perturbation theory, arXiv:1310.3626.
- [17] R. Sommer, A New way to set the energy scale in lattice gauge theories and its applications to the static force and alpha-s in SU(2) Yang-Mills theory, Nucl.Phys. B411 (1994) 839–854, [hep-lat/9310022].
- [18] M. Lüscher, Properties and uses of the Wilson flow in lattice QCD, JHEP 1008 (2010) 071, [arXiv:1006.4518].
- [19] M. Bruno and R. Sommer, On the -dependence of gluonic observables, PoS LATTICE2013 (2013) 321, [arXiv:1311.5585].