양자 시뮬레이션으로 기하학적 좌절 자석의 스케일링 이점

초록

본 연구는 D‑Wave 초전도 플럭스 큐비트 양자 어닐링 프로세서를 이용해, 원통형 경계조건을 갖는 1440 스핀 규모의 기하학적 좌절 격자 시스템을 시뮬레이션하였다. 초기 상태를 위상적 장애가 있는 형태로 준비함으로써 마이크로초 수준의 완화 시간을 관찰하고, 동일 조건에서 경로 적분 몬테카를로(PIMC) 시뮬레이션과 비교했을 때 시스템 크기와 온도가 낮아질수록 최대 10⁶ 배에 달하는 동적 속도 향상을 입증하였다. 이는 스투키오스틱 해밀토니안을 갖는 경우에도 양자 어닐링이 고전적 샘플링보다 실질적인 스케일링 이점을 제공한다는 실험적 증거이다.

상세 분석

이 논문은 양자 어닐링(QA)과 경로 적분 몬테카를로(PIMC) 사이의 동적 차이를 정량적으로 평가하기 위해, 2차원 사각‑팔각 격자 구조에 대한 TFIM(Transverse‑Field Ising Model) 해밀토니안을 구현하였다. 격자 내 각 플라크는 반강자성 결합이 홀수 개 존재해 완전 만족이 불가능한 ‘좌절’ 상태를 만들며, 이는 여섯 가지 고전적 기저 상태(의사스핀 방향)로 표현된다. 양자 플럭스(Γ)를 도입하면 네‑스핀 체인이 GHZ‑형 초포지션을 형성해 의사스핀 간의 강한 페르미온 상호작용을 유도, 결국 장거리 페르미온 순서를 발생시킨다.

실험적 핵심은 ‘위상적 장애’를 가진 초기 상태(시계/반시계 방향으로 의사스핀이 전체 원통형 경계에 걸쳐 꼬인 상태)를 준비함으로써, 시스템이 전역적인 ‘풀어내기’ 과정을 거쳐야만 평형에 도달하도록 만든 점이다. 이러한 초기 조건은 QA의 완화 속도를 마이크로초 수준으로 늦추어, 기존 장치의 1 µs 제어 한계 내에서 시간 해상도를 확보한다.

QA 측정은 ‘역 어닐링’ 프로토콜을 사용해 s = 1(플럭스 억제)에서 지정된 고전적 상태로 초기화한 뒤, s* < 1으로 급격히 전이하고 일정 시간(tₚ = 1–4 µs) 동안 시스템을 진화시킨다. 이후 다시 s = 1로 급속 퀀치하고 읽어내어, 다수 반복을 통해 평균 의사스핀 순서 매개변수 m(t)를 추출한다. PIMC는 연속시간 알고리즘으로 4‑스핀 체인 클러스터 업데이트를 적용해 동일한 초기 조건을 모사하고, Monte Carlo sweep를 시간 단위로 매핑한다.

두 시뮬레이션 모두 m(t)의 지수적 수렴을 관찰했으며, QA는 동일 정확도(오차 ≤ 0.03) 내에서 PIMC 대비 평균 10⁶ 배 빠른 수렴 시간을 보였다. 특히 온도 T가 낮고 시스템 폭 L이 클수록, 그리고 Γ/J 비율이 중간(≈0.7)일 때 이점이 극대화되었다. 저온·저Γ 영역에서는 하드웨어 내 불균일성(제조 공정에 의한 J_ij 변동)으로 인해 QA 순서가 억제되어 PIMC와 차이가 발생했으나, 이는 모델링 상의 ‘퀀치 보정’으로 보정 가능하였다.

결과적으로, 스투키오스틱 해밀토니안을 갖는 복잡한 좌절 시스템에서도 양자 어닐링이 고전적 샘플링보다 실질적인 동적 스케일링 이점을 제공함을 실험적으로 입증하였다. 이는 양자 어닐링이 단순 최적화 문제를 넘어, 물리적 평형 통계량을 빠르게 추정할 수 있는 ‘양자 시뮬레이션 엔진’으로 활용 가능함을 시사한다.