Efficient computation of real-time correlators using Pauli Propagation

기술 요약: 파울리 전파를 이용한 실시간 상관 함수의 효율적 계산

문제 정의
실시간 2점 시간 순서 상관 함수, $C(t) = \text{Tr}[\rho A(t)B]$,를 계산하는 것은 다체 물리(many-body physics)의 근본적인 과제로, 미시적인 해밀토니안의 세부 사항과 거시적인 실험적 관측량(예: 중성자 산란, 전도도) 사이를 잇는 가교 역할을 합니다. 전통적인 방법들은 상당한 병목 현상에 직면해 있습니다: 정확한 대각화(exact diagonalization)는 메모리에 의해 제한되고, 실시간 양자 몬테카를로(quantum Monte Carlo)는 부호 문제(sign problem)를 겪으며, 텐서 네트워크 방법(TEBD 또는 TDVP와 같은)은 시간 진화 과정 중 결합 차원(bond dimension)의 급격한 성장으로 인해 1차원을 넘어 확장하는 데 어려움을 겪습니다. 파울리 전파(Pauli propagation, PP)는 하이젠베르크 묘사에서 연산자를 직접 진화시킴으로써 유망한 대안을 제시하지만, 본질적으로 파울리 스트링(Pauli strings)의 지수적 증식에 의해 제한되며, 이는 절단(truncation)을 필연적으로 수반하여 편향(bias)을 도입하고 접근 가능한 시간 창을 제한합니다.

방법론
저자들은 짧은 시간 동안의 파울리 전파와 고급 신호 처리 및 바닥 상태 근사 기술을 결다하여, 직접 시뮬레이션 가능한 범위를 넘어 상관 함수를 확장하는 하이브리드 계산 전략을 제안합니다. 이 방법론은 세 가지 주요 구성 요소로 이루어집니다:

1. 절단된 파울리 전파 (Truncated Pauli Propagation): 시간 진화된 연산자 $A(t)$는 파울리 스트링의 기저로 확장됩니다, $A(t) = \sum_i a_i(t) P_i$. 파울리 스트링의 지수적 성장을 관리하기 위해 저자들은 두 가지 절단 전략을 채택합니다:

   • 계수 컷오프 ($\epsilon$): 계수 $|a_i| < \epsilon$인 스트링을 폐기합니다.
   • 가중치 컷오프 ($w$): 파울리 가중치(비항등 사이트의 수)가 $w$를 초과하는 스트링을 폐기합니다.
     진화는 트로터화(Trotterization)를 통해 수행되며, 각 단계에서 각 항은 해밀토니안과의 교환 관계에 따라 파울리 스트링을 분기(branching)시킵니다.

2. 클리포드 증강 행렬 곱 상태 (Clifford-Augmented Matrix Product States, CAMPS): 얽힌 바닥 상태(특히 2D의 경우)에 대해 기댓값 $\langle \Psi | A(t) | \Psi \rangle$을 평가하는 어려움을 해결하기 위해, 저자들은 CAMPS를 활용합니다. 이 방식은 행렬 곱 상태(MPS)에 클리포드 변환을 적용하여 파울리 합의 계수 구조를 보존하면서 필요한 결합 차원을 크게 줄입니다. 이를 통해 표준 MPS 방법으로는 계산이 매우 까다로운 복잡한 바닥 상태에 대한 기댓값을 효율적으로 평가할 수 있습니다.

3. 양의 준정부호 (Positive Semi-Definite, PSD) 확장: 상관 함수 $C_{AA}(t) = \text{Tr}[\rho A^\dagger(t)A(0)]$가 양의 함수임을 인식하여, 저자들은 이들의 시간 이산 토플리츠 행렬(Toeplitz matrices)이 양의 준정부호(PSD)라는 사실을 활용합니다. Carathéodory, Toeplitz, Pisarenko의 연구에 기반하여, 랭크 $r$인 PSD 토플리츠 행렬은 $2r$개의 노이즈 없는 데이터 포인트에 의해 유일하게 결정됩니다. 저자들은 짧은 시간의 PP 데이터를 "초고해상도(super-resolve)"하여, 추가적인 연산자 진화 없이도 상관 함수를 늦은 시간까지 확장하기 위해 저계수(low-rank) PSD 피팅을 사용합니다. 이는 많은 물리계가 적은 수의 특성 주파수에 의해 지배된다는 관찰에 의존하여 "파울리 폭발" 문제를 효과적으로 완화합니다.

