Significance of the dispersion force for ferroelectric switching in ZnO and related materials

기술 요약: ZnO 및 관련 물질의 강유전성 스위칭에 있어서 분산력의 중요성

문제 제기
우르츠 구조(Wurtzite-structured)를 가진 산화아연(ZnO)은 넓은 밴드갭을 가진 극성 물질이지만, 이 물질의 고유한 강유전성 스위칭 장벽은 실용적인 응용을 적용하기에는 지나치게 높다. Zn$_{0.5}$Mg$_{0.5}$O와 같이 치환된 물질들은 스위칭 가능한 분극을 나타내지만, 이러한 협동적 분극 스위칭을 지배하는 원자 수준의 메커니즘에 대해서는 여전히 논쟁이 진행 중이다. 특히, 평면 육방정계 중간상(h-ZnO 또는 구조 C)의 안정성은 논란의 여지가 있다. 메타스테이블(metastable)한 h-ZnO 나노 결정에 대한 실험적 관찰 결과는, 무한 결정 구조에서 C를 불안정한 전이 상태로 식별하는 많은 계산적 예측과 상충한다. 본 연구의 핵심 가설은 표준적인 계산 접근법들이 ZnO나 MgO와 같은 "경질 원자(hard atom)" 시스템에서 흔히 무시되는 반데르발스 분산력을 부적절하게 처리했을 가능성이 있으며, 이것이 중간상의 안정성에 관한 이론과 실험 사이의 불일치를 초导致할 수 있다는 것이다.

