イジング模型
統計力学において、イジング模型(英: Ising model、イジングモデルとも言う)とは二つの配位状態をとる格子点から構成され、最隣接する格子点のみの相互作用を考慮する格子模型。強磁性体の模型(モデル)であるとともに、二元合金、格子気体の模型としても用いられる。スピン系のモデルとしては非常に単純化されたモデルであるが、相転移現象を記述可能なモデルであり、多くの物理学者によって、研究されてきた[1]。また、この単純化された性質により、厳密な解析が可能であり、特に外部磁場の無い二次元イジング模型は、厳密解が得られる可解格子模型の一種である。1920年にドイツの物理学者ヴィルヘルム・レンツによって、提案された[2]。イジング模型の名は、レンツの博士課程の指導学生であり、その研究を行ったエルンスト・イジングの名前に因む[3]。1944年に、ラルス・オンサーガーが二次元イジング模型の厳密解を求め、相転移が起きることを示したが、この結果は、統計力学における金字塔の一つとされる[4]。
概要
ハイゼンベルク模型をより簡単化したもので、そのハミルトニアンは、
[math] H = - J \sum_{\left\langle i,j \right\rangle} \sigma_i \cdot \sigma_j [/math]
である。σi は(結晶)格子点 i 上のスピン。自由度は上向き (+1) と下向き (−1) のみである。J は最隣接スピン間の相互作用によるエネルギー(交換相互作用エネルギー)である。[math] \left\langle i,j \right\rangle [/math]は最隣接格子点のみ和を取ることを意味する。ここで、隣り合う格子点上のスピン同士が、お互い上向き (+1 × +1) または下向き (−1 × −1) の場合(スピンが平行)、−J/2、互いに逆向き(上向きと下向き:1 × −1または−1 × 1、スピン反平行)の場合、J/2 となる。
イジングの提案の段階で、一次元(格子系)での厳密な解は求められていて、有限温度での相転移を起こさないことが示されていた。その後、1944年にラルス・オンサーガーが二次元イジング模型の厳密解を求めた。これは相転移を起こし、この結果は、相転移現象の記述、理解のために大変重要な役割を果たしている。二次元の磁場の無い場合のこの模型の厳密解はオンサーガーの解法以外にもいくつかの方法が示されている.また,外部磁場が印加されたモデルの厳密解は得られていない.
三次元に関しての厳密解は現在求められていないが、共形ブートストラップを用いて解析的に臨界指数を求める試みがなされている[5] [6]。
厳密解が求められるのは、特殊な場合で多くの場合、平均場近似、繰り込み群、級数展開(低温展開、高温展開)の手法などと、これらを用いた数値計算手段を使って近似的に解かれる。
この模型(モデル)は、合金の規則‐不規則(秩序‐無秩序)転移や、異方性の大きな磁性の問題などに適用されている。
1次元モデル
相互作用が崩壊が α > 1 で [math] J_{ij} \sim |i-j|^{-\alpha} [/math] であれば、熱力学的極限が存在する[7]。
- 1 < α < 2 で強磁性の相互作用 [math] J_{ij} \sim |i-j|^{-\alpha} [/math] の場合、ダイソン(Dyson)は、階層を比較することにより、充分小さな温度で相転移があることを証明した[8]。
- 強磁性の相互作用 [math] J_{ij} \sim |i-j|^{-2} [/math] の場合、フローリッヒ(Fröhlich)とスペンサー(Spencer)は、(階層の場合と対照的に)充分小さな温度で相転移があることを示した[9]。
- α > 2 の相互作用 [math] J_{ij} \sim |i-j|^{-\alpha} [/math] の場合(このことは有限の範囲の相互作用を意味する)、自由エネルギー(free energy)が熱力学パラメータに対して解析的であるので、正の温度(つまり、有限の β)に対して相転移がない[7]。
- 近接相互作用の場合、イジング(E. Ising)は、モデルの完全解を示した。任意の正の温度(つまり、有限の β)で、自由エネルギーは熱力学的パラメータの中で解析的であり、省略された 2点相関函数は指数的に急速に減少する。温度 0 (つまり、β が無限大)では、第二種の相転移がある。自由エネルギーは無限大となり、領略された 2点スピンの相関函数は減少しない(定数のままである)。従って、T = 0 はこの場合の臨界温度であり、スケーリング公式を満す[10]。
イジングによる完全解
(周期的境界条件、または、自由境界条件)近接相互作用の場合、完全解が存在する。周期境界条件を持つ格子 L の上の 1次元イジングモデルのエネルギーは、
- [math]H(\sigma)=-J\sum_{i=1,\ldots,L} \sigma_i \sigma_{i+1} - h \sum_i \sigma_i[/math]
である。ここに J と h は、この単純化された場合には J は定数で近隣間の相互作用の強さを表し、h は格子に適用された定数の外場であるので、任意の数値で問題ない。従って、自由エネルギー(free energy)は、
- [math]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) [/math]
であり、スピン-スピン相関函数は、
- [math]\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 となる。
脚注
- ↑ Stephen G. Brush, Rev. Mod. Phys., 39, p.883 (1967)
- ↑ W. Lenz, "Beiträge zum Verständnis der magnetischen Eigenschaften in festen Körpern", Phys. ZS., 21, p. 613 (1920)
- ↑ E. Ising, "Beitrag zur theorie des ferromagnetismus," Z. Physik, 31, p. 253 (1925)
- ↑ Somendra M. Bhattacharjee et al., Curr.Sci., 69 p. 816 (1995)
- ↑ Sheer El-Showk et al., Rhys.Rev.D 86, 025022 (2012)
- ↑ Sheer El-Showk et al., J. Stat. Phys. 157, p. 869-914 (2014)
- ↑ 7.0 7.1 (1969) Statistical Mechanics:Rigorous Results.. New York: W.A. Benjamin Inc..
- ↑ Dyson, F.J. (1969). “Existence of a phase-transition in a one-dimensional Ising ferromagnet”. Comm. Math. Phys. 12: 91–107. Bibcode 1969CMaPh..12...91D. doi:10.1007/BF01645907.
- ↑ Fröhlich, J.; Spencer, T. (1982). “The phase transition in the one-dimensional Ising model with 1/r 2 interaction energy.”. Comm. Math. Phys. 84: 87–101. doi:10.1007/BF01208373 .
- ↑ Baxter, Rodney J. (1982), Exactly solved models in statistical mechanics*, London: Academic Press Inc. [Harcourt Brace Jovanovich Publishers], ISBN 978-0-12-083180-7, MR 690578
参考文献
- Lars Onsager: "Crystal statistics. I. A two-dimensional model with an order-disorder transition", Phys. Rev. Vol.65, pp.117-149 (1944).
- Stephen G. Brush, "History of Lentz-Ising Model," Rev. Mod. Phys., 39, p.883 (1967) doi:10.1103/RevModPhys.39.883
- Somendra M. Bhattacharjee, Avinash Khare, "Fifty Years of the Exact Solution of the Two-Dimensional Ising Model by Onsager," Curr.Sci., 69, p. 816 (1995) テンプレート:ArXiv
- Martin Niss, "History of the Lenz-Ising Model 1920–1950: From Ferromagnetic to Cooperative Phenomena," Arch. Hist. Exact Sci., 59, p. 267 doi:10.1007/s00407-004-0088-3
- 小口武彦 『磁性体の統計理論 』 裳華房 (1970) ISBN 978-4785323127
- I. A. Stepanov. Exact Solutions of the One-Dimensional, Two-Dimensional, and Three-Dimensional Ising Models. – Nano Science and Nano Technology: An Indian Journal. 2012. Vol. 6. No 3. 118 - 122. The paper is on the Journal’s website with a free access.