マクスウェルの方程式
マクスウェルの方程式(マクスウェルのほうていしき、英: Maxwell's equations)は、電磁場のふるまいを記述する古典電磁気学の基礎方程式である。マイケル・ファラデーが幾何学的考察から見出した電磁力に関する法則が1864年にジェームズ・クラーク・マクスウェルによって数学的形式として整理された[1]。マクスウェル-ヘルツの電磁方程式、電磁方程式などとも呼ばれ、マクスウェルはマックスウェルとも表記される。
真空中の電磁気学に限れば、マクスウェルの方程式の一般解は、ジェフィメンコ方程式として与えられる。
なお、電磁気学の単位系は、国際単位系に発展したMKSA単位系のほか、ガウス単位系などがあるが、以下では原則として、国際単位系を用いることとする。
Contents
4つの方程式
(微分形による)マクスウェルの方程式は、以下の4つの連立偏微分方程式である。
- [math]\begin{cases} \nabla\cdot\boldsymbol{B}(t,\boldsymbol{x})&=0\\ \nabla\times\boldsymbol{E}(t,\boldsymbol{x})+\dfrac{\partial\boldsymbol{B}(t,\boldsymbol{x})}{\partial t}&=0\\ \nabla\cdot\boldsymbol{D}(t,\boldsymbol{x})&=\rho(t,\boldsymbol{x})\\ \nabla\times\boldsymbol{H}(t,\boldsymbol{x})-\dfrac{\partial\boldsymbol{D}(t,\boldsymbol{x})}{\partial t}&=\boldsymbol{j}(t,\boldsymbol{x}) \end{cases}[/math]
ここで E は電場の強度、B は磁束密度、D は電束密度、H は磁場の強度を表す。
また ρ は電荷密度、j は電流密度を表す。記号「∇·」、「∇×」はそれぞれベクトル場の発散 (div) と回転 (rot) である。
次に、4つの個々の方程式(成分表示で8つの式、テンソル表示で2つの式)について説明する。
磁場の構造(磁束保存の式)
- [math]\nabla\cdot\boldsymbol{B}=0[/math](微分形の磁束保存の式)
積分形で表すと次の式になる。
- [math]\int_A\boldsymbol{B}\cdot\mathrm{d}\boldsymbol{A}=0[/math]
ここで B は磁束密度(単位はテスラ T )、dA は、領域の外側へ向かう方向と直交する閉じた曲面 A 上の微小な方形の領域である。
与えられるどんな体積要素についても、表面 A の外側の点のベクトル成分の総和が内側の点のベクトル成分の総和に等しくなり、構造的に見て磁力線が閉曲線でなければならないことを意味するためである。この式は電場の積分形と同様に、閉曲面上を積分したときにのみ意味がある。
これらの式は、磁気単極子(モノポール)が存在しないことを前提としており、もし磁気単極子が発見されたならば、上の式は次のように変更されなければならない。
- [math]\nabla\cdot\boldsymbol{B}=\rho_\mathrm{m}[/math]
ここで ρm は磁気単極子の磁荷密度である。
変化する磁場と電場(ファラデー-マクスウェルの式)
- [math]\nabla\times\boldsymbol{E}=-\frac{\partial\boldsymbol{B}}{\partial t}[/math](微分形のファラデー-マクスウェルの式)
この式を積分形で表すと次の式になる。
- [math]V=-\frac{\mathrm{d}\phi_{\boldsymbol{B}}}{\mathrm{d}t}[/math](レンツの法則)
ただし、
- [math]\phi_{\boldsymbol{B}}=\int_{A}\boldsymbol{B}\cdot\mathrm{d}\boldsymbol{\mathit{A}}[/math]
ここで φB は磁束保存の式で記述された面積 A を通過する磁束、V は面積 A の縁の周囲の起電力である。
磁束保存の式の説明で述べたように、閉じた曲面を通る磁束の総和は常に 0 であるため、この式は閉じていない曲面 A についてのみ働く。起電力はその曲面 A の縁に沿って測定されるが、閉じた曲面には縁がない。いくつかの電気工学の文献では、曲面 A の縁に巻かれたコイルの数 N を磁束の導関数の前に用いてこの積分形式を表現している。なお、式中の負号があるため、磁束密度の時間微分が正なら左回転に、負なら右回転になる。
電荷密度と電場(マクスウェル-ガウスの式)
- [math]\nabla\cdot\boldsymbol{D}=\rho[/math](微分形のマクスウェル-ガウスの式: D-H 対応)
ここで、ρ は、電荷密度(単位は C/m3)。D は電束密度(単位は C/m2)で、「線形な物質中」では D = εE の関係がある(E は電場の強度、ε は誘電率)。電場が非常に強くない限り、どんな物質も「線形」なものとして扱うことができる。上の式は、電束は電荷の存在するところで発生・消滅し、それ以外のところでは保存されることを意味している。真空の誘電率は ε0 と書かれ、次の式で表される。
- [math]\nabla\cdot\boldsymbol{E}=\frac{\rho}{\varepsilon_0}[/math](微分形のマクスウェル-ガウスの式: E-B対応)
また、εr = εε0 で定義される比誘電率などが用いられることもある。
マクスウェル-ガウスの式を積分形で表すと次の式になる。
- [math]\int_A\boldsymbol{E}\cdot\mathrm{d}\boldsymbol{A}=\frac{Q_{\rm encl}}{\varepsilon_0}[/math]
ここで dA は、電荷の外側へ向かう方向と直交する閉じた曲面 A 上の微小な方形の領域であり、Qencl はその閉曲面当たりの電荷である。この積分形は、閉曲面上を積分したときにのみ意味があり、ガウスの法則としてよく知られている。
電流・電場と磁場(アンペール-マクスウェルの式)
- [math]\nabla\times\boldsymbol{H}=\boldsymbol{j}+\frac{\partial\boldsymbol{D}}{\partial t}[/math](微分形のアンペール-マクスウェルの式: D-H対応)
ここで j は電流密度。H は磁場の強度(単位は A/m)で、「線形な物質中」で「磁場の強度が小さい範囲」において、B = μH の関係がある(B は磁束密度、μ は透磁率)。
真空中では透磁率 μ は真空の透磁率 μ0 = 4πテンプレート:E- W/Am で置き換えられる。したがって式は次のようになる。
- [math]\nabla\times\boldsymbol{B}=\mu_0\boldsymbol{j}+\mu_0\varepsilon_0\frac{\partial\boldsymbol{E}}{\partial t}[/math]
積分形は次のようになる。
- [math]\oint_s\boldsymbol{B}\cdot\mathrm{d}\boldsymbol{s}=\mu_0 I_{\rm{through}}-\mu_0\varepsilon_0\int_A\frac{\partial\boldsymbol{E}}{\partial t}\cdot\mathrm{d}\boldsymbol{A}[/math]
s は開曲面 A の縁となる曲線で、Ithrough は曲線 s で囲まれた曲面 A を通過する電流 Ithrough A = ∫A j · dA である。コンデンサや ∇ · j ≠ 0 となるほかの場所がなければ、右辺の第 2 項(変位電流)は一般に無視される。
それぞれの式の解釈
- 磁束保存の式
- 磁力線はどこかを起点とすることも終点とすることもできない、すなわち磁気単極子(モノポール)が存在しないことを示唆している。
- ファラデー-マクスウェルの式
- 磁場の時間変化があるところには電場が生じることを示すファラデーの電磁誘導の法則の定式化である。
- 電動機(モータ)や発電機など非常に多くの実用的な応用に関係している。
- ガウス-マクスウェルの式
- 電場の源は電荷であることを示し、クーロンの法則に相当する。
- アンペール-マクスウェルの式
- 電場の時間変化(変位電流という)と電流とで磁場が生じることを示す。
- これらの式は、電流によって磁場が生じるというアンペールの法則に電場の時間変化を加えたものである。
歴史的経緯
マクスウェルの方程式は、次の2つの組に分類されることが多い。
電磁場の拘束条件
第1の組は、 テンプレート:Numblk テンプレート:Numblk である。この式は電磁場の拘束条件を与える式である(ビアンキ恒等式)。
この式は E, B を電磁ポテンシャル φ, A により、 テンプレート:Numblk テンプレート:Numblk と表せば恒等的に満たすように出来る。
マクスウェル自身の原著論文『電磁場の動力学的理論』(1865年)や原著教科書『電気磁気論』(1873年)では上記のように表されていたが、1890年になってヘルツが改めて理論構成を考察し、上記2式から電磁ポテンシャルを消去しテンプレート:EquationNote, テンプレート:EquationNote を基本方程式とすることを要請した。このヘルツによる電磁ポテンシャルを消去した形をマクスウェルの方程式と見なすのが現在の主流となっている。この見かたでは テンプレート:EquationNote と テンプレート:EquationNote は電磁場の定義式と見なされる。
電磁場の運動方程式
第2の組は、 テンプレート:Numblk テンプレート:Numblk である。電荷、電流の分布が電磁場の源となっていることを表す式である(電磁場の運動方程式)。
電磁場の微分(左辺)が電荷、電流の分布(右辺)によって書かれており、電荷、電流の分布を与えると電磁場の形が分かる方程式になっている。
この式から、電荷、電流の分布には電気量保存則(連続の方程式)
[math]\frac{\partial\rho}{\partial t}+\nabla\cdot\boldsymbol{j}=0[/math]
が成り立つことが導かれる。
また、電磁場はローレンツ力
[math]\rho\boldsymbol{E}+\boldsymbol{j}\times\boldsymbol{B}[/math]
により電荷、電流の分布を変動させる。
それぞれの組は時間微分を片側に移し、
[math]\frac{\partial\boldsymbol{B}}{\partial t}=-\nabla\times\boldsymbol{E},~\nabla\cdot\boldsymbol{B}=0[/math]
[math]\frac{\partial\boldsymbol{D}}{\partial t}=\nabla\times\boldsymbol{H}+\boldsymbol{j},~\nabla\cdot\boldsymbol{D}=\rho[/math]
と変形すれば、時間発展の方程式とその初期条件と見ることができる。
また、E, B により記述した場合は式の組み合わせを変えて
[math]\nabla\cdot\boldsymbol{E}=\frac{\rho}{\epsilon},~\nabla\times\boldsymbol{E}=-\frac{\partial\boldsymbol{B}}{\partial t}[/math]
[math]\nabla\cdot\boldsymbol{B}=0,~\nabla\times\boldsymbol{B}=\mu\epsilon\frac{\partial\boldsymbol{E}}{\partial t}+\mu\boldsymbol{j}[/math]
として、ベクトル場の発散と回転を与える式と見ることができる。
これらの方程式系に整理されたことから、電場と磁場の統一(電磁場)、光が電磁波であることなどが導かれ、その時空論としての特殊相対性理論に至る。後年、アインシュタインは特殊相対性理論の起源はマクスウェルの電磁場方程式である旨を明言している。
マクスウェルが導出した方程式はベクトルの各成分をあたかも互いに独立な量であるかのように別々の文字で表して書かれており、現代の洗練された形式ではなかった。これを1884年にヘヴィサイドがベクトル解析の記法を適用して現在の見やすい形に書き改めた。しかも彼は既にそこで電磁ポテンシャルが消去出来ることを示して、方程式系を今日我々が知る形に整理していた。しかし、その意義は直ちには認められるに至らず、それとは独立に上記のヘルツの仕事がなされた。
ベクトル記法が一般化し始めるのは 1890年代半ばであって、ヘルツの論文ではまだそれを使っていない。いずれにせよ、このベクトル解析の記法の採用は場における様々な対称性を一目で見ることを可能にし、物理現象の理解に大いに役立った[2]。
E-B による表示
電束密度 D と磁場の強度 H は、電磁場の媒質中での振舞いを表現する量である。媒質の性質を他の部分で表現することにより、電場の強度 E と磁束密度 B により方程式を記述することが出来る。
真空中
媒質が存在しない真空中(自由空間中)においては、D, H はそれぞれ E, B と
[math]\boldsymbol{D}=\varepsilon_0\boldsymbol{E}[/math]
[math]\boldsymbol{H}=\frac{1}{\mu_0}\boldsymbol{B}[/math]
の関係にある。ここで ε0 は真空の誘電率、μ0 は真空の透磁率であり、これらは普遍定数である。このとき方程式は
[math]\nabla\cdot\boldsymbol{E}=\frac{1}{\varepsilon_0}\rho[/math]
[math]\nabla\times\boldsymbol{B}-\mu_0\varepsilon_0\frac{\partial\boldsymbol{E}}{\partial t}=\mu_0\boldsymbol{j}[/math]
となって E, B により記述される。
等方一様線型媒質中
線型媒質中においては、D, H は線型関係
[math]\boldsymbol{D}=\varepsilon\boldsymbol{E}[/math]
[math]\boldsymbol{B}=\mu\boldsymbol{H}[/math]
によって E, B と関係付けられる。ここで ε, μ はそれぞれその媒質の誘電率と透磁率であり、媒質の性質を特徴付ける物性値である。これらは一般にはテンソルであるが、等方的な媒質ではスカラーとなる。さらに一様な媒質であると考えれば、マクスウェルの方程式は
[math]\nabla\cdot\boldsymbol{E}=\frac{\rho}{\varepsilon}[/math]
[math]\nabla\times\boldsymbol{B}-\mu\varepsilon\frac{\partial\boldsymbol{E}}{\partial t}=\mu\boldsymbol{j}[/math]
と変形できて E, B により記述される。
一般の媒質中
一般の媒質中においては、D, H と E, B を関係付ける量として、誘電分極 P と磁化 M が
[math]\boldsymbol{D}=\varepsilon_0\boldsymbol{E}+\boldsymbol{P}[/math]
[math]\boldsymbol{B}=\mu_0(\boldsymbol{H}+\boldsymbol{M})[/math]
によって導入される。このとき、方程式は
[math]\nabla\cdot\boldsymbol{E}=\frac{1}{\varepsilon_0}(\rho-\nabla\cdot\boldsymbol{P})[/math]
[math]\nabla\times\boldsymbol{B}-\mu_0\varepsilon_0\frac{\partial\boldsymbol{E}}{\partial t}=\mu_0\left(\boldsymbol{j} +\frac{\partial\boldsymbol{P}}{\partial t}+\nabla\times\boldsymbol{M}\right)[/math]
となる。さらに分極電荷密度、分極電流密度、磁化電流密度を
[math]\rho_P=-\nabla\cdot\boldsymbol{P}[/math]
[math]\boldsymbol{j}_P=\frac{\partial\boldsymbol{P}}{\partial t}[/math]
[math]\boldsymbol{j}_M=\nabla\times\boldsymbol{M}[/math]
として導入すれば、方程式は
[math]\nabla\cdot\boldsymbol{E}=\frac{1}{\varepsilon_0}(\rho +\rho_P)=\frac{1}{\varepsilon_0}\rho_0[/math]
[math]\nabla\times\boldsymbol{B}-\mu_0\varepsilon_0\frac{\partial\boldsymbol{E}}{\partial t}=\mu_0\left(\boldsymbol{j} +\boldsymbol{j}_P+\boldsymbol{j}_M\right)=\mu_0\boldsymbol{j}_0[/math]
となる。これは真空中におけるマクスウェルの方程式と同じ形をしている。媒質は原子核や電子などの荷電粒子から構成されている。これらが真空中に分布しているものとして考えたときの電荷の分布が ρ0, j0 である。一方、媒質を構成する荷電粒子はマクロに見たとき、分子、或いは原子として束縛されている。電荷の分布を平均化した後の分布が ρ, j である。
ローレンツゲージでのマクスウェルの方程式
以下のローレンツ条件
- [math]\operatorname{div}\boldsymbol{A}+\frac{1}{c^2}\frac{\partial\phi}{\partial t}=0[/math]
における電磁ポテンシャル(ベクトルポテンシャル A と電磁ポテンシャル φ)を用いて、マクスウェルの方程式は以下の2組の方程式として表すことができる。
- [math]\begin{cases} \left(\triangle-\dfrac{1}{c^2}\dfrac{\partial^{2}}{\partial t^2}\right)\phi&=-\frac{\rho}{\varepsilon_0}\\ \left(\triangle-\dfrac{1}{c^2}\dfrac{\partial^{2}}{\partial t^2}\right)\boldsymbol{A}&=-\mu_0\boldsymbol{J} \end{cases}[/math]
いずれの式も左辺は線形の演算子であるところのダランベルシアンが作用しており、右辺は片やスカラー値の、片やベクトル値の連続関数である。ベクトルについては各々の成分について適用して考えることでスカラーの場合と同様に考えることができる。線形微分方程式に対してはグリーン関数法を考えることで解くことができる。すなわち、
[math]\left(\Delta-\frac{1}{c^{2}}\frac{\partial^{2}}{\partial t^{2}}\right)G({\bf x},t)=-\delta({\bf x},t)[/math]
の解となる関数(グリーン関数)[math]G({\bf x},t)[/math]を求めることで一般に
[math]\left(\Delta-\frac{1}{c^{2}}\frac{\partial^{2}}{\partial t^{2}}\right)f({\bf x},t)=-\rho({\bf x},t)[/math]
なる方程式に対して
[math]f({\bf x},t)=\int d^{3}x' dt'\ G({\bf x}-{\bf x'},t-t')\rho({\bf x'},t')[/math]
として求めることができる。このときのグリーン関数は先進グリーン関数と遅延グリーン関数の2つを得るが、物理的に意味のある遅延グリーン関数を採用することで遅延ポテンシャルを得ることができる。
遅延ポテンシャルを元に電場や磁場を計算するのが一般に運動している物体についての電磁場を検討する際に楽な方法であり、結果としてジェフィメンコ方程式を得ることになる。
電磁波の波動方程式
マクスウェルの方程式により、電磁波の伝搬を記述することができる[3]。真空または電荷分布がない絶縁体では、電場と磁場が波動方程式
- [math] \nabla^2 \boldsymbol{E}-\mu \epsilon \frac{\partial^2\boldsymbol{E}}{\partial t^2}=0[/math]
- [math] \nabla^2 \boldsymbol{B}-\mu \epsilon \frac{\partial^2\boldsymbol{B}}{\partial t^2}=0[/math]
を満たすことがマクスウェル方程式から示される。これは電磁場が媒質中を速さ
- [math]v=\frac{1}{\sqrt{\mu \epsilon}}[/math]
- [math]\eta=\sqrt{\frac{\mu \epsilon}{\mu_0 \epsilon_0}}=c\sqrt{\mu \epsilon}[/math]
を導入すれば、v は
- [math]v=\frac{c}{\eta}[/math]
とも表される。
ここで、真空の誘電率と真空の透磁率の各値から導かれる定数 c の値が光速度の値とほとんど一致することから、マクスウェルは光は電磁波ではないかという予測を行なったのである。その予測は1888年にハインリヒ・ヘルツによって実証される。ヘルツはマクスウェルの方程式の研究に貢献したので、マクスウェルの方程式はマクスウェル-ヘルツの(電磁)方程式と呼ばれることもある。
マクスウェルの方程式と特殊相対性理論
19世紀後半を通じて物理学者の大半は、マクスウェルの方程式において光速度が全ての観測者に対して不変であるという奇妙な予測のために、またそれがニュートン力学の運動法則と矛盾したために、これらの方程式が電磁場への近似的なものに過ぎないと考えた。しかし、1905年にアインシュタインが特殊相対性理論を提出したことによって、マクスウェルの方程式が正確で、ニュートン力学の方を修正すべきだったことが明確になった。これら電磁場の方程式は、特殊相対性理論と密接な関係にあり、ローレンツ変換に対する不変性(共変性)を満たす。磁場の方程式は、光速度に比べて小さい速度では、相対論的変換による電場の方程式の変形に結び付けられる。
電場と磁場による表現では、共変性が見にくいため、4元ポテンシャル Aμ を考える。
[math]A^\mu=(\phi/c,\boldsymbol{A}),~A_\mu=\eta_{\mu\nu}A^\nu=(\phi/c,-\boldsymbol{A})[/math]
但し、重複するギリシャ文字に対してはアインシュタインの縮約記法に従って和をとるものとし、計量テンソルは ημν = diag(1, −1, −1, −1) で与えるものとする。また、各ギリシャ文字は 0,1,2,3 の値を取り、0は時間成分、1,2,3は空間成分を表すものとする。特に時空の座標については (x0, x1, x2, x3) = (ct, x, y, z) である。
電磁ポテンシャルから構成される電磁場テンソル
[math]F_{\mu\nu}\equiv\partial_\mu A_\nu-\partial_\nu A_\mu=-F_{\nu\mu}[/math](テンプレート:EquationNote,テンプレート:EquationNoteに対応)
を導入する。電場、磁場との対応関係は
[math](F_{01},F_{02},F_{03})=(E_1/c, E_2/c, E_3/c),~(F_{32},F_{13},F_{21})=(B_1,B_2,B_3)[/math]
となる。
このとき、マクスウェル方程式はローレンツ変換に対しての共変性が明確な形式で、次のような2つの方程式にまとめられる。
[math]\partial_\rho F_{\mu\nu}+\partial_\mu F_{\nu\rho}+\partial_\nu F_{\rho\mu}=0[/math](テンプレート:EquationNote,テンプレート:EquationNoteに対応)
[math]\partial_{\mu}F^{\mu\nu}=\mu_0 j^\nu[/math](テンプレート:EquationNote,テンプレート:EquationNoteに対応)
但し、jμ は4元電流密度
[math]j^\mu=(c\rho,\boldsymbol{j})[/math]
である。このとき、電荷の保存則は
[math]\partial_\mu j^\mu=0[/math](3に対応)
と表される。なお、4元ポテンシャルで表現すると、マクスウェル方程式は次の一つの方程式にまとめられる。
[math]\Box A^\mu -\partial^\mu\partial_\nu A^\nu=-\mu_0 j^\mu[/math]
ここで、□はダランベール演算子を表す。
微分形式による表現
マクスウェルの方程式は多様体理論における微分形式によって簡明に表現することができる[4]。
まず電磁ポテンシャル Aμ により、1次微分形式
[math]A=A_\mu dx^\mu=\phi\,dt-A_x\,dx-A_y\,dy-A_z\,dz[/math]
を導入する。これに外微分を作用させることで2次微分形式
[math]\begin{align}F&\equiv dA=\tfrac{1}{2}(\partial_\mu A_\nu-\partial_\nu A_\mu)\,dx^\mu\wedge dx^\nu\\ &=\tfrac{1}{2}F_{\mu\nu}\, dx^\mu\wedge dx^\nu\\ &=E_x\,dt\wedge dx+E_y\,dt\wedge dy+E_z\,dt\wedge dz-B_x\, dy\wedge dz-B_y\, dz\wedge dx-B_z\, dx\wedge dy\end{align}[/math]
が定義される。 さらに F のホッジ双対として 2次微分形式
[math]\begin{align}H&\equiv \tfrac{1}{\mu_0}F^*=\tfrac{1}{4\mu_0}\epsilon_{\mu\nu\rho\sigma}F^{\mu\nu}dx^\rho\wedge dx^\sigma\\ &=\tfrac{1}{2}H_{\mu\nu}\,dx^\mu\wedge dx^\nu\\ &=H_x\,cdt\wedge dx+H_y\,cdt\wedge dy+H_z\,cdt\wedge dz+D_xc\,dy\wedge dz+D_yc\,dz\wedge dx+D_zc\,dx\wedge dy\end{align}[/math]
が定義される。
4元電流密度により1次微分形式
[math]J=j_\mu dx^\mu=\rho c^2dt-j_x dx-j_y dy-j_z dz[/math]
を導入し、これのホッジ双対により3次微分形式
[math]\begin{align}J^*&=\tfrac{1}{3!}\epsilon_{\mu\nu\rho\sigma}j^\mu dx^\nu\wedge dx^\rho\wedge dx^\sigma\\ &=\rho c\,dx \wedge dy\wedge dz-j_x\,cdt \wedge dy\wedge dz-j_y\,cdt\wedge dz\wedge dx-j_z\,cdt\wedge dx\wedge dy\end{align}[/math]
を定義すれば、外微分の作用により運動方程式(テンプレート:EquationNote,テンプレート:EquationNote)に対応して
[math]dH=J^*[/math]
となる。
外微分の性質 ddξ=0 から(テンプレート:EquationNote,テンプレート:EquationNote)に対応する
[math]dF=ddA=0[/math]
と、連続の方程式に対応する
[math]dJ^*=ddH=0[/math]
が得られる。
脚注
- 出典
参考文献
原論文
- Maxwell, J. C. (October 27, 1865). “A dynamical theory of the electromagnetic field [電磁場の動力学的理論]” (PDF). Phil. Trans. R. Soc. Lond. 155: 459-512. doi:10.1098/rstl.1865.0008. JSTOR 108892 .
書籍
- Lorentz, H.A. 『ローレンツ 電子論』 広重徹、1973。
- 広重, 徹 『物理学史II』 培風館〈新物理学シリーズ〉、1968-03。全国書誌番号:68001733。ISBN 978-4563024062。OCLC 673599647。
- 『場の古典論:電気力学, 特殊および一般相対性理論』 恒藤敏彦, 広重徹訳、東京図書〈ランダウ=リフシッツ理論物理学教程〉、1978-10、原書第6版。全国書誌番号:79000237。ISBN 978-4489011610。OCLC 841897028。
- 砂川, 重信 『理論電磁気学』 紀伊國屋書店、1999-09、第3版。全国書誌番号:99125994。ISBN 978-4314008549。OCLC 675159672。
- Jackson, J.D. 『電磁気学』上巻、西田稔訳、吉岡書店〈物理学叢書〉、2002-07、原書第3版。全国書誌番号:20301816。ISBN 978-4842703053。OCLC 123038116。
- (1989) Differential Forms with Applications to the Physical Sciences. Dover Publications. ISBN 0486661695.