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.

An Adaptive Inverse Iteration for Maxwell Eigenvalue Problem Based on Edge Elements.

Convergence Analysis of an Adaptive Edge Element Method for Maxwell's Equations

An Adaptive Edge Element Method and Its Convergence for a Saddle‐Point Problem from Magnetostatics

Edge Element Methods for Maxwell's Equations with Strong Convergence for Gauss' Laws.

A Meshless Method for Some Inverse Problems Associated with the Helmholtz Equation

An Adaptive High-Order Transient Algorithm to Solve Large-Scale Anisotropic Maxwell’s Equations

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

An Efficient Adaptive Finite Element Method for Eigenvalue Problems

Finite Element Formulations for Maxwell's Eigenvalue Problem Using Continuous Lagrangian Interpolations

An adaptive algorithm based on the shifted inverse iteration for the Steklov eigenvalue problem

An adaptive finite element method in reconstruction of coefficients in Maxwell's equations from limited observations

Optimized Operator-Splitting Methods in Numerical Integration of Maxwell's Equations

An Adaptive Edge Finite Element DtN Method for Maxwell's Equations in Biperiodic Structures

An adaptive finite element method in reconstruction of coefficients in Maxwell’s equations from limited observations

Finite Element Approximation of the Modified Maxwell's Stekloff Eigenvalues

Adaptive finite element approximation for steady-state Poisson-Nernst-Planck equations

Convergence of Adaptive Finite Element Approximations for Nonlinear Eigenvalue Problems

A General Adaptive Finite Element Eigen-Algorithm Stemming from Wittrick-Williams Algorithm

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

The Shifted-Inverse Iteration Based on the Multigrid Discretizations for Eigenvalue Problems