Ab initio methods using weakly interacting nucleons give a good description of condensed nuclear matter up to densities comparable to the nuclear saturation density. At higher densities palpable strong interactions between overlapping nucleons become important; we propose that the interactions will continuously switch over to follow a holographic model in this region. In order to implement this, we construct hybrid equations of state (EoSs) where various models are used for low density nuclear matter, and the holographic V-QCD model is used for non-perturbative high density nuclear matter as well as for quark matter. We carefully examine all existing constraints from astrophysics of compact stars and discuss their implications for the hybrid EoSs. Thanks to the stiffness of the V-QCD EoS for nuclear matter, we obtain a large family of viable hybrid EoSs passing the constraints. We find that quark matter cores in neutron stars are unstable due to the strongly first order deconfinement transition, and predict bounds on the tidal deformability as well as on the radius of neutron stars. By relying on universal relations, we also constrain characteristic peak frequencies of gravitational waves produced in neutron star mergers.

We use the holographic V-QCD models to analyse the physics of dense QCD and neutron stars. Accommodating lattice results for thermodynamics of QCD enables us to make generic predictions for the Equation of State (EoS) of the quark matter phase in the cold and dense regime. We demonstrate that the resulting pressure in V-QCD matches well with a family of neutron-star-matter EoSs that interpolate between state-of-the-art theoretical results for low and high density QCD. After implementing the astrophysical constraints, i.e., the largest known neutron star mass and the recent LIGO/Virgo results for the tidal deformability, we analyse the phase transition between the baryonic and quark matter phases. We find that the baryon density $n_B$ at the transition is at least 2.9 times the nuclear saturation density $n_s$. The transition is of strongly first order at low and intermediate densities, i.e., for $n_B/n_s \lesssim 7.5$.