# 1-2-3-flavor color superconductivity in compact stars

###### Abstract

We suggest a scenario where the three light quark flavors are sequentially deconfined under increasing pressure in cold asymmetric nuclear matter as, e.g., in neutron stars. The basis for our analysis is a chiral quark matter model of Nambu–Jona-Lasinio (NJL) type with diquark pairing in the single flavor color-spin-locking (CSL), two-flavor (2SC) and three-flavor color-flavor locking (CFL) channels, and a Dirac-Brueckner Hartree-Fock (DBHF) approach in the nuclear matter sector. We find that nucleon dissociation sets in at about the saturation density, , when the down-quark Fermi sea is populated (d-quark dripline) due to the flavor asymmetry imposed by -equilibrium and charge neutrality. At about u-quarks appear forming a two-flavor color superconducting (2SC) phase, while the s-quark Fermi sea is populated only at still higher baryon density. The hybrid star sequence has a maximum mass of 2.1 M. Two- and three-flavor quark matter phases are found only in gravitationally unstable hybrid star solutions.

Introduction - Recent results from observations of compact star properties provide constraints on the nuclear equation of state (EoS) [1]. In particular, the high masses of compact stars in low-mass X-ray binaries, e.g., 4U 1636-536 [2] and the indicated large radius km of the isolated neutron star RX J1856.5-3754 [3] point to a stiff EoS at high densities. Other limits on the stiffness of the EoS comes from heavy-ion collision data for kaon production and elliptic flow (see [1] for references). A key question is whether the phase transition to quark matter can occur inside compact stars [4] and if it is accompanied by observable signatures. Based on a DBHF approach in the nuclear matter sector and a chiral quark matter model with a vector meanfield interaction a recently developed class of hybrid EoS [5, 6] obtained stable hybrid stars with masses from 1.2 M up to 2.1 M in accordance with modern mass-radius constraints. Under the -equilibrium condition in compact stars the chemical potentials of quarks and electrons are related by and . The mass difference between the strange and the light quark flavors has two important consequences: (1) the down and strange quark densities are different so charge neutrality requires a finite electron density and, consequently, (2) .

When increasing the baryochemical potential, the d-quark chemical potential is the first to reach a critical value where the the partial density of free d-quarks becomes finite in a first order phase transition. Due to the finite value of the u-quark chemical potential is below the critical value and the s-quark density is zero due to the high s-quark mass. A single-flavor d-quark phase therefore forms. We discuss this phase here under the natural assumption that the neutralizing background is provided by nuclear matter, for details see [7].

Phase transition to quark matter: nucleon dissociation - The task to describe the phase transition from nuclear matter to quark matter as a dissociation of three-quark bound states into their quark constituents (Mott transition) has not yet been solved. Only some aspects have been considered yet, e.g., within a nonrelativistic potential model [8] or within the NJL model [9]. We will consider a chemical equilibrium of the type , which results in a mixed phase of nucleons and down quarks once the d-quark chemical potential exceeds the critical value. This scenario is analogous to the dissociation of nuclear clusters in the crust of neutron stars (neutron dripline) and the effect may therefore be called the d-quark dripline. We approximate the quark and nucleon components as subphases described by separate models, e.g., by the DBHF approach and a three-flavor quark model of NJL type [10, 11, 12, 13]. The partition function for the latter is given in path-integral representation by

(1) |

where is the diagonal quark chemical potential matrix and the current-quark mass matrix. For , , otherwise and are Gell-Mann matrices acting in flavor and color spaces, respectively. is the charge conjugation operator and . The coupling strength, , of the scalar quark-antiquark current-current interaction and the 3-momentum cutoff, , are fixed by low-energy QCD phenomenology (see table I of [14]). The relative coupling strengths of the spin-0 and spin-1 diquark currents, and , are essentially free parameters, but we restrict the discussion to the values from a Fierz transformation of the one-gluon exchange interaction, and , see [15]. The mass gaps, pairing gaps and quark matter EoS are obtained from the mean-field thermodynamic potential, , which is obtained from (1). In Fig. 1 the phase diagram and hybrid EoS are illustrated. The phase diagram maps the phases with lowest free energy in the plane of baryon chemical potential, , and electric charge chemical potential, . The EoS is the path between oppositely charged phases. Volume fractions are determined such that the mixture of the two phases is neutral. The equations for the mass gaps and pairing gaps are solved separately for the CSL [15] () and the 2SC/CFL [10] () phases, see [16] for a recent review. Each subphase is locally color neutral, and the mixture is charge neutral and -equilibrated [17]. The discussion of surface tension and charge screening effects goes beyond the present exploratory investigation. It shall be performed along the lines of Refs. [18, 19, 20].

