# Radial oscillations of color superconducting self-bound quark stars

## Abstract

We investigate the effect of the color-flavor locking pairing pattern on the adiabatic radial oscillations of pure self-bound quark stars using an equation of state in the framework of the MIT Bag model. We integrate the equations of relativistic radial oscillations to determine the fundamental and the first excited oscillation modes for several parameterizations of the equation of state. For low mass stars we find that the period of the fundamental mode is typically ms and has a small dependence on the parameters of the equation of state. For large mass stars the effect of color-flavor locking is related to the rise of the maximum mass with increasing . As for unpaired quark stars, the period of the fundamental mode becomes divergent at the maximum mass but now the divergence is shifted to large masses for large values of the pairing gap . As a consequence, the oscillation period is strongly affected by color superconductivity for stars with . We fit the period of the fundamental mode with appropriate analytical functions of the gravitational redshift of the star and the pairing gap . We further discuss the excitation and damping of the modes and their potential detectability during violent transient phenomena.

## I Introduction

Neutron stars’ oscillations have been investigated for over almost 50 years (since the work of Chandrasekhar in 1964 (1)) using various equations of state (EOS) for nuclear matter and quark matter. The first exhaustive compilation of radial modes for various zero temperature EOSs was presented by Glass & Lindblom in 1983 (2). Later studies by Väth & Chanmugan (3) complemented and corrected earlier studies and also considered the case of strange quark stars. Further work by Kokkotas & Ruoff (4) presented a new survey of the radial oscillation modes of neutron stars. Their study includes an extensive list of frequencies for the most common equations of state. Radial pulsations were also considered in the case of newborn neutron stars (protoneutron stars) by Gondek et al. (5). They found that the spectrum of the lowest modes of radial pulsations of proto-neutron stars is quite different from that of cold neutron stars. Generally, protoneutron stars are significantly softer with respect to the radial pulsations, than cold neutron stars, and this difference increases for higher modes and for lower stellar masses (5). These differences stem from different structure of proto-neutron stars, which in contrast to cold neutron stars have extended envelopes, inflated by thermal and trapped neutrinos effects.

Radial oscillations were also studied in the case of quark stars. Stars containing quark phases fall into two main classes: hybrid stars (where quark matter is restricted to the core) and strange stars (made up completely by quark matter). It is expected that both kinds of stars cannot exist simultaneously in Nature, but it is not know which one would be realized (if any). This depends on whether the energy per baryon of -equilibrated quark matter at zero pressure and zero temperature is less than the neutron mass (the so called “absolute stability” condition (6); (7)). Analysis made within the MIT bag model shows that there is a room in the parameter space for the existence of strange stars. Moreover, color superconductivity enlarges substantially the region of the parameter space where -stable quark matter has an energy per baryon smaller than the neutron mass (8); (9). As a consequence, a “color superconducting strange matter” is allowed for the same parameters that would otherwise produce unbound strange matter (8). On the other hand, within a Nambu-Jona-Lasinio (NJL) description of quark matter, the strange matter hypothesis is not favored, at least for the most accepted parameterizations of the EOS (10). Thus, stars containing quark phases are believed to be hybrid stars within the NJL model.

Several works have dealt with the problem of radial oscillations of strange and hybrid stars (3); (11); (12); (13); (14); (16); (15). Väth & Chanmugan (3) showed that the oscillation frequencies of strange stars have qualitatively a different dependence on the central density compared with the case of neutron stars, as the periods of all modes go to zero when the central density of the strange star approaches its smallest possible value. Benvenuto & Horvath (11) also presented calculations of radial oscillations of homogeneous strange stars, using a parameterized form of the equation of state. They showed that the particular form of the equation of state allows some simple and general scaling relations which may prove to be very useful for the search of these objects. Strange stars have been mostly studied in the framework of the MIT bag model, but some work has been made using other microphysical description of quark matter. In particular, Benvenuto and Lugones (12) studied the radial oscillations of strange stars in the Quark-Mass Density-Dependent model (17); (18) (QMDD model). Their results showed that oscillation periods are similar to those obtained within the MIT bag model. Recently, Anand et al. (13) studied the radial oscillations of magnetized strange stars in the QMDD model showing that the squares of the frequencies are always decreasing functions of the central density of the strange star. In a similar way Singh et al. (14) investigated the radial oscillations of rotating strange stars in strong magnetic fields in the QMDD model. They showed that the difference in frequency between rotating and non-rotating stars is larger for higher magnetic fields. The change is small for low mass stars but it increases with the mass of the star. This change of frequency is significant for the most massive stars whereas it is marginal for a 1.4 M star.

