CSIRO PUBLISHING Publications of the Astronomical Society of Australia, 2005, 22, 292–297 www.publish.csiro.au/journals/pasa
Strange Star Equations of State Revisited
Debora P. MenezesA,B,C and Don B. MelroseB
A Depto de Física - CFM - CP476, Universidade Federal de Santa Catarina, Florianópolis - SC, Brazil B School of Physics, University of Sydney, Sydney NSW 2006, Australia C Corresponding author. Email: [email protected], [email protected]
Received 2005 June 8, accepted 2005 September 8
Abstract: Motivated by recent suggestions that strange stars can be responsible for glitches and other obser- vational features of pulsars, we review some possible equations of state and their implications for models of neutron, hybrid, and strange stars. We consider the MIT bag model and also strange matter in the colour– flavour locked phase. The central energy densities for strange stars are higher than the central densities of ordinary neutron stars. Strange stars are bound by the strong force and so can also rotate much faster than neutron stars. These results are only weakly dependent on the model used for the quark matter. If just one of the existing mass-to-radius ratio constraint is valid, most neutron stars equations of state are ruled out, but all the strange stars equations of state presented in this work remain consistent with the constraint.
Keywords: stars: interiors — equation of state
1 Introduction when the neutrinos are still trapped (Panda et al. 2004a; Qualitatively, a neutron star is analogous to a white dwarf Menezes & Providência 2004a) and for the subsequent star, with the pressure due to degenerate neutrons rather (deleptonized) phase after the neutrinos escape (Menezes than degenerate electrons. General relativity is signifi- & Providência 2003, 2004b, 2004c; Panda et al. 2004b). cant for neutron stars, and the simplest models are found It is possible that the interior of a neutron-like star by solving the Tolman–Oppenheimer–Volkoff equations does not consist primarily of neutrons, but rather of the (Tolman 1939; Oppenheimer & Volkoff 1939), which are strange matter. Strange matter is composed of deconfined derived from Einstein’s equations in the Schwarzschild quarks, including up, down, and strange quarks, plus the metric for a static, spherical star composed of an ideal leptons necessary to ensure charge neutrality (Bodmer gas. The assumption that the neutrons in a neutron star 1971; Witten 1984). This possibility arises because at the can be treated as an ideal gas is not well justified: The high densities present in the interior of neutron stars, a effect of the strong force needs to be taken into account phase transition from hadronic to quark phase is possi- by replacing the equation of state (EOS) for an ideal gas ble (Glendenning 2000; Menezes & Providência 2003, by a more realistic EOS. Neutron stars can be regarded 2004b, 2004c). Approximately one-third of the quarks are as laboratories for testing different hypotheses relating strange quarks, with the exact fraction depending on the to the EOS through their influence on the properties of model used to describe the quark phase. For example, the neutron star models, including radius, mass, central den- strangeness content is slightly lower for the Nambu–Jona- sity, and Kepler frequency. The effect of many possible Lasinio model (Nambu & Jona-Lasinio 1961a, 1961b) EOS has been studied in this context (Glendenning 2000; than for the MIT bag model (Chodos et al. 1974), as can be Prakash et al. 1997). One class of EOS is based on the seen in Menezes & Providência (2003). It has been argued assumption that the neutron star is composed only of (Alcock et al. 1986) that strange matter is the true ground hadrons (Espíndola & Menezes 2002; Santos & Menezes state of all matter, and if this is the case then as soon as 2004) plus an essential small admixture of electrons. In the core of the star converts to the quark phase, the entire order to construct appropriate EOS for these hadronic star converts. This led to the suggestion that there may be stars, one must rely on models which describe nuclear no neutron stars, and that all neutron-like stars are in fact matter’s bulk properties. Another possibility is that the strange stars (Alcock et al. 1986). interior of the neutron star includes quarks; such mod- Further evidence in favour of strange stars is the els are known as quark stars (Ivanenko & Kurdgelaidze fact that some stars seem to rotate faster that what 1969) or hybrid stars (Glendenning 2000). In hybrid stars would be expected for a neutron star (Glendenning 2000; low-density regions are composed of hadronic matter, but Glendenning & Weber 1992). Neutron stars are gravita- in high-density regions a deconfinement of the quarks tionally bound with the pressure approaching zero as the from the hadrons occurs, leading to a quark phase. Many density approaches zero. Strange stars, on the other hand, different EOS have been built, including EOS for the are self-bound by the strong interaction and gravity is not phase immediately after the formation of the neutron star necessary for their stability. In this case the pressure goes
© Astronomical Society of Australia 2005 10.1071/AS05022 1323-3580/05/04292
Downloaded from https://www.cambridge.org/core. IP address: 170.106.35.76, on 01 Oct 2021 at 02:18:19, subject to the Cambridge Core terms of use, available at https://www.cambridge.org/core/terms. https://doi.org/10.1071/AS05022 Strange Star Equations of State Revisited 293
to zero at a finite density, so that there is a sharp change in an ideal Fermi gas. The energy density, the pressure, and density at the surface of a strange star. The restriction on the quark q density, respectively, are given by the rotational speed of a gravitationally bound star, such d3p that it does not break up due to the centrifugal accelera- ε = 3 × 2 p2 + m2 f + + f − + Bag (2π)3 q q q tion exceeding the gravitational acceleration at the equator q=u,d,s of the star, might not apply to a strange star and, in this (1) case, strange stars could rotate faster than neutron stars. 1 p4 P = dp f + + f − − Bag (2) It is however worth emphasizing that, although there are π2 q q p2 + m2 predictions that neutron or strange stars can rotate as fast q q as 1.5 kHz, such rapidly rotating objects have never been d3p observed. ρ = 3 × 2 (f + − f −) (3) q (2π)3 q q A strange star must have a thin layer on its surface dom- inated by the electrons, which are necessary to enforce where 3 stands for the number of colours, 2 for the spin charge neutrality. This layer could suspend a hadronic degeneracy, mq for the quark masses, and Bag represents crust, which would not be in contact with the stellar core the bag pressure. The distribution functions for the quarks (Alcock et al. 1986; Glendenning & Weber 1992). It has and anti-quarks are the Fermi–Dirac distributions been suggested that the presence of a crust provides a natu- fq± = 1/(1 + exp[( ∓ µq)/T ]), (4) ral explanation for pulsar glitches (Glendenning & Weber 1992), which are sudden changes in the rotation period of with µq the chemical potential for quarks (upper sign) the pulsar. Other authors have argued against the presence and anti-quarks (lower sign) of type q and energy of a crust (Xu 2003a; Usov 1998, 2001) and that strange = 2 + 2 p mq. In our calculations we assume the masses stars should be bare. Recently, Aksenov et al. (2004) iden- m = m = 5.0 MeV, m = 150.0 MeV,and adopt different tified some possible characteristics of the radiation that u d s values for the bag constant, Bag. would come from hot, bare-strange stars. Although our The foregoing equations apply at arbitrary tempera- main interest is in bare-strange stars, we investigate the ture. For T = 0, there are no antiparticles, the chemical extent to which the stellar properties are influenced by a potential is equal to the Fermi energy, and the distribu- hadronic crust. tion functions for the particles are the usual step functions In this paper we focus on the different possible mod- f + = θ(P2 − p2). For neutron stars, the typical temper- els for strange matter and their implications for models q Fq atures are somewhat lower than the Fermi temperatures, of strange stars. Besides the MIT bag model (Cho- and the approximation T = 0 is normally made. dos et al. 1974), normally used for strange stars, more sophisticated models are quark models with a colour– 2.2 Quark Matter in β-Equilibrium flavour locked phase (CFL; Shovkovy et al. 2003; Son & In a star with quark matter, we must impose both charge Stephanov 2000a, 200b; Buballa & Oertel 2002; Alford neutrality and equilibrium for the weak interactions, et al. 2001; Alford & Reddy 2003) and the Nambu–Jona- referred to as β-equilibrium (Glendenning 2000). In this Lasinio model (Nambu & Jona-Lasinio 1961a, 1961b). work we only consider the stage after deleptonization We consider only the MIT and the CFL phase models when the neutrinos have escaped. In this case the neutrino here. We also discuss the effects of the inclusion of a chemical potential is zero. For matter in β-equilibrium, we crust, not only in strange stars, for which the crust is must add the contribution of the leptons as free Fermi gases not in contact with the interior, but also for hadronic and (electrons and muons) to the energy and pressure. The hybrid stars for which it is in contact with the interior. relations between the chemical potentials of the different We calculate the Kepler frequencies using a prescription particles are given by the β-equilibrium conditions given in Lattimer & Prakash (2004) to test the hypothe- sis that strange stars can rotate faster than neutron stars. µs = µd = µu + µe µe = µµ (5) Finally, we also investigate how changes in the EOS affect the star properties within the temperature range of For charge neutrality, we must impose 10–15 MeV. 1 ρ + ρ = (2ρ − ρ − ρ ) e µ 3 u d s 2 Quark Matter For the electron and muon densities we have In this section we summarize the models we use to describe the properties of quark matter. d3p ρ = 2 (f + − f −) l = e, µ (6) l (2π)3 l l 2.1 MIT Bag Model where the distribution functions for the leptons are given The MIT bag model (Chodos et al. 1974) has been used by substituting q by l in Equation (4), with µl the chemi- extensively to describe quark matter. In its simplest form, cal potential for leptons of type l. At T = 0, Equation (6) the quarks are considered to be free inside a bag and the becomes = 3 π2 thermodynamic properties are derived by treating them as ρl kFl/3 . (7)
Downloaded from https://www.cambridge.org/core. IP address: 170.106.35.76, on 01 Oct 2021 at 02:18:19, subject to the Cambridge Core terms of use, available at https://www.cambridge.org/core/terms. https://doi.org/10.1071/AS05022 294 D. P. Menezes and D. B. Melrose
The pressure for the leptons is In the above expressions , the gap parameter, is taken to be 100 MeV (Alford & Reddy 2003). 1 p4dp P = (f + + f −). The electric charge density carried by the pion conden- l 2 l l (8) 3π 2 l p2 + m sate is given by l 4 2 mπ 2.3 Colour–Flavour Locked Quark Phase Q = fπµ 1 − (15) CFL e µ4 Recently many authors (Shovkovy et al. 2003; Son & e Stephanov 2000a, 2000b; Buballa & Oertel 2002; Alford For the case of strange stars, charge neutrality has to be et al. 2001) have discussed the possibility that the quark enforced. Hence, we take µe in the above expressions so matter is in a colour–superconducting phase, in which that QCFL vanishes. In writing down the thermodynamic quarks near the Fermi surface are paired, forming Cooper potential, we neglect the contribution due to the kaon con- 4 pairs which condense and break the colour gauge symme- densation, which is an effect of order ms and hence is small 2 2 try (Alford 2001).At sufficiently high density the favoured compared with the µq contribution to the thermody- ∼ phase is called CFL, in which quarks of all three colours namic potential for 100 MeV. We note that the CFL and all three flavours are allowed to pair. phase model is a variant of the MIT bag model in which In this section, we study the equation of state taking the quarks are allowed to pair. Hereafter, we refer to the into consideration a CFL quark-paired phase. We treat model where pairing is present as the CFL phase model. the quark matter as a Fermi sea of free quarks with an additional contribution to the pressure arising from the 3 Results formation of the CFL condensates. In order to describe the crust of hadronic, hybrid, and The CFL phase can be described with the thermody- strange stars we use the well known EOS calculated in namic potential (Alford & Reddy 2003) Baym et al. (1971) for very low densities. This is justified because all our codes only run up to subnuclear densities (µ ,µ ) = (µ ) + (µ ,µ ) + (µ ) CFL q e quarks q GB q e l e of the order of 0.01 fm−3 and the crust can, in principle, (9) bear much lower densities. with µ = µ /3 where µ is the neutron chemical poten- q n n Next we investigate the differences arising from dif- tial, and ferent EOS for strange matter. In Figure 1 we plot the 6 ν EOS obtained from the bag model for three values of the (µ ) = p2 p(p− µ ) quarks q 2 d q π 0 bag constants. One can see a clear discontinuity between the hadronic crust and the inner strange matter. This plot 3 ν + p2dp p2 + m2 − µ is similar to the one shown in Glendenning & Weber π2 s q 0 (1992). The differences arising from different bag values 2 2 appear only at reasonably large energy densities but do 3 µq − + Bag (10) affect the properties of the stars. In particular, the larger π2 the bag value, the larger the gap between the crust and −4 = with mu = md = 5 MeV, the core of the strange star. Note the units: 1 fm 3.5178 × 1014 gcm−3 = 3.1616 × 1035 dyne cm−2. m2 In Figure 2, we show the differences in stiffness in the ν = 2µ − µ2 + s (11) q q 3 EOS obtained from the bag and the CFL phase models within the range of energy densities where they are used. where GB(µq,µe) is a contribution from the Goldstone In previous studies of hybrid stars containing a core of bosons arising due to the chiral symmetry breaking in quark matter (Panda et al. 2004b) it was found that stable the CFL phase (Shovkovy et al. 2003; Son & Stephanov stars are sensitive to the Bag value. In a hybrid star with 2000a, 2000b; Buballa & Oertel 2002; Alford et al. 2001; the quark phase described by a CFL phase model, for a Alford & Reddy 2003) given gap constant, , a phase transition to a deconfined 2 2 CFL phase was found possible only for Bag greater than 1 2 2 mπ GB(µq,µe) =− fπµ 1 − (12) a critical value. For lower values pure quark matter stars e 2 1/4 2 µe were found. For = 100 MeV one should have Bag ≥ 1/4 = where 185 MeV.On the other hand, for Bag 190 MeV it was not possible to obtain stable stars with masses equal to (21 − 8ln2)µ2 32 most of known radio-pulsar masses. In the present work f 2 = q m2 = m (m + m ) π 2 π 2 2 s u d 36π π fπ we also vary the value of the bag constant to check its (13) influence on the properties of the strange stars. l(µe) is the negative of Equation (8), and the quark Given the EOS, the next step is to solve the Tolman– number densities are equal, i.e., Oppenheimer–Volkoff equations (Tolman 1939):
3 2 + + π 3 ν + 2 µq dP G (ε P)(M 4 r P) ρ = ρ = ρ = (14) =− (16) u d s π2 dr r (r − 2GM)
Downloaded from https://www.cambridge.org/core. IP address: 170.106.35.76, on 01 Oct 2021 at 02:18:19, subject to the Cambridge Core terms of use, available at https://www.cambridge.org/core/terms. https://doi.org/10.1071/AS05022 Strange Star Equations of State Revisited 295
10 2 1 0.1 1.5 0.01 ) 4