Role of Magnetic Interaction in Dense Plasma

Quasiparticle excitations and associated phenomena of energy and momentum transfer rates have been calculated in terms of the drag and the diffusion coefficients exposing clearly the dominance of the magnetic interaction over its electric counterpart. The results have been compared with the finite temperature results highlighting the similarities and dissimilarities in the two extreme regimes of temperature and density. Non-Fermi-liquid behavior of various physical quantities like neutrino mean free path and thermal relaxation time due to the inclusion of magnetic interaction has clearly been revealed. All the results presented in the current review are pertinent to the degenerate and ultradegenerate plasma.


Introduction
Understanding the properties of the hot and dense ultrarelativistic plasma has been at the forefront of contemporary research for the past few decades. Interests in this regime cover a broad area starting from the laboratory based heavy ion collisions (HIC) to the wider domain of naturally occurring astrophysical sites like neutron stars, supernovae, and white dwarfs and so forth.
In the present review we mainly focus on the properties of plasma with high chemical potential ( ) and zero or low temperature ( ). It is to be noted that the theoretical calculations so far involving both the quantum electrodynamics and quantum chromodynamics (QED/QCD) have been confined largely to the domain of high temperature particularly to address the issues related to heavy ion collisions. The main motivation for this has been the possibility of quark-hadron phase transition which might occur in HIC mimicking the conditions of microsecond old universe. Studies with dense plasma on the other hand find their applications mainly in astrophysics.
Several calculations have recently been performed exposing subtle departure of the characteristic behavior of the quasiparticle excitations in dense plasma than their high temperature counterparts. For example, the role of magnetic interaction in dense plasma, as we shall see, modifies the fermion self-energy in a nontrivial way involving fractional powers of the excitation energy [1] leading to phenomena not seen in the high temperature plasma. In fact the transverse or magnetic interaction seems to dominate over the corresponding electric or longitudinal interactions [1] in degenerate electron or quark matter. It is to be noted that such departures can only be significant in the relativistic plasma as in the nonrelativistic case magnetic interactions are = V/ suppressed [2][3][4][5]. The introduction of magnetic interaction changes some of the characteristic behaviours of dense plasma which show departure from the normal Fermi liquid case. In this review this particular aspect will be illuminated upon further.
Our discussion starts with the quasiparticle damping rate ( ). In the relativistic case ∼ ( − ) [6] unlike the known non-relativistic scenario where ∼ ( − ) 2 . Here, is the energy of the quasiparticle and is the chemical potential [6][7][8]. More importantly the infrared behavior of this quantity differs dramatically from the high temperature case [6][7][8].
In a plasma the other quantities of interest have been the 2 Advances in High Energy Physics drag ( ) and the diffusion coefficients (B) which eventually determine the equilibration time scale of the plasma. These two quantities are actually related to the energy loss and the momentum transfer due to the scattering of the constituent particles of the plasma. In high temperature case they are known to be connected via the Einstein relation. We address this issue to see the scaling behavior of these quantities in the domain of zero and high temperature plasma [2].
Recently the implications of the medium modified quasiparticle self-energies have been discussed by several authors showing their significance in the cooling behavior of neutron stars. The appearance of logarithmic terms in the specific heat, emissivity, and thermal relaxation time of the degenerate matter has been found to be of crucial importance. How the medium modified dispersion relation changes the cooling behavior of neutron star has been first exposed in [9] and later in [3,5,10].
The plan of the review is as follows. In Section 2 we have shown the behavior of quasiparticle damping rate in ultradegenerate plasma. In this section we compare the results of quasiparticle lifetime in the extreme limits of high temperature and high density. For the former it is known that the damping rate for the exchange of static gauge boson remains unscreened even when one uses hard thermal loop corrected propagator for the gauge bosons [11,12]. It has been shown in the literatures that a further resummation is required to resolve this issue by invoking Bloch-Nordseick techniques as discussed in great detail in [13,14]. For the degenerate case, as we shall see this problem does not appear. In this context however we shall expose the relative importance of the magnetic interaction over the electric part unlike the hot relativistic plasma.
In this review we also calculate the drag and the diffusion coefficients for both degenerate and ultradegenerate cases. Different properties of and B in the zero temperature and low temperature limit have been reported in Section 3. In Section 4 explicit analysis has been presented to exhibit how the nature of neutrino mean free path changes with the inclusion of the transverse mode. One distinguishing feature of the present work is to show the structure of the fermion self-energy near the Fermi surface. Of particular importance is the vanishing group velocity at the Fermi surface which in turn is responsible for the non-Fermi liquid (NFL) corrections. Section 5 describes the implication of the neutrino mean free path vis-a-vis cooling of the neutron star specially focusing on NFL corrections to the various physical quantities. We calculate the thermal relaxation time exposing the role of magnetic interaction further in Section 6 where we also discuss NFL corrections.
In the last two sections of this review we estimate the neutrino mean free path and thermal relaxation time of degenerate matter. At very high density nowadays it is well known that quark matter is expected to form a color superconducting color-flavour locked (CFL) phase [15]. In this phase all quark excitations are gapped, and the mean free path or the thermal relaxation time is modified by exponential factor involving gap parameter [16][17][18]. However, in order to highlight particular aspects of NFL behaviour of the relativistic plasma we deal with only ungapped quark matter. The extension of the present calculation for the gapped quark matter is straightforward.

