2021ApJ...908L..28N

The detection of the binary events GW170817 and GW190814 has provided invaluable constraints on the maximum mass of nonrotating configurations of neutron stars, M_TOV. However, the large differences in the neutron-star masses measured in GW170817 and GW190814 has also lead to significant tension between the predictions for such maximum masses, with GW170817 suggesting that M_TOV <= 2.3 M_☉, and GW190814 requiring M_TOV >= 2.5 M_☉ if the secondary was a (non- or slowly rotating) neutron star at merger. Using a genetic algorithm, we sample the multidimensional space of parameters spanned by gravitational-wave and astronomical observations associated with GW170817. Consistent with previous estimates, we find that all of the physical quantities are in agreement with the observations if the maximum mass is in the range of M_TOV=2.210_–0.123^+0.116M_☉ within a 2σ confidence level. By contrast, maximum masses with M_TOV >= 2.5 M_☉, not only require efficiencies in the gravitational-wave emission that are well above the numerical-relativity estimates, but they also lead to a significant underproduction of the ejected mass. Hence, the tension can be released by assuming that the secondary in GW190814 was a black hole at merger, although it could have been a rotating neutron star before.