Application of Bayesian statistics to the sector of decay constants in three-ﬂavour χ PT

. Abstract The sector of decay constants of the octet of light pseudoscalar mesons in the framework of ’resummed’ SU (3) chiral perturbation theory is investigated. A theoretical prediction for the decay constant of η -meson is compared to a range of available determinations. Compatibility of these determinations with the latest ﬁts of the SU (3) low energy coupling constants is discussed. Using a Bayesian statistical approach, constraints on the low energy coupling constants L r 4 and L r 5 , as well as higher order remainders to the decay constants F K and F η , are extracted from the most recent experimental and lattice QCD inputs for the values of the decay constants.


Introduction
Decay constants of the octet of light pseudoscalar mesons have a deep connection to the spontaneous breaking of chiral symmetry.Within the standard framework of SU (3) L × SU (3) R chiral perturbation theory (χPT) [1], the effective theory of quantum chromodynamics at low energies, the decay constants are directly connected to the renormalization and diagonalization of the kinetic part of the Lagrangian.
Starting from the effective generating functional e iZ ef f [v,a,s,p ] = DU e i d 4 x L ef f [U,v,a,s,p ] , (1) where ϕ a (x) are the pseudo-Goldstone boson fields collected in a matrix field U (x), one can obtain the connected n-point Green functions G a1...an P1...Pn (p 1 , . . .p n ) as on-shell residues of the Fourier transformed functional derivatives of Z ef f [v, a, s, p ] with respect to the axial vector sources a i .P 1 . . .P n are the pseudo-Goldstone bosons in the inand out-states with momenta p 1 , . . .p n .One then finds the following relation between the Green functions and the the elements of the scattering matrix A P1...Pn (p 1 , . . .p n ) [2] G a1...an P1...Pn (p 1 , . . .p n ) = F a1 P1 . . .F an Pn A P1...Pn (p 1 , . . .p n ).(3) a kolesar@ipnp.mff.cuni.czb jara.riha@email.czF ai Pi are the generalized decay constants which might in general include mixing terms.They correspond to the renormalization of the external legs of Feynman diagrams.Due to the mixing, it is necessary to distinguish between the on-shell particles in the in-and out-states P 1 . . .P n and the fields ϕ ai , coupled to the axial-vector currents.
The leading order (LO) effective Lagrangian has the form Here in the case when the scalar external sources are taken to be the quark masses.There is no mixing in the kinetic part of the effective Lagrangian at the leading order and thus it's straightforward to see that all the decay constants are equal to the low energy constant F 0 -the fundamental order parameter of the broken chiral symmetry.At the next-to-leading order (NLO), taking here the part containing the low-energy coupling constants (LECs) L 1 . . .L 10 π 0 -η mixing occurs in the kinetic part of the effective action as an isospin breaking effect, inversely proportional to the difference of light quark masses m d − m u .If the π 0 -η mixing is neglected, the only two terms contributing to the renormalization of the kinetic part are the ones proportional to the LECs L 4 and L 5 .
As can be seen, the sector of decay constants in the isospin limit can be viewed as the simplest self-contained subsystem of the theory, only involving two LECs at the leading order (F 0 , B 0 ) and two at the next-to-leading order (L 4 , L 5 ).Quite intriguingly, none of these four constants is known with very high certainty.As will be shown in more detail in what follows, at leading order a significant suppression of the order parameters, compared to the two-flavour values, is still possible or even probable, given the recent results from phenomenology [3] and lattice QCD [4].At next-to-leading order, the constant L r 4 is expected to be small due to its suppression in the limit of large number of colours [5], but if it is indeed the case is still unknown [3].Depending on the above, the value of L r 5 can also vary widely [3,6].The motivation of our work is to investigate this segment of χPT by Bayesian statistical methods, the question being how much can be told about the low-energy coupling constants just by restricting ourselves to this sector.We do not neglect the higher orders, but treat them as a source of statistical uncertainty, thus avoiding the large number of LECs appearing at next-to-next-leading order (NNLO) [7].The framework of 'resummed' χPT [8] is very well suited for such an approach.
Naturally, such a task might have been accomplished long ago, if the values of all the decay constants were known with sufficient precision.While that has indeed been the case for the pions and the kaons, the value of the decay constant of the η meson was considered to be very uncertain due to its strong mixing with η ′ .In fact, in χPT calculations the η decay constant has been usually treated using its chiral expansion, not as an independent observable [2,9].Our crucial input is thus a recent calculation of the η-η ′ sector on lattice QCD by the RQCD collaboration [10], which allows us to derive the SU (3) decay constant F η with some confidence.This work is a continuation of our initial inquiries [11] and [12], with the major new ingredients being the updated input for F η from lattice QCD and Bayesian statistical analysis, implemented in a numerical way.It can also be noted that the η decay constant has not been used as input for the purpose of extraction of the low-energy parameters of the SU (3) χPT until now.
The paper is organized in the following way -Section 2 provides a concise summary of our theoretical framework, while Section 3 introduces the sector of decay constants and connected phenomenology in a more detailed way.Our implementation of the Bayesian statistical analysis is outlined in Section 4, while the employed assumptions are discussed in Section 5. Section 6 then presents the results of the paper, which are subsequently summarized in Section 7.