Damping Rate
The propagation of a particle immersed in plasma gets modified through interaction with the medium. In case of ultradegenerate ( = 0, ̸ = 0) plasma these medium modified particles (quasiparticles) scatter with the other particles close to the Fermi surface. Information about the quasiparticle lifetime is obtained from the retarded propagator. The damping rate is expected to follow exponential decay in time ( , ) ∼ − − , so that | ( , )| 2 ∼ −Γ , where Γ = /2. The life time of the single-particle excitation is then ( ) = 1/Γ( ). The concept of quasiparticles becomes meaningful only if their life time is long or in other words the damping width is much smaller than its energy. In a plasma the damping rate of the quasiparticle can be written in terms of the imaginary part of the fermion self-energy or the scattering rate (see Figure 1). Mathematically one writes [19] where Σ is the fermion self-energy, → 0 + , and we have used = ( 0 , p). The aforementioned expression is true for both the high temperature and the zero temperature plasmas. The one-loop electron self-energy is given by In the previous equation ( − ) is the free fermion propagator and Δ ] ( ) is the photon propagator. Focusing on (2) we can see that it involves the boson propagator (Δ ] ) which is well known to diverge in the infrared regime. This is a generic feature of both the high temperature and the zero temperature plasmas. In case of hot QED or QCD plasmas one handles this problem by employing the hard thermal loop (HTL) prescription which was originally developed in [11,12]. In this method the exchanged momentum (q) integration is performed in two domains by introducing an intermediate cut-off parameter ( * ), below which the one-loop corrected dressed boson propagator has to be used and the bare propagator can be used to perform the integration above * (Braaten and Yuan's prescription) [20].
For electric interaction this is sufficient to remove the infrared singularity associated with the exchange of massless Advances in High Energy Physics 3 bosons (photons or gluons). The interaction mediated by the transverse bosons on the other hand still poses a problem at high temperature [13,14]. We shall discuss this further at the end of this section where we compare the zero temperature case with the corresponding high temperature calculations. It might be mentioned that to construct the one-loop corrected photon propagator one needs to evaluate the photon selfenergy in dense plasma [21]. In the ultradegenerate case, the dominant contribution to the photon self-energy comes from the loop momentum ∼ , and the external momentum for the soft modes is assumed to be q ∼ ( Figure 2). This, in principle, runs similar to the HTL approximation where the loop momentum is ∼ and the external momentum is ∼ . Likewise for the degenerate case one can construct hard dense loop (HDL) propagator where loop momentum is assumed to be ∼ [21].
The structure of the dressed photon propagator in the Coulomb gauge can be written as [19] with = ( −̂̂),̂= q /|q|, 0 = 0 = 00 = 0, and Δ , Δ are given by [19] Δ , are the longitudinal and the transverse components of the boson propagator, and Π , are the longitudinal and transverse parts of the self-energy [19]: where the Debye mass is = / . The bosonic spectral functions are the imaginary part of the propagators: where , in the HDL approximations are given by [6] ( 0 , ) To calculate in the ultradegenerate plasma the energy of the quasiparticle is therefore considered to be hard, that is, ∼ . Hence, for the quasiparticle with momentum close to the Fermi momentum the electron-photon vertex can be replaced by the bare one and the one-loop self-energy is dominated by the photon with soft momentum (q ∼ ). So, the calculation of imaginary part of Σ can be performed with the help of the free fermion spectral function ( − ) = 2 ( 0 ) ( 2 0 − 2 ) ( ( 0 ) = 0 /| 0 |) and the medium modified , [6]. Explicitly, In the previous equation ( 0 ) and ( 0 ) are the boson and fermion distribution functions. In case of the cold plasma these distribution functions can be replaced by (1 + ( 0 )) = Θ( 0 ) and ( 0 ) = Θ( − + 0 ), where Θ represents the step function. These theta functions, as we shall see, restrict the 4 Advances in High Energy Physics phase space of the 0 integration severely. This eventually with the help of the HDL propagator makes finite [13,14,21]. To proceed further an approximation in the region ( − )/ = ] ≪ 1 has been made. This is due to the fact that scatterings close to the Fermi surface are only relevant as mentioned earlier. To perform the integration the integration domain is divided into two sectors. In the soft sector ( < * ) with the HDL corrected spectral function mentioned in (8) one obtains [6] Only the leading order terms in ( − )/ have been retained in writing the previous equation. For the hard momentum region > * bare propagator is sufficient to calculate the damping rate: where max ≃ is the maximum momentum transfer that is allowed by kinematics. Two points are to be noted here. One that * dependent term appear only at the higher orders that is, the leading order result is independent of * . Second upon summation the * terms cancel out exactly. The total damping rate is then given by [6] where G 1 = ( 2 )/(24 ) and G 2 = ( 2 )/32. In the last equation the first term corresponds to the magnetic interaction, and the second term comes from the electric interaction revealing the fact that the magnetic interaction here dominates over the electric one. In this connection we can recall the results of the damping rate in case of the high temperature plasma with vanishing chemical potential. These are given by [13,14] where is the plasma frequency related to the Debye mass ( = / √ 3). The last equations have been obtained by using the oneloop corrected dressed propagator which makes finite. The transverse part, on the other hand, remains divergent that can be cured by another resummation as discussed in [13,14]. In the case of dense QED or QCD matter interestingly the second resummation is not required. This happens, as indicated above, due to the restrictions imposed by the Pauli blocking which severely cuts off the phase space. A comparative study between high and zero temperature damping rates has been listed in Table 1 [22].
In the nonrelativistic case the excitation energy dependence of can be derived from the simple phase-space factor without performing any rigorous calculation. A particle with energy > interacts with the particle of energy < .
The scattered particles are now in the energy state > and > . The total probability of the scattering process is proportional to ∫ ( + − − ) 3 3 . For − ≪ , the permissible regions of variation for the scattered particles are < < + − and 2 − < < . The integral over 3 and 3 can be computed with the approximation ≈ ≈ , that is, scattering occurs near the Fermi surface. The scattering probability is then ( − ) 2 [23]. Similar arguments can be made for the transverse sector which will render the damping rate of the particle ∼ ( − ) at the leading order for relativistic case.

