Linearized model collision operators for multiple ion species plasmas and gyrokinetic entropy balance equations