We use a high resolution collisionless simulation of a Virgo--like cluster in a $\Lambda$CDM cosmology to determine the velocity and clustering properties of the diffuse stellar component in the intracluster region at the present epoch. The simulated cluster builds up hierarchically and tidal interactions between member galaxies and the cluster potential produce a diffuse stellar component free-flying in the intracluster medium. Here we adopt an empirical scheme to identify tracers of the stellar component in the simulation and hence study its properties. We find that at $z=0$ the intracluster stellar light is mostly unrelaxed in velocity space and clustered in structures whose typical clustering radii are about 50 kpc at R=400--500 kpc from the cluster center, and predict the radial velocity distribution expected in spectroscopic follow-up surveys. Finally, we compare the spatial clustering in the simulation with the properties of the Virgo intracluster stellar population, as traced by ongoing intracluster planetary nebulae surveys in Virgo. The preliminary results indicate a substantial agreement with the observed clustering properties of the diffuse stellar population in Virgo.

We examine how the statistics of the quadrupoles of (projected) cluster masses can discriminate between flat cold dark matter (CDM) universes with or without a cosmological constant term. Even in the era of high precision cosmology that cosmic microwave background experiments should open soon, it is important to devise self consistency tests of cosmogonic theories tuned at the matter radiation decoupling epoch using data from the non--linear evolved universe. We build cluster catalogs from two large volume simulations of a ``tilted'' CDM model and a $\Lambda$CDM model with cosmic density parameter $\Omega_m=0.35$ and cosmological constant contribution $\Omega_\Lambda=0.65$. From the projected mass distribution of the clusters we work out the quadrupoles Q and examine their dependence on cluster mass and the cosmological model. We find that TCDM clusters have systematically larger quadrupoles than their \lcdm counterpart. The effect is mass dependent: massive clusters ($M\gsim 10^{15}\Msunh$) have quadrupoles differing by more than 30% in the two models, while for $M\lsim 4\times10^{14}\Msunh$ the difference rapidly drops to $\sim 1$%. Performing a K-S test of the Q distributions, we estimate that using just the 15 most massive clusters in the simulation volume ($360\Mpch$ a side) we can discriminate between TCDM and \lcdm at a confidence level better than 99.9%. In the volume probed by exhisting observations, there are potentially several hundred clusters with masses above the threshold for which the differences in the quadrupoles become relevant. Should weak lensing data become available for this whole set, a quadrupole analysis may be expected to discriminate among different values of $\Lambda$.

We use a cosmological simulation of the Local Group to make quantitative and speculative predictions for direct detection experiments. Cold dark matter (CDM) halos form via a complex series of mergers, accretion events and violent relaxation which precludes the formation of significant caustic features predicted by axially symmetric collapse. The halo density profiles are combined with observational constraints on the galactic mass distribution to constrain the local density of cold dark matter to lie in the range 0.18 <~ rho_CDM(R_solar)/GeV cm^-3 <~ 0.30. In velocity space, coherent streams of dark matter from tidally disrupted halos fill the halo and provide a tracer of the merging hierarchy. The particle velocities within triaxial CDM halos cannot be approximated by a simple Maxwellian distribution and is radially biased at the solar position. The detailed phase space structure within the solar system will depend on the early merger history of the progenitor halos and the importance of major mergers over accretion dominated growth. We follow the formation of a ``Draco'' sized dSph halo of mass 10^8M_solar with several million particles and high force accuracy. Its internal structure and substructure resembles that of galactic or cluster mass halos: the density profile has a singular central cusp and it contains thousands of sub-halos orbiting within its virial radius demonstrating a self-similar nature to collisionless dark matter sub-clustering. The singular cores of substructure halos always survive complete tidal disruption although mass loss is continuous and rapid. Extrapolating wildly to earth mass halos with velocity dispersion of 1 m s^-1 (roughly equal to the free streaming scale for neutralinos) we find that most of the dark matter may remain attached to bound subhalos. (Abridged)

For the last few years the Cold Dark Matter model (ticker: CDM), has been the dominant theory of structure formation. We briefly review the recent advancements and predictions of the model in the field of galaxy clusters. A new set of very high resolution simulations of galaxy clusters show that they have (1) density profiles with central slopes very close to -1.6 and (2) abundance of subhalos similar to the ones observed in real clusters. These results show a remarkably small cluster to cluster variation and a weak dependence from the particular CDM cosmology chosen (LCDM having ~40% less substructure than SCDM). While still a speculative theory with a high prediction/evidence ratio, subject to strong challenges from observational data and competition from other hierarchical theories, we give CDM a rating of ``market outperform'' and of ``long term BUY''.

