Abstract:This paper is concerned with a 1D Schrödinger scattering problem involving both oscillatory and evanescent regimes, separated by jump discontinuities in the potential function, to avoid "turning points". We derive 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. 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. We provide a complete error analysis of this hybrid method and illustrate our convergence results by numerical tests.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is about the numerical solution of the one - dimensional Schrödinger scattering problem in the semi - classical limit, especially involving the coupling of the oscillatory and decaying regions. Specifically:
- **Research Background**: When considering a one - dimensional Schrödinger equation with a piecewise - smooth potential function \(a(x)\), the solutions exhibit different behaviors in different regions. For the region where \(a(x)>0\), the solution is highly oscillatory; while for the region where \(a(x)<0\), the solution is exponentially growing or decaying. These different types of solutions pose challenges to numerical solutions.
- **Core of the Problem**: The paper aims to develop an efficient and accurate numerical method to deal with this multi - scale problem, especially in the semi - classical limit (i.e., a small parameter \(\varepsilon\ll1\)). Traditional methods require a very fine grid to accurately capture oscillations and boundary layers, which will lead to high computational costs. Therefore, the paper proposes a non - overlapping domain decomposition method (DDM), decomposes the original problem into sub - problems of oscillatory and decaying regions, and designs numerical methods suitable for these two regions respectively.
- **Specific Objectives**:
- For the oscillatory region, use the WKB - marching method based on the WKB approximation, which can obtain high - precision solutions on a coarse grid.
- For the decaying region, use the finite element method with WKB - basis functions (FEM with WKB - basis functions), which can also obtain high - precision solutions on a coarse grid.
- Provide a complete error analysis of the entire hybrid method and verify the convergence results through numerical experiments.
In summary, the goal of this paper is to develop and analyze an efficient numerical method for solving the one - dimensional Schrödinger equation containing oscillatory and decaying regions in the semi - classical limit, thereby avoiding the use of a very fine grid, reducing computational costs and improving the solution efficiency.
### Formula Display
The Schrödinger equation involved in the paper is:
\[ \varepsilon^2\psi''(x)+a(x)\psi(x) = 0,\quad x\in\mathbb{R}\]
where \(0 <\varepsilon\ll1\) is a small parameter and \(a(x)\) is a piecewise - smooth real function.
For the oscillatory region (\(a(x)>0\)), the wavelength is:
\[ \lambda(x)=\frac{2\pi\varepsilon}{\sqrt{a(x)}}\]
For the decaying region (\(a(x)<0\)), the wave function exhibits exponential growth or decay, and the typical thickness is:
\[ O\left(\varepsilon\sqrt{|a(x)|}\right)\]
The form of the WKB approximate solution is:
\[ \psi(x)\sim\exp\left(\frac{1}{\varepsilon}\sum_{p = 0}^{\infty}\varepsilon^p\phi_p(x)\right)\]
The WKB functions (the lowest three orders) in the oscillatory region are:
\[ \phi_{os}^{0}(x)=\exp\left(\pm\frac{i}{\varepsilon}\int_{0}^{x}\sqrt{a(y)}\,dy\right)\]
\[ \phi_{os}^{1}(x)=\exp\left(\pm\frac{i}{\varepsilon}\int_{0}^{x}\sqrt{a(y)}\,dy\right)\cdot\frac{1}{4\sqrt{a(x)}}\]
\[ \phi_{os}^{2}(x)=\exp\left(\pm\frac{i}{\varepsilon}\int_{0}^{x}[\sqrt{a(y)}-\varepsilon^2\beta(y)]\right)