Synthetic galaxy power spectra

Galaxy formation model

From this page you will be able to download the auto and cross power spectra of the 8000 synthetic galaxy samples presented in Zennaro et al., 2021b. Each galaxy sample is obtained from one of our four main simulations (Angulo et al., 2021; Contreras et al., 2020), at redshifts \(z=0,1\), at fixed number densities \(\bar{n} = \{0.0003, 0.001, 0.00316, 0.01\}\) \(h^3 \mathrm{Mpc}^{-3}\), and randomly varying the parameters of the galaxy formation model. To build our galaxy catalogues we employ the SHAMe model (Contreras et al., 2021b), both for stellar-mass selected samples (3 model parameters), and star-formation rate selected samples (5 model parameters).

Download and usage

All of our spectra are available as python dictionaries stored in HDF5 format. You can download them through these links:

sm_galaxy_samples.hdf5 21 M
sfr_galaxy_samples.hdf5 21 M

You can also directly download these files from the command line with 

wget https://bacco.dipc.org/SHAMeGalPk/sm_galaxy_samples.hdf5
wget https://bacco.dipc.org/SHAMeGalPk/sfr_galaxy_samples.hdf5

Once download, you can use these files loading them with your favorite HDF5 reading library and accessing the dictionary contained within. For example 

import deepdish as dd
sfr = dd.io.load('sfr_galaxy_samples.hdf5')
sm = dd.io.load('sm_galaxy_samples.hdf5')
print(sfr.keys())
will give 
dict_keys(['SHAMe_pars', 'atarget', 'auto', 'cosmology', 'cross', 'ndens', 'satellite_fraction'])

These keywords contain

  • SHAMe_pars: the choice of SHAMe parameters for each sample (\(\sigma_{M_*}, t_{\rm merger}, f_{\rm s}\) for SM-selected galaxies, \(\beta, \gamma, \log_{10}(M_1), \tau_0, \tau_{\rm s}\) for SFR-selected galaxies).
  • atarget: the expansion factor of the sample.
  • cosmology: the name of the cosmology (nenya, narya, vilya, power).
  • ndens: the number density of the sample in \(h^3 \mathrm{Mpc}^{-3}\).
  • satellite_fraction: the satellite fraction of the sample.
  • auto: a dictionary with wavenumbers k, galaxy auto power spectrum pk, number of modes in each wavenumber bin nmodes, statistical error assuming the simulation volume and accounting for fixing and pairing fnp_pk_error, shotnoise shotnoise, separation bins r, 2-point correlation function (from inverse FT of power spectrum) xi.
  • cross: a dictionary with wavenumbers k, galaxy matter cross power spectrum pk, number of modes in each wavenumber bin nmodes, statistical error assuming the simulation volume and accounting for fixing and pairing fnp_pk_error, shotnoise shotnoise.

If you use these spectra, please cite these two papers: