Numerical stability of the electromagnetic quasinormal and quasibound modes of Kerr black holes

The proper understanding of the electromagnetic counterpart of gravity waves emitters is of serious interest to the multimessenger astronomy. In this article, we study the numerical stability of the quasinormal modes (QNM) and quasibound modes (QBM) obtained as solutions of the Teukolsky Angular Equation and the Teukolsky Radial Equation with appropriate boundary conditions. We use the epsilon-method for the system featuring the confluent Heun functions to study the stability of the spectra with respect to changes in the radial variable. Ave find that the QNM. and QBM are stable in certain regions of the complex plane, just as expected, while the third "spurious" spectrum was found to be nithierically unstable and thus unphysical. This analysis shows the importance of understanding the numerical results in the framework of the theory of the confluent Heun functions, in order to be able to distinguish the physical spectra from the numerical artifacts.