Can N-body simulations reliably determine the structural properties of dark matter halos? Focussing on a Virgo-sized galaxy cluster, we increase the resolution of current ``high resolution simulations'' by almost an order of magnitude to examine the convergence of the important physical quantities. We have 4 million particles within the cluster and force resolution 0.5 kpc/h (0.05% of the virial radius). The central density profile has a logarithmic slope of -1.5, as found in lower resolution studies of the same halo, indicating that the profile has converged to the ``physical'' limit down to scales of a few kpc. Also the abundance of substructure is consistent with that derived from lower resolution runs; on the scales explored, the mass and circular velocity functions are close to power laws of exponents ~ -1.9 and -4. Overmerging appears to be globally unimportant for suhalos with circular velocities > 100 km/s. We can trace most of the cluster progenitors from z=3 to the present; the central object (the dark matter analog of a cD galaxy)is assembled between z=3 and 1 from the merging of a dozen halos with v_circ \sim 300 km/s. The mean circular velocity of the subhalos decreases by ~ 20% over 5 billion years, due to tidal mass loss. The velocity dispersions of halos and dark matter globally agree within 10%, but the halos are spatially anti-biased, and, in the very central region of the cluster, they show positive velocity bias; however, this effect appears to depend on numerical resolution.

We combine semi-analytical methods with a ultra-high resolution simulation of a galaxy cluster (of mass 2.3 10^14h-1Msolar, and 4 10^6 particles within its virial radius) formed in a standard CDM universe to study the spatial distribution and orbital properties of the present-day descendents of Lyman Break Galaxies (LBGs). At the present time only five (out of 12) of halos containing LBGs survive as separate entities inside the cluster virial radius. Their circular velocities are in the range 200 - 550 km/sec. Seven halos merged together to form the central object at the very center of the cluster. Using semi-analytical modeling of galaxy evolution we show that descendents of halos containing LBGs now host giant elliptical galaxies. Galaxy orbits are radial, with a pericenter to apocenter ratio of about 1:5. The orbital eccentricities of LBGs descendents are statistically indistinguishable from those of the average galaxy population inside the cluster, suggesting that the orbits of these galaxies are not significantly affected by dynamical friction decay after the formation of the cluster's main body. In this cluster, possibly due to its early formation time, the descendents of LBGs are contained within the central 60% of the cluster virial radius and have an orbital velocity dispersion lower than the global galaxy population, originating a mild luminosity segregation for the brightest cluster members. Mass estimates based only on LBGs descendents (especially including the central cD) reflect this bias in space and velocity and underestimate the total mass of this well virialized cluster by up to a factor of two compared to estimates using at least 20 cluster members.

We use numerical simulations to examine the substructure within galactic and cluster mass halos that form within a hierarchical universe. Clusters are easily reproduced with a steep mass spectrum of thousands of substructure clumps that closely matches observations. However, the survival of dark matter substructure also occurs on galactic scales, leading to the remarkable result that galaxy halos appear as scaled versions of galaxy clusters. The model predicts that the virialised extent of the Milky Way's halo should contain about 500 satellites with circular velocities larger than Draco and Ursa-Minor i.e. bound masses > 10^8Mo and tidally limited sizes > kpc. The substructure clumps are on orbits that take a large fraction of them through the stellar disk leading to significant resonant and impulsive heating. Their abundance and singular density profiles has important implications for the existence of old thin disks, cold stellar streams, gravitational lensing and indirect/direct detection experiments.

We examine the properties of dark matter halos within a rich galaxy cluster using a high resolution simulation that captures the cosmological context of a cold dark matter universe. The mass and force resolution permit the resolution of 150 halos with circular velocities larger than 80 kms within the cluster's virial radius of 2 Mpc. This enables an unprecedented study of the statistical properties of a large sample of dark matter halos evolving in a dense environment. The cumulative fraction of mass attached to these halos varies from 0% at 200 kpc, to 13% at the virial radius. Even at this resolution the overmerging problem persists; halos that pass within 200 kpc of the cluster center are tidally disrupted. Additional substructure is lost at earlier epochs within the massive progenitor halos. The median ratio of apocentric to pericentric radii is 6:1; the orbital distribution is close to isotropic, circular orbits are rare, radial orbits are common. The orbits of halos are unbiased with respect to both position within the cluster and with the orbits of the smooth dark matter background and no velocity bias is detected. The tidal radii of surviving halos are generally well-fit using the simple analytic prediction applied to their orbital pericenters. Halos within clusters have higher concentrations than those in the field. Within the cluster, halo density profiles can be modified by tidal forces and individual encounters with other halos that cause significant mass loss - ``galaxy harassment''. Mergers between halos do not occur inside the clusters virial radius.

We investigate the internal structure of cold dark matter halos using high resolution N-body simulations. As the mass and force resolution are increased, halo density profiles become steeper, asymptoting to a slope of $\sim -4/3$ in the central regions and may not have converged to a unique result. At our highest resolution we have nearly 3 million particles within the virial radius, $R_{200}$, and force softening that is $\sim 0.2% R_{200}$. This resolution has also allowed us to resolve a large part of the overmerging problem - we find over 1000 surviving dark halos orbiting within a single cluster potential. These data have given us unprecedented insights into the dynamics and structure of ``halos within halos'', allowing us for the first time, to compare the distribution of dark matter with observations of galaxies in clusters.

