tecktom.com

この章のまとめ

2026年10月15日9分で読めますAI生成教科書『虎の巻』第3部 第2章 等重率とミクロカノニカル

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第3部 第2章 §3.2.7 ★★ 等重率=相空間の体積で一様、状態数=体積を h3NN!h^{3N}N! で割った数、エントロピー=その対数、温度=そのエネルギー微分の逆数

孤立系について分かっているのはエネルギー EE(と VV, NN)だけで、相空間のどの点にいるかは分からない。その分からなさを、エネルギーの殻 E<H≤E+ΔEE<H\le E+\Delta E の中で相空間の体積 dΓ\mathrm{d}\Gamma について一様な分布(ミクロカノニカル分布)で表すのが等重率の原理である。dΓ\mathrm{d}\Gamma で一様に選ぶ理由は Liouville の定理で、dΓ\mathrm{d}\Gamma が流れで保たれるので、殻の中で一様な分布は時間で変わらない。自由度1では殻の上の測度が dt\mathrm{d} t に一致し、殻の平均は時間平均と厳密に等しい。自由度が2以上では、これはエルゴード性という力学の性質にかかっていて、一般には示されていないので、統計力学は等重率を原理として置く。

数えるます目の大きさは、1つの量子状態が相空間で占める体積 hh(自由度あたり)で、同じ種類の粒子を入れ替えた状態を数え直さないために N!N! で割る。こうして数えた状態数の対数に kBk_{\mathrm B} を掛けたものがエントロピー SS で、2つの系を触れさせると状態数の積を最大にするエネルギーの分け方がほぼ確実に実現し、∂S/∂E\partial S/\partial E がそろう。これを 1/T1/T と定め、体積についても同じ論法で P/T=∂S/∂VP/T=\partial S/\partial V を得た。理想気体では Φ∝VNE3N/2\Phi\propto V^NE^{3N/2} から E=32NkBTE=\frac32N k_{\mathrm B} T と PV=NkBTPV=N k_{\mathrm B} T が出る。高次元の体積は表面に集まるので、殻で数えても球で数えても、殻の厚みをどう選んでも、エントロピーは O(ln⁡N)\mathcal{O}(\ln N) しか変わらない。球の体積 Φ(E)\Phi(E) は自由度1では作用変数そのもので、断熱変化で変わらない。

次の章(第3部第3章)では、大きな系の一部分を取り出し、残りを熱浴とみなすことで、ミクロカノニカル分布からカノニカル分布 e−βH/Z\mathrm{e}^{-\beta H}/Z を導く。下の総合問題の (3) が、その最初の一歩になっている。

総合問題 独立な調和振動子の集まりを数える ★★★

