LAMMPS学习系列（11）

实例一

在用LAMMPS进行分子动力学模拟之前，最重要的一件事就是找到合适的相互作用势，那么如何判断所得到的经验势是否合适呢？今天就结合一个简单的实例介绍一下如何去判断：

直观判断：做个简单的动态弛豫模拟，看看体系是否稳定，结构是否变混乱，具体弛豫模拟可见我之前的推文：LAMMPS学习系列（1）；
能量判断（推荐）：计算相关的缺陷（空位，间隙子。。）形成能，层错能，迁移能等能量值或者弹性常数等力学值，然后与DFT计算结果进行对比；


一般来说，判断一个经验势好不好的标准就是它的计算值与DFT的计算值的吻合程度是否很高，但是DFT可以计算的能量值有很多，经验势往往可能只能与其中一部分值吻合的比较好，所以这也是在一定程度上限制了经验势的使用范围。例如上述文章，只显示了缺陷形成能以及迁移能这部分能量与DFT计算值吻合的很好，但是并没有提及到弹性常数之类的力学值，所以这个W-He经验势一般只用来做缺陷演化之类的模拟，很少应用于力学性能变化的研究，当然，没提及并不代表它在这方面吻合的不好，只是做之前需要自行验证。

接下来先介绍一个BCC结构W的空位形成能的计算，其计算流程如下：</