주기적 대각 구조를 가진 희소 행렬의 블록 인코딩

초록

본 논문은 주기적인 대각 성분을 갖는 희소(밴드) 행렬을 효율적으로 블록 인코딩하는 양자 회로를 제시한다. 선형 결합 유니터리(LCU) 프레임워크와 복소 지수 연산을 실수·허수 성분으로 분리하는 특수 유니터리를 활용해, 일반적인 밀집 행렬 대비 다항식 규모의 게이트 복잡도를 달성한다. 대각 행렬의 경우 O(n)·희소 밴드 행렬의 경우 O(poly(n))의 복잡도를 보이며, ADR(Advection‑Diffusion‑Reaction) 방정식 등 편미분 방정식 해결에 적용 가능함을 수치 실험으로 검증한다.

상세 분석

이 논문은 양자 알고리즘에서 핵심적인 서브루틴인 블록 인코딩을, 주기적인 대각 구조를 가진 희소 행렬에 특화시킨 방법론을 제시한다. 핵심 아이디어는 두 단계로 나뉜다. 첫 번째는 주기적인 위상 ω를 갖는 대각 유니터리 V(ω)=diag(1,e^{iω},e^{i2ω},…,e^{i(N‑1)ω})를 구현하는 것으로, 이는 각 비트에 대해 Z축 회전(P(θ)=diag(1,e^{iθ}))을 적용함으로써 O(n)개의 단일‑큐빗 위상 게이트만으로 구현 가능하다. 이 회로는 Hadamard 게이트와 결합되어 실수 행렬 C(ω)=diag(cos(kω))와 허수 행렬 S(ω)=diag(sin(kω))를 각각 ancilla 큐빗의 |0⟩, |1⟩ 상태에 매핑한다. 정리 1에 의해 U_C(ω)=(H⊗I)(I⊗V†)(|0⟩⟨0|⊗V(2ω)+|1⟩⟨1|⊗I)(H⊗I) 가 정확히 C(ω)의 블록 인코딩을 제공하며, 스케일링 팩터 α=1, ancilla 수 m=1이다.

두 번째 단계는 LCU 프레임워크를 이용해 주대각 성분이 주기적이고, 오프‑대각이 토플리츠 형태인 밴드 행렬 M을 표현한다. M은 α₀C(ω)+α₁L+α₂R+α₃I 로 분해되며, 여기서 L,R는 순환 이동 퍼뮤테이션이며 각각 O(N) 복잡도의 토플리(또는 제어 토플리) 회로로 구현된다. 상태 준비 연산 P_REP은 α_j 를 정규화된 양자 상태 |j⟩에 매핑하고, SELECT 연산은 |j⟩⟨j|⊗U_j 형태로 제어된 유니터리를 적용한다. 전체 블록 인코딩 회로는 ancilla 3개(대각 블록, LCU 가중치, 선택 제어)와 O(n²) 게이트 복잡도를 가진다. 특히 V(ω)의 제어 버전은 추가 토플리와 n개의 제어 위상 회전으로 구현 가능하며, L,R의 제어 버전은 기존 양자 가산기 설계를 활용해 O(n) 게이트로 구현한다.

복잡도 분석에서 저자들은 최악의 경우(밴드 행렬)에도 O(poly(n)) 수준을 유지하고, 순수 대각 행렬에 대해서는 O(n)으로 선형 스케일링한다는 점을 강조한다. 이는 일반적인 밀집 행렬에 대한 블록 인코딩이 요구하는 O(2^n) 혹은 O(poly(2^n)) 복잡도와 비교해 현저히 효율적이다. 또한, ω가 작을 경우 V(ω)를 근사해 상위 O(log(1/ε)) 비트만 사용함으로써 추가적인 게이트 절감이 가능함을 제시한다.

응용 측면에서는 ADR(Advection‑Diffusion‑Reaction) 방정식의 주기적 계수를 가진 연산자를 양자화하여, QSVT(Quantum Singular Value Transformation)와 같은 최적 스케일링 알고리즘에 바로 삽입할 수 있음을 보였다. 수치 실험에서는 4‑qubit 시스템을 대상으로 확률 p₀가 ω와 초기 상태에 따라 어떻게 진동하는지를 확인하고, 제시된 회로가 기대한 블록 인코딩을 정확히 수행함을 검증하였다.

전체적으로 이 논문은 (1) 주기적 대각 구조를 가진 행렬에 대한 특수한 유니터리 구현, (2) LCU 기반의 효율적인 가중치·선택 회로 설계, (3) 복잡도 절감과 실제 물리‑공학 문제 적용 가능성을 동시에 제공함으로써, 양자 선형대수와 양자 시뮬레이션 분야에 실용적인 도구를 추가한다는 점에서 큰 의의를 가진다.