======TAP 방정식======
복제본 트릭 없이 스핀 유리의 성질을 이해하기 위해 사울레스, 앤더슨, 그리고 팔머가 개발한 이론 틀.

====== 셰링턴-커크패트릭 모형 ======

[[물리:셰링턴-커크패트릭 모형|셰링턴-커크패트릭 모형]]에 대해서 [[신경망:공동_방법|공동 방법(cavity method)]]으로 유도할 수 있다.

기본적인 아이디어는, $N$개의 스핀(spin)으로 이루어진 계에서 $N+1$개의 스핀으로 이루어진 계로 변화시키면서도
그러한 두 계에서 얻는 관측 가능량(observables)이 열역학적 극한에서 차이가 없어야 한다는 점을 이용하는 것이다.

===== 유도 과정 =====
==== 국소장 ====

관측 가능량 중 하나로서 '국소적 자화량(local magnetization)'을 $m_0 = \langle s_0 \rangle$라고 표기한다면, 이는 한 곳에 대한 열적 평균(thermal average)을 내린 값이다. 이때, 원래 $N$개의 스핀으로 이루어진 계에 $s_0$를 추가하는 상황을 생각해보자. ($s_0$는 위치 $0$에 추가된 스핀을 의미하는 관습적인 표기법이다.)

$s_0$를 추가함에 따라 $N+1$개의 스핀들로 이루어진 계로 변화하였다. 추가된 스핀의 국소적 자화량은 $m_0=\langle s_0 \rangle$이다. 이번 글에서의 최종 목표는, 이러한 국소적 자화량이 TAP 방정식을 만족함을 보이는 것이다. 국소적 자화량에 대한 평균인 $\langle s_0 \rangle$을 얻기 위해서는, 국소적 스핀에 대한 분포함수 $P(s_0)$를 얻으면 된다. 이어지는 논의에서는 편의상 외부 자기장은 없다고 가정하여 $h_i=0$이라고 하자.

스핀 $s_0$가 추가되어 계의 크기가 $N+1$이 되었을 때의 확률분포는 다음처럼 기술된다:
\[ P \left(\vec{s}, s_0\right) \propto \exp\left[ -\beta H_N (\vec{s}) + \beta \sum_j J_{0j} s_0 s_j \right] = e^{\beta s_0 \tilde{h}_0} \exp\left[ -\beta H_N (\vec{s}) \right]. \]

여기에서 $\vec{s} \equiv (s_1, \ldots, s_N)$이고, $\tilde{h}_0$는 다른 스핀들이 $s_0$의 자리에 미치는 국소장으로서 다음처럼 정의되고
$$ \tilde{h}_0 \equiv \sum_{j=1}^N J_{0j} s_j.$$
공동 장(cavity field)라고도 불린다.

$s_0$ 없이 $N$개의 스핀 $\vec{s}$에 대한 분포 함수는 위 식에서 $s_0=0$을 대입하면 된다:
\[ P \left(\vec{s}, s_0=0\right) \propto \exp\left[-\beta H_N (\vec{s}) \right]. \]

/*
스핀 $s_0$가 추가되므로 그의 이웃인 스핀 $s_j$들의 상호작용을 고려해야 한다.
$s_0$와 $\tilde{h}_0 $의 결합 분포 $P(s_0,\tilde{h}_0)$를 구한 후에, 그를 통해 $s_0$와 $\tilde{h}_0$의 열적 평균 값을 구하면 될 것이다.
$s_0$에 의해 스핀 $N$개에서 $N+1$개로 변한다면, 에너지 $E^{(N)}$에서 $E^{(N+1)}$로 다음과 같이 변한다:

\begin{eqnarray*}
s_0 = 1, &~~& E^{(N+1)}=E^{(N)}-\beta \tilde{h}_0 \\
s_0 = -1, &~~& E^{(N+1)}=E^{(N)}+\beta \tilde{h}_0.
\end{eqnarray*}

즉 우리가 구하고자 하는 $P(s_0,\tilde{h}_0)$는 다음과 같이 표현될 것이다:

$$ P(s_0,\tilde{h}_0) \propto  \ e^{\beta \tilde{h}_0 s_0} P(\tilde{h}_0 \setminus s_0).$$

여기에서 $P(\tilde{h}_0 \setminus s_0)$는 $s_0$가 식에 포함되지 않은 확률 분포를 말한다. 더 자세하게 표현하자면, $P(\tilde{h}_0 \setminus s_0)$는 ($s_0$가 없을 때) $\tilde{h}_0$에 대한 확률 분포의 식이며, $e^{\beta \tilde{h}_0 s_0}$은 $N+1$의 계를 기술하기 위해 $P(\tilde{h}_0 \setminus s_0)$에 포함되었다.

이제 $ P(\tilde{h}_0 \setminus s_0)$의 식을 나타낼 수 있다면, 공동 방법의 기본적인 아이디어를 통해서 TAP 방정식을 얻는 것이 가능하다.
*/
==== 추가 스핀이 없을 때 국소장의 분포 ====

$s_0$가 없을 때 $\tilde{h}_0$의 확률분포를 $P\left(\tilde{h}_0 \setminus s_0\right)$라 하자.

만약 격자 구조가 아래와 같은 베테 격자 구조라면, $P\left(\tilde{h}_0 \setminus s_0\right)$에 대해 중심 극한 정리를 적용할 수 있다. 왜냐하면 위치 $0$의 $s_0$에 연결된 각각의 이웃 위치 $j$에 대한 $J_{0j} s_j$는 '$s_0$의 정보가 제외된다면' 서로 독립적이기 때문이다. 

{{:물리:tap_bethe_edited.png?300|}}

(위의 베테 격자 그림에서 $s_0$를 감싸는 공동(cavity)을 생각하고 $s_0$가 없는 상황을 떠올리면 보다 이해하기 쉽다. 그 상황에서 스핀의 개수는 $N$개 이다.)

이때 서로 다른 위치에 대한 상관관계가 약한 [[물리:셰링턴-커크패트릭_모형|셰링턴-커크패트릭 모형]]에 대해서도 이러한 가정을 도입할 수 있다고 하자.
즉, 독립적이고 동일하게 분포된(independent and identically distributed) $J_{0j} s_j$에 대해서, $s_0$가 식에 포함되지 않는 $\tilde{h}_0$의 분포는 ($N$이 증가함에 따라) 정규 분포를 따를 것이다. 즉 다음과 같은 형태로 표현 가능하다:

$$ P\left(\tilde{h}_0 \setminus s_0\right)=\frac{1}{\sqrt{2\pi V_0}} \exp \left[ -\frac{(\tilde{h}_0 - \langle \tilde{h}_0 \rangle_{\setminus 0})^2}{2V_0}\right].$$

이때 $h$는 공동 장의 평균으로서 다음처럼 정의되고:
\[\langle \tilde{h}_0 \rangle_{\setminus 0} = \sum_{j=1}^N J_{0i} \langle s_i \rangle_{\setminus 0}\]
$V_0$는 국소장의 분산이다. 아래 첨자 $\setminus 0$는 스핀 $s_0$를 제외한 $N$개의 스핀 $\vec{s}$를 가지는 계의 확률분포에 대해 평균하고 있음을 의미한다.
분산 $V_0$를 식으로 표현하면 아래와 같다.

