Die Magnetostatik ist ein Teilgebiet der Elektrodynamik. Sie behandelt magnetische Gleichfelder, also zeitlich konstante Magnetfelder.
In der Magnetostatik wird die räumliche Verteilung von Magnetfeldern in der Umgebung von Dauermagneten und von stationären Strömen (Konzept des Stromfadens) untersucht. Ein stationärer Strom ist beispielsweise Gleichstrom in einem elektrischen Leiter. Hierzu gehören neben den einzelnen magnetischen Eigenschaften der Stoffe wie Ferromagnetismus, Diamagnetismus etc. auch das Erdmagnetfeld. Außerdem beschreibt die Magnetostatik die Kraftwirkung derartig erzeugter Felder auf Magnete und Ströme. Hierzu gehört das Verhalten eines magnetischen Dipols in einem zeitlich konstanten Magnetfeld, beispielsweise das Verhalten einer (frei beweglichen) Magnetnadel im Erdmagnetfeld.
Die Grundbegriffe sind der Elektrostatik analog. Der positiven und negativen elektrischen Ladung entsprechen Nordpole und Südpole, quantitativ: positive und negative Polstärke. Allerdings können magnetische Pole im Gegensatz zu elektrischen Ladungen nicht isoliert werden, sondern treten in einem Körper immer zusammen auf.
Das grundlegende Kraftgesetz der Magnetostatik und somit des stationären Elektromagnetismus ist das Graßmann-Ampère’sche Kraftgesetz, welches wie folgt angegeben werden kann:[1]
![{\displaystyle {\vec {F}}_{m}^{1\rightarrow 2}={\frac {\mu _{0}I_{1}I_{2}}{4\pi }}\oint _{\ell _{1}}\oint _{\ell _{2}}\mathrm {d} {\vec {\ell }}_{2}\times \left(\mathrm {d} {\vec {\ell }}_{1}\times {\frac {{\vec {r}}_{2}-{\vec {r}}_{1}}{|{\vec {r}}_{2}-{\vec {r}}_{1}|^{3}}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/12bfc9a188f95fb12070faa976eeb668c60b5fb9)
In dieser Form beschreibt es die Kraft, welche von der geschlossenen Leiterschleife 1 auf die geschlossene Leiterschleife 2 wirkt. Für den Fall, dass die Leiterschleifen endliche Querschnittflächen besitzen ist folgende Form mit den Stromdichten
und
zu verwenden:
![{\displaystyle {\vec {F}}_{m}^{1\rightarrow 2}={\frac {\mu _{0}}{4\pi }}\int _{V_{1}}\int _{V_{2}}{\vec {\jmath }}_{2}({\vec {r}}_{2})\times \left({\vec {\jmath }}_{1}({\vec {r}}_{1})\times {\frac {{\vec {r}}_{2}-{\vec {r}}_{1}}{|{\vec {r}}_{2}-{\vec {r}}_{1}|^{3}}}\right)\mathrm {d} ^{3}r_{2}\mathrm {d} ^{3}r_{1}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d4cf10dc6ac3fc21beab6e7848248e6bdcf713ce)
Diese Formulierung versteht sich derart, dass die Kraft herrührend von der Stromschleife 1 (mit Volumen
und elektrischer Stromdichte
) auf die Stromschleife 2 (mit Volumen
und elektrischer Stromdichte
) wirkt. Die Volumen
and
sind nicht überlappend. Es wird weiterhin davon ausgegangen, dass es sich um geschlossene Stromschleifen mit endlichen Querschnittsflächen handelt. Diese Annahme impliziert, dass sich die Summe über alle Stromdichte-Vektoren entlang einer jeweiligen Schleife aufhebt:[1]
![{\displaystyle \int _{V_{i}}{\vec {\jmath }}_{i}({\vec {r}})\;\mathrm {d} ^{3}r_{i}={\vec {0}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/394b86b5cab92f5847bebbafc03b5f1bb346c7ee)
(Um dies zu veranschaulichen, betrachte eine gleichförmige Kreisbewegung: Die Summe aller Geschwindigkeitsvektoren in diesem Fall beträgt ebenfalls gleich dem Nullvektor.)
Unter Verwendung der Graßmann-Identität und der Eigenschaft von geschlossenen Kurvenintegralen im Kontext von Gradientenfeldern zeigt sich, dass der effektive Beitrag zur Kraft einfacher anzugeben ist. Zunächst findet sich mit der Graßmann-Identität:[1]
![{\displaystyle \mathrm {d} {\vec {\ell }}_{2}\times \left(\mathrm {d} {\vec {\ell }}_{1}\times {\frac {{\vec {r}}_{2}-{\vec {r}}_{1}}{|{\vec {r}}_{2}-{\vec {r}}_{1}|^{3}}}\right)=-\left(\mathrm {d} {\vec {\ell }}_{2}\cdot \mathrm {d} {\vec {\ell }}_{1}\right){\frac {{\vec {r}}_{2}-{\vec {r}}_{1}}{|{\vec {r}}_{2}-{\vec {r}}_{1}|^{3}}}+\mathrm {d} {\vec {\ell }}_{2}\left(\mathrm {d} {\vec {\ell }}_{1}\cdot {\frac {{\vec {r}}_{2}-{\vec {r}}_{1}}{|{\vec {r}}_{2}-{\vec {r}}_{1}|^{3}}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fe96e03d5912df6bbf5d56e10e492b311ab82ca9)
bzw.
![{\displaystyle {\vec {\jmath }}_{2}\times \left({\vec {\jmath }}_{1}\times {\frac {{\vec {r}}_{2}-{\vec {r}}_{1}}{|{\vec {r}}_{2}-{\vec {r}}_{1}|^{3}}}\right)=-\left({\vec {\jmath }}_{2}\cdot {\vec {\jmath }}_{1}\right){\frac {{\vec {r}}_{2}-{\vec {r}}_{1}}{|{\vec {r}}_{2}-{\vec {r}}_{1}|^{3}}}+{\vec {\jmath }}_{2}\left({\vec {\jmath }}_{1}\cdot {\frac {{\vec {r}}_{2}-{\vec {r}}_{1}}{|{\vec {r}}_{2}-{\vec {r}}_{1}|^{3}}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a8fdaf6f60637a1fc0274cbc511c98400edf85b3)
Die Eigenschaft der Integration entlang geschlossener Kurven führt dazu, dass sich jeweils letzterer Term zu Null ergibt. Folglich kann die Kraft auch einfacher beschrieben werden durch:[1]
![{\displaystyle {\vec {F}}_{m}^{1\rightarrow 2}=-{\frac {\mu _{0}I_{1}I_{2}}{4\pi }}\oint _{\ell _{1}}\oint _{\ell _{2}}\left(\mathrm {d} {\vec {\ell }}_{2}\cdot \mathrm {d} {\vec {\ell }}_{1}\right){\frac {{\vec {r}}_{2}-{\vec {r}}_{1}}{|{\vec {r}}_{2}-{\vec {r}}_{1}|^{3}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/37df525a1a93b7ab670c44b835c87a120a167ea3)
bzw.
![{\displaystyle {\vec {F}}_{m}^{1\rightarrow 2}=-{\frac {\mu _{0}}{4\pi }}\int _{V_{1}}\int _{V_{2}}\left({\vec {\jmath }}_{2}({\vec {r}}_{2})\cdot {\vec {\jmath }}_{1}({\vec {r}}_{1})\right){\frac {{\vec {r}}_{2}-{\vec {r}}_{1}}{|{\vec {r}}_{2}-{\vec {r}}_{1}|^{3}}}\mathrm {d} ^{3}r_{2}\mathrm {d} ^{3}r_{1}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ba10a56cd15aef7538f1a7763bffda4c9796a787)
In dieser Form des Kraftgesetzes wird deutlich, dass das 3. Newton'sche Axiom erfüllt ist (
). Die Formulierung mit Kreuzprodukt gibt diese Eigenschaft nicht explizit, sondern nur unter der Voraussetzung von geschlossenen Stromlinien. Die differentielle Form hat daher in diesem Sinn keine physikalische Aussagekraft[1], was jedoch kein Problem darstellt, da in realen Systemen elektrische Stromkreise in der Regel geschlossen sind. Die Begründung für die komplizierter anmutende Formulierung mit Kreuzprodukten liegt darin, dass so die Definition der magnetischen Flussdichte durch das Biot-Savart-Gesetz eleganter durchführbar ist. Bei solcher Formulierung des magnetostatischen Kraftgesetzes gibt es demnach den Freiheitsgrad eines Gradientenfeldes, welches nach der Integration keinen effektiven Beitrag liefert. In der ursprünglichen Formulierung von Ampere wurde z. B. eine andere „Eichung“ gewählt[2].
Weiterhin ist zu erwähnen, dass die Formulierung des Kraftgesetzes mit Stromdichten als allgemeiner anzusehen ist. Anhand dieser Formulierung schließt sich der Bogen vom Graßmann-Ampère’schen Kraftgesetzt zu den Maxwell-Gleichungen der Magnetostatik. Aus dem Kraftgesetzt lässt sich das Biot-Savart-Gesetz entnehmen, und schließlich lässt sich vom Biot-Savart-Gesetz auf die magnetostatischen Maxwell-Gleichungen für die magnetische Flussdichte schließen, in denen die elektrische Stromdichte enthalten ist und nicht etwa ein Linienstrom.
Aus letzterer Formulierung der magnetischen Kraft geht hervor, dass bei parallel gerichteten Stromdichten die Kraftwirkung anziehend ist, wohingegen bei antiparallel gerichteten Stromdichten eine abstoßende Kraft wirkt. Verglichen mit dem Coulomb’schen Kraftgesetz der Elektrostatik ist dies ein umgekehrtes Wirkungsprinzip. In der Elektrostatik ziehen sich ungleichnamige Ladungen an und gleichnamige stoßen sich ab.
Es zeigt sich, dass das Ampère’sche Kraftgesetz konservativ ist, was bedeutet, dass eine zugehörige potentielle Energie definiert werden kann. Zur Herleitung stellt man zunächst wie folgt um:
![{\displaystyle {\begin{aligned}{\vec {F}}_{m}^{1\rightarrow 2}&=-{\frac {\mu _{0}}{4\pi }}\int _{V_{1}}\int _{V_{2}}\left({\vec {\jmath }}_{2}({\vec {r}}_{2})\cdot {\vec {\jmath }}_{1}({\vec {r}}_{1})\right){\frac {{\vec {r}}_{2}-{\vec {r}}_{1}}{|{\vec {r}}_{2}-{\vec {r}}_{1}|^{3}}}\mathrm {d} ^{3}r_{2}\mathrm {d} ^{3}r_{1}\\&={\frac {\mu _{0}}{4\pi }}\int _{V_{1}}\int _{V_{2}}\left({\vec {\jmath }}_{2}({\vec {r}}_{2})\cdot {\vec {\jmath }}_{1}({\vec {r}}_{1})\right)\nabla _{{\vec {r}}_{2}}\left({\frac {1}{|{\vec {r}}_{2}-{\vec {r}}_{1}|}}\right)\mathrm {d} ^{3}r_{2}\mathrm {d} ^{3}r_{1}\\&={\frac {\mu _{0}}{4\pi }}\int _{V_{1}}\int _{V_{2}}{\frac {\partial }{\partial {\vec {r}}_{2}^{T}}}\left({\frac {{\vec {\jmath }}_{1}({\vec {r}}_{1})}{|{\vec {r}}_{2}-{\vec {r}}_{1}|}}\right)\cdot {\vec {\jmath }}_{2}({\vec {r}}_{2})\mathrm {d} ^{3}r_{2}\mathrm {d} ^{3}r_{1}\\&=\int _{V_{2}}{\frac {\partial }{\partial {\vec {r}}_{2}^{T}}}\left({\frac {\mu _{0}}{4\pi }}\int _{V_{1}}{\frac {{\vec {\jmath }}_{1}({\vec {r}}_{1})}{|{\vec {r}}_{2}-{\vec {r}}_{1}|}}\mathrm {d} ^{3}r_{1}\right)\cdot {\vec {\jmath }}_{2}({\vec {r}}_{2})\mathrm {d} ^{3}r_{2}\\&=\int _{V_{2}}\left.{\frac {\partial {\vec {A}}_{1}}{\partial {\vec {r}}^{T}}}\right|_{{\vec {r}}={\vec {r}}_{2}}\cdot {\vec {\jmath }}_{2}({\vec {r}}_{2})\mathrm {d} ^{3}r_{2}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e5e26a33214cb8cd66688a7b3d31fa0166e00581)
Hier wird das magnetische Vektorpotential
eingeführt:
![{\displaystyle {\vec {A}}_{1}({\vec {r}})={\frac {\mu _{0}}{4\pi }}\int _{V_{1}}{\frac {{\vec {\jmath }}_{1}({\vec {r}}')}{|{\vec {r}}-{\vec {r}}'|}}\mathrm {d} ^{3}r'}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ffba7b9e166e9226b219c4cd464a7cb9b1c0420c)
Wobei die Notation
der Jacobi-Matrix entspricht. Zur Bemessung der magnetischen Kraft an verschiedenen Orten führt man den Verschiebungsvektor
wie folgt ein, sodass sich die Ableitung vor das Integral ziehen lässt:
![{\displaystyle {\begin{aligned}{\vec {F}}_{m}^{1\rightarrow 2}({\vec {\xi }})&=\int _{V_{2}}\left.{\frac {\partial {\vec {A}}_{1}}{\partial {\vec {r}}^{T}}}\right|_{{\vec {r}}={\vec {r}}_{2}+{\vec {\xi }}}\cdot {\vec {\jmath }}_{2}({\vec {r}}_{2})\mathrm {d} ^{3}r_{2}\\&=\nabla _{\vec {\xi }}\int _{V_{2}}{\vec {A}}_{1}({\vec {r}}_{2}+{\vec {\xi }})\cdot {\vec {\jmath }}_{2}({\vec {r}}_{2})\mathrm {d} ^{3}r_{2}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c11586c30b63e1b919cf6d4bac1393684b62bb31)
Entsprechend der Gradienten-Beziehung von Kraft und potentieller Energie
, folgt daraus schließlich die potentielle Energie:
![{\displaystyle {\begin{aligned}E_{\mathrm {pot} }^{12}({\vec {\xi }})&=-\int _{V_{2}}{\vec {A}}_{1}({\vec {r}}_{2}+{\vec {\xi }})\cdot {\vec {\jmath }}_{2}({\vec {r}}_{2})\mathrm {d} ^{3}r_{2}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d4bb677b1287ada6150b40aa682191fc0d001e7d)
bzw. unabhängig von
:
![{\displaystyle {\begin{aligned}E_{\mathrm {pot} }^{12}&=-\int _{V_{2}}{\vec {\jmath }}_{2}({\vec {r}}_{2})\cdot {\vec {A}}_{1}({\vec {r}}_{2})\mathrm {d} ^{3}r_{2}=-{\frac {\mu _{0}}{4\pi }}\int _{V_{2}}\int _{V_{1}}{\frac {{\vec {\jmath }}_{1}({\vec {r}}_{1}){\vec {\jmath }}_{2}({\vec {r}}_{2})}{|{\vec {r}}_{2}-{\vec {r}}_{1}|}}\mathrm {d} ^{3}r_{2}\mathrm {d} ^{3}r_{1}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/682491d46b046586088a2f925686277b13dc7d7c)
Auch hier ist bis auf das negative Vorzeichen die Analogie zur potentiellen Energie in der Elektrostatik erkennbar.
Durch Definition der magnetischen Flussdichte gemäß des Biot-Savart-Gesetz:[1][3]
,
überführt sich die weiter oben angegebene Kraft wie folgt:[1][3]
![{\displaystyle {\vec {F}}_{m}^{1\rightarrow 2}={\frac {\mu _{0}}{4\pi }}\int _{V_{1}}\int _{V_{2}}{\vec {\jmath }}_{2}({\vec {r}}_{2})\times \left({\vec {\jmath }}_{1}({\vec {r}}_{1})\times {\frac {{\vec {r}}_{2}-{\vec {r}}_{1}}{|{\vec {r}}_{2}-{\vec {r}}_{1}|^{3}}}\right)\mathrm {d} ^{3}r_{2}\mathrm {d} ^{3}r_{1}\quad \rightarrow \quad {\vec {F}}_{m}^{1\rightarrow 2}=\int _{V_{2}}{\vec {\jmath }}_{2}({\vec {r}})\times {\vec {B}}_{1}({\vec {r}})\;\mathrm {d} ^{3}r}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a4ba850b7f6327dfbd6cd24de7b1fa7152a061f4)
Wird die Stromdichte als die Kombination von elektrischer Ladungsdichte
und einem Geschwindigkeitsfeld
gemäß
beschrieben, so findet sich hier der direkte Zusammenhang mit der Lorentzkraft im engeren Sinne (für elektrische Punktladung:
).
Die magnetostatische Kraft wirkend auf ein Volumen
mit Magnetisierung
welches sich in einem Magnetfeld
herrührend von einer externen Ursache befindet ist beschrieben durch:[3]
![{\displaystyle {\vec {F}}_{m}=\int _{V}[{\vec {M}}({\vec {r}})\cdot \nabla ]{\vec {B}}_{\mathrm {ext} }({\vec {r}})\mathrm {d} ^{3}r}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9a71ae4ad082b0b765b49c53034efdfc5ac308b7)
Explizit versteht sich diese Notation wie folgt:
![{\displaystyle [{\vec {M}}({\vec {r}})\cdot \nabla ]{\vec {B}}_{\mathrm {ext} }({\vec {r}})=\left[M_{x}({\vec {r}}){\frac {\partial }{\partial x}}+M_{y}({\vec {r}}){\frac {\partial }{\partial y}}+M_{z}({\vec {r}}){\frac {\partial }{\partial z}}\right]{\vec {B}}_{\mathrm {ext} }({\vec {r}})={\begin{bmatrix}M_{x}({\vec {r}}){\frac {\partial B_{x,\mathrm {ext} }}{\partial x}}+M_{y}({\vec {r}}){\frac {\partial B_{x,\mathrm {ext} }}{\partial y}}+M_{z}({\vec {r}}){\frac {\partial B_{x,\mathrm {ext} }}{\partial z}}\\M_{x}({\vec {r}}){\frac {\partial B_{y,\mathrm {ext} }}{\partial x}}+M_{y}({\vec {r}}){\frac {\partial B_{y,\mathrm {ext} }}{\partial y}}+M_{z}({\vec {r}}){\frac {\partial B_{y,\mathrm {ext} }}{\partial z}}\\M_{x}({\vec {r}}){\frac {\partial B_{z,\mathrm {ext} }}{\partial x}}+M_{y}({\vec {r}}){\frac {\partial B_{z,\mathrm {ext} }}{\partial y}}+M_{z}({\vec {r}}){\frac {\partial B_{z,\mathrm {ext} }}{\partial z}}\end{bmatrix}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4537bf0c3b0954ace6eb24c208ea874f23b73c9b)
Bei dieser Formulierung wird davon ausgegangen, dass innerhalb des Volumen die Rotation der magnetischen Flussdichte gleich Null ist (
). Sind elektrische Ströme innerhalb des Volumens
vorhanden, so genügt die hier angegebene Formulierung nicht mehr[3].
Wird ein magnetisiertes Volumen
mit Magnetisierung
von einer von extern eingeprägten Stromdichte
durchflossen, so ist die auf das magnetisierte Volumen wirkende Kraft beschrieben durch:[1][3]
,
wobei
.
Obwohl es keine isolierten magnetischen Ladungen (magnetische Monopole) gibt, können magnetostatische Effekte mit einer Analogie zur Elektrostatik veranschaulicht werden. Dies wird insbesondere in der Schulphysik benutzt: man betrachtet einen Stabmagneten der Länge l als zwei entgegensetzte magnetische Ladungen im Abstand l. Das Analogon zur elektrischen Ladung ist die magnetische Polstärke
. Die Polstärke ist so definiert, dass das magnetische Kraftgesetz (auch: magnetostatisches Kraftgesetz) analog zur Coulomb-Kraft formuliert werden kann:[Anm. 1]
![{\displaystyle F={\frac {1}{4\pi \mu _{0}}}{\frac {p_{1}\cdot p_{2}}{r^{2}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b0fb59e47fd0900be1e3f0de986d0ecdbcbced53)
F ist hierbei die magnetische Kraft, die zwischen zwei Magnetpolen der Polstärke
und
im Abstand
wirkt; μ0 ist die magnetische Feldkonstante. Die Polstärke ist von der gleichen Dimension wie der magnetische Fluss und wird somit in der Einheit Weber angegeben.[Anm. 1]
Aus der Definition folgt z. B. bei einem homogenen Feld mit bekannter Flussdichte B und Fläche A für die Kraft:
![{\displaystyle F={\frac {1}{4\pi \mu _{0}}}{\frac {\Phi _{1}\cdot \Phi _{2}}{r^{2}}}={\frac {1}{4\pi \mu _{0}}}{\frac {B_{1}A_{1}\cdot B_{2}A_{2}}{r^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/794e9d3dd04492242d2f6a0a82e6b87d5a264283)
Die Feldtheorie der Magnetostatik ist didaktisch auf im Wesentlichen zwei unterschiedliche Arten zu vermitteln. Der historische (oder induktive) Ansatz geht zunächst von der Naturbeobachtung von Kräften aus, deren Wirkungsprinzip folglich anhand von Experimenten tiefer untersucht, und schlussendlich in einem physikalischen Kraft-Gesetz auf eine möglichst einfache und allgemeine Weise mathematisch formuliert wird. Durch weitere mathematische Untersuchung des Kraftgesetzes im Sinne der Trennung von Ursache und Wirkung finden sich dann die Feldgleichungen (Maxwell-Gleichungen). Der deduktive Ansatz geht davon aus, dass die Maxwell-Gleichungen bereits bekannt sind.
Ausgegangen wird vom Graßmann-Ampere’schen Kraftgesetz:[1]
,
bei welchem in dieser Form die elektrische Stromdichte
als Ursache der magnetischen Kraft angesehen wird, welche auf die elektrische Stromdichte
wirkt. An dieser Stelle ist die Miteinbeziehung der Ladungserhaltung von wichtiger Bedeutung[1], welche sich in der Divergenzfreiheit der Stromdichten niederschlägt (
). wird die magnetische Flussdichte
entnommen. Diese Form der magnetischen Flussdichte wird als Biot-Savart-Gesetz bezeichnet:
![{\displaystyle {\vec {B}}_{1}({\vec {r}})={\frac {\mu _{0}}{4\pi }}\int _{V_{1}}{\vec {\jmath }}_{1}({\vec {r}}')\times {\frac {{\vec {r}}-{\vec {r}}'}{|{\vec {r}}-{\vec {r}}'|^{3}}}\mathrm {d} ^{3}r'}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c168fecd24ee26a5a209d9bdcc31b2a6861e3f90)
Somit schreibt sich das Graßmann-Ampere’sche Kraftgesetz wie folgt:
![{\displaystyle {\vec {F}}_{m}^{1\rightarrow 2}=\int _{V_{2}}{\vec {\jmath }}_{2}({\vec {r}})\times {\vec {B}}_{1}({\vec {r}})\mathrm {d} ^{3}r}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fae895fb0defdf585e63ee6c4cc756e57f42fa84)
Die weitere Untersuchung widmet sich nun dem Biot-Savart-Gesetz (zur Einfachheit verzichten wir auf den Index). Es zeigt sich zunächst, dass das Biot-Savart-Gesetz sich ebenfalls darstellen lässt durch:[1]
![{\displaystyle {\vec {B}}({\vec {r}})={\frac {\mu _{0}}{4\pi }}\nabla \times \int _{V}{\frac {{\vec {\jmath }}({\vec {r}}')}{|{\vec {r}}-{\vec {r}}'|}}\mathrm {d} ^{3}r'}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ee0e93c1a808dc5a588364aac5a1cdb020e0f042)
Dies impliziert bereits die Divergenzfreiheit der magnetischen Flussdichte:
.
Die Rotation der magnetischen Flussdichte findet sich ebenfalls:[1]
![{\displaystyle {\begin{aligned}\nabla \times {\vec {B}}({\vec {r}})&={\frac {\mu _{0}}{4\pi }}\nabla \times \nabla \times \int _{V}{\frac {{\vec {\jmath }}({\vec {r}}')}{|{\vec {r}}-{\vec {r}}'|}}\mathrm {d} ^{3}r'\\&={\frac {\mu _{0}}{4\pi }}\nabla \nabla \cdot \int _{V}{\frac {{\vec {\jmath }}({\vec {r}}')}{|{\vec {r}}-{\vec {r}}'|}}\mathrm {d} ^{3}r'-{\frac {\mu _{0}}{4\pi }}\Delta \int _{V}{\frac {{\vec {\jmath }}({\vec {r}}')}{|{\vec {r}}-{\vec {r}}'|}}\mathrm {d} ^{3}r'\\&={\frac {\mu _{0}}{4\pi }}\nabla \int _{V}{\frac {[\nabla _{{\vec {r}}'}\cdot {\vec {\jmath }}({\vec {r}}')]}{|{\vec {r}}-{\vec {r}}'|}}\mathrm {d} ^{3}r'-{\frac {\mu _{0}}{4\pi }}\int _{V}\Delta _{\vec {r}}\left({\frac {1}{|{\vec {r}}-{\vec {r}}'|}}\right){\vec {\jmath }}({\vec {r}}')\mathrm {d} ^{3}r'\\&={\frac {\mu _{0}}{4\pi }}\int _{V}(4\pi \delta ({\vec {r}}-{\vec {r}}')){\vec {\jmath }}({\vec {r}}')\mathrm {d} ^{3}r'\\&=\mu _{0}{\vec {\jmath }}({\vec {r}})\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9eda5915225c1634a1311cc5a00157360ebfc25c)
Die Rotation der magnetischen Flussdichte entspricht folglich:
.
Man beachte, dass in dieser Herleitung die Ladungserhaltung gemäß
, die Laplace-Identität
, und die Identität
verwendet wurden.
Zusammengefasst erhält man die Maxwellgleichungen der Magnetostatik in differentieller Form:
![{\displaystyle \nabla \cdot {\vec {B}}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/77d6b768c165d3d213f72b3204ae62d39bab1be2)
![{\displaystyle \nabla \times {\vec {B}}=\mu _{0}{\vec {\jmath }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d011d86a0332b209a419847d82e43ba02dd181bb)
Definiert man weiterhin das magnetische Vektorpotential:
,
sodass
, so wird aus der obigen Herleitung klar, dass die Divergenz des magnetischen Vektorpotentials ebenfalls verschwinden muss. Die Coulomb-Eichung des magnetischen Vektorpotentials, welche im nächsten Abschnitt bei der Herleitung der magnetischen Flussdichte aus den Maxwell-Gleichungen Verwendung findet, ist somit in dieser Herangehensweise eine Schlussfolgerung aus der Ladungserhaltung. Für das Vektorpotential gilt schließlich die Poisson-Gleichung
.
Für zeitlich konstante Felder „entkoppeln“ die Gleichungen für elektrische (E) und magnetische (B) Felder: setzt man in den Maxwellgleichungen alle Zeitableitungen gleich 0, so entstehen Gleichungen, die nicht gleichzeitig E und B enthalten. Die Phänomene der Magnetostatik lassen sich mit folgenden zwei reduzierten Maxwell-Gleichungen beschreiben:
![{\displaystyle \nabla \cdot {\vec {B}}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/77d6b768c165d3d213f72b3204ae62d39bab1be2)
![{\displaystyle \nabla \times {\vec {B}}=\mu _{0}{\vec {\jmath }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d011d86a0332b209a419847d82e43ba02dd181bb)
Beachtenswert ist, dass in dieser stationären Betrachtung (zeitlich konstant, statisch) die Divergenz der elektrische Stromdichte
stets gleich Null ist. Dieser Aspekt steht im Zusammenhang mit der Kontinuitätsgleichung
für die elektrische Stromdichte, wobei
die elektrische Volumenladungsdichte repräsentiert. D.h., im hier betrachteten statischen Fall gilt:
.
Man führt das Vektorpotential
als Hilfsfeld mit folgender Definition ein:
![{\displaystyle {\vec {B}}=\nabla \times {\vec {A}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7e7407e07bdb01a1048a72518d136819c32434f5)
Dadurch wird automatisch die Gleichung
erfüllt, da die Divergenz eines Rotationsfeldes identisch 0 ist
.
ist jedoch nicht eindeutig bestimmt, da
invariant ist unter einer Eichtransformation
mit
. D. h. die durch A und A’ festgelegten B-Felder sind identisch. Dies ergibt sich aus
,
da die Rotation des Gradienten eines Skalarfeldes verschwindet.
Setzt man
in die inhomogene Maxwellgleichung (obige Gleichung 2)
![{\displaystyle \mu _{0}{\vec {\jmath }}=\nabla \times \nabla \times {\vec {A}}=\nabla \left({\nabla \cdot {\vec {A}}}\right)-\Delta {\vec {A}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f2fea8434b53051f27372e89195be3071b353504)
ein (
ist der Laplace-Operator), so ergibt sich mit der Coulomb-Eichung
die besonders einfache Form:
![{\displaystyle \Delta {\vec {A}}=-\mu _{0}{\vec {\jmath }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5d3b6be94e9bdba556d1d2943597ff06248ce526)
Dies stellt für jede Komponente eine Poisson-Gleichung dar, die durch
![{\displaystyle {\vec {A}}({\vec {r}})={\frac {\mu _{0}}{4\pi }}\int _{\mathbb {R} ^{3}}\mathrm {d} ^{3}r'{\frac {{\vec {\jmath }}({\vec {r}}\,')}{|{\vec {r}}-{\vec {r}}\,'|}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3d824ca477b9c8a7a2366a1401b41b940d16e885)
gelöst wird.
Wendet man die Rotation auf A an so erhält man das Biot-Savart-Gesetz für das physikalisch relevante B-Feld
![{\displaystyle {\vec {B}}\left({\vec {r}}\right)=\nabla \times {\vec {A}}\left({\vec {r}}\right)={\frac {\mu _{0}}{4\pi }}\int _{\mathbb {R} ^{3}}{\nabla _{\vec {r}}\times {\frac {{\vec {\jmath }}\left({{\vec {r}}\,'}\right)}{\left|{{\vec {r}}-{\vec {r}}\,'}\right|}}}\mathrm {d} ^{3}r'={\frac {\mu _{0}}{4\pi }}\int _{\mathbb {R} ^{3}}{\left({\nabla _{\vec {r}}{\frac {1}{\left|{{\vec {r}}-{\vec {r}}\,'}\right|}}}\right)\times }{\vec {\jmath }}\left({{\vec {r}}\,'}\right)\mathrm {d} ^{3}r'={\frac {\mu _{0}}{4\pi }}\int _{\mathbb {R} ^{3}}{{\vec {\jmath }}\left({{\vec {r}}\,'}\right)\times {\frac {{\vec {r}}-{\vec {r}}\,'}{\left|{{\vec {r}}-{\vec {r}}\,'}\right|^{3}}}\mathrm {d} ^{3}r'}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8ad6f040bc064f0a0b50701b1833c2e89ea40bd9)
Für einen Stromfaden geht
zu
über:
![{\displaystyle {\vec {B}}\left({\vec {r}}\right)={\frac {\mu _{0}}{4\pi }}I\int _{s}{\mathrm {d} {\vec {s}}\,'\times {\frac {{\vec {r}}-{\vec {r}}\,'}{\left|{{\vec {r}}-{\vec {r}}\,'}\right|^{3}}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/50cbba141e90c8a3423644bfa2d2d73e226e8c43)
Allgemeiner stellt sich das Biot-Savart-Gesetz mit dem Faltungsprodukt „
“ (Faltungsintegral) dar. Die Verwendung der Faltungsalgebra vereinfacht die Darstellungen und ermöglicht einfachere Termumformungen. Das magnetische Vektorpotential und die magnetische Flussdichte stellen sich unter Verwendung des Faltungsproduktes wie folgt dar:
![{\displaystyle {\vec {A}}({\vec {r}})={\frac {\mu _{0}}{4\pi }}\left[{\frac {1}{|{\vec {r}}|}}\ast {\vec {\jmath }}\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a5255e43b3e285362f6f3c733a817c01c89b4f42)
![{\displaystyle {\vec {B}}({\vec {r}})={\frac {\mu _{0}}{4\pi }}\nabla \times \left[{\frac {1}{|{\vec {r}}|}}\ast {\vec {\jmath }}\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f95bdb80dbfd48928388a60ed4d3a0c58d54dcf7)
Durch korrekte Anwendung der Faltungsalgebra und zugehörigen Ableitungsregeln ist es zum Beispiel trivial zu zeigen, dass die Divergenz des magnetischen Vektorpotential verschwindet:
![{\displaystyle \nabla \cdot {\vec {A}}({\vec {r}})={\frac {\mu _{0}}{4\pi }}\nabla \cdot \left[{\frac {1}{|{\vec {r}}|}}\ast {\vec {\jmath }}\right]={\frac {\mu _{0}}{4\pi }}\left[{\frac {1}{|{\vec {r}}|}}\ast \nabla \cdot {\vec {\jmath }}\right]=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/333414ec627fc82000e1b9faa9c793378f9f5265)
(Man bemerke, dass die Divergenz der elektrischen Stromdichte aufgrund der Ladungserhaltung gleich Null ist.)
Eine alternative Formulierung, welche im Kontext von, z. B., Permanentmagneten von hohem Interesse ist, wird durch die Einführung der Magnetisierung
erreicht. Die Magnetisierung lässt sich in der statischen Betrachtung durch die Kontinuumgleichung für die elektrische Stromdichte einführen. Für die stationäre elektrische Stromdichte gilt die Gleichung
, welche durch Einführung der Magnetisierung durch
erfüllt ist. Für den Fall, dass die Magnetisierung im Bezug auf gebundene elektrische Ströme verwendet wird, muss weiterhin das folgende Integral verschwinden:[1]
.
Dies ist das hinreichende Integral-Kriterium für geschlossene, nicht-divergierende Stromlinien, wohingegen die differentielle Gleichung
eine notwendige Bedingung darstellt. Mit der Magnetisierung lassen sich die magnetostatischen Maxwell-Gleichungen wie folgt angeben:
![{\displaystyle \nabla \cdot {\vec {B}}=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/77d6b768c165d3d213f72b3204ae62d39bab1be2)
![{\displaystyle \nabla \times {\vec {B}}=\mu _{0}\nabla \times {\vec {M}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/74572b0e481db48538f1d8d8f563666cac22b6a5)
Analog zur Formulierung mit elektrischer Stromdichte lässt sich unter Berücksichtigung der Coulomb-Eichung nachfolgende Poisson-Gleichung für das magnetische Vektorpotential angeben:
![{\displaystyle \Delta {\vec {A}}=-\mu _{0}\nabla \times {\vec {M}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7f1bac3712f5bebc45fcf30b991454afecebed82)
An dieser Stelle bietet es sich an ein weiteres Vektorpotential
, der magnetische Hertz-Vektor (benannt nach Heinrich Hertz), gemäß
einzuführen. Dies ermöglicht es wiederum eine einfache Poisson-Gleichung zu formulieren:[4][5]
![{\displaystyle \Delta {\vec {A}}=-\mu _{0}\nabla \times {\vec {M}}\quad \rightarrow \quad \nabla \times \left(\Delta {\vec {\Pi }}_{m}+\mu _{0}{\vec {M}}\right)={\vec {0}}\quad \rightarrow \quad \Delta {\vec {\Pi }}_{m}=-\mu _{0}{\vec {M}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ef3d25afa90aa3786de6ef9827c09a54bcf29894)
Analog wie bei der Formulierung mit elektrischer Stromdichte lässt sich die Lösung mit folgendem Faltungs-Integral angeben:[4]
![{\displaystyle {\vec {\Pi }}_{m}({\vec {r}})={\frac {\mu _{0}}{4\pi }}\int _{\mathbb {R} ^{3}}\mathrm {d} ^{3}r'{\frac {{\vec {M}}({\vec {r}}')}{|{\vec {r}}-{\vec {r}}'|}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d13505c4cc42ac37d0a276406ba0f4b79ab93043)
Die zu diesem magnetischen Hertz-Vektor gehörigen Felder finden sich durch differentiation. Diese können wie folgt angegeben werden:
![{\displaystyle {\vec {A}}({\vec {r}})=\nabla \times {\vec {\Pi }}_{m}=-{\frac {\mu _{0}}{4\pi }}\int _{\mathbb {R} ^{3}}{\frac {{\vec {r}}-{\vec {r}}'}{|{\vec {r}}-{\vec {r}}'|^{3}}}\times {\vec {M}}({\vec {r}}')\mathrm {d} ^{3}r'}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6466a51006bbaf61769ce74f4154613b907df57d)
![{\displaystyle {\vec {B}}({\vec {r}})=\nabla \times \nabla \times {\vec {\Pi }}_{m}=\nabla \times {\vec {A}}={\frac {\mu _{0}}{4\pi }}\int _{\mathbb {R} ^{3}}\left[3{\frac {({\vec {r}}-{\vec {r}}')\otimes ({\vec {r}}-{\vec {r}}')}{|{\vec {r}}-{\vec {r}}'|^{5}}}-{\frac {I_{3}}{|{\vec {r}}-{\vec {r}}'|^{3}}}+{\frac {8\pi }{3}}I_{3}\delta ({\vec {r}}-{\vec {r}}')\right]{\vec {M}}({\vec {r}}')\mathrm {d} ^{3}r'}](https://wikimedia.org/api/rest_v1/media/math/render/svg/aa4a88df45bcceaab72560cb6479db7660c7f050)
![{\displaystyle \varphi _{m}({\vec {r}})=-\mu _{0}^{-1}\nabla \cdot {\vec {\Pi }}_{m}={\frac {1}{4\pi }}\int _{\mathbb {R} ^{3}}{\frac {{\vec {r}}-{\vec {r}}'}{|{\vec {r}}-{\vec {r}}'|^{3}}}\cdot {\vec {M}}({\vec {r}}')\mathrm {d} ^{3}r'}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6addfd6a457453748ed80096b1783b121e3a158d)
![{\displaystyle {\vec {H}}({\vec {r}})=\mu _{0}^{-1}\nabla \nabla \cdot {\vec {\Pi }}_{m}=-\nabla \varphi _{m}={\frac {1}{4\pi }}\int _{\mathbb {R} ^{3}}\left[3{\frac {({\vec {r}}-{\vec {r}}')\otimes ({\vec {r}}-{\vec {r}}')}{|{\vec {r}}-{\vec {r}}'|^{5}}}-{\frac {I_{3}}{|{\vec {r}}-{\vec {r}}'|^{3}}}-{\frac {4\pi }{3}}I_{3}\delta ({\vec {r}}-{\vec {r}}')\right]{\vec {M}}({\vec {r}}')\mathrm {d} ^{3}r'}](https://wikimedia.org/api/rest_v1/media/math/render/svg/96a812521c126a8127f0db84862d87f93180d40d)
Hier stellt „
“ das dyadische Produkt (äußeres Produkt) und
die Einheitsmatrix dar. Aus diesen Lösungen folgt ebenfalls die Beziehung von Magnetisierung
, magnetischer Feldstärke
und magnetischer Flussdichte
gemäß
. Ein Vorteil der Darstellung mit dyadischem Produkt ist die Möglichkeit der Separierung von Integrations-Kernel und Quelle. So wird offensichtlicher, dass sich alle obigen Formulierungen auch allgemeiner durch ein Faltungsprodukt darstellen lassen. Ein Vorteil der Verwendung der Faltungsalgebra liegt in der Übersichtlichkeit von Umformungen. So finden sich leichter alternative Darstellungen der obigen Integrale. Mit Faltungsprodukt (Faltungsintegral) „
“ stellen sich die Felder wie folgt dar:
![{\displaystyle {\vec {B}}({\vec {r}})={\frac {\mu _{0}}{4\pi }}\nabla \times \nabla \times \left[{\frac {1}{|{\vec {r}}|}}\ast {\vec {M}}\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bd57932be3fc3df00c956c6f4df09cd7a25f0ab0)
![{\displaystyle {\vec {H}}({\vec {r}})={\frac {1}{4\pi }}\nabla \nabla \cdot \left[{\frac {1}{|{\vec {r}}|}}\ast {\vec {M}}\right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d451701d49a83abc54c910297a6a91007ce80fd5)
Aufgrund der Kommutativität des Faltungsproduktes ist so offensichtlich, dass die Felder auf acht (bzw. 16 wenn man berücksichtigt, dass durch
von
auf
und umgekehrt geschlossen werden kann) unterschiedliche Weisen berechnet werden können. Die obigen Varianten sind solche, bei denen die Nabla-Operatoren ausschließlich auf die Greensche Funktion
angewandt werden. Eine weitere sehr übliche Variante findet sich durch Anwendung der inneren Nabla-Operationen auf die Magnetisierung:[1]
![{\displaystyle {\vec {B}}({\vec {r}})={\frac {\mu _{0}}{4\pi }}\nabla \times \nabla \times \left[{\frac {1}{|{\vec {r}}|}}\ast {\vec {M}}\right]={\frac {\mu _{0}}{4\pi }}\nabla \times \left[{\frac {1}{|{\vec {r}}|}}\ast \nabla \times {\vec {M}}\right]={\frac {\mu _{0}}{4\pi }}\nabla \times \int _{\mathbb {R} ^{3}}{\frac {[\nabla _{{\vec {r}}'}\times {\vec {M}}({\vec {r}}')]}{|{\vec {r}}-{\vec {r}}'|}}\mathrm {d} ^{3}r'}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9a5dfbccdddfefc10a667135a20ae3b3cebad657)
![{\displaystyle {\vec {H}}({\vec {r}})={\frac {1}{4\pi }}\nabla \nabla \cdot \left[{\frac {1}{|{\vec {r}}|}}\ast {\vec {M}}\right]={\frac {1}{4\pi }}\nabla \left[{\frac {1}{|{\vec {r}}|}}\ast \nabla \cdot {\vec {M}}\right]={\frac {1}{4\pi }}\nabla \int _{\mathbb {R} ^{3}}{\frac {[\nabla _{{\vec {r}}'}\cdot {\vec {M}}({\vec {r}}')]}{|{\vec {r}}-{\vec {r}}'|}}\mathrm {d} ^{3}r'}](https://wikimedia.org/api/rest_v1/media/math/render/svg/71f673a7885c45e04d0e267a0cd3a147e2901fba)
An dieser Stelle ist die fiktive magnetische Volumenladungsdichte
als Quelle der magnetischen Feldstärke ersichtlich. Ist die Magnetisierung
auf ein endliches Volumen
begrenzt so müssen Beiträge von Volumen und Oberfläche unterschieden werden:[1]
![{\displaystyle {\vec {B}}({\vec {r}})={\frac {\mu _{0}}{4\pi }}\nabla \times \int _{V}{\frac {[\nabla _{{\vec {r}}'}\times {\vec {M}}({\vec {r}}')]}{|{\vec {r}}-{\vec {r}}'|}}\mathrm {d} ^{3}r'+{\frac {\mu _{0}}{4\pi }}\nabla \times \oint _{\partial V}{\frac {[{\vec {M}}({\vec {r}}')\times {\vec {n}}']}{|{\vec {r}}-{\vec {r}}'|}}d^{2}r'}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4483cb533d9e9bbb9e2e52a487c5cb9040568e62)
![{\displaystyle {\vec {H}}({\vec {r}})={\frac {1}{4\pi }}\nabla \int _{V}{\frac {[\nabla _{{\vec {r}}'}\cdot {\vec {M}}({\vec {r}}')]}{|{\vec {r}}-{\vec {r}}'|}}\mathrm {d} ^{3}r'-{\frac {1}{4\pi }}\nabla \oint _{\partial V}{\frac {[{\vec {n}}'\cdot {\vec {M}}({\vec {r}}')]}{|{\vec {r}}-{\vec {r}}'|}}d^{2}r'}](https://wikimedia.org/api/rest_v1/media/math/render/svg/89b5349c05a8f7a856870f0fd79ae2e2719a16b1)
Die hier auftretenden Oberflächen-Terme berücksichtigen die Beiträge von fiktiven magnetischen Oberflächenladungen
und Oberflächenströmen
.
(Genau genommen: Beim Übergang von einem unendlich ausgedehnten Magnetisierungsvolumen zu einem endlichen Volumen entsteht eine Unstetigkeit an der Grenze zwischen Magnetisierungsvolumen und Vakuum. Somit ist die Magnetisierung allgemeiner mit Hilfe einer Indikatorfunktion
gemäß
zu beschreiben, wobei
mindestens einfach differenzierbar ist. Durch diese Beschreibung von
ist offensichtlich, dass eine Produktregel greift. Die Anwendung der Produktregel und durch den Einsatz der Delta-Distribution findet sich so der korrekte Übergang von der Beschreibung mit
zur Beschreibung mit endlichem Volumen
. In der Regel wird dies in Physik Büchern nicht derart ausführlich behandelt, was bei komplizierteren Problemstellungen jedoch zu Verwirrungen führen kann. Im Buch von J.D. Jackson[1] findet sich daher die Anmerkung: "Never combine the surface integral of
with (5.98)!". Wobei Gleichung (5.98) im Buch von J.D. Jackson der allgemeineren Formulierung mit
entspricht.)
Wie im Abschnitt zum Ampere'schen Kraftgesetz hergeleitet, ist die zum Ampere’schen Gesetz zugehörige potentielle Energie gegeben durch:
![{\displaystyle {\begin{aligned}E_{\mathrm {pot} }^{12}&=-{\frac {\mu _{0}}{4\pi }}\int _{V_{2}}\int _{V_{1}}{\frac {{\vec {\jmath }}_{1}({\vec {r}}_{1}){\vec {\jmath }}_{2}({\vec {r}}_{2})}{|{\vec {r}}_{2}-{\vec {r}}_{1}|}}\mathrm {d} ^{3}r_{2}\mathrm {d} ^{3}r_{1}=-\int _{V_{2}}{\vec {\jmath }}_{2}({\vec {r}}_{2})\cdot {\vec {A}}_{1}({\vec {r}}_{2})\mathrm {d} ^{3}r_{2}=-\int _{V_{1}}{\vec {\jmath }}_{1}({\vec {r}}_{1})\cdot {\vec {A}}_{2}({\vec {r}}_{1})\mathrm {d} ^{3}r_{2}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/95e7020670b0672f914691c0bb30601356faab47)
Bis auf das negative Vorzeichen der Energie ist die Analogie zur potentiellen Energie in der Elektrostatik erkennbar. Das negative Vorzeichen deutet an, dass parallel gerichtete Stromdichten sich gegenseitig anziehen und antiparallel gerichtete Stromdichten sich gegenseitig Abstoßen. Umgekehrt ist es in der Elektrostatik, wo gleichnamige Ladungen sich gegenseitig abstoßen und ungleichnamige Ladungen sich gegenseitig anziehen.
Ausgehend vom vorigen Abschnitt stellt sich die magnetostatische potentielle Energie einer geschlossenen Stromschleife mit Stromdichte
, befindlich in einem externen Magnetfeld mit magnetischem Vektorpotential
wie folgt dar:
![{\displaystyle E_{\mathrm {pot} }=-\int _{V}{\vec {\jmath }}({\vec {r}})\cdot {\vec {A}}_{\mathrm {ext} }({\vec {r}})\mathrm {d} ^{3}r}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f817248a5a95147a05eba2a6952441cdd9c98f70)
In der Magnetostatik ist die Formulierung mit Stromdichte und Magnetisierung als äquivalent anzusehen. Die Grundannahme für die Stromdichte ist
, welche durch die Einführung der Magnetisierung gemäß
erfüllt ist. Das Magnetisierungs-Vektorfeld
kann also als „Vektorpotential“ der Stromdichte
angesehen werden (analog zu
). Im Weiteren gilt es einen Ausdruck für die potentielle Energie in Abhängigkeit der Magnetisierung zu finden. Man beginnt mit:
![{\displaystyle {\begin{aligned}E_{\mathrm {pot} }&=-\int _{\mathbb {R} ^{3}}{\vec {\jmath }}({\vec {r}})\cdot {\vec {A}}_{\mathrm {ext} }({\vec {r}})\mathrm {d} ^{3}r\\&=-\int _{\mathbb {R} ^{3}}[\nabla \times {\vec {M}}]({\vec {r}})\cdot {\vec {A}}_{\mathrm {ext} }({\vec {r}})\mathrm {d} ^{3}r\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/351e19bd7b75aedc9613557f728d1c6c6a5e558d)
(Man beachte, dass hier die Formulierung mit
anstelle
wichtig ist. Dies ist dadurch gerechtfertigt, dass die Stromdichte
mathematisch stets so formuliert werden kann, dass diese außerhalb eines endlichen Volumens
gleich Null ist.) Durch Verwendung der Identität
, durch Anwendung des Gauß'schen Integralsatzes und der Beziehung
findet sich weiterhin:
![{\displaystyle {\begin{aligned}E_{\mathrm {pot} }&=-\int _{\mathbb {R} ^{3}}{\vec {M}}\cdot \nabla \times {\vec {A}}_{\mathrm {ext} }+\nabla \cdot ({\vec {M}}\times {\vec {A}}_{\mathrm {ext} })\,\mathrm {d} ^{3}r\\&=-\int _{\mathbb {R} ^{3}}{\vec {M}}\cdot \nabla \times {\vec {A}}_{\mathrm {ext} }\,\mathrm {d} ^{3}r+\oint _{\partial \mathbb {R} ^{3}}({\vec {M}}\times {\vec {A}}_{\mathrm {ext} })\cdot {\vec {n}}\,\mathrm {d} ^{2}r\\&=-\int _{\mathbb {R} ^{3}}{\vec {M}}({\vec {r}})\cdot {\vec {B}}_{\mathrm {ext} }({\vec {r}})\,\mathrm {d} ^{3}r\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/640382ceb738695d03745578227bd2fe89d14bd9)
Hier ist offensichtlich, dass das Oberflächenintegral verschwinden muss, da das magnetische Vektorpotential und die Magnetisierung im Unendlichen verschwinden. Dies zeigt, dass die Formulierungen mit Stromdichte und Magnetisierung identisch sind:[1][3][6]
![{\displaystyle {\begin{aligned}E_{\mathrm {pot} }&=-\int _{\mathbb {R} ^{3}}{\vec {\jmath }}({\vec {r}})\cdot {\vec {A}}_{\mathrm {ext} }({\vec {r}})\,\mathrm {d} ^{3}r=-\int _{\mathbb {R} ^{3}}{\vec {M}}({\vec {r}})\cdot {\vec {B}}_{\mathrm {ext} }({\vec {r}})\,\mathrm {d} ^{3}r\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/330cb4e201ffec9ebad44d1eee72e3f43beb46b3)
Die potentielle magnetostatische Energie eines magnetisierten Volumens
mit Magnetisierungs-Vektorfeld
, befindlich in einem von einer externen Ursache herrührenden Magnetfeldes
ist wie folgt anzugeben:[1][3][6]
![{\displaystyle E_{\mathrm {pot} }=-\int _{V}{\vec {M}}({\vec {r}})\cdot {\vec {B}}_{\mathrm {ext} }({\vec {r}})\;\mathrm {d} ^{3}r}](https://wikimedia.org/api/rest_v1/media/math/render/svg/525c96dc50be560de4232c9a4f93f7b1d47767db)
Im Fall der Interaktion zweier Dauermagnete, mit nicht-überlappenden endlichen Volumina
und
, mit Magnetisierungen
und
, ist das externe Magnetfeld formal bekannt und die potentielle magnetostatische Energie der Wechselwirkung ist wie folgt anzugeben:
![{\displaystyle E_{\mathrm {pot} }^{12}=-{\frac {\mu _{0}}{4\pi }}\int _{V_{1}}\int _{V_{2}}{\vec {M}}_{1}({\vec {r}}_{1})\cdot \left[3{\frac {({\vec {r}}_{1}-{\vec {r}}_{2})\otimes ({\vec {r}}_{1}-{\vec {r}}_{2})}{|{\vec {r}}_{1}-{\vec {r}}_{2}|^{5}}}-{\frac {I_{3}}{|{\vec {r}}_{1}-{\vec {r}}_{2}|^{3}}}+{\frac {8\pi }{3}}I_{3}\delta ({\vec {r}}_{1}-{\vec {r}}_{2})\right]\cdot {\vec {M}}_{2}({\vec {r}}_{2})\;\mathrm {d} ^{3}r_{2}\mathrm {d} ^{3}r_{1}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c9f02a34757abd0485dd3b7db6beff64b6d87ba5)
Entsprechend des 3. Newton'schen Axioms ist bei dieser Formulierung erkennbar, dass die potentielle Energie invariant gegenüber der Vertauschung der Indizes
und
ist.
Die Selbstenergie einer geschlossenen Stromschleife mit endlicher Querschnittsfläche ist gegeben durch:[1][7]
![{\displaystyle E_{\mathrm {pot} }^{\mathrm {self} }=-{\frac {1}{2}}\int _{V}{\vec {\jmath }}({\vec {r}})\cdot {\vec {A}}({\vec {r}})\;\mathrm {d} ^{3}r}](https://wikimedia.org/api/rest_v1/media/math/render/svg/de240e99847fe90c42622276ab788888ec4f7198)
Da es sich hier um die Selbstenergie handelt gilt
. Damit lässt sich alternative schreiben:
![{\displaystyle {\begin{aligned}E_{\mathrm {pot} }^{\mathrm {self} }&=-{\frac {1}{2}}\int _{V}{\vec {\jmath }}({\vec {r}})\cdot {\vec {A}}({\vec {r}})\;\mathrm {d} ^{3}r\\&=+{\frac {1}{2\mu _{0}}}\int _{\mathbb {R} ^{3}}\Delta {\vec {A}}({\vec {r}})\cdot {\vec {A}}({\vec {r}})\;\mathrm {d} ^{3}r\\&=-{\frac {1}{2\mu _{0}}}\int _{\mathbb {R} ^{3}}\sum _{k\in \{x,y,z\}}(\nabla A_{k})^{2}\mathrm {d} ^{3}r+{\frac {1}{2\mu _{0}}}\oint _{\partial \mathbb {R} ^{3}}\sum _{k\in \{x,y,z\}}(\nabla A_{k}A_{k}\cdot {\vec {n}})\;d^{2}r\\&=-{\frac {1}{2\mu _{0}}}\int _{\mathbb {R} ^{3}}\sum _{k\in \{x,y,z\}}(\nabla A_{k})^{2}\mathrm {d} ^{3}r\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b481e54e1967419465c7af1c833af25927b2119b)
Hier wurde die Produktregel
und der Gauß'sche Integralsatz verwendet. Der Oberflächenterm ergibt sich zu Null, da in unendlicher Entfernung von der Ursache das magnetischen Vektorpotential gegen Null strebt. Alternativ findet sich durch Verwendung der Identität
und unter Verwendung der verschwindenden Divergenz
ein weiterer Ausdruck für die Selbstenergie. Zunächst formen wir um:
![{\displaystyle {\begin{aligned}E_{\mathrm {pot} }^{\mathrm {self} }&=-{\frac {1}{2}}\int _{V}{\vec {\jmath }}({\vec {r}})\cdot {\vec {A}}({\vec {r}})\;\mathrm {d} ^{3}r\\&=+{\frac {1}{2\mu _{0}}}\int _{\mathbb {R} ^{3}}\Delta {\vec {A}}({\vec {r}})\cdot {\vec {A}}({\vec {r}})\;\mathrm {d} ^{3}r\\&=-{\frac {1}{2\mu _{0}}}\int _{\mathbb {R} ^{3}}(\nabla \times \nabla \times {\vec {A}})\cdot {\vec {A}}\;\mathrm {d} ^{3}r\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/38e873568ba87c58c6456d7665b6bb3b10a07305)
An dieser Stelle ist die Produktregel
anzuwenden. Wir erhalten:
![{\displaystyle {\begin{aligned}E_{\mathrm {pot} }^{\mathrm {self} }&=-{\frac {1}{2\mu _{0}}}\int _{\mathbb {R} ^{3}}(\nabla \times {\vec {A}})\cdot (\nabla \times {\vec {A}})+\nabla \cdot [(\nabla \times {\vec {A}})\times {\vec {A}}]\;\mathrm {d} ^{3}r\\&=-{\frac {1}{2\mu _{0}}}\int _{\mathbb {R} ^{3}}(\nabla \times {\vec {A}})\cdot (\nabla \times {\vec {A}})\,\mathrm {d} ^{3}r-{\frac {1}{2\mu _{0}}}\oint _{\partial \mathbb {R} ^{3}}[(\nabla \times {\vec {A}})\times {\vec {A}}]\cdot {\vec {n}}\;\mathrm {d} ^{2}r\\&=-{\frac {1}{2\mu _{0}}}\int _{\mathbb {R} ^{3}}(\nabla \times {\vec {A}})\cdot (\nabla \times {\vec {A}})\,\mathrm {d} ^{3}r\\&=-{\frac {1}{2\mu _{0}}}\int _{\mathbb {R} ^{3}}|{\vec {B}}|^{2}\mathrm {d} ^{3}r\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/84abf5fc4b8f89f19e4998f6c9008514a63f5b26)
Neben der obigen Produktregel wurde weiterhin der Gauß'sche Integralsatz verwendet. Der Oberflächen-Term ergibt sich zu Null unter der Annahme, dass das Vektorpotential im unendlichen gegen null konvergiert. Abschließend ist die Rotation des Vektorpotentials durch die magnetische Flussdichte ersetzbar (
).
Die potentielle magnetostatische Selbstenergie eines magnetisierten Volumens
mit Magnetisierungs-Vektorfeld
ist wie folgt anzugeben:
![{\displaystyle E_{\mathrm {pot} }^{\mathrm {self} }=-{\frac {1}{2}}\int _{V}{\vec {M}}({\vec {r}})\cdot {\vec {B}}({\vec {r}})\;\mathrm {d} ^{3}r}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2cd28a9a7a0e3621dd742d99116072adf77bdb74)
Der Unterschied zur Formulierung für externe Felder besteht darin, dass hier die magnetische Flussdichte
zugehörig zur Magnetisierung
ist (es gilt
). Durch Verwendung der im Abschnitt Feldtheorie angegebenen Beziehung
finden sich alternative Formulierungen. Zunächst lässt sich feststellen, dass sich die Magnetisierung durch
ausdrückbar ist. Mit dieser Erkenntnis findet sich:
![{\displaystyle E_{\mathrm {pot} }^{\mathrm {self} }=-{\frac {1}{2\mu _{0}}}\int _{\mathbb {R} ^{3}}|{\vec {B}}({\vec {r}})|^{2}\;\mathrm {d} ^{3}r+{\frac {1}{2}}\int _{\mathbb {R} ^{3}}{\vec {H}}({\vec {r}})\cdot {\vec {B}}({\vec {r}})\;\mathrm {d} ^{3}r}](https://wikimedia.org/api/rest_v1/media/math/render/svg/96a9f2ca873e7625c921a15be0e625997efa2ee4)
Durch weitere Untersuchung findet sich, dass der zweite Term gleich Null ist:[1][6][8]
![{\displaystyle \int _{\mathbb {R} ^{3}}{\vec {H}}({\vec {r}})\cdot {\vec {B}}({\vec {r}})\;\mathrm {d} ^{3}r=-\int _{\mathbb {R} ^{3}}\nabla \varphi _{m}\cdot {\vec {B}}\;\mathrm {d} ^{3}r=-\int _{\mathbb {R} ^{3}}\nabla \cdot (\varphi _{m}{\vec {B}})-\varphi _{m}\nabla \cdot {\vec {B}}\;\mathrm {d} ^{3}r=-\oint _{\partial \mathbb {R} ^{3}}(\varphi _{m}{\vec {B}})\cdot {\vec {n}}\;\mathrm {d} ^{2}r=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e9a0b264a1fae1736ee1acd132b099ece032d428)
Um dies zu zeigen ist die Produktregel
, der Gaußsche Integralsatz und die Divergenzfreiheit der magnetischen Flussdichte
zu verwenden. Das im letzten Schritt auftretende Hüllenintegral muss ebenfalls verschwinden, da die Felder
und
bei unendlicher Entfernung zu ihrer Ursache gegen Null streben. Folglich findet sich:
![{\displaystyle E_{\mathrm {pot} }^{\mathrm {self} }=-{\frac {1}{2}}\int _{V}{\vec {M}}({\vec {r}})\cdot {\vec {B}}({\vec {r}})\;\mathrm {d} ^{3}r=-{\frac {1}{2\mu _{0}}}\int _{\mathbb {R} ^{3}}|{\vec {B}}({\vec {r}})|^{2}\;\mathrm {d} ^{3}r}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b4c1c5c181575aab048525462ca29388d5bf8330)
In der hier aufgeführten Formulierung ist eine Feinheit bezüglich des Integrationsgebietes zu beachten. Bei letzterem Ausdruck muss über den kompletten Raum integriert werden, da die magnetische Flussdichte in der Regel in jedem Punkt
einen von Null verschiedenen Wert hat. Bei ersterem Ausdruck wird hier von einem lokal begrenztem Volumen
ausgegangen, sodass die Magnetisierung
außerhalb dieses Volumens nicht vorhanden ist. Alternativ könnte man diese räumliche Begrenzung dem Vektorfeld
durch Verwendung einer Indikatorfunktion überlassen, und allgemeiner das Integrationsgebiet
durch
ersetzen. Eine entsprechende Indikatorfunktion würde folglich innerhalb des Volumens
den Wert 1 und außerhalb des Volumens
den Wert 0 annehmen.
Im Kontext des Mikromagnetismus ist eine Grundannahme, dass das Magnetisierungs-Vektorfeld eine konstante Magnitude, die Sättigungs-Magnetisierung
, besitzt. Die Magnetisierung wird folglich durch
mit
beschrieben. Unter dieser Annahme findet sich:[8]
![{\displaystyle {\begin{aligned}E_{\mathrm {pot} }^{\mathrm {self} }&=-{\frac {1}{2\mu _{0}}}\int _{\mathbb {R} ^{3}}|{\vec {B}}({\vec {r}})|^{2}\;\mathrm {d} ^{3}r\\&=-{\frac {\mu _{0}}{2}}\int _{\mathbb {R} ^{3}}({\vec {H}}+{\vec {M}})\cdot ({\vec {H}}+{\vec {M}})\;\mathrm {d} ^{3}r\\&=-{\frac {\mu _{0}}{2}}\int _{\mathbb {R} ^{3}}|{\vec {H}}|^{2}+|{\vec {M}}|^{2}+2{\vec {H}}\cdot {\vec {M}}\;\mathrm {d} ^{3}r\\&=-{\frac {\mu _{0}}{2}}\int _{\mathbb {R} ^{3}}|{\vec {H}}|^{2}+|{\vec {M}}|^{2}+2{\vec {H}}\cdot ({\vec {B}}/\mu _{0}-{\vec {H}})\;\mathrm {d} ^{3}r\\&=+{\frac {\mu _{0}}{2}}\int _{\mathbb {R} ^{3}}|{\vec {H}}|^{2}\mathrm {d} ^{3}r-{\frac {\mu _{0}}{2}}\int _{\mathbb {R} ^{3}}|{\vec {M}}|^{2}\,\mathrm {d} ^{3}r\\&=+{\frac {\mu _{0}}{2}}\int _{\mathbb {R} ^{3}}|{\vec {H}}|^{2}\mathrm {d} ^{3}r-{\frac {\mu _{0}}{2}}M_{s}^{2}V\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/74c6215343ac30ea360c6ab2019cc36597874cb3)
Man beachte, dass sich im vierten Schritt das Volumenintegral über das Skalarprodukt
exakt zu Null ergibt (siehe voriger Abschnitt). Der letzte Schritt ist unter der Annahme des Mirkomagnetismus gerechtfertigt. Eine wichtige Aussage, welche anhand dieser Formulierung der magnetostatischen Selbstenergie gewonnen werden kann, wird als Pole-Avoidance Principle bezeichnet. Da es sich bei dem zweiten Term um eine Konstante handelt, nimmt dieser bei Energie-Minimierung keinen Einfluss. Der erste Term ist positiv definit. Dies hat zur Folge, dass bei Energie-Minimierung die magnetische Feldstärke
möglichst gegen Null strebt, was mit der Auslöschung von fiktiven magnetischen Polen einhergeht. Dies ist der Fall wenn das Magnetisierungs-Vektorfeld
einem geschlossenen Feldlinienbild folgt. Daher können in der magnetischen Mikrostruktur von Materialien wirbelartige Magnetisierungs-Vektorfelder auftreten.
Die Magnetostatik kann mit folgender Lagrange-Dichte formuliert werden:
![{\displaystyle {\mathcal {L}}=-{\frac {1}{2\mu _{0}}}\sum _{k\in \{x,y,z\}}(\nabla A_{k})^{2}+{\vec {\jmath }}\cdot {\vec {A}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cbfd341332b36e75a00437a87217dfb01c55e0e8)
Dabei ist der erste Term der Selbstenergie und der zweite Term der negierten potentiellen Energie zuzuordnen. Durch die Anwendung der Euler-Lagrange Gleichungen
![{\displaystyle {\frac {\partial {\mathcal {L}}}{\partial {\vec {A}}}}-\sum _{k\in \{x,y,z\}}{\frac {\partial }{\partial k}}{\frac {\partial {\mathcal {L}}}{\partial (\partial _{k}{\vec {A}})}}=\mathbf {0} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/37c8a611fc54b11a96a4d1ddd278565243e25d14)
lässt sich folglich die Poisson-Gleichung für das magnetische Vektorpotential
herleiten. In diesen Euler-Lagrange Gleichungen versteht sich der erste Ableitungsterm wie ein Gradient, sodass
.
Die partiellen Ableitungen im zweiten Term ergeben sich zunächst gemäß:
,
und folglich:
.
Somit ergibt sich aus den Euler-Lagrange Gleichungen die Poisson-Gleichung:
![{\displaystyle {\frac {\partial {\mathcal {L}}}{\partial {\vec {A}}}}-\sum _{k\in \{x,y,z\}}{\frac {\partial }{\partial k}}{\frac {\partial {\mathcal {L}}}{\partial (\partial _{k}{\vec {A}})}}={\vec {\jmath }}+{\frac {1}{\mu _{0}}}\Delta {\vec {A}}=\mathbf {0} \quad \rightarrow \quad \Delta {\vec {A}}=-\mu _{0}{\vec {\jmath }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/06f934a82d5701f3c26a549caee5786312695d47)
Eine äquivalente Lagrange-Dichte ist gegeben durch:
![{\displaystyle {\mathcal {L}}=-{\frac {1}{2\mu _{0}}}|\nabla \times {\vec {A}}|^{2}+{\vec {\jmath }}\cdot {\vec {A}}=-{\frac {1}{2\mu _{0}}}|{\vec {B}}|^{2}+{\vec {\jmath }}\cdot {\vec {A}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e76ded72536c12a38a22abe9ff5ab7206c6bdfe7)
|
|
|
Die Grundannahme ist eine
divergenzfreie Stromdichte, als
auch geschlossene Stromlinien
|
|
Diese Grundannahmen sind stets im
Kontext zu allen weiteren Gleichungen
zu sehen. Die Stromdichte Vektoren
verstehen sich als Tangentenvektoren an
die Stromlinien.
|
Ampere'sches Kraftgesetz in
Grassmann-Formulierung
|
|
Diese Formulierung des Kraftgesetzes ist üblich,
da es die direkte Verbindung mit der Lorentz-Kraft
explizit darstellt. Die vereinfachte Darstellung ist
äquivalent zu dieser Formulierung
|
Ampere'sches Kraftgesetz in
vereinfachter Darstellung
|
|
Das Ampere'sche Kraftgesetz ist konservativ
und erfüllt das 3. Newton'sche Axiom. Diese
Formulierung enthält die effektiven Beiträge der
Grassmann Formulierung in Kontext der Magnetostatik.
In dieser Formulierung ist explizit ersichtlich, dass parallel
fließende Ströme sich abstoßen.
|
Ampere'sches Potential
|
|
|
Biot-Savart Gesetz
|
|
|
Maxwell-Gleichungen
|
|
|
Poisson-Gleichung
|
|
|
Selbst-Energie
|
|
|
Lagrange-Dichte
|
|
|
- ↑ a b Man findet auch die Definition
. In diesem Fall hat die Polstärke die Dimension „Stromstärke × Länge“ und die Einheit A·m.
Im elektromagnetischen CGS-System gilt einfach:
, und die Polstärke hat die Dimension „√Kraft × Länge“.
- ↑ a b c d e f g h i j k l m n o p q r s t J. D. Jackson: Electrodynamics, Classical. In: digital Encyclopedia of Applied Physics. 15. April 2003.
- ↑ A.K.T. Assis, M.A. Bueno: Equivalence between Ampere and Grassmann's forces. In: IEEE Transactions on Magnetics. Band 32, Nr. 2, März 1996, ISSN 0018-9464, S. 431–436, doi:10.1109/20.486529.
- ↑ a b c d e f g Edward Purcell: Electricity and Magnetism. Cambridge University Press, 2011, ISBN 978-1-107-01360-5.
- ↑ a b E. A. Essex: Hertz vector potentials of electromagnetic theory. In: American Journal of Physics. Band 45, Nr. 11, 1. November 1977, ISSN 0002-9505, S. 1099–1101, doi:10.1119/1.10955.
- ↑ Hertzian electromagnetic potentials and associated gauge transformations. In: Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences. Band 231, Nr. 1185, 22. August 1955, ISSN 0080-4630, S. 250–263, doi:10.1098/rspa.1955.0170 (royalsocietypublishing.org [abgerufen am 19. Dezember 2023]).
- ↑ a b c J. M. D. Coey: Magnetism and Magnetic Materials. Cambridge University Press, Cambridge 2010, ISBN 978-0-521-81614-4, doi:10.1017/cbo9780511845000 (cambridge.org [abgerufen am 27. Dezember 2023]).
- ↑ Oleg D. Jefimenko: Electricity and Magnetism: An Introduction to the Theory of Electric and Magnetic Fields. 2. ed Auflage. Electret Scientific, Star City 1989.
- ↑ a b Amikam Aharoni: Introduction to the Theory of Ferromagnetism. Oxford University Press, 4. Januar 2001.