解析力学から様々な理論分野を組み立てる 第3部 第4章 §3.4.1 ★★★ 状態数から得た S ( E , V , N ) S(E,V,N) S ( E , V , N ) を E E E について解いた U ( S , V , N ) U(S,V,N) U ( S , V , N ) が出発点。1次同次性と凸性の2つの性質を持つ
熱力学の教科書には、内部エネルギー U U U ・Helmholtz の自由エネルギー F F F ・エンタルピー H e H_{\mathrm e} H e ・Gibbs エネルギー G G G ・グランドポテンシャル Ω \Omega Ω と、エネルギーの次元を持つ関数がいくつも並ぶ。それぞれに「自然な変数」があり、全微分の形と Maxwell 関係式を覚えることになっている。第1部第7章 では、U ( S , V ) U(S,V) U ( S , V ) のエントロピー S S S を温度 T = ∂ U / ∂ S T=\partial U/\partial S T = ∂ U / ∂ S に取り替える Legendre 変換として F = U − T S F=U-TS F = U − T S を作った。この章では、同じ操作を U ( S , V , N ) U(S,V,N) U ( S , V , N ) の3つの変数のどれを取り替えるかの組み合わせとして整理し、熱力学ポテンシャルがすべてその組み合わせの1つずつであることを示す。
この章で選ぶものは、ミクロな Hamiltonian H H H から状態数を数えて得られる1つの関数、基本関係式 U ( S , V , N ) U(S,V,N) U ( S , V , N ) である。当てる道具は、第1部第7章 の Legendre 変換を多変数に広げたものである。出てくるものは3つに分けて書く。3.4.3 節では4つのポテンシャルとその最小原理、Maxwell 関係式を出す。3.4.4 節では、統計力学の分配関数が状態数の Laplace 変換であり、粒子数が大きいと鞍点の値が Legendre 変換に一致することを示す。3.4.5 節では、凸でない自由エネルギーを Legendre 変換すると凸包に置き換わることが、気体と液体の共存(1次相転移)の記述になることを示す。3.4.6 節で授業の「全微分を取って項を移す」書き方と並べ、3.4.7 節で常磁性体を題材にした総合問題を解く。
ミクロな H H H から U ( S , V , N ) U(S,V,N) U ( S , V , N ) へ
第3部第2章 では、Hamiltonian H ( q , p ) H(q,p) H ( q , p ) を持つ N N N 粒子の系を体積 V V V の容器に入れ、エネルギーが幅 Δ E \Delta E Δ E の殻 [ E , E + Δ E ) [E,E+\Delta E) [ E , E + Δ E ) に入る状態の数 W ( E ) W(E) W ( E ) を数えて、エントロピーを S ( E , V , N ) = k B ln W ( E ) S(E,V,N)= k_{\mathrm B}\ln W(E) S ( E , V , N ) = k B ln W ( E ) と定めた。W ( E ) W(E) W ( E ) はふつう E E E について単調に増えるので、S = S ( E , V , N ) S=S(E,V,N) S = S ( E , V , N ) を E E E について解くことができる。その解を内部エネルギー U ( S , V , N ) U(S,V,N) U ( S , V , N ) と書き、これを基本関係式 と呼ぶ。この章で選ぶものはこの1つの関数で、熱力学のすべての量はこの関数の偏微分から出てくる。
U U U の1階偏微分に名前を付ける。
T = ∂ U ∂ S ∣ V , N , P = − ∂ U ∂ V ∣ S , N , μ = ∂ U ∂ N ∣ S , V , d U = T d S − P d V + μ d N . (3.4.1) T=\frac{\partial U}{\partial S}\Big|_{V,N},\qquad P=-\frac{\partial U}{\partial V}\Big|_{S,N},\qquad \mu=\frac{\partial U}{\partial N}\Big|_{S,V},\qquad
\mathrm{d} U=T\,\mathrm{d} S-P\,\mathrm{d} V+\mu\,\mathrm{d} N .\tag{3.4.1} T = ∂ S ∂ U V , N , P = − ∂ V ∂ U S , N , μ = ∂ N ∂ U S , V , d U = T d S − P d V + μ d N . ( 3.4.1 )
P P P にだけ負号が付くのは、体積を増やすと系が外へ仕事をしてエネルギーが減る向きを正の圧力と呼ぶ習慣による。( S , V , N ) (S,V,N) ( S , V , N ) は系を2つ並べると2倍になる量で示量変数 と呼び、それぞれの傾き ( T , − P , μ ) (T,-P,\mu) ( T , − P , μ ) は2つ並べても変わらない量で示強変数 と呼ぶ。示量変数と、それに対応する傾きの組を共役な組 と呼ぶ。
1次同次性と Euler の関係
短距離の相互作用しかない系では、同じ系を2つ並べてつなげても、境界での相互作用のエネルギーは全体に比べて小さい。したがって S , V , N S,V,N S , V , N をすべて λ \lambda λ 倍すると U U U も λ \lambda λ 倍になる。
U ( λ S , λ V , λ N ) = λ U ( S , V , N ) ( λ > 0 ) . (3.4.2) U(\lambda S,\lambda V,\lambda N)=\lambda\,U(S,V,N)\qquad(\lambda>0).\tag{3.4.2} U ( λ S , λV , λ N ) = λ U ( S , V , N ) ( λ > 0 ) . ( 3.4.2 )
この性質を1次同次性 と呼ぶ。両辺を λ \lambda λ で微分して λ = 1 \lambda=1 λ = 1 と置くと、左辺は S ∂ U / ∂ S + V ∂ U / ∂ V + N ∂ U / ∂ N S\,\partial U/\partial S+V\,\partial U/\partial V+N\,\partial U/\partial N S ∂ U / ∂ S + V ∂ U / ∂ V + N ∂ U / ∂ N になるので
U = T S − P V + μ N (3.4.3) U=TS-PV+\mu N\tag{3.4.3} U = T S − P V + μ N ( 3.4.3 )
を得る。これを Euler の関係 と呼ぶ。式 (3.4.3) の全微分を取り、式 (3.4.1) の d U \mathrm{d} U d U を引くと
S d T − V d P + N d μ = 0 (3.4.4) S\,\mathrm{d} T-V\,\mathrm{d} P+N\,\mathrm{d}\mu=0\tag{3.4.4} S d T − V d P + N d μ = 0 ( 3.4.4 )
が残る。これを Gibbs–Duhem の関係 と呼ぶ。3つの示強変数 ( T , P , μ ) (T,P,\mu) ( T , P , μ ) は独立に動けず、2つを決めると残りの1つが決まる。この関係は3.4.2 節で「3つの変数を全部取り替えると何も残らない」という形で戻ってくる。
凸性
U U U は ( S , V , N ) (S,V,N) ( S , V , N ) について下に凸である。理由は次のとおりである。系を壁で2つに仕切り、左に ( U 1 , V 1 , N 1 ) (U_1,V_1,N_1) ( U 1 , V 1 , N 1 ) 、右に ( U 2 , V 2 , N 2 ) (U_2,V_2,N_2) ( U 2 , V 2 , N 2 ) を入れる。壁を取り去ると、取り去る前に数えた状態はすべて取り去った後にも許されるので、状態数は減らない。したがってS ( U 1 + U 2 , V 1 + V 2 , N 1 + N 2 ) ≥ S ( U 1 , V 1 , N 1 ) + S ( U 2 , V 2 , N 2 ) S(U_1+U_2,V_1+V_2,N_1+N_2)\ge S(U_1,V_1,N_1)+S(U_2,V_2,N_2) S ( U 1 + U 2 , V 1 + V 2 , N 1 + N 2 ) ≥ S ( U 1 , V 1 , N 1 ) + S ( U 2 , V 2 , N 2 ) である。右辺を1次同次性で書き直すと、2点の中点で S S S が2点の値の平均以上になること、すなわち S S S が上に凸であることになる。S S S が U U U について増加関数で上に凸なら、その逆関数 U U U は S S S について下に凸である(V , N V,N V , N についても同じ議論が通る)。凸性は3.4.2 節で Legendre 変換が情報を落とさない条件として使い、3.4.5 節では、それが破れて見えるときに何が起こるかを調べる。
1次同次な関数は、原点を通る直線の上で1次式になるので、その方向には曲がっていない。U U U の ( S , V , N ) (S,V,N) ( S , V , N ) についての Hesse 行列は、( S , V , N ) (S,V,N) ( S , V , N ) 自身の方向に固有値 0 0 0 を持ち、残りの2方向で正の固有値を持つ。狭義に凸なのは、たとえば N N N を止めた ( S , V ) (S,V) ( S , V ) の2変数についてである。
例題 3.4.1 単原子理想気体の基本関係式 ★★★
オリジナル 単原子理想気体(粒子の質量 m m m )のエントロピーは、Planck 定数 h h h を含めてS = N k B [ ln { V N ( 4 π m U 3 N h 2 ) 3 / 2 } + 5 2 ] S=N k_{\mathrm B}\Bigl[\ln\Bigl\{\dfrac VN\Bigl(\dfrac{4\pi mU}{3Nh^2}\Bigr)^{3/2}\Bigr\}+\dfrac52\Bigr] S = N k B [ ln { N V ( 3 N h 2 4 π m U ) 3/2 } + 2 5 ] と書くことができる(Sackur–Tetrode の式。導出は第3部第2章 )。
(1) これを U U U について解き、U ( λ S , λ V , λ N ) = λ U U(\lambda S,\lambda V,\lambda N)=\lambda U U ( λ S , λV , λ N ) = λ U を確かめよ。
(2) T , P , μ T,P,\mu T , P , μ を求め、Euler の関係 (3.4.3) を確かめよ。
(3) μ \mu μ を T T T と P P P で表し、Gibbs–Duhem の関係 (3.4.4) を確かめよ。
解答 (1) S N k B − 5 2 = ln V N + 3 2 ln 4 π m U 3 N h 2 \dfrac{S}{N k_{\mathrm B}}-\dfrac52=\ln\dfrac VN+\dfrac32\ln\dfrac{4\pi mU}{3Nh^2} N k B S − 2 5 = ln N V + 2 3 ln 3 N h 2 4 π m U を U U U について解くと
U ( S , V , N ) = 3 h 2 N 5 / 3 4 π m V 2 / 3 exp ( 2 S 3 N k B − 5 3 ) . U(S,V,N)=\frac{3h^2N^{5/3}}{4\pi mV^{2/3}}\exp\Bigl(\frac{2S}{3N k_{\mathrm B}}-\frac53\Bigr). U ( S , V , N ) = 4 π m V 2/3 3 h 2 N 5/3 exp ( 3 N k B 2 S − 3 5 ) .
S , V , N S,V,N S , V , N を λ \lambda λ 倍すると、指数の中の S / N S/N S / N は変わらず、前の因子 N 5 / 3 V − 2 / 3 N^{5/3}V^{-2/3} N 5/3 V − 2/3 は λ 5 / 3 − 2 / 3 = λ \lambda^{5/3-2/3}=\lambda λ 5/3 − 2/3 = λ 倍になる。よって U U U は λ \lambda λ 倍になる。
(2) S S S で微分すると指数の係数 2 / ( 3 N k B ) 2/(3N k_{\mathrm B}) 2/ ( 3 N k B ) が出るので T = 2 U / ( 3 N k B ) T=2U/(3N k_{\mathrm B}) T = 2 U / ( 3 N k B ) 、すなわち U = 3 2 N k B T U=\frac32N k_{\mathrm B} T U = 2 3 N k B T 。V V V で微分すると V − 2 / 3 V^{-2/3} V − 2/3 から − 2 3 V -\frac{2}{3V} − 3 V 2 が出るので P = 2 U / ( 3 V ) P=2U/(3V) P = 2 U / ( 3 V ) 。2つを合わせると P V = N k B T PV=N k_{\mathrm B} T P V = N k B T 。N N N で微分すると、前の因子から 5 3 N \frac{5}{3N} 3 N 5 、指数から − 2 S 3 N 2 k B -\frac{2S}{3N^2 k_{\mathrm B}} − 3 N 2 k B 2 S が出るので μ = U ( 5 3 N − 2 S 3 N 2 k B ) \mu=U\bigl(\frac{5}{3N}-\frac{2S}{3N^2 k_{\mathrm B}}\bigr) μ = U ( 3 N 5 − 3 N 2 k B 2 S ) 。これらを並べると
T S − P V + μ N = 2 U S 3 N k B − 2 U 3 + 5 U 3 − 2 U S 3 N k B = U TS-PV+\mu N=\frac{2US}{3N k_{\mathrm B}}-\frac{2U}{3}+\frac{5U}{3}-\frac{2US}{3N k_{\mathrm B}}=U T S − P V + μ N = 3 N k B 2 U S − 3 2 U + 3 5 U − 3 N k B 2 U S = U
で、Euler の関係が成り立つ。
(3) (2) の μ \mu μ に U = 3 2 N k B T U=\frac32N k_{\mathrm B} T U = 2 3 N k B T と Sackur–Tetrode の S S S を入れて整理すると、熱的 de Broglie 波長 λ 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 を使って
μ = k B T ln ( N V λ T 3 ) = k B T ln P λ T 3 k B T \mu= k_{\mathrm B} T\ln\Bigl(\frac NV\lambda_T^3\Bigr)= k_{\mathrm B} T\ln\frac{P\lambda_T^3}{ k_{\mathrm B} T} μ = k B T ln ( V N λ T 3 ) = k B T ln k B T P λ T 3
となる。λ T 3 ∝ T − 3 / 2 \lambda_T^3\propto T^{-3/2} λ T 3 ∝ T − 3/2 なので μ = k B T [ ln P − 5 2 ln T + 定数 ] \mu= k_{\mathrm B} T\bigl[\ln P-\frac52\ln T+\text{定数}\bigr] μ = k B T [ ln P − 2 5 ln T + 定数 ] で、∂ μ ∂ P ∣ T = k B T / P = V / N \frac{\partial \mu}{\partial P}\big|_T= k_{\mathrm B} T/P=V/N ∂ P ∂ μ T = k B T / P = V / N 、∂ μ ∂ T ∣ P = μ / T − 5 2 k B = − S / N \frac{\partial \mu}{\partial T}\big|_P=\mu/T-\frac52 k_{\mathrm B}=-S/N ∂ T ∂ μ P = μ / T − 2 5 k B = − S / N (最後の等号は Sackur–Tetrode の式を μ \mu μ で書き直したもの)。したがって d μ = − ( S / N ) d T + ( V / N ) d P \mathrm{d}\mu=-(S/N)\mathrm{d} T+(V/N)\mathrm{d} P d μ = − ( S / N ) d T + ( V / N ) d P で、両辺に N N N を掛けると式 (3.4.4) になる。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com