High-order matrix method with delimited expansion domain

TL;DR

This paper introduces an advanced high-order matrix method using a mock-Chebyshev grid to accurately compute black hole quasinormal modes, especially in cases with metric discontinuities, improving stability and precision.

Contribution

It generalizes and enhances the matrix method for higher-order quasinormal mode calculations, addressing challenges posed by metric discontinuities.

Findings

The method achieves high accuracy in computing quasinormal modes.

It demonstrates improved convergence and speed over existing methods.

Effective handling of discontinuities in metric perturbations.

Abstract

Motivated by the substantial instability of the fundamental and high-overtone quasinormal modes, recent developments regarding the notion of black hole pseudospectrum call for numerical results with unprecedented precision. This work generalizes and improves the matrix method for black hole quasinormal modes to higher orders, specifically aiming at a class of perturbations to the metric featured by discontinuity intimately associated with the quasinormal mode structural instability. The approach is based on the mock-Chebyshev grid, which guarantees its convergence in the degree of the interpolant. In practice, solving for black hole quasinormal modes is a formidable task. The presence of discontinuity poses a further difficulty so that many well-known approaches cannot be employed straightforwardly. Compared with other viable methods, the modified matrix method is competent in speed and accuracy. Therefore, the method serves as a helpful gadget for relevant studies.