Spectral stability and perturbation results for kernel differentiation matrices on the sphere

We investigate the spectrum of differentiation matrices for certain operators on the sphere that are generated from collocation at a set of scattered points $X$ with positive definite and conditionally positive definite kernels. We focus on the case where these matrices are constructed from collocation using all the points in $X$ and from local subsets of points (or stencils) in $X$. The former case is often referred to as the global, Kansa, or pseudospectral method, while the latter is referred to as the local radial basis function (RBF) finite difference (RBF-FD) method. Both techniques are used extensively for numerically solving certain partial differential equations (PDEs) on spheres (and other domains). For time-dependent PDEs on spheres like the (surface) diffusion equation, the spectrum of the differentiation matrices and their stability under perturbations are central to understanding the temporal stability of the underlying numerical schemes. In the global case, we present a perturbation estimate for differentiation matrices which discretize operators that commute with the Laplace-Beltrami operator. In doing so, we demonstrate that if such an operator has negative (non-positive) spectrum, then the differentiation matrix does, too (i.e., it is Hurwitz stable). For conditionally positive definite kernels this is particularly challenging since the differentiation matrices are not necessarily diagonalizable. This perturbation theory is then used to obtain bounds on the spectra of the local RBF-FD differentiation matrices based on the conditionally positive definite surface spline kernels. Numerical results are presented to confirm the theoretical estimates.