We apply several statistical estimators to high-resolution N-body simulations of two currently viable cosmological models: a mixed dark matter model, having $\Omega_\nu=0.2$ contributed by two massive neutrinos (C+2\nuDM), and a Cold Dark Matter model with Cosmological Constant (\LambdaCDM) with $\Omega_0=0.3$ and h=0.7. Our aim is to compare simulated galaxy samples with the Perseus-Pisces redshift survey (PPS). We consider the n-point correlation functions (n=2-4), the N-count probability functions P_N, including the void probability function P_0, and the underdensity probability function U_\epsilon (where \epsilon fixes the underdensity threshold in percentage of the average). We find that P_0 (for which PPS and CfA2 data agree) and P_1 distinguish efficiently between the models, while U_\epsilon is only marginally discriminatory. On the contrary, the reduced skewness and kurtosis are, respectively, S_3\simeq 2.2 and S_4\simeq 6-7 in all cases, quite independent of the scale, in agreement with hierarchical scaling predictions and estimates based on redshift surveys. Among our results, we emphasize the remarkable agreement between PPS data and C+2\nuDM in all the tests performed. In contrast, the above \LambdaCDM model has serious difficulties in reproducing observational data if galaxies and matter overdensities are related in a simple way.

We use the void probability function to compare the redshift--space galaxy distribution in the Perseus--Pisces survey with artificial samples from N-body simulations of standard cold dark matter (CDM) and broken scale invariance (BSI) models. Galaxies are identified as residing in peaks of the evolved density field in such a way as to reproduce both the observed luminosity and two--point correlation functions. Using a similar approach, it was recently shown that the VPF can discriminate between CDM and a cold+hot dark matter (CHDM) model with cold/hot/baryon fractions = 0.6/0.3/0.1. Our main result is that both CDM (as expected from a previous analysis) and BSI fit observational data. The robustness of the result is checked against changing the observer's position in the simulations and the galaxy identification in the evolved density field. Therefore, while the void statistics is sensitive to the passage from CDM to CHDM (different spectrum and different nature of dark matter), it is not to the passage from CDM to BSI (different spectrum but same dark matter). On such a basis, we conjecture that the distribution of voids could be directly sensitive to the nature of dark matter, but scarcely sensitive to the shape of the transfer function.

We use the void probability statistics to study the redshift-space galaxy distribution as described by a volume-limited subsample of the Perseus-Pisces survey. We compare the results with the same analysis realized on artificial samples, extracted from high-resolution N-body simulations by reproducing the observational biases of the real data set. Simulations are run for the Cold+HotDM model (CHDM) and for unbiased and biased (b=1.5) CDM models in a 50 Mpc/h box. We identify galaxies as residing in peaks of the evolved density field. We fragment overmerged structures into individual galaxies so as to reproduce both the correct luminosity function (after assuming M/ L values for the resulting galaxy groups) and the two-point correlation function. Our main result is that a void-probability function (VPF) from the standard CHDM model with fractions 60% cold, 30% hot, 10% barions, exceeds the observational VPF with a high confidence level. CDM models produce smaller VPF independent of the biasing parameter. We verify the robustness of this result against changing the observer position in the simulations and the galaxy identification in the evolved density field.

We discuss the effects of redshift--space distortions on the estimate of high order galaxy correlation functions. We use both the Perseus--Pisces redshift survey and the results of high--resolution N--body simulations to explore the consequences of working in redshift space on the detection of deviations from hierarchical clustering. Both for real and simulated data, significant deviations from hierarchical clustering seem to be present in real space. Their behaviour is coherent with that expected from an initially biased galaxy field, once the displacement from the sites where galaxies formed, due to the nonlinear gravitational evolution, is taken into account. The passage to redshift space has the net effect of filtering out the higher powers required to fit the distribution in real space. This magnifies the distortions operated by nonlinear evolution on the initial distribution, erasing the residual narrow scale range (2--5~$h^{-1}$Mpc) where deviations from hierarchical clustering are still detectable in real space. We conclude that such deviations can hardly be estimated using data in redshift space.

We use high resolution PM N-body simulations to follow the development of non-linear clustering in a flat Universe, dominated by Cold + Hot Dark Matter (CHDM) with 60% of CDM, 30% of HDM and 10% of baryons; a simulation box of 100 Mpc a side ($h=0.5$) is used. We analyze two CHDM simulations with $b =1.5$ (COBE normalization). We also compare them with CDM simulations with $b=1.5$ and $b=1$. We evaluate high-order correlation functions and the void-probability-function (VPF). Correlation functions are obtained both from counts in cells and counts of neighbors. The analysis is made for DM particles and for galaxies, identified as massive halos in the evolved density field. We also check the effects of dynamical evolution and redshift space distortions. We find that clustering of DM particles exhibits deviations from the hierarchical scaling, which decrease somewhat in redshift space. Galaxies follow hierarchical scaling far more closely, with coefficients S_3=2.5 and S_4=7.5, in general agreement with observational results, just marginally affected by redshift distortions and on the choices for the initial spectra. The hierarchical scaling of galaxy clustering is confirmed by the VPF analysis and in all the cases considered a good agreement with observational results is obtained.