# WKB approximation

﻿
WKB approximation

In physics, the WKB (Wentzel–Kramers–Brillouin) approximation, also known as WKBJ (Wentzel–Kramers–Brillouin–Jeffreys) approximation, is the most familiar example of a semiclassical calculation in quantum mechanics in which the wavefunction is recast as an exponential function, semiclassically expanded, and then either the amplitude or the phase is taken to be slowly changing.

Brief history

This method is named after physicists Wentzel, Kramers, and Brillouin, who all developed it in 1926. In 1923, mathematician Harold Jeffreys had developed a general method of approximating linear, second-order differential equations, which includes the Schrödinger equation. But since the Schrödinger equation was developed two years later, and Wentzel, Kramers, and Brillouin were apparently unaware of this earlier work, Jeffreys is often neglected credit. Early texts in quantum mechanics contain any number of combinations of their initials, including WBK, BWK, WKBJ and BWKJ.

Earlier references to the method are: Carlini in 1817, Liouville in 1837, Green in 1837, Rayleigh in 1912 and Gans in 1915. Liouville and Green may be called the founders of the method, in 1837.
The important contribution of Wentzel, Kramers, Brillouin and Jeffreys to the method was the inclusion of the treatment of turning points, connecting the evanescent and oscillatory solutions at either side of the turning point. For example, this may occur in the Schrödinger equation, due to a potential energy hill.

WKB method

Generally, WKB theory is a method for approximating the solution of a differential equation whose highest derivative is multiplied by a small parameter ε. The method of approximation is as follows:

For a differential equation

