Równanie Diraca

Z testwiki
Przejdź do nawigacji Przejdź do wyszukiwania

Równanie Diraca – jedno z fundamentalnych równań w relatywistycznej mechanice kwantowej, sformułowane przez angielskiego fizyka Paula Diraca w 1928 roku[1], słuszne dla cząstek o dowolnie wielkich energiach (tzw. cząstek relatywistycznych) o spinie 1/2 (fermiony, np. elektrony, kwarki), swobodnych i oddziałujących z polem elektromagnetycznym. Istnienie spinu wynika z samego żądania relatywistycznej niezmienniczości równania ruchu cząstek. Odpowiada równaniu Pauliego, które także zawiera spin cząstek, ale wprowadza go w sposób fenomenologiczny, niejako sztuczny, a jedynie dlatego, by otrzymać zgodność z doświadczeniem Sterna-Gerlacha (rozszerzając formalizm nierelatywistycznego równania Schrödingera).

Równanie Diraca jest równaniem macierzowym – de facto stanowi ono układ 4 równań ze względu na fakt, iż symbole gamma (lub alfa, beta), występujące w tym równaniu, są macierzami 4×4.

Równania Diraca zapisuje się w postaci jawnie relatywistycznie niezmienniczej lub w tzw. obrazie Schrödingera. Ta ostatnia postać została najpierw wyprowadzona przez Diraca i jest stosowana ze względu na wygodę do wykonywania obliczeń, gdyż odróżnia współrzędne przestrzenne od współrzędnej czasowej.

Równanie Diraca zostało potwierdzone w odniesieniu do struktury subtelnej widma atomu wodoru, wykazując znakomitą zgodność z pomiarami. Przewiduje istnienie antycząstek. Niektóre jednak efekty, takie jak kreacja i anihilacja cząstek czy przesunięcie Lamba tłumaczy dopiero elektrodynamika kwantowa.

Macierze gamma γμ

Macierze gamma γμ to macierze zespolone 4×4 spełniające 16 reguł antykomutacyjnych w postaci

{γμ,γν}=2gμνI,

gdzie:

μ,ν=0,1,2,3,
{A,B}=AB+BA – tzw. antykomutator,
gμν – elementy tensora metrycznego czasoprzestrzeni gμν=(1000010000100001); np. g00=1,g11=1 itd.,
I=(1000010000100001) – macierz jednostkowa 4×4.

Powyższa reguła określająca macierze gamma wynika m.in. z wymagania, by spełnione było równanie Kleina-Gordona. Jest bardzo wiele sposobów wyboru tych macierzy, np. reprezentacja Pauliego-Diraca ma postać:

γ0=(I00I)
γi=(0σiσi0)
σi, i=1,2,3macierzami Pauliego, zaś I jest tu macierzą jednostkową 2×2.

Jawnie relatywistycznie niezmiennicza postać równania Diraca

Znaczenie jawnie niezmienniczej postaci

