Crossover in the Ordered Phase in the Non-Mermin-Wagner-Hohenberg Regime of Spin Models with Long-Range Coupling

기술 요약: 장거리 결합을 가진 스핀 모델의 비-메르민-와그너-호헨베르크(Non-Mermin-Wagner-Hohenberg) 영역 내 질서상에서의 크로스오버

문제 제기
$r^{-\sigma}$로 감쇠하는 장거리 상호작용을 갖는 고전적 스핀 모델은 단거리 상호작용 모델과 확연히 다른 복잡한 상 거동을 보인다. 메르민-와그너-호헨베르크(MWH) 정리는 저차원 시스템에서 단거리 상호작용이 존재할 때 진정한 장거리 질서(LRO)를 금지하지만, 장거리 상호작용의 존재는 "비-MWH 영역"($D < \sigma < 2D$, 여기서 $D$는 공간 차원)에서 진정한 LRO를 허용한다. 기존 연구들은 주로 질서상과 무질서상 사이의 전이 또는 특정 $\sigma$ 값에 집중해 왔다. 그러나 진정한 장거리 질서상 내부에서의 스케일링 거동에 대한 포괄적인 조사는 부족한 실정이다. 특히, 이 영역에서 자화(magnetization)가 유한한 값으로 단조 수렴하는지, 아니면 시스템 크기와 상호작용 범위에 의존하는 구별된 스케일링 영역을 보이는지는 여전히 불분명하다.

방법론
저자들은 1차원(1D) 및 2차원(2D) XY 모델과 하이젠베르크(Heisenberg) 모델을 장거리 결합 $J/|i-j|^\sigma$를 사용하여 연구하기 위해 단일 업데이트 메트로폴리스 알고리즘을 이용한 몬테카를로 시뮬레이션을 채택하였다.

• 모델: 시뮬레이션은 주기적 경계 조건을 가진 1D 체인과 2D 정사각형 격자에서 수행되었다. 해밀토니안은 $r^{-\sigma}$로 감쇠하는 상호작용을 포함한다.
• 시뮬레이션 파라미터: 시스템 크기($L$)는 1D의 경우 최대 $L=4000$, 2D의 경우 $L=120$까지 사용되었다. $10^5$회의 웜업 스윕(warm-up sweeps) 후 $5 \times 10^6$회의 샘플링 스윕을 수행하였다. 낮은 온도($T < 0.1$)에서의 메타스테빌리티(metastability)를 완화하기 위해 제안된 스핀 플립을 제한하는 특별한 조치를 취했다.
• 관측량: 연구에서는 에너지($E$), 자화($M$), 자기 감수율($\chi$), 그리고 자기 큐넌트(magnetic cumulants, Binder ratio, $U_n$)를 분석하였다. 결정적으로, 저자들은 전통적인 교차점(intersection) 방식보다 더 정밀하게 상 전이를 찾기 위해 역온도($\beta$)에 대한 큐넌트와 자화의 1차 도함수를 활용하였다.
• 스케일링 분석: 임계 지수($\nu, \gamma$)와 임계 온도($T_c$)를 추출하기 위해 유한 크기 스케일링(FSS)을 적용하였다. 저자들은 또한 서로 다른 질서 영역을 구분하기 위해 $\ln L$에 따른 자화의 의존성과 스핀-스핀 상관 함수 $G(r)$의 거동을 조사하였다.

주요 기여 및 결과

1. 구별된 스케일링 영역(EnLRO 및 ReLRO)의 식별:
   주요 발견은 진정한 장거리 질서상이 균일하지 않으며 두 가지 구별된 스케일링 영역으로 구성된다는 것이다:

   • 강화된 장거리 질서 (EnLRO): 이 영역에서 자화 $M$은 시스템 크기 $L$이 증가함에 따라 증가한다. 이는 더 작은 $\sigma$ 값(더 강한 장거리 상호작용)에서 발생한다.
   • 감소된 장거리 질서 (ReLRO): 이 영역에서 자화 $M$은 시스템 크기 $L$이 증가함에 따라 감소한다(비록 유한한 값을 유지하지만), 이는 더 큰 $\sigma$ 값(더 약한 장거리 상호작용)에서 발생한다.
   • 크로스오버: 이 영역들 사이의 크로스오버는 임계 $\sigma_c$에서 발생한다. 이 근처에서 $M$은 $L$에 대해 비단조적 거동(감소 후 증가)을 보인다.

