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.

The distribution of visible matter in the universe, such as galaxies and galaxy clusters, has its origin in the week fluctuations of density that existed at the epoch of recombination. The hierarchical distribution of the universe, with its galaxies, clusters and super-clusters of galaxies indicates the absence of a natural length scale. In the Newtonian formulation, numerical simulations of a one-dimensional system permit us to precisely follow the evolution of an ensemble of particles starting with an initial perturbation in the Hubble flow. The limitation of the investigation to one dimension removes the necessity to make approximations in calculating the gravitational field and, on the whole, the system dynamics. It is then possible to accurately follow the trajectories of particles for a long time. The simulations show the emergence of a self-similar hierarchical structure in both the phase space and the configuration space and invites the implementation of a multifractal analysis. Here, after showing that symmetry considerations leads to the construction of a family of equations of motion of the one-dimensional gravitational system, we apply four different methods for computing generalized dimensions $D_q$ of the distribution of particles in configuration space. We first employ the conventional box counting and correlation integral methods based on partitions of equal size and then the less familiar nearest-neighbor and k-neighbor methods based on partitions of equal mass. We show that the latter are superior for computing generalized dimensions for indices $q<-1$ which characterize regions of low density.

Concentrations of matter, such as galaxies and galactic clusters, originated as very small density fluctuations in the early universe. The existence of galaxy clusters and super-clusters suggests that a natural scale for the matter distribution may not exist. A point of controversy is whether the distribution is fractal and, if so, over what range of scales. One-dimensional models demonstrate that the important dynamics for cluster formation occur in the position-velocity plane. Here the development of scaling behavior and multifractal geometry is investigated for a family of one-dimensional models for three different, scale-free, initial conditions. The methodology employed includes: 1) The derivation of explicit solutions for the gravitational potential and field for a one-dimensional system with periodic boundary conditions (Ewald sums for one dimension); 2) The development of a procedure for obtaining scale-free initial conditions for the growing mode in phase space for an arbitrary power-law index; 3) The evaluation of power spectra, correlation functions, and generalized fractal dimensions at different stages of the system evolution. It is shown that a simple analytic representation of the power spectra captures the main features of the evolution, including the correct time dependence of the crossover from the linear to nonlinear regime and the transition from regular to fractal geometry. A possible physical mechanism for understanding the self-similar evolution is introduced. It is shown that hierarchical cluster formation depends both on the model and the initial power spectrum. Under special circumstances a simple relation between the power spectrum, correlation function, and correlation dimension in the highly nonlinear regime is confirmed.

We derive analytic solutions for the potential and field in a one-dimensional system of masses or charges with periodic boundary conditions, in other words Ewald sums for one dimension. We also provide a set of tools for exploring the system evolution and show that it's possible to construct an efficient algorithm for carrying out simulations. In the cosmological setting we show that two approaches for satisfying periodic boundary conditions, one overly specified and the other completely general, provide a nearly identical clustering evolution until the number of clusters becomes small, at which time the influence of any size-dependent boundary cannot be ignored. Finally we compare the results with other recent work with the hope of providing clarification over differences these issues have induced. We explain that modern formulations of physics require a well defined potential which is not available if the forces are screened directly.

We study a one dimensional model of gravitational instability in an Einstein-de Sitter universe. Scaling in both space and time results in an autonomous set of coupled Poisson-Vlasov equations for the field and phase space density, and the $N$-body problem. Using dynamical simulation, we find direct evidence of hierarchical clustering. A multi-fractal analysis reveals a bifractal geometry similar to that observed in the distribution of galaxies. To demonstrate the role of scaling, we compare the system to other one dimensional models recently employed to study structure formation. Finally we show that the model yields an estimate of the time of galaxy formation of the correct order.