tecktom.com

出てくるもの3:凸でない自由エネルギーと共通接線

2026年10月15日8分で読めますAI生成教科書『虎の巻』第3部 第4章 熱力学ポテンシャルは Legendre 変換
目次
  1. 安定性の条件
  2. Legendre 変換は凸包しか見ない
  3. の折れ目が1次相転移
  4. 木下の層

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第3部 第4章 §3.4.5 ★★★ Legendre 変換は凸でない部分を凸包に置き換える。これが気体と液体の共存で、G(P)G(P) の折れ目が1次相転移の体積の跳びになる

安定性の条件

3.4.2節で、F(T,V)F(T,V) は VV について下に凸であることを示した。2階偏微分で書くと

∂2F∂V2∣T=−∂P∂V∣T=1VκT≥0,∂2F∂T2∣V=−CVT≤0(3.4.13)\frac{\partial^2F}{\partial V^2}\Big|_T=-\frac{\partial P}{\partial V}\Big|_T=\frac{1}{V\kappa_T}\ge0,\qquad \frac{\partial^2F}{\partial T^2}\Big|_V=-\frac{C_V}{T}\le0\tag{3.4.13}

で、等温圧縮率 κT=−1V(∂V/∂P)T\kappa_T=-\frac1V(\partial V/\partial P)_T と定積熱容量 CVC_V が負でないことと同じである。κT<0\kappa_T<0 の状態では、外から少し押して体積が減ると圧力も下がるので、さらに押し込まれて潰れていく。凸性は熱力学的な安定性の条件である。

ところが、分子間の引力を平均場で取り込んだ van der Waals の状態方程式(由来は第3部第7章)では、低温で (∂P/∂V)T>0(\partial P/\partial V)_T>0 の区間が現れる。臨界点の値 Tc,PcT_{\mathrm c},P_{\mathrm c} と1分子あたりの体積 vcv_{\mathrm c} で測った換算量 T,P,vT,P,v では

P=8T3v−1−3v2,f(v)=−8T3ln⁡(3v−1)−3v(P=−dfdv)(3.4.14)P=\frac{8T}{3v-1}-\frac{3}{v^2},\qquad f(v)=-\frac{8T}{3}\ln(3v-1)-\frac{3}{v}\quad(P=-\frac{\mathrm{d} f}{\mathrm{d} v})\tag{3.4.14}

で、ff は1分子あたりの自由エネルギー(PcvcP_{\mathrm c}v_{\mathrm c} 単位。TT だけの関数の項は省いた)である。T<1T<1 では f(v)f(v) に上に凸な区間ができる。

Legendre 変換は凸包しか見ない

圧力 PP を外から決めた系の Gibbs エネルギーは g(P)=inf⁡v[f(v)+Pv]g(P)=\inf_v[f(v)+Pv] である。第1部第7章で見たとおり、sup⁡\sup(ここでは inf⁡\inf)で定義した Legendre 変換は、関数の凸でない部分を見ない。g(P)g(P) は、ff をその凸包 f∗∗f^{**}(ff の下に引ける凸関数のうち最大のもの)に置き換えても変わらない。ff の上に凸な区間では、凸包は2点 vl<vgv_{\mathrm l}<v_{\mathrm g} で ff に接する1本の直線、共通接線になる。共通接線の条件は

−f′(vl)=−f′(vg)=Pco,f(vl)+Pcovl=f(vg)+Pcovg(3.4.15)-f'(v_{\mathrm l})=-f'(v_{\mathrm g})=P_{\mathrm{co}},\qquad f(v_{\mathrm l})+P_{\mathrm{co}}v_{\mathrm l}=f(v_{\mathrm g})+P_{\mathrm{co}}v_{\mathrm g}\tag{3.4.15}

で、1つ目は2点の圧力が等しいこと、2つ目は2点の1分子あたりの Gibbs エネルギー(化学ポテンシャル)が等しいことである。2つ目を f(vl)−f(vg)=∫vlvgP dvf(v_{\mathrm l})-f(v_{\mathrm g})=\int_{v_{\mathrm l}}^{v_{\mathrm g}}P\,\mathrm{d} v と書き直すと ∫vlvg(P−Pco) dv=0\int_{v_{\mathrm l}}^{v_{\mathrm g}}(P-P_{\mathrm{co}})\,\mathrm{d} v=0 になり、PP–vv 図で水平線の上下の面積が等しいという Maxwell の等面積則になる。

凸包の直線部分は、体積 vlv_{\mathrm l} の液体と vgv_{\mathrm g} の気体が割合を変えて共存している状態を表す。系の一部を液体、残りを気体にすると、全体の自由エネルギーは2つの値の重み付き平均、すなわち2点を結ぶ線分の上の値になる。元の曲線 f(v)f(v) がこの線分より上にあるなら、一様な状態より2相に分かれた状態のほうが自由エネルギーが低い。3.4.3節の最小原理により、系は2相に分かれる。f∗∗f^{**} は、この分かれ方まで含めて最小にした自由エネルギーである。

図 図