Drag and Diffusion Coefficients
3.1. Zero Temperature Plasma. Next we review the drag and diffusion coefficients in zero temperature relativistic plasma and see how the magnetic interaction changes the behavior of these coefficients. These two are the important quantities to study the equilibration properties of the plasma as mentioned earlier and are related to the energy and the momentum transfer per scattering, respectively. can thus be defined as follows [24]: The energy loss −( / ) can be calculated by averaging over the interaction rate (Γ) times the energy transfer per scattering [25]: where Γ = /2. Inserting the energy exchange in the expression of to calculate / from the previous equation we obtain ( Figure 2) The basic assumption behind the calculation of lies in the fact that ( − )/ = ] ≪ 1. Taking into consideration the previous assumption in case of the hard photon exchange we get [2], The above expression is infrared divergent with algebraic dependence. Here, it would be worthwhile to comment that in case of the high temperature plasma the divergence is logarithmic ( 4 2 /24 ) ∫( / ). The divergence is even worse in case of the ultradegenerate plasma. But we show below that with the HDL propagator and Pauli blocking one Advances in High Energy Physics 5 Table 1: Comparison of high and zero temperature damping rates.
in high temperature plasma in zero temperature plasma Infrared behaviour IR divergent after the HTL resummation Finite with HDL resummation Dominance and contribute at the same order ( ( 2 )) Dominant contribution from can render finite. In the soft sector is given by the following expression: The integration domain ( ) above is restricted by the Θ functions: After explicit calculation, the electric and the magnetic contributions to the expressions of take the following form [2]: In the last two equations the leading order terms are independent of the cut-off parameter as they appear only at ( 7 ). In the region of hard photon momentum exchange we obtain [2] On addition of the soft and the hard sectors it can be seen that the expression of the drag coefficient becomes [2], where H 1 = ( 2 2 )/(48 ) and H 2 = ( 2 2 )/(96 ).
In the previous equation the first term corresponds to the magnetic interaction, and the second term has come from the longitudinal interaction. From (21) it is important to note that the hard sector involves even higher power. The entire contribution at the leading order then comes from the soft sector alone, and Braaten and Yuan's prescription is not required in case of the ultradegenerate plasma [2]. The final expression of reveals that with the magnetic interaction changes significantly from the Fermi liquid result. In the Fermi liquid theory, since the magnetic interaction is suppressed in comparison with the electric sector goes with ( − ) 3 . Here, we can recall the result of the drag coefficient in case of high temperature plasma and compare the same with the previous result. The coefficient at high temperature looks like [25][26][27][28][29]

