Jump to content

User:Gareth Owen/WKB approximation

From Wikipedia, the free encyclopedia

In physics, the WKB (Wentzel-Kramers-Brillouin) approximation, also known as WKBJ 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.

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.

Derivation

[edit]

We beginning with a one dimensional, time-independent wave equation in which the local wavenumber varies. Such an equation can typically be written as

,

where K is O(1) and is small.

We recast the wavefunction as the exponential of another function Φ (which is closely related to the action):

The function Φ must then satisfy

where Φ' indicates the derivative of Φ with respect to x. Now let us separate into real and imaginary parts by introducing the real functions A and B:

The amplitude of the wavefunction is then while its phase is . The governing equation implies that these functions must satisfy:

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

Next we want to find an asymptotic approximation to solve this. That means we expand each function as a power series in . From the equations we can already see that the power series must start with at least an order of to satisfy the real part of the equation.

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

Clearly then, the second of these equations tells us that either or must be identically zero. Since both functions are real, examination of the first equation tells us that if then and if , then

Oscillatory Regime

[edit]

In the former case, the leading order term

constitutes a rapid variation in phase. Thus, the solutions to the equation are predominately oscillatory in nature. For this case, we can calculate the next order correction. Taking the real and imaginary parts of term of order , and noting that we have

and

giving us

Thus, to the first two orders

Integrating directly, we have

Exponential Regime

[edit]

When (corresponding to an imaginary wavenumber, K) the leading order term is given by

leading to either exponential growth or decay. At the next order, noting that , we have

and

Therefore, and, as in the oscillatory regime,

Combining these two terms we obtain

It is clear that in neither regime is the approximation valid near , where the denominator becomes singular, as the assumption that is breaks down. In wave dynamics, particularly optics, the location where this happens is known as a caustic.

Across the Caustic

[edit]

It is apparent from the denominator, that both of these approximate solutions 'blow up' when the local wavenumber passes through zero, and cannot be valid. The approximate solutions that we have found are accurate away from this zero, but inaccurate near to it. We can find accurate approximate solutions near to this zero by approximating by its Taylor series

Let's label the zero of by . Now, if is near , we can write

,

where is the derivative of at

To first order, one finds

This differential equation is Airy equation, and the solution may be written in terms of Airy functions. (Alternatively, with some trickery, it may be transformed into a Bessel equation of fractional order.) The exact form of the solution depends on the sign of . The case when corresponds to the oscillatory regime being to the left of and the exponential regime to the right. In this case, the solution is given by,

When , the locations of the oscillatory and exponential regimes are reversed, and the solution is:

This solution is accurate near the zero, and should connect with the solutions . Thus, we should be able to determine the 2 coefficients and so that the solutions are identical in the region of overlap, where the are both accurate.

See also

[edit]

Perturbation methods, Quantum tunnelling, Airy Function

References

[edit]
  • Razavy, Moshen (2003). Quantum Theory of Tunneling. World Scientific. ISBN 9812380191.
  • Griffiths, David J. (2004). Introduction to Quantum Mechanics (2nd ed.). Prentice Hall. ISBN 0131118927.
  • Liboff, Richard L. (2003). Introductory Quantum Mechanics (4th ed.). Addison-Wesley. ISBN 0805387145.
  • Sakurai, J. J. (1993). Modern Quantum Mechanics. Addison-Wesley. ISBN 0201539292.
[edit]
  • The W.K.B. Approximation (Note that in this webpage, : there are two sets of labels for the equations.)