Bridging the gap by shaking superfluid matter

Bridging the gap by shaking superfluid matter


In cold compact stars, Cooper pairing between fermions in dense matter leads to the formation of a gap in their excitation spectrum and typically exponentially suppresses transport properties. However, we show here that weak Urca reactions become strongly enhanced and approach their ungapped level when the star undergoes density oscillations of sufficiently large amplitude. We study both the neutrino emissivity and the bulk viscosity due to direct Urca processes in hadronic, hyperonic and quark matter and discuss different superfluid and superconducting pairing patterns.

The dense matter in a compact star has certain transport properties, such as bulk viscosity and neutrino emissivity, that are dominated by beta (weak interaction) equilibration processes. However, this equilibration rate is expected to be exponentially suppressed as the star cools below the critical temperature for superconductivity/superfluidity, and Cooper pairing produces a gap in the relevant fermion excitation spectrum (1); (2). This suppression is relevant both in hadronic matter, where the critical temperature for proton superconductivity and neutron superfluidity is in the range (3) and in exotic phases such as quark matter, where critical temperatures as high as if flavor antisymmetric pairing channels are favored , and as low as when flavor singlet pairing is favored (4).

In this letter we show that the exponential suppression operates only at small amplitude and may be overcome in realistic situations. It is already known that beta-equilibration rates in normal matter are enhanced by suprathermal processes at high amplitude when equilibration is driven by higher order terms in the chemical potential imbalance (5); (6); (7). Here we report on a similar but more dramatic phenomenon in superfluid/superconducting phases of matter, arising from a threshold-like behavior with a rapid increase in available phase space when the typical energy in the equilibration processes approaches the gap.

The mechanism we discuss is generic and can be expected to operate in all situations where large perturbations drive the system out of beta-equilibrium. One possible application are unstable oscillations of rotating compact stars, like f- or r-modes (8). These modes are unstable due to gravitational wave emission and their amplitude grows exponentially until they are saturated by non-linear coupling to damped daughter modes (9) or damped by the supra-thermal bulk viscosity (10). Note, that even at amplitudes well below the level where bulk viscosity alone would saturate such modes, suprathermal effects can noticeably affect the thermal evolution of the star, via viscous heating and suprathermally enhanced neutrino emissivity (6). Another class of scenarios stems from large amplitude oscillations caused by singular events, like star quakes (11) or tidal forces preceding neutron star mergers (12). In this case the suprathermal bulk viscosity should control the amplitude of the oscillation and correspondingly the emitted radiation.

Here we will study the effect of large amplitude oscillations in the case of various gapped phases of dense matter. When, as in nuclear matter, Cooper pairs contain two particles of the same flavor, density oscillations cannot lead to pair breaking. (Thermal pair-breaking has been used (13); (14) to explain the rapid cooling observed in Cas A (15).) Nevertheless, oscillations displace the gapped Fermi seas, and at sufficiently large amplitude can bridge the gap, opening phase space for beta equilibration processes as illustrated in fig. 1. As we will show, this suprathermal effect, which begins at amplitudes far below the level where non-linear hydrodynamic effects, of higher order in the fluid velocity, would arise (5); (7), causes a huge increase in the corresponding rates, effectively bringing them up to their ungapped level.

Figure 1: Schematic illustration of the opening of phase space for Urca interactions in gapped nuclear matter when the deviation from beta equilibrium becomes sufficiently large relative to the pairing gaps ,.

Beta equilibration processes have a finite time scale, so a harmonic local density fluctuation with amplitude around the equilibrium value induces an oscillation in the displacement of the chemical potentials of the degenerate particles from their equilibrium. This oscillation leads both to heat generation (via dissipation due to bulk viscosity) and heat loss (from enhanced neutrino emission). Which of the two effects dominates depends among other factors on the amplitude of the oscillation (6). In (7) we found that over the entire range of physically reasonable density amplitudes the chemical potential oscillation is given by the linear, harmonic relation


where is a susceptibility characterizing the particular form of dense matter.

We are interested in the direct “Urca” processes and mediated by charged current interactions, which occur in hadronic matter at high densities where proton fraction (16). To study all relevant cases within a unified framework we use a notation where represents either an elementary negatively charged quark (down or strange) or a hadron that contains that quark. Similarly stands either for an elementary up quark or a hadron containing it; represents a charged lepton (an electron, or in dense neutron matter a muon) and is the corresponding neutrino. We specialize to Fermi liquid theory, where the general dispersion relation of the matter fields is


where labels the matter species and represents a gap in the particle spectrum arising from superfluidity or (color)-superconductivity. Non-Fermi liquid effects should not play a role in a gapped system (17). We neglect dependence of the gap on energy and momentum, but we include temperature and density dependence . The leptons can to very good approximation be described by a free dispersion relation , . Beta equilibrium with respect to the Urca processes enforces . Since compact stars contain ultradegenerate matter with we can restrict the analysis to leading order in where we find the general results for the net rate and the emissivity


