Przejdź do zawartości

Mechanika kwantowa/Wstęp do teorii promieniowania kwantów pola elektromagnetycznego

Z Wikibooks, biblioteki wolnych podręczników.
Mechanika kwantowa
Mechanika kwantowa
Wstęp do teorii promieniowania kwantów pola elektromagnetycznego

Licencja
Autor: Mirosław Makowiecki
Absolwent UMCS Fizyki Komputerowej Uniwersytetu Marii Curie-Skłodowskiej w Lublinie
Email: miroslaw(kropka)makowiecki(małpa)gmail(kropka)pl
Dotyczy: książki, do której należy ta strona, oraz w niej zawartych stron i w nich podstron, a także w nich kolumn, wraz z zawartościami.
Użytkownika książki, do której należy ta strona, oraz w niej zawartych stron i w nich podstron, a także w nich kolumn, wraz z zawartościami nie zwalnia z odpowiedzialności prawnoautorskiej nieprzeczytanie warunków licencjonowania.
Umowa prawna: Creative Commons: uznanie autorstwa, na tych samych warunkach, z możliwością obowiązywania dodatkowych ograniczeń.
Autor tej książki dołożył wszelką staranność, aby informacje zawarte w książce były poprawne i najwyższej jakości, jednakże nie udzielana jest żadna gwarancja, czy też rękojma. Autor nie jest odpowiedzialny za wykorzystanie informacji zawarte w książce, nawet jeśli wywołaby jakąś szkodę, straty w zyskach, zastoju w prowadzeniu firmy, przedsiębiorstwa lub spółki bądź utraty informacji, niezależnie czy autor (a nawet Wikibooks) został powiadomiony o możliwości wystąpienie szkód. Informacje zawarte w książce mogą być wykorzystane tylko na własną odpowiedzialność.


W tym rozdziale zapoznamy się z podstawowymi własnościami dla operatorów symbolizujących kwanty pola elektromagnetycznego, wyprowadzimy macierze spinowe, podamy równanie własne współrzędnej zetowej i kwadratu całkowitego momentu pędu kwantów pola elektromagnetycznego.

Całkowity moment pędu kwantu pola elektromagnetycznego a jego obroty

[edytuj]

Wiemy, że obroty są odpowiedzialne za zasadę zachowania momentu pędu. Jeśli znamy oś obrotu i jej kierunek, czyli znamy wektor n, który właśnie za nie jest odpowiedzialny, przy czym zakładamy, że zwrot tego wektora jest do góry względem pewnej płaszczyzny i prostopadły do niego, w której następuje obrót, jeśli obrót jest przeciwny z kierunkiem ruchu zegara, zatem macierz obrotu przedstawia się:

R(n,θ)=eθnl̂
(33.1)

Wektor wodzący tuż po obrocie względem operatora obrotu zdefiniowaną w punkcie (33.1) względem starego wektora zapisujemy:

r'=R(n,θ)r
(33.2)

Można też napisać transformację odwrotną do (33.2) transformująca nowy wektor tuż po obrocie na stary wektor wodzący przed obrotem względem wektora n:

r=R1(n,θ)r'
(33.3)
(Rys. 33.1) Pole wektorowe przed i po obrocie.

Transformacja funkcji falowej przedstawia się podobnie jak w przypadku współrzędnych wektora wodzącego (33.2), jeśli rozpatrujemy funkcje falowe w postaci wektora, to taką transformacje zapisujemy podobnie:

ψ'(r')=R(n,θ)ψ(r)
(33.4)

Jeśli dodatkowo założymy, że kwant pola elektromagnetycznego posiada dodatkowy moment pędu i nazwijmy ją spinowym momentem pędu kwantu promieniowania, to operator transformacji w analogii do (33.1) jest:

D̂(1)=eiθnŜ
(33.5)
  • gdzie Ŝ - jest to macierz spinowa:

Całkowity operator obrotu jest operatorem iloczynu dwóch operatorów, tzn. (33.1) oraz (33.5):

P̂=D̂(1)T̂=eiθ(nŜ)eiθnl̂=eiθn(l̂+Ŝ)
(33.6)

Zatem otrzymaliśmy (33.6), czyli macierz obrotu sumy operatorów momentu pędu orbitalnego i operatora momentu pędu spinowego:

P̂=eiθn(l̂+Ŝ)
(33.7)

Nazwijmy całkowitym operatorem momentu pędu jako suma operatora orbitalnego momentu pędu i operatora momentu pędu spinowego kwantu promieniowania:

