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.

Non-consistent approximations of self-adjoint eigenproblems: Application to the supercell method

Possibility of magnetic field induced Composite Fermi sea in neutral Graphene

The A Posteriori Error Estimates of the FE Approximation of Defective Eigenvalues for Non-Self-Adjoint Eigenvalue Problems

Convergence of supercell and superspace methods for computing spectra of quasiperiodic operators

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

Inverse Sensitivity Method to Self-Adjoint Systems with Repeated Eigenvalues

Discontinuous Galerkin Methods of the Non-Selfadjoint 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

Numerical Approximations of a Nonlinear Eigenvalue Problem and Applications to a Density Functional Model

A High Accuracy Post-processing Algorithm for the Eigenvalues of Elliptic Operators

Optimal Convergence Rate of Self-Consistent Field Iteration for Solving Eigenvector-dependent Nonlinear Eigenvalue Problems

Theory and numerics of subspace approximation of eigenvalue problems

The Lower/Upper Bound Property of Approximate Eigenvalues by Nonconforming Finite Element Methods for Elliptic Operators

Convergence of Adaptive Finite Element Approximations for Nonlinear Eigenvalue Problems

Generalised Soft Finite Element Method for Elliptic Eigenvalue Problems

Asymptotic Errors in the Superconvergence of Discontinuous Galerkin Methods Based on Upwind-Biased Fluxes for 1D Linear Hyperbolic Equations

A posteriori error estimates for Schr{ö}dinger operators discretized with linear combinations of atomic orbitals

A staggered discontinuous Galerkin method for a class of nonlinear elliptic equations

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

Numerical approaches to compute spectra of non-self adjoint operators in dimensions two and three

When Is A Semiclassical Approximation Self-consistent?