オリジナル質量 mm、角振動数 ω\omega の1次元調和振動子が NN 個あり、格子点ごとに置かれていて互いに区別できる(dΓ=∏j=1Ndqj dpj/h\mathrm{d}\Gamma=\prod_{j=1}^{N}\mathrm{d} q_j\,\mathrm{d} p_j/h、N!N! で割らない)。振動子どうしはエネルギーをやり取りできる程度に弱く結合しているが、全エネルギーへの結合の寄与は無視する。全エネルギーを EE とする。

  1. (1) 1つの振動子を作用角変数 (θj,Jj)(\theta_j,J_j) で書くと Hj=ωJjH_j=\omega J_j である。これを使って Φ(E)=1N!(Eℏω)N\Phi(E)=\frac{1}{N!}\Bigl(\frac{E}{\hbar\omega}\Bigr)^N を示せ。また、ミクロカノニカル分布が時間で変わらないことを Poisson 括弧で示せ。(使う道具:第1部第14章の作用角変数・第3部第1章の Liouville 方程式・3.2.2節)
  2. (2) NN が大きいとして S=kBln⁡ΦS= k_{\mathrm B}\ln\Phi を Stirling の公式で整理し、温度とエネルギーの関係を求めよ。(使う道具:3.2.3節・3.2.4節)
  3. (3) 1つの振動子のエネルギー E1E_1 の分布 P(E1)P(E_1) を求め、平均が E/NE/N であることを示せ。N=10N=10 で E1>E/2E_1>E/2 となる確率を求めよ。また、E=NkBTE=N k_{\mathrm B} T を保ったまま N→∞N\to\infty とすると P(E1)P(E_1) が e−E1/kBT/(kBT)\mathrm{e}^{-E_1/ k_{\mathrm B} T}/( k_{\mathrm B} T) に近づくことを示せ。(使う道具:3.2.4節)
  4. (4) すべての振動子の角振動数を ω\omega から 2ω2\omega へ、振動の周期に比べてゆっくり変える。EE、TT、SS はそれぞれ何倍になるか。(使う道具:第1部第14章の断熱不変量・3.2.5節)
  5. (5) 量子力学では N=3N=3 の振動子に MM 個のエネルギー量子を分ける方法の数が Wq=(M+22)W_{\mathrm q}=\binom{M+2}{2} である。M=4M=4 のときの WqW_{\mathrm q} と、古典の Wcl=g(E) ℏωW_{\mathrm{cl}}=g(E)\,\hbar\omega(E=(M+32)ℏωE=(M+\frac32)\hbar\omega)を比べ、差が MM によらないことを示せ。(使う道具:3.2.6節)
解答

(1) 作用角変数への変換は正準変換なので体積を保ち、dqj dpj=dθj dJj\mathrm{d} q_j\,\mathrm{d} p_j=\mathrm{d}\theta_j\,\mathrm{d} J_j である。θj\theta_j の積分は 2π2\pi を与えるので、εj=ωJj\varepsilon_j=\omega J_j と置くと

Φ(E)=∏j2πhω∫∑jεj≤E, εj≥0∏jdεj=(1ℏω)N⋅ENN!.\Phi(E)=\prod_j\frac{2\pi}{h\omega}\int_{\sum_j\varepsilon_j\le E,\ \varepsilon_j\ge0}\prod_j\mathrm{d}\varepsilon_j=\Bigl(\frac{1}{\hbar\omega}\Bigr)^N\cdot\frac{E^N}{N!}.

最後の EN/N!E^N/N! は NN 次元の単体 {εj≥0, ∑εj≤E}\{\varepsilon_j\ge0,\ \sum\varepsilon_j\le E\} の体積である(N=1,2,3N=1,2,3 で EE, E2/2E^2/2, E3/6E^3/6)。ミクロカノニカル分布は H=∑jHjH=\sum_jH_j の関数 F(H)F(H) なので、{H,F(H)}=F′(H){H,H}=0\{H,F(H)\}=F'(H)\{H,H\}=0 で、Liouville 方程式 (3.2.8) から ∂ρ/∂t=0\partial\rho/\partial t=0 である。

(2) ln⁡Φ=Nln⁡Eℏω−ln⁡N!≈Nln⁡ENℏω+N\ln\Phi=N\ln\frac{E}{\hbar\omega}-\ln N!\approx N\ln\frac{E}{N\hbar\omega}+N。S=kB[Nln⁡ENℏω+N]S= k_{\mathrm B}[N\ln\frac{E}{N\hbar\omega}+N] を EE で微分して 1T=NkBE\frac1T=\frac{N k_{\mathrm B}}{E}、すなわち E=NkBTE=N k_{\mathrm B} T。1つの振動子に kBTk_{\mathrm B} T(運動エネルギーと位置エネルギーに 12kBT\frac12 k_{\mathrm B} T ずつ)が割り当てられる。

(3) 式 (3.2.15) で、系1を振動子1つ(g1=1/(ℏω)g_1=1/(\hbar\omega)、一定)、系2を残りの N−1N-1 個(g2(E2)=E2N−2/[(N−2)!(ℏω)N−1]g_2(E_2)=E_2^{N-2}/[(N-2)!(\hbar\omega)^{N-1}])とすると

