Kwantowy oscylator harmoniczny

Kwantowy oscylator harmoniczny – układ fizyczny rozmiarów atomowych lub subatomowych (np. jon w sieci krystalicznej lub w cząsteczka gazu) wykonujący ruch drgający (oscylacyjny) pod wpływem siły proporcjonalnej do wychylenia od położenia równowagi. Właściwy opis ruchu wymaga zastosowania mechaniki kwantowej, co sprowadza się do znalezienia rozwiązań równania Schrödingera. Dowodem eksperymentalnym konieczności zastosowania mechaniki kwantowej do opisu właściwości mikroskopowych układów drgających jest np. nieciągłe widmo promieniowania emitowane przez drgające cząsteczki. Makroskopowym odpowiednikiem oscylatora kwantowego jest klasyczny oscylator harmoniczny, którym jest ciało makroskopowe o stosunkowo dużej masie, zawieszone np. na sprężynie i wykonujące drgania; do opisu jego ruchu wystarczająca jest mechanika klasyczna. Pojęcie oscylatora ma duże zastosowanie i znaczenie w wielu działach fizyki klasycznej i kwantowej.

Cząsteczka HCl jako oscylator kwantowy drgający na poziomie energii $E_{3}.$ Energia jest skwantowana, tzn. może przyjmować tylko skokowe wartości $E_{0},E_{1},\dots ,D_{0}$ jest energią dysocjacji, $r_{0}$ średnią odległością atomów, U energią potencjalną ich ruchu oscylacyjnego. Atom wodoru umieszczono w początku układu współrzędnych, aby pokazać zmiany średniej odległości atomów na krzywej.

Znaczenie oscylatora harmonicznego

Teoria oscylatora harmonicznego ma w fizyce doniosłe znaczenie. Jest tak dlatego, że wiele różnych układów fizycznych jest opisywanych równaniami o postaci identycznej z postacią pojedynczego oscylatora lub zespołu oscylatorów harmonicznych, które są słabo ze sobą sprzężone (czyli słabo ze sobą oddziałują). W pierwszym przybliżeniu zaniedbuje się oddziaływanie – wtedy układ jest matematycznie równoważny prostemu do opisu układowi niezależnie drgających oscylatorów harmonicznych, np. wszystkie cząstki wieloatomowe wykonują drgania, które z dobrym przybliżeniem można opisywać w ramach teorii oscylatorów harmonicznych.

Początki fizyki kwantowej wiążą się z pojęciem oscylatora kwantowego. Mianowicie Max Planck, próbując wyjaśnić zjawisko promieniowania termicznego ciał (tzw. zjawisko promieniowania ciała doskonale czarnego), przyjął, że cząstki materii emitujące i absorbujące promieniowanie zachowują się jak oscylatory. Planck założył, że energia oscylatora nie może być dowolna, lecz jest skwantowana, to znaczy może przyjmować tylko ściśle określone wartości. Założenie to nie posiadało wtedy żadnego uzasadnienia w znanych teoriach fizycznych.

Ważnym osiągnięciem było opisanie pola elektromagnetycznego jako pola kwantowego (tzw. drugie kwantowanie). W początkach rozwoju teorii kwantowej pole elektromagnetyczne traktowano jako pole klasyczne, będące źródłem potencjału $\phi ({\boldsymbol {r}},t)$ oddziałującego na cząstki naładowane, wstawianym np. do równania Schrödingera (nierelatywistycznego), a nawet do relatywistycznego równania Diraca. Jednak dokładny opis pól kwantowych wymaga potraktowania ich jako układów kwantowych, gdzie równania pola są analogiczne do równań oscylatorów kwantowych. Ze względu na to, że pola kwantowe mają charakter bozonowy lub fermionowy, odpowiadające im oscylatory określa się jako oscylatory bozonowe i fermionowe (patrz niżej).

Klasyczny oscylator harmoniczny

Klasyczny oscylator harmoniczny – to ciało o masie $m,$ na które działa siła proporcjonalna do wychylenia $x$ ciała od stanu równowagi i mająca przeciwny zwrot

$F=-kx,$

gdzie $k$ jest stałą wielkością (tzw. stałą sprężystości). Przykładem oscylatora harmonicznego jest ciało na sprężynie, wykonujące niewielkie drgania od położenia równowagi, co zapewnia słuszność założenia o proporcjonalności siły do wychylenia (dla dużych wychyleń założenie to nie byłoby słuszne). Układ drgający ma energię potencjalną:

$U(x)={\frac {1}{2}}kx^{2}={\frac {1}{2}}m\omega ^{2}x^{2},$

