Dirac-Gleichung

Die Dirac-Gleichung ist ein System von vier gekoppelten partiellen Differentialgleichungen für die vier Komponentenfunktionen des Dirac-Spinors ${\displaystyle \psi (x)}$. Die Variable ${\displaystyle x}$ steht hier für ${\displaystyle x=(x^{0},x^{1},x^{2},x^{3})\,,}$ worin der obere Index 0 die Zeit ${\displaystyle t=x^{0}}$ und die Indizes 1 bis 3 die Ortskoordinaten ${\displaystyle (x^{1},x^{2},x^{3})}$ bezeichnen.

In natürlichen Maßeinheiten mit ${\displaystyle c=1=\hbar }$ lautet die Dirac-Gleichung für ein ungeladenes Teilchen der Masse ${\displaystyle m}$

${\displaystyle \left[\mathrm {i} \sum _{\mu =0}^{3}\gamma ^{\mu }{\frac {\partial }{\partial x^{\mu }}}-m\right]\,\psi (x)=0\,.}$

Der Ausdruck in eckigen Klammern ist die Standardform eines Dirac-Operators.

Die konstanten Gamma- oder Dirac-Matrizen ${\displaystyle \gamma ^{0},\gamma ^{1},\gamma ^{2}}$ und ${\displaystyle \gamma ^{3}}$ wirken im Raum der vier Komponenten des Spinors und koppeln sie aneinander. Die Produkte von zwei Gamma-Matrizen haben die folgenden Eigenschaften:

${\displaystyle \gamma ^{\mu }\,\gamma ^{\nu }=-\gamma ^{\nu }\,\gamma ^{\mu }\quad (\,\mu ,\nu \in \{0,1,2,3\}\,,\mu \neq \nu ),}$

${\displaystyle \gamma ^{0}\gamma ^{0}=-\gamma ^{1}\gamma ^{1}=-\gamma ^{2}\gamma ^{2}=-\gamma ^{3}\gamma ^{3}=1\;.}$

Damit bilden sie eine Clifford- oder Dirac-Algebra. Wird der Dirac-Operator

${\displaystyle \mathrm {i} \sum _{\mu =0}^{3}\gamma ^{\mu }{\frac {\partial }{\partial x^{\mu }}}-m}$

auf beide Seiten der Dirac-Gleichung angewandt, entkoppeln die vier Differentialgleichungen und man erhält für jede Komponente von ${\displaystyle \psi }$ die Klein-Gordon-Gleichung:

${\displaystyle {\Bigl [}{\frac {\partial ^{2}}{(\partial {x^{0}})^{2}}}-{\frac {\partial ^{2}}{(\partial {x^{1}})^{2}}}-{\frac {\partial ^{2}}{(\partial {x^{2}})^{2}}}-{\frac {\partial ^{2}}{(\partial {x^{3}})^{2}}}+m^{2}{\Bigr ]}\,\psi (x)=0\;.}$

Die zweimalige Anwendung eines Dirac-Operators führt also auf die Klein-Gordon-Gleichung, weshalb die Dirac-Gleichung auch als die „Wurzel“ aus der Klein-Gordon-Gleichung angesehen wird. Für ein Teilchen in einem Impulseigenzustand ergibt die Klein-Gordon-Gleichung (in der Reihenfolge ihrer Terme) ${\displaystyle -E^{2}+{\vec {p}}^{\,2}+m^{2}=0}$, also die relativistische Energie-Impuls-Beziehung eines Teilchens der Masse ${\displaystyle m\,.}$

Jede irreduzible Darstellung der Dirac-Algebra besteht aus ${\displaystyle 4\times 4}$-Matrizen. In der Standard- oder Dirac-Darstellung haben sie die folgende Form (verschwindende Matrixelemente mit Wert Null sind dabei nicht angeschrieben):

${\displaystyle {\begin{array}{c c}\gamma ^{0}={\begin{pmatrix}1&&&\\&1&&\\&&-1&\\&&&-1\end{pmatrix}}\,,&\gamma ^{1}={\begin{pmatrix}&&&1\\&&1&\\&-1&&\\-1&&&\end{pmatrix}}\,,\\\,&\,\\\gamma ^{2}={\begin{pmatrix}&&&-\mathrm {i} \\&&\mathrm {i} &\\&\mathrm {i} &&\\-\mathrm {i} &&&\end{pmatrix}}\,,&\gamma ^{3}={\begin{pmatrix}&&1&\\&&&-1\\-1&&&\\&1&&\end{pmatrix}}\,.\end{array}}}$

