イジング模型

統計力学において**イジング模型**（イジングもけい、、イジングモデルとも言う）とは、二つの配位状態をとる格子点から構成され、最隣接する格子点のみの相互作用を考慮する格子模型である。二つの配位状態をスピンとする磁性体のモデルだが、二元合金、格子気体のモデルにも等価である。

スピン系のモデルとしては非常に単純化されたモデルであるが、相転移現象を記述可能なモデルであり、多くの物理学者によって研究されてきた。単純なモデルであるため厳密な解析が可能であり、特に外部磁場の無い二次元イジング模型は厳密解が得られる可解格子模型の一種である。

イジング模型は1920年にドイツの物理学者によって提案された。イジング模型という名前はレンツの博士課程の指導学生でありこの模型の研究を行っていたエルンスト・イジングに因んでいる。1944年にラルス・オンサーガーによって与えられた二次元イジング模型の厳密解は統計力学における金字塔の一つとされる。

## 概要
磁性体のモデルとして、 -次元空間の格子点に上向きと下向きの2状態をとるスピンが配置された格子模型を考える。
を  番目の格子点におけるスピンの状態を示す変数とし、が上向きのスピン、が下向きのスピンに対応するものとする。格子点の総数は  個とし、一つの格子点に最近接する格子点の数を  個とする。例えば、１次元格子では、2次元正方格子では 、3次元立方格子では である。

を２つの格子点間における交換相互作用、は格子点  における外部磁場とする。このとき、イジング模型のハミルトニアンは次式で与えられる。

: \mathcal{H} =  -  \sum_{\left\langle i,j \right\rangle}J_{ij} \, \sigma_i \cdot \sigma_j -  \sum_{i} h_{i} \, \sigma_i

第1項目は最近接する格子点におけるスピン間の相互作用のエネルギーを表す。記号は最近接する格子点のペアについての和であることを意味し、の和は個の項の和になる。
の場合を強磁性相互作用、の場合を反強磁性相互作用という。強磁性相互作用では最近接する格子点 のスピンのペアが同じ向きに揃い、となるとエネルギーは  だけ下がる。そのため、エネルギーが最も低い基底状態は全てのスピンの向きが揃った状態となる。一方、反強磁性相互作用では最近接する格子点のスピンのペアが異なる向きをとり、となるとエネルギーは  だけ下がる。第2項目は外部磁場に対するエネルギーを表す。格子点において、スピンの向き（符号）が外部磁場の向き（符号）と揃うと、エネルギーは  だけ下がる。

特に格子点上で交換相互作用と外部磁場を一定値とする一様なケースでは、イジング模型のハミルトニアンは

: \mathcal{H} =  - J \sum_{\left\langle i,j \right\rangle} \sigma_i \cdot \sigma_j -  h \sum_{i}  \sigma_i

となる。

統計力学において、温度の平衡状態での系の熱力学的な性質は分配関数  から求まる。分配関数は系の取りうる全ての状態についてのボルツマン因子の足し合わせで与えられる。 個の格子点をもつイジング模型においては、格子点のスピン変数がの値をとる個の状態が存在し、分配関数は

:
\begin{align}
Z ( \beta ,N ) &=\sum_{\sigma_1=\pm 1}\cdots\sum_{\sigma_N=\pm 1}e^{-\beta \mathcal{H}} \\
& =\sum_{\sigma_1=\pm 1}\cdots\sum_{\sigma_N=\pm 1}
\exp \biggl ( \beta \sum_{\left\langle i,j \right\rangle}
J_{ij} \sigma_i  \sigma_j  +\beta \sum_{i}h_{i} \sigma_i \biggr)
\end{align}


となる。分配関数から自由エネルギー、磁化、帯磁率が求まる。

一般に相互作用を含むモデルでは分配関数を求めることは困難であるが、交換相互作用と外部磁場を一様とする設定において、イジング模型では1次元のケース、外部磁場のない2次元のケースについては、厳密に分配関数を求めることが可能である。

エルンスト・イジングによる1925年の解析の段階で、一次元系での厳密な解は求められていて、有限温度での相転移を起こさないことが示されていた。その後、1944年にラルス・オンサーガーが二次元イジング模型の厳密解を求めた。これは相転移を起こし、この結果は相転移現象の記述と理解のために大変重要な役割を果たしている。オンサーガーの方法以外にも外部磁場のない二次元イジング模型の厳密解を求める方法がいくつか知られている。しかし、外部磁場のある場合の厳密解は得られていない。

三次元イジング模型の厳密解は知られていないが、共形ブートストラップを用いて解析的に臨界指数を求める試みがなされている
。

