Extracting the low-energy constant Lr0 at three flavors from pion-kaon scattering

Chaitra Kalmahalli Guruswamy,Ulf-G Mei?ner,2,3 and Chien-Yeah Seng

1 Helmholtz-Institut für Strahlen-und Kernphysik and Bethe Center for Theoretical Physics,Universit?t Bonn,D-53115 Bonn,Germany

2 Institute for Advanced Simulation,Institut für Kernphysik and Jülich Center for Hadron Physics,Forschungszentrum Jülich,D-52425 Jülich,Germany

3 Tbilisi State University,0186 Tbilisi,Georgia

Abstract Based on our analysis of the contributions from the connected and disconnected contraction diagrams to the pion-kaon scattering amplitude,we provide the first determination of the only free low-energy constant atO (p4),known as L0r,in SU (4|1) Partially-Quenched Chiral Perturbation theory using the data from the Extended Twisted Mass collaboration,L0 r(μ=Mρ)=0.77 (2 0)(2 5)(7)(7)(2) · 10-3.The theory uncertainties originate from the unphysical scattering length,the physical low-energy constants,the higher-order chiral corrections,the (lattice) meson masses and the pion decay constant,respectively.

Keywords:lattice QCD,partially quenched chiral perturbation theory,pion-kaon scattering

1.Introduction

Partially-Quenched Quantum Chromodynamics and its lowenergy effective field theory (EFT),Partially-Quenched Chiral Perturbation theory (PQChPT) [1–7],are powerful tools to assist the first-principles calculations of hadronic observables using lattice QCD.They were originally created to handle the so-called partially-quenched approximation in early lattice studies,where the ‘valence’ and ‘sea’ quark masses were made distinct in order to simplify the calculation of the fermion determinant.Nowadays such an approximation has been largely abandoned,but the devised theory frameworks have found their own ways to continue being useful.In particular,it was recently realized that since additional quark flavors are introduced in PQChPT,it allows an EFT description of each individual quark contraction diagram in the lattice simulation of a given physical observable,which is very useful for getting a better handle on the noisier and computationally-expensive contractions (the so-called ‘disconnected diagrams’).This idea was applied initially to the study of the hadronic vacuum polarization [8]and the pion scalar form factor [9],and was later extended to pion-pion scattering [10,11]and the parity-odd pion-nucleon coupling constant [12].

As in any EFT,the full predictive power of PQChPT to a given order is guaranteed only when all the low-energy constants(LECs)at that order are fixed.This is a non-trivial task since some of them are not constrained by any physical experiment and can only be determined through lattice simulations.In particular,in an extended flavor sector with Nf≥4,the Cayley–Hamilton relation for 3×3 matrices used in SU(3)ChPT cannot be applied anymore,which leads to the following extra term in the PQChPT Lagrangian at O (p4):

where ‘Str’ denotes the supertrace over the extended flavor space,U is the standard exponential representation of the pseudo-Nambu-Goldstone bosons,and L0is a new LEC that does not appear in ordinary two-flavor and threeflavor ChPT.Despite appearing at the next-to-leading order (NLO),this LEC contributes to static quantities such as the pion mass and decay constant only at the nextto-next-to-leading order (NNLO),when the valence and sea quark masses are kept distinct.Based on this,[13]has determined the renormalized LECLr0in the SU (4|2)PQChPT,which is equivalent to a two-flavor ChPT in computations of physical processes.The quoted result isL0r(μ=1 GeV)=1.0 (1.1) · 10-3,using dimensional regularization.

[10,11]provided an alternative determination ofLr0in SU(4|2) using the fact that it appears in separate contraction diagrams in the ππ scattering amplitude at NLO,and therefore can be obtained from the unphysical scattering length defined through an appropriate linear combination of contraction diagrams.The advantage of this new method is that,since only light quarks are involved in the procedure,the higher-order chiral corrections that scale generically asare expected to be small so the NLO fitting is more stable than NNLO,and of course the number of unknown LECs in the former is also smaller.Using the lattice data of connected ππ correlation functions by the Extended Twisted Mass (ETM) collaboration [11,14]reportedL0r(μ=1 GeV)=5.7 (1.9) · 10-3.The significant difference with the result from the NNLO fit in [13]is yet to be understood.Resolving such a discrepancy may improve our knowledge of the chiral dynamics at low energies,or even reveal some unexpected lattice systematics that could play important roles in precision physics.

