Higher Berry curvature, second Chern numbers and magnetoelectric coupling in crystalline insulators

기술 요약: 결정성 절연체에서의 높은 차수의 베리 곡률, 제2 천 수, 그리고 자기전기 결합

문제 정의
위상 불변량은 양자 다체 바닥 상태를 분류하는 데 필수적이다. 베리 곡률(2-form)과 그와 관련된 천 수(Chern number)가 브릴루앙 존의 0차원 바닥 상태(예: 블로흐 상태)를 성공적으로 특징짓는 반면, 최근의 이론적 발전은 공간적으로 확장된 시스템으로의 일반화를 제안한다. 이 프레임워크에서, $d$ 차원 공간의 단거리 얽힘 상태(short-range entangled states)의 가계(families)는 더 높은 차수의 미분 형식(예: 1차원 시스템을 위한 3-form)에 의해 특징지어진다. 구체적으로, 3차원 공간으로 매개되는 1차원 시스템의 경우, 관련 불변량은 더 높은 차수의 베리 곡률로부터 유도된 딕스미어-두디-카푸스틴-스포이네이코(Dixmier–Douady–Kapustin–Spodyneiko, DDKS) 수 $\nu$이다.

일반적인 상호작용 시스템 또는 자유 페르미온 시스템에 대해 이러한 고차 불변량을 명시적으로 양자화하여 계산하는 방법에는 상당한 간극이 존재한다. 4차원 디락 모델과 같은 특정 자유 페르미온 모델에 대해서는 해석적인 표현식이 존재하지만, 4차원에서 제2 천 수($ch_2$)를 계산하기 위한 일반적인 수치적 방법들은 격자 게이지 이론이나 실공간 자기장 응답에 의존한다. 이러한 표준 방법들은 "명시적 양자화(manifest quantization)"가 결여되어 있는 경우가 많은데, 즉 연속체 극한에서만 정수 값으로 수렴하거나 광범위한 메쉬 세분화(mesh refinement)를 요구하며, 구조 자체에 의해 정수를 반환하지는 못한다. 또한, 3차원 공간에서의 DDKS 수의 물리적 해석, 특히 궤도 자기전기 결합 및 Chern-Simons 액시온 각도($\theta_{CS}$)와의 관계는 여전히 미해결 과제로 남아 있다.

방법론
저자들은 4차원 Chern 절연체의 격자 모델을 3차원 브릴루앙 존($T^3$)에 의해 매개되는 변환 불변 무한 1차원(1D) 사슬의 가계로 재구성함으로써 이 문제들을 해결한다.

1. 모델 재구성: 4차원 디락 해밀토니안으로부터 시작하여, 저자들은 하나의 운동량 성분($k_4$)에 대해 부분 푸리에 변환을 수행하여 매개변수 $\lambda = (\lambda_1, \lambda_2, \lambda_3)$에 의존하는 1차원 해밀토니안 $\{H_\lambda\}$의 가계를 얻는다. 이 해밀토니안들은 단위 격자당 두 개의 궤도를 갖는 1차원 사슬을 기술한다.
2. 수치적 프레임워크 (iMPS): 1차원 사슬의 바닥 상태는 무한 행렬 곱 상태(infinite Matrix Product States, iMPS)로 근사된다. 저자들은 더 높은 차수의 베리 곡률을 계산하기 위해 Ref. [13]에서 제안된 알고리즘을 일반화하여 사용한다.
3. 이산적 고차 베리 곡률:
   • 파라미터 공간 $T^3$는 그리드로 이산화되고 사면체(tetrahedra)로 삼각 분할된다.
   • 각 정점(vertex)에 대해, iMPS 바닥 상태가 iDMRG를 통해 구해진다.
   • 인접한 iMPS 텐서 사이의 전이 함수(transition functions)는 혼합 전이 행렬(mixed transfer matrices)의 지배적인 고유벡터를 사용하여 계산된다.
   • 삼각형 면(faces)에 대해, 결합 차원(bond dimension)의 변화를 처리하기 위해 슈미트 스펙트럼(Schmidt spectrum)에 가중치를 둔 게이지 불변 위상 $\phi(\Delta_2)$가 정의된다.
   • 사면체에 대한 고차 베리 곡률 $F(\Delta_3)$는 이들의 경계에 대한 방향성이 있는 위상의 합으로 정의된다.
   • DDKS 수 $\nu$는 $T^3$의 전체 삼각 분할에 대해 $F(\Delta_3)$를 합산하여 얻어지며, 이는 명시적으로 양자화된 정수를 산출한다.
