A new algorithm for implementing the adaptive Monte Carlo method is given. It is used to solve the relativistic Boltzmann equations that describe the time evolution of a nonequilibrium electron-positron pair plasma containing high-energy photons and pairs. The collision kernels for the photons as well as pairs are constructed for Compton scattering, pair annihilation and creation, bremsstrahlung, and Bhabha & Moller scattering. For a homogeneous and isotropic plasma, analytical equilibrium solutions are obtained in terms of the initial conditions. For two non-equilibrium models, the time evolution of the photon and pair spectra is determined using the new method. The asymptotic numerical solutions are found to be in a good agreement with the analytical equilibrium states. Astrophysical applications of this scheme are discussed.

We study the emission and absorption spectra due to various photon and pair processes in a non-equilibrium pair plasma containing a significant density of photons. We present here some preliminary results from Monte-Carlo simulations. These investigations are likely to be useful in understanding the radiation and relaxation mechanisms in non-thermal gamma-ray sources in astrophysics.

Recent observations by the ROSAT X-ray satellite of the Vela supernova remnant have revealed, in addition to the previously identified compact nebula, a nearly circular emitting region with a radius of about 4 degrees. The Vela pulsar is slightly off the center of this circular region, consistent with its measured proper motion of about 100km/s and an age of about 10^4 yr. The emitting region is bounded by the main supernova shock. Just outside the shock, the X-ray image reveals several well-defined V-shaped features extending radially outwards. These features are most likely wakes produced by objects moving supersonically through the outside medium. The shapes and orientations of the wakes suggest that these objects have been ejected from the center of the supernova explosion. Their present positions indicate that they have been moving with a mean velocity of a few thousand km/s. We show that pre-existing objects, such as planets in orbit around the progenitor star, could not have been accelerated to sufficiently high velocities or would have been destroyed. Instead, we propose that the observed objects are fragments ejected during the formation of the neutron star. Fragmentation during gravitational collapse is a natural consequence of both convective and rotational instabilities.