Probing the Ground State of the Antiferromagnetic Heisenberg Model on the Kagome Lattice using Geometrically Informed Variational Quantum Eigensolver

기술 요약: 기하학적 정보 기반 변분 양자 고유값 솔버(VQE)를 이용한 카고메 격자 상의 반강자성 하이젠베르크 모델 바닥 상태 탐사

문제 정의
본 논문은 기하학적 좌절(geometric frustration)과 잠재적인 양자 스핀 액체 상으로 알려진 시스템인 카고메 격자 상의 반강자성 하이젠베르크 모델의 바닥 상태를 찾는 과제를 다룹니다. 이러한 $k$-로컬 해밀토니안의 바닥 상태를 결정하는 것은 QMA-hard 문제로, 힐베르트 공간의 지수적 성장으로 인해 대규모 시스템에 대한 고전적 시뮬레이션이 불가능합니다. 변분 양자 고유값 솔버(VQE)는 잡음이 있는 중간 규모 양자(NISQ) 장치에 적합한 하이브리드 양자-고전 접근 방식을 제공하지만, 깊은 회로, 결맞음 및 비결맞음 오류에 대한 취약성, 그래디언트가 사라지는 배런 플래토(barren plateaus) 현상, 그리고 표현력과 하드웨어 효율성을 모두 갖춘 안사츠(ansatz) 설계의 어려움과 같은 상당한 장애물에 직면해 있습니다. 구체적으로, 저자들은 복잡한 문제 특화적 축소(reduction)나 최적화 루프 중의 광범한 오류 완화에 의존하지 않고, 실제 양자 하드웨어에서 스핀-스핀 상관관계 및 정적 구조 인자(static structure factor)와 같은 바닥 상태의 특성을 규명하는 것을 목표로 합니다.

방법론
저자들은 단일 삼각형과 카고메 스타(여섯 개의 꼭짓점을 공유하는 삼각형)라는 두 가지 근본적인 기하학적 구조에 초점을 맞춘, 카고메 격자에 특화된 기하학적 정보 기반 VQE 접근 방식을 제안합니다.

1. 기하학적 정보 기반 안사츠 (Geometrically Informed Ansatz):
   핵심 혁신은 맞춤형의 얕은 하드웨어 효율적 안사츠입니다. 해밀토니안 변분 안사츠(HVA)가 NISQ 장치에 너무 깊을 수 있는 것과 달리, 이 안사츠는 네이티브 하드웨어 연결성에 제한된 매개변수화된 단일 큐비트 회전($R_y$)과 얽힘 게이트 패턴을 사용하여 구축됩니다(SWAP 게이트를 피함).

   • 유클리드 매개변수 공간: 이 설계의 핵심 특징은 푸비니-스터디(Fubini-Study) 메트릭 텐서가 분석적으로 대각이며 상수임을 증명했다는 점입니다. 이는 매개변수 공간이 자연스럽게 유클리드적이고, 특이점이 없으며, 과잉 매개변수화로부터 자유로움을 보장합니다. 결과적으로 양자 자연 그래디언트(quantum natural gradient)는 표준 그래디언트와 동일해지며, 이는 최적화 지형을 단순화합니다.
   • 하드웨어 호환성: 회로 깊이는 최소화되며(하나의 얽힘 레이어), $R_y$ 게이트의 배치를 교대로 배치함으로써 매개변수의 수를 줄였습니다.

2. 최적화 전략 (AQNGD/AGD):
   저자들은 적응형 그래디언트 하강법(AGD) 옵티마이저를 채택합니다. 안사츠가 대각이며 일정한 푸비니-스터디 메트릭을 생성하기 때문에, 적응형 양자 자연 그래디언트 하강법(AQNGD)은 AGD로 축소됩니다. 학습률은 아르미호 백트래킹 규칙(Armijo backtracking rule)을 사용하여 동적으로 적응됩니다. 이를 통해 초기 단계에서는 큰 보폭을 가져가 수렴을 가속화하면서도, 옵티마이저가 최솟값에 접근함에 따라 자동으로 보폭을 줄여 좁은 기하학적 구조에서의 정체를 방지할 수 있습니다. 이는 SPSA(Simultaneous Perturbation Stochastic Approximation)와 비교됩니다.

3. 오류 완화 (사후 처리):
   오류 완화는 반복적인 훈련 과정 중이 아니라, 최종 에너지 추정치와 구조적 특성을 정밀하게 하기 위해 VQE 최적화 이후에 적용됩니다.

   • 판독 오류 완화 (REM): 분할 응답 행렬 접근 방식을 사용하여 SPAM(상태 준비 및 측정) 오류를 역전시키며, 이는 전체 행렬 역행렬 계산보다 희소 측정(sparse measurements)에 더 효율적으로 확장됩니다.
   • 제로 노이즈 외삽 (ZNE): 저자들은 이차 외삽(및 베이지안 다항 회귀)을 사용하는 ZNE를 적용하여 결과를 제로 노이즈 한계로 외삽합니다. 그들은 ZNE가 엄격하게 레이리-리츠(Rayleigh-Ritz) 변분 원리를 따르지는 않지만, 정확도를 높이기 위해 사용되었다고 언급합니다.

