真空中のマクスウェル方程式(真空中のマクスウェル方程式(微分形)(Definition 3.1)[マクスウェル方程式と電磁波] )を SI 単位で書きます。ρ \rho ρ は電荷密度、j \boldsymbol{j} j は電流密度、c = 1 / ε 0 μ 0 c=1/\sqrt{\varepsilon_0\mu_0} c = 1/ ε 0 μ 0 です。
div E = ρ ε 0 , rot B = μ 0 j + 1 c 2 ∂ E ∂ t , div B = 0 , rot E = − ∂ B ∂ t . \begin{aligned}
\operatorname{div}\boldsymbol{E} &= \frac{\rho}{\varepsilon_0}, &\qquad
\operatorname{rot}\boldsymbol{B} &= \mu_0\boldsymbol{j}+\frac{1}{c^2}\frac{\partial \boldsymbol{E}}{\partial t},\\[4pt]
\operatorname{div}\boldsymbol{B} &= 0, &\qquad
\operatorname{rot}\boldsymbol{E} &= -\frac{\partial \boldsymbol{B}}{\partial t}.
\end{aligned} div E div B = ε 0 ρ , = 0 , rot B rot E = μ 0 j + c 2 1 ∂ t ∂ E , = − ∂ t ∂ B . 左列上と右列上、すなわち ρ , j \rho,\boldsymbol{j} ρ , j が現れる 2 本を非斉次方程式 、div B = 0 \operatorname{div}\boldsymbol{B}=0 div B = 0 と rot E = − ∂ t B \operatorname{rot}\boldsymbol{E}=-\partial_t\boldsymbol{B} rot E = − ∂ t B の 2 本を斉次方程式 と呼びます。この記事の主題は、斉次方程式の一般解を書き下すことです。
「回転がゼロなら勾配」「発散がゼロなら回転」という 2 つの事実を、証明つきで確認します。逆向き(rot grad f = 0 \operatorname{rot}\operatorname{grad}f=\boldsymbol{0} rot grad f = 0 と div rot A = 0 \operatorname{div}\operatorname{rot}\boldsymbol{A}=0 div rot A = 0 )は成分計算で直ちに従うので、問題は逆命題です。逆命題には領域の形についての仮定が要ります。
領域 Ω ⊂ R 3 \Omega\subset\mathbb{R}^3 Ω ⊂ R 3 が点 x 0 \boldsymbol{x}_0 x 0 について星形 であるとは、任意の x ∈ Ω \boldsymbol{x}\in\Omega x ∈ Ω に対して線分 { x 0 + s ( x − x 0 ) ∣ 0 ≤ s ≤ 1 } \{\boldsymbol{x}_0+s(\boldsymbol{x}-\boldsymbol{x}_0)\mid 0\le s\le 1\} { x 0 + s ( x − x 0 ) ∣ 0 ≤ s ≤ 1 } が Ω \Omega Ω に含まれることをいいます。以下では記号を簡単にするため x 0 = 0 \boldsymbol{x}_0=\boldsymbol{0} x 0 = 0 とします(一般の場合は平行移動するだけです)。全空間 R 3 \mathbb{R}^3 R 3 、球、直方体はすべて星形です。
Lemma 2.1 (回転がゼロなら勾配(ポアンカレの補題) )
Ω ⊂ R 3 \Omega\subset\mathbb{R}^3 Ω ⊂ R 3 を原点について星形な開集合とし、F : Ω → R 3 \boldsymbol{F}\colon\Omega\to\mathbb{R}^3 F : Ω → R 3 を C 1 C^1 C 1 級のベクトル場とする。Ω \Omega Ω 上で rot F = 0 \operatorname{rot}\boldsymbol{F}=\boldsymbol{0} rot F = 0 ならば、C 2 C^2 C 2 級のスカラー場 f : Ω → R f\colon\Omega\to\mathbb{R} f : Ω → R で F = grad f \boldsymbol{F}=\operatorname{grad}f F = grad f を満たすものが存在する。具体的には
f ( x ) = ∫ 0 1 x ⋅ F ( s x ) d s f(\boldsymbol{x})=\int_0^1 \boldsymbol{x}\cdot\boldsymbol{F}(s\boldsymbol{x})\,ds f ( x ) = ∫ 0 1 x ⋅ F ( s x ) d s がその 1 つである。
Proof(Lemma 2.1) f f f を上の式で定める(星形性より s x ∈ Ω s\boldsymbol{x}\in\Omega s x ∈ Ω なので、この積分は意味を持ちます)。l l l 成分を計算します。積分記号下の微分を行い、x ⋅ F ( s x ) = x j F j ( s x ) \boldsymbol{x}\cdot\boldsymbol{F}(s\boldsymbol{x})=x_jF_j(s\boldsymbol{x}) x ⋅ F ( s x ) = x j F j ( s x ) (j j j について和)を x l x_l x l で偏微分すると、積の微分法により
∂ ∂ x l [ x j F j ( s x ) ] = F l ( s x ) + s x j ( ∂ l F j ) ( s x ) . \frac{\partial}{\partial x_l}\bigl[x_jF_j(s\boldsymbol{x})\bigr]
= F_l(s\boldsymbol{x}) + s\,x_j\,(\partial_l F_j)(s\boldsymbol{x}). ∂ x l ∂ [ x j F j ( s x ) ] = F l ( s x ) + s x j ( ∂ l F j ) ( s x ) . 第 2 項は合成関数の微分(∂ ( s x ) l / ∂ x l = s \partial (s\boldsymbol{x})_l/\partial x_l = s ∂ ( s x ) l / ∂ x l = s )から出ます。ここで仮定 rot F = 0 \operatorname{rot}\boldsymbol{F}=\boldsymbol{0} rot F = 0 、すなわち ∂ l F j = ∂ j F l \partial_l F_j=\partial_j F_l ∂ l F j = ∂ j F l を使うと
∂ f ∂ x l ( x ) = ∫ 0 1 [ F l ( s x ) + s x j ( ∂ j F l ) ( s x ) ] d s . \frac{\partial f}{\partial x_l}(\boldsymbol{x})
=\int_0^1\Bigl[F_l(s\boldsymbol{x})+s\,x_j\,(\partial_j F_l)(s\boldsymbol{x})\Bigr]ds . ∂ x l ∂ f ( x ) = ∫ 0 1 [ F l ( s x ) + s x j ( ∂ j F l ) ( s x ) ] d s . 被積分関数は d d s [ s F l ( s x ) ] \dfrac{d}{ds}\bigl[s\,F_l(s\boldsymbol{x})\bigr] d s d [ s F l ( s x ) ] にほかなりません(実際、積の微分と連鎖律で d d s [ s F l ( s x ) ] = F l ( s x ) + s x j ( ∂ j F l ) ( s x ) \frac{d}{ds}[sF_l(s\boldsymbol{x})]=F_l(s\boldsymbol{x})+s\,x_j(\partial_jF_l)(s\boldsymbol{x}) d s d [ s F l ( s x )] = F l ( s x ) + s x j ( ∂ j F l ) ( s x ) )。したがって
∂ f ∂ x l ( x ) = [ s F l ( s x ) ] s = 0 s = 1 = F l ( x ) . \frac{\partial f}{\partial x_l}(\boldsymbol{x})=\Bigl[s\,F_l(s\boldsymbol{x})\Bigr]_{s=0}^{s=1}=F_l(\boldsymbol{x}). ∂ x l ∂ f ( x ) = [ s F l ( s x ) ] s = 0 s = 1 = F l ( x ) . l = 1 , 2 , 3 l=1,2,3 l = 1 , 2 , 3 で成り立つので grad f = F \operatorname{grad}f=\boldsymbol{F} grad f = F です。
∎ Lemma 2.2 (発散がゼロなら回転(ポアンカレの補題) )
Ω ⊂ R 3 \Omega\subset\mathbb{R}^3 Ω ⊂ R 3 を原点について星形な開集合とし、G : Ω → R 3 \boldsymbol{G}\colon\Omega\to\mathbb{R}^3 G : Ω → R 3 を C 1 C^1 C 1 級のベクトル場とする。Ω \Omega Ω 上で div G = 0 \operatorname{div}\boldsymbol{G}=0 div G = 0 ならば、C 2 C^2 C 2 級のベクトル場 A \boldsymbol{A} A で G = rot A \boldsymbol{G}=\operatorname{rot}\boldsymbol{A} G = rot A を満たすものが存在する。具体的には
A ( x ) = ∫ 0 1 s G ( s x ) × x d s \boldsymbol{A}(\boldsymbol{x})=\int_0^1 s\,\boldsymbol{G}(s\boldsymbol{x})\times\boldsymbol{x}\,ds A ( x ) = ∫ 0 1 s G ( s x ) × x d s がその 1 つである。
Proof(Lemma 2.2) 成分で書くと A i ( x ) = ∫ 0 1 s ε i j k G j ( s x ) x k d s A_i(\boldsymbol{x})=\int_0^1 s\,\varepsilon_{ijk}G_j(s\boldsymbol{x})x_k\,ds A i ( x ) = ∫ 0 1 s ε ij k G j ( s x ) x k d s です(ε i j k \varepsilon_{ijk} ε ij k はレヴィ=チヴィタ記号、繰り返し添字は和)。回転を取ります。
( rot A ) l = ε l m i ∂ m A i = ∫ 0 1 s ε l m i ε i j k ∂ m [ G j ( s x ) x k ] d s . (\operatorname{rot}\boldsymbol{A})_l=\varepsilon_{lmi}\partial_m A_i
=\int_0^1 s\,\varepsilon_{lmi}\varepsilon_{ijk}\,\partial_m\bigl[G_j(s\boldsymbol{x})x_k\bigr]ds . ( rot A ) l = ε l mi ∂ m A i = ∫ 0 1 s ε l mi ε ij k ∂ m [ G j ( s x ) x k ] d s . 括弧の中を積の微分法と連鎖律で展開すると
∂ m [ G j ( s x ) x k ] = s ( ∂ m G j ) ( s x ) x k + G j ( s x ) δ m k . \partial_m\bigl[G_j(s\boldsymbol{x})x_k\bigr]=s\,(\partial_mG_j)(s\boldsymbol{x})\,x_k+G_j(s\boldsymbol{x})\,\delta_{mk}. ∂ m [ G j ( s x ) x k ] = s ( ∂ m G j ) ( s x ) x k + G j ( s x ) δ mk . 一方 ε l m i ε i j k = ε i l m ε i j k = δ l j δ m k − δ l k δ m j \varepsilon_{lmi}\varepsilon_{ijk}=\varepsilon_{ilm}\varepsilon_{ijk}=\delta_{lj}\delta_{mk}-\delta_{lk}\delta_{mj} ε l mi ε ij k = ε i l m ε ij k = δ l j δ mk − δ l k δ mj です。これを代入して 2 つの項に分けます。
δ l j δ m k \delta_{lj}\delta_{mk} δ l j δ mk の項:s ( ∂ k G l ) ( s x ) x k + G l ( s x ) δ k k = s x ⋅ ( ∇ G l ) ( s x ) + 3 G l ( s x ) s\,(\partial_kG_l)(s\boldsymbol{x})x_k+G_l(s\boldsymbol{x})\delta_{kk}=s\,\boldsymbol{x}\cdot(\nabla G_l)(s\boldsymbol{x})+3G_l(s\boldsymbol{x}) s ( ∂ k G l ) ( s x ) x k + G l ( s x ) δ k k = s x ⋅ ( ∇ G l ) ( s x ) + 3 G l ( s x ) (δ k k = 3 \delta_{kk}=3 δ k k = 3 )。
δ l k δ m j \delta_{lk}\delta_{mj} δ l k δ mj の項:s ( ∂ j G j ) ( s x ) x l + G j ( s x ) δ j l = s ( div G ) ( s x ) x l + G l ( s x ) s\,(\partial_jG_j)(s\boldsymbol{x})x_l+G_j(s\boldsymbol{x})\delta_{jl}=s\,(\operatorname{div}\boldsymbol{G})(s\boldsymbol{x})\,x_l+G_l(s\boldsymbol{x}) s ( ∂ j G j ) ( s x ) x l + G j ( s x ) δ j l = s ( div G ) ( s x ) x l + G l ( s x ) 。ここで仮定 div G = 0 \operatorname{div}\boldsymbol{G}=0 div G = 0 を使うと第 1 項が消え、G l ( s x ) G_l(s\boldsymbol{x}) G l ( s x ) だけが残ります。
両者の差を取ると
( rot A ) l = ∫ 0 1 s [ s x ⋅ ( ∇ G l ) ( s x ) + 3 G l ( s x ) − G l ( s x ) ] d s = ∫ 0 1 [ s 2 x ⋅ ( ∇ G l ) ( s x ) + 2 s G l ( s x ) ] d s . (\operatorname{rot}\boldsymbol{A})_l=\int_0^1 s\Bigl[s\,\boldsymbol{x}\cdot(\nabla G_l)(s\boldsymbol{x})+3G_l(s\boldsymbol{x})-G_l(s\boldsymbol{x})\Bigr]ds
=\int_0^1\Bigl[s^2\,\boldsymbol{x}\cdot(\nabla G_l)(s\boldsymbol{x})+2s\,G_l(s\boldsymbol{x})\Bigr]ds . ( rot A ) l = ∫ 0 1 s [ s x ⋅ ( ∇ G l ) ( s x ) + 3 G l ( s x ) − G l ( s x ) ] d s = ∫ 0 1 [ s 2 x ⋅ ( ∇ G l ) ( s x ) + 2 s G l ( s x ) ] d s . 被積分関数は d d s [ s 2 G l ( s x ) ] \dfrac{d}{ds}\bigl[s^2G_l(s\boldsymbol{x})\bigr] d s d [ s 2 G l ( s x ) ] に一致します。よって
( rot A ) l = [ s 2 G l ( s x ) ] s = 0 s = 1 = G l ( x ) . (\operatorname{rot}\boldsymbol{A})_l=\Bigl[s^2G_l(s\boldsymbol{x})\Bigr]_{s=0}^{s=1}=G_l(\boldsymbol{x}). ( rot A ) l = [ s 2 G l ( s x ) ] s = 0 s = 1 = G l ( x ) . l = 1 , 2 , 3 l=1,2,3 l = 1 , 2 , 3 で成り立つので rot A = G \operatorname{rot}\boldsymbol{A}=\boldsymbol{G} rot A = G です。
∎ Theorem 3.1 (電磁ポテンシャルの存在 )
Ω ⊂ R 3 \Omega\subset\mathbb{R}^3 Ω ⊂ R 3 を星形な開集合、I ⊂ R I\subset\mathbb{R} I ⊂ R を開区間とする。Ω × I \Omega\times I Ω × I 上の C 2 C^2 C 2 級のベクトル場 E , B \boldsymbol{E},\boldsymbol{B} E , B が斉次マクスウェル方程式
div B = 0 , rot E = − ∂ B ∂ t \operatorname{div}\boldsymbol{B}=0,\qquad \operatorname{rot}\boldsymbol{E}=-\frac{\partial\boldsymbol{B}}{\partial t} div B = 0 , rot E = − ∂ t ∂ B を満たすとする。このとき Ω × I \Omega\times I Ω × I 上の C 1 C^1 C 1 級のスカラー場 φ \varphi φ とベクトル場 A \boldsymbol{A} A が存在して
B = rot A , E = − grad φ − ∂ A ∂ t \boldsymbol{B}=\operatorname{rot}\boldsymbol{A},\qquad
\boldsymbol{E}=-\operatorname{grad}\varphi-\frac{\partial\boldsymbol{A}}{\partial t} B = rot A , E = − grad φ − ∂ t ∂ A が Ω × I \Omega\times I Ω × I の全体で成り立つ。逆に、任意の C 2 C^2 C 2 級の ( φ , A ) (\varphi,\boldsymbol{A}) ( φ , A ) からこの 2 式で ( E , B ) (\boldsymbol{E},\boldsymbol{B}) ( E , B ) を定めれば、斉次マクスウェル方程式は恒等的に満たされる。
Proof(Theorem 3.1) (存在) 各時刻 t ∈ I t\in I t ∈ I を固定します。div B ( ⋅ , t ) = 0 \operatorname{div}\boldsymbol{B}(\cdot,t)=0 div B ( ⋅ , t ) = 0 なので、Lemma 2.2 により
A ( x , t ) = ∫ 0 1 s B ( s x , t ) × x d s \boldsymbol{A}(\boldsymbol{x},t)=\int_0^1 s\,\boldsymbol{B}(s\boldsymbol{x},t)\times\boldsymbol{x}\,ds A ( x , t ) = ∫ 0 1 s B ( s x , t ) × x d s とおけば rot A = B \operatorname{rot}\boldsymbol{A}=\boldsymbol{B} rot A = B です。この式は t t t を単なるパラメータとして含むだけなので、B \boldsymbol{B} B が t t t について C 2 C^2 C 2 級なら A \boldsymbol{A} A も t t t について微分できます。
次に E + ∂ t A \boldsymbol{E}+\partial_t\boldsymbol{A} E + ∂ t A を考えます。回転と時間微分は独立変数についての偏微分なので順序交換ができ、
rot ( E + ∂ A ∂ t ) = rot E + ∂ ∂ t ( rot A ) = − ∂ B ∂ t + ∂ B ∂ t = 0 . \operatorname{rot}\Bigl(\boldsymbol{E}+\frac{\partial\boldsymbol{A}}{\partial t}\Bigr)
=\operatorname{rot}\boldsymbol{E}+\frac{\partial}{\partial t}(\operatorname{rot}\boldsymbol{A})
=-\frac{\partial\boldsymbol{B}}{\partial t}+\frac{\partial\boldsymbol{B}}{\partial t}=\boldsymbol{0}. rot ( E + ∂ t ∂ A ) = rot E + ∂ t ∂ ( rot A ) = − ∂ t ∂ B + ∂ t ∂ B = 0 . 1 番目の等号で微分の順序交換、2 番目で仮定 rot E = − ∂ t B \operatorname{rot}\boldsymbol{E}=-\partial_t\boldsymbol{B} rot E = − ∂ t B と、いま作った rot A = B \operatorname{rot}\boldsymbol{A}=\boldsymbol{B} rot A = B を使いました。したがって各時刻で E + ∂ t A \boldsymbol{E}+\partial_t\boldsymbol{A} E + ∂ t A は回転がゼロなので、Lemma 2.1 により
φ ( x , t ) = − ∫ 0 1 x ⋅ ( E + ∂ A ∂ t ) ( s x , t ) d s \varphi(\boldsymbol{x},t)=-\int_0^1 \boldsymbol{x}\cdot\Bigl(\boldsymbol{E}+\frac{\partial\boldsymbol{A}}{\partial t}\Bigr)(s\boldsymbol{x},t)\,ds φ ( x , t ) = − ∫ 0 1 x ⋅ ( E + ∂ t ∂ A ) ( s x , t ) d s とおけば − grad φ = E + ∂ t A -\operatorname{grad}\varphi=\boldsymbol{E}+\partial_t\boldsymbol{A} − grad φ = E + ∂ t A 、すなわち E = − grad φ − ∂ t A \boldsymbol{E}=-\operatorname{grad}\varphi-\partial_t\boldsymbol{A} E = − grad φ − ∂ t A が得られます。
(逆) B = rot A \boldsymbol{B}=\operatorname{rot}\boldsymbol{A} B = rot A とすると、恒等式 div rot = 0 \operatorname{div}\operatorname{rot}=0 div rot = 0 から div B = 0 \operatorname{div}\boldsymbol{B}=0 div B = 0 。これは成分で書けば ∂ i ε i j k ∂ j A k = ε i j k ∂ i ∂ j A k \partial_i\varepsilon_{ijk}\partial_jA_k=\varepsilon_{ijk}\partial_i\partial_jA_k ∂ i ε ij k ∂ j A k = ε ij k ∂ i ∂ j A k で、ε i j k \varepsilon_{ijk} ε ij k が i , j i,j i , j の交換で符号を変え ∂ i ∂ j \partial_i\partial_j ∂ i ∂ j が変えないため和がゼロになる、ということです。また E = − grad φ − ∂ t A \boldsymbol{E}=-\operatorname{grad}\varphi-\partial_t\boldsymbol{A} E = − grad φ − ∂ t A の回転を取ると、恒等式 rot grad = 0 \operatorname{rot}\operatorname{grad}=\boldsymbol{0} rot grad = 0 (同じ理由で従います)により
rot E = − rot grad φ − ∂ ∂ t rot A = − ∂ B ∂ t \operatorname{rot}\boldsymbol{E}=-\operatorname{rot}\operatorname{grad}\varphi-\frac{\partial}{\partial t}\operatorname{rot}\boldsymbol{A}
=-\frac{\partial\boldsymbol{B}}{\partial t} rot E = − rot grad φ − ∂ t ∂ rot A = − ∂ t ∂ B となります。
∎ Definition 3.2 (電磁ポテンシャル )
Theorem 3.1 の A \boldsymbol{A} A をベクトルポテンシャル 、φ \varphi φ をスカラーポテンシャル と呼び、あわせて電磁ポテンシャル という。すなわち
B = rot A , E = − grad φ − ∂ A ∂ t \boldsymbol{B}=\operatorname{rot}\boldsymbol{A},\qquad
\boldsymbol{E}=-\operatorname{grad}\varphi-\frac{\partial\boldsymbol{A}}{\partial t} B = rot A , E = − grad φ − ∂ t ∂ A の関係にある ( φ , A ) (\varphi,\boldsymbol{A}) ( φ , A ) のことである。
flowchart TD
A["div B = 0"] --> B["B = rot A(ベクトルポテンシャル)"]
C["rot E = -∂B/∂t"] --> D["rot(E + ∂A/∂t)= 0"]
B --> D
D --> E["E = -grad φ - ∂A/∂t(スカラーポテンシャル)"] 斉次マクスウェル方程式からポテンシャルが出てくる筋道 Example 3.5 (ソレノイドの外側:磁場はゼロでもポテンシャルは消せない )
z z z 軸を中心軸とする半径 a a a の無限に長いソレノイドを考えます。定常電流と静磁場 で見たとおり(無限に長いソレノイド(Example 6.6)[Steady Currents and Magnetostatic Fields] )、内部は一様磁場 B = B e z \boldsymbol{B}=B\boldsymbol{e}_z B = B e z 、外部は B = 0 \boldsymbol{B}=\boldsymbol{0} B = 0 です。
円筒対称性から A = A ϕ ( r ) e ϕ \boldsymbol{A}=A_\phi(r)\boldsymbol{e}_\phi A = A ϕ ( r ) e ϕ の形を仮定します(ここでの ϕ \phi ϕ は方位角で、スカラーポテンシャル φ \varphi φ とは別の記号です)。半径 r r r の円周にストークスの定理(Theorem 5.3)[Steady Currents and Magnetostatic Fields] を適用すると
∮ r A ⋅ d l = 2 π r A ϕ ( r ) = ∫ B ⋅ d S = Φ ( r ) \oint_{r}\boldsymbol{A}\cdot d\boldsymbol{l}=2\pi r\,A_\phi(r)=\int \boldsymbol{B}\cdot d\boldsymbol{S}=\Phi(r) ∮ r A ⋅ d l = 2 π r A ϕ ( r ) = ∫ B ⋅ d S = Φ ( r ) なので、Φ ( r ) \Phi(r) Φ ( r ) をその円が囲む磁束として
A ϕ ( r ) = Φ ( r ) 2 π r = { B r 2 , r ≤ a , Φ 2 π r , r ≥ a , Φ = π a 2 B . A_\phi(r)=\frac{\Phi(r)}{2\pi r}=
\begin{cases}
\dfrac{Br}{2}, & r\le a,\\[8pt]
\dfrac{\Phi}{2\pi r}, & r\ge a,
\end{cases}
\qquad \Phi=\pi a^2B . A ϕ ( r ) = 2 π r Φ ( r ) = ⎩ ⎨ ⎧ 2 B r , 2 π r Φ , r ≤ a , r ≥ a , Φ = π a 2 B . 外側では B = 0 \boldsymbol{B}=\boldsymbol{0} B = 0 なのに A ≠ 0 \boldsymbol{A}\neq\boldsymbol{0} A = 0 です。もちろん外側で rot A = 0 \operatorname{rot}\boldsymbol{A}=\boldsymbol{0} rot A = 0 は成り立っています(円筒座標で ( rot A ) z = 1 r ∂ ∂ r ( r A ϕ ) = 1 r ∂ ∂ r ( Φ 2 π ) = 0 (\operatorname{rot}\boldsymbol{A})_z=\frac{1}{r}\frac{\partial}{\partial r}(rA_\phi)=\frac{1}{r}\frac{\partial}{\partial r}\bigl(\frac{\Phi}{2\pi}\bigr)=0 ( rot A ) z = r 1 ∂ r ∂ ( r A ϕ ) = r 1 ∂ r ∂ ( 2 π Φ ) = 0 )。しかし外部領域で A = grad χ \boldsymbol{A}=\operatorname{grad}\chi A = grad χ と書ける 1 価関数 χ \chi χ は存在しません。もし存在すれば閉曲線に沿った周回積分が 0 0 0 になるはずですが、実際には ∮ A ⋅ d l = Φ ≠ 0 \oint\boldsymbol{A}\cdot d\boldsymbol{l}=\Phi\neq0 ∮ A ⋅ d l = Φ = 0 だからです。
外部領域は z z z 軸のまわりに穴が開いていて星形ではなく、Lemma 2.1 の仮定が破れています。ポテンシャルの理論が領域の位相に敏感であることの、最も分かりやすい例です。この ∮ A ⋅ d l \oint\boldsymbol{A}\cdot d\boldsymbol{l} ∮ A ⋅ d l が量子力学で電子の干渉縞を動かす(アハラノフ=ボーム効果)ことが、後にポテンシャルの地位を決定づけました。
符号の組み合わせが重要です。A \boldsymbol{A} A には + grad χ +\operatorname{grad}\chi + grad χ 、φ \varphi φ には − ∂ t χ -\partial_t\chi − ∂ t χ 。この組でなければ E \boldsymbol{E} E が変わってしまいます。
Proposition 4.2 (場のゲージ不変性 )
Definition 4.1 のゲージ変換によって、( φ , A ) (\varphi,\boldsymbol{A}) ( φ , A ) から定まる場は変わらない。すなわち
rot A ′ = rot A , − grad φ ′ − ∂ A ′ ∂ t = − grad φ − ∂ A ∂ t . \operatorname{rot}\boldsymbol{A}'=\operatorname{rot}\boldsymbol{A},\qquad
-\operatorname{grad}\varphi'-\frac{\partial\boldsymbol{A}'}{\partial t}=-\operatorname{grad}\varphi-\frac{\partial\boldsymbol{A}}{\partial t}. rot A ′ = rot A , − grad φ ′ − ∂ t ∂ A ′ = − grad φ − ∂ t ∂ A . Proof(Proposition 4.2) 磁場については、恒等式 rot grad χ = 0 \operatorname{rot}\operatorname{grad}\chi=\boldsymbol{0} rot grad χ = 0 から
rot A ′ = rot A + rot grad χ = rot A . \operatorname{rot}\boldsymbol{A}'=\operatorname{rot}\boldsymbol{A}+\operatorname{rot}\operatorname{grad}\chi=\operatorname{rot}\boldsymbol{A}. rot A ′ = rot A + rot grad χ = rot A . 電場については、χ \chi χ が C 2 C^2 C 2 級なので grad \operatorname{grad} grad と ∂ t \partial_t ∂ t が交換すること(シュワルツの定理(Theorem 7.1)[多変数関数の微分と偏微分] )を使って
− grad φ ′ − ∂ A ′ ∂ t = − grad φ + grad ∂ χ ∂ t − ∂ A ∂ t − ∂ ∂ t grad χ = − grad φ − ∂ A ∂ t -\operatorname{grad}\varphi'-\frac{\partial\boldsymbol{A}'}{\partial t}
=-\operatorname{grad}\varphi+\operatorname{grad}\frac{\partial\chi}{\partial t}
-\frac{\partial\boldsymbol{A}}{\partial t}-\frac{\partial}{\partial t}\operatorname{grad}\chi
=-\operatorname{grad}\varphi-\frac{\partial\boldsymbol{A}}{\partial t} − grad φ ′ − ∂ t ∂ A ′ = − grad φ + grad ∂ t ∂ χ − ∂ t ∂ A − ∂ t ∂ grad χ = − grad φ − ∂ t ∂ A となり、grad ∂ t χ \operatorname{grad}\partial_t\chi grad ∂ t χ と ∂ t grad χ \partial_t\operatorname{grad}\chi ∂ t grad χ が打ち消し合います。
∎ Proposition 4.2 は「ゲージ変換をしても場は変わらない」と言っています。逆に「場が同じなら必ずゲージ変換で移り合う」でしょうか。これが成り立つからこそ、不定性を完全に把握したと言えます。
Theorem 4.3 (ゲージ自由度の完全性 )
Ω ⊂ R 3 \Omega\subset\mathbb{R}^3 Ω ⊂ R 3 を星形な開集合、I I I を開区間とする。Ω × I \Omega\times I Ω × I 上の 2 組の C 2 C^2 C 2 級ポテンシャル ( φ , A ) (\varphi,\boldsymbol{A}) ( φ , A ) と ( φ ′ , A ′ ) (\varphi',\boldsymbol{A}') ( φ ′ , A ′ ) が同一の場を与える、すなわち
rot A ′ = rot A , − grad φ ′ − ∂ A ′ ∂ t = − grad φ − ∂ A ∂ t \operatorname{rot}\boldsymbol{A}'=\operatorname{rot}\boldsymbol{A},\qquad
-\operatorname{grad}\varphi'-\frac{\partial\boldsymbol{A}'}{\partial t}=-\operatorname{grad}\varphi-\frac{\partial\boldsymbol{A}}{\partial t} rot A ′ = rot A , − grad φ ′ − ∂ t ∂ A ′ = − grad φ − ∂ t ∂ A が成り立つとする。このとき C 2 C^2 C 2 級のスカラー場 χ ( x , t ) \chi(\boldsymbol{x},t) χ ( x , t ) が存在して
A ′ = A + grad χ , φ ′ = φ − ∂ χ ∂ t \boldsymbol{A}'=\boldsymbol{A}+\operatorname{grad}\chi,\qquad \varphi'=\varphi-\frac{\partial\chi}{\partial t} A ′ = A + grad χ , φ ′ = φ − ∂ t ∂ χ が成り立つ。すなわち 2 組はゲージ変換で結ばれる。
Proof(Theorem 4.3) D = A ′ − A \boldsymbol{D}=\boldsymbol{A}'-\boldsymbol{A} D = A ′ − A 、u = φ ′ − φ u=\varphi'-\varphi u = φ ′ − φ とおきます。仮定の第 1 式から rot D = 0 \operatorname{rot}\boldsymbol{D}=\boldsymbol{0} rot D = 0 です。各時刻 t t t で Lemma 2.1 を適用すると、Ω \Omega Ω が星形なので
χ ( x , t ) = ∫ 0 1 x ⋅ D ( s x , t ) d s \chi(\boldsymbol{x},t)=\int_0^1\boldsymbol{x}\cdot\boldsymbol{D}(s\boldsymbol{x},t)\,ds χ ( x , t ) = ∫ 0 1 x ⋅ D ( s x , t ) d s が D = grad χ \boldsymbol{D}=\operatorname{grad}\chi D = grad χ を与えます。t t t は積分の中でパラメータとしてしか現れないので、χ \chi χ は t t t について微分可能です。
次に仮定の第 2 式に A ′ = A + grad χ \boldsymbol{A}'=\boldsymbol{A}+\operatorname{grad}\chi A ′ = A + grad χ を代入します。
− grad φ ′ − ∂ A ∂ t − ∂ ∂ t grad χ = − grad φ − ∂ A ∂ t -\operatorname{grad}\varphi'-\frac{\partial\boldsymbol{A}}{\partial t}-\frac{\partial}{\partial t}\operatorname{grad}\chi
=-\operatorname{grad}\varphi-\frac{\partial\boldsymbol{A}}{\partial t} − grad φ ′ − ∂ t ∂ A − ∂ t ∂ grad χ = − grad φ − ∂ t ∂ A より grad ( φ ′ − φ + ∂ t χ ) = 0 \operatorname{grad}\bigl(\varphi'-\varphi+\partial_t\chi\bigr)=\boldsymbol{0} grad ( φ ′ − φ + ∂ t χ ) = 0 。Ω \Omega Ω は星形なので連結であり、勾配がゼロの関数は x \boldsymbol{x} x に依りません。よってある関数 g ( t ) g(t) g ( t ) があって
φ ′ = φ − ∂ χ ∂ t + g ( t ) . \varphi'=\varphi-\frac{\partial\chi}{\partial t}+g(t). φ ′ = φ − ∂ t ∂ χ + g ( t ) . 最後に g g g を吸収します。h ( t ) h(t) h ( t ) を h ′ ( t ) = g ( t ) h'(t)=g(t) h ′ ( t ) = g ( t ) を満たす関数(g g g の原始関数、積分の基本定理 により連続な g g g に対して存在します)とし、χ ~ ( x , t ) = χ ( x , t ) − h ( t ) \tilde\chi(\boldsymbol{x},t)=\chi(\boldsymbol{x},t)-h(t) χ ~ ( x , t ) = χ ( x , t ) − h ( t ) とおきます。h h h は x \boldsymbol{x} x に依らないので grad χ ~ = grad χ = D \operatorname{grad}\tilde\chi=\operatorname{grad}\chi=\boldsymbol{D} grad χ ~ = grad χ = D は変わらず、一方
φ − ∂ χ ~ ∂ t = φ − ∂ χ ∂ t + h ′ ( t ) = φ − ∂ χ ∂ t + g ( t ) = φ ′ . \varphi-\frac{\partial\tilde\chi}{\partial t}=\varphi-\frac{\partial\chi}{\partial t}+h'(t)=\varphi-\frac{\partial\chi}{\partial t}+g(t)=\varphi' . φ − ∂ t ∂ χ ~ = φ − ∂ t ∂ χ + h ′ ( t ) = φ − ∂ t ∂ χ + g ( t ) = φ ′ . よって χ ~ \tilde\chi χ ~ が求めるゲージ関数です。
∎ Proposition 4.2 と Theorem 4.3 をあわせると、次のように言えます。ポテンシャルの全体は、同じ物理的な場を与えるもの同士でグループ分けされ、各グループ(ゲージ軌道)はちょうど 1 つのゲージ関数 χ \chi χ の分だけの自由度を持つ。 物理的に意味があるのはグループそのものであって、その中のどの代表元を選ぶかは私たちの都合です。
ポテンシャルの空間 (φ, A) 物理的な場の空間 (E, B) ゲージ軌道 場を計算する写像 ポテンシャルの空間はゲージ軌道に分かれ、各軌道が 1 つの物理的な場に対応する Example 4.4 (スカラーポテンシャルはいつでも消せる )
φ \varphi φ が物理量でないことを最もはっきり示す例です。任意の ( φ , A ) (\varphi,\boldsymbol{A}) ( φ , A ) に対して、時刻 t 0 t_0 t 0 を固定して
χ ( x , t ) = ∫ t 0 t φ ( x , t ′ ) d t ′ \chi(\boldsymbol{x},t)=\int_{t_0}^{t}\varphi(\boldsymbol{x},t')\,dt' χ ( x , t ) = ∫ t 0 t φ ( x , t ′ ) d t ′ とおきます。微分積分学の基本定理 (第一部(Theorem 5.2)[積分の基本定理と定積分] )から ∂ t χ = φ \partial_t\chi=\varphi ∂ t χ = φ なので、このゲージ変換の結果は
φ ′ = φ − ∂ χ ∂ t = 0 , A ′ = A + grad ∫ t 0 t φ d t ′ . \varphi'=\varphi-\frac{\partial\chi}{\partial t}=0,\qquad
\boldsymbol{A}'=\boldsymbol{A}+\operatorname{grad}\int_{t_0}^{t}\varphi\,dt' . φ ′ = φ − ∂ t ∂ χ = 0 , A ′ = A + grad ∫ t 0 t φ d t ′ . つまり、どんな電磁場でも φ ≡ 0 \varphi\equiv0 φ ≡ 0 となるポテンシャルの取り方が存在します(時間ゲージ 、あるいはワイルゲージ)。このとき電場は E = − ∂ t A ′ \boldsymbol{E}=-\partial_t\boldsymbol{A}' E = − ∂ t A ′ だけで表されます。静電気の問題ですら、φ \varphi φ を消して A \boldsymbol{A} A が時間とともに変化する記述に書き換えられるのです。「電位」に絶対的な意味がないことが、ここに極端な形で現れています。
Theorem 3.1 により斉次方程式は自動的に満たされたので、残る非斉次方程式をポテンシャルで書き直します。
ガウスの法則 div E = ρ / ε 0 \operatorname{div}\boldsymbol{E}=\rho/\varepsilon_0 div E = ρ / ε 0 に E = − grad φ − ∂ t A \boldsymbol{E}=-\operatorname{grad}\varphi-\partial_t\boldsymbol{A} E = − grad φ − ∂ t A を代入して
− ∇ 2 φ − ∂ ∂ t ( div A ) = ρ ε 0 . -\nabla^2\varphi-\frac{\partial}{\partial t}(\operatorname{div}\boldsymbol{A})=\frac{\rho}{\varepsilon_0}. − ∇ 2 φ − ∂ t ∂ ( div A ) = ε 0 ρ . アンペール=マクスウェルの法則には恒等式 rot rot A = grad ( div A ) − ∇ 2 A \operatorname{rot}\operatorname{rot}\boldsymbol{A}=\operatorname{grad}(\operatorname{div}\boldsymbol{A})-\nabla^2\boldsymbol{A} rot rot A = grad ( div A ) − ∇ 2 A を使います。
grad ( div A ) − ∇ 2 A = μ 0 j + 1 c 2 ∂ ∂ t ( − grad φ − ∂ A ∂ t ) \operatorname{grad}(\operatorname{div}\boldsymbol{A})-\nabla^2\boldsymbol{A}
=\mu_0\boldsymbol{j}+\frac{1}{c^2}\frac{\partial}{\partial t}\Bigl(-\operatorname{grad}\varphi-\frac{\partial\boldsymbol{A}}{\partial t}\Bigr) grad ( div A ) − ∇ 2 A = μ 0 j + c 2 1 ∂ t ∂ ( − grad φ − ∂ t ∂ A ) を整理すると
( ∇ 2 − 1 c 2 ∂ 2 ∂ t 2 ) A − grad ( div A + 1 c 2 ∂ φ ∂ t ) = − μ 0 j \Bigl(\nabla^2-\frac{1}{c^2}\frac{\partial^2}{\partial t^2}\Bigr)\boldsymbol{A}
-\operatorname{grad}\Bigl(\operatorname{div}\boldsymbol{A}+\frac{1}{c^2}\frac{\partial\varphi}{\partial t}\Bigr)=-\mu_0\boldsymbol{j} ( ∇ 2 − c 2 1 ∂ t 2 ∂ 2 ) A − grad ( div A + c 2 1 ∂ t ∂ φ ) = − μ 0 j となります。以下、ダランベール演算子 を
□ ≡ ∇ 2 − 1 c 2 ∂ 2 ∂ t 2 \Box\;\equiv\;\nabla^2-\frac{1}{c^2}\frac{\partial^2}{\partial t^2} □ ≡ ∇ 2 − c 2 1 ∂ t 2 ∂ 2 と書きます。2 本の方程式は φ \varphi φ と A \boldsymbol{A} A が絡み合っていて、このままでは扱いにくい。しかし両者に共通して現れる組み合わせ div A + c − 2 ∂ t φ \operatorname{div}\boldsymbol{A}+c^{-2}\partial_t\varphi div A + c − 2 ∂ t φ に目を付けると、道が開けます。Theorem 4.3 が保証するゲージ自由度は、まさにこの組み合わせを操作するために使えるのです。
Definition 5.1 (ローレンツゲージ条件 )
電磁ポテンシャル ( φ , A ) (\varphi,\boldsymbol{A}) ( φ , A ) が
div A + 1 c 2 ∂ φ ∂ t = 0 \operatorname{div}\boldsymbol{A}+\frac{1}{c^2}\frac{\partial\varphi}{\partial t}=0 div A + c 2 1 ∂ t ∂ φ = 0 を満たすとき、ローレンツゲージ にあるという。
Theorem 5.2 (ローレンツゲージはいつでも実現できる )
( φ , A ) (\varphi,\boldsymbol{A}) ( φ , A ) を任意の電磁ポテンシャルとする。非斉次波動方程式
□ χ = − ( div A + 1 c 2 ∂ φ ∂ t ) \Box\chi=-\Bigl(\operatorname{div}\boldsymbol{A}+\frac{1}{c^2}\frac{\partial\varphi}{\partial t}\Bigr) □ χ = − ( div A + c 2 1 ∂ t ∂ φ ) の解 χ \chi χ を用いてゲージ変換すれば、得られる ( φ ′ , A ′ ) (\varphi',\boldsymbol{A}') ( φ ′ , A ′ ) はローレンツゲージ条件を満たす。さらに、ローレンツゲージにあるポテンシャルどうしを結ぶゲージ変換は、□ χ = 0 \Box\chi=0 □ χ = 0 を満たす χ \chi χ によるものに限る。
Proof(Theorem 5.2) ゲージ変換 A ′ = A + grad χ \boldsymbol{A}'=\boldsymbol{A}+\operatorname{grad}\chi A ′ = A + grad χ 、φ ′ = φ − ∂ t χ \varphi'=\varphi-\partial_t\chi φ ′ = φ − ∂ t χ のもとで、ゲージ条件の左辺がどう変わるかを計算します。div grad = ∇ 2 \operatorname{div}\operatorname{grad}=\nabla^2 div grad = ∇ 2 に注意して
div A ′ + 1 c 2 ∂ φ ′ ∂ t = div A + ∇ 2 χ + 1 c 2 ∂ φ ∂ t − 1 c 2 ∂ 2 χ ∂ t 2 = ( div A + 1 c 2 ∂ φ ∂ t ) + □ χ . \operatorname{div}\boldsymbol{A}'+\frac{1}{c^2}\frac{\partial\varphi'}{\partial t}
=\operatorname{div}\boldsymbol{A}+\nabla^2\chi+\frac{1}{c^2}\frac{\partial\varphi}{\partial t}-\frac{1}{c^2}\frac{\partial^2\chi}{\partial t^2}
=\Bigl(\operatorname{div}\boldsymbol{A}+\frac{1}{c^2}\frac{\partial\varphi}{\partial t}\Bigr)+\Box\chi . div A ′ + c 2 1 ∂ t ∂ φ ′ = div A + ∇ 2 χ + c 2 1 ∂ t ∂ φ − c 2 1 ∂ t 2 ∂ 2 χ = ( div A + c 2 1 ∂ t ∂ φ ) + □ χ . したがって左辺をゼロにするには □ χ \Box\chi □ χ が元のゲージ条件の左辺の符号を変えたものに等しければよく、それが定理の主張する方程式です。この非斉次波動方程式は、右辺が与えられた関数であるとき(適切な滑らかさと減衰のもとで)遅延解
χ ( x , t ) = 1 4 π ∫ S ( x ′ , t − ∣ x − x ′ ∣ / c ) ∣ x − x ′ ∣ d 3 x ′ ( S = div A + c − 2 ∂ t φ ) \chi(\boldsymbol{x},t)=\frac{1}{4\pi}\int \frac{S\bigl(\boldsymbol{x}',\,t-|\boldsymbol{x}-\boldsymbol{x}'|/c\bigr)}{|\boldsymbol{x}-\boldsymbol{x}'|}\,d^3x'
\quad\bigl(S=\operatorname{div}\boldsymbol{A}+c^{-2}\partial_t\varphi\bigr) χ ( x , t ) = 4 π 1 ∫ ∣ x − x ′ ∣ S ( x ′ , t − ∣ x − x ′ ∣/ c ) d 3 x ′ ( S = div A + c − 2 ∂ t φ ) を持ちます。少なくとも 1 つ解が存在すればよいので、これで実現可能性が示せました。
後半も上の計算から直ちに従います。変換前後がどちらもローレンツゲージ条件を満たすなら、上の等式の両端がゼロなので □ χ = 0 \Box\chi=0 □ χ = 0 です。逆に □ χ = 0 \Box\chi=0 □ χ = 0 ならゲージ条件は保たれます。
∎ Corollary 5.3 (ローレンツゲージにおける波動方程式 )
ローレンツゲージ条件のもとで、マクスウェル方程式は
□ φ = − ρ ε 0 , □ A = − μ 0 j \Box\varphi=-\frac{\rho}{\varepsilon_0},\qquad
\Box\boldsymbol{A}=-\mu_0\boldsymbol{j} □ φ = − ε 0 ρ , □ A = − μ 0 j と同値である。すなわち φ \varphi φ と A \boldsymbol{A} A の 3 成分は、それぞれ独立な非斉次波動方程式に従う。
Proof(Corollary 5.3) §5.1 で得た 2 本の式に条件 Definition 5.1 を代入します。
第 1 式について。div A = − c − 2 ∂ t φ \operatorname{div}\boldsymbol{A}=-c^{-2}\partial_t\varphi div A = − c − 2 ∂ t φ を − ∇ 2 φ − ∂ t ( div A ) = ρ / ε 0 -\nabla^2\varphi-\partial_t(\operatorname{div}\boldsymbol{A})=\rho/\varepsilon_0 − ∇ 2 φ − ∂ t ( div A ) = ρ / ε 0 に入れると
− ∇ 2 φ + 1 c 2 ∂ 2 φ ∂ t 2 = ρ ε 0 ⟺ □ φ = − ρ ε 0 . -\nabla^2\varphi+\frac{1}{c^2}\frac{\partial^2\varphi}{\partial t^2}=\frac{\rho}{\varepsilon_0}
\quad\Longleftrightarrow\quad
\Box\varphi=-\frac{\rho}{\varepsilon_0}. − ∇ 2 φ + c 2 1 ∂ t 2 ∂ 2 φ = ε 0 ρ ⟺ □ φ = − ε 0 ρ . 第 2 式について。□ A − grad ( div A + c − 2 ∂ t φ ) = − μ 0 j \Box\boldsymbol{A}-\operatorname{grad}(\operatorname{div}\boldsymbol{A}+c^{-2}\partial_t\varphi)=-\mu_0\boldsymbol{j} □ A − grad ( div A + c − 2 ∂ t φ ) = − μ 0 j の第 2 項の括弧がゲージ条件によりゼロなので、□ A = − μ 0 j \Box\boldsymbol{A}=-\mu_0\boldsymbol{j} □ A = − μ 0 j です。
逆向きも同じ代入を逆にたどるだけです。
∎ これが電磁ポテンシャルを導入する実用上の最大の見返りです。6 成分・8 本の連立方程式が、同じ形をした 4 本の独立なスカラー波動方程式に化けました。しかもこの形から、ρ \rho ρ と j \boldsymbol{j} j が作る電磁場が速さ c c c で伝わることが即座に読み取れます。
Definition 5.5 (クーロンゲージ条件 )
電磁ポテンシャル ( φ , A ) (\varphi,\boldsymbol{A}) ( φ , A ) が
div A = 0 \operatorname{div}\boldsymbol{A}=0 div A = 0 を満たすとき、クーロンゲージ (横波ゲージ、放射ゲージ)にあるという。
Proposition 5.6 (クーロンゲージの実現と、そこでのスカラーポテンシャル )
任意の電磁ポテンシャル ( φ , A ) (\varphi,\boldsymbol{A}) ( φ , A ) に対し、ポアソン方程式 ∇ 2 χ = − div A \nabla^2\chi=-\operatorname{div}\boldsymbol{A} ∇ 2 χ = − div A の解 χ \chi χ によるゲージ変換は ( φ ′ , A ′ ) (\varphi',\boldsymbol{A}') ( φ ′ , A ′ ) をクーロンゲージに移す。さらにクーロンゲージのもとでスカラーポテンシャルは
∇ 2 φ = − ρ ε 0 \nabla^2\varphi=-\frac{\rho}{\varepsilon_0} ∇ 2 φ = − ε 0 ρ を満たし、全空間で φ → 0 \varphi\to0 φ → 0 (∣ x ∣ → ∞ |\boldsymbol{x}|\to\infty ∣ x ∣ → ∞ )を課せば
φ ( x , t ) = 1 4 π ε 0 ∫ ρ ( x ′ , t ) ∣ x − x ′ ∣ d 3 x ′ \varphi(\boldsymbol{x},t)=\frac{1}{4\pi\varepsilon_0}\int\frac{\rho(\boldsymbol{x}',t)}{|\boldsymbol{x}-\boldsymbol{x}'|}\,d^3x' φ ( x , t ) = 4 π ε 0 1 ∫ ∣ x − x ′ ∣ ρ ( x ′ , t ) d 3 x ′ と、同時刻 の電荷分布で書ける。
Proof(Proposition 5.6) ゲージ変換で div A ′ = div A + ∇ 2 χ \operatorname{div}\boldsymbol{A}'=\operatorname{div}\boldsymbol{A}+\nabla^2\chi div A ′ = div A + ∇ 2 χ ですから、これがゼロになる条件が ∇ 2 χ = − div A \nabla^2\chi=-\operatorname{div}\boldsymbol{A} ∇ 2 χ = − div A です。ポアソン方程式は(適切な条件下で)常に解を持ちます。
次に §5.1 の第 1 式 − ∇ 2 φ − ∂ t ( div A ) = ρ / ε 0 -\nabla^2\varphi-\partial_t(\operatorname{div}\boldsymbol{A})=\rho/\varepsilon_0 − ∇ 2 φ − ∂ t ( div A ) = ρ / ε 0 に div A = 0 \operatorname{div}\boldsymbol{A}=0 div A = 0 を代入すると ∇ 2 φ = − ρ / ε 0 \nabla^2\varphi=-\rho/\varepsilon_0 ∇ 2 φ = − ρ / ε 0 となります。この形は静電場とガウスの法則 で扱った静電ポテンシャルの方程式(ポアソン方程式(Theorem 6.5)[Electrostatic Fields and Gauss's Law] )とまったく同じで、時間 t t t はパラメータとしてしか入りません。境界条件 φ → 0 \varphi\to0 φ → 0 のもとでの解が上のクーロン積分です。
∎
ローレンツゲージ クーロンゲージ 条件 div A + c − 2 ∂ t φ = 0 \operatorname{div}\boldsymbol{A}+c^{-2}\partial_t\varphi=0 div A + c − 2 ∂ t φ = 0 div A = 0 \operatorname{div}\boldsymbol{A}=0 div A = 0 φ \varphi φ の方程式□ φ = − ρ / ε 0 \Box\varphi=-\rho/\varepsilon_0 □ φ = − ρ / ε 0 (波動)∇ 2 φ = − ρ / ε 0 \nabla^2\varphi=-\rho/\varepsilon_0 ∇ 2 φ = − ρ / ε 0 (瞬間的)A \boldsymbol{A} A の方程式□ A = − μ 0 j \Box\boldsymbol{A}=-\mu_0\boldsymbol{j} □ A = − μ 0 j □ A = − μ 0 j T \Box\boldsymbol{A}=-\mu_0\boldsymbol{j}_{\mathrm{T}} □ A = − μ 0 j T (Appendix 参照)ローレンツ変換 形を変えない 形を変える(慣性系ごとに条件が違う) 得意な場面 電磁波の放射、相対論的な計算 束縛電子系、量子光学、非相対論的な近似
Example 5.7 (静止した点電荷 )
原点に静止した点電荷 q q q を考えます。ρ ( x ) = q δ 3 ( x ) \rho(\boldsymbol{x})=q\delta^3(\boldsymbol{x}) ρ ( x ) = q δ 3 ( x ) 、j = 0 \boldsymbol{j}=\boldsymbol{0} j = 0 です。Proposition 5.6 の積分公式から
φ ( x ) = 1 4 π ε 0 ∫ q δ 3 ( x ′ ) ∣ x − x ′ ∣ d 3 x ′ = q 4 π ε 0 ∣ x ∣ , A = 0 \varphi(\boldsymbol{x})=\frac{1}{4\pi\varepsilon_0}\int\frac{q\delta^3(\boldsymbol{x}')}{|\boldsymbol{x}-\boldsymbol{x}'|}d^3x'=\frac{q}{4\pi\varepsilon_0|\boldsymbol{x}|},
\qquad \boldsymbol{A}=\boldsymbol{0} φ ( x ) = 4 π ε 0 1 ∫ ∣ x − x ′ ∣ q δ 3 ( x ′ ) d 3 x ′ = 4 π ε 0 ∣ x ∣ q , A = 0 です(A \boldsymbol{A} A については □ A = 0 \Box\boldsymbol{A}=\boldsymbol{0} □ A = 0 の解で無限遠でゼロかつ時間に依らないものを取りました)。この組はクーロンゲージ条件 div A = 0 \operatorname{div}\boldsymbol{A}=0 div A = 0 を満たし、同時に φ \varphi φ が時間に依らないのでローレンツゲージ条件 div A + c − 2 ∂ t φ = 0 \operatorname{div}\boldsymbol{A}+c^{-2}\partial_t\varphi=0 div A + c − 2 ∂ t φ = 0 も満たします。静的な問題では 2 つのゲージが一致するわけです。得られる場は
E = − grad φ = q 4 π ε 0 x ∣ x ∣ 3 , B = 0 \boldsymbol{E}=-\operatorname{grad}\varphi=\frac{q}{4\pi\varepsilon_0}\frac{\boldsymbol{x}}{|\boldsymbol{x}|^3},\qquad \boldsymbol{B}=\boldsymbol{0} E = − grad φ = 4 π ε 0 q ∣ x ∣ 3 x , B = 0 で、クーロンの法則が再現されました。
Example 5.8 (真空中の平面電磁波 )
ρ = 0 \rho=0 ρ = 0 、j = 0 \boldsymbol{j}=\boldsymbol{0} j = 0 の真空を考えます。Corollary 5.3 により □ φ = 0 \Box\varphi=0 □ φ = 0 、□ A = 0 \Box\boldsymbol{A}=\boldsymbol{0} □ A = 0 です。Example 4.4 と Theorem 5.2 の残留自由度を使うと、真空では φ ≡ 0 \varphi\equiv0 φ ≡ 0 かつ div A = 0 \operatorname{div}\boldsymbol{A}=0 div A = 0 を同時に課せます(放射ゲージ)。そこで
A ( x , t ) = A 0 ϵ cos ( k ⋅ x − ω t ) , ϵ ⋅ k = 0 , ∣ ϵ ∣ = 1 \boldsymbol{A}(\boldsymbol{x},t)=A_0\,\boldsymbol{\epsilon}\cos(\boldsymbol{k}\cdot\boldsymbol{x}-\omega t),\qquad \boldsymbol{\epsilon}\cdot\boldsymbol{k}=0,\ |\boldsymbol{\epsilon}|=1 A ( x , t ) = A 0 ϵ cos ( k ⋅ x − ω t ) , ϵ ⋅ k = 0 , ∣ ϵ ∣ = 1 とおきます。div A = − A 0 ( ϵ ⋅ k ) sin ( ⋯ ) = 0 \operatorname{div}\boldsymbol{A}=-A_0(\boldsymbol{\epsilon}\cdot\boldsymbol{k})\sin(\cdots)=0 div A = − A 0 ( ϵ ⋅ k ) sin ( ⋯ ) = 0 は条件 ϵ ⋅ k = 0 \boldsymbol{\epsilon}\cdot\boldsymbol{k}=0 ϵ ⋅ k = 0 から従い、□ A = 0 \Box\boldsymbol{A}=\boldsymbol{0} □ A = 0 は − ∣ k ∣ 2 + ω 2 / c 2 = 0 -|\boldsymbol{k}|^2+\omega^2/c^2=0 − ∣ k ∣ 2 + ω 2 / c 2 = 0 、すなわち ω = c ∣ k ∣ \omega=c|\boldsymbol{k}| ω = c ∣ k ∣ を要求します。場を計算します。
E = − ∂ A ∂ t = − A 0 ω ϵ sin ( k ⋅ x − ω t ) , \boldsymbol{E}=-\frac{\partial\boldsymbol{A}}{\partial t}=-A_0\omega\,\boldsymbol{\epsilon}\sin(\boldsymbol{k}\cdot\boldsymbol{x}-\omega t), E = − ∂ t ∂ A = − A 0 ω ϵ sin ( k ⋅ x − ω t ) , B = rot A = − A 0 sin ( k ⋅ x − ω t ) k × ϵ . \boldsymbol{B}=\operatorname{rot}\boldsymbol{A}=-A_0\sin(\boldsymbol{k}\cdot\boldsymbol{x}-\omega t)\,\boldsymbol{k}\times\boldsymbol{\epsilon}. B = rot A = − A 0 sin ( k ⋅ x − ω t ) k × ϵ . (2 番目は rot ( f ϵ ) = ( grad f ) × ϵ \operatorname{rot}(f\boldsymbol{\epsilon})=(\operatorname{grad}f)\times\boldsymbol{\epsilon} rot ( f ϵ ) = ( grad f ) × ϵ と grad cos ( k ⋅ x − ω t ) = − k sin ( ⋯ ) \operatorname{grad}\cos(\boldsymbol{k}\cdot\boldsymbol{x}-\omega t)=-\boldsymbol{k}\sin(\cdots) grad cos ( k ⋅ x − ω t ) = − k sin ( ⋯ ) から出ます。)大きさの比は
∣ E ∣ ∣ B ∣ = A 0 ω A 0 ∣ k ∣ = ω ∣ k ∣ = c \frac{|\boldsymbol{E}|}{|\boldsymbol{B}|}=\frac{A_0\omega}{A_0|\boldsymbol{k}|}=\frac{\omega}{|\boldsymbol{k}|}=c ∣ B ∣ ∣ E ∣ = A 0 ∣ k ∣ A 0 ω = ∣ k ∣ ω = c で、E ⊥ B ⊥ k \boldsymbol{E}\perp\boldsymbol{B}\perp\boldsymbol{k} E ⊥ B ⊥ k も見て取れます。マクスウェル方程式と電磁波 で場から直接導いた性質が、ポテンシャルからは 1 つのベクトル ϵ \boldsymbol{\epsilon} ϵ の選び方(偏光)として整理されました。ϵ \boldsymbol{\epsilon} ϵ は k \boldsymbol{k} k に直交する 2 次元平面の中で選べるので、独立な偏光は 2 種類。Remark 5.4 の数え上げと一致します。
ここまで φ \varphi φ と A \boldsymbol{A} A を別々の量として扱ってきましたが、両者は 1 つの対象の成分です。特殊相対論の記法を最小限だけ導入して、それを見ます。計量は η μ ν = d i a g ( + 1 , − 1 , − 1 , − 1 ) \eta_{\mu\nu}=\mathrm{diag}(+1,-1,-1,-1) η μν = diag ( + 1 , − 1 , − 1 , − 1 ) 、座標は x μ = ( c t , x ) x^\mu=(ct,\boldsymbol{x}) x μ = ( c t , x ) 、微分は ∂ μ = ∂ / ∂ x μ = ( c − 1 ∂ t , ∇ ) \partial_\mu=\partial/\partial x^\mu=(c^{-1}\partial_t,\nabla) ∂ μ = ∂ / ∂ x μ = ( c − 1 ∂ t , ∇ ) とします。
Definition 6.1 (4 元ポテンシャルと 4 元電流 )
A μ = ( φ c , A ) , J μ = ( c ρ , j ) A^\mu=\Bigl(\frac{\varphi}{c},\ \boldsymbol{A}\Bigr),\qquad
J^\mu=(c\rho,\ \boldsymbol{j}) A μ = ( c φ , A ) , J μ = ( c ρ , j ) をそれぞれ4 元ポテンシャル 、4 元電流密度 という。また
F μ ν = ∂ μ A ν − ∂ ν A μ F^{\mu\nu}=\partial^\mu A^\nu-\partial^\nu A^\mu F μν = ∂ μ A ν − ∂ ν A μ を電磁場テンソル という。
Theorem 6.2 (マクスウェル方程式の共変形 )
Definition 6.1 の記法のもとで、次が成り立つ。
斉次マクスウェル方程式は F μ ν F^{\mu\nu} F μν の定義から恒等的に従う。実際、∂ λ F μ ν + ∂ μ F ν λ + ∂ ν F λ μ = 0 \partial^\lambda F^{\mu\nu}+\partial^\mu F^{\nu\lambda}+\partial^\nu F^{\lambda\mu}=0 ∂ λ F μν + ∂ μ F ν λ + ∂ ν F λ μ = 0 が A μ A^\mu A μ の C 2 C^2 C 2 級性のみから成り立つ。
ゲージ変換は A μ → A μ − ∂ μ χ A^\mu\to A^\mu-\partial^\mu\chi A μ → A μ − ∂ μ χ と書け、F μ ν F^{\mu\nu} F μν を変えない。
ローレンツゲージ条件は ∂ μ A μ = 0 \partial_\mu A^\mu=0 ∂ μ A μ = 0 と書ける。
ローレンツゲージのもとで非斉次マクスウェル方程式は ∂ μ ∂ μ A ν = μ 0 J ν \partial_\mu\partial^\mu A^\nu=\mu_0 J^\nu ∂ μ ∂ μ A ν = μ 0 J ν となる。ここで ∂ μ ∂ μ = c − 2 ∂ t 2 − ∇ 2 = − □ \partial_\mu\partial^\mu=c^{-2}\partial_t^2-\nabla^2=-\Box ∂ μ ∂ μ = c − 2 ∂ t 2 − ∇ 2 = − □ である。
Proof(Theorem 6.2) 1. F μ ν = ∂ μ A ν − ∂ ν A μ F^{\mu\nu}=\partial^\mu A^\nu-\partial^\nu A^\mu F μν = ∂ μ A ν − ∂ ν A μ を巡回和に代入すると、∂ λ ∂ μ A ν \partial^\lambda\partial^\mu A^\nu ∂ λ ∂ μ A ν の形の項が 6 個現れ、偏微分の順序交換(シュワルツの定理)により符号違いで 2 個ずつ打ち消し合います。これは第 2 節で見た div rot = 0 \operatorname{div}\operatorname{rot}=0 div rot = 0 、rot grad = 0 \operatorname{rot}\operatorname{grad}=\boldsymbol{0} rot grad = 0 を 4 次元でまとめた形です。
2. ∂ μ = ( c − 1 ∂ t , − ∇ ) \partial^\mu=(c^{-1}\partial_t,-\nabla) ∂ μ = ( c − 1 ∂ t , − ∇ ) なので ∂ μ χ = ( c − 1 ∂ t χ , − grad χ ) \partial^\mu\chi=(c^{-1}\partial_t\chi,\,-\operatorname{grad}\chi) ∂ μ χ = ( c − 1 ∂ t χ , − grad χ ) です。よって
A μ − ∂ μ χ = ( φ c − 1 c ∂ χ ∂ t , A + grad χ ) = ( φ ′ c , A ′ ) A^\mu-\partial^\mu\chi=\Bigl(\frac{\varphi}{c}-\frac{1}{c}\frac{\partial\chi}{\partial t},\ \boldsymbol{A}+\operatorname{grad}\chi\Bigr)
=\Bigl(\frac{\varphi'}{c},\ \boldsymbol{A}'\Bigr) A μ − ∂ μ χ = ( c φ − c 1 ∂ t ∂ χ , A + grad χ ) = ( c φ ′ , A ′ ) となり、Definition 4.1 のゲージ変換に一致します。F μ ν F^{\mu\nu} F μν の変化分は − ( ∂ μ ∂ ν − ∂ ν ∂ μ ) χ = 0 -(\partial^\mu\partial^\nu-\partial^\nu\partial^\mu)\chi=0 − ( ∂ μ ∂ ν − ∂ ν ∂ μ ) χ = 0 です。
3. ∂ μ A μ = 1 c ∂ ∂ t ( φ c ) + div A = 1 c 2 ∂ φ ∂ t + div A \partial_\mu A^\mu=\frac{1}{c}\frac{\partial}{\partial t}\Bigl(\frac{\varphi}{c}\Bigr)+\operatorname{div}\boldsymbol{A}=\frac{1}{c^2}\frac{\partial\varphi}{\partial t}+\operatorname{div}\boldsymbol{A} ∂ μ A μ = c 1 ∂ t ∂ ( c φ ) + div A = c 2 1 ∂ t ∂ φ + div A で、Definition 5.1 と同じです。
4. ν = 0 \nu=0 ν = 0 成分を確かめます。Corollary 5.3 の □ φ = − ρ / ε 0 \Box\varphi=-\rho/\varepsilon_0 □ φ = − ρ / ε 0 を c c c で割ると □ ( φ / c ) = − ρ / ( c ε 0 ) \Box(\varphi/c)=-\rho/(c\varepsilon_0) □ ( φ / c ) = − ρ / ( c ε 0 ) 。∂ μ ∂ μ = − □ \partial_\mu\partial^\mu=-\Box ∂ μ ∂ μ = − □ なので ∂ μ ∂ μ A 0 = ρ / ( c ε 0 ) \partial_\mu\partial^\mu A^0=\rho/(c\varepsilon_0) ∂ μ ∂ μ A 0 = ρ / ( c ε 0 ) 。ここで c 2 = 1 / ( ε 0 μ 0 ) c^2=1/(\varepsilon_0\mu_0) c 2 = 1/ ( ε 0 μ 0 ) より 1 / ( c ε 0 ) = μ 0 c 1/(c\varepsilon_0)=\mu_0 c 1/ ( c ε 0 ) = μ 0 c なので、右辺は μ 0 c ρ = μ 0 J 0 \mu_0 c\rho=\mu_0J^0 μ 0 c ρ = μ 0 J 0 です。空間成分は □ A = − μ 0 j \Box\boldsymbol{A}=-\mu_0\boldsymbol{j} □ A = − μ 0 j をそのまま書き換えて ∂ μ ∂ μ A = μ 0 j \partial_\mu\partial^\mu\boldsymbol{A}=\mu_0\boldsymbol{j} ∂ μ ∂ μ A = μ 0 j です。
∎ Exercise 7.2 標準
任意の電磁ポテンシャルが φ ≡ 0 \varphi\equiv0 φ ≡ 0 となるようにゲージ変換できることを示し(Example 4.4 の復習)、その条件 φ = 0 \varphi=0 φ = 0 を保つゲージ変換をすべて決定してください。
Solution χ ( x , t ) = ∫ t 0 t φ ( x , t ′ ) d t ′ \chi(\boldsymbol{x},t)=\int_{t_0}^{t}\varphi(\boldsymbol{x},t')\,dt' χ ( x , t ) = ∫ t 0 t φ ( x , t ′ ) d t ′ とおくと、被積分関数が t ′ t' t ′ について連続なら微分積分学の基本定理により ∂ t χ ( x , t ) = φ ( x , t ) \partial_t\chi(\boldsymbol{x},t)=\varphi(\boldsymbol{x},t) ∂ t χ ( x , t ) = φ ( x , t ) です。よってゲージ変換後は φ ′ = φ − ∂ t χ = 0 \varphi'=\varphi-\partial_t\chi=0 φ ′ = φ − ∂ t χ = 0 となります。
残留自由度を求めます。φ ′ = 0 \varphi'=0 φ ′ = 0 からさらにゲージ関数 ψ \psi ψ で変換して φ ′ ′ = φ ′ − ∂ t ψ = − ∂ t ψ \varphi''=\varphi'-\partial_t\psi=-\partial_t\psi φ ′′ = φ ′ − ∂ t ψ = − ∂ t ψ 。これが 0 0 0 に等しい条件は ∂ t ψ = 0 \partial_t\psi=0 ∂ t ψ = 0 、すなわち ψ \psi ψ が時間に依らないことです。したがって φ ≡ 0 \varphi\equiv0 φ ≡ 0 を保つゲージ変換は
A → A + grad ψ ( x ) , ψ は t に依らない \boldsymbol{A}\to\boldsymbol{A}+\operatorname{grad}\psi(\boldsymbol{x}),\qquad \psi \text{ は } t \text{ に依らない} A → A + grad ψ ( x ) , ψ は t に依らない に限られます。時間ゲージを課してもなお、時間に依らないゲージ関数 1 つ分の自由度が残ることになります。
Exercise 7.3 標準
ローレンツゲージにおける波動方程式 Corollary 5.3 から、電荷保存則
∂ ρ ∂ t + div j = 0 \frac{\partial\rho}{\partial t}+\operatorname{div}\boldsymbol{j}=0 ∂ t ∂ ρ + div j = 0 が従うことを示してください。(ヒント:ゲージ条件の左辺に □ \Box □ を作用させます。)
Solution □ A = − μ 0 j \Box\boldsymbol{A}=-\mu_0\boldsymbol{j} □ A = − μ 0 j の両辺の発散を取ります。□ \Box □ は定数係数の微分作用素なので div \operatorname{div} div と交換し、
□ ( div A ) = − μ 0 div j . \Box(\operatorname{div}\boldsymbol{A})=-\mu_0\operatorname{div}\boldsymbol{j}. □ ( div A ) = − μ 0 div j . また □ φ = − ρ / ε 0 \Box\varphi=-\rho/\varepsilon_0 □ φ = − ρ / ε 0 の両辺を t t t で微分し c − 2 c^{-2} c − 2 を掛けると
□ ( 1 c 2 ∂ φ ∂ t ) = − 1 c 2 ε 0 ∂ ρ ∂ t = − μ 0 ∂ ρ ∂ t \Box\Bigl(\frac{1}{c^2}\frac{\partial\varphi}{\partial t}\Bigr)=-\frac{1}{c^2\varepsilon_0}\frac{\partial\rho}{\partial t}=-\mu_0\frac{\partial\rho}{\partial t} □ ( c 2 1 ∂ t ∂ φ ) = − c 2 ε 0 1 ∂ t ∂ ρ = − μ 0 ∂ t ∂ ρ (最後の等号で c 2 = 1 / ( ε 0 μ 0 ) c^2=1/(\varepsilon_0\mu_0) c 2 = 1/ ( ε 0 μ 0 ) 、すなわち 1 / ( c 2 ε 0 ) = μ 0 1/(c^2\varepsilon_0)=\mu_0 1/ ( c 2 ε 0 ) = μ 0 を使いました)。両者を足すと
□ ( div A + 1 c 2 ∂ φ ∂ t ) = − μ 0 ( div j + ∂ ρ ∂ t ) . \Box\Bigl(\operatorname{div}\boldsymbol{A}+\frac{1}{c^2}\frac{\partial\varphi}{\partial t}\Bigr)=-\mu_0\Bigl(\operatorname{div}\boldsymbol{j}+\frac{\partial\rho}{\partial t}\Bigr). □ ( div A + c 2 1 ∂ t ∂ φ ) = − μ 0 ( div j + ∂ t ∂ ρ ) . 左辺の括弧はローレンツゲージ条件 Definition 5.1 によりゼロなので、左辺は 0 0 0 です。μ 0 ≠ 0 \mu_0\neq0 μ 0 = 0 より div j + ∂ t ρ = 0 \operatorname{div}\boldsymbol{j}+\partial_t\rho=0 div j + ∂ t ρ = 0 が得られます。
電荷が保存しない状況ではローレンツゲージの波動方程式が矛盾する、ということでもあります。ゲージ不変性と電荷保存則が表裏の関係にあることの、初等的な現れです。
Exercise 7.4 難
Example 3.5 のソレノイド外部領域 Ω = { ( x , y , z ) ∣ x 2 + y 2 > a 2 } \Omega=\{(x,y,z)\mid x^2+y^2>a^2\} Ω = {( x , y , z ) ∣ x 2 + y 2 > a 2 } を考えます。
(1) Ω \Omega Ω 上で B = 0 \boldsymbol{B}=\boldsymbol{0} B = 0 、したがって rot A = 0 \operatorname{rot}\boldsymbol{A}=\boldsymbol{0} rot A = 0 であるにもかかわらず、A \boldsymbol{A} A をゲージ変換で恒等的に 0 \boldsymbol{0} 0 にはできないことを示してください。
(2) Ω \Omega Ω の点 1 つを固定し、そこを含む十分小さい球体 U ⊂ Ω U\subset\Omega U ⊂ Ω に制限すれば、U U U 上では A \boldsymbol{A} A をゲージ変換で 0 \boldsymbol{0} 0 にできることを示してください。(1) と矛盾しない理由を述べてください。
Solution (1) 背理法で示します。ゲージ関数 χ \chi χ があって A + grad χ = 0 \boldsymbol{A}+\operatorname{grad}\chi=\boldsymbol{0} A + grad χ = 0 、すなわち A = − grad χ \boldsymbol{A}=-\operatorname{grad}\chi A = − grad χ が Ω \Omega Ω 全体で成り立ったとします。χ \chi χ は Ω \Omega Ω 上で定義された 1 価の C 2 C^2 C 2 級関数です。半径 r > a r>a r > a の円 C C C (z = z= z = 一定、中心は z z z 軸上)に沿って線積分すると、勾配場の閉曲線に沿う積分は始点と終点が同じなので
∮ C A ⋅ d l = − ∮ C grad χ ⋅ d l = − [ χ ] 始点 終点 = 0. \oint_C\boldsymbol{A}\cdot d\boldsymbol{l}=-\oint_C\operatorname{grad}\chi\cdot d\boldsymbol{l}=-\bigl[\chi\bigr]_{\text{始点}}^{\text{終点}}=0 . ∮ C A ⋅ d l = − ∮ C grad χ ⋅ d l = − [ χ ] 始点 終点 = 0. ところが Example 3.5 の計算から
∮ C A ⋅ d l = 2 π r ⋅ Φ 2 π r = Φ = π a 2 B ≠ 0 \oint_C\boldsymbol{A}\cdot d\boldsymbol{l}=2\pi r\cdot\frac{\Phi}{2\pi r}=\Phi=\pi a^2B\neq0 ∮ C A ⋅ d l = 2 π r ⋅ 2 π r Φ = Φ = π a 2 B = 0 (B ≠ 0 B\neq0 B = 0 のとき)。矛盾です。したがってそのような χ \chi χ は存在しません。
(2) U U U は球体なので星形です。U U U 上で rot A = 0 \operatorname{rot}\boldsymbol{A}=\boldsymbol{0} rot A = 0 ですから、Lemma 2.1 により U U U 上で A = grad f \boldsymbol{A}=\operatorname{grad}f A = grad f となる f f f が存在します。χ = − f \chi=-f χ = − f とすれば U U U 上で A + grad χ = 0 \boldsymbol{A}+\operatorname{grad}\chi=\boldsymbol{0} A + grad χ = 0 です。
矛盾しない理由は、(2) で得た χ \chi χ が局所的にしか定義されていないことです。Ω \Omega Ω を球体で覆って各所で χ \chi χ を作っても、z z z 軸のまわりを 1 周する間に値が Φ \Phi Φ だけずれてしまい、Ω \Omega Ω 全体で 1 価な関数として貼り合わせられません。「回転がゼロなら勾配」という主張は局所的には常に正しく、大域的には領域の位相(ここでは縮められないループの存在)に依存する、というのがポアンカレの補題の正確な内容です。物理的には、A \boldsymbol{A} A の周回積分すなわち磁束が、外部領域だけを見ていても消せない情報として残ります。
砂川重信『理論電磁気学(第 3 版)』紀伊國屋書店、1999 — 電磁ポテンシャルとゲージ変換、遅延ポテンシャルを扱う章。
太田浩一『マクスウェル理論の基礎 ― 相対論と電磁気学』東京大学出版会、2002 — ポテンシャルの導入と 4 元形式の関係を丁寧に扱っています。
J. D. Jackson, Classical Electrodynamics , 3rd ed., Wiley, 1999 — 第 6 章(ゲージ変換、ローレンツゲージとクーロンゲージ、横成分電流)、第 11 章(特殊相対論と共変形式)。
L. D. ランダウ・E. M. リフシッツ『場の古典論(原書第 6 版)』東京図書 — 第 3 章「場の中の電荷」。4 元ポテンシャルとゲージ不変性を出発点に据えた構成です。
J. D. Jackson and L. B. Okun, “Historical roots of gauge invariance”, Reviews of Modern Physics 73 (2001), 663–680. DOI: 10.1103/RevModPhys.73.663 — ローレンツゲージの命名を含むゲージ不変性の歴史。
Y. Aharonov and D. Bohm, “Significance of Electromagnetic Potentials in the Quantum Theory”, Physical Review 115 (1959), 485–491. DOI: 10.1103/PhysRev.115.485
R. P. Feynman, R. B. Leighton, M. Sands, The Feynman Lectures on Physics , Vol. II, Chapter 15 “The Vector Potential” — オンライン版
横成分と縦成分への分解。 クーロンゲージで φ \varphi φ が瞬間的に決まってしまう(Proposition 5.6 )のに、なぜ因果性が壊れないのか。その仕組みを見ておきます。全空間で十分速く減衰する電流密度 j \boldsymbol{j} j は、ヘルムホルツ分解により
j = j L + j T , rot j L = 0 , div j T = 0 \boldsymbol{j}=\boldsymbol{j}_{\mathrm{L}}+\boldsymbol{j}_{\mathrm{T}},\qquad
\operatorname{rot}\boldsymbol{j}_{\mathrm{L}}=\boldsymbol{0},\quad \operatorname{div}\boldsymbol{j}_{\mathrm{T}}=0 j = j L + j T , rot j L = 0 , div j T = 0 と一意に分けられます。j L \boldsymbol{j}_{\mathrm{L}} j L を縦成分 、j T \boldsymbol{j}_{\mathrm{T}} j T を横成分 といいます。Lemma 2.1 により j L = grad ψ \boldsymbol{j}_{\mathrm{L}}=\operatorname{grad}\psi j L = grad ψ と書け、両辺の発散を取ると ∇ 2 ψ = div j \nabla^2\psi=\operatorname{div}\boldsymbol{j} ∇ 2 ψ = div j です。無限遠で減衰する解は
ψ ( x , t ) = − 1 4 π ∫ ( div ′ j ) ( x ′ , t ) ∣ x − x ′ ∣ d 3 x ′ \psi(\boldsymbol{x},t)=-\frac{1}{4\pi}\int\frac{(\operatorname{div}'\boldsymbol{j})(\boldsymbol{x}',t)}{|\boldsymbol{x}-\boldsymbol{x}'|}\,d^3x' ψ ( x , t ) = − 4 π 1 ∫ ∣ x − x ′ ∣ ( div ′ j ) ( x ′ , t ) d 3 x ′ です(∇ 2 ( − 1 4 π ∣ x − x ′ ∣ ) = δ 3 ( x − x ′ ) \nabla^2\bigl(-\frac{1}{4\pi|\boldsymbol{x}-\boldsymbol{x}'|}\bigr)=\delta^3(\boldsymbol{x}-\boldsymbol{x}') ∇ 2 ( − 4 π ∣ x − x ′ ∣ 1 ) = δ 3 ( x − x ′ ) を使いました)。
φ \varphi φ の時間変化がちょうど縦成分を作る。 Proposition 5.6 のクーロン積分を t t t で微分し、電荷保存則 ∂ t ρ = − div j \partial_t\rho=-\operatorname{div}\boldsymbol{j} ∂ t ρ = − div j を代入すると
ε 0 ∂ φ ∂ t = 1 4 π ∫ ∂ t ρ ( x ′ , t ) ∣ x − x ′ ∣ d 3 x ′ = − 1 4 π ∫ ( div ′ j ) ( x ′ , t ) ∣ x − x ′ ∣ d 3 x ′ = ψ ( x , t ) . \varepsilon_0\frac{\partial\varphi}{\partial t}
=\frac{1}{4\pi}\int\frac{\partial_t\rho(\boldsymbol{x}',t)}{|\boldsymbol{x}-\boldsymbol{x}'|}d^3x'
=-\frac{1}{4\pi}\int\frac{(\operatorname{div}'\boldsymbol{j})(\boldsymbol{x}',t)}{|\boldsymbol{x}-\boldsymbol{x}'|}d^3x'
=\psi(\boldsymbol{x},t). ε 0 ∂ t ∂ φ = 4 π 1 ∫ ∣ x − x ′ ∣ ∂ t ρ ( x ′ , t ) d 3 x ′ = − 4 π 1 ∫ ∣ x − x ′ ∣ ( div ′ j ) ( x ′ , t ) d 3 x ′ = ψ ( x , t ) . したがって ε 0 grad ∂ t φ = grad ψ = j L \varepsilon_0\operatorname{grad}\partial_t\varphi=\operatorname{grad}\psi=\boldsymbol{j}_{\mathrm{L}} ε 0 grad ∂ t φ = grad ψ = j L です。
結論。 §5.1 で得た一般の A \boldsymbol{A} A の方程式 □ A − grad ( div A + c − 2 ∂ t φ ) = − μ 0 j \Box\boldsymbol{A}-\operatorname{grad}\bigl(\operatorname{div}\boldsymbol{A}+c^{-2}\partial_t\varphi\bigr)=-\mu_0\boldsymbol{j} □ A − grad ( div A + c − 2 ∂ t φ ) = − μ 0 j に div A = 0 \operatorname{div}\boldsymbol{A}=0 div A = 0 を代入し、c − 2 = ε 0 μ 0 c^{-2}=\varepsilon_0\mu_0 c − 2 = ε 0 μ 0 を使うと
□ A = − μ 0 j + μ 0 ε 0 grad ∂ φ ∂ t = − μ 0 j + μ 0 j L = − μ 0 j T . \Box\boldsymbol{A}=-\mu_0\boldsymbol{j}+\mu_0\varepsilon_0\operatorname{grad}\frac{\partial\varphi}{\partial t}
=-\mu_0\boldsymbol{j}+\mu_0\boldsymbol{j}_{\mathrm{L}}
=-\mu_0\boldsymbol{j}_{\mathrm{T}} . □ A = − μ 0 j + μ 0 ε 0 grad ∂ t ∂ φ = − μ 0 j + μ 0 j L = − μ 0 j T . つまりクーロンゲージでは、A \boldsymbol{A} A は電流の横成分だけを源とする波動方程式に従います。φ \varphi φ の瞬間的な部分は A \boldsymbol{A} A の方程式から縦成分を差し引く形で正確に相殺されており、両者を組み合わせて作る E \boldsymbol{E} E と B \boldsymbol{B} B は光速でしか伝わりません。ゲージに依存する量が非因果的に見えても、ゲージ不変な量が因果的であればよい。Theorem 4.3 が保証する「物理はゲージ軌道の上でだけ定義される」という見方が、ここで具体的に働いています。