Langer Modification, Quantization condition and Barrier Penetration in Quantum Mechanics

TL;DR

This paper introduces a uniform asymptotic approximation method for solving Schrödinger equations, improving accuracy over WKB, explaining Langer modification, and accurately calculating wave functions, transmission coefficients, and quantization conditions.

Contribution

The paper presents a new analytical approximation method that enhances accuracy and provides clear insights into Langer modification in quantum mechanics.

Findings

The method accurately reproduces wave functions for exactly solvable potentials.

It yields precise transmission coefficients for potential barriers.

The approach clarifies the origin of Langer modification in quantum problems.

Abstract

The WKB approximation plays an essential role in the development of quantum mechanics and various important results have been obtained from it. In this paper, we introduce another method, {\it the so-called uniform asymptotic approximations}, which is an analytical approximation method to calculate the wave functions of the Schr\"odinger-like equations, and is applicable to various problems, including cases with poles (singularities) and multiple turning points. An distinguished feature of the method is that in each order of the approximations the upper bounds of the errors are given explicitly. By properly choosing the freedom introduced in the method, the errors can be minimized, which significantly improves the accuracy of the calculations. A byproduct of the method is to provide a very clear explanation of the Langer modification encountered in the studies of the hydrogen atom and harmonic oscillator. To further test our method, we calculate (analytically) the wave functions for several exactly solvable potentials of the Schr\"odinger equation, and then obtain the transmission coefficients of particles over potential barriers, as well as the quantization conditions for bound states. We find that such obtained results agree with the exact ones extremely well. Possible applications of the method to other fields are also discussed.

Tables

Table 1. Table 1: Some exactly solvable potentials and the choices of q ​ ( x ) 𝑞 𝑥 q(x) . The exact energy eignvalue spectra for these solvable potentials can be found in dong_wave_2011 .

Potentials	$V\left(x\right)$	$q\left(x\right)$
Hydrogen	$-\frac{e^2}{x}+\frac{{\hslash }^{2}l\left(l+1\right)}{2m{x}^2}$	$-\frac{1}{4x^2}$
Harmonic oscillator	$\frac{1}{2}m{\omega }^2{x}^2+\frac{{\hslash }^{2}l\left(l+1\right)}{2m{x}^2}$	$-\frac{1}{4x^2}$
Morse potential	$v_0{e}^{-2\alpha x}+v_1{e}^{-\alpha x}$	$0$
Pöschl-Teller potential	$\frac{v_0}{{\cosh }^2\left(\alpha x\right)}$	$\frac{{\alpha }^2}{4{\cosh }^2\left(\alpha x\right)}$
Eckart potential	$\frac{v_0}{{\sinh }^2\left(\alpha x\right)}+\frac{v_1}{\tanh \left(\alpha x\right)}$	$-\frac{{\alpha }^2}{4sin{h}^2\left(\alpha x\right)}$

Equations

$$\displaystyle\frac{d^{2}\Psi(x)}{dx^{2}}+\frac{2m}{\hbar^{2}}\Big{(}E-V(x)\Big{)}\Psi(x)=0,$$

$$\Psi (x)\simeq \frac{\hslash }{\sqrt{2|p(x)|}}\exp \left[\pm \frac{i}{\hslash }{\int }^{x}p(x^{\prime })d{x}^{\prime }\right],$$

$$\mathcal{Q}\equiv {\hslash }^2\left|\frac{3p^{\prime 2}}{4p^4}-\frac{p^{\prime \prime }}{2p^3}\right|\ll 1.$$

$$V_C(r)=\frac{{\hslash }^{2}l(l+1)}{2m{r}^2},$$

$$\frac{d^2\Psi (x)}{d{x}^2}=\left\{g(x)+q(x)\right\}\Psi (x),$$

$$g(x)=f(x)(x-x_0),$$

$$\displaystyle\Psi(x)=\left(\frac{\xi(x)}{g(x)}\right)^{1/4}\Bigg{(}a_{0}{\rm Ai}(\xi)+b_{0}{\rm Bi}(\xi)\Bigg{)},$$

$$\displaystyle\mathscr{H}(\xi)\equiv\frac{5}{24|\xi|^{3/2}}-\int_{x_{0}}^{x}\Bigg{\{}\frac{q(\tilde{x})}{g(\tilde{x})}-\frac{5g^{\prime 2}(\tilde{x})}{16g^{3}(\tilde{x})}+\frac{g^{\prime\prime}(\tilde{x})}{4g^{2}(\tilde{x})}\Bigg{\}}$$

$$\times \sqrt{|g(\tilde{x})|}d\tilde{x}.$$

$$\Psi (x)$$

$$\mathcal{I}(\zeta )$$

