3
# 1) You must use moltemplate.sh to create 3 files:
4
# system.data system.in.init system.in.settings
5
# (Follow the instructions in README_setup.sh, or run it using ./README_sh.)
6
# 2) You must equilibrate the system beforehand using "run.in.npt".
7
# This will create the file "system_after_npt.data" which this file reads.
8
# (Note: I have not verified that this equilibration protocol works well.)
10
# ------------------------------- Initialization Section --------------------
12
include "system.in.init"
14
# ------------------------------- Atom Definition Section -------------------
16
# Read the coordinates generated by an earlier NPT simulation
18
read_data "system_after_npt.data"
20
# OPLSAA atom charges are stored in a separate file.
23
include "system.in.charges"
25
# ------------------------------- Settings Section --------------------------
27
include "system.in.settings"
30
# (The "write_restart" and "read_restart" commands were buggy in 2012,
31
# but they should work also. I prefer "write_data" and "read_data".)
33
# ------------------------------- Settings Section --------------------------
35
include system.in.settings
37
# ------------------------------- Run Section -------------------------------
39
# -- simulation protocol --
43
dump 1 all custom 5000 traj_nvt.lammpstrj id mol type x y z ix iy iz
44
fix fxnvt all nvt temp 300.0 300.0 500.0 tchain 1
46
#thermo_modify flush yes
50
write_restart system_after_nvt.data