The analysis in [10,11]could be straightforwardly generalized to the three-flavor case.In particular,the socalled SU(4|1) PQChPT is of special interest because it is the simplest extension of the original three-flavor ChPT.Moreover,among all its renormalized LECs at O(p4),onlyLr0is undetermined,while all the others are identical to those in SU(3).This implies that the theory would be fully predictive atO (p4) onceL0ris fixed,and could then be used to aid the lattice studies of interesting hadronic processes involving strange quarks,such as Kπ→Kπ,ππ→and Kη→Kη scatterings,in particular channels where disconnected diagrams appear (e.g.the I=1/2 channel of Kπ scattering).

In this work,we present the first-ever numerical determination ofLr0in SU(4|1) based on the method outlined in our previous paper,[15].The main idea is that,by switching the relative sign between the two types of connected contraction diagrams that occur in the I=3/2 Kπ scattering,one obtains effectively an unphysical single-channel scattering amplitude Tβ,of which the scattering length depends onL0ratO (p4).Invoking the recent lattice data by the ETM collaboration [16],the unphysical scattering length is obtained through the usual Lüscher analysis of the discrete energy states [17],which consequently fixesLr0.This completes the SU(4|1) PQChPT Lagrangian at O(p4) and is potentially useful for all the purposes mentioned above.Besides,it serves as a prototype for future analyses of more complicated PQ-extensions of three-flavor ChPT.Such theories,after integrating out the strange quark,reduce to PQChPT with only two light flavors.This may give an independent check of value ofLr0in SU(4|2),and provide some hints towards the solution to the aforementioned discrepancy.

Figure 1.Connected contraction diagrams Ta and Tb.The thick line indicates an〈〉 contraction.

2.Formalism

As detailed in [15],the pertinent two contractions in Kπ scattering with total isospin I=3/2,as depicted in figure 1,are related to the amplitudes in the SU (4|1) PQChPT:

with s,t,u the conventional Mandelstam variables subject to the constraints+t+u=2 (MK2+Mπ2),u,d,s are the physical quarks and j denotes the additional valence quark(which comes together with a ghost quark)in PQChPT.The assigned quark masses aremu=md=mj=＜ms.For total isospin I=1/2,there is one additional scattering amplitude Tcthat is related to Tbthrough a simple crossing,Tc(s,t,u)≡Tb(u,t,s).Although such a crossing is analytically straightforward,from a lattice point of view Taand Tbare relatively easy to evaluate,whereas Tcinvolves a pair of quark propagators that start and end on the same temporal slice,and is exactly what we call a ‘disconnected diagram’.Such a diagram suffers from a low signal-to-noise ratio and represents a fundamental challenge in the first-principles study of the Kπ scattering in the I=1/2 channel.

From the above one may define three effective singlechannel scattering amplitudes:

where Tαand Tγcorrespond to the physical I=3/2 and I=1/2 scattering amplitudes,respectively,while Tβis an unphysical amplitude.The S-wave scattering lengths are defined through the threshold values of the amplitudes:

withs0=(MK+Mπ)2,t0=0 andu0=(MK-Mπ)2.Obviously,only the unphysical scattering lengthβa0can depend on the unphysical LECLr0.Its explicit expression at O (p4) reads:

Table 1.Basic parameters of the two lattice ensembles in [16],with lattice spacing a=0.0885(36)fm.The uncertainties are mainly from finite-volume effects.

Table 2.The unphysical scattering lengthevaluated from the ensembles A30.32 and A40.24.The uncertainties are mainly from thermal pollutions,which are treated using two different methods(E1 and E2).

Table 2.The unphysical scattering lengthevaluated from the ensembles A30.32 and A40.24.The uncertainties are mainly from thermal pollutions,which are treated using two different methods(E1 and E2).

μπK 0 a β Ensemble E1 E2 A30.32-0.0956(91)-0.0961(84)A40.24-0.1152(144)-0.1142(131)

3.Extraction of Lr0

The unphysical scattering lengthβa0is an essential input in our study.It is obtained from the analysis of the latticevolume-dependence of the discrete energy levels corresponding to the combination Ta–Tbthrough the standard Lüscher formula [17].In this work,we utilize the results in [16]for the Kπ system in the I=3/2 channel where the correlation functions Ca(τ) and Cb(τ) corresponding to the two contractions in figure 1 were separately calculated.The calculation was based on the Nf=2+1+1 twisted mass lattice QCD.For our analysis,we have considered the ensembles A30.32 and A40.24 for the determination of the discrete ground-state energies of the Kπ system.The basic parameters of each ensemble are summarized in table 1.

