Stationary Schrödinger equation in the semi-classical limit: numerical coupling of oscillatory and evanescent regions (Q681704)

This paper deals with the design, error analysis and numerics of an asymptotically correct scheme to numerically solve stationary, scalar Schrödinger equations \[ \epsilon^2\psi''(x)+a(x)\psi(x)=0,\quad x\in{\mathbb R}, \] in scattering situations, where \(0<\epsilon\ll 1\) and \(a\) denotes a piecewise (sufficiently) smooth function. This problem involves both oscillatory and evanescent regimes, separated by jump discontinuities in the potential function, to avoid ``returning points''. Quoting from the abstract, ``a non-overlapping domain decomposition method to split the original problem into sub-problems on these regions, both for the continuous and afterwards for the discrete problem'' is derived. ``Further, a hybrid WKB-based numerical method is designed for its efficient and accurate solution in the semi-classical limit: a WKB-marching method for the oscillatory regions and a FEM with WKB-basis functions in the evanescent regions.'' A complete error analysis of this hybrid method is provided. Numerical tests with a piecewise linear and a piecewise quadratic potential are given and illustrate the convergence analysis.