The statistics of gamma-ray pulsars

We use Monte Carlo methods to simulate the properties of a Galactic population of rotation-powered pulsars in the self-consistent outer gap model proposed by Zhang & Cheng, where the initial magnetic field and spatial and velocity distributions of the neutron stars at birth are obtained from the statistical results of radio pulsars. We obtain the distance, period, age, magnetic field, and gamma-ray flux distributions of the gamma-ray pulsars whose gamma-ray and radio fluxes are above the detectable threshold fluxes. Furthermore, we simulate the properties of Geminga-like pulsars and obtain the possible parametric region and the number of the Geminga-like pulsars. In our simulations, the different beaming effects of radio and gamma-ray beams are taken into account. We predict that there may be similar to 11 gamma-ray pulsars that are detectable at both radio and gamma-ray energies and similar to 37 Geminga-like pulsars in the Galactic plane (\b\ < 5 degrees) if the gamma-ray beaming fraction is 0.67 and the birthrate of the neutron stars is 1 per century. Compared to the flux distribution, individual spectrum, and distance distribution obtained by assuming OB star associations of our simulated Geminga-like pulsars and the unidentified EGRET sources at \b\ < 5 degrees, we conclude that the majority of the unidentified gamma-ray point sources near the Galactic plane may be Geminga-like pulsars. We also suggest that the most possible radio pulsars likely to be confirmed as gamma-ray pulsars in future are those with L-X/(E) over dot (sd) similar to 10(-3).