静止物体中の Maxwell の方程式 の履歴(No.9)

更新

電磁気学

静止物体†

物体は非常に多数の荷電粒子（電子および原子核）から構成されるが、 電子と原子核は空間的に重なって存在するため、 通常ほぼすべてが互いに打ち消し合って、遠くから見た際に電荷の存在を無視できる。

そこで、

  • 打ち消される電荷を無視し、
  • 残った電荷や電流のみを考えたい

例）１ｋｇの鉄球から、１兆個 = 10^{12} 個に１個の電子を取り出して、もう１つの鉄球に移す。 それら２つの鉄球を１ｍ 離して置くと４０トン以上の力で互いに引き合う。

普段、正の電荷と負の電荷がどれほど精密に打ち消し合っているかが分かる。

打ち消されずに残るもの†

原子イオン自由電子
構成核+電子核+(内殻)電子伝導電子
特徴動けない
中性		動けない
電荷		動ける
電荷
電荷 \rho 真電荷 \rho_e
分極電荷 \rho_d
電流 \bm i 分極電流 \bm i_d 伝導電流 \bm i_e
磁化電流 \bm i_m

このように分類すると、

  • 全電荷 \rho=\rho_e+\rho_d+ (打ち消されてゼロになるため無視する電荷)
  • 全電流 \bm i=\bm i_e+\bm i_d+\bm i_m+ (打ち消されてゼロになるため無視する電流)

と書ける。

本章では、真電荷、伝導電流を主に考える（イオンと自由電子に由来）

打ち消されてゼロになる成分（中性分子、内殻電子など）は完全に無視する

薄文字 はほぼ打ち消し合った電荷の消え残った影響
→ \varepsilon,\mu を介して間接的に取り扱う。

通常、静止物体中の電磁気学ではこの \rho_e,\bm i_e を単に \rho,\bm i と書く。 真空中の電磁気学で出てくる \rho,\bm i と定義が異なることに注意せよ。

分極†

通常、原子の 負電荷（電子雲） -e の中心は、 正電荷（原子核） +e の中心と一致している。

ここに外場 \bm E がかかり、両者の中心が \bm s だけずれたとしよう。

負電荷の中心を \bm x' 、正電荷の中心を \bm x'+\bm s とする。 すなわち \bm s は負電荷から正電荷へのベクトル。


　　　　　分極.png


このとき遠く離れた \bm x における電位は、

