[論文レビュー] Entropy of fully-packed rigid rods on generalized Husimi trees: a route to the square lattice limit

Nathann T. Rodrigues, Jürgen F. Stilck+1 more|arXiv (Cornell University)|Nov 30, 2021

Theoretical and Computational Physics参考文献 45被引用数 6

ひとこと要約

本稿では、サイズLが増加する正方形格子クラスタから構築された一般化されたフーシミ木を用いて、完全に詰め詰めの剛体k-merが正方形格子上に配置された系のサイト当たりエントロピーを体系的に推定する手法を提案する。Lの偶奇性に起因するL-パリティ効果とL→∞への外挿を用いることで、極めて高い精度の推定値が得られる：s(2)は正確解と0.03%の差異にとどまり、s(3)は3%の差異、s(4)については新たな推定値が得られ、k≥4では無秩序相である可能性が示唆される。

ABSTRACT

Although hard rigid rods ($k$-mers) defined on the square lattice have been widely studied in the literature, their entropy per site, $s(k)$, in the full-packing limit is only known exactly for dimers ($k=2$) and numerically for trimers ($k=3$). Here, we investigate this entropy for rods with $k \le 7$, by defining and solving them on Husimi lattices built with diagonal and regular square lattice clusters of effective lateral size $L$, where $L$ defines the level of approximation to the square lattice. Due to an $L$-parity effect, by increasing $L$ we obtain two systematic sequences of values for the entropies $s_L(k)$ for each type of cluster, whose extrapolations to $L ightarrow \infty$ provide estimates of these entropies for the square lattice. For dimers, our estimates for $s(2)$ differ from the exact result by only $0.03\%$, while that for $s(3)$ differs from best available estimates by $3\%$. In this paper, we also obtain a new estimate for $s(4)$. For larger $k$, we find that the extrapolated results from the Husimi tree calculations do not lie between the lower and upper bounds established in the literature for $s(k)$. In fact, we observe that, to obtain reliable estimates for these entropies, we should deal with levels $L$ that increase with $k$. However, it is very challenging computationally to advance to solve the problem for large values of $L$ and for large rods. In addition, the exact calculations on the generalized Husimi trees provide strong evidence for the fully packed phase to be disordered for $k\geq 4$, in contrast to the results for the Bethe lattice wherein it is nematic, thus providing evidence for a high density nematic-disordered transition in the system of $k$-mers with vacancies.

研究の動機と目的

• 完全に詰め詰めの剛体k-merが正方形格子上に配置された系のサイト当たりエントロピーs(k)を計算すること。この問題はk=2の場合にのみ正確に解かれており、k=3については数値的にのみ解かれた。
• kが大きい場合の正方形格子上での直接シミュレーションの計算不能性を克服するため、一般化されたフーシミ木を用いた階層的近似を用いること。
• k≥4の場合の完全に詰め詰め相の性質を調査すること—特に、それがネマチック相であるか無秩序相であるかを明らかにすること。
• フーシミ木の構築において有効なクラスターサイズLを増大させることで、正方形格子極限への体系的かつ段階的な到達ルートを確立すること。
• 特に理論的境界と計算スケーリングの課題を踏まえて、kが大きい場合におけるこの手法の信頼性を評価すること。

提案手法

• 有効な横方向サイズLを持つ正方形格子のクラスタから一般化されたフーシミ木を構築し、正方形格子の階層的近似を形成する。
• 有限Lにおける自由エネルギーとエントロピーを正確に計算するために、これらの木上にモデルを定義し、正確な転送行列技術を用いる。
• L-パリティ効果を活用する：偶数Lと奇数Lのそれぞれに対して、二つの異なるエントロピー推定値の系列sL(k)が得られ、L→∞への外挿が可能になる。
• L→∞極限における外挿値を、正方形格子上でのs(k)の推定値とし、Lが増大するにつれて精度が向上する。
• k≤7のk-merに対してこの手法を適用し、既知の境界と数値的推定値とを比較する。
• 特にk≥4の場合に注目して、フーシミ木上の基底状態の性質を分析することで、相の性質を調査する。

実験結果

リサーチクエスチョン

• RQ1 k>3の場合に、体系的な近似手法を用いて、完全に詰め詰めのk-merが正方形格子上に配置された系のサイト当たりエントロピーs(k)を高精度に推定できるか？
• RQ2 一般化されたフーシミ木におけるL-パリティ効果は、k≥4の場合のエントロピー推定値の収束性と信頼性にどのように影響するか？
• RQ3 k≥4の場合の完全に詰め詰め相の性質—ネマチック相か無秩序相か—は、フーシミ木上での正確解に基づいてどのように特定できるか？
• RQ4 外挿されたエントロピー値は、既存の理論的境界や数値的推定値とどのように一致するか？
• RQ5 kが大きくなることやLが高くなることに対応する際、この手法の計算的および概念的限界は何か？

主な発見

• 本手法によりs(2) = 0.29156が得られ、正確解と0.03%の差異にとどまり、二量体に対して高い精度が確認された。
• 三量体の場合、外挿されたs(3)は、利用可能な最高の数値的推定値と3%の差異にとどまり、強い一貫性が示された。
• s(4)について新たな推定値が提供され、s(4) ≈ 0.0876となったが、論文ではこの値が既知の理論的境界の外にあると指摘している。
• k≥4の場合、外挿されたエントロピー値は、文献に既存する下限と上限の間には位置しなかったため、より高いL値が必要であると示唆された。
• フーシミ木上での正確解は、k≥4の場合に完全に詰め詰め相が無秩序相であるという強い証拠を提供しており、ベーテ格子上ではネマチック相が得られるのとは対照的である。
• 結果は、特にkが大きい場合に、空孔を伴うk-mer系において高密度ネマチック相から無秩序相への転移が存在することを支持する。