Zero noise extrapolation of transport coefficients

기술 요약: 수송 계수의 제로 노이즈 외삽법 (Zero Noise Extrapolation)

문제 정의
제로 노이즈 외삽(Zero-noise extrapolation, ZNE)은 잡음이 있는 중간 규모 양자(NISQ) 컴퓨팅에서 표준적으로 사용되는 프로토콜로, 다양한 노이즈 강도($\gamma$)에서 관측량을 측정하고 이를 제로 노이즈 극한($\gamma \to 0$)으로 외삽하여 이상적인 시스템의 성질을 추정한다. ZNE는 노이즈 강도에 대해 해석적(analytically)으로 의존하는 관측량에는 효과적이지만, 보존량의 수송 계수(transport coefficients)를 추정할 때는 실패하는 것으로 알려져 있다. 이러한 실패는 노이즈가 일반적으로 보존 법칙을 위반하여, 임의의 양수 $\gamma$에서 수송을 소멸시키기 때문에 발생한다. 또한, 노이즈가 평균적으로 보존 법칙을 유지하더라도, 정지 상태(stationary state)를 변화시키거나 보편적인 유체역학적 보정(hydrodynamic corrections)으로 인해 비해석적(non-analytic) 의존성을 유발할 수 있다. 결과적으로, 표준 ZNE 방식은 확산 및 아확산(subdiffusion) 계수를 잘못 추정하는 경우가 빈번하다.

방법론
저자들은 GKLS(Gorini-Kossakowski-Lindblad-Sudarshan) 마스터 방정식을 통해 모델링된 노이즈 강도 $\gamma$를 갖는 1차원 비적분 스핀 체인(non-integrable spin chains)에서 ZNE를 조사한다. 이 연구는 해석적 유체역학 논증과 ITensor 텐서 네트워크 방법을 이용한 수치 시뮬레이션을 결합하여 수행된다.

1. 이율 가설 (Two-Rate Ansatz): 저자들은 $\gamma > 0$일 때, 보존량(스핀 밀도)의 푸리에 모드 감쇠율 $R(k, \gamma)$가 다음과 같은 형태를 따른다고 제안한다:
   $$R(k, \gamma) = D_{\text{eff}}(\gamma)k^z + R(0, \gamma)$$
   여기서 $z$는 동역학적 지수(dynamical exponent)이며, $R(0, \gamma)$는 총 전하의 감쇠율이다. 이를 통해 실제 수송이 억제되는 상황에서도 유효한 노이즈 확산 상수 $D_{\text{eff}}(\gamma)$를 정의할 수 있다.

2. 정지 상태 진화 처리: 노이즈가 시스템의 정지 상태(자화 밀도 변화)를 변화시키는 경우, 저자들은 $\gamma t$의 곱을 고정한 채 $\gamma \to 0$와 $t \to \infty$의 극한을 동시에 취함으로써 올바른 제로 노이즈 극한을 회복할 수 있다고 주장한다. 이는 외삽이 노이즈 유도 정지 상태가 아닌 특정 자화 밀도 $M^*$를 목표로 하도록 보장한다.

3. 유체역학적 보정: 논문은 유동 유체역학(fluctuating hydrodynamics)을 사용하여 $D_{\text{eff}}(\gamma)$의 노이즈 의존성을 도출한다. 저자들은 $D_{\text{eff}}(\gamma)$가 외삽된 값 $D_{\text{exp}}$, 해석적 보정 항 $f_a(\gamma)$, 그리고 전류-전류 상관 함수(current-current correlators)의 보편적인 장기 꼬리(long-time tails)에서 기인하는 비해석적 보정 항 $f_h(\gamma)$로 구성된다고 상정한다:
   $$D_{\text{eff}}(\gamma) = D_{\text{exp}} + f_a(\gamma) + f_h(\gamma)$$
   비해석적 항은 $\gamma^{\alpha-1}$로 스케일링되며, 여기서 $\alpha$는 무잡음 극한에서의 전류 상관 함수의 멱법칙 붕괴($\sim \tau^{-\alpha}$)에 의해 결정된다.

주요 결과
저자들은 두 가지 특정 모델에 대한 수치 시뮬레이션과 분석적 예측을 벤치마킹한다:

1. 소산적 하이젠베르크 체인 (확산 수송):

   • 다중 보존 법칙 (SU(2) 대칭): 스태거드 필드(staggered field)가 없는 경우, 시스템은 여러 개의 보존량을 가진다. 이 경우 전류-전 전류 상관 함수는 $\tau^{-3/2}$ ($\alpha = 3/2$)로 붕괴하며, 이는 $\sqrt{\gamma}$로 스케일링되는 비해석적 보정을 초래한다. 수치 데이터와 섭동 계산은 $D_{\text{eff}}(\gamma) \approx D_{\text{exp}} - c\sqrt{\gamma}$임을 확인해 준다.
   • 단일 보존 법칙 (스태거드 필드): 스태거드 필드가 SU(2) 대칭을 깨뜨려 단 하나의 보존량만 남기는 경우, 주요 비해석적 보정이 사라지거나(또는 부차적이 되어) 의존성이 해석적(선형 $\gamma$)이 된다.

2. 다이폴 보존 프랙톤 체인 (아확산 수송):

   • 저자들은 총 다이폴 모멘트를 보존하며 무잡음 극한에서 아확산 수송($z_{\text{true}} = 4$)을 보이는 체인을 분석한다.
   • 저자들은 유효 아확산 계수가 $\gamma^{1/2}$에 따라 제로 노이즈 값에 접근함을 입증한다. 이는 보편적인 유체역학적 보정이 아확산 영역에서도 비해석적 거동을 생성함을 확인시켜 준다.

3. 외삽 프로토콜: 연구는 $\gamma t$를 고정함으로써(즉, 자화 밀도를 고정함으로써) 정지 상태가 노이즈에 의존하는 경우에도 진정한 무잡음 확산 상수 $D_{\text{true}}(M^*)$를 성공적으로 회복할 수 있음을 검증한다.

의의 및 주장
본 논문은 표준 ZNE가 수송 계수를 잘못 추정하는 이유가 단순히 노이즈가 보존 법칙을 파괴하기 때문이 아니라, 유체역학에 대한 보편적인 비선형 보정이 노이즈 강도에 대한 비해석적 의존성을 유발하기 때문이라고 주장한다. 구체적으로:

• 표준 ZNE는 $\gamma$에 대한 해석성을 가정하지만, 다중 보존량(또는 아확산 시스템의 특정 대칭성)이 존재하는 시스템에서는 주요 보정 항이 비해석적(예: $\gamma^{1/2}$)이다.
• 유효 수송 계수는 다중 보존 법칙이 존재하는 경우 $\alpha < 2$인 $\gamma^{\alpha-1}$ 멱법칙을 통해 제로 노이즈 값에 접근한다.
• 저자들은 특정 자화 밀도에 대응하는 수송 계수를 정확히 식별하기 위해 $\gamma t$를 고정하는 수정된 외삽 프로토콜을 제공한다.

이 연구는 이러한 비해석적 유체역학적 꼬리를 고려하지 않으면 ZNE가 양자 시뮬레이터의 수송 계수를 체계적으로 오측정하게 될 것임을 보여주는 이론적 및 수치적 벤치마크 역할을 한다.