lammps案例：空位形成能的模拟与计算

大家好，我是小马老师。
在昨天的推文中，介绍了计算Cu间隙原子形成能的lammps模拟方法。
与间隙原子相对应的另一种缺陷是空位缺陷，本文介绍计算空位形成能的模拟方法。
空位形成能计算公式和间隙形成能类似：

主要模拟过程也差不多，间隙原子是额外插入了一个原子，而空位是指在原来的晶格中删除一个原子。
假如体系共有N个原子，删除一个原子后原子总数为（N-1），空位形成能等于删除一个原子后的能量Ef减去（N-1）个原子的能量。

具体in文件如下，代码来源于网络，我进行了修改和简化，对原作者表示感谢。代码已经注释，仅供参考。

#模型基本设置
variable ao equal 3.615
units   metal
dimension