tecktom.com

出てくるもの1:第二法則の3つの形

2026年10月15日11分で読めますAI生成教科書『虎の巻』第4部 第3章 第二法則は経路の D の非負性
目次
  1. 平衡から始めて平衡へ戻す:
  2. 終わりの分布から逆回す:Clausius の不等式
  3. 2つの選び方の関係
  4. 木下の層

AI 生成・木下が未勉強

相対エントロピーから情報と物理を組み立てる 第4部 第3章 §4.3.3 ★★★ 逆回しの初めの分布を平衡にとると β(⟨W⟩−ΔF)\beta(\langle W \rangle -\Delta F)、順回しの終わりの分布にとると Clausius の不等式、どちらも D≥0D\ge0

式 (4.3.7) の中で、選び方が残っているのは逆回しの初めの分布 p~0\tilde p_0 だけである。どの選び方でも定理 4.3.3 から ⟨Σ⟩≥0\langle \Sigma \rangle \ge0 が出るが、何について述べた不等式になるかは選び方で変わる。代表的な2つを見る。

平衡から始めて平衡へ戻す:⟨W⟩≥ΔF\langle W \rangle \ge\Delta F

順回しを平衡 p0=peqλ0p_0=p_{\mathrm{eq}}^{\lambda_0} から始め、逆回しも平衡 p~0=peqλn\tilde p_0=p_{\mathrm{eq}}^{\lambda_n} から始める。−ln⁡peqλ(x)=βEλ(x)−βFλ-\ln p_{\mathrm{eq}}^\lambda(x)=\beta E_\lambda(x)-\beta F_\lambda(Fλ=−kBTln⁡ZλF_\lambda=- k_{\mathrm B} T\ln Z_\lambda)を式 (4.3.7) に入れ、経路ごとの第一法則 Eλn(xn)−Eλ0(x0)=W+QE_{\lambda_n}(x_n)-E_{\lambda_0}(x_0)=W+Q を使うと

Σ[Γ]=β(Eλn(xn)−Eλ0(x0))−βΔF−βQ[Γ]=β(W[Γ]−ΔF),ΔF≔Fλn−Fλ0(4.3.13)\Sigma[\Gamma]=\beta\bigl(E_{\lambda_n}(x_n)-E_{\lambda_0}(x_0)\bigr)-\beta\Delta F-\beta Q[\Gamma]=\beta\bigl(W[\Gamma]-\Delta F\bigr),\qquad\Delta F\coloneqq F_{\lambda_n}-F_{\lambda_0}\tag{4.3.13}

になる。このとき、定理 4.3.3 と式 (4.3.9) はそれぞれ

β(⟨W⟩−ΔF)=D(P ∥ P~)≥0,⟨e−βW⟩=e−βΔF(4.3.14)\beta\bigl(\langle W \rangle -\Delta F\bigr)=D(P\,\|\,\tilde P)\ge0,\qquad\langle \mathrm{e}^{-\beta W} \rangle =\mathrm{e}^{-\beta\Delta F}\tag{4.3.14}

になる。前者は等温の操作の第二法則(系にする仕事は自由エネルギーの差以上、取り出せる仕事は自由エネルギーの減少以下)、後者は Jarzynski 等式1C. Jarzynski (1997), Phys. Rev. Lett. 78, 2690–2693。下の脚注へである。この選び方では、逆回しの過程も平衡から始まるので、順と逆が対等になり、Σ~[Γ~]=−Σ[Γ]\tilde\Sigma[\tilde\Gamma]=-\Sigma[\Gamma] が成り立って Crooks の関係 (4.3.11) が使える2G. E. Crooks (1999), Phys. Rev. E 60, 2721–2726。下の脚注へ。手順の初めと終わりで λ\lambda が同じ(λn=λ0\lambda_n=\lambda_0、サイクル)なら ΔF=0\Delta F=0 で、⟨W⟩≥0\langle W \rangle \ge0、すなわち1つの熱浴だけを使うサイクルで平均として仕事を取り出すことはできない。これは Kelvin の原理である。

⟨W⟩−ΔF\langle W \rangle -\Delta F は散逸した仕事と呼ばれ、それが kBTk_{\mathrm B} T 掛ける経路の DD に正確に等しい、というのが式 (4.3.14) の主張である。散逸した仕事は、逆回しの過程の経路を順回しの過程の経路と見分ける手がかりの量に比例する。

