“…11, for the neutron star model predicted with the empirical formula of ∆R Pb n , we only plot the stellar model whose central density is up to threefold the saturation density. In the same figure, we show the constraints on the neutron star mass and radius obtained from the various astronomical observations; the gravitational wave observations in the GW170817 event, i.e., the 1.4M neutron star radius is less than 13.6 km [8], whose constraint may become more severe by combining with the multimessenger observations and nuclear theory [71,72]; the x-ray observations via NICER for PSR J0030+0451 [16,17] and MSP J0740+6620 [18,19]; the observations of x-ray burst through the theoretical models [73]; and the identification of the magnetar quasi-periodic oscillations observed in GRB 200415A with the crustal torsional oscillations [74]. As the theoretical constraint, the top-left region can be excluded from the causality [75].…”