Structured Backward Error for the WKB method

Robert M. Corless,Nicolas Fillion

2024-12-01

Abstract:The classical WKB method (also known as the WKBJ method, the LG method, or the phase integral method) for solving singularly perturbed linear differential equations has never, as far as we know, been looked at from the structured backward error (BEA) point of view. This is somewhat surprising, because a simple computation shows that for some important problems, the WKB method gives the exact solution of a problem of the same structure that can be expressed in finitely many terms. This kind of analysis can be extremely useful in assessing the validity of a solution provided by the WKB method. In this paper we show how to do this and explore some of the consequences, which include a new iterative algorithm to improve the quality of the WKB solution. We also explore a new hybrid method where the potential is approximated by Chebyshev polynomials, which can be implemented in a few lines of Chebfun.

Numerical Analysis,Classical Analysis and ODEs

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: The performance of the classical WKB method in solving singularly perturbed linear differential equations has never been examined from the perspective of Structured Backward Error (SBE). The authors point out that for some important problems, the WKB method can give an exact solution to the same - structured problem that can be expressed in a finite number of terms. This analysis is very useful for evaluating the validity of the solutions provided by the WKB method.

Specifically, the main contributions of this paper include:
1. **Introducing Structured Backward Error Analysis**: Use structured backward error analysis to explain the success of the WKB method and the reasons for its failure.
2. **Proposing an Improved Algorithm**: Based on backward error analysis, a new iterative algorithm is proposed to improve the quality of the WKB solution.
3. **Exploring a Hybrid Method**: Combining the method of approximating the potential function with Chebyshev polynomials, a new hybrid WKB - Chebyshev method is proposed, and its implementation in Chebfun is demonstrated.

### Key Point Summary
- **Problem Background**: The WKB method is used to solve singularly perturbed linear differential equations of the form \(\varepsilon^2 y'' = Q(x)y\).
- **Innovation**: Analyze the WKB method from the perspective of structured backward error for the first time, and reveal the reasons why it can give exact solutions in some cases.
- **New Algorithms**: Propose an iterative algorithm and a hybrid method to improve the quality of the solutions of the WKB method.

### Mathematical Formula Representation
The mathematical formulas involved in the paper are represented in Markdown format as follows:

- The form of the WKB solution is:
\[
y_{\text{WKB}} = \exp\left(\frac{1}{\delta} S(x)\right)
\]
where \(S(x)\) is a possibly infinite series:
\[
S(x) = S_0(x) + \delta S_1(x) + \delta^2 S_2(x) + \cdots
\]

- For a differential equation of the form \(\varepsilon^2 y'' = Q(x)y\), the residual is defined as:
\[
r(y_n) = \varepsilon^2 y_n'' - Q(x) y_n
\]

- The differential equation satisfied by the approximate solution given by the WKB method is:
\[
\varepsilon^2 y'' = \tilde{Q}(x) y
\]
where \(\tilde{Q}(x)\) is the original potential \(Q(x)\) plus a correction term of \(O(\varepsilon^2)\):
\[
\tilde{Q}(x) = Q(x) + \varepsilon^2 Q_2(x)
\]

Through these analyses, the authors not only provide a deeper understanding of the WKB method, but also propose some practical new algorithms and methods to further improve the quality of the solutions of the WKB method.