삼각 격자 세슘·이트륨 화합물의 스핀 해밀토니안과 자기성

초록

본 연구는 DFT+DMFT 기반 Hubbard‑I 접근법을 이용해 Ce와 Yb 삼각 격자 화합물의 결정장 분할, 바닥 상태 Kramers 이중항, 그리고 근거리 및 원거리 교환 상호작용을 계산한다. 얻어진 스핀‑½ 모델을 36 스핀까지 정확 대각화(ED)한 결과, 기존에 양자 스핀 액정(QSL) 후보로 알려진 RbCeO₂와 YbZn₂GaO₅를 포함한 모든 시스템이 120° 반강자성 혹은 스트리프형 반강자성 등 순서된 자기 상태를 보인다. QSL을 안정화하려면 강한 반거리 반강자성 교환이 필요함을 제시한다.

상세 분석

이 논문은 삼각 격자에 위치한 Kramers 이온 Ce³⁺(4f¹)와 Yb³⁺(4f¹³)의 전자구조를 DFT와 동적 평균장 이론(DMFT)을 결합한 Hubbard‑I 근사법으로 정밀히 계산한다. 먼저, 실험 구조를 그대로 사용해 전자 밀도와 4f 궤도 점유를 자기적 비극성(고온) 상태에서 수렴시켰으며, 이를 통해 각 물질의 결정장(CF) 스펙트럼을 얻었다. Ce계에서는 SO 간격의 약 60%에 해당하는 큰 CF 분할이 관측돼 |M=±7/2⟩와 |M=±5/2⟩ 혼합이 눈에 띄게 나타났고, 이는 g‑텐서의 in‑plane 성분(g_ab)이 크게 우세함을 의미한다. 반면 YbZn₂GaO₅(YZGO)에서는 SO 간격이 1.37 eV로 크게 차이나며, CF 혼합이 거의 없고 g‑텐서가 거의 등방성에 가깝다.

교환 상호작용은 FT‑HI(force‑theorem Hubbard‑I) 방법으로 계산했으며, 여기에는 4f 전자의 가상 홉핑에 의한 초전도 교환과 고전적인 쌍극자‑쌍극자 상호작용이 모두 포함된다. 최근접 이웃(NN) 결합은 모두 반강자성(AF)이며, XXZ 형태의 이방성 ∆≥1을 보인다. 특히 Ce계에서는 J±±와 Jz± 같은 비대각 항이 비교적 큰 반면, Yb계에서는 이러한 비등방성 항이 작고, 대신 NNN 결합에서 쌍극자 상호작용이 지배적이다. 결과적으로 YZGO의 NNN 결합은 강한 ∆′≈6을 갖는 고이방성 형태이며, 절대값은 J에 비해 매우 작다(J′/J≈0.003).

이러한 파라미터를 바탕으로 스핀‑½ 해밀토니안을 구축하고, XDiag 라이브러리를 이용해 N=12, 32, 36 사이트 클러스터에 대해 정확 대각화를 수행했다. Ce계(KCeO₂, RbCeO₂)는 120° 반강자성(삼각형 AFM) 패턴을, KCeS₂는 스트리프형(Stripy⊥) 반강자성을 나타냈으며, YZGO 역시 120° AFM 구조를 보였다. 스핀 구조인자 S(k)와 저에너지 스펙트럼은 각각의 순서된 상태와 일치했으며, DSL(Dirac Spin Liquid) 특유의 X점에서의 무갭 싱글톤 모드가 현저히 억제돼 QSL이 실현되지 않음을 확인했다. 파라미터 공간을 확장해 NN 이방성 및 NNN 등방성 교환을 변형한 모델을 조사한 결과, QSL을 안정화하려면 NNN 반강자성 J′가 충분히 커야 하고, NN 이방성(∆, J±±)은 QSL을 억제하는 방향으로 작용한다는 일반적인 경향을 도출했다.

마지막으로, 계산된 교환 상수와 실험적 커리‑와이스 온도(θ)를 비교했을 때, Ce 산화물에서는 θ가 잘 재현되었으나 KCeS₂에서는 이론값이 현저히 작아 실제 교환이 다른 메커니즘(예: 직접 4f‑4f 중첩)으로 강화될 가능성을 제시한다. 전반적으로, 이 연구는 DFT+DMFT 기반 Hubbard‑I 접근법이 희귀‑희토류 삼각 격자 시스템의 실질적인 스핀 해밀토니안을 제공할 수 있음을 입증하고, QSL을 실현하려면 순수한 Heisenberg 삼각 격자 모델을 넘어선 비이상적인 교환 구조(특히 강한 NNN 반강자성)가 필요함을 강조한다.