Entropy-Driven Structural Phase Transition in Nb$_3$Cl$_8$ via Density Functional Theory and an Effective Model

기술 요약: DFT 및 유효 모델을 통한 Nb$_3$Cl$_8$의 엔트로피 주도 구조적 상전이

문제 제기
Nb$_3$Cl$_8$은 유효 삼각 격자 위에서 전형적인 플랫 밴드 클러스터 모트 절연체(flat-band cluster Mott insulator)로, 양자 스핀 액체(QSL) 상태를 구현할 수 있는 잠재적 후보 물질로 간여됩니다. 그러나 이 물질의 저온 물리학은 약 90 K에서 발생하는 1차 구조적 상전이에 의해 가려져 있습니다. 이 전이는 고온의 상자성 $\alpha$ 상을 저온의 비자성 $\beta$ 상으로 변형시키며, 이는 후보 QSL 영역을 억제합니다. 이 전이의 미시적 기원은 전하 불균등화(charge disproportionation) 또는 층간 다이머 싱글렛(interlayer dimer singlet) 형성이라는 상충하는 설명들 사이에서 논쟁이 되어 왔습니다. 또한, 열역학적 구동력과 외부 자극(예: 단축 압력)이 이 전이를 억제하는 메커니즘 역시 완전히 이해되지 않았습니다.

방법론
저자들은 격자, 전자, 스핀 기여를 분리하기 위해 제일원리 밀도 범함수 이론(DFT)과 확장된 허버드 모델(extended Hubbard model)을 결합한 통합 자유 에너지 프레임워크를 사용합니다.

1. DFT 계산: 비스핀 분극 DFT 계산(optB88-vdW 및 PBE 범함수를 사용한 VASP 활용)을 수행하여 $\alpha$ (공간군 $P\bar{3}m1$) 및 $\beta$ (Represented by $R\bar{3}m$) 상에 대한 완화된 구조, 바닥 상태 에너지 및 포논 스펙트럼을 얻었습니다. 두 구조 사이의 에너지 장벽을 추정하기 위해 선형 보간 경로를 구축했습니다.
2. 포논 열역학: 조화 근사(harmonic approximation) 내에서 밀도 범함수 섭동 이론(DFPT)을 사용하여 포논 자유 에너지($F_{\text{phonon}}$)와 엔트로피를 계산했습니다.
3. 유효 모델: 반 채워진 2$a_1$ 매니폴드의 모트 물리학을 다루기 위해, 클러스터 중심 궤도들의 워너(Wannier) 기저를 기반으로 확장된 허버드 모델을 구축했습니다. 파라미터(호핑 $t$, 온사이트 $U$, 및 격자 간 $V$ 쿨롱 상호작용)는 제약된 랜덤 위상 근사(cRPA)를 통해 도출되었습니다.
4. 이중 계수 보정(Double-Counting Correction): DFT 참조에서 이미 포착된 저에너지 2$a_1$ 기여를 빼기 위해 하트리 수준의 이중 계수 보정($E_{\text{DC}}$)을 구현하여, 일관된 총 자유 에너지($F_{\text{tot}}$)를 확보했습니다.
5. 스핀 열역학: 확장된 허버드 모델은 유효 하이젠베르크 모델로 축소되었습니다. $\alpha$ 상의 경우, 연결된 클러스터 고온 급수 전개(HTSE)를 사용하여 스핀 엔트로피를 계산했습니다. $\beta$ 상의 경우, 강한 층간 결합을 고려하여 단순화된 고립 다이머 모델을 채택했습니다.
6. 압력 효과: 응력-일(stress-work) 항을 포함하여, $c$축 단축 압력에 따른 효과를 평가하기 위해 $T=0$에서 일반화된 엔탈피 프레임워크($H_{\text{tot}}$)를 구축했습니다.

주요 결과

• 구조적 장벽: DFT 계산 결과, $\alpha$ 상과 $\beta$ 상 사이에 약 300 meV/2-site unit cell의 상당한 에너지 장벽이 존재함을 밝혔으며, 이는 실험적으로 관찰된 열적 이력 현상(thermal hysteresis)과 일치합니다.
• 열역학적 구동력: 전이는 내부 에너지와 엔트로피 사이의 경쟁에 의해 구동됩니다:
  • $\beta$ 상의 안정성: 저온 $\beta$ 상은 층간 다이머화에 의한 내부 에너지 이득에 의해 안정화됩니다. 이는 강한 층간 호핑($t^{(\text{inter})}_0 \approx -136.7$ meV)을 수반하며, 큰 반강자성 교환 상호작용($J \approx 70$ meV)과 스핀 싱글렛 형성을 유도하여 스핀 엔트로피를 억제합니다. 또한, $\beta$ 상의 더 강한 층간 결합은 포논 모드를 경화(harden)시켜 포논 엔트로피를 감소시킵니다.
  • $\alpha$ 상의 안정성: 고온 $\alpha$ 상은 엔트로피에 의해 안정화됩니다. 이 상은 더 부드러운 포논 모드(낮은 영점 에너지를 가지나 높은 엔트로피를 보유)를 가지며, 약한 층간 결합 덕분에 상자성 스핀 엔트로피를 유지합니다.
• 전이 온도: 계산된 총 헬름홀츠 자유 에너지의 교차점은 $T_c \approx 113$ K에서 발생하며, 이는 실험적 전이 온도인 $\sim$90 K와 반정량적으로 일치합니다. 계산된 전이 전후의 엔트로피 변화는 거시적 비열 측정값과 일치합니다.
• 압력 응답: $c$축 단축 압력 하에서의 일반화된 엔탈피 계산에 따르면, 구조적 반응의 차이로 인해 $\alpha$ 상이 $\beta$ 상에 비해 엔탈피 감소 측면에서 덜 압축적임을 나타냅니다. 계산 결과, 전이는 약 2.6 GPa의 임계 압력에서 완전히 억제되어 저온에서도 $\alpha$ 상을 안정화시킵니다.

의의 및 주장
본 논문은 Nb$_3$Cl$_8$의 구조적 전이에 대한 통합적인 열역학적 그림을 제공하며, 그 미시적 기원이 단일 메커니즘이 아닌 포논과 스핀 엔트로피의 공동 구동임을 밝혀냈다고 주장합니다. 격자, 스핀, 상관 관계 기여를 명시적으로 분리함으로써, 본 연구는 $c$축 단축 압력에 의한 전이 억제가 일반화된 엔탈피 지형에서 $\alpha$ 상을 안정화시키는 결과임을 입증했습니다. 이는 압력이 전이를 억제한다는 실험적 관찰에 대한 열역학적 설명을 제공하며, 저온에서 $\alpha$ 상의 후보 QSL 영역에 접근할 수 있는 실험적 경로를 제시할 수 있습니다. 저자들은 포논에 대한 조화 근사 및 압력 하에서의 아조화(anharmonicity) 및 유한 온도 재규격화의 무시와 같은 한계를 인정하며, 결과가 반정량적임을 명시하였습니다.