# Pulsar Death at an Advanced Age

## 1. Introduction

Radio emission from Rotation Powered Pulsars (RPPs) probably has its origin in the relativistic outflow of electron-positron pairs along the polar magnetic field lines of a dipole magnetic field frozen into the rotating neutron star (e.g., Arons 1992, Meszaros 1992). The evidence for dipole magnetic fields comes primarily from the electromagnetic theory of RPP spindown energy losses, which occur at the rate $`\dot{E}_R=k\mu ^2\mathrm{\Omega }_{}^4/c^3=I\mathrm{\Omega }_{}\dot{\mathrm{\Omega }}_{}`$ (Dyson 1971, Arons 1979, 1992). Here $`\mu `$ is the magnetic moment, $`\mathrm{\Omega }_{}=2\pi /P`$, $`P`$ is the rotation period, and $`k`$ is a function of any other parameters of significance, with magnitude on the order of unity. In the vacuum theory (Deutsch 1955), $`k=(2/3)\mathrm{sin}^2i`$, with $`i`$ the angle between the magnetic moment and the angular velocity. Theoretical work on the torques due to conduction currents steming back to Goldreich and Julian (1969), coupled to the approximate independence of spindown torques from observationally estimated values of $`i`$ (Lyne and Manchester 1988), suggest that in reality $`k`$ does not substantially depend on $`i`$. In the subsequent discussion, I assume $`k=4/9`$, the average of the vacuum value over the sphere. Application of this EM energy loss rate to the observations of RPPs’ periods and period derivatives yields $`\mu 10^{30}`$ cgs for “normal” RPPs, and $`\mu 10^{27}`$ cgs for millisecond RPPs.

The electromagnetic torque interpretation of pulsar spindown constrains only the exterior dipole moment of the magnetic field. However, Rankin (1990) and Kramer et al. (1998) have presented strong evidence in favor of a low altitude ($`rR_{}`$) dipole geometry for the site of the core component of pulsar radio emission. Arons (1993) gave evidence that spun up millisecond pulsars must have substantially dipolar fields at low altitute.

If electron-positron pair creation above the polar caps is important for radio emission, all observed pulsars must lie in the region of $`P\dot{P}`$ space where polar cap acceleration has sufficient vigor to lead to copious pair production. Yet, to date, all internally consistent theories of polar cap pair creation have required hypothesizing a large scale (e.g., quadrupole) component of the magnetic field with strength comparable to that of the dipole (Ruderman and Sutherland 1975, Arons and Scharlemann 1979, Barnard and Arons 1982, Gurevich and Istomin 1985). Such strong magnetic anomalies contradict the evidence in favor of an apparently dipolar low altitude geometry; the alteration of the magnetic geometry also ruins the internal consistency of many models’ electrodynamics.

Here I describe a low altitude polar cap acceleration theory which successfully associates pulsar “death” with the cessation of pair creation in an offset dipole low altitude magnetic field. The basic acceleration physics is that of a space charge limited relativistic particle beam accelerated along the field lines by the starvation electric field, as in the Arons and Scharlemann theory, but with the additional effect of inertial frame dragging, first pointed out by Muslimov and Tsygan (1990, 1992) and by Beskin (1990). If the dipole’s center is offset from the stellar center, the magnetic field at one pole becomes substantially stronger than it would be if the same magnetic dipole were star centered. The location of an individual pulsar’s pair creation death depends on the magnitude of the offset, thus yielding a “death valley” (Chen and Ruderman 1993) for the whole pulsar population. Finally, if thermal photon emission at temperature $`T10^5`$ K continues to great age ($`t>10^8`$ years), as is is expected in neutron star models with late time heating due to friction between the crust and core, the theory with dipole offsets predicts and accounts for pulsars with very long periods and great age.

## 2. Polar Acceleration

