[논문 리뷰] 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명|arXiv (Cornell University)|2021. 11. 30.

Theoretical and Computational Physics참고 문헌 45인용 수 6

한 줄 요약

이 논문은 점점 크기가 커지는 L 크기의 정사각형 격자 클러스터에서 구성된 일반화된 투스미 나무를 이용해, 정사각형 격자 위의 완전히 밀도가 높은 강성 있는 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일 경우에만 수치적으로 해결됨.
• 더 큰 k에 대해 정사각형 격자에서 직접 시뮬레이션을 수행하는 데 있어 계산적으로 불가능한 문제를 해결하기 위해, 일반화된 투스미 나무를 통한 계층적 근사 방법을 사용하는 것.
• k≥4일 경우의 완전히 밀도가 높은 상의 성질 — 특히 이 상이 정렬된(nematic) 상태인지 무질서한(disordered) 상태인지 여부를 조사하는 것.
• 투스미 나무 구성에서 효과적인 클러스터 크기 L을 증가시킴으로써 정사각형 격자 극한에 대한 체계적인 접근로드를 확립하는 것.
• 특히 계산 스케일링의 과제와 알려진 이론적 한계를 감안할 때, 더 큰 k에 대해 이 방법의 신뢰성을 평가하는 것.

제안 방법

• 정사각형 격자 클러스터에서 효과적인 횡방향 크기 L을 갖는 일반화된 투스미 나무를 구성하여, 정사각형 격자에 대한 계층적 근사로 사용한다.
• 유한한 L에 대해 정확한 전이 행렬 기법을 사용해 이 나무 위에 모델을 정의하고, 자유 에너지와 엔트로피를 계산한다.
• L의 기수 효과를 활용: 짝수와 홀수 L에 대해 두 개의 서로 다른 엔트로피 추정치 시퀀스 sL(k)가 도출되며, 이를 통해 L→∞로의 외삽이 가능해진다.
• L→∞ 극한에서의 외삽된 값들을 정사각형 격자에서의 s(k)에 대한 추정치로 사용하여, L이 증가할수록 정확도가 향상된다.
• k≤7에 대해 이 방법을 적용하고, 결과를 알려진 이론적 한계와 수치적 추정치와 비교한다.
• 특히 k≥4일 경우에 대해, 투스미 나무의 기저 상태 성질을 분석함으로써 상의 행동을 분석한다.

실험 결과

연구 질문

• RQ1 k>3에 대해 체계적인 근사 방법을 사용해 정사각형 격자 위 완전히 밀도가 높은 k-머스의 단위 면적당 엔트로피 s(k)를 정확하게 추정할 수 있는가?
• RQ2 일반화된 투스미 나무에서 관찰되는 L 기수 효과가 k≥4에 대해 엔트로피 추정치의 수렴성과 신뢰성에 어떤 영향을 미치는가?
• RQ3 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이 필요로 함을 시사한다.
• 투스미 나무에서의 정확한 해는 k≥4에 대해 완전히 밀도가 높은 상이 무질서한 상태임을 강력히 뒷받침하며, 베티 격자에서 관찰된 정렬된 상과는 대조된다.
• 결과는 특히 더 큰 k에 대해 흐린 빈자리가 존재하는 k-머스 시스템에서 고밀도 정렬상에서 무질서상으로의 전이가 존재함을 지지한다.