解析力学から様々な理論分野を組み立てる 第3部 第2章 §3.2.4 ★★★ 2つの系を触れさせると、状態数の積を最大にするエネルギーの分け方がほぼ確実に実現し、そこで ∂ S / ∂ E \partial S/\partial E ∂ S / ∂ E がそろう。これを 1 / T 1/T 1/ T と呼ぶ
2つの系を触れさせる
2つの系1と2を、エネルギーだけをやり取りできるように触れさせる。Hamiltonian は H = H 1 + H 2 + H i n t H=H_1+H_2+H_{\mathrm{int}} H = H 1 + H 2 + H int で、相互作用 H i n t H_{\mathrm{int}} H int はエネルギーの行き来を起こす程度には働くが、全エネルギーへの寄与は無視できるほど小さいとする(壁を通して熱だけを伝える状況)。全体を孤立系とみなしてミクロカノニカル分布を当てると、系1がエネルギー E 1 E_1 E 1 を持つ確率は、3.2.2 節の例題と同じ計算で
P ( E 1 ) = g 1 ( E 1 ) g 2 ( E − E 1 ) g ( E ) , g ( E ) = ∫ g 1 ( E 1 ) g 2 ( E − E 1 ) d E 1 (3.2.15) P(E_1)=\frac{g_1(E_1)\,g_2(E-E_1)}{g(E)},\qquad g(E)=\int g_1(E_1)\,g_2(E-E_1)\,\mathrm{d} E_1\tag{3.2.15} P ( E 1 ) = g ( E ) g 1 ( E 1 ) g 2 ( E − E 1 ) , g ( E ) = ∫ g 1 ( E 1 ) g 2 ( E − E 1 ) d E 1 ( 3.2.15 )
になる。全体の状態のうち「系1が E 1 E_1 E 1 、系2が E − E 1 E-E_1 E − E 1 」に当たるものの数は、それぞれの状態の数の積である。
N N N が大きい系では g 1 ( E 1 ) g_1(E_1) g 1 ( E 1 ) は E 1 E_1 E 1 の大きな冪で増え、g 2 ( E − E 1 ) g_2(E-E_1) g 2 ( E − E 1 ) は大きな冪で減るので、積は鋭いピークを持つ。ピークの位置は ln g 1 ( E 1 ) + ln g 2 ( E − E 1 ) \ln g_1(E_1)+\ln g_2(E-E_1) ln g 1 ( E 1 ) + ln g 2 ( E − E 1 ) の最大で決まり、
∂ ln g 1 ∂ E 1 = ∂ ln g 2 ∂ E 2 ( E 2 = E − E 1 ) (3.2.16) \frac{\partial \ln g_1}{\partial E_1}=\frac{\partial \ln g_2}{\partial E_2}\qquad(E_2=E-E_1)\tag{3.2.16} ∂ E 1 ∂ ln g 1 = ∂ E 2 ∂ ln g 2 ( E 2 = E − E 1 ) ( 3.2.16 )
を満たす。式 (3.2.6) のエントロピー S j = k B ln W j S_j= k_{\mathrm B}\ln W_j S j = k B ln W j で書けば ∂ S 1 / ∂ E 1 = ∂ S 2 / ∂ E 2 \partial S_1/\partial E_1=\partial S_2/\partial E_2 ∂ S 1 / ∂ E 1 = ∂ S 2 / ∂ E 2 である。触れさせた2つの系で等しくなる量は、熱力学では温度だった。そこで
1 T = ∂ S ∂ E (3.2.17) \frac{1}{T}=\frac{\partial S}{\partial E}\tag{3.2.17} T 1 = ∂ E ∂ S ( 3.2.17 )
を温度の定義にする。S S S が E E E の増加関数なら T > 0 T>0 T > 0 で、エネルギーは ∂ S / ∂ E \partial S/\partial E ∂ S / ∂ E の小さい系(温度の高い系)から大きい系(温度の低い系)へ移る。エネルギーが移ると S 1 + S 2 S_1+S_2 S 1 + S 2 が増え、ピークに着いたところで止まる。熱力学第二法則のこの形が、状態の数え上げから出てくる。
理想気体では g j ∝ E j 3 N j / 2 − 1 g_j\propto E_j^{3N_j/2-1} g j ∝ E j 3 N j /2 − 1 なので、式 (3.2.15) は x = E 1 / E x=E_1/E x = E 1 / E について
P ( x ) ∝ x 3 N 1 / 2 − 1 ( 1 − x ) 3 N 2 / 2 − 1 (3.2.18) P(x)\propto x^{3N_1/2-1}(1-x)^{3N_2/2-1}\tag{3.2.18} P ( x ) ∝ x 3 N 1 /2 − 1 ( 1 − x ) 3 N 2 /2 − 1 ( 3.2.18 )
となる(Beta 分布)。平均は N 1 / ( N 1 + N 2 ) N_1/(N_1+N_2) N 1 / ( N 1 + N 2 ) で、エネルギーは粒子の数に比例して分けられる。N 1 = N 2 = N N_1=N_2=N N 1 = N 2 = N のとき標準偏差は 1 / ( 2 3 N + 1 ) 1/(2\sqrt{3N+1}) 1/ ( 2 3 N + 1 ) で、N N N とともに N − 1 / 2 N^{-1/2} N − 1/2 で狭くなる(図 3.2.3 )。
図 3.2.3 同じ粒子数 N N N の2つの理想気体を触れさせたときの、系1の受け持つエネルギーの割合 x = E 1 / E x=E_1/E x = E 1 / E の分布(式 (3.2.18) )。点線(厳密)は N = 1 , 10 , 100 N=1,10,100 N = 1 , 10 , 100 の Beta 分布、細い実線(近似)は N = 100 N=100 N = 100 の分布と同じ平均 1 2 \frac12 2 1 ・標準偏差 1 / ( 2 301 ) 1/(2\sqrt{301}) 1/ ( 2 301 ) の正規分布で、点線とほとんど重なる。ピークは N N N とともに N − 1 / 2 N^{-1/2} N − 1/2 で狭くなる。
例題 3.2.5 エネルギーの分け方の鋭さ ★★
オリジナル 2つの理想気体(粒子数 N 1 N_1 N 1 , N 2 N_2 N 2 )を触れさせる。式 (3.2.18) を使って次に答えよ。
(1) N 1 = 100 N_1=100 N 1 = 100 , N 2 = 300 N_2=300 N 2 = 300 のとき、x = E 1 / E x=E_1/E x = E 1 / E の最も確からしい値(最頻値)と平均を求めよ。
(2) N 1 = N 2 = 100 N_1=N_2=100 N 1 = N 2 = 100 のときの x x x の標準偏差を求めよ。
(3) N 1 = N 2 = 6.02 × 10 23 N_1=N_2=6.02\times10^{23} N 1 = N 2 = 6.02 × 1 0 23 のときの x x x の標準偏差を求め、温度計で測れる程度の揺らぎかどうかを述べよ。
解答 (1) ln P = ( a − 1 ) ln x + ( b − 1 ) ln ( 1 − x ) \ln P=(a-1)\ln x+(b-1)\ln(1-x) ln P = ( a − 1 ) ln x + ( b − 1 ) ln ( 1 − x ) (a = 3 N 1 / 2 = 150 a=3N_1/2=150 a = 3 N 1 /2 = 150 , b = 3 N 2 / 2 = 450 b=3N_2/2=450 b = 3 N 2 /2 = 450 )を x x x で微分して 0 0 0 と置くと x ∗ = ( a − 1 ) / ( a + b − 2 ) = 149 / 598 ≈ 0.2492 x^*=(a-1)/(a+b-2)=149/598\approx0.2492 x ∗ = ( a − 1 ) / ( a + b − 2 ) = 149/598 ≈ 0.2492 。平均は Beta 分布の公式 a / ( a + b ) = 0.25 a/(a+b)=0.25 a / ( a + b ) = 0.25 で、N 1 / ( N 1 + N 2 ) N_1/(N_1+N_2) N 1 / ( N 1 + N 2 ) に等しい。最頻値と平均の差は O ( 1 / N ) \mathcal{O}(1/N) O ( 1/ N ) で、N N N が大きいと区別がなくなる。
(2) Beta 分布の分散は a b / [ ( a + b ) 2 ( a + b + 1 ) ] ab/[(a+b)^2(a+b+1)] ab / [( a + b ) 2 ( a + b + 1 )] で、a = b = 150 a=b=150 a = b = 150 なら 1 4 ⋅ 301 \frac{1}{4\cdot301} 4 ⋅ 301 1 。標準偏差は 1 / ( 2 301 ) ≈ 0.0288 1/(2\sqrt{301})\approx0.0288 1/ ( 2 301 ) ≈ 0.0288 。
(3) a = b = 3 2 N a=b=\frac32N a = b = 2 3 N なら標準偏差は 1 / ( 2 3 N + 1 ) ≈ 3.7 × 10 − 13 1/(2\sqrt{3N+1})\approx3.7\times10^{-13} 1/ ( 2 3 N + 1 ) ≈ 3.7 × 1 0 − 13 である。エネルギーの割合が 10 − 13 10^{-13} 1 0 − 13 の桁でしか揺れないので、巨視的な系ではエネルギーの分け方はピークの位置にほぼ確実に決まり、温度計で測れるような揺らぎは見えない。
理想気体の温度と圧力
式 (3.2.14) から、N N N が大きい理想気体では S = k B ln Φ S= k_{\mathrm B}\ln\Phi S = k B ln Φ と置いてよい(殻の W W W と体積の Φ \Phi Φ の違いは3.2.5 節)。E E E についての部分は 3 2 N k B ln E \frac32N k_{\mathrm B}\ln E 2 3 N k B ln E なので、
1 T = ∂ S ∂ E = 3 N k B 2 E ⟹ E = 3 2 N k B T (3.2.19) \frac1T=\frac{\partial S}{\partial E}=\frac{3N k_{\mathrm B}}{2E}\qquad\Longrightarrow\qquad E=\frac32N k_{\mathrm B} T\tag{3.2.19} T 1 = ∂ E ∂ S = 2 E 3 N k B ⟹ E = 2 3 N k B T ( 3.2.19 )
である。1粒子あたり、運動の3方向それぞれに 1 2 k B T \frac12 k_{\mathrm B} T 2 1 k B T のエネルギーがある(等分配則。一般の形は第3部第6章 )。
体積も同じ方法で扱える。2つの系の間の壁が動ける(ピストン)なら、体積の分け方 V 1 + V 2 = V V_1+V_2=V V 1 + V 2 = V についても状態数の積を最大にする分け方が実現し、∂ S 1 / ∂ V 1 = ∂ S 2 / ∂ V 2 \partial S_1/\partial V_1=\partial S_2/\partial V_2 ∂ S 1 / ∂ V 1 = ∂ S 2 / ∂ V 2 がそろう。力学的につり合ったとき2つの系で等しいのは圧力なので、P / T = ∂ S / ∂ V P/T=\partial S/\partial V P / T = ∂ S / ∂ V と置く。理想気体では S S S の V V V についての部分が N k B ln V N k_{\mathrm B}\ln V N k B ln V なので
P T = ∂ S ∂ V = N k B V ⟹ P V = N k B T (3.2.20) \frac{P}{T}=\frac{\partial S}{\partial V}=\frac{N k_{\mathrm B}}{V}\qquad\Longrightarrow\qquad PV=N k_{\mathrm B} T\tag{3.2.20} T P = ∂ V ∂ S = V N k B ⟹ P V = N k B T ( 3.2.20 )
を得る。理想気体の状態方程式は、Φ \Phi Φ が V N V^N V N に比例すること(各粒子が容器のどこにいてもよいこと)の言い換えである。S ( E , V , N ) S(E,V,N) S ( E , V , N ) を E E E と V V V で微分した量が温度と圧力になる、という組み立ては、熱力学の基本関係式 d E = T d S − P d V \mathrm{d} E=T\,\mathrm{d} S-P\,\mathrm{d} V d E = T d S − P d V をそのまま S S S について解いた形 d S = d E / T + ( P / T ) d V \mathrm{d} S=\mathrm{d} E/T+(P/T)\,\mathrm{d} V d S = d E / T + ( P / T ) d V である。自然な変数を E E E から T T T に替える操作は Legendre 変換で、第3部第4章 で扱う。
例題 3.2.6 アルゴンのエントロピーを数える ★★
オリジナル アルゴン(原子の質量 39.948 u 39.948\,\mathrm u 39.948 u 、1 u = 1.66054 × 10 − 27 k g 1\,\mathrm u=1.66054\times10^{-27}\,\mathrm{kg} 1 u = 1.66054 × 1 0 − 27 kg )の気体が、温度 300 K 300\,\mathrm K 300 K 、圧力 1.0 × 10 5 P a 1.0\times10^{5}\,\mathrm{Pa} 1.0 × 1 0 5 Pa にある。理想気体とみなす。
(1) 熱的波長 λ T = h / 2 π m k B T \lambda_T=h/\sqrt{2\pi m k_{\mathrm B} T} λ T = h / 2 π m k B T と、1原子あたりの体積 V / N = k B T / P V/N= k_{\mathrm B} T/P V / N = k B T / P を求めよ。
(2) 例題 3.2.4 の式で1原子あたりのエントロピー S / ( N k B ) S/(N k_{\mathrm B}) S / ( N k B ) と、1 mol あたりのエントロピーを求めよ。
(3) V / ( N λ T 3 ) V/(N\lambda_T^3) V / ( N λ T 3 ) の大きさから、相空間の体積で古典的に数えてよいかを判断せよ。
解答 (1) m = 39.948 × 1.66054 × 10 − 27 = 6.633 × 10 − 26 k g m=39.948\times1.66054\times10^{-27}=6.633\times10^{-26}\,\mathrm{kg} m = 39.948 × 1.66054 × 1 0 − 27 = 6.633 × 1 0 − 26 kg で、λ T = 6.626 × 10 − 34 / 2 π ⋅ 6.633 × 10 − 26 ⋅ 1.381 × 10 − 23 ⋅ 300 = 1.59 × 10 − 11 m \lambda_T=6.626\times10^{-34}/\sqrt{2\pi\cdot6.633\times10^{-26}\cdot1.381\times10^{-23}\cdot300}=1.59\times10^{-11}\,\mathrm m λ T = 6.626 × 1 0 − 34 / 2 π ⋅ 6.633 × 1 0 − 26 ⋅ 1.381 × 1 0 − 23 ⋅ 300 = 1.59 × 1 0 − 11 m 。V / N = 1.381 × 10 − 23 ⋅ 300 / ( 1.0 × 10 5 ) = 4.14 × 10 − 26 m 3 V/N=1.381\times10^{-23}\cdot300/(1.0\times10^5)=4.14\times10^{-26}\,\mathrm{m^3} V / N = 1.381 × 1 0 − 23 ⋅ 300/ ( 1.0 × 1 0 5 ) = 4.14 × 1 0 − 26 m 3 。
(2) V / ( N λ T 3 ) = 4.14 × 10 − 26 / ( 1.59 × 10 − 11 ) 3 = 1.02 × 10 7 V/(N\lambda_T^3)=4.14\times10^{-26}/(1.59\times10^{-11})^3=1.02\times10^{7} V / ( N λ T 3 ) = 4.14 × 1 0 − 26 / ( 1.59 × 1 0 − 11 ) 3 = 1.02 × 1 0 7 なので、S / ( N k B ) = ln ( 1.02 × 10 7 ) + 5 2 = 16.14 + 2.50 = 18.6 S/(N k_{\mathrm B})=\ln(1.02\times10^7)+\frac52=16.14+2.50=18.6 S / ( N k B ) = ln ( 1.02 × 1 0 7 ) + 2 5 = 16.14 + 2.50 = 18.6 。1 mol あたりは N A k B = 8.314 J m o l − 1 K − 1 N_{\mathrm A} k_{\mathrm B}=8.314\,\mathrm{J\,mol^{-1}K^{-1}} N A k B = 8.314 J mo l − 1 K − 1 を掛けて 155 J m o l − 1 K − 1 155\,\mathrm{J\,mol^{-1}K^{-1}} 155 J mo l − 1 K − 1 である。
(3) λ T 3 \lambda_T^3 λ T 3 は、1つの原子の位置と運動量の揺らぎが相空間で占める大きさ h 3 h^3 h 3 を、位置の体積に換算したものである。V / ( N λ T 3 ) ≈ 10 7 V/(N\lambda_T^3)\approx10^7 V / ( N λ T 3 ) ≈ 1 0 7 は、原子1つに 10 7 10^7 1 0 7 個ほどの状態(体積 h 3 h^3 h 3 のます目)が割り当てられていることを意味する。2つの原子が同じます目に入ることはまずないので、量子力学の粒子の統計(Bose・Fermi)の違いは表に出ず、N ! N! N ! で割る古典的な数え方でよい。V / ( N λ T 3 ) V/(N\lambda_T^3) V / ( N λ T 3 ) が 1 1 1 に近づく低温・高密度では、この数え方は破れる。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com