# KEK-CP-082KEK Preprint 98-217January 1999 Non-perturbative determination of quark masses in quenched lattice QCD with the Kogut-Susskind fermion action

## Abstract

We report results of quark masses in quenched lattice QCD with the Kogut-Susskind fermion action, employing the Reguralization Independent scheme (RI) of Martinelli et al. to non-perturbatively evaluate the renormalization factor relating the bare quark mass on the lattice to that in the continuum. Calculations are carried out at $`\beta =6.0`$, $`6.2`$, and $`6.4`$, from which we find $`m_{ud}^{\overline{\mathrm{MS}}}(2\mathrm{G}\mathrm{e}\mathrm{V})=4.23(29)\mathrm{MeV}`$ for the average up and down quark mass and, with the $`\varphi `$ meson mass as input, $`m_s^{\overline{\mathrm{MS}}}(2\mathrm{G}\mathrm{e}\mathrm{V})=129(12)\mathrm{MeV}`$ for the strange mass in the continuum limit. These values are about 20% larger than those obtained with the one-loop perturbative renormalization factor.

The values of quark masses are fundamental parameters of the Standard Model which are not directly accessible through experimental measurements. Lattice QCD allows their determination through a calculation of the functional relation between quark masses and hadron masses. For this reason a number of lattice QCD calculations have been carried out to evaluate quark masses, employing the Wilson, clover or Kogut-Susskind (KS) fermion action .

An important ingredient in these calculations is the renormalization factor relating the bare lattice quark mass to that in the continuum. While perturbation theory is often used to evaluate this factor, uncertainties due to higher order terms are quite significant in the range of the QCD coupling constant accessible in today’s numerical simulations. A non-perturbative determination of the renormalization factor is therefore necessary for a reliable calculation of quark masses, and effort in this direction has recently been pursued for the Wilson and clover fermion actions .

The need for a non-perturbative determination of the renormalization factor is even more urgent for the KS action since the one-loop correction  is as large as 50% in present simulations. In this article we report a study to meet this need  : we calculate the renormalization factor of bi-linear quark operators for the KS action non-perturbatively using the Reguralization Independent scheme (RI) of Ref.  developed for the Wilson/clover actions. The results for the scalar operator, combined with our previous calculation of bare quark masses , lead to a non-perturbative determination of the quark masses in the continuum limit.

In the RI sheme, the renormalization factor of a bi-linear operator $`𝒪`$ is obtained from the amputated Green function,

$$\mathrm{\Gamma }_𝒪(p)=S(p)^10|\varphi (p)𝒪\overline{\varphi }(p)|0S(p)^1$$
(1)   
where the quark two-point function is defined by $`S(p)=0|\varphi (p)\overline{\varphi }(p)|0`$. The quark field $`\varphi (p)`$ with momentum $`p`$ is related to the one-component KS field $`\chi (x)`$ by $`\varphi _A(p)=_y\mathrm{exp}(ipy)\chi (y+aA)`$, where $`y_\mu =2an_\mu `$, $`p_\mu =2\pi n_\mu /(aL)`$ with $`L/4n_\mu <L/4`$ and $`A_\mu =0,1`$.

Bi-linear operators have a form

$$𝒪=\underset{yABab}{}\overline{\varphi }_A^a(y)\overline{(\gamma _S\xi _F)}_{AB}U_{AB}^{ab}(y)\varphi _B^b(y)$$
(2)   
where $`\overline{(\gamma _S\xi _F)}`$ refers to Dirac ($`\gamma _S`$) and KS flavor ($`\xi _F`$) structure , and the indices $`a`$ and $`b`$ refer to color. The factor $`U_{AB}^{ab}(y)`$ is the product of gauge link variables along a minimum path from $`y+aA`$ to $`y+aB`$. We note that $`U_{AB}(y)`$ is absent for scalar and pseudo scalar operators as these operators are local.

The renormalization condition imposed on $`\mathrm{\Gamma }_𝒪(p)`$ is given by

