by far the study of neutron stars has been mainly focused on the relationship between the equation of state ( eos ) of nuclear matter and the observed maximum mass .
the connection has been achieved by solving the hydrostatic equilibrium equations based on general relativity .
the first generation observed masses exhibited an average value around 1.5 solar masses .
this value requires a soft eos that can be easily obtained by introducing new degrees of freedom like hyperons , kaons or quarks accompanied or not by a phase transition .
sometimes the softening was so large that the neutron star is predicted to collapse into a black hole , as for the sn1987a @xcite . in the new generation of observations
the masses are distributed within a large range , up to 2 solar masses , that requires a stiff eos , i.e. , hadronic matter without new degrees of freedom . since it is hard to imagine pure hadronic matter to sustain the high pressure predicted in the inner core , new scenarios have to be advanced to explain the coexistence , in the phenomenology of neutron stars , of low and high mass spectra .
recently @xcite it has been argued that the two observed classes of neutron stars might correspond to two different evolutionary scenarios of neutron stars . in one case , the hot and dense remnant of the supernova explosion rapidly evolves into a hybrid star , where the transition to a quark phase softens the nuclear matter so that @xmath0 ; in the other case a slow evolution could lead the neutron star to a large mass via a mass accreting from the coupling with a white dwarf . from this point of view the destiny of the remnant
is strongly affected by the initial conditions , i.e. density , temperature , leptonization degree etc .
for instance , if the mass of the remnant is below the mass threshold for quark nucleation the transition to the quark phase is forbidden @xcite . if the mass is slowly accreting the transition is allowed .
the role of temperature or other parameters defining the initial state of a new born neutron star has not yet been studied .    to investigate the possible phase transition to quark matter in neutron stars
, we need also to know the eos of quark matter .
although we have in hand the fundamental theory of strong interactions , i.e. , quantum chromodynamics ( qcd ) , we still do not know the true ground state .
it is now generally expected that quark matter is in the color - flavor locked phase ( cfl ) @xcite at extremely high densities when the finite current mass of strange quarks becomes unimportant . in the density range from nuclear saturation to cfl , there may exist a rich and varied landscape of phases , e.g. , the 2sc , g2sc , gcfl etc .
presently , however , these phases suffer from the so - called chromomagnetic instability problem for both the two- @xcite and three - flavor @xcite cases . on the other hand , experiments show that quarks become asymptotically free rather slowly @xcite .
therefore , in the present study we are dealing with the ordinary strange quark matter ( sqm ) @xcite .    the special problem in studying the eos of ordinary quark matter is to treat quark confinement in a proper way . in the conventional standard approach , an extra constant term
, the famous bag constant @xmath1 , is added to the energy density of the system , which provides a negative pressure to confine quarks within a finite volume , usually called a ` bag ' .
the quark mass is infinitely large outside the bag , and a finite constant within the bag
. a vast quantity of investigation have been performed within the framework of the bag model @xcite .    as is well known , however , particle masses vary with environment .
such masses are usually called effective masses .
effective masses of hadrons and quarks have been extensively discussed , e.g. , within the nambu - jona - lasinio model @xcite and within a quasi - particle model @xcite . in principle
, the density dependence of quark masses should be connected to the in - medium chiral condensates @xcite .
taking advantage of the density dependence , one can describe quark confinement without using the bag constant .
instead , the quark confinement is achieved by the density dependence of the quark masses derived from in - medium chiral condensates @xcite .
the two most important aspects in this model are the quark mass scaling  @xcite and the thermodynamic treatment  @xcite .
both aspects will be reviewed in this paper .    in the present contribution
, the transition from hadron phase ( hp ) to strange quark phase ( sqp ) in the inner core of a neutron star is investigated within the fully consistent nuclear and quark models . in the hadron sector
we adopt the equation of state from brueckner - bethe - goldstone ( bbg ) approach with three - body forces ( tbf ) @xcite .
this theory , being a completely microscopic approach , can easily incorporate degrees of freedom such as nucleon resonances [ @xmath2(1232 ) or @xmath3(1440 ) ] , which are expected to appear at higher hadron densities . it is found that the mixed hadron - quark phase can occur , for reasonable values of the confinement parameter , a little above the normal saturation density , and can undergo the transition to pure quark matter at about 5 - 6 times the saturation .
this result is quite different from the previous results from nambu - jona - lasinio ( njl ) model in which the mixed quark phase can not appear at neutron - star densities  @xcite .
afterwards , the influence of the mixed and quark phases on the structure of compact stars is discussed by solving the tolman - oppenheimer - volkov ( tov ) equation and extracting the mass - radius plots for neutron stars .
finally , it is shown that the transition to the deconfined phase turns out to be incompatible with the onset of kaon condensation .
sqm has been one of the hot topics in nuclear physics since witten s famous stability conjecture @xcite . in many studies ,
the quark confinement was treated adopting the bag mechanism @xcite .
an alternative approach to obtain confinement is based on the density dependence of quark masses @xcite .
this mechanism has been extensively applied to investigating the properties of sqm @xcite . in this section ,
we first give a short review of the two most important aspects , with the main inconsistency of the original model .
then we present a fully self - consistent thermodynamic treatment .
the properties of sqm will be given in the new treatment .
in the present paper , however , the main application of the new approach is the study of the phase transition in compact stars .
as mentioned above , the quark confinement in this model is achieved by the density dependence of quark masses .
therefore , the first important question is how to determine the quark mass scaling which can reasonably produce confinement .
originally , the interaction part of the quark masses was assumed to be inversely proportional to the density @xcite .
this linear scaling has been extensively applied to studying the properties of sqm .
there are also other mass scalings @xcite .
their main drawback is that they are pure parametrizations without any convincing derivation .
therefore , a cubic scaling was derived based on the in - medium chiral condensates and linear confinement at both zero @xcite and finite temperature @xcite .
this new scaling has been applied to investigate the viscosity of sqm and the damping time scale that is due to the coupling of the viscosity and r mode @xcite , the quark - diquark equation of state and compact star structure @xcite , the properties of strangelets versus the electric charge and strangeness @xcite , and the new solutions for cfl slets @xcite . in the present paper ,
we use the chirally determined quark mass scaling  @xcite to study the phase transition in neutron stars . for this
we need a completely self - consistent thermodynamic treatment of the eos of quark matter .
the thermodynamic treatment of the system with confinement via the density dependent quark masses has been a long story .
originally , the thermodynamic formulism was regarded the same as the constant - mass case @xcite . in this first treatment ,
the internal pressure can not be zero , and the properties of sqm were rather different from those in the bag model .
but it was later pointed out that the difference was caused by the incorrect thermodynamic treatment  @xcite .
it was found that an additional term is to be added to the pressure and energy expressions @xcite .
this second treatment makes it possible sqm to be self - bound .
however , two serious problems came out : one is the unreasonable vacuum limits , the other one is the discrepancy between the energy minimum and zero pressure .
it was shown that the added term in the pressure , due to the density dependence of quark masses , should not be appended to the energy . after discarding this term in the energy , while keeping it in the pressure ,
the two inconsistencies mentioned above were immediately removed @xcite .
this third treatment has recently been extended to finite temperature @xcite .
the thermodynamic formulism in ref .
@xcite was also adopted in ref .
@xcite , though a different quark mass scaling was used there .
a common feature of the last two thermodynamic treatments  @xcite , as well as other recent references using this model @xcite , is that they all regard the thermodynamic potential the same as in the fermi gas model .
because of the additional term , the pressure becomes obviously not equal to the minus thermodynamic potential density , contradicting the thermodynamic equality @xmath4 for a homogeneous system .
one can also easily check that the fundamental differentiation equality @xmath5 for homogeneous systems at zero temperature was not fulfilled in the mentioned references .    in the rest of this section
, we will present a fully self - consistent thermodynamic treatment of the confinement by density - dependent mass model ( cddm ) .
let us consider a quark model with three flavors . denoting the fermi momentum in the phase space by @xmath6 , the particle number densities
can then be expressed as @xmath7 and the energy density as @xmath8    equations ( [ nimod ] ) and ( [ emod ] ) are familiar expressions , where the summation index goes over all considered particle types . to let the model be valid for both particles and anti - particles , the particle number density , or accordingly ,
the fermi momentum , are formally assumed to be negative for anti - particles .
therefore , in the upper limit of the integration , the absolute value has to be taken .    if the particle masses @xmath9 are constant , the relation between the fermi momenta @xmath6 and the chemical potentials @xmath10 is @xmath11    as is well - known , however , the quark mass depends on density and temperature . in principle , the quark mass scaling should be determined from qcd , which is obviously impossible presently .
based on the in - medium chiral condensates , a cubic scaling law was derived at zero temperature @xcite , and it has been recently extended to finite temperature @xcite . at zero temperature , we have the simple cubic scaling @xmath12 where @xmath13 is the quark current mass , @xmath14 is the total baryon number density , the exponent of density is @xmath15 @xcite , and the constant @xmath16 is to be discussed a bit later .    in the following , we show that the density - dependence of particle masses will modify the fermi momentum , i.e.the relation in eq .
( [ nuimuifree ] ) for free - particle systems should be modified to include interactions .
in fact for the quark flavor @xmath17 we have @xmath18 since the quark masses are density dependent , the derivatives generate an additional term with respect to the free fermi gas model .
we get @xmath19 where @xmath20.\ ] ] the pressure is then given by @xmath21 with @xmath22  being the free - particle contribution : @xmath23.\end{aligned}\ ] ]    due to the additional term in the chemical potential , the pressure also has an extra term .
the inclusion of such a term guarantees that the hughenoltz - van hove theorem is fulfilled in the calculations .    in the quasi - particle model @xcite
, one also has an extra term .
because the quark masses there depend on chemical potentials , and the extra term is not used in the relation between the fermi momenta and chemical potentials , an effective bag constant has to be added to the energy expression to consider confinement @xcite .    in cddm quark model , however
, we do not need a bag constant any more .
quark confinement is achieved automatically by the density dependence of quark masses , or by the strong interaction between quarks .
in fact , the exponent @xmath24 in eq .
( [ mqt0 ] ) is derived from the linear confinement interaction @xcite .        in the present model ,
the parameters are : the electron mass @xmath25 mev , the quark current masses @xmath26 , @xmath27 , @xmath28 , and the confinement parameter @xmath16 . although the light - quark masses are not without controversy and remain under active investigations , they are anyway very small , and so we simply take @xmath29 the current mass of strange quarks is @xmath30 mev according to the latest version of the particle data group @xcite .
conventionally , the stability of strange quark matter ( sqm ) is judged by the minimum energy per baryon @xcite .
if it is less than 930 mev ( the mass of @xmath31fe divided by 56 ) , then sqm is absolutely stable .
if it is bigger than 930 mev but less than 939 mev ( the mass of nucleons ) , then it is metastable .
otherwise , if it is larger than 939 mev , sqm is un - stable . in case of two - flavor
quark matter , it should be no less than 930 mev , in order not to contradict standard nuclear physics .
this is the witten - bodmer hypothesis @xcite .    in fig .
[ dms0 ] , we show the different regimes in the @xmath32-@xmath28 plane . the area below the full line is forbidden where the energy per baryon of two - flavor quark matter is less than 930 mev .
above the dotted line , the energy per baryon of sqm is more than 939 mev , and thus sqm is unstable .
the area bounded by the dotted and dashed lines is the metastable region where the energy per baryon is between 930 mev and 939 mev .
only when the @xmath33 pair is in the range between the full and dashed lines , sqm can be absolutely stable , i.e. , its energy per baryon is less than 930 mev .
therefore , the range of @xmath16 values is very narrow for a chosen @xmath28 value , if the witten - bodmer hypothesis is correct . if we take the modest value @xmath34 mev , for example , then @xmath35 is in the range of 158160 mev .
the lower bound 158 mev is obtained by taking @xmath36 .
if @xmath26 and @xmath27 are given a small finite value , the lower bound can then be a little bit smaller , e.g. 156 mev @xcite .
unfortunately we presently do not have a definite conclusion on the stability of sqm , so we treat @xmath16 as a free parameter .
however , the first condition , i.e. , @xmath16 greater than about @xmath37 , should always be satisfied .
on the other hand , we can connect @xmath16 to the pion mass @xmath38 , pion decay constant @xmath39 , pion - nucleon sigma term @xmath40 , string tension @xmath41 , and the vacuum chiral condensate @xmath42 by @xcite @xmath43 from the known range of the vacuum condensate , we can have an upper bound @xmath44 .
therefore , @xmath35 should not be out of the range ( 156 , 270 ) mev .      as usually done , we consider sqm as a mixture of @xmath45 , @xmath46 , @xmath47 quarks , and electrons .
the relevant chemical potentials @xmath48 , @xmath49 , @xmath50 , and @xmath51 satisfy the weak - equilibrium condition @xmath52    because all particle masses do not depend on the density of electrons , i.e. , @xmath53 , eq .
( [ mui ] ) gives @xmath54 with @xmath55 for @xmath56 quarks , and @xmath57 for electrons .    in eq .
( [ muivsni ] ) , @xmath58 is the second term on the right hand side of eq .
( [ mqt0 ] ) , so we have @xmath59 the pressure is then obtained from eq .
( [ pressure ] ) as @xmath60    substituting these expressions into eqs .
( [ weak1 ] ) and ( [ weak2 ] ) , we have @xmath61 and @xmath62 we also have the baryon number density @xmath63 and the charge density @xmath64 the charge - neutrality condition requires @xmath65 .    for a given total baryon number density @xmath14
, we can obtain the respective @xmath66 , @xmath67 , @xmath68 , and @xmath69 by solving the four equations ( [ qmeq1 ] ) , ( [ qmeq2 ] ) , ( [ qmeq3 ] ) , ( [ qmeq4 ] ) .
the chemical potentials @xmath48 , @xmath49 , @xmath50 , and @xmath51 can then be calculated by eqs .
( [ muivsni ] ) and ( [ muevsne ] ) .
therefore , the energy density of the quark matter is a function of the baryon number density @xmath14 and the charge density @xmath70 , i.e. , @xmath71 or @xmath72 where the two partial derivatives , @xmath73 and @xmath74 are called the baryon chemical potential and charge chemical potential , respectively .
it can be easily shown that they are connected to the quark chemical potentials by @xmath75    in fact , according to the fundamental differentiation equality of thermodynamics , we have @xmath76 on the other hand , we have @xmath77 , @xmath78 , @xmath79 , and @xmath80 from eqs .
( [ weak1 ] ) , ( [ weak2 ] ) , ( [ qmeq3 ] ) , and ( [ qmeq4 ] ) . substituting these four equalities into eq .  ( [ deqthermo ] ) leads to @xmath81 .
comparison of this with eq .
( [ deqbaryonchem ] ) immediately gives eq .
( [ baryonchemmuud ] ) .
mev and @xmath82 mev . ,
width=309 ]    in fig .
[ nvsnb ] , the quark fractions , i.e. , @xmath83 , @xmath84 , @xmath85 , and the @xmath86 times the electron number divided by the total quark number , @xmath87 , have been shown versus the baryon number density for @xmath88 mev and @xmath82 mev .
it is seen that the fraction of up quarks is nearly always one third .
the fraction of down quarks increases rapidly with decreasing densities , while the fraction of strange quarks approaches to zero when the density decreases to a certain lower density .    .
, width=309 ]    in order to compare the relation between the fermi momentum and chemical potential , we plot , in fig .
[ nui - mu ] , the fermi momentum of up ( solid line ) , down ( dotted line ) , and strange ( dashed line ) quarks , respectively , as a function of the corresponding quark chemical potential , in both the present model ( lines with a solid circle ) and the previous model ( lines with an open circle ) .
it is very obvious that the difference is very large , especially at comparatively lower densities . in both models ,
the fermi momentum of up or down quarks is higher than that of strange quarks due to the fact that strange quarks are heavier than up or down quarks .
for the same chemical potential , however , the fermi momentum in the present model is generally bigger than that in the previous model , due to the quark mass - density - dependence which reflects the strong interaction between quarks .     in mev for the solid , dashed , dotted , and dash - dotted cures are ( 200,150 ) , ( 180,120 ) , ( 170,95 ) , and ( 160,80 ) , respectively .
the dash - dot - dot line is for the two - flavor quark matter at @xmath88 mev .
it is very obviously shown that the energy minimum , marked with a full triangle on each line , is located exactly at the same point of the zero pressure indicated by an open circle .
, width=309 ]    figure [ enb ] shows the energy per baryon of quark matter for different parameter sets in the present model .
each line has a minimum , corresponding to the lowest energy state ( marked with a solid triangle ) .
one can see that the pressure at this minimum is exactly zero .
so this special point is marked with an open circle as well . at the same time
, we also display the energy per baryon for the two - flavor quark matter by a dash - dot - dot line .
we see that the two - flavor quark matter is less stable than sqm . even for a smaller @xmath16 value , e.g. , @xmath89
, its energy will finally exceed that of sqm for a bigger @xmath16 with increasing densities .
mev , @xmath90 mev ) and in the bag model ( @xmath91 mev ) .
the former has a minimum at a lower density depending on the value of @xmath16 and @xmath28 , while the latter is always a monotonic function of density .
the velocity of sound in both models has also been given on the right axis .
, width=309 ]    in principle , the cddm quark model contains more physics than the simple bag model . to demonstrate this we plot , in fig .
[ comp ] , the baryon chemical potential in cddm with @xmath92 mev and @xmath90 mev , and in the bag model with @xmath91 mev . in cddm ,
the baryon chemical potential decreases with decreasing density to a certain value depending on @xmath16 and @xmath28 , then it increases very rapidly , i.e. , it saturates at a definite density marked with a bullet .
when the density is lower than the bullet , the derivative @xmath93 becomes negative , and so quark matter is unstable against phase separation and falls apart at lower densities . in the bag model
, however , the baryon chemical potential is always a monotonic function of density , which means that quark matter does not fall apart at any lower densities .
the velocity of sound has also been plotted on the right axis .
we observe that it is , in the bag model , nearly the same as that of a non - interacting fermi gas .
but in cddm , it decreases to zero with decreasing densities . at high densities
, it becomes asymptotically identical to the ultrarelativistic limit , as expected .
since the single particle energies depend on density and temperature via the quark masses , the thermal properties should be founded on the canonical ensemble , but , as is well known , the partition function is not easy to calculate . therefore a different statistical procedure is usually adopted , which is based on the quasiparticle assumption . according to that the energy density is written as @xmath94 where the fermi distribution function is @xmath95/t}}. \label{ass3}\ ] ] if antiparticles are included , the sum must be extended to antiparticles for which @xmath10 must be replaced by @xmath96 . from the landau definition of the single - particle energy extended to finite temperature , we have @xmath97 where @xmath98  is the dispersion relation of free particles .
the extra term @xmath99 can be added to the chemical potential , so defining @xmath100 accordingly , the net density of the particle type @xmath17 is @xmath101,$]or , explicitly , we have @xmath102/t}}\right .    \left.-\frac{1}{1+e^{[\varepsilon_i(p)+\mu^*_i]/t}}\right\ }    \frac{p^2\mbox{d}p}{2\pi^2}.\ ] ] inverting this equation , one determines @xmath103 as a function of @xmath104 so that the free energy density @xmath105 \label{fttot}\ ] ] with @xmath106 \nonumber\\ & & \phantom{g_i\int_0^{\infty}\bigg\ { }   \pm\frac{\mu_i^*}{1+e^{(\sqrt{p^2+m_i^2}\mp\mu_i^*)/t } }   \bigg\ } \frac{p^2\mbox{d}p}{2\pi^2 }   \label{fitmum}\end{aligned}\ ] ] will be a function of respective particle densities , instead of chemical potentials .
one can then determine the real chemical potentials and pressure , according to the well - known relations @xmath107 these quantities will completely describe the thermal equilibrium of pure quark phase and the transition to the quark - hadron mixed phase .
a more detailed analysis , where the thermodynamic formulism is developed to much more extent , is reported in appendix  [ appea ] . in ref .
@xcite , the @xmath10 should be implicitly understood as @xmath108 .
mev and @xmath90 mev .
, width=309 ]    in fig .
[ enbt ] , we give the energy per baryon of sqm as a function of both density and temperature .
the parameters used for this three - dimensional plot are @xmath109 mev and @xmath90 mev .
the brueckner - bethe - goldstone ( bbg ) theory is among the most advanced microscopic theories of nuclear matter . in the recent years
it was recognized that the three - body forces , which are expected to have a dominant role at high nuclear density , also affect the saturation point and , in fact , after including the three - body forces in the brueckner theory , the empirical saturation properties are reproduced quite well  @xcite . in the neutron star interior , where high baryonic density values are reached ,
processes like the excitation of nucleon - antinucleon pairs ( @xmath110-diagrams ) and nucleonic resonances ( together with the production of other hadrons ) sizeably influence the two - body nuclear interaction .
the former process involves the virtual excitation of negative energy states , which is absent from the standard brueckner theory , and thus it represents a pure relativistic effect . from the comparison with the dirac - brueckner theory it turns out that it is by far the most important relativistic effect @xcite .
this process together with nucleon resonances can be incorporated in the interaction as medium virtual excitations in tbfs .
one could guess that many body ( more that three ) forces are also important at high density as large as @xmath111 , but it is hard to imagine that pure baryon matter can exist at so high density . at lower density two and
three - body forces are dominant since the hole line expansion is , roughly speaking , an expansion in density powers .
the global effect of tbfs at high density is strongly repulsive , leading to a remarkable increase of the maximum mass in the study of the neutron star structure .
but also a correct estimate of the saturation point is important , since , as we will see below , in strongly asymmetric nuclear matter the threshold for the transition to mixed nucleon - quark phase can appear close to the saturation density .
therefore the corresponding eos could be used as input for transport - model simulations of heavy ion collisions , where strongly isospin - asymmetric systems are formed in central events .      the brueckner theory extended to tbf
is described elsewhere  @xcite . here
we simply give a brief review of the brueckner - hartree - fock ( bhf ) approximation at finite temperature @xmath112  @xcite .
the starting point is the reaction @xmath113-matrix , which satisfies the bbg equation , @xmath114 where @xmath115 , denotes the single particle momentum , the @xmath116-component of spin and isospin , respectively , and @xmath117 is the starting energy .
the @xmath113-matrix , the pauli operator @xmath118 and the single particle energies @xmath119 depend on the neutron and proton densities and temperature .
the interaction @xmath120 given by @xmath121 where @xmath122 is the bare two body force ( 2bf ) and @xmath123 is an effective 2bf derived by the average of the bare tbf on the third particle as follows @xmath124\displaystyle   \frac{1}{4 } { \rm tr}\sum_n \int { \rm d } { \vec r_3 } { \rm d } { \vec r_3^{\ \prime}}\phi^*_n(\vec r_3^{\ \prime } ) [ 1-\eta(r_{13}',t ) ] [ 1-\eta(r_{23}',t ) ] \\[6 mm ] \times \displaystyle w_3(\vec r_1^{\ \prime}\vec r_2^{\ \prime } \vec r_3^{\ \prime}|\vec r_1 \vec
r_2 \vec r_3 ) \phi_n(r_3 ) [ 1-\eta(r_{13},t)][1-\eta(r_{23},t ) ] .
\end{array } \label{eq : tbf}\ ] ]    since the defect function @xmath125 is directly determined by the solution of the bbg equation @xcite , @xmath123 must be calculated self - consistently with the @xmath113 matrix and the s.p
.  potential @xmath126 on the basis of bbg equation .
it is clear from eq .
( [ eq : tbf ] ) that the effective force rising from the tbf in nuclear medium is density and temperature dependent through the defect function .
a detailed description and justification of the method can be found in ref .
@xcite , including a discussion on the averaging procedure .
the validity of such a procedure has been numerically tested in the comparison between the bhf eos plus z - diagrams with @xmath127 meson exchange and the dirac - bhf eos , which are expected to be equal .
the calculation @xcite gives an impressive agreement between the two eos s , although the tbf due to the @xmath110-diagrams is averaged according to eq .
( [ eq : tbf ] ) . for @xmath122
we adopt the argonne @xmath128 two - body interaction  @xcite .
the tbf is constructed from the meson - exchange current approach  @xcite and contains virtual particle ( @xmath2 and @xmath129 ) excitations and , in addition , relativistic effects induced by the excitations of particle - antiparticle pairs .
this description of the interaction is not completely consistent since , in principle , the two and three body forces should be derived from the same meson parameters , but a recent calculation @xcite replacing argonne potential with bonn potential @xcite and the tbf built up with the bonn meson parameters substantially leads to the same results .
let us start with symmetric nuclear matter . in the bhf approximation
, the thermodynamic potential can be written @xmath130 , \\
w_l&=&-\frac{1}{2}\sum_k f_k(t ) u_k(t),\end{aligned}\ ] ] where @xmath131 is the thermodynamic potential for a system of independent particle with single particle spectrum @xmath132 , and the @xmath133 the sum of all linked cluster diagrams to the lowest order in the hole line expansion@xcite .
@xmath134 is the selfconsistent mean field , @xmath135 the finite temperature bhf approximation suffers from the same difficulty as any strongly interacting fermi system .
the difficulty is the same as in the cddm quark model , since in both cases the s.p .
spectrum is density and temperature dependent . whereas at @xmath136 the density ( or fermi momentum )
can be fixed , at @xmath137 the role of density is taken by chemical potential @xmath138 ( grand canonical ensemble ) , and one is forced to fix the chemical potential at each iteration due to the presence of the fermi distribution @xmath139/t)\}^{-1}$ ] . since this procedure does not converge @xcite , one should fix the density and invert the equation relating density and chemical potential , @xmath140 which is not a viable task .
the usual approximation is to drop out the derivative in the previous equation , which corresponds to the quasiparticle approximation above discussed within the cddm quark model .
so doing , the resulting @xmath141 looses its meaning of chemical potential . in this approximations
the energy and entropy densities are given by : @xmath142 , \\
s_n & = & -\frac{1}{v}\sum_k\big\{f_k(t)\ln f_k(t ) \nonumber\\ & & \phantom{-\sum_k\big\ { }    + \big[1-f_k(t)\big ] \ln\big[1-f_k(t)\big]\big\}.\end{aligned}\ ] ] after one calculates @xmath141 in terms of @xmath143 , the thermodynamics is developed from the free energy density @xmath144 . the free energy per particle , calculated from bhf approximation ,
is depicted in fig .
[ eoslast ] . due to the difficulty of extending the bhf code to very high temperature an extrapolation from the real numerical results to high @xmath112 has been performed adopting the so called frozen approximation based on @xmath112-independent single particle spectrum , i.e.  the latter
is frozen at @xmath136 .
this turns out to be a good approximation up to @xmath145 mev @xcite .
the relevant thermodynamical quantities , i.e.  chemical potentials and pressure , are derived from free energy as follows @xmath146 let us consider asymmetric nuclear matter with baryon density @xmath147 and asymmetry parameter @xmath148 , where @xmath149 ( @xmath150 ) is the neutron ( proton ) density .
the baryon chemical potentials can be expressed as @xmath151 where @xmath152 is the free energy density .    assuming the parabolic law for the latter , we get the simple expression for the chemical momentum isotopic shift @xmath153 where @xmath154 is the symmetry free energy density . the parabolic law is well satisfied at low density , but at high density additional terms of the @xmath155-expansion
must be considered .    in neutron star inner core nuclear matter
is supposed to be in @xmath155-equilibrium under the condition of charge neutrality . assuming that only electrons are present ( the muon contribution is negligible )
, the two preceding conditions require @xmath156 @xmath157 is the net charge density of nuclear matter .
it is zero for pure neutral nuclear matter . for a given set of @xmath158
, we can solve the chemical potentials @xmath159 , @xmath160 , and @xmath51 from the above equations . then all other quantities can be obtained for a fixed temperature . in other words ,
all thermodynamic quantities can be regarded as a function of the nucleon density @xmath143 , charge density @xmath157 , and temperature @xmath112 . at zero temperature , for example
, the energy density can be regarded as a function of @xmath143 and @xmath157 , i.e. , @xmath161 . with a similar approach as in the preceding section to obtain eq .
( [ baryonchemmuud ] ) , we can easily show that the baryon chemical and charge chemical potentials of nuclear matter can be expressed as @xmath162    the system turns out to be in a strongly isospin asymmetric state .
the isotherms of free energy and pressure of nuclear matter in @xmath155-equilibrium are shown in fig .
[ beos ] .
three values of the temperature are considered.,width=309 ]
let us study the nuclear matter , consisting of nucleons and electrons , in equilibrium with a gas of @xmath45 , @xmath46 , @xmath47 quarks and electrons . according to glendenning  @xcite
, we assume the total charge conservation , in addition to total baryon and energy conservation .
now we first consider the case of zero temperature , and then extend to finite temperature .
the conservation laws can be imposed by introducing the quark fraction @xmath163 defined as @xmath164 where @xmath165 is the total volume , @xmath166 is the volume occupied by quarks
. then the total baryon density is @xmath167 the total electric charge is @xmath168 and the total energy density is @xmath169 where @xmath143 , @xmath157 , and @xmath170 are , respectively , the baryon number density , electric charge density , and energy density of nuclear matter , while @xmath14 , @xmath70 , and @xmath171 are the corresponding quantities of quark matter .
@xmath170 is a function of @xmath143 and @xmath157 , @xmath171 is a function of @xmath14 and @xmath70 , i.e. , @xmath161 , @xmath71 . differentiating eq .
( [ etot ] ) , one obtains @xmath172 on the other hand , differentiating eqs .
( [ rhotot ] ) and ( [ qtot ] ) at a given pair of @xmath173 and @xmath174 , we have @xmath175    to minimize @xmath176 , we substitute eqs .
( [ drhoh ] ) and ( [ dqh ] ) into eq .
( [ de01 ] ) .
then setting @xmath177 , we find @xmath178 where @xmath179 the conditions in eq .
( [ gibbscond ] ) are nothing but the gibbs ones , i.e. , the baryon chemical potential , the charge chemical potential , and the pressure in nuclear and quark matter should be equal to each other to minimize the total energy of the mixed phase .    in the previous two sections , we have linked the baryon chemical potential and charge chemical potential to the respective constituent particle chemical potentials in eqs .
( [ baryonchemmuud ] ) and ( [ baryonchemnm ] ) . as application of these equalities
, we immediately see that the first two equations in ( [ gibbscond ] ) are equivalent to @xmath180    in general , all other chemical potentials in quark sector can be related to @xmath48 and @xmath49 , e.g. , @xmath78 , @xmath77 .
similarly , all chemical potentials in nuclear sector can be linked to @xmath159 and @xmath160 , e.g. , @xmath181 .
therefore , eq .  ( [ munpud ] )
means that we can choose either ( @xmath48,@xmath49 ) or ( @xmath159,@xmath160 ) as the two independent chemical potentials .
the latter can then be determined by solving the charge neutrality equation and the pressure balance equation for a given total baryon number or a given quark fraction .    at finite temperature , we similarly have the phase equilibrium condition @xmath182    the condition in eq .
( [ phasetcon03 ] ) only tells us that the temperature in nuclear and quark sectors are equal , so we have a common temperature @xmath112 . the chemical equilibrium condition in eq .
( [ phasetcon02 ] ) is equivalent to that in eq .
( [ munpud ] ) .
therefore , we still have only two independent chemical potentials . for a given total density @xmath173 at a fixed temperature @xmath112 , the two independent chemical potentials and the quark fraction @xmath163 can then be determined by solving the three equations in eqs .
( [ rhotot ] ) , ( [ qtot ] ) with @xmath183 , and ( [ phasetcon01 ] ) .
similar to the case at zero temperature , the lower critical density @xmath184 , which separates the nuclear and mixed phases , is defined by @xmath185 , while the critical density @xmath186 between the mixed and quark phases is determined by @xmath187 .
, width=309 ]    in fig .
[ mixe17095 ] , we display , for the parameter @xmath32 = 170 mev , the density dependence of the energy per baryon in pure nuclear matter , in pure quark matter , and in the mixed phase .
the quark fraction has also been depicted on the right y - axis . in fig .
[ mixp17095 ] the corresponding pressure is reported .
it is seen that the nuclear matter is the most favorite phase at lower densities , and the quark matter is the most stable phase at higher densities , while at intermediate densities , the mixed phase has the lowest energy .
( cures with an open circle ) and 20 mev ( cures with an solid circle ) . the nuclear density @xmath143 and the quark density @xmath14 in the mixed phase
are also shown on the right axis .
, width=309 ]    the quark baryon number density @xmath14 and nuclear density @xmath143 are also plotted on the right axis of fig .
[ mixp17095 ] .
we see that the quark density is always higher than the nuclear density .
the transition from hadron phase to mixed phase occurs at the density a bit less than @xmath188 @xmath189 , well below the saturation density .
but it is hard to observe in terrestrial laboratories , since the nuclear matter so far realized in exotic nuclei or heavy - ion collisions is much less neutron rich .
the transition from mixed phase to pure quark phase occurs at the total density @xmath190 @xmath189 , where the nuclear density is only @xmath191 @xmath189 .
the density range of mixed phase is only slightly depending on the temperature , at least in temperature range of interest in neutron stars .
the critical densities depend on the parameter @xmath16 . in the left panel of fig .
[ rhocd ] , both nuclear critical density ( the solid line , separating the pure nuclear phase and the mixed phase ) and the quark critical density ( the dashed line , delimiting the pure quark phase ) are displayed as a function of @xmath35 .
if @xmath192 mev , the two critical densities approach zero , and accordingly sqm is absolutely stable .
when 161.6 mev @xmath193 162.5 mev , mixed phase can exist at any lower densities .
only when @xmath194 mev , nuclear matter is more stable at lower densities
. if we have only two flavor quarks in the quark sector , the critical densities are usually higher . in the same panel
we also plot the lower critical density for the two - flavor case . because we know in our real world the two - flavor quark matter does not exist below the saturation density , @xmath35 should be on the right of the first full dot at @xmath195 mev ( the intersection of the dotted and dot - dashed lines ) in fig .
[ rhocd ] . on the other hand , to let sqm have a chance to appear below the saturation density , @xmath35 should be on the left of the second full dot where @xmath196 mev . in plotting fig .
[ mixe17095 ] and [ mixp17095 ] , we adopt @xmath197 mev . the temperature dependence of the lower and higher critical densities is also plotted in fig .  [ rhocd ] ( right panel ) .
the two lines at fixed @xmath16 mark the boundaries of the three phases .    .
the horizontal line is the nuclear saturation density .
the dotted line is the quark critical density in the case of two - flavor quark system .
right panel : phase diagram for three values of the confinement parameter and @xmath198 mev .
, width=309 ]
with the equation of state that has the mixed and/or quark phase derived in sec
.  [ phase ] , we are ready to study the structure of hybrid stars by solving the tolman - oppenheimer - volkov equation @xmath199 where @xmath200 mev@xmath201 is the gravitational constant , @xmath202 is the distance from the center of the star , @xmath203 and @xmath204 are the energy density and pressure at the radius @xmath202 , respectively .
the subsidiary condition is @xmath205 with @xmath206 being the mass within the radius @xmath202 .    at variance with pure nuclear or quark stars ,
a hybrid star contains pure quark matter in the core , pure nuclear matter near the outer part , and , in between , a mixed phase of the quark and nuclear matter . in this case , therefore , we must use the eos in whole density range .
mev for four values of the confinement parameter @xmath35 of 170 , 180 , 190 , 200 mev . ,
width=309 ]    the resulting gravitational mass for the hybrid star is plotted in fig .
[ mass ] , as a function of both radius and central density , for four values of the confinement parameter @xmath35 of 170 , 180 , 190 , 200 mev .
the main effect of the phase transition is , as expected , a large reduction of the maximum mass due to softening of the eos .
concerning the confinement parameter @xmath16 , we observe a slight decrease of the maximum mass when the @xmath16 value becomes smaller , from 1.74@xmath207 at @xmath35 = 200 mev down to 1.6 , 1.55 , 1.52 @xmath207 corresponding to @xmath35 = 190 , 180 , 170 mev .
this can be easily understood in this way : since quark phase is occurring at lower density for smaller values of @xmath16 , sometimes even less than the nuclear saturation density ( for example only 0.15 @xmath189 for @xmath35 = 170 mev ) , then the quark population of the star would be more numerous , since the stronger softening of the eos may only support less gravitational mass .
in particular , the most favorable case of @xmath35 = 180 and 190 mev , when the quark confinement appears at around 2@xmath208 consistent with the heavy ion experiment , the predicted maximum mass can match very well the s branch of neutron stars mentioned by haensel et al.with an error less than @xmath209 @xcite .        in order to study the temperature effect on the mass of hybrid stars
, we plot the star mass , vs the star radius or central density , at @xmath136 and @xmath210 mev in fig .
it is seen that the temperature influence on the maximum mass is very limited .
otherwise the effect is a quite strong increase of the ns radius for a fixed amount of gravitational mass .
but , larger is the mass smaller is the radius variation , similarly to a few other calculations @xcite .
the condensation of @xmath211 mesons in neutron stars is widely discussed in the literature ( see refs .
@xcite and references therein quoted ) . in dense matter
, the condensation of @xmath211 mesons is originated by the reaction @xmath212 if the effective mass of the @xmath211 drops below the chemical potential of the electron , this reaction becomes possible in dense matter , indicating the presence of kaon condensation . since almost all the studies of kaon medium properties @xcite
suggest the consistent picture that the attraction from nuclear matter would bring the @xmath211 mass down , so the threshold condition for the onset of @xmath211 condensation @xmath213 , which follows from eq .
( [ eq : kaon.1 ] ) , could be fulfilled in the center of neutron stars @xcite at @xmath214 .
however , the deconfinement phase transition from the hadronic phase to the quark phase occurs also at rather low density , which leads to the onset of mixed phase in neutron stars .
then it is very interesting to explore how the quark deconfinement affects the @xmath211 condensation threshold .
mev for three selected values of the confinement parameter @xmath35 of 180 , 200 , 225 mev .
, width=309 ]    we take the antikaon dispersion relation constrained by the heavy - ion data as empirical indication of an attractive antikaon potential in dense matter @xcite , and combine this with the bhf model of nuclear matter together with the above quark model without bag constant , to calculate the effective kaon mass in hybrid stars .
the result is illustrated in fig .
[ kaon ] , where the effective kaon mass and the electron chemical potential are shown in normal nuclear matter and especially in mixed phase for three selected values of the confinement parameter @xmath16 .
one sees that in normal nuclear matter , @xmath211 medium mass decreases with increasing density and meets the electron chemical potential at @xmath215 0.6 @xmath189 ( the solid bullet in fig .
14 ) , so the kaon condensation would be present for sure in this case ; however , once the quark phase sets in at the total density of 0.26 , 0.4 , 0.53 @xmath189 respectively for @xmath35 = 180 , 200 , 225 mev , the decreasing speed of the kaon effective mass in the matter slows down a little bit .
in addition , more conclusively , the electron population begin to decrease instead , for the reason that the electric charge neutrality can be achieved more efficiently through the charged quarks themselves . as a result
the threshold condition of @xmath211 condensation is much more difficult to satisfy , unless the confinement parameter @xmath16 is chosen to be extremely high ( at least @xmath216 mev ) when the presence of quark phase is pushed to very high densities , then one may expect that @xmath217 can finally equal to @xmath51 ( which is indicated with a circle )
. those high values of @xmath16 may even not be realistic , therefore we would like to conclude that the inclusion of quark phase may make the kaon condensation impossible in neutron stars , or at least hinder it very strongly .
in this paper the neutron star inner structure has been studied .
the lack of strong observational constraints demands for sophisticated models of the ns composition and interaction mechanisms . in this study
we only included baryons and quarks in equilibrium with leptons and kaons .
since in the ns interior high baryon densities are reached , n@xmath218 and nucleonic excitations are expected to play a major role .
their effects can be incorporated in a three body force .
the baryon eos in weak coupling equilibrium with electrons was derived within the bbg theory suitably extended so to include the three body force . in the quark sector
, we adopted the semiphenomenological cddm quark model , which exhibits a confinement mechanism alternative to the crude mit bag model .
furthermore , in contrast with the extension of the mit model , where the density dependence is introduced artificially @xcite , cddm quark model shapes its density dependence in agreement with chiral requirements @xcite .
hadrons and quarks in @xmath155-equilibrium were considered at zero and finite temperature . since typical temperatures of protoneutron stars are as high as @xmath219 mev and beyond , rigorous thermodynamic potentials have to be derived .
this was possible with cddm quark model , due to its simplicity , whereas some approximations are needed to calculate the hadron phase in the brueckner theory .
the latter is still a main drawback of the finite temperature microscopic theory of strongly interacting fermi systems .
the transition from the low density hadron phase to high density quark phase in beta equilibrium was studied under the glendenning hypothesis of total charge neutrality .
the gibbs construction enabled to follow the evolution of the mixed hadron - to - quark phase , varying the temperature @xmath112 and the confinement parameter @xmath16 .
the eos , in terms of pressure , energy density and chemical composition , was constructed as a function of @xmath112 and @xmath16 .
moreover the phase diagram @xmath112-@xmath16 was also depicted .
the transition density from hadron to mixed phase is strongly dependent on the confinement parameter and , for d small enough , it becomes lower than the nuclear saturation density .
however , at that point nuclear matter is so strongly neutron rich that it hard to imagine that it can be observed in terrestrial laboratory experiments . at high temperature
the transition density is comparatively smaller .
the tov equations have been solved in the above discussed model of neutron stars .
the transition to quark matter produces a strong reduction of the maximum mass , from @xmath220 to @xmath221 for the lowest value of the confinement parameter @xmath16 , but the corresponding radius is not changed , being in both cases about @xmath222 km .
the other @xmath223-@xmath224 configurations ly in a range @xmath225 km depending on the value of @xmath16 . with increasing @xmath16 , the maximum mass increases from @xmath226 to @xmath227 solar masses .
other quark models have been adopted to describe the ns transition to a deconfined phase ( see , for instance , refs .
@xcite ) . except for the nambu - jona lasinio model which exhibits instability of neutron stars @xcite
, the other models point to a softening of the eos of nuclear matter .
the mit model with density - dependent bag constant predicts a beta - stable quark phase quite similar to our model @xcite , but the comparison on the ns structure turns out to be difficult , since in that calculation the transition to quark phase is built on top of hyperonized nuclear matter .
the effect of temperature on the maximum mass is negligible up to @xmath210 mev , but as expected the heated system can support much less gravitational mass in the same volume for all @xmath223-@xmath224 configurations .
thus , an evolutionary study of isolated neutron stars would show a strong compression from the newborn phase to the long era phase as fast as the star cools down from @xmath219 mev to approximately zero .    beyond the confinement parameter ,
the transition to quark phase is interrelated to other properties , in particular the possible onset of kaon condensation . in our analysis
we concluded that the quark phase could not be compatible with the kaon condensation in neutron stars , or at least could hinder it very strongly .
the model developed in this paper misses some important aspects .
the first one is the effect of neutrinos .
neutrinos play the major role in new born neutron stars , when they are still trapped in the interior . in this case
the neutronization process is still hindered and the eos is that of symmetric nuclear matter . moreover ,
the onset of kaon condensation is shifted to higher densities @xcite and probably kaons have no more any chance to compete with quarks .
the second one is the inclusion of hyperons and their competition with other mechanisms such as kaon condensation . in refs .
@xcite the hybrid stars are studied by including hyperons in the hadron phase , which makes the hadron eos very soft and brings the maximum mass to @xmath228 .
the introduction of the quark phase rises up again the maximum mass to @xmath229 , a value which is consistent with our prediction without hyperons . on the other hand ,
the appearance of hyperons much depends on the threshold of kaon production .
therefore the interplay between hyperons and kaons turns out to be quite important and deserves additional investigation , as soon as more reliable empirical inputs will be available , especially on the hyperon - nucleon and hyperon - hyperon interaction .
one of us ( u.l . ) is indebted to dr m. baldo for several discussions on statistical thermodynamics .
g.x .  would like to thank support from nsfc ( 10675137 , 10375074 , 90203004 ) and kjcx2-yw - n2 , and also acknowledge warm hospitality at the dipartimento di fisica e astronomia ( universit di catania ) , the laboratori nazionali del sud ( infn , catania ) , and the mit center for theoretical physics ( mit - ctp ) ..
let s start from the fundamental thermodynamic differentiation relation @xmath230 where @xmath231 is the entropy density , and ( anti)particles are assumed to be uniformly distributed in a volume @xmath165 . by using the free energy density @xmath232
, it becomes @xmath233or , equivalently , @xmath234 because of the uniformity , the free energy density has nothing to do with the volume .
we thus have @xmath235    at finite temperature , we should consider both particles and anti - particles , and particles / anti - particles are not always located bellow the fermi energy .
therefore , the net particle number densities can be expressed as @xmath236 where the superscripts @xmath237  indicates particles , and @xmath238 signifies antiparticles : @xmath239    in eqs .
( [ nitmuspm ] ) and ( [ nitmum ] ) , @xmath108 are effective chemical potentials of respective particles .
if quark densities are not density and/or temperature dependent , @xmath108 are nothing but the actual chemical potentials . in our present case , however , quark masses depend on both density and temperature , in order to include the strong interaction between quarks . therefore , the real chemical potentials should be derived according to fundamental thermodynamic laws .      to determine the thermodynamic properties , we need to give a characteristic function . at zero temperature
, we use the energy density in eq .
( [ emod ] ) due to zero entropy .
now the temperature @xmath112 and the densities @xmath104 are chosen as the independent system variables , we should , therefore , choose the free energy in eq .
( [ fttot ] ) as the characteristic function .          defining @xmath247\ ] ] with @xmath248   p^2\mbox{d}p,\ ] ] then we can write @xmath249 . substituting this into eq .
( [ stmus ] ) and ( [ muimus ] ) , then applying eq .
( [ dnidmui ] ) and @xmath250 , we have @xmath251 and @xmath252 obviously , the free energy density can be given as @xmath253 the energy density is obtained by @xmath254 , giving @xmath255 and the pressure is obtained by substituting eq .
( [ muitmusexp ] ) into eq .
( [ ptexpr ] ) : @xmath256 eqs .
( [ stmusexp])-([ptmusexp ] ) here are in complete accordance with the eqs .
( 58)-(61 ) in @xcite if one regards the @xmath10 there as @xmath108 and @xmath257 as @xmath22 which were not explicitly stated .
the expression in eq .
( [ muitmusexp ] ) is special in the present paper because we need the real chemical potential to investigate the mixed phase .      in the above derivation , we choose volume @xmath165 , the temperature@xmath112 , and the particle number densities @xmath104 as the independent system variables . in this case , the free energy is the characteristic function from which we get the complete set of thermodynamic functions . for this purpose
we have defined the intermediate variables @xmath108 in eqs .
( [ nitmuspm ] ) and ( [ fitmum ] ) . because the quark matter we are considering is a strongly interacting system , the relations between the chemical potentials and the densities are , in principle , not the same as those of a free fermi gas .
however , with the ` effective ' chemical potentials @xmath108 , the densities and the free energy are really of the same form as those of a non - interacting fermi gas .
this is what the ` effective ' means .
the actual chemical potentials of each type of particles are determined from the fundamental thermodynamic equality ( [ dfexpr ] ) which results in eq .
( [ muitmusexp ] ) .
peng , nucl .
a _ 747 _ , 75 ( 2005 ) ; g.x .
peng , m. loewe , u. lombardo , and x.j .
wen , nucl .
b ( proc . suppl . ) * 133 * , 259 ( 2004 ) ; g.x .
peng , h.c .
chiang , p.z .
ning , u. lombardo , and m. loewe , int .
a * 18 * , 3151 ( 2003 ) ; g.x .
peng , u. lombardo , m. loewe , and h.c .
chiang , phys .
b * 548 * , 189 ( 2002 ) .
y. zhang , r.k .
su , phys .
c * 65 * , 035202 ( 2002 ) ; y. zhang , r. k. su , s. q. ying , and p. wang , europhys . lett . * 56 * , 361 ( 2001 ) ; y. zhang and r. k. su , phys . rev .
c * 67 * , 015202 ( 2003 ) .
p. wang , phys .
c * 62 * , 015204 ( 2000 ) .
wiringa , v.g.j .
stoks , and r. schiavilla , phys .
c * 51 * , 38 ( 1995 ) .
r. machleidt , adv .
nucl . phys .
* 19 * , 189 ( 1989 ) ; r. brockmann and r. machleidt , phys .
c * 42 * , 1965 ( 1990 ) ; r. machleidt , phys . rev .
c * 63 * , 024001 ( 2001 ) .
brown , c .- h .
lee , and m. rho , phys . rept . * 462 * , 1 ( 2008 ) . g.q .
lee , and g.e .
brown , phys .
lett . * 79 * , 5214 ( 1997 ) ; f. weber , prog .
. phys . * 54 * , 193 - 288 ( 2005 ) ; and quoted therein .