Ti/Zr/Hf-Bi/Sb 허니컴 단일층에서 양자 비정상 홀 효과: 제1원리 예측

초록

2차원 소재 연구의 발전은 새로운 위상학적 상의 기회를 지속적으로 드러내고 있습니다. 양자 비정상 홀 효과(QAHE)는 최근 자성 도핑 위상 절연체에서 입증되었으나 극저온에 국한되어 있습니다. 밀도범함수 이론을 사용해 전이금속 Ti, Zr, Hf와 전형소 Bi, Sb로 구성된 단일층 허니컴(MPn, M=Ti, Zr, Hf; Pn=Bi, Sb)이 고유한 QAHE를 가질 것을 예측했습니다. 연구 시스템 중 평면 TiBi와 HfBi가 각각 15 meV와 7 meV 밴드갭을 가진 QAH 절연체로 나타났으며, 이는 적당한 인장 변형 하에서입니다. 격자 상수와 버클링 거리 변화가 위상학적 상 전이를 유도하는 방식, 스핀-궤도 결합에 의한 밴드 반전 관찰, Wannier 유도 타이트바인딩 모델로 지그재그 나노리본에서 카이랄 모서리 모드 확인 등을 통해 실현 가능한 QAH 물질 목록을 확장하고 저에너지 스핀트로닉 소자 플랫폼을 제안합니다.

배경

2차원 위상 절연체(2D TI)는 절연 벌크를 유지하면서 스핀 분극된 무간극 모서리 상태를 보이며, 이는 양자 스핀 홀(QSH) 효과의 특징입니다[1]. 시간 반전 대칭(TRS)이 깨지면(고유 자성이나 외부 섭동에 의해) 나선형 모서리 채널이 카이랄 채널로 변환되어 양자 비정상 홀 효과(QAHE)가 발생합니다[3]. QAHE의 소산 없는 자기장 없는 전하 수송은 초저전력 스핀트로닉스에 매력적입니다[4], 그러나 실험적 구현은 희석 자성 위상 절연체에 국한되어 왔습니다[8]. 이론 연구는 강자성과 강한 스핀-궤도 결합(SOC)을 QSH 물질에 도입하면 위상학적 밴드갭을 열고 0이 아닌 체른 수를 생성할 수 있음을 보여줍니다[9,10]. IV족(Sn)과 V족(Bi, Sb) 단일층 등 많은 QSH 시스템이 광범위하게 연구되었으나, 전이금속 기반 허니컴의 QAHE 잠재력은 미개척 상태입니다. 본 연구는 Bi/Sb 허니컴에 Ti, Zr, Hf 치환이 고유 QAH 상을 안정화할 수 있는지 제1원리 계산으로 탐구함으로써 이 공백을 메웁니다. 최근 III-V 허니컴이 QSH 상을 가질 수 있고 전이금속 도핑이 자성을 유도할 수 있다는 예측[17,25-27]을 기반으로 합니다.

결과 및 논의

원시 허니컴 격자에서 Bi 원자 절반을 전이금속(Ti, Zr, Hf)으로 치환하면 MPn 단일층이 생성되며, 이는 낮은 버클링 또는 평면 구조를 취할 수 있습니다. 그림 1a는 M-Pn 격자의 평면도, 그림 1b, c는 버클링 및 평면 구조의 측면도, 그림 1d는 고대칭점을 가진 제1 브릴루앙 존을 보여줍니다.

a M-Sb/Bi 허니컴 결정 구조. b, c 버클링 및 평면 구조 측면도. d 고대칭점을 가진 제1 브릴루앙 존(BZ)

격자 상수를 변화시키고 원자 위치를 완화해 버클링 및 평면 상의 구조적 안정성과 변형 효과를 평가했습니다. TiBi, ZrBi, HfBi의 총에너지 곡선을 그림 2a-c에 도시했습니다. 낮은 버클링 MPn 구조가 평면 구조보다 일관되게 낮은 에너지를 가집니다. 그럼에도 QAHE는 주로 변형된 평면 영역에서 나타나며, 버클링 HfBi에서만 좁은 격자 상수 창이 QAH 상을 지지합니다(Fig. 2c 참조).

a TiBi, b ZrBi, c HfBi의 위상 다이어그램: 서로 다른 격자 상수에서 총 에너지. 다이어그램은 QAH(양자 비정상 홀 상), I(절연체), SM(반금속) 영역으로 나뉨. 파란 원과 빨간 삼각형은 각각 버클링과 평면 경우를 나타냄

표 1은 모든 MPn 시스템의 평형 격자 상수, 자기 모멘트, 밴드갭, 계산된 체른 수를 나열합니다. 0이 아닌 정수 체른 수는 QAH 절연체를 식별합니다: 평면 TiBi(C=+1, Δ=15 meV)와 평면 HfBi(C=-1, Δ=7 meV)가 가장 유망한 후보입니다. TiBi의 상 전이는 버클링 거리 감소로(Fig. 3), HfBi는 인장 변형으로 유도됩니다(Fig. 4).

