Maxwell-Gleichungen in der klassischen Elektrodynamik

In vorangegangenen Kapiteln wurden bereits die wesentlichsten Resultate der Elektrostatik und Elektrodynamik vorgestellt. Aus diesen kristallisierten sich insbesondere die sog. »Maxwell'schen Gleichungen« der Elektrodynamik heraus. Im Folgenden soll der umgekehrte Weg beschritten werden: Ausgehend von den Maxwell-Gleichungen wird gezeigt, wie sich aus ihnen die Gesetze der klassischen Elektrodynamik und Elektrostatik ergeben. Der Begriff »klassisch« bedeute hier »nicht quantenmechanisch«. Auf eine (speziell-) relativistische Formulierung der Maxwell-Gleichungen wird hingegen eingegangen.

Mikroskopische Maxwellgleichungen

Die Maxwell-Gleichungen geben Beziehungen zwischen folgenden physikalischen Größen an, die wir als bekannt voraussetzen, da sie bereits in den früheren Kapiteln über Elektrodynamik und Elektrostatik vorgestellt und diskutiert wurden:

• Ladung q: Neben der Masse m kann ein mechanischer Körper auch

eine Ladung besitzen. Ladungen können im Gegensatz zu Massen auch negativ sein;

• eine Ladungsdichte $\rho$ ist von der Dimension

»Ladung pro Volumen«, d.h. $q=\int d^{3}x\rho$, wobei wir hier über das Volumen des Körpers integriert haben;

• die Stromdichte $\overrightarrow{J}$ gibt Auskunft über die

Ladungsmenge, die durch eine Fläche je Zeiteinheit fließt, und hängt daher mit dem sog. (elektrischen) Strom I zusammen: $I=\int d^2\overrightarrow{a}\cdot \overrightarrow{J}$, wobei wir hier über z.B. die Querschnittsfläche eines Leiters integriert haben;

• In einem elektrisches Feld $\overrightarrow{E}$ erfahren geladene

Körper Kräfte. Diese Körper können dabei auch durchaus ruhen;

• In einem magnetischen Feld erfahren bewegte geladene Körper Kräfte.

In Abwesenheit von sog. »Dielektrika« oder »Magnetika« brauchen wir nicht zwischen der magnetische Induktion/ Flussdichte $\overrightarrow{B}$ und dem Magnetfeld (später mit $\overrightarrow{H}$ bezeichnet) unterscheiden;

Die Lichtgeschwindigkeit im Vakuum werde zudem mit c bezeichnet.

Wenn keine dielektrische oder (z.B. dia-, para-, ferro-)magnetische Materialien vorhanden sind, lauten die somit »mikroskopischen« Maxwell-Gleichungen (in ihrer differenziellen Formulierung):

• $\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{B}=0$: Diese Gleichung sagt

u.a. etwas darüber aus, dass es keine magnetischen Monopole gibt, während die folgende Gleichung zeigt, dass elektrische Monopole hingegen möglich sind.

• $\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{E}=4\pi \rho$: Hieraus lässt

sich nämlich u.a. das bekannte Coulomb-Gesetz der Elektrostatik herleiten. Eine Ladungsdichte ist also die Quelle eines elektrischen Feldes.

• $\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{B}=\frac{4\pi }{c}\overrightarrow{J}+\frac{1}{c}{\partial }_t\overrightarrow{E}$:

Diese Gleichung enthält das gleichermaßen bereits als bekannt vorausgesetzte Ampère'sches Durchflutungsgesetz und ist zudem noch um den sog. Maxwell'schen Verschiebungsstrom-Term erweitert. D.h. sowohl Ströme (bzw. wie hier: Stromdichten) als auch ein sich zeitlich änderndes elektrisches Feld bewirken ein magnetisches Wirbelfeld.

• $\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E}=-\frac{1}{c}{\partial }_t\overrightarrow{B}$:

Jedes sich ändernde Magnetfeld erzeugt ein elektrisches Wirbelfeld. Diese Gleichung beschreibt also Faradays Induktionsgesetz. Im Minuszeichen vor der (partiellen) Zeitableitung der magnetischen Induktion liegt die Ursache der sog. Lentz'schen Regel.

Makroskopische Maxwellgleichungen

Bei Anwesenheit von Dielektrika oder Magnetika werden vom elektrischen Feld $\overrightarrow{E}$ und von der magnetischen Flussdichte $\overrightarrow{B}$ noch die elektrische Flussdichte $\overrightarrow{D}$ bzw. das magnetisches Feld $\overrightarrow{H}$ unterschieden. Die somit makroskopischen Maxwell'schen Gleichungen lauten dann:

• $\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{B}=0$ (keine magnet. Monopole),

• $\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{D}=4\pi \rho$ (Coulomb-Gesetz),

• $\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{H}=\frac{4\pi }{c}\overrightarrow{J}+\frac{1}{c}{\partial }_t\overrightarrow{D}$

(Ampère'sches Durchflutungsgesetz, Verschiebungsstrom),

• $\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E}=-\frac{1}{c}{\partial }_t\overrightarrow{B}$

(Induktionsgesetz; Minuszeichen: Lentz'sche Regel).

Die elektrische Flussdichte hängt dabei mit dem elektrischen Feld über die elektrische Polarisation $\overrightarrow{P}$ zusammen:

$\overrightarrow{D}=\overrightarrow{E}+4\pi \overrightarrow{P}$.

Die magnetische Flussdichte wird mit dem magnetischen Feld über die magnetische Polarisation $\overrightarrow{M}$ verknüpft:

$\overrightarrow{B}=\overrightarrow{H}+4\pi \overrightarrow{M}$.

Man bezeichnet manchmal die Dielektrika bzw. Magnetika auch als polarisierbare Medien.

In sog. isotropen Dielektrika gilt

$\overrightarrow{D}=\epsilon \overrightarrow{E}\Rightarrow \overrightarrow{P}=\frac{1}{4\pi }(\epsilon -1)\overrightarrow{E}={\chi }_E\overrightarrow{E}$

und in isotropen Magnetika gilt entsprechend:

$\overrightarrow{B}=\mu \overrightarrow{H}\Rightarrow \overrightarrow{M}=\frac{1}{4\pi }(\mu -1)\overrightarrow{H}={\chi }_M\overrightarrow{H}$,

wobei wir zudem die elektrische Suszeptibilität

${\chi }_E=\frac{1}{4\pi }(\epsilon -1)$,

und die magnetische Suszeptibilität

${\chi }_M=\frac{1}{4\pi }(\mu -1)$

eingeführt haben.

Im nicht isotropen Fall existiert somit folgender Zusammenhang zwischen den mikroskopischen und den makroskopischen Maxwellgleichungen:

$\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{E}=4\pi \tilde{\rho }$, $\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{B}=\frac{4\pi }{c}\tilde{\overrightarrow{J}}+\frac{1}{c}{\partial }_t\overrightarrow{E}$ mit $\tilde{\rho }=\rho -\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{P}$, $\tilde{\overrightarrow{J}}=\overrightarrow{J}+c\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{M}+{\partial }_t\overrightarrow{P}$ .

Im isotropen Fall erhalten wir für diesen hingegen:

$\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{E}=4\pi \tilde{\rho }$, $\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{B}=\frac{4\pi \epsilon \mu }{c}\tilde{\overrightarrow{J}}+\frac{\epsilon \mu }{c}{\partial }_t\overrightarrow{E}$ mit $\tilde{\rho }=\frac{1}{\epsilon }\rho$, $\tilde{\overrightarrow{J}}=\frac{1}{\epsilon }\overrightarrow{J}$, der Lichtgeschwindigkeit im Medium $c_n=\frac{c}{n}$ und dem Brechungsindex $n=\sqrt{\epsilon \mu }$ .

Die homogenen Maxwellgleichungen, d.h. $\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{B}=0$ und $\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E}+\frac{1}{c}{\partial }_t\overrightarrow{B}=0$ (also jene ohne Quellterme wie Ladungs- bzw. Stromdichten), bleiben hingegen bei Anwesenheit polarisierbarer Medien unverändert.

Induktionsgesetz

Mit Hilfe des Stokes'schen Integralsatzes können wir die Maxwell-Gleichung $\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E}=-\frac{1}{c}{\partial }_t\overrightarrow{B}$ in eine Integralgleichung umwandeln:

$\underset{C=\partial A}{\int }d\overrightarrow{x}\cdot \overrightarrow{E}=\underset{A}{\int }d\overrightarrow{a}\cdot (\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E})=-\frac{1}{c}\underset{A}{\int }d\overrightarrow{a}\cdot {\partial }_t\overrightarrow{B}=-\frac{1}{c}\underset{A}{\frac{d}{dt}\int }d\overrightarrow{a}\cdot \overrightarrow{B}+\frac{1}{c}\underset{A}{\int }d\overrightarrow{a}\cdot \left[(\overrightarrow{v}\cdot \overrightarrow{\mathrm{\nabla }})\overrightarrow{B}\right]$,

weil $\frac{d}{dt}\overrightarrow{B}(t,\overrightarrow{x}\left(t\right))={\partial }_t\overrightarrow{B}+(\overrightarrow{v}\cdot \overrightarrow{\mathrm{\nabla }})\overrightarrow{B}$ mit $\overrightarrow{v}=\dot{\overrightarrow{x}}$ gilt. Mit den folgenden weiteren Umformungen (über gleiche Indizes werde dabei von 1 bis 3 summiert)

${[\overrightarrow{\mathrm{\nabla }}\times (\overrightarrow{v}\times \overrightarrow{B})]}_i={\epsilon }_{ikl}{\partial }_k{[\overrightarrow{v}\times \overrightarrow{B}]}_l={\epsilon }_{ikl}{\epsilon }_{lmn}{\partial }_k{v}_m{B}_n=({\delta }_{im}{\delta }_{kn}-{\delta }_{in}{\delta }_{kn}){\partial }_k{v}_m{B}_n=\overrightarrow{\mathrm{\nabla }}\cdot \left(v_i\overrightarrow{B}\right)-(\overrightarrow{v}\cdot \overrightarrow{\mathrm{\nabla }})B_i=-(\overrightarrow{v}\cdot \overrightarrow{\mathrm{\nabla }})B_i$,

weil $\overrightarrow{\mathrm{\nabla }}\cdot \left(v_i\overrightarrow{B}\right)=\left(\overrightarrow{\mathrm{\nabla }}{v}_i\right)\cdot \overrightarrow{B}+v_i(\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{B})=0$ mit $\overrightarrow{\mathrm{\nabla }}{v}_i=0$ und $\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{B}=0$ gilt,

erhalten wir schließlich durch erneutes Anwenden des Stokes'schen Satzes das Induktionsgesetz in Integralform:

