On the convergence of Krylov subspace methods for matrix Mittag-Leffler functions