© 2019 Author(s). The linearized model collision operator for multiple species plasmas given by Sugama et al. [Phys. Plasmas 16, 112503 (2009)] is improved to be properly applicable up to the highly collisional regime. The improved linearized model operator retains the conservation laws of particles, momentum, and energy, and it reproduces the same friction-flow relations as derived by the linearized Landau operator so that this model can be used to correctly evaluate neoclassical transport fluxes in all collisionality regimes. The adjointness relations and Boltzmann's H-theorem are exactly satisfied by the improved operator except in the case of collisions between unlike particle species with unequal temperatures where these relations and H-theorem still hold approximately because there is a large difference between the masses of the two species with significantly different temperatures. Even in the unequal-temperature case, the improved operator can also be modified so as to exactly satisfy the adjointness relations, while it causes the values of the friction coefficients to deviate from those given by the Landau operator. In addition, for application to gyrokinetic simulations of turbulent transport, the improved operator is transformed into the gyrophase-averaged form by keeping the finite gyroradius effect.