Sunyaev-Zeldovich effect from hydrodynamical simulations: Maps and low order statistics

We use moving mesh hydrodynamical simulations to make maps of the Sunyaev-Zeldovich effect. We present these maps for several cosmological models and explore their lowest moments. We find that the first moment, the mean Compton y parameter, is typically between $1-3\times {10}^{-6}$ for cluster abundance normalized models, the lower value corresponding to the high density models, and scales as ${\sigma }_8^{3-4}.$ Power spectra at $l\sim 1000$ have an amplitude between 1 and $3\times {10}^{-6}$ in the Rayleigh-Jeans regime, with a similar depedendence on the background cosmology. The amplitude of the power spectrum strongly depends on the power spectrum normalization and scales as ${\sigma }_8^{6-9}.$ On smaller scales $\mathrm{}(l>\mathrm{}1000\mathrm{})\mathrm{}$ the spectrum is dominated by halos below ${10}^{13}{M}_{\odot }$ and so is sensitive to the thermal history of galactic halo gas. On larger scales $\mathrm{}(l<\mathrm{}1000\mathrm{})\mathrm{}$ the power spectrum is less sensitive to nongravitational energy injection, but becomes very noisy, with large field to field variations in the power spectrum caused by rare bright sources in the maps, which dominate the spectrum over the large scale structure correlations. The cross-correlation power spectrum with weak lensing or projected galaxy distribution is significant and the cross-correlation coefficient is around 0.5 over a wide range of scales. A comparison with Press-Schechter predictions gives very good agreement for all statistics, indicating that there is no significant contribution to SZ from nonvirialized structures. The one point distribution function shows clear deviations from Gaussianity and can be well approximated as log-normal on small scales.