厳密解以外にも平均場近似や繰り込み群、級数展開（低温展開、高温展開）の手法などによる近似解が知られている。と、これらを用いた数値計算手段を使って近似的に解かれる。

この模型は、結晶表面のラフニング転移や合金の規則‐不規則（秩序‐無秩序）転移、異方性の大きな磁性の問題などに応用されている。

## 一般化
イジング模型は最近接する格子点以外にも任意の格子点間の相互作用を考慮する形に拡張することができる
。このとき、ハミルトニアン&#8459;は
: 
\begin{align}
\mathcal{H} 
&=  - \sum_{(i,j)} J_{ij} \, \sigma_i \cdot \sigma_j - \sum_{i} h_i \sigma_i \\
&=  - \frac{1}{2} \sum_{i=1}^{N}\sum_{j=1}^{N} 
J_{ij} \, \sigma_i \cdot \sigma_j - \sum_{i=1}^{N} h_i \sigma_i
\end{align}

となる。

より一般にイジング模型は、無向グラフ上で定義することができる。頂点を},頂点同士を繋ぐ辺をとする無向グラフにおいて、イジング模型のハミルトニアンは

: \mathcal{H} =  - \sum_{(i,j) \isin E} J_{ij} \, \sigma_i \cdot \sigma_j - \sum_{i\isin V}h_i \sigma_i


となる。

## 歴史
イジング模型はヴィルヘルム・レンツによって、磁性体のモデルとして考案され、1920年の論文の中で与えられた。1920年にレンツはロストック大学にいたが、1921年にハンブルク大学の教授に着任した。ハンブルク大学においてエルンスト・イジングはレンツの博士課程の学生であった。博士論文の中でレンツの考案したイジング模型について取り組み、博士論文は1924年に提出された。イジングは、一次元のケースについて、イジング模型の分配関数を厳密に求め、強磁性体の相転移が起きないことを示した。その結果をもって、3次元の場合についても相転移が起きないとする誤った結論に至った。イジングの博士論文の要約は1925年に論文として出版された。

1928年にヴェルナー・ハイゼンベルクはそれまでに得られていた量子力学の知見から、強磁性の相互作用の起源が量子力学的交換相互作用であるとする論文を出した。この中でハイゼンベルクはスピンがベクトル型相互作用するモデルを考案した。論文の序論において、イジングの結果も引用し、「他の困難についてはレンツとイジングが詳しく論じており、イジングは鎖状に隣り合う2つの原子の間に向きが揃った十分大きな力が働くという仮定も、強磁性を作り出すには十分でないことを示すことに成功した」と記している。

1936年にルドルフ・パイエルスは「強磁性のイジング模型について」という論文を出した。パイエルスは合金の秩序無秩序転移の議論を参考にしつつ、強磁性体の相転移について論じた。そして、十分低温であれば、2次元ではイジング模型は自発磁化が存在することを示した。その結果を基に3次元でも自発磁化が存在すると結論した。パイエルスの議論は、2次元格子において上向きのスピンが集まった領域と下向きのスピンが集まった領域の複数の集まりに分け、その境界（閉曲線）の長さから磁化の下限を評価するものである。この相転移の存在についての評価手法はパイエルスの議論と呼ばれる。なお、パイエルスの証明には不完全な部分があり、完全な証明は1964年にロバート・グリフィスによって与えられた。

1938年にジョン・G・カークウッドは分配関数を温度の逆数のべき乗で系統的に展開する方法を与えた。これはトルバルド・ティエレによって、研究されたキュムラント（半不変式）の性質に基づくものであった。

1941年にヘンリク・アンソニー・クラマースとグレゴリー・ワニエは2つの論文を出した。一つ目の論文で、彼らはイジング模型の分配関数が、ある種の行列の最大固有値から計算できることを示した。。この手法は統計力学において、転送行列の方法と呼ばれる。二つ目の論文でクラマースとワニエは2次元イジング模型について、相転移が起こるキューリー温度の値を求めた。クラマースとワニアは分配関数の高温と低温での級数展開について成り立つ対称性を導いた。その上で相転移が存在することを仮定し、転移温度を計算した。この対称性は2次元正方格子の表格子の温度と裏格子の温度についての対称性であり、クラマース＝ワニエ双対性と呼ばれる。

2次元イジング模型の厳密な結果を与えたのはラルス・オンサーガーである。オンサーガーは1942年のニューヨーク科学アカデミーの会合で外部磁場が無い場合の2次元イジング模型の厳密解を求めたことを報告した。そして、その結果は1944年に論文として出版された。

## 対称性
イジング模型はスピン反転対称性や副格子対称性と呼ばれる対称性をもつ。

