In this paper, we discuss the damping of density oscillations in dense nuclear matter in the temperature range relevant to neutron star mergers. This damping is due to bulk viscosity arising from the weak interaction ``Urca'' processes of neutron decay and electron capture. The nuclear matter is modelled in the relativistic density functional approach. The bulk viscosity reaches a resonant maximum close to the neutrino trapping temperature, then drops rapidly as temperature rises into the range where neutrinos are trapped in neutron stars. We investigate the bulk viscous dissipation timescales in a post-merger object and identify regimes where these timescales are as short as the characteristic timescale $\sim$10 ms, and, therefore, might affect the evolution of the post-merger object. Our analysis indicates that bulk viscous damping would be important at not too high temperatures of the order of a few MeV and densities up to a few times saturation density.

Supernovae and cooling neutron stars have long been used to constrain the properties of axions, such as their mass and interactions with nucleons and other Standard Model particles. We investigate the prospects of using neutron star mergers as a similar location where axions can be probed in the future. We examine the impact axions would have on mergers, considering both the possibility that they free-stream through the dense nuclear matter and the case where they are trapped. We calculate the mean free path of axions in merger conditions, and find that they would free-stream through the merger in all thermodynamic conditions. In contrast to previous calculations, we integrate over the entire phase space while using a relativistic treatment of the nucleons, assuming the matrix element is momentum-independent. In particular, we use a relativistic mean field theory to describe the nucleons, taking into account the precipitous decrease in the effective mass of the nucleons as density increases above nuclear saturation density. We find that within current constraints on the axion-neutron coupling, axions could cool nuclear matter on timescales relevant to neutron star mergers. Our results may be regarded as first steps aimed at understanding how axions affect merger simulations and potentially interface with observations.

We compute the mass, radius and tidal deformability of stars containing phase transitions from hadronic to quark phase(s). These quantities are computed for three types of hadronic envelopes: purely nuclear, hyperonic, and $\Delta$-resonance--hyperon admixed matter. We consider either a single first-order phase transition to a quark phase with a maximally stiff equation of state (EOS) or two sequential first-order phase transitions mimicking a transition from hadronic to a quark matter phase followed by a second phase transition to another quark phase. We explore the parameter space which produces low-mass twin and triplet configurations where equal-mass stars have substantially different radii and tidal deformabilities. We demonstrate that while for purely hadronic stiff EOS the obtained maximum mass is inconsistent with the upper limit on this quantity placed by GW170817, the inclusion of the hyperonic and $\Delta$-resonance degrees of freedom, as well as the deconfinement phase transition at sufficiently low density, produces a configuration of stars consistent with this limit. The obtained hybrid star configurations are in the mass range relevant for the interpretation of the GW170817 event. We compare our results for the tidal deformability with the limits inferred from GW170817 showing that the onset of non-nucleonic phases, such as $\Delta$-resonance--hyperon admixed phase and/or the quark phase(s), are favored by this data if the nuclear EOS is stiff. Also, we show that low-mass twins and especially triplets proliferate the number of combinations of possible types of stars that can undergo a merger event, the maximal number being six in the case of triplets. The prospects for uncovering the first-order phase transition(s) to and in quark matter via measurements of tidal deformabilities in merger events are discussed.

We study bulk viscosity arising from weak current Urca processes in dense baryonic matter at and beyond nuclear saturation density. We consider the temperature regime where neutrinos are trapped and therefore have nonzero chemical potential. We model the nuclear matter in a relativistic density functional approach, taking into account the trapped neutrino component. We find that the resonant maximum of the bulk viscosity would occur at or below the neutrino trapping temperature, so in the neutrino trapped regime the bulk viscosity decreases with temperature as $T^{-2}$, this decrease being interrupted by a drop to zero at a special temperature where the proton fraction becomes density-independent and the material scale-invariant. The bulk viscosity is larger for matter with lower lepton fraction, i.e., larger isospin asymmetry. We find that bulk viscosity in the neutrino-trapped regime is smaller by several orders of magnitude than in the neutrino-transparent regime, which implies that bulk viscosity in neutrino-trapped matter is probably not strong enough to affect significantly the evolution of neutron star mergers. This also implies weak damping of gravitational waves emitted by the oscillations of the postmerger remnant in the high-temperature, neutrino-trapped phase of evolution.

We calculate the bulk-viscous dissipation time for adiabatic density oscillations in nuclear matter at densities of 1-7 times nuclear saturation density and at temperatures ranging from 1 MeV, where corrections to previous low-temperature calculations become important, up to 10 MeV, where the assumption of neutrino transparency is no longer valid. Under these conditions, which are expected to occur in neutron star mergers, damping of density oscillations arises from beta equilibration via weak interactions. We find that for 1 kHz oscillations the shortest dissipation times are in the 5 to 20 ms range, depending on the equation of state, which means that bulk viscous damping could affect the dynamics of a neutron star merger. For higher frequencies the dissipation time can be even shorter.