In the case of hybrid stars Gupta et al. (15) studied the effect of a mixed quark-nuclear matter core on radial oscillations. They found that the effect of the mixed phase is to decrease the maximum mass of the stable neutron star and to cause a kink in radial oscillation frequencies at the onset of the mixed phase. This kink can be traced to a slight kink in the density profile as well as in the EOS for neutron stars having a mixed quark-hadron phase (15). Their results appears to be robust and independent of the EOS. Gondek and Zdunik (16) also studied radial pulsations of neutron stars and strange quark stars with a nuclear crust. They used neutron star models constructed using a realistic equation of state of dense matter and strange star models using a phenomenological Bag model of quark matter. They calculated the eigenfrequencies of the three lowest modes of linear, adiabatic pulsations and found an avoided crossing phenomenon that is strongly related to the changes of compressibility of the matter throughout the star.

In this work we shall focus on the radial oscillations of strange quark stars paying particular attention to the effect of color superconductivity, which for the best of our knowledge, has not been taken into account yet in the literature. In particular, we shall consider the effect of color-flavor locking within the frame of the MIT bag model. The paper is organized as follows: in Sec. II we present the equation of state of the CFL phase and find simple analytical expressions for the relevant thermodynamical quantities. In Sec. III we study the radial oscillations of CFL strange stars and we calculate the fundamental and first excited modes for different values of the parameters of the EOS. In Sec. IV we discuss our results.

## Ii Thermodynamics of the CFL phase

As mentioned before, we work here within the hypothesis that quark matter is absolutely stable and thus quark stars are completely made up by an almost symmetric mixture of up, down and strange quarks. At sufficiently large densities and low temperatures quark matter is a color superconductor, which is a degenerate Fermi gas of quarks with a condensate of Cooper pairs near the Fermi surface (19). Color superconducting quark matter can come in a multiplicity of different phases, based on different pairing patterns of the quarks. At asymptotically large densities, where the quark masses are negligibly small compared to the quark chemical potential, three-flavor quark matter is in the color-flavor locked (CFL) state (20). In this state quarks form Cooper pairs of different color and flavor where all quarks have the same Fermi momentum and electrons cannot be present (21). Color-flavor locking has a profound effect on the properties of quark matter, mainly on transport properties such as mean free paths, conductivities and viscosities. Concerning the equations of state, the effects enter as a term of order which is of a few percent for typical values of the color superconducting gap ( MeV) and the baryon chemical potential ( MeV). However, the effect is proportionally very large in the low pressure regime that affects the absolute stability of quark matter. Thus, pure self-bound quark-matter stars (strange stars) may exist for a wider range of parameters of the MIT Bag model equation of state (8). This affects considerably the mass-radius relationship of CFL stars, allowing for very large maximum masses (9); (22).

The equation of state for CFL quark matter can be obtained in the framework of the MIT bag model. To order , the thermodynamical potential reads (8)

(1) |

being the thermodynamical potential of a state of unpaired , and quarks in which all them have a common Fermi momentum , with chosen to minimize . The binding energy of the diquark condensate is included by subtracting a condensation term proportional to where the chemical potential is related to through , being the mass of the strange quark. Confinement is introduced through a phenomenological vacuum energy density or bag constant .

From the given above we can obtain the following expressions for the pressure and the energy density to order (8):

