ColRadPy icon indicating copy to clipboard operation
ColRadPy copied to clipboard

Use cubic interpolation of adf11 rates for solving ionization balances

Open dmkriete opened this issue 9 months ago • 0 comments

This switches the interpolation of adf11 rates from linear to cubic (the interpolation is still done in logarithmic space).

For adf11 rates this appears to give better results. See comparison below; the linear interpolation exhibits unphysical "bumps" in the curves that disappear when cubic interpolation is used. (technically the below comparison is technically energy balance, but ionization balance is a component of energy balance, so comparison is applicable to both).

Linear interpolation (old) c2+ thermalization with protium for Ti=1 0Te_old

Cubic interpolation (new) c2+ thermalization with protium for Ti=1 0Te_new

dmkriete avatar May 10 '24 17:05 dmkriete