Cell Natural Orbitals in Quantum Materials

기술 요약: 양자 물질에서의 셀 자연 궤도함수 (Cell Natural Orbitals)

문제 정의
상관된 양자 물질을 정확하게 모델링하기 위해서는 관련 저에너지 단일 입자 자유도를 설명하는 최소한의 기술이 필요하다. 운동량 공간 기술(블로흐 파동함수)은 분산과 대칭성을 포착하지만, 국소적 전자 상호작용을 다루기 위해서는 실공간 모델이 필수적이다. 그러나 대상 밴드의 분산, 대칭성, 그리고 양자 기하학(특히 단위 셀 내 전하 밀도 및 쌍극자 행렬 요소)을 충실히 재현하는 최소한의 국소 궤도함수 세트를 식별하는 것은 쉽지 않은 과제이다. 최대 국소화된 와니에 함수(MLWF)나 디스엔탱글먼트(disentanglement) 기법과 같은 기존 방법들은 종종 궤도함수의 수를 미리 지정해야 하거나 물리적으로 동기 부여된 초기 투영(projection)에 의존하며, 밴드 자체에 의해 결정되는 최소한의 국소 힐베르트 공간을 최적화하지 못한다. 또한, 위상적 장애(예: Chern 밴드에서의 Chern number) 또는 폐쇄된 원자 한계(obstructed atomic limits)는 표준적인 접근 방식으로는 대칭적이고 국소화된 기저가 대상 매니폴드를 생성하는 것을 방해할 수 있다.

방법론
저자들은 대상 블로흐 밴드에 대한 투영체로부터 유도되는, **셀 자연 궤도함수(Cell Natural Orbitals, CNOs)**라고 불리는 체계적인 국소 궤도함수 구축법을 소개한다. 핵심적인 수학적 대상은 **단위 셀 일입자 축약 밀도 행렬(Unit-Cell One-Particle Reduced Density Matrix, UC-1pRDM)**이며, 이를 $L$로 표기한다.

1. UC-1pRDM의 구축: $P$를 대상 블로히 밴드에 대한 투영체라 하고, $\Pi$를 힐베르트 공간을 하나의 참조 단위 셀로 제한하는 투영체라고 하자. UC-1pRDM은 $L = \Pi P \Pi$로 정의된다. 궤도함수 기저에서 행렬 요소는 $L_{\alpha\beta} = \langle c^\dagger_\beta c_\alpha \rangle$로 주어지며, 이는 대상 밴드들을 합산하고 단위 셀 내 궤도함수의 실공간 위치를 고려하는 임베딩 행 matrix $V_k$를 포함하여 계산된다.
2. 대각화: $L$의 고유벡터 $|\tau_a\rangle$는 CNOs가 된다. 그에 대응하는 고유값 $\lambda_a$는 대상 밴드 내에서의 각 CNO의 점유율을 나타낸다.
3. 물리적 해석:
   • 스펙트럼 가중치(Spectral Weight): $\lambda_a$는 각 국소 궤도함수의 중요도를 정량화한다. $\lambda_1 \approx 1$인 단일 궤도는 해당 밴드가 단일 궤도에 의해 잘 설명됨을 의미하며, 분수 점유율은 다중 궤도 특성을 나타낸다.
   • 엔탱글먼트(Entanglement): 고유값 $\lambda_a$는 국소 궤도함수와 밴드 부공간(subspace) 사이의 사상(map)에 대한 특이값(singular values)이다. 이들은 단위 셀 분할에 따른 엔탱글먼트 스펙트럼을 정의하며, 본 네이만 엔트로피 $S_{cell}$은 점유된 부공간이 셀 경계를 가로질러 인자화(factorize)되는 데 실패하는 정도를 측정한다.
   • 대칭성 및 위상: 특정 와이코프 위치(Wyckoff position) $q$를 중심으로 단위 셀 분할을 선택함으로써, CNO는 사이트 대칭군 $G_q$의 기약 표현(irreps)으로 변환된다. 이를 통해 CNO는 후보 원소 밴드 표현(Elementary Band Representations, EBRs)을 식별할 수 있게 한다.
   • 영점(Zeros): 포락선 함수(envelope functions) $s_{a,k}$(재조정된 블로흐 계수)는 브릴루앙 영역에서의 영점을 드러낸다. 위상적 장애(예: 0이 아닌 Chern number)는 제거할 수 없는 영점을 강제하며, 폐쇄된 원자 한계 또는 호환되지 않는 와이코프 선택은 대칭성에 의해 강제된 영점을 유발하며, 이는 분할 중심을 이동함으로써 제거될 수 있다.

주요 기여 및 결과

