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

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

David Blaschke, Fredrik Sandin and Thomas Klähn Institute for Theoretical Physics, University of Wroclaw, 50-204 Wroclaw, Poland
Bogoliubov Lab. for Theoretical Physics, JINR Dubna, 141980 Dubna, Russia
IFPA, Département AGO, Université de Liège, 4000 Liège, Belgium
Physics Division, Argonne National Laboratory, Argonne, IL 60439, USA

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].

Figure 1: Left panel: Phase diagram in the plane of baryon and charge chemical potential. The dash-dotted line denote the border between oppositely charged phases. Right panel: Pressure of matter in -equilibrium as a function of the density.

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


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].

Figure 2: Left panel: Compact star sequences. The phase structure of the core changes with increasing density, as indicated in the figure. Constraints on the mass come from 4U 1636 [2] and on the mass-radius relation from RX J1856 [3]. Right panel: Density profiles of two stars with masses  M and  M.

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.



  • [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.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description