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.

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

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

A Legendre spectral method based on a hybrid format and its error estimation for fourth-order eigenvalue problems

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

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

Spectral Galerkin Approximation and Error Analysis Based on a Mixed Scheme for Fourth‐Order Problems in Complex Regions

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

Efficient Spectral-Galerkin Method I. Direct Solvers of Second- and Fourth-Order Equations Using Legendre Polynomials

Boundary element approximation of Steklov eigenvalue problem for Helmholtz equation

The Spectral Element Method for the Steklov Eigenvalue Problem

A Fast Fourier-Galerkin Method Solving Boundary Integral Equations for the Helmholtz Equation with Exponential Convergence

Conforming element method and shifted-inverse two-grid dicretization for the fourth-order Steklov eigenvalue problem

A spectral Galerkin exponential Euler time-stepping scheme for parabolic SPDEs on two-dimensional domains with a C2-boundary

Galerkin-Chebyshev Spectral Method and Block Boundary Value Methods for Two-Dimensional Semilinear Parabolic Equations

Fourth-order time-stepping for stiff PDEs on the sphere

Solving Sturm–Liouville inverse problems by an orthogonalized enhanced boundary function method and a product formula for symmetric potential

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

An Iterative Method Based on Equation Decomposition for the Fourth‐order Singular Perturbation Problem

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

The weak Galerkin finite element method for the Steklov eigenvalue problem

Split-step Spectral Galerkin Method for the Two-Dimensional Nonlinear Space-Fractional Schrödinger Equation