$$(\mathcal{H},\mathcal{I})\to -{\int }^x\left(\frac{q}{g}-\frac{5g^{\prime 2}}{16g^3}+\frac{g^{\prime \prime }}{4g^2}\right)\sqrt{|g|}dx.$$

$$g(x)\sim a/x^2,$$

$$\left(\mathcal{H},\mathcal{I}\right)\to \underset{x\to 0}{\lim }\left(-\frac{\ln x}{4\sqrt{|a|}}-{\int }^x\frac{q}{\sqrt{|g|}}dx\right).$$

$$q(x)=-\frac{1}{4x^2},$$

$$g(x)=f(x)(x-x_1)(x-x_2),$$

$$\mathcal{I}(\zeta )$$

$$q(x)\simeq q_0+q_1(x-\mathrm{Re}(x_i)).$$

$$q_0={\frac{7f^{\prime 2}-6f{f}^{\prime \prime }}{32f^2}|}_{\mathrm{Re}(x_i)}.$$

$$q_0={\left(\frac{7g^{\prime \prime \prime 2}}{288g^{\prime \prime 2}}-\frac{g^{\prime \prime \prime \prime }}{32g^{\prime \prime }})\right|}_{x_m}.$$

$${\int }_{{\tilde{x}}_1}^{{\tilde{x}}_2}\sqrt{\frac{2m}{{\hslash }^2}\left(E-V(x)\right)}dx=\left(n+\frac{1}{2}\right)\pi ,$$

$${\int }_{x_b}^{{\tilde{x}}_0}\sqrt{\frac{2m}{{\hslash }^2}\left(E-V(x)\right)}dx=\left(n+\frac{3}{4}\right)\pi .$$

$$\frac{\pi {\zeta }_0^2}{2}={\int }_{x_1}^{x_2}\sqrt{-g(x)}=\left(n+\frac{1}{2}\right)\pi .$$

$${\int }_{x_0}^{x_b}\sqrt{-g(x)}=\left(n+\frac{3}{4}\right)\pi ,$$

$$V(x)=-\frac{e^2}{x}+\frac{{\hslash }^{2}l(l+1)}{2m{x}^2},$$

$$E_n=-\frac{m{e}^4}{2{\hslash }^2(n+l+1{)}^2}.$$

$$E_n=-\frac{m{e}^4}{2{\hslash }^2(n+1/2+\sqrt{l(l+1)}{)}^2},$$

$$q(x)\to -\frac{1}{4x^2}\mathrm{as}x\to 0,$$

$$q(x_m)={\left(\frac{7g^{\prime \prime \prime 2}}{288g^{\prime \prime 2}}-\frac{g^{\prime \prime \prime \prime }}{32g^{\prime \prime }})\right|}_{x_m}.$$

$$q(x)=-\frac{1}{4x^2}.$$

Full text

I Introduction

A fundamental interest in quantum mechanics (QM) is to derive various physical quantities from the wave function of the Schrödinger equation. Due to the complexity of the equation, it is extremely difficult to conduct analytical analyses, and various approximation methods have been proposed. Among them, the WKB approximation has played an essential role in the development of QM and been widely used in many fields of physics and chemistry berry_semiclassical_1972 ; friedrich_working_2004 ; price_semiclassical_2018 ; hyouguchi_divergencefree_2002 ; Karnakov_wkb ; froman_jwkb . In general, the one-dimensional Schrödinger equation reads,

[TABLE]

which describes a particle of mass $m$ moving with total energy $E$ in a potential $V(x)$. Then, using the WKB approximations, the wave function $\Psi(x)$ can be approximately written in the form

[TABLE]

where $p(x)=\sqrt{2m(E-V)}$ is the local momentum. The validity of the approximation is restricted to the regions where the WKB condition is fulfilled,

[TABLE]

In this way, one can treat the reduced Planck constant $\hbar$ as a small parameter and extend the leading-order solution (2) to high-orders.

However, it is well known that the above condition can be violated or not fulfilled completely in many cases berry_semiclassical_1972 . For example, the WKB condition is always violated around turning points ($p(x)=0$), at which both $\mathcal{Q}$ and the WKB wave function (2) diverge. In addition, the WKB condition is also violated around singular points (poles) of $p^{2}(x)$. For instance, for the radial Schrödinger equation, the effective potential $V(r)$ contains a centrifugal term

[TABLE]

which has a second-order pole at the origin $r=0$. Other typical potentials encountered in QM are presented in Table 1.

