Dirac-Gleichung

Die Dirac-Gleichung ist eine grundlegende Gleichung der relativistischen Quantenmechanik. Sie beschreibt die Eigenschaften und das Verhalten eines fundamentalen Fermions mit Spin 1/2 (zum Beispiel Elektron, Quark). Sie wurde 1928 von Paul Dirac entwickelt und erfüllt im Gegensatz zur Schrödingergleichung die Anforderungen der speziellen Relativitätstheorie.

Die Dirac-Gleichung ist eine partielle Differentialgleichung erster Ordnung sowohl in den drei Raumkoordinaten als auch in der Zeit, im Einklang mit der von der speziellen Relativitätstheorie geforderten Invarianz unter Lorentz-Transformationen. Im nichtrelativistischen Grenzfall (${\tfrac {v}{c}}\to 0$) geht sie in die Pauli-Gleichung über, die im Gegensatz zur Schrödingergleichung noch die Spin-Bahn-Kopplung und weitere Terme enthält. Jede Lösung der Dirac-Gleichung entspricht einem möglichen Zustand des betreffenden Teilchens, mit der Besonderheit, dass zur Darstellung dieses Zustands vier räumliche Wellenfunktionen nötig sind (s. Dirac-Spinor), statt zwei in der nichtrelativistischen Theorie mit Spin oder einer einzigen im Fall von spinlosen Teilchen. Für die von der Dirac-Gleichung beschriebenen Teilchen gilt:

• Für ein freies Teilchen ist die relativistische Energie-Impuls-Beziehung $E={\sqrt {m^{2}c^{4} p^{2}c^{2}}}$ erfüllt.

• Für ein Teilchen im elektrostatischen Feld einer Punktladung ergibt sich das Wasserstoffspektrum mit seiner Feinstruktur.

• Das Teilchen hat einen Eigendrehimpuls (Spin), der die Quantenzahl 1/2 hat und – weil dies in der klassischen Physik nicht vorkommt – nicht wie bei einem Kreisel auf die Rotation einer Massenverteilung zurückgehen kann.

• Trägt das Teilchen eine elektrische Ladung, so ist mit dem Spin stets auch ein magnetisches Dipolmoment verknüpft (→ Spinmagnetismus). Im Vergleich mit dem magnetischen Dipol, den das Teilchen durch eine Rotationsbewegung bei gleich großem Drehimpuls hervorrufen würde (→ orbitaler Magnetismus), hat das mit dem Spin verbundene Moment die doppelte Stärke (s. Anomales magnetisches Moment des Elektrons).

• Zu dem Teilchen existiert ein Antiteilchen (zum Elektron also ein sog. Positron) mit derselben Masse und demselben Spin, aber mit entgegengesetzter Ladung und magnetischem Moment.

Alle genannten Eigenschaften entsprechen den experimentellen Befunden. Zur Zeit der Entdeckung der Dirac-Gleichung 1928 waren die vier erstgenannten schon bekannt, nicht aber ihre gemeinsame Grundlage. Die letztgenannte Eigenschaft wurde durch die Dirac-Gleichung vorhergesagt, und der erste Nachweis eines Antiteilchens gelang 1932 Carl David Anderson (s. Positron).

Der in der Diracgleichung vorkommende Differentialoperator spielt auch in der Mathematik (Differentialgeometrie) eine große Rolle (Dirac-Operator).

Dirac-Gleichung eines ungeladenen Teilchens

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

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

$\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 $\gamma ^{0},\gamma ^{1},\gamma ^{2}$ und $\gamma ^{3}$ wirken im Raum der vier Komponenten des Spinors und koppeln sie aneinander. Die Produkte von zwei Gamma-Matrizen haben die folgenden Eigenschaften:

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

$\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

$\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 $\psi$ die Klein-Gordon-Gleichung:

${\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) $-E^{2} {\vec {p}}^{\,2} m^{2}=0$, also die relativistische Energie-Impuls-Beziehung eines Teilchens der Masse $m\,.$

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

