Astrophysical constraints on nuclear EOSs and coupling constants in RMF models

Abstract: Utilizing various astrophysical constraints on neutron star structures, we carry out a Bayesian analysis on the density-dependent behaviors of coupling constants in RMF models as well as the nuclear matter properties at supranuclear densities. The effective nucleon interactions in the isoscalar-scalar, isoscalar-vector, and isovector-vector channels are considered, where the corresponding coupling constants ($\alpha_S, \alpha_V, \alpha_{TV}$) are fixed by dividing entire density range into three regions with six independent parameters. In this work we focus on constraining the density-dependent point-coupling constants at supranuclear densities, while the coupling constants at subsaturation densities are derived from the covariant density functional DD-ME2. For those consistent with astrophysical observations, the coupling constants generally decrease with density and approach to small positive values at large enough densities, which qualitatively agrees with various RMF models. The posterior probability density functions and their correlations of the coupling constants and various nuclear matter properties are examined as well. At $1\sigma$ level, the constrained coupling constants at density $1.5n_0$ ($2.5n_0$) are $\alpha_S = 3.1^{+0.1}_{-0.05} (1.55^{+0.85}_{-0.2}) \times 10^{-4} \mathrm{MeV}^{-2}$, $\alpha_V = 2.3^{+0.1}_{-0.0} (1.3^{+0.55}_{-0.1}) \times 10^{-4} \mathrm{MeV}^{-2}$, and $\alpha_{TV} = 2.05^{+0}_{-0.4} (2.05^{+0}_{-0.5})\times 10^{-5} \mathrm{MeV}^{-2}$. At larger densities, we find the lower limit of $\alpha_{TV}$ is not well constrained, so that more extensive calculations with larger number of free parameters are necessary.