text:
# Absence of two-flavor color superconductivity in compact stars

## 1 Introduction

We are beginning to learn some interesting things about matter at very high densities. For experimental information, we have to rely on observations of compact stars, which consist of matter at nuclear density at the surface, increasing to an unknown maximum density at the core. So far these observations have not put strong constraints on the behavior of matter beyond nuclear density. Theoretically, several different possible phases have been conjectured for hyperdense matter of the kind we would expect to find in compact stars. Some of the most interesting possibilities are the color-superconducting phases expected if the density is high enough to transform hadrons into quark matter . Color superconductivity occurs because QCD predicts an attractive force between quarks that are antisymmetric in color, so we expect quarks near their Fermi surfaces to form Cooper pairs, which condense, breaking the color gauge symmetry. One pairing pattern, known to be favored at sufficiently high density, is the color-flavor locked (CFL) phase in which up-down, down-strange, and up-strange Cooper pairs all form, allowing quarks of all three colors and all three flavors to pair . At lower densities, where the strange quark mass disfavors pairing of the strange with the light quarks, previous work has suggested that the favored phase is a form of color superconductivity in which up and down quarks of two of the three colors participate in pairing. This “2SC” phase leaves the third color and strange quarks unpaired.111 The quark matter that we discuss in this paper always contains strange quarks, which may or may not pair with the light quarks. The phase where they do not, which we here call “2SC”, therefore contains unpaired strange quarks, and is the phase we have called “2SC+s” elsewhere .

In this paper we investigate what forms of superconducting quark matter are most likely to occur in compact stars, taking into account the constraints that must be imposed for a stable bulk phase. These are:

* Electromagnetic neutrality. Bulk matter must be neutral with respect to any Abelian gauge charge (broken or unbroken). If it were not, the energy would be nonextensive, growing faster than the volume as the size of the sample increases.
* Color neutrality. A macroscopic chunk of quark matter must be a color singlet. Color neutrality, the equality in the number density of red, green and blue color charges, is a prerequisite for this. Indeed, the authors of Ref.  have shown that so long as a macroscopic chunk of color superconductor is color neutral, carrying out the projection which imposes color singlet-ness has negligible effect on the free energy of the state. This result is analogous to the familiar fact from ordinary superconductivity that the projection which turns a BCS state, wherein particle number is formally indefinite, into a state with definite but very large particle number has no significant effect. Here, “macroscopic” means that the chunk of quark matter must be much larger than the inverse of the gap $`\mathrm{\Delta }`$ , which is of order $`10`$ to $`100`$ MeV in a color superconductor . We are interested in color superconducting regions of order kilometers in size. Thus to very good approximation all we need to worry about is color neutrality. Color singlet-ness follows without paying any further free energy price.

We shall see that when we impose these constraints, the free energy of the 2SC phase becomes large enough that it is unlikely to be found in nature. The CFL phase, by contrast, satisfies the neutrality constraints at a much smaller cost. (The free energy cost of imposing neutrality turns out to be less in the CFL phase than in unpaired quark matter.)

Early work on the 2SC phase ignored issues of color and electromagnetic neutrality, focussing instead on foundational issues like the size of the pairing gap. Calculations were performed in various toy models in which QCD with only two flavors of quarks was analyzed with the assumption of equal chemical potentials and equal number densities for up and down quarks. Subsequently, many authors have assumed that there is some range of values of the strange quark mass for which up-strange and down-strange pairing can be neglected, and so a stable bulk 2SC phase would occur. We argue here that this is not the case. If the strange quark mass precludes up-strange and down-strange pairing, then up-down pairing (i.e. the 2SC phase) is also precluded. Other authors have looked at the implications of maintaining electric and color neutrality in the 2SC phase , but the very existence of stable bulk 2SC quark matter has not previously been questioned.

In this paper we shall consider matter containing up ($`u`$), down ($`d`$) and strange ($`s`$) quarks, but no heavier quarks. We will take the $`u`$ and $`d`$ to be massless, and the strange quark to have an effective (“constituent”) mass $`M_s`$, which may depend on the density and phase, and is larger than its current mass $`m_s120\mathrm{MeV}`$. We shall assume that the quark number chemical potential $`\mu `$ is large compared to $`M_s`$, and work only to lowest nontrivial order in $`M_s^2/\mu ^2`$. This simplifies our analysis considerably, but we should remark that it is unlikely that $`\mu `$ is much larger than 500 MeV even in the center of a compact star, meaning that $`M_s^2/\mu ^2`$ may not be very small.

## 2 Enforcing neutrality

We shall limit our discussion to quark matter in compact stars older than a few minutes. This justifies our working at zero temperature, since the temperature is by then well below the expected value of all gaps and quasiparticle masses. It also means that there has been plenty of time to come to equilibrium under the weak interactions,222Other authors have considered color superconductivity in contexts where the weak interactions are not in equilibrium, as would for example be appropriate if it were somehow possible to realize color superconductivity in a heavy ion collision . so flavor symmetries are broken. It also means that the star is transparent to neutrinos, which leave the system, so lepton number is not conserved. The relevant symmetries are therefore just the color and electromagnetic gauge symmetries and the global symmetry related to baryon number conservation:

$$[SU(3)_{\mathrm{color}}]\times [U(1)_Q]\times U(1)_B.$$
(2.1)
In any color superconducting phase, this symmetry is broken down to some subgroup by condensation of quark pairs.

The baryon number of an isolated macroscopic body is fixed, but since baryon number is not a gauge symmetry, nothing prevents matter from having a constant nonzero baryon density in the infinite volume limit. There is no requirement for “baryon-number neutrality”, and we can treat the chemical potential for baryon number, $`3\mu `$, as a free parameter. A color superconductor is a BCS state in which baryon number is spontaneously broken by virtue of a condensate wherein $`qq0`$. A macroscopic color superconductor may therefore seem to have ill-defined baryon number, but just as in an ordinary superconductor or superfluid, and as shown explicitly in Ref. , projecting the BCS state onto a state of fixed large baryon number has no significant effect.

