On the pulse--width statistics in radio pulsars. I. Importance of the interpulse emission

TL;DR

This study uses Monte Carlo simulations and large pulsar data sets to analyze pulse-width statistics, interpulse emission rates, and magnetic inclination angles, revealing insights into pulsar beam geometry and population estimates.

Contribution

It introduces a comprehensive simulation approach combining new data to constrain pulsar emission models and magnetic inclination distributions, highlighting the importance of interpulse emission.

Findings

The parent period distribution is best modeled by a log-normal function.

The inclination angle distribution has two peaks near 0° and 90°.

The pulsar beam covers about 10% of the sky, suggesting many neutron stars are undetected.

Abstract

We performed Monte Carlo simulations of different properties of pulsar radio emission, such as: pulsar periods, pulse-widths, inclination angles and rates of occurrence of interpulse emission (IP). We used recently available large data sets of the pulsar periods P, the pulse profile widths W and the magnetic inclination angle alpha. We also compiled the largest ever database of pulsars with interpulse emission, divided into the double-pole (DP-IP) and the single-pole (SP-IP) cases. Their distribution on the P - Pdot diagram strongly suggests a secular alignment of the magnetic axis from the originally random orientation. We derived possible parent distribution functions of important pulsar parameters by means of the Kolmogorov-Smirnov significance test using the available data sets (P, W, alpha and IP), different models of pulsar radio beam rho = rho(P) as well as different trial distribution functions of pulsar period and the inclination angles. The best suited parent period distribution function is the log-normal distribution, although the gamma function distribution cannot be excluded. The strongest constraint on derived model distribution functions was the requirement that the numbers of interpulses were exactly (within 1sigma errors) at the observed level of occurrences. We found that a suitable model distribution function for the inclination angle is the complicated trigonometric function which has two local maxima, one near 0 deg and the other near 90 deg. The former and the latter implies the right rates of IP occurrence. It is very unlikely that the pulsar beam deviates significantly from the circular cross-section. We found that the upper limit for the average beaming factor fb describing a fraction of the full sphere (called also beaming fraction) covered by a pulsar beam is about 10%. This implies that the number of the neutron stars in the Galaxy might be underestimated.