Note that it is exactly because of this second-order pole that the WKB approximation fails to give correct results for hydrogen atoms and harmonic oscillators young_wentzelbrillouinkramers_1930 . This problem was studied by Langer several decades ago, and shown that it can be cured if one replaces $l(l+1)$ by $(l+1/2)^{2}$ in $V_{C}(r)$ langer_connection_1937 . This modification now is considered as a standard ingredient of the WKB method in QM berry_semiclassical_1972 . However, a rigorous and logically consistent derivation of this modification is still lacking.

Another situation that could violate the WKB condition is around the extreme point of $p^{2}(x)$ if $\mathcal{Q}\simeq|p^{\prime\prime}/2p^{3}|\sim\mathcal{O}(1)$. This can arise in the bound states with a potential well or in particle scattering with a potential barrier. In both cases the results from the WKB approximation becomes invalid if $\mathcal{Q}$ is not small enough at the extreme points.

To overcome these problems, various (approximate) methods have been proposed dong_wave_2011 , including the complex Zwaan29 ; Langer34 ; LL97 ; Berry90 ; DDP97 and uniform WKB methods langer_connection_1937 ; MG53 ; Dingle56 ; Alvarez04 ; DU14 . In the complex method, one first generalizes the above one-dimensional problem into a complex $x$ plane, and then joins the two parts $x>a$ and $x<a$ by a path in the complex plane, which is sufficiently far from the turning point $x=a$, so the WKB solutions are valid along the path. With the seminal work of Delabaere, Dillinger and Pham DDP97 , the long-standing problem of the method, the ambiguities of the unique matching between the two parts berry_semiclassical_1972 , was finally resolved. As a result, a simple and rigorous justification of the Zinn-Justin quantization condition and its solution in terms of the “multi-instanton expansion” for the symmetric double well problem ZJ02 , were derived. On the other hand, the uniform (in the sense of being smooth across the turning points) WKB approximation method is a generalization of the basic WKB approximations by employing the technique of comparison equations Slavyanov96 ; FF96 . In this approach, the solution is mapped on to the one of a simpler equation which, however, has the same disposition of turning points as the original. Without bypassing the turning points, as done in the complex method DDP97 , the uniform WKB approximations also provide the wave function in the neighborhood of the turning points. Using this method, Álvarez was also able to derive the “multi-instanton expansion” for the eigenvalues of the symmetric double well by using a Langer-Cherry uniform asymptotic expansion Alvarez04 . An excellent pedagogical introduction of this method was recently provided by Dunne and Ünsal with both the double well and the Sine-Gordon potentials as concrete examples DU14 , in which the relations among the uniform WKB approach, multi-instanton and resurgence theory were also explored.

The purpose of this paper is to present another approximation method, the so-called uniform asymptotic approximation method, developed systematically by Olver olver1975 ; olver1997_book . Among other things, an essential difference of this method from the above ones is that explicit bounds are constructed for the error terms associated with the approximations, although, similar to the uniform WKB method, it is also based on the comparison equations Slavyanov96 ; FF96 . In particular, we shall calculate accurately the wave function of the Schrödinger equation (1) with singular potentials such as those given in Table 1, although the method is quite general, and in principle can be applied to any problem of the Schrödinger equation (1), including the cases with poles and multiple turning points. This method has been already shown to be powerful and robust when applied to calculations of the mode functions and primordial power spectra in a variety of slow-roll inflationary models Habib:2002yi ; Habib:2004kc ; Wang10 ; Zhu:2014wfa ; AKN16 ; Wu:2017joj , adiabatic regularization of the power spectra in various inflationary models Alinea:2015pza ; Alinea:2016qlf , inflation with nonlinear dispersion relations Zhu:2013fha ; Zhu:2013upa ; Zhu:2014wda ; Zhu:2014aea ; Zhu:2016srz ; Qiao:2018dpp ; Ding:2019nwu ; Qiao:2019hkz , and loop quantum cosmology Zhu:2015ata ; Zhu:2015owa ; Zhu:2015xsa ; Li:2018vzr , as well as studying the parametric resonance during inflation and reheating Zhu:2018smk . As mentioned above, the major advantage of the method is that the errors in each order of approximations can be estimated and the upper bounds of the errors are always known. In particular, for certain models, it was found that the errors are no larger than $0.15\%$ up to the third-order of approximations Zhu:2016srz .

In the application of this method to hydrogen atoms and harmonic oscillators, we provide a rigorous derivation of the Langer modification and expressions for the quantization conditions and the barrier transmission coefficients. Applications of our method to some well-known examples are also presented, in order to further test it.

II Uniform asymptotic approximation method

II.1 Wave functions in the uniform asymptotic approximation

Let us first write the standard form of (1) in the form olver1997_book ; olver1975 ,

[TABLE]

