=====$n$개의 복제본에 대한 분배함수===== 아래 식의 좌변에서 $Z_n$을 감싸고 있는 꺽쇠는 무질서평균을 뜻하며, $\alpha, \beta, \gamma, \delta, \ldots$의 그리스 문자들은 복제본을 가리키는 인덱스들이다. \begin{align} \left[ Z^n \right] &= \int \left[ \prod_{i \beta} S_{i}^{\alpha} S_{j}^{\alpha} S_{i}^{\beta} S_{j}^{\beta} \right) \\ &= \sum_{ij} \left[ 2 \sum_{\alpha < \beta} S_{i}^{\alpha} S_{j}^{\alpha} S_{i}^{\beta} S_{j}^{\beta} + n \right] - \sum_{i=j} \left[ 2 \sum_{\alpha < \beta} S_{i}^{\alpha} S_{j}^{\alpha} S_{i}^{\beta} S_{j}^{\beta} + n \right]. \end{align} 여기에서 두 번째 항은 그 바로 윗줄의 $Y$ 표현식과 사실상 동일하다. 이를 좌변으로 넘겨 $2Y$라고 쓰도록 하자. 그리고 마지막 항은 $i=j$인 경우만 다루므로 $S_i^\alpha S_j^\alpha S_i^\beta S_j^\beta = 1$이다. 이를 $\alpha<\beta$에 대해 더하므로 답은 $n(n-1)/2$가 나온다. 이를 대입하자. \begin{equation} \therefore 2Y = \sum_{i,j} \left[ 2 \sum_{\alpha < \beta} S_{i}^{\alpha} S_{j}^{\alpha} S_{i}^{\beta} S_{j}^{\beta} + n \right] - \sum_i \left( n^2 - n + n \right) \end{equation} 이때 $\sum_{ij}$을 가진 앞의 항은 $N^2$에 비례하는 기여를 하고 $\sum_i$만을 가진 뒤의 항은 $N$에 비례하는 기여를 하는데, $N \gg 1$인 경우를 다루고 있으므로 뒤의 항은 무시할 수 있다. 그러므로 \begin{equation} Y \cong \sum_{i,j} \sum_{\alpha < \beta} S_{i}^{\alpha} S_{j}^{\alpha} S_{i}^{\beta} S_{j}^{\beta} + \frac{N^2}{2}n = \sum_{\alpha < \beta} \left( \sum_i S_{i}^{\alpha} S_{i}^{\beta} \right)^2 + \frac{N^2}{2}n \end{equation} 이다. ===스핀 2개의 결합항=== $W$에서도 인덱스 $i$와 $j$ 사이의 대소관계를 살펴보면 다음처럼 고쳐쓰게 된다: \begin{align} W &= \sum_{ij} \sum_{\alpha} S_{i}^{\alpha} S_{j}^{\alpha} - \sum_{i=j} \sum_{\alpha} S_{i}^{\alpha} S_{j}^{\alpha}. \end{align} 이 때 $\sum_{i>j} \sum_{\alpha} S_{i}^{\alpha} S_{j}^{\alpha}$ 는 $i$와 $j$ 인덱스만 바꾸면 $W$와 같으므로 \begin{equation} 2W = \sum_{i,j} \sum_{\alpha} S_{i}^{\alpha} S_{j}^{\alpha} - \sum_i \sum_\alpha \left( S_{i}^{\alpha} \right)^2 \end{equation} 이며 $\left( S_i^\alpha \right)^2 = 1$임은 자명하다. 여기에서도 앞의 항은 $N^2$에 비례하고 뒤의 항은 $N$에 비례하므로 뒤의 항은 무시한다. \begin{equation} W = \frac{1}{2}\sum_{\alpha} \left( \sum_i S_{i}^{\alpha} \right)^2 - \frac{N}{2}n \cong \frac{1}{2} \sum_{\alpha} \left( \sum_i S_{i}^{\alpha} \right)^2 \end{equation} 따라서 $N \gg 1$ 이면, \begin{align} \left[ Z^n \right] &\cong \text{Tr}_n \exp \left\{ \frac{\beta^2 J^2}{2N} \left[ \sum_{\alpha < \beta} \left( \sum_i S_{i}^{\alpha} S_{i}^{\beta} \right)^2 + \frac{N^2}{2}n \right] + \frac{\beta J_0}{2N} \sum_{\alpha} \left( \sum_i S_{i}^{\alpha} \right)^2 + \beta h \sum_{i} \sum_{\alpha} S_{i}^{\alpha} \right\} \\ &= \exp\left( \frac{\beta^2 J^2}{4}Nn \right) \text{Tr}_n \exp \left\{ \frac{\beta^2 J^2}{2N} \sum_{\alpha < \beta} \left( \sum_i S_{i}^{\alpha} S_{i}^{\beta} \right)^2 + \frac{\beta J_0}{2N}\sum_{\alpha}\left(\sum_i S^{\alpha}_i\right)^2+ \beta h \sum_{i} \sum_{\alpha} S_{i}^{\alpha} \right\} \end{align} 이다. ====변환==== 여기서 제곱항을 처리해주기 위해 [[수학:허바드-스트라토노비치 변환]]을 취하자 (해당 페이지에서 $a = \beta^2 J^2$과 $x = \frac{1}{\sqrt{N}}\sum_i S_i^{\alpha} S_i^{\beta}$를 대입해준 것과 같다). 이 과정에서 보조장(auxiliary field)인 $q_{\alpha \beta}$가 도입된다: \begin{equation} \exp \left[ \frac{\beta^2 J^2}{2N} \left( \sum_i S_{i}^{\alpha} S_{i}^{\beta} \right)^2 \right] = \sqrt{\frac{\beta^2 J^2 N}{2\pi}} \int_{-\infty}^{\infty} dq_{\alpha \beta} \ \exp \left[ -\frac{N \beta^2 J^2}{2} q_{\alpha \beta}^2 + \beta^2 J^2 q_{\alpha \beta} \left( \sum_i S_{i}^{\alpha} S_{i}^{\alpha} \right) \right]. \end{equation} 또다른 제곱항은 보조장 $m_{\alpha}$를 도입하여 마찬가지로 [[수학:허바드-스트라토노비치 변환]]을 취한다: \begin{equation} \exp \left[ \frac{\beta J_0}{2N} \left(\sum_i S_i^{\alpha} \right)^2 \right] = \sqrt{\frac{\beta J_0 N}{2 \pi}} \int_{-\infty}^{\infty} dm_{\alpha} \ \exp \left[ -\frac{N \beta J_0}{2}m_{\alpha}^2 + \beta J_0 m_{\alpha} \left( \sum_i S_{i}^{\beta} \right) \right]. \end{equation} 이 표현식들을 대입하고, 변환으로 인해 생기는 계수들을 생략한다면 다음처럼 적을 수 있다: \begin{equation} \left[ Z^n \right] ~\propto~ \exp \left( \frac{\beta^2 J^2}{4}Nn \right) \int \prod_{\alpha < \beta} dq_{\alpha \beta} \int \prod_{\alpha} dm_{\alpha} \exp \left( -\frac{N \beta^2 J^2}{2} \sum_{\alpha < \beta} q_{\alpha \beta}^2 - \frac{N \beta J_0}{2}\sum_{\alpha} m_{\alpha}^2 \right) \text{Tr}_n \exp \left( \beta^2 J^2 \sum_{\alpha < \beta} q_{\alpha \beta} \sum_i S_{i}^{\alpha} S_{i}^{\beta} + \beta \sum_{\alpha} \left( J_0 m_{\alpha} + h \right) \sum_i S_{i}^{\alpha} \right). \end{equation} 끄트머리의 지수 함수 안에 스핀 위치에 대한 합 $\sum_i$만이 들어있음에 주의하자. 일반적으로 \begin{equation} \text{Tr}_n \exp\left( \sum_i^N L_i \right) = \left( \text{Tr}^\prime e^L \right)^N \end{equation} 과 같은 식으로 분해해서 쓸 수 있다. 예를 들어 $N=2$이고 $n=1$이라면 \begin{align} \text{Tr}_n \exp \left( L_1 + L_2 \right) &= e^{L_1^{+}}e^{L_2^{+}} + e^{L_1^{+}}e^{L_2^{-}} + e^{L_1^{-}}e^{L_2^{+}} + e^{L_1^{-}}e^{L_2^{-}} \\ &= \left( e^{L_1^{+}} + e^{L_1^{-}} \right) \left( e^{L_2^{+}} + e^{L_2^{-}} \right) \end{align} 이라는 뜻이다. 혹은 $N=2$이고 $n=2$라면 \begin{align} \text{Tr}_n \exp \left( L_1 + L_2 \right) &= \left( e^{L_1^{++}} + e^{L_1^{+-}} + e^{L_1^{-+}} + e^{L_1^{--}} \right) \left( e^{L_2^{++}} + e^{L_2^{+-}} + e^{L_2^{-+}} + e^{L_2^{--}} \right) \end{align} 처럼 될 것이다. 물론 $\left( \text{Tr}^\prime e^L \right)^N = \exp \left( N\ln \text{Tr}^{\prime} e^{L} \right)$로도 쓸 수 있다. 여기서 $\text{Tr}^{\prime}$은 특정 위치의 스핀 하나에 대해 취하는 대각합이다 (단, 위치 $i$만을 고정했을 뿐, 다른 복제본들이 개입되었을 수는 있음). 지금의 경우 \begin{align} L \equiv \beta^2 J^2 \sum_{\alpha < \beta} q_{\alpha \beta} S^{\alpha} S^{\beta} + \beta \sum_{\alpha} \left( J_0 m_{\alpha} + h \right) S^{\alpha} \end{align} 가 된다. 다시 정리하면, 분배함수의 무질서평균은 아래와 같다: \begin{equation} \left[ Z^n \right] ~=~ \exp \left( \frac{N \beta^2 J^2 n}{4} \right) \int \prod_{\alpha < \beta} dq_{\alpha \beta} \int \prod_{\alpha} dm_{\alpha} \exp \left( -\frac{N \beta^2 J^2}{2} \sum_{\alpha < \beta} q_{\alpha \beta}^2 - \frac{N \beta J_0}{2}\sum_{\alpha} m_{\alpha}^2 \right) \exp \left( N\ln \text{Tr}^{\prime} e^{L} \right). \end{equation} [[수학:안장점 근사]]를 사용하면 $N \to \infty$에서는 적분되고 있는 지수함수 안의 표현식 \begin{equation} -\frac{\beta^2 J^2}{2} \sum_{\alpha < \beta} q_{\alpha \beta}^2 - \frac{\beta J_0}{2}\sum_{\alpha} m_{\alpha}^2 + \ln \text{Tr}^{\prime} e^{L} \end{equation} 가 최대가 되는 지점이 적분의 값을 결정할 것이라고 논할 수 있다. 이는 위 표현식을 $q_{\alpha \beta}$와 $m_\alpha$로 편미분한 값이 모두 0이 되게끔 하면 찾을 수 있으며, 그 결과는 아래와 같다. \begin{align} &q_{\alpha \beta} = \frac{1}{\beta^2 J^2} \frac{\partial}{\partial q_{\alpha \beta}} \ln \text{Tr}^{\prime} e^{L} = \frac{\text{Tr}^{\prime} S^{\alpha} S^{\beta} e^{L}}{\text{Tr}^{\prime} e^{L}} \\ &m_{\alpha} = \frac{1}{\beta J_0} \frac{\partial}{\partial m_{\alpha}} \ln \text{Tr}^{\prime} e^{L} = \frac{\text{Tr}^{\prime} S^{\alpha} e^{L}}{\text{Tr}^{\prime} e^{L}} \end{align} [[수학:안장점 근사]]를 사용해 적분을 수행한 결과는 위에서 찾아진 $q_{\alpha \beta}$와 $m_\alpha$를 피적분함수에 대입해 적분 밖으로 끄집어낸 결과와 비례한다. 적분으로 인해 얻어지는 비례상수들을 모두 모아서 $C^n$으로 적어놓으면 아래와 같이 쓸 수 있다: \begin{align} \left[ Z^n \right] \propto C^n \exp \left[ -\frac{N \beta^2 J^2}{2} \sum_{\alpha < \beta} q_{\alpha \beta}^2 - \frac{N \beta J_0}{2} \sum_{\alpha} m_{\alpha}^2 + N\ln \text{Tr}^{\prime} e^{L} + \frac{N}{4}\beta^2 J^2 n \right]. \end{align} ====자유에너지==== 이제 입자 하나당의 (무질서평균이 이루어진) 자유에너지 $\left[ f \right]$는 다음의 식을 만족한다. \begin{align} -\beta [f] =& \lim_{n \rightarrow 0} \frac{\left[ Z^n \right] - 1}{nN} = \lim_{n \rightarrow 0} \frac{1}{nN} \left[ -\frac{N \beta^2 J^2}{4} \sum_{\alpha \neq \beta} q_{\alpha \beta}^2 - \frac{N \beta J_0}{2} \sum_{\alpha} m_{\alpha}^2 + \frac{Nn}{4}\beta^2 J^2 + N\ln \text{Tr}^{\prime} e^{L} \right] \\ =& \lim_{n \rightarrow 0} \left[ -\frac{\beta^2 J^2}{4n} \sum_{\alpha \neq \beta} q_{\alpha \beta}^2 - \frac{\beta J_0}{2n} \sum_{\alpha} m_{\alpha}^2 + \frac{\beta^2 J^2}{4} + \frac{1}{n}\ln \text{Tr}^{\prime} e^{L} \right] \end{align} ====보조장의 의미==== 한 위치에서 복제본 사이의 상관관계를 보기 위해 $S_{i}^{\alpha} S_{i}^{\beta}$의 기대값을 계산하자. 먼저 열적평균(thermal average, $\left< \ldots \right>$)을 취한 다음 무질서평균을 취한다. \begin{align} \left[ \big\langle S_{i}^{\alpha} S_{i}^{\beta} \big\rangle \right] =& \left[ \frac{\text{Tr}_n \ S_{i}^{\alpha} S_{i}^{\beta} \exp \left( -\beta \sum_{\gamma} H^{\gamma} \right)}{\text{Tr}_n \ \exp \left( -\beta\sum_{\gamma} H^{\gamma} \right) } \right] \end{align} 여기에서 $H^\gamma = -\sum_{i \right] = \left[ \left< S_i^\alpha \right> \left< S_i^\beta \right> \right] = \left[ \left< S_i \right>^2 \right] = q \end{equation} 와 같은 형태가 된다. 스핀유리상(spin glass phase)에서 스핀들은 위아래를 가리킨 채로 굳어있게(frozen) 될 텐데 따라서 개별 스핀은 $\left< S_i \right> \neq 0$일 것이다. 무질서에 대해 평균을 내면 $m=\left[ \big\langle S_{i} \big\rangle \right]=0$이 될 테지만, $q = \left[ \left< S_i \right>^2 \right]$은 양수인 값을 가지는 것이 가능하다. 상자성 상(paramagnetic phase)에서는 $m=q=0$이고, 강자성 상(ferromagnetic phase)에서는 $m$과 $q$ 모두 0이 아닌 값을 가진다.