Limiting set of second order spectra
Lyonell Boulton
TL;DR
The paper develops and rigorously analyzes a second order spectrum method for detecting isolated eigenvalues of self-adjoint operators in gaps of the essential spectrum. By defining the second order spectrum $\mathrm{Spec}_2(M|\mathcal{L}_n)$ via zeros of $\det[P_n(M-z)^2|\mathcal{L}_n]$ and its Lipschitz scalar surrogate $\sigma_n(z)$, it constructs the limiting set $\Lambda$ and proves that, under the mild hypothesis (H), every isolated eigenvalue of finite multiplicity lies in $\Lambda$, while the spectra avoid spurious pollution away from the real axis. The paper also provides a concrete model with $M$ as a piecewise-constant multiplier where $\Lambda(M)=\mathbb{T}$, illustrating maximal pollution in linear Ritz methods and the non-polluting behavior of the second order approach. Finally, it offers practical criteria to verify (H) for unbounded operators, notably via a dominating operator with compact resolvent (e.g., the harmonic oscillator) and an explicit analysis for $M=-\partial_x^2+V$. This combination of theory and explicit examples demonstrates that second order spectra can achieve simultaneous non-pollution and accurate eigenvalue approximation in broad settings.
Abstract
M.Levitin and E.Shargorodsky purposed in a recent article, [math.SP/0212087], the use of the so called ``second order relative spectrum'', to find eigenvalues of self-adjoint operators in gaps of the essential spectrum. Let $M$ be a self-adjoint operator acting on the Hilbert space $H$. A complex number $z$ is in the second order spectrum of $M$ relative to a finite dimensional subspace $Λ\subset dom M^2$ if and only if the truncation to $Λ$ of $(M-z)^2$ is not invertible. It is remarkable that these sets seem to provide a general method to estimating eigenvalues free from the problems of spectral pollution present in most linear methods. In this notes we investigate rigorously various aspects of the use of second order spectrum to finding eigenvalues. Our main result shows that, under certain fairly mild hypothesis on $M$, the uniform limit of the second order spectra, as $Λ$ increases toward $H$, contains the isolated eigenvalues of $M$ of finite multiplicity. In applications the essential spectrum can be computed analytically, while precisely these eigenvalues are the ones that should be approximated numerically. Hence this method seems to combine non-pollution and approximation at a very high level of generality.
