HMFモデル

ハミルトニアン平均場モデルまたはHMFモデル[1] (英: Hamiltonian Mean Field model, HMF model) とは、物理学において多体系の研究に用いられる模型のひとつである。重力多体系を単純化した模型として宇宙物理学の観点から、また長距離相互作用のある系として統計力学の観点から研究されてきた。1次元系であるにもかかわらず相転移が起こる[2]ほか、violent relaxation、長期的に持続する metaequilibrium state、ゆっくりとした衝突緩和といった長距離相互作用系の興味深い性質を備えている[3]。

定義

HMF モデルは1次元 N {\displaystyle N} 粒子系のモデルであり、 θ i {\displaystyle \theta _{i}} , p i {\displaystyle p_{i}} をそれぞれ i {\displaystyle i} 番目の粒子の座標および運動量として、HMF モデルのハミルトニアンは

H = ∑ i = 1 N p i 2 2 + J 2 N ∑ i , j = 1 N [ 1 − cos ⁡ ( θ i − θ j ) ] {\displaystyle H=\sum _{i=1}^{N}{\frac {p_{i}^{2}}{2}}+{\frac {J}{2N}}\sum _{i,j=1}^{N}\left[1-\cos(\theta _{i}-\theta _{j})\right]}

により与えられる[2][4]。ただし座標 θ i {\displaystyle \theta _{i}} は θ i + 2 π {\displaystyle \theta _{i}+2\pi } と θ i {\displaystyle \theta _{i}} を同一視する。相互作用は J > 0 {\displaystyle J>0} のとき引力（強磁性）、 J < 0 {\displaystyle J<0} のとき斥力（反強磁性）である[2][5]。ただし反強磁性モデルは一様な状態がすべてのエネルギーで安定であり相転移が起きない[5]ため、以下では強磁性モデル ( J > 0 {\displaystyle J>0} ) について記述する。

しばしば秩序変数として1粒子あたりの磁化

M x = 1 N ∑ i = 1 N cos ⁡ θ i ,     M y = 1 N ∑ i = 1 N sin ⁡ θ i {\displaystyle M_{x}={\frac {1}{N}}\sum _{i=1}^{N}\cos \theta _{i},\ \ M_{y}={\frac {1}{N}}\sum _{i=1}^{N}\sin \theta _{i}}

が導入される[6][5][7]。これを用いると HMF モデルのハミルトニアンは

H = ∑ i = 1 N p i 2 2 + N J 2 ( 1 − M 2 ) {\displaystyle H=\sum _{i=1}^{N}{\frac {p_{i}^{2}}{2}}+{\frac {NJ}{2}}(1-M^{2})}

とも表せる[5]。

他の系との関係

統計力学において研究されるXY模型は、各格子点のスピンが2次元の単位ベクトル s = ( cos ⁡ θ i , sin ⁡ θ i ) {\displaystyle s=(\cos \theta _{i},\sin \theta _{i})} により表される模型であり、そのハミルトニアンは

H = − J ∑ ⟨ i , j ⟩ cos ⁡ ( θ i − θ j ) {\displaystyle H=-J\sum _{\langle i,j\rangle }\cos(\theta _{i}-\theta _{j})}

により与えられる[8]。ここに ⟨ i , j ⟩ {\displaystyle \langle i,j\rangle } は最近接格子の組に関する和であり、XY 模型では各スピンは最近接格子（1次元ならば隣接する2点）のスピンとのみ相互作用する。HMF モデルはその逆に各スピンが他のすべてのスピンと同じ強さで相互作用する（平均場近似）ものである[5][7][1]。

HMF モデルは1次元重力多体系において重力ポテンシャルのフーリエ級数表示を最低次で打ち切ったものに一致する[2]。すなわち、 N {\displaystyle N} 個の粒子の座標を θ i {\displaystyle \theta _{i}} とするとき、その系の重力ポテンシャル ψ {\displaystyle \psi } はポアソン方程式

∇ 2 ψ ( θ ) = k 2 ∑ i = 1 N [ δ ( θ − θ i ) − 1 2 π ] {\displaystyle \nabla ^{2}\psi (\theta )={\frac {k}{2}}\sum _{i=1}^{N}\left[\delta (\theta -\theta _{i})-{\frac {1}{2\pi }}\right]}

により定まる（ k {\displaystyle k} は定数、 δ ( x ) {\displaystyle \delta (x)} はディラックのデルタ関数）[2]。その解をフーリエ級数の形

ψ ( θ ) = k 2 ∑ i = 1 N ∑ n = 1 ∞ 1 − cos ⁡ n ( θ − θ i ) π n 2 {\displaystyle \psi (\theta )={\frac {k}{2}}\sum _{i=1}^{N}\sum _{n=1}^{\infty }{\frac {1-\cos n(\theta -\theta _{i})}{\pi n^{2}}}}

