## CAMB: transfer function for gamma at z~100

Use of Healpix, camb, CLASS, cosmomc, compilers, etc.
Kyungjin Ahn
Posts: 4
Joined: May 02 2008
Affiliation: Chosun University
Contact:

### CAMB: transfer function for gamma at z~100

Pardon me if this has been addressed before. I'm using January 2015 version of CAMB, and wanted to use the 4th column (photon density perturbation) of transfer function outputs at various redshifts. I set
transfer_kmax = 20
in the parameter file.

With redshifts selected as in the figure file(1000, 203, 101), the transfer function has discontinuous points at z=203 and 101 which looks suspicious. I ran it both with intel-fortran compiled and gfortran-compiled ones with the same result.

As a test I went back to September 2013 version, and the same outputs for the same redshifts have null values for the transfer function at k>~0.02, so it used to be worse.

Is it a glitch from error accumulation of the ODE solver used?

[/b]

Antony Lewis
Posts: 1636
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

### Re: CAMB: transfer function for gamma at z~100

By default it doesn't calculate the radiation multipoles accurately at late times, the glitches may be the sudden switch to the approximate solution. You can turn off DoLateRadTruncation and increase lAccuracyBoost if you want to get it more accurately, e.g. (in latest CAMB with python wrapper)

Code: Select all

options=&#91;&#123;&#125;,&#123;'DoLateRadTruncation'&#58;False&#125;,&#123;'DoLateRadTruncation'&#58;False, 'lAccuracyBoost'&#58;4&#125;&#93;
_, axs = plt.subplots&#40;1,3, figsize=&#40;12,4&#41;&#41;
for option, ax in zip&#40;options, axs&#41;&#58;
pars = camb.CAMBparams&#40;&#41;
pars.set_cosmology&#40;H0=67.5, ombh2=0.022, omch2=0.122&#41;
pars.InitPower.set_params&#40;ns=0.965&#41;
pars.WantCls = False
pars.NonLinear = model.NonLinear_none

pars.set_accuracy&#40;**option&#41;

results = camb.get_results&#40;pars&#41;

M = results.get_matter_transfer_data&#40;&#41;
for i in range&#40;3&#41;&#58;
tz = M.transfer_z&#40;'delta_photon',i&#41;
ax.plot&#40;M.q, tz*M.q**2&#41;
ax.set_xlim&#40;&#91;0,0.1&#93;&#41;
gives

Kyungjin Ahn
Posts: 4
Joined: May 02 2008
Affiliation: Chosun University
Contact:

### CAMB: transfer function for gamma at z~100

Thank you so much, Antony. Not having the proper gfortran version for python wrapper, I just followed your instruction inside the usual input file. Now the figure is what I have with