36
(− 1 2 From the previous result it is clear that there is no splitting between electric and magnetic modes, and in hot plasma both contribute at the same order. A comparative study between high temperature and zero temperature can be found in [22]: Another important quantity to study the equilibration property of the plasma is the momentum diffusion coefficient . The coefficient is related to the interaction rate via the following relation [24]: Decomposing into longitudinal ( ) and transverse components ( ) we get the following expression, , are the longitudinal, transverse squared momenta acquired by the particle through collision with the plasma. The longitudinal momentum diffusion coefficient is related to the drag coefficient via the Einstein relation in case of equilibrating plasma. Like longitudinal momentum diffusion coefficient ( = B, suppressing the index ) can be expressed as Here, ‖ is the longitudinal exchanged momentum transfer in a scattering and ‖ = cos . Explicitly B is then given by [2] B where B 1 = ( 2 3 )/(72 ) and B 2 = ( 2 3 )/128. In this case also the electric and magnetic modes split as the magnetic mode is ∼( − ) 3 whereas electric mode is ∼ ( − ) 4 . One can relate the leading order terms of and B from (22) and (27) via a common relation, namely, Einstein's relation as in case of high temperature plasma. In ultradegenerate plasma the relation can be written as B = 3 ( − ) /4. In the high temperature and zero chemical potential regime B has the following expression [25][26][27][28][29]: In case of high temperature plasma from (23) and (28), Einstein's relation becomes B = 2 . From (22), (23), (27), and (28) it can be seen that zero temperature and B show entirely different characteristics from the high temperature case. The logarithmic dependence in high temperature takes an algebraic form in the ultradegenerate limit [22].

Low Temperature Plasma.
In this section we extend our calculations to include finite temperature effects in the domain where | − | ∼ ≪ ≪ . Moreover, we also go beyond the leading order calculations as previously reported above. Finally we shall see how the leading order zero temperature results appear as a limiting case from our present analysis including higher order correction terms.
For the soft sector contribution we take / as an expansion parameter and = ( − )/ ∼ 1. We start with the following expression for the soft sector [4]: The previous equation can easily be derived from (16). After subtracting the energy independent part from the above equation we get, We substitute and 0 by introducing dimensionless variables and V in the above equation: where is the screening length in the magnetic sector and given by ( 2 0 /2) 1/3 . This fractional power in the transverse sector will eventually show up in different coefficients as we show below. The previous substitutions actually yields, After performing the 0 and integrations we obtain [4]: where In case of the longitudinal sector we substitute = and 0 = / or = and = / . Since screening length is different in electric and magnetic sectors the substitutions therefore involve different coefficients of Advances in High Energy Physics 7 and for the transverse and the longitudinal cases [1] as can be seen from the structure of , (8). For the leading term in the longitudinal sector one writes, The final expression for drag coefficient then becomes [4] From the previous expression it is evident that is polylogarithmic in nature and contains fractional powers in ( − ).
The fractional powers indicate the nonanalytic behavior of . Significant it is to note that here the subleading transverse term is greater than the leading order longitudinal one. Thus we see that here too the transverse contribution dominates over the longitudinal one. In the zero temperature limit the functions behave as ℎ ( ) → | | and ( ) → | |. Hence, in the limit becomes [4] In the last equation H 3 = (3 × 2 1/3 2 2 )/(72 7/3 ), and H 1 , H 2 are the same as before. We can infer from the previous observations that the inclusion of transverse interaction for a relativistic particle changes the nature of . In the nonrelativistic case considering only the electric interaction goes as ( − ) 3 , but with the inclusion of the transverse sector we see that anomalous fractional powers appear. Like one can also compute the longitudinal momentum diffusion coefficient in the low temperature limit. The expression for the longitudinal momentum diffusion coefficient as defined in (26) can be written as [4] where The final expression for B in the extreme zero temperature limit becomes [4] It is to be noted that, for the low temperature case, so far, no reference has been made to the case of hard momentum regime. This is justified as we have seen before that up to the order of our interest the hard sector does not contribute, and the entire contribution comes from the soft sector alone. Furthermore, it is worthwhile to mention that the NFL behaviour of these coefficients is related to the dynamical screening involving " " in the static limit when one evaluates Π . To see this one can recall the expressions for the polarization functions in (6). In the static limit ( 0 / ∼ 0) the functions become The " " that appears in the denominator of (5) from the Π above along with the Pauli blocking is responsible for this NFL behavior. In Figures 3 and 4 and B versus energy of the incoming fermion in the small temperature ( / ≪ 1) region ( = / is the Fermi temperature) have been plotted. From the figures, it is evident that with increasing / , both and B decrease. This nature is consistent with what one finds for the fermionic damping rate at small temperature [1].

Neutrino Mean Free Path
In this section we calculate the neutrino mean free path (MFP) for cold and warm QCD matters. When a new star is born following a supernova explosion, it starts emitting neutrinos via the direct or the modified Urca processes. Emission of these neutrinos is responsible for the initial cooling of the star so produced. For quark matter, the dominant contribution to the emission of neutrinos is given by the decay and the electron capture [9]. These reactions are named as "quark direct Urca" processes which have been studied in detail by Iwamoto [30,31]. For quark matter, the MFP was previously derived in [30] where the calculations were restricted to the leading order by assuming free Fermi gas interactions. In this context we recall two important works on the neutrino MFP in QED plasma. One is due to Tubbs and Schramm [32], and the other is done by Lamb and Pethick [33]. In [32], MFP was calculated in the neutranized core and just outside the core. On the other hand, it has been shown in [33] that neutrino degeneracy reduces the neutrino MFP which suggests that neutrino may flow out of the core rather slowly. Here, we revisit the calculation and go beyond the leading order results to incorporate the NFL corrections both for degenerate and nondegenerate neutrinos.
In the interior of a neutron star, there are two distinct phenomena for which the neutrino mean free path is calculated: one is absorption, and the other involves scattering of neutrinos [31]. The corresponding mean free paths are denoted as abs mean and scatt mean , and combine them to obtain the total mean free path [34]: In our model the Lagrangian density is described by [30] L ( ) = √ 2 ( ) J ( ) + H.C., where the weak coupling constant is ≃ 1.166 × 10 −11 in MeV units, and and J are the lepton and hadron charged weak currents, respectively. The weak currents are where is the Cabibbo angle (cos 2 ≃ 0.948) [35]. The mean free path is determined by the quark neutrino interaction in dense quark matter via weak processes. We consider the simplest decay reactions: the absorption process and the other is its inverse relation The neutrino mean free path is related to the total interaction rate due to neutrino emission averaged over the initial quark spins and summed over the final state phase space and spins. It is given by [31] where is the spin and color degeneracy, which in the present case is considered to be 6. Here, , , and are the energy, momentum, and distribution function for the corresponding particle. | | 2 is the squared invariant amplitude averaged over initial quark spin and summed over final spins of quark and electron as given by [31] | | 2 = 1 2 ∑ , , Here, we work with the two flavor system as the interaction involving strange quark is Cabibbo suppressed [9,36].

Degenerate Neutrinos.
In this section we consider the case of degenerate neutrinos. When the neutrino chemical potential ( ] ) is considered to be much larger than the temperature ( ), the neutrinos become degenerate. This is also termed as trapped neutrino matter. In this case both the direct (46) and inverse (47) processes can occur, and both terms in (48) under curly brackets are retained [31]. Consequently, the equilibrium condition becomes + ] = + . Neglecting the quark-quark interactions and by using (48) and (49), the mean free path becomes 1 abs, mean In the square bracket, the second term − ( ] − ] ) is due to the inverse process (47). One can neglect the mass effect on the neutrino MFP as the masses of , quark and electron are very small. To perform the momentum integration we define ≡ | + ] | = | + | as a variable. Following the same procedure as described by Iwamoto [31] one has where we denote the single particle energy ( ) as . For the free case / is the inverse quark velocity. It is well known that this slope of the dispersion relation changes in matter due to the scattering from the Fermi surface and excitation of the Dirac vacuum. The modified dispersion relation can be obtained by computing the on-shell one-loop self-energy. For quasiparticles with momenta close to the Fermi momentum, the one-loop self-energy is dominated by the soft gluon exchanges [8]. The quasiparticle energy satisfies the relation [8,21] = ( ) − Re Σ ( , ( )) , where we have retained only the real part of self-energy since the imaginary part of Σ turns out to be negligible compared to its real part [1,37]. For the detailed analysis we refer the reader to [11,12,21]. The analytical expressions for one-loop quark self-energy can be written as [1,8,[37][38][39][40] It exhibits a logarithmic singularity close to the Fermi surface, that is, when → . ( )/ can be obtained by differentiating the dispersion relation (55) with respect to . At leading order in / this yields where is the strong coupling constant, = ( 2 −1)/(2 ), and is the color factor. Using (57), (50), and (52) To derive (58), we use the chemical equilibrium condition ( ) + ( ) = ( ) + (]) neglecting the masses for quarks and electrons. For the phase-space integral we have [31,41] 10
The other contribution to the mean free path comes from the quark-neutrino scattering. The neutrino scattering process from degenerate quarks is given by for each quark component of flavor (= or ). The scattering mean free path of the neutrinos in degenerate case can be calculated similarly as evaluated by Lamb and Pethick in [33] for electron-neutrino scattering. Assuming / ≪ 1 and including the suitable modifications of the phase space, the mean free path is given by Here, is the quark mass. and are the vector and axial vector coupling constant given in Table II of [31]. The expression of (61) has been found to be in agreement with the results reported in [33] for dense and cold QED plasmas by making suitable changes for the color factors and by dropping the second square bracketed term. In (61) the constants 0 ≡ 4 2 2 / [32], and is the quark number density: where 2 is the quark spin degeneracy factor. Explicit form of Λ( ) can be written as [31,33] Λ ( ) = 4 3 Min ( ] , ) and = ] / if ] < , and = / ] if ] > .