Because we are only concerned with enforcing color neutrality, as opposed to color singlet-ness, we need only consider the the $`U(1)_3\times U(1)_8`$ subgroup of the color gauge symmetry generated by the Cartan subalgebra $`T_3=\mathrm{diag}(\frac{1}{2},\frac{1}{2},0)`$ and $`T_8=\mathrm{diag}(\frac{1}{3},\frac{1}{3},\frac{2}{3})`$ in color space, where we choose to label the colors as $`(r,g,b)`$. We introduce chemical (color-electrostatic) potentials $`\mu _3`$ and $`\mu _8`$ coupled to the color charges $`T_3`$ and $`T_8`$. Enforcing color neutrality means choosing $`\mu _3`$ and $`\mu _8`$ such that the $`T_3`$ and $`T_8`$ densities vanish. Enforcing $`T_3`$ and $`T_8`$ neutrality suffices to enforce equality in the number of red, green and blue quarks.

Electromagnetism is generated by $`Q=\mathrm{diag}(\frac{2}{3},\frac{1}{3},\frac{1}{3})`$ in flavor space, where we order the flavors as $`(u,d,s)`$. We will focus on the negative charge $`Q_e=Q`$, coupling it to the (negative) electrostatic potential $`\mu _e`$. ($`\mu _e>0`$ corresponds to a density of electrons; $`\mu _e<0`$ to positrons.) Enforcing electric neutrality means choosing $`\mu _e`$ so that $`Q_e`$ vanishes.

In a color superconducting phase, there is an expectation value for some diquark operator. Since a diquark cannot be electrically neutral and cannot be color neutral, some subgroup of $`U(1)_3\times U(1)_8\times U(1)_Q`$ is spontaneously broken. As we shall see, however, there is typically at least one linear combination of $`T_3`$, $`T_8`$ and $`Q`$ with respect to which the condensate is neutral, meaning that there is at least one $`U(1)`$ subgroup of $`U(1)_3\times U(1)_8\times U(1)_Q`$ that the condensate leaves unbroken. We now discuss unbroken and broken gauged charges in turn.

Unbroken gauged charges. Any uniform phase (we will not consider mixed phases, although they may be important in some contexts), must be neutral under unbroken gauged charges, to avoid the infrared-divergent energy cost of long-range electric fields. The corresponding chemical (i.e. electrostatic) potential is forced to the value of $`\mu _Q`$ that solves

$$Q=\frac{\mathrm{\Omega }}{\mu _Q}=0.$$
(2.2)
Choosing this value of $`\mu _Q`$ changes the contribution of any particles with $`Q`$-charge to the free energy of the system. The ensuing free energy cost is proportional to the volume, though, whereas the free energy cost of the electric fields that would arise if neutrality were not enforced grows with system size faster than the volume.

A macroscopic but finite chunk of matter, that has $`\mu _Q0`$ within it chosen so that it is neutral, has a $`Q`$-electric field (i.e. a gradient in $`\mu _Q`$) across its surface. This $`Q`$-electric field corresponds to a layer of polarization. (For analogous surface phenomena occurring at the interface between nuclear and quark matter, see Ref. .) We will not consider surfaces in this paper. For our purposes, it suffices to choose $`\mu _Q`$ to enforce $`Q`$-neutrality, thus obtaining matter that can sensibly be analyzed in the infinite volume limit.

Broken gauged charges. The argument requiring that any uniform phase be neutral applies to broken symmetries too. If a finite-sized sample has a net charge then there will be electric fields outside the sample that grow in strength with the size of the sample. This means that the sample must be neutral to have a good large-volume limit.

Thus for the purposes of imposing neutrality we can treat $`T_3`$, $`T_8`$ and $`Q`$ on the same footing, regardless of which linear combinations are broken and which are unbroken.

## 3 Free energy comparisons

In this section, we compute the free energy of noninteracting (and thus unpaired) quark matter, CFL quark matter, 2SC quark matter, and a variant of 2SC quark matter in which $`u`$ and $`s`$ quarks (rather than $`u`$ and $`d`$ quarks) pair. We work to order $`(M_s/\mu )^4`$ and $`(\mathrm{\Delta }/\mu )^2`$. We shall see that the free energies of these different phases must be compared in a regime where $`\mathrm{\Delta }M_s^2/\mu `$, meaning that an expansion of this sort is formally consistent. The analysis would be considerably more difficult were we to attempt to work to higher order in $`M_s`$, but this will not be necessary to make the qualitative points we wish to make. Quantitatively, we are interested in $`\mu 500`$ MeV and it is likely that $`120<M_s<500\mathrm{MeV}`$ while $`\mathrm{\Delta }10`$ to $`100\mathrm{MeV}`$.

We shall work in unitary gauge, where the color-direction of the quark pairing is chosen to be position-independent. It is always possible to make gauge transformations that make the color-direction of the quarks vary with position and transfer some color into the off-diagonal glue fields. This would change nothing, although it would make the analysis less transparent.

In our analysis, we do not include any contribution to the free energy from the nonabelian gauge bosons. In principle, these can carry charge and so one may worry that they contribute to maintaining (or upsetting) color neutrality. In unitary gauge, the color current carried by gluons arises from the gauge-connection part of the covariant derivative of the gauge field, $`J_{\mathrm{glue}}^\nu =g[A_\mu ,F^{\mu \nu }]`$. In order for the charge $`J_{\mathrm{glue}}^0`$ to be nonzero, there would have to be a nonzero color electric field $`F^{i0}`$, which is impossible in a conducting or broken (superconducting) phase.

### 3.1 Unpaired quark matter

We begin with noninteracting quark matter. If there really were no interactions, and therefore global symmetries only, there would be no reason for the quark matter to be neutral in any sense. So, we imagine turning on arbitrarily small gauge couplings, as this motivates the requirement that we impose electric and color neutrality and satisfy weak equilibrium, and ask what chemical potentials we must introduce in order to achieve neutrality. Given that quark masses are independent of color, we simply set $`\mu _3=\mu _8=0`$ and obtain quark matter in which the up, down and strange quarks are all separately color neutral. Imposing electrical neutrality is (a little) more challenging. Weak equilibrium relates the chemical potentials for the three flavors to just $`\mu `$ and $`\mu _e`$:

