First-principles study of bulk stacking, $J_{\rm eff}$ picture, magnetic Hamiltonian, $g$ factors, and structural distortions of $α$-RuCl$_3$

기술 요약: $\alpha$-RuCl$_3$의 벌크 적층, $J_{\text{eff}}$ 모델, 자기 해밀토니안, $g$ 인자 및 구조적 왜곡에 관한 제일원리 연구

문제 제기
$\alpha$-RuCl$_3$는 키타에프(Kitaev) 양자 스핀 액정 상태를 구현하기 위한 주요 후보 물질이지만, 그 근본적인 자기적 특성은 여전히 논쟁의 대상입니다. 주요 미해결 과제는 다음과 같습니다:

1. 벌크 구조: 최근 실험들은 저온 상이 $R\bar{3}$ 공간군에 속함을 시사하고 있으나(약 150 K에서 $C2/m$으로 전이), $R\bar{3}$ 구조가 에너지적으로 유리한 저온 상임을 확인하는 선행 이론 연구는 없었습니다.
2. 전자 구조: $\alpha$-RuCl$_3$에서 $J_{\text{eff}} = 1/2$ 모델의 타당성, 특히 넬(Néel) 벡터에 대한 각운동량 양자화 축의 방향에 대한 검증이 필요합니다.
3. 자기 해밀토니안: 기존 연구들은 자기 교환 매개변수(예: 근접 이웃 키타에프 교환 $K_1$의 부호)에서 상당한 불일치를 보고하며, 종종 2차 근접 이웃 상호작용을 간과합니다.
4. $g$-인자 이방성: $g$-인자의 이방성($g_{XY}$ 대 $g_Z$)에 대한 실험적 보고는 작은 이방성부터 큰 이방성까지 서로 모순됩니다. 또한, 이러한 인자를 계산하기 위한 이론적 방법(특히 원자 궤도 투영법)이 불충분할 수 있습니다.
5. 구조적 왜곡: RuCl$_6$ 팔면체의 특정 구조적 왜곡(삼각 왜곡 vs. Cl 삼각형의 상대적 비틀림)이 자기 이방성에 미치는 영향이 체계적으로 분리되어 연구된 바 없습니다.

방법론
저자들은 Quantum ESPRESSO 프레임워크 내에서 국소 밀도 근사(LDA) 및 회전 불변 LDA+U 방법을 사용하여 **제한된 밀도 범함수 이론(cDFT)**을 채택했습니다.

• 구조 완화: $C2/m$ 및 $R\bar{3}$ 적층 구성에 대해 실험적, 고정 셀, 가변 셀 완화 조건 하에서 계산을 수행하여 에너지적 안정성을 결정했습니다.
• $J_{\text{eff}}$ 분석: 파동 함수를 $J_{\text{eff}}$ 기저 상태로 투영하여 전자 구조를 분석했습니다. 결정적으로, $J_{\text{eff}}$ 모델을 테스트하기 위해 각운동량 양자화 축을 넬 벡터와 평행하게 정렬했습니다.
• 자기 해밀토니안 구축: 다양한 자기 구성(zigzag, Néel, stripy, ferromagnet, 그리고 특정 1-by-4 구성)으로부터의 총 에너지를 세 가지 모델에 피팅했습니다:
  1. 전통적인 키타에프-하이젠베르크(KH) 모델.
  2. 스핀 아이스(SI) 모델.
  3. 새로운 이방성 $J_1$-$J_2$-$J_3$ (A-$J_1J_2J_3$) 모델 (이방성 2차 및 3차 근접 이웃 상호작용 포함).
• $g$-인자 계산: 두 가지 방법을 비교했습니다:
  1. 원자 궤도 투영: 모멘트를 허바드 매니폴드(Hubbard manifold)에 투영.
  2. 워니어 보간을 이용한 모멘트 캔팅(Moment Canting): $J_{\text{eff}}$ 모멘트를 $1^\circ$ 기울이고, 변환 불변 워니어 보간(WannierBerri)을 사용하여 총 자기 모멘트의 변화를 계산.
• 왜곡 분석: 두 가지 특정 왜곡(삼각 왜곡: Ru-Cl-Z 각도 변화, 및 상대적 비틀림: XY 평면 내 상하부 Cl 삼각형의 회전)을 격리하여 분석했습니다.

주요 기여 및 결과

1. $R\bar{3}$ 적층 확인:
   본 연구는 다양한 허바드 $U$ 값과 완화 유형에 걸쳐 $R\text{ } \bar{3}$ 구조가 $C2/m$ 구조보다 약 5–6 meV/Ru만큼 낮음을 입증했습니다. 이는 실험적으로 관찰된 저온 $R\bar{3}$ 상에 대한 최초의 이론적 지지를 제공합니다.

