A Mathematical Theory of Correlated Error Accumulation in Quantum Circuits

기술적 요약: 양자 회로의 상관된 오차 누적에 관한 수학적 이론

문제 정의
무작위 벤치마킹(Randomised Benchmarking, RB)과 같은 기존의 양자 벤치마킹 기법은 모든 회로 인스턴스와 오차원을 통합하여 단일 평균 오차율을 추출한다. 이는 마르코프 가설(Markovian assumptions) 하에서 평균 게이트 충실도(average gate fidelity)를 추정하는 데는 효과적이지만, 오차 구조를 스칼라 값으로 압축해 버리기 때문에 상관된 오차(correlated errors)나 장거리 의존성(long-range dependencies)을 탐지할 수 없다. 결과적으로, 현재의 방법들은 오차가 독립적인 장치와 강력하게 상관된 장거리 오차를 가진 장치를 구분할 수 없으며, 후자의 경우 오차 불충실도(infidelity)가 회로 복잡도(큐비트 수 $n$)에 따라 어떻게 확장되는지를 밝혀내지 못한다. 이러한 한계는 동일한 RB 비율을 보이면서도 규모 확장 시 서로 다른 스케일링 동작을 보일 수 있는 장치들을 구별하는 능력을 저해한다.

방법론
본 논문은 임의적인 가정에 의존하지 않고 상관된 오차의 특성을 규명하기 위해 엄격한 다층적 수학적 프레임워크를 개발한다. 이 접근 방식은 세 가지 명시적인 공리(axioms)를 기반으로 한다:

1. A1 (정상성, Stationarity): 게이트 오차는 평균이 0인 약정상 과정(weakly stationary process)을 형성한다.
2. A2 (정규 변이, Regular Variation): 게이트 오차의 상관 함수 $C(k)$는 거듭제곱 법칙 $C(k) \sim L \cdot k^{-\beta}$에 따라 감쇠한다(느리게 변하는 인자가 포함될 수 있음). 이는 합산 가능한($\beta > 1$) 영역과 장거리 의존성($0 < \beta < 1$) 영역을 모두 포괄한다.
3. A3 (섭동 영역, Perturbative Regime): 누적된 위상 분산은 섭동 경계($E[\phi^2] < \pi^2/4$) 내에 머물러 있어, 불충실도가 0.5 미만임을 보장한다.

핵심 분석 도구는 **상관 오차 누적 정리(Correlated Error Accumulation Theorem, CEAT)**이다. 저자는 회로 크기 $n$에 대한 총 누적 공분산을 나타내는 누적 연산자 $T_n(C) = \sum_i \sum_j C(|i-j|)$를 정의한다. **오일러-마클로린 전개(Euler–Maclaurin expansion)**와 정규 변이(regular variation) 이론을 사용하여, 본 논문은 이 연산자의 점근적 스케일링을 유도한다.

데이터로부터 스케일링 지수를 추정하기 위해, 본 프레임워크는 모델 $f(n; \theta) = C_0 n^\alpha + R$에 대한 비선형 최소제곱법(Non-Linear Least Squares, NLS) 추정을 사용한다. 이 추정량의 통계적 특성(일치성, 점근적 정규성, 유한 표본 편향)은 Jennrich(1969)와 White(1980)가 확립한 조건을 통해 도출된다. 또한, 본 프레임워크는 장치의 상태(health regimes)를 분류하기 위한 최적의 결정 임계값을 결정하기 위해 베이즈 위험(Bayes risk) 분석을 활용한다.

주요 기여

1. 상관 오차 누적 정리 (CEAT): 상관 시퀀스 $C(k) \sim k^{-\beta}$에 대해, 누적 오차 연산자가 $T_n(C) \sim n^\alpha$로 스케일링되며, 이때 지수는 $\alpha = \max(1, 2-\beta)$임을 증명한다.
   • 만약 $0 < \beta < 1$ (장거리 의존성)이면, $\alpha = 2-\beta > 1$이 되어 초선형(super-linear) 오차 성장을 나타낸다.
   • 만약 $\beta > 1$ (합산 가능한 상관관계)이면, $\alpha = 1$이 되어 선형 성장을 나타낸다.
   • 경계 사례인 $\beta=1$은 로그 보정($n \ln n$)을 생성한다.
