Evidence for a maximum mass cut-off in the neutron star mass distribution and constraints on the equation of state. Issue 1 (26th April 2018)

• Record Type: Journal Article
• Authors: Alsing, Justin; Silva, Hector O; Berti, Emanuele
• Is Part Of: Monthly notices of the Royal Astronomical Society. Volume 478:Issue 1(2018)
• Page Start: 1377
• Page End: 1391
• Publication Date: 2018-04-26
• Subjects: equation of state -- stars: neutron; Astronomy -- Periodicals; Periodicals; 520.5
• DOI: 10.1093/mnras/sty1065
• Languages: English
• ISSNs: 0035-8711

Abstract

Abstract: We infer the mass distribution of neutron stars in binary systems using a flexible Gaussian mixture model and use Bayesian model selection to explore evidence for multimodality and a sharp cut-off in the mass distribution. We find overwhelming evidence for a bimodal distribution, in agreement with previous literature, and report for the first time positive evidence for a sharp cut-off at a maximum neutron star mass. We measure the maximum mass to be 2.0 M⊙ < m max < 2.2 M⊙ (68 per cent), 2.0 M⊙ < m max < 2.6 M⊙ (90 per cent), and evidence for a cut-off is robust against the choice of model for the mass distribution and to removing the most extreme (highest mass) neutron stars from the data set. If this sharp cut-off is interpreted as the maximum stable neutron star mass allowed by the equation of state of dense matter, our measurement puts constraints on the equation of state. For a set of realistic equations of state that support >2 M⊙ neutron stars, our inference of m max is able to distinguish between models at odds ratios of up to 12:1, whilst under a flexible piecewise polytropic equation-of-state model our maximum mass measurement improves constraints on the pressure at 3–7× the nuclear saturation density by ${\sim } 30\text{--}50\hbox{ per cent}$ compared to simply requiring m max > 2 M⊙ . We obtain a lower bound on the maximum sound speed attained inside the neutron star of $c_\mathrm{ s}^\mathrm{max} \gt 0.63c$ (99.8 per cent), ruling out $c_\mathrm{ s}^\mathrm{max} \lt c/\sqrt{3}$ at high significance. Our constraints on the maximum neutron star mass strengthen the case for neutron star–neutron star mergers as the primary source of short gamma-ray bursts.