We review the prospects for detecting quark matter in neutron star cores. We survey the proposed signatures and emphasize the importance of data from neutron star mergers, which provide access to dynamical properties that operate on short timescales that are not probed by other neutron star observables.

We contribute a publicly available set of tables and code to provide Equations of State (EoS) for matter at neutron star densities. Our EoSes are constrained only by input from hadron physics and fundamental principles, without feedback from neutron star observations, and so without relying on General Relativity. They can therefore be used to test General Relativity itself, as well as modified gravity theories, with neutron star observables, without logical circularity. We have adapted state of the art results from NN chiral potentials for the low--density limit, pQCD results for the asymptotically high-density EoS, and use monotonicity and causality as the only restrictions for intermediate densities, for the EoS sets to remain as model-independent as is feasible today.

We show that the commonly used criterion for beta equilibrium in neutrino-transparent dense nuclear matter becomes invalid as temperatures rise above 1 MeV. Such temperatures are attained in neutron star mergers. By numerically computing the relevant weak interaction rates we find that the correct criterion for beta equilibrium requires an isospin chemical potential that can be as large as 10-20 MeV, depending on the temperature at which neutrinos become trapped.

Compact stars may contain quark matter in their interiors at densities exceeding several times the nuclear saturation density. We explore models of such compact stars where there are two first-order phase transitions: the first from nuclear matter to a quark-matter phase, followed at higher density by another first-order transition to a different quark matter phase [e.g., from the two-flavor color superconducting (2SC) to the color-flavor-locked (CFL) phase). We show that this can give rise to two separate branches of hybrid stars, separated from each other and from the nuclear branch by instability regions and, therefore, to a new family of compact stars, denser than the ordinary hybrid stars. In a range of parameters, one may obtain twin hybrid stars (hybrid stars with the same masses but different radii) and even triplets where three stars, with inner cores of nuclear matter, 2SC matter, and CFL matter, respectively, all have the same mass but different radii.

Inferring the properties of dense matter is one of the most exciting prospects from the measurement of gravitational waves from neutron star mergers. However, it will require reliable numerical simulations that incorporate viscous dissipation and energy transport if these can play a significant role within the survival time of the post-merger object. We calculate timescales for typical forms of dissipation and find that thermal transport and shear viscosity will not be important unless neutrino trapping occurs, which requires temperatures above about 10 MeV and gradients over lengthscales of 0.1 km or less. On the other hand, if direct-Urca processes remain suppressed, leaving modified-Urca processes to establish flavor equilibrium, then bulk viscous dissipation could provide significant damping to density oscillations observed right after the merger. When comparing with data from a state-of-the-art merger simulation we find that the bulk viscosity takes values close to its resonant maximum in a typical neutron-star merger, motivating a more careful assessment of the role of bulk viscous dissipation in the gravitational-wave signal from merging neutron stars.

We investigate the stability of stars with a density discontinuity between a high-density core and a very low density mantle. Previous work on "strange dwarfs" suggested that such a discontinuity could stabilize stars that would have been classified as unstable by the conventional criteria based on extrema in the mass-radius relation. We investigate the stability of such stars by numerically solving the Sturm-Liouville equations for the lowest-energy modes of the star. We find that the conventional criteria are correct, and strange dwarfs are not stable.

In nuclear matter at neutron-star densities and temperatures, Cooper pairing leads to the formation of a gap in the nucleon excitation spectra resulting in exponentially strong Boltzmann suppression of many transport coefficients. Previous calculations have shown evidence that density oscillations of sufficiently large amplitude can overcome this suppression for flavor-changing $\beta$ processes, via the mechanism of "gap bridging". We address the simplifications made in that initial work, and show that gap bridging can counteract Boltzmann suppression of neutrino emissivity for the realistic case of modified Urca processes in matter with $^3P_2$ neutron pairing.

At very high densities, as for example in the core of a neutron star, matter may appear in the color-flavor locked (CFL) phase, which is a superfluid. This phase features topologically stable vortex solutions, which arise in a spinning superfluid as localized configurations carrying quanta of angular momentum. Despite the topological stability of these vortices they are not the lowest energy state of the system at neutron star densities and decay into triplets of semi-superfluid fluxtubes. In these proceedings we report on the progress of our numerical study in the Ginzburg-Landau approach, where we investigate lattices of semi-superfluid fluxtubes. The fluxtubes are obtained through controlled decay of global vortex configurations in the presence of a gauge field. Understanding the dynamics of semi-superfluid string configurations is important in the context of angular momentum transfer from a quark matter core of a neutron star beyond the core boundary, since not vortex-, but fluxtube pinning seems to be the relevant mechanism in this scenario.

The core region of a neutron star may feature quark matter in the color-flavor- locked (CFL) phase. The CFL condensate breaks the baryon number symmetry, such that the phenomenon of superfluidity arises. If the core of the star is rotating, vortices will form in the superfluid, carrying the quanta of angular momentum. In a previous study we have solved the question of stability of these vortices, where we found numerical proof of a conjectured instability, according to which superfluid vortices will decay into an arrangement of so-called semi-superfluid fluxtubes. Here we report first results of an extension of our framework that allows us to study multi-vortex dynamics. This will in turn enable us to investigate the structure of semi-superfluid string lattices, which could be relevant to study pinning phenomena at the boundary of the core.