2. 크로스오버 지점의 정량화:

   • 1D XY 모델: 낮은 온도($T=0.01$)에서 크로스오버는 $\sigma_c \approx 1.575$에서 발생한다. 크로스오버 온도는 $\sigma$에 따라 증가한다.
   • 2D XY 모델: 낮은 온도에서 크로스오버는 $\sigma_c \approx 3.2$에서 발생한다.
   • 하이젠베르크 모델: 1D 및 2D 하이젠베르크 모델에서도 각각 $\sigma_c \approx 1.575$ 및 $\sigma_c \approx 3.22$에서 유사한 크로스오버가 관찰되었다. 저자들은 이러한 유사성이 이 현상이 연속적인 회전 대칭성을 공유하기 때문에 발생하는, 장거리 상호작용을 가진 연속 스핀 모델의 일반적인 특징임을 시사한다고 제안한다.

3. 상관 함수의 거동:

   • ReLRO 영역에서 스핀-스핀 상관 함수는 0이 아닌 값으로 대수적으로 감쇠한다 ($G(r) \sim r^{-\kappa}$).
   • EnLRO 영역에서 감쇠는 0이 아닌 상수에 접근하는 지수 함수($G(r) \sim k_1 + k_2 e^{-r/\xi}$)로 더 잘 설명된다. 상관 길이 $\xi$는 $\sigma_c$ 근처에서 불연속적인 변화를 보인다.

4. 상 평면 및 임계 지수:
   $T-\sigma$ 평면 상에 1D 및 2D 모델에 대한 종합적인 상 평면이 구축되었다. 저자들은 질서-무질서 전이를 위한 임계 지수 $\nu$와 $\gamma$를 추출하였다. 예를 들어, $\sigma=1.1$인 1D 모델에서 $\nu \approx 1.525$인 반면, $\sigma=1.8$에서는 $\nu \approx 2.65$이다. 임계 온도 $T_c$는 $\sigma \to D$일 때 발산한다.

5. 크로스오버 메커니즘:
   논문은 스핀 정렬, 스핀파(spin-wave) 들뜸, 그리고 국소 결함 사이의 경쟁에 기반한 메커니즘을 제안한다:

   • 저온: $\sigma$가 감소하면 스핀파의 들뜸과 전파가 억제되어, 시스템을 ReLRO에서 EnLRO로 유도한다.
   • 고온: 국소 결함(체인이나 격자를 하위 시스템으로 분리하는 결함)은 스핀파의 전파를 차단한다. 시스템 크기가 증가함에 따라 결함이 스핀파의 교란보다 더 빠르게 증식하여, 효과적으로 질서 영역을 고립시키고 국소적 질서를 강화함으로써 EnLRO를 유도한다.

의의 및 주장
본 논문은 장거리 질서상 내에 이전에 간과되었던 구별된 스케일링 거동이 존재함을 식별함으로써, 비-MWH 영역에서의 고전적 스핀 모델에 대한 이해를 정교화한다고 주장한다. 저자들은 EnLRO와 ReLRO의 구분이 아직 확인되지 않은 새로운 질서 매개변수를 가진 전통적인 열역학적 상 전이가 아니라, 시스템의 유한 크기 효과에 대한 반응의 중요한 변화라고 단언한다. 저자들은 EnLRO-ReLRO 전이가 전통적인 상 전이는 아니지만, 유한 크기 영역에서 이러한 효과가 두드러지는 리드베리(Rydberg) 원자 배열과 같은 실험적 시스템에 가치 있는 참조점을 제공한다는 점을 강조한다. 이 연구는 이론적 모델과 장거리 상호작용 시스템의 현재 실험적 구현에서 관찰되는 유한 크기 제약 사이의 간극을 메운다.