イジング模型のソースを表示

[[File:2D ising model on lattice.svg|thumb|upright=1.4|２次元格子上のイジング模型]]
{{統計力学}}
[[統計力学]]において'''イジング模型'''（イジングもけい、{{lang-en-short|Ising model}}、イジングモデルとも言う）とは、二つの配位状態をとる格子点から構成され、最隣接する格子点のみの相互作用を考慮する[[格子模型]]である<ref name ="huang1991_ch14">[[#huang1991|K. Huang (1991), chapter14]]</ref>。二つの配位状態を[[スピン角運動量|スピン]]とする[[磁性体]]のモデルだが、[[二元合金]]、[[格子気体]]のモデルにも等価である<ref name ="huang1991_ch14"></ref>。

スピン系のモデルとしては非常に単純化されたモデルであるが、[[相転移]]現象を記述可能なモデルであり、多くの物理学者によって研究されてきた<ref name ="brush1967">[[#brush1967|Stephen G. Brush, ''Rev. Mod. Phys.'', '''39''', p.883 (1967)]]</ref><ref name ="niss2005">[[#niss2005|Martin Niss, ''Arch. Hist. Exact Sci.'', '''59''', p.267 (2005)]]</ref><ref name ="niss2009">[[#niss2009|Martin Niss, ''Arch. Hist. Exact Sci.'', '''63''', p.243 (2009)]]</ref><ref name ="niss2011">[[#niss2011|Martin Niss, ''Arch. Hist. Exact Sci.'', '''65''', p.625 (2011)]]</ref>。単純なモデルであるため厳密な解析が可能であり、特に外部磁場の無い二次元イジング模型は厳密解が得られる[[可解格子]]模型の一種である。

イジング模型は1920年にドイツの物理学者{{仮リンク|ヴィルヘルム・レンツ|en|Wilhelm Lenz}}によって提案された<ref name ="lenz1920">{{cite journal |author1=W. Lenz|authorlink1=ヴィルヘルム・レンツ||year=1920 |title=Beitrag zum Verständnis der magnetischen Erscheinungen in festen Körpern|journal=Physik. Z.|volume=21 |issue= |pages=613–615 |url=https://www.physik.uni-rostock.de/storages/uni-rostock/Alle_MNF/Physik/Historisches/Kalenderblaetter_Physik/KB_2013_06_Lenz/Lenz_1920.pdf}}</ref><ref name ="brush1967"></ref>。イジング模型という名前はレンツの博士課程の指導学生でありこの模型の研究を行っていた[[エルンスト・イジング]]に因んでいる<ref name ="ising1925">{{cite journal |author1=E. Ising|authorlink1=エルンスト・イジング||year=1925 |title=Beitrag zur Theorie des Ferromagnetismus|journal= Z. Physik|volume=31 |issue= |pages=253–258|doi=10.1007/BF02980577}}</ref><ref name ="brush1967"></ref>。1944年に[[ラルス・オンサーガー]]によって与えられた二次元イジング模型の厳密解は統計力学における金字塔の一つとされる<ref>Somendra M. Bhattacharjee ''et al.'',  ''Curr.Sci.'', '''69'''  p. 816 (1995)</ref>。

== 概要 ==
磁性体のモデルとして、{{mvar|d}} -次元空間の格子点に上向きと下向きの2状態をとる[[スピン角運動量|スピン]]が配置された格子模型を考える。
{{math|''&sigma;<sub>i</sub>''{{=}}&plusmn;1}}を {{mvar|i}} 番目の格子点におけるスピンの状態を示す変数とし、{{math|&plus;1}}が上向きのスピン、{{math|&minus;1}}が下向きのスピンに対応するものとする。格子点の総数は {{mvar|N}} 個とし、一つの格子点に最近接する格子点の数を {{mvar|z}} 個とする。例えば、１次元格子では{{math|''z'' {{=}}2}}、2次元正方格子では {{math|''z'' {{=}}4}}、3次元立方格子では {{math|''z'' {{=}}6}}である。

{{math|''J<sub>ij</sub>''}}を２つの格子点{{math|''i, j''}}間における[[交換相互作用]]、{{math|''h<sub>i</sub>''}}は格子点 {{mvar|i}} における外部[[磁場]]とする。このとき、イジング模型の[[ハミルトニアン]]は次式で与えられる<ref group="注">行列要素{{math|{{tilde|''J''}}''<sub>ij</sub>''}}が格子点{{math|''i, j''}}が最近接するときのみに {{math|''J<sub>ij</sub>''}} の値をとり、それ以外は{{math|0}}とすると、ハミルトニアンは
:<math> \mathcal{H} =  - \frac{1}{2} \sum_{i=1}^{N}\sum_{j=1}^{N} \tilde
J_{ij} \, \sigma_i \cdot \sigma_j - \sum_{i=1}^{N} h_i \sigma_i
</math>
:とも表せる。</ref>。

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

第1項目は最近接する格子点におけるスピン間の相互作用のエネルギーを表す。記号{{math|{{angbr|''i, j''}}}}は最近接する格子点のペアについての和であることを意味し、{{math|{{angbr|''i,j''}}}}の和は{{math|''zN''/2}}個の項の和になる。
{{math|''J<sub>ij</sub>'' >0}}の場合を強磁性相互作用、{{math|''J<sub>ij</sub>'' <0}}の場合を反強磁性相互作用という。強磁性相互作用では最近接する格子点 {{math|''i,j''}}のスピンのペアが同じ向きに揃い、{{math|''&sigma;<sub>i</sub>&middot;&sigma;<sub>j</sub>''{{=}}&plus;1}}となるとエネルギーは {{math|''J<sub>ij</sub>''}} だけ下がる。そのため、エネルギーが最も低い[[基底状態]]は全てのスピンの向きが揃った状態となる。一方、反強磁性相互作用では最近接する格子点のスピンのペアが異なる向きをとり、{{math|''&sigma;<sub>i</sub>&middot;&sigma;<sub>j</sub>''{{=}}&minus;1}}となるとエネルギーは {{math|{{abs|''J<sub>ij</sub>''}}}} だけ下がる。第2項目は外部磁場に対するエネルギーを表す。格子点{{math|''i''}}において、スピンの向き（符号）が外部磁場の向き（符号）と揃うと、エネルギーは {{math|{{abs|''h<sub>i</sub>''}}}} だけ下がる。

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

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

となる。

統計力学において、温度{{mvar|T}}の平衡状態での系の熱力学的な性質は[[分配関数]] {{mvar|Z}} から求まる。分配関数は系の取りうる全ての状態についての[[ボルツマン因子]]{{math|e<sup>&minus;''&beta;H</sup>''}}の足し合わせで与えられる。{{mvar|N}} 個の格子点をもつイジング模型においては、格子点のスピン変数が{{math|''&sigma;''{{=}}&plusmn;1}}の値をとる{{math|2<sup>''N''</sup>}}個の状態が存在し、分配関数は

:<math>
\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}
</math>

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

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

エルンスト・イジングによる1925年の解析の段階で、一次元系での厳密な解は求められていて、有限温度での[[相転移]]を起こさないことが示されていた<ref name ="ising1925"></ref>。その後、1944年に[[ラルス・オンサーガー]]が二次元イジング模型の厳密解を求めた<ref name ="onsager1944">{{Cite journal|author=Lars Onsager|authorlink1=ラルス・オンサーガー|title=Crystal statistics. I. A two-dimensional model with an order-disorder transition|journal=Phys. Rev.|year=1944|volume=65|pages=117-149|doi=10.1103/PhysRev.65.117}}</ref>。これは相転移を起こし、この結果は[[相転移]]現象の記述と理解のために大変重要な役割を果たしている。オンサーガーの方法以外にも外部磁場のない二次元イジング模型の厳密解を求める方法がいくつか知られている。しかし、外部磁場のある場合の厳密解は得られていない。

三次元イジング模型の厳密解は知られていないが、共形ブートストラップを用いて解析的に[[臨界指数]]を求める試みがなされている<ref>Sheer El-Showk ''et al.'',  ''Rhys.Rev.D'' '''86''', 025022 (2012)</ref>
<ref>Sheer El-Showk ''et al.'',  ''J. Stat. Phys.'' '''157''', p. 869-914 (2014)</ref>。

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

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

==一般化==
イジング模型は最近接する格子点以外にも任意の格子点間{{math|(''i'', ''j'')}}の相互作用を考慮する形に拡張することができる<ref name ="takahashi_nishimori2017_ch2">[[#nishimori1999|高橋、西森（2017）、第2章]]</ref><ref name ="goldenfeld1992_ch2">[[#goldenfeld1992|N.Goldenfeld (1992), chapter2]]</ref>
。このとき、ハミルトニアン&#8459;は
:<math> 
\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}
</math>
となる。

より一般にイジング模型は、[[無向グラフ]]上で定義することができる。頂点を{{math|''V''{{=}}{1,&hellip;, ''N''}}},頂点同士を繋ぐ辺を{{mvar|E}}とする無向グラフ{{math|''G''{{=}}(''V'', ''E'')}}において、イジング模型のハミルトニアンは

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

となる。

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

1928年に[[ヴェルナー・ハイゼンベルク]]はそれまでに得られていた量子力学の知見から、強磁性の相互作用の起源が量子力学的交換相互作用であるとする論文を出した<ref name ="heisenber1928">{{cite journal |author1=W. Heisenberg|authorlink1=ヴェルナー・ハイゼンベルク||year= |title=Zur Theorie des Ferromagnetismus|journal=1928 |volume=49  |issue= |pages=619-636|doi=10.1007/BF01328601}}</ref>。この中でハイゼンベルクはスピンがベクトル型相互作用するモデルを考案した。論文の序論において、イジングの結果も引用し、「他の困難についてはレンツとイジングが詳しく論じており、イジングは鎖状に隣り合う2つの原子の間に向きが揃った十分大きな力が働くという仮定も、強磁性を作り出すには十分でないことを示すことに成功した」と記している。

1936年に[[ルドルフ・パイエルス]]は「強磁性のイジング模型について」という論文を出した<ref>{{cite journal |author1=R. Peierls|authorlink1=ルドルフ・パイエルス|year=1936 |title=On Ising's model of ferromagnetism|journal=Proc. Camb. Soc. Philos.|volume=32 |issue=3 |pages=477–481 |doi=10.1017/S0305004100019174}}</ref>。パイエルスは合金の秩序無秩序転移の議論を参考にしつつ、強磁性体の相転移について論じた。そして、十分低温であれば、2次元ではイジング模型は自発磁化が存在することを示した。その結果を基に3次元でも自発磁化が存在すると結論した。パイエルスの議論は、2次元格子において上向きのスピンが集まった領域と下向きのスピンが集まった領域の複数の集まりに分け、その境界（閉曲線）の長さから磁化の下限を評価するものである。この相転移の存在についての評価手法はパイエルスの議論と呼ばれる。なお、パイエルスの証明には不完全な部分があり、完全な証明は1964年に[[ロバート・グリフィス]]によって与えられた<ref name ="griffiths1964>{{cite journal |author1=Robert B. Griffiths|authorlink1=ロバート・グリフィス||year=1964 |title=Peierls Proof of Spontaneous Magnetization in a Two-Dimensional Ising　Ferromagnet|journal=Phys. Rev. volume=136|issue=|pages=A437–|doi=10.1103/PhysRev.136.A437}}</ref>。

1938年に[[ジョン・G・カークウッド]]は分配関数を温度の逆数のべき乗で系統的に展開する方法を与えた<ref name ="kirwood1938>{{cite journal |author1=John G. Kirkwood |authorlink1=ジョン・G・カークウッド||year=1938|title=Order and Disorder in Binary Solid Solutions|journal=J. Chem. Phys.|volume=6 |issue= |pages=70-75 |doi=10.1063/1.1750205 }}</ref>。これは[[トルバルド・ティエレ]]によって、研究された[[キュムラント]]（半不変式）の性質に基づくものであった。

1941年に[[ヘンリク・アンソニー・クラマース]]と[[グレゴリー・ワニエ]]は2つの論文を出した。一つ目の論文で、彼らはイジング模型の分配関数が、ある種の行列の最大固有値から計算できることを示した。<ref name ="kramers_wannier1941a>{{cite journal |author1=H. A. Kramers|authorlink1=ヘンリク・アンソニー・クラマース|author2=G. H. Wannier|authorlink2=グレゴリー・ワニエ|year=1941 |title=Statitics of the Two-Dimensional Ferromagnet. Part I|journal=Physical Review |volume=60 |issue=3 |pages=252–262 |doi=10.1103/PhysRev.60.252}}</ref>。この手法は統計力学において、転送行列の方法と呼ばれる。二つ目の論文でクラマースとワニエは2次元イジング模型について、相転移が起こる[[キューリー温度]]の値を求めた<ref name ="kramers_wannier1941b>{{cite journal |author1=H. A. Kramers|authorlink1=ヘンリク・アンソニー・クラマース|author2=G. H. Wannier|authorlink2=グレゴリー・ワニエ|year=1941 |title=Statitics of the Two-Dimensional Ferromagnet. Part II|journal=Physical Review |volume=60 |issue=3 |pages=263–276 |doi=10.1103/PhysRev.60.263}}</ref>。クラマースとワニアは分配関数の高温と低温での級数展開について成り立つ対称性を導いた。その上で相転移が存在することを仮定し、転移温度を計算した。この対称性は2次元正方格子の表格子の温度と裏格子の温度についての対称性であり、[[クラマース＝ワニエ双対性]]と呼ばれる。

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

==対称性==
イジング模型はスピン反転対称性や副格子対称性と呼ばれる[[対称性 (物理学)|対称性]]をもつ<ref name ="takahashi_nishimori2017_ch2"></ref><ref name ="goldenfeld1992_ch2"></ref>。

===スピン反転対称性===
各格子点上のスピン変数{{mvar|&sigma;<sub>i</sub>}}の組をまとめて、{{math|{{mset|''&sigma;<sub>i</sub>''}}}}と表す。全ての格子点のスピン変数の向きを反転させる変換{{math|''&sigma;<sub>i</sub>''&rarr;&minus;''&sigma;<sub>i</sub>''}}を行うと、ハミルトニアンは

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

となり、これは外部磁場の向きの反転{{math|''h''&rarr;&minus;''h''}}と等価である。分配関数については、{{math|{{mset|''&sigma;<sub>i</sub>''}}}}の取りうる全ての状態についてのボルツマン因子{{math|e<sup>&minus;''&beta;&#8459;</sup>''}}の和と{{math|{{mset|&minus;''&sigma;<sub>i</sub>''}}}}の取りうる全ての状態についてのボルツマン因子{{math|e<sup>&minus;''&beta;&#8459;</sup>''}}の和は等価であり、

:<math> 
Z(h)=Z(-h)
</math> 

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

:<math> 
f(h)=f(-h)
</math> 

も成り立つ。これらの対称性をスピン反転対称性または{{math|Z<sub>2</sub>}}対称性という。

==1次元モデル==
相互作用の減衰が ''α'' > 1 で <math> J_{ij} \sim |i-j|^{-\alpha} </math> であれば、熱力学的極限が存在する<ref name="Ruelle">{{cite book|last1=Ruelle |title=Statistical Mechanics:Rigorous Results. |publisher=W.A. Benjamin Inc.|year=1969|location=New York}}</ref>。

* 1 < ''α'' < 2 で'''強磁性'''の相互作用 <math> J_{ij} \sim |i-j|^{-\alpha} </math> の場合について、ダイソン(Dyson)は階層を比較することにより充分小さな温度で相転移があることを証明した<ref>{{cite journal|last=Dyson|first=F.J.|title=Existence of a phase-transition in a one-dimensional Ising ferromagnet|journal=Comm. Math. Phys.|year=1969|volume=12|pages=91&ndash;107|doi=10.1007/BF01645907|bibcode = 1969CMaPh..12...91D }}</ref>。

* '''強磁性'''の相互作用  <math> J_{ij} \sim |i-j|^{-2} </math> の場合について、フレーリッヒ(Fröhlich)とスペンサー(Spencer)は（階層の場合と対照的に）充分小さな温度で相転移があることを示した<ref>{{cite journal|last1=Fröhlich|first1=J.|last2=Spencer|first2=T.|title=The phase transition in the one-dimensional Ising model with 1/r 2 interaction energy.|journal=Comm. Math. Phys.|year=1982|volume=84|url=http://www.springerlink.com/content/wu3782848714tt0l|doi=10.1007/BF01208373|pages=87–101}}</ref>。

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

* '''近接相互作用'''の場合についてはイジング(E. Ising)がモデルの完全解を示した。任意の正の温度（有限の ''β''）で、自由エネルギーは熱力学的パラメータの中で解析的であり、省略された 2点相関函数は指数的に急速に減少する。温度 0 （''β'' が無限大）では、第二種の相転移がある。自由エネルギーは無限大となり、領略された 2点スピンの相関函数は減少しない（定数のままである）。従って、''T'' = 0 はこの場合の臨界温度であり、スケーリング公式を満す<ref>{{citation | last1=Baxter | first1=Rodney J. | title=Exactly solved models in statistical mechanics* | url=http://tpsrv.anu.edu.au/Members/baxter/book | publisher=Academic Press Inc. [Harcourt Brace Jovanovich Publishers] | location=London | isbn=978-0-12-083180-7 | mr=690578 | year=1982}}</ref>。

===イジングによる完全解===
（周期的境界条件、または、自由境界条件）近接相互作用の場合、完全解が存在する。周期境界条件を持つ格子 ''L'' の上の1次元イジングモデルのエネルギーは、
:<math display="block">
\mathcal{H}(\sigma)
= -J\sum_{i=1,\ldots,L} \sigma_i \sigma_{i+1} - h \sum_i \sigma_i
</math>
である。ここに ''J'' と ''h'' は、この単純化された場合には ''J'' は定数で近隣間の相互作用の強さを表し、''h'' は格子に適用された定数の外場であるので、任意の数値で問題ない。従って、[[自由エネルギー]]は、
:<math display="block">
\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}
</math>
であり、スピン-スピン相関函数は、
:<math display="block">
\langle \sigma_i \sigma_j\rangle-\langle \sigma_i \rangle\langle\sigma_j\rangle
= C(\beta)e^{-c(\beta)|i-j|}
</math>

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

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

===二元合金===
2種類の金属原子{{math|''A'',''B''}}が格子点上に配置された二元合金の系を考える。格子点の総数を{{mvar|N}}とし、金属原子{{mvar|A}}の個数を{{mvar|N<sub>A</sub>}}、金属原子{{mvar|B}}の個数を{{mvar|N<sub>B</sub>}}とする。原子間の相互作用としては、最近接格子点に{{mvar|A}}同士が並んだ時に{{mvar|&varphi;<sub>AA</sub>}}、{{mvar|B}}同士が並んだ時に{{mvar|&varphi;<sub>BB</sub>}}、{{mvar|A}}と{{mvar|B}}が並んだ時に{{mvar|&varphi;<sub>AB</sub>}}だけのポテンシャルエネルギーをもつとする。また、{{mvar|N<sub>AA</sub>}}は{{mvar|A}}同士が最近接する格子点のペア数、{{mvar|N<sub>BB</sub>}}は{{mvar|B}}同士が最近接する格子点のペア数、{{mvar|N<sub>AB</sub>}}は{{mvar|A}}と{{mvar|B}}が最近接する格子点のペア数とする。系のポテンシャルエネルギーは

:<math>
E=\phi_{AA}N_{AA}+\phi_{BB}N_{BB}+\phi_{AB}N_{AB}
</math>

となる。{{math|''N<sub>AA</sub>''、''N<sub>BB</sub>''、''N<sub>AB</sub>''}}は独立でなく、
:<math> N_{A}+N_{B}=N </math>
:<math> zN_{A}=2N_{AA}+N_{AB} </math>
:<math> zN_{B}=2N_{BB}+N_{AB} </math>
の関係を満たす。ここで {{mvar|z}} は一つの格子点の最近接する格子点の数である。格子点 {{mvar|i}} における変数{{math|''&sigma;<sub>i</sub>''}}
を、金属{{mvar|A}}が占有しているときに{{math|''&sigma;<sub>i</sub>''{{=}}&plus;1}}、金属{{mvar|B}}が占有しているときに{{math|''&sigma;<sub>i</sub>''{{=}}&minus;1}}の値をとるものと定義する。このとき、

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

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

:<math>
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
</math>

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

とし、外部磁場{{mvar|h}}を
:<math>
h=  \frac{z}{4}\left (\phi_{AA}-\phi_{BB} \right ) 
</math>
とするイジング模型と定数項を除いて等価である。

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

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

である。相互作用の項の和は最近接する格子点のペア{{math|{{angbr|''i, j''}}}}についてとる。{{math|''J<sub>ij</sub>''}}は強磁性的（{{math|''J<sub>ij</sub>''>0}}）と反強磁性的（{{math|''J<sub>ij</sub>''<0}}）の両者の値を取りえる[[確率変数]]である。{{math|''J<sub>ij</sub>''}}の分布としては、[[確率密度関数]]が

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

である平均{{math|''J''<sub>0</sub>}}、分散{{math|''J''}}の[[ガウス分布]]や

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

と確率{{math|''p''}}で値{{math|''J'' (>0)}}をとり、確率{{math|1-''p''}}で値{{math|&minus;''J''}}をとる分布が用いられる。記号{{math|''&delta;''(''J<sub>ij</sub>'')}}は[[デルタ関数]]である。

一方、スピングラスのシェリントン・カークパトリック模型は、空間的にランダムな相互作用が全ての格子点のペア{{math|(''i, j'')}}についてわたる無限レンジであるとするモデルである<ref name ="nishimori1999_ch2"></ref>。このモデルでは、系のハミルトニアンはランダムに分布する相互作用を無限レンジとするイジング模型

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

である。確率変数{{math|''J<sub>ij</sub>''}}は確率密度関数が

:<math> 
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) } 
</math>

である平均{{math|''J''<sub>0</sub>/''N''}}、分散{{math|''J''/{{Sqrt|''N''}}}}のガウス分布に従う。

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

== 脚注 ==
{{脚注ヘルプ}}
=== 注釈 ===
{{reflist|group="注"}}
=== 出典 ===
{{reflist|30em}}

== 参考文献 ==
*{{Cite journal
|author=Stephen G. Brush
|authorlink1=
|title=History of Lentz-Ising Model
|journal=Rev. Mod. Phys.
|year=1967
|volume=39
|pages=883
|url=
|doi=10.1103/RevModPhys.39.883
|ref=brush1968}}
*{{Cite journal
|author=Martin Niss
|authorlink1=
|title=History of the Lenz-Ising Model 1920–1950: From Ferromagnetic to Cooperative Phenomena
|journal=Arch. Hist. Exact Sci.
|year=2005
|volume=59
|pages=267-318
|url=
|doi=10.1007/s00407-004-0088-3
|ref=niss2005}}
*{{cite journal 
|author1=Martin Niss|authorlink1=
||year=2009
|title=History of the Lenz–Ising Model 1950–1965: from irrelevance to relevanc
|journal=Arch. Hist. Exact Sci.
|volume=63
|issue=3 
|pages=243-287
|doi=10.1007/s00407-008-0039-5
|ref=niss2009}}
*{{cite journal 
|author1=Martin Niss|authorlink1=
|year=2011
|title=History of the Lenz-Ising model 1965-1971: the role of a simple model in understanding critical phenomena
|journal=Arch. Hist. Exact Sci.
|volume=65
|issue=6
|pages=625-658
|doi=10.1007/s00407-011-0086-1
|ref=niss2011}}
*{{Cite journal
|author1=Somendra M. Bhattacharjee
|authorlink1=
|author2=Avinash Khare
|title=Fifty Years of the Exact Solution of the Two-Dimensional Ising Model by Onsager,
|journal=
|year=1995
|volume=69
|pages=815
|url=
|doi=
|arXiv|cond-mat/9511003
|ref=bhattacharjee_khare1995}}
*{{Cite book |
|title= Statistical Mechanics
|author1=Kerson Huang
|authorlink1=Kerson Huang
|publisher=Wiley
|year=1991
|isbn=978-0471815181
|ref=huang1991}}
*{{Cite book |
|title= Lectures On Phase Transitions And The Renormalization Group
|author1= Nigel Goldenfeld
|authorlink1= Nigel Goldenfeld
|series= Frontiers in Physics
|publisher=Westview Press
|year=1992
|isbn=978-0201554090
|ref=goldenfeld1992}}
*{{Cite book |和書
|title= 磁性体の統計理論
|author1=小口武彦
|authorlink1=小口武彦
|series=物理学選書
|publisher=裳華房
|year=1970
|isbn=978-4785323127
|ref=oguchi1970}}
* 黒田耕嗣：「秩序・無秩序の世界：イジング・モデルと相転移のはなし」、丸善、ISBN 4-621-03178-3 (1987年6月30日)。
*{{Cite book |和書
|title= 統計力学
|author1=宮下精二
|authorlink1=宮下精二
|editor = [[植松恒夫]] (編集)、[[青山秀明]] (編集)、[[益川敏英]] (監修)
|series=基幹講座 物理学
|publisher=東京図書
|year=2020
|isbn=978-4489023446
|ref=miyashita2020}}
*{{Cite book |和書
|title=相転移・臨界現象とくりこみ群
|author1=高橋和孝
|authorlink1=高橋和孝
|author2=西森秀稔
|authorlink2=西森秀稔
|series=
|publisher=丸善出版
|year=2017
|isbn=978-4621301562
|ref=takahashi_nishimori2017}}
*{{Cite book |和書
|title=スピングラス理論と情報統計力学
|author1=西森秀稔
|authorlink1=西森秀稔
|series=新物理学選書
|publisher=岩波書店
|year=1999
|isbn=978-4000074148
|ref=nishimori1999}}

==関連記事==
*[[XY模型]]
*[[量子力学]]
*[[物性物理学]]

{{Normdaten}}
{{DEFAULTSORT:いしんくもけい}}
[[Category:磁気]]
[[Category:統計力学]]
[[Category:格子模型]]
[[Category:物理学のエポニム]]