We investigate the formation of virialized halos due to the gravitational collapse of collisionless matter using high-resolution N-body simulations. A variety of formation scenarios are studied, ranging from hierarchical clustering to monolithic radial collapse. The goal of these experiments was to study departures from the universal density profiles recently found to arise in cosmological settings. However, we found that even for models which exhibit quite a different formation history, the density and velocity dispersion profiles of the virialized halos are strikingly similar. Power law density profiles do not result even in models with initial power law profiles and no initial substructure or non-radial motions. Such initial conditions give rise to a radial orbit instability which leads to curved velocity dispersion and density profiles. The shapes of the density profiles in all our models are well parameterized by the profiles of halos formed in a generic cosmological setting. Our results show that the universality of dark halo density profiles does not depend crucially on hierarchical merging as has been suggested recently in the literature. Rather it arises because apparently different collapse histories produces a near universal angular momentum distribution among the halo particles. We conclude that the density and velocity dispersion profiles of virialized halos in an expanding universe are robust outcomes of gravitational collapse, nearly independent of the initial conditions and the formation history.

The gravitational lens effect of galaxy clusters can produce large arcs from source galaxies in their background. Typical source redshifts of ~ 1 require clusters at z ~ 0.3 for arcs to form efficiently. Given the cluster abundance at the present epoch, the fewer clusters exist at z ~ 0.3 the higher Omega_0 is, because the formation epoch of galaxy clusters strongly depends on Omega_0. In addition, at fixed Omega_0, clusters are less concentrated, and hence less efficient lenses, when the cosmological constant is positive, Omega_Lambda > 0. Numerical cluster simulations show that the expected number of arcs on the sky is indeed a sensitive function of Omega_0 and Omega_Lambda. The numerical results are compatible with the statistics of observed arcs only in a universe with low matter density, Omega_0 ~ 0.3, and zero cosmological constant. Other models fail by one or two orders of magnitude, rendering arc statistics a sensitive probe for cosmological parameters.

We use numerical simulations of galaxy clusters in different cosmologies to study their ability to form large arcs. The cosmological models are: Standard CDM (SCDM; Omega_0=1, Omega_Lambda=0); tauCDM with reduced small-scale power (parameters as SCDM, but with a smaller shape parameter of the power spectrum); open CDM (OCDM; Omega_0=0.3, Omega_Lambda=0); and spatially flat, low-density CDM (LambdaCDM; Omega_0=0.3, Omega_Lambda=0.7). All models are normalised to the local number density of rich clusters. Simulating gravitational lensing by these clusters, we compute optical depths for the formation of large arcs. For large arcs with length-to-width ratio >= 10, the optical depth is largest for OCDM. Relative to OCDM, the optical depth is lower by about an order of magnitude for LambdaCDM, and by about two orders of magnitude for S/tauCDM. These differences originate from the different epochs of cluster formation across the cosmological models, and from the non-linearity of the strong lensing effect. We conclude that only the OCDM model can reproduce the observed arc abundance well, while the other models fail to do so by orders of magnitude.

The formation of galaxy clusters in hierarchically clustering universes is investigated by means of high resolution N-body simulations. The simulations are performed using a newly developed multi-mass scheme which combines a PM code with a high resolution N-body code. Numerical effects due to time stepping and gravitational softening are investigated as well as the influence of the simulation box size and of the assumed boundary conditions. Special emphasis is laid on the formation process and the influence of various cosmological parameters. Cosmogonies with massive neutrinos are also considered. Differences between clusters in the same cosmological model seem to dominate over differences due differing background cosmogony. The cosmological model can alter the time evolution of cluster collapse, but the merging pattern remains fairly similar, e.g. number of mergers and mass ratio of mergers. The gross properties of a halo, such as its size and total angular momentum, also evolve in a similar manner for all cosmogonies and can be described using analytical models. It is shown that the density distribution of a halo shows a characteristic radial dependence which follows a power law with a slope of $\alpha=-1$ at small and $\alpha=-3$ at large radii, independent of the background cosmogony or the considered redshift. The shape of the density profiles follows the generic form proposed by Navarro et al. (1996) for all hierarchically clustering scenarios and retains very little information about the formation process or the cosmological model. Only the central matter concentration of a halo is correlated to the formation time and therefore to the corresponding cosmogony. We emphasise the role of non-radial motions of the halo particles in the evolution of the density profile.