Operator delocalization in disordered spin chains via exact MPO marginals

1. 연구 배경 및 문제 제기 (Problem)

• 양자 열화 및 스크램블링: 양자 열화 (thermalization) 는 엔트로피 증가와 함께 정보가 비국소적 자유도로 퍼져가는 '스크램블링 (scrambling)' 과정을 수반합니다. 이를 정량화하기 위해 기존에는 주로 **OTOC (Out-of-Time-Order Correlator)**가 사용되어 왔으나, 이는 특정 연산자 기저에 의존하지 않는 공간 - 시간 정보를 제공하지만, 연산자의 구체적인 구조적 복잡성 (예: 어떤 기저에서 얼마나 퍼져 있는지) 을 직접적으로 측정하기에는 한계가 있습니다.
• 무질서한 계의 동역학:
  • 앤더슨 국소화 (Anderson Localization): 상호작용이 없는 무질서한 계에서는 모든 단일 입자 상태가 국소화되어 연산자의 확산이 억제됩니다.
  • 다체 국소화 (MBL, Many-Body Localization): 상호작용이 있는 무질서한 계에서는 열화가 억제되지만, 로그 (logarithmic) 형태의 빛의 원뿔 (light cone) 을 따라 매우 느린 속도로 정보가 확산됩니다.
• 기존 방법론의 한계: 연산자의 복잡성을 측정하는 지표인 **비안정화 엔트로피 (Nonstabilizerness)**나 **스태빌라이저 레니 엔트로피 (SRE)**를 계산하려면 확률적 샘플링 (stochastic sampling) 이 필요하며, 시스템 크기가 커질수록 샘플 수가 지수적으로 증가하여 계산 비용이 매우 큽니다. 또한, 연산자의 공간적 확산 범위를 정밀하게 측정할 수 있는 효율적인 도구가 부족했습니다.

2. 방법론 (Methodology)

이 논문은 행렬 곱 연산자 (MPO, Matrix Product Operator) 프레임워크를 기반으로 두 가지 새로운 연산자 측정 지표를 도입하고 이를 정밀하게 계산하는 방법을 제시합니다.

• 새로운 지표 도입:
  1. 연산자 질량 (Operator Mass, $m(t)$): 시간 진화한 연산자를 파울리 기저 (Pauli basis) 로 전개했을 때, 항등원 (Identity) 이 아닌 파울리 행렬이 포함된 사이트의 평균 개수. 연산자가 얼마나 많은 사이트에 걸쳐 작용하는지 (내부 밀도) 를 나타냅니다.
  2. 연산자 길이 (Operator Length, $h(t)$): 파울리 스트링에서 항등원이 아닌 가장 오른쪽 사이트의 위치. 연산자의 공간적 지지 범위 (spatial support) 의 확장을 나타냅니다.
• 정밀 MPO 계산 기법:
  • 연산자를 힐베르트 - 슈미트 내적을 통해 파울리 스트링 기저로 전개합니다.
  • 이를 정규화된 파울리 기저를 사용하여 MPS (Matrix Product State) 형태로 매핑합니다 (Choi-Jamiołkowski 동형사상 활용).
  • 핵심 기여: 기존 확률적 샘플링 없이, MPO 의 텐서 네트워크 구조를 직접 활용하여 연산자 질량과 길이의 전체 확률 분포를 정밀하게 (exact) 그리고 효율적으로 계산합니다.
  • 이 방법은 연산자 엔트로피 (Operator Entanglement Entropy) 와 결합하여, MBL 위상에서의 연산자 성장 동역학을 추적합니다.

3. 주요 결과 (Key Results)

연구자들은 무질서한 XXZ 스핀-1/2 사슬 (disordered XXZ spin-1/2 chain) 을 대상으로 수치 시뮬레이션 (ED 및 MPS) 을 수행하고 다음과 같은 결과를 도출했습니다.