$\underset{C=\partial A}{\int }d\overrightarrow{x}\cdot \overrightarrow{E}=-\frac{1}{c}\underset{A}{\frac{d}{dt}\int }d\overrightarrow{a}\cdot \overrightarrow{B}-\frac{1}{c}\underset{A}{\int }d\overrightarrow{a}\cdot [\overrightarrow{\mathrm{\nabla }}\times (\overrightarrow{v}\times \overrightarrow{B})]=-\frac{1}{c}\underset{A}{\frac{d}{dt}\int }d\overrightarrow{a}\cdot \overrightarrow{B}-\frac{1}{c}\underset{C=\partial A}{\int }d\overrightarrow{x}\cdot (\overrightarrow{v}\times \overrightarrow{B})$.

Hierin führen wir die Induktionsspannung U und den magnetische Fluss ${\mathrm{\Phi }}_M$ ein:

$U=\underset{C=\partial A}{\int }d\overrightarrow{x}\cdot \overrightarrow{\overline{E}}=\underset{C=\partial A}{\int }d\overrightarrow{x}\cdot (\overrightarrow{E}+\frac{1}{c}(\overrightarrow{v}\times \overrightarrow{B}))=-\frac{1}{c}\underset{A}{\frac{d}{dt}\int }d\overrightarrow{a}\cdot \overrightarrow{B}=-\frac{1}{c}\frac{d}{dt}{\mathrm{\Phi }}_M$.

Multiplizieren wir das dort auftretende elektrische Feld $\overrightarrow{\overline{E}}$ mit einer Ladung q, so ergibt sich offensichtlich die bekannte Lorentzkraft auf eine mit der Geschwindigkeit $\overrightarrow{v}$ bewegten Ladung:

$\overrightarrow{F}=q\overrightarrow{\overline{E}}=q\overrightarrow{E}+\frac{q}{c}(\overrightarrow{v}\times \overrightarrow{B})$.

Wellengleichung für die elektrische Feldstärke

Für die elektrische Feldstärke lässt sich aus den mikroskopischen Maxwell-Gleichungen folgendermaßen eine Wellengleichung gewinnen:

Durch Einsetzen von

$\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E}=-\frac{1}{c}{\partial }_t\overrightarrow{B}$

in jene inhomogene Maxwell-Gleichung mit der Stromdichte,

$\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{B}=\frac{4\pi }{c}\overrightarrow{J}+\frac{1}{c}{\partial }_t\overrightarrow{E}$,

die zusätzlich noch partiell nach der Zeit abgeleitet wurde,

$-\frac{1}{c}{\partial }_t\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{B}=-\frac{4\pi }{c^2}{\partial }_t\overrightarrow{J}-\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}\overrightarrow{E}$,

d.h.

$-\frac{4\pi }{c^2}{\partial }_t\overrightarrow{J}-\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}\overrightarrow{E}=-\frac{1}{c}{\partial }_t\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{B}=\overrightarrow{\mathrm{\nabla }}\times (\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E})=\overrightarrow{\mathrm{\nabla }}\underset{4\pi \rho }{\underbrace{(\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{E})}}-{\overrightarrow{\mathrm{\nabla }}}^2\overrightarrow{E}$,

weil

${[\overrightarrow{\mathrm{\nabla }}\times (\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E})]}_i={\epsilon }_{ikl}{\partial }_k{[\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E}]}_l={\epsilon }_{ikl}{\epsilon }_{lmn}{\partial }_k{\partial }_m{E}_n$ $=({\delta }_{im}{\delta }_{kn}-{\delta }_{in}{\delta }_{km}){\partial }_k{\partial }_m{E}_n={\partial }_i(\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{E})-{\overrightarrow{\mathrm{\nabla }}}^2{E}_i$.

Die sich somit ergebende Wellengleichung lautet:

$(\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}-{\overrightarrow{\mathrm{\nabla }}}^2)\overrightarrow{E}=-4\pi \overrightarrow{\mathrm{\nabla }}\rho -\frac{4\pi }{c^2}{\partial }_t\overrightarrow{J}$.

Wenn weder Ladungen noch Ströme vorhanden sind, gilt:

$(\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}-{\overrightarrow{\mathrm{\nabla }}}^2)\overrightarrow{E}=0$.

Eine allgemeine Lösung der letzteren Wellengleichung ist

$\overrightarrow{E}(t,\overrightarrow{x})=\int \frac{d^{3}k}{{\left(2\pi \right)}^3}({\overrightarrow{E}}_0{e}^{-i\omega t+i\overrightarrow{k}\cdot \overrightarrow{x}}+{\overrightarrow{E}}_0^{\star }{e}^{+i\omega t-i\overrightarrow{k}\cdot \overrightarrow{x}})$,

wobei $\omega =kc$. Die Lösung $\overrightarrow{E}$ ist reell, während die Amplitude ${\overrightarrow{E}}_0$ aber im Allg. komplex sein kann.

Wellengleichung für das magnetische Feld

Für die magnetische Flussdichte ergibt sich aus den mikroskopischen Maxwellgleichungen wie beim elektrischen Feld eine Wellengleichung:

$\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{B}=\frac{4\pi }{c}\overrightarrow{J}+\frac{1}{c}{\partial }_t\overrightarrow{E}\Rightarrow$ $\frac{4\pi }{c}\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{J}+\frac{1}{c}{\partial }_t\underset{-\frac{1}{c}{\partial }_t\overrightarrow{B}}{\underbrace{\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E}}}=\overrightarrow{\mathrm{\nabla }}\times (\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{B})=\overrightarrow{\mathrm{\nabla }}\underset{0}{\underbrace{(\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{B})}}-{\overrightarrow{\mathrm{\nabla }}}^2\overrightarrow{B}$.

Diese lautet also:

$(\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}-{\overrightarrow{\mathrm{\nabla }}}^2)\overrightarrow{B}=\frac{4\pi }{c}\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{J}$.

Kontinuitätsgleichung

Auch die Kontinuitätsgleichung ist nur eine Folge der Maxwell-Gleichungen. Wir bilden die Divergenz von folgender Gleichung:

$\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{B}=\frac{4\pi }{c}\overrightarrow{J}+\frac{1}{c}{\partial }_t\overrightarrow{E}\Rightarrow$ $0=\overrightarrow{\mathrm{\nabla }}\cdot (\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{B})=\frac{4\pi }{c}\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{J}+\frac{1}{c}{\partial }_t\underset{4\pi \rho }{\underbrace{\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{E}}}$.

D.h.

${\partial }_t\rho +\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{J}=0$.

Anm.: Hingegen ergibt sich durch Divergenzbildung von $\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E}=-\frac{1}{c}{\partial }_t\overrightarrow{B}$ mit $0=\overrightarrow{\mathrm{\nabla }}\cdot (\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E})=-\frac{1}{c}{\partial }_t\underset{=0}{\underbrace{\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{B}}}$ nichts Neues. Letzteres gilt, weil

$\overrightarrow{\mathrm{\nabla }}\cdot (\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E})={\epsilon }_{ikl}{\partial }_i{\partial }_k{E}_l=\frac{1}{2}({\epsilon }_{ikl}+{\epsilon }_{ikl}){\partial }_i{\partial }_k{E}_l=\frac{1}{2}({\epsilon }_{ikl}{\partial }_i{\partial }_k{E}_l+{\epsilon }_{kil}{\partial }_k{\partial }_i{E}_l)$

$=\frac{1}{2}({\epsilon }_{ikl}+{\epsilon }_{kil}){\partial }_i{\partial }_k{E}_l=\frac{1}{2}({\epsilon }_{ikl}-{\epsilon }_{ikl}){\partial }_i{\partial }_k{E}_l=0$.

Leistungsdichte der Joule'schen Wärme

Für eine Punktladung q mit Stromdichte $\overrightarrow{J}(t,\overrightarrow{x})=\rho \left(\overrightarrow{x}\right)\dot{\overrightarrow{x}}$ und Ladungsdichte $\rho \left(\overrightarrow{x}\right)=q{\delta }^3(\overrightarrow{x}-{\overrightarrow{x}}^{\prime })$ im elektrischen Feld $\overrightarrow{E}$ lässt sich die Leistung der Joule'schen Wärme wie folgt angeben:

$-\frac{dW}{dt}=\overrightarrow{F}\cdot \dot{\overrightarrow{x}}=q\overrightarrow{E}\cdot \dot{\overrightarrow{x}}=\overrightarrow{E}\cdot q\dot{\overrightarrow{x}}=\int d^3{x}^{\prime }\overrightarrow{E}(t,{\overrightarrow{x}}^{\prime })\cdot q{\delta }^3(\overrightarrow{x}-{\overrightarrow{x}}^{\prime }){\dot{\overrightarrow{x}}}^{\prime }$ $=\int d^3{x}^{\prime }\overrightarrow{E}(t,{\overrightarrow{x}}^{\prime })\cdot \rho \left({\overrightarrow{x}}^{\prime }\right){\dot{\overrightarrow{x}}}^{\prime }=\int d^3{x}^{\prime }\overrightarrow{E}(t,{\overrightarrow{x}}^{\prime })\cdot \overrightarrow{J}\left({\overrightarrow{x}}^{\prime }\right)$.

Verwenden wir folgenden Zusammenhang zwischen elektrischen Strom I' und Stromdichte $\overrightarrow{J}$,

$I\int d{\overrightarrow{x}}^{\prime }\cdot ...=\int d^3{x}^{\prime }\overrightarrow{J}\left({\overrightarrow{x}}^{\prime }\right)...$,

dann erhalten wir für die als Joule'sche Wärme erzeugte Leistung:

$-\frac{dW}{dt}=\int d^3{x}^{\prime }\overrightarrow{E}(t,{\overrightarrow{x}}^{\prime })\cdot \overrightarrow{J}\left({\overrightarrow{x}}^{\prime }\right)=I\int d{\overrightarrow{x}}^{\prime }\cdot \overrightarrow{E}(t,{\overrightarrow{x}}^{\prime })=IU$

mit der Spannung

$U=\int d{\overrightarrow{x}}^{\prime }\cdot \overrightarrow{E}(t,{\overrightarrow{x}}^{\prime })$.

Die Leistungsdichte der Joule'schen Wärme ist somit $-\overrightarrow{J}\cdot \overrightarrow{E}$.

Poynting-Vektor

Die Leistungsdichte der Joule'schen Wärme $-\overrightarrow{J}\cdot \overrightarrow{E}$ lässt sich mittels $\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{H}=\frac{4\pi }{c}\overrightarrow{J}+\frac{1}{c}{\partial }_t\overrightarrow{D}$ folgendermaßen angeben:

$-\overrightarrow{J}\cdot \overrightarrow{E}=-\frac{c}{4\pi }(\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{H})\cdot \overrightarrow{E}-\frac{1}{4\pi }\overrightarrow{E}\cdot {\partial }_t\overrightarrow{D}$.

Verwenden wir hierin

