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.

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