Separation of Variables
Separation of Variables
Separation of Variables
The linear momenta of the particles can now be expressed where the partial derivatives with respect to r have been
in terms of the rates of change of x and X: replaced by complete derivatives because R depends only
on r. Multiply through by r2/RY and obtain
dx1 dX m1m2 dx
p1 = m1 = m1 +
dt dt m dt 2 2 d2R dR 2 2
2
− r + 2r + Vr − Λ Y = Er 2 (9d)
(5) 2 µ R dr 2 dr 2 µY
dx 2 dX m1m2 dx
p2 = m2 = m2 −
dt dt m dt
The functions Y satisfy the equation
Then, with
Λ2Y = −l ( l +1)Y (10)
1 1 1 mm mm
= + , or µ = 1 2 = 1 2 (6)
µ m1 m2 m1 + m2 m so eqn 9d becomes
Independent of θ, ϕ
it follows that because the Ys cancel
2 2 d2R dR 2 2l(l +1)
p12 p2
2
dX 1 dx
2
− r 2 + 2r + Vr + Y = Er 2 (9e)
+ 2 = 12 m + µ (7) 2 µ R dr dr 2 µY
2m1 2m2 dt 2 dt
By writing P = m(d X /dt ) for the linear momentum of the Step 3 Write an expression for the radial contribution
system as a whole and p = µ(d x /dt ), it follows that
Continuing the argument from step 2, rewrite eqn 9e
P p2 2 by dividing both sides by r2 and multiplying by R. The
E= + + V ( x ) (8a) result is
2m 2 µ
As in Part (b), the linear momenta of the particles are Step 2 Show that the equation separates
expressed in terms of the rates of change of the separation Write the overall wavefunction as the product ψtotal(X,x) =
x and the location of the centre of mass X: ψc.m.(X)ψ(x), where ψc.m.(X) is a function of only the
centre of mass coordinate and ψ(x) is a function of only
dx1 dX m1m2 dx
p1 = m1 = m1 + the separation of the atoms. The overall Schrödinger equa-
dt dt m dt
(14) tion then becomes
dx dX m1m2 dx
p2 = m2 2 = m2 −
dt dt m dt 2 ∂ 2ψ c.m. ( X )ψ ( x ) 2 ∂ 2ψ c.m. ( X )ψ ( x )
− −
2m ∂X 2 2µ ∂x 2 (19a)
Then, with
+ kf x ψ c.m. ( X )ψ ( x ) = Etotalψ c.m. ( X )ψ ( x )
1
2
2
1 1 1 mm mm
= + , or µ = 1 2 = 1 2 (15)
µ m1 m2 m1 + m2 m and therefore
it follows that 2 d 2ψ c.m. ( X ) 2 d 2ψ ( x )
− ψ (x ) − ψ ( X )
2 2
2m dX 2 2 µ c.m. dx 2 (19b)
p12 p2 dX 1 dx
+ 2 = 12 m + µ (16) + 12 kf x 2ψ c.m. ( X )ψ ( x ) = Etotalψ c.m. ( X )ψ (x )
2m1 2m2 dt 2 dt
By writing P = m(dX /dt ) for the linear momentum of the After division by ψtotal(X,x) = ψc.m.( X )ψ ( x ), this equation
molecule as a whole and p = µ(dx /dt ), the total energy becomes
becomes
Depends on X Depends on x
2 2
P p 2 d 2ψ c.m. ( X ) 2 d 2ψ ( x ) 1
Etotal = + + V ( x ) (17a) − − + 2 kf x 2 = Etotal
2m 2 µ 2mψ c.m. ( X ) 2 µψ ( x ) dx 2
dX 2
For a harmonic oscillator (19c)