Die beiden ersten Komponenten von ${\displaystyle \gamma _{0}}$ bilden also die zweikomponentige Einheitsmatrix, die beiden letzten Komponenten deren Negatives. Analog ergeben die beiden oberen Komponenten der zweiten, dritten bzw. vierten ${\displaystyle \gamma }$-Matrix die drei 2×2-Pauli-Matrizen ${\displaystyle \sigma _{j}}$ und die beiden letzten Komponenten von ${\displaystyle \gamma _{j}}$ deren Negatives. Letztere gehen im nichtrelativistischen Grenzfall wie ${\displaystyle {\tfrac {v}{c}}}$ gegen Null. Damit eignet sich diese Darstellung, die Standarddarstellung, besonders für die Behandlung langsam bewegter Elektronen. In der dazu mathematisch und physikalisch äquivalenten Weyl-Darstellung ist das Spinor-Transformationsverhalten bei Lorentztransformationen besonders einfach, in der ebenfalls äquivalenten Majorana-Darstellung ist die Dirac-Gleichung ein reelles Gleichungssystem. Weitere Darstellungen erhält man durch Äquivalenztransformationen.

Die vier Gamma-Matrizen lassen sich in symbolischer Schreibweise zu dem kontravarianten 4-Vektor

${\displaystyle \gamma ^{\mu }={\begin{pmatrix}\gamma ^{0}\\\gamma ^{1}\\\gamma ^{2}\\\gamma ^{3}\end{pmatrix}}}$

zusammenfassen. Dann hat der erste Term der Dirac-Gleichung die Form eines Skalarprodukts der Vektoren ${\displaystyle \gamma ^{\mu }}$ und ${\displaystyle {\frac {\partial }{\partial x^{\mu }}}}$. Dieses ist bei Lorentztransformation jedoch nicht invariant, denn ${\displaystyle \gamma ^{\mu }}$ bleibt konstant. Die Lorentzinvarianz der Dirac-Theorie ergibt sich erst dadurch, dass der Dirac-Operator auf einen Spinor ${\displaystyle \psi }$ wirkt, dessen vier Komponenten geeignet mittransformiert werden. Im Endergebnis geht damit eine Lösung ${\displaystyle \psi (x)}$ der Dirac-Gleichung durch Lorentztransformation in eine Lösung der entsprechend transformierten Dirac-Gleichung über.

Neben der eben beschriebenen Form im Ortsraum kann die Dirac-Gleichung auch im Impulsraum aufgeschrieben werden. Sie lautet dann

${\displaystyle {\Bigl [}\gamma ^{\mu }p_{\mu }-m{\Bigr ]}\,\psi (p)=0\,,}$

wobei zur Abkürzung die einsteinsche Summenkonvention benutzt wurde (die besagt, dass über gleiche Indizes summiert wird). In der noch weiter vereinfachten Feynman-Slash-Notation wird das Skalarprodukt mit den Gamma-Matrizen durch ein Slash-Symbol ausgedrückt. Es ergibt sich im Ortsraum

${\displaystyle {\Bigl [}\mathrm {i} \partial \!\!\!/\ -m{\Bigr ]}\,\psi (x)=0\,,}$

und im Impulsraum gilt

${\displaystyle {\Bigl [}p\!\!\!/\ -m{\Bigr ]}\,\psi (p)=0\,.}$

Wenn ${\displaystyle \psi (x)}$ die Dirac-Gleichung löst, dann löst auch der mit einer Phase ${\displaystyle \alpha }$ multiplizierte Spinor ${\displaystyle \mathrm {e} ^{\mathrm {i} q\alpha }\psi }$ die Dirac-Gleichung. Da alle physikalisch messbaren Größen mit jedem Faktor ${\displaystyle \psi }$ auch den konjugiert komplexen Faktor ${\displaystyle \psi ^{*}}$ enthalten, sind sie und die Dirac-Gleichung invariant unter dieser Phasentransformation des Dirac-Spinors ${\displaystyle \psi }$.