P(E1)=g1 g2(E−E1)g(E)=(N−1)(E−E1)N−2EN−1(0≤E1≤E).P(E_1)=\frac{g_1\,g_2(E-E_1)}{g(E)}=\frac{(N-1)(E-E_1)^{N-2}}{E^{N-1}}\qquad(0\le E_1\le E).

∫0EP dE1=1\int_0^EP\,\mathrm{d} E_1=1 で、平均は ∫0EE1P dE1=E/N\int_0^EE_1P\,\mathrm{d} E_1=E/N である。N=10N=10 で E1>E/2E_1>E/2 となる確率は ∫E/2EP dE1=(1/2)N−1=1/512≈0.00195\int_{E/2}^EP\,\mathrm{d} E_1=(1/2)^{N-1}=1/512\approx0.00195。(同じ E/N=kBTE/N= k_{\mathrm B} T で次の章の Boltzmann 因子を使うと e−E1/kBT\mathrm{e}^{-E_1/ k_{\mathrm B} T} の裾は e−5=0.00674\mathrm{e}^{-5}=0.00674 で、N=10N=10 ではまだ違う。)E=NkBTE=N k_{\mathrm B} T と置いて N→∞N\to\infty とすると

P(E1)=N−1NkBT(1−E1NkBT)N−2⟶1kBTe−E1/kBTP(E_1)=\frac{N-1}{N k_{\mathrm B} T}\Bigl(1-\frac{E_1}{N k_{\mathrm B} T}\Bigr)^{N-2}\longrightarrow\frac{1}{ k_{\mathrm B} T}\mathrm{e}^{-E_1/ k_{\mathrm B} T}

である。全体がミクロカノニカル分布でも、大きな系の中の1つの振動子のエネルギーは Boltzmann 因子 e−E1/kBT\mathrm{e}^{-E_1/ k_{\mathrm B} T} に従って分布する。これがカノニカル分布(第3部第3章)の出どころである。

(4) 各振動子の作用変数 JjJ_j は断熱不変量なので、ゆっくり変える間 Φ\Phi は変わらない。弱い結合でエネルギーが振動子の間を行き来する場合にも全体の Φ(E)\Phi(E) が断熱不変量になること(3.2.5節で予告した P. Hertz の定理)を使う。その証明は第3部第5章で行う。Φ=1N!(E/ℏω)N\Phi=\frac{1}{N!}(E/\hbar\omega)^N が一定なので E/ωE/\omega が一定で、EE は 22 倍になる。E=NkBTE=N k_{\mathrm B} T から TT も 22 倍、S=kBln⁡ΦS= k_{\mathrm B}\ln\Phi は変わらない(11 倍)。これは1つの振動子で E=ωJE=\omega J が ω\omega に比例して増えること(第1部第14章)を NN 個に広げたものである。

(5) Wq=(62)=15W_{\mathrm q}=\binom{6}{2}=15。古典では g(E)=dΦdE=EN−1(N−1)!(ℏω)Ng(E)=\frac{\mathrm{d} \Phi}{\mathrm{d} E}=\frac{E^{N-1}}{(N-1)!(\hbar\omega)^N} なので Wcl=12(M+32)2=(5.5)22=15.125W_{\mathrm{cl}}=\frac{1}{2}\bigl(M+\frac32\bigr)^2=\frac{(5.5)^2}{2}=15.125。一般の MM では

Wcl−Wq=M2+3M+942−M2+3M+22=18W_{\mathrm{cl}}-W_{\mathrm q}=\frac{M^2+3M+\frac94}{2}-\frac{M^2+3M+2}{2}=\frac18

で、MM によらない。状態の数そのものは M2M^2 で増えるので、相対的な差は MM が大きいほど小さくなる。零点エネルギー 32ℏω\frac32\hbar\omega を EE に含めて数えたことで、差が O(M)\mathcal{O}(M) でなく定数に収まっている。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