TIP4P水球蒸发

#模型基本设置

units real

atom_style full

boundary p p p

timestep   1.0

neighbor 2.0 bin

neigh_modify every 1 delay 0 check yes

#势函数类型设置

bond_style  harmonic

angle_style  harmonic

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

#读取模型文件

read_data H2O.data

#势函数参数设置

pair_coeff 1 1 0.1852 3.1589 #O-O

pair_coeff 2 2 0 0 #H-H

bond_coeff 1 0 0.9572 #H-O

angle_coeff 1 0 104.52 #H-O-H

fix 01 all shake 0.0001 50 0 b 1 a 1

#输出

thermo     1000

dump     1 all xyz 1000 dump.xyz

dump_modify 1 element O H

#蒸发过程

velocity   all create 273 4928459 mom yes rot yes dist gaussian 

fix         1 all npt temp 273 473 100 iso 1 1 1000

run    1000

unfix 1

fix 1 all nvt temp 473 473 100

run 100000