Only results that correspond to the physical combination Ca(τ)+Cb(τ) were displayed in [16]for obvious reasons.In this project,we acquired the unphysical S-wave scattering lengthβa0directly from the authors of that paper,who obtained its value through an unpublished analysis of the volume-dependence of the discrete energy levels extracted from the unphysical combination Cβ(τ)=Ca(τ)-Cb(τ)[18]:first,the energy shiftδEβ=EβπK-Mπ-MKwas obtained as a function of the lattice size L from the exponential behavior of Cβ(τ) at large Euclidean time τ.The scattering lengthβa0at infinite volume was then computed using the single-channel,Taylor-expanded Lüscher formula

as in equation (14) of [16].Here,μπKis the reduced mass of the Kπ system and c1,2are known coefficients.The final outcomes are summarized in table 2.The main uncertainty ofa0βis systematic,which comes from the unwanted timedependent contributions at finite τ (i.e.the ‘thermal pollutions’).Their effects were studied using two different methods labeled as E1 (weighting and shifting) and E2(dividing out the pollution),respectively [19].

To solve forLr0using equation(5),we further require the values of all the physical LECs.We took their values at μ=Mρ=770 MeV from [20]:

With all the above,we may now computeLr0straightforwardly.The outcome at μ=Mρreads:

where the uncertainties come fromβa0,the physical LECs,the higher-order ChPT corrections,the (lattice) meson masses and Fπ,respectively.In particular,the higher-order ChPT corrections are estimated by multiplying the central value with the usual chiral suppression factorWe see that the values ofL0rfrom all four determinations are consistent with each other within the error bars,so we may simply quote the number with the smallest theory uncertainty,namely the one from the ensemble A30.32 with method E2:

Finally,we comment on the relation between this result and the corresponding LEC in SU(4|2).In principle,relations between LECs in the two-flavor and three-flavor ChPT can be obtained by integrating out the strange quark in the latter.Possible PQChPT extensions of an ordinary two-flavor ChPT are SU (3|1),SU (4|2),SU (5|3)...,etc,but only SU (4|2)onwards possess an L0-dependence at tree-level as it requires at least four fermionic quarks.Similarly,possible PQChPT extensions of a three-flavor ChPT are SU (4|1),SU (5|2),SU (6|3)....In [15]we chose the simplest version which is SU (4|1).After integrating the strange quark,it reduces to SU (3|1) that does not depend on L0at tree-level.Therefore,it is not possible to discuss the matching between L0in two-and three-flavors based on the theory setup in [15].For that,one would have to repeat the calculations using a larger graded algebra,such as SU (5|2).This is of great interest because it may provide new insights into the apparent disagreement between the determination ofLr0at SU(4|2) from NLO and NNLO.However,it goes beyond the scope of this work and will be carried out in follow-up studies.

4.Summary

In this work,we have for the first time determined the unphysical LECL0r(Mρ) in the simplest PQ-extension of the three-flavor ChPT through an NLO analysis of contraction diagrams in Kπ scattering,L0r=0.77 (3 3) · 10-3.Utilizing the precise data from the ETM collaboration for Kπ scattering in the I=3/2 channel,we control the absolute uncertainty in this LEC to 3.3×10-4,which is better than the previous determinations of the corresponding LEC for two flavors in[13]that made use of an NNLO fitting,and in the NLO fitting to the ππ scattering amplitudes in [10,11]that depends on more unknown LECs.The major sources of uncertainty in this study are the systematic errors in the lattice extraction of the unphysical scattering lengthβa0,as well as the physical LECsL-r18.Our work completes the PQChPT Lagrangian at O(p4) and prepares it for future applications in studies of interesting hadronic observables involving strange quarks,in synergy with lattice QCD.

Acknowledgments

We are very grateful to Ferenc Pitler for making the ETM collaboration data available to us and for his detailed explanations concerning these.We thank Hans Bijnens for a useful communication.This work is supported in part by the DFG(Projektnummer 196 253 076—TRR 110)and the NSFC(Grant No.11 621 131 001) through the funds provided to the Sino-German CRC 110‘Symmetries and the Emergence of Structure in QCD’,by the Alexander von Humboldt Foundation through the Humboldt Research Fellowship,by the Chinese Academy of Sciences (CAS) through a President’s International Fellowship Initiative (PIFI) (Grant No.2018DM0034),by the VolkswagenStiftung(Grant No.93 562),and by the EU Horizon 2020 research and innovation programme,STRONG-2020 project under grant agreement No.824 093.