(2) |

(3) |

where

(4) |

In order to have the EOS in the form , we can invert Eq. (2) to find as a function of

(5) |

and we can write from Eqs. (2) and (3)

(6) |

If we need we may write

(7) |

with given by

(8) |

For the oscillation equation (see Sec. III) we need the adiabatic index . From Eq. (7) we have

(9) |

The derivative can be calculated from Eq. (8): . Thus, we have

(10) |

which allows to write the adiabatic index as function of using Eq. (5) or as a function of using Eq. (8).

Since the values of , and are not accurately known we shall consider them as free parameters in the equation of state. We emphasize that all the values of , and employed in this paper fall inside the stability windows presented in Fig. 2 of Ref. (8); i.e. we always obtain strange stars when integrating the stellar structure equations. Additionally, the parameters satisfy the stability condition given in Ref. (23).

In our calculations of the next section we shall also employ the MIT equation of state for a gas of unpaired massless quarks, i.e. . Notice that unpaired matter is the same that CFL matter with only if all the quark masses are zero. If is not zero, the Fermi momenta of , and quarks are different for unpaired quark matter, but they are still the same for CFL matter.

## Iii Radial Pulsations of CFL stars

In this section we shall study the radial oscillations of strange quark stars employing the equation of state of the previous section. We shall consider the unperturbed star to be composed of a perfect fluid, whose stress-energy tensor takes the form

(11) |

The generic background space-time of a static spherical star is expressed through the line element

(12) |

The Einstein equations in such a spacetime lead to the following set of stellar structure equations (Tolman-Oppenheimer-Volkoff equations)

(13) | |||

(14) | |||

(15) |

where is the gravitational mass inside the radius . The metric function has the boundary condition

(16) |

where is the radius of the star and its mass. With this condition the metric function will match smoothly to the Schwarzschild metric outside the star.

To obtain the equations that govern radial oscillations, both fluid and spacetime variables are perturbed in such a way that the spherical symmetry of the background body is not violated. These perturbations are inserted into the Einstein equations and into the energy, momentum and baryon number conservation equations and only the first-order terms are retained.

Several forms of the oscillation equation have been presented in the literature. The original form presented by Chandrasekhar (1) constitutes a Sturm-Lioville problem whose solution provides the eigenvalues and the eigenfunctions for the radial perturbations (see eq. (59) of Ref. (1)). The second order oscillation equation given by Chandrasekhar can be split into two first order equations. This has been done by Vath and Chanmugan (3) who derived a set of first order equations for the quantities and (see eqs. (9) and (10) of Ref. (3)). More recently, Gondek et al. (5) obtained a similar set of equations but for the relative radial displacement and the Lagrangian perturbation of the pressure (see eqs. (11) and (12) of Ref. (5)). All these sets of equations are equivalent. In this work, we adopted the equations of Gondek et al. (5), because this system is particularly suitable for numerical applications and the boundary condition at the star’s surface can be obtained by purely physical arguments. Another important advantage of this system of oscillation equations stems from the fact that they do not involve any derivatives of the adiabatic index, . Adopting the system of equations is

(17) |

(18) | |||

where is the eigenfrequency and the quantities and are assumed to have a harmonic time dependence .

To solve equations (17) and (III) one needs two boundary conditions. The condition of regularity at requires that for the coefficient of the term in Eq. (17) must vanish (3); (5); (16). Thus, we have

(19) |

Notice that the eigenfunctions can be normalized in order to have . The surface of the star is determined by the condition that for , one has . This implies that the Lagrangian perturbation in the pressure at the surface is zero. Therefore the second boundary condition is

(20) |