$$Z_𝒪^{\mathrm{RI}}(p)Z_\varphi (p)=\mathrm{Tr}[P_𝒪^{}\mathrm{\Gamma }_𝒪(p)]$$
(3)   
where $`P_𝒪^{}=\overline{(\gamma _S^{}\xi _F^{})}`$ is the projector onto the tree-level amputated Green function. The wave function renormalization factor $`Z_\varphi (p)`$ can be calculated by the condition $`Z_V(p)=1`$ for the conserved vector current corresponding to $`\overline{(\gamma _\mu I)}`$. Since the RI scheme explicitly uses the quarks in external states, gauge fixing is necessary. We employ the Landau gauge throughout the present work.

The relation between the bare operator on the lattice and the renormalized operator in the continuum takes the form,

$$𝒪_{\overline{\mathrm{MS}}}(\mu )=U_{\overline{\mathrm{MS}}}(\mu ,p)Z_{\mathrm{RI}}^{\overline{\mathrm{MS}}}(p)/Z_𝒪^{\mathrm{RI}}(p)𝒪$$
(4)   
where $`U_{\overline{\mathrm{MS}}}(\mu ,p)`$ is the renormalization-group running factor in the continuum from momentum scale $`p`$ to $`\mu `$. We adopt the naive dimensional regularization (NDR) with the modified minimum subtraction scheme ($`\overline{\mathrm{MS}}`$) in the continuum. The factor $`Z_{\mathrm{RI}}^{\overline{\mathrm{MS}}}(p)`$ provides matching from the $`\mathrm{RI}`$ scheme to the $`\overline{\mathrm{MS}}`$ scheme. These two factors are calculated perturbatively in the continuum. For our calculation of the quark mass we apply the relation (4) in the scalar channel in the chiral limit, i.e, $`1/Z_m^{\mathrm{RI}}=Z_S^{\mathrm{RI}}`$.

Our calculations are carried out in quenched QCD. Gauge configurations are generated with the standard plaquette action at $`\beta =6.0`$, $`6.2`$, and $`6.4`$ on an $`32^4`$ lattice. For each $`\beta `$ we choose three bare quark masses tabulated in Table I where the inverse lattice spacing $`1/a`$ is taken from our previous work .

We calculate Green function for $`15`$ momentum in the range $`0.038553(ap)^21.9277`$. Quark propagators are evaluated with a source in momentum eigenstate. We find that the use of such a source results in very small statistical errors of $`O(0.1\%)`$ in Green functions.

The RI method completely avoids the use of lattice perturbation theory. We do not have to introduce any ambiguous scale, such as $`q^{}`$ , to improve on one-loop results. An important practical issue, however, is whether the renormalization factor can be extracted from a momentum range $`\mathrm{\Lambda }_{\mathrm{QCD}}pO(1/a)`$ keeping under control the higher order effects in continuum perturbation theory, non-perturbative hadronization effects, and the discretization error on the lattice. These effects appear as $`p`$ dependence of the renormalization factor in (4), which should be absent if these effects are negligible.

In Fig. 1 we compare the scalar renormalization factor $`Z_S^{\mathrm{RI}}(p)`$ with that for pseudo scalar $`Z_P^{\mathrm{RI}}(p)`$ for three values of bare quark mass $`am`$ at $`\beta =6.0`$. From chiral symmetry of the KS fermion action, we naively expect a relation $`Z_S^{\mathrm{RI}}(p)=Z_P^{\mathrm{RI}}(p)`$ for all momenta $`p`$ in the chiral limit. Clearly this does not hold with our result toward small momenta, where $`Z_P^{\mathrm{RI}}(p)`$ rapidly increases as $`m0`$, while $`Z_S^{\mathrm{RI}}(p)`$ does not show such a trend.

To understand this result, we note that chiral symmetry of KS fermion leads to the following identities between the amputated Green function of the scalar $`\mathrm{\Gamma }_S(p)`$, pseudo scalar $`\mathrm{\Gamma }_P(p)`$, and the quark two-point function $`S(p)^1`$ :

