We study dynamical friction in interacting relativistic systems with arbitrary mean free paths and medium constituent masses. Our novel framework recovers the known limits of ideal gas and ideal fluid when the mean free path goes to infinity or zero, respectively, and allows for a smooth interpolation between these limits. We find that in an infinite system the drag force can be expressed as a sum of ideal-gas-like and ideal-fluid-like contributions leading to a finite friction even at subsonic velocities. This simple picture receives corrections in any finite system and the corrections become especially significant for a projectile moving at a velocity $v$ close to the speed of sound $v\approx c_s$. These corrections smoothen the ideal fluid discontinuity around the speed of sound and render the drag force a continuous function of velocity. We show that these corrections can be computed to a good approximation within effective theory of viscous fluid dynamics.

In this conference-proceedings contribution, we review recent advances in placing model-independent constraints on the properties of cold and dense QCD matter inside neutron stars. In addition to introducing new bounds for the Equation of State, we explain how these results may be used to make robust statements about the physical phase of strongly interacting matter in the centers of neutron stars of different masses. Our findings indicate that the existence of quark-matter cores inside massive neutron stars appears to be a very common feature of the allowed Equations of State, and should not be considered an exotic or unlikely scenario.

We combine results from astrophysical observations and theoretical ab-initio calculations to investigate the presence of quark matter (QM) inside neutron stars (NSs). We find a clear qualitative change in the material properties of NS matter at energy densities comparable to those where quark-gluon plasma is created in heavy-ion collisions. In the low-density phase, the system has characteristics closely resembling those of hadronic matter, while the high-density phase can be clearly identified with nearly conformal QM. We show that QM is never present inside 1.44-M_sol NSs, but resides inside maximally massive NSs barring very specific and extreme conditions. Finally, for the heaviest observed NSs of approximately two solar masses, the fate of QM is found to strongly depend on the behavior of the speed of sound in the matter: if this quantity does not strongly violate the conformal limit c_s^2 < 1/3, these stars host sizable quark cores of at least five kilometers in radius.

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

At high baryon chemical potential mu_B, the equation of state of QCD allows a weak-coupling expansion in the QCD coupling alpha_s. The result is currently known up to and including the full next-to-next-to-leading order (NNLO) alpha_s^2. Starting at this order, the computations are complicated by the modification of particle propagation in a dense medium, which necessitates non-perturbative treatment of the scale Sqrt(alpha_s)*mu_B. In this work, we apply a Hard-Thermal-Loop scheme for capturing the contributions of this scale to the weak-coupling expansion, and use it to determine the leading-logarithm contribution to NNNLO: alpha_s^3 (ln alpha_s)^2. This result is the first improvement to the equation of state of massless cold quark matter in 40 years. The new term is negligibly small, and thus significantly increases our confidence in the applicability of the weak-coupling expansion.

The LIGO/Virgo detection of gravitational waves originating from a neutron-star merger, GW170817, has recently provided new stringent limits on the tidal deformabilities of the stars involved in the collision. Combining this measurement with the existence of two-solar-mass stars, we generate a generic family of neutron-star-matter Equations of State (EoSs) that interpolate between state-of-the-art theoretical results at low and high baryon density. Comparing the results to ones obtained without the tidal-deformability constraint, we witness a dramatic reduction in the family of allowed EoSs. Based on our analysis, we conclude that the maximal radius of a 1.4-solar-mass neutron star is 13.6 km, and that smallest allowed tidal deformability of a similar-mass star is $\Lambda(1.4 M_\odot) = 120$.

We generalize the state-of-the-art perturbative Equation of State of cold quark matter to nonzero temperatures, needed in the description of neutron star mergers and core collapse processes. The new result is accurate to order g^5 in the gauge coupling, and is based on a novel framework for dealing with the infrared sensitive soft field modes of the theory. The zero Matsubara mode sector is treated using a dimensionally reduced effective theory, while the soft non-zero modes are resummed using the Hard Thermal Loop approximation. This combination of known effective descriptions offers unprecedented access to small but nonzero temperatures, both in and out of beta equilibrium.

In this review article, we argue that our current understanding of the thermodynamic properties of cold QCD matter, originating from first principles calculations at high and low densities, can be used to efficiently constrain the macroscopic properties of neutron stars. In particular, we demonstrate that combining state-of-the-art results from Chiral Effective Theory and perturbative QCD with the current bounds on neutron star masses, the Equation of State of neutron star matter can be obtained to an accuracy better than 30% at all densities.

In recent years, there have been several successful attempts to constrain the equation of state of neutron star matter using input from low-energy nuclear physics and observational data. We demonstrate that significant further restrictions can be placed by additionally requiring the pressure to approach that of deconfined quark matter at high densities. Remarkably, the new constraints turn out to be highly insensitive to the amount --- or even presence --- of quark matter inside the stars.

Lattice QCD studies of the thermodynamics of the hot quark-gluon plasma (QGP) demonstrate the importance of accounting for the interactions of quarks and gluons, if one wants to investigate the phase structure of strongly interacting matter. Motivated by this observation and using state-of-the-art results from perturbative QCD, we construct a simple effective equation of state for cold quark matter that consistently incorporates the effects of interactions and furthermore includes a built-in estimate of the inherent systematic uncertainties. This goes beyond the MIT bag model description in a crucial way, yet leads to an equation of state that is equally straightforward to use.

While QCD appears not to be accurately solvable in the regime of interest for neutron star physics, microscopic calculations are feasible at both low and very high densities. In this work, we propose using the most realistic calculations in these two regimes of nuclear physics and perturbative QCD, and construct equations of state by matching the results requiring thermodynamic consistency. We find that the resulting equations of state --- in contrast to several hadronic ones --- are able to reproduce current observational data on neutron stars without any fine tuning, and allow stable hybrid stars with masses up to 2.1M_{sun}. Using recent observations of star radii, we perform a maximum likelihood analysis to further constrain the equation of state, and in addition show that the effects of rotation on radii and masses should not be neglected in future precision studies.

We perform an O(alpha_s^2) perturbative calculation of the equation of state of cold but dense QCD matter with two massless and one massive quark flavor, finding that perturbation theory converges reasonably well for quark chemical potentials above 1 GeV. Using a running coupling constant and strange quark mass, and allowing for further non-perturbative effects, our results point to a narrow range where absolutely stable strange quark matter may exist. Absent stable strange quark matter, our findings suggest that quark matter in compact star cores becomes confined to hadrons only slightly above the density of atomic nuclei. Finally, we show that equations of state including quark matter lead to hybrid star masses up to M~2M_solar, in agreement with current observations. For strange stars, we find maximal masses of M~2.75M_solar and conclude that confirmed observations of compact stars with M>2M_solar would strongly favor the existence of stable strange quark matter.