${\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 $\gamma _{0}$ bilden also die zweikomponentige Einheitsmatrix, die beiden letzten Komponenten deren Negatives. Analog ergeben die beiden oberen Komponenten der zweiten, dritten bzw. vierten $\gamma$-Matrix die drei 2×2-Pauli-Matrizen $\sigma _{j}$ und die beiden letzten Komponenten von $\gamma _{j}$ deren Negatives. Letztere gehen im nichtrelativistischen Grenzfall wie ${\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

$\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 $\gamma ^{\mu }$ und ${\frac {\partial }{\partial x^{\mu }}}$. Dieses ist bei Lorentztransformation jedoch nicht invariant, denn $\gamma ^{\mu }$ bleibt konstant. Die Lorentzinvarianz der Dirac-Theorie ergibt sich erst dadurch, dass der Dirac-Operator auf einen Spinor $\psi$ wirkt, dessen vier Komponenten geeignet mittransformiert werden. Im Endergebnis geht damit eine Lösung $\psi (x)$ der Dirac-Gleichung durch Lorentztransformation in eine Lösung der entsprechend transformierten Dirac-Gleichung über.

Impulsraum und Slash-Notation

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

${\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

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

und im Impulsraum gilt

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

Eichinvarianz und elektromagnetische Wechselwirkung

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

Bei nichtkonstantem $\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, bei denen die Phase $\alpha (x)$ stetig-differenzierbar von Zeit und Ort abhängt,

$\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:

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

Die hier auftretenden vier Funktionen $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

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

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

$\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

$\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 $\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 $\phi$ und das Vektorpotential ${\vec {A}}$ der Elektrodynamik,

$(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

${\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 $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.

Schrödingerform

Nach Multiplikation mit $\gamma ^{0}$ kann man wegen $(\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,

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

$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 $\gamma$-Matrizen verschieden sind, $\alpha ^{k}=\gamma ^{0}\,\gamma ^{k},$ lassen sich ebenfalls kompakt mit Hilfe der Pauli-Matrizen durch Blöcke von 2×2-Matrizen $\sigma _{2\times 2}^{k}$ beschreiben:

$\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 $H_{\text{Dirac}}\,.$ Die möglichen Energien des Teilchens sind Eigenwerte dieses Hamiltonoperators.

Dabei zeigt die mathematische Untersuchung im Fall eines ungeladenen Teilchens ($q=0$), dass das Spektrum positive und negative Werte enthält, ebenso wie man aus der Energie-Impuls-Relation der Klein-Gordon-Gleichung $E^{2}-{\vec {p}}^{\,2}=m^{2}$ (in natürlichen Maßeinheiten mit $c=1=\hbar$) die positiven und negativen Energiewerte $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 und ist durch die Feynman-Stückelberg-Interpretation ersetzt. Sie deutet die Dirac-Gleichung als Gleichung für ein Quantenfeld $\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 $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 $g$ stimmt mit dem gemessenen Wert auf zehn Dezimalstellen überein.

Herleitung des gyromagnetischen Faktors

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

$\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 $\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:

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

Aus diesem Ansatz folgt:

$\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 $m$. Daher ist $\chi$ klein gegen $\varphi$ und ungefähr gleich

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

In die erste Zeile eingesetzt ergibt sich:

$\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

$(\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 $\varphi$ genügt daher der Pauli-Gleichung mit dem nichtklassischen Wert $g=2,$

$\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 $S^{k}={\tfrac {\sigma ^{k}}{2}}$ die Komponenten des Spin-Operators.

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

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

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

$\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 ${\vec {L}}$ und trägt nicht nur $-{\tfrac {q}{2m}}{\vec {L}}\cdot {\vec {B}}$ zur Energie bei. Der Faktor ${\tfrac {q}{2m}}$ ist das Magneton des Teilchens. In Drehimpulseigenzuständen ist ${\vec {L}}\cdot {\vec {B}}$ ein ganzzahliges Vielfaches der Magnetfeldstärke $|{\vec {B}}|$. Dagegen ergibt ${\vec {S}}\cdot {\vec {B}}$ ein halbzahliges Vielfaches, das erst nach Multiplikation mit $g$ ganzzahlig wird.

Realisierungen in Hochenergie- und Festkörperphysik

Die Dirac-Gleichung bildet (nach Quantisierung des zugehörigen klassischen Feldes) die Grundlage der relativistischen Quantenfeldtheorien der Hochenergiephysik. Erst seit wenigen Jahren 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) $m=0$ zu betrachten, und es ist zusätzlich die Lichtgeschwindigkeit $c$ durch die Grenzgeschwindigkeit $v_{F}$ des Elektronensystems, die sog. Fermi-Geschwindigkeit zu ersetzen. Als Konsequenz sind bei diesem System Energie $E$ und Impuls $p$ zueinander proportional ($E\propto p$), während sonst bei nichtrelativistischen Elektronen $E\propto p^{2}$ gilt. Darüber hinaus ergeben sich zahlreiche weitere Besonderheiten.

Siehe auch

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

Literatur

Artikel

• P. A. M. Dirac: The Quantum Theory of the Electron. In: Proceedings of the Royal Society of London. Series A. Band 117, Nr. 778, 1. Januar 1928, S. 610–624, doi:10.1098/rspa.1928.0023. 
• P. A. M. Dirac: The Quantum Theory of the Electron. Part II. In: Royal Society of London Proceedings Series A. Band 118, 1. Februar 1928, S. 351–361. 
• P. A. M. Dirac: A Theory of Electrons and Protons. In: Proceedings of the Royal Society of London. Series A, Containing Papers of a Mathematical and Physical Character. Band 126, Nr. 801, 1930, S. 360–365, JSTOR:95359. 
• Carl D. Anderson: The Positive Electron. In: Physical Review. Band 43, Nr. 6, 15. Februar 1933, S. 491–494, doi:10.1103/PhysRev.43.491. 

Bücher

• James Bjorken, Sidney Drell: Relativistische Quantenmechanik. Mannheim, Bibliographisches Institut, 1990. (BI Hochschultaschenbücher; 98/98a), ISBN 3-411-00098-8.
  Engl. Originalausgabe: Relativistic Quantum Mechanics. McGraw-Hill, New York 1964, ISBN 0-07-005493-2.
• James Bjorken, Sidney Drell: Relativistische Quantenfeldtheorie. (Dt. Übers.: J. Benecke, D. Maison, E. Riedel).
  Unveränd. Nachdruck: Mannheim, Zürich. BI-Wissenschaftsverlag, 1993. BI-Hochschultaschenbuch; 101, ISBN 3-411-00101-1.
  Engl. Originalausgabe: Relativistic Quantum Fields. McGraw-Hill, New York 1965, ISBN 0-07-005494-0.
• R. P. Feynman: Quantenelektrodynamik. 4. Auflage, ISBN 3-486-24337-3.
• Walter Greiner: Relativistische Quantenmechanik. Wellengleichungen. Band 6, ISBN 3-8171-1022-7.
• Franz Schwabl: Quantenmechanik für Fortgeschrittene (QM II). ISBN 978-3-540-25904-6.