---------- Forwarded message ----------
Date: Wed, 5 Mar 2003 11:54:19 +0100 (MET)
From: seraphimo.gmx.de
To: amber.heimdal.compchem.ucsf.edu
Subject: nasty npt md problem
hi everybody,
i have some deep troubles with my md simulations and wonder if somebody
could help.
the system i prepared is a polymer under pbc plus a tip3 water molecule in
order
to study diffusion behaviour. anyway, while i'm doing nvt md everything
works fine.
when i switch to npt md then there is a raise in density at a certain point
followed by
a drop and the simulation stops for ewald bomb reasons. obviously there is
some vlimit
exception. i used different cutoffs, different timesteps, shake ntc=3,
annealed for
more than two nanoseconds, minimized several times and always run into this
problem.
i attache my mdin file and would be glad if somebody could give me a hint.
cheers, mike.
&ctrl
imin=0,
ntx=5,
ntb=2
irest=1,
iwrap=1,
ntwx=100,
ntave=100,
cut=7.0,
ntc=3,
jfastw=0,
nstlim=200000,
dt=0.001,
temp0=300,
tempi=300,
ntt=1,
pres0=1.0,
taup=1.0,
comp=28.0,
&end
--
+++ GMX - Mail, Messaging & more http://www.gmx.net +++
Bitte lächeln! Fotogalerie online mit GMX ohne eigene Homepage!
Received on Wed Mar 05 2003 - 18:53:01 PST