$$ 
\begin{align}

V_0 &= \langle (\tilde{h}_0)^2\rangle_{\setminus 0} - \langle 
\tilde{h}_0 \rangle_{\setminus 0} ^2
\\
&=\sum_{j,k}J_{0j}J_{0k} \left(\langle s_j s_k\rangle_{\setminus 0} \ \ - \ \langle s_j\rangle_{\setminus 0} \langle s_k\rangle_{\setminus 0}  \right).
\end{align}$$

이때 결합 상수의 분포함수에 주목하자:

$$ P(J_{ij}) = \frac{1}{J}\sqrt{\frac{N}{2\pi}} \exp\left\{ - \frac{N}{2J^2}\left( J_{ij} - \frac{J_0}{N} \right)^2 \right\} .$$ 

위의 식으로부터 $J_{ij}$의 기댓값은 $J_0/N$이며 분산은 $J^2/N$임을 알 수 있다. 편의상 $J_0 \equiv 0$이고 $J\equiv 1$이라고 놓자. 그러면 $J_{ij}$의 크기는 $1/\sqrt{N}$ 정도이다.

또 [[신경망:공동_방법|공동 방법]]의 가정은 $N\to\infty$의 극한에서 아래의 비선형 감수율(nonlinear susceptibility)이 유한하다는 것인데
\[ \chi_\text{nl} = \frac{1}{N} \sum_{i<j} \left( \langle s_i s_j \rangle - \langle s_i \rangle \langle s_j \rangle \right)^2 \]
이는 $\langle s_i s_j \rangle_{\setminus 0} - \langle s_i \rangle_{\setminus 0} \langle s_j \rangle_{\setminus 0}$의 크기가 $1/\sqrt{N}$ 정도임을 의미한다.

그래서 $V_0$를 계산할 때, $N(N-1)$개의 서로 상관관계가 없고 무작위적으로 부호를 갖는 $j\ne k$의 항들의 합은 $N \to \infty$에서 $0$이 된다.

따라서 $N \to \infty$에서 $j=k$인 경우만이 $V_0$에 기여한다는 것을 알 수 있다. 그 결과는 다음과 같이 표현된다:

$$ V_0 \approx \sum_j J_{0j}^2 \left(1-\langle s_j \rangle^2_{\setminus 0} \right) \approx 1 - \frac{1}{N} \sum_i \langle s_i \rangle_{\setminus 0}^2 \approx 1 - \overline{\langle s_i \rangle_{\setminus 0}^2} = 1-q. $$
여기에서 $J_{ij} \sim 1/\sqrt{N}$임을 이용했고, 공간적인 대칭성을 이용해서 모든 지점에 대한 합 $\sum_i$를 한 지점에서의 무질서 평균으로 대체했다. $q \equiv \overline{\langle s_i \rangle_{\setminus 0}^2}$는 에드워즈-앤더슨 질서 변수(Edwards-Anderson order parameter)이다.

사실 마찬가지 논법으로
\[ \Big\langle \left( \tilde{h}_0 - \langle \tilde{h}_0 \rangle_{\setminus 0} \right)^4 \Big\rangle_{\setminus 0} = \sum_{ijkl} J_{0i} J_{0j} J_{0k} J_{0l} \langle (s_i - \langle s_i \rangle_{\setminus 0}) (s_j - \langle s_j \rangle_{\setminus 0}) (s_k - \langle s_k \rangle_{\setminus 0}) (s_l - \langle s_l \rangle_{\setminus 0}) \rangle_{\setminus 0} \]
을 고려해보면 $i=j$와 $k=l$, $i=k$와 $j=l$, 그리고 $i=l$과 $j=k$의 세 경우만 $N\to\infty$에서 살아남고 $i=j=k=l$은 살아남지 못하기 때문에 최종 결과는 $\langle \left( \tilde{h}_0 - \langle \tilde{h}_0 \rangle_{\setminus 0} \right)^4 \rangle_{\setminus 0} \approx 3(1-q)^2$인데 이는 [[수학:윅의_정리|윅의 정리]]와 비슷하다. 일반적으로
\[\langle (\tilde{h}_0)^{2p} \rangle_{\setminus 0} = (2p-1)!! (1-q)^p\]
로서 이는 다시 한번
\[ P\left(\tilde{h}_0 \setminus s_0\right) \propto \exp \left[ - \frac{(\tilde{h}_0-h)^2}{2(1-q)} \right] \]
임을 의미한다.

앞의 식을 활용한다면 아래처럼 표현할 수도 있다:
\[ P\left(\tilde{h}_0 \setminus s_0\right) = \int d\vec{\sigma} P\left(\vec{\sigma}, s_0=0 \right) \delta\left( \sum_{j=1}^N J_{0j} s_j - \tilde{h}_0 \right). \]

/*이와 같이 국소장의 분산 $V_0$을 얻었으므로 $ P(S_0) \propto \int d \tilde{h}_0 \ e^{\beta \tilde{h}_0 S_0} P(\tilde{h}_0 ＼S_0)$의 관계를 이용하여 $m_0$를 구하면 된다.
*/
==== 추가된 스핀의 기댓값 ====

이제 스핀 $s_0$를 추가하고 $N+1$개의 스핀을 가지는 계에 대해 $P(s_0)$를 구해보자. 이를 위해 먼저 $P\left( \tilde{h}_0, s_0 \right)$를 구한 다음 주변화(marginalization)를 시행할 것이다.
\begin{eqnarray*}
P\left( \tilde{h}_0, s_0 \right) &=& \int d\vec{s} P \left(\vec{s}, s_0 \right) \delta \left( \sum_{j=1}^N J_{0j} s_j - \tilde{h}_0 \right)\\
&=& \int d\vec{s} e^{\beta s_0 \tilde{h}_0} P \left(\vec{s}, s_0=0 \right) \delta \left( \sum_{j=1}^N J_{0j} s_j - \tilde{h}_0 \right)\\
&=& e^{\beta s_0 \tilde{h}_0} P \left(\tilde{h}_0 \setminus s_0 \right).
\end{eqnarray*}
여기에서 $J_{0j}$의 분산이 $1/N$임을 활용했다.

