arxiv: v2 [hep-ph] 5 Sep 2016

Transcription

1 Thermo-magnetic properties o the strong coupling in the local Nambu Jona-Lasinio model arxiv: v [hep-ph] 5 Sep 016 Alejandro Ayala 1,, C. A. Dominguez, L. A. Hernández, M. Loewe,3,4, Alredo Raya 5, J. C. Rojas 6, C. Villavicencio 7 1 Instituto de Ciencias Nucleares, Universidad Nacional Autónoma de México, Apartado Postal , México Distrito Federal 04510, Mexico. Centre or Theoretical and Mathematical Physics, and Department o Physics, University o Cape Town, Rondebosch 7700, South Arica. 3 Instituto de Física, Pontiicia Universidad Católica de Chile, Casilla 306, Santiago, Chile. 4 Centro Cientíico-Tecnológico de Valparaíso, Casilla 110-V, Valparaíso, Chile. 5 Instituto de Física y Matemáticas, Universidad Michoacana de San Nicolás de Hidalgo, Ediicio C-3, Ciudad Universitaria, Morelia, Michoacán 58040, Mexico. 6 Departamento de Física, Universidad Católica del Norte, Casilla 180, Antoagasta, Chile 7 Departamento de Ciencias Básicas, Facultad de Cienicas, Universidad del Bío-Bío, Casilla 447, Chillán, Chile. We study the thermo-magnetic properties o the strong coupling constant G and quark mass M entering the Nambu-Jona-Lasinio model. For this purpose, we compute the quark condensate and compare it to lattice QCD (LQCD) results to extract the behavior o G and M as unctions o the magnetic ield strength and temperature. We ind that at zero temperature, where the LQCD condensate is ound to monotonically increase with the ield strength, M also increases whereas G remains approximately constant. However, or temperatures above the chiral/deconinement phase transitions, where the LQCD condensate is ound to monotonically decrease with increasing ield, M and G also decrease monotonically. For inite temperatures, below the transition temperature, we ind that both G and M initially grow and then decrease with increasing ield strength. To study possible consequences o the extracted temperature and magnetic ield dependence o G and M, we compute the pressure and compare to LQCD results, inding an excellent qualitative agreement. In particular, we show that the transverse pressure, as a unction o the ield strength, is always negative or temperatures below the transition temperature whereas it starts o being positive and then becomes negative or temperatures above the transition temperature, also in agreement with LQCD results. We also show that or the longitudinal pressure to agree with LQCD calculations, the system should be described as a diamagnet. We argue that the turnover o M and G as unctions o temperature and ield strength is a key element that drives the behavior o the quark condensate going across the transition temperature and provides clues or a better understanding o the inverse magnetic catalysis phenomenon. PACS numbers: t, 1.38.Aw I. INTRODUCTION The study o the properties o strongly interacting matter at high temperature and/or density, under the inluence o magnetic ields, has become a research subject o growing interest over the last years. Examples o physical systems where such properties are relevant include heavy-ion collisions and compact stellar objects. One o the intriguing aspects o these properties, not yet well understood (at least in a consensual way), are the possible causes o the so called inverse magnetic catalysis (IMC) phenomenon ound by lattice QCD (LQCD) calculations [1 3]. Recall that IMC is characterized by a decreasing critical temperature (T c ) or the chiral/deconinement phase transition and a decreasing quark condensate above T c, with increasing ield strength. Dierent approaches have been explored in order to either ind or include IMC in QCD [4], some o these provide an explanation [5], and almost all suggest the need to include extra ingredients in terms o magneticinduced modiications o QCD properties. In particular, the modiication o the QCD coupling due to magnetic screening at low temperature and antiscreening at high temperature has been shown to be a plausible mechanism to explain IMC [6 9]. This also seems to be the reason why eective models without such modiications do not describe neither the behavior o the critical temperature nor the properties o the quark condensate at high temperature [10 17]. For recent reviews see [18, 19]. Deducing the detailed screening/antiscreening properties o the strong coupling as a unction o the ield strength is not a simple task since these properties belong or the largest portion o the parameter space to the non-perturbative domain. Nevertheless, it should be possible to extract general eatures o this coupling by resorting to combining inormation rom eective models and LQCD. The Nambu Jona-Lasinio(NJL) is one o such models. It has been extensively used to explore the chiral transition [0, 1]. In particular, the NJL model can be used to ormulate a simpliied version o the QCD gap equa-

