本文分享一个计算界面能的lammps案例。 主要的原理是先计算整个结构的总势能,然后从中间切开,形成两个界面,把两部分完全分离后计算新的势能,带入公式得到界面能。 本文介绍另外一种方法:在周期性边界ppp下计算体系的势能Ebulk,然后修改边界条件为mpp,计算新的势能Ecut。 边界改为m边界后,x方向的左右两侧原子不再产...