Three parametrizations DOPS1, DOPS2, and DOPS3 (named after the Department of Physics Shimla) of the relativistic mean field model have been proposed with the inclusion of all possible self and mixed interactions between the scalar-isoscalar ($\ensuremath{\sigma}$), vector-isoscalar ($\ensuremath{\omega}$), and vector-isovector ($\ensuremath{\rho}$) mesons up to quartic order. The generated parameter sets are in harmony with the finite and bulk nuclear matter properties. A set of equations of state (EOSs) composed of pure hadronic (nucleonic) matter and nucleonic with quark matter (hybrid EOSs) for superdense hadron-quark matter in $\ensuremath{\beta}$ equilibrium is obtained. The quark matter phase is calculated by using the three-flavor Nambu-Jona-Lasinio (NJL) model. The maximum mass of a nonrotating neutron star with DOPS1 parametrization is found to be around $2.6M\ensuremath{\bigodot}$ for the pure nucleonic matter, which satisfies the recent gravitational wave analysis of GW190814 [Abbott et al., Astrophys. J. Lett. 896, L44 (2020)] with possible maximum mass constraint indicating that the secondary component of GW190814 could be a nonrotating heaviest neutron star composed of pure nucleonic matter. EOSs computed with the DOPS2 and DOPS3 parametrizations satisfy the x-ray observational data [Steiner et al., Astrophys. J. 722, 33 (2010)] and the recent observations of GW170817 maximum mass constraint of a stable nonrotating neutron star in the range $2.01\ifmmode\pm\else\textpm\fi{}0.04--2.16\ifmmode\pm\else\textpm\fi{}0.03M\ensuremath{\bigodot}$ [Rezzolla et al., Astrophys. J. Lett. 852, L25 (2018)] and also in good agreement with constraints on mass and radius measurement for PSR $\mathrm{J}0740+6620$ (NICER) [Riley et al., Astrophys. J. Lett. 918, L27 (2021); Miller et al., Astrophys. J. Lett. 918, L28 (2021)]. The hybrid EOSs obtained with the NJL model also satisfy astrophysical constraints on the maximum mass of a neutron star from PSR J1614-2230 [Demorest et al., Nature (London) 467, 1081 (2010)]. We also present the results for dimensionless tidal deformability, $\mathrm{\ensuremath{\Lambda}}$ which are consistent with the waveform models analysis of GW170817.
Read full abstract