Anomalous entanglement scaling from eigenvector nonorthogonality in critical non-Hermitian free fermions

기술 요약: 고유벡터의 비직교성으로부터 기인하는 비정상적 얽힘 스케일링: 임계 비헤르미안 자유 페르미온 연구

문제 정의
얽힘 엔트로피는 양자 상(phase)과 임계점을 진단하는 보편적인 도구이다. 1차원 헤르미안 임계 계에서 부분계의 길이 $\ell$에 대한 본 네이만 얽힘 엔트로피 $S_1$은 $S_1 \sim (c/3) \ln \ell$와 같이 로그 함수적으로 스케일링되며, 여기서 계수는 기저가 되는 공형 장론(CFT)의 중심 전하 $c$에 의해 결정된다. 그러나 개방형 양자계, 이득-손실 과정, 그리고 사후 선택된 측정 궤적에서 발생하는 비헤르미안 시스템은 고유값(exceptional points) 및 비헤르미안 스킨 효과와 같은 독특한 현상을 보여준다. 기존 연구들은 음의 중심 전하를 갖거나 파라미터에 의존하는 계수를 갖는 등 다양한 거동을 보여주었으나, 이러한 이상 현상을 설명할 수 있는 통합된 메커니즘은 비헤르미안 자유 페르미온 정상 상태(steady states)에 대해 결여되어 있었다.

방법론
저자들은 임계 1차원 비헤르미안 자유 페르미온 체인, 특히 비헤르미안 SSH(Su–Schrieffer–Heiger) 모델에 초점을 맞추어 정상 상태를 조사한다. 정상 상태 $|\Psi_S\rangle$는 에너지의 허수 부분($\text{Im } E_\alpha$)이 가장 큰 우측 고유상태(right eigenstate) $|R_\alpha\rangle$를 선택하는 비유니타리 진화의 장기 극한으로 정의된다.

1. 모델 및 정상 상태 정의: 시스템은 비대칭 호핑(허수 게이지 장)을 가진 이차 형식 해밀토니안으로 기술된다. 정상 상태는 $\text{Im } \varepsilon_m > 0$인 모든 단일 입자 우측 고유 모드(right eigenmodes)를 채움으로써 형성되는 슬레이터 행렬식(Slater determinant)에 해당한다.
2. 얽임 계산: 레니(Rényi) 얽힘 엔트로피 $S_n$은 상관 행렬법을 사용하여 계산된다. 비헤르미안 시스템의 경우, 저자들은 상관 행렬이 일반적으로 직교하지 않는 점유된 우측 고유 벡터들로부터 구성되어야 함을 강조한다. 이를 위해 물리적 밀도 행렬 $\rho = |\Psi_S\rangle\langle\Psi_S|$를 정의하기 위한 직교화 단계(QR 분해 또는 슈어 분해를 통한)가 필요하다.
3. 해석적 유도: 저자들은 "점유 경계"(에너지의 허수 부분이 0이 되는 운동량 지점) 근처에서 상관 행렬의 저에너지 확장을 수행한다. 이들은 이 경계를 가로지르는 투영 연산자(projector)의 불연속성을 분석한다.
4. 수치적 검증: 다양한 파라미터 세트와 무질서 강도에 대해 얽힘 엔트로피를 계산하기 위해 대규모 격자 수치 계산($L \le 65536$)을 수행하였으며, 피팅된 로그 계수를 해석적 예측과 비교하였다.

주요 기여 및 결과

• 이상적 스케일링과 "허수" 디락 점(Imaginary Dirac Point):
  본 연구는 얽힘 엔트로피가 로그 스케일링($S_n \sim \zeta_n \ln \ell$)을 따르지만, 계수 $\zeta_n$이 헤르미안의 의미에서 보편적이지 않음을 밝혀냈다. 대신 $\zeta_n$은 시스템 파라미터에 따라 연속적으로 변하며, 단일 중심 전하로 재현될 수 없는 계열을 형성한다.
  저자들은 이 이상 현상의 기원을 "허수" 디락 점에서 찾았다. 이는 에너지의 허수 부분($\text{Im } E(k) = 0$)에서의 교차점으로, 이 지점에서 실수 부분의 밴드는 분리된 상태를 유지한다. 이 지점에서 시스템은 두 블로흐 상태 사이를 전이한다.

