We study the statistical properties of the spread complexity in the Krylov space of quantum systems driven across a quantum phase transition. Using the diabatic Magnus expansion, we map the evolution to an effective one-dimensional hopping model. For the transverse field Ising model, we establish an exact link between the growth of complexity and the Kibble-Zurek defect scaling: all cumulants of complexity exhibit the same power-law scaling as the defect density, with coefficients identical to the mean, and the full distribution asymptotically becomes Gaussian. These results yield general scaling arguments for the growth of complexity across arbitrary second-order quantum phase transitions.
논문 ID : 2510.13947제목 : Universal Growth of Krylov Complexity Across A Quantum Phase Transition저자 : András Grabarits, Adolfo del Campo분류 : quant-ph발표 시간 : 2025년 10월 17일논문 링크 : https://arxiv.org/abs/2510.13947 본 논문은 양자 상전이 과정에서 크릴로프 공간 내 확산 복잡도의 통계적 성질을 연구한다. 비단열 매그누스 전개를 사용하여 진화를 유효한 1차원 점프 모델로 매핑한다. 횡장 이싱 모델의 경우, 복잡도 성장과 Kibble-Zurek 결함 스케일링 간의 정확한 연결을 확립한다: 복잡도의 모든 누적량은 결함 밀도와 동일한 멱법칙 스케일링을 나타내며, 계수는 평균값과 동일하고, 완전한 분포는 점근적으로 가우스 분포로 수렴한다. 이러한 결과는 임의의 2차 양자 상전이에서의 복잡도 성장에 대한 일반적인 스케일링 논증을 제공한다.
해결해야 할 문제 :양자 다체 동역학에서 복잡도 성장은 열화, 혼돈 및 양자 진화의 계산 난이도를 설명하는 기초적 문제 시간 의존 설정에서 연산자 성장 및 크릴로프 복잡도의 특성화는 여전히 기본적으로 미개발 상태 복잡도와 같은 물리량이 비평형 보편성을 나타내는지 여부를 결정할 필요성 문제의 중요성 :크릴로프 부분공간 방법은 양자 복잡도 성장의 통용 프레임워크가 됨 양자 제어 및 변분 양자 알고리즘에서 실용적 가치 보유 복잡도 성장과 임계 양자 동역학의 보편적 특징 연결 기존 방법의 한계 :기존 시간 의존 설정의 일반 프레임워크는 비국소 시간의 Floquet 연산자를 사용하여 다체 시스템에 적용 어려움 복잡도 성장과 양자 상전이 보편성을 연결하는 이론적 프레임워크 부재 연구 동기 :시간 의존 설정에 적용 가능한 정확하고 통용적인 분석 프레임워크 개발 크릴로프 복잡도와 Kibble-Zurek 메커니즘 간의 연결 확립 임계점 근처의 복잡도 성장에 대한 보편적 프레임워크 제공 구동 양자 시스템의 시간 의존 설정에 적용 가능한 비단열 매그누스 연산자 기반의 크릴로프 부분공간 새로운 방법 제시 복잡도 성장과 Kibble-Zurek 결함 스케일링 간의 정확한 연결 확립 , 횡장 이싱 모델에서 모든 복잡도 누적량이 동일한 멱법칙 스케일링을 따름을 증명KZ 스케일링 구간 내 크릴로프 복잡도 통계가 보편적 가우스 분포로 수렴함을 증명 , 완전한 분포가 보편적 행동을 나타냄임의의 2차 양자 상전이에 적용 가능한 일반적 스케일링 논증 제공 , 결과를 더 광범위한 임계 시스템 범주로 일반화유한 속도로 양자 상전이를 통과하는 양자 시스템의 크릴로프 복잡도의 통계적 성질을 연구하고, 복잡도 성장과 상전이 동역학 간의 보편적 연결을 확립한다.
구동 시스템 H ( t ) = ∑ n E n ( t ) ∣ n ( t ) ⟩ ⟨ n ( t ) ∣ H(t) = \sum_n E_n(t)|n(t)\rangle\langle n(t)| H ( t ) = ∑ n E n ( t ) ∣ n ( t )⟩ ⟨ n ( t ) ∣ 에 대해, 비단열 시간 진화 연산자를 정의한다:
U ( t ) = U ( t ) U a d ( t ) † U(t) = U(t)U_{ad}(t)^\dagger U ( t ) = U ( t ) U a d ( t ) †
여기서 U a d = ∑ n ∣ n ( t ) ⟩ ⟨ n ( 0 ) ∣ U_{ad} = \sum_n |n(t)\rangle\langle n(0)| U a d = ∑ n ∣ n ( t )⟩ ⟨ n ( 0 ) ∣ 는 평행 수송 연산자이다. 비단열 매그누스 연산자는 다음과 같이 정의된다:
Ω ( t ) = i log ( U ( t ) ) \Omega(t) = i\log(U(t)) Ω ( t ) = i log ( U ( t ))
시간 진화 상태는 다음과 같이 표현될 수 있다:
∣ ψ ( t ) ⟩ = e − i Ω ( t ) ∣ G S ( t ) ⟩ = ∑ m = 0 ∞ ( − i ) m Ω m ( t ) m ! ∣ G S ( t ) ⟩ |\psi(t)\rangle = e^{-i\Omega(t)}|GS(t)\rangle = \sum_{m=0}^{\infty} \frac{(-i)^m \Omega^m(t)}{m!}|GS(t)\rangle ∣ ψ ( t )⟩ = e − i Ω ( t ) ∣ GS ( t )⟩ = ∑ m = 0 ∞ m ! ( − i ) m Ω m ( t ) ∣ GS ( t )⟩
해당하는 Lanczos 알고리즘은 다음과 같다:
∣ K n + 1 , t ⟩ b n + 1 , t = Ω ( t ) ∣ K n , t ⟩ − a n , t ∣ K n , t ⟩ − b n , t ∣ K n − 1 , t ⟩ |K_{n+1,t}\rangle b_{n+1,t} = \Omega(t)|K_{n,t}\rangle - a_{n,t}|K_{n,t}\rangle - b_{n,t}|K_{n-1,t}\rangle ∣ K n + 1 , t ⟩ b n + 1 , t = Ω ( t ) ∣ K n , t ⟩ − a n , t ∣ K n , t ⟩ − b n , t ∣ K n − 1 , t ⟩
여기서:
a n , t = ⟨ K n , t ∣ Ω ( t ) ∣ K n , t ⟩ a_{n,t} = \langle K_{n,t}|\Omega(t)|K_{n,t}\rangle a n , t = ⟨ K n , t ∣Ω ( t ) ∣ K n , t ⟩ b n , t = ⟨ K n − 1 , t ∣ Ω ( t ) ∣ K n , t ⟩ b_{n,t} = \langle K_{n-1,t}|\Omega(t)|K_{n,t}\rangle b n , t = ⟨ K n − 1 , t ∣Ω ( t ) ∣ K n , t ⟩ TFIM 해밀토니안:
H ^ ( t ) = − J ∑ j = 1 L [ σ j z σ j + 1 z + g ( t ) σ j x ] \hat{H}(t) = -J\sum_{j=1}^L [\sigma_j^z \sigma_{j+1}^z + g(t)\sigma_j^x] H ^ ( t ) = − J ∑ j = 1 L [ σ j z σ j + 1 z + g ( t ) σ j x ]
페르미온 표현에서, 각 운동량 모드는 독립적인 2에너지 준위 시스템에 대응된다:
H ^ = 2 ∑ k ψ ^ k † H k ( t ) ψ ^ k \hat{H} = 2\sum_k \hat{\psi}_k^\dagger H_k(t) \hat{\psi}_k H ^ = 2 ∑ k ψ ^ k † H k ( t ) ψ ^ k
여기서 H k ( t ) = ( g ( t ) − cos k ) τ z + sin k τ x H_k(t) = (g(t) - \cos k)\tau_z + \sin k \tau_x H k ( t ) = ( g ( t ) − cos k ) τ z + sin k τ x 이다.
시간 국소성 : 서로 다른 시간의 크릴로프 기저와 Lanczos 계수는 독립적으로 생성되어 시간 비국소성 문제를 회피한다정확한 해석 결과 : KZ 스케일링 구간 내에서 Lanczos 계수 및 크릴로프 파동함수의 정확한 표현식을 획득한다:
b n ∼ L 1 / 2 τ − 1 / 4 n , a n ∼ L τ − 1 / 2 b_n \sim L^{1/2}\tau^{-1/4}\sqrt{n}, \quad a_n \sim L\tau^{-1/2} b n ∼ L 1/2 τ − 1/4 n , a n ∼ L τ − 1/2 포아송 통계 : 크릴로프 파동함수는 주도 차수 근사에서 포아송 분포로 나타난다:
∣ ϕ n ( τ ) ∣ 2 ∼ ( 2 C L τ − 1 / 2 ) n e − 2 C L τ − 1 / 2 n ! |\phi_n(\tau)|^2 \sim \frac{(2CL\tau^{-1/2})^n e^{-2CL\tau^{-1/2}}}{n!} ∣ ϕ n ( τ ) ∣ 2 ∼ n ! ( 2 C L τ − 1/2 ) n e − 2 C L τ − 1/2 횡장 이싱 모델, 시스템 크기 L = 200 , 400 L = 200, 400 L = 200 , 400 유효 시스템 크기 L e f f = 46 , 48 L_{eff} = 46, 48 L e ff = 46 , 48 (Lanczos 알고리즘용) 구동 시간 범위: τ = 2 \tau = 2 τ = 2 에서 τ = 750 \tau = 750 τ = 750 크릴로프 복잡도의 처음 세 누적량: K 1 K_1 K 1 (평균값), K 2 K_2 K 2 (분산), K 3 K_3 K 3 (3차 중심 모멘트) 복잡도 분포의 가우스 적합 품질 KZ 스케일링 τ − 1 / 2 \tau^{-1/2} τ − 1/2 와의 일치도 선형 구동 프로토콜: g ( t ) = t / τ g(t) = t/\tau g ( t ) = t / τ 상자성 상에서 초기화: t = − g 0 τ t = -g_0\tau t = − g 0 τ g ( 0 ) = 0 g(0) = 0 g ( 0 ) = 0 에서 종료정확한 수치 시뮬레이션 계산 복잡도 제한: L ≈ 48 L \approx 48 L ≈ 48 누적량의 보편적 스케일링 :
모든 복잡도 누적량은 KZ 스케일링 구간 내에서 동일한 멱법칙을 따른다:
K q ≈ 2 C L τ − 1 / 2 K_q \approx 2CL\tau^{-1/2} K q ≈ 2 C L τ − 1/2 이는 결함 밀도의 KZ 스케일링 n ∼ τ − 1 / 2 n \sim \tau^{-1/2} n ∼ τ − 1/2 와 완전히 일치한다.가우스 분포 수렴 :
KZ 스케일링 구간 L τ − 1 / 2 ≫ 1 L\tau^{-1/2} \gg 1 L τ − 1/2 ≫ 1 내에서, 복잡도 분포는 가우스 분포로 수렴한다:
P K ( n ) ≈ exp [ − ( n − K 1 ) 2 / ( 2 K 2 ) ] 2 π K 2 P_K(n) \approx \frac{\exp[-(n-K_1)^2/(2K_2)]}{\sqrt{2\pi K_2}} P K ( n ) ≈ 2 π K 2 e x p [ − ( n − K 1 ) 2 / ( 2 K 2 )] Lanczos 계수 스케일링 :
수치 결과는 이론적으로 예측된 스케일링 행동을 검증하며, 모든 서로 다른 τ \tau τ 와 L L L 의 데이터는 적절한 재스케일링 후 단일 곡선으로 붕괴된다.임계점 근처에서 복잡도가 급격히 증가하여 양자 상전이와 복잡도 성장의 상호작용을 포착한다 임계점에서 멀어질 때 비보편적 진동 행동을 나타낸다 대칭 깨짐 상 내에서 점진적으로 점근 값으로 수렴한다 그림 1은 서로 다른 시스템 크기 및 구동 시간의 비대각 Lanczos 계수를 보여주며, 적절한 재스케일링 τ 1 / 4 L − 1 / 2 \tau^{1/4}L^{-1/2} τ 1/4 L − 1/2 후 완벽하게 단일 곡선으로 붕괴된다.
Parker 등이 제시한 보편적 연산자 성장 가정 양자 혼돈 및 얽힘 성장에서의 크릴로프 부분공간 방법 적용 시간 무관 설정에서의 복잡도 연구 2차 양자 상전이에서 결함 형성의 보편적 프레임워크 다양한 양자 시스템에서 KZ 스케일링 검증 결함 통계의 보편성 단열-충동 근사 이론 다체 시스템에서의 Landau-Zener 전이 비평형 양자 상전이의 보편성 복잡도와 상전이 동역학 간의 보편적 연결 확립 : 크릴로프 복잡도의 모든 누적량은 결함 밀도와 동일한 KZ 스케일링을 따른다복잡도 통계의 보편성 증명 : KZ 구간 내에서 복잡도 분포는 보편적 가우스 형태로 수렴한다일반적 프레임워크 제공 : 결과를 ( d − D ) (d-D) ( d − D ) 차원 임계면을 가진 d d d 차원 자유 페르미온 시스템으로 일반화한다일반적 임계 시스템의 경우, 복잡도 누적량은 다음을 따른다:
K q ∼ 2 C L d − D τ − α ( d − D ) K_q \sim 2CL^{d-D}\tau^{-\alpha(d-D)} K q ∼ 2 C L d − D τ − α ( d − D )
여기서 α \alpha α 와 d − D d-D d − D 는 각각 여기 확률과 임계면의 차원을 특징짓는다.
계산 복잡도 제한 : 매그누스 연산자의 직합 형태는 지수 증가하는 계산 비용을 초래하여 정확한 수치 시뮬레이션의 시스템 크기를 제한한다근사의 유효 범위 : 주도 차수 근사는 n ≪ L τ − 1 / 2 n \ll L\tau^{-1/2} n ≪ L τ − 1/2 일 때만 유효하다모델 특이성 : 상세 분석은 주로 TFIM을 대상으로 하며, 다른 모델은 추가 검증이 필요하다실험 검증 : 양자 시뮬레이터 및 양자 어닐링 장치에서 이론적 예측 검증상호작용 시스템 : 강상호작용 양자 시스템으로 일반화고차 상전이 : 1차 상전이 및 위상 상전이에서의 복잡도 행동 연구이론적 엄밀성 : 매그누스 연산자 구성에서 크릴로프 기저 생성까지 완전한 해석 프레임워크와 엄격한 수학적 유도 제공보편성 발견 : 복잡도 성장과 알려진 보편 현상(KZ 메커니즘) 간의 깊은 연결 확립으로 중요한 이론적 의의 보유방법론적 혁신성 : 비단열 매그누스 연산자 방법은 시간 의존 시스템에서 크릴로프 복잡도의 기술적 난제를 해결한다결과의 완전성 : 평균값 행동뿐 아니라 분포의 통계적 성질을 완전히 설명한다수치 검증 범위 제한 : 계산 복잡도 제한으로 인해 수치 검증은 주로 상대적으로 작은 시스템에 국한된다모델 의존성 : 일반적 논증이 제시되었지만 상세 검증은 주로 TFIM이라는 특수 모델에 기반한다실험 실현 가능성 : 이론적 예측의 실험 검증은 기술적 도전에 직면할 수 있다이론적 기여 : 양자 복잡도 이론과 비평형 양자 물리의 교차 분야에 중요한 기여방법론적 가치 : 제시된 분석 프레임워크는 다른 시간 의존 양자 시스템에 적용 가능하다실용적 전망 : 양자 계산 및 양자 시뮬레이션에서 잠재적 응용 가치 보유양자 상전이 연구 : 다양한 양자 상전이에서의 복잡도 행동 연구를 위한 이론적 도구 제공양자 시뮬레이션 : 냉원자 및 초전도 양자비트 시스템에서 이론적 예측 검증양자 알고리즘 설계 : 크릴로프 방법 기반 양자 알고리즘 최적화에 대한 지침 제공논문은 37편의 중요 문헌을 인용하며, 크릴로프 복잡도 이론, Kibble-Zurek 메커니즘, 양자 임계 동역학 등 관련 분야의 핵심 연구를 포함하여 연구에 견고한 이론적 기초를 제공한다.
본 논문은 양자 복잡도 이론과 비평형 양자 물리의 교차 분야에서 중요한 기여를 하였으며, 복잡도 성장과 양자 상전이 동역학 간의 보편적 연결을 확립하여 양자 시스템에서의 복잡도 증가 메커니즘을 이해하기 위한 새로운 이론적 관점을 제공한다.