양자 HHL 알고리즘으로 선형 방정식 풀기: 학부생을 위한 물리·수학 튜토리얼

초록

본 논문은 Harrow‑Hassidim‑Lloyd(HHL) 알고리즘의 물리적·수학적 원리를 학부 수준에서 이해하도록 돕는 교육용 튜토리얼을 제공한다. 선형 시스템을 양자 상태로 인코딩하고, 상태 준비, 양자 위상 추정(QPE), 앙상블 회전, 역 QPE 순으로 구현하는 과정을 상세히 설명한다. Qiskit 기반 구현 예시와 고전 시뮬레이션 비교를 통해 복잡도, 제한점, 향후 발전 방향을 논의한다.

상세 분석

이 논문은 HHL 알고리즘을 단순히 소개하는 수준을 넘어, 물리학과 컴퓨터 과학 전공 학부생이 직접 회로를 설계하고 실행할 수 있도록 단계별로 분해한다. 먼저 선형 시스템 (A\mathbf{x}=\mathbf{b}) 을 양자 상태 (|b\rangle) 와 Hermitian 행렬 (A) 로 변환하고, (|b\rangle)를 정규화하는 과정이 강조된다. Hermitian 변환이 필요할 경우 (\begin{pmatrix}0&A\A^\dagger&0\end{pmatrix}) 와 같은 블록 행렬을 도입함으로써 일반 행렬도 적용 가능하도록 설계한다.

알고리즘의 핵심은 양자 위상 추정(QPE)이다. 여기서는 (U=e^{iAt}) 를 제어 연산으로 구현하고, (n)개의 레지스터 (q_c) 에 위상 (\theta=\lambda t/2\pi) 를 저장한다. 논문은 QPE의 세 단계—균등 초위상, 제어 유니터리, 역 푸리에 변환(QFT†)—를 구체적인 게이트 시퀀스로 제시하고, 각 단계가 어떻게 고유값 (\lambda_i) 를 디지털화하는지 수식적으로 증명한다.

그 다음 앙상블(ancilla) 회전 단계에서는 고유값의 역수 (1/\lambda_i) 를 확률 진폭으로 변환한다. 회전 각도 (\arcsin(C/\tilde\lambda_i)) 를 선택함으로써, 측정 시 (|1\rangle_a) 가 관측될 확률이 (\propto 1/\lambda_i^2) 가 되며, 이는 최종 해 (|x\rangle\propto A^{-1}|b\rangle) 와 직접 연결된다.

역 QPE 단계는 레지스터 (q_c) 에 남아 있는 위상 정보를 제거해 원래 (|b\rangle) 와 동일한 형태의 해를 얻는다. 논문은 이 전체 흐름을 Qiskit 코드와 회로 다이어그램으로 구현하고, 2‑qubit, 4‑qubit 예시를 통해 실제 양자 상태의 진화 과정을 시각화한다.

복잡도 분석에서는 이상적인 상황에서 (O(\log N)) 시간에 해를 얻을 수 있지만, 실제 NISQ 디바이스에서는 게이트 오류, 디코히런스, 고유값 스케일링 문제 등으로 인해 제한이 존재함을 강조한다. 특히 고유값이 0에 가까워질 경우 회전 각도가 급격히 커져 성공 확률이 급감하고, 이를 완화하기 위한 사전 스케일링 및 정규화 기법을 제시한다.

마지막으로, 현재 구현 가능한 문제 규모가 매우 제한적이며, 오류 정정 코드와 더 깊은 회로 최적화가 필요함을 지적한다. 향후 연구 방향으로는 변분형 HHL, 하이브리드 양자‑고전 알고리즘, 그리고 양자 머신러닝과의 연계가 제시된다.