skypy.power_spectrum.camb

skypy.power_spectrum.camb(wavenumber, redshift, cosmology, A_s, n_s)

CAMB linear matter power spectrum. Compute the linear matter power spectrum on a two dimensional grid of redshift and wavenumber using CAMB [1].

Parameters

wavenumber(nk,) array_like

Array of wavenumbers in units of Mpc-1 at which to evaluate the linear matter power spectrum.

redshift(nz,) array_like

Array of redshifts at which to evaluate the linear matter power spectrum.

cosmologyastropy.cosmology.Cosmology

Cosmology object providing omega_matter, omega_baryon, Hubble parameter and CMB temperature in the present day

A_sfloat

Cosmology parameter, amplitude normalisation of curvature perturbation power spectrum

n_sfloat

Cosmology parameter, spectral index of scalar perturbation power spectrum

Returns

power_spectrum(nz, nk) array_like

Array of values for the linear matter power spectrum in Mpc3 evaluated at the input wavenumbers for the given primordial power spectrum parameters, cosmology. For nz redshifts and nk wavenumbers the returned array will have shape (nz, nk).

References

1

Lewis, A. and Challinor, A. and Lasenby, A. (2000), doi : 10.1086/309179.

Examples

>>> import numpy as np
>>> from astropy.cosmology import default_cosmology
>>> cosmology = default_cosmology.get()
>>> redshift = np.array([0, 1])
>>> wavenumber = np.array([1.e-2, 1.e-1, 1e0])
>>> A_s = 2.e-9
>>> n_s = 0.965
>>> power_spectrum = camb(wavenumber, redshift, cosmology, A_s, n_s)