Nondegenerate Neutrinos.
The mean free path has been also derived for nondegenerate neutrinos, that is, when ] ≪ . This is the case for untrapped neutrino matter. For this case the inverse process (47) does not contribute to the MFP. Hence we neglect the second term in the curly braces of (48). Here, only those fermions take part in the reaction whose momenta lie close to their respective Fermi surfaces. It is to be mentioned here, if quarks are treated as free, as discussed in [31,42,43], the matrix element vanishes since , quarks and electrons are collinear in momenta. The inclusion of strong interactions between quarks relaxes these kinematic restrictions resulting in a nonvanishing squared matrix amplitude. Since the neutrinos are produced thermally, we neglect the neutrino momentum in energymomentum conservation relation [31]. This is not the case for degenerate neutrinos where ] ≫ , and therefore such approximation is not valid there. By doing angular average over the direction of the outgoing neutrino, from (49) the squared matrix element is given by [9]: where = (1, V ) is the four velocity. To calculate (64) we have used the chemical equilibrium condition = + and also the relations derived from Fermi liquid theory [9] V = 1 − 2 ; = . (65) Putting | | 2 in (48) we have Neglecting the neutrino momentum in the neutrino momentum conserving function, the integrals can be decoupled into two parts. Following the procedure described by Iwamoto [31], the angular integral is given by Changing the variables to = ( − ) , = − ( − ) , and = − ( − ) and denoting the single particle energy ( ) as we have from (68) As the contribution dominates near the Fermi surfaces, extension of lower limit is reasonable [41,44]. The integration of (69) can be performed using (57) and following the procedure defined in [37,41,44] to give Using (66), (67), and (70), the mean free path at leading order in / is given by Advances in High Energy Physics 11 The first term is known from [31], and the additional terms are higher order corrections to the previous results derived in the present work. For the scattering of nondegenerate neutrinos in quark matter, the expression of mean free path was given by Iwamoto [31]. We incorporate the anomalous effect which enters through phase space modification giving rise to 1 scatt, mean Here, we have assumed / ≪ 1. The constants 0 and number density have been defined earlier. We now estimate the numerical values of the neutrino mean free paths. Here, ] is set to be equal to 3 and = 10 MeV [31,34]. For the quark chemical potential, following [9], we take ≃ 500 MeV corresponding to densities ≈ 6 0 , where 0 is the nuclear matter saturation density. The electron chemical potential is determined by using the charge neutrality and beta equilibrium conditions which yields = 11 Mev. The other parameters used are the same as [9].
From Figure 5, we find for degenerate neutrinos at the anomalous logarithmic terms reduce the value of the mean free path appreciably both in the low and high temperature regimes. Figure 6 shows that for nondegenerate neutrinos NFL correction is quite large at low temperature while at higher temperature it tends to merge. It is interesting to see from these two plots that NFL correction to the MFP in degenerate neutrinos is less than that of nondegenerate neutrinos [3]. This reduced mean free path is expected to influence the cooling of the compact stars.