where $g(x)+q(x)=-p^{2}(x)/\hbar^{2}$. Note that, for any given $p(x)$, here we introduce two functions $g(x)$ and $q(x)$, and to fix them uniquely, we require that the errors in each order of approximations be minimized. This is one of the major ingredients of the method. Then, such defined $g(x)$ in general can have zero points $g(x_{i})=0$, which are called turning points in the uniform asymptotic approximation. Except such points, $g(x)$ may also have other types of transition points, such as poles and extreme points. According to the theory of the second-order ordinary differential equations olver1997_book ; olver1975 , the wave function $\Psi(x)$ sensitively depends on the number and nature of turning points, poles and extreme points. Analyzing the corresponding error control function around poles and extreme points provides the main guidance on how to determine $g(x)$ and $q(x)$ olver1997_book ; olver1975 . In addition, around each of the turning points $x_{i}$, we require $|q(x)|\ll|g(x)/(x-x_{i})|$, while away from them we require $|q(x)|\ll|g(x)|$ olver1997_book ; olver1975 .

At the turning points, the WKB condition (3) is violated, and the wave function (2) becomes invalid. Generally the turning points can have different types, depending on the nature of the zeros of $g(y)$. The simplest case is the single turning point which denotes the simple and single zero of $g(y)$. Beside this, $g(y)$ may also have a pair of turning points, which could be both real and single, double, or even complex conjugated. In the following we are going to construct the uniform wave function for above two cases individually.

II.1.1 For single turning point

Let us first consider the single turning point, say, denoted by $x_{0}$. We can write the function $g(x)$ in the form of

[TABLE]

where $f(x)$ is a regular function. In quantum mechanics, this turning point can be thought as the critical point that separates the classically allowed and forbidden regions of a quantum system if $q(x)=0$. When we have to chose $q(x)\neq 0$, it will produce a small shift on this point. Around the turning point $x_{0}$, the wave function $\Psi(x)$ can be written as olver1997_book

[TABLE]

where ${\rm Ai}(\xi)$ and ${\rm Bi(\xi)}$ are the Airy functions, $a_{0}$ and $b_{0}$ are two integration constants, and $\xi(x)$ is chosen to be a monotonic function of $x$, which has the same sign of $g(x)$ with relations $\sqrt{|\xi|}d\xi=\sqrt{|g(x)|}dx$ and $\xi(x_{0})=0$. The errors are controlled by the error control function,

[TABLE]

II.1.2 A pair of turning points

For a pair of turning points $x_{1}$ and $x_{2}$, they could be: (1) both real and different $x_{1}\neq x_{2}$; (2) both real but equal $x_{1}=x_{2}$; or (3) complex conjugate $x_{1}=x_{2}^{*}$. In each case, between these points, $g(y)$ usually has one extreme.

If this extreme is a minimal point of $g(y)$, we can construct the wave function $\Psi(x)$ in terms of the parabolic cylinder functions $U(-\zeta_{0}^{2}/2,\sqrt{2}\zeta)$ and $\bar{U}(-\zeta_{0}^{2}/2\sqrt{2}\zeta)$ as olver1975 ,

[TABLE]

If the extreme is a maximal point, the wave function can be also constructed in terms of parabolic cylinder functions, but now in terms of $W(\zeta_{0}^{2}/2,\sqrt{2}\zeta)$ and $W(\zeta_{0}^{2}/2,-\sqrt{2}\zeta)$,

[TABLE]

In both cases, the variable $\zeta(x)$ is a monotonic function of $x$ and defined by $\sqrt{|\zeta^{2}-\zeta_{0}^{2}|}d\zeta=\sqrt{|g(x)|}dx$ with $\zeta(x_{1})=-|\zeta_{0}|$, $\zeta(x_{2})=|\zeta_{0}|$, and $\zeta_{0}^{2}=\pm(2/\pi)|\int_{x_{1}}^{x_{2}}\sqrt{|g(x)|}dx|$, where $``+"$ $(``-")$ corresponds to the case in which the two turning points $x_{1}$ and $x_{2}$ are real (complex conjugate). The associated error control function of the above two wave functions is,

[TABLE]

The wave functions given in (7), (9) and (10) are valid in the neighborhoods of the turning points $x_{i}$. The extension of them beyond these points crucially depends on the behaviors of the corresponding error control functions in the extended regions. In the following let us consider it for the case with a second-order pole, as shown in (4).

II.2 Second-order pole and Langer modification

As mentioned above, for the radial Schrödinger equation, the effective potential $V_{C}(r)$ given by (4) contains a second-order pole at the origin, at which we have $(\xi,\;\zeta)\to\pm\infty$, and

[TABLE]

Since near the second-order pole, $g(x)$ has the asymptotic behavior

[TABLE]

we find that the error control function takes the limit

[TABLE]