$`\mathrm{\Gamma }_S(p)={\displaystyle \frac{}{m}}S(p)^1`$ (5)   
$`\mathrm{\Gamma }_P(p)={\displaystyle \frac{1}{2m}}\left[\overline{(\gamma _5\xi _5)}S(p)^1+S(p)^1\overline{(\gamma _5\xi _5)}\right]`$ (6)   
We also find numerically that the quark two-point function can be well represented by

$$S(p)^1\underset{\mu }{}\overline{(\gamma _\mu I)}\mathrm{\Sigma }_\mu ^{}(p)iC_\mu (p)+M(p)$$
(7)   
with two real functions $`C_\mu (p)`$ and $`M(p)`$, where $`\mathrm{\Sigma }_\mu ^{}(p)=\mathrm{cos}(ap_\mu )i\overline{(\gamma _\mu \gamma _5\xi _\mu \xi _5)}\mathrm{sin}(ap_\mu )`$. From (6), (7), and (3) we obtain the relations between the renormalization factors and $`M(p)`$,

$`Z_S^{\mathrm{RI}}(p)Z_\varphi (p)=M(p)/m`$ (8)   
$`Z_P^{\mathrm{RI}}(p)Z_\varphi (p)=M(p)/m`$ (9)   
In Fig. 2 $`M(p)`$ in the chiral limit obtained by a linear extrapolation in $`m`$ is plotted. It rapidly dumps for large momenta, but largely increases toward small momenta. Combined with (9) this implies that $`Z_P^{\mathrm{RI}}(p)`$ diverges in the chiral limit for small momenta, which is consistent with the result in Fig. 1.

The function $`M(p)`$ is related to the chiral condensate as follows :

$$\varphi \overline{\varphi }=\underset{p}{}\mathrm{Tr}[S(p)]=\underset{p}{}\frac{M(p)}{_\mu C_\mu (p)^2+M(p)^2}$$
(10)   
A non-vanishing value of $`M(p)`$ for small momenta would lead to a non-zero value of the condensate. Therefore the divergence of $`Z_P^{\mathrm{RI}}(p)`$ near the chiral limit is a manifestation of spontaneous breakdown of chiral symmetry; it is a non-perturbative hadronization effect arising from the presence of massless Nambu-Goldstone boson in the pseudo scalar channel.

While we do not expect the pseudo scalar meson to affect the scalar renormalization factor $`Z_S^{\mathrm{RI}}(p)`$, as indeed observed in the small quark mass dependence seen in Fig. 1, the above result raises a warning that $`Z_S^{\mathrm{RI}}(p)`$ may still be contaminated by hadronization effects for small momenta.

In Fig. 3 we show the momentum dependence of $`Z_m(\mu ,p,1/a)U_{\overline{\mathrm{MS}}}(\mu ,p)Z_{\mathrm{RI}}^{\overline{\mathrm{MS}}}(p)Z_S^{\mathrm{RI}}(p)`$ which is the renormalization factor from the bare quark mass on the lattice to the renormalized quark mass at scale $`\mu `$ in the continuum. Here we set $`\mu =2\mathrm{G}\mathrm{e}\mathrm{V}`$ and use the three-loop formula  for $`U_{\overline{\mathrm{MS}}}`$ and $`Z_{\mathrm{RI}}^{\overline{\mathrm{MS}}}`$. While $`Z_m(\mu ,p,1/a)`$ should be independent of the quark momentum $`p`$, our results show a sizable momentum dependence which is almost linear in $`(ap)^2`$ for large momenta (filled symbols in Fig. 3).

For small momenta we consider that the momentum dependence arises from non-perturbative hadronization effects on the lattice and the higher order effects in continuum perturbation theory. It is very difficult to remove these effects from our results.

Toward large momenta, however, these effects are expected to disappear. The linear dependence on $`(ap)^2`$, which still remains, should arise from the discretization error on the lattice, i.e.,

$$Z_m(\mu ,p,1/a)=m^{\overline{\mathrm{MS}}}(\mu )/m+(ap)^2Z_\mathrm{H}+O(a^4)$$
(11)   
with the constant $`Z_\mathrm{H}`$ corresponding to the mixing to dimension 5 operators on the lattice.