ĵ=l̂+Ŝ
(33.8)

Całkowita macierz obrotu (33.7), korzystając z definicji całkowitego momentu pędu kwantu (33.8), jest przedstawiona:

P̂=eiθnĵ
(33.9)

Rozpatrując wartości własne operatorów momentu pędu, to całkowity moment pędu kwantu promieniowania jest zapisany wedle wzoru:

j=l+S
(33.10)

Wyznaczenie macierzy spinowych kwantów pola elektromagnetycznego

[edytuj]
(Rys. 33.2) Zmiana δ stałego pola pola przy obrocie układu o kąt δθ

W macierzy transformacji operatora całkowitego momentu pędu kwantu promieniowania napisaną według (33.7), jeśli w nim dokonamy operacji θδθ, i rozpatrując przy stałym polu wektorowym ψ(r)=const, wtedy wynik działania operatora momentu pędu na tą funkcję jest równa zero (ten operator jest w coś rodzaju liczenia pochodnej względem pewnych zmiennych przestrzennych według definicji operatora momentu pędu) i pozostaje nam tylko operator zależny od operatora spinu (33.5), a więc funkcję ψ możemy rozłożyć w szereg Taylora jednej zmiennej pomijając wyższe wyrazy inne niż liniowe:

ψk''(r)=ψk'(r)kiδθn(l̂+Ŝ)k'kψk(r)ψk''(r)=ψk'(r)kiδθ1nl̂ψk(r)0+
kiδθnŜk'kψk(r)
(33.11)

Przy zastosowania operatora (33.5) dla infitezymalnego kąta obrotu przy transformacji wektorowej funkcji własnej jakiegoś równania własnego otrzymujemy to samo równanie, co (33.11) przy stałej tej funkcji:

δψ'(r)k'=ψk''(r)ψk'(r)=kiδθnŜk'kψk(r)
(33.12)

Z drugiej jednak strony z rysunku obok możemy jednak zapisać przy pomocy wektora n, wzdłuż której następuje obrót o kierunku odwrotnym względem wskazówek zegara, jeśli zwrot tego wektora jest nad zegarem:

δψ=(ψ×n)δθ
(33.13)

Ponieważ wzory (33.12) i (33.13) przedstawiają to samo, dla tej samej współrzędnej wektora falowego, ale dla stałego pola wektorowego, przyrównujemy obie strony tychże równań do siebie:

kiδθnŜk'kψk(r)=(ψ×n)k'δθ(kinŜk'kψk(r))δθ=(ψ×n)k'δθ
(33.14)

I ostatecznie w ostatnim równaniu (33.14) dokonajmy operacji wymnożenia obustronnego przez stałą kreśloną Plancka, wtedy otrzymamy tożsamość, którą później wykorzystamy dla ściśle określonych k':

iknŜk'kψk=(ψ×n)k'ik(nŜk'k)ψk=(ψ×n)k'
(33.15)

Rozpatrujemy równanie (33.15) dla parametru k' równej jeden k'=1:

i[nx(Sx)11+ny(Sy)11+nz(Sz)11]ψ1+i[nx(Sx)12+ny(Sy)12+nz(Sz)12]ψ2+i[nx(Sx)13+ny(Sy)13+nz(Sz)13]ψ3=
=(ψ2nzψ3ny)
(33.16)

Porównując obie strony tożsamości (33.16) dochodzimy do wniosku, że odpowiednie współrzędne operatora spinu kwantu promieniowania, a właściwie jej niektóre elementy można policzyć z pierwszego z trzech równań, które będziemy rozważać:

{(Sx)11=(Sy)11=(Sz)11=0i[nx(Sx)12+ny(Sy)12+nz(Sz)12]=nzi[nx(Sx)13+nz(Sy)13+nz(Sz)13]=ny{(Sx)12=(Sy)12=0(Sz)12=i(Sx)13=(Sz)13=0(Sy)13=i
(33.17)

Gdy rozpatrzymy równanie (33.15) dla parametru k'=2, wtedy niektóre elementy macierzy współrzędnych operatora spinu można policzyć z drugiego z trzech równań, wtedy:

i[nx(Sx)21+ny(Sy)21+nz(Sz)21]ψ1+i[nx(Sx)22+ny(Sy)22+nz(Sz)22]ψ2+i[nx(Sx)23+ny(Sy)23+nz(Sz)23]ψ3=
=[ψ3nxψ1nz]
(33.18)

Porównując obie strony tożsamości (33.18), to dochodzimy do wniosku:

{i[nx(Sx)21+ny(Sy)21+nz(Sz)21]=nz(Sx)22=(Sy)22=(Sz)22=0i[nx(Sx)23+ny(Sy)23+nz(Sz)23]=nx{(Sx)21=(Sy)21=0(Sy)23=(Sz)23=0(Sz)21=i(Sx)23=i
(33.19)

Gdy rozpatrzymy równanie (33.15) dla parametru k'=3, wtedy niektóre elementy macierzy współrzędnych operatora spinu można policzyć z trzeciego z trzech równań, które będziemy rozważać:

i[nx(Sx)31+ny(Sy)31+nz(Sz)31]ψ1+i[nx(Sx)32+ny(Sy)32+nz(Sz)32]ψ2+i[nx(Sx)33+ny(Sy)33+nz(Sz)33]ψ3=
=[ψ1nyψ2nx]
(33.20)

Porównując obie strony tożsamości (33.20), dochodzimy więc do wniosku:

{i[nx(Sx)31+ny(Sy)31+nz(Sz)31]=nyi[nx(Sx)32+ny(Sy)32+nz(Sz)32]=nx(Sx)33=(Sy)33=(Sz)33=0{(Sx)31=(Sz)31=0(Sy)32=(Sz)32=0(Sy)31=i(Sx)32=i
(33.21)

Następnie wyznaczmy współrzędne operatora spinowe Ŝ na podstawie obliczeń dokonanych w punktach (33.17), (33.19) i (33.21):

Ŝx=(00000i0i0)
(33.22)
Ŝy=(00i000i00)
(33.23)
Ŝz=(0i0i00000)
(33.24)

Związki między współrzędnymi macierzy spinowych, które później udowodnimy, przedstawiamy podobnie jak dla zwykłych orbitalnych operatorów momentu pędu:

[Ŝx,Ŝy]=iŜz
(33.25)
[Ŝy,Ŝz]=iŜx
(33.26)
[Ŝz,Ŝx]=iŜy
(33.27)

Udowodnijmy komutator (33.25) na współrzędnych operatora momentu pędu spinowego:

[Ŝx,Ŝy]=ŜxŜyŜyŜx=2(00000i0i0)(00i000i00)2(00i000i00)(00000i0i0)=
=2(000100000)2(010000000)=2(010100000)=i2(0i0i00000)=iŜz
(33.28)

Następnie komutator (33.26) na współrzędnych operatora momentu pędu spinowego:

[Ŝy,Ŝz]=ŜyŜzŜzŜy=2(00i000i00)(0i0i00000)2(0i0i00000)(00i000i00)=
=2(000000010)2(000001000)=2(000001010)=i2(00000i00i)=iŜx
(33.29)

I ostatecznie udowodnijmy komutator (33.27) na współrzędnych operatora momentu pędu spinowego:

[Ŝz,Ŝx]=ŜzŜxŜxŜz=2(0i0i00000)(00000i0i0)2(00000i0i0)(0i0i00000)=
=2(001000000)(000000010)=2(001000010)=i2(00i0000i0)=iŜy
(33.30)

Komutatory (33.25), (33.26) i (33.27) można ogólnie zapisać tak jak w przypadku operatorów momentu pędu orbitalnego (7.13) według schematu:

[Ŝi,Ŝj]=iϵijkŜk
(33.31)

Teraz policzmy kwadrat macierzy operatora spinowego, który jak wiadomo jest sumą kwadratów współrzędnych operatora momentu pędu spinowego zdefiniowanych wedle (33.22)(Ŝx), (33.23)(Ŝy) i (33.24)(Ŝz):

Ŝ2=Ŝx2+Ŝy2+Ŝz2=2(00000i0i0)2+2(00i000i00)2+2(0i0i00000)2=
=2[(00000i0i0)(00000i0i0)+(00i000i00)(00i000i00)+(0i0i00000)(0i0i00000)]=
=2[(000010001)+(100000001)+(100010000)]=2(200020002)=22(100010001)=22E3x3

Również dobrze możemy napisać na podstawie ostatnich obliczeń:

Ŝ2=2S(S+1)E3x3
(33.32)
  • gdzie S=1, która jest liczbą kwantową spinowego operatora momentu pędu.

Wyznaczmy wartości własne operatora zetowego macierzy spinowej:

0=det[(0i0i00000)λ(100010001)]=|λi0iλ000λ|=λ|λiiλ|=
=λ(λ22)=λ(λ)(λ+)λ=λ=0λ=
(33.33)

ze względu na obliczenia (33.33) i tożsamości (33.32) dochodzimy do wniosku, że trzecia składowa spinu dla kwantu spełnia podobny warunek jak dla trzeciej składowej (zetowej) momentu pędu względem liczby kwantowej charakteryzującej wartości własne momentu pędu, czyli:

λ=mS=,0,mS=1,0,1
(33.34)

Liczba kwantowa kwadratu całkowitego momentu pędu jest równa S=1, a magnetyczna spinowa liczba magnetyczna jest liczbą mieszczącą się pomiędzy liczbami S=-1 i S=1, i ta liczba zmienia się co jeden, tak jak według przepisu (33.34), podobnie jest dla magnetycznej orbitalnej liczby kwantowej poprzez liczbę kwantową orbitalnego kwadratu momentu pędu (8.83).

Wektory spinowe, równania własne dla całkowitego momentu pędu kwantu

[edytuj]

Prowadźmy trójwymiarową przestrzeń spinową, którymi wektorami bazy w układzie kartezjańskim będą formalnie zapisane wektory spinowe χ, które jak można udowodnić są do siebie ortonormalne i są to kanoniczne wersory w rozważanej bazie:

χx=(100)
(33.35)
χy=(010)
(33.36)
χz=(001)
(33.37)

Wprowadźmy nowe wektory spinowe tensora rangi drugiej zdefiniowane przy pomocy (33.35), (33.36), a ich definicje są takie by ich norma była równa jeden, należy przy tym pamiętać, że te wektory spinowe rangi drugiej nie muszą być wcale do siebie prostopadłe, w przypadku pierwszym spinor o numerze jeden:

χ1=12(χx+iχy)=12((100)+i(010))=12(1i0)
(33.38)

Spinor o numerze zero piszemy bezpośrednio przy pomocy wektora ortonormalnego (33.37):

χ0=χz=(001)
(33.39)

i ostatecznie przy tych samych spinorach co (33.38), ale inaczej zdefiniowane przy oznaczeniu równym minus jeden:

χ1=12(χxiχy)=12((100)i(010))=12(1i0)
(33.40)

Zbudujmy operator spinowy ze znakiem plus, który jest sumą operatora spinowego iksowego (33.22) i igrekowanego (33.23), pomnożonej przez jednostkę urojoną, całość podzielona przez stałą kreśloną Plancka:

Ŝ+=1(Ŝx+iŜy)=(00000i0i0)+i(00i000i00)=(00100i1i0)
(33.41)

Dalej zbudujmy operator spinowy ze znakiem minus, który jest różnicą operatora spinowego iksowego (33.22) i igrekowanego (33.23), pomnożonej przez jednostkę urojoną, całość podzielona przez stałą kreśloną Plancka:

Ŝ=1(ŜxiŜy)=(00000i0i0)i(00i000i00)=(00100i1i0)
(33.42)

Ostatecznie zbudujmy operator spinowy ze znakiem zero, który jest zetowym operatorem spinowym (33.24) podzielonej przez stałą Plancka.

Ŝ0=1Ŝz=(0i0i00000)
(33.43)

Również można udowodnić, że zachodzą związki na wektorach spinowych drugiego rzędu na macierzach spinowych (33.41), (33.42) i (33.43)

Ŝ+χk=(Sk)(S+k+1)χk+1
(33.44)
Ŝχk=(S+k)(Sk+1)χk1
(33.45)
Ŝ0χk=kχk
(33.46)
  • gdzie S=1, zaś k=+1,0,1

Udowodnijmy zależności (33.44), (33.45) i (33.46) dla trzech możliwych k i jednego S:

Ŝ+χ1=(00100i1i0)12(1i0)=12(002)=2(001)=(1(1))(1+(1))(001)
Ŝ+χ0=(00100i1i0)(001)=(1i0)=(10)(1+0+1)12(1i0)
Ŝχ1=(00100i1i0)12(1i0)=12(002)=2(001)=(1+1)(11+1)(001)


Ŝχ0=(00100i1i0)(001)=(1i0)=(1+0)(10+1)12(1i0)

Ŝ0χ1=(0i0i00000)12(1i0)=12(1i0)=12(1i0)


Ŝ0χ0=(0i0i00000)(001)=(000)=0(001)


Ŝ0χ1=(0i0i00000)12(1i0)=12(1i0)=112(1i0)

Wiadomo, że jeśli zachodzi:l̂zYlm=mYlm, to powinno również zachodzić inne równanie własne, ale wynikającego ze wspomnianego równania:

ĵz(l̂Ylm)=m(l̂Ylm)
(33.47)

A teraz przejdźmy do dowodu ostatniego stwierdzenia wykorzystując przy tym, że całkowity operator momentu pędu promieniowania elektromagnetycznego jest wyrażony wzorem (33.10), to operator całkowitego momentu pędu wyraża się bardzo podobnym wzorem, zatem wspomniany ostatnio wzór można udowodnić wedle przekształceń:

ĵz(l̂Ylm)=(l̂z+Ŝz)(l̂xYlml̂yYlml̂zYlm)=(l̂zl̂xYlml̂zl̂yYlml̂zl̂zYlm)+(010i00000)(l̂xYlml̂yYlml̂zYlm)=
=(l̂xl̂zYlm+il̂yYlml̂yl̂zYlmil̂xYlml̂yl̂zYlm+0)+(il̂yYlmil̂xYlm0)=(l̂xl̂zYlml̂yl̂zYlml̂yl̂zYlm)=m(l̂xYlml̂yYlml̂zYlm)=m(l̂Ylm)
(33.48)
  • gdzie m jest to wartość własna operatora l̂z, a zarazem operatora ĵz

Jeśli zachodzi warunek l̂2Ylm=l(l+1)Ylm, to można udowodnić:

ĵ2(l̂Yjm)=j(j+1)(l̂Yjm)
(33.49)

a zarazem zachodzi j=l, co później udowodnimy intuicyjnie. Aby udowodnić równanie własne (33.49), musimy wiedzieć, że: [l̂,Ŝ]=0, co jest trywialnym dowodem, ze względu, że współrzędne operatora momentu pędu spinowego są macierzami, współrzędne operatora momentu pędu orbitalnego są operatorami, coś w rodzaju różniczkowania, dlatego te operatory są przemienne, zatem przejdźmy do głównego nurtu dowodu wykorzystując przy tym, że operator całkowitego momentu pędu promieniowania elektromagnetycznego jest wyrażony wzorem (33.10):

ĵ2(l̂Ylm)=(l̂+Ŝ)2(l̂Ylm)=l̂2(l̂Ylm)+Ŝ2(l̂Ylm)+2l̂Ŝ(l̂Ylm)=(l(l+1)2+S(S+1)2)(l̂Ylm)+2l̂Ŝ(l̂Ylm)
(33.50)

Obliczmy wyrażenie pomocnicze: l̂Ŝ(l̂Yjm), występujące w dowodzie (33.50), by dalej przeprowadzić ten nas powyższy dowód:

l̂Ŝ(l̂Ylm)=(l̂xŜx+l̂yŜy+l̂zŜz)(l̂Ylm)=[l̂x(00000i0i0)+l̂y(00i000i00)+l̂z(0i0i00000)](l̂Ylm)=
=(0il̂zil̂yil̂z0il̂xil̂yil̂x0)(l̂xl̂yl̂z)Ylm=(i(l̂yl̂zl̂zl̂z)i(l̂zl̂xl̂xl̂z)i(l̂xl̂yl̂yl̂x))Ylm=i2(il̂xil̂yil̂z)Ylm=2l̂Ylm
(33.51)

Teraz powróćmy do dowodu (33.50) na wartość własną całkowitego momentu pędu promieniowania elektromagnetycznego, korzystając przy tym z obliczeń pomocniczych (33.51):

ĵ2(l̂Ylm)=(l(l+1)+S(S+1)2)2(l̂Ylm)=(l(l+1)+22)2(l̂Ylm)=l(l+1)2(l̂Ylm)
(33.52)

Z definicji wartość własna kwadratu momentu pędu całkowitego przyjmuje się jako j(j+1)2, a momentu pędu orbitalnego l(l+1)2. Dochodzimy do wniosku, że w obliczeniach (33.52) musimy przyjąć j=l, ponieważ przy składaniu orbitalnego momentu pędu ze spinem, którego liczba kwantowa jest w trzech postaciach, że j=l-1,l,l+1, ale pole elektromagnetyczne wybiera jedną z nich, tzn. j=l i dlatego musimy napisać dla (33.47) i (33.52) w postaci równań własnych:

ĵz(l̂Yjm)=m(l̂Yjm)
(33.53)
ĵ2(l̂Yjm)=j(j+1)2(l̂Yjm)
(33.54)