NAMD minimisation

From: Hua Wong (wong_at_ebgm.jussieu.fr)
Date: Thu Nov 27 2003 - 07:29:27 CST

Hello. I have trouble with setting up a minimisation under NAMD.
Here is the NAMD script I use....

# NAMD configuration file for PROT

# molecular system
structure PROT.psf

# force field
paratypecharmm on
parameters ../par_all22_prot.inp
exclude scaled1-4
1-4scaling 1.0

# approximations
switching on
switchdist 8
cutoff 12
pairlistdist 13.5
margin 0
stepspercycle 20

#integrator
numsteps 10000
timestep 1

#output
outputenergies 1
outputtiming 100
binaryoutput no
imdon yes
imdfreq 1
imdport 3111

# molecular system
coordinates PROT.pdb

#output
outputname PROT_out
#dcdfreq 1000

# constraints/Fixed
fixedAtoms on
fixedAtomsFile PROT.pdb
fixedAtomsCol B

#protocol
temperature 400
langevin no

rescaleFreq 50
rescaleTemp 400

#script

#minimize 400

run 20000

and the results I obtain

ETITLE: TS BOND ANGLE DIHED
IMPRP ELECT VDW BOUNDARY MISC
KINETIC TOTAL TEMP TOTAL2 TOTAL3
TEMPAVG

ENERGY: 0 2685892.8098 -99999999.9999 -99999999.9999
-99999999.9999 -13233.6249 99999999.9999 0.0000
0.0000 -99999999.9999 -99999999.9999 -99999999.9999 -99999999.9999
-99999999.9999 -99999999.9999

ENERGY: 1 -99999999.9999 -99999999.9999 -99999999.9999
-99999999.9999 -99999999.9999 -99999999.9999 0.0000
0.0000 -99999999.9999 -99999999.9999 -99999999.9999 -99999999.9999
-99999999.9999 -99999999.9999

ENERGY: 2 -99999999.9999 -99999999.9999 -99999999.9999
-99999999.9999 -99999999.9999 -99999999.9999 0.0000
0.0000 -99999999.9999 -99999999.9999 -99999999.9999 -99999999.9999
-99999999.9999 -99999999.9999

I am a bit puzzled because someone else made a script that is identical
and it seems to work.....

Any hint?

Thanks

Wong Hua

This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:37:10 CST