비평형 임계 급랭 및 어닐링 동역학 연구

초록

본 논문은 1차원 장거리 상호작용 이징 모델(σ=0.75)을 대상으로, 무작위 초기 상태에서 임계점으로 급랭한 경우와 정렬된 초기 상태에서 임계점으로 급히 어닐링한 경우의 비평형 동역학을 대규모 몬테카를로 시뮬레이션으로 조사한다. 임계 온도와 정적·동적 임계 지수를 전력법칙을 이용해 추정하고, 유한 범위 스케일링을 통해 무한 상호작용 한계의 임계 온도를 구한다. 정적 지수는 RG 예측과 일치하지만, 동적 지수 z는 급랭과 어닐링에서 서로 다른 값을 보이며, 이는 σ≈1에서 발생하는 Kosterlitz‑Thouless 전이와 사용된 Metropolis 알고리즘에 기인한다는 결론을 제시한다.

상세 분석

이 연구는 장거리 상호작용을 갖는 1차원 이징 모델에 대해 비평형 임계 현상을 정량적으로 분석한 점에서 의미가 크다. 먼저, 상호작용 형태를 $J(r)\sim r^{-(d+\sigma)}$ 로 두고 σ=0.75를 선택함으로써, 기존에 알려진 임계 차원 $d_c=2\sigma$와 비교했을 때 1차원에서도 비평형 임계 현상이 나타날 수 있는 파라미터 영역을 탐색한다. 급랭(quench) 실험에서는 무한 온도에서 완전히 무작위인 초기 스핀 배열을 사용하고, 시스템을 즉시 임계 온도 $T_c$ 로 내려놓는다. 이때 관측되는 물리량—예를 들어 초기 마그네틱 모멘트 $M(t)$, 자기상관함수 $C(t)$, 그리고 제2 모멘트 $U(t)$—은 시간에 대한 전력법칙 $M(t)\sim t^{-\beta/\nu z}$ 등으로 스케일링한다. 이러한 전력법칙을 통해 $\beta/\nu$, $\gamma/\nu$, $1/\nu$, 그리고 동적 지수 $z$를 동시에 추정한다.

반면 어닐링(annealing) 실험에서는 초기 상태를 완전 정렬된(모든 스핀이 +1) 상태로 두고, 동일하게 $T_c$ 로 급히 온도를 변환한다. 여기서는 초기 마그네틱 모멘트가 1에서 시작하므로, 초기 감쇠 속도와 장기적인 성장/감소 패턴이 급랭 경우와 다르게 나타난다. 특히, 초기 슬립 지수 $\theta$가 두 실험에서 서로 다른 값을 보이며, 이는 초기 조건에 대한 민감도가 동적 지수에 크게 영향을 미친다는 점을 시사한다.

정적 지수들의 경우, 저자들은 유한 범위 스케일링(finite‑range scaling) 기법을 적용해 서로 다른 상호작용 절단 거리 $r_{\text{cut}}$ 에 대해 측정된 $T_c(r_{\text{cut}})$ 를 $r_{\text{cut}}\to\infty$ 로 외삽하였다. 이 과정에서 얻어진 $T_c$ 값은 기존 RG 예측 및 이전 수치 연구와 일치함을 확인한다.

동적 지수 $z$에 대한 가장 흥미로운 결과는 급랭과 어닐링에서 서로 다른 값을 얻었다는 점이다. 급랭에서는 $z_{\text{quench}}$ 가 RG 예측보다 크게 나타났으며, 이는 σ=1 근처에서 발생하는 Kosterlitz‑Thouless 전이가 장거리 상호작용의 유효 차원을 변형시켜 동적 스케일링에 비정상적인 영향을 미친 것으로 해석된다. 반면 어닐링에서는 $z_{\text{anneal}}$ 가 RG 예측치와 근접했으며, 이는 정렬된 초기 상태가 시스템을 보다 “균일한” 동역학 경로로 이끌어 전통적인 동적 스케일링 법칙을 따르게 만든 것으로 보인다. 또한, 사용된 Metropolis 업데이트가 비보존적이며, 전통적인 모델 A 동역학과 차이를 보이기 때문에 $z$와 초기 슬립 지수 $\theta$가 이론적 RG 값과 약간의 편차를 보이는 것이 자연스럽다.

전체적으로, 이 논문은 비평형 초기 조건이 장거리 상호작용 이징 모델의 동적 지수에 미치는 영향을 체계적으로 조사했으며, 특히 급랭과 어닐링 사이의 차이를 통해 초기 상태와 상호작용 범위가 동적 스케일링에 미치는 복합적인 역할을 밝히는 데 성공했다. 이러한 결과는 장거리 상호작용을 갖는 물리계(예: 유전체, 초전도체, 그리고 일부 양자 스핀 체인)에서 비평형 현상을 이해하고, 실험적 프로토콜을 설계하는 데 중요한 지침을 제공한다.