• 체계적인 기저 선택: 저자들은 궤도함수의 수나 특정 원자적 성격에 대한 사전 가정 없이 최적의 국소 궤도함수 기저를 결정하는 방법을 제공한다. CNO는 블로흐 파동함수와 실공간 분할의 선택에 의해 유일하게 결정된다.
• 폐쇄된 원자 한계(Obstructed Atomic Insulators)에 대한 적용: 저자들은 1H-TMD 모델을 사용하여, 원자 사이트(1a 와이코프 위치)를 중심으로 단위 셀 분할을 선택하면 분수 점유율과 포락선 함수의 대칭성 강제 영점이 발생하여 밴드의 폐쇄적 특성을 반영함을 보여준다. 분할 중심을 와니에 중심(1c 와이코프 위치)으로 옮기면 $\lambda_1 \approx 1$인 지배적인 CNO와 영점이 없는 상태가 되어, 최소한의 국소적 설명을 올바르게 식별한다.
• 위상 밴드(Topological Bands)에 대한 적용: 스핀 없는 BHZ 모델에서, CNO 점유율은 위상 전이 과정을 따라 매끄럽게 변하지만, 포락선 함수는 위상 단계($C=1$)에서 제거할 수 없는 영점을 보이며, 이는 실공간에서의 와니에 장애(Wannier obstruction)를 나타낸다.
• Twisted Bilayer WSe$_2$의 국소 모델 구축: 저자들은 다양한 비틀림 각도($1^\circ$ ~ $5^\circ$)에 대해 비틀린 이중층 WSe$_2$ (tWSe$_2$)에 이 방법을 적용한다.
  • 이들은 상위 두 모아레 밴드가 비틀림 각도와 분할의 선택(MM, MX, 또는 XM 적층 사이트)에 따라 서로 다른 국소적 설명이 필요함을 확인하였다.
  • 위상 영역($\theta > 1.51^\circ$)에서는 유한한 총 Chern number로 인해 단일 국소 궤도가 밴드를 포괄할 수 없다.
  • 총 스펙트럼 가중치를 최대화하고 완비성 조건(중첩 행렬 $S_k$의 비제로 행렬식)을 만족하는 CNO 부분 집합(예: MM, MX, XM 궤도를 결합한 다중 중심 집합)을 선택함으로써, 이들은 타이트 바인딩 모델을 구축한다.
  • 단일 샷 와니에화(Single-Shot Wannierization): 이러한 CNO 부분 집합은 단일 샷 와니에화 절차를 위한 시험 상태(trial states) 역할을 한다. 결과적으로 얻어진 타이트 바인딩 모델은 밴드 분산뿐만 아니라 연속체 모델의 양자 기하학(Berry curvature)까지 재현한다.
  • 트레이드오프(Trade-offs): 저자들은 MM 중심 모델(한 사이트 위의 3개 궤도)과 다중 중심 모델(서로 다른 사이트 위의 궤도)을 비교한다. 다중 중심 모델은 더 국소화된 와니에 함수(더 짧은 범위의 호핑)를 생성하지만, 대칭성을 보존하는 다체 처리를 위해 더 큰 클러스터가 필요하다. 반면, MM 중심 모델은 $C_3$ 대칭성을 더 쉽게 보존하지만 더 긴 범위의 상호작용을 가진다.

의의 및 주장
본 논문은 CNO가 대상 밴드의 다중 궤도 특성에 대한 직접적이고 파동함수 유도적인 척도를 제공한다고 주장한다. 화학적 성질이나 대칭성 고유값으로부터 궤도를 추론하는 전통적인 방식과 반대로, CNO를 통해 국소적 자유도를 밴드 매니폴드로부터 직접 추론할 수 있다.

이 연구의 의의는 다음과 같다:

1. 최적성: CNO는 양자 기하학을 보존하면서 대상 밴드를 포괄하는 데 필요한 궤도함수의 수를 최소화하는 최적의 시험 상태를 제공한다.
2. 진단 능력: 점유 스펙트럼과 포락선 함수의 영점은 위상적 장애와 폐쇄된 원자 한계를 구분하는 진단 도구로서, 제거 가능한(대칭성) 영점과 제거할 수 없는(위상적) 영점을 구별한다.
3. 상호작용 모델의 기초: CNO는 단위 셀 상관 행렬을 대각화하므로, 밴드 투영된 상호작용의 자연스러운 기저를 제공한다. 점유율 계층 구조 $\lambda_a$는 상호작용 척도를 정렬하며, 임의의 선택 없이 다중 궤도 허바드 모델을 구축하고 양자 임베딩 방법(예: DMFT)을 위한 임퓨리티 공간을 정의하는 경로를 제시한다.
4. 일반성: 이 방법은 격자 모델과 연속체 모델(평면파 기저로부터 유도된 실공간 그리드 사용) 모두에 적용 가능하며, 밴드 투영체만을 입력값으로 요구하므로 제일원리 계산과 호환된다.

저자들은 결론적으로, 이러한 관점이 제일원리 파동함수 및 연속체 모델과 상관된 위상 양자 물질의 유효 이론에 필요한 최소한의 국소적 자유도 사이를 연결한다고 밝힌다.