これまでのところをまとめておく。
マクスウェル方程式
\begin{eqnarray}
\nabla\cdot \boldsymbol{D} &=& \rho \quad \tag{1}\\
\nabla\cdot\boldsymbol{B} &=& 0 \quad \tag{2}\\
\nabla\times\boldsymbol{E} + \color{red}{\frac{\partial \boldsymbol{B}}{\partial t} }&=& \boldsymbol{0} \quad \tag{3}\\
\nabla\times\boldsymbol{H} – \color{red}{\frac{\partial \boldsymbol{D}}{\partial t} }&=& \boldsymbol{J} \quad \tag{4}\end{eqnarray}
時間変動しない電磁場の場合は赤字の部分が消えて
\begin{eqnarray}
\nabla\cdot \boldsymbol{D} &=& \rho \quad \tag{1}\\
\nabla\cdot\boldsymbol{B} &=& 0 \quad \tag{2}\\
\nabla\times\boldsymbol{E} &=& \boldsymbol{0} \quad \tag{3a}\\
\nabla\times\boldsymbol{H} &=& \boldsymbol{J} \quad \tag{4a}\end{eqnarray}
静電場の基本方程式
$\boldsymbol{D} = \varepsilon_0 \boldsymbol{E}$ を使って $\boldsymbol{E}$ だけの式にして，$(\mbox{1})$ 式と $(\mbox{3a})$ 式は
$$\nabla\cdot \boldsymbol{E} = \frac{\rho}{\varepsilon_0},
\quad \nabla\times\boldsymbol{E} = \boldsymbol{0}$$
静磁場の基本方程式
$\boldsymbol{H} = \frac{1}{\mu_0} \boldsymbol{B}$ を使って $\boldsymbol{B}$ だけの式にして，$(\mbox{2})$ 式と $(\mbox{4a})$ 式は
$$ \nabla\cdot\boldsymbol{B} = 0, \quad \nabla\times\boldsymbol{B}= \mu_0\,\boldsymbol{J} = \frac{\boldsymbol{J} }{\varepsilon_0 c^2}$$
$\displaystyle \varepsilon_0 \mu_0 = \frac{1}{c^2}$ となる理屈はこのへんにまとめておいた。
電磁ポテンシャル
$$\nabla\times\boldsymbol{E} = \boldsymbol{0} \quad\Rightarrow\quad \boldsymbol{E} \equiv – \nabla\phi$$
$$\nabla\cdot\boldsymbol{B} = 0 \quad\Rightarrow\quad \boldsymbol{B} \equiv \nabla\times\boldsymbol{A}$$
というようにスカラーポテンシャル（静電ポテンシャル）\(\phi\) とベクトルポテンシャル \(\boldsymbol{A}\) を導入できた。
\(\phi\) と \(\boldsymbol{A}\) をまとめて電磁ポテンシャルと呼ぼう。また，今後，静電場（時間変動しない電場）でない場合についても述べるので，一時「静電」ポテンシャルとも呼んでいた \(\phi\) の呼称を「スカラーポテンシャル」に統一する。
ゲージ条件
ベクトルポテンシャル \(\boldsymbol{A}\) に対して，クーロンゲージ条件
$$\nabla\cdot\boldsymbol{A} = 0$$
を課すと…
ポアソン方程式
電磁ポテンシャルは以下のようなポアソン方程式に従うことがわかり，
$$\nabla^2 \phi = – \frac{\rho}{\varepsilon_0}, \quad \nabla^2 \boldsymbol{A} = -\frac{\boldsymbol{J}}{\varepsilon_0c^2}$$
ただちに完全な解が以下のように求まるのであった。
\begin{eqnarray}
\phi(\boldsymbol{r}) &=& \frac{1}{4\pi\varepsilon_0} \iiint \frac{\rho(\boldsymbol{r}’)}{|\boldsymbol{r} – \boldsymbol{r}’|} dV’
\end{eqnarray}
$$ \boldsymbol{A} (\boldsymbol{r})= \frac{1}{4\pi \varepsilon_0 c^2} \iiint \frac{\boldsymbol{J}(\boldsymbol{r}’)}{|\boldsymbol{r} – \boldsymbol{r}’|} \,dV’$$
電荷密度から直接静電場を求める式
\begin{eqnarray}
\boldsymbol{E} (\boldsymbol{r}) &=& – \nabla\phi (\boldsymbol{r})
= \frac{1}{4\pi \varepsilon_0} \iiint\frac{\rho(\boldsymbol{r}’) (\boldsymbol{r} – \boldsymbol{r}’) }{|\boldsymbol{r} – \boldsymbol{r}’|^3}dV’
\end{eqnarray}
この積分が（比較的）簡単にできるいくつかの例を以下にまとめる。
点電荷の電荷密度と電場
位置 $\boldsymbol{r}_1$ においた点電荷 $q_1$ の電荷密度は，ディラックのデルタ関数を使って
$$\rho(\boldsymbol{r}) = q_1 \, \delta^3(\boldsymbol{r} – \boldsymbol{r}_1)$$
電場は
$$\boldsymbol{E} = \frac{q_1}{4\pi \varepsilon_0} \frac{ \boldsymbol{r} – \boldsymbol{r}_1}{|\boldsymbol{r} – \boldsymbol{r}_1|^3} $$
特に点電荷が原点にある場合は $\boldsymbol{r}_1 = \boldsymbol{0}$ とすればいいし，複数の点電荷の場合は重ね合わせの原理から
$$\boldsymbol{E} = \sum_i \frac{q_i}{4\pi \varepsilon_0} \frac{ \boldsymbol{r} – \boldsymbol{r}_i}{|\boldsymbol{r} – \boldsymbol{r}_i|^3} $$
電気双極子による電場
電気双極子とは正電荷 $q$ と負電荷 $-q$ が微小距離 $\boldsymbol{d}$ で対になっている状態。電気双極モーメント $\boldsymbol{p} \equiv q \boldsymbol{d}$ による遠方の電場は
$$\boldsymbol{E}= \frac{1}{4\pi\varepsilon_0}\left\{ 3 \frac{\boldsymbol{r}\cdot \boldsymbol{p} }{r^5} \boldsymbol{r} – \frac{\boldsymbol{p}}{r^3}\right\}$$
軸対称な電荷分布による電場
$z$ 軸について軸対称な電荷分布を表す電荷密度は
$$\rho(\boldsymbol{r}) = \rho(\varrho), \quad\boldsymbol{\varrho}\equiv (x, y, 0), \ \ \varrho\equiv \sqrt{x^2+y^2}$$
電場は
$$\boldsymbol{E} = \frac{Q_{\varrho}}{2\pi \varepsilon_0} \frac{\boldsymbol{\varrho}}{{\varrho}^2}$$
ここで $Q_{\varrho}$ は，$z$ 軸を中心とし，底面積が $\pi \varrho^2$ の単位高さの円柱内の全電荷である。特に $z$ 軸上の一様な線電荷密度 $\lambda$ の場合は，$Q_{\varrho} \rightarrow \lambda$ とすればよい。
面対称な電荷分布による電場
$x = 0$ の $yz$ 平面について面対称な電荷分布を表す電荷密度は
$$\rho(\boldsymbol{r}) = \rho(|x|)$$
電場は
$$ E_x = \frac{Q_{|x|}}{2 \epsilon_0} \frac{x}{|x|}, \quad E_y = E_z = 0$$
$Q_{|x|}$ は $yz$ 面が単位面積，$x$ 方向が $-x$ から $x$ までの直方体内の全電荷。特に，$yz$ 面上の一様面電荷密度 $\sigma$ の場合は，$Q_{|x|} \rightarrow \sigma$ とすればよい。
球対称な電荷分布による電場
原点を中心とした球対称な電荷分布を表す電荷密度は
$$\rho(\boldsymbol{r}) = \rho(r), \quad r \equiv \sqrt{x^2 + y^2 + z^2}$$
電場は
$$\boldsymbol{E} = \frac{Q_r}{4\pi \varepsilon_0} \frac{\boldsymbol{r}}{r^3}$$
$Q_r$ は原点を中心とする半径 $r$ の球内の全電荷。球対称な電荷分布の場合は，この $Q_r$ をもつ点電荷を原点においたときの電場と等価であることがわかるのであった。
電流密度から直接静磁場を求める式：ビオ – サバールの法則
\begin{eqnarray}
\boldsymbol{B} (\boldsymbol{r}) = \nabla\times \boldsymbol{A} (\boldsymbol{r})&=& \frac{1}{4\pi \varepsilon_0 c^2} \iiint \frac{\boldsymbol{J}(\boldsymbol{r}’)\times (\boldsymbol{r} – \boldsymbol{r}’)}{|\boldsymbol{r} – \boldsymbol{r}’|^3} \,dV’
\end{eqnarray}
この積分が（比較的）簡単にできるいくつかの例を以下にまとめる。
直線電流による磁場
$z$ 軸上の直線電流 $\boldsymbol{I} = (0, 0, I)$ を表す電流密度 $\boldsymbol{J}$ はディラックのデルタ関数を使って
$$\boldsymbol{J} = (0, 0, J_z) = \boldsymbol{I}\delta(x) \delta(y)$$
磁場は
$$\begin{eqnarray}
\boldsymbol{B} &=& \frac{1}{2\pi \varepsilon_0 c^2} \frac{\boldsymbol{I}\times\boldsymbol{\varrho}}{\varrho^2}
\end{eqnarray}, \quad \boldsymbol{\varrho} = (x, y, 0), \ \varrho = \sqrt{x^2+y^2}$$
円電流による磁場
$xy$ 平面上の原点を中心とした半径 $a$ の円周上を流れる電流 $I$ を表す電流密度は
\begin{eqnarray}
\boldsymbol{J} &=& (J_x, J_y, 0) \\
&=& (-I \sin\phi , I \cos\phi , 0)\times \delta(\varrho-a) \,\delta(z)
\end{eqnarray}
ビオ・サバールの法則に入れて磁場を求めるが，一般に積分ができるわけではなく，あえて書くと $z$ 軸上 $\boldsymbol{r} = (0, 0, z)$ の磁場は簡単に積分できて以下のようになる。
\begin{eqnarray}
B_x(0,0,z) &=& 0\\
B_y(0,0,z)&=& 0 \\
B_z(0,0,z) &=& \frac{I a^2}{2 \varepsilon_0 c^2\left\{z^2 + a^2\right\}^{\frac{3}{2}}}
\end{eqnarray}
微小な円電流・磁気双極子による磁場
微小な円電流がつくる遠方の磁場は
$$a^2 \ll r^2 = x^2 + y^2 + z^2$$
とすると，磁気双極モーメント $\boldsymbol{\mu} \equiv (0, 0, I\pi a^2)$ を持つ磁気双極子がつくる磁場となり，
$$\boldsymbol{B} = \frac{1}{4\pi\varepsilon_0 c^2} \left\{3 \frac{\boldsymbol{r}\cdot\boldsymbol{\mu} }{r^5} \boldsymbol{r} – \frac{\boldsymbol{\mu}}{r^3} \right\}$$
ソレノイドを流れる電流による磁場
$z$ 軸を中心とし，半径 $a$ で $z$ 方向に単位長さあたり $n$ 回巻きの十分に長いソレノイドを電流 $I$ が流れる。これを表す電流密度は（ほんとは螺旋状に巻いてあるところを，円電流を$z$ 方向に単位長さあたり $n$ 個重ねて表すことにして）
\begin{eqnarray}
\boldsymbol{J} &=& (J_x, J_y, 0) \\
&=& (-n I \sin\phi , n I \cos\phi , 0)\times \delta(\sqrt{x^2 + y^2}-a)
\end{eqnarray}
磁場は
$$B_x = B_y = 0, \quad
B_z =\left\{\begin{array}{ll}\frac{n I}{\varepsilon_0 c^2} & (\sqrt{x^2 + y^2}< a)\\ \ \\0 & (\sqrt{x^2 + y^2}>a)\end{array}\right.$$
ソレノイドの外部 $\sqrt{x^2 + y^2}>a$ では $B_z$ も含めて $\boldsymbol{B} = \boldsymbol{0}$，ソレノイドの内部 $\sqrt{x^2 + y^2}<a$ では，ソレノイドに沿って一定の大きさの $B_z$ のみが存在する。