図 4.3.2 は、βε\beta\varepsilon を 00 から ln⁡3\ln3 へ nn 回に分けて上げる操作で、⟨Σ⟩=β(⟨W⟩−ΔF)\langle \Sigma \rangle =\beta(\langle W \rangle -\Delta F) が nn とともに減る様子である。n=1n=1 の値 0.14380.1438 は例題 4.3.2 の値で、nn を増やすとほぼ 1/n1/n に比例して 00 に近づき、n=400n=400 で n⟨Σ⟩≈0.206n\langle \Sigma \rangle \approx0.206 である。ゆっくり動かすほど、各刻みでの λ\lambda の変化が小さくなり、系がその都度 Gibbs 分布に近い所にとどまるので、順と逆の過程の経路が見分けにくくなる。n→∞n\to\infty の極限が準静的な操作である。

図

図 4.3.2 2準位の βε\beta\varepsilon を 00 から ln⁡3\ln3 へ nn 回に分けて上げる操作(各回の後に1回跳ぶ。遷移は例題 4.3.1 と同じ)。上の線は経路の DD =散逸した仕事 /kBT/ k_{\mathrm B} T、下の線は終わりの分布と終わりの Gibbs 分布の DD(4.3.4 節の下からの評価)。どちらも厳密な計算(scripts/verify_ch4-03.py)。

終わりの分布から逆回す:Clausius の不等式

初めの分布 p0p_0 は平衡でなくてもよいとし、逆回しの初めの分布を順回しの終わりの分布 p~0=pn\tilde p_0=p_n にとる。式 (4.3.7) を平均すると

⟨Σ⟩=H(pn)−H(p0)−β⟨Q⟩=D(P ∥ P~)≥0(4.3.15)\langle \Sigma \rangle =H(p_n)-H(p_0)-\beta\langle Q \rangle =D(P\,\|\,\tilde P)\ge0\tag{4.3.15}

になる。kBk_{\mathrm B} を掛けると kBΔH−⟨Q⟩/T≥0k_{\mathrm B}\Delta H-\langle Q \rangle /T\ge0 で、系のエントロピーの増加は、系がもらった熱を温度で割ったもの以上である。これは Clausius の不等式で、系のエントロピーとして Shannon エントロピー(に kBk_{\mathrm B} を掛けたもの)を使えば、初めと終わりが平衡でなくても成り立つ。第2部第4章 では、系と熱浴を合わせた全体の可逆な運動から Σ=ΔH−βQ=I(S;B)+D(ρB ∥ ρBcan)\Sigma=\Delta H-\beta Q=I(S;B)+D(\rho_B\,\|\,\rho_B^{\mathrm{can}}) を出した。そこでは増えたエントロピーが「系と熱浴の相関」と「熱浴の平衡からのずれ」に分かれたが、この章では熱浴をあらわに書かず、同じ量を経路の確率の DD として得た。

式 (4.3.15) を仕事で書き直すことができる。非平衡の自由エネルギー F[p]≔⟨Eλ⟩p−kBT H(p)F[p]\coloneqq\langle E_\lambda \rangle _p- k_{\mathrm B} T\,H(p)(第2部第2章)を使い、⟨Q⟩=⟨ΔE⟩−⟨W⟩\langle Q \rangle =\langle \Delta E \rangle -\langle W \rangle を入れると

⟨W⟩−(F[pn]−F[p0])=kBT D(P ∥ P~)≥0(4.3.16)\langle W \rangle -\bigl(F[p_n]-F[p_0]\bigr)= k_{\mathrm B} T\,D(P\,\|\,\tilde P)\ge0\tag{4.3.16}

になる(F[pn]F[p_n] は λn\lambda_n、F[p0]F[p_0] は λ0\lambda_0 のエネルギーで計算する)。初めと終わりが平衡でないときの第二法則で、Esposito と Van den Broeck が示した形である(第2部第2章で引用した)。

2つの選び方の関係

p~0\tilde p_0 の2つの選び方は、xnx_n の分布だけが違い、xnx_n を与えたときの残りの経路の条件付き分布は同じである。そこで DD の連鎖律(第1部第3章)を xnx_n で分けて当てると、平衡から始めた場合に

β(⟨W⟩−ΔF)=D(P ∥ P~p~0=pn)⏟全体のエントロピー生成+D(pn ∥ peqλn)(4.3.17)\beta\bigl(\langle W \rangle -\Delta F\bigr)=\underbrace{D\bigl(P\,\|\,\tilde P_{\tilde p_0=p_n}\bigr)}_{\text{全体のエントロピー生成}}+D\bigl(p_n\,\|\,p_{\mathrm{eq}}^{\lambda_n}\bigr)\tag{4.3.17}

