Solutions Quantum Mechanics 2 PDF
Solutions Quantum Mechanics 2 PDF
Solutions Quantum Mechanics 2 PDF
Similarly
Z Z ! Z
X X X
dx ψ ∗ Qψ = dx a∗i u∗i Q aj u j =
∗
ai aj dxu∗i Quj
i j ij
X Z X X
= a∗i aj qj dx u∗i uj = |ai |2 qi = Pi qi
ij i i
which is by definition the expectation value of Q.
2.22 Show that for any three operators A, B and C, the Jacobi identity holds:
[A, [B, C]] + [B, [C, A]] + [C, [A, B]] = 0. (2.1)
Soln: We expand the inner commutator and then use the rule for expanding the commutator of a
product:
[A, [B, C]] + [B, [C, A]] + [C, [A, B]] = [A, BC] − [A, CB] + [B, CA] − [B, AC] + [C, AB] − [C, BA]
= [A, B]C + B[A, C] − [A, C]B − C[A, B] + [B, C]A + C[B, A] − [B, A]C
− A[B, C] + [C, A]B + A[C, B] − [C, B]A − B[C, A]
= 2 ([A, B]C − C[A, B] − [A, C]B + B[A, C] + [B, C]A − A[B, C])
= 2 ABC − BAC − CAB + CBA − ACB + CAB + BAC − BCA + BCA
− CBA − ABC + ACB = 0
3.9∗ By expressing the annihilation operator A of the harmonic oscillator in the momentum rep-
resentation, obtain hp|0i. Check that your expression agrees with that obtained from the Fourier
transform of r
1 −x2 /4ℓ2 h̄
hx|0i = e , where ℓ ≡ . (3.1)
(2πℓ2 )1/4 2mω
Soln: In the momentum representation x = ih̄∂/∂p so [x, p] = ih̄∂p/∂p = ih̄. Thus from Problem
3.7
x ℓ ℓp h̄ ∂
A= +i p =i +
2ℓ h̄ h̄ 2ℓ ∂p
ℓp h̄ ∂u0 2 2 2
0 = Au0 ⇒ u0 = − ⇒ u0 (p) ∝ e−p ℓ /h̄
h̄ 2ℓ ∂p
1 In the most elegant formulation of qantum mechanics, this last result is the basic postulate of the theory, and
one derives other rules for the physical interpretation of the qn , an etc. from it – see J. von Neumann, Mathematical
Foundations of Quantum Mechanics.
2
3.10 Show that for any two N × N matrices A, B, trace([A, B]) = 0. Comment on this result in
the light of the results of Problem 3.7 and the canonical commutation relation [x, p] = ih̄.
Soln: X X X X X
trace(AB) = (AB)jj = Ajk Bkj ; trace(BA) = (BA)kk = Bkj Ajk = Ajk Bkj
j jk k kj kj
These two expressions differ only in the order Pof the summations over j and k. We can reverse
one order and thus show them to be equal iff jk |Ajk Bkj |2 is finite, which it will be for finite-
dimensional matrices with finite elements. Hence x and p can only be represented by infinite matrices
with elements that are not square summable, (which the matrices of Problem 3.7 are not).
3.11∗ A Fermi oscillator has Hamiltonian H = f † f , where f is an operator that satisfies
f2 = 0 ; f f † + f † f = 1. (3.2)
2
Show that H = H, and thus find the eigenvalues of H. If the ket |0i satisfies H|0i = 0 with
h0|0i = 1, what are the kets (a) |ai ≡ f |0i, and (b) |bi ≡ f † |0i?
In quantum field theory the vacuum is pictured as an assembly of oscillators, one for each
possible value of the momentum of each particle type. A boson is an excitation of a harmonic
oscillator, while a fermion in an excitation of a Fermi oscillator. Explain the connection between
the spectrum of f † f and the Pauli principle.
Soln:
H 2 = f † f f † f = f † (1 − f † f )f = f † f = H
Since eigenvalues have to satisfy any equations satisfied by their operators, the eigenvalues of H
must satisfy λ2 = λ, which restricts them to the numbers 0 and 1. The Fermi exclusion principle
says there can be no more than one particle in a single-particle state, so each such state is a Fermi
oscillator that is either excited once or not at all.
stationary, all its energy is invested in the spring and the mass suddenly gets poorer. We expect
quantum mechanics to yield a probability distribution for the changes in E that is similar to the
average of the classical changes over phase. Of course classical mechanics averages probabilities for
states of even and odd parity because it doesn’t recognise the parity of a state.
3.13∗ P is the probability that at the end of the experiment described in Problem 3.12, the
oscillator is in its second excited state. Show that when f = 12 , P = 0.144 as follows. First show
that the annihilation operator of the original oscillator
A = 21 (f −1 + f )A′ + (f −1 − f )A′† ,
(3.3)
where A′ and A′† are the annihilation and creation operators P of the final oscillator. Then writing
the ground-state ket of the original oscillator as a sum |0i = n cn |n′ i over the energy eigenkets
of the final oscillator, impose the condition A|0i = 0. Finally use the normalisation of |0i and the
orthogonality of the |n′ i. What value do you get for the probability of the oscillator remaining in
the ground state?
Show that at the end of the experiment the expectation value of the energy is 0.2656h̄ω. Explain
physically why this is less than the original ground-state energy 12 h̄ω.
This example contains the physics behind the inflationary origin of the Universe: gravity ex-
plosively enlarges the vacuum, which is an infinite collection of harmonic oscillators (Problem 3.11).
Excitations of these oscillators correspond to elementary particles. Before inflation the vacuum is
unexcited so every oscillator is in its ground state. At the end of inflation, there is non-negligible
probability of many oscillators being excited and each excitation implies the existence of a newly
created particle.
Soln: From Problem 3.6 we have
mωx + ip mf 2 ωx + ip
A≡ √ A′ ≡ p
2mh̄ω 2mh̄f 2 ω
x iℓ fx iℓ
= + p = + p
2ℓ h̄ 2ℓ f h̄
Hence
f 2iℓ 1 f
A′ + A′† = x A′ − A′† = f p so A = (A′ + A′† ) + (A′ − A′† )
ℓ f h̄ 2f 2
X
0 = A|0i = 21 (f −1 + f )ck A′ |k ′ i + (f −1 − f )ck A′† |k ′ i
k
Xn
−1
√ √ o
= 1
2 (f + f ) kck |k − 1′ i + (f −1 − f ) k + 1ck |k + 1′ i
k
Multiply through by hn′ |:
√ √
0 = (f −1 + f ) n + 1cn+1 + (f −1 − f ) ncn−1 ,
which is a recurrence relation from which all non-zero cn can be determined in√ terms of c0 . Put
c0 = 1 and solve for the cn . Then evaluate S ≡ |cn |2 and renormalise: cn →Pcn / S.
The probability of remaining in the ground state is |c0 |2 = 0.8. hEi = n |cn |2 (n + 21 )h̄f 2 ω. It
is less than the original energy because of the chance that energy is in the spring when the stiffness
is reduced.
3.14∗ In terms of the usual ladder operators A, A† , a Hamiltonian can be written
H = µA† A + λ(A + A† ). (3.4)
What restrictions on the values of the numbers µ and λ follow from the requirement for H to be
Hermitian?
Show that for a suitably chosen operator B, H can be rewritten
H = µB † B + constant. (3.5)
†
where [B, B ] = 1. Hence determine the spectrum of H.
4
We have to choose the plus sign in order to achieve consistency with our previously established value
of Emin ; thus finally p
En = − 12 ǫ + (n + 21 ) ǫ2 − λ2
3.17∗ This problem is all classical emag, but it gives physical insight into quantum physics. It is
hard to do without a command of Cartesian tensor notation. A point charge Q is placed at the
origin in the magnetic field generated by a spatially confined current distribution. Given that
Q r
E= (3.8)
4πǫ0 r3
and B = ∇ × A with ∇ · A = 0, show that the field’s momentum
Z
P ≡ ǫ0 d3 x E × B = QA(0). (3.9)
Write down the relation between the particle’s position and momentum and interpret this relation
physically in light of the result you have just obtained.
Hint: write E = −(Q/4πǫ0 )∇r−1 and B = ∇ × A, expand the vector triple product and
integrate each of the resulting terms by parts so as to exploit in one ∇R· A = 0 and in
H the other
∇2 r−1 = −4πδ 3 (r). The tensor form of Gauss’s theorem states that d3 x ∇i T = d2 Si T no
matter how many indices the tensor T may carry.
Soln: In tensor notation with ∂i ≡ ∂/∂xi
Q
Z
P=− d3 x ∇r−1 × (∇ × A)
4π
becomes
Q
Z X
Pi = − d3 x ǫijk ∂j r−1 ǫklm ∂l Am
4π
jklm
Q
Z X
=− d3 x ǫkij ǫklm ∂j r−1 ∂l Am
4π
jklm
Q
Z X
=− d3 x (δil δjm − δim δjl )∂j r−1 ∂l Am
4π
jlm
Q
Z X
=− d3 x (∂j r−1 ∂i Aj − ∂j r−1 ∂j Ai )
4π j
In the first integral we use the divergence theorem to move ∂j from r−1 to ∂i Aj and in the second
integral we move the ∂j from Aj to ∂j r−1 . We argue that the surface integrals over some bounding
sphere of radius R that arise in this process vanish as R → ∞ because their integrands have explicit
R−2 scaling, so they will tend to zero as R → ∞ so long as A → 0 no matter how slowly. After this
has been done we have Z
Q
Z
Pi = − − d3 x (r−1 ∂j ∂i Aj + d3 x ∂j ∂j r−1 Ai
4π
The first integral vanishes because ∂j Aj = ∇ · A = 0 and in the second integral we note that
∂j ∂j r−1 = ∇2 r−1 = −4πδ 3 (x), so evaluation of the integral is trivial and yields the required expres-
sion.
The physical interpretation is that when we accelerate a charge, we have to inject momentum
into the emag field that moves with the charge. So the conserved momentum p associated with the
particle’s motion includes both the momentum in the particle and that in the field, which we have
just shown to be QA(x), where x is the particle’s location. Thus p = mẋ + QA. The particle’s
kinetic energy is H = 21 mẋ2 = (p − QA)2 /2m.
3.18∗ From equation (3.58) show that the the normalised wavefunction of a particle of mass m
that is in the nth Landau level of a uniform magnetic field B is
2 2
rn e−r /4rB
e −inφ
hx|ni = √ n+1
, (3.10)
2(n+1)/2 n! π rB
6
p
where rB = h̄/QB. Hence show that the expectation of the particle’s gyration radius is
√ rB
hrin ≡ hn|r|ni = 2{(n + 12 )(n − 21 ) × · · · × 12 } . (3.11)
n!
Show further that
δ ln hrin 1
≃ (3.12)
δn √ 2n
and thus show that in the limit of large n, hri ∝ E, where E is the energy of the level. Show that
this result is in accordance with the correspondence principle. R
2 2
by evaluating dr r r2n e−r /2rB (n+r2 /2rB 2
R
Soln: We obtain the required expression R ∞ forz hx|ni ) dφ
with the help of the definition z! = 0 dx x e−x . Then
2 2 √
dr r2 r2n e−r /2rB 2π 23/2 rB ∞ dz n+1/2 −z
R
2rB
Z
hrin = 2n+2 = z e = (n + 12 )!
2n+1 n! π rB n! 0 2 n!
√
where z ≡ r/ 2rB .
(n + 23 )!n! n + 32 3
1 + 2n
δ lnhrin 1
= hrin+1 − hrin = ln 1 = ln = ln 1 ≃
δn (n + 2 )!(n + 1)! n+1 1+ n 2n
Summing this expression over several increments we have that the overall change ∆ ln(hrin ) =
∆ ln(n1/2 ), so for n ≫ 1, hrin ∼ n1/2 ∼ E 1/2 . Classically
2
mv 2
QB
= QvB ⇒ E = 21 mv 2 = 21 m r2
r m
so r ∝ E 1/2 as in QM.
3.19 A particle of charge Q is confined to move in the xy plane, with electrostatic potential φ = 0
and vector potential A satisfying
∇ × A = (0, 0, B). (3.13)
Consider the operators ρx , ρy , Rx and Ry , defined by
1
ρ= êz × (p − QA) and R = r − ρ, (3.14)
QB
where r and p are the usual position and momentum operators, and êz is the unit vector along B.
Show that the only non-zero commutators formed from the x- and y-components of these are
2 2
[ρx , ρy ] = irB and [Rx , Ry ] = −irB , (3.15)
2
where rB = h̄/QB.
The operators a, a† , b and b† are defined via
1 1
a= √ (ρx + iρy ) and b= √ (Ry + iRx ). (3.16)
2 rB 2 rB
Evaluate [a, a†] and [b, b† ]. Show that for suitably defined ω, the Hamiltonian can be written
H = h̄ω a† a + 12 .
(3.17)
Given that there exists a unique state |ψi satisfying
a|ψi = b|ψi = 0, (3.18)
what conclusions can be drawn about the allowed energies of the Hamiltonian and their degeneracies?
What is the physical interpretation of these results?
7
Soln: We have
1 1
ρx = − (py − QAy ) ; ρy = (px − QAx )
QB QB
so
1 1
[ρx , ρy ] = − [(py − QAy ), (px − QAx )] = ([py , Ax ] + [Ay , px ])
(QB)2 QB 2
−ih̄ ∂Ax ∂Ay ih̄ 2
= − = = irB
QB 2 ∂y ∂x QB
Also
1 1 ih̄ 2
[ρx , x] = − [py , x] = 0 ; [ρx , y] = − [py , y] = = irB
QB QB QB
1 2
[ρy , x] = [px , x] = −irB ; [ρy , y] = 0
QB
So
2
[ρx , Rx ] = [ρx , x − ρx ] = 0 ; [ρx , Ry ] = [ρx , y − ρy ] = irB (1 − 1) = 0
and similarly [ρy , Ri ] = 0. Finally
2
[Rx , Ry ] = [x − ρx , y − ρy ] = −[ρx , y] − [x, ρy ] + [ρx , ρy ] = −irB
Now
1 1
[a, a† ] = 2 [ρx + iρy , ρx − iρy ] = − 2 2i[ρx , ρy ] = 1
2rB 2rB
1 1
[b, b† ] = 2 [Ry + iRx , Ry − iRx ] = 2 2i[Rx , Ry ] = 1
2rB 2rB
1 1
a† a = 2 (ρx − iρy )(ρx + iρy ) = 2 (ρ2x + ρ2y + i[ρx , ρy ])
2rB 2rB
2
1 |p − QA|2
1 |p − QA| 2
= 2 2
− rB = − 12
2rB (QB) 2h̄ QB
But
|p − QA|2 QB
H= = (a† a + 21 )h̄ = (a† a + 21 )h̄ω,
2m m
where ω ≡ QB/m.
Given a|ψi = 0, we have H|ψi = (a† a + 12 )h̄ω|ψi = 12 h̄ω|ψi so |ψi is a stationary state of energy
1 †
2 h̄ω. Given any stationary state |Ei of energy |Ei we have that a |Ei is a stationary state of energy
E + h̄ω:
Ha† |Ei = h̄ωa† (aa† + 21 )|Ei = h̄ωa† (a† a + [a, a† ] + 21 )|Ei = a† (H + h̄ω)|Ei = (E + h̄ω)a† |Ei,
so the energies are h̄ω × ( 21 , 32 , 52 , . . .).
Since [ρi , Rj ] = 0 we have [a, b] = [a, b† ] = 0 etc, and it follows that [b† , H] = 0. Hence b† |Ei is
a stationary state of energy E.
To show that it is distinct from |Ei, we define the number operator N ≡ b† b. N |ψi = 0 because
b|ψi = 0. Given that N |ni = n|ni, we have that
N (b† |ni) = b† (bb† |ni) = b† (b† b + [b, b† ])|ni = b† (N + 1)|ni = (n + 1)b† |ni
Thus by repeatedly applying b† to |ψi we can create distinct states of energy 12 h̄ω, and by applying
it to a† |ψi we can make distinct states of energy 32 h̄ω, etc. Hence all energy levels are infinitely
degenerate.
Physically the energy levels are degenerate because there are an infinite number of possible
locations of the gyrocentre of a particle that is gyrating with a given energy.
3.20∗ Equation (2.87) gives the probability current density associated with a wavefunction. Show
that the flux given by this expression changes when we make a gauge change ψ → ψ ′ = eiQΛ/h̄ ψ.
Why is this state of affairs physically unacceptable?
8
1 1 1
(hψ|xihx|p|ψi + hx|ψi(hx|p|ψi)∗ ) = ψ 2 (−ih̄∇ψ) + ψ ∗ (−ih̄∇ψ)∗
ℜ (hψ|xihx|p|ψi) =
m 2m 2m
which agrees with the definition of J. The principle of minimal coupling requires us to replace p by
p − QA, so we redefine J to be
1
J(x) = ℜ (hψ|xihx|(p − QA)|ψi) .
m
From H we recognise (p − QA)/m as the particle’s velocity, and clearly this is what determines the
current. Converting back to the position representation, we find
1
J= (hψ|xihx|(p − QA)|ψi + hx|ψi(hx|(p − QA)|ψi)∗ )
2m
1
= (ψ ∗ (−ih̄∇ψ − QAψ) + ψ(−ih̄∇ψ − QAψ)∗ )
2m
1
−ih̄(ψ ∗ ∇ψ − ψ∇ψ ∗ ) − 2QA|ψ|2
=
2m
In amplitude-phase format we now have
|ψ|2
J= (h̄∇θ − QA).
m
This is clearly invariant when we add QΛ/h̄ to θ and Q∇Λ to A.
In plane polars,
∂ 1 ∂
∇= , ,
∂r r ∂φ
so we obtain the required expression for J of the Landau level when we take |ψ| and θ from (3.10)
and use A = 12 Brêφ .
9
The figure shows the current density as a function of radius for the first six Landau levels. The
ground state is
√ in a class by itself. For the excited states the characteristic radius r increases with
n roughly as n as predicted by classical physics, while the peak magnitude of the flow is almost
independent of n because J ∝ ω|ψ 2 | and ω ∼ constant and |ψ | 2 ≃ 1/r as expected from classical
mechanics.
3.21∗ Determine the probability current density associated with the nth Landau ground-state
wavefunction (3.67) (which for n = 4 is shown on the front cover). Use your result to explain in as
much detail as you can why this state can be interpreted as a superposition of states in which the
electron gyrates around different gyrocentres. Hint: adapt equation (3.21).
2
Why is the energy of a gyrating electron incremented if we multiply the wavefunction e−(mω/4h̄)r
n n n n
by v = (x − iy) but not if we multiply it by u = (x + iy) ?
Soln: The wavefunction of the nth state in the ground-state level differs from that of the first state
2
of the nth excited level only in having un instead of v n in front of e−uv/4rB . Consequently these
th
wavefunctions have the same amplitudes but where the phase of the n excited state is θ = −nφ,
that of the nth ground state is θ = nφ. Since J ∝ (h̄∇θ − QA) it follows that the probability current
is given by equation (3.21) with a minus sign in front of the n:
2 2
h̄r2n−1 e−r /2rB r2
J(r) = − n+1 2n+2 2 − n êφ ,
2 πn!mrB 2rB
√
J now changes sign at r = 2nrB : at small r it flows in the opposite sense to the current of the nth
excited level (and the opposite sense to that of classical gyrations), but flows in the same sense at
large r. This behaviour
√ arises naturally if you arrange clockwise vortices of radius ∼ rB around a
circle of radius 2nrB . So √ whereas in the nth excited level the particle is moving counter-clockwise
around a path of radius√∼ nrB , in this state it is moving on a circle of radius ∼ rB that is offset
from the centre by ∼ nrB . The direction of the offset is completely uncertain, so we have to
imagine a series of offset circles. These clearly generate an anti-clockwise flow around the origin:
When the ground-state ψ(r) is multiplied by either un or v n , a factor rn e±inφ is added. This
factor moves outwards the radius rm at which |ψ|2 peaks in the same way regardless of whether
we multiply by un or v n . But our expressions for J(r) show that when we multiply by un we shift
rm out to the radius near which J = 0 so there is not much kinetic energy. Had we multiplied by
v n , J and the kinetic energy would have been large at rm . This asymmetry between the effects of
multiplying by einφ or e−inφ arises because A always points in the direction êφ .
3.22∗ In classical electromagnetism the magnetic moment of a planar loop of wire that has area
A, normal n̂ and carries a current I is defined to be
µ = IAn̂. (3.22)
Use this formula and equation (3.21) to show that the magnetic moment of a charge Q that is in
a Landau level of a magnetic field B has magnitude µ = E/B, where E is the energy of the level.
Rederive this formula from classical mechanics.
Soln: Eq. (3.21) gives the probability current density J. We multiply it by Q to get the electrical
current density QJ. The current dI = QJ(r)dr through the annulus of radius r contributes magnetic
moment dµ = πr2 dI. Summing over annuli we obtain
Z ∞ Z ∞
r2
2 Qh̄ 2n+1 −r 2 /2rB
2
µ = Qπ dr r J(r) = − n+1 2n+2 dr r e n+ 2
0 2 n!mrB 0 2rB
Z ∞
h̄Q h̄Q h̄Q
=− dz z n e−z (n + z) = − [n n! + (n + 1)!] = (n + 12 )
2n!m 0 2n!m m
10
4.3∗ We have shown that [vi , Jj ] = i k ǫijk vk for any operator whose components vi form a
P
vector. The expectation value of this operator relation in any state |ψi is then hψ|[vi , Jj ]|ψi =
i k ǫijk hψ|vk |ψi. Check that with U (α) = e−iα·J this relation is consistent under a further rotation
P
|ψi → |ψ ′ i = U (α)|ψi by evaluating both sides separately.
Soln: Under the further rotation the LHS → hψ|U † [vi , Jj ]U |ψi. Now
U † [vi , Jj ]U = U † vi Jj U − U † Jj vi U = (U † vi U )(U † Jj U ) − (U † Jj U )(U † vi U )
X X
= [Rik vk , Rjl Jl ] = Rik Rjl [vk , Jl ].
kl kl
Similar |ψi → U |ψi on the RHS yields
X
i Rkm ǫijk hψ|vm |ψi.
km
We now multiply each side by Ris Rjt and sum over i and j. On the LHS this operation yields
[vs , Jt ]. On the right it yields
X X
i Ris Rjt Rkm ǫijk hψ|vm |ψi = i ǫstm hψ|vm |ψi,
ijkm m
By considering the form taken by R for infinitesimal φ calculate from √ R the matrix J z that appears
√
in R(φ) = exp(−iJ Jz φ). Introduce new coordinates u1 ≡ (−x+iy)/ 2, u2 = z and u3 ≡ (x+iy)/ 2.
Write down the matrix M that appears in u = M ·x [where x ≡ (x, y, z)] and show that it is unitary.
Then show that
J z′ ≡ M · J z · M† . (4.2)
is identical with Sz in the set of spin-one Pauli analogues
0 1 0 0 −i 0 1 0 0
1 1
Sx = √ 1 0 1 , Sy = √ i 0 −i , Sz = 0 0 0 . (4.3)
2 2
0 1 0 0 i 0 0 0 −1
Write down the matrix J x whose exponential generates rotations around the x axis, calculate J x′
by analogy with equation (4.2) and check that your result agrees with Sx in the set (4.3). Explain
as fully as you can the meaning of these calculations.
Soln: For an infinitesimal rotation angle δφ to first order in δφ we have
1 −δφ 0
1 − iJ
Jz δφ = R(δφ) = δφ 1 0
0 0 1
comparing coefficients of δφ we find
0 −1 0
Jz = i 1 0 0
0 0 0
In components u = M · x reads
u1 −1 i √0 x
1
u2 = √ 0 0 2y
2
u3 1 i 0 z
so M is the matrix above. We show that M is unitary by calculating the product MM† . Now we
have
−1 i √0 0 −i 0 −1 0 1
J z′ = 21 0 0 2 i 0 0 −i √0 −i
1 i 0 0 0 0 0 2 0
−1 i √0 −1 0 −1 2 0 0
= 12 0 0 2 −i 0 i = 21 0 0 0
1 i 0 0 0 0 0 0 −2
Similarly, we have
0 0 0
J x = i 0 0 −1
0 1 0
so
−1 i √0 0 0 0 −1 0 1
J x′ = 12 0 0 2 0 0 −i −i √0 −i
1 i 0 0 i 0 0 2 0
√
−1 i √0 0 0√ 0 0
√ 2 √0
= 12 0 0 2 0 −i 2 0 = 12 2 √0 2
1 i 0 1 0 1 0 2 0
These results show that the only difference between the generators of rotations of ordinary 3d
vectors and the spin-1 representations of the angular-momentum operators, is that for conventional
vectors we use a different coordinate system than we do for spin-1 amplitudes. Apart from this, the
three amplitudes for the spin of a spin-1 particle to point in various directions are equivalent to the
components of a vector, and they transform among themselves when the particle is reoriented for
the same reason that the rotation of a vector changes its Cartesian components.
12
4.5 Determine the commutator [J Jx′ , J z′ ] of the generators used in Problem 4.4. Show that it is
′ ′
equal to −iJ
Jy , where J y is identical with Sy in the set (4.3).
Soln: This is straightforward extension of the previous problem.
4.6∗ Show that if α and β are non-parallel vectors, α is not invariant under the combined rotation
R(α)R(β). Hence show that RT (β)RT (α)R(β)R(α) is not the identity operation. Explain the
physical significance of this result.
Soln: R(α)α = α because a rotation leaves its axis invariant. But the only vectors that are
invariant under R(β) are multiples of the rotation axis β. So R(β)α is not parallel to α.
If RT (β)RT (α)R(β)R(α) were the identity, we would have
RT (β)RT (α)R(β)R(α)α = α ⇒ R(β)R(α)α = R(α)R(β)α ⇒ R(β)α = R(α)(R(β)α)
which would imply that R(β)α is invariant under R(α). Consequently we would have R(β)α = α.
But this is true only if α is parallel to β. So our original hypothesis that RT (β)RT (α)R(β)R(α) = I
is wrong. This demonstrates that when you rotate about two non-parallel axes and then do the
reverse rotations in the same order, you always finish with a non-trivial rotation.
4.7∗ In this problem you derive the wavefunction
hx|pi = eip·x/h̄ (4.4)
of a state of well defined momentum from the properties of the translation operator U (a). The state
|ki is one of well-defined momentum h̄k. How would you characterise the state |k′ i ≡ U (a)|ki? Show
that the wavefunctions of these states are related by uk′ (x) = e−ia·k uk (x) and uk′ (x) = uk (x − a).
Hence obtain equation (4.4).
Soln: U (a)|ki is the result of translating a state of well-defined momentum by k. Moving to the
position representation
uk′ (x) = hx|U (a)|ki = hk|U † (a)|xi∗ = hk|x − ai∗ = uk (x − a)
Also
hx|U (a)|ki = hx|e−ia·p/h̄ |ki = e−ia·k hx|ki = e−ia·k uk (x)
Putting these results together we have uk (x − a) = e−ia·kuk (x). Setting a = x we find uk (x) =
eik·x uk (0), as required.
4.8 By expanding the anticommutator on the left and then applying the third rule of the set (2.22),
show that any three operators satisfy the identity
[{A, B}, C] = {A, [B, C]} + {[A, C], B}. (4.5)
Soln:
[{A, B}, C] = [AB, C] + [BA, C] = [A, C]B + A[B, C] + [B, C]A + B[A, C]
= {A, [B, C]} + {[A, C], B}.
4.9 Let P be the parity operator and S an arbitrary scalar operator. Explain why P and S must
commute.
Soln: A scalar is something that is unaffected by orthogonal transformations of three-dimensional
space. These transformations are made up of the rotations and reflections. The parity operator
reflects states, so it is a member of the group of transformations under which a scalar is invariant.
Consequently S commutes with P .
4.10 In this problem we consider discrete transformations other than that associated with parity.
Let S be a linear transformation on ordinary three-dimensional space that effects a a reflection in a
plane. Let S be the associated operator on kets. Explain the physical relationship between the kets
|ψi and |ψ ′ i ≡ S|ψi. Explain why we can write
Shψ|x|ψi = hψ|S † xS|ψi. (4.6)
What are the possible eigenvalues of S?
13
Given that S reflects in the plane through the origin with unit normal n̂, show, by means of a
diagram or otherwise, that its matrix is given by
Sij = δij − 2ni nj . (4.7)
√
Determine the form of this matrix in the case that n = (1, −1, 0)/ 2. Show that in this case
Sx = yS and give an alternative expression for Sy.
Show that a potential of the form
p
V (x) = f (R) + λxy, where R ≡ x2 + y 2 (4.8)
satisfies V (Sx) = V (x) and explain the geometrical significance of this equation. Show that [S, V ] =
0. Given that E is an eigenvalue of H = p2 /2m + V that has a unique eigenket |Ei, what equation
does |Ei satisfy in addition to H|Ei = E|Ei?
Soln: |ψ ′ i ≡ S|ψi is the state in which the amplitude to be at the reflected point x′ = Sx is equal
to the amplitude to be at x. The mathematical statement of this fact is
hx′ |ψ ′ i = hSx|Ŝ|ψi = hx|ψi, (4.9)
Equation (4.6) states that the expectation value of x for the reflected state is the reflection of the
expectation value of x for the original state. Since x = SSx for all x, S 2 = I and the eigenvalues of
S have to be ±1.
S reverses the component of its target vector that is perpendicular to the plane. Thus we get
x′ by taking from x twice (n̂ · x)n̂:
X X X X
x′ = Sx = x − 2(n̂ · x)n → Sij xj = δij xj − 2( nj xj )ni = (δij − 2ni nj )xj . (4.10)
j j j j
√
For the concrete example given, n̂ = (1, −1, 0)/ 2, so
1 0 0 1 −1 0 0 1 0
Sij = 0 1 0 − −1 1 0 = 1 0 0 (4.11)
0 0 1 0 0 0 0 0 0
Since S(x, y, z) = (y, x, z), we now have
hx, y, z|Sy|ψi = hy, x, z|y|ψi = xhy, x, z|ψi, (4.12)
where the last equality follows because y, being Hermitian, can be considered to act backwards. But
xhy, x, z|ψi = xhx, y, z|S|ψi = hx, y, z|xS|ψi.
Consequently, Sy = xS. Similarly, Sx = yS.
Swapping the x and y coordinates of x clearly leaves invariant both R and xy and therefore
V (Sx) = V (x). Geometrically a contour plot of V is symmetrical around the line x = y, so the
potentials at mirrored points are identical.
hx, y, z|SV |ψi = hy, x, z|V |ψi = V (y, x, z)hy, x, z|ψi
Similarly
hx, y, z|V S|ψi = V (x, y, z)hx, y, z|S|ψi = V (x, y, z)hy, x, z|ψi.
But V (x, y, z) = V (y, x, z) so for and |ψi these two expressions are equal, and it follows that
[V, S] = 0.
14
Figure 5.0 The real part of the wavefunction when a free particle of energy E is scattered by a classically forbidden
square barrier barrier (top) and a potential well (bottom). The upper panel is for a barrier of height V0 = E/0.7 and
half width a such that 2mEa2 /h̄2 = 1. The lower panel is for a well of depth V0 = E/0.2 and half width a such that
2mEa2 /h̄2 = 9. In both panels (2mE/h̄2 )1/2 = 40.
5.7 Reproduce the plots shown in Figure 5.0 of the wavefunctions of particles that are scattered
by a square barrier and a square potential well. Give physical interpretations of as many features
of the plots as you can.
Soln: For the top plot one has to plot the real part of
B sin(−kx + φ) − B ′ sin(−kx + φ′ ) for x < −a
ψ(x) = cos(Kx) + A sin(Kx) for |x| < a
B sin(kx + φ) + B ′ sin(kx + φ′ ) for x > a
p
where k = 40, a = 1/k, K = k 1/0.7 − 1, φ = arctan((k/K) coth(Ka))−ka, φ′ = arctan((k/K) tanh(Ka))−
′
ka, B = cosh(Ka)/ sin(ka+φ), B ′ = −Bei(φ −φ) andp A = B ′ sin(ka+φ′ )/ sinh(Ka). The correspond-
ing equations for the lower panel are a = 3/k, K = k 1 + 1/0.2, φ = − arctan((k/K) cot(Ka))−ka,
′
φ′ = arctan((k/K) tan(Ka)) − ka, B = cos(Ka)/ sin(ka + φ), B ′ = −Bei(φ −φ) and A = B ′ sin(ka +
′
φ )/ sin(Ka).
The dotted lines mark the discontinuities in the potential energy. In the barrier we see the
exponential decay of the wavefunction through the barrier, which flows into rather a small-amplitude
transmitted wave on the right. The profile is close to an exponential because the wave reflected by
the right-hand discontinuity is of very low amplitude compared to that transmitted by the left-
hand discontinuity. This results in A ≃ 1. In the lower panel the shorter wavelength in the well is
apparent, indicating that he particle speeds up as it enters the well. Its high speed is also responsible
to the wave having lower amplitude in the well than on the right, although the amplitude reduction
is modest because in the well particles move in both directions and the flux carried away on the
right is only the difference in the fluxes carried each way in the well. In both panels the wave on the
left represents particles moving in both senses, but this fact is not apparent from the real part of ψ
alone.
15
5.12 An electron moves along an infinite chain of potential wells. For sufficiently low energies
we can assume that the set {|ni} is complete, where |ni is the state of definitely being in the nth
well. By analogy with our analysis of the NH3 molecule we assume that for all n the only non-
vanishing matrix elements of the Hamiltonian are E ≡ hn|H|ni and A ≡ hn ± 1|H|ni. Give physical
interpretations of the numbers A and E.
Explain why we can write
X∞
H= E|nihn| + A (|nihn + 1| + |n + 1ihn|) . (5.1)
n=−∞
P
Writing an energy eigenket |Ei = n an |ni show that
am (E − E) − A (am+1 + am−1 ) = 0. (5.2)
Obtain solutions of these equations in which am ∝ eikm and thus find the corresponding energies
Ek . Why is there an upper limit on the values of k that need be considered?
Initially the electron is in the state
1
|ψi = √ (|Ek i + |Ek+∆ i) , (5.3)
2
where 0 < k ≪ 1 and 0 < ∆ ≪ k. Describe the electron’s subsequent motion in as much detail as
you can.
Soln: E is the energy of an electron in a well. A is the amplitude for tunnelling through to a
neighbouring well. Since the set {|ni} may be considered complete, any operator can be written
P
mn Qmn |mihn|, where Qmn = hm|Q|ni are the operator’s matrix elements. In the case of H most
of the matrix elements vanish and we are left with the given expression. When |Ei is expanded in
the |ni basis, the tise is
X ∞ ∞
X X
H|Ei = E|Ei = E an |ni = {E|nihn| + A (|nihn + 1| + |n + 1ihn|)} am |mi
n=−∞ n=−∞ m
X∞
= Ean |ni + A (an+1 |ni + an |n + 1i)
n=−∞
The required recurrence relation is obtained when this eqn is multiplied on the left by hm|. Substi-
tuting in the trial solution am = Keikm gives
0 = eikm (Ek − E) − A eik(m+1) + eik(m−1)
Dividing through by eikm we find
0 = (Ek − E) − 2A cos(k)
ikm
P
Thus for any k |ki = K m e |mi is an energy eigenstate with Ek = E + 2A cos(k). However,
replacing k with k ′ = k + 2rπ changes neither the energy nor the amplitude to be in any well, so
it does not provide a physically distinct state. Hence we consider only states with k in the range
(0, 2π).
Motion is going to be generated by quantum interference between the states |Ek i and |Ek+∆ i.
The energy difference δE can be estimated by differentiating our expression for Ek :
δE ≃ −2A sin(k)∆
At general time t the state |ψi is
e−iEk t/h̄ X ikm e−iEk t/h̄ X ikm
|ψi = √ e 1 + ei∆m−iδEt/h̄ |mi = √ e 1 + ei∆(m+2A sin(k)t/h̄) |mi
2 m
2 m
16
At some points in the chain, the exponential on the extreme right will cancel on the 1 and the
amplitude to be there will vanish. At other points the exponential will be unity so there is a large
amplitude to be at these wells. As t increases the points of maximum and minimum amplitude will
migrate down the chain as the peaks and troughs of a wave. The speed of this wave is 2A sin(k)/h̄
wells per unit time.
5.13∗ In this problem you investigate the interaction of ammonia molecules with electromagnetic
waves in an ammonia maser. Let |+i be the state in which the N atom lies above the plane of the H
atoms and |−i be the state in which the N lies below the plane. Then when there is an oscillating
electric field E cos ωt directed perpendicular to the plane of the hydrogen atoms, the Hamiltonian in
the |±i basis becomes
E + qEs cos ωt −A
H= . (5.4)
−A E − qEs cos ωt
Transform this Hamiltonian from the |±i basis
√ to the basis provided by the states of well-defined
parity |ei and |oi (where |ei = (|+i + |−i)/ 2, etc). Writing
|ψi = ae (t)e−iEe t/h̄ |ei + ao (t)e−iEo t/h̄ |oi, (5.5)
show that the equations of motion of the expansion coefficients are
dae
= −iΩao (t) ei(ω−ω0 )t + e−i(ω+ω0 )t
dt (5.6)
dao
= −iΩae (t) ei(ω+ω0 )t + e−i(ω−ω0 )t ,
dt
where Ω ≡ qEs/2h̄ and ω0 = (Eo − Ee )/h̄. Explain why in the case of a maser the exponentials
involving ω + ω0 a can be neglected so the equations of motion become
dae dao
= −iΩao (t)ei(ω−ω0 )t ; = −iΩae (t)e−i(ω−ω0 )t . (5.7)
dt dt
Solve the equations by multiplying the first equation by e−i(ω−ω0 )t and differentiating the result.
Explain how the solution describes the decay of a population of molecules that are initially all in
the higher energy level. Compare your solution to the result of setting ω = ω0 in (5.7).
Soln: We have
he|H|ei = 21 (h+| + h−|) H (|+i + |−i)
1
= 2 (h+|H|+i + h−|H|−i + h−|H|+i + h+|H|−i)
= E − A = Ee
1
ho|H|oi = 2 (h+| − h−|) H (|+i − |−i)
1
= 2 (h+|H|+i + h−|H|−i − h−|H|+i − h+|H|−i)
= E + A = Eo
1
ho|H|ei = he|H|oi = 2 (h+| + h−|) H (|+i − |−i)
1
= (h+|H|+i − h−|H|−i + h−|H|+i − h+|H|−i)
2
= qEs cos(ωt)
Now we use the tdse to calculate the evolution of |ψi = ae e−iEe t/h̄ |ei + ao e−iEo t/h̄ |oi:
∂|ψi
ih̄ = ih̄ȧe e−iEe t/h̄ |ei+ae Ee e−iEe t/h̄ |ei+ih̄ȧo e−iEo t/h̄ |oi+ao Eo e−iEo t/h̄ |oi = ae e−iEe t/h̄ H|ei+ao e−iEo t/h̄ H|oi
∂t
We now multiply through by first he| and then ho|. After dividing through by some exponential
factors to simplify, we get
ih̄ȧe + ae Ee = ae he|H|ei + ao ei(Ee −Eo )t/h̄ he|H|oi
ih̄ȧo + ao Eo = ae ei(Eo −Ee )t/h̄ ho|H|ei + ao ho|H|oi
17
From the condition ae (0) = 0, the required solution is ae (t) ∝ (eiω+ t − eiω− t ) and the constant of
proportionality follows from the second initial condition, so finally
−Ω
ae (t) = p (eiω+ t − eiω− t ) (∗)
(ω − ω0 )2 + 4Ω2
The probability oscillates between the odd and even states. First the oscillating field stimulates
emission of radiation and decay from |oi to |ei. Later the field excites molecules in the ground state
to move back up to the first-excited state |oi.
If we solve the original equations (1) exactly on resonance (ω = ω0 ), the relevant solution is
ae (t) = 12 (e−iΩt − eiΩt ),
which is what our general solution (∗) reduces to as ω → ω0 .
5.14 238 U decays by α emission with a mean lifetime of 6.4 Gyr. Take the nucleus to have a
diameter ∼ 10−14 m and suppose that the α particle has been bouncing around within it at speed
∼ c/3. Modelling the potential barrier that confines the α particle to be a square one of height V0
and width 2a, give an order-of-magnitude estimate of W = (2mV0 a2 /h̄2 )1/2 . Given that the energy
released by the decay is ∼ 4 MeV and the atomic number of uranium is Z = 92, estimate the width
of the barrier through which the α particle has to tunnel. Hence give a very rough estimate of the
barrier’s typical height. Outline numerical work that would lead to an improved estimate of the
structure of the barrier.
Soln: The time between impacts with the wall is ∼ 10−14 /108 s, so in 6.4 Gyr ≃ 2 × 1017 s the α
makes ∼ 2 × 1039 escape attempts and the probability of any one attempt succeeding is ∼ 5 × 10−40 .
This probability is sin2 (∆φ) ≃ (4k/K)2 e−4W . Taking 4k/K to be order unity, we conclude that
W = 14 ln(1039 ) = 22.5.
The distance 2a from the nucleus at which the α becomes classically allowed is that at which
its energy of electrostatic repulsion from the thorium nucleus (Z = 92 − 2 = 90) is 4 MeV: 2a =
2 × 90e2 /(4πǫ0 4 × 106 e) = 6.5 × 10−14 m. If the barrier were square rather than wedge shaped, it
would not be so thick, so let’s adopt 2a = 3×10−14 m. Then V0 = (22.5h̄/a)2 /8mp ≃ 1.6×10−12 J ≃
10 MeV.
This calculation has two major problems: (a) the assumption of a square barrier instead of
a skew and sharply peaked one, and (b) the use of a one-dimensional rather than spherical model
given that α becomes allowed at a radius that is significantly greater than the radius of the nucleus.
Problem (a) is probably the most serious and could be addressed by numerically solving the one-d
tise as in Problem 5.15. To address problem (b) we would solve the radial wavefunction for a
particle with vanishing angular momentum.
18
5.15∗ Particles of mass m and momentum h̄k at x < −a move in the potential
0 for x < −a
(
V (x) = V0 12 [1 + sin(πx/2a)] for |x| < a (5.8)
1 for x > a,
where V0 < h̄2 k 2 /2m. Numerically reproduce the reflection probabilities plotted Figure 5.20 as
follows. Let ψi ≡ ψ(xj ) be the value of the wavefunction at xj = j∆, where ∆ is a small increment
in the x coordinate. From the tise show that
ψj ≃ (2 − ∆2 k 2 )ψj+1 − ψj+2 , (5.9)
p
where k ≡ 2m(E − V )/h̄. Determine ψj at the two grid points with the largest values of x from
a suitable boundary condition, and use the recurrence relation (5.9) to determine ψj at all other
grid points. By matching the values of ψ at the points with the smallest values of x to a sum of
sinusoidal waves, determine the probabilities required for the figure. Be sure to check the accuracy
of your code when V0 = 0, and in the general case explicitly check that your results are consistent
with equal fluxes of particles towards and away from the origin.
Equation (11.40) gives an analytical approximation for ψ in the case that there is negligible
reflection. Compute this approximate form of ψ and compare it with your numerical results for
larger values of a.
Soln:
We discretise the tise
h̄2 d2 ψ h̄2 ψj+1 + ψj−1 − 2ψj
− + V ψ = Eψ by − + Vj ψj = Eψj
2m dx2 2m ∆2
which readily yields the required recurrence relation. At the right-hand boundary we require a pure
outgoing wave, so ψj = exp(ijK∆) gives ψ at the two last grid points. From the recurrence relation
we obtain ψ elsewhere. At the left boundary we solve for A+ and A− the equations
A+ exp(i0k∆) + A− exp(−i0k∆) = ψ0
A+ exp(i1k∆) + A− exp(−i1k∆) = ψ1
The transmission probability is (K/k)/|A+ |2 . The code must reproduce the result of Problem 5.4
in the appropriate limit.
5.16∗ In this problem we obtain an analytic estimate of the energy difference between the even-
and odd-parity states of a double square well. Show that for large θ, coth θ − tanh θ ≃ 4e−2θ . Next
letting δk be the difference between the k values that solve
s p
2 coth W 2 − (ka)2 even parity
W
tan [rπ − k(b − a)] 2
−1= p (5.10a)
(ka) tanh W 2 − (ka)2 odd parity,
where r
2mV0 a2
W ≡ (5.10b)
h̄2
19
In the case of interest the right side of the original equation is close to unity, so we can simplify the
last equation by using s
W2
tan [rπ − k(b − a)] −1≃1
(ka)2
With the help of the identity s2 θ = 1 + tan2 θ we obtain the required relation. We now approximate
the left side for W ≫ ka. This yields
W √ 2
(b − a)δk ≃ −4e−2W 1−(ka/W ) ($)
ka
Since E = h̄2 k 2 /2m, δE/E = 2δk/k and
2mEa2 h̄2
(ka/W )2 = × = E/V0 .
h̄2 2mV0 a2
The required relation follows when we use these relations in ($).
6.3 Given that the state |ABi of a compound system can be Pwritten as a product |Ai|Bi of states
of the individual systems, show that when |ABi is written as ij cij |A; ii|B; ji in terms of arbitrary
basis vectors for the
P subsystems, every column of the matrix cij is a multiple of the leftmost column.
Soln: If |Ai = i ai |A; ii and similarly for |Bi, cij = ai bj , so the ratio of corresponding terms in
the j th and k th column is cij /cik = bj /bk is independent of i as required.
6.6 Show that when the Hadamard operator UH is applied to every qubit of an n-qubit register
that is initially in a member |mi of the computational basis, the resulting state is
n
2X −1
1
|ψi = ax |xi, (6.1)
2n/2 x=0
where ax = 1 for all x if m = 0, but exactly half the ax = 1 and the other half the ax = −1 for any
other choice of m. Hence show that
1 X |0i if all ax = 1
UH ax |xi = (6.2)
2 n/2 |mi 6 = |0i if half the ax = 1 and the other ax = −1.
x
20
Soln: We are given that the register is in a member of the computational basis, so each of its
th
qubits is in either√the state |0i or the state |1i. The Hadamard √ operator puts the r qubit into the
state (|0i + |1i)/ 2 if it was in the state |0i and (|0i − |1i)/ 2 if it was in the state |1i. So if |mi =
1
|0i = |0i × · · · |0i the Hadamard operators produce U n |0i = 2n/2 (|0i + |1i)(|0i + |1i) · · · (|0i + |1i).
When we multiply this up we get every possible sequence of |0i and |1i exactly once and all added
together.
Each original qubit that is |1i gives rise to a factor (|0i − |1i). Multiply all the factors (|0i + |1i)
together to get all positive terms. Then multiply in the first factor (|0i − |1i) and half the resulting
terms will be positive and half negative from the factor −|1i. When we multiply by the next factor
(|0i − |1i) half of the existing terms will change sign and half will retain their sign. Hence we will
still have half positive and half negative terms. The required result now follows by induction.
6.7 Show that the trace of every Hermitian operator is real.
Soln: Since the trace is invariant under unitary transformations, it is equal to the sum of the
operator’s eigenvalues. But these are all real.
6.8 Let ρ be the density operator of a two-state system. Explain why ρ can be assumed to have
the matrix representation
a c
ρ= , (6.3)
c∗ b
where a and b are real numbers. Let E0 and E1 > E0 be the eigenenergies of this system and |0i and
|1i the corresponding stationary states. Show from the equation of motion of ρ that in the energy
representation a and b are time-independent while c(t) = c(0)eiωt with ω = (E1 − E0 )/h̄.
Determine the √ values of a, b and c(t) for the case that initially the system is in the state
|ψi = (|0i + |1i)/ 2. Given that the parities of |0i and |1i are even and odd respectively, find the
time evolution of the expectation value x in terms of the matrix element h0|x|1i. Interpret your
result physically.
Soln: The density operator is Hermitian and any 2 × 2 Hermitian matrix can be written thus.
With this form of ρ and in the E representation, the equation of motion is
ȧ ċ E0 0 a c a c E0 0 0 −c
ih̄ ∗ = − = (E1 − E0 )
ċ ḃ 0 E1 c∗ b c∗ b 0 E1 c∗ 0
So a and b are constant and c satisfies ċ = iωc, and the required form of c(t) follows trivially.
For the given initial condition
1 1 1 1
ρ(0) = 2 (|0i + |1i)(h0| + h1|) ⇒ ρ(0) = 2
1 1
We get ρ(t) from our solutions to the equations of motion of a, b and c. From the given parity
condition we know that the diagonal elements of the matrix for x vanish, so
eiωt
1 0 x01
1
x = 2 Tr = x01 eiωt + x∗01 e−iωt = 2|x01 | cos(ωt + φ01 ),
e−iωt 1 x∗01 0
where φ01 = arg(x01 ). We have recovered a result familiar from our study of the harmonic oscillator:
oscillations occur at the frequency associated with the energy difference between adjacent stationary
states as a result of interference between the amplitudes of these states.
6.9 In this problem we consider an alternative interpretation of the density operator. Any quantum
state can be expanded in the energy basis as
N
X √
|ψ; φi ≡ pn eiφn |ni, (6.4)
n=1
where φn is real and pn is the probability that a measurement of energy will return En . Suppose
we know the values of the pn but not the values of the phases φn . Then the density operator is
Z 2π N
d φ
ρ= |ψ; φihψ; φ|. (6.5)
0 (2π)N
21
P
Show that this expression reduces to n pn |nihn|. Contrast the physical assumptions made in this
derivation of ρ with those made in §6.3.
Clearly |ψ; φi can be expanded in some other basis {|qr i} as
Xp
|ψ; φi ≡ Pr eiηr |qr i, (6.6)
r
where Pr is the probability of obtaining qr on a measurement of the observable Q and the ηr (φ)
are unknown phases. Why does this second expansion not lead to the erroneous conclusion that ρ
is necessarily diagonal in the {|qr i} representation?
Soln: In §6.3 we assumed that the system was in an (unknown) eigenstate |ni, which is in practice
usually taken to be a stationary state. This is an improbable assumption. Here we assume that
the system is in a generic and non-stationary state. So this physical interpretation of the density
operator is superior to our earlier one.
In equation (6.5) we assume that all values of the phases φn are equally likely. Assuming that
all values of the phases ηr are equally likely will lead to ρ being diagonal in the {|qr i} representaton.
By substituting X
|ni = |hqr |ni| eiχnr |qr i
r
into (6.4), we can show that
!
1 X r pn
ηr (φ) = ln |hqr |ni| ei(φn +χnr ) .
i n
Pr
Thus the dependence of ηr on the φn is highly non-linear, and a uniform distribution in the φn will
in general not correspond to a uniform distribution in the ηr .
6.10∗ Show that when the density operator takes the form ρ = |ψihψ|, the expression Q = Tr Qρ for
the expectation value of an observable can be reduced to hψ|Q|ψi. Explain the physical significance
of this result. For the given form of the density operator, show that the equation of motion of ρ
yields
∂|ψi
|φihψ| = |ψihφ| where |φi ≡ ih̄ − H|ψi. (6.7)
∂t
Show from this equation that |φi = a|ψi, where a is real. Hence determine the time evolution of |ψi
given the at t = 0 |ψi = |Ei is an eigenket of H. Explain why ρ does not depend on the phase of
|ψi and relate this fact to the presence of a in your solution for |ψ, ti.
Soln: X
Tr(Qρ) = hn|Q|ψihψ|ni
n
We choose a basis that |ψi is a member. Then there is only one non-vanishing term in the sum,
when |ni = |ψi, and the right side reduces to hψ|Q|ψi as required. This result shows that density
operators recover standard experimental predictions when the system is in a pure state.
Differentiating the given ρ we have
dρ ∂|ψi ∂hψ| 1
= hψ| + |ψi = (H|ψihψ| − |ψihψ|H)
dt ∂t ∂t ih̄
Gathering the terms proportional to hψ| on the left and those proportional to |ψi on the right we
obtain the required expression. Now
|φihψ| = |ψihφ| ⇒ |φihψ|φi = |ψihφ|φi,
which establishes that |φi ∝ |ψi. We define a as the constant of proportionality. Using |φi = a|ψi
in |φihψ| = |ψihφ| we learn that a = a∗ so a is real.
Returning to the definition of |φi we now have
∂|ψi
ih̄ = (H − a)|ψi.
∂t
22
This differs from the tdse in having the term in a. If |ψi is an eigenfunction of H, we find that its
time dependence is |ψ, ti = |ψ, 0ie−i(E−a)t/h̄ rather than the expected result |ψ, ti = |ψ, 0ie−iEt/h̄ .
We cannot determine a from the density-matrix formalism because ρ is invariant under the trans-
formation |ψi → e−iχ |ψi, where χ is any real number.
P
6.11 The density operator is defined to be ρ = α pα |αihα|, where pα is the probability
P that the
system is in the state α. Given an arbitrary basis {|ii} and the expansions |αi = i aαi |ii, calculate
the matrix elements ρij = hi|ρ|ji of ρ. Show that the diagonal elements ρii are non-negative real
numbers and interpret them as probabilities.
Soln: X X
ρ= pα aαi a∗αj ρij |iihj| ⇒ ρij = pα aαi a∗αj
αij α
The diagonal element ρii = α pα |aαi |2 is the product of two non-negative real numbers so it’s real
P
non-negative itself.
Consider the observable Q = |iihi| which is diagonal in the given basis and has eigenvalues 1
and 0: you measure unity if the system is in the state |ii and zero otherwise. ρii = Tr(Qρ) = Q is
the expectation value of Q. Consequently it’s the probability that if you measure Q you force the
system into the state |ii. Loosely speaking we can consider ρii to be the probability of finding the
system in the state |ii.
P
6.12 Consider the density operator ρ = ij ρij |iihj| of a system that is in a pure state. Show
that every row of the matrix ρij is a multiple of the first row and every column is a multiple of the
first column. Given that these relations between the rows and columns of a density matrix hold,
show that the system is in a pure state. Hint: exploit the real, non-negativity of ρ11 established in
Problem 6.11 and the Hermiticity of ρ.
Soln: Since the system is in a pure state |ψi = i ai |ii, ρ = |ψihψ|. So ρij = ai a∗j . The ratio of
P
The result follows on substituting for hn|[Q, H]|ni in the previous equation.
6.14 Find the probability distribution (p1 , . . . , pn ) for n possible outcomes that maximises the
Shannon entropy. Hint: use a Lagrange multiplier.
Soln: Bearing in mind the requirement for the pi to be normalised, we set to zero each derivative
∂/∂pj of !
X X
S= pi ln pi + λ pi − 1
i i
and obtain
0 = ln pj − 1 + λ ⇒ pj = e1−λ .
Thus the pi are all equal, and from the normalisation condition they must equal 1/n.
6.15 Use Lagrange multipliers λ and β to P extremise the Shannon
P entropy of the probability dis-
tribution {pi } subject to the constraints (i) i pi = 1 and (ii) i pi Ei = U . Explain the physical
significance of your result.
Soln: We set to zero each derivative ∂/∂pj of
! !
X X X
S= pi ln pi + λ pi − 1 + β pi Ei − U
i i i
and obtain
0 = ln pj − 1 + λ + βEj ⇒ pj = e1−λ e−βEj .
We determine λ from the condition i pi = 1, so eλ−1 is the partition function Z. We determine
P
β from the condition that the expectation of the energy is U : β is the inverse temperature 1/kB T
because the smaller it is the higher U is. This calculation establishes that the Gibbs distribution
maximises the Shannon entropy of the system subject to the given constraints.
6.16 A composite system is formed from uncorrelated subsystem A and subsystem B, both in
impure states. The numbers {pAi } are the probabilities of the members of the complete set of states
{|A; ii} for subsystem A, while the numbers {pBi } are the probabilities of the complete set of states
{|B; ii} for subsystem B. Show that the Shannon entropy of the composite system is the sum of the
Shannon entropies of its subsystems. What is the relevance of this result for thermodynamics?
Soln: The state |AB; iji of the composite system has probability pij = pAi pBj so the Shannon
entropy of the whole system is
X X X X X
s= pij ln(pij ) = pAi pBj ln(pAi ) + pAi pBj ln(pBj ) = pAi ln(pAi ) + pBj ln(pBj )
ij ij ij i j
= sA + sB
P
where we’ve used the normalisation conditions i pAi = 1, etc. This result establishes that entropy
is an “extensive” thermodynamic quantity: the entropy of two litres of water at 20◦ C is twice the
entropy of two one-litre bottles of water at the same temperature.
6.17 The |0i state of a qubit has energy 0, while the |1i state has energy ǫ. Show that when the
qubit is in thermodynamic equilibrium at temperature T = 1/(kB β) the internal energy of the qubit
is ǫ
U = βǫ . (6.10)
e +1
Show that when βǫ ≪ 1, U ≃ 21 ǫ, while for βǫ ≫ 1, U ≃ ǫe−βǫ . Interpret these results physically
and sketch the specific heat C = ∂U/∂T as a function of T .
Soln: The partition function is
Z = 1 + e−βǫ
so the internal energy is
∂ ln Z ǫe−βǫ ǫ
U =− = −βǫ
= βǫ .
∂β 1+e e +1
When βǫ ≪ 1 we can approximate eβǫ ≃ 1 + βǫ and then binomial-expand the denominator in U :
ǫ
U≃ ≃ 21 ǫ(1 − βǫ + · · ·)
2 + βǫ
24
This argument and the result it leads to is wrong. Obtain the correct result by integrating by parts
dθ sin θ dφ (f ∗ L+ g), where f and g are arbitrary functions of θ and φ. What is the fallacy in
R R
the given argument?
Soln: Z
∂g ∂g
Z Z Z
dθ sin θ dφ (f ∗ L+ g) = dθ sin θ dφ f ∗ eiφ + i cot θ
∂θ ∂φ
∂g ∂g
Z Z Z Z
= dφ eiφ dθ sin θf ∗ + i dθ cos θ dφ f ∗ eiφ
∂θ ∂φ
∗
∂(sin θf )
Z Z
= dφ eiφ [sin θ f ∗ g] − dθ g
∂θ
∂(f ∗ eiφ )
Z Z
+ i dθ cos θ [f ∗ eiφ g] − dφ g
∂φ
The square brackets vanish so long f, g are periodic in φ. Differentiating out the products we get
∂f ∗
Z Z Z Z Z
dθ sin θ dφ (f ∗ L+ g) = − dφ eiφ dθ sin θg + dθ cos θgf ∗
∂θ
∂f ∗
Z Z Z
− i dθ cos θ dφ eiφ g + i dφ eiφ gf ∗
∂φ
The two integrals containing f ∗ g cancel as required leaving us with
∗
∂f ∗
Z
∂f
Z Z Z Z Z
dθ sin θ dφ (f ∗ L+ g) = − dθ sin θ dφ geiφ + i cot θ = dθ sin θ dφ g(L− f )∗
∂θ ∂φ
where ∂ ∂
L− = −e−iφ − i cot θ .
∂θ ∂φ
The fallacy is the proposition that ∂/∂θ is anti-Hermitian: the inclusion of the factor sin θ in the
integral prevents this being so.
7.4∗ By writing h̄2 L2 = (x × p) · (x × p) = ijklm ǫijk xj pk ǫilm xl pm show that
P
h̄2 L2 1
p2 = + 2 (r · p)2 − ih̄r · p .
2
(7.2)
r r
By showing that p · r̂ − r̂ · p = −2ih̄/r, obtain r · p = rpr + ih̄. Hence obtain
h̄2 L2
p2 = p2r +
. (7.3)
r2
Give a physical interpretation of one over 2m times this equation.
Soln: From the formula for the product of two epsilon symbols we have
X
h̄2 L2 = (δjl δkm − δjm δkl )xj pk xl pm
jklm
X
= xj pk xj pk − xj pk xk pj .
jk
The first term is X X X
xj pk xj pk = xj (xj pk + [pk , xj ])pk = xj (xj pk − ih̄δjk )pk
jk jk jk
= r2 p2 − ih̄r · p.
27
7.9 Consider a stationary state |E, li of a free particle of mass m that has angular-momentum
quantum number l. Show that Hl |E, li = E|E, li, where
l(l + 1)h̄2
1 2
Hl ≡ pr + . (7.4)
2m r2
Give a physical interpretation of the two terms in the big bracket. Show that Hl = A†l Al , where
1 (l + 1)h̄
Al ≡ √ ipr − . (7.5)
2m r
Show that [Al , A†l ] = Hl+1 − Hl . What is the state Al |E, li? Show that for E > 0 there is no upper
bound on the angular momentum. Interpret this result physically.
Soln: The Hamiltonian is
p2 h̄2 L2
1
H= = p2r + 2
2m 2m r
For kets |E, li that are simultaneous eigenkets of l, this expression immediately reduces to the
required expression. The second term in the bracket is the kinetic energy of tangential motion
28
(classically 12 mvt2 ) and the first term is the radial kinetic energy (classically 12 mvr2 ). For the given
form of Al
1 (l + 1)h̄ (l + 1)h̄
A†l Al = −ipr − ipr −
2m r r
2 2
1 (l + 1) h̄
= p2r + + i(l + 1)h̄[p r , r −1
]
2m r2
But [pr , r−1 ] = −r−2 [pr , r] = ih̄r−2 so
(l + 1)2 h̄2 (l + 1)h̄2
1
A†l Al = p2r + − = Hl
2m r2 r2
as required.
Clearly,
(l + 1)h̄2
† 1 (l + 1)h̄ (l + 1)h̄ (l + 1)ih̄
[Al , Al ] = ipr − , −ipr − =− [pr , r−1 ] =
2m r r m mr2
But Hl+1 − Hl = {(l + 1)(l + 2) − l(l + 1)}h̄2 /(2mr2 ) = (l + 1)h̄2 /mr2 .
Now
Al Hl = Al A†l Al = (A†l Al + [Al , A†l ])Al = (Hl + Hl+1 − Hl )Al = Hl+1 Al
So if we multiply both sides of E|E, li = Hl |E, li by Al , we get
EAl |E, li = Al Hl |E, li = Hl+1 Al |E, li,
which establishes that Al |E, li ∝ |E, l + 1i is a state with the same energy but more angular
momentum. To see whether there is an upper limit on the angular momentum, we evaluate
2
|Al |E, li| = hE, l|A†l Al |E, li = E > 0
so there is no limit to the angular momentum. Physically, when there is no confining potential, with
a given energy the particle can move at a given speed along a path that is as far from the origin as
it pleases. Hence its angular momentum is not constrained by its energy.
7.10 Write down the expression for the commutator [σi , σj ] of two Pauli matrices. Show that the
anticommutator of two Pauli matrices is
{σi , σj } = 2δij . (7.6)
Soln: Since σi = 2si is twice the spin operator of a spin-half particle, from the angular-momentum
commutation relations we have [σi , σj ] = 2iǫijk σk .
Calculating explicitly
2 0 1 0 1 1 0
σx = =
1 0 1 0 0 1
and similarly σy2 = σz2 = I so {σi , σi } = 2I as required. Further
0 1 0 −i 0 −i 0 1 i 0 −i 0
{σx , σy } = + = + =0
1 0 i 0 i 0 1 0 0 −i 0 i
and similarly for {σy , σz } and {σx , σz }.
7.13 Show that a classical top with spin angular momentum S which is subject to a torque G =
µS × B/|S| precesses at angular velocity ω = µB/|S|. Explain the relevance of this calculation to
nuclear magnetic resonance in general and equation (7.137b) in particular.
Soln: Since torque is equal to rate of change of angular momentum, we have
dS
= µS × B/|S|.
dt
29
When a vector S is rotated with angular velocity ω its rate of change is (Box 4.3)
dS
= ω × S.
dt
Comparing these equations we see that the torque rotates S at a rate ω = −µB/|S|.
The relevance to NMR is that classically the magnetic field would apply a torque µp S × B/|S|
to the proton, and for the proton |S| = 21 h̄, so classical physics predicts precession of the proton’s
spin at the angular frequency given by equation (7.137b).
7.14 For a spin-half particle at rest, the rotation operator J is equal to the spin operator S. Use
the result of Problem 7.10 to show that in this case the rotation operator U (α) ≡ exp(−iα · J) is
α α
U (α) = I cos − iα̂ · σ sin , (7.7)
2 2
where α̂ is the unit vector parallel to α. Comment on the value this gives for U (α) when α = 2π.
Soln:
(α · σ/2)2 (α · σ/2)3 (α · σ/2)4
U (α) = exp(− 21 iα · σ) = I − iα · σ/2 − +i + + ···
2! 3! 4!
2n
Using α = αα̂ and (α̂ · σ) = I, and gathering real and imaginary parts, we have
(α/2)2 (α/2)4 (α/2)3
U (α) = 1 − + + · · · I − i (α/2) − + · · · (α̂ · σ).
2! 4! 3!
The required result now follows from the power series for cos θ and sin θ. For α = 2π we have rotated
the system right around but the value of the wavefunction has changed sign because U = −I.
7.18∗ Repeat the analysis of Problem 7.15 for spin-one particles coming on filters aligned succes-
sively along +z, 45◦ from z towards x [i.e. along (1,0,1)], and along x.
Use classical electromagnetic theory to determine the outcome in the case that the spin-one
particles were photons and the filters were polaroid. Why do you get a different answer? Which
answer is correct?
Soln: We√adapt the calculation of Problem 7.15 by replacing the matrix for Jx by that for n · J =
(Jx + Jz )/ 2. So if now (a, b, c) is | + ni in the usual basis, we have
b
−1/2 1
a = 2 − √2
2 2 0 a a
12 0 1
2
b = b ⇒
1 −1/2 b
0 −2 c c c= √
2
2+ 2
√
The normalisation yields b = 12 , so a = 21 /(2 − 2) and the required probability is the square of this,
√
0.25/(6 − 4 2) ≃ 0.73. So the probability of getting through all three filters is 31 × (0.73)2 ≃ 0.177.
In electromagnetism just one of two polarisations gets through the first filter, so we must say
that a photon has a probability of half of passing the first filter. Then we resolve
√ its E field along the
direction of the second filter and find that the amplitude of E falls by 1/ 2 on passing the second
filter, so half the energy and therefore photons that pass the first filter pass the second. Of these
just a half pass the third filter. Hence in total 81 = 0.125 of the photons get right through.
Although photons are spin-one particles, there are two major difference between the two cases.
Most obviously, polaroid selects for linear polarisation rather than circular polarisation, and a photon
with well-defined angular momentum is circularly polarised. The other difference is that a photon
can be in the state | + zi or | − zi but not the state |0zi, where the z-axis is parallel to the photon’s
motion. This fact arises because emag waves are transverse so they do not drive motion in the
direction of propagation k; an angular momentum vector perpendicular to k would require motion
along k. Our theory does not allow for this case because it is non-relativistic, whereas a photon,
having zero rest mass, is an inherently relativistic object; we cannot transform to a frame in which
a photon is at rest so all three directions would be equivalent.
7.19∗ Show that l excitations can be divided amongst the x, y or z oscillators of a three-dimensional
harmonic oscillator in ( 12 l + 1)(l + 1) ways. Verify in the case l = 4 that this agrees with the number
of states of well defined angular momentum and the given energy.
30
Soln: If we assign nx of the l excitations to the x oscillator, we can assign 0, 1, . . . , l−nx excitations
to the y oscillator [(l − nx + 1) possibilities], and the remaining excitations go to z. So the number
of ways is
X l Xl Xl
S≡ (l − nx + 1) = (l + 1) − nx = (l + 1)2 − 21 l(l + 1) = (l + 1)( 21 l + 1)
nx =0 nx =0 nx =1
In the case of 4 excitations, the possible values of l are 4, 2 and 0, so the number of states is
(2 ∗ 4 + 1) + (2 ∗ 2 + 1) + 1 = 15, which is indeed equal to (4 + 1) ∗ (2 + 1).
7.20∗ Let
1 (l + 1)h̄
Al ≡ √ ipr − + mωr . (7.8)
2mh̄ω r
be the ladder operator of the three-dimensional harmonic oscillator and |E, li be the oscillator’s
stationary state of energy E and angular-momentum√ quantum number l. Show that if we write
Al |E, li = α− |E − h̄ω, l + 1i, then α− = L − l, where L is the angular-momentum quantum
number of a circular orbit of energy E. Show similarly that if A†l−1 |E, li = α+ |E + h̄ω, l − 1i, then
√
α+ = L − l + 2.
Soln: Taking the mod-square of each side of Al |E, li = α− |E − h̄ω, l + 1i we find
Hl E
|α− |2 = hE, l|A†l Al |E, li = hE, L| − (l + 32 ) |E, li = − (l + 23 ).
h̄ω h̄ω
3
In the case l = L, |α− |2 = 0, so L = (E/h̄ω) − 2 and therefore |α− |2 = L − l as required. We can
choose the phase of α− at our convenience.
Similarly
α2+ = hE, l|Al−1 A†l−1 |E, li = hE, l|(A†l−1 Al−1 + [Al−1 , A†l−1 ])|E, li
Hl−1 1 Hl − Hl−1 E 1
= hE, l| − (l + 2 ) + + 1 |E, li = −l+ 2 =L−l+2
h̄ω h̄ω h̄ω
7.21∗ Show that the probability distribution in radius of a particle that orbits in the three-
dimensional harmonic-oscillator
p potential on a circular orbit with angular-momentum quantum
number l peaks at r/ℓ = 2(l + 1), where
r
h̄
ℓ≡ . (7.9)
2mω
Derive the corresponding classical result.
Soln: The radial wavefunctions of circular orbits are annihilated by Al , so Al |E, li = 0. In the
position representation this is
∂ 1 l+1 r
+ − + 2 u(r) = 0
∂r r r 2ℓ
Using the integrating factor,
Z
l r
= r−l exp r2 /4ℓ2 ,
exp dr − + 2 (7.10)
r 2ℓ
2 2 2
/2ℓ2
to solve the equation, we have u ∝ rl e−r /4ℓ . The radial distribution is P (r) ∝ r2 |u|2 = r2(l+1) e−r .
Differentiating to find the maximum, we have
√
2(l + 1)r2l+1 − r2(l+1) r/ℓ2 = 0 ⇒ r = 2(l + 1)1/2 a
For the classical result we have
mv 2 lh̄
mrv = lh̄ and = mω 2 r ⇒ r = v/ω =
r mrω
31
(1) (2)
S+ S− = (Sx(1) + iSy(1) )(Sx(2) − iSy(2) ) = Sx(1) Sx(2) + Sy(1) Sy(2) + i(Sy(1) Sx(2) − Sx(1) Sy(2) )
Similarly,
(1) (2)
S− S+ = Sx(1) Sx(2) + Sy(1) Sy(2) − i(Sy(1) Sx(2) − Sx(1) Sy(2) )
Adding these expressions we obtain the desired relation.
We choose to orient the z-axis along a. Then H becomes
K 1 (1) (2) (1) (2) (1) (2) (1) (2)
H= (S S + S S ) + S S − 3S S . (7.13)
a 2 + − − + z z z z
The eigenkets of S 2 and Sz are the three spin-one kets |1, 1i, |1, 0i and |1, −1i and the single spin-zero
ket |0, 0i. We multiply each of these kets in turn by H:
K 1 (1) (2) (1) (2)
H|1, 1i = H|+i|+i = 2 (S+ S− + S− S+ ) − 2Sz(1) Sz(2) |+i|+i
a
K
= − |1, 1i
2a
(i)
which uses the fact that S+ |+i = 0. Similarly H|1, −1i = H|−i|−i = −(K/2a)|1, −1i.
1 K 1 (1) (2) (1) (2) (1) (2)
H|1, 0i = H √ (|+i|−i + |−i|+i) = √ (S + S − + S − S + ) − 2S z S z (|+i|−i + |−i|+i)
2 2a 2
K 1 K
=√ 2 + 1 (|+i|−i + |−i|+i) = a |1, 0i
2a
32
7.25∗ In this problem you show that many matrix elements of the position operator x vanish when
states of well defined l, m are used as basis states.
P These results will lead to selection rules for electric
dipole radiation. First show that [L2 , xi ] = i jk ǫjik (Lj xk + xk Lj ). Then show that L · x = 0 and
using this result derive
X
[L2 , [L2 , xi ]] = i ǫjik Lj [L2 , xk ] + [L2 , xk ]Lj = 2(L2 xi + xi L2 ).
(7.14)
jk
By squeezing this equation between angular-momentum eigenstates hl, m| and |l′ , m′ i show that
0 = (β − β ′ )2 − 2(β + β ′ ) hl, m|xi |l′ , m′ i,
where β ≡ l(l + 1) and β ′ = l′ (l′ + 1). By equating the factor in front of hl, m|xi |l′ , m′ i to zero,
and treating the resulting equation as a quadratic equation for β given β ′ , show that hl, m|xi |l′ , m′ i
must vanish unless l + l′ = 0 or l = l′ ± 1. Explain why the matrix element must also vanish when
l = l′ = 0.
Soln: X X X
[L2j , xi ] = (Lj [Lj , xi ] + [Lj , xi ]Lj ) = i ǫjik (Lj xk + xk Lj )
j j jk
X X X
h̄L · x = ǫijk xj pk xi = ǫijk (xj xi pk + xj [pk , xi ]) = ǫijk (xj xi pk − ih̄xj δki )
ijk ijk ijk
P
Both terms on the right side of this expression involve ik ǫijk Sik where Sik = Ski so they vanish
by Problem 7.1. Hence x · L = 0 as in classical physics.
Now X X
[L , [L2 , xi ]] = i
2
ǫjik [L2 , (Lj xk + xk Lj )] = i ǫjik (Lj [L2 , xk ] + [L2 , xk ]Lj )
jk jk
X
=− ǫjik ǫlkm (Lj {Ll xm + xm Ll } + {Ll xm + xm Ll }Lj )
jklm
X
=− (δjm δil − δjl δim )(Lj {Ll xm + xm Ll } + {Ll xm + xm Ll }Lj )
jlm
X
=− (Lj {Li xj + xj Li } + {Li xj + xj Li }Lj − Lj {Lj xi + xi Lj } − {Lj xi + xi Lj }Lj )
j
X X
=− (Lj Li xj + xj Li Lj ) − L2 xi − (Lj xi Lj + Lj xi Lj ) − xi L2
j j
33
where to obtain the last line we have identified occurrences of L · x and x · L. Now
X X X
Lj Li xj = (Lj xj Li + Lj [Li , xj ]) = i ǫijk Lj xk
j j jk
P P
Similarly, j xj Li Lj = i jk ǫjik xk Lj . Moreover
X X X
Lj xi Lj = ([Lj , xi ]Lj + xi Lj Lj ) = i ǫjik xk Lj + xi L2
j j jk
X X
= (Lj [xi , Lj ] + Lj Lj xi ) = i ǫijk Lj xk + L2 xi
j jk
Assembling these results we find
X X
[L2 , [L2 , xi ]] = − i ǫijk [Lj , xk ] − L2 xi − i ǫjik [xk , Lj ] − xi L2 − L2 xi − xi L2
jk jk
2 2
= 2(L xi + xi L )
as required. The relevant matrix element is
hlm|[L2 , [L2 , xi ]]|l′ m′ i = hlm|(L2 L2 xi − 2L2 xi L2 + xi L2 L2 )|l′ m′ i = 2hlm|(L2 xi + xi L2 )|l′ m′ i
which implies
β 2 hlm|xi |l′ m′ i − 2βhlm|xi |l′ m′ iβ ′ + hlm|xi |l′ m′ iβ ′2 = 2βhlm|xi |l′ m′ i + 2hlm|xi |l′ m′ iβ ′
Taking out the common factor we obtain the required result.
The quadratic for β(β ′ ) is
β 2 − 2(β ′ + 1)β + β ′ (β ′ − 2) = 0
so p p
β = β ′ + 1 ± (β ′ + 1)2 − β ′ (β ′ − 2) = β ′ + 1 ± 4β ′ + 1
p
= l′ (l′ + 1) + 1 ± 4l′2 + 4l′ + 1 = l′ (l′ + 1) + 1 ± (2l′ + 1)
= l′2 + 3l′ + 2 or l′2 − l′
We now have two quadratic equations to solve
l2 + l − (l′2 + 3l′ + 2) = 0 ⇒ l = 12 [−1 ± (2l′ + 3)]
l2 + l − (l′2 − l′ ) = 0 ⇒ l = 21 [−1 ± (2l′ − 1)]
Since l, l ≥ 0, the only acceptable solutions are l + l′ = 0 and l = l′ ± 1 as required. However, when
′
l = l′ = 0 the two states have the same (even) parity so the matrix element vanishes by the proof
given in eq (4.41) of the book.
8.8 Tritium, 3 H, is highly radioactive and decays with a half-life of 12.3 years to 3 He by the emission
of an electron from its nucleus. The electron departs with 16 keV of kinetic energy. Explain why
its departure can be treated as sudden in the sense that the electron of the original tritium atom
barely moves while the ejected electron leaves.
Calculate the probability that the newly-formed 3 He atom is in an excited state. Hint: evaluate
h1, 0, 0; Z = 2|1, 0, 0; Z = 1i.
Soln: The binding energy of H is just 13.6 eV and by the virial p theorem its kinetic energy is half
this, so the speed of the ejected electron is larger by a factor 16000/6.8 ≃ 48.5 ≫ 1. Hence the
orbital electron barely moves in the time required for the ejected electron to get clear of the atom.
34
After the decay, the orbital electron is still in the ground state of H. The amplitude for it to be
in the ground state of the new Hamiltonian is h100; Z = 2|100; Z = 1i. In the position representation
this is
3/2 Z
1 4 2
h100; Z = 2|100; Z = 1i = d3 x e−2r/a0 Y00 e−r/a0 Y00
2 a0 a0
23/2 4
Z Z
=4 3 dr r2 e−3r/a0 = 3 23/2 dx x2 e−x (x ≡ 3r/a0 )
a0 3
4 3/2
= 3 2 2!
3
So the probability of being in an excited state
64 × 8
P = 1 − h100; Z = 2|100; Z = 1i2 = 1 − = 0.298
272
solution of this equation changes when l′ changes from l as a result of β increasing from zero to a
small number, show that
2β
∆=1+ 2 + O(β 2 ). (8.7)
4l − 1
Explain the physical significance of this result.
Soln: The given eqn is a quadratic in l′ :
p
−1 ± 1 + 4l(l + 1) − 4β
q
′2 ′ ′
l +l −l(l +1)+β = 0 ⇒ l = ⇒ l′ + 12 = (l + 21 )2 − β, (8.8)
2
where we’ve chosen the root that makes l′ > 0.
Squaring up this equation, we have
(l′ + 21 )2 = (l + 12 )2 − β ⇒ (l′ + ∆ + 21 )2 = (l + 23 )2 − β
Taking the first eqn from the second yields
∆2 + 2(l′ + 12 )∆ = (l + 23 )2 − (l + 12 )2 = 2(l + 1)
This is a quadratic equation for ∆, which is solved by ∆ = 1 when l′ = l. We are interested in
the small change δ∆ in this solution when l′ changes by a small amount δl′ . Differentiating the
equation, we have
2∆δl′
2∆δ∆ + 2∆δl′ + (2l′ + 1)δ∆ = 0 ⇒ δ∆ = −
2∆ + 2l′ + 1
′
Into this we put ∆ = 1, l = l, and by binomial expansion of (8.8)
β
δl′ = −
2l + 1
and have finally
−2β
δ∆ =
(2l + 1)(2l + 3)
Eq (8.51) gives the energy of a circular orbit as
Z02 e2
E=− ,
8πǫ0 a0 (l′ (l) + k + 1)2
with k the number of nodes in the radial wavefunction. This differs from Rydberg’s formula in that
(l′ (l) + k + 1) is not an integer n. Crucially l′ (l) + k does not stay the same if k in decreased by unity
and l increased by unity – in fact these changes (which correspond to shifting to a more circular
orbit) cause l′ (l) + k to increase slightly and therefore E to decrease slightly: on a more circular
orbit, the electron is more effectively screened from the nucleus. So in the presence of screening the
degeneracy in H under which at the same E there are states of different angular momentum is lifted
by screening.
8.13 Show that Ehrenfest’s theorem yields equation (8.70) with B = 0 as the classical equation of
motion of the vector S that is implied by the spin–orbit Hamiltonian (8.71).
Soln: Bearing in mind that L commutes with S and the scalars r and dΦ/dr, Ehrenfest’s theorem
yields
d hSi i dΦ eh̄2 X
ih̄ = h[Si , HSO ]i = − h[Si , Sj Lj ]i
dt dr 2rm2e c2 j
dΦ eh̄2 X
=− h[Si , Sj ]Lj + Sj [Si , Lj ]i
dr 2rm2e c2 j
dΦ eh̄2 X
=− iǫijk hSk Lj i
dr 2rm2e c2
jk
37
P
Cancelling ih̄ from each side and using jk ǫijk Sk Lj = −(S × L)i on the right we obtain
d hSi dΦ eh̄
= hS × Li
dt dr 2rm2e c2
which agrees with (8.70) for Q = −e and B = 0.
9.7∗ The Hamiltonian of a two-state system can be written
A1 + B1 ǫ B2 ǫ
H= , (9.1)
B2 ǫ A2
where all quantities are real and ǫ is a small parameter. To first order in ǫ, what are the allowed
energies in the cases (a) A1 6= A2 , and (b) A1 = A2 ?
Obtain the exact eigenvalues and recover the results of perturbation theory by expanding in
powers of ǫ.
Soln: When A1 6= A2 , the eigenvectors of H0 are (1, 0) and (0, 1) so to first-order in ǫ the perturbed
energies are the diagonal elements of H, namely A1 + B1 ǫ and A2 .
When A1 = A2 the unperturbed Hamiltonian is degenerate and degenerate perturbation theory
applies: we diagonalise the perturbation
B1 ǫ B2 ǫ B1 B2
H1 = =ǫ
B2 ǫ 0 B2 0
The eigenvalues λ of the last matrix satisfy
q
λ2 − B1 λ − B22 = 0 ⇒ λ = 12 B1 ± B12 + 4B22
So
for j = 32 , m = 32
2
2
EB /(µB B) = mgL = 3for j = 23 , m = 12
1
for j = 21 , m = 12
3
with the values for negative m being minus the values for positive m.
(b) We have | 32 , 32 i = |+i|11i so h 23 , 23 |Sz | 32 23 i = 12 and EB /(µB B) = m + hψ|Sz |ψi = 3
2 + 1
2 =2
in agreement with the Landé factor. Similarly
q q
| 32 , 12 i = 13 |−i|11i + 23 |+i|10i ⇒ h 23 , 12 |Sz | 32 , 12 i = 13 (− 12 ) + 32 12 = 61
1
so EB /(µB B) = 2 + 61 = 32 Finally
q q
| 12 , 21 i = 23 |−i|11i − 13 |+i|10i ⇒ h 12 , 12 |Sz | 21 , 12 i = 23 (− 12 ) + 11
23 = − 16
1 1 1
so EB /(µB B) = 2 − 6 = 3
2 2
9.11∗ Using the result proved in Problem 9.10, show that the trial wavefunction ψb = e−b r /2
yields −8/(3π)R as an estimate of hydrogen’s ground-state energy, where R is the Rydberg constant.
2 2 2 2
Soln: With ψ = e−b r /2 , dψ/dr = −b2 re−b r /2 , so
2 4Z ,Z
h̄ b e2
Z
4 −b2 r 2 −b2 r 2 2 2
hHi = dr r e − dr re dr r2 e−b r
2m 4πǫ0
2 Z , Z
h̄ e2 1
Z
4 −x2 −x2 2
= dx x e − dx xe dx x2 e−x
2mb 4πǫ0 b2 b3
Now " #∞
2
e−x
Z
−x2
dx xe = = 21
−2
0
"
−x2
# ∞ √
xe π
Z Z
2
−x2
dx x2 e−x = 1
+ 2 dx e =
−2 4
" #0 ∞ √
3 −x2
x e 3 π
Z Z
2 2
dx x4 e−x = + 23 dx x2 e−x =
−2 8
0
so √ ,√
2
e2 1 3h̄2 b2 e2 b
h̄ 3 π π
hHi = − 2 2 3
= − 3/2
2mb 8 4πǫ0 b 4b 4m 2π ǫ0
At the stationary point of hHi b = me2 /(3π 3/2 ǫ0 h̄2 ). Plugging this into hHi we find
2 2
3h̄2 m2 e4 e2 me2 8 m e 8
hHi = − = − = R
4m 9π 3 ǫ20 h̄4 2π 3/2 ǫ0 3π 3/2 ǫ0 h̄2 3π 2 4πǫ0 3π
39
9.13∗ A particle travelling with momentum p = h̄k > 0 from −∞ encounters the steep-sided
potential well V (x) = −V0 < 0 for |x| < a. Use the Fermi golden rule to show that the probability
that a particle will be reflected by the well is
V2
Preflect ≃ 0 2 sin2 (2ka),
4E
where E = p2 /2m. Show that in the limit E ≫ V0 this result is consistent with the exact reflection
probability derived in Problem 5.10. Hint: adopt periodic boundary conditions so the wavefunctions
of the in and out states can be normalised.
Soln: We consider a length L of the x axis where L ≫ a and k = 2nπ/L, where n ≫ 1 is an
integer. Then correctly normalised wavefunctions of the in and out states are
1 1
ψin (x) = √ eikx ; ψout (x) = √ e−ikx
L L
The required matrix element is
1 L/2
Z a
sin(2ka)
Z
ikx ikx
dx e V (x)e = −V0 dx e2ikx = −V0
L −L/2 −a Lk
so the rate of transitions from the in to the out state is
2π 2π sin2 (2ka)
Ṗ = g(E)hout|V |ini = g(E)V02
h̄ h̄ L2 k 2
2 2 2
Now we need the density of states g(E). E = p /2m = h̄ k /2m is just kinetic energy. Eliminating
k in favour of n, we have
L √
n= 2mE
2πh̄
As n increases by one, we get one extra state to scatter into, so
r
dn L 2m
g= = .
dE 4πh̄ E
Substituting this value into our scattering rate we find
V02 2m sin2 (2ka)
r
Ṗ = 2
2h̄ E Lk 2
This vanishes as L → ∞ because the fraction of the available space that is occupied by the scattering
potential
p is ∼ 1/L. If it is not scattered, the particle covers distance L in a time τ = L/v =
L/ 2E/m. So the probability that it is scattered on a single encounter is
V02 m sin2 (2ka) V2
Ṗ τ = 2 2
= 0 2 sin2 (2ka)
Eh̄ k 4E
Equation (5.78) gives the reflection probability as
(K/k − k/K)2 sin2 (2Ka)
P =
(K/k + k/K)2 sin2 (2Ka) + 4 cos2 (2Ka)
When V0 ≪ E, K − k = 2mV0 ≪ k 2 , so we approximate Ka with ka and then the reflection
2 2
probability becomes
2 2
K − k2
P ≃ sin2 (2ka)
2kK
which agrees with the value we obtained from Fermi’s rule.
40
9.14∗ Show that the number states g(E) dE d2 Ω with energy in (E, E + dE) and momentum in
the solid angle d2 Ω around p = h̄k of a particle of mass m that moves freely subject to periodic
boundary conditions on the walls of a cubical box of side length L is
3 3/2
L m √
g(E) dE d2 Ω = 2E dE dΩ2 . (9.2)
2π h̄3
Hence show from Fermi’s golden rule that the cross section for elastic scattering of such particles by
a weak potential V (x) from momentum h̄k into the solid angle d2 Ω around momentum h̄k′ is
2
m2
Z
2 3 i(k−k′ )·x
dσ = 4 d Ω d xe V (x) . (9.3)
2
(2π) h̄
Explain in what sense the potential has to be “weak” for this Born approximation to the scattering
cross section to be valid.
Soln: We have kx = 2nx π/L, where nx is an integer, and similarly for ky , kz . So each state
occupies volume (2π/L)3 in k-space. So the number of states in the volume element k 2 dkd2 Ω is
3
2 L
g(E)dEd Ω = k 2 dkd2 Ω
2π
Using k 2 = 2mE/h̄2 to eliminate k we obtain the required expression.
From the Fermi rule the probability of making the transition k → k′ is
3
2π 2 2 2π L dk 2
Ṗ = g(E)d Ω|hout|V |ini| = k2 d Ω|hout|V |ini|2
h̄ h̄ 2π dE
The matrix element is
1
Z
′
hout|V |ini = 3 d3 x e−ik ·x V (x)eik·x
L
Now the cross section dσ is defined by Ṗ = dσ × incoming flux = (v/L3 )dσ = (h̄k/mL3 )dσ. Putting
everything together, we find
Z 2 3
h̄k 1 3 −ik′ ·x ik·x 2π L dk 2
k2
dσ = d x e V (x)e d Ω
mL3 L6 h̄ 2π dE
Z 2
mk dk/dE 3 −ik′ ·x ik·x
⇒ dσ = 2 2 d xe V (x)e .
(2π) h̄
Eliminating k with h̄2 k dk = mdE we obtain the desired expression.
The Born approximation is valid providing the unperturbed wavefunction is a reasonable ap-
proximation to the true wavefunction throughout the scattering potential. That is, we must be able
to neglect “shadowing” by the scattering potential.
10.5∗ Assume that a LiH molecule comprises a Li+ ion electrostatically bound to an H− ion, and
that in the molecule’s ground state the kinetic energies of the ions can be neglected. Let the centres
of the two ions be separated by a distance b and calculate the resulting electrostatic binding energy
under the assumption that they attract like point charges. Given that the ionisation energy of Li
is 0.40R and using the result of Problem 10.4, show that the molecule has less energy than that of
well separated hydrogen and lithium atoms for b < 4.4a0 . Does this calculation suggest that LiH is
a stable molecule? Is it safe to neglect the kinetic energies of the ions within the molecule?
Soln: When the LI and H are well separated, the energy required to strip an electron from the Li
and park it on the H− is E = (0.4 + 1 − 0.955)R = 0.445R. Now we recover some of this energy
by letting the Li+ and H− fall towards each other. When they have reached distance b the energy
released is
e2 a0
= 2R
4πǫ0 b b
This energy equals our original outlay when b = (2/0.445)a0 = 4.49a0, which establishes the required
proposition.
41
Since the wavefunction has to be antisymmetric overall, for the l = 1 state it must be symmetric
in the spins of the electrons, so the total spin has to be 1.
10.7∗ Suppose we have three spin-one gyros in a box. Express the state |0, 0i of the box in which
it has no angular momentum as a linear combination of the states |1, mi|1, m′ i|1, m′′ i in which the
individual gyros have well-defined angular momenta. Hint: start with just two gyros in the box,
giving states |j, mi of the box, and argue that only for a single value of j will it be possible to get
|0, 0i by adding the third gyro; use results from Problem 10.6.
Explain the relevance of your result to the fact that the ground state of nitrogen has l = 0.
Deduce the value of the total electron spin of an N atom.
Soln: Since when we add gyros with spins j1 and j2 the resulting j satisfies |j1 − j2 | ≤ j ≤ j1 + j2 ,
we will be able to construct the state |0, 0i on adding the third gyro to the box, only if the box has
j = 1 before adding the last gyro. From Problem 10.6 we have that
1
|0, 0i = √ (|1, −1i|1, 1i − |1, 0i|1, 0i + |1, 1i|1, −1i),
3
where we can consider the first ket in each product is for the combination of 2 gyros and the second
ket is for the third gyro. We use Problem 10.6 again to express the kets of the 2-gyro box as linear
combinations of the kets of individual gyros:
1 1 1
|0, 0i = √ √ |1, −1i|1, 0i − |1, 0i|1, −1i |1, 1i − √ |1, −1i|1, 1i − |1, 1i|1, −1i |1, 0i
3 2 2
1
+ √ |1, 0i|1, 1i − |1, 1i|1, 0i |1, −1i ,
2
1
= √ |1, −1i|1, 0i|1, 1i − |1, 0i|1, −1i|1, 1i − |1, −1i|1, 1i|1, 0i
6
+ |1, 1i|1, −1i|1, 0i + |1, 0i|1, 1i|1, 0i − |1, 1i|1, 0i|1, 0i
This state is totally antisymmetric under exchange of the m values of the gyros.
When we interpret the gyros as electrons and move to the position representation we find that
the wavefunction of the valence electrons is a totally antisymmetric function of their coordinates,
x1 , x2 , x3 . Hence the electrons do an excellent job of keeping out of each other’s way, and this will
be the ground state. To be totally antisymmetric overall, the state must be symmetric in the spin
labels of the electrons, so the spin states will be |+i|+i|+i and the states obtained from this by
application of J− . Thus the total spin will be s = 32 .
10.8∗ Consider a system made of three spin-half particles with individual spin states |±i. Write
down a linear combination of states such as |+i|+i|−i (with two spins up and one down) that is
symmetric under any exchange of spin eigenvalues ±. Write down three other totally symmetric
states and say what total spin your states correspond to.
Show that it is not possible to construct a linear combination of products of |±i which is totally
antisymmetric.
What consequences do these results have for the structure of atoms such as nitrogen that have
three valence electrons?
Soln: There are just three of these product states to consider because there are just three places
to put the single minus sign. The sum of these states is obviously totally symmetric:
1
|ψi = √ |+i|+i|−i + |+i|−i|+i + |−i|+i|+i
3
Three other totally symmetric state are clearly |+i|+i|+i and what you get from this ket and the
one given by everywhere interchanging + and −. These four kets are the kets | 32 , mi.
A totally antisymmetric state would have to be constructed from the same three basis kets used
above, so we write it as
|ψ ′ i = a|+i|+i|−i + b|+i|−i|+i + c|−i|+i|+i
On swapping the spins of the first and the third particles, the first and third kets would interchange,
and this would have to generate a change of sign. So a = −c and b = 0. Similarly, by swapping the
43
spins on the first and second particles, we can show that a = 0. Hence |ψi = 0, and we have shown
that no nonzero ket has the required symmetry.
States that satisfy the exchange principle can be constructed by multiplying a spatial wave-
function that is totally antisymmetric in its arguments by a totally symmetric spin function. Such
states have maximum total spin. In contrast to the situation with helium, conforming states can-
not be analogously constructed by multiplying a symmetric wavefunction by an antisymmetric spin
function.