[gmx-users] How to overcome Fatal error of bonded interactions Gromacs 4.0.3

Claus Valka lastexile7gr at yahoo.de
Wed Mar 4 12:12:19 CET 2009


the system I'm referring at is an already equilibrated one. Before, I had already made steepest descents. 

You are correct that I have non-bonded interactions at distances that are more than cut-off+table ext distance. I used to put table ext very bin in contrast to my box. My system box was something less than 5nm before and table-ext was 5nm at the time. Now that I simulate a much bigger system than before (in one dimension 6 times bigger), I realized that I had to increase also this parameter.

Since this parameter affects only the pairs interactions I do not see the reason why my simulations are affected by it, even if it is already big enough (5nm in a 30nm box). My pairs interactions are (should) not (be) more that 1nm far away.

I had extensively tested this parameter in gromacs 3, and I saw that it affects the pairs interactions only. So I decided to increase it very much and not bother with it.

Is there any difference in newest versions? Does it affect anything more?

Thank you,

--- Berk Hess <gmx3 at hotmail.com> schrieb am Mi, 4.3.2009:
Von: Berk Hess <gmx3 at hotmail.com>
Betreff: RE: [gmx-users] How to overcome Fatal error of bonded interactions Gromacs 4.0.3
An: "Discussion list for GROMACS users" <gmx-users at gromacs.org>
Datum: Mittwoch, 4. März 2009, 9:23

#yiv174728022 #yiv779874946 .hmmessage P
#yiv174728022 #yiv779874946 {



Your solution seems to indicate that your system is unstable.
My guess is that your system was not well equilibrated and at some point
some non-bonded interactions were at a distance of more than cut-off+tab_ext.
That would make your system explode. In that case increasing tab_ext to 2 nm
should also suffice.


Date: Tue, 3 Mar 2009 08:33:52 -0800
From: lastexile7gr at yahoo.de
To: gmx-users at gromacs.org
Subject: [gmx-users] How to overcome Fatal error of bonded interactions	Gromacs 4.0.3


this has to do with my personal experience of using Gromacs version 4.0.3. So this is not an
 official response.

The story:

While I was trying to run an NVT simulation in a monoclinic box I encountered the following problem:

Fatal error:
25 of the 336080 bonded interactions could not be calculated because some atoms involved moved further apart than the multi-body cut-off distance (0.895 nm) or the two-body cut-off distance (0.895 nm), see option -rdd, for pairs and tabulated bonds also see option -ddcheck

and the:
Not all bonded interactions have been properly assigned to the domain decomposition cells

Trying to put aside the fact that my system had any bug I found in the net the

In order to overcome this error I increased the value of
table-extension          = 30
in my mdp file in such a way, as to be bigger from my biggest dimension I had in my system. The ran proceeded (until now that I'm writing this e-mail) without any problem.

This value has to do with the cut-off of pairs1-4 LJ interactions and is calculated in nm.


See all the ways you can stay connected to friends and family 
gmx-users mailing list    gmx-users at gromacs.org
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the 
www interface or send it to gmx-users-request at gromacs.org.
Can't post? Read http://www.gromacs.org/mailing_lists/users.php

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20090304/b9c5e88a/attachment.html>

More information about the gromacs.org_gmx-users mailing list