In quantum mechanics, the Wentzel-Kramers-Brillouin or simply the WKB
approximation is a technique to approximate the wave function ψ(x) of
the one-dimensional time-independent Schrödinger equation. It is an example
of a semiclassical approximation, because it tries to find a
balance between classical and quantum physics.
In classical mechanics, a particle travelling in a potential V(x)
along a path x(t) has a total energy E as follows, which we
rearrange:
E=21mx˙2+V(x)⟹m2(x′)2=2m(E−V(x))
The left-hand side of the rearranged version is simply the momentum squared,
so we define the magnitude of the momentum p(x) accordingly:
p(x)=2m(E−V(x))
Note that this is under the assumption that E>V,
which is always true in classical mechanics,
but not necessarily in quantum mechanics.
We rewrite the Schrödinger equation:
0=dx2d2ψ+ℏ22m(E−V)ψ=dx2d2ψ+ℏ2p2ψ
If V(x) were constant, and by extension p(x) too, then the solution
is easy:
ψ(x)=ψ(0)exp(±ipx/ℏ)
This form is reminiscent of the generator of translations. In practice,
V(x) and p(x) vary with x, but we can still salvage this solution
by assuming that V(x) varies slowly compared to the wavelength
λ(x)=2π/k(x), where k(x)=p(x)/ℏ is the
wavenumber. The solution then takes the following form:
ψ(x)=ψ(0)exp(±ℏi∫0xχ(ξ)dξ)
χ(ξ) is an unknown function, which intuitively should be related
to p(x). The purpose of the integral is to accumulate the change of
χ from the initial point 0 to the current position x.
Let us write this as an indefinite integral for convenience:
ψ(x)=ψ(0)exp(±ℏi(∫χ(x)dx−C))
Where C=∫χ(x)dx∣x=0 is the initial point of the definite integral.
For simplicity, we absorb the constant C into ψ(0).
We can now clearly see that:
ψ′(x)=±ℏiχ(x)ψ(x)⟹χ(x)=±iℏψ(x)ψ′(x)
Next, we insert this ansatz for ψ(x) into the Schrödinger equation
to get:
Dividing out ψ and rearranging gives us the following, which is
still exact:
±iℏχ′=p2−χ2
Next, we expand this as a power series of ℏ. This is why it is
called semiclassical: so far we have been using full quantum mechanics,
but now we are treating ℏ as a parameter which controls the
strength of quantum effects:
χ(x)=χ0(x)+iℏχ1(x)+i2ℏ2χ2(x)+...
The heart of the WKB approximation is its assumption that quantum effects are
sufficiently weak (i.e. ℏ is small enough) that we only need to
consider the first two terms, or, more specifically, that we only go up to
ℏ, not ℏ2 or higher. Inserting the first two terms of this
expansion into the equation:
±iℏχ0′=p2−χ02−2iℏχ0χ1
Where we have discarded all terms containing ℏ2. At order
ℏ0, we then get the expected classical result for χ0(x):
0=p2−χ02⟹χ0(x)=p(x)
While at order ℏ, we get the following quantum-mechanical
correction:
±iℏχ0′=−2iℏχ0χ1⟹χ1(x)=∓21χ0(x)χ0′(x)
Therefore, our approximated wave function ψ(x) currently looks like
this:
ψ(x)≈ψ(0)exp(±ℏi∫χ0(x)dx)exp(±∫χ1(x)dx)
We can reduce the latter exponential using integration by substitution:
In the WKB approximation for E>V, the solution ψ(x) is thus
given by:
ψ(x)≈p(x)Aexp(±ℏi∫p(x)dx)
What if E<V? In classical mechanics, this is just not allowed; a ball
cannot simply go through a potential bump without the necessary energy.
On the other hand, in quantum physics, particles can tunnel through barriers.
Luckily, the only thing we need to change for the WKB approximation
is to let the momentum take imaginary values:
p(x)=2m(E−V(x))=i2m(V(x)−E)
And then take the absolute value in the appropriate place in front of ψ(x):
ψ(x)≈∣p(x)∣Aexp(±ℏi∫p(x)dx)
In the classical region (E>V), the wave function oscillates, and
in the quantum-physical region (E<V) it is exponential.
Note that for E≈V the approximation breaks down,
because of the appearance of p(x) in the denominator.