In order to make our analytic solutions in (7), (9) and (10) also valid at the origin, we must require that the error control function given by (14) are finite at the pole $x=0$. Using the asymptotic behavior of $g(x)$ given in (13), it is straightforward to show that the right-hand side of (14) can be made finite at the pole as long as Zhu:2013upa

[TABLE]

which is nothing but exactly the Langer modification. Thus, the latter is simply the result of imposing the condition that the error control function be finite at the pole in the uniform asymptotic approximation method.

II.3 Extreme point and the elimination of the error term

The extreme points of $g(x)$ in general arise from quantum system with a potential well or barrier in the region between two turning points. These extreme points are the same as the bottom of the well or the top of the potential if one chooses $q(x)=0$. As mentioned above, the existence of the extreme points will make the WKB approximation invalid if $\mathcal{Q}$ is not small enough at the extreme. Such cases rise when the potential wells or barriers are sharply peaked, for which $g(x)$ has two coalescent turning points $x_{1}$ and $x_{2}$. To be more specific, let us write $g(x)$ in the form

[TABLE]

where $f(x)$ is a finite and regular function with $f(x_{i})\not=0$. Then, we expect the dominant contribution to the integral of (II.1.2) arises from the lower limit. Therefore, we can formally expand the error control function $\mathscr{I}(\zeta)$ about the turning points and find 111 It should be noted that the integration of the form, $I=\int{dxg(x)\exp\left[if(x)\right]}$, by using the method of stationary phase, is well-established, see, for example, Miller70 , where there are two roots, $x_{1}$ and $x_{2}$, of the equation $f^{\prime}(x)=0$. When $x_{2}-x_{1}$ is small, the method leads to the solution, $I=g_{0}\exp(if_{0})2\pi\left|2/f^{\prime\prime\prime}_{0}\right|^{1/3}Ai[f_{0}(2/f^{\prime\prime\prime}_{0})^{1/3}]$, where $x_{0}$ is the value of $x$ for which $f^{\prime\prime}(x_{0})=0$, $g_{0}\equiv g(x_{0})$ and so on. However, to our current purpose, we find that the expression of (17) is more suitable.,

[TABLE]

where ${\rm Re}(x_{i})$ denotes the real part of the turning point $x_{i}$. Note that in deriving the above expression we had ignored the small corrections. We also note that when the turning points $x_{1}$ and $x_{2}$ are close to each other, the dominant contributions to the error control function come from the $\ln|x_{2}-x_{1}|$ term, which becomes divergent at $x_{1}=x_{2}$. It is somehow surprising that such dominant contributions seemingly had never been noticed before, and later they will play an essential role in determining the extension of the wave functions to the regions near the extreme point.

With the knowledge of the error control function (17), we are now in a position to eliminate the dominant error term in (17) by properly choosing $q(x)$ in the second term of (17). To achieve this, we expand $q(x)$ at one of the turning points 222 It should be noted that in general the expansion should be carried out in terms of $x_{i}$ in the complex $x$ plane LL97 . However, we find that for the analysis of the error control function $\mathscr{I}(\zeta)$ defined by (17) near the turning point, the expansion alone the real axis is sufficient. In particular, it is involved only with the choice of the zeroth-order term $q_{0}$, as can be seen from (19) and (20).,

[TABLE]

Then the elimination of the first error term in (17) requires

[TABLE]

This represents one of the important conditions for the choice of the function $q(x)$. Additional requirements include that $|q(x)|$ be negligible in comparison with $|g(x)|$ in the regions that is away from the extreme and turning points. Then, we expect that the right hand side of (19) does not contain $q_{0}$ and should be independent of the nature of the two turning points $x_{1}$ and $x_{2}$. For this requirement, one can relate $q_{0}$ to the derivatives of the function $g(x)$ at the extreme point $x_{m}$ via the relation,

[TABLE]

III Improved quantization conditions and potential barrier transmission coefficients

With the above considerations, now we are at the position to generalize the wave functions near poles and extreme points. These wave functions then can be utilized to derive the quantization conditions for bound states or the quantum transmission coefficients of a particle through a potential barrier. In this section, we are going to discuss quantization conditions and potential barrier transmission individually by using the quantum mechanical wave functions derived in the above section with the corresponding $q(x)$ calculated for specific potentials.

III.1 Improved Quantization Condition

One of the most important topics in quantum mechanics is to determine the bound states and energy eigenvalue spectra for a quantum mechanical system. As mentioned in the introduction, the conventional WKB method fails to describe the correct behaviors of wave functions around the poles, extreme points, and turning points of $E-V(x)$ in the Schrödinger equation (1). Especially at the turning point, the WKB wave function (2) becomes divergent. One approach to avoid this singularity is the complex method in which the WKB approximation has been extended to the complex $x$ plane. In this way, the two WKB wave functions in the classically allowed and forbidden regions can be connected without the knowledge of the wave functions around the turning points. With this procedure, when a particle is trapped in a potential well $V(x)$, the corresponding energy eigenvalue spectra can be obtained by using the conventional WKB quantization condition,

