Deconfined quark matter may be present in the cores of massive neutron stars, and is in addition expected to be produced in binary neutron star mergers. While the bulk thermodynamic properties of the medium may undergo only moderate changes at the transition density, transport properties are expected to be qualitatively different in the nuclear and quark matter phases, possibly leading to observable effects. Motivated by recent advances in neutron star observations, we perform the first-ever systematic study of the bulk and shear viscosities as well as the thermal and electrical conductivities in dense unpaired quark matter using holography. Working in the bottom-up V-QCD model and comparing our results to the predictions of a top-down D3-D7 system, we arrive at results that are in qualitative disagreement with the available leading-order predictions from perturbative QCD. Our findings highlight the differing transport properties of weakly and strongly interacting systems, and call for caution in the use of the perturbative results in neutron-star applications.

The detection of gravitational waves from a neutron star merger has opened up the possibility of detecting the presence or creation of deconfined quark matter using the gravitational wave signal. To investigate this possibility, we construct a family of neutron star matter equations of state at nonzero density and temperature by combining state-of-the-art nuclear matter equations of state with holographic equations of state for strongly interacting quark matter. The emerging picture consistently points toward a strong first order deconfinement transition, with a temperature-dependent critical density and latent heat that we quantitatively examine. Recent neutron star mass measurements are further used to discriminate between the different equations of state obtained, leaving a tightly constrained family of preferred equations of state.

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.

We investigate a simple holographic model for cold and dense deconfined QCD matter consisting of three quark flavors. Varying the single free parameter of the model and utilizing a Chiral Effective Theory equation of state (EoS) for nuclear matter, we find four different compact star solutions: traditional neutron stars, strange quark stars, as well as two non-standard solutions we refer to as hybrid stars of the second and third kind (HS2 and HS3). The HS2s are composed of a nuclear matter core and a crust made of stable strange quark matter, while the HS3s have both a quark mantle and a nuclear crust on top of a nuclear matter core. For all types of stars constructed, we determine not only their mass-radius relations, but also tidal deformabilities, Love numbers, as well as moments of inertia and the mass distribution. We find that there exists a range of parameter values in our model, for which the novel hybrid stars have properties in very good agreement with all existing bounds on the stationary properties of compact stars. In particular, the tidal deformabilities of these solutions are smaller than those of ordinary neutron stars of the same mass, implying that they provide an excellent fit to the recent gravitational wave data GW170817 of LIGO and Virgo. The assumptions underlying the viability of the different star types, in particular those corresponding to absolutely stable quark matter, are finally discussed at some length.

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$.

According to common lore, Equations of State of field theories with gravity duals tend to be soft, with speeds of sound either below or around the conformal value of $v_s=1/\sqrt{3}$. This has important consequences in particular for the physics of compact stars, where the detection of two solar mass neutron stars has been shown to require very stiff equations of state. In this paper, we show that no speed limit exists for holographic models at finite density, explicitly constructing examples where the speed of sound becomes arbitrarily close to that of light. This opens up the possibility of building hybrid stars that contain quark matter obeying a holographic equation of state in their cores.

It has been conjectured that the speed of sound in holographic models with UV fixed points has an upper bound set by the value of the quantity in conformal field theory. If true, this would set stringent constraints for the presence of strongly coupled quark matter in the cores of physical neutron stars, as the existence of two-solar-mass stars appears to demand a very stiff Equation of State. In this article, we present a family of counterexamples to the speed of sound conjecture, consisting of strongly coupled theories at finite density. The theories we consider include ${\cal N}=4$ super Yang-Mills at finite R-charge density and non-zero gaugino masses, while the holographic duals are Einstein-Maxwell theories with a minimally coupled scalar in a charged black hole geometry. We show that for a small breaking of conformal invariance, the speed of sound approaches the conformal value from above at large chemical potentials.

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.

We use a top-down holographic model for strongly interacting quark matter to study the properties of neutron stars. When the corresponding Equation of State (EoS) is matched with state-of-the-art results for dense nuclear matter, we consistently observe a first order phase transition at densities between two and seven times the nuclear saturation density. Solving the Tolman-Oppenheimer-Volkov equations with the resulting hybrid EoSs, we find maximal stellar masses in the excess of two solar masses, albeit somewhat smaller than those obtained with simple extrapolations of the nuclear matter EoSs. Our calculation predicts that no quark matter exists inside neutron stars.

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.

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.

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.

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.

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.