andreab at physiol.unizh.ch
Fri Jan 21 17:21:56 CET 2005
Hallo to all the Gromacs users
I would like to ask some general questions about the meaning of the Fmax
value that you get after the energy minimization procedure. The
question is simple I want to simulate the beaviour of a protein (about
390 aa) in water solution; using the alghoritms like steepest descent
and l-bfgs I am not able to get an Fmax below 2000 no constraints were
settled (checked the tpr file and no CONS were found), instead I'm not
able to use cg at all because the program grompp complains about this
error: ERROR can not do conjugate gradients with constraints (29796).
I have such a kind of problem only when I introduce water because when I
am working with the protein alone I can go down with Fmax up to 50.
My questions are:
i) What is a reasonable value for Fmax to set for a MD simulation in
water? Is it there a general good Fmax rang of values or its changing
accordingly to the simulated system?
ii) When I convert the pdb file with pdb2gmx if I use option 0 I
always get extra charge to balance, which instead I don't get if I'm
using option 1 (G43b1 GROMOS96 43b1 Vacuum Forcefield). Is it ok to use
option 1 as forcefield, then solvate the previous output with water and
go for an energy minimization run?
iii) any suggestions on how to write a working em.mdp script with a
reasonable Fmax value?
I already posted some questions about this em problem and I got some
good suggestions but I was not able to solve it.
Thank you for your attention and great help
More information about the gromacs.org_gmx-users