Quantum chaos and late-time equipartition of symmetry-resolved Krylov complexity

기술 요약: 양자 카오스와 대칭성 분해된 크릴로프 복잡도의 후기 시간 등분배

문제 제기
크릴로프 복잡도($C_K$)는 양자 다체계의 연산자 성장(operator growth)을 진단하는 도구로 부상하였으며, 일반적으로 카오스 시스템에서 초기 지수적 성장 이후 선형적 단계와 최종적인 포화 단계로 특징지어진다. 그러나 해밀토니안 $H$와 초기 시드 연산자 $O$가 보존 전하 $P$와 교환(commute)하는 경우, 후기 시간 포화값에 대한 해석은 정확한 대칭성의 존재로 인해 복잡해진다. $H$와 $O$가 동일한 대칭성을 공유할 때, 초기 시간대에는 전체 크릴로프 복잡도가 (초기 연산자의 투영에 의해 결정되는 가중치에 따라) 대칭성별 분해된 복잡도의 가중 평균과 관련이 있음이 확립되어 있다. 그러나 포화 이후의 후기 시간에서 전체 복잡도와 분해된 복잡도 사이의 관계는 여전히 불분명한 상태로 남아 있었다. 단순한 기대치는 포화값이 전체 힐베르트 공간 차원의 제곱($D^2$)에 비례할 것이라는 것이었으나, 대칭성의 존재가 이러한 스케일링을 변화시켜, 만약 적절히 분해되지 않는다면 크릴로프 복잡도가 양자 카오스의 진단 도구로서 갖는 능력을 가릴 수 있다.

방법론
저자들은 유한 차원 카오스 시스템에서 대칭성 분해된 크릴로프 복잡도의 후기 시간 거동을 조사하기 위해 해석적 유도와 수치적 시뮬레이션의 결합을 사용한다.

• 해석적 프레임워크: 본 연구는 리우빌리안 슈퍼 연산자 $\mathcal{L} = [H, \cdot]$에 의해 생성되는 크릴로프 공간 형식을 활용한다. 힐베르트 공간 $\mathcal{H}$를 차원이 $d_q$인 대칭 섹터 $\mathcal{H}_q$들로 분해함으로써, 저자들은 각 섹터 내에서 접근 가능한 크릴로프 공간의 차원($K_q$)과 전체 미분해 공간($K$)을 분석한다. 저자들은 일반적인 카오스 해밀토니안에 대해 (보편적인 제로 모드를 제외하고) 에너지 차이가 비퇴화적이며, 시드 연산자가 모든 대칭 허용 리우빌리안 주파수와 비제로 중첩을 가진다고 가정한다. 이러한 가정 하에, 저자들은 크릴로프 공간의 차원을 유도하고, 크릴로프 파동함수가 접근 가능한 크릴로프 체인 전체에 균일하게 분포한다는 후기 시간 델로컬라이제이션(delocalization) 가설을 적용한다.
• 수치적 구현: 저자들은 크릴로프 복잡도를 계산하기 위해 란초스 알고리즘(수치적 안정성을 위해 아르놀디 반복법 사용)을 구현한다. 저자들은 네 가지 모델에 대해 해석적 예측을 테스트한다:
  1. 실수 SYK 모델: $Z_2$ 페르미온 패리티 대칭성을 가진 마요라나 페르미온 모델.
  2. 복소 SYK 모델: 전역 $U(1)$ 전하 보존을 가진 복소 페르미온 모델.
  3. 보존 스핀 모델: $Z_2 \times Z_2$ 대칭성을 가진 카오스 이차 스핀 모델.
  4. 혼합 필드 이싱 체인: 반사 대칭성을 가진 비적분적 스핀 체인.
     또한, 카오스적 행동과 대조하기 위해 이싱 체인의 적분 가능한 극한을 분석한다.

