A fidelity metric for quantum annealing benchmarked by extreme scaling quantum Monte-Carlo simulations

기술 요약: 극한 스케일링 양자 몬테카를로 시뮬레이션에 의해 벤치마킹된 양자 어닐링의 충실도 지표

문제 정의
현재 양자 어닐링(QA)의 벤치마크는 주로 실용적인 측면에 치중되어 있으며, 특정 최적화 문제(통상적으로 QUTO로 정식화됨)를 해결하는 속도나 성공률만을 측정합니다. 저자들은 이러한 접근 방식이 결과가 특정 문제 인스턴스와 참조로 사용되는 고전 알고리즘에 크게 의존하기 때문에 "본질적으로 취약하다"고 주장합니다. 또한, 이러한 방법론은 어닐링 과정의 내재적 한계(예: 랜다우-제너 전이를 유발하는 지수적으로 소멸하는 에너지 갭)와 외재적 하드웨어 한계(예: 노이즈, 결맞음, 유한 온도)를 구분하지 못합니다. 결과적으로, 최종 솔루션의 품질과는 독립적으로 양자 어닐링 과정 자체의 충실도를 정량화할 수 있는 견고한 지표가 부족한 실정입니다.

방법론
본 논문은 새로운 충실도 지표로서, 어닐링 과정 중 상태 방정식(equation of state) $E(h_x) \equiv \langle \hat{H} \rangle$의 정확도($\epsilon$)를 제안합니다. 이 지표는 게이트 기반 양자 컴퓨팅의 "게이트당 충실도(fidelity-per-gate)"와 유사합니다. 실험적으로 상태 방정식을 얻으려면 어닐링 스케줄의 여러 지점에서 종방향 에너지($\langle \hat{H}_z \rangle$)와 횡방향 자기장 에너지($\langle \hat{H}_x \rangle$)를 모두 측정해야 합니다.

$\epsilon$의 벤치마크 값을 설정하고 고전 시뮬레이션의 한계를 테스트하기 위해, 저자들은 리드베리 원자 시스템에서 양자 몬테카를로(QMC) 기법을 사용하는 두 가지 변형된 시뮬레이드 양자 어닐링(SQA)을 채택했습니다:

1. 변분 양자 몬테카를로 (VMC):

   • 안사츠(Ansatz): 저자들은 시뮬레이드 열적 어닐링에서 영감을 얻은 "열적 안사츠"를 사용합니다. 파동함수는 $\psi(n) \propto e^{-E_{\text{eff}}(n)/2}$로 모델링되며, 여기서 유효 에너지 $E_{\text{eff}}(n)$은 수정된 고전적 에너지 범함수입니다.
   • 변형: 저자들은 (1) 단순 열적 안사츠($E_{\text{eff}} = \beta_{\text{eff}} F(n)$)와 (2) 대칭성 깨짐을 포착하기 위한 격자형 열적 안사츠($E_{\text{eff}} = \beta_{\text{eff}} F(n) + \sum \beta_k N_k(n)$)를 테스트합니다.
   • 규모: 이 접근 방식은 계산적으로 가벼워 단일 CPU에서 최대 100,000,000개의 원자($N=10^8$)까지 시뮬레이션이 가능합니다.
   • 정확도: 약 $\epsilon \sim 10^{-2} - 10^{-3}$의 정확도를 달성합니다.

2. 그린 함수 몬테카를로 (GFMC):

   • 방법: VMC 이후 단계의 방법으로, 참된 바닥 상태로 투영하기 위해 허수 시간(imaginary time) 동안 변분 안사츠를 진화시킵니다. 이는 편향되지 않았으며 횡방향 이싱 모델(transverse field Ising models)에 대해 사인 문제(sign problem)가 없습니다.
   • 규모: 최대 100,000개의 원자($N=10^5$)까지 시뮬레이션을 수행합니다.
   • 정확도: $\epsilon \sim 10^{-4}$의 더 높은 정밀도를 달성합니다.

저자들은 이 방법들을 사용하여 정사각형 및 삼각형 격자에서의 리드베리 원자와 실제 QUBO 인스턴스("fallen angel problem", 금융 리스크 관리 관련)에 대한 상태 방정식을 계산합니다.

주요 결과

• 확장성 및 정밀도: 열적 안사츠를 사용한 VMC 접근 방식은 상대 오차 $\delta E/E_0 < 10^{-2}$로 최대 $10^8$개의 원자를 성공적으로 시뮬레이션했습니다. GFMC 접근 방식은 $N=10^5$에 대해 $\epsilon \sim 10^{-4}$에 도달했습니다.
• 오차 국소화: 상태 방정식의 상대 오차는 에너지 갭이 닫히는 위치(상전이 지점)에서 정확히 정점을 찍으며, 이는 랜다우-제너 전이가 발생하는 영역과 일치합니다. 이러한 임계점 이외의 영역에서는 오차가 현저히 낮습니다.
• 실험과의 비교: 저자들은 현재의 리드베리 원자 실험 플랫폼이 이러한 고전적 벤치마크보다 훨씬 낮은 정확도로 작동하고 있다고 추정합니다. 고전 시뮬레이션은 시스템 크기와 정밀도 모두에서 현재의 실험 플랫폼을 수 차수(orders of magnitude) 앞섭니다.
• QUBO 벤치마킹: "fallen angel" QUBO 인스턴스의 경우, 단순 열적 안사츠는 $\approx 10^{-3}$의 상대 오차를 나타냈습니다. 저자들은 최종 솔루션 오차($\delta F/F$)는 샷(shot) 수를 늘림으로써 임의로 줄일 수 있는 반면, 상태 방정식 오차($\delta E/E_0$)는 하드웨어나 알고리즘의 역량을 보여주는 더 본질적인 척도라고 언급합니다.
• 열적 vs 양자: 본 연구는 정확도 $\epsilon \sim 10^{-2} - 10^{-3}$ 범위 내에서 양자 어닐러는 사실상 그에 대응하는 고전적 열적 모델과 구별할 수 없음을 시사하며, 이는 현재의 실험적 노이즈 수준이 진정한 양자 우위를 가리고 있을 수 있음을 의미합니다.

의의 및 주장
본 논문은 상태 방정식의 정확도를 측정하는 것이 다음과 같은 이유로 QA 벤치마킹의 중요한 진전이라고 주장합니다:

1. 이상적인 단열 진화로부터의 편차를 일으키는 지배적인 메커니즘(내재적 갭 문제와 외재적 노이즈의 구분)을 진단하는 직접적인 도구를 제공합니다.
2. 선택된 특정 QUBO 인스턴스의 난이도에 의존하지 않는 하드웨어 불가지론적(hardware-agnostic) 척도를 제공합니다.
3. 미래 하드웨어에 엄격한 제약을 설정합니다: 양자 우위를 입증하려면 실험 플랫폼은 이러한 고전적 벤치마크(특히 대규모 시스템에 대해 $\epsilon < 10^{-4}$)를 능가하는 상태 방정식 정확도를 달성해야 합니다.

저자들은 "기만적일 정도로 단순한" 열적 변분 안사츠가 QA를 시뮬레이션하는 강력한 도구임을 결론지으며, 양자 우위를 이해하는 길은 양자 요동이 가장 정교한 고전적 열적 시뮬레이션보다 뚜렷한 이점을 제공하는 영역을 식별하는 데 있음을 시사합니다. 저자들은 QA가 현재 고전적 휴리스틱보다 우월하다고 주장하는 것이 아니라, 향로의 가설을 엄격하게 테스트할 수 있는 도구를 제공하고자 합니다.