Although finding numerically the quasinormal modes of a nonrotating black hole is a well-studied question, the physics of the problem is often hidden behind complicated numerical procedures aimed at avoiding the direct solution of the spectral system in this case. In this article, we use the exact analytical solutions of the Regge-Wheeler equation and the Teukolsky radial equation, written in terms of confluent Heun functions. In both cases, we obtain the quasinormal modes numerically from spectral condition written in terms of the Heun functions. The frequencies are compared with ones already published by Andersson and other authors. A new method of studying the branch cuts in the solutions is presented -the epsilon-method. In particular, we prove that the mode n = 8 is not algebraically special and find its value with more than 6 firm figures of precision for the first time. The stability of that mode is explored using the ǫ method, and the results show that this new method provides a natural way of studying the behavior of the modes around the branch cut points.
QUASI-NORMAL MODES OF BLACK HOLESThe study of quasinormal modes (QNMs) of a black hole (BH) has long history [1][2][3][4][5][6][7]. The reason behind this interest is that the QNMs offer a direct way of studying the key features of the physics of compact massive objects, without the complications of the full 3D general relativistic simulations. For example, by comparing the theoretically obtained gravitational QNMs with the frequencies of the gravitational waves, one can confirm or refute the nature of the central engines of many astrophysical objects, since those modes differ for the different types of objects -black holes, superspinars (naked singularities), neutron stars, black hole mimickers etc. [8][9][10][11][12][13].To find the QNMs, one needs to solve the second-order linear differential equations describing the linearized perturbations of the metric: the Regge-Wheeler equation (RWE) and the Zerilli equation for the Schwarzschild metric or the Teukolsky radial equation (TRE) for the Kerr metric and to impose the appropriate boundary conditions -the so-called black hole boundary conditions (waves going simultaneously into the horizon and into infinity) [1,3]. Additionally, one requires a regularity condition for the angular part of the solutions. And then, one needs to solve a connected problem with two complex spectral parameters -the frequency ω and the separation constant E (E = l(l + 1) -real for a nonrotating BH, with l the angular momentum of the perturbation). This Because of the complexity of the differential equations, until now, those equations were solved either approximately or numerically meeting an essential difficulty [1]. The indirect approaches like the continued fractions method have some limitations and are not directly related with the physics of the problem. The RWE, the Zerilli equation and TRE, however, can be solved analytically in terms of confluent Heun functions, as done for the first time in [16][17][18][19]. Imposin...