HOME

TheInfoList



OR:

In
mathematical physics Mathematical physics refers to the development of mathematical methods for application to problems in physics. The '' Journal of Mathematical Physics'' defines the field as "the application of mathematics to problems in physics and the developme ...
, the WKB approximation or WKB method is a method for finding approximate solutions to linear differential equations with spatially varying coefficients. It is typically used for a semiclassical calculation in
quantum mechanics Quantum mechanics is a fundamental theory in physics that provides a description of the physical properties of nature at the scale of atoms and subatomic particles. It is the foundation of all quantum physics including quantum chemistr ...
in which the wavefunction is recast as an exponential function, semiclassically expanded, and then either the amplitude or the phase is taken to be changing slowly. The name is an initialism for Wentzel–Kramers–Brillouin. It is also known as the LG or Liouville–Green method. Other often-used letter combinations include JWKB and WKBJ, where the "J" stands for Jeffreys.


Brief history

This method is named after physicists
Gregor Wentzel Gregor Wentzel (17 February 1898 – 12 August 1978) was a German physicist known for development of quantum mechanics. Wentzel, Hendrik Kramers, and Léon Brillouin developed the Wentzel–Kramers–Brillouin approximation in 1926. In his early ...
,
Hendrik Anthony Kramers Hendrik Anthony "Hans" Kramers (17 December 1894 – 24 April 1952) was a Dutch physicist who worked with Niels Bohr to understand how electromagnetic waves interact with matter and made important contributions to quantum mechanics and statistical ...
, and
Léon Brillouin Léon Nicolas Brillouin (; August 7, 1889 – October 4, 1969) was a French physicist. He made contributions to quantum mechanics, radio wave propagation in the atmosphere, solid state physics, and information theory. Early life Brillouin ...
, who all developed it in 1926. In 1923, mathematician
Harold Jeffreys Sir Harold Jeffreys, FRS (22 April 1891 – 18 March 1989) was a British mathematician, statistician, geophysicist, and astronomer. His book, ''Theory of Probability'', which was first published in 1939, played an important role in the revival ...
had developed a general method of approximating solutions to linear, second-order differential equations, a class that includes the
Schrödinger equation The Schrödinger equation is a linear partial differential equation that governs the wave function of a quantum-mechanical system. It is a key result in quantum mechanics, and its discovery was a significant landmark in the development of th ...
. The Schrödinger equation itself was not developed until two years later, and Wentzel, Kramers, and Brillouin were apparently unaware of this earlier work, so Jeffreys is often neglected credit. Early texts in quantum mechanics contain any number of combinations of their initials, including WBK, BWK, WKBJ, JWKB and BWKJ. An authoritative discussion and critical survey has been given by Robert B. Dingle. Earlier appearances of essentially equivalent methods are:
Francesco Carlini Francesco Carlini (January 7, 1783 – August 29, 1862) was an Italian astronomer. Born in Milan, he became director of the Brera Astronomical Observatory there in 1832. He published ''Nuove tavole de moti apparenti del sole'' in 1832. In 1810 ...
in 1817,
Joseph Liouville Joseph Liouville (; ; 24 March 1809 – 8 September 1882) was a French mathematician and engineer. Life and work He was born in Saint-Omer in France on 24 March 1809. His parents were Claude-Joseph Liouville (an army officer) and Thérèse ...
in 1837, George Green in 1837,
Lord Rayleigh John William Strutt, 3rd Baron Rayleigh, (; 12 November 1842 – 30 June 1919) was an English mathematician and physicist who made extensive contributions to science. He spent all of his academic career at the University of Cambridge. Am ...
in 1912 and
Richard Gans Richard is a male given name. It originates, via Old French, from Old Frankish and is a compound of the words descending from Proto-Germanic ''*rīk-'' 'ruler, leader, king' and ''*hardu-'' 'strong, brave, hardy', and it therefore means 'stron ...
in 1915. Liouville and Green may be said to have founded the method in 1837, and it is also commonly referred to as the Liouville–Green or LG method. The important contribution of Jeffreys, Wentzel, Kramers, and Brillouin 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.


Formulation

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 : \varepsilon \frac + a(x)\frac + \cdots + k(x)\frac + m(x)y= 0, assume a solution of the form of an
asymptotic series In mathematics, an asymptotic expansion, asymptotic series or Poincaré expansion (after Henri Poincaré) is a formal series of functions which has the property that truncating the series after a finite number of terms provides an approximation to ...
expansion : y(x) \sim \exp\left frac\sum_^\delta^nS_n(x)\right/math> in the limit . The asymptotic scaling of in terms of will be determined by the equation – see the example below. Substituting the above
ansatz In physics and mathematics, an ansatz (; , meaning: "initial placement of a tool at a work piece", plural Ansätze ; ) is an educated guess or an additional assumption made to help solve a problem, and which may later be verified to be part of the ...
into the differential equation and cancelling out the exponential terms allows one to solve for an arbitrary number of terms in the expansion. WKB theory is a special case of
multiple scale analysis In mathematics and physics, multiple-scale analysis (also called the method of multiple scales) comprises techniques used to construct uniformly valid approximations to the solutions of perturbation problems, both for small as well as large value ...
.


An example

This example comes from the text of
Carl M. Bender Carl M. Bender (born 1943) is an American Applied mathematics, applied mathematician and Mathematical physics, mathematical physicist. He currently holds the Wilfred R. and Ann Lee Konneker Distinguished Professorship of Physics at Washington Uni ...
and Steven Orszag. Consider the second-order homogeneous linear differential equation : \epsilon^2 \frac = Q(x) y, where Q(x) \neq 0. Substituting :y(x) = \exp\left frac\sum_^\infty \delta^nS_n(x)\right/math> results in the equation :\epsilon^2\left frac\left(\sum_^\infty \delta^nS_n'\right)^2 + \frac\sum_^\delta^nS_n''\right= Q(x). To leading order (assuming, for the moment, the series will be asymptotically consistent), the above can be approximated as :\fracS_0'^2 + \fracS_0'S_1' + \fracS_0'' = Q(x). In the limit , the dominant balance is given by :\fracS_0'^2 \sim Q(x). So is proportional to ''ε''. Setting them equal and comparing powers yields :\epsilon^0: \quad S_0'^2 = Q(x), which can be recognized as the
Eikonal equation An eikonal equation (from Greek εἰκών, image) is a non-linear first-order partial differential equation that is encountered in problems of wave propagation. The classical eikonal equation in geometric optics is a differential equation of ...
, with solution :S_0(x) = \pm \int_^x \sqrt\,dx'. Considering first-order powers of fixes :\epsilon^1: \quad 2S_0'S_1' + S_0'' = 0. This has the solution :S_1(x) = -\frac\ln Q(x) + k_1, where is an arbitrary constant. We now have a pair of approximations to the system (a pair, because can take two signs); the first-order WKB-approximation will be a linear combination of the two: :y(x) \approx c_1Q^(x)\exp\left frac\int_^x\sqrt\,dt\right+ c_2Q^(x)\exp\left \frac\int_^x\sqrt\,dt\right Higher-order terms can be obtained by looking at equations for higher powers of . Explicitly, : 2S_0'S_n' + S''_ + \sum_^S'_jS'_ = 0 for ≥ 2.


Precision of the asymptotic series

The asymptotic series for is usually a
divergent series In mathematics, a divergent series is an infinite series that is not convergent, meaning that the infinite sequence of the partial sums of the series does not have a finite limit. If a series converges, the individual terms of the series mus ...
, whose general term starts to increase after a certain value . 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 = Q(x) y, with <0 an analytic function, the value n_\max and the magnitude of the last term can be estimated as follows: :n_\max \approx 2\epsilon^ \left, \int_^ \sqrt\,dz \ , :\delta^S_(x_0) \approx \sqrt \exp n_\max where x_0 is the point at which y(x_0) needs to be evaluated and x_ is the (complex) turning point where Q(x_)=0, closest to x=x_0. The number can be interpreted as the number of oscillations between x_0 and the closest turning point. If \epsilon^Q(x) is a slowly changing function, :\epsilon\left, \frac \ \ll Q^2 ,^ the number will be large, and the minimum error of the asymptotic series will be exponentially small.


Application to the Schrödinger equation

The above example may be applied specifically to the one-dimensional, time-independent
Schrödinger equation The Schrödinger equation is a linear partial differential equation that governs the wave function of a quantum-mechanical system. It is a key result in quantum mechanics, and its discovery was a significant landmark in the development of th ...
, :-\frac \frac \Psi(x) + V(x) \Psi(x) = E \Psi(x), which can be rewritten as :\frac \Psi(x) = \frac \left( V(x) - E \right) \Psi(x).


Approximation away from the turning points

The wavefunction can be rewritten as the exponential of another function (closely related to the
action Action may refer to: * Action (narrative), a literary mode * Action fiction, a type of genre fiction * Action game, a genre of video game Film * Action film, a genre of film * ''Action'' (1921 film), a film by John Ford * ''Action'' (1980 fil ...
), which could be complex, :\Psi(x) = e^, so that :\Phi''(x) + \left Phi'(x)\right2 = \frac \left( V(x) - E \right), where ' indicates the derivative of with respect to ''x''. This derivative ' can be separated into real and imaginary parts by introducing the real functions ''A'' and ''B'', :\Phi'(x) = A(x) + i B(x). The amplitude of the wavefunction is then :\exp\left int_^x A(x')\,dx'\right while the phase is :\int_^x B(x')\,dx'. The real and imaginary parts of the Schrödinger equation then become :A'(x) + A(x)^2 - B(x)^2 = \frac \left( V(x) - E \right), :B'(x) + 2 A(x) B(x) = 0. Next, the semiclassical approximation is used. This means that each function is expanded as a power series in . From the above equations, it can be seen that the power series must start with at least an order of 1/ to satisfy the real part of the equation: :A(x) = \frac \sum_^\infty \hbar^n A_n(x), :B(x) = \frac \sum_^\infty \hbar^n B_n(x). To the zeroth order in this expansion, the conditions on ''A'' and ''B'' can be written, :A_0(x)^2 - B_0(x)^2 = 2m \left( V(x) - E \right), :A_0(x) B_0(x) = 0 \;. The first derivatives and were discarded, because they include factors of order 1/, higher than the dominant −2. Then, if the amplitude varies sufficiently slowly as compared to the phase (A_0(x) = 0), it follows that :B_0(x) = \pm \sqrt, 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(x) \approx C_0 \frac.
On the other hand, if it is the phase that varies slowly (as compared to the amplitude), (B_0(x) = 0) then :A_0(x) = \pm \sqrt, which is only valid when the potential energy is greater than the total energy (the regime in which
quantum tunneling In physics, a quantum (plural quanta) is the minimum amount of any physical entity (physical property) involved in an interaction. The fundamental notion that a physical property can be "quantized" is referred to as "the hypothesis of quantizati ...
occurs). Finding the next order of the expansion yields, as in the example of the previous section, In the classically allowed region, namely the region where V(x) the integrand in the exponent is imaginary and the approximate wave function is oscillatory. In the classically forbidden region V(x)>E, the solutions are growing or decaying. It is evident in the denominator that both of these approximate solutions become singular near the classical turning points, where , and cannot be valid. (The turning points are the points where the classical particle changes direction.)


Behavior near the turning points

We now consider the behavior of the wave function near the turning points. For this, we need a different method. Near the first turning points, 1, the term \frac\left(V(x)-E\right) can be expanded in a power series, :\frac\left(V(x)-E\right) = U_1 \cdot (x - x_1) + U_2 \cdot (x - x_1)^2 + \cdots\;. To first order, one finds :\frac \Psi(x) = U_1 \cdot (x - x_1) \cdot \Psi(x). This differential equation is known as the
Airy equation In the physical sciences, the Airy function (or Airy function of the first kind) is a special function named after the British astronomer George Biddell Airy (1801–1892). The function and the related function , are linearly independent solutio ...
, and the solution may be written in terms of
Airy function In the physical sciences, the Airy function (or Airy function of the first kind) is a special function named after the British astronomer George Biddell Airy (1801–1892). The function and the related function , are linearly independent solut ...
s, :\Psi(x) = C_A \textrm\left( \sqrt \cdot (x - x_1) \right) + C_B \textrm\left( \sqrt \cdot (x - x_1) \right). Although for any fixed value of \hbar, the wave function is bounded near the turning points, the wave function will be peaked there, as can be seen in the images above. As \hbar gets smaller, the height of the wave function at the turning points grows.


The matching conditions

It now remains to construct a global (approximate) solution to the Schrödinger equation. For the wave function to be square-integrable, we must take only the exponentially decaying solution in the two classically forbidden regions. These must then "connect" properly through the turning points to the classically allowed region. For most values of ''E'', this matching procedure will not work: The function obtained by connecting the solution near +\infty to the classically allowed region will not agree with the function obtained by connecting the solution near -\infty to the classically allowed region. The requirement that the two functions agree imposes a condition on the energy ''E'', which will give an approximation to the exact quantum energy levels. Given the two 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 the Airy function to connect them. Thus, a relationship between C_0,\theta and C_,C_ can be found. This relationship is obtained using known asymptotic of the Airy function. The relationship can be found to be as follows (often referred to as the "connection formulas"): : C_ = + \frac C_0 \cos, : C_ = - \frac C_0 \sin. Now the global (approximate) solutions can be constructed. The same can be done at the other turning points; assume there is just another one, 2. The expression there, however, will appear different than the one determined above at 1 by a difference in the argument of these trigonometric functions. The matching condition, needed to get a single-valued, square-integrable approximate solution, takes the following form: ::\int_^ \sqrt\,dx = (n+1/2)\pi \hbar , where x_1,x_2 are the turning points of the potential discussed, where the integrand vanishes. Here ''n'' is a non-negative integer. This condition can also be rewritten as saying that ::The area enclosed by the classical energy curve is 2\pi\hbar(n+1/2). Either way, the condition on the energy is a version of the Bohr–Sommerfeld quantization condition, with a " Maslov correction" equal to 1/2. It is possible to show that after piecing together the approximations in the various regions, one obtains a good approximation to the actual eigenfunction. In particular, the Maslov-corrected Bohr–Sommerfeld energies are good approximations to the actual eigenvalues of the Schrödinger operator. Specifically, the error in the energies is small compared to the typical spacing of the quantum energy levels. Thus, although the "old quantum theory" of Bohr and Sommerfeld was ultimately replaced by the Schrödinger equation, some vestige of that theory remains, as an approximation to the eigenvalues of the appropriate Schrödinger operator.


The probability density

One can then compute the probability density associated to the approximate wave function. The probability that the quantum particle will be found in the classically forbidden region is small. In the classically allowed region, meanwhile, the probability the quantum particle will be found in a given interval is approximately the ''fraction of time the classical particle spends in that interval'' over one period of motion. Since the classical particle's velocity goes to zero at the turning points, it spends more time near the turning points than in other classically allowed regions. This observation accounts for the peak in the wave function (and its probability density) near the turning points. Applications of the WKB method to Schrödinger equations with a large variety of potentials and comparison with perturbation methods and path integrals are treated in Müller-Kirsten.Harald J.W. Müller-Kirsten, Introduction to Quantum Mechanics: Schrödinger Equation and Path Integral, 2nd ed. (World Scientific, 2012).


See also

*
Instanton An instanton (or pseudoparticle) is a notion appearing in theoretical and mathematical physics. An instanton is a classical solution to equations of motion with a finite, non-zero action, either in quantum mechanics or in quantum field theory. Mo ...
*
Airy function In the physical sciences, the Airy function (or Airy function of the first kind) is a special function named after the British astronomer George Biddell Airy (1801–1892). The function and the related function , are linearly independent solut ...
* Field electron emission * Langer correction *
Maslov index In mathematics, the Lagrangian Grassmannian is the smooth manifold of Lagrangian subspaces of a real symplectic vector space ''V''. Its dimension is ''n''(''n'' + 1) (where the dimension of ''V'' is ''2n''). It may be identified with the homogeneous ...
*
Method of steepest descent In mathematics, the method of steepest descent or saddle-point method is an extension of Laplace's method for approximating an integral, where one deforms a contour integral in the complex plane to pass near a stationary point (saddle point), in r ...
*
Method of matched asymptotic expansions In mathematics, the method of matched asymptotic expansions is a common approach to finding an accurate approximation to the solution to an equation, or system of equations. It is particularly used when solving singularly perturbed differential equ ...
*
Old quantum theory The old quantum theory is a collection of results from the years 1900–1925 which predate modern quantum mechanics. The theory was never complete or self-consistent, but was rather a set of heuristic corrections to classical mechanics. The theory ...
*
Einstein–Brillouin–Keller method The Einstein–Brillouin–Keller method (EBK) is a semiclassical method (named after Albert Einstein, Léon Brillouin, and Joseph B. Keller) used to compute eigenvalues in quantum-mechanical systems. EBK quantization is an improvement from Boh ...
*
Perturbation methods In mathematics and applied mathematics, perturbation theory comprises methods for finding an approximate solution to a problem, by starting from the exact solution of a related, simpler problem. A critical feature of the technique is a middle ...
*
Quantum tunneling In physics, a quantum (plural quanta) is the minimum amount of any physical entity (physical property) involved in an interaction. The fundamental notion that a physical property can be "quantized" is referred to as "the hypothesis of quantizati ...
* Slowly varying envelope approximation *
Supersymmetric WKB approximation In physics, the supersymmetric WKB (SWKB) approximation is an extension of the WKB approximation that uses principles from supersymmetric quantum mechanics to provide estimations on energy eigenvalues in quantum-mechanical systems. Using the supe ...


References


Modern references

* * * * * * * *


Historical references

* * * * * * * * *


External links

* {{cite web, first=Richard , last=Fitzpatrick, url=http://farside.ph.utexas.edu/teaching/jk1/lectures/node70.html, title= The W.K.B. Approximation, year=2002 (An application of the WKB approximation to the scattering of radio waves from the ionosphere.) Approximations Theoretical physics Asymptotic analysis Mathematical physics