Die Kugelfunktionen in der IT

Die kartesischen Koordinaten lassen sich durch Transformationen in Polarkoordinaten umwandeln. Die Darstellung eines Punktes bezüglich der zwei Systeme sind dann durch die folgenden Beziehungen miteinander verknüpft:

$$X \ = \ R\, sin(\theta)\, cos(\phi)$$ (9.4)

$$Y \ = \ R\, sin(\theta)\, sin(\phi)$$ (9.5)

$$Z \ = \ R\, cos(\theta)$$ (9.6)

$$R^2 \ = \ X^2\, + \, Y^2\, + \, Z^2$$ (9.7)

$$\phi \ = \ arctan \, \left(\, \displaystyle\frac {\, Y\, }{X}\, \right) \quad$$ (9.8)

$$\theta \ = \ arctan \, \left(\, \displaystyle\frac {\sqrt{\, X^2\, + Y^2\, }}{Z}\, \right) \quad$$ (9.9)

Jede infinitesimale Änderungen der Variabeln $(\, X, Y, Z\, )$ der kartesischen Koordinaten führt zu einer afinen Transfomation der Polarkoordinaten und damit zur Änderungen der Variabeln $(\, R, \theta, \phi\, )$. Nun werden wir die Matrix dieser Transformation an Hand der Gleichungen (9.4), (9.5), und (9.6) ableiten. Differenzieren wir diese Gleichungen, so stellen wir fest, daß folgende Differentialquotienten bei der Ableitungen der kartesischen Variabeln die Komponenten der gesuchte Transformationmatrix ergeben:

$$\fbox {$\rule[-4mm]{0cm}{1cm}R_{1k} \ = \ \left( \ 0 \ , \ \ (\di... ...artial R}Y)\ , \ (\displaystyle\frac {\partial}{\partial R}Z)\, \right)\quad $}$$ (9.10)

$$\fbox {$\rule[-4mm]{0cm}{1cm}\theta_{1k} \ = \ \left( \ 0 \ , \ \... ...e\frac {1}{R}\displaystyle\frac {\partial}{\partial \theta}Z)\, \right)\quad $}$$ (9.11)

$$\fbox {$\rule[-4mm]{0cm}{1cm}\phi_{1k} \ = \ \left( \ 0 \ , \ \ (... ...R sin(\theta)} \displaystyle\frac {\partial}{\partial \phi}Z)\, \right)\quad $}$$ (9.12)

Nach Ausführung der Differentiation erhalten wir:

$$\fbox {$\rule[-4mm]{0cm}{1cm}R_{1k} \ = \ \left( \ 0 \ , \ \ (sin... ... cos(\phi))\ , \ (sin(\theta)\, sin(\phi))\ , \ (cos(\theta))\, \right)\quad $}$$ (9.13)

$$\fbox {$\rule[-4mm]{0cm}{1cm}\theta_{1k} \ = \ \left( \ 0 \ , \ \... ...cos(\phi))\ , \ (cos(\theta)\, sin(\phi))\ , \ -(sin(\theta))\, \right)\quad $}$$ (9.14)

$$\fbox {$\rule[-4mm]{0cm}{1cm}\phi_{1k} \ = \ \left( \ 0 \ , \ \ -(sin(\phi))\ , \ (cos(\phi))\ , \ (\, 0 \, ))\, \right)\quad $}$$ (9.15)

Die drei Transformationmatrix sind unter sich wie folgt verknüpft:

$$\fbox {$\rule[-4mm]{0cm}{1cm}R_{jk} \ . \ R_{kj} \ = \ \delta_j^k \quad $}$$ (9.16)

$$\fbox {$\rule[-4mm]{0cm}{1cm}\theta_{jk} \ . \ \theta_{kj} \ = \ \delta_j^k \quad $}$$ (9.17)

$$\fbox {$\rule[-4mm]{0cm}{1cm}\phi_{jk} \ . \ \phi_{kj} \ = \ \delta_j^k \quad $}$$ (9.18)

Untereinander weisen sie folgende Eigenschaften auf:

$$\fbox {$\rule[-4mm]{0cm}{1cm}R_{jk} \ . \ \theta_{jk} \ = \ \phi_... ... \fbox {$\rule[-4mm]{0cm}{1cm}\theta_{jk} \ . \ R_{jk} \ = \ \phi_{kj} \quad $}$$ (9.19)

$$\fbox {$\rule[-4mm]{0cm}{1cm}\theta_{jk} \ . \ \phi_{jk} \ = \ R_... ... \fbox {$\rule[-4mm]{0cm}{1cm}\phi_{jk} \ . \ \theta_{jk} \ = \ R_{kj} \quad $}$$ (9.20)

$$\fbox {$\rule[-4mm]{0cm}{1cm}\phi_{jk} \ . \ R_{jk} \ = \ \theta_... ... \fbox {$\rule[-4mm]{0cm}{1cm}R_{jk} \ . \ \phi_{jk} \ = \ \theta_{kj} \quad $}$$ (9.21)

Bezeichen wir den infinitesimalen polaren Viererortstensor mit $dY_{jk}$, so können wir für seine Komponente schreiben:

$$\fbox {$\rule[-4mm]{0cm}{1cm}dY_{1k} \ = \ \left( \, ic\, dT \ , ... ...\, R\, d\theta\, \ , \ \phi_{1k}\, R\, sin(\theta)\ d\phi\, )\, \right)\quad $}$$ (9.22)

Die Überführung der kartesischen in Polarkoordinaten wurde realisiert, in dem wir jeden Komponenten des Linienelements $(dR, \, R\, d\theta, R\, sin(\theta)\, d\phi)$ eine explizite orthogonale Transformationen zugeordnet haben. Die so entstandene Raumelemente sind Komponenten der polaren Vierertensoren mit den Einheitsvektoren der kartesischen Koordinaten. Wir weisen an dieser Stelle besonders darauf hin, daß im Gegensatz zur kartesischen Koordinaten $(X\, , Y\, , Z)$ wir mit fünf Variabln $(dR\, ,d\theta\, ,d\phi\, ,R\, ,sin(\theta))$ der Polarkoordinaten zu tun haben. Für das Abstandquadrat erhalten wir unter Berücksichigung von (9.16), (9.17) und (9.18):

$$\fbox {$\rule[-4mm]{0cm}{1cm}dS^2 \ = \ (\, -c^2\, dT^2\, + \, (\... ...^2 \, + \, R^2\, d\theta^2\, + \, R^2\, sin(\theta)^2\, d\phi^2\, )\, )\quad $}$$ (9.23)

Analog zu (9.22) kann man die Vierergradienten der Polarkoordinaten ermitteln. Man findet sie:

$$\fbox {$\rule[-4mm]{0cm}{1cm}\partial_{1k}\ = \ \left(\, i\displa... ...sin(\theta)} \displaystyle\frac {\partial}{\partial \phi}\, )\, \right)\quad $}$$ (9.24)

Dieses Ergebnis in (9.23) kommt der Vorstellung eines orthogonalen polaren Ortsvektors mit Einheitsvektoren $(\vec{e_r}\, , \vec{e_\theta}\, , \vec{e_\phi})$ nahe. Diese Äquivalenz zu den kartesischen Koordinaten bleibt aufrechtenthalten, solange die Berechnungen (die Tensorprodukte) keine Vierergradienten enthalten. Hierzu betrachten wir die vereinfachte Gleichung (6.18):

$$\fbox {$\rule[-4mm]{0cm}{1cm}(\, -i {\bf div} \vec{D^-}_{(\vec{R})}\, \ , \ \ {\bf rot} \displaystyle\frac {\vec{H}^+}{c}\, ) \quad $}$$ (9.25)

Man findet oft in der Literatur Integralgleichungen, die zur Definition^9.1$^)$ (9.25) der div und rot der Vektoren beitragen: Nach Gauß - Ostrogradskischen Formel ist das Volumenintegral der div gleich dem Vektorfluß der Oberfläche des Bereiches:

$$\fbox {$\rule[-4mm]{0cm}{1cm}{\bf div} \vec{D^-}_{(R,\theta,\phi)... ...playstyle\oint \vec{D^-}_{(R,\theta,\phi)} . d\vec{F}\, }{\triangle V} \quad $}$$ (9.26)

und nach Stockesschen Formel ist die Zirkulation des Vektorfeldes $$\frac {\vec{H}^+}{c}$$ längs der Kontur $(\vec{S})$ einer Fläche $(\vec{F})$ gleich dem Integral der Normalkomponente der rot über die Fläche selbst, d.h. gleich dem Wirbelfluß durch die Fläche:

$$\fbox {$\rule[-4mm]{0cm}{1cm}{\bf rot} \vec{H^+}_{(R,\theta,\phi)... ...playstyle\oint \vec{H^+}_{(R,\theta,\phi)} . d\vec{S}\, }{\triangle F} \quad $}$$ (9.27)