The discovery of pulsars as heavy as 2 solar masses has led astrophysicists to rethink the core compositions of neutron stars, ruling out many models for the nuclear equations of state (EoS). We explore the hybrid stars that occur when hadronic matter is treated in a relativistic mean-field approximation and quark matter is modeled by three-flavor local and non-local Nambu Jona-Lasinio (NJL) models with repulsive vector interactions. The NJL models typically yield equations of state that feature a first order transition to quark matter. Assuming that the quark-hadron surface tension is high enough to disfavour mixed phases, and restricting to EoSes that allow stars to reach 2 solar masses, we find that the appearance of the quark matter core either destabilizes the star immediately (this is typical for non-local NJL models) or leads to a very short hybrid star branch in the mass-radius relation (this is typical for local NJL models). Using the Constant-Sound-Speed parametrization we can see that the reason for the near-absence of hybrid stars is that the transition pressure is fairly high and the transition is strongly first order.

We describe two aspects of the physics of hybrid stars that have a sharp interface between a core of quark matter and a mantle of nuclear matter. Firstly, we analyze the mass-radius relation. We describe a generic "Constant Speed of Sound" (CSS) parameterization of the quark matter equation of state (EoS), in which the speed of sound is independent of density. In terms of the three parameters of the CSS EoS we obtain the phase diagram of possible forms of the hybrid star mass-radius relation, and we show how observational constraints on the maximum mass and typical radius of neutron stars can be expressed as constraints on the CSS parameters. Secondly, we propose a mechanism for the damping of density oscillations, including r-modes, in hybrid stars with a sharp interface. The dissipation arises from the periodic conversion between quark matter and nuclear matter induced by the pressure oscillations in the star. We find the damping grows nonlinearly with the amplitude of the oscillation and is powerful enough to saturate an r-mode at very low saturation amplitude, of order $10^{-10}$, which is compatible with currently-available observations of neutron star spin frequencies and temperatures.

Superfluid vortices in the color-flavor-locked (CFL) phase of dense quark matter are known to be energetically disfavored relative to well-separated triplets of "semi-superfluid" color flux tubes. However, the short-range interaction (metastable versus unstable) has not been established. In this paper we perform numerical calculations using the effective theory of the condensate field, mapping the regions in the parameter space of coupling constants where the vortices are metastable versus unstable. For the case of zero gauge coupling we analytically identify a candidate for the unstable mode, and show that it agrees well with the results of the numerical calculations. We find that in the region of the parameter space that seems likely to correspond to real-world CFL quark matter the vortices are unstable, indicating that if such matter exists in neutron star cores it is very likely to contain semi-superfluid color flux tubes rather than superfluid vortices.

We propose a mechanism for the damping of density oscillations in multicomponent compact stars. The mechanism is the periodic conversion between different phases, i.e., the movement of the interface between them, induced by pressure oscillations in the star. The damping grows nonlinearly with the amplitude of the oscillation. We study in detail the case of r-modes in a hybrid star with a sharp interface, and we find that this mechanism is powerful enough to saturate the r-mode at very low saturation amplitude, of order $10^{-10}$, and is therefore likely to be the dominant r-mode saturation mechanism in hybrid stars with a sharp interface.

We discuss the "constant speed of sound" (CSS) parametrization of the equation of state of high-density matter and its application to the field correlator method (FCM) model of quark matter. We show how observational constraints on the maximum mass and typical radius of neutron stars are expressed as constraints on the CSS parameters. We find that the observation of a $2\,M_{\odot}$ star already severely constrains the CSS parameters, and is particularly difficult to accommodate if the squared speed of sound in the high-density phase is assumed to be around $1/3$ or less. We show that the FCM equation of state can be accurately represented by the CSS parametrization, which assumes a sharp transition to a high-density phase with density-independent speed of sound. We display the mapping between the FCM and CSS parameters, and see that FCM only allows equations of state in a restricted subspace of the CSS parameters.

The cores of compact stars reach the highest densities in nature and therefore could consist of novel phases of matter. We demonstrate via a detailed analysis of pulsar evolution that precise pulsar timing data can constrain the star's composition, through unstable global oscillations (r-modes) whose damping is determined by microscopic properties of the interior. If not efficiently damped, these modes emit gravitational waves that quickly spin down a millisecond pulsar. As a first application of this general method, we find that ungapped interacting quark matter is consistent with both the observed radio and x-ray data, whereas for ordinary nuclear matter some additional enhanced damping mechanism is required.

Background: The two-flavor color superconducting (2SC) phase of quark matter is a possible constituent of the core of neutron stars. To assess its impact on the observable behavior of the star one must analyze transport properties, which in 2SC matter are controlled by the scattering of gapless fermionic modes by each other and possibly also by color-magnetic flux tubes. Purpose: We determine the electrical and thermal conductivities and the shear viscosity of 2SC matter. Methods: We use a variational formulation of transport theory, treating the strong and electromagnetic interactions via a weak coupling expansion. Results: We provide the leading order scaling of the transport coefficients with temperature and chemical potential as well as accurate fits to our numerical results. We also find that the scattering of fermions by color-magnetic flux tubes is insignificant for thermal conductivity, but may contribute to the electrical conductivity and shear viscosity in the limit of very low temperature or high magnetic field. We also estimate the transport coefficients in unpaired quark matter. Conclusions: Our calculation has set the stage for exploration of possible signatures of the presence of 2SC quark matter in neutron stars.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

Astrophysical observations of compact stars provide, in addition to collider experiments, the other big source of information on matter under extreme conditions. The largest and most precise data set about neutron stars is the timing data of radio pulsars. We show how this unique data can be used to learn about the ultra-dense matter in the compact star interior. The method relies on astro-seismology based on special global oscillation modes (r-modes) that emit gravitational waves. They would prevent pulsars from spinning with their observed high frequencies, unless the damping of these modes, determined by the microscopic properties of matter, can prevent this. We show that for each form of matter there is a distinct region in a frequency/spindown-rate diagram where r-modes can be present. We find that stars containing ungapped quark matter are consistent with both the observed radio and x-ray data, whereas, even when taking into account the considerable uncertainties, neutron star models with standard viscous damping are inconsistent with both data sets and additional damping mechanisms would be required.

We highlight the progress, current status, and open challenges of QCD-driven physics, in theory and in experiment. We discuss how the strong interaction is intimately connected to a broad sweep of physical problems, in settings ranging from astrophysics and cosmology to strongly-coupled, complex systems in particle and condensed-matter physics, as well as to searches for physics beyond the Standard Model. We also discuss how success in describing the strong interaction impacts other fields, and, in turn, how such subjects can impact studies of the strong interaction. In the course of the work we offer a perspective on the many research streams which flow into and out of QCD, as well as a vision for future developments.

Neutron stars undergoing r-mode oscillation emit gravitational radiation that might be detected on earth. For known millisecond pulsars the observed spindown rate imposes an upper limit on the possible gravitational wave signal of these sources. Taking into account the physics of r-mode evolution, we show that only sources spinning at frequencies above a few hundred Hertz can be unstable to r-modes, and we derive a more stringent universal r-mode spindown limit on their gravitational wave signal, exploiting the fact that the r-mode saturation amplitude is insensitive to the structural properties of individual sources. We find that this refined bound limits the gravitational wave strain from millisecond pulsars to values below the detection sensitivity of next-generation detectors. Young sources are therefore a more promising option for the detection of gravitational waves emitted by r-modes and to probe the interior composition of compact stars in the near future.

We study the mass-radius curve of hybrid stars, assuming a single first-order phase transition between nuclear and quark matter, with a sharp interface between the quark matter core and nuclear matter mantle. We use a generic parameterization of the quark matter equation of state, which has a constant, i.e. density-independent, speed of sound ("CSS"). We argue that this parameterization provides a framework for comparison and empirical testing of models of quark matter. We obtain the phase diagram of possible forms of the hybrid star mass-radius relation, where the control parameters are the transition pressure, energy density discontinuity, and the quark matter speed of sound. We find that this diagram is sensitive to the quark matter parameters but fairly insensitive to details of the nuclear matter equation of state. We calculate the maximum hybrid star mass as a function of the parameters of the quark matter EoS, and find that there are reasonable values of those parameters that give rise to hybrid stars with mass above $2\,M_\odot$.

The rotation frequencies of young pulsars are systematically below their theoretical Kepler limit. R-modes have been suggested as a possible explanation for this observation. With the help of semi-analytic expressions that make it possible to assess the uncertainties of the r-mode scenario due to the impact of uncertainties in underlying microphysics, we perform a quantitative analysis of the spin-down and the emitted gravitational waves of young pulsars. We find that the frequency to which r-modes spin down a young neutron star is surprisingly insensitive both to the microscopic details and the saturation amplitude. Comparing our result to astrophysical data, we show that for a range of sufficiently large saturation amplitudes r-modes provide a viable spindown scenario and that all observed young pulsars are very likely already outside the r-mode instability region. Therefore the most promising sources for gravitational wave detection are unobserved neutron stars associated with recent supernovae, and we find that advanced LIGO should be able to see several of them. We find the remarkable result that the gravitational wave strain amplitude is completely independent of both the r-mode saturation amplitude and the microphysics, and depends on the saturation mechanism only within some tens of per cent. However, the gravitational wave frequency depends on the amplitude and we provide the required expected timing parameter ranges to look for promising sources in future searches.

We study the spindown of pulsars due to gravitational wave emission and show that r-modes in neutron stars provide a quantitative explanation for the observed low rotation frequencies of young pulsars if the r-mode saturation amplitude is sufficiently large. In contrast for such large saturation amplitudes old hadronic millisecond pulsars would spin down much faster than observed. We discuss resolutions of this apparent discrepancy that could make this mechanism consistent with the observational data.

