[gmx-users] pressure scaling
Berk Hess
gmx3 at hotmail.com
Tue Aug 3 11:11:49 CEST 2004
>
>Before using this values of rvdw, rcoulomb, rlist and nstlist I had
>rcoulomb=rlist=nstlist=1.3 and only SR energies, and when I turned rlist to
>1.1 I get SR and LR energies.
OK.
But this should not have a large effect on the simulation.
Having all cut-off's at 1.3 is more accurate, but is
computationally more expensive.
Something else:
Are you using constraints?
If you use constraints you should be able to use a time step
of at least 0.002 and maybe even up to 0.005.
Many other things could be wrong.
The most probable cause is always an not well equilibrated
starting structure.
But there could also be an error in your topology.
What is your tau_t?
The coupling time for temperature bath should be a least
a factor of two smaller than tau_p.
Berk.
_________________________________________________________________
Hotmail en Messenger on the move
http://www.msn.nl/communicatie/smsdiensten/hotmailsmsv2/
More information about the gromacs.org_gmx-users
mailing list