이제 $P(s_0)$는 다음과 같이 간단한 [[수학:윅의_정리|가우스 적분]]으로 표현된다:
/*
$$
\begin{align}
P(s_0) & \propto \int d \tilde{h}_0 \ \exp[{\beta \tilde{h}_0 s_0}] P(\tilde{h}_0 \setminus s_0) \\
& = \int d\tilde{h}_0 \frac{1}{\sqrt{2\pi V_0}} \exp[{\beta \tilde{h}_0 s_0}] \exp \left[ -\frac{(\tilde{h}_0 - \langle \tilde{h}_0\rangle _{\setminus0})^2}{2V_0}\right]\\
& = \frac{1}{\sqrt{2\pi V_0}} \int  d\tilde{h}_0 \exp\left[-\frac{(\tilde{h}_0 - \langle \tilde{h}_0\rangle _{\setminus0})^2}{2V_0} +\beta \tilde{h}_0 s_0 \right] \\

&= \frac{1}{\sqrt{2\pi V_0}} \int  d\tilde{h}_0 \exp\Biggl[ -\frac{1}{2V_0}\left(\tilde{h}_0^2 -2\tilde{h}_0\langle \tilde{h_0}\rangle_{\setminus0} +\langle \tilde{h_0}\rangle_{\setminus0}^2 -2V_0\beta s_0 \tilde{h}_0 \right) \Biggr]\\

&=\frac{1}{\sqrt{2\pi V_0}} \int  d\tilde{h}_0 \exp\Biggl[ -\frac{1}{2V_0}\Bigl(\tilde{h}_0^2 -2\tilde{h}_0(\langle \tilde{h_0}\rangle_{\setminus0}+V_0\beta s_0)  \\
&\qquad \qquad +(\tilde{h}_0\langle \tilde{h_0}\rangle_{\setminus0}+V_0\beta s_0) ^2-(\tilde{h}_0\langle \tilde{h_0}\rangle_{\setminus0}+V_0\beta s_0) ^2 +\langle \tilde{h_0}\rangle_{\setminus0}^2 \Bigr) \Biggr]\\

&=\frac{1}{\sqrt{2\pi V_0}} \int  d\tilde{h}_0 \exp\Biggl[ -\frac{1}{2V_0}\Bigl( \tilde{h}_0 -(\langle \tilde{h_0}\rangle_{\setminus0}+V_0\beta s_0)\Bigr)^2\Biggr]\\
&\qquad \qquad \qquad × \exp\Biggl[-\frac{1}{2V_0}\Bigl(-( \langle \tilde{h_0}\rangle_{\setminus0}+V_0\beta s_0) ^2 +\langle \tilde{h_0}\rangle_{\setminus0}^2 \Bigr)\Biggr]\\

&=\frac{1}{\sqrt{2\pi V_0}}\exp\Biggl[-\frac{1}{2V_0}\Bigl(-( \langle \tilde{h_0}\rangle_{\setminus0}+V_0\beta s_0) ^2 +\langle \tilde{h_0}\rangle_{\setminus0}^2 \Bigr)\Biggr] \sqrt{2\pi V_0}) \\

&=\exp\Biggl[-\frac{1}{2V_0}\Bigl(-( \langle \tilde{h_0}\rangle_{\setminus0}+V_0\beta s_0) ^2 +\langle \tilde{h_0}\rangle_{\setminus0}^2 \Bigr)\Biggr] \\

&=\exp\Biggl[ \frac{1}{2V_0}\Bigl( 2 \langle \tilde{h_0}\rangle_{\setminus0}V_0\beta s_0  +(V_0\beta s_0)^2  \Bigr)\Biggr] \\

&=\exp\left[\frac{1}{2V_0}(V_0\beta)^2 \right]\exp\Bigl[    \beta \langle \tilde{h_0}\rangle_{\setminus0} s_0\Bigr] 
\end{align}$$
*/

\begin{eqnarray*}
P(s_0) &\propto& \int d \tilde{h}_0 \ e^{\beta \tilde{h}_0 s_0} P(\tilde{h}_0 \setminus s_0) \\
&=& \int d\tilde{h}_0 \frac{1}{\sqrt{2\pi V_0}} e^{\beta \tilde{h}_0 s_0} \exp \left[ -\frac{(\tilde{h}_0 - \langle \tilde{h}_0\rangle _{\setminus0})^2}{2V_0}\right]\\
&=& \exp\left[\frac{1}{2V_0}(V_0\beta)^2 \right]\exp\Bigl[    \beta \langle \tilde{h_0}\rangle_{\setminus0} s_0\Bigr].
\end{eqnarray*}
중간에 $s_0^2=1$을 활용했다.
이제 $P(s_0)$가 $\exp\Bigl[    \beta \langle \tilde{h_0}\rangle_{\setminus0} s_0\Bigr] $에 비례한다는 것을 알았으므로, 아래와 같이 $m_0$룰 계산하자.

$$
\begin{align}