に表示するとき、最低次の n = 1 {\displaystyle n=1} の項のみを残す近似が HMF モデルである[2]。 θ i = θ j {\displaystyle \theta _{i}=\theta _{j}} での特異性を持たず系のサイズが有限であるため、HMF モデルは重力多体系特有の困難のないごく単純化した模型とみなすことができる[9]。

熱力学的性質

熱力学極限（粒子数 N → ∞ {\displaystyle N\to \infty } の極限）で系の統計分布は一体分布関数で記述できるようになる[10]。この1体分布関数は、粒子密度

ρ = A e β cos ⁡ θ {\displaystyle \rho =Ae^{\beta \cos \theta }}

（ A {\displaystyle A} , B {\displaystyle B} は定数）を導き、 B = 0 {\displaystyle B=0} ならば粒子は一様に分布し、 B ≠ 0 {\displaystyle B\neq 0} ならばクラスターを形成する[11]。定数 B {\displaystyle B} はセルフコンシステント条件

B = J 2 I 1 ( β B ) I 0 ( β B ) {\displaystyle B={\frac {J}{2}}{\frac {I_{1}(\beta B)}{I_{0}(\beta B)}}}

（ I n {\displaystyle I_{n}} は修正ベッセル関数）を満足し、これが非自明解 B ≠ 0 {\displaystyle B\neq 0} を持つ条件が

T < T c := J 2 {\displaystyle T<T_{c}:={\frac {J}{2}}}

と求まる[11]。すなわち、温度 T c {\displaystyle T_{c}} より低温のときにのみ熱平衡状態としてクラスター状態が可能である[11]。そして熱力学的安定性の要求から、低温側ではクラスター状態が安定であり、一方高温側では一様状態が安定であることが結論される[12]。この転移は2次相転移である[11]。

あるいは、この結果は以下の統計力学的な考察に基づいて導出することもできる。一般に長距離相互作用する系ではミクロカノニカルアンサンブルとカノニカルアンサンブルが等価ではなく異なる結果を導く可能性があるが、HMF モデルの場合には両者は等価である[13]。その1粒子あたりの自由エネルギー ϕ ( β ) {\displaystyle \phi (\beta )} ( β {\displaystyle \beta } は逆温度) は、極限 N → ∞ {\displaystyle N\to \infty } に対して

ϕ ( β ) = β 2 − 1 2 ln ⁡ 2 π + 1 2 ln ⁡ β + inf x ≥ 0 [ β x 2 2 − ln ⁡ I 0 ( β x ) ] {\displaystyle \phi (\beta )={\frac {\beta }{2}}-{\frac {1}{2}}\ln 2\pi +{\frac {1}{2}}\ln \beta +\inf _{x\geq 0}\left[{\frac {\beta x^{2}}{2}}-\ln I_{0}(\beta x)\right]}

と求まる（適当な規格化を施した）[5]。自由エネルギーの極値条件として方程式

x = I 1 ( x ) I 0 ( x ) {\displaystyle x={\frac {I_{1}(x)}{I_{0}(x)}}}

が非自明解を持つか、という条件が得られ、上の考察が再現される[14]。短距離相互作用する1次元系では自発的対称性の破れによる相転移が起きないことが保証されており（マーミン・ワグナーの定理（英語版））、この結果は HMF モデルの熱力学的性質において長距離相互作用が本質的であることを示している[2]。

歴史

HMF モデルの原型となるモデルは1992年に小西哲郎と金子邦彦によって導入された[15][16]。彼らは非線型力学（カオス理論）の観点からこの模型に興味を持ち、ある離散的なシンプレクティック写像の反復により粒子群がクラスターを形成し有限時間の後にクラスターが散開することを観察した[15]。1993年に稲垣省五と小西はその連続時間における対応物について運動論的方程式を用いて研究した[17]が、これが現在 HMF モデルとして知られているものである。また稲垣と小西はこのモデルにおけるクラスターの形成過程は自己重力系のジーンズ不安定性に相当する不安定性によるものであると指摘した。続く論文で稲垣は HMF モデルの熱平衡状態の熱力学的安定性を論じ[18]、1996年には運動学的方程式を用いて HMF モデルのボルツマンエントロピーが減少しないことを示した[19]。なお稲垣はこの模型を modified Konishi-Kaneko model と呼称している[18][20]。同時期に Christophe Pichon とドナルド・リンデン＝ベルは棒渦巻銀河におけるバー構造の観点から同じ模型について調べていた[3]。

