2018ApJ...859...90Z

Within the parameter space of the equation of state (EOS) of dense neutron-rich matter limited by existing constraints mainly from terrestrial nuclear experiments, we investigate how the neutron star maximum mass M_max > 2.01 ± 0.04 M_☉, radius 10.62 km < R_1.4 < 12.83 km and tidal deformability Λ_1.4 <= 800 of canonical neutron stars together constrain the EOS of dense neutron-rich nucleonic matter. While the 3D parameter space of K_sym (curvature of nuclear symmetry energy), J_sym, and J₀ (skewness of the symmetry energy and EOS of symmetric nuclear matter, respectively) is narrowed down significantly by the observational constraints, more data are needed to pin down the individual values of K_sym, J_sym, and J₀. The J₀ largely controls the maximum mass of neutron stars. While the EOS with J₀ = 0 is sufficiently stiff to support neutron stars as massive as 2.37 M_☉, supporting the hypothetical ones as massive as 2.74 M_☉ (composite mass of GW170817) requires J₀ to be larger than its currently known maximum value of about 400 MeV and beyond the causality limit. The upper limit on the tidal deformability of Λ_1.4 = 800 from the recent observation of GW170817 is found to provide upper limits on some EOS parameters consistent with but far less restrictive than the existing constraints of other observables studied.