Octet baryon magnetic moments at next-to-next-to-leading order in covariant chiral perturbation theory

We calculate the octet baryon magnetic moments in covariant baryon chiral perturbation theory with the extended-on-mass-shell renormalization scheme up to next-to-next-to-leading order. At this order, there are nine low-energy constants, which cannot be uniquely determined by the seven experimental data alone. We propose two strategies to circumvent this problem. First, we assume that chiral perturbation theory has a certain convergence rate and use this as one additional constraint to fix the low-energy constants by fitting to the experimental data. Second, we fit to lattice QCD simulations to determine the low-energy constants. We then compare the resulting predictions of the light and strange quark mass dependence of the octet baryon magnetic moments by the three mostly studied formulations of baryon chiral perturbation theory, namely, the extended-on-mass-shell, the infrared, and the heavy baryon approaches. It is shown that once more precise lattice data become available, one will learn more about the convergence pattern of baryon chiral perturbation theory.


Introduction
SU(3) flavor symmetry and its breaking play an important role in our understanding of the strong interaction in the nonperturbative regime. In the limit of an exact SU(3) flavour symmetry, one can relate the magnetic moments of the octet baryons and the 0 transition to those of the proton and the neutron via the celebrated Coleman-Glashow formulae a e-mail: lisheng.geng@buaa.edu.cn [1]. Nonetheless, in nature SU(3) flavor symmetry is broken. This must be properly taken into account in order to improve the description of the baryon magnetic moments by inducing a realistic SU(3)-breaking mechanism. Chiral perturbation theory (ChPT), the low-energy effective field theory of QCD (see e.g. [2][3][4][5][6]), provides an appropriate framework to tackle this problem in a systematic fashion. However, it was noticed long ago that the leading order (LO) chiral corrections are large and tend to worsen the results, as exemplified, e.g., in [7][8][9][10]. This issue has often been used to question the validity of the SU(3) baryon ChPT altogether, see e.g. [11].
It was shown in Ref. [12] that one can achieve, however, an order by order improvement in the description of the octet baryon magnetic moments with the extended-on-mass-shell (EOMS) formulation of baryon ChPT [13,14]. 1 Although it seems that the puzzle has been solved, a natural question is what happens at the next-to-next-to-leading order (NNLO). Because of the increased number of unknown low-energy constants (LECs), a clear answer to this question has not yet been provided.
In the last two decades several calculations of the octet baryon magnetic moments in heavy baryon (HB) ChPT up to NNLO have been performed both with [9,18,19] and without [20] the inclusion of the baryon decuplet. It was shown in Ref. [20] that at NNLO the convergence of the HBChPT is quite good, contrary to the pattern exhibited at NLO. One should note that, however, in Ref. [20] the contributions of the two NNLO LECs b D ,F 6 are absorbed into the two LO ones b D,F 6 . This is legitimate as long as one works at the physical quark masses, as the primed LECs merely amount to a quark mass dependent shift of the unprimed ones that can not be disentangled.
Despite of all these studies, it remains unclear whether the convergence pattern of the HB [9], the infrared (IR) [10], and the EOMS formulation [12] observed in the description of the octet baryon magnetic moments, one of the cleanest observables, is particular to NLO and accidental, where no unknown LECs contribute, or it might be a more genuine feature of different formulations. Given that ChPT plays an indispensable role in our understanding of low energy strong interaction physics, it is of utmost importance to clarify this puzzling situation.
In this work we address this question by performing a study of the octet baryon magnetic moments at NNLO using the EOMS renormalization scheme. Since at this order the number of LECs is larger than that of the data, we first use convergence as a criterion to constrain the two LO LECs, b D,F 6 , which cannot be distinguished from the two NNLO LECs b D ,F 6 at the physical point. We then predict the light quark mass dependence of the octet baryon magnetic moments and contrast them with the results of the HB and IR formulations, and the state of the art lattice QCD simulations. Second, we fit the LECs to the lattice QCD data and then predict the strange quark mass dependence of the magnetic moments. It is shown that depending on how one determines the LECs, the predicted dependence is rather different, which could be investigated in more detail by future lattice QCD simulations. This paper is organized as follows. In Sect. 2, we present the theoretical framework. Results and discussions are given in Sect. 3, followed by a short summary in Sect. 4.

