What does perturbative QCD really have to say about neutron stars

The implications of perturbative QCD (pQCD) calculations on neutron stars are carefully examined. While pQCD calculations above baryon chemical potentials $\mu_B\sim2.4$ GeV demonstrate the potential of ruling out a wide range of neutron star equations of state (EOSs), such constraints only affect the most massive neutron stars in the vicinity of the Tolman-Oppenheimer-Volkoff (TOV) limit, resulting in constraints that are orthogonal to current or expected astrophysical bounds. In the most constraining scenario, pQCD considerations favor low values of the squared speed sound $C_s$ at high $\mu_B$ relevant for the most massive neutron stars, but leave predictions of the radii and tidal deformabilities almost unchanged. Such considerations become irrelevant if the maximum speed of sound squared inside neutron stars does not exceed about $C_{s,\mathrm{max}}\sim0.5$, or if pQCD breaks down below $\mu_B\simeq2.9$ GeV. Furthermore, the large pQCD uncertainties preclude any meaningful bounds on the neutron star EOS at the moment. Interestingly, if pQCD predictions for the pressure at around $\mu_B\simeq2.5$ GeV are refined and found to be low ($\lesssim 1.5$ GeV/fm$^3$), evidence for a soft neutron star inner core EOS in combination with the existence of two-solar-mass pulsars would indicate the presence of color superconductivity beyond neutron star densities. We point out that two-solar-mass pulsars place robust upper bounds on this non-perturbative effect and require the pairing gap to be less than $\Delta_\mathrm{CFL}\lesssim500$ MeV at $\mu_B\simeq2.5$ GeV.