-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathmsld_prod.inp
160 lines (132 loc) · 4.96 KB
/
msld_prod.inp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
* CHARMM input file for Multi-Site lambda-dynamics
* MSLD references:
* J.L. Knight and C.L. Brooks III. Multi-site lambda-dynamics:
* free energy method for simulating Structure-Activity Relationship
* studies. J. Chem. Theory Comput. 7, 2728-2739 (2011).
* J.L. Knight and C.L. Brooks III. Applying efficient
* implicit non-geometric constraints in free energy simulations.
* J. Comput. Chem. 32, 3423-3432 (2011).
* Stream file for MSLD created by process.py
* Written by Ryan Hayes 2016-06-18
*
!set fnex = 5.5
set temp = 298.15
stream "variablesprod.inp"
bomblev -1
stream "prep/@{SYSNAME}.inp"
! Set nonbonded options
stream ../nbond.str
faster on
! system "mkdir -p res dcd"
! bomlev -6
! domdec gpu only dlb off ndir 1 @nnodes 1
! mini sd nstep 100 nprint 10 step 0.005 !! Do min, Steepest-Decent
! bomlev 0
!---------------------------------------------------------------------------
! Write out structure after minimization
!---------------------------------------------------------------------------
! write psf card form name aftermin.psf
! * after minimization psf
! *
! write coor pdb form name minimized.pdb
! * after minimization pdb
! *
read coor card name prep/minimized.crd
!---------------------------------------------------------------------------
! Setup of dynamics run
!---------------------------------------------------------------------------
if @nreps .eq. 1 then
set myrep = 0
else
!! Initiate BP-REX MSLD
set dph = 0.5
calc sph = 7.0 - @dph*@ncentral
repd fast nrep @nreps EXLM FREQuency 1000 MSPH SPH @sph DPH @dph unit 27
!! Direct output of replicas to their individual files
set myrep = ?myrep
endif
open write form unit 28 name output_@{itt}_@myrep
outu 28
stream "../nbshift/vb.inp"
shake fast bonh param
scalar fbeta set 0.1 sele all end
domdec gpu only dlb off ndir 1 @nnodes 1
calc ittm1 = @itt - 1
!! Production run of 1000 ps / 1 ns
open unit 21 write unform name "dcd/@{SYSNAME}_prod@{ITT}.dcd_@MYREP"
open unit 22 write form name "res/@{SYSNAME}_prod@{ITT}.res_@MYREP"
! open unit 23 read form name "res/@{SYSNAME}_prod@{ITTM1}.res_@MYREP"
open unit 24 write file name "res/@{SYSNAME}_prod@{ITT}.lmd_@MYREP"
if nreps .gt. 1 then
open write form unit 27 name "@{SYSNAME}_prod@{ITT}.rex_@MYREP"
endif
if @itt .eq. 1 then
set resttoken = start
set iurest = -1
else
set resttoken = restart
open unit 23 read form name "res/@{SYSNAME}_prod@{ITTM1}.res_@MYREP"
set iurest = 23
endif
calc seed = @seed + @itt
calc pmass = ?natom * 0.12
dynamics cpt @resttoken -
blade prmc iprs 100 pref 1 prdv 100 -
iseed @seed @seed @seed @seed -
timestep 0.002 - !!timestep interval
nstep @nsteps - !!no. of steps (== 2 ns)
nprint 1000 - !!freq of printout
iprfrq 1000 - !!freq of calc avg/rms energy
nsavc @nsavc - !!freq of writing coordinates
nsavl 10 - !!freq of writing lambda histograms
isvfrq @nsteps - ! 10000 - !!freq of writing out restart files
iunread @iurest - !!read restart file
iuncrd 21 - !!write out coordinates
iunwri 22 - !!write out restart file
iunldm 24 - !!write out lambda histograms and biasing potential
firstt @temp - !!initial temp
finalt @temp - !!final temp
tstruc @temp - !!temp at which starting structure has been equilibrated
tbath @temp - !!temp of langevin bath
pconstant - !! use constant pressure (pressure.doc recommendations)
pmass @pmass - !! pressure "mass" 2% of mass or 400.0 amu
pref 1.0 - !! pressure in atm
pgamma 20.0 - !! pressure collision frequency
hoover -
reft @temp -
tmass 1000 -
ichecw 0 - !!do not scale velocities to final temp (i.e. equilibrate)
ihtfrq 0 - !!freq of heating
ieqfrq 0 - !!freq of scaling/assigning velocity(X)
iasors 1 - !!assign velocities during heating (0 will be scale velocities)
iasvel 1 - !!using gaussian distribution to assign velocities
iscvel 0 -
inbfrq -1 -
ilbfrq 0 -
imgfrq -1 -
ntrfrq 0 - ! 500 -
echeck -1 !!energy tolerance check before crashing
!----------------------------------------------------------------------------
! Write out structure after simulation
!----------------------------------------------------------------------------
NOE
RESET
END
cons harm clear
if @itt .eq. 1 then
write psf card name "@{SYSNAME}_prod@{ITT}.psf_@MYREP"
* after dynamics psf
*
endif
write coor pdb card name "@{SYSNAME}_prod@{ITT}.pdb_@MYREP"
* after dynamics pdb
*
write coor card name "@{SYSNAME}_prod@{ITT}.crd_@MYREP"
* after dynamics crd
*
!---------------------------------------------------------------------------
! Analyze l-dynamics output
!---------------------------------------------------------------------------
open unit 44 read file name "res/@{SYSNAME}_prod@{ITT}.lmd_@MYREP"
traj lamb print ctlo 0.95 cthi 0.99 first 44 nunit 1
stop