Evolution of bias in different cosmological models

We study the evolution of the halo-halo correlation function and small-scale (approximate to 0.2-7 h(-1) Mpc) bias in four cosmological models (Lambda CDM, OCDM, tau CDM, and SCDM) using very high resolution n-body simulations with a dynamical range of similar to 10,000-32,000 (force resolution of approximate to 2-4 h(-1) kpc and particle mass of approximate to 10(9) h(-1) M.). The high force and mass resolution allows dark matter halos to survive in the tidal fields of high-density regions and thus prevents the ambiguities related with the "overmerging problem." This allows us to estimate for the first time the evolution of the correlation function and bias at small (down to similar to 100 h(-1) kpc) scales. We find that at all epochs the two-point correlation function of galaxy-size halos xi(hh) is well approximated by a power law with slope approximate to 1.6-1.8. The difference between the shape of xi(hh) and the shape of the correlation function of matter results in the scale-dependent bias at: scales less than or similar to 7 h(-1) Mpc, which we find to be a generic prediction of the hierarchical models, independent of the epoch and of the model details. The bias evolves rapidly from a high value of similar to 2-5 at z similar to 3-7 to the antibias of b similar to 0.5-1 at small less than or similar to 5 h(-1) Mpc scales at z = 0. Another generic prediction is that the comoving amplitude of the correlation function for halos above a certain mass evolves non-monotonically: it decreases from an initially high value at z similar to 3-7, and very slowly increases at z less than or similar to 1. We find that our results agree well with existing clustering data at different redshifts, indicating the general success of the hierarchical models of structure formation in which galaxies form inside the host DM halos. Particularly, we find an excellent agreement in both slope and the amplitude between xi(hh)(z = 0) in our ACDM(60) simulation and the galaxy correlation function measured using the Automatic Plate Measuring Facility galaxy survey. At high redshifts, the observed clustering of the Lyman-break galaxies is also well reproduced by the models. We find good agreement at z greater than or similar to 2 between our results and predictions of the analytical models of bias evolution. This indicates that we have a solid understanding of the nature of the bias and of the processes that drive its evolution at these epochs. We argue, however, that at lower redshifts the evolution of the bias is driven by dynamical processes inside the nonlinear high-density regions such as galaxy clusters and groups. These processes do not depend on cosmology and tend to erase the differences in clustering properties of halos that exist between cosmological models at high z.