Resummed χPT
We use an approach to chiral perturbation theory, dubbed 'resummed' χPT [2,8,13], which was proposed as a way to accommodate the possibility of an irregular convergence of the chiral expansion.Such a scenario might occur if some of the leading order LECs (F 0 or B 0 ) were suppressed to a sufficient degree, so that the leading order was not dominant in the chiral expansion.In such a case the chiral series should be handled carefully, as unexpectedly large higher orders might result from reordering of the expansion.In our case, we will assume a large range of possible values of the leading order constants, so various scenarios are naturally possible.
Let us summarized the procedure in a few points: • We use the standard χPT Lagrangian (4, 6), based on the usual power counting m q ∼ O(p 2 ) [14].• Expansions of quantities related linearly to Green functions of QCD currents are trusted ("safe observables").We assumed that for these expansions the NNLO and higher order terms are reasonably small, though not necessary negligible.Leading order terms are not required to be dominant.• The expansions are expressed explicitly to next-toleading order, all higher order contribution are summed into higher order remainders.Thus for an observable A the 'resummed' chiral expansion has the form • These higher order remainders will not be neglected, but estimated and treated as sources of error.In general, they might have a non-trivial analytical structure, though this is not the case for the decay constants.All higher order LECs are effectively contained in the remainders, the large number of NNLO constants is thus traded off for a relatively smaller number of remainders.

