제3 양자화와 4n 곱하기 4n 행렬을 이용한 개방 페르미계 마스터 방정식 해법

이 논문은 개방 양자 시스템의 동역학을 기술하는 Lindblad 마스터 방정식을, 특히 2n개의 마요라나 연산자로 구성된 이차 페르미계에 대해 완전 해석적으로 푸는 방법을 제시한다. 저자는 먼저 시스템의 Hamiltonian H와 Lindblad 연산자 Lμ를 w·H·w와 lμ·w 형태로 표현한다. 여기서 wj는 반반대칭 마요라나 연산자이며, {wj,wk}=2δjk 를 만족한다. 이러한 설정 하에 Liouvillian L̂은 adjoint Fock 공간 K(차원 4n) 위에서 a‑fermion 연산자 ĉj, ĉj† 로 구성된 quadratic form 으로 재작성된다. 구체적으로 L̂=−4i ĉ†·H·ĉ+2ĉ†·M·ĉ+2ĉ†·(M−Mᵀ)·ĉ† 로, M은 Lindblad 연산자들의 복소 계수로 만든 Hermitian 행렬이다. 다음 단계에서는 a‑fermion을 실수 마요라나 연산자 âr 로 변환해 L̂을 4n 차원의 복소 반대칭 행렬 A 로 표현한다(L̂=â·A·â−A₀·1). A는 일반적인 반대칭 행렬이므로 고유값은 ±βj 쌍으로 나타난다. 저자는 A를 정규형으로 대각화하고, 새로운 정상 모드(bj, bj′)를 정의한다. 이때 Liouvillian은 ‑2∑j βj bj′bj 형태가 되며, 이는 ‘제3 양자화’라 부르는 표준 형태이다. βj는 ‘rapidity’라 불리며, 실부분이 양수이면 시스템은 안정적이다. 또한 βj들의 합은 2 tr M 와 동일한 sum rule 를 만족한다(∑βj=tr M). NESS는 모든 bj가 진공 상태인 |NESS⟩ 로, Liouvillian의 영 고유값에 대응한다. 전체 스펙트럼은 2^{2n}개의 이진 점유수 νj∈{0,1} 로 라벨링된 고유벡터와 λν=−2∑j βj νj 로 완전히 기술된다. 따라서 물리량(예: 두점 상관함수, 전류, 스핀 밀도 등)은 βj와 대응하는 고유벡터만 알면 즉시 계산할 수 있다. 이론적 틀을 검증하기 위해 저자는 근접 이웃 XY 스핀 ½ 체인에 적용한다. XY Hamiltonian은 마요라나 연산자로 변환 가능하고, 체인의 양 끝에 서로 다른 온도와 자기장을 가진 Lindblad 바스를 연결한다. 이 설정은 열 및 스핀 전송 문제를 비평형 정상 상태(NESS)에서 다루게 만든다. 해석적 대각화와 수치적 시뮬레이션을 통해 다음과 같은 주요 결과를 얻는다. 첫째, 무질서가 없는 균일 체인에서는 전도율이 체인 길이 L에 대해 L⁰에 비례하는 ballistic 전도가 나타난다. 이는 전류가 체인 전체에 걸쳐 일정함을 의미한다. 둘째, Liouvillian의 스펙트럼 갭 Δ는 L^{−3} 로 스케일링되며, 이는 긴 체인에서 이완 시간이 L³에 비례함을 보여준다. 셋째, 긴 체인에서 에너지와 스핀 프로파일은 바스 접촉점 근처에서 지수적으로 감소하는 evanescent 모드 형태를 보이며, 이는 전송이 바스와 체인 사이에서 급격히 소멸함을 의미한다. 넷째, 마요라나 결합 상수를 무작위로 변동시키는 디스오더를 도입하면 βj의 실부분이 모두 양의 최소값을 갖게 되어 전도율이 급격히 감소하고, 전반적으로 절연 거동을 보인다. 이러한 절연 현상은 온도와 외부 자기장에 관계없이 나타난다. 제3 양자화 방법은 기존의 직접적인 수치 해법(예: 차원 2^{2n}인 밀도 행렬 직접 대각화)보다 훨씬 큰 시스템(N≈10⁴)까지 적용 가능하게 하며, 비평형 양자 통계역학, 양자 열전달, 그리고 양자 정보 전송 등 다양한 분야에 활용될 수 있다. 또한, 시간 의존성 Hamiltonian이나 비Markovian 확장에도 직관적으로 일반화될 수 있음을 언급한다. 논문은 이와 같은 일반적인 프레임워크가 quadratic open fermionic 시스템을 완전하게 해석하고, 물리적 현상을 정량적으로 예측하는 데 강력한 도구가 됨을 강조한다.