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.

High-Accuracy Numerical Approximations to Several Singularly Perturbed Problems and Singular Integral Equations by Enriched Spectral Galerkin Methods

Enriched Spectral Methods and Applications to Problems with Weakly Singular Solutions

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

An Efficient and Accurate Numerical Method for the Spectral Fractional Laplacian Equation

Thin plate spline Galerkin scheme for numerically solving nonlinear weakly singular Fredholm integral equations

Numerical Solution of Singularly Perturbed Problems via both Galerkin and Subdomain Galerkin methods

Iterative Solvers for a Spectral Galerkin Approach to Elliptic Partial Differential Equations with Fuzzy Coefficients

Convergence analysis of a spectral-Galerkin-type search extension method for finding multiple solutions to semilinear problems

Efficient and Accurate Spectral Method Using Generalized Jacobi Functions for Solving Riesz Fractional Differential Equations

Enrichment terms for line singularities in second-order elliptic boundary value problems in 3D heterogeneous media: Application to heat conduction problems

Understanding the ultraspherical spectral method

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

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

Computing Nearly Singular Solutions Using Pseudo-Spectral Methods

Error estimation for non-linear singularly perturbed reaction–diffusion parabolic problems via element-free Galerkin method

Efficient spectral–Galerkin methods for fractional partial differential equations with variable coefficients

A Highly Accurate and Efficient Genocchi-based Spectral Technique Applied to Singular Fractional Order Boundary Value Problems

Space‐time Spectral Method for Two‐dimensional Semilinear Parabolic Equations

A parameter robust numerical method for a singularly perturbed Volterra equation in security technologies

An Exponential Convergence Approximation to Singularly Perturbed Problems by Log Orthogonal Functions

Pseudo-Spectral Galerkin Method Using Shifted Vieta-Fibonacci Polynomials for Stochastic Models: Existence, Stability, and Numerical Validation