氧化石墨烯+水体系平衡模拟

#模型基本设置

boundary  p p p

units    real

atom_style  full

#势函数参数设置

bond_style  harmonic

angle_style  harmonic

dihedral_style  opls

improper_style  harmonic

#tip4p

pair_style  lj/cut/tip4p/long 1 2 1 1 0.1546 12.0

kspace_style  pppm/tip4p 1.0e-4

pair_modify   shift yes mix arithmetic

special_bonds  lj/coul 0.0 0.0 0.5 angle yes dihedral yes

#读取模型文件

read_data  data.lammps

#按分子方式进行分组

group    H2O type 1

group    H2O include molecule

group    GO type 3

group    GO include molecule

#系综设置

fix    mynve all nve

fix    myber all temp/berendsen 300 300 100

#轨迹文件输出设置

dump            mydmp all atom 1000 dump.lammpstrj

#固定水分子键角

fix    s1 H2O shake 1.0e-4 200 0 b 1 a 1

#模拟步长

timestep   1.0

thermo   1000

#启动运算

run    50000