35 lines
901 B
Plaintext
35 lines
901 B
Plaintext
# 3d Lennard-Jones melt - MDI driver script
|
|
|
|
units lj
|
|
atom_style atomic
|
|
atom_modify map array
|
|
|
|
read_data data.snapshot
|
|
|
|
velocity all create 1.44 87287 loop geom
|
|
|
|
pair_style lj/cut 2.5
|
|
pair_coeff 1 1 1.0 1.0 2.5
|
|
|
|
neighbor 0.3 bin
|
|
neigh_modify delay 0 every 1 check yes
|
|
|
|
fix 1 all nve
|
|
|
|
fix 2 all mdi/qm add no every 100 virial yes
|
|
|
|
compute 1 all pressure NULL virial
|
|
variable evirial equal (f_2[1]+f_2[2]+f_2[3])/3
|
|
|
|
thermo_style custom step temp pe c_1 c_1[1] c_1[2] c_1[3]
|
|
|
|
thermo 100
|
|
|
|
dump 1 all custom 100 dump.snapshot.driver &
|
|
id type x y z f_2[1] f_2[2] f_2[3]
|
|
dump_modify 1 sort id
|
|
|
|
run 300 pre no post no every 100 &
|
|
"print 'QM eng = $(f_2/atoms)'" &
|
|
"print 'QM virial = $(v_evirial) $(f_2[1]) $(f_2[2]) $(f_2[3])'"
|