MIT5 61F13 Lecture8
MIT5 61F13 Lecture8
MIT5 61F13 Lecture8
Last time
functions)
mm
* μ = 1 2
reduced mass
m1 + m2
d2x k
From F = ma we get 2 = − x [we get x(t) from this, not from ψ(x)]
dt m
1/2
⎛ 2E ⎞
turning point x± (E) = ± ⎜ from E = V ( x± (E))
⎝ k ⎟⎠
ν,ω,τ
T (t),T (kinetic energy)
Today
2
� = T� + V� = p̂ + 1 kx̂ 2
H
2μ 2
� 2 ∂ 1 2
=− + kx̂
2 μ ∂x 2 2
where α = ( kμ)
1/2
ξ = α1/2 x � (ξ is dimensionless,
1/2
� 2 ∂ 2 � 2 ∂ 2 � ⎛ k ⎞ ∂2 � ∂ 2
= α 2= ⎜ ⎟ = ω
2 μ ∂x 2
2 μ ∂ξ 2 ⎝ μ⎠ ∂ξ 2 2 ∂ξ 2
⎡ ∂ ∂ ∂ξ ⎤
⎢ because =
∂x ⎥
∂x ∂ξ ax
⎢⎣ ⎥
α1/2 ⎦
1 2 1 ⎛ 1 ⎞
2
kx = k ⎜ ⎟ ξ
2 2 ⎝ α ⎠
1/2
1⎛ k⎞ 1
= ⎜ ⎟ �ξ 2 = �ωξ 2
2 ⎝ μ⎠ 2
�=�
ω ⎡ ∂ 2 ⎤ ( much simpler form)
H ⎢ − 2 + ξ2 ⎥
2 ⎣ ∂ξ ⎦
� ψ(ξ) = Eψ ( ξ )
H
⎡ ∂2 2E ⎤ (entire differential equation, except
0 = ⎢− 2 + ξ2 − ψ
⎣ ∂ξ �ω ⎦⎥ ψ, is dimensionless)
One can convert this into the Hermite differential equation by making the substitution
ψ ( ξ ) = e− ξ f ( ξ ) and finding a new differential equation for f(ξ). The reason for doing this is
2
2
that e− ξ 2 ensures that ψ → 0 as |ξ| → ∞. Note that letting ξ2 → ∞ means that 2E/�ω is
2
negligible with respect to ξ2. What is the solution to the differential equation if we ignore the
2E
term?
�ω
This is a very clean form of the Schrödinger equation because all of the k,μ-specific factors are
absorbed into a dimensionless ξ variable. Why would we want this?
The Hermite polynomials (in integer powers of ξ) are solutions to the differential equation
d2Hn dH n
2
− 2ξ + 2nH n = 0 . (Hermite equation)
dξ dξ
There are very convenient “recursion relations” that relate the Hn–1 to the Hn, etc.
dH n
= 2nH n−1 ( ξ ) .
dξ
This looks like the effect of p̂ on Hn (because the leading term in Hn is ξn).
1
ξH n (ξ) =H (ξ) + nH n−1 (ξ).
2 n+1
This looks like the effect of x̂ on Hn. We will use this second recursion relation to compute
integrals of the form
∫ ψ v ( x̂ n p̂ m ) ψ v+f dx
(we will postpone the actual evaluation until next lecture) and we derive the “selection rule” for
nonzero integrals
f = n + m, n + m – 2, … – (n + m)
d n − ξ2
H n ( ξ ) = ( −1) en ξ2
e
dξ n
1/2
⎛ kμ ⎞ Gaussian
α=⎜ 2⎟
⎝ !� ⎠ envelope
ensures
ξ = α1/2x ψ→0
1/4 assx → o
±∞
⎛α⎞ 1
ψ v (x) = ⎜ ⎟ H v (ξ ) e −ξ2 2
( 2 v v!) ⎝ π ⎠ a x
1/2
!#########"
vibrational Nv Hermite polynomials
quantum
number Normalization
v = 0, 1, 2, …
* Normalized
* ψv(±∞) = 0
* ψv(0) = 0 for odd-v (odd function)
dψ v
* = 0 for even-v (even function), to be worked out
dx x=0
Ev = �ω(v + ½)
Semi-Classical Picture – applicable to all 1-dimensional problems for insight, prediction, and
obtaining ψ(x) without solving any differential equation.
p(x) is momentum in classical mechanics but just a convenient function in quantum mechanics.
h
Quantum: d e Broglie λ =
p
valid not just for free particle or a piecewise constant potential
=
h
Semi-classical: λ(x) ≡
pclassical (x)
* outer lobes (near x+ and x–) are largest [see McQuarrie, page 226, Fig. 5.10 right side]
* envelope within classically allowed region resembles what you expect from classical
mechanics
dx
ψ * (x)ψ(x)dx ∝ (slow speed↔high probability)
v classical
1
vclassical = pclassical μ = {2 μ[ E − V (x)]}
1/2
μ
To get the proportionality constant, consider the fraction of time the particle is
found between x and x + dx:
time ( x → x + dx ) ⎧⎪ probability of finding particle moving ⎫⎪ dx vclassical ( x )
= ⎨ ⎬=
time ( x− to x+ ) ⎪⎩ to right between x and x + dx ⎪⎭ τ/2
1/2
1 2π ⎛ μ⎞
τ= = = 2π ⎜ ⎟ ⎡v
⎣ = ω / 2π, ω = [k / μ] ⎤⎦
1/2
ν ω ⎝ k⎠
−1/2 1/2
dx dx 1 ⎛ μ⎞ ⎡ k / 2π 2 ⎤
ψ * (x)ψ(x)dx = = ⎜ ⎟ =⎢ ⎥ dx
( vclassical ) ( τ 2 ) ⎧2
1/2
⎫ π⎝ k⎠ ⎣ E − V (x) ⎦
⎨ [ E − V (t)]⎬
⎩μ ⎭
1/2
⎡ k / 2π 2 ⎤
ψ * (x)ψ(x) = ⎢ ⎥ gives the classical average of ψ * (x)ψ(x) near x
⎣ E − V (x) ⎦
(but not the phase). [To get the classical envelope, assume that the maximum
value of ψ*ψ is twice the average value. [This is always a good approximation
for a rapidly oscillating always positive function.] Thus the envelope of
⎡ 2k / π 2 ⎤ ⎤
1/2
ψ∗(x)ψ(x) envelope is ⎢ ⎥ ⎥.
⎣ E − V (x) ⎦ ⎥⎦
* node spacing
h
(1) recall λ = , use classical ideas to qualitatively locate nodes,
p(x)
(2) nodes are closest together when p is largest (near x = 0), envelope and
node spacing allows you to sketch ψ∗(x)ψ(x) without solving a differential
equation.
Δ x 1
(3) or compute a “phase integral” = (want to find value of Δx that is
λ(x) 2
equal to λ/2, the distance between nodes)
2 x2
h ∫x1
phase integral pE (x)dx tells us how many nodes there are between x1 and x2
at energy E. This is the same as knowing how many bound energy levels lie at or
below E.
Ev
= ωr ( v + 1 / 2 ) − ωx ( v + 1 / 2 )2 (~ means cm–1 units)
hc
“anharmonicity”, comes mostly
from x3 and x4 terms in V(x)
Spectral intensities
vf
2
I fi ∝ ∫ dx ψ μ( x ) ψ vi
*
hν
v f
electric dipole moment vi
(HCl vs. H2, N2, O2, Cl2)
(think of radio antenna)
Δv = 0 Δv = ±1 Δv = 0, ±2
pure rotation vibrational vibrational
spectrum fundamental overtone
μ
+ μ1 ∫ dx ψ v*f xψ vi + 2 ∫ dx ψ v*f x 2ψ vi
�j j jj 2
see recursion
relationship
2 ψ ψ
xψ v v+1 v−1
selection rule: Δv = ±1
1
ξ2 H n = ξH n+1 + nξH n−1
2
1⎛ 1 ⎞ ⎛1 ⎞
= ⎜ H n+2 + (n + 1)H n ⎟ + n ⎜ H n + (n − 1)H n−2 ⎟
2⎝ 2 ⎠ ⎝2 ⎠
Δ v = 0, ±2
Next time: a†, a treatment
For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.