tecktom.com

この章のまとめ

2026年10月15日11分で読めますAI生成教科書『虎の巻』第3部 第4章 熱力学ポテンシャルは Legendre 変換

AI 生成・木下が未勉強

解析力学から様々な理論分野を組み立てる 第3部 第4章 §3.4.7 ★★★ 熱力学ポテンシャルは U(S,V,N)U(S,V,N) のどの変数を取り替えたかで決まり、分配関数はその Legendre 変換を Laplace 変換として実行している

この章で選んだものは、状態数から得た基本関係式 U(S,V,N)U(S,V,N) である。UU は1次同次で、Euler の関係 U=TS−PV+μNU=TS-PV+\mu N と Gibbs–Duhem の関係を満たし、示量変数について下に凸である。当てた道具は第1部第7章の Legendre 変換で、3つの示量変数のどれを傾きに取り替えるかで、FF・HeH_{\mathrm e}・GG・Ω\Omega が1つずつ決まる。取り替えた変数の傾きは −xa-x_a に、残した変数の傾きはそのまま引き継がれ、2階偏微分の対称性から Maxwell 関係式が出る。示量変数を全部取り替えると 00 しか残らず、1つだけ残すと G=μNG=\mu N、Ω=−PV\Omega=-PV のようにその変数の1次式になる。inf⁡\inf は熱浴や圧力の環境に接した系の最小原理そのものである。統計力学の分配関数は、状態数の Laplace 変換で、粒子数が大きいと被積分関数の最大値で決まり、その最大値が Legendre 変換になる。アンサンブルを取り替えることと、熱力学ポテンシャルを取り替えることは、同じ操作である。凸でない自由エネルギーは、Legendre 変換の段階で凸包に置き換わり、共通接線が気体と液体の共存を、G(P)G(P) の折れ目が1次相転移の体積の跳びを表す。次の章(第3部第5章)では、エントロピーを一定に保つ断熱変化を、第1部第14章の断熱不変量と結びつける。

この章の道具

  • 基本関係式 U(S,V,N)U(S,V,N):T=∂U/∂ST=\partial U/\partial S、P=−∂U/∂VP=-\partial U/\partial V、μ=∂U/∂N\mu=\partial U/\partial N。1次同次から U=TS−PV+μNU=TS-PV+\mu N、S dT−V dP+N dμ=0S\,\mathrm{d} T-V\,\mathrm{d} P+N\,\mathrm{d}\mu=0。
  • ポテンシャル=取り替える変数の組:F=inf⁡S[U−TS]F=\inf_S[U-TS]、He=inf⁡V[U+PV]H_{\mathrm e}=\inf_V[U+PV]、G=inf⁡S,V[U−TS+PV]=μNG=\inf_{S,V}[U-TS+PV]=\mu N、Ω=inf⁡S,N[U−TS−μN]=−PV\Omega=\inf_{S,N}[U-TS-\mu N]=-PV。
  • 偏微分の規則:∂Φ/∂ya=−xa\partial\Phi/\partial y_a=-x_a、∂Φ/∂xb=yb\partial\Phi/\partial x_b=y_b。混合偏微分の対称性が Maxwell 関係式。
  • 鞍点:−kBTln⁡∫e(S(E)−E/T)/kBdE/ΔE=inf⁡E[E−TS(E)]+O(kBTln⁡N)- k_{\mathrm B} T\ln\int\mathrm{e}^{(S(E)-E/T)/ k_{\mathrm B}}\mathrm{d} E/\Delta E=\inf_E[E-TS(E)]+\mathcal{O}( k_{\mathrm B} T\ln N)。Ξ→Ω\Xi\to\Omega、圧力アンサンブル →G\to G も同じ。
  • 凸でない ff の Legendre 変換は凸包 f∗∗f^{**} しか見ない:共通接線・等面積則・GG の折れ目=1次相転移。

総合問題 常磁性体の自由エネルギーを磁場と磁化で取り替える ★★★

