CMB Emulator
We also build an emulator for the CMB power spectra (TT, TE and EE), up to \(\ell=2500\). The priors of the cosmological parameters are roughly centred on the recent Planck 2018 results and are as follows:
| \(\boldsymbol{\theta}\) | Distribution | Minimum | Scale | Fiducial |
|---|---|---|---|---|
| \(\sigma_8\) | Uniform | 0.7 | 0.2 | 0.8 |
| \(\Omega_{\text{cdm}}\) | Uniform | 0.20 | 0.15 | 0.2 |
| \(\Omega_b\) | Uniform | 0.04 | 0.02 | 0.04 |
| \(h\) | Uniform | 0.62 | 0.12 | 0.7 |
| \(n_s\) | Uniform | 0.90 | 0.2 | 1.0 |
A minimal example of how we can use the CMB emulator is as follows:
from jax_cosmo.emulator import EMUCMBdata, prediction_cmb_cls
from cmbrun.cmbcls import get_config
cfg = get_config('planck')
emudata = EMUCMBdata(cfg)
cosmology = np.array([0.8120, 0.265, 0.04938, 0.6732, 0.96605])
# using the emulator to predict the power spectra
mean_tt_jax = prediction_cmb_cls(cosmology, emudata, 'tt')
mean_te_jax = prediction_cmb_cls(cosmology, emudata, 'te')
mean_ee_jax = prediction_cmb_cls(cosmology, emudata, 'ee')