Dear Grzegorz:

Under no conditions should any of the tip5p geometry change (for the standard 
tip5p model). 
If you find that this is happening, then that is certainly an error. You can 
check if you like by analyzing 
your trajectory. However, flexible bonds will allow the distance from the 
arginine N to the arginine
H to vary, which my allow a closer approach of the arginine H to the tip5p 
dummy site.

Did you verify that a water box (no protein) simulates without error?

Did you post a redmine issue with .mdp , .gro , and .top files?

Chris.

-- original message --

Dear Chris,
Authors answered me very quickly. They did not have such a problem, but 
I still don't know the details of their input. They used gromacs-3.3, so 
I decided to give the old one a try. I did some tests with 3.3.4. 
Although the same problem occurred during steep minimization, some 
interesting things popped out. When I tried to grompp the system with 
plain tip5p for cg minimization, it failed:
"ERROR: can not do Conjugate Gradients with constraints (8484)", even 
though I did not set any constraints. The error is the same for tip4p 
unless you use flexible model, which tip5p does not have - the water had 
to be constrained then. I guess that treatment of virtual sites in 
gromacs-3.3 has something to do with this. I noticed, that constraints 
make simulations with tip5p more stable. It should not happen, that the 
LP virtual atoms are pulled further then the defined 0.7 from the 
oxygen, right? I will keep you updated.
Best Regards,
Grzegorz

--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Reply via email to