Bei nichtkonstantem ${\displaystyle \alpha }$ ergibt das eine zusätzliche U(1)-Eichinvarianz, und die partiellen Ableitungen müssen durch sog. kovariante Ableitungen ersetzt werden: Aus der Forderung der Invarianz unter allen Phasentransformationen, die stetig-differenzierbar von Zeit und Ort abhängen,

${\displaystyle \psi ^{\prime }(x)=\mathrm {e} ^{\mathrm {i} q\alpha (x)}\psi (x)\,,}$ ergibt sich die Notwendigkeit, die

partiellen Ableitungen in der Dirac-Gleichung durch die kovariante Ableitung zu ersetzen:

${\displaystyle D_{\mu }={\frac {\partial }{\partial x^{\mu }}}+\mathrm {i} qA_{\mu }\;.}$

Die hier auftretenden vier Funktionen ${\displaystyle A_{\mu }}$ bilden in der Physik das sog. Viererpotential oder Eichfeld. Mathematisch handelt es sich um eine Konnexion oder einen Zusammenhang. Definiert man das transformierte Eichfeld durch

${\displaystyle A_{\mu }^{\prime }=A_{\mu }-{\frac {\partial \alpha (x)}{\partial x^{\mu }}}\,,}$

dann löst ${\displaystyle \psi }$ die Dirac-Gleichung mit dem Eichfeld ${\displaystyle A_{\mu }}$

${\displaystyle \left(\mathrm {i} \sum _{\mu =0}^{3}\gamma ^{\mu }{\frac {\partial }{\partial x^{\mu }}}-q\sum _{\mu =0}^{3}\gamma ^{\mu }A_{\mu }-m\right)\psi (x)=0}$

oder in Slash-Notation

${\displaystyle \left(\mathrm {i} \partial \!\!\!/\,-qA\!\!\!/\,-m\right)\psi (x)=0}$

genau dann, wenn der transformierte Dirac-Spinor die Dirac-Gleichung mit dem transformierten Eichfeld erfüllt. Transformationen, deren Parameter so wie hier die Phase ${\displaystyle \alpha (x)}$ beliebig von Zeit und Ort abhängen dürfen, heißen in der Physik lokale Eichtransformationen.

Bei dem Eichfeld handelt es sich um das skalare Potential ${\displaystyle \phi }$ und das Vektorpotential ${\displaystyle {\vec {A}}}$ der Elektrodynamik,

${\displaystyle (A_{0},A_{1},A_{2},A_{3})=(\phi ,-A_{x},-A_{y},-A_{z})\,.}$

Wenn man sie wie angegeben transformiert, bleiben die elektrische und magnetische Feldstärke

${\displaystyle {\vec {E}}=-{\vec {\nabla }}\phi -{\frac {\partial {\vec {A}}}{\partial t}}\quad {\mbox{sowie}}\quad {\vec {B}}={\vec {\nabla }}\times {\vec {A}}}$

und alle anderen messbaren Größen unverändert.

Die Dirac-Gleichung mit kovarianter Ableitung und die Elektrodynamik sind invariant unter beliebigen zeit- und ortsabhängigen Transformationen der Phase des Dirac-Spinors. Der Parameter ${\displaystyle q}$ in der kovarianten Ableitung bestimmt die Stärke der Ankopplung der elektromagnetischen Potentiale an den Dirac-Spinor. Er entspricht dabei genau der elektrischen Ladung des Teilchens.

Die Ersetzung der partiellen Ableitungen in der Dirac-Gleichung durch eine kovariante Ableitung koppelt die elektromagnetischen Potentiale an den Dirac-Spinor. Man spricht dabei von sog. minimaler Kopplung im Gegensatz zu einem Kopplungsterm wie „magnetische Feldstärke mal Dirac-Spinor“, der auch eichinvariant wäre, aber nicht zur Ergänzung einer Ableitung zu einer kovarianten Ableitung erforderlich ist.

Nach Multiplikation mit ${\displaystyle \gamma ^{0}}$ kann man wegen ${\displaystyle (\gamma ^{0})^{2}=1}$ in der Dirac-Gleichung nach der Zeitableitung auflösen und die Dirac-Gleichung in die Form einer Schrödinger-Gleichung bringen,

