Hamiltonian spectra in quantum computers through the generalized eigenvalue method

기술 요약: 일반화된 고유값 방법을 통한 양자 컴퓨터의 해밀토니안 스펙트럼 분석

문제 정의
해밀토니안의 에너지 스펙트럼을 결정하는 것은 양자 물리학의 근본적인 목표이며, 특히 섭동론적 방법이 실패하는 강하게 결합된 장론(strongly coupled field theories)에서 매우 중요하다. 고전적 격자 장론은 유클리드 상관 함수(Euclidean correlators)의 몬테카를로 평가를 활용하여 스펙트럼을 추출하지만, 이 접근 방식은 "양날의 검"과 같은 문제를 안고 있다. 즉, 허수 시간(imaginary time)이 커질수록 바닥 상태가 지배적이 되지만, 들뜬 상태(excited states)는 지수적으로 억제되어 신호 대 잡음비(signal-to-noise ratio)의 저하로 인해 추출이 어려워진다.

양자 컴퓨터는 실시간 상관 함수(real-time correlators)를 자연스럽게 생성할 수 있다는 점에서 뚜렷한 이점을 제공하며, 여기서는 모든 상태가 감쇠되지 않는 진동하는 가중치와 함께 기여한다. 그러나 이러한 진동 신호로부터 에너지 고유값을 추출하는 작업은 통상적으로 푸리에 변환(Fourier transforms)에 의존한다. 푸리에 변환의 해상도는 전체 진화 시간($1/T$)의 역수로 제한된다. 현재의 양자 하드웨어는 유한한 결맞음 시간(coherence times)과 축적되는 게이트/트로터 오차(gate/Trotter errors)로 인해, 푸리에 분석을 통해 거의 퇴행적인(nearly degenerate) 레벨을 구별하거나 링잉 아티팩트(ringing artifacts)를 실제 피크와 구분하기 위해 필요한 긴 진화 시간을 유지하기 어렵다.

방법론
저자들은 양자 컴퓨터에서 생성된 실시간 상관 함수에 일반화된 고유값 문제(Generalized Eigenvalue Problem, GEVP) — 격자 장론의 표준 기술 — 를 적용할 것을 제안한다. 핵심 방법론은 다음과 같다:

1. 연산자 기저 및 상관 함수 행렬: 단일 상관 함수를 분석하는 대신, $C_{ij}(t) = \langle \Omega | O_i(t) O_j^\dagger(0) | \Omega \rangle$ 형태의 상관 함수 행렬을 구성한다. 여기서 $O_i$는 불완전한 보간 연산자(interpolating operators)의 집합이며, $|\Omega\rangle$는 근사적인 바닥 상태이다.
2. 이론적 프레임워크: 저자들은 행렬 곱 $C^{-1}(0)C(t)$의 고유값 $\lambda_k(t)$가 $e^{-iE_k t}$처럼 거동함을 입증한다. 여기서 $E_k$는 에너지 고유값이다. 유효 에너지(effective energy)는 $E_k^{\text{eff}}(t) = \frac{i}{t} \log(\lambda_k(t))$를 통해 추출된다. 이론적 분석에 따르면, 연산자 기저가 관련 고유 상태의 부분 공간을 충분히 포괄한다면, 정확한 에너지에 대한 보정 항은 $1/t$ (또는 파괴적 간섭으로 인해 그보다 빠르게) 감소한다.
3. 양자 구현: 단일 보조 큐비트(ancilla qubit)를 사용하여 양자 컴퓨터에서 파울리 연산자의 비대각 상관 함수를 측정하는 구체적인 알고리즘을 상세히 기술한다. 절차는 다음과 같다:
   • 근사 바닥 상태 $|\Omega\rangle$와 보조 큐비트를 준비한다.
   • 연산자 $O_j$를 포함하는 제어 연산을 통해 시스템과 보조 큐비트를 얽히게 한다.
   • 시스템을 시간 진화시킨다.
   • 상관 함수의 실수부와 허수부를 재구성하기 위해 특정 관측량 $X \otimes O_i$ 및 $Y \otimes O_i$를 측정한다.