図 3.4.3 換算温度 T=0.9T=0.9 の van der Waals 気体。左:f(v)f(v) に PcovP_{\mathrm{co}}v を足して傾けた図(Pco=0.647P_{\mathrm{co}}=0.647)。共通接線は水平線になり、vl=0.603v_{\mathrm l}=0.603 と vg=2.349v_{\mathrm g}=2.349 で曲線に接する。その間で曲線が直線より上にあるので、凸包は直線部分を持つ。右:g(P)=f(v)+Pvg(P)=f(v)+Pv を vv を動かして描いた曲線。液体の枝と気体の枝が PcoP_{\mathrm{co}} で交わり、g(P)=inf⁡v[f+Pv]g(P)=\inf_v[f+Pv] は2つの枝の低いほうをたどって PcoP_{\mathrm{co}} で折れる。灰色は (∂P/∂v)>0(\partial P/\partial v)>0 の不安定な枝。どの曲線も近似を含まない

G(P)G(P) の折れ目が1次相転移

図 3.4.3 の右のように、g(P)=inf⁡v[f(v)+Pv]g(P)=\inf_v[f(v)+Pv] は、液体の枝と気体の枝の低いほうをたどる。2つの枝は PcoP_{\mathrm{co}} で交わり、そこで gg は折れる。式 (3.4.7) により ∂g/∂P=v\partial g/\partial P=v なので、折れ目の左右で傾き、すなわち1分子あたりの体積が vgv_{\mathrm g} から vlv_{\mathrm l} へ跳ぶ。Gibbs エネルギーの1階導関数が不連続になる転移を 1次相転移と呼ぶ。第1部第7章の例題では (x2−1)2(x^2-1)^2 の Legendre 変換に折れ目ができることを見たが、それと同じ数学が、気体と液体の転移の体積の跳びを記述している。

逆に、凸な関数の Legendre 変換は滑らかである。ff が狭義に凸なら gg は PP について微分可能で、体積は圧力の連続関数になる。T>1T>1 では ff は全域で凸なので、気体と液体の区別はなく、連続的につながる。T→1T\to1 で vlv_{\mathrm l} と vgv_{\mathrm g} が一致して跳びが 00 になる点が臨界点で、臨界点の近くの振る舞いを自由エネルギーの展開で調べるのが第3部第10章の Landau 理論である。

例題 3.4.6 van der Waals 気体の共存線 ★★

オリジナル換算した van der Waals の式 (3.4.14) を使う。

  1. (1) (∂P/∂v)T=0(\partial P/\partial v)_T=0 が 4Tv3=(3v−1)24Tv^3=(3v-1)^2 と同じであることを示し、(∂P/∂v)T>0(\partial P/\partial v)_T>0 の区間が現れるのは T<1T<1 のときに限ることを示せ。T=0.9T=0.9 でその区間の端の vv と、そこでの PP を数値で求めよ。
  2. (2) T=0.9T=0.9 で共通接線の条件 (3.4.15) を数値で解き、Pco,vl,vgP_{\mathrm{co}},v_{\mathrm l},v_{\mathrm g} を求めよ。等面積則が成り立つことも確かめよ。
  3. (3) 転移での1分子あたりの体積の跳びを求めよ。
解答

(1) ∂P/∂v=−24T(3v−1)2+6v3=6[(3v−1)2−4Tv3]v3(3v−1)2\partial P/\partial v=-\dfrac{24T}{(3v-1)^2}+\dfrac{6}{v^3}=\dfrac{6\bigl[(3v-1)^2-4Tv^3\bigr]}{v^3(3v-1)^2} で、分母は v>1/3v>1/3 で正なので、符号は (3v−1)2−4Tv3(3v-1)^2-4Tv^3 で決まる。∂P/∂v>0\partial P/\partial v>0 は T<h(v)≔(3v−1)2/(4v3)T<h(v)\coloneqq(3v-1)^2/(4v^3) と同じである。h′(v)=3(3v−1)(1−v)4v4h'(v)=\dfrac{3(3v-1)(1-v)}{4v^4} なので、hh は v=1v=1 で最大値 h(1)=1h(1)=1 を取る。したがって T<1T<1 のときに限り、h(v)>Th(v)>T となる vv の区間が現れる。T=0.9T=0.9 で 4Tv3=(3v−1)24Tv^3=(3v-1)^2 を数値で解くと v=0.719v=0.719 と 1.5291.529、そこでの圧力は P=0.420P=0.420(極小)と 0.7240.724(極大)である。

(2) 3つの未知数に3つの式(P(vl)=PcoP(v_{\mathrm l})=P_{\mathrm{co}}、P(vg)=PcoP(v_{\mathrm g})=P_{\mathrm{co}}、f(vl)+Pcovl=f(vg)+Pcovgf(v_{\mathrm l})+P_{\mathrm{co}}v_{\mathrm l}=f(v_{\mathrm g})+P_{\mathrm{co}}v_{\mathrm g})を Newton 法で解くと

Pco=0.647,vl=0.603,vg=2.349.P_{\mathrm{co}}=0.647,\qquad v_{\mathrm l}=0.603,\qquad v_{\mathrm g}=2.349 .

PcoP_{\mathrm{co}} は (1) の極小 0.4200.420 と極大 0.7240.724 の間にある。∫vlvg(P−Pco) dv\int_{v_{\mathrm l}}^{v_{\mathrm g}}(P-P_{\mathrm{co}})\,\mathrm{d} v を数値積分すると 00(10−1110^{-11} 以下)になり、等面積則が成り立つ。

(3) vg−vl=1.745v_{\mathrm g}-v_{\mathrm l}=1.745(臨界体積の単位)で、気体は液体の 3.893.89 倍の体積を持つ。この跳びは g(P)g(P) の PcoP_{\mathrm{co}} での傾きの差である。

木下の勉強ステータス

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

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

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