Decay constants
Decay constants of the light pseudoscalar meson nonet, consisting of the pions, kaons, η a η ′ , can be introduced in terms of the QCD axial-vector currents where A a µ = qγ µ γ 5 λ a q.The pion and kaon decay constants take a straightforward form in the isospin limit and their values are very well established from either experimental data or lattice QCD calculations [4,15].In contrast, η nad η ′ decay constants were not very well known, until quite recently, due to significant mixing.A lot of theoretical and phenomenological work has thus been devoted to the η-η ′ sector, see e.g.[16][17][18][19][20][21][22][23][24].This list, far from exhaustive, includes investigations of the sector in the U (3) L × U (3) R large N c framework as well as phenomenological studies, aiming to extract the values of the decay constants and related mixing angles from experimental inputs.The results of the phenomenological studies span quite a range of values, some of which are not compatible with others (see [10] for a detailed overview).
Masses of the η and η ′ mesons in a scheme with a single mixing angle were obtained in lattice QCD simulations around a decade ago [25][26][27].However, until recently, to our knowledge only the EMT collaboration [28] has attempted to calculate the full sector of mixing parameters, which they did in the quark flavour basis.Finally, as already mentioned, a comprehensive study of the sector, which goes down to the physical pion mass, has now been published by the RQCD collaboration [10].
The SU (3) decay constant F η , which is the point of our interest, is defined identically to F 8 η in (8) and can therefore be related to the mixing parameters in the U (3) octet-singlet basis For the purpose of this work, our main input will be the recent lattice QCD determination of F 8 η by the RQCD collaboration [10] For comparison, we will also use two model dependent results from phenomenology Here, EGMS15 [22] is a more recent determination which is representative of lower values of this observable, better compatible with (10).On the other hand, EF05 [19] lies on the opposite end of the spectrum and is an example of a very high value of F 8 η .Reported uncertainties are quite low in both cases and thus these results are essentially incompatible with each other.
In the framework of SU (3) L × SU (3) R chiral perturbation theory [1], the chiral expansion of the pseudoscalar meson octet in the isospin limit can be written in the following way [11] This form is obtained directly from the generation functional of two-point Green functions in the logic of 'resummed' approach to χPT [8].A strict form of the chiral expansion is used, where the original parameters of the Lagrangian are retained, thus avoiding any reordering of the series.F 2 P δ F P are the sum of all higher orders, the higher order remainders, which are not neglected.These effectively contain all low-energy coupling constants at higher orders.It should be also noted that in this case the remainders are real constants with no analytical structure and no scale dependence.
Chiral logarithms are denoted as where m P are the pseudoscalar masses at leading order.In particular As can be seen from (13)(14)(15), chiral expansions of the decay constants up to next-to-leading order does indeed depend only on the two leading-order and two next-toleading order LECs -F 0 , B 0 and L r 4 , L r 5 , respectively.The sector of decay constants can thus be considered as a simple, self-contained system, which can be investigated on its own.
For convenience, we introduce a reparametrization of the chiral order parameters F 0 and B 0 where Σ 0 is the three-flavour chiral condensate and M π is the physical pion mass.Such a reparametrization is convenient as the parameters X and Z are restricted to the range (0, 1).Furthermore, the so-called paramagnetic inequality [29] puts an upper bound in the form of the twoflavour LO LECs: where the two-flavour parameters are defined analogously to the three-flavour ones in (19).
Standard approach to the chiral perturbation series usually assumes values of X and Z reasonably close to one, with the leading order dominating the expansion.On the other hand, Z = 0 would correspond to a restoration of chiral symmetry, while X = 0 to a scenario with a vanishing chiral condensate, which also implies Y = 0.
The most recent NNLO standard χPT fit [3] provides two different sets for the NLO LECs.It's based on a large number of inputs, including πK and ππ scattering lengths, K l4 form factors and pion scalar and vector form factors.It also uses the ratio F K /F π (but not F η ).Overall, it uses 16 input observables to fit 8+34 NLO and NNLO parameters.The main fit (BE14) fixes L r 4 by hand, in order to ensure the expected suppression in the large N c limit [5].FF14 (free fit) releases this constraint.Their results for L r 4 and L r 5 are (at µ = 770 MeV): and As can be seen, the obtained values are quite different.The difference is less pronounced for the LO LECs: and For comparison, quite different values were obtained in [30] by constructing SU (3) amplitudes in a Large N c framework.F 0 was found to be very large (88.1±4.1 MeV) and L r 4 compatible with zero ((−0.05 ± 0.22) × 10 −3 ).The Flavour Lattice Averaging Group [4,6] cites several lattice QCD determinations of L r 4 and L r 5 .The last report [4] highlights the results by HPQCD [31]: and MILC [32] The leading-order LECs have also been recently calculated on lattice by the χQCD collaboration [33].Though the results have not been fully published yet, the work has been cited by the Flavour Lattice Averaging Group [4] with a favorable rating.While there are several other older determinations, for example by the MILC collaboration [32,34,35] or based on RBC/UKQCD [36], χQCD provides the first highly-rated calculation of these parameters in more than a decade, as far as we are aware of.The results were quoted by FLAG in the following form: We will use these values as alternative inputs for the leading-order LECs.The purpose of this work is twofold -first, we will show that the 'resummed' χPT framework leads to a simple, but robust prediction for F η .Then we will use the values of F 8 η (10-12) as an input and use Bayesian statistical inference to obtain constraints on the higher order remainders δ F K , δ Fη and the NLO LECs L r 4 and L r 5 .We will compare these results with the two versions of the fit [3] (BE14 and FF14) and lattice QCD values (HPQCD 13A [31] and MILC 10 [32]) and thus check the compatibility of the various values of F η and NLO LECs.

