Justin A. Lemkul wrote:
Guido Polles wrote:
Hi everybody,
Sorry about this, but I'm pretty new to gromacs and I couldn't find
any clue about my problem in the archives or in the manual.
I'm using gromacs 4.5.1, and my system is a rather big system (about
12x12x12) that includes a proteon, an atp molecule, water and some
ions.
I am using ffamber99sb force field. Now, I post the log for my mdrun.
Initializing Domain Decomposition on 2 nodes
Dynamic load balancing: auto
Will sort the charge groups at every domain (re)decomposition
Initial maximum inter charge-group distances:
two-body bonded interactions: 14.813 nm, LJ-14, atoms 25682 25685
multi-body bonded interactions: 14.813 nm, Angle, atoms 25685 25686
Minimum cell size due to bonded interactions: 16.294 nm
Using 0 separate PME nodes
Scaling the initial minimum size with 1/0.8 (option -dds) = 1.25
Optimizing the DD grid for 2 cells with a minimum initial size of
20.368 nm
The maximum allowed number of cells is: X 0 Y 0 Z 0
-------------------------------------------------------
Program mdrun, VERSION 4.5.1
Source code file: domdec.c, line: 6428
Fatal error:
There is no domain decomposition for 2 nodes that is compatible with
the given box and a minimum cell size of 20.368 nm
Change the number of nodes or mdrun option -rdd or -dds
Look in the log file for details on the domain decomposition
For more information and tips for troubleshooting, please check the
GROMACS
website at http://www.gromacs.org/Documentation/Errors
Of course the main point here is this
Initial maximum inter charge-group distances:
two-body bonded interactions: 14.813 nm, LJ-14, atoms 25682 25685
multi-body bonded interactions: 14.813 nm, Angle, atoms 25685 25686
but they're CL ions! Why should they have any bonded interaction or
angle interaction??
They shouldn't. Check your [bonds] and [angles] directives for these
atom numbers. Apparently mdrun is finding bonds to these ions in your
topology. Were these ions added with genion? Present in the initial
structure processed by pdb2gmx? Manually added in some location?
Actually, I made some tests, and:
1- if I solvate the protein+atp+ions (I want the ions to be in
specific positions) system I get that 14nm bond error when I run
mdrun, but grompp let it go.
OK, so you want ions in specific positions. Did you put them there? If
so, was it before or after you generated the topology (presumably with
pdb2gmx)?
Also relevant here, but I forgot to ask: are you using distance restraints to
keep these ions in place? If so, that's your problem. If the restraints can't
be placed in the same DD cell, you have to use particle decomposition (mdrun
-pd), which is somewhat slower.
-Justin
--
========================================
Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
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 [email protected]
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 [email protected].
Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
