エネルギー保存・運動量保存 のバックアップ(No.6)

更新


電磁気学

点電荷の系

電磁ポテンシャルの導入 では、 \bm i(\bm x,t),\rho(\bm x,t) を与えて、そこから生まれる \bm E(\bm x,t),\bm B(\bm x,t) を求める問題について考えた。

ここでは、 \bm E(\bm x,t),\bm B(\bm x,t) ばかりでなく \bm i(\bm x,t),\rho(\bm x,t) も未知であり 電荷が電磁場と相互作用しながら共に時間発展していく様子を記述することを考える。

系は n 個の点電荷からなり、それぞれ

  • 電荷 e_k
  • 質量 m_k
  • 位置座標 \bm z_k(t)

を持つとする。

電荷密度

1つの点電荷 e が作る電荷密度を考えよう。

\bm x=\bm z(t) に電荷があり、 それ以外で電荷密度はゼロであるから、 電荷密度は Dirac のデルタ関数を用いて、

\rho(\bm x,t)=e\delta^3(\bm x-\bm z(t))

と表せる。この式から時間と共に点電荷が移動し、 電荷密度が変化する様子を思い浮かべてみよう。

ここで、 \delta^3(\bm x)=\delta(x)\delta(y)\delta(z) であり、 \delta(x)

\delta(x)=\begin{cases}0&x\ne 0\\+\infty&x=0\end{cases}

\delta(-x)=\delta(x)

\int_{-\infty}^\infty\delta(x)dx=1

を満たす。

電流密度

点電荷 e の電流密度は次のように書ける。

\bm i(\bm x,t)=e\dot{\bm z}(t)\delta^3(\bm x-\bm z(t))

と言われても、

電流密度が点電荷の存在する位置 \bm z(t) 以外でゼロになること、
大きさが電荷 e およびその速度 \dot{\bm z}(t) に比例すること、
は納得できるが、他に係数が必要ないのかちょっと心配になるかもしれない。

この点は、電荷保存則が成り立っていることを確認すれば納得できる。

&math( \frac{\PD}{\PD t}\rho(\bm x,t) &=-e\dot{\bm z}(t)\cdot\GRAD\delta^3(\bm x-\bm z(t))\\ &=-e\sum_{a=x,y,z}\dot z_a(t)\frac{\PD}{\PD a}\delta^3(\bm x-\bm z(t)) );

&math( \DIV \bm i(\bm x,t) &=e\,\DIV\{\dot{\bm z}(t)\delta^3(\bm x-\bm z(t))\}\\ &=\sum_{a=x,y,z}\frac{\PD}{\PD a}\dot z_a(t)\delta^3(\bm x-\bm z(t))\\ &=\sum_{a=x,y,z}\dot z_a(t)\frac{\PD}{\PD a}\delta^3(\bm x-\bm z(t)) );

&math( \therefore \frac{\PD}{\PD t}\rho(\bm x,t)+\DIV \bm i(\bm x,t)=0 );

として、上記の電流密度の式が正しいことが分かる。

ローレンツ力

点電荷 e に働く力は、

&math( \bm f(t)&=m\ddot{\bm z}(t)\\ &=\int\Big[\rho(\bm x,t)\bm E(\bm x,t)+\bm i(\bm x,t)\times\bm B(\bm x,t)\Big]d^3x\\ &=\int\Big[e\delta^3(\bm x-\bm z(t))\bm E(\bm x,t)+e\delta^3(\bm x-\bm z(t))\dot{\bm z}(t)\times\bm B(\bm x,t)\Big]d^3x\\ &=e\bm E(\bm z(t),t)+e\dot{\bm z}(t)\times\bm B(\bm z(t),t) );

と表せる。

解くべき方程式

以上を合わせると、

  • Maxwell 方程式
  • 電荷密度、電流密度
  • 運動方程式

をすべて連立させることで、系の時間発展を記述できる。