Bayesian statistical analysis
We use a statistical approach based on the Bayes' theorem [8,37] where P (X i |data) is the probability density function (PDF) of an explored set of theoretical parameters X i having a specific value given some experimental data.
In the case of independent experimental inputs, P (data|X i ) is the known probability density of obtaining the observed values of the observables O k in a set of experiments with uncertainties σ k under the assumption that the true values of X i are known, typically given as a normal distribution (28) are prior probability distributions of X i .We use them to implement theoretical assumptions, available experimental information and uncertainties connected with our parameters.
Traditionally, the prior has been understood as a degree of subjective belief.However, in our view, one does not necessarily needs to 'believe' in the validity of the prior in a scientific context, which we think can then be more appropriately interpreted as the quantification of available information and beyond that, the assumptions entering the analysis.Naturally, predictions might depend on the assumptions used.The Bayesian formalism allows us to straightforwardly implement a variety of assumptions and explore their consequences, which we consider to be an important feature of this approach.
In our case we have three observables in the form of the three decay constants F π , F K , F η .We will consider the ratios of quark masses as known and thus we are left with the following free theoretical parameters: -leading order: Z, Y -next-to-leading order: L r 4 , L r 5 -higher orders: δ Fπ , δ F K , δ Fη .
As discussed in the next section, we will use several assumptions about these parameters, which will determine the prior distributions.
Our implementation of the Bayesian statistical analysis is numerical.It consists of two steps -first we numerically generate a large ensemble of theoretical predictions O th k (X i ) for the decay constants, depending on the free parameters, and then we calculate the probability density functions (28), effectively using Monte Carlo integration.

