5. in 6. predavanje iz Kvantne mehanike
Table of Contents
1. Simetrije - unitarni operatorji
1.1. Translacija
Grafično smo risali realen del \( \psi(x) \), nato pa smo narisali identično funkcijo \( \tilde{\psi}(x) \), ki pa je premaknjena za \( s \) in velja
\[ \tilde{\psi}(x) = \psi(x - s) \]
Za unitaren operator je značilno, da se ohranja skalarni produkt. Translacijo zapišemo
\begin{align*} U(s) \psi(x) &= \tilde{\psi}(x) = \psi(x - s) \\ &= \psi(x) + s \frac{\partial }{\partial x} \psi(x) + \frac{s ^2}{2!} \frac{\partial ^2 }{\partial x ^2} \psi(x) + \ldots \\ &= \exp \left\{ -s \frac{\partial }{\partial x} \right\}\psi(x) \end{align*}V zadnji vrstici smo dobili zametek operatorja, vendar ga bomo sedaj nadomestili z gibalno količino \( p_x = - \mathrm{i} \hbar \frac{\partial }{\partial x} \)
\begin{equation} \label{eq:1} U(s) = \exp \left\{ - s \frac{\partial }{\partial x} \right\} = \exp \left\{ - \mathrm{i} \frac{s p_x}{\hbar} \right\} \end{equation}Operator translacije \ref{eq:1} lahko zapišemo tudi v treh dimenzijah in sicer
\[ U(\vec{s}) = \exp \left\{ - \mathrm{i} \frac{\vec{s} \cdot \vec{p}}{\hbar} \right\} \]
kjer je produkt \( \vec{s} \cdot \vec{p} \) generator premika.
1.2. Rotacija, vrtenje
Imamo valovno funkcijo \( \psi \left( \vec{r} \right) \), ki se nahaja na ravnini z vektorjem normale \( \vec{n} \). Funkcijo ravninsko zasučemo za kot \( \phi \) tako, da ostane na enaki razdalji od izhodišča. Zasukano funkcijo označimo z \( \psi \left( \vec{r}_1 \right) \), kjer je \( \vec{r}_1 = \vec{r} - \mathrm{d} \vec{r} \).
Rotacija je infinitizimalno majhna, kar pomeni, da nam bo v pomoč tudi limita, ko gre \( \mathrm{d} \phi \to 0 \). Novi radij vektor \( \vec{r}_1 \) je vsota vektorjev \( \vec{r} \) in \( \mathrm{d} \vec{r} \), kjer sta tedva vektorja zaradi infinitizimalnosti pravokotnotna (glej sliko). Velja torej \( \mathrm{d} \vec{r} = \mathrm{d} \phi \vec{n} \times \vec{r} \).
Unitarna opracija bo po razvoju po Taylorju
\begin{align*} U \left( \mathrm{d} \phi \vec{n} \right) \psi(r) &= \psi \left( \vec{r} - \mathrm{d} \vec{r} \right) \\ &= \psi \left( \vec{r} \right) - \mathrm{i} \mathrm{d} \phi \frac{\left( \vec{n} \times \vec{r} \right) \vec{p}}{\hbar} \psi \left( \vec{r} \right) + o \left( \mathrm{d} \phi ^2 \right) \\ &= \left( I - \mathrm{i} \mathrm{d} \phi \frac{\left( \vec{n} \times \vec{r} \right) \cdot \vec{p}}{\hbar} \right) \psi \left( \vec{r} \right) + o \left( \mathrm{d} \phi ^2 \right) \\ &= \left( I - \mathrm{i} \mathrm{d} \phi \frac{\vec{n} \cdot \left( \vec{r} \times \vec{p} \right)}{\hbar} \right) \psi \left( \vec{r} \right) + \ldots \end{align*}Definiramo operator vrtilne količine \( \vec{L} = \vec{r} \times \vec{p} \). Izpeljani operator
\[ U \left( \mathrm{d} \phi \vec{n} \right) = \left( I - \mathrm{i} \mathrm{d} \phi \frac{\vec{n} \cdot \vec{L}}{\hbar} \right) \psi + \ldots \]
še ni uporaben, saj obravnavamo zgolj majhen premik
Namesto tega predpostavimo, da naredimo \( N \) majhnih premikov, kjer gre \( N \) proti neskončno. Torej
\[ U \left( \phi \vec{n} \right) = \left. \left( I - \mathrm{i} \frac{\phi}{N} \frac{\vec{n} \cdot \vec{L}}{\hbar} \right)^N \right|_{N \to \infty} = \exp \left\{ - \mathrm{i} \phi \frac{\vec{n} \cdot \vec{L}}{\hbar} \right\} \]
Dobili smo torej operator rotacije
\[ U \left( \vec{\phi} \right) = \exp \left\{ - \mathrm{i} \frac{\vec{\phi} \cdot \vec{L}}{\hbar} \right\} \]
kjer je \( \vec{\phi} = \phi \cdot \vec{n} \) ter je \( \vec{\phi} \cdot \vec{L} \) generator.
1.3. Inverzija prostora
Obravnavamo operacijo in iščemo operator \( \Pi \), da velja
\[ \Pi \psi \left( \vec{r} \right) = \psi \left( - \vec{r} \right) \]
V eni dimenziji za operator inverzije velja
\begin{align*} \psi(x) &\to \psi(- x) \\ \frac{\mathrm{d} }{\mathrm{d} x} &\to - \frac{\mathrm{d} }{\mathrm{d} x}, \quad p \to -p \\ \frac{\mathrm{d} ^2 }{\mathrm{d} x ^2} &\to \frac{\mathrm{d} ^2 }{\mathrm{d} x ^2}, \quad p ^2 \to p ^2 \end{align*}Naj imamo sod potencial \( V(x) = V(- x) \). Hamiltonian je po definiciji \( H = \frac{ p ^2}{2m} + V \). Zaradi zgoraj izpeljanih lastnosti velja
\[ \Pi H = H \Pi \implies \left[ \Pi H \right] = 0 \]
Torej bodo rešitve enačbe
\[ H \psi(x) = E \psi(x), \]
na kateri smo opravili inverzijo prostora
\[ \Pi H \psi = H \Pi \psi(x) = H \psi(- x) = E \psi(-x) \]
enake
\[ \psi_{\pm} (x) = \frac{1}{\sqrt{2}} \left( \psi(x) \pm \psi(- x) \right), \]
saj je rešitev soda. Če energija \( E \) ni degenerirana, to pomeni, da \( \psi \) soda ali liha.
1.4. Obrat časa (time reversal) Sakurai_QM_book
Obravnavajmo najprej klasično mehanski primer. Imejmo delec, ki se giblje po trajektoriji \( \vec{r} (t) \) s hitrostjo \( \vec{v} \). V primeru, če se začne gibati s hitrostjo \( - \vec{v} \), se bo gibal nazaj po isti trajektoriji. To nam nakažejo tudi gibalne enačbe, saj za obrat časa \( t \to - t \) in posledično \( \vec{r} (t) \to \vec{r} (- t) \) velja
\begin{align*} m \oddv{r} &= \vec{F} (\vec{r}) \\ m \oddv{r} (t) &= \vec{F} (\vec{r} (t)) && t \to - t \\ m \oddv{r} (-t) &= \vec{F} \left( \vec{r} (- t) \right) \end{align*}Velja \( \frac{\mathrm{d} ^2 }{\mathrm{d} t ^2} = \frac{\mathrm{d} ^2 }{\mathrm{d} t ^2} \).
Za valovno enačbo
\[ \frac{\partial ^2 u}{\partial x ^2} = \frac{1}{c ^2} \frac{\partial ^2 }{\partial t ^2} u(x, t) \]
katere rešitev je
\[ u = \sin \left( k x - \omega t \right) = \sin \left( k (x - vt) \right) = \sin \left( k (x + vt) \right) \]
pa je obrat časa zaznaven. Pri klasični mehaniki, če zanemarimo entropijo, ne ločimo med obratom časa.
Obravnavajmo Schrödingerjevo enačbo, kjer je Hamiltonian neodvisen od časa
\[ \mathrm{i} \hbar \frac{\partial }{\partial t} \psi(x, t) = H \psi(x, t) \]
\( \psi(x, - t) \) ne reši enačbe, medtem kot \( \psi^{\ast} (x, - t) \) pa je. To lahko hitro preverimo, saj
\[ \mathrm{i} \hbar \frac{\partial }{\partial t} \psi^{\ast} (x, - t) = H \psi^{\ast} (x, - t) \]
Uvedemo operator konjugacije \( K \), ki \( z \in \mathbb{C} \) priredi njegovo konjugirano vrednosti:
\[ K z = z^{\ast} \]
V splošnem primeru, če \( z \) sledi še kakšna druga zadevščina \( \left\langle \cdot \right\rangle \), ki je lahko kompleksno konjugirana, potem
\begin{equation} \label{eq:3} K z \left\langle \cdot \right\rangle = z^{\ast} K \left\langle \cdot \right\rangle \end{equation}Operator obrata časa v kvantni mehaniki je ravno \( K \). Imamo valovno funkcijo \( \psi(x, t) = \exp \left\{ \mathrm{i} \frac{px}{\hbar} - \mathrm{i} \omega t \right\} \), na katerom delujemo z operatorjem obrata časa
\[K \psi = \psi^{\ast} = \exp \left\{ - \mathrm{i} \frac{p x}{\hbar} + \mathrm{i} \omega t \right\} \\ \]
Če delujemo z operatorjem časa na valovno funkcijo \( \psi(x, - t) = \exp \left\{ \mathrm{i} \frac{px}{\hbar} + \mathrm{i} \omega t \right\} \), pa dobimo
\[ K \psi(x, - t) = \exp \left\{ - \mathrm{i} \frac{px}{\hbar} - \mathrm{i} \omega t \right\}, \]
kjer je \( \omega t = \frac{E t}{\hbar} \). Tukaj tudi pridemo do spoznanja, da operator obrata časa ne obrne časa nazaj, kakor da bi film nazaj vrteli, ampak obrne hitrost!
\begin{align*} \vec{p} & t\to - \vec{p} \\ \vec{L} &\to - \vec{L} \end{align*}[vstavi sliko primera]
Poglejmo si lastna stacionarna stanja
\[ H \psi = E \psi, \]
na katero apliciramo operator obrata časa
\[ K H \psi = H K \psi = H \psi^{\ast} = E \psi^{\ast} \]
Rešitve te enačbe so
\[ \psi_{R, I} = \frac{1}{\sqrt{2}} e^{\mathrm{i} \delta} \left( \psi \pm \psi^{\ast} \right) \]
kjer je \( \psi_R \in \mathbb{R} \) in \( \psi_I = \frac{\mathrm{i}}{\sqrt{2}} \left( \psi - \psi^{\ast} \right) \in \mathrm{i} \mathbb{R} \). V primeru, če \( E \) ni degenerirana, potem \( \psi \in \mathbb{R} \).
1.5. Umeritvena transformacija (ang., nem. gauge)
Umeritvena transformacija ima sledeč učinek
\[ \psi(x, t) \to e^{\mathrm{i} \delta} \psi(x, t) = \tilde{\psi}, \]
kjer je \( \delta \in \mathbb{R} \) konstanta. Če imamo umeritveno transformacijo na vse elemente v našem prostoru, ji pravimo, da je globalna umeritvena transformacija
\[ \phi_n \to e^{\mathrm{i} \delta} \phi_n = \tilde{\phi} \]
Za umeritvene transformacije velja, da ohranjajo skalarni produkt
\[ \left\langle \phi \middle| A \middle| \psi \right\rangle = \left\langle \tilde{\phi} \middle| A \middle| \tilde{\psi} \right\rangle = e^{- \mathrm{i} \delta} e^{\mathrm{i} \delta} \left\langle \phi \middle| A \middle| \psi \right\rangle \]
2. Vrtilna količina
Srečali smo jo (prvič) pri rotacijski simetriji
\[ \vec{L} = \vec{r} \times \vec{p}. \]
Za začetek nas zanima, če je operator vrtilne količine hermitski. Pred nadaljevanjem samo spomnimo na nekaj identitet
\begin{align} x^{\dagger} &= x \\ p ^{\dagger} &= p \\ \left( AB \right) ^{\dagger} &= A^{\dagger} B^{\dagger} \end{align}Vrtilno količino zapišemo kot
\[ L^{\dagger} = \left( \sum\limits_{\beta, \gamma}^{} \epsilon_{\alpha\beta\gamma} x_{\beta} p_{\gamma} \right) ^{\dagger} = \sum\limits_{\beta \gamma}^{} \epsilon_{\alpha \beta\gamma} p_{\gamma} ^{\dagger} x_{\beta} ^{\dagger} = \sum\limits_{\beta \gamma}^{} \epsilon_{\alpha \beta \gamma} x_{\beta} p_{\gamma} = L \]
kjer smo upoštevali nekaj glede komutacije \( x \) in \( p \).
Velikost operatorja vrtilne količine je enaka
\[ L ^2 = \vec{L} \cdot \vec{L} = \sum\limits_{\alpha}^{} L_{\alpha} ^2 \]
2.1. (Glavne) lastnosti vrtilne količine
2.1.1. Uvod
Spomnimo
\begin{align*} \left[ x_{\alpha}, x_{\beta} \right] &= 0 \\ \left[ x_{\alpha}, p_{\beta} \right] &= \mathrm{i} \hbar \delta_{\alpha\beta} \end{align*}Komutacija vrtilnih količin v smeri \( x \) in \( y \) je enaka
\begin{align*} \left[ L_x, L_y \right] &= \left[ y p_z - z p_y, z p_x - x p_z \right] \\ &= \left[ y p_z, z p_x \right] - \underbrace{\left[ p_y z, zp_x \right]}_{= 0} - \underbrace{\left[ yp_z, xp_z \right]}_{= 0} + \left[ p_y z, x p_z \right] \\ &= y p_x \left[ p_z, z \right] + x p_y \left[ z, p_z \right] \\ &= \mathrm{i} \hbar \left( x p_y - y p_x \right) = \mathrm{i} \hbar L_z \end{align*}kjer je \( L_z = xp_y - yp_x \).
V splošnem torej velja
\begin{equation} \label{eq:2} \left[ L_{\alpha}, L_{\beta} \right] = \mathrm{i} \hbar \epsilon_{\alpha \beta \gamma} L_{\gamma} \end{equation}Enačba \ref{eq:2} je zelo pomembna zveza.
Velja tudi
\[ \left[ L_{\alpha}, A_{\beta} \right] = \mathrm{i} \hbar \epsilon_{\alpha \beta \gamma} A_{\gamma}, \]
če je \( \vec{A} \) vektorski operator. Vektorski operator je operator, ki se pri rotaciji prostora, transformira kot vektor (aka uporabiš rotacijsko matriko). Primere takih vektorski operatorjev je so vrtilna količina \( \vec{L} \), radij vektor \( \vec{r} \) in gibalne količine \( \vec{p} \).
2.1.2. z \( \vec{L} \) komutirajoči operatorji
Iščemo operatorje, ki komutirajo z vrtilno količino, saj v tem primeru lahko uporabimo (neke?) lastne vrednosti. Matematično zapisano
\[ \left[ \vec{L}, A \right] = 0 \]
Rotacijski operator na \( A \) zapišemo kot
\[ U \left( \vec{\phi} \right) A = A U \left( \vec{\phi} \right), \]
kjer smo to zapisali za isto lastnost kot \ref{eq:3}. Če za vsako rotacijo \( \phi \) velja \( \left[ U, A \right] = 0 \), potem je \( A \) skalar. Posledično velja \( \left[ \vec{L}, A \right] = 0 \).
Količine \( A = 0 \), \( \vec{r}, \vec{r} = r ^2 \), \( \vec{p} \cdot \vec{p} = p ^2 \), \( V \left( \vec{r} \right) = V(r) \), \( H = \frac{p ^2}{2m} + V(x) \) in \( L ^2 \) komutirajo z vsako komponento.
2.1.3. Lestvični operator
Uvedemo lestvični operator (ang. ladder operator), ki je definiran kot
\[ L_{\pm} = L_x \pm \mathrm{i} L_y \]
Hitro lahko preverimo, da velja \( \left( L_{+} \right)^{\dagger} = L_- \) in \( \left[ L_{\pm}, L ^2 \right] = 0 \). Lestvična operatorja nam bosta zviševala/zniževala stanja podobna kot sta to počela \( a^{\dagger} \) in \( a \). Takrat smo uvedli tudi številski operator \( n \). Pri vrtilni količini bo vlogo številskega operatorja prevzel \( L_z \).
\begin{align*} \left[ L_z, L_{\pm} \right] &= \left[ L_z, L_x \pm \mathrm{i} L_y \right] \\ &= \left[ L_z, L_x \right] \pm \mathrm{i} \left[ L_z, L_y \right] \\ &= \mathrm{i} \hbar L_y \pm \mathrm{i} \left( - \mathrm{i} \hbar \right) L_x \\ &= \hbar \left( \pm L_x + \mathrm{i} L_y \right) = - \hbar L_{\pm} \end{align*}Izbrali smo si komponento \( L_z \) zaradi dogovora, saj je problem simetričen na vse 3 osi. Če hitro povzamemo, velja
\[ \left[ L_z, L_{\pm} \right] = \pm \hbar L_{\pm} \]
Pogledalo smo si \( a^{\dagger} a \), zato bomo tukaj tudi pogledali rezultat analogen temu:
\begin{equation} \label{eq:5} \begin{aligned} L_{\pm} L_{mp} &= \left( L_x \pm \mathrm{i} L_y \right) \left( L_x \mp \mathrm{i} L_y \right) \\ &= L_x ^2 + L_y ^2 \pm \mathrm{i} L_y L_x \mp \mathrm{i} L_x L_y \\ &= L ^2 - L_z ^2 \pm \hbar L_z \end{aligned} \end{equation}Komutirajoče zveze
\[ \left[ L_+, L_- \right] = 2 \hbar L_z \]
ne bomo potrebovali according to Ramšak. Spomnimo, da se rezultat razlikuje od komutiranja anihilacijskega in kreacijskega operatorja, ki pa je bil
\[ \left[ a^{\dagger}, a \right] = -1 \]
2.1.4. Lastne vrednosti vrtilne količine
Operator vrtilne količine je definiran kot
\[ \vec{L} = \vec{r} \times \left( - \mathrm{i} \hbar \nabla \right), \]
preko katerega bi lahko iskali s pomočjo Schrödingerjeve enačbe lastne vrednosti. To smo (baje) že delali in bomo namesto tega iskale lastne vrednosti s pomočjo operatorjev.
Uporabili bomo identiteti
\[ \left[ L ^2, L_z \right] = 0, \quad \left[ L_z, L_+ \right] = \hbar L_+ = L_z L_+ - L_+ L_z . \]
Naj bo \( m \in \mathbb{R} \) lastna vrednost vrtilne količine
\begin{equation} \label{eq:4} L_z \left| m \right\rangle = m \hbar \left| m \right\rangle \end{equation}Na lastno stanje vplivamo z operatorjem \( L_+ \) in iz zgornje identitete sledi
\[ L_z L_+ \left| m \right\rangle = \left( L_+ L_z + \hbar L_+ \right) \left| m \right\rangle = \hbar (m + 1)L_+ \left| m \right\rangle, \]
kjer smo upoštevali \ref{eq:4}.
Zanimajo nas tudi lastnosti lastne vrednosti \( \lambda \), za katere velja
\[ L ^2 \left| m \right\rangle = \lambda \left| m \right\rangle. \]
Operator \( L ^2 = \vec{L} \cdot \vec{L}\) je produkt dveh vektorjev, zato postopamo malo drugače pri lastnih stanjih
\begin{align*} \left\langle m \middle| L ^2 \middle| m \right\rangle &= \left\langle m \middle| \sum\limits_{\alpha}^{} L_{\alpha} ^2 \middle| m \right\rangle \\ &= \sum\limits_{\alpha}^{} \left\langle L_{\alpha}^{\dagger} m \middle| L_{\alpha} m \right\rangle \ge 0 \end{align*}Iz tega sledi, da ker je \( \left\langle m \middle|L ^2 \middle| m \right\rangle \ge 0 \), preko enačbe
\[ \left\langle m \middle| L ^2 \middle| m \right\rangle = \lambda \left\langle m \middle| m \right\rangle, \]
tudi \( \lambda \ge 0 \), saj je vedno \( \left\langle m \middle| m \right\rangle \ge 0 \)
Kakšna pa je \( \lambda \)? Zaradi komutacije velja zveza
\[ L ^2 L_+ \left| m \right\rangle = L_+ L ^2 \left| m \right\rangle = \lambda L_+ \left| m \right\rangle, \]
iz česar skupaj z zgornjim spoznanjem sledi
\[ \lambda = l (l + 1) \hbar ^2 \ge 0 \]
Če si narišemo graf funkcije \( \lambda(l) \), dobimo parabolo značilno za kvadratno funkcijo. Zaradi pogoja, da je večja ali enaka \( 0 \), sledi, da so primerne lastne vrednosti \( l \) tudi večje ali enake \( 0 \), saj z njimi zadostimo vse možne vrednosti \( \lambda \).
Velja(?)
\[ \left\langle L_+ \psi_{lm} \middle| L_+ \middle| \psi_{lm} \right\rangle = \left\langle lm \middle| L_- L_+ \middle| lm \right\rangle \ge 0 \]
Preko zveze \ref{eq:5} se nam izraz preobrazi v
\begin{align*} \left\langle lm \middle| L_- L_+\middle| lm \right\rangle & \ge 0 \\ \left\langle lm \middle| \left( L ^2 - L_z ^2 - \hbar L_z \right) \middle| L_z \right\rangle & \ge 0 \\ \left\langle lm \middle| \left| \left( l (l + 1) \hbar ^2 - m ^2 \hbar ^2 - \hbar ^2 \eta \right) \right| \middle| lm \right\rangle & \ge 0 \\ \left( l(l + 1) - m(m + 1) \right) \hbar ^2 \left\langle lm \middle| lm \right\rangle &\ge 0 \end{align*}Ker je \( \left\langle lm \middle| lm \right\rangle \) vedno večji od \( 0 \), to pomeni, da je tudi funkcija predfaktor večji od \( 0 \).
Iz neenakosti
\[ l(l + 1) - m(m + 1) \ge 0 \implies f(l) = l(l + 1) \ge f(m) = m(m + 1) \]
velja \( l \ge m \) za \( m \ge 0 \). Zvezo lahko zapišemo tudi kot \( 0 \le m \le l \).
Kaj pa v primeru, da je \( m \le 0 \)? Potem zapišemo \( m = - \left| m \right| \) in sledi
\begin{align*} l(l + 1) & \ge - \left| m \right| \left( - \left| m \right| - 1\right) = \left| m \right| \left( \left| m \right| + 1 \right) \\ \left| m \right| &\ge l \end{align*}oziroma z drugimi besedami, \( m \) je spodaj omejen z \( -l \) in ko združimo skupaj, dobimo
\[ - l \le m \le l \]
Do tega lahko pridemo analogno kakor za pozitivne \( m \), samo da uporabimo zvezo \( \left\langle L_- lm \middle| L_- \middle| lm \right\rangle \)
Lastno stanje zapišemo
\[ L_+ \left| lm \right\rangle = c \left| l, m + 1 \right\rangle \]
Potem sledi
\[ \left\langle L_+ \psi_{lm} \middle| L_+ \psi_{lm} \right\rangle = \left| c \right| ^2 \left\langle \psi_{l, m + 1}, \middle| \psi_{l, m + 1} \right\rangle \]
in nadalje
\begin{equation} \label{eq:6} L_{\pm} \left| lm \right\rangle = \hbar \sqrt{l (l + 1) - m (m \pm 1)} \left| l, m \pm 1 \right\rangle \end{equation}Odgovoriti moramo še vprašanje, kaj so lastne vrednosti \( l \)? Gledamo najvišje možno stanje \( \left| ll \right\rangle \), ki ga nižamo z lestvičnim operatorjem
\begin{align*} L_- \left| ll \right\rangle &= c_{l, l - 1} \left| l, l - 1 \right\rangle \\ L_-^2 \left| ll \right\rangle &= c_{l-1, l - 2} \left| l, l - 2 \right\rangle \\ & \vdots \\ L_- ^k \left| ll \right\rangle &= c_{l - k + 1, l - k} \left| l, l - k \right\rangle \end{align*}Najmanjše stanje, do katerega lahko pridemo brez divergence zaradi omejitve \( -l \le m \le m \), je
\[ l - k = - l \implies \ 2l = k \in \mathbb{N}_{0} \]
Iz tega sledi, da \( l = \frac{k}{2} = 0, \frac{1}{2}, 1, \frac{3}{2}, 2, \ldots \).
\[ l = \begin{cases} 0, 1, 2, \ldots &; k \text{ sod} \\ \frac{1}{2}, \frac{3}{2}, \frac{5}{2}, \ldots &; k \text{ lih} \end{cases} \quad\text{ in } - l \le m \le l \]
Rešitev Schrödingerjeve enačbe, na katero smo delovali z operatorjem
\[ L_z = - \mathrm{i} \hbar \frac{\partial }{\partial \phi} = - \mathrm{i} \hbar \left( x \frac{\partial }{\partial y} - y \frac{\partial }{\partial y} \right). \]
Enačba je tako
\[ -\mathrm{i} \hbar \frac{\partial }{\partial \phi} \psi_m (\phi) = m \hbar \psi_m (\phi) \]
in njena rešitev je
\[ \psi_m (\phi) = C e^{\mathrm{i} m \phi}, \]
ki je zvezna funkcija. Ker je \( L_z \) operator v odvisnosti od \( \phi \), pomeni, da se sučemo v ravnini okoli osi \( z \). Zveznost naše rešitve \( \psi \), da je periodična
\[ \psi_m (\phi + 2 \pi) = \psi_m (\phi) \]
Uporabimo nastavek rešitve in zveznost, da dobimo identiteto
\[ e^{i m (\phi + 2 \pi)} = e^{\mathrm{i} m \phi} \implies \ e^{\mathrm{i} 2 \pi m} = 1, \]
iz katere sledi, da je \( m \in \mathbb{Z} \).
Za lihe \( k \) pa \( l \) nima celih vrednosti, kar pomeni, da smo screwed - to predstavlja spin, ki pa ne reši Schrödingerjeve enačbe.
Kot zanimivost za \( l = \frac{1}{2} \) lahko najdemo lastne funkcije, ki rešijo
\[ L ^2 f (\theta, \phi) = \frac{1}{2} \left( \frac{1}{2} + 1 \right) \hbar ^2 f(\theta, \phi), \]
kjer je
\[ f(\theta, \phi) = C \begin{cases} \sin \frac{\theta}{2} e^{\mathrm{i} \frac{\phi}{2}} &; m = \frac{1}{2}\\ \cos \frac{\theta}{2} e^{- \mathrm{i} \frac{\phi}{2}} & ; m = - \frac{1}{2} \end{cases} \]
Rešitev torej obstaja, vendar pa ni zvezna.
2.1.5. Sferične funkcije \( Y_l^m (\theta, \phi) \)
Sferične funkcije ni tako enostavno zgenerirati kakor harmonične funkcije. Nekatere sferične funkcije so
\begin{align*} Y_0^0 &= \frac{1}{\sqrt{4\pi}} \\ Y_1^0 &= \sqrt{\frac{3}{8 \pi}} \cos \theta \\ Y_1^\mathbf{1} &= \sqrt{\frac{3}{8 \pi}} \sin \theta e^{\mathrm{i} \phi} \in \mathbb{C} \end{align*}Obstaja semi-konvencija, da če je indeks \( m \) zapis kot supskript predstavlja sferične funkcije v kompleksni ravnini. V primeru, da sta oba indeksa spodaj zapisana, potem je sferična funkcija realna
\[ Y_{ml} = \sqrt{\frac{3}{8 \pi}} \sin \theta \cos \theta \in \mathbb{R} \]
3. Zapis operatorja z matriko
To poglavje ne spada pod vrtilno količino, samo uporabili smo jo za primer. Na isti način bi lahko izpeljali tudi za harmonični oscilator.
Operator lahko predstavimo z matriko. Gledamo tako valovno funkcijo, da opazujemo samo kotni del
\[ \psi = \sum\limits_{l = 0}^{\infty} \sum\limits_{m = -l}^l c_{lm} \left| lm \right\rangle, \quad \left\langle \vec{r} \middle| lm \right\rangle = Y_l^m (\theta, \phi) \]
Ne želimo se ukvarjati z reševanjem diferencialne enačbe, ampak z operatorji. Lestvični operator izrazimo v neki bazi
\[ L_+ \sum\limits_{l', l, m', m}^{} \left| l' m' \right\rangle \left( L_+ \right)_{l' l m ' m} \left\langle lm \right| \]
Zanima nas, kakšna bo matrika tega operatorja. Opazujmo najprej za \( l = 1 \):
\[ \left( L_+ \right)_{11m' m} = \left\langle 1 m' \middle| L_+ \middle| 1m \right\rangle \]
V nadaljevanju uporabimo enakost \ref{eq:8} in se preobrazi v
\[ \left( L_+ \right) = \hbar \sqrt{1(1 + 1) - m (m + 1)} \left\langle 1 m' \middle| 1, m + 1 \right\rangle = \hbar \sqrt{1(1 + 1) - m (m + 1)} \delta_{m', m + 1} \]
Navada je, da se v vektorskem zapisu piše stanja z višjim \( m \) na vrhu. To si lahko zapomnimo intuitivno tako, da je stanje z največjim \( m \) usmerjeno navpično navzgor, medtem ko stanje z najmanjšim \( m \) pa navpično navzdol. Koeficienti so
\[ c_{lm} = \begin{pmatrix} c_{11} \\ c_{10} \\ c_{1,-1} \end{pmatrix} \]
Koeficiente \( L_+ \psi \), ki jih označimo z \( d_{lm} \) zapišemo v matrični obliki kot
\[ \begin{pmatrix} d_{11} \\ d_{10} \\ d_{1, -1} \end{pmatrix} = \begin{pmatrix} 0 & \sqrt{2} & 0 \\ 0 & 0 & \sqrt{2} \\ 0 & 0 & 0 \end{pmatrix} \begin{pmatrix} c_{11} \\ c_{10} \\ c_{1, -1} \end{pmatrix} \hbar \]
Elementi v stolpcu predstavljajo vrednosti \( m' \), elementi v vrstici pa vrednosti \( m \). Zaradi \( \delta_{m', m + 1} \) je diagonala ničelna.
Torej so matrike lestvičnih operatorjev
\[ L_+ = \hbar \sqrt{2} \begin{pmatrix} 0 & 1 & 0 \\ 0 & 0 & 1 \\ 0 & 0 & 0 \end{pmatrix} = L_x + \mathrm{i} L_y; \quad L_- = \hbar \sqrt{2} \begin{pmatrix} 0 & 0 & 0 \\ 1 & 0 & 0 \\ 0 & 1 & 0 \end{pmatrix} \]
Preko definicije \( L_+ \) in \( L_- \) lahko zapišemo tudi
\begin{align*} L_x &= (L_+ + L_- )\frac{1}{2} = \frac{\hbar}{\sqrt{2}} \begin{pmatrix} 0 & 1 & 0 \\ 1 & 0 & 1 \\ 0 & 1 & 0 \end{pmatrix} \\ L_y &= \frac{\hbar}{\sqrt{2}} \begin{pmatrix} 0 & -\mathrm{i} & 0 \\ \mathrm{i} & 0 & -\mathrm{i} \\ 0 & \mathrm{i} & 0 \end{pmatrix} \\ L_z &= \hbar \begin{pmatrix} 1 & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & -1 \end{pmatrix} \end{align*}4. Centralni potencial
Klasično smo imeli potencial \( V(r) \) in potem je bila sila
\[ \vec{F} = - \frac{\partial V}{\partial r} \frac{\vec{r}}{r}, \]
iz česar sledi, da je vrtilna količina \( \vec{L} = \mathrm{konst} \). V kvantni mehaniki pa imamo Hamiltonian \( H = \frac{p ^2}{2m} + V(r) \), ki komutira
\[ \left[ H, L_z \right] = 0, \ \left[ H, L ^2 \right] = 0 \]
Hamiltonian lahko prevedemo iz kartezičnih tudi v sferične koordinate
\[ H = - \frac{\hbar ^2 \nabla ^2}{2m} + V(r) = - \frac{\hbar}{2m} \left( \frac{1}{r ^2} \frac{\partial }{\partial r} r ^2 \frac{\partial }{\partial r} \right) + \frac{L ^2}{2m r ^2} + V(r) \]
Problem je v tem, da ko mi začne padati [flomaster], si ga ne upam loviti, ker imam potem roke čisto črne.
Drugi del si je precej enostavno zapomniti, vendar dokazovali pa ne bomo. Uporabimo nastavek
\[ \Psi (r, \theta, \phi) = \psi(r) Y_l^m (\theta, \phi) \]
v Schrödingerjevi enačbi s Hamiltonianom v sferičnih koordinatah
\[ \left[ - \frac{\hbar ^2}{2m} \left( \frac{1}{r ^2} \frac{\partial }{\partial r} r ^2 \frac{\partial }{\partial r} \right) + \underbrace{\left( \frac{l(l + 1) \hbar ^2}{2m r ^2} + V(r) \right)}_{V_{ef}} \right] \psi Y_l ^m = E \psi Y_l ^m \]
Uporabimo še nastavek za radialno funkcijo \( \psi(r) = \frac{u(r)}{r} \).
\[ \left[ - \frac{\hbar ^2}{2m} \left( \frac{1}{r ^2} \frac{\partial }{\partial r} r ^2 \frac{\partial }{\partial r} \right)+ \left( \frac{l(l+ 1) \hbar ^2}{2m r ^2} + V(r) \right) \right] \frac{u(r)}{r} = E \frac{u}{r} \]
Prvi seštevanec se poenostavi
\[ \frac{1}{r ^2}\frac{\partial }{\partial r} r ^2 \frac{\partial }{\partial r} \left( \frac{u}{r} \right) = \frac{1}{r ^2} \frac{\partial }{\partial r} r ^2 \left( \frac{u'}{r} - \frac{u}{r ^2} \right) = \frac{1}{r ^2} \frac{\partial }{\partial r} \left( r u' - u \right) = \frac{u''}{r} = \frac{1}{r} \frac{\partial ^2 u}{\partial r ^2} \]
To pomeni, da rešujemo enačbo oblike
\begin{equation} \label{eq:7} -\frac{\hbar ^2}{2m} \frac{\mathrm{d} ^2 }{\mathrm{d} r ^2} u(r) + \left( \frac{l(l + 1) \hbar ^2}{2m r ^2} + V(r) \right) u(r) = E u(r) \end{equation}4.1. Lastnosti rešitev \( u(r) \)
Za razliko od matematikov si ne bomo pogledali vseh vrednosti za \( r \), ampak samo tiste, ki so nam smiselne, ker za več nimamo časa.
Prvi primer, ki ga pogledamo, je \( r\to 0 \), kjer predpostavimo \( \lim_{r \to 0} v ^2 V(r) = 0 \). To za nas pomeni, da v enačbi \ref{eq:7} zanemarimo \( V(r) \), saj pada počasneje kot \( \frac{1}{r ^2} \). Koeficienti se pokrajšajo s koeficienti skritimi v \( E \).
Ostane nam zgolj
\[ \frac{\mathrm{d} ^2 }{\mathrm{d} r ^2} u = \frac{l (l + 1)}{r ^2} u. \]
Uporabimo nastavek \( u = C r^{\lambda} \) in dobimo enačbo
\[ \lambda (\lambda - 1) r ^{\lambda - 2} = \frac{l(l + 1)}{r ^2} r^{\lambda}, \]
katere rešitve so
\[ \lambda_{1, 2} = l + 1, - 1. \]
Rešitve uporabimo v nastavku in je splošna rešitev diferencialne enačbe
\[ u(r) = C_l r^{l + 1} + D \frac{1}{r^l}, l = 0,1, 2, \ldots \]
Vrnemo se nazaj k nastavku za sferične rešitve
\begin{align*} \left\langle \Psi \middle| \Psi \right\rangle &= 1\\ &= \int\limits_0^{\infty} \left| \psi \right| ^2 r ^2 \, \mathrm{d} r \underbrace{\int\limits_0^{2\pi} \, \mathrm{d} \phi \int\limits_0^{\pi} \sin \theta \, \mathrm{d} \theta \left| Y_l^m \right| ^2}_1 \\ &= \int\limits_0^{\infty} \left| u \right|^2 \, \mathrm{d} r \end{align*}Za \( l > 0 \) sledi, da
\[ D_l ^2 \int\limits_0^{\infty} \frac{1}{r^{2l}} \, \mathrm{d} r \to \infty \implies \ D_l = 0 \]
Če pa je \( l = 0 \), pa je naš nastavek
\[ u(r) = C_0 r + C_0. \]
Radialna funkcija je
\[ \psi(r) = \frac{u}{r} = C_0 + \frac{D_0}{r} \]
Uporabimo Gaussov zakon
\[ \nabla \cdot \vec{E} = \frac{\rho \left( \vec{r} \right)}{\epsilon_0} , \]
definicijo potenciala \( \vec{E} = - \nabla \phi \) in imamo točkast naboj \( \rho \left( \dot{r} \right) = e \delta ^3 \left( \vec{r} \right) \). Poissonova enačba je
\[ - \nabla ^2 \frac{e}{4 \pi \epsilon_0 r} = \frac{e \delta ^3 \left( \vec{r} \right)}{\epsilon_0} \]
Isto logiko uoprabimo na naši radialni funkciji
\[ \nabla ^2 \psi = - 4 \pi D_0 \delta ^3 \left( \vec{r} \right) = 0 \implies \ D_0 0 \]
Zaključimo, da je radialna funkcija za majhne \( r \) enaka
\[ \psi(r) = C_l r^l; \ l = 0,1,2 , \ldots \]
Poglejmo si še, kako je z velikimi razdaljami, ko \( r \to \infty \). Gledamo primer, ko je \( E > 0 \) in potencial \( V \) pada proti 0 “dovolj hitro”.
Rešujemo torej enačbo
\[ - \frac{\hbar ^2}{2m} \frac{\mathrm{d} ^2 }{\mathrm{d} r ^2} u = E u, \]
ki ima rešitev
\[ u(r) = C_+ e^{\mathrm{i} k r} + C_- e^{-\mathrm{i} k r} = r \begin{cases} j_l \sim \sin \left( kr - \frac{\pi}{2} l \right)\\ n_l \sim \cos \left( kr - \frac{\pi}{2} l \right) \end{cases} \]
Energija je \( E_{\gamma} = \frac{p ^2}{2m} \) za \( p = \hbar k \)
Tretje možnost je \( E_n < 0 \). Potem bo rešitev
\[ u = D_- e^{- \kappa_n r} + D_+ e^{\kappa _n r} \]
Ko gre \( r \to \infty \), mora biti integral
\[ \int\limits_0^{\infty} \left| u \right| ^2 \, \mathrm{d} r < \infty \]
omejen. Iz tega sledi, da je \( D_+ = 0 \).