[논문 리뷰] Exploring Lee-Yang and Fisher Zeros in the 2D Ising model through multipoint Padé approximants

Simran Singh, Massimo Cipressi외 1명|arXiv (Cornell University)|2023. 12. 05.

Theoretical and Computational Physics인용 수 4

한 줄 요약

이 논문은 유한체적 몬테카를로 시뮬레이션을 사용하여 2차원 이징 모형에서 리-양 및 파이저 영역을 다항식 패드에 응용하는 새로운 방법을 제시한다. 여러 자기장 및 역온도 값에서 분할함수의 누적량을 분석함으로써 복소 평면상의 영역과 극을 정확히 특정하여 임계 지수를 정밀하게 추정한다: $\beta\delta = 1.879(16)$ 및 $\nu = 1.014(60)$, $\beta_c = 0.4405(5)$로, 이는 방법의 위상전이 연구에 대한 강건성을 확인한다.

ABSTRACT

We present a numerical calculation of the Lee-Yang and Fisher zeros of the 2D Ising model using multi-point Padé approximants. We perform simulations for the 2D Ising model with ferromagnetic couplings both in the absence and in the presence of a magnetic field using a cluster spin-flip algorithm. We show that it is possible to extract genuine signature of Lee Yang and Fisher zeros of the theory through the poles of magnetization and specific heat, using multi-point Padé method. We extract the poles of magnetization using Padé approximants and compare their scaling with known results. We verify the circle theorem associated to the well known behaviour of Lee Yang zeros. We present our finite volume scaling analysis of the zeros done at $T=T_c$ for a few lattice sizes, extracting to a good precision the (combination of) critical exponents $βδ$. The computation at the critical temperature is performed after the latter has been determined via the study of Fisher zeros, thus extracting both $β_c$ and the critical exponent $ν$. Results already exist for extracting the critical exponents for the Ising model in 2 and 3 dimensions making use of Fisher and Lee Yang zeros. In this work, multi-point Padé is shown to be competitive with this respect and thus a powerful tool to study phase transitions.

연구 동기 및 목표

• 유한체적 시스템에서 외부 자기장 $H$와 역온도 $\beta$의 다중값에서 분할함수의 주요 누적량을 재수열하는 다중점 패드 근사법을 사용하여 리-양 및 파이저 영역을 추출하는 방법을 개발하고 검증한다.
• 몬테카를로 시뮬레이션에서 발생하는 통계적 노이즈가 존재하는 상황에서 다중점 패드 방법의 강건성을 시험한다.
• 2차원 이징 모형에서 유한체적 스케일링을 통해 영역의 위치를 분석하여 임계 지수 $\beta\delta$, $\nu$, $\beta_c$를 추출한다.
• 리-양 영역에 대한 원의 정리가 성립하는지 확인하고, 파이저 영역 분석을 통해 $T_c$가 자가일관적으로 결정되는지 검증한다.
• 다중점 패드 방법을 랜드스케일 모델에서 임계 현상 분석을 위한 전통적 유한체적 스케일링 대비 경쟁 가능한 대안으로 정립한다.

제안 방법

• 외부 자기장 $H$와 역온도 $\beta$의 여러 값에서 평가된 분할함수의 주요 누적량을 재수열하는 데 다중점 패드 근사법을 사용한다.
• 자기화는 다중점 패드를 통해 $H$의 유리함수로 근사하고, 그 극을 분석하여 복소 $H$ 평면상의 리-양 영역을 특정한다.
• 비슷하게 열용량은 $\beta$의 유리함수로 근사하여 복소 $\beta$ 평면상의 파이저 영역을 특정한다.
• 유한체적 스케일링을 리-양 및 파이저 영역 중 가장 가까운 영역의 위치에 적용하여 임계 지수 $\beta\delta$, $\nu$, $\beta_c$를 추출한다.
• 다양한 크기 $L$의 정사각형 격자에서 $T_c$에서 클러스터 스핀-반전 알고리즘을 사용하여 시뮬레이션을 수행하며, 각 $H$ 또는 $\beta$에 대해 약 650,000개의 구성표를 확보하여 정확도를 확보한다.
• 오차 분석은 $\chi^2/dof$와 신뢰구간을 사용하여 임계 지수 스케일링 법칙에 대한 피팅의 신뢰성을 평가한다.

실험 결과

연구 질문

• RQ1 다중점 패드 근사법은 통계적 노이즈가 있는 유한체적 몬테카를로 데이터로부터 리-양 및 파이저 영역을 신뢰성 있게 추출할 수 있는가?
• RQ2 자기화 및 열용량 근사식의 극은 2차원 이징 모형에서 알려진 리-양 영역에 대한 원의 정리를 재현하는가?
• RQ3 가장 가까운 리-양 및 파이저 영역의 유한체적 스케일링을 통해 임계 지수 $\beta\delta$, $\nu$, $\beta_c$의 정밀한 추정치를 도출할 수 있는가?
• RQ4 파이저 영역 분석을 통해 임계 온도 $T_c$가 자가일관적으로 결정되며, 기존 값과 일치하는가?
• RQ5 임계 현상 분석에서 전통적 유한체적 스케일링 대비 다중점 패드 방법의 정밀도와 강건성은 어떻게 비교되는가?

주요 결과

• 다중점 패드 방법은 자기화 근사식의 극으로 진정한 리-양 영역을 성공적으로 식별하며, 모든 안정된 극이 허수 $H$축 위에 위치하여 원의 정리를 확인한다.
• 가장 가까운 리-양 영역의 허수부에 대한 유한체적 스케일링은 $\beta\delta = 1.879(16)$를 도출하며, 정확한 값 1.875와 뛰어난 일치를 보인다.
• 파이저 영역의 유한체적 스케일링은 $\nu = 1.014(60)$ 및 $\beta_c = 0.4405(5)$를 도출하며, 정확한 결과와 일치하고 $T_c$의 자가일관적 결정을 검증한다.
• 통계적 노이즈가 존재하더라도 강건성을 보이며, 유한 통계 조건에서도 고정밀 결과를 확보함으로써, QCD와 같은 더 복잡한 시스템에도 적용 가능함을 시사한다.
• 임계 지수 $\beta\delta$는 상대 오차 약 $0.8\%$로 추정되어, 유한체적 스케일링 분석에서 높은 정밀도를 입증한다.
• 결과는 다중점 패드 방법이 랜드스케일 시뮬레이션에서 임계 지수를 추출하는 데 있어 전통적 유한체적 스케일링 대비 경쟁력 있고 신뢰할 수 있는 대안임을 확인한다.