-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathremdff_template.namd
132 lines (118 loc) · 3.42 KB
/
remdff_template.namd
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
############################################################################
#cr
#cr (C) Copyright 1995-2009 The Board of Trustees of the
#cr University of Illinois
#cr All Rights Reserved
#cr
############################################################################
############################################################################
# RCS INFORMATION:
#
# $RCSfile: remdff_template.namd,v $
# $Author: ryanmcgreevy $ $Locker: $ $State: Exp $
# $Revision: 1.2 $ $Date: 2017/05/23 15:16:12 $
#
############################################################################
#
# DESCRIPTION:
#
# NAMD template file for MDFF simulations.
#
# Authors: Leonardo Trabuco <[email protected]>
# Elizabeth Villa <[email protected]>
# Kwok-Yan Chan <[email protected]>
#
############################################################################
#firsttimestep 0
# Force-Field Parameters
exclude scaled1-4
1-4scaling $SCALING_1_4
switching on
dielectric $DIEL
margin $MARGIN
if {$GBISON} {
gbis on
alphacutoff 14.
switchdist 15.
cutoff 16.
pairlistdist 17.
ionconcentration 0.1
solventDielectric 80.0
sasa on
} else {
cutoff 10.
switchdist 9.
pairlistdist 11.
}
# Integrator Parameters
timestep 1.0 ;# 1fs/step
rigidBonds water ;# TIP3P
nonbondedFreq 2
fullElectFrequency 4
stepspercycle 20
# Constant Temperature Control
if {$ITEMP == $FTEMP} {
langevin on ;# do langevin dynamics
langevinDamping 5 ;# damping coefficient (gamma) of 5/ps
langevinTemp $FTEMP
langevinHydrogen off ;# don't couple langevin bath to hydrogens
} else {
reassignFreq 1000
if {$FTEMP > $ITEMP} {
reassignIncr 10
} else {
reassignIncr -10
}
reassignTemp $ITEMP
reassignHold $FTEMP
}
# Output
#outputName $OUTPUTNAME
#restartfreq 1000
#dcdfreq 1000
#outputEnergies 100
#outputPressure 100
# gridForces (to apply forces from 3D potential)
if {$GRIDON} {
mgridForce on
for {set i 0} {$i < [llength $GRIDFILE]} {incr i} {
mgridForceFile $i [lindex $GRIDPDB $i]
mgridForceCol $i O
mgridForceChargeCol $i B
mgridForcePotFile $i [lindex $GRIDFILE $i]
mgridForceScale $i [lindex $GSCALE $i] [lindex $GSCALE $i] [lindex $GSCALE $i]
}
}
# extraBonds (for secondary structure restraints)
if {$EXTRAB != 0} {
extraBonds yes
foreach fil $EXTRAB {
extraBondsFile $fil
}
}
# Positional restraints
if {$CONSPDB != 0} {
Constraints yes
ConsExp 2
ConsRef $CONSPDB
ConsKFile $CONSPDB
ConskCol $CONSCOL
}
# Fixed atoms
if {$FIXPDB != 0} {
fixedAtoms yes
fixedAtomsFile $FIXPDB
fixedAtomsCol $FIXCOL
}
#############################################################
## EXECUTION SCRIPT ##
#############################################################
#if {$MS != 0} {
# minimize $MS
# reinitvels $ITEMP
#}
#if {$ITEMP == $FTEMP} {
# run $TS
#} else {
# run [expr abs($FTEMP-$ITEMP)*100]
#}