Table of contents
Back to the main page 1. Compute Background Evolution 2. Compute Recombination History 3. Compute Perturbations 4. Compute Powerspectrum
Compute Powerspectrum
Here you can select the cosmologically parameters and numerically compute the CMB and matter powerspectrum in $\Lambda$CDM and make some plots. We also show a comparison to Planck 2018 data. 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 ~10sec. 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 powerspectra 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(1a)$:  
$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 ($320$):  
$n_\ell$:  
(Precision) Number of neutrino multipoles ($320$):  
$n_\ell$:  
(Precision) Minimum wavenumber:  
$k_{\rm min}\eta_0$:  
(Precision) Maximum wavenumber:  
$k_{\rm max}\eta_0$:  
(Precision) Number of $k$values per loginterval ($10200$):  
$n$:  
(Precision) Spacing between points to store while solving:  
$\Delta x$:  
Powerspectrum parameters 

Primordial amplitude:  
$10^{10}A_s$:  
Spectral index:  
$n_s$:  
Pivot scale:  
$k_{\rm pivot}$ (1/Mpc):  
Maximum $\ell$ to compute powerspectrum up to:  
$\ell_{\rm max}$:  
(Precision) Number of points per oscillation to sample Bessel functions on:  
$n$:  
(Precision) Number of points per oscillation to sample $C_\ell$ integrand on:  
$n$:  
(Precision) Number of samples per oscillation for LOS integration:  
$n$:  
(Precision) Number of samples in $\log a$ for LOS integration:  
$n$:  
Compute temperature $C_\ell$'s:  
Compute Emode polarization $C_\ell$'s:  
Compute correlation function (powerspectrum we always compute):  
Comments (if something is not working / suggestions) max 200 characters:  
(can take a few seconds to load, just be patient!)