• 고유벡터 비직교성의 역할:
  헤르미안 임계점에서 교차하는 고유벡터들은 서로 직교하여 점유 투영 연산자의 최대 불연속성을 유도한다. 반면 비헤르미안 경우, 허수 디락 점에서의 우측 고유벡터들은 비직교적이다. 이 비직교성은 점유된 스피너(spinor)들의 내적에 의해 정의되는 중첩각(overlap angle) $\theta$에 의해 정량화된다.
  이러한 비직교성은 점유 투영 연산자의 불연속성을 1에서 $\cos \theta$로 약화시킨다. 이 "약화된 불연속성"은 얽힘 스펙트럼에 갭을 생성하고 로그 계수 $\zeta_n$을 감소시킨다.

• 폐쇄형 계수(Closed-Form Coefficients):
  저자들은 중첩각 $\theta$의 함수로서 레니 계수 $\zeta_n(\theta)$에 대한 폐쇄형 식을 유도하였다.
  $$\zeta_n = \frac{2}{\pi^2} \int_{E_\theta}^\infty d\epsilon \frac{\cos \theta \tanh(\epsilon/2)}{\sqrt{\tanh^2(\epsilon/2) - \sin^2 \theta}} h_n\left(\frac{1}{1+e^{-\epsilon}}\right)$$
  여기서 $E_\theta = 2 \text{arctanh}|\sin \theta|$이다. 이 해석적 결과는 전체 갭리스(gapless) 영역과 모든 레니 지수에 대해 격자 수치 계산과 매우 우수한 정량적 일치를 보인다.

• 가법성(Additivity) 및 일반화:
  이 프레임워크는 임의의 1차원 임계 비헤르미안 정상 상태로 일반화된다. 전체 로그 계수는 다음과 같이 가법적임을 보였다:
  $$\zeta_n = \sum_{i=1}^m \zeta_n(\theta_i) + N_F \frac{1 + 1/n}{12}$$
  여기서 $m$은 허수 디락 점의 개수(중첩각 $\theta_i$를 가짐)이고, $N_F$는 단일 밴드 페르미 점의 개수(표준 키랄 페르미온 값을 기여함)이다. 이 공식은 $\text{Re } E(k) = 0$의 교차로 정의되는 비헤르미안 바닥 상태에도 적용된다.

• 무질서에 대한 견고성:
  중요한 발견 중 하나는 약한 실수 온사이트 무질서(real onsite disorder)에 대한 이 로그 스케일링의 견고성이다. 헤르미안 무질서 자유 페르미온이 일반적으로 앤더슨 국소화(area-law) 상이나 랜덤 싱글렛 고정점으로 흐르는 것과 달리, 비헤르미안 정상 상태는 로그 스케일링을 유지한다. 약한 무질서는 유효 중첩각 $\theta_{\text{eff}}$를 재규격화하여 얽힘 계수를 증가시키지만, 스케일링 형태 자체를 파괴하지는 않는다. 저자들은 이를 각 복소 에너지가 단일 전파 채널에 의해 지원되기 때문에 복소 에너지 스펙트럼에서 불순물 후방 산란(backscattering)이 존재하지 않는다는 점에 기인한다고 설명한다.

의의
본 논문은 임계 비헤르미안 자유 페르미온 정상 상태에서의 얽힘 스케일링에 대한 일반적인 메커니즘을 제공한다. 연구는 로그 계수의 비보편성이 경계 효과나 특정 모델의 특수한 특징이 아니라, 허수 디락 점에서의 고유벡터 비직교성에서 기인하는 벌크 효과임을 입증하였다. 이는 비헤르미안 임계성을 헤르미안 임계성(고정된 $c$) 및 다른 비헤르미안 시나리오(예: 측정 유도 전이)와 구별 짓는다. 결과적으로 고유벡터 번들의 기하학적 성질(중첩각)을 얽힘 엔트로피와 직접 연결함으로써 얽힘에 대한 통합된 분류 체계를 제시하며, 이러한 스케일링이 헤르미안 무질서 시스템의 거동과 극명하게 대조되는 방식으로 약한 무질서에 대해 견고함을 보여준다.