The freezing phase transition in hard core lattice gases on triangular lattice with exclusion up to seventh next-nearest neighbor

技術要約：三角格子上のハードコア格子ガスにおける凍結相転移

問題提起
ハードコア格子ガス（HCLG）モデルは、密度増加に伴うエントロピー駆動型の相転移、特に流体から固体への凍結転移を研究するための最小限の枠組みとして機能する。$k$-NNモデルでは、粒子は自身から$k$番目の次近接隣接サイトまでの占有を排除する。正方格子およびハニカム格子における$k$-NNモデルの相図は比較的よく確立されているが、三角格子における$k \ge 4$の場合については、依然として理解が進んでいない。三角格子に関する先行研究は、$k \le 3$に限定されていた。$k=4$および$k=5$については、テンソル繰り込み群（TRG）法を用いた結果、転移の性質（連続か不連続か）に関して決定的な結論が出ず、臨界パラメータについても広い範囲の値が示された。また、大きな排除体積を持つHCLGモデルの標準的なモンテカルロ・シミュレーションは、局所的な移動のみを許可する場合、充填率が最大に近づくとメタステーブルな状態（準安定状態）にトラップされるため、平衡化の問題に直面することが多い。

手法
著者らは、ストリップ・クラスター・ワン・ランダウ（SCWL）アルゴリズムを用いて、三角格子上の$1 \le k \le 1 \le 7$の$k$-NNモデルを調査している。この手法は、フラットヒストグラム・アプローチと棄却のないクラスター移動を組み合わせることで、密度状態 $g(N, L)$（ここで $N$ は粒子数、$L$ は線形システムサイズ）を効率的にサンプリングする。

主要なアルゴルリズムの特徴は以下の通りである：

• クラスター移動: アルゴリズムは、格子を横断するストリップ上で構成を更新する。特定の行（$k$に依存する方向）を選択し、排除サイト間の開いたセグメントを特定する。セグメント内の粒子は、現在の密度状態の推定値から導出された確率に従って、除去および再配置される。
• フラットヒストグラム・プロトコル: シミュレーションでは、訪問された粒子数のヒストグラムを平坦化するためにワン・ランダウ・プロトコルを採用する。修正係数 $f$ は、ヒストグラムが十分に平坦（$H_{min} > 0.8 H_{max}$）になるまで、反復的に減少させる（$f=1$ から $f_{min} = 2^{-22}$ まで）。
• 熱力学的解析: $g(N, L)$ が決定された後、グランドカノニカル分配関数、圧力、および平均観測量が計算される。相転移の性質は以下の手法を用いて分析される：
  • 分配関数ゼロ（PFZ）: 複素化学ポテンシャル平面（$z = e^\mu$）におけるゼロの軌跡を調べる。
  • エントロピーの非凸性（NC）: エントロピー $S(N) = \ln g(N, L)$ を解析し、一次転移を示す非凸領域を特定する。
  • 有限サイズスケーリング: 二次元における一次転移に適したスケーリング則（$\nu = 1/2$）を用いて、臨界指数およびパラメータを熱力学的極限（$L \to \infty$）へと外挿する。

主な貢献と結果

1. ベンチマークと検証:

   • $k=1$（ハードヘキサゴンモデル）において、SCWLによる臨界化学ポテンシャル（$\mu_c \approx 2.4064$）および臨界密度（$\rho_c \approx 0.826$）の結果は厳密解とよく一致しており、アルゴリズムの正確性を検証している。
   • $k=2$ については、臨界密度（$\rho_c = 0.7419(5)$）の改良された推定値を提供しており、これは転送行列の結果と密接に一致し、一部の先行モンテカルロ推定とは異なる。

2. $4 \le k \le 7$ における新知見:

   • 単一相転移: 大きな排除体積に対する予想に反して、検討されたすべての $k$ 値（$4, 5, 6, 7$）において、系は低密度の無秩序流体相から高密度の副格子秩序相への「単一の」相転移を起こす。
   • 不連続な性質: 転移は $k \ge 4$ で不連続（一次）であることが確立された。これは以下の証拠による：
     • PFZが円を形成し、角度 $\pi/2$ で実軸を一点で挟み込むこと。
     • 数値的に得られたエントロピーに非凸領域が存在すること。
   • 臨界パラメータ: 著者らは、$k=4$ から $7$までの臨界化学ポテンシャル（$\mu_c$）、共存密度（$\rho_f$および$\rho_s$）、および臨界圧力（$P_c$）の正確な推定値を提供している。
     • 例 ($k=4$): $\mu_c = 2.8696(2)$, $\rho_f = 0.7404(2)$, $\rho_s = 0.9067(2)$。
     • 例 ($k=5$): $\mu_c = 4.720(3)$, $\rho_f = 0.916(3)$, $\rho_s = 0.988(3)$。

3. 副格子秩序:

   • 高密度相は、特定の副格子の占有に対応する。副格子の数（$N_{sub}$）は $k$ とともに増加する：$k=4$ で $N_{sub}=9$、$k=5$ で $12$、$k=6$で$13$、$k=7$で$16$ である。

意義とTRGとの比較
本論文は、$k=4$ および $k=5$ におけるSCWLの結果と、従来のTRG研究との間の重大な相違を強調している。

• TRGの限界: 従来のTRG研究は $k=4$ に対して連続転移を示唆しており、提供された臨界化学ポテンシャルの範囲（$2.65 < \mu_c < 2.7$）は、SCWLの結果（$\mu_c \approx 2.87$）と重ならない。$k=5$ について、TRGは $4.4 < \mu_c < 4.5$ という範囲を示唆したが、SCWLは $\mu_c \approx 4.72$ を見出している。
• SCWLの優位性: 著者らは、高密度状態への効率的なアクセスと、局所移動アルゴリズムやTRGにおいてこれらの特定の領域で問題となる平衡化の問題を克服できる能力により、SCWLアルゴリズムが優れた精度を持つことを挙げている。

結論
本研究は、三角格子上の $k$-NN モデルにおいて、$4 \le k \le 7$ の場合、系は単一の不連続な相転移を示すことを確立した。SCWLアルゴリズムは、大きな排除体積を持つハードコア格子ガスの相図および臨界パラメータを決定するための堅牢なツールであり、この領域においてTRG法よりも精密な結果を提供する。著者らは、連続体極限（ハードスフェアガス）では、非常に大きな $k$ に対して中間相（ヘキサティック相など）が現れることが予想されるものの、今回研究した計算可能な範囲の $k \le 7$ ではそのような相は観察されないと述べている。