Abstract:In this article, we study numerical approximation of eigenvalue problems of the Schrödinger operator $-\Delta u + \frac{c^2}{|x|^2}u$. There are three stages in our investigation: We start from a ball of any dimension, in which case the exact solution in the radial direction can be expressed by Bessel functions of fractional degrees. This knowledge helps us to design two novel spectral methods by modifying the polynomial basis to fit the singularities of the eigenfunctions. At the second stage, we move to circular sectors in the two-dimensional setting. Again the radial direction can be expressed by Bessel functions of fractional degrees. Only in the tangential direction some modifications are needed from stage one. At the final stage, we extend the idea to arbitrary polygonal domains. We propose a mortar spectral element approach: a polygonal domain is decomposed into several subdomains with each singular corner including the origin covered by a circular sector, in which origin and corner singularities are handled similarly as in the former stages, and the remaining domains are either a standard quadrilateral/triangle or a quadrilateral/triangle with a circular edge, in which the traditional polynomial based spectral method is applied. All subdomains are linked by mortar elements (note that we may have hanging nodes). In all three stages, exponential convergence rates are achieved. Numerical experiments indicate that our new methods are superior to standard polynomial based spectral (or spectral element) methods and $hp$-adaptive methods. Our study offers a new and effective way to handle eigenvalue problems of the Schrödinger operator including the Laplacian operator on polygonal domains with reentrant corners.

Bernstein spectral method for quasinormal modes and other eigenvalue problems

Efficient Chebyshev Spectral Method for Solving Linear Elliptic Pdes Using Quasi-Inverse Technique

A Bernstein Polynomial Collocation Method for the Solution of Elliptic Boundary Value Problems

Bernstein collocation technique for a class of Sturm-Liouville problems

Overdamped modes in Schwarzschild-de Sitter and a Mathematica package for the numerical computation of quasinormal modes

Spectral Methods for Solution of Differential and Functional Equations

Computing eigenfunctions and eigenvalues of boundary value problems with the orthogonal spectral renormalization method

Efficient Spectral and Spectral Element Methods for Eigenvalue Problems of Schrödinger Equations with an Inverse Square Potential

Pseudo-Spectral Galerkin Method Using Shifted Vieta-Fibonacci Polynomials for Stochastic Models: Existence, Stability, and Numerical Validation

Natural polynomials for Kerr quasi-normal modes

Hyperboloidal Method for Quasinormal Modes of Non-Relativistic Operators

Bernstein approximation and beyond: proofs by means of elementary probability theory

A New Spectral Method Using Nonstandard Singular Basis Functions for Time-Fractional Differential Equations

Numerical solution of fractional Fredholm integro-differential equations by spectral method with fractional basis functions

Spectral Jacobi approximations for Boussinesq systems

Fast computation of spectral densities for generalized eigenvalue problems

Efficient Spectral Methods for Quasi-Equilibrium Closure Approximations of Symmetric Problems on Unit Circle and Sphere

Quantum eigenvalue processing

Mathematical foundations of spectral methods for time-dependent PDEs