解析力学から様々な理論分野を組み立てる 第1部 第7章 §1.7.5 ★★★ F ( T , V ) = inf S [ U − T S ] = − U ∗ ( T ) F(T,V)=\inf_S[U-TS]=-U^*(T) F ( T , V ) = inf S [ U − T S ] = − U ∗ ( T ) 。体積は取り替えない変数で、L → H L\to H L → H と同じ図に載る
同じ図に載せる
熱力学では、粒子数 N N N を止めた系の内部エネルギーを U ( S , V ) U(S,V) U ( S , V ) と書き、温度と圧力を
T = ∂ U ∂ S ∣ V , P = − ∂ U ∂ V ∣ S T=\frac{\partial U}{\partial S}\Big|_V,\qquad P=-\frac{\partial U}{\partial V}\Big|_S T = ∂ S ∂ U V , P = − ∂ V ∂ U S
で定める。安定な系では T > 0 T>0 T > 0 で、U U U は S S S について下に凸、∂ 2 U / ∂ S 2 ∣ V = T / C V > 0 \partial^2U/\partial S^2|_V=T/C_V>0 ∂ 2 U / ∂ S 2 ∣ V = T / C V > 0 (C V C_V C V は定積熱容量)である。S S S は実験で直接制御できないが、T T T は熱浴に接触させて一定に保てるので、S S S を T T T に取り替えたい。U U U の S S S についての傾きが T T T なので、これは1.7.2 節の Legendre 変換で x → S x\to S x → S 、p → T p\to T p → T とし、V V V を取り替えない変数としたものである。
熱力学の習慣では、切片に負号を付けずにそのまま使う。傾き T T T の接線の切片は U − T S U-TS U − T S で、これを Helmholtz の自由エネルギー と呼ぶ。
F ( T , V ) ≔ inf S [ U ( S , V ) − T S ] = − U ∗ ( T ) ( U ∗ は S についての Legendre 変換、 V は止める ) . (1.7.13) F(T,V)\coloneqq\inf_S\bigl[U(S,V)-TS\bigr]=-U^*(T)\qquad(U^*\text{ は }S\text{ についての Legendre 変換、}V\text{ は止める}).\tag{1.7.13} F ( T , V ) : = S inf [ U ( S , V ) − T S ] = − U ∗ ( T ) ( U ∗ は S についての Legendre 変換、 V は止める ) . ( 1.7.13 )
inf \inf inf と sup \sup sup は負号で入れ替わる(inf S [ U − T S ] = − sup S [ T S − U ] \inf_S[U-TS]=-\sup_S[TS-U] inf S [ U − T S ] = − sup S [ T S − U ] )。F F F は T T T について上に凸(− U ∗ -U^* − U ∗ だから)になる。
図 1.7.4 は、L → H L\to H L → H と U → F U\to F U → F を同じ図に描いたものである。左は L L L を q ˙ \dot q q ˙ の関数として描き、傾き p p p の接線の切片が − H -H − H になる。右は U U U を S S S の関数として描き、傾き T T T の接線の切片が F F F になる。2つは同じ作図で、違いは切片に負号を付けるかどうかの約束だけである。H H H は L L L の Legendre 変換そのもの、F F F は U U U の Legendre 変換の − 1 -1 − 1 倍である。
図 1.7.4 同じ作図の2つの例。左:相対論的な粒子の L = − m c 2 1 − q ˙ 2 / c 2 L=-mc^2\sqrt{1-\dot q^2/c^2} L = − m c 2 1 − q ˙ 2 / c 2 (V = 0 V=0 V = 0 )。q ˙ = 0.6 c \dot q=0.6c q ˙ = 0.6 c での傾きは p = 0.75 m c p=0.75mc p = 0.75 m c 、切片は − 1.25 m c 2 = − H -1.25mc^2=-H − 1.25 m c 2 = − H 。右:体積を止めた単原子理想気体の U ∝ e 2 S / ( 3 N k B ) U\propto\mathrm{e}^{2S/(3N k_{\mathrm B})} U ∝ e 2 S / ( 3 N k B ) (単位を 3 N k B / 2 = 1 3N k_{\mathrm B}/2=1 3 N k B /2 = 1 に選んだ)。S = 0.5 S=0.5 S = 0.5 での傾き T T T の接線の切片が F = U − T S F=U-TS F = U − T S 。曲線は実線、接線は破線(どちらも近似を含まない)
左の図の L L L は ∣ q ˙ ∣ < c \lvert \dot q\rvert <c ∣ q ˙ ∣ < c でしか定義されないが、傾き p = m q ˙ / 1 − q ˙ 2 / c 2 p=m\dot q/\sqrt{1-\dot q^2/c^2} p = m q ˙ / 1 − q ˙ 2 / c 2 は q ˙ → ± c \dot q\to\pm c q ˙ → ± c で ± ∞ \pm\infty ± ∞ に向かうので、どんな p p p にも接線を引くことができる。その結果 H = p 2 c 2 + m 2 c 4 H=\sqrt{p^2c^2+m^2c^4} H = p 2 c 2 + m 2 c 4 はすべての p p p で定義され、p → ± ∞ p\to\pm\infty p → ± ∞ で傾きが ± c \pm c ± c に近づく。L L L の定義域の端が、H H H の傾きの上限として現れている(相対論的な粒子の L L L の由来は第2部第3章 )。
取り替えない変数と2階微分
規則 (1.7.10) と同じ計算を F = U − T S F=U-TS F = U − T S (最小点で ∂ U / ∂ S = T \partial U/\partial S=T ∂ U / ∂ S = T )に当てると
∂ F ∂ T ∣ V = − S , ∂ F ∂ V ∣ T = ∂ U ∂ V ∣ S = − P (1.7.14) \frac{\partial F}{\partial T}\Big|_V=-S,\qquad \frac{\partial F}{\partial V}\Big|_T=\frac{\partial U}{\partial V}\Big|_S=-P\tag{1.7.14} ∂ T ∂ F V = − S , ∂ V ∂ F T = ∂ V ∂ U S = − P ( 1.7.14 )
を得る。1つ目は定理 1.7.2 の ( f ∗ ) ′ = x (f^*)'=x ( f ∗ ) ′ = x に負号が付いたもの、2つ目は取り替えない変数 V V V の偏微分が引き継がれたものである。F = − U ∗ F=-U^* F = − U ∗ と約束したので、H H H の場合と違って2つ目には負号が付かない。2階微分は定理 1.7.2 から
∂ 2 F ∂ T 2 ∣ V = − 1 ∂ 2 U / ∂ S 2 ∣ V = − C V T < 0 \frac{\partial^2F}{\partial T^2}\Big|_V=-\frac{1}{\partial^2U/\partial S^2|_V}=-\frac{C_V}{T}<0 ∂ T 2 ∂ 2 F V = − ∂ 2 U / ∂ S 2 ∣ V 1 = − T C V < 0
で、U U U が S S S について凸であることと、F F F が T T T について上に凸であることは同じ条件 C V > 0 C_V>0 C V > 0 の言い換えになる。
Young の不等式 (1.7.6) は、熱力学では U ( S , V ) − F ( T , V ) ≥ T S U(S,V)-F(T,V)\ge TS U ( S , V ) − F ( T , V ) ≥ T S 、すなわち任意の S S S について U ( S , V ) − T S ≥ F ( T , V ) U(S,V)-TS\ge F(T,V) U ( S , V ) − T S ≥ F ( T , V ) と読むことができる。温度 T T T の熱浴に接した系では、エントロピーをどんな値にしても U − T S U-TS U − T S は F F F を下回らず、平衡(∂ U / ∂ S = T \partial U/\partial S=T ∂ U / ∂ S = T )で等号になる。自由エネルギー最小の原理は、この不等式の等号条件である。
未定乗数との関係
U ( S , V ) − T S U(S,V)-TS U ( S , V ) − T S を S S S について最小にする問題は、U U U を「S S S の値を決めた」という拘束のもとで扱う代わりに、拘束を外して − T S -TS − T S を足した形になっている。第1部第6章 の未定乗数の言葉では、T T T は拘束 S = 一定 S=\text{一定} S = 一定 に付く乗数で、「拘束の値 S S S を1単位変えたときの最適値 U U U の変化」∂ U / ∂ S \partial U/\partial S ∂ U / ∂ S に等しい。Legendre 変換は、拘束の値(S S S )とその乗数(T T T )のどちらを変数として持つかを入れ替える操作と読むこともできる。
体積も取り替えると、エンタルピー H e = U + P V H_{\mathrm e}=U+PV H e = U + P V (V → − P V\to-P V → − P )や Gibbs エネルギー G = F + P V G=F+PV G = F + P V が同じ手続きで得られる。どの変数を取り替えるかの組み合わせと、統計力学の分配関数との関係は第3部第4章 で扱う。
例題 1.7.5 単原子理想気体の自由エネルギー ★★★
オリジナル 単原子理想気体(粒子数 N N N )のエントロピーは、定数 s 0 s_0 s 0 を含めてS ( U , V ) = N k B [ ln ( V / N ) + 3 2 ln ( U / N ) + s 0 ] S(U,V)=N k_{\mathrm B}\bigl[\ln(V/N)+\frac32\ln(U/N)+s_0\bigr] S ( U , V ) = N k B [ ln ( V / N ) + 2 3 ln ( U / N ) + s 0 ] と書くことができる。
(1) これを U U U について解いて U ( S , V ) U(S,V) U ( S , V ) を求め、T = ∂ U / ∂ S T=\partial U/\partial S T = ∂ U / ∂ S から U = 3 2 N k B T U=\frac32N k_{\mathrm B} T U = 2 3 N k B T を示せ。
(2) F ( T , V ) = U − T S F(T,V)=U-TS F ( T , V ) = U − T S を T T T と V V V で表せ。
(3) 式 (1.7.14) から P P P を求め、∂ 2 F / ∂ T 2 = − C V / T \partial^2F/\partial T^2=-C_V/T ∂ 2 F / ∂ T 2 = − C V / T を確かめよ。
解答 (1) 3 2 ln ( U / N ) = S / ( N k B ) − ln ( V / N ) − s 0 \frac32\ln(U/N)=S/(N k_{\mathrm B})-\ln(V/N)-s_0 2 3 ln ( U / N ) = S / ( N k B ) − ln ( V / N ) − s 0 より
U ( S , V ) = N ( N V ) 2 / 3 exp [ 2 3 ( S N k B − s 0 ) ] . U(S,V)=N\Bigl(\frac NV\Bigr)^{2/3}\exp\Bigl[\frac{2}{3}\Bigl(\frac{S}{N k_{\mathrm B}}-s_0\Bigr)\Bigr]. U ( S , V ) = N ( V N ) 2/3 exp [ 3 2 ( N k B S − s 0 ) ] .
S S S で微分すると指数の係数 2 / ( 3 N k B ) 2/(3N k_{\mathrm B}) 2/ ( 3 N k B ) が出るので T = ∂ U / ∂ S = 2 U / ( 3 N k B ) T=\partial U/\partial S=2U/(3N k_{\mathrm B}) T = ∂ U / ∂ S = 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 。∂ 2 U / ∂ S 2 = 2 3 N k B T > 0 \partial^2U/\partial S^2=\frac{2}{3N k_{\mathrm B}}\,T>0 ∂ 2 U / ∂ S 2 = 3 N k B 2 T > 0 なので U U U は S S S について下に凸で、C V = 3 2 N k B C_V=\frac32N k_{\mathrm B} C V = 2 3 N k B と合わせて T / C V T/C_V T / C V に等しい。
(2) U = 3 2 N k B T U=\frac32N k_{\mathrm B} T U = 2 3 N k B T を S S S の式に入れて S = N k B [ ln ( V / N ) + 3 2 ln ( 3 2 k B T ) + s 0 ] S=N k_{\mathrm B}\bigl[\ln(V/N)+\frac32\ln(\frac32 k_{\mathrm B} T)+s_0\bigr] S = N k B [ ln ( V / N ) + 2 3 ln ( 2 3 k B T ) + s 0 ] 。よって
F = 3 2 N k B T − N k B T [ ln V N + 3 2 ln ( 3 2 k B T ) + s 0 ] . F=\frac32N k_{\mathrm B} T-N k_{\mathrm B} T\Bigl[\ln\frac VN+\frac32\ln\Bigl(\frac32 k_{\mathrm B} T\Bigr)+s_0\Bigr]. F = 2 3 N k B T − N k B T [ ln N V + 2 3 ln ( 2 3 k B T ) + s 0 ] .
(3) P = − ∂ F / ∂ V ∣ T = N k B T / V P=-\partial F/\partial V|_T=N k_{\mathrm B} T/V P = − ∂ F / ∂ V ∣ T = N k B T / V で、理想気体の状態方程式が出る。∂ F / ∂ T = 3 2 N k B − N k B [ ⋯ ] − 3 2 N k B = − S \partial F/\partial T=\frac32N k_{\mathrm B}-N k_{\mathrm B}[\cdots]-\frac32N k_{\mathrm B}=-S ∂ F / ∂ T = 2 3 N k B − N k B [ ⋯ ] − 2 3 N k B = − S で式 (1.7.14) の1つ目と合う。もう一度微分すると ∂ 2 F / ∂ T 2 = − ∂ S / ∂ T ∣ V = − 3 2 N k B / T = − C V / T \partial^2F/\partial T^2=-\partial S/\partial T|_V=-\frac32N k_{\mathrm B}/T=-C_V/T ∂ 2 F / ∂ T 2 = − ∂ S / ∂ T ∣ V = − 2 3 N k B / T = − C V / T で、定理 1.7.2 の逆数の関係∂ 2 F ∂ T 2 = − 1 / ∂ 2 U ∂ S 2 \frac{\partial^2F}{\partial T^2}=-1\big/\frac{\partial^2U}{\partial S^2} ∂ T 2 ∂ 2 F = − 1 / ∂ S 2 ∂ 2 U と一致する。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com