4. 비교 분석: 저자들은 동일한 모델에 대해 표준 격자 게이지 이론 접근법(Ref. [34])을 사용하여 제2 천 수 $ch_2$를 계산하고, 그 결과 나타나는 위상 도표를 비교한다.
5. 자기전기 결합: 연구는 DDKS 수와 3차원 궤도 자기전기 결합을 지배하는 Chern-Simons 액시온 각도 $\theta_{CS}$ 사이의 관계를 조사한다. 저자들은 펌핑 사이클(pumping cycle)을 따라 제2 천 형식(second Chern form)으로부터 유도된 $\theta_{CS}$의 변화와 "고차 베리 위상(higher Berry phase)" $\theta_{HBP}$(iMPS 구축으로부터 유도됨)의 변화를 비교한다.

주요 기여 및 결과

• $ch_2$의 명시적 양자화 계산: 저자들은 iMPS 기반의 고차 베리 곡률 알고리즘을 통해 계산된 DDKS 수가 4D Chern 절연체의 제2 천 수 $ch_2$와 정확히 일치하는 위상 도표를 생성함을 입증한다. 이 일치는 질량 영역($m$)에 관계없이 유지되며, 표준 격자 게이지 이론 계산이 심각한 유한 그리드 오차를 겪는 영역에서도 마찬가지이다. iMPS 방법은 거친 그리드(예: $4^3$)에서도 정확한 정수를 산출하는 반면, 격자 게이지 방법은 수렴을 위해 훨씬 더 미세한 그리드($10^4$ 포인트)를 요구한다.
• Künneth 공식을 통한 이론적 연결: 논문은 Künneth 공식을 사용하여 4D 제2 천 수와 3D DDKS 수 사이의 엄격한 연결을 확립한다. $T^4$를 $S^1 \times T^3$로 분해함으로써, 저자들은 $T^4$ 상의 제2 Chern class가 $H^1(S^1, \mathbb{Z})$의 생성자와 $\nu$에 의해 측정되는 $H^3(T^3, \mathbb{Z})$의 클래스의 곱에 대응함을 보여준다. 따라서 이 특정 분해에 대해 $ch_2 = \nu$이다.
• 4D 디락 모델에서의 자기전기 결합: 4D 디락 모델에 대해, 저자들은 펌핑 사이클을 따르는 Chern-Simons 액시온 각도의 변화 $\Delta\theta_{CS}$가 고차 베리 위상 $\Delta\theta_{HBP}$의 변화에 의해 정확하게 재현됨을 발견했다. 이는 첫 번째 Chern 수가 0인 슬라이스들의 경우, 고차 베리 위상이 자기전기 응답과 직접적인 대응 관계를 가짐을 시사한다.
• Hopf 절연체에서의 반례: 저자들은 비자명한 Hopf 불변량을 갖는 2-밴드 모델인 Hopf 절연체에서 이 대응 관계를 테스트한다. 이 경우, Chern-Simons 각도는 Hopf 불변수에 의해 양자화되지만, 계산된 고차 베리 위상 $\theta_{HBP}$는 양자화되지 않으며 $\theta_{CS}$와도 일치하지 않는다. 저자들은 이러한 불일치를 Hopf 불변의 "섬세한(delicate)" 성질 때문이라고 설명하는데, 즉 Hopf 불변은 운동량 공간에서의 2-밴드 모델에 대해서만 엄격하게 잘 정의되며, 이 조건이 열역학적 극한의 실공간 iMPS 기술에서는 유지되지 않을 수 있다는 것이다.

의의 및 주장
본 논문은 고차 베리 곡률과 iMPS를 사용하여 4D 격자 모델의 제2 천 수를 명시적으로 양자화하여 계산한 첫 사례라고 주장한다. 이는 고차 베리 곡률이 공간적으로 확장된 상태들의 가계에 대한 고차 위상 불변량을 탐지하는 도구가 될 수 있다는 이론적 프레임워크를 검증한다.

저자들은 이 결과를 양자 기하학 및 위상적 성질의 맥락에 배치하며, 고차 베리 위상과 자기전기 결합 사이의 대응 관계가 4D 디락 모델에서는 성립하지만, Hopf 절연체의 사례에서 나타나듯 보편적이지는 않다는 점에 주목한다. 저자들은 이 방법이 단일 입자 블로흐 기술이 실패하는 상호작용 시스템에 특히 적합하다고 강조하지만, 현재는 자유 페르미온 벤치마크에 집중하고 있다. 이 연구는 합성 차원(synthetic dimensions)과 준주기적 구동(quasi-periodic driving)이 고차원 위상 현상을 구현하고 탐구하는 실험적 플랫폼에서 활용될 수 있음을 강조한다.

논문은 고차 베리 곡률과 3차원에서의 물리적 응답 사이의 관계가 여전히 열린 문제, 특히 $\theta_{CS}$와의 보편적 연결성에 관한 문제가 남아 있다고 언급하며 겸허하게 결론을 맺는다. 이는 고차 베리 위상이 표준적인 Chern-Simons 응답과는 구별되는 기하학적 정보를 포착할 수 있음을 시사하며, 비선형 수송 및 바닥 상태 물리학에서의 역할을 규명하기 위한 추가 연구의 필요성을 제시한다.