Cooling of Neutron Star
To analyze the cooling of the compact star [45][46][47], the specific heat capacity of the quark matter core needs to be taken into consideration along with the emissivity via the cooling equation [9,36] where is the internal energy, V is the specific heat at constant volume and baryon number, is the emissivity, is time, and is the temperature, and we have assumed that there is no surface emission. Emissivity under different scenarios both for quark and nuclear matters has been studied extensively. For the specific heat however earlier literatures were concerned only with the kinetic energy contribution without paying much attention to the QCD interaction. Recently several authors have shown that the interaction can lead to a significant correction to V . It has been shown that this interaction mediated by the gluons can give rise to a leading ln −1 term in the specific heat [48,49]. This logarithmic enhancement comes from longrange chromomagnetic fields which eventually leads to non-Fermi-liquid behavior. The calculation of V including NFL  correction has been done in [48,49] which we quote here for further use: On the other hand, the neutrino emissivity of degenerate quark matter was first derived by Iwamoto [31]. Latter Schäfer and Schwenzer [9] calculated the neutrino emissivity including the corrections for NFL effects. The calculation has been made by considering the energy loss due to neutrino emission from the core of the neutron star. Alternatively, one can also determine the neutrino emissivity from the mean free path by evaluating the integral [ (75) This is the same as given in [9]. In this context, it is to be mentioned that without such non-Fermi liquid effects we have ∼ 6 and V ∼ . In this case the temperature scales turn out to be ∝ 1/ 1/4 [9]. With logarithmic corrections included, no simple analytic solution is possible, and we solve (73) numerically. We observe from Figure 7 that the cooling of neutron star is significantly faster with NFL corrections as compared to the Fermi liquid result [9].

Thermal Relaxation Time
In this section we address the issue of thermal conduction that follows the initial cooling by the neutrino emission. This phenomenon is intimately connected to the thermal relaxation mechanism. A new born star emits a large amount of neutrinos to cool the core, but the crust remains hot. Hence, a temperature gradient is set up between the crust and the core. Immediately after the neutrino emission discussed earlier in the paper thermal energy gradually flows from outer crust to the inner core by heat conduction which, alternatively, might be viewed as the propagation of the cooling waves from the center towards the surface leading to thermalization. We in this section estimate the thermalization time scale ( ) and reveal the NFL behavior of the quantity in the context of the neutron star crust. In the latter the electrons constitute an almost ideal Fermi gas and scatter between themselves. The kinetic theory definition of is In case of the strongly degenerate electron gas the electron thermal conductivity ( ) can be expressed in terms of the thermal current as follows: where ] is the total effective collision frequency. In the neutron star crust the main components that contribute to different transport coefficients are electrons and ions. Hence, the total collision frequency is the sum of the partial collision frequencies of electron-ion (] ) and electron-electron (] ). Evidently is related to and via the following expression [50,51]: Thus the derivation of involves the calculation of ] . To derive ] we start with the Boltzmann equation which describes the kinetics of the particles: where p is the momentum of the quasiparticle, F is the external force, k p is the velocity of the heat carrier, and is the distribution function of electrons. The right hand side of the Boltzmann equation is the collision integral describing the rate of scattering of fermions. In the presence of a weak stationary temperature gradient Boltzmann equation takes the following form: The Fermi-Dirac distribution functions deviate from equilibrium distribution functions due to the presence of a weak temperature gradient, which we write as is the particle energy, is the chemical potential, and is the temperature. The term with Φ measures the deviation from equilibrium. The collision integral can be written as follows:  In the previous equation ], ] are the degeneracy factors and | | 2 is the squared matrix element for the 2 → 2 scattering process. The ± sign is meant for stimulated emission or Pauli blocking. Since we consider only the electron-electron scattering, the negative sign will remain in the phase space factor. We write the previous equation in the form | ⟩ = |Φ⟩, where | ⟩ = ( − )V and is the integral operator. The thermal conductivity is then given by the maximum of the following equation [52,53]: ⟨⋅ | ⋅⟩ denotes an inner product; for Ψ = Φ the quantity ⟨ |Ψ⟩ 2 / ⟨Ψ| |Ψ⟩ is minimal with the minimal value . With the help of this definition one can write the first bracketed term in the denominator is the thermal current . In principle Φ should be determined by minimizing (84) using the variational principle. But in the present scenario we consider the simplest trial function as [50,52] Ψ ∝ ( − ) V .
The term in the bracket can be averaged over the axis keeping and fixed using the previous trial function, where, is the azimuthal angle between V and V . After averaging we obtain [52] For small energy transfer the electron-electron scattering squared matrix element is given by [52] | | 2 = 32 4 [ 1 The medium modified photon propagator in the previous equation contains the polarization functions Π ( 0 , ) and Π ( 0 , ), which describe plasma screening of interparticle interaction by longitudinal and transverse plasma perturbations, respectively. Now, we first analyze the denominator of (85). The denominator is the thermal current as defined earlier and is given by Here, the degeneracy factor for electrons is ] = 2.
The fermionic dispersion relation gets modified due to the inclusion of the fermion self-energy in the presence of the medium. Hence, energy momentum relation changes from the vacuum. For the momentum integration in the phasespace integral in (85) we will use the medium modified dispersion relation given in (55). For this one needs to know the fermion self-energy. The fermion self-energy for the QCD matter has already been quoted in (56). In case of electrons at low temperature with NLO correction it is given by the following [1]: where 2 = 2 /2 and is chosen to be ( − ). The phasespace correction due to the medium modified dispersion relation can now be written as [10] ≃ 1 + where = / . The expression for −1 now becomes where P = 3 4 /4 5 and In the previous equation we have used the following expression [10,52]: while writing the previous equation we have considered the medium modified phase-space factor only for the particle with momentum . One can in principle take the correction for both the particles with momentum and , but that will contribute at higher order in coupling constant. The final expression for the electron thermal conductivity now takes the following form [10]: Unlike the Fermi-liquid result where ∝ 1/ , here the temperature dependence is nonanalytical and anomalous in nature. This is reminiscent of other quantities involving ultradegenerate plasma presented in the previous sections [1,[3][4][5].
involves fractional powers in ( / ) coming from the medium modified phase space factor too.
For the estimation of relaxation time the other quantity which we require is the specific heat. For the degenerate electron gas it can be written as [48,49] With (95) The thermal relaxation time up to NLO terms contains some anomalous fractional powers originated from the transverse interaction. This in turn changes the temperature dependence of nontrivially. The departure from the Fermi liquid result ( ∝ 1/ 2 ) is evident from the final expression of . In Figures 8, 9, 10, and 11 we have plotted and with using (95) and (97). From the plots it can be seen that the inclusion of both the medium modified propagator and decreases the value of and . It shows strong deviation from the Fermi liquid in both cases. In Figures 9  and 11 it has been shown that reduces both the thermal conductivity and the thermal relaxation time. This has a serious implication on the total electron conductivity . In [50] it has been shown that magnetic interaction decreases which in turn increases the electron-electron collision frequency. Thus to the total electron thermal conductivity electronelectron scattering dominates over electron-ion scattering. The phase space correction due to the medium modified electron dispersion relation further enhances the electronelectron collision frequency [10].

Conclusion
Investigation of the quasiparticle excitation in ultradegenerate plasma has been the cardinal focus of the present review. The main interesting feature which has been exposed here is the role of magnetic interaction mediated by the transverse gauge bosons leading to the phenomena very different from its high temperature counterpart. In particular the NFL behaviour of various quantities like V , , and total mean involving excitations near the Fermi surface is a very special characteristic behaviour of degenerate or ultradegenerate plasma. We also show clearly how the dynamical screening leads to finite damping rate emanating from the restricted phase space driven by the Pauli blocking. Besides the damping rate we also present results for the drag and the diffusion coefficients, that is, the energy and the momentum relaxation of the quasiparticle excitation in such a plasma revealing further the domination of the magnetic interaction over the electric interactions in ultradegenerate plasma. This separation between the longitudinal and the transverse sectors is a very special feature of the cold and dense plasma which is not seen in the finite temperature case. Another aspect of the present work has been the modification of the quasiparticle dispersion relation which changes the phase-space in a nontrivial way leading to modification of and mean . Further investigations along this line are pertinent to understand various phenomena of astrophysical interests. Some of these extensions have recently been undertaken by various authors showing the importance of this emerging area particularly due to the possibility of experiments in the domain of high chemical potential and low temperature in addition to its relevance in the astrophysical context.