$\overrightarrow{\mathrm{\nabla }}\cdot (\overrightarrow{E}\times \overrightarrow{H})={\epsilon }_{jkl}{\partial }_j\left(E_k{H}_l\right)={\epsilon }_{jkl}{H}_l{\partial }_j{E}_k+{\epsilon }_{jkl}{E}_k{\partial }_j{H}_l=\overrightarrow{H}\cdot \underset{=-\frac{1}{c}{\partial }_t\overrightarrow{B}}{\underbrace{(\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E})}}-\overrightarrow{E}\cdot (\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{H})$,

dann resultiert daraus

$-\overrightarrow{J}\cdot \overrightarrow{E}=\frac{c}{4\pi }\overrightarrow{\mathrm{\nabla }}\cdot \underset{=:\frac{4\pi }{c}\overrightarrow{S}}{\underbrace{(\overrightarrow{E}\times \overrightarrow{H})}}\cdot \overrightarrow{E}+\frac{1}{4\pi }\underset{=:{\partial }_{t}w}{\underbrace{(\overrightarrow{E}\cdot {\partial }_t\overrightarrow{D}+\overrightarrow{H}\cdot {\partial }_t\overrightarrow{B})}}$,

d.h.

${\partial }_{t}w+\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{S}=-\overrightarrow{J}\cdot \overrightarrow{E}$

mit dem Poynting-Vektor $\overrightarrow{S}=\frac{c}{4\pi }(\overrightarrow{E}\times \overrightarrow{H})$ und der Energiedichte w.

Poynting-Vektor für Felder mit harmonischer Zeitabhängigkeit

Wir nehmen Felder mit harmonischer Zeitabhängigkeit $e^{-i\omega t}$ an:

$\overrightarrow{E}(t,\overrightarrow{x})=\mathrm{\Re }\left(\overrightarrow{E}\left(\overrightarrow{x}\right)e^{-i\omega t}\right)=\frac{1}{2}[\overrightarrow{E}\left(\overrightarrow{x}\right)e^{-i\omega t}+{\overrightarrow{E}}^{\star }\left(\overrightarrow{x}\right)e^{i\omega t}]$, $\overrightarrow{H}(t,\overrightarrow{x})=\mathrm{\Re }\left(\overrightarrow{H}\left(\overrightarrow{x}\right)e^{-i\omega t}\right)=\frac{1}{2}[\overrightarrow{H}\left(\overrightarrow{x}\right)e^{-i\omega t}+{\overrightarrow{H}}^{\star }\left(\overrightarrow{x}\right)e^{i\omega t}]$.

Mit diesen lässt sich der Poynting-Vektor angeben:

$\overrightarrow{S}(t,\overrightarrow{x})=\frac{c}{4\pi }\overrightarrow{E}(t,\overrightarrow{x})\times \overrightarrow{H}(t,\overrightarrow{x})$

$=\frac{c}{16\pi }[\overrightarrow{E}\left(\overrightarrow{x}\right)\times \overrightarrow{H}\left(\overrightarrow{x}\right)e^{-i2\omega t}+{\overrightarrow{E}}^{\star }\left(\overrightarrow{x}\right)\times {\overrightarrow{H}}^{\star }\left(\overrightarrow{x}\right)e^{i2\omega t}+{\overrightarrow{E}}^{\star }\left(\overrightarrow{x}\right)\times \overrightarrow{H}\left(\overrightarrow{x}\right)+\overrightarrow{E}\left(\overrightarrow{x}\right)\times {\overrightarrow{H}}^{\star }\left(\overrightarrow{x}\right)]$.

Hieraus resultiert der zeitliche Mittelwert des Poynting-Vektors:

$⟨\overrightarrow{S}(t,\overrightarrow{x})⟩=\underset{T\to \mathrm{\infty }}{\lim }\frac{1}{2T}\underset{-T}{\stackrel{T}{\int }}dt\overrightarrow{S}(t,\overrightarrow{x})=\frac{c}{16\pi }[{\overrightarrow{E}}^{\star }\left(\overrightarrow{x}\right)\times \overrightarrow{H}\left(\overrightarrow{x}\right)+\overrightarrow{E}\left(\overrightarrow{x}\right)\times {\overrightarrow{H}}^{\star }\left(\overrightarrow{x}\right)]=\frac{c}{8\pi }\mathrm{\Re }(\overrightarrow{E}\left(\overrightarrow{x}\right)\times {\overrightarrow{H}}^{\star }\left(\overrightarrow{x}\right))$,

weil

$\frac{1}{2T}\underset{-T}{\stackrel{T}{\int }}dt{e}^{\pm i2\omega t}=\pm \frac{1}{4i\omega T}(e^{\pm i2\omega T}-e^{\mp i2\omega T})=\frac{\sin 2\omega T}{2\omega T}\underset{T\to \mathrm{\infty }}{⟶}0$,

$\mathrm{\Re }\left(z\right)=\frac{1}{2}(z+z^{\star })$.

Wenn wir in

$\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E}=-\frac{1}{c}{\partial }_t\overrightarrow{B}$

die zuvor bereits präsentierte Darstellung von $\overrightarrow{E}$ und $\overrightarrow{B}$ mit $\overrightarrow{E}\left(\overrightarrow{x}\right)={\overrightarrow{E}}_0{e}^{i\overrightarrow{k}\cdot \overrightarrow{x}}$ bzw. $\overrightarrow{B}\left(\overrightarrow{x}\right)={\overrightarrow{B}}_0{e}^{i\overrightarrow{k}\cdot \overrightarrow{x}}$ einsetzen, wobei $\overrightarrow{B}=\mu \overrightarrow{H}$ gilt, dann erhalten wir:

$\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E}(t,\overrightarrow{x})=\frac{1}{2}[\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E}\left(\overrightarrow{x}\right)e^{-i\omega t}+\overrightarrow{\mathrm{\nabla }}\times {\overrightarrow{E}}^{\star }\left(\overrightarrow{x}\right)e^{i\omega t}]=\frac{1}{2}[\overrightarrow{\mathrm{\nabla }}\times {\overrightarrow{E}}_0{e}^{i\overrightarrow{k}\cdot \overrightarrow{x}}{e}^{-i\omega t}+\overrightarrow{\mathrm{\nabla }}\times {\overrightarrow{E}}_0^{\star }{e}^{-i\overrightarrow{k}\cdot \overrightarrow{x}}{e}^{i\omega t}]=$

$=\frac{1}{2}[i\overrightarrow{k}\times {\overrightarrow{E}}_0{e}^{i\overrightarrow{k}\cdot \overrightarrow{x}}{e}^{-i\omega t}-i\overrightarrow{k}\times {\overrightarrow{E}}_0^{\star }{e}^{-i\overrightarrow{k}\cdot \overrightarrow{x}}{e}^{i\omega t}]=-\frac{1}{c}{\partial }_t\overrightarrow{B}(t,\overrightarrow{x})=\frac{1}{2}[-\overrightarrow{B}\left(\overrightarrow{x}\right)\frac{1}{c}{\partial }_t{e}^{-i\omega t}-{\overrightarrow{B}}^{\star }\left(\overrightarrow{x}\right)\frac{1}{c}{\partial }_t{e}^{i\omega t}]=$

$=\frac{1}{2}[\overrightarrow{B}\left(\overrightarrow{x}\right)i\frac{\omega }{c}{e}^{-i\omega t}-{\overrightarrow{B}}^{\star }\left(\overrightarrow{x}\right)i\frac{\omega }{c}{e}^{i\omega t}]=\frac{1}{2}[i\frac{\omega }{c}{\overrightarrow{B}}_0{e}^{i\overrightarrow{k}\cdot \overrightarrow{x}}{e}^{-i\omega t}-i\frac{\omega }{c}{\overrightarrow{B}}_0^{\star }{e}^{-i\overrightarrow{k}\cdot \overrightarrow{x}}{e}^{i\omega t}]\Rightarrow$

$\frac{\omega }{c}{\overrightarrow{B}}_0=\overrightarrow{k}\times {\overrightarrow{E}}_0\Rightarrow \overrightarrow{B}\left(\overrightarrow{x}\right)=\frac{c}{\omega }\overrightarrow{k}\times \overrightarrow{E}\left(\overrightarrow{x}\right)$ .

Bei verschwindenden Ladungen und Strömen (d.h. im Vakuum, aber bei Anwesenheit polarisierbarer, isotroper Medien) gilt zudem:

$\overrightarrow{D}=\epsilon \overrightarrow{E},\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{D}=0\Rightarrow \overrightarrow{k}\cdot \overrightarrow{E}\left(\overrightarrow{x}\right)=0$.

Somit erhalten wir

$\overrightarrow{E}\times {\overrightarrow{B}}^{\star }=\frac{c}{\omega }\overrightarrow{E}\times (\overrightarrow{k}\times {\overrightarrow{E}}^{\star })=\frac{c}{\omega }({\left|\overrightarrow{E}\right|}^2\overrightarrow{k}-\underset{=0}{\underbrace{(\overrightarrow{k}\cdot \overrightarrow{E})}}{\overrightarrow{E}}^{\star })=\frac{c}{\omega }\overrightarrow{k}{\left|\overrightarrow{E}\right|}^2$

und daraus wiederum

$⟨\overrightarrow{S}(t,\overrightarrow{x})⟩=\frac{c}{8\pi }\mathrm{\Re }(\overrightarrow{E}\left(\overrightarrow{x}\right)\times {\overrightarrow{H}}^{\star }\left(\overrightarrow{x}\right))=\frac{c^2}{8\pi \mu \omega }\overrightarrow{k}{\left|\overrightarrow{E}\right|}^2\underset{\left|\overrightarrow{k}\right|=\frac{\omega }{c_n},n=\sqrt{\epsilon \mu },c_n=c/n}{=}\frac{c_n}{8\pi }\epsilon {\left|\overrightarrow{E}\right|}^2\frac{\overrightarrow{k}}{\left|\overrightarrow{k}\right|}$.

Jetzt berechnen wir noch die mittlere Energie,

$⟨u(t,\overrightarrow{x})⟩=\underset{T\to \mathrm{\infty }}{\lim }\frac{1}{2T}\underset{-T}{\stackrel{T}{\int }}dtu(t,\overrightarrow{x}),u=\frac{1}{8\pi }(\overrightarrow{E}\cdot \overrightarrow{D}+\overrightarrow{H}\cdot \overrightarrow{B})$.

Hierzu benötigen wir

$\overrightarrow{E}(t,\overrightarrow{x})\cdot \overrightarrow{D}(t,\overrightarrow{x})=\frac{1}{4}[\overrightarrow{E}\left(\overrightarrow{x}\right)\cdot \overrightarrow{D}\left(\overrightarrow{x}\right)e^{-i2\omega t}+{\overrightarrow{E}}^{\star }\left(\overrightarrow{x}\right)\cdot {\overrightarrow{D}}^{\star }\left(\overrightarrow{x}\right)e^{i2\omega t}+{\overrightarrow{E}}^{\star }\left(\overrightarrow{x}\right)\cdot \overrightarrow{D}\left(\overrightarrow{x}\right)+\overrightarrow{E}\left(\overrightarrow{x}\right)\cdot {\overrightarrow{D}}^{\star }\left(\overrightarrow{x}\right)]$,