:$epsilon frac\left\{d^ny\right\}\left\{dx^n\right\} + a\left(x\right)frac\left\{d^\left\{n-1\right\}y\right\}\left\{dx^\left\{n-1 + cdots + k\left(x\right)frac\left\{dy\right\}\left\{dx\right\} + m\left(x\right)y= 0$

assume a solution of the form of an asymptotic series expansion

:$y\left(x\right) sim expleft \left[frac\left\{1\right\}\left\{delta\right\}sum_\left\{n=0\right\}^\left\{infty\right\}delta^nS_n\left(x\right) ight\right]$

In the limit $delta ightarrow 0$. Substitution of the above ansatz into the differential equation and canceling out the exponential terms allows one to solve for an arbitrary number of terms $S_n\left(x\right)$ in the expansion. WKB Theory is a special case of Multiple Scale Analysis.

An example

Consider the second-order homogeneous linear differential equation

:$epsilon^2 frac\left\{d^2 y\right\}\left\{dx^2\right\} = Q\left(x\right) y$

where $Q\left(x\right) eq 0$. Plugging in

:$y\left(x\right) = expleft \left[frac\left\{1\right\}\left\{delta\right\}sum_\left\{n=0\right\}^\left\{infty\right\}delta^nS_n\left(x\right) ight\right]$

results in the equation

:$epsilon^2left \left[frac\left\{1\right\}\left\{delta^2\right\}left\left(sum_\left\{n=0\right\}^\left\{infty\right\}delta^nS_n\text{'} ight\right)^2 + frac\left\{1\right\}\left\{delta\right\}sum_\left\{n=0\right\}^\left\{infty\right\}delta^nS_n" ight\right] = Q\left(x\right)$

To leading order, (assuming, for the moment, the series will be asymptotically consistent) the above can be approximated as

:$frac\left\{epsilon^2\right\}\left\{delta^2\right\}S_0\text{'}^2 + frac\left\{2epsilon^2\right\}\left\{delta\right\}S_0\text{'}S_1\text{'} + frac\left\{epsilon^2\right\}\left\{delta\right\}S_0" = Q\left(x\right)$

In the limit $delta ightarrow 0$, the dominant balance is given by

:$frac\left\{epsilon^2\right\}\left\{delta^2\right\}S_0\text{'}^2 sim Q\left(x\right)$

So δ is proportional to ε. Setting them equal and comparing powers renders

:$epsilon^0: ; ; ; S_0\text{'}^2 = Q\left(x\right)$

Which can be recognized as the Eikonal equation, with solution

:$S_0\left(x\right) = pm int_\left\{x_0\right\}^\left\{x\right\}sqrt\left\{Q\left(t\right)\right\},dt$

Looking at first-order powers of $epsilon$ gives

:$epsilon^1: ; ; ; 2S_0\text{'}S_1\text{'} + S_0" = 0$

Which is the unidimensional transport equation, which has the solution

:$S_1\left(x\right) = -frac\left\{1\right\}\left\{4\right\}logleft\left(Q\left(x\right) ight\right) + k_1.,$

And $k_1$ is an arbitrary constant. We now have a pair of approximations to the system (a pair because $S_0$ can take two signs); the first-order WKB-approximation will be a linear combination of the two:

:$y\left(x\right) approx c_1Q^\left\{-frac\left\{1\right\}\left\{4\left(x\right)expleft \left[frac\left\{1\right\}\left\{epsilon\right\}int_\left\{x_0\right\}^xsqrt\left\{Q\left(t\right)\right\}dt ight\right] + c_2Q^\left\{-frac\left\{1\right\}\left\{4\left(x\right)expleft \left[-frac\left\{1\right\}\left\{epsilon\right\}int_\left\{x_0\right\}^xsqrt\left\{Q\left(t\right)\right\}dt ight\right]$

Higher-order terms can be obtained by looking at equations for higher powers of ε. Explicitly

:$2S_0\text{'}S_n\text{'} + S"_\left\{n-1\right\} + sum_\left\{j=1\right\}^\left\{n-1\right\}S\text{'}_jS\text{'}_\left\{n-j\right\} = 0$

for $n>2$. This example comes from Bender and Orszag's textbook (see references).

Application to Schrödinger equation

The one dimensional, time-independent Schrödinger equation is:$-frac\left\{hbar^2\right\}\left\{2m\right\} frac\left\{mathrm\left\{d\right\}^2\right\}\left\{mathrm\left\{d\right\}x^2\right\} Psi\left(x\right) + V\left(x\right) Psi\left(x\right) = E Psi\left(x\right)$,

which can be rewritten as

:$frac\left\{mathrm\left\{d\right\}^2\right\}\left\{mathrm\left\{d\right\}x^2\right\} Psi\left(x\right) = frac\left\{2m\right\}\left\{hbar^2\right\} left\left( V\left(x\right) - E ight\right) Psi\left(x\right)$.

The wavefunction can be rewritten as the exponential of another function Φ (which is closely related to the action):

:$Psi\left(x\right) = e^\left\{Phi\left(x\right)\right\}, !$

so that

:$Phi"\left(x\right) + left \left[Phi\text{'}\left(x\right) ight\right] ^2 = frac\left\{2m\right\}\left\{hbar^2\right\} left\left( V\left(x\right) - E ight\right),$

where Φ' indicates the derivative of Φ with respect to "x". The derivative $Phi\text{'}\left(x\right)$ can be separated into real and imaginary parts by introducing the real functions "A" and "B":

:$Phi\text{'}\left(x\right) = A\left(x\right) + i B\left(x\right). ;$

The amplitude of the wavefunction is then $e^\left\{A\left(x\right)\right\}$ while the phase is $B\left(x\right)$. The Schrödinger equation implies that these functions must satisfy:

:$A\text{'}\left(x\right) + A\left(x\right)^2 - B\left(x\right)^2 = frac\left\{2m\right\}\left\{hbar^2\right\} left\left( V\left(x\right) - E ight\right)$

and therefore, since the right hand side of the differential equation for Φ is real,

:$B\text{'}\left(x\right) + 2 A\left(x\right) B\left(x\right) = 0. ;$

Next, the semiclassical approximation is invoked. This means that each function is expanded as a power series in $hbar$. From the equations it can be seen that the power series must start with at least an order of $hbar^\left\{-1\right\}$ to satisfy the real part of the equation. In order to achieve a good classical limit, it is necessary to start with as high a power of Planck's constant as possible.

:$A\left(x\right) = frac\left\{1\right\}\left\{hbar\right\} sum_\left\{n=0\right\}^infty hbar^n A_n\left(x\right)$

:$B\left(x\right) = frac\left\{1\right\}\left\{hbar\right\} sum_\left\{n=0\right\}^infty hbar^n B_n\left(x\right)$

To first order in this expansion, the conditions on "A" and "B" can be written.

:$A_0\left(x\right)^2 - B_0\left(x\right)^2 = 2m left\left( V\left(x\right) - E ight\right)$

:$A_0\left(x\right) B_0\left(x\right) = 0 ;$

If the amplitude varies sufficiently slowly as compared to the phase ($A_0\left(x\right) = 0$), it follows that

:$B_0\left(x\right) = pm sqrt\left\{ 2m left\left( E - V\left(x\right) ight\right) \right\},$

which is only valid when the total energy is greater than the potential energy, as is always the case in classical motion. After the same procedure on the next order of the expansion it follows that

:$Psi\left(x\right) approx C_0 frac\left\{ e^\left\{i int mathrm\left\{d\right\}x sqrt\left\{frac\left\{2m\right\}\left\{hbar^2\right\} left\left( E - V\left(x\right) ight\right)\right\} + heta\right\} \right\}\left\{sqrt \left[4\right] \left\{frac\left\{2m\right\}\left\{hbar^2\right\} left\left( E - V\left(x\right) ight\right)$

On the other hand, if it is the phase varies that varies slowly (as compared to the amplitude), ($B_0\left(x\right) = 0$) then

:$A_0\left(x\right) = pm sqrt\left\{ 2m left\left( V\left(x\right) - E ight\right) \right\}$

which is only valid when the potential energy is greater than the total energy (the regime in which quantum tunneling occurs). Grinding out the next order of the expansion yields

:$Psi\left(x\right) approx frac\left\{ C_\left\{+\right\} e^\left\{+int mathrm\left\{d\right\}x sqrt\left\{frac\left\{2m\right\}\left\{hbar^2\right\} left\left( V\left(x\right) - E ight\right) + C_\left\{-\right\} e^\left\{-int mathrm\left\{d\right\}x sqrt\left\{frac\left\{2m\right\}\left\{hbar^2\right\} left\left( V\left(x\right) - E ight\right)\right\}\left\{sqrt \left[4\right] \left\{frac\left\{2m\right\}\left\{hbar^2\right\} left\left( V\left(x\right) - E ight\right).$

It is apparent from the denominator, that both of these approximate solutions 'blow up' near the classical turning point where $E = V\left(x\right)$ and cannot be valid. These are the approximate solutions away from the potential hill and beneath the potential hill. Away from the potential hill, the particle acts similarly to a free wave—the phase is oscillating. Beneath the potential hill, the particle undergoes exponential changes in amplitude.

To complete the derivation, the approximate solutions must be found everywhere and their coefficients matched to make a global approximate solution. The approximate solution near the classical turning points $E=V\left(x\right)$ is yet to be found.

For a classical turning point $x_1$ and close to $E=V\left(x_1\right)$, $frac\left\{2m\right\}\left\{hbar^2\right\}left\left(V\left(x\right)-E ight\right)$ can be expanded in a power series.

:$frac\left\{2m\right\}\left\{hbar^2\right\}left\left(V\left(x\right)-E ight\right) = U_1 \left(x - x_1\right) + U_2 \left(x - x_1\right)^2 + cdots$

To first order, one finds

:$frac\left\{mathrm\left\{d\right\}^2\right\}\left\{mathrm\left\{d\right\}x^2\right\} Psi\left(x\right) = U_1 \left(x - x_1\right) Psi\left(x\right).$

This differential equation is known as the Airy equation, and the solution may be written in terms of Airy functions.:$Psi\left(x\right) = C_A extrm\left\{Ai\right\}left\left( sqrt \left[3\right] \left\{U_1\right\} \left(x - x_1\right) ight\right) + C_B extrm\left\{Bi\right\}left\left( sqrt \left[3\right] \left\{U_1\right\} \left(x - x_1\right) ight\right).$

This solution should connect the far away and beneath solutions. Given the 2 coefficients on one side of the classical turning point, the 2 coefficients on the other side of the classical turning point can be determined by using this local solution to connect them. Thus, a relationship between $C_0, heta$ and $C_\left\{+\right\},C_\left\{-\right\}$ can be found.

Fortunately the Airy functions will asymptote into sine, cosine and exponential functions in the proper limits. The relationship can be found to be as follows (often referred to as "connection formulas"):

:

Now the global (approximate) solutions can be constructed.

Precision of the asymptotic series

The asymptotic series for $y\left(x\right)$ is usually a divergent series whose general term $delta ^n S_n\left(x\right)$ starts to increase after a certain value $n=n_max$. Therefore the smallest error achieved by the WKB method is at best of the order of the last included term. For the equation:$epsilon^2 frac\left\{d^2 y\right\}\left\{dx^2\right\} = Q\left(x\right) y$with $Q\left(x\right)<0$ an analytic function, the value $n_max$ and the magnitude of the last term can be estimated as follows (see Winitzki 2005),

:$n_max approx 2epsilon^\left\{-1\right\} left| int_\left\{x_0\right\}^\left\{x_\left\{* dzsqrt\left\{-Q\left(z\right)\right\} ight| ,$:$delta^\left\{n_max\right\}S_\left\{n_max\right\}\left(x_0\right) approx sqrt\left\{frac\left\{2pi\right\}\left\{n_max exp \left[-n_max\right] ,$

where $x_0$ is the point at which $y\left(x_0\right)$ needs to be evaluated and $x_\left\{*\right\}$ is the (complex) turning point where $Q\left(x_\left\{*\right\}\right)=0$, closest to $x=x_0$. The number $n_max$ can be interpreted as the number of oscillations between $x_0$ and the closest turning point. If $epsilon^\left\{-1\right\}Q\left(x\right)$ is a slowly-changing function, :$epsilonleft| frac\left\{dQ\right\}\left\{dx\right\} ight| ll Q^2 ,$ the number $n_max$ will be large, and the minimum error of the asymptotic series will be exponentially small.

ee also

* Airy Function
* Langer correction
* Method of steepest descent / Laplace Method
* Perturbation methods
* Quantum tunneling
* Old quantum theory

References

Modern references

*cite book | author=Razavy, Moshen | title=Quantum Theory of Tunneling | publisher=World Scientific | year=2003 | id=ISBN 981-238-019-1
*cite book | author=Griffiths, David J. | title=Introduction to Quantum Mechanics (2nd ed.) | publisher=Prentice Hall |year=2004 |id=ISBN 0-13-111892-7
*cite book | author=Liboff, Richard L. | title=Introductory Quantum Mechanics (4th ed.) | publisher=Addison-Wesley |year=2003 |id=ISBN 0-8053-8714-5
*cite book | author=Sakurai, J. J. | title=Modern Quantum Mechanics | publisher=Addison-Wesley |year=1993 |id=ISBN 0-201-53929-2
*cite book | author=Bender, Carl; Orszag, Steven | title=Advanced Mathematical Methods for Scientists and Engineers | publisher=McGraw-Hill | year=1978 | id=ISBN 0-07-004452-X
*cite book | author=Olver, Frank J. W. | title=Asymptotics and Special Functions | publisher=Academic Press | year=1974 | id=ISBN 0-12-525850-X
*cite journal | author=Winitzki, Sergei | year=2005 | title=Cosmological particle production and the precision of the WKB approximation | journal=Physical Review D | volume=72 | pages=104011 | doi=10.1103/PhysRevD.72.104011

Historical references

*cite book | author=Carlini, Francesco | year=1817 | title=Richerche sulla convergenza della serie che serva aal soluzione del problema di Keplero | publisher=Milano
*cite journal | author=Liouville, Joseph | year=1837 | title=Sur le développement des fonctions et séries..| journal=Journal de Mathématiques Pures et Appliquées | volume=1 | pages=16–35
*cite journal | author=Green, George | year=1837 | title=On the motion of waves in a variable canal of small depth and width | journal=Transactions of the Cambridge Philosophical Society | volume=6 | pages=457–462
*cite journal | author=Rayleigh, Lord (John William Strutt) | year=1912 | title=On the propagation of waves through a stratified medium, with special reference to the question of reflection | journal=Proceedings of the Royal Society London, Series A | volume=86 | pages=207–226 | doi=10.1098/rspa.1912.0014
*cite journal | author=Gans, Richard | year=1915 | title=Fortplantzung des Lichts durch ein inhomogenes Medium | journal=Annalen der Physik | volume=47 | pages=709–736
*cite journal | author=Jeffreys, Harold | year=1924 | title=On certain approximate solutions of linear differential equations of the second order | journal=Proceedings of the London Mathematical Society | volume=23 | pages=428–436 | doi=10.1112/plms/s2-23.1.428
*cite journal | author=Brillouin, Léon | year=1926 | title=La mécanique ondulatoire de Schrödinger: une méthode générale de resolution par approximations successives | journal=Comptes Rendus de l'Academie des Sciences | volume=183 | pages=24–26
*cite journal | author=Kramers, Hendrik A. | year=1926 | title=Wellenmechanik und halbzählige Quantisierung | journal=Zeitschrift der Physik | volume=39 |pages=828–840 | doi=10.1007/BF01451751
*cite journal | author=Wentzel, Gregor | year=1926 | title=Eine Verallgemeinerung der Quantenbedingungen für die Zwecke der Wellenmechanik | journal=Zeitschrift der Physik | volume=38 | pages=518–529 | doi=10.1007/BF01397171

* Richard Fitzpatrick, " [http://farside.ph.utexas.edu/teaching/jk1/lectures/node70.html The W.K.B. Approximation] " (2002). (An application of the WKB approximation to the scattering of radio waves from the ionosphere.)
* [http://www.asymptotus.com Free WKB library for Microsoft Visual C v6 for some special functions]

Wikimedia Foundation. 2010.

### Look at other dictionaries:

• WKB (disambiguation) — WKB may mean: *WKB approximation in mathematics is an asymptotic approximation to the solution of a differential equation. *Well known binary in Geographical Information Systems is a markup language for representing vector geometry objects on a… …   Wikipedia

• WKB-Methode — Die semiklassische WKB Näherung aus der Quantenmechanik, die nach Gregor Wentzel, Hendrik Anthony Kramers und Leon Brillouin benannt ist, liefert eine Näherung der Lösung der eindimensionalen, stationären Schrödingergleichung. Die Lösung der… …   Deutsch Wikipedia

• WKB-Verfahren — Die semiklassische WKB Näherung aus der Quantenmechanik, die nach Gregor Wentzel, Hendrik Anthony Kramers und Leon Brillouin benannt ist, liefert eine Näherung der Lösung der eindimensionalen, stationären Schrödingergleichung. Die Lösung der… …   Deutsch Wikipedia

• WKB-Näherung — Die semiklassische WKB Näherung aus der Quantenmechanik, die nach Gregor Wentzel, Hendrik Anthony Kramers und Leon Brillouin benannt ist, liefert eine Näherung der Lösung der eindimensionalen, stationären Schrödingergleichung. Die Näherung… …   Deutsch Wikipedia

• Wkb — Cette page d’homonymie répertorie les différents sujets et articles partageant un même nom. {{{image}}}   Sigles d une seule lettre   Sigles de deux lettres > Sigles de trois lettres …   Wikipédia en Français

• Approximation BKW — Pour les articles homonymes, voir BKW. En physique, l approximation BKW Brillouin[1] Kramers [2] Wentzel[3] est une m …   Wikipédia en Français

• WKB — Cette page d’homonymie répertorie les différents sujets et articles partageant un même nom.   Sigles d’une seule lettre   Sigles de deux lettres > Sigles de trois lettres   Sigles de quatre lettres …   Wikipédia en Français

• WKB — Wentzell Krammers Brillouin Approximation Contributor: CASI …   NASA Acronyms

• Aproximación WKB — En Física, la aproximación WKB es el ejemplo más familiar de un cálculo semi clásico (ver antigua teoría cuántica) en el cual la función de onda es redactada como una función exponencial, semiclásicamente expandida, y entonces o la amplitud o la… …   Wikipedia Español

• Eikonal approximation — The eikonal approximation is a method of approximation useful in wave scattering equations within the realms of quantum mechanics, optics, quantum electrodynamics, and partial wave expansionInformal descriptionThe main advantage the eikonal… …   Wikipedia