Evaluation of real-space second Chern number using the kernel polynomial method

기술 요약: 커널 다항식 방법을 이용한 실공간 제2 천 르 수(Second Chern Number)의 평가

문제 정의
4차원(4D) 천 절연체에서 나타나는 제2 천 르 수($C_2$)나 6차원(6D) 시스템에서 나타나는 제3 천 르 수($C_3$)와 같은 $2n$ 차원 시스템의 위상적 상들은 전통적으로 운동량 공간(momentum space)에서 정의된다. 이러한 정의는 번역 불변성(translational invariance)에 의존하는데, 이는 무질서(disorder)나 준주기성(quasiperiodicity)이 포함된 실제적인 시스템에서는 흔히 결여되는 특성이다. 이를 해결하기 위해 실공간 공식들이 존재하지만, 직접적인 행렬 대각화에 의존하는 표준적인 접근 방식은 약 $10^4$개의 자유도를 가진 시스템 크기로 제한된다. 이러한 크기 제약은 차원이 높아질수록 힐베르트 공간의 차원이 지수적으로 증가($L^{2n}$)하기 때문에 특히 높은 차원의 시스템에서 문제가 되며, 위상 불변량의 양자화가 명확하게 나타나기 위해 필요한 충분히 큰 시스템 크기를 확보하는 것을 어렵게 만든다. 슈퍼셀 근사법이나 키타에프(Kitaev)의 공식을 사용한 기존의 실공간 $C_2$ 계산 시도는 $4^4$ 또는 $8^4$의 시스템 크기에 국한되었으며, 유한 크기 효과(finite-size effects)로 인해 완전한 양자화를 달eric하지 못하는 경우가 많았다.

방법론
저자들은 훨씬 더 큰 시스템 크기에 대해 실공간 천 르 수를 평가하기 위해 **커널 다항식 방법(Kernel Polynomial Method, KPM)**을 채택하였다. 이들의 접근 방식의 핵심은 다음과 같다:

1. 모델: 저자들은 격자 상수 $a=1$인 초입방 격자(hypercubic lattice) 상에 이산화된 $2n$ 차원 윌슨-디락(Wilson-Dirac) 격자 모델을 활용한다. 해밀토니안은 클리포드 대수(Clifford algebra)를 만족하는 감마 행렬들을 사용하여 정의된다.
2. 실공간 공식: 저자들은 $n$차 천 르 수를 위한 일반적인 보편적 위상 마커(universal topological marker) 공식을 채택한다:
   $$C_n = \frac{(2\pi i)^n}{n!} \epsilon_{j_1 j_2 \dots j_{2n}} \text{Tr} (P X_{j_1} P X_{j_2} \dots P X_{j_{2n}} P)$$
   여기서 $P$는 점유된 밴드(occupied bands)로의 투영 연산자이며, $X_j$는 좌표 연산자이다.
3. KPM 구현: 비용이 많이 드는 대각화를 피하기 위해, 계단 함수(step function)를 $M=256$ 차수의 체비쇼프 다항식(Chebyshev polynomials)을 사용하여 확장함으로써 $P$를 근사한다. 트레이스(trace)는 $R=5$개의 무작위 위상 벡터를 사용하는 확률적 트레이스 근사(stochastic trace approximation)를 통해 추정된다. 수치적 안정성을 보장하기 위해 해밀토니안은 $[-1, 1]$ 구간으로 재스케일링된다.
4. 무질서: 안데르슨형 무질서(Anderson-type disorder)는 $[-W, W]$ 범위 내에서 균등 분포된 무작위 온사이트 에너지(on-site energies)를 통해 도입된다.

주요 기여 및 결과

• 깨끗한 시스템에서의 4D 제2 천 르 수 ($C_2$):
  저자들은 최대 $30^4$ 사이트를 가진 4D 시스템에 대해 실공간 계산을 수행하였다 (본문에는 304개 사이트로 표기되어 있으나, $30^4$ 표기는 한 변의 길이 $L=30$을 의미함).

  • 양자화: 시스템 크기($L$)가 4에서 20으로 증가함에 따라, 수치적으로 계산된 $C_2$는 이론적인 운동량 공간 예측값($C_2 = -3$ for 특정 질량 매개변수)으로 수렴하였다.
  • 수렴 속도: $\ln(|C_2 + 3|)$와 시스템 크기 $L$ 사이의 편차 분석은 지수적 감소를 보여주었으며, 이는 유한 크기 효과가 더 큰 시스템에서는 무시할 수 있는 수준이 됨을 확인시켜 주었다.

• 무질서한 시스템에서의 4D 제2 천 르 수 ($C_2$):

  • 견고성(Robustness): 이 방법은 약한 무질서 하에서 위상적 상의 안정성을 성공적으로 포착하였다. 양자화된 $C_2$ 값은 적당한 무질서 강도($W$)에서도 안정적으로 유지되었다.
  • 상도표(Phase Diagram): 무질서 평균화된 $C_2$ 상도표는 자기 일관적 본 근사(self-consistent Born approximation)로부터의 예측과 잘 일치하였다.
  • 임계 거동: 임계점 근처(벌크 갭이 작은 지점)에서는 아주 약한 무질서조차 위상 마커를 양자화된 값에서 벗어나게 했으나, 갭이 더 큰 시스템은 더 강한 무질서에 대해서도 견고함을 유지했다.

• 6D 제3 천 르 수 ($C_3$) 탐색적 계산:
  저자들은 실공간 제3 천 르 수를 계산하기 위해 이 방법을 6D 시스템으로 확장하였다.

  • 질적 일치: $L=5$까지의 시스템 크기에 대해, 수치적 결과는 이론적인 운동량 공간 예측과 질적으로 일치하였으며, 상전이 위치를 정확히 식별하였다.
  • 한계: 완전한 양자화는 달성되지 않았다. 저자들은 이를 상당한 유한 크기 효과와 극심한 계산 복잡성 때문이라고 설명한다. $C_3$ 계산은 4D의 경우 24개의 항이 필요한 것에 비해 720개의 6D 레비-치비타 텐서 항을 합산해야 하므로 자원 요구량이 급격히 증가한다.

의의 및 주장
본 논문은 고차원 위상 상의 실공간 특성화에 있어 "진일보한 성과"를 나타낸다고 주장한다. 주요 의의는 커널 다항식 방법이 이전의 실공간 접근 방식($8^4$)보다 수십 배 더 큰 시스템 크기($30^4$)에서 실공간 제2 천 르 수를 계산할 수 있음을 입증했다는 데 있다. 이 규모는 4D 시스템에서 양자화된 값으로의 지수적 수렴을 관찰하기에 충분하며, 무질서한 위상 절연체를 분석하는 데 있어 이 방법의 유용성을 검증한다.

6D 확장과 관련하여, 저자들은 현재의 계산적 한계로 인해 완전한 양자화가 불가능함을 인정하며 신중한 어조를 유지하고 있다. 그들은 6D에서의 양자화된 $C_3$ 값을 얻기 위해 필요한 더 큰 시스템 크기에 접근하려면 향후 텐서 네트워크 기술과의 통합이 필요할 수 있다고 제안한다. 이 연구는 새로운 실험 설업을 제안하는 것이 아니라, 무질서가 존재하는 기존의 이론적 모델을 분석하기 위한 강력한 수치적 도구를 제공한다.