${\displaystyle \mathrm {i} {\frac {\partial }{\partial t}}\,\psi =H_{\text{Dirac}}\psi }$

${\displaystyle H_{\text{Dirac}}=\sum _{k=1}^{3}\alpha ^{k}\left(-\mathrm {i} {\frac {\partial }{\partial x^{k}}}-qA_{k}\right)+\gamma ^{0}m+q\phi \;.}$

Die hier auftretenden 4×4-Matrizen, die leicht von den entsprechenden ${\displaystyle \gamma }$-Matrizen verschieden sind, ${\displaystyle \alpha ^{k}=\gamma ^{0}\,\gamma ^{k},}$ lassen sich ebenfalls kompakt mit Hilfe der Pauli-Matrizen durch Blöcke von 2×2-Matrizen ${\displaystyle \sigma _{2\times 2}^{k}}$ beschreiben:

${\displaystyle \alpha ^{k}={\begin{pmatrix}&\sigma _{2\times 2}^{k}\\\sigma _{2\times 2}^{k}&\end{pmatrix}}\,,\,k\in \{1,2,3\}\,,\quad \gamma ^{0}={\begin{pmatrix}1_{2\times 2}&\\&-1_{2\times 2}\end{pmatrix}}\,.}$

Der Differentialoperator auf der rechten Seite der Schrödinger-Gleichung ist der zur Dirac-Gleichung gehörige Hamiltonoperator ${\displaystyle H_{\text{Dirac}}\,.}$ Die möglichen Energien des Teilchens sind Eigenwerte dieses Hamiltonoperators.

Dabei zeigt die mathematische Untersuchung im Fall eines ungeladenen Teilchens (${\displaystyle q=0}$), dass das Spektrum positive und negative Werte enthält, ebenso wie man aus der Energie-Impuls-Relation der Klein-Gordon-Gleichung ${\displaystyle E^{2}-{\vec {p}}^{\,2}=m^{2}}$ (in natürlichen Maßeinheiten mit ${\displaystyle c=1=\hbar }$) die positiven und negativen Energiewerte ${\displaystyle E=\pm {\sqrt {m^{2}+{\vec {p}}^{\,2}}}}$ erhält.

Da Teilchen mit negativer Energie nie beobachtet wurden und da eine Welt mit Teilchen, deren Energien nach oben und nach unten unbeschränkt ist, instabil wäre, postulierte Dirac, dass das Vakuum ein Dirac-See sei, in dem jeder denkbare Zustand negativer Energie schon besetzt sei, sodass weitere Elektronen nur positive Energien annehmen könnten. Füge man diesem Dirac-See genügend Energie, mindestens die Ruheenergie zweier Elektronen, hinzu, so könne man einem See-Elektron positive Energie verleihen und das entstehende Loch verhielte sich wie ein Zustand mit der restlichen, ebenfalls positiven Energie und der fehlenden, entgegengesetzten Ladung. So sagte Dirac die Existenz von Antiteilchen und die Paarerzeugung von Elektron-Positron-Paaren voraus, die ein Jahr später beobachtet wurden.

Die Vorstellung eines Dirac-Sees gilt allerdings heute als unhaltbar[3] und ist durch die Feynman-Stückelberg-Interpretation ersetzt. Sie deutet die Dirac-Gleichung als Gleichung für ein Quantenfeld ${\displaystyle \psi (x)}$, das ist mathematisch ein Operator, der in den quantenmechanischen Zuständen Teilchen oder Antiteilchen erzeugt oder vernichtet. Die Erzeugung und Vernichtung von Teilchen während der Wechselwirkung des Elektrons mit dem Proton führt in der Quantenelektrodynamik zu einer kleinen Verschiebung der Energien verschiedener Zustände des Wasserstoffatoms, die ohne diese Erzeugungs- und Vernichtungsvorgänge gleiche Energie hätten. Die berechnete Größe dieser Lamb-Verschiebung stimmt innerhalb der Messgenauigkeit von sechs Stellen mit dem gemessenen Wert überein.

