[论文解读] 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的正方晶格簇构建广义Husimi树并求解该模型，估算正方晶格上完全填充刚性k-簇的每格点熵。利用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-簇的每格点熵s(k)，该问题仅对k=2有精确解，对k=3仅有数值解。
• 通过广义Husimi树的分层近似方法，克服在正方晶格上对较大k值直接模拟的计算不可行性。
• 研究完全填充相的性质——特别是当k≥4时是否为取向有序（向列相）或无序相。
• 通过在Husimi树构建中逐步增大有效簇尺寸L，建立通往正方晶格极限的系统性路径。
• 评估该方法在较大k值下的可靠性，尤其结合已知理论边界与计算规模增长的挑战。

提出的方法

• 从具有有效横向尺寸L的正方晶格簇构建广义Husimi树，形成对正方晶格的分层近似。
• 在这些树上使用精确转移矩阵技术定义模型，计算有限L下的自由能与熵。
• 利用L的奇偶性效应：对偶数L和奇数L分别产生两组不同的熵估计序列sL(k)，从而实现外推至L→∞。
• 将L→∞极限下的外推值用作正方晶格上s(k)的估计值，且随着L增大，精度不断提高。
• 将该方法应用于k≤7的k-簇，与已知边界和数值估计进行比较。
• 通过分析Husimi树上基态的性质，特别是对k≥4的情况，研究相行为。

实验结果

研究问题

• RQ1 能否使用系统性近似方法准确估算正方晶格上完全填充k-簇的每格点熵s(k)，其中k>3？
• RQ2 广义Husimi树中的L奇偶性效应如何影响k≥4时熵估计的收敛性与可靠性？
• RQ3 基于Husimi树上的精确解，k≥4时完全填充相的性质是向列相还是无序相？
• RQ4 外推得到的熵值与现有理论边界及s(k)的数值估计相比如何？
• RQ5 当应用于更大k值和更高L时，该方法在计算与概念层面存在哪些局限性？

主要发现

• 该方法得到s(2) = 0.29156，与精确结果仅相差0.03%，证实对二聚体具有极高精度。
• 对三聚体而言，外推得到的s(3)与目前最佳数值估计相差3%，表明结果高度一致。
• 为s(4)提供了新估计值，s(4) ≈ 0.0876，但论文指出该值位于已知理论边界的外部。
• 对于k≥4，外推得到的熵值未落在文献中已建立的下限与上限之间，表明需进一步增大L才能获得可靠估计。
• Husimi树上的精确解强烈表明，当k≥4时完全填充相为无序相，与Bethe晶格上发现的向列相形成对比。
• 结果支持在存在空位的k-簇系统中存在高密度向列相到无序相的转变，尤其在k较大时。