Study of polar cap relativstic particle acceleration in the 1970’s had led to the conclusion that acceleration of a space charge limited particle beam from the stellar surface with energy/particle high enough to emit magnetically convertible curvature gamma rays occurs because of curvature of the magnetic field (Scharlemann et al. 1978, Arons and Scharlemann 1979). In a curved $`B`$ field, matching of the beam charge density to the Goldreich-Julian density occurs only at the surface. Along field lines which curve toward the rotation axis (“favorably curved” field lines), the beam fails to short out the vacuum above the surface, Therefore, particles accelerate along $`B`$ through a potential drop $`\mathrm{\Delta }\mathrm{\Phi }_{}=\mathrm{\Delta }\mathrm{\Phi }_{SAF}\mathrm{\Phi }_{\mathrm{pole}}(R_{}/\rho _B)10^2P^{1/2}\mathrm{\Phi }_{\mathrm{pole}}.`$ The numerical value assumes field lines have dipolar radius of curvature $`\rho _B\sqrt{R_{}c/\mathrm{\Omega }_{}}`$. Here $`\mathrm{\Phi }_{\mathrm{pole}}\mathrm{\Omega }_{}^2\mu /c^2=1.09\times 10^{13}(I_{45}/k)^{1/2}(\dot{P}_{15}/P^3)^{1/2}`$ Volts, with $`\dot{P}_{15}\dot{P}/10^{15}\mathrm{s}/\mathrm{s}`$ and $`I_{45}=I/10^{45}`$ g-cm<sup>2</sup>. Particles drop through the potential $`\mathrm{\Delta }\mathrm{\Phi }_{SAF}`$ over a length $`L_{}R_{}`$.

Curvature gamma rays have typical energy $`\epsilon _c(\mathrm{}c/\rho _B)(e\mathrm{\Delta }\mathrm{\Phi }_{}/mc^2)^3\mathrm{\Phi }_{\mathrm{pole}}^3/\rho _B^4`$, while the optical depth for pair creation, due to one photon conversion of gamma rays emitted by electrons (or positrons) can be shown to be $`\tau =\mathrm{\Lambda }\mathrm{exp}[a(mc^2/\epsilon _c)(B_q/B_{})(\rho _B/R_{})]`$ (Arons and Scharlemann 1979, Luo 1996, Bjornsson 1996), where $`a`$ is a pure number (typically $`30`$) and $`\mathrm{\Lambda }`$ is a combination of the basic parameters which is quite large ($`\mathrm{ln}\mathrm{\Lambda }20`$). A reasonable theoretical definition of the death line is $`\tau =1`$. Using $`B_{}=2(\mathrm{\Phi }_{\mathrm{pole}}/R_{})(c/\mathrm{\Omega }_{}R_{})^2`$, $`\mathrm{\Delta }\mathrm{\Phi }_{SAF}`$ and setting $`\tau `$ equal to unity yields the death line, expressed as $`\mathrm{\Phi }_{\mathrm{death}}(P)`$ such that stars $`\mathrm{\Phi }_{\mathrm{pole}}<\mathrm{\Phi }_{\mathrm{death}}`$ do not make pairs. This death line, appears as the dashed line in Figure 1. This figure shows clearly that the large dynamic range in $`\mathrm{\Phi }_{pole},P`$ space made available by the cataloging of millisecond pulsars falsifies this theory, even if one invokes non-dipolar radii of curvature to move the position of the death line vertically in the diagram - the scaling with period flatly disagrees with the shape of the boundary of pulsar radio emission in the $`\mathrm{\Phi }_{pole},P`$ diagram. These results imply either that something else governs the low altitude acceleration which leads to pair creation, or that pair creation is not important to radio emission.

