lammps案例：SiC结合能的模拟

大家好，我是小马老师。

本文介绍lammps模拟结合能的案例代码。

分散的原子之所以能构成稳定的晶体，是因为原子之间存在一种将原子或分子结合在一起的力。

原子在结合成稳定晶体的过程中，所释放出来的能量，称为晶体的结合能(binding energy)，也称为内聚能（cohesive energy）。

纯金属的内聚能可参考：
Fe内聚能的模拟代码

多种原子构成的晶体可以按照同样的原理计算出结合能。

以sic为例，对应的in文件为：

units        metal
boundary     p p p 
atom_style   full
read_data   sic.data   
pair_style    tersoff
pair_coeff    *