<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Verdana
}
--></style>
</head>
<body class='hmmessage'>
Hi,<br><br>This is not a decision, but a limitation of the grompp code.<br>I have fixed this in the git release.<br><br>Simply editing the charge group numbers by hand to make them all different<br>should work.<br><br>Berk<br><br>&gt; Date: Mon, 3 May 2010 10:22:54 -0400<br>&gt; From: jalemkul@vt.edu<br>&gt; To: gmx-users@gromacs.org<br>&gt; Subject: Re: [gmx-users] cgnr in top file<br>&gt; <br>&gt; <br>&gt; <br>&gt; Erik Marklund wrote:<br>&gt; &gt; True. I guess gromacs groups together the terminal hydrogens, overriding <br>&gt; &gt; the original force field parameters. Can't say why that decision was <br>&gt; &gt; made, but it will influence your simulation very little.<br>&gt; <br>&gt; I think there was a bug report some time ago involving virtual sites on the <br>&gt; termini of all-atom force fields causing a crash when the terminus crossed a <br>&gt; periodic boundary.  Could be related to that, although I can't remember the <br>&gt; exact progress behind the fix.<br>&gt; <br>&gt; -Justin<br>&gt; <br>&gt; &gt; <br>&gt; &gt; Erik<br>&gt; &gt; <br>&gt; &gt; leila karami skrev:<br>&gt; &gt;&gt; Dear Erik<br>&gt; &gt;&gt;  <br>&gt; &gt;&gt; you said Amber03 was not parameterized with charge groups and last <br>&gt; &gt;&gt; column in your rtp entry holds the chargegroup numbers. but<br>&gt; &gt;&gt; cgnr column in my top file does not looks just like it. in case of H1, <br>&gt; &gt;&gt; H2 and H3.<br>&gt; &gt;&gt;  <br>&gt; &gt;&gt; /     H1    amber99_17   0.16420     2<br>&gt; &gt;&gt; //     H2    amber99_17   0.16420     3<br>&gt; &gt;&gt; //     H3    amber99_17   0.16420     4<br>&gt; &gt;&gt; //    // /. that's ok. but why<br>&gt; &gt;&gt;  <br>&gt; &gt;&gt; /and my top file is as follow :<br>&gt; &gt;&gt; //  //                                                        //    <br>&gt; &gt;&gt; nr       type          resnr  residue  atom   cgnr   charge     mass  <br>&gt; &gt;&gt; //      //2 amber99_17      1   NGLY     H1      2     0.1642      <br>&gt; &gt;&gt; 1.008  //      3 amber99_17      1   NGLY     H2      2     <br>&gt; &gt;&gt; 0.1642      1.008  //      4 amber99_17      1   NGLY     H3      <br>&gt; &gt;&gt; 2     0.1642      1.008   //   /<br>&gt; &gt;&gt; thanks for your attention        <br>&gt; &gt; <br>&gt; &gt; <br>&gt; <br>&gt; -- <br>&gt; ========================================<br>&gt; <br>&gt; Justin A. Lemkul<br>&gt; Ph.D. Candidate<br>&gt; ICTAS Doctoral Scholar<br>&gt; MILES-IGERT Trainee<br>&gt; Department of Biochemistry<br>&gt; Virginia Tech<br>&gt; Blacksburg, VA<br>&gt; jalemkul[at]vt.edu | (540) 231-9080<br>&gt; http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin<br>&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/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/mailing_lists/users.php<br>                                               <br /><hr />New Windows 7: Find the right PC for you. <a href='http://windows.microsoft.com/shop' target='_new'>Learn more.</a></body>
</html>