전이금속 칼코게나이드의 라시바 스핀분할: 층별 궤도 비대칭과 내재 전계의 기원

초록

**
본 연구는 2차원 전이금속 칼코게나이드(TMD)에서 역전 및 미러 대칭이 동시에 깨질 때 나타나는 라시바 스핀분할의 미시적 기원을 첫 원리 계산으로 규명한다. 단층과 이층 구조의 전계 의존 라시바 파라미터 λₙᴿ와 내재 궤도 전계 Eₙ⁰를 추출하고, 궤도 편극 불균형 ξₙ을 새로운 주문 매개변수로 도입해 스핀 순서를 설명한다. 결과는 이층이 단층보다 라시바 분할이 작으며, 이는 내부 전기편극과 층간 혼성의 경쟁에 기인함을 보여준다.

**

상세 분석

**
본 논문은 2D TMD 소재에서 라시바 스핀‑궤도 결합을 미시적으로 이해하기 위해 DFT‑GGA(PBE)+D3‑BJ 기반 전산을 수행하고, 전자 구조를 44개의 Wannier 함수(전이금속 d와 칼코게나이드 p 궤도)로 전환하였다. 단층(Monolayer, ML)에서는 전기장 E_ext가 없을 경우 미러 대칭 m_z가 보존되어 선형 라시바 항이 금지되고, 외부 전계가 가해질 때만 α_R∝E_ext가 나타난다. 반면 R‑형 이층(H‑stack이 아닌 R‑stack)에서는 구조적으로 I와 m_z가 모두 깨져, 외부 전계 없이도 내재 전계 E_n⁰가 존재한다. 라시바 해밀토니안 H_R=α_nᴿ(k_xσ_y−k_yσ_x)에서 α_nᴿ=λ_nᴿ(E_ext+E_n⁰)로 정의하고, k와 E_ext에 대한 선형 구간에서 ∆E=α_nᴿk를 직접 피팅함으로써 λ_nᴿ와 E_n⁰를 추출하였다.

주요 결과는 다음과 같다. (1) λ_nᴿ는 밴드마다 크게 다르며, B‑밴드(결합형)에서 양의 값이, A‑밴드(반결합형)에서는 음의 값이 관측된다. (2) 전계 없이도 B‑밴드와 A‑밴드가 서로 반대 방향의 스핀 텍스처를 보이며, E_ext=−E_n⁰일 때 스핀 분할이 소멸한다. (3) 무거운 칼코게나이드(Se, Te)로 갈수록 λ_nᴿ와 |E_n⁰|가 증가하는데, 이는 원자 전기적 편극성(이온성 감소)과 SOC 강도가 동시에 작용하기 때문이다. 반면 전이금속(Mo, W)의 무게 변화는 λ_nᴿ에 미미한 영향을 미치며, 이는 SOC와 궤도 편극성 사이의 상쇄 효과로 해석된다.

핵심적인 새로운 개념으로 도입된 궤도 편극 불균형 ξ_n은 각 밴드의 전하 분포를 z축에 대해 비대칭적으로 정량화한다. 단층에서는 ξ_n^ML = (∫{z_M}^{∞}ρ−∫{−∞}^{z_M}ρ)/C 로 정의하고, 이층에서는 두 단층의 ξ_n^L1, ξ_n^L2를 평균해 ξ_n을 얻는다. ξ_n은 E_ext에 따라 급격히 변하는 스텝을 보이며, 이 스텝이 바로 ∆E가 0이 되는 전계와 일치한다. 따라서 ξ_n은 라시바 스핀 분할을 실험적으로 추적할 수 있는 물리량으로서, 스핀 순서(양/음)와 내재 전계의 부호를 동시에 결정한다.

또한, 퍼트된 Wannier 기반 모델을 통해 d_{xz/yz}와 p_{x/y} 궤도가 1차 라시바 과정에 핵심 매개체임을 확인하였다. 이들 궤도는 층간 하이브리다이제이션에 민감하게 반응하며, 내부 전기 편극과 결합해 라시바 파라미터를 조절한다. 결과적으로, 이층에서 관측되는 라시바 분할이 단층보다 작아지는 현상은 (i) 내부 전계가 이미 존재해 외부 전계에 대한 민감도가 감소하고, (ii) 층간 혼성으로 인해 반결합 A‑밴드와 결합 B‑밴드의 궤도 구성이 서로 상쇄되기 때문임을 설명한다.

이러한 정량적·정성적 분석은 라시바 스핀트로닉스 구현을 위한 소재 설계에 직접적인 가이드를 제공한다. 예를 들어, 큰 λ_nᴿ와 |E_n⁰|를 원한다면 무거운 칼코게나이드(Se, Te)와 R‑형 이층 구조를 선택하고, 외부 전계는 E_n⁰와 반대 부호로 조절해 스핀 텍스처를 역전시킬 수 있다. 또한, ξ_n을 실험적으로 측정(예: 전하 분포를 탐지하는 STEM‑EELS)함으로써 라시바 효과를 비파괴적으로 검증할 수 있다.

**