[TABLE]

where $\tilde{x}_{1}$ and $\tilde{x}_{2}$ are two zeros of $E-V(x)$ (or classical turning points of (1)) and $n$ denotes the quantum number. It is worth noting that $\tilde{x}_{1}$ and $\tilde{x}_{2}$ is slightly different from the turning points $x_{1}$ and $x_{2}$ of $g(y)$ with $q(x)\neq 0$. For one single turning point system, where a particle is trapped in a classically allowed region with a boundary at $x_{b}$, the WKB quantization gives,

[TABLE]

However, in the uniform asymptotic approximations, if $g(x)$ has two turning points $x_{1}$ and $x_{2}$, the wave function (9) leads to the quantization condition,

[TABLE]

When $g(x)$ has only one turning point $x_{0}$, the wave function (7) leads to the quantization condition,

[TABLE]

where $x_{b}$ is the boundary of the classically allowed region. The above quantization condition is different from that derived in the WKB approximation except for the case with $q(x)=0$. For the quantum mechanical system that $q(x)$ has to be nonzero, the above quantization condition can provide a significant improvement on the calculations of the energy eigenvalue spectra for a quantum mechanical system. In the following, for applications of the above quantization condition, we consider several representative exactly solvable systems, including hydrogen atoms, harmonic oscillators in D dimensions, Morse potential, Pöschl-Teller potential, and Eckart potential, as given in Table. 1.

III.1.1 Hydrogen atoms

For hydrogen atoms, the electron moves in the atom with a Coulomb potential and a centrifugal term dong_wave_2011 ,

[TABLE]

where $e$ denotes the electron charge. The Schrödinger equation (1) of the hydrogen atom can be solved exactly, and from the exact solution one can easily obtain the exact energy eigenvalue spectra as dong_wave_2011

[TABLE]

If we apply the WKB quantization condition (21), the WKB energy eigenvalue spectra read

[TABLE]

which does not match the exact result. A standard procedure to solve this problem in quantum mechanics can be implemented using the Langer modifciation by simply replacing $l(l+1)$ by $(l+1/2)^{2}$ in the Schrödinger equation langer_connection_1937 .

In the uniform asymptotic approximation, the above problem can be resolved by properly choosing $q(x)$ to eliminate the divergent errors of the approximate solution arising from the existence of the second-order pole and the extreme point. It can be seen that the Schrödinger equation (1) for the hydrogen atom has a second order pole, one extreme point, and two turning points. Then, according to the analysis in the above section, $q(x)$ has to be chosen such that

[TABLE]

To meet these two conditions, $q(x)$ must be

[TABLE]

With this choice, the two turning points $x_{1}$ and $x_{2}$ can be solved from $g(x)=0$, yielding

[TABLE]

Substituting them into the improved quantization condition (23), one immediately obtains

[TABLE]

which precisely recovers the exact result. The eigenstate for each quantum number then is given by the wave function in (9) with $q(x)$ given in (30), which is uniformly valid in the whole region $x\in(0,+\infty)$. It is worth noting that the above analysis can be easily extended to the hydrogen-like atoms.

III.1.2 Harmonic Oscillators

The potential for the harmonic oscillator in D dimensions is given by dong_wave_2011

[TABLE]

The exact solution of the Schrödinger equation (1) with the above potential leads to dong_wave_2011

[TABLE]

If we perform the WKB quantization condition (21), the WKB energy eigenvalue spectra read

[TABLE]

Similar to the hydrogen atom, it does not match the exact result.

To do the analysis in the uniform asymptotic approximation, we observe that Schrödinger equation (1) for the harmonic oscillator in D dimensions has one second-order pole, one extreme point, and two turning points. This is very similar to the case in the hydrogen atom. Then the choice of $q(x)$ has to meet the two conditions in (28) and (29), which leads to

[TABLE]

With this choice, the improved quantization condition (23) then gives

[TABLE]

which again recovers the exact result. Similarly the eigenstate for each quantum number is given by the wave function in (9) which is valid in the whole region of interest.

III.1.3 Morse potential

The Morse potential is given by dong_wave_2011

[TABLE]

where $v_{0}$, $v_{1}$, and $\alpha$ are three parameters in the potential. The exact solution of the Schrödinger equation (1) with the Morse potential leads to the following energy eigenvalue spectra dong_wave_2011 ,

[TABLE]

