Hierarchical clustering and the baryon distribution in galaxy clusters

The baryon fraction of galaxy clusters in numerical simulations is found to be dependant on the cluster formation method. In all cases, the gas is anti-biased compared with the dark matter. However, clusters formed hierarchically are found to be more depleted in baryons than clusters formed non-hierarchically. There is a depletion of 10 to 15% for hierarchically formed clusters while the depletion is less than 10% for those clusters formed non-hierarchically. This difference is dependent on the mass of the clusters. The mean baryon enrichment profile for the hierarchically formed clusters shows an appreciable baryon enhancement around the virial radius not seen in the clusters formed without substructure. If this phenomenon applies to real clusters, it implies that determinations of the baryon fractions in clusters of galaxies require data extending beyond the virial radius of the clusters in order to achieve an unbiased value.Comment: 13 pages including 2 tables and 2 figures. Submitted to MNRA

Hydra: An Adaptive--Mesh Implementation of PPPM--SPH

We present an implementation of Smoothed Particle Hydrodynamics (SPH) in an adaptive-mesh PPPM algorithm. The code evolves a mixture of purely gravitational particles and gas particles. The code retains the desirable properties of previous PPPM--SPH implementations; speed under light clustering, naturally periodic boundary conditions and accurate pairwise forces. Under heavy clustering the cycle time of the new code is only 2--3 times slower than for a uniform particle distribution, overcoming the principal disadvantage of previous implementations\dash a dramatic loss of efficiency as clustering develops. A 1000 step simulation with 65,536 particles (half dark, half gas) runs in one day on a Sun Sparc10 workstation. The choice of time integration scheme is investigated in detail. A simple single-step Predictor--Corrector type integrator is most efficient. A method for generating an initial distribution of particles by allowing a a uniform temperature gas of SPH particles to relax within a periodic box is presented. The average SPH density that results varies by $\sim\pm1.3$\%. We present a modified form of the Layzer--Irvine equation which includes the thermal contribution of the gas together with radiative cooling. Tests of sound waves, shocks, spherical infall and collapse are presented. Appropriate timestep constraints sufficient to ensure both energy and entropy conservation are discussed. A cluster simulation, repeating Thomas andComment: 29 pp, uuencoded Postscrip

Smoothed Particle Hydrodynamics in cosmology: a comparative study of implementations

We analyse the performance of twelve different implementations of Smoothed Particle Hydrodynamics (SPH) using seven tests designed to isolate key hydrodynamic elements of cosmological simulations which are known to cause the SPH algorithm problems. In order, we consider a shock tube, spherical adiabatic collapse, cooling flow model, drag, a cosmological simulation, rotating cloud-collapse and disc stability. In the implementations special attention is given to the way in which force symmetry is enforced in the equations of motion. We study in detail how the hydrodynamics are affected by different implementations of the artificial viscosity including those with a shear-correction modification. We present an improved first-order smoothing-length update algorithm that is designed to remove instabilities that are present in the Hernquist and Katz (1989) algorithm. For all tests we find that the artificial viscosity is the most important factor distinguishing the results from the various implementations. The second most important factor is the way force symmetry is achieved in the equation of motion. Most results favour a kernel symmetrization approach. The exact method by which SPH pressure forces are included has comparatively little effect on the results. Combining the equation of motion presented in Thomas and Couchman (1992) with a modification of the Monaghan and Gingold (1983) artificial viscosity leads to an SPH scheme that is both fast and reliable.Comment: 30 pages, 26 figures and 9 tables included. Submitted to MNRAS. Postscript version available at ftp://phobos.astro.uwo.ca/pub/etittley/papers/sphtest.ps.g

The effect of radiative cooling on scaling laws of X-ray groups and clusters

We have performed cosmological simulations in a ΛCDM cosmology with and without radiative cooling in order to study the effect of cooling on the cluster scaling laws. Our simulations consist of 4.1 million particles each of gas and dark matter within a box size of 100 h-1 Mpc, and the run with cooling is the largest of its kind to have been evolved to z = 0. Our cluster catalogs both consist of over 400 objects and are complete in mass down to ~1013 h-1 M☉. We contrast the emission-weighted temperature-mass (Tew-M) and bolometric luminosity-temperature (Lbol-Tew) relations for the simulations at z = 0. We find that radiative cooling increases the temperature of intracluster gas and decreases its total luminosity, in agreement with the results of Pearce et al. Furthermore, the temperature dependence of these effects flattens the slope of the Tew-M relation and steepens the slope of the Lbol-Tew relation. Inclusion of radiative cooling in the simulations is sufficient to reproduce the observed X-ray scaling relations without requiring excessive nongravitational energy injection