주요 결과
연구는 실제 양자 하드웨어(IBM Oslo, IBM Kyoto, IBM Torino, IBM Algiers)에서 수행되었으며 상태 벡터 시뮬레이션을 통해 검증되었습니다.

• 수렴 및 최적화:

  • 3-큐비트 삼각형 시스템에서, AGD로 최적화된 VQE는 단 19회의 반복만으로 높은 충실도(시뮬레이션에서 99.89%)를 보이며 바닥 상태 에너지($E \approx -3$)로 매끄럽게 수렴했습니다. 최종 매개변수는 이론적 정확값과 밀접하게 일치했습니다.
  • 12-큐비트 카고메 스타의 경우, AGD는 SPSA보다 우수한 수렴 성능을 보였으며, 바닥 상태 에너지($E \approx -18$)에 도달하기 위해 더 적은 반복 횟수를 필요로 했습니다. 시뮬레이션에서의 최종 충실도는 96.37%였습니다.
  • 아르미호 규칙을 통한 단계 크기 적응은 비용 지형을 효과적으로 탐색하는 것으로 나타났으며, 큰 단계로 시작하여 수렴이 접근함에 따라 정밀하게 조정되었습니다.

• 구조적 특성의 노이즈 탄력성:

  • 중요한 발견은 스핀-스핀 상관관계와 정적 구조 인자($S(q)$)가 노이즈에 탄력적이라는 점입니다. 오류 완화 없이도, 실험적 $S(q)$ 맵은 바닥 상태의 질적 특징(예: 강자성 질서가 없음을 나타내는 제로 모멘텀에서의 억제)을 유지했습니다.
  • 노이즈가 있는 실험적 $S(q)$와 정확한 이론적 $S(q)$ 사이의 피어슨 상관 계수는 여전히 매우 높았습니다(>99.9%).
  • 시뮬레이션 실험은 정적 구조 인자의 노이즈 수준에 대한 "계단형" 의존성을 보여주었습니다: 상관 구조는 임계 노이즈 임계값까지 견고하게 유지되다가, 그 이후에는 급격히 저하되며 이는 상전이와 유사한 모습을 보입니다. 이는 에너지가 노이즈에 따라 단조 감소하는 것과는 대조적입니다.

• 오류 완화 효능:

  • REM과 ZNE를 최적화 후에 적용하는 것은 에너지 추정치를 크게 개선했습니다. 카고메 스타의 경우, 가공되지 않은 에너지는 $\approx -14.9$였으나, REM만 적용해도 $\approx -17.3$로 개선되었고, ZNE(이차)를 적용했을 때 정확한 값인 $-18$에 더 가까워졌습니다.
  • 저자들은 REM과 ZNE를 결합했을 때 때때로 바닥 에너지를 초과(overshooting)하여 변분 경계(variational bound)를 위반하는 현상을 관찰했는데, 이는 ZNE의 비변분적 성질과 일치하는 동작입니다. ZNE 단독 사용이 테스트된 구성 중 가장 높은 정확도를 제공하는 경우가 많았습니다.

의의 및 주장
본 논문은 기하학적으로 정보가 담긴 안사츠 설계가 자연스럽게 유클리드적 매개변수 공간을 보장함으로써, 복잡하고 문제 특화적인 축소 없이도 현재의 NISQ 하드웨어에서 효율적인 VQE 훈련을 가능하게 한다고 주장합니다.

• 강건성: 본 연구는 노이즈가 있는 양자 하드웨어로부터 충분한 정확도로 스핀-스핀 상관관계 및 정적 구조 인자와 같은 의미 있는 물리적 특성을 추출할 수 있으며, 이를 통해 오류 완화 과정 중의 최적화 루프 없이도 바닥 상태를 규명할 수 있음을 입증했습니다.
• 최적화 효율성: 이 특정 하드웨어 효율적 안사츠에 대해 적응형 단계 크기를 가진 AGD를 사용하는 것이 SPSA보다 더 효율적임을 입증하였으며, 반복 횟수 측면에서 더 빠른 수렴을 달ach했습니다.
• 확장성: 깊은 회로와 훈련 중의 무거운 오류 완화를 피함으로써, 이 접근 방식은 좌절된 양자 자석을 탐사하기 위한 확장 가능한 경로를 제공합니다. 저자들은 이 프레임워크가 더 큰 카고메 격자 및 다른 좌절된 양자 시스템으로 확장될 수 있다고 제안합니다.
• 겸허함: 저자들은 자신의 방법이 작은 기본 단위(삼각형 및 스타)를 성공적으로 규명했지만, 무한 카고메 격자의 일반적인 바닥 상태는 여전히 문헌상의 미해결 과제임을 인정합니다. 그들은 이 문제를 해결했다고 주장하는 것이 아니라, 현재의 장치로 이러한 시스템을 조사하기 위한 실행 가능한 방법론을 제시한 것입니다. 또한 ZNE가 에너지 추정치를 개선하지만 변분 원리를 보존하지 않으며, REM과 ZNE의 결합은 주의 깊은 해석이 필요하다는 점을 명시했습니다.