방법론
저자들은 ZnO와 Zn$_{0.5}$Mg$_{0.5}$O의 협동적 강유전성 스위칭을 위한 퍼텐셜 에너지 표면(PES)을 매핑하기 위해 종합적인 ab initio 및 밀도범함수 이론(DFT) 방법들을 사용하였다. 스위칭 경로는 두 가지 우르츠 구조(A와 A') 사이를 평면 육방정계 중간상(C)과 잠재적 전이 상태(B)를 거쳐 이동하는 과정으로 모델링되었다.

계산 접근 방식은 다음과 같다:

1. 고수준 파동함수 방법: 하트리-포크(HF), 2차 뫼러-플렉스(MP2) 섭동 이론, 그리고 랜덤 위상 근사(RPA). 이 방법들은 전자 상관관계와 분산을 다르게 처리하는 벤치마크를 제공하기 위해 선택되었으며, MP2와 RPA는 HF보다 체계적인 개선을 제공한다.
2. 밀도범함수 이론 (DFT): PBESol, PBE0, HSE06, r2SCAN을 포함한 다양한 범함수.
3. 분산 보정: 표준 범함수에서 장거리 분산의 누락 문제를 해결하기 위해, HF, PBESol, PBE0, HSE0 się HSE06에 대해서는 D3(BJ)를 적용하였고, r2SCAN에는 rVV10을 적용하였으며, 기타 보정법(D2, D3, ULG, optB88)을 테스트하였다.
4. 구조 분석: 4-원자 단위 격자(Zn$_2$O$_2$ 및 ZnMgO$_2$)에 대해 계산을 수행하였다. 반응 좌표는 c-축을 따른 양이온의 분율 변위로 정의되었다. 전체 구조 최적화와 특정 격자 벡터(평면 구조에 적용된 우르츠 격자 벡터, C$_v$로 명명됨)에 대한 구속 최적화를 모두 수행하였다.
5. 수렴성: MP2 및 RPA의 정량적 정확성을 확보하기 위해, 특히 k-포인트 그리드(최대 12 $\times$ 12 $\times$ 8)와 평면파 컷오프(최대 700 eV)에 대한 광범위한 테스트를 수행하였다.

주요 결과

• h-ZnO (ZnO)의 안정성: MP2 및 RPA 계산의 합의된 결과는 평면 육방정계 구조 C가 무한 ZnO 결정에서 메타스테이블한 상이 아니라 전이 상태라는 것이다. 일부 DFT 범함수(예: PBE0/D3, HSE06/D3)에서의 분산 보정이 C를 인위적으로 안정화시켜 메타스테이블하게 보이게 만들기도 하지만, 고수준 방법들은 여전히 C가 불안정함을 보여준다. 이는 실험적으로 관찰되는 안정적인 h-ZnO 나노 결정이 내부적 벌크 안정성보다는 유한 크기 효과나 용매 효과에 의한 것임을 시사한다.
• h-Zn$_{0.5}$Mg$_{0.5}$O의 안정성: 대조적으로, Mg가 치환된 시스템의 경우, MP2와 RPA는 평면 육방정계 상 C가 적어도 메타스테이블하며, 어떤 경우에는 우르츠 상 A보다 더 안정적이라고 일관되게 예측한다. 그러나 표준 분산 보정(D3(BJ), rVV10)을 사용하는 DFT 방법들은 MP2 및 RPA에 비해 C의 안정성을 크게 과대평가하는 경傾向이 있다. 오직 r2SCAN/rVV10 접근법만이 A가 C보다 약간 더 안정적이라고 예측하여, 고수준 벤치마크와 더 가깝게 일치한다.
• 분산의 역할: 분산력은 우르츠 격자 벡터에 구속된 가상의 평면 육방정계 구조(C$_v$)에서 가장 유의미한 것으로 나타났다. 최적화된 구조의 경우, 분산 보정은 일반적으로 A에 비해 C를 안정화시키지만, 이 효과의 크기는 서로 다른 DFT 범함수 간에 매우 다양하게 나타난다.
• 방법론 성능:
  • MP2 및 RPA는 이온, 공유 및 분산 상호작용에 대한 자가 일관적이고 균형 잡힌 설명을 제공하며, 가장 신뢰할 수 있는 벤치마크 역할을 한다.
  • 분산이 없는 DFT는 ZnO에 대해서는 MP2/RPA와 질적으로 일치하지만, Zn$_{0.5}$Mg$_{0.5}$O에 대해서는 실패한다.
  • 경험적 분산 보정 (특히 D3(BJ) 및 rVV10을 사용하는 DFT)은 종종 평면상을 과도하게 안정화시켜 고수준 이론과 질적인 불일치를 야기한다. ULG 보정과 optB88 범함수는 D3(BJ)보다 MP2/RPA 결과와 더 나은 일치도를 보였다.
• 격자 상수: 분산 보정은 에너지론에는 결정적이지만, 평형 격자 상수를 예측하는 데는 미미한 역할을 하며, 대부분의 방법에서 실험값과의 편차는 작게 유지되었다.

의의 및 주장
본 논문은 "경질 원자"로 구성되어 분산이 전통적으로 무시되는 물질이라 할지라도, 우르츠 구조 물질의 분극 스위칭과 강유전성 거동을 기술하는 데 있어 정확한 분산력 처리가 필수적이라고 주장한다. 본 연구는 다음을 입증한다:

1. 인지된 h-ZnO의 안정성은 분산 보정의 선택에 따라 크게 영향을 받으며, 고수준 이론(MP2/RPA)으로부터 h-ZnO가 벌크에서 불안정하다는 합의를 도출함으로써, 환경적 요인에 의한 것임을 밝혀 실험적 나노 결정 데이터와의 충돌을 해결하였다.
2. Zn$_{0.5}$Mg$_{0.5}$O의 경우, 스위칭 장벽의 감소와 중간상의 안정화는 실제 현상이지만, 그 정량적 기술은 계산 방법에 매우 민감하다. 경험적 분산 보정을 사용하는 표준 DFT는 이 시스템에 대해 신뢰할 수 없는 질적 예측을 내놓는 경우가 많다.
3. MP2와 RPA는 이러한 시스템에서 진정한 메타스테빌리티와 전이 상태를 구분할 수 있는 필요한 정량적 정확성을 제공하며, 분산을 무시하기 위한 "경질 원자" 가정이 강유전성 스위칭 메커니즘을 이해하는 데 불충분함을 강조한다.

저자들은 분산이 스위칭 에너지론에 중요한 기여를 하지만, 그것이 Mg 치환에 의해 관찰되는 장벽 감소를 단독으로 설명하는 것은 아니라고 결론짓는다. 즉, 치환에 의해 유도된 이온 및 공유 결합 변화 사이의 상호작용 또한 결정적이며, 이를 위해서는 전자 상관관계를 일관되게 다루는 방법이 필요하다.