Quantum linear solvers for quantum chemistry: prospects of exponential quantum advantage

기술 요약: 양자 화학을 위한 양자 선형 솔버: 지수적 양자 우위의 전망

문제 정의
Harrow-Hassidem-Lloyd(HHL) 알고리즘과 같은 양자 선형 솔버(QLS)는 선형 방정식 시스템을 해결하는 데 있어 이론적으로 지수적 양자 우위를 제공할 잠재력을 가지고 있다. 양자 화학에서 선형화된 결합 클러스터(LCC) 방정식은 선형 시스템($A\vec{x} = \vec{b}$)으로 재구성될 수 있으므로, QLS의 자연스러운 후보가 된다. 그러나 이러한 이점의 실현 여부는 행렬 $A$의 조건수($\kappa$)가 시스템 크기에 따라 어떻게 스케일링되는지에 달려 있다. 기존 문헌에 따르면 유리한(다항 로그, polylogarithmic) $\kappa$ 스케일링은 드물며, 많은 응용 분야에서 다항식 또는 지수적 성장을 보인다. 또한, 이전의 HHL을 이용한 양자 화학 적용은 주로 단일 참조 선형화 결합 클러스터(SRLCC) 방법에 국한되었으며, 이는 결합 끊어짐 및 형성을 설명하는 데 필수적인 강한 정적 상관관계를 포착하는 데 실패한다. 본 연구가 다루는 핵심 과제는 다음과 같다: (1) 강한 상관관계를 처리하기 위해 다중 참조 시스템으로 QLS 프레임워크를 확장하는 것, 그리고 (2) LCC 행렬의 조건수와 희소성(sparsity)이 켤리 Conjugate Gradient(CG) 알고리즘과 같은 고전적 방법 대비 지수적 양자를 보장할 만큼 유리하게 스케일링되는지를 엄밀하게 평가하는 것이다.

방법론
저자들은 이론적 유도, 스케일링 분석, 그리고 수치적 개념 증명 시뮬레이션을 결합한 종합적인 프레임워크를 제안한다:

1. 이론적 확장: 저자들은 내부 계약 다중 참조 선형화 결합 클러스터(icMRLCC) 방법에 기반하여 기존의 QLS-SRLCC 프레임워크를 다중 참조 형식으로 확장한다. 이들은 icMRLCCSD(단일 및 이중 들뜸)에 대한 작동 방정식을 유도하며, 진폭 방정식이 선형 시스템 $A\vec{x} = \vec{b}$로 공식화될 수 있음을 입증한다. 여기서 해 벡터 $\vec{x}$는 들뜸 진폭을 포함하며, 상관 에너지는 입력 상태 $|b\rangle$와 출력 상태 $|x\rangle$ 사이의 중첩 측정(Hong-Ou-Mandel 회로 사용)을 통해 추출된다.
2. 스케일링 분석: 양자 우위의 타당성을 평가하기 위해, 저자들은 시스템 크기($N$)에 따른 $A$ 행렬의 조건수($\kappa$)와 희소성($s$)의 스케일링을 분석한다. 이를 위해 세 가지 상호 보완적인 진단 도구를 사용한다:
   • 명시적 계산: 가상 스핀 궤도(virtual spin-orbitals)의 수를 변화시키며 소규모 모델 시스템(LiH, H4, BeH$_2$, BeH$_2^+$)에 대해 $\kappa$를 직접 계산한다.
   • 대각 비율 프록시(Diagonal Ratio Proxy): $A$의 최대 대각 성분과 최소 대각 성분의 비율($d_{max}/d_{min}$)을 사용하여 $\kappa$의 상한을 유도하며, 이는 대각 우세(diagonally dominant) 및 비대각 우세 사례 모두에서 검증된다.
   • 에지 생성 추측(Edge Spawning Conjecture): 그래프 인접 행렬에 관한 데이터 기반 추측을 채택한다. 비제로(non-zero) 행렬 요소를 그래프의 에지로 매핑함으로써, 행렬 패턴의 "확산성(diffuseness)"을 분석하여 $\kappa$ 스케일링 거동을 추론한다.
