From: Norman Geist (norman.geist_at_uni-greifswald.de)
Date: Fri Aug 05 2011 - 05:56:04 CDT
Hi,
where is the error?
Mit freundlichen Grüßen
Norman Geist.
-----Ursprüngliche Nachricht-----
Von: owner-namd-l_at_ks.uiuc.edu [mailto:owner-namd-l_at_ks.uiuc.edu] Im Auftrag
von Narasimhan LOGANATHAN
Gesendet: Freitag, 5. August 2011 12:41
An: NAMD
Betreff: namd-l: water equilibration
Hi all,
I am a newbie to NAMD. Right now i am trying to equilibrate water box.
I am not using the CHARMM forcefield. Hence i defined my own parameter file.
However when i tried to equilibrate or minimize the water i am getting the
following error.
However i have defined the splitPatch hydrogen in the file.
Could someone kindly help me to overcome this.
Herewith, i attach my files for information
## Minimization and Equilibration of Water Box ##
structure water.psf
coordinates water.pdb
set temperature 300 ;# Tcl variable
set outputname water-equi ;# Tcl variable
## starting from scratch ##
temperature $temperature ;# intialize velocities randomly using
Boltzmann Distribution
firsttimestep 0
## Input and ForceField Parameters ##
paratypeCharmm off
parameters clayff-namd-1.inp
if {0} {
cellBasisVector1 31.212 0.0 0.0
cellBasisVector2 0.0 27.000 0.0
cellBasisVector3 0.0 0.0 30.000
cellOrigin 0.0 0.0 0.0
}
exclude none
## 1-4scaling 1.0 ##
cutoff 12.0
switching off
## swtichdist 10.0 ## to be used when switching is on ##
pairlistdist 14.0 ;# cutoff + 2 in general
wrapAll on
## Integrator Parameters ##
stepspercycle 10 ;# redefine pairlists every this many
timesteps
timestep 1.0 ;# 1 fs/step
rigidBonds water ;# to keep the water molecules rigid
nonbondedFreq 1 ;# no. of timesteps used to calculate nonbonded
interactions
fullElectFrequency 1 ;# no. of timesteps used to calculate
electrostatic interactions
splitPatch hydrogen
## Constant Temperature Control ##
langevin on ;# Uses langevin dynamics for controlling T
(or) P
langevinDamping 2 ;# Damping Coefficient in pico seconds
langevinTemp $temperature ;# random noise at this level
langevinHydrogen no ;# no coupling of hydrogens
## Constant Pressure Control is done by the use of the langevin
dynamics (see manual) ##
## Particle Mesh Ewald ##
PME yes
PMEGridSpacing 1.0 ;# Alternatively we can specify the PMEgrid
Manually in case direction ##
#################################################################
### Output Options ###
#################################################################
outputname $outputname
restartfreq 500
dcdFreq 500
xstFreq 200
outputEnergies 100
outputPressure 100
###################################################################
### IMD Settings ###
###################################################################
if {0} {
IMDon on
IMDport 3000 ;# Enter this port number in VMD ##
IMDfreq 1 ;# send every 1 frame ##
IMDwait no ;# wait for VMD to connect before running ##
}
###################################################################
### Execution Script ###
###################################################################
## Minimization ##
if {0} {
minimize 100
reinitvels $temperature
}
run 100000
Regards
Naresh
This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:57:32 CST