[论文解读] 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 more|ArXiv.org|Oct 21, 2009

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.

研究动机与目标

• 评估多种电子结构方法（力场、紧束缚、DFT和QMC）在描述硅体系势能面（PES）方面的准确性。
• 评估力场在预测扩散系数等动力学性质方面的可靠性，这些性质依赖于能垒高度和PES的拓扑结构。
• 确定广泛使用的力场是否忠实再现了硅体系（特别是缺陷和非晶结构）的构型态密度（C-DOS）和势阱分布。
• 比较多种力场（Tersoff、Stillinger-Weber、EDIP、Lenosky）的性能，并识别其PES表征中的系统性误差。
• 为大规模硅体系模拟（尤其是涉及动力学或亚稳态性质时）提供方法选择的基准。

提出的方法

• 采用量子蒙特卡罗（QMC）方法，特别是变分蒙特卡罗（VMC）和扩散蒙特卡罗（DMC），作为硅体系能量景观的准精确参考。
• 采用最小值跳跃法（MHM）对势能面进行采样，并计算最小值跳跃态密度（MH-DOS），以估算单位能量区间内的局部最小值数量。
• 从MHM数据中推导出构型态密度（C-DOS），表示不同能量层级上稳定构型的真实分布。
• 比较多种电子结构方法：DFT（LDA和PBE）、紧束缚（Lenosky T-B）以及四种力场（Tersoff、Stillinger-Weber、EDIP、Lenosky），所有方法均使用相同的调用序列以确保公平比较。
• 研究采用1000个原子的硅团簇，并进行长时间的MHM运行（最多1,000,000个构型），以确保高能量非晶和剪切结构的收敛与可靠采样。
• 通过优化试探波函数，最小化DMC中的固定节点误差，所得DMC结果用作准确性评估的基准。

实验结果

研究问题

• RQ1 与准精确的QMC参考相比，标准DFT泛函（LDA、PBE）在描述硅体系能垒高度方面的准确性如何？
• RQ2 广泛使用的经典力场（Tersoff、Stillinger-Weber、EDIP、Lenosky）在多大程度上再现了硅的真实构型态密度（C-DOS）？
• RQ3 为何使用不同力场对复杂硅体系（如缺陷和非晶相）进行模拟会得到显著不同的结果？
• RQ4 力场是否正确描述了势能面上局部最小值和鞍点的分布，特别是在非周期性和缺陷体系中？
• RQ5 DFT与力场生成的能量景观在预测硅团簇全局最小值和亚稳态方面，其表现如何比较？

主要发现

• DFT，特别是LDA泛函，能为硅重排过程提供高度准确的能垒高度，其误差与QMC相当，因此DFT在动力学模拟中具有可靠性。
• 最小值跳跃法（MHM）显示，Tersoff、Stillinger-Weber和EDIP等力场生成的MH-DOS和C-DOS与QMC相比存在显著差异，表明其势能面存在重大误差。
• 除基于MEAM的Lenosky力场外，所有力场在非晶和缺陷构型中均表现出过于崎岖的能量景观，伴随大量虚假的局部最小值和鞍点。
• 从力场计算得到的C-DOS存在显著偏差：Tersoff和EDIP在约10 eV处产生单一峰，而Stillinger-Weber和Lenosky T-B在约35 eV处出现第二个峰，这是由于对高能剪切结构的采样所致，表明收敛性差且采样不真实。
• 不同力场的MH-DOS在高能区显著发散，证实其C-DOS不可靠，因此在可能的情况下，使用单一力场的模拟结果应通过DFT进行交叉验证。
• 研究结论认为，由于力场对PES的描述较差，其不适用于硅体系的精确动力学模拟，而DFT为大规模模拟提供了稳健且准确的替代方案。