If we perform the WKB quantization condition (21), the WKB energy eigenvalue spectra takes exactly the same form as the exact result. This is in contrast to the cases for the hydrogen atom and the harmonic oscillator in D dimensions.

The analysis in the uniform asymptotic approximation can give the same result as well. Observing that the Morse potential has one extreme point between two turning points, application of the criterion (20) immediately leads to

[TABLE]

For this case once we choose $q(x)=0$, and then the improved quantization condition (23) recovers the WKB quantization condition (21).

III.1.4 Pöschl-Teller potential

The Pöschl-Teller potential is given by dong_wave_2011

[TABLE]

where $v_{0}<0$ and $\alpha$ are two parameters describing the potential. The exact solution of the Schrödinger equation (1) with the Pöschl-Teller potential leads to the exact energy eigenvalue spectra dong_wave_2011 ,

[TABLE]

If we perform the WKB quantization condition (21), the WKB energy eigenvalue spectra read

[TABLE]

Obviously it does not match the exact result.

Considering that the Pöschl-Teller potential only contains one extreme point between the two turning points, we can use the the criterion (20) to determine the choice of $q(x)$. An immediately calculation shows that

[TABLE]

As a result, a proper choice for $q(x)$ is

[TABLE]

With this choice, the improved quantization condition (23) gives

[TABLE]

which again is precisely the same as the exact result. Similarly the eigenstate for each quantum number is given by the wave function in (9) which is valid in the whole region of interest.

III.1.5 Eckart potential

The Eckart potential takes the form dong_wave_2011

[TABLE]

where $v_{0}$, $v_{1}$, and $\alpha$ are three parameters describing the potential. The exact solution of the Schrödinger equation (1) with the Eckart potential leads to the exact energy eigenvalue spectra, which reads dong_wave_2011

[TABLE]

If we perform the WKB quantization condition (21), the WKB energy eigenvalue spectra read

[TABLE]

We observe that it does not match the exact result.

Since the Eckart potential contains one second-order pole at $x=0$ and one extreme point between the two turning points, we need to use the two criterions (15) and (20) to determine the choice of $q(x)$. A straightforward calculation shows that a proper choice for $q(x)$ to meet both conditions is

[TABLE]

With this choice, the improved quantization condition (23) gives

[TABLE]

which again is precisely the same as the exact result. Similarly the eigenstate for each quantum number is given by the wave function in (9), which is valid in the whole region of interest.

In summary, it can be seen that from the above five examples that, while the standard WKB quantization condition fails to predict correct energy eigenvalues $E_{n}$ except for the Morse potential, the quantization condition (23) yields precisely the exact results for all these potentials. We note that the exact energy eigenvalue spectra for these potentials can also be obtained from alternative approaches, for examples, by using the supersymmetry quantum mechanics cooper_supersymmetry_2001 or the proper quantization rules developed in qiang_proper_2010 ; dong_wave_2011 ; serrano_qiang_2010 . In contrast to these approaches for solvable potentials, we would like to mention here that the improved quantization condition (23) is not only applicable to the solvable potentials, but also valid for more general potentials that do not admit exact solutions. More importantly, the bound states for a general potential can also be obtained from the wave functions (9) which are derived from our approach.

III.2 Potential Barrier Transmission

In the above subsection, we have studied the bound state and energy eigenvalue problems when a particle moves in a potential well. In this subsection, let us move on to the scattering problem for a particle passing through a potential barrier. As mentioned in the above, the standard WKB approximation cannot give the correct behaviors of the wave function about the turning points. When the energy of the particle is below the peak of the barrier, the Schrödinger equation (1) normally has two real turning points $\tilde{x}_{1}$ and $\tilde{x}_{2}$. For this case, the complex method of the WKB approximation leads to an approximate transmission coefficient,

[TABLE]

where $\tilde{x}_{1}$ and $\tilde{x}_{2}$ are two turning points defined by $E-V(x)=0$.

The transmission coefficient in the uniform asymptotic approximation can be directly obtained from the wave function given in (10) with known $q(x)$. This treatment is universal and does not depend on whether the energy of the particle is above or below the peak of the barrier. Using the wave function (10), the improved transmission coefficient reads

[TABLE]

We note that $\zeta_{0}^{2}$ is positive when $x_{1}$ and $x_{2}$ are real and negative when $x_{1}$ and $x_{2}$ are complex conjugated. It is worth emphasizing again that (54) is valid in either case in which the energy $E$ is above or below the peak of the potential barriers.

