MONTE CARLO SIMULATION OF NEUTRON STAR MASSES

Abstract

We perform a Monte Carlo simulation to investigate the distribution and evolution of 66 mass measured pulsars. We get the best fits to the distribution at 1.35 ± 0.27M⨀ (1σ confidence level). In addition, we notice bimodal distributions in 1.34 ± 0.15M⨀ and 1.48 ± 0.53M⨀, this can be led to the idea that radio pulsars in binary systems have recycled. Thus we divide the data according to the characteristic spin period into two groups, millisecond pulsars (MSPs), P Spin ≤ 20 ms and less recycled pulsars P Spin ≥ 20 ms , respectively. We show that the distributions of MSPs at 1.42 ± 0.36M⨀, and 1.32 ± 0.18M⨀ for less recycled pulsars. As such, the mass of MSPs are heavier than those in less recycled pulsars by ~ 0.1M⨀, since they accreting material from their companions. On the other hand, the formation of heavier pulsars from the accretion induced collapse of accreting white dwarfs, must be invoked.