Influence of ligand field and correlation on the electronic structure of NiO and CoO from DFT+DMFT calculations

전하 자기일관성 DFT+DMFT 계산을 사용하여 본 연구는 결정 구조, 리간드장, 그리고 다양한 상관 강도(산소 2p 상관 포함) 가 상자성 NiO 와 CoO 의 전자 구조 및 스펙트럼 함수에 미치는 영향을 조사한다.

핵심

작은 에너지 넘치는 댄서들 (전자) 이 무대 위를 움직이는 세상을 상상해 보세요. 대부분의 물질에서 이 댄서들은 콘서트장의 군중처럼 자유롭게 움직입니다. 하지만 전이금속 산화물 (산화 니켈과 산화 코발트 등) 이라는 특수한 물질에서는 댄서들이 너무 빽빽하게 모여 서로에 매우 민감하게 반응하기 때문에, 마치 긴밀하게 결속된 집단처럼 행동합니다. 한 명이 움직이면 나머지 모두가 즉각 반응합니다. 이를 "강한 상관관계"라고 하며, 이로 인해 이러한 물질은 까다롭고 매혹적인 방식으로 행동합니다. 때로는 전기를 차단하는 부도체처럼 행동하기도 하고, 도체처럼 행동하지도 않습니다.

제공된 논문은 산화 니켈 (NiO) 과 산화 코발트 (CoO) 라는 두 가지 특정 물질에서 이러한 댄서들이 정확히 어떻게 행동하는지 심층적으로 탐구한 것입니다. 연구자들은 이러한 물질들이 왜 그런 방식으로 행동하는지 파악하고, 컴퓨터 시뮬레이션을 통해 그 행동을 정확하게 예측하는 방법을 찾고자 했습니다.

다음은 그들의 발견을 간단한 비유로 정리한 내용입니다:

1. 문제: "나쁜 지도"

과학자들은 이러한 물질을 매핑하는 표준적인 방법으로 DFT(밀도범함수이론) 를 사용합니다. DFT 를 열린 고속도로를 운전할 때 일반적으로 아주 잘 작동하는 GPS 라고 생각하세요. 하지만 이러한 빽빽한 "댄싱" 물질의 경우, GPS 는 고장 났습니다. 댄서들이 자유롭게 움직일 수 있다고 (금속성) 예측하지만, 실제로는 제자리에 묶여 있습니다 (부도체성). 지도가 잘못된 이유는 댄서들이 끊임없이 서로 부딪히고 반응한다는 사실을 무시했기 때문입니다.

2. 해결책: "자기 수정" 기능이 있는 더 나은 GPS

이를 해결하기 위해 저자들은 DFT+DMFT라는 초고급 방법을 사용했습니다.

• DMFT(동역학적 평균장 이론): 시뮬레이션에 모든 댄서를 지켜보며 이웃들과 실시간으로 어떻게 반응하는지 고려하는 "안무가"를 추가한다고 상상해 보세요. 이렇게 하면 "금속 대 부도체" 문제가 해결됩니다.
• SIC(자기 상호작용 보정): 연구자들은 댄서들이 금속 원자 (니켈 또는 코발트) 일 뿐만 아니라, 그들을 둘러싼 산소 원자들도 춤의 일부라는 사실을 깨달았습니다. 시뮬레이션이 산소 댄서들의 에너지를 중복 계산하지 않도록 산소 댄서들을 위한 특별한 "자기 수정" 규칙을 추가했습니다. 이 조합 (DFT+DMFT+SIC) 은 실제 실험 결과와 거의 완벽하게 일치하는 지도를 만들어냈습니다.

3. 두 개의 무대: 팔면체 대 사면체

연구자들은 두 가지 다른 "무대" (결정 구조) 에서 이러한 물질을 테스트했습니다.

• 염화나트륨형 (RS): 댄서들은 팔면체 (6 개의 꼭짓점이 있는 공과 같은 모양) 로 배열되어 있습니다. 이는 이러한 물질들의 안정된 자연스러운 형태입니다.
• 아연-blende 형 (ZB): 댄서들은 사면체 (4 개의 꼭짓점이 있는 피라미드와 같은 모양) 로 배열되어 있습니다. 이는 연구자들이 댄서들이 다른 환경에 어떻게 반응하는지 보기 위해 인위적으로 만든 불안정한 형태입니다.