オリジナル磁場 BB の中に、互いに相互作用しない NN 個のスピン 1/21/2 を置く。各スピンの磁気モーメントは磁場の向きに +μB+\mu_{\mathrm B} か −μB-\mu_{\mathrm B} を取り、エネルギーはそれぞれ ∓μBB\mp\mu_{\mathrm B}B である。全磁気モーメントを MM、x=μBB/(kBT)x=\mu_{\mathrm B}B/( k_{\mathrm B} T)、m=M/(NμB)m=M/(N\mu_{\mathrm B}) と書く。

  1. (1) (使う道具:第3部第3章のカノニカル分布)分配関数 ZZ と自由エネルギー F(T,B)=−kBTln⁡ZF(T,B)=- k_{\mathrm B} T\ln Z を求めよ。
  2. (2) (使う道具:3.4.2節の偏微分の規則)dF=−S dT−M dB\mathrm{d} F=-S\,\mathrm{d} T-M\,\mathrm{d} B として、M(T,B)M(T,B) と S(T,B)S(T,B) を求めよ。
  3. (3) (使う道具:第1部第7章・3.4.2節)BB を MM に取り替えた F~(T,M)=sup⁡B[F(T,B)+MB]\tilde F(T,M)=\sup_B[F(T,B)+MB] を求め、F~=−NkBT hb(1+m2)\tilde F=-N k_{\mathrm B} T\,h_{\mathrm b}\bigl(\frac{1+m}{2}\bigr)(hb(p)=−pln⁡p−(1−p)ln⁡(1−p)h_{\mathrm b}(p)=-p\ln p-(1-p)\ln(1-p) は二値エントロピー)となることと、∂F~/∂M=B\partial\tilde F/\partial M=B を示せ。
  4. (4) (使う道具:3.4.3節の Maxwell 関係式)(∂S/∂B)T=(∂M/∂T)B(\partial S/\partial B)_T=(\partial M/\partial T)_B を、(2) の結果で確かめよ。
  5. (5) (使う道具:(2)・3.4.3節)温度 1.0 K1.0\,\mathrm K、磁場 1.0 T1.0\,\mathrm T で磁化させたあと、熱を出入りさせずに(エントロピー一定で)磁場を 0.010 T0.010\,\mathrm T まで下げる。最後の温度を求めよ(スピン間の相互作用は無視する)。はじめの状態の S/(NkB)S/(N k_{\mathrm B}) も求めよ。μB=9.274×10−24 J T−1\mu_{\mathrm B}=9.274\times10^{-24}\,\mathrm{J\,T^{-1}}、kB=1.381×10−23 J K−1k_{\mathrm B}=1.381\times10^{-23}\,\mathrm{J\,K^{-1}}。
  6. (6) (使う道具:3.4.4節の鞍点)ZZ を上向きのスピンの数 nn についての和 Z=∑n(Nn)ex(2n−N)Z=\sum_n\binom Nn\mathrm{e}^{x(2n-N)} と見て、N=100N=100、x=0.5x=0.5 で ln⁡Z\ln Z と最大項の対数を比べよ。N→∞N\to\infty では、最大項の対数を NN で割ったものが sup⁡m[hb(1+m2)+xm]=ln⁡(2cosh⁡x)\sup_m\bigl[h_{\mathrm b}\bigl(\frac{1+m}2\bigr)+xm\bigr]=\ln(2\cosh x) に近づくことを示せ。
解答

(1) 1つのスピンの分配関数は ex+e−x=2cosh⁡x\mathrm{e}^{x}+\mathrm{e}^{-x}=2\cosh x で、スピンは独立なので Z=(2cosh⁡x)NZ=(2\cosh x)^N。F(T,B)=−NkBTln⁡(2cosh⁡x)F(T,B)=-N k_{\mathrm B} T\ln(2\cosh x)。

(2) ∂x∂B=μBkBT\frac{\partial x}{\partial B}=\dfrac{\mu_{\mathrm B}}{ k_{\mathrm B} T} を使って M=−∂F∂B∣T=NμBtanh⁡xM=-\frac{\partial F}{\partial B}\big|_T=N\mu_{\mathrm B}\tanh x、すなわち m=tanh⁡xm=\tanh x。∂x∂T=−xT\frac{\partial x}{\partial T}=-\dfrac xT を使って