To calculate the transmission coefficient (54) for a specific form of the potential, we first need to determine the choice of $q(x)$. In this subsection, instead of calculating various potentials in details, let us consider the Pöschl-Teller potential barrier with a positive $v_{0}$ as an example. We note that this potential barrier has been applied to the calculation of the black hole quasi-normal modes berti_quasinormal_2009 ; konoplya_quasinormal_2011 and the recent studying of the primordial perturbation during quantum bounce in the loop quantum cosmology Wu:2018sbr ; Wu:2018sbr1 ; Wu:2018sbr2 ; Wu:2018sbr3 . The choice of $q(x)$ is the same as in the case for the Pöschl-Teller potential well. Then, the transmission coefficient (54) reads

[TABLE]

while the (standard) WKB approximation gives gottfried

[TABLE]

It is worth noting that the choice of $q(x)$ is only possible when $8mv_{0}/\hbar^{2}>\alpha^{2}$ for positive $v_{0}$. When $8mv_{0}/\hbar^{2}<\alpha^{2}$, as mentioned above, the choice of $q(x)$ changes the properties of turning points significantly, which makes the approximation not applicable. In Fig. 1, we present the transmission coefficient (55), the WKB transmission coefficient (56) and the exact result, from which one can see that the transmission coefficient (55) fits the exact result extremely well.

IV Summary and Outlook

In this paper, a new analytical approximation method to solve the Schrödinger equation, the so-called uniform asymptotic approximation method, developed systematically by Olver olver1975 ; olver1997_book , is presented and applied to several cases of interest in quantum mechanics. This new method makes use of the technique of uniform approximations to map the original Schrödinger equation to a simpler equation for which an approximate solution can be found analytically. One of the major advantage of the method is that the errors in each order of approximations can be estimated and their upper bounds are always known explicitly. In particular, for certain models, it was found that the errors are no larger than $0.15\%$ up to the third-order of approximations Zhu:2016srz .

To illustrate the above, in this paper we have given the general solutions from the new method in the region where one single turning point or two turning points are encountered. One of the advantages of the new method is that by carefully choosing the function $q(x)$, the errors of the approximate solutions can be well-controlled. In general, the function $q(x)$ has to be chosen to satisfy two conditions: (a) near the turning points $x_{i}$, $|q(x)|\ll|g(x)/(x-x_{i})|$; and (b) away from the turning points, $|q(x)|\ll|g(x)|$. Moreover, when one of the the boundaries of the interesting region, for example $x=0$, becomes a second-order pole, $q(x)$ has to be fixed correspondingly to $-\frac{1}{4x^{2}}$ in order to make the error control function finite at the pole. One of the results from this particular choice of the function $q(x)$ is that the quantization conditions for bound states from the exact solutions can be recovered from our approximate solutions for several quantum mechanical systems. Moreover, when the form of the error control function can be analytically expressed in terms of $q(x)$ and $g(x)$, the function $q(x)$ can also be chosen to minimize the errors at the tuning point or at the extreme point. With the new method, we have also derived the transmission coefficients for particle scattering over a potential barrier, which significantly improve the accuracy of the results obtained from the standard WKB approximations.

It must be noted that in this article we have mainly compared the uniform asymptotic approximation method introduced in this paper with the conventional (standard) WKB approximations. Other methods, such as the complex Zwaan29 ; Langer34 ; LL97 ; Berry90 ; DDP97 and uniform WKB langer_connection_1937 ; MG53 ; Dingle56 ; Alvarez04 ; DU14 , have also achieved great successes, and high (even exponential) precisions can be obtained. However, comparing the uniform asymptotic approximation method with them directly is not an easy task. In particular, the precisions of these methods, including the asymptotic approximation method, depend on the models and properties of the potentials $V(x)$ appearing in (1), and for a given model, one method might achieve much better precisions than others, but in other models the situation can be completely different. In the uniform WKB approximations and the current method, the analysis and precisions of a given problem also depend on the specific choice of the comparison equations RY64 ; Slavyanov96 ; FF96 . Considering the scope of this paper, we shall leave this issue to another occasion.

The new analytical approximation method presented in this paper is also useful to future investigations on the high-order approximations, extension to the nonlinear Schrödinger equation, and applications to various quantum systems with general potentials. Our formulas are general and simple to use, and can also be applied to other interesting cases, for example, Hawking radiation and quasi-normal modes of black holes konoplya_quasinormal_2011 ; berti_quasinormal_2009 , primordial perturbations during inflation and reheating bassett_inflation_2006 , and Schwinger vacuum pair productions due to laser pulses dumlu_stokes_2010 .

Acknowledgements

This work is partially supported by the National Natural Science Foundation of China under the Grants Nos. 11675143, 11675145, and 11975203; the Zhejiang Provincial Natural Science Foundation of China under Grant No. LY20A050002; and the Fundamental Research Funds for the Provincial Universities of Zhejiang, China under Grants No. RF-A2019015. B.-F. L is partially supported by NSF grant PHY-1454832.