High Resolution Microwave Full Sky Simulalations
SZ Halo Catalogs
Cosmological parameters are based closely on the five year WMAP data release.
nsinit = 0.96
hubble0 = 0.71
wde = -1.0
sigma8 = 0.80
omegam = 0.264
omegal = 0.736 = 1-omegam
omegab = 0.044
YHe = 0.24
XH = 0.76
Tcmb = 2.726
The TPM N-body code was run with 1024^3 particles in a 1000 Mpc/h box. Halos were identified using a
friends-of-friends algorithm and those belonging to a light cone spanning an octant for redshifts 0 < z < 3 have
been saved.
TPM N-body simulation
L = 1000 Mpc/h
Np = 1024^3
mp = 6.82E10 Msolar/h = 9.61E10 Msolar
FoF halo finder
bfof = 0.2 times mean interparticle spacing
Nmin = 293 particles
Mmin = 2.00E13 Msolar/h = 2.82E13 Msolar
Nhalo = 1414339
SZ maps are available as binary FITS files with the minimum header. Maps are in Healpix format with Nside = 8192
and Npix = 12*Nside^2. Full sky maps are constructed by replicating the octant as follow:
Northern hemisphere
theta_N = theta_O
phi_N = phi_O + n*pi/2 where n = 0,1,2,3
Southern hemisphere
theta_S = pi - theta_N
phi_S = 2*pi - phi_N
Sky_sz_nu.fits are distortion maps, in units of Delta T/T, at various GHz frequencies. These maps are for the
redshift range 0 < z < 10. They include relativistic corrections.
Sky_tsz.fits and Sky_ksz.fits are the Compton y and Doppler b maps, respectively. These maps are for the same
redshift range, but do not include relativistic corrections.
Sky_sz_148.fits
Sky_sz_219.fits
Sky_sz_277.fits
Sky_sz_30.fits
Sky_sz_90.fits
Sky_sz_350.fits
Sky_tsz.fits
Sky_ksz.fits
SZ halo catalogs are available in ascii and binary formats. There are 1414339 halos in the first octant for
redshifts 0 < z < 3. For each halo, there are 48 numbers as described below.
Here is a sample codelet to read in the binary data file
integer(4), parameter :: nhalo = 1414339
real(4), dimension(48,nhalo) :: halo
open(1,file='halo_sz.binary',form='binary')
do i=1,nhalo
read(1) halo(:,i)
enddo
close(1)
Halo array
halo( 1,i) = redshift
halo( 2,i) = ra in degrees
halo( 3,i) = dec in degrees
halo( 4:6,i) = comoving position of halo potential minimum in Mpc
halo( 7:9,i) = proper peculiar velocity in km/s
halo( 10,i) = Mfof in Msolar
halo( 11,i) = Mvir in Msolar
halo( 12,i) = Mgas_vir in Msolar
halo( 13,i) = Rvir in proper Mpc
halo( 14,i) = Integrated TSZ within Rvir in arcminute^2
halo( 15,i) = Integrated KSZ within Rvir in arcminute^2
halo(16:21,i) = Integrated SZ at 148,219,277,30,90,350 GHz
halo(22:32,i) = Same as 11-21 but for within R200
halo(33:43,i) = Same as 11-21 but for within R500
halo( 44,i) = Mstar in Msolar
halo( 45,i) = central gas density in Msolar/Mpc^3
halo( 46,i) = central temperature in keV
halo( 47,i) = central pressure in Msolar/Mpc/Gyr^2
halo( 48,i) = central potential in (Mpc/Gyr)^2