The pion mass and decay constant at three loops in two-flavour chiral perturbation theory

A calculation of the pion mass and decay constant at NNNLO in two-flavour chiral perturbation theory is presented. The results are cross-checked by using both the exponential and square root parametrizations of the Goldstone matrix field, as well as by comparing to the known leading log coefficients of the two quantities. A small numerical study of the quark mass dependence is performed, and for a physical quark mass there is good agreement with lower order results.


JHEP11(2017)181
general method to NNLO is described in detail in [5]. We extend it to one order higher in the expansion.
The motivation behind this work is twofold. The expressions themselves are of intrinsic interest but in so-called hard-pion ChPT it was argued that mass logarithms could be calculated also for pions with hard momenta. This was checked at two-loop order in ChPT [9]. At three-loop order it was found that the chiral mass logarithm does not agree with the prediction of [9] in [10]. This work is a first step towards checking the results of [10] and possibly being able to correct and extend the arguments of [9] towards a full proof.

Mass in the O(N ) φ 4 model
The Lagrangian for the O(N ) φ 4 model is given by φ is a vector of N real fields φ a and f is the external current. We have indicated here the higher order terms with l i , c i and d i as well and have used the equations of motion (or field redefinitions) to discard the higher order terms in the coupling to the external current f . Up to two-loop order the counterterms are given by 1 The mass is defined as the pole of the two point function, see e.g. the discussion in [6], as

JHEP11(2017)181
where Σ is the sum of one-particle irreducible diagrams. The relevant diagrams are shown in figure 1 when neglecting those involving vertices with more than four legs. p 2 corresponds to terms with 1, p 4 with l i , p 6 with c i and p 8 with d i in (2.1). The physical mass M φ is given by the solution of We write the mass as and the self-energy as Where we used the subscript 4,6,8 for NLO, NNLO and NNNLO respectively. With this expansion we can solve for the mass perturbatively. We evaluate diagrams at p 2 = M 2 φ as a perturbative expansion away from M 2 . This is why derivatives of the self-energy show up. Taking into account that here ∂ 2 Σ 4 /(∂p 2 ) 2 = 0 we get (2.8) All nonlocal divergences cancel as they should and we get a finite result by setting

JHEP11(2017)181
The full result for the mass at three-loop order is This result can be checked in a number of ways. The nonlocal divergences cancelled as they should. The terms leading in N can be derived using a gap equation similar to what was done for the nonlinear sigma model in [11,12]. The renormalization group equations are known to five loop order [13], these can be used to check the 1/ε terms in (2.9). All checks are satisfied.

Chiral perturbation theory
The effective Lagrangian in ChPT is expanded in powers of p 2 as L = L 0 + L 4 + L 6 + L 8 + . . . .
The fields l µ , r µ , s and p are the usual N f × N f external fields of ChPT. F and B are the two low-energy-constants (LECs) at leading order for the two-flavour case. The nextorder Lagrangian was classified in [2]. The NNLO Lagrangian can be found in [14]. The NNNLO Lagrangian L 8 is at present not known, but there will be one combination of p 8 LECs contributing to the mass and another to the decay constant, we will call these combinations r M 8 and r F 8 , respectively. The NLO and NNLO low-energy-constants (LECs) are conventionally denoted as l i and c i , respectively. The divergent parts needed to one- [2] and two-loop order [15] are known JHEP11(2017)181 in general and the equivalent formulas to (2.2) can be found there. For later convenience we introduce the lowest order pion mass where 2m = m u + m d . In the remainder we will work in the isospin limit with m u = m d .

The calculation and checks
The diagrams contributing are shown in figure 1. The Feynman diagrams are programmed in FORM [16]. The derivatives w.r.t. p 2 needed are obtained by taking the derivative diagram by diagram at this stage. Then the expressions are rewritten in integrals. These are reduced to a set of master integrals. This is done using integration-by-parts and Lorentz invariance identities through a Laporta algorithm. We have used the program Reduze [17] for this. The resulting master integrals are all known to the order in ε required and we quote them in appendix A.
A three-loop calculation needs a large number of checks. We have checked that the nonlocal divergences cancel, that we reproduce the known two-loop and leading logarithm results. As a final check we use two different parametrizations for u in terms of a traceless Hermitian 2 × 2 matrix Φ, namely the exponential, u = exp i √ 2F Φ, and a square root, The diagrams are quite different in these two parametrizations but the final result must of course be the same. All checks are satisfied by our results. Since essentially the same programs were used for φ 4 the checks discussed in section 2 are another partial check on our main results.

