We report on a series of tests of agreement between three types of N-body simulations: PM, P$^3$M, and Tree codes. We find good agreement in both the individual and the statistical properties only on scales larger than the mean interparticle separation. As a result, we question most numerical results at and below below galaxy scales, either concerning primordial dark matter or baryonic matter coupled to it by gravitation. Mass resolution (specifically, sufficient to have the mean interparticle separation comparable to the smoothing scales) appears be necessary for results to converge.

We explore some of the effects that discreteness and two-body scattering may have on N-body simulations with ``realistic'' cosmological initial conditions. We use an identical subset of particles from the initial conditions for a $128^3$ Particle-Mesh (PM) calculation as the initial conditions for a variety P$^3$M and Tree code runs. We investigate the effect of mass resolution (the mean interparticle separation) since most ``high resolution'' codes only have high resolution in gravitational force. The phase-insensitive two--point statistics, such as the power spectrum (autocorrelation) are somewhat affected by these variations, but phase-sensitive statistics show greater differences. Results converge at the mean interparticle separation scale of the lowest mass-resolution code. As more particles are added, but the force resolution is held constant, the P$^3$M and the Tree runs agree more and more strongly with each other and with the PM run which had the same initial conditions. This shows high particle density is necessary for correct time evolution, since many different results cannot all be correct. However, they do not so converge to a PM run which continued the fluctuations to small scales. Our results show that ignoring them is a major source of error on comoving scales of the missing wavelengths. This can be resolved by putting in a high particle density. Since the codes never agree well on scales below the mean comoving interparticle separation, we find little justification for quantitative predictions on this scale. Some measures vary by 50%, but others can be off by a factor of three or more. Our results suggest possible problems with the density of galaxy halos, formation of early generation objects such as QSO absorber clouds, etc.

Two-body scattering and other discreteness effects are unimportant in cosmological gravitational clustering in most scenarios, since the dark matter has a small particle mass. The collective field should determine evolution: Two-body scattering in simulations violates the Poisson-Vlasov equations. We test this in PM, P$^3$M, Tree, and NGPM codes, noting that a collisionless code will preserve the one-dimensional character of plane wave collapse. We find collisionality vanishing as the softening parameter approaches the mean interparticle separation. Solutions for the problem are suggested, involving greater computer power, PM-based nested grid codes, and a more conservative approach to resolution claims.

In this study we use simulations of 128$^3$ particles to study the ellipticity and orientation of clusters of galaxies in N-body simulations of differing power-law initial spectra (P(k) \propto k^n ,n = +1, 0, -1, -2$), and density parameters ($\Omega_0 = 0.2$ to 1.0). Furthermore, unlike most theoretical studies we mimic most observers by removing all particles which lie at distances greater than 2 1/h Mpc from the cluster center of mass. We computed the axial ratio and the principal axes using the inertia tensor of each cluster. The mean ellipticity of clusters increases strongly with increasing $n$. We also find that clusters tend to become more spherical at smaller radii. We compared the orientation of a cluster to the orientation of neighboring clusters as a function of distance (correlation). In addition, we considered whether a cluster's major axis tends to lie along the line connecting it to a neighboring cluster, as a function of distance (alignment). Both alignments and correlations were computed in three dimensions and in projection to mimic observational surveys. Our results show that significant alignments exist for all spectra at small separations ($D < 15 h^{-1}$ Mpc) but drops off at larger distance in a strongly $n-$dependent way.

I describe a nested-grid particle-mesh (NGPM) code designed to study gravitational instability in three-dimensions. The code is based upon a standard PM code. Within the parent grid I am able to define smaller sub-grids allowing us to substantially extend the dynamical range in mass and length. I treat the fields on the parent grid as background fields and utilize a one-way interactive meshing. Waves on the coarse parent grid are allowed to enter and exit the subgrid, but waves from the subgrid are precluded from effecting the dynamics of the parent grid. On the parent grid the potential is computed using a standard multiple Fourier transform technique. On the subgrid I use a Fourier transform technique to compute the subgrid potential at high resolution. I impose quasi-isolated boundary conditions on the subgrid using the standard method for generating isolated boundary conditions, but rather than using the isolated Green function I use the Ewald method to compute a Green function on the subgrid which possesses the full periodicity of the parent grid. I present a detailed discussion of my methodology and a series of code tests.

In this work we demonstrate the ability of the Minimal Spanning Tree to duplicate the information contained within a percolation analysis for a point dataset. We show how to construct the percolation properties from the Minimal Spanning Tree, finding roughly an order of magnitude improvement in the computer time required. We apply these statistics to Particle-Mesh simulations of large-scale structure formation. We consider purely scale-free Gaussian initial conditions ($P(k) \propto k^n$, with $n = -2, -1, 0 \ \& +1$) in a critical density universe. We find in general the mass of the percolating cluster is a much better quantity by which to judge the onset of percolation than the length of the percolating cluster.

Davidsen et al. (1991) have argued that the failure to detect uv photons from the dark matter DM) in cluster A665 excludes the decaying neutrino hypothesis. Sciama et al. (1993) argued that because of high central concentration the DM in that cluster must be baryonic. We study the DM profile in clusters of galaxies simulated using the Harrison--Zel'dovich spectrum of density fluctuations, and an amplitude previously derived from numerical simulations (Melott 1984b; Anninos et al. 1991) and in agreement with microwave background fluctuations (Smoot et al. 1992). We find that with this amplitude normalization cluster neutrino DM densities are comparable to observed cluster DM values. We conclude that given this normalization, the cluster DM should be at least largely composed of neutrinos. The constraint of Davidsen et al. can be somewhat weakened by the presence of baryonic DM; but it cannot be eliminated given our assumptions.