We have simulated the formation of an X-ray cluster in a cold dark matter universe using 12 different codes. The codes span the range of numerical techniques and implementations currently in use, including SPH and grid methods with fixed, deformable or multilevel meshes. The goal of this comparison is to assess the reliability of cosmological gas dynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at different epochs, global properties such as mass, temperature and X-ray luminosity, and radial profiles of various dynamical and thermodynamical quantities. On the whole, the agreement among the various simulations is gratifying although a number of discrepancies exist. Agreement is best for properties of the dark matter and worst for the total X-ray luminosity. Even in this case, simulations that adequately resolve the core radius of the gas distribution predict total X-ray luminosities that agree to within a factor of two. Other quantities are reproduced to much higher accuracy. For example, the temperature and gas mass fraction within the virial radius agree to about 10%, and the ratio of specific kinetic to thermal energies of the gas agree to about 5%. Various factors contribute to the spread in calculated cluster properties, including differences in the internal timing of the simulations. Based on the overall consistency of results, we discuss a number of general properties of the cluster we have modelled.

High-resolution N-body simulations are used to investigate systematic trends in the mass profiles and total masses of clusters as derived from 3 simple estimators: (1) the weak gravitational lensing shear field under the assumption of an isothermal cluster potential, (2) the dynamical mass obtained from the measured velocity dispersion under the assumption of an isothermal cluster potential, and (3) the classical virial estimator. The clusters consist of order 2.5e+05 particles of mass $m_p \simeq 10^{10} \Msun$, have triaxial mass distributions, and significant substructure exists within their virial radii. Not surprisingly, the level of agreement between the mass profiles obtained from the various estimators and the actual mass profiles is found to be scale-dependent. The virial estimator yields a good measurement of the total cluster mass, though it is systematically underestimated by of order 10%. This result suggests that, at least in the limit of ideal data, the virial estimator is quite robust to deviations from pure spherical symmetry and the presence of substructure. The dynamical mass estimate based upon a measurement of the cluster velocity dispersion and an assumption of an isothermal potential yields a poor measurement of the total mass. The weak lensing estimate yields a very good measurement of the total mass, provided the mean shear used to determine the equivalent cluster velocity dispersion is computed from an average of the lensing signal over the entire cluster (i.e. the mean shear is computed interior to the virial radius). [abridged]

The growth of structure from scale-free initial conditions is one of the most important tests of cosmological simulation methods, providing a realistically complex problem in which numerical results can be compared to rigorous analytic scaling laws. Previous studies of this problem have incorporated gravitational dynamics and adiabatic gas dynamics, but radiative cooling, an essential element of the physics of galaxy formation, normally violates the conditions necessary for self-similar evolution. We show that for any specified value of the initial power spectrum index n [where P(k)~k^n], there is a family of power-law cooling functions that preserves self-similarity by ensuring that the cooling time in an object of the characteristic mass M_* is a fixed fraction t_C of the Hubble time. We perform hydrodynamic numerical simulations with an Einstein-de Sitter cosmology, a baryon fraction of 5%, Gaussian initial conditions, two different power spectrum indices, and four values of t_C for each index, ranging from no cooling to strong cooling. We restrict the numerical simulations to two dimensions in order to allow exploration of a wide parameter space with adequate dynamic range. In all cases, the simulations are remarkably successful at reproducing the analytically predicted scalings of the mass function of dissipated objects and the gas temperature distributions and cooled gas fractions in collapsed systems. While similar success with 3-D simulations must still be demonstrated, our results have encouraging implications for numerical studies of galaxy formation, indicating that simulations with resolution comparable to that in many current studies can accurately follow the collapse and dissipation of baryons into the dense, cold systems where star formation is likely to occur.

This paper presents an alternative formulation of the ASPH algorithm for evolving anisotropic smoothing kernels, in which the geometric approach of Shapiro et al. (1996; Paper I) is replaced by an approach involving a local transformation of coordinates to those in which the underlying anisotropic volume changes appear to be isotropic. The ASPH method is presented in 2D and 3D, including a number of details not previously included in Paper I, some of which represent either advances or different choices with respect to Paper I. Among the advances included here are an asynchronous time-integration scheme with different time steps for different particles and the generalization of the ASPH method to 3D. The shock-tracking algorithm described in Paper I for locally adapting the artificial viscosity to restrict viscous heating just to particles encountering shocks, is not included here. Instead, we adopt a different interpolation kernel for use with the artificial viscosity, which has the effect of spatially localizing effects of the artificial viscosity. This version of the ASPH method in 2D and 3D is then applied to a series of 1D, 2D, and 3D test problems, and the results are compared to those of standard SPH applied to the same problems. These include the problem of cosmological pancake collapse, the Riemann shock tube, cylindrical and spherical Sedov blast waves, the collision of two strong shocks, and problems involving shearing disks intended to test the angular momentum conservation properties of the method. These results further support the idea that ASPH has significantly better resolving power than standard SPH for a wide range of problems, including that of cosmological structure formation. (Abridged)

