Actually I solved the problem, <div>I did a mistake in the bond conversion. Actually I found the right formula:</div><div><br></div><div>k_gromacs = k_charmm*<b>418.4</b>*2</div><div><br></div><div>rising form the conversion from angstrom to nanometer....</div>
<div><br></div><div><br></div><div><br></div><div><br><div class="gmail_quote">---------- Messaggio inoltrato ----------<br>Da: <b class="gmail_sendername">francesco oteri</b> <span dir="ltr"><<a href="mailto:francesco.oteri@gmail.com">francesco.oteri@gmail.com</a>></span><br>
Date: 08 marzo 2012 10:49<br>Oggetto: charm to gromacs<br>A: Discussion list for GROMACS users <<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>><br><br><br>Dear gromacs users, <div>I am trying to convert a charmm style parameter file into gromacs. So I al trying to understand how converting the force constants from CHARMM to gromacs.</div>
<div>Using the CHARMM27 force field as a reference, I've noticed that to convert the bond constant the following formula is supposed to be used:</div>
<div><br></div><div>k_gromacs = k_charmm*<b>4184</b>*2</div><div>In fact if I try to convert the following line from par_all27_prot_lipid.prm charmm file:</div><div>CEL1 CEL1 440.000 1.340 ! butene, yin,adm jr., 12/95<br clear="all">
<div><br></div><div>I obtain the right result ( 440.000*4184*2. = 368192.0)</div><div>charmm27/ffbonded.itp CEL1 CEL1 1 <a href="tel:0.134%20%C2%A0%20368192.0" value="+3901343681920" target="_blank">0.134 368192.0</a> </div>
<div><br></div><div><br></div><div>But to convert the angle constants I've to use the following formula:</div>
<div><br></div><div>k_gromacs = k_charmm*<b>4.184</b>*2.</div><div><br></div><div>how demostrated by the conversion of the following line:</div><div><div><div><br></div><div>par_all27_prot_lipid.prm CEL1 CEL1 CTL2 48.00 123.50 ! from 2-butene, yin,adm jr., 12/95</div>
</div><div><br></div><div>into</div><div><br></div><div>charmm27/ffbonded.itp CEL1 CEL1 CTL2 5 123.50 401.664 0.0 0.0</div></div><div><br></div><div><br></div><div><br></div><div>My question is: </div>
<div><br></div><div>Why in the former case I need the constant 4184 while in the latter the constant 4.184 is supposed to be used? Am I miss something?</div><br>
</div>
</div><br><br clear="all"><div><br></div>-- <br>Cordiali saluti, Dr.Oteri Francesco<br>
</div>