paper

Bernstein spectral method for quasinormal modes and other eigenvalue problems

arXiv:2003.06232 · doi:10.1140/epjc/s10052-023-12350-9

Abstract

Spectral methods are now common in the solution of ordinary differential eigenvalue problems in a wide variety of fields, such as in the computation of black hole quasinormal modes. Most of these spectral codes are based on standard Chebyshev, Fourier, or some other orthogonal basis functions. In this work we highlight the usefulness of a relatively unknown set of non-orthogonal basis functions, known as Bernstein polynomials, and their advantages for handling boundary conditions in ordinary differential eigenvalue problems. We also report on a new user-friendly package, called \texttt{SpectralBP}, that implements Berstein-polynomial-based pseudospectral routines for eigenvalue problems. We demonstrate the functionalities of the package by applying it to a number of model problems in quantum mechanics and to the problem of computing scalar and gravitational quasinormal modes in a Schwarzschild background. We validate our code against some known results and achieve excellent agreement. Compared to continued-fraction or series methods, global approximation methods are particularly well-suited for computing purely imaginary modes such as the algebraically special modes for Schwarzschild gravitational perturbations.

arXiv comment46 pages, 9 figures added references, changed format, reorganized sections to enhance flow of discussions, improved the subsection discussing the subtleties in computing the ASM, strengthened discussion on the QHO, edits to graphs so that they are readable in black and white, fixed some typos

References in corpus (14)

Cited by in corpus (12)