求助：如何用vasp估算体系的键能

XYL123

各位老师好，学生想估算一下石墨烯、氮化硼以及联苯烯内，C-C键、B-N键的键能。
采用的方法是①孤立原子的能量减去体系的总能，得到的能量除以键的个数，得到每个键的平均键能。
②用断键后的体系总能减去初始体系的总能，得到的能量就是断裂的几个键的总能，在平均一下。比如我采用石墨烯2x2的超胞，在MS中把晶格常数扩大，调整原子之间的距离，然后进行优化。（PS：这样操作，放入vasp里优化后结构会跑散，与我构建的模型完全不同。只能进行自洽计算算总能了）
下面是我的方法二的设想建模