Exploring the Pseudo-modes of Schrödinger Operators with Complex Potentials: A Focus on Resolvent Norm Estimates and Spectral Stability
Sameh Gana
TL;DR
The paper analyzes the ε-pseudospectrum of the non-self-adjoint complex harmonic oscillator $H_c=-\dfrac{d^{2}}{dx^{2}}+cx^{2}$ with $Re(c)>0$, $Im(c)\neq0$, focusing on resolvent-norm behavior along curves $z_{\eta}=b\eta+c\eta^{p}$ and, in particular, the border case $p=\tfrac{1}{3}$. Employing semiclassical scaling and microlocal factorization, it proves that for $p=\tfrac{1}{3}$ the resolvent norm satisfies $\|(H_c-zI)^{-1}\|=O(|z|^{-1/3})$ as $|z|\to\infty$ along $z_{\eta}$, indicating spectral stability in this regime; this contrasts with known blow-up results for $\tfrac{1}{3}<p<3$. The numerical portion, using Chebyshev collocation and pseudo-spectral techniques, corroborates the theoretical predictions and demonstrates pronounced pseudospectral growth far from the stability line, signaling high-energy instability under perturbations. Overall, the work clarifies the border between stability and instability in the pseudospectrum of complex potentials and provides a methodology for validating spectral sensitivity via numerical pseudo-spectra.
Abstract
This paper aims to investigate the pseudo-modes of the one-dimensional Schrödinger operator with complex potentials, focusing on the behavior of the resolvent norm along specific curves in the complex plane and assessing the stability of the spectrum under small perturbations. The study builds upon previous work of E.B. Davies, L.S. Boulton, and N. Trefethen, specifically examining the resolvent norm of the complex harmonic oscillator along curves of the form $z_{η}= bη+ cη^{p} $ where $ b > 0$, $ \frac{1}{3}< p <3 $ independent of $η> 0$. The present work narrows the focus to the case where $p = \frac{1}{3}$. Numerical computations of pseudo-eigenvalues are performed to verify spectral instability.