발견: 원형 무대와 사각형 무대에서 댄서가 다르게 느끼는 것처럼, 전자는 무대 모양에 따라 에너지 준위를 다르게 나눕니다.

• 팔면체 (RS) 무대에서는 에너지 준위가 한 가지 방식으로 나뉩니다.
• 사면체 (ZB) 무대에서는 준위가 뒤집힙니다.
  이 연구는 팔면체 무대가 더 안정적임을 확인했습니다. 댄서들 (전자) 이 더 낮은 에너지이고 더 편안한 자리에 정착할 수 있기 때문입니다. 사면체 무대는 너무 많은 댄서들을 "반결합" (불편하고 고에너지) 자리에 밀어 넣기 때문에 구조가 불안정해집니다. 이것이 왜 산화 니켈이 자연적으로 팔면체 형태로 발견되지만 사면체 형태로는 발견되지 않는지 설명해 줍니다.

4. "위성"과 "갭"

주요 목표 중 하나는 댄서가 움직이기에 필요한 에너지 점프인 밴드 갭을 측정하는 것이었습니다.

• 실험: 실제 실험 결과 약 5~6 전자볼트 (eV) 의 갭이 관찰되었습니다.
• 구 시뮬레이션: 특별한 보정 없이 컴퓨터는 약 2.5~3 eV 의 작은 갭을 예측했는데, 이는 틀렸습니다.
• 신 시뮬레이션: 산소 원자에 대한 "자기 수정" (SIC) 을 사용함으로써 연구자들의 컴퓨터 모델은 5.1 eV의 갭을 예측했습니다. 이는 현실과 완벽한 일치입니다!

또한 "위성 피크" (데이터 내의 고에너지 돌기) 를 살펴보았습니다. 그들은 모델이 주요 갭을 개선했지만, 위성 피크는 여전히 완벽하게 고정하기에는 다소 까다로웠음을 발견했지만, 전체적인 그림은 이전보다 훨씬 더 명확해졌습니다.

5. 니켈 대 코발트: "한 걸음" 차이

니켈과 코발트는 주기율표에서 이웃입니다. 니켈은 코발트보다 전자 하나를 더 가지고 있습니다.

• 산화 니켈: 매우 안정적이고 고에너지인 "장 - 리 (Zhang-Rice)" 상태 (특수하게 긴밀하게 결합된 댄서 쌍) 를 만드는 특정 수의 댄서를 가지고 있습니다.
• 산화 코발트: 댄서가 하나 적습니다. 이로 인해 "구멍" (빈 자리) 이 생깁니다. 연구자들은 이 missing 댄서가 산화 코발트 물질의 행동을 약간 다르게 만든다는 것을 발견했는데, 에너지 밴드 가장자리에 더 날카롭고 강렬한 피크를 생성합니다. 마치 붐비는 방에서 한 사람을 빼면 나머지 사람들이 빈 공간에 더 강하게 반응하고 이동하는 것과 같습니다.

요약

간단히 말해, 이 논문은 산화 니켈과 산화 코발트에서 전자가 어떻게 춤추는지에 대한 완벽하게 정확한 컴퓨터 시뮬레이션을 구축하는 것에 관한 것입니다.

1. 올바른 답을 얻으려면 금속뿐만 아니라 산소 원자를 고려해야 함을 증명했습니다.
2. 결정의 모양 (팔면체 대 사면체) 이 전자가 에너지를 나누는 방식을 바꾸어 왜 어떤 모양은 안정적이고 어떤 것은 그렇지 않은지 설명했습니다.
3. 실제 실험 결과와 거의 완벽한 일치를 이루어, 이 방법이 이러한 복잡한 물질을 이해하는 신뢰할 수 있는 도구임을 입증했습니다.

저자들은 이 방법이 촉매 (물 분해 또는 연료 제조) 나 에너지 저장과 같은 더 나은 물질을 설계하려는 과학자들에게 강력한 도구라고 결론지었습니다. 전자의 "춤"을 이해하는 것이 그 잠재력을 unlocking 하는 열쇠이기 때문입니다.

기술 요약

기술적 요약: DFT+DMFT 계산을 통한 NiO 및 CoO 의 전자 구조에 대한 리간드 장과 상관 효과의 영향

