Certified Fidelity Susceptibility from Classical Shadows

기술 요약: 클래식 섀도우를 이용한 인증된 충실도 민감도(Certified Fidelity Susceptibility)

문제 정의
충실도 민감도($\chi_F$)는 양자 상전이를 탐지하고 양자 메트롤로지를 수행하는 데 있어 중요한 프로브 역할을 하며, 질서 매개변수(order parameter)에 대한 사전 지식 없이도 탐지 능력을 제공한다. 그러나 이를 양자 장치에서 직접 평가하는 것은 두 가지 주요 과제로 인해 저해된다:

1. 스펙트럼 요구 사항: 표준 섭동 전개(perturbative expansion)는 전체 들뜬 상태 스펙트럼과 전이 행렬 요소(transition matrix elements)를 필요로 하는데, 이는 일반적인 바닥 상태 준비 시나리오에서는 접근 불가능하다.
2. 노이즈 증폭: 인접한 바닥 상태 사이의 충실도를 통해 $\chi_F$를 추정하는 유한 차분(finite-difference) 방식은 통계적 노이즈를 증폭시킨다.
   클래식 섀도우 토모그래피(classical shadow tomography)가 적은 측정 횟수로 많은 관측량을 추정할 수 있는 경로를 제공하지만, 참조 바닥 상태의 단일 복사(single-copy) 측정으로부터 $\chi_F$를 직접 재구성하는 것은 여전히 미해결 과제로 남아 있다. 변분적 접근법(variational approaches)이나 양자 특이값 변환(QSVT)과 같은 기존 방법들은 종종 깊은 회로를 요구하거나 실제 환경에서의 보증(guarantee)이 부족하다.

방법론
저자는 **크릴로프 부공간 방법(Krylov-subspace methods)**과 충실도 민감도의 **분해자 기반 재정식화(resolvent-based reformulation)**를 통합하는 프레임워크를 제안한다. 핵심 전략은 비국소적 연산자 문제를 클래식 모멘트 재구성 문제로 변환하는 것이다.

1. 분해자 재정식화 (Resolvent Reformulation):
   전체 스펙트럼에 대해 $\chi_F$를 전개하는 대신, 저자는 연결된 섭동 연산자 $B = H_I - \langle \Psi_0 | H_I | \Psi_0 \rangle I$와 상태 $|\phi\rangle = B|\Psi_0\rangle$를 정의한다. 충실도 민감도는 다음과 같이 재작성된다:
   $$\chi_F = \langle \phi | (K^+)^2 | \phi \rangle$$
   여기서 $K = H - E_0 I$이고 $K^+$는 무어-펜로즈 유사 역행렬(Moore–Penrose pseudoinverse)이다. 이 정식화는 들뜬 상태에 대한 명시적 의존성을 피하면서도 비국소적 연산자 $(K^+)^2$를 유지한다.

2. 크리롤 부공간 근사 (Krylov Subspace Approximation):
   비국소적 연산자를 처리하기 위해, 이 방법은 벡터 $|y\rangle = K^+|\phi\rangle$를 $n$차 크릴로프 부공간 $\mathcal{V}_n = \text{span}\{K|\phi\rangle, \dots, K^n|\phi\rangle\}$ 내에서 근사한다. 정확한 민감도는 $\|y\|^2$이다. 근사치 $\chi_F^{(n)}$는 $\mathcal{V}_n$ 위로의 $|y\rangle$의 직교 투영(orthogonal projection)의 제곱 노름으로 정의된다.

이는 다음과 같은 **단조 증가하는 하한 계층(monotone hierarchy of lower bounds)**을 도출한다:
$$0 \le \chi_F^{(1)} \le \chi_F^{(2)} \le \dots \le \chi_F$$
$n$차 근사치는 응답 모멘트 $M_k = \langle \phi | K^k | \phi \rangle$에 의해 완전히 결정된다. 구체적으로, $\chi_F^{(n)} = \mathbf{m}_n^\dagger \Gamma_n^{-1} \mathbf{m}_n$이며, 여기서 $\mathbf{m}_n$은 모멘트 벡터이고 $\Gamma_n$은 $M_0$부터 $M_{2n}$까지의 모멘트로 구성된 그람(Gram, 또는 Hankel) 행렬이다.

