
Neutron star mass formula with nuclear saturation parameters for asymmetric nuclear matter
Hajime Sotani1, 2, ∗and Shinsuke Ota3
1Astrophysical Big Bang Laboratory, RIKEN, Saitama 351-0198, Japan
2Interdisciplinary Theoretical & Mathematical Science Program (iTHEMS), RIKEN, Saitama 351-0198, Japan
3Research Center for Nuclear Physics (RCNP), Osaka University, Ibaragi, Osaka 567-0047, Japan
(Dated: October 24, 2022)
Low-mass neutron stars are directly associated with the nuclear saturation parameters because their central
density is definitely low. We have already found a suitable combination of nuclear saturation parameters for
expressing the neutron star mass and gravitational redshift, i.e., η≡(K0L2)1/3with the incompressibility for
symmetric nuclear matter, K0, and the density-dependent nuclear symmetry energy, L. In this study, we newly
find another suitable combination given by ητ≡(−KτL5)1/6with the isospin dependence of incompressibility
for asymmetric nuclear matter, Kτ, and derive the empirical relations for the neutron star mass and gravitational
redshift as a function of ητand the normalized central number density. With these empirical relations, one
can evaluate the mass and gravitational redshift of the neutron star, whose central number density is less than
threefold the saturation density, within ∼10% accuracy, and the radius within a few % accuracy. In addition,
we discuss the neutron star mass and radius constraints from the terrestrial experiments, using the empirical
relations, together with those from the astronomical observations. Furthermore, we find a tight correlation
between ητand η. With this correlation, we derive the constraint on Kτas −348 ≤Kτ≤ −237 MeV,
assuming that L= 60 ±20 and K0= 240 ±20 MeV.
PACS numbers: 04.40.Dg, 26.60.+c, 21.65.Ef
I. INTRODUCTION
Neutron star is a massive remnant left after a supernova explosion, which happens at the last moment of the massive star’s life.
The density inside the star becomes much higher than the standard nuclear density, ρ0= 2.7×1014 g/cm3, and the gravitational
and magnetic fields inside/around the star are significantly stronger than those observed in our solar system [1]. The neutron
star mass and radius strongly depend on the equation of state (EOS) for neutron star matter under the β-equilibrium. The mass
of a neutron star model with a higher central density generally becomes larger, even though the EOS is not fixed yet. So, the
astronomical observations of neutron stars or their activities tell us the information about the EOS for a relatively higher density
region, while the terrestrial experiments tell us that for a lower density region (e.g., Fig. 2 in Ref. [2]).
In practice, the discovery of 2Mneutron stars [3–6] has excluded some soft EOSs, with which the expected maximum mass
is less than the observed mass. This discovery simultaneously reveals the problem that most of the EOSs with hyperon are too
soft to construct the 2Mneutron star, i.e., the so-called hyperon puzzle. Meanwhile, the gravitational wave event, GW170817
[7], tells us the information on the tidal deformability of the neutron star, which leads to the constraint that a 1.4Mneutron
star radius should be less than 13.6 km [8]. We note that the constraint on neutron star radius may become more stringent in
view of the existing multi-messenger observational data [9, 10]. The light bending due to the strong gravitational field induced
by the neutron star is also one of the important phenomena to see the neutron star properties. That is, the pulsar light curve
from the rotating neutron star would be modified due to this relativistic effect and one could get the neutron star properties
by carefully observing it (e.g., [11–16]). Based on this idea, the Neutron star Interior Composition Explorer (NICER) is now
operating on an International Space Station (ISS) and it has already announced the constraint on two neutron stars properties,
i.e., PSR J0030+0451 [17, 18] and PSR J0740+6620 [19, 20]. Furthermore, the direct detection of the gravitational waves from
the neutron star in the future may enable us to extract the neutron star properties (e.g., [21–29]).
On the other hand, the EOS in a lower density region is also gradually constrained through terrestrial nuclear experiments,
but still, there are large uncertainties in EOS parameters (or in neutron star properties) constrained from terrestrial experiments.
For instance, the fiducial value of the density-dependent nuclear symmetry energy Lis L'60 ±20 MeV [30, 31], while the
constraints of Lobtained recently seem to be significantly larger than the fiducial value [32, 33]. This is because one has to
usually transform the experimental constraint to the EOS parameters, even if the information determined via experiments is
associated with some aspects of nuclear EOS. Then, one can eventually discuss the neutron star mass and radius as a solution of
the Tolman-Oppenheimer-Volkoff (TOV) equation. Anyway, the terrestrial experiments are definitely crucial for understanding
the neutron star EOS as well as the astronomical observation of neutron stars.
∗sotani@yukawa.kyoto-u.ac.jp
arXiv:2210.11651v1 [nucl-th] 21 Oct 2022