보조 큐비트 플뢰켓으로 양자 알고리즘 가속

초록

본 논문은 보조 큐비트를 이용한 주기적 플뢰켓(키크) 연산을 도입해, 전통적인 양자 어닐링의 최소 에너지 갭 제한을 극복하고 시뮬레이션 시간을 절반으로 단축시키면서 정확도는 향상시키는 방법을 제시한다. 짧은 거리와 무한 장거리의 횡자장 이징 모델, 그리고 양자 화학에서의 수소 분자 모델을 대상으로 수치 시뮬레이션과 시간 평균 해밀토니안 이론을 통해 효과를 검증하였다.

상세 분석

이 연구는 양자 어닐링이 시스템 규모가 커질수록 최소 여기갭(Δ_min)이 급격히 감소하면서 어닐링 속도를 제한한다는 근본적인 문제에 주목한다. 저자들은 이 한계를 극복하기 위해 보조(ancillary) 큐비트 집합 A와 주 시스템 큐비트 집합 S를 결합한 확대 힐베르트 공간 S⊗A를 정의하고, 항상 작동하는 시스템 유니터리 U_S(t)와 동시에 작동하는 엔탱글링 플뢰켓 유니터리 U_SA(Θ̂_A, Ŝ)=exp(−i Θ̂_A⊗Ŝ)를 도입한다. 여기서 Θ̂_A는 보조 큐비트에 적용되는 단일 파라미터 θ를 갖는 Pauli 연산의 선형 결합이며, Ŝ는 시스템 전체에 걸친 집합 연산 Σ_{l=1}^{N_S}σ_S^l 로 정의된다. 핵심 아이디어는 Θ̂_A와 Ŝ의 선택을 통해 보조 자유도(게이지 DOF)를 ‘조정’함으로써, 시스템 서브스페이스에서의 에너지 풍경을 인위적으로 평탄화하거나 원하는 방향으로 변형시켜, 비아디아틱 전이를 억제하고 목표 바닥 상태에 더 빠르게 도달하도록 하는 것이다.

플뢰켓은 디지털-아날로그 혼합 모델에서 구현되며, δ(t−KΔt_K) 형태의 디랙 임펄스가 일정 간격으로 시스템 해밀토니안에 삽입된다. 이때 최적의 킥 각 θ_opt은 시간 평균 해밀토니안(마그누스 전개 1차 항) 분석을 통해 도출되며, 식 (2)에서 제시된 바와 같이 θ_opt≈(1/N_A N_K)·(1−E_T^S/E_mix^S(τ)) 로 근사된다. 여기서 E_T^S는 실제 바닥 상태 에너지, E_mix^S(τ)는 믹서 해밀토니안에 의한 평균 에너지이다. 이 식은 시스템 해밀토니안의 상세 구조에 크게 의존하지 않으며, 믹서와 킥 해밀토니안이 지배적인 경우에 보편적으로 적용 가능함을 저자들은 수치 실험으로 확인한다.

구체적인 실험으로는 (1) 개방 사슬의 근접 이웃 TFIM(NN‑TFIM), (2) 완전 연결 장거리 TFIM(ILR‑TFIM), (3) 브라위‑키타에프 변환 후 두 스핀으로 인코딩된 수소 분자 모델을 선택하였다. NN‑TFIM에서는 θ_0^XX와 θ_0^Z를 각각 문제와 믹서 강도로 설정하고, 초기 큰 횡자장(θ_0^Z) 하에서 작은 θ_opt≈0.004의 플뢰켓을 적용하면, Δ_min보다 빠른 어닐링 속도에서도 시스템 에너지가 바닥 상태에 빠르게 수렴한다는 것을 보여준다. 킥 강도가 과도하면 비아디아틱 전이가 발생해 최종 에너지가 실제 바닥보다 높아지지만, 적절히 조절된 θ_opt은 이러한 과잉을 방지한다. 또한 킥 간격 Δt_K를 충분히 짧게(θ_0^Z Δt_K≪1, θ_0^XX Δt_K≪1) 유지하면 고주파 킥이 지속적으로 에너지 풍경을 보정해 주어, 연속적인 약한 킥(θ→0, N_K→∞) 경우에도 비슷한 가속 효과를 얻을 수 있다.

ILR‑TFIM에서는 연결도가 N_S(N_S−1)/2 로 급증함에도 불구하고, 동일한 θ_opt와 킥 스케줄을 적용했을 때 성능 저하가 거의 없으며, 이는 제안된 프로토콜이 시스템 규모와 연결성에 대해 강인함을 의미한다. 마지막으로 수소 분자 모델에서는 파리티 보존 특성 때문에 전통적인 Z‑기반 믹서가 작동하지 않으며, X‑기반 믹서를 사용하고, 보조 연산자를 Z_A 로 선택한 뒤 동일한 플뢰켓을 적용하였다. 결과적으로 바닥 상태 에너지(≈0 a.u.)를 정확히 재현하면서도 어닐링 시간은 절반 이하로 단축되었다.

전체적으로 이 논문은 (i) 보조 큐비트를 이용한 단일 파라미터 플뢰켓 설계, (ii) 시간 평균 해밀토니안을 통한 최적 각도 추정, (iii) 다양한 모델에 대한 수치 검증이라는 세 축을 통해, 디지털‑아날로그 양자 시뮬레이션에서 비아디아틱 효과를 제어하고 어닐링 속도를 크게 향상시킬 수 있음을 입증한다. 특히 최적 킥 각도가 시스템 크기에 거의 독립적이라는 점은 실험적 구현 시 스케일업에 큰 장점을 제공한다.