コンテンツにスキップ

DLR 条件と Gibbs 測度:仕様・存在・極値測度と Dobrushin の一意性定理

対象読者: 本テーマ第 1 章(圧力と熱力学的極限)と、条件付き期待値・確率測度の弱収束(Prokhorov の定理)を学んだ修士課程の人。Gibbs 測度の知識は仮定しません。

前提となる記事:熱力学的極限:圧力の存在・凸性・Gibbs の変分原理と Cramér の定理、条件付き期待値:σ-加法族で条件づけ、ラドン・ニコディムで存在を保証する、確率測度の弱収束:Portmanteau・Prokhorov・Skorokhod の表現定理

生 Markdown

磁石の中には、ほぼ無限に多い小さな磁石(スピン)が並んでいます。有限個ならエネルギーから確率を決める公式がありますが、無限個ではその公式がそのままでは使えません。この記事では、「ある範囲の外を決めたとき、範囲の中はどんな確率で決まるか」という局所的な規則だけから、無限個の系の確率を定義します。読み終えると、その確率がいつ存在するか、いつただ一つに決まるか(決まらないときが相転移です)を、証明つきで説明できるようになります。

  • 無限個のスピンの同時分布は、Boltzmann 因子 e−H/Ze^{-H}/Z では直接書けません(HH も ZZ も発散する)。代わりに「有限領域 Λ\Lambda の外を固定したときの Λ\Lambda 内の条件付き分布」を指定します。この指定の整合的な族が仕様 γ=(γΛ)\gamma=(\gamma_\Lambda) で、それを条件付き分布にもつ確率測度が Gibbs 測度です(DLR 条件 μγΛ=μ\mu\gamma_\Lambda=\mu)。
  • 存在(定理 5.1): 仕様が連続関数を連続関数に写す(準局所的)なら、有限体積の分布の弱収束極限は Gibbs 測度です。スピンが有限集合に値をとれば配位空間がコンパクトで緊密性は自動的に成り立ち、Gibbs 測度は必ず存在します。
  • 凸構造(定理 6.2): Gibbs 測度全体 G(γ)\mathcal{G}(\gamma) は凸でコンパクトです。その端点(極値 Gibbs 測度。混ぜ合わせで作れない測度のこと)は、ちょうど尾部 σ-加法族(有限個のスピンを変えても変わらない事象の全体)の上で自明(確率が 00 か 11)なものです。異なる極値 Gibbs 測度は互いに特異(確率 1 を持つ事象が互いに交わらない)で、物理的には別々の「純粋な相」です。
  • Dobrushin の一意性定理(定理 7.3): 1 点の条件付き分布が、遠くのスピンを 1 個変えたときに全変動距離でどれだけ動くかを並べた行列 CC の行和が 11 未満なら、Gibbs 測度はただ一つです。Ising 模型では 2dtanh⁡(β∣J∣)<12d\tanh(\beta|J|)<1(外場は任意)で一意性が出ます。
  • 相転移は「∣G(γ)∣≥2|\mathcal{G}(\gamma)|\ge2」として定義でき、次章 の Peierls の議論はこれを低温で示します。

前章 では、圧力 P(Φ)=lim⁡∣Λ∣−1log⁡ZΛP(\Phi)=\lim|\Lambda|^{-1}\log Z_\Lambda が存在し、境界条件によらないことを示しました。圧力は熱力学量ですが、相転移の本質は状態、すなわち無限系の確率測度の側にあります。たとえば低温の強磁性体では、境界をすべて ++ にして体積を大きくしたときと、すべて −- にしたときで、原点のスピンの期待値が異なる値に収束します。境界の影響が無限遠まで残るのです。

この現象を記述するには、無限体積の確率測度そのものを扱う必要があります。Dobrushin(1968)と Lanford–Ruelle(1969)は、Boltzmann 分布の代わりに条件付き分布を指定するという考えでこれを定式化しました(DLR 条件)。これは確率過程を有限次元分布で指定する Kolmogorov の方法(Kolmogorov の拡張定理)と似ていますが、決定的な違いがあります。有限次元分布の整合的な族は確率測度を一意に決めますが、条件付き分布の整合的な族は確率測度を一意に決めるとは限りません。この非一意性が相転移です。

前章と同じく SS は有限集合で離散位相を入れ、Ω=SZd\Omega=S^{\mathbb{Z}^d} に直積位相を入れます。Λ⊂Zd\Lambda\subset\mathbb{Z}^d について FΛ\mathcal{F}_\Lambda を座標 ωx\omega_x(x∈Λx\in\Lambda)が生成する σ-加法族、F=FZd\mathcal{F}=\mathcal{F}_{\mathbb{Z}^d}、TΛ=FZd∖Λ\mathcal{T}_\Lambda=\mathcal{F}_{\mathbb{Z}^d\setminus\Lambda} とし、

T=⋂Λ⋐ZdTΛ\mathcal{T}=\bigcap_{\Lambda\Subset\mathbb{Z}^d}\mathcal{T}_\Lambda

を尾部 σ-加法族と呼びます。尾部事象は「有限個のスピンを変えても変わらない事象」です。たとえば {lim⁡n∣Λn∣−1∑x∈Λnωx\{\lim_n|\Lambda_n|^{-1}\sum_{x\in\Lambda_n}\omega_x が存在して正}\}(S⊂RS\subset\mathbb{R})は尾部事象です。

FΛ\mathcal{F}_\Lambda-可測な関数(Λ\Lambda は有限)を局所関数と呼びます。

補題 2.1(配位空間の性質)✓

(i) ∥x∥=max⁡i∣xi∣\|x\|=\max_i|x_i| とし、ω≠η\omega\ne\eta に対して d(ω,η)=2−kd(\omega,\eta)=2^{-k}、k=min⁡{∥x∥:ωx≠ηx}k=\min\{\|x\|:\omega_x\ne\eta_x\} とおくと、dd は直積位相を定める距離で、(Ω,d)(\Omega,d) はコンパクトである。

(ii) F\mathcal{F} は Ω\Omega のボレル σ-加法族に等しい。

(iii) 局所関数は連続で、C(Ω)C(\Omega) の中で一様ノルムについて稠密である。

証明(補題 2.1)

(i) d(ω,η)<2−kd(\omega,\eta)<2^{-k} は「∥x∥≤k\|x\|\le k で ωx=ηx\omega_x=\eta_x」と同値なので、dd の開球は箱 Bk={∥x∥≤k}B_k=\{\|x\|\le k\} 上の値を固定した集合(柱集合)で、これは直積位相の基本開集合です。逆に有限個の座標を固定した集合は、それを含む BkB_k の上の値を固定した柱集合の有限和です。三角不等式は強い形 d(ω,ζ)≤max⁡(d(ω,η),d(η,ζ))d(\omega,\zeta)\le\max(d(\omega,\eta),d(\eta,\zeta)) で成り立ちます。コンパクト性は点列で示します。(ωn)(\omega^n) に対し、SS は有限なので ωB1n\omega^n_{B_1} が一定になる部分列をとり、その中で ωB2n\omega^n_{B_2} が一定になる部分列をとり……と続け、対角線の部分列をとれば、それは各 BkB_k 上で最終的に一定で、収束します。

(ii) 柱集合は可算個しかなく、dd の開集合は開球(柱集合)の和集合なので、可算個の柱集合の和です。よって開集合は F\mathcal{F} に属します。逆に座標写像 ω↦ωx\omega\mapsto\omega_x は連続なので F\mathcal{F} はボレル集合族に含まれます。