która jest tym większa, im większe jest rozciągnięcie sprężyny ($\omega ={\sqrt {k/m}}$ jest częstotliwością kołową ruchu drgającego). Energia całkowita układu jest sumą energii kinetycznej i potencjalnej

$E={\frac {p^{2}}{2m}}+U(x),$

gdzie $p=mv$ oznacza pęd ciała drgającego, zajmujące położenie $x.$ Całkowita energia układu drgającego harmonicznie nie ulega zmianie w czasie, mimo że energia potencjalna zamienia się cyklicznie w energię kinetyczną i odwrotnie, kinetyczna przechodzi w potencjalną.

Kwantowy oscylator harmoniczny – przypadek stałej energii drgań

Porównanie ruchu oscylatora harmonicznego klasycznego i kwantowego. (A–B) Oscylator klasyczny – to cząstka masywna (reprezentowana przez kulkę na sprężynie), której ruch można dobrze opisać za pomocą równań Newtona. (C–H) Oscylator kwantowy – to cząstka mikroskopowa, której ruch można poprawnie opisać jedynie za pomocą równania Schrödingera. Na rys. C–H pokazano niektóre w wielu możliwych funkcji falowych $\Psi (x,t),$ stanowiących rozwiązania równania Schrödingera. Na osi pionowej odłożono część rzeczywistą (kolor niebieski) oraz część urojoną (kolor czerwony) funkcji falowej, uzależnione od położenia cząstki $x,$ odłożonego na osi poziomej. Rys. C,D,E,F, ale nie G,H przedstawiają stany stacjonarne (o stałej energii). H przedstawia stan koherentny – stan kwantowy, który przybliża ruch oscylatora klasycznego.

W mechanice kwantowej do opisu ruchu układów fizycznych stosuje się zamiast równania Newtona równanie Schrödingera. Konkretna jego postać zależy od opisywanej sytuacji fizycznej. Jedną z metod znalezienia postaci równania Schrödingera jest tzw. metoda kwantowania, polegająca na zamianie w równaniach ruchu mechaniki klasycznej wektora pędu ciała $p$ na operator pędu ${\hat {p}}.$ Współrzędne położenia ciała, np. $x,$ pozostawia się przy tym bez zmian (formalnie każdej współrzędnej nadaje się nazwę operatora położenia). Słuszność tej metody uzasadnia fakt, że otrzymane za jej pomocą równania dają przewidywania zgodne z wynikami eksperymentów. W przypadku ruchu jednowymiarowego operator pędu ma tylko jedną składową i ma postać:

${\hat {p}}=-\mathrm {i} \hbar {\frac {\mathrm {d} }{\mathrm {d} x}}.$

Ponieważ poszukiwany jest opis stanu układu w zależności od współrzędnych $x,$ dlatego trzeba znaleźć jawną postać równania Schrödingera w reprezentacji położeniowej, przy czym dla uproszczenia założymy, że energia układu jest niezmienna. (Podobnie zakłada się, rozwiązując zagadnienie poziomów energetycznych atomu wodoru). Jest to uzasadnione, jeżeli układ drgający pozostaje dłuższy czas w izolacji od otoczenia. Dlatego stosuje się równanie Schrödingera niezależne od czasu:

${\hat {H}}\psi (x)=E\ \psi (x).$

gdzie $E$ – energia układu, $\psi (x)$ – składnik funkcji falowej niezależny od czasu, opisujący stan układu, ${\hat {H}}$ – operator Hamiltona. Postać operatora ${\hat {H}}$ otrzymuje się, zastępując w wyrażeniu na energię całkowitą $E$ oscylatora klasycznego (patrz wyżej), pęd $p$ operatorem pędu ${\hat {p}}{:}$

${\hat {H}}={\frac {{\hat {p}}^{2}}{2m}}+{\frac {1}{2}}m\omega ^{2}x^{2}.$

Podstawiając jawną postać operatora wartości pędu, otrzymuje się ostatecznie:

${\hat {H}}={\frac {-\hbar ^{2}}{2m}}{\frac {\mathrm {d} ^{2}}{\mathrm {d} x^{2}}}+{\frac {1}{2}}m\omega ^{2}x^{2}.$

Równanie Schrödingera bez czasu przyjmuje więc postać:

${\frac {-\hbar ^{2}}{2m}}{\frac {d^{2}}{dx^{2}}}\psi (x)+{\frac {1}{2}}m\omega ^{2}x^{2}\psi (x)=E\psi (x).$

Rozwiązanie tego równania daje zbiór możliwych stanów stacjonarnych

