On their paper Puchwein et al. show their results for \(T_0\) and \(\gamma\), figures are below:

Comparing to my results, seems like there are differences in \(\gamma\) for the Puchwein model.

Optical Depth

Ionization Fraction

The H ionization fractions seem to match, so why are the optical depths so different for z>4.5?