[gmx-users] Energy minimization

[Markus O Kaukonen]

Dear All,

My understanding is that Gromacs with constraints does not work with
conjugate gradient minimisation.

To my opinion:
on the other hand one needs constrain water molecules (SPCwater, oplsFF)
(ie bond angle and bond length) in order to
get correct charging energies (and solvation energies)
of ions in a water box. I made this conclusion after some
testing and comparing to literature values
(darden-JCP-v109-pp10921-10935-y1998).

So with systems involving water solvent one should use
constrained water, steepest descent and simulated annealing.
To my experiense the error increases with increasing system size when not
constraining water molecules.

Please correct me if I'm wrong. My test calculation were only qualitive
with energy only (not free energy).

Cheers, Markus


-- 
-----------www= http://www.iki.fi/markus.kaukonen---------------
Markus.Kaukonen at helsinki.fi   \ HOME->   Viinirinne 3 F 12
Helsinki Bioenergetics Group   \_______  02630 Espoo, Finland
Institute of Biotechnology             \ tel:+358-(0)9-5127122/h
Biocenter 3, P.O.BOX 65 (Viikinkaari 1) \____
FIN-00014 University of Helsinki              \______ H+ ______
tel: +358-(0)9-19159752/o, fax: +358-(0)9-19159920             \
--- Rikos ei kannata, eika maatalous Suomessa. (Paimio 1998) ---