We use a lattice formulation to study the spectra of the Dirac and the Dirac–Kähler operators on the 2-sphere. This lattice formulation uses differentiation matrices which yield exact values for the derivatives of polynomials, preserving the Leibniz rule in subspaces of polynomials of low degree and therefore, this formulation can be used to study the fermion–boson symmetry on the lattice. In this context, we find that the free Dirac and Dirac–Kähler operators on the 2-sphere exhibit fermionic as well as bosonic-like eigensolutions for which the corresponding eigenvalues and the number of states are computed. In the Dirac case these solutions appear in doublets, except for the bosonic mode with zero eigenvalue, indicating the possible existence of a supersymmetry of the square of the Dirac operator.