View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Xiamen University Institutional Repository CHIN.PHYS.LETT. Vol. 25, No. 12 (2008) 4233 Hot Nuclear Matter Equation of State and Finite Temperature Kaon Condensation ∗ LI Ang(o[)1∗∗, ZUO Wei()2, Schulze H. J.3, Lombardo U.3 1Department of Physics and Institute of Theoretical Physics and Astrophysics, Xiamen University, Xiamen 361005 2 Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000 3 INFN-LNS, Via Santa Sofia 44, I-95123 Catania, and Department of Physics and Astrophysics, Catania University, Via Santa Sofia 64, I-95123, Italy (Received 7 June 2008) We perform a systematic calculation of the equation of state of asymmetric nuclear matter at finite temperature within the framework of the Brueckner{Hartree{Fock approach with a microscopic three-body force. When applying it to the study of hot kaon condensed matter, we find that the thermal effect is more profound in comparison with normal matter, in particular around the threshold density. Also, the increase of temperature makes the equation of state slightly stiffer through suppression of kaon condensation. PACS: 21. 65. +f, 26. 60. +c, 24. 10. Cn, 97. 60. Jd, 13. 75. Cs One of the challenging problems in nuclear physics the strangeness content of the proton, which is still is to elucidate the behaviour of nuclear matter under quite controversial.[3−5] Neutrino-trapping may shift high-density and/or high-temperature environments. the kaon threshold to higher densities as compared to Despite the importance of the structure and properties the neutrino-free matter.[6] Also, the presence of other of β-stable matter at extreme densities of several times strange particles (for instance, hyperons or more gen- −3 of normal nuclear matter density (ρ0 ∼ 0:17 fm ), as erally strange-quark degrees of freedom) was found yet, its internal constitution and the equation of state to push the onset of kaon condensation to higher (EoS) are not yet known with certainty. Strangeness densities, even out of the physically relevant density may occur in the form of hadrons (such as hyperons or regime.[17] a K− meson condensation) or in the form of strange Overall most of the previous investigations have quarks. The existence of these strange phases may been carried out for cold matter, and have underesti- have important consequences for the structure of com- mated the dependence of kaon condensation on tem- pact stars, for the cooling dynamic of a remnant star perature. That is not perfectly justified because ther- after supernova explosion and the subsequent forma- mal effects are found to significantly affect the prop- tion of a neutron star. erties of proto-neutron stars (PNS).[18] Therefore in With respect to kaons, Kaplan and Nelson this Letter we extend our previous work[15;16] to hot [1] showed in tree order with SU(3)L × SU(3)R chi- matter. Our objective here is to investigate the im- ral Lagrangian that at a density around 3ρ0, the pact of a kaon condensation on NS matter at finite ground state of baryonic matter might contain a Bose{ temperature. Einstein condensation of negatively charged kaons. In the kaon condensed phase of neutron star mat- Following that, kaon condensation has prompted ex- ter, the energy density consists of three parts of con- tensive investigations and discussions[2−16] on its im- tributions, i.e. plications for astrophysical phenomena in a neutron " = " + " + " ; (1) star (NS). An improved mechanism for kaon conden- NN lep KN sations going beyond the tree order which is consis- where "NN is the baryonic part, "lep denotes the con- tent with kaon{nuclear interactions was proposed by tribution of leptons, and KN is the contribution from Brown, Kubodera, Rho and Thorsson.[2] This work the kaon{nucleon interaction. as well as others that followed it confirmed that the In the present work, we employ the Brueckner{ K K critical density ρc lies in the range 2ρ0 < ρc < Hartree{Fock (BHF) approach for asymmetric nuclear [18−21] 4ρ0, more or less in accord with the first prediction matter at finite temperature to calculate the of Ref. [1]. After a kaon condensation is present, baryonic contribution to the EoS of stellar matter. proton abundance increases dramatically, and even The essential ingredient of the BHF approach is the in- antileptons are allowed to exist. Nevertheless, the teraction matrix G, which satisfies the self-consistent critical density depends sensitively on the value of equations ∗Supported by the National Natural Science Foundation of China under Grant No 10575119, the Knowledge Innovation Project of Chinese Academy of Sciences under Grant No KJCX3-SYW-N2, and the Asia-Europe Link project (CN/ASIA-LINK/008(94791)) of the European Commission. ∗∗Email: [email protected] c 2008 Chinese Physical Society and IOP Publishing Ltd 4234 LI Ang et al. Vol. 25 the correlations at T 6= 0 are assumed to be essen- tially the same as at T = 0. This means that the s.p. X 1 G(ρ, xp; !) = vNN + vNN potential Ui(k) at finite temperature for the compo- ! − (k1) − (k2) k1k2 nent i is approximated by the one calculated at T = 0. · jk1k2i[1 − n(k1)][1 − n(k2)]hk1k2j Then the free energy density has the following simpli- fied expression: · G(ρ, xp; !); (2) X X k2 1 and f = n (k) + U (k) − T s ; (6) i 2m 2 i i i i X 0 0 0 0 k U(k; ρ)=Re n(k )hkk jG[ρ; e(k)+e(k )]jkk ia; (3) 0 k ≤kF where the entropy density si for component i is ob- tained from the free gas model with a spectrum where ! is the starting energy; xp = ρp/ρ is the of i(k) (see, for example, Ref. [18]). It turns out proton fraction, i.e. ρp and ρ the proton and the that the assumed independence is valid to a good total baryon density, respectively; (k) is the single accuracy,[18;19] at least for not too high temperature. particle energy given by (k) = ¯h2k2=(2m) + U(k). The realistic nucleon-nucleon (NN) interaction vNN adopted in the present calculation is the Argonne v18 two-body force[22] supplemented with the microscopic three-body force (TBF) given in the following. At finite temperature, n(k) is a Fermi distribution. For a given density and temperature, Eqs. (2) and (3) have to be solved self-consistently with the following equation for the auxiliary chemical potentialµ ~, X X h (k) − µ~ i−1 ρ = n(k; T ) = exp + 1 : (4) T k k The microscopic TBF adopted in the present calcula- tions is constructed from the meson-exchange current approach, which has been reduced to an effective two- Fig. 1. Finite temperature EoS for symmetric (left-hand body interaction in r-space as panel) and purely neutron (right-hand panel) matter. The upper panels show the internal energy and the lower pan- 0 0 eff hr1r2jV3 (T )jr1r2i els the free energy per particle as a function of the nucleon density. The temperatures vary from 0 to 50 MeV in steps 1 X Z = Tr n(k; T ) dr dr0 φ∗ (r0 ) of 10 MeV. 4 3 3 n 3 k In Fig.1 we display the EoSs obtained following 0 0 × (1 − η(r13;T ))(1 − η(r23;T )) the above-discussed procedure, for symmetric nuclear 0 0 0 × W3(r1r2r3jr1r2r3)φn(r3) matter and purely neutron matter. In the upper pan- els we display the internal energy per particle, whereas × (1 − η(r13;T ))(1 − η(r23;T )); (5) the lower panels show the free energy as a function of where the trace is taken with respect to the spin and the baryon density, for several values of temperature isospin of the third nucleon. The function η(r; T ) between 0 and 50 MeV. We notice that the free en- is the defect function. Since the defect function is ergy of symmetric matter is a monotonically decreas- directly determined by the solution of the Bethe{ ing function of temperature. At T = 0 the free energy Goldstone (BG) equation,[23] i.e. Eq. (2), and depends coincides with the internal energy and the correspond- on temperature. The transformation of the TBF to ing curve is just the usual nuclear matter saturation the above effective interaction entails a self-consistent curve. In contrast, the internal energy is an increasing coupling between the TBF and the Brueckner proce- function of temperature. The effect is less pronounced dure of solving the BG equation. It is obvious from for pure neutron matter due to the larger Fermi en- eff Eq. (5) that the temperature dependence of V3 (T ) ergy of the neutrons at given density. We notice the comes from the Fermi distribution n(k; T ) and the de- results of this microscopic TBF are always stiffer than fect function η(r; T ), therefore it depends implicitly those of the phenomenological Urbana TBF.[18] on temperature via density. As to the thermal kaon-nucleon part, we apply the To save computational time and simplify the nu- new formulation constructed by Tatsumi et al. and merical procedure, in the following we introduce the treat fluctuations around the condensate within the so-called frozen correlations approximation,[18] i.e., framework of chiral symmetry. A detailed description No. 12 LI Ang et al. 4235 of the procedure is given in Refs. [27,28]. In addition and charge neutrality conditions leads to the following to the classic kaon contribution, the thermodynamic three coupled equations:[4;5;27;28] potential of the thermal kaon contribution is intro- ρ duced as follows: 0 = f 2 sin θ m2 + (a x + a + 2a )m K 1 p 2 3 s f 2 Z d3p th th −βE+(p) −βE−(p) ρ(1 + xp) 2 @(ΩK =V ) ΩK =TV ln 1−e 1−e ; (7) − 2µ − µ cos θ + ; (2π)3 K 4f 2 K @θ (11) and the kaon{nucleon energy density obtained in this way is (θ is the amplitude of the condensation) and θ θ µ = µ − µ = 4(1 − 2x )F sec2 " = f 2 − µ2 sin2 θ=2 + 2m2 sin2 K n p p sym 2 KN K K 2 2 θ Z d3p − 2a1ms tan ; (12) + [E (p)f (E (p)) 2 (2π)3 − B − + E+(p)fB(E+(p))] + µK QK ; (8) QK + Qe + Qµ = Qp(= xpρ); (13) [4;5] where E± are the zero-momentum kaonic energies where µ represents the chemical potential of the neg- and the kaonic charge density QK is given by atively charged particles.