• 상호작용 유무에 따른 극명한 차이:
  • 비상호작용계 ($\Delta = 0$, 앤더슨 국소화): 연산자 질량, 길이, 엔트로피가 모두 매우 빠르게 포화 (saturation) 됩니다. 이는 정보가 확산되지 않고 국소화됨을 의미합니다.
  • 상호작용계 ($\Delta \neq 0$, MBL): 약한 상호작용만으로도 모든 양 (질량, 길이, 엔트로피) 이 시간에 대해 로그 ($\ln t$) 형태로 성장합니다. 이는 MBL 에서 관찰되는 로그 빛의 원뿔 (logarithmic light cone) 과 일치합니다.
• 로그 성장의 물리적 기작:
  • 이 로그 성장은 $\ell$-bit 모델 (국소화된 적분 운동량) 로 설명됩니다. 국소화된 $\ell$-bits 간의 상호작용이 거리에 따라 지수적으로 약해지기 때문에 ($J_{jl} \sim e^{-\kappa r}$), 연산자가 거리 $r$만큼 퍼지기 위해서는 $t \sim e^{r/\kappa}$의 시간이 소요되어 $r(t) \sim \ln t$ 관계를 유도합니다.
  • 상호작용이 약하더라도 ($\Delta > 0$), 앤더슨 국소화는 파괴되며 비국소화가 시작됨을 확인했습니다.
• 시스템 크기 및 무질서 강도 의존성:
  • 시간 평균된 연산자 길이와 질량을 시스템 크기 $L$로 스케일링 ($h/\ln L$) 하면, 강한 무질서 영역에서 데이터가 하나의 보편적 곡선으로 수렴 (collapse) 합니다. 이는 MBL 위상의 로그 확산 특성을 강력하게 지지합니다.
  • 중간 정도의 무질서 영역에서는 비단조적인 거동을 보이며, 이는 시스템이 적분 가능 (integrable) 한 지점과 국소화 지점 사이를 전이함을 시사합니다.

4. 의의 및 기여 (Significance)

• 계산적 효율성: 확률적 샘플링 없이 MPO 를 통해 연산자 질량과 길이의 전체 분포를 정밀하게 계산할 수 있는 방법을 최초로 제시했습니다. 이는 시스템 크기에 대해 다항식 시간 복잡도로 계산이 가능하여 대규모 계의 장시간 동역학을 연구할 수 있게 합니다.
• MBL 동역학의 새로운 통찰: 연산자 '길이'와 '질량'이라는 물리적으로 직관적인 지표를 통해 MBL 의 느린 스크램블링 과정을 명확하게 규명했습니다. 특히, 연산자 엔트로피와 공간적 확산이 동일한 로그 스케일로 성장함을 보여줌으로써 MBL 의 내부 복잡성과 공간적 확장의 밀접한 연관성을 입증했습니다.
• 실험적 검증 가능성: 제안된 프로토콜은 양자 시뮬레이션 플랫폼 (초전도 큐비트, 트랩드 이온 등) 에서 실행 가능합니다.
  • Choi-Jamiołkowski 동형사상을 이용한 시스템 - 안실라 (ancilla) 상태 준비.
  • 클래식 쉐도우 (Classical Shadows) 기법과 벨 측정 (Bell measurements) 을 결합하여 연산자 한계값 (marginals) 을 효율적으로 추정.
  • 이는 실험실 환경에서 MBL 현상과 연산자 확산을 직접 관측할 수 있는 길을 열었습니다.

5. 결론

이 논문은 무질서한 상호작용 스핀 사슬에서 연산자의 비국소화를 정량화하기 위해 연산자 길이 (Operator Length) 개념을 도입하고, 이를 정밀 MPO 기법을 통해 계산 가능한 지표로 확립했습니다. 연구 결과는 상호작용이 있는 MBL 계에서 연산자가 로그적으로 느리게 확산됨을 확인했으며, 이는 기존 OTOC 연구 결과와 일치하면서도 연산자의 내부 구조와 공간적 범위를 더 구체적으로 보여줍니다. 또한, 제안된 방법은 이론적 계산의 효율성을 높였을 뿐만 아니라, 향후 양자 시뮬레이션을 통한 실험적 검증을 위한 구체적인 로드맵을 제공한다는 점에서 중요한 의의를 가집니다.