Die Erzeugung und Vernichtung von Teilchen während der Wechselwirkung des Elektrons mit einem Magnetfeld ändert auch den Dirac-Wert ${\displaystyle g=2}$ des gyromagnetischen Faktors. Sie bewirkt ein sogenanntes anomales magnetisches Moment, von dem man auch als g-2-Anomalie spricht. Der in der Quantenelektrodynamik berechnete Wert von ${\displaystyle g}$ stimmt mit dem gemessenen Wert auf zehn Dezimalstellen überein.

Ausgehend von der Schrödingerform der Dirac-Gleichung für ein Teilchen im elektromagnetischen Feld wird der Dirac-Spinor in zwei Zweierspinoren aufgespalten.

${\displaystyle \mathrm {i} \partial _{t}{\begin{pmatrix}\varphi _{1}\\\varphi _{2}\end{pmatrix}}={\begin{pmatrix}\sigma ^{k}\pi _{k}\varphi _{2}\\\sigma ^{k}\pi _{k}\varphi _{1}\end{pmatrix}}+q\phi {\begin{pmatrix}\varphi _{1}\\\varphi _{2}\end{pmatrix}}+m{\begin{pmatrix}\varphi _{1}\\-\varphi _{2}\end{pmatrix}}}$ mit ${\displaystyle \pi _{k}=-\mathrm {i} \partial _{k}-qA_{k}\;.}$

Unter der Annahme, dass sich das Teilchen nur langsam bewegt, sodass seine Energie nur wenig größer als seine Ruheenergie ist, kann die schnelle Zeitentwicklung, die von der Ruheenergie herrührt, abgespalten werden:

${\displaystyle {\begin{pmatrix}\varphi _{1}\\\varphi _{2}\end{pmatrix}}=\mathrm {e} ^{-\mathrm {i} mt}{\begin{pmatrix}\varphi \\\chi \end{pmatrix}}\;.}$

Aus diesem Ansatz folgt:

${\displaystyle \mathrm {i} \partial _{t}{\begin{pmatrix}\varphi \\\chi \end{pmatrix}}={\begin{pmatrix}\sigma ^{k}\pi _{k}\chi \\\sigma ^{k}\pi _{k}\varphi \end{pmatrix}}+q\phi {\begin{pmatrix}\varphi \\\chi \end{pmatrix}}+{\begin{pmatrix}0\\-2m\chi \end{pmatrix}}\,.}$

In der zweiten Zeile sind nach Annahme sowohl die Zeitableitung als auch die kinetischen Energien und die elektrostatische Energie klein gegenüber der Ruheenergie ${\displaystyle m}$. Daher ist ${\displaystyle \chi }$ klein gegen ${\displaystyle \varphi }$ und ungefähr gleich

${\displaystyle \chi \approx {\frac {\sigma ^{k}\pi _{k}}{2m}}\varphi }$.

In die erste Zeile eingesetzt ergibt sich:

${\displaystyle \mathrm {i} \partial _{t}\varphi ={\frac {(\sigma ^{k}\pi _{k})^{2}}{2m}}\varphi +q\phi \varphi \;.}$

Für das Produkt der Pauli-Matrizen erhält man

${\displaystyle (\sigma ^{k}\pi _{k})^{2}=(\delta ^{ij}+\mathrm {i} \varepsilon ^{ijk}\sigma _{k})\pi _{i}\pi _{j}={\vec {\pi }}^{2}-q\,\sigma ^{k}B_{k}.}$

Der Spinor ${\displaystyle \varphi }$ genügt daher der Pauli-Gleichung mit dem nichtklassischen Wert ${\displaystyle g=2,}$

${\displaystyle \mathrm {i} \partial _{t}\varphi ={\frac {{\vec {\pi }}^{2}}{2m}}\varphi +q\phi \varphi -g{\frac {q}{2m}}S^{k}B_{k}\varphi .}$

Dabei sind ${\displaystyle S^{k}={\tfrac {\sigma ^{k}}{2}}}$ die Komponenten des Spin-Operators.

Im homogenen Magnetfeld gilt ${\displaystyle \phi =0,{\vec {A}}={\tfrac {1}{2}}{\vec {B}}\times {\vec {x}}}$ und mit ${\displaystyle p_{k}=-\mathrm {i} \partial _{k}}$

${\displaystyle ({\vec {p}}-q{\vec {A}})^{2}={\vec {p}}^{2}-q{\vec {L}}\cdot {\vec {B}},}$

