<table cellspacing="0" cellpadding="0" border="0" ><tr><td valign="top" style="font: inherit;"><DIV>Dear GMX users,</DIV>
<DIV>&nbsp; Recently, I want to caculate the free energy for K+-----&gt;NA+ in bulk with GMX4.5.4 using charmm 27 force field. &nbsp;I change ion.itp</DIV>
<DIV>[ moleculetype ]<BR>; molname&nbsp;nrexcl<BR>K&nbsp;&nbsp;1</DIV>
<DIV>[ atoms ]<BR>; id&nbsp;at type&nbsp;&nbsp;&nbsp; res nr &nbsp;residu name&nbsp;at name&nbsp; cg nr&nbsp;&nbsp; charge&nbsp;&nbsp; </DIV>
<DIV>1&nbsp;POT&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp; K&nbsp;&nbsp;&nbsp;&nbsp; K&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp; </DIV>
<DIV>&nbsp;</DIV>
<DIV>to</DIV>
<DIV>&nbsp;</DIV>
<DIV>[ moleculetype ]<BR>; molname&nbsp;nrexcl<BR>K&nbsp;&nbsp;1</DIV>
<DIV>[ atoms ]<BR>; id&nbsp;at type&nbsp;&nbsp;&nbsp; res nr &nbsp;residu name&nbsp;at name&nbsp; cg nr&nbsp;&nbsp; charge&nbsp;&nbsp; mass&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; typeB&nbsp;&nbsp;&nbsp; chargeB&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; massB<BR>1&nbsp;POT&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp; K&nbsp;&nbsp;&nbsp;&nbsp; K&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp; 39.102000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; SOD&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 22.989770 </DIV>
<DIV>#ifdef POSRES_K<BR>[ position_restraints ]<BR>;&nbsp; i funct&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; fcx&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; fcy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; fcz<BR>&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1000000&nbsp;&nbsp;&nbsp; 1000000&nbsp;&nbsp;&nbsp;&nbsp; 1000000<BR>#endif</DIV>
<DIV>&nbsp;</DIV>
<DIV>and include it into .top file.</DIV>
<DIV>&nbsp;</DIV>
<DIV>And I set .mdp file</DIV>
<DIV>; Free energy control stuff<BR>free-energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = yes<BR>init-lambda&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.0<BR>delta-lambda&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<BR>foreign_lambda&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.05<BR>sc-alpha&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.5<BR>sc-power&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.0<BR>sc-sigma&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.3<BR>nstdhdl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 10<BR>separate-dhdl-file&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
 yes<BR>dhdl-derivatives&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = yes<BR>dh_hist_size&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<BR>dh_hist_spacing&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.1<BR>couple-moltype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = K<BR>couple-lambda0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = vdw-q<BR>couple-lambda1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = vdw-q<BR>couple-intramol&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = no</DIV>
<DIV>&nbsp;</DIV>
<DIV>however when I do EM, there are two warnings and stop.</DIV>
<DIV>WARNING 1 [file min.mdp, line 329]:<BR>&nbsp; The lambda=0 and lambda=1 states for coupling are identical</DIV>
<DIV>Generated 28 of the 28 non-bonded parameter combinations<BR>Generating 1-4 interactions: fudge = 1<BR>Generated 28 of the 28 1-4 parameter combinations</DIV>
<DIV>WARNING 2 [file ion.itp, line 22]:<BR>&nbsp; Some parameters for bonded interaction involving perturbed atoms are<BR>&nbsp; specified explicitly in state A, but not B - copying A to B</DIV>
<DIV>Excluding 3 bonded neighbours molecule type 'gra'<BR>Excluding 2 bonded neighbours molecule type 'SOL'<BR>Excluding 1 bonded neighbours molecule type 'K'<BR>Coupling 1 copies of molecule type 'K'</DIV>
<DIV>NOTE 1 [file topol.top, line 16]:<BR>&nbsp; System has non-zero total charge: 1.000000e+00<BR>&nbsp; </DIV>
<DIV><BR>Analysing residue names:<BR>There are:&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Other residues<BR>There are:&nbsp;&nbsp; 531&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Water residues<BR>There are:&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ion residues<BR>Analysing residues not classified as Protein/DNA/RNA/Water and splitting into groups...<BR>Analysing residues not classified as Protein/DNA/RNA/Water and splitting into groups...<BR>Number of degrees of freedom in T-Coupling group rest is 3660.00<BR>Largest charge group radii for Van der Waals: 0.040, 0.040 nm<BR>Largest charge group radii for Coulomb:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.079, 0.079 nm<BR>Calculating fourier grid dimensions for X Y Z<BR>Using a fourier grid of 20x20x30, spacing 0.115 0.115 0.115<BR>Estimate for the relative computational load of the PME mesh part: 0.36<BR>This run will generate roughly 22 Mb of data</DIV>
<DIV>There was 1 note</DIV>
<DIV>There were 2 warnings</DIV>
<DIV>-------------------------------------------------------<BR>Program grompp, VERSION 4.5.4<BR>Source code file: grompp.c, line: 1582</DIV>
<DIV>Fatal error:<BR>Too many warnings (2), grompp terminated.<BR>If you are sure all warnings are harmless, use the -maxwarn option.<BR>For more information and tips for troubleshooting, please check the GROMACS<BR>website at <A href="http://www.gromacs.org/Documentation/Errors">http://www.gromacs.org/Documentation/Errors</A><BR>-------------------------------------------------------</DIV>
<DIV>"Your Proposal is Accepted" (Men In Black)</DIV>
<DIV>&nbsp;</DIV>
<DIV>It seems that </DIV>
<DIV>[ moleculetype ]<BR>; molname&nbsp;nrexcl<BR>K&nbsp;&nbsp;1</DIV>
<DIV>[ atoms ]<BR>; id&nbsp;at type&nbsp;&nbsp;&nbsp; res nr &nbsp;residu name&nbsp;at name&nbsp; cg nr&nbsp;&nbsp; charge&nbsp;&nbsp; mass&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; typeB&nbsp;&nbsp;&nbsp; chargeB&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; massB<BR>1&nbsp;POT&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp; K&nbsp;&nbsp;&nbsp;&nbsp; K&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp; 39.102000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; SOD&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 22.989770 </DIV>
<DIV>do not take effect.</DIV>
<DIV>&nbsp;</DIV>
<DIV>Thanks in advance!</DIV>
<DIV>&nbsp;</DIV>
<DIV>Zhongjin He</DIV></td></tr></table>