Table of contents


Back to the main page 1. Compute Background Evolution 2. Compute Recombination History 3. Compute Perturbations 4. Compute Powerspectrum


Compute Perturbations


Here you can select the cosmologically parameters and numerically compute the evolution of the perturbations in $\Lambda$CDM and make some plots. We assume massless neutrinos. NB: The server this runs on have old CPUs and old compilers so the runtime is unfortunately slow. The precision settings are therefore not for perfect (subpercent) accuracy, but fiducial settings are still fairly accurate and runs in ~5sec. Accuracy can be improved, but if settings are improved too much the code might not finish in time (I automatically kill any runs that takes more than a minute or so) so just be aware of that. The results using the fiducial parameters below can also be found here. Important assumptions: 1) for evolution of perturbations and power-spectra the equations implemented assumes $\Omega_k = 0$ (i.e. the curvature is just implemented in the background), 2) neutrinos are always assumed to be massless, 3) for recombination we only use the Saha approximation for Helium and the Peebles for hydrogen.

Background parameters

Hubble parameter ($H_0\equiv 100 h{\rm km/s/ Mpc}$):
$h$:
Baryon density:
$\Omega_b$:
CDM density:
$\Omega_{\rm CDM}$:
Curvature density ($\Omega_k \equiv -k/H_0^2$):
$\Omega_k$:
Temperature of the CMB today ($\Omega_\gamma \propto T_{\rm CMB}^4 / h^2$):
$T_{\rm CMB}$ (K):
Effective number of neutrinos ($\Omega_\nu \equiv N_{\rm eff}\frac{7}{8}\left(\frac{4}{11}\right)^{4/3}\Omega_\gamma$):
$N_{\rm eff}$:
Dark energy equation of state (CPL parametrisation) $w = w_0 + w_a(1-a)$:
$w_0$:
$w_a$:

Recombination parameters

Primordial helium (mass) abundance:
$Y_p$:
Include Reionization:
Reionization (H and He) redshift:
$z_{\rm reion}$:
Reionization (H and He) redshift width:
$\Delta z_{\rm reion}$:
Include He+ Reionization:
He+ reionization redshift:
$z_{\rm HeReion}$:
He+ reionization redshift width:
$\Delta z_{\rm HeReion}$:

Perturbations parameters

Include polarization:
Include neutrino perturbations:
(Precision) Number of photon multipoles ($3-20$):
$n_\ell$:
(Precision) Number of neutrino multipoles ($3-20$):
$n_\ell$:
(Precision) Minimum wavenumber:
$k_{\rm min}\eta_0$:
(Precision) Maximum wavenumber:
$k_{\rm max}\eta_0$:
(Precision) Number of $k$-values per log-interval ($10-200$):
$n$:
(Precision) Spacing between points to store while solving:
$\Delta x$:
Comments (if something is not working / suggestions) max 200 characters:

(can take a few seconds to load, just be patient!)