Studying the frequency and temperature evolution of a compact star can give us valuable information about the microscopic properties of the matter inside the star. In this paper we study the effect of dissipative reheating of a neutron star due to r-mode oscillations on its temperature evolution. We find that there is still an impact of an r-mode phase on the temperature long after the star has left the instability region and the r-mode is damped completely. With accurate temperature measurements it may be possible to detect this trace of a previous r-mode phase in observed pulsars.

If the surface tension of quark matter is low enough, quark matter is not self bound. At sufficiently low pressure and temperature, it will take the form of a crystal of positively charged strangelets in a neutralizing background of electrons. In this case there will exist, in addition to the usual family of strange stars, a family of low-mass large-radius objects analogous to white dwarfs, which we call "strangelet dwarfs". Using a generic parametrization of the equation of state of quark matter, we calculate the mass-radius relationship of these objects.

We analyze the viscous damping of r-mode oscillations of compact stars, taking into account non-linear viscous effects in the large-amplitude regime. The qualitatively different cases of hadronic stars, strange quark stars, and hybrid stars are studied. We calculate the viscous damping times of r-modes, obtaining numerical results and also general approximate analytic expressions that explicitly exhibit the dependence on the parameters that are relevant for a future spindown evolution calculation. The strongly enhanced damping of large amplitude oscillations leads to damping times that are considerably lower than those obtained when the amplitude dependence of the viscosity is neglected. Consequently, large-amplitude viscous damping competes with the gravitational instability at all physical frequencies and could stop the r-mode growth in case this is not done before by non-linear hydrodynamic mechanisms.

We study the viscous damping of r-modes of compact stars and analyze in detail the regions where small amplitude modes are unstable to the emission of gravitational radiation. We present general expressions for the viscous damping times for arbitrary forms of interacting dense matter and derive general semi-analytic results for the boundary of the instability region. These results show that many aspects, like in particular the physically important minima of the instability boundary, are surprisingly insensitive to detailed microscopic properties of the considered form of matter. Our general expressions are applied to the cases of hadronic stars, strange stars, and hybrid stars, and we focus on equations of state that are compatible with the recent measurement of a heavy compact star. We find that hybrid stars with a sufficiently small core can "masquerade" as neutron stars and feature an instability region that is indistinguishable from that of a neutron star, whereas neutron stars with a core density high enough to allow direct Urca reactions feature a notch on the right side of the instability region.

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.

Pulsar spin frequencies and their time evolution are an important source of information on compact stars and their internal composition. Oscillations of the star can reduce the rotational energy via the emission of gravitational waves. In particular unstable oscillation modes, like r-modes, are relevant since their amplitude becomes large and can lead to a fast spin-down of young stars if they are saturated by a non-linear saturation mechanism. We present a novel mechanism based on the pronounced large-amplitude enhancement of the bulk viscosity of dense matter. We show that the enhanced damping due to non-linear bulk viscosity can saturate r-modes of neutron stars at amplitudes appropriate for an efficient spin-down.

The recent measurement of the Shapiro delay in the radio pulsar PSR J1614-2230 yielded a mass of 1.97 +/- 0.04 M_sun, making it the most massive pulsar known to date. Its mass is high enough that, even without an accompanying measurement of the stellar radius, it has a strong impact on our understanding of nuclear matter, gamma-ray bursts, and the generation of gravitational waves from coalescing neutron stars. This single high mass value indicates that a transition to quark matter in neutron-star cores can occur at densities comparable to the nuclear saturation density only if the quarks are strongly interacting and are color superconducting. We further show that a high maximum neutron-star mass is required if short duration gamma-ray bursts are powered by coalescing neutron stars and, therefore, this mechanism becomes viable in the light of the recent measurement. Finally, we argue that the low-frequency (<= 500 Hz) gravitational waves emitted during the final stages of neutron-star coalescence encode the properties of the equation of state because neutron stars consistent with this measurement cannot be centrally condensed. This will facilitate the measurement of the neutron star equation of state with Advanced LIGO/Virgo.

We study the bulk viscosity of dense matter, taking into account non-linear effects which arise in the large amplitude "supra-thermal" region where the deviation $\mu_\Delta$ of the chemical potentials from chemical equilibrium fulfills $\mu_\Delta>T$. This regime is relevant to unstable modes such as r-modes, which grow in amplitude until saturated by non-linear effects. We study the damping due to direct and modified Urca processes in hadronic matter, and due to nonleptonic weak interactions in strange quark matter. We give general results valid for an arbitrary equation of state of dense matter and find that the viscosity can be strongly enhanced by supra-thermal effects. Our study confirms previous results on quark matter and shows that the non-linear enhancement is even stronger in the case of hadronic matter. Our results can be applied to calculations of the r-mode-induced spin-down of fast-rotating neutron stars, where the spin-down time will depend on the saturation amplitude of the r-mode

