磁石の中には、ほぼ無限に多い小さな磁石(スピン)が並んでいます。有限個ならエネルギーから確率を決める公式がありますが、無限個ではその公式がそのままでは使えません。この記事では、「ある範囲の外を決めたとき、範囲の中はどんな確率で決まるか」という局所的な規則だけから、無限個の系の確率を定義します。読み終えると、その確率がいつ存在するか、いつただ一つに決まるか(決まらないときが相転移です)を、証明つきで説明できるようになります。
無限個のスピンの同時分布は、Boltzmann 因子 e − H / Z e^{-H}/Z e − H / Z では直接書けません(H H H も Z Z Z も発散する)。代わりに「有限領域 Λ \Lambda Λ の外を固定したときの Λ \Lambda Λ 内の条件付き分布」を指定します。この指定の整合的な族が仕様 γ = ( γ Λ ) \gamma=(\gamma_\Lambda) γ = ( γ Λ ) で、それを条件付き分布にもつ確率測度が Gibbs 測度 です(DLR 条件 μ γ Λ = μ \mu\gamma_\Lambda=\mu μ γ Λ = μ )。
存在 (定理 5.1 ): 仕様が連続関数を連続関数に写す(準局所的)なら、有限体積の分布の弱収束極限は Gibbs 測度です。スピンが有限集合に値をとれば配位空間がコンパクトで緊密性は自動的に成り立ち、Gibbs 測度は必ず存在します。
凸構造 (定理 6.2 ): Gibbs 測度全体 G ( γ ) \mathcal{G}(\gamma) G ( γ ) は凸でコンパクトです。その端点(極値 Gibbs 測度。混ぜ合わせで作れない測度のこと)は、ちょうど尾部 σ-加法族(有限個のスピンを変えても変わらない事象の全体)の上で自明 (確率が 0 0 0 か 1 1 1 )なものです。異なる極値 Gibbs 測度は互いに特異 (確率 1 を持つ事象が互いに交わらない)で、物理的には別々の「純粋な相」です。
Dobrushin の一意性定理 (定理 7.3 ): 1 点の条件付き分布が、遠くのスピンを 1 個変えたときに全変動距離でどれだけ動くかを並べた行列 C C C の行和が 1 1 1 未満なら、Gibbs 測度はただ一つです。Ising 模型では 2 d tanh ( β ∣ J ∣ ) < 1 2d\tanh(\beta|J|)<1 2 d tanh ( β ∣ J ∣ ) < 1 (外場は任意)で一意性が出ます。
相転移は「∣ G ( γ ) ∣ ≥ 2 |\mathcal{G}(\gamma)|\ge2 ∣ G ( γ ) ∣ ≥ 2 」として定義でき、次章 の Peierls の議論はこれを低温で示します。
この記事で使う事前知識
条件付き期待値 : 一部の情報(σ-加法族)を与えたとき、残りを平均して得られる確率変数。→ 条件付き期待値
弱収束 : 連続関数の積分がすべて収束する、という確率測度の収束。→ 確率測度の弱収束
緊密性と Prokhorov の定理 : 確率測度の族が緊密なら、弱収束する部分列がとれる。→ 確率測度の弱収束
Radon–Nikodym の密度 : ν ≪ μ \nu\ll\mu ν ≪ μ のとき ν \nu ν を μ \mu μ の重みつき測度として書く関数。→ 条件付き期待値
相互作用 Φ A \Phi_A Φ A : 有限集合 A A A ごとに、A A A 上のスピンだけに依存するエネルギーを与える関数。→ 相互作用
前章 では、圧力 P ( Φ ) = lim ∣ Λ ∣ − 1 log Z Λ P(\Phi)=\lim|\Lambda|^{-1}\log Z_\Lambda P ( Φ ) = lim ∣Λ ∣ − 1 log Z Λ が存在し、境界条件によらないことを示しました。圧力は熱力学量ですが、相転移の本質は状態 、すなわち無限系の確率測度の側にあります。たとえば低温の強磁性体では、境界をすべて + + + にして体積を大きくしたときと、すべて − - − にしたときで、原点のスピンの期待値が異なる値に収束します。境界の影響が無限遠まで残るのです。
この現象を記述するには、無限体積の確率測度そのものを扱う必要があります。Dobrushin(1968)と Lanford–Ruelle(1969)は、Boltzmann 分布の代わりに条件付き分布 を指定するという考えでこれを定式化しました(DLR 条件)。これは確率過程を有限次元分布で指定する Kolmogorov の方法(Kolmogorov の拡張定理 )と似ていますが、決定的な違いがあります。有限次元分布の整合的な族は確率測度を一意に 決めますが、条件付き分布の整合的な族は確率測度を一意に決めるとは限りません。この非一意性が相転移です。
前章と同じく S S S は有限集合で離散位相を入れ、Ω = S Z d \Omega=S^{\mathbb{Z}^d} Ω = S Z d に直積位相を入れます。Λ ⊂ Z d \Lambda\subset\mathbb{Z}^d Λ ⊂ Z d について F Λ \mathcal{F}_\Lambda F Λ を座標 ω x \omega_x ω x (x ∈ Λ x\in\Lambda x ∈ Λ )が生成する σ-加法族、F = F Z d \mathcal{F}=\mathcal{F}_{\mathbb{Z}^d} F = F Z d 、T Λ = F Z d ∖ Λ \mathcal{T}_\Lambda=\mathcal{F}_{\mathbb{Z}^d\setminus\Lambda} T Λ = F Z d ∖ Λ とし、
T = ⋂ Λ ⋐ Z d T Λ \mathcal{T}=\bigcap_{\Lambda\Subset\mathbb{Z}^d}\mathcal{T}_\Lambda T = Λ ⋐ Z d ⋂ T Λ
を尾部 σ-加法族 と呼びます。尾部事象は「有限個のスピンを変えても変わらない事象」です。たとえば { lim n ∣ Λ n ∣ − 1 ∑ x ∈ Λ n ω x \{\lim_n|\Lambda_n|^{-1}\sum_{x\in\Lambda_n}\omega_x { lim n ∣ Λ n ∣ − 1 ∑ x ∈ Λ n ω x が存在して正} \} } (S ⊂ R S\subset\mathbb{R} S ⊂ R )は尾部事象です。
F Λ \mathcal{F}_\Lambda F Λ -可測な関数(Λ \Lambda Λ は有限)を局所関数 と呼びます。
補題 2.1 (配位空間の性質 ) ✓
(i) ∥ x ∥ = max i ∣ x i ∣ \|x\|=\max_i|x_i| ∥ x ∥ = max i ∣ x i ∣ とし、ω ≠ η \omega\ne\eta ω = η に対して d ( ω , η ) = 2 − k d(\omega,\eta)=2^{-k} d ( ω , η ) = 2 − k 、k = min { ∥ x ∥ : ω x ≠ η x } k=\min\{\|x\|:\omega_x\ne\eta_x\} k = min { ∥ x ∥ : ω x = η x } とおくと、d d d は直積位相を定める距離で、( Ω , d ) (\Omega,d) ( Ω , d ) はコンパクトである。
(ii) F \mathcal{F} F は Ω \Omega Ω のボレル σ-加法族に等しい。
(iii) 局所関数は連続で、C ( Ω ) C(\Omega) C ( Ω ) の中で一様ノルムについて稠密である。
証明(補題 2.1) (i) d ( ω , η ) < 2 − k d(\omega,\eta)<2^{-k} d ( ω , η ) < 2 − k は「∥ x ∥ ≤ k \|x\|\le k ∥ x ∥ ≤ k で ω x = η x \omega_x=\eta_x ω x = η x 」と同値なので、d d d の開球は箱 B k = { ∥ x ∥ ≤ k } B_k=\{\|x\|\le k\} B k = { ∥ x ∥ ≤ k } 上の値を固定した集合(柱集合)で、これは直積位相の基本開集合です。逆に有限個の座標を固定した集合は、それを含む B k B_k B k の上の値を固定した柱集合の有限和です。三角不等式は強い形 d ( ω , ζ ) ≤ max ( d ( ω , η ) , d ( η , ζ ) ) d(\omega,\zeta)\le\max(d(\omega,\eta),d(\eta,\zeta)) d ( ω , ζ ) ≤ max ( d ( ω , η ) , d ( η , ζ )) で成り立ちます。コンパクト性は点列で示します。( ω n ) (\omega^n) ( ω n ) に対し、S S S は有限なので ω B 1 n \omega^n_{B_1} ω B 1 n が一定になる部分列をとり、その中で ω B 2 n \omega^n_{B_2} ω B 2 n が一定になる部分列をとり……と続け、対角線の部分列をとれば、それは各 B k B_k B k 上で最終的に一定で、収束します。
(ii) 柱集合は可算個しかなく、d d d の開集合は開球(柱集合)の和集合なので、可算個の柱集合の和です。よって開集合は F \mathcal{F} F に属します。逆に座標写像 ω ↦ ω x \omega\mapsto\omega_x ω ↦ ω x は連続なので F \mathcal{F} F はボレル集合族に含まれます。
(iii) 局所関数は有限個の座標の関数で、座標は連続なので連続です。f ∈ C ( Ω ) f\in C(\Omega) f ∈ C ( Ω ) と ε > 0 \varepsilon>0 ε > 0 について、一様連続性(補題 2.2[Arzelà–Ascoli と Stone–Weierstrass] )から k k k を、d ( ω , η ) < 2 − k d(\omega,\eta)<2^{-k} d ( ω , η ) < 2 − k なら ∣ f ( ω ) − f ( η ) ∣ < ε |f(\omega)-f(\eta)|<\varepsilon ∣ f ( ω ) − f ( η ) ∣ < ε となるようにとります。η 0 ∈ Ω \eta^0\in\Omega η 0 ∈ Ω を固定して f k ( ω ) = f ( ω B k η B k c 0 ) f_k(\omega)=f(\omega_{B_k}\eta^0_{B_k^c}) f k ( ω ) = f ( ω B k η B k c 0 ) とおくと、f k f_k f k は局所関数で ∥ f − f k ∥ ∞ ≤ ε \|f-f_k\|_\infty\le\varepsilon ∥ f − f k ∥ ∞ ≤ ε です。
∎
確率測度 μ \mu μ は ( Ω , F ) (\Omega,\mathcal{F}) ( Ω , F ) 上で考え、μ ( f ) = ∫ f d μ \mu(f)=\int f\,d\mu μ ( f ) = ∫ f d μ と書きます。局所関数(あるいは連続関数)の積分が一致する 2 つの確率測度は等しいので(補題 2.2[確率測度の弱収束] と (iii))、確率測度は局所関数の期待値で決まります。
この節のねらいは、「外を固定したときの内側の確率」を一般的な形で扱うための言葉をそろえることです。たとえば Ising 模型で 1 点 x x x だけを再標本化するなら、隣のスピンの値(外)を与えると、x x x のスピンが + 1 +1 + 1 になる確率は式 ( 1 + tanh β m ) / 2 (1+\tanh\beta m)/2 ( 1 + tanh β m ) /2 (m m m は隣のスピンと外場から決まる局所場)で決まります(第 7 節で導きます)。この「外の配位から内側の確率分布を返す対応」が確率核で、以下の定義はそれを一般化したものです。
Ω \Omega Ω 上の確率核 π \pi π を G ⊂ F \mathcal{G}\subset\mathcal{F} G ⊂ F から F \mathcal{F} F への写像 ( B , ω ) ↦ π ( B ∣ ω ) (B,\omega)\mapsto\pi(B\mid\omega) ( B , ω ) ↦ π ( B ∣ ω ) で、各 ω \omega ω について π ( ⋅ ∣ ω ) \pi(\cdot\mid\omega) π ( ⋅ ∣ ω ) が確率測度、各 B ∈ F B\in\mathcal{F} B ∈ F について ω ↦ π ( B ∣ ω ) \omega\mapsto\pi(B\mid\omega) ω ↦ π ( B ∣ ω ) が G \mathcal{G} G -可測なものとします。有界可測な f f f に π f ( ω ) = ∫ f ( η ) π ( d η ∣ ω ) \pi f(\omega)=\int f(\eta)\,\pi(d\eta\mid\omega) π f ( ω ) = ∫ f ( η ) π ( d η ∣ ω ) 、確率測度 μ \mu μ に μ π ( B ) = ∫ π ( B ∣ ω ) μ ( d ω ) \mu\pi(B)=\int\pi(B\mid\omega)\,\mu(d\omega) μ π ( B ) = ∫ π ( B ∣ ω ) μ ( d ω ) を対応させます。μ π ( f ) = μ ( π f ) \mu\pi(f)=\mu(\pi f) μ π ( f ) = μ ( π f ) は、f f f が指示関数なら定義そのもので、単関数を経て有界可測関数に広がります。核の合成は ( π π ′ ) ( B ∣ ω ) = ∫ π ′ ( B ∣ η ) π ( d η ∣ ω ) (\pi\pi')(B\mid\omega)=\int\pi'(B\mid\eta)\,\pi(d\eta\mid\omega) ( π π ′ ) ( B ∣ ω ) = ∫ π ′ ( B ∣ η ) π ( d η ∣ ω ) です。
定義 3.1 (仕様 ) ↗
Λ ⋐ Z d \Lambda\Subset\mathbb{Z}^d Λ ⋐ Z d ごとに T Λ \mathcal{T}_\Lambda T Λ から F \mathcal{F} F への確率核 γ Λ \gamma_\Lambda γ Λ が与えられ、次を満たすとき、族 γ = ( γ Λ ) Λ ⋐ Z d \gamma=(\gamma_\Lambda)_{\Lambda\Subset\mathbb{Z}^d} γ = ( γ Λ ) Λ ⋐ Z d を仕様 という。
(固有性 )B ∈ T Λ B\in\mathcal{T}_\Lambda B ∈ T Λ ならば γ Λ ( B ∣ ω ) = 1 B ( ω ) \gamma_\Lambda(B\mid\omega)=\mathbf{1}_B(\omega) γ Λ ( B ∣ ω ) = 1 B ( ω ) 。
(整合性 )Λ ⊂ Δ \Lambda\subset\Delta Λ ⊂ Δ ならば γ Δ γ Λ = γ Δ \gamma_\Delta\gamma_\Lambda=\gamma_\Delta γ Δ γ Λ = γ Δ 。
さらに、すべての Λ \Lambda Λ と f ∈ C ( Ω ) f\in C(\Omega) f ∈ C ( Ω ) について γ Λ f ∈ C ( Ω ) \gamma_\Lambda f\in C(\Omega) γ Λ f ∈ C ( Ω ) のとき、γ \gamma γ は準局所的 であるという。
固有性は「Λ \Lambda Λ の外は固定されたまま」を、整合性は「大きな領域で再標本化した後に小さな領域で再標本化しても分布は変わらない」を意味します。
補題 3.2 (固有性からの取り出し ) ✓
仕様 γ \gamma γ と有界 T Λ \mathcal{T}_\Lambda T Λ -可測な g g g 、有界可測な f f f について γ Λ ( g f ) = g γ Λ f \gamma_\Lambda(gf)=g\,\gamma_\Lambda f γ Λ ( g f ) = g γ Λ f 。
証明(補題 3.2) g = 1 B g=\mathbf{1}_B g = 1 B (B ∈ T Λ B\in\mathcal{T}_\Lambda B ∈ T Λ )のとき、ω ∈ B \omega\in B ω ∈ B なら γ Λ ( B c ∣ ω ) = 0 \gamma_\Lambda(B^c\mid\omega)=0 γ Λ ( B c ∣ ω ) = 0 なので γ Λ ( 1 B f ) ( ω ) = γ Λ f ( ω ) \gamma_\Lambda(\mathbf{1}_Bf)(\omega)=\gamma_\Lambda f(\omega) γ Λ ( 1 B f ) ( ω ) = γ Λ f ( ω ) 、ω ∉ B \omega\notin B ω ∈ / B なら γ Λ ( B ∣ ω ) = 0 \gamma_\Lambda(B\mid\omega)=0 γ Λ ( B ∣ ω ) = 0 なので γ Λ ( 1 B f ) ( ω ) = 0 \gamma_\Lambda(\mathbf{1}_Bf)(\omega)=0 γ Λ ( 1 B f ) ( ω ) = 0 です。線形性で単関数に、一様近似(∣ γ Λ h ∣ ≤ ∥ h ∥ ∞ |\gamma_\Lambda h|\le\|h\|_\infty ∣ γ Λ h ∣ ≤ ∥ h ∥ ∞ )で有界可測関数に広がります。
∎
物理で使うのは相互作用から作る仕様です。前章と違い、ここでは並進不変性を仮定しません。
定義 3.3 (Gibbs 仕様 ) ↗
相互作用 Φ = ( Φ A ) \Phi=(\Phi_A) Φ = ( Φ A ) (定義 2.1[熱力学的極限] )が sup x ∑ A ∋ x ∥ Φ A ∥ ∞ < ∞ \sup_{x}\sum_{A\ni x}\|\Phi_A\|_\infty<\infty sup x ∑ A ∋ x ∥ Φ A ∥ ∞ < ∞ を満たすとする。Λ ⋐ Z d \Lambda\Subset\mathbb{Z}^d Λ ⋐ 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})}, H Λ ( ω ) = A ∩ Λ = ∅ ∑ Φ A ( ω ) , Z Λ ω = ξ ∈ Ω Λ ∑ e − H Λ ( ξ ω Λ c ) , γ Λ Φ ( f ∣ ω ) = 1 Z Λ ω ∑ ξ ∈ Ω Λ 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}) γ Λ Φ ( f ∣ ω ) = Z Λ ω 1 ξ ∈ Ω Λ ∑ e − H Λ ( ξ ω Λ c ) f ( ξ ω Λ c ) とおく。γ Φ = ( γ Λ Φ ) \gamma^\Phi=(\gamma^\Phi_\Lambda) γ Φ = ( γ Λ Φ ) を Φ \Phi Φ の Gibbs 仕様 という。
H Λ H_\Lambda H Λ の級数は ∑ 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 ∑ A ∩ Λ = ∅ ∥ Φ A ∥ ∞ ≤ ∣Λ∣ sup x ∑ A ∋ x ∥ Φ A ∥ ∞ で絶対収束します。境界条件 ω Λ c \omega_{\Lambda^c} ω Λ c との相互作用も含めたエネルギーであることに注意してください(前章の H Λ τ H_\Lambda^\tau H Λ τ )。
命題 3.4 (Gibbs 仕様は準局所的な仕様 ) ✓
γ Φ \gamma^\Phi γ Φ は準局所的な仕様である。
証明(命題 3.4) 核であること・固有性. γ Λ Φ ( ⋅ ∣ ω ) \gamma^\Phi_\Lambda(\cdot\mid\omega) γ Λ Φ ( ⋅ ∣ ω ) は有限個の点 ξ ω Λ c \xi\omega_{\Lambda^c} ξ ω Λ c に重みを置いた確率測度で、重みは ω Λ c \omega_{\Lambda^c} ω Λ c だけに依存するので T Λ \mathcal{T}_\Lambda T Λ -可測です。B ∈ T Λ B\in\mathcal{T}_\Lambda B ∈ T Λ なら 1 B ( ξ ω Λ c ) = 1 B ( ω ) \mathbf{1}_B(\xi\omega_{\Lambda^c})=\mathbf{1}_B(\omega) 1 B ( ξ ω Λ c ) = 1 B ( ω ) なので固有性が成り立ちます。
整合性. Λ ⊂ Δ \Lambda\subset\Delta Λ ⊂ Δ とし、Δ \Delta Δ 上の配位を ( α , ζ ) (\alpha,\zeta) ( α , ζ ) (α ∈ Ω Λ \alpha\in\Omega_\Lambda α ∈ Ω Λ 、ζ ∈ Ω Δ ∖ Λ \zeta\in\Omega_{\Delta\setminus\Lambda} ζ ∈ Ω Δ ∖ Λ )と書きます。H Δ − H Λ = ∑ Φ A H_\Delta-H_\Lambda=\sum\Phi_A H Δ − H Λ = ∑ Φ A (和は A ∩ Δ ≠ ∅ A\cap\Delta\ne\emptyset A ∩ Δ = ∅ 、A ∩ Λ = ∅ A\cap\Lambda=\emptyset A ∩ Λ = ∅ )は α \alpha α に依存しないので、これを K ( ζ ω Δ c ) K(\zeta\omega_{\Delta^c}) K ( ζ ω Δ 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})}. e − H Δ ( α ζ ω Δ c ) = e − K ( ζ ω Δ c ) e − H Λ ( α ζ ω Δ c ) . g = γ Λ Φ f g=\gamma^\Phi_\Lambda f g = γ Λ Φ f は α \alpha α に依存しないので
γ Δ Φ g ( ω ) = 1 Z Δ ω ∑ ζ e − K ( ∑ α e − H Λ ( α ζ ω Δ c ) ) 1 Z Λ ζ ω Δ c ∑ α ′ e − H Λ ( α ′ ζ ω Δ c ) f ( α ′ ζ ω Δ c ) = 1 Z Δ ω ∑ ζ ∑ α ′ e − K e − 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} γ Δ Φ g ( ω ) = Z Δ ω 1 ζ ∑ e − K ( α ∑ e − H Λ ( α ζ ω Δ c ) ) Z Λ ζ ω Δ c 1 α ′ ∑ e − H Λ ( α ′ ζ ω Δ c ) f ( α ′ ζ ω Δ c ) = Z Δ ω 1 ζ ∑ α ′ ∑ e − K e − H Λ ( α ′ ζ ω Δ c ) f ( α ′ ζ ω Δ c ) = γ Δ Φ f ( ω ) . 2 行目では括弧内の和が Z Λ ζ ω Δ c Z_\Lambda^{\zeta\omega_{\Delta^c}} Z Λ ζ ω Δ c に等しいことを使いました。
準局所性. ω ↦ H Λ ( ξ ω Λ c ) \omega\mapsto H_\Lambda(\xi\omega_{\Lambda^c}) ω ↦ H Λ ( ξ ω Λ c ) は連続関数(局所関数)の一様収束級数なので連続です。f f f が連続なら γ Λ Φ f \gamma^\Phi_\Lambda f γ Λ Φ f は連続関数の有限和を正の連続関数で割ったものなので連続です。
∎
例 3.5 (相互作用のない場合 )
Φ = 0 \Phi=0 Φ = 0 なら γ Λ ( ⋅ ∣ ω ) \gamma_\Lambda(\cdot\mid\omega) γ Λ ( ⋅ ∣ ω ) は「Λ \Lambda Λ 上は一様分布、外は ω \omega ω のまま」です。後で見るように、この仕様の Gibbs 測度は一様分布の無限直積だけです(演習 8.1 )。
定義 4.1 (Gibbs 測度(DLR 条件) ) ↗
仕様 γ \gamma γ について、Ω \Omega Ω 上の確率測度 μ \mu μ がすべての Λ ⋐ Z d \Lambda\Subset\mathbb{Z}^d Λ ⋐ Z d で
μ γ Λ = μ \mu\gamma_\Lambda=\mu μ γ Λ = μ を満たすとき、μ \mu μ を γ \gamma γ の Gibbs 測度 という。これを DLR 方程式 (Dobrushin–Lanford–Ruelle)と呼び、Gibbs 測度全体を G ( γ ) \mathcal{G}(\gamma) G ( γ ) と書く。γ = γ Φ \gamma=\gamma^\Phi γ = γ Φ のときは G ( Φ ) \mathcal{G}(\Phi) G ( Φ ) とも書く。
命題 4.2 (DLR 方程式と条件付き期待値 ) ✓
確率測度 μ \mu μ と Λ ⋐ Z d \Lambda\Subset\mathbb{Z}^d Λ ⋐ Z d について、次は同値である。
(a) μ γ Λ = μ \mu\gamma_\Lambda=\mu μ γ Λ = μ 。
(b) すべての有界可測な f f f について E μ [ f ∣ T Λ ] = γ Λ f E_\mu[f\mid\mathcal{T}_\Lambda]=\gamma_\Lambda f E μ [ f ∣ T Λ ] = γ Λ f (μ \mu μ -a.s.)。
証明(命題 4.2) (a)⇒(b): γ Λ f \gamma_\Lambda f γ Λ f は T Λ \mathcal{T}_\Lambda T Λ -可測で有界です。B ∈ T Λ B\in\mathcal{T}_\Lambda B ∈ T Λ について、補題 3.2 と (a) から
μ ( 1 B γ Λ f ) = μ ( γ Λ ( 1 B f ) ) = μ γ Λ ( 1 B f ) = μ ( 1 B f ) . \mu(\mathbf{1}_B\gamma_\Lambda f)=\mu(\gamma_\Lambda(\mathbf{1}_Bf))=\mu\gamma_\Lambda(\mathbf{1}_Bf)=\mu(\mathbf{1}_Bf). μ ( 1 B γ Λ f ) = μ ( γ Λ ( 1 B f )) = μ γ Λ ( 1 B f ) = μ ( 1 B f ) . これは条件付き期待値の定義(定義 3.1[条件付き期待値] )の条件そのものです。
(b)⇒(a): B = Ω B=\Omega B = Ω とすれば μ ( γ Λ f ) = μ ( f ) \mu(\gamma_\Lambda f)=\mu(f) μ ( γ Λ f ) = μ ( f ) 、すなわち μ γ Λ ( f ) = μ ( f ) \mu\gamma_\Lambda(f)=\mu(f) μ γ Λ ( f ) = μ ( f ) です。
∎
Gibbs 仕様では、(b) は「Λ \Lambda Λ の外の配位 ω Λ c \omega_{\Lambda^c} ω Λ c を与えたときの ω Λ \omega_\Lambda ω Λ の条件付き分布が e − H Λ ( ξ ω Λ c ) / Z Λ ω e^{-H_\Lambda(\xi\omega_{\Lambda^c})}/Z_\Lambda^\omega e − H Λ ( ξ ω Λ c ) / Z Λ ω である」という意味です。条件付き期待値 は μ \mu μ -a.s. にしか決まらないので、「すべての境界条件 ω \omega ω で」条件付き分布を指定する仕様は、測度より多くの情報を含んでいます。
定理 5.1 (Gibbs 測度の存在 ) ✓
γ \gamma γ を準局所的な仕様とし、Λ 1 ⊂ Λ 2 ⊂ ⋯ \Lambda_1\subset\Lambda_2\subset\cdots Λ 1 ⊂ Λ 2 ⊂ ⋯ を ⋃ n Λ n = Z d \bigcup_n\Lambda_n=\mathbb{Z}^d ⋃ n Λ n = Z d となる有限集合の列、ω ( n ) ∈ Ω \omega^{(n)}\in\Omega ω ( n ) ∈ Ω を任意の境界条件とする。μ n = γ Λ n ( ⋅ ∣ ω ( n ) ) \mu_n=\gamma_{\Lambda_n}(\cdot\mid\omega^{(n)}) μ n = γ Λ n ( ⋅ ∣ ω ( n ) ) とおく。
(i) ( μ n ) (\mu_n) ( μ n ) は弱収束する部分列をもつ。
(ii) ( μ n ) (\mu_n) ( μ n ) の部分列の弱収束極限はすべて G ( γ ) \mathcal{G}(\gamma) G ( γ ) に属する。
特に G ( γ ) ≠ ∅ \mathcal{G}(\gamma)\ne\emptyset G ( γ ) = ∅ である。
証明(定理 5.1) (i) Ω \Omega Ω はコンパクト距離空間(補題 2.1 )なので、確率測度の族はすべて緊密です(定義 3.1[確率測度の弱収束] で K = Ω K=\Omega K = Ω )。Prokhorov の定理(定理 3.3[確率測度の弱収束] )から弱収束する部分列があります。
(ii) μ n k ⇒ μ \mu_{n_k}\Rightarrow\mu μ n k ⇒ μ とし、Λ ⋐ Z d \Lambda\Subset\mathbb{Z}^d Λ ⋐ Z d を固定します。Λ \Lambda Λ は有限なので、ある k 0 k_0 k 0 から先 Λ ⊂ Λ n k \Lambda\subset\Lambda_{n_k} Λ ⊂ Λ n k です。そのとき整合性から
μ n k γ Λ ( f ) = ∫ γ Λ f d γ Λ n k ( ⋅ ∣ ω ( n k ) ) = ( γ Λ n k γ Λ f ) ( ω ( n k ) ) = ( γ Λ n k f ) ( ω ( n k ) ) = μ n k ( 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). μ n k γ Λ ( f ) = ∫ γ Λ f d γ Λ n k ( ⋅ ∣ ω ( n k ) ) = ( γ Λ n k γ Λ f ) ( ω ( n k ) ) = ( γ Λ n k f ) ( ω ( n k ) ) = μ n k ( f ) . f ∈ C ( Ω ) f\in C(\Omega) f ∈ C ( Ω ) なら準局所性から γ Λ f ∈ C ( Ω ) \gamma_\Lambda f\in C(\Omega) γ Λ f ∈ C ( Ω ) なので、弱収束の定義(定義 2.1[確率測度の弱収束] )から
μ ( γ Λ f ) = lim k μ n k ( γ Λ f ) = lim k μ n k ( f ) = μ ( f ) . \mu(\gamma_\Lambda f)=\lim_k\mu_{n_k}(\gamma_\Lambda f)=\lim_k\mu_{n_k}(f)=\mu(f). μ ( γ Λ f ) = k lim μ n k ( γ Λ f ) = k lim μ n k ( f ) = μ ( f ) . 2 つの確率測度 μ γ Λ \mu\gamma_\Lambda μ γ Λ と μ \mu μ がすべての連続関数で一致するので(補題 2.2[確率測度の弱収束] )μ γ Λ = μ \mu\gamma_\Lambda=\mu μ γ Λ = μ です。
∎
この証明でコンパクト性を使ったのは (i) の緊密性だけです。スピンが R \mathbb{R} R に値をとる格子場(たとえば格子 ϕ 4 \phi^4 ϕ 4 模型)では Ω = R Z d \Omega=\mathbb{R}^{\mathbb{Z}^d} Ω = R Z d はコンパクトでなく、境界条件と相互作用について一様な積率評価 を示して緊密性を確保する必要があります(Ruelle の超安定性評価など)。構成的場の理論で無限体積極限を作る議論(相関不等式 の単調性による方法など)は、この緊密性を別の手段で得るものと見ることができます。
系 5.2 (Gibbs 測度の集合は凸でコンパクト ) ✓
準局所的な仕様 γ \gamma γ について、G ( γ ) \mathcal{G}(\gamma) G ( γ ) は凸で、弱収束について閉じている(したがって点列コンパクトである)。
証明(系 5.2) μ γ Λ \mu\gamma_\Lambda μ γ Λ は μ \mu μ について線形なので、DLR 方程式は凸結合で保たれます。μ n ∈ G ( γ ) \mu_n\in\mathcal{G}(\gamma) μ n ∈ G ( γ ) 、μ n ⇒ μ \mu_n\Rightarrow\mu μ n ⇒ μ なら、定理 5.1 (ii) の証明と同じく、f ∈ C ( Ω ) f\in C(\Omega) f ∈ C ( Ω ) について μ ( γ Λ f ) = lim μ n ( γ Λ f ) = lim μ n ( f ) = μ ( f ) \mu(\gamma_\Lambda f)=\lim\mu_n(\gamma_\Lambda f)=\lim\mu_n(f)=\mu(f) μ ( γ Λ f ) = lim μ n ( γ Λ f ) = lim μ n ( f ) = μ ( f ) で、μ ∈ G ( γ ) \mu\in\mathcal{G}(\gamma) μ ∈ G ( γ ) 。点列コンパクト性は Prokhorov の定理から従います。
∎
この節のねらいは、「純粋な相」(混ぜ物でない状態)を数学の言葉で捉えることです。直観では、+ + + 相と − - − 相を半分ずつ混ぜた状態は、どちらの相にいるかという「無限遠の情報」について確率 1 / 2 1/2 1/2 の不確かさを持ちます。つまり混合であることと、尾部事象の確率が 0 0 0 か 1 1 1 でないことは同じ現象で、これを次の補題と定理で正確にします。
凸集合 G ( γ ) \mathcal{G}(\gamma) G ( γ ) の元 μ \mu μ が極値的 とは、μ = t ν 1 + ( 1 − t ) ν 2 \mu=t\nu_1+(1-t)\nu_2 μ = t ν 1 + ( 1 − t ) ν 2 (0 < t < 1 0<t<1 0 < t < 1 、ν i ∈ G ( γ ) \nu_i\in\mathcal{G}(\gamma) ν i ∈ G ( γ ) )ならば ν 1 = ν 2 = μ \nu_1=\nu_2=\mu ν 1 = ν 2 = μ となることです。極値 Gibbs 測度全体を ex G ( γ ) \operatorname{ex}\mathcal{G}(\gamma) ex G ( γ ) と書きます。
次の補題は、極値性と尾部自明性を結ぶ道具です。ある Gibbs 測度が別の Gibbs 測度を「含む」とき、その重みの関数は有限個のスピンの値に左右されない、つまり尾部の情報だけで決まる、と述べています。
補題 6.1 (Gibbs 測度どうしの密度は尾部可測 ) ✓
μ , ν ∈ G ( γ ) \mu,\nu\in\mathcal{G}(\gamma) μ , ν ∈ G ( γ ) で ν ≪ μ \nu\ll\mu ν ≪ μ ならば、密度 d ν / d μ d\nu/d\mu d ν / d μ は μ \mu μ -a.s. で T \mathcal{T} T -可測な関数に等しい。
証明(補題 6.1) h = d ν / d μ h=d\nu/d\mu h = d ν / d μ (Radon–Nikodym の定理、定理 2.1[条件付き期待値] )とし、Λ ⋐ Z d \Lambda\Subset\mathbb{Z}^d Λ ⋐ Z d を固定して h Λ = E μ [ h ∣ T Λ ] h_\Lambda=E_\mu[h\mid\mathcal{T}_\Lambda] h Λ = E μ [ h ∣ T Λ ] とおきます。有界可測な f f f について、ν , μ ∈ G ( γ ) \nu,\mu\in\mathcal{G}(\gamma) ν , μ ∈ G ( γ ) 、γ Λ f \gamma_\Lambda f γ Λ f の T Λ \mathcal{T}_\Lambda T Λ -可測性、補題 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). ν ( f ) = ν ( γ Λ f ) = μ ( h γ Λ f ) = μ ( h Λ γ Λ f ) = μ ( γ Λ ( h Λ f ) ) = μ ( h Λ f ) . 3 番目の等号は既知量の取り出し(定理 4.4[条件付き期待値] 、γ Λ f \gamma_\Lambda f γ Λ f は有界 T Λ \mathcal{T}_\Lambda T Λ -可測)、最後は μ γ Λ = μ \mu\gamma_\Lambda=\mu μ γ Λ = μ です。4 番目と最後の等号では h Λ f h_\Lambda f h Λ f が有界とは限りませんが、f ≥ 0 f\ge0 f ≥ 0 としておけば、補題 3.2 と μ ( γ Λ g ) = μ ( g ) \mu(\gamma_\Lambda g)=\mu(g) μ ( γ Λ g ) = μ ( g ) は g = ( h Λ ∧ M ) f g=(h_\Lambda\wedge M)f g = ( h Λ ∧ M ) f に対して成り立ち、M → ∞ M\to\infty M → ∞ で単調収束定理(定理 2.1[確率変数と期待値] )により g = h Λ f g=h_\Lambda f g = h Λ f に広がります(一般の有界 f f f は正負に分けます)。よって h Λ h_\Lambda h Λ も ν \nu ν の μ \mu μ に関する密度で、密度の一意性から h = h Λ h=h_\Lambda h = h Λ (μ \mu μ -a.s.)です。
Λ n ↑ Z d \Lambda_n\uparrow\mathbb{Z}^d Λ n ↑ Z d をとり、各 h Λ n h_{\Lambda_n} h Λ n の T Λ n \mathcal{T}_{\Lambda_n} T Λ n -可測な版を固定して h ~ = lim sup n h Λ n \tilde h=\limsup_nh_{\Lambda_n} h ~ = lim sup n h Λ n とおきます。T Λ n \mathcal{T}_{\Lambda_n} T Λ n は n n n について減少するので、m m m を固定すると n ≥ m n\ge m n ≥ m の h Λ n h_{\Lambda_n} h Λ n はすべて T Λ m \mathcal{T}_{\Lambda_m} T Λ m -可測で、h ~ \tilde h h ~ は T Λ m \mathcal{T}_{\Lambda_m} T Λ m -可測です。任意の有限 Λ \Lambda Λ はある Λ m \Lambda_m Λ m に含まれ T Λ m ⊂ T Λ \mathcal{T}_{\Lambda_m}\subset\mathcal{T}_\Lambda T Λ m ⊂ T Λ なので、h ~ \tilde h h ~ は T \mathcal{T} T -可測です。可算個の零集合を除いて h Λ n = h h_{\Lambda_n}=h h Λ n = h なので h ~ = h \tilde h=h h ~ = h (μ \mu μ -a.s.)です。
∎
定理 6.2 (極値 Gibbs 測度の特徴づけ ) ✓
μ ∈ G ( γ ) \mu\in\mathcal{G}(\gamma) μ ∈ G ( γ ) について、次は同値である。
(a) μ \mu μ は G ( γ ) \mathcal{G}(\gamma) G ( γ ) の極値点である。
(b) μ \mu μ は尾部自明である。すなわちすべての B ∈ T B\in\mathcal{T} B ∈ T で μ ( B ) ∈ { 0 , 1 } \mu(B)\in\{0,1\} μ ( B ) ∈ { 0 , 1 } 。
証明(定理 6.2) (b)⇒(a). μ = t ν 1 + ( 1 − t ) ν 2 \mu=t\nu_1+(1-t)\nu_2 μ = t ν 1 + ( 1 − t ) ν 2 (0 < t < 1 0<t<1 0 < t < 1 、ν i ∈ G ( γ ) \nu_i\in\mathcal{G}(\gamma) ν i ∈ G ( γ ) )とします。ν 1 ≤ t − 1 μ \nu_1\le t^{-1}\mu ν 1 ≤ t − 1 μ なので ν 1 ≪ μ \nu_1\ll\mu ν 1 ≪ μ で、補題 6.1 から密度 h h h は T \mathcal{T} T -可測です。尾部自明な測度に関して T \mathcal{T} T -可測な関数は a.s. 定数です({ h > c } ∈ T \{h>c\}\in\mathcal{T} { h > c } ∈ T の確率が 0 0 0 か 1 1 1 なので、c ∗ = sup { c : μ ( h > c ) = 1 } c^*=\sup\{c:\mu(h>c)=1\} c ∗ = sup { c : μ ( h > c ) = 1 } とすれば h = c ∗ h=c^* h = c ∗ a.s.)。μ ( h ) = ν 1 ( Ω ) = 1 \mu(h)=\nu_1(\Omega)=1 μ ( h ) = ν 1 ( Ω ) = 1 から h = 1 h=1 h = 1 、すなわち ν 1 = μ \nu_1=\mu ν 1 = μ 。同様に ν 2 = μ \nu_2=\mu ν 2 = μ です。
(a)⇒(b). B ∈ T B\in\mathcal{T} B ∈ T 、0 < μ ( B ) < 1 0<\mu(B)<1 0 < μ ( B ) < 1 と仮定し、ν = μ ( ⋅ ∣ B ) \nu=\mu(\cdot\mid B) ν = μ ( ⋅ ∣ B ) 、ν ′ = μ ( ⋅ ∣ B c ) \nu'=\mu(\cdot\mid B^c) ν ′ = μ ( ⋅ ∣ B c ) とおきます。B ∈ T Λ B\in\mathcal{T}_\Lambda B ∈ T Λ なので 補題 3.2 から、有界可測な f f f について
ν ( γ Λ f ) = μ ( 1 B γ Λ f ) μ ( B ) = μ ( γ Λ ( 1 B f ) ) μ ( B ) = μ ( 1 B f ) μ ( 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). ν ( γ Λ f ) = μ ( B ) μ ( 1 B γ Λ f ) = μ ( B ) μ ( γ Λ ( 1 B f )) = μ ( B ) μ ( 1 B f ) = ν ( f ) . よって ν ∈ G ( γ ) \nu\in\mathcal{G}(\gamma) ν ∈ G ( γ ) 、同様に ν ′ ∈ G ( γ ) \nu'\in\mathcal{G}(\gamma) ν ′ ∈ G ( γ ) で、μ = μ ( B ) ν + μ ( B c ) ν ′ \mu=\mu(B)\nu+\mu(B^c)\nu' μ = μ ( B ) ν + μ ( B c ) ν ′ です。ν ( B ) = 1 ≠ μ ( B ) \nu(B)=1\ne\mu(B) ν ( B ) = 1 = μ ( B ) なので ν ≠ μ \nu\ne\mu ν = μ となり、μ \mu μ は極値的でありません。
∎
系 6.3 (尾部での値が Gibbs 測度を決める ) ✓
(i) μ , ν ∈ G ( γ ) \mu,\nu\in\mathcal{G}(\gamma) μ , ν ∈ G ( γ ) が T \mathcal{T} T の上で一致すれば μ = ν \mu=\nu μ = ν 。
(ii) μ ≠ ν \mu\ne\nu μ = ν がともに極値 Gibbs 測度ならば、μ ( B ) = 1 \mu(B)=1 μ ( B ) = 1 、ν ( B ) = 0 \nu(B)=0 ν ( B ) = 0 となる B ∈ T B\in\mathcal{T} B ∈ T がある。特に μ \mu μ と ν \nu ν は互いに特異である。
証明(系 6.3) (i) ρ = ( μ + ν ) / 2 ∈ G ( γ ) \rho=(\mu+\nu)/2\in\mathcal{G}(\gamma) ρ = ( μ + ν ) /2 ∈ G ( γ ) について μ , ν ≪ ρ \mu,\nu\ll\rho μ , ν ≪ ρ で、補題 6.1 から密度 h , h ′ h,h' h , h ′ は T \mathcal{T} T -可測にとれます。B ∈ T B\in\mathcal{T} B ∈ T で ρ ( h 1 B ) = μ ( B ) = ν ( B ) = ρ ( h ′ 1 B ) \rho(h\mathbf{1}_B)=\mu(B)=\nu(B)=\rho(h'\mathbf{1}_B) ρ ( h 1 B ) = μ ( B ) = ν ( B ) = ρ ( h ′ 1 B ) なので、B = { h > h ′ } B=\{h>h'\} B = { h > h ′ } 、{ h < h ′ } \{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 ) = ν ( B ) となる B ∈ T B\in\mathcal{T} B ∈ T があり、定理 6.2 から値は 0 0 0 か 1 1 1 です。必要なら B B B を補集合に替えます。
∎
極値 Gibbs 測度は、尾部の情報(「無限遠での境界条件」)が確定した状態であり、物理的には純粋な相 と解釈されます。たとえば低温の Ising 模型の + + + 相と − - − 相は、尾部事象「空間平均磁化が正」で区別されます。
この節のねらいは、Gibbs 測度がただ一つに決まるための、確かめやすい十分条件を与えることです。考え方は、スピンが遠くのスピンから受ける影響の強さを数 C x y C_{xy} C x y で測り、その総和が 1 1 1 未満なら、境界条件の違いが遠くへ行くほど薄まって消える、というものです。
S S S 上の確率測度 p , q p,q p , q の全変動距離 を
∥ p − q ∥ = max B ⊂ S ∣ p ( B ) − q ( B ) ∣ = 1 2 ∑ 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)| ∥ p − q ∥ = B ⊂ S max ∣ p ( B ) − q ( B ) ∣ = 2 1 s ∈ S ∑ ∣ p ( s ) − q ( s ) ∣
とします。2 番目の等号は、最大値が B = { p > q } B=\{p>q\} B = { p > q } でとられることからわかります。仕様 γ \gamma γ の 1 点集合 { x } \{x\} { x } の核から、x x x のスピンの条件付き分布 γ x ( ⋅ ∣ ω ) = γ { x } ( ω x ∈ ⋅ ∣ ω ) \gamma_x(\cdot\mid\omega)=\gamma_{\{x\}}(\omega_x\in\cdot\mid\omega) γ x ( ⋅ ∣ ω ) = γ { x } ( ω x ∈ ⋅ ∣ ω ) が決まります。固有性から γ { x } ( ⋅ ∣ ω ) \gamma_{\{x\}}(\cdot\mid\omega) γ { x } ( ⋅ ∣ ω ) は { η : η x c = ω x c } \{\eta:\eta_{x^c}=\omega_{x^c}\} { η : η x c = ω x c } に集中し、γ x ( ⋅ ∣ ω ) \gamma_x(\cdot\mid\omega) γ x ( ⋅ ∣ ω ) は ω x c \omega_{x^c} ω x c だけに依存します。したがって
γ { x } f ( ω ) = ∑ s ∈ S γ x ( s ∣ ω ) f ( s ω x c ) \gamma_{\{x\}}f(\omega)=\sum_{s\in S}\gamma_x(s\mid\omega)\,f(s\,\omega_{x^c}) γ { x } f ( ω ) = s ∈ S ∑ γ x ( s ∣ ω ) f ( s ω x c )
です(s ω x c s\,\omega_{x^c} s ω x c は x x x で s s s 、他で ω \omega ω に等しい配位)。
定義 7.1 (Dobrushin の相互作用行列 ) ↗
x ≠ y x\ne y x = y について
C x y = 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\}, C x y = sup { ∥ γ x ( ⋅ ∣ ω ) − γ x ( ⋅ ∣ η ) ∥ : ω , η ∈ Ω , ω z = η z ( z = y ) } , C x x = 0 C_{xx}=0 C xx = 0 とおく。α ( γ ) = sup x ∑ y C x y \alpha(\gamma)=\sup_x\sum_yC_{xy} α ( γ ) = sup x ∑ y C x y を Dobrushin の定数 という。
C x y C_{xy} C x y は「y y y のスピン 1 個を変えたとき、x x x の条件付き分布がどれだけ動くか」の最大値です。関数 f f f の側では、y y y のスピン 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)\} δ y ( f ) = sup { ∣ f ( ω ) − f ( η ) ∣ : ω z = η z ( z = y )} を考え、
D = { f ∈ C ( Ω ) : ∑ y δ y ( f ) < ∞ } \mathcal{D}=\Bigl\{f\in C(\Omega):\ \sum_y\delta_y(f)<\infty\Bigr\} D = { f ∈ C ( Ω ) : y ∑ δ y ( f ) < ∞ }
とします。局所関数は D \mathcal{D} D に属します(有限個の y y y でしか δ y ( f ) ≠ 0 \delta_y(f)\ne0 δ y ( f ) = 0 にならない)。
補題 7.2 (変動の評価 ) ✓
(i) f ∈ C ( Ω ) f\in C(\Omega) f ∈ C ( Ω ) ならば sup f − inf f ≤ ∑ y δ y ( f ) \sup f-\inf f\le\sum_y\delta_y(f) sup f − inf f ≤ ∑ y δ y ( f ) 。
(ii) γ \gamma γ が準局所的な仕様で f ∈ D f\in\mathcal{D} f ∈ D 、x ∈ Z d x\in\mathbb{Z}^d x ∈ Z d ならば、δ x ( γ { x } f ) = 0 \delta_x(\gamma_{\{x\}}f)=0 δ x ( γ { x } f ) = 0 、y ≠ x y\ne x y = x について
δ y ( γ { x } f ) ≤ δ y ( f ) + C x y δ x ( f ) . \delta_y(\gamma_{\{x\}}f)\le\delta_y(f)+C_{xy}\,\delta_x(f). δ y ( γ { x } f ) ≤ δ y ( f ) + C x y δ x ( f ) . 特に ∑ y C x y < ∞ \sum_yC_{xy}<\infty ∑ y C x y < ∞ なら γ { x } f ∈ D \gamma_{\{x\}}f\in\mathcal{D} γ { x } f ∈ D 。
証明(補題 7.2) (i) ω , η ∈ Ω \omega,\eta\in\Omega ω , η ∈ Ω と Λ n ↑ Z d \Lambda_n\uparrow\mathbb{Z}^d Λ n ↑ Z d について、ω \omega ω から Λ n \Lambda_n Λ n の点を 1 個ずつ η \eta η の値に変えていくと ∣ f ( ω ) − f ( η Λ n ω Λ n c ) ∣ ≤ ∑ y ∈ Λ n δ y ( f ) |f(\omega)-f(\eta_{\Lambda_n}\omega_{\Lambda_n^c})|\le\sum_{y\in\Lambda_n}\delta_y(f) ∣ f ( ω ) − f ( η Λ n ω Λ n c ) ∣ ≤ ∑ y ∈ Λ n δ y ( f ) 。η Λ n ω Λ n c → η \eta_{\Lambda_n}\omega_{\Lambda_n^c}\to\eta η Λ n ω Λ n c → η なので、f f f の連続性から ∣ f ( ω ) − f ( η ) ∣ ≤ ∑ y δ y ( f ) |f(\omega)-f(\eta)|\le\sum_y\delta_y(f) ∣ f ( ω ) − f ( η ) ∣ ≤ ∑ y δ y ( f ) 。
(ii) γ { x } f ( ω ) \gamma_{\{x\}}f(\omega) γ { x } f ( ω ) は ω x c \omega_{x^c} ω x c だけに依存するので δ x ( γ { x } f ) = 0 \delta_x(\gamma_{\{x\}}f)=0 δ x ( γ { x } f ) = 0 。y ≠ x y\ne x y = x 、ω z = η z \omega_z=\eta_z ω z = η z (z ≠ y z\ne y z = y )とし、p = γ x ( ⋅ ∣ ω ) p=\gamma_x(\cdot\mid\omega) p = γ x ( ⋅ ∣ ω ) 、q = γ x ( ⋅ ∣ η ) q=\gamma_x(\cdot\mid\eta) q = γ x ( ⋅ ∣ η ) とおくと
γ { x } f ( ω ) − γ { x } f ( η ) = ∑ s p ( s ) ( f ( s ω x c ) − f ( s η x c ) ) + ∑ s ( p ( s ) − q ( s ) ) f ( s η x c ) . \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}). γ { x } f ( ω ) − γ { x } f ( η ) = s ∑ p ( s ) ( f ( s ω x c ) − f ( s η x c ) ) + s ∑ ( p ( s ) − q ( s ) ) f ( s η x c ) . 第 1 項の絶対値は δ y ( f ) \delta_y(f) δ y ( f ) 以下です(s ω x c s\omega_{x^c} s ω x c と s η x c s\eta_{x^c} s η x c は y y y でしか違わない)。第 2 項では、g ( s ) = f ( s η x c ) g(s)=f(s\eta_{x^c}) g ( s ) = f ( s η x c ) の最大値と最小値の中点を c c c とすれば ∑ ( p − q ) = 0 \sum(p-q)=0 ∑ ( p − q ) = 0 から
∣ ∑ s ( p ( s ) − q ( s ) ) ( g ( s ) − c ) ∣ ≤ ∑ s ∣ p ( s ) − q ( s ) ∣ ⋅ max g − min g 2 ≤ ∥ p − q ∥ δ x ( f ) ≤ C x y δ 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). s ∑ ( p ( s ) − q ( s )) ( g ( s ) − c ) ≤ s ∑ ∣ p ( s ) − q ( s ) ∣ ⋅ 2 max g − min g ≤ ∥ p − q ∥ δ x ( f ) ≤ C x y δ x ( f ) . 最後に、準局所性から γ { x } f \gamma_{\{x\}}f γ { x } f は連続で、変動の和は ∑ y δ y ( f ) + δ x ( f ) ∑ y C x y < ∞ \sum_y\delta_y(f)+\delta_x(f)\sum_yC_{xy}<\infty ∑ y δ y ( f ) + δ x ( f ) ∑ y C x y < ∞ です。
∎
定理 7.3 (Dobrushin の一意性定理 ) ✓
γ \gamma γ を準局所的な仕様とし、α ( γ ) < 1 \alpha(\gamma)<1 α ( γ ) < 1 とする。このとき G ( γ ) \mathcal{G}(\gamma) G ( γ ) はちょうど 1 点からなる。
証明の方針は次のとおりです。2 つの Gibbs 測度 μ , ν \mu,\nu μ , ν の期待値の差を、各座標の変動 δ y ( f ) \delta_y(f) δ y ( f ) の重みつき和で上から押さえます(この重みが「許容的」な b b b です)。1 点を再標本化するたびに、その点の重みが、それまでの重みの最大値の α \alpha α 倍以下に小さくなるので、これを繰り返せば重みを限りなく 0 0 0 に近づけられ、差が 0 0 0 になります。
証明(定理 7.3) 存在は 定理 5.1 なので、μ , ν ∈ G ( γ ) \mu,\nu\in\mathcal{G}(\gamma) μ , ν ∈ G ( γ ) について μ = ν \mu=\nu μ = ν を示します。α = α ( γ ) \alpha=\alpha(\gamma) α = α ( γ ) とおきます。有界な関数 b : Z d → [ 0 , ∞ ) b:\mathbb{Z}^d\to[0,\infty) b : Z d → [ 0 , ∞ ) が
∣ μ ( f ) − ν ( f ) ∣ ≤ ∑ y b y δ y ( f ) ( すべての f ∈ D ) |\mu(f)-\nu(f)|\le\sum_yb_y\,\delta_y(f)\qquad(\text{すべての }f\in\mathcal{D}) ∣ μ ( f ) − ν ( f ) ∣ ≤ y ∑ b y δ y ( f ) ( すべての f ∈ D ) を満たすとき、b b b を許容的 と呼びます。補題 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) ∣ μ ( f ) − ν ( f ) ∣ ≤ sup f − inf f ≤ ∑ y δ y ( f ) なので、定数 1 1 1 は許容的です。
1 点の改良. b b b が許容的で x ∈ Z d x\in\mathbb{Z}^d x ∈ Z d とし、b ′ b' b ′ を b x ′ = ∑ y C x y b y b'_x=\sum_yC_{xy}b_y b x ′ = ∑ y C x y b y 、b z ′ = b z b'_z=b_z b z ′ = b z (z ≠ x z\ne x z = x )で定めます。f ∈ D f\in\mathcal{D} f ∈ D について、DLR 方程式 μ γ { x } = μ \mu\gamma_{\{x\}}=\mu μ γ { x } = μ 、ν γ { x } = ν \nu\gamma_{\{x\}}=\nu ν γ { x } = ν と 補題 7.2 (ii) から
∣ μ ( f ) − ν ( f ) ∣ = ∣ μ ( γ { x } f ) − ν ( γ { x } f ) ∣ ≤ ∑ y ≠ x b y ( δ y ( f ) + C x y δ x ( f ) ) = ∑ y b y ′ δ 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). ∣ μ ( f ) − ν ( f ) ∣ = ∣ μ ( γ { x } f ) − ν ( γ { x } f ) ∣ ≤ y = x ∑ b y ( δ y ( f ) + C x y δ x ( f ) ) = y ∑ b y ′ δ y ( f ) . よって b ′ b' b ′ も許容的です。
有限集合での改良. β 0 = sup y b y \beta_0=\sup_yb_y β 0 = sup y b y とし、有限集合 Λ = { x 1 , … , x N } \Lambda=\{x_1,\dots,x_N\} Λ = { x 1 , … , x N } の点で順に 1 点の改良を行います。改良された値は ∑ y C x i y b y ≤ α β 0 ≤ β 0 \sum_yC_{x_iy}b_y\le\alpha\beta_0\le\beta_0 ∑ y C x i y b y ≤ α β 0 ≤ β 0 なので、各段階で上限は β 0 \beta_0 β 0 以下に保たれ、改良済みの座標は α β 0 \alpha\beta_0 α β 0 以下、改良していない座標は β 0 \beta_0 β 0 以下です。したがって、f ∈ D f\in\mathcal{D} f ∈ 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). ∣ μ ( f ) − ν ( f ) ∣ ≤ α β 0 y ∈ Λ ∑ δ y ( f ) + β 0 y ∈ / Λ ∑ δ y ( f ) . Λ ↑ Z d \Lambda\uparrow\mathbb{Z}^d Λ ↑ Z d とすれば第 2 項は 0 0 0 に収束するので、定数 α β 0 \alpha\beta_0 α β 0 も許容的です。
結論. 定数 1 1 1 から始めて繰り返せば、すべての k k k で定数 α k \alpha^k α k が許容的で、∣ μ ( f ) − ν ( f ) ∣ ≤ α k ∑ y δ y ( f ) → 0 |\mu(f)-\nu(f)|\le\alpha^k\sum_y\delta_y(f)\to0 ∣ μ ( f ) − ν ( f ) ∣ ≤ α k ∑ y δ y ( f ) → 0 。局所関数は D \mathcal{D} 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 osc Φ A = sup Φ A − inf Φ A とする。
sup x ∑ A ∋ x ( ∣ A ∣ − 1 ) osc Φ A < 1 2 \sup_x\sum_{A\ni x}(|A|-1)\operatorname{osc}\Phi_A<\frac12 x sup A ∋ x ∑ ( ∣ A ∣ − 1 ) osc Φ A < 2 1 ならば ∣ G ( Φ ) ∣ = 1 |\mathcal{G}(\Phi)|=1 ∣ G ( Φ ) ∣ = 1 。
証明(系 7.4) x ≠ y x\ne y x = y 、ω z = η z \omega_z=\eta_z ω z = η z (z ≠ y z\ne y z = y )とします。γ x ( s ∣ ω ) ∝ e − H { x } ( s ω x c ) \gamma_x(s\mid\omega)\propto e^{-H_{\{x\}}(s\omega_{x^c})} γ x ( s ∣ ω ) ∝ e − H { x } ( s ω x c ) で、u ( s ) = H { x } ( s η x c ) − H { x } ( s ω x c ) u(s)=H_{\{x\}}(s\eta_{x^c})-H_{\{x\}}(s\omega_{x^c}) u ( s ) = H { x } ( s η x c ) − H { x } ( s ω x c ) には x , y x,y x , y の両方を含む A A A だけが寄与するので、∣ u ( s ) ∣ ≤ r x y : = ∑ A ∋ x , y osc Φ A |u(s)|\le r_{xy}:=\sum_{A\ni x,y}\operatorname{osc}\Phi_A ∣ u ( s ) ∣ ≤ r x y := ∑ A ∋ x , y osc Φ A です。p = γ x ( ⋅ ∣ ω ) p=\gamma_x(\cdot\mid\omega) p = γ x ( ⋅ ∣ ω ) 、q = γ x ( ⋅ ∣ η ) q=\gamma_x(\cdot\mid\eta) q = γ x ( ⋅ ∣ η ) とおくと q ( s ) = p ( s ) e − u ( s ) / ∑ t p ( t ) e − u ( t ) ≥ e − 2 r x y p ( s ) q(s)=p(s)e^{-u(s)}/\sum_tp(t)e^{-u(t)}\ge e^{-2r_{xy}}p(s) q ( s ) = p ( s ) e − u ( s ) / ∑ t p ( t ) e − u ( t ) ≥ e − 2 r x y p ( s ) なので
∥ p − q ∥ = ∑ s ( p ( s ) − q ( s ) ) + ≤ ∑ s p ( s ) ( 1 − e − 2 r x y ) ≤ 2 r x y . \|p-q\|=\sum_s(p(s)-q(s))^+\le\sum_sp(s)(1-e^{-2r_{xy}})\le2r_{xy}. ∥ p − q ∥ = s ∑ ( p ( s ) − q ( s ) ) + ≤ s ∑ p ( s ) ( 1 − e − 2 r x y ) ≤ 2 r x y . よって ∑ y C x y ≤ 2 ∑ y ≠ x ∑ A ∋ x , y osc Φ 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 ∑ y C x y ≤ 2 ∑ y = x ∑ A ∋ x , y osc Φ A = 2 ∑ A ∋ x ( ∣ A ∣ − 1 ) osc Φ A < 1 で、定理 7.3 を使います。
∎
命題 7.5 (Ising 模型の高温での一意性 ) ✓
S = { − 1 , 1 } S=\{-1,1\} S = { − 1 , 1 } とし、最近接 Ising 模型の相互作用 Φ { x , y } ( ω ) = − β J ω x ω y \Phi_{\{x,y\}}(\omega)=-\beta J\omega_x\omega_y Φ { x , y } ( ω ) = − β J ω x ω y (∣ x − y ∣ = 1 |x-y|=1 ∣ x − y ∣ = 1 )、Φ { x } ( ω ) = − β h ω x \Phi_{\{x\}}(\omega)=-\beta h\omega_x Φ { x } ( ω ) = − β h ω x (β > 0 \beta\gt0 β > 0 、J , h ∈ R J,h\in\mathbb{R} J , h ∈ R )、それ以外の Φ A = 0 \Phi_A=0 Φ A = 0 を考える。2 d tanh ( β ∣ J ∣ ) < 1 2d\tanh(\beta|J|)\lt1 2 d tanh ( β ∣ J ∣ ) < 1 ならば、h h h によらず ∣ G ( Φ ) ∣ = 1 |\mathcal{G}(\Phi)|=1 ∣ G ( Φ ) ∣ = 1 である。
証明(命題 7.5) x x x を固定し、局所場を m ( ω ) = J ∑ y ∼ x ω y + h m(\omega)=J\sum_{y\sim x}\omega_y+h m ( ω ) = J ∑ y ∼ x ω y + h とします。H { x } ( s ω x c ) = − β s m ( ω ) + c ( ω ) H_{\{x\}}(s\omega_{x^c})=-\beta sm(\omega)+c(\omega) H { x } ( s ω x c ) = − β s m ( ω ) + c ( ω ) (c ( ω ) c(\omega) c ( ω ) は s s s によらない)なので、定義 3.3 から
γ x ( + 1 ∣ ω ) = e β m e β m + e − β m = 1 + tanh β m 2 , γ x ( − 1 ∣ ω ) = 1 − tanh β m 2 \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 ( + 1 ∣ ω ) = e β m + e − β m e β m = 2 1 + tanh β m , γ x ( − 1 ∣ ω ) = 2 1 − tanh β m です。γ x ( ⋅ ∣ ω ) \gamma_x(\cdot\mid\omega) γ x ( ⋅ ∣ ω ) は ω \omega ω の x x x の隣接点での値だけに依存するので、y ≁ x y\not\sim x y ∼ x なら C x y = 0 C_{xy}=0 C x y = 0 です。y ∼ x y\sim x y ∼ x で ω , η \omega,\eta ω , η が y y y でだけ異なるとき、局所場の差は ± 2 J \pm2J ± 2 J なので、2 点上の確率測度の全変動距離 ∥ p − q ∥ = ∑ s ( p ( s ) − q ( s ) ) + = ∣ p ( + 1 ) − q ( + 1 ) ∣ \|p-q\|=\sum_s(p(s)-q(s))^+=|p(+1)-q(+1)| ∥ p − q ∥ = ∑ s ( p ( s ) − q ( s ) ) + = ∣ p ( + 1 ) − q ( + 1 ) ∣ から
C x y ≤ sup m ∈ R 1 2 ∣ 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| C x y ≤ m ∈ R sup 2 1 tanh ( β ( m + 2∣ J ∣ )) − tanh ( β m ) です(J < 0 J\lt0 J < 0 のときは m m m を m − 2 ∣ J ∣ m-2|J| m − 2∣ J ∣ に置き換える)。a = 2 β ∣ J ∣ ≥ 0 a=2\beta|J|\ge0 a = 2 β ∣ J ∣ ≥ 0 、g ( t ) = tanh ( t + a ) − tanh t ≥ 0 g(t)=\tanh(t+a)-\tanh t\ge0 g ( t ) = tanh ( t + a ) − tanh t ≥ 0 とおくと g ′ ( t ) = cosh − 2 ( t + a ) − cosh − 2 t g'(t)=\cosh^{-2}(t+a)-\cosh^{-2}t g ′ ( t ) = cosh − 2 ( t + a ) − cosh − 2 t は ∣ t + a ∣ < ∣ t ∣ |t+a|\lt|t| ∣ t + a ∣ < ∣ t ∣ 、すなわち t < − a / 2 t\lt-a/2 t < − a /2 で正、t > − a / 2 t\gt-a/2 t > − a /2 で負なので、g g g の最大値は g ( − a / 2 ) = 2 tanh ( a / 2 ) g(-a/2)=2\tanh(a/2) g ( − a /2 ) = 2 tanh ( a /2 ) です。よって C x y ≤ tanh ( β ∣ J ∣ ) C_{xy}\le\tanh(\beta|J|) C x y ≤ tanh ( β ∣ J ∣ ) で、隣接点は 2 d 2d 2 d 個なので α ( γ Φ ) ≤ 2 d tanh ( β ∣ J ∣ ) < 1 \alpha(\gamma^\Phi)\le2d\tanh(\beta|J|)\lt1 α ( γ Φ ) ≤ 2 d tanh ( β ∣ J ∣ ) < 1 です。γ Φ \gamma^\Phi γ Φ は準局所的な仕様(命題 3.4 。sup x ∑ A ∋ x ∥ Φ A ∥ ∞ = β ( 2 d ∣ J ∣ + ∣ h ∣ ) < ∞ \sup_x\sum_{A\ni x}\|\Phi_A\|_\infty=\beta(2d|J|+|h|)\lt\infty sup x ∑ A ∋ x ∥ Φ A ∥ ∞ = β ( 2 d ∣ J ∣ + ∣ h ∣ ) < ∞ )なので、定理 7.3 から ∣ G ( Φ ) ∣ = 1 |\mathcal{G}(\Phi)|=1 ∣ G ( Φ ) ∣ = 1 です。
∎
例 7.6 (Ising 模型の高温での一意性 )
最近接 Ising 模型(逆温度 β \beta β 、結合 J J J 、外場 h h h )では、x x x の局所場を m = J ∑ y ∼ x ω y + h m=J\sum_{y\sim x}\omega_y+h m = J ∑ y ∼ x ω y + h とすると γ x ( + 1 ∣ ω ) = ( 1 + tanh β m ) / 2 \gamma_x(+1\mid\omega)=(1+\tanh\beta m)/2 γ x ( + 1 ∣ ω ) = ( 1 + tanh β m ) /2 です。隣の y y y のスピンを反転すると m m m は m ± 2 J m\pm2J m ± 2 J に変わるので
C x y ≤ sup m ∈ R 1 2 ∣ tanh ( β ( m + 2 J ) ) − 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|) C x y ≤ m ∈ R sup 2 1 tanh ( β ( m + 2 J )) − tanh ( β m ) = tanh ( β ∣ J ∣ ) (上限は m = − J m=-J m = − J でとられる)。隣接点は 2 d 2d 2 d 個なので α ≤ 2 d tanh ( β ∣ J ∣ ) \alpha\le2d\tanh(\beta|J|) α ≤ 2 d tanh ( β ∣ J ∣ ) で、
2 d tanh ( β ∣ J ∣ ) < 1 2d\tanh(\beta|J|)<1 2 d tanh ( β ∣ J ∣ ) < 1 なら外場 h h h によらず Gibbs 測度はただ一つです(命題 7.5 )。d = 2 d=2 d = 2 では β ∣ J ∣ < artanh ( 1 / 4 ) ≈ 0.255 \beta|J|<\operatorname{artanh}(1/4)\approx0.255 β ∣ J ∣ < artanh ( 1/4 ) ≈ 0.255 です(正確な臨界値は β c J = 1 2 log ( 1 + 2 ) ≈ 0.441 \beta_cJ=\frac12\log(1+\sqrt2)\approx0.441 β c J = 2 1 log ( 1 + 2 ) ≈ 0.441 )。
演習 8.1 易
Φ = 0 \Phi=0 Φ = 0 の Gibbs 仕様(例 3.5 )について、G \mathcal{G} G が S S S 上の一様分布の無限直積だけからなることを、(1) DLR 方程式から直接、(2) 定理 7.3 から、の 2 通りで示してください。
解答 (1) μ ∈ G \mu\in\mathcal{G} μ ∈ G 、Λ ⋐ Z d \Lambda\Subset\mathbb{Z}^d Λ ⋐ 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|} μ ( ω Λ = ξ ) = μ γ Λ ( ω Λ = ξ ) = ∫ ∣ S ∣ − ∣Λ∣ d μ = ∣ S ∣ − ∣Λ∣ 。有限次元分布が一様分布の直積と一致するので、柱集合(π-系)の上で一致し、測度の一意性(定理 4.5[測度空間と Carathéodory の拡張定理] )から μ \mu μ は一様分布の無限直積です。逆にこの直積が DLR 方程式を満たすことは、Λ \Lambda Λ 上の一様分布と外の分布が独立であることから従います。
(2) γ x ( ⋅ ∣ ω ) \gamma_x(\cdot\mid\omega) γ x ( ⋅ ∣ ω ) は ω \omega ω によらない一様分布なので C x y = 0 C_{xy}=0 C x y = 0 、α = 0 < 1 \alpha=0<1 α = 0 < 1 です。
演習 8.2 標準
μ 1 ≠ μ 2 \mu_1\ne\mu_2 μ 1 = μ 2 を極値 Gibbs 測度とし、μ = 1 2 ( μ 1 + μ 2 ) \mu=\frac12(\mu_1+\mu_2) μ = 2 1 ( μ 1 + μ 2 ) とします。系 6.3 の尾部事象 B B B について μ ( B ) = 1 / 2 \mu(B)=1/2 μ ( B ) = 1/2 を示し、μ \mu μ が局所関数の相関について「遠くで無相関」にならない例を作ってください。具体的には S = { − 1 , 1 } S=\{-1,1\} S = { − 1 , 1 } で μ 1 ( ω x ) = m = − μ 2 ( ω x ) \mu_1(\omega_x)=m=-\mu_2(\omega_x) μ 1 ( ω x ) = m = − μ 2 ( ω x ) (m > 0 m>0 m > 0 、すべての x x x )のとき、μ ( ω 0 ω x ) − μ ( ω 0 ) μ ( ω x ) ≥ m 2 \mu(\omega_0\omega_x)-\mu(\omega_0)\mu(\omega_x)\ge m^2 μ ( ω 0 ω x ) − μ ( ω 0 ) μ ( ω x ) ≥ m 2 を示してください(極値測度では μ i ( ω 0 ω x ) ≥ μ i ( ω 0 ) μ i ( ω x ) \mu_i(\omega_0\omega_x)\ge\mu_i(\omega_0)\mu_i(\omega_x) μ i ( ω 0 ω x ) ≥ μ i ( ω 0 ) μ i ( ω x ) が成り立つと仮定してよい)。
解答 μ 1 ( B ) = 1 \mu_1(B)=1 μ 1 ( B ) = 1 、μ 2 ( B ) = 0 \mu_2(B)=0 μ 2 ( B ) = 0 なので μ ( B ) = 1 / 2 \mu(B)=1/2 μ ( B ) = 1/2 。相関については μ ( ω 0 ) = μ ( ω x ) = 0 \mu(\omega_0)=\mu(\omega_x)=0 μ ( ω 0 ) = μ ( ω x ) = 0 で、
μ ( ω 0 ω x ) = 1 2 ( μ 1 ( ω 0 ω x ) + μ 2 ( ω 0 ω x ) ) ≥ 1 2 ( m 2 + m 2 ) = m 2 . \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. μ ( ω 0 ω x ) = 2 1 ( μ 1 ( ω 0 ω x ) + μ 2 ( ω 0 ω x ) ) ≥ 2 1 ( m 2 + m 2 ) = m 2 . したがって共分散は ∣ x ∣ |x| ∣ x ∣ によらず m 2 m^2 m 2 以上で、∣ x ∣ → ∞ |x|\to\infty ∣ x ∣ → ∞ でも 0 0 0 に近づきません。混合状態では、遠く離れたスピンが「どちらの相にいるか」という尾部の情報を通じて相関します。仮定した不等式は強磁性 Ising 模型の ± \pm ± 状態では GKS 不等式(定理 3.2[相関不等式] )の極限として成り立ちます。
演習 8.3 標準
S = { − 1 , 1 } S=\{-1,1\} S = { − 1 , 1 } 、Φ { x , y } = − J x y ω x ω y \Phi_{\{x,y\}}=-J_{xy}\omega_x\omega_y Φ { x , y } = − J x y ω x ω y (J x y = J y x J_{xy}=J_{yx} J x y = J y x 、sup x ∑ y ∣ J x y ∣ < ∞ \sup_x\sum_y|J_{xy}|<\infty sup x ∑ y ∣ J x y ∣ < ∞ )と外場 Φ { x } = − h x ω x \Phi_{\{x\}}=-h_x\omega_x Φ { x } = − h x ω x からなる Ising 模型(逆温度は J , h J,h J , h に含める)について、C x y ≤ tanh ∣ J x y ∣ C_{xy}\le\tanh|J_{xy}| C x y ≤ tanh ∣ J x y ∣ を示し、sup x ∑ y tanh ∣ J x y ∣ < 1 \sup_x\sum_y\tanh|J_{xy}|<1 sup x ∑ y tanh ∣ J x y ∣ < 1 なら Gibbs 測度が一意であることを導いてください。系 7.4 の条件と比べてどちらが強いですか。
解答 例 7.6 と同じく、x x x の局所場 m = ∑ y J x y ω y + h x m=\sum_yJ_{xy}\omega_y+h_x m = ∑ y J x y ω y + h x について γ x ( + 1 ∣ ω ) = ( 1 + tanh m ) / 2 \gamma_x(+1\mid\omega)=(1+\tanh m)/2 γ x ( + 1 ∣ ω ) = ( 1 + tanh m ) /2 で、y y y のスピンを反転すると m m m は ± 2 J x y \pm2J_{xy} ± 2 J x y だけ変わります。tanh \tanh tanh の差は sup m 1 2 ∣ tanh ( m + 2 a ) − tanh m ∣ = tanh ∣ a ∣ \sup_m\frac12|\tanh(m+2a)-\tanh m|=\tanh|a| sup m 2 1 ∣ tanh ( m + 2 a ) − tanh m ∣ = tanh ∣ a ∣ なので C x y ≤ tanh ∣ J x y ∣ C_{xy}\le\tanh|J_{xy}| C x y ≤ tanh ∣ J x y ∣ 。あとは 定理 7.3 です。系 7.4 の条件は osc Φ { x , y } = 2 ∣ J x y ∣ \operatorname{osc}\Phi_{\{x,y\}}=2|J_{xy}| osc Φ { x , y } = 2∣ J x y ∣ から sup x ∑ y 2 ∣ J x y ∣ < 1 / 2 \sup_x\sum_y2|J_{xy}|<1/2 sup x ∑ y 2∣ J x y ∣ < 1/2 、すなわち ∑ y ∣ J x y ∣ < 1 / 4 \sum_y|J_{xy}|<1/4 ∑ y ∣ J x y ∣ < 1/4 で、tanh a ≤ a \tanh a\le a tanh a ≤ a なのでこれは今の条件 ∑ y tanh ∣ J x y ∣ < 1 \sum_y\tanh|J_{xy}|<1 ∑ y tanh ∣ J x y ∣ < 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 方程式を扱う章。
執筆: Claude Code(Opus 5.5)、Claude Code(Sonnet 5.5)・監査 3 回(Claude Code)・最終 2026-10-12