tecktom.com

この章で選ぶもの:基本関係式 U(S,V,N)U(S,V,N)

2026年10月15日9分で読めますAI生成教科書『虎の巻』第3部 第4章 熱力学ポテンシャルは Legendre 変換
目次
  1. ミクロな から へ
  2. 1次同次性と Euler の関係
  3. 凸性
  4. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第3部 第4章 §3.4.1 ★★★ 状態数から得た S(E,V,N)S(E,V,N) を EE について解いた U(S,V,N)U(S,V,N) が出発点。1次同次性と凸性の2つの性質を持つ

熱力学の教科書には、内部エネルギー UU・Helmholtz の自由エネルギー FF・エンタルピー HeH_{\mathrm e}・Gibbs エネルギー GG・グランドポテンシャル Ω\Omega と、エネルギーの次元を持つ関数がいくつも並ぶ。それぞれに「自然な変数」があり、全微分の形と Maxwell 関係式を覚えることになっている。第1部第7章では、U(S,V)U(S,V) のエントロピー SS を温度 T=∂U/∂ST=\partial U/\partial S に取り替える Legendre 変換として F=U−TSF=U-TS を作った。この章では、同じ操作を U(S,V,N)U(S,V,N) の3つの変数のどれを取り替えるかの組み合わせとして整理し、熱力学ポテンシャルがすべてその組み合わせの1つずつであることを示す。

この章で選ぶものは、ミクロな Hamiltonian HH から状態数を数えて得られる1つの関数、基本関係式 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節で常磁性体を題材にした総合問題を解く。

ミクロな HH から U(S,V,N)U(S,V,N) へ

第3部第2章では、Hamiltonian H(q,p)H(q,p) を持つ NN 粒子の系を体積 VV の容器に入れ、エネルギーが幅 ΔE\Delta E の殻 [E,E+ΔE)[E,E+\Delta E) に入る状態の数 W(E)W(E) を数えて、エントロピーを S(E,V,N)=kBln⁡W(E)S(E,V,N)= k_{\mathrm B}\ln W(E) と定めた。W(E)W(E) はふつう EE について単調に増えるので、S=S(E,V,N)S=S(E,V,N) を EE について解くことができる。その解を内部エネルギー U(S,V,N)U(S,V,N) と書き、これを基本関係式と呼ぶ。この章で選ぶものはこの1つの関数で、熱力学のすべての量はこの関数の偏微分から出てくる。

UU の1階偏微分に名前を付ける。

T=∂U∂S∣V,N,P=−∂U∂V∣S,N,μ=∂U∂N∣S,V,dU=T dS−P dV+μ dN.(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}

PP にだけ負号が付くのは、体積を増やすと系が外へ仕事をしてエネルギーが減る向きを正の圧力と呼ぶ習慣による。(S,V,N)(S,V,N) は系を2つ並べると2倍になる量で示量変数と呼び、それぞれの傾き (T,−P,μ)(T,-P,\mu) は2つ並べても変わらない量で示強変数と呼ぶ。示量変数と、それに対応する傾きの組を共役な組と呼ぶ。

1次同次性と Euler の関係

短距離の相互作用しかない系では、同じ系を2つ並べてつなげても、境界での相互作用のエネルギーは全体に比べて小さい。したがって S,V,NS,V,N をすべて λ\lambda 倍すると UU も λ\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}

この性質を1次同次性と呼ぶ。両辺を λ\lambda で微分して λ=1\lambda=1 と置くと、左辺は S ∂U/∂S+V ∂U/∂V+N ∂U/∂NS\,\partial U/\partial S+V\,\partial U/\partial V+N\,\partial U/\partial N になるので

U=TS−PV+μN(3.4.3)U=TS-PV+\mu N\tag{3.4.3}

を得る。これを Euler の関係と呼ぶ。式 (3.4.3) の全微分を取り、式 (3.4.1) の dU\mathrm{d} U を引くと

S dT−V dP+N dμ=0(3.4.4)S\,\mathrm{d} T-V\,\mathrm{d} P+N\,\mathrm{d}\mu=0\tag{3.4.4}

が残る。これを Gibbs–Duhem の関係と呼ぶ。3つの示強変数 (T,P,μ)(T,P,\mu) は独立に動けず、2つを決めると残りの1つが決まる。この関係は3.4.2節で「3つの変数を全部取り替えると何も残らない」という形で戻ってくる。

凸性

UU は (S,V,N)(S,V,N) について下に凸である。理由は次のとおりである。系を壁で2つに仕切り、左に (U1,V1,N1)(U_1,V_1,N_1)、右に (U2,V2,N2)(U_2,V_2,N_2) を入れる。壁を取り去ると、取り去る前に数えた状態はすべて取り去った後にも許されるので、状態数は減らない。したがってS(U1+U2,V1+V2,N1+N2)≥S(U1,V1,N1)+S(U2,V2,N2)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) である。右辺を1次同次性で書き直すと、2点の中点で SS が2点の値の平均以上になること、すなわち SS が上に凸であることになる。SS が UU について増加関数で上に凸なら、その逆関数 UU は SS について下に凸である(V,NV,N についても同じ議論が通る)。凸性は3.4.2節で Legendre 変換が情報を落とさない条件として使い、3.4.5節では、それが破れて見えるときに何が起こるかを調べる。