Muslimov and Tsygan (1990, 1992) uncovered a previously overlooked effect on the acceleration of the non-neutral beam from the stellar surface. Stellar rotation drags the inertial frame into rotation, at the angular velocity $`\omega _{LT}=(2GI/R_{}^3c^2)\mathrm{\Omega }_{}(R_{}/r)^3`$. Therefore, the electric field required to bring a charged particle into corotation is $`𝐄_{co}=(1/c)[(𝛀_{}\omega _{LT})\times 𝐫]\times 𝐁`$; the magnetic field rotates with respect to local inertial space, not inertial space at infinity. The charge density required to support this local corotation electric field is $`\eta _R=[(𝛀_{}\omega _{LT})𝐁]/2\pi c=[𝛀_{}𝐁/2\pi c][1\kappa _g(R_{}/r)^3],`$ where $`\kappa _g=2GI/R_{}^3c^2=0.17(I_{45}/R_{10}^3)`$. Relativistic space charge limited flow from the surface has a beam charge density $`\eta _b=(𝛀_{}𝐁_{}/2\pi c)(1\kappa _g)(B/B_{})`$. Above the surface, this charge density is too small to short out $`E_{}`$ on all polar field lines, not just the favorably curved part of a polar flux tube, thus providing a theoretical basis for polar cap acceleration models to be in accord with the observed rough symmetry of radio emission with respect to the magnetic axis (Lyne and Manchester 1988). One can graphically describe this general relativistic origin of electrical starvation as the consequence of the field lines rotating faster with respect to inertial space as the radius increases, at the angular speed $`\mathrm{\Omega }_{}\omega _{LT}(r)=\mathrm{\Omega }_{}[1\kappa _g(R_{}/r)^3]`$. The constraint of relativistic flow along $`B`$ allows the beam to provide only a charge density sufficient to support corotation at the angular speed $`\mathrm{\Omega }_{}(1\kappa _g)`$. The difference not surprisingly leads to an accelerating potential drop $`\mathrm{\Delta }\mathrm{\Phi }_{}\kappa _g\mathrm{\Phi }_{pole}[1(R_{}/r)^3].`$ For normal ($`P1`$ sec) pulsars with dipole fields, the effect of dragging of inertial frames on the beam’s acceleration yields curvature gamma ray energies 1000 times greater than occur in the Arons and Scharlemann pair creation theory; for MSPs, the theories yield comparable results, although of course the symmetry of the beam with respect to the magnetic axis differs.

## 3. Death Lines and Death Valley

When curvature emission is the only source of gamma rays, the death line for a star centered dipole is shown in Figure 1. Dragging of inertial frames clearly improves the agreement between the boundary of pair activity in the $`\mathrm{\Phi }P`$ diagram and the region where pulsars occur, but the discrepancy is still too large - something else is missing. If the field geometry must be locally dipolar at low altitude, then the only ingredients still not included are 1) offset of the dipole from the stellar center and 2) additional gamma ray emission and absorption processes. The simplest dipole offset has the magnetic field of a point dipole, with the center of the dipole displaced from the stellar center by an offset vector parallel to $`\mu `$. This has the effect of increasing the magnetic field at one pole to strength $`B_{}=2\mu /(R_{}\delta )^3`$, with a resulting drastic increase in the gamma ray opacity, while leaving the accelerating potential unaltered. The results with curvature radiation as the only gamma ray emission process show that dipole offsets do allow such a curvature radiation theory of pulsar death to survive the challenge of modern observations, although at the price of displacements of the dipole center from the stellar center comparable to moving the dipole’s center to the base of the crust. Such a model still does not account for PSR J2144-3933 (Young et al. 1999), however.

Curvature emission is not the only means of converting beam energy to gamma rays. ROSAT observations have revealed the long sought thermal X-rays from neutron star surfaces (Becker and Trümper 1997). Resonant Compton scattering creates magnetically convertible gamma rays at a spatial rate $`(dN_\gamma /ds)_{rC}T/\mathrm{\Gamma }^2`$, (e.g., Luo 1996) where $`T`$ is the temperature of the cooling neutron star (polar cap heating is unimportant in death valley) and $`\mathrm{\Gamma }=e\mathrm{\Phi }/m_\pm c^2`$ is the Lorentz factor of an electron or positron in the beam. Compton scattering thus can become a significant source of gamma rays in stars with small accelerating potentials. Compton scattering thus may contribute significantly for stars with low overall voltage.

This expectation is correct, if internal heating (e.g., Umeda et al. 1993) keeps the surface temperature above $`10^5`$ K at spindown ages in excess of $`10^{7.5}`$ years. In this case, resonant Compton scattering of thermal photons by a polar electron beam does extend death valley to include all the observed pulsars, with offsets required in the lowest voltage RPPs of order 60% to 70%, as is shown in the right panel of Figure 1. Note that this theory predicted (Arons 1998) the existence of RPPs with large periods ($`P>5`$ seconds) and unusually low voltage $`\mathrm{\Phi }_{pole}<10^{12.5}`$ Volts \[great age, $`t=170(10^{12}V/\mathrm{\Phi }_{pole})^2(10^s/P)^2`$ Myr.\]

Indeed, recent observations (Young et al. 1999) have found a a $`P=8.51`$ second pulsar, PSR J2144-3933, which falls below the death line for a star centered dipole (even including Compton scattering), but lies comfortably within death valley, when both Compton scatering and dipole offsets are included.

## 4. Conclusion

