2024ApJ...974..244T

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 ${{\rm{c}}}_{{\rm{s}}}^{2}\leqslant 0.781$. In this work, we demonstrate that this upper bound is corroborated by the reconstructed equation of state (EOS; modeled with a nonparametric method) for ultradense matter. The reconstruction integrates multimessenger observation for neutron stars, in particular, the latest radius measurements for PSR J0437-4715 (${11.36}_{-0.63}^{+0.95}$ km), PSR J0030+0451 (${11.71}_{-0.83}^{+0.88}$ km, in the ST+PDT model), and PSR J0740+6620 (${12.49}_{-0.88}^{+1.28}$ km) by NICER have been adopted. The result shows in all cases, the ${{\rm{c}}}_{{\rm{s}}}^{2}\leqslant 0.781$ upper limit for EOS will naturally yield the properties of matter near the center of the massive neutron star consistent with the causality-driven constraint from pQCD, where, in practice, the density in implementing the pQCD likelihood (n_L) is applied at $10{n}_{s}$ (where ${n}_{s}$ is the nuclear saturation density). We also note that there is a strong correlation for the maximum c_s² with n_L, and ${{\rm{c}}}_{{\rm{s}}}^{2}\leqslant 0.781$ is somehow violated when n_L = n_c,TOV. The result indicates that a higher n_L, 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. In addition, the latest joint constraint for R_1.4, R_2.0, R_1.4 - R_2.0, and M_TOV are ${11.94}_{-0.68}^{+0.77}$ km, ${11.99}_{-0.67}^{+0.88}$ km, $-{0.1}_{-0.27}^{+0.42}$ km, and ${2.24}_{-0.10}^{+0.13}{M}_{\odot }$ (at 90% credible level), respectively.