also

$⟨\overrightarrow{E}(t,\overrightarrow{x})\cdot \overrightarrow{D}(t,\overrightarrow{x})⟩=\frac{1}{4}[{\overrightarrow{E}}^{\star }\left(\overrightarrow{x}\right)\cdot \overrightarrow{D}\left(\overrightarrow{x}\right)+\overrightarrow{E}\left(\overrightarrow{x}\right)\cdot {\overrightarrow{D}}^{\star }\left(\overrightarrow{x}\right)]=\frac{1}{2}\epsilon {\left|\overrightarrow{E}\left(\overrightarrow{x}\right)\right|}^2$

und dazu analog:

$⟨\overrightarrow{H}(t,\overrightarrow{x})\cdot \overrightarrow{B}(t,\overrightarrow{x})⟩=\frac{1}{4}[{\overrightarrow{H}}^{\star }\left(\overrightarrow{x}\right)\cdot \overrightarrow{B}\left(\overrightarrow{x}\right)+\overrightarrow{H}\left(\overrightarrow{x}\right)\cdot {\overrightarrow{B}}^{\star }\left(\overrightarrow{x}\right)]=\frac{1}{2\mu }{\left|\overrightarrow{B}\left(\overrightarrow{x}\right)\right|}^2$.

Wegen

$\overrightarrow{B}=\frac{c}{\omega }\overrightarrow{k}\times \overrightarrow{E},\overrightarrow{k}\cdot \overrightarrow{E}=0\Rightarrow$ ${\left|\overrightarrow{B}\right|}^2=\overrightarrow{B}\cdot {\overrightarrow{B}}^{\star }={\left(\frac{c}{\omega }\right)}^2(\overrightarrow{k}\times \overrightarrow{E})\cdot (\overrightarrow{k}\times {\overrightarrow{E}}^{\star })={\left(\frac{c}{\omega }\right)}^2\overrightarrow{k}\cdot (\overrightarrow{E}\times (\overrightarrow{k}\times {\overrightarrow{E}}^{\star }))={\left(\frac{c}{\omega }\right)}^2\overrightarrow{k}\cdot ({\left|\overrightarrow{E}\right|}^2\overrightarrow{k}-\underset{=0}{\underbrace{(\overrightarrow{k}\cdot \overrightarrow{E})}}{\overrightarrow{E}}^{\star })={\left(\frac{c}{\omega }\right)}^2{\overrightarrow{k}}^2{\left|\overrightarrow{E}\right|}^2$

erhalten wir somit folgenden Ausdruck für die mittlere Energie:

$⟨u(t,\overrightarrow{x})⟩=\frac{1}{2}\epsilon {\left|\overrightarrow{E}\left(\overrightarrow{x}\right)\right|}^2+\frac{1}{2\mu }{\left|\overrightarrow{B}\left(\overrightarrow{x}\right)\right|}^2=\frac{1}{16\pi }(\epsilon {\left|\overrightarrow{E}\left(\overrightarrow{x}\right)\right|}^2+\frac{1}{\mu }{\left(\frac{c}{\omega }\right)}^2{\overrightarrow{k}}^2{\left|\overrightarrow{E}\right|}^2)\underset{\left|\overrightarrow{k}\right|=\frac{\omega }{c_n},n=\sqrt{\epsilon \mu },c_n=c/n}{=\frac{\epsilon }{8\pi }{\left|\overrightarrow{E}\left(\overrightarrow{x}\right)\right|}^2}$.

Der Zusammenhang zwischen dem Mittelwert des Poyntingvektors und der elektromagnetischen Energie lautet daher

$⟨\overrightarrow{S}(t,\overrightarrow{x})⟩=\frac{c_n}{8\pi }\epsilon {\left|\overrightarrow{E}\right|}^2\frac{\overrightarrow{k}}{\left|\overrightarrow{k}\right|}=c_n⟨u(t,\overrightarrow{x})⟩\frac{\overrightarrow{k}}{\left|\overrightarrow{k}\right|}$

mit der Lichtgeschwindigkeit $c_n=c/n$ in einem Medium mit dem Brechungsindex $n=\sqrt{\epsilon \mu }$.

Polarisation

Im Vakuum, d.h. bei Abwesenheit von Quellen wie Ladungen und Ströme, sind bereits ebene Wellen Lösungen der Maxwell-Gleichungen:

$\overrightarrow{E}(t,\overrightarrow{x})=\mathrm{\Re }\left(\overrightarrow{E}\left(\overrightarrow{x}\right)e^{-i\omega t}\right)=\frac{1}{2}[\overrightarrow{E}\left(\overrightarrow{x}\right)e^{-i\omega t}+{\overrightarrow{E}}^{\star }\left(\overrightarrow{x}\right)e^{i\omega t}]$, $\overrightarrow{E}\left(\overrightarrow{x}\right)={\overrightarrow{E}}_0{e}^{i\overrightarrow{k}\cdot \overrightarrow{x}}$

Zudem gilt für die Wellen die Transversalität im Vakuum:

$\overrightarrow{k}\cdot {\overrightarrow{E}}_0=0\Rightarrow {\overrightarrow{E}}_0=E_1{\hat{e}}_1+E_2{\hat{e}}_2,\overrightarrow{k}\cdot {\hat{e}}_i=0,{\hat{e}}_i^2=1(i=1,2),{\hat{e}}_1\cdot {\hat{e}}_2=0,\Rightarrow {\hat{e}}_1\times {\hat{e}}_2=\frac{\overrightarrow{k}}{k}$

${\overrightarrow{E}}_0$ ist hierbei im Allg. komplex: $E_j=\left|E_j\right|e^{i{\phi }_j},(i=1,2)$.

Mit letzterer Gleichung und mit Hilfe der Euler-Formel $e^{ix}=\cos x+i\sin x$ ($\Rightarrow e^{\pm i\frac{\pi }{2}}=\pm i$) können wir die ebene Welle umschreiben:

$\overrightarrow{E}(t,\overrightarrow{x})=\mathrm{\Re }\left(\overrightarrow{E}\left(\overrightarrow{x}\right)e^{-i\omega t}\right)=\mathrm{\Re }\left({\overrightarrow{E}}_0{e}^{-i\omega t+i\overrightarrow{k}\cdot \overrightarrow{x}}\right)$ $=\mathrm{\Re }({\hat{e}}_1\left|E_1\right|e^{-i\omega t+i\overrightarrow{k}\cdot \overrightarrow{x}+i{\phi }_1}+{\hat{e}}_2\left|E_2\right|e^{-i\omega t+i\overrightarrow{k}\cdot \overrightarrow{x}+i{\phi }_2})$ $={\hat{e}}_1\left|E_1\right|\cos (\omega t-\overrightarrow{k}\cdot \overrightarrow{x}-{\phi }_1)+{\hat{e}}_2\left|E_2\right|\cos (\omega t-\overrightarrow{k}\cdot \overrightarrow{x}-{\phi }_2)$.

Je nach der gegenseitigen Beziehung der Phasen ${\phi }_1$ bzw. ${\phi }_2$ werden gerne folgende Formen der sog. Polarisation unterschieden:

• Lineare Polarisation: $E_1$ und $E_2$

sind in Phase, d.h. ${\phi }_1={\phi }_2$, z.B. ${\phi }_1={\phi }_2=0$.

• Zirkulare Polarisation: $\left|E_1\right|=\left|E_2\right|=E_0$,

$|{\phi }_1-{\phi }_2|=\frac{\pi }{2}$, z.B. ${\phi }_1=0,{\phi }_2=\pm \frac{\pi }{2}$:

$\overrightarrow{E}(t,\overrightarrow{x})={\hat{e}}_1{E}_0\cos (\omega t-\overrightarrow{k}\cdot \overrightarrow{x})\pm {\hat{e}}_2{E}_0\sin (\omega t-\overrightarrow{k}\cdot \overrightarrow{x})$ erfüllt die Kreisgleichung ${({\hat{e}}_1\cdot \overrightarrow{E}(t,\overrightarrow{x}))}^2+{({\hat{e}}_2\cdot \overrightarrow{E}(t,\overrightarrow{x}))}^2=E_0^2=const.$;

${\overrightarrow{E}}_0=E_1{\hat{e}}_1+E_2{\hat{e}}_2=E_0({\hat{e}}_1\pm i{\hat{e}}_2)$, woraus die Definition der zirkularen Einheitsvektoren folgt:

${\widehat{ϵ}}_{\pm }=\frac{1}{\sqrt{2}}({\hat{e}}_1\pm i{\hat{e}}_2),{ϵ}_3=\frac{\overrightarrow{k}}{k}\Rightarrow {\widehat{ϵ}}_{\pm }\cdot {ϵ}_3=0,{\widehat{ϵ}}_{\pm }\cdot {\widehat{ϵ}}_{\pm }^{\star }={\widehat{ϵ}}_{\pm }\cdot {\widehat{ϵ}}_{\mp }=1,{\widehat{ϵ}}_{\pm }\cdot {\widehat{ϵ}}_{\mp }^{\star }={\widehat{ϵ}}_{\pm }\cdot {\widehat{ϵ}}_{\pm }=0$.

Eine Zerlegung von ${\overrightarrow{E}}_0$ ist daher auch bezüglich der ${\widehat{ϵ}}_{\pm }$ möglich: ${\overrightarrow{E}}_0=E_{+}{\widehat{ϵ}}_{+}+E_{-}{\widehat{ϵ}}_{-}$.

Die sog. Helizität ist der Drehsinn der Welle parallel zum Wellenvektor $\overrightarrow{k}:{\widehat{ϵ}}_{\pm }$ gehören zu den Helizitäten $\pm 1$ .

Impuls des Strahlungsfeldes

Ausgehend von der Lorentz-Kraft können wir eine Größe gewinnen, die etwas über den Impuls eines Strahlungsfeldes ausdrückt:

$\frac{1}{c}{\partial }_t{\overrightarrow{P}}_M=\overrightarrow{F}=\int d^{3}x(\underset{=\frac{1}{4\pi }\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{D}}{\underbrace{\rho (t,\overrightarrow{x})}}\overrightarrow{E}(t,\overrightarrow{x})+\frac{1}{c}\underset{=\frac{c}{4\pi }\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{H}-\frac{1}{4\pi }{\partial }_t\overrightarrow{D}}{\underbrace{\overrightarrow{J}(t,\overrightarrow{x})}}\times \overrightarrow{B}(t,\overrightarrow{x}))$ $=\frac{1}{4\pi }\int d^{3}x(\overrightarrow{E}(\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{D})+(\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{H})\times \overrightarrow{B}-\frac{1}{c}{\partial }_t\overrightarrow{D}\times \overrightarrow{B})$.