S=−∂F∂T∣B=NkBln⁡(2cosh⁡x)+NkBTtanh⁡x ∂x∂T=NkB[ln⁡(2cosh⁡x)−xtanh⁡x].S=-\frac{\partial F}{\partial T}\Big|_B=N k_{\mathrm B}\ln(2\cosh x)+N k_{\mathrm B} T\tanh x\,\frac{\partial x}{\partial T}=N k_{\mathrm B}\bigl[\ln(2\cosh x)-x\tanh x\bigr].

ここで FF は磁場 BB(示強変数)を変数に持つので、3.4.2節の言葉では、全磁気モーメント MM をすでに BB に取り替えたポテンシャルにあたる。−M dB-M\,\mathrm{d} B の負号は、取り替えた変数の傾きが −xa-x_a になる規則と同じである。

(3) sup⁡B[F+MB]\sup_B[F+MB] の最大点の条件は ∂F/∂B+M=0\partial F/\partial B+M=0、すなわち (2) から m=tanh⁡xm=\tanh x、B=kBTμBartanh⁡mB=\dfrac{ k_{\mathrm B} T}{\mu_{\mathrm B}}\operatorname{artanh}m。そこでの値は

F~=F+MB=−NkBT[ln⁡(2cosh⁡x)−xtanh⁡x]=−TS.\tilde F=F+MB=-N k_{\mathrm B} T\bigl[\ln(2\cosh x)-x\tanh x\bigr]=-TS .

cosh⁡x=1/1−m2\cosh x=1/\sqrt{1-m^2}、x=12ln⁡1+m1−mx=\frac12\ln\frac{1+m}{1-m} を入れて整理すると

ln⁡(2cosh⁡x)−xtanh⁡x=ln⁡2−1+m2ln⁡(1+m)−1−m2ln⁡(1−m)=hb(1+m2)\ln(2\cosh x)-x\tanh x=\ln2-\frac{1+m}{2}\ln(1+m)-\frac{1-m}{2}\ln(1-m)=h_{\mathrm b}\Bigl(\frac{1+m}2\Bigr)

で、F~=−NkBT hb(1+m2)\tilde F=-N k_{\mathrm B} T\,h_{\mathrm b}\bigl(\frac{1+m}2\bigr)。最後の等号は、p=(1+m)/2p=(1+m)/2 を hbh_{\mathrm b} に入れて展開すると確かめることができる。F~\tilde F を MM で微分すると hbh_{\mathrm b} の導関数 ddmhb=−artanh⁡m\frac{\mathrm{d}}{\mathrm{d} m}h_{\mathrm b}=-\operatorname{artanh}m から ∂F~∂M=kBTμBartanh⁡m=B\frac{\partial \tilde F}{\partial M}=\dfrac{ k_{\mathrm B} T}{\mu_{\mathrm B}}\operatorname{artanh}m=B で、逆向きの変換の傾きが BB に戻る。磁化を決めた系の自由エネルギーは −T×-T\times(磁化 mm と両立するスピンの並べ方の数の対数)で、エネルギーの項を含まない。相互作用のないスピンでは、磁化の値を決めた時点でエネルギーは −MB-MB に決まり、それを F~\tilde F から外して FF に移したのが Legendre 変換である。

(4) (2) の SS を xx で微分すると ddx[ln⁡(2cosh⁡x)−xtanh⁡x]=tanh⁡x−tanh⁡x−xcosh⁡2x=−xcosh⁡2x\frac{\mathrm{d} }{\mathrm{d} x}\bigl[\ln(2\cosh x)-x\tanh x\bigr]=\tanh x-\tanh x-\frac{x}{\cosh^2x}=-\frac{x}{\cosh^2x}。∂x∂B=xB\frac{\partial x}{\partial B}=\frac xB なので∂S∂B∣T=−NkBx2Bcosh⁡2x=−NμBxTcosh⁡2x\frac{\partial S}{\partial B}\big|_T=-\dfrac{N k_{\mathrm B} x^2}{B\cosh^2x}=-\dfrac{N\mu_{\mathrm B}x}{T\cosh^2x}。一方 ∂M∂T∣B=NμB1cosh⁡2x⋅(−xT)\frac{\partial M}{\partial T}\big|_B=N\mu_{\mathrm B}\dfrac{1}{\cosh^2x}\cdot\Bigl(-\dfrac xT\Bigr) で、2つは等しい。F(T,B)F(T,B) の混合偏微分の対称性で、式 (3.4.9) の2本目の (V,P)(V,P) を (B,−M)(B,-M) に読み替えたものである。