Assumptions
For the LO LECs F 0 and B 0 , we use similar theoretical constraints as in [37], which define our priors for these parameters.Their approximate range then is where the explicit form of Y max , derived in [8], is Here δ Mπ and δ M K are higher order remainders for the chiral expansions of the pseudoscalar masses, which we treat analogously to the remainders of the decay constants (see (40) below).
In order to calculate the prior distributions, we consider X and Z as the primary variables: P (Y, Z|data) dY dZ ≡ ≡ P (X, Z|data)| X→ZY dY ZdZ = = P (X, Z|data) dXdZ.(34) One might naturally ask, why not use Y and Z as the primary variables, as these are the free parameters in our case.That would mean using uniform distributions in the range 0 < Y < Y max and 0 < Z < 1 as the a priori assumption.However, this leads to a quickly rising probability distribution for X = ZY towards zero.Quite clearly, a very small chiral condensate is not a reasonable initial expectation.On the other hand, starting with uniform distributions for 0 < X < 1, 0 < Z < 1 and adding Y < Y max ensures a relatively flat prior for X and a vanishing distribution for Z at Z = 0.That we find reasonable, as it excludes the scenario with unbroken chiral symmetry and thus a world without the pseudo-Goldstone bosons.Then by including the paramagnetic inequality (20) we obtain the set (30)(31)(32).These assumptions lead to probability distributions for the priors depicted in Figure 1.
For the purpose of obtaining constraints on the NLO LECs L r 4 and L r 5 , we will use the determination of Y from η → 3π decays [37] as an additional assumption This value was obtained from a Bayesian analysis of the η → 3π decay widths of the two decay channels and the Dalitz parameter a in the charged channel.The tendency towards Y > 1 is ultimately tied to the very large overall experimental decay rate compared to the simple estimate at leading order, given that the isospin violating parameter R is now known with a fairly good precision from lattice QCD [4].However, while this value is higher than the result of the fits BE14 and FF14 (23)(24).it is compatible and does provide us with a reasonable uncertainty range.As can be seen in Figure 2, this input effectively excludes very low values of Y , which correspond to a significantly suppressed chiral condensate.In other words, such a scenario can be understood to be excluded by the phenomenology of the η → 3π decays.
As an alternative, we will also use the recent result by the χQCD Collaboration (27), expressed in the form: In comparison with our main input (35), χQCD21 also excludes high values of Y and fixes Z in a narrow range at a fairly low value.The prior distributions can be seen in Figure 3.
The three sets of priors introduced above should be understood as proceeding from more conservative to more restricted.The set (30)(31)(32)(33)(34), is based on very general consideration rooted in QCD.The additional inputs (35) or alternatively (36) then supplement an assumption about the value of the leading order LECs.Here, (35) is a more conservative one, effectively only excluding very low values of Y, based on results from η → 3π decays.It is compatible with all values quoted in Section 3. On the other hand, (36) uses very specific values from a recent lattice QCD study, which might be in tension with some other determinations.Our main results will be based on the more conservative assumption of the first two sets of priors, while the third one will be used to explore the consequence of assuming a particular value of the chiral order parameters, which are not yet firmly established, though.
Furthermore, while the Bayes theorem (28) separates the probability distributions P (data|X i ), which depend on the theoretical predictions O th k (X i ), and the priors P (X i ), it's quite desirable to implement the priors on the level of theoretical predictions, thus effectively incorporate the priors P (X i ) into P (data|X i ), so one can examine the theoretical predictions including a realistic set of assumptions.Thus we implement the relations (30)(31)(32)(33)(34), i.e. our default prior, when numerically generating the theoretical predictions O th k (X i ).As for the NLO LECs L r 4 and L r 5 , where our goal is to extract constraints, we limit them to the range (at µ = 770 MeV) which we implement as a uniform distribution.The choice of the uniform distribution signifies the lack of preference for any particular value in the allowed range, which we think is appropriate given the range of values for these parameters available in the literature.Hence we do not use any particular value as our prior and the results are therefore not directly dependent on previous analyses.The allowed range was chosen in a way to cover the values of all recent determinations we are aware of, including NNLO χPT [3] and lattice QCD [4,6].The upper bound is high enough to contain the PDFs of our main results.The distribution is cut off in particular cases where a lower bound for L r 5 is obtained (see later), but in such instances we find reasonable to stick with a more conservative result.
We find no credible reason to assume L r 5 smaller than zero.Our basic assumption is that L r 5 is positive, which is consistent with available determinations, which are all clearly larger than zero.
The situation is more subtle concerning L r 4 .As commented above, this constant is suppressed in the large N c limit and some results indeed find its value close to zero or even slightly negative [4,30].From the paramagnetic inequalities (20) it follows that there is a critical value of L r 4 [8,38].For the value of r we use (see (41) below) we find: Hence we use a lower bound L r 4 > −0.5 × 10 −3 .We estimate the higher order remainders statistically, based on general arguments about the convergence of the chiral series [8].Our initial ansatz is implement this by normal distributions, therefore the remainders are limited only statistically, not by any upper bound.However, our initial analysis will provide us with constraints on the higher order remainders obtained from the data for the decay constants.Subsequently, we will reuse these constraints as priors for the determination of NLO LECs L r 4 and L r 5 , thus effectively shifting the initial ansatz (40).
We use the lattice QCD average [4] for the value of the strange-to-light quark mass ratio r r = 27.23 ± 0.10. (41) Finally, the inputs for the pion and kaon decay constants are [15] F π = 92.32 ± 0.09 MeV, F K = 110.10± 0.21 MeV. (42) We use inputs from PDG [15] for the masses of the particles as well, with the experimental uncertainties being negligible compared to other sources of error.

