Impact of molecular orbital localization on quantum computational resources for Hamiltonian simulation: A benchmark study of hydrogen chain systems

기술 요약: 해밀토니안 시뮬레이션을 위한 양자 컴퓨팅 자원에 대한 분자 궤도 국소화의 영향

문제 정의
양자 시뮬레이션 알고리즘(예: Trotterization, qubitization) 및 압축 전략(예: tensor hypercontraction, QCELS)을 최적화하는 데는 상당한 연구가 집중되어 왔으나, 분자 궤도 기저(basis)의 선택이 양자 컴퓨팅 자원에 어떻게 영향을 미치는지에 대해서는 상대적으로 관심이 적었습니다. 전통적인 양자 화학에서 하트리-포크(Hartree–Fock, HF) 이론으로부터 유도된 정규 분자 궤도(Canonical Molecular Orbitals, CMOs)는 비국소화(delocalized)되어 있어, 2차 양자화된 해밀토니안에서 "all-to-all" 결합을 초래합니다. 이러한 비국소화는 시뮬레이션을 위해 방대한 수의 양자 게이트를 필요로 합니다. 반면, 국소 분자 궤도(Localized Molecular Orbitals, LMOs)는 공간적 국소성을 활용하기 위해 고전적 전자 구조 방법론에서 널리 사용되며, 여기서 멀리 떨어진 궤도 간의 상호작용은 무시할 수 있는 수준입니다. 본 연구가 다루는 핵심 질문은 궤도 국소화가 트로터 분해(Trotter decomposition) 기반 접근 방식에 대해 이와 유사하게 복잡도와 자원 요구량을 줄일 수 있는지 여부입니다.

방법론
저자들은 1.0 Å의 원자 간 거리를 가진 $n$개의 원자(8에서 100개 사이)로 구성된 1차원 선형 수소 사슬 시스템($H_n$)을 사용하여 벤치마크 연구를 수행했습니다. 연구에서는 두 가지 궤도 표현을 비교했습니다:

1. CMOs: 표준 하트리-포크 정규 궤도.
2. LMOs: 상대적 공간 위치에 따라 재정렬된 파이프-메제이(Pipek–Mezey) 국소 궤도.

해밀토니안은 조던-위그너(Jordan–Wigner, JW) 변환을 사용하여 구축되었습니다. 저자들은 항의 개수를 줄이기 위해 두 가지 서로 다른 해밀토니안 절단(truncation) 전략을 조사했습니다:

• 연산자 계수 기반 절단 ($c_{thre}$): 절대값이 임계값 미만인 항을 제거함.
• 연산자 국소성 기반 절단 ($k_{max}$): 연산자 국소성 $k$(파울리 문자열 내의 identity가 아닌 파울리 연산자 수)가 임계값을 초과하는 항을 제거함. 2-전자 항의 경우, 스핀 편향을 피하기 위해 하나의 $\alpha$ 전자와 하나의 $\beta$ 전자의 여기를 기준으로 국소성을 정의했습니다.

정확도는 미절단 참조 계산(작은 사슬의 경우 Full-CI, 큰 사슬의 경우 DMRG)에 대한 바닥 상태 에너지 오차($\Delta E$) 및 충실도($F$)를 통해 평가되었습니다. 높은 충실도($F \geq 0.99$)를 달ach하기 위해 필요한 임계값을 추정하기 위해, 절단 임계값과 충실도 사이의 관계에 대한 경험적 공식이 피팅되었습니다. 이렇게 추정된 임계값들은 $H_{100}$까지의 사슬에 대한 양자 게이트 수(Clifford, CNOT, 그리고 non-Clifford $R_z$ 게이트)를 투영하는 데 사용되었습니다.

주요 기여 및 결과

1. 기저별 최적 절단 전략:

   • CMO 기반 전개에서는 계수 기반 절단($c_{thre}$)이 효과적입니다. 느슨한 임계값에서도 바닥 상태 에너지 오차와 충실도가 이상적인 값에 가깝게 유지되는데, 이는 비국소화된 CMO의 특성상 파동 함수를 기술하기 위해 고국소성 항들이 필수적이기 때문입니다.
   • LMO 기반 전개에서는 국소성 기반 절단($k_{max}$)이 더 우수합니다. 계수 기반 절단은 LMO에 대해 즉각적으로 실패하며 정확도를 급격히 저하시킵니다. 그러나 국소성 기반 절단은 낮은 $k_{max}$ 값에서도 높은 충실도를 유지하는데, 이는 공간적으로 분리된 LMO 기반의 궤도들이 약한 상관관계를 갖기 때문입니다.