1次同次な関数は、原点を通る直線の上で1次式になるので、その方向には曲がっていない。UU の (S,V,N)(S,V,N) についての Hesse 行列は、(S,V,N)(S,V,N) 自身の方向に固有値 00 を持ち、残りの2方向で正の固有値を持つ。狭義に凸なのは、たとえば NN を止めた (S,V)(S,V) の2変数についてである。

例題 3.4.1 単原子理想気体の基本関係式 ★★★

オリジナル単原子理想気体(粒子の質量 mm)のエントロピーは、Planck 定数 hh を含めてS=NkB[ln⁡{VN(4πmU3Nh2)3/2}+52]S=N k_{\mathrm B}\Bigl[\ln\Bigl\{\dfrac VN\Bigl(\dfrac{4\pi mU}{3Nh^2}\Bigr)^{3/2}\Bigr\}+\dfrac52\Bigr]と書くことができる(Sackur–Tetrode の式。導出は第3部第2章)。

  1. (1) これを UU について解き、U(λS,λV,λN)=λUU(\lambda S,\lambda V,\lambda N)=\lambda U を確かめよ。
  2. (2) T,P,μT,P,\mu を求め、Euler の関係 (3.4.3) を確かめよ。
  3. (3) μ\mu を TT と PP で表し、Gibbs–Duhem の関係 (3.4.4) を確かめよ。
解答

(1) SNkB−52=ln⁡VN+32ln⁡4πmU3Nh2\dfrac{S}{N k_{\mathrm B}}-\dfrac52=\ln\dfrac VN+\dfrac32\ln\dfrac{4\pi mU}{3Nh^2} を UU について解くと

U(S,V,N)=3h2N5/34πmV2/3exp⁡(2S3NkB−53).U(S,V,N)=\frac{3h^2N^{5/3}}{4\pi mV^{2/3}}\exp\Bigl(\frac{2S}{3N k_{\mathrm B}}-\frac53\Bigr).

S,V,NS,V,N を λ\lambda 倍すると、指数の中の S/NS/N は変わらず、前の因子 N5/3V−2/3N^{5/3}V^{-2/3} は λ5/3−2/3=λ\lambda^{5/3-2/3}=\lambda 倍になる。よって UU は λ\lambda 倍になる。

(2) SS で微分すると指数の係数 2/(3NkB)2/(3N k_{\mathrm B}) が出るので T=2U/(3NkB)T=2U/(3N k_{\mathrm B})、すなわち U=32NkBTU=\frac32N k_{\mathrm B} T。VV で微分すると V−2/3V^{-2/3} から −23V-\frac{2}{3V} が出るので P=2U/(3V)P=2U/(3V)。2つを合わせると PV=NkBTPV=N k_{\mathrm B} T。NN で微分すると、前の因子から 53N\frac{5}{3N}、指数から −2S3N2kB-\frac{2S}{3N^2 k_{\mathrm B}} が出るので μ=U(53N−2S3N2kB)\mu=U\bigl(\frac{5}{3N}-\frac{2S}{3N^2 k_{\mathrm B}}\bigr)。これらを並べると

TS−PV+μN=2US3NkB−2U3+5U3−2US3NkB=UTS-PV+\mu N=\frac{2US}{3N k_{\mathrm B}}-\frac{2U}{3}+\frac{5U}{3}-\frac{2US}{3N k_{\mathrm B}}=U

で、Euler の関係が成り立つ。

(3) (2) の μ\mu に U=32NkBTU=\frac32N k_{\mathrm B} T と Sackur–Tetrode の SS を入れて整理すると、熱的 de Broglie 波長 λT=h/2πmkBT\lambda_T=h/\sqrt{2\pi m k_{\mathrm B} T} を使って

μ=kBTln⁡(NVλT3)=kBTln⁡PλT3kBT\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}

となる。λT3∝T−3/2\lambda_T^3\propto T^{-3/2} なので μ=kBT[ln⁡P−52ln⁡T+定数]\mu= k_{\mathrm B} T\bigl[\ln P-\frac52\ln T+\text{定数}\bigr] で、∂μ∂P∣T=kBT/P=V/N\frac{\partial \mu}{\partial P}\big|_T= k_{\mathrm B} T/P=V/N、∂μ∂T∣P=μ/T−52kB=−S/N\frac{\partial \mu}{\partial T}\big|_P=\mu/T-\frac52 k_{\mathrm B}=-S/N(最後の等号は Sackur–Tetrode の式を μ\mu で書き直したもの)。したがって dμ=−(S/N)dT+(V/N)dP\mathrm{d}\mu=-(S/N)\mathrm{d} T+(V/N)\mathrm{d} P で、両辺に NN を掛けると式 (3.4.4) になる。

木下の勉強ステータス

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

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

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