In order to numerically solve the oscillation equations we proceed as follows. First, we integrate the Tolman-Oppenheimer-Volkoff equations for each set of the parameters of the equation of state (, and ) in order to obtain the coefficients of the oscillation equations for a given central pressure. Then we solve the oscillation equations by means of the shooting method: we start the numerical integration of Eqs. (17) and (III) for a trial value of and a given set of initial values of and which satisfy at the center the boundary condition given above. The equations are integrated outwards trying to match the boundary condition at the star’s surface. After each integration, the trial value of is corrected in order to improve the matching of the surface boundary condition until the desired precision is achieved. The discrete values of for which Eq. (20) is satisfied are the eigenfrequencies of the radial perturbations. Our code was able to reproduce the results of Vath & Chanmugam (3) and of Kokkotas & Ruoff (4).

Our results are shown in Figs. 14 for different values of the parameters , and of the equation of state (falling inside the stability windows presented in Fig. 2 of Ref. (8)). In Figs. 1 and 2 we show the period of the fundamental and the first excited modes of CFL strange stars as a function of the stellar mass and the central energy density . The effect of color superconductivity has two aspects. By one hand, for a fixed mass of the star the period of the fundamental mode is smaller as increases. On the other hand, the maximum mass of the star increases significantly for large values of . As a consequence, the oscillation period is largely affected by color superconductivity for stars with masses near the maximum mass. In Figs. 3 and 4 we show the period of the fundamental and the first excited modes as a function of the gravitational redshift at the surface of the compact star. The effect of is similar to the observed in the plots of versus (c.f. Figs. 1 and 2). For comparison, the calculations were also performed for hadronic stars described by the Bethe-Johnson EOS (24) and for quark stars without color superconductivity. For low mass stars there is a large difference between the oscillation periods of hadronic stars and quark stars, as already known from previous calculations without color superconductivity (see lower panels of Figs. 1 and 2). The difference is also large for high mass stars but this is due to the difference between the maximum mass in different models. A similar behavior is found in the plots of versus (see Figs. 3 and 4).

Notice that while other EOS parameters such as or can have similar effects than on the oscillation periods, the effect of color superconductivity is stronger for MeV or larger. In fact, and only enter in the equation of state through the parameter defined in Eq. (4) and they have opposite effects. However, since is multiplied by and by , the effect of is proportionally smaller that the effect of . In the case of , we notice that the period of the fundamental mode of stars with is rather independent of the value of (and also ) but changes up to a factor of in the here-considered range of . It is also interesting to note that, for equations of state of the form , it has been shown that there exists a simple scaling law relating the oscillation periods corresponding to different values of . The scaling law is where the oscillation period corresponds to and corresponds to (see (11) for more details). Thus, for that equation of state, there is shifting of the periods for different values of . In our case, a similar shifting is observed in the left upper panels of Figs. 1 and 2 but it is not governed by a simple scaling law like the one given above because of the presence the condensation term in Eq. (7), which is not a constant but depends on the energy density .

For applications it is convenient to have analytic fittings of the results. The period of the fundamental mode can be expressed as a function of the redshift of the star and the pairing gap as:

(21) |

where is in milliseconds and / (100 MeV). The coefficients are given in Table 1 for different values of the bag constant . This fitting is better than a in the range 0 MeV 150 MeV and . The maximum gravitational redshift of the star as a function of and can be obtained through

(22) |

with in MeV fm and in MeV. The coefficients are given in Table 2. This fitting is better than a for 0 MeV 150 MeV and .

B = 60 MeV fm | B = 80 MeV fm | B = 100 MeV fm | B = 120 MeV fm | |
---|---|---|---|---|

6.04684 | 5.21480 | 4.59254 | 4.12120 | |

-4.39575 | -2.20430 | -1.86426 | -1.55666 | |

5.06630 | 3.81115 | 3.39918 | 3.09221 | |

-9.53275 | -7.82398 | -6.88127 | -6.26093 | |

2.10216 | 1.20291 | 8.61653 | 6.49367 | |

-2.85725 | -1.73731 | -1.82305 | -1.77994 | |

-1.25543 | -1.81158 | -1.85871 | -1.92614 | |

7.81002 | 5.28224 | 4.17669 | 3.44168 | |

-1.59488 | -4.63609 | -3.34299 | -1.82102 | |