(5) (2) の SS は x=μBB/(kBT)x=\mu_{\mathrm B}B/( k_{\mathrm B} T) だけの関数なので、SS を一定に保つ変化では B/TB/T が一定に保たれる。したがって

Tf=TiBfBi=1.0 K×0.0101.0=0.010 K.T_{\mathrm f}=T_{\mathrm i}\frac{B_{\mathrm f}}{B_{\mathrm i}}=1.0\,\mathrm K\times\frac{0.010}{1.0}=0.010\,\mathrm K .

はじめの状態では xi=9.274×10−24×1.01.381×10−23×1.0=0.672x_{\mathrm i}=\dfrac{9.274\times10^{-24}\times1.0}{1.381\times10^{-23}\times1.0}=0.672 で、S/(NkB)=ln⁡(2cosh⁡0.672)−0.672tanh⁡0.672=0.510S/(N k_{\mathrm B})=\ln(2\cosh0.672)-0.672\tanh0.672=0.510。磁場のない無秩序な状態の ln⁡2=0.693\ln2=0.693 の 0.7360.736 倍で、磁場でスピンをそろえた分だけエントロピーが下がっている。磁場を下げてもエントロピーを保つには、スピンのそろい方 m=tanh⁡xm=\tanh x を保つ必要があり、そのために温度が磁場と同じ割合で下がる。(4) の Maxwell 関係式 (∂S/∂B)T<0(\partial S/\partial B)_T<0 が、この冷却の向き(等温で磁場をかけるとエントロピーが減る)を表している。実際の物質ではスピン間の相互作用や内部の磁場があるので、温度は Bf→0B_{\mathrm f}\to0 でも 00 にはならない。

(6) ln⁡Z=100ln⁡(2cosh⁡0.5)=81.33\ln Z=100\ln(2\cosh0.5)=81.33。和の各項の対数 ln⁡(100n)+0.5(2n−100)\ln\binom{100}{n}+0.5(2n-100) を n=0,…,100n=0,\dots,100 で比べると、最大は n=73n=73(m=0.46m=0.46、tanh⁡0.5=0.462\tanh0.5=0.462 に最も近い)で、その値は 78.9178.91。差は 2.412.41 で、ln⁡Z\ln Z の 3.0%3.0\% である。この差は最大項のまわりの山の幅の対数で、二項分布の幅から 12ln⁡(2πNp(1−p))=2.41\frac12\ln\bigl(2\pi Np(1-p)\bigr)=2.41(p=(1+tanh⁡x)/2p=(1+\tanh x)/2)と見積もることができる。N→∞N\to\infty では Stirling の公式により 1Nln⁡(Nn)→hb(n/N)\frac1N\ln\binom Nn\to h_{\mathrm b}(n/N) なので、最大項の対数を NN で割ったものは sup⁡m[hb(1+m2)+xm]\sup_m[h_{\mathrm b}(\frac{1+m}2)+xm] に近づく。最大点の条件 −artanh⁡m+x=0-\operatorname{artanh}m+x=0 から m=tanh⁡xm=\tanh x で、その値は (3) の計算を逆にたどると hb(1+m2)+xm=ln⁡(2cosh⁡x)−xtanh⁡x+xtanh⁡x=ln⁡(2cosh⁡x)h_{\mathrm b}(\frac{1+m}{2})+xm=\ln(2\cosh x)-x\tanh x+x\tanh x=\ln(2\cosh x) である。これは (1) の厳密な 1Nln⁡Z\frac1N\ln Z に等しい。和を最大項で置き換えることが、mm についての Legendre 変換((3) の逆向き)と同じ計算になっている。

木下の勉強ステータス

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

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

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