Diameter truncated operator evolution

기술 요약: 직경 절단 연산자 진화 (Diameter-Truncated Operator Evolution, DTOE)

문제 정의
평형에서 벗어난 다체 양자계(out-of-equilibrium many-body quantum systems)의 역학을 시뮬레이션하는 것은 복잡성의 급격한 증가로 인해 근본적으로 어렵습니다. 이러한 계가 진화함에 따라 얽힘(entanglement), 비안정성(non-stabilizerness), 비가우시안성(non-Gaussianity)이 증가하며, 이로 인해 초기에 국소적이었던 연산자들이 매우 비국소적이고 "스크램블링(scrambled)"된 상태가 됩니다. 이러한 성장은 초기 파동함수나 가능한 모든 관측량을 효율적으로 포착할 수 없게 만들어, 전통적인 텐서 네트워크 시뮬레이션을 매우 짧은 시간대 너머로 수행하는 것을 불가능하게 만듭니다. 후기 열화(late-time thermalization)는 유체역학적 기술을 허용하지만, 단기 및 중기 단계에서 평형으로 접근하는 과정을 정확하게 계산하는 것은 여전히 어려운 과제입니다.

방법론
저자들은 제한된 공간적 범위를 가진 연산자의 부분 공간으로 시뮬레이션을 제한하도록 설계된 연산자 역학 시뮬레이션 수치 방법인 **직경 절단 연산자 진화(DTOE)**를 소개합니다.

• 핵심 개념: 연산자의 가중치(Pauli 항의 개수)를 기준으로 연산자를 절단하는 기존의 소산 보조 연산자 진화(DAOE)나 Pauli 전파(Pauli Propagation) 방식과 달리, DTOH는 **직경(diameter)**을 기준으로 절단합니다. 직경은 격자 상에서 연산자가 비자명하게 지지(support)되는 실제 공간 영역의 크기로 정의됩니다.
• 구현: 이 방법은 1+1D 브릭워크 양자 회로(brickwork quantum circuits) 프레임워크 내에서 정식화됩니다. 연산자를 이중 힐베르트 공간으로 벡터화하는 "폴딩(folded)" 표기법을 활용합니다. 시간 진화는 직경 $d$를 가진 연산자 구간을 매핑하는 전달 행렬(transfer matrices)을 적용함으로써 계산됩니다.
• 절단 방식: 하드 절단(hard truncation)이 적용됩니다. 직경이 고정된 컷오프 $d$를 초과하는 모든 연산자는 시뮬레이션에서 엄격히 제외됩니다. 상관 함수는 이러한 전달 행렬들을 체이닝하여 로컬 연산자 사이의 경로를 생성함으로써 재구성됩니다.
• 이론적 동기: 이 접근법은 연산자의 대부분의 노름(norm) 흐름이 매우 비국소적인 궤적으로 흘러가더라도, 국소적인 이점 상관 함수(two-point correlation functions)는 종종 유계된 상수 직경 영역 내에 머무는 "이례적인(atypical)" 궤적들에 의해 지배된다는 관찰에 기초합니다. 일반적인 회로의 경우, 시공간 상에 (통상 임계 속도 $v_c$ 이상에서) 상관 관계가 상수 $O(1)$ 직경의 연산자들에 의해 지배되는 "쉬운(easy)" 상(phase)이 존재합니다. 이 영역에서 DTOE는 $O(t)$의 메모리 스케일링과 함께 다항 시간 복잡도($O(t^2)$)를 제공합니다.

주요 기여 및 결과
본 논문은 고성능 컴퓨팅 자원을 사용하지 않고도 가적분(integrable) 및 비가적분(non-integrable) 계 모두에 대한 광범위한 수치 시뮬레이션을 통해 DTOE의 효능을 입증합니다.