$\psi _{n}(x)={\frac {1}{\sqrt {2^{n}\,n!}}}\left({\frac {m\omega }{\pi \hbar }}\right)^{1/4}\mathrm {e} ^{-{\frac {m\omega x^{2}}{2\hbar }}}H_{n}\left({\sqrt {\frac {m\omega }{\hbar }}}x\right),\quad n=0,1,2,\dots$

Energia potencjalna oscylatora (niebieska parabola) i kilka pierwszych stanów własnych $\psi _{n}(x)$ w zależności od wychylenia $x.$ Wykresy umieszczono na wysokościach odpowiadających energiom $E_{n}$ tych stanów. Zaznaczono skok wartości energii $\hbar \omega .$ Energia stanu podstawowego jest większa od zera. Widać też, że amplituda ruchu wzrasta z energią.

Gęstości prawdopodobieństwa $P_{n}(x)=|\psi _{n}(x)|^{2}$ odpowiadające funkcjom falowym $\psi _{n}(x).$ Z wykresów widać, że szansa znalezienia oscylatora w pobliżu maksymalnego wychylenia $x$ od położenia równowagi jest największa. Odpowiada to sytuacji oscylatora klasycznego, który w punktach amplitudy przebywa najdłużej. U dołu – gęstość prawdopodobieństwa dla stanu podstawowego: widać, że w tym stanie oscylator najczęściej znajduje się w położeniu równowagi. Zaznaczono wartości energii $E_{n}(x)$ odpowiadające rozkładom $P_{n}(x).$

gdzie:

$H_{n}(x)=(-1)^{n}\mathrm {e} ^{x^{2}}{\frac {\mathrm {d} ^{n}}{\mathrm {d} x^{n}}}\mathrm {e} ^{-x^{2}}$

oznaczają wielomiany Hermite’a, gdzie np.

${\begin{aligned}H_{0}(x)&=1,\\H_{1}(x)&=2x,\\H_{2}(x)&=(2x)^{2}-2=4x^{2}-2,\\H_{3}(x)&=(2x)^{3}-6(2x)=8x^{3}-12x,\\H_{4}(x)&=(2x)^{4}-12(2x)^{2}+12=16x^{4}-48x^{2}+12,\\\ldots \end{aligned}}$

Postać funkcji falowych $\psi _{0}(x),\dots \psi _{7}(x)$ pokazuje rysunek obok.

Stanom $\psi _{n}(x)$ odpowiadają energie oscylatora:

$E_{n}=\hbar \omega \left(n+{\frac {1}{2}}\right),\quad n=0,1,2\ldots$

Układ kwantowy drgający harmonicznie przyjmuje tylko wyróżnione wartości energii, czym różni się od układu klasycznego (makroskopowego) – ten ostatni może drgać, mając dowolną wartość energii. Ponieważ drgające układy mikroskopowe faktycznie przyjmują dyskretne poziomy energii, widoczne się staje, że teoria Schrödingera dostarcza właściwego ich opisu.

Różnica między kolejnymi poziomami energii jest stała i wynosi $\hbar \omega .$ Animacja u góry pokazuje, że poziomy rzeczywistej cząsteczki HCl dla większych energii $E_{n}$ stopniowo zagęszczają się. Dla większych energii wzrasta amplituda i drgania przestają być harmoniczne, co było zakładane wcześniej. Opis takiego ruchu wymagałby dodania do hamiltonianu dodatkowego wyrazu, odpowiadającego za nieharmoniczny składnik siły. Drugie spostrzeżenie: najmniejsza energia drgań nie jest zerowa, gdyż $E_{0}=1/2\hbar \omega .$ Jest to tzw. energia drgań zerowych, która nie jest znana fizyce klasycznej. Istnienie tej energii oznacza, że układ kwantowy nigdy nie może być w absolutnym spoczynku.

Fundamentalne relacje komutacyjne

Można pokazać, że powyżej zdefiniowane operatory położenia i pędu spełniają fundamentalną relację komutacyjną

$[{\hat {x}},{\hat {p}}]=\mathrm {i} \hbar$

gdzie $[{\hat {x}},{\hat {p}}]={\hat {x}}{\hat {p}}-{\hat {p}}{\hat {x}}$ – tzw. komutator operatorów, ${\hat {x}}\equiv x$ – operator położenia.

