tecktom.com

条件付き極限定理:Gibbs 分布が出てくる

2026年10月15日8分で読めますAI生成教科書『虎の巻』第1部 第8章 大偏差・Sanov・Stein の補題
目次
  1. 条件をつけた世界の経験分布
  2. Gibbs 分布が出る理由
  3. 木下の層

AI 生成・木下が未勉強

相対エントロピーから情報と物理を組み立てる 第1部 第8章 §1.8.4 ★★★ めったに起きないことが起きたと知ると、1個の標本は I 射影の分布に従う。微視的カノニカルからカノニカルへの橋

条件をつけた世界の経験分布

Sanov の定理は、めったに起きないことが「最もありそうな起こり方」1つでほぼ起きると言った。これを確率の言葉に直したものが次の定理である。

定理 1.8.7(条件付き極限定理 ★★★)

EE を凸な閉集合で内部の閉包に等しいものとし、q∉Eq\notin E、p∗p^* を qq の EE への I 射影とする。X1,…,XnX_1,\dots,X_n が独立に qq に従うとき、条件 pXn∈Ep_{X^n}\in E のもとで、

  1. 任意の δ>0\delta>0 について qn(D(pXn∥p∗)>δ ∣ pXn∈E)→0q^n\bigl(D(p_{X^n}\|p^*)>\delta\ \big|\ p_{X^n}\in E\bigr)\to0(n→∞n\to\infty)。
  2. 1個の標本の条件付き分布は p∗p^* に近づく:qn(X1=x ∣ pXn∈E)→p∗(x)q^n\bigl(X_1=x\ \big|\ p_{X^n}\in E\bigr)\to p^*(x)。
証明

(1) 分子と分母を別々に評価する。分子は、EE の型のうち D(p∥p∗)>δD(p\|p^*)>\delta のものの確率である。そういう型について Pythagoras の不等式 (1.8.8) からD(p∥q)≥D(p∥p∗)+D(p∗∥q)>D(p∗∥q)+δD(p\|q)\ge D(p\|p^*)+D(p^*\|q)>D(p^*\|q)+\delta なので、Sanov の上界 (1.8.6) と同じ数え方で、分子は (n+1)∣X∣e−n(D(p∗∥q)+δ)(n+1)^{\lvert \mathcal X\rvert }\mathrm{e}^{-n(D(p^*\|q)+\delta)} 以下になる。分母は Sanov の定理 (1.8.7) から、nn が十分大きければ e−n(D(p∗∥q)+δ/2)\mathrm{e}^{-n(D(p^*\|q)+\delta/2)} 以上になる。比は (n+1)∣X∣e−nδ/2→0(n+1)^{\lvert \mathcal X\rvert }\mathrm{e}^{-n\delta/2}\to0。

(2) 条件 pXn∈Ep_{X^n}\in E は系列の並べ替えで変わらないので、条件のもとで X1,…,XnX_1,\dots,X_n はどれも同じ分布に従う。したがってqn(X1=x ∣ E)=1n∑iqn(Xi=x ∣ E)=⟨pXn(x) ∣ E⟩q^n(X_1=x\,|\,E)=\frac1n\sum_iq^n(X_i=x\,|\,E)=\langle p_{X^n}(x)\,|\,E \rangle で、条件付きの経験分布の平均に等しい。(1) と Pinsker の不等式 ∑x∣p(x)−p∗(x)∣≤2D(p∥p∗)\sum_x\lvert p(x)-p^*(x)\rvert \le\sqrt{2D(p\|p^*)} から、経験分布は条件のもとで p∗p^* の近くに集まるので、その平均も p∗(x)p^*(x) に近づく。

第2部第1章(定理 2.1.2)は、この定理を最大エントロピー原理の側から使い、有限の個数のずれまで評価する。この定理は、有限の X\mathcal X では型の数え上げだけで示すことができる。一般の可測空間で、I 射影が存在しない場合も含めて示したのは I. Csiszár である1I. Csiszár (1984), Ann. Probab. 12, 768–793。下の脚注へ。定理の (2) は1個の標本についての主張だが、全エネルギーをちょうど決める条件(EE が超平面のとき)でも、超平面の上の型を数えれば同じ結論が得られる。図 1.8.3 はこの場合の数値である。同じ論法で、固定した kk 個の標本 X1,…,XkX_1,\dots,X_k が条件のもとで漸近的に独立に p∗p^* に従うことも示すことができる。

Gibbs 分布が出る理由

EE を平均を決めた集合 Ea={p:⟨f⟩p≥a}E_a=\{p:\langle f \rangle _p\ge a\}(a>μa>\mu)に選ぶと、1.8.3節で見たとおり I 射影は指数型の分布 qλ∗q_{\lambda^*} である。ここで物理の言葉に読み替える。X\mathcal X を1個の粒子の微視的な状態の集合、f(x)=ε(x)f(x)=\varepsilon (x) をその状態のエネルギー、qq を一様分布とする。系列 XnX^n は nn 個の粒子の微視的な状態で、qnq^n はすべての微視的な状態を等しく重みづける。全エネルギーを ∑iε(Xi)≈nu\sum_i\varepsilon (X_i)\approx nu に限ることは、微視的カノニカル集団をとることにあたる。条件付き極限定理によれば、このとき1個の粒子の状態は

