Spin evolution and mass distribution of the Galactic Binary Neutron Stars

TL;DR

This study models the spin evolution and mass distribution of Galactic binary neutron stars, explaining observed properties and predicting a significant fraction of high-mass mergers detectable via gravitational waves.

Contribution

It introduces an improved binary star evolution model that accurately reproduces observed BNS properties and predicts the mass distribution of merging BNSs, including high-mass systems.

Findings

Low-mass BNSs can spin up to millisecond periods via accretion.

Most high-mass BNSs have shorter radio lifetimes, reducing detection probability.

Approximately 19-22% of merging BNSs have masses exceeding 3 solar masses.

Abstract

Binary neutron stars (BNSs) detected in the Milky Way have the total masses distributing narrowly around $\sim2.6-2.7M_\odot$, while the BNS merger GW190425 detected via gravitational wave has a significantly larger mass ($\sim3.4M_\odot$). This difference is not well understood, yet. In this paper, we investigate the BNS spin evolution via an improved binary star evolution model and its effects on the BNS observability, with implementation of various relevant astrophysical processes. We find that the first-born neutron star component in low-mass BNSs can be spun up to millisecond pulsars by the accretion of Roche-lobe overflow from its companion and its radio lifetime can be comparable to the Hubble time. However, most high-mass BNSs have substantially shorter radio lifetime than the low-mass BNSs, and thus smaller probability being detected via radio emission. Adopting the star formation and metal enrichment history of the Milky Way given by observations, we obtain the survived Galactic BNSs with pulsar components from our population synthesis model and find that their distributions on the diagrams of spin period versus spin-period-time-derivative ($P-\dot{P}$) and orbital period versus eccentricity ($P_{\rm orb}-e$) can well match those of the observed Galactic BNSs. The total mass distribution of the observed Galactic BNSs can also be matched by the model. A significant fraction ($\sim19\%-22\%$) of merging BNSs at redshift $z\sim0$ have masses $\gtrsim3M_\odot$, which seems compatible with the GW observations. Future radio observations may detect many more Galactic BNSs, which will put strong constraint on the spin evolution of BNSs during their formation processes.