Powyższa relacja komutacyjna oraz analogiczne reguły dla operatorów położenia i pędu w kierunkach y oraz z stanowią fundamentalne warunki, na których zbudowana jest mechanika kwantowa. Relacje te zostały najpierw odkryte przez Heisneberga w 1925 r. i stanowią podstawę sformułowania mechaniki kwantowej w wersji zwanej mechaniką macierzową. Schrödinger niezależnie doszedł do własnego sformułowania mechaniki kwantowej w 1926 r., przyjmując powyżej zdefiniowane operatory położenia i pędu, które dają identyczny warunek komutacyjny. Później pokazano, iż sformułowania Heisenberga i Schrödingera są równoważne.

Bozonowy oscylator harmoniczny

Rozwiązanie równania Schrödingera oscylatora metodą bezpośrednią jest bardzo złożone. Powyżej został podany jedynie wynik. Jednak można uprościć poszukiwanie rozwiązania, stosując tzw. metodę algebraiczną. Metoda ta polega na zastąpieniu operatorów $x,{\hat {p}}$ operatorami anihilacji ${\hat {a}}$ oraz kreacji ${\hat {a}}^{\dagger }.$

Operatory kreacji i anihilacji

Operatory anihilacji i kreacji definiuje się następująco:

${\hat {a}}={\frac {1}{\sqrt {2}}}\left({\sqrt {\frac {m\omega }{\hbar }}}x+\mathrm {i} {\frac {1}{\sqrt {m\omega \hbar }}}{\hat {p}}\right),$

${\hat {a}}^{\dagger }={\frac {1}{\sqrt {2}}}\left({\sqrt {\frac {m\omega }{\hbar }}}x-\mathrm {i} {\frac {1}{\sqrt {m\omega \hbar }}}{\hat {p}}\right).$

Operatory położenia $x$ i pędu ${\hat {p}}$ wyrażone przez te operatory mają postać:

$x={\sqrt {\frac {\hbar }{2m\omega }}}({\hat {a}}+{\hat {a}}^{\dagger }),$

${\hat {p}}=\mathrm {i} {\sqrt {\frac {m\omega \hbar }{2}}}({\hat {a}}^{\dagger }-{\hat {a}}).$

Użyteczność metody algebraicznej bierze się stąd, że operatory ${\hat {a}}$ oraz ${\hat {a}}^{\dagger }$ dają proste reguły komutacjne:

$[{\hat {a}},{\hat {a}}^{\dagger }]=1,$

$[{\hat {N}},{\hat {a}}]=-{\hat {a}},$

$[{\hat {N}},{\hat {a}}^{\dagger }]={\hat {a}}^{\dagger },$

gdzie ${\hat {N}}={\hat {a}}^{\dagger }{\hat {a}}$ to tzw. operator liczby cząstek, przy czym przez komutator dowolnych operatorów ${\hat {A}},{\hat {B}}$ rozumie się wyrażenie $[{\hat {A}},{\hat {B}}]={\hat {A}}{\hat {B}}-{\hat {B}}{\hat {A}}.$

Dzięki temu złożone przekształcenia zostają zastąpione prostszymi manipulacjami na symbolach. (Powyższe reguły komutacyjne są konsekwencją powyżej podanej fundamentalnej reguły komutayjnej pomiędzy operatorami położenia i pędu).

Operator Hamiltona wyrażony przez te operatory przyjmuje postać:

${\hat {H}}={\frac {1}{2}}\hbar \omega ({\hat {a}}^{\dagger }{\hat {a}}+{\hat {a}}{\hat {a}}^{\dagger })=\hbar \omega \left({\hat {a}}^{\dagger }{\hat {a}}+{\frac {1}{2}}\right),$

przy czym ostatni wzór uzyskuje się, wykorzystując własność $[{\hat {a}},{\hat {a}}^{\dagger }]=1.$ Niech $|n\rangle \equiv \psi _{n}$ (w zapisie Diraca) oznacza stan własny oscylatora o energii $E_{n}.$ Ponieważ szukane są stany stacjonarne oscylatora, to należy rozwiązać równanie Schrödingera bez czasu:

${\hat {H}}|n\rangle =\hbar \omega \left({\hat {a}}^{\dagger }{\hat {a}}+{\frac {1}{2}}\right)|n\rangle =E_{n}|n\rangle .$

Aby to zrobić, zostanie najpierw pokazane, jak operatory kreacji i anihilacji działają na stany oscylatora.

Działanie operatorów kreacji i anihilacji na stany własne oscylatora

Mnożąc powyższe równanie z lewej strony przez ${\hat {a}}^{\dagger },$ otrzymuje się:

$\hbar \omega \left({\hat {a}}^{\dagger }{\hat {a}}^{\dagger }{\hat {a}}+{\frac {1}{2}}{\hat {a}}^{\dagger }\right)|n\rangle =E_{n}{\hat {a}}^{\dagger }|n\rangle .$

