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.

Rigorous estimation for the difference quotients of multiple eigenvalues

An Inexact Cayley Transform Method for Inverse Eigenvalue Problems with Multiple Eigenvalues

Rounding Errors of Partial Derivatives of Simple Eigenvalues of the Quadratic Eigenvalue Problem.

An Ulm-like Cayley Transform Method for Inverse Eigenvalue Problems with Multiple Eigenvalues

Estimates of the Gaps Between Consecutive Eigenvalues of Laplacian

Numerical optimisation of Dirac eigenvalues

2D Eigenvalue Problem II: Rayleigh Quotient Iteration and Applications

Convergence Analysis of Newton-Like Methods for Inverse Eigenvalue Problems with Multiple Eigenvalues

Fast computation of spectral densities for generalized eigenvalue problems

A numerical approach for investigating multiple eigenpairs of a quasilinear elliptic system

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

Isoperimetric Inequalities and the Gap Between the First and Second Eigenvalues of an Euclidean Domain

Asymptotically Exact a Posteriori Error Estimates of Eigenvalues by the Crouzeix-Raviart Element and Enriched Crouzeix-Raviart Element

Fastest quotient iteration with variational principles for self-adjoint eigenvalue problems

Optimal Domains for Elliptic Eigenvalue Problems with Rough Coefficients

Convergence of Q-Ritz Pairs, Refined Q-Ritz Vectors and Q-Rayleigh-ritz Method for Quadratic Eigenvalue Problems

First moment of Hecke eigenvalues at the integers represented by binary quadratic forms

General estimate of the first eigenvalue on manifolds

Computing eigenvalues of the Laplacian on rough domains