221A Lecture Notes

WKB Method
1

HamiltonJacobi Equation

We start from the Schrodinger equation for a single particle in a potential

h
2 ~ 2
i
h (~x, t) =
+ V (~x) (~x, t).
t
2m
"

(1)

We can rewrite this equation by using (~x, t) = eiS(~x,t)/h :

"

i
h ~2
1 ~ 2
S
=
(S)
( S) + V .

t
2m
2m

(2)

Assuming 6= 0, this leads to an equation

S
1 ~ 2
i
h ~2
=
(S)
( S) + V.
t
2m
2m

(3)

Now taking the formal limit h

0, it becomes the same as the classical
HamilgonJacobi equation

1 ~ 2
S
=
(S) + V.
t
2m

(4)

It is often said that the h

0 limit is indeed the classical mechanics because
it reduces to the HamiltonJacobi equation.

Classical Limit

The fact that the exponent S satisfies the HamiltonJacobi equation may
not mean the classical behavior to you right away. The way to see it is by
forming a wave packet, and ask the question how does the wave packet moves
as a whole. To simplify the discussion, let us consentrate on the case where
the Hamiltonian does not depend explicitly on time, and we can separate the
variable as
x, E) Et.
S(~x, t) = S(~
(5)
1

The inverse Legendre transformation says

x, E)
S(~
t=
.
(6)
E
This condition can be viewed as an implicit equation to determine the position of the particle ~x as a function of time t.
To form a wave packet, we need to put together waves of slighly different
energies. We can write
(~x, t) =

dEg(E)ei(S(~x,E)Et)/h ,

(7)

The question here is where the wave packet is. At most places, the averaging
over E makes the phase oscillate very fast because of the small h
, and there
is a strong cancellation. Only at special points where the phase factor is
stationary with respect to E, the wave function is sizable. Therefore, the
position of the wave packet is determined by the stationary condition

x, E)
S(~

t = 0.
(8)
S(~x, E) Et =
E
E
This is nothing but Eq. (6). In other words, the wave packet is sizable only
where the position satisfies the classical equation of motion. Therefore, the
wave packet follows the classical equation of motion.

h
Expansion

Now that we have rederived classical mechanics from quantum mechanics,

our interest now turns to quantum mechanics itself. We will use the fact
that h
is indeed small to find a useful approximation method in quantum
mechanics.
The Eq. (3) is exact as long as 6= 0. One can use this equation, and
consider an expansion
S(~x, t) = S0 + h
S1 + h
2 S2 + .

(9)

This is an expansion in h
, and hence called h
-expansion or semi-classical
expansion. Plugging in the expansion into Eq. (3), we find
S0
1 ~
=
(S0 )2 + V,
t
2m
i
S1
1 h ~2
~ 0 )(S
~ 1) ,
i S0 + 2(S

=
t
2m

(10)
(11)

and similarly for higher terms in h

. The leading equation has only S0 , and
it is exactly the same as HamiltonJacobi equation. Once you solve these
equations starting from S0 , S1 , etc, you have solved the wave function in
a systematic expansion in h
.
There is one important difference between the classical limit which
reduces to the HamiltonJacobi equation and the idea of expansion in h
.
In the classical limit, S is of course real because it does not make sense
otherwise. However, in the sense of an expansion in h
, we dont know if S0
has to be real. In fact, in many interesting cases, S0 turns out to be complex .
We will see these examples below.

WKB Approximation

WKB Approximation, due to Wentzel, Kramers, and Brillouin, keeps terms

up to O(
h) in the h
expansion. It is used mostly for the time-independent
case, or in other words, for an eigenstate of energy E. Then the wave function
has the ordinary time dependence eiEt/h . We also restrict ourselves to the
one-dimensional problem. In terms of S, it corresponds to
S(~x, t) = S(~x) Et.

(12)

Therefore only S0 has the time dependence S0 (x, t) = S0 (x) Et, while
higher order terms Si = Si (x) for i 6= 0 do not depend on time.
The lowest order term S0 satisfies the HamiltonJacobi equation (see
Eq. (10)
1
(S 0 )2 + V (x).
(13)
E=
2m 0
The differential equation can be solve immediately as
S0 (x) =

2m(E V (x0 ))dx0 =

p(x0 )dx0

(14)

up to an integration constant which can be determined only after imposing

aq
boundary condition on the wave function. We used the notation p(x) =
2m(E V (x)) because it is nothing but the momentum of the particle
in the classical sense. Once we know S0 , we can also solve for S1 . Starting
from Eq. (11), and using S1 /t = 0, we find
2S00 S10 = iS000 ,
3

(15)

which has a solution

S1 (x) = i

i
S000 (x0 ) 0
dx = log p(x) + constant.
0
0
2S0 (x )
2

(16)

Therefore the general solution to the Schrodinger equation up to this order

is
(x, t) = ei(S0 (x)+hS1 (x))/h eiEt/h


1
i Z xq
0
0
= c
2m(E V (x ))dx eiEt/h ,
exp
p(x)1/2
h

(17)

and the overall constant c is of course undetermined from this analysis.

This solution makes it immediately clear that this approximation breaks
down when p(x) goes to zero. Or in other words, the approximation is bad
where the classical particle stops and turns because of the potential. Such
points are called classical turning points.

4.1

Validity of the WKB Approximation

The approximation to stop with S1 in the h

expansion is valid only when S1
is much smaller than S0 . Or in other words, if the term with h
in Eq. (2) is
much smaller than the other terms. In particular, we require
~ 2|  h
~ 2 S|.
|(S)
|

(18)

In the one-diimensional time-independent case discussed above, this is

p(x)2  h
|p0 (x)|.

(19)

Using the definition of p(x) = 2m(E V (x)), we find




h
dV (x)/dx 


  1.
 2(E V (x))p(x) 

(20)

Again we find the same conclusion: the WKB approximation breaks down
close to the classical turning point V (x) = E (e.g., p(x) = 0).
For example, take a harmonic oscillator V (x) = 21 m 2 x2 . The validity
condition Eq. (20) can be rewritten as

3


1
2 2

8 E m x   (
h)2 m 2 x2 .

(21)

This inequality is always satisfied exactly at the origin x = 0, but once away
from the origin, it is impossible to satisfy unless E  h
. In this sense, we
are indeed in the classical regime. However, even for a large E  h
, the
1
approximation is not valid close to the classical turning points E = 2 m 2 x2 .
Here is the surprise. The validity condition Eq. (20) may be satisfied
even in the region where the particle cannot enter classicaly E < V (x). For
example with the harmonic
q oscillator again, the validity condition is always
satisfied for large x  2E/m 2 for any value of E. In other words, the
WKB approximation is good away from the classical turning points both
where a classical particle exists and where a classical particle cannot exist.
This is why the WKB approximation is not really a classical limit. It applies
also where physics is truly quantum mechanical.
In the classically forbidden region, the solution Eq. (17) needs to be modified to
(x, t) = ei(S0 (x)+hS1 (x))/h eiEt/h


1 Z xq
1
0
0
2m(V (x ) E)dx eiEt/h ,(22)
exp
=c
(2m(V (x) E))1/4
h

by following the same steps as in the classically allowed region.

4.2

Matching

WKB approximation can be good both in the region E > V (x) and the region
E < V (x) but cannot be good in between the regions close to the classical
turning point E = V (xc ). In order to utilize the WKB approximation to
work out wave functions, we need to somehow overcome this limitation. The
standard method is to expand around xc and solve for the wave function
exactly. Then you can match on to the WKB solutions away from xc to
determine the entire wave function.
The common method is to approximate the potential around the classical
turning point xc by a linear one:
V (x) = V (xc ) + V 0 (xc )(x xc ) + O(x xc )2 ,

(23)

and ignore the second order term. By definition V (xc ) = E. Schrodinger

equation around this point is therefore
h
2 d2
h
2 d2
+
V
(x)

E
+ V 0 (xc )(x xc ) = 0.

2m dx2
2m dx2
!

(24)

Using the new variable

!1/3

u=

2m dV
(xc )
h
2 dx

(x xc ),

(25)

the differential equation simplifes drastically to

d2
u = 0.
du2
!

(26)

The solution to this equation is known as the Airy function



1Z
1 3
dt cos t + ut .
Ai(u) =
0
3

(27)

This can be checked as follows. By acting the differential operator in Eq. (26)
on the definition of the Airy function, we find
!


1Z d
1 3
d2

u
Ai(u)
=

dt
sin
t
+
ut
.
du2
0
dt
3

(28)

The boundary term at t = 0 obviously vanishes. The behavior at t = is

trickier. The point is that the argument of the sin grows as t3 and oscillates
more and more rapidly as t . Therefore for any infinitesimal interval of
large t, the oscillation basically cancels the integrand except for a left-over
that goes down as 1/t2 . Therefore the boundary term for t can also
be dropped.
It can be shown using the steepest descent method that the asymptotic
behavior of the Airy function smoothly matches to the WKB solutions. The
asymptotic behavior is

1/2

exp 23 u3/2
u0
.
Ai(u)  u1/2



1
cos
u
u
+
u

0
3
4
u
1
2

(29)

Note first that, for u  0, the asymptotic behavior is

1/2

Ai(u) =

h
(2mV 0 /
h2 )1/3

2m(E V )

1 Z xq

cos
2m(E V )dx0 +
.
h
xc
4


(30)

Here weqused the linear expansion V (x) = E +V 0 (xc )(xxc ) to relate powers
of u to 2m(E V ). This expression is consistent with the WKB solution
for a particular choice of the overall constant. Similarly, for u  0,
1/2

(2mV 0 /
h2 )1/3
1 h
q
Ai(u) =
2 2m(V E)

4.3

1 Z xq
exp
2m(V (x0 ) E)dx0 . (31)
h
xc


Bound States

Consider the following situation. In the region I (x < a), E < V and it
is classically forbidden. In the region II (a < x < b), E > V and the
the particle is classically allowed. But again in the region III (x > b) it
is classically forbidden. In this case, we expect bound states with discrete
energy levels.
From the matching using the Airy function at x = b, we want the wave
function in the classicaly allowed region II to be as in Eq. (30),
1/2

(x) = cb

h
(2mV 0 (b)/
h2 )1/3

2m(E V )

1 Z xq
2m(E V (x0 ))dx0 +
cos
.
h
b
4


(32)
On the other hand from the matching at x = a, the wave function in the
region II must be
1/2

(x) = ca

h
(2mV 0 (a)/
h2 )1/3

2m(E V )

cos

1 Z xq
2m(E V (x0 ))dx0
.
h
a
4


(33)
cb and ca are normalization constants. The minus sign in front of /4 is
there because the classically allowed region is to the right of the turning
point and hence we need to use Ai(u) instead of Ai(u) for matching. The
two behaviors obtained from matching at both sides must be the same in
order to have a consistent wave function for the entire region II. Because we
do not know the coefficients ca,b , which can in particular differ in signs, we
have to require
1 Z xq

1 Z xq

2m(E V (x0 ))dx0 + =

2m(E V (x0 ))dx0 n. (34)
h
b
4
h
a
4
7

Further simplifying the equation we obtain

Z

2m(E V (x))dx = n +

h.
2


(35)

This is reminiscent of BohrSommerfeld quantization condition that appeared at the early stage of quantum mechanics,
I

pdq = nh

(36)

with h = 2
h except that there is an extra contribution 1/2 here. Remember
that the BohrSommerfeld quantization condition was an ad hoc requirement
which happened to give correct energy spectra for the hydrogen-like atoms.
Nonetheless, this is an important equation because it puts constraints on
the allowed values of E, and hence we find discrete energy levels. Because
the WKB approximation is better for larger E, the constraint Eq. (35) is
expected to give correct energy levels at least for highly excited states.
Let us apply this constraint Eq. (35) for a harmonic oscillator. The
condition then is


Z xmax
1
2
2
2
2mE m x dx = n +

h.
(37)
2
xmin
The l.h.s. is an elementary integral

2mE
m 2

and we find

1
E = n+
h
.
2


(38)

This is precisely the energy levels of a harmonic oscillator. Of course, obtaining the exact energy levels is an accident for a harmonic oscillator. In
general, we expect to get approximately correct energy levels for high n. Indeed, if we apply the same constraint for the hydrogen-like atoms (for l = 0
and applying the WKB method to the radial wave function), we would obtain the correct results if we drop 1/2 (this is what BohrSommerfeld) did,
but dont if we keep 1/2 as obtained by the WKB method.
What we obtain from the WKB method is not just the energy levels,
but more details such as the approximate wave function. To simplify the
notation, let us choose a unit system such that m = = h
= 1. Then the
wave function in the region II for the harmonic oscillator is

1
cos
(2E x2 )1/4

Z

2E

2E x02 dx0

.
4


(39)

Doing the standard integral, we find

x2

1
x

. (40)
+ +
1
cos E sin1
2
1/4
(2E x )
2
2E
4
2E
2E

Because E = n + 1/2 in this unit, it further simplifies to

4.4

1
(2E x2 )1/4

 

cos E sin1
 

sin E sin1

x
2E
x
2E

+
+

x
2E

x
2E

x2
2E



x2
2E



n even
(41)
n odd

Tunneling

To study a tunneling process, we again have three regions, a classically allowed region I x < a, where the particle initially exist, a classically forbidden
region II a < x < b, and a classically allowed region III x > b where the
particle tunnels to. We follow the same matching procedure as in the bound
state example. A difference is that the matching at x = b uses a different
version of the Airy function which rises exponentially away from x = b. This
version of the Airy function Bi is explained in the note on steepest descent
method. The bottom line is that the matching at x = a is done as
1/2

'

h
(2mV 0 (a)/
h2 )1/3

2m(E V )

1 Z a q
2m(E V )dx0 +
cos
,
h
x
4


1/2

1 h
(2mV 0 (a)/
h2 )1/3
q
'
2
2m(V E)

exp

1 Z x q
2m(V E)dx0 .
h
a

(x < a)

(x > a)

The matching at x = b is, on the other hand,

1/2

i h
(2mV 0 (b)/
h2 )1/3
q
' C
2
2m(V E)

1/2

i h
(2mV 0 (b)/
h2 )1/3
' C q
2
2m(E V )

+1 Z b q
2m(V E)dx0 ,
exp
h
x
1 Z xq

sin
2m(E V )dx0 +
.
h
b
4


(x < b)

(x > b)

The overall normalization factor C is determined by the requirement that

the behavior of the wave function is consistent for a < x < b between two
9

matching procedures. We therefore find

V 0 (a)
C = i
V 0 (b)

!1/6

1 Z bq
exp
2m(V (x) E)dx .
h
a

(42)

Comparing two classically allowed regions, and taking advantage of a further

normalization change, the matching reduces to
1
1 Z a q

'
cos
2m(E V )dx0 +
,
(x < a)
1/4
(2m(E V ))
h
x
4
!
 Z xq

1 Z bq
1
1
1

0
exp
sin
.
'
2m(V (x) E)dx
2m(E V )dx +
2
h
a
(2m(E V ))1/4
h
b
4


In other words, the amplitude in the region x > b due to tunneling from the
region x < a is suppressed by
!

1
1 Z bq
exp
2m(V (x) E)dx .
2
h
a

(43)

Usually people refer to Gamovs transmission coefficient

1 Z bq
2m(V (x) E)dx
exp 2
h
a

(44)

as a suppression factor for the tunneling rate (square of the amplitude).

10