## Brazilian Journal of Physics

*On-line version* ISSN 1678-4448

### Braz. J. Phys. vol.37 no.1a São Paulo Mar. 2007

#### http://dx.doi.org/10.1590/S0103-97332007000100010

**Cold and dense perturbative Yukawa theory with massive fermions **

**L. F. Palhares; E. S. Fraga **

Instituto de Física, Universidade Federal do Rio de Janeiro C.P. 68528, Rio de Janeiro, RJ 21941-972, Brazil

**ABSTRACT**

Recent results for the two-loop thermodynamic potential of QCD at finite density have shown that nonzero quark mass corrections to the pressure are relevant and can dramatically affect the structure of compact stars. Motivated by these findings, we consider a simple toy model - cold and dense Yukawa theory - to study the effects of finite fermion masses on the pressure. The role of renormalization group running of the coupling and mass is also discussed. Results within this simple model might be useful in the description of condensates in the core of neutron stars.

**Keywords:** Finite-temperature field theory; Yukawa theory; Equation of state

The role of finite quark masses in QCD thermodynamics has received increasing attention in the last few years. Since their effects might be relevant in the description of phenomena associated with the chiral and the deconfinement transitions, more realistic investigations of the phase structure of strong interactions should incorporate them.

For finite temperature and vanishing chemical potential, lattice simulations can now investigate the effect of nonzero quark masses more reliably, and preliminary results for the corrections to the behavior of the quark condensate and the renormalized Polyakov loop with temperature have been obtained [1]. Although none of these quantities is a true order parameter under such conditions, they still provide useful information on chiral symmetry breakdown and confinement within the framework of effective field theories [2].

In the case of cold and dense QCD, it was generally believed that the effects of nonzero quark masses on the equation of state were negligible, thereby yielding only minor corrections to the mass-radius diagram of compact stars. At the time of the first investigations of (strange) quark star structure [3], the seminal computations of the thermodynamic potential for cold QCD at finite chemical potential in perturbation theory to ~ had already been done [4-6]. Nevertheless, properties of quark stars were then computed making use of the bag model description of the equation of state with corrections only up to order a_{s} from perturbative QCD, a practice that remains in most treatments to date [7]. In the massless case, however, first-order corrections cancel out in the equation of state, so that one ends up with a free gas of quarks modified only by a bag constant. Including a finite quark mass for the strange quark, for instance, brings a correction that was first estimated to modify the equation of state by less than 5% and was essentially ignored.

In fact, mass, as well as color superconductivity gap, contributions to the pressure are supressed by two powers of the chemical potential as compared to zero-mass interacting quark gas terms. Therefore, assuming a critical chemical potential for the chiral transition of the order of a few hundred MeV, naively those terms should not matter. However, recent results for the thermodynamic potential to two loops have shown that corrections are sizable, and may dramatically affect the structure of compact stars [8]. Moreover, the situation in which mass (as well as gap) effects are significant corresponds to the critical region for chiral symmetry breakdown in the phase diagram of QCD. Hence, not only the value of the critical chemical potential will be affected, but also the nature of the chiral transition. In particular, if the latter is strongly first-order there might be a new class of compact stars, smaller and denser, with a deconfined quark matter core [9]. Of course, contributions due to color superconductivity [10] as well as chiral condensation [11] will also affect this picture.

In this paper, we consider a simple toy model - cold and dense Yukawa theory - to study the influence of fermion masses on the pressure. The extremely simplified formal structure allows for a detailed investigation of the features that come about in higher-order contributions to the equation of state, some of which will be preserved in the case of cold and dense full QCD. Here, we present a two-loop calculation of the pressure with massive fermions in the modified minimal subtraction () renormalization scheme, and discuss qualitatively the role to be played by the renormalization group running of the coupling and the fermion mass. Higher-order corrections and a thorough analysis of renormalization group effects will be presented elsewhere [12].

The use of corrections ~ with a modern definition of the running coupling constant to model the non-ideality in the equation of state for cold, dense QCD with massless quarks [9] has shown that the renormalization group flow plays a key role, a feature that is probably even more important in the case of massive quarks [8]. In fact, as will be discussed below, the introduction of a nonzero mass for the fermions brings the first logarithmic contribution to the equation of state to two-loop order, whereas in the massless case such corrections appear only at three loops.

Although our framework is assumedly very simplified, our results could be useful not only to provide insight for in-medium QCD phenomena, but also in many different applications, since Yukawa-type couplings appear in several contexts as will be discussed later.

In what follows, we consider a gas of massive fermions whose interaction is mediated by a real scalar field, f, with an interaction Lagrangian of the Yukawa form:

where *g* is the coupling constant.