3. 클래식 섀도우 추정 (Classical Shadow Estimation):
   각 모멘트 $M_k$는 바닥 상태 밀도 행렬 $\rho_0$의 선형 범함수(linear functional)이므로, 단일 세트의 무작위 단일 복사 측정(클래식 섀도우)으로부터 동시에 추정될 수 있다. 이 프로토콜은 스냅샷을 생성하기 위해 독립적인 무작위 단일 큐비트 파울리 측정을 사용한다.

저자는 추정된 모멘트들이 균등 오차 경계 $\eta$를 만족할 때, 실제 충실도 민감도가 (확률 $1-\delta$로) 다음을 만족함을 보장하는 인증된 하한(certified lower bound) $B^{(n)}$을 도출한다:
$$\chi_F \ge \Omega^{-2} \left[ \tilde{\chi}_F^{(n)} - L_n \eta - R_n \eta^2 \right]$$
여기서 $\tilde{\chi}_F^{(n)}$는 노이즈가 있는 모멘트로 재구성된 값이고, $L_n, R_n$은 추정된 그람 행렬의 조건수(condition number)와 오차 $\eta$에 의존하는 패널티 항이다.

주요 결과 및 이론적 보증

• 단조성 및 인증: 이 방법은 정확한 값으로 단조 증가하며 인증되는 하한 시퀀스를 제공한다. 인증은 섀도우 데이터의 통계적 특성에 기초하여 엄격하게 이루어진다.
• 기하급수적 수렴 (Geometric Convergence): 오차 $\chi_F - \chi_F^{(n)}$는 크릴로프 차수 $n$에 따라 기하급수적으로 수렴한다. 수렴 속도는 전체 해밀토니안 스펙트럼이 아니라 응답 스펙트럼의 조건수($\kappa_\phi$)에 의해 결정되며, 이는 스펙트럼이 희소한 응답을 가진 시스템에서 빠른 수렴을 가능하게 할 잠재력이 있다.
• 샘플 복잡도: $2n+1$개의 모멘트를 동시에 추정하는 데 필요한 샷(shot)의 수는 모멘트 수에 대해 로그 스케일($O(\log n)$)로 증가하여, 적절한 크기의 크릴로프 차수에 대해 효율적이다.
• 강건성 (Robustness): 이 프레임워크는 불완전한 바닥 상태 준비를 다룰 수 있도록 확장되었으며, 수정된 오차 예산 내에서 인증이 유효하게 유지된다.

응용 분야
저자는 두 가지 즉각적인 응용 사례를 제시한다:

1. 양자 메트롤로지: 순수 상태에 대해 양자 피셔 정보(QFI)는 $F_\lambda = 4\chi_F$이므로, $\chi_F$에 대한 인증된 하한은 QFI에 대한 통계적으로 인증된 하한을 직접적으로 제공한다.
2. 선형 정적 민감도 (Linear Static Susceptibilities): 이 구성은 정적 외부 장에 대한 응답을 정량화하는 선형 정적 민감도($\chi_O$)로 일반화된다. 이는 이동된(shifted) 크릴로프 부공간과 모멘트를 필요로 하며, 유사한 단조 증가 하한 계층을 생성한다.

의의 및 주장
저자는 이 연구가 바닥 상태의 무작위 단일 복사 측정만을 사용하여 충실도 민감도를 인증하는 첫 번째 방법임을 주장하며, 이를 통해 들뜬 상태의 스펙트럼이나 깊은 회로의 필요성을 제거했다.

• 운용 가치: 단조 증가하며 기하급수적으로 수렴하는 하한 계층을 제공함으로써, 이 방법은 근접 양자 장치(NISQ)에서 양자 상전이를 탐지하고 메트롤로지적 유용성을 평가하는 실용적인 도구를 제공한다.
• 효율성: 이 접근법은 양자 추정 문제를 클래식 모멘트 재구성 문제로 축소하여, 클래식 섀도우 토모그래피의 효율성을 활용한다.
• 향후 방향: 저자는 현재 작업이 특정 함수($1/x^2$ 및 $1/x$)에 집중되어 있지만, 프레임워크가 일반적인 함수 $f(x)$ 및 혼합 상태(mixed states)로 확장될 수 있음을 언급했다. 또한 샘플 복잡도와 측정 앙상블(Pauli vs. Clifford) 선택 사이의 트레이드오프를 향후 분석 과제로 꼽았다.

결론적으로, 이 논문은 크릴로프 기법을 분해자 정식화 및 클래식 섀도우와 통합하는 것이 양자 하드웨어 상의 양자 시스템을 특성화하기 위한 강력하고 확장 가능하며 인증 가능한 경로를 제공한다고 밝히고 있다.