[논문 리뷰] Dynamic critical behavior of cluster algorithms for 2D Ashkin-Teller and Potts models

Jesús Salas|arXiv (Cornell University)|2000. 04. 19.

Theoretical and Computational Physics참고 문헌 27인용 수 3

한 줄 요약

이 논문은 2차원 아쉬킨–틀러 모델에 대한 스웨ン센–워우와 관련된 클러스터 알고리즘의 동적 임계 행동을 조사하며, 자가상관 시간에 대한 리–소칼 하한 경계가 정확하지는 않지만 거의 정확하다는 것을 발견한다. 이 비율인 통합 자가상관 시간과 비열의 비율은 체계 크기의 로그 또는 작은 거듭제곱(0.05–0.12)으로 증가한다. 에너지에 대한 지수 및 통합 자가상관 시간이 동일한 동적 임계 지수로 스케일링됨을 보여, 거의 최적의 동적 스케일링 행동임을 시사한다.

ABSTRACT

We study the dynamic critical behavior of two algorithms: the Swendsen-Wang algorithm for the two-dimensional Potts model with q=2,3,4 and a Swendsen-Wang-type algorithm for the two-dimensional symmetric Ashkin-Teller model on the self-dual curve. We find that the Li--Sokal bound on the autocorrelation time τ_{{ m int},{\cal E}} \geq const imes C_H is almost, but not quite sharp. The ratio τ_{{ m int},{\cal E}}/C_H appears to tend to infinity either as a logarithm or as a small power (0.05 \ltapprox p \ltapprox 0.12). We also show that the exponential autocorrelation time τ_{{ m exp},{\cal E}} is proportional to the integrated autocorrelation time τ_{{ m int},{\cal E}}.

연구 동기 및 목표

• 스웨닝슨–워우 알고리즘이 일부 모델에서는 잘 작동하지만 다른 모델에서는 그렇지 않은 이유를 이해하기 위해.
• 자기상관 시간에 대한 리–소칼 하한 경계의 날카움을 테스트하기 위해, 이는 동적 임계 지수와 정적 임계 지수를 연결한다.
• 2차원 클러스터 알고리즘에서 에너지에 대한 지수 및 통합 자가상관 시간이 동일한 동적 임계 지수로 스케일링되는지 확인하기 위해.
• 비국소 몬테카를로 알고리즘의 역학에서 느린 모드의 성격을 조사하기 위해.

제안 방법

• 스웨닝슨–워우 알고리즘과 관련된 알고리즘을 사용하여 2차원 q-상 페츠 모델(q=2,3,4)과 자기 dual 곡선 상의 대칭 아쉬킨–틀러 모델의 수치적 몬테카를로 시뮬레이션을 수행한다.
• 긴 시뮬레이션에서 얻은 고정밀도 데이터를 이용해 에너지 관측치의 통합 및 지수 자가상관 시간을 추정한다.
• 체계 크기 L에 따른 자가상관 시간의 스케일링을 분석하여 유한 체계 스케일링 기법을 적용해 동적 임계 지수를 추출한다.
• 자기상관 함수가 순수 지수 함수로 작동하는지 테스트하기 위해 스케일링 가정 ρ_AA(t;L) ≈ |t|^{-p_A} h_A[t/τ_exp,A; ξ(L)/L] 을 적용한다.
• 로그 및 거듭제곱 형태의 피팅을 사용하여, τ_int,ℰ / C_H 비율과 리–소칼 경계의 날카움을 평가한다.
• 다양한 격자 크기에서 ρ_ℰℰ(t) vs. t/τ_int,ℰ 의 데이터 콜라프스를 확인함으로써, τ_int,ℰ ≈ τ_exp,ℰ 의 가정을 검증한다.

실험 결과

연구 질문

• RQ1 2차원 시스템에서 클러스터 알고리즘에 대한 리–소칼 자가상관 시간 하한 경계는 날카운가?
• RQ2 스웨닝슨–워우 알고리즘에서 에너지에 대한 지수 및 통합 자가상관 시간이 동일한 동적 임계 지수로 스케일링되는가?
• RQ3 비율 τ_int,ℰ / C_H 는 체계 크기 L에 따라 어떻게 스케일링되며, 로그 또는 작은 거듭제곱으로 증가하는가?
• RQ4 2차원 페츠 및 아쉬킨–틀러 모델의 역학에서 가장 느린 모드의 성격은 무엇인가?
• RQ5 왜 리–소칼 경계는 2차원에서는 거의 정확하지만 고차원에서는 명백히 정확하지 않은가?

주요 결과

• 리–소칼 경계 τ_int,ℰ ≥ const × C_H 는 정확하지 않지만, 비율 τ_int,ℰ / C_H 는 체계 크기 L에 따라 로그 또는 작은 거듭제곱(0.05 ≤ p ≤ 0.12)으로 증가한다.
• 2차원 페츠 모델(q=2,3,4)과 아쉬킨–틀러 모델에서 에너지의 동적 임계 지수 z_int,ℰ 는 z_int,ℰ ≈ z_exp,ℰ 와 일치하여, 지수 및 통합 자가상관 시간이 동일한 지수로 스케일링됨을 나타낸다.
• 2차원 이징 모델(q=2)의 경우 비율 τ_int,ℰ / C_H 는 로그로 증가하여 알고리즘의 동적 스케일링이 거의 최적임을 시사한다.
• 스케일링 τ_int,ℰ / C_H ∼ L^p 에서의 거듭제곱 p 는 이징 모델에서 약 0.05에서 4상 페츠 모델에서 0.12로 증가하여, 다양한 모델 간의 연속적인 추세를 보인다.
• L ≥ 64 에서 ρ_ℰℰ(t) vs. t/τ_int,ℰ 의 데이터 콜라프스는 뛰어나며, 에너지 자가상관 함수가 순수 지수 함수로 작동한다는 가정을 지지한다.
• 이 연구는 d=2 에서는 에너지와 결합 점유율이 가장 느린 모드에 속하지만, d>2 에서는 그렇지 않다는 점을 시사하며, 리–소칼 경계의 날카움 차이를 설명한다.