2. 안정성 및 보편성 정리: 지수 $\alpha$가 가법적 섭동(주요 상관관계보다 빠르게 감쇠하는 경우)과 곱셈적 느리게 변하는 인자, 그리고 유한한 절단(finite truncations)에 대해 강건함을 입증한다. 이 결과는 특정 물리적 메커니즘이 아닌, 상관 함수의 정규 변이 클래스에만 의존하는 보편적인 것으로 나타났다.
3. 추정량의 통계적 특성 규명: NLS 추정량 $\hat{\alpha}$가 일치성을 가지며 점근적으로 정규 분포를 따름을 증명한다. 결정적으로, 이 추정량이 $O(n^{\beta-1})$로 감소하는 **유한 표본 상향 편향(finite-sample upward bias)**을 정량화한다. 이를 통해 특정 추정 정확도를 달성하기 위해 필요한 최소 큐비트 수($N(\beta, \epsilon)$)를 계산할 수 있다.
4. GHZ 회로를 통한 물리적 연결: 저자는 GHZ 상태 준비 회로가 누적된 상관 오차를 직접 조사하는 프로브 역할을 함을 입증한다. 명시된 공리 하에서, $n$-큐비트 GHZ 상태의 불충실도는 $T_n(C)/4$에 비례한다. 교대 급수(alternating series)를 통해 섭동 경계 $|C_{GHZ} - \phi^2/4| \leq \phi^4/48$이 증명되었으며, 이는 GHZ 회로의 진단적 유효성을 검증한다.
5. 의사결정론적 프레임워크: $n=2\text{--}10$ 범위에서 확률적(건강한) 상태와 임계(상관된) 상태를 구분하기 위한 베이즈 최적 임계값($\alpha^* \approx 1.57$)을 도출하며, 유한 데이터셋의 경우 $\alpha=1$이라는 점근적 임계값이 최적이 아님을 언급한다.

결과

• 분석적 검증: 일반적인 거듭제곱 법칙 상관관계에 대해 CEAT가 분석적으로 증명되었다. 두 항 전개식 $T_n(C) \approx A n^{2-\beta} + B n$은 $n \geq 1000$에서 상대 오차가 $0.01\%$ 미만임을 보였다.
• 수치적 검증: 2000회의 시행을 포함한 몬테카를로 시뮬레이션은 NLS 추정량의 점근적 정규성과 예측된 상향 편향을 확인하였다. 편향은 $\beta=1$ 근처에서 가장 심각한 것으로 나타났다.
• 하드웨어 검증 맥락: 본 논문에는 실험 데이터가 포함되어 있지 않으나, 세션 7개에 걸쳐 세 대의 IBM Heron r2 프로세서($ibm\_fez$, $ibm\_kingston$, $ibm\_marrakesh$)의 데이터를 적용한 동반 실험 논문 [6]을 참조한다.
  • 해당 연구에서 거듭제곱 법칙 피팅은 모든 세션에서 $R^2 > 0.97$을 기록했다.
  • 장치들은 뚜렷한 영역을 보였다: $ibm\_fez$와 $ibm\_kingston$은 $\alpha < 1.569$ (확률적/건강함)를 보였으나, $ibm\_marrakesh$는 한 세션에서 $\alpha = 3.78$ (임계/결맞음)을 보였으며, 이는 이후의 유지보수 이벤트 및 T2 회복과 상관관계가 있었다.
  • 노이즈 분광학 추정치($\mu$)와 $\alpha = 1+\mu$ 관계를 통해 $\hat{\alpha}$를 비교하는 전향적 예측 테스트(Prediction P1)를 $ibm\_fez$에 수행하였다. 8/8 세션에서 모델이 기각되지 않았으며, 이는 해당 장치에 대한 단일 메커니즘 가설을 지지한다.
  • $ibm\_marrakesh$에 대한 검증 테스트(Prediction P2)는 테스트 통계량 $Z=7.29$를 산출하여 단일 메커니즘 모델을 결정적으로 기각하였으며, 이는 Born 근사를 넘어서는 결맞은 ZZ 결합 또는 다른 메커니즘의 존재를 나타낸다.

의의 및 주장
본 논문은 회로 불충실도의 스케일링 지수와 게이트 오차의 상관 감쇠 지수를 연결하는 최초의 엄격한 수학적 이론을 제공한다고 주장한다. 본 논문은 현상학적인 거듭제곱 법칙 피팅과 달리, 제1원리로부터 추정량의 통계적 특성을 유도하고 이론이 성립하는 조건(A1–A3)을 확립함으로써 차별화된다.

저자는 이 프레임워크가 **조건부적(conditional)**임을 강조한다: 수학적 정리는 무조건적으로 증명되었으나, 이를 특정 하드웨어에 적용하는 것은 공리의 경험적 타당성에 달려 있다. 본 논문은 K–R 파라미터화가 CEAT의 적용이지 유도가 아님을 명시한다. 또한, A2 공리는 본 텍스트 내에서 근본적인 욕조 역학(bath dynamics)으로부터 유도된 것이 아니라 가정된 것이다. A2의 유도는 열린 양자계 이론(Nakajima–Zwanzig 방정식 및 Tauberian 정리)을 사용하는 향후 "Level 2" 프로그램의 과제로 제시되었다.

이 연구의 의의는 단순히 평균 오차율을 측정하는 것을 넘어 노이즈 메커니즘(상관된 오차 vs 독립적인 오차)을 드러내는 진단 도구를 제공한다는 데 있다. 스케일링 지수 $\alpha$를 식별함으로써, 연구자들은 표준 벤치마킹 하에서는 유사해 보일 수 있지만 근본적으로 다른 오차 누적 동작을 가진 장치들을 구분할 수 있으며, 이는 양자 하드웨어의 확장성에 매우 중요하다. 논문은 마지막으로 프레임워크의 가정을 체계적으로 테스트하고 실패 모드를 식별하기 위한 "Ten-Row Falsification Protocol"을 제시하며, 이 이론을 블랙박스 휴리스틱이 아닌 검증 가능한 과학적 모델로 포지셔닝하며 마무리한다.