We present numerical simulations of the electron-positron plasma creation process in a simple neutron star magnetosphere. We have developed a set of cascade `kernels', which represent the endpoint of the pair cascades resulting from monoenergetic photon seeds. We explore two popular models by convolving these kernels with the seed photon distributions produced by curvature radiation and by inverse Compton scattering. We find that the pair plasma in either case is well-described in its rest frame by a relativistic Maxwellian distribution with temperature near mc^2/k_B. We present cascade multiplicities and efficiencies for a range of seed particle energies and stellar magnetic fields. We find that the efficiencies and multiplicities of pair creation are often lower than has been assumed in previous work.

The pulsar magnetosphere contains a strongly magnetized, relativistic plasma. We need to understand the physics of that plasma if we want to connect the data to the models. Our group in Socorro is mixing theory and observations in order to study the radio-loud pulsar plasma. In this paper we report on several aspects of our current work.