arxiv: v2 [hep-lat] 8 Sep 2017

Transcription

1 A finite-density transition line for QCD with 695 MeV dynamical fermions Jeff Greensite and Roman Höllwieser, 3 Physics and Astronomy Department, San Francisco State University, San Francisco, CA 943, USA Department of Physics, New Meico State University, Las Cruces, NM, 883-, USA 3 Institute of Atomic and Subatomic Physics, Vienna University of Technology, Operngasse 9, 4 Vienna, Austria (Dated: September, 7) arxiv:78.83v [hep-lat] 8 Sep 7 We apply the relative weights method to SU(3) gauge theory with staggered fermions of mass 695 MeV at a set of temperatures in the range 9 T 6, to obtain an effective Polyakov line action at each temperature. We then apply a mean-field method look for phase transitions in the effective theory at finite densities. The result is a transition line in the plane of temperature and chemical potential, with an endpoint at high temperature, as epected, but also a second endpoint at a lower temperature. We cannot rule out the possibilities that a transition line reappears at temperatures lower than the range investigated, or that the second endpoint is absent for light quarks. I. INTRODUCTION The effective Polyakov line action (PLA) of a lattice gauge theory is the theory which results from integrating out all of the degrees of freedom of the theory, subject to the condition that the Polyakov lines are held fied, and it is hoped that this effective theory is more tractable than the underlying lattice gauge theory (LGT) when confronting the sign problem at finite density. The general idea was pioneered in [], and the derivation of the PLA from the underlying LGT has been pursued by various methods, e.g. [ 4]. The relative weights method [5] is a simple numerical technique for finding the derivative of the PLA in any direction in the space of Polyakov line holonomies. Given some ansatz for the PLA, depending on some set of parameters, we can use the relative weights method to determine those parameters. Then, given the PLA at some fied temperature T, we can apply a mean field method to search for phase transitions at finite chemical potential µ. This is the strategy which we have outlined in some detail in [6], where some preliminary results for finite densities were presented. The relative weights method has strengths and weaknesses; on the positive side the approach is not tied to either a strong coupling or hopping parameter epansion, and the non-holomorphic character of the fermion action is irrelevant. The main weakness is that the validity of the results depends on a good choice of ansatz for the PLA. We have suggested, for eploratory work, an ansatz for the PLA inspired first by the success of the relative weights method applied to pure gauge theories [5], and secondly by the form of the PLA obtained for heavy-dense quarks. at a lower temperature. Whether a transition line reappears at still lower temperatures, outside the range we have investigated, or whether the second transition point disappears for lighter quarks, or whether this second transition point is instead indicative of some deficiency in our ansatz for the PLA, remains to be seen. In the net section we briefly review the relative weights method and associated mean field technique at finite density, referring to our previous work [6] for some of the technical details. Section 3 contains our results, followed by conclusions in section 4. A result for a recently introduced observable ξ/ξ nd [], which is sensitive to ecitations above the lowest lying ecitation, is presented in an appendi. II. RELATIVE WEIGHTS It is simplest to work in temporal gauge, where we can fi the timelike links to the identity everywhere ecept on one timeslice, say at t =, on the periodic lattice. In this gauge the timelike links at t = are the Polyakov line holonomies, which are held fied, and the PLA S P for an SU(N) gauge theory is defined by [ ] ep S P [U,U ] { = DU (,)DU k DψDψ δ[u U (,)] } e S L, where S L is the lattice action for an SU(N) gauge theory with dynamical fermions, and we also define the Polyakov line () In this article we follow up on the work in ref. [6] to obtain a tentative transition line in the µ T plane for SU(3) gauge theory with dynamical staggered unrooted quarks of mass 695 MeV. It is generally believed that this line has an endpoint at high temperatures, and this is what we find. A second, unepected finding is that there is also an endpoint In fact some conjectured QCD phase diagrams do contain a second endpoint (see, e.g. [7, 8]), but this is based on the idea of quark hadron continuity at N f = 3 flavors [9], and it is not clear that a similar argument would apply in our case, with unrooted staggered fermions.

