Spin SYK-dot chains: entanglement bound and resource estimates towards a quantum advantage

기술 요약: 스핀 SYK-도트 체인: 양자 우위를 향한 얽힘 경계 및 자원 추정

문제 정의
무작위적인 전방향(all-to-all) 4체 상호작용을 특징으로 하는 사크데브-예-키타예프(SYK) 모델은 강상관 양자 물질, 최대 혼돈(maximal chaos), 그리고 홀로그래피 이중성을 나타내는 전형적인 모델이다. 그러나 그 기저 상태는 부피 법칙(volume law)의 얽힘을 따르기 때문에, 텐서 네트워크를 이용한 고전적 시뮬레이션 시 지수적으로 높은 비용($\chi \sim 2^{N/2}$)이 발생하며, 이러한 얽힘을 생성하기 위해 필요한 회로 깊이로 인해 양자 하드웨어에서 준비하는 데에도 어려움이 있다. 비용을 완화하기 위해 희소(sparse) 변형 모델들이 탐구되어 왔으나, 이들은 종종 부피 법칙 얽힘을 유지하거나 제한된 연결성을 가진 하드웨어에서 과도한 SWAP 게이트 오버헤드를 겪는다. 핵심 과제는 오류 수정 없이도 SYK 유형의 모델에 대해 양자 하드웨어가 고전적 시뮬레이션을 능가할 수 있는 영역을 식별하는 동시에, 강력한 고전적 기준선(baseline)을 유지할 수 있을 만큼 모델이 다루기 쉬운 상태를 찾는 것이다.

방법론
저자들은 구조화된 희소화 SYK 모델인 "스핀 SYK-도트 체인"을 제안하고 분석한다. 이 시스템은 밀집된 무작위 전방향 4-바디 파울리 상호작용(SYK 도트)을 포함하는 $p$개의 큐비트를 가진 $q$개의 "도트(dots)"로 구성된다. 이 도트들은 경계 큐비트 사이의 강도 $K$를 가진 단일 비무작위 $ZZ$ 결합에 의해 체인 형태로 연결된다.

본 연구는 두 가지 접근 방식을 병행한다:

1. 고전적 시뮬레이션: 저자들은 허수 시간 행렬 곱 상태(imaginary-time Matrix Product State, MPS) 알고리즘을 활용한다. 결정적으로, 각 $p$-큐비트 도트를 하나의 블록으로 취급하여 밀집된 도트 해밀토니안의 정확한 트로터 전파자(Trotter propagator)를 단일 $2^p \times 2^p$ 연산자로 적용한다. 이 "블록 MPS" 방식은 체인 구조를 활용하여 (고정된 $\chi$에 대해) 결합 차원을 $q$에 대해 선형적으로, $p$에 대해 다항식적으로 유지하며, 이는 기존의 게이트 단위 MPS 업데이트와 대조된다. 또한 소규모 시스템에 대해서는 정확 대각화(Exact Diagonalization, ED)와 비교하며, 변분 양자 허수 시간 진화(VarQITE) 안사츠(ansatz)와도 비교한다.
2. 양자 자원 추정: 양자 프로세서(특히 IBM heavy-hex 아키텍처에 매핑)에서의 단일 트로터 단계에 필요한 2-큐비트 게이트 깊이와 개수를 추정한다. 저자들은 $p$와 $q$에 따른 회로 깊이의 스케일링을 분석하여, SYK-도트 체인을 밀집 SYK 및 이전에 제안된 희소 SYK 모델과 비교한다.

주요 기여 및 결과