1. 섭동된 이중 유니타리 회로 (Kicked Ising Model):

   • 저자들은 $J=B=\pi/4+\eta$ 및 $h=1.2$인 키크드 이징 모델(kicked Ising model)을 시뮬레이션하였으며, 섭동 강도 $\eta \in [0, 0.5]$를 변화시켰습니다.
   • 무한 온도 이점 상관 함수($C_{XX}$)를 $t=100$ 회로 레이어까지 계산하였습니다.
   • 결과: $\eta = 0.24588$일 때, 낮은 직경 근사($d=1, 2$)는 정확한 데이터에서 벗어났으나, 더 높은 직경($d \ge 9$)에서는 $10^{-2}$에서 $10^{-3}$ 수준의 상대 오차를 달 achievement했습니다. 이 방법은 매우 비섭동적인 영역에서도 효과적이었으며, 일반적으로 약 25회의 플로케(Floquet) 사이클만을 수행할 수 있는 현재의 디지털 양자 컴퓨터 시뮬레이션보다 뛰어난 성능을 보였습니다.

2. Heisenberg XXZ 모델:

   • 상관 폭 $\sigma(t) \sim t^\alpha$의 스케일링을 위한 유체역학적 수송 지수($\alpha$)를 추출하기 위해 DTOE를 XXZ 회로에 적용하였습니다.
   • 탄성(Ballistic) 영역 ($|\Delta| < 1$): $d=10$인 DTOE는 $\alpha = 1.00 \pm 0.01$을 정확하게 추출하였습니다. 저자들은 방법이 확산(diffusion)으로 드리프트하기 전까지 정확도를 유지하는 시간 척도 $t^*$가 직경 $d$에 대해 지수적으로 스케일링됨($t^* \sim O(\exp(d))$)을 관찰했습니다.
   • 초확산(Superdiffusive) 영역 ($|\Delta| = 1$): 등방성 지점에서, 이 방법은 긴 준입자 스트링(quasiparticle strings)을 포착하는 데 따르는 알려진 어려움에도 불구하고, 이론적으로 예측된 지수 $\alpha = 2/3$ (구체적으로 $d=10$일 때 $\alpha \in [0.6663, 0.6684]$)를 포착했습니다.
   • 확산(Diffusive) 영역 ($|\Delta| > 1$): 방법은 낮은 직경에서도 확산 지수 $\alpha = 1/2$를 정확하게 재현하였는데, 이는 이 영역에서 수송이 짧은 준입자 스트링에 의해 지배되기 때문입니다.

의의 및 주장
본 논문은 DTOE가 가중치 기반 방법들에 비해 더 "가벼운(leaner)" 절단 프로토콜을 제공한다고 주장합니다. 이론적 분석에 따르면, 크기가 $L$인 격자에서 직경 $d$를 가진 연산자의 수는 $O(L)$로 스케일링되는 반면, 가중치 $d$를 가진 연산자의 수는 $O(L^d)$로 스케일링됩니다. 따라서 큰 $L$에 대해 DTOE와 가중치 기반 절단 간의 유지되는 연산자 비율은 0으로 수렴하며, 이는 DTOE가 본질적으로 더 효율적임을 시사합니다.

저자들은 DTOE가 국소 상관 관계가 유계된 직경 궤적들에 의해 지배되는 계에 대해 물리적으로 잘 동기화되어 있다고 주장합니다. 그들은 DTOE가 현재의 양자 하드웨어 능력을 넘어서는 시간대에서 수송 특성과 상관 함수를 성공적으로 추출한다는 점을 강조합니다. 또한, 1+1D에서의 결과가 제시되었으나, 기존 문헌은 2D 일반 회로에서 시공간 격자 전체가 "쉬운" 상에 놓일 수 있음을 시사하며, 이는 DTOE가 고차원에서 특히 강력할 수 있음을 의미합니다(다만, 이 일반화는 본 연구에서 명시적으로 테스트되지 않았습니다). 마지막으로, DAOE가 특정 소산 슈퍼오퍼레이터를 가진 MPO 방법으로 프레임화되는 것과 달리, DTOE를 유사하게 프레임화하는 것은 향후 과제로 남겨두는 것이 도전 과제임을 언급하였습니다.