2 P = N TrU. In this article we use the standard SU(3) Wilson action for the gauge field, and unrooted staggered fermions as the dynamical matter fields. Given the PLA at µ =, the PLA at finite µ is obtained by the simple replacement S µ P [U,U ]=Sµ= P [e N t µ U,e N t µ U ], () where N t is the lattice etension in the time direction, with µ in lattice units. Let us consider some path through the space of Polyakov line holonomies U (λ) parametrized by λ. The relative weights method allows us to compute the derivative ds P /dλ anywhere along the path. Let U and U represent the Polyakov line holonomy field at parameter λ ± λ respectively. Defining S L,S L as the lattice actions with the Polyakov line holonomies fied to U,U respectively, and S P = S P [U ] S P[U ], then in temporal gauge we have by definition e S DUk P Dφ e S L = DUk Dφ e S L DUk Dφ ep[s L = S L ]es L DUk Dφ e S L = ep[s L S L], (3) where... indicates that the epectation value is to be taken in the probability measure e S L DUk Dφ e S L. (4) The epectation value in the last line of (3) can be calculated by standard lattice Monte Carlo, only holding fied timelike links at t =. From this calculation we find the derivative ds P /dλ S P / λ. The SU() and SU(3) gauge groups are special in the sense that the trace P determines the holonomy U. So we epand P = a k e ik, (5) k and compute the derivatives S P / a k with respect to a set of specific Fourier components a k, keeping the remaining components fied to values taken from a thermalized configuration. Our ansatz for the PLA is e S P = det[+he TrU ] p det[+he TrU ]p [ ] ep P P y K( y), (6) y where p = for unrooted staggered fermions, and p = N f for Wilson fermions, where N f is the number of flavors. The determinant factors det[+he TrU ] = +he TrU + h e TrU + h 3 e 3, (7) det[+he TrU ] = +he TrU + h e TrU + h 3 e 3 (8) are motivated by the PLA for heavy dense quarks [, ] in which h=(κ) N t, with κ the hopping parameter for Wilson fermions, or κ = /m for staggered fermions. In our ansatz, h becomes a fit parameter. The part of the action involving the kernel K( y) is motivated by previous successful treatments [5] of pure gauge theory and gauge-higgs theory. All in all, using where,y P K( y)p y = a k a k K(k), (9) k K( y)= L 3 k K(k)e k ( y) on an L 3 three-space volume, and choosing p=, we have () S P [U ]= k { a k a k K(k)+ log(+he Tr[U ]+h e Tr[U ]+h3 e 3 ) } +log(+he Tr[U ]+h e Tr[U ]+h3 e 3 ). () The Fourier transform K(k) of K( y) is determined numerically at µ =. For k=, ( ) SP L 3 a R a =α { = K()α+ (3h+3h ) L 3 } Q + c.c, () where Q = +3hP + 3h P + h 3. (3) If h, then dropping terms of O(h ) and higher the deriva-