I have shown that polar pair creation based on acceleration of a steadily flowing, space charged limited non-neutral beam in a locally dipolar magnetic geometry at low altitude is consistent with pulsar radio emission throughout the $`P\dot{P}`$ diagram, provided 1) the effect of dragging of inertial frames is included in estimates of the starvation electric field; 2) the dipole center is strongly offset from the stellar center in older stars, perhaps as much as $`0.70.8R_{}`$; and 3) inverse Compton emission of thermal photons from a neutron star cooling slower than exponentially at ages in excess of $`10^6`$ years plays an important role in the emission of magnetically convertible gamma rays. Stars not deep in death valley have copious pair outflow on all polar field lines. Such outflow shorts out the “outer gaps” proposed as the dynamical basis for gamma ray emission in the outer magnetosphere (Cheng, Ho and Ruderman 1986, Romani 1996). Either polar cap gamma ray emission (e.g. Zhang and Harding 2000) or possible outer magnetsophere emission from a dissipative return current layer (Arons 1981) remain as candidates for the gamma ray emission observed in RPPs.

The development of new diagnostics of the low altitude magneic field, gamma ray observations sensitive to low altitude emission, and optical and UV observations of thermal emission from old, nearby RPPs will eventually provide tests of these ideas.

## 5. Acknowledgments

My research on pulsars is supported in part by NSF grant AST 9528271 and NASA grant NAG 5-3073, and in part by the generosity of California’s taxpayers.

## References

Arons, J. 1979, Space Sci. Rev., 24, 437

. 1981, in Proc. IAU Symp. No. 95 ‘Pulsars’, W. Sieber and R. Wielebinski, eds. (Dordrecht: Reidel), 69

. 1992, Proc. IAU Colloq. No. 128, ‘The Magnetospheric Structure and Emission Mechanisms of Raio Pulsars’, T.H. Hankins, J.M. Rankin and J. A. Gil, eds. (Zielona Gora: Pedagogical University Press), 56

. 1993, Ap.J., 408, 160

. 1998, in ‘Neutron Stars and Pulsars: Thirty Years After the Discovery’, N. Shibazaki et al., eds (Tokyo: Universal Academy Press), 339 (astro-ph/9802198)

Arons, J., and Scharlemann, E.T. 1979, Ap.J., 231, 854

Beskin, V.S. 1990, Pis’ma Ast. Zh., 16, 665 (Sov. Ast. - Letters, 16)

Barnard, J.J., and Arons, J. 1982, Ap.J., 254, 713

Becker, W., and Trümper, J. 1997, A&A, 326, 682

Bjornsson, C.-I. 1996, Ap.J., 471, 321

Chen, K., and Ruderman, M.A. 1993, Ap.J., 402, 264

Cheng, K.S., Ho., C., and Ruderman, M. 1986, Ap.J., 300, 500

Deutsch, A.J. 1955, Ann. Ap., 18, 1

Dyson, F.J. 1971, “Neutron Stars and Pulsars: Fermi Lectures 1970” (Rome: Academia Nazionale dei Lincei), 25-26

Goldreich, P. and Julian, W. 1969, Ap.J., 157, 869

Gurevich, A.V., and Istomin, Ya.N. 1985, Zh.Eksp.Teor.Fiz., 89,3 (Soviet Physis - JETP, 62, 1)

Kramer, M., et al. 1998, Ap.J., 501, 270

Luo, Q. 1996, Ap.J., 468, 338

Lyne, A.G., and Manchester, R.N. 1988, MNRAS, 234, 477

Meszaros, P. 1992, “High Energy Radiation from Magnetized Neutron Stars” (Chicago: University of Chicago Press)

Muslimov, A., and Tsygan, A.I. 1990, Ast. Zh., 67, 263 (Soviet Ast., 34, 133)

. 1992, MNRAS, 255, 61

Rankin, J.M. 1990, Ap.J., 352, 247

Romani, R. 1996, Ap.J., 470, 469

Ruderman, M.A., and Sutherland, P.G. 1975, Ap.J., 196, 51

Scharlemann, E.T., Arons, J., and Fawley, W.M. 1978, Ap.J., 222, 297

Umeda, H., et al. 1993, Ap.J., 408, 186

Young, M. et al. 1999, Nature, 400, 848

Zhang, B., and Harding, A.K. 2000, Ap.J., in press (astro-ph/9911028)