경계 구동 XXZ 스핀 체인의 비평형 정상상태와 전류 통계

초록

본 논문은 Lindblad 마스터 방정식으로 기술되는 개방형 Heisenberg XXZ 스핀 체인의 비평형 정상상태(NESS)를 정확히 구하고, 그 상태에서의 스핀 프로파일, 전류 및 전류 전수 통계(full counting statistics)를 계산한다. 양자 적분 가능성, Lax 연산자와 Cholesky 형태의 밀도 연산자를 이용해 최대 구동(maximally driven) 경우의 해를 얻으며, 약한 시스템‑밥 결합 하에서 전류 누적량의 누적 생성함수(cumulant generating function)를 전개한다.

상세 분석

이 연구는 두 개의 무한한 열원/열극에 의해 경계가 구동되는 1차원 XXZ 스핀 체인을 대상으로 한다. 시스템‑밥 상호작용을 Born‑Markov 및 회전파 근사 하에 평균화하여 Lindblad 형태의 마스터 방정식(2)을 도입하고, 고정점 ρ∞(3)을 찾는다. 핵심은 양자 군 Uq(sl₂) 대수와 그에 대응하는 Lax 연산자 L(φ,s) (8)를 이용해 전역적인 양자 적분 가능성을 확보하는 것이다. Sutherland 방정식(9)은 Lax 연산자와 Hamiltonian hXXZ 사이의 교환 관계를 제공하며, 이를 통해 경계 조건을 만족하는 스펙트럼 파라미터 φ와 연속 스핀 파라미터 s를 결정한다(14‑15). 특히 최대 구동(ε) 상황에서 φ=0, s가 순수 허수값을 갖는 해가 유일하게 존재한다는 점이 중요한 결과다.

정상상태 밀도 연산자는 ρ∞=SS† 형태의 Cholesky 분해(12‑13)로 표현되며, 여기서 S는 Lax 연산자를 n번 곱한 후 최저 상태에 투사한 것이다. 이 구조 덕분에 관측량 ⟨O⟩=tr(Oρ∞)/tr ρ∞(17)를 효율적으로 계산할 수 있다. 특히 스핀 프로파일 ⟨σz_j⟩와 전류 ⟨j_z⟩에 대해, 물리적 공간에 대한 부분 트레이스를 이용해 전이 행렬 T, V, W를 정의하고(19‑21), 이들 사이의 대수적 관계(23)와 연속 방정식으로부터 전류 식(24)를 도출한다.

양자 이방성 Δ=cosγ에 따라 보조 공간의 차원이 유리수 γ=πl/m 형태일 때는 T와 V가 유한 차원(m+1)으로 축소돼 정확한 고유값 분석이 가능하다. 예를 들어 Δ=½(γ=π/3)에서는 스핀 프로파일이 평탄하고 전류가 배럴리스트(∝ε⁴)임을 확인한다. 반면 |Δ|>1(easy‑axis) 영역에서는 전류가 지수적으로 감소해 절연성질을 보인다.

전류 전수 통계는 Lindblad 점프 연산자를 ±1 스핀 변화를 일으키는 두 종류로 분리하고, 카운팅 필드 χ를 도입해 수정된 점프 초연산자 D_jump^χ(29)를 정의한다. 이때 가장 큰 실부를 갖는 고윳값 λ(χ)가 누적 생성함수가 된다(30‑31). 약한 결합 ε에 대해 퍼트urbation 전개(34)를 수행하면, 0차에서는 평형형 ρ(0)=2^{-n}∏_j(1+ν(χ)σz_j) 형태가 나오고, 1차에서 λ^(1)(χ)와 ν(χ)를 구한다(38). 대칭 구동(μ) 경우에는 ν=0, ρ(0)=2^{-n}I가 되며 λ^(1)(χ)=−1+cosχ−iμ sinχ가 얻어진다. 따라서 짝수·홀수 차수 누적량이 각각 ε/2, εμ/2로 단순화된다.

XXZ 체인에 대해 3차 교정 λ^(3)(χ)를 구하기 위해 최대 구동 해의 파라미터 미분으로 정의된 연산자 Z를 도입하고,