Cosmology is entering an era of percent level precision due to current large observational surveys. This precision in observation is now demanding more accuracy from numerical methods and cosmological simulations. In this paper, we study the accuracy
of $N$-body numerical simulations and their dependence on changes in the initial conditions and in the simulation algorithms. For this purpose, we use a series of cosmological $N$-body simulations with varying initial conditions. We test the influence of the initial conditions, namely the pre-initial configuration (preIC), the order of the Lagrangian perturbation theory (LPT), and the initial redshift, on the statistics associated with the large scale structures of the universe such as the halo mass function, the density power spectrum, and the maximal extent of the large scale structures. We find that glass or grid pre-initial conditions give similar results at $zlesssim 2$. However, the initial excess of power in the glass initial conditions yields a subtle difference in the power spectra and the mass function at high redshifts. The LPT order used to generate the ICs of the simulations is found to play a crucial role. First-order LPT (1LPT) simulations underestimate the number of massive haloes with respect to second-order (2LPT) ones, typically by 2% at $10^{14} h^{-1} M_odot$ for an initial redshift of 23, and the small-scale power with an underestimation of 6% near the Nyquist frequency for $z_mathrm{ini} = 23$. Moreover, at higher redshifts, the high-mass end of the mass function is significantly underestimated in 1LPT simulations. On the other hand, when the LPT order is fixed, the starting redshift has a systematic impact on the low-mass end of the halo mass function.
Large galaxy redshift surveys have long been used to constrain cosmological models and structure formation scenarios. In particular, the largest structures discovered observationally are thought to carry critical information on the amplitude of large
-scale density fluctuations or homogeneity of the universe, and have often challenged the standard cosmological framework. The Sloan Great Wall (SGW) recently found in the Sloan Digital Sky Survey (SDSS) region casts doubt on the concordance cosmological model with a cosmological constant (i.e. the flat LCDM model). Here we show that the existence of the SGW is perfectly consistent with the LCDM model, a result that only our very large cosmological N-body simulation (the Horizon Run 2, HR2) could supply. In addition, we report on the discovery of a void complex in the SDSS much larger than the SGW, and show that such size of the largest void is also predicted in the LCDM paradigm. Our results demonstrate that an initially homogeneous isotropic universe with primordial Gaussian random phase density fluctuations growing in accordance with the General Relativity, can explain the richness and size of the observed large-scale structures in the SDSS. Using the HR2 simulation we predict that a future galaxy redshift survey about four times deeper or with 3 magnitude fainter limit than the SDSS should reveal a largest structure of bright galaxies about twice as big as the SGW.
In support of the new Sloan III survey, which will measure the baryon oscillation scale using the luminous red galaxies (LRGs), we have run the largest N-body simulation to date using $4120^3 = 69.9$ billion particles, and covering a volume of $(6.59
2 h^{-1} {rm Gpc})^3$. This is over 2000 times the volume of the Millennium Run, and corner-to-corner stretches all the way to the horizon of the visible universe. LRG galaxies are selected by finding the most massive gravitationally bound, cold dark matter subhalos, not subject to tidal disruption, a technique that correctly reproduces the 3D topology of the LRG galaxies in the Sloan Survey.We have measured the covariance function, power spectrum, and the 3D topology of the LRG galaxy distribution in our simulation and made 32 mock surveys along the past lightcone to simulate the Sloan III survey. Our large N-body simulation is used to accurately measure the non-linear systematic effects such as gravitational evolution, redshift space distortion, past light cone space gradient, and galaxy biasing, and to calibrate the baryon oscillation scale and the genus topology. For example, we predict from our mock surveys that the baryon acoustic oscillation peak scale can be measured with the cosmic variance-dominated uncertainty of about 5% when the SDSS-III sample is divided into three equal volume shells, or about 2.6% when a thicker shell with $0.4<z<0.6$ is used. We find that one needs to correct the scale for the systematic effects amounting up to 5.2% to use it to constrain the linear theories. And the uncertainty in the amplitude of the genus curve is expected to be about 1% at 15 $h^{-1}$Mpc scale. We are making the simulation and mock surveys publicly available.
