# ランダムクラスター表示と浸透：FK 表示・Edwards–Sokal 結合・FKG による比較と長距離模型

> Potts・Ising 模型を辺の浸透模型に書き換える FK（ランダムクラスター）表示と Edwards–Sokal 結合を証明し、FKG 不等式で独立浸透と比較して高温・低温の判定を導く。長距離浸透と 1/r² Ising の不連続性（Aizenman–Newman・ACCN）の骨格も述べる。
> https://rikai.mugen-giken.com/mathematics/gibbs-measures/random-cluster-percolation

## 0. この記事の要点

磁石の模型で「隣どうしが揃いやすい」という性質を、格子の辺を「つなぐ」「つながない」とランダムに決める問題に言い換えると、磁石の性質が「大きなつながりができるか」という問いに変わります。この記事は、その言い換えを正確に述べ、道や輪郭の数え上げだけで高温・低温での磁化の有無を判定し、さらに遠くまで届く結合のある 1 次元の模型にまで進みます。

- **FK 表示**（<Ref to="thm-edwards-sokal" />）: $q$ 状態 Potts 模型（各点が $q$ 色のどれかをとり、隣り合う点が同じ色だと有利になる模型。$q=2$ が Ising 模型）は、辺を確率的に開閉する**ランダムクラスター模型** $\varphi_{p,q}(\omega)\propto\prod_ep_e^{\omega_e}(1-p_e)^{1-\omega_e}\,q^{k(\omega)}$（$k$ はクラスターの個数）と、**Edwards–Sokal 結合**で結ばれます。スピンの相関は連結確率になり、Ising 模型では $\langle\sigma_x\sigma_y\rangle=\varphi(x\leftrightarrow y)$、$+$ 境界条件の磁化は「境界とつながる確率」です。
- **FKG**（<Ref to="thm-rc-fkg" />）: $q\ge1$ ならランダムクラスター測度は [FKG の束条件](/physics/constructive-qft/correlation-inequalities#def-fkg-lattice-condition)を満たします。鍵はクラスターの個数の**超モジュラ性**（辺が多い配置ほど、辺を 1 本足したときのクラスター数の減り方が小さい、という性質）$k(\omega\vee\eta)+k(\omega\wedge\eta)\ge k(\omega)+k(\eta)$ です。
- **比較不等式**（<Ref to="thm-comparison" />）: $q\ge1$ では $\varphi_{p,q}$ は独立浸透 $P_p$ で上から、$P_{p'}$（$p'=p/(p+q(1-p))$）で下から確率的に抑えられます。これだけで、Ising 模型の高温での磁化の消失（経路の数え上げ）と低温での[自発磁化](/mathematics/gibbs-measures/ising-peierls#def-spontaneous-magnetization)（[Peierls の章](/mathematics/gibbs-measures/ising-peierls) の輪郭の数え上げ）が得られます。
- **長距離模型**: 1 次元で $\sum_rr\,p_r<\infty$ なら浸透は起きません（<Ref to="prop-lr-no-percolation" />）。境界の $p_r\sim\beta/r^2$ では、Aizenman–Newman と Aizenman–Chayes–Chayes–Newman（ACCN）が「浸透確率 $\theta$ は $0$ か $\theta^2\beta\ge1$」という二分律を示し、浸透転移と $1/r^2$ Ising 模型の磁化が臨界点で**不連続**であることを導きました（主張と骨格）。

<Aside type="tip" title="この記事で使う事前知識">

- **最近接 Ising 模型**: 格子の各点にスピン $\pm1$ が乗り、隣り合う点が同じ向きだとエネルギーが下がる模型。→ [最近接 Ising 模型](/mathematics/gibbs-measures/ising-peierls#def-ising)
- **自発磁化**: 外場を $0$ にした無限体積の極限で、$+$ 境界条件のもとでも残る磁化。→ [自発磁化](/mathematics/gibbs-measures/ising-peierls#def-spontaneous-magnetization)
- **Peierls の議論**: 低温で、スピンの向きが変わる境目（輪郭）の数え上げから自発磁化を示す方法。→ [Ising 模型と Peierls の議論](/mathematics/gibbs-measures/ising-peierls)
- **FKG の束条件**: 測度の重みが $\varphi(\omega\vee\eta)\varphi(\omega\wedge\eta)\ge\varphi(\omega)\varphi(\eta)$ を満たすという条件。→ [FKG の束条件](/physics/constructive-qft/correlation-inequalities#def-fkg-lattice-condition)
- **FKG 不等式**: 束条件を満たす測度では、増加する 2 つの量が正の相関をもつという不等式。→ [相関不等式](/physics/constructive-qft/correlation-inequalities)

</Aside>

## 1. 動機

Ising 模型の相関不等式（[Griffiths・FKG・Lebowitz](/physics/constructive-qft/correlation-inequalities)）や、[ランダムカレント表示](/physics/constructive-qft/random-current) は、スピンの相関を組合せ的な対象（カレントの経路）の言葉に書き換えて評価する方法でした。Fortuin と Kasteleyn（1969–72）は、別の書き換えを見つけました。各辺の Boltzmann 因子を「辺を開く」「閉じる」の和に分解すると、Potts 模型は**辺の浸透模型**に化け、スピンの長距離秩序は「無限に大きなクラスターがあること」に対応します。

この表示の利点は 3 つあります。(1) スピンの値の対称性によらない議論になり、Ising だけでなく $q$ 状態 Potts 模型を一様に扱えます。(2) 独立浸透との比較という強力な道具が使えます。(3) Swendsen–Wang のクラスター更新（クラスターごとにまとめてスピンを塗り替えるモンテカルロ法。[格子モンテカルロ法](/physics/lattice-field-theory/lattice-monte-carlo)）は Edwards–Sokal 結合そのものです。本章では、この表示を証明し、比較不等式から相転移の上下からの評価を導き、最後に長距離模型での不連続転移の話題を紹介します。

## 2. ランダムクラスター模型

$G=(V,E)$ を有限グラフ、$W\subset V$ を（空でもよい）**境界**の頂点の集合とします。辺の配置 $\omega\in\{0,1\}^E$ で $\omega_e=1$ の辺を**開**、$0$ を**閉**と呼び、開いた辺の集合も $\omega$ と書きます。$x\leftrightarrow y$ は「$x$ と $y$ が開いた辺の道でつながる」を、$x\leftrightarrow W$ は「$x$ が $W$ のある点とつながる」を表します。$k(\omega)$ はグラフ $(V,\omega)$ の連結成分（クラスター、孤立点を含む）の個数、$k_W(\omega)$ はそのうち $W$ と交わらないものの個数です。

<Definition id="def-random-cluster" title="ランダムクラスター測度">
$p=(p_e)_{e\in E}$、$p_e\in(0,1)$、$q>0$ について
$$
\varphi^W_{p,q}(\omega)=\frac1{Z^W_{p,q}}\prod_{e\in E}p_e^{\omega_e}(1-p_e)^{1-\omega_e}\;q^{k_W(\omega)}
$$
を**ランダムクラスター測度**（FK 測度）という。$W=\emptyset$ のとき（$k_\emptyset=k$）を**自由境界条件**、$W\ne\emptyset$ のときを**配線（wired）境界条件**と呼ぶ。$q=1$ では独立な Bernoulli 浸透 $P_p$ である。
</Definition>

配線境界条件の名前は、$W$ と交わるクラスターをすべて 1 つに数える（$W$ の点を外で配線してつないだと見る）ことから来ています。$q>1$ ではクラスターが多い配置が優遇され、辺が閉じやすくなります。

**Potts 模型.** $q\in\{2,3,\dots\}$、$J_e\ge0$、$\beta>0$ とし、$\sigma\in\{1,\dots,q\}^V$ で $W$ 上は $\sigma\equiv1$ のものに
$$
\pi^W(\sigma)\propto\prod_{e=\{x,y\}\in E}\exp\bigl(\beta J_e(\mathbf{1}\{\sigma_x=\sigma_y\}-1)\bigr)
$$
を与えたものを、境界 $W$ で色 $1$ に固定した **Potts 模型**と呼びます。$q=2$ で色 $1,2$ を $+1,-1$ と読めば、$\sigma_x\sigma_y=2\cdot\mathbf{1}\{\sigma_x=\sigma_y\}-1$ なので、結合 $J_e/2$ の Ising 模型です。

## 3. Edwards–Sokal 結合と FK 表示

$p_e=1-e^{-\beta J_e}$ とおくと、各辺の因子は
$$
e^{\beta J_e(\mathbf{1}\{\sigma_x=\sigma_y\}-1)}=(1-p_e)+p_e\mathbf{1}\{\sigma_x=\sigma_y\}
$$
と分解できます（一致すれば $1$、しなければ $1-p_e=e^{-\beta J_e}$）。右辺の 2 項を「辺を閉じる」「辺を開く（ただし両端の色が等しいとき）」と読むのが FK 表示の出発点です。

<Theorem id="thm-edwards-sokal" title="Edwards–Sokal 結合">
$q\in\{2,3,\dots\}$、$p_e=1-e^{-\beta J_e}\in(0,1)$ とし、$W$ 上で $\sigma\equiv1$ となる $(\sigma,\omega)$ に
$$
\mathbb{P}(\sigma,\omega)=\frac1{\mathcal{Z}}\prod_{e=\{x,y\}\in E}\Bigl((1-p_e)\mathbf{1}\{\omega_e=0\}+p_e\mathbf{1}\{\omega_e=1\}\mathbf{1}\{\sigma_x=\sigma_y\}\Bigr)
$$
を与える。このとき次が成り立つ。

(i) $\sigma$ の周辺分布は Potts 模型 $\pi^W$ である。

(ii) $\omega$ の周辺分布はランダムクラスター測度 $\varphi^W_{p,q}$ である。

(iii) $\omega$ を与えたとき、$\sigma$ の条件付き分布は「$W$ と交わるクラスターは色 $1$、他のクラスターには独立に一様な色を塗り、クラスター内は同色」である。

(iv) $\sigma$ を与えたとき、$\omega_e$ は独立で、$\sigma_x=\sigma_y$ の辺では確率 $p_e$ で開、$\sigma_x\ne\sigma_y$ の辺は閉である。
</Theorem>

<Proof of="thm-edwards-sokal">
**(i)** $\sigma$ を固定して $\omega$ について和をとると、辺ごとの和 $\sum_{\omega_e\in\{0,1\}}$ に分かれ、各辺で $(1-p_e)+p_e\mathbf{1}\{\sigma_x=\sigma_y\}=e^{\beta J_e(\mathbf{1}\{\sigma_x=\sigma_y\}-1)}$。

**(ii)** $\omega$ を固定すると、$\mathbb{P}(\sigma,\omega)\ne0$ となるのは開いた辺の両端が同色の $\sigma$、すなわちクラスターごとに一定の $\sigma$ で、そのとき値は $\prod_ep_e^{\omega_e}(1-p_e)^{1-\omega_e}/\mathcal{Z}$ です。$W$ と交わるクラスターの色は $1$ に決まり、他の $k_W(\omega)$ 個のクラスターの色は自由なので、そのような $\sigma$ は $q^{k_W(\omega)}$ 個あります。

**(iii)** (ii) の計算から、$\omega$ を固定したとき許される $\sigma$ はすべて同じ重みをもち、それはちょうど述べた塗り方です。

**(iv)** $\sigma$ を固定すると $\mathbb{P}(\sigma,\omega)$ は辺ごとの因子の積で、$\sigma_x=\sigma_y$ の辺の因子は $\omega_e=0,1$ に対して $1-p_e,p_e$、$\sigma_x\ne\sigma_y$ の辺では $1-p_e,0$ です。
</Proof>

<Corollary id="cor-correlations-connectivity" title="相関は連結確率">
上の設定で、$x,y\in V$ について
$$
\pi^\emptyset(\sigma_x=\sigma_y)=\frac1q+\Bigl(1-\frac1q\Bigr)\varphi^\emptyset_{p,q}(x\leftrightarrow y),\qquad
\pi^W(\sigma_x=1)=\frac1q+\Bigl(1-\frac1q\Bigr)\varphi^W_{p,q}(x\leftrightarrow W).
$$
特に $q=2$（結合 $J_e/2$ の Ising 模型）では
$$
\langle\sigma_x\sigma_y\rangle^{\text{自由}}=\varphi^\emptyset_{p,2}(x\leftrightarrow y),\qquad\langle\sigma_x\rangle^{W,+}=\varphi^W_{p,2}(x\leftrightarrow W).
$$
</Corollary>

<Proof of="cor-correlations-connectivity">
(iii) から、$x\leftrightarrow y$ なら $\sigma_x=\sigma_y$、そうでなければ $\sigma_x,\sigma_y$ の少なくとも一方は独立に一様な色なので一致する確率は $1/q$ です（一方が $W$ につながり色 $1$ でも同じ）。$\mathbb{P}$ について平均すれば第 1 式です。第 2 式も同様で、$x\leftrightarrow W$ なら $\sigma_x=1$、そうでなければ確率 $1/q$。Ising では $\langle\sigma_x\sigma_y\rangle=2\pi(\sigma_x=\sigma_y)-1$、$\langle\sigma_x\rangle=2\pi(\sigma_x=1)-1$ です。
</Proof>

**Ising 模型の $+$ 境界条件.** $\mathbb{Z}^d$ の最近接 Ising 模型（逆温度 $\beta$、結合 $J>0$、$h=0$）の $\Lambda\Subset\mathbb{Z}^d$ での $+$ 境界条件の分布 $\mu^+_\Lambda$（[Peierls の章](/mathematics/gibbs-measures/ising-peierls) の <Ref to="mathematics/gibbs-measures/ising-peierls#def-ising" />）は、頂点 $V=\Lambda\cup\partial^{\mathrm{out}}\Lambda$（$\partial^{\mathrm{out}}\Lambda$ は $\Lambda$ の外で $\Lambda$ に隣接する点）、辺 $E_\Lambda$（少なくとも一方の端点が $\Lambda$ にある最近接対）、$W=\partial^{\mathrm{out}}\Lambda$ の上の $q=2$ の Potts 模型（結合 $2J$）です。よって $p=1-e^{-2\beta J}$ として
$$
\langle\sigma_0\rangle^+_\Lambda=\varphi^{W}_{p,2}(0\leftrightarrow\partial^{\mathrm{out}}\Lambda)
$$
です。自発磁化は「原点が遠くの境界とつながる確率」、すなわち配線ランダムクラスター模型の**浸透確率**です。

## 4. FKG 不等式と比較

$\{0,1\}^E$ に成分ごとの順序を入れ、$\omega\vee\eta$、$\omega\wedge\eta$ を成分ごとの最大・最小（開いた辺の集合としては和・共通部分）とします。関数 $f$ が**増加**とは $\omega\le\eta\Rightarrow f(\omega)\le f(\eta)$、事象が増加とはその指示関数が増加であることです。$\{x\leftrightarrow y\}$ は増加事象です。

次の補題は、辺を足したときにクラスターの個数がどれだけ減るかを比べるもので、FKG の束条件を確かめるのに使う唯一の入力です。直感的には、辺がすでに多い配置では新しい辺が別々のクラスターをつなぐ機会が少なく、クラスター数の減り方も小さくなります。下の図は 4 頂点の四角形での例で、実線が開いた辺です。

<Figure caption="四角形の 4 辺のうち、ω は上と右の 2 辺、η は下と左の 2 辺を開いたもの。クラスター数 k を数えると、k(ω∨η)+k(ω∧η)=1+4=5 ≥ 2+2=k(ω)+k(η) となる。">
<svg viewBox="0 0 400 115" width="100%" role="img" aria-label="四角形の辺の開閉による 4 つの配置ω、η、ω∨η、ω∧η と、それぞれのクラスター数 2、2、1、4">
  <g>
    <line x1="20" y1="15" x2="80" y2="15" stroke="var(--sl-color-accent)" stroke-width="4" />
    <line x1="80" y1="15" x2="80" y2="75" stroke="var(--sl-color-accent)" stroke-width="4" />
    <line x1="20" y1="75" x2="80" y2="75" stroke="currentColor" stroke-width="1" stroke-dasharray="3 4" opacity="0.4" />
    <line x1="20" y1="15" x2="20" y2="75" stroke="currentColor" stroke-width="1" stroke-dasharray="3 4" opacity="0.4" />
    <circle cx="20" cy="15" r="5" fill="currentColor" />
    <circle cx="80" cy="15" r="5" fill="currentColor" />
    <circle cx="80" cy="75" r="5" fill="currentColor" />
    <circle cx="20" cy="75" r="5" fill="currentColor" />
    <text x="50" y="104" text-anchor="middle" font-size="13" fill="currentColor">ω：k=2</text>
  </g>
  <g>
    <line x1="120" y1="15" x2="180" y2="15" stroke="currentColor" stroke-width="1" stroke-dasharray="3 4" opacity="0.4" />
    <line x1="180" y1="15" x2="180" y2="75" stroke="currentColor" stroke-width="1" stroke-dasharray="3 4" opacity="0.4" />
    <line x1="120" y1="75" x2="180" y2="75" stroke="var(--sl-color-accent)" stroke-width="4" />
    <line x1="120" y1="15" x2="120" y2="75" stroke="var(--sl-color-accent)" stroke-width="4" />
    <circle cx="120" cy="15" r="5" fill="currentColor" />
    <circle cx="180" cy="15" r="5" fill="currentColor" />
    <circle cx="180" cy="75" r="5" fill="currentColor" />
    <circle cx="120" cy="75" r="5" fill="currentColor" />
    <text x="150" y="104" text-anchor="middle" font-size="13" fill="currentColor">η：k=2</text>
  </g>
  <g>
    <line x1="220" y1="15" x2="280" y2="15" stroke="var(--sl-color-accent)" stroke-width="4" />
    <line x1="280" y1="15" x2="280" y2="75" stroke="var(--sl-color-accent)" stroke-width="4" />
    <line x1="220" y1="75" x2="280" y2="75" stroke="var(--sl-color-accent)" stroke-width="4" />
    <line x1="220" y1="15" x2="220" y2="75" stroke="var(--sl-color-accent)" stroke-width="4" />
    <circle cx="220" cy="15" r="5" fill="currentColor" />
    <circle cx="280" cy="15" r="5" fill="currentColor" />
    <circle cx="280" cy="75" r="5" fill="currentColor" />
    <circle cx="220" cy="75" r="5" fill="currentColor" />
    <text x="250" y="104" text-anchor="middle" font-size="13" fill="currentColor">ω∨η：k=1</text>
  </g>
  <g>
    <line x1="320" y1="15" x2="380" y2="15" stroke="currentColor" stroke-width="1" stroke-dasharray="3 4" opacity="0.4" />
    <line x1="380" y1="15" x2="380" y2="75" stroke="currentColor" stroke-width="1" stroke-dasharray="3 4" opacity="0.4" />
    <line x1="320" y1="75" x2="380" y2="75" stroke="currentColor" stroke-width="1" stroke-dasharray="3 4" opacity="0.4" />
    <line x1="320" y1="15" x2="320" y2="75" stroke="currentColor" stroke-width="1" stroke-dasharray="3 4" opacity="0.4" />
    <circle cx="320" cy="15" r="5" fill="currentColor" />
    <circle cx="380" cy="15" r="5" fill="currentColor" />
    <circle cx="380" cy="75" r="5" fill="currentColor" />
    <circle cx="320" cy="75" r="5" fill="currentColor" />
    <text x="350" y="104" text-anchor="middle" font-size="13" fill="currentColor">ω∧η：k=4</text>
  </g>
</svg>
</Figure>

<Lemma id="lem-supermodular" title="クラスターの個数の超モジュラ性">
すべての $\omega,\eta\in\{0,1\}^E$ と $W\subset V$ について
$$
k_W(\omega\vee\eta)+k_W(\omega\wedge\eta)\ge k_W(\omega)+k_W(\eta).
$$
また、$W$ と交わるクラスターの個数 $N_W(\omega)=k(\omega)-k_W(\omega)$ は $\omega$ について減少である。
</Lemma>

<Proof of="lem-supermodular">
**$W=\emptyset$ の場合.** 辺 $e=\{x,y\}$ を $\omega$（$e\notin\omega$）に加えると、$x,y$ が $\omega$ で別のクラスターなら $k$ は $1$ 減り、同じなら変わりません。$\omega\subset\eta$、$e\notin\eta$ のとき、$\eta$ で $x,y$ が別のクラスターなら $\omega$ でも別なので（$\omega$ のクラスターは $\eta$ のクラスターを細分する）、
$$
k(\omega)-k(\omega\cup\{e\})\ge k(\eta)-k(\eta\cup\{e\}).
$$
有限集合 $F$（$F\cap\eta=\emptyset$）の辺を 1 本ずつ加えてこれを足し合わせると、$k(\omega)-k(\omega\cup F)\ge k(\eta)-k(\eta\cup F)$。$\omega\wedge\eta\subset\eta$、$F=\omega\setminus\eta$ に使うと、$(\omega\wedge\eta)\cup F=\omega$、$\eta\cup F=\omega\vee\eta$ なので
$$
k(\omega\wedge\eta)-k(\omega)\ge k(\eta)-k(\omega\vee\eta),
$$
これが主張です。

**一般の $W\ne\emptyset$.** $W$ の点を結ぶ新しい辺の集合 $F_W$（$W$ 上の全域木）を加えたグラフを考えると、$W$ と交わるクラスターはちょうど 1 つにまとまるので $k_W(\omega)=k(\omega\cup F_W)-1$ です。$\omega\mapsto\omega\cup F_W$ は $\vee$ と $\wedge$ を保つので、$W=\emptyset$ の場合に帰着します。

**$N_W$ の単調性.** 辺を加えるとクラスターは合併するだけなので、$W$ と交わるクラスターの個数は増えません。
</Proof>

<Theorem id="thm-rc-fkg" title="ランダムクラスター測度の FKG 不等式">
$q\ge1$ ならば、$\varphi^W_{p,q}$ は FKG の束条件（<Ref to="physics/constructive-qft/correlation-inequalities#def-fkg-lattice-condition" />）を満たす。したがって増加関数 $f,g$ について $\varphi^W_{p,q}(fg)\ge\varphi^W_{p,q}(f)\varphi^W_{p,q}(g)$。
</Theorem>

<Proof of="thm-rc-fkg">
$\varphi=\varphi^W_{p,q}$ はすべての点で正です（$p_e\in(0,1)$）。$|\omega\vee\eta|+|\omega\wedge\eta|=|\omega|+|\eta|$ は辺ごとに成り立つので、$p_e$ の因子は束条件の両辺で打ち消し合い、
$$
\frac{\varphi(\omega\vee\eta)\varphi(\omega\wedge\eta)}{\varphi(\omega)\varphi(\eta)}=q^{k_W(\omega\vee\eta)+k_W(\omega\wedge\eta)-k_W(\omega)-k_W(\eta)}\ge1
$$
です（<Ref to="lem-supermodular" />、$q\ge1$）。あとは FKG 不等式（<Ref to="physics/constructive-qft/correlation-inequalities#thm-fkg" />）です。
</Proof>

測度 $\mu,\nu$ について、すべての増加関数 $f$ で $\mu(f)\le\nu(f)$ のとき $\mu\preceq\nu$（$\nu$ が $\mu$ を**確率的に優越する**）と書きます。比較には次の単純な観察を繰り返し使います。

<Lemma id="lem-tilt" title="増加関数で傾けると大きくなる">
$\mu$ が FKG の束条件を満たし、$g>0$ が増加関数ならば、$\mu_g=g\mu/\mu(g)$ について $\mu\preceq\mu_g$。$g>0$ が減少関数なら $\mu_g\preceq\mu$。
</Lemma>

<Proof of="lem-tilt">
増加関数 $f$ について、FKG 不等式から $\mu(fg)\ge\mu(f)\mu(g)$、すなわち $\mu_g(f)\ge\mu(f)$。$g$ が減少なら FKG 不等式の後半（<Ref to="physics/constructive-qft/correlation-inequalities#thm-fkg" />）で不等号が逆になります。
</Proof>

次の定理のねらいは、ランダムクラスター測度を、扱いやすい独立浸透で上下から挟むことです。上側の $P_p$ は重み $q^{k_W}$ を外したもの、下側の $P_{p'}$ は、他の辺の状況によらず 1 本の辺が開く条件付き確率の最小値 $p'$ をすべての辺に使ったものと読めます（演習の 1 辺の計算がその様子です）。

<Theorem id="thm-comparison" title="比較不等式">
$q\ge1$ とし、$p'_e=\dfrac{p_e}{p_e+q(1-p_e)}$ とおく。

(i) $P_{p'}\preceq\varphi^W_{p,q}\preceq P_p$（$P_p$ は独立な Bernoulli 浸透）。

(ii) $\varphi^\emptyset_{p,q}\preceq\varphi^W_{p,q}$（自由 $\preceq$ 配線）。

(iii) $p_e\le\tilde p_e$（すべての $e$）ならば $\varphi^W_{p,q}\preceq\varphi^W_{\tilde p,q}$。
</Theorem>

<Proof of="thm-comparison">
**(i) 上側.** $\varphi^W_{p,q}$ は $P_p$ を $g=q^{k_W}$ で傾けたものです。辺を加えると $k_W$ は増えない（クラスターが合併するだけ）ので $g$ は減少で、$P_p$ は束条件を満たす（$q=1$ の場合）ので <Ref to="lem-tilt" /> から $\varphi^W_{p,q}\preceq P_p$。

**(i) 下側.** $p_e/(1-p_e)=q\,p'_e/(1-p'_e)$ から
$$
\prod_ep_e^{\omega_e}(1-p_e)^{1-\omega_e}q^{k_W(\omega)}=\Bigl(\prod_e\frac{1-p_e}{1-p'_e}\Bigr)\prod_e{p'_e}^{\omega_e}(1-p'_e)^{1-\omega_e}\;q^{|\omega|+k_W(\omega)}.
$$
辺を 1 本加えると $|\omega|$ は $1$ 増え、$k_W$ は高々 $1$ 減るので $|\omega|+k_W(\omega)$ は増加です。よって $\varphi^W_{p,q}$ は $P_{p'}$ を増加関数で傾けたもので、$P_{p'}\preceq\varphi^W_{p,q}$。

**(ii)** $q^{k_W}=q^{k}q^{-N_W}$ で、$N_W$ は減少（<Ref to="lem-supermodular" />）なので $q^{-N_W}$ は増加です。$\varphi^W_{p,q}$ は $\varphi^\emptyset_{p,q}$（<Ref to="thm-rc-fkg" /> で束条件を満たす）を $q^{-N_W}$ で傾けたものです。

**(iii)** $\varphi^W_{\tilde p,q}$ は $\varphi^W_{p,q}$ を $\prod_e\bigl(\frac{\tilde p_e(1-p_e)}{p_e(1-\tilde p_e)}\bigr)^{\omega_e}$ で傾けたもので、各因子の底は $1$ 以上なのでこれは増加です。
</Proof>

## 5. Ising 模型への応用

$\mathbb{Z}^d$ の最近接 Ising 模型（$h=0$）で $p=1-e^{-2\beta J}$ とおきます。§3 の表示と <Ref to="thm-comparison" /> (i) から、すべての有限集合 $\Lambda\ni0$ について
$$
P_{p'}\bigl(0\leftrightarrow\partial^{\mathrm{out}}\Lambda\bigr)\le\langle\sigma_0\rangle^+_\Lambda\le P_p\bigl(0\leftrightarrow\partial^{\mathrm{out}}\Lambda\bigr),\qquad p'=\frac{p}{2-p},
$$
が成り立ちます（$P$ はグラフ $(V,E_\Lambda)$ 上の独立浸透）。

次の命題は、高温で磁化が消えることを、原点から出る長い開いた道の数え上げで示します。同じ点を二度通らない道（自己回避路）の個数は、1 歩ごとに高々 $2d-1$ 通りの選び方しかないので、$(2d-1)p<1$ なら長い道が開いている確率は急速に小さくなります。

<Proposition id="prop-high-temperature" title="経路の数え上げによる高温評価">
$(2d-1)(1-e^{-2\beta J})<1$ ならば、$\Lambda_n=\{-n,\dots,n\}^d$ について $\langle\sigma_0\rangle^+_{\Lambda_n}\le\frac{2d}{2d-1}\bigl((2d-1)p\bigr)^{n+1}\to0$ で、自発磁化 $m^*(\beta)$ は $0$ である。
</Proposition>

<Proof of="prop-high-temperature">
$0\leftrightarrow\partial^{\mathrm{out}}\Lambda_n$ なら、$0$ から出て $\Lambda_n$ の外に達する開いた自己回避路があり、その長さは $n+1$ 以上です。その最初の $n+1$ 歩は長さ $n+1$ の自己回避路で、その個数は $2d(2d-1)^n$ 以下（最初の 1 歩は $2d$ 通り、以後は直前の点に戻らないので $2d-1$ 通り以下）、各路がすべて開である確率は $p^{n+1}$ です。和の上界をとれば主張の評価になります。$m^*=\lim_n\langle\sigma_0\rangle^+_{\Lambda_n}$（<Ref to="mathematics/gibbs-measures/ising-peierls#thm-plus-state" />）です。
</Proof>

$d=2$ ではこの条件は $\beta J<\frac12\log\frac32\approx0.203$ で、Dobrushin の判定 $\beta J<\operatorname{artanh}\frac14\approx0.255$（<Ref to="mathematics/gibbs-measures/dlr-gibbs-measures#ex-ising-dobrushin" />）より少し弱い結果です。経路の数え上げを独立浸透の臨界値 $p_c(\mathbb{Z}^d)$ の知識に置き換えれば、$p<p_c$ で磁化が消えることがわかります。

次の命題は、低温で磁化が残ることを、Peierls の議論と同じ輪郭の数え上げで示します。原点が境界とつながらないとき、原点のクラスターは閉じた辺の閉曲線で囲まれているので、その長さごとに個数を数えて、すべて閉じている確率を足し上げます。

<Proposition id="prop-low-temperature-fk" title="輪郭の数え上げによる低温評価">
$d=2$ とし、$u=e^{-2\beta J}$、$y=9\bigl(\frac{2u}{1+u}\bigr)^2<1$ とおく。正方形 $\Lambda=\{-n,\dots,n\}^2$ について
$$
\langle\sigma_0\rangle^+_\Lambda\ge1-\sum_{j\ge2}2j\,y^j .
$$
特に $\beta J\ge\frac12\log11$ ならば $\langle\sigma_0\rangle^+_\Lambda\ge11/18$ で、$m^*(\beta)>0$。
</Proposition>

<Proof of="prop-low-temperature-fk">
$1-p'=\frac{2(1-p)}{2-p}=\frac{2u}{1+u}$ です。$P_{p'}$ のもとで $0\not\leftrightarrow\partial^{\mathrm{out}}\Lambda$ とし、$C$ を $0$ の開いたクラスターとすると $C\subset\Lambda$ で、$C$ から外へ出る辺はすべて閉です。<Ref to="mathematics/gibbs-measures/ising-peierls#thm-peierls" /> の証明 (1) と同じく $C$ の穴を埋めた $B$ をとると、$B$ は有限・連結で補集合も連結、$0\in B\subset\Lambda$ で、$\partial_eB$ の辺はすべて $C$ から $C$ の外への辺なので閉です。独立性から $P_{p'}(\partial_eB\text{ がすべて閉})=(1-p')^{|\partial_eB|}$ で、同証明 (3) の数え上げ（長さ $\ell$ の輪郭は $\ell\,3^\ell$ 個以下、$\ell$ は $4$ 以上の偶数）から
$$
P_{p'}(0\not\leftrightarrow\partial^{\mathrm{out}}\Lambda)\le\sum_{\ell\ge4,\ \ell\text{ 偶}}\ell\,3^\ell(1-p')^\ell=\sum_{j\ge2}2j\,y^j .
$$
上の下界の不等式と合わせれば主張です。$\beta J\ge\frac12\log11$ なら $u\le1/11$、$\frac{2u}{1+u}\le\frac16$、$y\le\frac14$ で、和は $\frac7{18}$ 以下です。
</Proof>

Peierls の議論（<Ref to="mathematics/gibbs-measures/ising-peierls#thm-peierls" />）はスピンの反転という Ising に固有の写像を使いましたが、この証明は独立浸透との比較だけを使うので、$q\ge2$ の Potts 模型にも同じ形で通ります（$p'=p/(p+q(1-p))$ にするだけ）。

## 6. 長距離浸透と 1/r² 模型

$\mathbb{Z}$ の任意の 2 点 $x<y$ を辺で結び、辺 $\{x,y\}$ を確率 $p_{y-x}$ で独立に開く**長距離浸透**を考えます（$p_r\in[0,1)$）。$C(0)$ を原点の開いたクラスター、$\theta=P(|C(0)|=\infty)$ を**浸透確率**とします。1 次元の最近接浸透（$p_r=0$、$r\ge2$）では $p_1<1$ なら閉じた辺が必ず現れて $\theta=0$ ですが、長い辺は閉じた辺を飛び越えます。どのくらい長い辺があれば浸透が起きるかが問題です。

<Lemma id="lem-zero-one" title="Kolmogorov の 0-1 法則">
$(\xi_i)_{i\in\mathbb{N}}$ を独立な確率変数の列とし、事象 $A$ がすべての $n$ について $\sigma(\xi_i:i>n)$ に属するとする。このとき $P(A)\in\{0,1\}$。
</Lemma>

<Proof of="lem-zero-one">
$\mathcal{A}_n=\sigma(\xi_1,\dots,\xi_n)$ とすると、独立性から $\mathcal{A}_n$ と $\sigma(\xi_i:i>n)$ は独立で、$A$ は $\bigcup_n\mathcal{A}_n$ のすべての元と独立です。$\bigcup_n\mathcal{A}_n$ は π-系で $\sigma(\xi_i:i\in\mathbb{N})\ni A$ を生成するので、2 つの測度 $B\mapsto P(A\cap B)$ と $B\mapsto P(A)P(B)$ は π-系の上で一致し、測度の一意性（<Ref to="mathematics/real-analysis/abstract-measures#thm-uniqueness-measures" />）から生成された σ-加法族全体で一致します。$B=A$ として $P(A)=P(A)^2$。
</Proof>

次の命題の方針は、ある点 $k+\frac12$ をまたぐ開いた辺が 1 本もない場所が、原点の左右にそれぞれ無限個あれば、原点のクラスターはその間に閉じ込められる、というものです。そのような場所が無限個あることを、0-1 法則と確率の下からの評価で示します。

<Proposition id="prop-lr-no-percolation" title="短い辺が支配的なら浸透しない">
$\sum_{r\ge1}r\,p_r<\infty$ ならば $\theta=0$。
</Proposition>

<Proof of="prop-lr-no-percolation">
$k\in\mathbb{Z}$ について $A_k$ を「$x\le k<y$ となる開いた辺 $\{x,y\}$ がない」事象とします。長さ $r$ の辺で点 $k+\frac12$ をまたぐものは $r$ 本なので、独立性から
$$
P(A_k)=\prod_{r\ge1}(1-p_r)^r=:\delta .
$$
$\sum_rr\,p_r<\infty$ から $p_r\to0$ で、$p_r\le\frac12$ では $-\log(1-p_r)\le2p_r$ なので $\sum_rr(-\log(1-p_r))<\infty$、各因子は正なので $\delta>0$ です。

辺に番号を付けて独立な列と見ます。有限個の辺は有限個の $A_k$ にしか影響しないので、事象 $A^+=\{A_k\text{ が無限個の }k\ge0\text{ で起きる}\}$ は <Ref to="lem-zero-one" /> の仮定を満たします。また
$$
P(A^+)=\lim_{n\to\infty}P\Bigl(\bigcup_{k\ge n}A_k\Bigr)\ge\limsup_kP(A_k)=\delta>0
$$
なので $P(A^+)=1$。同様に負の側でも、$A_k$ が無限個の $k<0$ で起きる確率は $1$ です。両方が起きれば、$C(0)$ は $k_-<0\le k_+$ となる $A_{k_\pm}$ により $\{k_-+1,\dots,k_+\}$ に閉じ込められて有限です。
</Proof>

<Corollary id="cor-lr-ising" title="短距離的な 1 次元 Ising 模型に自発磁化はない">
$\mathbb{Z}$ 上の外場 $0$ の Ising 模型で結合を $J_{xy}=J(|x-y|)\ge0$、$\sum_rr\,J(r)<\infty$ とする。このとき任意の $\beta$ で $\lim_n\langle\sigma_0\rangle^+_{\{-n,\dots,n\}}=0$。
</Corollary>

<Proof of="cor-lr-ising">
§3 の表示を長距離の結合に使います。$+$ 境界条件は、$\Lambda=\{-n,\dots,n\}$ の外の点をすべて $W$ に含めることに当たり、外の点どうしの結合は分布に影響しないので、外の点を 1 点 $w$ にまとめ、$x\in\Lambda$ と $w$ を結ぶ辺を「$x$ から外への辺のどれかが開」と読み替えられます。<Ref to="thm-comparison" /> (i) から、$p_r=1-e^{-2\beta J(r)}\le2\beta J(r)$ の独立長距離浸透について
$$
\langle\sigma_0\rangle^+_\Lambda\le P_p\bigl(0\text{ が }\Lambda\text{ の外の点とつながる}\bigr).
$$
右辺の事象は $n$ について減少し、その共通部分は $\{C(0)\text{ は有界でない}\}=\{|C(0)|=\infty\}$ なので、右辺は $\theta$ に収束します。$\sum_rr\,p_r\le2\beta\sum_rr\,J(r)<\infty$ なので <Ref to="prop-lr-no-percolation" /> から $\theta=0$。（無限個の辺をもつグラフへの比較不等式の適用は、$\Lambda$ の外への辺の確率を有限個の辺にまとめれば有限グラフの場合に帰着します。）
</Proof>

この境界にあるのが $p_r\sim\beta/r^2$ です。$\sum_rr\cdot r^{-2}$ は対数発散するので <Ref to="prop-lr-no-percolation" /> は使えず、実際 $\beta$ と短い辺の確率しだいで浸透が起こりえます。

条件 $\theta^2\beta^+\ge1$ に $\theta$ の二乗が現れるのは、近くにある 2 つの区間の両方が無限クラスターに属する必要があり、その割合がそれぞれ $\theta$ だからです（骨格は後の注意で述べます）。この定理の証明は本章では与えず、主張と骨格だけを述べます。

<Theorem id="thm-accn" title="1/r² 浸透の二分律（Aizenman–Newman, ACCN の主張）">
1 次元の長距離浸透で $\sup_rp_r<1$、$\beta^+=\limsup_{r\to\infty}r^2p_r<\infty$ とする。

(i) $\beta^+\le1$ ならば、短い辺の確率がどれほど $1$ に近くても $\theta=0$（Aizenman–Newman 1986）。

(ii) $\theta=0$ または $\theta^2\beta^+\ge1$（Aizenman–Newman 1986）。

(iii) $\liminf_rr^2p_r>1$ のとき、$p_1$ が十分 $1$ に近ければ $\theta>0$（Newman–Schulman 1986）。

したがって、$p_r$ をパラメータについて連続かつ単調に増やしていく族では、浸透確率は臨界点で $0$ から $1/\sqrt{\beta^+}$ 以上へ**跳ぶ**。Aizenman–Chayes–Chayes–Newman（1988）は、同じ二分律をランダムクラスター表示を通じて結合 $J(r)\sim J^+/r^2$ の 1 次元 Ising・Potts 模型に移し、自発磁化が臨界点で不連続である（**Thouless 効果**）ことを示した。
</Theorem>

<Remark id="rem-accn-skeleton" title="証明の骨格と所在">
証明は本章では与えません（原論文を参照）。骨格は次のとおりです。

1. **ブロックの繰り込み.** 長さ $L$ の区間 2 つが距離 $L$ 程度離れているとき、その間の辺の本数の期待値は $\sum_{x,y}\beta/(y-x)^2$ の程度で、$L$ によらない定数になります（$1/r^2$ がスケール不変であることの表れ）。各区間で無限クラスターに属する点の割合が $\theta$ なら、無限クラスターどうしを結ぶ辺の数はおよそ $\theta^2\beta$ の定数倍です。
2. **連結の伝播と切断.** $\theta^2\beta^+<1$ だと、スケールを $L\to2L\to4L\to\cdots$ と上げるごとに「隣のブロックの無限クラスターとつながりそこなう」確率が一定以上残り、それが無限回起きて無限クラスターが切れる、という矛盾を導きます。この評価には、有限体積の連結確率についての再帰的な不等式を使います。
3. **Ising への移植.** $q\ge1$ のランダムクラスター測度は独立浸透の比較（<Ref to="thm-comparison" />）と FKG（<Ref to="thm-rc-fkg" />）を満たすので、独立性の代わりにこれらの単調性を使って同じ議論が通ります。

Dyson（1969）は $J(r)\sim r^{-s}$（$1<s<2$）の 1 次元 Ising 模型に相転移があることを示し、Fröhlich–Spencer（1982）は $s=2$ で相転移と Thouless 効果を示しました。$s>2$ では <Ref to="cor-lr-ising" /> により相転移はありません。
</Remark>

## 7. 演習

<Exercise id="exr-single-edge" difficulty="易">
辺 $e=\{x,y\}$ 以外の辺の配置 $\xi$ を固定したとき、$\varphi^\emptyset_{p,q}(\omega_e=1\mid\xi)$ は、$\xi$ で $x\leftrightarrow y$ なら $p_e$、そうでなければ $p_e/(p_e+q(1-p_e))$ であることを示してください。これと <Ref to="thm-comparison" /> (i) の $p'$ との関係を説明してください。
<Solution>
$\omega_e=1$ と $0$ の重みの比は $\frac{p_e}{1-p_e}q^{k(\xi\cup\{e\})-k(\xi)}$ です。$\xi$ で $x\leftrightarrow y$ なら $k$ は変わらず比は $p_e/(1-p_e)$、条件付き確率は $p_e$。そうでなければ $k$ が $1$ 減り、比は $p_e/(q(1-p_e))$、条件付き確率は $p_e/(p_e+q(1-p_e))=p'_e$ です。$q\ge1$ では $p'_e\le p_e$ なので、各辺の条件付き確率はつねに $[p'_e,p_e]$ に入ります。これは比較 $P_{p'}\preceq\varphi\preceq P_p$ の「1 辺ずつの」版で、Holley の不等式を使えばここから <Ref to="thm-comparison" /> (i) を導くこともできます。
</Solution>
</Exercise>

<Exercise id="exr-potts-two-sites" difficulty="易">
2 頂点 1 辺のグラフで、$q$ 状態 Potts 模型の $\pi(\sigma_x=\sigma_y)$ を直接計算し、<Ref to="cor-correlations-connectivity" /> の右辺と一致することを確かめてください。
<Solution>
直接計算では、一致する $q$ 通りの配位の重みが $1$、一致しない $q(q-1)$ 通りの重みが $e^{-\beta J}=1-p$ なので $\pi(\sigma_x=\sigma_y)=\frac{q}{q+q(q-1)(1-p)}=\frac1{1+(q-1)(1-p)}$。ランダムクラスター側では $\varphi(\omega_e=1)=\frac{pq}{pq+(1-p)q^2}=\frac{p}{p+q(1-p)}$ で、
$$
\frac1q+\Bigl(1-\frac1q\Bigr)\frac{p}{p+q(1-p)}=\frac{p+q(1-p)+(q-1)p}{q(p+q(1-p))}=\frac{q}{q(p+q(1-p))}=\frac1{1+(q-1)(1-p)}.
$$
</Solution>
</Exercise>

<Exercise id="exr-swendsen-wang" difficulty="標準">
<Ref to="thm-edwards-sokal" /> の (iii)(iv) を交互に使う Markov 連鎖（$\sigma$ から (iv) で $\omega$ を引き、$\omega$ から (iii) で新しい $\sigma$ を引く）が Potts 模型 $\pi^W$ を定常分布にもつことを示してください（**Swendsen–Wang 法**）。
<Solution>
$\sigma\sim\pi^W$ とすると、(i) から $\sigma$ は結合分布 $\mathbb{P}$ の $\sigma$-周辺分布に従い、(iv) で引いた $\omega$ との組 $(\sigma,\omega)$ は $\mathbb{P}$ に従います。したがって $\omega$ は $\mathbb{P}$ の $\omega$-周辺分布に従い、(iii) で引いた $\sigma'$ との組 $(\sigma',\omega)$ も $\mathbb{P}$ に従います。よって $\sigma'\sim\pi^W$ です。1 ステップでクラスター全体の色を一度に変えられるので、臨界点の近くでも 1 点ずつの更新より速く混ざることが知られています。
</Solution>
</Exercise>

## 参考文献

- G. Grimmett, *The Random-Cluster Model*, Springer, 2006 — 第 1 章（Edwards–Sokal 結合）、第 2・3 章（FKG と比較不等式）。
- H. Duminil-Copin, "Lectures on the Ising and Potts models on the hypercubic lattice", 2017（講義録）. [arXiv:1707.00520](https://arxiv.org/abs/1707.00520)
- R. G. Edwards, A. D. Sokal, "Generalization of the Fortuin–Kasteleyn–Swendsen–Wang representation and Monte Carlo algorithm", *Physical Review D* 38 (1988), 2009–2012.
- M. Aizenman, C. M. Newman, "Discontinuity of the percolation density in one dimensional 1/|x−y|² percolation models", *Communications in Mathematical Physics* 107 (1986), 611–647.
- M. Aizenman, J. T. Chayes, L. Chayes, C. M. Newman, "Discontinuity of the magnetization in one-dimensional 1/|x−y|² Ising and Potts models", *Journal of Statistical Physics* 50 (1988), 1–40.
- C. M. Newman, L. S. Schulman, "One dimensional 1/|j−i|^s percolation models: the existence of a transition for s≤2", *Communications in Mathematical Physics* 104 (1986), 547–571.
