カノニカル分布の導出 のバックアップ差分(No.3)

更新


  • 追加された行はこの色です。
  • 削除された行はこの色です。
* 概要 [#e3a482b6]

納得のいく「カノニカル分布の導出法」がネット上でなかなか見つけられなかったため、自分でやってみました。

でもまだ完全には納得していません。。。

** 目次 [#ja7a6f13]
#contents

&mathjax();

* カノニカル分布の導出1:状態密度を使う [#r56d6657]

注目する系 &math(\mathrm S); (System) が温度 &math(T); の熱浴 &math(\mathrm R); (Reserver) と平衡状態にあるとします。

また、熱浴側のエネルギー &math(E_\mathrm{R}); が &math(E\le E_\mathrm{R}<E+\delta E); となるような熱浴の微視的状態の数を
&math(W_{\mathrm R}(E)\delta E); と書けるとします。すなわち &math(W_{\mathrm R}); は熱浴の状態密度です。

この熱浴に対して、

>1. 「エネルギーが &math([E, E+\delta E)); の範囲にある」という巨視的状態の持つエントロピーは~
 &math(S_{\mathrm R}(E,\delta E)=k_B\log\{W_{\mathrm R}(E)\delta E\});  (ボルツマンのエントロピー)
>
>2. 閉じた系に対する熱力学の等式:~
 &math(dE=TdS);
>
>3. 熱浴の、エネルギー &math(E); に対する確率密度関数 &math(p_\mathrm{R}); は次式で与えられる。~
 &math(
p_\mathrm{R}(E)=\frac{W_{\mathrm R}(E)}{\int_{-\infty}^\infty W_{\mathrm R}(E')dE'}\propto W_{\mathrm R}(E)
);  (等分配の法則)
);  (等確率の原理)

が成り立つことを前提としてカノニカル分布を導出します。

1., 2. から、温度 &math(T); が統計力学的な量である &math(W_\mathrm{R}(E)); により次のように定義されることが分かる。

 &math(
T(E)&\equiv\left(\frac{dS}{dE}\right)^{-1}\\
&=\left[k_B\frac{d}{dE}\Big\{\log \Big(W_\mathrm{R}(E)\delta E\Big)\Big\}\right]^{-1}\\
&=\left[k_B\frac{d}{dE}\Big\{\log W_\mathrm{R}(E)+\cancel{\log\delta E}\Big\}\right]^{-1}\\
&=\frac{1}{k_B}\left[\frac{d\log W_\mathrm{R}(E)}{dE}\right]^{-1}\\
);

このとき、小さなエネルギーの変化 &math(-\Delta E); に対して、

 &math(
\frac{1}{k_BT}=\frac{d\log W_\mathrm{R}(E)}{dE}
= \frac{\log W_{\mathrm R}(E-\Delta E)-\log W_{\mathrm R}(E)}{-\Delta E}
);

 &math(
\exp\left[-\frac{\Delta E}{k_BT}\right]=\frac{W_{\mathrm R}(E-\Delta E)}{W_{\mathrm R}(E)}
);

さらに 3. を用いると、

 &math(
p_\mathrm{R}(E-\Delta E)=p_\mathrm{R}(E)\exp\left[-\frac{\Delta E}{k_BT}\right]
);

ここで、系 &math(\mathrm S);と熱浴 &math(\mathrm R); とを合わせた全エネルギーを &math(E); とすれば、~
注目する系が &math(\Delta E); を取る確率 &math(p_{\mathrm S}(\Delta E)); は、~
熱浴が &math(E-\Delta E); を取る確率 &math(p_\mathrm{R}(E-\Delta E)); 
と等しくなります。

すなわち、

 &math(
p_{\mathrm S}(\Delta E)=p_{\mathrm S}(0)\exp\left[-\frac{\Delta E}{k_BT}\right]
);

** 「$\Delta E$ が小さい」の意味 [#gcc2e687]

「&math(\Delta E); が小さい」という上記の条件を物理学的な言葉で記述するなら、「系が熱浴とエネルギーをやりとりしても熱浴の温度が変化しない」という条件に読み替えられるため、系が熱浴に対して小さい限り一般に、

 &math(
p_{\mathrm S}(E)=p_{\mathrm S}(0)\exp\left[-\frac{E}{k_BT}\right]
);

や、

 &math(
p_{\mathrm S}(E+\Delta E)=p_{\mathrm S}(E)\exp\left[-\frac{\Delta E}{k_BT}\right]
);

が成り立つことになります。

一方、「系が熱浴とエネルギーをやりとりしても熱浴の温度が変化しない」の数学的意味は、

 &math(W_{\mathrm R}(E+\Delta E)=W_{\mathrm R}(E)\left(e^{\frac{1}{k_BT}}\right)^{\Delta E});

の右辺の &math(\left(e^{\frac{1}{k_BT}}\right)); が定数ということなので、
「熱浴の状態密度が &math(E); の近傍 &math(\Delta E); でエネルギーとともにほぼ指数関数的に増加する」
であると理解できます。

** 熱平衡に於いて系全体で温度が等しくなる理由 [#d2315946]

注目する系 &math(\mathrm S); のエネルギーが &math(E\sim E+\delta E); の範囲に入るとき、
熱浴 &math(\mathrm R); はエネルギー &math(E_\mathrm{total}-E-\delta E\sim E_\mathrm{total}-E); 
の範囲に入るから、そのような巨視的状態を実現する系全体の状態数は &math(\mathrm S,\mathrm R); 
の状態数を掛け合わせて、

 &math(W_\mathrm{S}(E)\delta E\cdot W_\mathrm{R}(E_\mathrm{total}-E)\delta E);

と表せます。

等確率の原理により、熱平衡に於いてはこの状態数を最大化するように系のエネルギー &math(E); 
が分配されます。すなわち &math(E); に対する微係数はゼロで、

 &math(\frac{d}{dE}\Big\{W_\mathrm{S}(E)\cdot W_\mathrm{R}(E_\mathrm{total}-E)\Big\}=W'_\mathrm{S}(E)\cdot W_\mathrm{R}(E_\mathrm{total}-E)-W_\mathrm{S}(E)\cdot W'_\mathrm{R}(E_\mathrm{total}-E)=0);

 &math(W'_\mathrm{S}(E)\cdot W_\mathrm{R}(E_\mathrm{total}-E)=W_\mathrm{S}(E)\cdot W'_\mathrm{R}(E_\mathrm{total}-E));

 &math(\frac{W'_\mathrm{S}(E)}{W_\mathrm{S}(E)}=\frac{W'_\mathrm{R}(E_\mathrm{total}-E)}{W_\mathrm{R}(E_\mathrm{total}-E)});

 &math(\left .\frac{d\log W_\mathrm{S}}{dE_\mathrm{S}}\right |_E=
\left .\frac{d\log W_\mathrm{R}}{dE_\mathrm{R}}\right |_{E_\mathrm{total}-E});

両辺に &math(k_B); をかけて逆数を取れば、

 &math(T_\mathrm{S}(E_\mathrm{S})=T_\mathrm{R}(E_\mathrm{R}));

が熱平衡の条件であることが分かります。

同様の議論を、全系を &math(\mathrm{S1},\mathrm{S2},\mathrm R); の三つに分けて行えば、
&math(E_1,E_2); を変数として

 &math(W_\mathrm{S1}(E_1)\cdot W_\mathrm{S2}(E_2)\cdot W_\mathrm{R}(E_\mathrm{total}-E1-E2));

を最大化する条件から、

 &math(T_\mathrm{S1}(E_\mathrm{1})=T_\mathrm{S2}(E_\mathrm{2})=T_\mathrm{R}(E_\mathrm{R}));

が得られ、系の各部がすべて同じ温度を持つことが熱平衡の条件であることが分かります。

すなわち、温度は示強性の状態変数です。

** カノニカル分布はどこから出てきたか [#z19e7fa8]

上記の議論を見返すと、冒頭の前提条件のうち 1., 2. についてはそれぞれ、
「エントロピー」、「温度」の用語を定義しているに過ぎず、
物理的な「原理」を表すものではありません。

すると、「実現確率は状態密度に比例する」
「熱平衡において巨視的状態の持つ微視的状態数は最大化される」なども含めて、
すべての源流は「等確率の原理」であることに気づきます。

話の流れは、

「等確率の原理」~
  ↓~
「実現確率は状態密度に比例する」~
  ↓~
「熱平衡において巨視的状態の持つ微視的状態数は最大化される」~
  ↓~
「熱平衡において系の各部は同じ『温度』を持つ」~
  ↓~
「温度は示強性の状態変数である」~
  ↓~
「温度とエントロピーの定義から『熱力学の等式』は自明」~
  ↓~
「カノニカル分布が導かれる」

となっていることを読み取れるはずです。

** 無数の調和振動子からなる系 [#x7b12d10]

https://twitter.com/kikumaco/status/879971741845434368

を参考にして、同じ固有振動数 &math(\omega); を持つ &math(M); 個の調和振動子からなる熱浴を考える。

全エネルギーを &math(E=N\hbar\omega); とすると、それぞれの調和振動子へエネルギーを分配する方法の数は &math(N); 個の玉の間に &math(M); 本の区切りを入れる方法の数と等しいため(正確には &math(M-1); 本であるが &math(M\gg 1); なので同一視する)、

 &math({}_{M+N}C_M=\frac{(M+N)!}{M!N!});

この状態数が &math(\hbar\omega); 毎に存在するため、状態密度は

 &math(W(E)=\frac{1}{\hbar\omega}\frac{(M+N)!}{M!N!}); ただし &math(N(E)=E/\hbar\omega);

「エネルギーが &math(E\sim E+\delta E); ただし &math(\delta E\gg\hbar\omega); 
の範囲にある」という状態を巨視的状態とすると、対応するボルツマンのエントロピーは、
スターリングの公式を使って

 &math(
S(E)=k_B\log W(E)\delta E=(M+N)k_B\log(M+N) - Mk_B\log M-Nk_B\log N+k_B\log\frac{\delta E}{\hbar\omega}
);

温度は、

 &math(
T(E)&=\left(\frac{dS}{dE}\right)^{-1}=\left(\frac{1}{\hbar\omega}\frac{dS}{ dN}\right)^{-1}\\
&=\frac{\hbar\omega}{k_B}\bigl\{\log(M+N)-\log N\bigr\}^{-1}\\
&=\frac{\hbar\omega}{k_B}\bigl\{\log\frac{M+N}{N}\bigr\}^{-1}\\
);

この熱浴と相互作用する相対的に小さな系が熱浴からエネルギー &math(\varepsilon=n\hbar\omega); 
ただし &math(n\ll N); だけ奪うと、

 &math(
W(E-\varepsilon)&=\frac{1}{\hbar\omega}\frac{(M+N-n)!}{M!(N-n)!}\\
&\sim\left(\frac{N}{M+N}\right)^n\frac{1}{\hbar\omega}\frac{(M+N)!}{M!N!}\\
&=\left(\frac{N}{M+N}\right)^nW(E)\\
);

となるが、両辺の対数を取ると、

 &math(
\log W(E-\varepsilon)=n\log \frac{N}{M+N}+\log W(E)
);

 &math(
k_B\log W(E-\varepsilon)-k_B\log W(E)=(n\hbar\omega) \frac{k_B}{\hbar\omega}\log \frac{N}{M+N}
);

 &math(
\frac{k_B\log W(E-\varepsilon)-k_B\log W(E)}{-\varepsilon}=\frac{dS}{dE}=-\frac{k_B}{\hbar\omega}\log \frac{N}{M+N}=\frac{k_B}{\hbar\omega}\log \frac{M+N}{N}=1/T
);

となって、&math(n\ll N); の条件は、エントロピーの増加が線形と見なせる範囲であるという条件、
また、&math(N-n); を &math(N); と同一視できるという条件、となっていることが分かる。

** 納得のいかない点 [#ec74646e]

完全に離散的なままだと「エネルギーによる微分」が扱えないので、
離散的な状態を微分可能な程度にならした「状態密度」を導入して議論したのですが、
もっと離散的な系に密接に対応づけられないのかどうか、
まだ完全には納得がいっていません。
このような議論だと巨視的状態には &math(\hbar\omega); よりずっと大きい &math(\delta E); 
のエネルギー幅を持たせなければならなくなってしまい、

注目する系 &math(\mathrm S); がある正確なエネルギー &math(E_k); を取る確率、

を議論できません。

離散的記述の中に「状態密度」あるいは「エントロピー」がエネルギーで微分可能であることをどううまく入れ込むべきか、まだ完全には納得がいっていません。

* 質問・コメント [#l774a821]

#article_kcaptcha


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