We construct a family of models with negative gravitational mass in the context of Newtonian gravity. We focus in particular on a model that reproduces the features of the so-called Dirac-Milne universe, a matter-antimatter symmetric universe that was recently proposed as an alternative cosmological scenario [A. Benoit-Levy and G. Chardin, A&A 537, A78 (2012)]. We perform one-dimensional N-body simulations of these negative-mass models for an expanding universe and study the associated formation of gravitational structures. The similarities and differences with the standard cosmological model are highlighted and discussed.

Numerical simulations of self-gravitating systems are generally based on N-body codes, which solve the equations of motion of a large number of interacting particles. This approach suffers from poor statistical sampling in regions of low density. In contrast, Vlasov codes, by meshing the entire phase space, can reach higher accuracy irrespective of the density. Here, we performed one-dimensional Vlasov simulations of a long-standing cosmological problem, namely the fractal properties of an expanding Einstein-de Sitter universe in Newtonian gravity. The N-body results were confirmed for high-density regions and extended to regions of low matter density, where the N-body approach usually fails.