주요 결과
저자들은 1D 및 2D 하이젠버그 모델에 대해 이 결합된 접근 방식의 효능을 입증했습니다:

• 1D 강자성 하이젠버그 사슬 ($N=64$): 저자들은 공격적인 절단(시간 궤적이 정확한 결과와 질적으로만 일치하는 경우)을 적용하더라도, 푸리에 변환된 스펙트럼이 올바른 물리적 특징(예: 마그논 분산)을 유지함을 보여줍니다. ad-hoc 감쇠 함수나 PSD 확장을 사용하면 원시 시간 궤적에서의 장기 역학 정확도 손실에도 불구하고 스펙트럼 특징을 회복할 수 있습니다.
• 1D 반강자성 하이젠버그 사슬 ($N=64$): DMRG로 준비된 바닥 상태를 사용하여, 저자들은 PP를 TDVP와 비교합니다. 공격적인 PP 절단은 MPS-TDVP에 비해 상당한 주파수 넓어짐(broadening)과 해상도 감소를 초래하지만, 이 방법은 스피논 연속체(spinon continuum)의 일반적인 특징을 성공적으로 포착합니다. 이 연구는 얽힌 바닥 상태의 경우, 병목 현상이 연산자 진화에서 기댓값 평가로 이동함을 강조합니다.
• 2D 반강자성 하이젠버그 격자 ($8 \times 8$): 이는 표준 TNS 방법이 너무 비싸지는 영역을 나타냅니다. PP를 CAMPS(결합 차원 $D=240$)와 결격-1 PSD 확장과 결합함으로써, 저자들은 단 7번의 타임 스텝만으로 동적 스핀 구조 인자 $S(k, \omega)$를 계산합니다. 결과 스펙트럼은 선형 스핀파 이론(linear spin-wave theory)과 잘 일치하며, $\Gamma$에서의 골드스톤 모드, $X$로의 분산, 그리고 $M$에서의 연화(softening)를 정확히 식별합니다.
• 확장성: 저자들은 2D의 경우, 제한 사항이 (단기 확장에 의해 피할 수 있는) 늦은 시간에서의 무제한적인 연산자 성장이라기보다, 바닥 상태의 준비와 기댓값의 평가에 있음을 언급합니다. 또한 부록에서 이 방법이 $100 \times 100$ 강자성 격자의 마그논 스펙트럼을 회복할 수 있음을 추가로 보여줍니다.

의의 및 주장
이 논문은 결합된 접근 방식이 클래식 방법이 상호작용하는 양자 다체계의 실시간 역학을 신뢰성 있게 탐사할 수 있는 영역을 넓힌다고 주장합니다.

• 지수적 성장 극복: 이 방법은 연산자의 지수적 증식을 피하기 위해 직접 시뮬레이션을 관리 가능한 수준의 짧은 시간 창으로 제한한 뒤, 물리적 제약(양의 성질 및 저계수 주파수 지배)을 사용하여 데이터를 확장합니다.
• 고차원에서의 효율성: CAMPS를 통합함으로써, 저자들은 평형 상관 함수를 표준 MPS 시간 진화와 관련된 가파른 결합 차원 성장 없이도 2D 시스템에 대해 계산 가능한 경계에서 계산할 수 있음을 보여줍니다.
• 스펙트럼 충실도: 결과는 시스템의 역학이 제한된 수의 주파수에 의해 지배되는 경우, 극도로 제한된 시간 데이터(2D의 경우 단 7개 타임 스텝)로부터 높은 해상도의 스펙트럼 정보를 추출할 수 있음을 나타냅니다.

저자들은 이 방법의 한계를 겸허히 인정합니다: 저계수 가정은 연속적인 주파수를 가진 시스템(예: 1D 스피논 연속체)에서는 실패하며, 방법의 성공 여부는 절단이 기저에 깔린 주파수를 식별할 만큼 충분히 정확한지에 달려 있습니다. 저자들은 이 작업을 1D에서의 텐서 네트워크를 대체하는 것이 아니라, 현재의 최첨단 방법들이 어려움을 겪는 대규모 고차원 격자의 동적 상관 함수를 계산하기 위한 유망하고 실행 가능한 경로로 위치시킵니다.