Although we are interested here in the limit in which the temperature vanishes, it is convenient to work initially with finite temperature and density. Within the imaginary-time framework [13], the time dimension is compactified and associated with temperature (*T*). Furthermore, imposition of periodicity or anti-periodicity of the fields in the imaginary time reflects the spin-statistics theorem, characterizing, respectively, bosons (*B*) or fermions (*F*). As a consequence, only specific discrete Fourier modes depending on the statistical nature of the field are allowed. Therefore, as is customary in finite-temperature field theory, integrals over the zeroth four-momentum component are mapped into discrete sums over Matsubara frequencies, denoted by = 2*n*p*T* and = (2*n*+1)p*T*, with n integer. To take finite density effects into account, the chemical potential µ is included, representing, in practice, a shift in the zeroth component of the fermionic four-momentum.

From the partition function written in terms of the euclidean action *S _{E}*,

*Z*(

*T*,µ) =

*Tr*exp(-

*S*), one derives the perturbative series for the pressure

_{E}*P*:

where *V* is the volume of the system, *Z*_{0} is the partition function of the free theory and S_{I} represents the interaction action. In our case, Wick's theorem implies that only even powers in the above expansion survive. Therefore, the perturbative method results in a power series of a_{Y} º *g*^{2}/4p [15]. Up to *O*(a_{Y}), the first non-trivial contributions to the pressure are given by the free massive gas term, *P*_{0}, and the ''exchange diagram'', *P*_{1}.

Using the Feynman rules in the imaginary time formalism at finite temperature [13], one can compute the pressure of a gas of interacting massive fermions order-by-order in perturbation theory. The free gas contribution for fermions of mass *m* has the form:

where the trace is performed over the Dirac structure, *P*^{µ} = (+µ,**p**) is the fermionic four-momentum, *p _{f}* = is the Fermi momentum, and we use the following notation for the sum-integrals:

The *O*(a_{Y}) correction is given by the exchange term [13,14]:

where b = 1/*T* and *K*^{µ} = (,**k**) is the bosonic four-momentum.

Resorting to a particular representation of the Dirac delta function [14], one can rewrite the Matsubara sums in (5) in a decoupled form and solve them independently via equivalent contour integrations. We obtain, performing all the sum-integrals explicitly, the following renormalized result in the limit *T* ® 0:

where L is the renormalization scale in the scheme.

It is interesting to notice that the inclusion of mass for the fermions brings the presence of logarithmic corrections one order down in a_{Y}. In massless QCD, for instance, one has ~ a_{s} corrections at two loops and ~ {, ln a_{s}, ln(L/µ)} at three loops [9]. However, in the massive case one finds not only ~ a_{s} terms, but also a contribution ~ *m*^{2}a_{s}ln(L/µ) at two loops, as can be seen above. In our case, this feature will be important when we incorporate renormalization group running of a_{Y} and *m*.

In Fig. 1 we illustrate the effect of modifying the mass on the total pressure to *O*(a_{Y}), *P* = *P*_{0}+*P*_{1}, normalized by the pressure of the free gas, as a function of the fermionic chemical potential. The choice of range for µ, and accordingly for the masses, are inspired by the scales found in the case of QCD [8]. In the same vein, the coupling is fixed to a_{Y} = 0.3. It is clear from the figure that mass corrections bring significant changes to the pressure, even in the absence of renormalization group (RG) running for the coupling and the mass. On the other hand, to perform a RG analysis it is more convenient to plot *P* vs. µ and look at the small µ region [12], given the sign of the beta function for this theory.

Fig. 2 displays the dependence on the renormalization scale L. The values chosen are motivated by the ones which appear in QCD, as before. Although the effects of varying L appear to be relatively small, it would be premature to conclude that this feature will remain after implementing the RG flow. In fact, the results presented in Fig. 2 most probably underestimate the scale dependence of the full correction, since not only the coupling but also the mass will run with L. In the Yukawa theory, in contrast to QCD, the effect will become larger as we increase the chemical potential. For fixed coupling, larger values of L yield larger modifications in the pressure. However, after the inclusion of RG running, this behavior can be mantained, as should be the case here, or become the opposite, as is the case in QCD, depending on the sign of the beta function. Since the L-dependence comes from the term ~ *m*^{2}a_{Y}ln(L/*m*) in (6), there will be a competition between the behavior of the renormalization scale L and that of m and a_{Y} as functions of µ.

Recently, finite fermion mass corrections were shown to affect significantly the cold and dense equation of state for QCD [8], despite the previous disbelief in its relevance. In this work, we investigated further these effects in the context of the Yukawa theory. The first non-trivial correction to the pressure of a free gas of massive fermions at finite density was calculated. Even at two-loop order, as discussed before, mass effects bring into play logarithmic corrections originated in the subtraction scheme. As usual, they bring about a non-physical dependence on the renormalization scale L, since one has to cut the perturbative series at some order. Higher-order computations in this framework are in progress [12], and will give a better handle on the choice of this scale, which in our case should be a function of µ and *m*. On the other hand, one can also choose the scale in a phenomenological way in a given model, imposing physical constraints to the equation of state, as was done in Ref. [9] to model the non-ideality of QCD at finite density with massless quarks.

A detailed renormalization group flux analysis is an essential step towards understanding the role of fermion masses in the equation of state at different energy scales and is also under current development [12]. These results may show that the modifications remain stable through the perturbative expansion and contribute to establish the importance of nonzero fermion mass corrections in high-density strongly interacting media. In the framework of effective theories, in which Yukawa-type couplings appear frequently, corrections due to finite fermion masses might yield appreciable deviations from the conventional results and affect different physical systems. If one also considers temperature effects, several applications are found from the realms of the electroweak phase transition in the early universe to the formation of condensates in the core of neutron stars.

**Acknowledgments **

The authors would like to thank R.D. Pisarski, J. Schaffner-Bielich and C. Villavicencio for fruitful discussions. This work was partially supported by CAPES, CNPq, FAPERJ and FUJB/UFRJ.

[1] C. Bernard *et al.* [MILC Collaboration], Phys. Rev. D **71**, 034504 (2005); [ Links ]Z. Fodor and S. D. Katz, JHEP **0404**, 050 (2004); [ Links ]P. Petreczky and K. Petrov, Phys. Rev. D **70**, 054503 (2004); [ Links ]Y. Aoki, Z. Fodor, S. D. Katz, and K. K. Szabo, JHEP **0601**, 089 (2006); [ Links ]P. de Forcrand and O. Philipsen, hep-lat/0607017. [ Links ]

[2] A. Dumitru, D. Roder, and J. Ruppert, Phys. Rev. D **70**, 074001 (2004); [ Links ]A. Mocsy, J. Phys. G **31**, S1203 (2005); [ Links ]J. Braun, B. Klein, H. J. Pirner, and A. H. Rezaeian, Phys. Rev. D **73**, 074010 (2006). [ Links ]

[3] E. Witten, Phys. Rev. D **30**, 272 (1984); [ Links ]E. Farhi and R. L. Jaffe, Phys. Rev. D **30**, 2379 (1984); [ Links ]C. Alcock, E. Farhi, and A. Olinto, Astrophys. J. **310**, 261 (1986); [ Links ]P. Haensel, J. L. Zdunik, and R. Schaeffer, Astron. Astrophys. **160**, 121 (1986). [ Links ]

[4] B. A. Freedman and L. D. McLerran, Phys. Rev. D **16**, 1130 (1977); * [ Links ]ibid.*, D **16**, 1147 (1977); *ibid.*, **D16**, 1169 (1977); *ibid.*, D **17**, 1109 (1978).

[5] V. Baluni, Phys. Rev. D **17**, 2092 (1978). [ Links ]

[6] T. Toimela, Int. J. Theor. Phys. **24**, 901 (1985) [Erratum-ibid. **26**, 1021 (1987)]. [ Links ]

[7] N. K. Glendenning, *Compact Stars - Nuclear Physics, Particle Physics, and General Relativity* (Springer, New York, 2000); [ Links ]H. Heiselberg and M. Hjorth-Jensen, Phys. Rept. **328**, 237 (2000). [ Links ]

[8] E. S. Fraga and P. Romatschke, Phys. Rev. D **71**, 105014 (2005); [ Links ]E. S. Fraga, Nucl. Phys. A **774**, 819 (2006). [ Links ]

[9] E. S. Fraga, R. D. Pisarski, and J. Schaffner-Bielich, Phys. Rev. D **63**, 121702 (2001); [ Links ]E. S. Fraga, R. D. Pisarski, and J. Schaffner-Bielich, Nucl. Phys. A **702**, 217 (2002). [ Links ]

[10] M. Alford and S. Reddy, Phys. Rev. D **67**, 074024 (2003); [ Links ]M. Alford, M. Braby, M. W. Paris, and S. Reddy, Astrophys. J. **629**, 969 (2005). [ Links ]

[11] M. Buballa and M. Oertel, Nucl. Phys. A **703**, 770 (2002); [ Links ]S. B. Ruster, V. Werth, M. Buballa, I. A. Shovkovy, and D. H. Rischke, Phys. Rev. D **72**, 034004 (2005); [ Links ]D. Blaschke, S. Fredriksson, H. Grigorian, A. M. Oztas, and F. Sandin, Phys. Rev. D **72**, 065020 (2005). [ Links ]

[12] L. F. Palhares and E. S. Fraga, work in progress. [ Links ]

[13] J. I. Kapusta, *Finite-temperature field theory* (Cambridge University Press, 1989). [ Links ]

[14] J. I. Kapusta, Nucl. Phys. B **148**, 461 (1979). [ Links ]

[15] Since we are concerned only with the zero-temperature limit, there are no odd powers of g coming from resummed contributions of the zero Matsubara mode for bosons in the perturbative series.

Received on 29 September, 2006