が成り立つ。散逸した仕事のうち、第2項は操作が終わった瞬間にまだ平衡に達していない分(自由エネルギーの余り。この後放っておけば熱として失われるが、うまく使えば取り出せる)で、第1項が操作の最中に実際に失われた分である。式 (4.3.17) は式 (4.3.16) と同じ内容である。

例題 4.3.3 序章の緩和を経路で読み直す ★★

オリジナル例題 4.3.1 の2準位で βε=ln⁡3\beta\varepsilon =\ln3 に固定し(λ\lambda を動かさない)、系を確実に e に置いた分布 p0=(0,1)p_0=(0,1) から1回だけ跳ばせる。

  1. (1) p1p_1 と ⟨Q⟩\langle Q \rangle を求め、式 (4.3.15) で ⟨Σ⟩\langle \Sigma \rangle を求めよ。
  2. (2) 経路 (x0,x1)(x_0,x_1) を全部書き、D(P ∥ P~)D(P\,\|\,\tilde P)(p~0=p1\tilde p_0=p_1)を直接計算して (1) と一致することを確かめよ。
  3. (3) (1) の値が D(p0 ∥ peq)−D(p1 ∥ peq)D(p_0\,\|\,p_{\mathrm{eq}})-D(p_1\,\|\,p_{\mathrm{eq}}) に等しいことを示せ。
解答

(1) e から g へ 12\frac12、e にとどまる 12\frac12 なので p1=(12,12)p_1=(\frac12,\frac12)。系がもらう熱は g へ落ちる経路だけで −ε-\varepsilon、⟨Q⟩=−12ε\langle Q \rangle =-\frac12\varepsilon、−β⟨Q⟩=12ln⁡3-\beta\langle Q \rangle =\frac12\ln3。H(p0)=0H(p_0)=0、H(p1)=ln⁡2H(p_1)=\ln2 なので ⟨Σ⟩=ln⁡2+12ln⁡3≈1.2425\langle \Sigma \rangle =\ln2+\frac12\ln3\approx1.2425 nat ≈1.79\approx1.79 ビット。

(2) e から出る経路は2本。(e,g)(\mathrm e,\mathrm g):P=12P=\frac12、P~=p1(g) K(e ∣ g)=12⋅16=112\tilde P=p_1(\mathrm g)\,K(\mathrm e\,|\,\mathrm g)=\frac12\cdot\frac16=\frac1{12}、ln⁡PP~=ln⁡6\ln\frac{P}{\tilde P}=\ln6。(e,e)(\mathrm e,\mathrm e):P=12P=\frac12、P~=12⋅12=14\tilde P=\frac12\cdot\frac12=\frac14、ln⁡2\ln2。平均は 12ln⁡6+12ln⁡2=12ln⁡12=ln⁡2+12ln⁡3\frac12\ln6+\frac12\ln2=\frac12\ln12=\ln2+\frac12\ln3 で一致する。

(3) peq=(34,14)p_{\mathrm{eq}}=(\frac34,\frac14) で D(p0 ∥ peq)=ln⁡4D(p_0\,\|\,p_{\mathrm{eq}})=\ln4、D(p1 ∥ peq)=12ln⁡43D(p_1\,\|\,p_{\mathrm{eq}})=\frac12\ln\frac43。差は ln⁡4−12ln⁡4+12ln⁡3=ln⁡2+12ln⁡3\ln4-\frac12\ln4+\frac12\ln3=\ln2+\frac12\ln3。一般に λ\lambda を固定すると W=0W=0、ΔF=0\Delta F=0 なので式 (4.3.16) から kBT⟨Σ⟩=F[p0]−F[p1]k_{\mathrm B} T\langle \Sigma \rangle =F[p_0]-F[p_1]、これに F[p]−Feq=kBT D(p ∥ peq)F[p]-F_{\mathrm{eq}}= k_{\mathrm B} T\,D(p\,\|\,p_{\mathrm{eq}}) を入れたものである。第1部第1章 の総合問題 (5) で「失われた DD =生じたエントロピー」と読んだ 1.2421.242 が、経路の確率の DD として正確に得られた。

1 C. Jarzynski (1997), Phys. Rev. Lett. 78, 2690–2693。

2 G. E. Crooks (1999), Phys. Rev. E 60, 2721–2726。

木下の勉強ステータス

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

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

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