In this talk I will discuss how to find and compute the eigenvalues of Dirac operators in their spectral gaps. In order to do so in an optimal way, the delicate study of the domains of critical Dirac operators is important and necessary. The results are concerned with variational methods, spectral theory and the development of optimal algorithms to compute the eigenvalues in a robust and efficient manner.