$`\mu _u`$ $`=`$ $`\mu \frac{2}{3}\mu _e`$
$`\mu _d`$ $`=`$ $`\mu +\frac{1}{3}\mu _e`$
$`\mu _s`$ $`=`$ $`\mu +\frac{1}{3}\mu _e,`$ (3.3)
independent of color because $`\mu _3=\mu _8=0`$. The free energy is minimized by filling Fermi seas for each quark flavor, up to Fermi momenta given by:

$`p_F^u`$ $`=`$ $`\mu _u`$
$`p_F^d`$ $`=`$ $`\mu _d`$
$`p_F^s`$ $`=`$ $`\sqrt{\mu _s^2M_s^2},`$ (3.4)
since the energy of a strange quark with momentum $`p_F`$ is $`\sqrt{p_F^2+M_s^2}`$. The free energy of unpaired quark matter is then given by

$$\begin{array}{ccc}\hfill \mathrm{\Omega }_{\mathrm{unpaired}}(\mu ,\mu _e,M_s)& =& \frac{3}{\pi ^2}_0^{p_F^u}(p\mu _u)p^2𝑑p+\frac{3}{\pi ^2}_0^{p_F^d}(p\mu _d)p^2𝑑p\hfill \\ & +& \frac{3}{\pi ^2}_0^{p_F^s}(\sqrt{p^2+M_s^2}\mu _s)p^2𝑑p+\frac{1}{\pi ^2}_0^{\mu _e}(p\mu _e)p^2𝑑p.\hfill \end{array}$$

The last term, which is just $`\mu _e^4/12\pi ^2`$, is the free energy of the electrons that are present because $`\mu _e>0`$. We now fix $`\mu _e`$ by requiring electrical neutrality

$$\frac{\mathrm{\Omega }_{\mathrm{unpaired}}}{\mu _e}=0.$$
(3.5)
To lowest nontrivial order in $`M_s`$, this yields

$$\mu _e=\frac{M_s^2}{4\mu },$$
(3.6)
meaning that in electrically neutral noninteracting quark matter the Fermi momenta are given by

$`p_F^d`$ $`=`$ $`\mu +{\displaystyle \frac{M_s^2}{12\mu }}=p_F^u+{\displaystyle \frac{M_s^2}{4\mu }}`$
$`p_F^u`$ $`=`$ $`\mu {\displaystyle \frac{M_s^2}{6\mu }}`$
$`p_F^s`$ $`=`$ $`\mu {\displaystyle \frac{5M_s^2}{12\mu }}=p_F^u{\displaystyle \frac{M_s^2}{4\mu }},`$ (3.7)
to lowest order. Once (3.5) is satisfied, $`\mu _e`$ does not occur linearly in $`\mathrm{\Omega }`$. This means that to see the first effects of $`\mu _e`$, we need to work to order $`\mu _e^2M_s^4`$. To this order, upon substituting (3.6) into (3.1) we find

$$\mathrm{\Omega }_{\mathrm{unpaired}}^{\mathrm{neutral}}=\frac{3\mu ^4}{4\pi ^2}+\frac{3M_s^2\mu ^2}{4\pi ^2}\frac{712\mathrm{log}(M_s/2\mu )}{32\pi ^2}M_s^4.$$
(3.8)
Several points are worth noting before we proceed:

* The electronic contribution to $`\mathrm{\Omega }`$ is of order $`\mu _e^4M_s^8`$, and therefore does not appear in (3.8). We therefore neglect it from the beginning in our analyses of paired phases below.
* Below, we will mostly ignore interactions, but we will include the effects of BCS pairing, which gives rise to gaps $`\mathrm{\Delta }`$ for some or all of the quarks. Interactions of course result in other changes to the free energy. For example, the perturbative (in the QCD coupling constant) corrections to the free energy are well-studied . The leading effect is a change in the coefficient of the $`\mu ^4`$ term in the free energy. These effects are the same for all the phases of quark matter we consider, and thus cancel in the comparisons we shall make. (More precisely, any differences between the free energies of the different color superconducting phases introduced by perturbative QCD interactions are perturbative corrections to the $`\mathrm{\Delta }^2\mu ^2`$ effects we consider explicitly.) These effects would matter only if we compared the free energy of any of the quark matter phases we analyze to that of hadronic matter. Similarly, we leave the bag constant out since it also only matters in comparisons between hadronic and quark matter. We shall not attempt such comparisons.
* Note that the effect of the strange quark mass, combined with the requirement of electric neutrality, is to push $`p_F^d`$ up and $`p_F^s`$ down relative to $`p_F^u`$. It is clear that, at some point, this will cause the CFL phase in which all flavors pair with each other to be disfavored relative to unpaired quark matter. The question we seek to answer is whether at that point pairing of up with down remains possible.

### 3.2 Color-flavor locked quark matter

In color-flavor locked quark matter , quarks form a condensate in which