주요 기여 및 결과
본 논문은 포화 영역($t_{sat} \lesssim t \ll t_{rec}$)에서 전체 크릴로프 복잡도가 대칭 섹터들에 대해 가산적(additive)임을 입증한다:
$$C_K(t) \simeq \sum_q C_K^{(q)}(t)$$
결정적으로, 각 섹터의 분율적 기여는 (초기 시간 영역에서와 같이) 초기 투영 가중치에 의해 결정되거나 균등하게 배분되는 것이 아니라, 동적으로 접근 가능한 연산자 공간의 차원에 의해 결정된다. 구체적으로, 후기 시간 분율 복잡도는 다음을 따른다:
$$\frac{C_K^{(q)}}{C_K} \simeq \frac{d_q(d_q - 1)}{\sum_{q'} d_{q'}(d_{q'} - 1)} \approx \frac{d_q^2}{\sum_{q'} d_{q'}^2}$$
여기서 $d_q$는 섹터 $q$의 힐베르트 공간 차원이다. 이 현상은 **차원 가중 등분배(dimension-weighted equipartition)**라고 명명되었다.

• 수치적 검증:
  • 실수 SYK 모델 ($N=16$)에서, 두 패리티 섹터는 동일한 차원($d_+ = d_-$)을 가지며, 이는 $d_q^2$ 스케일링과 일치하는 균등한 분할($50\%$)을 나타낸다.
  • 복소 SYK 모델 ($N=7, 8$)에서, 전하 섹터들은 서로 다른 차원($d_q = \binom{N}{q}$)을 가진다. 수치 결과는 더 큰 $d_q$를 가진 섹터가 전체 복잡도에 더 큰 분율을 기여함을 보여주며, 이는 $d_q^2$ 가중치를 정확히 따른다.
  • 보존 스핀 모델 ($N=6$)은 네 개의 동일한 섹터를 가지며 균등한 분할($25\%$)을 확인시켜 준다.
  • 혼합 필드 이싱 체인 ($N=8$)은 서로 다른 차원을 가진 반사 섹터($d_+ = 136, d_- = 120$)를 특징으로 한다. 포화값은 $50/50$에서 벗어나 예측된 비율 $d_+^2 / (d_+^2 + d_-^2) \approx 0.56$과 일치한다.
• 적분 가능성 대조: 이싱 체인의 적분 가능한 극한에서, 저자들은 리우빌리안 퇴화와 선택 규칙이 실제 크릴로프 차원($K_q$)을 일반적인 카오스 경계($d_q(d_q-1)+1$)보다 훨씬 낮게 감소시킴을 보여준다. 결과적으로, 후기 시간 분율 복잡도는 $d_q^2$ 스케일링이 아닌 실제 크릴로프 차원의 비율 $(K_q-1)/(K_{full}-1)$에 따라 포화된다. 이는 차원 가중 법칙이 크릴로프 공간이 최대한으로 탐색되는 카오스 역학에 특유한 것임을 강조한다.

의의
본 논문은 카오스적 연산자 성장의 진단으로서 크릴로프 복잡도의 포화값을 올바르게 해석하기 위해 정확한 대칭성을 분해하는 것이 필수적이라고 주장한다. 대칭성을 분해하지 않으면, 포화값이 전체 힐베르트 공간 차원의 제곱($D^2$)에 비례할 것이라고 잘못 가정할 수 있다. 대신, 저자들은 카오스 시스템의 경우 포화가 섹터 차원의 제곱의 합($\sum d_q^2$)에 비례하며, 이는 섹터가 불균등할 경우 $D^2$보다 매개변수적으로 작음을 보여준다.

이 연구는 카오스 시스템에서 전체 크릴로프 복잡도와 대칭성 분해된 크릴로프 복잡도 사이의 "단순하고 일반적인 관계"를 제공하며, 초기 시간 영역(초기 상태 투영에 의해 결정됨)과 후기 시간 영역(접근 가능한 연산자 공간 차원에 의해 결정됨)을 구분한다. 이는 크릴로프 복잡도 포화의 해석을 명확히 하여, 포화값이 대칭성에 의한 블록 구조에 의해 가려지지 않고 동적으로 접근 가능한 부분 공간 내에서의 실제 연산자 성장 정도를 반영하도록 보장한다. 이러한 결과는 대칭성 분해가 카오스 역학의 더 정밀한 특성화를 제공하는 정교화 도구 역할을 함을 시사한다.