Hierin verwenden wir

$\frac{1}{c}{\partial }_t(\overrightarrow{D}\times \overrightarrow{B})=\frac{1}{c}{\partial }_t\overrightarrow{D}\times \overrightarrow{B}+\overrightarrow{D}\times \underset{=-\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E}}{\underbrace{\frac{1}{c}{\partial }_t\overrightarrow{B}}}=\frac{1}{c}{\partial }_t\overrightarrow{D}\times \overrightarrow{B}-\overrightarrow{D}\times (\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E})$,

sodass

$\frac{1}{c}{\partial }_t{\overrightarrow{P}}_M+\frac{1}{c}{\partial }_t\underset{={\overrightarrow{P}}_F}{\underbrace{\int d^{3}x\frac{1}{4\pi }\overrightarrow{D}\times \overrightarrow{B}}}$ $=\frac{1}{4\pi }\int d^{3}x(\overrightarrow{E}(\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{D})+\overrightarrow{H}\underset{=0}{\underbrace{(\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{B})}}-\overrightarrow{B}\times (\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{H})-\overrightarrow{D}\times (\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E}))$, $\overrightarrow{D}\times (\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E})=\sum _{k,l,m,n}{\epsilon }_{ikl}{\epsilon }_{lmn}{D}_k{\partial }_m{E}_n=\sum _{k,l,m,n}({\delta }_{im}{\delta }_{kn}-{\delta }_{in}{\delta }_{km})D_k{\partial }_m{E}_n=\overrightarrow{D}\cdot {\partial }_i\overrightarrow{E}-(\overrightarrow{D}\cdot \overrightarrow{\mathrm{\nabla }})E_i$, $(\overrightarrow{D}\cdot \overrightarrow{\mathrm{\nabla }})E_i=\overrightarrow{\mathrm{\nabla }}\cdot \left(\overrightarrow{D}{E}_i\right)-E_i(\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{D})$.

Außerdem nehmen wir einen linearen Zusammenhang zwischen $\overrightarrow{E}$ und $\overrightarrow{D}$ (wie z.B. isotrope Dielektrika: $\underset{\_}{\epsilon }=\epsilon \underset{\_}{1}$) an:

$\overrightarrow{D}=\underset{\_}{\epsilon }\overrightarrow{E}\Rightarrow$ ${\partial }_i(\overrightarrow{D}\cdot \overrightarrow{E})={\partial }_i({\overrightarrow{D}}^T\cdot \overrightarrow{E})={\partial }_i\left({\overrightarrow{E}}^T{\underset{\_}{\epsilon }}^T\overrightarrow{E}\right)$ $={\left({\partial }_i\overrightarrow{E}\right)}^T{\underset{\_}{\epsilon }}^T\overrightarrow{E}+{\overrightarrow{E}}^T{\underset{\_}{\epsilon }}^T{\partial }_i\overrightarrow{E}={\left({\partial }_i\overrightarrow{E}\right)}^T{\left({\overrightarrow{E}}^T\underset{\_}{\epsilon }\right)}^T+\overrightarrow{D}\cdot {\partial }_i\overrightarrow{E}$;

Skalarprodukt: $\overrightarrow{a}\cdot \overrightarrow{b}={\overrightarrow{a}}^T\overrightarrow{b}=\overrightarrow{b}\cdot \overrightarrow{a}={\overrightarrow{b}}^T\overrightarrow{a}\Rightarrow$

Wenn $\underset{\_}{\epsilon }$ symmetrisch ist, d.h.

${\underset{\_}{\epsilon }}^T=\underset{\_}{\epsilon }$

${\left({\partial }_i\overrightarrow{E}\right)}^T{\left({\overrightarrow{E}}^T\underset{\_}{\epsilon }\right)}^T={\overrightarrow{E}}^T\underset{\_}{\epsilon }{\partial }_i\overrightarrow{E}={\left({\underset{\_}{\epsilon }}^T\overrightarrow{E}\right)}^T{\partial }_i\overrightarrow{E}={\left(\underset{\_}{\epsilon }\overrightarrow{E}\right)}^T{\partial }_i\overrightarrow{E}=\overrightarrow{D}\cdot {\partial }_i\overrightarrow{E}\Rightarrow$

${\partial }_i(\overrightarrow{D}\cdot \overrightarrow{E})=2\overrightarrow{D}\cdot {\partial }_i\overrightarrow{E}$

Somit erhalten wir :

${[\overrightarrow{D}\times (\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E})]}_i=\sum _k{\partial }_k(\frac{1}{2}{\delta }_{ki}\overrightarrow{D}\cdot \overrightarrow{E}-E_i{D}_k)+\overrightarrow{E}(\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{D})$.

Ganz analog hierzu nehmen wir einen linearen Zusammenhang zwischen $\overrightarrow{B}$ und $\overrightarrow{H}$ (wie z.B. isotrope Magnetika: $\underset{\_}{\mu }=\mu \underset{\_}{1}$) an:

$\overrightarrow{B}=\underset{\_}{\mu }\overrightarrow{H}$.

$\underset{\_}{\mu }$ sei symmetrisch: ${\underset{\_}{\mu }}^T=\underset{\_}{\mu }$:

${[\overrightarrow{B}\times (\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{H})]}_i=\sum _k{\partial }_k(\frac{1}{2}{\delta }_{ki}\overrightarrow{B}\cdot \overrightarrow{H}-B_i{H}_k)+\overrightarrow{H}(\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{B})$.

Somit ergibt sich folgender Zusammenhang,

$\frac{1}{c}{\partial }_t{\left[{\overrightarrow{P}}_M\right]}_i+\frac{1}{c}{\partial }_t{\left[{\overrightarrow{P}}_F\right]}_i=\sum _k\int d^{3}x{\partial }_k{T}^{ik}$,

zwischen dem mechanischen Impuls ${\overrightarrow{P}}_M$ und dem Impuls des elektromagnetischen Feldes ${\overrightarrow{P}}_F=\int d^{3}x\frac{1}{4\pi }\overrightarrow{D}\times \overrightarrow{B}$.

Hierin kommt der Maxwell'sche Spannungstensor $T^{ik}={\delta }_{ik}w-\frac{1}{4\pi }(E_i{D}_k+H_i{B}_k)$ mit der Energiedichte $w=\frac{1}{8\pi }(\overrightarrow{D}\cdot \overrightarrow{E}+\overrightarrow{B}\cdot \overrightarrow{H})$ vor.

Elektromagnetische Potentiale

Die elektromagnetischen Feldvektoren lassen sich mittels sog. Potentiale darstellen, wie z.B. die magnetische Flussdichte als Rotation eines Vektorpotentials $\overrightarrow{A}$:

$\overrightarrow{B}=\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{A}$,

weil dies automatisch die Maxwell-Gleichung $\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{B}=0$ erfüllt. Für das elektrische Feld schließen wir hieraus auf die Existenz eines Skalarpotentials $A_0$:

$\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E}=-\frac{1}{c}{\partial }_t\overrightarrow{B}\Rightarrow 0=\overrightarrow{\mathrm{\nabla }}\times (\overrightarrow{E}+\frac{1}{c}{\partial }_t\overrightarrow{A})=-\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{\mathrm{\nabla }}{A}_0$,

weil:

${[\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{\mathrm{\nabla }}{A}_0]}_i={\epsilon }_{ikl}{\partial }_k{\partial }_l{A}_0=\frac{1}{2}({\epsilon }_{ikl}+{\epsilon }_{ikl}){\partial }_k{\partial }_l{A}_0=\frac{1}{2}({\epsilon }_{ikl}{\partial }_k{\partial }_l{A}_0+{\epsilon }_{ilk}{\partial }_l{\partial }_k{A}_0)$ $=\frac{1}{2}({\epsilon }_{ikl}+{\epsilon }_{ilk}){\partial }_k{\partial }_l{A}_0=\frac{1}{2}({\epsilon }_{ikl}-{\epsilon }_{ikl}){\partial }_k{\partial }_l{A}_0=0$.

Das elektrische Feld lässt sich also wie folgt mit Hilfe des Skalarpotentials $A_0$ und des Vektorpotentials $\overrightarrow{A}$ angeben:

$\overrightarrow{E}=-\overrightarrow{\mathrm{\nabla }}{A}_0-\frac{1}{c}{\partial }_t\overrightarrow{A}$.

Wellengleichungen für elektromagnetische Potentiale

Genau so wie die elektromagnetischen Felder Wellengleichungen erfüllen, ist dies auch für die elektromagnetischen Potentiale der Fall, aus denen sie berechnet werden können.

(1) Wellengleichung für das Vektorpotential:

$\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{B}=\frac{4\pi }{c}\overrightarrow{J}+\frac{1}{c}{\partial }_t\overrightarrow{E},\overrightarrow{B}=\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{A},\overrightarrow{E}=-\overrightarrow{\mathrm{\nabla }}{A}_0-\frac{1}{c}{\partial }_t\overrightarrow{A}\Rightarrow$ $\overrightarrow{\mathrm{\nabla }}(\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A})-{\overrightarrow{\mathrm{\nabla }}}^2\overrightarrow{A}=\overrightarrow{\mathrm{\nabla }}\times (\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{A})=\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{B}$ $=\frac{4\pi }{c}\overrightarrow{J}+\frac{1}{c}{\partial }_t\overrightarrow{E}=\frac{4\pi }{c}\overrightarrow{J}-\overrightarrow{\mathrm{\nabla }}\frac{1}{c}{\partial }_t{A}_0-\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}\overrightarrow{A}$.

D.h. wir erhalten die Wellengleichung

$(\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}-{\overrightarrow{\mathrm{\nabla }}}^2)\overrightarrow{A}=\frac{4\pi }{c}\overrightarrow{J}-\overrightarrow{\mathrm{\nabla }}(\frac{1}{c}{\partial }_t{A}_0+\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A})$.

Anm.: Die Schlussfolgerungen aus folgenden analogen Operationen sind hingegen eher trivial:

$\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E}=-\frac{1}{c}{\partial }_t\overrightarrow{B},\overrightarrow{B}=\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{A},\overrightarrow{E}=-\overrightarrow{\mathrm{\nabla }}{A}_0-\frac{1}{c}{\partial }_t\overrightarrow{A}\Rightarrow$ $-\overrightarrow{\mathrm{\nabla }}\times \frac{1}{c}{\partial }_t\overrightarrow{A}=-\frac{1}{c}{\partial }_t\overrightarrow{B}=\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{E}=-\underset{=0}{\underbrace{\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{\mathrm{\nabla }}{A}_0}}-\overrightarrow{\mathrm{\nabla }}\times \frac{1}{c}{\partial }_t\overrightarrow{A}=-\overrightarrow{\mathrm{\nabla }}\times \frac{1}{c}{\partial }_t\overrightarrow{A}$.