### スピン反転対称性
各格子点上のスピン変数の組をまとめて、と表す。全ての格子点のスピン変数の向きを反転させる変換を行うと、ハミルトニアンは

: 
\mathcal{H} (h,\{-\sigma_i\})=- J\sum_{(i,j)}  \, \sigma_i \cdot \sigma_j - (-h)\sum_{i}\sigma_i
=\mathcal{H} (-h,\{\sigma_i\})
 

となり、これは外部磁場の向きの反転と等価である。分配関数については、の取りうる全ての状態についてのボルツマン因子の和との取りうる全ての状態についてのボルツマン因子の和は等価であり、

: 
Z(h)=Z(-h)
 

が成りたつ。その結果、単位スピン当たりの自由エネルギーについても

: 
f(h)=f(-h)
 

も成り立つ。これらの対称性をスピン反転対称性または対称性という。

## 1次元モデル
相互作用の減衰が *α* > 1 で  J_{ij} \sim |i-j|^{-\alpha}  であれば、熱力学的極限が存在する。

* 1  J_{ij} \sim |i-j|^{-\alpha}  の場合について、ダイソン(Dyson)は階層を比較することにより充分小さな温度で相転移があることを証明した。

* **強磁性**の相互作用   J_{ij} \sim |i-j|^{-2}  の場合について、フレーリッヒ(Fröhlich)とスペンサー(Spencer)は（階層の場合と対照的に）充分小さな温度で相転移があることを示した。

**α* > 2 の相互作用  J_{ij} \sim |i-j|^{-\alpha}  の場合（このことは有限の範囲の相互作用を意味する）においては、自由エネルギー(free energy)が熱力学パラメータに対して解析的であるので、正の温度（有限の *β*）に対して相転移がない。

* **近接相互作用**の場合についてはイジング(E. Ising)がモデルの完全解を示した。任意の正の温度（有限の *β*）で、自由エネルギーは熱力学的パラメータの中で解析的であり、省略された 2点相関函数は指数的に急速に減少する。温度 0 （*β* が無限大）では、第二種の相転移がある。自由エネルギーは無限大となり、領略された 2点スピンの相関函数は減少しない（定数のままである）。従って、*T* = 0 はこの場合の臨界温度であり、スケーリング公式を満す。

### イジングによる完全解
（周期的境界条件、または、自由境界条件）近接相互作用の場合、完全解が存在する。周期境界条件を持つ格子 *L* の上の1次元イジングモデルのエネルギーは、
:
\mathcal{H}(\sigma)
= -J\sum_{i=1,\ldots,L} \sigma_i \sigma_{i+1} - h \sum_i \sigma_i

である。ここに *J* と *h* は、この単純化された場合には *J* は定数で近隣間の相互作用の強さを表し、*h* は格子に適用された定数の外場であるので、任意の数値で問題ない。従って、自由エネルギーは、
:
\begin{align}
f(\beta, h)
&= -\lim_{L\to \infty} \frac{1}{\beta L} \ln (Z(\beta)) \\
&= -\frac{1}{\beta} \ln\left(e^{\beta J} \cosh \beta h+\sqrt{e^{2\beta J}(\sinh\beta h)^2+e^{-2\beta J}}\right)
\end{align}

であり、スピン-スピン相関函数は、
:
\langle \sigma_i \sigma_j\rangle-\langle \sigma_i \rangle\langle\sigma_j\rangle
= C(\beta)e^{-c(\beta)|i-j|}


である。ここに *C*(*β*) と *c*(*β*) は *T* > 0 の正の値の函数である。しかし、*T* → 0 とすると、逆の相関の長さ *c*(*β*) は 0 となる。

## 応用
イジング模型は強磁性体や反強磁性体のモデルではあるが、二元合金や格子気体のモデルとも等価である
。また、イジング模型は不規則磁性体の秩序相であるスピングラスのモデルにも用いられる。スピングラスでは、強磁性と反強磁性の相互作用が空間的にランダムに入り混じったイジング模型が用いられる。スピングラス理論における解析手法は、ニューラルネットワーク（神経回路網）における連想記憶の理論や組合せ最適化問題にも適用されており、これらの分野においてもイジング模型が応用されている。

### 二元合金
2種類の金属原子が格子点上に配置された二元合金の系を考える。格子点の総数をとし、金属原子の個数を、金属原子の個数をとする。原子間の相互作用としては、最近接格子点に同士が並んだ時に、同士が並んだ時に、とが並んだ時にだけのポテンシャルエネルギーをもつとする。また、は同士が最近接する格子点のペア数、は同士が最近接する格子点のペア数、はとが最近接する格子点のペア数とする。系のポテンシャルエネルギーは

