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.

Exponential Collocation Methods Based on Continuous Finite Element Approximations for Efficiently Solving the Cubic Schrödinger Equation

Exponential Collocation Methods for the Cubic Schrödinger Equation

Arbitrary-order Energy-Preserving Exponential Integrators for the Cubic Schrödinger Equation

Scattering for the Focusing Energy-Subcritical Nonlinear Schrödinger Equation

An Extended Subequation Rational Expansion Method and Solutions of (2+1)-Dimensional Cubic Nonlinear Schrödinger Equation

Symplectic and Multi-Symplectic Wavelet Collocation Methods for Two-Dimensional Schrödinger Equations

A new general algebraic method with symbolic computation to construct new exact analytical solution for a (2 + 1)-dimensional cubic nonlinear Schrödinger equation

Exponential Collocation Methods for Conservative or Dissipative Systems.

Strong convergence rate of finite difference approximations for stochastic cubic Schrödinger equations

Accurate and Efficient Numerical Methods for the Nonlinear Schrödinger Equation with Dirac Delta Potential

A Time Two-Mesh Compact Difference Method for the One-Dimensional Nonlinear Schrödinger Equation

Novel Conservative Methods for Schrödinger Equations with Variable Coefficients over Long Time

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

Mass- and Energy-Preserving Exponential Runge-Kutta Methods for the Nonlinear Schrodinger Equation

Large Global Solutions for the Energy-Critical Nonlinear Schrödinger Equation.

A discrete least squares collocation method for two-dimensional nonlinear time-dependent partial differential equations

An Efficient Spectral-Collocation Difference Method for Two-Dimensional Schrödinger Equation with Neumann Boundary Conditions

Numerical Analysis and Computation of the Finite Volume Element Method for the Nonlinear Coupled Time-Fractional Schrödinger Equations

An Efficient Numerical Approach to Solve Schrödinger Equations with Space Fractional Derivative

Large global solutions for nonlinear Schr"odinger equations III, energy-supercritical cases

Numerical Solutions of Schrödinger Equations in ℝ3