2019ApJ...887...48B

This paper presents a new analysis of the thermal emission from the neutron star (NS) surface to constrain the dense matter equation of state. We employ an empirical parameterization of the equation of state with a Markov Chain Monte Carlo approach to consistently fit the spectra of quiescent low-mass X-ray binaries in globular clusters with well-measured distances. Despite previous analyses predicting low NS radii, we show that it is possible to reconcile the astrophysical data with nuclear physics knowledge with or without including a prior on the slope of the symmetry energy L_sym. With this empirical parameterization of the equation of state, we obtain radii of the order of about 12 km without worsening the fit statistic. More importantly, we obtain the following values for the slope of the symmetry energy, its curvature K_sym, and the isoscalar skewness parameter Q_sat: L_sym=37.2_–8.9^+9.2 MeV, K_sym=-85_–70⁺⁸² MeV, and Q_sat=318_–366⁺⁶⁷³ MeV. These are the first measurements of the empirical parameters K_sym and Q_sat. Their values are only weakly impacted by our assumptions, such as the distances or the number of free empirical parameters, provided the latter are taken within a reasonable range. We also study the weak sensitivity of our results to the set of sources analyzed, and we identify a group of sources that dominates the constraints. The resulting masses and radii obtained from this empirical parameterization are also compared to other measurements from electromagnetic observations of NSs and gravitational wave signals from the NS-NS merger GW170817.