The angular correlation function \wth of faint galaxies is affected both by nonlinear gravitational evolution and by magnification bias due to gravitational lensing. We compute the resulting \wth for different cosmological models and show how its shape and redshift evolution depend on \Omega and \Lambda. For galaxies at redshift greater than 1 (R magnitude fainter than about 24), magnification bias can significantly enhance or suppress \wth, depending on the slope of the number-magnitude relation. We show for example how it changes the ratio of \wth for two galaxy samples with different number-count slopes.

The formation history of rich clusters is investigated using a hybrid N-body simulation in which high spatial and mass resolution can be achieved self-consistently within a small region of a very large volume. The evolution of three massive clusters is studied via mass accretion, spherically-averaged density profiles, three-dimensional and projected shapes, and degree of substructure. Each cluster consists of order $4\times 10^5$ particles at the present epoch and in the case that rich cluster evolution is well-described by a 1-parameter family, the simulations have sufficient resolution to demonstrate this. At $z=0$ the clusters have similar masses, $M(r \le 1.5h^{-1} Mpc) \sim 2\times 10^{15} h^{-1} M_\odot$, and similar spherically-averaged density profiles, however markedly different formation histories are observed. No single, dominant pattern is apparent in the time variation of the mass accretion rate, the cluster shape, or the degree of substructure. Although not a statistically large sample, these objects suggest that the detailed formation history of rich clusters cannot be characterized by a simple 1-parameter family. These results suggest that the use of observations of rich clusters over a wide range of redshifts to constrain cosmological parameters may not be entirely straightforward.

Weak gravitational lensing from large-scale structure enhances and reduces the fluxes from extragalactic point sources with an rms amplitude of order 15%. In cosmic microwave background (CMB) experiments, sources exceeding some flux threshold are removed, which means that lensing will modulate the brightness map of the remaining unresolved sources. Since this mean brightness is of order 100 micro-Kelvin at 30 GHz for a reasonable flux cut, one might be concerned that this modulation could cause substantial problems for future CMB experiments. We present a detailed calculation of this effect and, fortunately, find that its power spectrum is always smaller than the normal point source power spectrum. Thus although this effect should be taken into account when analysing future high-precision CMB measurements, it will not substantially reduce the accuracy with which cosmological parameters can be measured.

We investigate the properties of hybrid gravitational/hydrodynamical simulations, examining both the numerics and the general physical properties of gravitationally driven, hierarchical collapse in a mixed baryonic/dark matter fluid. We demonstrate that, under certain restrictions, such simulations converge with increasing resolution to a consistent solution. The dark matter achieves convergence provided that the relevant scales dominating nonlinear collapse are resolved. If the gas has a minimum temperature (as expected when intergalactic gas is heated by photoionization due to the ultraviolet background) and the corresponding Jeans mass is resolved, then the baryons also converge. However, if there is no minimum baryonic collapse mass or if this scale is not resolved, then the baryon results err in a systematic fashion. In such a case, as resolution is increased the baryon distribution tends toward a higher density, more tightly bound state. We attribute this to the fact that under hierarchical structure formation on all scales there is always an earlier generation of smaller scale collapses, causing shocks which irreversibly alter the state of the baryon gas. In a simulation with finite resolution we miss such earlier generation collapses, unless a physical scale is introduced below which structure formation is suppressed in the baryons. We also find that the baryon/dark matter ratio follows a characteristic pattern, such that collapsed structures possess a baryon enriched core (enriched by factors of 2 or more over the universal average) which is embedded within a dark matter halo, even without accounting for radiative cooling of the gas. The dark matter is unaffected by changing the baryon distribution (at least in the dark matter dominated case investigated here).