2 tion by means o the Schwinger-Dyson technique where the dynamically generated mass M is constant (momentum independent) and the interaction is given by a ourermion contact term whose strength is controlled by a coupling G. A pertinent question is whether it is possible to extract inormation on the behavior o the coupling G as a unction o the magnetic ield strength and the temperature by combining the NJL model with LQCD data or the quark condensate in the presence o a magnetic ield, and whether this inormation can be used to get clues on the microscopical origins o IMC. In this work we take this approach. We use the NJL model to extract the behavior ogand M as unctions o the magnetic ield or dierent temperatures using LQCD data or the quark condensate []. Notice that i there is a link between the ading-out o the condensate, as a unction o eb, above the critical temperature or the chiral/deconinement transition, and the thermo-magnetic dependent coupling, then the latter should also decrease with the ield intensity. This behavior would signal that a decreasing coupling contributes to a less intense bound between quark-antiquark pairs above the critical temperature, as the ield intensity grows. In this work we show that this is the case. The work is organized as ollows: In Sec II we set up the ramework writing the expression or the gap equation and the quark propagator obtained rom the NJL model in the constant (momentum independent) mass approximation. These equations, together with the expression or the quark condensate given in terms o the quark propagator, provide the set o equations that allow inding the behavior o M and G as unctions o the ield strength, or dierent temperatures. In Sec. III we include the eects o the magnetic ield by means o Schwinger s proper time method. Since the NJL model is not renormalizable, in order to ind the behavior o M and G as unctions o the magnetic ield, we irst separate the vacuum contribution rom the thermo-magnetic one. This procedure has been shown to give reliable results or matter and magnetic ield induced properties in the NJL model []. In Sec. IV we compute the thermomagnetic dependence o M using as input the LQCD behavior o the average light quark condensate as a unction o the ield strength. We also ind the value o G rom the gap equation. Since the light quarks have dierent charges, we ind the values or G and M by averaging over the light-quark lavors. We also compute the thermo-magnetic contribution to the pressure and show that above T c the transverse pressure starts o being positive or small ield strengths. Below T c this pressure starts rom zero and then becomes negative as the ield strength increases, in agreement with LQCD calculations. We argue that this result goes in line with the idea that above (below) T c quarks are brought together (pushed apart) and this makes the coupling become weaker (stronger) due to asymptotic reedom. We also explore the sensitivity o the results or the pressure to variations o the vacuum parameters and ind that there is no signiicant dependence. We inally summarize and conclude in Sec. V. II. GAP EQUATION AND QUARK CONDENSATE IN THE NJL MODEL The NJL model is deined by means o the Lagrangian density L = ψ(i m)ψ +G( ( ψψ ) + ( ψiγ 5 τψ ) ), (1) where τ are the Pauli matrices in isospin space, and ψ is a quark ield. On general grounds, in the mean ield approximation, and ater a bosonization process, the Lagrangian can be rewritten as a vacuum term plus a ree ermion Lagrangian with a dressed mass, namely L MF = σ 4G + ψ(i M)ψ, () where σ = 4G ψψ and M = m+σ. Here we do not consider pion condensation eects, so the only contribution comes rom the sigma meson in the bosonization procedure. The value o the mean ield is determined through the gap equation, obtained by minimizing the eective potential with respect to the mean ield [0, 1] M m = 4G d 4 p Tr[iS(p)], (3) (π) 4 with the trace reerring to color and Lorentz indices. We notice that the quark condensate ψψ is given by d 4 p ψψ = Tr[iS(p)]. (4) (π) 4 In the absence o thermo-magnetic eects, the propagator is given by S(p) = p +M 0 p M0 (5) +iǫ. We now proceed to include magnetic ield and temperature eects in the model. III. THERMO-MAGNETIC EFFECTS Eqs. (3) and (4) represent the two independent equations providing inormation on the thermo-magnetic behavior o the coupling G and the dynamically generated mass M, ater using LQCD results or the quark condensates []. To account or the magnetic ield, we emphasize that the above described bosonization does not aect the orm o the gap equation (3) nor the condensate (4), and the eect o the magnetic ield is relected in the dressing o

