(no subject)

From: Daniel Hornburg (Hleinad_at_gmx.de)
Date: Tue Mar 16 2010 - 07:17:02 CDT

Sorry I think there was a problem with my mail account resulting in not
sending the parts of the mail were I added my script. Hope this time it




I am using NAMD for md and I am relatively new. I have a question concerning
my equilibrium script. When I am watching the resulting trajectory in VMD
(there is no error or something like that produced, all seems to be fine)
after a short movement (probably the minimization see #1) the motion of the
whole system stops for a while. Then it continues with movements at the
point where constrains are decreased (probably, #2).

I think, that the whole loop with the reassignTemp $TEMP is not working
instead the whole system stays at 0.0 kelvin for this loop. I got that
script from another person who is more skilled than me but not available
right now. Anyway I am asking myself what the benefit of this whole loop was
intended to be and I hesitate to change something right on my own without

Hope you might help me. Have great day and Thank you!

Best wishes,



constraints on

consref crd/mwt26ra.pdb

conskfile crd/mwt26ra.pdb

conskcol O

constraintScaling 5


exclude scaled1-4

1-4scaling 1.0

switching on

switchdist 8.0

cutoff 12.0

pairlistdist 14.0


rigidBonds all


fullElectFrequency 1

nonbondedFreq 1


langevin on

langevinTemp 0

langevinDamping 5


cellBasisVector1 104.0 0.0 0.0

cellBasisVector2 0.0 60.0 0.0

cellBasisVector3 0.0 0.0 60.0


PME on

PMEGridSpacing 1.0


timestep 1.00

temperature 0.0


#1 minimation which I can see in vmd/trjfile

minimize 100000


#2 loop wich might not work

for { set TEMP 30 } { $TEMP < 300 } { incr TEMP 20 } {

   reassignTemp $TEMP

   run 10000



#here I think the movment starts again

langevinTemp 295


run 200000


constraintScaling 4

run 200000


constraintScaling 3

run 200000


constraintScaling 2

run 200000


constraintScaling 1

run 200000


constraintScaling 0

run 200000

This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:53:54 CST