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.

The Spectral Element Method for the Steklov Eigenvalue Problem

The a Posteriori Error Estimates of the Spectral Element Method for the Steklov Eigenvalue Problem

Spectral Method with the Tensor-Product Nodal Basis for the Steklov Eigenvalue Problem

An Efficient Spectral-Galerkin Method for a New Steklov Eigenvalue Problem in Inverse Scattering

A Stochastic Finite Element Method for Real Eigenvalue Problems

Conforming Finite Element Approximations for A Fourth-Order Steklov Eigenvalue Problem

The weak Galerkin finite element method for the Steklov eigenvalue problem

The a Priori and a Posteriori Error Estimates of DG Method for the Steklov Eigenvalue Problem in Inverse Scattering

A Priori and a Posteriori Error Estimates for a Virtual Element Method for the Non-Self-adjoint Steklov Eigenvalue Problem

Boundary element approximation of Steklov eigenvalue problem for Helmholtz equation

The Adaptive Finite Element Method for the Steklov Eigenvalue Problem in Inverse Scattering

A Highly Efficient Spectral-Galerkin Method Based on Tensor Product for Fourth-Order Steklov Equation with Boundary Eigenvalue

Local A Priori/A Posteriori Error Estimates of Conforming Finite Elements Approximation for Steklov Eigenvalue Problems

An adaptive algorithm for Steklov eigenvalues of the Lamé operator in linear elasticity

A Multigrid Correction Scheme for a New Steklov Eigenvalue Problem in Inverse Scattering

Two-grid Discretizations and a Local Finite Element Scheme for a Non-Selfadjoint Stekloff Eigenvalue Problem

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

Discontinuous Galerkin Methods of the Non-Selfadjoint Steklov Eigenvalue Problem in Inverse Scattering

Matlab Experiments on Extrapolation of The Nonconforming Crouzeix-Raviart Element for Steklov Eigenvalue Problem

A two-grid method of the non-conforming Crouzeix–Raviart element for the Steklov eigenvalue problem

Finite Element Approximation of the Modified Maxwell's Stekloff Eigenvalues