Towards mitigation of apparent tension between nuclear physics and astrophysical observations by improved modeling of neutron star matter

Observations of neutron stars (NSs) by the LIGO-Virgo and NICER collaborations have provided reasonably precise measurements of their various macroscopic properties. In this paper, we employ a Bayesian framework to combine them and place improved joint constraints on the properties of NS equation of state (EoS). We use a hybrid EoS formulation that employs a parabolic expansion-based nuclear empirical parameterization around the nuclear saturation density augmented by a generic 3-segment piecewise polytrope model at higher densities. Within the $90 \%$ credible level this parameterization predicts $R_{1.4} = 12.57_{-0.92}^{+0.73}$ km and $\Lambda_{1.4} = 550_{-225}^{+223}$ for the radius and dimensionless tidal deformability, respectively, of a $1.4 M_{\odot}$ NS. Finally, we show how the construction of the full NS EoS based solely on the nuclear empirical parameters at saturation density leads to certain tension with the astrophysical data, and how the hybrid approach provides a resolution to it.