[논문 리뷰] The energy landscape of silicon systems and its description by force fields, tight binding schemes, density functional methods and Quantum Monte Carlo methods

S. Alireza Ghasemi, Maximilian Amsler외 8명|ArXiv.org|2009. 10. 21.

Advanced Chemical Physics Studies인용 수 4

한 줄 요약

이 연구는 힘장, 타이트 본드, 밀도함수이론(DFT), 양자 몽테카를로(QMC) 방법이 실리콘 시스템의 에너지 표면을 기술하는 데 있어 정확도를 평가한다. QMC를 준정확한 기준으로 사용하여, DFT, 특히 LDA가 장벽 높이와 위치 에너지 표면을 신뢰할 만하게 제공하는 것으로 밝혀졌으며, 힘장은 다수의 임의의 최소점과鞍점(사슬점)을 포함한 과도하게 거친 표면을 생성하여 실리콘의 동역학적 시뮬레이션에 사용하는 데에 부적합함을 밝혀냈다.

ABSTRACT

The accuracy of the energy landscape of silicon systems obtained from various density functional methods, a tight binding scheme and force fields is studied. Quantum Monte Carlo results serve as quasi exact reference values. In addition to the well known accuracy of DFT methods for geometric ground states and metastable configurations we find that DFT methods give a similar accuracy for transition states and thus a good overall description of the energy landscape. On the other hand, force fields give a very poor description of the landscape that are in most cases too rugged and contain many fake local minima and saddle points or ones that have the wrong height.

연구 동기 및 목표

• 실리콘 시스템의 위치 에너지 표면(PES)을 기술하는 데 있어 힘장, 타이트 본드, DFT, QMC 등의 전자 구조 방법의 정확도를 평가하는 것.
• 장벽 높이와 PES의 기하학적 구조에 따라 영향을 받는 확산 계수 등의 동역학적 성질을 예측하는 데 있어 힘장의 신뢰성을 평가하는 것.
• 특히 결함과 비정질 구조에서의 구성도 상태(C-DOS)와 최소점 분포를 충실하게 재현하는지 여부를 평가하는 것.
• 여러 힘장(Tersoff, Stillinger-Weber, EDIP, Lenosky)의 성능을 비교하고, 그들의 PES 표현에 대한 체계적 오류를 규명하는 것.
• 특히 동역학적 또는 준안정 상태 성질에 관심이 있는 대규모 실리콘 시뮬레이션에서 적절한 방법을 선택하기 위한 기준을 제공하는 것.

제안 방법

• 양자 몽테카를로(QMC) 방법, 특히 변분 몽테카를로(VMC)와 확산 몽테카를로(DMC)는 실리콘 시스템의 에너지 표면에 대한 준정확한 기준으로 사용된다.
• 최소점 이동 방법(Minima Hopping Method, MHM)을 사용하여 위치 에너지 표면을 샘플링하고, 최소점 이동 밀도 상태(MH-DOS)를 계산한다. 이는 에너지 간격당 국소 최소점의 수를 추정한다.
• 구성도 밀도 상태(C-DOS)는 MHM 데이터에서 유도되며, 에너지 수준에 따른 안정한 구조의 진짜 분포를 나타낸다.
• 여러 전자 구조 방법을 비교한다: DFT(LDA 및 PBE), 타이트 본드(Lenosky T-B), 그리고 네 가지 힘장(Tersoff, Stillinger-Weber, EDIP, Lenosky)이며, 공정한 비교를 위해 동일한 호출 순서로 구현된다.
• 연구는 1000원자 실리콘 클러스터를 사용하고, 고에너지 비정질 및 비틀림 구조를 충분히 수렴시키기 위해 장기간의 MHM 실행(최대 1,000,000개의 구조)을 수행한다.
• DMC의 고정 노드 오차는 시험 파동함수를 최적화하여 최소화하고, 그 결과로 얻은 DMC 결과를 정확도 평가의 기준으로 사용한다.

실험 결과

연구 질문

• RQ1 표준 DFT 기능성(LDA, PBE)이 QMC 준정확 기준과 비교하여 실리콘 시스템의 장벽 높이를 얼마나 정확하게 기술하는가?
• RQ2 일반적으로 사용되는 고전적 힘장(Tersoff, Stillinger-Weber, EDIP, Lenosky)이 실리콘의 진정한 구성도 밀도 상태(C-DOS)를 어느 정도 정확하게 재현하는가?
• RQ3 다양한 힘장이 결함과 비정질 상과 같은 복잡한 실리콘 시스템에서 상당히 다른 결과를 낳는 이유는 무엇인가?
• RQ4 힘장은 비정기적 및 결함이 있는 시스템에서 위치 에너지 표면의 국소 최소점과 사슬점 분포를 올바르게 표현하는가?
• RQ5 DFT와 힘장이 실리콘 클러스터에서 전역 최소점과 준안정 상태를 예측하는 데 있어 에너지 표면의 능력은 어떻게 비교되는가?

주요 결과

• DFT, 특히 LDA 기능성이 실리콘 재배열 과정의 장벽 높이를 매우 정확하게 제공하며, QMC와 유사한 오차 수준을 보이므로 동역학적 시뮬레이션에 신뢰할 수 있다.
• 최소점 이동 방법(MHM)은 Tersoff, Stillinger-Weber, EDIP 등의 힘장이 QMC와 비교해 극명하게 다른 MH-DOS와 C-DOS를 생성하는 것으로 밝혀져, 그들의 위치 에너지 표면에 심각한 오류가 있음을 시사한다.
• Lenosky 잠재력 기반 MEAM를 제외한 모든 힘장은 비정질 및 결함이 있는 구성에서 다수의 임의의 국소 최소점과 사슬점을 포함한 과도하게 거친 에너지 표면을 생성한다.
• 힘장에서 계산된 C-DOS에는 중대한 차이가 있다: Tersoff와 EDIP는 약 10 eV 부근에 단일 피크를 보이며, Stillinger-Weber와 Lenosky T-B는 고에너지 비틀림 구조를 샘플링하면서 약 35 eV 부근에 두 번째 피크를 보여, 수렴 불량과 물리적으로 비합리적인 샘플링을 시사한다.
• 다른 힘장에서 유도된 MH-DOS는 특히 고에너지 영역에서 크게 다름을 보이며, 이는 그들의 C-DOS가 신뢰할 수 없음을 확인하고, 가능한 한 DFT와 함께 단일 힘장으로의 시뮬레이션을 교차 검증해야 함을 시사한다.
• 연구는 힘장이 PES를 잘못 기술함으로써 실리콘의 정확한 동역학적 시뮬레이션에 부적합하다고 결론내리며, DFT는 대규모 시뮬레이션에 대해 강력하고 정확한 대안을 제공한다.