相対エントロピーから情報と物理を組み立てる 第1部 第8章 §1.8.4 ★★★ めったに起きないことが起きたと知ると、1個の標本は I 射影の分布に従う。微視的カノニカルからカノニカルへの橋
条件をつけた世界の経験分布
Sanov の定理は、めったに起きないことが「最もありそうな起こり方」1つでほぼ起きると言った。これを確率の言葉に直したものが次の定理である。
定理 1.8.7(条件付き極限定理 ★★★)
E を凸な閉集合で内部の閉包に等しいものとし、q∈/E、p∗ を q の E への I 射影とする。X1,…,Xn が独立に q に従うとき、条件 pXn∈E のもとで、
- 任意の δ>0 について qn(D(pXn∥p∗)>δ pXn∈E)→0(n→∞)。
- 1個の標本の条件付き分布は p∗ に近づく:qn(X1=x pXn∈E)→p∗(x)。
証明
(1) 分子と分母を別々に評価する。分子は、E の型のうち D(p∥p∗)>δ のものの確率である。そういう型について Pythagoras の不等式 (1.8.8) からD(p∥q)≥D(p∥p∗)+D(p∗∥q)>D(p∗∥q)+δ なので、Sanov の上界 (1.8.6) と同じ数え方で、分子は (n+1)∣X∣e−n(D(p∗∥q)+δ) 以下になる。分母は Sanov の定理 (1.8.7) から、n が十分大きければ e−n(D(p∗∥q)+δ/2) 以上になる。比は (n+1)∣X∣e−nδ/2→0。
(2) 条件 pXn∈E は系列の並べ替えで変わらないので、条件のもとで X1,…,Xn はどれも同じ分布に従う。したがってqn(X1=x∣E)=n1∑iqn(Xi=x∣E)=⟨pXn(x)∣E⟩ で、条件付きの経験分布の平均に等しい。(1) と Pinsker の不等式 ∑x∣p(x)−p∗(x)∣≤2D(p∥p∗) から、経験分布は条件のもとで p∗ の近くに集まるので、その平均も p∗(x) に近づく。
第2部第1章(定理 2.1.2)は、この定理を最大エントロピー原理の側から使い、有限の個数のずれまで評価する。この定理は、有限の X では型の数え上げだけで示すことができる。一般の可測空間で、I 射影が存在しない場合も含めて示したのは I. Csiszár である1I. Csiszár (1984), Ann. Probab. 12, 768–793。。定理の (2) は1個の標本についての主張だが、全エネルギーをちょうど決める条件(E が超平面のとき)でも、超平面の上の型を数えれば同じ結論が得られる。図 1.8.3 はこの場合の数値である。同じ論法で、固定した k 個の標本 X1,…,Xk が条件のもとで漸近的に独立に p∗ に従うことも示すことができる。
Gibbs 分布が出る理由
E を平均を決めた集合 Ea={p:⟨f⟩p≥a}(a>μ)に選ぶと、1.8.3節で見たとおり I 射影は指数型の分布 qλ∗ である。ここで物理の言葉に読み替える。X を1個の粒子の微視的な状態の集合、f(x)=ε(x) をその状態のエネルギー、q を一様分布とする。系列 Xn は n 個の粒子の微視的な状態で、qn はすべての微視的な状態を等しく重みづける。全エネルギーを ∑iε(Xi)≈nu に限ることは、微視的カノニカル集団をとることにあたる。条件付き極限定理によれば、このとき1個の粒子の状態は
p∗(x)=Z(β)e−βε(x),Z(β)=x∑e−βε(x),⟨ε⟩p∗=u(1.8.15)
に従う。β=−λ∗ で、u が一様分布の平均エネルギーより小さいときは λ∗<0、すなわち β>0 である(1.8.3節は a>μ で書いたが、a<μ でも符号を変えて同じ議論ができる)。「大きな系全体のエネルギーを決めると、その一部分はカノニカル分布に従う」という統計力学の標準的な議論は、普通は熱浴の状態数を展開して示す。ここではそれが「エネルギーを決めた条件のもとで最もありそうな経験分布は、エネルギーの平均を決めて D(p∥一様) を最小にする分布、すなわちエントロピー H(p) を最大にする分布である」と言い換えられた。D(p∥一様)=ln∣X∣−H(p) なので、D の最小とエントロピーの最大は同じことである。この「最大エントロピー」の側からの読み方は 第2部第1章 で推論の原理として取り上げ直す。
近似の裏付け
型 B(極限) | 小さいもの 1/n→0(粒子数 n→∞。微視的カノニカルとカノニカルの等価性) | 証明 証明済み(有限の状態空間・独立な粒子で、条件の集合が凸で内部の閉包に等しいとき。Van Campenhout–Cover 1981・Csiszár 1984) | 行き先 近似の数理 第6章

図 1.8.3 エネルギー 0,1,2(単位 ε0)の3準位の粒子 n 個で、全エネルギーを n/2 に決めたときに粒子1個が基底状態にある確率。丸印は状態の数え上げによる厳密な値、点線は例題 1.8.5 の Gibbs 分布の値。
例題 1.8.5 3準位の粒子の集まり ★★
オリジナル1個の粒子のエネルギーが 0,ε0,2ε0 の3つの状態をとる。n 個の粒子の微視的な状態を等しく重みづけ、全エネルギーを nε0/2 に決める。
- 条件付き極限定理から、粒子1個の状態の分布 p∗ を求めよ。z=e−βε0 とおいて z と βε0 を数値で出せ。
- n 個の粒子の状態でエネルギーの和が E(単位 ε0)になるものの数は、(1+t+t2)n の tE の係数 cn(E) である。粒子1個が基底状態にある厳密な条件付き確率を cn で表せ。n=200 で 0.6158 になることを使い、(1) と比べよ。
解答
(1) p∗∝(1,z,z2) で、平均エネルギーの条件は 1+z+z2z+2z2=21、すなわち 3z2+z−1=0。正の根は z=6−1+13=0.434 で、βε0=−lnz=0.834。分布は p∗=1+z+z2(1,z,z2)=(0.616,0.268,0.116)。
(2) 粒子1が基底状態なら、残り n−1 個でエネルギー n/2 を分け合うので、その数は cn−1(n/2)。したがってqn(X1=0∣⋅)=cn−1(n/2)/cn(n/2)。n=200 の値 0.6158 は (1) の 0.6162 と 4×10−4 だけ違い、図 1.8.3 のとおり n とともに近づく。全エネルギーを決めた微視的カノニカルな数え上げと、温度 β を決めたカノニカル分布が、粒子数が多いと1粒子の分布について同じ答えになることが数値で確かめられた。
1 I. Csiszár (1984), Ann. Probab. 12, 768–793。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com