Rachunek zaburzeń dla stanów zdegenerowanych
Tak jak poprzednio poszukujemy rozwiązania równania ϕ
ϕ E Hˆ = ,
znając ścisłe rozwiązania równania o zbliżonym hamiltonianie Hˆ(0). Teraz jednak stan niezaburzony jest zdegenerowany. Przyjmujemy, że mamy n stanów opisywanych funkcjami falowymi ϕ1(0),ϕ(20),Kϕ(n0), o energii E(0). Hamiltonian Hˆ i poszukiwane funkcje własne
ϕ
przestawiamy w postaci. , , 2 , 1 ,
ˆ , ˆ
ˆ H(0) H(1) (0) (1) i n H = + ϕi =ϕi +ϕi = K
Stany niezaburzoneϕ(i0) są zdegenerowane. Zaburzenie może spowodować usunięcie degeneracji, wówczas zaburzone energie odpowiadające różnym stanom będą różne. Stany zaburzone ϕi, mogą również pozostać zdegenerowane, tzn. poprawki do energii będą takie same dla wszystkich stanów. Możliwy jest także sytuacja pośrednia – zaburzenie usuwa degenerację częściowo.
Zakładamy, jak poprzednio, że Hˆ(1) jest dużo mniejsze niż Hˆ(0), ϕ(i1) jest dużo mniejsze niż ϕ(i0) oraz E(i1) jest dużo mniejsze od E(0).
Równanie, które chcemy rozwiązać, przyjmuje postać1
(
Hˆ(0)+Hˆ(1)) (ϕ(i0) +ϕ(i1)) (
= E(0) +E(i1))(
ϕ(i0) +ϕ(i1))
.
Pamiętając, że Hˆ(0)ϕ(i0) = E(0)ϕ(i0) i pomijając wyrazy „kwadratowo małe”
dostajemy
i i i i
i H E E
Hˆ(0)ϕ(1)+ ˆ(1)ϕ(0) = (0)ϕ(1) + (1)ϕ(0).
Funkcje tworzące lewą i prawą stronę równania mnożymy teraz skalarnie przez funkcję ϕ(j0)i dostajemy
(
ϕ(j0),Hˆ(0)ϕ(i1) +Hˆ(1)ϕ(i1))
=(
ϕ(j0),E(0)ϕ(i1) +E(i1)ϕ(i0))
.Pamiętając, że iloczyn skalarny jest linowy w drugim argumencie mamy
(
ϕ(0j),Hˆ(0)ϕ(i1)) (
+ ϕ(j0),Hˆ(1)ϕ(i0))
=(
ϕ(0j),E(0)ϕ(i1)) (
+ ϕ(j0),E(i1)ϕ(i0))
.Prawą stronę równania zapisujemy jako
(
ϕ(j0),E(0)ϕ(i1)) (
+ ϕ(j0),E(i1)ϕ(i0))
=E(0)(
ϕ(j0),ϕ(i1))
+E(i1)(
ϕ(j0),ϕ(i0))
. Ponieważ operator Hˆ(0) jest hermitowski, więc(
ϕ(j0),Hˆ(0)ϕ(i1)) (
= Hˆ(0)ϕ(j0),ϕ(i1))
=E(0)(
ϕ(j0),ϕ(i1))
.Uwzględniając jeszcze, że funkcje falowe ϕ(i0)są unormowane i wzajemnie ortogonalne tzn.
(
ϕ(i0),ϕ(j0))
=δij, otrzymujemy(
j i) ( j H i )
E (
j i )
Ei ij
E(0) ϕ(0),ϕ(1) + ϕ(0), ˆ(1)ϕ(0) = (0) ϕ(0),ϕ(1) + (1)δ . Równanie na poprawki do energii niezaburzonej przyjmuje postać.
(
ϕ(j0),Hˆ(1)ϕ(i0))
−E(i1)δij =0. Widzimy, że jeśli i≠ j, mamy(
ϕ(j0),Hˆ(1)ϕ(i0))
=0.Czyli wyprowadzone równanie ma sens tylko wtedy, gdy niediagonalne wyrazy macierzy
(
ϕ(j0),Hˆ(1)ϕ(i0))
znikają tzn.(
ϕ(j0),Hˆ(1)ϕ(i0))
=0 dla i≠ j. Wówczas(
i i)
i H
E(1) = ϕ(0), ˆ(1)ϕ(0) .
Otrzymaliśmy więc wyrażenie bardzo zbliżone do formuły odpowiadającej przypadkowi stanów niezdegenerowanych. Bardzo podobnie też znajdujemy poprawki do funkcji falowych. Warto też podkreślić, że warunek
(
ϕ(j0),Hˆ(1)ϕ(i0))
=0 dla i≠ j nie rozstrzyga, czy zaburzenie usuwa degenerację czy nie.Efekt skończonych rozmiarów jądra
Obliczymy teraz przesunięcia wzbudzonych poziomów atomu wodoru, które, jak wiadomo, są zdegenerowane. Wyprowadzony poprzednio Hamiltonian zaburzający równy jest
≥
<
−
− −
=
. ,
0
, 1,
2 2
3
ˆ 3
2
2 ) 1 (
R r
R r r
R r e R
H
Rozważmy n-ty poziom, który zdegenerowany jest n2 krotnie ze względu na liczby kwantowe l i m. Wyprowadzony wzór na poprawkę do energii ma zastosowanie, jeśli
(
ϕnlm,Hˆ(1)ϕnl'm')
=0 dla l≠l' lub m≠m'. Ponieważ Hˆ(1) zależy tylko od zmiennej radialnej, więc ortogonalność harmonik sferycznych gwarantuje spełnienie tego warunku. A zatem poprawkę do energii stanu (n,l,m)obliczamy ze wzoru
(
, ˆ)
( ) ˆ(1) ( )* 3 )
1 ( )
1
( H d r r H r
Enlm nlm nlm nlm r nlm r ϕ ϕ
ϕ
ϕ =
∫
= .
Przedstawiając funkcję falową atomu wodoru w postaci
) , ( ) ( ) , ,
( θ ϕ θ ϕ
ϕnlm r =Rnl r Ylm ,
mamy
) ˆ (
) ( )
, ( ) ,
( (1)
0 2
* 2 )
1
( d Y Y drr R r H R r
Enlm =
∫
Ω lm θ ϕ lm θ ϕ∫
∞ nl nl .Ponieważ harmoniki sferyczne są unormowane tzn.
∫
d2ΩYlm*(θ,ϕ)Ylm(θ,ϕ)=1, znajdujemy) ˆ (
) ( (1)
0 2 )
1
( drr R r H R r
Enl =∞
∫
nl nl .Poprawka nie zależy od liczby kwantowej m, więc zburzenie spowodowane skończonymi rozmiarami jądra atomowego nie powoduje usunięcia degeneracji poziomów atomu ze względu na liczbę kwantową m.
Efekt Zeemana
Przyjmując, jak poprzednio, że pole magnetyczne jest jednorodne i skierowane wzdłuż osi z, hamiltonian zaburzający ma postać
z e
cL m
H eB ˆ
2
ˆ(1) = .
Warunek
(
ϕnlm,Hˆ(1)ϕnl'm')
=0 dla l≠l' lub m≠m' spełniony jest trywialnie, bowiem ϕnl'm' jest funkcja własnąLˆz(z wartością własną m'h). A zatem(
(1) ' ') (
' ')
' ' '2 , ˆ
2
, ˆ ll mm
e m
nl z nlm e m
nl
nlm m
c m L eB
c m
H ϕ eB ϕ ϕ δ δ
ϕ = = h .
Poprawkę do energii znajdujemy jako
( ) ( )
mhc m L eB
c m H eB
E
e nlm
z nlm e nlm
nlm nlm
2 , ˆ
2 , ˆ(1)
) 1
( = ϕ ϕ = ϕ ϕ = .
Widzimy, że poprawka zależy tylko od liczby kwantowej m, więc efekt Zeemana usuwa degenerację ze względu na liczbę kwantową m, pozostawiając degenerację ze względu na liczbę kwantową l.
Efekt Starka
Przyjmując, jak poprzednio, że pole elektryczne jest jednorodne i skierowane wzdłuż osi z, hamiltonian zaburzający ma postać
z e Hˆ(1) = E . Obliczamy
(
, ˆ)
( ) ˆ ( ) ( ) ' '( )* 3 '
' ) 1 (
* 3 '
' ) 1
( d r r H r e d rz r r
H nlm nlm nlm nlm nlm
nlm
r r
r
r ϕ ϕ ϕ
ϕ ϕ
ϕ =
∫
= E∫
.Przedstawiając funkcję falową w postaci ϕnlm(r,θ,ϕ)=Rnl(r)Ylm(θ,ϕ) i pamiętając, że z=rcosθ, znajdujemy
(
, ˆ)
cos ( , ) ( , ) ( ) '( )0 3 '
'
* 2
' ' ) 1
( e d Y Y drr R r R r
H nlm lm lm nl nl
nlm ϕ =
∫
Ω θ θ ϕ θ ϕ∫
∞ϕ E .
Jak widzimy, w ogólności
(
ϕnlm,Hˆ(1)ϕnl'm')
≠0 dla l≠l' lub m≠m',więc wyprowadzone formuły na poprawkę do energii nie stosują się do efektu Starka.
Przypadek: (ϕ(j0),Hˆ(1)ϕ(i0))
≠0.
Fakt, że
(
ϕ(j0),Hˆ(1)ϕ(i0))
≠ 0 , jest sygnałem, rachunek zaburzeń prowadzimy dla niewłaściwie wybranych stanów zdegenerowanych. Należy pamiętać, że zamiast stanów (n0)2 ) 0 ( 1
) 0
( ,ϕ , ϕ
ϕ K o energii E(0) możemy rozważać ich kombinacje liniowe tzn. stany ϕ~(10),ϕ~(20),Kϕ~(n0) dane wzorem
∑
=
=
n
j
j j i i
1 ) 0 ( )
0 (
~ α ϕ
ϕ
(αij są współczynnikami liczbowymi), które również mają energię E(0). Jeśli zaburzenie usuwa degenerację, to zaburzone stany niezdegenerowane mogą odpowiadać kombinacji liniowej niezaburzonych stanów zdegenerowanych.
Rachunek zaburzeń więc prowadzimy dla stanów (n0) 2
) 0 ( 1
) 0 (
, ~ , ~
~ ϕ ϕ
ϕ K i sam
rachunek zdeterminuje wartości współczynników αij. Tak zatem poszukujemy rozwiązania równania
ϕ ϕ E Hˆ = ,
znając ścisłe rozwiązania równania o zbliżonym hamiltonianie Hˆ(0). Stan niezaburzony jest zdegenerowany. Przyjmujemy, że mamy n stanów opisywanych funkcjami falowymi ϕ~(10),ϕ~(20),Kϕ~(n0) o energii E(0). Hamiltonian
Hˆ i poszukiwane funkcje własne
ϕ
przestawiamy w postaci . , , 2 , 1~ , ˆ ,
ˆ
ˆ H(0) H(1) (0) (1) i n H = + ϕi =ϕi +ϕi = K
Zakładamy, jak poprzednio, że Hˆ(1) jest dużo mniejsze niż Hˆ(0), ϕ(i1) jest
dużo mniejsze niż (i0)
ϕ~ oraz E(i1) jest dużo mniejsze od E(0).
Równanie, które chcemy rozwiązać, przyjmuje postać
(
H(0) H(1)) ( (i0) (i1)) (
E(0) E(i1))(
(i0) (i1))
~ ˆ ~
ˆ + ϕ +ϕ = + ϕ +ϕ .
Pamiętając, że Hˆ(0)ϕ~(i0) =E(0)ϕ~(i0) i pomijając wyrazy „kwadratowo małe”
dostajemy
i i i i
i H E E
Hˆ(0)ϕ(1)+ ˆ(1)ϕ~(0) = (0)ϕ(1)+ (1)ϕ~(0).
Funkcje tworzące lewą i prawą stronę równania mnożymy teraz skalarnie przez funkcję ϕ(j0)i dostajemy
(
ϕj ,Hˆ ϕi +Hˆ ϕ~i)
=(
ϕj ,E ϕi +Ei ϕ~i)
Pamiętając, że iloczyn skalarny jest linowy w drugim argumencie, mamy
(
(j0) H(0) (i1)) (
(j0) H(1) (i0)) ( (j0) E(0) (i1)) (
(j0) E(i1) (i0))
, ~
~ , , ˆ
, ˆ ϕ ϕ ϕ ϕ ϕ ϕ ϕ
ϕ + = + .
Prawą stronę równania zapisujemy jako
(
ϕ(j0),E(0)ϕ(i1)) (
+ ϕ(j0),E(i1)ϕ(i0))
=E(0)(
ϕ(j0),ϕ(i1))
+E(i1)(
ϕ(j0),ϕ~(i0))
. Ponieważ operator Hˆ(0) jest hermitowski, więc(
ϕ(j0),Hˆ(0)ϕ(i1)) (
= Hˆ(0)ϕ(j0),ϕ(i1))
=E(0)(
ϕ(j0),ϕ(i1))
. Dostajemy zatem(
(j0) H(1) (i0))
E(i1)(
(j0) (i0))
,~ ˆ ~
, ϕ ϕ ϕ
ϕ = .
Zakładając, że funkcje falowe ϕ(i0)są unormowane i wzajemnie ortogonalne, tzn.
(
ϕ(i0),ϕ(j0))
=δij, obliczamy(
(j0) (i0))
,ϕ~
ϕ oraz
(
(j0) H(1) (i0))
ˆ ~
, ϕ
ϕ :
( ) ( )
ijjk n
k k i k
j n
k k i n
k
k k i j
i
j ϕ ϕ α ϕ α ϕ ϕ α δ α
ϕ = = =
=
∑ ∑ ∑
=
=
= 1
) 0 ( ) 0 ( 1
1
) 0 ( )
0 ( )
0 ( ) 0
( ,~ , , ,
( ) ∑ ∑ ( )
=
=
=
=
n
k
k i k j
n
k
k k i j
i
j H H H
1
) 0 ( ) 1 ( ) 0 ( 1
) 0 ( )
1 ( ) 0 ( )
0 ( ) 1 ( ) 0
( , ˆ ϕ~ ϕ , ˆ α ϕ ϕ , ˆ ϕ α
ϕ .
Po wprowadzeniu oznaczenia W jk ≡
(
ϕ(j0),Hˆ(1)ϕ(k0))
dostajemy) 0
1 ( 1
=
∑
−=
j i i n
k
k i
jk E
W α α ,
co można przepisać w postaci
0 )
(
1
) 1
( =
∑
−=
k i n
k
jk i
jk E
W δ α .
Przy określonym i, jest to jednorodne równanie na n współczynników αik, .
, 2 ,
1 n
k = K Równanie jednorodne AX =0, gdzie A jest macierzą a X szukanym wektorem, ma rozwiązanie, jeśli detA=0. A zatem równanie na poprawki do energii E przybiera postać (1)
0
)
det( W − E
(1)I =
,Efekt Starka dla pierwszego stanu wzbudzonego atomu wodoru
Przyjmując, jak poprzednio, że pole elektryczne jest jednorodne i skierowane wzdłuż osi z, hamiltonian zaburzający ma postać
z e Hˆ(1) = E .
Funkcje falowe niezaburzonego pierwszego stanu wzbudzonego atomu wodoru mają zwykłą postać ϕnlm(r,θ,ϕ)=Rnl(r)Ylm(θ,ϕ), przy czym
B B
a r
B B
a r
B B
a e r a r
R
a e r a
r R
2 3
21 3
2 3
20 3
2 3 ) 1 (
2 2 ) 1 (
−
−
=
−
=
θ ϕ
ϕ π θ
π θ ϕ
θ ϕ π θ
e i
Y Y Y
±
± =
=
=
8 sin ) 3
, (
4 cos ) 3
, (
4 ) 1 , (
1 1 10 00
m
Obliczamy teraz
(
ϕ2lm,Hˆ(1)ϕ2l'm')
, pamiętając, że z=rcosθ. Wprowadzamy oznaczenia:ϕ
1 ≡ϕ
200,ϕ
2 ≡ϕ
21−1,ϕ
3 ≡ϕ
210,ϕ
4 ≡ϕ
211 oraz(
i j)
ij H
W ≡ ϕ, ˆ(1)ϕ . Ponieważ Hˆ(1) nie zależy od kąta azymutalnego, więc ortogonaloność funkcji
4 3 2
1,
ϕ
,ϕ
,ϕ
ϕ
ze względu na liczbę m sprawia, że43 0
34 42 24 32 23 41 14 21
12 =W =W =W =W =W =W =W =W =W =
W .
Łatwo też zauważyć, że znikają wszystkie diagonalne wyrazy macierzy W ij
44 0
33 22
11=W =W =W =
W ,
bowiem
0 cos
sin cos cos
0
3
0 0
2
∫
∫
=∫
=π =π π
θ θ
θ θ
θ θ
θ d d
d .
Jedyne nieznikające wyrazy macierzy W to ij W i 13 W , przy czym 31 W13 =W31, ponieważ funkcje
ϕ
1,ϕ
3 są rzeczywiste. Obliczamy. 3 1 2
cos 3 sin
) ( ) ( )
, ( ) , ( cos
/ 4
2 2
21 20 0
3 10
* 00 2
13
a
r e a
r e r
dr d
d e
r R r R r dr Y
Y d
e W
B E
E E
−
=
−
=
Ω
=
−
∞
∞
∫
∫
∫
∫
∫
π π
θ θ
θ ϕ
ϕ θ ϕ θ θ
Równanie na zaburzenia energii ma postać
(
( ) ( ))
( ) 0) ( 0
0 0
0
0 0
0
0 0
0 0
0
0 0
2 ) 1 ( 2 13 2 ) 1 ( ) 1 (
) 1 ( 13
) 1 (
13 )
1 (
) 1 ( ) 1 ( 31
) 1 (
13 )
1 (
=
−
=
−
−
−
−
=
−
−
−
−
E W E
E E W
E W E
E E
W
E W E
.
Rozwiązaniami są E(1) =0 oraz E(1) = ±W13 = m3eEaB. A zatem, pierwszy poziom wzbudzony atomu wodoru o energii
2 4
E =−Ry, gdzie 2
4
2h e
Ry≡ me jest stałą Rydberga, rozszczepia się w polu elektrycznym na trzy poziomy E2 −3eEaB, E2
oraz E2+3eEaB . Poziomowi E2 odpowiada kombinacja liniowa stanów
1 21 2 ≡
ϕ
−ϕ
iϕ
4 ≡ϕ
211 zaś poziomom E2±3eEaB kombinacja liniowaϕ
1 ≡ϕ
200i