  It is shown that Milne models (a subclass of FLRW spacetimes with negative
spatial curvature) are nonlinearly stable in the set of solutions to the
Einstein-Vlasov-Maxwell system, describing universes with ensembles of
collisionless self-gravitating, charged particles. The system contains various
slowly decaying borderline terms in the mutually coupled equations describing
the propagation of particles and Maxwell fields. The effects of those terms are
controlled using a suitable hierarchy based on the energy density of the matter
fields.