m_0 &= \sum_{s_0=\pm 1} \frac{1}{[\exp[\beta \langle \tilde{h_0}\rangle_{\setminus0} \ + \exp[-\beta \langle \tilde{h_0}\rangle_{\setminus0}]]} s_0\exp\bigl[ \beta \langle \tilde{h_0}\rangle_{\setminus0} s_0\bigr] \\
&=\tanh[\beta \langle \tilde{h_0}\rangle_{\setminus0}]
\end{align}
$$

여기에서 $\langle \tilde{h_0}\rangle_{\setminus0}$라고 함은 스핀 $N$개에 대한 $\tilde{h}_0$의 평균값이다.
그러므로 $\langle \tilde{h_0}\rangle_{\setminus0}$의 식을, 스핀 $N+1$개에 대한 $\langle \tilde{h_0}\rangle$에 대해서 얻으면 $m_0$를 온전히 표현할 수 있다.

$\langle \tilde{h_0}\rangle_{\setminus0}$이 아닌 $\langle \tilde{h_0}\rangle$ 로서 '$s_0$를 제외 시키지 않은' 경우의 $\tilde{h}_0$의 기대값은 다음과 같이 쓸 수 있다.

$$ \langle \tilde{h}_0 \rangle =
\text{Tr}_{s_0} \int d\tilde{h}_0 \tilde{h}_0 P(s_0,\tilde{h}_0) $$

앞서 살펴본 대로 $P(s_0,\tilde{h}_0) \propto e^{\beta\tilde{h}_0 s_0}P(\tilde{h}_0 \setminus s_0)$를 이용하면 
또 한번의 [[수학:윅의_정리|가우스 적분]]에 의해서 $\langle \tilde{h_0}\rangle$와 $\langle \tilde{h_0}\rangle_{\setminus0}$에 대한 관계식을 얻는다.
/*
$$
\begin{align}
\langle \tilde{h}_0 \rangle &\propto \frac{1}{\sqrt{2\pi V_0}} \ \text{Tr}_{s_0} \int d\tilde{h}_0 \tilde{h}_0\exp \Biggl[ -\frac{(\tilde{h}_0 - \langle \tilde{h}_0\rangle _{\setminus0})^2}{2V_0} +\beta\tilde{h}_0 s_0 \Biggr]\\

&=  \frac{1}{\sqrt{2\pi V_0}} \text{Tr}_{s_0}  \int  d\tilde{h}_0 \tilde{h}_0 \exp\Biggl[ -\frac{1}{2V_0}\Bigl( \tilde{h}_0 -(\langle \tilde{h_0}\rangle_{\setminus0}+V_0\beta s_0)\Bigr)^2\Biggr]\\
&\qquad \qquad \qquad × \exp\Biggl[-\frac{1}{2V_0}\Bigl(-( \langle \tilde{h_0}\rangle_{\setminus0}+V_0\beta s_0) ^2 +\langle \tilde{h_0}\rangle_{\setminus0}^2 \Bigr)\Biggr]\\

&=\frac{1}{\sqrt{2\pi V_0}} \exp\left[\frac{1}{2V_0}(V_0\beta)^2 \right] × \\
&\text{Tr}_{s_0} \exp\Bigl[    \beta \langle \tilde{h_0}\rangle_{\setminus0} s_0\Bigr] \int  d\tilde{h}_0 \tilde{h}_0 \exp\Biggl[ -\frac{1}{2V_0}\Bigl( \tilde{h}_0 -(\langle \tilde{h_0}\rangle_{\setminus0}+V_0\beta s_0)\Bigr)^2\Biggr]\\


&=\frac{1}{\sqrt{2\pi V_0}} \exp\left[\frac{1}{2V_0}(V_0\beta)^2 \right] ×\\
&\text{Tr}_{s_0}\exp\Bigl[    \beta \langle \tilde{h_0}\rangle_{\setminus0} s_0\Bigr]  \int  d\tilde{h}_0 \Bigl(  \langle \tilde{h_0}\rangle_{\setminus0}+V_0\beta s_0 \Bigr) \exp\Biggl[ -\frac{1}{2V_0}\Bigl( \tilde{h}_0 -(\langle \tilde{h_0}\rangle_{\setminus0}+V_0\beta s_0)\Bigr)^2\Biggr]\\

&= \exp\left[\frac{1}{2V_0}(V_0\beta)^2 \right]\text{Tr}_{s_0} \exp\Bigl[    \beta \langle \tilde{h_0}\rangle_{\setminus0} s_0\Bigr]  \Bigl(\langle \tilde{h_0} \rangle_{\setminus0}+V_0\beta s_0 \Bigr) \\

&=  \langle \tilde{h_0} \rangle_{\setminus0}+V_0\beta \langle s_0 \rangle

\end{align}
$$
*/
\begin{eqnarray*}
\langle \tilde{h}_0 \rangle &\propto& \frac{1}{\sqrt{2\pi V_0}} \ \text{Tr}_{s_0} \int d\tilde{h}_0 \tilde{h}_0\exp \Biggl[ -\frac{(\tilde{h}_0 - \langle \tilde{h}_0\rangle _{\setminus0})^2}{2V_0} +\beta\tilde{h}_0 s_0 \Biggr]\\
&=& \exp\left[ \frac{1}{2V_0} (\beta V_0)^2 \right] \text{Tr}_{s_0} \left( \langle \tilde{h}_0\rangle _{\setminus0} + \beta s_0 V_0 \right) \exp \left( \beta s_0 \langle \tilde{h}_0\rangle _{\setminus0} \right)\\
&=&  \langle \tilde{h_0} \rangle_{\setminus0}+V_0\beta \langle s_0 \rangle\\
&=& \langle \tilde{h_0} \rangle_{\setminus0}+V_0\beta m_0.
\end{eqnarray*}
/*
$$ 
\begin{align}
\therefore  \langle \tilde{h_0} \rangle &= \langle \tilde{h_0} \rangle_{\setminus0}+V_0\beta \langle s_0 \rangle \\
&= \langle \tilde{h_0} \rangle_{\setminus0}+V_0\beta m_0
\end{align}$$
*/
위의 유도 과정에서는 앞서서 계산한 적분 결과를 곧바로 사용하였고, 
$P(s_0) \propto \exp\left[\frac{1}{2V_0}(V_0\beta)^2 \right]\exp\Bigl[    \beta \langle \tilde{h_0}\rangle_{\setminus0} s_0\Bigr] $임을 적용하여 관계식을 등호로 바꿀 수 있었다.

이때 '$s_0$를 제외 시키지 않은' 경우에 해당하는 국소장의 기댓값은 $\langle \tilde{h_0}\rangle=\sum_j J_{0j} m_j$이므로
앞서 얻었던 식인 $m_0 = \tanh[\beta \langle \tilde{h_0}\rangle_{\setminus0}]$은 결론적으로 아래와 같다.

\begin{align}
m_0 &=  \tanh\left[\beta \left( \langle \tilde{h_0} \rangle - V_0\beta m_0 \right)\right] \\
&=  \tanh\left[\beta \left( \sum_j J_{0j} m_j -\beta \sum_j J_{0j}^2 (1-m_j^2)m_0 \right)\right]
\end{align}

위의 유도 결과는 TAP 방정식이라고 불린다. $s_0$가 포함 되었으므로, 스핀 $N+1$개에 대한 계를 기술한다.

$$ \\ $$

===== TAP 방정식의 가정 =====

TAP 방정식은 고온의 영역에서만 옳은 결과를 준다는 한계가 있다.

우리가 유도 과정을 거칠 때에 사용하였던 주된 가정은 $\sum_{j,k} (\langle s_j s_k\rangle - \langle s_j \rangle \langle s_k \rangle)$이 
$j \ne k$ 일 때는 $N \to \infty$에서 크게 기여하지 못한다는 가정이었다.

그러나 저온의 영역에서는  $N \to \infty$에서도 $j,k$에 대해서 $(\langle s_j s_k\rangle - \langle s_j \rangle \langle s_k \rangle)$가 유한한 값을 가지므로
TAP 방정식이 정확히 옳은 결과를 주지 못한다.

===== 온사거 반응 장 =====

TAP 방정식인 다음의 식을 다시 살펴보면

$$ m_0 =  \tanh\left[\beta \left( \sum_j J_{0j} m_j -\beta \sum_j J_{0j}^2 (1-m_j ^2)m_0 \right)\right]$$

$-\beta \sum_j J_{0j}^2 (1-m_j ^2)m_0   $의 항은 마치 추가로 붙게 된 항처럼 보인다.
사실 이 항은 상쇄의 역할을 하며 중요한 의미를 갖는다.

자화량 $m_0$는 $j$의 위치에 내부장(internal field) $J_{0j}m_0$만큼의 영향을 준다.
그에 따라, 위치 $j$의 자화량 $m_j$는 $\chi_{jj}J_{0j}m_0$만큼 변하게 된다. $\chi_{jj}$는 다음과 같다.
$$ \chi_{jj} = \frac{\partial m_j}{\partial h_j} \Biggl|_{h_j \to 0}=\beta(1-m_j^2)$$
$\chi$가 자기 감수율(magnetic susceptibility)이며, 그 값은 앞서 분산 $V_0$를 구할 때와 마찬가지의 방식으로 구할 수 있다.

이때, 이렇게 변화한 $m_j$의 값에 의해서 $m_0$가 다시 영향을 받게 되는 값은 다음과 같다:
$$ J_{0j} \times \chi_{jj}J_{0j}m_0 = \beta J_{0j}^2(1-m_j^2 )m_0.$$
그런데 이것을 발생시킨 $m_j$의 변화는 $m_0$에 의한 것이므로 이러한 항의 효과는 제거해주어야 한다.

TAP 방정식에서 $-\beta \sum_j J_{0j}^2 (1-m_j ^2)m_0   $의 항이 해당 역할을 적절하게 수행하며
온사거 반응 장(reaction field)으로 불린다.

===== 복제 대칭 해의 유도 =====

====자기화====

우리가 [[물리:복제_대칭_해|복제 대칭 해]]에서 살펴본 결과는 TAP 방정식을 이용해서도 유도가 가능하다.
우선 상호작용인 $J_{0j}$를 강자성에 대한 항과 난수로 표현되는 항으로 분리하여 아래와 같이 나타내자.

$$ J_{0j} = \frac{J_0}{N} + \frac{J}{\sqrt{N}}z_{0j} $$

여기에서 $z_{0j}$는 평균이 $0$이며 분산이 $1$인 정규 분포에서 추출한 난수이다.
앞서 얻었던 $ \langle \tilde{h_0} \rangle = \langle \tilde{h_0} \rangle_{\setminus0}+V_0\beta m_0$의 식에 이를 대입하면 다음과 같다.

$$\langle \tilde{h_0} \rangle_{\setminus0} = \frac{J_0}{N}\sum_j m_j + \frac{J}{\sqrt{N}}\sum_j z_{0j}m_j - V_0 \beta m_0 $$

$m$은 강자성 질서 맺음 변수라고 하면, 첫번째 항은 $J_0 m$과 같다.
세번째 항의 역할은 온사거의 반응 장에서 살펴본 공동 보정(cavity correction)의 효과를 주므로
두번째 항의 $ z_{0j}m_j$가 독립적인 담금질(quenched) 무작위 변수인 것으로 가정할 수 있도록 한다.

보다 자세히 설명하면, $\langle \tilde{h_0} \rangle_{\setminus0}$를 표현하는 식의 의미는 '공동 장'의 평균으로서 $s_0$가 제거된 경우를 보는 것이므로 
한 위치 $j$로 부터 발생하는 기여분은 또 다른 $j$로 부터 발생하는 기여분과 서로 간섭하지 않는다.
따라서 '중심 극한 정리'를 적용할 수 있고, $\sum_j z_{0j}m_j$가 가우스 분포를 따른다고 가정할 수 있다.
즉 평균은 $0$이고, 그에 따라 분산은 다음과 같이 표현된다.

$$ \sum_j \sum_k [z_{0j}z_{0k}]m_j m_k = \sum_j m_j^2 = N q$$

그러므로 두번째 항은 정규 분포에서 추출된 난수 $z$를 이용하여 $\sqrt{Nq}z$로서 표현이 가능하다.
이러한 결과들을 가지고서, $m_0 = \tanh \left(\beta \langle \tilde{h_0}\rangle_{\setminus0} \right)$를 $z$의 분포에 대해서 평균을 내리면

$$ m= \int \text{D}z \tanh (\beta J_0 m + \beta J\sqrt{q}z) $$
의 결과를 얻고, 이는 [[물리:복제_대칭_해|복제 대칭 해]]에 대해서 설명된 식과 같다.

==== 에드워즈-앤더슨 질서 변수 ====
$N+1$개의 스핀이 있는 계에 대해 다음 질서 변수를 계산해보자:
\[q = \overline{\langle s_0 \rangle^2} = \overline{\tanh \left(\beta h \right)^2} \]
여기에서 $h \equiv \langle \tilde{h}_0 \rangle_{\setminus 0} = \sum_{i=1}^N J_{0i} \langle s_i \rangle_{\setminus 0}$는 $s_0$ 없이 $N$개의 스핀만 있을 때의 공동 장의 평균을 간단히 표기한 것이다. 무질서 평균을 취하면, 결합상수들의 평균이 $0$이므로 아래 결과를 얻는다:
\[ \overline{h} = \sum_i \overline{J_{0i} \langle s_i \rangle_{\setminus 0}} = 0.\]
제곱을 해서 무질서 평균을 취하면 아래 이중합에서 $i=j$인 항들만 주로 기여하고 나머지는 자기들끼리 상쇄되므로
\[ \overline{h^2} = \sum_{i,j} \overline{J_{0i} J_{0j} \langle s_i \rangle_{\setminus 0} \langle s_j \rangle_{\setminus 0}} \approx \frac{1}{N} \sum_i \overline{\langle s_i \rangle^2_{\setminus 0}} = q. \]
이렇게 통계적 모멘트들을 따져보면 무질서의 샘플들로부터 보았을 때 $h$가 분산 $q$인 가우스 분포를 따름을 알 수 있다. 따라서
\[ q = \overline{\tanh \left(\beta h \right)^2} = \int \frac{dh}{\sqrt{2\pi q}} \exp\left( -\frac{h^2}{2q} \right) \tanh^2 \left(\beta h\right). \]
==== 비선형 감수율 ====

다음의 양을 계산하기 위해
\[ \chi_\text{nl} = \frac{\beta^2}{N} \sum_{i,j} \overline{\left( \langle s_i s_j \rangle - \langle s_i \rangle \langle s_j \rangle \right)^2} \]
국소적인 감수율인 다음의 양을 고려하도록 하자:
\[ \chi_{ij} = \beta \left( \langle s_i s_j \rangle - \langle s_i \rangle \langle s_j \rangle \right). \]
[[물리:요동-흩어지기 정리|]]에 의하면 이 양은 $\chi_{ij} = dm_i / dh_j = \chi_{ji}$로도 쓸 수 있으며, $h_i$는 $s_i$의 위치에 아주 작게 걸리는 국소적인 외부 장이다.

계산을 위해 두 개의 스핀 $s_0$와 $s_{0'}$을 계에 추가해보자. $N+2$개의 스핀을 가지는 계의 확률분포는
\[ P\left( \vec{s}, s_0, s_{0'}\right) = \exp\left[ -\beta H_N \left(\vec{s}\right) + \beta \tilde{h}_0 s_0 + \beta \tilde{h}_{0'} s_{0'} + \beta J_{00'} s_0 s_{0'} \right] \]
로서, $\tilde{h}_0 \equiv \sum_{i=1}^N J_{0i} s_i$이고 $\tilde{h}_{0'} \equiv \sum_{i=1}^N J_{0'i} s_i$이다.
$\tilde{h}_0$와 $\tilde{h}_{0'}$의 결합확률이 다음처럼 주어진다:
\[ P\left( \tilde{h}_0, \tilde{h}_{0'} \right) \propto \exp\left[ - \frac{(\tilde{h}_0 - h)^2}{2(1-q)} - \frac{(\tilde{h}_{0'} - h')^2}{2(1-q)} + \epsilon (\tilde{h}_0 - h) (\tilde{h}_{0'} - h')\right]. \]
표기를 간단히 하기 위해 원래의 스핀 $N$개만 있을 때에 공동 장을 평균한 값들을 $h \equiv \langle \tilde{h}_0 \rangle_{\setminus \left\{0,0'\right\}}$과 $h' \equiv \langle \tilde{h}_{0'} \rangle_{\setminus \left\{0,0'\right\}}$으로 적었다.

일반적으로 상관함수
\[ \Big\langle \left( \tilde{h}_0 - \langle \tilde{h}_0 \rangle \right) \left( \tilde{h}_{0'} - \langle \tilde{h}_{0'} \rangle \right) \Big\rangle 
= \sum_{i,j} J_{0i} J_{0'j} \left( \langle s_i s_j \rangle - \langle s_i \rangle \langle s_j \rangle \right) \]
가 가지는 크기는 $i=j$인 $N$개의 항들에 의해 주로 결정되고, 각각의 항의 크기는 $1/N$ 정도인데 부호가 제각각이므로 결국 이 상관함수의 크기는 약 $1/\sqrt{N}$ 정도이다 (크기와 부호는 샘플마다 다르다). 따라서 결합확률 안에서 $\tilde{h}_0$와 $\tilde{h}_{0'}$의 상관도를 나타내는 $\epsilon$은 $1/\sqrt{N}$ 정도의 크기를 가지는 값이다.

추가된 두 개의 스핀까지 고려한 확률분포 함수는 아래와 같다:
\[ P\left( \tilde{h}_0, \tilde{h}_{0'}, s_0, s_{0'} \right) \propto P\left( \tilde{h}_0, \tilde{h}_{0'} \right) \exp\left[ \beta \tilde{h}_0 s_0 + \beta \tilde{h}_{0'} s_{0'} + \beta J_{00'} s_0 s_{0'} \right]. \]
[[수학:윅의_정리|가우스 함수의 적분]]을 통해 $\tilde{h}_0$와 $\tilde{h}_{0'}$을 주변화(marginalization)한다:
\begin{eqnarray*}
P\left( s_0, s_{0'} \right) &\propto& \exp\left\{ \frac{\beta^2(1-q) \left[ 1-\epsilon(1-q)s_0 s_{0'} \right]}{1-\epsilon^2(1-q)^2} + \beta h s_0 + \beta h' s_{0'} + \beta J_{00'} s_0 s_{0'} \right\}
\end{eqnarray*}
이 확률분포를 통해 계산했을 때 다음의 결과를 얻는다:
\begin{eqnarray*}
\chi_{00} &=& \beta \left( 1 - \langle s_0 \rangle^2 \right) = \beta \text{ sech}^2 \beta h + O\left(1/\sqrt{N}\right)\\
\chi_{00'} &=& \beta \left( \langle s_0 s_{0'} \rangle - \langle s_0 \rangle \langle s_{0'} \rangle \right) = \underbrace{\beta^2 J_{00'} \text{ sech}^2 \beta h \text{ sech}^2 \beta h' + \beta^3 \epsilon (1-q)^2 \text{ sech}^2 \beta h \text{ sech}^2 \beta h'}_{O\left(1/\sqrt{N}\right)} + O\left( 1/N \right).
\end{eqnarray*}
무질서 평균까지 취했을 때 모든 스핀이 대칭적인 상황에 있으므로 $s_0$가 $s_i$ 역할을 하게끔 놓으면
\begin{eqnarray*}
\chi_\text{nl} &=& \frac{1}{N} \sum_{i,j} \overline{ \chi_{ij}^2 }\\
&\approx& \sum_{i=0}^N \overline{ \chi_{i0}^2 }\\
&=& \overline{\chi_{00}^2} + \sum_{i=1}^N \overline{ \chi_{i0}^2 }.
\end{eqnarray*}

$m_0$가 추가되었을 때 $m_i$의 변화는 다음처럼 기술되고
\begin{eqnarray*}
\delta m_i &=& \sum_j \left( \frac{\partial m_i}{\partial h_j} \right) \delta h_j\\
&\approx& \sum_j \chi_{ij} J_{j0} m_0,
\end{eqnarray*}
따라서
\[ \chi_{i0} \approx \frac{dm_i}{dh_0} = \sum_j \chi_{ij} J_{j0} \chi_{00}. \]
이제 $\chi_\text{nl}$을 적어보자:
\begin{eqnarray*}
\chi_\text{nl} &=& \overline{\chi_{00}^2} + \sum_{i=1}^N \overline{\chi_{i0}^2}\\
&\approx& \overline{\chi_{00}^2} + \sum_{i=1}^N \overline{ \left( \sum_j \chi_{ij} J_{j0} \chi_{00} \right)^2}\\
&=& \overline{\chi_{00}^2} + \sum_{i=1}^N \overline{ \chi_{00}^2 \sum_{j,k} \chi_{ij} \chi_{ik} J_{j0} J_{k0} }\\
&\approx& \overline{\chi_{00}^2} + \sum_{i=1}^N \overline{ \chi_{00}^2 \sum_{j} \chi_{ij}^2 \frac{1}{N}}\\
&=& \overline{\chi_{00}^2} + \overline{ \chi_{00}^2 \sum_{j} \chi_{ij}^2}\\
&\approx& \overline{\chi_{00}^2} \left( 1 + \sum_{j} \overline{ \chi_{ij}^2} \right)\\
&\approx& \overline{\chi_{00}^2} \left( 1 + \chi_\text{nl} \right).
\end{eqnarray*}
따라서
\[ \chi_\text{nl} \approx \frac{\overline{\chi_{00}^2}}{1-\overline{\chi_{00}^2}} \approx \frac{\beta^2 \overline{\text{sech}^4 \beta h}}{1-\beta^2 \overline{\text{sech}^4 \beta h}}. \]
즉 [[신경망:공동_방법|공동 방법]]에서 $\chi_{\text{nl}}$이 유한하다는 가정은 $\beta^2 \overline{\text{sech}^4 \beta h}<1$을 의미하고, 이는 [[물리:드알메이다-사울레스_선|드알메이다-사울레스 선]]으로 나뉘어진 고온 영역을 가리킨다.

======구면 $p$-스핀 유리 모형======
스핀 $N$개를 가진 [[물리:구면_p-스핀_유리_모형|구면 $p$-스핀 유리 모형]]의 해밀토니안을 다음처럼 적자:
$$H = - \sum_{i_1<\ldots<i_p} J_{i_1 \ldots i_p} s_{i_1} \cdots s_{i_p} - \frac{r}{2} \sum_{i=1}^N s_i^2.$$
이때 $r$은 구면 조건을 만족하기 위해 조정해야 하는 맺음변수이다.
결합 상수들은 평균이 $0$이고 분산이 $p!/(2N^{p-1})$인 정규 분포로부터 추출된 담금질(quenched) 난수들이다.

=====추가된 스핀=====
이제 $(N+1)$번째 스핀 $s_0$를 추가하자. 이 스핀이 가지고 들어오는 결합 상수는 $J_{0i_2 \ldots i_p}$이다. 스핀 $s_0$에 작용하는 국소장은 다음과 같다:
$$\tilde{h}_0 = \sum_{i_2<\ldots<i_p} J_{0 i_2 \ldots i_p} s_{i_2} \cdots s_{i_p}.$$

스핀 $s_0$가 없을 때의 국소장인 $\tilde{h}_0\backslash s_0$의 분포를 알아보자. 스핀 $s_0$가 없으므로 독립적인 난수의 합이어서 정규 분포로 근사할 수 있다:
$$P \left(\tilde{h}_0\backslash s_0\right) = \frac{1}{\sqrt{2\pi V_0}} \exp\left[ - \frac{\left(\tilde{h}_0 - \langle \tilde{h}_0 \rangle_{\backslash 0}\right)^2}{2V_0} \right].$$

결합 상수와 스핀 값들이 모두 독립적이므로 이 분포의 평균은 아래와 같다:
$$\langle \tilde{h}_0 \rangle_{\backslash 0} = \sum_{i_2<\ldots<i_p} J_{0 i_2 \ldots i_p} \langle s_{i_2}\rangle_{\backslash 0} \cdots \langle s_{i_p} \rangle_{\backslash 0}.$$
[[수학:인자_그래프|인자 그래프]]에서 계산한 방법대로 $\langle s_i \rangle_{\backslash 0}$을 $m_i \equiv \langle s_i \rangle$에 온사거 보정항(Onsager correction term)을 더함으로써 표현하면,
\begin{eqnarray*}
\langle \tilde{h}_0 \rangle_{\backslash 0} &\approx& \sum_{i_2<\ldots<i_p} J_{0 i_2 \ldots i_p} m_{i_2} \cdots m_{i_p} - \beta (1-q) m_0 \frac{p(p-1)}{2} q^{p-2}.
\end{eqnarray*}

분산은 다음처럼 계산되는데 이 중에서 대각항($i_1=k_1, \ldots, i_p=k_p$)만 살아남고 나머지 교차항들은 상쇄되어 평균적으로 0이라고 근사하면
\begin{eqnarray*}
V_0 &=& \sum_{i_2<\ldots<i_p} J_{0 i_2 \ldots i_p} \left( s_{i_2} \cdots s_{i_p} - \langle s_{i_2}\rangle_{\backslash 0} \cdots \langle s_{i_p} \rangle_{\backslash 0} \right) \sum_{k_2<\ldots<k_p} J_{0 k_2 \ldots k_p} \left( s_{k_2} \cdots s_{k_p} - \langle s_{k_2}\rangle_{\backslash 0} \cdots \langle s_{k_p} \rangle_{\backslash 0} \right)\\
&=& \frac{1}{\left[(p-1)!\right]^2} \sum_{i_2,\ldots,i_p} J_{0 i_2 \ldots i_p} \left( s_{i_2} \cdots s_{i_p} - \langle s_{i_2}\rangle_{\backslash 0} \cdots \langle s_{i_p} \rangle_{\backslash 0} \right) \sum_{k_2,\ldots,k_p} J_{0 k_2 \ldots k_p} \left( s_{k_2} \cdots s_{k_p} - \langle s_{k_2}\rangle_{\backslash 0} \cdots \langle s_{k_p} \rangle_{\backslash 0} \right)\\
&\approx& \frac{1}{(p-1)!} \sum_{i_2,\ldots,i_p} J_{0 i_2 \ldots i_p}^2 \left(1 - \langle s_{i_2}\rangle_{\backslash 0}^2 \cdots \langle s_{i_p} \rangle_{\backslash 0}^2 \right)\\
&\approx& \frac{1}{(p-1)!} \left( \frac{p!}{2N^{p-1}} \right) \left(1 - q^{p-1} \right) N^{p-1}\\
&=& \frac{p}{2} \left(1 - q^{p-1} \right).
\end{eqnarray*}

이제 $s_0$를 포함한 결합 분포 함수를 구하고
\begin{eqnarray*}
P\left( s_0, \tilde{h}_0 \right) &\propto& \exp\left[ -\frac{\beta r}{2}s_0^2 + \beta \tilde{h}_0 s_0 - \frac{\left(\tilde{h}_0 - \langle \tilde{h}_0 \rangle_{\backslash 0}\right)^2}{2V_0} \right]\\
\end{eqnarray*}
$h_0$에 대해 주변화를 시행한 다음,
\begin{eqnarray*}
P\left( s_0 \right) &\propto& \int d\tilde{h}_0 \exp\left[ -\frac{\beta r}{2}s_0^2 + \beta \tilde{h}_0 s_0 - \frac{\left(\tilde{h}_0 - \langle \tilde{h}_0 \rangle_{\backslash 0}\right)^2}{2V_0} \right]\\
&\propto& \exp\left[ - \frac{\beta r}{2} s_0^2 + \beta \langle \tilde{h}_0 \rangle_{\backslash 0} s_0 + \frac12 \beta^2 s_0^2 V_0 \right],
\end{eqnarray*}
이로부터 $s_0$의 평균값과 (단, $r>\beta V_0$)
\begin{eqnarray*}
\langle s_0 \rangle &=& \frac{\int ds_0 s_0 P\left(s_0\right)}{\int ds_0 P\left(s_0\right)} = \frac{\langle \tilde{h}_0 \rangle_{\backslash 0}}{r - \beta V_0}.
\end{eqnarray*}
분산을 계산한다:
\begin{eqnarray*}
\langle \left( s_0 - \langle s_0 \rangle \right)^2 \rangle &=& \frac{\int ds_0 \left( s_0- \langle s_0 \rangle \right)^2 P\left(s_0\right)}{\int ds_0 P\left(s_0\right)} = \frac{1}{\beta \left(r - \beta V_0\right)} = 1-q.
\end{eqnarray*}
마지막의 표현식은 구면 조건과 $q \equiv N^{-1} \sum_i m_i^2$의 정의로부터 나온다. 이것을 이용하면 $s_0$가 실제로 추가되었을 때 그 평균값인 $m_0 \equiv \langle s_0 \rangle$에 대한 또다른 표현식은,
\begin{eqnarray*}
\frac{m_0}{\beta(1-q)} &=& \langle \tilde{h}_0 \rangle_{\backslash 0} \approx \sum_{i_2<\ldots<i_p} J_{0 i_2 \ldots i_p} m_{i_2} \cdots m_{i_p} - \beta (1-q) m_0 \frac{p(p-1)}{2} q^{p-2}.
\end{eqnarray*}

$\partial q/\partial m_i = 2m_i/N$에 유의하면, 위의 식은 아래와 같은 자유 에너지 밀도를 미분함으로써 얻을 수 있다:
$$f = -\frac{1}{N} \sum_{i_1<\ldots<i_p} J_{i_1 \ldots i_p} m_{i_1} \cdots m_{i_p} - \frac{1}{2\beta} \ln(1-q) - \frac{\beta}{4} \left[ (p-1) q^p - pq^{p-1} +1 \right].$$

/*
=====다른 스핀과의 관계=====
스핀 $0$가 아닌 다른 스핀 $i$를 포함하여 식을 적어보자.
그러기 위해 국소장을 다음처럼 $i$를 포함한 항과 포함하지 않은 항으로 나누어 쓰자:
\begin{eqnarray*}
\tilde{h}_0 &=& \sum_{i_2<\ldots<i_p} J_{0 i_2 \ldots i_p} s_{i_2} \cdots s_{i_p}\\
&=& s_i \sum_{i_3<\ldots<i_p} J_{0 i i_3 \ldots i_p} \cdots s_{i_p} + \sum_{i_2<\ldots<i_p\\i \notin \{i_2,\ldots,i_p\}} J_{0 i_2 \ldots i_p} s_{i_2} \cdots s_{i_p}\\
&=& s_i \tilde{h}_i + \tilde{h}_i'.
\end{eqnarray*}
스핀 $s_0$가 없을 때 $\tilde{h}_i$의 평균과 분산을 각각 $\langle \tilde{h}_i \rangle_{\backslash 0}$과 $\Bigl< \left( \tilde{h}_i - \langle \tilde{h}_i \rangle_{\backslash 0} \right)^2 \Bigr>_{\backslash 0}$라고 하자. [[수학:인자_그래프|인자 그래프]]에서 다룬 바와 같이
\begin{eqnarray*}
\Bigl< \left( \tilde{h}_i - \langle \tilde{h}_i \rangle_{\backslash 0} \right)^2 \Bigr>_{\backslash 0} &=& \langle \tilde{h}_i^2 \rangle_{\backslash 0} - \langle \tilde{h}_i \rangle_{\backslash 0}^2\\
&\approx& \frac{1}{\left[(p-2)!\right]^2} (p-2)! \sum_{i_3,\ldots,i_p} J_{0 i i_3 \ldots i_p}^2 \left( \langle s_{i_3}^2 \rangle_{\backslash 0} \cdots \langle s_{i_p}^2 \rangle_{\backslash 0} - \langle s_{i_3} \rangle_{\backslash 0}^2 \cdots \langle s_{i_p}\rangle_{\backslash 0}^2 \right)\\
&\approx& \frac{1}{(p-2)!} \sum_{i_3,\ldots,i_p} \frac{p!}{2N^{p-1}} \left( 1 - q^{p-2} \right)\\
&=& \frac{1}{(p-2)!} N^{p-2} \frac{p!}{2N^{p-1}} \left( 1 - q^{p-2} \right)\\
&=& \frac{p(p-1)}{2N} \left( 1 - q^{p-2} \right).
\end{eqnarray*}
마찬가지로 스핀 $s_0$가 없을 때 $\tilde{h}_i'$의 평균과 분산을 각각 $\langle \tilde{h}_i' \rangle_{\backslash 0}$과 $\Bigl< \left( \tilde{h}_i' - \langle \tilde{h}_i' \rangle_{\backslash 0} \right)^2 \Bigr>_{\backslash 0}$라고 하자. 셈에서 스핀 $i$가 빠지기는 하지만 $N$이 매우 크다면 이 계산은 사실상 위에서 했던 $V_0$의 계산과 같아진다. 따라서
$$\Bigl< \left( \tilde{h}_i' - \langle \tilde{h}_i' \rangle_{\backslash 0} \right)^2 \Bigr>_{\backslash 0} = \frac{p}{2} \left( 1- q^{p-1} \right).$$
끝으로, 스핀 $s_0$가 없을 때에 $s_i$는 평균 $\langle s_i \rangle_{\backslash 0}$와 분산 $1-q$를 가진다.

$s_0$를 계에 추가하고 $s_0$, $s_i$, $\tilde{h}_i$, 그리고 $\tilde{h}_i'$의 결합 확률을 구하자. 앞서와 마찬가지로 $s_0$가 없을 때의 확률분포를 먼저 적고 나서 $s_0$로 인한 볼츠만 인자 $\exp\left( -\beta r s_0^2/2 + \beta h s_0\right)$를 곱하는 식으로 할 것이다. $s_0$가 없을 때, $\tilde{h}_0$를 구성하는 $s_i$, $\tilde{h}_i$, 그리고 $\tilde{h}_i'$들은 서로 독립적인 난수라고 가정한다.
\begin{eqnarray*}
P\left(s_0, s_i, \tilde{h}_i, \tilde{h}_i' \right) &\propto& \exp\left[ -\frac{\beta r}{2} s_0^2 + \beta s_0 \left(s_i \tilde{h}_i + \tilde{h}_i' \right) - \frac{\left( s_i - \langle s_i \rangle_{\backslash 0} \right)^2}{2(1-q)} - N \frac{\left( \tilde{h}_i - \langle \tilde{h}_i \rangle_{\backslash 0} \right)^2}{p(p-1)\left( 1-q^{p-1} \right)} - \frac{\left( \tilde{h}_i' - \langle \tilde{h}_i' \rangle_{\backslash 0} \right)^2}{p\left( 1-q^{p-1} \right)} \right].
\end{eqnarray*}
주변화를 거치면 아래와 같은 적분 결과를 얻고
\begin{eqnarray*}
P\left(s_0, s_i \right) &=& \int d\tilde{h}_i d\tilde{h}_i' P\left(s_0, s_i, \tilde{h}_i, \tilde{h}_i' \right)\\
&\propto& \exp\left[ -\frac{\beta r}{2}s_0^2+\frac{p\beta^2}{4}\left(1-q^{p-1}\right) s_0^2 + \beta s_0 \left( s_i \langle \tilde{h}_i \rangle_{\backslash 0} + \langle \tilde{h}_i' \rangle_{\backslash 0} \right)-\frac{\left(s_i-\langle s_i \rangle_{\backslash 0}\right)^2}{2(1-q)} -\frac{1}{4N}\beta^2 (1-p)p \left(1-q^{p-1}\right) s_0^2 s_i^2 \right]\\
&\approx& \exp\left[ -\frac{s_0^2}{2(1-q)} + \beta s_0 \left( s_i \langle \tilde{h}_i \rangle_{\backslash 0} + \langle \tilde{h}_i' \rangle_{\backslash 0} \right)-\frac{\left(s_i-\langle s_i \rangle_{\backslash 0}\right)^2}{2(1-q)} \right],
\end{eqnarray*}
지수 중 $N^{-1}$에 비례하는 항은 $N\to\infty$에서 무시했다. 이제 주변화로 $P(s_i) = \int ds_0 P(s_0, s_i)$를 얻은 다음 $\langle s_i \rangle$을 $\tilde{h}_i \sim O\left(N^{-1/2}\right)$의 일차항까지 계산해보면,
\begin{eqnarray*}
\langle s_i \rangle &=& \frac{\int ds_i ~s_i P\left( s_i \right)}{\int ds_i P\left( s_i \right)} = \frac{\langle s_i \rangle_{\backslash 0} + \beta^2 (1-q)^2 \langle \tilde{h}_i \rangle_{\backslash 0} \langle \tilde{h}_i' \rangle_{\backslash 0}}{1-\beta^2(1-q)^2 \langle \tilde{h}_i \rangle_{\backslash 0}^2} \approx \langle s_i \rangle_{\backslash 0} + \beta^2 (1-q)^2 \langle \tilde{h}_i \rangle_{\backslash 0} \langle \tilde{h}_i' \rangle_{\backslash 0}\\
&=& \langle s_i \rangle_{\backslash 0} + \beta (1-q) \langle s_0 \rangle \langle \tilde{h}_i' \rangle_{\backslash 0}\\
&\approx& \langle s_i \rangle_{\backslash 0} + \beta (1-q) \langle s_0 \rangle \sum_{i_3<\ldots<i_p} J_{0 i i_3 \ldots i_p} \langle s_{i_3} \rangle_{\backslash 0} \cdots \langle s_{i_p} \rangle_{\backslash 0}.
\end{eqnarray*}
*/
======함께 보기======
  * [[신경망:공동_방법|공동 방법]]
  * [[수학:인자_그래프|인자 그래프]]
  * [[물리:드알메이다-사울레스_선|드알메이다-사울레스 선]]
  * [[물리:구면_p-스핀_유리_모형|구면 $p$-스핀 유리 모형]]

====== 참고문헌 ======
  * Hidetoshi Nishimori, //Statistical Physics of Spin Glasses and Information Processing: An Introduction//, 2001.
  * Gino Del Ferraro, Chuang Wang, Dani Martí, and Marc Mézard, //Cavity Method: Message Passing from a Physics Perspective// [[https://doi.org/10.48550/arXiv.1409.3048|arXiv:1409.3049]] (2014).
  * Marc Mézard, Giorgio Parisi, and Miguel Ángel Virasoro, //Spin Glass Theory and Beyond//, 1987.
  * A. Barrat, //The p-spin spherical spin glass model//, [[https://arxiv.org/pdf/cond-mat/9701031|cond-mat/9701031]] (1997).
  * T. Rizzo, //TAP complexity, the cavity method and supersymmetry//, [[https://doi.org/10.1088/0305-4470/38/15/005|J. Phys. A: Math. Gen. 38 3287 (2005)]].