We describe a gamma ray pulsar code which computes the high energy photon emissivities from vacuum gaps in the outer magnetosphere, after the model outlined by Cheng, Ho and Ruderman (1986) and Ho (1989). Pair‐production due to photon‐photon interactions and radiation processes including curvature, synchrotron and inverse Compton processes are computed with an iterative scheme which converges to self‐consistent photon and particle distributions for a sampling of locations in the outer magnetosphere. We follow the photons from these distributions as they propagate through the pulsar magnetosphere toward a distant observer. We include the effects of relativistic aberration, time‐of‐flight delays and reabsorption by photon‐photon pair‐production to determine an intensity map of the high energy pulsar emission on the sky. Using data from radio and optical observations to constrain the geometry of the magnetosphere as well as the possible observer viewing angles, we derive light curves and phase dependent spectra which can be directly compared to data from the Compton Observatory. Observations for Crab, Vela and the recently identified gamma ray pulsars Geminga, PSR1706−44 aNd PSR 1509−58 will provide important tests of our model calculations, help us to improve our picture of the relevant physics at work in pulsar magnetospheres and allow us to comment on the implications for future pulsar discoveries.

This content is only available via PDF.
You do not currently have access to this content.