上一篇文章中介绍了npt系综的设置方法,本文给出一个lammps模拟氢气压缩的案例,熟悉npt系综的设置方法。
模拟温度300k,压强100GPa,在该条件下,氢气应该为固态。
在MS中建立H2模型,共包含1000个H原子,导出H.car,使用msi2lmp软件转化为H.data文件。
H原子之间作用力使用LJ势描述,H-H键作用力使用harmonic势,具体参数在H.data文件中。
模拟源代码为:
#模拟参数设置
variable t equal 300
variable p equal 1000000
units metal
boundary p p p
atom_style full
neighbor 0.2 bin
timestep 0.0001
neigh_modify ev