第一章 なぜ永久機関は実現不可能なのか
指数分布の理由
2006/08/21  

前節で、エネルギーの最も自然な分布は、ボルツマン分布

Ni 比例 exp(−Ei/kT)
に従うということを述べた。 なぜエネルギーの分布はこのような指数関数となるのだろうか。

まず、問題を確認しよう。
「多数の分子にエネルギーを配るときに、最も確率が高くなるのはどのように配分したときか」

最初に、問題をうんと単純化して分子が2個の場合を考える。 分配すべき全エネルギーをEとしよう。 このとき、1個の分子の持つエネルギーがE1となる確率をPr(E1)とすれば

Pr(E1) = 1/E
となるであろう。 なぜなら、分子の持つエネルギーは0〜Eの間で均一だからである。

なぜ指数分布2次元

次に、分子が3個あった場合を考える。 3つの分子の持つエネルギーを E1, E2, E3 として、この3つをそれぞれX,Y,Z軸にあてはめた(立体の)グラフを書いてみよう。 全エネルギー一定という条件は、グラフ上の正三角形の面として表現されている。 E1の軸に着目して三角形を見ると、E1がEに近づくほど、それに比例して三角形の切り口の長さが短くなる、つまり確率が下がることがわかる。 三角形をE1軸に投影した面積(三角形の切り口の長さの積分)は E^2 / 2 であることを考慮すると、求める確率は

Pr(E1) = (E - E1) / 三角形の投影面積 = (E - E1) * 2 / E^2
となる。

なぜ指数分布3次元

分子が4個になると、図示するのは難しい。(4次元なので) 上の分子3個で用いた様な図を、4個目の分子が持つエネルギーの大きさ順に並べて想像するしかない。 ここでは、着目する分子のエネルギーE1の順に図が並んでいるものと考える。 全エネルギー一定という条件は、今度は個々の図の正三角形を重ね合わせた三角錐の体積として表現されている。 エネルギーがE1であるような切り口の図に着目すると、E1がEに近づくほど、三角形の面積が小さくなる。 面積は長さの2乗に比例するので、確率は (E - E1) ^ 2 に比例して小さくなる。 三角錐の体積は E^3 / 3 となるので、求める確率は

Pr(E1) = (E - E1) ^ 2 / 三角錐の体積 = (E - E1) ^ 2 * 3 / E^3

なぜ指数分布4次元

ここまで来ると、分子がもっと増えたときの確率の予想が付くだろう。 分子がn+1個あった場合の確率は

Pr(E1) = (E - E1) ^ (n - 1) * n / E^n
となる。
(直線の長さ -> 三角形の面積 -> 三角錐の体積 -> 4次元三角錐の体積(?) -> ... は、N次の積分を考えればよい)
Pr(E1)
= (E - E1) ^ (n - 1) * n / E^n
= (E - E1) ^ (n - 1) * 1 / E^(n-1) * n/E  -- 後の項を分けて(n-1)という因子を出した
= (1 - E1/E) ^ (n - 1) * n/E
ここで、分子数が非常に大きくなったとき lim(n→∞) を考える。 自然対数 e の定義から考えると ※
lim(n→∞) [ (1 - E1/E) ^ (n - 1) * n/E ]
= Exp( - E1 / E ) * lim(n→∞) n/E
ここで全エネルギーを分子数で割った平均値、E/n を A と置くと、
lim(n→∞) n/E = 1/A
これを用いて
上式 = Exp( - E1 / E ) * 1/A
結果はエネルギーE1 に対する指数関数となっていることが確認された。

以上が指数分布の理由なのだが、それでも実感が湧かないという方のために数式抜きの説明を付け足そう。 例えば人生をハードル越えの連続と見なして、生まれ落ちたときには全員が等しくスタート地点に並んでいるものだとしよう。 仮にスタート地点に1万人が並んでいたとする。 1つ目のハードルを越えるときに、そのうちの半数が脱落するものとしよう。 無事ハードルを越えた人数は5000人である。 2個目のハードルを越えるときにも、残りの半数が脱落するものとする。 すると、2個目を越えた人数は2500人となる。 3個目で1250人、4個目で625人、5個目で312(又は313人、割り切れないので)、以下半分の人数、半分の人数が残ることになる。 この様子を天から眺めれば、1つ目のハードルに残っている人数が最も多く、ハードルの数が増すにつれて徐々に人数が減っていく分布が見えることになる。 これがまさしく指数分布である。 ハードルの高さは、いわばエネルギーの段差に相当する。 人生のたとえを借りるなら、ハードル一段の違いは持てる財産の違いのようなものである。 (何とも切ない話ではあるが。) 上の例では1つのハードルで半数が脱落するものとしたが、ここで脱落者の割合を80%、90%といった具合に低く抑え、その代わりハードルの数をたくさんに増やす。 その極限として、スタート地点の人数が最も多く、スタートから離れるに従って人数が漸次減少するような分布となるのである。

※ 自然対数の定義は・・・
e = lim(n→∞) (1 + 1/n) ^ n

自然対数を底とする指数関数は
Exp( a )
 = lim(n→∞) (1 + 1/n) ^ N
ここで an = N と置いて、N→∞ を考えると
 = lim(N→∞) (1 + 1/(N/a)) ^ N
 = lim(N→∞) (1 + a/N) ^ N
この式の形と上の式を見比べてみよう。


ここに述べた内容は、いわゆる教科書的な導出ではなく、直感的に把握するための方法である。 大抵の教科書には「ラグランジュの未定乗数法」と「スターリングの公式」を用いた導出が示されている。 私も教科書流の方が正道とは思う。 ただ、私自身のことを申し上げると、教科書の数式を追ってみて正しいことは分かっても、それで納得した気にはなれなかった。 それゆえ、ここでは私自身が納得した直感的な方法を挙げることにした。 興味のある方は、正統な統計力学の教科書を参照して欲しい。
ページ先頭に戻る▲