양자 어닐링으로 고전 임계 현상 재현

초록

본 논문은 D‑Wave 양자 어닐러를 이용해 2차원 피일드‑업 도미노(PUD) 모델의 상전이와 임계 현상을 정확히 재현한다. 에너지 스케일을 조절해 효과 온도를 제어하고, Binder 누적량과 유한 크기 스케일링을 적용해 임계 온도와 임계 지수를 추정한다. 실험 결과는 정확 해와 일치하며, 양자 어닐링이 고전 통계 물리 시뮬레이션에 유효한 도구임을 입증한다.

상세 분석

이 연구는 양자 어닐링(Quantum Annealing, QA)이 고전 통계 물리 모델, 특히 임계 현상이 나타나는 시스템을 시뮬레이션하는 데 실용적이고 정확한 방법이 될 수 있음을 실증한다. 저자들은 먼저 피일드‑업 도미노(Piled‑Up Domino, PUD) 모델을 선택했는데, 이 모델은 수직 결합 상수의 비율 s에 따라 완전 반강제(Villain) 모델과 2D 페리마그넷 이징 모델 사이를 연속적으로 연결한다. s = 0에서는 전통적인 페리마그넷 이징 모델이, s = 1에서는 완전 반강제 모델이 되며, 중간값에서는 부분적인 좌절(frustration)이 도입된다. PUD 모델은 전이 온도 T_c(s)가 정확히 알려진 드문 경우이며, 페리마그넷‑파라마그넷, 반강제‑파라마그넷, 그리고 반강제‑페리마그넷 사이의 세 가지 상이 s‑T 평면에 명확히 구분된다.

양자 어닐러가 실제로는 고정된 장치 온도 T_sampler 에서 Boltzmann 분포를 샘플링한다는 점을 이용해, 입력 해밀토니안에 스케일링 인자 J를 곱함으로써 효과 온도 T_eff ∝ 1/J를 선형적으로 조절한다. 이는 기존에 물리적 온도를 직접 제어해야 했던 Monte Carlo 방법과 달리, 하드웨어 설정을 바꾸지 않고도 온도 스캔이 가능하게 한다는 큰 장점을 제공한다. 저자들은 D‑Wave Advantage2 프로토타입 2.6을 사용해 L = 6, 8, 10, 12 (총 144 스핀) 크기의 토로이달 격자에 대해 10 000번씩 샘플링하고, 평균 자화 ⟨|m|⟩와 교번 자화 ⟨|m_AFM|⟩를 계산해 전반적인 상도표를 재구성했다. 결과는 정확 해와 거의 일치했으며, 특히 임계선이 명확히 드러났다.

임계점을 정밀하게 찾기 위해 Binder 누적량 U = 1 − ⟨m⁴⟩/(3⟨m²⟩²) 를 시스템 크기별로 계산하고, 서로 다른 L에 대한 U(J⁻¹) 곡선이 교차하는 지점을 임계 스케일 J_c⁻¹(s) 로 정의했다. 교차점은 약간의 산란을 보였지만, 이는 제한된 시스템 크기와 장치 노이즈에 기인한다. 중요한 점은 J_c⁻¹(s)와 정확 해에서 얻은 T_c(s) 사이에 선형 관계가 성립한다는 것이며, 이는 앞서 제시한 T_eff ∝ 1/J 가 실험적으로 검증된 결과이다.

또한, 저자들은 자화 민감도 χ 와 열용량 C_V 를 각각 χ/β = N(⟨m²⟩ − ⟨m⟩²)와 C_V/β² = N(⟨E²⟩ − ⟨E⟩²) 로 정의하고, 시스템 크기별 피크 위치와 높이를 분석했다. 피크가 시스템이 커질수록 더 뾰족해지고, 피크 위치가 낮은 J⁻¹(즉, 높은 온도) 쪽으로 이동하는 전형적인 유한 크기 스케일링 현상이 관찰되었다. 이를 바탕으로 χ와 C_V의 스케일링 법칙 χ ∝ L^{γ/ν} · f(tL^{1/ν}) 와 C_V ∝ L^{α/ν} · g(tL^{1/ν}) 를 적용해 임계 지수 ν, γ, α 를 추정했으며, 얻어진 값은 2D 이징 모델의 알려진 지수와 통계적으로 일치했다.

전체적으로 이 논문은 (1) 양자 어닐러가 효과 온도 제어를 통해 전통적인 Monte Carlo와 동등한 온도 스캔을 수행할 수 있음을, (2) Binder 누적량과 유한 크기 스케일링을 이용해 임계 온도와 임계 지수를 정확히 추정할 수 있음을, (3) 장치의 비이상적인 샘플링(노이즈, 편향 등)에도 불구하고 고전 통계 물리의 핵심 현상을 재현할 수 있음을 입증한다. 이는 양자 어닐링이 고전 통계 물리 시뮬레이션, 특히 좌절이 강한 시스템의 상전이 연구에 새로운 도구로 자리매김할 가능성을 시사한다.