landau fermi - liquid theory has seen considerable success in describing a wide variety of fermionic many - particle systems , such as liquid helium-3 , electrons in metals , nuclei and nuclear matter @xcite .
it gives an effective description of the low - lying elementary excitations ( quasiparticles ) at low temperatures , i.e. , at a temperature @xmath1 such that the system can be considered degenerate ( @xmath2 , where @xmath3 denotes the chemical potential ) but still in the normal phase above any symmetry breaking phase transition ( @xmath4 ) , for example , to a magnetic or superconducting phase .
nevertheless , the theory is not only important for the description of the above `` normal '' ( fermi - liquid ) phase , but is also vital to correctly describe the emergence of a possible ordered phase @xcite . indeed , according to the bardeen - cooper - schrieffer ( bcs ) theory the onset of superconductivity is to be viewed as an instability of the fermi liquid under an attractive interaction , which results in the formation of cooper pairs @xcite . the same is true for magnetism in the case of repulsive interactions .
in particular , bcs theory implies that the high - density and low - temperature region of the phase diagram of quantum chromodynamics ( qcd ) contains a color superconducting phase of quarks @xcite .
in fact , since quarks carry color , flavor and spin quantum numbers , many distinct superconducting phases are possible and are characterized by the various symmetries of the cooper - pair wavefunction .
these phases of cold and dense qcd might occur in the core of neutron stars where matter is compressed to several times the nuclear density @xmath5 .
the presence of such a superconducting phase is expected to have observable consequences on the cooling and magnetic fields of neutron stars @xcite .
a particularly interesting phase is the so - called 2sc phase in which only two flavors of quarks are paired , while the unpaired quarks of the third flavor remain a fermi liquid .
this phase further illustrates the importance of understanding the quark - gluon plasma , i.e. , the normal phase of quarks .    fortunately , due to asymptotic freedom , the coupling constant @xmath6 of qcd becomes small at large chemical potential @xmath3 , such that a systematic study of the high - density and low - temperature region of the qcd phase diagram is possible using perturbation theory . at high densities
the dominant interaction between quarks is that of one - gluon exchange , where the long - range behavior of the gluons is screened due to the quark - gluon plasma .
while the electric gluons ( longitudinal ) are screened by a debye mass @xmath7 , the magnetic ( transverse ) gluons are only dynamically screened @xcite .
this residual long - range behavior of the magnetic gluons dominates the low - temperature behavior of the system .
examples of the effects of magnetic gluons can be seen in the non - standard scaling of the critical temperature of the superconducting phase @xmath8 @xcite , which depends exponentially only on @xmath9 instead of showing the expected @xmath10 behavior , or of the life - time of elementary excitations of ( grand - canonical ) energy @xmath11 near the fermi surface @xmath12 @xcite that signal a ( marginal ) non - fermi - liquid behavior at zero temperature .    in the context of a possible deconfined phase of quarks inside neutron stars ,
the fermi system is expected to be imbalanced due to the different chemical potentials of the various quarks .
this deconfined phase may contain @xmath13 , @xmath14 , and @xmath15 quarks with different densities as a consequence of their different masses and charges . in cold atomic physics
, an analogous population imbalance has been realized experimentally between two spin states @xcite , which resulted in frustration of the pairing between particles .
this will cause very different behavior from the normal bcs case and is currently a very hot topic in a wide variety of fields @xcite . in particular , the imbalance between different flavors of quarks is expected to significantly alter the properties of the qcd phase diagram , such as the 2sc phase mentioned above or the color - flavor locked phase where all the colors and flavors are paired @xcite . in the present paper ,
we therefore generalize the quasiparticle properties based on fermi - liquid theory to an imbalanced system .    to do so , the effect of the dressed gluons on the thermodynamic potential , specific heat and effective mass of the quasiparticles is determined at low temperatures .
the thermodynamic potential will be calculated for a two - flavor quark system in the limit of a large number of quark flavors ( large-@xmath0 ) .
this approximation corresponds to the random phase approximation ( rpa ) , which has been quite successful in condensed - matter systems such as the interacting electron gas in metals @xcite . using the framework of landau fermi - liquid theory for relativistic systems , derived in section [ sec - theory ] ,
the effective mass of the quasiparticles can be determined from the thermodynamic potential by considering the specific heat .
the result and its implications on the applicability of fermi - liquid theory is discussed .
as we will show , the logarithmic dependence of the temperature , which is due to the transverse gluons , shows the breakdown of fermi - liquid theory at low temperatures .
furthermore , even though within our approximations there is no interaction between the different flavors of quarks in the limit of weak coupling , the rpa correction still presents some mixing between the different flavors , which stems from the long - wavelength screening of the longitudinal gluons . above and throughout the article natural units
are used , i.e. , units such that @xmath16 . other conventions and technical details
of the calculation can be found in the appendices .
to describe the normal state of the quark - gluon plasma at nonzero temperatures , which is a strongly interacting gas of quarks , landau fermi - liquid theory can be used .
this theory takes as a starting point the non - interacting fermi gas and switches on the interaction adiabatically .
as long as the temperature of the system is much higher than the critical temperature @xmath17 for superconductivity , no bound states ( cooper pairs ) will form and each state of the non - interacting fermi gas is transformed into a state of the interacting gas .
therefore , the excitations of such a fermi liquid remain of a fermionic nature . for simplicity
, fermi - liquid theory will be introduced here for a system at zero temperature even though the quark - gluon plasma is formally a marginal fermi liquid in that case .
however , the introduced concepts turn out to be valid for nonzero temperatures much lower than the fermi temperature @xmath18 , which is expected to be around @xmath19 for quarks in the core of neutron stars .
the basis of landau fermi - liquid theory is to consider the effect of a small change of the ground - state fermi distribution on the thermodynamic potential density , which can be written as @xmath20 here @xmath21 is a small change in the ground - state momentum distribution of species @xmath22 , @xmath23 is the particle density , @xmath24 is the energy of a quasiparticle that for our purposes only depends on the magnitude of the vector @xmath25 , and @xmath26 is the chemical potential . in the case of electrons
the index @xmath22 specifies the spin of the electron , while for the case of quarks it specifies the spin , color and flavor of the quark .
hence , the quasiparticle energy and the effective interaction between quasiparticles are defined as @xmath27    an important quantity in fermi - liquid theory is the effective mass @xmath28 of a quasiparticle of species @xmath22 near its fermi surface .
it is defined in terms of the group velocity @xmath29 evaluated at the fermi momentum @xmath30 @xmath31_{k = k_\sigma}^{-1},\ ] ] such that upon linearizing the quasiparticle energy around the fermi surface we obtain @xmath32 with @xmath33 as the fermi energy of species @xmath22 . in the non - interacting limit ,
the effective mass reduces to @xmath34 and @xmath35 for non - relativistic and relativistic dispersions , respectively .    using the fact that the pressure is an invariant under galilean or lorentz transformations ( see appendix .  [ app - invariance ] ) , of which
the effect can be written as a change in the distribution as shown in fig .
[ fig : boosted_fermi_sphere ] , eq .
( [ thermpotential ] ) should give zero for this particular choice of @xmath36 and hence relate @xmath24 to @xmath37 .
in other words , by performing an infinitesimal galilean or lorentz transformation on the fermi sphere , the effective mass can be related to the effective interaction at the fermi surface @xcite .    ]    to this end , consider an infinitesimal change in the fermi surface from @xmath30 to @xmath38 @xmath39 where the change in the fermi surface solely depends on the direction in momentum space denoted by @xmath40 .
the corresponding change in the thermodynamic potential density can be written as    @xmath41+\frac{1}{2}\sum_{\sigma,\sigma'}\iint\frac{\mathrm{d}\hat{\mathbf{k}}}{(2\pi)^3}\frac{\mathrm{d}\hat{\mathbf{k}}'}{(2\pi)^3}\intop_{k_\sigma}^{k_\sigma+\delta k_\sigma}\mathrm{d}kk^2\intop_{k_{\sigma'}}^{k_{\sigma'}+\delta
k_{\sigma'}}\mathrm{d}k'k'^2f_{\sigma\sigma'}(\mathbf{k},\mathbf{k}')-\sum_\sigma n_\sigma\delta\mu_\sigma\nonumber\\ = & \sum_{\sigma}\frac{k_\sigma^3}{4\pi^2m_\sigma^*}\int\frac{\mathrm{d}\hat{\mathbf{k}}}{4\pi}\delta k_\sigma^2(\hat{\mathbf{k}})+\sum_{\sigma,\sigma'}\frac{k_\sigma^2k_{\sigma'}^2}{2(2\pi^2)^2}\iint\frac{\mathrm{d}\hat{\mathbf{k}}}{4\pi}\frac{\mathrm{d}\hat{\mathbf{k}}'}{4\pi}f_{\sigma\sigma'}(k_\sigma , k_{\sigma'},\theta)\delta k_\sigma(\hat{\mathbf{k}})\delta k_{\sigma'}(\hat{\mathbf{k}}')-\sum_\sigma n_\sigma\delta\mu_\sigma.\end{aligned}\ ] ]    in the second line the integrals were expanded to second order in @xmath42 , and @xmath43 was defined as the effective interaction at the fermi surfaces between two species , with @xmath44 the angle between the directions of @xmath25 and @xmath45 .
next we expand @xmath46 in terms of spherical harmonics and @xmath47 in terms of legendre polynomials as @xmath48 for an imbalanced system , the difference in density @xmath23 results in a difference in @xmath26 , which can cause @xmath49 to be different from each other even if their non - interacting bare masses @xmath50 are the same .
for each species , the particle density obeys @xmath51 . inserting these expansions and using the orthogonality of the spherical harmonics , eq .
( [ thermpotexpand ] ) becomes @xmath52 in the balanced case , all fermi momenta are the same , @xmath53 , @xmath54 , then eq .  ( [ reducedthermpotexpansion ] ) reduces to @xmath55,\ ] ] where the intensive quantities @xmath56 , @xmath3 , @xmath57 and @xmath58 are the same for all species , and @xmath59 is the total density , i.e. , the density summed over all species .
the ( dimensionless ) landau fermi - liquid parameters @xmath58 are defined as @xmath60_{k_f}^{-1}=\sum_{\sigma'}\frac{k_fm^*f_{l,\sigma\sigma'}(k_f)}{2\pi^2}.\ ] ] note that @xmath58 is species independent , but the inter - species and intra - species interactions can be different , so the subscript @xmath22 is shown explicitly in @xmath61 but not in @xmath58 . from the above it is possible to derive the equation for the effective mass , as we first show for the simpler balanced case .    consider a change in the distribution due to a lorentz transformation as shown in fig .
[ fig : boosted_fermi_sphere ] .
the shape of the fermi sphere is determined from the equation @xmath62 , where @xmath63 is the four - momentum , @xmath64 is the four - velocity of the fermi sphere with the lorentz factor @xmath65 .
the minus sign in the previous expression and in @xmath66 are a consequence of our choice of the metric @xmath67 . in the rest frame , where @xmath68 , the above reduces to the condition @xmath69 .
similarly , the shape of the fermi surface can be determined from the lorentz transformation of the momentum to a frame moving with velocity @xmath70 , @xmath71 which for an infinitesimal lorentz transformation reduces to @xmath72 and in the non - relativistic limit to @xmath73 .
note that the deformation of the fermi sphere due to the lorentz contraction shows up only in higher - order terms of @xmath74 .
first , consider the non - relativistic case in which an infinitesimal galilean transformation shifts all momenta from @xmath25 to @xmath75 .
this results in a change of the fermi surface according to @xmath76 , such that @xmath77 .
furthermore , the total energy of the system transforms as @xmath78 , where @xmath79 is the total mass of the system , which induces a change in the chemical potential @xmath80 . inserting this into eq .
( [ landaufermithermpot ] ) gives the well - known landau effective - mass relation in the balanced case @xcite @xmath81    secondly , for a relativistic system the momenta are shifted by @xmath72 such that @xmath82 .
since the total energy of the system transforms as @xmath83 , the change in the chemical potential is @xmath84 . inserting these into eq .
( [ landaufermithermpot ] ) gives the relativistic landau effective - mass relation @xcite @xmath85 note that in the right - hand side @xmath86 also depends on @xmath3 , because the fermi momentum @xmath57 is a , in general complicated , function of @xmath3 .
now for the more general imbalanced case . since the lorentz transformation is applied uniformly to each species , according to the symmetry of the formula , we can expect that the contribution from each species in eq .
( [ reducedthermpotexpansion ] ) should vanish , namely , @xmath87 then we get the general expression for the effective mass , @xmath88,\ ] ] where @xmath89 is the density of states at the fermi energy of species @xmath22 .
the arguments in @xmath90 emphasize that , in the imbalanced case , the interaction may depend not only on the angle between the momenta but also on the absolute value of the fermi momentum of each species .
in fact , the same expression of @xmath28 can be obtained more strictly by considering the addition of a single particle to the system and comparing the energy increase in two different frames .
we can introduce the average effective mass , which is of practical importance , as @xmath91 with @xmath92 the total number of species .
correspondingly , the average chemical potential is @xmath93 . in order to get a similar expression as in eq .
( [ rel - landaumass ] ) , we introduce the relative weight of each species as @xmath94 such that @xmath95 , and then generalize the landau effective - mass relation with average values as @xmath96,\ ] ] where the double - sum term , which can be defined again as @xmath86 , plays the role of the landau parameter in the balanced case as in eq .
( [ rel - landaumass ] ) , therefore we obtain now @xmath97 which , obviously , reduces to eq .
( [ landaufermiparameter ] ) as the system becomes balanced .    in the following sections the effective mass will be derived from a microscopic rpa calculation . in particular , in sec .
[ sec - thermpotential ] it will be determined from the specific heat .
the specific heat ( per volume ) for low temperatures can be expressed in terms of the effective masses by @xmath98_{v , n}\nonumber\\ & \simeq\sum_\sigma\frac{k_\sigma m^*_\sigma}{6}t=\sum_\sigma\frac{\pi^2}{3}n_\sigma(0)t,\end{aligned}\ ] ] where it was used that @xmath99 and for the temperature derivative the sommerfeld expansion was used , c.f .
( [ eq : ap_sommerfeld_exp ] ) .
this result is a trivial generalization of the single species or balanced case , since the effective mass already contains the effect of interactions , such that the contribution to the heat capacity of a single quasiparticle is additive .
our starting point is the partition function @xmath100 containing the euclidean qcd lagrangian density fixed in a linear gauge @xmath101 , @xmath102 where @xmath103 is the quark field of flavor @xmath104 ( the color @xmath105 and spin @xmath15 degrees of freedom are not shown explicitly ) with mass @xmath106 and chemical potential @xmath107 , @xmath108 are the gluon fields , @xmath109 are the ghost fields , @xmath6 is the qcd coupling constant , @xmath110 is a gauge fixing parameter , @xmath111 are the fine - structure constants of the color @xmath112 group , the covariant derivative @xmath113 and the antisymmetric gluon field tensor @xmath114 . for large densities ,
the fermionic degrees of freedom become increasingly important , such that the large-@xmath0 limit could give insight into the behavior of quarks at high densities .
therefore , consider the large-@xmath0 limit with a fixed rescaled coupling @xmath115 , @xmath116 integrating out the fermions and the ghosts gives @xmath117\\ & -\mathrm{tr}_c\ln\left(\partial_{\mu}f_\mu\delta^{ab}+\mathrm{g}n_f^{-\frac{1}{2}}f^{abc}a_\mu^cf_\mu\right)\\ & + \frac{1}{2}a_\mu^ad^{-1}_{0,\mu\nu}\delta^{ab}a_\nu^b\\ & + \frac{1}{2}\mathrm{g}n_f^{-\frac{1}{2}}\left(\partial_{\mu}a_\nu^a-\partial_{\nu}a_\mu^a\right)f^{abc}a_\mu^ba_\nu^c\\ & + \frac{1}{4}\mathrm{g}n_f^{-1}f^{abc}f^{ade}a_\mu^ba_\nu^ca_\mu^da_\nu^e,\end{aligned}\ ] ] where the subscripts @xmath118 explicitly indicate that the trace should also be taken over color , flavor and spin space .
when we expand the first line in @xmath0 , the first - order contribution will vanish due to conservation of color , i.e. , @xmath119=0 $ ] .
only the second - order expansion in @xmath0 of the first line will give a contribution .
thus expansion in powers of @xmath0 gives @xmath120 where we defined the polarization tensor @xmath121 from the contribution of the quark with flavor @xmath104 as @xmath122 which is shown diagrammatically in fig .
[ fig - polarizationtensor ] .
the polarization tensor @xmath121 is diagonal in color space , as can be seen explicitly form eq .
( [ eq : def_polarization_tensor ] ) . upon integrating out the gluons
the thermodynamic potential density is @xmath123\nonumber\\ & \left.+\frac{n_g}{2}\mathrm{tr}\ln\left(1+d_0\sum_f\pi_f\right)\right\},\label{eq : therm_pot}\end{aligned}\ ] ] where @xmath124 is the number of colors , @xmath125 is the number of gluons , and we use the qcd values @xmath126 and @xmath127 in the following .
the first two terms are the ideal fermi and bose gas contributions to the thermodynamic potential density , where the term @xmath128 cancels the two unphysical degrees of freedom from @xmath129 .
the last term is the rpa correction ( ring sum ) to the thermodynamic potential density , as is shown in fig  [ fig - thermopotential ] .
since we are mainly interested in the temperature dependence of @xmath130 and the fact that the above definition contains ( divergent ) zero - temperature contributions , the @xmath131 expression will be subtracted , i.e. , we consider @xmath132 . in sec .
[ sec - weakcoupling ] an example is considered to check the expression for the landau effective mass analytically and in sec .
[ sec - thermpotential ] the above thermodynamic potential density is calculated numerically for the case of a two - flavor imbalanced quark system .     to the first order . ]     or rpa approximation where the double gluon line signifies the dressed gluon propagator . ]
in the limit of weak coupling the interaction contribution to the thermodynamic potential density is given solely by the exchange diagram @xmath133 using the definition of the free gluon propagator in the lorentz gauge , this can be written as @xmath134,\end{aligned}\ ] ] where the fermi and bose distributions are defined as @xmath135}+1},\quad n_b^s(\mathbf{p})=\frac{1}{e^{\beta[s\epsilon^g(\mathbf{p})]}-1},\ ] ] respectively , the explicit expression for @xmath136 is shown in eq .
( [ polarizationtensor ] ) , the matsubara sum over the bosonic frequencies @xmath137 was performed , @xmath138 denotes the free quark dispersion while @xmath139 is the free gluon dispersion , and the function @xmath140 was defined as @xmath141 which will soon be shown to play the role of the interaction between two quarks ( @xmath142 ) , two anti - quarks ( @xmath143 ) or a quark and an antiquark ( @xmath144 or @xmath145 ) mediated by the emission or absorption of a gluon ( depending on the sign of @xmath146 relative to @xmath147 and @xmath148 ) .
the quasiparticle energy and effective interaction can be obtained by varying the thermodynamic potential with respect to the fermi distribution . in the flavor - imbalanced case ,
the distributions depend only on the flavor but not on the color and spin indices . to distinguish between the different contributions for each particle species
it can be checked that the multiplication factor @xmath149 could be written as @xmath150 , where @xmath15 and @xmath151 are just dummy indices denoting the spin degrees of freedom , which are helpful to arrive at the following expressions .
note that @xmath152 is independent of the spin index , and in fact the spin degrees of freedom has already been summed over in the expression for @xmath140 .
explicitly we find for these quantities    @xmath153-\sum_{s'}\frac{g^2}{4}\sum_{c'}t_{cc'}^at_{c'c}^a\sum_{s_2,s_3=\pm1}\int\frac{\mathrm{d}^3p}{(2\pi)^3}\mathcal{f}^f_{+,s_2,s_3}(\mathbf{k},\mathbf{q})\nonumber\\ & \qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\times[(1+n_b^{s_3}(\mathbf{k}-\mathbf{q})(1-n_f^{s_2}(\mathbf{q}))+n_b^{s_3}(\mathbf{k}-\mathbf{q}))n_f^{s_2}(\mathbf{q})],\label{quasiparticle - energy}\\ & f_{\{c , f , s\};\{c',f',s'\}}(\mathbf{k},\mathbf{k}')\equiv\frac{\delta^2\omega}{\delta n_{c , f , s}^+(\mathbf{k})\delta n_{c',f',s'}^+(\mathbf{k}')}=\frac{g^2}{4}t_{cc'}^at_{c'c}^a\delta_{ff'}\sum_{s_3=\pm1}\mathcal{f}^f_{+,+,s_{3}}(\mathbf{k},\mathbf{k}'),\label{quasiparticle - int}\end{aligned}\ ] ]    where @xmath154 and @xmath155 , according to the definitions in appendix  [ app - euclidean ] .    in the quasiparticle energy eq .
( [ quasiparticle - energy ] ) the first term is the noninteracting part , while the second term corresponds to the quark self - energy , as is shown in appendix  [ app - quark - selfenergy ] .
the angular averaged interaction parameters are defined by the coefficients of the legendre polynomial expansions of the effective interaction eq .
( [ quasiparticle - int ] ) evaluated on the fermi sphere , by denoting @xmath156 ,    @xmath157    which gives for the zeroth and first interaction parameters @xmath158,\nonumber\\ f_{1;f}&=\frac{g^2n_gn_s}{8n_c}\frac{3}{\mu_{0f}^2}\frac{\lambda^2 + 2m_f^2}{2k_f^2}\nonumber\\ & \qquad\qquad\times\left[1+\frac{\lambda^2 + 2k_f^2}{4k_f^2}\ln\left(\frac{\lambda^2}{\lambda^2 + 4k_f^2}\right)\right],\label{weakcouplandauparam}\end{aligned}\ ] ] where @xmath159 is the number of spin degrees of freedom , and a regulatory gluon mass @xmath160 was introduced , i.e. , @xmath161 .
we see that in this approximation the interaction is color- and spin - independent . for the sake of simplicity , therefore , we from now on omit the unnecessary color and spin indices , and use @xmath162 and @xmath163 for the fermi momentum and the non - interacting fermi energy of the quark with flavor @xmath104 .
the effective mass of the balanced case , eq .
( [ rel - landaumass ] ) , can be obtained from the interaction parameters eq .
( [ weakcouplandauparam ] ) , and eq .  ( [ landaufermiparameter ] ) , by expanding to the first order in the coupling constant @xmath164 where it was used that @xmath165 as follows from eq .
( [ quasiparticle - energy ] ) . for the imbalanced case
, a similar result can be obtained for each species according to eq .
( [ imbalancemass ] ) , @xmath166 where the kronecker delta @xmath167 in the interaction simplifies the expression considerably , such that we obtain the similar expression as in eq .
( [ effectivemassexpand ] ) .    as a special case
, we can use the above results to discuss a two - flavor imbalanced system , where the average effective mass and chemical potential can be defined as @xmath168 and @xmath169 with the subscripts @xmath170 and @xmath171 for the majority and minority flavors , respectively .
the imbalance can be quantified as @xmath172 . for small imbalance
, the factors @xmath173 can be expanded in @xmath174 with respect to the value in the balanced case , namely @xmath175 .
note that , the linear coefficients of the expansion can , in general , be different for each flavor .
we obtain the expression of @xmath56 for small @xmath174 @xmath176,\end{aligned}\ ] ] where the @xmath177 term is kept since the linear term may vanish for the symmetric case , namely @xmath178 , and we have replaced @xmath179 with @xmath3 because the difference is in higher order of @xmath6 .
compared with eq .
( [ averageimbalancemass ] ) , the effective landau parameter for the average effective mass is @xmath180.\end{aligned}\ ] ] similarly , we can obtain the effective mass difference @xmath181 as a function of @xmath174 @xmath182.\end{aligned}\ ] ] it is quite natural to find @xmath183 in the leading order .    the consistency of the above results will now be verified to the lowest order in the coupling constant using the quark self - energy . starting from the dispersion and self - energy of eq .
( [ quasiparticle - energy ] ) , the effective mass in eq .
( [ effectivemassexpand ] ) or eq .
( [ imbalancemassexpand ] ) can also be derived in a different manner .
the quasiparticle pole of the quark propagator is @xmath184 where @xmath185 is the renormalized positive - energy projected self - energy .
the effective mass corresponding to this pole is most easily defined from the fermi velocity of the quasiparticle @xmath186_{p = k_f}.\ ] ] inserting the definition of the quasiparticle pole into the above gives , using @xmath187 , @xmath188_{p = k_f}^{-1}+\mathcal{o}(g^4)\nonumber\\ = & \mu_f-\left[\frac{\partial\epsilon_{0f}(\mathbf{p})\sigma_{c , f , s}^+(\epsilon_{0f}(\mathbf{p})-\mu_f,\mathbf{p})}{\partial\epsilon_{0f}(\mathbf{p})}\right]_{p = k_f}+\mathcal{o}(g^4).\label{effectivemassnew}\end{aligned}\ ] ] using that the explicit form of the self - energy in the zero - temperature limit can be written as ( c.f .
( [ quarkselfenergy ] ) ) @xmath189 and the relation eq .
( [ eq : ap_eff_int_transf ] ) between the vector derivatives of the effective interaction , the derivative in eq .
( [ effectivemassnew ] ) can be rewritten as    @xmath190_{p = k_f}=\frac{\mu_{0f}}{k_f}\left[\sum_{c',f',s'}\int\frac{\mathrm{d}^3q}{(2\pi)^3}\hat{\mathbf{p}}\cdot\frac{\partial\epsilon_{0f}(\mathbf{p})f_{\{c , f , s\};\{c',f',s'\}}(\mathbf{p},\mathbf{q})}{\partial\mathbf{p}}n_{f'}^+(\mathbf{q})\right]_{p = k_f}\\ = & -\frac{\mu_{0f}}{k_f}\int\frac{\mathrm{d}q}{2\pi^2}q^2\epsilon_{0f}(\mathbf{q})\delta(q - k_{f})\frac{1}{3}\left[3\int\frac{\mathrm{d}\hat{\mathbf{q}}}{4\pi}\hat{\mathbf{p}}\cdot\hat{\mathbf{q}}\sum_{c',f',s'}f_{\{c , f , s\};\{c',f',s'\}}(\mathbf{p},\mathbf{q})\right]_{p = k_f } = -\frac{k_f\mu_{0f}^2}{2\pi^2}\frac{f_{1;\{c , f , s\}}(k_f)}{3},\end{aligned}\ ] ]    where partial integration was used and @xmath191 at zero temperature . therefore eq .
( [ effectivemassnew ] ) is identical to eq .
( [ imbalancemassexpand ] ) , which shows that the phenominological landau argument is indeed consistent with the microscopic diagrammatic calculation .
before we start the calculation of the temperature dependence of the large-@xmath0 thermodynamic potential , it is useful to examine the dressed gluon propagator , which is a crucial ingredient of the rpa theory . the most physical gauge to study the propagator is the coulomb gauge , since its form results from considering linear response @xcite .
this is also by far the most used gauge in condensed - matter theory . in this gauge
the gluon propagator is @xcite @xmath192 where @xmath193 is the three - dimensional transverse projector defined in appendix  [ app - green ] , @xmath194 and @xmath195 are related to the longitudinal and transverse projections of @xmath196 including the contributions of various species as well as the vacuum .
more details can be found in appendix  [ app - gluon - selfenergy ] .
consider the spectral functions of the transverse and longitudinal propagator , which are defined by @xmath197.\ ] ] in general the spectral function depends on the gauge , however , the positions of the poles are gauge independent , and the coulomb gauge has the additional property @xmath198 for @xmath199 as required of a physical spectral function .
for the balanced case , the spectral functions are shown in fig .  [
fig : rhotl ] for several values of @xmath3 and in fig .
[ fig : rhotltemp ] for several values of @xmath1 .
the latter includes a small-@xmath1 correction to @xmath194 and @xmath195 , see eq .
( [ eq : ap_pi_temperature_correction]).-correction the spectral functions will no longer satisfy @xmath200 . ] the transverse and longitudinal plasmon modes and the large contribution due to the decay of the gluon into the particle - hole continuum ( @xmath201 ) are clearly visible .
note that massless quarks are used in the limit of high density , c.f .
appendix  [ app - gluon - selfenergy ] . however , a small but finite quark mass is necessary to renormalize the real part of the polarization tensor of the vacuum , as shown in eq .
( [ eq : ap_renormalised_vacuum ] ) .
since the vacuum contribution plays no significant role in the following calculation , we simply take this nonzero quark mass @xmath202 in the vacuum term the same for different flavors .    for the two - flavor imbalanced case ,
it is clear from eq .
( [ eq : therm_pot ] ) that in the leading - order correction each flavor contributes separately to the polarization tensor . because of the above setting , there is no mass imbalance in the present system , therefore it is enough to consider only positive @xmath174 due to the symmetry . following the notation in sec .
[ sec - weakcoupling ] , all the above results can be easily generalized to such an imbalanced system by using @xmath203 to replace @xmath204 , and @xmath205 for @xmath206 . as @xmath207
, the system reduces to the balanced case .
now the low temperature condition requires @xmath208 , such that @xmath174 can not be too close to @xmath209 , namely the extremely imbalanced case .
the spectral function for a two - flavor imbalanced system is shown for various @xmath174 in fig .
[ fig : rhotlh ] . comparing with fig .
[ fig : rhotltemp ] , we see that increasing @xmath174 has a similar effect as increasing @xmath1 .    the dispersion relations of the modes can be found by solving @xmath210=0.\label{eq : re_invprop_zero}\ ] ] in the transverse case there is a single solution @xmath211 , however , in the longitudinal case there are two solutions @xmath212 , as shown in fig .
[ fig : dispersion ] .
note that @xmath213 is not a real propagating mode because the imaginary part in the region @xmath201 is large due to particle - hole creation processes , as explained in appendix  [ app - gluon - selfenergy ] .
the two plasmon modes @xmath214 and @xmath215 approach the so - called plasma frequency @xmath216 as @xmath217 .
the plasma frequency of the transverse and longitudinal mode can be found by expanding the inverse propagators for small @xmath218 and small @xmath219 , which in the zero - temperature limit gives for both cases @xmath220 whose solution for small coupling constant yields @xmath221 @xcite .
here @xmath222 is the gluon thermal mass whose expression is obtained in the hard dense and hard thermal loop approximation as shown in eq .
( [ thermalmass ] ) . in the limit @xmath223 all solutions reduce to @xmath224 ,
( [ eq : ap_fg_large_pomega_behavior ] ) .    furthermore , for the limit @xmath225 the vacuum becomes increasingly dominant such that eq .
( [ eq : re_invprop_zero ] ) has a zero at large @xmath226 called the landau pole , @xmath227 @xcite , c.f .
( [ eq : ap_fg_large_pomega_behavior ] ) , which , however , plays no role for our purpose as we are interested in the low - temperature behavior of the theory that is hardly influenced by the high - energy behavior of the gluon propagator .
the dispersion relations for the case with nonzero @xmath1 or @xmath174 is similar but just with a little higher @xmath228 .
in fact , as a generalization of eq .
( [ thermalmass ] ) , the thermal mass of a two - flavor imbalanced quark system reads @xmath229,\ ] ] where @xmath174 thus plays the same role as @xmath230 .
the full ideal - gas contribution to the thermodynamic potential density with the zero - temperature contribution subtracted is ( c.f .
appendix  [ app - thermpotential ] ) @xmath231 the first part is the contribution of an ideal massless fermi gas , while the second part is the stefan - boltzmann law of an ideal bose gas .
the rpa correction to @xmath130 , the last term in eq .
( [ eq : therm_pot ] ) , can be written as @xmath232 where the lorentz gauge is used , while for the coulomb gauge the second identity is not valid but the last result is still the same , which is a consequence of the gauge invariance of the thermodynamic potential density . using contour deformations to carry out the matsubara sum , as shown in appendix  [ app - matsubara ] , we obtain @xmath233\\ & \times\{\im[\ln\tilde{f}(\omega_+,\mathbf{q})]+2\im[\ln\tilde{g}(\omega_+,\mathbf{q})]\},\end{aligned}\ ] ] where @xmath234 , and similar for @xmath235 by replacing @xmath194 with @xmath195 . the temperature dependence can be obtained by subtracting the zero - temperature contribution .
since @xmath194 and @xmath195 contain corrections of order @xmath236 , the leading - order correction in @xmath237 comes from two parts , @xmath238\mathrm{sgn}[\im(\tilde{f})^0]\right.\\ & \left.+2\arctan\frac{\im(\tilde{g}^0)}{\re(\tilde{g}^0)}+2\pi\theta[-\re(\tilde{g}^0)]\mathrm{sgn}[\im(\tilde{g}^0)]\right)\\ & \left.\frac{1}{2}[\im(\ln\tilde{f}^t)-\im(\ln\tilde{f}^0)+2\im(\ln\tilde{g}^t)-2\im(\ln\tilde{g}^0)]\vphantom{\frac{1}{2}}\right ] , \end{aligned}\ ] ] where the superscript @xmath1 or @xmath239 means the corresponding terms are taken at nonzero @xmath1 or @xmath131
. we will refer to the first part as the @xmath240 term and the second as the non-@xmath240 term .
the leading correction from the non-@xmath240 term can be shown to be proportional to @xmath236 and is not of great interest in our study , since we will concentrate on the anomalous and dominant @xmath1 dependence , which is a consequence of the @xmath240 term . in the @xmath240 term
, the arctangent terms can be interpreted as contributions due to production and decay of thermal gluons , because of their dependence on the imaginary part of the gluon self - energies @xmath194 and @xmath195 , while the theta function terms are interpreted as a correction to the ideal gas law due to thermal plasmon modes .
the frequency integral over the theta function can be performed explicitly . in both the transverse and longitudinal case
the sign of the imaginary part is positive when the real part is negative , such that after the frequency integration the result of the integral is proportional to @xmath241 where @xmath242 and @xmath243 signifies the region where @xmath244 . in the limit of low temperature
this integral will go to a constant .
the theta - function contribution can thus be neglected since it is of higher order in the temperature than is of interest to us here .
next we perform the integrals over the arctangents , whose structure at @xmath245 is shown in fig .
[ fig : arctan ] .
note that , for the study of the leading - order @xmath1 corrections , it is not necessary to include the @xmath236 term in @xmath194 and @xmath195 since the integral with @xmath246 at low @xmath1 is already in the order of @xmath236 .
the dominant contribution for small temperatures ( @xmath2 ) comes from the frequency integration over the domain @xmath247 $ ] for the case @xmath248 , which is due to particle - hole creation .
for the two - flavor balanced case ( @xmath249 , @xmath250 ) , it was found numerically that the integral in the limit of small temperatures is @xmath251 for the electric ( longitudinal ) gluons and @xmath252 for the magnetic ( transverse ) gluons , with @xmath253 , @xmath254 , @xmath255 , and @xmath256 .
these results are quite close to the analytic results for the leading - order correction terms at small @xmath6 and @xmath1 obtained in ref .
@xcite , where @xmath257/6\approx0.4460 $ ] , @xmath258 , @xmath259/18\approx0.4032 $ ] , and @xmath260 .
for the imbalanced case , the above result for the transverse gluons can be generalized to @xmath261,\ ] ] where the contributions from both flavors with chemical potentials @xmath262 are additive .
the longitudinal part needs some further discussion . since the static long - wavelength longitudinal modes are screened ( c.f .
appendix  [ app - gluon - selfenergy ] ) , the two chemical potentials contributing to @xmath263 can not be separated even in the lowest - order term .
this is different from the result obtained in weak coupling , as shown in eq .
( [ quasiparticle - int ] ) , where no interaction between the two flavors is involved .
therefore it is not surprising to find that a simple generalization of the balanced case , as @xmath264 , does not fit well with the numerical results . to obtain a reasonable ansatz for the longitudinal part
, we integrate the corresponding arctangent term of the imbalanced case in the low temperature limit up to @xmath265 to obtain @xmath266,\ ] ] where we see that the factor @xmath267 is a mixture effect of the two flavors .
this expression fits the numerical results very well with an error of only @xmath268 .
furthermore , just based on the numerical data , we find another ansatz which fits the results even better @xmath269 where the mixture effect is shown implicitly in the logarithm , since the denominator becomes dimensionless by canceling with the average chemical potential @xmath270 in the numerator .    the specific heat ( per volume ) at fixed volume and particle number is @xcite @xmath271 ^ 2}{(\partial n_f/\partial\mu_f)_t},\ ] ] where the entropy density @xmath272 and particle number density @xmath273 . in the low - temperature limit the second term can be neglected ,
therefore the specific heat can be obtained , using eqs .
( [ ct ] ) and ( [ cl ] ) , as    @xmath274,\ ] ]    where the specific heat @xmath275 of an ideal gas has been subtracted .
the effective mass and the first landau parameter can be determined in the high - density limit by comparing eq .
( [ cv ] ) and eq .
( [ heatcapacity ] ) , @xmath276,\end{aligned}\ ] ] where we used @xmath277 in the massless limit at high density .
the first term in the effective mass is due to the ideal - gas specific heat .
however , as pointed out before , the factor @xmath278 shows the mixing of the two flavors , unlike the weak - coupling results obtained in sec .
[ sec - weakcoupling ] , e.g. , eq .
( [ imbalancemassexpand ] ) .
this is because , with the rpa correction , we incorporate the sum of an infinite chain of gluon self - energies , which incorporates interactions between different quark flavors . since the longitudinal gluon is screened in the static long - wavelength limit , the mixing is even present in the low-@xmath1 and small-@xmath6 limit .
the numerical fit in eq  ( [ clnew ] ) also provides another expression for the effective mass , @xmath279.\end{aligned}\ ] ] as expected , all the above results return to the balanced case as @xmath207 .
as mentioned previously , the contribution from the non-@xmath240 term is not included , which acts as a constant shift on the @xmath280 and @xmath281 factors . in fig .
[ fig : mass ] the change of effective mass due to the interaction is given as function of imbalance at various temperatures .
finally , we emphasize again that @xmath174 should not be too close to @xmath209 even though the divergence from @xmath282 is suppressed by the prefactor @xmath283 , because in such an extremely imbalanced case the condition @xmath284 is not satisfied for the minority flavor .
we have calculated the thermodynamic potential perturbatively in the large-@xmath0 limit and the effective mass of the quarks is determined by using fermi - liquid theory for an imbalanced cold dense quark system .
the temperature dependence is obtained by using the gluon self - energy , from which the contributions from transverse and longitudinal gluons are explicitly shown .
for the two - flavor imbalanced quark system , the effective mass is obtained both analytically within the weak - coupling limit , and numerically within the rpa approximation .
we find that , in contrast to the weak - coupling result , where the effective mass of each flavor is independent of each other due to the lack of an inter - flavor interaction , the effective mass obtained from the rpa calculation depends on the chemical potentials of both flavors , which is a consequence of the static screening of the electric gluons at long wavelength .    from the rpa results , the logarithmic dependence on temperature of the specific heat and effective mass signals a breakdown of fermi - liquid theory at zero - temperature .
non - fermi - liquid behavior , due to unscreened long - range magnetic interactions , was already discussed several decades ago for the case of the electron gas @xcite . in the large-@xmath0
limit the qcd thermodynamic potential is essentially the same as that of qed , apart from group theory factors .
the non - abelian effects of qcd only show up if gluon self - interaction corrections are included .
the logarithmic behavior in the balanced case has previously been seen in analytic and numerical calculations of the qcd thermodynamic potential and specific heat using the large-@xmath0 limit , dimensional reduction , and hard - dense loop qcd perturbation theory @xcite .
we expect that such a logarithmic dependence shall have important effects on the imbalanced qcd phase diagram at low temperatures and directly influence the properties of quark matter in the core of neutron stars . to better understand these effects , it is necessary to go further , such as including the gluon self - interaction and extending our discussion to the three - flavor imbalanced case . we hope progress along these directions will be achieved in the near future .
this work is supported by the stichting voor fundamenteel onderzoek der materie ( fom ) and the nederlandse organisatie voor wetenschaplijk onderzoek ( nwo ) .
the four - momentum vectors in euclidean space - time will be written with capital letters @xmath285 , while for the wick - rotated case ( @xmath286 ) the roman capital letters @xmath287 are used .
three momentum vectors are written in bold face @xmath288 and its length as @xmath218 . and
@xmath289 are sometimes written as @xmath290 for short .
the euclidean gamma matrices in the standard representation are @xmath291 @xmath292 all the above entries are @xmath293 matrices and @xmath294 are the pauli spin matrices @xmath295 the gamma matrices obey the following relations @xmath296 the charge conjugation matrix is @xmath297 which satisfies @xmath298    the eight generators of the fundamental representation of @xmath299 are taken to be @xmath300 @xmath301 @xmath302 @xmath303 and have been normalized according to @xmath304=\frac{1}{2}\delta^{ab}$ ] . in general , the product of the generators is @xmath305 furthermore , one will frequently encounter the following group - theory factors @xmath306 @xmath307 in the above expressions @xmath124 is the number of colors and @xmath308 the number of gluons .
fourier transforms are normalized as @xmath309 where @xmath310 are the fermionic matsubara frequencies , @xmath311 is the imaginary time phase space volume and @xmath312 and @xmath313 .
the dirac equation in euclidean space is @xmath314 where @xmath315 . in momentum space
the dirac equation reads @xmath316 the eigenvalues of the matrix @xmath317 are @xmath318 , since @xmath319 .
the eigenspinors corresponding to these eigenvalues are @xmath320 note that @xmath317 is not a hermitian matrix and that @xmath321 satisfies the same equation as @xmath322 .
however , they can also be viewed as the the eigenspinors of a hermitian matrix @xmath323 the positive and negative energy eigenspinors have the following form @xmath324 where @xmath325 and @xmath326 .
they are orthonormal in the sense that @xmath327 and satisfy the completeness relation @xmath328=1.\ ] ]    using the above eigenspinors the positive and negative energy projectors can be defined as @xmath329 which in the massless case reduce to @xmath330 the helicity projection operators project the spin along the momentum of the particle and read @xmath331 where @xmath332 is the spin operator .
note that @xmath333=0 $ ] .
the green s functions of quarks and gluons in euclidean space are defined by @xmath334 in momentum space , the free quark propagator is @xmath335 in the lorentz gauge @xmath336 , the momentum - space free gluon propagator is @xmath337=\frac{p_{\mu\nu}}{q^2}+\frac{\xi}{q^2}\frac{q_\mu q_\nu}{q^2},\ ] ] while in the coulomb gauge @xmath338 it is @xcite @xmath339 where @xmath110 is a gauge - fixing parameter and the projectors can be written as @xmath340 note that the free gluon propagator in the lorentz gauge satisfies @xmath341 , such that in the landau gauge ( @xmath342 ) it is purely four - momentum transverse , i.e. , @xmath343 .
the free gluon propagator in the coulomb gauge satisfies @xmath344 , such that in the landau gauge it is three - momentum transverse .
these propagators are diagonal in color space .
in this section the lorentz transformation properties of some quantities are summarized , such as the thermodynamic potential , the distribution function and the effective interaction .
the lorentz transformation to a frame moving with velocity @xmath74 is @xmath345 such that the four - momentum @xmath346 of a particle transforms as @xmath347 not to be confused with the dirac matrices @xmath348 , the @xmath349 used in this section is the lorentz factor @xmath65 .
the invariance of the thermodynamic potential density under lorentz transformations can be shown using the stress - energy tensor .
consider the change in the thermodynamic potential density under a lorentz transformation from the rest frame to a frame moving with velocity @xmath74 .
since an interacting gas of quarks in the rest frame is specified only by an energy density @xmath350 and a pressure @xmath351 , the stress - energy tensor is diagonal and of the form @xmath352 such that under a lorentz transformation the thermodynamic potential density @xmath353 changes as @xmath354 the transformation of the chemical potential ( @xmath355 ) follows from the lorentz boosted total energy of the system @xmath356 , where @xmath357 is the total momentum of the system in the rest frame , which results in @xmath358 .
the transformation of the density is due to a lorentz contraction in the volume @xmath359 .
thus the total change in @xmath360 is @xmath361 , such that under a lorentz transformation @xmath362 . for an alternative derivation which uses that the pressure transforms the same way as a force per area , see ref .
@xcite .
consider the lowest - order correction to the free thermodynamic potential density @xmath364 note that the distribution function @xmath365 is a lorentz invariant , which can be easily derived from the fact that the number of particles in a volume @xmath366 of phase space is invariant under lorentz transformations @xcite , i.e. , @xmath367 where the tilde signifies the lorentz transformed quantity .
subsequently , it is possible to derive a transformation law for @xmath363 by using that the distribution , the thermodynamic potential density and @xmath368 are lorentz invariant .
it follows that @xmath369 should be lorentz invariant , giving @xmath370 where @xmath371 .
expand @xmath372 to the lowest order in @xmath74 , using @xmath373 , @xmath374 , and @xmath375 ,    @xmath376\left[1-\mathbf{v}\cdot\frac{\mathbf{p}'}{\epsilon_{0\sigma'}(\mathbf{p}')}\right]+\mathcal{o}(v^2)\\ = & \tilde{f}_{\sigma\sigma'}(\mathbf{p},\mathbf{p}')-\epsilon_{0\sigma}(\mathbf{p})\mathbf{v}\cdot\frac{\partial f_{\sigma\sigma'}(\mathbf{p},\mathbf{p}')}{\partial\mathbf{p}}-\epsilon_{0\sigma'}(\mathbf{p}')\mathbf{v}\cdot\frac{\partial f_{\sigma\sigma'}(\mathbf{p},\mathbf{p}')}{\partial\mathbf{p}'}-f_{\sigma\sigma'}(\mathbf{p},\mathbf{p}')\left[\mathbf{v}\cdot\frac{\partial\epsilon_{0\sigma}(\mathbf{p})}{\partial\mathbf{p}}-\mathbf{v}\cdot\frac{\partial\epsilon_{0\sigma'}(\mathbf{p}')}{\partial\mathbf{p}'}\right]+\mathcal{o}(v^2)\\ = & \tilde{f}_{\sigma\sigma'}(\mathbf{p},\mathbf{p}')-\mathbf{v}\cdot\left[\frac{\partial\epsilon_{0\sigma}(\mathbf{p})f_{\sigma\sigma'}(\mathbf{p},\mathbf{p}')}{\partial\mathbf{p}}+\frac{\partial\epsilon_{0\sigma'}(\mathbf{p}')f_{\sigma\sigma'}(\mathbf{p},\mathbf{p}')}{\partial\mathbf{p}'}\right]+\mathcal{o}(v^2).\end{aligned}\ ] ]    if it is assumed that to the lowest order the interaction does not depend on any distribution functions , i.e. , @xmath377 , the above implies @xmath378 note that the above derivation is similar to that given in ref .
in this section it is summarized how to turn the minkowski quantum field theory ( qft ) of quantum chromodynamics into a euclidean statistical field theory ( sft ) suitable for studying the dynamical properties of a many - particle system .
the starting point is the gauge - fixed path integral for qcd @xmath379 where @xmath380 are the gluon fields , @xmath381 and @xmath382 the quark fields , @xmath383 and @xmath384 the ghost fields and the qcd lagrangian in minkowski space fixed in a linear gauge ( @xmath385 ) is given by @xmath386\nonumber\\ & -\frac{1}{2}a_a^\mu(\partial_{\mu}\partial_\nu-\partial^2\eta_{\mu\nu})a_a^\nu\nonumber\\ & -\bar{\eta}^a\partial_\mu\partial^\mu\eta^a - gf^{abc}(\bar{\eta}^a\partial^\mu\eta^b)a_\mu^c\nonumber\\ & -\frac{1}{2}g(\partial^\mu a_a^\nu-\partial^\nu a_a^\mu)f^{abc}\eta_{\mu\sigma}\eta_{\nu\rho}a_b^\sigma a_c^\rho\nonumber\\ & -\frac{1}{4}g^2f^{abc}f^{ade}\eta_{\mu\sigma}\eta_{\nu\rho}a_a^\mu a_c^\nu a_d^\sigma a_e^\rho+\frac{1}{2\xi}(f^\mu a_\mu^a)^2,\end{aligned}\ ] ] where @xmath103 is the quark field with flavor @xmath104 , and the summation over color and spin indices are shown implicitly .
the metric was chosen to be @xmath387 and the minkowski gamma matrices in the standard representation are @xmath388 the gamma matrices satisfy @xmath389    by performing a wick rotation the above quantum field theory can be transformed into a statistical field theory . a wick rotation amounts to taking an analytic continuation from real time to imaginary time ( @xmath390 ) , which turns the minkowski metric @xmath391 into the euclidean metric @xmath392 . to do this consistently the zeroth component of all four - vectors
need to change accordingly .
the procedure is most easily understood by considering the length of the position four - vector @xmath393 where it is seen that the minus sign of the minkowski metric is absorbed in the definition of the euclidean four - vectors @xmath394 . in euclidean space
no distinction is made between upper and lower indices .
a simple way to obtain the euclidean form of a vector is to multiply the contravariant vector by the matrix @xmath395 or the covariant vector by @xmath396 and set @xmath390 .
note that the spatial components do not change .
for example , the euclidean position four - vector , the four - divergence and the zeroth gamma matrix are in terms of their minkowski definitions eq .
( [ eq : ap_minkowski_gamma ] ) @xmath397    generalizations to tensors is straightforward and follows for instance from the example @xmath398 . using the above procedure to find the euclidean versions of all tensors , the partition function is easily found from eqs .
( [ eq : ap_minkowski_path_int ] ) and ( [ eq : ap_minkowski_qcd_lagrangian ] ) by setting the tensors to their euclidean versions , taking @xmath399 and @xmath400 .
the partition function is @xmath401 where the euclidean action is defined as @xmath402 with the euclidean lagrangian @xmath403\\ & + \frac{1}{2}a_\mu^a(\partial_\mu\partial_\nu-\partial^2\delta_{\mu\nu})a_\nu^a\\ & + \bar{\eta}^a\partial^2\eta^a+gf^{abc}(\bar{\eta}^a\partial_\mu\eta^b)a_\mu^c\\ & + \frac{1}{2}g(\partial_\mu a_\nu^a-\partial_\nu a_\mu^a)f^{abc}a_\mu^ba_\nu^c\\ & + \frac{1}{4}g^2f^{abc}f^{ade}a_\mu^ba_\nu^ca_\mu^da_\nu^e-\frac{1}{2\xi}(\partial_\mu a_\mu^a)^2.\end{aligned}\ ] ] in the above the tensors are all euclidean but the index @xmath11 has been dropped for convenience , the conjugate fields are now defined as @xmath404 and @xmath405 and the chemical potential has been added as the lagrange multiplier of the density @xmath406 . additionally , to complete the connection between qft and sft , the time integration domain is changed to @xmath407 $ ] , where @xmath408 is the inverse temperature @xmath1 of the system . due to the definition of the partition function as a trace over all states , the bosonic ( fermionic ) fields
are required to obey symmetric ( anti - symmetric ) boundary conditions , namely @xmath409 .
at nonzero temperature one needs to calculate matsubara summations which usually can be done using contour integration . in the following an expression
is derived for such summations and the interpretation of the result is discussed .
consider to this end the sum over bosonic matsubara frequencies ( @xmath410 ) of the function @xmath411 @xmath412 where the contour @xmath413 is given in fig .
[ fig : ap_bose_matsubara_contours ] .     along the imaginary axis is used to rewrite the frequency summation as a contour integral .
the contours @xmath414 are the positive and negative energy semicircular contours obtained by adding arcs at infinity to @xmath413.[fig : ap_bose_matsubara_contours ] ]    if @xmath415 when @xmath416 then it is possible to close the contour by adding the arcs of @xmath414 , in which case @xmath417 . in general @xmath418
will only have poles or branch cuts on the real axis , such that the contours can be contracted along the real axis , which gives @xmath419\mathrm{d}\omega\nonumber\\ & -\frac{1}{\pi}\int_0^\infty[1+n_b(\omega)]\im[f(-\omega+i0)]\mathrm{d}\omega,\label{eq : ap_matsubara_bosons}\end{aligned}\ ] ] where @xmath420\equiv f(\omega+i0)-f(\omega - i0)$ ] and it was used that @xmath421 for the contour @xmath422 , such that it is more convenient to take @xmath423 and use @xmath424 . the first line can be interpreted as due to thermal gluons and the second line due to thermal ( anti)gluons and a vacuum contribution .
the fact that the gluon is its own antiparticle will be reflected by @xmath425=-\im[f(\omega+i0)]$ ] .
let us evaluate the above for the specific case @xmath426 $ ] which satisfies @xmath427=-\im[g(\omega+i0)]$ ] .
expanding the logarithm in terms of its real and imaginary parts @xmath428=&\ln|g(z)|+i\arctan\left(\frac{\im[g(z)]}{\re[g(z)]}\right)\\ & + i\pi\theta(-\re[g(z)])\mathrm{sgn}(\im[g(z)]),\end{aligned}\ ] ] then eq .
( [ eq : ap_matsubara_bosons ] ) can be written as @xmath429\\ = & \frac{1}{\pi}\int_0^\infty\mathrm{d}\omega(1 + 2n_{b}(\omega))\left[\arctan\left(\frac{\im\left[g(\omega_+)\right]}{\re\left[g(\omega_+)\right]}\right)\right.\\ & \hphantom{\times}\left.+\pi\theta(-\re[g(\omega_+)])\mathrm{sgn}(\im[g(\omega_+)])\vphantom{\arctan\left(\frac{\im\left[g(\omega_+)\right]}{\re\left[g(\omega_+)\right]}\right)}\right],\end{aligned}\ ] ] which will be used in the calculation of the rpa correction to the thermodynamic potential .    for fermionic matsubara frequencies
a similar derivation can be done , but now the chemical potential is included by writing @xmath430 @xmath431 where the contour @xmath413 is given in fig .
[ fig : ap_fermi_matsubara_contours ] and @xmath432 has poles at the fermionic matsubara frequencies with residue @xmath433 . if @xmath434 when @xmath416 , then it is possible to close the contour as given in fig .
[ fig : ap_fermi_matsubara_contours ] by adding the arcs of @xmath414 and the lower- and upper - boundaries of @xmath435 .
the full matsubara sum can thus be written as a contour integral over @xmath436 . again
if @xmath418 only has poles and branch cuts on the real axis , contracting the contours along the real axis gives @xmath437\mathrm{d}\omega\nonumber\\ & -\frac{1}{\pi}\int_0^\mu[1-n(\mu-\omega)]\im[f(\omega+i0)]\mathrm{d}\omega\nonumber\\ & -\frac{1}{\pi}\int_0^\infty[1-n(\omega+\mu)]\im[f(-\omega+i0)]\mathrm{d}\omega.\end{aligned}\ ] ] where in the first - to - last line it was used that @xmath438 for the contour @xmath439 , such that it is more convenient to write @xmath440 .
similarly for @xmath441 , where @xmath421 , we take @xmath423 in the integral and write @xmath442 . from the above it
is clear that the first line corresponds to particles above the fermi sphere , the second line to those in the fermi sphere and the last line is due to anti - particles and contains a vacuum contribution .
note that the above formula must be used with care , since one has to check if the above simplifications are valid on a case - to - case basis .
but for the fermionic case , where @xmath435 , the contour between @xmath438 , is related to processes inside the fermi sphere.[fig : ap_fermi_matsubara_contours ] ]
the gluon self - energy , or the so - called polarization tensor , is given to the lowest order by @xmath443\nonumber\\ & \equiv\delta^{ab}\pi_{f,\mu\nu}.\end{aligned}\ ] ] the polarization tensor is purely transverse , i.e. , @xmath444 , which can be easily seen using @xmath445 and the cyclicity of the trace .
expanding the polarization tensor in terms of the longitudinal and transverse parts relative to the three - momentum @xmath288 gives @xmath446 where @xmath447 and @xmath448 can be obtained by @xmath449 noticing that even in the imbalanced case , the contribution from each flavor simply adds up to the total polarization tensor , therefore , throughout this section , we will not explicitly write out the subscript @xmath104 , and all the formulae presented here are applicable to any flavor .
however , in the main text , we use @xmath196 , @xmath194 and @xmath195 as the summation over all flavors if the subscript @xmath104 is not explicitly shown .    after performing the matsubara sum
, the individual processes involved can be identified    @xmath450}{i\omega_q-\epsilon_0(\mathbf{p}+\mathbf{q})+\epsilon_0(\mathbf{p})}[n(\mathbf{p}+\mathbf{q})(1-n(\mathbf{p}))-(1-n(\mathbf{p}+\mathbf{q}))n(\mathbf{p})]\right.\nonumber\\ & \qquad\qquad\quad\left.+\frac{\sum_{ss'}[\bar{v}_s(\mathbf{p})\gamma_{\mu}u_{s'}(\mathbf{p}+\mathbf{q})\bar{u}_{s'}(\mathbf{p}+\mathbf{q})\gamma_{\nu}v_s(\mathbf{p})]}{i\omega_q-\epsilon_0(\mathbf{p}+\mathbf{q})-\epsilon_0(\mathbf{p})}[n(\mathbf{p}+\mathbf{q})\bar{n}(\mathbf{p})-(1-n(\mathbf{p}+\mathbf{q}))(1-\bar{n}(\mathbf{p}))]\right.\nonumber\\ & \qquad\qquad\quad\left.+\frac{\sum_{ss'}[\bar{u}_s(\mathbf{p})\gamma_{\mu}v_{s'}(\mathbf{p}+\mathbf{q})\bar{v}_{s'}(\mathbf{p}+\mathbf{q})\gamma_{\nu}u_s(\mathbf{p})]}{i\omega_q+\epsilon_0(\mathbf{p}+\mathbf{q})+\epsilon_0(\mathbf{p})}[(1-\bar{n}(\mathbf{p}+\mathbf{q}))(1-n(\mathbf{p}))-\bar{n}(\mathbf{p}+\mathbf{q})n(\mathbf{p})]\right.\nonumber\\ & \qquad\qquad\quad\left.+\frac{\sum_{ss'}[\bar{v}_s(\mathbf{p})\gamma_{\mu}v_{s'}(\mathbf{p}+\mathbf{q})\bar{v}_{s'}(\mathbf{p}+\mathbf{q})\gamma_{\nu}v_s(\mathbf{p})]}{i\omega_q+\epsilon_0(\mathbf{p}+\mathbf{q})-\epsilon_0(\mathbf{p})}[(1-\bar{n}(\mathbf{p}+\mathbf{q}))\bar{n}(\mathbf{p})-\bar{n}(\mathbf{p}+\mathbf{q})(1-\bar{n}(\mathbf{p}))]\right\}\nonumber\\ = & \frac{g^2}{2}\sum_{s_1,s_2=\pm1}\int\frac{\mathrm{d}^3p}{(2\pi)^3}\frac{\mathrm{tr}[p_e^{s_1}(\mathbf{p})\gamma_0\gamma_{\mu}p_e^{s_2}(\mathbf{p}+\mathbf{q})\gamma_0\gamma_\nu]}{i\omega_q - s_2\epsilon_0(\mathbf{p}+\mathbf{q})+s_1\epsilon_0(\mathbf{p})}[n^{s_2}(\mathbf{p}+\mathbf{q})(1-n^{s_1}(\mathbf{p}))-(1-n^{s_2}(\mathbf{p}+\mathbf{q}))n^{s_1}(\mathbf{p})],\end{aligned}\ ] ]    where @xmath451}$ ] , @xmath452 is the fermi distribution of the corresponding flavor and @xmath453 + 1}=1-n^-(\mathbf{p})$ ] is the fermi distribution for the corresponding anti - particle .
the first line is related to particle - hole creation and annihilation , the second and third line to particle - antiparticle ( and hole - antihole ) production and annihilation , and the last line to antiparticle - antiparticlehole creation and annihilation .
note that , as expected , the first line vanishes in the @xmath454 limit and the last line vanishes in the @xmath455 limit .
the second and third line contain infinite vacuum ( @xmath456 ) contributions , which need to be renormalized .
the gluon self - energy will be decomposed in terms of a matter part , a renormalized vacuum and an infinite vacuum contribution according to @xmath457 where the various functions have been defined as @xmath458,\\ \pi_{\mu\nu}^{\textrm{inf.vac.}}&\equiv\frac{p_{\mu\nu}}{3}q^2\left[\frac{\pi_{\lambda\lambda}(q,0,0)}{q^2}\right]_{q^2\rightarrow0}.\end{aligned}\ ] ] for the vacuum expressions it is possible to extract the projection matrix @xmath459 , since the polarization tensor is purely transverse and for @xmath456 there is no preferential frame such that the tensor is built up out of only two possible quantities @xmath460 and @xmath461 . to include the dynamical properties of the vacuum expression
a renormalization procedure is necessary . in the above the renormalized vacuum expression
was obtained by extracting the infinite contributions from the vacuum polarization tensor at the renomalization point @xmath462 .
the reason for extracting the factor @xmath463 from the vacuum @xmath464 before setting @xmath462 is because in the vacuum the dressed gluons are massless , i.e. , @xmath465 , c.f .
( [ eq : ap_renormalised_vacuum ] ) .
therefore @xmath462 is still a pole for the vacuum expression , at which the residue is @xmath209 after the renormalization . in conclusion , the following full renormalized polarization tensor will be used @xmath466    in the high - density limit ( @xmath467 ) , a good approximation is to evaluate the matter part for the case of massless quarks .
thus for massless quarks in the zero - temperature limit the matter parts of @xmath468 and @xmath469 can be calculated explicitly    @xmath470\ln\left[1-\frac{2s_1s_2\mu}{q - is_1\omega_q}\right]\right.\nonumber\\ & \qquad\qquad\qquad\qquad\qquad\qquad\qquad\left.+\sum_{s_2=\pm1}\left[s_2\frac{\mu(3q^2 - 4\mu^2)}{24q}-\frac{\mu^2}{4}\frac{i\omega_q}{q}\right]\ln\left[\frac{1-\frac{2s_2\mu}{q - i\omega_q}}{1+\frac{2s_2\mu}{q+i\omega_q}}\right]\right\},\label{eq : ap_pi_mat_00}\\ \pi_{\mu\mu}^{\mathrm{mat}}(t=0)=&\frac{g^2}{\pi^2}\left[\frac{\mu^2}{2}-\frac{q^2}{8q}\left\{\mu\sum_{s_1=\pm1}s_1\ln\left[\frac{1+\frac{2s_1\mu}{q - i\omega_q}}{1-\frac{2s_1\mu}{q+i\omega_q}}\right]+\frac{1}{2}\sum_{s_1,s_2=\pm1}(q - is_1\omega_q)\ln\left[1+\frac{2s_1s_2\mu}{q - is_1\omega_q}\right]\right\}\right].\label{eq : ap_pi_mat_mumu}\end{aligned}\ ] ]    these expressions give the same result as can be found in refs .
@xcite . performing an analytic continuation to real time ( @xmath471 ) ,
the above expressions have branch cuts from @xmath472 and @xmath473 , the origins of which can be found from eq .
( [ polarizationtensor ] ) and are due to particle - hole processes and finite-@xmath474 contributions to particle - antiparticle production and annihilation , respectively . because lorentz invariance is broken due to the presence of the gas , @xmath475 is a function of @xmath476 and @xmath218 separately since it can be a function of @xmath477 and @xmath478 , where @xmath479 defines the rest frame of the system @xcite
.    the renormalized vacuum part @xmath480 has been found in ref .
@xcite @xmath481\nonumber\\ \simeq&p_{\mu\nu}\frac{g^2q^2\left[\frac{5}{3}+\ln\left(\frac{m^2}{q^2}\right)\right]}{24\pi^2}.\label{eq : ap_renormalised_vacuum}\end{aligned}\ ] ] where it was assumed that the quarks have equal masses and in the last line it was expanded for small masses . in real time
this expression has a branch cut for @xmath482 , which is due to particle - antiparticle production and annihilation .    using a sommerfeld expansion temperature corrections to the gluon self - energy
can be obtained .
the sommerfeld expansion can be summarized as @xmath483 the derivation is analogous to that of the non - relativistic sommerfeld expansion @xcite .
one has to keep in mind that the above is not a complete expansion for small temperatures since also the chemical potential depends on temperature . in this manner
the @xmath484 correction to the gluon self - energy is found to be    @xmath485,\nonumber\\ \pi_{\mu\mu}^{\mathrm{mat}}-\pi_{\mu\mu}^{\mathrm{mat}}(t=0)=&-\frac{g^2}{2}\frac{1}{6}t^2\frac{8\mu^2[3(i\omega_q)^2-(2\mu - q)(2\mu+q)]}{[(i\omega_q)^2-(2\mu - q)^2][(i\omega_q)^2-(2\mu+q)^2]}.\label{eq : ap_pi_temperature_correction}\end{aligned}\ ] ]    using eqs .
( [ fg]-[eq : ap_pi_temperature_correction ] ) the behavior of @xmath194 and @xmath195 can be found in several limits . keeping the ratio @xmath486 fixed and expanding up to zeroth order in @xmath219 , the hard dense and hard thermal loop ( hdl / htl ) expressions are re - obtained @xcite @xmath487\right|_{x = q/\omega},\\ \lim_{\omega\rightarrow0}&g(\omega , q = x\omega)\\=&m_g^2\left.\left[\frac{1}{x^2}-\frac{1}{2x}\left(1-\frac{1}{x^2}\right)\ln\left(\frac{1-x}{1+x}\right)\right]\right|_{x = q/\omega},\end{aligned}\ ] ] where @xmath488 is the gluon thermal mass .
the reason that the above limit returns the hdl and htl expressions , is due to the fact that the hdl / htl approximation ( @xmath489 and @xmath490 ) takes into account only the low - energy processes around the fermi surface , namely particle - hole processes .
some other useful limits are @xmath491 where it can be seen that in the static long - wavelength limit ( @xmath492 ) the electric gluons ( @xmath194 ) are screened , while magnetic gluons ( @xmath195 ) are dynamically screened .
furthermore , the large momenta and frequency behavior of the matter and vacuum parts are separately seen to be @xmath493}{24\pi^2}.\label{eq : ap_fg_large_pomega_behavior}\end{aligned}\ ] ] the vacuum expressions are proportional to @xmath494 since in the vacuum the gluons remain massless and thus @xmath495 is still a pole of the propagator .
the logarithmic @xmath494 dependence is usually absorbed into the vertex and subsequently interpreted as the varying of the coupling constant with the energy scale @xmath226 .
consider first the well - known ideal fermi gas term of a single species @xmath496\\ & = -\frac{1}{\beta}\sum_{\omega_n}\int\frac{\mathrm{d}^3p}{(2\pi)^3}2\ln[-(i\omega_n+\mu_0)^2+\epsilon_0 ^ 2(\mathbf{p})].\end{aligned}\ ] ] using eq .
( [ eq : ap_matsubara_fermions ] ) with @xmath497 $ ] and using the principle value logarithm with a branch cut on the negative real axis @xmath498&=\lim_{\eta\downarrow0}\im[\ln(-\omega^2+\epsilon_0 ^ 2(\mathbf{p})-2i\omega\eta)]\\ & = \pi\theta[\omega^2-\epsilon_0 ^ 2(\mathbf{p})]\mathrm{sgn}(-\omega),\end{aligned}\ ] ] the ideal fermi gas contribution becomes @xmath499\\ & \qquad\qquad\qquad\times\theta[\omega^2-\epsilon_0 ^ 2(\mathbf{p})]\mathrm{d}\omega\\ = & -2\beta^{-1}\int\frac{\mathrm{d}^3p}{(2\pi)^3}\left\{\ln\left[1+e^{-\beta(\epsilon_0(\mathbf{p})-\mu_0)}\right]\right.\\ & \qquad\left.+\ln\left[1+e^{-\beta(\epsilon_0(\mathbf{p})+\mu_0)}\right]\right\}+2\int\frac{\mathrm{d}^3p}{(2\pi)^3}\int_{\epsilon_0 ^ 2(\mathbf{p})}^{\infty}d\omega.\end{aligned}\ ] ] in the above it was assumed that the arc at infinity vanishes , which is the case if the time - ordering in the path - integral is taken into account properly by multiplying the integrand by @xmath500 and taking the limit @xmath501 .
the infinite vacuum contribution is clearly visible and should be subtracted , giving for massless fermions @xmath502.\ ] ]    the ideal bose gas follows from @xmath503,\ ] ] which will lead to the same expression in both the lorentz and coulomb gauges . in the following , the coulomb gauge is taken ( @xmath504 ) , @xmath505 - 2\mathrm{tr}\ln q^2-\mathrm{tr}\ln\frac{1}{\xi}\right)\\ & = \frac{1}{2v\beta}\sum_{\omega_n , q}\left(\ln\left[\frac{q^2}{\xi}q^2(q^2)^2\right]-2\ln q^2-\ln\frac{1}{\xi}\right)\\ & = \frac{1}{v\beta}\sum_{\omega_n , q}\ln[-(i\omega_n)^2+q^2].\end{aligned}\ ] ] the unphysical degrees of freedom clearly drop out due to the ghost contribution @xmath506 . using eq .
( [ eq : ap_matsubara_bosons ] ) with @xmath507 , the ideal bose gas contribution is @xmath508-\int\frac{\mathrm{d}^3q}{(2\pi)^3}\int_q^\infty\mathrm{d}\omega.\ ] ] again the vacuum term is clearly present and will be subtracted , giving the stefan - boltzmann law @xmath509
for completeness also the quark self - energy will be derived in the lorentz gauge to lowest order .
@xmath510.\ ] ] performing the matsubara sum gives    @xmath511}{i\omega_p+\mu-\epsilon_0(\mathbf{q})-\epsilon^g(\mathbf{p}-\mathbf{q})}[(1+n_b(\mathbf{p}-\mathbf{q}))(1-n(\mathbf{q}))+n_b(\mathbf{p}-\mathbf{q})n(\mathbf{q})]\right.\nonumber\\ & \qquad\qquad\quad\left.+\frac{\sum_s[\gamma_\mu u_s(\mathbf{q})\bar{u}_s(\mathbf{q})\gamma_\mu]}{i\omega_p+\mu-\epsilon_0(\mathbf{q})+\epsilon^g(\mathbf{p}-\mathbf{q})}[n_b(\mathbf{p}-\mathbf{q})(1-n(\mathbf{q}))+(1+n_b(\mathbf{p}-\mathbf{q}))n(\mathbf{q})]\right.\nonumber\\ & \qquad\qquad\quad\left.+\frac{\sum_s[\gamma_\mu v_s(\mathbf{q})\bar{v}_s(\mathbf{q})\gamma_\mu]}{i\omega_p+\mu+\epsilon_0(\mathbf{q})-\epsilon^g(\mathbf{p}-\mathbf{q})}[(1+n_b(\mathbf{p}-\mathbf{q}))\bar{n}(\mathbf{q})+n_b(\mathbf{p}-\mathbf{q})(1-\bar{n}(\mathbf{q}))]\right.\nonumber\\ & \qquad\qquad\quad\left.+\frac{\sum_s[\gamma_\mu v_s(\mathbf{q})\bar{v}_s(\mathbf{q})\gamma_\mu]}{i\omega_p+\mu+\epsilon_0(\mathbf{q})+\epsilon^g(\mathbf{p}-\mathbf{q})}[n_b(\mathbf{p}-\mathbf{q})\bar{n}(\mathbf{q})+(1+n_b(\mathbf{p}-\mathbf{q}))(1-\bar{n}(\mathbf{q}))]\right\},\nonumber\\ \delta\sigma(i\omega_p,\mathbf{p})=&-g^2\frac{n_g}{4n_c}\sum_{s_2,s_3=\pm1}\int\frac{\mathrm{d}^3q}{(2\pi)^3}\left[\frac{s_3}{\epsilon^g(\mathbf{p}-\mathbf{q})}\frac{\gamma_{\mu}p_e^{s_2}(\mathbf{q})\gamma_0\gamma_\mu}{i\omega_p+\mu - s_2\epsilon_0(\mathbf{q})-s_3\epsilon^g(\mathbf{p}-\mathbf{q})}\right]\nonumber\\ & \qquad\qquad\qquad\qquad\qquad\times[(1+n_b^{s_3}(\mathbf{p}-\mathbf{q}))(1-n^{s_2}(\mathbf{q}))+n_b^{s_3}(\mathbf{p}-\mathbf{q})n^{s_2}(\mathbf{q})-i_{s_2,s_3}].\label{eq : ap_quark_self_en}\end{aligned}\ ] ]    in the last line the indicator function @xmath512 is exactly the vacuum contribution such that @xmath513 is the vacuum subtracted self - energy , i.e. , @xmath514 . to first order the self - energy for a single quark and antiquark
can be found by using the free quark energy projectors and helicity projectors @xmath515.\ ] ] after projecting and evaluating the self - energy at @xmath516 the function in the first square brackets in eq .
( [ eq : ap_quark_self_en ] ) reduces to @xmath517 defined in eq .
( [ functionf ] ) . in the zero - temperature limit
the self - energy for a single quark then reduces to @xmath518 note that the final result is independent of helicity @xmath15 .