p∗(x)=e−βε(x)Z(β),Z(β)=∑xe−βε(x),⟨ε⟩p∗=u(1.8.15)p^*(x)=\frac{\mathrm{e}^{-\beta\varepsilon (x)}}{Z(\beta)},\qquad Z(\beta)=\sum_x\mathrm{e}^{-\beta\varepsilon (x)},\qquad\langle \varepsilon \rangle _{p^*}=u\tag{1.8.15}

に従う。β=−λ∗\beta=-\lambda^* で、uu が一様分布の平均エネルギーより小さいときは λ∗<0\lambda^*<0、すなわち β>0\beta>0 である(1.8.3節は a>μa>\mu で書いたが、a<μa<\mu でも符号を変えて同じ議論ができる)。「大きな系全体のエネルギーを決めると、その一部分はカノニカル分布に従う」という統計力学の標準的な議論は、普通は熱浴の状態数を展開して示す。ここではそれが「エネルギーを決めた条件のもとで最もありそうな経験分布は、エネルギーの平均を決めて D(p∥一様)D(p\|\text{一様}) を最小にする分布、すなわちエントロピー H(p)H(p) を最大にする分布である」と言い換えられた。D(p∥一様)=ln⁡∣X∣−H(p)D(p\|\text{一様})=\ln\lvert \mathcal X\rvert -H(p) なので、DD の最小とエントロピーの最大は同じことである。この「最大エントロピー」の側からの読み方は 第2部第1章 で推論の原理として取り上げ直す。

近似の裏付け

型 B(極限) | 小さいもの 1/n→01/n\to0(粒子数 n→∞n\to\infty。微視的カノニカルとカノニカルの等価性) | 証明 証明済み(有限の状態空間・独立な粒子で、条件の集合が凸で内部の閉包に等しいとき。Van Campenhout–Cover 1981・Csiszár 1984) | 行き先 近似の数理 第6章

図

図 1.8.3 エネルギー 0,1,20,1,2(単位 ε0\varepsilon _0)の3準位の粒子 nn 個で、全エネルギーを n/2n/2 に決めたときに粒子1個が基底状態にある確率。丸印は状態の数え上げによる厳密な値、点線は例題 1.8.5 の Gibbs 分布の値。

例題 1.8.5 3準位の粒子の集まり ★★

オリジナル1個の粒子のエネルギーが 0,ε0,2ε00,\varepsilon _0,2\varepsilon _0 の3つの状態をとる。nn 個の粒子の微視的な状態を等しく重みづけ、全エネルギーを nε0/2n\varepsilon _0/2 に決める。

  1. 条件付き極限定理から、粒子1個の状態の分布 p∗p^* を求めよ。z=e−βε0z=\mathrm{e}^{-\beta\varepsilon _0} とおいて zz と βε0\beta\varepsilon _0 を数値で出せ。
  2. nn 個の粒子の状態でエネルギーの和が EE(単位 ε0\varepsilon _0)になるものの数は、(1+t+t2)n(1+t+t^2)^n の tEt^E の係数 cn(E)c_n(E) である。粒子1個が基底状態にある厳密な条件付き確率を cnc_n で表せ。n=200n=200 で 0.61580.6158 になることを使い、(1) と比べよ。
解答

(1) p∗∝(1,z,z2)p^*\propto(1,z,z^2) で、平均エネルギーの条件は z+2z21+z+z2=12\frac{z+2z^2}{1+z+z^2}=\frac12、すなわち 3z2+z−1=03z^2+z-1=0。正の根は z=−1+136=0.434z=\frac{-1+\sqrt{13}}{6}=0.434 で、βε0=−ln⁡z=0.834\beta\varepsilon _0=-\ln z=0.834。分布は p∗=(1,z,z2)1+z+z2=(0.616, 0.268, 0.116)p^*=\frac{(1,z,z^2)}{1+z+z^2}=(0.616,\,0.268,\,0.116)。

(2) 粒子1が基底状態なら、残り n−1n-1 個でエネルギー n/2n/2 を分け合うので、その数は cn−1(n/2)c_{n-1}(n/2)。したがってqn(X1=0 ∣ ⋅)=cn−1(n/2)/cn(n/2)q^n(X_1=0\,|\,\cdot)=c_{n-1}(n/2)/c_n(n/2)。n=200n=200 の値 0.61580.6158 は (1) の 0.61620.6162 と 4×10−44\times10^{-4} だけ違い、図 1.8.3 のとおり nn とともに近づく。全エネルギーを決めた微視的カノニカルな数え上げと、温度 β\beta を決めたカノニカル分布が、粒子数が多いと1粒子の分布について同じ答えになることが数値で確かめられた。

ほかの巻では

1 I. Csiszár (1984), Ann. Probab. 12, 768–793。

木下の勉強ステータス

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

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

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