버클링 거리 변화에 의한 상 전이. a a=4.6 Å에서 TiBi 위상 다이어그램. 화살표는 전이 경로 표시. b-f 버클링 거리(δ)를 0.44에서 0.4 Å로 줄일 때 밴드 구조 전이. 전이는 δ=0.41 Å에서 발생

격자 상수 변화에 의한 상 전이. a 버클링 HfBi 위상 다이어그램. 화살표는 전이 경로. b-h 격자 상수를 4.7에서 5.1 Å로 증가시킬 때 밴드 구조 전이

평형 기하학에서의 전자 밴드 구조는 그림 5에 제시되어 있습니다. 평면 기하학에서 TiBi는 명확한 밴드 반전과 15 meV 갭을 보이며, 평면 HfBi는 C=-1로 7 meV 갭을 보입니다. 그러나 버클링 HfBi는 높은 체른 수(C=-3)를 가진 반금속입니다. 버클링 ZrBi, TiSb 등 다른 구성은 자명한 절연체로 남습니다.

평형 격자 상수에서 M-Pn(M=Ti, Zr, Hf; Pn=Sb, Bi) 전자 밴드 구조: a 평면, b 버클링 경우. 빨간/파란 원은 +s_z, -s_z 기여 표시

SOC와 자성의 역할을 분석하기 위해 평면 TiBi(a=4.76 Å)에 초점을 맞춥니다. 그림 6은 비자성 계산이 금속 상태를 보이나, 강자성 정렬(단위 세포당 자기 모멘트 ≈1.05 μ_B, Ti 주도)이 반금속 스펙트럼을 생성함을 보여줍니다. SOC 포함 시 15 meV 갭이 열리며, 밴드 반전이 SOC 주도이며 QAHE를 담당함을 확인합니다.

평면 TiBi 박막 a=4.76 Å에서 전자 밴드 구조: 비자성 계산 (a) SOC 없이, (c) SOC 포함; 강자성 계산 (b) SOC 없이, (d) SOC 포함. 빨간/파란 원은 (c) 비자성, (d) 강자성 계산에서 +s_z, -s_z 기여 표시

모서리 상태는 위상학적 상의 필수적 지문입니다. Wannier 유도 타이트바인딩 해밀토니안을 사용해 127 Å 폭 지그재그 HfBi 나노리본을 구성했습니다(Fig. 7). 밴드 구조는 페르미 레벨을 홀수 번(|C|=1) 가로지르는 카이랄 모서리 모드를 보여주며 평면 HfBi의 QAH 특성을 확인합니다.

버클링 HfBi 지그재그 나노리본 모서리를 따른 밴드 구조(a=4.9 Å, 폭 127 Å). 파란(빨간) 원은 왼쪽(오른쪽) 모서리 기여 표시. 주황색 채워진 영역은 벌크 밴드 표시

포논 계산에서 자유 단일층의 작은 허수 주파수가 나타나 실험적 구현을 위해 기판이 필요할 수 있음을 시사합니다. 또한 2×2 초세포 내에서 강자성 배열이 버클링 경우 반강자성보다 에너지적으로 유리하며, 평면 시스템은 FM과 AFM 상태가 근축퇴됨을 보입니다.

결론

제1원리 연구는 Ti, Zr, Hf를 Bi 또는 Sb 허니컴 격자에 치환하면 고유 QAHE를 유도할 수 있음을 입증했습니다. 버클링 기하학이 에너지적으로 선호되나, 적당한 인장 변형이나 버클링 거리 감소가 0이 아닌 체른 수를 가진 평면 구성으로의 전이를 유도합니다. 평면 TiBi와 HfBi는 각각 15 meV와 7 meV 밴드갭을 가진 견고한 QAH 절연체로 나타납니다. 이 물질들은 성장하는 2차원 QAH 후보 목록에 추가되며 저전력 스핀트로닉 소자를 위한 유망한 플랫폼입니다.

실험 방법

모든 계산은 VASP 5.3에 구현된 일반화 기울기 근사(GGA)와 투영자 증강파(PAW) 포텐셜을 사용한 밀도범함수 이론 프레임워크 내에서 수행되었습니다. 평면파 컷오프는 350 eV로 설정했고, 잔여 힘이 5×10⁻³ eV/Å 미만이 될 때까지 구조를 완화했습니다. 전자적 자기일관성 루프 수렴은 10⁻⁶ eV에서 달성했습니다. z축 20 Å 진공 간격으로 층간 상호작용을 방지했고, 24×24×1 감마 중심 Monkhorst-Pack 격자로 2D 브릴루앙 존을 샘플링했습니다. WANNIER90으로 최대 국소화 Wannier 함수를 생성해 타이트바인딩 모델 구성과 모서리 상태 계산을 가능케 했습니다. 위상학적 불변량은 Z2Pack 패키지로 하이브리드 Wannier 전하 중심을 추적해 체른 수를 얻었습니다.