A numerical investigation of the steady states of the spherically symmetric Einstein–Vlasov–Maxwell system

We construct, by numerical means, static solutions of the spherically symmetric Einstein-Vlasov-Maxwell system and investigate various features of the solutions. This extends a previous investigation (Andreasson and Rein 2007 Class. Quantum Grav. 24 1809) of the chargeless case. We study the possible shapes of the energy density profile as a function of the area radius when the electric charge of an individual particle is varied as a parameter. We find profiles which are multi-peaked, where the peaks are separated either by vacuum or a thin atmosphere, and we find that for a sufficiently large charge parameter the solutions break down at a finite radius. Furthermore, we investigate the inequality