• 얽힘 예산(Entanglement Budget): 논문은 체인의 임의의 절단(cut)에 대한 얽힘 엔트로피의 엄격한 경계를 도출한다. 임의의 순수 상태에 대해, 얽힘은 $\lfloor p/2 \rfloor \ln 2$ (최대 도트 내 얽힘)와 최대 도트 간 얽힘의 합으로 제한된다. 기저 상태에서 도트 간 얽힘은 정적이며 주로 결합 $K$에 의해 결정된다. 이 "얽힘 예산"은 문제의 난이도가 도트 크기 $p$에 의해 설정되며, 도트의 수 $q$는 얽힘을 지수적으로 증가시키지 않으면서 시스템 크기를 선형적으로 확장함을 의미한다.
• 회로 깊이 독립성: 전체 시스템 크기 $N$에 따라 회로 깊이가 급격히 증가하는 밀집 SYK 또는 무작위 희소 모델과 달리, SYK-도트 체인의 단일 트로터 단계에서의 2-큐비트 깊이는 $p$에 의해서만 결정되며($\sim p^3$으로 스케일링), $q$와는 독립적이다(도트가 하드웨어의 서로 다른 영역에 매핑된 경우).
• 양자 우위를 위한 자원 추정: 블록 MPS의 고전적 비용($p$에 대해 지수적으로 증가하지만 $q$에 대해 선형적으로 증가)과 양자 게이트 깊이를 결합하여, 저자들은 양자 우위가 나타날 수 있는 조건을 추정한다. 이들은 $p \approx 15$에서 $21$사이의 도트 크기(가정:$q=5$ 및 완만한 도트 간 결합 차원)에서 기저 상태 준비를 위한 양자 우위가 나타날 수 있음을 발견했다.
• NISQ의 한계: 이 임계 크기($p \ge 15$)에서, 단일 트로터 단계는 충실도(fidelity) $1/e$를 유지하기 위해 $10^{-6}$ 차수의 2-큐비트 오류율을 요구한다. 현재의 NISQ 하드웨어는 $10^{-3}$에서 $10^{-4}$ 수준의 오류율로 작동하며, 이는 $p \le 5$에서 $8$까지만 충분하다. 이 더 작은 영역에서는 고전적 시뮬레이션이 여전히 저렴하게 수행된다. 따라서 저자들은 양자 하드웨어가 고전적 시뮬레이션을 능가할 수 있는 영역이 현재의 NISQ 장치의 범위를 벗어나 있다고 결론짓는다.
• 검증:
  • 정확도: 블록 MPS는 소규모 시스템($p=6, q=2$)에서 정확 대각화와 일치한다.
  • 수렴성: $N=40$ ($p=8, q=5$)에서, 기저 상태 에너지는 결합 차원 $\chi = 768$에서 수렴하며, 이는 정확한 차원인 $2^{20}$보다 훨씬 낮은 수준으로 블록 접근법의 효율성을 확인시켜 준다.
  • 얽힘 구조: 시뮬레이션은 도트 간 얽힘이 작으며 $K$에 의해 결정됨을 확인한다. 강한 결합($|K|=2$)에서도 도트 간 엔트로피는 도트 내부 규모에 비해 작다.
  • 변분 회로: $(p,q)=(5,2)$에서 고전적으로 훈련된 변분 회로는 정확한 값의 3.6% 이내의 에너지를 달성한다. 노이즈 시뮬레이션은 충실도가 $e^{-\bar{n}}$ ($\bar{n}$은 평균 게이트 오류 수)에 따라 감소함을 보여주며, 이는 자원 추정에 사용된 오류 가정을 검증한다.
  • 연산자 확산(Operator Spreading): 링크를 가로지르는 OTOC(Out-of-time-order correlators)는 $|K|$가 커짐에 따라 급격히 줄어드는 시간 척도에서 감쇠하며, 이는 기저 상태는 통제 가능하지만 실시간 역학(real-time dynamics)은 빠르게 체인을 가로지를 수 있음을 나타낸다.

의의 및 주장
본 논문은 현재의 하드웨어로 양자 우위를 입증한다고 주장하지 않는다. 대신, 그 의의는 양자 우위가 나타날 수 있는 이론적 영역을 위치시키는 데 있다. 저자들은 SYK-도트 체인이 다음과 같은 특성을 제공한다고 주장한다:

1. 도트 크기 $p$에 대해서는 지수적으로 증가하지만 도트 수 $q$에 대해서는 선형적으로 증가하는 고전적 비용을 갖는 모델을 제공한다.
2. 기저 상태의 "난이도"가 $p$와 $K$를 통해 조절 가능하다는 것을 보여준다.
3. 현재의 NISQ 장치가 도달할 수 없는 $10^{-6}$ 수준의 매우 낮은 오류율이 필요함을 보여줌으로써, 고전적 실행 가능성과 양자 하드웨어 능력 사이의 간극을 증명한다.

이 연구는 미래의 하드웨어 개발을 위한 벤치마크 역할을 하며, 밀집된 혼돈 모델에 대한 양자 우위는 아마도 현재의 NISQ 장치보다 2개 차수 더 나은 오류율(특히 도트 크기가 15 큐비트를 초과하는 경우)을 요구할 것임을 시사한다. 저자들은 본 결과가 기저 상태에 국한된 것이라고 언급하며, 긴 시간의 실시간 진화 과정에서는 도트 간 얽힘이 성장하여 자원 추정치를 변화시킬 수 있다고 덧붙였다.