1995年に Antoni & Ruffo は強磁性体の統計力学の文脈で、XY模型の平均場長距離相互作用の極限に相当するモデルについて研究し、それを Hamiltonian mean-field X-Y model と命名した[7]。またその論文の中でこの模型が稲垣、小西らによって既に調べられていたものと同じものであると指摘している[21]。Dauxios らは2002年に HMF モデルの相転移の性質についてより詳しく調べている[22]。なお結合定数が 1 / N {\displaystyle 1/N} に比例してスケールするのではない場合、熱力学極限でエネルギーの示量性が成立せず熱力学的振る舞いが異なるが、この場合の HMF モデルの熱力学については Tamarit & Anteneodo (2000), Toral (2004) らによって調べられた[13][23][24]。その他に HMF モデルの quasi-stationary state (QSS) の性質[25][26][27][28]、リンデン＝ベル統計の応用[29][30][31]、対応する量子多体系の性質[32][33]といった研究が行われている。

脚注

  1. ^ a b 後藤 & 山口 2005, p. 1.
  2. ^ a b c d e f g Levin et al. 2014, p. 36.
  3. ^ a b Chavanis, Vatteville & Bouchet 2005, p. 62.
  4. ^ Campa, Dauxois & Ruffo 2009, pp. 87–88.
  5. ^ a b c d e f Campa, Dauxois & Ruffo 2009, p. 88.
  6. ^ Levin et al. 2014, p. 37.
  7. ^ a b c Antoni & Ruffo 1995, p. 2361.
  8. ^ 西森秀稔『相転移・臨界現象の統計物理学』 35巻、培風館〈新物理学シリーズ〉、2005年、16頁。ISBN 978-4-563-02435-2。 
  9. ^ 稲垣 1993, pp. 137–138.
  10. ^ Chavanis, Vatteville & Bouchet 2005, pp. 64–65.
  11. ^ a b c d Chavanis, Vatteville & Bouchet 2005, p. 65.
  12. ^ Chavanis, Vatteville & Bouchet 2005, p. 68.
  13. ^ a b Campa, Dauxois & Ruffo 2009, p. 89.
  14. ^ Campa, Dauxois & Ruffo 2009, pp. 88–89.
  15. ^ a b Konishi, Tetsuro; Kaneko, Kunihiko (1992). “Clustered motion in symplectic coupled map systems”. Journal of Physics A: Mathematical and General 25 (23): 6283-6296. Bibcode: 1992JPhA...25.6283K. doi:10.1088/0305-4470/25/23/023. 
  16. ^ Chavanis, Vatteville & Bouchet 2005, p. 61.
  17. ^ Inagaki, Shogo; Konishi, Tetsuro (1993). “Dynamical Stability of a Simple Model Similar to Self-Gravitating Systems”. Publications of the Astronomical Society of Japan 43: 733-735. Bibcode: 1993PASJ...45..733I. 
  18. ^ a b Inagaki, S. (1993). “Thermodynamic Stability of Modified Konishi-Kaneko System”. Progress of Theoretical Physics 90 (3): 577–584. Bibcode: 1993PThPh..90..577I. doi:10.1143/ptp/90.3.577. ISSN 0033-068X. 
  19. ^ Inagaki, S. (1996). “Kinetic Equation for the Modified Konishi-Kaneko System”. Progress of Theoretical Physics 96 (6): 1307–1309. Bibcode: 1996PThPh..96.1307I. doi:10.1143/PTP.96.1307. ISSN 0033-068X. 
  20. ^ 稲垣 1993.
  21. ^ Antoni & Ruffo 1995, p. 2362.
  22. ^ Dauxois, Thierry; Latora, Vito; Rapisarda, Andrea; Ruffo, Stefano; Torcini, Alessandro (2002). The Hamiltonian Mean Field Model: From Dynamics to Statistical Mechanics and Back. 602. pp. 458–487. arXiv:cond-mat/0208456. Bibcode: 2002LNP...602..458D. doi:10.1007/3-540-45835-2_16. ISSN 0075-8450. 
  23. ^ Tamarit, Francisco; Anteneodo, Celia (2000). “Rotators with Long-Range Interactions: Connection with the Mean-Field Approximation”. Physical Review Letters 84 (2): 208–211. arXiv:cond-mat/9911030. Bibcode: 2000PhRvL..84..208T. doi:10.1103/PhysRevLett.84.208. ISSN 0031-9007. 
  24. ^ Toral, Raúl (2004). “On the Nonextensivity of the Long Range X-Y Model”. Journal of Statistical Physics 114 (5/6): 1393–1398. arXiv:cond-mat/0304018. Bibcode: 2004JSP...114.1393T. doi:10.1023/B:JOSS.0000013963.16180.a3. ISSN 0022-4715. 
  25. ^ Levin et al. 2014, pp. 38–43.
  26. ^ Morita, Hidetoshi; Kaneko, Kunihiko (2006). “Collective Oscillation in a Hamiltonian System”. Physical Review Letters 96 (5). Bibcode: 2006PhRvL..96e0602M. doi:10.1103/PhysRevLett.96.050602. ISSN 0031-9007. 
  27. ^ Pakter, Renato; Levin, Yan (2013). “Topology of Collisionless Relaxation”. Physical Review Letters 110 (14). arXiv:1304.5252. Bibcode: 2013PhRvL.110n0601P. doi:10.1103/PhysRevLett.110.140601. ISSN 0031-9007. 
  28. ^ Campa, Alessandro; Chavanis, Pierre-Henri (2013). “Caloric curves fitted by polytropic distributions in the HMF model”. The European Physical Journal B 86 (4). arXiv:1210.4082. Bibcode: 2013EPJB...86..170C. doi:10.1140/epjb/e2013-30947-0. ISSN 1434-6028. 
  29. ^ Levin et al. 2014, pp. 43–44.
  30. ^ Chavanis, P. H. (2006). “Lynden-Bell and Tsallis distributions for the HMF model”. The European Physical Journal B 53 (4): 487–501. arXiv:cond-mat/0604234. Bibcode: 2006EPJB...53..487C. doi:10.1140/epjb/e2006-00405-5. ISSN 1434-6028. 
  31. ^ Assllani, Mallbor; Fanelli, Duccio; Turchi, Alessio; Carletti, Timoteo; Leoncini, Xavier (2012). “Statistical theory of quasistationary states beyond the single water-bag case study”. Physical Review E 85 (2). arXiv:1109.5934. Bibcode: 2012PhRvE..85b1148A. doi:10.1103/PhysRevE.85.021148. ISSN 1539-3755. 
  32. ^ Chavanis, Pierre-Henri (2011). “The quantum HMF model: I. Fermions”. Journal of Statistical Mechanics: Theory and Experiment 2011 (08): P08002. Bibcode: 2011JSMTE..08..002C. doi:10.1088/1742-5468/2011/08/P08002. ISSN 1742-5468. 
  33. ^ Chavanis, Pierre-Henri (2011). “The quantum HMF model: II. Bosons”. Journal of Statistical Mechanics: Theory and Experiment 2011 (08): P08003. Bibcode: 2011JSMTE..08..003C. doi:10.1088/1742-5468/2011/08/P08003. ISSN 1742-5468. 