We argue that if color-superconducting quark matter exists in the core of a neutron star, it may contain a high density of flux tubes, carrying flux that is mostly color-magnetic, with a small admixture of ordinary magnetic flux. We focus on the two-flavor color-superconducting ("2SC") phase, and assume that the flux tubes are energetically stable, although this has not yet been demonstrated. The density of flux tubes depends on the nature of the transition to the color-superconducting phase, and could be within an order of magnitude of the density of magnetic flux tubes that would be found if the core were superconducting nuclear matter. We calculate the cross-section for Aharonov-Bohm scattering of gapless fermions off the flux tubes, and the associated collision time and frictional force on a moving flux tube. We discuss the other forces on the flux tube, and find that if we take in to account only the forces that arise within the 2SC core region then the timescale for expulsion of the color flux tubes from the 2SC core is of order 10^10 years.

Color-flavor locked (CFL) quark matter at high densities is a color superconductor, which spontaneously breaks baryon number and chiral symmetry. Its low-energy thermodynamic and transport properties are therefore dominated by the H (superfluid) boson, and the octet of pseudoscalar pseudo-Goldstone bosons of which the neutral kaon is the lightest. We study the CFL-K^0 phase, in which the stress induced by the strange quark mass causes the kaons to condense, and there is an additional ultra-light "K^0" Goldstone boson arising from the spontaneous breaking of isospin. We compute the bulk viscosity of matter in the CFL-K^0 phase, which arises from the beta-equilibration processes K^0<->H+H and K^0+H<->H. We find that the bulk viscosity varies as T^7, unlike the CFL phase where it is exponentially Boltzmann-suppressed by the kaon's energy gap. However, in the temperature range of relevance for r-mode damping in compact stars, the bulk viscosity in the CFL-K^0 phase turns out to be even smaller than in the uncondensed CFL phase, which already has a bulk viscosity much smaller than all other known color-superconducting quark phases.

It has recently been pointed out that if the surface tension of quark matter is low enough, the surface of a strange star will be a crust consisting of a crystal of charged strangelets in a neutralizing background of electrons. This affects the behavior of the surface, and must be taken into account in efforts to observationally rule out strange stars. We calculate the thickness of this ``mixed phase'' crust, taking into account the effects of surface tension and Debye screening of electric charge. Our calculation uses a generic parametrization of the equation of state of quark matter. For a reasonable range of quark matter equations of state, and surface tension of order a few MeV/fm^2, we find that the preferred crystal structure always involves spherical strangelets, not rods or slabs of quark matter. We find that for a star of radius 10 km and mass 1.5 Msolar, the strangelet-crystal crust can be from zero to hundreds of meters thick, the thickness being greater when the strange quark is heavier, and the surface tension is smaller. For smaller quark stars the crust will be even thicker.

Matter at high density and low temperature is expected to be a color superconductor, which is a degenerate Fermi gas of quarks with a condensate of Cooper pairs near the Fermi surface that induces color Meissner effects. At the highest densities, where the QCD coupling is weak, rigorous calculations are possible, and the ground state is a particularly symmetric state, the color-flavor locked (CFL) phase. The CFL phase is a superfluid, an electromagnetic insulator, and breaks chiral symmetry. The effective theory of the low-energy excitations in the CFL phase is known and can be used, even at more moderate densities, to describe its physical properties. At lower densities the CFL phase may be disfavored by stresses that seek to separate the Fermi surfaces of the different flavors, and comparison with the competing alternative phases, which may break translation and/or rotation invariance, is done using phenomenological models. We review the calculations that underlie these results, and then discuss transport properties of several color-superconducting phases and their consequences for signatures of color superconductivity in neutron stars.

We study the behavior of Goldstone bosons in color-flavor-locked (CFL) quark matter at nonzero temperature. Chiral symmetry breaking in this phase of cold and dense matter gives rise to pseudo-Goldstone bosons, the lightest of these being the charged and neutral kaons K^+ and K^0. At zero temperature, Bose-Einstein condensation of the kaons occurs. Since all fermions are gapped, this kaon condensed CFL phase can, for energies below the fermionic energy gap, be described by an effective theory for the bosonic modes. We use this effective theory to investigate the melting of the condensate: we determine the temperature-dependent kaon masses self-consistently using the two-particle irreducible effective action, and we compute the transition temperature for Bose-Einstein condensation. Our results are important for studies of transport properties of the kaon condensed CFL phase, such as bulk viscosity.

The bulk viscosity of three-flavor color-superconducting quark matter originating from the nonleptonic process u+s <-> u+d is computed. It is assumed that up and down quarks form Cooper pairs while the strange quark remains unpaired (2SC phase). A general derivation of the rate of strangeness production is presented, involving contributions from a multitude of different subprocesses, including subprocesses that involve different numbers of gapped quarks as well as creation and annihilation of particles in the condensate. The rate is then used to compute the bulk viscosity as a function of the temperature, for an external oscillation frequency typical of a compact star r-mode. We find that, for temperatures far below the critical temperature T_c for 2SC pairing, the bulk viscosity of color-superconducting quark matter is suppressed relative to that of unpaired quark matter, but for T >~ 10^(-3) T_c the color-superconducting quark matter has a higher bulk viscosity. This is potentially relevant for the suppression of r-mode instabilities early in the life of a compact star.