This relation implies that, if we take a continuum extrapolation of $`Z_m(\mu ,p,1/a)m`$ at a fixed physical momentum $`p`$, the discretization error in $`Z_m`$ is removed. This procedure also removes the $`a^2`$ discretization error in the lattice bare quark mass $`m`$ itself reflecting that in hadron masses.

The momentum $`p`$ should be chosen in the region where the linear dependence on $`(ap)^2`$ is confirmed in our results. This region starts from a similar value of $`p^23\mathrm{G}\mathrm{e}\mathrm{V}^2`$ for the three $`\beta `$ values, and extends to $`p^21.9/a^2`$, the highest momentum measured. Hence we are able to use only a rather narrow range $`3\mathrm{G}\mathrm{e}\mathrm{V}^2<p^2<6.6\mathrm{GeV}^2`$, the upper bound dictated by the value of $`1.9/a^2`$ for the largest lattice spacing at $`\beta =6.0`$.

In Fig. 4 we show the continuum extrapolation for the averaged up and down quark mass at $`\mu =2\mathrm{G}\mathrm{e}\mathrm{V}`$. Filled circles are obtained for $`p=1.8\mathrm{GeV}`$ and squares for $`p=2.6\mathrm{GeV}`$ for which the value of $`Z_m`$ is obtained by a linear fit in $`(ap)^2`$ employing the filled points in Fig. 3. The bare quark mass  is determined by a linear extrapolation of pseudo scalar meson mass squared in the Nambu-Goldstone channel $`\overline{(\gamma _5\xi _5)}`$ and that of vector meson mass in the $`VT`$ channel $`\overline{(\gamma _i\xi _i)}`$ to the physical point of $`\pi `$ and $`\rho `$ meson masses.

We observe that the continuum extrapolation completely removes the momentum dependence of the quark mass at finite lattice spacings. Furthermore the values are substantially larger than those obtained with one-loop perturbation theory (open circles for $`q^{}=1/a`$ and squares for $`q^{}=\pi /a`$ ). Making a linear extrapolation in $`a^2`$, our final result in the continuum limit is

$$m_{ud}^{\overline{\mathrm{MS}}}(2\mathrm{G}\mathrm{e}\mathrm{V})=4.23(29)\mathrm{MeV}.$$
(12)   
where we adopt the value for $`p=2.6\mathrm{GeV}`$ since this is the largest momentum accessible and the momentum dependence is negligible. This value is about $`20\%`$ larger than the perturbative estimates : $`3.46(23)\mathrm{MeV}`$ for $`q^{}=1/a`$ and $`3.36(22)\mathrm{MeV}`$ for $`q^{}=\pi /a`$. We collect the values of renormalization factor and quark masses in Table II and III.

Applying our renormalization factor to the strange quark mass, we obtain

$`m_s^{\overline{\mathrm{MS}}}(2\mathrm{G}\mathrm{e}\mathrm{V})=`$ $`106.0(7.1)\mathrm{MeV}`$ for $`m_K`$ (13)   
$`=`$ $`129(12)\mathrm{MeV}`$ for $`m_\varphi `$ (14)   
where we use $`K`$ or $`\varphi `$ meson mass to determine the bare strange mass. Results from perturbative estimation are given in Table IV.

The CP-PACS Collaboration recently reported the results  $`m_{ud}^{\overline{\mathrm{MS}}}(2\mathrm{G}\mathrm{e}\mathrm{V})=4.6(2)\mathrm{MeV}`$, $`m_s^{\overline{\mathrm{MS}}}(2\mathrm{G}\mathrm{e}\mathrm{V})=115(2)\mathrm{MeV}(m_K)`$ and $`143(6)\mathrm{MeV}(m_\varphi )`$ from a large-scale precision simulation of hadron masses with the Wilson action. Our values are 10% smaller, which may be due to the use of one-loop perturbative renormalization factor in the CP-PACS analysis.

This work is supported by the Supercomputer Project No.32 (FY1998) of High Energy Accelerator Research Organization (KEK), and also in part by the Grants-in-Aid of the Ministry of Education (Nos. 08640404, 09304029, 10640246, 10640248, 10740107, 10740125). S.K. and S.T. are supported by the JSPS Research Fellowship.