FEM-DtN-SIM Method for Computing Resonances of Schr\"{o}dinger Operators

Abstract

The study of resonances of the Schr\"{o}dinger operator has a long-standing tradition in mathematical physics. Extensive theoretical investigations have explored the proximity of resonances to the real axis, their distribution, and bounds on the counting functions. However, computational results beyond one dimension remain scarce due to the nonlinearity of the problem and the unbounded nature of the domain. We propose a novel approach that integrates finite elements, Dirichlet-to-Neumann (DtN) mapping, and the spectral indicator method. The DtN mapping, imposed on the boundary of a truncated computational domain, enforces the outgoing condition. Finite elements allow for the efficient handling of complicated potential functions. The spectral indicator method effectively computes (complex) eigenvalues of the resulting nonlinear algebraic system without introducing spectral pollution. The viability of this approach is demonstrated through a range of numerical examples.