<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
</head>
<body class='hmmessage'>
Hi,<br><br>The choice of charge is significant even with PME, because of the way<br>Gromacs makes the neighborlists.<br><br>As I mailed before, use the -nochargegrp option of pdb2gmx with the Charmm27 ff.<br><br>Berk<br><br>&gt; Date: Fri, 17 Sep 2010 11:46:36 +0200<br>&gt; From: qiaobf@gmail.com<br>&gt; To: gmx-users@gromacs.org<br>&gt; Subject: Re: [gmx-users] ARG Charmm gmx 4.5.1<br>&gt; <br>&gt; <br>&gt; You  should change the charge group number. Every charge group include<br>&gt; several, like 3-4, atoms. Actually, the definition of charge group is<br>&gt; not significant if you use PME for the coulomb interaction.<br>&gt; <br>&gt; <br>&gt; nahren manuel wrote:<br>&gt; &gt; Dear Gromacs Users,<br>&gt; &gt;<br>&gt; &gt; I am using plain cutoff for my 12-mer protein.<br>&gt; &gt; The grompp reports ARG to have a big charge group. this was also highlighted in the following mail<br>&gt; &gt; http://www.mail-archive.com/gmx-users@gromacs.org/msg32098.html<br>&gt; &gt;<br>&gt; &gt; I was just think if changing the charges on these atoms would help,<br>&gt; &gt; from<br>&gt; &gt; 13    CT2      1    ARG     CD     4        0.2     12.011   ; qtot 1.2<br>&gt; &gt; 14     HA      1    ARG    HD1     4       0.09      1.008   ; qtot 1.29<br>&gt; &gt; 15     HA      1    ARG    HD2     4       0.09      1.008   ; qtot 1.38<br>&gt; &gt; 16    NC2      1    ARG     NE     4       -0.7     14.007   ; qtot 0.68<br>&gt; &gt; 17     HC      1    ARG     HE     4       0.44      1.008   ; qtot 1.12<br>&gt; &gt; 18      C      1    ARG     CZ     4       0.64     12.011   ; qtot 1.76<br>&gt; &gt; 19    NC2      1    ARG    NH1     4       -0.8     14.007   ; qtot 0.96<br>&gt; &gt; 20     HC      1    ARG   HH11     4       0.46      1.008   ; qtot 1.42<br>&gt; &gt; 21     HC      1    ARG   HH12     4       0.46      1.008   ; qtot 1.88<br>&gt; &gt; 22    NC2      1    ARG    NH2     4       -0.8     14.007   ; qtot 1.08<br>&gt; &gt; 23     HC      1    ARG   HH21     4       0.46      1.008   ; qtot 1.54<br>&gt; &gt; 24     HC      1    ARG   HH22     4       0.46      1.008   ; qtot 2<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt; to<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt;         <br>&gt; &gt;         <br>&gt; &gt;         <br>&gt; &gt;         <br>&gt; &gt;         <br>&gt; &gt;                 <br>&gt; &gt;    CD    0.18<br>&gt; &gt;   HD1    0.06<br>&gt; &gt;   HD2    0.06<br>&gt; &gt;    NE    -0.7<br>&gt; &gt;    HE    0.4<br>&gt; &gt;    CZ    0.6<br>&gt; &gt;   NH1    -0.8<br>&gt; &gt;  HH11    0.5<br>&gt; &gt;  HH12    0.5<br>&gt; &gt;   NH2    -0.8<br>&gt; &gt;  HH21    0.5<br>&gt; &gt;  HH22    0.5<br>&gt; &gt;<br>&gt; &gt; The above transformation of charges seems reasonable.<br>&gt; &gt;<br>&gt; &gt; Would like to know if this is okay...<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt; Best,<br>&gt; &gt; nahren<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt;<br>&gt; &gt;       <br>&gt; &gt;   <br>&gt; <br>&gt; -- <br>&gt; gmx-users mailing list    gmx-users@gromacs.org<br>&gt; http://lists.gromacs.org/mailman/listinfo/gmx-users<br>&gt; Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting!<br>&gt; Please don't post (un)subscribe requests to the list. Use the <br>&gt; www interface or send it to gmx-users-request@gromacs.org.<br>&gt; Can't post? Read http://www.gromacs.org/Support/Mailing_Lists<br>                                               </body>
</html>