Bei Anwesenheit isotroper Magnetika bzw. Dielektrika müssen wir statt $\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{B}=\frac{4\pi }{c}\overrightarrow{J}+\frac{1}{c}{\partial }_t\overrightarrow{E}$ die Gleichung $\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{B}=\frac{4\pi \epsilon \mu }{c}\tilde{\overrightarrow{J}}+\frac{\epsilon \mu }{c}{\partial }_t\overrightarrow{E}$ mit $\tilde{\overrightarrow{J}}=\frac{1}{\epsilon }\overrightarrow{J}$ verwenden und erhalten daher:

$(\frac{1}{c_n^2}\underset{t}{\overset{2}{\partial }}-{\overrightarrow{\mathrm{\nabla }}}^2)\overrightarrow{A}=\frac{4\pi n}{c_n\epsilon }\overrightarrow{J}-\overrightarrow{\mathrm{\nabla }}(\frac{n}{c_n}{\partial }_t{A}_0+\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A})$

mit der Lichtgeschwindigkeit im Medium $c_n=\frac{c}{n}$ und dem Brechungsindex $n=\sqrt{\epsilon \mu }$.

(2) Wellengleichung für das skalare Potential:

$\overrightarrow{E}=-\overrightarrow{\mathrm{\nabla }}{A}_0-\frac{1}{c}{\partial }_t\overrightarrow{A}\Rightarrow$ $-{\overrightarrow{\mathrm{\nabla }}}^2{A}_0-\frac{1}{c}{\partial }_t\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A}=\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{E}=4\pi \rho$.

Bei Anwesenheit isotroper Magnetika bzw. Dielektrika müssen wir statt $\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{E}=4\pi \rho$ die Gleichung $\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{E}=\frac{4\pi }{\epsilon }\rho$ verwenden:

$-{\overrightarrow{\mathrm{\nabla }}}^2{A}_0-\frac{1}{c}{\partial }_t\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A}=\frac{4\pi }{\epsilon }\rho$.

Eichfreiheit

Bei Abwesenheit von Ladungen und Strömen gilt ja

$(\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}-{\overrightarrow{\mathrm{\nabla }}}^2)\overrightarrow{E}=0$,

wobei die allgemeine Lösung

$\overrightarrow{E}(t,\overrightarrow{x})=\int \frac{d^{3}k}{{\left(2\pi \right)}^3}({\overrightarrow{E}}_0{e}^{-i\omega t+i\overrightarrow{k}\cdot \overrightarrow{x}}+{\overrightarrow{E}}_0^{\star }{e}^{+i\omega t-i\overrightarrow{k}\cdot \overrightarrow{x}})$

mit $\omega =kc$ lautet (die Lösung $\overrightarrow{E}$ ist reell, aber die Amplitude ${\overrightarrow{E}}_0$ wird im Allg. komplex sein).

Analoges gilt für $\overrightarrow{B}$:

$(\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}-{\overrightarrow{\mathrm{\nabla }}}^2)\overrightarrow{B}=0$.

Hieraus folgt auch der Spezialfall ebener Wellen:

$\overrightarrow{E}(t,\overrightarrow{x})={\overrightarrow{E}}_0{e}^{-i\omega t+i\overrightarrow{k}\cdot \overrightarrow{x}}$, $\overrightarrow{B}(t,\overrightarrow{x})={\overrightarrow{B}}_0{e}^{-i\omega t+i\overrightarrow{k}\cdot \overrightarrow{x}}$.

Mit $\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{B}=0$ und $\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{E}=0$ (d.h. es sind keine Ladung vorhanden) folgt hieraus:

$\overrightarrow{k}\cdot {\overrightarrow{B}}_0=0$

und

$\overrightarrow{k}\cdot {\overrightarrow{E}}_0=0$,

d.h. die Wellen sind transversal. $A_0$ und $\overrightarrow{A}$ besitzen zusammen vier Komponenten. Wegen der beiden Transversalitätsbedingung sind also zwei Komponenten zuviel, um die Welle beschreiben zu wollen: Es sind vier Komponenten vorhanden, es werden aber wegen der Transversalität der elektromagnetischen Wellen nur zwei benötigt. Die überschüssigen Freiheitsgrade werden daher »weggeeicht« .

Eine Transformation

$A_0\to A_0-\frac{1}{c}\frac{\partial }{\partial t}\mathrm{\Lambda },\overrightarrow{A}\to \overrightarrow{A}+\overrightarrow{\mathrm{\nabla }}\mathrm{\Lambda }$

ändert $\overrightarrow{E}$ und $\overrightarrow{B}$ nicht (und heißt »Eichtransformation«):

$\overrightarrow{E}=-\overrightarrow{\mathrm{\nabla }}{A}_0-\frac{1}{c}{\partial }_t\overrightarrow{A}\to -\overrightarrow{\mathrm{\nabla }}{A}_0-\frac{1}{c}{\partial }_t\overrightarrow{A}+\overrightarrow{\mathrm{\nabla }}\frac{1}{c}\frac{\partial }{\partial t}\mathrm{\Lambda }-\frac{1}{c}{\partial }_t\overrightarrow{\mathrm{\nabla }}\mathrm{\Lambda }=\overrightarrow{E}$, $\overrightarrow{B}=\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{A}\to \overrightarrow{\mathrm{\nabla }}\times \overrightarrow{A}+\underset{=0}{\underbrace{\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{\mathrm{\nabla }}\mathrm{\Lambda }}}=\overrightarrow{B}$,

Letzteres weil $rotgrad\equiv 0$. Folgender Term aus den Wellengleichungen verhält sich unter dieser Transformation hingegen wie:

$\frac{1}{c}{\partial }_t{A}_0+\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A}\to \frac{1}{c}{\partial }_t{A}_0+\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A}-(\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}-{\overrightarrow{\mathrm{\nabla }}}^2)\mathrm{\Lambda }$.

Die bekanntesten Beispiele für Eichungen sind:

• Lorenz-Eichung: Wähle $\mathrm{\Lambda }$ so, dass unter

der Transformation $\frac{1}{c}{\partial }_t{A}_0+\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A}\to 0$ gilt, d.h.

$(\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}-{\overrightarrow{\mathrm{\nabla }}}^2)\mathrm{\Lambda }=\frac{1}{c}{\partial }_t{A}_0+\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A}$.

• Coulomb-Eichung: Wähle $\mathrm{\Lambda }$ so, dass unter der Transformation

$\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A}\to 0$ gilt, d.h.

$(\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}-{\overrightarrow{\mathrm{\nabla }}}^2)\mathrm{\Lambda }=\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A}$.

In den folgenden beiden Kapiteln möchten wir gerne diese beiden Sorten von Eichungen auf die mikroskopischen Wellengleichungen anwenden.

Mikroskopische Wellengleichungen in Coulomb-Eichung

Die mikroskopische Wellengleichung für Skalar- und Vektorpotential möchten wir hier in der Coulomb-Eichung, d.h.

$\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A}=0$,

betrachten. Für das Vektorpotential gilt dann:

$(\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}-{\overrightarrow{\mathrm{\nabla }}}^2)\overrightarrow{A}=\frac{4\pi }{c}\overrightarrow{J}-\overrightarrow{\mathrm{\nabla }}(\frac{1}{c}{\partial }_t{A}_0+\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A})\Rightarrow$ $(\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}-{\overrightarrow{\mathrm{\nabla }}}^2)\overrightarrow{A}=\frac{4\pi }{c}\overrightarrow{J}-\overrightarrow{\mathrm{\nabla }}\frac{1}{c}{\partial }_t{A}_0=:\frac{4\pi }{c}{\overrightarrow{J}}_{\perp }$.

Für das Skalarpotential erhalten wir die Gleichung:

$-{\overrightarrow{\mathrm{\nabla }}}^2{A}_0-\frac{1}{c}{\partial }_t\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A}=4\pi \rho \Rightarrow {\overrightarrow{\mathrm{\nabla }}}^2{A}_0=-4\pi \rho$,

woraus wir wie folgt das Skalarpotential bestimmen können:

${\overrightarrow{\mathrm{\nabla }}}^2{A}_0(t,\overrightarrow{x})=-4\pi \rho (t,\overrightarrow{x})=\int d^3{x}^{\prime }\rho (t,{\overrightarrow{x}}^{\prime })\underset{{\overrightarrow{\mathrm{\nabla }}}^2\frac{1}{|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|}}{\underbrace{(-4\pi {\delta }^3(\overrightarrow{x}-{\overrightarrow{x}}^{\prime }))}}$ $=\int d^3{x}^{\prime }\rho (t,{\overrightarrow{x}}^{\prime }){\overrightarrow{\mathrm{\nabla }}}^2\frac{1}{|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|}={\overrightarrow{\mathrm{\nabla }}}^2\int d^3{x}^{\prime }\frac{\rho (t,{\overrightarrow{x}}^{\prime })}{|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|}$,

d.h.

$A_0(t,\overrightarrow{x})=\int d^3{x}^{\prime }\frac{\rho (t,{\overrightarrow{x}}^{\prime })}{|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|}$.

Hieraus und aus der Kontinuitätsgleichung ${\partial }_t\rho +\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{J}=0$ folgt zum einen

${\partial }_t{A}_0(t,\overrightarrow{x})=\int d^3{x}^{\prime }\frac{{\partial }_t\rho (t,{\overrightarrow{x}}^{\prime })}{|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|}=-\int d^3{x}^{\prime }\frac{{\overrightarrow{\mathrm{\nabla }}}^{\prime }\cdot \overrightarrow{J}(t,{\overrightarrow{x}}^{\prime })}{|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|}$

und die Transversalität von ${\overrightarrow{J}}_{\perp }$:

$\frac{4\pi }{c}{\overrightarrow{J}}_{\perp }=\frac{4\pi }{c}\overrightarrow{J}-\overrightarrow{\mathrm{\nabla }}\frac{1}{c}{\partial }_t{A}_0\Rightarrow$ $\frac{4\pi }{c}\overrightarrow{\mathrm{\nabla }}\times {\overrightarrow{J}}_{\perp }=\frac{4\pi }{c}\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{J}$, $\frac{4\pi }{c}\overrightarrow{\mathrm{\nabla }}\cdot {\overrightarrow{J}}_{\perp }=\frac{4\pi }{c}\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{J}-\frac{1}{c}{\partial }_t\underset{-4\pi \rho }{\underbrace{{\overrightarrow{\mathrm{\nabla }}}^2{A}_0}}=\frac{4\pi }{c}(\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{J}+{\partial }_t\rho )=0$.

Die Stromdichte können wir daher in einen transversalen und einen longitudinalen Anteil zerlegen:

$\overrightarrow{J}={\overrightarrow{J}}_{\perp }+{\overrightarrow{J}}_{\parallel },\overrightarrow{\mathrm{\nabla }}\cdot {\overrightarrow{J}}_{\perp }=0,\overrightarrow{\mathrm{\nabla }}\times {\overrightarrow{J}}_{\parallel }:=0\Rightarrow$ $\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{J}=\overrightarrow{\mathrm{\nabla }}\cdot {\overrightarrow{J}}_{\parallel },\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{J}=\overrightarrow{\mathrm{\nabla }}\times {\overrightarrow{J}}_{\perp }$

mit dem longitudinalen Anteil

${\overrightarrow{J}}_{\parallel }=\frac{1}{4\pi }\overrightarrow{\mathrm{\nabla }}{\partial }_t{A}_0(t,\overrightarrow{x})$.

Mikroskopische Wellengleichungen in Lorenz-Eichung

Die mikroskopische Wellengleichung für Skalar- und Vektorpotential möchten wir hier in der Lorenz-Eichung,

$\frac{1}{c}{\partial }_t{A}_0+\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A}=0$,

diskutieren. Für das Vektorpotential erhalten wir

$(\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}-{\overrightarrow{\mathrm{\nabla }}}^2)\overrightarrow{A}=\frac{4\pi }{c}\overrightarrow{J}-\overrightarrow{\mathrm{\nabla }}(\frac{1}{c}{\partial }_t{A}_0-\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A})\Rightarrow$ $(\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}-{\overrightarrow{\mathrm{\nabla }}}^2)\overrightarrow{A}=\frac{4\pi }{c}\overrightarrow{J}$,

während für das Skalarpotential Folgendes gilt:

$4\pi \rho =-{\overrightarrow{\mathrm{\nabla }}}^2{A}_0-\frac{1}{c}{\partial }_t\underset{=-\frac{1}{c}{\partial }_t{A}_0}{\underbrace{\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A}}}=(\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}-{\overrightarrow{\mathrm{\nabla }}}^2)A_0$.

Zusammengefasst zu einer Gleichung mittels (der. sog. »Vierervektoren«) $A^{\mu }=(A_0,\overrightarrow{A})$, $J^{\mu }=(J_0,\overrightarrow{J})=(c\rho ,\overrightarrow{J})$, $x=(x_0,\overrightarrow{x})=(ct,\overrightarrow{x})$, wobei hier und in Zukunft griechische Indizes immer Werte von 0 bis 3 annehmen sollen, $\mu =0,1,2,3$, (während wir vereinbaren, dass römische Indizes nur die Werte von 1 bis 3 durchlaufen) erhalten wir

$◻A^{\mu }=(\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}-{\overrightarrow{\mathrm{\nabla }}}^2)A^{\mu }=\frac{4\pi }{c}{J}^{\mu }$.

Mit der »Vierer-Ableitung«,

${\partial }_{\mu }=({\partial }_0,\overrightarrow{\mathrm{\nabla }})=(\frac{1}{c}{\partial }_t,\overrightarrow{\mathrm{\nabla }})$,

können wir auch die Lorenz-Eichbedingung folgendermaßen schreiben:

$\underset{\mu =0}{\sum ^4}{\partial }_{\mu }{A}^{\mu }=\frac{1}{c}{\partial }_t{A}_0+\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A}=0$.

Lösungen der mikroskopischen Wellengleichungen

Mit Hilfe der zuvor bereits vorgestellten Greensfunktion der Elektrodynamik können wir jetzt sogar eine Lösung der Wellengleichungen in Lorenz-Eichung ($\frac{1}{c}{\partial }_t{A}_0+\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A}=0$), d.h.

$◻A^{\mu }=(\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}-{\overrightarrow{\mathrm{\nabla }}}^2)A^{\mu }=\frac{4\pi }{c}{J}^{\mu }$,

bzw. in Coulomb-Eichung ($\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{A}=0$), d.h.

$◻\overrightarrow{A}=(\frac{1}{c^2}\underset{t}{\overset{2}{\partial }}-{\overrightarrow{\mathrm{\nabla }}}^2)\overrightarrow{A}=\frac{4\pi }{c}\overrightarrow{J}-\overrightarrow{\mathrm{\nabla }}\frac{1}{c}{\partial }_t{A}_0=:\frac{4\pi }{c}{\overrightarrow{J}}_{\perp }$,

${\overrightarrow{\mathrm{\nabla }}}^2{A}_0=-4\pi \rho$ (mit der bereits bekannten Lösung $A_0(t,\overrightarrow{x})=\int d^3{x}^{\prime }\frac{\rho (t,{\overrightarrow{x}}^{\prime })}{|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|}$)

angeben. Für die Greensfunktion gilt ja:

$◻G(x-x^{\prime })=4\pi {\delta }^4(x-x^{\prime })$

.

Am Beispiel der Lorenz-Eichung folgt hieraus:

$◻A^{\mu }\left(x\right)=\frac{4\pi }{c}{J}^{\mu }\left(x\right)=\frac{1}{c}\int d^4{x}^{\prime }{J}^{\mu }\left(x^{\prime }\right)4\pi {\delta }^4(x-x^{\prime })=\frac{1}{c}\int d^4{x}^{\prime }{J}^{\mu }\left(x^{\prime }\right)◻G(x-x^{\prime })\Rightarrow$ $A^{\mu }\left(x\right)=\frac{1}{c}\int d^4{x}^{\prime }G(x-x^{\prime })J^{\mu }\left(x^{\prime }\right)$

.

Im nächsten Kapitel werden wir diese formale Lösung verbunden mit der bereits bekannten Gestalt der Greensfunktion verwenden, um Strahlungssysteme zu betrachten.

Strahlungssysteme

Strahlungsquellen, d.h. Ladungen und Ströme, sollen (als Spezialfall einer Fourier-Zerlegung) eine harmonische Zeitabhängigkeit besitzen:

$\rho (t,\overrightarrow{x})=\rho \left(\overrightarrow{x}\right)e^{-i\omega t}$, $\overrightarrow{J}(t,\overrightarrow{x})=\overrightarrow{J}\left(\overrightarrow{x}\right)e^{-i\omega t}$,

d.h.

$J^{\mu }(t,\overrightarrow{x})=J^{\mu }\left(\overrightarrow{x}\right)e^{-i\omega t}$.

In $A^{\mu }\left(x\right)=\frac{1}{c}\int d^4{x}^{\prime }G(x-x^{\prime })J^{\mu }\left(x^{\prime }\right)$ verwenden wir die retardierte Greensfunktion $G_{+}(x_0,\overrightarrow{x})=\mathrm{\Theta }\left(x_0\right)\frac{1}{r}\delta (x_0-r)$, weil nur bei ihr die Wirkung auf der Ursache folgt (und nicht umgekehrt, wie dies bei der avancierten Greensfunktion der Fall ist):

$A^{\mu }\left(x\right)=\frac{1}{c}\int d^4{x}^{\prime }G(x-x^{\prime })J^{\mu }\left(x^{\prime }\right)=$ $\int d^3{x}^{\prime }\int d{t}^{\prime }\frac{J^{\mu }\left({\overrightarrow{x}}^{\prime }\right)e^{-i\omega t^{\prime }}}{|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|}\underset{=\frac{1}{c}\delta (t-t^{\prime }-|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|/c)}{\underbrace{\delta (c(t-t^{\prime })-|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|)}}\mathrm{\Theta }(t-t^{\prime })$ $=\frac{1}{c}\int d^3{x}^{\prime }\frac{J^{\mu }\left({\overrightarrow{x}}^{\prime }\right)e^{i\frac{\omega }{c}|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|}}{|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|}\underset{=1}{\underbrace{\mathrm{\Theta }\left(\stackrel{\ge 0}{\overbrace{\frac{1}{c}|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|}}\right)}}{e}^{-i\omega t}=A^{\mu }\left(\overrightarrow{x}\right)e^{-i\omega t}$.

Somit gilt bei $\overrightarrow{A}(t,\overrightarrow{x})=\overrightarrow{A}\left(\overrightarrow{x}\right)e^{-i\omega t}$:

$\overrightarrow{A}\left(\overrightarrow{x}\right)=\frac{1}{c}\int d^3{x}^{\prime }\frac{\overrightarrow{J}\left({\overrightarrow{x}}^{\prime }\right)e^{ik|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|}}{|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|}$

mit $k=\frac{\omega }{c}$.

Die lineare Ausdehnung d der Quelle sei viel kleiner als die Wellenlänge $\lambda$ der von einer zeitlich veränderlichen Ladungs- und Stromverteilung erzeugten Strahlung: $d\ll \lambda$ (bzw. wegen $k=\frac{2\pi }{\lambda }:kd\ll 1$). Der Beobachter befinde sich im Abstand r von der Strahlungsquelle. Wir betrachten im Folgenden zwei Fälle:

(1) Nahzone: $d\ll r\ll \lambda$, $k=\frac{2\pi }{\lambda }\Rightarrow kr\ll 1\Rightarrow e^{ik|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|}\approx 1$.

Wir können $\frac{1}{|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|}$ in Kugelfunktionen $Y_{lm}\left(\hat{x}\right)=Y_{lm}(\vartheta ,\phi )$ entwickeln:

$\frac{1}{|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|}=\underset{l=0}{\sum ^{\mathrm{\infty }}}\underset{m=-l}{\sum ^l}\frac{4\pi }{2l+1}\frac{r_{<}^l}{r_{>}^{l+1}}{Y}_{lm}\left(\hat{x}\right)Y_{lm}^{\star }\left({\hat{x}}^{\prime }\right)$,

$r_{<}=r^{\prime }$, $r_{>}=r$ , so dass sich in der Nahzone für das Vektorpotential

$\overrightarrow{A}\left(\overrightarrow{x}\right)\underset{kr\ll 1}{\approx }\frac{1}{c}\underset{l=0}{\sum ^{\mathrm{\infty }}}\underset{m=-l}{\sum ^l}\frac{4\pi }{2l+1}\frac{Y_{lm}\left(\hat{x}\right)}{r^{l+1}}\underset{={\mu }_{lm}}{\underbrace{\int d^3{x}^{\prime }\overrightarrow{J}\left({\overrightarrow{x}}^{\prime }\right)r^{\prime l}{Y}_{lm}^{\star }\left({\hat{x}}^{\prime }\right)}}$

mit dem Multipolmoment der Stromdichte,

${\mu }_{lm}=\int d^3{x}^{\prime }\overrightarrow{J}\left({\overrightarrow{x}}^{\prime }\right)r^{\prime l}{Y}_{lm}^{\star }\left({\hat{x}}^{\prime }\right)$,

ergibt.

(2) Fernzone: $d\ll \lambda \ll r$, $k=\frac{2\pi }{\lambda }\Rightarrow kr\gg 1$.

Wir können somit $|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|$ nach Taylor entwickeln:

$|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|\underset{\overrightarrow{n}=\frac{\overrightarrow{x}}{r}}{=}r\sqrt{1+{\left(\frac{r^{\prime }}{r}\right)}^2-2\frac{{\overrightarrow{x}}^{\prime }\cdot \overrightarrow{n}}{r}}\underset{r\gg r^{\prime }\approx d}{\approx }r(1-\frac{{\overrightarrow{x}}^{\prime }\cdot \overrightarrow{n}}{r})\Rightarrow$ $\frac{1}{|\overrightarrow{x}-{\overrightarrow{x}}^{\prime }|}\approx \frac{1}{r}(1+\underset{\underset{r\to \mathrm{\infty }}{⟶}0}{\underbrace{\frac{{\overrightarrow{x}}^{\prime }\cdot \overrightarrow{n}}{r}}})\approx \frac{1}{r}$

.

$e^{-ik\overrightarrow{n}\cdot {\overrightarrow{x}}^{\prime }}=\underset{l=1}{\sum ^{\mathrm{\infty }}}\frac{{(-ik)}^l}{l!}{(\overrightarrow{n}\cdot {\overrightarrow{x}}^{\prime })}^l$ konvergiert zudem schnell, da $r^{\prime }\approx d$ und $kd\ll 1$. In der Fernzone erhalten wir daher für das Vektorpotential:

$\overrightarrow{A}\left(\overrightarrow{x}\right)\underset{kr\gg 1}{\approx }\frac{e^{ikr}}{cr}\int d^3{x}^{\prime }\overrightarrow{J}\left({\overrightarrow{x}}^{\prime }\right)e^{-ik\overrightarrow{n}\cdot {\overrightarrow{x}}^{\prime }}\approx \frac{e^{ikr}}{cr}\underset{l=1}{\sum ^{\mathrm{\infty }}}\frac{{(-ik)}^l}{l!}\int d^3{x}^{\prime }\overrightarrow{J}\left({\overrightarrow{x}}^{\prime }\right){(\overrightarrow{n}\cdot {\overrightarrow{x}}^{\prime })}^l$

.

Hierin ist $\frac{e^{ikr}}{r}$ eine auslaufende Kugelwelle.

In niedrigster Ordnung dieser Entwicklung erhalten wir einen elektrischen Dipol. Für die Dipolnäherung (die sowohl in der Nah- als auch der Fernzone gültig ist) gilt wegen $r^{\prime }\approx d$, $kd\ll 1\Rightarrow {(\overrightarrow{n}\cdot {\overrightarrow{x}}^{\prime })}^l\approx 0(l\ne 0)$:

$\overrightarrow{A}\left(\overrightarrow{x}\right)\approx \frac{e^{ikr}}{cr}\int d^3{x}^{\prime }\overrightarrow{J}\left({\overrightarrow{x}}^{\prime }\right)$

.

Das Integral über die Stromdichte werden wir umzuformen versuchen:

$\int d^3{x}^{\prime }{J}_l\left({\overrightarrow{x}}^{\prime }\right)=\int d^3{x}^{\prime }{\hat{e}}_l\cdot \overrightarrow{J}\left({\overrightarrow{x}}^{\prime }\right)=\int d^3{x}^{\prime }\overrightarrow{J}\left({\overrightarrow{x}}^{\prime }\right)\cdot {\overrightarrow{\mathrm{\nabla }}}^{\prime }{x}_l^{\prime }$

$=\int d^3{x}^{\prime }{\overrightarrow{\mathrm{\nabla }}}^{\prime }\cdot \left(\overrightarrow{J}\left({\overrightarrow{x}}^{\prime }\right)x_l^{\prime }\right)-\int d^3{x}^{\prime }{x}_l^{\prime }{\overrightarrow{\mathrm{\nabla }}}^{\prime }\cdot \overrightarrow{J}\left({\overrightarrow{x}}^{\prime }\right)$

.

Hierin verwenden wir den Gauß'schen Integralsatz und die Annahme, dass im Unendlichen die Ströme verschwinden:

$\underset{{ℝ}^3}{\int }{d}^3{x}^{\prime }{\overrightarrow{\mathrm{\nabla }}}^{\prime }\cdot \left(\overrightarrow{J}\left({\overrightarrow{x}}^{\prime }\right)x_l^{\prime }\right)=\underset{\partial {ℝ}^3}{\int }d\overrightarrow{a}\cdot \overrightarrow{J}\left({\overrightarrow{x}}^{\prime }\right)x_l^{\prime }⟶0$

.

Dann machen wir noch von der Kontinuitätsgleichung Gebrauch:

$\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{J}\left(\overrightarrow{x}\right)e^{-i\omega t}=\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{J}(t,\overrightarrow{x})=$ $-{\partial }_t\rho (t,\overrightarrow{x})=-{\partial }_t\left(\rho \left(\overrightarrow{x}\right)e^{-i\omega t}\right)=i\omega \rho \left(\overrightarrow{x}\right)e^{-i\omega t}$

,

d.h. $\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{J}\left(\overrightarrow{x}\right)=i\omega \rho \left(\overrightarrow{x}\right)$.

Somit ergibt sich schließlich

$\int d^3{x}^{\prime }\overrightarrow{J}\left({\overrightarrow{x}}^{\prime }\right)=-i\omega \int d^3{x}^{\prime }\rho \left({\overrightarrow{x}}^{\prime }\right){\overrightarrow{x}}^{\prime }=-i\omega \overrightarrow{p}$

,

wobei wir das sog. Dipolmoment

$\overrightarrow{p}=\int d^3{x}^{\prime }\rho \left({\overrightarrow{x}}^{\prime }\right){\overrightarrow{x}}^{\prime }$

(das offensichtlich konstant in t und $\overrightarrow{x}$ ist) eingeführt haben. Mit der Definition

$k=\frac{2\pi }{\lambda }=\frac{\omega }{c}$

resultiert daher für das Vektorpotential:

$\overrightarrow{A}\left(\overrightarrow{x}\right)=-ik\overrightarrow{p}\frac{e^{ikr}}{r}$

.

Hieraus können wir wieder auf gewohnte Weise die magnetische Flussdichte bestimmen:

$\overrightarrow{B}=\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{A}=k^2(\overrightarrow{n}\times \overrightarrow{p})\frac{e^{ikr}}{r}(1-\frac{1}{ikr})$

,

wegen $\overrightarrow{\mathrm{\nabla }}r=\frac{\overrightarrow{x}}{r}=\overrightarrow{n}$, $\overrightarrow{\mathrm{\nabla }}\frac{1}{r}=-\frac{\overrightarrow{x}}{r^3}$.

Weit entfernt von den Strahlungsquellen, d.h. $\overrightarrow{J}\approx \overrightarrow{0}$, und bei Abwesenheit von Dielektrika oder Magnetika, d.h. u.a. $\overrightarrow{B}=\overrightarrow{H}$, sowie einer gleichermaßen harmonische Zeitabhängigkeit von $\overrightarrow{E}$, gilt:

$\overrightarrow{E}=\frac{i}{k}\overrightarrow{\mathrm{\nabla }}\times \overrightarrow{B}=k^2(\overrightarrow{n}\times \overrightarrow{p})\times \overrightarrow{n}\frac{e^{ikr}}{r}+(3\overrightarrow{n}(\overrightarrow{n}\cdot \overrightarrow{p})-\overrightarrow{p})(\frac{1}{r^3}-\frac{ik}{r^2})e^{ikr}$

• In der Fernzone, d.h. $kr\gg 1\Rightarrow \frac{1}{kr}\ll 1$,

wo jedoch $\frac{e^{ikr}}{kr}\ne 0$ gilt, da $e^{ikr}$ schneller wächst als $\frac{1}{kr}$ fällt, erhalten wir daher

$\overrightarrow{B}=k^2(\overrightarrow{n}\times \overrightarrow{p})\frac{e^{ikr}}{r}$, $\overrightarrow{E}=k^2(\overrightarrow{n}\times \overrightarrow{p})\times \overrightarrow{n}\frac{e^{ikr}}{r}=\overrightarrow{B}\times \overrightarrow{n}$.

Für die in ein Raumwinkelelement $d\mathrm{\Omega }$ abgestrahlte Leistung

$dP=(\overrightarrow{S}\cdot \overrightarrow{n})r^{2}d\mathrm{\Omega }$

mit dem Poynting-Vektor

$\overrightarrow{S}=\frac{c}{8\pi }\mathrm{\Re }(\overrightarrow{E}\left(\overrightarrow{x}\right)\times {\overrightarrow{B}}^{\star }\left(\overrightarrow{x}\right))$

ergibt sich mittels

$\overrightarrow{n}\cdot [((\overrightarrow{n}\times \overrightarrow{p})\times \overrightarrow{n})\times (\overrightarrow{n}\times \overrightarrow{p})]={[(\overrightarrow{n}\times \overrightarrow{p})\times \overrightarrow{n}]}^2$

und $\overrightarrow{p}\cdot \overrightarrow{n}=\left|\overrightarrow{p}\right|\cos \vartheta$:

$\frac{dP}{d\mathrm{\Omega }}=\frac{c}{8\pi }{k}^4{\left|\overrightarrow{p}\right|}^2{\sin }^2\vartheta$.

Dies über den gesamten Raumwinkel integriert, resultiert in der gesamten abgestrahlten Leistung

$P=\int d\mathrm{\Omega }\frac{dP}{d\mathrm{\Omega }}=\frac{c}{8\pi }{k}^4{\left|\overrightarrow{p}\right|}^{2}2\pi \underset{-1}{\stackrel{1}{\int }}d\cos \vartheta (1-{\cos }^2\vartheta )=\frac{c{k}^4}{3}{\left|\overrightarrow{p}\right|}^2$.

• In der Nahzone, d.h. $kr\ll 1\Rightarrow \frac{1}{kr}\gg 1,e^{ikr}\approx 1$,

erhalten wir hingegen die üblichen Dipol-Formeln, die wir auch noch in der Elektrostatik herleiten werden:

$\overrightarrow{B}=ik(\overrightarrow{n}\times \overrightarrow{p})\frac{1}{r^2}$, $\overrightarrow{E}=(3\overrightarrow{n}(\overrightarrow{n}\cdot \overrightarrow{p})-\overrightarrow{p})\frac{1}{r^3}$.

Hierin überwiegt das statische elektrische Dipolfeld $\overrightarrow{E}\sim \frac{1}{r^3}$ gegenüber $\overrightarrow{B}\sim \frac{1}{r^2}$.

Referenzen

1. F. Scheck, Theoretische Physik 1: Mechanik. Von den Newton'schen Gesetzen zum deterministischen Chaos (Springer, 2007).
2. A. Sommerfeld, Band III: Elektrodynamik (Harri Deutsch, 1988).
3. A. Lindner, Grundkurs Theoretische Physik (Teubner, 1997).
4. J.D. Jackson, Klassische Elektrodynamik (de Gruyter, 1985).