電磁ポテンシャルの導入 のバックアップ差分(No.9)
更新- バックアップ一覧
- 現在との差分 を表示
- ソース を表示
- バックアップ を表示
- 電磁気学/電磁ポテンシャルの導入 へ行く。
- 追加された行はこの色です。
- 削除された行はこの色です。
[[電磁気学]]
#contents
* 電荷密度と電流密度を与えて Maxwell 方程式を解く問題を考える [#c9b93e5b]
Maxwell 方程式
&math(\left\{\begin{array}{c@{\ }l@{\ \ \ \ \ \ \ \ \ \ \ \ }l}\displaystyle {\rm rot}\bm E+\frac{\partial \bm B}{\partial t}&=\bm 0&\MARU{1}\vspace{2mm}\\{\rm div} \bm B &= 0&\MARU{2}\vspace{2mm}\\\displaystyle\frac{1}{\mu_0}{\rm rot} \bm B-\varepsilon_0 \frac{\partial \bm E}{\partial t} &= \bm i&\MARU{3}\vspace{2mm}\\\varepsilon_0{\rm div} \bm E &= \rho&\MARU{4}\\\end{array}\right .);
に、
- 電荷密度 &math(\rho(\bm x,t));
- 電流密度 &math(\bm i(\bm x,t));
を与えて解き、
- 電場 &math(\bm E(\bm x,t));
- 磁束密度 &math(\bm B(\bm x,t));
を求める問題を考える。
電磁気学の典型的な問題の1つである
* 方程式が多すぎる? [#rcf08c0d]
ベクトルの &math(x,\, y,\,z); 成分を1つ1つ独立変数と考えれば
- 求める変数は6個 … &math(E_x,E_y,E_z,B_x,B_y,B_z);
- 与えられた式は8個 … &math(3+1+3+1=8);
条件が多すぎて解がないのでは?
→ 実は &math(\MARU{2}); と &math(\MARU{4}); は独立の条件になっていないため大丈夫
* ② は ① に含まれている [#tbdb1b7e]
&math(\MARU{1}); の div を取る
&math({\rm div}\,{\rm rot}\bm E+{\rm div}\frac{\partial}{\partial t}\bm B={\rm div} \bm \,0);
ベクトル解析の公式より、任意の &math(\bm E); に対して &math({\rm div}\,{\rm rot}\bm E=0); が成り立つから、
&math(\frac{\partial}{\partial t}\big({\rm div}\bm B\big)=0);
すなわち、&math({\rm div}\bm B); の値は時間によって変化しない。
したがって、ある時刻、例えば &math(t=0); において &math({\rm div}\bm B(\bm x,0)=0); が成り立っていれば、&math(\MARU{2}); は無くても &math(\MARU{1}); のみから任意の時刻 &math(t); において &math({\rm div}\bm B(\bm x,t)=0); が保証される。
すなわち、&math({\rm div}\bm B(\bm x,0)=0); となる「境界条件」を与える限り、&math(\MARU{2}); は &math(\MARU{1}); に含まれており、独立な条件にはならない。
* ④ は ③ に含まれている [#tbdb1b7e]
同様に &math(\MARU{3}); の div を取る
&math(\frac{1}{\mu_0}{\rm div}\,{\rm rot} \bm B-\varepsilon_0{\rm div}\frac{\partial \bm E}{\partial t} = {\rm div}\bm i);
&math(-\frac{\partial}{\partial t}\varepsilon_0{\rm div}\bm E = {\rm div}\bm i);
ここで、条件として与える &math(\rho,\,\bm i); は電荷保存則を満たしていなければならないため、
&math({\rm div}\bm i=-\frac{\partial}{\partial t}\rho);
したがって、
&math(\frac{\partial}{\partial t}\varepsilon_0{\rm div}\bm E=\frac{\partial}{\partial t}\rho);
&math(\frac{\partial}{\partial t}\big(\varepsilon_0{\rm div}\bm E-\rho\big)=0);
すなわち、&math(\varepsilon_0{\rm div}\bm E-\rho); は時間によって変化しない。
したがって、ある時刻、例えば &math(t=0); において &math(\varepsilon_0{\rm div}\bm E(\bm x,0)=\rho(\bm x,0)); が成り立っていれば、&math(\MARU{4}); は無くても &math(\MARU{3}); のみから任意の時刻 &math(t); において &math(\varepsilon_0{\rm div}\bm E(\bm x,t)=\rho(\bm x,t)); が保証される。
すなわち、&math(\varepsilon_0{\rm div}\bm E(\bm x,0)=\rho(\bm x,0)); となる「境界条件」を与える限り、&math(\MARU{4}); は &math(\MARU{3}); に含まれており、独立な条件にはならない。
* 独立な条件の数 [#y6577d4e]
&math(\MARU{1}\sim\MARU{4}); のうち、&math(\MARU{2}, \MARU{4}); を除くと、独立な条件は &math(\MARU{1}, \MARU{3}); のみとなり、成分で考えれば
- 求める変数は6個 … &math(E_x,E_y,E_z,B_x,B_y,B_z);
- 与えられた式は6個 … &math(3+3=6);
で、ちょうど解ける計算になる。
* 電磁ポテンシャルとは [#b08b2f1b]
&math(\MARU{1}, \MARU{3}); の式は
&math(\left\{\begin{array}{c@{\ }l@{\ \ \ \ \ \ \ \ \ \ \ \ }l}\displaystyle {\rm rot}\bm E+\frac{\partial \bm B}{\partial t}&=\bm 0&\MARU{1}\vspace{2mm}\\\displaystyle\frac{1}{\mu_0}{\rm rot} \bm B-\varepsilon_0 \frac{\partial \bm E}{\partial t} &= \bm i&\MARU{3}\vspace{2mm}\end{array}\right .);
のように &math(\bm E, \bm B); の各成分が複雑に入り組んでいるため、解くのが難しい。
「電磁ポテンシャル」を導入することで問題を見通しの良い形に直せる。
ここで、「電磁ポテンシャル」と「静電ポテンシャル」の2つの用語を使い分けていることに注意せよ。
- 電磁ポテンシャル … 電磁場が時間に依存して変化する時に使用する
-- ベクトルポテンシャル
-- スカラーポテンシャル
- 静電ポテンシャル … 電磁場が時間に依存しない時に使用する~
(電磁ポテンシャルのスカラーポテンシャルから時間に依存する項を除いた形)
* 静電ポテンシャルの復習 [#y4833a49]
時間に依存しない電場を考えるとき、静電ポテンシャルを使うことで問題を簡単に解くことができた。
静電ポテンシャル:&math(\phi(\bm x));
2点 &math(A, B); が与えられたとき、2点間の電位差(電場を &math(AB); 間で線積分した値)は点 &math(A, B); における &math(\phi(\bm x)); の値のみを使って、
&math(\Delta \phi_{BA}=\int_{\bm x_A}^{\bm x_B}\bm E\cdot d\bm r=\phi(\bm x_B)-\phi(\bm x_A));
と書ける。あるいは同じ意味であるが、
&math(\bm E(\bm x)=-{\rm grad} \phi(\bm x));
と書ける。
すなわち電場 &math(\bm E); は静電ポテンシャル &math(\phi); の傾きである。
というのが静電ポテンシャルの定義であった。
このような静電ポテンシャルが定義可能であることは、
上記 Maxwell 方程式の &math(\MARU{1}); から時間に依存する項を除いた
&math({\rm rot}\bm E(\bm x)=0);
という式から導かれる。
&attachref(静電ポテンシャル.png,,50%);
2点 &math(A, B); 間に2つの経路 &math(R_1, R_2); を考える。
&math(R_1); を行って &math(R_2); を帰る閉曲線を &math(C=R_1-R_2);、&math(C); に囲まれる面積を &math(S); とすると、
&math(\int_{R_1}\bm E\cdot d\bm r-\int_{R_2}\bm E\cdot d\bm r=\int_{C}\bm E\cdot d\bm r=\oint_S{\rm rot}\bm E\cdot\bm n dS=0);
となり、任意の経路に沿った線積分が同じ値を取ることが分かり、これがポテンシャルの定義可能性を示す。
** 静電ポテンシャルを用いて Maxwell 方程式を書き直す [#n140b745]
&math(\bm E(\bm x)=-{\rm grad} \phi(\bm x)); を &math(\MARU 4); へ代入すれば、
&math(-{\rm div}\,{\rm grad}\phi(\bm x)=-\bigtriangleup\phi(\bm x)=\rho(\bm x)/\varepsilon_0);
の形のラプラス方程式を得る。
以上、静電ポテンシャル &math(\rho); を導入することにより、
&math({\rm rot}\bm E(\bm x)=0);~
&math({\rm div}\bm E(\bm x)=\rho(\bm x)/\varepsilon_0);
の複雑な連立方程式から &math(\bm E(\bm x)); を決める難しい問題が、
&math(-\bigtriangleup\phi(\bm x)=\rho(\bm x)/\varepsilon_0);
を解いて
&math(\bm E(\bm x)=-{\rm grad} \phi(\bm x));
へ代入する、易しい問題に書き換えられたことに注目せよ。
* 電磁ポテンシャルの導入 [#q3201fd3]
上で見たとおり、動電場においては電場の線積分で静電ポテンシャルを定義できない。
時間に依存する電磁場を考えるときは、静電ポテンシャルのかわりに電磁ポテンシャルを使うことになる。
** ベクトルポテンシャル [#m5bb4434]
まず、&math(\MARU{2}); よりベクトルポテンシャル &math(\bm A(\bm x,t)); を導入し、これを用いて磁束密度を
&math(\bm B(\bm x,t)={\rm rot} \bm A(\bm x,t)\hspace{2cm}\MARU{5});
と書き表す。これにより &math(\MARU{2}); は自動的に満たされる。
&math({\rm div}\,{\rm rot} \bm A=0);
** ベクトルポテンシャルの意味 [#r35960b0]
静電ポテンシャルが2点間の電場の線積分から得られたのと比較しながら理解しよう。
閉曲線 &math(C); を考える。&math(C); を縁とする2つの曲面 &math(S_1, S_2); を取り、それらを合わせてできる閉曲面を &math(S=S_1-S_2);、&math(S); に囲まれる体積を &math(V); とする。
&attachref(ベクトルポテンシャル.png,,50%);
&math(S_1, S_2); を貫く磁束を比較すると、
&math(\int_{S_1}\bm B\cdot\bm ndS-\int_{S_2}\bm B\cdot\bm ndS=\int_{S}\bm B\cdot\bm ndS=\int_{V}{\rm div}\bm Bd^3x=0);
となり、閉曲線 &math(C); を貫く磁束 &math(N_C); は面の取り方によらず一意に定義される。
&math(N_C=\int_{S_\mathrm{any}}\bm B\cdot\bm ndS);
ここに &math(\bm B(\bm x,t)={\rm rot} \bm A(\bm x,t)); を代入すると、
&math(N_C=\int_{S}{\rm rot}\bm A\cdot\bm ndS=\oint_{C}\bm A\cdot d\bm r);
となり、&math(N_C); を &math(C); 上の &math(\bm A); の値のみを含む式で表せることになる。
2点 &math(A,B); 間の電位差 &math(\phi_\mathrm{AB}); が、2点 &math(A,B); の上での静電ポテンシャル &math(\phi); の値のみを使って &math(\phi_\mathrm{AB}=\phi(\bm x_B)-\phi(\bm x_A)); と表せたのと、
閉曲線 &math(C); を貫く磁束 &math(N_C); が、閉曲線 &math(C); 上でのベクトルポテンシャル &math(\bm A); の値のみを使って &math(N_C=\oint_{C}\bm A\cdot d\bm r); と表せたのと、
を対比させて理解したい。
** スカラーポテンシャル [#da84e85b]
&math(\MARU{1}); に &math(\bm B(\bm x,t)={\rm rot} \bm A(\bm x,t)); を代入すると、
&math({\rm rot}\bm E+\frac{\partial}{\partial t}{\rm rot}\bm A=\bm 0);
&math(\ROT\left(\bm E+\frac{\partial}{\partial t}\bm A\right)=\bm 0);
したがって、静電ポテンシャルの時の議論と同様に、
&math(\bm E+\frac{\partial}{\partial t}\bm A=-\GRAD\phi);
と置けて、
&math(\bm E=-\GRAD\phi-\frac{\partial}{\partial t}\bm A\hspace{2cm}\MARU{6});
である。これにより &math(\MARU{1}); は自動的に満たされる。
第2項がなければ静電ポテンシャルの式と同じであるが、動的な場では電磁誘導の項である第2項が現れる事に注意せよ。
すなわち、
&math(\Delta \phi_{BA}=\int_{\bm x_A}^{\bm x_B}\left(\bm E+\frac{\partial}{\partial t}\bm A\right)\cdot d\bm r=\phi(\bm x_B)-\phi(\bm x_A));
* 電磁ポテンシャルを用いて Maxwell 方程式を書き直す [#c1c983bc]
電磁ポテンシャルを導入したことで、&math(\MARU{1}, \MARU{2}); は自動的に満たされるため、&math(\MARU{3}, \MARU{4}); を &math(\bm A, \phi); の式に直すと次式を得る(詳細は省く)。
&math(\left\{\begin{array}{rl@{\ \ \ \ \ \ \ \ }l}\displaystyle {\rm grad}\left({\rm div} \bm A+\frac{1}{c^2}\frac{\partial\phi}{\partial t}\right)-\left(\bigtriangleup-\frac{1}{c^2}\frac{\partial^2}{\partial t^2}\right)\bm A&=\mu_0\bm i&\MARU{3}^\prime\\\displaystyle-\frac{\partial}{\partial t}\bigg( {\rm div} \bm A \ \,\phantom{+\frac{1}{c^2}\frac{\partial \phi}{\partial t}}\bigg)-\left(\bigtriangleup\ \,\phantom{-\frac{1}{c^2}\frac{\partial^2}{\partial t^2}}\right)\,\phi&\displaystyle=\frac{1}{\varepsilon_0}\rho&\MARU{4}^\prime\end{array}\right .);
ここで &math(\begin{cases}\displaystyle\bigtriangleup=\DIV\GRAD=\frac{\PD^2}{\PD x^2}+\frac{\PD^2}{\PD y^2}+\frac{\PD^2}{\PD z^2}\\\displaystyle\frac{1}{c^2}=\mu_0\varepsilon_0\end{cases}); であり、
&math(\bigtriangleup\bm A=\begin{pmatrix}\bigtriangleup A_x\\\bigtriangleup A_y\\\bigtriangleup A_z\end{pmatrix}); である。
さらに、ラプラシアンの4次元版であるダランベールの演算子(ダランベルシアン)を~
&math(\square=\left(\bigtriangleup-\frac{1}{c^2}\frac{\PD^2}{\PD t^2}\right)); として導入して、
&math(\left\{\begin{array}{rl@{\ \ \ \ \ \ \ \ }l}\displaystyle\GRAD\left(\DIV \bm A+\frac{1}{c^2}\frac{\PD\phi}{\PD t}\right)-\square\bm A&=\mu_0\bm i&\MARU{3}'\\\displaystyle\frac{\partial}{\partial t}\left(\DIV\bm A+\frac{1}{c^2}\frac{\PD \phi}{\PD t}\right)-\square\phi&\displaystyle=\frac{1}{\varepsilon_0}\rho&\MARU{4}'\end{array}\right .);
と書くこともできる。
この形は &math(\bm A); と &math(\phi); の各成分が複雑に混じり合った難しい問題になっている。
(特に第一項が)
* ゲージ変換 [#fde7eaa7]
電磁ポテンシャルの定義
&math(\bm B=\ROT \bm A\hspace{4cm}\MARU{5});~
&math(\bm E=-\GRAD\phi-\frac{\PD}{\PD t}\bm A\hspace{2cm}\MARU{6});
∴ &math(\bm A, \phi); を決めれば &math(\bm E, \bm B); は一意に定まる。
しかし、その逆は真ではない。~
すなわち、同じ &math(\bm E, \bm B); を与える &math(\bm A, \phi); が複数存在する。
簡単に思いつくのは、&math(\bm A_0,\phi_0); を定数として~
&math(\bm A'(\bm x,t)=\bm A(\bm x,t)+\bm A_0);~
&math(\phi'(\bm x,t)=\phi(\bm x,t)+\phi_0);~
→ 他にもたくさんある
任意の &math(\chi(\bm x,t)); に対して &math(\ROT\GRAD\chi=0); より、
&math(\bm A'(\bm x,t)=\bm A(\bm x,t)+\GRAD\chi(\bm x,t)); は &math(\bm A); と同じ &math(\bm B); を与える。
&math(\bm E=-\GRAD\phi'-\frac{\PD}{\PD t}\bm A');~
&math(=-\GRAD\phi'-\frac{\PD}{\PD t}\bm A-\frac{\PD}{\PD t}\GRAD\chi);~
&math(=-\GRAD(\phi'+\frac{\PD}{\PD t}\chi)-\frac{\PD}{\PD t}\bm A);~
より、
&math(\bm A'(\bm x,t)=\bm A(\bm x,t)+\GRAD\chi(\bm x,t));~
&math(\phi'(\bm x,t)=\phi(\bm x,t)-\frac{\PD}{\PD t}\chi(\bm x,t));
は、&math(\bm A,\phi); と同じ &math(\bm E,\bm E); を与える。
同じ &math(\bm E,\bm E); を与えるように、ある関数 &math(\chi); を使って
&math(\bm A,\phi); を &math(\bm A',\phi'); にする変換を「ゲージ変換」と言う。
物理的には、&math(\bm A,\phi); と &math(\bm A',\phi'); は同じ状態を表す。
したがって、Maxwell 方程式の解が1つ &math(\bm A,\phi); 求まれば、
それをゲージ変換した &math(\bm A',\phi'); も同じ方程式の解になる。
(同じ &math(\bm E,\bm E); を与えるのだから当然)
* Lorentz ゲージ [#zdf59729]
Maxwell 方程式を解いても、ただ1つだけの &math(\bm A,\phi); が求まるわけではなく、
解の集合は多くの関数を含むことが分かった。
そのような Maxwell 方程式の解の集合の中には、~
&math(\DIV \bm A+\frac{1}{c^2}\frac{\PD\phi}{\PD t}=0\hspace{1cm}\MARU{7});~
を満たす解が必ず存在することを証明できる。この条件は Lorentz 条件と呼ばれる。
そのような解を「Lorentz ゲージにおける解」と呼び、しばしば &math(\bm A_L, \phi_L); と書く。
(Maxwell 方程式&math(\MARU{3}',\MARU{4}');の解) ⊇ (Lorentz ゲージにおける解)
という関係である。
** Lorentz ゲージにおける解の存在 [#c0066b71]
&math(\bm A,\phi); が Maxwell 方程式の解であるとき、
&math(\bigtriangleup\chi-\frac{1}{c^2}\frac{\PD^2\chi}{\PD t^2}=-\left(\DIV \bm A+\frac{1}{c^2}\frac{\PD\phi}{\PD t}\right));
の解 &math(\chi_L); は必ず存在する。この &math(\chi_L); を用いて &math(\bm A,\phi\leftarrow\bm A_L,\phi_L); の変換を行えば、
&math(\DIV \bm A_L+\frac{1}{c^2}\frac{\PD\phi_L}{\PD t}=\DIV\bm A+\frac{1}{c^2}\frac{\PD\phi}{\PD t}+\bigtriangleup\chi_L-\frac{1}{c^2}\frac{\PD^2}{\PD t^2}\chi_L=0);~
となって、Lorentz ゲージにおける解を生成できる。
* Lorentz ゲージにおける Maxwell 方程式 [#sa3f7ffa]
&math(\MARU{3}',\MARU{4}'); の解の中に、&math(\MARU{7}); も満たす解があることが分かった。
実は、&math(\MARU{3}',\MARU{4}'); に &math(\MARU{7}); も条件として加えることで、
解空間は狭くなるにもかかわらず、問題が解きやすくなる場合が多い。
Lorentz ゲージにおいて、Maxwell 方程式は次の形となる。
&math(\left\{\begin{array}{rl@{\ \ \ \ \ \ \ \ }l}\displaystyle-\left(\bigtriangleup-\frac{1}{c^2}\frac{\PD^2}{\PD t^2}\right)\bm A_L&=\mu_0\bm i&\MARU{3}''\\\displaystyle-\left(\bigtriangleup-\frac{1}{c^2}\frac{\PD^2}{\PD t^2}\right)\phi_L&\displaystyle=\frac{1}{\varepsilon_0}\rho&\MARU{4}''\\\displaystyle\DIV \bm A_L+\frac{1}{c^2}\frac{\PD\phi_L}{\PD t}&=0&\MARU{7}\\\displaystyle\bm E&\displaystyle=-\GRAD\phi_L-\frac{\PD\bm A}{\PD t}&\MARU{6}\\\bm B&=\ROT \bm A_L&\MARU{5}\\\end{array}\right .);
あるいは、
&math(\left\{\begin{array}{rl@{\ \ \ \ \ \ \ \ }l}\displaystyle-\square\bm A&=\mu_0\bm i&\MARU{3}''\\\displaystyle-\square\phi&\displaystyle=\frac{1}{\varepsilon_0}\rho&\MARU{4}''\\\displaystyle\DIV \bm A+\frac{1}{c^2}\frac{\PD\phi}{\PD t}&=0&\MARU{7}\\\displaystyle\bm E&\displaystyle=-\GRAD\phi_L-\frac{\PD\bm A}{\PD t}&\MARU{6}\\\bm B&=\ROT \bm A_L&\MARU{5}\\\end{array}\right .);
この形を &math(x,y,z); 成分に分けると非常に対称性がよく、数学的にも解きやすい。
&math(\square=\frac{\PD^2}{\PD x^2}+\frac{\PD^2}{\PD y^2}+\frac{\PD^2}{\PD z^2}-\frac{1}{c^2}\frac{\PD^2}{\PD y^2});
&math(\DIV \bm A+\frac{1}{c^2}\frac{\PD\phi}{\PD t}=\frac{\PD A_x}{\PD x}+\frac{\PD A_y}{\PD y}+\frac{\PD A_z}{\PD z}+\frac{1}{c^2}\frac{\PD\phi}{\PD t});
* Lorentz ゲージにおける解の任意性 [#q9ba8749]
[[Lorentz ゲージにおける解の存在>#c0066b71]] での議論から、Lorentz ゲージの解 &math(\bm A_L,\phi_L); を、斉次方程式
&math(\bigtriangleup\chi_0-\frac{1}{c^2}\frac{\PD^2}{\PD t^2}\chi_0=0);
の解 &math(\chi_0); を用いてゲージ変換しても、やはり Lorentz 条件を満たすことが容易に分かる。
すなわち、Lorentz ゲージにおける解には &math(\chi_0); を用いたゲージ変換に相当する自由度が存在するため、解は一意には決まらない。
物理的な問題を解く上では、無数に存在する解のうち1つが見つかれば十分である。
Counter: 120300 (from 2010/06/03),
today: 26,
yesterday: 0