2004ApJ...604..317Z

We reexamine the outer-gap size by taking the geometry of the dipole magnetic field into account. Furthermore, we also consider that instead of taking the gap size at half of the light cylinder radius to represent the entire outer gap, it is more appropriate to average the entire outer-gap size over the distance. When these two factors are considered, the derived outer-gap size f(P,B,<r>(α)) is a function not only of the period P and magnetic field B of the neutron star but also of the average radial distance to the neutron star, <r>, which depends on the magnetic inclination angle α. We use this new outer-gap model to study the γ-ray luminosity of pulsars, which is given by L_γ=f³(P,B,<r>(α))L_sd, where L_sdis the pulsar spin-down power, and to study the death lines of γ-ray emission of the pulsars. Our model can predict the γ-ray luminosity of an individual pulsar if its P, B, and α are known. Since different pulsars have different α, this explains why some γ-ray pulsars have very similar P and B but very different γ-ray luminosities. In determining the death line of γ-ray pulsars, we have used a new criterion based on a concrete physical property, i.e., that the fractional size of the outer gap at the null-charge surface for a given pulsar cannot be larger than unity. In an estimate of the fractional size of the outer gap, two possible X-ray fields are considered: (1) X-rays produced by neutron star cooling and polar-cap heating, and (2) X-rays produced by the bombardment of relativistic particles from the outer gap onto the stellar surface (the outer gap is called a ``self-sustained outer gap''). Since it is very difficult to measure α in general, we use a Monte Carlo method to simulate the properties of γ-ray pulsars in our Galaxy. We find that this new outer-gap model predicts many more weak γ-ray pulsars, which have a typical age between 0.3 and 3 Myr, than does the old model. For all simulated γ-ray pulsars with self-sustained outer gaps, the γ-ray luminosity L_γ satisfies L_γ∝L^δ_sd, where the value of δ depends on the sensitivity of the γ-ray detector. For EGRET, δ∼0.38, whereas δ∼0.46 for GLAST. For γ-ray pulsars with L_sd≲L^crit_sd, δ∼1, and L^crit_sd=1.5x10³⁴P^1/3 ergs/s is determined by f(<r≳r_L)=1. These results are roughly consistent with the observed luminosity of γ-ray pulsars. These predictions are very different from those of the previous outer-gap model, which predicts a very flat relation between L_γand L_sd.