3. 실행 시간 및 자원 추정: 저자들은 고전적 CG와 양자 솔버(HHL 및 최적에 가까운 Childs-Kothari-Somma(CKS) 알고리즘)를 비교하는 실행 시간 비율($R$)을 도출한다. 또한, 두 큐비트 게이트 수를 추정하고, 고전적 전처리 비용(행렬 $A$와 $\vec{b}$의 구성) 및 후처리 비용에 대해 논의한다.
4. 수치 시뮬레이션: Qiskit을 사용하여 HHL 알고리즘을 통한 LiH, H4, BeH$_2$의 퍼텐셜 에너지 곡선(PEC)에 대한 개념 증명 시뮬레이션을 수행한다. 이 시뮬레이션은 HHL-icMRLCCSD 결과를 고전적 방법(CASSCF, icMRCI, CCSD, CCSD(T)) 및 Full Configuration Interaction(FCI)과 벤치마킹한다.

주요 기여 및 결과

• 다중 참조 프레임워크: 본 논문은 icMRLCC 프레임워크를 통해 LCC 방정식을 QLS에 적합한 선형 시스템으로 성공적으로 공식화하였으며, 이를 통해 기존 단일 참조 방식의 한계를 극복하고 강한 상관관계 처리를 가능하게 하였다.
• 조건수 스케일링: 세 가지 진단 접근 방식 모두 시스템 크기(특히 가상 궤도의 수)에 따라 조건수 $\kappa$가 다항 로그(polylogarithmic) 방식으로 스케일링된다는 일관된 예측을 보여준다. 명시적 계산 결과는 다항 로그 트렌드와 일치하는 "계단식(staircase-like)" 성장을 보여주며, $d_{max}/d_{min}$ 비율이 이 거동을 밀접하게 추적한다.
• 희소성 및 우위: LCC 행렬의 희소성 $s$는 시스템 크기에 따라 하위 선형(sub-linear)으로 스케일링됨($s \sim N^{2/E}$, 여기서 $E$는 들뜸 랭크)을 확인하였다.
• 양자 우위 전망: 다항 로그 $\kappa$ 스케일링과 하위 선형 희소성 스케일링을 결합하여, 저자들은 다음과 같은 실행 시간 비율을 도출한다:
  • HHL 알고리즘의 경우, 들뜸 랭크 $E > 2$ (예: LCCSDT 이상)에서 지수적 우위가 실현된다.
  • 최적에 가까운 CKS 알고리즘의 경우, LCCSD 수준($E=2$)에서도 지수적 우위 달성이 가능하다.
• 수치적 정확도: 모델 시스템(최대 4개 원자)에 대한 개념 증명 시뮬레이션 결과, 벤치마크 고전 방법 대비 오차가 0.009%를 초과하지 않는 범위 내에서 바닥 상태 에너지를 회복하였다. HHL-icMRLCCSD 방법은 강한 상관관계 영역(예: H4 및 BeH$_2$의 늘어난 결합)에서 안정성을 보이며, 이 영역에서 단일 참조 방법보다 뛰어난 성능을 나타낸다.
• 자원 분석: 본 연구는 상태 준비 및 HHL 유니터리의 게이트 수를 상세히 기술한다. 양자 알고리즘 단계가 잠재적 가속을 제공하지만, 고전적 전처리 비용(적분 및 $A$ 행렬 구성)이 전체 워크플로우에서 여전히 중요한 병목 현상으로 남아 있음을 지적한다.

의의 및 주장
본 논문은 icMRLCC 프레임워크를 통해 강한 상관관계를 가진 양자 화학 문제에 대한 QLS의 첫 번째 일반적 처리를 제공한다고 주장한다. 본 연구의 주요 의의는 양자 우위를 결정하는 데 있어 결정적이지만 종종 간과되는 요소인 조건수 스케일링을 엄밀하게 평가했다는 점에 있다. LCC 문제에서 $\kappa$가 유리하게(다항 로그 방식으로) 스케일링됨을 입증함으로써, 본 연구는 특히 CKS와 같은 최적에 가까운 솔버를 사용할 때 지수적 양자 우위의 가능성을 뒷받침한다.

저자들은 QLS-LCC를 고정밀 활성 공간 방법인 양자 위상 추정(QPE)의 대체재가 아닌 보완재로 위치시킨다. 이들은 QPE가 화학적으로 유의미한 활성 공간(정적 상관관계를 포착)을 처리하고, QLS-LCC가 더 큰 외부 궤도 공간(동적 상관관계를 포착)을 로그 단위의 큐비트 요구량으로 처리하는 하이브리드 워크플로우를 제안한다. 결론적으로, 양자 알고리즘 단계는 유망하지만, 엔드 투 엔드(end-to-end)의 우위를 실현하기 위해서는 현재의 고전적 병목 현상을 대체할 효율적인 양자 전처리 기술의 개발이 필요하다고 밝힌다.