양자 자코비‑데이비드슨 방법: 빠른 수렴을 위한 새로운 하이브리드 알고리즘

초록

**
본 논문은 기존 양자 데이비드슨(QD) 방법보다 빠른 수렴을 보이는 양자 자코비‑데이비드슨(QJD)과 샘플 기반 변형(SBQJD)을 제안한다. LCU와 뉴턴‑유사 교정 방정식을 활용해 서브스페이스를 효율적으로 확장하고, 8‑큐빗 대각우세 행렬, 12‑큐빗 1차원 이징 모델, 10‑큐빗 물(H₂O) 해밀토니안에 대한 정확한 시뮬레이션을 통해 측정 횟수와 반복 횟수에서 QD 대비 현저히 적은 자원을 요구함을 입증한다.

**

상세 분석

**
본 연구는 고전적인 자코비‑데이비드슨(JD) 알고리즘의 핵심 아이디어를 양자 컴퓨팅에 그대로 옮겨온 점이 가장 큰 혁신이다. JD는 잔차 벡터와 직교하는 교정 벡터 |t⟩ 를 구해 서브스페이스를 확장하는데, 이때 교정 방정식은 (I‑|rv⟩⟨rv|)(H‑E′I)(I‑|rv⟩⟨rv|)|t⟩=‑|r⟩ 로 정의된다. 논문은 이 방정식을 양자 회로 수준에서 구현하기 위해 두 가지 형태를 제시한다. 첫 번째는 정확한 역연산 (H‑E′I)⁻¹ 을 LCU(Linear Combination of Unitaries) 기법으로 근사하는 방법이며, 두 번째는 대각 원소만을 이용한 전처리 M≈Diag(H)‑E′I 를 이용해 M⁻¹ 을 적용하는 프리컨디셔닝이다. LCU 회로는 ancilla 비트를 ⌈log₂m⌉개 사용해 A=∑α_iU_i 를 구현하고, 성공 확률은 1/s² ( s=∑|α_i| ) 로 제어된다. 여기서 amplitude amplification을 적용하면 전체 복잡도가 O(s·polylog(s/p)) 으로 감소한다.

또한 저자는 교정 연산을 뉴턴 스텝으로 해석함으로써, 수렴 속도가 1차 방법인 QD의 선형 수렴을 넘어 최소 2차 수렴을 보장한다는 이론적 근거를 제시한다. 이는 잔차 벡터만을 이용하는 QD와 달리, 교정 벡터가 실제 에너지 표면의 기울기를 반영하므로 작은 잔차에서도 급격히 오차가 감소한다는 점에서 실용적이다.

샘플 기반 양자 대각화(SQDiag)와의 결합(SBQJD)은 또 다른 중요한 기여이다. SQDiag은 측정 결과에서 가장 빈번히 나타나는 베이스 상태들을 샘플링해 서브스페이스를 구성하고, 그 서브스페이스 내에서 고전적인 대각화를 수행한다. 이를 QJD와 연계하면, 초기 레퍼런스 상태가 충분히 좋은 경우 교정 벡터를 매우 적은 수의 샘플로도 정확히 추정할 수 있어 전체 회로 깊이가 크게 감소한다.

실험적으로는 8‑큐빗 대각우세 행렬, 12‑큐빗 1차원 이징 모델, 10‑큐빗 물 분자 해밀토니안을 대상으로 정확한 시뮬레이션을 수행했다. 모든 사례에서 QJD와 SBQJD는 QD 대비 평균 3‑5배 적은 반복 횟수와 2‑4배 적은 파울리 측정 샷을 요구했으며, 최종 에너지 오차는 화학 정확도(1 mHa) 이하로 수렴했다. 특히 대각우세 행렬에서는 프리컨디셔닝 M⁻¹이 거의 정확한 역연산을 제공하므로, 교정 벡터 계산 비용이 크게 낮아졌다.

자원 분석 측면에서 QJD는 각 반복마다 역연산을 포함하므로 QPU 사이클이 QD보다 무겁지만, 전체 반복 수가 크게 감소해 총 실행 시간은 오히려 짧다. 또한, fault‑tolerant 양자 컴퓨터가 보편화될 경우, 역연산을 위한 양자 신호 처리(예: Hamiltonian simulation + QPE)와 결합해 실용적인 고정밀 전자구조 계산에 바로 적용 가능할 것으로 기대된다.

**