Abstract:This paper is dedicated to the efficient numerical computation of solutions to the 1D stationary Schrödinger equation in the highly oscillatory regime. We compute an approximate solution based on the well-known WKB-ansatz, which relies on an asymptotic expansion w.r.t. the small parameter $\varepsilon$. Assuming that the coefficient in the equation is analytic, we derive an explicit error estimate for the truncated WKB series, in terms of $\varepsilon$ and the truncation order $N$. For any fixed $\varepsilon$, this allows to determine the optimal truncation order $N_{opt}$ which turns out to be proportional to $\varepsilon^{-1}$. When chosen this way, the resulting error of the optimally truncated WKB series behaves like $\mathcal{O}(\exp(-r/\varepsilon))$, with some parameter $r>0$. The theoretical results established in this paper are confirmed by several numerical examples.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the efficient numerical computation of the one - dimensional stationary Schrödinger equation in a highly oscillatory regime. Specifically, the authors focus on how to effectively solve the one - dimensional stationary Schrödinger equation in the following form when the parameter $ \varepsilon $ is very small: \[ \begin{cases} \varepsilon^2 \varphi''(x) + a(x) \varphi(x) = 0, & x \in I := [\xi, \eta], \\ \varphi(\xi) = \varphi_0, \\ \varepsilon \varphi'(\xi) = \varphi_1, \end{cases} \] where $ 0 < \varepsilon \ll 1 $ is a very small parameter, and $ a(x) \geq a_0 > 0 $ is a real - valued function, which is related to the potential energy in quantum mechanics problems. Since $ \varepsilon $ is very small, the solution $ \varphi(x) $ becomes highly oscillatory, especially in the semi - classical limit $ \varepsilon \to 0 $. To solve this problem, the authors construct an approximate solution based on the well - known WKB approximation method (Wentzel - Kramers - Brillouin approximation) by using the asymptotic expansion with respect to the small parameter $ \varepsilon $. They derive the explicit error estimate of the truncated WKB series and determine the optimal truncation order $ N_{\text{opt}} $, which is proportional to $ \varepsilon^{-1} $. After choosing this optimal truncation order, the error behavior of the truncated WKB series is $ O(\exp(-r/\varepsilon)) $, where $ r > 0 $ is a certain parameter. In addition, the authors verify the validity of the theoretical results through numerical experiments. Their goal is to develop an efficient numerical method so that when $ \varepsilon $ is very small, the highly oscillatory one - dimensional stationary Schrödinger equation can be solved with high precision.

Optimally truncated WKB approximation for the 1D stationary Schrödinger equation in the highly oscillatory regime

Optimally truncated WKB approximation for the highly oscillatory stationary 1D Schrödinger equation

WKB-based third order method for the highly oscillatory 1D stationary Schrödinger equation

High-order WKB-based Method For The 1D Stationary Schrödinger Equation In The Semi-classical Limit

Uniformly accurate numerical schemes for highly oscillatory Klein-Gordon and nonlinear Schrödinger equations

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

Optimal Error Estimate of the Extended-WKB Approximation to the High Frequency Wave-Type Equation in the Semi-classical Regime

On numerical solutions of the time-dependent Schrödinger equation

Reducibility of 1-D Schrödinger Equation with Unbounded Oscillation Perturbations

Asymptotics for the cubic 1D NLS with a slowly decaying potential

Absorbing Boundary Conditions for Variable Potential Schrödinger Equations via Titchmarsh-Weyl Theory

Minimizing travelling waves for the one-dimensional nonlinear Schrödinger equation with non-zero condition at infinity

Optimal error bounds on time-splitting methods for the nonlinear Schrödinger equation with low regularity potential and nonlinearity

Standing wave solutions of a quasilinear Schrödinger equation in the small frequency limit

High-order asymptotic expansion for the nonlinear Klein-Gordon equation in the non-relativistic limit regime

Uniform error bound of an exponential wave integrator for the long-time dynamics of the nonlinear Schrödinger equation with wave operator

Stability and Error Analysis for a Second-Order Fast Approximation of the One-dimensional Schrödinger Equation under Absorbing Boundary Conditions

Tenth-Order Accurate Numerical Method for Solving the Time-Dependent Schrödinger Equation

Numerological analysis of the WKB approximation in large order

An Efficient Second-Order Finite Difference Method for the One-Dimensional Schrödinger Equation with Absorbing Boundary Conditions

Explicit summation of the constituent WKB series and new approximate wave functions