We compute the two-point angular correlation function $w(\theta)$ for a sample of $\sim$ 1700 galaxies to a magnitude-limit equivalent to $R \sim 29.5$ using a catalog derived from the Hubble Deep Field images. A non zero value of $w(\theta)$ is measured down to $R=29.0$. The amplitude of $w(\theta)$ at the bright magnitude limit ($R \sim 26$) is consistent with previous ground-based observations. At fainter magnitudes the clustering amplitude continues to decrease but at a slower rate than that predicted by the power law $w(1'')\propto 10^{-0.27R}$ observed for shallower samples. The observed $w(\theta)$ over the magnitude range $20 <R< 29$ is consistent with linear evolution of the clustering of a galaxy population which at present has a correlation length $r_0$ of about 4 h^(-1) Mpc, close to that of local IRAS galaxies. We also investigate the impact that magnification bias induced by weak gravitational lensing may have on our results. Although the observed amplitude of $w(\theta)$ can differ from the true amplitude by up to 30%, this effect is not large enough to affect our conclusions. Finally, by using a color-selected sample, we examine whether the expected effects of magnification bias can be used for an independent determination of cosmological parameters in deep images. We conclude that the amplitude of the effect can be large and in some cases even produce an upturn of the amplitude of the correlation with limiting magnitude. However, we find that it is not possible to detect the effects of magnification bias on $w(\theta)$ from images alone. If redshift information becomes available, it is possible to measure the effects of magnification bias directly and thus constrain the density parameter $\Omega_0$ and the bias factor b.

Weak gravitational lensing by large scale structure affects the number counts of faint galaxies through the ``magnification bias'' and thus affects the measurement of the angular two-point correlation function $\w $. At faint magnitudes the clustering amplitude will decrease differently with limiting magnitude than expected from Limber's equation. The amplitude will hit a minimum and then rise with limiting magnitude. This behavior occurs because $\w$ due to clustering decreases with distance, while the ``magnification bias'' due to weak lensing increases with distance. The apparent magnitude $m_{min}$ at which the magnification bias starts to dominate the observed clustering is model and color dependent. It is given by $\omega(m=m_{min},\theta=5^\prime) \approx (1\ -\ 2)\times 10^{-3}(5s-2)^2 \Omega_0^2 \sigma_8^2$, where $s$ is the logarithmic slope of the number counts. Already published measurements of $\w$ at $R=25$ may be strongly influenced by the ``magnification bias''. An experiment using the ratio of blue and red number counts across the sky can be designed such that the effects of the ``true'' clustering is minimized. The magnification bias is a measurement of the clustering of the mass. This weak lensing experiment does not require measuring shapes and position angles of galaxies. I derive a revised Limber's Equation including the effects of magnification bias.

A single very deep $10'$ field has been imaged to $r \sim 26$. There are 2682 galaxies with $23 \le r \le 26$ with a size significantly larger than the seeing disk in the field. After correcting for telescope aberrations, possible guiding errors, and signal degradation due to seeing, I find a polarization signal of amplitude $|p|=2.7\%$, with an uncertainty of 1.2\% (95\% confidence limit). For the 1773 galaxies with $23 \le r \le 25$, the measured amplitude is 2.4\% $\pm$ 1.2\% (95\% confidence limit). If this observed polarization is due to gravitational lensing by large scale structure, cosmological inferences can be made. The polarization amplitude would indicate that $\Omega_0 \times \sigma_8 \approx 1$. The amplitude probability distribution is a Rayleigh distribution with variance $\sigma_p^2$ so the probability $P(1/2 < |p|/\sigma_p <2) = 75\%$, and $P(1/3 < |p|/\sigma_p < 3) = 93.5\%$. The measurement of the polarization amplitude in a single field can be used to constrain cosmological models since $\sigma_p \propto \Omega_0 \times \sigma_8$.

Weak lensing is the distortion (polarization) of images of distant objects, such as high redshift galaxies, by gravitational fields in the limit where the distortion is small. Gravitational potential fluctuations due to large scale structure cause correlated distortions of the images of high redshift galaxies. These distortions are observable with current large telescopes and instrumentation. In a Friedmann-Robertson-Walker (FRW) metric I calculate the weak lensing pattern due to large scale structure for arbitrary $\Omega_0$ and zero cosmological constant $\Lambda$. For a given cosmological model, specified by $\Omega_0$ and a power spectrum of density fluctuations, I calculate the statistical properties of the polarization field for an arbitrary redshift source distribution in a simple closed form. It is shown that for low redshift $z$ of the sources, the polarization amplitude is proportional to $\Omega_0$, while at higher redshift the polarization measures the value of $\Omega(z)$, where $z$ is the characteristic source redshift. Moreover, the statistics of the polarization field are a direct measure of the power spectrum of density fluctuations.