r-mode oscillations have been shown to have a significant potential to constrain the composition of fast spinning neutron stars. Due to their high rotation rates, millisecond pulsars (MSPs) provide a unique platform to constrain the properties of such oscillations, if their surface temperatures can be inferred. We present the results of our investigations of archival X-ray data of a number of MSPs, including recent XMM-Newton observations of PSR J1810+1744 and PSR J2241-5236. Using the neutron star atmosphere (NSA) model and taking into account various uncertainties, we present new bounds on the surface temperature of these sources. We were then able to significantly strengthen previous bounds on the amplitude of the r-mode oscillations in millisecond pulsars and find values as low as $\alpha \lesssim 10^{-9}$. This is by now about three orders of magnitude below what standard saturation mechanisms in neutron stars could provide, which requires very strong dissipation in the interior, strongly pointing towards a structurally complex or exotic composition of these sources. At such low temperatures sources could even be outside of the instability region, and taking into account the various uncertainties we obtain for an observed surface temperature a simple frequency bound below which r-modes are excluded in slower spinning pulsars.

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.

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.

r-mode astroseismology provides a unique way to study the internal composition of compact stars. Due to their precise timing, recycled millisecond radio pulsars present a particularly promising class of sources. Although their thermal properties are still poorly constrained, X-ray data is very useful for astroseismology since r-modes could strongly heat a star. Using known and new upper bounds on the temperatures and luminosities of several non-accreting millisecond radio pulsars we derive bounds on the r-mode amplitude as low as $\alpha\lesssim10^{-8}$ and discuss the impact on scenarios for their internal composition.

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.

Next generation gravitational wave detectors will start taking data in the near future. Here we discuss the chances to detect the continuous emission from r-mode oscillations in compact stars and study which properties of compact stars we can infer from such novel data. In particular we show that the combination of the gravitational wave data with electromagnetic multi-messenger observations could give us detailed insight into compact star properties, ranging from precise mass-radius measurements to the determination of the equation of state and the phase structure of dense matter.

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.

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.

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.

Long-range interactions lead to non-Fermi liquid effects in dense matter. We show that, in contrast to other material properties, their effect on the bulk viscosity of quark matter is significant since they shift its resonant maximum and can thereby change the viscosity by many orders of magnitude. This is of importance for the damping of oscillations of compact stars, like in particular unstable r-modes, and the quest to detect signatures of deconfined matter in astrophysical observations. We find that, in contrast to neutron stars with standard damping mechanisms, compact stars that contain ungapped quark matter are consistent with the observed data on low mass x-ray binaries.

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.

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.

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 study neutrino emission from a normal, ungapped, quark phase in the core of a compact star. Neutrino emission from non-interacting quark matter leads to an emissivity that scales as $\epsilon\sim T^7$. We show that the emissivity is enhanced by a combination of Fermi liquid and non-Fermi liquid effects. Fermi liquid effects lead to an emissivity that scales as $\epsilon\sim \alpha_s T^6$, as originally shown by Iwamoto. We demonstrate that non-Fermi liquid effects further enhance the rate, leading to $\epsilon\sim \alpha_s^3 T^6\log(m/T)^2$, where $m$ is the electric screening scale and $m\gg T$ under the conditions found in compact stars. We show, however, that combined with non-Fermi liquid effects in the specific heat the enhancement in the emissivity only leads to a modest reduction in the temperature of the star at late times. Our results confirm existing bounds on the presence of ungapped quark matter in compact stars. We also discuss neutrino emission from superconducting phases with ungapped fermionic excitations.