リエナール・ヴィーヘルト・ポテンシャル

リエナール・ヴィーヘルト・ポテンシャル（英: Liénard–Wiechert potentials）は点電荷の運動によって生じる古典的な電磁場を記述する、ローレンツ・ゲージにおけるベクトル・ポテンシャルとスカラー・ポテンシャルの総称である。名前は提案者であるアルフレド＝マリー・リエナール（英語版）とエミール・ヴィーヘルトに因む。 リエナール・ヴィーヘルト・ポテンシャルはマクスウェルの方程式から直接導かれ、点電荷の任意の運動に対する時間変化する電磁場を完全に、相対論的に正しく記述する。しかしながら、場を古典的に扱うため量子力学的な効果は記述できない。 波によって表される電磁輻射はリエナール＝ヴィーヘルト・ポテンシャルから得ることができる。 リエナール・ヴィーヘルト・ポテンシャルの表式は、一部を1898年にアルフレド＝マリー・リエナールが、1900年から1900年初頭にかけてエミール・ヴィーヘルトがそれぞれ独立に与えた。 古典電磁気学を規範として、アルベルト・アインシュタインは特殊相対性理論の発見に成功した。 電磁波の運動と伝播を調べることで、相対論的な空間と時間の記述を導いた。 リエナール・ヴィーヘルトの公式は運動する相対論的な粒子系のより詳細な解析をするための出発点として重要な役割を果たす。 巨視的で互いに独立に運動する粒子に対しては、リエナール・ヴィーヘルトの公式による記述は正確だが、粒子の運動が量子論的になる領域においては正確ではなくなる。 量子力学では粒子の電磁放射について制限が加わる。粒子の放射現象に関する古典的な記述は、実験結果と明らかに食い違ってしまう。例えば原子を構成する電子は古典論が示すような放射現象を起こさず、原子は安定に存在できる。このことは電子のエネルギー状態が量子化されることによって説明される。 放射を理解するには電磁場を量子化する必要があり、量子電磁力学として20世紀の後半に構築された。 リエナール・ヴィーヘルト・ポテンシャルは運動する点電荷が作る電磁ポテンシャルである。電荷 q {\displaystyle q} の位置を r s {\displaystyle {\boldsymbol {r}}_{\mathrm {s} }} 速度を v s {\displaystyle {\boldsymbol {v}}_{\mathrm {s} }} とするとスカラー・ポテンシャル φ {\displaystyle \varphi } とベクトル・ポテンシャル A {\displaystyle {\boldsymbol {A}}} は次のように表される。 φ ( r , t ) = 1 4 π ϵ 0 ( q ( 1 − n ⋅ β s ) | r − r s | ) t r , {\displaystyle \varphi ({\boldsymbol {r}},t)={\frac {1}{4\pi \epsilon _{0}}}\left({\frac {q}{(1-{\boldsymbol {n}}\cdot {\boldsymbol {\beta }}_{\mathrm {s} })|{\boldsymbol {r}}-{\boldsymbol {r}}_{\mathrm {s} }|}}\right)_{t_{\mathrm {r} }}\,,} A ( r , t ) = μ 0 c 4 π ( q β s ( 1 − n ⋅ β s ) | r − r s | ) t r = β s ( t r ) c φ ( r , t ) . {\displaystyle {\boldsymbol {A}}({\boldsymbol {r}},t)={\frac {\mu _{0}c}{4\pi }}\left({\frac {q{\boldsymbol {\beta }}_{\mathrm {s} }}{(1-{\boldsymbol {n}}\cdot {\boldsymbol {\beta }}_{\mathrm {s} })|{\boldsymbol {r}}-{\boldsymbol {r}}_{\mathrm {s} }|}}\right)_{t_{\mathrm {r} }}={\frac {{\boldsymbol {\beta }}_{\mathrm {s} }(t_{\mathrm {r} })}{c}}\varphi ({\boldsymbol {r}},t)\,.} ここで β s ( t ) {\displaystyle {\boldsymbol {\beta }}_{\mathrm {s} }(t)} および n {\displaystyle {\boldsymbol {n}}} は次のように定義される。 β s ( t ) = v s ( t ) c , {\displaystyle {\boldsymbol {\beta }}_{\mathrm {s} }(t)={\frac {{\boldsymbol {v}}_{\mathrm {s} }(t)}{c}}\,,} n = r − r s | r − r s | . {\displaystyle {\boldsymbol {n}}={\frac {{\boldsymbol {r}}-{\boldsymbol {r}}_{\mathrm {s} }}{|{\boldsymbol {r}}-{\boldsymbol {r}}_{\mathrm {s} }|}}\,.} 速度因子 β s {\displaystyle {\boldsymbol {\beta }}_{\mathrm {s} }} は光速を単位とする点電荷の速度を表し、 ベクトル n {\displaystyle {\boldsymbol {n}}} は点電荷の位置 r s {\displaystyle {\boldsymbol {r}}_{\mathrm {s} }} を原点とし位置 r {\displaystyle {\boldsymbol {r}}} へ向かう単位ベクトルを表している。 電磁ポテンシャルの定義から直接、電場および磁場を求めることができる。 E = − ∇ φ − ∂ A ∂ t , {\displaystyle {\boldsymbol {E}}=-\nabla \varphi -{\dfrac {\partial {\boldsymbol {A}}}{\partial t}}\,,} B = ∇ × A . {\displaystyle {\boldsymbol {B}}=\nabla \times {\boldsymbol {A}}\,.} 具体的な電磁ポテンシャルとしてリエナール・ヴィーヘルト・ポテンシャルを与えると次の共変でない電磁場の表式が得られる。 E ( r , t ) = 1 4 π ϵ 0 ( q ( n − β ) γ 2 ( 1 − n ⋅ β ) 3 | r − r s | 2 + q n × ( ( n − β ) × β ˙ ) c ( 1 − n ⋅ β ) 3 | r − r s | ) t r , {\displaystyle {\boldsymbol {E}}({\boldsymbol {r}},t)={\frac {1}{4\pi \epsilon _{0}}}\left({\frac {q({\boldsymbol {n}}-{\boldsymbol {\beta }})}{\gamma ^{2}(1-{\boldsymbol {n}}\cdot {\boldsymbol {\beta }})^{3}|{\boldsymbol {r}}-{\boldsymbol {r}}_{\mathrm {s} }|^{2}}}+{\frac {q{\boldsymbol {n}}\times {\big (}({\boldsymbol {n}}-{\boldsymbol {\beta }})\times {\dot {\boldsymbol {\beta }}}{\big )}}{c(1-{\boldsymbol {n}}\cdot {\boldsymbol {\beta }})^{3}|{\boldsymbol {r}}-{\boldsymbol {r}}_{\mathrm {s} }|}}\right)_{t_{\mathrm {r} }}\,,} B ( r , t ) = μ 0 4 π ( q c ( β × n ) γ 2 ( 1 − n ⋅ β ) 3 | r − r s | 2 + q n × ( n × ( ( n − β ) × β ˙ ) ) ( 1 − n ⋅ β ) 3 | r − r s | ) t r = n ( t r ) c × E ( r , t ) . {\displaystyle {\boldsymbol {B}}({\boldsymbol {r}},t)={\frac {\mu _{0}}{4\pi }}\left({\frac {qc({\boldsymbol {\beta }}\times {\boldsymbol {n}})}{\gamma ^{2}(1-{\boldsymbol {n}}\cdot {\boldsymbol {\beta }})^{3}|{\boldsymbol {r}}-{\boldsymbol {r}}_{\mathrm {s} }|^{2}}}+{\frac {q{\boldsymbol {n}}\times {\Big (}{\boldsymbol {n}}\times {\big (}({\boldsymbol {n}}-{\boldsymbol {\beta }})\times {\dot {\boldsymbol {\beta }}}{\big )}{\Big )}}{(1-{\boldsymbol {n}}\cdot {\boldsymbol {\beta }})^{3}|{\boldsymbol {r}}-{\boldsymbol {r}}_{\mathrm {s} }|}}\right)_{t_{\mathrm {r} }}={\frac {{\boldsymbol {n}}(t_{\mathrm {r} })}{c}}\times {\boldsymbol {E}}({\boldsymbol {r}},t)\,.} ここで γ ( t ) = 1 1 − | β ( t ) | 2 {\displaystyle \gamma (t)={\frac {1}{\sqrt {1-|{\boldsymbol {\beta }}(t)|^{2}}}}} はローレンツ因子である。 もし電荷が一定の速度 c β {\displaystyle c{\boldsymbol {\beta }}} で運動する場合、速度因子の時間微分 β ˙ ( t ) = d β ( t ) d t {\displaystyle {\dot {\boldsymbol {\beta }}}(t)={\frac {d{\boldsymbol {\beta }}(t)}{dt}}} はゼロになるので、電場は n − β {\displaystyle {\boldsymbol {n}}-{\boldsymbol {\beta }}} の項だけが残る。このとき電場の向きは n − β {\displaystyle {\boldsymbol {n}}-{\boldsymbol {\beta }}} によって決まる。 電場の初項は、特に速度因子がゼロ β ( t ) = 0 {\displaystyle {\boldsymbol {\beta }}(t)=0} の場合には、電荷 q {\displaystyle q} の点電荷がつくる静電場に一致し、電荷のつくる電磁場の静的な部分を表している。 第二項は運動する電荷が放射する電磁波に対応する。位置 r {\displaystyle {\boldsymbol {r}}} で電場 E ( r , t ) {\displaystyle {\boldsymbol {E}}({\boldsymbol {r}},t)} を観測するとき、 電荷と観測点を結ぶ方向に直交するように電荷が加速されると第二項の輻射項が観測される。 この輻射項の電磁場の向きは遅延時間における電荷の位置を向いている。 外界に電磁場のソースがない境界条件の下で、非斉次な波動方程式における電磁ポテンシャルの遅延解は、ローレンツ・ゲージを採用すれば、電荷密度 ρ ( r , t ) {\displaystyle \rho ({\boldsymbol {r}},t)} および電流密度 J ( r , t ) {\displaystyle {\boldsymbol {J}}({\boldsymbol {r}},t)} をソースとして、 φ ( r , t ) = 1 4 π ϵ 0 ∫ ρ ( r ′ , t r ′ ) | r − r ′ | d 3 r ′ {\displaystyle \varphi ({\boldsymbol {r}},t)={\frac {1}{4\pi \epsilon _{0}}}\int {\frac {\rho ({\boldsymbol {r}}',t_{\mathrm {r} }')}{|{\boldsymbol {r}}-{\boldsymbol {r}}'|}}d^{3}{\boldsymbol {r}}'} および A ( r , t ) = μ 0 4 π ∫ J ( r ′ , t r ′ ) | r − r ′ | d 3 r ′ {\displaystyle {\boldsymbol {A}}({\boldsymbol {r}},t)={\frac {\mu _{0}}{4\pi }}\int {\frac {{\boldsymbol {J}}({\boldsymbol {r}}',t_{\mathrm {r} }')}{|{\boldsymbol {r}}-{\boldsymbol {r}}'|}}d^{3}{\boldsymbol {r}}'} となる。ここで t r ′ = t − 1 c | r − r ′ | {\displaystyle t_{\mathrm {r} }'=t-{\frac {1}{c}}|{\boldsymbol {r}}-{\boldsymbol {r}}'|} は遅延時間 (retarded time) である。 ソースとなる点電荷の運動の軌跡は時間の関数 r s ( t ′ ) {\displaystyle {\boldsymbol {r}}_{\mathrm {s} }(t')} として与えられる。電荷密度および電流密度は以下のものを考える。 ρ ( r ′ , t ′ ) = q δ 3 ( r ′ − r s ( t ′ ) ) , {\displaystyle \rho ({\boldsymbol {r}}',t')=q\delta ^{3}({\boldsymbol {r'}}-{\boldsymbol {r}}_{\mathrm {s} }(t'))\,,} J ( r ′ , t ′ ) = q v s ( t ′ ) δ 3 ( r ′ − r s ( t ′ ) ) . {\displaystyle {\boldsymbol {J}}({\boldsymbol {r}}',t')=q{\boldsymbol {v}}_{\mathrm {s} }(t')\delta ^{3}({\boldsymbol {r'}}-{\boldsymbol {r}}_{\mathrm {s} }(t'))\,.} ここで δ 3 {\displaystyle \delta ^{3}} は 3 次元のディラックのデルタ関数であり、 v s ( t ′ ) {\displaystyle {\boldsymbol {v}}_{\mathrm {s} }(t')} は点電荷の速度である。 電荷密度と電流密度を具体的な形に書き直せば、以下の電磁ポテンシャルが与えられる。 φ ( r , t ) = 1 4 π ϵ 0 ∫ q δ 3 ( r ′ − r s ( t r ′ ) ) | r − r ′ | d 3 r ′ {\displaystyle \varphi ({\boldsymbol {r}},t)={\frac {1}{4\pi \epsilon _{0}}}\int {\frac {q\delta ^{3}({\boldsymbol {r'}}-{\boldsymbol {r}}_{\mathrm {s} }(t_{\mathrm {r} }'))}{|{\boldsymbol {r}}-{\boldsymbol {r}}'|}}d^{3}{\boldsymbol {r}}'} A ( r , t ) = μ 0 4 π ∫ q v s ( t r ′ ) δ 3 ( r ′ − r s ( t r ′ ) ) | r − r ′ | d 3 r ′ {\displaystyle {\boldsymbol {A}}({\boldsymbol {r}},t)={\frac {\mu _{0}}{4\pi }}\int {\frac {q{\boldsymbol {v}}_{\mathrm {s} }(t_{\mathrm {r} }')\delta ^{3}({\boldsymbol {r'}}-{\boldsymbol {r}}_{\mathrm {s} }(t_{r}'))}{|{\boldsymbol {r}}-{\boldsymbol {r}}'|}}d^{3}{\boldsymbol {r}}'} これらの積分を簡単にするために、デルタ関数 δ ( t ′ − t r ′ ) {\displaystyle \delta (t'-t_{\mathrm {r} }')} を用いて t r ′ {\displaystyle t_{\mathrm {r} }'} を t ′ {\displaystyle t'} で置き換え、 t ′ {\displaystyle t'} の積分に直せば次のようになる。 φ ( r , t ) = 1 4 π ϵ 0 ∬ q δ 3 ( r ′ − r s ( t ′ ) ) | r − r ′ | δ ( t ′ − t r ′ ) d t ′ d 3 r ′ , {\displaystyle \varphi ({\boldsymbol {r}},t)={\frac {1}{4\pi \epsilon _{0}}}\iint {\frac {q\delta ^{3}({\boldsymbol {r'}}-{\boldsymbol {r}}_{\mathrm {s} }(t'))}{|{\boldsymbol {r}}-{\boldsymbol {r}}'|}}\delta (t'-t_{\mathrm {r} }')\,dt'\,d^{3}{\boldsymbol {r}}'\,,} A ( r , t ) = μ 0 4 π ∬ q v s ( t ′ ) δ 3 ( r ′ − r s ( t ′ ) ) | r − r ′ | δ ( t ′ − t r ′ ) d t ′ d 3 r ′ . {\displaystyle {\boldsymbol {A}}({\boldsymbol {r}},t)={\frac {\mu _{0}}{4\pi }}\iint {\frac {q{\boldsymbol {v}}_{\mathrm {s} }(t')\delta ^{3}({\boldsymbol {r'}}-{\boldsymbol {r}}_{\mathrm {s} }(t'))}{|{\boldsymbol {r}}-{\boldsymbol {r}}'|}}\delta (t'-t_{\mathrm {r} }')\,dt'\,d^{3}{\boldsymbol {r}}'\,.} また、積分の順序を入れ替えれば以下の式を得る。 φ ( r , t ) = 1 4 π ϵ 0 ∬ δ ( t ′ − t r ′ ) | r − r ′ | q δ 3 ( r ′ − r s ( t ′ ) ) d 3 r ′ d t ′ , {\displaystyle \varphi ({\boldsymbol {r}},t)={\frac {1}{4\pi \epsilon _{0}}}\iint {\frac {\delta (t'-t_{\mathrm {r} }')}{|{\boldsymbol {r}}-{\boldsymbol {r}}'|}}q\delta ^{3}({\boldsymbol {r'}}-{\boldsymbol {r}}_{\mathrm {s} }(t'))\,d^{3}{\boldsymbol {r}}'dt'\,,} A ( r , t ) = μ 0 4 π ∬ δ ( t ′ − t r ′ ) | r − r ′ | q v s ( t ′ ) δ 3 ( r ′ − r s ( t ′ ) ) d 3 r ′ d t ′ . {\displaystyle {\boldsymbol {A}}({\boldsymbol {r}},t)={\frac {\mu _{0}}{4\pi }}\iint {\frac {\delta (t'-t_{\mathrm {r} }')}{|{\boldsymbol {r}}-{\boldsymbol {r}}'|}}q{\boldsymbol {v}}_{\mathrm {s} }(t')\delta ^{3}({\boldsymbol {r'}}-{\boldsymbol {r}}_{\mathrm {s} }(t'))\,d^{3}{\boldsymbol {r}}'dt'\,.} デルタ関数によって r ′ = r s ( t ′ ) {\displaystyle {\boldsymbol {r}}'={\boldsymbol {r}}_{\mathrm {s} }(t')} となるような対角成分だけが取り出され、内側の積分が簡単になる。 t r ′ {\displaystyle t_{\mathrm {r} }'} は r ′ {\displaystyle {\boldsymbol {r}}'} の関数であることに注意し、 t r = t r ( r s ( t ′ ) , t ′ ) {\displaystyle t_{\mathrm {r} }=t_{\mathrm {r} }({\boldsymbol {r}}_{s}(t'),t')} と書き直して積分を行う。 φ ( r , t ) = 1 4 π ϵ 0 ∫ q δ ( t ′ − t r ′ ) | r − r s ( t ′ ) | d t ′ , {\displaystyle \varphi ({\boldsymbol {r}},t)={\frac {1}{4\pi \epsilon _{0}}}\int q{\frac {\delta (t'-t_{\mathrm {r} }')}{|{\boldsymbol {r}}-{\boldsymbol {r}}_{\mathrm {s} }(t')|}}dt'\,,} A ( r , t ) = μ 0 4 π ∫ q v s ( t ′ ) δ ( t ′ − t r ′ ) | r − r s ( t ′ ) | d t ′ . {\displaystyle {\boldsymbol {A}}({\boldsymbol {r}},t)={\frac {\mu _{0}}{4\pi }}\int q{\boldsymbol {v}}_{\mathrm {s} }(t'){\frac {\delta (t'-t_{\mathrm {r} }')}{|{\boldsymbol {r}}-{\boldsymbol {r}}_{\mathrm {s} }(t')|}}\,dt'\,.} 遅延時間 t r ′ {\displaystyle t_{\mathrm {r} }'} およびソースの位置 r s ( t ′ ) {\displaystyle {\boldsymbol {r}}_{\mathrm {s} }(t')} は ( r , t ) {\displaystyle ({\boldsymbol {r}},t)} の関数であり、従って t ′ {\displaystyle t'} に依存する。 この積分を計算するために、以下の恒等式を利用する。 δ ( f ( t ′ ) ) = ∑ i δ ( t ′ − t i ) | f ′ ( t i ) | {\displaystyle \delta (f(t'))=\sum _{i}{\frac {\delta (t'-t_{i})}{|f'(t_{i})|}}} ここで t i {\displaystyle t_{i}} はそれぞれ関数 f ( t ′ ) {\displaystyle f(t')} がゼロとなるような点である。 任意の与えられた時間と空間の座標 ( r , t ) {\displaystyle ({\boldsymbol {r}},t)} とソースの軌跡 r s ( t ′ ) {\displaystyle {\boldsymbol {r}}_{\mathrm {s} }(t')} に対して、 遅延時間 t r {\displaystyle t_{\mathrm {r} }} は唯一つに定まるので、デルタ関数は次のように簡単にできる。 δ ( t ′ − t r ′ ) = δ ( t ′ − t r ) ∂ ∂ t ′ ( t ′ − t r ′ ) | t ′ = t r = δ ( t ′ − t r ) ∂ ∂ t ′ ( t ′ − ( t − 1 c | r − r s ( t ′ ) | ) ) | t ′ = t r = δ ( t ′ − t r ) 1 + 1 c ( r − r s ( t ′ ) ) / | r − r s ( t ′ ) | ⋅ ( − v s ( t ′ ) ) | t ′ = t r = δ ( t ′ − t r ) 1 − β s ⋅ ( r − r s ) / | r − r s | {\displaystyle {\begin{aligned}\delta (t'-t_{\mathrm {r} }')&={\frac {\delta (t'-t_{\mathrm {r} })}{{\frac {\partial }{\partial t'}}(t'-t_{\mathrm {r} }')|_{t'=t_{\mathrm {r} }}}}\\&={\frac {\delta (t'-t_{\mathrm {r} })}{{\frac {\partial }{\partial t'}}(t'-(t-{\frac {1}{c}}|{\boldsymbol {r}}-{\boldsymbol {r}}_{\mathrm {s} }(t')|))|_{t'=t_{\mathrm {r} }}}}\\&={\frac {\delta (t'-t_{\mathrm {r} })}{1+{\frac {1}{c}}({\boldsymbol {r}}-{\boldsymbol {r}}_{\mathrm {s} }(t'))/|{\boldsymbol {r}}-{\boldsymbol {r}}_{\mathrm {s} }(t')|\cdot (-{\boldsymbol {v}}_{\mathrm {s} }(t'))|_{t'=t_{r}}}}\\&={\frac {\delta (t'-t_{\mathrm {r} })}{1-{\boldsymbol {\beta }}_{\mathrm {s} }\cdot ({\boldsymbol {r}}-{\boldsymbol {r}}_{\mathrm {s} })/|{\boldsymbol {r}}-{\boldsymbol {r}}_{\mathrm {s} }|}}\end{aligned}}} ここで β s = v s / c {\displaystyle {\boldsymbol {\beta }}_{\mathrm {s} }={\boldsymbol {v}}_{\mathrm {s} }/c} , β s {\displaystyle {\boldsymbol {\beta }}_{\mathrm {s} }} および r s {\displaystyle {\boldsymbol {r}}_{\mathrm {s} }} は遅延時間 t r {\displaystyle t_{\mathrm {r} }} の関数であり、また恒等式 | x | ′ = x