문제 제기
NiO 와 CoO 와 같은 전이 금속 산화물 (TMO) 은 전자 상관, 리간드 장, 결정 구조, 그리고 혼성화의 상호작용에 의해 구동되는 복잡한 전자적 거동을 보입니다. 국소 근사 (LDA) 나 일반화 기울기 근사 (GGA) 를 포함한 표준 밀도 범함수 이론 (DFT) 은 이러한 강상관 계의 절연체 바닥상태와 정확한 여기 스펙트럼을 재현하지 못하므로, 더 정교한 방법이 필요합니다. 구체적으로, 결정 구조 (염화나트륨형 대 아연-blende 형), 조성 (Ni 대 Co), 그리고 전이 금속 (TM) $3d$ 및 산소 $2p$ 오비탈에서의 상관 처리가 어떻게 스펙트럼 함수와 밴드 갭에 영향을 미치는지 체계적으로 이해할 필요가 있습니다. 이전 연구들은 종종 특정 상호작용 매개변수에 초점을 맞추거나, 전하 이동 절연체를 정확하게 기술하는 데 필수적인 산소 리간드에서의 상관을 간과해 왔습니다.

방법론
저자들은 산소 오비탈에 대한 자기 상호작용 보정 (SIC) 체계가 추가된 (DFT+sicDMFT) 밀도 범함수 이론과 동적 평균 장 이론 (DFT+DMFT) 의 전하 자기 일관적 결합을 사용합니다.

• 계산 프레임워크: DFT 구성 요소는 노름 보존赝퍼텐셜을 사용하는 혼합 기저赝퍼텐셜 코드 (MBPP) 를 활용합니다. 상관 서브스페이스는 TM 이온에 중심을 둔 원자 유사 $d$-오비탈로 정의되며, 단축된 Kohn-Sham 기저 내에서 분해됩니다.
• 상관 처리: DMFT 임피urity 문제는 CTHYB(TRIQS 패키지) 인 연속 시간 양자 몬테카를로 (QMC) 솔버를 사용하여 해결됩니다. TM $3d$ 서브스페이스 내의 전자 - 전자 상호작용은 허바드 매개변수 $U$와 훈드 교환 $J$를 가진 회전 불변 슬레이터 해밀토니안을 사용하여 모델링됩니다.
  • NiO 의 경우: $U = 10$ eV, $J = 1$ eV.
  • CoO 의 경우: $U = 9$ eV ($U=6$ eV 로도 테스트됨), $J = 1$ eV ($J=0.9$ eV 로도 테스트됨).
• 리간드 상관: TM 상태와 혼성화된 O($2p$) 오비탈의 상관을 고려하기 위해, 저자들은 산소赝퍼텐셜에 SIC 체계를 적용합니다. 매개변수는 $\alpha=0.8$로 설정되며, O($2s$) 에 대해서는 완전 보정 ($w_l=1.0$), O($2p$) 에 대해서는 80% 보정 ($w_l=0.8$) 을 적용합니다.
• 연구 대상 시스템: 상자성 NiO 와 CoO 는 입방정 염화나트륨형 (RS) 과 아연-blende 형 (ZB) 구조 모두에서 분석됩니다. RS 는 안정한 상이지만, ZB 는 리간드 장 효과 (팔면체 대 사면체 배위) 를 분리하기 위한 모델 시스템으로 연구됩니다.
• 분석: 실수 주파수 스펙트럼 함수는 해석적 연속 (최대 엔트로피 방법) 을 통해 유도됩니다. 연구는 총, 원소 투영 (TM 및 O), 그리고 국소 오비탈 분해 스펙트럼 함수에 초점을 맞추며, 특히 중심 (COGs) 과 밴드 갭 특성을 분석합니다.

주요 결과

1. SIC 와 상호작용 강도 ($U$) 의 역할:

   • RS-CoO 의 경우, $U$를 6 eV 에서 9 eV 로 증가시키면 저에너지 $d^7\underline{L}$ 피크 (여기서 $\underline{L}$은 리간드 홀) 가 $-1.0$ eV 에서 $-1.6$ eV 로 이동합니다.
   • O($2p$)-SIC 를 적용하면 이 피크가 반대 방향으로 $-0.6$ eV 로 이동합니다. $U=9$ eV 와 O($2p$)-SIC 의 조합은 O($2p$) 밴드에 대해 $-3.1$ eV 인 스펙트럼 모양과 피크 위치를 제공하여 UPS 실험 데이터와 잘 일치합니다.
   • 결정적으로, 충분히 큰 $U$에서만 유의미한 밴드 갭이 열립니다. O($2p$)-SIC 를 포함하면 비점유 상부 허바드 밴드 (UHB) 가 더 높은 에너지로 이동하여 RS-CoO 에 대해 계산된 갭이 5.1 eV 가 됩니다. 이 값은 5–6 eV 범위의 실험적 범위 내에 있으며, O($2p$)-SIC 없이 수행된 이전 DFT+DMFT 연구들이 보고한 2.5–3.8 eV 까지 낮은 갭 값에 비해 현저한 개선입니다.