In these expressions is the effective coupling of the fields to the weak current, is a function depending on the density and the equation of state


which is identical for both quantities, and the parenthesis vanishes for a free massless dispersion relation. Aside from the neglected temperature corrections any modification of the dispersion relations due to interactions or masses opens phase space and yields a non-zero result. Eqs. (3) and (4) also contain a threshold function , a characteristic temperature and amplitude dependence, and a modification function given for singlet gaps by the dimensionless integrals


where , , and . These functions reflect the phase space available to Urca reactions. They are normalized so that and are symmetric in the gap parameters . In equilibrium they are pure reduction factors that can become extremely small (1). However, out of equilibrium these modification functions exhibit a suprathermal enhancement for (5); (6); (7) and can then greatly exceed one. They depend on the equation of state only via dimensionless ratios and are the same for all direct Urca processes (1). In the ungapped case they have an analytic polynomial form (6).

The oscillation period of a compact star is much smaller than its evolution time scale, so the relevant quantity is the averaged emissivity which takes the same form eq. (4) with the averaged modification function depending on . The bulk viscosity of dense matter is (7)


Inserting the rate eq. (3) and the low-amplitude form of the chemical potential oscillation eq. (1) gives


with an analogous modification function


In hyperon and quark matter the viscosity is dominated by neutral current rather than Urca processes, but for those we also expect a high-amplitude enhancement (7).

Figure 2: The modification functions eqs. (4) & (8) as a function of oscillation amplitude , for different pairing patterns with maximum gap : symmetric case where only one particle is gapped, e.g. , (solid lines); intermediate cases (long dashed lines); (short dashed lines); symmetric case (dotted lines). Gap ranges from to . Left panel: for the averaged neutrino emissivity. Right panel: eq. (9) for the bulk viscosity.
Figure 3: Enhancement via suprathermal oscillations of the neutrino emissivity and the bulk viscosity of hadronic matter with a proton gap at and a neutron gap at .
or or
for a free gas for a free gas
free gas:
Table 1: Parameters for direct Urca processes in dense matter, with Fermi constant , Cabbibo angle , vector and axial couplings and (1), proton fraction , symmetry energy (16), in-medium hadron masses and quark interaction parameter (18); (7). The quark expression is to leading order in and all results to leading order in , and .

To see the gap-bridging enhancement of beta-equilibration rates for oscillations of suprathermal amplitude, without specifying a particular form of matter, we show in fig. 2 the modification functions and . We vary the maximum gap and the ratio of the two gaps. For strongly gapped particles , the larger gap tends to dominate, e.g. is more suppressed than . Once the amplitude becomes comparable to , the Urca reactions can bridge the gap(s) so that the modification functions rise steeply, and quickly reach their ungapped levels. Interestingly, for matter where the two gaps are the same (dotted lines) this rise starts already for and eventually merges into the solid curve for a single gap . In contrast, in the asymmetric case where one gap is much larger than the other, the effects of the smaller gap are only visible at large amplitudes and become negligible when it is more than an order of magnitude smaller.

Fig. 2 shows that suprathermal enhancement in gapped matter is considerably bigger than in ungapped matter (the line) (7). This enhancement could be realistically achieved in compact star oscillations: is related to the density amplitude by the susceptibility eq. (1). For hadronic matter with an APR equation of state (19), increases from to (2); for quark matter it decreases from to . Thus for amplitudes the amplitude can indeed become large enough to bridge gaps (7).

To apply these general results to explicit processes requires dispersion relations and effective couplings for the relevant particles. Here we consider direct Urca processes in hadronic (16); (20); (2), hyperonic (21); (2) and quark matter (22); (23). Parameter values are given in tab. 1. Inserted in the general results eqs. (4) and (8) they reproduce previous expressions in the subthermal case . In the ideal gas approximation our results simplify to


where is the temperature in units of and the values of and are given in tab. 1.

As an explicit example that takes into account the complications in gapped interacting matter we consider APR hadronic matter (19) with the in-medium hadron masses given in (24). For the density dependence of the proton gap we use a Gaussian fit to data from (3), with maximum centered at . For illustrative purposes we assume the neutrons also have a gap (in reality it is expected to be ) with the density dependence proposed to explain the Cas A cooling data in (14), i.e. with maximum at . The gaps have a BCS temperature dependence, but this has only a modest effect in small regions where . Our results for neutrino emissivity and bulk viscosity at are shown in fig. 3. At infinitesimal amplitudes (thick lines) there is enormous suppression by the proton gap, and moderate suppression by the smaller neutron gap. As the density oscillation amplitude increases, it first, at , overcomes the suppression of emissivity and viscosity due to the neutron gap, increasing them at densities potentially relevant for direct Urca processes by many orders of magnitude. Then, at , even the large proton gap at lower density is bridged, boosting the results by up to ! As in previous results for ungapped matter (7), these effects are expected to be even more pronounced for the realistic case of modified Urca reactions (1); (6), that are actually present at moderate densities.