2. 정교화된 $J_{\text{eff}}$ 모델:
   전도대 최솟값 상태는 각운동량 양자화 축을 넬 벡터와 평행하게 설정했을 때만 거의 순수한 $J_{\text{eff}} = 1/2$ 및 $m_{\text{eff}} = -1/2$ 특성을 보임을 보여줍니다. 이 관점은 이전에 간과되었으나, 상태의 궤도적 성격을 올바르게 식별하는 데 필수적입니다. $J_{\text{eff}}$ 모멘트의 크기는 하이브리드화와 왜곡로 인해 이상적인 $1/2$보다는 약간 작지만, 여러 구성에서 견고하게 나타났습니다.

3. 자기 해밀토니안 및 교환 매개변수:

• 전통적인 KH 및 스핀 아이스 모델은 특히 Néel 구성에 대해 DFT 총 에너지를 정확하게 재현하는 데 실패합니다.
• 제안된 A-$J_1J_2J_3$ 모델(및 단순화된 A-$J_1J_2$-I-$J_3$ 변형)은 DFT 에너지를 정확하게 재현합니다.
• 2차 근접 이웃 상호작용이 매우 중요합니다. 2차 근접 이웃 교환($K_2, \Gamma_2, \Gamma'_2$)의 이방성 성분은 무시할 수 없으며, 그 크기는 종종 고려되는 1차 근접 이웃의 이방성 항($\Gamma'_1$)보다 큽니다. 이러한 항들을 포함함으로써 기존 모델 대비 총 에너지의 제곱평균제곱근(RMS) 오차를 3.6배 줄였습니다.

4. $g$-인자 이방성 및 계산 방법:

• 모멘트 캔팅 방법과 워니어 보간을 결합한 결과, $g_{XY} \approx 2.22$ 및 $g_Z \approx 1.84$의 $g$-인자를 얻었습니다. 이 값들은 작은 이방성을 나타내며, 이는 최근의 실험 보고(예: Loidl 등)와 일치하며 큰 이방성을 주장했던 이전의 보고들과는 상반됩니다.
• 원자 궤도 투영법은 $g$-인자 크기를 현저히 과소평가하며($g_{XY} \approx 1.54$), 총 자기 모멘트의 무시할 수 없는 부분을 놓칩니다. 저자들은 정확한 궤도 자기 모멘트 계산을 위해서는 워니어 보간 접근법이 필요하다고 결론지었습니다.

5. 구조적 왜곡의 역할:
   왜곡을 격리하여 분석한 결과, 상하부 Cl 삼각형의 상대적 비틀림(XY 평면 내)이 자기 이방성과 올바른 바닥 상태 모멘트 방향($\theta \sim 60^\circ$)을 결정하는 지배적인 요인임을 발견했습니다. 삼각 왜곡만으로는 실험적 바닥 상태를 재현할 수 없습니다. 이 비틀림 왜곡은 외부 섭동 없이 벌크 물질에 존재하는 것으로, 이전의 이론 모델에서는 간과되어 왔습니다.

의의
본 논문은 $\alpha$-RuCl$_3$의 구조 및 자기적 묘사에 대한 오랜 모순을 해결함으로써 이에 대한 포괄적인 제일원리 이해를 제공한다고 주장합니다.

• $R\bar{3}$ 저온 구조를 이론적으로 검증했습니다.
• $J_{\text{eff}}$ 모델을 올바르게 해석하기 위해 양자화 축을 넬 벡터와 정렬해야 함을 확립했습니다.
• 자기 모델링에는 이방성 2차 근접 이웃 교환의 포함이 필수적임을 입증했습니다.
• $g$-인자 계산 방법론을 교정하여, 원자 궤도 투영법은 불충분하며 실제 이방성은 작다는 것을 보여주었습니다.
• Cl 삼각형의 상대적 비틀림이 자기 이방성을 지배하는 핵심적인, 이전에 간과되었던 구조적 매개변수임을 식별하였으며, 이는 키타에프 물질의 자기성을 엔지니어링할 수 있는 새로운 경로를 제시합니다.

저자들은 본 연구가 $\alpha$-RuCl$_3$의 자기성에 대한 다른 관점을 제공하며, 향후 이론 모델과 고압 실험은 정확한 해석을 위해 이러한 특정 구조적 왜곡과 교환 항들을 포함해야 한다고 결론지었습니다.