[gmx-users] Charge system simulation problem

Justin Lemkul jalemkul at vt.edu
Mon Nov 26 15:58:23 CET 2018



On 11/26/18 12:01 AM, Karpurmanjari Kakati wrote:
> Dear All,
>
> My System contains, while doing NVT of the same, I got the following Note
>
>
>
>
>
>
> *NOTE 1 [file sa.top, line 80]:  System has non-zero total charge:
> 12.936997  Total charge should normally be an integer. See
> http://www.gromacs.org/Documentation/Floating_Point_Arithmetic
> <http://www.gromacs.org/Documentation/Floating_Point_Arithmetic>  for
> discussion on how close it should be to an integer.*
>
>
> This is possibly because of  the partial charges in the molecules in the
> ITP fie.

Something isn't parametrized correctly because you need to have an 
integer charge. Your total charge is not physical.

-Justin

> My doubt is can't I have a system with non zero and a decimal charge value.
>
>
> Also I got the following warning:
>
>
>
>
>
>
> *WARNING 1 [file sa.top, line 80]:  You are using Ewald electrostatics in a
> system with net charge. This can  lead to severe artifacts, such as ions
> moving into regions with low  dielectric, due to the uniform background
> charge. We suggest to  neutralize your system with counter ions, possibly
> in combination with a  physiological salt concentration.*
>
> Can Someone help me with this?
>
> Also this :
>
>
>
>
>
>
>
> *Velocities were taken from a Maxwell distribution at 300 KRemoving all
> charge groups because cutoff-scheme=VerletAnalysing residue names:There
> are:     1      Other residuesThere are: 100000      Water residues*
>
>
> Removing all charge group because the cut off scheme is verlet.
>
> If That happens then will it show the behavior I want to see. Or shall I
> put any other scheme on that ?
>
> Please help

-- 
==================================================

Justin A. Lemkul, Ph.D.
Assistant Professor
Office: 301 Fralin Hall
Lab: 303 Engel Hall

Virginia Tech Department of Biochemistry
340 West Campus Dr.
Blacksburg, VA 24061

jalemkul at vt.edu | (540) 231-3129
http://www.thelemkullab.com

==================================================



More information about the gromacs.org_gmx-users mailing list