inclusion of the neutralizing density $\phi$ is also necessary in order to arrive at finite energy solutions for (\[RVMFP\]) with a single species of ion.

Over the past twenty-five years significant progress has been made in the analysis of (\[RVM\]), specifically, the global existence of weak solutions (which also holds for the non-relativistic system (VM); see [@DPL]) and the determination of conditions which ensure global existence of classical solutions (originally discovered in [@GlStr], and later in [@BGP] and [@KlStaf] using different methods) for the Cauchy problem. Additionally, a wide array of results have been obtained regarding electrostatic simplifications of (\[RVM\]) - the Vlasov-Poisson and relativistic Vlasov-Poisson systems, obtained by taking the limit as $c \to \infty$ [@Sch86] and $B \equiv 0$, respectively. These models do not include magnetic effects, and the electric field is given by an elliptic equation rather than a hyperbolic PDE. This simplification has led to a great deal of progress concerning the electrostatic systems, including theorems regarding the well-posedness of solutions [@LP; @VPSSA3DGeneral; @Pfaff; @SchVP]. General references on kinetic equations of plasma dynamics, such as (\[RVM\]) and (\[RVMFP\]), include [@Glassey] and [@VKF].

Independent of recent advances, many of the most basic existence and regularity questions remain 