对SPPARKS了解不算深,稍微折腾过一点点。
SPPARKS是支持off-lattice 的 standard Metropolis MC atomic relaxiation的,但不清楚Fe-Co的这种随机交换有没有包含在标准算法里面。如果没有的话,强行进行格点交换的barrier估计会大于5 eV,仅仅靠随机移动很难越过。
如果引入空位作为扩散载体的话barrier应该会降低不少,但估计效率还是堪忧。
@yang元祐 最近看LAMMPS手册,发现LAMMPS也支持MC模拟,和你的需求看起来挺匹配: