Michael Fowler 1/21/08
The WKB (Wentzel, Kramers, Brillouin) approximation is,
in sense to be made clear below, a
quasi-classical method for solving the one-dimensional (and effectively
one-dimensional, such as radial) time-independent Schrödinger equation. The nontrivial step in the method is the
connection formulas (see below), that problem was first solved by Lord Rayleigh
(Proc. Roy. Soc. A, 86, 1912, 207) and as Jeffries notes (Mathematical Physics, p 526) “it has been rediscovered by several
later writers” presumable referring to W, K and B. By the way, the English call it the Jeffries approximation,
or, if feeling sufficiently ecumenical, the WKBJ approximation. (In this lecture, we only consider bound
states: the most famous application of WKB, α-decay,
was covered in detail in the undergraduate quantum mechanics course based on
We’ll follow the development in Landau and Lifshitz, who consider this all sufficiently obvious that they don’t mention any of these people. In fact, they call it
Consider a particle moving along in a slowly varying one-dimensional potential. By “slowly varying” we mean here that in any small region the wave function is well approximated by a plane wave, and that the wavelength only changes over distances long compared with a wavelength. We’re also assuming for the moment that the particle has positive kinetic energy in the region. Under these conditions, it’s easy to see the general form of the solution to the time independent Schrödinger equation
Very approximately, will look like where p(x) is the “local momentum” we define classically by
and A(x) is slowly varying compared with the phase factor.
Clearly this is a semiclassical limit: has to be sufficiently small that there are many oscillations in the typical distance over which the potential varies.
To handle this a little more precisely, we emphasize the rapid phase variation in this semiclassical limit by writing the wave function
and writing Schrödinger ’s equation for .
Schrödinger’s equation written for the phase function is:
And, since we’re assuming the system is close to classical, it makes sense to expand as a series in (following Landau and Lifshitz):
The zeroth order approximation is
and fixing the sign of p by
we conclude that
(As we discussed in the lecture on path integrals, in the classical limit one path dominates, and the phase of the wave function is times the classical action S along that path. In the present case, , we’ve already factored out the Et since we’re dealing here with the time-independent wave function.)
From the Schrödinger equation , it is evident that this approximate solution is only valid if we can ignore that first term. That is to say, we must have
But in leading approximation and , so the condition is
This just means the change in wavelength over a distance of one wavelength must be small.
Obviously, this cannot always be the case: if the particle is confined by an attractive potential, at the edge of the classically allowed region, that is, where , p is zero and the wavelength is infinite. The approximation is only good well away from that point, to which we shall return shortly.
The second term in theexpansion of the phase, satisfies
So the wave function to this order is:
(Recall we fixed the sign of p to be positive.)
To interpret the factor, consider the first term, a wave moving to the right. Since p is real, the exponential has modulus unity, and the local amplitude squared is proportional to 1/p, that is, 1/v , where v is the velocity of the particle. This is simple to understand physically: the probability of finding the particle in any given small interval is proportional to the time it spends there, hence inversely proportional to its speed.
We turn now to the wave function in the classically forbidden region,
Here p is of course pure imaginary, but the same formal phase solution of the Schrodinger equation works, again provided that the particle is well away from the points where E = V(x).
The wave function is:
Let us assume we’re dealing with a one-dimensional potential, and the classically allowed region is (I’m just following Landau’s notation here.) Clearly, in the forbidden region to the right, only the first term in the above equation for appears, and for only the second term. Furthermore, in the “inside” (classically allowed) region, the wave function has the oscillating form discussed earlier.
But how do we connect the three regions together? We make an assumption: we take it that the potential varies sufficiently smoothly that it’s a good approximation to take it to be linear in the vicinity of the classical turning points. That is to say, we assume that a linear potential is a sufficiently good approximation out to the point where the short wavelength (or decay length for tunneling regions) description is adequate.
Therefore, near x = a, we take the potential to be
(so F0 would be the force) and then approximate the wave function by the known exact solution for a linear everywhere potential: the Airy function.
It is known that for the Airy function, the solution having the form
to the right becomes
for (The derivation of this “connection” is given in my notes here.)
At b, the same argument gives
For these two expressions to be consistent, we must have
where the latter integral is over a complete cycle of the classical motion.
Here n is the number of zeros of the wave function: this is the quantization condition.
The time for a complete classical circuit is
Note that is the area of the classical path in phase space, so we see each state has an element of phase space From this, we can figure out the approximate energy splitting between levels in the quasi-classical limit: the change in the integral with energy corresponding to one level must be That is,
This is just saying that if the particle emits one photon and drops to the next level, the frequency of the photon emitted is just the orbital frequency of the particle, a very natural conclusion in the quasi-classical limit.
In the above analysis for a particle confined to one dimension, the connection formulas can be understood with a simple picture: the wave function “spills over” into the forbidden regime, and its twisting there counts as an extra of phase change, so in the lowest state the total phase change in the allowed region need only be . In the radial case, assuming the potential is well behaved at the origin, the wave function goes to zero there. A bound state will still spill over beyond the classical turning point at r0, say, but clearly there must be a total phase change of in the allowed region for the lowest state, since there can be no spill over to negative r.
The general formula will be
the series terminating if and when the turning point reaches infinity.
Warning: actually, some potentials, including the Coulomb potential and the centrifugal barrier for , are in fact singular at These cases require special treatment.