Korzystając z komutatora $[{\hat {a}},{\hat {a}}^{\dagger }]=1,$ dostaniemy:

$\hbar \omega \left(({\hat {a}}^{\dagger }{\hat {a}}{\hat {a}}^{\dagger }-{\hat {a}}^{\dagger })+{\frac {1}{2}}{\hat {a}}^{\dagger }\right)|n\rangle =E_{n}{\hat {a}}^{\dagger }|n\rangle$

i stąd:

$\hbar \omega \left({\hat {a}}^{\dagger }{\hat {a}}+{\frac {1}{2}}\right)({\hat {a}}^{\dagger }|n\rangle )=(E_{n}+\hbar \omega )({\hat {a}}^{\dagger }|n\rangle ),$

czyli

${\hat {H}}\,({\hat {a}}^{\dagger }|n\rangle )=(E_{n}+\hbar \omega )({\hat {a}}^{\dagger }|n\rangle ).$

Wynika stąd, że stany $({\hat {a}}^{\dagger }|n\rangle )$ są stanami własnymi operatora Hamiltona, którym odpowiadają wartości własne $(E_{n}+\hbar \omega ).$ Inaczej patrząc na ten wynik, można powiedzieć, że operator ${\hat {a}}^{\dagger },$ działając na dowolny stan $|n\rangle ,$ tworzy stan o energii powiększonej o kwant $\hbar \omega$ względem energii $E_{n}$ stanu, na który działa. Stąd nazwa tego operatora – operator kreacji.

Podobnie, mnożąc równanie Schrödingera przez operator anihilacji ${\hat {a}},$ otrzymamy

${\hat {H}}\,({\hat {a}}|n\rangle )=(E_{n}-\hbar \omega )({\hat {a}}|n\rangle ),$

co oznacza, że operator anihilacji ${\hat {a}},$ działając na dowolny stan $|n\rangle ,$ tworzy stan o energii mniejszej o $\hbar \omega$ od energii $E_{n}$ stanu, na który działa; stanowi $a|n\rangle$ odpowiada bowiem energia $(E_{n}-\hbar \omega ).$

Dokładne obliczenia pokazują, że działanie operatorów kreacji i anihilacji na stany własne jest następujące:

${\begin{aligned}{\hat {a}}\left|n\right\rangle &={\sqrt {n}}\,\left|n-1\right\rangle ,\\{\hat {a}}^{\dagger }\left|n\right\rangle &={\sqrt {n+1}}\,\left|n+1\right\rangle .\end{aligned}}$

Stan zerowy oscylatora i energie własne

Aby znaleźć najniższy możliwy stan oscylatora zauważmy, że operator anihilacji, działając wielokrotnie na dany stan wyjściowy, będzie tworzył stany o coraz mniejszej energii. Ponieważ energia oscylacji nie może być mniejsza od zera, więc trzeba przyjąć, że istnieje stan najniższy, $|0\rangle ,$ taki że działanie operatorem anihilacji na ten stan daje zero:

${\hat {a}}|0\rangle =0,$

przy czym jeżeli w obliczeniach stan $|0\rangle$ zostanie wyzerowany, to działając następnie jakimkolwiek operatorem, otrzyma się nadal zero, czyli:

${\hat {H}}\,({\hat {a}}|0\rangle )=(E_{0}-\hbar \omega )({\hat {a}}|0\rangle )=0.$

Działając operatorem Hamiltona na stan zerowy, otrzyma się:

${\hat {H}}|0\rangle =\hbar \omega \left({\hat {a}}^{\dagger }{\hat {a}}+{\frac {1}{2}}\right)|0\rangle ={\frac {1}{2}}\hbar \omega {\hat {a}}^{\dagger }{\hat {a}}|0\rangle +{\frac {1}{2}}\hbar \omega |0\rangle =\ldots ={\frac {1}{2}}\hbar \omega |0\rangle ,$

co oznacza, że energia stanu zerowego wynosi $E_{0}={\frac {1}{2}}\hbar \omega .$

Ponieważ $E_{n+1}=E_{n}+\hbar \omega ,$ to otrzymuje się wartości dowolnych energii własnych:

$E_{n}=\hbar \omega \left(n+{\frac {1}{2}}\right),\quad n=0,1,\dots$

Z powyższego wzoru widać, że oscylator bozonowy może przyjmować dowolną energię. Ilość kwantów energii w danym stanie bozonowym nie jest więc niczym ograniczona.