2. 구조적 및 리간드 장 효과 (RS 대 ZB):

   • 계산은 결정장 (CF) 이론 예측을 확인합니다. RS(팔면체) 에서는 $t_{2g}$ 상태가 $e_g$보다 에너지가 낮습니다. ZB(사면체) 에서는 계층이 반전됩니다 ($t_2$보다 $e$가 낮음).
   • 적분된 스펙트럼 함수를 통한 점유 분석은 고스핀 배치를 확인합니다: RS-NiO 의 경우 $t_{2g}^6 e_g^2$, RS-CoO 의 경우 $t_{2g}^5 e_g^2$입니다. ZB 구조의 경우 점유는 NiO 에 대해 $e^4 t_2^4$, CoO 에 대해 $e^4 t_2^3$입니다.
   • ZB 구조는 반결합 상태의 점유로 인해 덜 안정한 것으로 나타납니다. ZB 에서 $t_2$ 오비탈은 반결합성이며, 이러한 상태에 존재하는 3 개 (CoO) 와 4 개 (NiO) 의 전자는 RS 상에 비해 구조를 불안정하게 만듭니다.

3. 스펙트럼 특징 및 차폐:

   • CoO 스펙트럼은 $d$-껍질 내 추가 홀로 인해 NiO 에 비해 더 날카로운 피크와 증가된 저에너지 스펙트럼 무게를 보입니다.
   • CoO 의 가전자대 가장자리 근처에 있는 날카로운 피크는 Li 도핑된 NiO 에서 논의된 $d^8\underline{Z}$ 거동과 유사한 $d^7\underline{Z}$ 상태 (여기서 $\underline{Z}$는 Zhang-Rice 홀을 나타냄) 로 해석됩니다.
   • 점유된 O($2p$) 와 비점유 TM($3d$) 상태 사이의 에너지 차이로 정의된 전하 이동 에너지 ($\Delta$) 는 구조와 원소에 따라 4.3–4.9 eV 로 계산되며, 실험적 경향 (NiO < CoO) 과 일치하는 계층을 보입니다.

의의 및 주장
이 논문은 DFT+sicDMFT 접근법이 강상관 TMO 의 전자 구조를 기술하는 강력한 프레임워크를 제공한다고 주장합니다. SIC 를 통해 산소 $2p$ 오비탈의 상관 효과를 명시적으로 포함함으로써, 이 방법은 정적 평균장 접근법 (DFT+U 등) 이나 표준 DFT+DMFT 가 종종 정확하게 포착하지 못하는 상자성 NiO 와 CoO 의 실험적 여기 스펙트럼 및 밴드 갭 값을 성공적으로 재현합니다.

이 연구는 다음을 보여줍니다:

• 허바드 $U$ 매개변수와 리간드 상관 사이의 상호작용은 실험과의 정량적 일치에 중요합니다.
• 이 방법론은 스펙트럼 함수에 대한 결정 구조 (리간드 장) 와 조성의 효과를 체계적으로 분리할 수 있습니다.
• 계산된 스펙트럼 특징, 예를 들어 밴드 갭과 전하 이동 에너지는 전자 구조가 정확하게 모델링된다면 촉매 및 에너지 변환 응용을 위한 TMO 최적화를 지원할 수 있는 신뢰할 수 있는 기술자 역할을 합니다.

저자들은 겸손한 입장을 유지하며, ZB-NiO 는 안정한 벌크 상이 아니지만, 리간드 장 효과를 체계적으로 비교할 수 있도록 모델 시스템으로서 연구되었다고 명시합니다. 이 작업은 NiO 에 대한 이전 발견들을 바탕으로 CoO 에 대한 접근법을 검증하고 구조적 변이에 대한 분석을 확장합니다.