(iii) 局所関数は有限個の座標の関数で、座標は連続なので連続です。f∈C(Ω)f\in C(\Omega) と ε>0\varepsilon>0 について、一様連続性(補題 2.2[Arzelà–Ascoli と Stone–Weierstrass])から kk を、d(ω,η)<2−kd(\omega,\eta)<2^{-k} なら ∣f(ω)−f(η)∣<ε|f(\omega)-f(\eta)|<\varepsilon となるようにとります。η0∈Ω\eta^0\in\Omega を固定して fk(ω)=f(ωBkηBkc0)f_k(\omega)=f(\omega_{B_k}\eta^0_{B_k^c}) とおくと、fkf_k は局所関数で ∥f−fk∥∞≤ε\|f-f_k\|_\infty\le\varepsilon です。

確率測度 μ\mu は (Ω,F)(\Omega,\mathcal{F}) 上で考え、μ(f)=∫f dμ\mu(f)=\int f\,d\mu と書きます。局所関数(あるいは連続関数)の積分が一致する 2 つの確率測度は等しいので(補題 2.2[確率測度の弱収束] と (iii))、確率測度は局所関数の期待値で決まります。

この節のねらいは、「外を固定したときの内側の確率」を一般的な形で扱うための言葉をそろえることです。たとえば Ising 模型で 1 点 xx だけを再標本化するなら、隣のスピンの値(外)を与えると、xx のスピンが +1+1 になる確率は式 (1+tanh⁡βm)/2(1+\tanh\beta m)/2(mm は隣のスピンと外場から決まる局所場)で決まります(第 7 節で導きます)。この「外の配位から内側の確率分布を返す対応」が確率核で、以下の定義はそれを一般化したものです。