Wyrażenie stanów własnych za pomocą operatora kreacji

Stany własne $|n\rangle$ można wyrazić za pomocą operatora kreacji ${\hat {a}}^{\dagger }$

$|n\rangle ={\frac {1}{\sqrt {n!}}}({\hat {a}}^{\dagger })^{n}|0\rangle .$

Dowód:

${\begin{aligned}\langle n|{\hat {a}}{\hat {a}}^{\dagger }|n\rangle &=\langle n|\left([{\hat {a}},{\hat {a}}^{\dagger }]+{\hat {a}}^{\dagger }{\hat {a}}\right)|n\rangle =\langle n|\left({\hat {N}}+1\right)|n\rangle =n+1\\\Rightarrow {\hat {a}}^{\dagger }|n\rangle &={\sqrt {n+1}}|n+1\rangle \\\Rightarrow |n\rangle &={\frac {{\hat {a}}^{\dagger }}{\sqrt {n}}}|n-1\rangle ={\frac {\left({\hat {a}}^{\dagger }\right)^{2}}{\sqrt {n(n-1)}}}|n-2\rangle =\ldots ={\frac {\left({\hat {a}}^{\dagger }\right)^{n}}{\sqrt {n!}}}|0\rangle .\end{aligned}}$

Funkcja własna stanu zerowego

Postać stanu próżni ma fundamentalne znaczenie, gdyż dopiero znając ten stan, można dokonać obliczeń innych stanów.

Postać stanu $|0\rangle$ w reprezentacji położeniowej wyznacza się, przedstawiając operator anihilacji w jawnej postaci, tj. podstawiając ${\hat {p}}=-\mathrm {i} \hbar {\frac {\mathrm {d} }{\mathrm {d} x}}$ do wyrażenia na ten operator:

${\hat {a}}={\frac {1}{\sqrt {2}}}\left({\sqrt {\frac {m\omega }{\hbar }}}x+{\sqrt {\frac {\hbar }{m\omega }}}{\frac {\mathrm {d} }{\mathrm {d} x}}\right).$

Podstawiając stałą pomocniczą $\xi ={\sqrt {\frac {\hbar }{m\omega }}},$ operator anihilacji przyjmie postać

${\hat {a}}={\frac {1}{\sqrt {2}}}\left({\frac {1}{\xi }}x+\xi {\frac {\mathrm {d} }{\mathrm {d} x}}\right).$

Stan $|0\rangle$ wyrażony w bazie położeniowej jest pewną funkcją zmiennej $x,$ tj. pewną funkcją $\psi _{0}(x)\equiv \langle x|0\rangle ,$ ponieważ operator anihilacji, działając na stan $|0\rangle ,$ ma go zerować, to musi być spełnione równanie

${\hat {a}}[\psi _{0}(x)]=0,$

czyli

${\frac {1}{\sqrt {2}}}\left({\frac {1}{\xi }}x+\xi {\frac {\mathrm {d} }{\mathrm {d} x}}\right)\psi _{0}(x)=0.$

Jest to równanie różniczkowe 1-go stopnia. Po znalezieniu rozwiązania i podstawieniu z powrotem wyrażenie na $\xi ,$ otrzymuje się:

$\psi _{0}(x)=C_{0}\exp \left(-{\frac {m\omega }{\hbar }}{\frac {x^{2}}{2}}\right),$

gdzie $C_{0}$ – stała normalizacyjna. Funkcja ta jest funkcją wykładniczą, symetrycznie zanikającą w nieskończonościach, mającą maksimum dla $x=0.$ Oznacza to, że dla energii drgań zerowych największe jest prawdopodobieństwo znalezienia oscylatora w stanie równowagi (porównaj wykresy gęstości prawdopodobieństw $|\psi _{n}(x)|^{2}$ umieszczone w poprzednim rozdziale).

Funkcje własne stanów wzbudzonych

Za pomocą operatora kreacji ${\hat {a}}^{\dagger }$ można teraz obliczyć funkcje falowe stanów wzbudzonych:

$\psi _{1}={\hat {a}}^{\dagger }\psi _{0},$

$\psi _{2}={\frac {1}{\sqrt {2}}}{\hat {a}}^{\dagger }\psi _{1}={\frac {1}{\sqrt {2}}}({\hat {a}}^{\dagger })^{2}\psi _{0},$

$\psi _{n}(x)={\frac {1}{\sqrt {n!}}}({\hat {a}}^{\dagger })^{n}\psi _{0}.$