In Fig. 2 we show compact star sequences that correspond to the hybrid EoS described above and a -equilibrated nuclear DBHF EoS. We also plot density profiles of two stars with masses M and M. It has been suggested that matter at high densities can be probed with analyses of the cooling behavior [21, 22, 23, 24] and the stability of rapidly rotating stars against r-modes [25, 26]. Since the mixed d-quark CSL and nuclear matter phase extends up to the crust-core boundary and therefore could have significant consequences, we suggest further investigations to clarify the impact of this new phase using the probing tools mentioned above.

Conclusions - In this contribution we have suggested a new quark-nuclear hybrid equation of state for compact star applications that fulfills modern mass-radius constraints. Due to isospin asymmetry, down-quarks may “drip out” from nucleons and form a single-flavor color superconducting (CSL) phase coexisting with nuclear matter already at the crust-core boundary in compact stars. The CSL phase has interesting cooling and transport properties that are in accordance with constraints from the thermal and rotational evolution of compact stars [27]. It remains to be investigated whether this new compact star composition could lead to unambiguous observational consequences. Eventually, it could contribute to an explanation of the puzzling superburst phenomenon [28].

Acknowledgements - We thank J. Berdermann for providing results for the CSL phase and C. Fuchs for the DBHF EoS. D.B. is supported in part by the Polish Ministry of National Education (MENiSW). F.S. acknowledges support from the Belgian fund for scientific research (FNRS). T.K. is grateful for support from the Department of Energy, Office of Nuclear Physics, contract no. DE-AC02-06CH11357. D.B. and F.S. thank the European Science Foundation (ESF) for support from the CompStar programme.

## References

## References

- [1] Klähn T et al. 2006, Phys. Rev. C 74 035802.
- [2] Barret D, Olive J F and Miller M C 2005, Mon. Not. Roy. Astron. Soc. 361 855.
- [3] Trümper J E, Burwitz V, Haberl F and Zavlin V E 2004, Nucl. Phys. Proc. Suppl. 132 560.
- [4] Alford M et al. 2007, Nature 445 E7.
- [5] Klähn T et al. 2006, Phys. Lett. B 654 170.
- [6] Blaschke D B et al. 2007, Phys. Rev. C 75 065804.
- [7] Blaschke D, Sandin F, Klähn T and Berdermann J 2008, arXiv:0807.0414 [nucl-th].
- [8] Röpke G, Blaschke D and Schulz H 1986, Phys. Rev. D 34 3499.
- [9] Lawley S, Bentz W and Thomas A W 2006, J. Phys. G 32 667.
- [10] Blaschke D et al. 2005, Phys. Rev. D 72 065020.
- [11] Rüster S B et al. 2005, Phys. Rev. D 72 034004.
- [12] Abuki H and Kunihiro T 2006, Nucl. Phys. A 768 118.
- [13] Warringa H J, Boer D and Andersen J O 2005, Phys. Rev. D 72 014015.
- [14] Grigorian H 2007, Phys. Part. Nucl. Lett. 4 382.
- [15] Aguilera D N, Blaschke D, Buballa M and Yudichev V L 2005, Phys. Rev. D 72 034008
- [16] Alford M G, Schmitt A, Rajagopal K and Schäfer T 2007, arXiv:0709.4635 [hep-ph].
- [17] Glendenning N K 1992, Phys. Rev. D 46 1274.
- [18] Alford M G, Rajagopal K, Reddy S and Wilczek F 2001, Phys. Rev. D 64, 074017.
- [19] Reddy S and Rupak G 2005, Phys. Rev. C 71, 025201.
- [20] Voskresensky D N,Yasuhira M and Tatsumi T 2002, Phys. Lett. B 541, 93.
- [21] Blaschke D and Grigorian H 2007, Prog. Part. Nucl. Phys. 59 139.
- [22] Popov S, Grigorian H, Turolla R and Blaschke D 2006, Astron. Astrophys. 448 327.
- [23] Popov S, Grigorian H and Blaschke D 2006, Phys. Rev. C 74 025803.
- [24] Blaschke D, Klähn T and Voskresensky D N 2000, Astrophys. J. 533 406
- [25] Madsen J 1999, Phys. Rev. Lett. 85 10.
- [26] Drago A, Pagliara G and Parenti I 2007, Astrophys. J 678 L117.
- [27] Blaschke D B and Berdermann J 2007, AIP Conf. Proc. 964 290.
- [28] Page D and Cumming A 2005, Astrophys. J. 635 L157.