Ω\Omega 上の確率核 π\pi を G⊂F\mathcal{G}\subset\mathcal{F} から F\mathcal{F} への写像 (B,ω)↦π(B∣ω)(B,\omega)\mapsto\pi(B\mid\omega) で、各 ω\omega について π(⋅∣ω)\pi(\cdot\mid\omega) が確率測度、各 B∈FB\in\mathcal{F} について ω↦π(B∣ω)\omega\mapsto\pi(B\mid\omega) が G\mathcal{G}-可測なものとします。有界可測な ff に πf(ω)=∫f(η) π(dη∣ω)\pi f(\omega)=\int f(\eta)\,\pi(d\eta\mid\omega)、確率測度 μ\mu に μπ(B)=∫π(B∣ω) μ(dω)\mu\pi(B)=\int\pi(B\mid\omega)\,\mu(d\omega) を対応させます。μπ(f)=μ(πf)\mu\pi(f)=\mu(\pi f) は、ff が指示関数なら定義そのもので、単関数を経て有界可測関数に広がります。核の合成は (ππ′)(B∣ω)=∫π′(B∣η) π(dη∣ω)(\pi\pi')(B\mid\omega)=\int\pi'(B\mid\eta)\,\pi(d\eta\mid\omega) です。

定義 3.1(仕様)↗

Λ⋐Zd\Lambda\Subset\mathbb{Z}^d ごとに TΛ\mathcal{T}_\Lambda から F\mathcal{F} への確率核 γΛ\gamma_\Lambda が与えられ、次を満たすとき、族 γ=(γΛ)Λ⋐Zd\gamma=(\gamma_\Lambda)_{\Lambda\Subset\mathbb{Z}^d} を仕様という。

  1. (固有性)B∈TΛB\in\mathcal{T}_\Lambda ならば γΛ(B∣ω)=1B(ω)\gamma_\Lambda(B\mid\omega)=\mathbf{1}_B(\omega)。
  2. (整合性)Λ⊂Δ\Lambda\subset\Delta ならば γΔγΛ=γΔ\gamma_\Delta\gamma_\Lambda=\gamma_\Delta。

さらに、すべての Λ\Lambda と f∈C(Ω)f\in C(\Omega) について γΛf∈C(Ω)\gamma_\Lambda f\in C(\Omega) のとき、γ\gamma は準局所的であるという。

固有性は「Λ\Lambda の外は固定されたまま」を、整合性は「大きな領域で再標本化した後に小さな領域で再標本化しても分布は変わらない」を意味します。

補題 3.2(固有性からの取り出し)✓

仕様 γ\gamma と有界 TΛ\mathcal{T}_\Lambda-可測な gg、有界可測な ff について γΛ(gf)=g γΛf\gamma_\Lambda(gf)=g\,\gamma_\Lambda f。

証明(補題 3.2)

g=1Bg=\mathbf{1}_B(B∈TΛB\in\mathcal{T}_\Lambda)のとき、ω∈B\omega\in B なら γΛ(Bc∣ω)=0\gamma_\Lambda(B^c\mid\omega)=0 なので γΛ(1Bf)(ω)=γΛf(ω)\gamma_\Lambda(\mathbf{1}_Bf)(\omega)=\gamma_\Lambda f(\omega)、ω∉B\omega\notin B なら γΛ(B∣ω)=0\gamma_\Lambda(B\mid\omega)=0 なので γΛ(1Bf)(ω)=0\gamma_\Lambda(\mathbf{1}_Bf)(\omega)=0 です。線形性で単関数に、一様近似(∣γΛh∣≤∥h∥∞|\gamma_\Lambda h|\le\|h\|_\infty)で有界可測関数に広がります。

物理で使うのは相互作用から作る仕様です。前章と違い、ここでは並進不変性を仮定しません。

定義 3.3(Gibbs 仕様)↗

相互作用 Φ=(ΦA)\Phi=(\Phi_A)(定義 2.1[熱力学的極限])が sup⁡x∑A∋x∥ΦA∥∞<∞\sup_{x}\sum_{A\ni x}\|\Phi_A\|_\infty<\infty を満たすとする。Λ⋐Zd\Lambda\Subset\mathbb{Z}^d と ω∈Ω\omega\in\Omega について

HΛ(ω)=∑A∩Λ≠∅ΦA(ω),ZΛω=∑ξ∈ΩΛe−HΛ(ξωΛc),H_\Lambda(\omega)=\sum_{A\cap\Lambda\ne\emptyset}\Phi_A(\omega),\qquad Z_\Lambda^\omega=\sum_{\xi\in\Omega_\Lambda}e^{-H_\Lambda(\xi\omega_{\Lambda^c})},γΛΦ(f∣ω)=1ZΛω∑ξ∈ΩΛe−HΛ(ξωΛc)f(ξωΛc)\gamma^\Phi_\Lambda(f\mid\omega)=\frac{1}{Z_\Lambda^\omega}\sum_{\xi\in\Omega_\Lambda}e^{-H_\Lambda(\xi\omega_{\Lambda^c})}f(\xi\omega_{\Lambda^c})

とおく。γΦ=(γΛΦ)\gamma^\Phi=(\gamma^\Phi_\Lambda) を Φ\Phi の Gibbs 仕様という。

HΛH_\Lambda の級数は ∑A∩Λ≠∅∥ΦA∥∞≤∣Λ∣sup⁡x∑A∋x∥ΦA∥∞\sum_{A\cap\Lambda\ne\emptyset}\|\Phi_A\|_\infty\le|\Lambda|\sup_x\sum_{A\ni x}\|\Phi_A\|_\infty で絶対収束します。境界条件 ωΛc\omega_{\Lambda^c} との相互作用も含めたエネルギーであることに注意してください(前章の HΛτH_\Lambda^\tau)。

命題 3.4(Gibbs 仕様は準局所的な仕様)✓

γΦ\gamma^\Phi は準局所的な仕様である。

証明(命題 3.4)

核であること・固有性. γΛΦ(⋅∣ω)\gamma^\Phi_\Lambda(\cdot\mid\omega) は有限個の点 ξωΛc\xi\omega_{\Lambda^c} に重みを置いた確率測度で、重みは ωΛc\omega_{\Lambda^c} だけに依存するので TΛ\mathcal{T}_\Lambda-可測です。B∈TΛB\in\mathcal{T}_\Lambda なら 1B(ξωΛc)=1B(ω)\mathbf{1}_B(\xi\omega_{\Lambda^c})=\mathbf{1}_B(\omega) なので固有性が成り立ちます。

整合性. Λ⊂Δ\Lambda\subset\Delta とし、Δ\Delta 上の配位を (α,ζ)(\alpha,\zeta)(α∈ΩΛ\alpha\in\Omega_\Lambda、ζ∈ΩΔ∖Λ\zeta\in\Omega_{\Delta\setminus\Lambda})と書きます。HΔ−HΛ=∑ΦAH_\Delta-H_\Lambda=\sum\Phi_A(和は A∩Δ≠∅A\cap\Delta\ne\emptyset、A∩Λ=∅A\cap\Lambda=\emptyset)は α\alpha に依存しないので、これを K(ζωΔc)K(\zeta\omega_{\Delta^c}) と書けば

e−HΔ(αζωΔc)=e−K(ζωΔc)e−HΛ(αζωΔc).e^{-H_\Delta(\alpha\zeta\omega_{\Delta^c})}=e^{-K(\zeta\omega_{\Delta^c})}e^{-H_\Lambda(\alpha\zeta\omega_{\Delta^c})}.

g=γΛΦfg=\gamma^\Phi_\Lambda f は α\alpha に依存しないので

γΔΦg(ω)=1ZΔω∑ζe−K(∑αe−HΛ(αζωΔc))1ZΛζωΔc∑α′e−HΛ(α′ζωΔc)f(α′ζωΔc)=1ZΔω∑ζ∑α′e−Ke−HΛ(α′ζωΔc)f(α′ζωΔc)=γΔΦf(ω).\begin{aligned} \gamma^\Phi_\Delta g(\omega) &=\frac1{Z^\omega_\Delta}\sum_\zeta e^{-K}\Bigl(\sum_\alpha e^{-H_\Lambda(\alpha\zeta\omega_{\Delta^c})}\Bigr)\frac{1}{Z_\Lambda^{\zeta\omega_{\Delta^c}}}\sum_{\alpha'}e^{-H_\Lambda(\alpha'\zeta\omega_{\Delta^c})}f(\alpha'\zeta\omega_{\Delta^c})\\ &=\frac1{Z^\omega_\Delta}\sum_\zeta\sum_{\alpha'}e^{-K}e^{-H_\Lambda(\alpha'\zeta\omega_{\Delta^c})}f(\alpha'\zeta\omega_{\Delta^c})=\gamma^\Phi_\Delta f(\omega). \end{aligned}

2 行目では括弧内の和が ZΛζωΔcZ_\Lambda^{\zeta\omega_{\Delta^c}} に等しいことを使いました。

準局所性. ω↦HΛ(ξωΛc)\omega\mapsto H_\Lambda(\xi\omega_{\Lambda^c}) は連続関数(局所関数)の一様収束級数なので連続です。ff が連続なら γΛΦf\gamma^\Phi_\Lambda f は連続関数の有限和を正の連続関数で割ったものなので連続です。

例 3.5(相互作用のない場合)

Φ=0\Phi=0 なら γΛ(⋅∣ω)\gamma_\Lambda(\cdot\mid\omega) は「Λ\Lambda 上は一様分布、外は ω\omega のまま」です。後で見るように、この仕様の Gibbs 測度は一様分布の無限直積だけです(演習 8.1)。

定義 4.1(Gibbs 測度(DLR 条件))↗

仕様 γ\gamma について、Ω\Omega 上の確率測度 μ\mu がすべての Λ⋐Zd\Lambda\Subset\mathbb{Z}^d で

μγΛ=μ\mu\gamma_\Lambda=\mu

を満たすとき、μ\mu を γ\gamma の Gibbs 測度という。これを DLR 方程式(Dobrushin–Lanford–Ruelle)と呼び、Gibbs 測度全体を G(γ)\mathcal{G}(\gamma) と書く。γ=γΦ\gamma=\gamma^\Phi のときは G(Φ)\mathcal{G}(\Phi) とも書く。

命題 4.2(DLR 方程式と条件付き期待値)✓

確率測度 μ\mu と Λ⋐Zd\Lambda\Subset\mathbb{Z}^d について、次は同値である。

(a) μγΛ=μ\mu\gamma_\Lambda=\mu。

(b) すべての有界可測な ff について Eμ[f∣TΛ]=γΛfE_\mu[f\mid\mathcal{T}_\Lambda]=\gamma_\Lambda f(μ\mu-a.s.)。

証明(命題 4.2)

(a)⇒(b): γΛf\gamma_\Lambda f は TΛ\mathcal{T}_\Lambda-可測で有界です。B∈TΛB\in\mathcal{T}_\Lambda について、補題 3.2 と (a) から

μ(1BγΛf)=μ(γΛ(1Bf))=μγΛ(1Bf)=μ(1Bf).\mu(\mathbf{1}_B\gamma_\Lambda f)=\mu(\gamma_\Lambda(\mathbf{1}_Bf))=\mu\gamma_\Lambda(\mathbf{1}_Bf)=\mu(\mathbf{1}_Bf).

これは条件付き期待値の定義(定義 3.1[条件付き期待値])の条件そのものです。

(b)⇒(a): B=ΩB=\Omega とすれば μ(γΛf)=μ(f)\mu(\gamma_\Lambda f)=\mu(f)、すなわち μγΛ(f)=μ(f)\mu\gamma_\Lambda(f)=\mu(f) です。

Gibbs 仕様では、(b) は「Λ\Lambda の外の配位 ωΛc\omega_{\Lambda^c} を与えたときの ωΛ\omega_\Lambda の条件付き分布が e−HΛ(ξωΛc)/ZΛωe^{-H_\Lambda(\xi\omega_{\Lambda^c})}/Z_\Lambda^\omega である」という意味です。条件付き期待値は μ\mu-a.s. にしか決まらないので、「すべての境界条件 ω\omega で」条件付き分布を指定する仕様は、測度より多くの情報を含んでいます。

注意 4.3(1 点の仕様で足りること)

Gibbs 仕様では、すべての 1 点集合 Λ={x}\Lambda=\{x\} について DLR 方程式が成り立てば、すべての有限集合で成り立つことが知られています(Georgii の教科書 第 1 章)。本章では使いません。

定理 5.1(Gibbs 測度の存在)✓

γ\gamma を準局所的な仕様とし、Λ1⊂Λ2⊂⋯\Lambda_1\subset\Lambda_2\subset\cdots を ⋃nΛn=Zd\bigcup_n\Lambda_n=\mathbb{Z}^d となる有限集合の列、ω(n)∈Ω\omega^{(n)}\in\Omega を任意の境界条件とする。μn=γΛn(⋅∣ω(n))\mu_n=\gamma_{\Lambda_n}(\cdot\mid\omega^{(n)}) とおく。

(i) (μn)(\mu_n) は弱収束する部分列をもつ。

(ii) (μn)(\mu_n) の部分列の弱収束極限はすべて G(γ)\mathcal{G}(\gamma) に属する。

特に G(γ)≠∅\mathcal{G}(\gamma)\ne\emptyset である。

証明(定理 5.1)

(i) Ω\Omega はコンパクト距離空間(補題 2.1)なので、確率測度の族はすべて緊密です(定義 3.1[確率測度の弱収束] で K=ΩK=\Omega)。Prokhorov の定理(定理 3.3[確率測度の弱収束])から弱収束する部分列があります。

(ii) μnk⇒μ\mu_{n_k}\Rightarrow\mu とし、Λ⋐Zd\Lambda\Subset\mathbb{Z}^d を固定します。Λ\Lambda は有限なので、ある k0k_0 から先 Λ⊂Λnk\Lambda\subset\Lambda_{n_k} です。そのとき整合性から

μnkγΛ(f)=∫γΛf dγΛnk(⋅∣ω(nk))=(γΛnkγΛf)(ω(nk))=(γΛnkf)(ω(nk))=μnk(f).\mu_{n_k}\gamma_\Lambda(f)=\int\gamma_\Lambda f\,d\gamma_{\Lambda_{n_k}}(\cdot\mid\omega^{(n_k)})=(\gamma_{\Lambda_{n_k}}\gamma_\Lambda f)(\omega^{(n_k)})=(\gamma_{\Lambda_{n_k}}f)(\omega^{(n_k)})=\mu_{n_k}(f).

f∈C(Ω)f\in C(\Omega) なら準局所性から γΛf∈C(Ω)\gamma_\Lambda f\in C(\Omega) なので、弱収束の定義(定義 2.1[確率測度の弱収束])から

μ(γΛf)=lim⁡kμnk(γΛf)=lim⁡kμnk(f)=μ(f).\mu(\gamma_\Lambda f)=\lim_k\mu_{n_k}(\gamma_\Lambda f)=\lim_k\mu_{n_k}(f)=\mu(f).

2 つの確率測度 μγΛ\mu\gamma_\Lambda と μ\mu がすべての連続関数で一致するので(補題 2.2[確率測度の弱収束])μγΛ=μ\mu\gamma_\Lambda=\mu です。

この証明でコンパクト性を使ったのは (i) の緊密性だけです。スピンが R\mathbb{R} に値をとる格子場(たとえば格子 ϕ4\phi^4 模型)では Ω=RZd\Omega=\mathbb{R}^{\mathbb{Z}^d} はコンパクトでなく、境界条件と相互作用について一様な積率評価を示して緊密性を確保する必要があります(Ruelle の超安定性評価など)。構成的場の理論で無限体積極限を作る議論(相関不等式 の単調性による方法など)は、この緊密性を別の手段で得るものと見ることができます。

系 5.2(Gibbs 測度の集合は凸でコンパクト)✓

準局所的な仕様 γ\gamma について、G(γ)\mathcal{G}(\gamma) は凸で、弱収束について閉じている(したがって点列コンパクトである)。

証明(系 5.2)

μγΛ\mu\gamma_\Lambda は μ\mu について線形なので、DLR 方程式は凸結合で保たれます。μn∈G(γ)\mu_n\in\mathcal{G}(\gamma)、μn⇒μ\mu_n\Rightarrow\mu なら、定理 5.1 (ii) の証明と同じく、f∈C(Ω)f\in C(\Omega) について μ(γΛf)=lim⁡μn(γΛf)=lim⁡μn(f)=μ(f)\mu(\gamma_\Lambda f)=\lim\mu_n(\gamma_\Lambda f)=\lim\mu_n(f)=\mu(f) で、μ∈G(γ)\mu\in\mathcal{G}(\gamma)。点列コンパクト性は Prokhorov の定理から従います。

この節のねらいは、「純粋な相」(混ぜ物でない状態)を数学の言葉で捉えることです。直観では、++ 相と −- 相を半分ずつ混ぜた状態は、どちらの相にいるかという「無限遠の情報」について確率 1/21/2 の不確かさを持ちます。つまり混合であることと、尾部事象の確率が 00 か 11 でないことは同じ現象で、これを次の補題と定理で正確にします。

凸集合 G(γ)\mathcal{G}(\gamma) の元 μ\mu が極値的とは、μ=tν1+(1−t)ν2\mu=t\nu_1+(1-t)\nu_2(0<t<10<t<1、νi∈G(γ)\nu_i\in\mathcal{G}(\gamma))ならば ν1=ν2=μ\nu_1=\nu_2=\mu となることです。極値 Gibbs 測度全体を ex⁡G(γ)\operatorname{ex}\mathcal{G}(\gamma) と書きます。

次の補題は、極値性と尾部自明性を結ぶ道具です。ある Gibbs 測度が別の Gibbs 測度を「含む」とき、その重みの関数は有限個のスピンの値に左右されない、つまり尾部の情報だけで決まる、と述べています。

補題 6.1(Gibbs 測度どうしの密度は尾部可測)✓

μ,ν∈G(γ)\mu,\nu\in\mathcal{G}(\gamma) で ν≪μ\nu\ll\mu ならば、密度 dν/dμd\nu/d\mu は μ\mu-a.s. で T\mathcal{T}-可測な関数に等しい。

証明(補題 6.1)

h=dν/dμh=d\nu/d\mu(Radon–Nikodym の定理、定理 2.1[条件付き期待値])とし、Λ⋐Zd\Lambda\Subset\mathbb{Z}^d を固定して hΛ=Eμ[h∣TΛ]h_\Lambda=E_\mu[h\mid\mathcal{T}_\Lambda] とおきます。有界可測な ff について、ν,μ∈G(γ)\nu,\mu\in\mathcal{G}(\gamma)、γΛf\gamma_\Lambda f の TΛ\mathcal{T}_\Lambda-可測性、補題 3.2 を順に使うと

ν(f)=ν(γΛf)=μ(h γΛf)=μ(hΛ γΛf)=μ(γΛ(hΛf))=μ(hΛf).\nu(f)=\nu(\gamma_\Lambda f)=\mu(h\,\gamma_\Lambda f)=\mu(h_\Lambda\,\gamma_\Lambda f)=\mu\bigl(\gamma_\Lambda(h_\Lambda f)\bigr)=\mu(h_\Lambda f).

3 番目の等号は既知量の取り出し(定理 4.4[条件付き期待値]、γΛf\gamma_\Lambda f は有界 TΛ\mathcal{T}_\Lambda-可測)、最後は μγΛ=μ\mu\gamma_\Lambda=\mu です。4 番目と最後の等号では hΛfh_\Lambda f が有界とは限りませんが、f≥0f\ge0 としておけば、補題 3.2 と μ(γΛg)=μ(g)\mu(\gamma_\Lambda g)=\mu(g) は g=(hΛ∧M)fg=(h_\Lambda\wedge M)f に対して成り立ち、M→∞M\to\infty で単調収束定理(定理 2.1[確率変数と期待値])により g=hΛfg=h_\Lambda f に広がります(一般の有界 ff は正負に分けます)。よって hΛh_\Lambda も ν\nu の μ\mu に関する密度で、密度の一意性から h=hΛh=h_\Lambda(μ\mu-a.s.)です。

Λn↑Zd\Lambda_n\uparrow\mathbb{Z}^d をとり、各 hΛnh_{\Lambda_n} の TΛn\mathcal{T}_{\Lambda_n}-可測な版を固定して h~=lim sup⁡nhΛn\tilde h=\limsup_nh_{\Lambda_n} とおきます。TΛn\mathcal{T}_{\Lambda_n} は nn について減少するので、mm を固定すると n≥mn\ge m の hΛnh_{\Lambda_n} はすべて TΛm\mathcal{T}_{\Lambda_m}-可測で、h~\tilde h は TΛm\mathcal{T}_{\Lambda_m}-可測です。任意の有限 Λ\Lambda はある Λm\Lambda_m に含まれ TΛm⊂TΛ\mathcal{T}_{\Lambda_m}\subset\mathcal{T}_\Lambda なので、h~\tilde h は T\mathcal{T}-可測です。可算個の零集合を除いて hΛn=hh_{\Lambda_n}=h なので h~=h\tilde h=h(μ\mu-a.s.)です。

定理 6.2(極値 Gibbs 測度の特徴づけ)✓

μ∈G(γ)\mu\in\mathcal{G}(\gamma) について、次は同値である。

(a) μ\mu は G(γ)\mathcal{G}(\gamma) の極値点である。

(b) μ\mu は尾部自明である。すなわちすべての B∈TB\in\mathcal{T} で μ(B)∈{0,1}\mu(B)\in\{0,1\}。

証明(定理 6.2)

(b)⇒(a). μ=tν1+(1−t)ν2\mu=t\nu_1+(1-t)\nu_2(0<t<10<t<1、νi∈G(γ)\nu_i\in\mathcal{G}(\gamma))とします。ν1≤t−1μ\nu_1\le t^{-1}\mu なので ν1≪μ\nu_1\ll\mu で、補題 6.1 から密度 hh は T\mathcal{T}-可測です。尾部自明な測度に関して T\mathcal{T}-可測な関数は a.s. 定数です({h>c}∈T\{h>c\}\in\mathcal{T} の確率が 00 か 11 なので、c∗=sup⁡{c:μ(h>c)=1}c^*=\sup\{c:\mu(h>c)=1\} とすれば h=c∗h=c^* a.s.)。μ(h)=ν1(Ω)=1\mu(h)=\nu_1(\Omega)=1 から h=1h=1、すなわち ν1=μ\nu_1=\mu。同様に ν2=μ\nu_2=\mu です。

(a)⇒(b). B∈TB\in\mathcal{T}、0<μ(B)<10<\mu(B)<1 と仮定し、ν=μ(⋅∣B)\nu=\mu(\cdot\mid B)、ν′=μ(⋅∣Bc)\nu'=\mu(\cdot\mid B^c) とおきます。B∈TΛB\in\mathcal{T}_\Lambda なので 補題 3.2 から、有界可測な ff について

ν(γΛf)=μ(1BγΛf)μ(B)=μ(γΛ(1Bf))μ(B)=μ(1Bf)μ(B)=ν(f).\nu(\gamma_\Lambda f)=\frac{\mu(\mathbf{1}_B\gamma_\Lambda f)}{\mu(B)}=\frac{\mu(\gamma_\Lambda(\mathbf{1}_Bf))}{\mu(B)}=\frac{\mu(\mathbf{1}_Bf)}{\mu(B)}=\nu(f).

よって ν∈G(γ)\nu\in\mathcal{G}(\gamma)、同様に ν′∈G(γ)\nu'\in\mathcal{G}(\gamma) で、μ=μ(B)ν+μ(Bc)ν′\mu=\mu(B)\nu+\mu(B^c)\nu' です。ν(B)=1≠μ(B)\nu(B)=1\ne\mu(B) なので ν≠μ\nu\ne\mu となり、μ\mu は極値的でありません。

系 6.3(尾部での値が Gibbs 測度を決める)✓

(i) μ,ν∈G(γ)\mu,\nu\in\mathcal{G}(\gamma) が T\mathcal{T} の上で一致すれば μ=ν\mu=\nu。

(ii) μ≠ν\mu\ne\nu がともに極値 Gibbs 測度ならば、μ(B)=1\mu(B)=1、ν(B)=0\nu(B)=0 となる B∈TB\in\mathcal{T} がある。特に μ\mu と ν\nu は互いに特異である。

証明(系 6.3)

(i) ρ=(μ+ν)/2∈G(γ)\rho=(\mu+\nu)/2\in\mathcal{G}(\gamma) について μ,ν≪ρ\mu,\nu\ll\rho で、補題 6.1 から密度 h,h′h,h' は T\mathcal{T}-可測にとれます。B∈TB\in\mathcal{T} で ρ(h1B)=μ(B)=ν(B)=ρ(h′1B)\rho(h\mathbf{1}_B)=\mu(B)=\nu(B)=\rho(h'\mathbf{1}_B) なので、B={h>h′}B=\{h>h'\}、{h<h′}\{h<h'\} とおけば h=h′h=h'(ρ\rho-a.s.)、よって μ=ν\mu=\nu。

