大家好,我是小马老师!
最近在咨询过程中,发现不少朋友对势函数混合势写法不是很清楚,以致运行出错。
比较常见的一个错误是:All pair coeffs are not set,如下图所示:
错误提示为“不是所有的原子都设置了势函数”,表示有的原子可能没有设置势参数。
出现这种错误主要有两种原因:
第一,因为原子数量多,确实少写了某种原子的势参数,这种情况仔细检查,找出未设置势的原子,补充设置即可。
第二,设置了全部原子的势函数,但是设置方法出错,这种情况再怎么检查原子势也不能解决问题。
因为表面上看,所有的势已经全部被设置。
下面以一个简单的Cu-Al界面模型为例,给出这种错误的解决方法。
界面模型如下图所示,下层为Cu,原子类别1,上层为Al,原子类型2,Cu和Al均使用meam势,cu-Al之间的势使用lj势。