On 10/15/12 11:45 AM, benjfitz wrote:
As a further check I ran a neat solvent box and the simulation ends abruptly
at 1.07 us.
These results suggest that your physical model is unstable. A linear model of a
triatomic species is not stable. A more appropriate approach is to use virtual
sites as in the following tutorial:
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin/gmx-tutorials/vsites/index.html
Also note that your cutoffs should not be dictated by your box size, rather the
opposite should be true. The cutoffs are an essential part of the force field
model and changing them haphazardly can easily invalidate that force field model.
-Justin
--
========================================
Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin
========================================
--
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