Równania Diraca zapisane w postaci jawnie relatywistycznie niezmienniczej to taka postać równania Diraca, która formalnie nie odróżnia czasu od współrzędnych przestrzennych, ale: (1) traktuje czas i współrzędne przestrzenne położenia jako współrzędne czterowektora położenia cząstki w czasoprzestrzeni (2) nie wyróżnia pochodnej po czasie w stosunku do pochodnych po współrzędnych przestrzennych (pochodna po czasie jest elementem czterogradientu, którego pozostałymi trzema elementami są pochodne po współrzędnych przestrzennych). Równanie tak zapisane ma identyczną postać w dowolnym układzie inercjalnym (z jedyną zmianą, że zamiast współrzędnych xν pojawią się współrzędne x'ν, właściwe dla innego układu).

Równanie cząstki swobodnej

W zapisie jawnie relatywistycznie niezmienniczym równanie Diraca dla cząstki swobodnej ma postać

(iγμμmc)Ψ(xν)=0,

gdzie:

μ=0,1,2,3,
xν=(x0=ct,x1,x2,x3) – współrzędne punktu w czasoprzestrzeni,
μ=xμ – element czterogradientu μ=,μ=xμ=(1ct,)=(tc,),
γμmacierze gamma Diraca, tj. γ0,γ1,γ2,γ3,
m – masa cząstki (tzw. masa spoczynkowa),
Ψ(xν)funkcja falowa o 4 składowych zespolonych, tzw. bispinor Diraca,
ijednostka urojona,
stała Plancka podzielona przez 2π,
cprędkość światła.

Równanie cząstki oddziałującej z polem elektromagnetycznym

Jeżeli cząstka nie jest swobodna, ale oddziałuje z zewnętrznym polem elektromagnetycznym, to równanie Diraca przyjmuje postać

(γμ(iμqAμ)mc)Ψ=0,

gdzie:

q – ładunek cząstki,
Aμ=(ϕ0,A)potencjał wektorowy pola zapisany jako czterowektor kowariantny.

Formalnie równanie to można otrzymać z równania Diraca cząstki swobodnej dokonując podstawienia (tzw. reguły Jordana)

iμiμqAμ.

Funkcja falowa Ψ(xν)

Funkcja falowa Ψ(xν), zwana bispinorem Diraca, jest funkcją o 4 składowych zespolonych; zapisuje się ją w postaci kolumny

Ψ(xν)=[ψ1ψ2ψ3ψ4](xν),

przy czym xν𝐑4 oznacza położenie cząstki w czasoprzestrzeni. Nazwa bi-spinor oznacza podwójny spinor. Spinor występuje w równaniu Pauliego, gdzie jest funkcją falową o 2 składnikach, opisujących 2 składowe spinowe (w równaniu Schrödingera funkcja falowa jest 1-składnikowa).

Interpretacja składowych bispinora

Jeżeli pęd jest skierowany w kierunku osi z, to dwie górne składowe bispinora są funkcjami falowymi cząstki:

  • jedna z nich opisuje składową spinu w kierunku zgodnym z wektorem zewnętrznego pola magnetycznego,
  • druga w kierunku przeciwnym.

Dwie dolne składowe odpowiadają analogicznym stanom spinowym antycząstki.

Dla innego skierowania pędu interpretacja taka nie jest jednak właściwaSzablon:Odn.

Bispinor Ψ hermitowsko sprzężony

Definiuje się bispinor Ψ hermitowsko sprzężony do bispinora Ψ – przedstawia on wektor w postaci wiesza, którego elementami są sprzężenia zespolone składowych bispinora (przy czym oznacza sprzężenie hermitowskie)

Ψ=[ψ1*,ψ2*,ψ3*,ψ4*].

Gęstość prawdopodobieństwa w teorii Diraca

Gęstość prawdopodobieństwa definiuje się analogicznie jak w teorii Schrödingera

ρ=ΨΨ.

W definicji gęstości prawdopodobieństwa dla równania Diraca istotna jest kolejność czynników: Ψ musi być przed Ψ, gdyż występuje tu mnożenie wektorów w postaci wiersza i kolumny, i tylko dla takiej kolejności mnożenie da w wyniku skalar. (W analogicznym wyrażeniu na gęstość prawdopodobieństwa dla równania Schrödingera funkcja falowa jest skalarem, stąd kolejność mnożenia nie ma znaczenia).

Wykonując obliczenia otrzymamy

ρ=|ψ1|2+|ψ2|2+|ψ3|2+|ψ4|2.

Wielkość ρ(xν) oznacza, że prawdopodobieństwo znalezienia cząstki materii w położeniu xν𝐑4 jest sumą prawdopodobieństw znalezienia jej w postaci cząstki w stanach spinowych w górę lub w dół, lub w postaci antycząstki w stanach spinowych w górę lub w dół.

Trzeci rodzaj bispinora – bispinor Ψ¯

Prócz bispinorów Ψ oraz Ψ definiuje się bispinor Ψ¯ w postaci

Ψ¯=Ψγ0=[ψ1*,ψ2*,ψ3*,ψ4*].

Powyższy bispinor jest używany do wyrażenia prądów prawdopodobieństwa, odpowiadających relatywistycznie niezmienniczej postaci równania Diraca.

Równanie Diraca w obrazie Schrödingera

Obraz Schrödingera

Równanie Schrödingera ma postać

H^Ψ(r,t)=itΨ(r,t),

gdzie:

H^=p22m

jest operatorem Hamiltona zależnym tylko od współrzędnych przestrzennych, zaś po prawej stronie równania występuje pochodna cząstkowa po czasie.

Dowolne równanie mechaniki kwantowej można zapisać w analogicznej postaci, tj. takiej że z jednej strony równania mamy operator Hamiltona, a z drugiej operator pochodnej czasowej. Taki zapis nazywa się obrazem Schrödingera (lub postacią Schrödingera).

Równanie Diraca w obrazie Schrödingera

Równanie Diraca można przekształcić do postaci w obrazie Schrödingera, wprowadzając macierze alfa i beta

αi=γ0γi,
β=γ0.

Mnożąc obustronnie równanie Diraca podane w postaci jawnie relatywistycznie niezmienniczej przez macierz βγ0, otrzymuje się równanie

(cαp+mc2β)Ψ(r,t)=itΨ(r,t),

gdzie:

c – prędkość światła,
α=(α1,α2,α3) – wektor utworzony z macierzy alfa,
p=(p^1,p^2,p^3) – wektorowy operator pędu,
m – masa cząstki,
Ψ(r,t) – czteroskładnikowa funkcja falowa Diraca.

Operator

H^=[cαp+mc2β]

jest więc operatorem Hamiltona swobodnego, relatywistycznego fermionu o spinie 1/2, analogicznym do operatora Hamiltona cząstki swobodnej w równaniu Schrödingera. W równaniu Diraca operator Hamiltona ma postać operatora macierzowego 4×4, podczas gdy w równaniu Schrödingera wyraża się przez pojedynczy operator (1×1).

Równanie Diraca zapisane w obrazie Schrödingera nie jest jawnie relatywistycznie niezmiennicze, gdyż współrzędna czasowa jest tu wyróżniona. Zapis taki jest jednak wygodny do wykonywania obliczeń w konkretnym układzie odniesienia.

Rozwiązanie równania Diraca dla cząstki swobodnej

Gdy cząstka jest swobodna, to funkcja falowa nie powinna zależeć od współrzędnych, czyli Ψx=Ψy=Ψz=0, co formalnie oznacza, że p=0 i równanie Diraca przyjmuje postaćSzablon:Odn

mc2βΨ(t)=itΨ(t).

Rozwiązania tego równania mają postać

ψA=ei(mc2/)t[ψ1(0)ψ2(0)],
ψB=e+i(mc2/)t[ψ3(0)ψ4(0)].

Pierwsze odpowiada cząstce (np. elektronowi) o energii E=mc2, drugie antycząstce (np. pozytonowi) także o energii E=mc2Szablon:Odn.

Równanie Diraca dla cząstki w polu elektromagnetycznym

Jeżeli cząstka ma ładunek q i oddziałuje z polem elektromagnetycznym o potencjale skalarnym ϕ i potencjale wektorowym A, to operator Hamiltona w równaniu Diraca, zapisanym w obrazie Schrödingera, otrzymuje się, stosując podstawiania (tzw. reguły Jordana)

EEqϕ,
ppqA.

Operator Hamiltona przyjmuje postać

H^=[cα(pqA)+mc2β+qϕ].

Pole traktuje się tu jako klasyczne pole Maxwella, tj. nie poddane tzw. procesowi drugiego kwantowania. Oznacza to, że nie uwzględnia się tu faktu, iż pole elektromagnetyczne występuje de facto w postaci kwantów energii, fotonów. Pominięcie tego jest uzasadnione wtedy, gdy pole ma dużą energię wobec energii cząstki.

Operator spinu

Pokażemy, że operator spinu wynika w naturalny sposób z równania Diraca, tj. z samego faktu, iż równanie to ma postać relatywistycznie niezmienniczą. M.in. z tej racji równanie Diraca stanowi „klejnot fizyki”. (Dla porównania: Pauli wprowadził operator spinu w sposób fenomenologiczny, tj. zmodyfikował jedynie równania Schrödingera tak, by uzyskać zgodność opisu z wynikami eksperymentów, gdzie ujawnia się spin cząstek).

Macierze sigma Diraca σiD

σiD=(σi00σi),
σi, i=x,y,z lub i=1,2,3 macierzami Pauliego, zaś 0 jest macierzą zerową 2×2.

Macierz te mają wymiar 4×4. Przy czym zachodzą równości

σzD=i2[αx,αy],
σyD=i2[αz,αx],
σxD=i2[αy,αz].

Komutator hamiltonianu z operatorem momentu pędu

Obliczamy komutator hamiltonianu H^ cząstki swobodnej z operatorem momentu pędu L^, gdzie:

H^=cαp+mc2β.

Np. dla składowej L^x

L^x=yp^zzp^y

otrzymamy

[L^x,H^]=cα{(yp^zzp^y)pp(yp^zzp^y)}=ic(αyp^zαzp^y)0.

Oznacza to, że moment pędu nie komutuje z hamiltonianem, nie jest więc zachowany (nie jest stałą ruchu).

Definicja operatora spinu Diraca

Wektorowy operator spinu S=[S^x,S^y,S^z] definiuje się, żądając (1) operator całkowitego momentu pędu cząstki J (tj. suma operatora spinu Diraca i operatora orbitalnego momentu pędu L) musi komutować z hamiltonianem równania Diraca dla cząstki swobodnej (jeżeli bowiem cząstka jest swobodna, to jej całkowity moment pędu musi być zachowany) (2) operator spinu Diraca musi spełniać odpowiednie reguły komutacyjne (dokładniej warunek ten omówiono niżej – patrz sekcja „Komutatory operatorów spinu S^x,S^y,S^z”).

Mamy więc

J=L+S

oraz

(1) [J,H^]=0,
(2) [S^k,S^l]=iϵklmS^m, k,l,m=x,y,z,

gdzie ϵklmtensor zupełnie antysymetryczny.

Warunki (1) i (2) są spełnione, jeżeli składowe operatora spinu S=[S^x,S^y,S^z] mają postać

S^i=12σiD, i=x,y,z,

czyli:

S^x=12σxD=12(0001001001001000),
S^y=12σyD=12(000i00i00i00i000),
S^z=12σzD=12(1000010000100001).

Składowe operator spinu są więc w reprezentacji macierzowej macierzami 4×4, w odróżnieniu od składowych operatora spinu Pauliego, które są macierzami 2×2

S^i=12σi, i=x,y,z.

Np.

S^z=12σz=12(1001).

Sens fizyczny operatora spinu Diraca

Sens fizyczny każdej ze składowych operatora spinu jest analogiczny. Np. operator spinu Diraca Sz odpowiada pomiarowi składowej z-owej spinu cząstki – zgodnej z kierunkiem osi z lub przeciwnej do kierunku tej osi, oraz pomiarowi składowej spinu antycząstki zgodnej i przeciwnej do osi z. (Dla porównania, operator spinu Pauliego Sz odpowiada tylko pomiarowi składowej z-owej spinu cząstki; równanie Pauliego nie przewiduje bowiem istnienia antycząstek.)

Kwadrat operatora spinu S2

Kwadrat operatora spinu Diraca ma postać:

S2=SS=S^x2+S^y2+S^z2.

Podstawiając wyrażenia na operatory S^x,S^y,S^z, otrzymuje się:

S2=324(1001)=324I,

gdzie 1,0 są macierzami 2×2, odpowiednio jednostkową i zerową, zaś I – macierz jednostkowa 4×4.

Pierwiastek ze średniej wartości S2 operatora S2 określa wartość mierzonego spinu, przy czym

S2=ΨS2Ψ=324|Ψ|2.

Ponieważ wektor stanu jest z założenia unormowany, to |Ψ|2=1. Stąd:

S2=32.

Powyższy wynik jest zgodny z ogólnym wzorem na długość wektora momentu pędu o liczbie spinowej s

S2=s(s+1),

przy czym dla s=12 otrzymuje się wcześniej podany wynik.

Tak więc pomiar spinu na cząstce Diraca daje zawsze wartość spinu S2=32, przy czym mierzy się spin cząstki albo antycząstki.

Komutatory operatorów spinu S^x,S^y,S^z

Z pomiarów wynika, że jest możliwe zmierzenie tylko jednej spośród trzech składowych wektora spinu. Z tej racji na operatory spinu S^x,S^y,S^z nakłada się reguły komutacyjne identyczne jak reguły komutacyjne operatorów momentu pędu L^x,L^y,L^z czy operatorów spinu Pauliego:

[S^x,S^y]=iS^z,
[S^z,S^x]=iS^y,
[S^y,S^z]=iS^x.

Operatory te nie komutują ze sobą (tzn. komutatory są 0), co odpowiada faktom eksperymentalnym, iż jest możliwe jednoczesne zmierzenie tylko jednej ze składowych spinu.

Komutatory operatorów S2 oraz S^x,S^y,S^z

Z pomiarów wynika, że jest możliwe zmierzenie jednoczesne jednej spośród trzech składowych wektora spinu oraz całkowitej wartości spinu. Z tej racji na operatory spinu S^x,S^y,S^z muszą komutować z operatorem S2.

Podane wyżej operatory spełniają te reguły, gdyż operator S2 wyraża się przez macierz jednostkową, a w związku z tym komutuje z dowolną ze składowych spinu, np.

[S2,S^x]=0.

Komutatory operatorów S,L oraz hamiltonianu

„Stożki wektorowe” momentów pędu: całkowitego J (fiolet), orbitalnego L (niebieski) i spinowego S (zielony). Stożki powstają na skutek nieoznaczoności kwantowej składowych tych momentów

(1) Operatory S,L komutują ze sobą, tj.

[S,L]=0,

co oznacza, że jest możliwe zmierzenie jednoczesne wartości spinu oraz momentu pędu (operatory te działają w innych przestrzeniach Hilberta).

(2) Operatory S,L nie komutują z osobna z operatorem Hamiltona cząstki swobodnej

[S,H^]0,
[L,H^]0,

ale suma tych operatorów J=L+S komutuje, tj.

[J,H^]=0.

Oznacza to, że moment pędu orbitalny i spinowy cząstki swobodnej mogą zmieniać się w czasie, ale tak, że ich suma jest stała, przy czym każdy z wektorów z osobna może przyjąć w miarę dowolne położenie w przestrzeni – wektory te osobno nie są zachowane, bo nie komutują z hamiltonianem. Pokazane na rysunku stożki wektorowe uwidaczniają dobrze tę zależność: jeżeli wektor momentu pędu wykonuje precesję po stożku niebieskim, to wektor spinu musi odpowiednio zmienić swoje położenie na stożku zielonym tak, by sumaryczny wektor pozostał na stożku fioletowym.

Prawdopodobieństwa pomiaru spinu S^x,S^y,S^z

Aby obliczyć prawdopodobieństwo otrzymania w eksperymencie np. antycząstki ze spinem skierowanym w kierunku osi z rozkłada się bi-spinor Diraca (o postaci takiej, że odpowiada stanowi cząstki) w bazie wektorów własnych |1,|2,|3,|4 operatora spinu S^z, gdzie

|1=[1000],|2=[0100],|3=[0010],|4=[0001]

(tzw. notacja Diraca), przy czym:

|1 – wektor własny operatora S^z odpowiadający pomiarowi spinu cząstki w kierunku +z,
|2 – wektor własny operatora S^z odpowiadający pomiarowi spinu cząstki w kierunku z,
|3 – wektor własny operatora S^z odpowiadający pomiarowi spinu antycząstki w kierunku +z itd.

Wtedy

Ψ=ψ1|1+ψ2|2+ψ3|3+ψ4|4, oraz np.
|ψ3|2 – prawdopodobieństwo otrzymania wartości rzutu spinu w kierunku +z dla antycząstki itp.

Analogicznie oblicza się prawdopodobieństwa uzyskania rzutów spinu przy pomiarze w kierunkach x oraz y (przy czym teraz trzeba rozłożyć bi-spinor Diraca w bazach wektorów własnych operatorów S^x,S^y).

Średnia wartość pomiaru spinu S^z

Średnią wartość pomiaru spinu S^z na cząstce opisanej stanem Ψ oblicza się ze wzoru

S^z=ΨS^zΨ=[ψ1*,ψ2*,ψ3*,ψ4*]12(1000010000100001)[ψ1ψ2ψ3ψ4]=12(|ψ1|2|ψ2|2+|ψ3|2|ψ4|2),

przy czym minusy odpowiadają skierowaniu spinu cząstki i antycząstki w kierunku z, a plusy w kierunku +z.

Operator spinu cząstki w polu elektromagnetycznym centralnym

Jeżeli cząstka naładowana znajduje się w polu elektromagnetycznym centralnym (jak jest np. w przypadku atomu wodoru), to operator spinu jest identyczny z operatorem spinu cząstki swobodnej, tj. S=[S^x,S^y,S^z]

oraz

S^i=12σiD, i=x,y,z.

W polu centralnym bowiem całkowity moment pędu cząstki jest stały (jest to analogiczne do prawa zachowania momentu pędu w polu centralnym, znanym z fizyki klasycznej).

Fermiony Majorany

Cząstki spełniające równanie Diraca są fermionami. Jednak teoretycznie mogą istnieć inne fermiony, które nie spełniają równania Diraca – są to tzw. cząstki Majorany.

Lagranżjan Diraca

Równanie Diraca i sprzężone równanie Diraca można otrzymać dokonując wariacji działania

𝒮=t1t2dtd3x

w której gęstość lagranżjanu dana jest wzorem

=icψγμμψmc2ψψ

Wariując działanie względem ψ otrzyma się równanie Diraca. Wariując działanie względem ψ otrzyma się sprzężone równanie Diraca.

Zobacz też

Przypisy

Szablon:Przypisy

Bibliografia

Linki zewnętrzne

Szablon:Szablon nawigacyjny Szablon:Równania różniczkowe

Szablon:Kontrola autorytatywna