Unified neutron star EOSs and neutron star structures in RMF models

TL;DR

This study systematically compares various covariant density functionals in modeling neutron star matter, analyzing their EOSs, structures, and observational consistency, highlighting the importance of accurate measurements for constraining nuclear physics models.

Contribution

It provides a comprehensive comparison of multiple covariant density functionals in neutron star EOSs and structures, emphasizing the impact of symmetry energy parameters and observational constraints.

Findings

EOSs are similar at low and certain mid densities but vary elsewhere.

Larger symmetry energy slope L increases neutron drip density and related parameters.

Most functionals predict neutron stars exceeding 2 solar masses, with some matching observational radii.

Abstract

In the framework of Thomas-Fermi approximation, we study systematically the EOSs and microscopic structures of neutron star matter in a vast density range with $n_\mathrm{b}\approx 10^{-10}$-2 $\mathrm{fm}^{-3}$, where various covariant density functionals are adopted, i.e., those with nonlinear self couplings (NL3, PK1, TM1, GM1, MTVTC) and density-dependent couplings (DD-LZ1, DDME-X, PKDD, DD-ME2, DD2, TW99). It is found that the EOSs generally coincide with each other at $n_\mathrm{b}\lesssim 10^{-4}$ fm${}^{-3}$ and 0.1 fm${}^{-3}\lesssim n_\mathrm{b} \lesssim 0.3$ fm${}^{-3}$, while in other density regions they are sensitive to the effective interactions between nucleons. By adopting functionals with larger slope of symmetry energy $L$, the curvature parameter $K_\mathrm{sym}$ and neutron drip density generally increase, while the droplet size, proton number of nucleus, core-crust transition density, and onset density of non-spherical nuclei decrease. All functionals predict neutron stars with maximum masses exceeding the two-solar-mass limit, while those of DD2, DD-LZ1, DD-ME2, and DDME-X predict optimum neutron star radii according to the observational constraints. Nevertheless, the corresponding skewness coefficients $J$ are much lager than expected, while only the functionals MTVTC and TW99 meet the start-of-art constraints on $J$. More accurate measurements on the radius of PSR J0740+6620 and the maximum mass of neutron stars are thus essential to identify the functional that satisfies all constraints from nuclear physics and astrophysical observations. Approximate linear correlations between neutron stars' radii at $M=1.4 M_{\odot}$ and $2 M_{\odot}$, the slope $L$ and curvature parameter $K_\mathrm{sym}$ of symmetry energy are observed as well, which is mainly attributed to the curvature-slope correlations in the functionals adopted here.