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.

We investigate structure formation in a one dimensional model of a matter-dominated universe using a quasi-newtonian formulation. In addition to dark matter, luminous matter is introduced to examine the potential bias in the distributions. We use multifractal analysis techniques to identify structures, including clusters and voids. Both dark matter and luminous matter exhibit fractal geometry as the universe evolves over a finite range. We present the results for the generalized dimensions computed on various scales for each matter distribution. We compare and contrast the fractal dimensions of two types of matter for the first time and show how dynamical considerations cause them to differ.

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.