CAMB: replacing reionization history

Use of Healpix, camb, CLASS, cosmomc, compilers, etc.
Post Reply
Shek Yeung
Posts: 2
Joined: February 18 2009
Affiliation: CUHK

CAMB: replacing reionization history

Post by Shek Yeung » February 18 2009

I'm doing research about a model with non-tanh reionization history. I already changed the function Reionization_xe in reionization.f90, but now I'm not sure about the value of ReionHist%tau_start. The comment says it's the time when xe is 'very small', how should I change it if the reionization history is not tanh anymore?

Another question is about the calculation of actual_opt_depth. I've read the subroutine inithermo but I don't understand how it's calculated. What is the theory behind this value? Any recommended readings?

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

Re: CAMB: replacing reionization history

Post by Antony Lewis » February 20 2009

Any value where your x_e is very small should be OK (e.g. 1e-4).

The "actual" optical depth is the calculated optical depth due to reionization (which may differ from input value since that is mapped approximately onto a reionization redshift). It's the integral using the difference of the total x_e to the x_e there would be with no reionization (so is well defined, not depending on an arbitrary high-z cut-off, and zero when no reionization even though there is finite scattering due to the ~2e-4 ionization fraction left from recombination).

Post Reply