:
E=\phi_{AA}N_{AA}+\phi_{BB}N_{BB}+\phi_{AB}N_{AB}


となる。は独立でなく、
: N_{A}+N_{B}=N 
: zN_{A}=2N_{AA}+N_{AB} 
: zN_{B}=2N_{BB}+N_{AB} 
の関係を満たす。ここで  は一つの格子点の最近接する格子点の数である。格子点  における変数
を、金属が占有しているときに、金属が占有しているときにの値をとるものと定義する。このとき、

: \sum_{i}\sigma_{i}=N_{A}-N_{B} 
: \sum_{\left\langle i,j \right\rangle}\sigma_{i} \cdot \sigma_{j}=N_{AA}+N_{BB}-N_{AB} 

であるから、系のポテンシャルエネルギーは

:
E= - \frac{1}{4} \left (2\phi_{AB}-\phi_{AA}-\phi_{BB} \right ) \sum_{\left\langle i,j \right\rangle} \sigma_i \cdot \sigma_j 
+ \frac{z}{4}\left (\phi_{AA}-\phi_{BB} \right ) \sum_{i}  \sigma_i
+ \frac{z}{4}\left (2\phi_{AB}+\phi_{AA}+\phi_{BB} \right )N


と書き表せる。これは交換相互作用を
:
J=  \frac{1}{4} \left (2\phi_{AB}-\phi_{AA}-\phi_{BB} \right ) 


とし、外部磁場を
:
h=  \frac{z}{4}\left (\phi_{AA}-\phi_{BB} \right ) 

とするイジング模型と定数項を除いて等価である。

### スピングラス
常磁性体金属に微量の磁性元素を添加した磁性希薄合金では、スピングラスと呼ばれる磁気的秩序相が存在する。エドワーズ・アンダーソン模型では、正負の値を取りえる磁気的相互作用が空間的にランダムに分布した不規則磁性体としてスピングラスを扱う。このモデルでは、系のハミルトニアンはランダムな磁気的相互作用を持つイジング模型

: \mathcal{H} = - \sum_{\left\langle i,j \right\rangle} J_{ij} \, \sigma_i \cdot \sigma_j 

である。相互作用の項の和は最近接する格子点のペアについてとる。は強磁性的（）と反強磁性的（）の両者の値を取りえる確率変数である。の分布としては、確率密度関数が

: 
P(J_{ij})=\frac{1}{\sqrt{2\pi J^2}}
\exp{\left (-\frac{1}{2J^2}\biggl(J_{ij}-J_{0}\biggr)^2 \right) } 


である平均、分散のガウス分布や

: 
P(J_{ij})=p \, \delta(J_{ij}-J)+(1-p) \, \delta(J_{ij}+J)


と確率で値をとり、確率で値をとる分布が用いられる。記号はデルタ関数である。

一方、スピングラスのシェリントン・カークパトリック模型は、空間的にランダムな相互作用が全ての格子点のペアについてわたる無限レンジであるとするモデルである。このモデルでは、系のハミルトニアンはランダムに分布する相互作用を無限レンジとするイジング模型

: \mathcal{H} = - \sum_{\left ( i,j \right)} J_{ij} \, \sigma_i \cdot \sigma_j 

である。確率変数は確率密度関数が

: 
P(J_{ij})=\left(\frac{N}{2\pi J^2}\right)^{\frac{1}{2}} \cdot
\exp{\left (-\frac{N}{2J^2}\biggl(J_{ij}-\frac{J_{0}}{N}\biggr)^2 \right) } 


である平均、分散のガウス分布に従う。

### 組合せ最適化問題
組合せ最適化問題では、与えられた制約条件の下、コスト関数を最小化する組合せの解を探索する。一般に要素数が増大すると、組合せ数が指数関数的に増大する組合せ爆発が生じ、解の探索は困難になる。シミュレーティド・アニーリングや量子アニーリングの手法では、組合せ最適化問題をイジング模型の問題に帰着させ、解の候補の探索を行い、近似的な解を与える。このとき、イジング模型のエネルギーを最小とするスピンの配位状態が解となる。これらの手法を実装したハードウェアをイジングマシンという。

## 脚注

### 注釈

### 出典

## 参考文献
*
*
*
*
*
*
*
*
* 黒田耕嗣：「秩序・無秩序の世界：イジング・モデルと相転移のはなし」、丸善、ISBN 4-621-03178-3 (1987年6月30日)。
*
*
*
* 「イジング模型の動的臨界指数の下限z ≥ 2を厳密に証明:100年来の未解決問題に対しての理論的進展」(東京大学、2025年5月28日)

## 関連記事
*XY模型
*量子力学
*物性物理学