[gmx-users] System has non-zero total charge: 1.999878e+00
zakiya
zakiya at tju.edu.cn
Sat Sep 12 11:04:53 CEST 2009
Hi Gromacs users,
I have some problems about the total charge of system using Gromacs.
I wanted to perform a protein denaturation simulation induced by urea. At first, I added urea molecules into the box including one CI2 protein. Then I performed the energy minimization step. However, there’s a note that “NOTE 1 [file ci2.top, line 3998]: System has non-zero total charge: 1.999878e+00 “.
My question is that why the total charge of system is not an integer. The whole charge of CI2 protein is 2.000e+00. I got the ITP file of urea molecule from the “share” file in the Gromacs library directory. The total charge is zero. So this phenomena puzzled me.
The [atoms] field in the ITP file of urea molecule is as follows:
[ atoms ]
; nr type resnr residu atom cgnr charge mass
#ifdef Boek
1 C 1 UREA C1 1 0.38 12.0110
2 O 1 UREA O2 1 -0.38 15.9994
3 NT 1 UREA N3 2 -0.83 14.0067
4 H 1 UREA H4 2 0.415 1.0080
5 H 1 UREA H5 2 0.415 1.0080
6 NT 1 UREA N6 3 -0.83 14.0067
7 H 1 UREA H7 3 0.415 1.0080
8 H 1 UREA H8 3 0.415 1.0080
#else
1 C 1 UREA C1 1 0.683 12.0110
2 O 1 UREA O2 1 -0.683 15.9994
3 NT 1 UREA N3 2 -0.622 14.0067
4 H 1 UREA H4 2 0.346 1.0080
5 H 1 UREA H5 2 0.276 1.0080
6 NT 1 UREA N6 3 -0.622 14.0067
7 H 1 UREA H7 3 0.346 1.0080
8 H 1 UREA H8 3 0.276 1.0080
#endif
Any help/suggestions will be greatly appreciated!
Thank you!!
Na Zhang
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20090912/9d5b10e9/attachment.html>
More information about the gromacs.org_gmx-users
mailing list