3 3 the quark propagator. For this we resort to Schwinger proper time representation o the two-point unction ds is(p) = 0 cos(q Bs) eis(p p [ (cos(q Bs)+γ 1 γ sin(q Bs))(M + p ) tan(q Bs) q Bs M +iǫ) ] p, cos(q Bs) (6) where q is the absolute value o the quark charge (i.e. q u = e /3 and q d = e /3), and we have chosen the homogeneous magnetic ield to point in the ẑ direction, namely B = Bẑ. This coniguration can be obtained rom an external vector potential which we choose in the so called symmetric gauge We have also deined A µ = B (0, y,x,0). (7) p µ (0,p 1,p,0), p µ (p 0,0,0,p 3 ), p p 1 +p, p p 0 p 3. (8) Notice that since the magnetic ield breaks Lorentz invariance, the propagator involves a non-local, albeit path independent phase. However, by taking the trace over a closed one-loop diagram, as is required or the calculation o the condensate, this phase does not contribute and thus we ignore it in the sequel. Using Eq. (6) to take the trace in Eq. (4), we obtain ψψ = 4N c M 1 0 d 4 p (π) 4 dse is(p p tan(q Bs) q Bs M +iǫ), (9) where in order to account or the dierent quark charges q we have averaged over quark-lavors. The integration over the transverse momentum components can be carried out, leading to d p (π) e i tan(q Bs) q B p = q B 4πi 1 tan(q Bs). (10) In order to introduce a inite temperature, within the Matsubara ormalism, we transorm the integrals to Euclidean space by means o d p (π) it + dp3 (π), (11) where the integral over the zeroth component o the ermion momentum has been discretized. We also perorm the change o variable s = iτ. Thereore, the expression or the quark condensate in Eq. (9) becomes ψψ = N c M 1 + q B π T + + dp 3 (π) dτ tanh(q Bτ) e τ( ω n +ω), (1) where we have introduced the ermion Matsubara requencies ω n = (n+1)πt and ω = p 3 +M. Since the NJL model is non-renormalizable the integral above needs to be regularized. This can be done in dierent ways, the simplest being the introduction o an ultraviolet cut-o. This is tantamount to introducing a regulator as the lower limit cuto in the proper-time representation. The parameter represents such regulator. The sum in Eq. (1) can be expressed in terms o Jacobi s ϑ 3 (z,x) unction, deined as ϑ 3 (z,x) = + exp(iπxn +iπzn), (13) whereby + e τ(n+1) π T = e π τt ϑ 3 (πiτt,4πiτt ). (14) For our purposes it is useul to invoke the inversion ormula ( i z ) ( ) ϑ 3 (z,x) = x exp π z ϑ 3 ix x, 1, (15) x which leads to ψψ = N cm 1 4π Using ϑ 3 ( 1, ϑ 3 ( 1, i 4τπT ) i 4πτT = + q B dτ τ tanh(q Bτ) e τm ). (16) + = 1+ ( 1) n exp ( n 4τT n=1 ( 1) n exp ) ( n 4τT ),(17) wenoticethatthet = 0termcorrespondston = 0inthe above expression. Thereore, the vacuum contribution is obtained rom the n = 0 term in the limit where q B 0. Adding and subtracting 1 in the integrand we get ψψ = N cm 4π { + 1 q B + q B dτ τ e τm dτ τ e τm ( 1) n n=1 ] q Bτ tanh(q Bτ) 1 [ e dτ e τm n } 4τT, (18) τ tanh(q Bτ)