&math( \left\{\begin{array}{l} \displaystyle\ROT\bm E(\bm x,t)+\frac{\PD}{\PD t}\bm B(\bm x,t)=\bm 0\vspace{2mm}\\ \displaystyle\DIV\bm B(\bm x,t)=0\vspace{2mm}\\ \displaystyle\frac{1}{\mu_0}\ROT \bm B(\bm x,t)-\varepsilon_0\frac{\PD}{\PD t}\bm E(\bm x,t)=\bm i(\bm x,t)\vspace{2mm}\\ \displaystyle\varepsilon_0\DIV\bm E(\bm x,t)=\rho(\bm x,t)\vspace{2mm}\\ \displaystyle\rho(\bm x,t)=\sum_{k=1}^n e_k\delta^3(\bm x-\bm z_k(t))\vspace{2mm}\\ \displaystyle\bm i(\bm x,t)=\sum_{k=1}^n e_k\dot{\bm z}_k(t)\delta^3(\bm x-\bm z_k(t))\vspace{2mm}\\ m_1\ddot{\bm z}_1(t)=e_1\Big[\bm E(\bm z_1(t),t)+\dot{\bm z}_1(t)\times\bm B(\bm z_1(t),t)\Big]\vspace{2mm}\\ m_2\ddot{\bm z}_2(t)=e_2\Big[\bm E(\bm z_2(t),t)+\dot{\bm z}_2(t)\times\bm B(\bm z_2(t),t)\Big]\vspace{2mm}\\ m_3\ddot{\bm z}_3(t)=e_3\Big[\bm E(\bm z_3(t),t)+\dot{\bm z}_3(t)\times\bm B(\bm z_3(t),t)\Big]\vspace{2mm}\\ \vdots\\ \end{array}\right. );

実際にこれを解くのは難しいが・・・

ガリレイの相対性について

ニュートン方程式では、

  • 静止した人から見た物理現象と、
  • 等速直線運動をしている人から見た物理現象は、

どちらも同じ方程式、つまりニュートン方程式を満たすことになる。
このことをガリレイの相対性と呼ぶ。
→ 絶対的な静止系は存在せず、静止・運動は相対的であるということ

上記の方程式に対して同様の相対性が成り立つだろうか?

ある座標で見て電荷が一点で静止していれば、その周りには球対象な電場のみが存在する。

同じ状況を等速直線運動する人から見れば、 電荷は移動しており、その周りには電場だけではなく磁場も発生しているはずである。

このことは、磁場の大きさ・方向、あるいは有無自体が観測者の運動に依存することを 示しており、直感的な意味でのガリレイの相対性が失われていることを示している。

このあたりをしっかり考えていくことで、特殊相対性理論にたどり着くことになる。

エネルギー保存則

遠隔相互作用では:

  • 静電ポテンシャル \frac{1}{4\pi\varepsilon_0}\frac{e_1e_2}{r^2}
  • コイルの電流 \frac{L}{2}I^2
  • コンデンサの電荷 \frac{Q^2}{2C}
  • ジュール熱 \int_0^\infty IV\,dt
  • 光の運ぶエネルギー (太陽定数)=1366\,\mathrm{W/m^2}

というように、個々の問題に対して別々の方法でエネルギーを測っていた。

近接相互作用の考え方で:

  • 電磁場の存在する空間の持つエネルギー
  • 電磁場により運ばれるエネルギー

を、局所的な電磁場の値のみによって定義しようというのが以下の目標である。

これらは同じエネルギーを別々の形で表わそうという物なので、 ある意味排他的である。すなわち、両方でエネルギーを記述してしまうと、 同じエネルギーを2度測ってしまうことになる。

基本

エネルギー保存則を式で書けば:

  \frac{d}{dt}(エネルギー)=0

あるいは、

  \frac{d}{dt}(エネルギー)=(単位時間あたりに入ってくるエネルギー)

となるはず。

一般に、このような式は Newton 方程式に速度を掛けることで得られる。

&math( &\bm f=m\ddot{\bm z}\\ &\bm f\cdot\dot{\bm z}=m\ddot{\bm z}\cdot\dot{\bm z}= \frac{d}{dt}\Big(\frac{1}{2}m|\dot{\bm z}|^2\Big)\\ );

この左辺 \bm f\cdot\dot{\bm z} は系に「加えられた仕事」を表わし、
右辺は運動エネルギー \frac{1}{2}m|\dot{\bm z}|^2 の時間変化を表わす。

Lorentz 力に適用する

&math( &m_k\dot{\bm z}_k(t)\cdot\ddot{\bm z}_k(t)=\frac{d}{dt}\Big(\frac{1}{2}m_k|\dot{\bm z}_k(t)|^2\Big)\\ &=e_k\dot{\bm z}_k(t)\cdot\bm E(\bm z_k(t),t)+e_k\underbrace{\dot{\bm z}_k(t)\cdot\big[\dot{\bm z}_k(t)\times\bm B(\bm z_k(t),t)\big]}_{\displaystyle \bm a\cdot(\bm a\times\bm b) の形なのでゼロ} );

すべての粒子について足し合わせると、

&math( &\frac{d}{dt}\left(\sum_{k=1}^n\frac{1}{2}m_k|\dot{\bm z}_k(t)|^2\right)\\ &=\sum_{k=1}^n e_k\dot{\bm z}_k(t)\cdot\bm E(\bm z_k(t),t)\\ &=\sum_{k=1}^ne_k\dot{\bm z}_k(t)\cdot\int d^3x\ \delta^3(\bm x-\bm z_k(t))\bm E(\bm x,t)\\ &=\int d^3x\ \left(\sum_{k=1}^ne_k\dot{\bm z}_k(t)\delta^3(\bm x-\bm z_k(t))\right)\cdot\bm E(\bm x,t)\\ &=\int d^3x\ \bm i(\bm x,t)\cdot\bm E(\bm x,t) );

Maxwell 方程式を用いて \bm i(\bm x,t) を書き直すと、

&math( &=\int d^3x\ \left(\frac{1}{\mu_0}\ROT B-\varepsilon_0\frac{\PD\bm E}{\PD t}\right)\cdot\bm E(\bm x,t)\\ );

Maxwell 方程式より \frac{1}{\mu_0}\bm B\cdot\underbrace{\left(\frac{\PD\bm B}{\PD t}+\ROT\bm E\right)}_{=\bm 0}=0 を加えると、

&math( &=\int d^3x\ \left[\frac{1}{\mu_0}\ROT\bm B\cdot\bm E-\frac{1}{\mu_0}\bm B\cdot\ROT\bm E-\frac{1}{\mu_0}\bm B\cdot\frac{\PD\bm B}{\PD t}-\varepsilon_0\bm E\cdot\frac{\PD\bm E}{\PD t}\right]\\ );

ベクトル解析の公式から、 \DIV(\bm E\times\bm B)=\bm B\cdot\ROT\bm E-\bm E\cdot\ROT\bm B なので、

&math( &=\int d^3x\ \left[

  • \frac{1}{\mu_0}\DIV(\bm E\times\bm B)
  • \frac{\PD}{\PD t}\left\{\frac{1}{\mu_0}|\bm B|^2-\varepsilon_0|\bm E|^2\right\}\right]\\ &=\int_S \underbrace{\frac{1}{\mu_0}(\bm E\times\bm B)}_{=\,\bm S(\bm x,t)}\cdot\bm ndS
  • \frac{d}{dt}\int d^3x\ \underbrace{\left\{\frac{1}{\mu_0}|\bm B|^2-\varepsilon_0|\bm E|^2\right\}}_{=\,w(\bm x,t)}\\ );

結果として、

&math( \frac{d}{dt}\Bigg[\underbrace{\sum_{k=1}^n\frac{1}{2}m_k|\dot{\bm z}_k(t)|^2}_{全運動エネルギー}+\underbrace{\rule[-2.8ex]{0pt}{1ex}\int d^3x\,w(\bm x,t)}_{電磁場のエネルギー}\Bigg]= \underbrace{\rule[-2.8ex]{0pt}{0pt}\int_S\bm S(\bm x,t)\cdot\bm ndS}_{界面から入ってくる量} );

の形を得る。

ここで、

  • 電磁場のエネルギー密度
    &math( w(\bm x,t)&=\frac{1}{2}\bigg[\varepsilon_0|\bm E|^2+\frac{1}{\mu_0}|\bm B|^2\bigg]\\ &=\frac{1}{2}\bigg[\bm E\cdot\bm D+\bm B\cdot\bm H\bigg] );
  • 電磁場によるエネルギーフラックス密度 (ポインティングベクトル)
    &math( \bm S(\bm x,t)&=\frac{1}{\mu_0}\bm E\times\bm B\\ &=\bm E\times\bm H );

ただし、 \bm D=\varepsilon\bm E,\bm H=\frac{1}{\mu_0}B を導入することにより係数を消した。

energy.png

運動量保存則

光が運動量を運ぶことは有名であり、これを利用した光子ヨットなども検討されている。
すなわち、電磁場はそれ自体が運動量を持ち、その値は物質系や外界との相互作用で刻々と変化する。

そこで、

 (系の全運動量)=(質点系の運動量)+(電磁場の運動量)

と考え、全運動量の時間変化について考えることにしよう。

基本

ニュートン方程式を

\bm f=m\ddot{\bm z}=\frac{d}{dt}\underbrace{\Big(\,m\dot{\bm z}\,\Big)}_{運動量}

のように変形すると、 運動量の時間変化は外力に等しいことを示している。

ローレンツ力の式から、系の全運動量の時間変化を考えよう。

&math( &\sum_{k}m_k\ddot{\bm z}_k(t)=\frac{d}{dt}\Big(\sum_km_k\dot{\bm z}_k(t)\Big)\\ &=\sum_k \Big[e_k\bm E(\bm z_k,t)+e_k\bm \dot{\bm z}_k(t)\times\bm B(\bm z_k,t)\Big]\\ &=\sum_k\int d^3x \Big[e_k\delta(\bm x-\bm z_k(t))\bm E(\bm x,t)+e_k\bm \dot{\bm z}_k(t)\times\delta(\bm x-\bm z_k(t))\bm B(\bm x,t)\Big]\\ &=\int d^3x \Big[\Big\{\sum_ke_k\delta(\bm x-\bm z_k(t))\Big\}\bm E(\bm x,t)+\Big\{\sum_ke_k\bm \dot{\bm z}_k(t)\delta(\bm x-\bm z_k(t))\Big\}\times\bm B(\bm x,t)\Big]\\ &=\int d^3x \Big[\rho(\bm x,t)\bm E(\bm x,t)+\bm i(\bm x,t)\times\bm B(\bm x,t)\Big]\\ );

Maxwell 方程式を用いて \rho,\bm i を書き換えると、

&math( &=\int d^3x \Big[\Big(\varepsilon_0\DIV \bm E\Big)\bm E+\Big(\frac{1}{\mu_0}\ROT\bm B-\varepsilon_0\frac{\PD \bm E}{\PD t}\Big)\times\bm B\Big]\\ );

さらに、

&math( \frac{\PD}{\PD t}\Big(\varepsilon_0\bm E\times\bm B\Big) &=\varepsilon_0\frac{\PD\bm E}{\PD t}\times \bm B+\varepsilon_0\bm E\times\frac{\PD\bm B}{\PD t}\\ &=\varepsilon_0\frac{\PD\bm E}{\PD t}\times \bm B-\varepsilon_0\bm E\times\ROT \bm E\\ );

を積分して両辺に加えれば、

&math( &\frac{d}{dt}\Big(全運動量\Big)+\frac{d}{dt}\Big[\underbrace{1/c^2\rule[-2.4ex]{0pt}{1ex}}_{\displaystyle=\varepsilon_0\mu_0}\int d^3x\underbrace{\Big(\frac{1}{\mu_0}\bm E\times\bm B\Big)}_{\displaystyle =\bm S}\Big]\\ &=\int d^3x\ \Big[\varepsilon_0(\DIV \bm E)\bm E+\varepsilon_0(\ROT \bm E)\times\bm E+\\ &\hspace{18mm}\underbrace{\frac{1}{\mu_0}(\DIV \bm B)\bm B}_{\displaystyle=0}+ \frac{1}{\mu_0}(\ROT \bm B)\times\bm B\Big] );

を得る。ただし、式の対称性を高めるために常にゼロとなる項 \frac{1}{\mu_0}(\DIV \bm B)\bm B を右辺に加えた。

右辺は対称性は良い物の、まだごちゃごちゃしていて分かりづらい。ここを何とかしよう。

$(\mathrm{div}\,\bm X)\bm X+(\mathrm{rot}\,\bm X)\times\bm X$ という量について

次のような2階のテンソル量を考える。
(2階のテンソルとは、ベクトルと同様に、座標変換に対して決まった法則に従い形を変える行列量のことである)

&math( T^e&=\varepsilon_0\begin{pmatrix} E_x^2-\frac{1}{2}|\bm E|^2&E_xE_y&E_xE_z\\ E_yE_x&E_y^2-\frac{1}{2}|\bm E|^2&E_yE_z\\ E_zE_x&E_zE_y&E_z^2-\frac{1}{2}|\bm E|^2\\ \end{pmatrix}\\ &=\varepsilon_0\begin{pmatrix} E_x^2&E_xE_y&E_xE_z\\ E_yE_x&E_y^2&E_yE_z\\ E_zE_x&E_zE_y&E_z^2\\ \end{pmatrix}-\frac{\varepsilon_0}{2}|\bm E|^2I\\ &=\varepsilon_0\Big[

 \begin{pmatrix}E_x\\E_y\\E_z\end{pmatrix}
 \begin{pmatrix}E_x&E_y&E_z\end{pmatrix}
  • \frac{1}{2}\begin{pmatrix}E_x&E_y&E_z\end{pmatrix}\begin{pmatrix}E_x\\E_y\\E_z\end{pmatrix}I\Big]\\ &=\varepsilon_0\Big[\bm E\bm E^T-\frac{1}{2}(\bm E^T\bm E) I\Big]\\ );

直交行列 R により \bm E\to \bm E'=R\bm E の形の座標変換が施されるとき、 T^e の変換規則は次のようになり、この量が2階のテンソル量であることを確認できる。

&math( T^e\to {T^e}'&=\varepsilon_0\Big[(R\bm E)(R\bm E)^T-\frac{1}{2}\{(R\bm E)^T(R\bm E)\} I\Big]\\ &=\varepsilon_0\Big[R\bm E\bm E^TR^T-\frac{1}{2}(\bm E^T\underbrace{R^TR}_{=I}\bm E) \underbrace{R^TR}_{=I}\Big]\\ &=\varepsilon_0R\Big[\bm E\bm E^T-\frac{1}{2}(\bm E^T\bm E) I\Big]R^T\\ &=R\,T^eR^T\\ );

ここで、

&math(\DIV \bm E=\begin{pmatrix}\frac{\PD}{\PD x}&\frac{\PD}{\PD y}&\frac{\PD}{\PD z}\end{pmatrix} \begin{pmatrix}E_x\\E_y\\E_z\end{pmatrix});

であったのと同様に、

&math( \DIV T^e &=\begin{pmatrix}\frac{\PD}{\PD x}&\frac{\PD}{\PD y}&\frac{\PD}{\PD z}\end{pmatrix} \begin{pmatrix}T^e_{xx}&&T^e_{xy}&&T^e_{xz}\\T^e_{yx}&&T^e_{yy}&&T^e_{yz}\\T^e_{zx}&&T^e_{zy}&&T^e_{zz}\\\end{pmatrix}\\ &=\begin{pmatrix}\rule[-4ex]{0pt}{9ex} \big(\frac{\PD T^e_{xx}}{\PD x}+\frac{\PD T^e_{yx}}{\PD y}+\frac{\PD T^e_{zx}}{\PD z}\big)& \big(\frac{\PD T^e_{xx}}{\PD x}+\frac{\PD T^e_{yx}}{\PD y}+\frac{\PD T^e_{zx}}{\PD z}\big)& \big(\frac{\PD T^e_{xx}}{\PD x}+\frac{\PD T^e_{yx}}{\PD y}+\frac{\PD T^e_{zx}}{\PD z}\big) \end{pmatrix} );

を定義すると、その成分を計算することにより

&math( \DIV T^e=\varepsilon_0(\DIV\bm E)\bm E+\varepsilon_0(\ROT\bm E)\times\bm E );

となることを確かめられる。

同様に、

&math( T^m=\frac{1}{\mu_0}\Big[\bm B\bm B^T-\frac{1}{2}(\bm B^T\bm B) I\Big]\\ );

と置く時、 T=T^e+T^m を Maxwell の応力テンソルと呼ぶ。

電磁場の運動量・電磁場の及ぼす応力

この T を用いることにより、

&math( \frac{d}{dt}\Big[\underbrace{\sum_km_k\dot{\bm z}_k(t)}_{粒子の運動量}+ \int d^3x\underbrace{\rule[-3ex]{0pt}{1pt}\Big\{\frac{1}{c^2}\bm S(\bm x,t)\Big\}}_{運動量密度}\Big]\\ =\int d^3x\,\DIV T(\bm x,t)=\underbrace{\int_S T(\bm x,t)\bm ndS\rule[-3ex]{0pt}{1pt}}_{\displaystyle=\bm F(t)} );

のように書き表せる。

  • 左辺2項目は電磁場の持つ運動量で、その単位体積あたりの密度は
    ポインティングベクトルを使って \bm S(\bm x,t)/c^2 と表せる
  • 右辺は外界から系に加わる力で、界面部分の電磁場の値のみを用いて書かれている
  • 外界から系に加わる力は必ずしも界面に垂直とは限らず、 界面に水平な成分も持っている。これは、 T\bm n が必ずしも \bm n と平行にならないことから分かる。

図で表わすと

moment.png

すべての質点(点電荷)を含む閉曲面 S と、 S に囲まれる体積 V を考える。

全質点の運動量の他に、電磁場の存在する空間自体が運動量を持つというのが上での結論。

微小体積 d^3x の持つ運動量は、体積に運動量密度を掛けて d\bm p=\frac{1}{c^2}\bm S(\bm x,t)d^3x と表せる。

系に及ぼされる外力は界面 S における Maxwell 応力を積分することにより求められる。

微小面積 dS にかかる応力は、法線ベクトルに応力テンソル T を掛けて求められ、 d\bm f=T\bm ndS と表せる。

この描像では個々の質点がうけるローレンツ力は、全運動量の変化に影響しない物として無視できる。

近接相互作用の視点から見れば、個々の質点が感じるローレンツ力は、 まさに質点の存在する点に存在する電磁場との相互作用であり、 系内部での運動量のやりとりにしかならない(全運動量に寄与しない)。 これに対し、系外部との運動量のやりとりは界面 S 上でのみで起きていることになる。

例:一様な電磁場

\bm E=(0\ 0\ E_z) のとき、

&math(T=T^e=\varepsilon_0\begin{pmatrix}

  • \frac{1}{2}E_z^2&&\\ &-\frac{1}{2}E_z^2&\\ &&\frac{1}{2}E_z^2 \end{pmatrix});

となる。

maxwell-stress.png

面が電場に垂直である場合:

\bm n=(0\ 0\ \pm1)\parallel\bm E の時、 T\bm n=\frac{\varepsilon_0}{2}E_z^2\bm n

力は常に \bm n と同じ向きになるため、 微小領域における電場に垂直な界面には引っ張り(外向き)の応力がかかることになる。

面が電場に平行である場合:

\bm n=(\pm 1\ 0\ 0)\perp\bm E あるいは \bm n=(0\ \pm 1\ 0)\perp\bm E の時、 T\bm n=-\frac{\varepsilon_0}{2}E_z^2\bm n

力は常に \bm n と逆の向きになるため、 微小領域における電場に平行な界面には圧縮(内向き)の応力がかかることになる。

空間にそのような影響を与える理由として、電気力線には長さ方向に縮もうとするほうこうに、 また、線間をなるべく広く取ろうとする方向に力を生じる性質があるとする古典的な解釈が 役に立つ。

electrical flux line.png

電気力線が線に沿って縮もうとすると空間は伸ばされ、戦艦を広げようとすると空間は圧縮される、 作用・反作用の関係になっている。

質問・コメント





Counter: 37490 (from 2010/06/03), today: 2, yesterday: 0