Theoretical framework
The octet baryon magnetic moments are defined via baryon matrix elements of the electromagnetic current J μ as follows whereū B and u B are Dirac spinors, m is the baryon mass, and F B 1 and F B 2 denote the Dirac and Pauli form factor, respectively. The four-momentum transfer is defined as q = p f − p i and t = q 2 . At t = 0, F B 2 (0) is the so-called anomalous magnetic moment, κ B , and the magnetic moment is μ B = κ B + q B , with q B the charge of the baryon.
In ChPT, one can calculate the baryon magnetic moments order by order, i.e., where the numbers in the superscripts are the chiral order, defined as n ChPT = 4L − 2N M − N B + k kV k for a properly renormalized diagram with L loops, N M (N B ) meson (baryon) propagators, and V k vertices from the k th order Lagrangian. Because of the large non-zero baryon mass in the chiral limit, this power counting is broken in a naive application of the MS regularization scheme [3]. To recover the power counting, several approaches have been proposed, such as the HB method [21,22], the IR approach [23] and the EOMS scheme [13,14]. In recent years, it has been shown that the EOMS scheme has some advantages because it satisfies all symmetry and analyticity constrains and converges relatively faster in certain cases, see e.g. Ref. [24] for a short review. The Feynman diagrams needed to calculate μ B up to NNLO are shown in Fig. 1. The LO contributions are provided by the following Lagrangian [25], where where m 0 denotes the baryon mass in the chiral limit, In what follows, we work in the isospin limit m l = (m u + m d )/2. Further, B 0 = | 0|qq|0 |/F 2 φ . The axial vector couplings D and F are determined from hyperon decays to be D = 0.8 and F = 0.46, and these values will be taken throughout.
The explicit expressions of the LO and NLO results can be found in Ref. [12]. In the following, we focus on the NNLO contributions. At this order one has to include oneloop diagrams with exactly one vertex from L (2) M B as well as additional tree contributions from L (4) M B . The fourth order contribution to μ B is given as (see Fig. 1): The terms contributing to μ (4,d) B collect the tree contributions with exactly one insertion from the following chiral Lagrangian [20,26], where b D 6 , b F 6 , α 1,2,3,4 and β 1 are LECs. At this order, we also have to consider double derivative operators at the meson-baryon vertex with the photon hooking on to the meson loop, see Fig. 1h. The corresponding terms of the dimension two Lagrangian read [20,26] where b 9,10,11 are LECs. They are estimated via resonance saturation in Ref. [20] and re-evaluated at m 0 = 0.94 GeV, yielding b 9 = 0.43 GeV −1 , b 10 = 0.86 GeV −1 , and b 11 = 0.45 GeV −1 . We will call these values of the LECs set I. An improved determination has been given in Ref. [10], the corresponding values are b 9 = 1.36 GeV −1 , b 10 = 1.24 GeV −1 , and b 11 = 0.46 GeV −1 . We will refer to these as set II. At NNLO, one also needs the LO chiral corrections to the baryon masses, which are provided by the following Lagrangian [25]: The two LECs b D and b F are fixed by fitting to the octet baryon mass splittings at O( The NNLO tree-level contributions can be obtained rather straightforwardly. The results are shown in the Appendix. Following the EOMS prescription to restore the power counting, we obtain the following NNLO loop results, where the coefficients C are tabulated in the Appendix and λ is the renormalization scale. We have checked that our results agree with those of Ref. [20] in the heavy mass limit up to analytical terms.

Results and discussions
At LO and NLO, the two LECs b D,F 6 can be determined by fitting to the seven experimental data μ p,n, , + , − , − , 0 . An extensive discussion of the EOMS results in comparison with the HB and IR results is given in Ref. [12]. At NNLO, however, there are nine LECs, two from the LO contribution and seven from the NNLO contribution. As a result, the experimental data alone can not uniquely determine all the nine LECs. In Ref. [20], the two NNLO LECs b D ,F 6 are absorbed by the two LO LECs, b D,F 6 , while in the present work we keep explicitly b D ,F 6 . Note that at the physical point, only the combinationsb D,F are relevant. In our numerical analysis, the decay constant and the chiral limit value of the baryon masses are chosen to be F φ = 0.108 GeV, an average of the pion, kaon, and eta decay constants, and m 0 = 0.94 GeV, following the argument in Ref. [20]. The dimensional regularization scale is set at λ = 1.0 GeV. We also have performed calculations allowing for slight variation of m 0 or λ about these values, e.g., λ = 0.9 ∼ 1.1 GeV and m 0 = 0.8 ∼ 1.1 GeV, and found that such changes have negligible effects on our results. It should be noted that in principle at the order we are working, one should use the chiral limit values for D, F, and F φ . Using the physical/average values for them will lead to differences of higher chiral order. We checked that using the chiral limit values for D, F, and F φ , such as those determined in Ref. [27,28], affect little our results and conclusions.  Table 1 Contributions of different chiral orders of the HB, IR, and EOMS schemes up to NNLO with the LECs of Table 2 Baryons EOMS IR HB

Using convergence to constrain the low-energy constants
We assume that BChPT has a reasonable convergence rate in the u, d, s three flavor sector, namely, higher order contributions are suppressed compared to lower order ones, in terms of M K / ChPT ≈ 0.5. Under this assumption, we can use the convergence criterion and the experimental data to constrain the nine LECs. More specifically, we define the convergence rate (C R ) as When fitting, we can set the convergence rate to a particular value, e.g., 0.5, and search for a combination of LECs that can satisfy such a requirement. In the left panel of Fig. 2, we show how one can fix the range of b D,F 6 by requiring C R ≤ 0.6 in the EOMS scheme. It can be seen that indeed there exist some combinations of b D,F 6 which can satisfy the requirement, namely the green shaded area. In the right panel of Fig. 2, we show the best convergence rate achievable in the EOMS scheme, in comparison with those in the HB and IR schemes. One can see that only in the EOMS scheme a convergence rate about 0.6 can be achieved, 2 while the convergence rates in the HB and IR schemes are relatively larger. This is consistent with the findings at NLO [12]. On the other  Table 2, while the shaded bands represent higher order contributions not considered (see text for details). The half filled square indicates that we have chosen a positive sign for μ 0 hand, from a typical convergence pattern of the three renormalization schemes shown in Table 1 with the corresponding LECs given in Table 2, one notices that the contributions of different chiral orders in the HB ChPT are more moderate, though not as small as one naively expects. It should be noted that in the numerical study we have taken the values given above for the LECs D, F, b D , b F , b 9 , b 10 , and b 11 . In principle, there are some uncertainties. These can originate either from the data used to fix them, or from the chiral orders at which they are determined, or from the validity of the assumption adopted, such as resonance saturation for the case of b 9 , b 10 , and b 11 . Unfortunately, there is no easy way to quantify these uncertainties. We have checked that using set II instead of set I, we find some quantitative differences but the overall trends are not affected.
One way to test the LECs determined above and also to distinguish different formulations of BChPT is to study the light quark mass dependence of the magnetic moments. Fixing the strange quark mass to its physical value with the LO ChPT relation 2B 0 m s = (2M 2 K − M 2 π )| phys. , we show the pion mass dependence of the magnetic moments in Fig. 3. The lines are obtained with the LECs tabulated in Table 2, while the bands denote higher order contributions not considered in the NNLO study. They are obtained according to Ref. [29]: where δμ (i) B are the uncertainties of chiral order i, Q = M π / ChPT , M π = 0.138 GeV, and ChPT ∼ 1 GeV is the chiral symmetry breaking scale. In addition, the requirement [29] is also satisfied. It is clear that these contributions reflect partly the convergence pattern, i.e., the slower the convergence, the larger the higher order contributions.
One can see that the three schemes display rather different pion mass dependence. For comparison, the state of the art lattice QCD results [30][31][32] are shown as well. Note that such a comparison is only meant to be qualitative, since    . 6 Predicted dependence of the octet baryon magnetic moments on the strange quark mass with the LECs determined by fitting to the lattice QCD data of Ref. [30] at O( p 4 ) with the EOMS, HB, and IR BChPT. The u/d quark mass is fixed at their physical value while the strange quark mass is proportional to 2M 2 K − M 2 π according to leading order ChPT. The bands represent higher order contributions Fig. 7 The same as Fig. 6, but with the LECs given in Table 2 these simulations are performed with the strange quark mass close to its physical value but not exactly at the physical point and furthermore current lattice QCD simulations still contain systematic uncertainties not under complete control. It is clear that the EOMS results are in better agreement with the lattice data, as corroborated by the unweighted χ 2 = (μ th − μ lattice ) 2 between the results of each scheme and the lattice QCD data [30][31][32] shown in Table 3. One should note that only at relatively large pion masses, e.g., M π > 200 or 300 MeV, one can distinguish the results from different formulations of BChPT using lattice QCD simulations.

Fitting to the lattice QCD data
Now we take a more practical attitude, forgetting about the convergence constraint, and determine the relevant LECs by fitting to the lattice QCD data directly. Among the three lattice QCD studies we considered, the one of Ref. [30] features the largest number of simulation points. Therefore, we fit to these data to determine the LECs. At O( p 3 ), one has only two LECs. Their values from the best fit together with the correspondingχ 2 are tabulated in Table 4. The predicted pion mass dependence is shown in Fig. 4, where the lattice data from Refs. [31,32] are also shown. One can see that only the EOMS formulation can describe the lattice QCD data reasonably well, consistent with the finding in the SU(2) sector [33]. Furthermore, as shown in Table 4, the LECs of the EOMS formulation determined from the fit to the lattice QCD data are similar to those determined by the experimental data [12].
The situation becomes different, however, if we fit to the lattice data with the O( p 4 ) BChPT results, all the three formulations can describe the lattice data with similar quality, as shown in Fig. 5 and Table 5. On the other hand, it seems that with the LECs determined from the best fit, the predicted μ 0 by the IR formulation is in better agreement with the experimental value (if a positive sign is taken). Nevertheless, one should note that the convergence pattern, particularly those of the IR and HB schemes, is destroyed (see Table 6).
One way to distinguish the different formulations in the present case is to study the strange quark mass dependence of the magnetic moments, shown in Figs. 6, and 7. One can see that depending on how one determines the LECs, either by fitting to the lattice QCD data or to the experimental data with the convergence constraint, the dependences on the strange quark mass are quite different. 3 It is clear that one needs more lattice data with varying strange quark mass to check which scenario is more realistic.

Summary
We studied the octet baryon magnetic moments in baryon chiral perturbation theory with the extended-on-mass-shell renormalization scheme up to next-to-next-to-leading order. We determined the low-energy constants following two different strategies, either by fitting to the experimental data with convergence as a further constraint or by fitting to lattice QCD data directly. It was shown that in the first case the extended-on-mass shell formulation seems to describe better the lattice QCD data, while in the second case, although all three formulations of baryon chiral perturbation theory can describe the lattice QCD data, they predict rather different strange quark mass dependence. Clearly more lattice QCD simulations are needed to better understand the situation and the convergence pattern of baryon chiral perturbation theory. and C (4, f ) The wave function renormalization coefficients are invariant under SU(2) transformations, therefore we give only the values of the different multiplets and C (4,h) pπ = 2(b 10 + b 11 ), C (4,h) pK = b 9 + 4b 11 , C (4,h) 11 , C (4,h) 0 π = −2b 10 + 2b 11 , and and C (4, j)