8.21011 | 5.2234 | 5.08038 | 4.77180 |

4.60362 | 3.62306 | -8.89339 | 7.93274 | 1.74139 | -5.88166 |

## Iv Discussion

In this work we have presented a study of the radial pulsational properties of strange quark stars, paying particular attention to the effect of color superconductivity. We have shown that the effect of color flavor locking into the oscillation periods is different for low mass stars (with ) and for large mass stars (with ).

For low mass stars the period of the fundamental mode is typically ms. For a fixed , is almost independent on the value of and but decreases up to a factor of as goes from 0 to MeV.

For large mass stars the effect of color flavor locking is related to the rise of the maximum mass with increasing . As for unpaired quark stars, becomes divergent at the maximum mass but now the divergence is shifted to large masses for large values of . As a consequence, the oscillation period is strongly affected by color superconductivity for .

It is well know that the periods of radial oscillations of strange stars behave very differently from those of hadronic stars (3), specially for low mass stars. This difference is amplified by color superconductivity because the oscillation periods tend to be smaller as the pairing gap increases. Also, as for unpaired strange stars, the oscillation periods go to zero when the central energy density moves towards the smallest possible value for which the pressure is zero (see Eq. (6)). This can be understood by noticing that low-mass quark stars are very well described as non-relativistic constant-density spheres, which have oscillation periods proportional to the adiabatic index (for the fundamental mode we have (24)). Since the pressure tends to zero everywhere for a quark star with , the adiabatic index tends to infinity, diverges, and the oscillation period tends to zero. This behavior is clear in the right upper panels of Figs. 1 and 2.

In addition to knowing the spectrum of pulsations for a given internal composition, it is important to realize whether the resulting modes are able to survive for a sufficiently long time provided they are excited in any astrophysically realistic situation. In principle, radial pulsations could be excited during a variety of catastrophic events taking place during the compact star’s life. Starquakes and stellar collisions are potential mechanisms if enough energy is transferred to vibrational modes. However, it has been shown that radial pulsations of quark stars are quickly damped due to the enormous bulk viscosity of hot unpaired quark matter (25); (26). The volume oscillation forces the system out of chemical equilibrium with respect to the non-leptonic process and the semi-leptonic processes and . In general, semi-leptonic processes are slower than the non-leptonic one, and the most effective damping reaction in unpaired quark matter is (25); (27). Since and quarks are essentially massless while MeV, the forward rate of cannot keep equal to the reverse rate and the system cannot keep in equilibrium during oscillations. This leads to irreversible processes and damping in unpaired quark matter. For a typical stellar oscillation time of 10 s, high amplitude oscillations are damped in fractions of a second due to the non-leptonic process (28). As a consequence, the detection of any signal related to radial stellar vibrations looks unlikely if quark matter is unpaired.

Still, in the CFL phase the contributions to the bulk viscosity from the above processes are exponentially suppressed (29). The thermodynamic and hydrodynamic properties are rather determined by the massless superfluid phonons and thermally excited light pseudo-Nambu-Goldstone bosons. The contribution to bulk viscosity from phonons alone () has been calculated in (30) (see also (31) for a more recent analysis). For MeV and a typical oscillation period ms, the resulting transport coefficient is several orders of magnitude smaller than for unpaired quark matter.

However, depending on the poorly known value for , the dominant contribution to the bulk viscosity may come from weak equilibrium processes involving the neutral kaon and the bosons , e.g. and (32); (29). For oscillations with a timescale of milliseconds and temperatures above a few MeV, the bulk viscosity can become larger than for unpaired quark matter (see (29) and references therein). However, is much less than that of unpaired quark matter for low temperatures. For example, for MeV and MeV, is more than six orders of magnitude smaller than for unpaired matter. Since the dissipation time scale is (26); (28) we expect that radial oscillations of CFL stars would last for several seconds and are potentially observable after a catastrophic event involving the quark star.

