2004ApJ...604..775G

We present results of a pulsar population synthesis study that incorporates a number of recent developments and some significant improvements over our previous study. We have included the results of the Parkes multibeam pulsar survey in our select group of nine radio surveys, doubling our sample of radio pulsars. More realistic geometries for the radio and γ-ray beams are included in our Monte Carlo computer code, which simulates the characteristics of the Galactic population of radio and γ-ray pulsars. We adopted with some modifications the radio-beam geometry of Arzoumanian, Chernoff, and Cordes. For the γ-ray beam, we have assumed the slot gap geometry described in the work of Muslimov and Harding. To account for the shape of the distribution of radio pulsars in the P{dot}-P diagram, we continue to find that decay of the magnetic field on a timescale of 2.8 Myr is needed. With all nine surveys, our model predicts that EGRET should have seen seven radio-quiet (below the sensitivity of these radio surveys) and 19 radio-loud γ-ray pulsars. AGILE (nominal sensitivity map) is expected to detect 13 radio-quiet and 37 radio-loud γ-ray pulsars, while GLAST, with greater sensitivity, is expected to detect 276 radio-quiet and 344 radio-loud γ-ray pulsars. When the Parkes multibeam pulsar survey is excluded, the ratio of radio-loud to radio-quiet γ-ray pulsars decreases, especially for GLAST. The decrease for EGRET is 45%, implying that some fraction of EGRET unidentified sources are radio-loud γ-ray pulsars. In the radio geometry adopted, short-period pulsars are core dominated. Unlike the EGRET γ-ray pulsars, our model predicts that when two γ-ray peaks appear in the pulse profile, a dominant radio core peak appears in between the γ-ray peaks. Our findings suggest that further improvements are required in describing both the radio and γ-ray geometries.