静磁場のソースを表示

'''静磁場'''（せいじば、{{lang-en|static magnetic field}}）とは、時間的に変動しない磁場のことである
<ref name="pano1">
W.K.H. パノフスキー (著), M. フィリップス (著),林 忠四郎 (翻訳), 西田 稔 (翻訳);「新版 電磁気学〈上〉」吉岡書店; POD版 (2002/09)
</ref>
<ref name="mizoguchi">
[http://books.google.tt/books/about/電磁気学.html?id=69NMAAAACAAJ [★]溝口 正  (著) 「電磁気学―SI UNITS 」裳華房 (2001/03)](特にP188付近)
</ref>
<ref name="takeyama9">
[★]竹山 説三 (著) 「電磁気学現象理論」丸善出版; 3版 (1949)
</ref>
<ref name="sandf">[★]P.P.シルベスタ(著),
R.L.フェラーリ（著）,
本間 利久（著）,
田中 康博（著）;
「有限要素法による電磁界解析 (Information & computing (26)) 」
サイエンス社 (1988/09) <br />
[原書]Peter P. Silvester (Author), Ronald L. Ferrari (Author);
「Finite Elements for Electrical Engineers」Cambridge University Press; 3版 (1996/9/5)
[https://books.google.co.jp/books?id=QXWFEIesd14C&pg=PA478&lpg=PA478&dq=Finite+Elements+for+Electrical+Engineers&source=bl&ots=e6xjqZBZer&sig=IGgeInNnlFg47hAzj7ax8WwN1ZQ&hl=ja&sa=X&ei=jFqzUd7gE87TkgWt6oCYCQ&ved=0CGMQ6AEwBw#v=onepage&q=Finite%20Elements%20for%20Electrical%20Engineers&f=false]
</ref>
<ref name="saito">{{PDFlink|[http://www.saito-lab.jp/papers/Conference/2003MAG-03-55.pdf 斎藤兆古,坂本禎智,藤原耕二；「磁気回路法と有限要素法の理論的関係」]}} 電気学会マグネティックス研究会資料 MAG-03 号：54-64 ページ：5-10 発行年：2003年03月31日</ref>
<ref name="miyata">[http://www.akita-nct.ac.jp/yamamoto/study/thesis/2005/senkoka_paper_miyata 静磁場解析のための二次要素を用いる有限要素法の研究]
</ref><ref name="yugen_2d">[http://computation.cside.com/mag-index.html 有限要素法による磁場解析]</ref>
<ref name="akita2">[http://www.akita-nct.jp/~yamamoto/study/electromagnetics/cylindrical_magnet/html/node3.html#SECTION00032000000000000000 磁場が満たす偏微分方程式（秋田高専講義録）]</ref>
<!--
-->
<ref name="sonodaumi">園田英徳;「大学院生のための基礎物理学」講談社（2011/09/29）ISBN 978-4-06-153277-9</ref>
<ref name="hira3">平川浩正;「電磁気学（新物理学シリーズ２）」培風館 (1986/04) ISBN 9784563024024
</ref>
<ref name="mori4">守末 利弥「数値電磁気学のためのゲージ理論」森北出版 (1996/04) ISBN 978-4-627-71600-1
</ref>
<ref name="yori5">依田 潔 (著) ;「Mathematicaによる電磁界シミュレーション入門 - ＰＯＤ版 (計算電気・電子工学シリーズ)」森北出版(2012/2/24) ISBN 978-4-627-71529-5
</ref>
<ref name="tanaka6">
{{PDFlink|[http://www-het.phys.sci.osaka-u.ac.jp/~tanaka/emIIge/2009/section3_4.pdf 電磁気学II([[大阪大学]] 田中実教授の講義録）]}}
</ref>
<ref name="rikyo">[http://www.ne.rikkyo.ac.jp/~ieki/em/jan11.html][[立教大学]]<span>講義ノート</span>

[http://www.ne.rikkyo.ac.jp/~ieki/em/text.html 本講義の参考文献欄]
等、至る所に[[文字化け]]があるが、たとえば、[http://lab.kiki-verb.com/mojibakeratta/ もじばけらった]等によって解読できる。
</ref>
<ref name="rikadai">
{{PDFlink|[http://www.rs.tus.ac.jp/nikuni/elemag/elemagnotes12.pdf [[東京理科大学]]講義ノート]}}
{{PDFlink|[http://www.rs.tus.ac.jp/nikuni/elemag/elemag12.pdf および、講義スライド]}}
</ref>
<ref name="kagosima7">鹿児島 誠一 (著) ;「電磁気学 (パリティ物理学コース)」丸善(1997/01) ISBN 978-4-621-04277-9</ref>
<ref name="goto8">後藤 憲一(著), 山崎修一郎(著)  ;「詳解 電磁気学演習」共立出版 (1970/12) ISBN 978-4-320-03022-0</ref>
<ref name="maeda10">前田 和茂 (著) ,小林 俊雄(著);「ビジュアルアプローチ電磁気学」森北出版 (2009/12/5) ISBN 978-4-627-16221-1</ref>
<ref name="yousho1">Julius Adams Stratton;「Electromagnetic Theory」Wiley-IEEE Press; 1版 (2007/1/22) ISBN 978-0-470-13153-4</ref>
<ref name="yousyo2">A. Pramanik (著) ;「Electro-Magnetism: Theory and Applications 」Prentice-Hall of India Pvt.Ltd (2004/8/15)
[https://books.google.co.jp/books?id=gnEEwy12S5cC&pg=PT232&dq=uniform+magnetization&hl=ja&sa=X&ei=GdkFUaaMJKP2mAWa9YDABg&ved=0CGYQ6AEwCDgU]
</ref>
<ref name="yousyo3">Bo-nan Jiang(著) ;「The Least-Squares Finite Element Method: Theory and Applications in Computational Fluid Dynamics and Electromagnetics (Scientific Computation) 」Springer; 1998版 (1998/6/22) [https://books.google.co.jp/books?id=bp1csE6-UBoC&pg=PA356&dq=%22vector+potential%22++magnetostatic&hl=ja&sa=X&ei=SRsjUZfjKcmNmQXm5YHQBA&ved=0CD4Q6AEwAg#v=onepage&q=%22vector%20potential%22%20%20magnetostatic&f=false]
</ref>
<ref name="akita1">[http://akita-nct.jp/yamamoto/lecture/2007/p1/index.html 電気磁気学特論（秋田高専講義録）]</ref>
<ref name="takaha">高橋 康人 (早稲田大学):{{PDFlink|[http://www.cc.kyushu-u.ac.jp/publish/kohobkno/genkoVol6No2/15-takahashi2.pdf "高速多重極法を組み込んだ磁気モーメント法による 磁性体解析に関する基礎的検討"]}}</ref>
<ref name="dendai1">東京電機大学 (編集)「入門 電磁気学」東京電機大学出版局 (2006/03) ISBN 9784501004200</ref>
<ref name="densan">早川 義晴 (著) 「電気教科書 電験三種合格ガイド」翔泳社 (2011/2/25) ISBN 9784798126623</ref>
<ref name="jikikairo">
{{PDFlink|[http://www.daido-electronics.co.jp/qa/magnet_qa/documents/magnet_qa10.pdf 磁気回路と電気回路]}}</ref>
<ref name="tuboi">坪井 一洋 (著) 「システムと微分方程式」三恵社 (2011/5/22) ISBN 9784883618248</ref>
<ref name="fuku">[http://irobutsu.a.la9.jp/mybook/ykwkrEM/magnets_f.html 複数磁石による静磁場(簡易シミュレータ)]</ref>
<ref name="kindai">{{PDFlink|[http://www.rist.kindai.ac.jp/~kiguchi/phystext/chap_21.pdf 近畿大学 講義ノート]}}</ref>。

本記事では、静磁気学(Magnetostatics)の視点から、静磁場について述べる。

== 定常な電流が作り出す静磁場の一般論 ==
本節では、真空中に定常な（つまり時刻tに依存しない）電流密度が作り出す磁束密度について、一般に成り立つ事柄について述べる。ただし、時間的な変動の影響はもちろんのこと、これ以外にも、電場や強制電荷、分極電荷の影響は排除されているものとする。
本記事では、専ら体積電流密度を中心に扱い、線電流近似については、例えば<ref name="todoriki">
等々力 二郎：[https://doi.org/10.1585/jspf1958.57.318 「有限の太さの矩形断面ヘリカルコイルの磁場の計算」]核融合研究,57(1987)318</ref>
<ref name="watabe1">[https://doi.org/10.1585/jspf1958.35.235 渡辺 二太:「空心円筒コイル群によって作られる磁場の計算」] 核融合研究 Vol.35 (1976) No.3 P235-242 </ref>
<ref name="watabe2">渡辺 二太:[https://doi.org/10.1585/jspf1958.63.482 「多様な形状のコイルに対する磁場計算法」] 核融合研究 Vol.63 (1990) No.6 P482-507</ref>等に委ねることとする。


真空中に定常な（つまり時刻tに依存しない）電流密度
<math>\boldsymbol{i}(\mathbf{r})</math>が与えられた
とする。このとき、<math>\boldsymbol{i}(\mathbf{r})</math>
は、以下の磁気ベクトルポテンシャル
<math>\mathbf{A}_{\boldsymbol{i}}(\mathbf{r})</math>を空間内に作り出す。

:<math>\mathbf{A}_{\boldsymbol{i}}(\mathbf{r})
=\frac{\mu_{0}}{4\pi}
{\int}_{\mathbf{s}\in\mathbb{R}^{3}}
\left(\frac{\boldsymbol{i}(\mathbf{s})}
{|\mathbf{r}-\mathbf{s}|}
\right)
\ d^{3}\mathbf{s}
</math>    （1-1）
となる。

<math>\mathbf{B}_{\boldsymbol{i}}=\operatorname{rot}[\mathbf{A}_{\boldsymbol{i}}]</math>を考え併せると、
<math>\boldsymbol{i}(\mathbf{r})</math>が直接的に作り出す磁束密度<math>\mathbf{B}_{\boldsymbol{i}}</math>は、

:<math>\mathbf{B}_{\boldsymbol{i}}(\mathbf{r})
=\frac{\mu_{0}}{4 \pi }
\int_{\mathbf{s}\in\mathbb{R}^{3}}
\left(\frac{
\boldsymbol{i}(\mathbf{s}) \times (\mathbf{r} - \mathbf{s})
}{|\mathbf{r} - \mathbf{s}|^3}\right) d^3
\mathbf{s}</math>    （1-2）

となる。これは、すなわち、[[ビオ・サバールの法則]]である。

上記の<math>\mathbf{B}_{\boldsymbol{i}}(\mathbf{r})</math>に対し、
新たな場<math>{\mathbf{H}}_{\boldsymbol{i}}</math>を、
:<math>\mathbf{H}_{\boldsymbol{i}}(\mathbf{r})
:=\frac{1}{{\mu}_{0}}\mathbf{B}_{\boldsymbol{i}}(\mathbf{r})
</math>   (1-3)
と定義する。この場のことを、「電流密度iが作り出す磁場」と呼ぶ。ここでμ<sub>0</sub>は、真空の透磁率である。
尚、定義の上では、「電流密度iが作り出す磁場」<math>{\mathbf{H}}_{\boldsymbol{i}}</math>は、
透磁率がμの場所でも、:<math>\mathbf{H}_{\boldsymbol{i}}(\mathbf{r})
:=\frac{1}{{\mu}_{0}}\mathbf{B}_{\boldsymbol{i}}(\mathbf{r})
</math>であることに特に注意されたい。

式(1-2),式(1-3)より、
:<math>\mathbf{H}_{\boldsymbol{i}}(\mathbf{r})
=\frac{1}{4 \pi}
\int_{\mathbf{s}\in\mathbb{R}^{3}}
\left(\frac{
\boldsymbol{i}(\mathbf{s}) \times (\mathbf{r} - \mathbf{s})
}{|\mathbf{r} - \mathbf{s}|^3}\right) d^3
\mathbf{s}</math>   (1-4)

である。これに、回転微分を作用させると、
:<math>\operatorname{rot}_{\mathbf{r}}[\mathbf{H}_{\boldsymbol{i}}]= \boldsymbol{i}</math>   (1-5)
が得られる。実際、ベクトル解析の公式より、
<blockquote>

:<math>
\begin{align}
\operatorname{rot}_{\mathbf{r}}\left[
\left(\frac{
\boldsymbol{i}(\mathbf{s}) \times (\mathbf{r} - \mathbf{s})
}{|\mathbf{r} - \mathbf{s}|^3}\right)
\right]
&= 4 \pi \left(
\operatorname{div}_{\mathbf{r}}
\left[
\frac{
(\mathbf{r} - \mathbf{s})
}{|\mathbf{r} - \mathbf{s}|^3}
\right]
\right)\boldsymbol{i}(\mathbf{s}) \\
&= 4\pi{\delta}^{3}(\mathbf{r}-\mathbf{s})
\boldsymbol{i}(\mathbf{s})
\qquad \text{(1-6)}
\end{align}
</math>

従って、

:<math>\begin{align}
\operatorname{rot}_{\mathbf{r}}
[\mathbf{H}_{\boldsymbol{i}}(\mathbf{r})]
&= \frac{1}{4\pi}\int_{\mathbf{s}\in\mathbb{R}^{3}}
\operatorname{rot}_{\mathbf{r}}
\left[
\left(\frac{
\boldsymbol{i}(\mathbf{s}) \times (\mathbf{r} - \mathbf{s})
}{|\mathbf{r} - \mathbf{s}|^3}\right)
\right]
d^3 \mathbf{s} \\
&= {\int}_{\mathbf{s}\in\mathbb{R}^{3}}
{\delta}^{3}(\mathbf{r}-\mathbf{s})
\boldsymbol{i}(\mathbf{s})
d^3\mathbf{s} \\
&= \boldsymbol{i}(\mathbf{r})
\qquad \text{(1-7)}
\end{align}
</math>

が得られる。ここでδ<sup>3</sup>は、3変数の[[δ関数]]（ディラックのデルタ）を意味する。
</blockquote>

== 時間的に定常な磁化が作り出す静磁場の一般論 ==
本節では、定常な（つまり時刻tに依存しない）[[磁化]]が作り出す磁束密度について、一般に成り立つ事柄について述べる。ただし、時間的な変動の影響はもちろんのこと、これ以外にも、電場や強制電荷、分極電荷の影響は排除されているものとする。


=== 時間的に定常な磁化が作り出す磁気ベクトルポテンシャル ===
空間内の領域<math>\Omega</math>に物質が置かれ、前記物質が、定常な（つまり時刻tに依存しない）[[磁化]]<math>\mathbf{M}(\mathbf{r})</math>を帯びているとする。このとき、磁化ベクトルは、「単位体積当たりの磁気モーメントの密度」を表すものであるため、磁化の定義より、

:物質内の各点<math>\mathbf{s}\in\Omega</math>それぞれに、それぞれ
::<math>\mathbf{M}(\mathbf{s}){d}^{3}\mathbf{s}</math>  （2-1-1）
:で与えられる磁気モーメントが配置されている
と考えることができる<ref name="mizoguchi" />。

まず、時間的に定常な磁化が作り出す磁気ベクトルポテンシャルについて考えよう。
<BLOCKQUOTE>
(1)'''原点に置かれた'''磁気モーメント'''m'''は、空間上の位置'''r'''に作り出す磁気ベクトルポテンシャルは、
::<math>\mathbf{a}_{0}(\mathbf{r})=\frac{\mu_{0}}{4\pi}\left(\frac{\mathbf{m}\times \mathbf{r}}
{|\mathbf{r}|}\right)</math>  （2-1-2）
:::である。<br />
(2)従って、"(1)"を平行移動すれば、'''位置<math>\mathbf{s}\in\mathbb{R}^{3}</math>に置かれた'''磁気モーメント'''m'''は、空間上の位置'''r'''に作り出す磁気ベクトルポテンシャルは、
::<math>\mathbf{a}_{s}(\mathbf{r})=\frac{\mu_{0}}{4\pi}\left(\frac{\mathbf{m}\times (\mathbf{r}-\mathbf{s})}
{|\mathbf{r}-\mathbf{s}|}\right)</math>  （2-1-3）
:::である。
</BLOCKQUOTE>

従って、上記の磁気モーメント<math>\mathbf{M}(\mathbf{s})d\mathbf{s}</math>それぞれは、磁気ベクトルポテンシャル

:<math>\frac{\mu_{0}}{4\pi}\left(\frac{(\mathbf{M}(\mathbf{s})d\mathbf{s})\times (\mathbf{r}-\mathbf{s})}
{|\mathbf{r}-\mathbf{s}|}\right)
</math>  （2-1-4）

を作り出す。

上記の磁気ベクトルポテンシャルそれぞれを、全ての<math>\mathbf{s}\in\Omega</math>に渡って足し合わせると、

:<math>\mathbf{A}_{M}(\mathbf{r})
=\frac{\mu_{0}}{4\pi}
{\int}_{\mathbf{s}\in\Omega}
\left(\frac{(\mathbf{M}(\mathbf{s}))\times (\mathbf{r}-{\mathbf{s}})}
{|\mathbf{r}-\mathbf{s}|}
\right)
\ {d}^{3}\mathbf{s}
</math>  （2-1-5）

を得る。即ち、物質の磁化<math>\mathbf{M}(\mathbf{r})</math>
は、上記の磁気ベクトルポテンシャル
<math>\mathbf{A}_{M}</math>を空間内に作り出す<ref name="mizoguchi" />
<ref group="注釈">
尚、磁化が<math>\Omega</math>の外では0であるからといって、<math>\Omega</math>の外で<math>\mathbf{A}_{M}</math>が0となるとは限らない（大概の場合は<math>\Omega</math>の外でも<math>\mathbf{A}_{M}</math>は、0ではない）ことに注意されたい。</ref>
<ref group="注釈">
なお、磁化は、<math>\Omega</math>の外では0なので、上式は、
:<math>\mathbf{A}_{M}(\mathbf{r})
=\frac{\mu_{0}}{4\pi}
{\int}_{\mathbf{s}\in\mathbb{R}^{3}}
\left(\frac{(\mathbf{M}(\mathbf{s}))\times (\mathbf{r}-\mathbf {s})}
{|\mathbf{r}-\mathbf{s}|}
\right)
\ {d}^{3}\mathbf{s}
</math>
としてもよい（紛らわしいのであまり推奨しない）。
</ref>

=== 時間的に定常な磁化が作り出す磁束密度 ===
次に、時間的に定常な磁化が作り出す磁束密度について考えよう。磁気ベクトルポテンシャルの回転微分をとれば、磁束密度が得られる。
<BLOCKQUOTE>
(1)'''原点に置かれた'''磁気モーメント'''m'''は、空間上の位置'''r'''に作り出す磁気ベクトルポテンシャルは、前述の通り、
::<math>\mathbf{a}_{0}(\mathbf{r})=\frac{\mu_{0}}{4\pi}\left(\frac{\mathbf{m}\times \mathbf{r}}
{|\mathbf{r}|}\right)</math>  （2-2-1）
:である。これの回転微分をとることで、'''原点に置かれた'''磁気モーメント'''m'''が、空間上の位置'''r'''に作り出す磁束密度は、
::<math>\mathbf{b}_{0}(\mathbf{r})=\operatorname{rot}_{\mathbf{r}}\left[\frac{\mu_{0}}{4\pi}\left(\frac{\mathbf{m}\times \mathbf{r}}
{|\mathbf{r}|}\right)\right]=
\frac{{\mu}_{0}}{4\pi}\operatorname{grad}_{\mathbf{r}}\left[\frac{\langle\mathbf{m}|\mathbf{r}\rangle}{{|\mathbf{r}|}^{3}}\right]+
\frac{{\mu}_{0}}{4\pi}\operatorname{div}_{\mathbf{r}}\left[\frac{\mathbf{r}}{|\mathbf{r}|}\right]\mathbf{m}
</math>  （2-2-2）
:であることが判る。ここで、<math><\ |\ ></math>は、内積を表す。<!-- &lt;, &gt; ではなく \langle, \rangle を使うべきである --><br />
(2)従って、"(1)"を平行移動すれば、'''位置<math>\mathbf{s}\in\mathbb{R}^{3}</math>に置かれた'''磁気モーメント'''m'''が、空間上の位置'''r'''
に作り出す磁束密度が、

::<math>\mathbf{b}_{s}(\mathbf{r})=\operatorname{rot}_{\mathbf{r}}\left[\frac{\mu_{0}}{4\pi}\left(\frac{\mathbf{m}\times (\mathbf{r}-\mathbf{s})}
{|\mathbf{r}-\mathbf{s}|}\right)\right]=
\frac{{\mu}_{0}}{4\pi}\ \operatorname{grad}_{\mathbf{r}}\left[\frac{\langle\mathbf{m}|(\mathbf{r}-\mathbf{s})\rangle}{{|\mathbf{r}-\mathbf{s}|}^{3}}\right]+
\frac{{\mu}_{0}}{4\pi}\ \operatorname{div}_{\mathbf{r}}\left[\frac{\mathbf{r}-\mathbf{s}}{|\mathbf{r}-\mathbf{s}|}\right]\mathbf{m}
</math>  （2-2-3）
:であることが判る。<br />
</BLOCKQUOTE>

従って、上記の磁気モーメント<math>\mathbf{M}(\mathbf{s})d\mathbf{s}</math>それぞれは、磁束密度

:<math>\frac{{\mu}_{0}}{4\pi}\operatorname{grad}_{\mathbf{r}}\left[\frac{\langle\mathbf{M}(\mathbf{s}){d}^{3}\mathbf{s} |(\mathbf{r}-\mathbf{s})\rangle}{{|\mathbf{r}-\mathbf{s}|}^{3}}\right]+
\frac{{\mu}_{0}}{4\pi}\operatorname{div}_{\mathbf{r}}\left[\frac{\mathbf{r}-\mathbf{s}}{|\mathbf{r}-\mathbf{s}|}\right]\mathbf{M}(\mathbf{s}){d}^{3}\mathbf{s}
</math>  （2-2-4）

を作り出す。

上記の磁気ベクトルポテンシャルそれぞれを、全ての<math>\mathbf{s}\in\Omega</math>に渡って足し合わせると、

:<math>\mathbf{B}_{M}(\mathbf{r})
=\frac{\mu_{0}}{4\pi}
{\int}_{\mathbf{s}\in\Omega}
\operatorname{grad}_{\mathbf{r}}\left[\frac{\langle\mathbf{M}(\mathbf{s})|(\mathbf{r}-\mathbf{s})\rangle}{{|\mathbf{r}-\mathbf{s}|}^{3}}\right]+
\operatorname{div}_{\mathbf{r}}\left[\frac{\mathbf{r}-\mathbf{s}}{|\mathbf{r}-\mathbf{s}|}\right]\mathbf{M}(\mathbf{s})
\ {d}^{3}\mathbf{s}
</math>  （2-2-5）

を得る。即ち、物質の磁化<math>\mathbf{M}(\mathbf{r})</math>
は、上記の磁気ベクトルポテンシャル
<math>\mathbf{B}_{M}</math>を空間内に作り出す
<ref name="mizoguchi" />
<ref group="注釈">
尚、磁化が<math>\Omega</math>の外では0であるからといって、<math>\Omega</math>の外で<math>\mathbf{B}_{M}</math>が0となるとは限らない（大概の場合は<math>\Omega</math>の外でも<math>\mathbf{B}_{M}</math>は、0ではない）ことに注意されたい。</ref>
<ref group="注釈">
なお、磁化は、<math>\Omega</math>の外では0なので、上式は、
:<math>\mathbf{B}_{M}(\mathbf{r})
=\frac{\mu_{0}}{4\pi}
\operatorname{grad}_{\mathbf{r}}\left[\frac{\langle\mathbf{M}(\mathbf{s})|(\mathbf{r}-\mathbf{s})\rangle}{{|\mathbf{r}-\mathbf{s}|}^{3}}\right]+
\operatorname{div}_{\mathbf{r}}\left[\frac{\mathbf{r}-\mathbf{s}}{|\mathbf{r}-\mathbf{s}|}\right]\mathbf{M}(\mathbf{s})
\ {d}^{3}\mathbf{s}
</math>
としてもよい（紛らわしいのであまり推奨しない）。
</ref>

さて、ベクトル解析の公式から、
:<math>\mathbf{M}(\mathbf{s})\ =\ {\int}_{\mathbf{s}\in\Omega}\ \operatorname{div}_{\mathbf{r}}\left[\frac{\mathbf{r}-\mathbf{s}}{|\mathbf{r}-\mathbf{s}|}\right]\mathbf{M}(\mathbf{s})
\ {d}^{3}s</math>  （2-2-6）
が判る。<ref group="注釈">
なお、磁化は、<math>\Omega</math>の外では0なので、上式は、
:<math>\mathbf{M}(\mathbf{s})\ =\ {\int}_{\mathbf{s}\in\mathbf{R}^{3}}\ \operatorname{div}_{\mathbf{r}}\left[\frac{\mathbf{r}-\mathbf{s}}{|\mathbf{r}-\mathbf{s}|}\right]\mathbf{M}(\mathbf{s})\ {d}^{3}\mathbf{s}</math>
でもある。</ref>

従って、上記の磁気モーメント<math>\mathbf{M}(\mathbf{s})d\mathbf{s}</math>それぞれが作り出す磁束密度は、

:<math>\frac{{\mu}_{0}}{4\pi}\operatorname{grad}_{\mathbf{r}}\left[\frac{\langle\mathbf{M}(\mathbf{s}){d}^{3}\mathbf{s} |(\mathbf{r}-\mathbf{s})\rangle}{{|\mathbf{r}-\mathbf{s}|}^{3}}\right]
+{\mu}_{0}\delta(\mathbf{r}-\mathbf{s})\mathbf{M}(\mathbf{s}){d}^{3}\mathbf{s}
</math>  （2-2-7）

と書けることが判る。従って、

:<math>\mathbf{B}_{M}(\mathbf{r})
=\frac{\mu_{0}}{4\pi}
{\int}_{\mathbf{s}\in\Omega}
\operatorname{grad}_{\mathbf{r}}\left[\frac{\langle\mathbf{M}(\mathbf{s})|(\mathbf{r}-\mathbf{s})\rangle}{{|\mathbf{r}-\mathbf{s}|}^{3}}\right]\ {d}^{3}{s}
+\mu_{0}\mathbf{M}(\mathbf{r})</math>  （2-2-8）

を得る。<ref group="注釈">
なお、磁化は、<math>\Omega</math>の外では0なので、上式は、
:<math>\mathbf{B}_{M}(\mathbf{r})
=\frac{\mu_{0}}{4\pi}
{\int}_{\mathbf{s}\in\mathbf{R}^{3}}
\operatorname{grad}_{\mathbf{r}}\left[\frac{\langle\mathbf{M}(\mathbf{s})|(\mathbf{r}-\mathbf{s})\rangle}{{|\mathbf{r}-\mathbf{s}|}^{3}}\right]\ {d}^{3}{s}
+\mu_{0}\mathbf{M}(\mathbf{r})</math>
としてもよい（紛らわしいのであまり推奨しない）。
</ref>

今、新たな場<math>\mathbf{H}_{M}</math>を、
:<math>\mathbf{H}_{M}(\mathbf{r})
:=\frac{1}{4\pi}{\int}_{\mathbf{s}\in\Omega}
\operatorname{grad}_{\mathbf{r}}\left[\frac{\langle\mathbf{M}(\mathbf{s})|(\mathbf{r}-\mathbf{s})\rangle}{{|\mathbf{r}-\mathbf{s}|}^{3}}\right]\ {d}^{3}{s}
</math>  （2-2-9）

と定義すると
<ref group="注釈">
なお、磁化は、<math>\Omega</math>の外では0なので、上式は、
:<math>\mathbf{H}_{M}(\mathbf{r})
=\frac{1}{4\pi}
{\int}_{\mathbf{s}\in\Omega}
\operatorname{grad}_{\mathbf{r}}\left[\frac{\langle\mathbf{M}(\mathbf{s})|(\mathbf{r}-\mathbf{s})\rangle}{{|\mathbf{r}-\mathbf{s}|}^{3}}\right]\ {d}^{3}{s}
</math>
としてもよい（紛らわしいのであまり推奨しない）。
</ref>、
:<math>\mathbf{B}_{M}(\mathbf{r})={\mu}_{0}(\mathbf{H}_{M}(\mathbf{r})+\mathbf{M}(\mathbf{r}))</math>  （2-2-10）

が得られる。

さて、以上の議論から「'''物質の磁化Mが既知である場合に<u>限れば</u>、<u>その磁化Mが作り出す磁束密度B<sub>M</sub></u>を計算する術が得られた'''」
ことになる。然しながら、物質の磁化Mが既知でない場合には、上述の関係式のみからは、'''B<sub>M</sub>'''も'''M'''も、判らない。上述の関係式は、
一つの拘束条件を与えているに過ぎないのである。（だが、極めて重要な関係式ではある）。

例えば、磁化'''M<sub>ini</sub>'''を帯びた鉄心が空間におかれていたとき、外部からの磁束密度<math>\mathbf{B}_{\text{ext}}(\mathbf{r})</math>
が与えられたとき、鉄心の磁化は、元々の磁化'''M<sub>ini</sub>'''と、外部からの磁束密度<math>\mathbf{B}_{\text{ext}}(\mathbf{r})</math>
の影響で、元々の磁化'''M<sub>ini</sub>'''とは異なる新たな磁化'''M<sub>con</sub>'''を得ることになる。仮に、この'''M<sub>con</sub>'''
が判れば、全系の磁束密度が計算できるのだが、

:（難所）元々磁化'''M<sub>ini</sub>'''を帯びている物質に、外部から磁束密度<math>\mathbf{B}_{\text{ext}}(\mathbf{r})</math>を印加したとき、最終的に、物質がどのような'''M<sub>con</sub>'''を得るか？

が、実のところは難しい。そこで、一般には。B-H曲線等の実測結果と、上記の拘束条件を考え合わせ、数値計算によって磁束密度や磁化が計算されるのである。但し、線形物質に関して言えば、上の難所は比較的簡単である。このような特殊な物質に関する問題については、次章で述べることにする。

=== 磁荷密度の導入 ===
空間内の領域<math>\Omega</math>に物質が置かれ、前記物質が、定常な（つまり時刻tに依存しない）[[磁化]]<math>\mathbf{M}(\mathbf{r})</math>を帯びているとする。

先に定義した、（磁化による）磁場
:<math>\mathbf{H}_{M}(\mathbf{r})
=\frac{1}{4\pi}
{\int}_{\mathbf{s}\in\Omega}
\operatorname{grad}_{\mathbf{r}}\left[
\frac{\langle\mathbf{M}(\mathbf{s})|(\mathbf{r}-\mathbf{s})\rangle}{{|\mathbf{r}-\mathbf{s}|}^{3}}
\right]
\ {d}^{3}{s}
</math>  （2-3-1）
の原因が、[[磁荷]]であると考えた場合に、それは、どのようなものであるのかを検討しよう。

ベクトル解析の公式より、
:<math>
\begin{align}
\operatorname{div}_{s}\left[\frac{\mathbf{M}(\mathbf{s})}{|\mathbf{r}-\mathbf{s}|} \right]
&= \left\langle
\operatorname{grad}_{\mathbf{s}}
\left[
\frac{1}{|\mathbf{r}-\mathbf{s}|}
\right]
\left|\right.
\mathbf{M}(\mathbf{s})
\right\rangle+
\frac{
\operatorname{div}_{\mathbf{s}}[\mathbf{M}]}
{|\mathbf{r}-\mathbf{s}|} \\
&= \left\langle
\frac{\mathbf{r}-\mathbf{s}}{|\mathbf{r}-\mathbf{s}|^{3}}
\left|\right.
\mathbf{M}(\mathbf{s})
\right\rangle+
\frac{
\operatorname{div}_{\mathbf{s}}[\mathbf{M}]}
{|\mathbf{r}-\mathbf{s}|} \\
&= \frac{\langle\mathbf{M}(\mathbf{s})|(\mathbf{r}-\mathbf{s})\rangle}{{|\mathbf{r}-\mathbf{s}|}^{3}}+
\frac{
\operatorname{div}_{\mathbf{s}}[\mathbf{M}]}
{|\mathbf{r}-\mathbf{s}|}
\qquad \text{(2-3-2)}
\end{align}
</math>

従って、
:<math>
\frac{\langle\mathbf{M}(\mathbf{s})|(\mathbf{r}-\mathbf{s})\rangle}{{|\mathbf{r}-\mathbf{s}|}^{3}}
=\operatorname{div}_{s}\left[\frac{\mathbf{M}(\mathbf{s})}{|\mathbf{r}-\mathbf{s}|} \right]
 - \frac{
\operatorname{div}_{\mathbf{s}}[\mathbf{M}]}
{|\mathbf{r}-\mathbf{s}|}
</math>  （2-3-3）

今、スカラー値関数<math>{\phi}_{M}</math>を、
:<math>{\phi}_{M}(\mathbf{r}):=
-{\int}_{\mathbf{s}\in\Omega}
\frac{\langle\mathbf{M}(\mathbf{s})|(\mathbf{r}-\mathbf{s})\rangle}{{|\mathbf{r}-\mathbf{s}|}^{3}}
\ {d}^{3}\Omega
</math>  （2-3-4）

と定める（磁化が定める磁位、あるいは、磁化が定めるスカラーポテンシャルと言う）と、

:<math>{\phi}_{M}(\mathbf{r})=
-{\int}_{\mathbf{s}\in\Omega}
\operatorname{div}_{s}\left[\frac{\mathbf{M}(\mathbf{s})}{|\mathbf{r}-\mathbf{s}|} \right]
\ {d}^{3}\Omega
+{\int}_{\mathbf{s}\in\Omega}
\frac{
\operatorname{div}_{\mathbf{s}}[\mathbf{M}]}
{|\mathbf{r}-\mathbf{s}|}
\ {d}^{3}\Omega
</math>  （2-3-5）

右辺第一項に、[[発散定理|ガウスの発散定理]]を適用すると、
:<math>
{\int}_{\mathbf{s}\in\Omega}
\operatorname{div}_{s}\left[\frac{\mathbf{M}(\mathbf{s})}{|\mathbf{r}-\mathbf{s}|} \right]
\ {d}^{3}\Omega
={\int}_{\mathbf{s}\in\partial\Omega}
\frac{\mathbf{M}(\mathbf{s})}{|\mathbf{r}-\mathbf{s}|}
\ {d}^{2}(\partial\Omega)
</math>  （2-3-6）

さらに、ベクトル解析の公式を適用すると、
:<math>
{\int}_{\mathbf{s}\in\partial\Omega}
\frac{\mathbf{M}(\mathbf{s})}{|\mathbf{r}-\mathbf{s}|}
\ {d}^{2}(\partial\Omega)
=
{\int}_{\mathbf{s}\in\partial\Omega}
\frac{\langle\mathbf{M}(\mathbf{s})|\mathbf{n}_{\partial\Omega}\rangle}{|\mathbf{r}-\mathbf{s}|}
\ |{d}^{2}(\partial\Omega)|
</math>  （2-3-7）

従って、
:<math>{\phi}_{M}(\mathbf{r})=
-{\int}_{\mathbf{s}\in\Omega}
\operatorname{div}_{s}\left[\frac{\mathbf{M}(\mathbf{s})}{|\mathbf{r}-\mathbf{s}|} \right]
\ {d}^{3}\Omega\ +
{\int}_{\mathbf{s}\in\partial\Omega}
\frac{\langle\mathbf{M}(\mathbf{s})|\mathbf{n}_{\partial\Omega}\rangle}{|\mathbf{r}-\mathbf{s}|}
\ |{d}^{2}(\partial\Omega)|
</math>  （2-3-8）
ここで、<math>\partial\Omega</math>は、領域<math>\Omega</math>の境界面を意味する。
また、<math>\mathbf{n}_{\partial\Omega}</math>は、<math>\partial\Omega</math>の法線ベクトルを意味する。

今、体積磁荷密度と、表面磁荷密度を、
*<math>{\rho}_{M}:=-\operatorname{div}_{\mathbf{r}}[\mathbf{M}]</math> （体積磁荷密度）  （2-3-10）
*<math>{\sigma}_{M}:=<\mathbf{M}\ |\ \mathbf{n}_{\partial\Omega}></math> （表面磁荷密度）  （2-3-11）
により定めると、

:<math>{\phi}_{M}(\mathbf{r})=
{\int}_{\mathbf{s}\in\Omega}
\frac{{\rho}_{M}(\mathbf{s})}{|\mathbf{r}-\mathbf{s}|}
\ {d}^{3}\Omega\ +
{\int}_{\mathbf{s}\in\partial\Omega}
\frac{
{\sigma}_{M}(\mathbf{s})
}{|\mathbf{r}-\mathbf{s}|}
\ |{d}^{2}(\partial\Omega)|
</math>  （2-3-12）

を得る。

一方、<math>{\phi}_{M}(\mathbf{r})</math>の定義により、
:<math>{\mathbf{H}}_{M}(\mathbf{r})=-\operatorname{grad}_{\mathbf{r}}[{\phi}_{M}(\mathbf{r})]</math>  （2-3-13）

であるため、
:<math>{\mathbf{H}}_{M}(\mathbf{r})=
{\int}_{\mathbf{s}\in\Omega}
\operatorname{grad}_{\mathbf{r}}
\left[
\frac{{-\rho}_{M}(\mathbf{s})}{|\mathbf{r}-\mathbf{s}|}
\right]
\ {d}^{3}\Omega\ +
{\int}_{\mathbf{s}\in\partial\Omega}
\operatorname{grad}_{\mathbf{r}}
\left[
\frac{
{-\sigma}_{M}(\mathbf{s})
}{|\mathbf{r}-\mathbf{s}|}
\right]
\ |{d}^{2}(\partial\Omega)|
</math>  （2-3-14）
である。

=== 磁化電流の導入 ===
時間的に定常な磁化が作り出す磁気ベクトルポテンシャルを、別の側面から考察してみることにしよう。ここでは、「磁束密度の原因は、電流に帰される」という思想に従い、
だとすれば、「磁化と等価な効果を発揮する電流」がどのようなものかを検討することにする。

前記<math>\mathbf{A}_{M}(\mathbf{r})</math>に、
[[ベクトル解析]]の公式を適用すると、

:<math>\begin{align}
\mathbf{A}_{M}(\mathbf{r})
&= \frac{\mu_{0}}{4\pi}
{\int}_{\mathbf{s}\in\Omega}
\left(\frac{(\mathbf{M}(\mathbf{s}))\times (\mathbf{r}-\mathbf{s})}
{|\mathbf{r}-\mathbf{s}|}
\right)
\ {d}^{3}\mathbf{s} \\
&= \frac{\mu_{0}}{4\pi}
{\int}_{\mathbf{s}\in\Omega}
\left(\frac{\operatorname{rot}_{\mathbf{s}}[\mathbf{M}(\mathbf{s})]}
{|\mathbf{r}-\mathbf{s}|}
\right)
\ {d}^{3}\mathbf{s}
+ \frac{\mu_{0}}{4\pi}
{\int}_{\mathbf{s}\in\partial\Omega}
\left(
\frac{\mathbf{M}(\mathbf{s})\times \mathbf{n}_{\partial\Omega}}{|\mathbf{r}-\mathbf{s}|}
\right)
\ |{d}^{2}(\partial\Omega)|
\qquad \text{(2-4-1)}
\end{align}
</math>

となることが判る。ここで、<math>\operatorname{rot}_{\mathbf{s}}</math>
は、変数<math>\mathbf{s}</math>についての回転微分を意味する。<math>\partial\Omega</math>は、領域<math>\Omega</math>の境界を意味する。
上式の右辺第二項の積分において、「<math>\partial\Omega</math>の面素の絶対値」を意味し
所謂普通の面積分ではない（この「面積分」は、「ベクトル場からベクトルを作る」特殊なもの）ので注意が必要である。
この積分において、<math>\mathbf{n}_{\partial\Omega}</math>は、<math>\partial\Omega</math>
の法線ベクトルを意味する。
<blockquote>
実際、スカラー倍の回転微分の公式より、

:<math>\begin{align}
\operatorname{rot}_{\mathbf{s}}\left[\frac{\mathbf{M}(\mathbf{s})}{|\mathbf{r}-\mathbf{s}|} \right]
&= \left(\operatorname{grad}_{\mathbf{s}}\left[\frac{\mathbf{M}(\mathbf{s})}{|\mathbf{r}-\mathbf{s}|} \right]\right)\times\mathbf{M}(\mathbf{s})
+ \frac{\operatorname{rot}_{\mathbf{s}}\left[ \mathbf{M}(\mathbf{s})\right]}{|\mathbf{r}-\mathbf{s}|} \\
&= \frac{(\mathbf{r}-\mathbf{s})\times \mathbf{M}(\mathbf{s})}{{|\mathbf{r}-\mathbf{s}|}^{3}}
+\frac{\operatorname{rot}_{\mathbf{s}}\left[ \mathbf{M}(\mathbf{s})\right]}{|\mathbf{r}-\mathbf{s}|}
\qquad \text{(2-4-2)}
\end{align}
</math>
従って、
:<math>
\frac{\mathbf{M}(\mathbf{s})\times (\mathbf{r}-\mathbf{s}) }{{|\mathbf{r}-\mathbf{s}|}^{3}}
=\frac{\operatorname{rot}_{\mathbf{s}}\left[ \mathbf{M}(\mathbf{s})\right]}{|\mathbf{r}-\mathbf{s}|}
-\operatorname{rot}_{\mathbf{s}}\left[\frac{\mathbf{M}(\mathbf{s})}{|\mathbf{r}-\mathbf{s}|} \right]
</math>  （2-4-3）

従って、
:<math>\mathbf{A}_{M}(\mathbf{r})
=\frac{\mu_{0}}{4\pi}
{\int}_{\mathbf{s}\in\Omega}
\frac{\operatorname{rot}_{\mathbf{s}}\left[ \mathbf{M}(\mathbf{s})\right]}{|\mathbf{r}-\mathbf{s}|}\ {d}^{3}\mathbf{s}\ 
-\frac{\mu_{0}}{4\pi}
{\int}_{\mathbf{s}\in\Omega}
\operatorname{rot}_{\mathbf{s}}\left[\frac{\mathbf{M}(\mathbf{s})}{|\mathbf{r}-\mathbf{s}|} \right] 
\ {d}^{3}\mathbf{s}\ 
</math>　　（2-4-4）

を得る。さらに、上式の右辺第二項の積分にベクトル解析の公式を適用すると、
:<math>-\frac{\mu_{0}}{4\pi}
{\int}_{\mathbf{s}\in\Omega}
\operatorname{rot}_{\mathbf{s}}\left[\frac{\mathbf{M}(\mathbf{s})}{|\mathbf{r}-\mathbf{s}|} \right] 
\ {d}^{3}\mathbf{s}\ 
=\frac{\mu_{0}}{4\pi}
{\int}_{\mathbf{s}\in\partial\Omega}
\left(\frac{\mathbf{M}(\mathbf{s})}{|\mathbf{r}-\mathbf{s}|} \right) 
\times {d}^{2}\mathbf{s}\ 
</math>　　（2-4-5）
を得る。さらに、右辺にベクトル解析の公式を適用すると、
:<math>
\frac{\mu_{0}}{4\pi}
{\int}_{\mathbf{s}\in\partial\Omega}
\left(\frac{\mathbf{M}(\mathbf{s})}{|\mathbf{r}-\mathbf{s}|} \right) \times {d}^{2}\mathbf{s}\ 
=
\frac{\mu_{0}}{4\pi}
{\int}_{\mathbf{s}\in\partial\Omega}
\left(
\frac{\mathbf{M}(\mathbf{s})\times \mathbf{n}_{\partial\Omega}}{|\mathbf{r}-\mathbf{s}|} 
\right)
\ |{d}^{2}(\partial\Omega)|\ 
</math>  （2-4-6）

が得られる。以上から示すべき式が証明された。
</blockquote>

今、<math>\boldsymbol{i}_{M}(\mathbf{r}),\ {\mathbf{K}}_{M}</math>を、
*<math>\boldsymbol{i}_{M}(\mathbf{r}) :=
\operatorname{rot}[\mathbf{M}(\mathbf{r})]
</math> （体積磁化電流密度）  （2-4-7）
*<math>{\mathbf{K}}_{M}(\mathbf{r}) :=
\mathbf{M}(\mathbf{r})\times \mathbf{n}_{\partial\Omega}(\mathbf{r})</math> (表面磁化電流密度)  （2-4-8）

と置くと、結局、
:<math>\mathbf{A}_{M}(\mathbf{r})
=\frac{\mu_{0}}{4\pi}
{\int}_{\mathbf{s}\in\Omega}
\left(\frac{\boldsymbol{i}_{M}(\mathbf{s})}{|\mathbf{r}-\mathbf{s}|} 
\right)
\ {d}^{3}\mathbf{s}
+\ \frac{\mu_{0}}{4\pi}
{\int}_{\mathbf{s}\in\partial\Omega} 
\left(
\frac{\mathbf{K}_{M}(\mathbf{s})}{|\mathbf{r}-\mathbf{s}|} 
\right)\ |{d}^{2}(\partial\Omega)|\ 
</math>  （2-4-9）

となる。

両辺の回転微分を取ると、
:<math>\mathbf{B}_{M}(\mathbf{r})
=\frac{\mu_{0}}{4 \pi }
\int_{\mathbf{s}\in\mathbb{R}^{3}} 
\left(\frac{
\boldsymbol{i}_{M}(\mathbf{s}) \times (\mathbf{r} - \mathbf{s})
}{|\mathbf{r} - \mathbf{s}|^3}\right) d^3 
\mathbf{s}\ +\ 
\frac{\mu_{0}}{4\pi}
{\int}_{\mathbf{s}\in\partial\Omega} 
\left(\frac{\mathbf{K}_{M}(\mathbf{s}) \times (\mathbf{r} - \mathbf{s})}{|\mathbf{r} - \mathbf{s}|^3}\right) 
\ |{d}^{2}(\partial\Omega)|</math>  （2-4-10）

が、判る。この見方は、特に、磁化が一様な場合（より一般にはrot[M]=0）といった特殊な場合に、特に威力を発揮する。


== 強制電流と磁化の両方が既知のとき(一般論) ==
=== 全系の磁束密度 ===
強制電流i<sub>fc</sub> (forced current)と磁化Mの両方が既知とし、これら以外に磁束密度の原因となるものがないとした場合、
全系の磁束密度'''B'''<sub>tot</sub>
は、強制電流i<sub>fc</sub>に起因する磁束密度の成分（式(1-2)から求められる）
'''B'''<sub>fc</sub>と、磁化Mが作り出す磁束密度の成分（式(2-2-8)から求められる）
'''B'''<sub>M</sub>によって、

:'''B'''<sub>tot</sub>='''B'''<sub>fc</sub>+'''B'''<sub>M</sub> (3-1-1)

と表される。
<blockquote>
ここで、よく注意しておかないといけないことは、「強制電流i<sub>fc</sub>と磁化Mの両方が<u>既知</u>」という言葉の意味であるが、
仮に、強制電流i<sub>fc</sub>（や、他の磁化）がない状態での磁化'''M<sub>ini</sub>'''が既知したとして、
(磁化'''M<sub>ini</sub>'''が作り出す磁束密度'''B'''<sub>Mini</sub> がどうなるかは(2-2-8)から計算できるが)

:'''B'''<sub>tot</sub>='''B'''<sub>fc</sub>+'''B'''<sub>Mini</sub> (間違った式)

は、よほど特殊な場合を除き成り立たない。要は間違いである。
</blockquote>
式(3-1-1)より'''正確に書くならば、「外場の影響等により磁化が変化した後の磁化'''
'''M<sub>con</sub>'''
'''が作り出す磁束密度'''」
'''B'''<sub>Mcon</sub>を用いて、

:'''B'''<sub>tot</sub>='''B'''<sub>fc</sub>+'''B'''<sub>Mcon</sub> (正しい式)

である。結局以下の難所は残ったままである。

:（難所）元々磁化'''M<sub>ini</sub>'''を帯びている物質に、外部から磁束密度<math>\mathbf{B}_{\text{ext}}(\mathbf{r})</math>を印加したとき、最終的に、物質がどのような'''M<sub>con</sub>'''を得るか？

この問題が、実のところは難しい。つまり、強制電流が作った磁束密度や、磁化自身が作り出す磁束密度により、物質の磁化が最初の磁化から変化してしまうという問題がある。
一般には。B-H曲線等の実測結果と、上記の拘束条件を考え合わせ、数値計算によって磁束密度や磁化が計算されるのである。但し、線形物質に関して言えば、上の難所は比較的簡単である。このような特殊な物質に関する問題については、次章で述べることにする。

再び式(3-1-1)について考えよう。ベクトル解析の公式から
:<math>\operatorname{div}[\mathbf{B}_{fc}]=0</math>   (3-1-2a)
:<math>\operatorname{div}[\mathbf{B}_{M}]=0</math>   (3-1-2b)
:<math>\operatorname{div}[\mathbf{B}_{\text{tot}}]=0</math>   (3-1-2c)

であることが判る。即ち、磁束保存の式が満たされることが判る。

=== B=μ<sub>0</sub>(H+M)について ===
前節同様に、強制電流i<sub>fc</sub>と磁化Mの両方が既知とし、これら以外に磁束密度の原因となるものがないとした場合
について考える。式(3-1-1)のB<sub>tot</sub>に対し、新たな場H<sub>tot</sub>を、

:<math>\mathbf{H}_{\text{tot}}(\mathbf{r}):=\frac{\mathbf{B}_{\text{tot}}(\mathbf{r})
}{{\mu}_{0}}-\mathbf{M}(\mathbf{r})</math> （3-2-1）

と定める。式(3-2-1)は、所謂”B=μ<sub>0</sub>(H+M)”に他ならない。


さらに、強制電流<math>\boldsymbol{i}_{fc}</math>が作る磁束密度<math>\mathbf{B}_{fc}(\mathbf{r})</math>は、式(1-2)に
より、（'''周辺の物質の有無にかかわらず'''）定まるが、式(1-3)に倣い
「強制電流<math>\boldsymbol{i}_{fc}</math>が作る磁場<math>\mathbf{H}_{fc}(\mathbf{r})</math>」を、
:<math>\mathbf{H}_{fc}(\mathbf{r}):=\frac{1}{{\mu}_{0}}\mathbf{B}_{fc}(\mathbf{r})</math> （3-2-2）
によって定め、磁化による磁場H<sub>M</sub>を式（2-2-10）のように定めると、
式（3-2-1）と式（3-2-2）、式（2-2-10）より、

:<math>\begin{align}
\mathbf{H}_{\text{tot}}(\mathbf{r}) &=
\frac{\mathbf{B}_{\text{tot}}(\mathbf{r})}{{\mu}_{0}}-\mathbf{M}(\mathbf{r}) \\
&=
\frac{
\mathbf{B}_{fc}(\mathbf{r})+\mathbf{B}_{M}(\mathbf{r})
}{{\mu}_{0}}-\mathbf{M}(\mathbf{r}) \\
&= \mathbf{H}_{fc}(\mathbf{r})+(\mathbf{H}_{M}(\mathbf{r})+\mathbf{M}(\mathbf{r}))-\mathbf{M}(\mathbf{r}) \\
&= \mathbf{H}_{fc}(\mathbf{r})+\mathbf{H}_{M}(\mathbf{r})
\qquad \text{(3-2-3)}
\end{align}</math>
が得られる。従って、
:<math>\operatorname{rot}[\mathbf{H}_{\text{tot}}(\mathbf{r})]
=\operatorname{rot}[\mathbf{H}_{fc}(\mathbf{r})]
+\operatorname{rot}[\mathbf{H}_{M}(\mathbf{r})]
=\boldsymbol{i}_{fc}</math>  （3-2-4）
が得られる。これは、即ち（静磁場の）アンペールの法則である。
<blockquote>
式（3-2-4）を示そう。式（2-2-9） 及び以下の式(3-2-5)より、
:rot grad=0 (3-2-5)

<math>\mathbf{H}_{M}(\mathbf{r})</math>に対し、回転微分を作用させると、
:<math>\operatorname{rot}[\mathbf{H}_{M}(\mathbf{r})]
=\frac{1}{4\pi}
{\int}_{\mathbf{s}\in\Omega}
\operatorname{rot}_{\mathbf{r}}\left[
\operatorname{grad}_{\mathbf{r}}\left[
\frac{\langle\mathbf{M}(\mathbf{s})|(\mathbf{r}-\mathbf{s})\rangle}{{|\mathbf{r}-\mathbf{s}|}^{3}}
\right]
\right]
\ {d}^{3}{s}\ =0
</math> (3-2-6)
が判る。従って、磁化に起因する磁場<math>\mathbf{H}_{M}(\mathbf{r})</math>は、
:<math>\operatorname{rot}[\mathbf{H}_{M}(\mathbf{r})]=0</math> (3-2-7)
を満たす。


</blockquote>
式(3-2-7)と式(1-5)より、式（3-2-4）が示された。

=== 透磁率の導入 ===
透磁率なる概念を、新たに導入する。
即ち、3次正方行列に値を取る行列値関数<math>{\mu}(\mathbf{r})</math>（透磁率）を用いて、

:<math>\mathbf{B}_{\text{tot}}(\mathbf{r})={\mu}(\mathbf{r})\mathbf{H}_{\text{tot}}(\mathbf{r})</math> (3-3-1)

と書けるものとする。この透磁率は、磁化の概念の”すり替え”に過ぎない概念である。
尚、通常は、透磁率は、スカラー値関数と考えてよい場合が多いのだが、その場合は、
「対角成分が全部同じ値で、それ以外の成分が0の行列値関数と、スカラー値関数が同一視できる」ことを思い起こせばいい。

さらに、式(3-3-1) で見たような、透磁率（透磁率テンソル）を用いて、
:<math>\mathbf{B}_{\text{tot}}(\mathbf{r})={\mu}(\mathbf{r})\mathbf{H}_{\text{tot}}(\mathbf{r})</math>   (3-3-2)
と書き表せ、かつ、μが全点で正則行列（逆行列を持つ）とする。(そう考えても、”あまり”一般性を失わない）逆に言えば、
:<math>\mathbf{H}_{\text{tot}}(\mathbf{r})={\nu}(\mathbf{r})\mathbf{B}_{\text{tot}}(\mathbf{r})</math>   (3-3-3)
である。ここで、

:ν:=μ<sup>-1</sup>  (3-3-4)

は、各点でμの逆行列を与えるような行列値関数である。即ち、

:μ('''s''') ν('''s''')=μ('''s''') μ<sup>-1</sup>('''s''')=E<sub>3</sub>   (3-3-5)

を充たすような行列値関数である。（ここでE<sub>3</sub>は3次の単位行列)である。


これまでの議論では、電流素片や磁気双極子モーメントが作り出した磁束密度/磁場について論じてきたが、
式(3-1-2),式（3-2-4）より、このようにして作られた全系の磁場が、磁束保存の式と、アンペールの法則を
充たすことが判った。
ここからは、逆に、磁束保存の式(3-3-6)と、アンペールの法則(3-3-7)、即ち、

:<math>\operatorname{div}[\mathbf{B}_{fc}]=0</math>   (3-3-6a)
:<math>\operatorname{div}[\mathbf{B}_{M}]=0</math>   (3-3-6b)
:<math>\operatorname{div}[\mathbf{B}_{\text{tot}}]=0</math>   (3-3-6c)
:<math>\operatorname{rot}[\mathbf{H}_{\text{tot}}]=\boldsymbol{i}_{fc}</math>    (3-3-7)

を出発点とし、全系の磁気ベクトルポテンシャル<math>\mathbf{A}_{\text{tot}}</math>が充たす微分方程式を導出する。

まず、磁束保存の式  (3-3-6c)は、強制電流に起因する成分 (3-3-6a)、磁化に起因する成分 (3-3-6b)
それぞれについて成り立つため、（ポアンカレの補助定理より）それぞれがベクトルポテンシャルを持ち、

:<math>\operatorname{rot}[\mathbf{A}_{fc}]=\mathbf{B}_{fc}</math>   (3-3-8a)
:<math>\operatorname{rot}[\mathbf{A}_{M}]=\mathbf{B}_{M}</math>   (3-3-8b)

をみたすようなベクトル場<math>\mathbf{A}_{fc}</math>と、<math>\mathbf{A}_{M}</math>が存在する。
このようなベクトル場<math>\mathbf{A}_{fc}</math>と、<math>\mathbf{A}_{M}</math>は、ゲージ不定性を除き一意に定まるが、
本記事では、強制電流に起因する成分<math>\mathbf{A}_{fc}</math>は、式(1-1)のものを採用し、
と、磁化に起因する成分<math>\mathbf{A}_{M}</math>式(2-1-5）を採用することにする。

全系のベクトルポテンシャルを、
:<math>\mathbf{A}_{\text{tot}}:=\mathbf{A}_{M}+\mathbf{A}_{fc}</math>   (3-3-9)
と定めると、
:<math>\operatorname{rot}[\mathbf{A}_{\text{tot}}]=\mathbf{B}_{\text{tot}}</math>   (3-3-8c)
が成り立つ。これは、(3-3-6c)の、磁束保存の式をベクトルポテンシャルを用いて書いたものに他ならない。

さらに、式(3-3-8c)に、式 (3-3-2)と、アンペールの法則(3-3-7)を考え併せると、
:<math>
\operatorname{rot}[{\nu}\operatorname{rot}[\mathbf{A}]]
- \boldsymbol{i}_{fc}=0
</math>  (3-3-10)

が得られ、未知のベクトル場BやHが消え、未知のベクトル場はAのみとなる。これを
「ベクトルポテンシャルによる静磁場の方程式」（静磁場の支配方程式）という。

=== 静磁場の境界条件 ===
物質の境界において、式(3-3-6a)、式 (3-3-7)に、ガウスの発散定理や、[[ケルビンストークスの定理]]を適用すると、
磁場および磁束密度の境界条件が得られる。


=== 静磁場のエネルギーと停留値関数 ===
前節の前提条件において、
全系の磁気ベクトルポテンシャル<math>\mathbf{A}_{\text{tot}}</math>は、以下の汎関数Fの停留関数となることが判る。<ref name="sandf" /><ref name="akita2" />

:<math>\begin{align}
F[\mathbf{A}] &:= {\int}_{\mathbf{s}\in\mathbb{R}^{3}}
\frac{1}{2}
<\ \mathbf{H}(\mathbf{s})
\ |\ \mathbf{B}(\mathbf{s})>
-
<\ \boldsymbol{i}_{fc}(\mathbf{s})|\ \mathbf{A}_{\text{tot}}(\mathbf{s})>
\ d\mathbf{s} \\
&= {\int}_{\mathbf{s}\in\mathbb{R}^{3}}\frac{1}{2}
<\ {\nu}(\mathbf{s})(\operatorname{rot}_{\mathbf{s}}[\mathbf{A}(\mathbf{s})])
\ |\ (\operatorname{rot}_{\mathbf{s}}[\mathbf{A}(\mathbf{s})])>
-
<\ \boldsymbol{i}_{fc}(\mathbf{s})|\ \mathbf{A}(\mathbf{s})>
\ d\mathbf{s}
\qquad \text{(3-5-1)}
\end{align}</math>

次元解析をすると、式(3-5-1)は、エネルギーの次元を持つことが判るが、実際に、式(3-5-1)は、全系のエネルギーとなっている。
ここで、<math><\ |\ ></math>は、内積を表す。νの定義は、式(3-3-4)に記載のとおりである。


式(3-5-1)が、本当に、「磁気ベクトルポテンシャルの停留汎関数」ことを検証しよう。
Aに対し、微小な摂動δA<ref name="dd" group="注釈">ここでいうδは、摂動あるいは変分を表し、ディラックのδと紛らわしいが全然別物である。</ref>
を与えた際の第一変分δF<ref name="dd" group="注釈" />、即ち、

:δF=F[A+δA]-F[A]    (3-5-2)

を求める。

F[A+δA]は、以下の被積分関数を、全空間でｓについて積分したものである。
:<math>
\frac{1}{2}
<\ {\nu}(\mathbf{s})(\operatorname{rot}_{\mathbf{s}}[\mathbf{A}(\mathbf{s})+\delta\mathbf{A}(\mathbf{s})])
\ |\ (\operatorname{rot}_{\mathbf{s}}[\mathbf{A}(\mathbf{s})+\delta\mathbf{A}(\mathbf{s})])>
\ -
<\ \boldsymbol{i}_{fc}(\mathbf{s})|\ \mathbf{A}(\mathbf{s})+\delta\mathbf{A}(\mathbf{s})>
</math>    (3-5-3)

一方、F[A]は、以下の被積分関数を、全空間でｓについて積分したものである。
:<math>
\frac{1}{2}
<\ {\nu}(\mathbf{s})(\operatorname{rot}_{\mathbf{s}}[\mathbf{A}(\mathbf{s})])
\ |\ (\operatorname{rot}_{\mathbf{s}}[\mathbf{A}(\mathbf{s})])>
\ -
<\ \boldsymbol{i}_{fc}(\mathbf{s})|\ \mathbf{A}(\mathbf{s})>
</math>    (3-5-4)

したがって、δFは、以下の(3-5-5)を、全空間でｓについて積分したものである。
:<math>
<\ {\nu}(\mathbf{s})(\operatorname{rot}_{\mathbf{s}}[\mathbf{A}(\mathbf{s})])
\ |\ (\operatorname{rot}_{\mathbf{s}}[\delta\mathbf{A}(\mathbf{s})])>+
\frac{1}{2}<\ {\nu}(\mathbf{s})(\operatorname{rot}_{\mathbf{s}}[\delta\mathbf{A}(\mathbf{s})])
\ |\ (\operatorname{rot}_{\mathbf{s}}[\delta\mathbf{A}(\mathbf{s})])>
\ -
<\ \boldsymbol{i}_{fc}(\mathbf{s})|\ \delta\mathbf{A}(\mathbf{s})>
</math>    (3-5-5)

式(3-5-5)から、二次の微小項を無視すると、Fの第”一”変分
:<math>\delta F[A]={\int}_{\mathbf{s}\in\mathbb{R}^{3}}
<\ {\nu}(\mathbf{s})(\operatorname{rot}_{\mathbf{s}}[\mathbf{A}(\mathbf{s})])
\ |\ (\operatorname{rot}_{\mathbf{s}}[\delta\mathbf{A}(\mathbf{s})])>
\ -
<\ \boldsymbol{i}_{fc}(\mathbf{s})|\ \delta\mathbf{A}(\mathbf{s})>
\ {d}^{3}\mathbf{s}
</math>     (3-5-6)

を得る。
式(3-5-6) に、ベクトル解析の公式を適用する。
:<math>
\operatorname{div}[\mathbf{X}\times\mathbf{Y}]
=<\ \mathbf{Y}\ |\ \operatorname{rot}[\mathbf{X}]\ >-
<\ \mathbf{X}\ |\ \operatorname{rot}[\mathbf{Y}]\ >
</math>  (3-5-7)

において
:<math>\begin{align} 
\mathbf{X} &:= {\nu}(\mathbf{s})\operatorname{rot}[\mathbf{A}(\mathbf{s})] & \text{(3-5-8)} \\ 
\mathbf{Y} &:= \delta\mathbf{A} & \text{(3-5-9)} 
\end{align}
</math>

を代入すると、
:<math>\begin{align}
\delta F[\mathbf{A}] &=
{\int}_{\mathbf{s}\in\mathbb{R}^{3}}
-\operatorname{div}_{\mathbf{s}}[({\nu}(\mathbf{s})
\operatorname{rot}_{\mathbf{s}}[\mathbf{A}]
)\times\delta\mathbf{A}]
\ {d}^{3}\mathbf{s} \\
&+
{\int}_{\mathbf{s}\in\mathbb{R}^{3}}
<(\operatorname{rot}_{\mathbf{s}}[({\nu}(\mathbf{s}))
\operatorname{rot}_{\mathbf{s}}[\mathbf{A}(\mathbf{s})]\ ]
- \boldsymbol{i}_{fc}(\mathbf{s})
\ 
|\ 
\delta\mathbf{A}(\mathbf{s})
>
\ {d}^{3}\mathbf{s}
\qquad \text{(3-5-10)}
\end{align}</math>

式(3-5-10)の第一項については、式(3-4-6)と、「真空中ではdiv[H]=0であること」を考え併せると、
結局物質Ω内の効果しか寄与しないことが判る。さらにガウスの発散定理を考慮すると、式(3-5-10)の第一項は、

:<math>
{\int}_{\mathbf{s}\in\Omega}
-\operatorname{div}_{\mathbf{s}}[({\nu}(\mathbf{s}) 
\operatorname{rot}_{\mathbf{s}}[\mathbf{A}]
)\times\delta\mathbf{A}\ ] 
\ {d}^{3}\mathbf{s}
= 
-{\int}_{\mathbf{s}\in\partial\Omega}
({\nu}(\mathbf{s}) 
\operatorname{rot}_{\mathbf{s}}[\mathbf{A}]
)\times\delta\mathbf{A} 
\ {d}^{2}(\partial\Omega)
</math>  (3-5-11)

となる。従って、式(3-5-10)の第一項が、任意の摂動δAに対して0になるためには、ノイマン条件即ち、
:<math>(\operatorname{rot}[\mathbf{A}])\times{n}_{\partial\Omega}\ =\ 0</math>  (3-5-12)

が満たされればよい。さらに、式(3-5-10)の第二項が、任意の摂動δAに対して0になるためには、
:<math>
\operatorname{rot}[{\nu}\operatorname{rot}[\mathbf{A}]]
- \boldsymbol{i}_{fc}=0
</math>  (3-5-12)

であればよい。これは、「ベクトルポテンシャルによる静磁場の方程式」（静磁場の支配方程式）に(式(3-3-10))他ならない。以上から、ノイマン条件下での、静磁場の支配方程式の解は、汎関数Fの停留関数となること判る。

== ベクトル解析に関する補足 ==
静磁場の解析には、静磁場の解析には、やや高度なベクトル解析の知識
(例えば 藤本<ref name="fujimoto">
藤本 淳夫 (著) ベクトル解析 (現代数学レクチャーズ C- 1) 培風館 (1979/01) </ref>参照）
が要求される。同じ式変形が何度も繰り返されるのを避ける目的から、
一般的な知識は既知とする立場から、本記事の内容理解に必要と思われる事柄に限り、簡単に説明する。
が要求される。同じ式変形が何度も繰り返されるのを避ける目的から、
一般的な知識は既知とする立場から、本記事の内容理解に必要と思われる事柄に限り、簡単に説明する。
=== ベクトル場の代数演算と微分作用素 ===
ベクトル場に代数演算を施したものに、微分作用素を作用させた場合に成り立つ公式について、本記事で用い、かつ、あまり本に載っていないものについて、簡潔にまとめる。

==== (1) <F|∇>について ====
F=(f<sub>1</sub>,f<sub>2</sub>,f<sub>3</sub>)をベクトル場とする。このとき、
:<math>
\left\langle\ \mathbf{F}\ |\ \nabla \ \right\rangle :=
{f}_{1}{\partial \over \partial {x}_{1}}  +  {f}_{2}{\partial \over \partial {x}_{2}} + {f}_{3}{\partial \over \partial {x}_{3}}
</math>
と定義する。ここで、∇は、
:<math>\nabla := \mathbf{\hat{x}} {\partial \over \partial x}  + \mathbf{\hat{y}} {\partial \over \partial y} + \mathbf{\hat{z}} {\partial \over \partial z}</math>
を意味する。<F|∇>のことを、”F・∇”と書くこともある。

Gを、ベクトル場としたとき、<F|∇>をGに作用させると、
:<math>
\left\langle\ \mathbf{F}\ |\ \nabla \ \right\rangle \mathbf{G}=
{f}_{1}{\partial{g}_{1} \over \partial {x}_{1}}  +
{f}_{2}{\partial{g}_{2} \over \partial {x}_{2}} +
{f}_{3}{\partial{g}_{3} \over \partial {x}_{3}}
=(J[\mathbf{G}])\cdot\mathbf{F}
</math>

が成り立つ。ここで、<math>J[\mathbf{G}]</math>は、Gの[[ヤコビ行列]]（ヤコビアンではない）を意味する。

==== (2)rot とスカラー倍、ベクトル積 ====
F,Gを、ベクトル場、fを、スカラー値関数とする。このとき、以下が成り立つ。
:<math>
\operatorname{rot}[f\mathbf{F}]=
\operatorname{grad}[f]\times\mathbf{F}+f\cdot \operatorname{rot}[\mathbf{F}]</math>

:<math>
\operatorname{rot}[\mathbf{F}\times\mathbf{G}]=
\left\langle\ \mathbf{G}\ |\ \nabla \ \right\rangle \mathbf{F}
-\left\langle\ \mathbf{F}\ |\ \nabla \ \right\rangle \mathbf{G}
+\mathbf{F}\operatorname{div}[\mathbf{G}]
-\mathbf{G}\operatorname{div}[\mathbf{F}]
</math>

==== (3)div とスカラー倍、ベクトル積 ====
F,Gを、ベクトル場、fを、スカラー値関数とする。このとき、
:<math>
\operatorname{div}
\left[
f\mathbf{F}
\right]
=\left\langle
\mathbf{F}|\operatorname{grad}[f]
\right\rangle+
f\operatorname{div}[\mathbf{F}]
</math>

=== ディラックのデルタと、重積分 ===
一般のn変数関数（関数はベクトル値関数であってよい)fの定義域をΩとしたとき、
:<math>f(\mathbf{r})={\int}_{\mathbf{s}\in \Omega} f(\mathbf{s}){\delta}^{n}(\mathbf{r}-\mathbf{s})\ {d}^{n}s\ =f*{\delta}^{n}</math>
であることが知られている。ここで、上式の"*"は、[[合成積]]([[乗算]]ではない)である。また、δ<sup>n</sup>は、n変数の[[δ関数]]である。

=== 発散微分とディラックのデルタ ===
'''原点(r=0)を除いて'''、
::<math>\operatorname{div}_{\mathbf{r}}\left[\frac{\mathbf{r}}{|\mathbf{r}|}\right]=\frac{3}{|\mathbf{r}|^{3}}-\frac{3{|\mathbf{r}|}^{2}}{{|\mathbf{r}|}^{5}}=0</math>
:であり、原点(r=0)を中心とする、球体BLに対し、[[ガウスの発散定理]]を用いると、
:::<math>{\int}_{\mathbf{r}\in BL} \operatorname{div}_{\mathbf{r}}\left[\frac{\mathbf{r}}{|\mathbf{r}|}\right]\ d^{3}r =4\pi</math>
:となるので、結局、
::<math>\operatorname{div}_{\mathbf{r}}\left[\frac{\mathbf{r}}{|\mathbf{r}|}\right]=4\pi\delta(\mathbf{r})</math>
であることが判る。(ここで、δは、3変数の[[δ関数]]である。<br />
以上の議論を平行移動させると、
:<math>\operatorname{div}_{\mathbf{r}}\left[\frac{\mathbf{r}-\mathbf{s}}{|\mathbf{r}-\mathbf{s}|}\right]=4\pi\delta(\mathbf{r}-\mathbf{s})</math>
であることが判る。

=== 分数関数の偏微分 ===
:<math>\begin{align}
\frac{\partial }{\partial \mathbf{s}_{i}}
\left[\frac{1}{|\mathbf{r}-\mathbf{s}|}\right]
&=
\frac{\partial }{\partial {s}_{i}}
\left[\frac{1}{{(({r}_{i}-{s}_{i})^{2}+\text{const.})}^{1/2}}\right] \\
&=
\frac{\partial }{\partial {s}_{i}}
\left[{(({r}_{i}-{s}_{i})^{2}+\text{const.})}^{-(1/2)}\right]
\end{align}
</math>

=== 面積分に関する補足 ===
ここでは、様々なベクトル場の面積分についてまとめる。
「絶対値による面積分」「外積面積分」という用語は一般的な用語ではないが、他に
適切な表現がないため、この場限りでそのような言い方をする。
本記事内での定義は、それぞれ以下の通り。

Iを<math>\mathbb{R}^{2}</math>の閉長方形、<br />
<math>\varphi:I\to \mathbb{R}^{3}</math>は、
Iの近傍で、殆ど至る所区分的に滑らかかつ、非退化であり、かつ、Iの内部で[[単射]]なベクトル値関数、<br />
<math>S:=\varphi[I]</math>を、<math>\mathbb{R}^{3}</math>内の曲面片とする。<br />
また、Xを、Sの近傍で定義された、区分的に滑らかなベクトル場とする。このとき、<br />

==== (1)普通の面積分 ====
:定義は、以下の通りである。
::<math>
{\int}_{\mathbf{s}\in S} \mathbf{X}(\mathbf{s}) \ {d}^{2}S
:={\int}_{({u}_{1},{u}_{2})\in I}
\langle \mathbf{X}(\varphi({u}_{1},{u}_{2}))\ |
\ \left(
\frac{\partial \varphi }{\partial{u}_{1}}
\times
\frac{\partial \varphi }{\partial{u}_{2}}
\right)
\rangle
\ d{u}_{1}d{u}_{2}</math>

:上式の右辺は、<math>({u}_{1},{u}_{2})</math>についてのスカラー値関数
:::<math>
\langle \mathbf{X}(\varphi({u}_{1},{u}_{2}))\ |
\ \left(
\frac{\partial \varphi }{\partial{u}_{1}}
\times
\frac{\partial \varphi }{\partial{u}_{2}}
\right)
\rangle
</math>
:を、区間I上で[[重積分]]したものを意味する。今、Sの単位法線ベクトル<math>{\mathbf{n}}_{S}</math>を、
::<math>
{\mathbf{n}}_{S}:=
\frac{1}{
\left|
\frac{\partial \varphi }{\partial{u}_{1}}
\times
\frac{\partial \varphi }{\partial{u}_{2}}
\right|
}
\left(
\frac{\partial \varphi }{\partial{u}_{1}}
\times
\frac{\partial \varphi }{\partial{u}_{2}}
\right)
</math>
:と定めると、

::<math>
{\int}_{\mathbf{s}\in S} \mathbf{X}(\mathbf{s}) \ {d}^{2}S
={\int}_{({u}_{1},{u}_{2})\in I} \langle \mathbf{X}(\varphi({u}_{1},{u}_{2}))\ |\ {\mathbf{n}}_{S} \rangle
\left|
\frac{\partial \varphi }{\partial{u}_{1}}
\times
\frac{\partial \varphi }{\partial{u}_{2}}
\right|
\ d{u}_{1}d{u}_{2}</math>
:である。

==== (2)絶対値による面積分 ====
:定義は、以下の通りである。
::<math>
{\int}_{\mathbf{s}\in S} \mathbf{X}(\mathbf{s}) \ |{d}^{2}S|
:={\int}_{({u}_{1},{u}_{2})\in I}
\mathbf{X}(\varphi({u}_{1},{u}_{2}))\ 
\left|\frac{\partial \varphi }{\partial{u}_{1}} \times \frac{\partial \varphi }{\partial{u}_{2}}\right| 
\ d{u}_{1}d{u}_{2}</math>

:本記事では、絶対値による面積分の場合は、<math>\ |{d}^{2}S|</math>のように、面素に絶対値記号をつけることにする。
:右辺は、<math>({u}_{1},{u}_{2})</math>についてのベクトル値関数
::<math>\mathbf{X}(\varphi({u}_{1},{u}_{2}))\ 
\left|\frac{\partial \varphi }{\partial{u}_{1}} \times \frac{\partial \varphi }{\partial{u}_{2}}\right|</math>
:を、成分ごとに区間I上で[[重積分]]したものを意味する。即ち、絶対値による面積分は「関数の面積分を各成分ごとにやる」というのと同じ意味である。

==== (3)外積面積分 ====
:定義は、以下の通りである。
::<math>
{\int}_{\mathbf{s}\in S} \mathbf{X}(\mathbf{s}) \times {d}^{2}S
:={\int}_{({u}_{1},{u}_{2})\in I} 
\left(
\mathbf{X}(
\varphi(
{u}_{1},{u}_{2})
)\right)\times
\left(
\frac{\partial \varphi }{\partial{u}_{1}} \times \frac{\partial \varphi }{\partial{u}_{2}}
\right)\ 
d{u}_{1}d{u}_{2}
</math>

:本記事では、外積面積分の場合は、<math>\times\ {d}^{2}S</math>のように、面素にの前に、×をつけることにする。
:右辺は、<math>({u}_{1},{u}_{2})</math>についてのベクトル値関数
::<math>\mathbf{X}(\varphi({u}_{1},{u}_{2}))\times
\left(\frac{\partial \varphi }{\partial{u}_{1}} \times \frac{\partial \varphi }{\partial{u}_{1}}\right)</math>

:を、成分ごとに区間I上で[[重積分]]したものを意味する。即ち、(2)との関係でいえば、
::<math>
{\int}_{\mathbf{s}\in S} \mathbf{X}(\mathbf{s}) \times {d}^{2}S
={\int}_{\mathbf{s}\in S} \mathbf{X}(\mathbf{s}) \times \mathbf{n}_{S}|{d}^{2}S|
</math>
となる。

== 脚注 ==
{{脚注ヘルプ}}
=== 注釈 ===
<references group="注釈" />

=== 参考文献 ===
{{参照方法|section=1|date=2014年1月22日 (水) 04:35 (UTC)}}
{{Reflist|2}}

== 関連項目 ==
* [[磁場]]
* [[静電気学]]
* [[ベクトルポテンシャル]]
* [[遅延ポテンシャル]]
* [[ジェフィメンコ方程式]]
* [[アンペールの法則]]
* [[ビオ・サバールの法則]]
* [[クーロンの法則]]
* [[電場|電場の強さ]]('''E''')、[[磁束密度]]('''B''')、[[電束密度]]('''D''')
* [[E-B対応とE-H対応]]
* [[地磁気]]
* [[ホール素子]]
* [[SQUID]]
* [[ファラデー効果]]
* [[マクスウェルの方程式]]
* [[ローレンツ力]]
* [[フレミングの法則]]

== 外部リンク ==
* {{PDFlink|[http://park.itc.u-tokyo.ac.jp/kato-yusuke-lab/nagai/note_0519_jika.pdf 磁化Mと磁場Hについて]}}
* {{PDFlink|[http://maildbs.c.u-tokyo.ac.jp/~torii/lectures/EM/EM_A_p25-34.pdf 磁場]}}
* {{PDFlink|[http://www.if.t.u-tokyo.ac.jp/~mita/EM/TD20101220-kaito.pdf 電気磁気学演習解答例]}}

{{DEFAULTSORT:せいしは}}
[[Category:電磁気学]]
[[Category:物理量]]
[[Category:磁気]]
[[Category:ポテンシャルエネルギー]]