(ii) (i) から μ(B)≠ν(B)\mu(B)\ne\nu(B) となる B∈TB\in\mathcal{T} があり、定理 6.2 から値は 00 か 11 です。必要なら BB を補集合に替えます。

極値 Gibbs 測度は、尾部の情報(「無限遠での境界条件」)が確定した状態であり、物理的には純粋な相と解釈されます。たとえば低温の Ising 模型の ++ 相と −- 相は、尾部事象「空間平均磁化が正」で区別されます。

注意 6.4(極値分解と境界条件からの再構成)

本章では証明しませんが、次の 2 つが成り立ちます(Georgii の教科書 第 7 章)。(1) 各 μ∈G(γ)\mu\in\mathcal{G}(\gamma) は、ex⁡G(γ)\operatorname{ex}\mathcal{G}(\gamma) 上のただ一つの確率測度 wμw_\mu による混合 μ=∫ν wμ(dν)\mu=\int\nu\,w_\mu(d\nu) で書ける(極値分解)。(2) μ\mu が極値的なら、μ\mu-a.e. の ω\omega について γΛn(⋅∣ω)⇒μ\gamma_{\Lambda_n}(\cdot\mid\omega)\Rightarrow\mu。(2) は逆向きマルチンゲールの収束定理と、尾部自明性から出ます。極値 Gibbs 測度は「典型的な境界条件から体積を広げて得られる測度」です。