Fig. 3 shows that, if the star is cold enough that Urca processes are blocked by pairing gaps throughout the star, density oscillations could provide the dominant contribution to those processes, by bridging the pairing gaps, starting in the regions where the rate-controlling gap is the smallest. The smaller this minimum gap, the lower the amplitude required to bridge it. In hyperonic (21); (2) or quark matter (25); (26) there are processes which are only suppressed by MeV, that could be bridged by oscillations with amplitude as small as .

We are grateful to C. Horowitz and J. Read for helpful discussions. This work was initiated during the workshop INT-11-2b: “Astrophysical Transients: Multi-messenger Probes of Nuclear Physics”. We thank the Institute for Nuclear Theory at the University of Washington for its hospitality and the U.S. Department of Energy for support under contracts #DE-FG02-91ER40628, #DE-FG02-05ER41375, #DE-FG02-00ER41132 and by the Topical Collaboration “Neutrinos and nucleosynthesis in hot dense matter”.


  1. D. Yakovlev, A. Kaminker, O. Y. Gnedin, and P. Haensel, Phys.Rept. 354, 1 (2001), astro-ph/0012122.
  2. P. Haensel, K. P. Levenfish, and D. G. Yakovlev, Astron. Astrophys. 357, 1157 (2000), astro-ph/0004183; Astron. Astrophys. 381, 1080 (2002), astro-ph/0110575.
  3. D. J. Dean and M. Hjorth-Jensen, Rev. Mod. Phys. 75, 607 (2003), nucl-th/0210033.
  4. M. G. Alford, A. Schmitt, K. Rajagopal, and T. Schafer, Rev. Mod. Phys. 80, 1455 (2008), 0709.4635.
  5. J. Madsen, Phys. Rev. D46, 3290 (1992).
  6. A. Reisenegger, Astrophys. J. 442, 749 (1995), astro-ph/9410035.
  7. M. G. Alford, S. Mahmoodifar, and K. Schwenzer, J. Phys. G37, 125202 (2010), 1005.3769; AIP Conf. Proc. 1343, 580 (2011), 1012.4834.
  8. N. Stergioulas, Living Rev.Rel. 6, 3 (2003), gr-qc/0302034.
  9. P. Arras et al., Astrophys. J. 591, 1129 (2003), astro-ph/0202345.
  10. M. G. Alford, S. Mahmoodifar, and K. Schwenzer, (2011), 1103.3521
  11. L. M. Franco, B. Link, and R. I. Epstein, Astrophys. J. 543, 987 (2000), arXiv:astro-ph/9911105.
  12. D. Tsang, J. S. Read, T. Hinderer, A. L. Piro, and R. Bondarescu, (2011), 1110.0467.
  13. D. Page, M. Prakash, J. M. Lattimer, and A. W. Steiner, Phys. Rev. Lett. 106, 081101 (2011), 1011.6142.
  14. P. S. Shternin et al., Mon. Not. Roy. Astron. Soc. 412, L108 (2011), 1012.0045.
  15. C. O. Heinke and W. C. G. Ho, Astrophys. J. 719, L167 (2010), 1007.4719.
  16. J. M. Lattimer, M. Prakash, C. J. Pethick, and P. Haensel, Phys. Rev. Lett. 66, 2701 (1991).
  17. T. Schafer and K. Schwenzer, Phys. Rev. D70, 114037 (2004), astro-ph/0410395.
  18. M. Alford, M. Braby, M. W. Paris, and S. Reddy, Astrophys. J. 629, 969 (2005), nucl-th/0411016.
  19. A. Akmal, V. R. Pandharipande, and D. G. Ravenhall, Phys. Rev. C58, 1804 (1998), nucl-th/9804027.
  20. P. Haensel and R. Schaeffer, Phys. Rev. D45, 4708 (1992).
  21. M. Prakash, M. Prakash, J. M. Lattimer, and C. J. Pethick, Astrophys. J. Lett. 390, L77 (1992).
  22. N. Iwamoto, Phys. Rev. Lett. 44, 1637 (1980).
  23. B. A. Sa’d, I. A. Shovkovy, and D. H. Rischke, Phys. Rev. D75, 125004 (2007), astro-ph/0703016.
  24. D. Page, J. M. Lattimer, M. Prakash, and A. W. Steiner, Astrophys. J. Suppl. 155, 623 (2004), astro-ph/0403657.
  25. A. Schmitt, I. A. Shovkovy, and Q. Wang, Phys. Rev. D73, 034012 (2006), hep-ph/0510347.
  26. X. Wang and I. A. Shovkovy, Phys.Rev. D82, 085007 (2010), 1006.1293.
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 minumum 40 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