Upper Limit of Sound Speed in Nuclear Matter: A Harmonious Interplay of Transport Calculation and Perturbative QCD Constraint

Very recently, it has been shown that there is an upper bound on the squared sound speed of nuclear matter from the transport, which reads cₛ² ≤ 0.781 {2024arXiv240214085H}. In this work, we demonstrate that this upper bound is corroborated by the reconstructed equation of state (EOS) for ultra-dense matter. The reconstruction integrates multi-messenger observation for neutron stars (NSs), as well as the theoretical constraints, including that from chiral effective field theory and perturbative quantum chromodynamics (pQCD). In particular, the latest radius measurements for PSR J0437-4715 (11.36+0.95-0.63 km) and PSR J0030+0451 (11.71+0.88-0.83 km, in the ST+PDT model) by NICER have been adopted. The result shows in all cases, the cₛ² ≤ 0.781 upper limit for EOS will naturally yield the properties of matter near center of the massive neutron star consistent with the causality-driven constraint from pQCD, where in practice, the constraint is applied at ten nuclear saturation density (nL=10nₛ). We also note that there is a strong correlation for the maximum square of sound speed cₛ² with nL, and cₛ² ≤ 0.781 is somehow violated when nL = nc,TOV. The result indicates that a higher density in implementing the pQCD constraint, even considering the uncertainties from statistics, is more natural. Moreover, the remarkable agreement between the outcomes derived from these two distinct and independent constraints (i.e., the transport calculation and pQCD boundary) lends strong support to their validity.