参考文献

  • 稲垣省五「Stability of Modified Konishi-Kaneko System〔邦文〕」『物性研究』第61巻第2号、物性研究刊行会、1993年11月、138-141頁、hdl:2433/95196、ISSN 0525-2997、NDLJP:10937334、2024年7月18日閲覧。 
  • Antoni, Mickael; Ruffo, Stefano (1995). “Clustering and relaxation in Hamiltonian long-range dynamics”. Physical Review E 52 (3): 2361–2374. Bibcode: 1995PhRvE..52.2361A. doi:10.1103/PhysRevE.52.2361. ISSN 1063-651X. 
  • Chavanis, Pierre-Henri; Vatteville, J.; Bouchet, Freddy (2005). “Dynamics and thermodynamics of a simple model similar to self-gravitating systems: the HMF model”. European Physical Journal B: Condensed Matter and Complex Systems (Springer-Verlag) 46 (1): 61-99. arXiv:cond-mat/0408117. Bibcode: 2005EPJB...46...61C. doi:10.1140/epjb/e2005-00234-0. 
  • 後藤振一郎, 山口義幸 (2005年). “<お探しのページは見つかりません> 長距離力相互作用を有する非線形格子ハミルトン系におけるマクロ変数のカノニカル期待値への漸近的振舞”. 2021年4月22日閲覧。[リンク切れ]
  • Campa, Alessandro; Dauxois, Thierry; Ruffo, Stefano (2009). “Statistical mechanics and dynamics of solvable models with long-range interactions”. Physics Reports 480 (3-6): 57–159. arXiv:0907.0323. Bibcode: 2009PhR...480...57C. doi:10.1016/j.physrep.2009.07.001. ISSN 03701573. 
  • Levin, Yan; Pakter, Renato; Rizzato, Felipe B.; Teles, Tarcísio N.; Benetti, Fernanda P.C. (2014). “Nonequilibrium statistical mechanics of systems with long-range interactions”. Physics Reports 535 (1): 1–60. arXiv:1310.1078. Bibcode: 2014PhR...535....1L. doi:10.1016/j.physrep.2013.10.001. ISSN 03701573.