JHEP11(2017)181 5 The pion mass and decay constant
The physical mass is defined as the pole of the two-point function (2.4) with Σ(p 2 ) the self-energy. The physical pion mass M 2 π is then found as the solution of The decay constant is defined through the relation and is calculated using diagrams of the same topology as those in Σ (the only difference is that one of the external legs corresponds to the axial current). For the decay constant one also needs to calculate the wave function renormalization factor Z defined as the residue of the propagator in (2.4), i.e.
π . The physical pion mass and decay constant can be written in expanded form (5.4)

Mass
We can solve (5.1) perturbatively and obtain here we used the fact that ∂ 2 Σ 4 /(∂p 2 ) 2 = 0. In order to obtain a finite result we need the subtraction r M 8 is the combination of p 8 LECs that contributes to the mass.

JHEP11(2017)181
This is a single scale problem and only logarithms of the mass scale show up, the expression is thus fairly compact. We use the abbreviations The results can be written in the form The coefficients are where ζ k is the Riemann-Zeta function. The leading log coefficient a M 33 agrees with [11,12]. The LECs l r i are well known but the c r i less well. In [5,15] combinations of the p 6 LECs appearing at p 6 in ππ-scattering, the mass and decay constant were defined, r 1 , . . . , r 6 , r M , r F , and numerical estimates using resonance saturation were done in [5].

Decay constant
For the decay constant everything is analogous except that we need to evaluate the diagrams with one leg replaced by the axial current and take into account the wave function renormalization factor Z. Denoting the sum of one-particle-irreducible diagrams of the axial current as A(p 2 = M 2 π ) = A 4 (p 2 ) + A 6 (p 2 ) + A 8 (p 2 ) the expression for the decay constant is (normalized to 1 at lowest order) Putting in the expanded expressions for Σ and A and using ∂ 2 Σ 4 /(∂p 2 ) 2 = ∂A 4 /∂p 2 = 0, we obtain with all right hand sides evaluated at p 2 = M 2 .

JHEP11(2017)181
The full results for the coefficients are, with r q F 8 = (16π 2 ) 3 r r F 8 , Also here the leading log coefficient, a F 33 , agrees with [12]. We can similarly to the previous subsection rewrite the results in terms of the combinations r q i i = 1, . . . , 6, M, F . The rewriting is not fully possible here.

Numerical study: mass dependence
Now that the analytic forms of the mass and decay constant have been obtained, we can do a first numerical analysis of the mass dependence. We present only results for one choice of input parameters to give an impression of the size of the NNNLO correction.

JHEP11(2017)181
The expansions given in the previous section correspond to an expansion expressed in terms of the lowest order mass M and decay constant F in the form There are many ways to rewrite this expansion but the second most standard version is the inverse, namely where ξ = M 2 π /(16π 2 F 2 π ) and L π = log(M 2 π /µ 2 ). The analytic expressions of the a F,M ij are in the previous section and the b F,M ij can be found in appendix C. These are often referred to as the x-and ξ-expansion, see e.g. [18]. As

Conclusions
In this paper we presented the calculation of the NNNLO contributions to the pion mass and decay constant in the isospin limit of two-flavour ChPT. We also calculated the mass in the O (N )

A Master integrals
Below, the master integrals needed for the calculation 2 are listed. To simplify the expressions of the master integrals below, we define C ε, M 2 = (4π) ε Γ(1 + ε)M −2ε .
(A.1) Also, we denote the external momentum as p with p 2 = M 2 , as is relevant for the quantities considered here.

A.1 One loop
We only need one one-loop integral, the well-known tadpole:

A.2 Two loops
All needed two-loop integrals can be reduced to the equal mass on-shell sunset or products of tadpoles. The sunset result is, see e.g. [20] 1