この節のねらいは、Gibbs 測度がただ一つに決まるための、確かめやすい十分条件を与えることです。考え方は、スピンが遠くのスピンから受ける影響の強さを数 CxyC_{xy} で測り、その総和が 11 未満なら、境界条件の違いが遠くへ行くほど薄まって消える、というものです。

SS 上の確率測度 p,qp,q の全変動距離を

∥p−q∥=max⁡B⊂S∣p(B)−q(B)∣=12∑s∈S∣p(s)−q(s)∣\|p-q\|=\max_{B\subset S}|p(B)-q(B)|=\frac12\sum_{s\in S}|p(s)-q(s)|

とします。2 番目の等号は、最大値が B={p>q}B=\{p>q\} でとられることからわかります。仕様 γ\gamma の 1 点集合 {x}\{x\} の核から、xx のスピンの条件付き分布 γx(⋅∣ω)=γ{x}(ωx∈⋅∣ω)\gamma_x(\cdot\mid\omega)=\gamma_{\{x\}}(\omega_x\in\cdot\mid\omega) が決まります。固有性から γ{x}(⋅∣ω)\gamma_{\{x\}}(\cdot\mid\omega) は {η:ηxc=ωxc}\{\eta:\eta_{x^c}=\omega_{x^c}\} に集中し、γx(⋅∣ω)\gamma_x(\cdot\mid\omega) は ωxc\omega_{x^c} だけに依存します。したがって

γ{x}f(ω)=∑s∈Sγx(s∣ω) f(s ωxc)\gamma_{\{x\}}f(\omega)=\sum_{s\in S}\gamma_x(s\mid\omega)\,f(s\,\omega_{x^c})

です(s ωxcs\,\omega_{x^c} は xx で ss、他で ω\omega に等しい配位)。

定義 7.1(Dobrushin の相互作用行列)↗

x≠yx\ne y について

Cxy=sup⁡{∥γx(⋅∣ω)−γx(⋅∣η)∥: ω,η∈Ω, ωz=ηz (z≠y)},C_{xy}=\sup\bigl\{\|\gamma_x(\cdot\mid\omega)-\gamma_x(\cdot\mid\eta)\|:\ \omega,\eta\in\Omega,\ \omega_z=\eta_z\ (z\ne y)\bigr\},

Cxx=0C_{xx}=0 とおく。α(γ)=sup⁡x∑yCxy\alpha(\gamma)=\sup_x\sum_yC_{xy} を Dobrushin の定数という。

CxyC_{xy} は「yy のスピン 1 個を変えたとき、xx の条件付き分布がどれだけ動くか」の最大値です。関数 ff の側では、yy のスピン 1 個を変えたときの変動 δy(f)=sup⁡{∣f(ω)−f(η)∣:ωz=ηz (z≠y)}\delta_y(f)=\sup\{|f(\omega)-f(\eta)|:\omega_z=\eta_z\ (z\ne y)\} を考え、

D={f∈C(Ω): ∑yδy(f)<∞}\mathcal{D}=\Bigl\{f\in C(\Omega):\ \sum_y\delta_y(f)<\infty\Bigr\}

とします。局所関数は D\mathcal{D} に属します(有限個の yy でしか δy(f)≠0\delta_y(f)\ne0 にならない)。

