양자 메트로폴리스 헤이스팅스 알고리즘 회로 설계

초록

본 논문은 메트로폴리스‑헤이스팅스(MH) 알고리즘의 양자화에 대해, 전통적인 전이 확률을 코히런트하게 계산하는 대신 제안‑수락 로직을 그대로 이용하는 새로운 Szegedy 양자 워크 구조를 제시한다. 제안된 회로는 고정된 수의 제안·수락 오라클 호출만으로 구현 가능하며, 기존 방법에 비해 논리 큐 수를 크게 절감한다. 또한 이중(dual) 커널을 이용한 양자 워크가 원래 MH 커널과 동일한 스펙트럼 갭을 유지함을 증명해, 이론적인 2차 속도 향상이 실제 시뮬레이션에서도 유지됨을 보인다.

상세 분석

이 논문은 Szegedy의 양자 워크 프레임워크를 메트로폴리스‑헤이스팅스(MH) 마코프 체인에 적용하면서, 기존의 “전이 확률을 코히런트하게 계산”하는 접근법이 요구하는 높은 양자 연산 비용과 큐 수 문제를 근본적으로 해결한다는 점에서 의미가 크다. 전통적인 Szegedy 양자 워크는 마코프 커널 (P) 를 직접 양자화해 (|x\rangle|y\rangle) 형태의 이중 레지스터에 전이 확률 (P(x,y)) 를 부호화한다. 이 과정에서 각 전이 확률을 정밀히 계산해야 하며, 이를 위해 가역 연산을 수행하는 회로가 필요하고, 그 회로는 전이 확률의 복잡도에 비례해 보조 큐를 다수 필요로 한다. 특히 MH 알고리즘은 제안 단계와 수락 단계로 구성되며, 거부 확률을 직접 계산해야 하는 “샘플 포겟팅” 단계가 비가역적이라 양자화가 까다롭다.

저자들은 이러한 문제를 회피하기 위해 두 가지 핵심 아이디어를 도입한다. 첫 번째는 “이중 커널(dual kernel)” 개념이다. 원래의 상태 공간 (S) 대신 상태 쌍 ((x,y)) 로 이루어진 에지 공간을 마코프 체인의 새로운 상태 공간으로 정의한다. 이 공간에서 한 단계는 (i) 현재 에지의 헤드(목표 상태)를 새로운 제안 상태로 교체하고, (ii) 현재 에지를 뒤집을지 말지를 수락 확률에 따라 결정한다. 중요한 점은 이 과정에서 “거부”가 존재하지 않으며, 모든 연산이 유니터리하게 구현될 수 있다는 것이다. 따라서 전이 확률을 직접 코히런트하게 계산할 필요가 없고, 제안 오라클 (O_T)와 수락 오라클 (O_A)만을 일정 횟수 호출하면 된다.

두 번째는 회로 설계 최적화이다. 저자들은 4개의 상태 레지스터와 3개의 보조 레지스터만을 사용해 전체 양자 워크 연산을 구현한다. 상태 레지스터는 각각 (\lceil\log_2 n\rceil) 비트를 차지하므로, 전체 큐 수는 (4\lceil\log_2 n\rceil+3) 로, 기존 방법이 요구하던 (O(\text{poly}(n))) 수준의 큐보다 훨씬 적다. 특히 메트로폴리스‑조정 랭게뱅 알고리즘(MALA) 같은 연속 공간 샘플링에 대해서도, 좌표를 이산화해 동일한 구조를 적용함으로써 선형적인 큐 스케일링을 달성한다.

스펙트럼 분석 측면에서도 중요한 결과를 제시한다. 이중 커널의 전이 행렬은 원래 MH 커널과 동일한 스펙트럼 갭 (\delta) 를 갖는다. Szegedy 양자 워크의 일반적인 특성에 따라, 양자 워크 연산자의 스펙트럼 갭은 (\Omega(\sqrt{\delta})) 로 확대된다. 저자들은 실제 회로 시뮬레이션을 통해, MALA에 대한 양자 워크 연산자 (W) 가 이론적 최소 갭 (\cos^{-1}(\sqrt{1-\delta/2})) 를 초과하는 갭을 보임을 확인했다. 이는 양자 워크가 원래 마코프 체인의 믹싱 타임을 제곱근 수준으로 가속할 수 있음을 실증한다.

또한 기존 연구와의 비교를 통해, 제안된 방법이 Childs 등(2020)과 Lemieux 등(2021)의 회로보다 큐 효율성이 현저히 우수함을 표로 정리한다. 특히 Childs 방법은 각 좌표를 (d) 비트로 이산화하고, 여러 보조 레지스터를 추가로 필요로 하여 전체 큐 수가 (O(pd)) 수준이다. Lemieux 방법은 무관계형 스핀 시스템에 특화된 유니터리 인코딩을 사용하지만, 여전히 (N) 비트 규모의 무브 레지스터가 필요하다. 반면 본 논문의 회로는 이진 인코딩을 사용해 (m=\lceil\log_2 n\rceil) 비트만으로 모든 상태를 표현하고, 고정된 3개의 보조 레지스터만 추가한다.

결과적으로, 이 논문은 “전이 확률을 직접 코히런트하게 계산하지 않아도 되는” 양자 MH 워크 구현 방식을 제시함으로써, 근본적인 큐 제한을 완화하고, 근접한 실제 양자 하드웨어에서도 적용 가능한 설계를 제공한다. 이는 양자 가속을 기대하는 MCMC 시뮬레이션, 특히 분자 동역학, 베이즈 추정, 머신러닝 분야에 큰 파급 효과를 가질 것으로 기대된다.