相対エントロピーから情報と物理を組み立てる 第2部 第2章 §2.2.2 ★★★ Gibbs 分布の族に対して相対エントロピーは直角三角形のように分かれ、自由エネルギーの余りが「同じエネルギーでの最大エントロピーからの不足」と「温度の違い」の和になる
凸性
序章では非負性とデータ処理不等式を使った。この章ではもう1つ、凸性を使う。
命題 2.2.2(F[p] の凸性 ★★)
0<λ<1 と2つの分布 p1=p2 について
F[λp1+(1−λ)p2]<λF[p1]+(1−λ)F[p2].
証明
平均エネルギー ⟨E⟩p は p について1次なので、左右で等しい。残りは −H(p) が狭義に凸であることで、関数 xlnx が x>0 で2階導関数 1/x>0 を持つことから、各成分ごとに狭義の凸不等式が成り立つ。p1=p2 なら少なくとも1つの成分で狭義の不等号になる。
凸関数の極小は最小で、狭義の凸関数の最小を与える点はただ1つである。したがって F[p] を下げる方向に分布を少しずつ動かす計算は、途中で別の極小に止まることがない。2.2.4節で試す分布の族を狭めると、この性質は失われる。族が凸集合でないと、F を族の上に制限した関数に極小が複数現れうるからである。命題を言い換えると、2つの分布を混ぜると自由エネルギーは平均より下がる。下がる分は kBT と混合のエントロピー H(λp1+(1−λ)p2)−λH(p1)−(1−λ)H(p2) の積で、混合のエントロピーは λ=1/2 のとき ln2 を超えない。
Pythagoras の分解
温度 T の Gibbs 分布を pβ と書き、β を動かした族 {pβ} を考える。任意の分布 p の平均エネルギーを Up=⟨E⟩p とし、同じ平均エネルギーを持つ Gibbs 分布を pβ∗ とする。pβ の平均エネルギー U(β)=−∂βlnZ は β について狭義に減少する(2.2.5節で示すように ∂βU=−Varβ(E)<0)ので、Up が最小と最大のエネルギーの間にあれば β∗ はただ1つ決まる。β∗ は負でもよい。上の準位ほど多く占める分布は、負の温度の Gibbs 分布と同じ平均エネルギーを持つ。
定理 2.2.3(Pythagoras の分解 ★★★)
⟨E⟩p=⟨E⟩pβ∗ なら、すべての β について
D(p∥pβ)=D(p∥pβ∗)+D(pβ∗∥pβ),(2.2.7)
さらに
D(p∥pβ∗)=H(pβ∗)−H(p)(2.2.8)
が成り立つ。
証明
lnpβ(i)pβ∗(i)=−(β∗−β)Ei−lnZ(β∗)+lnZ(β) は Ei の1次式なので、p で平均しても pβ∗ で平均しても同じ値になる。よって
D(p∥pβ)−D(p∥pβ∗)=i∑p(i)lnpβ(i)pβ∗(i)=i∑pβ∗(i)lnpβ(i)pβ∗(i)=D(pβ∗∥pβ).
式 (2.2.8) も同じ理由で、lnpβ∗(i)=−β∗Ei−lnZ(β∗) が Ei の1次式なので −∑ip(i)lnpβ∗(i)=−∑ipβ∗(i)lnpβ∗(i)=H(pβ∗) となり、D(p∥pβ∗)=−H(p)+H(pβ∗) を得る。
式 (2.2.8) と非負性から H(p)≤H(pβ∗) がすべての p について成り立つ。これは平均エネルギーを決めたときエントロピーが Gibbs 分布で最大になる、という前の章の最大エントロピー原理そのものである。式 (2.2.7) に kBT を掛け、式 (2.2.3) と組み合わせると
F[p]−Feq=同じ平均エネルギーでのエントロピーの不足T[kBH(pβ∗)−kBH(p)]+温度の違いkBTD(pβ∗∥pβ)(2.2.9)
となる。第1項は、平均エネルギーを変えずにエントロピーを最大にすれば消える部分で、前の章の最大エントロピー原理が扱った部分である。第2項は、温度 T∗=1/(kBβ∗) の平衡状態が温度 T の熱浴と接したときの余りで、熱のやりとりで消える部分である。最大エントロピー原理と自由エネルギー最小の原理は、1つの恒等式 (2.2.7) の2つの項として並ぶ。

図 2.2.1 Pythagoras の分解 (2.2.7) の模式図。平均エネルギーが Up に等しい分布の集合(縦の線)と Gibbs 分布の族(曲線)は pβ∗ で交わり、相対エントロピーについて「直交」している。距離ではなく相対エントロピーの和として成り立つ関係で、図の長さは目安である。
この「直交」は Euclid 空間の直交とは意味が違う。相対エントロピーは距離ではなく、2乗した距離に近い量(D(pθ∥pθ+δ)≈21J(θ)δ2)なので、2乗の距離が足し算になる関係が直角三角形にあたる。式 (2.2.7) が近似でなく等式で成り立つのは、Gibbs 分布の対数がエネルギーの1次式であることによる。この性質を持つ分布の族を指数型分布族と呼び、統計学の最尤推定でも同じ分解が成り立つ(第1部第4章)。
例題 2.2.2 3準位系で分解を確かめる ★★
オリジナルエネルギー 0,ε,2ε の3準位系が、βε=ln2 の熱浴と接している。系の分布が p=(1/2,0,1/2) のとき、次を求めよ。
- (1) 熱浴の Gibbs 分布 pβ と、p と同じ平均エネルギーを持つ Gibbs 分布 pβ∗。
- (2) 式 (2.2.7) の3つの相対エントロピーを求め、等式を確かめよ。F[p]−Feq のうち温度の違いの分は何割か。
解答
(1) e−βε=1/2 なので pβ∝(1,21,41)、すなわち pβ=(4,2,1)/7。p の平均エネルギーは 21⋅0+21⋅2ε=ε である。x=e−β∗ε とおくと Gibbs 分布の平均エネルギーは ε(x+2x2)/(1+x+x2) で、これを ε とおくと x+2x2=1+x+x2、すなわち x2=1。x>0 から x=1、β∗=0 で、pβ∗=(1/3,1/3,1/3) は一様分布である。無限に高い温度の Gibbs 分布が、p と同じ平均エネルギーを持つ。
(2) D(p∥pβ∗)=ln3−H(p)=ln3−ln2=ln23≈0.405。D(pβ∗∥pβ)=31(ln4/71/3+ln2/71/3+ln1/71/3)=ln37−31ln8=ln67≈0.154。D(p∥pβ)=21ln4/71/2+21ln1/71/2=21ln1649=ln47≈0.560。ln23+ln67=ln1221=ln47 で等式が成り立つ。温度の違いの分は 0.154/0.560≈0.28、すなわち約3割で、残りの約7割は同じ平均エネルギーでのエントロピーの不足である。p は真ん中の準位を空にしているので、平均エネルギーを保ったまま真ん中の準位へ確率を移すだけで、余りの約7割は消える。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com