2023ApJ...949...11J

We have investigated the systematic differences introduced when performing a Bayesian-inference analysis of the equation of state (EOS) of neutron stars employing either variable- or constant-likelihood functions. The former has the advantage of retaining the full information on the distributions of the measurements, making exhaustive usage of the data. The latter, on the other hand, has the advantage of a much simpler implementation and reduced computational costs. In both approaches, the EOSs have identical priors and have been built using the sound speed parameterization method so as to satisfy the constraints from X-ray and gravitational waves observations, as well as those from chiral effective theory and perturbative quantum chromodynamics. In all cases, the two approaches lead to very similar results and the 90% confidence levels essentially overlap. Some differences do appear, but in regions where the probability density is extremely small and are mostly due to the sharp cutoff on the binary tidal deformability $\tilde{{\rm{\Lambda }}}\leqslant 720$ set in the constant-likelihood approach. Our analysis has also produced two additional results. First, an inverse correlation between the normalized central number density, n_c,TOV/n_s, and the radius of a maximally massive star, R_TOV. Second, and most importantly, it has confirmed the relation between the chirp mass and the binary tidal deformability. The importance of this result is that it relates ${{ \mathcal M }}_{\mathrm{chirp}}$, which is measured very accurately, and $\tilde{{\rm{\Lambda }}}$, which contains important information on the EOS. Hence, when ${{ \mathcal M }}_{\mathrm{chirp}}$ is measured in future detections, our relation can be used to set tight constraints on $\tilde{{\rm{\Lambda }}}$.