lammps蒙特卡洛MC模拟案例

大家好，我是小马老师。
本文介绍lammps和蒙特卡洛(MC)相结合的模拟方法。
MC是一种随机模拟方法，可以模拟一定区域内的随机原子交换过程。在lammps中，和MC相关的模拟命令有三个：fix gcmc、fix atom/swap、fix sgcmc。
本文主要介绍fix sgcmc命令。sgcmc命令可应用到合金的模拟中，如短程有序模拟。
fix sgcmc语法格式：

fix ID group-ID sgcmc every_nsteps swap_fraction temperature deltamu ...

• ID, group-ID：fix命令的id和原子组的id sgcmc：该命令的关键词 every_nsteps：每隔every_nsteps步执行一次MC循环
  swap_fraction：每一个MC循环周期内执行的交换步数的比例
  temperature：MC循环中的温度，通常这个温度与MD模拟的温度相同 deltamu：化学势差

在一个lammps模拟中，mc和md可同时存在，间隔运行，所以设置fix sgcmc命令的同时，可以同时设置系综（nve、nvt或者npt）。
如：

fix mc all sgcmc 50