Linearly interpolating between 2 orientations represented by axis angle results in non-constant angular velocity. Spherical interpolation however doesn't suffer from this problem.

Why is that? Hopefully someone could point me to mathematical derivation. Thus far I am able to derive axis angle to matrix representation and back, but haven't found explanation for this.

Lerp vs slerp demo can be seen here. Demo has been made from this article