Complete quasinormal modes of type-D black holes

Quasinormal mode (QNM) spectra of black holes exhibit two open problems [Conf. Proc. C 0405132, 145 (2004); Classical Quantum Gravity 26, 163001 (2009) ]: (i) the discontinuity in highly damped QNMs between Schwarzschild and Kerr solutions as $a\ensuremath{\rightarrow}0$, and (ii) the unexplained spectral proximity between QNMs and algebraically special (AS) frequencies, particularly the anomalous multiplet splitting for Kerr $\ensuremath{\ell}=2$, $m\ensuremath{\ge}0$ modes. We develop a novel method to compute complete QNM spectra for type-D black holes, solving both problems and establishing a mathematical framework for boundary value problems of dissipative systems. Using analytic continuation of radial eigenvalue equations, our method eliminates the dependence on auxiliary parameters in the connection formulas for confluent Heun solutions. This breakthrough overcomes the longstanding challenge of calculating QNMs that cross or lie on the negative imaginary axis (NIA). For Schwarzschild and Kerr spacetimes ($0\ensuremath{\le}n\ensuremath{\le}41$, $2\ensuremath{\le}\ensuremath{\ell}\ensuremath{\le}16$), we present complete spectra validated through scattering amplitudes with errors $<{10}^{\ensuremath{-}10}$. The results provide definitive solutions to both open problems: (i) The inability of conventional methods to compute QNMs crossing or residing on the NIA leads to apparent discontinuities in Kerr spectra as $a\ensuremath{\rightarrow}0$. (ii) When a QNM exactly coincides with the AS frequency, an additional QNM (unconventional mode) appears nearby. For the Kerr case with $\ensuremath{\ell}=2$, overtone sequences from both unconventional and AS modes exhibit precisely $2\ensuremath{\ell}+1$ branches, without multiplets or supersymmetry breaking at AS frequencies. Our method confirms Leung's conjecture of high-$\ensuremath{\ell}$ unconventional mode deviations from the NIA through the first $\ensuremath{\ell}=3$ calculation. Moreover, this paradigm surpasses the state-of-the-art Cook's continued fraction and isomonodromic methods in computational efficiency.