skypy.power_spectrum.transfer_no_wiggles

skypy.power_spectrum.transfer_no_wiggles(wavenumber, A_s, n_s, cosmology)

Eisenstein & Hu transfer function without wiggles. Eisenstein & Hu fitting formula for the transfer function without baryon acoustic oscillation wiggles. This is described in [1] and [2].

Parameters

wavenumber(nk,) array_like

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

A_s, n_s: float

Amplitude and spectral index of primordial scalar fluctuations.

cosmologyastropy.cosmology.Cosmology

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

Returns

transfer(nk, ) array_like

Transfer function evaluated at the given wavenumbers for the input primordial power spectrum parameters A_s and n_s, cosmology and kwmap normalization.

References

1

Eisenstein D. J., Hu W., ApJ, 496, 605 (1998)

2

Eisenstein D. J., Hu W., ApJ, 511, 5 (1999)

3

Komatsu et al., ApJS, 180, 330 (2009)

Examples

This returns the transfer function without wiggles for the Astropy default cosmology at a given array of wavenumbers:

>>> import numpy as np
>>> from astropy.cosmology import default_cosmology
>>> wavenumber = np.logspace(-3, 1, num=5, base=10.0)
>>> A_s, n_s = 2.1982e-09, 0.969453
>>> cosmology = default_cosmology.get()
>>> transfer_no_wiggles(wavenumber, A_s, n_s, cosmology)
array([9.91959695e-01, 7.84518347e-01, 1.32327555e-01, 4.60773671e-03,
   8.78447096e-05])