4 4 ψψ 1/3 0 M 0 G 0 m Tc NJL (GeV) (GeV) (GeV) (GeV) (GeV) (GeV) TABLE I. Two sets o values or the vacuum regulator, condensate ψψ 0 and dynamically generated mass M 0 stemming rom requiring that the pion mass and the pion decay constant computed in the NJL model attain their physical values. Shown also are the corresponding vacuum values or the coupling constant G 0, current quark mass m and the critical temperature or eb = 0. where we can identiy the vacuum condensate as given by the expression ψψ 0 = N cm 0 4π 0 dτ τ e τm 0, (19) whereas the thermo-magnetic contribution is given by ψψ B,T = N cm 1 { 4π q B [ ] dτ q Bτ τ e τm tanh(q Bτ) 1 + ( 1) n n=1 0 e n } dτ e τm 4τT. (0) τ tanh(q Bτ) ThequantityM M(B,T)inEq.(0)issuchthatwhen B,T 0, M M 0. It turns out that the integrals in Eq. (0) are inite as the lower limit o integration goes to zero. This means that the thermo-magnetic eects are independent otheregulator andwehaveconsequently set this lower limit to zero in Eq. (0). Also, the irst integral in this equation can be computed analytically, with the result [ ] { dτ q Bτ M 0 τ e τm tanh(q Bτ) 1 = q B 1 ( ) ( ) ( ) M M M ln ln q B q B q B ]} + lnγ [1+ M. (1) q B Equation (19) can be used to ix the vacuum values o the quark condensate and the dynamically generated mass M 0 rom a choice o. Two sets o consistent choices, that reproduce the physical values o the pion mass and o the pion decay constant [] are shown in Table 1. Shown also are the corresponding vacuum values or the coupling constant G 0 and the current quark mass m. Beore proceeding to use Eqs. (0) and (1) to ind the thermo-magnetic behavior o the dynamically generated mass M and the coupling G and the consequences or the pressure, it is important to notice that the critical temperature given by the model depends on the choice om 0 and doesnotcoincidewith the correspondingvalue reported by lattice. It is thereore necessary to scale the values o the model temperatures to make them correspond to the physical values. The simplest choice is a linear scaling such that ( ) T T NJL NJL c = T, () T c where T represents the physical value o the temperature and T c and Tc NJL are the physical and model critical temperatures, respectively. We use as the condition to ind the critical temperature the vanishing o the derivative o the thermal piece o the quark condensate with respect to the temperature, namely where ψψ 0,T = N cm π d dt ψψ 0,T = 0, (3) 0 dp p 1.(4) p +M e p +M /T +1 A straightorward calculation in the approximation where Tc NJL and M are o the same order gives T NJL c.38m. (5) To have a better estimate o Tc NJL, one needs to evaluate the above equation at an appropriatevalue o M. We observe that at the critical temperature, the dynamically generated mass drops rom its vacuum value to about hal o it, namely to M 0 /. Using this as the working criterium and the reported LQCD value or the critical temperature T c = GeV, we obtain the corresponding values or the model critical temperature, which we also show in Table 1. IV. THERMO-MAGNETIC BEHAVIOR OF THE DYNAMICALLY GENERATED MASSES AND COUPLINGS To establish how well the solutions or M(B,T) describe the condensates as unctions o eb and T one considers Fig. 1. The igure shows the LQCD average quark condensate (Σ u + Σ d )/ compared to the equivalent quantity 1+ ψψ B,T / ψψ 0 computed within the model, using Eq. (0) or ψψ B,T and one o the values

5 5 FIG. 1. Comparison between the lattice QCD results rom Re. [] or the average quark condensate and the model calculation, as a unction o the magnetic ield. Curves (a), (b), (c) and (d) correspond to T = 0, 113, 130 and 176MeV, respectively. The model describes better the lattice results or lower temperatures. in Table 1 or ψψ 0. Equation (0) is a transcendental equation with none, one or multiple solutions or M, depending on the values ot and B. The procedure we ollowto ind the reported FIG.. Thermo-magnetic average quark mass M(B, T) as a unction o the ield strength eb or the temperatures T = 0MeV, 113MeV, 130MeV and 176MeV computed using the irst set o values in Table 1. FIG. 3. Thermo-magnetic average coupling G(B, T) as a unction o the ield strength eb or the temperatures T = 0MeV, 113MeV, 130MeV and 176MeV computed using the irst set o values in Table 1.

6 6 value om is to averagethe multiple solutionsin the case that there is more than one or to deine as the solution the value o M that provides the closest distance between the lattice value and the model. From Fig. 1 we notice that our description o the LQCDresults is better or low temperatures. Figure 1 has been prepared using the irst set o values in Table 1. The behavior o the thermo-magnetic average mass M(B,T) and coupling G(B,T) as unctions o the ield strength are depicted in Figs. and 3, respectively. Notice that or T = 0, the mass increases monotonically with the magnetic ield. However, there is a turn-over behavior or intermediate values o T where, as unctions o eb the masses start o increasing to then decrease as the ield strength increases. For the largest temperature, which is above the transition temperature, the mass becomes a monotonically decreasing unction o the ield strength. A similar behavior is observed or the coupling. Although G(B,0) is practically constant, or the temperature above the transition temperature the coupling becomes a monotonically decreasing unction o the ield strength. For intermediate temperatures there is also a a turn-over behavior where as unctions o eb the couplings start o increasing and then decrease as the ield strength increases. To test the sensitivity o the results to a change in the vacuum parameters, Figs. 4 and 5 show the behavior o the M(B,T) and G(B,T) as unctions o the ield strength or dierent temperatures, when using the second set o values in Table 1 or the calculation. Notice that the results are qualitative and quantitatively similar to the ones obtained rom the irst set o values in Table 1. In order to study one o the consequences o the behavior o the mass and coupling, we proceed to compute the thermo-magnetic contribution to the pressure. Notice that the magnetic ield induces a dierence between the pressure in the directions parallel and perpendicular to the ield; a magnetization in the ormer direction is absent, while in the latter it contributes. We call the irst kind o pressure longitudinal, that is, directed along the ẑ axis, whereas we call the second kind o pressure transverse, that is, directed along the ˆx, ŷ directions. We consider only the renormalized contribution to the pressure in the so called Φ-scheme [3]. In the mean ield approximation, the longitudinal contribution to the pressure can be written as P z = V e, (6) where V e, is the eective potential. Thereore, using Eqs. (3) and (4), P z can be written as P z = (M(B,T) m) 4G i d 4 p Tr (π) 4 ln(is 1 ), (7) whereas the magnetization M is given by M = Ve (8) (eb)ẑ, rom where the transverse pressure can be computed as [3] P x,y = P z +e B M. (9) To compute M, we observe that M(B,T) has a mild dependence on eb. Thereore, we only consider the terms coming rom the explicit dependence o eb o the eective potential and o G(B,T), which are by ar the dominant contributions. Notice that or the computation o the pressure we use M(B,T) and G(B,T), namely, the averagemass and FIG. 4. Thermo-magnetic average quark mass M(B, T) as a unction o the ield strength eb or the temperatures T = 0MeV, 113MeV, 130MeV and 176MeV computed using the second set o values in Table 1. FIG. 5. Thermo-magnetic average coupling G(B, T) as a unction o the ield strength eb or the temperatures T = 0MeV, 113MeV, 130MeV and 176MeV computed using the second set o values in Table 1.

7 7 FIG. 6. Longitudinal and transverse pressures as unctions o the ield strength eb or T = 113MeV, computed using the irst set o vacuum values in Table 1. The longitudinal (transverse) pressure is a monotonically increasing (decreasing) unction that or this temperature starts o rom zero and grows (decreases) towards positive (negative) values as the ield strength increases. FIG. 7. Longitudinal and transverse pressures as unctions o the ield strength eb or T = 176MeV, computed using the irst set o vacuum values in Table 1. The longitudinal (transverse) pressure is a monotonically increasing (decreasing) unction that or this temperature starts o rom positive values and grows (decreases) towards positive (negative) values as the ield strength increases. coupling, respectively. Thereore, the pressure and magnetization are correspondingly also computed as an average over the light lavors. Figures 6 and 7 show the longitudinal and transverse pressures as unctions o the ield strength, or T = 113 MeVandT = 176MeV,respectively,computedusingthe irst set ovacuum parametersin Table 1. Notice that or T = 113 MeV, that is, or a temperature below T c, these pressures start o rom zero and have opposite behaviors; the longitudinal pressure is a monotonically increasing unction towards positive values whereas the transverse pressure is a monotonically decreasing unction towards negative values. For the case o T = 176 MeV, that is or a temperature above T c, both pressures start o rom positive values. However, there is an interval o ield strengths where the transverse pressure is positive to then change sign and become negative. This results are in agreement with the indings o Re. [3]. To test the sensitivity o the results to the vacuum parameters, Figs. 8 and 9 show the same pressures computed using the second set o parameters in Table 1. The results are equivalent. It is important to notice that the calculation describing the longitudinal pressure in Fig. 6 Fig. 9 agree with the LQCD indings [3], provided that e B M = ebm, (30) where M represents the magnitude o the magnetization vector. This means, that the magnetization is overall opposite to the magnetic ield which in turn means that the system is well described in the model as possessing diamagnetic properties, both, below and above the critical temperature. Finally, in order to study the pressure behavior reerred to purely thermal eects, Figs. 10 and 11 show P P z (B,T) P z (0,T), computed or the irst and second set o vacuum values in Table 1, respectively. In all cases P is well described by a monotonically and positive deinite increasing unction o eb. This behavior is also in agreement with LQCD calculations [3]. We observe that the rate o change shows a turn over behavior as the temperature increases. For T = 0 the rate o increase is small, becoming aster or intermediate temperatures to then decrease or the highest temperature. V. SUMMARY AND CONCLUSIONS In this paper we studied the thermo-magnetic behavior o the coupling constant and mass in the NJL model. We used the gap equation and the quark condensate expressions obtained rom the model together with LQCD results or the light-quark condensates in the presence o a magnetic ield []. Although we obtained the behavior o the couplings G(B,T) and M(B,T) as a unction o the temperature and the magnetic ield strength, we did not attempt a detailed magnetic ield dependence study o a list o observables, neither to look or a parametrization o the coupling as a unction o the temperature and/or the ield strength. A study o this sort has been recently carried out in Re. [3]. Our results show that or temperatures above the transition temperature, the couplings are monotonically decreasing unctions o the ield strength. This means that at these temperatures the melting o the quark condensates is accompanied by a corresponding decrease in the strength o the interaction that binds these quarks. For temperatures close to, but below the transition

8 8 FIG. 8. Longitudinal and transverse pressures as unctions o the ield strength eb or T = 113MeV, computed using the second set o vacuum values in Table 1. The results are equivalent to the ones obtained using the irst set o vacuum values in Table 1. FIG. 9. Longitudinal and transverse pressures as unctions o the ield strength eb or T = 176MeV, computed using the second set o vacuum values in Table 1. The results are equivalent to the ones obtained using the irst set o vacuum values in Table 1. temperature, we ind a turn over behavior o the couplings. As the ield strength starts increasing, the couplings increase. However, or intermediate values o the ield strength the couplings decrease. This signals that as the temperature decreases below, but close to the transition temperature, the strength o the coupling increases. This increase is accompanied by a corresponding increase in the value o the condensate, as shown by LQCD calculations. Nevertheless, this increase is not sustained, since or stronger ields the couplings decrease, as do the LQCD computed condensates. The results or G(B,T) strengthen the picture advocated in Res. [8] where the behavior o the condensate as a unction o the magnetic ield is directly linked to the properties o the strong coupling constant at high and low temperatures. FIG. 10. P P z(b,t) P z(0,t), computed or the irst set o vacuum values in Table 1 or T = 0, 113, 130 and 176 MeV. We also computed the thermo-magnetic contribution to the longitudinal and transverse pressures. We ound thatbelowt c, the transversepressureasaunction othe magnetic ield, decreases towards negative values starting o rom zero. However, or temperatures above the transition temperature, although the transverse pressure still decreases as a unction o the ield strength, it starts o rom positive values. This turnover behavior o the transverse pressure means that above T c particles are pulled closer together, at least or small values o the magnetic ield. The act that at the same time the coupling decreases can be viewed as signaling that the strength o the bound o the condensate is smaller, due to asymptotic reedom and this can be responsible or the decrease o the condensate as the magnetic ield strength is turned on. Last but not least, we ound that in order or the computation o the longitudinal pressure to agree with LQCD calculations, the system should be described as a diamagnet at inite temperature, that is, with its magnetization opposite to the magnetic ield direction, both above and below the critical temperature. Overall, the results suggest that IMC, as described by the thermo-magnetic behavior o the quark condensate, can be linked to the properties o the coupling constant as a unction o the magnetic ield in a wide range o temperatures. ACKNOWLEDGMENTS A.A. and L.A.H. acknowledge valuable discussions with N. Scoccola. This work has been supported in part by UNAM-DGAPA-PAPIIT grant number IN101515, by Consejo Nacional de Ciencia y Tecnología grant number 56494, by CIC-UMSNH (México) grant number 4., National Research Foundation (South Arica), the Harry Oppenheimer Memorial Trust OMT Re. 04/0 and by FONDECYT (Chile) grant numbers , , and CV acknowledges sup-

9 port rom the group Física de Altas Energías at UBB. L.A.H. acknowledges the University o Cape Town and the National Research Foundation (South Arica) or unding assistance. M.L acknowledges support rom Proyecto Basal (Chile) FB FIG. 11. P P z(b,t) P z(0,t), computed or the second set o vacuum values in Table 1 or T = 0, 113, 130 and 176 MeV. [1] G. Bali, F. Bruckmann, G. Endrodi, Z. Fodor, S. Katz,et al., JHEP 10, 044 (01). [] G. Bali, F. Bruckmann, G. Endrodi, Z. Fodor, S. Katz, et al., Phys. Rev. D 86, (01). [3] G. Bali, F. Bruckmann, G. Endrodi, S. Katz, and A. Shaer, JHEP 1408, 177 (014). [4] E. S. Fraga, J. Noronha and L. F. Palhares, Phys. Rev. D 87, (013); F. Bruckmann, G. Endrodi and T. G. Kovacs, J. High Energy Phys. 1304, 11 (013); G. S. Bali, F. Bruckmann, G. Endrodi, F. Gruber, A. Schaeer, J. High Energy Phys. 1304, 130 (013); E. S. Fraga, B. W. Mintz, J. Schaner-Bielich, Phys. Lett. B 731, 154 (014); G. Endrödi, J. High Energy Phys. 1507, 173 (015); J. O. Andersen, W. R. Naylor, A. Tranberg, J. High Energy Phys. 150, 04 (015); A. Ayala, C. A. Dominguez, L. A. Hernández, M. Loewe, J. C. Rojas, C. Villavicencio, Phys. Rev. D 9, (015); R. Rougemont, R. Critelli, J. Noronha, arxiv: [hep-th]. [5] K. Fukushima and Y. Hidaka, Phys. Rev. Lett. 110, (013); J. Chao, P. Chu, M. Huang, Phys.Rev. D 88, (013); M. Ferreira, P. Costa, D. P. Menezes, C. Providencia and N. N. Scoccola, Phys. Rev. D 89, (014); Sh. Fayazbakhsh and N. Sadooghi, Phys. Rev. D 90, (014); N. Mueller and J. M. Pawlowski, Phys. Rev. D 91, 11, (015). [6] A. Ayala, M. Loewe, A. J. Mizher, R. Zamora, Phys. Rev. D 90, (014); J. Braun, W. A. Mian, S. Rechenberger, arxiv: [hep-ph]; A. Ayala, M. Loewe, R. Zamora, Phys. Rev. D 91, (015); A. Ayala, C. A. Dominguez, L. A. Hernández, M. Loewe, R. Zamora, Phys. Rev. D 9, (015); [7] R. L. S. Farias, K. P. Gomes, G. Krein and M. B. Pinto, Phys. Rev. C 90, 0503 (014); M. Ferreira, P. Costa, O. Lourenço, T. Frederico, C. Providência, Phys. Rev. D 89, (014). [8] A. Ayala, J. J. Cobos-Martínez, M. Loewe, M. E. Tejeda- Yeomans, R. Zamora, Phys. Rev. D 91, (015); A. Ayala, C. A. Dominguez, L. A. Hernandez, M. Loewe, R. Zamora, Phys. Lett. B 759, (016). [9] E. J. Ferrer, V. de la Incera, X. J. Wen, Phys. Rev. D 91, (015); A. Ahmad and A. Raya, J. Phys. G 43, (016). [10] J. K. Boomsma and D. Boer, Phys. Rev. D 81, (010). [11] M. Loewe, C. Villavicencio, R. Zamora, Phys. Rev. D 89, (014). [1] N. Agasian and S. Fedorov, Phys. Lett. B663, 445 (008). [13] E. S. Fraga and A. J. Mizher, Phys. Rev. D 78, (008). [14] A. J. Mizher, M. N. Chernodub, and E. S. Fraga, Phys. Rev. D 8, (010). [15] J. O. Andersen and R. Khan, Phys. Rev. D 85, (01). [16] D. Blaschke, S. Fredriksson, H. Grigorian, A. M. Oztas, and F. Sandin, Phys. Rev. D 7, (005). [17] I. General, D. Gomez-Dumm, and N. Scoccola, Phys. Lett. B506, 67 (001). [18] J. O. Andersen, W. R. Naylor, A. Tranberg, Rev. Mod. Phys. 88, (016). [19] V. A. Miransky, I. A. Shovkovy, Phys. Rept. 576, 1-09 (015). [0] M. Buballa, Phys. Rept. 407, 015 (005). [1] S. Klevansky, Rev. Mod. Phys. 64, 649 (199). [] N. Scoccola, private communication. [3] R.L.S.Farias, V.S.Timoteo, S.S.Avancini, M.B.Pinto, G. Krein, arxiv: [hep-ph].