Numerical solutions (Green's functions) for intermediate-type distortions and Mathematica code
Spectrum created by energy release in the early Universe, before
recombination, creates distortions which are a superposition of
μ-type, y-type and intermediate-type distortions, shown in the above
figure. The final spectrum can thus be constructed from the templates,
once energy injection rate as a function of redshift is known. This
package contains the spectral distortions (Green's functions) spaced at dy=0.001 for y<1 and dy=0.01
for y>1 covering a range 0.001 < y < 10.
Also included is a mathematica
code which can combine these numerical solutions of Kompaneets
equation (Green's functions) for user-defined rate of
energy injection as a function of redshift. Silk damping, particle
decay and annihilation examples are also included. More details can
README file and the companion paper
See also Appendix of arXiv:1303.7212 for the algorithm and how to use the Green's functions provided to write your own code.