양자 컴퓨터로 구현한 다중체 혼돈의 동역학 시뮬레이션

초록

본 논문은 91개의 초전도 트랜스몬 큐비트를 이용해 최대 혼돈성을 갖는 듀얼 유니터리(dual‑unitary) 회로를 구현하고, 텐서‑네트워크 기반 오류 완화(TEM) 기법으로 무한 온도 자동 상관함수를 정확히 재현한다. 듀얼 유니터리 조건을 벗어난 비정칙 회로에 대해서는 고전 텐서‑네트워크 시뮬레이션과 비교해 양자 장치의 신뢰성을 검증한다.

상세 분석

이 연구는 디지털 양자 시뮬레이션이 연속시간 해밀토니안 분해에 의존하지 않고, 자체적으로 이산시간 동역학을 구현할 수 있다는 점에 착안한다. 특히, 공간과 시간 양쪽 모두에서 유니터리성을 유지하는 듀얼 유니터리 회로는 특정 상관함수와 OTOC(Out‑of‑Time‑Order Correlator)의 정확한 해를 제공한다는 이론적 특성을 활용한다. 저자들은 이론적으로 풀 수 있는 킥드 이징 모델을 선택하고, J와 b를 π/4로 고정해 듀얼 유니터리 포인트를 만들었다. h 파라미터를 변화시켜 비정칙(비‑클리포드) 회로로 확장함으로써, 양자 혼돈의 전이와 엔탱글먼트 성장률을 실험적으로 탐색한다.

양자 하드웨어는 IBM Quantum의 91‑큐비트 초전도 프로세서이며, 두‑큐비트 게이트는 ECR(Excitation‑Cross‑Resonance)으로 구현된다. 각 레이어별로 파울리 트와이링을 적용해 효과적인 파울리 채널을 추정하고, 기존 노이즈 학습 기법에 듀얼 유니터리 회로의 클리포드 포인트(h=0)를 이용해 모델의 자유도를 추가 보정한다. 이렇게 구축된 노이즈 모델을 텐서‑네트워크 형태로 역변환(M′)하여, 측정된 연산자 O를 변형된 연산자 O′=M′†(O)로 바꾸어 노이즈가 제거된 기대값을 추정한다. 이 과정은 정보완전(Informationally Complete) 측정을 위해 각 큐비트를 X, Y, Z 기저 중 무작위로 선택하고, 신호 큐비트는 X에 80 % 비중을 두는 비대칭 샘플링을 적용한다.

실험 결과는 듀얼 유니터리 포인트에서 무한 온도 자동 상관함수 C_n(t)=cos(2h)^t (t=n) 를 정확히 재현한다. 비‑클리포드 파라미터 h>0에서도 TEM을 적용하면 이론적 지수 감쇠율과 거의 일치한다. 시스템 규모가 51, 71, 91 큐비트로 증가함에 따라 회로 깊이도 늘어나지만, 오차 완화 후에도 평균 오차는 10⁻² 이하로 유지된다. 다만, 깊이가 깊어질수록 노이즈 모델의 누적 오차가 미세하게 나타나며, 이는 파울리 채널 추정의 불완전성, 시간 의존적 드리프트, 혹은 레이어 간 상관성 무시 등에 기인한다.

듀얼 유니터리 조건을 벗어난 b≠π/4 상황에서는 정확한 해가 존재하지 않으므로, 저자들은 고전 텐서‑네트워크 시뮬레이션(Heisenberg 및 Schrödinger 그림)과 비교한다. Heisenberg 방식(χ=500)과 실험값은 전반적으로 좋은 일치를 보였으며, Schrödinger 방식(χ=1500)은 큰 차이를 보였다. 이는 Heisenberg 텐서‑네트워크가 연산자 전파에 더 적합하고, 양자 회로의 광범위한 얽힘 구조를 효율적으로 포착하기 때문으로 해석된다.

이 논문은 (1) 듀얼 유니터리 회로가 디지털 양자 시뮬레이션의 벤치마크로 활용될 수 있음을, (2) 대규모 노이즈 모델링과 TEM이 프리‑포트볼트 양자 프로세서에서 비‑클리포드 동역학을 신뢰성 있게 측정하도록 만든다는 점을, (3) 고전 텐서‑네트워크와의 정량적 비교를 통해 양자 장치의 실용적 한계를 명확히 제시한다는 점에서 의미가 크다. 향후 연구는 더 복잡한 비평형 현상(예: 측정‑유도 전이, 시간결정성 파동)과 다중 파라미터 최적화를 결합해, 오류 완화와 하드웨어 설계가 공동으로 작동하는 양자 시뮬레이션 프레임워크를 구축하는 방향으로 나아갈 수 있다.