補題 7.2(変動の評価)✓

(i) f∈C(Ω)f\in C(\Omega) ならば sup⁡f−inf⁡f≤∑yδy(f)\sup f-\inf f\le\sum_y\delta_y(f)。

(ii) γ\gamma が準局所的な仕様で f∈Df\in\mathcal{D}、x∈Zdx\in\mathbb{Z}^d ならば、δx(γ{x}f)=0\delta_x(\gamma_{\{x\}}f)=0、y≠xy\ne x について

δy(γ{x}f)≤δy(f)+Cxy δx(f).\delta_y(\gamma_{\{x\}}f)\le\delta_y(f)+C_{xy}\,\delta_x(f).

特に ∑yCxy<∞\sum_yC_{xy}<\infty なら γ{x}f∈D\gamma_{\{x\}}f\in\mathcal{D}。

証明(補題 7.2)

(i) ω,η∈Ω\omega,\eta\in\Omega と Λn↑Zd\Lambda_n\uparrow\mathbb{Z}^d について、ω\omega から Λn\Lambda_n の点を 1 個ずつ η\eta の値に変えていくと ∣f(ω)−f(ηΛnωΛnc)∣≤∑y∈Λnδy(f)|f(\omega)-f(\eta_{\Lambda_n}\omega_{\Lambda_n^c})|\le\sum_{y\in\Lambda_n}\delta_y(f)。ηΛnωΛnc→η\eta_{\Lambda_n}\omega_{\Lambda_n^c}\to\eta なので、ff の連続性から ∣f(ω)−f(η)∣≤∑yδy(f)|f(\omega)-f(\eta)|\le\sum_y\delta_y(f)。

(ii) γ{x}f(ω)\gamma_{\{x\}}f(\omega) は ωxc\omega_{x^c} だけに依存するので δx(γ{x}f)=0\delta_x(\gamma_{\{x\}}f)=0。y≠xy\ne x、ωz=ηz\omega_z=\eta_z(z≠yz\ne y)とし、p=γx(⋅∣ω)p=\gamma_x(\cdot\mid\omega)、q=γx(⋅∣η)q=\gamma_x(\cdot\mid\eta) とおくと

γ{x}f(ω)−γ{x}f(η)=∑sp(s)(f(sωxc)−f(sηxc))+∑s(p(s)−q(s))f(sηxc).\gamma_{\{x\}}f(\omega)-\gamma_{\{x\}}f(\eta)=\sum_sp(s)\bigl(f(s\omega_{x^c})-f(s\eta_{x^c})\bigr)+\sum_s\bigl(p(s)-q(s)\bigr)f(s\eta_{x^c}).

第 1 項の絶対値は δy(f)\delta_y(f) 以下です(sωxcs\omega_{x^c} と sηxcs\eta_{x^c} は yy でしか違わない)。第 2 項では、g(s)=f(sηxc)g(s)=f(s\eta_{x^c}) の最大値と最小値の中点を cc とすれば ∑(p−q)=0\sum(p-q)=0 から

∣∑s(p(s)−q(s))(g(s)−c)∣≤∑s∣p(s)−q(s)∣⋅max⁡g−min⁡g2≤∥p−q∥ δx(f)≤Cxyδx(f).\Bigl|\sum_s(p(s)-q(s))(g(s)-c)\Bigr|\le\sum_s|p(s)-q(s)|\cdot\frac{\max g-\min g}{2}\le\|p-q\|\,\delta_x(f)\le C_{xy}\delta_x(f).

最後に、準局所性から γ{x}f\gamma_{\{x\}}f は連続で、変動の和は ∑yδy(f)+δx(f)∑yCxy<∞\sum_y\delta_y(f)+\delta_x(f)\sum_yC_{xy}<\infty です。

定理 7.3(Dobrushin の一意性定理)✓

γ\gamma を準局所的な仕様とし、α(γ)<1\alpha(\gamma)<1 とする。このとき G(γ)\mathcal{G}(\gamma) はちょうど 1 点からなる。

証明の方針は次のとおりです。2 つの Gibbs 測度 μ,ν\mu,\nu の期待値の差を、各座標の変動 δy(f)\delta_y(f) の重みつき和で上から押さえます(この重みが「許容的」な bb です)。1 点を再標本化するたびに、その点の重みが、それまでの重みの最大値の α\alpha 倍以下に小さくなるので、これを繰り返せば重みを限りなく 00 に近づけられ、差が 00 になります。

証明(定理 7.3)

存在は 定理 5.1 なので、μ,ν∈G(γ)\mu,\nu\in\mathcal{G}(\gamma) について μ=ν\mu=\nu を示します。α=α(γ)\alpha=\alpha(\gamma) とおきます。有界な関数 b:Zd→[0,∞)b:\mathbb{Z}^d\to[0,\infty) が

∣μ(f)−ν(f)∣≤∑yby δy(f)(すべての f∈D)|\mu(f)-\nu(f)|\le\sum_yb_y\,\delta_y(f)\qquad(\text{すべての }f\in\mathcal{D})

を満たすとき、bb を許容的と呼びます。補題 7.2 (i) から ∣μ(f)−ν(f)∣≤sup⁡f−inf⁡f≤∑yδy(f)|\mu(f)-\nu(f)|\le\sup f-\inf f\le\sum_y\delta_y(f) なので、定数 11 は許容的です。

1 点の改良. bb が許容的で x∈Zdx\in\mathbb{Z}^d とし、b′b' を bx′=∑yCxybyb'_x=\sum_yC_{xy}b_y、bz′=bzb'_z=b_z(z≠xz\ne x)で定めます。f∈Df\in\mathcal{D} について、DLR 方程式 μγ{x}=μ\mu\gamma_{\{x\}}=\mu、νγ{x}=ν\nu\gamma_{\{x\}}=\nu と 補題 7.2 (ii) から

∣μ(f)−ν(f)∣=∣μ(γ{x}f)−ν(γ{x}f)∣≤∑y≠xby(δy(f)+Cxyδx(f))=∑yby′δy(f).|\mu(f)-\nu(f)|=|\mu(\gamma_{\{x\}}f)-\nu(\gamma_{\{x\}}f)|\le\sum_{y\ne x}b_y\bigl(\delta_y(f)+C_{xy}\delta_x(f)\bigr)=\sum_yb'_y\delta_y(f).

よって b′b' も許容的です。

有限集合での改良. β0=sup⁡yby\beta_0=\sup_yb_y とし、有限集合 Λ={x1,…,xN}\Lambda=\{x_1,\dots,x_N\} の点で順に 1 点の改良を行います。改良された値は ∑yCxiyby≤αβ0≤β0\sum_yC_{x_iy}b_y\le\alpha\beta_0\le\beta_0 なので、各段階で上限は β0\beta_0 以下に保たれ、改良済みの座標は αβ0\alpha\beta_0 以下、改良していない座標は β0\beta_0 以下です。したがって、f∈Df\in\mathcal{D} について

∣μ(f)−ν(f)∣≤αβ0∑y∈Λδy(f)+β0∑y∉Λδy(f).|\mu(f)-\nu(f)|\le\alpha\beta_0\sum_{y\in\Lambda}\delta_y(f)+\beta_0\sum_{y\notin\Lambda}\delta_y(f).

Λ↑Zd\Lambda\uparrow\mathbb{Z}^d とすれば第 2 項は 00 に収束するので、定数 αβ0\alpha\beta_0 も許容的です。

結論. 定数 11 から始めて繰り返せば、すべての kk で定数 αk\alpha^k が許容的で、∣μ(f)−ν(f)∣≤αk∑yδy(f)→0|\mu(f)-\nu(f)|\le\alpha^k\sum_y\delta_y(f)\to0。局所関数は D\mathcal{D} に属するので μ\mu と ν\nu は局所関数の期待値が一致し、補題 2.1 (iii) と 補題 2.2[確率測度の弱収束] から μ=ν\mu=\nu です。