3 3 tive simplifies to L 3 ( SP a R ) = K()α+ 6h, (4) a =α where the R superscript refers to the real part. For k, again dropping terms of order h and higher, ( ) SP L 3 a R = K(k)α. (5) k a k =α The derivatives on left-hand sides of (4) and (5) are determined by the method of relative weights at a variety of α. By plotting those results vs. α and fitting the data to a straight line, K(k) is determined from the slope. The h parameter can in principle be determined by the y-intercept of the K() vs. α data, but a better method is to choose h by requiring that P computed in the PLA at µ = agrees with P computed in the underlying LGT. K(R) β=5.7, ma=.6. R line fit poly fit.66/r 4 FIG.. Comparison of the position space K(R) vs. R obtained using the double-line fit and the polynomial fit. The deviation of the two procedures indicates a finite-size effect beyond R 5 (see tet). K(k) β=5.7, ma=.6 data fit fit k L (a) ous work [5, 6], we fit the data for K(k) at k by a double straight-line fit where K f it (k)={ c c k L k L k d d k L k L k, (6) k L = 3 i= sin (k i /) (7) is the lattice momentum, and carry out the inverse Fourier transform, but taking K() from the data rather than the fit. The second method is to fit the data to a third-order polynomial K(k) β=5.7, ma=.6 data polynomial fit k L (b) FIG.. Fits to the relative weights results for K(k) vs k L, obtained at β = 5.7, ma=.6, N t = 6: (a) double line fit (6); (b) third order polynomial fit (8). Since K(k) is determined at only a small set of k values on a 6 3 volume, it is necessary to fit the data to some analytic epression in order to carry out the inverse Fourier transform to K( y). We have tried two procedures: First, as in previ- K f it (k)=c + c k L + c k L + c 3k 3 L. (8) These two fits to the same set of relative weights data for K(k), obtained at β = 5.7, ma=.6, N t = 6, are illustrated in Fig.. On taking the inverse Fourier transform, both methods agree reasonably well for the position-space K(R) at moderate values of R, in this instance up to R 4 on a 6 3 volume, and in this interval K(R) falls off as /R 4, as seen from a straight-line fit on the log-log plot of Fig.. At larger R the two fits disagree, with K(R) either going negative and then oscillating (double line procedure), or flattening out (polynomial procedure). We take this as an indication of a finite size problem, which should be studied further by going to larger lattice volumes. Unfortunately, as we are dealing with dynamical fermions, we are limited by computer resources, since carrying out the relative weights calculation with dynamical fermions is very costly on the LGT side. We therefore take as a guiding principle the fact that K(R) cannot continue to fall with a power law on a large lattice, because this would inevitably result in a power law falloff, rather than an eponential falloff, in the connected Polyakov correlator. On these grounds we impose a simple cutoff, K(R) = for R > R cut 4 6, at the value of R cut where the inverse Fourier transform of the polynomial fit begins to deviate from

4 4 the /R 4 falloff. With K(R) determined by the procedure just outlined, we find h by simulating the PLA in eq. () at µ = with a series of trial values of the h parameter, until P computed in the PLA agrees with the value computed in the LGT up to error bars. This is not entirely straightforward. As we found in our earlier work [6], the highly non-local term containing K( y) in the action leads to metastable states which, in a lattice simulation of the PLA, depend on the initialization and persist for thousands of iterations. A cold start, which sets P() = initially, generally leads to the system in the deconfined phase, with a large epectation value P, even at h =. This is in strong disagreement with the LGT. Instead we initialize the system at P() =, which then stays in the confined phase. Having obtained h, we can then go on to compute the Polyakov line correlator at µ = G( y )= P()P(y), (9) in the PLA, and compare with the corresponding correlator obtained in the underlying LGT. A. Mean field at µ The PLA still has a sign problem at µ, which we deal with via a mean field method []. We summarize here only the essential points; a detailed derivation may be found in the cited references. Starting from the partition function of the effective theory Z = du D (µ,tru,tru )e S S = 9 K( y)tru TrU y, () where D (µ,tru,tru ) is the product of the determinant factors (7) and (8), we rewrite S = J (vtru + utru ) uvj V +a Tr[U ]Tr[U ]+E, () where V = L 3 is the lattice volume, and we have defined E = (TrU u)(tru y v) K( y), 9 (y) J = 9 K(), a = K(). () 9 Parameters u and v are to be chosen such that E can be treated as a perturbation, to be ignored as a first approimation. In particular, E = when u= TrU, v= TrU, (3) and these conditions turn out to be equivalent to the stationarity of the mean field free energy. The mean field approimation is obtained, at leading order, by dropping E, in which case the partition function factorizes, and can be solved analytically as a function of u and v. After some manipulations (cf. []), one finds the mean field approimations u,v to TrU and TrU respectively, by solving the pair of equations u G G A = and v G =, (4) G B where A=J v, B=J u. The epression G(A,B) is given by ( G(A,B)=D µ, A, ) [ D s i j B I [ ] AB],(5) det s= where D s i j is the i, j-th component of a matri of differential operators { D s i j = Di, j+s s D i+ s, j s<, ( i j B) i j D i j = ( j i, (6) A) i< j The mean field free energy density f m f and fermion number density n are f m f T = J uv logg(a,b) n= G G µ. (7) The stationarity conditions (4) may have more than one solution, and here it is important to take account of the eistence of very long lived metastable states in the PLA. The state at µ = which corresponds to the LGT is the one obtained by initializing at P =, and in the mean field analysis this is actually not the state of lowest free energy (its stability in a Monte Carlo simulation is no doubt related to the highly nonlocal couplings in the PLA). By analogy, at finite µ we look for solutions of (4) by starting the search at u = v =, regardless of whether another solution eists at a slightly lower f m f. III. RESULTS Our LGT numerical simulations were carried out in SU(3) lattice gauge theory with unrooted staggered fermions on a lattice. For scale setting we have taken the lattice spacing from the Necco-Sommer epression [3] a(β)=(.5 fm)ep[ (β 6) (β 6).448(β 6) ]. (8) We take the quark mass in lattice units to be ma =.6 at β = 5.7. This corresponds to a mass of m = 695 MeV, and temperature T = 88 MeV in physical units. We keep the physical mass and the etension N t = 6 in the time direction fied, and vary the temperature by varying the lattice spacing, i.e. by varying β.

5 5. LGT PLA LGT PLA. <P + P(r)>.. <P + P(r)> r r (a) β = 5.7, ma=.6 (b) β = 5.75, ma=.536 LGT PLA LGT PLA. <P + P(r)> <P + P(r)> r r (c) β = 5.8, ma=.48 (d) β = 5.85, ma=.435 FIG. 3. Comparison of Polyakov line correlators P()P (r) vs. r at µ = and m=695 MeV obtained from simulations of the effective PLA and underlying LGT. The LGT lattice volume is 6 3 6, and different temperatures are obtained by varying the lattice coupling: (a) T = 88 MeV, (b) T = MeV, (c) T = 35 MeV, (d) T = 6 MeV. An important check is to compare the Polyakov line correlators G(R) obtained in the LGT and the PLA at µ =. Typical results are shown in Fig. 3 for β = 5.7,5.75,5.8,5.85, corresponding to temperatures T = 88,,35,6 MeV respectively. Correlators in the LGT and the PLA are in very good agreement. Given the PLA at µ =, the Polyakov line epectation values TrU and TrU at finite µ are calculated by the mean field method outlined above, with a sample of our results displayed in Figs. 4. A discontinuity in a plot of TrU, TrU vs. µ is the sign of a transition at finite density, and conversely the absence of any discontinuity indicates the absence of any transition. When a transition occurs at some value of chemical potential µ, then there is a second transition at some µ > µ. However, while the first transition occurs at some relatively low density (in lattice units) on the order of n., the second transition always occurs at a density close to the saturation value, which for staggered fermions is n = 3. An eample of density n vs. µ at β = 5.75 is shown in Fig. 5; transitions occur at the sharp jumps in density. Since the saturation value is a lattice artifact, we do not attach much physical significance to the transition at µ = µ. Our full set of results is shown in Table I. Note the comparison between P(LGT), which is the Polyakov line computed in the LGT at µ = (and used to determine h by numerical simulation of the PLA), with P(mfd), which is the corresponding estimate for the Polyakov line obtained from a mean field approimation to the PLA at µ =. The very close agreement, in what amounts to a spin system in D=3dimensions, can be attributed to the highly non-local nature of the PLA. Mean field treatments work best when each degree of freedom is coupled to many other degrees of freedom. For nearestneighbor couplings, this generally means that the treatment is only accurate in higher dimensions. But it seems that the relatively long range of K(R), which results in each Polyakov line being coupled to many other lines, serves the same purpose as high dimensions in a nearest-neighbor theory. This close agreement of P(LGT) with P(mfd) gives us some confidence in the mean field treatment also at µ >. We should also note the good agreement found in ref. [] between Langevin simulations at finite density in a spin system of this type, with the mean field treatment we have discussed. From the µ vs. T data shown in the table, we plot a transition line in the µ T plane for staggered, unrooted quarks of mass 695 MeV, Fig. 6. This figure is the main result of our paper, and holds for the temperature range 9 T 6 MeV

6 6.4. <TrU> <TrU + >.8.6 <TrU> <TrU + >.4 <TrU>, <TrU + > <TrU>, <TrU + > (a) β = 5.6, ma=.767 (b) β = 5.7, ma=.6 <TrU>, <TrU + > <TrU> <TrU + > <TrU>, <TrU + > <TrU> <TrU + > (c) β = 5.75, ma=.536 (d) β = 5.8, ma=.48 FIG. 4. Eamples of mean-field calculations of TrU and TrU at finite µ. There are phase transitions in the two mid-range temperatures, but no transitions at the highest and lowest temperatures. number density is a critical endpoint. The fact that there is a critical endpoint at high temperature was epected. What was unepected is the eistence of a second critical endpoint at a lower temperature of T 58 MeV. We cannot rule out the possibility that a high density transition reappears at some temperature lower than the lowest temperature (9 MeV) that we have considered. Or perhaps the second critical endpoint goes away for light quarks. These possibilities we reserve for later investigation FIG. 5. Mean field result for fermion number density per lattice site, vs., at β = 5.75 corresponding to T = MeV. Two transitions are visible. Note that the transition at higher chemical potential occurs close to the saturation value, which is n = 3 for staggered fermions. Since saturation is a lattice artifact, this second transition is probably unphysical. that we have investigated. We see that the phase transition line eists to an upper temperature of T MeV, where there IV. CONCLUSIONS We have found a first-order phase transition line for SU(3) gauge theory with dynamical unrooted staggered fermions of mass 695 MeV, by the method of relative weights combined with mean field theory, in the plane of chemical potential µ and temperature T. This line has two critical endpoints, at T = 58 MeV and T = MeV. We offer this result with reservations. There is certainly a degree of arbitrariness in our approach, particularly in the choice of ansatz for the PLA. We have used a product of local determinants for the µ-dependent part of the PLA, and a non-local bilinear form for the µ-independent part. This form

7 7 β T[MeV] a[fm] ma P (LGT) P (mfd) R cut h µ /T µ /T TABLE I. Lattice coupling β for the LGT and parameters h, R cut in the corresponding PLA, with the LGT on a lattice and dynamical staggered fermions with m q = 695 MeV. Also shown are the temperature T and lattice spacing a in physical units, the quark mass ma in lattice units, the LGT and mean field results (at µ = ) for the Polyakov lines, and chemical potentials µ /T, µ /T of the finite density transitions determined by the mean field method described in the tet. 3 5 phase transition line critical endpoints JG s research is supported by the U.S. Department of Energy under Grant No. DE-SC368. RH s research is supported by the Erwin Schrödinger Fellowship program of the Austrian Science Fund FWF ( Fonds zur Förderung der wissenschaftlichen Forschung ) under Contract No. J345-N7. T [MeV] 5 Appendi: The ξ/ξ nd observable µ [MeV] FIG. 6. The phase transition line in the µ T plane for staggered unrooted fermions of mass m = 695 MeV, obtained by the method described in the tet. has given us an ecellent match to the Polyakov line correlator at µ =, computed in the LGT. But there is no guarantee that a more complicated form, involving e.g. Polyakov lines in higher representations, trilinear couplings, etc., is not required at high densities. It would be interesting to probe the eistence and possible importance of such terms, supplementing the method of relative weights and mean field with, e.g., the method of inverse Monte Carlo [4, 4] and/or the approach of Bergner et al. [3]. It would also be interesting to see what happens to the second critical endpoint in a simulation with lighter quarks, or whether the epected transition line reappears, for m = 695 MeV, at temperatures below 9 MeV. We reserve these questions for later study. Gauge theories have a complicated spectrum, and correlators at intermediate distances should be sensitive to more than just the lowest-lying ecitation. Caselle and Nada [] have introduced an interesting observable which tests for the presence, in an effective Polyakov line action, of a spectrum of ecitations contributing to two-point correlators. We briefly review their idea in this section, as implemented at µ =. Denoting P = P,y,z, define the average of Polyakov lines in a plane on an L 3 lattice volume P(z)= L P,y,z, y and consider the connected correlator G(τ)= P()P (τ) P (A.) (A.) We etract the correlation length ξ from a best fit of G(τ) to the form G(τ) c (e τ/ξ + e (L τ)/ξ) (A.3) which is a form we epect to be true asymptotically for sufficiently large τ, L τ. But suppose that in fact G(τ) is more accurately epressed as a sum of eponentials, reflecting the ACKNOWLEDGMENTS We adhere in this appendi to the notation of [], but G(τ) should not be confused with the usual Polyakov line correlator defined in eq. (9).

8 8 eistence of ecited states, i.e. for τ L/ G(τ) c i e τ/ξ i i (A.4) Let G (τ) be G(τ) in the L limit. Following [] we consider ξ nd = τ= τ G (τ) τ= G (τ) (A.5) If we approimate the sums over τ by integrals, and assuming that G (τ) has the form (A.) for all τ, we find that ξnd = i c i ξi 3 i c i ξ i (A.6) Now suppose that the sum is dominated by the first term, with all other terms negligible. In that case ξ defined in (A.3) equals ξ, and we have ξ/ξ nd =. In the Ising model in D = 3 dimensions, the deviation from unity is on the order of a few percent. Conversely, if the ratio ξ/ξ nd is greater than one, then this is evidence that there is a spectrum of ecited states, not just the lowest lying ecitation, that make a non-negligible contribution to the sum. For SU() pure gauge theory, lattice Monte Carlo simulations well away from the deconfinement transition show a much more substantial deviation, on the order 4-5% []. G(τ) τ data fit FIG. 7. G(τ) defined in (A.) vs. τ for the PLA corresponding to the LGT at β = 5.65 and a lattice volume. On the lattice we have lattice periodicity, and we cannot take the L limit to get G (τ). The proposal is instead to approimate (A.5) by ξnd = ( τ ma t= τ G(τ)+ τ=τ ma + τ G(τ ma )ep ( τ ma t= G(τ)+ τ=τ ma + G(τ ma)ep τ τ ma ξ τ τ ma ξ ) ) (A.7) where ξ is determined from (A.3). There are, of course, possible sources of systematic error in the choice of τ ma and the etraction of ξ. But for a PLA derived for a system of dynamical fermions there is an additional ambiguity. Polyakov lines at a separation less than the stringbreaking scale represent a quark-antiquark system joined by a flu tube, and the contributions to G(τ) come from a spectrum of string-like ecitations. Beyond the string-breaking scale, the Polyakov lines represent two bound states, namely a static quark + light antiquark, and a static antiquark + light quark. In this case the contributions to the connected correlator are associated with hadron echange. So in this case the ratio ξ/ξ nd is picking up contributions from quite different regimes. In view of this, we have calculated the ratio ξ/ξ nd at β = 5.65, where the epectation value of the Polyakov line is very small, and the contributions to ξ/ξ nd are coming mainly from the string-like spectrum, at least on the comparatively small lattice volume of 6 4 used in our simulations. We obtain ξ =.96() from the fit (A.3), with data and fit shown in Fig. 7 shown, and using (A.7) with τ ma = 5 we finally obtain ξ ξ nd =.7±.3 (A.8) which is comparable to some of the results for pure SU() lattice gauge theory quoted in ref. []. The fit to a single correlation length (A.3) is epected to be very accurate at large τ, and the influence of higher ecitations would only be evident at smaller τ. Note that in Fig. 7 the fit (A.3) is in fact quite accurate for all τ >, so in this case the deviation of the ratio from unity must be attributed mainly to the data point at τ =. We must stress again that this result is obtained on a relatively small lattice volume, and is subject to the caveats already mentioned. [] M. Fromm, J. Langelage, S. Lottini, and O. Philipsen, JHEP, 4 (), arxiv: [] C. Gattringer, Nucl.Phys. B85, 4 (), arxiv:4.53. [3] G. Bergner, J. Langelage, and O. Philipsen, JHEP, (5), arxiv:55.. [4] C. Wozar, T. Kaestner, A. Wipf, and T. Heinzl, Phys.Rev. D76, 854 (7), arxiv: [5] J. Greensite and K. Langfeld, Phys. Rev. D9, 457 (4), arxiv: J. Greensite and K. Langfeld, Phys.Rev. D88, 7453 (3), arxiv: [6] J. Greensite and R. Höllwieser, Phys. Rev. D94, 454 (6), arxiv: [7] K. Fukushima and T. Hatsuda, Rept. Prog. Phys. 74, 4 (), arxiv: [8] C. Sasaki, Nucl. Phys. A83, 649C (9), arxiv: [9] T. Schäfer and F. Wilczek, Phys. Rev. Lett. 8, 3956 (999), arxiv:hep-ph/ [] M. Caselle and A. Nada, (7), arxiv: [] I. Bender et al., Nucl.Phys.Proc.Suppl. 6, 33 (99).

9 9 T. C. Blum, J. E. Hetrick, and D. Toussaint, Phys.Rev.Lett. 76, 9 (996), arxiv:hep-lat/959. J. Engels, O. Kaczmarek, F. Karsch, and E. Laermann, Nucl.Phys. B558, 37 (999), arxiv:hep-lat/9933. R. De Pietri, A. Feo, E. Seiler, and I.-O. Stamatescu, Phys.Rev. D76, 45 (7), arxiv: [] J. Greensite and K. Splittorff, Phys.Rev. D86, 745 (), arxiv:6.59. J. Greensite, Phys. Rev. D9, 457 (4), arxiv: [3] S. Necco and R. Sommer, Nucl. Phys. B6, 38 (), arxiv:hep-lat/88. [4] B. Bahrampour, B. Wellegehausen, and L. von Smekal, PoS LATTICE6, 7 (6), arxiv:6.85.