Prediction for F η
We will employ several ways of dealing with the system of equations (13)(14)(15).At the first stage, it is possible to eliminate F 0 , L r 4 and L r 5 by simple algebraic manipulations and thus we obtain a single equation The equation depends, beyond the remainders δ F P , only on a single parameter Y and the dependence is very weak, as already noted in [8] and [11].A histogram of 10 6 numerically generated theoretical predictions is depicted in Figure 4, where the default assumptions (30-34) (illustrated in Fig. 1) and (40-42) were used.A Gaussian fit leads to a value This is an improved prediction over [11] and lies in between the values of EGMS15 (11) and EF05 (12), discussed above, while still being compatible with RQCD21 (10).As noted, this result depends only very weakly on the value of Y and thus the choice of the prior.E.g., adding the most restricting assumption (36), illustrated in Fig. 3, leads to an almost identical prediction F η = 117.7 ± 9.3 MeV (χQCD21). (45)

Higher order remainders
Next, given the weak dependence of (43) on Y , we can use RQCD21 (10), EGMS15 (11) and EF05 ( 12) as alternative inputs for F η and employ the Bayesian statistical approach to extract information about the remainders.A contour plot with confidence levels can be found in Figure 5, which leads to where ρ is the correlation These values are compatible with the prior assumption (40).We can also compare these results with the NNLO contributions for F K obtained in [3] F K /F π = 1 + 0.176 + 0.023 (BE14), (49) F K /F π = 1 + 0.121 + 0.077 (FF14). (50) As can be seen, both are positive, while EF05 (48) implies a negative remainder δ F K .It should be noted, however, that the work [3] uses a different form of the chiral expansion and thus this can only be taken as an indication that lower values of F η might be better compatible with the fits BE14/FF14.In the following, we will use the relation (43) as an additional constraint, thus effectively implementing the results (46-48).In this way we obtain improved priors for the higher order remainders, compared to the initial ansatz (40).

Extraction of L r 5
As a second step, we can algebraically eliminate F 0 and L r 4 by using equation (13), which leads to a system of two equations for F K and F η , now depending on Y , L r 5 (µ) and the remainders δ F P .We numerically generated 10 8 theoretical predictions for the kaon and eta decay constants (at µ = 770 MeV), shown in Figure 6, in comparison with the data (10) and (42).Once again, the default priors (30)(31)(32)(33)(34) has been implemented here, along with the assumptions (40-42).
Our first task is to verify the general compatibility of our ensemble of theoretical predictions with the data.As can be seen in Figure 6, this is indeed quite clearly the Horizontal lines -data from [10,15].
case, as the predictions are compatible with the data in the whole range of values.The data for F K are a little better compatible with higher values of L r 5 , while the low value of F η (RQCD21) slightly prefers lower values of L r 5 .However, it's quite evident that without additional information no values of L r 5 can be excluded at statistically significant levels.For this reason we need to employ an additional assumption about the values of the LO LECs, either (35) or (36).
In the first case, using the priors based on the additional input Y = 1.44 ± 0.32 (35) and the relation (43), depicted in Fig. 2 and Fig. 5, we obtain the following constraints on L r 5 by employing the Bayesian analysis.Figure 7 shows our main result, the probability density function for L r 5 using RQCD21 (10), in comparison with only using F K as an input.Quite clearly, incorporating F η into the analysis has a strong influence.
By approximating with a normal distribution, or alternatively putting a 2σ CL bound, we find for all the alternative inputs for F η L r 5 > 0.78 • 10 −3 at 2σ CL.
In the case of RQCD21 and EGMS15, the obtained values of L r 5 are compatible with both fits BE14/FF14 (21-22) and lattice QCD calculations (25)(26).However, for a high value of F η from EF05 (12), we obtain a lower bound for L r 5 , which is incompatible with the value from the fit FF14 (22) -L r 5 = (0.5 ± 0.07) • 10 −3 .Alternatively, using the lattice QCD input for the LO LECs (36), depicted on Fig. 3  really significant.It might seem surprising that dramatically restricting Y to a more narrow range (compare Fig. 3 vs Fig. 2) actually leads to a slightly larger uncertainty, but that is a result of a weaker of F η on L r 5 at smaller values of Y .In other words, a larger value of Y is correlated more strongly with smaller values of L r 5 .