4. 비교 벤치마크: 이 방법은 깁스 현상(Gibbs phenomena)을 완화하기 위한 댐핑 인자(damping factor)를 활용한 상관 함수 행렬 트레이스의 이산 푸리에 변환(DFT)과 벤치마킹된다.

주요 기여 및 결과
본 논문은 이 접근 방식을 물리적으로 유의미한 장론(예: 점근적 자유도를 가진 이론)과 유사한 특징을 공유하는 퍼지 $\sigma$-모델(fuzzy $\sigma$-model, 하이젠베르크 콤과 동일)을 사용하여 검증한다. 연구에는 세 가지 뚜렷한 테스트 환경이 사용되었다:

• 정확한 시뮬레이션 ($L=8$): 정확한 시간 진화(노이즈나 트로터 오차 없음)를 사용하여, 세 가지 단계적으로 복잡해지는 연산자 집합을 통해 8개, 36개, 276개의 에너지 레벨을 성공적으로 분리하였다. 결과는 최대 진화 시간 $T_{\text{max}} = 5.6$ 내에서 낮은 에너지 상태들에 대해 1% 미만의 정밀도를 달お성하였다. 또한, 이 방법은 모델의 $O(3)$ 대칭성에 의해 결정되는 퇴행 패턴을 정확하게 재현하였다.
• 게이트 기반 시뮬레이션 ($L=4$): 저자들은 트로터 오차의 영향을 분석하였다. 다양한 시간 단계($\Delta t$)로부터의 결과를 $\Delta t \to 0$으로 외삽함으로써, 이 방법이 견고하며 편차가 단계 크기에 따라 이차적으로 스케일링됨을 확인하였다.
• 양자 하드웨어 ($L=4$ on IonQ Forte-Enterprise): 실제 양자 하드웨어에서 이 방법이 실행되었다. 단일 보조 큐비트와 최소한의 연산자 집합을 사용하여 4개의 낮은 에너지 레벨을 분리하였다. 추출된 에너지는 정확한 대각화(exact diagonalization) 결과와 몇 퍼센트 이내의 오차로 일치하였다. 잔여 불일치는 GEVP 방법 자체의 실패라기보다는, 거친 시간 단계($\Delta t = 0.4$)로 인한 트로터 오차와 하드웨어 노이즈/결맞음 해제(decoherence)의 결합 때문인 것으로 분석되었다.

의의 및 주장
본 논문은 GEVP 접근 방식이 스펙트럼 추출에 있어 푸리에 변환 방법보다 "실질적으로 더 효율적"이라고 주장한다. 구체적으로, 이 방법은 푸리에 변환이 피크와 아티팩트를 구분하지 못하는 시간 범위에서도 낮은 에너지 레벨을 분리할 수 있게 해준다.

저자들은 이 방법이 완벽한 연산자나 정확한 바닥 상태를 요구하지 않는다는 점을 강조한다. 대신, 관심 있는 상태들을 대략적으로 포괄하는 연산자 기저와 좋은 중첩(overlap)을 가진 근사 바닥 상태만을 필요로 한다. 비록 여기서는 필요한 상태들을 구축하기 위해 강한 결합 한계(strong-coupling limits)를 사용하였으나, 저자들은 이 방법이 일반적이며 변분적(variational) 또는 단열적(adiabatic) 상태 준비 방식에 모두 적용 가능하다고 언급하였다.

이 연구의 주요 의의는 이 기술이 유클리드 접근법을 가로막는 "부호 문제(sign problem)"가 존재하는 이론들, 특히 유한 밀도 게이지 이론 및 실시간 수송 현상에서의 분광학(spectroscopy)을 가능하게 할 잠재력에 있다. 본 연구는 실시간 양자 시뮬레이션이 주파수 영역 분석에 필요한 과도하게 긴 진화 시간 없이도 정밀한 스펙트럼 데이터를 산출할 수 있음을 보여주는 개념 증명(proof-of-principle) 역할을 한다.