Efficient Spectral Methods for Some Singular Eigenvalue Problems

Abstract

We propose and analyze some efficient spectral/spectral element methods to solve singular eigenvalue problems related to the Schrodinger operator with an inverse-power potential. For the Schrodinger eigenvalue problem $$-\Delta u +V(x)u=\lambda u$$ with a regular potential $$V(x)=c_1|x|^{-1}$$ , we first design an efficient spectral method on a ball of any dimension by adopting the Sobloev-orthogonal basis functions with respect to the Laplacian operator to overwhelm the homogeneous inverse potential and to eliminate the singularity of the eigenfunctions. Then we extend this spectral method to arbitrary polygonal domains by the mortar element method with each corner covered by a circular sector and origin covered by a circular disc. Furthermore, for the Schrodinger eigenvalue problem with a singular potential $$V(x)=c_3|x|^{-3}$$ , we devise a novel spectral method by modifying the former Sobloev-orthogonal bases to fit the stronger singularity. As in the case of $$|x|^{-1}$$ potential, this approach can be extended to arbitrary polygonal domains by the mortar element method as well. Finally, for the singular elliptic eigenvalue problem $$-\frac{\partial ^2}{\partial x^2}u-\frac{1}{x^2}\frac{\partial ^2}{\partial y^2}u =\lambda u$$ on rectangles, we propose a novel spectral method by using tensorial bases composed of the $$L^2$$ - and $$H^1$$ -simultaneously orthogonal functions in the y-direction and the Sobolev-orthogonal functions with respect to the Schrodinger operator with an inverse-square potential in the x-direction. Numerical experiments indicate that all our methods possess exponential orders of convergence, and are superior to the existing polynomial based spectral/spectral element methods and hp-adaptive methods.