2. 절단 임계값의 스케일링:

   • $F \geq 0.99$를 유지하기 위해, CMO의 계수 임계값은 $c_{thre} \propto n^{-2.07}$로 스케일링됩니다.
   • LMO의 경우, 요구되는 국소성 임계값 $k_{max}$는 시스템 크기에 따라 로그-다항식(polylogarithmically)으로 스케일링됩니다, 구체적으로 $k_{max} \propto (\log n)^{1.86}$입니다.

3. 양자 게이트 자원 추정:

   • CMO 기저: 해밀토니안 시뮬레이션에 필요한 양자 게이트 수는 시스템 크기에 따라 다항식적으로 성장하며, 구체적으로 $O(n^4)$에서 $O(n^5)$입니다. 절단은 계수를 줄이지만, 다항식 스케일링 차수는 변화시키지 못합니다.
   • LMO 기저: 국소성 기반 절단과 결합될 때, 게이트 수는 로그-다항식적으로 성장합니다, 구체적으로 $O((\log n)^{12})$에서 $O((\log n)^{14})$입니다.
   • CNOT 게이트: CMO 기저에서는 절단 후에도 많은 고국소성 항들이 남아 있어 높은 CNOT 수를 초래합니다. LMO 기저에서는 절단이 해밀토니안을 저국소성 항들로 제한하여 CNOT 게이트의 성장을 억제합니다. $n \geq 50$인 경우, LMO 기반 접근 방식은 점군 대칭성(point-group symmetry)의 파괴로 인해 초기에는 더 많은 해밀토니안 항을 가짐에도 불구하고, CMO 기반 접근 방식보다 적은 CNOT 게이트를 요구합니다.

4. 인코딩 비교:

   • 연구는 또한 조던-위그너(JW)와 브라비-키타예프(Bravyi–Kitaev, BK) 인코딩을 비교했습니다. BK는 CMO에 대한 CNOT 스케일링을 $O(n^5)$에서 $O(n^4 \log n)$으로 줄이지만, 결과적인 게이트 수는 LMO 기반 접근 방식과 유사한 수준입니다. 또한, BK는 파울리 문자열 내의 가변적인 Pauli-X/Y 항의 수로 인해 단일 큐비트 Clifford 게이트의 수를 증가시킵니다.

의의 및 주장
본 논문은 궤도 표현의 선택과 그에 따른 절단 전략이 대규모 양자 화학 시뮬레이션의 실현 가능성을 결정하는 결정적인 요소라고 주장합니다. 저자들은 트로터 기반 해밀토니안 시뮬레이션에서 LMO와 국소성 기반 절단을 함께 사용하는 것이 전통적인 CMO 기반 접근 방식에 비해 상당한 지수적 이점(exponential advantage)을 제공한다고 결론짓습니다. 특히, 게이트 요구량의 다항식 스케일링에서 로그-다항식 스케일링으로의 전환은 LMO가 미래의 결함 허용(fault-tolerant) 양자 컴퓨터에서 훨씬 더 큰 분자 시스템을 시뮬레이션할 수 있게 할 것임을 시사합니다.

저자들은 본 연구가 연산자 국소성이 공간적 거리와 직접적으로 상관관계가 있는 1D 수소 사슬에 집중되어 있지만, 그 결과가 기저 설계 전략에 유용한 가이드라인을 제공한다고 언급했습니다. 그들은 복잡한 3D 분자의 경우, 엄격한 국소성 기반 절단보다는 거리 기반 절단이 더 적절할 수 있다고 제안하지만, 궤도 국소성을 활용하는 일반적인 원칙은 계산 비용을 줄이기 위한 유망한 방안으로 남아 있습니다. 또한, 연구는 상태 준비(state preparation)가 잠재적인 병목 구간이 될 수 있음을 강조하며, 완전한 엔드 투 엔드 비용 평가를 위해서는 상태 준비에 대한 추가적인 자원 추정이 필요하다고 명시했습니다.