“…Furthermore, these solutions have been utilized to calculate QNMs of the Teukolsky equation describing the Schwarzschild black hole [38], including its continuous spectrum [39], as well as for the central engine of Gamma-ray bursts (GRB) and cosmic jet of the Kerr black hole [40,41]. Following Fiziev's work, numerous researchers have employed Heun class functions to study source-free perturbation equations in various spacetimes, such as QNMs [42] of Dirac field in 2+1 dimensional GW background [43], QNMs and the reflection coefficient [44] of massless fields in Kerr-Newman-de Sitter BHs [45], and QNMs of the massive scalar field in Kerr-AdS 5 BHs [46]. Moreover, Cook et al [47] converted the modes equations of Kerr BHs into the confluent Heun equation and then solved the radial equation using the continued fraction method, and the angular equation using the spectral method, resulting in improved accuracy for the QNMs.…”