Unsere Beweislage wird sich aber nicht auf diese Gleichungen stützen. Lediglich werden wir sie zum Vergleich heranziehen. Wir errechnen zuerst das Volumen und das Flächenelement aus den Komponenten in (9.22). Man erhält:

$$\fbox {$\rule[-4mm]{0cm}{1cm}dV_1^k \ = \ dR\, R\, d\theta\, R\, ... ..._{jk}\, ) \ = \ - \ R^2\, dR\, d\theta\, sin(\theta)\, d\phi \delta_1^k\quad $}$$ (9.28)

$$\fbox {$\rule[-4mm]{0cm}{1cm}d\vec{F} \ = \ (\, (\, R^2\, sin(\th... ...d\phi\, )\, \theta_{1k}\, \ , (\, R\, dR\, d\theta\, )\, \phi_{1k}\, ) \quad $}$$ (9.29)

Nun errechnen wir das Tensorprodukt $\partial_{1k}. \acute{A_{kj}}$ aus: In ${\bf div} \vec{D^-}_{(R,\theta,\phi)}$ können außer $(dR\, ,d\theta\, ,d\phi)$ noch $(R^2$ und $sin(\theta))$ und in ${\bf rot}\vec{H^+}_{(R,\theta,\phi)}$ die Größen $(R$ und $sin(\theta))$ in die Berechnungen eingehen. Bei der Ableitungen werden die Transformationenmatrix konstant gehalten und nur ihre Faktoren, des Flächen

$$F_{(R,\theta,\phi)}=(\, R^2\, sin(\theta)\, ,R\, sin(\theta)\, ,R\, )$$ (9.30)

bzw. des Linienelements

$$S_{(R,\theta,\phi)}= (1\, , R\, , R\, sin(\theta))$$ (9.31)

differenziert. Führen wir die Differenzation durch, so finden wir für den polaren $\vec{D}_{(R,\theta,\phi)}$ und axialen Vektor $\vec{H}_{(R,\theta,\phi)}$

$$\left( \displaystyle\frac {1}{F_R}\displaystyle\frac {\partial}{\... ... sin(\theta)} \displaystyle\frac {\partial}{\partial \phi}(S_RH^+_R)\theta_{1k}$$ (9.32)

$$\displaystyle\frac {1}{S_\theta}\displaystyle\frac {\partial}{\pa... ...eta)} \displaystyle\frac {\partial}{\partial \phi}(S_\theta H^+_\theta ) R_{1k}$$ (9.33)

$$-\displaystyle\frac {1}{S_\phi}\displaystyle\frac {\partial}{\par... ...(\theta)} \displaystyle\frac {\partial}{\partial \phi}(F_\phi D^-_\phi )\right)$$ (9.34)

Berücksichtigen wir diesen Umstand bei der Ableitungen nach $(R\, ,\theta$ und $\phi)$, so finden wir für div eines polaren Vektors $({\bf rot}\, \vec{D}=0)$ im Fall der Oberflächenabhängigkeit den Matrix $(R_{kj}\, ,\theta_{kj}\, , \phi_{kj})$ übereinstimmend mit den Termen in den runden Klammern von (9.32), (9.33) und (9.34):

$$\fbox {$\rule[-4mm]{0cm}{1cm}{\bf div}\vec{D^-}_{(R,\theta,\phi)}... ...ac {1}{R sin(\theta)} \displaystyle\frac {\partial}{\partial \phi}D^-_{\phi} $}$$ (9.35)

Entsprechende Formel erhalten wir für die rot eines axialen Vektors $({\bf div}\, \vec{H}=0)$ , falls wir die genannten Matrixen als linear (9.32), (9.33), (9.34) betrachten:

$$\fbox {$\rule[-4mm]{0cm}{1cm}{\bf rot} (\vec{H^+}_{(R,\theta,\phi... ...style\frac {\partial}{\partial \theta}H^+_R\, ) \phi_{1k} \end{array}\right) $}$$ (9.36)

Äquivalent zu div kann man den Laplace Operator ableiten. Analog zu (9.33) erhalten wir für das Drehimpulsquadrat:

$$\fbox {$\rule[-4mm]{0cm}{1cm}(\hbar {\bf grad})^2 = \hbar^2 \left... ...c {1}{R sin(\theta)} \displaystyle\frac {\partial}{\partial \phi})^2 \right) $}$$ (9.37)

Diese, sind die Differentialgleichungen der Kugelfuntionen. Sie werden in radialen und winkelabhängigen Kugelfunktionen unterteilt.