wenn man Terme vernachlässigt, die quadratisch in ${\displaystyle {\vec {B}}}$ sind. Dann besagt die Pauli-Gleichung

${\displaystyle \mathrm {i} \partial _{t}\varphi ={\frac {{\vec {p}}^{2}}{2m}}\varphi -{\frac {q}{2m}}({\vec {L}}+g{\vec {S}})\cdot {\vec {B}}\varphi .}$

Das Magnetfeld koppelt folglich nicht nur an den Bahndrehimpuls ${\displaystyle {\vec {L}}}$ und trägt nicht nur ${\displaystyle -{\tfrac {q}{2m}}{\vec {L}}\cdot {\vec {B}}}$ zur Energie bei. Der Faktor ${\displaystyle {\tfrac {q}{2m}}}$ ist das Magneton des Teilchens. In Drehimpulseigenzuständen ist ${\displaystyle {\vec {L}}\cdot {\vec {B}}}$ ein ganzzahliges Vielfaches der Magnetfeldstärke ${\displaystyle |{\vec {B}}|}$. Dagegen ergibt ${\displaystyle {\vec {S}}\cdot {\vec {B}}}$ ein halbzahliges Vielfaches, das erst nach Multiplikation mit ${\displaystyle g}$ ganzzahlig wird.[4]

Die Dirac-Gleichung bildet (nach Quantisierung des zugehörigen klassischen Feldes)[5] die Grundlage der relativistischen Quantenfeldtheorien der Hochenergiephysik. Erst seit wenigen Jahren[6] weiß man, dass auch bei nichtrelativistischen Energien Realisierungen existieren, nämlich bei Graphenen, das sind Schichtsysteme, die mit Graphit zusammenhängen. Und zwar braucht man hier nur den Grenzwert verschwindender Masse (sog. chiraler Limes) ${\displaystyle m=0}$ zu betrachten, und es ist zusätzlich die Lichtgeschwindigkeit ${\displaystyle c}$ durch die Grenzgeschwindigkeit ${\displaystyle v_{F}}$ des Elektronensystems, die sog. Fermi-Geschwindigkeit zu ersetzen. Als Konsequenz sind bei diesem System Energie ${\displaystyle E}$ und Impuls ${\displaystyle p}$ zueinander proportional (${\displaystyle E\propto p}$), während sonst bei nichtrelativistischen Elektronen ${\displaystyle E\propto p^{2}}$ gilt. Darüber hinaus ergeben sich zahlreiche weitere Besonderheiten.[6]

• Nichtlineare Dirac-Gleichung
• Dirac–Kähler-Gleichung

1. PAM Dirac: The Quantum Theory of the Electron. In: Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character. A, Nr. 778, 1928, S. 610–624, doi:10.1098/rspa.1928.0023 (Online).
2. C. D. Anderson: The Positive Electron. In: Physical Review. Band 43, Nr. 6, 1933, S. 491–494, doi:10.1103/PhysRev.43.491 (Online).
3. J. Schwinger: A Report on Quantum Electrodynamics. In: The Physicist’s Conception of Nature. Reidel, Dordrecht 1973, S. 415.
4. Bei isolierten Atomen oder Ionen muss man den Gesamt-Bahndrehimpuls und den Gesamt-Spindrehimpuls des Atoms bzw. Ions zu einem Gesamtdrehimpuls J (= L+S) addieren und erhält den sog. Landé-Faktor g(L, S; J). Dieser ist 1 bei reinem Gesamt-Bahndrehimpuls und 2 bei reinem Gesamt-Spindrehimpuls und hat sonst von 1 und 2 verschiedene Werte. Wenn ferner die betroffenen Atome in einen Festkörper eingebaut sind, erhält man Zusatzbeiträge, die ${\displaystyle g}$ wesentlich verändern können. Der Ferro- und Paramagnetismus typischer Repräsentanten ferromagnetischer oder paramagnetischer fester Körper bzw. paramagnetischer Moleküle ist trotzdem meist überwiegend Spinmagnetismus, weil experimentell sehr oft ${\displaystyle g\sim 2}$ gemessen wird.
5. Oft spricht man von zweiter Quantisierung.
6. Siehe den Artikel Graphen.