The bulk viscosities of two color-superconducting phases, the color-flavor locked (CFL) phase and the 2SC phase, are computed and compared to the result for unpaired quark matter. In the case of the CFL phase, processes involving kaons and the superfluid mode give the largest contribution to the bulk viscosity since all fermionic modes are gapped. In the case of the 2SC phase, ungapped fermionic modes are present and the process u+d <-> u+s provides the dominant contribution. In both cases, the bulk viscosity can become larger than that of the unpaired phase for sufficiently large temperatures (T >~ 1 MeV for CFL, T >~ 0.1 MeV for 2SC). Bulk viscosity (as well as shear viscosity) is important for the damping of r-modes in compact stars and thus can potentially be used as an indirect signal for the presence or absence of color-superconducting quark matter.

We calculate the bulk viscosity of color-superconducting quark matter in the color-flavor-locked (CFL) phase. We assume that the lightest bosons are the superfluid mode H and the kaons K0 and K+, and that there is no kaon condensate. We calculate the rate of strangeness-equilibrating processes that convert kaons into superfluid modes, and the resultant bulk viscosity. We find that for oscillations with a timescale of milliseconds, at temperatures T << 1 MeV, the CFL bulk viscosity is much less than that of unpaired quark matter, but at higher temperatures the bulk viscosity of CFL matter can become larger.

Ozel, in a recent reanalysis of EXO 0748-676 observational data (astro-ph/0605106), concluded that quark matter probably does not exist in the center of compact stars. We show that the data is actually consistent with the presence of quark matter in compact stars.

At ultra-high density, matter is expected to form a degenerate Fermi gas of quarks in which there is a condensate of Cooper pairs of quarks near the Fermi surface: color superconductivity. In this chapter we review some of the underlying physics, and discuss outstanding questions about the phase structure of ultra-dense quark matter. We then focus on describing recent results on the crystalline color superconducting phase that may be the preferred form of cold, dense but not asymptotically dense, three-flavor quark matter. The gap parameter and free energy for this phase have recently been evaluated within a Ginzburg-Landau approximation for many candidate crystal structures. We describe the two that are most favorable. The robustness of these phases results in their being favored over wide ranges of density. However, it also implies that the Ginzburg-Landau approximation is not quantitatively reliable. We describe qualitative insights into what makes a crystal structure favorable which can be used to winnow the possibilities. We close with a look ahead at the calculations that remain to be done in order to make quantitative contact with observations of compact stars.

We construct strangelets, taking into account electrostatic effects, including Debye screening, and arbitrary surface tension sigma of the interface between vacuum and quark matter. We find that there is a critical surface tension sigma_crit below which large strangelets are unstable to fragmentation and below which quark star surfaces will fragment into a crystalline crust made of charged strangelets immersed in an electron gas. We derive a model-independent relationship between sigma_crit and two parameters that characterize any quark matter equation of state. For reasonable model equations of state, we find sigma_crit typically of order a few MeV/fm^2. If sigma <= sigma_crit, the size-distribution of strangelets in cosmic rays could feature a peak corresponding to the stable strangelets that we construct.

We show that a hybrid (nuclear + quark matter) star can have a mass-radius relationship very similar to that predicted for a star made of purely nucleonic matter. We show this for a generic parameterization of the quark matter equation of state, and also for an MIT bag model, each including a phenomenological correction based on gluonic corrections to the equation of state. We obtain hybrid stars as heavy as 2 M_solar for reasonable values of the bag model parameters. For nuclear matter, we use the equation of state calculated by Akmal, Pandharipande, and Ravenhall using many-body techniques. Both mixed and homogeneous phases of nuclear and quark matter are considered.

The gapless color-flavor locked (gCFL) phase is the second-densest phase of matter in the QCD phase diagram, making it a plausible constituent of the core of neutron stars. We show that even a relatively small region of gCFL matter in a star will dominate both the heat capacity C_V and the heat loss by neutrino emission L_\nu. The gCFL phase is characterized by an unusual quasiparticle dispersion relation that makes both its specific heat c_V and its neutrino emissivity epsilon_\nu parametrically larger than in any other phase of nuclear or quark matter. During the epoch in which the cooling of the star is dominated by direct Urca neutrino emission, the presence of a gCFL region does not strongly alter the cooling history because the enhancements of C_V and L_\nu cancel against each other. At late times, however, the cooling is dominated by photon emission from the surface, so L_\nu is irrelevant, and the anomalously large heat capacity of the gCFL region keeps the star warm. The temperature drops with time as T\sim t^{-1.4} rather than the canonical T\sim t^{-5}. This provides a unique and potentially observable signature of gCFL quark matter.

