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,
6
# or run the file as a script using ./README_setup.sh)
7
# 2) You must equilibrate the system beforehand using "run.in.npt".
8
# This will create the file "system_after_npt.data" which this file reads.
9
# (Note: I have not verified that this equilibration protocol works well.)
11
# ------------------------------- Initialization Section --------------------
13
include "system.in.init"
15
# ------------------------------- Atom Definition Section -------------------
17
# Read the coordinates generated by an earlier NPT simulation
19
read_data "system_after_npt.data"
21
# ------------------------------- Settings Section --------------------------
23
include "system.in.settings"
26
# (The "write_restart" and "read_restart" commands were buggy in 2012,
27
# but they should work also. I prefer "write_data" and "read_data".)
29
# ------------------------------- Settings Section --------------------------
31
include system.in.settings
33
# ------------------------------- Run Section -------------------------------
35
# -- simulation protocol --
39
dump 1 all custom 5000 traj_nvt.lammpstrj id mol type x y z ix iy iz
40
fix fxnvt all nvt temp 300.0 300.0 500.0 tchain 1
42
#thermo_modify flush yes
46
write_data system_after_nvt.data