Pulsars are highly magnetised rotating neutron stars, emitting in a broad electromagnetic energy range. These objects were discovered more than 55 years ago and are astrophysical laboratories for studying physics at extreme conditions. Reproducing the observed pulsar population helps refine our understanding of their formation and evolution scenarios, as well as their radiation processes and geometry. In this paper, we improve our previous population synthesis by focusing on both the radio and gamma -ray pulsar populations, investigating the impact of the Galactic gravitational potential and of the radio emission death line. To elucidate the necessity of a death line, we implemented our refined initial distributions of the spin period and spacial position at birth. This approach allowed us to elevate the sophistication of our simulations to the most recent state-of-the-art approaches. The motion of each individual pulsar was tracked in the Galactic potential by a fourth-order symplectic integration scheme. Our pulsar population synthesis took into account the secular evolution of the force-free magnetosphere and magnetic field decay simultaneously and self-consistently. Each pulsar was evolved from birth to the present time. The radio and gamma -ray emission locations were modelled by the polar cap geometry and striped wind model, respectively. By simulating ten million pulsars, we found that including a death line allows us to better reproduce the observational trend. However, when simulating one million pulsars, we obtained an even more realistic $P- P $ diagram, whether or not a death line was included. This suggests that the ages of the detected pulsars might be overestimated and so, it sets the need for a death line in pulsar population studies into question. Kolmogorov-Smirnov tests confirm the statistical similarity between the observed and simulated $P- P $ diagram. Additionally, simulations with increased gamma -ray telescope sensitivities hint at a significant contribution coming from the gamma -ray pulsars to the GeV excess in the Galactic centre.
Read full abstract