According to quantum chromodynamics (QCD), matter at ultra-high densities will take the form of a color-superconducting quark liquid, in which there is a condensate of Cooper pairs of quarks near the Fermi surface. I present a review of the physics of color superconductivity. I give particular attention to the recently proposed gapless CFL (gCFL) phase, which has unusual properties such as quasiquarks with a near-quadratic dispersion relation, and which may well be the favored phase of quark matter in the density range relevant to compact stars. I also discuss the effects of color superconductivity on the mass-radius relationship of compact stars, showing that one would have to fix the bag constant by other measurements in order to see the effects of color superconductivity. An additional parameter in the quark matter equation of state connected with perturbative corrections allows quark matter to imitate nuclear matter over the relevant density range, so that hybrid stars can show a mass-radius relationship very similar to that of nuclear matter, and their masses can reach 1.9 M_solar.

The densest predicted state of matter is colour-superconducting quark matter, in which quarks near the Fermi surface form a condensate of Cooper pairs. This form of matter may well exist in the core of compact stars, and the search for signatures of its presence is an ongoing enterprise. Using a bag model of quark matter, I discuss the effects of colour superconductivity on the mass-radius relationship of compact stars, showing that colour superconducting quark matter can occur in compact stars at values of the bag constant where ordinary quark matter would not be allowed. The resultant ``hybrid'' stars with colour superconducting quark matter interior and nuclear matter surface have masses in the range 1.3-1.6 Msolar and radii 8-11 km. Once perturbative corrections are included, quark matter can show a mass-radius relationship very similar to that of nuclear matter, and the mass of a hybrid star can reach $1.8 \Msolar$.

We study compact stars that contain quark matter. We look at the effect of color superconductivity in the quark matter on the nuclear-quark matter transition density, mass-radius relationship, and the density discontinuity at the boundary between nuclear and quark matter. We find that color superconducting quark matter will occur in compact stars at values of the bag constant where ordinary quark matter would not be allowed. We are able to construct ``hybrid'' stars with a color superconducting quark matter interior and nuclear matter surface, with masses in the range 1.3-1.6 M_solar and radii 8-11 km. Our results are consistent with recent mass-radius limits based on absorption lines from EXO0748-676.

At nuclear matter density, electrically neutral strongly interacting matter in weak equilibrium is made of neutrons, protons and electrons. At sufficiently high density, such matter is made of up, down and strange quarks in the color-flavor locked phase, with no electrons. As a function of increasing density (or, perhaps, increasing depth in a compact star) other phases may intervene between these two phases which are guaranteed to be present. The simplest possibility, however, is a single first order phase transition between CFL and nuclear matter. Such a transition, in space, could take place either through a mixed phase region or at a single sharp interface with electron-free CFL and electron-rich nuclear matter in stable contact. Here we construct a model for such an interface. It is characterized by a region of separated charge, similar to an inversion layer at a metal-insulator boundary. On the CFL side, the charged boundary layer is dominated by a condensate of negative kaons. We then consider the energetics of the mixed phase alternative. We find that the mixed phase will occur only if the nuclear-CFL surface tension is significantly smaller than dimensional analysis would indicate.

In any context in which color superconductivity arises in nature, it is likely to involve pairing between species of quarks with differing chemical potentials. For suitable values of the differences between chemical potentials, Cooper pairs with nonzero total momentum are favored, as was first realized by Larkin, Ovchinnikov, Fulde and Ferrell (LOFF). Condensates of this sort spontaneously break translational and rotational invariance, leading to gaps which vary periodically in a crystalline pattern. Unlike the original LOFF state, these crystalline quark matter condensates include both spin zero and spin one Cooper pairs. We explore the range of parameters for which crystalline color superconductivity arises in the QCD phase diagram. If in some shell within the quark matter core of a neutron star (or within a strange quark star) the quark number densities are such that crystalline color superconductivity arises, rotational vortices may be pinned in this shell, making it a locus for glitch phenomena.

After a brief review of the phenomena expected in cold dense quark matter, color superconductivity and color-flavor locking, we sketch some implications of recent developments in our understanding of cold dense quark matter for the physics of compact stars. We give a more detailed summary of our recent work on crystalline color superconductivity and the consequent realization that (some) pulsar glitches may originate in quark matter.

We discuss the Meissner effect for a color superconductor formed by cold dense quark matter. Though color and ordinary electromagnetism are broken in a color superconductor, there is a linear combination of the photon and a gluon that remains massless. Consequently, a color superconducting region may be penetrated by an external magnetic field. We show that at most a small fraction of the magnetic field is expelled, and if the screening distance is the smallest length scale in the problem there is no expulsion at all. We calculate the behavior of the magnetic field for a spherical geometry relevant for compact stars. If a neutron star contains a quark matter core, this core is a color superconductor. Our results demonstrate that such cores admit magnetic fields without restricting them to quantized flux tubes. Such magnetic fields within color superconducting neutron star cores are stable on time scales longer than the age of the universe, even if the spin period of the neutron star is changing.