Combined Constraints on the Equation of State of Dense Neutron-rich Matter from Terrestrial Nuclear Experiments and Observations of Neutron Stars

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-circle dot, radius 10.62 km < R-1.4< 12.83 km and tidal deformability Lambda(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(0 )(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(0). The J(0) largely controls the maximum mass of neutron stars. While the EOS with J(0) = 0 is sufficiently stiff to support neutron stars as massive as 2.37 M-circle dot, supporting the hypothetical ones as massive as 2.74 M-0 (composite mass of GW170817) requires J(0) 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 Lambda(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.