証明の構造は、Markov 連鎖のエルゴード定理を全変動距離の縮小で示す議論(Markov 連鎖と Poisson 過程 の 補題 3.1[Markov 連鎖と Poisson 過程])と同じ発想です。1 点ずつ再標本化するたびに、2 つの測度の「違い」が係数 α\alpha 倍に縮みます。

系 7.4(Gibbs 仕様での判定)✓

Φ\Phi を 定義 3.3 の相互作用とし、osc⁡ΦA=sup⁡ΦA−inf⁡ΦA\operatorname{osc}\Phi_A=\sup\Phi_A-\inf\Phi_A とする。

sup⁡x∑A∋x(∣A∣−1)osc⁡ΦA<12\sup_x\sum_{A\ni x}(|A|-1)\operatorname{osc}\Phi_A<\frac12

ならば ∣G(Φ)∣=1|\mathcal{G}(\Phi)|=1。

証明(系 7.4)

x≠yx\ne y、ωz=ηz\omega_z=\eta_z(z≠yz\ne y)とします。γx(s∣ω)∝e−H{x}(sωxc)\gamma_x(s\mid\omega)\propto e^{-H_{\{x\}}(s\omega_{x^c})} で、u(s)=H{x}(sηxc)−H{x}(sωxc)u(s)=H_{\{x\}}(s\eta_{x^c})-H_{\{x\}}(s\omega_{x^c}) には x,yx,y の両方を含む AA だけが寄与するので、∣u(s)∣≤rxy:=∑A∋x,yosc⁡ΦA|u(s)|\le r_{xy}:=\sum_{A\ni x,y}\operatorname{osc}\Phi_A です。p=γx(⋅∣ω)p=\gamma_x(\cdot\mid\omega)、q=γx(⋅∣η)q=\gamma_x(\cdot\mid\eta) とおくと q(s)=p(s)e−u(s)/∑tp(t)e−u(t)≥e−2rxyp(s)q(s)=p(s)e^{-u(s)}/\sum_tp(t)e^{-u(t)}\ge e^{-2r_{xy}}p(s) なので

∥p−q∥=∑s(p(s)−q(s))+≤∑sp(s)(1−e−2rxy)≤2rxy.\|p-q\|=\sum_s(p(s)-q(s))^+\le\sum_sp(s)(1-e^{-2r_{xy}})\le2r_{xy}.

よって ∑yCxy≤2∑y≠x∑A∋x,yosc⁡ΦA=2∑A∋x(∣A∣−1)osc⁡ΦA<1\sum_yC_{xy}\le2\sum_{y\ne x}\sum_{A\ni x,y}\operatorname{osc}\Phi_A=2\sum_{A\ni x}(|A|-1)\operatorname{osc}\Phi_A<1 で、定理 7.3 を使います。

命題 7.5(Ising 模型の高温での一意性)✓

S={−1,1}S=\{-1,1\} とし、最近接 Ising 模型の相互作用 Φ{x,y}(ω)=−βJωxωy\Phi_{\{x,y\}}(\omega)=-\beta J\omega_x\omega_y(∣x−y∣=1|x-y|=1)、Φ{x}(ω)=−βhωx\Phi_{\{x\}}(\omega)=-\beta h\omega_x(β>0\beta\gt0、J,h∈RJ,h\in\mathbb{R})、それ以外の ΦA=0\Phi_A=0 を考える。2dtanh⁡(β∣J∣)<12d\tanh(\beta|J|)\lt1 ならば、hh によらず ∣G(Φ)∣=1|\mathcal{G}(\Phi)|=1 である。

証明(命題 7.5)

xx を固定し、局所場を m(ω)=J∑y∼xωy+hm(\omega)=J\sum_{y\sim x}\omega_y+h とします。H{x}(sωxc)=−βsm(ω)+c(ω)H_{\{x\}}(s\omega_{x^c})=-\beta sm(\omega)+c(\omega)(c(ω)c(\omega) は ss によらない)なので、定義 3.3 から

γx(+1∣ω)=eβmeβm+e−βm=1+tanh⁡βm2,γx(−1∣ω)=1−tanh⁡βm2\gamma_x(+1\mid\omega)=\frac{e^{\beta m}}{e^{\beta m}+e^{-\beta m}}=\frac{1+\tanh\beta m}2,\qquad\gamma_x(-1\mid\omega)=\frac{1-\tanh\beta m}2

です。γx(⋅∣ω)\gamma_x(\cdot\mid\omega) は ω\omega の xx の隣接点での値だけに依存するので、y≁xy\not\sim x なら Cxy=0C_{xy}=0 です。y∼xy\sim x で ω,η\omega,\eta が yy でだけ異なるとき、局所場の差は ±2J\pm2J なので、2 点上の確率測度の全変動距離 ∥p−q∥=∑s(p(s)−q(s))+=∣p(+1)−q(+1)∣\|p-q\|=\sum_s(p(s)-q(s))^+=|p(+1)-q(+1)| から

Cxy≤sup⁡m∈R12∣tanh⁡(β(m+2∣J∣))−tanh⁡(βm)∣C_{xy}\le\sup_{m\in\mathbb{R}}\tfrac12\bigl|\tanh(\beta(m+2|J|))-\tanh(\beta m)\bigr|

です(J<0J\lt0 のときは mm を m−2∣J∣m-2|J| に置き換える)。a=2β∣J∣≥0a=2\beta|J|\ge0、g(t)=tanh⁡(t+a)−tanh⁡t≥0g(t)=\tanh(t+a)-\tanh t\ge0 とおくと g′(t)=cosh⁡−2(t+a)−cosh⁡−2tg'(t)=\cosh^{-2}(t+a)-\cosh^{-2}t は ∣t+a∣<∣t∣|t+a|\lt|t|、すなわち t<−a/2t\lt-a/2 で正、t>−a/2t\gt-a/2 で負なので、gg の最大値は g(−a/2)=2tanh⁡(a/2)g(-a/2)=2\tanh(a/2) です。よって Cxy≤tanh⁡(β∣J∣)C_{xy}\le\tanh(\beta|J|) で、隣接点は 2d2d 個なので α(γΦ)≤2dtanh⁡(β∣J∣)<1\alpha(\gamma^\Phi)\le2d\tanh(\beta|J|)\lt1 です。γΦ\gamma^\Phi は準局所的な仕様(命題 3.4。sup⁡x∑A∋x∥ΦA∥∞=β(2d∣J∣+∣h∣)<∞\sup_x\sum_{A\ni x}\|\Phi_A\|_\infty=\beta(2d|J|+|h|)\lt\infty)なので、定理 7.3 から ∣G(Φ)∣=1|\mathcal{G}(\Phi)|=1 です。

例 7.6(Ising 模型の高温での一意性)

最近接 Ising 模型(逆温度 β\beta、結合 JJ、外場 hh)では、xx の局所場を m=J∑y∼xωy+hm=J\sum_{y\sim x}\omega_y+h とすると γx(+1∣ω)=(1+tanh⁡βm)/2\gamma_x(+1\mid\omega)=(1+\tanh\beta m)/2 です。隣の yy のスピンを反転すると mm は m±2Jm\pm2J に変わるので

Cxy≤sup⁡m∈R12∣tanh⁡(β(m+2J))−tanh⁡(βm)∣=tanh⁡(β∣J∣)C_{xy}\le\sup_{m\in\mathbb{R}}\tfrac12\bigl|\tanh(\beta(m+2J))-\tanh(\beta m)\bigr|=\tanh(\beta|J|)

(上限は m=−Jm=-J でとられる)。隣接点は 2d2d 個なので α≤2dtanh⁡(β∣J∣)\alpha\le2d\tanh(\beta|J|) で、

2dtanh⁡(β∣J∣)<12d\tanh(\beta|J|)<1