Extraction of L r 4
As the last option, we will try to extract information on L r 4 .We will essentially repeat the procedure from the last subsection, but in this case, we will use the equation ( 14) to eliminate L r 5 , which gives us a system of two equations for F π and F η , the free variables being Z, Y , L r 4 (µ) and the remainders δ F P .Once again, we numerically generated 10 8 theoretical predictions for F π and F η , shown in Figure 9, using the default priors (30-34) (along with (40-42)).As can be seen, while the dependence on L r 4 is markedly different for the two decay constants, our theoretical predictions are compatible with the data in the whole range of values.As in the previous case, we need to employ additional information, i.e. the priors (35) or (36).
First, using the more conservative choice of priors for the statistical analysis based on (35) (Fig. 2, Fig. 5), we obtained the following probability density functions for L r 4 .Figure 10 depicts the full result using RQCD21 (10) and also the distribution given solely by F π .
Interestingly, in this case the strongest constraint is generated by the chiral expansion of F π (13) and adding F η into the analysis does not make a marked change.As can be seen, varying the input for F η does not have a significant impact and all results are compatible with both the fits BE14/FF14 (21)(22) and lattice QCD calculations (25)(26).
Using the alternative priors for the leading order LECs from χQCD21 (36) (Fig. 3), we obtain F η = 117.5 ± 9.4 MeV = (1.28 ± 0.10)F π , which is compatible with recent determinations [10,19,22].Utilizing these determinations as inputs for F η , we have applied Bayesian statistical inference to extract the values of next-to-leading order low-energy constants L r 4 , L r 5 and higher order remainders δ F K and δ Fη .L r 5 was assumed to be positive, while L r 4 > L r 4 (crit) = −0.5 × 10 −3 .By using the most recent lattice QCD data from the RQCD Collaboration [10], which provided us with the best estimate F η = (1.123± 0.035)F π , we have obtained our main result (at µ = 770 MeV): These results have used conservative estimates for the priors of the low-energy constants at the leading order (F 0 and B 0 ).Alternatively, we have used a recent computation of the leading order LECs F 0 and B 0 by the χQCD Collaboration [33] as an additional input.Though the work has not been fully published yet, it has been cited by the Flavour Lattice Averaging Group [4].In this case we have obtained: As our main conclusion, all these values are compatible within uncertainties with the most recent standard χPT fits BE14 and FF14 [3], as well as the lattice QCD computations cited by the FLAG review [4].So quite clearly, while we independently confirm the generic range of values available in the literature, an additional source of information needs to be found in order to pin down the values of the low energy constants more precisely.
However, when testing inputs for F η from phenomenology, we have found some tension if a high value of F η = (1.38 ± 0.05)F π (EF05) [19] was assumed.This lead to a negative sign of the remainder δ F K , while both fits BE14 and FF14 in [3] have positive NNLO contributions for F K .In addition, such a high value of F η produced a lower bound L r 5 < 0.78 • 10 −3 (at 2σ CL), which is incompatible with the value from the fit FF14 (L r 5 = (0.5±0.07)•10 −3 ).