$$q_a^\alpha C\gamma _5q_b^\beta \mathrm{\Delta }\left(ϵ^{\alpha \beta 1}ϵ_{ab1}+ϵ^{\alpha \beta 2}ϵ_{ab2}+ϵ^{\alpha \beta 3}ϵ_{ab3}\right),$$
(3.9)
where $`\alpha ,\beta `$ are color indices $`(r,g,b)`$ and $`a,b`$ are flavor indices $`(u,d,s)`$. This form follows from the QCD interaction, which is attractive in the color-antisymmetric channel, and from requiring that rotational invariance be preserved, which leads to a spin-antisymmetric state. The flavor structure is then forced to be antisymmetric.333We neglect the additional condensate which is symmetric in color and flavor, because although it must be nonzero it is much smaller than the antisymmetric condensate (3.9, and because its presence does not modify the symmetries of CFL quark matter . Each of the three terms in (3.9) satisfy these antisymmetry requirements; summing the three of them allows all nine quarks to pair, maximizing the pairing energy. In this state, $`rd`$ and $`gu`$ quarks pair yielding two quasiparticles with gap $`\mathrm{\Delta }`$, $`rs`$ and $`bu`$ quarks pair yielding two quasiparticles with gap $`\mathrm{\Delta }`$, $`bd`$ and $`gs`$ quarks pair yielding two quasiparticles with gap $`\mathrm{\Delta }`$, and $`ru`$, $`gd`$ and $`bs`$ quarks pair yielding two quasiparticles with gap $`\mathrm{\Delta }`$ and one with gap $`2\mathrm{\Delta }`$. The $`U(1)_3\times U(1)_8\times U(1)_Q`$ symmetry is broken to $`U(1)_{\stackrel{~}{Q}}`$, where $`\stackrel{~}{Q}=QT_3\frac{1}{2}T_8`$ is the generator of the unbroken symmetry.

To order $`\mathrm{\Delta }^2`$, the free energy of the CFL phase can be described as follows. One begins with the free energy of a (fictional) state of unpaired quark matter in which all quarks that are “going to pair” have a common Fermi momentum $`p_F^{\mathrm{common}}`$, with $`p_F^{\mathrm{common}}`$ chosen to minimize the free energy of this fictional unpaired state. The binding energy of the diquark condensate is included by subtracting $`\mathrm{\Delta }^2\mu ^2/4\pi ^2`$ for every quasiparticle with gap $`\mathrm{\Delta }`$. For all the phases we study, $`\mu _e,\mu _3,\mu _8`$ are of order $`M_s^2/\mu `$, so terms such as $`\mu _e\mu \mathrm{\Delta }^2`$, $`\mu _e^2\mathrm{\Delta }^2`$, etc, can be consistently neglected. Thus,

$$\mathrm{\Omega }_{\mathrm{CFL}}=\frac{1}{\pi ^2}\underset{i=1}{\overset{9}{}}_0^{p_{F,i}^{\mathrm{common}}}\left(\sqrt{p^2+M_i^2}\mu _i\right)p^2𝑑p\frac{3}{\pi ^2}\mathrm{\Delta }^2\mu ^2,$$
(3.10)
where the sum runs over all nine quarks, with $`M_i=0`$ for the up and down quarks and $`M_i=M_s`$ for the strange quarks. The chemical potential $`\mu _i=\mu Q\mu _e+T_3\mu _3+T_8\mu _8`$ for each quark is specified by its electric and color charges. Quarks that pair with each other have the same $`p_{F,i}^{\mathrm{common}}`$: it is in this sense that the parameter is “common”. Since not all nine quarks pair with each other, it appears that there can be four different values of $`p_{F,i}^{\mathrm{common}}`$ for $`(ru,gd,bs)`$, $`(rd,gu)`$, $`(rs,bu)`$, and $`(gs,bd)`$, respectively.444We thank Sanjay Reddy for pointing this out. Each of these four parameters should be determined by minimizing the free energy with respect to it. We shall see below that in neutral CFL quark matter, it turns out that all the $`p_{F,i}^{\mathrm{common}}`$ take on the same value. This form for the free energy of a state with BCS pairing between different species goes back to the work of Clogston and Chandrasekhar , and was derived in the CFL context in Refs. . Note that it is only valid to order $`\mathrm{\Delta }^2`$ and $`M_s^4`$, and that it follows from substituting back into a more general expression the value of $`\mathrm{\Delta }`$ that solves the gap equation; in other words it has already been minimized with respect to $`\mathrm{\Delta }`$ .

We emphasize that, to this order, the nature of the interaction that generates $`\mathrm{\Delta }`$ does not matter. The free energy is given by this prescription regardless of whether the pairing is due to a point-like four-fermion interaction, as in NJL models, or due to the exchange of a gluon, as in QCD at asymptotically high energies . Of course, the strength and form of the interaction determine the value of $`\mathrm{\Delta }`$, which we shall keep as a free parameter.

For massless quarks, evaluating $`p_F^{\mathrm{common}}`$ yields simply the average of the chemical potentials of the quarks that pair. Evaluation of $`p_F^{\mathrm{common}}`$ for pairing between massless and strange quarks yields the result that the lowest order effect of $`M_s`$ on $`p_F^{\mathrm{common}}`$ is to weight the strange quarks in the average as if their chemical potential were depressed by $`M_s^2/2\mu `$. That is, evaluating the parameters $`p_{F,i}^{\mathrm{common}}`$ by minimizing the free energy with respect to them yields

$`p_{F,(ru,gd,bs)}^{\mathrm{common}}`$ $`=`$ $`\mu {\displaystyle \frac{M_s^2}{6\mu }}`$
$`p_{F,(rd,gu)}^{\mathrm{common}}`$ $`=`$ $`\mu \frac{1}{6}\mu _e+\frac{1}{3}\mu _8`$
$`p_{F,(rs,bu)}^{\mathrm{common}}`$ $`=`$ $`\mu \frac{1}{6}\mu _e+\frac{1}{4}\mu _3\frac{1}{6}\mu _8{\displaystyle \frac{M_s^2}{4\mu }}`$
$`p_{F,(gs,bd)}^{\mathrm{common}}`$ $`=`$ $`\mu +\frac{1}{3}\mu _e\frac{1}{4}\mu _3\frac{1}{6}\mu _8{\displaystyle \frac{M_s^2}{4\mu }},`$ (3.11)
to lowest order. This is sufficient to obtain $`\mathrm{\Omega }_{\mathrm{CFL}}`$ to order $`M_s^4`$. Note that the quark number densities in the CFL phase are not simply given by terms proportional to $`(p_F^{\mathrm{common}})^3`$: $`\mathrm{\Omega }_{\mathrm{CFL}}/\mu `$ receives a contribution also from the $`\mathrm{\Delta }^2\mu ^2`$ term. And, of course, in a paired state there is no sharp Fermi surface anyway as the single particle density of states is smeared out by of order $`\mathrm{\Delta }`$. Thus the fictional unpaired quark matter used in the construction of $`\mathrm{\Omega }_{\mathrm{CFL}}`$ is fictional in two senses: (i) the Fermi momenta are not the same as in the unpaired quark matter of Section 3.1; (ii) the CFL state has no Fermi momentum and the parameters $`p_{F,i}^{\mathrm{common}}`$ that arise in its description do not fully specify the quark number density in the CFL state.

To the order at which we are working,

$`\mathrm{\Omega }_{\mathrm{CFL}}={\displaystyle \frac{3\mu ^4}{4\pi ^2}}`$ $`+`$ $`{\displaystyle \frac{3M_s^2\mu ^2}{4\pi ^2}}{\displaystyle \frac{1}{96\pi ^2}}[(16\mu _e^2+12\mu _3^2+16\mu _8^224\mu _3\mu _e16\mu _e\mu _8)\mu ^2`$
$`+`$ $`(16\mu _88\mu _e)M_s^2\mu +7M_s^436M_s^4\mathrm{log}(M_s/2\mu )]{\displaystyle \frac{3\mathrm{\Delta }^2\mu ^2}{\pi ^2}}.`$
(3.12)
Hence,

$$\begin{array}{ccc}\hfill \frac{\mathrm{\Omega }_{\mathrm{CFL}}}{\mu _e}& =& \frac{1}{12\pi ^2}\left(M_s^2\mu +3\mu _3\mu ^2+2\mu _8\mu ^24\mu _e\mu ^2\right)\hfill \\ \hfill \frac{\mathrm{\Omega }_{\mathrm{CFL}}}{\mu _3}& =& \frac{\mu ^2}{4\pi ^2}(\mu _e\mu _3)\hfill \\ \hfill \frac{\mathrm{\Omega }_{\mathrm{CFL}}}{\mu _8}& =& \frac{1}{6\pi ^2}\left(M_s^2\mu 2\mu _8\mu ^2+\mu _e\mu ^2\right).\hfill \end{array}$$
(3.13)
In discussing the neutrality of the CFL phase, it is convenient to use the basis of broken ($`X,Y`$) and unbroken ($`\stackrel{~}{Q}`$) charges,

$$\begin{array}{cccccc}\hfill \stackrel{~}{Q}& =& QT_3\frac{1}{2}T_8\hfill & \hfill \mu _{\stackrel{~}{Q}}& =& \frac{4}{9}(\mu _Q\mu _3\frac{1}{2}\mu _8)\hfill \\ \hfill X& =& Q+T_34T_8\hfill & \hfill \mu _X& =& \frac{1}{18}(\mu _Q+\mu _34\mu _8)\hfill \\ \hfill Y& =& QT_3\hfill & \hfill \mu _Y& =& \frac{1}{2}(\mu _Q\mu _3),\hfill \end{array}$$
(3.14)
where $`\mu _Q=\mu _e`$. Then, we see from (3.13) that, to the order we are working, electric and color neutrality require

$`\mu _X`$ $`=`$ $`{\displaystyle \frac{M_s^2}{9\mu }}`$
$`\mu _Y`$ $`=`$ $`0.`$ (3.15)
The role of the nonzero $`\mu _X`$ is to enforce $`X`$-neutrality, and it does so by forcing the apparently different $`p_{F,i}^{\mathrm{common}}`$’s in (3.11) to all be equal, as was assumed without proof in Ref. . To this order, we see from (3.13) that $`\mathrm{\Omega }_{\mathrm{CFL}}`$ does not depend on $`\mu _{\stackrel{~}{Q}}`$, so that this combination of chemical potentials cannot be specified by the requirement of neutrality. This means that CFL quark matter is a $`\stackrel{~}{Q}`$-insulator . As we discuss below, this result holds even at higher order. Even though CFL quark matter is a $`\stackrel{~}{Q}`$-insulator, meaning that it is neutral for a range of values of $`\mu _{\stackrel{~}{Q}}`$, in order to obtain neutral matter we must ensure that, once we go to higher order, no electrons are present. Electrons give contributions to $`\mathrm{\Omega }`$ of the form $`\mu _e^4=(\mu _{\stackrel{~}{Q}}+\mu _X+\mu _Y)^4`$, and enforcing neutrality therefore fixes $`\mu _{\stackrel{~}{Q}}=\mu _X`$.555An electron quasiparticle excitation above the CFL ground state should be thought of as carrying $`\stackrel{~}{Q}`$ charge. This can be described by saying that if an electron is injected into a lump of CFL material from outside it, the $`X`$-charge of the electron will be screened by the CFL condensate, leaving the electron excitation propagating within the CFL matter carrying only $`\stackrel{~}{Q}`$-charge. As this example makes clear, though, in order to enforce electrical neutrality it is the total $`Q`$-charge of the electron that matters, regardless of the fact that the electron carries its $`\stackrel{~}{Q}`$ with it while losing its $`X`$ to the condensate. Thus, when we enforce neutrality at higher order by ensuring that no electrons are present, we do so by setting $`\mu _e=0`$ and thus $`\mu _{\stackrel{~}{Q}}=\mu _X0`$.

To order $`M_s^4`$, the free energy for electric and color neutral CFL quark matter is given by

$`\mathrm{\Omega }_{\mathrm{CFL}}^{\mathrm{neutral}}`$ $`=`$ $`{\displaystyle \frac{3\mu ^4}{4\pi ^2}}+{\displaystyle \frac{3M_s^2\mu ^2}{4\pi ^2}}{\displaystyle \frac{112\mathrm{log}(M_s/2\mu )}{32\pi ^2}}M_s^4{\displaystyle \frac{3\mathrm{\Delta }^2\mu ^2}{\pi ^2}}`$ (3.16)
$`=`$ $`\mathrm{\Omega }_{\mathrm{unpaired}}^{\mathrm{neutral}}+{\displaystyle \frac{3M_s^448\mathrm{\Delta }^2\mu ^2}{16\pi ^2}},`$
We conclude that CFL quark matter is favored over unpaired quark matter as long as the interactions are strong enough to generate a gap in the CFL phase satisfying

$$\mathrm{\Delta }>\frac{M_s^2}{4\mu }.$$
(3.17)
Since $`\mathrm{\Delta }`$ is of order $`M_s^2`$ in this criterion, our strategy of evaluating $`\mathrm{\Omega }`$ to order $`M_s^4`$ and to order $`\mathrm{\Delta }^2`$ is consistent.

In the calculation above, we have neglected the fact that when $`M_s0`$, the gaps in the CFL phase are $`M_s`$-dependent and, furthermore, are non-degenerate . For example, the gap for the quasiparticles obtained upon pairing the $`rd`$ and $`gu`$ quarks differs from that obtained upon pairing the $`rs`$ and $`bu`$ quarks. The leading $`M_s`$ dependence of the latter is $`\mathrm{\Delta }(M_s)\mathrm{\Delta }(0)(1cM_s^2/\mu ^2)`$, for some coefficient $`c`$ . These effects are therefore of order $`\mathrm{\Delta }(0)^2M_s^2`$ in the free energy and it is consistent to neglect them. Although their neglect is consistent, they are important because they introduce new $`\mu _X`$ and $`\mu _Y`$ dependence into the free energy. At a high enough order that the differences among gaps matters, the fact that different Cooper pairs in the CFL phase have different $`X`$ and $`Y`$ charges matters. (At this order, the expression (3.10) is itself incomplete.) It is worth noting that the $`M_s`$-dependent differences among gaps cannot introduce any dependence on $`\mu _{\stackrel{~}{Q}}`$ because, by the definition of $`\stackrel{~}{Q}`$, each Cooper pair in the CFL phase has $`\stackrel{~}{Q}=0`$. (For example, $`rd`$ with $`\stackrel{~}{Q}=1`$ pairs only with $`gu`$ with $`\stackrel{~}{Q}=+1`$.) The CFL phase is $`\stackrel{~}{Q}`$-neutral and a $`\stackrel{~}{Q}`$-insulator, even at higher order in $`M_s`$. On the other hand, CFL Cooper pairs do have $`X`$ and $`Y`$ charges. The total $`X`$ and $`Y`$ charge of the condensate vanishes only to the extent that all the gaps are the same. This means that imposing neutrality with respect to the broken gauge charges of the CFL phase will be less trivial at higher order than it is at the order we are working. This warrants investigation.

We have also neglected the possibility that when $`M_s0`$, the left-handed CFL condensate and the right-handed CFL condensate may rotate relative to one another in flavor space . By neglecting this kaon condensation, we have neglected an effect that, if it occurs, lowers the free energy of the CFL phase at order $`M_s^4`$  meaning that in this instance our neglect is not consistent. If we take the quantitative results for the correction to $`\mathrm{\Omega }_{\mathrm{CFL}}`$ due to maximal kaon condensation given in Ref. , we find that the $`4`$ in the denominator of (3.17) increases but remains smaller than $`5`$. We shall see below that in the comparison between CFL and 2SC, there is a much larger uncertainty to be concerned with.

### 3.3 2SC quark matter

In 2SC quark matter, quarks form a condensate in which

$$q_a^\alpha C\gamma _5q_b^\beta \mathrm{\Delta }ϵ^{\alpha \beta 3}ϵ_{ab3}.$$
(3.18)
In this state, $`rd`$ and $`gu`$ quarks pair yielding two quasiparticles with gap $`\mathrm{\Delta }`$, and $`ru`$ and $`gd`$ quarks pair, yielding two quasiparticles with gap $`\mathrm{\Delta }`$. Five quarks are left unpaired. The $`U(1)_3`$ symmetry (and indeed a color $`SU(2)`$ symmetry of which it is a subgroup) is left unbroken, as is $`U(1)_{\stackrel{~}{Q}}`$, where $`\stackrel{~}{Q}=Q\frac{1}{2}T_8`$.

To the order we are working, the free energy of the 2SC phase is

$$\begin{array}{ccc}\hfill \mathrm{\Omega }_{2\mathrm{S}\mathrm{C}}& =& \frac{1}{\pi ^2}\underset{i=1}{\overset{4}{}}_0^{p_F^{\mathrm{common}}}(p\mu _i)p^2𝑑p\hfill \\ & +& \frac{1}{\pi ^2}\underset{i=5}{\overset{9}{}}_0^{\sqrt{\mu _i^2M_i^2}}(\sqrt{p^2+M_i^2}\mu _i)p^2𝑑p\hfill \\ & & \frac{1}{\pi ^2}\mathrm{\Delta }^2\mu ^2,\hfill \end{array}$$
(3.19)
where the first four quarks are those that pair while the last five are those that don’t. This time, $`p_F^{\mathrm{common}}`$ is just the average of the chemical potentials of the quarks that pair:

$$p_F^{\mathrm{common}}=\mu \frac{1}{6}\mu _e+\frac{1}{3}\mu _8.$$
(3.20)
To the order at which we are working,

$$\begin{array}{ccc}\hfill \frac{\mathrm{\Omega }_{2\mathrm{S}\mathrm{C}}}{\mu _e}& =& \frac{M_s^2\mu 2\mu _e\mu ^2}{2\pi ^2},\hfill \\ \hfill \frac{\mathrm{\Omega }_{2\mathrm{S}\mathrm{C}}}{\mu _3}& =& \frac{\mu _3\mu ^2}{2\pi ^2},\hfill \\ \hfill \frac{\mathrm{\Omega }_{2\mathrm{S}\mathrm{C}}}{\mu _8}& =& \frac{2\mu _8\mu ^2}{\pi ^2},\hfill \end{array}$$
(3.21)
meaning that we find electric and color neutrality with $`\mu _3=\mu _8=0`$ and $`\mu _e=M_s^2/2\mu `$. Note that $`Q`$ is a linear combination of the unbroken charge $`\stackrel{~}{Q}`$ and the orthogonal broken charge $`XT_8\frac{1}{2}Q`$. It is more natural to describe the 2SC phase in terms of $`\stackrel{~}{Q}`$ and $`X`$ charges, in which case the neutrality condition becomes $`\mu _{\stackrel{~}{Q}}=2M_s^2/(5\mu )`$ and $`\mu _X=M_s^2/(5\mu )`$. The 2SC phase contains ungapped $`\stackrel{~}{Q}`$-carrying modes, and is therefore a $`\stackrel{~}{Q}`$-conductor.

To order $`M_s^4`$, the free energy for electric and color neutral 2SC quark matter is given by

$`\mathrm{\Omega }_{2\mathrm{S}\mathrm{C}}^{\mathrm{neutral}}`$ $`=`$ $`{\displaystyle \frac{3\mu ^4}{4\pi ^2}}+{\displaystyle \frac{3M_s^2\mu ^2}{4\pi ^2}}{\displaystyle \frac{512\mathrm{log}(M_s/2\mu )}{32\pi ^2}}M_s^4{\displaystyle \frac{\mathrm{\Delta }^2\mu ^2}{\pi ^2}}`$ (3.22)
$`=`$ $`\mathrm{\Omega }_{\mathrm{unpaired}}^{\mathrm{neutral}}+{\displaystyle \frac{M_s^416\mathrm{\Delta }^2\mu ^2}{16\pi ^2}}.`$
### 3.4 2SCus quark matter

There is an obvious variant of the 2SC phase which we must consider, before we make the comparison between 2SC and CFL. In the 2SC phase above, it is the $`u`$ and $`d`$ quarks that pair. Given (3.7), it seems as likely for $`u`$ and $`s`$ to pair. In the resulting 2SCus quark matter, quarks form a condensate in which

$$q_a^\alpha C\gamma _5q_b^\beta \mathrm{\Delta }ϵ^{\alpha \beta 2}ϵ_{ab2}.$$
(3.23)
In this state, $`rs`$ and $`bu`$ quarks pair yielding two quasiparticles with gap $`\mathrm{\Delta }`$, and $`ru`$ and $`bs`$ quarks pair, yielding two quasiparticles with gap $`\mathrm{\Delta }`$. Five quarks are left unpaired. The color $`U(1)`$ symmetry generated by $`\mathrm{diag}(\frac{1}{2},0,\frac{1}{2})=\frac{1}{2}T_3+\frac{3}{4}T_8`$ (and indeed a color $`SU(2)`$ symmetry of which this is a subgroup) is left unbroken, as is $`U(1)_{\stackrel{~}{Q}}`$, where $`\stackrel{~}{Q}=QT_3\frac{1}{2}T_8`$.

To the order we are working, the free energy of the 2SCus phase is

$`\mathrm{\Omega }_{2\mathrm{S}\mathrm{C}\mathrm{u}\mathrm{s}}`$ $`=`$ $`{\displaystyle \frac{1}{\pi ^2}}{\displaystyle \underset{i=1}{\overset{4}{}}}{\displaystyle _0^{p_F^{\mathrm{common}}}}(\sqrt{p^2+M_i^2}\mu _i)p^2𝑑p`$ (3.24)
$`+`$ $`{\displaystyle \frac{1}{\pi ^2}}{\displaystyle \underset{i=5}{\overset{9}{}}}{\displaystyle _0^{\sqrt{\mu _i^2M_i^2}}}(\sqrt{p^2+M_i^2}\mu _i)p^2𝑑p{\displaystyle \frac{1}{\pi ^2}}\mathrm{\Delta }^2\mu ^2,`$
where the first four quarks are those that pair while the last five are those that don’t. This time,

$$p_F^{\mathrm{common}}=\mu \frac{1}{6}\mu _e+\frac{1}{3}\mu _8\frac{M_s^2}{4\mu }.$$
(3.25)
To the order at which we are working,

$`{\displaystyle \frac{\mathrm{\Omega }_{2\mathrm{S}\mathrm{C}\mathrm{u}\mathrm{s}}}{\mu _e}}`$ $`=`$ $`{\displaystyle \frac{\mu _e\mu ^2}{\pi ^2}}`$
$`{\displaystyle \frac{\mathrm{\Omega }_{2\mathrm{S}\mathrm{C}\mathrm{u}\mathrm{s}}}{\mu _3}}`$ $`=`$ $`{\displaystyle \frac{\mu _3\mu ^2}{2\pi ^2}}`$
$`{\displaystyle \frac{\mathrm{\Omega }_{2\mathrm{S}\mathrm{C}\mathrm{u}\mathrm{s}}}{\mu _8}}`$ $`=`$ $`{\displaystyle \frac{2\mu _8\mu ^2}{\pi ^2}},`$ (3.26)
meaning that we find electric and color neutrality with $`\mu _3=\mu _8=\mu _e=0`$. (At higher order, all would become nonzero.) The 2SCus phase is a $`\stackrel{~}{Q}`$-conductor.

To order $`M_s^4`$, the free energy for electric and color neutral 2SCus quark matter is identical to that for electric and color neutral 2SC quark matter, given in (3.22).

### 3.5 2SC vs. CFL Comparison

The first comparison to make is that between all the phases we have discussed (unpaired, CFL, 2SC, 2SCus) upon making the assumption that $`M_s`$ is the same in all phases and $`\mathrm{\Delta }`$ is the same in all the paired phases. This comparison is shown in Fig. 1.

If Fig. 1 were the end of the story, we would conclude that for $`\mathrm{\Delta }<M_s^2/4\mu `$ unpaired quark matter is favored, whereas for $`\mathrm{\Delta }>M_s^2/4\mu `$ CFL quark matter is favored. Precisely at $`\mathrm{\Delta }=M_s^2/4\mu `$, unpaired, CFL, 2SC and 2SCus quark matter all have the same free energy. There is therefore (just barely) no 2SC window. Given the fragility of this conclusion, let us enumerate the effects we know we have left out:

* Effects that are of order $`M_s^6`$ and higher. We have discussed some of these above. We do not know how they all add up, but we can consistently neglect them.
* As mentioned above, we have neglected the possibility of kaon condensation in CFL quark matter. This can only lower the CFL free energy. It therefore works against opening a 2SC window.
* For a given interaction, the gap $`\mathrm{\Delta }`$ is not the same in the 2SC and CFL phases. If the interaction is single-gluon exchange as at asymptotically high densities, the CFL gap is smaller than the 2SC gap by a factor of $`2^{1/3}0.79`$ . If the interaction is a point-like interaction with the quantum numbers of single-gluon exchange, the ratio of the CFL gap to the 2SC gap is $`0.75`$ . The fact that the 2SC gap is somewhat bigger than the CFL gap tends to open up a 2SC window. This effect is comparable in magnitude (and opposite in sign) to that of kaon condensation.
* The biggest effect comes from the fact that $`M_s`$ is not the same in the CFL and 2SC phases. The free energy of all the phases we have discussed includes a common $`3/(4\pi ^2)M_s^2\mu ^2`$ term, which plays no role in the comparison of Fig. 1. However when comparing phases with different $`M_s`$ it will parametrically dominate the $`M_s^4`$ and $`\mathrm{\Delta }^2`$ effects on which we have focussed to this point. We now discuss this effect in more detail.

A full treatment of how the constituent mass $`M_s`$ and gaps $`\mathrm{\Delta }`$ vary between different phases would require analysis of the coupled gap equations in some model Hamiltonian, an exercise that goes beyond the general considerations we employ in this paper. Such an analysis has recently been done for several NJL models by Buballa and Oertel , albeit without enforcing electric and color neutrality. They find that $`M_s`$ is considerably smaller in the CFL phase than in the 2SC phase. They do not consider the 2SCus phase, but it is reasonable to assume that $`M_s`$ in this phase will be comparable to that in the 2SC phase.

In Fig. 2, we illustrate the effect of a relatively small difference in the value of $`M_s`$ between phases. We plot the relative free energy of the CFL, 2SC, and unpaired phases, taking ($`M_s=250\mathrm{MeV}`$, $`\mathrm{\Delta }=30\mathrm{MeV}`$) in the CFL phase and ($`M_s=275\mathrm{MeV}`$, $`\mathrm{\Delta }=40\mathrm{MeV}`$) in the unpaired and 2SC phases. (In the models of Ref. , the change in $`M_s`$ between the CFL and 2SC phases is significantly larger, of order $`30\%`$, and there is no analogue of our unpaired quark matter.) We see in Fig. 2 that the fragility of the conclusion of Fig. 1 is illusory: the 2SC phase is firmly excluded if $`M_s`$ varies as Buballa and Oertel find. Because it is an order $`M_s^2`$ effect, the greatest uncertainty in our conclusions arises from our uncertainty in how much $`M_s`$ changes between different phases. Our analysis has also neglected $`\mu `$-independent contributions to the free energy that arise from the binding energy of the $`\overline{s}s`$ condensate, and so favor the 2SC phase where $`M_s`$ is larger. These are subleading relative to the $`\mu ^2M_s^2`$ term.

## 4 Conclusions

We find that, taking into account the constraints imposed by color neutrality, electrical neutrality, and weak equilibrium, the two-flavor color superconducting (2SC) phase is not expected to be found in compact stars. Our argument is model independent, and is based on an expansion in $`M_s/\mu `$ and $`\mathrm{\Delta }/\mu `$. Since it is possible that in lower density quark matter $`M_s`$ is not much smaller than $`\mu `$, it would be of considerable interest to perform the full calculation for a specific model, solving coupled gap equations under the constraint of electrical and color neutrality, to see whether the 2SC phase vanishes as expected. An NJL model such as that of Buballa and Oertel  would be a natural starting point.

Our results indicate that two qualitative possibilities remain viable for the phase diagram of quark matter that is electrically and color neutral. The first is a single phase transition directly from hadronic matter to color-flavor locked quark matter. The second possibility is that there may be a window of intermediate densities in which we find quark matter that is, at the level of the analysis of this paper, unpaired. That is, as a function of increasing $`\mu `$ we first go from hadronic matter to unpaired quark matter and only at a higher $`\mu `$ make the transition to CFL quark matter.

If there is a window of “unpaired” three-flavor quark matter, there will certainly be pairing in this window: all that we have shown is that there will be no BCS pairing between quarks of different flavors. One possible pattern of pairing is the formation of $`uu`$, $`dd`$ and $`ss`$ condensates. These must be either $`J=1`$ or symmetric in color, and are therefore much smaller than the $`J=0`$ color-antisymmetric condensates we have investigated in this paper. The gaps in these phases may be as large as of order 1 MeV , or could be much smaller . In all the condensates we have discussed to this point, the Cooper pairs are made of quarks with equal and opposite momenta. Another possibility for pairing in the “unpaired” quark matter is crystalline color superconductivity , which involves pairing between quarks whose momenta do not add to zero, as first considered in a condensed matter physics context in Refs. . The unpaired quark matter of (3.7) is susceptible to the formation of a crystalline color superconducting condensate constructed from pairs of quarks with different flavors, both of which have momenta near their respective unpaired Fermi surfaces. Our work demonstrates that the electrically neutral unpaired quark matter of Section 3.1 is the correct starting point for an analysis of crystalline color superconductivity. Crystalline color superconductivity need not coexist with the 2SC phase, as previously thought.

Acknowledgements

We have had many illuminating discussions with Sanjay Reddy. We thank him for pointing out an error in the derivation of $`\mathrm{\Omega }_{\mathrm{CFL}}^{\mathrm{neutral}}`$ in an earlier draft. KR is grateful for the support and hospitality of the Institute for Theoretical Physics at Santa Barbara. The research of KR is supported in part by the DOE under cooperative research agreement DE-FC02-94ER40818 and by the NSF under Grant No. PHY99-07949. The research of MGA is supported in part by the UK PPARC.