Do obliczenia stanów wzbudzonych wystarczy znaleźć wynik działania potęgi operatora kreacji na stan zerowy. Operator ${\hat {a}}^{\dagger }$ w jawnej postaci uzyskuje się analogicznie jak w przypadku operatora ${\hat {a}}{:}$

${\hat {a}}^{\dagger }={\frac {1}{\sqrt {2}}}\left({\frac {1}{\xi }}x-\xi {\frac {\mathrm {d} }{\mathrm {d} x}}\right).$

Powyższy operator różniczkowy, działając n-krotnie na funkcję wykładniczą $\psi _{0}=C_{0}\exp \left(-{\frac {m\omega }{\hbar }}{\frac {x^{2}}{2}}\right),$ reprodukuje ten sam czynnik wykładniczy, pomnożony przez wielomian $n$-tego rzędu względem $x.$

Ostatecznie otrzymamy:

$\psi _{n}(x)={\frac {1}{\sqrt {n!}}}({\hat {a}}^{\dagger })^{n}\psi _{0}(x)=C_{n}H_{n}\left({\frac {x}{\xi }}\right)\psi _{0}(x),$

gdzie:

$H_{n}(x)=(-1)^{n}\mathrm {e} ^{x^{2}}{\frac {\mathrm {d} ^{n}}{\mathrm {d} x^{n}}}\mathrm {e} ^{-x^{2}}$

są wielomianami Hermite’a, $C_{n}$ jest stałą normalizacyjną, $\xi ={\sqrt {\frac {\hbar }{m\omega }}}.$

Algebra Heisenberga

Powyżej zdefiniowane operatory tworzą grupę operatorów

${\hat {X}}_{i}=\left\{{\hat {I}},{\hat {a}},{\hat {a}}^{\dagger },{\hat {N}}\right\}.$

Grupa ta rozpina algebrę Heisenberga (która jest jedną z algebr Liego) o następujących komutatorach:

$[{\hat {a}},{\hat {a}}^{\dagger }]=1,$

$[{\hat {a}},{\hat {a}}]=[{\hat {a}}^{\dagger },{\hat {a}}^{\dagger }]=0,$

$[{\hat {N}},{\hat {a}}]=-{\hat {a}},$

$[{\hat {N}},{\hat {a}}^{\dagger }]={\hat {a}}^{\dagger },$

$[{\hat {I}},{\hat {X}}_{i}]=0.$

Fermionowy oscylator harmoniczny

(1) Przypomnijmy, że bozonowy oscylator harmoniczny opisuje hamiltonian:

${\hat {H}}={\frac {1}{2}}\hbar \omega ({\hat {a}}^{\dagger }{\hat {a}}+{\hat {a}}{\hat {a}}^{\dagger }).$

(2) Analogicznie definiuje się fermionowy oscylator harmoniczny – opisuje go hamiltonian:

${\hat {H}}={\frac {1}{2}}\hbar \omega ({\hat {c}}^{\dagger }{\hat {c}}-{\hat {c}}{\hat {c}}^{\dagger }),$

przy czym spełnione są związki:

$\{{\hat {c}},{\hat {c}}^{\dagger }\}=1,$

$\{{\hat {c}},{\hat {c}}\}=\{{\hat {c}}^{\dagger },{\hat {c}}^{\dagger }\}=0,$

gdzie:

${\hat {c}},{\hat {c}}^{\dagger }$ – operatory anihilacji i kreacji pojedynczego fermionu (o częstotliwości $\omega$),

$\{..,..\}$ – antykomutator operatorów, tj.

$\{{\hat {c}},{\hat {c}}^{\dagger }\}={\hat {c}}{\hat {c}}^{\dagger }+{\hat {c}}^{\dagger }{\hat {c}},$

$\{{\hat {c}},{\hat {c}}\}={\hat {c}}{\hat {c}}+{\hat {c}}{\hat {c}}.$

Operator liczby fermionów ma postać

${\hat {N}}={\hat {c}}^{\dagger }{\hat {c}}.$

Powyżej zdefiniowane operatory tworzą grupę operatorów

${\hat {X}}_{i}=\{{\hat {I}},{\hat {c}},{\hat {c}}^{\dagger },N={\hat {c}}^{\dagger }{\hat {c}}\}.$

Grupa ta rozpina algebrę gradowaną o następujących komutatorach

$[{\hat {N}},{\hat {c}}]=-{\hat {c}},$

$[{\hat {N}},{\hat {c}}^{\dagger }]={\hat {c}}^{\dagger },$

$[{\hat {I}},{\hat {X}}_{i}]=0.$

(3) Hamiltonian fermionowy można przekształcić do postaci:

${\hat {H}}=\hbar \omega ({\hat {c}}^{\dagger }{\hat {c}}-{\frac {1}{2}})=\hbar \omega {\hat {c}}^{\dagger }{\hat {c}}+E_{0},$

gdzie $E_{0}=-{\frac {1}{2}}\hbar \omega$ jest energią stanu podstawowego.

(4) Zakaz Pauliego

Reguła komutacyjna $\{{\hat {c}}^{\dagger },{\hat {c}}^{\dagger }\}=0$ wyraża zakaz Pauliego:

• fermionowy oscylator harmoniczny istnieje tylko w stanie próżni $|0\rangle$ lub w pierwszym stanie wzbudzonym $|1\rangle ={\hat {c}}^{\dagger }|0\rangle ,$
• drugi stan wzbudzony $|2\rangle =({\hat {c}}^{\dagger })^{2}|0\rangle$ nie istnieje, bo z reguł antykomutacyjnych wynika, iż $({\hat {c}}^{\dagger })^{2}=0$ (czyli $|2\rangle =({\hat {c}}^{\dagger })^{2}|0\rangle =0$).

Operator Hamiltona – to operator energii. Dozwolonym stanom oscylatora odpowiadają wartości własne operatora Hamiltona:

• $E_{0}=-{\frac {1}{2}}\hbar \omega$ – dla stanu $|0\rangle ,$
• $E_{1}=+{\frac {1}{2}}\hbar \omega$ – dla stanu $|1\rangle .$

Kwantowa teoria pola

Dokładnego opisu pól fizycznych dostarcza kwantowa teoria pola. Kwantowanie pól fizycznych polega na zastąpieniu wielkości polowych (skalarnych czy wektorowych) operatorami. Przy tym pola dzieli się na bozonowe (o spinie całkowitym) i fermionowe (o spinie połówkowym). Pola opisywane są za pomocą hamiltonianów analogicznych do hamiltonianu oscylatora harmonicznego bozonowego lub fermionowego.

Stan pola bozonowego (np. przypisanego fotonom, mające spin 1) opisuje się jako sumę wzbudzeń wielu oscylatorów bozonowych, z których każdy ma inną częstotliwość drgań i właściwą sobie energię (przy czym energie te są skwantowane, tj. mogą przyjmować dyskretne wartości $E_{n}(\omega )=\hbar \omega \left(n+{\frac {1}{2}}\right),\quad n=0,1,\dots ,$ podobnie jak oscylator bozonowy). Częstotliwości mogą zaś przyjmować wartości dodatnie, continuum.

Stan pola fermionowego (np. przypisanego elektronowi, mające spin 1/2) opisuje się jako sumę wzbudzeń wielu oscylatorów fermionowych, które dla danej częstotliwości mogą przyjmować tylko dwa stany energii $E_{0}(\omega )=-{\frac {1}{2}}\hbar \omega$ oraz $E_{1}(\omega )=+{\frac {1}{2}}\hbar \omega .$ Częstotliwości mogą zaś przyjmować wartości dodatnie, continuum.

Supersymetria

Z połączenia hamiltonianu bozonowego i fermionowego tworzy się hamiltonian

${\hat {H}}={\frac {1}{2}}\hbar \omega (\{{\hat {a}}^{\dagger },{\hat {a}}\}+[{\hat {c}}^{\dagger },{\hat {c}}]).$

Hamiltonian ten ma łącznie dodatkową symetrię, zwaną supersymetrią – miesza ona bozonowe stopnie swobody z fermionowymi. Symetria ta jest generowana przez operatory: ${\hat {Q}}={\sqrt {2\omega }}{\hat {a}}^{\dagger }{\hat {c}}$ oraz ${\hat {Q}}^{\dagger }={\sqrt {2\omega }}{\hat {c}}^{\dagger }{\hat {a}},$ które spełniają relację:

$\{{\hat {Q}},{\hat {Q}}^{\dagger }\}=2{\hat {H}}.$

Ta własność jest podstawą konstrukcji supersymetrycznej teorii pola.

Zobacz też

• kwantowa teoria pola
• kwantowa teoria pola grawitacyjnego
• kwantowa teoria pola w zakrzywionej czasoprzestrzeni
• oscylator anharmoniczny

Bibliografia

• R.L. Liboff: Wstęp do mechaniki kwantowej. Warszawa: PWN, 1987, s. 164–180.
• Claude Cohen-Tannoudji, Bernard Diu, Frank Laloë, Quantum Mechanics, Vol. I, 1991. Wiley, New-York, ISBN 0-471-16433-X, s. 481–541.