24.10.2014 Views

My PhD thesis - Condensed Matter Theory - Imperial College London

My PhD thesis - Condensed Matter Theory - Imperial College London

My PhD thesis - Condensed Matter Theory - Imperial College London

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

APPENDIX B. THE CUSP CONDITIONS<br />

may be written as<br />

[− 2<br />

Ĥ({r i }) = ∑ ( 1<br />

∇ 2 i + V (r i ) +<br />

2m<br />

i>2 e 4πɛ 0 2<br />

]<br />

+<br />

[− 2<br />

∇ 2 R + 2V (R) +<br />

4m e<br />

e2<br />

∑<br />

j≠i, j>2<br />

[− 2<br />

m e<br />

∇ 2 r +<br />

1<br />

r ij<br />

+<br />

)]<br />

2<br />

|R − r i |<br />

]<br />

.<br />

e2<br />

4πɛ 0 r + O[r2 ]<br />

(B.5)<br />

For any eigenstate of Ĥ, the divergence in the Coulomb energy must cancel out<br />

exactly; the only other term which can give rise to a matching divergence is<br />

− 2<br />

m e<br />

∇ 2 r.<br />

(B.6)<br />

Therefore, when two electrons are close together, the Schrödinger equation is entirely<br />

dominated by these two terms and is approximately separable; the part of the wave<br />

function depending on r is a solution of the reduced Schrödinger equation<br />

)<br />

(− 2<br />

∇ 2 r +<br />

e2<br />

Ψ r (r) = EΨ r (r).<br />

m e 4πɛ 0 r<br />

(B.7)<br />

This hydrogenic equation has the well-known solution in spherical polar coordinates<br />

Ψ r (r) = ∑ l, m<br />

f lm (r)Y lm (θ, φ)<br />

where the Y lm are spherical harmonics and f lm satisfies the radial equation<br />

(<br />

− 2 d<br />

r 2 df ) ( )<br />

lm e<br />

2<br />

l(l + 1)2<br />

+ + f<br />

m e r 2 dr dr 4πɛ 0 r m e r 2 lm = Ef lm .<br />

Substituting the power series solution<br />

f lm = r p<br />

∞<br />

∑<br />

j=0<br />

a lmj r j<br />

(B.8)<br />

(B.9)<br />

(B.10)<br />

into equation (B.9) gives<br />

r p<br />

∞<br />

∑<br />

j=0<br />

( )<br />

<br />

2 [ ]<br />

a lmj l(l + 1) − (p + j)(p + j + 1) r j−2 +<br />

e2<br />

r j−1 − Er j = 0 (B.11)<br />

m e 4πɛ 0<br />

so that 1 p = l and<br />

( )<br />

me e 2 1<br />

a lm1 = a lm0<br />

4πɛ 0 2 2(l + 1) . (B.12)<br />

1 The alternative result, p = −l − 1, is rejected on physical grounds.<br />

193

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!