Three-dimensional numerical model of the formation of large-scale structure in the Universe

A. A. Klypin,S. F. Shandarin
DOI: https://doi.org/10.1093/mnras/204.3.891
IF: 4.8
1983-10-01
Monthly Notices of the Royal Astronomical Society
Abstract:The first results of numerical fully three-dimensional simulations of formation and evolution of the large-scale structure of the Universe are presented. The simulations were carried out in the framework of the adiabatic scenario of galaxy formation. The model contains 323 = 32 768 collisionless particles interacting only gravitationally. Equal mass particles are moving in a collective gravitational field which is smoothed at small scales. Evolution of perturbations is followed in an expanding cosmological model. The numerical procedure starts at zstart ⋍ 5.25, when (δρ/ρ)r0 ⋍ 0.22 and $$r_0=5\enspace h^{-1}_{100}\enspace \text {Mpc}$$ is our cell length. Previous evolution was computed in accord with the approximate non-linear solution by Zeldovich. The initial perturbations of density as well as velocity were taken as a Fourier sum with random coefficients ranging from $$\lambda_\text {min}=40\enspace h^{-1}_{100}\enspace \text {Mpc}$$ up to $$\lambda_\text {max}=160\enspace h^{-1}_{100}\enspace \text {Mpc}$$. The latter is the size of the cube where particles travel. The total number of cells to solve the Poisson equation is 323, thus the shortest wavelength in the initial perturbations contains eight cells. The structure formation begins from pancakes, then evolves to the network structure, then to formation of clumps connected by strings and afterwards to formation of huge isolated clumps. Because of statistical initial perturbations one can observe different phases of the process in different places simultaneously. The two-point correlation function was calculated. The model correlation function is shown to agree reasonably over the range $$1\gtrsim\xi\gtrsim0.05$$ with the observational one. We were mainly interested in the scale $$5\enspace h^{-1}_{100}\enspace \text {Mpc}\lesssim r\lesssim160\enspace h^{-1}_{100}\enspace \text {Mpc}$$; the smaller as well as larger scales were beyond the opportunity of our simulations. The mass function of clumps agree reasonably with that of Abell clusters.
astronomy & astrophysics
What problem does this paper attempt to address?