\phi(\bm x)=\frac{1}{4\pi\varepsilon_0}\left(\frac{+e}{|\bm x-(\bm x'+\bm s)|}+\frac{-e}{|\bm x-\bm x'|}\right)

= \frac{1}{4\pi\varepsilon_0}e\bm s\cdot\GRAD\!_{\bm x'}\frac{1}{|\bm x-\bm x'|}

電気双極子能率 \bm p\equiv e\bm s を導入すると、

=\frac{1}{4\pi\varepsilon_0}\bm p\cdot\GRAD\!_{\bm x'}\frac{1}{|\bm x-\bm x'|}

=\frac{1}{4\pi\varepsilon_0}\bm p\cdot\frac{1}{|\bm x-\bm x'|^2}\frac{\bm x-\bm x'}{|\bm x-\bm x'|}

を得る。

ここで、 \frac{1}{|\bm x-\bm x'|} は \bm x' をパラメータに含む \bm x の関数であり、その値は \bm x' に対して中心対称。 \bm x' からの距離に対して 1/r で減衰する関数である。

\GRAD\!_{\bm x'}\frac{1}{|\bm x-\bm x'|}=\begin{pmatrix}\displaystyle\frac{\PD}{\PD x'_x}\frac{1}{|\bm x-\bm x'|}\vspace{2mm}\\ \displaystyle\frac{\PD}{\PD x'_y}\frac{1}{|\bm x-\bm x'|}\vspace{2mm}\\ \displaystyle\frac{\PD}{\PD x'_z}\frac{1}{|\bm x-\bm x'}\\ \end{pmatrix}

はそれを \bm x' に対して微分したベクトル値関数。 grad は勾配を表わしていて、傾きの最も大きい方向に上る向きのベクトルになる。

ここでは \bm x' に対する勾配なので、 \bm x に対する勾配とは逆になり、 下図のように \bm x' を中心にして外へ向くベクトル場となる。

grad 1 over r - 2d.png grad 1 over r - 3d.png

実際の電位分布はこれと方向ベクトルの内積をとった物となるため、 例えば \bm s が +x 方向であれば次の形になる。 右は電位分布を微分して得られる「電場」。

dielectric polarization-1.jpg dielectric polarization-2.png

電場の形は 棒磁石の周りにできる磁力線 をイメージするとよい。 事実、「非常に小さな棒磁石」は磁気双極子そのものとなる。

具体的な形を４章でもう少し詳しく見ることになる。

Mathematica ソース: fileMathematica.pdf

分極の連続分布†

電気双極子能率を持つ分子や原子が空間的に敷き詰められており、 局所的な分極率を体積に比例する形で、

\bm p_{\rm local}=\bm p(\bm x)d^3x

と書けるとする。この分布により生じる電位は、

\rho_d(\bm x)=\frac{1}{4\pi\varepsilon_0}\int_V \bm p(\bm x')\cdot\GRAD\!_{\bm x'}\frac{1}{|\bm x-\bm x'|}\,d^3x'

次式を用いて変形すると、

\DIV\!_{\bm x'}\frac{\bm p(\bm x')}{|\bm x-\bm x'|}=\frac{\DIV\!_{\bm x'}\bm p(\bm x')}{|\bm x-\bm x'|}+\bm p(\bm x')\cdot\GRAD\!_{\bm x'}\frac{1}{|\bm x-\bm x'|}

\rho_d(\bm x)=\frac{1}{4\pi\varepsilon_0}\int_V \DIV\!_{\bm x'}\frac{\bm p(\bm x')}{|\bm x-\bm x'|}\,d^3x'-\frac{1}{4\pi\varepsilon_0}\int_V \frac{\DIV\!_{\bm x'}\bm p(\bm x')}{|\bm x-\bm x'|}\,d^3x'

=\frac{1}{4\pi\varepsilon_0}\int_S \frac{\bm p(\bm x')}{|\bm x-\bm x'|}\cdot\bm n\,dS-\frac{1}{4\pi\varepsilon_0}\int_V \frac{\DIV\!_{\bm x'}\bm p(\bm x')}{|\bm x-\bm x'|}\,d^3x'

物質が系の外にはみ出していない限り、系の境界面 S において、 \bm p(\bm x')=\bm o となるから、この第１項は消える。

=\frac{1}{4\pi\varepsilon_0}\int_V \frac{-\DIV\!_{\bm x'}\bm p(\bm x')}{|\bm x-\bm x'|}\,d^3x'

この式を、真電荷 \rho_e に対する電位分布の式

\phi_e(\bm x)=\frac{1}{4\pi\varepsilon_0}\int_V \frac{\rho_e(\bm x')}{|\bm x-\bm x'|}\,d^3x'

と比較すると、

\rho_d(\bm x)=-\DIV\!_{\bm x'}\bm p(\bm x')

という関係があることが分かる。

これは、 \bm p は負電荷から正電荷に向けてのベクトルであるから、 \bm p に沸き出しがあれば、中心部には負電荷が残されるという意味として理解できる。

個々の電荷はほとんど動かない（ \bm s だけしか）のであるが、 その空間的偏りにより電荷分布が生じることになる。

電荷保存則†

真電荷としてカウントされる伝導電子やホールが、原子からイオンを残して自由に動けるようになった電荷であるのに対して、分極電荷を生じる電荷は原子核と電子が強固に結びついている。

分極電荷を生じる荷電粒子がイオン化して真電荷となったり、その逆の過程が起きない限り、 真電荷の総量と、分極電荷の総量とは独立に保存する。すなわち、

  • \DIV\bm i_e+\frac{\PD\rho_e}{\PD t}=0
  • \DIV\bm i_d+\frac{\PD\rho_d}{\PD t}=0

特に第２式は、

\DIV\bm i_d-\frac{\PD\DIV\bm p(\bm x)}{\PD t}=0

\DIV\left(\bm i_d-\frac{\PD\bm p(\bm x)}{\PD t}\right)=0

と書け、これはある関数 \bm M(\bm x) を用いて、

\bm i_d-\frac{\PD\bm p}{\PD t}=-\ROT\bm M

の関係があることを示している。

以降では、これまでの \bm i_d を \bm i_{d+m} と書き、 改めて \bm i_d を分極の時間変化によるものとして、

\bm i_d=\frac{\PD\bm p}{\PD t}

と定義する。 \bm i_m=-\ROT\bm M と置けば、上式は

\bm i_{d+m}=\bm i_d+\bm i_m

と書ける。

磁化電流†

電流が流れても電荷保存則に影響がないというのはどういう事であろうか？

電流に発散がないとき、そのような電流は必ずループの形で流れるため、電荷分布を変化させない。 そのため電荷保存則には現れないのである。

また我々は「動ける電荷」である真電荷が作る電流を伝導電流 \bm i_e としてすでに取り入れた。残りは「動けない電荷」だけであるが、動けない電荷が電流を作ることが可能であろうか？

ここで考える磁化電流 \bm i_m とは、空間的に移動しない電荷が、局所領域で回転運動することにより生じる電流である。そのような電流は微少な磁気双極子を形成する。

例えば、そのような電荷が空間のある領域に一定密度で分布するとき、内部を流れる電流の影響は互いに打ち消し合うが、領域の縁を流れる電流は打ち消されず残ることになる。結果的に、個々の電荷は大きく移動しないにもかかわらず、あたかも領域の縁を一周する電流が流れているかのような磁場を発生させることになる。

磁化電流 \bm i_m は、上で見たようにある関数 \bm M の回転で書くことができる。

\bm i_m=\ROT\bm M

この \bm M を磁化ベクトルと呼び、 \bm M を用いると局所体積あたりの磁気双極子能を \bm M d^3x と書き表せる。

物質中の Maxwell 方程式†

真空中の Maxwell 方程式の電荷および電流を

  • 全電荷 \rho=\rho_e+\rho_d+ (打ち消されてゼロになるため無視する電荷)
  • 全電流 \bm i=\bm i_e+\bm i_d+\bm i_m+ (打ち消されてゼロになるため無視する電流)

と書き直し括弧内を無視することにより、

\ROT\bm E(\bm x,t)+\frac{\PD\bm B(\bm x,t)}{\PD t}=\bm o

\DIV\bm B(\bm x,t)=0

\frac{1}{\mu_0}\ROT\bm B(\bm x,t)-\varepsilon_0\frac{\PD\bm E(\bm x,t)}{\PD t}=\bm i_e+\bm i_d+\bm i_m

\varepsilon_0\DIV\bm E(\bm x,t)=\rho_e+\rho_d

を得る。さらに、 \rho_d,\bm i_d,\bm i_m に具体的な表式を代入すれば第３、４式は、

\ROT\left\{\frac{1}{\mu_0}\bm B(\bm x,t)-\ROT\bm M(\bm x,t)\right\}-\frac{\PD}{\PD t}\left\{\varepsilon_0\bm E(\bm x,t)+\bm P(\bm x,t)\right\}=\bm i_e

\DIV\{\varepsilon_0\bm E(\bm x,t)+\bm P(\bm x,t)\}=\rho_e

と変形できる。

そこで、

電即密度： \bm D(\bm x,t)=\varepsilon_0\bm E(\bm x,t)+\bm P(\bm x,t)

磁場の強さ： \bm H(\bm x,t)=\frac{1}{\mu_0}\bm B(\bm x,t)-\ROT\bm M(\bm x,t)

と書く事にすると、物質中での Maxwell 方程式として次式を得る。

\ROT\bm E(\bm x,t)+\frac{\PD\bm B(\bm x,t)}{\PD t}=\bm o

\DIV\bm B(\bm x,t)=0

\ROT\bm H(\bm x,t)-\frac{\PD\bm D(\bm x,t)}{\PD t}=\bm i_e

\DIV\bm D(\bm x,t)=\rho_e

この形は、真空の Maxwell 方程式と非常に似た形になっているものの、 実際には非常に大きな違いを持っている。

  • ここで右辺に現れる電荷や電流は、真空の式では考慮に入れていた全成分の中から 真電荷および伝導電流のみを抜き出した値になっている。
  • 空間的に平均されて消えてしまう成分は、完全に省略されたためこの式には現れてこない。
  • 分極による影響は \bm D(\bm x,t) に押し込められている。
    • \bm D(\bm x,t) は見た目の電荷量（真電荷の量）から予想される電場に相当する（次元は異なるが）
    • 実際の電場 \bm E(\bm x,t) は分極や磁化の影響を含むため、 \bm D(\bm x,t) とは異なる値になる。
  • 磁化電流による影響は \bm H(\bm x,t) に押し込められている。
    • \bm H(\bm x,t) は見た目の電流量（伝導電流の量）から予想される磁束密度に相当する（次元は異なるが）
    • 実際の磁束密度 \bm B(\bm x,t) は分極や磁化の影響を含むため、 \bm H(\bm x,t) とは異なる値になる。

一般には \bm E(\bm x,t) と \bm D(\bm x,t) 、 \bm B(\bm x,t) と \bm H(\bm x,t) は完全に独立な物理量で、大きさはもちろん方向も異なる。

したがって、 \rho_e や \bm i_e が与えられただけでは変数が多すぎて解を求められない。 この方程式を解くには Maxwell 方程式に加えて、 外場に対して物質の分極や磁化がどのように応答するかの知識が必要である。

素直な系では†

実験によると、強強誘電体や強磁性体、あるいは非等方性物質を除き、 微弱な電界、磁界の元で次のような簡単な関係が成り立つことが認められる。

電即密度： \bm D(\bm x,t)=\varepsilon\bm E(\bm x,t)

磁場の強さ： \bm H(\bm x,t)=\frac{1}{\mu}\bm B(\bm x,t)

\varepsilon,\mu は単なる比例定数であるから、このとき \bm D//\bm E 、 \bm H//\bm B である。

特に真空では \varepsilon=\varepsilon_0 、 \mu=\mu_0 であるから、

\varepsilon=\varepsilon^*\varepsilon_0

\mu=\mu^*\mu_0

のように、無次元量である比誘電率 \varepsilon^* 、比透磁率 \mu^* を抜き出すこともしばしば行われる。

上で示した \bm D,\bm E および \bm H,\bm B の関係は、物性論から得られる物であり、電磁気学からの帰結ではないことに注意せよ。

比誘電率 \varepsilon^* は物質により 1\sim 10 程度で１より大きな値を取る。一般には外場の周波数や強度にも依存する。

比透磁率 \mu^* は物質により１より大きくも、小さくもなる。常磁性体では \mu^*-1= 10^{-6}\sim 10^{-4} である。反磁性体では \mu^*=10^{-5} 程度まで小さくなる。

Counter: 27492 (from 2010/06/03), today: 1, yesterday: 6