Griffiths First Part
Griffiths First Part
Quantum
Mechanics
SECOND EDITION
DAVID J. GRIFFITHS
Fundamental Equations
Schrodinger equation:
9*
ih\342\200\224 = //*
dt
Hamiltonian operator:
2m
Momentum operator:
p = -ihV
^->\342\200\242\302\253>-(\302\245:
Canonical commutator:
[x, p] = ih
Angular momentum:
Pauli matrices:
'0 P '0 -i'
a, = =
I
a~
A 0 ^=/ 0 0
Fundamental Constants
Charge
of proton: e = 1.60218 x 10\"I9C
Hydrogen Atom
Fine structureconstant: or =
1/137.036
47t\302\243ohc
a =
fl2
4jT\342\202\254o h
Bohr radius: 5.29177 x 10-11 m
mee2 amec
m ee4
Bohr energies: En
=
-4 (n= 1,2,3,...)
2(47te 0)\302\2762
9 9
h2 a~mec~
Binding energy: -EX = = 13.6057 eV
2m ea2 2
1
Ground state: e -rja
y/jta^
l
Rydbergformula: X
~
ny nj
Quantum Mechanics
Second Edition
David J. Griffiths
Reed College
^^^KiiH PearsonEducation
International
Editor-in-Chief. Science: John Challice
Senior Editor: Erik Fahlgren
Associate Editor: Christian Bolting
Editorial Assistant: Andrew Sobel
Vice President and Director of Production and Manufacturing, ESM: David W. Riccardi
Production Editor: Beth Lew
Director of Creative Services:Paul Belfanli
__ \302\251
2005, 1995 Pearson Education, Inc.
Pearson Prentice Hall
Prentice
Hall Pearson Education, Inc.
Upper Saddle River, NJ 07458
All rights reserved, No part of this book may be reproduced in any form or by any means, without
10 9 8 7
ISBN D-13-nil7S-T
If you purchased this book wilhin the United States or Canada you should be aware that il has been
wrongfully imported without the approval of the Publisher or the Author.
PREFACE vii
PARTI THEORY
1.3 Probability 5
1.4 Normalization 12
1.5 Momentum 15
1.6 The Uncertainty Principle 18
3 FORMALISM 93
3.1 Hilbert Space 93
3.2 Observables96
3.3 Eigenfunctions of a Hermiiian Operator 100
iii
3.4 Generalized Statistical Interpretation106
3.5 The
Uncertainty Principle 110
3.6 Dirac Notation 118
PART II APPLICATIONS
11 SCATTERING 394
11.1 Introduction 394
11.2 Partial Wave Analysis 399
11.3 Phase Shifts 405
408
11.4 The Born Approximation
12 AFTERWORD 420
INDEX 459
PREFACE
UnlikeNewton'smechanics, Maxwell's or
electrodynamics, or Einstein's relativity,
quantum theory was not even created\342\200\224or definitively packaged\342\200\224by one
individual, and it retains to this day some of the scars of its exhilarating
but traumatic
youth. There is no general consensus as to what its fundamental principles are, how
it should be taught, or what it really \"means.\" Every competent physicist can \"do\"
quantum mechanics,but the stories tell ourselves about what we are doing are
we
questions are saved for the end. I do not believe one can intelligently discuss what
quantum mechanics means until one has a firm sense of quantum what
realistic problems. Accordingly, this book is divided into two parts;' Part I covers
the basic theory, and Part II assembles arsenal of approximation schemes, withan
separate, it is not necessary to study the material in the order presented h ere. Some
'This structure was inspired by David Park's classic text, Introduction to the Quantum Theory,
3rd ed.. McGraw-Hill, New York (1992).
vii
instructors, for example, may wish to treat
time-independent perturbation theory
and radically counterintuitive way of thinking about the world. That, indeed, is
They
digging. may develop blisters at first, but I still is the most efficient think this
and exciting way to learn.) any rate, I can assure you that there is no deep
At
mathematics in this book, and if you run into something unfamiliar, and you don't
find my explanation adequate, by all means ask someone about it, or look it up.
There are many good booksonmathematical methods\342\200\224I
particularly recommend
Mary Boas, Mathematical Methods in the Physical Sciences, 2nd ed., Wiley, New
(No stars at all means fast food: OK if you're hungry, but not very nourishing.)
Most of the one-star problems appear at the end of the relevant section; most of
the three-starproblemsare at the end of the chapter. A solution manual is available
(to instructors only)
from the publisher.
In preparing the second edition I have tried to retain as much as possible the
spirit of the first.The only wholesale change is Chapter 3, which was much too
long and
completely
diverting; rewritten, with the background material
it has been
vectorspaces
on finite-dimensional with which most students at this level (a subject
in Chapter 2 (and fixed the awkward definition of raising and lowering operators
for the harmonic oscillator).In laterchapters I have made as few changes as I
could, even preservingthe numberingof problems equations, where possible. and
10, and a new section on partial wave phase shiftsin 11). Inevitably, the Chapter
second edition is a bit longer the first, which I regret, but I hope it is cleaner
than
and moreaccessible.
I have benefited and advice of
from the comments
colleagues, who many
read the originalmanuscript, pointed out weaknesses (or errors) the first edition, in
(Baylor), David Boness (Seattle), Burt Brody (Bard), Carter (Drew), Edward Ash
International), Robert Reynolds (Reed), Keith Riles (Michigan), Mark Semon (Bates),
HerschelSnodgrass Clark), John Taylor (Lewis and (Colorado), Stavros Theodor-
akis (Cyprus), A. S. Tremsin (Berkeley), Dan Velleman (Amherst), Nicholas
Wheeler (Reed), Scott Willenbrock (Illinois), William Wootters (Williams),Sam
Wurzel (Brown), and Jens Zorn (Michigan).
Introduction to
Quantum Mechanics
PARTI THEORY
CHAPTER 1
\342\200\2243V/3.V,
and Newton's law reads mdrxjdt1 = \342\200\224dV/dx.) This, together with
determines x{t).
9vj/ fi1 dH ,, ,
axil A]
at 2m
1
Magnetic forces are an exception, but let's not worry about ihem just yet. By the way. we shall
assume throughout this book that the motion is nonrelalivislic (,i\302\273
<SC c).
1
2 Chapter 1 The Wave Function
m
* r{X,t)
o n
X
x(t)
of a specified force.
h = =
\342\200\224
1.054572 x 10~34J s. [1.2]
2tt
But what exactly is this \"wave function,\" and what does it do for you once you've
got it? After all, a particle, by its nature, is localized at a point, whereas the wave
function (as its name suggests) is spread out in space(it'sa function of x, for any
Probability is the area under the graph of \\ty |2. For the wave function in Figure 1.2,
you
would be quite likely to find the particle in the vicinity of point A, where |4>|2
is large, and relatively mh likely to find it near point B.
^For first-hand
a delightful account of the origins of the Schrodinger equation see the article by
*M2
B C x
FIGURE 1.2: A typical wave function. The shaded area represents the probability of
rinding
the particle between a and b. The particlewouldbe relatively likely to be found
near A, and unlikely to befoundnear B.
the particle (to wit: its wave function), still cannot predict with certainty
you
the
outcome of a simple experimentto measure its
position\342\200\224all quantum
mechanics has to offer is statistical information about the possible results. This
indeterminacy has been profoundly disturbing to physicists and philosophers alike,
and it is natural to wonder whether it is a fact of nature, or a defectin the
theory.
Suppose
I do measure the position of the particle, and I find it to be at point
C.4 Question: Where was the particlejust before I made the measurement? There
are three plausible answers to thisquestion,
and they serve to characterize the main
schools of thought regardingquantum indeterminacy:
at C, and
yet quantum mechanics was unable to tellus so. To the realist,
indeterminacy
is not a fact of nature, but a reflection of our ignorance. As
d'Espagnat put
it, \"the position of the particle was never indeterminate, but was merely unknown
to the experimenter.\"5Evidently 4> is not the whole story\342\200\224some additional
4Of course, no measuring instrument is perfectly precise: what I mean is that the particle was
found in the vicinity of C. to within the tolerance of the equipment.
\342\200\242'Bernardd'Espagnat, \"The Quantum Theory and Reality\" (Scientific American. November 1979.
p. 165).
4 Chapter 1 The Wave Function
word) to worry about something that cannot, by its nature, be tested. Pauli said:
\"Oneshouldno more rack one's brain about the problem of whether somethingone
cannotknow anything about exists all the same, than about the ancientquestion of
how
many angels are able to sit on the pointofa needle.\"7 For decades this was the
\"fall-back\" position of most physicists: They'd try to sell you the orthodox answer,
but if you were persistent they'd retreatto theagnostic response,
and terminate the
conversation.
one particular number, and thereby in a sense createsthe specific result, limited
\"Quoted in a lovely article by N. David Mennin. \"Is the moon there when nobody looks?\"
(Physics Today. April 1985. p. 38).
tM2
C x
FIGURE 1.3: Collapse of the wave function: graph of |>P|2 immediately after a
measurement has found the particle at point C.
account
interpretation for the fact that the second measurement is bound to yield
the value C? Evidently the first measurement radically alters the wave function,
so that it is now sharply peaked about C (Figure1.3). We say that the wave
function
collapses, upon spike at the point C (it soonspreads
measurement, to a
out
1.3 PROBABILITY
JThe role of measurement in quantum mechanics is so critical and so bizarre that you may
well be wondering what precisely constitutes a measurement. Does it have to do with the interaction
between a microscopic(quanlum) system and a macroscopic (classical) measuring apparatus (as Bohr
insisted), or is it characterized by the leaving of a permanent \"record\" (as Heisenberg claimed), or does
it involve the intervention of a conscious \"observer\" (as Wigner proposed)? I'll return to this thorny
issue in the Afterword: for the moment let's lake the naive view: A measurement is the kind of thing
that a scientist does in the laboratory, with rulers, stopwatches, Geiger counters, and so on.
6 Chapter 1 The Wave Function
N(U)
= 1,
N(15) = 1,
N(16) = 3,
N(22)= 2,
= 2,
N(24)
N(25) = 5,
while N(17), for instance,is zero. The total number of people in the room is
00
tf = X>c/). [1.4]
7=0
(In the example, of course, N = 14.) Figure 1.4 is a histogram of the data. The
is the probability that this person's age would be 15?Answer: One chance in
14, since there are 14 possiblechoices,all equally likely, of whom only one has
that particular age. If P(j) is the probability of getting age j, then P(14) =
1/14, P(15) = 1/14,P(16)= 3/14, and so on. In general,
=
N(j)
P(j) [1.5]
N
N(J) I
j i i J I L
10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 j
FIGURE 1.4: Histogram showing the number of people, NO), with age ;, for the
distributionin Section1.3.1.
Section1.3: Probability 7
Notice that the probability of getting either 14 or 15is the sum of the individual
probabilities (in this case, 1/7). In particular, the sum of all the probabilities is
1\342\200\224you're
certain to get some age:
oo
= i. [1.6]
En/)
7=0
0\ = \342\200\224jv\342\200\224
=
1^, j nj). [1.7]
7=0
Notice that there with the average age or the medianage\342\200\224in
need not be anyone
this example nobody happens 21 or 23. In quantum mechanics to be
the average
is usually the quantity of interest; in that context it has come to be called the
expectation value. It'sa misleading term, since it suggests that this is the outcome
you would be most likely
to get if you made a single measurement {that would
be the most probable value, not the average value)\342\200\224but I'm afraid we're stuck
with it.
Question 5. What is the average of the squares of the ages? Answer: You
could get 142 = 196, with probability 1/14, or 152= 225, with
probability 1/14,
00
7=0
In general,the average value of some function of j is given by
[1.9]
8 Chapter 1 The Wave Function
N(j) A Nil) i
i i i i i i i i i i i i i i i i i i i i i
> >
123456789 10/ 123456789 10/
FIGURE 1.5: Two histograms with the same median, same average, and same most
probable value, but different standard deviations.
Aj = j-{j), [l.io]
and
compute the average of Aj. Trouble is, of course, that
you get zero, since, by
the nature of the average,Aj is as often
negative as positive:
(aj) =
J2u
- u))pu) = J2 Jpw - u) \302\243
p(j)
= U)
- (j) = o.
(Note that (j) is constant\342\200\224it does not change as you go from one member of
the sampleto another\342\200\224so it can be taken outside the summation.) To avoid this
irritating problem you might decide to average the absolute value of Aj. But
absolutevaluesare nasty to work with; instead, we get around the sign problem
by squaring before averaging:
a2 = ((Aj)2). [1.11]
Section1.3: Probability 9
root of the average of the square of the deviation average\342\200\224gulp!) is from the
called the standard deviation. The latter the customarymeasureof the spread
is
about (j).
There is a useful little theorem on variances:
a2 = ((Aj)2)
=
\302\243(A/)2P(./)
=
\302\243(./
-
(j))2P(j)
= -
\302\243
j2P(j) 2U)^jP(j) + U)2 J2 P{j)
= U2)-2U)(j) + U)2
=
U2)-U)2.
= ~
\302\260
y/(J2) U)2- [1.12]
In practice, this is a much faster way to a: Simply calculate (/2) and (j)2,
get
subtract, and take the square root. Incidentally, I warned you a moment ago that
{j~) is not, in general, equal to (j)~. Since a~ is plainly nonnegative
(from its
(r)>0*)2, [i.i3]
that can on only certain isolated values (in the example,j had
take to be an
lasted six hours, we'll take intervals of a secondor less,to be on the safe side.
Technically, Thus
we're talking about infinitesimalintervals.)
= /I P(x)dx,
Pah= P(x [1.15]
Ja
and the rules we deduced for discrete distributionstranslatein the obvious way:
p(x)dx, [1.16]
/+oo-oc
xp(x)dx, [1.17]
-oo
/+oo
f(x)p(x)dx, [1.18]
/+0C -oo
Example1.1 Suppose
I
drop a rock off a cliff of height h. As it falls. I snap a
million photographs, at random intervals. On each pictureI measure the dssUince
the rock has fallen. Question: What is the average of all these distances' That rs
Solution: The rock starts out at rest, and picks up speed as it falls; it spends more
time near the top, so the average distancemust be less than h/2. Ignoring air
resistance, the distance x at time t is
=
1 9
*(') \342\200\242
J*'
l0A statistician will complain that I am confusing the average of a finite sample (a million, in
this case) with the \"true\" average the whole continuum). This can
(over be an awkward problem for
the experimentalist, especially when the sample size is small, but here I am only concerned, of course,
with the true average, lo which the sample average is presumably a good approximation.
Section 1.3: Probability 11
P(x)\"
dt _ dx
_
1
jg dx.
T~~giy2h~27n=x
1
p(x) = , (0 < x < h)
2Vhx
1 =
1
= l.
dx
I 0 ly/hx 2Vh M
The
average distance (Equation 1.17) is
1 /2 h
r3/2
(b) Determine Ay for each j, and use Equation 1.11 to computethe standard
deviation.
(c) Use your results in (a) and (b) to check Equation 1.12.
Problem 1.2
(a) Find the standard deviation of the distribution in Example 1.1.
(b) What is the probability that a photograph, selected at random, would showa
distancex more than one standard deviation away from the average?
p(x)
= Ae-k{x~a)\\
1.4 NORMALIZATION
We return now to the statistical interpretation of the wave function (Equation 1.3),
which says that \\ty(x, t)\\2 is the probabilitydensity for
finding the particle at point
x, at time t. It follows (Equation1.16) that the integral of |^|2 must be 1 (the
particle's got to be somewhere):
= l. [1.20]
/+oo
-co
\\V(x,t)\\2dx
is called normalizing the wave function. For some solutions to the Schrodinger
How do I knowthat it will stay normalized, as time goes on, and ^ evolves? (You
can't keep ^normalizing the wave function, for then becomes a function of t, A
+0O +oc a
l*(*, t)\\2d \\V(x,t)\\2dx. [1.21]
dt;_ oo -oo dt
(Note that the integral is a function only of t, so I usea total derivative (d/dt)
in the first expression, but the integrand is a functionof x as well as r, so it's a
partial derivative (d/dt) in the secondone.)By the product rule,
9 9 9* 8V*
\342\200\224
|vl/|2
= \342\200\224
(vi/*vi/)
= ijr + vl/. [1.22]
dt dt dt dt
Now the Schrodinger equationsays that
9^ iti d2V i
[1.23]
dt 2m dx- n
9** ih d2V* / *
[1-24]
dt 2m dx1 ft
so
d f+00 w
\342\226\240 m\"> J
^ /\"t*9^ 9Vr / +00
\342\200\224
|vl/(.r,r)|~</.r
= \342\200\224 \342\200\224
* [1.26]
af / 2m (**- 3jc 3x \342\200\22400
J.qo \\
d f+\302\260\302\260
\342\200\224
\\V(x,t)\\2dx
= 0, [1.27]
dt J_oo
(d) What is the probability of finding the particle to the left of al Check your
,2A good mathematician can supply you with pathological counterexamples, but they do not arise
in physics; for us the wave function always goes to zero at infinity.
Section 1.5: Momentum 15
in x. What
which er represents the \"spread\" is the probability that the particle
would be found outsidethisrange?
1.5 MOMENTUM
For a particle in state vj>, the expectation value of x is
x\\V(x,t)\\2dx [1.28]
-00
/+oo
What exactly does this mean? It emphatically does not mean that if you measure
the position of one particle over and over again,j x\\^\\2dx
is the average of the
results you'll get. On the contrary:The first measurement (whose outcome is
indeterminate) collapse will the the value
wave functionobtained, to a spike at actually
and the subsequent measurements (if they're performed quickly) will simplyrepeat
that same result. Rather, (x) is the average of measurements performedon particles
all the state ^, which means that either you must some way of returning the
in find
wholeensembleof particles, in the same state ^, and measure the positionsof each
a shelf, each containing a particle in the state ^ (relative tothecenter the bottle). of
A graduate student with a ruler is assignedto each and at a signal they all bottle,
perfect agreement, but the more bottles use, the closer we ought to come.)) In we
To keep things from gelling too cluttered. I'll suppress the limits of integration.
The Wave Function
(I used the fact that dx/dx = 1, and threw away the boundary term, on the ground
that vj>
goes to zero at ( + ) infinity.) Performing another integration by parts, on
the second term, we conclude:
d(x)
=
/ V* \342\200\224dx [1.31]
dt m J cix
What are we to make of this result? Note that we're talking about the
\"velocity\" of the expectation value of A', which is not the same thing as the velocityof
theparticle. Nothing
we have seen so far would enable us to calculatethe velocity
It's
of a particle. not even clear what velocity means in quantum mechanics: If the
particledoesn't have a determinate position (prior to measurement), neither does it
have a well-definedvelocity. All we could reasonably probability of
ask for is the
Equation 1.31 tells us, then, how to calculate (v) directly from ^.
Actually, it is customary to work with momentum (p = mv), rather than
velocity:
[1.33]
14
The product rule says that
d.f, ,.dg df
ax ax dx
Under the integral sign. then, you can peel a derivative off one factor in a product, and slap it onto the
other one\342\200\224it'll cost you a minus sign, and you'll pick up a boundary term.
Section 1.5: Momentum 17
= [1.34]
(a-)
J ty*(x)Vd.x,
(p)= [1.35]
Jv*(j^)vdx.
We say that the operator15 x \"represents\" position, and the operator(/?//)(3/3.v)
\"represents\" momentum, in quantum mechanics; to calculate expectation values we
\"sandwich\"the appropriateoperator between ^* and ^, and integrate.
That's cute, but about other quantities? what The fact is, all classical
dynamicalvariables can be expressed in terms of position and momentum.Kineticenergy,
for example, is
l
rT = -mv 2 = P2 \342\200\224,
2 2m
L = r x m\\ = r x p
(Qix.p)) = [1.36]
J**Q(x,~^Vdx.
a2*
(T) = I1371
W*a^
Equation
1.36 is a recipe for computing the expectationvalue of any dynamical
'-''An \"operator\" is an instruction lo do something lo the function that follows it. The position
operator lells you lo multiply by .v: Ihe momentum operator tells you to differentiate with respect lo
.v (and multiply the result by
\342\200\224
ih). In this book all operators will be derivatives
(d/clt, ch/clt~.
a-/i)xciy. etc.) or multipliers (2. i. x~. etc.). or combinations of these.
18 Chapter 1 The Wave Function
Equations1.32 (or the first part of 1.33) and 1.38 are instances of Ehrenfest's
theorem, which tells us that expectation values obey classical laws.
anything, but what about quantum mechanics? Show that the wave function picks
up a time-dependent phase factor: exp(\342\200\224iV^t/h).
What effect does this have on
the expectation of a dynamical value variable?
Imagine that
you're holding rope, and you generatea
one end of a very long
wave
by shaking it up and down rhythmically (Figure 1.7). If someone asked you
\"Precisely where is that wave?\" you'd probably think he a little bit nutty: The was
(Figure 1.8), you'd get a relatively narrow bump traveling down the line. Thistime
the first
question (Where precisely is the wave?) is a sensibleone, the second and
(What is its wavelength?) seems nutty\342\200\224it isn't even vaguely periodic, so how
can you assign a wavelength to it? Of course, you can draw intermediate cases, in
which the wave is fairly well localizedand the wavelength is fairly well defined,
but there is an inescapable trade-offhere:The more precise a wave's position is,
the less preciseis its wavelength,
and vice versa.16 A theorem in Fourier analysis
makes all thisrigorous, but for the moment I am only concerned with the qualitative
argument.
That's why a piccolo player must be right on pitch, whereas a double-bass player can afford to
wear garden gloves. For the piccolo, a sixty-fourth note contains many full cycles, and the frequency
(we're working in the time domain now, instead of space) is well defined, whereas for the bass, at a
much lower register, the sixty-fourth note contains only a few cycles, and all you hear is a general sort
of \"oomph,\" with no very clear pitch.
Section 1.6: The Uncertainty Principle 19
50 x (feet)
position.
AH
* 10 20 30 40 50 x (feet)
h lit ft
P = = [1.39]
l
Thusa spreadin wavelength corresponds to a spread in momentum, and our general
observation now
says that the more precisely determined a particle's position is,
is
the lessprecisely its momentum. Quantitatively,
[1-40]
l7I\"ll prove this in due course. Many authors lake the de Broglie formula as an axiom, from
which they then deduce the association of momentum with the operator (h/i)(B/dx). Although this is
a conceplually cleaner approach, il involves diverting mathematical complications lhal I would rather
* Problem 1.9 A
particle of mass m is in the state
V(x,t) = Ae-al(mx2/ti)+i'\\
(a) Find A.
(b) For what
potential energy function V(x) does ^ satisfy the Schrodinger
equation?
(c) Calculate the
expectation values of x, x~, p, and p .
(d) Find ax and
ap.
Is their product consistent with the uncertainty principle?
1,5,9,...).
(b) What is the most probable digit? What is the median digit? What is the
average value?
bounces perfectly off the pins at either end, so that if you give it a flick it is
equally likely to cometo restat any angle between 0 and tt.
(a) What is the probability density, p(0)? Hint: p(6)d6 is the probabilitythat
the needle will come to rest between 9 and (0+d0). Graphp(0) asa function
of 6, from \342\200\224tt/2
to 3tt/2. (Of course, part of this interval is excluded,so p
is zero there.) Make sure that the total probability is 1.
Further Problems for Chapter 1 21
Problem 1.12 We consider the same device as the previous problem, but this time
we are interested in the .*-coordinate of the needle point\342\200\224that is, the \"shadow,\"
or \"projection,\" of the needle on the horizontalline.
(a) What is the probability density p(.r)? Graph p(x) as a functionof ,v, from
to
\342\200\2242r +2/-, where r is the length of the needle.Makesurethe total
probability
is 1. Hint: p(x)dx is the probability that the projection lies between
x and (x + dx). know (from You Problem 1.11) the probability that 9 is in
a given range;the questionis, what interval dx corresponds to the
interval dOl
**J\302\260\302\261
=
j(a.t)-J(b.t),
dt
where
ift ( dV
J(x.t) = \342\200\224
4/
[V^^-V*
\342\200\224
2\302\273r V
2/7/ \\ 3,v
dx dx
(b) Find the probability current for the wave function in Problem 1.9. (This is
not a very pithy example, I'm afraid; we'll encounter more substantial ones
in due course.)
The Wave Function
'+00
Pit) \\V(x,t)\\2dx = e-'tT.
/-t-oo
-00
V = Vo-iT,
dt n
(b) Solve for P(t), and find the lifetime of the particle in terms of F.
d f\302\260\302\260
\342\200\224
/ ** vj/, dx = 0
dt y_oo
for any
two (normalizable) solutions to the Schrodinger equation, ty\\ and vi/2.
f A(a2-x2]
I o. otherwise.
(c) What is the expectation value of p (at time t = 0)? (Note that you cannot
(g)
Find the uncertainty in p (crp).
(h) Check that your results are consistentwith the uncertainty principle.
p- =
\342\200\224
3
-kBT
2/77 2
_ h
v/3/77/V/jr
this problemis
The purpose of to
anticipate which systems will have to be treated
quantum mechanically,and which can safely be described classically.
(a) Solids. The lattice spacingin a typicalsolidis around d = 0.3 nm. Find the
temperature below the free18 electrons in which a solid are quantum
mechanical. Below what temperature are the nuclei in a solid quantum mechanical?
(Use sodium as a typical case.) Moral: The free electronsin a solid are
l8In a solid Ihe inner electrons are attached to a particular nucleus,and for them the relevant
size would be the radius of the atom. But the outermost electrons are not attached, and for them the
relevant distance is the lattice spacing. This problem pertains to the outer electrons.
CHAPTER 2
TIME-INDEPENDENT
SCHRODINGER EQUATION
In Chapter 1 we talked a lot about the wave function, and how you use it to
calculate various quantities of interest.Thetimehascome to stop procrastinating,
and confront what is, logically, the prior question:How do you get ^(x, t) in the
first place? We need to solve the Schrodinger equation,
9vl/ h1 92*
/ft\342\200\224=
-
\342\200\224\342\200\224+ V\302\245. [2.1]
dt 2m dx-
for a specified potential1 V(x, t). In this chapter (and most of this book) I shall
assume that V is independent of t. In that case the Schrodinger equationcan be
solved by
the method of separation of variables (the physicist's first line of attack
on any partial differential equation): We look for solutions that are simpleproducts,
'It is tiresome to keep saying \"potential energy function.\" so most people just call V the
\"potential.\" even though this invites occasional confusion with electric potential, which is actually
potential energy per unit charge.
24
Section 2.1: Stationary States 25
subsetof all solutions in this way. But hang on, because the solutions we do obtain
turn out to be of great interest. Moreover(as is typically the case with separation
of variables) we will be able at the end to patch together the separable solutions
in such a way
as to construct the most general solution.
For separable solutions we have
equationreads
(ordinary derivatives, now), and the Schrodinger
h2
d2f
2m dx-
h2 1 d2ir
\342\200\224
in + V. [2.3]
<p dt 2m \\f/
dx2
Now, the left side is a functionof t alone, and the right side is a function of
x alone.2 The only way this can possibly be true is if both sides are in fact
argument, so if it's new to you,be sure to pause and think it through.) For reasons
that will appearin a moment, we shall call the separation constant E. Then
.\342\200\236ld<p
//2--7-
= \302\243.
<p
dt
or
d<p iE
~dt=~ h
[2.4]
and
h2 1 d2xj/
+ v --= E.
2m \\j/ dx2
or
h2 d2i/
2m dx1
Ef. [2-5]
ofvariableshas
Separation turned a
partial differential equation into two
ordinary differential equations (Equations 2.4 and 2.5). The first of these(Equation
2.4)
-Note that this would not be true if V were a function of t as well as .v.
Time-Independent Schrodinger Equation
(p(t)
= e-[Et^. [2.6]
(Q(X, p))
=
^*Q L l\302\261\\^djCt [2.9]
j
Evety expectation value is constant in time; we might as well dropthe factor <p(t)
\"the wave function,\" but this is sloppy language can be dangerous, and that it is
important to remember that the true wave function always carries that exponential
time-dependent factor.) In particular, {x) is constant, and hence (Equation 1.33)
(p) = 0. Nothingever happensin a stationary
state.
P2
H(x,p) = ^-
2m
+ V(x). [2.10]
Hf = Ejfr, [2.12]
and hence
So the variance of H is
But remember, if a = 0, then every memberof the sample must share the same
value (the distribution has zero spread). Conclusion:A separable solution has the
property that every measurement of the total energyis certainto return the value
E. (That's why I chose that letter for the separation constant.)
3. The
general solution is a linear combination of separable solutions. As
we're about to discover, the time-independent Schrodinger equation (Equation2.5)
yields an infinite collection of solutions (ijfi(x), xj/jix), foix),...), each with
its associated value of the separation constant (E\\, Ei, thus there is a \302\2433,...);
4Whenever confusion might arise. I'll put a \"hat\" C) on the operator, to distinguish it from the
dynamical variable it
represents.
5A linear combination of the functions f\\ (z). /2(2) is an expression of the form
is itself a solution. Once we have found the separable solutions, then, we can
immediatelyconstructa much more general solution, of the form
CO
\302\245(*. t)
= [2.15]
J^cMx)e-',E\"'/h-
11=1
put
it into more elegant language, but the main point is this:Once solved you've
independent) potential V(x), and the starting wave function^(.v,0); your job is
to find the wave function, ^(x, t), for any subsequent time t. To do this you must
solvethe (time-dependent)
Schrodinger equation (Equation 2.1). The strategy6 is
first to solve the time-in dependent Schrodinger equation (Equation 2.5); this yields,
in general, an infinite set of solutions (\\j/\\ (x), \\j/2(x), 1^3CO,.. each
\342\226\240), with its own
associated energy (\302\2431, Ei, \302\2433,...). To fit ^(a-,0) you write down the general
linearcombination of these solutions:
00
vI/(x,0) = [2-16]
]Tc,,iMa-):
choice of the .
constants c\\, ci, C3, ... To construct W(x, t) you simply tack onto
each term its characteristic
timedependence, exp(\342\200\224/\302\243,,///2):
00 CO
-iE\342\200\236t/tt
V(x. 0-- ,if\342\200\236(x)e r%(-V 0.
\302\273=1 \302\273=1
\"Occasionally you can solve the time-dependent Schrodinger equation without recourse to
separation of variables\342\200\224-see. for instance. Problems 2.49 and 2.50. But such cases are extremely rare.
Section 2.1: Stationary States 29
are stationarystates, sense that in the all probabilities and expectation values are
independentof time, property but this is emphatically not shared by the general
solution(Equation2.17); the energies are different, for different stationary states,
and the exponentialsdo not cancel, when you calculate |^|2.
Example 2.1 Suppose a particle starts out in a linear combination of just two
stationary states:
real.) What is the wave function *I>(a\\ t) at subsequent times? Find the probability
density, and describeits motion.
Solution: The first part is easy:
vl/(A% t)
= c\\f\\{x)e-iE^fh + c2 i/o (*)<?\"'\" \302\2432'//!,
is that if you've got one that is not, it can always be expressed as a linear
combination of solutions (with the same energy) that are. So you might as
well stick to i/r's that are real. Hint: If i/(x) satisfies Equation 2.5, for a
given E, so toodoes its complex conjugate, and hence also the real linear
combinations (i/r + \\f/*) and i(if/
\342\200\224
if/*).
*Problem 2.2 Show that E must exceed the minimum value of V(x), for every
normalizable solution to the time-independent Schrodinger equation.What is the
d2i/ 2m
dxl tr
if E < Vmin, then \\j/ and its second derivative always have the same sign\342\200\224argue
Suppose
0. ifO<A'<\302\253,
V(x)
otherwise
[2.19]
oo.
bumpers\342\200\224it just keeps bouncing back and forth forever. (This potential is
artificial, of course, but I urge you to treat it with respect. Despiteits simplicity\342\200\224or
frequently.)
V(x)i
Outside the well, ij/(x) = 0 (the probabilityof finding the particle there is
zero). Inside the well, where V = 0, the time-independent Schrodinger equation
(Equation 2.5) reads
h2 d2ir
=
y Ex//, [2.20]
2m dx
or
d2^ 1^ , u r
=
^2m\302\243
mil
\342\200\224=-= \342\200\224k-d/, where /: \342\200\224-\342\200\224. [2.21]
dx1 h
where A and B are arbitrary constants. Typically, these constants are fixed by the
^(0) = \\/r(a)
= 0, [2.23]
t/K0)
= Asin0 + ficos0 = fi,
so B = 0, and hence
i/r(x)
= Asmkx. [2.24]
that
But k = 0 is no good (again, that would imply yj/(x) = 0), and the negative
solutions give nothing new, since sin(\342\200\2240)
= \342\200\224
sin(0) and we can absorb the
minus sign into A. Sothe distinct solutions are
mt
=
k\342\200\236 \342\200\224. with n = 1, 2, 3. ... [2.26]
a
Time-Independent Schrodinger Equation
A
Vi(*)| V|/2(X) \302\2453W-
FIGURE 2.2: The first three stationary states oftheinfinite square well (Equation 2.28).
[2.27]
In radical contrast to the classical case, a quantum particle in the infinite square
well cannot have
just any old energy\342\200\224it
has to be one of these special allowed
values.8To find A, we normalize \\}/:
This only determinesthe magnitudeof A, but it is simplest to pick the positive real
root: A =
*J2/a (the phase of A carries no physical significance anyway).
Inside
(X)
= J~ [2.28]
tyn
v!sin(T*)-
As
promised,
the time-independent Schrodinger equation has delivered an
infinite set of solutions (one for each positive integer n). The first few of theseare
plottedin Figure
2.2. They look just like the standing waves on a stringof lengtha;
\\J/\\, which carries the lowest energy, is called the ground state, the others,whose
energies increase in proportion to n2, are called excited states. As a collection, the
functions ^\342\200\236(x) have some interesting and important properties:
1. They are alternately even and odd, with respect to the center of the well:
t/^i is even, \\j/2 is odd, 1//3 is even, and so on.9
8Nolice lhal the quantization of energy emerged as a rather technical consequence of the
boundaryconditions on solutions to the time-independent Schrodinger equation.
9To make this symmetry more apparent, some authors center the well at the origin (running it
from \342\200\224
a to The even
+\302\253)\342\200\242 functions are then cosines, and the odd ones are sines.SeeProblem 2.36.
\\
[jfm(xrfa(x)dx
= 0, [2.29]
whenever m ^ n. Proof:
1 rr \\ fm+n
= -1 cos (m-n Ttx \342\200\224
cos Ttx
\\~|
dx
,
a Jo I V a J \\
a
/J
\342\200\224
. (m n \\ 1 . (m + n \\\\
\342\200\224\342\200\224 \342\200\224
sin Ttx sin Ttx \\
(m
\342\200\224
n)Tt
n)it \\
\\ aa )/ (in +
(in + n)Tt \\ a ) \\
sin[(w \342\200\224
n)rt] sin[(w + n)7t]
= 0.
il (m \342\200\224
n) (m + n)
Note that this argument does not work if m = n. (Can you spot the point at which
it fails?) In that case nonnalization tells us that the integral is 1. In fact, we can
combine orthogonalityand normalization
intoa single statement:10
_ J 0, if iw
#\302\253;
[231]
1. \\fm=n.
5\"\"'-j
In this case the Vr's-are real, so the * on rjrm is unnecessary, but for future purposes it's a good
idea to get in the habit of putting it there.
Time-Independent Schrodinger Equation
studied advanced calculus you will recognize Equation 2.32 is nothing but the that
Fourier series for f(x), the fact that \"any\" function can be expanded in this
and
way
is sometimes called Dirichlet's theorem.11
The coefficients can
c\342\200\236 be evaluated\342\200\224for a given f(x)\342\200\224by
a method I call
oo . oo
,\342\200\236\302\253./(*)*
= = =
cm. [2.33]
\302\243>,\342\200\236\342\200\236
/ E,/\302\253-v\302\273V,W*
/1=1 /1 = 1
(Notice how the Kronecker delta kills every term in the sum except the one for
which n = m.) Thus the nth coefficient in the expansion of f(x) is12
=
\342\226\240n / fnO x)*f(x)dx. [2.34]
These four properties are extremely powerful,and they are not peculiar to the
infinite square well. The first is true whenever the potential itself is a
symmetric
function;the secondis universal, regardless of the shape of the potential.13
Orthogonality
is also quite general\342\200\224I'll
show you the proof in Chapter 3. Completeness
holds for all the
potentials you are likely to encounter, but the proofs tend to be
nasty and laborious; I'm afraid most physicists simply assume completeness, and
hope for the best.
The stationary states (Equation 2.18) of the infinitesquarewell are evidently
''
See, for example, Mary Boas, Mathematical Methods in the Physical Sciences, 2d cd. (New
York: John Wiley, 1983), p. 313; ,/'(.v) can even have a finite number of finite discontinuities.
'-It doesn't matter whether you use m or n as the \"dummy index\" here (as long as you are
consistent on the two sides of the equation, of course):whatever letter you use. it just stands for \"any
positive integer.\"
1*
See. for example, John L. Powell and Bemd Crasemann, Quantum Mechanics (Addison-
Wesley, Reading. MA, 1961), p. 126.
Section 2.2: The Infinite Square Well 35
CO
=
^-,0) ^,,^).
/1=1
of the
The completeness \\j/'s (confirmed in this case by Dirichlet's
theorem)
guarantees that I can always express ^(.v,0) in this way, and their orthonormality
licensesthe use of Fourier's trick to determine the actual coefficients:
=
c\342\200\236 J- sin *(-v. 0) dx. [2.37]
(\342\200\224Jc)
dynamical quantities of interest, using the procedures in Chapter 1. And this same
ritual applies to any potential\342\200\224the only things that change are the functional form
of the 1^'s and the equation for the allowed energies.
Example 2.2 A
particle in the infinite square well has the initial wave function
for some constant A (see Figure 2.3). Outside the well, of course, ^ = 0. Find
V(x.t).
A ^(x, 0)
1= = |A|2
r|\302\245(jr,0)|2</jc fax2(a-x)2dx = \\A\\2^
M
Jo JO
SO
A= /?\302\260.
V a5
\342\200\236=^r,n(-v)/lo,(a_.v)rf,
a Jo V \\ a / a3 y
\342\226\240
2v/l5 T
fa . nm v
f 2 . //77T \\
= \342\200\224=\342\200\224
<:/ .v sin I \342\200\224x dx \342\200\224
I x sin I \342\200\224x dx
/
a* v I \\ J
L Jo a Jo a
_
\342\200\224
\342\200\236/ \\2a . (nit \\ (nitx/a)- 2 /tut \\
\342\200\224
2 I .v sin \342\200\224a*
J
cos
=\342\200\224 \342\200\224x
J
V/77r/ \\ a / {nit ay \\ a /
2v^l5 o3 3(/77T)2 -2
: =
C0S(/77T) + a1 COS(/77T)+ CT =- COS(0)
a- /77T (nity (nity
4</l5 r
-
-=; [COS(0) COS(/77T)]
(mty
if /7 is even.
10,8vT5/(/7tt)3. if /7 is odd.
a 7 *-? V fl /
V
x
\\izJ 773
//=1.3.5...
[2.38]
\\
^ / CO \\
(CO
oc oc
= fm(x)*fn(x)dx
E EC'\302\273C\302\273
/
m=\\ n=\\
CO CO CO
= =
J21lc\302\273>c\"8\302\273\"' J2\\c\"\\2-
lt=l 111=1 \302\273=1
(Again, the Kronecker delta picks out the term m = n in the summation over m.)
Moreover, the expectationvalue of the energy
must be
CO
[2.39]
n=\\
so
(H) = V*HVdx = H dx
J J (]Tc\342\200\236,Vo\342\200\236)* (][>\"V0,)
= =
E Ec'\302\273c\302\273\302\243\302\273 *\302\243*\302\273dx
/
E ic\302\273i2\302\243\302\273-
*
,,= 1 \\ J n=TT.5....n
The expectation value of the energy,in thisexample, is
480/?2
_ 5/32
y^ J_
7T4ma2 fr1^ n4 ma2
,,=1.3.5....
special case of the general theorem in Problem 2.2, but this time do it by explicitly
solving the Schrodinger equation, and showing that you cannot meet the boundary
conditions.)
^Problem 2.4 Calculate (x), (x~), {p), (p2), ax, and ap, for the /ith stationary state
of the infinite
square well. Check that the uncertainty principle is satisfied. Which
*Problem 2.5 A
particle
in the infinite square well has as its initial wave function
1
You can look up the series
1 7T6
J_ + + j_
+\"\"~_
lfl 3^ 56 %()
and
111 1 1 h \342\226\240
\342\226\240
jr4
= \342\200\224
\342\226\240
H 34 54 96
in math tables, under \"Sums of Reciprocal Powers\"or \"Riemann Zela Function.\"
Section 2.2: The Infinite Square Well 39
(a) Normalize ^(x,0). (That is, find A. This is very easy,if you exploit the
orthonormalityof yfri. Recall that, having
%f/\\ and normalized ^ at t = 0,
you can rest assured stays normalized\342\200\224if
that it you doubt this, check it
explicitly after doingpart (b).)
(b) Find ty(x,t) and \\ty(x,t)\\2. Express the latter as a sinusoidalfunctionof
time,as in
Example 2.1. To simplify the result, let a>
=
jr-ft/2ma-.
(c) Compute (x). Notice that it oscillates in time.What is the angular frequency
of the oscillation? What is the amplitude oftheoscillation? your (If amplitude
is greater than a/2, go directly to jail.)
(d) Compute (p). (As Peter Lorre would say, \"Do it ze kveek
vay, Johnny!\
(e) If you measured the energy of this particle, what values might you get, and
what is the probabilityof gettingeachof them? Find the expectation value
of H. How doesit comparewith E\\ and \302\2432?
Problem 2.6 Although the overall phase constant of the wave function is of no
physical significance (it cancels out whenever you calculate a measurable quantity),
the relative phase of the coefficients in Equation 2.17 does matter. For example,
suppose we change the relative phase of \\f/\\ and \\j/2 in Problem 2.5:
where (p is some constant. Find ^(x, t), \\ty(x,t)\\-, and (,r), and compareyour
results with what you got before. Study the special cases <p = tt/2 and $> = tt.
of this
(For a graphicalexploration problem
see the applet in footnote 7.)
*Problem 2.7 A
particle
in the infinite square well has the initial wave function15
Ax< -x
\302\260 - a/2<
nx v o); = f
J A(a-x). a/2<x <a.
1
is no
\342\200\242'There restriction in principle on the shape of the starting wave function, as long
as il is normalizable. In parlicular, ty(.v.O) need nol have a conlinuous derivative\342\200\224in facl, il
doesn't even have to be a conlinuous function. However,
you iry if 10 calculate (H) using
f ty(.v. 0)*//ty(.v-. 0)dx in such a case,you may encounter technical difficulties, because the second
derivative of #(x. 0) is ill-defined. Il works in Problem 2.9 because the discontinuities occur at Ihe end
points, where ihe wave function is zero anyway. In Problem 2.48 you'll see how to manage cases like
Problem 2.7.
40 Chapter2 Time-Independent Schrodinger Equation
(c) What is the probability that a measurement of the energy would yieldthe
value E\\ ?
Problem2.8 A
particle of mass m in the infinite square well (of width a) starts
(a) What is its initial wave function, vJ/(a',0)? (Assume it is real. Don't forget
to normalizeit.)
(b) What is the probability that a measurement of the energy would yieldthe
value
7T2/72/2ma2l
Problem 2.9 For the wave function in Example 2.2, find the expectation value of
H, at time t = 0, the \"oldfashioned\"way:
(H)
= I V(x. 0)*HV(x, 0) dx.
Comparethe resultobtained in
Example 2.3. Note: Because (H) is independent of
time, there is no loss of
generality in using t = 0.
d2x
Fv = i =
\342\200\224kx m --^
dt-
where
co = J- m [2.41]
V
V(x)= ^kx2;
[2.42]
V(X)A
subtract V(xo) (you can add a constant to V(a\") with impunity, since that doesn't
V(x) =
^V\"(x0)(x-x0)2,
which describessimple harmonic oscillation the point xo), with an effective
(about
(it is customaryto eliminate spring constant in the favor of the classical frequency,
using Equation2.41). we have seen, it suffices As to solve the time-independent
Schrodingerequation:
tr 1
2
2
= Ef. [2.44]
d2xj/
\342\226\240-
+ -mco x f
2m ax1 pj-
2
l6Note that V\"(.vo) > 0. since by assumption xq is a minimum. Only in the rare case V\"(.vo) = 0
is the oscillation not even approximately simple harmonic.
Time-Independent Schrodinger Equation
In the literature find two you will entirely different approaches to this problem.
The first is a straightforward \"brute force\" solution to the differential equation,
using the power series method; it has the virtue that the same strategy can be
appliedto many
other potentials (in fact, we'll to use it in Chapter 4 treat the
stage.
\342\200\224
[p2 + {mcox)2]^
= Ef. [2.45]
2/7/
where p = {h/ i)d/dx is, of course, the momentum operator. The basic idea is to
factor the Hamiltonian,
H= \342\200\224[p2 + {mcox)2]. [2.46]
2m
Here, however, it's not quite so simple,becausep x are operators, and and
operators do not, in general, commute {xp is not the same as px). Still, this does
^=^(^+^) [2.47]
(the factorin front is just there to make the final results look nicer).
Well, what is the product a-ci+7
1
a~a+ = {ip + mcox){\342\200\224ip + mcox)
2ft mco
= \342\200\224
[p2 + {mcox)2
- imco{xp - px)].
Inmco
17We'll encounter some of the same strategies in the theory of angular momentum (Chapter 4).
and the technique generalizes to a broad class of potentials in super-symmetric quantum mechanics
(see. for example. Richard W. Robinett. Quantum Mechanics. (Oxford U.P., New York. 1997). Section
14.4).
Section 2.3: The Harmonic Oscillator 43
there's
As anticipated, an extra term, involving (xp
\342\200\224
px). We call this the
commutator of X and p; it is a measure of how badly they fail to commute. In general,
the commutator of operators A and B (written with square brackets) is
[A, B] = AB
- BA. [2.48]
In this notation,
1 '
CI-CI+ = IP'
-)
+ (ma>x)~]
->
- \342\200\224
[a\\ p]. [2.49]
ol.
2ft m co 2ft
We need to figure
out the commutator of a- and p. Warning: Operators are
notoriouslyslipperyto work with in the abstract, and you are bound to make
mistakes unless you give them a \"test function,\"f(x), to act on. At the end you
can throw away the test function, and you'll beleftwith an equation involving the
=
ft d
\342\200\236
ft d
=
n df
( X\342\200\224 df \\
This lovely
and ubiquitous result is known as the canonical commutation
relation.18
1 1
a-a+ = -\342\200\224//
ftCD
+ -.
2
[2.52]
or
H = ftcD I
\342\200\224
[2.53]
a-a+
Evidently the Hamiltonian does not factor perfectly\342\200\224there's that extra \342\200\2241/2
on the
right. Notice that the ordering of a+ and a- is important here; the same argument,
with a+ on the left, yields
1
\342\200\236
1
a+a- = \342\200\224H . [2.54]
ftCD 2
In particular,
[a-,a+] = 1. [2.55]
18In a deep sense all of the mysteries of quantum mechanics can be traced to the fact that position
and momentum do not commute. Indeed, some authors take the canonical commutation relation as an
1
H = hco I
a+a- + [2.56]
In terms of a\302\261, then, the Schrodinger equation19 for the harmonic oscillator takes
the form
= [2.57]
hco
\302\261-)f Ef
(\302\253\302\261*F\302\261i)
H(a+\\j/)
= hco I
a+ci- + - 1
(a+xf/) = hco I a+a-a+ + -a+ J
if/
=
hcoa+ lci-d+ + -) \\f/
=
a+ hco I
a+a- + 1+ -1^
= a+(H + hco)\\j/
=
a+(E + ftco)\\j/ = (E + hco)(a+\\j/).
(I used Equation
2.55 to replace a~a+ by a+a- + 1, in the secondline. Notice
that whereas the ordering of a+ and a- does matter, the ordering of a\302\261
and
any constants\342\200\224such as h, co, and E\342\200\224does not; an operator commutes with any
constant.)
By the same token,a-\\f/
\342\200\224
is a solution with energy (E hco):
= hco -
\342\200\224 = hcoci- I a+ci- -
\342\200\224
I
] (fl-VO
I
H(ci-\\j/) a-a+ \\f/
= a- hco I a-a+
- 1- - =
a-(H
\342\200\224 = a-(E
\342\200\224
= (E
- hoj)(a-i/).
Here, then, is a wonderful machine for generating new solutions, with higher and
lower energies\342\200\224if
we could find one solution,
just to get started! We call a\302\261
ladder operators, because they allow us to climb up and down in energy;a+ is
the
raising operator, and a- the lowering operator. The \"ladder\" of states is
19I'm getting tired of writing ''time-independent Schrodinger equation,\" so when it's clear from
the context which one I mean. I'll just call it the \"Schrodinger equation.\"
Section 2.3: The Harmonic Oscillator 45
But wait! What if I apply the lowering operator repeatedly? Eventually I'm
going to reacha state with
energy less than zero, which (according to the general
theoremin Problem 2.2) does not exist! At some point the machine must fail.
How can that happen? We know that ci-\\j/ is a new solution to the Schrodinger
equation, but there is no guaranteethatit will be normalizable\342\200\224it might be zero,
or its square-integral might be infinite.In practice it is the former: There occurs a
\"lowest rung\" (call it \\j/o) such that
a-i,0 = 0. [2.58]
ltd \\
= 0.
, [h-r + m(0X )fo
s/lhmco \\
ax J
Time-Independent Schrodinger Equation
or
dxf/Q
mco
Xf0-
dx h
so
_ into
A
2
xj/0(x)
= Ae 2ft .
\\ =
\\A\\2 e-'\"\302\273x2/hdx
(\302\260\302\260 =
\\A\\\\[^-
J-oo V m cc
[2.59]
To determinethe energy
of this state we plug it into the Schrodinger equation (in
the form of Equation 2.57), hco(a+ci- + 1/2)^0 = Eofa, and exploitthe fact that
= 0:
a-xj/o
=
\302\243o [2.60]
\\hco.
[2.61]
where is
A\342\200\236 the normalization constant. the raising operator
By applying
(repeatedly) to t/tq, then, we can (in principle) construct all21 the stationary states of
20In the case of the harmonic oscillator it is customary, for some reason, to depart from the usual
-lNote that we obtain all the (normalizable) solutions by this procedure. For it\" there were some
other solution, we could generatefrom it a second ladder, by repeated application of the raising and
lowering operators.But the bottom rung of this new ladder would have to satisfy Equation 2.58, and
since that leads inexorably to Equation 2.59, the bottom rungs would be the same, and hence the two
ladders would in fact be identical.
Section 2.3: The Harmonic Oscillator 47
Solution: UsingEquation2.61,
'
_>mo i
1/4
( d
i/i (x) =
A\\ \\ (mo)\\
= e u'
A\\a+\\j/0 -^3-+m(i)X )\\\342\200\224T
[2.62]
/ma>\\' lima) _w<jdl.2
=AiU v\342\200\224~
so, as it happens, = 1.
A\\
I wouldn't want to calculate xf/so this way (applying the raising operator fifty
times!), but never mind: In principle Equation 2.61 does the job\342\200\224except
for the
normalization.
You can even get the normalization algebraically,but it takes some fancy
a+\\lfn
=
fn+\\,
c\342\200\236 a- fn = dn VOi\342\200\224i [2-63]
= [2.64]
f*(a\302\261g)dx (aTf)*gdx.
/00
-co \302\253/\342\200\224CO
= max
/ f*(ct\302\261g)dx / f*(Th\342\200\224+ )gdx..
J-cc y/lhmo) J-oo V
dx J
-~0f course, the integrals must exist, and this means that /(.v) and g(x) must go to zero at
\302\261oo.
Time-Independent Scbrodinger Equation
/00 -co
\302\261g)dx
y/2hm(D J-oQ l_\\
dx /OO
-oc
QED
In particular,
/\302\273CO
* = *
(fl ^/i) (A ^/i)
\302\261
^A\"
(flTfl VOi) VOi ^A' \342\200\242
/OO
\302\261
/ \302\261
-co J-oc
so
/\302\2730C
=
|T/r\342\200\236+I|2c/A\" (/1 +
1)/ ./-co l^nl2^-
J-oc /OO -co
\302\273oo /\302\273CO /\302\273CO
= = n
(\302\253_^,)*(\302\253-^)^ K|2 / |^,-i|2\302\253/jr llM2*/*.
/OO
-OO \302\253/\342\200\224oo
/
\302\253/\342\200\224OO
and hence
a+f\342\200\236
= Jn + 11/0,+1, fl-^n = Vnfn-i- [2.66]
Thus
1 1 ,
1 1 1 1
* = =
(fl+)Vo*
71\302\260^1 VJT2 V4 V4\342\200\2423
\342\200\242
2
= [2.67]
^/1 -j== (0+^^0-
{a-\\jrm)*{a-ifn)dx= / (a+a-f,\342\200\236)*$\342\200\236
dx
/OO
-oo J\342\200\224oc
CO
fmfndx.
='\"/ oo
Example 2.5 Find the expectation value of the potential energy in the /7 th state
of the harmonic oscillator.
Solution:
1 tt\\ 1 o/ 1
(V) = { -mco~x~) = -ma>~ I \\js*x-\\j/n
d: x.
2\"
h lima)
=
(C1++CI-): p /./ \342\200\224-\342\200\224(\302\253+\342\200\224fl_). [2.69]
2m (o
So
ft\302\256
f I*
(V) *: (a+)~ + (a+a-) + (fl-fl+) + (ci-Y \\j/n
dx.
-T
Time-Independent Schrodinger Equation
the same goes for (a-)2\\f/n, which is proportional to i/n-2- So those terms drop
(V)
=
ha)
\342\200\224
(n + n + 1) =
1
-ha> In +
(
- .
\\\\
J
As it
happens, the expectation value of the potential energy is exactly half the
total (the other half, of course, is kinetic). This is a peculiarity of the harmonic
oscillator,as we'll see later on.
\342\231\246Problem 2.10
(c) Check the orthogonality of t/tq, \\j/\\, and 1/0, by explicit integration. Hint: If
you exploit the even-ness and odd-ness of the functions, there is really only
one integral left to do.
\342\231\246Problem 2.11
(a) Compute (x), (/?). (.v2), and (p2), for the states\\f/o (Equation 2.59) and \\f/\\
(c) Compute (T) (the average kinetic energy) and {V) (the average potential
energy) for these stales. (No new integrationallowed!)Is theirsumwhat you
would expect?
*Problem 2.12 Find (a*),(p), (x2), (p2), and (T), for the /7th stationary state of the
harmonic oscillator,usingthemethodof Example
2.5. Check that the uncertainty
principle is satisfied.
Problem 2.13 A
particle in the harmonic oscillator potential starts out in the state
*(.r.O) = A[3^o(*)+4TM.v)].
(a) Find A.
(c) Find (x) and (p). Don't get too excited if they oscillate at the classical
frequency; what would it have been had I specified faix), instead of \\j/\\ (x)?
Check that Ehrenfest's theorem (Equation 1.38) holds for this wave function.
(d) If you measured the energy of this particle, what values might you get, and
with what probabilities?
Problem 2.14 A
particle is in the harmonic oscillator
ground state of the with
classical frequency a>, when suddenly the spring constant quadruples, so to' = 2a>,
without initially changing the wave function (of course, will now evolve vj/
differently,
because the Hamillonian has changed). What is the probability a that
measurement of the energy would still return the value ftco/21 is the probability What
2.3.2 AnalyticMethod
We return now to the Schrodinger equation for the harmonic oscillator,
ft1 d2xj/ V 1 , ,
+ -m(D2x2xj/
= Ef. [2.70]
2m dx2 2
and solveit directly, by the series method. Things look a little cleaner if we
I mm
Imco
III [2.71]
Vfiv;
in terms of the
\302\243 Schrodinger equation reads
d2yj/
i^-K [2.72]
_ 2E
K [2.73]
ftco
Our problemis to solve Equation 2.72, and in the process obtain the \"allowed\"
values of K (and henceof E).
To begin with, note that at very large \302\243
(which is to say, at very large x), \302\2432
1r(%)
% Ae~t2/2 + Be+*1/2. [2.75]
Time-Independent Schrodinger Equation
The B term is clearly not normalizable (it blows up as |jc| -> oo);the physically
^(S) -\342\226\272
( )^2/2< at large \302\243. [2.76]
in hopes that what remains, /?(\302\243), has a simpler functional form than i/r(\302\243)
itself.23
df (dh \\ _ti/2
and
dH
/^_2^ + (f._1)/V-,V2!
df- \\df- \"s^
I propose to look for solutions to Equation 2.78 in the form of power series
in \302\243.
h<&)
= fl0 + fl|? +02^2 + =
\342\200\242
\342\226\240-
J2ajtJ- [2J9]
/=0
that
-\342\200\242'Note some approximations to motivate
although we invoked Equation 2.77. what
follows is exact. of stripping off the asymptotic
The device behavior is the standard first step in
the power series method for solving differential equations-\342\200\224see, for example. Boas (footnote 11).
Chapter 12.
-4This is known as the Frobenius method for solving a differential equation. According to
Taylor's iheorem. any reasonably well-behaved function can be expressed as a power scries, so
Equation 2.79 ordinarily involves no loss of generality. For conditions on the applicability of the
method, see Boas (footnote 11) or George B. Arfken and Hans-Jurgen Weber, Mathematical Methods
oo
It follows (from the uniquenessof power series expansions25) that the coefficient
of each power of \302\243
must vanish,
(2j + 1-K)
a;+->
J
= a;.J [2.81]
(./ + 1)0+2)
=
\302\2533
(3-K)
^ =
(7-/0 = (7-/0(3-/0 \"\342\200\242\342\200\242\342\200\242
\342\200\224g\342\200\224-.. -^-^ 120 fl,<
where
hcvcn(Z) = +
\302\2530 aii1 + \302\2534?4 H
= +
---
/*odd(\302\243) fli\302\243+fl3\302\2433+tf5\302\2435
two arbitrary constants (ao and a\\)\342\200\224which is just what we would expect, for a
second-orderdifferentialequation.
However, not all the solutions so obtained are normalizable. For at very large
2
c'j+2 \302\253
-aj,
%
C
a,- .
7
072)!
normalizable solutions the power series must terminate. There must occur some
\"highest\" j (call it /?), such that the recursion formula spitsout an+2
= 0 (this will
truncate either the series /zCvcn or the series /7 odd; the one must be zero from
other
the start: ci\\
= 0 if n is even, and ciq
= 0 if n is odd). For physically acceptable
solutions,then, Equation
2.81 requires that
K = 2h + 1.
=
E\342\200\236 (n + -)h(o. for n =0.1,2 [2.83]
Thus we recover, by a completely different method, the fundamental quantization
condition we found algebraically in Equation 2.61.
It seems at first rather surprising that the quantization of energy should
emerge from a technical detail in the powerseriessolution to the Schrodinger
equation, but let's look at it from a different perspective.Equation2.70 has
solutions, of course, for any value of E (in fact, it has two linearly independent
solutions for every E). But almost all of thesesolutions blow
up exponentially at
-\"It's no surprise lhal the ill-behaved solutions are still contained in Equation 2.81: this recursion
relation is equivalent to the Schrodinger equation, so it's got to include both the asymptotic forms we
found in Equation 2.75.
27
It is possible to set this up on a computer, and discover the allowed energies\"experimentally.\"
You might call it the wag the dog method: When the tail wags, you know you've just passed over an
allowed value. See Problems2.54-2.56.
Section2.3: The Harmonic Oscillator 55
2 :
1.5 '-
>v
X 0.5
1 i 1
-4 -2 2
4$
-0.5
-1
\\
-1.5
-2
(a)
ji
1.5
/o.5 \"
1 1 >
-4 ./^1 -2 2 ^
-0.5 n
\\
-1 ':
-1.5
-2
(b)
FIGURE 2.6: Solutionsto the Schrodinger equation for (a) E = 0.49 fico, and
(b)E = 0.51/1\302\253.
-
-20? j)
\302\253y+2 fly- [2.84]
(7 + 1)(7+2)
If n = 0, there is only one term in the series (we must pick a \\ 0 to kill /20Cid,
/?()(\302\243) =flo.
and hence
= $~12
V'o(\302\243) a0e
Time-Independent Schrodmger Equation
M$) = <M$.
and hence
= a^e-^~'2
tfr,($)
h2(t)=a0(\\-2$2).
and
and so on. (Compare Problem 2.10, where this last resultwas obtained by algebraic
means.)
In general, will
//\342\200\236(\302\243)
be a polynomial of degree 11in involving
\302\243, even powers
only, if n is an even integer, and odd powers only, if n is an odd integer. Apart
from the overall factor (\302\253o or a\\) they are the so-called Hermite polynomials,
The
#\342\200\236(\302\243).29 first few of them are listed in Table 2.1. By tradition, the arbitrary
multiplicative factor is chosen so that the coefficient of the highest power of \302\243
is 2\". With this convention, the normalized30 stationary states for the harmonic
oscillatorare
1/4
fn(x) =
mco 1
-f-/2
:#\342\200\236($)<? [2.85]
They are identical (of course) to the oneswe obtained algebraically in Equation 2.67.
polynomials, H\342\200\236(%).
= 1,
H0-
=
Hy- \302\276.
== 2,
H2 4^2_
== -
#3 8\302\2763 1\302\276
= -
HA-- 16\302\2764 48\302\2762+ 12,
=
H5- 32\302\2765-160^- f120\302\276.
30I shall not work out the normalization constant here: if you are interested in knowing how it is
done, see for example Leonard Schiff, Quantum Mechanics, 3rd ed., McGraw-Hill, New York (1968).
Section 13.
Section 2.3: The Harmonic Oscillator 57
In Figure2.7(a) I have plotted \\j/n(x) for the first few /?'s. The quantum
oscillator is strikingly different from its classical counterpart\342\200\224not only are the
* *Problem 2.17 In this problem we explore some of the more useful theorems(stated
without
proof) involving Hermite polynomials.
Hn^)
= e-*\\ [2.86]
(-\\)\"e^(J^j
Use it to derive H$ and //4.
(b) The following recursion relation gives you Hn+\\ in terms of the two
preceding
Hermite polynomials:
= -
//\342\200\236+,(\302\243) 2\302\243//\342\200\236(\302\243)
2/7#\342\200\236_,(\302\243). [2.87]
Use it, together with your answer in (a), to obtain//5 and H^.
-
'The parallel is perhaps more direct if you interpret the classical distribution as an ensemble of
oscillators all with the same energy, but with random starting times.
2 Time-Independent Schrodinger Equation
/7 = 3
/7 = 2
/7=1
/7 = 0
(a)
|Vl00(x)i20.24r
0.20-
0.16-
0.12
0.06
-
0.04 -
(b)
FIGURE 2.7: (a) The first four stationary states of the harmonic oscillator. This
material is used by permission of John Wiley 8c Sons, Inc.; Stephen Gasiorowicz,
Quantum Physics, John Wiley
8c Sons, Inc., 1974. (b) Graph of |^iool2> w^h l^e
(c) If
you differentiate an nth-order polynomial, you get a polynomialof order
(n
\342\200\224
1). For the Hermite polynomials, in fact,
^=2,7//,,,^). [2.88]
2z%)\\ or, to put it another way, it is the coefficientof z\"/n\\ in the Taylor
series expansion for this function:
=
\302\253-=2+2*
5^/^). [2.89]
/7=0
We turn next to what should have been the simplest case of all: the free particle
velocity,
but in quantum mechanics the problem is surprisingly subtle and tricky. The
Schrodinger
time-independent equation reads
h2 d2i/
= Ef. [2.90]
2m dx\\
or
d~\\l/ i y/2mE
\342\200\224L=
-k2f, vj\\ierek=- . [2.91]
dx- h
vl/(x, t)
= Aeik{x-^n + Be-iHx+^n. [2.93]
Now, any function of x and t that depends on these variables in the special
combination(x + vt) (for some constant v) represents a wave of fixed profile,
travelingin the :pf-direction, at
speed v. A fixed point on the waveform (for
Schrodinger
Time-Independent Equation
Since every point on the waveform is moving along the same velocity, its with
represents a wave traveling to the right, and the second represents a (of the same wave
energy) going to the left. By the since they only differ by the sign in front of
way,
k, we
might as well write
=
, y/lmE f
k > 0 => traveling to the right, \342\200\236
__.,
k \302\261 . with .. , ,
traveling to the left.
{ , r\302\253. [2.951 J
/j [
k < 0 =>
wavelength
= 2n/\\k\\, and,
is k accordingto the de Broglie (Equation formula 1.39),
^quantum
-*'*' _
\342\200\224
~
\342\200\224
w
E
\302\253 \342\226\240
nan
l~-\"'J
On the other hand, the classical speed of a free particle with energy E is given by
E = (l/2)mir (pure kinetic,since V = 0), so
_=
2E _
= [2.98J
^classical \\/
m ^'quantum-
In the case of the free particle, then, the separablesolutionsdo not represent
physically realizable states. A free particle cannot exist in a stationary state; or,
to put
it another way, there is no such thing as a free particle a definite with
energy.
Section 2.4: The Free Particle 61
But that doesn't meanthe separable solutions are of no use to us, for they
play a mathematical role that is entirely independent of their physical
interpretation. The general solution to the time-dependent Schrodinger equation is stilla
linear combination of separable solutions (only this time it's an integral over the
continuous variable k, instead of a sum over the discrete index n):
[2.100]
(The quantity l/y/ln is factored out for convenience; what plays the role of the
coefficient in
c\342\200\236 Equation 2.17 is the combination
(\\/y/2n)<p(k)dk.) Now this
wave function can be normalized (for appropriate 0(fc)). But it necessarily carries
a range of fc's, and hence a range of energies and speeds. call it a wave We
packet.32
In the generic quantum problem, we are given ty(x, 0), and we are asked to
find y(x,t). For a free
particle the solution takes the form of Equation 2.100;
the only question is how to determine 0(fc) so as to match the initial wave
function:
0) =
ikx
* (jr. 0 (*)<?'\" dk. [2.101]
-co
[2.102]
transform of
F(k) only (theis in the sign of the exponent).Thereis, of
difference
course, some restriction on the allowable functions: The integrals have to exist}*
For our purposes this is guaranteed by the physical requirement that ^(.v. 0) itself
-^Sinusoidal waves extend out to infinity, and they are not normalizable. But superpositions of
such waves lead to interference, which allows for localization and normalizability.
\342\226\240,3The
necessary and sufficient condition on f(x) is that
/^, \\f(x)\\2dx be finite. (In that
case f^ -3C \\F(k)\\-dk is also finite, and in fact the two integrals are equal.) SecArfken (footnote 24).
Section 15.5.
Time-Independent Scbrodinger Equation
[2.103]
_
_
f A, if \342\200\224
a < x < a,
[ 0, otherwise,
pa
|vI/(a\\0)|2c/a- = \\A\\2 / dx = 2a\\A\\2 => A =
/oo-oo J\342\200\224a
11 1
1 Ca 1 ,-ikx
=
If00 sin(fctf) ;a... m2m \342\200\236
*(jc.r) \342\200\224=
/
\342\200\224^-\342\200\224e'^-^0 dk. [2.104]
Tty/la J-oo k
<t>{k)
Section 2.4: The FreeParticle 63
u|\302\245(x,0l2
V(X, 0) m
V2a
/a/n
(a) (b)
FIGURE 2.9: Example2.6, forsmalla. (a) Graph of W(x, 0). (b) Graph of 4>{k).
it's yfotf, since the fc's cancelled out (Figure 2.9(b)). This is an example of the
=
'TTsinffcfl)
*(*)
ka
(Figure 2.10(b)). This time it's got a momentum well-defined but an ill-defined
position.
Time-Independent Scbrodinger Equation
V(x. 0)
V2a
(a) (b)
FIGURE 2.10: Example 2.6, for large a. (a) Graph of *(x, 0). (b) Graph of $(&).
I return now to the paradox noted earlier: the fact that the separable solution
*I>it(..Y, f) in Equation 2.94 travels at the \"wrong\" speed for the particle it
ostensibly represents. Strictly when we discovered
speaking, the problem evaporated
that Wk is not state.
a physically Nevertheless, it is of interestto
realizable
discover how velocity is contained in the free particle
information about wave function
of the waves, can be greater less than, or equal to, the velocity of the ripples
than,
that
go to make it up. For waves on a string,the group is the same as the velocity
noticed when you toss a rock into a pond (if you concentrate ona particularripple,
you
will see it build up from the rear, move through the group, and fade
forward
away at the front, while the groupasa propagates out at half the speed). What
whole
I need to is that
show for the wave function of a free particlein mechanics quantum
the group velocity is twice the phase velocity\342\200\224just right to represent the classical
particle speed.
The problem, then, is to determine the
group velocity of a wave packet with
the general form
1 P+OO
t) = -== dk.
*(*,
s/2tt
/
J-oo
0(Jt)e'(*jr_fttf)
(In our case co = (hk2/2m), but what I have to say now applies to any kind
of wave packet, regardlessof its dispersion formula for a> as a relation\342\200\224the
function of k.) Let us assume 0(fc) is narrowly peaked about some particular
that
Since the integrand is negligible except in the vicinity of ko, may we as well
Taylor-expand the function co(k) about that point, and keep only the leading terms:
where co'0
is the derivative of co with respect to k, at the point ko.
V 2TX J-oo
At t = 0,
1 f+oa
*(x, 0) = -= / 0(\302\243o + s)e'(kQ+s)x ds,
\\I2ti -co
J\342\200\224oo
1 r+oo
vl/(.v. t) = -\342\200\224e'\"<-a>n'+M\302\273'> /
0(^o + synko+sKx-^n ds
\\j2ti J\342\200\224oo
- [2.105]
vJ/(a\\ t) = g-'^-^o)^^ co'0t.0).
Apart
from the phase factor in front (which won't affect |vj/|2in any event) the
wave packet evidently moves along at a speed (o'Q:
dco
= [2.106]
\"group -77-
not the phase velocity of the stationary states, that matches the classical particle
Velocity:
^classical\342\200\224 = [2.10oJ
ugroup 2l)phase-
Problem 2.18 Show that [Aelkx+ Be~lkx]and [C coskx + D sin kx] are equivalent
ways of writing
the same function of x, and determine the constants C and D in
termsof A and B, and vice versa. Comment: In quantum mechanics, when V = 0,
the exponentials represent traveling waves, and are most convenient in discussing
the free particle, whereas sines and cosines correspond to standing waves, which
Problem2.19 Find the probability current, J (Problem 1.14) for the free particle
wave function Equation 2.94. Which direction does the probability current flow?
f(x)
= sin(n-7rx/a) + bn cosijmx/a)}.
2^[\302\253/i
00
n=\342\200\224OG
What is in
c\342\200\236, terms of and
a\342\200\236 bnl
vl>(.T.0)
= Ae~alxl,
(d) Discuss the limiting cases {a very large, and a very small).
*Problem 2.22 The gaussian wave packet. A free particle has the initial wave
function
vl>(.v,0)= Ae-\"x\\
where A and a are constants (a is real and positive).
(a) NormalizeV(x,0).
(b) Find W(x,t). Hint: Integrals of the form
e-\302\253,x-+bx) dx
/+OC
-oo
/2a\\i/4 e-ax2/[l+ainal/m)\\
\342\200\224
vl>(.v.D= =\342\226\240
(c) Find 1^(.1-, t)\\^. Express your answer in terms of the quantity
I a
w =
Y
1 +(2hat/m)2'
Sketch |^|2 (as a function of x) at t = 0, and again for some very large t.
(e) Doesthe uncertainty principle hold? At what time t does the system come
closestto the uncertainty
limit?
68 Chapter 2 Time-Independent Schrodinger Equation
(over whereas
\302\253.), for the second it is an integral (over k). is the physical What
rise to tworather different kinds of motion.If V{x)rises than the particle's higher
total energy (E) on either side (Figure2.12(a)), the particle is \"stuck\" in the
then
other hand, E exceeds V(x) on one side(or both), the particle comes in from then
to infinity (Figure 2.12(b)). (It can't get trapped in the potentialunlessthereis some
mechanism, such as friction, to dissipate energy, but again, we're not about talking
that.) We call this a scattering state. Some potentials admit only bound states (for
instance, the harmonic oscillator); some allow only scattering states (a potential
hill no
with
dips in it, for example); some permit both kinds,depending on the
energy
of the particle.
The solutions to the Schrodingerequationcorrespond
two kinds of to
precisely
bound and scattering states. The distinction is even cleanerin the domain, quantum
because the phenomenon of tunneling (which we'll come to allows the shortly)
particle to \"leak\" through any finite potential barrier,so the thing that matters only
In \"real life\" most potentials go to zero at infinity, in which case the criterion
simplifies even further:
IE
E > 0 => scattering state.
Because the infinite square well and harmonic oscillator potentials go to infinity as
x \342\200\224>
+oo, they admit bound states only; because the free particlepotentialis zero
Section 2.5: The Delta-Function Potential 69
(a)
V(x)k
(c)
FIGURE 2.12: (a) A bound state, (b) Scattering states, (c) A classical bound state, but
a quantum scattering state.
34
If you are irritatingly observant, you may have noticed that the general theorem requiring
E > Vmjn (Problem 2.2) does not really apply to scattering states, since they are not normalizable
anyway. If this bothers you, try solving the Schrodinger equation with E < 0. for the free particle, and
Time-Independent Schrodinger Equation
8(x)
with = 1. [2.111]
^H\302\260oo.lltl]- \302\243%*)\302\253**
f(x)8(x
- a) dx = f(a) / 8(x - a) dx =
f(a). [2.113]
/+0O
-oo J\342\200\224oo
serves to \"pick out\" the value of /(.*) at the point a. course, the integral need (Of
not go from to
\342\200\224oo -f-oo: all that matters is that the domain of integrationinclude
the
point a, so a e to a + e would
\342\200\224
do, for e > 0.) any
where a is somepositiveconstant.36
This is an artificial potential, to be sure (so was
the infinite squarewell), it's but
delightfully simple to work with, and illuminates
the basic theory a minimum with of analytical clutter. The Schrodinger equation
for the delta-function reads well
h2 d2i/ Y
-
a8(x)1r = Ejr; [2.115]
2m dx2
it yieldsboth bound states (E < 0) and scattering states (E > 0).
We'lllook first at the bound states. In the region x < 0, V(x) = 0, so
2m E 2
d2f
= ron*!
zt$=k$, [2.116]
dx1
\342\200\224T
ft1
where
K = ^\342\200\224r \342\200\242 [2.117]
n
Equation 2.116 is
i/(x) = Ae-KX+BeK\\ [2.118]
F exp(-Kx) + Gexp(/c.r); this time it's the secondterm that blows up (as x ->
+00), so
yj/{x)
= Fe-KX. (x > 0). [2.120]
:%Thedelta function itself carries units of Mkngth (see Equation 2.111). so a has the dimensions
energy x length.
Time-Independent Schrodinger Equation
FIGURE 2.14: Bound state wave function for the delta-function potential
(Equation 2.122).
\\}/{x) is plotted in Figure 2.14. The second boundary condition tells us nothing;
this is (like the infinite squarewell) the exceptional casewhere V is infinite at the
join, and it's clear from the graph that this function has a kink at x = 0. Moreover,
up to this point the delta function has not come intothe story at all. Evidently the
delta function must determine the discontinuityin the derivativeof \\j/, at x = 0.
I'll show you now how this works, and as a by-product we'll see why d\\f//dx is
ordinarily continuous.
The idea is to integrate the Schrodingerequation, from
\342\200\224
e to +e, and then
take the limit as e \342\200\224>
0:
-\342\200\224
I V(x)xj,(x)dx = E I f{x)dx. [2.123]
-^jdx+l
The first integral is nothing but d\\[//dx, evaluated at the two end points;the last
integralis zero, in the limit e -\302\2730, since it's the area of a sliver with
vanishing
width and finite height. Thus
dir
A
'djj_ = lim =
-^- lim / V(x)xj/(x)dx. [2.124]
dx ~dx +\342\202\254
dx
-J
)
/r 6->o
J-\342\202\254
d\\J/ 2ma
A \342\226\240^(0). [2.125]
.dx J h2
df/dx
= -BKe~K\\ for >
(\342\200\236y 0). so d\\j//dx\\, = \342\200\224
Bk,
h-K- ma~
E = [2.127]
2m
\\i/(x)\\2 dx = 2\\B\\2
e~lKX clx = \342\200\224 = 1,
/
/+oo-oo ./0 \302\253
B = Jk = y/ma
[2.128]
Evidently the delta-function well, regardless of its \"strength\" a, has exactly one
bound state:
ma~
nx)=\"/mae-'na^t'2; E = \" [2.129]
h 2/72
What about scattering states, with E > 0? For .v < 0 the Schrodinger equation
reads
d2f 2mE
, /2 ,
~dx~5 n-
where
=
\\j2mE
k [2.130]
h
and this time we cannot rule out either term, since neither of them blows up.
Similarly, for .v > 0,
yj/ (x)
= Feikx + Ge~ikx. [2.132]
F + G= A + B. [2.133]
2ma
ik(F -G- A + B) = (A + fl), [2.134]
or, morecompactly,
met
F- G =
A(l + 2ij8)
- 5(1 - 2/\302\243). where 0 = [2.135]
Having imposed both boundary conditions,we are left with two equations
(Equations 2.133 and 2.135) in four unknowns (A, B, F, and G)\342\200\224five, if you
count k. Normalization won't help\342\200\224this
isn't a normalizable state. Perhaps we'd
better pause, then, and examinethe physicalsignificance these various of
constants. Recall that exp(ikx) gives rise (when coupled with the time-dependent
factor to
exp(\342\200\224i.Et/ft))
a wave function propagating to the right, and exp(\342\200\224ikx)
leads to a wave propagating to the left. It follows A (in Equation 2.131) is the that
amplitude of a wave
coming in from the left, B is the amplitude of a wave
to
returning
the left, F (Equation 2.132) is the amplitude of a traveling off to the wave
say, from the left. In that case the amplitude of the wave coming in from the right
will be zero:
G = 0, (for scattering from the left); [2.136]
A is the amplitude of the incident B is the wave, amplitude of the reflected wave,
and F is the amplitudeof the transmitted wave.
Solving Equations 2.133 and
2.135 for B and F, find we
ifi 1
B = A, F = [2.137]
l-ifi \\-ifi
(If you want to study scattering from the right, set A = 0; then G is the incident
amplitude, F is thereflectedamplitude, and B is the transmitted amplitude.)
Ae ikx Fe ikx
\342\226\240ikx -ikx
B& Ge
|vl> |2, so the relative27 probability that an incident particle will be reflected back is
\\B\\2 P2
& =
= 7TTT = t2-138]
T-T-^T-
\\A\\2 1+^2
R is called the reflection coefficient. (If you have a beam of particles,it tells
you
the fraction of the incoming number that will bounce back.) Meanwhile,the
is
probabilityof transmission given by the transmission coefficient
IFI2 1
T = = [2'1391
W TW
Of course, the sum of these probabilities should be 1\342\200\224and it is:
R + T=\\. [2.140]
Notice that R and T are functions of /3, and hence (Equations 2.130and 2.135)
of E:
1 1
R = = -. T = [2.141]
1 + (2h2E/ma2)' 1+ (ma2/2h2E)
seems reasonable).
This very tidy, but there is a stickymatter of principle
is all that we cannot
altogether ignore: These scattering wave functions are not normalizable, so they
don't actually represent possible particle states. But we know what the resolution to
this problem is: We must form normalizable linear combinations of the stationary
states,just as we did for the free particle\342\200\224true physical particles are represented
by the resulting wave packets. Though straightforward
in messy
principle, this is a
business in practice, and this point at it is best to turn the problem over to a
computer.38 Meanwhile, since it is impossible to create a normalizable free-particle
wave function without involving a range of energies, R and T shouldbe interpreted
as the
approximate reflection and transmission probabilities for particles in the
vicinity
of E.
-7This is not a normalizable wave function, so the absolute probability of finding the particle
al a particular location is not well defined: nevertheless, the ratio of probabilities for the incident and
reflected waves is meaningful. More on this in the next paragraph.
38Numerical studies of wave packets scattering off wells and barriers reveal extraordinarily rich
su'ucture. The classic analysis is A. Goldberg,H.M. Schey. and J. L. Schwartz, Am. J. Pliys.35, 177
(1967); more recent work can be found on the Web.
Time-Independent Schrodinger Equation
= a5(x)
_._iV(x)
x FIGURE 2.16: The delta-function barrier.
and flies off to infinity) using stationcuy states. After all, \\j/ (in Equations 2.131
and 2.132) is simply a sinusoidalfunction,
complex, time-independent, extending
suppose, the fact that by taking linear combinations of statesspread over all space,
and with
essentially trivial time dependence, we can construct wave functions that
are concentrated about a (moving) point, with quite elaborate behaviorin time (see
Problem2.43).
As
long relevant equations on the table, let's lookbriefly
as we've got the at
change the sign of a. Thiskills the bound state, of course (Problem 2.2). On the
other hand, the reflectionand transmission coefficients, which depend only on a\",
are unchanged.Strangeto say, the particle is just as likely to pass through the barrier
as to cross over the well! Classically, of course, a particle cannot make it over an
possibility that the particle will bounce back\342\200\224though I wouldn't advise driving
off a cliff in the hope that quantum mechanics will save you (see Problem 2.35).
- tt) dx.
(b) /0\302\260\302\260[cos(3x)
+ 2]8(x
(c) /+1exp(|.t|+3)5(A--2)dT.
Section2.5: The Delta-Function Potential 77
Problem 2.24 Delta functions live under integralsigns,and two expressions (D\\ (x)
and Di{x)) involving delta functionsare saidto be equalif
\342\200\242
+ 00 C+OQ
/(.t)D, (x) dx = / f(x)D2(x) dx,
-CO J\342\200\224OQ
\302\243<\302\243
^
\342\200\242Hi
(In the rare case where it actually matters, we define 0(0) to be 1/2.) Show
that
d9/dx
= 8(x).
(P1)
= Ona/h)2.
Equation 2.144 to mean the average value of the finite integral, as L \342\200\224>\342\200\242
oo).
The source of the problem is that the delta function doesn't meet the requirement
(square-integrability) for Plancherel'stheorem(see footnote 33). In spite of this,
Equation 2.144 can be extremelyuseful,if handled with care.
\342\200\224 \342\200\224
Vq, for a<x<a,
V(x)
=
0. for |.v| > a,
[2.145]
h1 d1* d2ir 9
=
E$.
\342\200\224\342\200\224\342\200\224T
or -\342\200\224=K-1r,
2m ax- dx-
where
k =
y/\342\200\2242mE
[2.146]
is real and positive. The general solution is \\j/(x) = A exp(\342\200\224kx) + B exp(/CA'), but
the first term blows up (as a\" \342\200\224>
\342\200\224oo), so the physically admissible solution (as
before\342\200\224see
Equation 2.119) is
(Equation 2.145).
Section 2.6: The FiniteSquare Well 79
In the region
\342\200\224
a < x < a, V(x) = and
\342\200\224Vo, the Schrodinger equation reads
2m dx1
- - Vof = E\\j/, or
dx
\342\200\2243-
=
-1-f,
where
I. MH. li
[2..48]
Although E is negative, for bound states, it must be greater Vq, by the than
\342\200\224
old theorem E > Vm\\n (Problem 2.2); so / is also real and positive. The general
solution is39
\\f/(x)
= C sin(/.r) -f D cos(/x), for \342\200\224
a < x < a, [2.149]
where C and D are constants. arbitrary Finally, in the region .v > a the potential
is again zero; the generalsolutionis \\jr{x)
=
Fexp(-Kx) + Gexp(/c.t), but the
second term blows up (as x 00), so \342\200\224> we are left with
ones in Problem 2.29. The cosineis even (and the sine is odd), so I'm looking for
solutionsof the form
Mx)
= D cos(/.v). for 0 < x < a. [2.151]
\\j/(\342\200\224x).
for a < 0.
/r=/tan(/fl). [2.154]
can.
\342\200\242^You if you like, write the general solution in exponential form (C'e'lx + D'e~'x). This
leads to the same final result, but since the potential is symmetric we know the solutions will be either
even or odd. and the sine/cosine notation allows us to exploit this directly.
Time-Independent Schrodinger Equation
tan z = - 1. [2.156]
yjizu/z)2
This is a transcendental equation for z (and hence for E) as a function of zo
(which is a measure of the \"size\"of the well).It can be solved numerically, using
a computer, or graphically, by plotting tanz and yJ(zo/z)2
\342\200\224
1 on the same grid,
and looking for points of intersection (seeFigure 2.18). Two limiting cases are of
special interest:
1. Wide, deep well. If zo is very large, the intersections occur just slightly
below zn
=
hjt/2, with n odd; it follows that
rrn-n-
En + V0
=
9 \342\200\242 [2.157]
2m(2a)
But E + Vb is the energy above the bottomof the well, and on the right side
we have precisely the infinite square well energies, for a well of width 2a (see
Equation 2.27)\342\200\224or rather, half of them, since this n is odd. (The other ones,of
course, come from the odd wave functions, as you'll discover in Problem 2.29.)So
the finite
square well goes over to the infinite square well, as Vq
\342\200\224>
oo; however,
for any finite Vq there are only a finite number of bound states.
2. Shallow,narrow well. As zo decreases, there are fewer and fewer bound
states, until finally (for zo < the lowest odd state disappears)only one
tt/2, where
remains. It is
interesting to
note, however, that there is always one boundstate,no
matter how \"weak\" the well becomes.
Section 2.6: TheFiniteSquare Well 81
k=?-\342\200\224. [2.159]
n
f(x)
= C sin(Zx) + D cos(/a), for (-a < x < a), [2.160]
where, as before,
, = 2^\302\261W. [2.m
n
ik[Ae~ika
- Beika] =
l[C cos(/a) + D sin(/a)] [2.164]
continuity
of f(x) at +a yields
C sin(/\302\253) + D cos(/\302\253)
= Feika. [2.165]
/[Ccos(/fl)
- Dsin(/fl)] = ikFeika. [2.166]
40We could look for even and odd functions, as we did in the case of bound states, but the
scaltering problem is inherenlly asymmetric, since the waves come in from one side only, and ihc
exponential notation (representing traveling waves) is more natural in this context.
Time-Independent Schrodinger Equation
We can use two of these to eliminateC and D, and solve the remaining two for
B and F (seeProblem2.32):
= i
.sin(2/a) 2 2
\342\200\224
B (/
\342\200\224tt-.\342\200\224 k )F. [2.167]
2kl
-2ika
F = [2.168]
cos(2/a)-/^psin(2/fl)
coefficient
The transmission (T
=
|F|2/|A|2), expressed in terms of the
original variables, is given by
V2
T~x = 1+ [2.169]
4^^0)^(1^^^)-
Notice that T = 1 (the well becomes \"transparent'\")whenever the sine is zero,
which is to say, when
where n is any integer. The energies for perfect transmission,then, are given by
n~7T~n~
E\302\273+ Vo = [2.171]
2m(2fl)2'
which happento be precisely the allowed energies for the infinite square well. T
is plottedin Figure2.19,asa function of energy.41
*Problem 2.29 Analyze the odd bound state wave functions for the finite square
well. Derive the transcendental equation for the allowedenergies, and solve it
graphically. Examine the two limiting cases. Is there always an odd bound state?
41 This remarkable has been observed in the laboratory, in the form of the Ramsauer-
phenomenon
Townsend effect. For an illuminating discussion see Richard W. Robineii, Quantum Mechanics. Oxford
U.P.. 1997. Section 12.4.1.
Section 2.6: TheFiniteSquare Well 83
Problem 2.31 The Dirac delta functioncanbe thought of as the limiting case of a
rectangle of area 1, as theheightgoes infinity and the width goes to zero. Show
to
that your answer is consistent with Equation 2.129. Also show Equation 2.169 that
Problem2.32 Derive Equations 2.167 and 2.168. Hint: Use Equations 2.165 and
2.166 to solve for C and D in terms of F:
k
C= sin(la) + i- cos(la) eikaF\\ D = cos(la) \342\200\224
i- sin(/tf) eikaF.
Plug these back into Equations 2.163and 2.164. Obtain the transmission coefficient,
and confirm Equation 2.169.
E < V0,42
... . ( 0, if jc <0.
V(X)=\\ V0. ifx>0.
(a) Calculate the reflection coefficient, for the case E < Vq, and comment on
the answer.
(c) For a potential such as this, which does not go back to zero to the right of
the barrier,the transmission coefficient is not simply |F|2/|A|2 (with A the
nV(X)
-V0
for E > Vq. Hint: You can figure it out using Equation 2.98, or\342\200\224more
elegantly,
but less informatively\342\200\224from the probability current (Problem 2.19).
What is T, for E < V01
(d) For E > Vb, calculate the transmission coefficient for the step potential, and
check that T + R = 1.
Problem 2.35 A
particle of mass m and kinetic energy E > 0 approaches an
(a) What is the probability that it will \"reflect\" back, if E = V0/3? Hint: This
is just like Problem 2.34, except that the step now goes down, insteadof up.
(b) I drew the figure so as to make you think of a car approaching a cliff, but
obviously the probabilityof \"bouncing
back\" from the edge of a cliff is far
smaller than what you got in (a)\342\200\224unless you're Bugs Bunny. Explain why
a cliff. Hint:
this potential does not correctlyrepresent In Figure 2.20 the
potential energy of the car dropsdiscontinuously .t = 0;
\342\200\224
to Vo, as it passes
potential
energy, V = 0 outside to
from around 12 MeV (million electron volts) \342\200\224
<
\342\200\224a x < +a), V(x) = oo (otherwise). Check that your allowed energies are
consistent with mine (Equation 2.27), and confirmthat your \\f/'$ can be obtained
from mine (Equation 2.28) by the substitution x \342\200\224>\342\226\240
(x + a)/2 (and appropriate
renormalization). Sketch your first three solutions,and compareFigure 2.2. Note
that the width of the well is now 2a.
Problem 2.37 A particle in the infinite squarewell (Equation 2.19) has the initial
wave function
*(.*. 0) = A
sin3(7TA/fl) (0 < .v < a).
Determine A, find W(x. t), and calculate (x), as a function of time. What is the
expectation value of the energy?Hint: sin\" 6 and cos\" 6 can be reduced, by repeated
application of the trigonometric sum formulas, to linear combinationsof sin(w#)
and cos(w0), with m = 0. 1. 2, ... . /?.
(Equation 2.19). Suddenly the well expands to twice its original size\342\200\224the right
(a) What is the most probable result? What is the probability of getting that
result?
(b) What is the next most probable result, and what is its probability?
(c) What is the expectation value of the energy? Hint: If you
find yourself
confronted with an infinite series, try
another method.
Problem 2.39
(a) Show that the wave functionof a particle in the infinite square well returns
to its original form after a quantum revival time T = Anur/nti. That is:
vl/(jr. T) = W(x. 0) for any state (not just a stationarystate).
(b) What is the classical revival time, for a particle of energy E bouncing back
(c) For what energy are the two revival times equal'?43
4- The fact that the classical and quantum revival times bear no obvious relation to one another
(and the quantum one doesn't even depend on the energy) is a curious paradox: see Daniel Slyer.
Am.J.Phys.69,56(2001).
86 Chapter
2 Time-Independent Schrodinger Equation
unco \\ _wm1.2
vI/(a,0)
= AI I -2 /\342\200\224a-I <T^A ,
\\ \"no ,.2
1 +2 / e~^x ,
(unco a)
* ^Problem 2.43 InProblem2.22 analyzed you the stationary gaussian free particle
wave packet. Now solvethe sameproblem for the traveling gaussian wave packet,
starting with the initial function wave
vl/(A-,Q) = A<rflA'V/v
[
oo, for |jc| > a.
Treat the even and odd wave functions separately. Don't bother to normalize them.
Find the allowedenergies if necessary). How (graphically,
do they compare with
the corresponding energiesin the absenceof the delta function?
Explain why the
odd solutions are not affected the delta function. Comment by on the limiting cases
a 0
\342\200\224> and a oo.
\342\200\224>
representing motion to the right, and the othermotiontotheleft. we have never But
Suppose there are two solutions, and \\j/2, with the same energy E. Multiply the
\\j/\\
Schrodingerequation for
\\f/\\ by i/^, and the Schrodinger equation for \\f/2 by \\f/\\,
and subtract, to show that {^/id^\\/dx \342\200\224
\\f/\\d\\//2/dx) is a constant. Use the fact
that for normalizablesolutions\\j/ 0 at
\342\200\224> + oo to demonstrate that this constantis
in fact zero. Conclude that xj/i is a multiple of \\j/\\,
and hence that the two solutions
are not distinct.
L) = t^(a')-) Find the stationary states (with appropriate normalization) and the
corresponding allowed energies. Note that there are two independent solutions for
allowed! Consider the \"double square well\" potential (Figure 2.21). Supposethe
44If two solutions differ only by a multiplicative constant (so that, once normalized, they differ
only by a phase factor e\"^). they represent the same physical slate, and in this sense they are not distinct
solutions. Technically, by \"distinct\" I mean \"linearly independent.\"
kV(x)
-vn
(a) Sketch the ground state wave function \\j/i and the first excited state t/t2,
(i) for the case b = 0, (ii) for b % a, and (iii) for b ^> a.
(b) Qualitatively,how do the correspondingenergies and Ei) (E\\ vary, as b goes
from 0 to oo? Sketch and Ejib) on the same graph. E\\(b)
(c) The double well is a very primitive one-dimensional model for the potential
experienced by an electron in a diatomic(the two wells represent molecule
the attractive force of the nuclei). are free to move,
If the nuclei will they
(b), does the electron tend to draw the nuclei together, or push them apart?
(Of course, there is alsothe internuclear repulsion
to consider, but that's a
separate problem.)
summing the series in Equation 2.39, but 1 you (in footnote 15) not to
warned try
it
you could have done it using integration by parts, but the Dirac delta function
affords a much cleanerway to handle such anomalies.
(a) Calculate the first derivative of ^(.v.O) (in Problem2.7),and express the
(c) Evaluate the integral f W(x.0)*HV(x,0)dx, and check that you get the
same answer as before.
Further Problems for Chapter 2 89
* * ^Problem2.49
(a) Show that
(b) Find \\^(x, t)\\2, and describe the motion of the wave packet.
(c) Compute (x) and (p), and check that Ehrenfest's theorem (Equation 1.38) is
satisfied.
V(x,t) = -aS(x-vt),
(a) Show that the time-dependent Schrodinger equation admits the exact solution
vl/(x< t)
= \"^\"\"c-'\302\253o\\x-vr\\/H2c-i\\(E+(\\/2)mv2)t-mvx\\/h
h
where E = \342\200\224ma~/2tr
is the bound-state energy of the stationary delta
function. Hint: Plugit in and check it! Use the result of Problem 2.24(b).
(b) Find the expectationvalueof the Hamiltonian in this state, and comment on
the result.
V(x) = secrr(fl.Y),
m
where a is a positiveconstant, and \"sech\" stands for the hyperbolic secant.
\\j/q(x) = A sech(rtA').
and find its energy. Normalize \\J/-q, and sketch its graph.
46This rare exampleof an exact closed-form solution to the time-dependent Schrodinger equation
was discovered by Schrodinger himself, in 1926.
Time-Independent Schrodinger Equation
\342\200\224
/ik atenh(ax)\\
fkti) = A I
,-;\342\200\242\"
I eikx,
V ik+a J
= as usual) solvesthe
(where k y/2niE/ti, Schrodingerequation for any
fk(x)
% Ae,kx. for large negative x.
This represents,then,a wave
coming in from the left with no
reflected
accompanying
wave (i.e., no term exp(\342\200\224/Aw)).
What is the asymptotic form
of fkix) at large positivex? What are R and T, for this potential?
Comment: This is a famous example of a reflectionless potential\342\200\224every incident
a pretty obvious way to arbitrary localized potentials (Figure 2.22). To the left
(Region I), V(x) = 0, so
/2mF
fix) = Ae,kx + Be~iL\\ where k = -\342\200\224-.
n
[2.173]
f(x)
= Cfix) + Dgix).
Vix) A
Aeikx peikx
Be-'kx ^/ \\ Ge~ikx
Regions II and III). Two of these can be used to eliminate C and D, and the other
two can be \"solved\" for B and F in terms of A and G:
The four coefficients S{j, which depend on k (and hence on E), constitute a 2 x 2
matrix S, called the scattering matrix (or S-matrix, for short).The S-matrixtells
you
the outgoing amplitudes (B and F) in terms of the incomingamplitudes (A
and G):
B\\ fSu Sl2\\(A\\
[2A75]
C So
In the typicalcaseof scattering from the left, G = 0, so the reflectionand
transmission coefficients are
2 |F|2
\\B\\2 I/\"!2
= |52i I2 [2.176]
G=0
For scattering from the right, A =0, and
(Equation 2.114)'.
(b) Construct the S-matrix for the finite square well (Equation2.145).Hint:This
requires no new work, if you carefully exploit the symmetry of the problem.
* * ^Problem 2.53 The transfer matrix. The S-matrix (Problem 2.52) tells you the
outgoing amplitudes (B and F) in terms of the incomingamplitudes and (A
G)\342\200\224Equation 2.175. For some purposes it is more convenient to with the work
transfer matrix, M, which gives you the amplitudes the right of the potential
to
(F and
G) in terms of those to the left and B): (A
\\Gj \\M2l
(a) Find the four elements of the Af-matrix, in terms of the elements of the
S-matrix, and vice versa.Express /?/, 7/, /?,-, and 7)- (Equations 2.176 and
2.177) in terms of elementsof the M-matrix.
^
Mh /W.
FIGURE 2.23: A
potential consisting of two isolated pieces (Problem 2.53).
(c) Construct the M-matrix for scattering from a single delta-function potential
at point a:
V (.t)
= -a8(x-a).
(d) By the method of part (b), find the M-matrix for scattering from the double
deltafunction
V(x) = -a[8(x + a) + 8(x
-
a)].
What is the transmissioncoefficient for this potential?
Problem 2.54 Find the ground state energy of the harmonicoscillator, five to
significant digits, by the \"wag-the-dog\" method. That is, solve Equation 2.72
numerically, varying K until you get a wave function goes to zero at large that In
\302\243.
u'[0] == 0), u[x], {x, 10\"8,10), -> 10000]], {x, a, b), MaxSteps
wave function does. Now K = 1.1, and note that the tail flipsover.
try
Somewhere
tighter and tighter. As you do so, you want to adjust a, b, c, and d, to zero may
Problem 2.55 Find the three excited first state energies (to five significant digits)
for the harmonicoscillator, wagging by
the dog (Problem 2.54). For the first (and
third) excitedstate will need to set you u[0] == 0. \302\253'[0]
== 1.
the infinite square well, by wagging the dog. Hint: Refer to Problem 2.54, making
appropriatechangesto the differential equation. This time the condition you are
looking for is u{\\)
= 0.
CHAPTER 3
FORMALISM
states). The purpose of this chapter is to recast the in a more powerful theory form,
with that in mind. There is not much here that is genuinely new; the idea, rather,
is to make coherentsense of what we have already discovered in particular cases.
Quantumtheory is basedon two constructs: wave functions and operators. The
state of a system is represented by
its wave function, observables are represented
by operators. Mathematically, wave functions satisfy the defining conditions for
abstract vectors, and operatorsact on them as linear transformations. So the
natural language of quantum mechanicsis linearalgebra.1
But a form of linear algebra
it is not, I suspect,which you are immediately with
ai
|or> a =
-\302\273\342\200\242
[3.1]
W/
If you have never studied linear algebra, you should read the Appendix before continuing.
93
Formalism
= + 4- -
\342\200\242
\342\226\240
+ [3.2]
(ct\\P) a\\b\\ a$b2 a*NbN.
converge, in which case the inner productdoesnotexist, any argument involving and
The collection of all functions of constitutes a vector space, but for our .v
purposes it is too
much
large. To represent a possible physical state, the wave
= 1.
\\V\\2dx
I
The set of all square-integrable functions, on a specified interval,2
2For us.the limits (a and b) will almost always be \302\261oo. but we might as well keep things more
integrable functions is only one example of a Hilberl space\342\200\224indeed, every finite-dimensional vector
space is trivially a Hilberl space. But since Lj is Ihe arena of quantum mechanics, it's what physicists
generally mean when they say \"Hilbert space.\" By Ihe way. ihe word complete here means lhal any
Cauchy sequence of functions in Hilberl space converges lo a function that is also in the space: it has no
\"holes\" in il. jusl as the set of all real numbers has no holes (by contrast, ihe space of all polynomials,
for example, like ihe set of all rational numbers, certainly does have holes in it). The completeness
of a space has nothing lo do wilh ihe completeness (same word, unfortunately) of a set of functions,
which is the properly lhal any oilier function can be expressed as a linear combination of them.
Section 3.1: Hilbert Space 95
We define the inner product of two functions, f(x) and g(x), as follows:
If / and g are both square-integrable (that is, if they are both in Hilbert space),
their inner productis to exist guaranteed (the integral in Equation 3.6 converges to
a finite number).4This from the follows integral Schwarz inequality:5
\\g{x)\\2dx. [3.7]
\\j f{xTg{x)dx\\<JJ,\\f{x)\\2dxj
You can check for yourself that Equation 3.6 satisfiesall the conditions
for an inner
= [3-8]
<*l/> (f\\8)*-
(fm\\f,,)=Sm\342\200\236. [3.10]
Finally,a set of functions is complete if any other function (in Hilbert space) can
be expressed as a linearcombination them: of
CO
f(x) = [3.11]
^2cnf\342\200\236{x).
/7=1
4ln Chapter 2 we were obligedon occasion to work with functions that were not nqrmalizable.
Such functions lie outside Hilbert space,and wc are going to have to handle them with special care, as
you will see shorlly. For the moment. 1 shall assumethat all Ihe functions we encounter are in Hilbert
space.
a proof,
\342\200\242'For see F. Riesz and B. Sz.-Nagy.Functional Analysis (Unger, New York, 1955),
Section 21. In -a finite dimensional vector space the Schwarz inequality, |(a|y3)|- < (a|a)(/0|/3), is
easy to prove (see Problem A.5). But that proof assumes the existence of the inner products, which is
precisely what we are trying to estabiisli here.
6What about a function that is zero
everywhere except at a few isolated points? The integral
(Equation 3.9) would vanish,
though the function
still evenitself does not. If this bothers you, you
should have been a math major. In physics such pathological functions do not occur, but in any case, in
Hilbert space two functions that have the same square integral are considered equivalent. Technically,
vectors in Hilbert space represent equivalence classes of functions.
96 Chapter 3 Formalism
Problem 3.1
(b) Show that the integral in Equation 3.6 satisfies the conditionsfor an inner
A.2).
product(Section
* Problem 3.2
(a) For what range of v is the function f(x) = xv in Hilbert space, on the
interval (0. 1)? Assume v is real, but not necessarilypositive.
(b) For the specific case v =
1/2, is f(x) in this Hilbert space? What about
3.2 OBSERVABLES
of
But the complexconjugate an inner product reverses the order (Equation 3.8), so
(f\\Qf)
= (Qf\\f) for all/(.v). [3.16]
(f\\Qg)= (Qf\\z)
for all fix) and all g(x). [3.17]
definition (Equation 3.16), as you prove in Problem 3.3. So use whichever will
you like. The essentialpointis a hermitian operator can be applied either to that
the member
first of an inner product or to the second, the same result, and with
tide h oo C\302\260\302\260
gdx = (pf\\g).
f\302\260\302\260 /hdf\\*
if\\pg)= r~dX
dx
= -f*g
i -do
+ [\"-f) [3.19]
J_oo i J_00 \\
i dx)
space into a function outside it (see Problem 3.2(b)). and in this case the domain of the operator may
have to be restricted.
,s
is not quite true. A.s 1 mention in Chapter
Actually, this 1. there exist pathological functions
that are but do not go to zero at infinity.
square-integrablc However, such functions do not arise in
physics, and if you are worried about it we will simply restrict the domain of our operators to exclude
them. On finite intervals, though, you really do have to be more careful with the boundary terms,
and an operator that is hermitian on (\342\200\224oo.
oo) may not be hermitian on (0. co) or (\342\200\224n, n). If you're
wondering about the infinite square well, it's safest to think of those wave functions as residing on the
infinite line\342\200\224they just happen to be zero outside (0. a).
Formalism
Notice how the complex conjugation of i compensates for the minus sign picked
up
from integration by parts\342\200\224the operator d/dx (without the i) is not hermitian,
a possible
and it does not represent observable.
*Problem 3.3 Show that if {h\\Qh) = (Qh\\h) for all functions h (in Hilbert space),
men (f\\Qg) = (Qf\\g) f\302\260r all
/ and g (i-e-> tne two definitionsof
\"hermitian\"\342\200\224Equations
3.16 and 3.17\342\200\224are equivalent). Hint: First let h = f + g, and
then let h =
f + ig.
Problem 3.4
(d) Show that the position operator (.v = x) and the hamiltonianoperator(H =
are hermitian.
3.2.2 DeterminateStates
Ordinarily,
when you measure an observable Q on an ensemble of identically
prepared systems, state ^, you do not get the same
all in the same result each
I'm talking about competent measurements, of course\342\200\224it's always possible to make a mistake.
and simply get the wrong answer, but that's not the fault of quantum mechanics.
Section 3.2: Observahles 99
-
o1 = (\302\2532
-
(<2\302\2732>
=
(*|(G <?)2*>= (\302\2532
~
q)V\\{Q
~
<?)*> = 0. [3.21]
(Of course,if every
measurement gives q, their average is also q: (Q) = q. I also
QV=qV. [3.22]
This is the eigenvalueequation for the
operator Q; ^ is an eigenfunction of Q,
and q is the
corresponding eigenvalue. Thus
Measurement
of Q on such is certain to yield the eigenvalue,q.
a state
Note that the eigenvalue is a number (not an operator or a function). You can
of//).
d
Q= i \342\200\224
< [3.25]
dq>
where 0 is the usual polar coordinatein two dimensions. (This operator might arise
in a physical contextif we were studying the bead-on-a-ring; see Problem 2.46.)
Is Q hermitian?Findits eigenfunctions and eigenvalues.
100 Chapter 3 Formalism
Solution: Here we are working with functions /(0.) on the finite interval 0 < 0 <
2tt, and stipulate that
i-^-fW)=
d<$>
\302\242/(0)- [3-27]
Equation3.26restricts the
possible values of the q:
Problem 3.6 Consider the operator Q = d2/d02, where (as in Example 3.1)
0 is the azimuthal angle in polar coordinates, and the functions are subject to
Equation 3.26. Is Q hermitian?Findits eigenfunctions
and eigenvalues. What is
the spectrum of <2? Is the spectrumdegenerate?
If the spectrumis continuous (i.e., eigenvalues fill out an entire range) then the
the eigenfunctionsare not normalizable, they do not represent possible wave and
theory (the eigenvectors of a hermitian matrix). I'll treat the discretecasefirst, and
of a hermitianoperator
Mathematically,the normalizableeigenfunctions have two
important properties:
(/1(2/) = ((2/1/)
ci(f\\f)=cj*(f\\f)
iq is a number, so it comes outside the integral, and becausethe first function
in the inner product is complex conjugated (Equation 3.6), so toois the q on
the right). But (/1/) cannot be zero (/(.y) = 0 is not a legal eigenfunction),
so q = q*, and hence q is real. QED
Theorem 2: Eigenfunctions
belonging to distinct eigenvalues are
orthogonal.
Proof: Suppose
q'{f\\g)=q*{f\\g)
(again, the inner productsexist eigenfunctions are in because the Hilbert
space by assumption). But q is real (from Theorem 1),soif q' ^ q
it must
,()U is here ihai we assume ihe eigenfunctions arc in Hilbert space\342\200\224otherwise the inner product
might not exist at all.
102 Chapter 3 Formalism
That's why
the stationary states of the infinite square well, for example,or the
harmonicoscillator, are
orthogonal\342\200\224they are eigenfunctions of the Hamiltonian
with distinct eigenvalues. But thisproperty is not peculiar to them, or even to the
Hamiltonian\342\200\224the same holds for determinate states of any observable.
Unfortunately, Theorem2 tells us
nothing about degenerate states (q' = q).
However, if more) two
eigenfunctions (or share the same eigenvalue, any
linear combination of them is itself an eigenfunction, with the same eigenvalue
(Problem 3.7(a)), use the Gram-Schmidt orthogonalization
and we can
procedure (Problem A.4) to construct orthogonal eigenfunctions witliin each degenerate
subspace.It is neveralmost
necessary to do this explicitly (thank God!), but it can
always
be done in principle. So even in the presence of degeneracythe
can
eigenfunctions be chosen to be orthogonal, and in setting up the formalismof quantum
mechanics we shall assume that this has already been done. licenses the use That
have a third fundamental property: They span the space vector can be (every
them.12
Problem 3.7
At
(a) Suppose that f{x) and g(x) are two eigenfunctions of an operator Q, with
the same eigenvalue q. Show that any linear combination of / and g is itself
an
eigenfunction of Q, with eigenvalue q.
(b) Check that f(x)
= exp(.v) and g(x) = exp(\342\200\224a)
are eigenfunctions of the
P. A. M. Dirac. 77k' Principle!; of Quantum Mechanics. Oxford University Press, New York
(1958).
'-In some specific cases completeness is provable (wc know that the stationary states of the
infinite square well, for example, arc complete,because of Dirichlef s theorem). It is a little
awkward lo call something an \"axiom\" that is provable in some cases, but 1 don't know a better way to
handle it.
Section 3.3: Eigenfunctions of a Hermitian Operator 103
Problem 3.8
(a) Check that the eigenvalues of the hermitian operator in Example 3.1 are real.
Show that the eigenfunctions (for distinct eigenvalues) are orthogonal.
may not exist. Nevertheless,there is a sense which the three essential properties in
(reality, orthogonality, and completeness) s till hold. think it's best to approach I
= [3.30]
T-^-fP(.x)
i ax pfP(x).
is
The generalsolution
= Ae*x>h.
fp{x)
/\302\2730O
fp(x)fp(x)dx
= \\A\\2 / eHP-W dx = \\A\\227rh8(p - p'). [3.31]
/OO
-OO J \342\200\224
OO
= [3.32]
fp(x) -jL=eipxlK
then
(fp'\\fP)=&(p-p'). [3-33]
which is strikingly reminiscent of true
orthonormality (Equation 3.10)\342\200\224the indices
are now continuous variables, and the Kronecker delta has becomea Dirac delta,
but otherwise it looks just the same. I'll call Equation 3.33 Dirac orthonormality.
Formalism
Most important, the eigenfunctions are complete, with thesum (in Equation3.11)
replaced by
an integral: Any (square-integrable) function f{x) can be written in
the form
1
poo
c{p) fp{x) dp
= -==
/ dp.
c{p)eiPx'h [3.34]
/OC-oo \\JlMn J-oo
trick:
/\302\2730O
= c(p)8(p
- p')dp = c(p'). [3.35]
c(p)(fP'\\fp)dp /
/OC-oo J-oc
of
The eigenfunctions momentum (Equation 3.32) are sinusoidal, with
wavelength
X = ^. [3.36]
P
This is the old de
Broglie formula (Equation 1.39), which I promised to prove at
= [3.37]
xgy(x) ygy(x).
13What about the eigenfunctions with notire&l eigenvalues? These are not merely non-
normalizable\342\200\224they actually blow up at \302\261oo. Functions in what 1 called the \"suburbs\" of Hilbert space
(the entire area is sometimes called a \"rigged Hilbert space\":see. for example,
metropolitan Leslie
Ballentine's Quantum Mechanics: A Modern Development, World Scientific, 1998) have the property
that although they have no (finite) inner product with themselves, they do admit inner products with all
members of Hilbert space. This is not true for eigenfunctions of p with nonreal eigenvalues. In
I showed
particular. that the momentum operator is hermitian for functions in Hilbert space, but the argument
depended on dropping the boundary term (in Equation 3.19). That term is still zero if g is an
of p with a real eigenvalue (as long
eigenfunction as f is in Hilbert space), but not if the eigenvalue has an
imaginary part. In this sense any complex number is an eigenvalue of the operator p, but only real numbers
are eigenvalues of the hermitian operator /3\342\200\224theothers lie outside the space over which p is hermitian.
Section 33: Eigenfunctions of a Hermitian Operator 105
= A8(x -y).
gy(x)
This time the eigenvaluehastobe real; eigenfunctions the are not square-integrable,
but again they admit Dirac orthonormality:
/\302\27300
= \\A\\2 8{x
-
y')8(x - y)dx = -
y'). [3.38]
/00
-co
gp(x)gy(x)dx
\" /
J\342\200\224oo
\\A\\28(y
If we pick A = 1, so
= [3-40]
{gy>\\gy) &(y-y')-
with
c{y)
= /(>') [3.42]
(trivial, in this case, but you can get it from Fourier's trick if you insist).
Problem 3.9
(a) Cite a Hamiltonianfrom Chapter2 (other than the harmonic oscillator) that
has only a discrete spectrum.
(c) Cite a Hamiltonian (other from Chapter 2 than the finite square well) that
has both a discrete and a continuous part
to its spectrum.
Problem 3.10 Is the ground state of the infinite square well an eigenfunction of
momentum? If so, what is its momentum? If not, why not?
In Chapter 1 I showed you how to calculate the probability that a particle would be
found in a particular location, and how to determine the expectation value of any
observable quantity. In Chapter2 learned how to find the possible
you outcomes
of an energy measurementand their probabilities.I now in a position to state am
the generalized statistical interpretation, subsumes all of this and enables which
you to figure out the possible results any measurement, and their probabilities.
of
Together with the Schrodingerequation tells you how the wave function (which
of
probability getting the particular eigenvalue qn associated with the orthonormalized
eigenfunction f\342\200\236(x) is
where =
|c\342\200\236|2. c\342\200\236 (f\342\200\236\\V). [3.43]
V{x,t) = [3-45]
J2c'iMx)-
n
In the case of continuous specira ihe collapse is lo a narrow range aboul the measured value,
depending on the precision of the measuring device.
Section 3.4: Generalized StatisticalInterpretation 107
I'll
(Forsimplicity, assume that the spectrum is discrete; it's easy to generalize this
argument
to the continuous case.) Because the eigenfunctions are orthonormal, the
coefficientsare given by Fourier's trick:15
=
-\302\273 (/i.l*>= f\342\200\2360c)*nx,t)dx. [3.46]
j
Qualitatively, tells
c\342\200\236 you \"how much is
/\342\200\236 contained in vj>,\" and given that a
measurement has to retum one of the eigenvalues of Q, it seems reasonable that
and sure enough, this follows from the normalizationof the wave function:
= =
Jl J2 c*,,c\"8\"'\" J2 c'*c\302\253= Jl ic\302\273|2- [3.48]
n' \"
=
<<2)
E>\"U2- [3.49]
Indeed,
l5Notice thai the lime dependence\342\200\224which is not al issue here\342\200\224is carried by ihe coefficients;
to make this explicit, we should really wrile c\342\200\236(.t).
l6Again, I am
scrupulously avoiding the all-loo-common assertion \"\\cn\\2 is ihe probabilily lhal
the
particle the stale f\342\200\236.\"
is in This is nonsense. The particle is in the stale *, period. Rather, \\cn |2 is
the probabilily that a measurement of Q would yield the value qn. It is true that such a measurement
will collapse ihe stale lo the eigenfunclion so one could
/\342\200\236, corrcclly say \"|c-;;l** is ihe probabilily that
a particle which is now in the stale * u'/7/ be in the stale J'\342\200\236 subsequent lo a measurement of Q\" ...
but that's a completely different assertion.
108 Chapter 3 Formalism
but Qf\342\200\236 so
=q\342\200\236f\342\200\236,
(Q) = = =
[3.51]
^2^2c*,cHq\342\200\236(fn'\\f\342\200\236) ^2^2c*,cHq\342\200\2368\342\200\236'\342\200\236 J^<?/i|cw|2.
1 f\302\260\302\260
c(p)
=
(fpm
= -==
/ e-\"}X/h*(x.t)dx. [3.53]
Thisis such an important quantity that we give it a special name and symbol:the
momentumspace wave function, <\302\243(/?, r). It is essentially the Fourier transform
of the (position space) wave function *(x,f)\342\200\224which, by Plancherel's theorem,
is its inverse Fourier transform:
v!/(a\\ t)
= \342\200\224==
e'px/n<$>{p, t)dp. [3.55]
/
V2tt/7 7-00
of
measurement momentum would yield a result in the range dp is
Mp,t)\\2dp. [3.56]
Example
3.4 A
particle of mass m is bound in the delta function well V(x) =
(where E = \342\200\224ma2/2fi2).
The momentum space wave function is therefore
1 Vma^-iEt/h 'E\"
H ~e
4>0>,r) =
f\302\260\302\260
e-iPx/ne-ma\\x\\/n2 dx = y-
J-z-co * P2 + Pq
(I looked up
the integral). The probability, then, is
oo
2 I*00 i . -
i PPO -l
3
= + tan
p2 + pI
\342\200\224i^dp
(pl + ps> \\P0/
A\302\2730
1 1
\"
_ = 0.0908
4 2Jr
(again, I looked up
the integral).
Problem 3.11 Find the momentum-space wave function, <\302\243(/?, f), for a particle in
the ground stateof the harmonic oscillator. is What the probability (to 2
significant
digits) that a measurement of p on a particle in this state would a value
yield
outside the classical range (for the same energy)?
Hint: Look in a math table
under \"NormalDistribution\" or Function\" \"Error for the numerical part\342\200\224or
use
Mathematica.
n a
(x) = I<S>*(-^iz\\<S>dp. [3.57]
^P,
in position space;
/**fiRs)*'x'
(Q(x.p))= \342\200\242
h 9
[3.58]
a\\ = ((A
- (A))*\\(A
-
(A)m = (/|/>,
where / = (A
\342\200\224
(A))^. Likewise, for any other observable, B,
<7B
=
(g\\g)' where 8 = & - WW-
Therefore (invoking
the Schwarz inequality, Equation 3.7),
=
(f\\f)(g\\g)
> \\(f\\g)\\2- [3-59]
\302\26014
n2
\\z\\2
=
[Re(z)]2 + [Im(z)]2 > [Im(z)]2 = [3.60]
2i
Therefore, letting z
=
(f\\g),
o\\a\\ >
- . [3.61]
(J;[(f\\g) (g\\f)]j
But
(f\\g)
= ((A - (A))V\\(B
-
(B)m = (vl/|(A - (A))(B - <*))\302\245)
=
{V\\(AB
- A(B) - B(A) + (A)(B)m
= (V\\ABV)
-
(B)(V\\AV)
- (A)(V\\BV) + {A){B)-(V\\V)
= + (A)(B)
(AB)-(B)(A)-(A)(B)
= (AB)-(A)(B).
Similarly,
(g\\f) = (BA)-(A)(B),
Section 3.5: The Uncertainty Principle 111
so
a a a a a a
</!*)
\"
(g\\f) = (AB) - (BA) = {[A, B])t
where
of the
is the commutator two
operators (Equation 2.48). Conclusion:
1
[3.62]
<rlcri>[-([lB])J
(:
[x. p] = ih.
So
> = ,
\302\260l*\\
(\302\261itl) g)
or, since standarddeviationsare by their nature positive,
VxOp
> -. [3.63]
That's the originalHeisenberg uncertainty principle, but we now see that it is just
one applicationof a more general theorem. much
l7More precisely,the commutator of two hermitian operators is itself anti-hermitian (Q* = \342\200\224Q),
l8This corresponds to the fact that noncommuting matrices cannot be simultaneously diagonalized
(that is, they both be brought
cannot to diagonal form by the same similarity transformation), whereas
For example, in the hydrogenatom (as we shall see in Chapter 4) the Hamilto-
nian, the magnitudeof the momentum, angular
and the z component of angular
momentumare mutually compatible observables, and we will construct
simultaneous
eigenfunctions of all three, labeled by their respective eigenvalues.Butthere is
no eigenfunction of
position that is also an eigcnfunction of momentum, because
these operatorsareincompatible.
Note that the uncertainty principle is not an extra assumption in quantum
the
particle no longer has the position you got in the first measurement.19 The
problem, then, is that the second measurement rendersthe outcomeof the first
disturbing the state of the particle (the second collapse wouldn'tchange anything,
in that case). But this is only possible, in general, if the observables two are
compatible.
*Problem 3.13
[x\".p]
= ihtvcn-1.
IJNiels Bohr was at pains to track down the mechanism by which the measurement of .v (for
instance) destroys the previously existing value of p. The crux of the mailer is lhal in order to determine
the position of a particle you have to poke it with something\342\200\224shine light on it. say. Bui these photons
imparl to the particle a momentum you cannot control. You now know the position, but you no longer
know ihe momentum. His famous debates with Einstein include many delightful examples, showing
in detail how experimental constraints enforce the uncertainly principle. For an inspired account see
Bohr's article in Albert Einstein: Philosopher-Scientist, edited by P. A. Schilpp, Tudor, New York
(1949).
Section 3.5: The Uncertainty Principle 113
*Problem 3.14 Provethe famous \"(your name) uncertainty principle,\" relating the
uncertainty in position (A
= x) to the uncertainty in energy (B = p2/2m+ V):
h
oxoh >\342\200\224\\(p)l
2/77
Meanwhile, in Equation 3.60 I threw away the real part of z\\ equality results if
* = ia{x~ [3-67]
(7 Tx ~{p))
<jr))**
examples
we encountered earlier were gaussians.20
-^Note thai it is only the dependence of * on .v that is at issue here\342\200\224the \"constants\" -4. a. (x).
and (p) may all be functions of lime, and for thai mailer * may evolve away from ihe minimal form.
All I'm asserting is lhat if, at some instant, the wave function is gaussian in .v. then (al thai instant) the
h
AxAp>-\\ [3.69]
AtAE>-. [3.70]
principle, and in the course of that derivation to persuade you that it is really an
altogether differentbeast,whosesuperficial resemblance to the position-momentum
uncertainty principle is actually quite misleading.
After all, position, momentum, and energy are all dynamical variables\342\200\224
measurable characteristics of the system, at any given time. But time itself is
not a
dynamical variable (not, at any rate, in a nonrelativistictheory): You don't
go out and measure the \"time\"of a particle,as you might its position or its energy.
Time is the independentvariable,of the dynamical quantities are which
In
functions.
particular, the At in the energy-time uncertainty principle not the
is
standarddeviation of a collection of time measurements; roughly speaking (I'll
make thismoreprecise a moment) it is the time it takes the
in to change system
substantially.
\"'Many casual applications of the uncertainly principle are actually based (often inadvertently) on
a completely different\342\200\224and sometimes quite unjustified\342\200\224measure of \"uncertainty.\" Conversely, some
perfectly rigorous arguments use other definitions of \"uncertainly.'' See Jan Hilgevoord. Am. J. Phxs.
70, 983 (2002).
Section 3.5: The Uncertainty Principle 115
d
=
d * /9* dQ
IT QV)+[V
V + *
*<G> *<*'\302\253*>= dt
Now, the Schrodinger equation says
9*
ih \342\200\224
= HV
dt
d 1-* 1 * \342\200\236 I dQ
-r(Q) = -tt{HV\\QV) + -mQHV) + -f-
dt in in \\ at
a a
But H is hermitian, so (HV\\QV) = {y\\HQV), and hence
[3.71]
This is an interesting and useful result in its own right (see Problems3.17 and
3.31). In the typical case where the operator does not dependexplicitly on time,22
Now, suppose we
pick A = H and B = Q, in the generalized uncertainty
principle (Equation 3.62), and assume that Q does not depend explicitlyon t:
n>
\342\200\224
aHaQ
(W-tef^HfH\342\204\242)2-
\342\200\224Operators that depend explicitly on t are quite rare, so almost always 'dQ/'dt = 0. As an
example of explicit lime dependence, consider the potcnlial energy of a harmonic oscillator whose
spring constant is changing (perhaps the temperature is rising, so the spring becomes more flexible):
Q = (l/2)w[ft>(/)]V.
Formalism
Then
~
AEAt>~,
2 [3.74]
and that's the energy-time uncertainty principle. But notice what is meant by Af,
here: Since
d(Q)
OQ = Ar,
dt
Af
represents the amount of time expectation value of Q to change
it takes the by
very gradual; or, to put it the otherway around, if any observable changes rapidly,
the \"uncertainty\"in the energy must be large.
Example 3.5 In the extreme case of a stationary state, for which the energy is
uniquely determined, all expectationvaluesare constant in (AE
= 0 time =>\342\226\240
states\342\200\224say:
E2 -
2ah^ {x)in{x)cos(-
E\\
|*(*, t)\\2
= a2d/r, (.0)2 + b2tyi(x))2 + t
AE At = 27ih,
Example3.6 How
long does it take a free-particle wave packet to passby a
particular
point (Figure 3.1)? Qualitatively (an exact version is exploredin Problem 3.19),
Ar =
Ax/v
= mAx/p, but E = p2/2m, so AE =
pAp/m. Therefore,
pAp mAx
AE At = -\342\200\224 =
AxAp,
m p
which is >ft/2 by
the position-momentum uncertainty principle.
is
\342\200\242^This sometimes called the \"Mandelstam-Tamm\" formulation of the energy-time uncertainty
principle. For a review of alternative approaches see Paul Busch, Found. Phys. 20, 1 (1990).
Section3.5: The
Uncertainty Principle 117
Ax-
a kind of bell-shaped curve centered at 1232 MeV/c2, a width of about 120 with
MeV/c2 (Figure 3.2). does the rest energy (mc2) sometimes come out higher
Why
than 1232, and sometimes lower? Is this experimental error? No, for
whereas h/2
= 3 x 10-22 MeV sec. So the spreadin m is about as small as the
uncertainty principle allows\342\200\224a particle with so short a lifetime just doesn't have
a very
well-defined mass.24
24
Actually, Example 3.7 is a bit of a fraud. You can't measure 10-23 secon a stop-watch, and in
practice the lifetime of such a short-lived particle is inferred from the width of the mass plot, usingthe
uncertainty principle as input. However, the point is valid, even if the logic is backwards. Moreover,
if you assume the A is about the same size as a proton (~10-'*5 m), then 10--3 sec is roughly the
time it takes light to cross the particle, and it's hard to imagine that the lifetime could be much less
than that.
118 Chapter 3 Formalism
particle. In every case, however, At is the time it takes forthe systemto undergo
\"substantial\" change.
It is often said that the uncertainty principle means energy is not strictly
Problem 3.18 Testthe energy-time uncertainty principle for the wave function in
Problem 2.5 and the observablex, calculating by 07/, ax, and d{x)/dt exactly.
Problem 3.19 Test the energy-timeuncertainty principle for the free particle wave
packet in Problem 2.43 and the observable
\342\200\236v,
by calculating an, ax, and d(x)/dt
exactly.
Problem 3.20 Show that the energy-time uncertainty principle reduces to the \"your
name\" uncertainty principle (Problem 3.14), when the observable in question is x.
3.6 DIRACNOTATION
Imagine
an ordinary vector A in two dimensions (Figure 3.3(a)). How wouldyou
describethis vector to someone? The most convenient way is to set up
cartesian axes, x and >\\ and specify the components of A: Ax = f
\342\200\242
A, = ./
\342\200\242
A
Ay
(Figure 3.3(b)). Of course, your sister might have drawn a different set of axes,
x' and y\\
and she would report different components: Ax = V
\342\226\240
A. A'y
= j' \342\226\240
A
(Figure 3.3(c)). But it's all the same vector\342\200\224we're simply expressing it with
respect to two different bases ({/\\ j} and {?',./'}). The vector itself lives \"out there
(with \\x) standing for the eigenfunction of .v with eigenvalue .v),25 whereas the
momentum
space wavefunction <t>(p.r) is the expansion of \\S) in the basis of
momentum eigenfunctions:
<5>{p,t)
= (p\\S{t)) [3.76]
(with \\n) standing for the /7th eigenfunction of H)\342\200\224Equation 3.46. But it's all the
same state; the functionsW and \302\242,and the collection of coefficients {c,,}, contain
exactly the same information\342\200\224they
are simply three different ways of describing
the same vector:
I
9ipx/h
\302\245(*./>
=
t)S(x-y)dy $>(p.t) dp
/#(>'. -I Jhih
= [3.78]
J2^e-iE\"r/h^n(x).
-5I don't want to call it gx (Equation 3.39). because that is its form in the position basis, and
the whole point here is to free ourselves from any particular basis. Indeed, when I first defined Hilbert
space as the set of squarc-intcgrable functions\342\200\224over x\342\200\224that was already too restrictive, committing
us to a specificrepresentation (the position basis). 1 want now to think of it as an abstract vector space,
whose members can be expressed with respect to any basis you like.
26In it would be
position space /;>(.v) (Equation 3.32).
Chapter 3 Formalism
=
\\fi) Q\\a). [3.79]
Just as vectors are represented, with respect to a particular basis {\\e\342\200\236)},21 by their
components,
elements28
{em\\Q\\en)=Q,m. [3.81]
J>.ki.) = J>\302\253fik,i).
[3.82]
n n
=
\302\243M*\302\273\302\273k\302\253) \302\243fl\302\253<*i\302\253lfiki.).
[3.83]
n n
and hence
\302\273
Example3.8 Imagine
a system in which there are just two linearly independent
states:29
|1>= and |2)
=
Q (\302\260).
Til assume llic basis is discrete: otherwise n becomes a continuous index and the sums arc
replacedby integrals.
-sThis terminology is inspired, obviously, by the linile-dimcnsional case, but the \"matrix\" will
now typically have an infinite (maybe even uncountable) number of elements.
2
Technically, the \"equals\" signs here mean \"is represented by.\" but I don't think any confusion
will arise if we adopt the customary informal notation.
Section 3.6: Dirac Notation 121
-c o-
where g and h are real constants. If the system starts out (at t = 0) in state |1),
what is its state at time tl
ih \342\200\224
\\S)
= H\\S). [3.85]
dt
As always, we begin by solving the time-w dependent Schrodinger equation:
that is, we look for the eigenvectors and eigenvaluesof H. The characteristic
a (h-E
det . = (/7 - E)2 - g2
= 0 h
=\302\273
- E = T8 => E\302\261 = h \302\261
g.
\302\243
8 h E\\
\\
= {h\302\2618) =* ha +sP =
= {h\302\261g)a^f} a.
\302\261
C 0 CO CO
1 ( 1 \\
=
l*\302\261>
J2\\\302\261K
K(0)) = = + k-)).
Q -^(k+)
Finally,
we tack on the standard time-dependence exp(\342\200\224iE\342\200\236t/h):
=
\\4(t)) -l=[e-i{h+8)t'n\\4.+)+e-i(h-g)t/n\\4-)]
V2
122 Chapter 3 Formalism
2 c-Whl^+c\302\256/*
= le-iht/h fe-.s.,.>+e>-8'/n\\_-
-igt/h
~ e h cos(gtfh)
e-igt/H eigtltt J -i sin(gf//J)
If you doubt this result, by all means check it: Doesit satisfy the
time-dependent
Dirac proposedto chop the bracket notation for the inner product, (a|/3), into
two pieces, he called bra,
which (a|, and ket, \\/3) (I don't know what happened to
the c). The latter is a vector, but what exactly is the former? It's a linear function
of vectors,in the sense that when its right) it yieldsa (complex)
it hits a vector (to
number\342\200\224the inner product. (When an operator hits a vector, it delivers another
vector; when a bra hits a vector, it delivers a number.) In a functionspace,the bra
(/1
-/' [\342\200\242\342\200\242\342\200\242]dx,
(a\\
=
{a*a*_ ... a*) [3.88]
The collection of all bras constitutes another vector space\342\200\224the
so-called dual
space.
The license to treat bras as separateentitiesin their own right allows for
some powerful and pretty notation(though I shall not exploit it in this book). For
example, if |ck)
is a normalized vector, the operator
P = \\a)(ct\\ [3.89]
Section 3.6: Dirac Notation 123
picks out the portionof any other vector that \"lies along\" |a):
P\\fi) = (ct\\P)\\a);
(\302\253\302\253k\302\273)-=8\302\273\302\273. [3.90]
then
= l
\302\243h,)<e\342\200\236l
[3.91]
n
(the identity operator). For if we let this operator act on any vector |a), we recover
the expansion of \\a) in the {\\en)} basis:
= \\a). [3.92]
J2\\en)(e\342\200\236\\a)
n
(e-\\ez>)=8(z-z'), [3.93]
then
\\ez)(ez\\dz=\\. [3.94]
J
|a)=i|l)-2|2)-/|3), |/J)=/|1)+2|3).
(a) Construct {a\\ and (/31 (in terms of the dual basis (1|, (2|, (3|).
(b) Find (a\\P) and (0\\a), and confirm that (0\\a) = (a\\0)*.
(c) Find all nine matrix elements of the operator A = \\a)(fi\\, in this basis, and
construct the matrix A. Is it hermitian?
124 Chapter 3 Formalism
energy. Find its eigenvalues and eigenvectors (as linearcombinations 11) of and
(ii
= 1,2.3,...).
Q\\e\342\200\236)=q\342\200\236\\e\342\200\236)
Q\\a) = ]T//Hk\342\200\236
)(<?\342\200\236! I\302\253>,
for any
vector \\a).
(Problem A.4) to orthonormalize the functions 1, .v, x2, and x3, on the interval <
\342\200\2241
x < 1. You may recognize the results\342\200\224they are (apart from the normalization)30
Legendre polynomials (Table 4.1).
(b) Show that the commutatorof two hermitian operators is anti-hermitian. How
about the commutator of two fl/m-hermitian operators?
?'Legendre didn'l know what the best convention would be: he pickedthe overall factor so that
all his functions would ao to 1 at x = 1. and we're stuck with his unfortunate choice.
Further Problems for Chapter 3 125
eigenstates 0| and 02, with eigenvalues h\\ and bi. The eigenstates are related by
(a) Observable A is measured, and the value a\\ is obtained. What is the state of
the system (immediately)after thismeasurement?
(b)
If B is now measured, what are the possible results, and what are their
probabilities?
**Problem3.28 Find the momentum-space wave function \302\242,,(p. t) for the nth
stationary state of the infinite square well. Graph \\<&\\(p,t)\\2and \\<&2(p< t)\\2, as
to calculate the expectation value of p2. Compare your answer to Problem 2.4.
1
_ il7txl'\\
\342\226\240.e -nX < x < nX.
q/(A\\0)= \\/2n~A.
0, otherwise,
where n is some positive integer. This functionis purely sinusoidal (with wavelength
X) on the interval \342\200\224nX< x < nX, but it still carries a range of momenta,because
the oscillations do not continue out to infinity. Find the momentum space wave
function \302\242(/7,0). Sketch the graphs of |^(a-,0)|2 and |\302\253J>(/?, 0)|2, and determine
their widths, wx and wp (the distance between zeros on either side of the main
peak).Note what happens to each width as n oo.
->\342\200\242
Using wx and wp as estimates
of A.t and A/?,
check that the uncertainty principle is satisfied. Warning: If you try
calculating ap, you're in for a rude suiprise. Can you diagnosethe problem?
(d) Use <\302\243(/?, 0) to calculate (/?), {p2), and ap (at time t = 0).
(e) ChecktheHeisenberg uncertainty principle for this state.
[3-%i
i<\">=2<r>-(*9'
2(T)
= [3.97]
Lx^-\\.
This is calledthe virial theorem. Use it to prove that (T) = (V) for stationary
statesof theharmonicoscillator, and check that this is consistent with the results
you got in Problems2.11 2.12. and
(lA/2)Wri(jr) + fo(*)].
**Problem 3.33 Findthe elements matrix {n\\x\\n') and {n\\p\\n') in the (orthonormal)
basis of stationarystates the harmonic for oscillator (Equation 2.67). You already
elements
calculated the \"diagonal\" =
(\302\273 ri) in Problem 2.12; use the same
technique for the general case. the corresponding (infinite) matrices, and
Construct X P.
(\302\2731*1\302\273')
= + \342\200\242 [3-98]
yj^ (Vn7*,,.,,'-! V\302\273~5\342\200\236',\342\200\236-i)
largest possible value of (p) in such a state? If it assumes maximal this value at
time t = 0, what is W(x, r)?
a-\\a) = a\\a)
(the eigenvaluea canbe any complex number).
CO
|Qf) = ]Pc\342\200\236|/?).
//=0
\\n)^e-iE\"'/h\\n),
and show that \\a(t)) remains an eigenstate of a-, but the eigenvalue evolves
in time:
a(t) = e~itota.
(f) Is the groundstate(|\302\253= 0)) itself a coherent state? If so, what is the
eigenvalue?
aAaB '
^^C)
At At
^I>~((02+(\302\243>)2), [3-99]
yy At At Al
what is |4(/))?
what is \\S{t))l
For interesting commentary and references, see R. R. Puri, Phys. Rev. A 49, 2178 (1994).
Further Problems for Chapter 3 129
the matrix
/1 0 0N
H = Ao> 0 2 0
\\0 0 2,
/0 1 0\\ /2 0 0N
A = X 1 0
0]
, B = /x 0 0 1
\\0 0 2/ V0 1 0,
of H,
(a) Find the eigenvalues and (normalized)eigenvectors A, and B.
values might you get, and what is the probability of each?Answer the same
* ^Problem 3.39
fix+x0)
= ^^/11 f(x)
(where xq is any
constant distance). For this reason, p/h is called the
generator of translations in space. Note: The exponential of an operator is defined
by
the power series expansion: e& = 1 + Q + (1/2)(22+ d/3!)fi3+ ....
(b) If vj/(.v, r) satisfies the (time-dependent) Schrodinger equation, show that
(c) Show that the expectation value of a dynamical variable Q(x, p, t), at time
Use this to recover Equation 3.71. Hint: Let to = dt, and expandto first
order in dt.
* ^Problem 3.40
for a free particle, and solve it. Answer: exp(\342\200\224ip2t/2mh) 0(p, 0).
independent of time.
(c) Calculate (p) and {p2) evaluating the appropriate by integrals involving 0,
and compare your answers to Problem2.43.
(d) Show that {H) = (p)2/2m + {H)q (where the subscript 0 denotes the
stationary gaussian), and comment on this result.
= <*(.*.
{Q(t))= (*(A-.r>|j2l*(.v.0> 0)11/-1 QUmx.0)).
where U = exp{\342\200\224iHr/ti). This says that you can calculate expectation values of Q either by
sandwiching Q between ^(.v, t)* and *(.v.f), as we have always done (letting the wave functions carry
the time dependence), or else by sandwiching U~]QU between *(.v.0)* and *(.v.0). letting the
operator carry the lime dependence. The former is called the Schrodinger picture, and the lallcr ihe
Heisenberg picture.
CHAPTER 4
QUANTUM MECHANICS IN
THREE DIMENSIONS
4.1 SCHRODINGER
EQUATION IN SPHERICAL COORDINATES
says
=
ili\342\200\224 HV; [4.1]
at
{-mv2 + V =
\302\261-{p2.+p1x + p1) + V
2 2m
h B fl B Ad
Px -^77-- Pv -^77-. Pc \"* 7 7-. [4.2]
i 3.v i 3v / az
'Where confusion might otherwise occur I have been putting \"hats\" on operators, to distinguish
them from the corresponding classical obscrvables. I don't think there will be much occasion for
ambiguily in this chapter, and the hats gel lo be cumbersome,so I am going lo leave them off from
now on.
131
132 Chapter 4 Quantum Mechanicsin Three Dimensions
or
[4.3]
9* n2 ,
=
ih\342\200\224 [4.4]
dt 2m
where
a2 a2 92
V2 = [4.5]
9^+9v^+9?
is the
Laplacian, in cartesian coordinates.
The potential energy V and the wave function ^ are now functions of r=
(.v. v, z) and t. The probability of finding the particlein the infinitesimal
volume
|2 j3
\\V\\-d*r= 1 [4.6]
/
with the integral taken is
over all space.If the potential independent
of time, there
^\342\200\236(r.t)
=
^,,(^-1 E\"'/fi, [4.7]
h1 ,
\342\200\224
V^ + Vf = Eyjt. [4.8]
2m
^{T,t) = [4.9]
J2c\"^n(T)e-iE\"'/r\\
*Problem4.1
(a) Work out all of the canonical commutation relations for componentsofthe
r
operators and p: [a-, v], [x, py], [x, px], [pY. pz], and so on. Answer:
in, Pj]
= -lPh
rj] = ihSij, [/7./7] = [ph pj] = 0, [4.10]
component.) Hint: First check that Equation 3.71 is valid in three dimensions.
oxaPx
> h/2, a,alh. > h/2. aza,h> h/2, [4.12]
, 1 a /?3\\ 1 3 / 9 \\ 1 ( 92 \\
9r
/\342\226\240-
\\ 9/-/ r*sm9d9\\ 99 J r2sin-9 \\30^/
then,
In sphericalcoordinates, the
time-independent Schrodinger equation reads
r
h2 1 3 / 23^\\
\"
1 3 / . 1
/32^
+
<hfr\\
+
2m 7-lu- V ~dr~) ~r~l^9d9 \\*m0~dO ) ~rh*n1~9\\WJ]
+ V$ = E\\lf. [4.14]
-In principle, this can be obtained by change of variables from the cartesian expression
(Equation 4.5). However, there-arc much more efficient ways of getting it; see. for instance, M. Boas,
Mathematical Methods in the Physical Sciences. 2nd ed.. (Wiley. New York. 1983), Chapter 10.
Section 9.
134 Chapter4 Quantum Mechanics in Three Dimensions
FIGURE 4.1: Spherical coordinates: radius r, polarangle 9, and azimuthal angle \302\242.
//2 Y d R 9 R 92y
sin0\342\200\224 +
2/77 P-dr [! dr ) r2 sin 0 90 M J ,-2 sin2 0 902 J
+VRY=ERY.
-> it-T.
Dividing by RY and multiplying by -2////-//7
[1 d ( -ydR\\ 2m/-2 )
+ ^1^(,^ + -^1=0.
Y sin 0 90 39 J sin2 0 902
The term in the first curly bracket depends only on /\\ whereas the remainder
1 d / -,dR\\ 2m r2
[4.16]
1 f 1 9 / . JY\\ 1 92y)
\342\200\224
sir\302\2730 ) + -^-- =-/(/ + 1). [4.17]
Y sin 0 90 90/ sin20 90:
(c) What is the degeneracy of \302\24314, and why is this case interesting?
it becomes:
9 ( dY\\ d2Y 9
[418]
sin^r^)+^=-'('+i)sin^
You
might recognize this equation\342\200\224it
occurs in the solution to Laplace's equation
in classical electrodynamics. As always, we try separation of variables:
f 1 o ) 1 dr$
sin0\342\200\224 sin0\342\200\224 + /(/+1)51112 0+--^=0.
0 dO \\ d9 J
<$> d<$>-
The first term is a function only of 9, and the secondis a function only of 0, so
each must be a constant.This time4 I'll call the separation constant m2:
'
1
+ 1) sin29 = m2:
\342\200\242 \342\226\240
0
sin
\302\253d
9\342\200\224
d9
(
\\
J\302\256'
sin0\342\200\224\342\200\224
d9
+ /(/ [4.20]
,\302\243/2*=-m'. [4.21]
O dfr
\"'Again, there is no loss of generality here, sinceat this stage in could be any complex number:
in a moment, though, we discover
will that m must in fact be an integer. Beware: The letter 111 is
now doing double duty, as mass and as a separation constant. There is no graceful way to avoid this
confusion, since both uses are standard. Some authors now switch to M or /i for mass, but I hate to
change notation in mid-stream, and I don*t think confusion will arise, as long as you are aware of the
problem.
136 Chapter 4 Quantum Mechanicsin ThreeDimensions
is
The 0 equation easy:
*
d2\302\256
^,,,
= -m-$
...2*
=*.
. _Jm*
= e,m<p.
\302\242(0) [4.22]
[Actually, there are /wo solutions: exp(zw0) and exp(\342\200\224?'m0), but we'll cover the
latter by allowing m to run negative. There could also be a constant factor in
front, but we might as well absorb that into 0. Incidentally,in electrodynamics we
would write the azimuthal function (\302\242) in terms of sines and cosines, instead of
exponentials,becauseelectricpotentials must be real. In quantum mechanics there
is no such constraint,and the exponentials are a lot easier to work with.] Now,
when 0 advances by 2rt, we return to the same point in space (see Figure4.1),so
it is natural to require that5
The 9 equation,
+ [1(1+ 1) sin2 9
-
m2]0 = 0. [4.25]
sin9^ (sin 0^)
0(0) = A
Pj\" (cos 9), [4.26]
I\"1'
/ d \\
= (1 - a-2)I'\"I/2 \342\200\224
P,(x). [4.27]
P/\"(jr) f
J
[4-28]
M-Mil1*2-1*-
is more
\342\226\240\"'This slippery than it looks. After all. the probability density (|4>|2) is single-valued
regardless of m. In Section 4.3 we'll obtain the condition on in by an entirely different\342\200\224and more
compelling\342\200\224argument.
^Notice that
Pj~'\"
\342\200\224
Pj\".
Some authors adopt a different sign convention for negative values
of in: see Boas (footnote 2), p. 505.
Section 4.1: Schrodinger Equation in
Spherical Coordinates 137
For example,
1 ( d \\ ~, , 1 -,
and so on. The first few Legendre polynomials are listed in Table 4.1. As the name
a factor of Vl \342\200\224
x1:
9 r-
2
^(.v) = ]f1
-(3a* 1) = 3(1-a-2).
(l-.v2)(\302\243) 2
/ = 0, 1. 2. m = -l. -/+1 0.
\342\226\2401. 1, ... , /- 1. /. [4.29]
But wait! Equation 4.25 is a second-order differential equation: It shouldhave two
linearly independent solutions, for any old values of / and m. Where are all the
138 Chapter 4 Quantum Mechanics in ThreeDimensions
r0
-
p\302\260- i1 pV =
^(3 cos29-1) p,'(^
-e&
P\\
= sin 9 p!
= 15 sin 9(1- cos2 9)
'i.:
= cos 9 p}
= 15 sin2 $ cos0 ^(^
P\\]
8-** -^
Pl = 3s\\rr9 rh \\sin 9(5cos20-1) ->Z
P\\
= 3 sin B cos e P{h | (5 cos3 0 - 3 cos 0) pV(W) CO ^>
(a) (b)
=
dyT r1^medrdOd(f), [4.30]
nlTt nrt
poo
I \\R\\2r2dr=[ and / / |r|2 sin0</0 d<t>
= 1. [4.31]
Jo Jo Jo
'/2
COS0 ^ \" = ~ 3
(5 9
C\302\260S3 ^
C\302\260S
y?-(i)\"2 (ifc)
= +
'/2
sin 0 (5 cos2 0-
y^1 (^) 1)*\302\261'*
= I ~\342\200\224
I sin- 6> cos Qe x-\"''
^-(4)^(3^-1) Yf-
'/2
Yfl
=+ sin & cos Be**
(\302\247|)
=
/(2/+l)(/-|m |)!
V\"\"*
,-mfl4 ,\342\200\236
1^(0, 0) g/A/ +
P/\"(cos0). [4.321
47T (/ |w|)!
where e = (\342\200\2241)'\" for m > 0 and e = 1 for /// < 0. As we shall prove later on,
they are automatically orthogonal, so
pin pir
/
' &]*[YI?(0,0)]
/ [Yi\"{0, sin0</0<ty
=
5//'5mm', [4.33]
JO JO
In Table 4.3 I have listedthe first few spherical harmonics. For historical reasons, /
is calledthe azimuthal quantumnumber, and /// the magnetic quantum number.
^Problem 4.3 Use Equations 4.27, 4.28,and 4.32, to construct Fj and 7-,1.Check
that they are normalizedand orthogonal.
0(0) = Aln[tan(0/2)]
8Thc normalization factor is derived in Problem 4.54: e (which is always 1 or -1) is chosen
for consistency with the notation we will be using in the theory of angular momentum: it is reasonably
i\" lv\302\273'\\*
=
y,-'\" (-\\)\"'{y;\
140 Chapter 4 Quantum Mechanics in Three Dimensions
\302\273//'\342\226\240 [4.34]
\302\243\342\204\242w*=(tt
Notice that the angular wave function, Y(9, \302\242),is the same for all
part of the
spherically symmetric potentials; the actual shape of the potential, V{r), affectsonly
the radial part of the wave function, R(r), which is determined by Equation 4.16:
'
odR'
d_ (
-^[V(r)-E]R
= l(l + l)R. [4.35]
dt n~
and hence
=
h2 1(1 + 1)
Veff V + [4.38]
2/7/ r2
Those w's are masses,of course\342\200\224ihc separation consiant m does not appear in the radial
equation.
Section 4.1: Schrodinger Equation in SphericalCoordinates 141
0, if r < a:
V(r) =
oo, if r > a.
[4.40]
[4.41]
dr2 L
>-2 J
where
s/lmE
k_ [4.42]
d2u
y
= \342\200\224
k^u =$> u(r) = A sin(A-r) + B cos(kr).
dr
But remember, theactual radial wave function is R(r) = w('*)/'\". and [cos(kr)]/r
blows up as r -> 0. So10 must we choose B = 0. The boundary condition then
requires sin(/rfl)
= 0, and hence ka n.
= rm, for someinteger The allowed energies
are evidently
n-n-n-
E\342\200\2360
= (n = 1.2.3, ...), [4.43]
2m a2
I//- at the origin is normalizable (becauseof the r2 in Equation 4.31). For a more compelling proof
that B = 0. sec R. Shankar. Principles of Quantum Mechanics (Plenum. New York. 1980).p. 351.
142 Chapter4 Quantum Mechanics in Three Dimensions
Normalizing u(r) yields A = -J2ja\\ tacking on the angular part (trivial, in this instance,
since Yq{0, \302\242)
=
1/V4~7r), we conclude that
1 sin(fljrr/fl)
= \342\200\242
Vr\342\200\236oo -j= [4.44]
V27Ttf r
[Notice that the stationary states are labeled by three quantum numbers, /7,/, and
tn: tynimO'i #> 0). The energy, however, depends only on and /: n \302\243\342\200\236/.]
where ji(x) is thespherical Bessel function of order /, and /?/(a) is the spherical
Neumannfunction /. of order They are defined as follows:
For example,
\\a-<1y/ a \\xdxj a3
3 sin a \342\200\224
3a cos x \342\200\224
x2 sin x
1-a2/2 + a4/4! ),
1 A A2
7o(.v)%l; ./iCO%-: /2(.v)
* -;
/10(.0\302\276\342\200\224:
a 3 15
etc. Notice that Bessel functions are finite at the origin, but Neumann functions
blowup at the origin. Accordingly, we must have Bi = 0, and hence
TABLE 4.4: The first few spherical Bessel and Neumann functions,
7/z (-v) and /i/(.v); asymptotic forms for small .v.
sinx COS
./0=
.V
=
\302\2730
x
./3= 1 \\ . 3 (3 1 \\
ji
J~ \342\200\224r sin.v ~ cos,v ~
>h=-\\ \342\200\224$ cos .v y sin .v
V.v3 W x2 Vx3 x) A-
-
Jl
2'n
(2/+ 1)!A
I
' ,^.(2^
i
1, for.v\302\253i.
2'/! x/+1
functions are oscillatory (see Figure 4.2); each one has infinite number an of zeros.
V8 / v10 \\ 12 V .14 X
-.31-
But (unfortunately for us) they are not located at nice sensible points (such as n.
or mj, or
something): they have to be computed numerically.11 At any rate, the
k = -/3,,/. [4.49]
a
where /3,,/ is the \302\253thzero of the /th spherical Bessel function. The allowed energies,
then, are given by
=
E,\342\200\236
/^Ai [4.50]
2m ci-
9. =
fnlmir. \302\242) A\342\200\236,M[}\342\200\236ir/a)Y;\"{0. \302\242). [4.51]
Problem 4.7
(b) Expand the sines and cosines to obtain approximate formulas for /ii(.v) and
ni(x), valid when .v <$C 1. Confirm that they blow up at the origin.
Problem4.8
(a) Check that Arj\\(kr) satisfies the radial equation with V(r)
= 0 and /=1.
Abramowiiz and Siegun. eds.. Handbook of Mathematical Functions. (Dover. New York. 1965).
Chapter 10. provides an extensive listing.
Section 4.2: The Hydrogen Atom 145
* *Problem4.9 A
particle of mass m is placed in a.finite spherical well:
if r < a;
V(r) =
-Vb.
\342\226\240
0, ifr>a.
Find the groundstate, by solving the radial equation with / = 0. Show that there
The
hydrogen atom consists of a heavy, essentially motionless proton(we may as
well put it at the origin), of chargee, together with a much lighter electron (charge
that
\342\200\224e) orbits around it, bound by the mutual attraction of oppositecharges(see
Figure 4.3). From Coulomb's law, the potential energy (in SI units)is
e2 1
V(r)
=
[4.52]
solutions this time\342\200\224we'll work them out in detail, by the method we used in the
analyticalsolutiontotheharmonicoscillator. (If any step in this process is unclear,
you may wish to refer backto Section2.3.2 for a more complete explanation.)
-e
(electron)
sJ-lmE
K = [4.54]
h
me-
p= Kr, and po \342\200\224
[4.55]
27T\342\202\254qITK
so that
d2u po /(/+1)
1 T t [4.56]
dp*
d2u =
u.
dp7
up (as p
\342\200\224
but ep blows \342\200\224>
oo), so B 0. Evidently,
~ Ae~p. [4.58]
u(p)
for largep. On the other hand, as p -\302\2730 the centrifugal term dominates;12
approximately,
then:
d2u /(/+1)
dp- p-
'-This argument does nol apply when I = 0 (although the conclusion. Equation 4.59, is in fact
valid for that case loo). But never mind: All I am Irving lo do is provide some motivation lor a change
of variables (Equation 4.60).
Section 4.2: The Hydrogen Atom 147
=
u(p) Cpl+* +Dp~!.
~ /+1
\302\253(p) Cp [4.59]
for small p.
The next
step is to peel off the asymptotic behavior, introducing the new
function v(p):
u(p)
= pl+le-pv(p), [4.60]
in the hope v(p) that will turn out to be simpler than u(p). The first indications
are not auspicious:
du dv~
l -pp
\342\200\224
=
p'e (/ + 1- p)v + p \342\200\224
dp dp.
and
d~v
- dv - 2(/ +
+ 2(/ + 1 p)\342\200\224 + [p0 1 )]v = 0. [4.61]
p\342\200\2247
dp- ap
v(p)
= \342\226\240 [4.62]
Ylc.iPi
7=0
by term:
oc 00
dv x\342\200\224v .. , -c-^ .
dp
u u
d~V t\342\200\224v ; I
= + ikv+ip'-1.
^ E^./
Insertingthese into
Equation 4.61, we have
OO oc
Y,j(j
+ DoW+2(/ +1) j^u + u<--j+\\pj
./=0 j=0
i=0 7=0
or:
= [ 2(,/+/+1)-/00 ]
W
{(; +1)0- +*+ 2)}^
[4.63]
This recursion formula determines the coefficients, and hence the functionu(p):
We start with co (this becomes an overall constant, to be fixed eventually by
normalization), and Equation 4.63 gives us c\\\\ putting this back in, we obtain ci,
and so on.13
Now let's see what the coefficients look like for large j (this corresponds to
large p, where the higher powers dominate). In this regimethe recursion formula
says14
2' -~ 2
Cj+[ ~
~
JuTT)Cj JTiCj-
'-'You might wonder why ] didn't use the series method directly on u(p)\342\200\224why factor out the
asymptotic behavior before applying this procedure? Well, the reason for peeling off p'+l is largely
aeslhelic: Wilhoul ihis. the sequence would begin wilh a long siring of zeros(the firsl nonzero coefficient
being r/+|): by factoring out p/+l we obtain a scries thai starts out wilh The
p\302\260. e~p factor is more
critical\342\200\224if you don't pull lhal out. you gel a three-term recursion formula, involving fy+2* <*/+i>
and
Cj (try ii!) and lhal is enormously more difficult to work wilh.
14You might ask why 1 don't drop ihe 1 in j + 1\342\200\224-afterall, I am ignoring 2(/+ 1)- po in ihe
numerator, and 2/ +2 in the denominator. In ihis approximation il would be fine to drop ihe 1 as well,
but keeping il makes the argument a lillle cleaner.Try doing il without the 1. and you'll see what I
mean.
Section 4.2: The Hydrogen Atom 149
V
Cj = \342\200\224cq. [4.64]
so
00
2J
=
V(p) =C0y \342\200\224pi CQ(TP,
7
7=0
and hence
u(p) = cQpl+[ep, [4.65]
equation\342\200\224they just don't happen to be the ones we're interested in,because they
aren't normalizable.) There is only one way out of this dilemma: The series must
terminate. Theremust occursome maximal
integer, ./max, such that
cOm,v+D=0. [4.66]
20\"max+/+0-/00 = 0.
Defining
n = ymax +/ + 1 [4.67]
P0
= In. [4.68]
?
n ^
k me4
E = -^- = ,,,-,- [4.69
2m
S7t2\342\202\254^h2pl
This is the famous Bohr formula\342\200\224by any measure the most important result in
all of quantum mechanics. Bohr obtained it in 1913 by
a serendipitous mixture
of inapplicable physics and premature
classical quantum theory (the Sehrodinger
equation did not come until 1924).
CombiningEquations4.55 and 4.68, we find that
me' 1 1
K=( [4.71]
\\47teo^'2 / n an'
where
a =
47re0^ = 0.529x 10\"10 m [4.72]
me\"
P
= [4.73]
an
The spatial wave functionsfor hydrogen are labeled by three quantum numbers (n,
I, and m):
fnlm(r,0, \302\242)
= Rnl(r) lf(0, 0), [4.74]
2(j+l + l-n)
Cj+1
(./ + 1)04-2/ +
,7' [4.76]
2)6
2~l
m
E\\ =- = -13.6 eV. [4.77]
2ft2 \\47T\342\202\2540i
is traditional
\342\226\2405It to write the Bohr radius with a subscript: oq. But this: is cumbersome and
unnecessary, so 1 prefer to leave the subscript off.
Section 4.2: The Hydrogen Atom 151
[4.79]
a
it,
Normalizing in accordance with Equation 4.31:
i2
f
/
id i2
\\Rw\\~rdr
2 j lC0l
= -^=- If e \342\200\224Irla
~r/a
r~dr
2 7
= 11^
c0
- - 11,
k a2 Jo 4
so co \342\200\224
%1-sfa. Meanwhile,
=
Y\302\256
ljy/4n, and hence the ground state of hydrogen
is
1
= -r/a
^1()0 (r, 0, \302\242) [4.80]
VjTtf,3
If n = 2 the energyis
=
-13.6 eV
\342\200\224\342\200\224 =
E2 -3.4 eV; [4.81]
this is the first excitedstate\342\200\224or rather, states, since we can have either I \342\200\224
0 (in
ci \342\200\224
(using
\342\200\224Co j = 0),.. and =
\302\2422 0 (using ;' = 1).,
so v(p) \342\200\224
cq(1
\342\200\224
p), and therefore
[Notice that the expansion coefficients {cj} are completelydifferent for different
[4.83]
4a2
Problem 4.11.)
152 Chapter 4 Quantum Mechanics in ThreeDimensions
For arbitrary n, the possible, values of I (consistent with Equation 4.67) are
/= 0,1,2,...,^-1, [4.84]
/j-i
d(n)
=
J] (21 + 1) = a2; [4.85]
The
polynomial v.(p) (defined by the recursion formula, Equation 4.76)isa function
well known to applied mathematicians; apart from normalization, it can be written
as
- [4-86]
v(p) ^SLi(2a)>
where
=
*\342\226\240*
IWi^\302\247^'m O' [*->(2rH *\342\204\242 \302\253\342\226\240 [4.89]
16As usual, there, are rival normalization conventions in the literature; 1 have adopted the most
nearly standard one.
17If you want to see how the normalization factor is calculated, study (for example), L. Schiff,
Quantum Mechanics, 2nd ed., (McGraw-Hill, New York, 1968), page 93.
Section 4.2: The Hydrogen Atom 153
Lp-= 1
Lr-= -x- +- 1
= --x2- - Ax + 2
L2
L6 = --X\302\253--36^ + 450x4
- 2400x3+ 5400.\302\2762 -4320x + 720
polynomials,
Lp (x).
j 0 _ 11 L2 = 2
L,0-
Lf=-jc+l
L2 = -6x+18
=
\302\243.0 .^-4^+2 L|=12x2-96x+144
L]-l L3 = 6
Z.J
= -2x + 4 Lj= -24x + 96
- 18
1^ = 3x2 18x + L|=60x2-600x+1200
case of the spherical well, you may recall, the energies dependalsoonI
(Equation 4.50). The wave functions are mutually orthogonal:
/ ^nim ifti'i'm*
>'2 sin # dr d$ # = 8ntI>8n>&mri
\342\226\240 [4,90]
and (for n, ^ nr) from the fact that they are eigenfunctionsof H with distinct
eigenvalues.
^Problem 4.10 Work out the radialwave functions ./\302\275> ^31, and R32, using the
recursion formula (Equation 4.76). Don'tbotherto normalize them.
154 Chapter 4 Quantum Mechanics in Three Dimensions
TABLE 4.7: The first few radial wave functions for hydrogen,
RiQ=2a~3/2 exp{-r/d)
\342\200\224
#2i =-f=a-312 exp (-r/2a)
V24 a
n _ 2 a,,-3/2
/tan
^30
\342\200\224\342\226\240
-r\342\200\224\342\226\240
V27 3 a 27
i'-h^i&h^-^ \\a
^4^0474(^1\302\276\302\276^^
- *
-
^4]
ltiV3 -3/204f4&)2)i-p^)
^--ih^i^uw^^
W= \342\200\224l\342\200\224crm(\342\200\224) exp (-r/4a)
768^35
^Problem 4.11
^Problem 4.12
(a) Using Equation 4.88, work out the first four Laguerre polynomials.
(b) Using Equations 4.86, 4.87, and 4.88, find v(p), for the case n = 5, I = 2,
FIGURE4.4: Graphs
of the first few hydrogen radial wave functions, HM/(r).
*Problem4.13
hydrogen? (The answer is not zero!) Hint: First you must figure out theprobability that
Problem4.15 A
hydrogen atom starts out in the following linear combinationof
the stationary states n \342\200\224
2,1
\342\200\224
1, m..= 1 and n \342\200\224 /
2,
= 1, m \342\200\224
\342\200\2241:
1
W(r,0)
= \342\200\224Ofon+^21-1).
:x/2
FIGURE 4.6: Surfaces of constant|^|2 for the first few hydrogen wave functions.
Reprinted by permissipn from Siegmund Brandt and Hans Dieter Dahmen, ThePicture
Bookof Quantum Mechanics, 3rd ed.. Springer, New York (2001).
In
principle,
if you put a hydrogen atom into some stationary state ^nim, it should
stay there forever. However, if you tickle it slightly collision with another
(by atom,
say, or by shining light on it), the electron may undergo a transition to some other
- ~
Ey = Ei -Ef = -13.6 eV
[ \\
n.
)
. [4.91]
n;f
EY
= hv. [4.92]
[4.93]
where 2
R= -^- (
-^
] = 1.097 x 107m\"1 [4.94]
is the
known
Rydberg as constant. Equation 4.93 is the Rydberg formula forthe
spectrum hydrogen; of it was discovered empirically in the nineteenth century, and
the greatest triumph of Bohr's theory was its ability to account for this result\342\200\224and
region; they
constitute the Balmer series. Transitions to rif = 3 (the Paschen
series)are in the
infrared; atid so on (see. Figure 4.7), (At room temperature,most
hydrogen atoms are in the obtain the emissionspectrum
ground state; to you
must
first populate the various excited states; typically this is done passing an by
electric
^Problem 4.16 A
hydrogenic atom consists of a single electron orbiting a nucleus
with Z protons (Z = 1 would be hydrogenitself,Z \342\200\224
2 is ionized helium, Z = 3..
By its nature, this involves a unie-ifependent interaction, and the details will have to wait for
Chapter 9; for our present purposes the actual mechanism involved is immaterial.
19The photon is- a quantum of electromagnetic radiation; it's a relativistic objectif thereever was
one, and therefore outside the scope of rionrelalivistic quantum mechanics. It will be useful in a lew
places to speak of photons, and to invoke the Planck formula for their energy, but please bear in mind
that this is external to the theory we are developing..
Section 4.2: The HydrogenAtom 159
n
\342\226\240
OQ
0
6
5
-1.0 4
Y I Y 'v
3
-2.0 Paschen
series
-3.0 u u u u u
-4.0 Balmer
series
-5.0
-6.0
>
-7.0h
S3
CD
-8.0
C
LU
-9.0h
-10.0
-11.0
-12.0
-13.0
v u u u ,r u
1
\342\200\242\342\200\224
\342\200\22414.0 Lyman series
is doubly ionized lithium, and so on). Determine the Bohr energies En(Z), the
binding energy E\\(Z), the Bohr the Rydberg constantR(Z)
radius a(Z), and for a
values.) Where in the electromagnetic spectrum would the series fall, Lyman
for
(b) What is the \"Bohr radius,\" ag, for this system?Work out the actual number.
much energy (in Joules) would be released? What would the wavelength of
the emitted photon (or, more likely, graviton)
be? (Express your answer in
light years\342\200\224is the remarkable answer20 a coincidence?)
As we have seen, the stationary states of the hydrogen atom are labeledby three
quantum numbers: n, I, and m. The principal quantum number determines
(\302\253.)
the
Lx = =
ypz-zpy, Ly zpx\342\200\224 xpz, Lz=xj?y-ypx. [4.96]
Px \342\200\224>
\342\200\224ihd/dx, py
\342\200\224>
\342\200\224ihd/Byi pt
-> \342\200\224ihd/dz- In the following section we'll
obtain the eigenvalues of the angularmomentum operators by a purely algebraic
technique reminiscent of the one we used in Chapter get the allowed energies 2 to
of the harmonicoscillator; all it is based on the clever exploitation of commutation
relations. we will turn
After that to the more difficult problem of determining the
eigerifunctions.
4.3.1 Eigenvalues
[Lx, Ly]
=
[ypz
-
zpy, zpx
- xpj
= [yPz, ZPx]
~
[yPz> *Pz\\
-
kpy, zpx\\ +[ZPy, *Pzl [4.-97]
From the eanonieal commutation relations (Equation 4.10) we know that the only
operators here thatfail to commute are x with px, y with py, and z with
pz. So
[LXi Ly]
= ypxipz, z] + xpylzi pz] = ih{xpy- ypx)
=
ihL-z. [4.98]
Of course, we could have startedoutwith [Ly, Lz] or [Lz, Lx], but there is no need
to calculate these separately\342\200\224-we can get them immediately by cyclic permutation
of the indices(x -> v, y
-> z, z-> x):
relations
These are the fundamentalcommutation for angular momentum;
~ '
vLsOLy >~\\(LZ)\\. [4.100]
[L2, LA
=
[L2X, Lx] + [L2, Lx] + [L% Lx]
= Ly[Ly,Lx]+ [Ly, Lx}Ly+Lz[Lz, LJ+ [LZJ LX]LZ
= Ly(-ihLz) + {-ihLz)Ly + Lz{ihLy) + {ihLy)Lz
= .0.
(I used Equation 3.64 to the commutators; simplify note also that any operator
commutes with itself) It follows,of course, that L2 also commutes with Ly and
Lz-
[L2,Z^] = 0, [L2,Ly]=0, [L^LJ^O, [4.102]
162 Chapter 4 Quantum Mechanics in Three Dimensions
or, morecompactly,
[L2,L]= 0. [4.103]
L+ = Lr \302\261iL.y [4.105]
SO
= + .
[LZ,L\302\261]
\302\261HL [4.1Q6]
And, of course,
= o.
[LM\302\261] [4.107]
I claim that if f is an
eigenfunction of L2 and Lz, so also is L\302\261/': Equation 4.107
says
L2(L\302\261/)
=
L\302\261(L2/)
= L\302\261(Xf)
=
HL\302\261f), [4.108]
says:
= + = \302\261hL + L^{iif)
Lz{L\302\261f) {LzL\302\261-L\302\261Lz)f L\302\261LJ + f
= ({j,\302\261K](L+f% [4.109]
so L + f is an
eigenfunction of Lz with the new eigenvalue fi \342\226\240+
/z. We call L+
the \"raising\" of
operator, because it increasesthe eigenvalue Lz by h, and L_ the
\"lowering*' operator, because it lowersthe eigenvalue h. by
Figure 4.8). To ascend the ladder we apply the raisingoperator, and to descend, the
lowering operator. But this process cannotgo on forever: Eventually we're going
Section 4.3: Angular Momentum 163
L+ft = 0. [4.110]
Let hi be the eigenvalue of Lz at this top rung (the appropriatenessof the letter
\"/'\" will appear in a moment):
Now,
L+LT = (Lx
+
iLy)(Lx ^iLy)
= L2 + L2^ i(LxLy - LyLx)
I?
\342\200\224 -
l}z^i(ihLz),
L2 = + hLz. [4.112]
L\302\261LT L2Z^
It follows that
and henee
L-fh = 0. [4.114]
Let hi be the eigenvalue of Lz at this bottomrung:
Lzfb
=
Hlfbi L2fh = kfh. [4.115]
I>fp
\342\200\242=
(L+L- + L2- hLz)fb = (0 + h2t - h2l)fb
=
h2J(l -\\)h,
and therefore
=
X /z27(7-l). [4.116]
= HH(l + rin
iS-f? 1)jf; LJ\342\204\242
=
hmfi [4.118]
where
- 1, /,
1 = 0, 1/2, 1,-3/2, . m = \342\200\224
I, +
\342\200\2241 1, I [4.119]
Fora given value of /, there are 2/ -f 1 different values of m (i.e., 11 + 1 \"rungs\"
on the \"ladder\.
Some peoplelike to illustrate this result with the diagram in Figure 4.9
(drawn for the arrows are supposedto represent
case I = 2). The
possible angular
Vfr and
\342\200\224
2.45), their z components are the allowed values of 1, m
\342\200\224
(\342\200\2242, 0, 1,
2). Notice that the magnitude of the vectors(theraditisof the sphere) is greater
uncertainty principle (Equation 4.100)says that's impossible. \"Well, all right, but
defined value, then Lx and Ly do not. It is misleading even to draw the vectors
in Figure4.9\342\200\224-at best they should be smeared out around the latitude lines, to
indicate that Lx and Ly are indeterminate.
I hope you're impressed: By purely algebraic means, starting with the
fundamental commutation relations for angular momentum (Equation 4.99), we have
determinedthe eigenvalues L2 and of
seeing the eigenfunctions Lz\342\200\224without ever
I should warn you that this is a much messierbusiness. lust you know where so
L2 and Lz are nothing but the old spherical harmonics, came upon by which we
a quite different route in Section 4.1.2 (that's chose the letters I and m, why
I
of course). I can now tell you why the spherical harmonicsare orthogonal:
And
They
are
eigenfunctions of hermitian operators (L2 and Lz) belonging to distinct
eigenvalues (Theorem 2, Section 3.3.1),
*Problem 4.1\302\276 The raising and lowering operators change the value of m by one
unit:
L\302\261f{n={Af)f^\302\261\\ [4.120]
Lx and
Ly are observahles, you may assume they
are herimtian . but
\342\226\240..
prove it
if you like); then use Equation 4.112.Answer:
Af = hy/lQ
+ l)-m(m\302\261\\) =
h^(l Tm)(I \302\261m + 1). [4,121]
Note what happens at the top and bottomof the ladder (i.e., when you apply L_|_
to fj or L_ to /\"').
4.19
^Problem
(a) Starting with the canonical commutation relations for position and momentum
(Equation4.10), work out the following commutators:
(c) Evaluate the commutators[Lz,r2] and {LA, p2] (where, of course, r2 \342\200\224
x2 + +
y\302\260- z2 and p2 = p2x+ p$ + /\302\276).
(d) Show that the Hamiltonian H = (p2/2m) + V commutes with all three
components of L, provided that V\" depends only on r. (Thus H, L2, and Lz
* ^Problem 4.20
(a) Prove that for a particlein a potential V(r) the rate of change of the
expectation value of the orbital angular momentum L is equal to the expectation
value of the torque:
= (N),
|(L)
where
N = rx(-VV).
4.3.2 Eigenfunctions
V = +0
r\342\200\224 + 0\342\200\224 ; [4.123]
dr r dO r sin0 90
meanwhile, r = rr, so
h ~ 9 \342\226\240 * \342\226\240
\342\200\236 -9 ~ * 1 9
L = r(r x r) f- y
(r x 61)\342\200\224+ (r
K
x 0)
v
_
JBr ;W sine 90.
\"^GeorgeArfken and
Hans-Jurgen Weber, Mathematical Methods for Physicists,5th ed.,
Academic Press, Orlando (2000), Section 2,5.
168 Chapter 4 Quantum Mechanicsin Three Dimensions
\302\273
/ --a -: i a
L=-
i 0--0- [4.124]
'
\\
3<9 sin 0 30 /
The unit vectors 9 and 0 canbe resolved into their cartesian components:
\342\200\224
0 \342\200\224(sin0)f-f- (cos0)J, [4.126]
Thus
. 3
(.-sin0i +cos0;) \342\200\224
1 3
-(eos# cos 0i + cos9 sin0 J \342\200\224
sin 9k)
\342\200\224
sin# 30 _
Evidently
/
-\302\276 3 3 \\
= - - -
cos0eot#
\342\200\224-
Lx I sin0-\342\200\224 , [4.127]
z
\\ 69 \302\242\302\242/
h 3
\342\200\224
Ly
\342\200\224 I +CDS0- sin0:eot0 [4.128]
z' V 3n 30,
and
[4-129]
h 3 , 3
L+ \342\200\224
Lx + iLi, = (\342\200\224sin0
i eos0)\342\200\224-
\302\261
\342\200\224
(cos0 i
\302\261 sm0)cot#\342\200\224
i d$ 30 _
But eos 0 i
\302\261 sin0
= e*+ i.
'*,, so
' A
3 3
L+ =
,
+
\342\200\224: \342\200\224
\302\261he\302\261ix> J: COt\302\243 [4.130]
.36\302\273 30 ,
/ 32 32
X+L_ - -^
3 , 3
1 3/.9
\342\200\224 \342\200\224 _L 3*_'
L2 = -h' sin (9 + [4.132]
sin e de \\ m sin2 0 d(j>2
1 3/..3sin.0
2 i-KfJ.
= -A\" \342\200\224
+ _1 3^ = fc2/(/ + l)//!
Z//f jf
sin 9(9
6\302\273 V 96 sin2 302
6\302\273 _
But this is preciselythe \"angular equation\" (Equation 4.18). And it's also an
eigenfunction of Lz, with the eigenvalue mh:
=
*<zf?
j^fr=*mfr,
1 -2 9
+V -f Vifr
= Exfr.
2mr2 dr \\j/
integer values (Equation 4.29). You might suppose that the half-integersolutions
are spurious, but it turns out that they are of profound importance, as we shall see
^Problem4.21
(a) Derive Equation 4.131 from Equation 4A30,Hint: Use a test function;
otherwise you're likely to drop some terms.
^Problem 4.22
(b) Use the result of (a), togetherwith Equation 4.130 and the fact that LzYf
=
Problem 4.24 Two particles of mass m are attachedto theends of a massless rigid
rod of length a. The systemis freeto rotatein three dimensions about the center
(but the center point itself is fixed).
(a) Show that the allowed energies of this rigid rotor are
= = 1, 2,.,.
En \342\200\224-\342\200\224Y\342\200\224,for
\302\253 0,
(b) What are the normalized eigenfunctions for this system? What is the
4*4 SPIN
In classical mechanics, a
rigid object admits two kinds of angular momentum:
orbital (L = r x p), associated with the motion of the center of mass, and spin
(S = 1(d),. with associated motion about the center of mass. For example,theearth
has orbital angular momentum attributable to its annual revolution around the sun,
and
spin angular momentum coming from its daily rotation about the north-south
axis. In the classical context this distinction is largely a matter of convenience, for
when you come right down to it, S is nothingbut the sumtotal the \"orbitaT' of
angular momenta of all the rocks and dirt clods go to make up the earth, asthat
therefore, we use the same word). It doesn't pay to press this too far: The analogy
electron (as far as we know) is a structureless point and its spin angular particle,
It follows(as before)
that the eigenvectors of 52 and satisfy27
51\302\276
S2\\sm)
= h2s(s + l)\\sm); Sz\\sm} = hm\\s m); [4.135]
a contrary
\342\226\240For interpretation, see Hans C. Ohanian, \"What is Spin?\", Am. J. Phys. 54, 500
(1986).
26
We as postulates for the theory of spin; the analogous formulas for orbital
shall take these
angularmomentum were derived from the known
(Equation 4.99) form of the operators (Equation 4.96). In
a more sophisticatedtreatment they can both be obtained from rotational invariance in three dimensions
(see, for example, LeslieE. Ballentihe, Quantum Mechanics: A Modern Development, World Scientific,
Singapore (1998), Section 3.3). Indeed, thesefundamental commutation relations apply \\a all forms of
angular mpmentum, whether spin, orbital, or the combined angular momentum of a compositesystem.
Which Could include some spin and some orbital.
Because the eigenstates of spin are use the \"ket\"
uqL notation
functwnstfor them. (I I will
could have done the same in Section 4.3, of Yfl, but in that context the function
writing \\lm) in place
notation seems more natural.) By the way, I'm running out of letters, so I'll use m for the eigenvalue
of SV, just as I did for Lv: (some authors write M\\ and ms at this stage, just to be absolutely clear).
172 in
Chapter 4 Quantum Mechanics Three Dimensions
and *
= + V-miml 1) (m \302\261
1)), [4.136]
S\302\261\\sm) \\s
hy/sis
1 3
= 05
.*\342\226\240 -, lr -,-..; m = -s, -5+1, ... , s - 1, s. [4.137]
0; electrons have spih 1/2;photons spin 1; deltas.have spin 3/2; gravitons have
have
spin 2; and so oh. contrast, the orbitalBy angular momentum quantum number
I (for electron an in a hydrogen atom, say) Can take on (integer) value you any
please, and will change from one to another the system is perturbed. But s when
is fixed, for given particle, and this makes the theory of spin comparatively
any
simple.28
Problem 4.25 If the electron were a classical solid sphere, with radius
e2
re = \342\200\224
[4.138]
y
4jr6o7wcz
attributable to energy stored in its electric field, via the Einstein formula E = mc2),
and its angular momentum i$ (l/2)h, thenhow fast (in m/s) would a point on the
\"equator\"be moving?Doesthis model make sense? (Actually, the radius of the
electronis known experimentally to be much less than rc, but this only makes
mattersworse.)
28Indeed, in a mathematical sense, spin 1/2 is the simplest possible nontrivial quantum system, for
it admits just two basis states. In place of an infinite-dimensional Hilbcrt space, with all its subtleties
and complications,we find ourselves working in an ordinary 2-dimensional vector space; in place
of unfamiliar differential equations and fancy we are confronted
functions, with 2x2 matrices and
2-component vectors. For this reason, some authors begin-quantum mechanics with the study of spin.
(An outstanding example is John S. Townsend, A Modern Approach to Quantum Mechanics, University
Books, Sausalito, C.A, 2000.) But the price of mathematical simplicity is conceptual absti'action, and I
X+=L [4-140]
X-= V [4-141]
S2 =
-
:: flffl-H).
\342\200\224
= a The second
so c (3/4)/\302\2762 and 0. equation says
' \302\260
'c d\\ /0\\ 3\342\200\2362/0\\
/<T
so d = 0 and / \342\200\224
(3/4)\302\276 . Conclusion:
* = ^14\302\276
l*Co ?)\342\226\240
174 Chapter 4 Quantum Mechanics in ThreeDimensions
Similarly,
h o - --x-
*
PzX+
=
s?x- [4.144]
^x+,
so
'0 1\" '0 0^
S-U=fi S_ = A..
1
[4.146]
,0 0, o,
+ S_) and -
Now S\302\261
=
Sx \302\261
iSy,
so Sx = (1/2)(5+
Sy
=
(1/2/)(5+ 5_), and
hence
^ - 2 or
11 bar\"2U. o
[4.147]
Since Sjt,. Sy.,, and Sz all carry a factor of hj'1,it is tidier to write S = (hj2)a,
where
'Q -f 1 0
'^s(l o)> \302\260>
j 0 7' ^^ 10 -1 [4.148]
These are the famous Pauli spin matrices. Notice that Sx, Sy, S2, and S are all
hermitian (as they should be, since they represent observables). On the otherhand,
S+ and S- are not hermitian\342\200\224evidently they are not observable.
The eigenspinors of Sz are (of course):
H'
X+
=
eigenvalue + X- eigenvalue J. [4.149]
If you measure Sz
on a particle in the general state x (Equation4.139), you could
only possibilities,
\\af + \\b\\2
= l [4.150]
29People often say that |a.|2 is the \"probability that the particle is in the spin-up state,\" but this
is sloppy language; what they mean is that if you measured Sz, la]2 is the probability you'd get h/2.
See footnote 16 in Chapter \"3.
Section 4,4: Spin 175
But what if, instead, you chose to measureSx1What are the possible results,
and what are their respective probabilities?According to the generalized
statistical
interpretation, we heed to know the eigenvalues and eigenspinors of Sx. The
characteristic equation is
-k /J/2 'W ft
= o r
=>\342\226\240 = =*3t=\302\261
ft/2 -X
Not surprisingly, the possible values for Sx are the same as those for Sz. The
= + 'a'
;(?.oca-* tew.'
so ^ +
\342\200\224 a. Evidently the (normalized) eigenspinors of Sx are
1
/J_\\
ft\" 72 ft'
yW
eigenvalue + *\302\253 eigenvalue [4.151]
1
\342\226\240
72'
As the eigenvectors of a hermitian matrix, they span the space; the generic spinor
X (Equation 4.139) can be expressed as a linear combinationof them:
Example
4.2 Suppose a spin-1/2 particle is in the state
Wl+T
=
X
'2
76 V ,
5( h\\ l{ h\\ ft
+ =
6 (+2J 6 N ?
which we could also have obtained more directly:
particle in the state x+- asks, \"What is the z-componentof that particle's
If someone
spin angular momentum?\", answer unambiguously: +H/2. For a
we could
measurement of SV is certain asksinstead,
to return that value. But if our interrogator \"What
you telling me that you don't know the true state of that particle?'' On the contrary;
I know precisely what the state of the particle is: x+. \"Well, then, how come you
can't tell me what the ^-componentof its spin is?\" Because it simply does not
have a particular x -componentof spin.Indeed, it cannot, for if both Sx and Sz
of its spin; let's say he gets the +/I/2. \"Aha!\" (he shouts in triumph), \"You
value
lied! This particlehas a well-defined perfectly value of Sx: h/2,\" Well, sure\342\200\224it
does now, but that doesn't prove it had that value,priorto measurement. your
\"You
\"But
5\302\276.. I was extremely careful not to disturb the particle I measured Sx\"
when
Very well, if you don't believe me, check it out: Sz, and see
Measurewhat you
get. (Of course, he get +/i/2,
may which will be embarfassihg to my if ease\342\200\224but
we repeat this whole scenario over and over, half the time he get will
\342\200\224h/2.)
almost impossible to convey to anyone who has not studied quantum mechanics
(if you don't, you probably haven't understood the problem), cqme back to the
Problem 4.26
(b) Show that the Pauli spin matrices (Equation 4.148) satisfy the productrule
ffj<rk
=
8jk + i
2^ ejkim* [4-!53]
l
where the indices stand for x, y, or z, and is
\342\202\254jm
the Levi-Civita symbol:
+ 1 if jkl = 123, 231,.or 312;-1 if jkl
= 132, 213, or 321; 0 otherwise.
-\"CO-
constantA.
(a) Determinethe normalization
(b) Find the expectation values of Sx., Sy, and Sz,
(c) Findthe\"uncertainties\"
crsr, a$v,
and
0\302\276. (Note:
These sigmas are standard
deviations, not Pauli matrices!)
(d) Confirm that your results are consistent with all three uncertainty principles
(Equation4.100 and its cyclic permutations-\342\200\224only with S in place of L, of
course).
=
(S2).
*Problem 4.29
(c) If you measuredSy,. what values might you get, arid with what probabilities?
Problem 4,31 Construct the spin matrices (S*, Sy, and Sz) for a particle of
5+,. and SL on each of these states. Followthe procedure in used the text for
spin 1/2.
A
spinning charged particle constitutes a magnetic dipole. Its magnetic dipole
moment,ji, is proportional
to its spin angular momentum, S:
\"See. for example,D. Griffiths, Introduction to Electrodynamics, 3rd ed. (Prentice Hall, Upper
SaddleRiver, NJ, 1999), page 252. Classically, the gyromagnetic ratio of an object whose charge--and
mass are identically distributed is q/1m, where q is the charge and m is the mass. For reasonsthai are
fully explained only in relatiyistic quantum theory, the gyromagnetic ratio of the electron is (almost)
exactly twice
\342\200\224
\"the classical value: y \342\200\224
e/m.
Section 4.4: Spin 179
H = -\\lB, [4.157]
B = Bok. [4.159]
The Hamiltonian (Equation 4.158), in matrix form, is
H =-yfcS, [4.1601
=-*\302\247\302\243 _0,).
\302\243
Evidently the energy is lowest when the dipole moment is parallel to the field\342\200\224just
as it would be classically.
Since the Hamiltonian the generalsolution
is time-independent, to the time-
ih%\302\243-=Hx, [4T62]
at
can be expressedin terms of the stationary states:
n
iyBot/2
iyB{)t/2
31
Griffiths (footnote 30), page 281.
32If the particle is allowed to move, there will also be kinetic energy to consider; moreover, it
will be subject to the Lorentz force (q\\ x B), which is not derivable from a potential energy function,
arid hence does,not fit the Schrodinger equation as: we have formulated it so far. I'll show you later
on how to handle this (Problem 4.59), but for the moment let's just assume that the particle is free to
rotate, but otherwise stationary.
180 Chapter 4 Quantum Mechanics in ThreeDimensions
X-(0) = (f) ,
(of course, \\a\\2 + |&|2 = 1). With no essential loss of generality33 I'll write a \342\200\224
:Cos'(a/2) and b = sin (a/2), where a is a fixed angle physical significance whose
,, /cos(ar/2)ez>B\302\260'/2\\ rA..\342\200\236,
To get a feel for what is happening here* let's calculate the expectation value of S,
as a functionof time:
(\302\276)
=
X (OfS^X (?) = (cos(a/2)^*<\"/2 sin^)^^/2)
ft /0 l\\ /cos(a/2)e'>B^/2\\
X [4.164]
2
(l oj ^/2),-^^/2)
\342\200\224
-sina cos[yB[)t).
Similarly,
(Sy)
= X(.()%X(!) = --sinasin()/5oO, [4.165]
and
Evidently (S) is tiltedat a constantangle ql to the z-axis, and precesses about the
field at the Larmor frequency
=
co yB0, [4.167]
theorem (in the form derived in Problem 4.20) guarantees that (S) evolves
according to
the classical laws. But it's nice to see how this works out in a specific
context.
does
\342\226\240^This assume that a and b are refill you Can work out 'the general case if you like, but all
34See, for instance, Tfve Feynman Lectures on Physics (Addison-Wesley, Reading, 1964),Volume
II, Section 34-3; Of course, in the classical case it is the; angular momentum vector itself, not just its
expectation value, that precesses around the magnetic field.
Section 4.4: Spin 181
z\302\273
F = VQi-.B). [4.168]
h(x, y,z)
= -axi + (Bo + az)k, [4.169]
where Bo is a strong uniform field and the constanta describes a small deviation
from homogeneity. (Actually, what we'd like is just the z component^ but
unfortunately
that's impossible\342\200\224it would violate the electromagnetic law V B = 0; like
it or not, an x component comes along for the ride.) The forceon these atoms is
^Griffiths, (footnote 30), page 258. Note that F is the negative gradient, of the energy
(Equation 4.157).
Zn Spin up
Spin down
Magnet
But because of the Larmor precession about Bq, & oscillates rapidly,
and
=
Fz y<xSz, [4.170]
beautifully demonstrating the quantization of angular momentum. (If you use silver
atoms,for example, all the inner electrons: are paired, in such a. way that their
Hamiltohian starts out zero, turns on for a time (as the particle passes through T
0, for t < 0,
\342\226\240H
(*)
=
-y(Bo + aZ)Sz, forO<f <T, [4.171]
0, for t > T.
(I ignore the of B,
pesky x cornporient which\342\200\224-lot reasons indicated above\342\200\224is
irrelevant to the problem,) Suppose the atom has spin 1/2, and startsout in the
state
X ft)
=
ax+ + H-, for t < 0.
(for i > T), The two terms now carry momentum in the z direction (see
Equation 3.32); the spin-up component has momentum
ayTh
Pz = -^-s [4.174]
and it moves in the plus-z direction; the spin-down component has the
opposite momentum, and it moves in the minus-z direction. Thus the beam splits
in two, as before. (Note that Equation 4.174 is consistent with the earlier result
you Want to
prepare a beam of atoms in a given spin configuration, you pass an
B \342\200\224
BQCOs(a)t)k,
sur srn(a)f) .
\\2cd J
(d) What is the minimum field (Bq) required to force a completeflipin Sx?
where the first arrow refers to the electron arid the second to the proton.Question:
What is the total angular momentum of the atom? Let
3^l put them in the ground stale so there won't be any orbital angular momentum to worry about.
39More precisely, each particle is in a linear combination of spin up; and spin down, and the
simply add:
szxiX2 =
(si{) + sf )xiX2= (sPxi)xi + xiisfhi)
=
(hmiXx)X2 + X] (AOT2X2) = ft(m.i + wii)xtxa,
(note that S^ acts only on xi? and S(2^ acts only on xi\\ this notation may not be
elegant, but it. does the job). So m (the quantum number for the composite system)
is just mi + mi;
ff: m l;
f|: m 0;
If: m 0;
44,^ m -1.
At first glance, this doesn't look right: m is supposed to advancein integer
steps,
from \342\200\224
$ to +s, so it appears that s \342\200\224
1 \342\200\224but there is an \"extra\" state
with m = 0. One way
to untangle this problem is to apply the lowering operator,
S_ = s!P + 5l2) to the state ft,- usin\302\247 Equation 4.146:
=.(* ^) t + r i)
(\302\276
- mt + u>.
r in) -tt 1
<
110) =\302\261tu + m. 5 = 1 (triplet). [4-177]
..11-1) = 44
(Asa check, try applying the lowering operator to |1 0); what should you get? See
= i(U
- 5 = 0 (singlet). [4.178]
{|00) W}
(If you apply the raising or lowering operator to this state,you'll get zero. See
Problem 4.34(b).)
186 Chapter 4 Quantum Mechanicsin ThreeDimensions
I claim, then, that the combination of two spin-1/2 particles- can carry a total
S2- = _ . s(2)_
(S(1) + S(2))
.
(S(1) + SC2>) (S(l)f + (vS(2))2 + 2S(1) [4>179fj
s(1) \342\200\242
smm)
= (sp t)(sf $.+ ism txsf ^) + $n t)(^2) i).
h2
=
-(24t-H)-
Similarly,
= it)-
sP .s^(w
^-(2 n-
It follows that
,2 i fc2
fe2 i
. S^IOO) = - -
3^2
[4.181]
-j jj=C24rt
S^ -2 n + ;t)
t4-
-^100).
Returning
to Equation 4.179 (and using Equation 4.142), we conclude that
0 - (
/sn2 3h2 H2\\
= 2/72|10),
9
S2\\\\ 0}
\342\200\224
+
\342\200\224
+2\342\200\224 ] |10) [4.182]
Problem 4,34(6)*)
Section 4.4: Spin, 187
What we have just done (combining spin 1/2 with spin 1/2 to get spin1 and
spin 0) is the simplest example of a larger problem:If you combine spin s\\ with
spin s.2.,. what total spins s can you get?4\"The answer41 is that you get every spin
from (s-i + 52) down to (s\\ \342\200\224
.\302\276)\342\200\224of
(\302\276
\342\200\224
i'l), if >
\302\2532 integer
\302\2431\342\200\224in steps:
[4.184]
(Roughly speaking, the highest total spin occurswhen the individual spins are
\\s m)
= 1*1 \302\2731 \302\253a>
> 1\302\276 [4.185]
\302\243 C%%n
m.\\ -\\-m2=m
(because the z components add, the only composite states that contribute are those
for which mi + m.2 \342\200\224
in): Equations 4.177 and 4,178 are special cases of this
general form, 1/2 (I usedthe informalnotationt =
with si = s% =
\\)-> It ^
\342\200\224
1\302\276(-5)).)-
The constants Cs,n^m2m ar? called Clebsch-Gotdan coefficients. A few
of the simplest cases are listed in Table 4.8.42For example, shaded column the of
the 2x1 table tells us that
|30}
= -L|2 1)|1 -1) + |20)|10}
y\302\247
+
-^|2- 1)|11).
In particular,if two
particles (of spin 2 and spin 1) are at rest in a box, and the total
I say spins, for simplicity, but either one (or both) could just as Well be orbital angular
momentum (for which, however, we would use the letter I).
41For a must look in a more advanced instance, ClaudeCohen-Tannoudji,
text; see, for
proof you
Bernard Diu, and Franck Laloe,: Quantum Mechanics, (Wiley, York, 1977), Vol. 2, Chapter X.
New
1/5).Notice that the probabilities add up to 1 (the sum of the squaresof any column
\\si Ml)\\s2m1)
= V C%$ \\sm)^ [4.186]
^Problem 4.34
(a) Apply S_ to 110) (Equation4-177), and confirm that you get V2h\\l -1).
(b) Apply S\302\261 to )00} (Equation 4.178), and confirm that you get zero,
appropriate eigenvalue.
Problem 4.36
probability of each?
Problem 4.37 Determine the commutatorof S2 with s|1} (where S = S<1}+ S(2)).
Generalize your result to show that
obvious inference from Equation4.187 that the sum S^ +SC9- does commute with
S2, which is a specialcase something we of already knew (see Equation 4.103).
EA = (n+3/2)h(D. [4.189]
**Problem4.40
(a) Prove the three-dimensional vi'rial theorem:
(b) Apply the virial theoremto thecaseof hydrogen, and show that
(T) = (V)
=
En/2. [4.192]
* ^Problem
\342\226\240*\342\226\240\342\226\2404.41 [Attempt this problem only if you are familiar with vector
V-J = -|-|*|2, at
[4.194]
h
:re-r/asm0$.
L = m / (r x J) d?r.
Use this to calculate Lz for the state #2.11, and commentonthe result.
* *. ^Problem 4.42 The (time independent) momentum space wave function in three
(a) Findthe momentum space wave function for the ground state of hydrogen
(Equation4.80).Hint: Use spherical coordinates, setting the polar axis along
the direction of p. Do the 0
integral first. Answer:
1 / o \\ 3/2 1
**) = [4197]
;f(t) [i + iaP,nn-
(d) What is the expectation value of the kineticenergy in this state? Express
your answer as a multiple of Ei, and check that it is consistent with the
virial theorem (Equation4.191).
Problem4.43
(VO for hydrogen in the state n
\342\200\224
(a) Construct the spatial wave function 3,
I = 2, m \342\200\224
1. Express your answer as a function of r, 0, </>, and a (the
Bohr radius) only\342\200\224no
other variables (p, z, etc.) or functions (Y, v, etc.), or
constants (A, cq, etc.), or derivatives, allowed {n is okay, and e, and 2, etc.).
(b) Check that this wave function is properly normalized, by carrying out the
appropriateintegrals over r, 0, and \302\242.
(c)
Find the expectation value of rs in this state. For what range of s (positive
and negative) is the result finite?
Problem4.44
for hydrogen in the state = 3, m = 3.
\342\200\224
(a) Construct the wave function n 4,-. /
Express your answer as a functionof the spherical coordinates r, 6.., and <fi.
(b) Find the expectation value of r in this state.(As always, look up any nontrivial
integrals.)
(c) If you could somehow measure the observableL2 + L2, on an atom in this
state, what value (or values) could you get, and what is the probabilityof
each?
(a) First calculatethe exact answer, assuming the wave function (Equation 4,80)
is correct all the way
down to r = 0. Let b be the radiusof the nucleus.
Further Problems for Chapter 4 193
(b) Expand your result as a powerseries in the small number e = 2b/d, and
show that the lowest-order term is the cubic: P \302\253*
(4/3)(fr/a)3. This should
be a suitable approximation, provided that h <^a (which it is).
(c) Alternatively, we might assume that ^.(r) is essentially constant over the
(tiny) volume of the nucleus,so that P -\302\276
(4/3)nb3\\f(0)\\^. Check that you
get the same answer this way.
Problem4.46
(a) Use the recursion formula (Equation 4.76) to confirm that when I n \342\200\224
\342\200\224 1
p ,
^n(f?-i)
, _ a!
\342\200\224 vn\342\200\224^0\342\200\224rjna
Nnr e ,
(c) Show that the \"uncertainty\" in r (oy) is (r)l\\/2n + 1 for such states. Note
that the fractional spread in r decreases, with increasing n (in this sensethe
system\"begins to; look classical,\" with identifiable circular \"orbits,\" for large
n). Sketchtheradial wave functions for several values of n, to illustrate this
point.
pairs {n,,n/} that yield the same X. For example, {6851, 6409} and {15283, 11687}
will do it, but you're not allowed to use those!
(b) If you measured Sz on this electron, values what could you get, and what is
the probabilityof each? is the expectation What value of S^?
(c) If you measured S'x: on this; electron, what values could you get, and what is
(d) If you measured Sy on this electron, what values could you get, and what is
the probability of each? What is the expectationvalue of Sy?
configuration (Equation 4.178). Let .\302\276 be the component of the spin angular momentum
of particlenumber 1 in the direction defined by the unit vector a. Similarly, let
angular momentum in the directionb. Show
\342\226\240
$1 be the component of 2's that
= [4.198]
(S^sf)
-jcosB,
* * ^Problem 4.51
coefficients
(a) Work out the Clebsch-Gofdari for the case s\\ = 1/2, \342\200\224
.\302\276
Is m) - A\\\\ 1)1\302\276 (m
-
i)> + B\\\302\261 (-^))1\302\276 (m + |)),
refer back to Equation 4.136 and the line before Equation 4.147. Answer:
2*2 + l V
2^2+1
Problem 4.52 Find the matrix representing Sx for a particleof spin3/2 (using, as
* * ^Problem 4.53 Work out the spin matrices for arbitrary spin s, generalizing spin 1/2
(Equations 4.145 and 4.147), spin 1 (Problem
4,31), and
spin 3/2 (Problem 4.52).
Answer:
*s 0 0 0\\
0 s --1 0 0
S,= h 0 0 s-2 0
^ 0 0 ~s)
/0 hs 0 0 0 0 \\
h 0 ^-1 0 0 0
0 bS-\\ 0 bs-2 o 0
h 0
,or \342\200\224
\342\200\224
0 ^-2 0 0 0
2 ; fc .
\342\226\240I T '. '. I I
0 0 0 0 0
b~s+i
\\o
0 0 0 b-s+i o )
/0 -ihs 0 0 w . .
0 0 \\
ibs 0 -ih -1 0 0 0
0 ibs~] 0 -th -2 0
\342\200\242\342\200\242\342\200\242
0
S, =
H
0 0 ibs~ 2 0 0 0
2 \342\200\242 \342\226\240 \342\200\242 \342\200\242\342\226\240
;\342\226\240 ;
0 0 0 0 . . . 0 -ik-s+i
0 0 0 i-fe-j+1 0
W
where
= + ,/)Cy+l-7)
bj V^
*. * ^Problem 4.54 Work out the normalization factor for the spherical harmonics, as
follows. From Section 4.1.2 we know that
=
Yff B?t?m*PFiGQs6y,
relation giving B^
in terms of
*Z?zm.
Solve it by induction on m to get B\342\204\242
up to
an overall constant,CQ). Finally, use the resultof Problem4.22 fix the constant. to
You may find the following formula for the derivativeof an associated Legendre
function useful:
dPm
- x1) \342\200\224l- = - [4.199]
(1
dx Vl-x2^1 l mxP'r.
position state
R2l
(7i/3 3fx+ + VV3YU-) \342\226\240
(a) If you measured the orbital angular momentum squared (L2), what values
Let J L
\342\226\240= + S be the total angular momentum.
(e) If you
measured J2, what values might you get, and what is the probability
of each?
(h) If you measured boththe z component of thespin and the distance from the
origin (note that these are compatibleobservables), what is the probability
density for finding the particle with spin up and at radius r?
* ^Problem
*\342\226\240 4.56
X, = ei\302\253x^2x [4-200]
(b) Construct the (2 x 2) representing matrix rotation by 180\302\260 about the jc-axis,
and show that it converts \"spin up\" (x+) into \"spin down\" ^s
(x\342\200\224)\302\273 y\302\260u
would expect.
(c) Construct the matrix representing rotation by 90Qaboutthe y-axis, and check
(d) Construct the matrix representing rotation by ,360\302\260 about the z-axis. If the
answer is not quite what you expected, discuss its implications.
Let a be some convenient constant with the dimensions of length (the Bohr radius,
say, if we're talking about hydrogen), and define the operators
Pi = -
-L
=\342\226\240
q\\ + ;
-\302\276[px (fi/a^y] ;
[* (a2/H)py]
12 s 7!
~ =
[px. + (hia2)y\\ \342\226\240
[x (a2(hyPy\\
r P2
7i
[4.203]
dt
thereis no
problem; the classical Hamiltonian45 is
1
H = --(p
- qAY + 9
q<pt [4.204]
2m
[4.205]
45.
See, for example,H. Goldstein, C. P. Poole, and J. L. Safko, ClassicalMechanics, 3rd ed.,
(Prentice Hall, Upper Saddle River, Nj, 2002), page 342.
Further Problems for Chapter 4 199
(a) Showthat
~^ = -<(p-?A)}. [4,206]
at m
m^P- at
= q(E) +
^- {(p xB-Bxp))-%x
2m- m
B)). [4.207]
(c) In particular, if the fields E and B are uniform over the volume of the wave
packet, show that
sothe expectation value of (v) moves according to the Lorentz force law, as
we would expect from Ehrenfesf s theorem.
A: = - A and = Kz2,
-y(xJ cp
(b) Find the allowed energies^ for a particleof mass m and charge q, in these
fields. Answer:
^Problem
.\342\226\240* 4.61 [Refer to Problem 4.59 for background.]In classical
the
electrodynamics
potentials A and p are hot uniquely determined;47the physicalquantities
are the fields, E and B.
9A
=
\302\242/ 0-.-, A' = A + VA [4.210]
at
(where A is an arbitrary real functionof position and time) yield the same
fields as tp and A. Equation 4.210 is called a. gauge transformation, and the
theory is saidto be gauge invariant..
W'=e'>A/% [4.211]
48That is to say, (r)* d(r)/dt, etc. are unchanged. BecauseA depends on position, (p) (with p
represented by the operator (ft/i)V) does change, but as we found in Equation 4.206, p does not
represent the mechanical momentum (mv) in this context (in lagrangian mechanics it is so-called
canonical momentum).