なら外場 hh によらず Gibbs 測度はただ一つです(命題 7.5)。d=2d=2 では β∣J∣<artanh⁡(1/4)≈0.255\beta|J|<\operatorname{artanh}(1/4)\approx0.255 です(正確な臨界値は βcJ=12log⁡(1+2)≈0.441\beta_cJ=\frac12\log(1+\sqrt2)\approx0.441)。

注意 7.7(Dobrushin 条件のもとでの他の帰結)

Dobrushin 条件のもとでは、一意性だけでなく、相関の減衰(∣μ(fg)−μ(f)μ(g)∣|\mu(fg)-\mu(f)\mu(g)| の評価)や、相互作用についての Gibbs 測度の滑らかさも同じ行列 CC で制御できます。また、スピン系の Glauber 力学の指数的な緩和や、関数不等式の章 で扱う対数 Sobolev 不等式も、Dobrushin 型の条件から導かれます(Stroock–Zegarlinski)。

演習 8.1易

Φ=0\Phi=0 の Gibbs 仕様(例 3.5)について、G\mathcal{G} が SS 上の一様分布の無限直積だけからなることを、(1) DLR 方程式から直接、(2) 定理 7.3 から、の 2 通りで示してください。

解答

(1) μ∈G\mu\in\mathcal{G}、Λ⋐Zd\Lambda\Subset\mathbb{Z}^d、ξ∈ΩΛ\xi\in\Omega_\Lambda について μ(ωΛ=ξ)=μγΛ(ωΛ=ξ)=∫∣S∣−∣Λ∣ dμ=∣S∣−∣Λ∣\mu(\omega_\Lambda=\xi)=\mu\gamma_\Lambda(\omega_\Lambda=\xi)=\int|S|^{-|\Lambda|}\,d\mu=|S|^{-|\Lambda|}。有限次元分布が一様分布の直積と一致するので、柱集合(π-系)の上で一致し、測度の一意性(定理 4.5[測度空間と Carathéodory の拡張定理])から μ\mu は一様分布の無限直積です。逆にこの直積が DLR 方程式を満たすことは、Λ\Lambda 上の一様分布と外の分布が独立であることから従います。

(2) γx(⋅∣ω)\gamma_x(\cdot\mid\omega) は ω\omega によらない一様分布なので Cxy=0C_{xy}=0、α=0<1\alpha=0<1 です。

演習 8.2標準

μ1≠μ2\mu_1\ne\mu_2 を極値 Gibbs 測度とし、μ=12(μ1+μ2)\mu=\frac12(\mu_1+\mu_2) とします。系 6.3 の尾部事象 BB について μ(B)=1/2\mu(B)=1/2 を示し、μ\mu が局所関数の相関について「遠くで無相関」にならない例を作ってください。具体的には S={−1,1}S=\{-1,1\} で μ1(ωx)=m=−μ2(ωx)\mu_1(\omega_x)=m=-\mu_2(\omega_x)(m>0m>0、すべての xx)のとき、μ(ω0ωx)−μ(ω0)μ(ωx)≥m2\mu(\omega_0\omega_x)-\mu(\omega_0)\mu(\omega_x)\ge m^2 を示してください(極値測度では μi(ω0ωx)≥μi(ω0)μi(ωx)\mu_i(\omega_0\omega_x)\ge\mu_i(\omega_0)\mu_i(\omega_x) が成り立つと仮定してよい)。

解答

μ1(B)=1\mu_1(B)=1、μ2(B)=0\mu_2(B)=0 なので μ(B)=1/2\mu(B)=1/2。相関については μ(ω0)=μ(ωx)=0\mu(\omega_0)=\mu(\omega_x)=0 で、

μ(ω0ωx)=12(μ1(ω0ωx)+μ2(ω0ωx))≥12(m2+m2)=m2.\mu(\omega_0\omega_x)=\tfrac12\bigl(\mu_1(\omega_0\omega_x)+\mu_2(\omega_0\omega_x)\bigr)\ge\tfrac12(m^2+m^2)=m^2.

したがって共分散は ∣x∣|x| によらず m2m^2 以上で、∣x∣→∞|x|\to\infty でも 00 に近づきません。混合状態では、遠く離れたスピンが「どちらの相にいるか」という尾部の情報を通じて相関します。仮定した不等式は強磁性 Ising 模型の ±\pm 状態では GKS 不等式(定理 3.2[相関不等式])の極限として成り立ちます。

演習 8.3標準

S={−1,1}S=\{-1,1\}、Φ{x,y}=−Jxyωxωy\Phi_{\{x,y\}}=-J_{xy}\omega_x\omega_y(Jxy=JyxJ_{xy}=J_{yx}、sup⁡x∑y∣Jxy∣<∞\sup_x\sum_y|J_{xy}|<\infty)と外場 Φ{x}=−hxωx\Phi_{\{x\}}=-h_x\omega_x からなる Ising 模型(逆温度は J,hJ,h に含める)について、Cxy≤tanh⁡∣Jxy∣C_{xy}\le\tanh|J_{xy}| を示し、sup⁡x∑ytanh⁡∣Jxy∣<1\sup_x\sum_y\tanh|J_{xy}|<1 なら Gibbs 測度が一意であることを導いてください。系 7.4 の条件と比べてどちらが強いですか。

解答

例 7.6 と同じく、xx の局所場 m=∑yJxyωy+hxm=\sum_yJ_{xy}\omega_y+h_x について γx(+1∣ω)=(1+tanh⁡m)/2\gamma_x(+1\mid\omega)=(1+\tanh m)/2 で、yy のスピンを反転すると mm は ±2Jxy\pm2J_{xy} だけ変わります。tanh⁡\tanh の差は sup⁡m12∣tanh⁡(m+2a)−tanh⁡m∣=tanh⁡∣a∣\sup_m\frac12|\tanh(m+2a)-\tanh m|=\tanh|a| なので Cxy≤tanh⁡∣Jxy∣C_{xy}\le\tanh|J_{xy}|。あとは 定理 7.3 です。系 7.4 の条件は osc⁡Φ{x,y}=2∣Jxy∣\operatorname{osc}\Phi_{\{x,y\}}=2|J_{xy}| から sup⁡x∑y2∣Jxy∣<1/2\sup_x\sum_y2|J_{xy}|<1/2、すなわち ∑y∣Jxy∣<1/4\sum_y|J_{xy}|<1/4 で、tanh⁡a≤a\tanh a\le a なのでこれは今の条件 ∑ytanh⁡∣Jxy∣<1\sum_y\tanh|J_{xy}|<1 より強い(狭い)条件です。

  • H.-O. Georgii, Gibbs Measures and Phase Transitions, 2nd ed., de Gruyter, 2011 — 第 1・2 章(仕様と Gibbs 測度)、第 4 章(存在)、第 7 章(極値分解)、第 8 章(Dobrushin の一意性)。
  • S. Friedli, Y. Velenik, Statistical Mechanics of Lattice Systems: A Concrete Mathematical Introduction, Cambridge University Press, 2017 — 第 6 章。
  • R. L. Dobrushin, “The description of a random field by means of conditional probabilities and conditions of its regularity”, Theory of Probability and its Applications 13 (1968), 197–224.
  • O. E. Lanford, D. Ruelle, “Observables at infinity and states with short range correlations in statistical mechanics”, Communications in Mathematical Physics 13 (1969), 194–215.
  • B. Simon, The Statistical Mechanics of Lattice Gases, Vol. I, Princeton University Press, 1993 — Gibbs 状態と DLR 方程式を扱う章。

このテーマの 100 本ノック(50 問)

執筆: Claude Code(Opus 5.5)、Claude Code(Sonnet 5.5)・監査 3 回(Claude Code)・最終 2026-10-12

この記事の誤りを報告する ・運営: 夢現技研合同会社 ・料金 ・利用条件 ・特定商取引法に基づく表記

© 2026 夢現技研合同会社 ・コード例は MIT ライセンスです。