Ground State Energy via Adiabatic Evolution and Phase Measurement for a Molecular Hamiltonian on an Ion-Trap Quantum Computer

기술 요약: 이온 트랩 양자 컴퓨터 상의 분자 해밀토니안에 대한 단열 진화 및 위상 측정을 통한 바닥 상태 에너지 측정

문제 정의
분자 바닥 상태 에너지(GSE)를 정확하게 추정하는 것은 양자 컴퓨팅의 주요 응용 분야이지만, 정확한 양자 상태 준비, 고정밀 관측량 추정, 그리고 하드웨어 노이즈 완화의 어려움으로 인해 실제 구현에 제약이 따릅니다. 고전적인 방법인 하트리-포크(Hartree–Fock, HF)나 결합 클러스터(coupled cluster) 방식은 강하게 상관된 시스템에서 필수적인 전자 상관 효과를 포착하는 데 실패할 수 있으며, 더 정교한 고전적 방법들은 지수적으로 스케일링되는 문제를 가집니다. 양자 알고리즘, 특히 단열 상태 준비(Adiabatic State Preparation, ASP)와 양자 위상 추정(Quantum Phase Estimation, QPE)은 체계적인 개선을 제공하지만, 파라미터 최적화, 배런 플래토(barren plateaus), 그리고 제어된 $U^{2^k}$ 연산에 필요한 높은 회로 깊이 등의 과제에 직면해 있습니다. 또한, 소프트웨어를 통해 완화 가능한 오류와 하드웨어 수준의 감소가 필요한 오류를 구분하는 것은 이 분야의 발전에 매우 중요합니다.

방법론
저자들은 $H_3^+$ 분자의 GSE(6개의 큐비트에 인코딩됨)를 추정하기 위해 Quantinuum의 H1-1 이온 트랩 양자 컴퓨터에서 순수 양자 기반의 ab initio 프로토콜을 구현했습니다. 워크플로우는 세 가지 주요 구성 요소로 이루어집니다:

1. 단열 상태 준비 (ASP): 이산화 오차(discretization errors)를 유발하는 표준 Trotter-Suzuki 분해 대신, 저자들은 TETRIS(Time Evolution Through Random Independent Sampling) 알고리즘을 채택했습니다. 이 방법은 무작위 유니터리를 생성하여 평균적으로 이산화 오차 없이 시간 진화 연산자를 정확하게 재현합니다. 시스템은 단순한 시작 해밀토니안($H_Z$)의 바닥 상태에서 타겟 분자 해밀토니안($H$)으로 진화합니다. 효율성을 높이기 위해, 저자들은 보존량 항(입자 수 연산자)을 추가하여 해밀토니안의 상호작용 노름($\mu_I$)을 최적화하였으며, 회로 깊이를 줄이기 위해 패리티 기반 및 점유수 기반 연산자 축소를 포함한 고전적 전처리 기술을 적용했습니다.
2. 에너지 추정: 저자들은 노이즈에 강한 변형된 반복 양자 위상 추정(Iterative Quantum Phase Estimation, iQPE)을 활용했습니다. 기대값을 직접 샘플링하는 대신, 이 알고리즘은 Hadamard 테스트를 사용하여 중첩값 $\rho(s) = \text{Im} \langle \psi(T) | e^{is(E_{\text{trial}} - H)} | \psi(T) \rangle$의 허수부를 측정합니다. 두 개의 시도 에너지($E_{\text{trial}}^\pm$)를 선택하고 고정된 중심 시간 $s$에서 $\rho(s)$를 측정함으로써, 아크탄젠트 공식을 통해 GSE를 추출합니다. 이 접근 방식은 진폭 정보를 활용하며 탈분극 노이즈(depolarizing noise)에 강합니다.
3. 오류 완화 및 후처리: 실험에는 두 가지 고전적 사후 선택(post-selection) 기술이 사용되었습니다:
   • 패리티 사후 선택 (Parity Post-selection): 물리적 큐비트의 전체 패리티가 예상되는 섹터에서 벗어나 누설(leakage)이나 오류를 나타내는 경우 측정 결과를 폐기합니다.
   • 하트리-포크 (HF) 투영: 비-보조(non-ancilla) 큐비트가 초기 HF 상태로 돌아오지 않은 샷(shot)들에 대해 보조 큐비트의 기대값을 0으로 설정합니다. 이는 노이즈가 없는 환경에서 그러한 기여가 평균적으로 0이 된다는 이론적 기대에 근거합니다.

주요 결과

• 에너지 추정: 실험은 정확한 GSE로부터 $-25.5 \pm 5.3$ mHa의 편차를 갖는 최종 에너지 추정치를 산출했습니다. 이는 고전적 하트리-포크 추정치가 $52.8$ mHa의 편차를 보인 것보다 개선된 결과입니다. 양자 결과는 HF 에너지로부터 5 표준 편차 이상 통계적으로 구별됩니다.
• 노이즈 분석: 저자들은 서로 다른 노이즈 소스의 영향을 격리하기 위해 광범형 고전 시뮬레이션을 수행했습니다:
  • 결맞음 및 비결맞음 노이즈: 현실적인 결맞음(예: $Z$-회전) 및 비결맞음(탈분극) 오류율을 포함한 시뮬레이션은 화학적 정확도(약 $\sim 2$ mHa 이내)에 근접한 에너지 편차를 예측했습니다. 이는 알고리즘이 이러한 표준 노이즈 유형에 내성이 있음을 시사합니다.
  • 누설 오류 (Leakage Errors): 시뮬레이션에 누설 오류(계산 서브스페이스 밖으로의 전이)를 도입했을 때, 상당한 계통적 편향이 나타났습니다. 시뮬레이션에서의 편향 크기는 누설율이 포함되었을 때만 실험적 편차($\approx -25.5$ mHa)와 일치했습니다.
• 누설 징후: 측정 통계 분석 결과, 노이즈가 없는 시뮬레이션과 비교했을 때 하드웨어 데이터에서 물리적 큐비트와 보조 큐비트 모두에 대해 $|1\rangle$ 상태를 측정하는 경향이 강한 편향을 보였습니다. 이러한 불균형은 누설 이벤트, 특히 과도한 셔틀링(shuttling)과 게이팅을 거치는 보조 큐비트가 누설에 매우 취약하다는 점과 일치합니다.

의의 및 주장
본 논문은 확장 불가능한 작업을 고전 컴퓨터로 오프로딩하지 않는, 실용적이고 확장 가능한 양자 화학 알고리즘을 입증합니다. 주요 기여는 현재의 이온 트랩 하드웨어에서 화학적 정확도 달성을 가로막는 주된 장애물이 결맞음 또는 비결맞음 노이즈가 아니라 누설 오류임을 식별한 것입니다.

저자들은 누설 오류가 없다면, 자신들의 실험 설정과 노이즈 내성 iQPE 프로토콜이 샷 노이즈를 포함하더라도 화학적 정확도(2배 이내의 인자)에 근접한 결과를 얻었을 것이라고 주장합니다. 결과적으로, 본 연구는 이온 트랩 장치에서의 향후 양자 화학 응용의 발전이 단순히 표준 게이트 불확실성을 줄이는 것뿐만 아니라, 하드웨어 설계(예: 보조 큐비트를 위한 누설 탐지 코드) 및 알고리즘적 적응을 통해 누설 억제를 목표로 하는 데 달려 있음을 강조합니다.