Additionally, we must keep in mind that the spherical symmetry is broken if the star is subject to rotation causing a non-zero varying quadrupole moment that generates gravitational radiation (33). This energy sink has a characteristic timescale P yr, where P is the period of rotation in seconds. Thus, while the persistence of primordial radial pulsations should not be expected, some fingerprints of radial oscillation modes could emerge in future observations of violent transient phenomena.

## V Acknowledgements

C. Vásquez Flores acknowledges the financial support received from UFABC and FAPESP. G. Lugones acknowledges the financial support received from FAPESP and CNPq.

### References

- S. Chandrasekhar, Phys. Rev. Lett. 12, 114 (1964).
- E. N. Glass & L. Lindblom, Astrophys. J. Supp 53, 93 (1983).
- H. M. Väth, & G. Chanmugam, A&A 260, 250 (1992).
- K. D. Kokkotas & J. Ruoff, Astronomy and Astrophysics 366, 565 (2001).
- D. Gondek, P. Haensel & J. L. Zdunik, A&A 325, 217 (1997).
- E. Witten, Phys. Rev. D 30, 272 (1984).
- E. Farhi and R.L. Jaffe, Phys. Rev. D 30, 2379 (1984).
- G. Lugones & J. E. Horvath, Phys. Rev. D 66, 074107 (2002).
- G. Lugones and J. E. Horvath, Astronomy and Astrophysics 403, 173 (2003).
- M. Buballa, Phys. Rept. 407, 205 (2005).
- O. G. Benvenuto and J. E. Horvath, Mon. Not. Roy. Astron. Soc., vol. 250, 679 (1991).
- O. G. Benvenuto & G. Lugones, Int. Jour. Mod. Phys. D 7, 29 (1998).
- J.D. Anand, N.C. Devi, V.K. Gupta & S. Singh,The Astrophysical Journal538,870 (2000).
- S. Singh, N. Chandrika Devi, V. K. Gupta, Asha Gupta and J. D. Anand, J. Phys. G: Nucl. Part. Phys. 28, 2525 (2002).
- V. K. Gupta, V. Tuli & A. Goyal, Astrophys. J. 579, 374 (2002).
- D. Gondek & J.L. Zdunik, A&A 344,117 (1999).
- O. G. Benvenuto & G. Lugones, Phys. Rev. D 51, 1989 (1995).
- O. G. Benvenuto & G. Lugones, Phys. Rev. D 52, 1276 (1995).
- Bailin D. and Love A., Phys. Rep. 107, 325 (1984).
- M. G. Alford, K. Rajagopal, and F. Wilczek, Nucl. Phys. B537, 443â458 (1999).
- K. Rajagopal & F. Wilczeck, Phys. Rev. Lett. 86, 3492 (2001).
- J. E. Horvath and G. Lugones, Astronomy and Astrophysics 422, L1 (2004).
- M. Alford, C. Kouvaris, K. Rajagopal, Phys. Rev. Lett. 92, 222001 (2004)
- S. L. Shapiro & S. A. Teukolsky, Black Holes, White Dwarfs and Neutron Stars: The Physics of Compact Objects. Wiley Eds.(1983).
- Q. D. Wang & T. Lu, Phys. Lett B 148, 211 (1984).
- R. Sawyer, Phys. Lett B 233, 412 (1989).
- J. Madsen, Phys. Rev. D47, 325 (1993).
- J. Madsen, Phys. Rev. D 46, 3290 (1992).
- M. G. Alford and A. Schmitt, AIP Conf. Proc. 964, 256 (2007).
- C. Manuel and F. J. Llanes-Estrada, J. Cosmol. Astropart. Phys. 08, 001 (2007).
- M. A. Escobedo, M. Mannarelli and C. Manuel, Phys. Rev. A 79, 063623 (2009)
- M. G. Alford, M. Braby, S. Reddy, and T. Schäfer, Phys. Rev. C 75, 055209 (2007).
- W. Chau, Astrophys. J. 147, 665 (1967).