大家好,我是小马老师。
在昨天的推文中,介绍了计算Cu间隙原子形成能的lammps模拟方法。
与间隙原子相对应的另一种缺陷是空位缺陷,本文介绍计算空位形成能的模拟方法。
空位形成能计算公式和间隙形成能类似:
主要模拟过程也差不多,间隙原子是额外插入了一个原子,而空位是指在原来的晶格中删除一个原子。
假如体系共有N个原子,删除一个原子后原子总数为(N-1),空位形成能等于删除一个原子后的能量Ef减去(N-1)个原子的能量。
具体in文件如下,代码来源于网络,我进行了修改和简化,对原作者表示感谢。代码已经注释,仅供参考。
#模型基本设置
variable ao equal 3.615
units metal
dimension